James P. McNamara

Storage as a Metric of Catchment Comparison

James P. McNamara et al. — Thu, 05 Jan 2012 15:34:42 PST

The volume of water stored within a catchment, and its partitioning among groundwater, soil moisture, snowpack, vegetation, and surface water are the variables that ultimately characterize the state of the hydrologic system. Accordingly, storage may provide useful metrics for catchment comparison. Unfortunately, measuring and predicting the amount of water present in a catchment is seldom done; tracking the dynamics of these stores is even rarer. Storage moderates fluxes and exerts critical controls on a wide range of hydrologic and biologic functions of a catchment. While understanding runoff generation and other processes by which catchments release water will always be central to hydrologic science, it is equally essential to understand how catchments retain water. We have initiated a catchment comparison exercise to begin assessing the value of viewing catchments from the storage perspective. The exercise is based on existing data from five watersheds, no common experimental design, and no integrated modelling efforts. Rather, storage was estimated independently for each site. This briefing presents some initial results of the exercise, poses questions about the definitions and importance of storage and the storage perspective, and suggests future directions for ongoing activities.

Small Soil Storage Capacity Limits Benefit of Winter Snowpack to Upland Vegetation

T. J. Smith et al. — Thu, 05 Jan 2012 15:34:39 PST

In the western United States, the mountain snowpack is an important natural reservoir that extends spring and summer water delivery to downstream users and ecosystems. The importance of winter snow accumulation to upland ecosystems is not as clearly defined. This study investigates the relative contribution of winter precipitation to upland spring and summer soil moisture storage and availability in a semi-arid mountainous watershed. At this site, coarse soil textures and shallow soil depths limit soil storage capacity to 6–16 cm. Winter precipitation exceeds soil storage capacity by 2.5 times. Accordingly, soil moisture profiles at most locations in the watershed reach field capacity in early winter. With soil storage near capacity, water released by snowmelt primarily contributes to deep drainage and makes a limited contribution to the soil moisture reservoir. Water that is retained by the soil after the snowpack melts is lost to evapotranspiration in as little as 10 days. In contrast, spring precipitation extends moist soil conditions by up to 90 days into the warm season, when ecological water demand is highest. These field observations suggest that changes in spring precipitation, not winter snowpack, may have the greater impact on upland ecosystems in this environment. Furthermore, because winter precipitation is in excess compared to the soil storage capacity, soil moisture availability may be fairly insensitive to climate change-induced transitions from snow to rain.

Aspect Influences on Soil Water Retention and Storage

I. J. Geroy et al. — Thu, 05 Jan 2012 15:34:36 PST

Many catchment hydrologic and ecologic processes are impacted by the storage capacity of soil water, which is dictated by the profile thickness and water retention properties of soil. Soil water retention properties are primarily controlled by soil texture, which in turn varies spatially in response to microclimate-induced differences in insolation, wetness and temperature. All of these variables can be strongly differentiated by slope aspect. In this study, we compare quantitative measures of soil water retention capacity for two opposing slopes in a semi-arid catchment in southwest Idaho, USA. Undisturbed soil cores from north and south aspects were subjected to a progressive drainage experiment to estimate the soil water retention curve for each sample location. The relatively large sample size (35) supported statistical analysis of slope scale differences in soil water retention between opposing aspects. Soils on the north aspect retain as much as 25% more water at any given soil water pressure than samples from the south aspect slope. Soil porosity, soil organic matter and silt content were all greater on the north aspect, and each contributed to greater soil water retention. These results, along with the observation that soils on north aspect slopes tend to be deeper, indicate that north aspect slopes can store more water from the wet winter months into the dry summer in this region, an observation with potential implications on ecological function and landscape evolution.

A Simplified Approach for Estimating Soil Carbon and Nitrogen Stocks in Semi-Arid Complex Terrain

Melvin L. Kunkel et al. — Wed, 28 Sep 2011 16:04:37 PDT

