Three-Dimensional Numerical Investigation of Pore Water Pressure and Deformation of Pumped Aquifer Systems

Excessive groundwater withdrawal has caused severe land subsidence worldwide. The pore water pressure and the deformation of pumped hydrostratigraphic units are complex. A fully coupled three-dimensional numerical simulation was carried out for different pumping plans in this paper. When groundwater is pumped from a confined aquifer, the great compaction occurs in the pumped aquifer and its upper and lower adjacent aquitard units. Land subsidence is smaller and the area affected by land subsidence is greater when groundwater is pumped from the deeper confined aquifer. The pore water pressure in the pumped confined aquifer changes immediately with pumpage. In the adjacent aquitard units, however, the pore water pressure increases in the early pumping time and decreases in the early recharging time. The decrease in the pore water pressure vertically spreads from the interface between aquitard and pumped aquifer to the other surface of the aquitard. The pumped aquifer compacts and rebounds immediately with pumping and non-pumping or recharging actions, while the compaction and rebounding of the aquitard units clearly lag behind. The compaction of the adjacent aquitard unit first occurs near the interface between aquitard and pumped aquifer units, and the compaction zone spreads outward as the pumping goes on. The aquitards may expand vertically within some zones. Due to the inelastic deformation of soil skeleton, different pumping plans result in different land subsidence. For the same net pumpage, maximal land subsidence and horizontal displacement are the smallest for constant discharge and the greatest for recharge-discharge cycle.     

Variation in the rate of land subsidence induced by groundwater extraction and its effect on the response pattern of soil microbial communities

Excessive extraction of groundwater leads to (irreversible) changes in the physical soil properties, causing land subsidence associated with soil compaction to occur. Using a combined image processing and field approach, we examined: (1) how variation in the land subsidence rate induces different soil compaction degrees; and (2) the response patterns of microbial communities to such variations. By using Sentinel Synthetic Aperture Radar image processing, we selected three locations that exhibited different land subsidence rates, including high (HSR), moderate (MSR), and low (LSR). Then, soil sampling was undertaken within these representative locations. Indicators of soil compaction, including total porosity, air-filled porosity, water-filled porosity, and bulk density, were measured. The soil microbial community was determined using qPCR and sequencing. The highest and lowest values for bulk density were observed in the HSR–MSR and LSR zones, respectively. The greatest values of total porosity and macropore volume were displayed in the LSR zone compared to other zones. Bacterial abundance in the LSR zone was significantly greater than that in the HSR and MSR zones. The relative abundances of bacterial taxa indirectly demonstrated that the anaerobic phyla were significantly increased (by 10–13%), and the aerobic phyla decreased (by 30–40%) in the HSR zone compared to the LSR zone. This result demonstrates that the aerobes declined as larger volumes of the soil became more anaerobic. Indeed, the increased abundance of anaerobes was not able to compensate for the larger decrease in the abundance of aerobes. Our work showed that at the increased rates of land subsidence, the abundance distribution of the microbial community critically declined. These findings highlight the critical impacts of increasing the land subsidence rate on the emergence of high soil compaction degrees, which can significantly affect the resilience thresholds of the microbial communities in dryland soils.     

Land Subsidence Near Oil and Gas Fields,Houston, Texasa

Subsidence profiles across 29 oil and gas fields in the 12,200-km² Houston, Texas, regional subsidence area, which is caused by decline of ground-water level, suggest that the contribution of petroleum withdrawal to local land subsidence is small. Despite large volumes of petroleum production, subsidence at most fields was not increased by oil and gas withdrawal. Local increases of subsidence were detected at only six fields—Alco-Mag, Chocolate Bayou, Goose Creek, Hastings, Mykawa, and South Houston. With the exception of the 1-m subsidence from 1917 to 1925 at Goose Creek, differential subsidence across oil and gas fields was smaller by a factor of two or more than subsidence caused by aquifer compaction. At four fields—Barbers Hill, Cedar Bayou, Humble, and Pierce Junction—subsidence was substantially less than in the surrounding area. Except for Cedar Bayou, these fields are associated with shallow salt domes that partly occupy the aquifer system; for the three fields, subsidence during the periods of record came to less than half the subsidence in the surrounding area. In addition to land subsidence, faults with an aggregate length of more than 240 km (150 mi) have offset the land surface in historical time. Natural geologic deformation, ground-water pumping, and petroleum withdrawal have all been considered as potential causes of the historical offset across these faults. The minor amount of localized land subsidence associated with oil and gas fields, however, suggests that petroleum withdrawal is not a major cause of the historical faulting, at least by a differential compaction mechanism.     

