skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The effect of weathering in runoff-to-groundwater partitioning in the Island of Hawai'i: Perspectives for landscape evolution
We investigated how runoff-to-groundwater partitioning changes as a function of substrate age and degree of regolith development in the Island of Hawai’i, by modeling watershed-scale hydrodynamic properties for a series of volcanic catchments of different substrate age developed under different climates. In the younger catchments, rainfall infiltrates directly into the groundwater system and surface runoff is minimal, consisting of ephemeral streams flowing on the scale of hours to days. The older catchments show increasing surface runoff, with deeper incision and perennial discharge. We hypothesize that watershed-scale hydrodynamic properties change as a function of their weathering history—the convolution of time and climate: as surfaces age and become increasingly weathered, hydraulic conductivity is reduced, leading to increased runoff-to-recharge ratios. To test this relationship, we calculated both saturated hydraulic conductivity (k) and aquifer thickness (D) using recession flow analysis. We show that the average k in the younger catchments can be between 3 to 6 orders of magnitude larger than in older catchments, whereas modeled D increases with age. Ephemeral streams with zero baseflow at daily timescales cannot be evaluated using the same method. Instead, we calculated the recession constant for two contiguous catchments developed on young ash or lava deposits of different ages. Increasing bedrock age results in slower recession response in these ephemeral streams, which is consistent with decreasing hydraulic conductivity. Our results highlight the role of the weathering history in determining the evolution of watershed-scale hydrologic properties in volcanic catchments.  more » « less
Award ID(s):
2141989
PAR ID:
10563497
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Earth and Planetary Science Letters
Volume:
635
Issue:
C
ISSN:
0012-821X
Page Range / eLocation ID:
118687
Subject(s) / Keyword(s):
hydrology weathering runoff incision
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The western U.S. is experiencing increasing rain to snow ratios due to climate change, and scientists are uncertain how changing recharge patterns will affect future groundwater‐surface water connection. We examined how watershed topography and streambed hydraulic conductivity impact groundwater age and stream discharge at eight sites along a headwater stream within the Manitou Experimental Forest, CO USA. To do so, we measured: (a) continuous stream and groundwater discharge/level and specific conductivity from April to November 2021; (b) biweekly stream and groundwater chemistry; (c) groundwater chlorofluorocarbons and tritium in spring and fall; (d) streambed hydraulic conductivity; and (e) local slope. We used the chemistry data to calculate fluorite saturation states that were used to inform end‐member mixing analysis of streamflow source. We then combined chlorofluorocarbon and tritium data to estimate the age composition of riparian groundwater. Our data suggest that future stream drying is more probable where local slope is steep and streambed hydraulic conductivity is high. In these areas, groundwater source shifted seasonally, as indicated by age increases, and we observed a high fraction of groundwater in streamflow, primarily interflow from adjacent hillslopes. In contrast, where local slope is flat and streambed hydraulic conductivity is low, streamflow is more likely to persist as groundwater age was seasonally constant and buffered by storage in alluvial sediments. Groundwater age and streamflow paired with characterization of watershed topography and subsurface characteristics enabled identification of likely controls on future stream drying patterns. 
    more » « less
  2. Abstract Groundwater flow direction within the critical zone of headwater catchments is often assumed to mimic land surface topographic gradients. However, groundwater hydraulic gradients are also influenced by subsurface permeability contrasts, which can result in variability in flow direction and magnitude. In this study, we investigated the relationship between shallow groundwater flow direction, surface topography, and the subsurface topography of low permeability units in a headwater catchment at the Hubbard Brook Experimental Forest (HBEF), NH. We continuously monitored shallow groundwater levels in the solum throughout several seasons in a well network (20 wells of 0.18–1.1 m depth) within the upper hillslopes of Watershed 3 of the HBEF. Water levels were also monitored in four deeper wells, screened from 2.4 to 6.9 m depth within glacial drift of the C horizon. We conducted slug tests across the well network to determine the saturated hydraulic conductivity (Ksat) of the materials surrounding each well. Results showed that under higher water table regimes, groundwater flow direction mimics surface topography, but under lower water table regimes, flow direction can deviate as much as 56 degrees from surface topography. Under these lower water table conditions, groundwater flow direction instead followed the topography of the top of the C horizon. The interquartile range ofKsatwithin the C horizon was two orders of magnitude lower than within the solum. Overall, our results suggest that the land surface topography and the top of the C horizon acted as end members defining the upper and lower bounds of flow direction variability. This suggests that temporal dynamics of groundwater flow direction should be considered when calculating hydrologic fluxes in critical zone and runoff generation studies of headwater catchments that are underlain by glacial drift. 
    more » « less
  3. The hydrology of alpine and subalpine areas in the Colorado Front Range (USA) is evolving, driven by warming and by the alteration of precipitation patterns, the timing of snowmelt, and other components of the hydrologic budget. Field measurements of soil hydraulic conductivity and moisture along 30-m transects (n = 13) of representative soils developed in surficial deposits and falling head slug tests of shallow groundwater in till demonstrate that hydraulic conductivity in the soil is comparable to hydraulic conductivity values in the shallow aquifer. Soil hydraulic conductivity values were variable (medians ranged from 5.6 × 10−7 to 4.96 × 10−5 m s−1) and increased in alpine areas underlain by periglacial deposits. Hydraulic conductivities measured by a modified Hvorslev technique in test wells ranged from 4.86 × 10−7 to 1.77 × 10−4 m s−1 in subalpine till. The results suggest a gradient from higher hydraulic conductivity in alpine zones, where short travel paths through periglacial deposits support ephemeral streams and wetlands, to lower hydraulic conductivity in the till-mantled subalpine zone. In drier downstream areas, streambed infiltration contributes substantially to near-channel groundwater. As summer temperatures and evapotranspiration (ET) increase and snowmelt occur earlier, alpine soils are likely to become more vulnerable to drought, and groundwater levels in the critical zone may lower, affecting the connectivity between late-melting snow, meltwater streams, and the areas they affect downstream. 
    more » « less
  4. Abstract In land surface models (LSMs), the hydraulic properties of the subsurface are commonly estimated according to the texture of soils at the Earth's surface. This approach ignores macropores, fracture flow, heterogeneity, and the effects of variable distribution of water in the subsurface oneffectivewatershed‐scale hydraulic variables. Using hydrograph recession analysis, we empirically constrain estimates of watershed‐scale effective hydraulic conductivities (K) and effective drainable aquifer storages (S) of all reference watersheds in the conterminous United States for which sufficient streamflow data are available (n = 1,561). Then, we use machine learning methods to model these properties across the entire conterminous United States. Model validation results in high confidence for estimates of log(K) (r2 > 0.89; 1% < bias < 9%) and reasonable confidence forS(r2 > 0.83; −70% < bias < −18%). Our estimates of effectiveKare, on average, two orders of magnitude higher than comparable soil‐texture‐based estimates of averageK, confirming the importance of soil structure and preferential flow pathways at the watershed scale. Our estimates of effectiveScompare favorably with recent global estimates of mobile groundwater and are spatially heterogeneous (5–3,355 mm). Because estimates ofSare much lower than the global maximums generally used in LSMs (e.g., 5,000 mm in Noah‐MP), they may serve both to limit model spin‐up time and to constrain model parameters to more realistic values. These results represent the first attempt to constrain estimates of watershed‐scale effective hydraulic variables that are necessary for the implementation of LSMs for the entire conterminous United States. 
    more » « less
  5. Abstract Rock weathering impacts atmospheric CO2levels with silicate rock dissolution removing CO2,and carbonate dissolution, pyrite oxidation, and organic rock carbon oxidation producing CO2. Glacierization impacts the hydrology and geomorphology of catchments and glacier retreat due to warming can increase runoff and initiate landscape succession. To investigate the impact of these changes on catchment scale weathering CO2balances, we report monthly samples of solute chemistry and continuous discharge records for a sequence of glacierized watersheds draining into Kachemak Bay, Alaska. We partition solute and acid sources and estimate inorganic weathering CO2balances using an inverse geochemical mixing model. Furthermore, we investigated how solutes vary with discharge conditions utilizing a concentration‐runoff framework. We develop an analogous fraction‐runoff framework which allows us to investigate changes in weathering contributions at different flows. Fraction‐runoff relationships suggest kinetic limitations on all reactions in glacierized catchments, and only silicate weathering in less glacierized catchments. Using forest cover as a proxy for landscape age and stability, multiple linear regression shows that faster reactions (pyrite oxidation) contribute less to the solute load with increasing forest cover, whereas silicate weathering (slow reaction kinetics) contributes more. Overall, in glacierized catchments, we find elevated weathering fluxes at high runoff despite significant dilution effects. This makes flux estimates that account for dilution more important in glacierized catchments. Our findings quantify how glaciers modify the inorganic weathering CO2balance of catchments through hydrologic and geomorphic forcings, and support the previous hypothesis that deglaciation will be accompanied by a shift in inorganic weathering CO2balances. 
    more » « less