Anthropogenic and natural ground deformation in the Hengill geothermal area, Iceland

We investigate crustal deformation due to the extraction of water and steam from a high-enthalpy geothermal reservoir; a common occurrence, yet not well understood. The cause of this deformation can be a change in pressure or in temperature in the reservoir, both of which can be caused by extraction or injection of geothermal fluids. Our study area, the Hengill mountains in SW Iceland, is an active volcanic center and a plate triple junction that hosts two power plants producing geothermal energy. This combination of natural and anthropogenic processes causes a complex displacement field at the surface. We analyze geodetic data—Global Navigation Satellite System and Interferometric Synthetic Aperture Radar—to obtain the surface velocity field, which we then simulate using an inverse modeling approach. We focus on the deformation around the geothermal power plants but need to model the regional tectonic and volcanic deformation as well, because the signals are overlapping. We find that plate motion and a deep contracting body can explain the broad scale signal in the area. Local deformation near the two power plants, Hellisheidi and Nesjavellir, can be explained by extraction of geothermal fluids. We estimate reservoirs extending from 0.6 to 3.0 km depth at Hellisheidi, and 1.0 to 3.0 km depth at Nesjavellir for observed pressure decrease rates of 0.25 MPa/yr and 0.1 MPa/yr, respectively. We find that the main cause for the subsidence in the geothermal area is the observed pressure drawdown.

[1]  P. Dobson,et al.  A 3D hydrogeological and geomechanical model of an Enhanced Geothermal System at The Geysers, California , 2014 .

[2]  Jonny Rutqvist,et al.  Monitoring deformation at the Geysers Geothermal Field, California using C‐band and X‐band interferometric synthetic aperture radar , 2013 .

[3]  F. Sigmundsson,et al.  The 1994–1995 seismicity and deformation at the Hengill triple junction, Iceland: Triggering of earthquakes by minor magma injection in a zone of horizontal shear stress , 1997 .

[4]  E. C. Robertson.  THERMAL PROPERTIES OF ROCKS , 1988 .

[5]  Andrew Hooper,et al.  A multi‐temporal InSAR method incorporating both persistent scatterer and small baseline approaches , 2008 .

[6]  Robert Tibshirani,et al.  Bootstrap Methods for Standard Errors, Confidence Intervals, and Other Measures of Statistical Accuracy , 1986 .

[7]  James L. Beck,et al.  Bayesian inversion for finite fault earthquake source models I—theory and algorithm , 2013 .

[8]  G. Gunnarsson,et al.  Model Simulations of the Hengill Area, Southwestern Iceland , 2011 .

[9]  P. LaFemina,et al.  A complex earthquake sequence captured by the continuous GPS network in SW Iceland , 2009 .

[10]  Sigurjón Jónsson,et al.  Plate boundary deformation and man-made subsidence around geothermal fields on the Reykjanes Peninsula, Iceland , 2010 .

[11]  Tim J. Wright,et al.  Reassessing the 2006 Guerrero slow‐slip event, Mexico: Implications for large earthquakes in the Guerrero Gap , 2015 .

[12]  Mark Simons,et al.  Deformation and seismicity in the Coso geothermal area, Inyo County, California: Observations and modeling using satellite radar interferometry , 2000 .

[13]  P. Einarsson.  Earthquakes and present-day tectonism in Iceland , 1991 .

[14]  S. Hurwitz,et al.  Numerical models of caldera deformation: Effects of multiphase and multicomponent hydrothermal fluid flow , 2009 .

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

[16]  T. Wright,et al.  Statistical comparison of InSAR tropospheric correction techniques , 2015 .

[17]  Freysteinn Sigmundsson,et al.  Fault slip distribution of two June 2000 M W 6.5 earthquakes in South Iceland estimated from joint inversion of InSAR and GPS measurements , 2003 .

[18]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[19]  Walter H. F. Smith,et al.  Free software helps map and display data , 1991 .

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

