Soil frost effects on soil water and runoff dynamics along a boreal transect: 2. Simulations

A physically based SVAT‐model was tested with soil and snow physical measurements, as well as runoff data from an 8600 m² catchment in northern Sweden in order to quantify the influence of soil frost on spring snowmelt runoff in a moderately sloped, boreal forest. The model was run as an array of connected profiles cascading to the brook. For three winter seasons (1995–98) it was able to predict the onset and total accumulation of the runoff with satisfactory accuracy. Surface runoff was identified as only a minor fraction of the total runoff occurring during short periods in connection with ice blocking of the water‐conducting pores. Little surface runoff, though, does not mean that soil frost is unimportant for spring runoff. Simulations without frost routines systematically underestimated the total accumulated runoff. The possibility of major frost effects appearing in response to specific combinations of weather conditions were also tested. Different scenarios of critical initial conditions for the winter, e.g. high water saturation and delayed snow accumulation leading to an increased frost penetration, were tested. These showed that under special circumstances there is potential for increased spring runoff due to soil frost. Copyright © 2001 John Wiley & Sons, Ltd.     

Snowmelt-generated runoff and soil erosion in Fife,Scotland

The hydrology and contrasting erosional responses of two snowmelt events on arable farmland in Fife, Scotland, are compared. Snowmelt-generated runoff in January 1993 caused widespread soil erosion across eastern Scotland. Gullying was exemplified by three sites in Fife, where thaw of a drifted snowpack was augmented by rainfall to produce a larger erosive response than meteorological data alone would have predicted. Up to 127 m³ of soil was lost from individual gullies in fields sown to winter cereals. In February 1996 snowfall of comparable depth again covered the field area, but a more uniform snowpack, slower thaw, greater crop cover and lower rainfall during the thaw phase combined to lessen the impact of erosion. These case studies demonstrate the complexity of the erosion/runoff relationship for rain on snow events, in which erosional severity depends not just on snow depth but on snow distribution, thaw rate and the amount and timing of rainfall during the thaw phase. © 1998 John Wiley & Sons, Ltd.     

Responses of soil moisture at different topographic positions to rainfall events along a steep subtropical forested hillslope

Understanding the dynamic response of soil moisture to rainfall is crucial for describing hydrological processes at the hillslope scale. However, because of sparse monitoring coupled with the complexity of water movement and steep topography, the findings of rainfall-related soil moisture dynamics have not always been consistent, indicating a need for systematic investigations of soil moisture dynamics and infiltration patterns following rainfall inputs at multiple topographic positions along a hillslope. This study aimed to examine the nature of these responses by characterizing and quantifying the response amplitude, rate and time for 37 large rainfall events at 25 combinations of topographic positions and soil depths along a steep forested hillslope. Our results showed that soil moisture responses under different rainfall patterns could be attributed to one or the other rainfall characteristics, such as rainfall intensity and amount. However, soil moisture dynamics at different hillslope positions after rainfall varied widely due to the controls of soil properties, topography, and non-equilibrium flow. Preferential flow was more evident under dry initial soil conditions than under wet initial soil conditions. Findings of this study reveal that the dynamic response patterns of soil moisture to rainfall do not always follow topographic controls, which can improve our understanding of water cycling related to the infiltration process at the hillslope scale, and support water resources management in subtropical mountain ecosystems.     

Soil moisture monitoring on a steep hillside

In this study, a soil monitoring system for a hillslope with steep relief and shallow soil depth was designed and installed to represent efficiently the spatial and temporal features of soil moisture. The study was conducted on a mountainous hillslope of the Sulmachun catchment (northeastern South Korea). The positions of soil moisture sensors were determined through a sequential procedure including intensive geomorphologic surveying of the study area, surface and subsurface terrain analysis, and inverse surveying. Using 26 sensors, soil moisture data from 11 locations were measured and recorded at hourly intervals over 380 h from 6 to 22 November 2003. Soil moisture response patterns were captured for a few consecutive rainfall events. The monitoring results are discussed in the context of soil moisture variations with terrain attributes. The immediate recharge and fast recession after a peak are the primary features of soil moisture in the upper zone. Stability and significant storage increase are distinct characteristics of soil moisture in the buffer zone and the flow path zone respectively. Spatial distribution of temporal soil moisture variations can be characterized in terms of recession, stability and recharge depending upon the topographic classification of a hillslope for this approach. Copyright © 2007 John Wiley & Sons, Ltd.     