Storage coefficient revisited: is purely vertical strain a good assumption?

The storage coefficient that is used ubiquitously today was first defined by the analytical work of Theis and Jacob over a half-century ago. Inherent within this definition is the restriction of purely vertical compression of the aquifer during a reduction in pressure. The assumption is revisited and quantitatively evaluated by comparing numerical results using both one- and three-dimensional strain models in the presence of three-dimensional flow. Results indicate that (1) calculated hydraulic head values are nearly identical for both models; (2) the release of water from storage in terms of volume strain is nearly identical for both models and that the location of maximum production moves outward from the well as a function of time; (3) the vertical strain components are markedly different with at least 50% of the total volume of water pumped originating from horizontal strain (and increasing to as much as 70%); and (4) for the one-dimensional strain model to yield the necessary quantity of water to the pumped well, the resulting vertical compaction (land subsidence) is as much as four times greater and vertical strain is as much as 60% greater than the three-dimensional strain model. Results indicate that small changes in porosity resulting from horizontal strain can yield extremely large quantities of water to the pumping well. This study suggests that the assumption of purely vertical strain used in the definition of the storage coefficient is not valid.     

Simulating pumping-induced regional land subsidence with the use of InSAR and field data in the Toluca Valley, Mexico

A multidisciplinary approach is presented here for quantifying land subsidence in a heavily pumped aquifer system with complex stratigraphy. The methodology consists in incorporating Terzaghi’s 1D instantaneous compaction principle into a 3D groundwater flow model that is then applied and calibrated to reproduce observed hydraulic heads and compaction for the Toluca Valley, Mexico. Differential Interferometric Synthetic Aperture Radar (D-InSAR), a generated 3D-geological model, extensometers, monitoring wells, and available literature are used to constrain the model. The D-InSAR measured subsidence, extensometers, and numerical simulations of subsidence agree relatively well. Simulations show that since regional subsidence began in the mid 1960s there has been up to 2 m of subsidence in the industrial corridor, where heavy pumping and thick clay layers are found. This study shows that an approach using various sources of data is useful in estimating and constraining the vertical component of the inelastic skeletal specific storage.     

The value of subsidence data in ground water model calibration

The accurate estimation of aquifer parameters such as transmissivity and specific storage is often an important objective during a ground water modeling investigation or aquifer resource evaluation. Parameter estimation is often accomplished with changes in hydraulic head data as the key and most abundant type of observation. The availability and accessibility of global positioning system and interferometric synthetic aperture radar data in heavily pumped alluvial basins can provide important subsidence observations that can greatly aid parameter estimation. The aim of this investigation is to evaluate the value of spatial and temporal subsidence data for automatically estimating parameters with and without observation error using UCODE-2005 and MODFLOW-2000. A synthetic conceptual model (24 separate cases) containing seven transmissivity zones and three zones each for elastic and inelastic skeletal specific storage was used to simulate subsidence and drawdown in an aquifer with variably thick interbeds with delayed drainage. Five pumping wells of variable rates were used to stress the system for up to 15 years. Calibration results indicate that (1) the inverse of the square of the observation values is a reasonable way to weight the observations, (2) spatially abundant subsidence data typically produce superior parameter estimates under constant pumping even with observation error, (3) only a small number of subsidence observations are required to achieve accurate parameter estimates, and (4) for seasonal pumping, accurate parameter estimates for elastic skeletal specific storage values are largely dependent on the quantity of temporal observational data and less on the quantity of available spatial data.     

