Brian A. Ebel

Near-surface hydrologic response for a steep, unchanneled catchment near Coos Bay, Oregon: 2. Physics-based simulations

The comprehensive physics-based hydrologic-response model InHM was used to simulate 3D variably-saturated flow and solute transport for three controlled sprinkling experiments at the Coos Bay 1 (CB1) experimental catchment in the Oregon Coast Range. The InHM-simulated hydrologic-response was evaluated against observed discharge, pressure head, total head, soil-water content, and deuterium concentration records. Runoff generation, tensiometric/piezometric response in the soil, pore-water pressure generation, and solute (tracer) transport were all simulated well, based on statistical and graphical model performance evaluation. The InHM simulations reported herein indicate that the 3D geometry and hydraulic characteristics of the layered geologic interfaces at CB1 can control the development of saturation and pore-water pressures at the soil-saprolite interface. The weathered bedrock piezometric response and runoff contribution were not simulated well with InHM in this study, most likely as a result of the uncertainty in the weathered bedrock layer geometry and fractured-rock hydraulic properties that preclude accurate fracture flow representation. Sensitivity analyses for the CB1 boundary-value problem indicate that: (i) hysteretic unsaturated flow in the CB1 soil is important for accurate hydrologic-response simulation, (ii) using an impermeable boundary condition to represent layered geologic interfaces leads to large errors in simulated magnitudes of runoff generation and pore-water pressure development, and (iii) field-based retention curve measurements can dramatically improve variably-saturated hydrologic-response simulation at sites with steep soil-water retention curves. The near-surface CB1 simulations reported herein demonstrate that physics-based models like InHM are useful for characterizing detailed spatio-temporal hydrologic-response, developing process-based concepts, and identifying information shortfalls for the next generation of field experiments. The field-based observations and hydrologic-response simulations from CB1 highlight the challenges in characterizing/simulating fractured bedrock flow at small catchments, which has important consequences for hydrologic response and landslide initiation.

Near-surface hydrologic response for a steep, unchanneled catchment near Coos Bay, Oregon: 1. sprinkling experiments

Sprinkling systems are frequently used to simulate rainfall for process-based investigations of near-surface hydrologic response without measuring or accounting for spatial variability. Data analyses from three sprinkling experiments at the Coos Bay 1 experimental catchment (CB1) demonstrate considerable spatial variability in sprinkling. Furthermore, simulated rainfall from sprinklers was found to be more heterogeneous than natural storms at CB1. Water balance calculations and evapotranspiration estimates indicate that evaporation of airborne droplets is a significant portion of applied sprinkling rates, although still less than the amount blown off the field site by strong winds. Incorporation of spatial variability in sprinkling input and soil-water storage did not significantly change water balance calculations. Saturation patterns within the near-surface soil profile and the timing of tensiometric response are affected by sprinkling heterogeneity. Pore-water pressure and saturation development at the soil-saprolite interface are primarily controlled by convergent surface / subsurface topography and bedrock fracture flow, but are also sensitive to sprinkling spatial variations. The analyses presented herein suggest that incorporating spatial variability in sprinkling rates is important when conducting hydrologic-response modeling of sprinkler experiments. This paper is the first-part of a two-part series focused on CB1. The data analyses in this paper are used to parameterize comprehensive physics-based hydrologic-response simulations of three CB1 sprinkling experiments reported in the companion paper.

Assessing the detail needed to capture rainfall-runoff dynamics with physics-based hydrologic response simulation