[21]  Richard M. Goldstein,et al.  Rapid subsidence over oil fields measured by SAR interferometry , 1998 .

[22]  Richard G. Gordon,et al.  Geologically current plate motions , 2010 .

[23]  Xiaoyong Zhan,et al.  Predicting subsidence at Wairakei and Ohaaki geothermal fields, New Zealand , 2000 .

[24]  R. Reilinger,et al.  Geodetic evidence for subsidence due to groundwater withdrawal in many parts of the United States of America , 1984 .

[25]  Jessica R. Murray,et al.  Modeling crustal deformation near active faults and volcanic centers: a catalog of deformation models and modeling approaches , 2013 .

[26]  Ari Tryggvason,et al.  Three-dimensional imaging of the P- and S-wave velocity structure and earthquake locations beneath Southwest Iceland , 2002 .

[27]  S. Jakobsdóttir.  Seismicity in Iceland: 1994–2007 , 2008, Jökull.

[28]  P. Einarsson.  Plate boundaries, rifts and transforms in Iceland , 2008, Jökull.

[29]  Freysteinn Sigmundsson,et al.  Surface effects of faulting and deformation resulting from magma accumulation at the Hengill triple junction, SW Iceland, 1994-1998 , 2002 .

[30]  Michael Heinert,et al.  Crustal deformation measured by GPS in the South Iceland Seismic Zone due to two large earthquakes in June 2000 , 2001 .

[31]  Freysteinn Sigmundsson,et al.  Crustal deformation near Hengill volcano, Iceland 1993–1998: Coupling between magmatic activity and faulting inferred from elastic modeling of satellite radar interferograms , 2000 .

[32]  Paul Segall,et al.  Subsidence at The Geysers Geothermal Field, N. California from a comparison of GPS and leveling surveys , 1997 .

[33]  Freysteinn Sigmundsson,et al.  The 2008 May 29 earthquake doublet in SW Iceland , 2010 .

[34]  Remko Scharroo,et al.  Generic Mapping Tools: Improved Version Released , 2013 .

[35]  J. Biggs,et al.  The influence of cooling, crystallisation and re-melting on the interpretation of geodetic signals in volcanic systems , 2014 .

[36]  Dante Fratta,et al.  Time-series analysis of surface deformation at Brady Hot Springs geothermal field (Nevada) using interferometric synthetic aperture radar , 2016 .

[37]  James H. Dieterich,et al.  Deformation from Inflation of a Dipping Finite Prolate Spheroid in an Elastic Half‐Space as a Model for Volcanic Stressing , 1988 .

[38]  V. Drouin.  Constraints on deformation processes in Iceland from space geodesy: seasonal load variations, plate spreading, volcanoes and geothermal fields , 2016 .

[39]  T. Herring,et al.  Introduction to GAMIT/GLOBK , 2006 .

[40]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[41]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[42]  J. D. Eshelby.  The determination of the elastic field of an ellipsoidal inclusion, and related problems , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[43]  S. Hickman,et al.  Subsidence rates at the southern Salton Sea consistent with reservoir depletion , 2016 .

[44]  S. Hurwitz,et al.  Hydrothermal fluid flow and deformation in large calderas: Inferences from numerical simulations , 2007 .

[45]  Rupture in the South Iceland Seismic Zone forced by magmatic intrusion in the Hengill area. , 2003 .

[46]  K. Mogi.  Relations between the Eruptions of Various Volcanoes and the Deformations of the Ground Surfaces around them , 1958 .

[47]  Wei Jiang,et al.  Glacial rebound and plate spreading: Results from the first countrywide GPS observations in Iceland , 2009 .

[48]  G. Axelsson,et al.  7.02 – The Physics of Geothermal Energy , 2012 .

[49]  H. Zebker,et al.  Fault Slip Distribution of the 1999 Mw 7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements , 2002 .

[50]  Paul Segall,et al.  Stress and subsidence resulting from subsurface fluid withdrawal in the epicentral region of the 1983 Coalinga earthquake , 1985 .