Sleepers River Research Watershed Research Articles

Controls on decadal, annual, and seasonal concentration‐discharge relationships in the Sleepers River Research Watershed, Vermont, northeastern United States

AbstractPast studies on concentration‐discharge (C‐Q) relationships have focused on short‐term or low‐temporal resolution data. While advancing understanding of catchment processes, these studies provided limited insight on catchment response over time or to climate change. Using 15 solutes from 1992 to 2015 at Sleepers River Research Watershed, Vermont, we compared C‐Q relationships over decades, years, and seasons to elucidate controls on stream chemical variation. We applied end‐member mixing analysis (EMMA) to identify solute sources and flow path routing. EMMA identified three end‐members: near‐surface runoff (NSR), riparian groundwater, and hillslope hollow groundwater. Shifting mixing proportions of these end‐members accounted for the temporal variability of conservative (no chemical reaction en route from source to stream) solutes in streamflow. For example, an increase in NSR fraction, typical of high flow, caused flushing (increased concentrations) of NO3−, DOC, Al, and Fe, which were greatest in NSR, dilution of specific conductance and base cation, SO42−, Si, Sr, Ba, and Mn concentrations, which were greatest in the two groundwater end‐members. This behaviour is reflected in the b‐coefficient of the C‐Q relation (C = aQb), which indicates the strength of dilution (b < −0.1) and flushing (b > 0.1) effects. For conservative solutes, the b‐coefficient decreased significantly (p < 0.01) with an increase in the groundwater to NSR concentration ratio. Solutes that are conservative and have relatively constant concentrations in end‐members over time showed consistent annual C‐Q patterns over years and decades. Furthermore, the strength of dilution or flushing was stronger during the snowmelt period, when the NSR fraction peaked, than during the dormant and growing seasons. With shorter snowmelt periods and snow to rain shifts, the flushing or dilution power of snowmelt runoff will weaken and alter catchment response to climate change. These insights provide more tools for the interpretation of catchment processes and responses to climate change.

Hydrology and biogeochemistry datasets from Sleepers River Research Watershed, Danville, Vermont, USA

AbstractThe Sleepers River Research Watershed (SRRW) in Danville, Vermont, USA, was established in 1957 by the U.S Department of Agriculture, Agricultural Research Service (ARS). The 111‐km2 area was selected as a representative mixed‐land‐use upland landscape in the northeastern USA, with the objective to study fundamental hydrological processes. The U.S. Geological Survey (USGS) currently operates SRRW and is the fourth federal agency to lead it. SRRW has hosted historic research on snowmelt modelling and streamflow generation processes, and USGS has augmented this hydrological foundation with biogeochemical cycling research. This data note describes five freely available data releases in stable repositories, extending to the late 2010s, from the indicated start date: aqueous chemistry and isotopes (1991); precipitation, water temperature, and streamflow (1991); groundwater levels (1991); ground frost depth (1983); and snow depth and water equivalent (1960). USGS research at SRRW is expected to continue, and we plan to update these releases every 1–3 years.

Drivers of Dissolved Organic Carbon Mobilization From Forested Headwater Catchments: A Multi Scaled Approach

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.

Soil Aggregates as a Source of Dissolved Organic Carbon to Streams: An Experimental Study on the Effect of Solution Chemistry on Water Extractable Carbon

Over the past two decades, headwater streams of the northern hemisphere have shown increased amounts of dissolved organic carbon (DOC), coinciding with decreased acid deposition. The exact nature of the mechanistic link between precipitation composition and stream water DOC is still widely debated in the literature. We hypothesize that soil aggregates are the main source of stream water DOC and that DOC release is greater in organic rich, riparian soils versus hillslope soils. To test these hypotheses, we collected soils from two main landscape positions (hillslope and riparian zones) from the acid-impacted Sleepers River Research Watershed in northeastern Vermont. We performed aqueous soil extracts with solutions of different ionic strength (IS) and composition to simulate changes in soil solution. We monitored dynamic changes in soil particle size, aggregate architecture and composition, leachate DOC concentrations, dissolved organic matter (DOM) characteristics with fluorescence spectroscopy and trends in bioavailability. In low IS solutions, extractable DOC concentrations were significantly higher, particle size (by laser diffraction) was significantly smaller and organic material was separated from mineral particles in scanning electron microscope observations. Furthermore, higher DOC concentrations were found in Na+ compared to Ca2+ solutions of the same IS. These effects are attributed to aggregate dispersion due to expanding diffuse double layers in decreased IS solutions and to decreased bridging by divalent cations. Landscape position impacted quality but not quantity of released DOC. Overall, these results indicate that soil aggregates might be one important link between Critical Zone inputs (i.e. precipitation) and exports in streams.

