[1] Spatial and temporal trends in stream chemistry were investigated in a large (1600 km 2 ) alpine watershed in the southern Rocky Mountains of Colorado to help understand mechanisms of streamflow generation. We observed linear increases of concentrations of chemical constituents in streamflow as accumulated drainage area increased along the main channel of Saguache Creek. We tested two conceptual models of streamflow generation against our stream chemistry observations. One model is essentially two-dimensional and treats streamflow generation at the large watershed scale as the aggregation of runoff responses from individual hillslopes, primarily surface and shallow subsurface flow paths. Alternatively, a fully three-dimensional conceptual model treats streamflow generation as being controlled by a distribution of large-scale groundwater flow paths as well as surface and shallow subsurface flow paths. The structure and magnitude of groundwater contributions in streamflow as a function of increasing scale provided a key distinction between these two conceptual models. End-member mixing analysis and measurements of hydraulic head gradients in streambeds were used to quantify basin-scale groundwater contributions to streamflow with increasing spatial scale in the Saguache Creek watershed. Our data show that groundwater contributions are important in streamflow generation at all scales and, more importantly, that groundwater contributions to streamflow do increase with increasing watershed scale. These results favor the three-dimensional conceptual model in which long groundwater flow paths provide a streamflow generation process at large scales that is not operative at smaller scales. This finding indicates that large watersheds may be more than simply the aggregation of hillslopes and small catchments.

Soil moisture distributions are expected to be closely tied to ecosystem processes in water-limited environments of the southwest United States. Nevertheless, few studies have addressed how soil moisture varies across grassland to forest transitions frequently observed in semiarid mountain settings. In this study, we quantify the vegetation controls on surface soil moisture by sampling a range of different ecosystems present in the Valles Caldera, New Mexico. Soil and atmospheric variables were measured during a 2-week field campaign conducted in late July to early August 2005 during the North American monsoon. Field observations were supplemented by a network of continuous instruments used to assess conditions prior to and after the sampling campaign. Results reveal that soil moisture responds directly to summer precipitation events and is mediated by plant interception, which differs across the grassland-forest continuum. The nature of the spatial and temporal variations in soil moisture changes across the different sampled ecosystems: wetlands, riparian forests, grasslands, ponderosa, deciduous and mixed conifer forests. In particular, statistical analyses of soil moisture distributions indicate that distinct regimes (e.g. probability density functions) exist along the semiarid vegetation gradient, which may not be revealed through simple metrics such as the ecosystem average. Ecosystem differences are further elucidated through comparison of the spatial variations in each vegetation type, indicating higher variability in wetland and grassland sites.

[1] Residence times provide vital information on hydrological, geochemical, and ecological processes in watersheds. The common perception is that mean residence times in watersheds are very short, on the order of days to years. However, there is growing concern that longer residence times of centuries to millennia are not being captured by traditional surface water age-dating methods. We hypothesize that if mean residence times are biased short, then weathering rates calculated from mean residence times will be forced unrealistically high to match observed solute concentrations. We test this hypothesis by calculating weathering rates from springs based upon residence times estimated using three different age-dating methods. Observed solute concentrations require unrealistically large weathering rates if typical short residence times are employed as compared to rates derived from longer residence times. Residence time distributions in watersheds have substantially longer tails than previously thought, with implications for age-dating methods and their interpretation to infer process behavior.

No abstract

Estimates of groundwater circulation depths based on field data are lacking. These data are critical to inform and refine hydrogeologic models of mountainous watersheds, and to quantify depth and time dependencies of weathering processes in watersheds. Here we test two competing hypotheses on the role of geology and geologic setting in deep groundwater circulation and the role of deep groundwater in the geochemical evolution of streams and springs. We test these hypotheses in two mountainous watersheds that have distinctly different geologic settings (one crystalline, metamorphic bedrock and the other volcanic bedrock). Estimated circulation depths for springs in both watersheds range from 0.6 to 1.6 km and may be as great as 2.5 km. These estimated groundwater circulation depths are much deeper than commonly modeled depths suggesting that we may be forcing groundwater flow paths too shallow in models. In addition, the spatial relationships of groundwater circulation depths are different between the two watersheds. Groundwater circulation depths in the crystalline bedrock watershed increase with decreasing elevation indicative of topography‐driven groundwater flow. This relationship is not present in the volcanic bedrock watershed suggesting that both the source of fracturing (tectonic versus volcanic) and increased primary porosity in the volcanic bedrock play a role in deep groundwater circulation. The results from the crystalline bedrock watershed also indicate that relatively deep groundwater circulation can occur at local scales in headwater drainages less than 9.0 km2 and at larger fractions than commonly perceived. Deep groundwater is a primary control on streamflow processes and solute concentrations in both watersheds.

