Sinking Tide Gauge Revealed by Space-borne InSAR : Implications for Sea Level Acceleration at Pohang , South Korea

Vertical land motion at tide gauges influences sea level rise acceleration; this must be addressed for interpreting reliable sea level projections. In recent years, tide gauge records for the Eastern coast of Korea have revealed rapid increases in sea level rise compared with the global mean. Pohang Tide Gauge Station has shown a +3.1 cm/year sea level rise since 2013. This study aims to estimate the vertical land motion that influences relative sea level rise observations at Pohang by applying a multi-track Persistent Scatter Interferometric Synthetic Aperture Radar (PS-InSAR) time-series analysis to Sentinel-1 SAR data acquired during 2015–2017. The results, which were obtained at a high spatial resolution (10 m), indicate vertical ground motion of −2.55 cm/year at the Pohang Tide Gauge Station; this was validated by data from a collocated global positioning system (GPS) station. The subtraction of InSAR-derived subsidence rates from sea level rise at the Pohang Tide Gauge Station is 6 mm/year; thus, vertical land motion significantly dominates the sea level acceleration. Natural hazards related to the sea level rise are primarily assessed by relative sea level changes obtained from tide gauges; therefore, tide gauge records should be reviewed for rapid vertical land motion along the vulnerable coastal areas.


Introduction
Sea level rise (SLR) is a global phenomenon that threatens low-lying coasts and islands [1].Global climate change and anthropogenic activities are considered vital factors for global and regional sea level rise [2].Until now, long-term global mean sea level rise and regional sea level rise have mainly been estimated using sea level measurements provided by tide gauges (TG) installed around the world [3,4].Tide gauges provide multi-decadal local sea level records of historical sea level variations which are utilized to project future sea level rise [5].The spatial and temporal variations of the tide gauges are used to investigate rapid sea level accelerations and the potential impacts on coastal areas [6][7][8].In addition, tide gauges are often utilized to calibrate and validate the sea level estimates from satellite altimeters.Since the early 1990s, satellite altimeters such as TOPEX/Poseidon and Jason-1, 2, and 3 have been employed in global mean sea level studies (GMSL) and are used to evaluate the model projections of SLR [9][10][11][12].In addition to the T/P and Jason satellites, the successors of these kinds of satellites, such as ERS-1 & 2, SARAL/Altika, Cryosat, HY-2A, and Sentinel-3 a/b have provided global and regional sea level records for the past 25 years that allow us to detect spatiotemporal variations [13].These long-term altimeter records are calibrated with the tide gauges to provide accurate sea level records.However, oftentimes, inconsistencies between tide gauge records among neighboring tide gauge stations arise from coastal processes, climate variability, tectonics, and subsidence effects [6,7].In particular, in delta regions or basins, the relative sea level is often dominated by subsidence due to both long-term glacial isostatic adjustments (GIA) and recent compaction of deltaic sediments due to natural and anthropogenic processes.Anthropogenic subsidence caused by groundwater pumping, oil/gas extraction, and sediment loading also perturbs the sea level records.As tide gauges measure relative SLR with reference to a tidal benchmark grounded on the adjacent coast, even minor ground motions will significantly affect sea level height measurements.Previous studies have reported on vulnerable sites for coastal subsidence and short-term sea level rise, including Hampton Roads, Virginia [14], New York City [15], Taiwan [16], and North Carolina [17].Hence, the estimation and correction of the local ground motion at and around the tide gauges is important to minimize uncertainties in the local and global sea level changes observed [18][19][20][21][22][23].In recent decades, long-term glacial adjustments were quantified using historical variations and climate models for the reconstruction of sea level records [24,25].However, for vertical ground motion at tide gauge sites, few attempts have been made to attribute the dominance of the local subsidence rate, as tide gauges around the world are not always equipped with global positioning system (GPS) instruments [26,27].
According to the Korea Hydrographic and Oceanographic Agency's report (KHOA), the mean SLR along the entire coast of Korea during 1971-2010 was 2.48 mm/year, which is greater than the global mean SLR (i.e., 2 ± 0.3 mm/year) [28].The South Korean peninsula is located in the Southeastern part of Asia; several tide gauge stations along the coast continuously record sea level changes (Figure 1).Among these stations, Pohang, Ulsan, Mokpo, Busan, and Mukho have been recording for more than 50 years [29].A number of tide gauge stations on the East coast of Korea have shown inconsistent relative SLR in recent years [28,29].Among them, the Pohang Tide Gauge Station has observed the highest relative SLR since the year 2010.The Pohang Tide Gauge Station is located along the breakwater structure at Guhang Port, Pohang (Figure 1).It is thought that the Pohang tide gauge is considerably impacted by surface displacements, which may influence sea level records.Some tide gauge sites use GPS to measure reliable absolute ground motion in the geocentric reference frame [8,30].However, although this approach is widely accepted, very few tide gauges are equipped with permanent GPS stations.Recently, other approaches have been presented to estimate vertical land motion using short-term satellite altimetry data and long-term sea level change patterns recorded by tide gauges [19,31].Consequently, the use of space-borne interferometric Synthetic Aperture Radar (InSAR) techniques has gained importance in recent years [20].InSAR techniques are widely accepted to estimate Earth's surface displacements [32].Time-series InSAR analysis methods can estimate relative displacements at mm-level accuracy; for example, Persistent Scatter Interferometry (PSI) [33,34] and Small Baseline Subset (SBAS) [35].Unlike GPS, which measures precise vertical displacement at specific locations in the geocentric frame, time-series InSAR estimates relative displacements with high spatial and temporal resolution.Wöppelmann,et al. [36] explored the predominance of line-of-sight (LOS) land motion on sea level change in coastal areas using the PS-InSAR technique.However, there remains a need to estimate precise vertical and horizontal ground motion at tide gauges to apply corrections to the measured sea level changes.
In this study, we aimed to estimate vertical land motion at the Pohang Tide Gauge Station by incorporating a multi-track PS-InSAR time-series analysis in order to evaluate the driver for the overestimated sea level rise.The reliability of the vertical velocity was validated against GPS observations collocated at the tide gauge station.The results confirm that vertical land motion at the tide gauge station significantly influences the observed sea level changes.The difference between the rate of vertical ground motion and sea level rise is used to reconstruct a corrected rate of sea level rise since 2013.