Does Stream Water Composition at Sleepers River in Vermont Reflect Dynamic Changes in Soils During Recovery From Acidification?

Stream water pH and composition are widely used to monitor ongoing recovery from the deposition of strong anthropogenic acids in many forested headwater catchments in the northeastern US. However, stream water composition is a function of highly complex and coupled processes, flowpaths, and variations in soil and bedrock composition. Spatial heterogeneity is especially pronounced in headwater catchments with steep topography, potentially limiting stream water composition as an indicator of changes in critical zone (CZ) dynamics during system recovery. To investigate the link between catchment characteristics, landscape position, and stream water composition we used long-term data (1991-2015) from the Sleepers River Research Watershed (SRRW) in northeastern Vermont. We investigated trends with time in stream water and trends with time, depth, and landscape position (upslope, midslope, and riparian zone) in groundwater (GW) and soil solution. We further determined soil elemental composition and mineralogy on archived (1996) and modern (2017) soil samples to assess changes in composition with time. SRRW is inherently well-buffered by calcite in bedrock and till, but soils had become acidified and are now recovering from acidification. Although base cations, especially Ca, decrease progressively with time in GW, riparian soils have become more enriched in Ca, due to a mixture of lateral and vertical transfers. At the same time stream water Ca fluxes increased over the past two decades, likely due to the leaching of (transient) legacy Ca from riparian zones and increased water fluxes. The stream water response therefore reflects the dynamic changes in soil chemistry, flow routing and water inputs.

Approaches to stream solute load estimation for solutes with varying dynamics from five diverse small watersheds

AbstractEstimating streamwater solute loads is a central objective of many water‐quality monitoring and research studies, as loads are used to compare with atmospheric inputs, to infer biogeochemical processes, and to assess whether water quality is improving or degrading. In this study, we evaluate loads and associated errors to determine the best load estimation technique among three methods (a period‐weighted approach, the regression‐model method, and the composite method) based on a solute's concentration dynamics and sampling frequency. We evaluated a broad range of varying concentration dynamics with stream flow and season using four dissolved solutes (sulfate, silica, nitrate, and dissolved organic carbon) at five diverse small watersheds (Sleepers River Research Watershed, VT; Hubbard Brook Experimental Forest, NH; Biscuit Brook Watershed, NY; Panola Mountain Research Watershed, GA; and Río Mameyes Watershed, PR) with fairly high‐frequency sampling during a 10‐ to 11‐yr period. Data sets with three different sampling frequencies were derived from the full data set at each site (weekly plus storm/snowmelt events, weekly, and monthly) and errors in loads were assessed for the study period, annually, and monthly. For solutes that had a moderate to strong concentration–discharge relation, the composite method performed best, unless the autocorrelation of the model residuals was <0.2, in which case the regression‐model method was most appropriate. For solutes that had a nonexistent or weak concentration–discharge relation (model R2 < about 0.3), the period‐weighted approach was most appropriate. The lowest errors in loads were achieved for solutes with the strongest concentration–discharge relations. Sample and regression model diagnostics could be used to approximate overall accuracies and annual precisions. For the period‐weighed approach, errors were lower when the variance in concentrations was lower, the degree of autocorrelation in the concentrations was higher, and sampling frequency was higher. The period‐weighted approach was most sensitive to sampling frequency. For the regression‐model and composite methods, errors were lower when the variance in model residuals was lower. For the composite method, errors were lower when the autocorrelation in the residuals was higher. Guidelines to determine the best load estimation method based on solute concentration–discharge dynamics and diagnostics are presented, and should be applicable to other studies.

Water's Way at Sleepers River watershed – revisiting flow generation in a post‐glacial landscape, Vermont USA