We investigated soil carbon (C) and nitrogen (N) distribution and developed a model, using readily available geospatial data, to predict that distribution across a mountainous, semi-arid, watershed in southwestern Idaho (USA). Soil core samples were collected and analyzed from 133 locations at 6 depths (n=798), revealing that aspect dramatically influences the distribution of C and N, with north-facing slopes exhibiting up to 5 times more C and N than adjacent southfacing aspects. These differences are superimposed upon an elevation (precipitation) gradient, with soil C and N contents increasing by nearly a factor of 10 from the bottom (1100 m elevation) to the top (1900 m elevation) of the watershed. Among the variables evaluated, vegetation cover, as represented by a Normalized Difference Vegetation Index (NDVI), is the strongest, positively correlated, predictor of C; potential insolation (incoming solar radiation) is a strong, negatively correlated, secondary predictor. Approximately 62% (as R²) of the variance in the C data is explained using NDVI and potential insolation, compared with an R² of 0.54 for a model using NDVI alone. Soil N is similarly correlated to NDVI and insolation. We hypothesize that the correlations between soil C and N and slope, aspect and elevation reflect, in part, the inhibiting influence of insolation on semi-arid ecosystem productivity via water limitation. Based on these identified relationships, two modeling techniques (multiple linear regression and cokriging) were applied to predict the spatial distribution of soil C and N across the watershed. Both methods produce similar distributions, successfully capturing observed trends with aspect and elevation. This easily applied approach may be applicable to other semi-arid systems at larger scales.

Seasonal Recharge Components in an Urban/Agricultural Mountain Front Aquifer System Using Noble Gas Thermometry

Michael J. Thoma et al. — Wed, 28 Sep 2011 16:04:35 PDT

Thirteen noble gas samples were collected from eleven wells and two mountain springs in the Treasure Valley, Idaho, USA to derive recharge temperatures using noble gas thermometry. One common assumption with noble gas thermometry is that recharge temperatures are roughly equal to the mean annual surface temperature. When water table depths are shallow or variable, or infiltration is seasonal recharge temperatures may be significantly different from the mean annual surface temperature. Water table depths throughout the study area were used to estimate recharge source temperatures using an infiltration-weighted recharge temperature model which takes into account a time-variable water table. This model was applied to six different seasonally-dependent recharge scenarios. The modeled recharge temperatures for all scenarios showed a strong dependence of recharge temperature on mean annual depth to water. Temperature results from the different recharge scenarios ranged from near the mean annual surface temperature to as much as 6 °C warmer. This compared well to noble gas derived recharge temperatures from the valley wells which ranged from 5 °C below to 7.4 °C above the mean annual surface temperature of the valley. Cooler temperatures suggest an influence of recharge through the adjacent mountain block while warmer temperatures suggest an influence from summer irrigation.

Modeling the Water and Energy Balance of Vegetated Areas with Snow Accumulation

T. J. Kelleners et al. — Thu, 14 Jul 2011 09:38:34 PDT

The ability to quantify soil–atmosphere water and energy exchange is important in understanding agricultural and natural ecosystems, as well as the earth's climate. We developed a one-dimensional vertical model that calculates solar radiation, canopy energy balance, surface energy balance, snowpack dynamics, soil water flow, and snow–soil–bedrock heat exchange, including soil water freezing. The processes are loosely coupled (solved sequentially) to limit the computational burden. The model was applied to describe water and energy dynamics for a northeast-facing mountain slope in the Dry Creek Experimental Watershed near Boise, ID. Calibration was achieved by optimizing the saturated soil hydraulic conductivity. Validation results showed that the model can successfully calculate seasonal dynamics in snow height, soil water content, and soil temperature. Both the calibration and validation years confirmed earlier results that evapotranspiration on the northeast-facing slope consumes approximately 60% of yearly precipitation, while deep percolation from the soil profile constitutes about 40% of yearly precipitation.

Simulated Soil Water Storage Effects on Streamflow Generation in a Mountainous Snowmelt Environment, Idaho, USA

M. S. Seyfried et al. — Fri, 08 Jul 2011 07:44:31 PDT