Study Area
The Pohang Tide Gauge Station is located on the Southeastern coast of the Korean peninsula, in the Pohang Basin (Figure 1).The coastline borders the East Sea, a marginal sea between the Korean peninsula, Japan, and Russia.The Pohang Basin has witnessed extensive land reclamation over past years, as the city is one of the economic zones of Korea.In particular, the Guhang Port area, where the tide gauge is grounded, is built on the reclaimed land of the Hyongsan River.The Pohang Tide Gauge Station has measured sea level records continuously from 1977-2018 (KHOA).Based on the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4), Pohang has shown increasing trends of extreme precipitation; however, its relationship with low-frequency climate events such as the El Niño-Southern Oscillation (ENSO) still need to be investigated [37].Seasonal analysis of long-term meteorological and hydrological variables suggests that the high intensity of precipitation along the East coast is mainly due to typhoons and extreme rain events [38].In recent years, a permanent GPS station (POHA) has been collocated at the Pohang Tide Gauge Station; it is maintained by KHOA.

Study Area
The Pohang Tide Gauge Station is located on the Southeastern coast of the Korean peninsula, in the Pohang Basin (Figure 1).The coastline borders the East Sea, a marginal sea between the Korean peninsula, Japan, and Russia.The Pohang Basin has witnessed extensive land reclamation over past years, as the city is one of the economic zones of Korea.In particular, the Guhang Port area, where the tide gauge is grounded, is built on the reclaimed land of the Hyongsan River.The Pohang Tide Gauge Station has measured sea level records continuously from 1977-2018 (KHOA).Based on the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4), Pohang has shown increasing trends of extreme precipitation; however, its relationship with low-frequency climate events such as the El Niño-Southern Oscillation (ENSO) still need to be investigated [37].Seasonal analysis of long-term meteorological and hydrological variables suggests that the high intensity of precipitation along the East coast is mainly due to typhoons and extreme rain events [38].In recent years, a permanent GPS station (POHA) has been collocated at the Pohang Tide Gauge Station; it is maintained by KHOA.

