Subsidence activity maps derived from DInSAR data: Orihuela case study

A new methodology is proposed to produce subsidence activity maps based on the geostatistical analysis of persistent scatterer interferometry (PSI) data. PSI displacement measurements are interpolated based on conditional Sequential Gaussian Simulation (SGS) to calculate multiple equiprobable realizations of subsidence. The result from this process is a series of interpolated subsidence values, with an estimation of the spatial variability and a confidence level on the interpolation. These maps complement the PSI displacement map, improving the identification of wide subsiding areas at a regional scale. At a local scale, they can be used to identify buildings susceptible to suffer subsidence related damages. In order to do so, it is necessary to calculate the maximum differential settlement and the maximum angular distortion for each building of the study area. Based on PSI-derived parameters those buildings in which the serviceability limit state has been exceeded, and where in situ forensic analysis should be made, can be automatically identified. This methodology has been tested in the city of Orihuela (SE Spain) for the study of historical buildings damaged during the last two decades by subsidence due to aquifer overexploitation. The qualitative evaluation of the results from the methodology carried out in buildings where damages have been reported shows a success rate of 100%.

Efectos del cambio climático y usos del suelo sobre los recursos hídricos de la cuenca del río Tordera (Barcelona, España)

Se presentan los efectos del cambio global en la cuenca del río Tordera (España) para el periodo 2000-2050, escenarios climáticos A2 (medio-alto) definidos por el Panel Intergubernamental del Cambio Climático (IPCC, 200) y escenarios socioeconómicos (cambios previstos en la cuenca) denominados estable y tendencial. Los efectos sobre los recursos hídricos se han analizado de forma conjunta superficial-subterránea mediante una metodológica de tipo acoplado. Para establecer los impactos futuros sobre los recursos hídricos se ha seleccionado el Modelo de Circulación Global ECHAM5 (Max Planck Institute). Los resultados obtenidos indican una disminución de la precipitación del 11.3% y un aumento de la temperatura de 1ºC, respecto a los valores históricos de la zona. De acuerdo a la proyección futura (2050) sobre cambios en los recursos hídricos, la escorrentía superficial obtenida mediante simulación con el código HEC-HMS 3.4 experimenta una reducción del 31.8% respecto al valor histórico y la recarga natural, estimada mediante VISUAL-Balan, se reduce en un 11.7%. El balance en el acuífero deltaico simulado mediante MODFLOW 2009.1 Pro muestra igualmente una disminución de los parámetros del balance. Los cambios del uso del suelo previstos de acuerdo a la legislación vigente (escenarios socioeconómicos) no conducen a la generación de un impacto apreciable en los recursos hídricos; según los escenarios definidos la variación de precipitación y temperatura son los parámetros fundamentales del cambio previsto.

Twenty-year advanced DInSAR analysis of severe land subsidence: The Alto Guadalentín Basin (Spain) case study

A twenty-year period of severe land subsidence evolution in the Alto Guadalentín Basin (southeast Spain) is monitored using multi-sensor SAR images, processed by advanced differential interferometric synthetic aperture radar (DInSAR) techniques. The SAR images used in this study consist of four datasets acquired by ERS-1/2, ENVISAT, ALOS and COSMO-SkyMed satellites between 1992 and 2012. The integration of ground surface displacement maps retrieved for different time periods allows us to quantify up to 2.50 m of cumulated displacements that occurred between 1992 and 2012 in the Alto Guadalentín Basin. DInSAR results were locally compared with global positioning system (GPS) data available for two continuous stations located in the study area, demonstrating the high consistency of local vertical motion measurements between the two different surveying techniques. An average absolute error of 4.6 ± 4 mm for the ALOS data and of 4.8 ± 3.5 mm for the COSMO-SkyMed data confirmed the reliability of the analysis. The spatial analysis of DInSAR ground surface displacement reveals a direct correlation with the thickness of the compressible alluvial deposits. Detected ground subsidence in the past 20 years is most likely a consequence of a 100–200 m groundwater level drop that has persisted since the 1970s due to the overexploitation of the Alto Guadalentín aquifer system. The negative gradient of the pore pressure is responsible for the extremely slow consolidation of a very thick (> 100 m) layer of fine-grained silt and clay layers with low vertical hydraulic permeability (approximately 50 mm/h) wherein the maximum settlement has still not been reached.

Regional subsidence modelling in Murcia city (SE Spain) using 1-D vertical finite element analysis and 2-D interpolation of ground surface displacements