[1] Concept development simulation with distributed, physics‐based models provides a quantitative approach for investigating runoff generation processes across environmental conditions. Disparities within data sets employed to design and parameterize boundary value problems used in heuristic simulation inevitably introduce various levels of bias. The objective was to evaluate the impact of boundary value problem complexity on process representation for different runoff generation mechanisms. The comprehensive physics‐based hydrologic response model InHM has been employed to generate base case simulations for four well‐characterized catchments. The C3 and CB catchments are located within steep, forested environments dominated by subsurface stormflow; the TW and R5 catchments are located in gently sloping rangeland environments dominated by Dunne and Horton overland flows. Observational details are well captured within all four of the base case simulations, but the characterization of soil depth, permeability, rainfall intensity, and evapotranspiration differs for each. These differences are investigated through the conversion of each base case into a reduced case scenario, all sharing the same level of complexity. Evaluation of how individual boundary value problem characteristics impact simulated runoff generation processes is facilitated by quantitative analysis of integrated and distributed responses at high spatial and temporal resolution. Generally, the base case reduction causes moderate changes in discharge and runoff patterns, with the dominant process remaining unchanged. Moderate differences between the base and reduced cases highlight the importance of detailed field observations for parameterizing and evaluating physics‐based models. Overall, similarities between the base and reduced cases indicate that the simpler boundary value problems may be useful for concept development simulation to investigate fundamental controls on the spectrum of runoff generation mechanisms.

Relations between soil hydraulic properties and burn severity

Wildfire can affect soil hydraulic properties, often resulting in reduced infiltration. The magnitude of change in infiltration varies depending on the burn severity. Quantitative approaches to link burn severity with changes in infiltration are lacking. This study uses controlled laboratory measurements to determine relations between a remotely sensed burn severity metric (dNBR, change in normalised burn ratio) and soil hydraulic properties (SHPs). SHPs were measured on soil cores collected from an area burned by the 2013 Black Forest fire in Colorado, USA. Six sites with the same soil type were selected across a range of burn severities, and 10 random soil cores were collected from each site within a 30-m diameter circle. Cumulative infiltration measurements were made in the laboratory using a tension infiltrometer to determine field-saturated hydraulic conductivity, Kfs, and sorptivity, S. These measurements were correlated with dNBR for values ranging from 124 (low severity) to 886 (high severity). SHPs were related to dNBR by inverse functions for specific conditions of water repellency (at the time of sampling) and soil texture. Both functions had a threshold value for dNBR between 124 and 420, where Kfs and S were unchanged and equal to values for soil unaffected by fire. For dNBRs >~420, the Kfs was an exponentially decreasing function of dNBR and S was a linearly decreasing function of dNBR. These initial quantitative empirical relations provide a first step to link SHPs to burn severity, and can be used in quantitative infiltration models to predict post-wildfire infiltration and resulting runoff.

Aspect-dependent soil saturation and insight into debris-flow initiation during extreme rainfall in the Colorado Front Range

Hydrologic processes during extreme rainfall events are poorly characterized because of the rarity of measurements. Improved understanding of hydrologic controls on natural hazards is needed because of the potential for substantial risk during extreme precipitation events. We present field measurements of the degree of soil saturation and estimates of available soil-water storage during the September 2013 Colorado extreme rainfall event at burned (wildfire in 2010) and unburned hillslopes with north- and south-facing slope aspects. Soil saturation was more strongly correlated with slope aspect than with recent fire history; south-facing hillslopes became fully saturated while north-facing hillslopes did not. Our results suggest multiple explanations for why aspect-dependent hydrologic controls favor saturation development on south-facing slopes, causing reductions in effective stress and triggering of slope failures during extreme rainfall. Aspect-dependent hydrologic behavior may result from (1) a larger gravel and stone fraction, and hence lower soil-water storage capacity, on south-facing slopes, and (2) lower weathered-bedrock permeability on south-facing slopes, because of lower tree density and associated deep roots penetrating bedrock as well as less intense weathering, inhibiting soil drainage.

Hydrologic and Isotopic Modeling of Alpine Lake Waiau, Mauna Kea, Hawai'i

ABSTRACT Analysis of hydrologic, meteorologic, and isotopic data collected over 3 yr quantifies and explains the enormous variability and isotopic enrichment (δ18O = +16.9, δD = +50.0) of alpine Lake Waiau, a culturally and ecologically significant perched lake near the summit of Mauna Kea, Hawai‘i. Further, a simple one-dimensional hydrologic model was developed that couples standard water budget modeling with modeling of δD and δ18O isotopic composition to provide daily predictions of lake volume and chemistry. Data analysis and modeling show that winter storms are the primary source of water for the lake, adding a distinctively light isotopic signature appropriate for high-altitude precipitation. Evaporation at the windy, dry summit is the primary loss mechanism for most of the year, greatly enriching the lake in heavy isotopes.