Methodology: SAR Data Collection and Time-Series PS-InSAR
InSAR is a widely accepted technique for estimating the Earth's surface displacements at submillimeter accuracy and a high spatial resolution [34].We used 49 SAR scenes in ascending tracks and 42 scenes in descending tracks acquired by the C-band Sentinel-1 A/B from 2015-2017.Sentinel-1 SAR

Methodology: SAR Data Collection and Time-Series PS-InSAR
InSAR is a widely accepted technique for estimating the Earth's surface displacements at sub-millimeter accuracy and a high spatial resolution [34].We used 49 SAR scenes in ascending tracks and 42 scenes in descending tracks acquired by the C-band Sentinel-1 A/B from 2015-2017.Sentinel-1 SAR sensors have the ability to cover a large spatial area (i.e., 250 km swath) and have a revisit cycle of 6 days.Details of data acquisition are shown in Table 1.A stack of multi-look interferograms (2 range × 1 Azimuth looks) resulting in a spatial resolution of 10 × 20 m was generated.Sentinel-1 Terrain Observation by Progressive Scans (TOPS) InSAR processing was performed to generate interferograms with a single master scene using the ISCE software [39,40].SAR scenes acquired on 19 February 2017 and 14 April 2017 were chosen as the master scenes for ascending and descending stacks, respectively.Master scenes were selected based on the minimum perpendicular and temporal baselines, as shown in Figure 2 [41].A Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) with 30-m spatial resolution was used to remove the topographic phase effect from the interferometric phase.In order to estimate the time-series linear trend of sea level change, we obtained hourly sea level records of the Pohang Tide Gauge Station from KHOA.The sea level records were corrected for local datum over the period.Sea-level trend analyses were performed in order to decompose trends, seasonality, and residual components from observed sea level records with 95% confidence intervals to the linear trend.The 95% confidence intervals on the linear trend were calculated by the standard statistical method.Sea level records are available from 1977-2017; however, the hourly SLR data provided by KHOA only show a rapid acceleration in sea level records since 2013.Therefore, in this study, we focused on recent records (i.e., 2013-2017).In addition, continuous POHA station GPS observations along the N-S, E-W, and vertical directions since 2013 were provided by KHOA (Figure 3).The data show that the Pohang Tide Gauge Station has remained stable for horizontal movement, but is significantly impacted by vertical displacement.Although GPS observations are discontinuous during 2015 and 2017, the displacement pattern shows a linear trend over the study period.Hence, the vertical displacement trend from GPS observations was utilized for comparison with InSAR-derived vertical displacements.Time-series PSI processing was applied using the Stanford Method of Persistent Scatterer (StaMPS) to estimate ground displacement over the study period [34,41].In this method, initial PS candidates are selected based on the amplitude dispersion index criterion.In this step, if the ratio between the standard deviation and mean amplitude in time is smaller than a threshold value of 0.4, then such pixels are selected as initial PS candidates to form the reference network for further processing.Following that, the phase values of these PS candidates are filtered with neighboring pixels assuming that the deformation, orbital error, and atmospheric artifacts are spatially correlated by bandpass filtering [34].Then, persistent scatterers are iteratively selected as target pixels found to be phase stable over a period using a measure called temporal coherence.Residual PS candidates having temporal coherence of <0.3 are filtered out at this stage.The resulting PS pixels have minimized spatial and temporal decorrelation noise levels.
The StaMPS PS-InSAR algorithm incorporates advanced time-series processing steps to reduce the decorrelation of noise in the differential interferograms.Furthermore, residual topographic errors that are spatially uncorrelated but correlated with perpendicular baseline information were removed by the least square method before phase unwrapping.The wrapped phase of the PS pixels was unwrapped by 3D (phase in space and time field) unwrapping using the SNAPHU algorithm [42].The thresholds were carefully selected for filter grid size, and an unwrapped grid size of 5 × 5 was used to control phase unwrapping errors.Spatial filtering was then applied for re-estimation of atmospheric and look-angle errors.Subsequently, the time-series surface displacement rate was estimated with respect to a reference point that presumed a nearly-stable ground location.
In order to estimate the time-series linear trend of sea level change, we obtained hourly sea level records of the Pohang Tide Gauge Station from KHOA.The sea level records were corrected for local datum over the period.Sea-level trend analyses were performed in order to decompose trends, seasonality, and residual components from observed sea level records with 95% confidence intervals to the linear trend.The 95% confidence intervals on the linear trend were calculated by the standard statistical method.Sea level records are available from 1977-2017; however, the hourly SLR data provided by KHOA only show a rapid acceleration in sea level records since 2013.Therefore, in this study, we focused on recent records (i.e., 2013-2017).In addition, continuous POHA station GPS observations along the N-S, E-W, and vertical directions since 2013 were provided by KHOA (Figure 3).The data show that the Pohang Tide Gauge Station has remained stable for horizontal movement, but is significantly impacted by vertical displacement.Although GPS observations are discontinuous during 2015 and 2017, the displacement pattern shows a linear trend over the study period.Hence, the vertical displacement trend from GPS observations was utilized for comparison with InSAR-derived vertical displacements.