AbstractThe Sleepers River Research Watershed (SRRW) in Vermont, USA, has been the site of active hydrologic research since 1959 and was the setting where Dunne and Black demonstrated the importance and controls of saturation‐excess overland flow (SOF) on streamflow generation. Here, we review the early studies from the SRRW and show how they guided our conceptual approach to hydrologic research at the SRRW during the most recent 25 years. In so doing, we chronicle a shift in the field from early studies that relied exclusively on hydrometric measurements to today's studies that include chemical and isotopic approaches to further elucidate streamflow generation mechanisms. Highlights of this evolution in hydrologic understanding include the following: (i) confirmation of the importance of SOF to streamflow generation, and at larger scales than first imagined; (ii) stored catchment water dominates stream response, except under unusual conditions such as deep frozen ground; (iii) hydrometric, chemical and isotopic approaches to hydrograph separation yield consistent and complementary results; (iv) nitrate and sulfate isotopic compositions specific to atmospheric inputs constrain new water contributions to streamflow; and (v) convergent areas, or ‘hillslope hollows’, contribute disproportionately to event hydrographs. We conclude by summarizing some remaining challenges that lead us to a vision for the future of research at the SRRW to address fundamental questions in the catchment sciences. Copyright © 2014 John Wiley & Sons, Ltd.

Responses of stream nitrate and DOC loadings to hydrological forcing and climate change in an upland forest of the northeastern United States

In coming decades, higher annual temperatures, increased growing season length, and increased dormant season precipitation are expected across the northeastern United States in response to anthropogenic forcing of global climate. We synthesized long‐term stream hydrochemical data from the Sleepers River Research Watershed in Vermont, United States, to explore the relationship of catchment wetness to stream nitrate and DOC loadings. We modeled changes in growing season length and precipitation patterns to simulate future climate scenarios and to assess how stream nutrient loadings respond to climate change. Model results for the 2070–2099 time period suggest that stream nutrient loadings during both the dormant and growing seasons will respond to climate change. During a warmer climate, growing season stream fluxes (runoff +20%, nitrate +57%, and DOC +58%) increase as more precipitation (+28%) and quick flow (+39%) occur during a longer growing season (+43 days). During the dormant season, stream water and nutrient loadings decrease. Net annual stream runoff (+8%) and DOC loading (+9%) increases are commensurate with the magnitude of the average increase of net annual precipitation (+7%). Net annual stream water and DOC loadings are primarily affected by increased dormant season precipitation. In contrast, decreased annual loading of stream nitrate (−2%) reflects a larger effect of growing season controls on stream nitrate and the effects of lengthened growing seasons in a warmer climate. Our findings suggest that leaching of nitrate and DOC from catchment soils will be affected by anthropogenic climate forcing, thereby affecting the timing and magnitude of annual stream loadings in the northeastern United States.

Seasonal and event variations in δ34S values of stream sulfate in a Vermont forested catchment: Implications for sulfur sources and cycling

Stable sulfur (S) isotope ratios can be used to identify the sources of sulfate contributing to streamwater. We collected weekly and high-flow stream samples for S isotopic analysis of sulfate through the entire water year 2003 plus the snowmelt period of 2004. The study area was the 41-ha forested W-9 catchment at Sleepers River Research Watershed, Vermont, a site known to produce sulfate from weathering of sulfide minerals in the bedrock. The δ 34S values of streamwater sulfate followed an annual sinusoidal pattern ranging from about 6.5‰ in early spring to about 10‰ in early fall. During high-flow events, δ 34S values typically decreased by 1 to 3‰ from the prevailing seasonal value. The isotopic evidence suggests that stream sulfate concentrations are controlled by: (1) an overall dominance of bedrock-derived sulfate ( δ 34S ~ 6–14‰); (2) contributions of pedogenic sulfate ( δ 34S ~ 5–6‰) during snowmelt and storms with progressively diminishing contributions during base flow recession; and (3) minor effects of dissimilatory bacterial sulfate reduction and subsequent reoxidation of sulfides. Bedrock should not be overlooked as a source of S in catchment sulfate budgets.

Carbon isotope fractionation of dissolved inorganic carbon (DIC) due to outgassing of carbon dioxide from a headwater stream