Although soil processes affect the timing and amount of streamflow generated from snowmelt, they are often overlooked in estimations of snowmelt-generated streamflow in the western USA. The use of a soil water balance modelling approach to incorporate the effects of soil processes, in particular soil water storage, on the timing and amount of snowmelt generated streamflow, was investigated. The study was conducted in the Reynolds Mountain East (RME) watershed, a 38 ha, snowmelt-dominated watershed in southwest Idaho. Snowmelt or rainfall inputs to the soil were determined using a well established snow accumulation and melt model (Isnobal). The soil water balance model was first evaluated at a point scale, using periodic soil water content measurements made over two years at 14 sites. In general, the simulated soil water profiles were in agreement with measurements (P < 0·05) as further indicated by high R² values (mostly > 0·85), y-intercept values near 0, slopes near 1 and low average differences between measured and modelled values. In addition, observed soil water dynamics were generally consistent with critical model assumptions. Spatially distributed simulations over the watershed for the same two years indicate that streamflow initiation and cessation are closely linked to the overall watershed soil water storage capacity, which acts as a threshold. When soil water storage was below the threshold, streamflow was insensitive to snowmelt inputs, but once the threshold was crossed, the streamflow response was very rapid. At these times there was a relatively high degree of spatial continuity of satiated soils within the watershed. Incorporation of soil water storage effects may improve estimation of the timing and amount of streamflow generated from mountainous watersheds dominated by snowmelt.

The Impact of a Shrinking Cryosphere on the Form of Arctic Alluvial Channels

James P. McNamara et al. — Fri, 08 Jul 2011 07:44:29 PDT

A shrinking cryosphere has important implications for the geomorphology of alluvial channels in arctic catchments because ice and permafrost alter the driving and resisting forces responsible for shaping local channel cross sections. For example, bedfast ice in shallow channels can suppress bedload transport during snowmelt events causing a reduction in the frequency of geomorphically effective flows, while cap ice in deeper channels can alter stream hydraulics. The impact of these local controls on catchment-scale geomorphic patterns, however, is less certain. A case study of one arctic river in northern Alaska, USA, reveals hydraulic geometry anomalies that can be reasonably explained by the influence of ice and permafrost on local channel form. However, there is a critical lack of data from rivers wholly contained in the continuous permafrost zone preventing comprehensive understanding of how arctic rivers might respond in a shrinking cryosphere.

Bedrock Infiltration and Mountain Block Recharge Accounting Using Chloride Mass Balance

Pam Aishlin et al. — Fri, 08 Jul 2011 07:44:27 PDT

Mountain front catchment net groundwater recharge (NR) represents the upper end of mountain block recharge (MBR) groundwater flow paths. Using environmental chloride in precipitation, streamflow and groundwater, we apply chloride mass balance (CMB) to estimate NR at multiple catchment scales within the 27 km² Dry Creek Experimental Watershed (DCEW) on the Boise Front, southwestern Idaho. The estimate for average annual precipitation partitioning to NR is approximately 14% for DCEW. In contrast, as much as 44% of annual precipitation routes to NR in ephemeral headwater catchments. NR in headwater catchments is likely routed to downgradient springs, baseflow, and MBR, while downgradient streamflow losses contribute further to MBR. A key assumption in the CMB approach is that the change in stored chloride during the study period is zero. We found that this assumption is violated in some individual years, but that a 5-year integration period is sufficient.

The Effect of Various Soil Hydraulic Property Estimates on Soil Moisture Simulations

Molly M. Gribb et al. — Wed, 22 Jun 2011 12:31:37 PDT

Models of water movement in unsaturated soils require accurate representations of the soil moisture retention and hydraulic conductivity curves; however, commonly used laboratory methods and pedotransfer functions (PTFs) are rarely verified against field conditions. In this study, we investigated the effects of using soil hydraulic property information obtained from different measurement and estimation techniques on one-dimensional model predictions of soil moisture content. Pairs of time domain reflectometry waveguides and tensiometers were installed at two depths in the side of a soil pit face to obtain in situ measurements. Undisturbed soil samples were taken near the instruments and subjected to particle size analysis, multistep outflow (MSO), and falling-head permeability tests to obtain estimates of the soil moisture retention curves. Three scaling methods were then applied to improve the fit of the various estimates to the field data. We found that soil hydraulic property estimates obtained from inverse methods lead to the best simulations of soil moisture dynamics, and that laboratory MSO tests or commonly used PTFs perform poorly. These laboratory and PTF estimates can be dramatically improved, however, by simply constraining the range of possible moisture contents to the minimum and maximum measured in the field. It appears that this method of scaling PTF results can be used to obtain soil hydraulic property inputs of sufficient accuracy for plot-scale modeling efforts without requiring expensive laboratory or in situ tests.

