Understanding and predicting catchment responses to a regional disturbance is difficult because catchments are spatially heterogeneous systems that exhibit unique moderating characteristics. Changes in precipitation composition in the Northeastern U.S. is one prominent example, where reduction in wet and dry deposition is hypothesized to have caused increased dissolved organic carbon (DOC) export from many northern hemisphere forested catchments; however, findings from different locations contradict each other. Using shifts in acid deposition as a test case, we illustrate an iterative “process and pattern” approach to investigate the role of catchment characteristics in modulating the steam DOC response. We use a novel dataset that integrates regional and catchment-scale atmospheric deposition data, catchment characteristics and co-located stream Q and stream chemistry data. We use these data to investigate opportunities and limitations of a pattern-to-process approach where we explore regional patterns of reduced acid deposition, catchment characteristics and stream DOC response and specific soil processes at select locations. For pattern investigation, we quantify long-term trends of flow-adjusted DOC concentrations in stream water, along with wet deposition trends in sulfate, for USGS headwater catchments using Seasonal Kendall tests and then compare trend results to catchment attributes. Our investigation of climatic, topographic, and hydrologic catchment attributes vs. directionality of DOC trends suggests soil depth and catchment connectivity as possible modulating factors for DOC concentrations. This informed our process-to-pattern investigation, in which we experimentally simulated increased and decreased acid deposition on soil cores from catchments of contrasting long-term DOC response [Sleepers River Research Watershed (SRRW) for long-term increases in DOC and the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO) for long-term decreases in DOC]. SRRW soils generally released more DOC than SSHCZO soils and losses into recovery solutions were higher. Scanning electron microscope imaging indicates a significant DOC contribution from destabilizing soil aggregates mostly from hydrologically disconnected landscape positions. Results from this work illustrate the value of an iterative process and pattern approach to understand catchment-scale response to regional disturbance and suggest opportunities for further investigations.
Despite the absence of tectonic activity, cratonic environments are characterized by strongly variable, and in places significant, rock weathering rates. This is shown here through an exploration of the weathering rates in two inter-tropical river basins from the Atlantic Central Africa: the Ogooué and Mbei River basins, Gabon. We analyzed the elemental and strontium isotope composition of 24 water samples collected throughout these basins. Based on the determination of the major element sources we estimate that the Ogooué and Mbei rivers total dissolved solids (TDS) mainly derive from silicate chemical weathering. The chemical composition of the dissolved load and the area-normalized solute fluxes at the outlet of the Ogooué are similar to those of other West African rivers (e.g., Niger, Nyong, or Congo). However, chemical weathering rates ( rate expressed as the release rate of the sum of cations by silicate chemical weathering) span the entire range of chemical weathering intensities hitherto recorded in worldwide cratonic environments. In the Ogooué-Mbei systems, three regions can be distinguished: (i) the Eastern sub-basins draining the Plateaux Batéké underlain by quartz-rich sandstones exhibit the lowest rates, (ii) the Northern sub-basins and the Mbei sub-basins, which drain the southern edge of the tectonically quiescent South Cameroon Plateau, show intermediate rates and (iii) the Southern sub-basins characterized by steeper slopes record the highest rates. In region (ii), higher DOC concentrations are associated with enrichment of elements expected to form insoluble hydrolysates in natural waters (e.g., Fe, Al, Th, REEs) suggesting enhanced transport of these elements in the colloidal phase. In region (iii), we suggest that a combination of mantle-induced dynamic uplift and lithospheric destabilization affecting the rim of the Congo Cuvette induces slow base level lowering thereby enhancing soil erosion, exhumation of fresh primary minerals, and thus weathering rates. The study points out that erosion of lateritic covers in cratonic areas can significantly enhance chemical weathering rates by bringing fresh minerals in contact with meteoric water. The heterogeneity of weathering rates amongst cratonic regions thus need to be considered for reconstructing the global, long-term carbon cycle and its control on Earth climate.
Hyporheic zones act as critical ecological links between terrestrial and aquatic systems where redox-sensitive metals of iron (Fe) and manganese (Mn) significantly impact nutrient cycling and water quality. However, the geochemical controls on the release and speciation of Fe(II) and Mn(II) in these biogeochemical hotspots are still poorly understood. Here we conducted batch incubation experiments and analyzed Fe K-edge extended X-ray absorption fine structure (EXAFS) spectroscopy data using sediment samples from a hyporheic zone of the East River floodplain in Colorado to understand the production, release and speciation of Fe(II) and Mn(II) in groundwater. Our results indicate that the production and release of Fe(II) and Mn(II) vary with sediment reducing conditions and subsurface positions, and the rates were determined either by a zero- or first-order rate equation. The sediments with higher Fe(II) production did not necessarily result in higher release of dissolved Fe(II), and ≥97% Fe(II) is accumulated in solid phase. We found that the majority of Fe(II) exists as siderite (FeCO3), Fe(II)-natural organic matter (NOM) complexes and ferrosmectite, and the equilibrium concentrations of dissolved Fe(II) are controlled primarily by siderite solubility, and enhanced greatly by formation of strong Fe(II)-NOM complexes as dominant aqueous Fe(II) species. By contract, dissolved Mn(II) increases slowly and linearly, and an equilibrium concentration was not reached during the incubation period, and the roles of rhodochrosite (MnCO3) and Mn(II)-NOM complexes are insignificant. Furthermore, we reviewed and calibrated the literature reported binding constants (log K) of Fe(II)-NOM complexes which successfully predicted our experimental data. This work reveals that siderite and dissolved NOM are the controlling phases in release and speciation of dissolved Fe(II), and the finding is expected to be applicable in many hyporheic zones and subsurface environments with similar geochemical conditions.
Climate warming in alpine regions is changing patterns of water storage, a primary control on alpine plant ecology, biogeochemistry, and water supplies to lower elevations. There is an outstanding need to determine how the interacting drivers of precipitation and the critical zone (CZ) dictate the spatial pattern and time evolution of soil water storage. In this study, we developed an analytical framework that combines intensive hydrologic measurements and extensive remotely-sensed observations with statistical modeling to identify areas with similar temporal trends in soil water storage within, and predict their relationships across, a 0.26 km2 alpine catchment in the Colorado Rocky Mountains, U.S.A. Repeat measurements of soil moisture were used to drive an unsupervised clustering algorithm, which identified six unique groups of locations ranging from predominantly dry to persistently very wet within the catchment. We then explored relationships between these hydrologic groups and multiple CZ-related indices, including snow depth, plant productivity, macro- (102->103 m) and microtopography (<100-102 m), and hydrological flow paths. Finally, we used a supervised machine learning random forest algorithm to map each of the six hydrologic groups across the catchment based on distributed CZ properties and evaluated their aggregate relationships at the catchment scale. Our analysis indicated that ~40–50% of the catchment is hydrologically connected to the stream channel, lending insight into the portions of the catchment that likely dominate stream water and solute fluxes. This research expands our understanding of patch-to-catchment-scale physical controls on hydrologic and biogeochemical processes, as well as their relationships across space and time, which will inform predictive models aimed at determining future changes to alpine ecosystems.
Concentration-discharge (C-Q) relations can provide insight into the dynamic behavior of the Critical Zone (CZ), as C-Q relations integrate the spatial distribution and timing of watershed hydrogeochemical processes. This study blends geomorphologic analysis, C-Q relations and reactive-transport modeling using a rich dataset from an elevation gradient of eight watersheds in the Southern Sierra Nevada, California. We found that the CZ structure exerts a strong control on the C-Q relations, and on the hydrogeochemical behavior of headwater watersheds. Watersheds with thin regolith, a large stream network, and limited water storage have fast mean transit times along subsurface flow lines, and show limited seasonal variability in ionic concentrations in streamflow (i.e., chemostatic behavior). In contrast, watersheds with thicker regolith, a small stream network and more water storage have longer transit times along subsurface flow lines, and exhibit greater chemical variability (i.e., chemodynamic behavior). Independent estimates of mean transit times and water storage from other isotopic, hydrologic and geophysical studies were consistent with results from modeling C-Q relations. The stream chemistry and its variability were controlled by lateral flow within the regolith, and no mixing with deep groundwater was needed to explain the observed chemical variability. This study opens the possibility to estimate water-storage capacity and mean transit times, and thus drought resistance in watersheds, by using quantitative modeling of C-Q relations.