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 The compaction measurements of Quaternary and Tertiary Gulf Coast aquifer system sediments in the Houston-Galveston region (TX) show spatially variable compression of 0.08 to 8.49 mm/yr because of geohistorical overburden pressure when groundwater levels in the aquifer system were stable after about the year 2000. An aquifer-system creep equation is developed for evaluating this variable compression, with a thickness-weighted average creep coefficient based on Taylor's (1942) secondary consolidation theory. The temporal variation of aquifer system creep can be neglected in a short-term observation period (such as a decade) after a long-term creep period (such as over 1,000 years) in geohistory. The creep coefficient of the Gulf Coast aquifer system is found to be in a range of 8.74 × 10−5 to 3.94 × 10−3 (dimensionless), with an average of 1.38 × 10−3. Moreover, for silty clay or clay-dominant aquitards in the Gulf Coast aquifer system the creep coefficient value varies in the range of 2.21 × 10−4 to 3.94 × 10−3, which is consistent with values found by Mesri (1973) for most soils, which vary in the range of creep coefficient, 1 × 10−4 to 5 × 10−3. Land subsidence due to secondary consolidation of the Gulf Coast aquifer system is estimated to be 0.04 to 4.33 m in the 20th century and is projected to be 0.01 to 0.64 m in the 21st century at the 13 borehole extensometer locations in the Houston-Galveston region. The significant creep should be considered in the relative sea level rise, in addition to tectonic subsidence and primary consolidation.more » « less
An official website of the United States government

