Estimation of small surface displacements in the Upper Rhine Graben area from a combined analysis of PS-InSAR, levelling and GNSS data

The intra-plate deformation of the Upper Rhine Graben (URG) located in Central Europe is investigated using geodetic measurement techniques. We present a new approach to calculate a combined velocity field from InSAR, levelling and GNSS measurements. As the expected tectonic movements in the URG area are small (less than 1mma−1), the best possible solutions for linear velocity rates from single-technique analyses are estimated in a first step. Second, we combine the velocity rates obtained from InSAR (line of sight velocity rates in ascending and descending image geometries), levelling (vertical velocity rates) and GNSS (horizontal velocity rates) using least-squares adjustment (LSA). Focusing on the Northern URG area, we analyse SAR data on four different image stacks (ERS ascending, ERS descending, Envisat ascending, Envisat descending) using the Persistent Scatterer (PS) approach. The linear velocity rates in ascending and descending image geometries, respectively, are estimated in an LSA from joint time-series analysis of ERS and Envisat data. Vertical velocity rates from levelling are obtained from a consistent adjustment of more than 40 000 measured height differences using a kinematic displacement model. Horizontal velocity rates in east and north direction are calculated from a time-series analysis of daily coordinate estimates at 76 permanently operating GNSS sites in the URG region. As the locations, at which the measurement data of PS-InSAR, levelling and GNSS reside, do not coincide, spatial interpolation is needed during several steps of the rigorous processing. We use Ordinary Kriging to interpolate from a given set of data points to the locations of interest with a special focus on the modeling and propagation of errors. The final 3-D velocity field is calculated at a 200 m grid, which carries values only close to the location of PS points, resulting in a mean horizontal and vertical precision of 0.30 and 0.13 mm a−1, respectively. The vertical component of the combined velocity field shows a significant subsidence of about 0.5 mm a−1 in the northern part of the graben coinciding with a well-known quaternary basin structure. Horizontal displacement rates of up to 0.8 mm a−1 in southeast direction are observed outside the graben, in reasonable alignment with the average direction of maximum horizontal stress. Within the graben, the velocity directions rotate toward east in the non-subsiding part, while an opposite trend is observed in the subsiding part of the graben. The complexities of the observed velocity field are compatible to the geomechanical situation in our investigation area which is characterized by a transition from a restraining to a releasing bend setting. Glacial isostatic adjustment is another potential source influencing the observed velocity field, as well as anthropogenic signals due to mining, oil exploration and groundwater usage that have been identified in some places.

[1]  Z. Altamimi,et al.  ITRF2005 : A new release of the International Terrestrial Reference Frame based on time series of station positions and Earth Orientation Parameters , 2007 .

[2]  Recomputation of vertical crustal motions near Palmdale, California, 1959–1975 , 1983 .

[3]  F. Cotton,et al.  GPS network monitors the Western Alps' deformation over a five-year period: 1993–1998 , 2002 .

[4]  Hermann Baehr,et al.  Orbital Effects in Spaceborne Synthetic Aperture Radar Interferometry , 2013 .

[5]  F. V. Leijen,et al.  Persistent Scatterer Interferometry based on geodetic estimation theory , 2014 .

[6]  J. Catalão,et al.  INTEGRATION OF INSAR AND GPS FOR VERTICAL DEFORMATION MONITORING : A CASE STUDY IN FAIAL AND PICO ISLANDS , 2010 .

[7]  T. Buchmann,et al.  Contemporary kinematics of the Upper Rhine Graben: A 3D finite element approach , 2007 .

[8]  A. Schenk,et al.  Towards A Fusion Of SAR-Interferometry, GNSS And Precise Levelling In The Upper Rhine Graben Area, Southwest Germany , 2013 .

[9]  E. Calais,et al.  Geodetic Measurements of Crustal Deformation in the Western Mediterranean and Europe , 2004 .

[10]  R. Hanssen,et al.  Radar Interferometry with Public Domain Tools , 2004 .

[11]  T. Fuhrmann,et al.  Vertical displacement rates in the Upper Rhine Graben area derived from precise leveling , 2014, Journal of Geodesy.

[12]  Ramon F. Hanssen,et al.  Radar target type classification and validation , 2013, 2013 IEEE International Geoscience and Remote Sensing Symposium - IGARSS.

[13]  S. Schmid,et al.  Evolution of the European Cenozoic Rift System: interaction of the Alpine and Pyrenean orogens with their foreland lithosphere , 2004 .

[14]  A. Hooper,et al.  Recent advances in SAR interferometry time series analysis for measuring crustal deformation , 2012 .

[15]  W. Peltier,et al.  Constraining Models of Postglacial Rebound Using Space Geodesy , 2008 .

