Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico

Jonathan P. Perkins; Mong-Han Huang; Corina Cerovski-Darriau; William H. Schulz; Collin Cronkite-Ratcliff; Kelli Wadsworth Baxstrom; Mason Muir Einbund; Jane Willenbring; Ben Leshchinsky; Samuel Johnstone

doi:10.1029/2025jf008847

USGS · 70280631

Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico

Jonathan P. Perkins·Mong-Han Huang·Corina Cerovski-Darriau·William H. Schulz·Collin Cronkite-Ratcliff·Kelli Wadsworth Baxstrom·Mason Muir Einbund·Jane Willenbring·Ben Leshchinsky·Samuel Johnstone

View original source ↗Export metadata ↓

Abstract

Landscapes in humid climates often experience rapid chemical and mechanical weathering that can change groundwater fluxes, drainage pathways, and frictional strength via alterations in subsurface mineral composition and matrix porosity. These changes to hillslope critical zone structure and strength should control landslide susceptibility, size, and failure style. Here we test the impact of differential weathering on landslide failure processes by exploiting a natural experiment along the dioritic Utuado batholith in Puerto Rico's Cordillera Central, where a sharp weathering gradient follows an upstream-propagating knickzone into an uplifted shoreline platform that was affected by Hurricane Maria in 2017. Through a combination of landslide mapping, field geophysical and geotechnical surveys, topographic analysis, calibrated stream profile inversion, 3D back-analysis of hillslope strength, and debris-flow inundation assessments, we show that differential weathering induces changes to critical zone structure and consequently landslide susceptibility and hazard. Our mapping reveals that the hurricane-induced landslides preferentially cluster along less weathered hillslopes despite a marked increase in frictional strength. Jointly inverted seismic refraction and electrical resistivity profiles show that the frictionally weak, heavily weathered hillslopes are highly porous and dry, which we posit acts to stabilize slopes through limiting deeper hydrologic connectivity that can produce shallow water tables with destabilizing lateral seepage forces during storms. Debris-flow inundation mapping and analysis reveals that major contributors to debris-flow growth, channel steepness and landslide coalescence, scale linearly with hillslope friction angle. This leads to higher-mobility debris flows along knickzone reaches, which disproportionately affects communities at the base of the plateau.

Explore related subjects

JGR Earth Surface Article

Source-reported bounding extent: 18.243967032677233° to 18.32926547571553° latitude; -66.72485429203331° to -66.5607322183898° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Jonathan P. Perkins, Mong-Han Huang, Corina Cerovski-Darriau, William H. Schulz, Collin Cronkite-Ratcliff, Kelli Wadsworth Baxstrom, Mason Muir Einbund, Jane Willenbring, Ben Leshchinsky, Samuel Johnstone. 2026-09-18. Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico. https://doi.org/10.1029/2025jf008847

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2025

Modeling the influence of upper and lower shoreface dynamics on barrier island evolution

Barrier island resilience to climate impacts depends on sediment redistribution between the subaqueous shoreface and subaerial barrier during sea-level rise and storms. However, autogenic interactions between the upper and lower shoreface and their influence on the subaerial barrier are poorly characterized. Here, we explore the influences of various shoreface components on barrier morphology using a model of barrier and shoreface evolution under sea-level rise, the Articulated Barrier Shoreface (ABSF) Model. This reduced-complexity model divides the shoreface into upper and lower shoreface panels that respond independently to sea-level rise and deviations from the equilibrium slope. We couple the ABSF with the Lorenzo-Trueba & Ashton, 2014, https://doi.org/10.1002/2013jf002941 model (LTA), a barrier island evolution model driven by overwash and sea-level rise. Through this coupled framework, we examine the influences of upper and lower shoreface slopes, their respective depths, and sensitivity to wave climate on long-term barrier evolution. Results show that the relative depths of the upper and lower shoreface toes influence barrier response to rising seas, alongside overwash flux and closure depth. Notably, the lower shoreface response to sea-level change lags that of the upper shoreface over decades, diminishing the resilience of the barrier over centennial timescales by slowing the overall barrier response. In fact, the ABSF model predicts barriers will drown faster and more than predicted with a linear shoreface. Results highlight the shoreface as an important sediment reservoir for barrier islands and that differences in upper and lower shoreface responses can reduce barrier resilience to sea-level rise due to limited lower shoreface sediment accessibility.

JGR Earth Surface↗

USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗

USGS2025

Numerical simulation of sound-side barrier-island inundation and breaching during Hurricane Dorian (2019)

Hurricane-induced morphological changes and associated community hazards along sandy, barrier-island coastlines have been studied primarily from the perspective of ocean-side attack by storm-driven ocean surge and large waves. Thus, our understanding of long-term barrier island morphological change focuses on beach erosion, overwash, and inlet formation. In contrast, outwash events with inundation from the sound side, such as one that occurred in Cape Lookout National Seashore, North Carolina, USA during Hurricane Dorian (September 2019), are understudied. Studying such events can improve understanding of barrier island response and stability for a broader range of conditions. Here, we model the hydrodynamics and morphological evolution of a barrier island using a coupled wave-current-sediment transport modeling system. Wind-driven surge in Pamlico Sound led to overtopping from the sound side, which eroded outwash channels and transported sediment seaward into the nearshore. Simulations reproduce the channel features observed with aerial imagery and provide information not available from the remote-sensing observations, including channel depths (>2 m) and the fate of the eroded sand. We found that >99% of the eroded sand was deposited in the nearshore, within 1,000 m of the shoreline in depths <10 m, suggesting that the deposited sediment remains available for littoral transport and beach recovery. Simulations with combinations of coarse or fine sediment and vegetated or unvegetated landcover indicate that channel position did not vary with grain size or vegetation, while volume of erosion and channel morphology were more responsive to variations in grain size and less responsive to presence of vegetation.

North Carolina↗