Interbasin groundwater flow (IGF) can play a significant role in the generation and geochemical evolution of streamflow. However, it is exceedingly difficult to identify IGF and to determine the location and quantity of water that is exchanged between watersheds. How does IGF affect landscape/watershed geomorphic evolution? Can geomorphic metrics be used to identify the presence of IGF? We examine these questions in two adjacent sedimentary watersheds in northern New Mexico using a combination of geomorphic/landscape metrics, springflow residence times, and spatial geochemical patterns. IGF is expressed geomorphically in the landscape placement of springs and flow direction and shape of stream channels. Springs emerge preferentially on one side of stream valleys where landscape incision has intercepted IGF flow paths. Stream channels grow toward the IGF source and show little bifurcation. In addition, radiocarbon residence times of springs decrease and the geochemical composition of springs changes as the connection to IGF is lost.

Spatial and temporal trends in stream chemistry have been used to provide insights into the scale dependencies of streamflow generation processes in small catchments. However, these scale dependencies have not been thoroughly investigated at large watershed scales (defined as drainage areas greater than 1000 km2). Quantifying these scale dependencies is critical to understanding how large watersheds will respond to future perturbations; e.g., the long‐term streamflow response to climate change and/or changes in land‐cover and land‐use. Here we investigate the spatial and temporal scaling relationships of all dominant streamflow generation processes in a large alpine watershed in the southern Rocky Mountains of Colorado. Observations in the watershed indicate that dominant streamflow processes are spatially and temporally variable. The relative strengths of dominant streamflow mechanisms vary as a function of internal watershed structure (i.e., spatial variability in topographic relief, soil development, groundwater flowpath development, and stream network structure) and external forcing such as timing and character of precipitation. This behavior coupled with previous observations that streamflow from the watershed contained a significant component of basin‐scale groundwater, suggests that similar large watersheds may have internal buffering, at least initially, against rapid change.

Abstract:Twelve modified passive capillary samplers (M-PCAPS) were installed in remote locations within a large, alpine watershed located in the southern Rocky Mountains of Colorado to collect samples of infiltration during the snowmelt and summer rainfall seasons. These samples were collected in order to provide better constraints on the isotopic composition of soil-water endmembers in the watershed. The seasonally integrated stable isotope composition (υ 18 O and υ 2 H) of soil-meltwater collected with M-PCAPS installed at shallow soil depths <10 cm was similar to the seasonally integrated isotopic composition of bulk snow taken at the soil surface. However, meltwater which infiltrated to depths >20 cm evolved along an isotopic enrichment line similar to the trendline described by the evolution of fresh snow to surface runoff from snowmelt in the watershed. Coincident changes in geochemistry were also observed at depth suggesting that the isotopic and geochemical composition of deep infiltration may be very different from that obtained by surface and/or shallow-subsurface measurements. The M-PCAPS design was also used to estimate downward fluxes of meltwater during the snowmelt season. Shallow and deep infiltration averaged 8Ð4 and 4Ð7 cm of event water or 54 and 33% of the measured snow water equivalent (SWE), respectively. Finally, dominant shallow-subsurface runoff processes occurring during snowmelt could be identified using geochemical data obtained with the M-PCAPS design. One soil regime was dominated by a combination of slow matrix flow in the shallow soil profile and fast preferential flow at depth through a layer of platy, volcanic rocks. The other soil regime lacked the rock layer and was dominated by slow matrix flow. Based on these results, the M-PCAPS design appears to be a useful, robust methodology to quantify soil-water fluxes during the snowmelt season and to sample the stable isotopic and geochemical composition of soil-meltwater endmembers in remote watersheds.