Proc. IAHS, 372, 335341, 2015 proc-iahs.net/372/335/2015/ doi:10.5194/piahs-372-335-2015 Author(s) 2015. CC Attribution 3.0 License.
O. Sarychikhina and E. Glowacka
Center for Scientic Research and Higher Education of Ensenada (CICESE), Ensenada, Mexico
Correspondence to: O. Sarychikhina ([email protected])
Published: 12 November 2015
Abstract. Ground deformation in Mexicali Valley, Baja California, Mexico, the southern part of the Mexicali-Imperial valley, is inuenced by active tectonics and human activity, mainly that of geothermal uid extraction in the Cerro Prieto Geothermal Field. Signicant ground deformation, mainly subsidence ( 18 cm yr1), and
related ground ssures cause severe damage to local infrastructure.
The technique of Differential Synthetic Aperture Radar Interferometry (DInSAR) has been demonstrated to be a very effective remote sensing tool for accurately measuring the spatial and temporal evolution of ground displacements over broad areas. In present study ERS-1/2 SAR and ENVISAT ASAR images acquired between 1993 and 2010 were used to perform a historical analysis of aseismic ground deformation in Mexicali Valley, in an attempt to evaluate its spatio-temporal evolution and improve the understanding of its dynamic. For this purpose, the conventional 2-pass DInSAR was used to generate interferograms which were used in stacking procedure to produce maps of annual aseismic ground deformation rates for different periods. Differential inter-ferograms that included strong co-seismic deformation signals were not included in the stacking and analysis.
The changes in the ground deformation pattern and rate were identied. The main changes occur between 2000 and 2005 and include increasing deformation rate in the recharge zone and decreasing deformation rate in the western part of the CPGF production zone. We suggested that these changes are mainly caused by production development in the Cerro Prieto Geothermal Field.
1 Introduction
Any process involving large-scale extraction of groundwater, including geothermal energy development, can be accompanied by ground deformation. Ground deformation related to geothermal energy production has been observed and studied in many geothermal elds all over the world (e.g. Mossop and Segall, 1997; Carnec and Fabriol, 1999; Fialko and Simons, 2000; Hole et al., 2007; Allis et al., 2009; Keiding et al., 2009). Ground deformation associated with geothermal energy production can be reduced or controled through monitoring combined with aquifer management (Lowe, 2012).
Differential synthetic aperture radar interferometry (DIn-
SAR) is a space-based geodetic technique widely applied for
large-scale monitoring (spatially continuous) of ground deformation at a low cost (Massonnet and Rabaute, 1993; Zebker et al., 1994; Massonnet and Feigl, 1998), and particularly for subsidence regardless of its cause (Massonnet et al., 1997; Fielding et al., 1998; Avallone et al., 1999; Carnec and Fabriol, 1999; Wright and Stow, 1999; Carnec and Dela-court, 2000). Compared with ground-based geodetic techniques, such as precise leveling and GPS surveys, which provide accurate and precise information but only at a number of measuring points in a deforming surface, DInSAR allows detecting, measuring and monitoring ground deformation over large areas ( 104 km2) with ne spatial resolution
( 102 m2) and high accuracy ( 1 cm) (Gabriel et al., 1989;
Brgmann et al., 2000; Hanssen, 2001). The systematic cov-
Published by Copernicus Publications on behalf of the International Association of Hydrological Sciences.
Preventionandmitigationofnaturalandanthropogenichazardsduetolandsubsidence
Open Access
Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley (Baja California, Mexico) from 1993 to 2010, using differential SAR interferometry
336 O. Sarychikhina and E. Glowacka: Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley
erage and frequent observations (e.g. approximately monthly repeat-cycle for ERS and ASAR) provided by SAR sensors allows the extraction of information on the spatial and temporal evolution of ground deformation at more frequent intervals than traditional techniques (Massonnet and Feigl, 1998;Baer et al., 2002; Toms et al., 2005; Herrera et al., 2007).The availability since 1992 of a large archive of satellite acquisitions allows an analysis of historical deformation.
In the present study ERS-1/2 SAR and ENVISAT ASAR images acquired between 1993 and 2010 were used to perform an analysis of historical ground deformation caused by geothermal energy development in the Cerro Prieto Geothermal Field (CPGF) in Mexicali Valley, in an attempt to evaluate its spatio-temporal evolution and improve the understanding of its dynamic.
2 Study area
Ground deformation in Mexicali Valley, Baja California, Mexico, the southern part of the Mexicali-Imperial Valley, is inuenced by active tectonics and human activity. Mexicali Valley is located in a complex tectonic environment, at the tectonic boundary between the Pacic and North American plates. This area is characterized by high tectonic seismicity, heat ow and surface deformation, related to the tec-tonic regime of the area. Besides the tectonic deformation, extraction of uids in the CPGF which is the oldest and largest Mexican geothermal eld in operation and the second largest in the world produces deformation of large magnitude (Glowacka et al., 1999). Signicant ground deformation (mainly subsidence) and related ground ssures cause severe damage to local infrastructure such as roads, irrigation canals and other facilities. A detailed plan-view map of the study area is presented in Fig. 1.
Several monitoring techniques have been employed in the study area including repeat ground-based geodetic surveys, geotechnical instruments measurements and DInSAR.Ground-based geodetic surveys (precise leveling and GPS) in Mexicali Valley began in the 1960s (Velasco, 1963) and have been carried out up to the present, with varying frequency, precision, coverage, and density (Garca, 1978; Grannell et al., 1979; de la Pea, 1981; Wyman, 1983; Lira and Arellano, 1997; Glowacka et al., 1999, 2012). Since 1996, the geotechnical instruments network has operated in Mexicali Valley (REDECVAM) for continuous recording of deformation phenomena (Nava and Glowacka, 1999; Glowacka et al., 2002, 2010a, b; Sarychikhina et al., 2015).
Analysis of leveling data and CPGF extraction history (Glowacka et al., 1999), and modelling of the tectonic and anthropogenic components of subsidence (Glowacka et al., 2005) suggested that the current deformation rate is mainly related to the uid extraction, while the eld observations (Gonzalez et al., 1998; Glowacka et al., 1999, 2010a, c;Surez-Vidal et al., 2008) and geotechnical instruments data
Figure 1. Plan-view map of the study area. The black thick dotted line frames the limits of the Cerro Prieto Geothermal Field. The borders of the evaporation pond and the Cerro Prieto volcano (CPV) are also shown (thin gray lines). A striped orange rectangle indicates the extraction area of CPIV which started the operation since 2000. Solid red lines are well-known surface traces of tectonic faults. CPF = Cerro Prieto Fault, GF = Guerrero Fault,
IF = Imperial Fault, MF = Morelia Fault, SF = Saltillo fault. Dotted
red lines are proposed surface fault traces based on mapped ssure zones from Gonzlez et al. (1998), Glowacka et al. (2006, 2010a, c), Lira (2006) and Surez-Vidal et al. (2007, 2008). SF is continuation of Saltillo fault as proposed by Surez-Vidal et al. (2008). Black square marks the location of the reference, considered stable, region. Black thick line corresponds to the prole A-A illustrated in Fig. 3. Modied from Glowacka et al. (2010c).
(Glowacka et al., 1999, 2010a, c) reveal that the geometry of the subsiding area is controlled by tectonic faults.
Ground deformation in the Mexicali Valley has also been measured and monitored using the DInSAR technique (Carnec and Fabriol, 1999; Hanssen, 2001; Sarychikhina et al., 2011, 2015). Carnec and Fabriol (1999) and Hanssen (2001) used ERS1/2 images acquired during 1993 97 and 199597, respectively. Signicant local subsidence has been detected. Due to the high degree of temporal and spatial decorrelation, reliable data was obtained only over the mainly desert area of CPGF. Although agricultural activity in the area limited the investigation, interferometric monitoring revealed that the ground deformation is associated with the geothermal uid withdrawal and was in agreement with the leveling data. More recent studies by Sarychikhina et al. (2011, 2014, 2015) and Sarychikhina and Glowacka (2015) used Envisat ASAR images acquired during 200310 to evaluate the spatial pattern and rate of land subsidence. The Differential Interferograms Stacking (DIS) procedure was applied to estimate the annual ground deformation rate for different periods. The results were compared to 199497 and 19972006 leveling data, revealing signi-cant changes in ground deformation pattern and rate.
Proc. IAHS, 372, 335341, 2015 proc-iahs.net/372/335/2015/
ing and descending tracks (for 2007), along with the known vertical displacement from ground-based measurements including leveling surveys, suggest that the observed LOS displacements may be interpreted as reecting mostly vertical surface displacement. However, Fig. 3, which presents the annual LOS displacement rates along the prole (A-A in Fig. 1) across the subsidence zone, shows important differences ( 2 cm of magnitude) between ascending and de
scending datasets in the central part of the CPGF. This indicates the inuence of a horizontal displacement component, probably related to horizontal displacement at and beyond the margin of the subsidence bowl, as proposed by Segall (1989).
The negative rate values in Fig. 2 indicate an increase in the LOS range that corresponds to ground subsidence. The subsiding area boundaries appear to correlate with faults, as can be seen in Fig. 2 where they are superimposed onto the maps of LOS displacement rate. The rst region of subsidence, with larger negative-valued LOS displacement rate for all analysed periods (Fig. 3), is located in the area between the eastern limits of the CPGF and the Saltillo fault, which was proposed as a recharge zone in previous studies (Glowacka et al., 1999, 2005; Sarychikhina, 2003). The maximum LOS displacement rate in this region increased from about 12.5 to 17 cm yr1 during 199397, and
200507 and 2009, respectively. The largest negative-valued LOS displacement rate (17 cm yr1) corresponds to about
18 cm yr1 of ground subsidence. The area of this subsidence bowl is also evidently increasing over time. This increasing
LOS displacement (subsidence) rate and areal extend of subsidence is probably related to geothermal uid production causing continued uid-pressure declines in the area of extraction in the CPGF that extend to the recharge zone. The 199497 leveling survey data zone (Glowacka et al., 1999, 2005; Sarychikhina et al., 2011) projected to the LOS direction shows a maximum LOS displacement rate of about
8 cm yr1 in the recharge which is not in accord with our results. This discordance is probably due to limited coverage of the leveling network.
The second region of high subsidence rate is located in the CPGF production zone. Because the ascending and descending pass datasets present important differences in this zone, we compared the data only from the same pass.
For the descending pass data, the maximum LOS displacement rate occurring in the eastern part of the CPGF production zone is about 10 cm yr1 for all compared periods
(199397, 19982000, 200507). In the western and central parts of the CPGF production zone the LOS displacement rate is lower in 200507 than in 199397 and 19982000. It could be related to decreased production in the western part of CPGF (CP I). The LOS displacement rates for individual years during 200507 period are similar.
For ascending pass data, we found important differences in LOS displacement rate and pattern between the 2004 and the 2007, 2009 periods. However, due to the fact that the 2004
proc-iahs.net/372/335/2015/ Proc. IAHS, 372, 335341, 2015
O. Sarychikhina and E. Glowacka: Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley 337
Here we extend these studies by using aditional SAR images (ERS 1/2 SAR) and the DIS procedure to obtain ground deformation rate for different periods from 1993 to 2009.
3 Methodology for deformation rate estimation
The conventional 2-pass DInSAR method using ERS 1/2 SAR and ENVISAT ASAR images was used to produce interferograms. The DIS procedure was used to estimate ground deformation rates in the study area for different time periods. Here, the DIS procedure, which consisted in averaging of multiple differential interferograms, has been applied in order to overcome limitations of conventional DInSAR, such as low coherence over long temporal separations and phase distortions introduced by atmospheric effects.
A total of 27 ERS 1/2 SAR and 52 ENVIAST ASAR images were acquired from the European Space Agency (ESA), as part of ESA CAT-1 project (ID C1P3508). ERS 1/2 SAR images were from descending pass (track 84, frame 2961) and covered the period 19932000. ENVISAT ASAR images were from ascending (track 306, frame 639) and descending (track 84, frame 2961) passes and covered the period 200309.
Data processing was performed using the GAMMA
commercial software, developed by Wegmller and Werner (1997) by rst performing two-pass interferometry and unwrapping of differential interferograms, and then the stacking of multiple differential interferograms and geocoding. The selection of interferograms for stacking was based on their quality in terms of high coherence and low noise level.
Unfortunately, it was impossible to form the annual stacks of coherent short temporal baseline (35140 days) differential interferograms with reasonable coherence level from ERS 1/2 dataset because of important gaps in data acquisition and geometrical decorrelation of some short time period differential interferograms. For this reason, the pairs from ERS 1/2 were grouped into two time periods (stacks): 1993 97 and 19982000.
Because here we attempt to present results of only aseismic deformation the differential interferograms that presented strong co-seismic deformation signals were not included in the stacking and analysis. It was impossible to generate an aseismic LOS displacement rate map for 2008 due to a high level of moderate sized seismicity during this year.
4 Historical deformation analysis
Figure 2 illustrates the maps of annual LOS displacement rate for different time periods obtained using the DIS method. All maps of LOS displacement rate show a roughly NESW oriented elliptical-shaped pattern of deformation with two bowls exhibiting high deformation rates. The similarity between maps of LOS displacement rate from ascend-
338 O. Sarychikhina and E. Glowacka: Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley
Figure 2. Maps of LOS displacement rate (cm yr1) showing the time evolution of the ground deformation eld in and around the Cerro Prieto Geothermal Field in Mexicali Valley. Maps were obtained using the multiple Differential Interferograms Stacking (DIS) procedure.
The numbers in the right lower corner indicate the number of differential interferograms used in stack. Negative displacement rates indicate an increase in the LOS range (ground subsidence). Black square marks the location of the reference, considered stable, region. Tectonic faults (solid and dotted red lines), CPGF limits (black thick dotted line), evaporation pond and CPV (thin gray lines) are superimposed onto the maps. The maps are in UTM coordinates, the marks are every 5 km.
LOS displacement rate was estimated using only 2 differential interferograms, we think that this difference is due to possible errors (e.g. atmospheric artifacts) in the 2004 dataset.There is no signicant difference between the 2007 and 2009LOS displacement rates.
5 Conclusions
A comparison of the annual aseismic LOS displacement rates for 199397, 19982000, 200407, and 2009 periods obtained using multiple Differential Interferograms Stacking (DIS) method reveals the important changes in spatial pattern and rate of ground deformation in the study area.The main changes occur between 2000 and 2005 and in-
Proc. IAHS, 372, 335341, 2015 proc-iahs.net/372/335/2015/
O. Sarychikhina and E. Glowacka: Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley 339
Figure 3. Comparison between LOS displacement rates derived from multiple Differential Interferograms Stacking (DIS) procedure for different periods along the prole A-A (location shown in Fig. 1). Solid red lines indicate location of tectonic faults which cross the prole. Faults notation as in Fig. 1.
clude increasing deformation rate in the recharge zone and decreasing deformation rate in the western part of the CPGF production zone. The observed changes may be caused by changes in production in the CPGF, an increased production in the eastern part of the eld due to the newest power plant (CP IV), which started operating in 2000, and decreased production in the western part of CPGF (CP I). Since 2000 there was not considerable changes in the total CPGF production, so based on the results of this study, we can state that it took a maximum 5 years for ground deformation to stabilize after the perturbation caused by changes in the CPGF production in 2000.
During 19932000 and 200509 periods, the spatial pattern and rate of deformation could be considered constant (or without signicant changes).
Acknowledgements. Data from the European Space Agencys ENVISAT satellite was used to generate the interferometric data. The data were obtained as part of ESA Cat-1 Project (ID C1P3508).
This research was sponsored in part by CONACYT, project 105907 and CICESE internal funds.
References
Allis, R. G., Bromley, C. J., and Currie, S.: Update on subsidence at the WairakeiTauhara Geothermal System, Geothermics, 38, 169180, 2009.
Avallone, A., Briole, P., Delacourt, C., Zollo, A., and Beauducel, F.: Subsidence at Campi Flegrei (Italy) detected by SAR interferometry, Geophys. Res. Lett., 26, 23032306, 1999.
Baer, G., Schattner, U., Wachs,D., Sandwell, D., Wdowinski, S., and Frydman, S.: The lowest place on Earth is subsiding An In-SAR (interferometric synthetic aperture radar) perspective, GSA Bulletin, 114, 1223, 2002.
Brgmann, R., Rosen, P. A., and Fielding, E. J.: Synthetic aperture radar interferometry to measure Earths surface topography and its deformation, Ann. Rev. Earth Planet. Sci., 28, 169209, 2000.Carnec, C. and Delacourt, C.: Three years of mining subsidence monitored by SAR interferometry, near Gardanne, France, J.Appl. Geophys., 43, 4354, 2000.
Carnec, C. and Fabriol, H.: Monitoring and modeling land subsidence at the Cerro Prieto geothermal eld, Baja California, Mexico, using SAR interferometry, Geophys. Res. Lett., 26, 1211 1214, 1999.de la Pea, L. A.: Results from the rst order leveling surveys carried out in the Mexicali Valley and at the Cerro Prieto eld, Baja California, Proceedings, in: Proceedings of the 3rd Symposium on the Cerro Prieto Geothermal Field, Baja California, Mexico, 2426 March 1981, San Francisco, California, USA, 281291, 1981.
Fialko, Y. and Simons, M.: Deformation and seismicity in the Coso geothermal area, Inyo County, California: Observations and modeling using satellite radar interferometry, J. Geophys.Res., 105, 2178121794, 2000.
Fielding, E. J., Blom, R. G., and Goldstein, R. M.: Rapid subsidence over oil elds measured by SAR interferometry, Geophys. Res.Lett., 27, 32153218, 1998.
Gabriel, A. G., Goldstein, R. M., and Zebker, H. A.: Mapping small elevation changes over large areas: Differential radar interferometry, J. Geophys. Res., 94, 91839191, 1989.
Garca, J. R.: Estudios de nivelacin de primer orden en Cerro Prieto, in: Proceedings of the 1st Symposium on the Cerro Prieto Geothermal Field, Baja California, Mxico, 2022 September 1978, San Diego, California, USA, 148150, 1978.Glowacka, E., Gonzlez, J., and Fabriol, H.: Recent vertical deformation in Mexicali Valley and its relationship with tectonics, seismicity, and the exploitation of the Cerro Prieto geothermal eld, Mexico, Pure Appl. Geophys., 156, 591614, 1999.Glowacka, E., Nava, F. A., de Cossio, G. D., Wong, V., and Farfan, F.: Fault slip, seismicity, and deformation in Mexicali Valley, Baja California, Mexico, after the M 7.1 1999 Hector Mine earthquake, Bull. Seismol. Soc. Am., 92, 12901299, 2002.Glowacka, E., Sarychikhina, O., and Nava, F. A.: Subsidence and stress change in the Cerro Prieto Geothermal Field, B.C., Mexico, Pure Appl. Geophys., 162, 20952110, 2005.
Glowacka, E., Sarychikhina, O., Surez, F., Mendoza, R., and Nava,F. A.: Estudio geolgico para denir la zona de hundimiento con el n de relocalizacin del canal Nuevo Delta en el Valle de Mexicali, Informe Tcnico, CICESE, Mxico, 2006.
Glowacka, E., Sarychikhina, O., Nava, F. A., Suarez, F., Ramrez, J., Guzman, M., Robles, B., Farfan, F., and De Cossio Batani, G. D.:
proc-iahs.net/372/335/2015/ Proc. IAHS, 372, 335341, 2015
340 O. Sarychikhina and E. Glowacka: Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley
Continuous monitoring techniques of fault displacement caused by geothermal uid extraction in the Cerro Prieto Geothermal Field (Baja California, Mexico), in: Land Subsidence, Associated Hazards and the Role of Natural Resources, edited by: Corren-Freyre D., Cerca, M., and Galowey, D., IAHS Publ., 339, 326332, 2010a.
Glowacka, E., Sarychikhina, O., Surez, F., Nava, F. A., and Mellors, R.: Anthropogenic subsidence in Mexicali Valley, Baja California, Mexico, and slip on the Saltillo fault, Environ. Earth Sci., 59, 15151524, 2010b.
Glowacka, E., Sarychikhina, O., Surez, F., Nava, F. A., Farfan,F., De Cossio Batani, G. D., and Garcia Arthur, M. A.: Anthropogenic subsidence in the Mexicali Valley, B.C., Mexico, caused by the uid extraction in the Cerro Prieto geothermal Field and the role of faults, in: Proceedings of the World Geothermal Congress (WGC) 2010, 2529 April 2010, Bali, Indonesia, 2010c.
Glowacka, E., Sarychikhina, O., Robles, B., Surez, F., Ramrez,J., and Nava, F. A.: Estudio geolgico para denir la lnea de hundimiento cero y monitorear la subsidencia de los mdulos 10, 11 y 12 en el Valle de Mexicali, en el distrito de riego 014, Rio Colorado, B.C., Reporte Tcnico Final, Convenio con CONAGUA, CICESE, Mexico, 2012.
Gonzlez, J., Glowacka, E., Surez, F., Quiones, J. G., Guzmn,M., Castro, J. M., Rivera, F., and Flix, M. G.: Movimiento reciente de la Falla Imperial, Mexicali, B. C, Ciencia para todos Divulgare, Universidad Autnoma de Baja California, 6, 415, 1998.
Grannell, R. B., Tarnman, D. W., Clover, R. C., Leggewie, R. M.,
Aronstam, P. S., Kroll, R. C., and Eppink, J.: Precision gravity studies at Cerro Prieto, Proceedings, Second Symposium on the Cerro Prieto, in: Proceedings of the 2nd Symposium on the Cerro Prieto Geothermal Field, Baja California, Mexico, 1719 October 1979, Mexicali, Baja California, Mexico, 329331, 1979.Hanssen, R. F.: Radar Interferometry: Data Interpretation and Error
Analysis, Kluwer Academic Publishers, Dordrecht, the Netherlands, 2001.
Herrera, G., Toms, R., Lopez-Sanchez, J. M., Delgado, J., Mallorqui, J. J., Duque, S., and Mulas, J.: Advanced DInSAR analysis on mining areas: La Union case study (Murcia, SE Spain), Eng. Geol., 90, 148159, 2007.
Hole, J., Bromley, C., Stevens, N., and Wadge, G.: Subsidence in the geothermal elds of the Taupo Volcanic Zone, New Zealand from 1996 to 2005 measured by InSAR, J. Volcanol. Geoth. Res., 166, 125146, 2007.
Keiding, M., Hooper, A., rnadttir, T., Jnsson, S., and Decriem,J.: Natural and man-made deformation around geothermal elds on the Reykjanes peninsula, SW Iceland, in: Proceedings, of the Workshop Fringe 2009, 30 November4 December 2009, ESA, Frascati, Italy, 4 pp., 2009.
Lira, H.: Caractersticas del sismo del 23 de Mayo de 2006, Informe RE-023/2006, Comisin Federal de Electricidad, Residencia de Estudios, Mxico, 2006.
Lira, H. and Arellano, J. F.: Resultados de la nivelacin de precisin realizada en 1997, en el campo geotrmico Cerro Prieto, Informe Tcnico RE 07/97, Comisin Federal de Electricidad, Residencia de Estudios, Mxico, 1997.
Lowe, M.: Subsidence in sedimentary basins due to groundwater withdrawal for geothermal energy development, Utah Geo-
logical Survey, Utah Department of Natural Resources, Open File Report 601, 13 pp., available at: http://files.geology.utah.gov/online/ofr/ofr-601.pdf
Web End =http://les.geology.utah. http://files.geology.utah.gov/online/ofr/ofr-601.pdf
Web End =gov/online/ofr/ofr-601.pdf (last access: 12 October 2015), 2012. Massonnet, D. and Feigl, K. L.: Radar interferometry and its application to changes in the Earths surface, Rev. Geophys., 36, 441500, 1998.
Massonnet, D. and Rabaute, T.: Radar interferometry: limits and potential, IEEE T. Geosci. Remote, 31, 455464, 1993. Massonnet, D., Holzer, T., and Vadon, H.: Land subsidence caused by the East Mesa Geothermal eld, California, observed using SAR interferometry, Geophys. Res. Lett., 24, 901904, 1997. Mossop, A. and Segall, P.: Subsidence at the Geysers geothermal eld, N. California from a comparison of GPS and leveling surveys, Geophys. Res. Lett., 24, 18391842, 1997.
Nava, F. A. and Glowacka, E.: Fault slip triggering, healing, and viscoelastic after working in sediments in the Mexicali-Imperial Valley, Pure Appl. Geophys., 156, 615629, 1999. Sarychikhina, O.: Modelacin de subsidencia en el campo geotrmico Cerro Prieto, MS thesis, CICESE, Mexico, 101 pp., 2003. Sarychikhina, O. and Glowacka, E.: Application of DInSAR Stacking Method for Monitoring of Surface Deformation Due to Geothermal Fluids Extraction in the Cerro Prieto Geothermal Field, Baja California, Mexico, in: Proceedings of the World Geothermal Congress (WGC) 2015, 1925 April 2015, Melbourne, Australia, 2015.
Sarychikhina, O., Glowacka, E., Mellors, R., and Surez-Vidal, F.: Land subsidence in the Cerro Prieto Geothermal Field, Baja California, Mexico, from 1994 to 2005. An integrated analysis of DInSAR, leveling and geological data, J. Volcanol. Geoth. Res., 204, 7690, 2011.
Sarychikhina, O., Glowacka, E., Robles, B., and Mojarro, J.: Spatio-temporal evolution of anthropogenic deformation around Cerro Prieto geothermal eld in the Mexicali Valley, B.C., Mexico, between 1993 and 2009 from DInSAR and leveling, EGU General Assembly 2014, EGU2014-14709, 27 April2 May 2014, Vienna, Austria, 2014.
Sarychikhina, O., Glowacka, E., Robles, B., Nava, F. A., and Guzmn, M.: Estimation of seismic and aseismic deformation in Mexicali Valley, Baja California, Mexico, in the 2006-2009 period, using precise leveling, DInSAR, geotechnical instruments data, and modeling, Pure Appl. Geophys., doi:http://dx.doi.org/10.1007/s00024-015-1067-0
Web End =10.1007/s00024- http://dx.doi.org/10.1007/s00024-015-1067-0
Web End =015-1067-0 , in press, 2015.
Segall, P.: Earthquakes triggered by uid extraction, Geology, 17,942946, 1989.
Surez-Vidal, F., Mungua-Orozco, L., Gonzlez-Escobar, M.,
Gonzlez-Garca, J., and Glowacka, E.: Surface rupture of the Morelia fault near the Cerro Prieto Geothermal Field, Mexicali, Baja California, Mexico, during the Mw 5.4 earthquake of24 May 2006, Seismol. Res. Lett., 78, 394399, 2007. Surez-Vidal, F., Mendoza-Borunda, R., Nafarrete-Zamarripa, L.,
Rmirez, J., and Glowacka, E.: Shape and dimensions of the Cerro Prieto pull-apart basin, Mexicali, Baja California, Mxico, based on the regional seismic record and surface structures, Int. Geol. Rev., 50, 636649, 2008.
Toms, R., Mrquez, Y., Lopez-Sanches, J. M., Delgado, J., Blanco,P., Mallorqu, J. J., Martnez, M., Herrera, G., and Mulas, J.: Mapping ground subsidence induced by aquifer overexploitation using advanced Differential SAR Interferometry: Vega Media of
Proc. IAHS, 372, 335341, 2015 proc-iahs.net/372/335/2015/
O. Sarychikhina and E. Glowacka: Spatio-temporal evolution of aseismic ground deformation in the Mexicali Valley 341
the Segura River (SE Spain) case study, Remote Sens. Environ., 98, 269283, 2005.
Velasco, J.: Levantamiento Gravimtrico, Zona Geotrmica de Mexicali, Baja California, Internal report, Consejo de Recursos Naturales no Renovables, Mexico, 1963.
Wegmller, U. and Werner, C. L.: Gamma SAR processor and interferometry software, in: Proceedings of the 3rd ERS Scientic Symposium, 1421 March 1997, Florence, Italy, 16861692, 1997.
Wright, P. A. and Stow, R. J.: Detecting mining subsidence from space, Int. J. Remote Sens., 20, 11831188, 1999.
Wyman, R. M.: Potential Modeling of Gravity and Leveling Data over Cerro Prieto Geothermal Field, MS thesis, California State University, Long Beach, CA, USA, 1983.
Zebker, H. A., Rosen, P. A., Goldstein, R. M., Gabriel, A., and Werner, C. L.: On the derivation of coseismic displacement elds using differential radar interferometry: The Landers earthquake,J. Geophys. Res., 99, 1961719634, 1994.
proc-iahs.net/372/335/2015/ Proc. IAHS, 372, 335341, 2015
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright Copernicus GmbH 2015
Abstract
Ground deformation in Mexicali Valley, Baja California, Mexico, the southern part of the Mexicali-Imperial valley, is influenced by active tectonics and human activity, mainly that of geothermal fluid extraction in the Cerro Prieto Geothermal Field. Significant ground deformation, mainly subsidence (~ 18 cm yr<sup>-1</sup>), and related ground fissures cause severe damage to local infrastructure. <br><br> The technique of Differential Synthetic Aperture Radar Interferometry (DInSAR) has been demonstrated to be a very effective remote sensing tool for accurately measuring the spatial and temporal evolution of ground displacements over broad areas. In present study ERS-1/2 SAR and ENVISAT ASAR images acquired between 1993 and 2010 were used to perform a historical analysis of aseismic ground deformation in Mexicali Valley, in an attempt to evaluate its spatio-temporal evolution and improve the understanding of its dynamic. For this purpose, the conventional 2-pass DInSAR was used to generate interferograms which were used in stacking procedure to produce maps of annual aseismic ground deformation rates for different periods. Differential interferograms that included strong co-seismic deformation signals were not included in the stacking and analysis. <br><br> The changes in the ground deformation pattern and rate were identified. The main changes occur between 2000 and 2005 and include increasing deformation rate in the recharge zone and decreasing deformation rate in the western part of the CPGF production zone. We suggested that these changes are mainly caused by production development in the Cerro Prieto Geothermal Field.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer