Bedrock vadose zone water storage (i.e., rock moisture) dynamics are sparsely observed but potentially key to understanding drought responses. Exploiting a borehole network at a Mediterranean blue oak savanna site-Rancho Venada-we document how water storage capacity in a deeply weathered bedrock profile regulates woody plant water availability and groundwater recharge. The site is in the Northern California Coast Range within steeply dipping turbidites. In a wet year (water year 2019; 647 mm of precipitation), rock moisture was quickly replenished to a characteristic storage capacity, recharging groundwater that emerged at springs to generate streamflow. In the subsequent rainless summer growing season, rock moisture was depleted by about 93 mm. In two drought years that followed (212 and 121 mm of precipitation) the total amount of rock moisture gained each winter was about 54 and 20 mm, respectively, and declines were observed exceeding these amounts, resulting in progressively lower rock moisture content. Oaks, which are rooted into bedrock, demonstrated signs of water stress in drought, including reduced transpiration rates and extremely low water potentials. In the 2020-2021 drought, precipitation did not exceed storage capacity, resulting in variable belowground water storage, increased plant water stress, and no recharge or runoff. Rock moisture deficits (rather than soil moisture deficits) explain these responses.

Unlike streamflow, which can be sampled in aggregate at the catchment outlet, evapotranspiration (ET) is spatially dispersed, challenging large-scale age estimation. Here, we introduce an approach for constraining the age of ET via mass balance and present the minimum flux-weighted age of ET across the continental US using distributed, publicly available water flux datasets. The lower-bound constraint on ET age can be calculated by assuming that ET is preferentially sourced from the most recent precipitation through a last-in, first-out algorithm. From 2012-2017, ET was at least several months old across large areas of the western continental US, including in Mediterranean and (semi-)arid climate zones and shrub and evergreen needleleaf plant communities. The primary limitation of this approach is that it provides only a minimum flux-weighted average age to satisfy mass balance of outgoing fluxes; true ET fluxes are composed of distributions of ages and may be composed of much older water. The primary advantage of the approach is that flux timeseries of precipitation and ET are sufficient to constrain ET age, and model parameterization is unnecessary. ET ages can be used to validate tracer-aided and modeling approaches and inform studies of biogeochemistry, water-rock interactions, and plant water sourcing under drought.

Understanding how soil thickness and bedrock weathering vary across ridge and valley topography is needed to constrain the flowpaths of water and sediment within a landscape. Here, we investigate how soil and weathered bedrock properties vary across a ridge-valley system in the Northern California Coast Ranges where topography varies with slope aspect such that north facing slopes, which are more densely vegetated, are steeper. In this study, we use seismic refraction surveys to extend observations made in boreholes and soil pits to the hillslope scale and identify that while soils are thicker on north facing slopes, the thickness of weathered bedrock does not vary with slope aspect. We estimate the porosity of the weathered bedrock and find that it is several times the annual rainfall, indicating that water storage is not limited by the available pore space, but rather the amount of precipitation delivered. Bedding-parallel and bedding-perpendicular seismic refraction surveys reveal weathering profiles that are thickest upslope and taper downslope to channels. We do not find a clear linear scaling relationship between depth to bedrock and hillslope length, which may be due to local variation in incision rate or bedrock hydraulic conductivity. Together, these findings, which suggest that the aspect-independent weathering profile structure is a legacy of past climate and vegetation conditions and that weathering varies strongly with hillslope position, have implications for hydrologic processes across this landscape.

Water age and flow pathways should be related; however, it is still generally unclear how integrated catchment runoff generation mechanisms result in streamflow age distributions at the outlet. Here, we combine field observations of runoff generation at the Dry Creek catchment with StorAge Selection (SAS) age models to explore the relationship between streamwater age and runoff pathways. Dry Creek is a 3.5 km2 catchment in the Northern California Coast Ranges with a Mediterranean climate, and, despite an average rainfall of ~1,800 mm/yr, is an oak savannah due to the limited water storage capacity. Runoff lag to peak—after initial seasonal wet-up—is rapid (~1-2 hours), and total annual streamflow consists predominantly of saturation overland flow, based on field mapping of saturated extents and an inferred runoff threshold for the expansion of saturation extent beyond the geomorphic channel. SAS modeling based on daily isotope sampling reveals that streamflow is typically older than one day. Because streamflow is mostly overland flow, this means that a significant portion of overland flow must not be event-rain but instead derive from older, non-event groundwater returning to the surface, consistent with field observations of exfiltrating head gradients, return flow through macropores, and extensive saturation days after storm events. We conclude that even in a landscape with widespread overland flow, runoff pathways may be longer and slower than anticipated. Our findings have implications for the assumptions built into widely used hydrograph separation inferences, namely, the assumption that overland flow consists of new (event) water.

Predicting the behavior of overland flow with analytical solutions to the kinematic wave equation is appealing due to its relative ease of implementation. Such simple solutions, however, have largely been constrained to applications on simple planar hillslopes. This study presents analytical solutions to the kinematic wave equation for hillslopes with modest topographic curvature that causes divergence or convergence of runoff flowpaths. The solution averages flow depths along changing hillslope contours whose lengths vary according hillslope width function, and results in a one-dimensional approximation to the two-dimensional flow field. The solutions are tested against both two-dimensional numerical solutions to the kinematic wave equation (in ParFlow) and against experiments that use rainfall simulation on machined hillslopes with defined curvature properties. Excellent agreement between numerical, experimental and analytical solutions is found in all cases. The solutions show that curvature drives large changes in maximum flow rate qmax and time of concentration tc, predictions frequently used in engineering hydrologic design and analysis.