AbstractThe stable isotopic composition of dissolved inorganic carbon (δ13C‐DIC) was investigated as a potential tracer of streamflow generation processes at the Sleepers River Research Watershed, Vermont, USA. Downstream sampling showed δ13C‐DIC increased between 3–5‰ from the stream source to the outlet weir approximately 0·5 km downstream, concomitant with increasing pH and decreasing PCO2. An increase in δ13C‐DIC of 2·4 ± 0·1‰ per log unit decrease of excess PCO2 (stream PCO2 normalized to atmospheric PCO2) was observed from downstream transect data collected during snowmelt. Isotopic fractionation of DIC due to CO2 outgassing rather than exchange with atmospheric CO2 may be the primary cause of increased δ13C‐DIC values downstream when PCO2 of surface freshwater exceeds twice the atmospheric CO2 concentration. Although CO2 outgassing caused a general increase in stream δ13C‐DIC values, points of localized groundwater seepage into the stream were identified by decreases in δ13C‐DIC and increases in DIC concentration of the stream water superimposed upon the general downstream trend. In addition, comparison between snowmelt, early spring and summer seasons showed that DIC is flushed from shallow groundwater flowpaths during snowmelt and is replaced by a greater proportion of DIC derived from soil CO2 during the early spring growing season. Thus, in spite of effects from CO2 outgassing, δ13C of DIC can be a useful indicator of groundwater additions to headwater streams and a tracer of carbon dynamics in catchments. Copyright © 2007 John Wiley & Sons, Ltd.

Mercury and Organic Carbon Dynamics During Runoff Episodes from a Northeastern USA Watershed

Mercury and organic carbon concentrations vary dynamically in streamwater at the Sleepers River Research Watershed in Vermont, USA. Total mercury (THg) concentrations ranged from 0.53 to 93.8 ng/L during a 3-year period of study. The highest mercury (Hg) concentrations occurred slightly before peak flows and were associated with the highest organic carbon (OC) concentrations. Dissolved Hg (DHg) was the dominant form in the upland catchments; particulate Hg (PHg) dominated in the lowland catchments. The concentration of hydrophobic acid (HPOA), the major component of dissolved organic carbon (DOC), explained 41–98% of the variability of DHg concentration while DOC flux explained 68–85% of the variability in DHg flux, indicating both quality and quantity of the DOC substantially influenced the transport and fate of DHg. Particulate organic carbon (POC) concentrations explained 50% of the PHg variability, indicating that POC is an important transport mechanism for PHg. Despite available sources of DHg and wetlands in the upland catchments, dissolved methylmercury (DmeHg) concentrations in streamwaters were below detection limit (0.04 ng/L). PHg and particulate methylmercury (PmeHg) had a strong positive correlation (r 2 = 0.84, p < 0.0001), suggesting a common source; likely in-stream or near-stream POC eroded or re-suspended during spring snowmelt and summer storms. Ratios of PmeHg to THg were low and fairly constant despite an apparent higher methylmercury (meHg) production potential in the summer. Methylmercury production in soils and stream sediments was below detection during snowmelt in April and highest in stream sediments (compared to forest and wetland soils) sampled in July. Using the watershed approach, the correlation of the percent of wetland cover to TmeHg concentrations in streamwater indicates that poorly drained wetland soils are a source of meHg and the relatively high concentrations found in stream surface sediments in July indicate these zones are a meHg sink.

A carbon-based method for estimating the wetness of forest surface soil horizons

The degree of wetness in forest surface soils has an effect on chemical and biological processes but is not easily measured. The high spatial variability in carbon (C) concentration creates high variability in water-holding capacity, and gravimetric water content is not informative. Local hydrology can create patchiness in soil moisture, with saturated soils often found near well-drained ones. When sampling to measure such factors as nitrification potential, it would be advantageous to have a simple metric that reflects the relative wetness of the soil. The relationship between C concentration (range 51.5–520.8 g·kg–1) and gravimetric water content was found to be linear for a set of 113 H- and A-horizon samples assumed to be at field capacity. The wetness ratio is defined as the actual water content of a sample divided by the water content predicted by the least squares regression equation based on C concentration (soil water content (kg·kg–1) = 0.080 + 0.0057 soil C concentration (g·kg–1)). Soil moisture retention curves were developed for a small number of samples in the range of 0 to about –10 kPa and showed that the equation predicted that water would be held at relatively high potential. In samples taken from 10 watersheds in the northeastern USA, wetness ratios between 1.25 and 3.1 were associated with soils identified in the field as ranging from wet to boglike. A median ratio of 0.49 was found in a watershed sampled after an extended dry period. At the Sleepers River Research Watershed, high wetness ratios were associated with a high soil calcium concentration, presumably from enriched groundwater. The ratio should be a useful measurement in watershed studies.