Modeling Runoff Generation in a Small Snow-Dominated Mountainous Catchment

T. J. Kelleners et al. — Mon, 13 Jun 2011 16:24:44 PDT

Snowmelt in mountainous areas is an important contributor to river water flows in the western United States. We developed a distributed model that calculates solar radiation, canopy energy balance, surface energy balance, snow pack dynamics, soil water flow, snow–soil–bedrock heat exchange, soil water freezing, and lateral surface and subsurface water flow. The model was applied to describe runoff generation in a subcatchment of the Dry Creek Experimental Watershed near Boise, ID. Calibration was achieved by optimizing the soil water field capacity (a trigger for lateral subsurface flow), lateral saturated soil hydraulic conductivity, and vertical saturated hydraulic conducitvity of the bedrock. Validation results show that the model can successfully calculate snow dynamics, soil water content, and soil temperature. Modeled streamflow for the validation period was underestimated by 53%. The timing of the streamflow was captured reasonably well (modeling efficiency was 0.48 for the validation period). The model calculations suggest that 50 to 53% of the yearly incoming precipitation in the subcatchment is consumed by evapotranspiration. The model results further suggest that 34 to 36% of the incoming precipitation is transformed into deep percolation into the bedrock, while only 11 to 16% is transformed into streamflow.

Improvement of Distributed Snowmelt Energy Balance Modeling with MODIS-Based NDSI-Derived Fractional Snow-Covered Area Data

Joel W. Homan et al. — Wed, 01 Jun 2011 10:54:18 PDT

Describing the spatial variability of heterogeneous snowpacks at a watershed or mountain-front scale is important for improvements in large-scale snowmelt modelling. Snowmelt depletion curves, which relate fractional decreases in snowcovered area (SCA) against normalized decreases in snow water equivalent (SWE), are a common approach to scale-up snowmelt models. Unfortunately, the kinds of ground-based observations that are used to develop depletion curves are expensive to gather and impractical for large areas. We describe an approach incorporating remotely sensed fractional SCA (FSCA) data with coinciding daily snowmelt SWE outputs during ablation to quantify the shape of a depletion curve. We joined melt estimates from the Utah Energy Balance Snow Accumulation and Melt Model (UEB) with FSCA data calculated from a normalized difference snow index snow algorithm using NASA’s moderate resolution imaging spectroradiometer (MODIS) visible (0Ð545–0Ð565 μm) and shortwave infrared (1•628–1•652 μm) reflectance data. We tested the approach at three 500 m2 study sites, one in central Idaho and the other two on the North Slope in the Alaskan arctic. The UEB-MODIS-derived depletion curves were evaluated against depletion curves derived from ground-based snow surveys. Comparisons showed strong agreement between the independent estimates.

Comparison of In-Channel Mobile–Immobile Zone Exchange During Instantaneous and Constant Rate Stream Tracer Additions: Implications for Design and Interpretation of Non-Conservative Tracer Experiments

Michael N. Gooseff et al. — Thu, 31 Mar 2011 15:44:38 PDT

The stream tracer experiment, including field tracer application and subsequent analysis of solute transport and storage, is an important tool in stream hydrology and ecology. However, there have been few comparisons of tracer dynamics between the commonly applied methods of instantaneous (IA) and constant rate (CRA) tracer additions. To determine whether there are fundamental differences between the two addition techniques due to surface storage zone loading and flushing during experiments, we compare longitudinal distributions of tracer dynamics of stream in-channel dead zones during IA and CRA experiments. Back-to-back IA and CRA additions were carried out in two morphologically distinct tundra stream reaches in Alaska. Dead zone tracer time series are determined by an aggregate of upstream transport and individual dead zone residence time distributions (RTDs). The dead zone breakthrough curves for both tracer addition techniques were not consistent, neither were aggregate RTDs observed in each dead zone. Flushing patterns of tracer from dead zones reveal that stream flushing after IA additions was slower than after CRA additions. However, whole-stream RTDs were similar between IA and CRA techniques in each reach. The implications of these findings are important to design and interpretation of IA and CRA stream tracer experiments, particularly those with reactive solutes whose transformations may depend on solute concentration. Thus, IA and CRA experiments may yield differing conclusions about non-conservative transport in streams because of the inherent differences in loading of transient storage zones between these two addition techniques, and potential differences in biogeochemical processing that may occur as a consequence.

