Abstract Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first‐order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated‐soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage. 
                        more » 
                        « less   
                    
                            
                            Does soil erosion rejuvenate the soil phosphorus inventory?
                        
                    
    
            Phosphorus (P) is an essential nutrient for life. Deficits in soil P reduce primary production and alter biodiversity. A soil P paradigm based on studies of soils that form on flat topography, where erosion rates are minimal, indicates P is supplied to soil mainly as apatite from the underlying parent material and over time is lost via weathering or transformed into labile and less-bioavailable secondary forms. However, little is systematically known about P transformation and bioavailability on eroding hillslopes, which make up the majority of Earth's surface. By linking soil residence time to P fractions in soils and parent material, we show that the traditional concept of P transformation as a function of time has limited applicability to hillslope soils of the western Southern Alps (New Zealand) and Northern Sierra Nevada (USA). Instead, the P inventory of eroding soils at these sites is dominated by secondary P forms across a range of soil residence times, an observation consistent with previously published soil P data. The findings for hillslope soils contrast with those from minimally eroding soils used in chronosequence studies, where the soil P paradigm originated, because chronosequences are often located on landforms where parent materials are less chemically altered and therefore richer in apatite P compared to soils on hillslopes, which are generally underlain by pre-weathered parent material (e.g., saprolite). The geomorphic history of the soil parent material is the likely cause of soil P inventory differences for eroding hillslope soils versus geomorphically stable chronosequence soils. Additionally, plants and dust seem to play an important role in vertically redistributing P in hillslope soils. Given the dominance of secondary soil P in hillslope soils, limits to ecosystem development caused by an undersupply of bio-available P may be more relevant to hillslopes than previously thought. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 1015454
- PAR ID:
- 10635926
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Geoderma
- Volume:
- 332
- Issue:
- C
- ISSN:
- 0016-7061
- Page Range / eLocation ID:
- 45 to 59
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Chemical weathering influences many aspects of the Earth system, including biogeochemical cycling, climate, and ecosystem function. Physical erosion influences chemical weathering rates by setting the supply of fresh minerals to the critical zone. Vegetation also influences chemical weathering rates, both by physical processes that expose mineral surfaces and via production of acids that contribute to mineral dissolution. However, the role of vegetation in setting surface process rates in different landscapes is unclear. Here we use 10Be and geochemical mass balance to quantify soil production, physical erosion, and chemical weathering rates in a landscape where a migrating drainage divide separates catchments with an order-of magnitude contrast in erosion rates and where vegetation spans temperate rainforest, tussock grassland, and unvegetated alpine ecosystems in the western Southern Alps of New Zealand. Soil production, physical erosion, and chemical weathering rates are significantly higher on the rapidly eroding versus the slowly eroding side of the drainage divide. However, chemical weathering intensity does not vary significantly across the divide or as a function of vegetation type. Soil production rates are correlated with ridgetop curvature, and ridgetops are more convex on the rapidly eroding side of the divide, where soil mineral residence times are lowest. Hence our findings suggest fluvially-driven erosion rates control soil production and soil chemical weathering rates by influencing the relationship between hillslope topography and mineral residence times. In the western Southern Alps, soil production and chemical weathering rates are more strongly mediated by physical rock breakdown driven by landscape response to tectonics, than by vegetation.more » « less
- 
            Abstract Process‐based modeling of soil water movement with the Richards equation requires the description of soil hydraulic material properties, which are highly uncertain and heterogeneous at all scales. This limits the applicability of the Richards equation at larger scales beyond the patch scale. The experimental capabilities of the three hillslopes of the Landscape Evolution Observatory (LEO) at Biosphere 2 provide a unique opportunity to observe the heterogeneity of hydraulic material properties at the hillslope scale. We performed a gravity flow experiment where through constant irrigation the water content increases until the hydraulic conductivity matches the irrigation flux above. The dense water content sensor network at LEO then allows mapping of the heterogeneity of hydraulic conductivity at a meter scale resolution. The experiment revealed spatial structures within the hillslopes, mainly a vertical trend with the lowest hydraulic conductivity close to the surface. However, the variation between neighboring sensors is high, showing that the heterogeneity cannot be fully resolved even at LEO. By representing the heterogeneity in models through Miller scaling we showed the impact on hillslope discharge. For the hillslope with the smallest heterogeneity, representing the dominant structures was sufficient. However, for the two hillslopes with the larger overall heterogeneity, adding further details of the local heterogeneity did impact the discharge further. This highlights the limitations of the Richards equation, which requires the heterogeneous field of material properties, at the hillslope scale and shows the relevance to improving our understanding of effective parameters to be able to apply the process‐based model to larger scales.more » « less
- 
            Abstract Debris flows are powered by sediment supplied from steep hillslopes where soils are often patchy and interrupted by bare‐bedrock cliffs. The role of patchy soils and cliffs in supplying sediment to channels remains unclear, particularly surrounding wildfire disturbances that heighten debris‐flow hazards by increasing sediment supply to channels. Here, we examine how variation in soil cover on hillslopes affects sediment sizes in channels surrounding the 2020 El Dorado wildfire, which burned debris‐flow prone slopes in the San Bernardino Mountains, California. We focus on six headwater catchments (<0.1 km2) where hillslope sources ranged from a continuous soil mantle to 95% bare‐bedrock cliffs. At each site, we measured sediment grain size distributions at the same channel locations before and immediately following the wildfire. We compared results to a mixing model that accounts for three distinct hillslope sediment sources distinguished by local slope thresholds. We find that channel sediment in fully soil‐mantled catchments reflects hillslope soils (D50 = 0.1–0.2 cm) both before and after the wildfire. In steeper catchments with cliffs, channel sediment is consistently coarse prior to fire (D50 = 6–32 cm) and reflects bedrock fracture spacing, despite cliffs representing anywhere from 5% to 95% of the sediment source area. Following the fire, channel sediment size reduces most (5‐ to 20‐fold) in catchments where hillslope sources are predominantly soil covered but with patches of cliffs. The abrupt fining of channel sediment is thought to facilitate postfire debris‐flow initiation, and our results imply that this effect is greatest where bare‐bedrock cliffs are present but not dominant. A patchwork of bare‐bedrock cliffs is common in steeplands where hillslopes respond to channel incision by landsliding. We show how local slope thresholds applied to such terrain aid in estimating sediment supply conditions before two destructive debris flows that eventually nucleated in these study catchments in 2022.more » « less
- 
            Hydrologic behavior and soil properties across forested landscapes with complex topography exhibit high variability. The interaction of groundwater with spatially distinct soils produces and transports solutes across catchments, however, the spatiotemporal relationships between groundwater dynamics and soil solute fluxes are difficult to directly evaluate. While whole-catchment export of solutes by shallow subsurface flow represents an integration of soil environments and conditions but many studies compartmentalize soil solute fluxes as hillslope vs. riparian, deep vs. shallow, or as individual soil horizon contributions. This potentially obscures and underestimates the hillslope variation and magnitude of solute fluxes and soil development across the landscape. This study determined the spatial variation and of shallow soil base cation fluxes associated with weathering reactions (Ca, Mg, and Na), soil elemental depletion, and soil saturation dynamics in upland soils within a small, forested watershed at the Hubbard Brook Experimental Forest, NH. Base cation fluxes were calculated using a combination of ion-exchange resins placed in shallow groundwater wells (0.3 – 1 m depth) located across hillslope transects (ridges to lower backslopes) and measurements of groundwater levels. Groundwater levels were also used to create metrics of annual soil saturation. Base cation fluxes were positively correlated with soil saturation frequency and were greatest in soil profiles where primary minerals were most depleted of base cations (i.e., highly weathered). Spatial differences in soil saturation across the catchment were strongly related to topographic properties of the upslope drainage area and are interpreted to result from spatial variations in transient groundwater dynamics. Results from this work suggest that the structure of a catchment defines the spatial architecture of base cation fluxes, likely reflecting the mediation of subsurface stormflow dynamics on soil development. Furthermore, this work highlights the importance of further compartmentalizing solute fluxes along hillslopes, where certain areas may disproportionately contribute solutes to the whole catchment. Refining catchment controls on base cation generation and transport could be an important tool for opening the black box of catchment elemental cycling.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    