Tracing Sources of Streamwater Sulfate During Snowmelt Using S and O Isotope Ratios of Sulfate and 35S Activity

The biogeochemical cycling of sulfur (S) was studied during the 2000 snowmelt at Sleepers River Research Watershed in northeastern Vermont, USA using a hydrochemical and multi-isotope approach. The snowpack and 10 streams of varying size and land use were sampled for analysis of anions, dissolved organic carbon (DOC), 35S activity, and δ34S and δ18O values of sulfate. At one of the streams, δ18O values of water also were measured. Apportionment of sulfur derived from atmospheric and mineral sources based on their distinct δ34S values was possible for 7 of the 10 streams. Although mineral S generally dominated, atmospheric-derived S contributions exceeded 50% in several of the streams at peak snowmelt and averaged 41% overall. However, most of this atmospheric sulfur was not from the melting snowpack; the direct contribution of atmospheric sulfate to streamwater sulfate was constrained by 35S mass balance to a maximum of 7%. Rather, the main source of atmospheric sulfur in streamwater was atmospheric sulfate deposited months to years earlier that had microbially cycled through the soil organic sulfur pool. This atmospheric/pedospheric sulfate (pedogenic sulfate formed from atmospheric sulfate) source is revealed by δ18O values of streamwater sulfate that remained constant and significantly lower than those of atmospheric sulfate throughout the melt period, as well as streamwater 35S ages of hundreds of days. Our results indicate that the response of streamwater sulfate to changes in atmospheric deposition will be mediated by sulfate retention in the soil.

Tracing sources of nitrate in snowmelt runoff using a high‐resolution isotopic technique

The denitrifier method to determine the dual isotopic composition (δ15N and δ18O) of nitrate is well suited for studies of nitrogen contributions to streams during runoff events. This method requires only 70 nmol of NO3− and enables high throughput of samples. We studied nitrate sources to a headwater stream during snowmelt by generating a high‐temporal resolution dataset at the Sleepers River Research Watershed in Vermont, USA. In the earliest phase of runoff, stream NO3− concentrations were highest and stream discharge, NO3− concentrations, and δ18O of NO3− generally tracked one another during diurnal melting. The isotopic composition of stream NO3− varied in‐between atmospheric and groundwater NO3− end members indicating a direct contribution of atmospherically‐derived NO3− from the snow pack to the stream. During the middle to late phases of snowmelt, the source shifted toward soil NO3− entering the stream via shallow subsurface flow paths.

Controls on old and new water contributions to stream flow at some nested catchments in Vermont, USA

AbstractFactors controlling the partitioning of old and new water contributions to stream flow were investigated for three events in four catchments (three of which were nested) at Sleepers River Research Watershed in Danville, Vermont. In the 1993 snowmelt period, two‐component isotopic hydrograph separations showed that new water (meltwater) inputs to the stream ranged widely from 41 to 74%, and increased with catchment size (41 to 11 125 ha) (with one exception) and with open land cover (0–73%). Peak dissolved organic carbon concentrations and relative alkalinity dilution in stream water ranked in the same order among catchments as the new water fractions, suggesting that new water followed shallow flow paths. During the 1994 snowmelt, despite similar timing and magnitude of melt inputs, the new‐water contribution to stream flow ranged only from 30 to 36% in the four catchments. We conclude that the uncommonly high and variable new water fractions in streamwater during the 1993 melt were caused by direct runoff of meltwater over frozen ground, which was prevalent in open land areas during the 1993 winter. In a high‐intensity summer rainstorm in 1993, new water fractions were smaller relative to the 1993 snowmelt, ranging from 28 to 46%, but they ranked in the identical catchment order. Reconciliation of the contrasting patterns of new–old water partitioning in the three events appears to require an explanation that invokes multiple processes and effects, including: topographically controlled increase in surface‐saturated area with increasing catchment size; direct runoff over frozen ground; low infiltration in agriculturally compacted soils; differences in soil transmissivity, which may be more relevant under dry antecedent conditions. These data highlight some of the difficulties faced by catchment hydrologists in formulating a theory of runoff generation at varying basin scales. Copyright © 2002 John Wiley &amp; Sons, Ltd.

