Search USGS⌕ Search

Geology topics

Francis K. Rengers

Publications and source records attributed to Francis K. Rengers.

72 records · Page 4Linked to original sources

Incorporating spatially heterogeneous infiltration capacity into hydrologic models with applications for simulating post‐wildfire 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 s ) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of K s , and the spatial correlation length scale associated with K s . The proposed method for determining 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 s measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on 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.

California↗

Overview of the geologic effects of the November 14, 2016, Mw 7.8 Kaikoura, New Zealand, earthquake

The November 14, 2016, Kaikoura, New Zealand, earthquake (moment magnitude [Mw] 7.8) triggered more than 10,000 landslides over an area of about 12,000 square kilometers in the northeastern part of the South Island of New Zealand. In collaboration with GNS Science (the Institute of Geological and Nuclear Science Limited), we conducted ground and helicopter reconnaissance of the affected areas and assisted in rapid hazard evaluation. The majority of the triggered landslides were shallow- to moderate-depth (1–10 meters), highly disrupted falls and slides in rock and debris from Lower Cretaceous graywacke sandstone in the Seaward Kaikoura Range. Deeper, more coherent landslides in weak Upper Cretaceous to Neogene sedimentary rock also were numerous in the gentler topography south and inland (west) of the Seaward Kaikoura Range. The principal ground-failure hazards from the earthquake were the hundreds of valley-blocking landslides, many of which impounded lakes and ponds that posed potential downstream flooding hazards. Both large and small landslides also blocked road and rail corridors in many locations, including the main north-south highway (State Highway 1), which was still closed in October 2017. As part of our investigation, we compared post-earthquake field observations to the output of models used to estimate near-real-time landslide probabilities following earthquakes. The models generally over-predicted landslide occurrence and thus need further refinement.

Kaikoura↗

Improving near‐real‐time coseismic landslide models: Lessons learned from the 2016 Kaikōura, New Zealand, earthquake

The U.S. Geological Survey (USGS) is developing near‐real‐time global earthquake‐triggered‐landslide products to augment the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) system. The 14 November 2016 M w "> M w Mw 7.8 Kaikōura, New Zealand, earthquake provided a test case for evaluating the performance and near‐real‐time response applicability of three published global seismically induced landslide models. All three models obtain shaking estimates from the USGS ShakeMap, which is updated and sometimes changes significantly in the hours to days after an earthquake. The Kaikōura earthquake is a particularly valuable event that helps us better understand how changes to the ShakeMap affect the landslide models because the ShakeMap evolved significantly over several weeks as multifault rupture and seismic data were incorporated. We used the detailed landslide inventory available for this event for qualitative landslide model evaluation. We found that once a point source was replaced with an approximate rupture extent in ShakeMap, the landslide models were all successful at roughly identifying the area of highest hazard. This is notable, given that the models are relatively simple, coarse in resolution, and are based solely on input proxies that are globally available. However, all of the models dramatically overpredicted the hazard level, which indicates that improvements can be made. Subsequent updates to the ShakeMap resulted in improvements to model performance by some metrics and declining performance by others. In all cases, details of the ShakeMap strongly controlled the spatial pattern, even when those details were erroneous, such as the inclusion of a fault segment that did not rupture. If maps of landslide hazard are to be used effectively for rapid response, then we need to understand and clearly communicate the control that ShakeMap has over the models and how that typically evolves with time and is (or is not) reflected in reported uncertainties.

Bulletin of the Seismological Society of America↗

Landslides triggered by the 14 November 2016 Mw 7.8 Kaikōura Earthquake, New Zealand

