Search for: All records

This Review synthesizes progress and outlines a new framework for understanding how land surface hazards interact and propagate as sediment cascades across Earth’s surface, influenced by interactions among the atmosphere, biosphere, hydrosphere, and solid Earth. Recent research highlights a gap in understanding these interactions on human timescales, given rapid climatic change and urban expansion into hazard-prone zones. We review how surface processes such as coseismic landslides and post-fire debris flows form a complex sequence of events that exacerbate hazard susceptibility. Moreover, innovations in modeling, remote sensing, and critical zone science can offer new opportunities for quantifying cascading hazards. Looking forward, societal resilience can increase by transforming our understanding of cascading hazards through advances in integrating data into comprehensive models that link across Earth systems. 

Measurements of cosmic-ray-produced beryllium-10, neon-21, and helium-3 in quartz in a soil profile from a forested landscape in the Oregon Coast Range show that the cosmogenic noble gases 21Ne and 3He are depleted relative to 10Be in the shallow subsurface. The noble gases are mobile in mineral grains via thermally activated diffusion and 10Be is not, implying that noble gas depletion is the result of surface heating by wildfires and subsequent mixing of partially degassed quartz downward into the soil. Cosmogenic noble gas depletion by wildfire heating of soils is a potential means of estimating wildfire intensity and/or frequency over pre-observational timescales. 

This data release contains two debris-flow inventories summarizing observations from burned and unburned areas in the western Cascade Range of Oregon (OR). The burned inventory focuses on debris flows that occurred during the first two years after the 2020 Archie Creek, Holiday Farm, Beachie Creek/Lionshead, and Riverside fires (OR_field_observations.csv). The unburned inventory (1995-2022) focuses on debris flows in the same areas (excluding the Riverside Fire). The inventories are derived from field observations (OR_field_observations.csv) and aerial imagery (OR_imagery_observations.csv). They include mapped debris-flow initiation locations, descriptions of the inferred initiation process, other notable site characteristics, and rainfall data. Locations of debris flows observed after wildfires are also linked to USGS postfire debris-flow hazard assessments (USGS, 2022; Staley and others, 2017; Thomas and others 2023). Rainfall characteristics for each debris flow in the inventory are derived from the closest rainfall gage to an observed debris flow (gage_locations.csv). Peak rainfall rates during the known time window of debris-flow initiation are reported for durations of 15 minutes, 30 minutes, 60 minutes, 12 hours, 24 hours, 36 hours, and 48 hours. More detailed explanations of the headers for each of these csv files can be found within the README_csvname.txt file. References: Landslide Hazards Program. (n.d.). Emergency assessment of post-fire debris-flow hazards. U.S. Geological Survey. https://landslides.usgs.gov/hazards/postfire_debrisflow Staley, D. M., Negri, J. A., Kean, J. W., Laber, J. L., Tillery, A. C., and Youberg, A. M., 2017, Prediction of spatially explicit rainfall intensity–duration thresholds for post-fire debris-flow generation in the western United States. Geomorphology, 278, 149–162. https://doi.org/10.1016/j.geomorph.2016.10.019 Thomas, M. A., Kean, J. W., McCoy, S. W., Lindsay, D. N., Kostelnik, J., Cavagnaro, D. B., Rengers, F. K., East, A. E., Schwartz, J. Y., Smith, D. P., and Collins, B. D., 2023, Postfire hydrologic response along the Central California (USA) coast: insights for the emergency assessment of postfire debris-flow hazards. Landslides, 20, 2421-2436. https://doi.org/10.1007/s10346-023-02106-7 