利用GPS与InSAR研究西安现今地面沉降与地裂缝时空演化特征

西安孕育严重的地面沉降及地裂缝灾害,严重制约着城市的现代化发展,本文采用GPS精密定位和InSAR遥感差分技术对西安地面沉降和地裂缝进行变形监测与分析,获取了西安地面沉降与地裂缝整体变形现状的珍贵信息,通过对这些变形信息的研究分析,揭示了西安现今地面沉降与地裂缝时空演化特征和机理：随着停止或限采地下水,西安地面沉降量级由20世纪90年代中期的最大年沉降速率20～30 cm/a减少到不足10 cm/a,且超过60%的沉降区域的年沉降速率已由90年代中期的5~8 cm/a减少到不足2 cm/a;原有的沉降中心大部分已不存在或大大减小;地裂缝在时空活动与分布上与地面沉降存在明显的关联性;地面沉降和地裂缝随着西安高新区的建设向南、西南、东南逐步扩展.     

Microbiological response to well pumping

To better constrain sampling strategies for observing biologically sensitive parameters in ground water, we vigorously pumped for 120 h a lightly pumped well completed in a confined glacial aquifer while observing how various physical and chemical parameters evolve in the water produced. The parameters commonly monitored when sampling a well stabilized within about an hour, after 5 wellbore volumes were produced; these parameters include temperature, pH, dissolved oxygen, oxidation-reduction potential (Eh), and electrical conductivity. The concentrations of ferrous iron, sulfide, and sulfate and various biological or biologically sensitive parameters, including the concentrations of dissolved hydrogen and methane, direct cell counts, and the microbial community profile, in contrast, required more than 8 h or 36 well volumes to stabilize. We interpret this result to mean that the zone of influence of the wellbore on biologic processes in the aquifer extends beyond the commonly recognized zone where physical properties are affected. A second period of adjustment of these biologically sensitive parameters began after about 50 h of pumping, following displacement of 230 wellbore volumes, and continued to the end of the experiment. During this period, the cell density and the composition of the microbial community suspended in the water samples changed. This finding indicates that the microbial community in and near the wellbore changed in response to pumping and the changes affected aspects of the composition of water produced from the well. The study demonstrates the importance of allowing adequate pumping time when sampling ground water for the analysis of biologically sensitive parameters.     

南亚热带典型调水型水库—广东大镜山水库的富营养化特征分析

调水型水库是一种以抽水入库为主要来水水源的水体,是沿海地区重要的供水水源地.为了解这类水库的富营养化特点,于2005年全年每月2次对地处我国南亚热带地区(广东珠海市)的大镜山水库进行采样和监测.监测和测定指标主要包括氮、磷营养盐浓度、水温、透明度及叶绿素a浓度等,结合水库水文数据对水库富营养化特征和主要的影响因素进行分析.结果表明,2005年,大镜山水库的富营养化状态TSI_M指数在45-53之间,水库处于中富营养状态,多数时间处于富营养状态.水体富营养化主要参数表现出明显的季节变化,即叶绿素a浓度和富营养化状态指数在早春和晚秋出现两个峰值,明显地与温带富营养化水体在夏季出现单个峰值的特征不同.调水入库增加了水库营养盐负荷的同时,也在很大程度上影响了水库水动力学过程,与夏季的集中强降雨一起成为影响该水库富营养化的关键因素,这些因素改变了浮游植物群落对营养盐的直接响应,导致叶绿素a浓度与总磷、总氮浓度之间呈弱相关关系,降水和调水量在时间上相对配置重要性决定了叶绿素a浓度与营养盐浓度的关系.     

Sinkholes Due to Groundwater Withdrawal in Tazerbo Wellfield,SE Libya