The 14 November 2016 M w "> M w Mw 7.8 Kaikōura earthquake generated more than 10,000 landslides over a total area of about 10 , 000    km 2 "> 10 , 000 km 2 10,000 km2 , with the majority concentrated in a smaller area of about 3600    km 2 "> 3600 km 2 3600 km2 . The largest landslide triggered by the earthquake had an approximate volume of 20 ( ± 2 )    M   m 3 "> 20 ( ± 2 ) M m 3 20(±2) M m3 , with a runout distance of about 2.7 km, forming a dam on the Hapuku River. In this article, we present version 1.0 of the landslide inventory we have created for this event. We use the inventory presented in this article to identify and discuss some of the controls on the spatial distribution of landslides triggered by the Kaikōura earthquake. Our main findings are (1) the number of medium to large landslides (source area ≥ 10 , 000    m 2 "> ≥ 10 , 000 m 2 ≥10,000 m2 ) triggered by the Kaikōura earthquake is smaller than for similar‐sized landslides triggered by similar magnitude earthquakes in New Zealand; (2) seven of the largest eight landslides (from 5 to 20    M   m 3 "> 20 M m 3 20 M m3 ) occurred on faults that ruptured to the surface during the earthquake; (3) the average landslide density within 200 m of a mapped surface fault rupture is three times that at a distance of 2500 m or more from a mapped surface fault rupture; (4) the “distance to fault” predictor variable, when used as a proxy for ground‐motion intensity, and when combined with slope angle, geology, and elevation variables, has more power in predicting landslide probability than the modeled peak ground acceleration or peak ground velocity; and (5) for the same slope angles, the coastal slopes have landslide point densities that are an order of magnitude greater than those in similar materials on the inland slopes, but their source areas are significantly smaller.

Bulletin of the Seismological Society of America↗

The evolution of a colluvial hollow to a fluvial channel with periodic steps following two transformational disturbances: A wildfire and a historic flood

The transition of a colluvial hollow to a fluvial channel with discrete steps was observed after two landscape-scale disturbances. The first disturbance, a high-severity wildfire, changed the catchment hydrology to favor overland flow, which incised a colluvial hollow, creating a channel in the same location. This incised channel became armored with cobbles and boulders following repeated post-wildfire overland flow events. Three years after the fire, a record rainstorm produced regional flooding and generated sufficient fluvial erosion and sorting to produce a fluvial channel with periodically spaced steps. An analysis of the step spacing shows that after the flood, newly formed steps retained a similar spacing to the topographic roughness spacing in the original colluvial hollow (prior to channelization). This suggests that despite a distinct change in channel form roughness and bedform morphology, the endogenous roughness periodicity was conserved. Variations in sediment erodibility helped to create the emergent steps as the largest particles ( >D84) remained immobile, becoming step features, and downstream soil was easily winnowed away.

Geomorphology↗

Estimating post-fire debris-flow hazards prior to wildfire using a statistical analysis of historical distributions of fire severity from remote sensing data

Following wildfire, mountainous areas of the western United States are susceptible to debris flow during intense rainfall. Convective storms that can generate debris flows in recently burned areas may occur during or immediately after the wildfire, leaving insufficient time for development and implementation of risk mitigation strategies. We present a method for estimating post-fire debris-flow hazards prior to wildfire using historical data to define the range of potential fire severities for a given location based on the statistical distribution of severity metrics obtained from remote sensing. Estimates of debris-flow likelihood, magnitude, and triggering rainfall threshold based upon the statistically simulated fire severity data provide hazard predictions consistent with those calculated from fire severity data collected after wildfire. Simulated fire severity data also produce hazard estimates that replicate observed debris-flow occurrence, rainfall conditions, and magnitude at a monitored site in the San Gabriel Mountains of southern California. Future applications of this method should rely upon a range of potential fire severity scenarios for improved pre-fire estimates of debris-flow hazard. The method presented here is also applicable to modeling other post-fire hazards, such as flooding and erosion risk, and for quantifying trends in observed fire severity in a changing climate.

California, Colorado, Washington↗

Quantifying post-wildfire hillslope erosion with lidar