Vertical Ground Motion at Pohang Tide Gauge Station
The mean velocity map obtained from the PSI analysis is presented in Figure 4.These PSI results indicate ground motion along the satellite LOS.PS pixels were identified in the urban areas of Pohang, and especially on the Pohang tide gauge.The time-series PSI results reveal that the vicinity of the Pohang Tide Gauge Station is affected by significant relative ground motion away from the satellite.Mean LOS displacement rates of approximately −1.85 cm/year and −1.92 cm/year were observed in ascending (Figure 4a) and descending (Figure 4b) passes, respectively.Other than the Pohang Tide Gauge Station and its neighborhood, most of the study area appeared stable, with minimum downward ground motion and no significant uplift based on PSI results.Time-series profiles show a linear LOS displacement trend at the tide gauge location over the study period for both ascending and descending tracks (Figure 4c).It is worth noting that the LOS displacements at the tide gauge station for both passes

Vertical Ground Motion at Pohang Tide Gauge Station
The mean velocity map obtained from the PSI analysis is presented in Figure 4.These PSI results indicate ground motion along the satellite LOS.PS pixels were identified in the urban areas of Pohang, and especially on the Pohang tide gauge.The time-series PSI results reveal that the vicinity of the Pohang Tide Gauge Station is affected by significant relative ground motion away from the satellite.Mean LOS displacement rates of approximately −1.85 cm/year and −1.92 cm/year were observed in ascending (Figure 4a) and descending (Figure 4b) passes, respectively.Other than the Pohang Tide Gauge Station and its neighborhood, most of the study area appeared stable, with minimum downward ground motion and no significant uplift based on PSI results.Time-series profiles show a linear LOS displacement trend at the tide gauge location over the study period for both ascending and descending tracks (Figure 4c).It is worth noting that the LOS displacements at the tide gauge station for both passes showed a similar and constant linear trend.Using the area that overlapped in both the spatial and temporal coverage of Sentinel-1 for the ascending/descending data, PSI results with respect to a common reference point were interpolated by cubic interpolation on the uniform grid before 2D decomposition.With the availability of LOS displacement for the two different geometries, the vertical (dver) and horizontal displacement components (dhor) were reconstructed, as presented in Equation 1 [43]: Using the area that overlapped in both the spatial and temporal coverage of Sentinel-1 for the ascending/descending data, PSI results with respect to a common reference point were interpolated by cubic interpolation on the uniform grid before 2D decomposition.With the availability of LOS displacement for the two different geometries, the vertical (d ver ) and horizontal displacement components (d hor ) were reconstructed, as presented in Equation ( 1) [43]: where d asc and d dsc are displacement along LOS in ascending and descending mode, respectively; θ asc and θ dsc are incidence angle in ascending and descending mode, respectively; α asc , α dsc are the satellite heading angle for ascending and descending passes; and d ver and d hor are vertical displacement and horizontal displacement, respectively, along the East-West direction.
The projections of the n,e,u vector components of the LOS velocities are presented in Figure 5.The Azimuth look direction (ALD) is perpendicular to the satellite heading and comprises the combined contributions of the E-W and N-S components.In the across-track SAR imaging system, sensitivity towards N-S components is much smaller than sensitivity towards E-W components.The decomposed vertical and horizontal displacements were obtained according to the incidence angles of ascending and descending passes.Hence, we assumed that the derived horizontal component represents the E-W velocity component along the descending Azimuth look direction.The decomposed vertical and horizontal components for the overlapping period (May 2015-December 2017) are presented in Figure 6a,b, respectively.It is clear that ground motions detected over the study area by PSI analysis have primarily undergone vertical movement with minimum to zero horizontal movements.The velocities of ground motion in the vertical and horizontal directions are approximately −2.55 cm/year and +0.01 cm/year, respectively.
where dasc and ddsc are displacement along LOS in ascending and descending mode, respectively; θasc and θdsc are incidence angle in ascending and descending mode, respectively;  ,  are the satellite heading angle for ascending and descending passes; and dver and dhor are vertical displacement and horizontal displacement, respectively, along the East-West direction.
The projections of the n,e,u vector components of the LOS velocities are presented in Figure 5.The Azimuth look direction (ALD) is perpendicular to the satellite heading and comprises the combined contributions of the E-W and N-S components.In the across-track SAR imaging system, sensitivity towards N-S components is much smaller than sensitivity towards E-W components.The decomposed vertical and horizontal displacements were obtained according to the incidence angles of ascending and descending passes.Hence, we assumed that the derived horizontal component represents the E-W velocity component along the descending Azimuth look direction.The decomposed vertical and horizontal components for the overlapping period (05/2015-12/2017) are presented in Figure 6a,b, respectively.It is clear that ground motions detected over the study area by PSI analysis have primarily undergone vertical movement with minimum to zero horizontal movements.The velocities of ground motion in the vertical and horizontal directions are approximately −2.55 cm/year and +0.01 cm/year, respectively.Temporal profiles for vertical ground motions estimated from POHA GPS and PS-InSAR analyses and the non-seasonal sea level trend computed from the Pohang tide gauge are presented in Figure 7.The vertical and horizontal components of the continuous GPS station (POHA) collocated with the Pohang Tide Gauge Station were estimated.The vertical velocity at POHA was found to be −2.67 cm/year from 2013-2016 with respect to the 2008 International Terrestrial Reference Frame (ITRF) coordinate system (Figure 7), for which the horizontal movements of GPS stations are stable.Temporal profiles for vertical ground motions estimated from POHA GPS and PS-InSAR analyses and the non-seasonal sea level trend computed from the Pohang tide gauge are presented in Figure 7.The vertical and horizontal components of the continuous GPS station (POHA) collocated with the Pohang Tide Gauge Station were estimated.The vertical velocity at POHA was found to be −2.67 cm/year from 2013-2016 with respect to the 2008 International Terrestrial Reference Frame (ITRF) coordinate system (Figure 7), for which the horizontal movements of GPS stations are stable.