The desert of eastern Libya forms one of the most arid regions of the Sahara. The Great Man‐Made River Project (GMRP) was established. It transports millions of cubic meters of water a day from desert wellfields to the coastal cities, where over 80% of the population lives. The Tazerbo Wellfield is one of the wellfields designed within the GMRP, delivering water to the eastern coast of Libya through an underground pipe network. Tazerbo Wellfield consists of 108 production wells; each well was designed to pump 100 L/s. The planned total groundwater withdrawal from all wells is 1 million m³/d. The deep sandstone aquifer (Nubian sandstone) is covered by a thick mudstone‐siltstone aquitard and is being heavily pumped. The aquifer and fine‐grained sediments of the aquitard may be compacted resulting in land subsidence as a result of high exploitation. Local sinkholes have developed in the area of Tazerbo since the start of the pumping from the wellfield in 2004. These sinkholes have been caused mainly by lowering of the piezometric heads due to the withdrawal of groundwater. In this study, a hydrogeological investigation is presented about the effect of large groundwater pumping from the Nubian sandstone aquifer in Tazerbo Wellfield, SE Libya, based on physical parameters for 108 production wells and 23 observation wells.     

Hydrologic trade-offs in conjunctive use management

An aquifer, in a stream/aquifer system, acts as a storage reservoir for groundwater. Groundwater pumping creates stream depletion that recharges the aquifer. As wells in the aquifer are moved away from the stream, the aquifer acts to filter out annual fluctuations in pumping; with distance the stream depletion tends to become equal to the total pumping averaged as an annual rate, with only a small fluctuation. This is true for both a single well and an ensemble of wells. A typical growing season in much of the western United States is 3 to 4 months. An ensemble of irrigation wells spread more or less uniformly across an aquifer several miles wide, pumping during the growing season, will deplete the stream by approximately one-third of the total amount of water pumped during the growing season. The remaining two-thirds of stream depletion occurs outside the growing season. Furthermore, it takes more than a decade of pumping for an ensemble of wells to reach a steady-state condition in which the impact on the stream is the same in succeeding years. After a decade or more of pumping, the depletion is nearly constant through the year, with only a small seasonal fluctuation: ±10%. Conversely, stream depletion following shutting down the pumping from an ensemble of wells takes more than a decade to fully recover from the prior pumping. Effectively managing a conjunctive groundwater and surface water system requires integrating the entire system into a single management institution with a long-term outlook.     

The Influence of Clay Zones on Land Subsidence from Groundwater Pumping

The objective of this article is to analyze the influence of clay zones on subsidence from groundwater pumping. Finite element analyses were conducted on a sand‐only aquifer and a sand aquifer with two clay zones located at different distances from the well face. A model that accounts for recoverable and nonrecoverable strains was used to simulate the sand and clay. This model couples the groundwater flow with the stress‐deformation response of the aquifer materials. Each aquifer was pumped from a single well for a period of 6 months, and then the groundwater level was lowered gradually to an elevation below the elevation of the clay zones and kept there for 10 years. The groundwater level was then raised gradually back to the original elevation over a period of 10 years. The results of the analyses show that the ground surface subsidence profile is strongly influenced by the presence of the clays zones. The ground surface sags where these clay zones are present resulting in a wavy ground surface profile. Subsidence continued when pumping is stopped, albeit at a much slower rate than during pumping, and when the groundwater level is below the elevation of the clay zones. Clay zones further away from the well face lag the subsidence of clay zones nearer the well face because of lower changes in hydrostatic head. Sags in ground surface subsidence profile from groundwater pumping are indicators of the presence of low hydraulic conductive geological materials.     

Evaluation of irrigation water use efficiency using remote sensing in the middle reach of the Heihe river,in the semi‐arid Northwestern China