Estimating source regions for snowmelt runoff in a Rocky Mountain basin: tests of a data‐based conceptual modeling approach

In many mountain basins, river discharge measurements are located far away from runoff source areas. This study tests whether a basic snowmelt runoff conceptual model can be used to estimate relative contributions of different elevation zones to basin‐scale discharge in the Cache la Poudre, a snowmelt‐dominated Rocky Mountain river. Model tests evaluate scenarios that vary model configuration, input variables, and parameter values to determine how these factors affect discharge simulation and the distribution of runoff generation with elevation. Results show that the model simulates basin discharge well (NSCE and R >0.90) when input precipitation and temperature are distributed with different lapse rates, with a rain‐snow threshold parameter between 0 and 3.3 °C, and with a melt rate parameter between 2 and 4 mm °C^?1 d^?1 because these variables and parameters can have compensating interactions with each other and with the runoff coefficient parameter. Only the hydrograph recession parameter can be uniquely defined with this model structure. These non‐unique model scenarios with different configurations, input variables, and parameter values all indicate that the majority of basin discharge comes from elevations above 2900 m, or less than 25% of the basin total area, with a steep increase in runoff generation above 2600 m. However, the simulations produce unrealistically low runoff ratios for elevations above 3000 m, highlighting the need for additional measurements of snow and discharge at under‐sampled elevations to evaluate the accuracy of simulated snow and runoff patterns. Copyright © 2013 John Wiley & Sons, Ltd.     

The role of overland flow and subsurface flow on the spatial distribution of soil moisture in the topsoil

Many investigations show relationships between topographical factors and the spatial distribution of soil moisture in catchments. However, few quantitative analyses have been carried out to elucidate the role of different hydrological processes in the spatial distribution of topsoil moisture in catchments. A spatially distributed rainfall—runoff model was used to investigate contributions of subsurface matric flow, macropore flow and surface runoff to the spatial distribution of soil moisture in a cultivated catchment. The model results show that lateral subsurface flow in the soil matrix or in macropores has a minor effect on the spatial distribution of soil moisture. Only when a perched groundwater table is maintained long enough, which is only possible if the subsurface is completely impermeable, may a spatial distribution in moisture content occur along the slope. Surface runoff, producing accumulations of soil moisture in flat flow paths of agricultural origin (field boundaries), was demonstrated to cause significant spatial variations in soil moisture within a short period after rainfall (<2 days). When significant amounts of surface runoff are produced, wetter moisture conditions will be generated at locations with larger upstream contributing areas. Copyright © 2001 John Wiley & Sons, Ltd.     

Hydrologic connectivity and threshold behavior of hillslopes with fragipans and soil pipe networks

Many concepts have been proposed to explain hydrologic connectivity of hillslopes with streams. Hydrologic connectivity is most often defined by qualitative assessment of spatial patterns in perched water tables or soil moisture on hillslopes without a direct linkage to water flow from hillslopes to streams. This form of hydrologic connectivity may not explain the hydrologic response of catchments that have network(s) of preferential flow paths, for example, soil pipes, which can provide intrinsic connectivity between hillslopes and streams. Duplex soils are known for developing perched water tables on hillslopes and fostering lateral flows, but the connectivity of localized perched water tables on hillslopes with soil pipes has not been fully established. The objectives of this study were to characterize pipeflow dynamics during storm events, the relationships between perched water tables on hillslopes and pipeflows, and their threshold behaviour. Two well‐characterized catchments in loess soil with a fragipan were selected for study because they contain multiple, laterally extensive (over 100 m) soil pipe networks. Hillslopes were instrumented with shallow wells adjacent to the soil pipes, and the wells and pipe collapse features were equipped with pressure transducers. Perched water tables developed on hillslopes during a wetting up period (October–December) and became well connected spatially across hillslope positions throughout the high flow period (January–March). The water table was not spatially connected on hillslopes during the drying out (April–June) and low flow (July–September) periods. Even when perched water tables were not well‐connected, water flowing through soil pipes provided hydrologic connectivity between upper hillslopes and catchment outlets. Correlations between soil pipeflow and perched water tables depended on the size and location of soil pipes. The threshold relationship between available soil‐moisture index plus storm precipitation and pipeflow was dependent on the season and strongest during dry periods and not high‐flow seasons. This study demonstrated that soil pipes serve as a catchment backbone of preferential flow paths that provide intrinsic connectivity between upper hillslopes and streams.     