Following a wildfire, flooding and debris- flow hazards are common and pose a threat to human life and infrastructure in steep burned terrain. Wildfire enhances both water runoff and soil erosion, which ultimately shape the debris flow potential. The erosional processes that route excess sediment from hillslopes to debris-flow channels in recently burned areas, however, are poorly constrained. In this study we examined erosional processes through repeat terrestrial lidar surveys in a steep mountainous watershed that experienced a high-severity burn in the 2016 San Gabriel complex fire. Three lidar surveys were conducted during a wet winter (2016-2017) on a hillslope plot. We used geomorphometric techniques to better contextualize erosion observations in areas with rills and between rills (interrill areas). A challenge was effectively differentiating DEM pixels that were in the constantly evolving rill network as well as those outside the rill network. By applying a series of DEM filtering processes we found that it was possible to efficiently identify the small-scale rill networks. Our results challenge previously held beliefs about sediment erosion on burned hillslopes, suggesting that prior estimates made without access to high resolution topography likely underestimated the role of interrill erosion.

Conference Paper↗

Debris flow initiation by runoff in a recently burned basin: Is grain-by-grain sediment bulking or en masse failure to blame?

Postwildfire debris flows are frequently triggered by runoff following high-intensity rainfall, but the physical mechanisms by which water-dominated flows transition to debris flows are poorly understood relative to debris flow initiation from shallow landslides. In this study, we combined a numerical model with high-resolution hydrologic and geomorphic data sets to test two different hypotheses for debris flow initiation during a rainfall event that produced numerous debris flows within a recently burned drainage basin. Based on simulations, large volumes of sediment eroded from the hillslopes were redeposited within the channel network throughout the storm, leading to the initiation of numerous debris flows as a result of the mass failure of sediment dams that built up within the channel. More generally, results provide a quantitative framework for assessing the potential of runoff-generated debris flows based on sediment supply and hydrologic conditions.

California↗

Can thermoluminescence be used to determine soil heating from a wildfire?

The Silverado wildfire occurred from September 12 to 20, 2014, burning 960 acres in Orange County, California. Soil samples from within the burn area were obtained and the thermoluminescence (TL) properties of those samples were compared against a control sample to understand wildfire heating. We performed a series of experiments investigating the degree to which the control differed from the wildfire soil samples. This work showed that soil heated by a wildfire had a distinctly different glow curve shape than the unburned soil sample. Moreover, it was possible to see changes in the TL signal as a function of soil depth in wildfire-heated samples. Our experiments suggest that minimal soil heating occurred below approximately 10 cm. Estimates of wildfire temperatures, however, were nuanced.

California↗

Observed and simulated hydrologic response for a first-order catchment during extreme rainfall 3 years after wildfire disturbance

Hydrologic response to extreme rainfall in disturbed landscapes is poorly understood because of the paucity of measurements. A unique opportunity presented itself when extreme rainfall in September 2013 fell on a headwater catchment (i.e., <1 ha) in Colorado, USA that had previously been burned by a wildfire in 2010. We compared measurements of soil-hydraulic properties, soil saturation from subsurface sensors, and estimated peak runoff during the extreme rainfall with numerical simulations of runoff generation and subsurface hydrologic response during this event. The simulations were used to explore differences in runoff generation between the wildfire-affected headwater catchment, a simulated unburned case, and for uniform versus spatially variable parameterizations of soil-hydraulic properties that affect infiltration and runoff generation in burned landscapes. Despite 3 years of elapsed time since the 2010 wildfire, observations and simulations pointed to substantial surface runoff generation in the wildfire-affected headwater catchment by the infiltration-excess mechanism while no surface runoff was generated in the unburned case. The surface runoff generation was the result of incomplete recovery of soil-hydraulic properties in the burned area, suggesting recovery takes longer than 3 years. Moreover, spatially variable soil-hydraulic property parameterizations produced longer duration but lower peak-flow infiltration-excess runoff, compared to uniform parameterization, which may have important hillslope sediment export and geomorphologic implications during long duration, extreme rainfall. The majority of the simulated surface runoff in the spatially variable cases came from connected near-channel contributing areas, which was a substantially smaller contributing area than the uniform simulations.