[16]  B. Kampes Displacement parameter estimation using permanent scatterer interferometry , 2005 .

[17]  G. Hein,et al.  Height Changes in the Rhenish Massif: Determination and Analysis , 1983 .

[18]  Xavier Collilieux,et al.  ITRF2008 plate motion model , 2011 .

[19]  R. Egli,et al.  A modified least-squares collocation method for the determination of crustal deformation: first results in the Swiss Alps , 2007 .

[20]  P. Cross,et al.  Geodetic determination of tectonic deformation in central Greece from 1900 to 1988 , 1991, Nature.

[21]  D. Giardini,et al.  Active Normal Faulting in the Upper Rhine Graben and Paleoseismic Identification of the 1356 Basel Earthquake , 2001, Science.

[22]  J. Ritter,et al.  The 2005 Heidelberg and Speyer earthquakes and their relationship to active tectonics in the central Upper Rhine Graben , 2009 .

[23]  R. Goldstein,et al.  Mapping small elevation changes over large areas: Differential radar interferometry , 1989 .

[24]  T. Wright,et al.  Constraining crustal velocity fields with InSAR for Eastern Turkey: Limits to the block‐like behavior of Eastern Anatolia , 2014 .

[25]  Hans-Georg Scherneck,et al.  Recent results based on continuous GPS observations of the GIA process in Fennoscandia from BIFROST , 2010 .

[26]  Simon D. P. Williams,et al.  CATS: GPS coordinate time series analysis software , 2008 .

[27]  Ramon F. Hanssen,et al.  Target Characterization and Interpretation of Deformation Using Persistent Scatterer Interferometry and Polarimetry , 2011 .

[28]  J. Nocquet Present-day kinematics of the Mediterranean: A comprehensive overview of GPS results , 2012 .

[29]  C. C. Wang,et al.  Robust Estimating Three-Dimensional Ground Motions from Fusion of InSAR and GPS Measurements , 2011, 2011 International Symposium on Image and Data Fusion.

[30]  K. Bonjer,et al.  The stress field in the Rhine Graben area inferred from earthquake focal mechanisms and estimation of frictional parameters , 1997 .

[31]  J. Nocquet,et al.  Crustal velocity field of western Europe from permanent GPS array solutions, 1996–2001 , 2003 .

[32]  D. Hayashi,et al.  Finite element modelling of the pull-apart formation: implication for tectonics of Bengo Co pull-apart basin, southern Tibet , 2010 .

[33]  Peter Teunissen,et al.  Least-squares prediction in linear models with integer unknowns , 2007 .

[34]  A. Demoulin,et al.  Recent crustal movements in the southern Black Forest (western Germany) , 1998 .

[35]  G. Peters Active tectonics in the Upper Rhine Graben Integration of paleoseismology, geomorphology and geomechanical modeling , 2007 .

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

[37]  Enrico Serpelloni,et al.  Vertical GPS ground motion rates in the Euro‐Mediterranean region: New evidence of velocity gradients at different spatial scales along the Nubia‐Eurasia plate boundary , 2013 .

[38]  R. L. Hardy,et al.  Solvability and multiquadric analysis as applied to investigations of vertical crustal movements , 1979 .

[39]  P. A. Ziegler European Cenozoic rift system , 1992 .

[40]  Ramon Hanssen,et al.  Surface Deformation Of The Whole Netherlands After PSI Analysis , 2012 .

[41]  D. Argus,et al.  Tests of the rigid-plate hypothesis and bounds on intraplate deformation using geodetic data from very long baseline interferometry , 1996 .

[42]  Helmut Moritz,et al.  Least‐squares collocation , 1978 .

[43]  C. Deutsch,et al.  Teacher's Aide Variogram Interpretation and Modeling , 2001 .

[44]  J. Behrmann,et al.  Anatomy and kinematics of oblique continental rifting revealed: A three-dimensional case study of the southeast Upper Rhine graben (Germany) , 2003 .

[45]  H. Zebker,et al.  Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcán Alcedo, Galápagos , 2007 .

[46]  Jean-Mathieu Nocquet,et al.  Geodetic constraints on glacial isostatic adjustment in Europe , 2005 .

[47]  Karl Zippelt,et al.  Recent surface displacements in the Upper Rhine Graben — Preliminary results from geodetic networks , 2013 .

[48]  B. Hoffers,et al.  Stress pattern and strain release in the Alpine foreland , 1981 .

[49]  R. Hanssen Radar Interferometry: Data Interpretation and Error Analysis , 2001 .

[50]  W. Peltier,et al.  Space geodesy constrains ice age terminal deglaciation: The global ICE‐6G_C (VM5a) model , 2015 .

[51]  K. McClay,et al.  4D analogue modelling of transtensional pull-apart basins , 2009 .