Geology topics
Benjamin B. Mirus
Publications and source records attributed to Benjamin B. Mirus.
The Galapagos archipelago: A natural laboratory to examine sharp hydroclimatic, geologic and anthropogenic gradients
Poor understanding of the water cycle in tropical ecosystems has the potential to exacerbate water shortages and water crises in the region. We suggest that the Galápagos Islands provide an excellent proxy to regions across the tropics as a result of sharp hydroclimatic, anthropogenic, and pedohydrologic gradients across the archipelago. Hydroclimatic and pedohydrologic gradients are found across different elevations on single islands, as well as across the archipelago, whereas anthropogenic gradients reflect land use and land cover change across islands as population and growth in tourism have affected individual islands differently. This article highlights specific opportunities to further examine our understanding of the interactions between water and critical zone processes in tropical ecosystems, making connections between the Galápagos archipelago and much of the understudied tropics. The Galápagos archipelago offers a natural laboratory through which we can examine current threats to freshwater security as well as the dynamics of coupled natural and human systems.
Testing the suitability of geologic frameworks for extrapolating hydraulic properties across regional scales
The suitability of geologic frameworks for extrapolating hydraulic conductivity ( K ) to length scales commensurate with hydraulic data is difficult to assess. A novel method is presented for evaluating assumed relations between K and geologic interpretations for regional-scale groundwater modeling. The approach relies on simultaneous interpretation of multiple aquifer tests using alternative geologic frameworks of variable complexity, where each framework is incorporated as prior information that assumes homogeneous K within each model unit. This approach is tested at Pahute Mesa within the Nevada National Security Site (USA), where observed drawdowns from eight aquifer tests in complex, highly faulted volcanic rocks provide the necessary hydraulic constraints. The investigated volume encompasses 40 mi 3 (167 km 3 ) where drawdowns traversed major fault structures and were detected more than 2 mi (3.2 km) from pumping wells. Complexity of the five frameworks assessed ranges from an undifferentiated mass of rock with a single unit to 14 distinct geologic units. Results show that only four geologic units can be justified as hydraulically unique for this location. The approach qualitatively evaluates the consistency of hydraulic property estimates within extents of investigation and effects of geologic frameworks on extrapolation. Distributions of transmissivity are similar within the investigated extents irrespective of the geologic framework. In contrast, the extrapolation of hydraulic properties beyond the volume investigated with interfering aquifer tests is strongly affected by the complexity of a given framework. Testing at Pahute Mesa illustrates how this method can be employed to determine the appropriate level of geologic complexity for large-scale groundwater modeling.
Identifying long term empirical relationships between storm characteristics and episodic groundwater recharge
Shallow aquifers are an important source of water resources and provide base flow to streams; yet actual rates of groundwater recharge are difficult to estimate. While climate change is predicted to increase the frequency and magnitude of extreme precipitation events, the resulting impact on groundwater recharge remains poorly understood. We quantify empirical relations between precipitation characteristics and episodic groundwater recharge for a wide variety of geographic and land use types across North Carolina. We extract storm duration, magnitude, average rate, and hourly weighted intensity from long-term precipitation records over periods of 12–35 years at 10 locations. Using time series of water table fluctuations from nearby monitoring wells, we estimate relative recharge to precipitation ratios (RPR) to identify statistical trends. Increased RPR correlates with increased storm duration, whereas RPR decreases with increasing magnitude, average rate, and intensity of precipitation. Agricultural and urban areas exhibit the greatest decrease in RPR due to increasing storm magnitude, average rate, and intensity, while naturally vegetated areas exhibit a larger increase in RPR with increased storm duration. Though RPR is generally higher during the winter than the summer, this seasonal effect is magnified in the Appalachian and Piedmont regions. These statistical trends provide valuable insights into the likely consequences of climate and land use change for water resources in subtropical climates. If, as predicted, growing seasons lengthen and the intensity of storms increases with a warming climate, decreased recharge in Appalachia, the Piedmont, and rapidly growing urban areas of the American Southeast could further limit groundwater availability.
An overview of current applications, challenges, and future trends in distributed process-based models in hydrology
Process-based hydrological models have a long history dating back to the 1960s. Criticized by some as over-parameterized, overly complex, and difficult to use, a more nuanced view is that these tools are necessary in many situations and, in a certain class of problems, they are the most appropriate type of hydrological model. This is especially the case in situations where knowledge of flow paths or distributed state variables and/or preservation of physical constraints is important. Examples of this include: spatiotemporal variability of soil moisture, groundwater flow and runoff generation, sediment and contaminant transport, or when feedbacks among various Earth’s system processes or understanding the impacts of climate non-stationarity are of primary concern. These are situations where process-based models excel and other models are unverifiable. This article presents this pragmatic view in the context of existing literature to justify the approach where applicable and necessary. We review how improvements in data availability, computational resources and algorithms have made detailed hydrological simulations a reality. Avenues for the future of process-based hydrological models are presented suggesting their use as virtual laboratories, for design purposes, and with a powerful treatment of uncertainty.
Evaluating the importance of characterizing soil structure and horizons in parameterizing a hydrologic process model
Incorporating the influence of soil structure and horizons into parameterizations of distributed surface water/groundwater models remains a challenge. Often, only a single soil unit is employed, and soil-hydraulic properties are assigned based on textural classification, without evaluating the potential impact of these simplifications. This study uses a distributed physics-based model to assess the influence of soil horizons and structure on effective parameterization. This paper tests the viability of two established and widely used hydrogeologic methods for simulating runoff and variably saturated flow through layered soils: (1) accounting for vertical heterogeneity by combining hydrostratigraphic units with contrasting hydraulic properties into homogeneous, anisotropic units and (2) use of established pedotransfer functions based on soil texture alone to estimate water retention and conductivity, without accounting for the influence of pedon structures and hysteresis. The viability of this latter method for capturing the seasonal transition from runoff-dominated to evapotranspiration-dominated regimes is also tested here. For cases tested here, event-based simulations using simplified vertical heterogeneity did not capture the state-dependent anisotropy and complex combinations of runoff generation mechanisms resulting from permeability contrasts in layered hillslopes with complex topography. Continuous simulations using pedotransfer functions that do not account for the influence of soil structure and hysteresis generally over-predicted runoff, leading to propagation of substantial water balance errors. Analysis suggests that identifying a dominant hydropedological unit provides the most acceptable simplification of subsurface layering and that modified pedotransfer functions with steeper soil-water retention curves might adequately capture the influence of soil structure and hysteresis on hydrologic response in headwater catchments.
Measurement of unsaturated hydraulic properties and evaluation of property-transfer models for deep sedimentary interbeds, Idaho National Laboratory, Idaho
Operations at the Idaho National Laboratory (INL) have the potential to contaminate the underlying Eastern Snake River Plain (ESRP) aquifer. Methods to quantitatively characterize unsaturated flow and recharge to the ESRP aquifer are needed to inform water-resources management decisions at INL. In particular, hydraulic properties are needed to parameterize distributed hydrologic models of unsaturated flow and transport at INL, but these properties are often difficult and costly to obtain for large areas. The unsaturated zone overlying the ESRP aquifer consists of alternating sequences of thick fractured volcanic rocks that can rapidly transmit water flow and thinner sedimentary interbeds that transmit water much more slowly. Consequently, the sedimentary interbeds are of considerable interest because they primarily restrict the vertical movement of water through the unsaturated zone. Previous efforts by the U.S. Geological Survey (USGS) have included extensive laboratory characterization of the sedimentary interbeds and regression analyses to develop property-transfer models, which relate readily available physical properties of the sedimentary interbeds (bulk density, median particle diameter, and uniformity coefficient) to water retention and unsaturated hydraulic conductivity curves. During 2013–14, the USGS, in cooperation with the U.S. Department of Energy, focused on further characterization of the sedimentary interbeds below the future site of the proposed Remote Handled Low-Level Waste (RHLLW) facility, which is intended for the long-term storage of low-level radioactive waste. Twelve core samples from the sedimentary interbeds from a borehole near the proposed facility were collected for laboratory analysis of hydraulic properties, which also allowed further testing of the property-transfer modeling approach. For each core sample, the steady-state centrifuge method was used to measure relations between matric potential, saturation, and conductivity. These laboratory measurements were compared to water-retention and unsaturated hydraulic conductivity parameters estimated using the established property-transfer models. For each core sample obtained, the agreement between measured and estimated hydraulic parameters was evaluated quantitatively using the Pearson correlation coefficient (r). The highest correlation is for saturated hydraulic conductivity ( K sat ) with an r value of 0.922. The saturated water content (q sat ) also exhibits a strong linear correlation with an r value of 0.892. The curve shape parameter (λ) has a value of 0.731, whereas the curve scaling parameter (y o ) has the lowest r value of 0.528. The r values demonstrate that model predictions correspond well to the laboratory measured properties for most parameters, which supports the value of extending this approach for quantifying unsaturated hydraulic properties at various sites throughout INL.
How runoff begins (and ends): characterizing hydrologic response at the catchment scale
Improved understanding of the complex dynamics associated with spatially and temporally variable runoff response is needed to better understand the hydrology component of interdisciplinary problems. The objective of this study was to quantitatively characterize the environmental controls on runoff generation for the range of different streamflow-generation mechanisms illustrated in the classic Dunne diagram. The comprehensive physics-based model of coupled surface-subsurface flow, InHM, is employed in a heuristic mode. InHM has been employed previously to successfully simulate the observed hydrologic response at four diverse, well-characterized catchments, which provides the foundation for this study. The C3 and CB catchments are located within steep, forested terrain; the TW and R5 catchments are located in gently sloping rangeland. The InHM boundary-value problems for these four catchments provide the corner-stones for alternative simulation scenarios designed to address the question of how runoff begins (and ends). Simulated rainfall-runoff events are used to systematically explore the impact of soil-hydraulic properties and rainfall characteristics. This approach facilitates quantitative analysis of both integrated and distributed hydrologic responses at high-spatial and temporal resolution over the wide range of environmental conditions represented by the four catchments. The results from 140 unique simulation scenarios illustrate how rainfall intensity/depth, subsurface permeability contrasts, characteristic curve shapes, and topography provide important controls on the hydrologic-response dynamics. The processes by which runoff begins (and ends) are shown, in large part, to be defined by the relative rates of rainfall, infiltration, lateral flow convergence, and storage dynamics within the variably saturated soil layers.
Balancing practicality and hydrologic realism: a parsimonious approach for simulating rapid groundwater recharge via unsaturated-zone preferential flow
The impact of preferential flow on recharge and contaminant transport poses a considerable challenge to water-resources management. Typical hydrologic models require extensive site characterization, but can underestimate fluxes when preferential flow is significant. A recently developed source-responsive model incorporates film-flow theory with conservation of mass to estimate unsaturated-zone preferential fluxes with readily available data. The term source-responsive describes the sensitivity of preferential flow in response to water availability at the source of input. We present the first rigorous tests of a parsimonious formulation for simulating water table fluctuations using two case studies, both in arid regions with thick unsaturated zones of fractured volcanic rock. Diffuse flow theory cannot adequately capture the observed water table responses at both sites; the source-responsive model is a viable alternative. We treat the active area fraction of preferential flow paths as a scaled function of water inputs at the land surface then calibrate the macropore density to fit observed water table rises. Unlike previous applications, we allow the characteristic film-flow velocity to vary, reflecting the lag time between source and deep water table responses. Analysis of model performance and parameter sensitivity for the two case studies underscores the importance of identifying thresholds for initiation of film flow in unsaturated rocks, and suggests that this parsimonious approach is potentially of great practical value.
Advanced methods for modeling water-levels and estimating drawdowns with SeriesSEE, an Excel add-in
Water-level modeling is used for multiple-well aquifer tests to reliably differentiate pumping responses from natural water-level changes in wells, or “environmental fluctuations.” Synthetic water levels are created during water-level modeling and represent the summation of multiple component fluctuations, including those caused by environmental forcing and pumping. Pumping signals are modeled by transforming step-wise pumping records into water-level changes by using superimposed Theis functions. Water-levels can be modeled robustly with this Theis-transform approach because environmental fluctuations and pumping signals are simulated simultaneously. Water-level modeling with Theis transforms has been implemented in the program SeriesSEE, which is a Microsoft® Excel add-in. Moving average, Theis, pneumatic-lag, and gamma functions transform time series of measured values into water-level model components in SeriesSEE. Earth tides and step transforms are additional computed water-level model components. Water-level models are calibrated by minimizing a sum-of-squares objective function where singular value decomposition and Tikhonov regularization stabilize results. Drawdown estimates from a water-level model are the summation of all Theis transforms minus residual differences between synthetic and measured water levels. The accuracy of drawdown estimates is limited primarily by noise in the data sets, not the Theis-transform approach. Drawdowns much smaller than environmental fluctuations have been detected across major fault structures, at distances of more than 1 mile from the pumping well, and with limited pre-pumping and recovery data at sites across the United States. In addition to water-level modeling, utilities exist in SeriesSEE for viewing, cleaning, manipulating, and analyzing time-series data.
Practical estimates of field-saturated hydraulic conductivity of bedrock outcrops using a modified bottomless bucket method
The bottomless bucket (BB) approach (Nimmo et al., 2009a) is a cost-effective method for rapidly characterizing field-saturated hydraulic conductivity Kfs of soils and alluvial deposits. This practical approach is of particular value for quantifying infiltration rates in remote areas with limited accessibility. A similar approach for bedrock outcrops is also of great value for improving quantitative understanding of infiltration and recharge in rugged terrain. We develop a simple modification to the BB method for application to bedrock outcrops, which uses a non-toxic, quick-drying silicone gel to seal the BB to the bedrock. These modifications to the field method require only minor changes to the analytical solution for calculating Kfs on soils. We investigate the reproducibility of the method with laboratory experiments on a previously studied calcarenite rock and conduct a sensitivity analysis to quantify uncertainty in our predictions. We apply the BB method on both bedrock and soil for sites on Pahute Mesa, which is located in a remote area of the Nevada National Security Site. The bedrock BB tests may require monitoring over several hours to days, depending on infiltration rates, which necessitates a cover to prevent evaporative losses. Our field and laboratory results compare well to Kfs values inferred from independent reports, which suggests the modified BB method can provide useful estimates and facilitate simple hypothesis testing. The ease with which the bedrock BB method can be deployed should facilitate more rapid in-situ data collection than is possible with alternative methods for quantitative characterization of infiltration into bedrock.
Assessing controls on perched saturated zones beneath the Idaho Nuclear Technology and Engineering Center, Idaho
Waste byproducts associated with operations at the Idaho Nuclear Technology and Engineering Center (INTEC) have the potential to contaminate the eastern Snake River Plain (ESRP) aquifer. Recharge to the ESRP aquifer is controlled largely by the alternating stratigraphy of fractured volcanic rocks and sedimentary interbeds within the overlying vadose zone and by the availability of water at the surface. Beneath the INTEC facilities, localized zones of saturation perched on the sedimentary interbeds are of particular concern because they may facilitate accelerated transport of contaminants. The sources and timing of natural and anthropogenic recharge to the perched zones are poorly understood. Simple approaches for quantitative characterization of this complex, variably saturated flow system are needed to assess potential scenarios for contaminant transport under alternative remediation strategies. During 2009-2011, the U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, employed data analysis and numerical simulations with a recently developed model of preferential flow to evaluate the sources and quantity of recharge to the perched zones. Piezometer, tensiometer, temperature, precipitation, and stream-discharge data were analyzed, with particular focus on the possibility of contributions to the perched zones from snowmelt and flow in the neighboring Big Lost River (BLR). Analysis of the timing and magnitude of subsurface dynamics indicate that streamflow provides local recharge to the shallow, intermediate, and deep perched saturated zones within 150 m of the BLR; at greater distances from the BLR the influence of streamflow on recharge is unclear. Perched water-level dynamics in most wells analyzed are consistent with findings from previous geochemical analyses, which suggest that a combination of annual snowmelt and anthropogenic sources (for example, leaky pipes and drainage ditches) contribute to recharge of shallow and intermediate perched zones throughout much of INTEC. The source-responsive fluxes model was parameterized to simulate recharge via preferential flow associated with intermittent episodes of streamflow in the BLR. The simulations correspond reasonably well to the observed hydrologic response within the shallow perched zone. Good model performance indicates that source-responsive flow through a limited number of connected fractures contributes substantially to the perched-zone dynamics. The agreement between simulated and observed perched-zone dynamics suggest that the source-responsive fluxes model can provide a valuable tool for quantifying rapid preferential flow processes that may result from different land management scenarios.