Subsidence is a hazard that may have natural or anthropogenic origin causing important economic losses. The area of Murcia city (SE Spain) has been affected by subsidence due to groundwater overexploitation since the year 1992. The main observed historical piezometric level declines occurred in the periods 1982–1984, 1992–1995 and 2004–2008 and showed a close correlation with the temporal evolution of ground displacements. Since 2008, the pressure recovery in the aquifer has led to an uplift of the ground surface that has been detected by the extensometers. In the present work an elastic hydro-mechanical finite element code has been used to compute the subsidence time series for 24 geotechnical boreholes, prescribing the measured groundwater table evolution. The achieved results have been compared with the displacements estimated through an advanced DInSAR technique and measured by the extensometers. These spatio-temporal comparisons have showed that, in spite of the limited geomechanical data available, the model has turned out to satisfactorily reproduce the subsidence phenomenon affecting Murcia City. The model will allow the prediction of future induced deformations and the consequences of any piezometric level variation in the study area.

A hydrological–economic model for sustainable groundwater use in sparse-data drylands: Application to the Amtoudi Oasis in southern Morocco, northern Sahara

A hydrological–economic model is introduced to describe the dynamics of groundwater-dependent economics (agriculture and tourism) for sustainable use in sparse-data drylands. The Amtoudi Oasis, a remote area in southern Morocco, in the northern Sahara attractive for tourism and with evidence of groundwater degradation, was chosen to show the model operation. Governing system variables were identified and put into action through System Dynamics (SD) modeling causal diagrams to program basic formulations into a model having two modules coupled by the nexus ‘pumping’: (1) the hydrological module represents the net groundwater balance (G) dynamics; and (2) the economic module reproduces the variation in the consumers of water, both the population and tourists. The model was operated under similar influx of tourists and different scenarios of water availability, such as the wet 2009–2010 and the average 2010–2011 hydrological years. The rise in international tourism is identified as the main driving force reducing emigration and introducing new social habits in the population, in particular concerning water consumption. Urban water allotment (PU) was doubled for less than a 100-inhabitant net increase in recent decades. The water allocation for agriculture (PI), the largest consumer of water, had remained constant for decades. Despite that the 2-year monitoring period is not long enough to draw long-term conclusions, groundwater imbalance was reflected by net aquifer recharge (R) less than PI + PU (G < 0) in the average year 2010–2011, with net lateral inflow from adjacent Cambrian formations being the largest recharge component. R is expected to be much less than PI + PU in recurrent dry spells. Some low-technology actions are tentatively proposed to mitigate groundwater degradation, such as: wastewater capture, treatment, and reuse for irrigation; storm-water harvesting for irrigation; and active maintenance of the irrigation system to improve its efficiency.

Vinasse application to sugar cane fields. Effect on the unsaturated zone and groundwater at Valle del Cauca (Colombia)

Extensive application of vinasse, a subproduct from sugar cane plantations for bioethanol production, is currently taking place as a source of nutrients that forms part of agricultural management in different agroclimatic regions. Liquid vinasse composition is characterised by high variability of organic compounds and major ions, acid pH (4.7), high TDS concentration (117,416–599,400 mg L− 1) and elevated EC (14,350–64,099 μS cm− 1). A large-scale sugar cane field application is taking place in Valle del Cauca (Colombia), where monitoring of soil, unsaturated zone and the aquifer underneath has been made since 2006 to evaluate possible impacts on three experimental plots. For this assessment, monitoring wells and piezometers were installed to determine groundwater flow and water samples were collected for chemical analysis. In the unsaturated zone, tensiometers were installed at different depths to determine flow patterns, while suction lysimeters were used for water sample chemical determinations. The findings show that in the sandy loam plot (Hacienda Real), the unsaturated zone is characterised by low water retention, showing a high transport capacity, while the other two plots of silty composition presented temporal saturation due to La Niña event (2010–2011). The strong La Niña effect on aquifer recharge which would dilute the infiltrated water during the monitoring period and, on the other hand dissolution of possible precipitated salts bringing them back into solution may occur. A slight increase in the concentration of major ions was observed in groundwater (~ 5% of TDS), which can be attributed to a combination of factors: vinasse dilution produced by water input and hydrochemical processes along with nutrient removal produced by sugar cane uptake. This fact may make the aquifer vulnerable to contamination.

Subsidence damage assessment of a Gothic church using differential interferometry and field data

The Santas Justa and Rufina Gothic church (fourteenth century) has suffered several physical, mechanical, chemical, and biochemical types of pathologies along its history: rock alveolization, efflorescence, biological activity, and capillary ascent of groundwater. However, during the last two decades, a new phenomenon has seriously affected the church: ground subsidence caused by aquifer overexploitation. Subsidence is a process that affects the whole Vega Baja of the Segura River basin and consists of gradual sinking in the ground surface caused by soil consolidation due to a pore pressure decrease. This phenomenon has been studied by differential synthetic aperture radar interferometry techniques, which illustrate settlements up to 100 mm for the 1993–2009 period for the whole Orihuela city. Although no differential synthetic aperture radar interferometry information is available for the church due to the loss of interferometric coherence, the spatial analysis of nearby deformation combined with fieldwork has advanced the current understanding on the mechanisms that affect the Santas Justa and Rufina church. These results show the potential interest and the limitations of using this remote sensing technique as a complementary tool for the forensic analysis of building structures.