Discussion
With the continued melting of ice sheets and glaciers in recent years, due to ocean and atmospheric warming, sea level is rising at an increasing rate.Based on tide gauge records from the past century and from satellite altimetry, the global mean sea level rises by 1.7 ± 0.3 mm/year and 3.3 mm/year, respectively [28].Nevertheless, local sea level change at the Pohang Tide Gauge Station deviates significantly from regional and global sea level change, and from sea level records at nearby tide gauge stations (e.g., Ulsan, Mukho).This is because such sea level measurements may include large contributions from ground motion.Vertical ground motion estimates obtained by various approaches have been reported for various sites; for example, Alexandria, Egypt [36] and Manila, and the Philippines [23].However, these studies estimated the ground motion in the LOS directions at tide gauge stations.

Discussion
With the continued melting of ice sheets and glaciers in recent years, due to ocean and atmospheric warming, sea level is rising at an increasing rate.Based on tide gauge records from the past century and from satellite altimetry, the global mean sea level rises by 1.7 ± 0.3 mm/year and 3.3 mm/year, respectively [28].Nevertheless, local sea level change at the Pohang Tide Gauge Station deviates significantly from regional and global sea level change, and from sea level records at nearby tide gauge stations (e.g., Ulsan, Mukho).This is because such sea level measurements may include large contributions from ground motion.Vertical ground motion estimates obtained by various approaches have been reported for various sites; for example, Alexandria, Egypt [36] and Manila, and the Philippines [23].However, these studies estimated the ground motion in the LOS directions at tide gauge stations.
Time-series PSI results reveal vulnerable areas of the Pohang coast; that is, those experiencing vertical (downward) ground motions (Figure 7b).In particular, the breakwater where the Pohang Tide Gauge Station is located is subsiding at a rate of 2.55 cm/year, with InSAR displacements from both ascending and descending tracks showing similar displacement patterns (i.e., a linear subsidence trend).Furthermore, the reliability of the PSI ground motion trend is in good agreement with a continuous GPS vertical velocity profile from a station (POHA) collocated with the Pohang Tide Gauge Station (Figure 7).The Guhang Port area where the tide gauge is grounded was built on the reclaimed land of the Hyongsan River; ground motion at the reclaimed land is primarily induced by the compaction of underlying clay soils [44].
After removing seasonality from the observed tide gauge observations, the non-seasonal sea level trend showed significant SLR of about +3.1 cm/year at the Pohang Tide Gauge Station (Figure 7c).The slope of the SLR has been large in recent years compared with that of long-term sea level rise [19,28].It is worth remarking that the sea level acceleration trend is possibly the combination of original SLR and local ground motion.Therefore, we infer that the rapid increase in SLR is most likely the result of vertical ground motion along the breakwater.This is substantiated by the vertical ground motion observed from the POHA GPS measurements (Figure 7).The difference between the linear sea level trend (+3.1 cm/year) and the estimated vertical ground motion (−2.55 cm/year) since 2013 (i.e., ~6 mm/year) reflects regional sea level rise around the Korean peninsula.However, 6 mm/year is still higher than the global mean SLR.According to KHOA [28], the SLR at the Pohang Tide Gauge Station was 5.82 mm/year from 1977-2015, hence, the SLR rate of ~6 mm/year is the considered the long-term rate.However, our PSI results represent short-term SLR during the InSAR acquisition period (2015-2017); therefore, our estimated vertical ground motion and subsequent reconstructed SLR rate (6 mm/year) apply only during the 2015-2017 period.We believe that this SLR rate of 6 mm/year reflects the influence of subsidence due to long-term natural or anthropogenic processes and seasonal climate variability [45,46].
Although tide gauges play a vital role in observing relative sea level changes, the observations are limited to land motion, and by sparse distributions around the world.Several studies have provided the sea level reconstruction for post-glacial rebound and other natural processes [31,47].Nevertheless, until now, rapid land motion at tide gauges caused by anthropogenic processes have not been considered for reconstructions of rapid SLR.Without a doubt, contemporary GPS measurements at the tide gauge stations can provide continuous and reliable information on ground motion rates.Unfortunately, not all tide gauge stations are collocated with permanent GPS stations.Where GPS stations do exist, they may not provide continuous observations.Therefore, the vertical ground motion of the tide gauge stations is not clear enough to reconstruct the original sea level rise.In recent years, InSAR data with high temporal and spatial resolutions for both ascending and descending data have become freely available.InSAR products are useful for the continuous monitoring of vertical ground motion at tide gauge stations and in their vicinities.