Water Resources Research↗

Constraining the relative importance of raindrop- and flow-driven sediment transport mechanisms in postwildfire environments and implications for recovery time scales

Mountain watersheds recently burned by wildfire often experience greater amounts of runoff and increased rates of sediment transport relative to similar unburned areas. Given the sedimentation and debris flow threats caused by increases in erosion, more work is needed to better understand the physical mechanisms responsible for the observed increase in sediment transport in burned environments and the time scale over which a heightened geomorphic response can be expected. In this study, we quantified the relative importance of different hillslope erosion mechanisms during two postwildfire rainstorms at a drainage basin in Southern California by combining terrestrial laser scanner-derived maps of topographic change, field measurements, and numerical modeling of overland flow and sediment transport. Numerous debris flows were initiated by runoff at our study area during a long-duration storm of relatively modest intensity. Despite the presence of a well-developed rill network, numerical model results suggest that the majority of eroded hillslope sediment during this long-duration rainstorm was transported by raindrop-induced sediment transport processes, highlighting the importance of raindrop-driven processes in supplying channels with potential debris flow material. We also used the numerical model to explore relationships between postwildfire storm characteristics, vegetation cover, soil infiltration capacity, and the total volume of eroded sediment from a synthetic hillslope for different end-member erosion regimes. This study adds to our understanding of sediment transport in steep, postwildfire landscapes and shows how data from field monitoring can be combined with numerical modeling of sediment transport to isolate the processes leading to increased erosion in burned areas.

California↗

The influence of vegetation cover on debris-flow density during an extreme rainfall in the northern Colorado Front Range

We explored regional influences on debris-flow initiation throughout the Colorado Front Range (Colorado, USA) by exploiting a unique data set of more than 1100 debris flows that initiated during a 5 day rainstorm in 2013. Using geospatial data, we examined the influence of rain, hillslope angle, hillslope aspect, and vegetation density on debris-flow initiation. In particular we used a greenness index to differentiate areas of high tree density from grass and bare soil. The data demonstrated an overwhelming propensity for debris-flow initiation on south-facing hillslopes. However, when the debris-flow density was analyzed with respect to total rainfall and greenness we found that most debris flows occurred in areas of high rainfall and low tree density, regardless of hillslope aspect. These results indicate that present-day tree density exerts a stronger influence on debris-flow initiation locations than aspect-driven variations in soil and bedrock properties that developed over longer time scales.

Geology↗

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 9 to 13 September 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.

Geophysical Research Letters↗

Amplification of postwildfire peak flow by debris

In burned steeplands, the peak depth and discharge of postwildfire runoff can substantially increase from the addition of debris. Yet methods to estimate the increase over water flow are lacking. We quantified the potential amplification of peak stage and discharge using video observations of postwildfire runoff, compiled data on postwildfire peak flow ( Q p ), and a physically based model. Comparison of flood and debris flow data with similar distributions in drainage area ( A ) and rainfall intensity ( I ) showed that the median runoff coefficient ( C = Q p / AI ) of debris flows is 50 times greater than that of floods. The striking increase in Q p can be explained using a fully predictive model that describes the additional flow resistance caused by the emergence of coarse-grained surge fronts. The model provides estimates of the amplification of peak depth, discharge, and shear stress needed for assessing postwildfire hazards and constraining models of bedrock incision.

Geophysical Research Letters↗

Model simulations of flood and debris flow timing in steep catchments after wildfire