Imaging Land Subsidence Induced by Groundwater Extraction in Beijing (China) Using Satellite Radar Interferometry

Beijing is one of the most water-stressed cities in the world. Due to over-exploitation of groundwater, the Beijing region has been suffering from land subsidence since 1935. In this study, the Small Baseline InSAR technique has been employed to process Envisat ASAR images acquired between 2003 and 2010 and TerraSAR-X stripmap images collected from 2010 to 2011 to investigate land subsidence in the Beijing region. The maximum subsidence is seen in the eastern part of Beijing with a rate greater than 100 mm/year. Comparisons between InSAR and GPS derived subsidence rates show an RMS difference of 2.94 mm/year with a mean of 2.41 ± 1.84 mm/year. In addition, a high correlation was observed between InSAR subsidence rate maps derived from two different datasets (i.e., Envisat and TerraSAR-X). These demonstrate once again that InSAR is a powerful tool for monitoring land subsidence. InSAR derived subsidence rate maps have allowed for a comprehensive spatio-temporal analysis to identify the main triggering factors of land subsidence. Some interesting relationships in terms of land subsidence were found with groundwater level, active faults, accumulated soft soil thickness and different aquifer types. Furthermore, a relationship with the distances to pumping wells was also recognized in this work.

Uso do zoneamento ambiental do Médio Rio Grande na sub-bacia hidrográfica do Ribeirão Bocaina - Passos/MG

Este trabalho é parte integrante do Projeto Grande Minas - União pelas Águas, que realizou o Zoneamento Ambiental das Sub-bacias Hidrográficas dos Afluentes Mineiros do Médio Rio Grande. Após a obtenção dos produtos finais do zoneamento, entra-se agora, na fase de “Implementação”, ou seja, de sua aplicação direta para contribuir no processo de gestão dos recursos hídricos das sub-bacias hidrográficas da área de estudo. Neste trabalho estuda-se, especificamente, a sub-bacia hidrográfica do Ribeirão Bocaina, considerada uma das 34 sub-bacias que envolvem a Bacia Hidrográfica dos Afluentes Mineiros do Médio Rio Grande e situa-se no município de Passos-MG, possuindo uma área de 457,9 km². O trabalho busca dar uma contribuição para a gestão dos recursos hídricos nesta sub-bacia, que vem sofrendo com a degradação do seu recurso hídrico. O objetivo é disponibilizar instrumentos cartográficos para subsidiar a gestão dos recursos hídricos e, com isto, contribuir com a sua preservação e uso sustentável. A sistemática metodológica envolve uma análise dos aspectos legais da área de estudo; a construção de um Banco de Dados Digital georreferenciado com informações sobre os meios físico, biótico e socioeconômico da sub-bacia e a produção de uma carta derivada, de cunho interpretativo e de fácil leitura, que possa ser utilizada diretamente pelo gestor público nas tomadas de decisões. A base digital produzida conta com uma série de mapas digitais, dentre eles: mapas climáticos, de solos, declividades, geomorfológico, geológico, de sistemas aquíferos, hidrográfico, uso e ocupação do solo. Com relação à carta derivada, destaca os terrenos que apresentam maior predisposição à alterações diretas nos recursos hídricos.

ROCK OUTCROPS IN THE CANADIAN SHIELD: AN INVESTIGATION OF CONTAMINANT TRANSPORT FROM SURFACE SOURCES IN FRACTURED ROCK AQUIFERS

