The structure of the critical zone is a product of feedbacks between hydrologic, climatic, biotic, and chemical processes. Ample research within snow-dominated systems has shown that aspect-dependent solar radiation inputs can produce striking differences in vegetation composition, topography, and soil depth between opposing hillslopes. However, more research is needed to understand the role of microclimates on critical zone development within rain-dominated systems, especially below the soil and into weathered bedrock. To address this need, we characterized the critical zone of a north-facing and south-facing slope within a first-order headwater catchment located in central coastal California. We combined terrain analysis of vegetation distribution and topography with field-based soil pit characterization, geophysical surveys and hydrologic measurements between slope-aspects. We observed thicker soil profiles, higher shallow soil moisture, and denser vegetation on north facing slopes, which matched previously documented observations in snow-dominated sites. However, average topographic gradient and saprolite thickness were uniform across our study hillslopes, which did not match common observations from the literature. These results suggest dominant processes for critical zone evolution are not necessarily transferable across regions. Thus, there is a continued need to expand critical zone research, especially in rain-dominated systems. Here, we present four non-exclusive, testable hypotheses of mechanisms that may explain these unexpected similarities in slope and saprolite thickness between hillslopes with opposing aspects. Specifically, we propose both past and present ecohydrologic processes must be taken into account to understand what shaped the present day critical zone.

In steep wildfire-burned terrains, intense rainfall can produce large volumes of runoff that can trigger highly destructive debris flows. The ability to accurately characterize and forecast debris-flow hazards in burned terrains, however, remains limited. Here, we augment the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) to simulate both overland and channelized flows and assess postfire debris-flow hazards over a regional domain. We perform hindcast simulations using high-resolution weather radar-derived precipitation and reanalysis data to drive non-burned baseline and burn scar sensitivity experiments. Our simulations focus on January 2021 when an atmospheric river triggered numerous debris flows within a wildfire burn scar in Big Sur – one of which destroyed California’s famous Highway 1. Compared to the baseline, our burn scar simulation yields dramatic increases in total and peak discharge, and shorter lags between rainfall onset and peak discharge. At Rat Creek, where Highway 1 was destroyed, discharge volume increases eight-fold and peak discharge triples relative to the baseline. For all catchments within the burn scar, we find that the median catchment-area normalized discharge volume increases nine-fold after incorporating burn scar characteristics, while the 95th percentile volume increases 13-fold. Catchments with anomalously high hazard levels correspond well with post-event debris flow observations. Our results demonstrate that WRF-Hydro provides a compelling new physics-based tool to investigate and potentially forecast postfire hydrologic hazards at regional scales.