Spectral albedo/reflectance of littered forest snow during the melt season

AbstractDespite the importance of litter on forest floor albedo and brightness, previous studies have not documented forest floor albedo or litter cover in any detail. Our objective was to describe the seasonal influence of litter on spectral albedos and nadir reflectances of a forest snowpack in a mixed‐hardwood stand in the Sleepers River Research Watershed (SRRW) in Danville, Vermont (37°39′ N, 119°2′ W). Experimental measurements in a nearby open area at the Snow Research Station of the SRRW nearly duplicated the spectral trend observed in the forest. Spectral albedo and nadir reflectance measurements in the visible and near infrared (350–2500 nm) transitioned from a gently curved shape through the visible range (for finer‐grained, lightly littered snow) to one having a peak in the red/near‐infrared (near 760 nm) as the snowmelt season progressed (for coarser‐grained, more heavily littered snow). The snowpack became optically thin as surface litter reached high percentages. A point‐in‐time digital photographic survey of the late‐lying snowpacks of three forest stands and the open showed that median litter cover percentages in the coniferous, deciduous, mixed‐forest, and an open area were 17·5, 6·1, 1·2, and 0·04 respectively. A Kruskal–Wallis ANOVA on ranks and pairwise comparisons using Dunn's test indicated that the litter covers of the three forest stands were significantly different with &gt;95% confidence. The snowpack was relatively shallow (&lt;1 m), as is typical for this area of Vermont. From a remote‐sensing standpoint, and since shallow snow and increased grain size also lower the visible albedo, we can expect that snowpack litter will cause decreased albedo earlier in the snowmelt season, at deeper snow depths, and will tend to shift the maximum albedo peak to the red/NIR range as the melt season progresses. Published in 2001 by John Wiley &amp; Sons, Ltd.

The effect of frozen soil on snowmelt runoff at Sleepers River, Vermont

Soil frost depth has been monitored at the Sleepers River Research Watershed in northeastern Vermont since 1984. Soil frost develops every winter, particularly in open fields, but its depth varies greatly from year to year in inverse relation to snow depth. During the 15 years of record at a benchmark mid-elevation open site, the annual maximum frost depth varied from 70 to 390 mm. We empirically tested the hypothesis that frozen soil prevents infiltration and recharge, thereby causing an increased runoff ratio (streamflow/(rain + snowmelt)) during the snowmelt hydrograph rise and a decreased runoff ratio during snowmelt recession. The hypothesis was not supported at the 111 km 2 W-5 catchment; there was no significant correlation of the runoff ratio with the seasonal maximum frost depth for either the pre-peak or post-peak period. In an analysis of four events, however, the presence of frost promoted a large and somewhat quicker response to rainfall relative to the no-frost condition, although snow cover caused a much greater time-to-peak regardless of frost status. For six years of flow and frost depth measured at the 59 ha agricultural basin W-2, the hypothesis appeared to be supported. The enhancement of runoff due to soil frost is evident on small plots and in extreme events, such as rain on frozen snow-free soil. In the northeastern USA rind eastern Canada, the effect is often masked in larger catchments by several confounding factors, including storage of meltwater in the snowpack, variability in snowmelt timing due to elevational and aspect differences, interspersed forested land where frost may be absent, and the timing of soil thawing relative to the runoff peak.

Air permeability and capillary rise as measures of the pore structure of snow: an experimental and theoretical study

