Woody plant encroachment is a global phenomenon, observed in many of the world's drylands. In those with shallow soils overlying karst geology, rock moisture can be an important source of water for the encroaching woody plants. This source can be particularly important for trees to maintain basic physiological functions during extended droughts, which are becoming more frequent and intense owing to climate change. However, our understanding of rock moisture dynamics in karst drylands undergoing woody plant encroachment is still limited because of the scarcity of direct measurements. In this study, we evaluated soil and rock moisture dynamics at a semiarid site in the Edwards Plateau region of Texas. Our measurements over the course of three years showed that in shallow upslope terrain, the dynamic water storage in bedrock was roughly twice that of soil, while in downslope terrain, the dynamic storage was largely restricted to the soil layer. Most of the bedrock storage gains occurred during the first year, after two major storm events of approximately 95 mm, and that storage was gradually depleted during the following two years, when precipitation was below average. Importantly, in upslope terrain we found substantially larger water storage under woody plants, which suggests that they not only can access and utilize rock moisture but also play a role in enhancing bedrock water storage capacity. These interconnected abilities can help woody plants survive extended droughts---a factor crucial for understanding their persistence and proliferation in the shallow soils of the Edwards Plateau and similar karst regions.

Across diverse biomes and climate types, plants use water stored in bedrock to sustain transpiration. Bedrock water storage ($S_{bedrock}$, mm), in addition to soil moisture, thus plays an important role in water cycling and should be accounted for in the context of surface energy balances and streamflow generation. Yet, the extent to which bedrock water storage impacts hydrologic partitioning and influences latent heat fluxes has yet to be quantified at large scales. This is particularly important in Mediterranean climates, where the majority of precipitation is offset from energy delivery and plants must rely on water retained from the wet season to support summer growth. Here we present a simple water balance approach and random forest model to quantify the role of $S_{bedrock}$ on controlling hydrologic partitioning and land surface energy budgets. Specifically, we track evapotranspiration in excess of precipitation and mapped soil water storage capacity ($S_{soil}$, mm) across the western US in the context of Budyko’s water partitioning framework. Our findings indicate that $S_{bedrock}$ is necessary to sustain plant growth in forests in the Sierra Nevada — some of the most productive forests on Earth — as early as April every year, which is counter to the current conventional thought that bedrock is exclusively used late in the dry season under extremely dry conditions. We show that the average latent heat flux used in evapotranspiration of $S_{bedrock}$ can exceed 100 $W/m^{2}$ during the dry season and the proportion of water that returns to the atmosphere would decrease dramatically without access to $S_{bedrock}$.

Accurate observation of hillslope groundwater storage and instantaneous recharge remains difficult due to limited monitoring and the complexity of mountainous landscapes. We introduce a novel storage-discharge method to estimate hillslope recharge and the recharge ratio---the fraction of precipitation that recharges groundwater. The method, which relies on streamflow data, is corroborated by independent measurements of water storage dynamics inside the Rivendell experimental hillslope at the Eel River Critical Zone Observatory, California USA. We find that along-hillslope patterns in bedrock weathering and plant-driven storage dynamics govern the seasonal evolution of recharge ratios. Thinner weathering profiles and smaller root-zone storage deficits near-channel are replenished before larger ridge-top deficits. Consequently, precipitation progressively activates groundwater from channel to divide, with an attendant increase in recharge ratios throughout the wet season. Our novel approach and process observations offer valuable insights into controls on groundwater recharge, enhancing our understanding of a critical flux in the hydrologic cycle.

The water storage capacity of the root zone determines whether plants survive dry periods and controls the partitioning of precipitation into streamflow and evapotranspiration. It is currently thought that top-down, climatic factors are the primary control on this capacity via their interaction with plant rooting adaptations. However, it remains unclear to what extent bottom-up, geologic factors can provide an additional constraint on storage capacity. Here we use a machine learning approach to identify regions with lower than climatically expected apparent storage capacity. We find that in seasonally dry California these regions overlap with particular geologic substrates. We hypothesize that these patterns reflect diverse mechanisms by which substrate can limit storage capacity, and highlight case studies consistent with limited weathered bedrock extent (melange in the Northern Coast Range), toxicity (ultramafic substrates in the Klamath-Siskiyou region), nutrient limitation (phosphorus-poor plutons in the southern Sierra Nevada), and low porosity capable of retaining water (volcanic formations in the southern Cascades). The observation that at regional scales climate alone does not ‘size’ the root zone has implications for the parameterization of storage capacity in models of plant dynamics (and the interrelated carbon and water cycles), and also underscores the importance of geology in considerations of climate-change induced biome migration and habitat suitability.

