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
Land Subsidence Due to Creep of the Gulf Coast Aquifer System in the Houston-Galveston Region
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
- PAR ID:
- 10422560
- Date Published:
- Journal Name:
- Environmental & Engineering Geoscience
- Volume:
- 28
- Issue:
- 3
- ISSN:
- 1078-7275
- Page Range / eLocation ID:
- 237 to 254
- 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
-
Abstract Relative sea level rise at tide gauge Galveston Pier 21, Texas, is the combination of absolute sea level rise and land subsidence. We estimate subsidence rates of 3.53 mm/a during 1909–1937, 6.08 mm/a during 1937–1983, and 3.51 mm/a since 1983. Subsidence attributed to aquifer-system compaction accompanying groundwater extraction contributed as much as 85% of the 0.7 m relative sea level rise since 1909, and an additional 1.9 m is projected by 2100, with contributions from land subsidence declining from 30 to 10% over the projection interval. We estimate a uniform absolute sea level rise rate of 1.10 mm ± 0.19/a in the Gulf of Mexico during 1909–1992 and its acceleration of 0.270 mm/a2at Galveston Pier 21 since 1992. This acceleration is 87% of the value for the highest scenario of global mean sea level rise. Results indicate that evaluating this extreme scenario would be valid for resource-management and flood-hazard-mitigation strategies for coastal communities in the Gulf of Mexico, especially those affected by subsidence.more » « less
-
Abstract We present the first investigation of subsidence due to sediment compaction and consolidation in two laboratory‐scale river delta experiments. Spatial and temporal trends in subsidence rates in the experimental setting may elucidate behavior which cannot be directly observed at sufficiently long timescales, except for in reduced scale models such as the ones studied. We compare subsidence between a control experiment using steady boundary conditions, and an otherwise identical experiment which has been treated with a proxy for highly compressible marsh deposits. Both experiments have non‐negligible compactional subsidence rates across the delta‐top, comparable in magnitude to our boundary condition relative sea level rise rate of 250 μm/hr. Subsidence in the control experiment (on average 54 μm/hr) is concentrated in the lowest elevation (<10 mm above sea level) areas near the coast and is likely related to creep induced by a rising water table near the shoreface. The treatment experiment exhibits larger (on average 126 μm/hr) and more spatially variable subsidence rates controlled mostly by compaction of recent marsh deposits within one channel depth (∼10 mm) of the sediment surface. These rates compare favorably with field and modeling based subsidence measurements both in relative magnitude and location. We find that subsidence “hot spots” may be relatively ephemeral on longer timescales, but average subsidence across the entire delta can be variable even at our shortest measurement window. This suggests that subsidence rates over a short time frame may exceed thresholds for marsh platform drowning, even if the long term trend does not.more » « less
-
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
An official website of the United States government