Irrigation of agricultural oases is the main water consumer in semi‐arid and arid regions of Northwestern China. The accurate estimation of evapotranspiration (ET) on the oases is extremely important for evaluating water use efficiency so as to reasonably allocate water resources, particularly in semi‐arid and arid areas. In this study, we integrated the soil moisture information into surface energy balance system (SEBS) for improving irrigated crop water consumption estimation. The new approach fed with the moderate resolution imaging spectro‐radiometer images mapped spatiotemporal ET on the oasis in the middle reach of the Heihe river. The daily ET outputs of the new approach were compared with those of the original SEBS using the eddy correlation observations, and the results demonstrate that the modified SEBS remedied the shortcoming of general overestimating ET without regard to soil water stress. Meanwhile, the crop planting structure and leaf area index spatiotemporal distribution in the studied region were derived from the high‐resolution Chinese satellite HJ‐1/CCD images for helping analyse the pattern of the monthly ET (ET_monthly). The results show that the spatiotemporal variation of ET_monthly is closely related to artificial irrigation and crop growth. Further evaluation of current irrigation water use efficiency was conducted on both irrigation district scale and the whole middle reach of the Heihe river. The results reveal that the average fraction of consumed water on irrigation district scale is 57% in 2012. The current irrigation water system is irrational because only 52% of the total irrigated amount was used to fulfil plant ET requirement and the rest of the irrigation water recharged into groundwater in the oasis in 2012. However, in view of the whole middle reach of the Heihe river, the irrigation water use efficiency could reach to 66% in 2012. But pumping groundwater for reused irrigation wastes mostly energy instead of water. An improved irrigation water allocation system according to actual ET requirement is needed to increase irrigation efficiency per cubic meter water resource in an effort to save both water and energy. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of the Land Use and Check Dams on Flood in Upper Catchment of Fuping Hydrological Station by Hydrological Modeling

Flood peaks and volumes had been detected a downward trend in Fuping hydrological station. To quantify the effects of check dams on flood peaks and volumes, a hydrological model integrating land use was established. The model performed well in flood processes simulation, and the Nash efficiency of the model was 0.72. Then the model was used to identify the comprehensive effects of land use and land cover change on flood processes by comparing the simulation results of the selected flood events under 1980 and 2000 land use and land cover conditions. 24.5, 37.7 and 51.3% decrease in flood peaks for flood events of greater than 10 years, 5~10 years and less than 5 years return periods, respectively, and 16.3, 27.9 and 28.5% decrease in flood volumes for the three groups flood events of different return periods. Contributions of land use change and check dams to decrease in flood peak and volume were simulated, and it was found that peak discharge and volume for each flood event responded differently to the two factors. The results in this study can provide valuable information on design flood calculation in the basin under land use and land cover change.     

Does Pumping Volume Affect the Concentration of Micropollutants in Groundwater Samples?

Information on concentrations of micropollutants (such as pharmaceuticals, pesticides, and industrial chemicals) in most highly dynamic riverbank filtration (RBF) systems is lacking, in contrast to data on standard chemical parameters. Sampling protocols have thus far been based on the stabilization of standard chemical parameters in relatively pristine environments. To determine whether groundwater samples for micropollutant analysis can be taken at a similar pumping volume as samples for testing standard chemical parameters in both environments, three groundwater monitoring wells in an RBF system were sampled at two points in time (after pumping of 3 well volumes and after pumping of 15 well volumes). Micropollutant concentrations were not significantly different between the two sampling points; therefore, appropriate samples can be drawn after pumping 3 well volumes. For a specific microbiological parameter (leucin incorporation), a statistically significant difference was found.     

Stream-aquifer interactions: evaluation of depletion volume and residual effects from ground water pumping

Numerical modeling techniques were used to simulate stream-aquifer interactions from seasonal ground water pumping. We used stream-aquifer models in which a shallow stream penetrates the top of an aquifer that discharges ground water to the stream as base flow. Because of the pumping, the volume of base flow discharged to the stream was reduced, and as the pumping continued, infiltration from the stream to the aquifer was induced. Both base-flow reduction and stream infiltration contributed to total stream depletion. We analyzed the depletion rates and volumes of the reduced base flow and induced stream infiltration during pumping and postpumping periods. Our results suggested that for a shallow penetrating stream with a low streambed conductance, base-flow reduction accounts for a significant percentage of the total stream depletion. Its residual effects in postpumping can last very long and may continue into the next pumping season for areas where recharge is nominal. In contrast, the contribution of the induced stream infiltration to the total stream depletion is much smaller, and its effects often become negligible shortly after pumping was stopped. For areas where surface recharge replenishes the aquifer, the residual effects of base-flow reduction and thus its depletion volume will be significantly reduced. A stream of large conductance has a high hydraulic connection to the aquifer, but the relationship between stream conductance and stream depletion is not linear.     