Use of soil moisture data and curve number method for estimating runoff in the Loess Plateau of China

The Soil Conservation Service curve number (CN) method commonly uses three discrete levels of soil antecedent moisture condition (AMC), defined by the 5‐day antecedent rainfall depth, to describe soil moisture prior to a runoff event. However, this way may not adequately represent soil water conditions of fields and watersheds in the Loess Plateau of China. The objectives of this study were: (1) to determine the effective soil moisture depth to which the CN is most related; (2) to evaluate a discrete and a linear relationship between AMC and soil moisture; and (3) to develop an equation between CN and soil moisture to predict runoff better for the climatic and soil conditions of the Loess Plateau of China. The dataset consisted of 10 years of rainfall, runoff and soil moisture measurements from four experimental plots cropped with millet, pasture and potatoes. Results indicate that the standard CN method underestimated runoff depths for 85 of the 98 observed plot‐runoff events, with a model efficiency E of only 0·243. For our experimental conditions, the discrete and linear approaches improved runoff estimation, but still underestimated most runoff events, with E values of 0·428 and 0·445 respectively. Based on the measured CN values and soil moisture values in the top 15 cm of the soil, a non‐linear equation was developed that predicted runoff better with an E value of 0·779. This modified CN equation was the most appropriate for runoff prediction in the study area, but may need adjustments for local conditions in the Loess Plateau of China. Copyright © 2006 John Wiley & Sons, Ltd.     

Configuration of hydrological processes in a steep hillslope through causality analysis of soil moisture measurements

The vertical and lateral profiles of temporal variations in soil moisture are important for understanding the hydrological process along hillside transects. In this study, relationships among measured soil moistures were explored to configure the hydrological contributions of different flowpaths. All the measured soil moistures included a common stochastic structure because rainfall, the hydrometeological driver, is the main factor that determines the soil moisture response feature, and the infiltration process through the topsoil at a shallow depth is also common in all measured soil moisture histories. Therefore, the relationships between the measured series are also affected by both rainfall and topsoil infiltration. The common stochastic structure of the soil moisture series was removed via a prewhitening procedure. A systematic analysis procedure is presented to delineate the exclusive causal relationships among multiple soil moisture measurements. A monitoring system based on multiplexed time domain reflectometry was used to obtain soil moisture time series along two transects on a steep hillslope during the rainy season. The application of the proposed method for monitoring points in two adjacent locations provided 8, 12, 14, and 13, 16, 22 causal relationships for vertical, lateral in parallel, and diagonal directions, respectively, along the two transects. The point‐based contributions of the internal flowpath can be evaluated as the correlation is normalized in the context of inflow and outflow. The hydrological processes in the soil layer, vertical flow, lateral flow, downslope recharge, and return flow were quantified, and the relative importance of each hydrological component was determined to improve our understanding of the hydrological processes along the two transects of the study area. Copyright © 2011 John Wiley & Sons, Ltd.     

Lateral flow thresholds for aspen forested hillslopes on the Western Boreal Plain,Alberta, Canada

