Search for: All records

Abstract This study explores the impact of deep (5–40 m) critical zone (CZ) structure on vegetation distribution in a semiarid snow‐dominated climate. Utilizing seismic refraction surveys, we identified a significant negative correlation between seismically derived saprolite thickness and light detecting and ranging‐derived vegetation heights (R= −0.66). We argue that CZ structure, specifically shallow fractured bedrock under valley bottoms, provides moisture near the surface where trees are established—suggesting the trees are situated in locations with access to nutrients and water. This work provides a unique spatially exhaustive perspective and adds to growing evidence that in addition to other factors such as slope, aspect, and climate, deep CZ structure plays a vital role in ecosystem development. 

Abstract Within Earth's critical zone, weathering processes influence landscape evolution and hillslope hydrology by creating porosity in bedrock, transforming it into saprolite and eventually soil. In situ weathering processes drive much of this transformation while preserving the rock fabric of the parent material. Inherited rock fabric in regolith makes the critical zone anisotropic, affecting its mechanical and hydrological properties. Therefore, quantifying and studying anisotropy is an important part of characterising the critical zone, yet doing so remains challenging. Seismic methods can be used to detect rock fabric and infer mechanical and hydrologic conductivity anisotropy across landscapes. We present a novel way of measuring seismic anisotropy in the critical zone using Rayleigh and Love surface waves. This method leverages multi‐component surface seismic data to create a high‐resolution model of seismic anisotropy, which we compare with a nuclear magnetic resonance log measured in a nearby borehole. The two geophysical data sets show that seismic anisotropy and porosity develop at similar depths in weathered bedrock and both reach their maximum values in saprolite, implying that in situ weathering enhances anisotropy while concurrently generating porosity in the critical zone. We bolster our findings with in situ measurements of seismic and hydrologic conductivity anisotropy made in a 3 m deep soil excavation. Our study offers a fresh perspective on the importance of rock fabric in the development and function of the critical zone and sheds new insights into how weathering processes operate. 

ABSTRACT To accurately predict earth system response to global change, we must be able to predict the responses of important properties of that system, such as the depths over which plant roots are distributed. In 2008, H. J. Schenk proposed a model for the depth distribution of plant roots based on a simple hydrological scheme and the assumptions that plants will take up the shallowest water available first and will distribute their roots in proportion to long‐term mean uptake at each depth. Here, we derive an analytical solution to the Schenk model under an idealised climate (in which infiltration events are treated as a marked Poisson process), explore properties of the result and compare with data. The solution suggests that in very humid and arid climates, the soil wetting and drying cycles induced by root water uptake are generally confined to a characteristic depth below the surface. This depth depends on the typical magnitude of rainfall events (most strongly so in arid climates), the typical total transpiration demand between rainfall events (most strongly in humid climates) and the plant‐available water holding capacity of the soil. Root water uptake (and thus predicted root density) in very humid and arid landscapes decreases exponentially with depth at a rate determined by this characteristic depth. However, in a mesic climate, soils may be wet or dry to greater depths below the near‐surface, and the duration spent in each state increases with depth. Consequently, root water uptake and root density in mesic climates more closely resemble a power law distribution. When the aridity index is exactly 1, the characteristic depth diverges and the mean rooting depth approaches infinity. This suggests that the most skewed root depth distributions might occur in mesic environments. We compared this model to another analytical solution and a compiled database of root distributions (159 combined locations). For a larger comparison dataset, we also compared 99th percentile rooting depth to rooting depths modeled by two other frameworks and a database of observed rooting depths (1271 combined locations). Results demonstrate that the analytical formulation of the Schenk model performs well as a shallow bound on rooting depths and captures something of the nonexponential form of root distributions, and its error is similar to or less than that of other modeling frameworks. Errors may be partly explained by the deviation of real climate from the idealisations used to obtain an analytical solution (exponentially distributed infiltration events and no seasonality). 

Abstract Topography is a key control on runoff generation, as topographic slope affects hydraulic gradients and curvature affects water flow paths. Simultaneously, runoff generation shapes topography through erosion, affecting landscape morphology over long timescales. Previous modeling efforts suggest that subsurface hydrological properties, relative to climate, are key mediators of this relationship. Specifically, when subsurface transmissivity and water storage capacity are low, (a) saturated areas and storm runoff should be larger and more variable, and (b) hillslopes shorter and lower relief, assuming other geomorphic factors are held constant. However, it remains uncertain whether subsurface properties can exert such strong controls on emergent properties in real landscapes. We compared emergent hydrological function and topography in two watersheds with very similar climatic and tectonic history, but very different subsurface properties due to contrasting bedrock lithology. We found that hillslopes were systematically shorter and saturated areas more dynamic at the lower transmissivity site. To test whether these features could be the result of coevolution between topography, hydrological function, and subsurface properties, we estimated all parameters of a coupled groundwater‐landscape evolution model for each site. Limitations were revealed in the model's ability to reproduce aspects of morphology and hydrologic behavior, however, model results suggested differences in hillslope length and variably saturated area between the sites could be explained by differences in subsurface properties, and not by differences in geomorphic process rates alone. This work demonstrates one way subsurface hydrology can profoundly affect landscape evolution. 

