Benjamin B. Mirus

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.

Elucidating the role of vegetation in the initiation of rainfall-induced shallow landslides: Insights from an extreme rainfall event in the Colorado Front Range

More than 1100 debris flows were mobilized from shallow landslides during a rainstorm from September 9-13, 2013 in the Colorado Front Range, with the vast majority initiating on sparsely vegetated, south-facing terrain. To investigate the physical processes responsible for the observed aspect control, we made measurements of soil properties on a densely forested north-facing hillslope and a grassland-dominated south-facing hillslope in the Colorado Front Range and performed numerical modeling of transient changes in soil pore water pressure throughout the rainstorm. Using the numerical model, we quantitatively assessed interactions among vegetation, rainfall interception, subsurface hydrology, and slope stability. Results suggest that apparent cohesion supplied by roots was responsible for the observed connection between debris flow initiation and slope aspect. Results suggest that future climate-driven modifications to forest structure could substantially influence landslide hazards throughout the Front Range and similar water-limited environments where vegetation communities may be more susceptible to small variations in climate.

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.

Advanced methods for modeling water-levels and estimating drawdowns with SeriesSEE, an Excel add-in

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

Hydrologic Impacts of Landslide Disturbances: Implications for Remobilization and Hazard Persistence

Landslides typically alter hillslope topography, but may also change the hydrologic connectivity and subsurface water-storage dynamics. In settings where mobile materials are not completely evacuated from steep slopes, influences of landslide disturbances on hillslope hydrology and susceptibility to subsequent failures remain poorly characterized. Since landslides often recur at the site of previous failures, we examine differences between a stable vegetated hillslope (VH) and a recent landslide (LS). These neighboring hillslopes exhibit similar topography and are situated on steep landslide-prone coastal bluffs of glacial deposits along the northeastern shore of Puget Sound, Washington. Our control hillslope, VH, is mantled by a heterogeneous colluvium, supporting a dense forest. In early 2013, our test hillslope, LS, also supported a forest before a landslide substantially altered the topography and disturbed the hillslope. In 2015, we observed a clay-rich landslide deposit at LS with sparse vegetation and limited root reinforcement, soil structures, and macropores. Our characterization of the sites also found matrix porosity and hydraulic conductivity are both lower at LS. Continuous monitoring during 2015–2016 revealed reduced effective precipitation at VH (due to canopy interception), an earlier seasonal transition to near-saturated conditions at LS, and longer persistence of positive pore pressures and slower drainage at LS (both seasonally and between major storm events). These differences, along with episodic, complex slope failures at LS support the hypothesis that, despite a reduced average slope, other disturbances introduced by landsliding may promote the hydrologic conditions leading to slope instability, thus contributing to the persistence of landslide hazards.

Integrating real-time subsurface hydrologic monitoring with empirical rainfall thresholds to improve landslide early warning

Early warning for rainfall-induced shallow landsliding can help reduce fatalities and economic losses. Although these commonly occurring landslides are typically triggered by subsurface hydrological processes, most early warning criteria rely exclusively on empirical rainfall thresholds and other indirect proxies for subsurface wetness. We explore the utility of explicitly accounting for antecedent wetness by integrating real-time subsurface hydrologic measurements into landslide early warning criteria. Our efforts build on previous progress with rainfall thresholds, monitoring, and numerical modeling along the landslide-prone railway corridor between Everett and Seattle, Washington, USA. We propose a modification to a previously established recent versus antecedent (RA) cumulative rainfall thresholds by replacing the antecedent 15-day rainfall component with an average saturation observed over the same timeframe. We calculate this antecedent saturation with real-time telemetered measurements from five volumetric water content probes installed in the shallow subsurface within a steep vegetated hillslope. Our hybrid rainfall versus saturation (RS) threshold still relies on the same recent 3-day rainfall component as the existing RA thresholds, to facilitate ready integration with quantitative precipitation forecasts. During the 2015–2017 monitoring period, this RS hybrid approach has an increase of true positives and a decrease of false positives and false negatives relative to the previous RA rainfall-only thresholds. We also demonstrate that alternative hybrid threshold formats could be even more accurate, which suggests that further development and testing during future landslide seasons is needed. The positive results confirm that accounting for antecedent wetness conditions with direct subsurface hydrologic measurements can improve thresholds for alert systems and early warning of rainfall-induced shallow landsliding.

Variability in soil-water retention properties and implications for physics-based simulation of landslide early warning criteria

Rainfall-induced shallow landsliding is a persistent hazard to human life and property. Despite the observed connection between infiltration through the unsaturated zone and shallow landslide initiation, there is considerable uncertainty in how estimates of unsaturated soil-water retention properties affect slope stability assessment. This source of uncertainty is critical to evaluating the utility of physics-based hydrologic modeling as a tool for landslide early warning. We employ a numerical model of variably saturated groundwater flow parameterized with an ensemble of texture-, laboratory-, and field-based estimates of soil-water retention properties for an extensively monitored landslide-prone site in the San Francisco Bay Area, CA, USA. Simulations of soil-water content, pore-water pressure, and the resultant factor of safety show considerable variability across and within these different parameter estimation techniques. In particular, we demonstrate that with the same permeability structure imposed across all simulations, the variability in soil-water retention properties strongly influences predictions of positive pore-water pressure coincident with widespread shallow landsliding. We also find that the ensemble of soil-water retention properties imposes an order-of-magnitude and nearly two-fold variability in seasonal and event-scale landslide susceptibility, respectively. Despite the reduced factor of safety uncertainty during wet conditions, parameters that control the dry end of the soil-water retention function markedly impact the ability of a hydrologic model to capture soil-water content dynamics observed in the field. These results suggest that variability in soil-water retention properties should be considered for objective physics-based simulation of landslide early warning criteria.

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.

Effect of Hydraulic Hysteresis on Stability of Infinite Slopes under Steady Infiltration

AbstractHydraulic hysteresis, including capillary soil water retention (SWR), air entrapment SWR, and hydraulic conductivity, is a common phenomenon in unsaturated soils. However, the influence of ...

Incorporating spatially heterogeneous infiltration capacity into hydrologic models with applications for simulating post-wildfire debris flow initiation: Post-wildfire runoff and debris flow initiation

Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity (Kₛ) to derive an expression for an aerially averaged saturated hydraulic conductivity (Ke*) that depends on the rainfall rate, the statistical properties of Kₛ, and the spatial correlation length scale associated with Kₛ. The proposed method for determining Ke* is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale Kₛ measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on Ke* to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.