To predict the long‐term sustainability of water resources on the Boreal Plain region of northern Alberta, it is critical to understand when hillslopes generate runoff and connect with surface waters. The sub‐humid climate (P ≤ ET) and deep glacial sediments of this region result in large available soil storage capacity relative to moisture surpluses or deficits, leading to threshold‐dependent rainfall‐runoff relationships. Rainfall simulation experiments were conducted using large magnitude and high intensity applications to examine the thresholds in precipitation and soil moisture that are necessary to generate lateral flow from hillslope runoff plots representative of Luvisolic soils and an aspen canopy. Two adjacent plots (areas of 2·95 and 3·4 m²) of contrasting antecedent moisture conditions were examined; one had tree root uptake excluded for two months to increase soil moisture content, while the second plot allowed tree uptake over the growing season resulting in drier soils. Vertical flow as drainage and soil moisture storage dominated the water balances of both plots. Greater lateral flow occurred from the plot with higher antecedent moisture content. Results indicate that a minimum of 15–20 mm of rainfall is required to generate lateral flow, and only after the soils have been wetted to a depth of 0·75 m (C‐horizon). The depth and intensity of rainfall events that generated runoff > 1 mm have return periods of 25 years or greater and, when combined with the need for wet antecendent conditions, indicate that lateral flow generation on these hillslopes will occur infrequently. Copyright © 2008 John Wiley & Sons, Ltd.     

High-frequency monitoring of the occurrence of preferential flow on hillslopes and its relationship with rainfall features,soil moisture and landscape

ABSTRACT

High-frequency monitoring was conducted to quantify the frequency and controlling factors of preferential flow (PF) in a monsoon-influenced sub-humid mountainous catchment (6.48 km²) of Northern China. Rainfall was measured using nine bucket raingauges. Soil moisture probes were set up at 12 sites to observe the PF. Overall, 129 rainfall events were identified during the years 2014–2016. The average PF occurrence was 41%, which increased to 71% during heavy rainfall events (>20 mm) revealing a strong influence of the amount and intensity of rainfall. The study also revealed that the PF increased with antecedent soil moisture. Soil moisture was much higher on flat sites compared to sloping sites, providing evidence that the topography has a strong influence on rainfall infiltration and runoff which, subsequently, influence soil moisture variation and the occurrence of PF. Our findings provide valuable insights into the hydrological processes for studies in regions with similar environmental conditions.     

Soil temperature and soil moisture dynamics in winter and spring under heavy snowfall conditions in North-Eastern Japan

Warm winters and high precipitation in north-eastern Japan generate snow covers of more than three meters depth and densities of up to 0.55 g cm⁻³. Under these conditions, rain/snow ratio and snowmelt have increased significantly in the last decade under increasing warm winters. This study aims at understanding the effect of rain-on-snow and snowmelt on soil moisture under thick snow covers in mid-winter, taking into account that snowmelt in spring is an important source of water for forests and agriculture. The study combines three components of the Hydrosphere (precipitation, snow cover and soil moisture) in order to trace water mobility in winter, since soil temperatures remained positive in winter at nearly 0.3°C. The results showed that soil moisture increased after snowmelt and especially after rain-on-snow events in mid-winter 2018/2019. Rain-on-snow events were firstly buffered by fresh snow, increasing the snow water equivalent (SWE), followed by water soil infiltration once the water storage capacity of the snowpack was reached. The largest increase of soil moisture was 2.35 vol%. Early snowmelt increased soil moisture with rates between 0.02 and 0.035 vol% hr⁻¹ while, rain-on-snow events infiltrated snow and soil faster than snowmelt and resulted in rates of up to 1.06 vol% hr⁻¹. These results showed the strong connection of rain, snow and soil in winter and introduce possible hydrological scenarios in the forest ecosystems of the heavy snowfall regions of north-eastern Japan. Effects of rain-on-snow events and snowmelt on soil moisture were estimated for the period 2012–2018. Rain/snow ratio showed that only 30% of the total precipitation in the winter season 2011/2012 was rain events while it was 50% for the winter 2018/2019. Increasing climate warming and weakening of the Siberian winter monsoons will probably increase rain/snow ratio and the number of rain-on-snow events in the near future.     

Hillslope‐scale probabilistic characterization of soil moisture dynamics and average water balance

Probabilistic water balance modelling provides a useful framework for investigating the interactions between soil, vegetation, and the atmosphere. It has been used to estimate temporal soil moisture dynamics and ecohydrological responses at a point. This study combines a nonlinear rainfall–runoff theory with probabilistic water balance model to represent varied source area runoff as a function of rainfall depth and a runoff coefficient at hillslope scale. Analytical solutions of the soil‐moisture probability density function and average water balance model are then developed. Based on a sensitivity analysis of soil moisture dynamics, we show that when varied source area runoff is incorporated, mean soil moisture is always lower and total runoff higher, compared with the original probabilistic water balance model. The increased runoff from the inclusion of varied source area runoff is mainly because of a reduction in leakage when the index of dryness is less than one and evapotranspiration when the index of dryness is greater than one. Inclusion of varied source area runoff in the model means that the actual evapotranspiration is limited by less available water (i.e. water limit), which is stricter than Budyko’s and Milly’s water limit. Application of the model to a catchment located in Western Australia showed that the method can predict annual value of actual evapotranspiration and streamflow accurately. Copyright © 2012 John Wiley & Sons, Ltd.     