Meta‐analysis of field‐saturated hydraulic conductivity recovery following wildland fire: Applications for hydrologic model parameterization and resilience assessment

Hydrologic recovery after wildfire is critical for restoring the ecosystem services of protecting of human lives and infrastructure from hazards and delivering water supply of sufficient quality and quantity. Recovery of soil-hydraulic properties, such as field-saturated hydraulic conductivity (Kfs), is a key factor for assessing the duration of watershed-scale flash flood and debris flow risks after wildfire. Despite the crucial role of Kfs in parameterizing numerical hydrologic models to predict the magnitude of postwildfire run-off and erosion, existing quantitative relations to predict Kfs recovery with time since wildfire are lacking. Here, we conduct meta-analyses of 5 datasets from the literature that measure or estimate Kfs with time since wildfire for longer than 3-year duration. The meta-analyses focus on fitting 2 quantitative relations (linear and non-linear logistic) to explain trends in Kfs temporal recovery. The 2 relations adequately described temporal recovery except for 1 site where macropore flow dominated infiltration and Kfs recovery. This work also suggests that Kfs can have low hydrologic resistance (large postfire changes), and moderate to high hydrologic stability (recovery time relative to disturbance recurrence interval) and resilience (recovery of hydrologic function and provision of ecosystem services). Future Kfs relations could more explicitly incorporate processes such as soil-water repellency, ground cover and soil structure regeneration, macropore recovery, and vegetation regrowth.

Disturbance Hydrology: Preparing for an Increasingly Disturbed Future

This special issue is the result of several fruitful conference sessions on disturbance hydrology, which started at the 2013 AGU Fall Meeting in San Francisco and have continued every year since. The stimulating presentations and discussions surrounding those sessions have focused on understanding both the disruption of hydrologic functioning following discrete disturbances, as well as the subsequent recovery or change within the affected watershed system. Whereas some hydrologic disturbances are directly linked to anthropogenic activities, such as resource extraction, the contributions to this special issue focus primarily on those with indirect or less pronounced human involvement, such as bark-beetle infestation, wildfire, and other natural hazards. However, human activities are enhancing the severity and frequency of these seemingly natural disturbances, thereby contributing to acute hydrologic problems and hazards. Major research challenges for our increasingly disturbed planet include the lack of continuous pre and postdisturbance monitoring, hydrologic impacts that vary spatially and temporally based on environmental and hydroclimatic conditions, and the preponderance of overlapping or compounding disturbance sequences. In addition, a conceptual framework for characterizing commonalities and differences among hydrologic disturbances is still in its infancy. In this introduction to the special issue, we advance the fusion of concepts and terminology from ecology and hydrology to begin filling this gap. We briefly explore some preliminary approaches for comparing different disturbances and their hydrologic impacts, which provides a starting point for further dialogue and research progress.

Hydraulic Property and Soil Textural Classification Measurements for Rainier Mesa, Nevada Test Site, Nevada

This report presents particle size analysis, field-saturated hydraulic conductivity measurements, and qualitative descriptions of surficial materials at selected locations at Rainier Mesa, Nevada. Measurements and sample collection were conducted in the Rainier Mesa area, including unconsolidated sediments on top of the mesa, an ephemeral wash channel near the mesa edge, and dry U12n tunnel pond sediments below the mesa. Particle size analysis used a combination of sieving and optical diffraction techniques. Field-saturated hydraulic conductivity measurements employed a single-ring infiltrometer with analytical formulas that correct for falling head and spreading outside the ring domain. These measurements may prove useful to current and future efforts at Rainier Mesa aimed at understanding infiltration and its effect on water fluxes and radionuclide transport in the unsaturated zone.

Postwildfire measurement of soil physical and hydraulic properties at selected sampling sites in the 2011 Las Conchas wildfire burn scar, Jemez Mountains, north-central New Mexico

..........................................................................................................................................................