Abstract Controls on the physical and chemical architecture of the subsurface critical zone are somewhat controversial, with multiple hypotheses proposed to account for variations in the depth of weathering between sites, and with landscape position at a site. In the Piedmont region of the Mid‐Atlantic US weathering of crystalline bedrock has been observed to extend tens of meters below the surface and groundwater in a'bow‐tie’ shape – i.e. weathering extends to lower elevations below ridges than below channels. The chemical and physical structure of a hillslope transect in the Maryland Piedmont was explored with a 45 m borehole in the ridge, as well as shallow bedrock boreholes at the toe of the slope and valley. Chemical weathering fronts were characterized using elemental abundances and mineralogical analysis. The ridge borehole did not extend deeper than the chemically and physically weathered rock. Surface and borehole geophysics and density measurements were used to characterize the weathered rock and saprolite. Na and Ca results suggest that plagioclase feldspar weathering is similar between samples collected from 45 m under the ridge and 2.2 m under the valley bottom. A narrow Fe oxidation garnet weathering front co‐insides with the transition from weathered bedrock to saprolite, suggesting that this reaction may generate initial saprolite porosity. Muscovite weathering co‐occurs with complete depletion of plagioclase a few meters above the Fe oxidation front. These nested weathering fronts in the saprolite appear to follow a subdued version of the surface topography. The location and shape of the nested saprolite weathering fronts may be controlled by the feedback between the transport of reactants and solutes and reaction‐generated porosity, consistent with the conceptual “valve” hypothesis. Differing dominant control mechanisms on deep bedrock weathering and saprolite initiating reactions may explain the thickness and structure of the critical zone at our site. 

Abstract Features of landscape morphology—including slope, curvature, and drainage dissection—are important controls on runoff generation in upland landscapes. Over long timescales, runoff plays an essential role in shaping these same features through surface erosion. This feedback between erosion and runoff generation suggests that modeling long‐term landscape evolution together with dynamic runoff generation could provide insight into hydrological function. Here we examine the emergence of variable source area runoff generation in a new coupled hydro‐geomorphic model that accounts for water balance partitioning between surface flow, subsurface flow, and evapotranspiration as landscapes evolve over millions of years. We derive a minimal set of dimensionless numbers that provide insight into how hydrologic and geomorphic parameters together affect landscapes. Across the parameter space we investigated, model results collapsed to a single inverse relationship between the dimensionless relief and the ratio of catchment quickflow to discharge. Furthermore, we found an inverse relationship between the Hillslope number, which describes topographic relief relative to aquifer thickness, and the proportion of the landscape that was variably saturated. While the model generally produces fluvial topography visually similar to simpler landscape evolution models, certain parameter combinations produce wide valley bottom wetlands and non‐dendritic, trellis‐like drainage networks, which may reflect real conditions in some landscapes where aquifer gradients become decoupled from topography. With these results, we demonstrate the power of hydro‐geomorphic models for generating new insights into hydrological processes, and also suggest that subsurface hydrology may be integral for modeling aspects of long‐term landscape evolution. 

Abstract Sulfate is a potential pollutant and important nutrient linked with the nitrogen, carbon, and phosphorus cycles. The importance of different anthropogenic sulfate sources in suburban streams (septic systems, fertilizer, road salt, and infrastructure) is uncertain, and the temporal dynamics of stream export sparsely documented. We study sources and export dynamics of sulfate in suburban and forested headwater catchments. Stream baseflow discharge and sulfate concentrations were strongly positively correlated in both watersheds with the highest values in spring. Suburban concentrations and fluxes (2.48–7.5 mg/L or 25.8–78.1 μM, 16.6 kg/ha/yr) were consistently higher than forested (0.56–2.78 mg/L or 5.8–28.9 μM, 5 kg/ha/yr). Following precipitation, sulfate concentrations in both forested and suburban streams increased to concentrations above pre‐storm values and remained high after peak discharge. These dynamics suggest that both catchments have a large pool of sulfate that can be mobilized under wet conditions. Ridge‐top forest soil samples contained 210 kg/ha stored, extractable sulfate. Current atmospheric sulfate deposition rates (5–7 kg/ha/yr) are approximately in balance with sulfate export in the forested stream. In the suburban watershed, we estimated septic fields contribute up to 11 kg/ha/yr (about half from surfactants) and lawn care up to 4.3 kg/ha/yr and are the most likely sources of elevated stream sulfate. Sulfate sulfur (4.9–5.8‰ forested; 6.1–7.0‰ suburban) and oxygen isotope values (0.7–2.0‰ forested; −0.1–4.1‰ suburban) are consistent with this interpretation, but do not provide strong corroboration due to large variation and overlap in estimated source values. 

Data, model output, and scripts supporting the manuscript: Litwin, D. G., & Harman, C. J. (2024) Evidence of subsurface control on the coevolution of hillslope morphology and runoff generation. Water Resources Research, 60, e2024WR037301. https://doi.org/10.1029/2024WR037301