Land cover effects on soil infiltration capacity measured using plot scale rainfall simulation in steep tropical lowlands of Central Panama

Understanding land use/land cover (LULC) effects on tropical soil infiltration is crucial for maximizing watershed scale hydro-ecosystem services and informing land managers. This paper reports results from a multiyear investigation of LULC effects on soil bulk infiltration in steep, humid tropical, and lowland catchments. A rainfall simulator applied water at measured rates on 2 × 6 m plots producing infiltration through structured, granulated, and macroporous Ferralsols in Panama's central lowlands. Time-lapse electrical resistivity tomography (ERT) helped to visualize infiltration depth and bulk velocity. A space-for-time substitution methodology allowed a land-use history investigation by considering the following: (a) a continuously heavy-grazed cattle pasture, (b) a rotationally grazed traditional cattle pasture, (c) a 4-year-old (y.o.) silvopastoral system with nonnative improved pasture grasses and managed intensive rotational grazing, (d) a 7 y.o. teak (Tectona grandis) plantation, (e) an approximately 10 y.o. secondary succession forest, (f) a 12 y.o. coffee plantation (Coffea canephora), (g) an approximately 30 y.o. secondary succession forest, and (h) a >100 y.o. secondary succession forest. Within a land cover, unique plot sites totalled two at (a), (c), (d), (e), and (g); three at (b); and one at (f) and (h). Our observations confirmed measured infiltration scale dependency by comparing our 12 m² plot-scale measurements against 8.9 cm diameter core-scale measurements collected by others from nearby sites. Preferential flow pathways (PFPs) significantly increased soil infiltration capacity, particularly in forests greater than or equal to 10 y.o. Time-lapse ERT observations revealed shallower rapid bulk infiltration and increased rapid lateral subsurface flow in pasture land covers when compared with forest land covers and highlighted how much subsurface flow pathways can vary within the Ferralsol soil class. Results suggest that LULC effects on PFPs are the dominant mechanism by which LULC affects throughfall partitioning, runoff generation, and flow pathways.     

Estimation of mean water transit time on a steep hillslope in South Korea using soil moisture measurements and deuterium excess

A method for estimating daily mean transit time (DMTT) within a soil layer was proposed using field measurements of soil moisture. Vertical profiles of soil moisture time series were used for storage estimation. Water fluxes were evaluated through matrix and bypass flow. Variations in soil moisture and soil thickness were used to evaluate matrix flow. Exponential decay in depth of macropores was also used for bypass flow approximation. DMTT evaluation was compared to results obtained from a stable water isotope model using two years of data acquired on a steep granite hillslope in the Sulmachun watershed, South Korea. Various uncertainties in transit time evaluation such as model structure, non‐stationary assumption and data acquisition of existing approaches can be accounted for in the proposed methodology, and the flowpath contribution can be further configured in conjunction with hydrometric measurements. Probability density functions of isotope analyses were partially explained by transit time distributions that were based on soil moisture measurements. Supplementary sensitivity analyses for uncertainty configurations indicate that matrix flow is the primary process in determining transit time distribution while the impact of bypass flow is minor. The feasibility of a DMTT approach over isotope‐based methodologies highlights not only the strength of this proposed method, both in cost and time, but also its further application potential for existing soil moisture measurements. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

Soil moisture and rainfall intensity thresholds for runoff generation in southwestern Wisconsin agricultural watersheds