In forested, seasonally dry watersheds, winter rains commonly replenish water storage deficits in the vadose zone before recharging underlying hillslope groundwater systems that sustain streamflow. However, the relative inaccessibility of the subsurface has hindered efforts to include the role of storage deficits, primarily generated by plant-water uptake, in moderating groundwater recharge. Here, we compare groundwater recharge inferred from the storage-discharge relationship with independent, distributed estimates of vadose zone storage deficits across 12 undisturbed California watersheds, thereby tracking the evolution of the deficit-recharge relationship without intensive field instrumentation. We find accrued deficits during the dry season (driven by evapotranspiration) alone insufficiently explain differences in the wet season partitioning of rainfall due to the inter-storm behavior of the deficit during the subsequent wet season. Tracking the deficit at the storm event-scale within the wet season reveals a characteristic response in groundwater to increasing rainfall not captured in the seasonal analysis, and may improve estimates of the rainfall required to generate recharge and streamflow on a per-storm basis. Our findings demonstrate the potential for existing public datasets to better capture water partitioning within the subsurface using a combined deficit-recharge approach, though our analysis is currently limited to basins with select characteristics.

Understanding how soil thickness and bedrock weathering vary across ridge and valley topography is needed to constrain the flowpaths of water and sediment production within a landscape. Here, we investigate saprolite and weathered bedrock properties across a ridge-valley system in the Northern California Coast Ranges, USA, where topography varies with slope aspect such that north facing slopes have thicker soils and are more densely vegetated than south facing slopes. We use active source seismic refraction surveys to extend observations made in boreholes to the hillslope scale. Seismic velocity models across several ridges capture a high velocity gradient zone (from 1000 to 2500 m/s) located ~4-13 m below ridgetops, that coincides with transitions in material strength and chemical depletion observed in boreholes. Comparing this transition depth across multiple north and south-facing slopes, we find that the thickness of saprolite does not vary with slope aspects. Additionally, seismic survey lines perpendicular and parallel to bedding planes reveal weathering profiles that thicken upslope and taper downslope to channels. Using a rock physics model incorporating seismic velocity, we estimate the total porosity of the saprolite and find that inherited fractures contribute a substantial amount of pore space in the upper 6 m, and the lateral porosity structure varies strongly with hillslope position. The aspect-independent weathering structure suggests the contemporary critical zone structure at Rancho Venada is a legacy of past climate and vegetation conditions.

Quantifying the volume of water that is stored in the subsurface is critical to studies of water availability to ecosystems, slope stability, and water-rock interactions. In a variety of settings, water is stored in fractured and weathered bedrock as rock moisture. However, few techniques are available to measure rock moisture in unsaturated rock, making direct estimates of water storage dynamics difficult to obtain. Here, we use borehole nuclear magnetic resonance (NMR) at two sites in seasonally dry California to quantify dynamic rock moisture storage. We show strong agreement between NMR estimates of dynamic storage and estimates derived from neutron logging and mass balance techniques. The depths of dynamic storage are up to 9 m and likely reflect the depth extent of root water uptake. To our knowledge, these data are the first to quantify the volume and depths of dynamic water storage in the bedrock vadose zone via NMR.

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.

A better understanding of how vegetation influences alluvial channels could improve (a) assessments of channel stability and flood risks, (b) applications of vegetation as a river management tool, and (c) predictions of channel responses to climate change and other human impacts. We take advantage of a natural field experiment in the semi-arid to arid Henry Mountains, Utah, USA: Large spatial differences in bed and bank vegetation are found along some alluvial channels due to localized perennial springs caused by aquicludes in the underlying bedrock. Airborne LiDAR topography and flood modeling are used to constrain channel morphology, vegetation density, and flow velocity at different flood discharges for three spring-fed reaches along intermittently-flowing streams. The spatial distribution of vegetation quantitatively influences both the magnitude and direction of channel adjustment. Reaches with abundant bed vegetation are significantly wider (by an average of ≈ 50%), with shallower flows and lower velocities, than reaches with little bed vegetation. Reaches with dense channel bank vegetation are ≈ 25% narrower and ≈ 25% deeper than sparsely-vegetated reaches. We interpret that sediment grain size influences the spatial distribution of vegetation within spring reaches, but that bank vegetation may be more important than grain size for “threshold” width adjustments. Widths, depths and velocities are fairly insensitive to whether local hydraulic roughness is parameterized in terms of local vegetation density or is assumed spatially constant, suggesting that the underlying “bare earth” topography of the channel bed, banks and floodplain exerts more control on local flow than does local vegetation density.

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.