Debris flows are a typical hazard on steep slopes after wildfire, but unlike debris flows that mobilize from landslides, most post-wildfire debris flows are generated from water runoff. The majority of existing debris-flow modeling has focused on landslide-triggered debris flows. In this study we explore the potential for using process-based rainfall-runoff models to simulate the timing of water flow and runoff-generated debris flows in recently burned areas. Two different spatially distributed hydrologic models with differing levels of complexity were used: the full shallow water equations and the kinematic wave approximation. Model parameter values were calibrated in two different watersheds, spanning two orders of magnitude in drainage area. These watersheds were affected by the 2009 Station Fire in the San Gabriel Mountains, CA, USA. Input data for the numerical models were constrained by time series of soil moisture, flow stage, and rainfall collected at field sites, as well as high-resolution lidar-derived digital elevation models. The calibrated parameters were used to model a third watershed in the burn area, and the results show a good match with observed timing of flow peaks. The calibrated roughness parameter (Manning's $n$) was generally higher when using the kinematic wave approximation relative to the shallow water equations, and decreased with increasing spatial scale. The calibrated effective watershed hydraulic conductivity was low for both models, even for storms occurring several months after the fire, suggesting that wildfire-induced changes to soil-water infiltration were retained throughout that time. Overall the two model simulations were quite similar suggesting that a kinematic wave model, which is simpler and more computationally efficient, is a suitable approach for predicting flood and debris flow timing in steep, burned watersheds.

Water Resources Research↗

Episodic bedrock erosion by gully-head migration, Colorado High Plains, USA

This study explores the frequency of bedrock exposure in a soil-mantled low-relief (i.e. non-mountainous) landscape. In the High Plains of eastern Colorado, gully headcuts are among the few erosional features that will incise through the soil mantle to expose bedrock. We measured the last time of bedrock exposure using optically stimulated luminescence dating of alluvial sediment overlying bedrock in gully headcuts. Our dating suggests that headcuts in adjacent gullies expose bedrock asynchronously, and therefore, the headcuts are unlikely to have been triggered by a base-level drop in the trunk stream. This finding supports the hypothesis that headcuts can develop locally in gullies as a result of focused scour in locations where hydraulic stress during a flash flood is sufficiently high, and/or ground cover is sufficiently weak, to generate a scour hole that undermines vegetation. Alluvium dating also reveals that gullies have been a persistent part of this landscape since the early Holocene.

Colorado↗

Illuminating wildfire erosion and deposition patterns with repeat terrestrial lidar

Erosion following a wildfire is much greater than background erosion in forests because of wildfire-induced changes to soil erodibility and water infiltration. While many previous studies have documented post-wildfire erosion with point and small plot-scale measurements, the spatial distribution of post-fire erosion patterns at the watershed scale remains largely unexplored. In this study lidar surveys were collected periodically in a small, first-order drainage basin over a period of 2 years following a wildfire. The study site was relatively steep with slopes ranging from 17° to > 30°. During the study period, several different types of rain storms occurred on the site including low-intensity frontal storms (2.4 mm h −1 ) and high-intensity convective thunderstorms (79 mm h −1 ). These storms were the dominant drivers of erosion. Erosion resulting from dry ravel and debris flows was notably absent at the site. Successive lidar surveys were subtracted from one another to obtain digital maps of topographic change between surveys. The results show an evolution in geomorphic response, such that the erosional response after rain storms was strongly influenced by the previous erosional events and pre-fire site morphology. Hillslope and channel roughness increased over time, and the watershed armored as coarse cobbles and boulders were exposed. The erosional response was spatially nonuniform; shallow erosion from hillslopes (87% of the study area) contributed 3 times more sediment volume than erosion from convergent areas (13% of the study area). However, the total normalized erosion depth (volume/area) was highest in convergent areas. From a detailed understanding of the spatial locations of erosion, we made inferences regarding the processes driving erosion. It appears that hillslope erosion is controlled by rain splash (for detachment) and overland flow (for transport and quasi-channelized erosion), with the sites of highest erosion corresponding to locations with the lowest roughness. By contrast, in convergent areas we found erosion caused by overland flow. Soil erosion was locally interrupted by immobile objects such as boulders, bedrock, or tree trunks, resulting in a patchy erosion network with increasing roughness over time.

Colorado↗

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.

Geology↗