Conclusions
Despite the importance of sea level measurement accuracy, unexpected local factors can deteriorate sea level records and result in unreliable interpretations.The Pohang Tide Gauge Station has been recording sea level rises that are higher than those of neighboring stations and of the global mean.Here, we have presented an estimation of vertical ground motion at the Pohang Tide Gauge Station using multi-track time-series InSAR analysis; the results verify that ground motion is a vital factor in the measured sea level change, resulting in an overestimated sea level trend (+3.1 cm/year since 2013).The estimated vertical ground motion (−2.55 cm/year) was validated against the ground motion rate from GPS observations (−2.67 cm/year) at a collocated station.The multi-track PS-InSAR technique can be applied to any tide gauge located in a coastal area to determine the influence of

Figure 1 .
Figure 1.Study area map of Guhang Port, Pohang city (aerial imagery of Pohang city obtained from the National Geographic Information Institute (NGII).The red dot denotes the location of the Pohang Tide Gauge Station and the POHA GPS station.The inset graph shows a time-series profile of monthly sea level rise from 1977-2017 (KHOA).

Figure 1 .
Figure 1.Study area map of Guhang Port, Pohang city (aerial imagery of Pohang city obtained from the National Geographic Information Institute (NGII).The red dot denotes the location of the Pohang Tide Gauge Station and the POHA GPS station.The inset graph shows a time-series profile of monthly sea level rise from 1977-2017 (KHOA).

Figure 3 .
Figure 3. Global positioning system (GPS) displacements measured at the POHA station from 2013 to 2016 (KHOA).(a) North-South, (b) East-West, and (c) vertical displacements.Horizontal components are measured in degrees and vertical displacement is measured in meters.

Figure 3 .
Figure 3. Global positioning system (GPS) displacements measured at the POHA station from 2013 to 2016 (KHOA).(a) North-South, (b) East-West, and (c) vertical displacements.Horizontal components are measured in degrees and vertical displacement is measured in meters.

Figure 4 .
Figure 4. Line-of-sight (LOS) velocities from time-series Interferometric Synthetic Aperture Radar (InSAR) results obtained from (a) ascending (November 2014-December 2017) and (b) descending (May 2015-May 2018) tracks.Negative values indicate displacement away from the satellite.(c) Time-series InSAR displacement at Pohang Tide Gauge Station for ascending (blue) and descending (green) tracks.

Figure 5 .
Figure 5. Projection of Interferometric Synthetic Aperture Radar (InSAR) line-of-sight (LOS) velocity decomposition into 2D velocity components along the vertical and horizontal directions (the Azimuth look direction for the descending heading).

Figure 5 .
Figure 5. Projection of Interferometric Synthetic Aperture Radar (InSAR) line-of-sight (LOS) velocity decomposition into 2D velocity components along the vertical and horizontal directions (the Azimuth look direction for the descending heading).

Figure 6 .
Figure 6.Interferometric Synthetic Aperture Radar (InSAR)-derived ground motion velocities from May 2015 to December 2017.(a) Vertical ground motion component and (b) horizontal ground motion component along the look direction (descending).

Figure 7 .
Figure 7. Temporal profiles of (a) absolute vertical ground motion from global positioning system (GPS) data, (b) relative vertical ground motion from Interferometric Synthetic Aperture Radar (InSAR) displacements, and (c) sea level records observed at Pohang Tide Gauge Station since 2013.

Figure 7 .
Figure 7. Temporal profiles of (a) absolute vertical ground motion from global positioning system (GPS) data, (b) relative vertical ground motion from Interferometric Synthetic Aperture Radar (InSAR) displacements, and (c) sea level records observed at Pohang Tide Gauge Station since 2013.

Table 1 .
Sentinel-1 A/B SAR scenes used in this study.