Measuring Thaw Depth Beneath Peat-Lined Arctic Streams Using Ground-Penetrating Radar

John H. Bradford et al. — Tue, 29 Mar 2011 15:31:55 PDT

In arctic streams, depth of thaw beneath the stream channel is likely a significant parameter controlling hyporheic zone hydrology and biogeochemical cycling. As part of an interdisciplinary study of this system, we conducted a field investigation to test the effectiveness of imaging substream permafrost using ground-penetrating radar (GPR). We investigated three sites characterized by low-energy water flow, organic material lining the streambeds, and water depths ranging from 0·2 to 2 m. We acquired data using a 200 MHz pulsed radar system with the antennas mounted in the bottom of a small rubber boat that was pulled across the stream while triggering the radar at a constant rate. We achieved excellent results at all three sites, with a clear continuous image of the permafrost boundary both peripheral to and beneath the stream. Our results demonstrate that GPR can be an effective tool for measuring substream thaw depth.

Association of Ice and River Channel Morphology Determined Using Ground-Penetrating Radar in the Kuparuk River, Alaska

Heather Best et al. — Tue, 29 Mar 2011 15:31:53 PDT

We collected ground-penetrating radar data at 10 sites along the Kuparuk River and its main tributary, the Toolik River, to detect unfrozen water beneath river ice. We used 250 MHz and 500 MHz antennas to image both the ice-water interface and the river channel in late April 2001, when daily high temperatures were consistently below freezing and river ice had attained its maximum seasonal thickness. The presence of water below the river ice appears as a strong, horizontal reflection observed in the radar data and is confirmed by drill hole data. A downstream transition occurs from ice that is frozen to the bed, called bedfast ice, to ice that is floating on unfrozen water, called floating ice. This transition in ice type corresponds to a downstream change in channel size that was detected in previously conducted hydraulic geometry surveys of the Kuparuk River. We propose a conceptual model wherein the downstream transition from bedfast ice to floating ice is responsible for an observed step change in channel size due to enhanced bank erosion in large channels by floating ice.

Modeling the Spatially Varying Water Balance Processes in a Semi- Arid Mountainous Watershed of Idaho

Benjamin T. Stratton et al. — Tue, 29 Mar 2011 11:59:34 PDT

Mountainous watersheds in semi-arid regions are complex hydrologic systems. To critically evaluate the hydrological processes, high resolution spatio-temporal information is necessary. Also, calibrating and validating a watershed-scale model is necessary to enable our understanding of the water balance components in the gauged watersheds. The distributed Soil Water Assessment Tool (SWAT) hydrologic model was applied to a research watershed, the Dry Creek Experimental Watershed (DCEW), near Boise Idaho to investigate its water balance components both temporally and spatially. Daily streamflow data from four streamflow gauges were used for calibration and validation of the model. Monthly estimates of streamflow during the calibration phase by SWAT produced satisfactory results with a Nash Sutcliffe coefficient of model efficiency 0.79. Since it is a continuous simulation model, as opposed to an event-based model, it demonstrated the limited ability in capturing both streamflow and soil moisture for selected rain-on-snow events during the validation period between 2005-07. Our implementation of SWAT showed that seasonal and annual water balance partitioning of precipitation into evapotranspiration, streamflow, soil moisture and drainage was not only possible but closely followed the trends of a typical semiarid watershed in the intermountain west. This study highlights the necessity for better techniques to precisely identify and drive the model with commonly observed climatic inversion-related snowmelt or rain-on-snow weather events. Estimation of key parameters pertaining to soil (e.g., available water content and saturated hydraulic conductivity), snow (e.g., lapse rates, melting) and vegetation (e.g., leaf area index and maximum canopy index) using additional field observations in the watershed is critical for better prediction.