The goal of this study was to improve understanding of the factors that influence runoff generation during non‐frozen ground periods in small agricultural watersheds in southwestern Wisconsin where the landscapes are controlled by dolostone bedrock in order to provide agricultural producers with a manure management tool. Six small watersheds (ranging from 6 to 17 ha) within two southwestern Wisconsin farm sites (Discovery Farms Program (DFP) and Pioneer Farm (PF)) were instrumented, and surface runoff was continuously monitored from 2004 to 2007. The soils in all watersheds were formed in deep (~1 m) loessial sites. A direct‐plant management strategy and corn‐soybean crop rotation were utilized within watersheds at DFP. A conventional tillage system (chisel plow in the fall followed by soil finisher in the spring) and a corn‐oat‐alfalfa crop rotation were utilized within watersheds at PF. At PF, the amount of precipitation leaving the landscape as surface runoff (1.8%) was two times greater compared to DFP (0.9%), indicating that the direct‐plant management system was better at retaining precipitation than the chisel plow/soil finisher system. Using breakpoint regression analysis, a non‐linear response in runoff generation with antecedent soil moisture (ASM) was observed with a threshold ASM of 0.39 cm³cm^?3 (approximately 80% of total porosity) for all six watersheds. Below this threshold, runoff coefficients were near zero. Above this threshold, runoff coefficients increased with ASM. A non‐linear response in runoff generation with maximum 30 min rainfall intensity (I30) was also observed, and threshold I30 values increased as ASM decreased and as crop cover increased. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of soil moisture dynamics beneath olive trees

Olive cultivation is a widespread land use in Mediterranean climates. The proper implementation of soil and water conservation practices in groves requires detailed knowledge of the governing hydrological processes. In this work topsoil moisture dynamics under wet and dry conditions and across a small catchment was investigated in the inter row (IR) and directly under the olive tree canopies (UC). We do this using a sensor network (11 stations) and a simple bucket model which was calibrated (June, 2011–2012) and validated (June, 2012–2013). During most of the year the normalized soil moisture contents (s) were greater in the IR than under UC, with an average normalized soil moisture difference of 0.12. The difference between UC and IR normalized soil moisture followed a seasonal pattern, reaching a maximum near 0.30 during spring. An analysis of the normalized soil moisture probability density functions (pdfs) was bimodal, showing characteristic dominant wet and dry soil moisture states, with the highest probability densities for the dry state. Overall the spatial variability of soil moisture was lower UC than in the IR. This was a result of the soil moisture buffering capacity of the canopy with respect to rainfall and evaporation, in addition to observed differences in soil properties. Hourly soil moisture data were successfully modelled (R² > 0.85), both UC and in the IR, yet with the inclusion of a simple formulation for canopy interception for the former. The results provide insight into how olive trees change hydrological processes in their neighbourhood. Copyright © 2016 John Wiley & Sons, Ltd.     

Soil moisture response to snowmelt timing in mixed‐conifer subalpine forests

Western US forest ecosystems and downstream water supplies are reliant on seasonal snowmelt. Complex feedbacks govern forest–snow interactions in which forests influence the distribution of snow and the timing of snowmelt but are also sensitive to snow water availability. Notwithstanding, few studies have investigated the influence of forest structure on snow distribution, snowmelt and soil moisture response. Using a multi‐year record from co‐located observations of snow depth and soil moisture, we evaluated the influence of forest‐canopy position on snow accumulation and snow depth depletion, and associated controls on the timing of soil moisture response at Boulder Creek, Colorado, Jemez River Basin, New Mexico, and the Wolverton Basin, California. Forest‐canopy controls on snow accumulation led to 12–42 cm greater peak snow depths in open versus under‐canopy positions. Differences in accumulation and melt across sites resulted in earlier snow disappearance in open positions at Jemez and earlier snow disappearance in under‐canopy positions at Boulder and Wolverton sites. Irrespective of net snow accumulation, we found that peak annual soil moisture was nearly synchronous with the date of snow disappearance at all sites with an average deviation of 12, 3 and 22 days at Jemez, Boulder and Wolverton sites, respectively. Interestingly, sites in the Sierra Nevada showed peak soil moisture prior to snow disappearance at both our intensive study site and the nearby snow telemetry stations. Our results imply that the duration of soil water stress may increase as regional warming or forest disturbance lead to earlier snow disappearance and soil moisture recession in subalpine forests. Copyright © 2014 John Wiley & Sons, Ltd.