Identification and quantification of streamflow trends on the Canadian Prairies

Abstract

Trend analysis was performed on streamflow data for a collection of stations on the Canadian Prairies, in terms of spring and summer runoff volumes, peak flow rates and peak flow occurrences, as well as an annual volume measure, for analysis periods of 1966–2005, 1971–2005, and 1976–2005. The Mann-Kendall statistical test for trend and bootstrap resampling were used to identify the trends and to determine the field significance of the trends. Partial correlation analysis was used to identify relationships between hydrological variables that exhibit a significant trend and meteorological variables that exhibit a significant trend. Noteworthy results include decreasing trends in the spring snowmelt runoff event volume and peak flow, decreasing trends (earlier occurrence) in the spring snowmelt runoff event peak date and decreasing trends in the seasonal (1 March–31 October) runoff volume. These trends can be attributed to a combination of reductions in snowfall and increases in temperatures during the winter months.     

Regional impacts of climate change on irrigation water demands

This paper presents an approach to model the expected impacts of climate change on irrigation water demand in a reservoir command area. A statistical downscaling model and an evapotranspiration model are used with a general circulation model (GCM) output to predict the anticipated change in the monthly irrigation water requirement of a crop. Specifically, we quantify the likely changes in irrigation water demands at a location in the command area, as a response to the projected changes in precipitation and evapotranspiration at that location. Statistical downscaling with a canonical correlation analysis is carried out to develop the future scenarios of meteorological variables (rainfall, relative humidity (RH), wind speed (U₂), radiation, maximum (Tmax) and minimum (Tmin) temperatures) starting with simulations provided by a GCM for a specified emission scenario. The medium resolution Model for Interdisciplinary Research on Climate GCM is used with the A1B scenario, to assess the likely changes in irrigation demands for paddy, sugarcane, permanent garden and semidry crops over the command area of Bhadra reservoir, India. Results from the downscaling model suggest that the monthly rainfall is likely to increase in the reservoir command area. RH, Tmax and Tmin are also projected to increase with small changes in U₂. Consequently, the reference evapotranspiration, modeled by the Penman–Monteith equation, is predicted to increase. The irrigation requirements are assessed on monthly scale at nine selected locations encompassing the Bhadra reservoir command area. The irrigation requirements are projected to increase, in most cases, suggesting that the effect of projected increase in rainfall on the irrigation demands is offset by the effect due to projected increase/change in other meteorological variables (viz., Tmax and Tmin, solar radiation, RH and U₂). The irrigation demand assessment study carried out at a river basin will be useful for future irrigation management systems. Copyright © 2012 John Wiley & Sons, Ltd.     

Integrated water resource management in a major canal command in eastern India

The present rice‐dominated cropping system in the Hirakud canal command (eastern India) is under severe threat due to imbalance between irrigation water supply and demand. The canal water supply, which is the only source of irrigation, only meets 54% of the demand at 90% probability of exceedance (PE). In order to mitigate the irrigation water deficit from canal water, groundwater is considered as a supplemental source. Quasi‐three‐dimensional groundwater flow simulation modelling was, therefore, carried out by using Visual MODFLOW to detect the change in hydraulic head due to transient pumping stresses. The simulation model was calibrated and validated satisfactorily. Sensitivity analysis of the model parameters shows that groundwater recharge is most sensitive followed by aquifer hydraulic conductivity at almost all the sites of the command area, whereas the model is comparatively less sensitive to specific storage and specific yield. Enhanced pumping scenarios showed that groundwater extraction can be increased up to 50 times of the existing pumping without causing any adverse effect to the aquifer but the aquifer does not permit to exploit water in order to fulfill the irrigation water demand even at 10% PE. Hence, it is imperative to develop an optimal land and water resources management plan of the command area. Copyright © 2011 John Wiley & Sons, Ltd.