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.
Coastal deltaic aquifers are vulnerable to degradation from seawater intrusion, geogenic and anthropogenic contamination, and groundwater abstraction. The distribution and transport of contaminants are highly dependent on the subsurface sedimentary architecture, such as the presence of channelized features that preferentially conduct flow. Surface deposition changes in response to sea‐level rise (SLR) and sediment supply, but it remains unclear how these surface changes affect the distribution and transport of groundwater solutes in aquifers. Here, we explore the influence of SLR and sediment supply on aquifer heterogeneity and resulting effects on contaminant transport. We use realizations of subsurface heterogeneity generated by a process‐based numerical model, DeltaRCM, which simulates the evolution of a deltaic aquifer with different input sand fractions and rates of SLR. We simulate groundwater flow and solute transport through these deposits in three contamination scenarios: (a) vertical transport from widespread contamination at the land surface, (b) vertical transport from river water infiltration, and (c) lateral seawater intrusion. The simulations show that the vulnerability of deltaic aquifers to seawater intrusion correlates to sand fraction, while vertical transport of contaminants, such as widespread shallow contamination and river water infiltration, is influenced by channel stacking patterns. This analysis provides new insights into the connection between the depositional system properties and vulnerability to different modes of groundwater contamination. It also illustrates how vulnerability may vary locally within a delta due to depositional differences. Results suggest that groundwater management strategies may be improved by considering surface features, location within the delta, and the external forcings during aquifer deposition.
more » « less- NSF-PAR ID:
- 10376051
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 9
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract River deltas are densely populated regions of the world with vulnerable groundwater reserves. Contamination of these groundwater aquifers via saline water intrusion and pollutant transport is a growing threat due to both anthropogenic and climate changes. The arrangement and composition of subsurface sediment is known to have a significant impact on aquifer contamination; however, developing accurate depictions of the subsurface is challenging. In this work, we explore the relationship between surface and subsurface properties and identify the metrics most sensitive to different forcing conditions. To do so, we simulate river delta evolution with the rule‐based numerical model, DeltaRCM, and test the influence of input sand fraction and steady sea level rise (SLR) on delta evolution. From the model outputs, we measure a variety of surface and subsurface metrics chosen based on their applicability to imagery and modeling results. The Kullback‐Leibler (KL) divergence is then used to quantitatively gauge which metrics are most indicative of the imposed forcings. Both qualitative observations and the KL divergence analysis suggest that estimates of subsurface connectivity can be constrained using surface information. In particular, more variable shoreline roughness values and higher surface wetted fraction values correspond to increased subsurface connectivity. These findings complement traditional methods of estimating subsurface structure in river‐dominated delta systems and represent a step toward the identification of a direct link between surface observations and subsurface form.
-
more » « less
Abstract Climate change is raising sea levels across the globe. On river deltas, sea‐level rise (SLR) may result in land loss, saline intrusion into groundwater aquifers, and other problems that adversely impact coastal communities. There is significant uncertainty surrounding future SLR trajectories and magnitudes, even over decadal timescales. Given this uncertainty, numerical modeling is needed to explore how different SLR projections may impact river delta evolution. In this work, we apply the pyDeltaRCM numerical model to simulate 350 years of deltaic evolution under three different SLR trajectories: steady rise, an abrupt change in SLR rate, and a gradual acceleration of SLR. For each SLR trajectory, we test a set of six final SLR magnitudes between 5 and 40 mm/yr, in addition to control runs with no SLR. We find that both surface channel dynamics as well as aspects of the subsurface change in response to higher rates of SLR, even over centennial timescales. In particular, increased channel mobility due to SLR corresponds to higher sand connectivity in the subsurface. Both the trajectory and magnitude of SLR change influence the evolution of the delta surface, which in turn modifies the structure of the subsurface. We identify correlations between surface and subsurface properties, and find that inferences of subsurface structure from the current surface configuration should be limited to time spans over which the sea level forcing is approximately steady. As a result, this work improves our ability to predict future delta evolution and subsurface connectivity as sea levels continue to rise.
-
more » « less
Abstract A density‐dependent, variably saturated groundwater flow and solute transport model was used to investigate the influence of swash motions on subsurface flow and moisture dynamics in beach aquifers with heterogeneous distributions of hydraulic conductivity (
K ) and capillarity. The numerical simulations were performed within a Monte Carlo framework using field measurements conducted in the swash zone of a sandy beach. Our results show that heterogeneous capillarity causes spatially variable capillary rise above the groundwater table. In response to swash motions, heterogeneity creates capillary barriers that result in pockets of elevated moisture content beneath the swash zone. These moisture hotspots persist within the unsaturated zone even at ebb tide when the swash motions recede seaward. Heterogeneous capillarity also results in highly tortuous preferential flow paths and alters the flow rates from the sand surface to the water table. HeterogeneousK greatly enhances the seawater infiltration into the swash zone and modulates its spatial distribution along the beach surface. Due to heterogeneousK and capillarity, complex mixing patterns emerge. Both strain‐dominated and vorticity‐dominated flow regions develop and dissipate as tides and waves move across the beach surface. Complex mixing patterns of seawater percolating from the swash zone surface to the water table, with localized areas of high and low mixing intensities, are further demonstrated by analysis of dilution index. Our findings reveal the influence of geologic heterogeneity on swash zone moisture and flow dynamics, which may have important implications for sediment transport and chemical processing in beach aquifers. -
more » « less
Abstract Groundwater is the primary source of water in the Bengal Delta but contamination threatens this vital resource. In deltaic environments, heterogeneous sedimentary architecture controls groundwater flow; therefore, characterizing subsurface structure is a critical step in predicting groundwater contamination. Here, we show that surface information can improve the characterization of the nature and geometry of subsurface features, thus improving the predictions of groundwater flow. We selected three locations in the Bengal Delta with distinct surface river network characteristics—the lower delta with straighter tidal channels, the mid‐delta with meandering and braided channels, and the inactive delta with transitional sinuous channels. We used surface information, including channel widths, depths, and sinuosity, to create models of the subsurface with object‐based geostatistical simulations. We collected an extensive set of lithologic data and filled in gaps with newly drilled boreholes. Our results show that densely distributed lithologic data from active lower and mid‐delta are consistent with the object‐based models generated from surface information. In the inactive delta, metrics from object‐based models derived from surface geometries are not consistent with subsurface data. We further simulated groundwater flow and solute transport through the object‐based models and compared these with simulated flow through lithologic models based only on variograms. Substantial differences in flow and transport through the different geologic models show that geometric structure derived from surface information strongly influences groundwater flow and solute transport. Land surface features in active deltas are therefore a valuable source of information for improving the evaluation of groundwater vulnerability to contamination.