Profiles of Temporal Thaw Depths Beneath Two Arctic Stream Types Using Ground-Penetrating Radar

Troy R. Brosten et al. — Tue, 29 Mar 2011 11:59:32 PDT

Thaw depths beneath arctic streams may have significant impact on the seasonal development of hyporheic zone hydraulics. To investigate thaw progression over the 2004 summer season we acquired a series of ground-penetrating radar (GPR) profiles at five sites from May–September, using 100, 200 and 400 MHz antennas. We selected sites with the objective of including stream reaches that span a range of geomorphologic conditions on Alaska's North Slope. Thaw depths interpreted from GPR data were constrained by both recorded subsurface temperature profiles and by pressing a metal probe through the active layer to the point of refusal. We found that low-energy stream environments react much more slowly to seasonal solar input and maintain thaw thicknesses longer throughout the late season whereas thaw depths increase rapidly within high-energy streams at the beginning of the season and decrease over the late season period.

Estimating 3D Variation in Active-Layer Thickness Beneath Arctic Streams Using Ground-Penetrating Radar

Troy R. Brosten et al. — Tue, 29 Mar 2011 11:59:31 PDT

We acquired three-dimensional (3D) ground-penetrating radar (GPR) data across three stream sites on the North Slope, AK, in August 2005, to investigate the dependence of thaw depth on channel morphology. Data were migrated with mean velocities derived from multi-offset GPR profiles collected across a stream section within each of the 3D survey areas. GPR data interpretations from the alluvial-lined stream site illustrate greater thaw depths beneath riffle and gravel bar features relative to neighboring pool features. The peat-lined stream sites indicate the opposite; greater thaw depths beneath pools and shallower thaw beneath the connecting runs. Results provide detailed 3D geometry of active-layer thaw depths that can support hydrological studies seeking to quantify transport and biogeochemical processes that occur within the hyporheic zone.

Multi-Offset GPR Methods for Hyporheic Zone Investigations

T. R. Brosten et al. — Tue, 29 Mar 2011 11:59:28 PDT

Porosity of stream sediments has a direct effect on hyporheic exchange patterns and rates. Improved estimates of porosity heterogeneity will yield enhanced simulation of hyporheic exchange processes. Ground-penetrating radar (GPR) velocity measurements are strongly controlled by water content thus accurate measures of GPR velocity in saturated sediments provides estimates of porosity beneath stream channels using petrophysical relationships. Imaging the substream system using surface based reflection measurements is particularly challenging due to large velocity gradients that occur at the transition from open water to saturated sediments. The continuous multi-offset method improves the quality of subsurface images through stacking and provides measurements of vertical and lateral velocity distributions. We applied the continuous multi-offset method to stream sites on the North Slope, Alaska and the Sawtooth Mountains near Boise, Idaho, USA. From the continuous multi-offset data, we measure velocity using reflection tomography then estimate water content and porosity using the Topp equation. These values provide detailed measurements for improved stream channel hydraulic and thermal modelling.

Correction of Electronic Record for Weighing Bucket Precipitation Gauge Measurements

Anurag Nayak et al. — Tue, 08 Mar 2011 15:04:05 PST

Electronic sensors generate valuable streams of forcing and validation data for hydrologic models but are often subject to noise which must be removed as part of model input and testing database development. We developed an automated precipitation correction program (APCP) for weighing bucket precipitation gauge records, which are subject to several types of mechanical and electronic noise and discontinuities, including gauge maintenance, missing data, wind vibration, and sensor drift. Corrected cumulative water year precipitation from APCP did not exhibit an error bias and matched measured water year total precipitation within 2.1% for 58 station years tested. Removal of low-amplitude periodic noise was especially important for developing accurate instantaneous precipitation records at subdaily time steps. Model flexibility for use with other data types is demonstrated through application to time domain reflectometry soil moisture content data, which are also frequently subject to substantial noise.