An official website of the United States government
Abstract Both hydrological and geophysical data can be used to calibrate hillslope hydrologic models. However, these data often reflect hydrological dynamics occurring at disparate spatial scales. Their use as sole objectives in model calibrations may thus result in different optimum hydraulic parameters and hydrologic model behavior. This is especially true for mountain hillslopes where the subsurface is often heterogeneous and the representative elementary volume can be on the scale of several m3. This study explores differences in hydraulic parameters and hillslope‐scale storage and flux dynamics of models calibrated with different hydrological and geophysical data. Soil water content, groundwater level, and two time‐lapse electrical resistivity tomography (ERT) data sets (transfer resistance and inverted resistivity) from two mountain hillslopes in Wyoming, USA, are used to calibrate physics‐based surface–subsurface hydrologic models of the hillslopes. Calibrations are performed using each data set independently and all data together resulting in five calibrated parameter sets at each site. Model predicted hillslope runoff and internal hydrological dynamics vary significantly depending on the calibration data set. Results indicate that water content calibration data yield models that overestimate near‐surface water storage in mountain hillslopes. Groundwater level calibration data yield models that more reasonably represent hillslope‐scale storage and flux dynamics. Additionally, ERT calibration data yield models with reasonable hillslope runoff predictions but relatively poor predictions of internal hillslope dynamics. These observations highlight the importance of carefully selecting data for hydrologic model calibration in mountain environments. Poor selection of calibration data may yield models with limited predictive capability depending on modeling goals and model complexity.
more »
« less
Hydrogeophysical Inversion of Time‐Lapse ERT Data to Determine Hillslope Subsurface Hydraulic Properties
Abstract Time‐lapse electrical resistivity tomography (ERT) data are increasingly used to inform the hydrologic dynamics of mountainous environments at the hillslope scale. Despite their popularity and recent advancements in hydrogeophysical inversion methods, few studies have shown how time‐lapse ERT data can be used to determine hydraulic parameters of subsurface water flow models. This study uses synthetic and field‐collected, hillslope‐scale, time‐lapse ERT data to determine subsurface hydraulic properties of a two‐layer, physics‐based, 2‐D vertical flow model with predefined layer and boundary locations. Uncoupled and coupled hydrogeophysical inversion methods are combined with a fine‐earth fraction optimization scheme to reduce the number of parameters needing calibration and interpret the influence of the hydraulic parameters on the hydrologic model predictions. Inversions of synthetic ERT data recover the prescribed fine‐earth fraction bulk density to within 0.1 g cm−3. Field‐collected ERT data from a mountain hillslope result in hydrologic model dynamics that are consistent with previous studies and measured water content data but struggle to capture measured groundwater levels. The uncoupled hydrogeophysical inversion method is more sensitive to changes in hydraulic parameter values of the lower hydrologic model layer than the coupled hydrogeophysical inversion method. Time series of minimum objective function value simulations indicate that periodically collected ERT data may recover hydraulic parameters to a similar level of uncertainty as daily ERT data. Using simple hydrologic model domains within hydrogeophysical inversions shows promise for providing reasonable hydrologic predictions while maintaining relatively simple calibration schemes and should be explored further in future studies.
more »
« less
- Award ID(s):
- 1818550
- PAR ID:
- 10419210
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 4
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The belowground architecture of the critical zone (CZ) consists of soil and rock in various stages of weathering and wetness that acts as a medium for biological growth, mediates chemical reactions, and controls partitioning of hydrologic fluxes. Hydrogeophysical imaging provides unique insights into the geometries and properties of earth materials that are present in the CZ and beyond the reach of direct observation beside sparse wellbores. An improved understanding of CZ architecture can be achieved by leveraging the geophysical measurements of the subsurface. Creating categorical models of the CZ is valuable for driving hydrologic models and comparing belowground architectures between different sites to interpret weathering processes. The CZ architecture is revealed through a novel comparison of hillslopes by applying facies classification in the elastic-electric domain driven by surface-based hydrogeophysical measurements. Three pairs of hillslopes grouped according to common geologic substrates — granite, volcanic extrusive, and glacially altered — are classified by five different hydrofacies classes to reveal the relative wetness and weathering states. The hydrofacies classifications are robust to the choice of initial mean values used in the classification and noncontemporaneous timing of geophysical data acquisition. These results will lead to improved interdisciplinary models of CZ processes at various scales and to an increased ability to predict the hydrologic timing and partitioning. Beyond the hillslope scale, this enhanced capability to compare CZ architecture can also be exploited at the catchment scale with implications for improved understanding of the link between rock weathering, hydrochemical fluxes, and landscape morphology.more » « less
-
Abstract. Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.more » « less
-
Abstract The hydrologic dynamics and geomorphic evolution of watersheds are intimately coupled—runoff generation and water storage are controlled by topography and properties of the surface and subsurface, while also affecting the evolution of those properties over geologic time. However, the large disparity between their timescales has made it difficult to examine interdependent controls on emergent hydrogeomorphic properties, such as hillslope length, drainage density, and extent of surface saturation. In this study, we develop a new model coupling hydrology and landscape evolution to explore how runoff generation affects long‐term catchment evolution, and analyze numerical results using a nondimensional scaling framework. We focus on hydrologic processes dominating in humid climates where storm runoff primarily arises from shallow subsurface flow and from precipitation on saturated areas. The model solves hydraulic groundwater equations to predict the water‐table elevation given prescribed, constant groundwater recharge. Water in excess of the subsurface capacity for transport becomes overland flow, which generates shear stress on the surface and may detach and transport sediment. This affects the landscape form that in turn affects runoff generation. We show that (a) four dimensionless parameters describe the possible steady state landscapes that coevolve under steady recharge; (b) hillslope length increases with increasing transmissivity relative to the recharge rate; (c) three topographic metrics—steepness index, Laplacian curvature, and topographic index—together provide a basis for interpreting landscapes that have coevolved with runoff generated via shallow subsurface flow. Finally we discuss the possibilities and limitations for quantitative comparisons between the model results and real landscapes.more » « less
-
Abstract Hydrologic modeling has been a useful approach for analyzing water partitioning in catchment systems. It will play an essential role in studying the responses of watersheds under projected climate changes. Numerous studies have shown it is critical to include subsurface heterogeneity in the hydrologic modeling to correctly simulate various water fluxes and processes in the hydrologic system. In this study, we test the idea of incorporating geophysics‐obtained subsurface critical zone (CZ) structures in the hydrologic modeling of a mountainous headwater catchment. The CZ structure is extracted from a three‐dimensional seismic velocity model developed from a series of two‐dimensional velocity sections inverted from seismic travel time measurements. Comparing different subsurface models shows that geophysics‐informed hydrologic modeling better fits the field observations, including streamflow discharge and soil moisture measurements. The results also show that this new hydrologic modeling approach could quantify many key hydrologic fluxes in the catchment, including streamflow, deep infiltration, and subsurface water storage. Estimations of these fluxes from numerical simulations generally have low uncertainties and are consistent with estimations from other methods. In particular, it is straightforward to calculate many hydraulic fluxes or states that may not be measured directly in the field or separated from field observations. Examples include quickflow/subsurface lateral flow, soil/rock moisture, and deep infiltration. Thus, this study provides a useful approach for studying the hydraulic fluxes and processes in the deep subsurface (e.g., weathered bedrock), which needs to be better represented in many earth system models.more » « less
