A total of 14 extensometers were installed in Houston-Galveston Region, Texas, USA, at 12 locations to record compaction. The earliest extensometer began to record compaction in 1973. Records from three of the extensometers installed at Baytown (Shallow and Deep) and Pasadena exhibit anomalous subsidence from 2009 to 2017. The maximum compaction occurred around 2014 with Baytown Shallow recording 164 mm, Baytown Deep 72 mm, and Pasadena 135 mm. The anomalous subsidence exhibits features not related to primary consolidation subsidence (PCS) and secondary consolidation subsidence (SCS) of the Gulf Coast Aquifer System (GCAS). Groundwater level records at the extensometer locations indicate that the anomalous subsidence is not related to groundwater exploitation and creep of the GCAS in this region. Analysis of compaction data for the three sites indicates that the subsidence is partially elastic. Salt dome growth/evolution resulting in activation/reactivation of subsurface and surface faults is proposed as the mechanism responsible for the anomalous subsidence.
more »
« less
Coupled Stress-Dependent Groundwater Flow-Deformation Model to Predict Land Subsidence in Basins with Highly Compressible Deposits
In this study, a stress-dependent groundwater model, MODFLOW-SD, has been developed and coupled with the nonlinear subsidence model, NDIS, to predict vertical deformation occurring in basins with highly compressible deposits. The MODFLOW-SD is a modified version of MODFLOW (the USGS Modular Three-Dimensional Groundwater Flow Model) with two new packages, NONK and NONS, to update hydraulic conductivity and skeletal specific storage due to change in effective stress. The NDIS package was developed based on Darcy–Gersevanov Law and bulk flux to model land subsidence. Results of sample simulations run for a conceptual model showed that hydraulic heads calculated by MODFLOW significantly overestimated for confining units and slightly underestimated for aquifer ones. Moreover, it showed that applied stress due to pumping changed initially homogeneous layers to be heterogeneous ones. Comparison of vertical deformations calculated by NDIS andMODFLOW-SUB showed that neglecting horizontal strain and stress-dependency of aquifer parameters can overestimate future subsidence. Furthermore, compared to the SUB (Subsidence and Aquifer-System Compaction) package, NDIS is more likely to provide a more accurate compaction model for a complex aquifer system with vertically variable compression (Cc), recompression (Cr), and hydraulic conductivity change (Ck) indices.
more »
« less
- Award ID(s):
- 1832065
- PAR ID:
- 10168032
- Date Published:
- Journal Name:
- Hydrology
- Volume:
- 6
- Issue:
- 78
- ISSN:
- 2306-5338
- Page Range / eLocation ID:
- 1-17
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
As much as 3.05 m of land subsidence was observed in 1979 in the Houston-Galveston region as a result primarily of inelastic compaction of aquitards in the Chicot and Evangeline aquifers between 1937 and 1979. The preconsolidation pressure heads for aquitards within these two aquifers were continuously updated in response to lowering groundwater levels, which in turn was caused by continuously increasing groundwater withdrawal rates from 0.57 to 4.28 million m3/day. This land subsidence occurred without any management of changes in groundwater levels. However, the management of recovering groundwater levels from 1979 to 2000 successfully decreased inelastic compaction from about 40 mm/yr in the early 1980s to zero around 2000 through decreasing groundwater withdrawal rates from 4.3 to 3.0 million m3/day. The inelastic consolidation that had existed for about 63 years roughly from 1937 to 2000 caused a land subsidence hazard in this region. Some rebounding of the land surface was achieved from groundwater level recovering management. It is found in this paper that subsidence of 0.08 to 8.49 mm/yr owing to a pseudo-constant secondary consolidation rate emerged or tended to emerge at 13 borehole extensometer station locations while the groundwater levels in the two aquifers were being managed. It is considered to remain stable in trend since 2000. The subsidence due to the secondary consolidation is beyond the control of any groundwater level change management schemes because it is caused by geo-historical overburden pressure on the two aquifers. The compaction measurements collected from the 13 extensometers since 1971 not only successfully corroborate the need for groundwater level change management in controlling land subsidence but also yield the first empirical findings of the occurrence of secondary consolidation subsidence in the Quaternary and Tertiary aquifer systems in the Houston-Galveston region.more » « less
-
We investigated how runoff-to-groundwater partitioning changes as a function of substrate age and degree of regolith development in the Island of Hawai’i, by modeling watershed-scale hydrodynamic properties for a series of volcanic catchments of different substrate age developed under different climates. In the younger catchments, rainfall infiltrates directly into the groundwater system and surface runoff is minimal, consisting of ephemeral streams flowing on the scale of hours to days. The older catchments show increasing surface runoff, with deeper incision and perennial discharge. We hypothesize that watershed-scale hydrodynamic properties change as a function of their weathering history—the convolution of time and climate: as surfaces age and become increasingly weathered, hydraulic conductivity is reduced, leading to increased runoff-to-recharge ratios. To test this relationship, we calculated both saturated hydraulic conductivity (k) and aquifer thickness (D) using recession flow analysis. We show that the average k in the younger catchments can be between 3 to 6 orders of magnitude larger than in older catchments, whereas modeled D increases with age. Ephemeral streams with zero baseflow at daily timescales cannot be evaluated using the same method. Instead, we calculated the recession constant for two contiguous catchments developed on young ash or lava deposits of different ages. Increasing bedrock age results in slower recession response in these ephemeral streams, which is consistent with decreasing hydraulic conductivity. Our results highlight the role of the weathering history in determining the evolution of watershed-scale hydrologic properties in volcanic catchments.more » « less
-
Abstract Many irrigated agricultural areas seek to prolong the lifetime of their groundwater resources by reducing pumping. However, it is unclear how lagged responses, such as reduced groundwater recharge caused by more efficient irrigation, may impact the long‐term effectiveness of conservation initiatives. Here, we use a variably saturated, simplified surrogate groundwater model to: (a) analyze aquifer responses to pumping reductions, (b) quantify time lags between reductions and groundwater level responses, and (c) identify the physical controls on lagged responses. We explore a range of plausible model parameters for an area of the High Plains aquifer (USA) where stakeholder‐driven conservation has slowed groundwater depletion. We identify two types of lagged responses that reduce the long‐term effectiveness of groundwater conservation, recharge‐dominated and lateral‐flow‐dominated, with vertical hydraulic conductivity (KZ) the major controlling variable. When highKZallows percolation to reach the aquifer, more efficient irrigation reduces groundwater recharge. By contrast, when lowKZimpedes vertical flow, short term changes in recharge are negligible, but pumping reductions alter the lateral flow between the groundwater conservation area and the surrounding regions (lateral‐flow‐dominated response). For the modeled area, we found that a pumping reduction of 30% resulted in median usable lifetime extensions of 20 or 25 years, depending on the dominant lagged response mechanism (recharge‐ vs. lateral‐flow‐dominated). These estimates are far shorter than estimates that do not account for lagged responses. Results indicate that conservation‐based pumping reductions can extend aquifer lifetimes, but lagged responses can create a sizable difference between the initially perceived and actual long‐term effectiveness.more » « less
-
Abstract Understanding subsurface structure and groundwater flow in deltaic aquifers is essential for evaluating the vulnerability of groundwater resources in delta systems. Deltaic aquifers contain coarse‐grained paleochannels that preserve a record of former surface river channels as well as fine‐grained floodplain deposits. The distribution of these deposits and how they are interconnected control groundwater flow and contaminant transport. In this work, we link depositional environments of deltaic aquifers to stratigraphic (static) and flow and transport (dynamic) connectivity metrics. Numerical models of deltaic stratigraphy were generated using a reduced‐complexity numerical model (DeltaRCM) with different input sand fractions (ISF) and rates of sea‐level rise (SLR). The groundwater flow and advective transport behavior of these deltas were simulated using MODFLOW and MODPATH. By comparing the static and dynamic metrics calculated from these numerical models, we show that groundwater behavior can be predicted by particular aspects of the subsurface architecture, and that horizontal and vertical connectivity display different characteristics. We also evaluate relationships between connectivity metrics and two environmental controls on delta evolution: ISF and SLR rate. The results show that geologic setting strongly influences both static and dynamic connectivity in different directions. These results provide insights into quantitatively differentiated subsurface hydraulic behavior between deltas formed under different external forcing (ISF and SLR rate) and they are a potential link in using information from delta surface networks and depositional history to predict vulnerability to aquifer contamination.more » « less
An official website of the United States government