Flow, recharge and transport dynamics in fractured rock aquifers with low lying rock outcrops is a largely unexplored area of study in hydrogeology. The purpose of this thesis is to examine these topics in an agricultural area in Eastern Ontario. The study consists of a regional scale groundwater quality study, an infiltration experiment that considers bacteria transport from the ground surface to a well, and a numerical modelling study that tests the parameters that affect surface infiltration of a tracer from a rock outcrop to a deeper horizontal fracture. In the water quality study, approximately 65% of the samples contained total coliform, 16% contained E. coli, and 1% contained nitrate-N at greater than 5 mg/L. Occurrence of E. coli increased when considering seasonality, where wells were drilled on rock outcrops, and for shallow well intervals. Nitrate-N did not occur above the Guidelines for Canadian Drinking Water Quality (Health Canada, 2012) of 10 mg/L. Rapid arrival times were observed in the infiltration study for both the microspheres (30 minutes) and a dye tracer (45 minutes) in a well approximately 6.0 m in horizontal and 2.8 m in vertical distance from the tracer source. Transport velocities were approximately 38.9 m/day for the dye tracer and 115.2 m/day for the colloidal tracer. Results of the model runs indicate that overburden can provide an effective protective layer to transport in fractures, that high groundwater velocities occur in larger fracture apertures and higher gradients dilute tracer concentrations, and that lower groundwater velocities occur with smaller fracture apertures and lower gradients result in elevated tracer concentrations. Lower rainfall rates, larger fracture apertures, early tracer time, larger gradients, and lower water levels maintained unsaturated conditions for longer time periods such that tracer transport was delayed until saturated conditions were attained. The overall heterogeneity of this aquifer environment creates a source water protection conundrum where the water quality is generally good, while transport can occur very quickly in proximity to rock outcrops and in areas with limited overburden.

Characterization of residual soil used for infiltration of reclaimed water

The effects of climate change and the growing demand for water for domestic, industrial, agricultural and recreational activities have been led the use of treated wastewater (reclaimed water) for such applications. The artificial recharge of aquifers with treated wastewater can be an alternative way for restoring underground water volumes that can be used for satisfying some activities, particularly in water shortage areas or where their quality is not suitable for use. After a two-year monitoring period in the Vila Fernando WWTP (Guarda, Portugal), the characteristics of the treated effluent suggest that it could be used for infiltration for aquifer recharge. A multi-criteria analysis based on GIS was developed for site location of infiltration sites. The procedure has involved the combination of six thematic maps and environmental, technical and economic criteria, over an area of 6687.1 ha. About 6.4 ha were selected for suitable sites for infiltration and one of these sites (Quinta de Gonçalo Martins, Guarda) was selected for collecting soil samples. The characterization of the soil indicates that is favorable to the infiltration of treated wastewater for artificial recharge of aquifers.

Metodologias utilizadas no estudo do escoamento em aquíferos cársicos e o caso prático da captação do Olho de Mira (Maciço Calcário Estremenho)

Tese de mestrado em Geologia Aplicada, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2016

Acompanhamento da construção de uma captação de água subterrânea

Tese de mestrado, Geologia Aplicada (Hidrogeologia) Universidade de Lisboa, Faculdade de Ciências, 2016

Isotope Study of Moisture Sources, Recharge Areas and Groundwater Flowpaths within the Eastern Batinah Coastal Plain, Sultanate of Oman

Oxygen and hydrogen isotope analyses of rainfall samples collected on the eastern Batinah coastal plain of northern Oman between 1995 and 1998 indicate two different principal water vapor sources for precipitation in the area: a northern, Mediterranean source and a southern, Indian Ocean source. As a result, two new local meteoric water lines were defined for the study area. Isotopic analyses of groundwater samples from over 200 springs and wells indicate that the main source of water to the Batinah coastal alluvial aquifer is high-altitude rainfall from the adjacent Jabal Akhdar Mountains, originating from a combination of northern and southern moisture sources. The groundwater recharged at high-altitude forms two plumes of water which is depleted in the heavy isotopes 18O and 2H and stretches from the mountains across the coastal plain to the sea, thereby retaining a chemical homogeneity horizontally and vertically down to a depth exceeding 300 m. In contrast, in areas adjacent to these two plumes the alluvial aquifer is geochemically stratified. Near the coast, saline intrusion results in abrupt changes in chloride concentrations and isotope values.

(Table 1) Gravity and pressure within the Sleipner area relative to ones at the reference station in 2002 and 2005

At Sleipner, CO2 is being separated from natural gas and injected into an underground saline aquifer for environmental purposes. Uncertainty in the aquifer temperature leads to uncertainty in the in situ density of CO2. In this study, gravity measurements were made over the injection site in 2002 and 2005 on top of 30 concrete benchmarks on the seafloor in order to constrain the in situ CO2 density. The gravity measurements have a repeatability of 4.3 µGal for 2003 and 3.5 µGal for 2005. The resulting time-lapse uncertainty is 5.3 µGal. Unexpected benchmark motions due to local sediment scouring contribute to the uncertainty. Forward gravity models are calculated based on both 3D seismic data and reservoir simulation models. The time-lapse gravity observations best fit a high temperature forward model based on the time-lapse 3D seismics, suggesting that the average in situ CO2 density is about to 530kg/m**3. Uncertainty in determining the average density is estimated to be ±65 kg/m**3 (95% confidence), however, this does not include uncertainties in the modeling. Additional seismic surveys and future gravity measurements will put better constraints on the CO2 density and continue to map out the CO2 flow.