Integrated space geodesy for mapping land deformation over Choushui River Fluvial Plain, Taiwan

ABSTRACT Groundwater is an important part of the precious water resources. As the fresh surface water resources become scarcer because of climate change, population growth, and industrial activities, more and more groundwater has been extracted to meet the demands of various water uses (e.g. municipal, industrial, and agricultural). Excessive groundwater extraction leads to severe ground subsidence which compromises the safety of surface and underground infrastructures. Modelling the effects of groundwater extraction is vital to the management and sustainable use of groundwater. However, results of such modelling have to be validated with inputs such as the field survey of ground subsidence. Levelling and continuous global positioning system (GPS) receiver networks are routinely used to collect these field measurements. Unfortunately, these techniques have limitations in terms of areal coverage and density of survey marks and, as a result, subsidence hot spots can be easily missed out. In order to provide a comprehensive picture of subsidence to aid geotechnical modelling and to assess the effectiveness of measures used to mitigate ground subsidence, satellite imaging radar interferometry techniques (interferometric synthetic aperture radar (InSAR) can be used to complement other deformation monitoring techniques. In this study, 20 Advanced Land Observing Satellite (ALOS) Phased Array L-band Synthetic Aperture Radar (PALSAR) images acquired from 31 December 2006 to 26 February 2011 were used to map the land displacement over the Choushui River Fluvial Plain (CRFP), Taiwan. The GPS measurements acquired at 10 continuously operating reference stations (CORS) were used to refine the orbit error in the each differential interferogram obtained from each radar image pair. The displacement time series over the distributed scatterers and the persistent scatterers were analysed. Several subsidence bowls were identified in CRFP. A quantitative comparison was conducted to compare the radar measurements to the GPS measurements over 36 GPS CORS stations. Good agreement between both measurements was observed with coefficient of determination (R2) of 0.97, absolute mean difference of 3.2 mm year−1, and standard deviation of 4 mm year−1. The InSAR-measured Line-of-Sight displacement and GPS-measured horizontal displacement were integrated to derive the vertical displacement map. Two displacement maps were generated using two ALOS-2 PALSAR-2 pairs acquired between 2015 and 2016. Similar subsidence patterns were found in the two maps compared to the 2006–2011 displacement rate map, suggesting the land over the same region might have continued to fall.

[1]  Xiaojing Li,et al.  Land subsidence characteristics of Bandung Basin as revealed by ENVISAT ASAR and ALOS PALSAR interferometry , 2014 .

[2]  Stefano Tebaldini,et al.  On the Exploitation of Target Statistics for SAR Interferometry Applications , 2008, IEEE Transactions on Geoscience and Remote Sensing.

[3]  L. Mervart,et al.  Bernese GPS Software Version 5.0 , 2007 .

[4]  Zheyuan Du,et al.  Satellite-based Estimates of Ground Subsidence in Ordos Basin, China , 2016 .

[5]  T. Lee,et al.  High-resolution sequence stratigraphic analysis of Late Quaternary deposits of the Changhua Coastal Plain in the frontal arc-continent collision belt of Central Taiwan , 2010 .

[6]  D. Agnew,et al.  The complete (3‐D) surface displacement field in the epicentral area of the 1999 MW7.1 Hector Mine Earthquake, California, from space geodetic observations , 2001 .

[7]  L. Ge,et al.  Assessments of land subsidence in the Gippsland Basin of Australia using ALOS PALSAR data , 2015 .

[8]  B. Morris,et al.  The impact of urbanisation on groundwater quality (project summary report) , 1994 .

[9]  R. Nikolaidis.  Observation of geodetic and seismic deformation with the Global Positioning System , 2002 .

[10]  Kui Zhang,et al.  Mapping land subsidence in Jakarta, Indonesia using persistent scatterer interferometry (PSI) technique with ALOS PALSAR , 2012, Int. J. Appl. Earth Obs. Geoinformation.

[11]  V. B. H. Ketelaar,et al.  Satellite Radar Interferometry , 2009 .

[12]  Estimating coefficients of volume compressibility from compression of strata and piezometric changes in a multiaquifer system in west Taiwan , 2004 .

[13]  Manabu Watanabe,et al.  PALSAR Characterization and Initial Calibration , 2006, 2006 IEEE International Symposium on Geoscience and Remote Sensing.

[14]  Chris Rizos,et al.  Crustal deformation in Australia measured by satellite radar interferometry using ALOS/PALSAR imagery , 2009 .

[15]  Yuei-An Liou,et al.  Observing Land Subsidence and Revealing the Factors That Influence It Using a Multi-Sensor Approach in Yunlin County, Taiwan , 2015, Remote. Sens..

[16]  Timothy H. Dixon,et al.  Space geodetic imaging of rapid ground subsidence in Mexico City , 2008 .

[17]  Y. Hsu,et al.  Coseismic displacements and slip distribution from GPS and leveling observations for the 2006 Peinan earthquake (Mw 6.1) in southeastern Taiwan , 2009 .

[18]  P. Rosen,et al.  On the derivation of coseismic displacement fields using differential radar interferometry: The Landers earthquake , 1994, Proceedings of IGARSS '94 - 1994 IEEE International Geoscience and Remote Sensing Symposium.

[19]  Fabio Rocca,et al.  Permanent scatterers in SAR interferometry , 2001, IEEE Trans. Geosci. Remote. Sens..

[20]  Chris Rizos,et al.  Mapping accumulated mine subsidence using small stack of SAR differential interferograms in the Southern coalfield of New South Wales, Australia , 2010 .

[21]  L. Ge,et al.  Monitoring ground surface deformation over the North China Plain using coherent ALOS PALSAR differential interferograms , 2013, Journal of Geodesy.

[22]  K. Feigl,et al.  The displacement field of the Landers earthquake mapped by radar interferometry , 1993, Nature.

[23]  K. Feigl,et al.  Discrimination of geophysical phenomena in satellite radar interferograms , 1995 .

[24]  Chris Rizos,et al.  Mine subsidence monitoring using multi-source satellite SAR images , 2007 .

[25]  Xiaojing Li,et al.  Monitoring ground deformation in Beijing, China with persistent scatterer SAR interferometry , 2012, Journal of Geodesy.

[26]  Combination with precise leveling and PSInSAR observations to quantify pumping-induced land subsidence in central Taiwan , 2015 .

[27]  T. Wright,et al.  Toward mapping surface deformation in three dimensions using InSAR , 2004 .

[28]  V. Ketelaar.  Satellite Radar Interferometry: Subsidence Monitoring Techniques , 2009 .

[29]  D. Sandwell.  BIHARMONIC SPLINE INTERPOLATION OF GEOS-3 AND SEASAT ALTIMETER DATA , 1987 .

[30]  Jyr‐Ching Hu,et al.  Assessments of serious anthropogenic land subsidence in Yunlin County of central Taiwan from 1996 to 1999 by Persistent Scatterers InSAR , 2012 .

[31]  Chris Rizos,et al.  Multi-path PALSAR interferometric observations of the 2008 magnitude 8.0 Wenchuan Earthquake , 2010 .

[32]  Xiaojing Li,et al.  Estimating horizontal and vertical movements due to underground mining using ALOS PALSAR , 2012 .

[33]  J. Saastamoinen.  Contributions to the theory of atmospheric refraction , 1972 .

[34]  Y. Hsu,et al.  Correlation between groundwater level and altitude variations in land subsidence area of the Choshuichi Alluvial Fan, Taiwan , 2010 .

[35]  Chung-Pai Chang,et al.  Monitoring severe aquifer-system compaction and land subsidence in Taiwan using multiple sensors: Yunlin, the southern Choushui River Alluvial Fan , 2010 .

[36]  Chung-Pai Chang,et al.  Surface deformation from persistent scatterers SAR interferometry and fusion with leveling data: A case study over the Choushui River Alluvial Fan, Taiwan , 2011 .

[37]  Land subsidence using absolute and relative gravimetry: a case study in central Taiwan , 2010 .

[38]  M. Antonellini,et al.  Salt water intrusion in the coastal aquifer of the southern Po Plain, Italy , 2008 .

[39]  E. Chaussard,et al.  Potential and limits of InSAR to characterize interseismic deformation independently of GPS data: Application to the southern San Andreas Fault system , 2016 .

[40]  Claudio Prati,et al.  A New Algorithm for Processing Interferometric Data-Stacks: SqueeSAR , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[41]  Frederick R. Forst,et al.  On robust estimation of the location parameter , 1980 .

[42]  W. Hung,et al.  Modeling aquifer-system compaction and predicting land subsidence in central Taiwan , 2012 .

[43]  V. Kumar,et al.  Glacier surface velocity estimation using SAR interferometry technique applying ascending and descending passes in Himalayas , 2011, Int. J. Appl. Earth Obs. Geoinformation.

[44]  Jyh-Jong Liao,et al.  Characterization of land subsidence in the Choshui River alluvial fan, Taiwan , 2004 .

[45]  R. Goldstein,et al.  Satellite Radar Interferometry for Monitoring Ice Sheet Motion: Application to an Antarctic Ice Stream , 1993, Science.