Air permeability and capillary pressure are macroscopic snow properties that are influenced by the pore structure of the snow cover. Formulas for predicting fluid transport, species elution, and acoustive wave propagation require parameterization of one or both of these properties. We report paired measurements of permeability and capillary rise from snow samples at field sites in Hanover, New Hampshire, and Sleepers River Research Watershed, Danville, Vermont. We augment these data with laboratory tests on sieved snow and glass beads. Our measurements demonstrate a linear relationship between permeability and the ratio of porosity and the square of capillary rise, which we corroborate theoretically using a simple conduit model of the pore space. We propose that scatter in the data results, in part, from the effect of crystal shape on air flow and imbibition contact angle. Since the early measurements and classification schemes of Bader in 1939, many investigators have expanded the database of permeability observations for a wide range of snow types. We summarize these data and report our own recent observations from the New England sites and from an additional site in Manitoba, Canada. Our measurements are in the high range of reported values. However, after normalizing our data by the square of grain diameter, they follow the empirical function of Shimizu fairly closely. This agreement supports our measurements, and demonstrates the usefulness of Shimizu's function for snow types other than the relatively dense, fine-grained snow used in his analysis. Our normalized permeability data for low density snow, as well as the Shimizu function, are below theoretical predictions for suspensions of spheres and infinite cylinders. By extending the model for spheres to oblate spheroids and discs, we estimate permeability that is in closer agreement with our data. We suggest that a decrease in surface-to-volume ratio as snow ages may account for a relative increase in normalized permeability. Copyright © 1999 John Wiley & Sons, Ltd.

Cross‐Well Slug Testing in Unconfined Aquifers: A Case Study from the Sleepers River Watershed, Vermont

AbstractNormally, slug test measurements are limited to the well in which the water level is perturbed. Consequently, it is often difficult to obtain reliable estimates of hydraulic properties, particularly if the aquifer is anisotropic or if there is a wellbore skin. In this investigation, we use partially penetrating stress and observation wells to evaluate specific storage, radial hydraulic conductivity and anisotropy of the aquifer, and the hydraulic conductivity of the borehole skin. The study site is located in the W9 subbasin of the Sleepers River Research Watershed, Vermont. At the site, ∼3 m of saturated till are partially penetrated by a stress well located in the center of the unconfined aquifer and six observation wells located above, below, and at the depth of the stress well at radial distances of 1.2 and 2.4 m. The observation wells were shut in with inflatable packers.The semianalytical solution of Butler (1995) was used to conduct a sensitivity analysis and to interpret slug test results. The sensitivity analysis indicates that the response of the stress well is primarily sensitive to radial hydraulic conductivity, less sensitiive to anisotropy and the conductivity of the borehole skin, and nearly insensitive to specific storage. In contrast, the responses of the observation wells are sensitive to all four parameters. Interpretation of the field data was facilitated by generating type curves in a manner analogous to the method of Cooper et al. (1967). Because the value of radial hydraulic conductivty is obtained from a match point, the number of unknowns is reduced to three. The estimated values of radial hydraulic conductivity and specific storage are comparable to those derived from the methods of Bouwer and Rice (1976) and Cooper et al. (1967). The values and skin conductivity, however, could not have been obtained without the use of observation wells.

Effects of Subbasin Size on Topographic Characteristics and Simulated Flow Paths in Sleepers River Watershed, Vermont

The effects of subbasin size on topographic characteristics and simulated flow paths were determined for the 111.5‐km2 Sleepers River Research Watershed in Vermont using the watershed model TOPMODEL. Topography is parameterized in TOPMODEL as the spatial and statistical distribution of the index ln (a/tan B), where In is the Napierian logarithm, a is the upslope area per unit contour length, and tan B is the slope gradient. The mean, variance, and skew of the ln (a/tan B) distribution were computed for several sets of nested subbasins (0.05 to 111.5 km2)) along streams in the watershed and used as input to TOPMODEL. In general, the statistics of the ln (a/tan B) distribution and the simulated percentage of overland flow in total streamflow increased rapidly for some nested subbasins and decreased rapidly for others as subbasin size increased from 0.05 to 1 km2, generally increased up to a subbasin size of 5 km2, and remained relatively constant at a subbasin size greater than 5 km2. Differences in simulated flow paths among subbasins of all sizes (0.05 to 111.5 km2) were caused by differences in the statistics of the ln (a/tan B) distribution, not by differences in the explicit spatial arrangement of ln (a/tan B) values within the subbasins. Analysis of streamflow chemistry data from the Neversink River watershed in southeastern New York supports the hypothesis that subbasin size affects flow‐path characteristics.