Abstract. Estimation of erosion rate is an important component of landscapeevolution studies, particularly in settings where transience or spatialvariability in uplift or erosion generates diverse landform morphologies.While bedrock rivers are often used to constrain the timing and magnitude of changes in baselevel lowering, hilltop curvature (or convexity), CHT, provides an additional opportunity to map variations in erosion rate given that average slope angle becomes insensitive to erosion rate owing to threshold slope processes. CHT measurement techniques applied in prior studies (e.g., polynomial functions), however, tend to be computationallyexpensive when they rely on high-resolution topographic data such as lidar,limiting the spatial extent of hillslope geomorphic studies to small studyregions. Alternative techniques such as spectral tools like continuouswavelet transforms present an opportunity to rapidly document trends inhilltop convexity across expansive areas. Here, we demonstrate howcontinuous wavelet transforms (CWTs) can be used to calculate the Laplacianof elevation, which we utilize to estimate erosion rate in three catchmentsof the Oregon Coast Range that exhibit varying slope angle, slope length,and hilltop convexity, implying differential erosion. We observe thatCHT values calculated with the CWT are similar to those obtained from2D polynomial functions. Consistent with recent studies, we find thaterosion rates estimated with CHT from both CWTs and 2D polynomialfunctions are consistent with erosion rates constrained with cosmogenicradionuclides from stream sediments. Importantly, our CWT approachcalculates curvature at least 103 times more quickly than 2Dpolynomials. This efficiency advantage of the CWT increases with domainsize. As such, continuous wavelet transforms provide a compelling approachto rapidly quantify regional variations in erosion rate as well aslithology, structure, and hillslope sediment transport processes, which areencoded in hillslope morphology. Finally, we test the accuracy of CWT and 2Dpolynomial techniques by constructing a series of synthetic hillslopesgenerated by a theoretical nonlinear transport model that exhibit a range oferosion rates and topographic noise characteristics. Notably, we find thatneither CWTs nor 2D polynomials reproduce the theoretically prescribedCHT value for hillslopes experiencing moderate to fast erosion rates,even when no topographic noise is added. Rather, CHT is systematicallyunderestimated, producing a power law relationship between erosion rate andCHT that can be attributed to the increasing prominence of planarhillslopes that narrow the zone of hilltop convexity as erosion rateincreases. As such, we recommend careful consideration of measurement lengthscale when applying CHT to estimate erosion rate in moderate tofast-eroding landscapes, where curvature measurement techniques may be prone to systematic underestimation. 

Abstract We deployed a network of 68 three-component geophones on the slow-moving Two Towers earthflow in northern California. We compute horizontal-to-vertical spectral ratios (HVSRs) from the ambient seismic field. The HVSRs have two prominent peaks, one near 1.23 Hz and another between 4 and 8 Hz at most stations. The 1.23 Hz resonance is a property of the background noise field and may be due to a velocity contrast at a few hundred meters depth. We interpret the higher frequency peaks as being related to slide deposits and invert the spectral ratios for shallow velocity structure using in situ thickness measurements as a priori constraints on the inversion. The thickness of the shallowest, low-velocity layer is systematically larger than landslide thicknesses inferred from inclinometer data acquired since 2013. Given constraints from field observations and boreholes, the inversion may reflect the thickness of deposits of an older slide that is larger in spatial extent and depth than the currently active slide. Because the HVSR peaks measured at Two Towers are caused by shallow slide deposits and represent frequencies that will experience amplification during earthquakes, the depth of the actively sliding mass may be less relevant for assessing potential slide volume and associated hazard than the thicknesses determined by our inversions. More generally, our results underscore the utility of combining both geotechnical measurements and subsurface imaging for landslide characterization and hazard assessment. 

Abstract Debris flows pose persistent hazards and shape high‐relief landscapes in diverse physiographic settings, but predicting the spatiotemporal occurrence of debris flows in postglacial topography remains challenging. To evaluate the debris flow process in high‐relief postglacial terrain, we conducted a geomorphic investigation to characterize geologic, glacial, volcanic, and land use contributions to landslide initiation across Southeast Alaska. To evaluate controls on landslide (esp. debris flow) occurrence in Sitka, we used field observation, geomorphic mapping, landslide characteristics as documented in the Tongass National Forest inventory, and a novel application of the shallow landslide model SHALSTAB to postglacial terrain. A complex geomorphic history of glaciation and volcanic activity provides a template for spatially heterogeneous landslide occurrence. Landslide density across the region is highly variable, but debris flow density is high on south‐ or southeast‐facing hillslopes where volcanic tephra soils are present and/or where timber harvest has occurred since 1900. High landslide density along the western coast of Baranof and Kruzof islands coincides with deposition of glacial sediment and thick tephra and exposure to extreme rainfall from atmospheric rivers on south‐facing aspects but the relative contributions of these controls are unclear. Timber harvest has also been identified as an important control on landslide occurrence in the region. Focusing on a subset of geo‐referenced landslides near Sitka, we used the SHALSTAB shallow landslide initiation model, which has been frequently applied in non‐glacial terrain, to identify areas of high landslide potential in steep, convergent terrain. In a validation against mapped landslide polygons, the model significantly outperformed random guessing, with area under the curve (AUC) = 0.709 on a performance classification curve of true positives vs. false positives. This successful application of SHALSTAB demonstrates practical utility for hazards analysis in postglacial landscapes to mitigate risk to people and infrastructure.