A TEST OF THE RADIATIVE ENERGY BALANCE OF THE SHAW MODEL FOR SNOWCOVER

Snow and ice present interesting challenges to hydrologists. Simulating the radiative balance over snow, which is an important part of surface–atmosphere interactions, is particularly challenging because of the decay in albedo over time and the difficulty in estimating surface temperature and incoming long-wave radiation fluxes. Few models are available that include a comprehensive energy and water balance for cold season conditions. The simultaneous heat and water model (SHAW) is a detailed, physical process model of a vertical, one-dimensional canopy–snow–residue–soil system which integrates the detailed physics of heat and water transfer through a plant canopy, snow, residue and soil into one simultaneous solution. Detailed provisions for metamorphosis of the snowpack are included. The SHAW model was applied to data for one winter/spring season (November to May) on a ploughed field in Minnesota without prior calibration to test the performance of the radiation components. Maximum snow depth during this period was 30 cm. For the nearly 100 days of snowcover, the model accounted for 69% of the variation in net solar radiation, 66% of the variation in incoming long-wave radiation, 87% of the variation in emitted long-wave radiation, 26% of the variation in net long-wave radiation and 55% of the variation in net radiation balance. Mean absolute error in simulated values ranged from 10 W m⁻² for emitted long-wave radiation to 27 W m⁻² for the entire net radiation balance. Mean bias error ranged from 8 W m⁻² for emitted long-wave radiation to −16 W m⁻² for the entire net radiation balance. When the entire 170 days of simulation, which included periods without snowcover, were included in the analysis, the variation in observed values increased greatly. As a result, the variation in observed values accounted for by the model increased to 97, 71, 93, 56 and 94%, respectively, while the mean absolute and mean bias errors in simulated values remained nearly the same. Model modifications and parameter adjustments necessary to improve winter-time simulation were investigated. Simulation results suggest that the SHAW model may be a useful tool in simulating the interactive influences of radiative transfer at the surface–atmosphere interface.     

The effects of forest litter on snow energy budget in the Tianshan Mountains,China

Forest litter exerts an impact on the energy budget of snow surfaces, which lie beneath forest canopies. In this study, we measured shortwave and longwave radiation levels, as well as quantities of Asian spruce (Picea schrenkinan ) forest litter, over 3 snow study plots that representing an open environment, 20% forest canopy openness (20% FCO), and 80% forest canopy openness (80% FCO). The fractional litter coverage (lc ) was obtained through the binarization of digital photographs of forest litter. The effects of forest litter on snow surface albedo (α ), snow surface temperature (T _s ), upward shortwave and longwave radiation (K _↑ and L _↑), and sensible heat flux (H ) were then analyzed. According to our results, the energy budget over snow surface influenced by forest litter principally due to forest litter forcing α decrease and T _s increase. The effects of forest litter on the energy budget increased with time and lc . We found that forest litter exerted the most significant impact on K _↑ and L _↑ at daytime during the latter stages of the snowmelt period. The influence of forest litter on H was more apparent on windy days. The presence of forest litter increased gains in shortwave radiation and losses in longwave radiation and decreased gains in H . Compared to the simulated energy (K _↑ + L _↑ + H ) over a snow surface without litter, the calculated energy decreased by ?13.4 W/m² and increased by 9.0 W/m², respectively, at the 20% FCO and 80% FCO sites during the latter stages of the snowmelt period. Overall, forest litter facilitated snow surface energy gains at the 80% FCO site and impeded them at the 20% FCO site during the latter stages of the snowmelt period.     

Influences of forest canopy on snowpack accumulation and isotope ratios

The stable water isotopes, ²H and ¹⁸O, can be useful environmental tracers for quantifying snow contributions to streams and aquifers, but characterizing the isotopic signatures of bulk snowpacks is challenging because they can be highly variable across the catchment landscape. In this study, we investigate one major source of isotopic heterogeneity in snowpacks: the influence of canopy cover. We measured amounts and isotopic compositions of bulk snowpack, throughfall, and open precipitation during seven campaigns in mid-winter 2018 along forest-grassland transects at three different elevations (1196, 1297, and 1434 m above sea level) in a pre-Alpine catchment in Switzerland. Snowpack storages under forest canopies were 67 to 93% less than in adjacent open grasslands. On average, the water isotope ratios were higher in the snowpacks under forest canopy than in open grasslands (by 13.4 ‰ in δ²H and 2.3 ‰ in δ¹⁸O). This isotopic difference mirrored the higher isotope values in throughfall compared with open snowfall (by 13.5 ‰ in δ²H and 2.2 ‰ in δ¹⁸O). Although this may suggest that most of the isotopic differences in snowpacks under forests versus in open grasslands were attributable to canopy interception effects, the temporal evolution of snowpack isotope ratios indicated preferential effluxes of lighter isotopes as energy inputs increased and the snowpack ripened and melted. Understanding these effects of forest canopy on bulk snowpack snow water equivalent and isotopic composition are useful when using isotopes to infer snowmelt processes in landscapes with varying forest cover.     

Hydrological and hydrochemical response of a small canadian shield catchment to late winter rain-on-snow events

A study was undertaken during the winter of 1990–1991 in a small (3.7 ha) Canadian Shield catchment to examine the hydrological and hydrochemical response during rain-on-snow events. The results are presented of two large (37.9 and 34.6 mm) rain-on-snow events occurring in early and late March 1991. Peak and total runoff and the groundwater response from the two events are significantly different. Hydrological data indicate that these differences can be attributed to a combination of meteorological (temperature) and physical conditions (antecedent snowpack ripeness, soil moisture and groundwater levels). An immature snowpack (low temperature and density) combined with low antecedent soil moisture conditions significantly reduced the magnitude of the net hydrological input and runoff from the catchment during the early March event, whereas a more mature snowpack and high antecedent soil moisture conditions led to a large runoff event during late March. During both rain-on-snow events a significant portion of the pre-event snowpack chemical load was lost. Based on the maximum snowpack chemical load measured before the events, the two large rain-on-snow events and a brief mid-March warm period during which there were two much smaller rain-on-snow events removed 78% of the hydrogen ion and 63% of the sulphate and nitrate load from the snowpack, while only reducing snowpack water equivalence by 7%. A two-component (rain and snowmelt) isotopic (δ¹⁸O SMOW %₀) separation of snowmelt lysimeter water during the two events indicated that snowmelt was an important (50 and 65%, respectively) water source available for infiltration and runoff at the snow-soil interface. Considering the high hydrogen ion loadings to the catchment during these two events (3.3 and 3.0 mequiv.m^?2, respectively) streamflow pH was not significantly reduced due to an increase in the discharge of well-buffered groundwater. A two-component isotopic hydrograph separation of peak stream discharge during the 2–3 March event indicated that 75% of the total flow was groundwater. In mid-latitude acid-sensitive catchments, winter rain-on-snow events are an important hydrological occurrence due to their ability to elute much of the chemical load (H⁺, SO₄, NO₃) from the snowpack before the onset of spring melt when the maximum annual hydrological input typically occurs.     

Vegetation water sources in California's Sierra Nevada (USA) are young and change over time,a multi-isotope (δ18O, δ2H, 3H) tracer approach.

Sierra Nevada forests transpire a significant amount of California's water resources, sparking interest in applying forest management to improve California's water supply. Determining the source water of evapotranspiration enables forest managers to make informed decisions. To this end, a significant interest in critical zone science is to develop new methods to work across time scales to predict subsurface water storage and use. In this study, forest vegetation accessed young water and switched sources depending on availability, suggesting that forest drought vulnerability may depend on the range of water sources available (rain, snowmelt and deeply stored water). This finding also suggests that changes in transpiration rates may have immediate effects on water sources in close proximity to vegetation, and delayed effects on storage and runoff. New δ¹⁸O, δ²H and ³H data were used to track precipitation, runoff, evapotranspiration and storage through the critical zone seasonally, including seasons where evapotranspiration and snowmelt were in phase (winter snowmelt) and out of phase (seasonally dry summer). The main source of this headwater catchment's runoff is derived from its meadow saturated zone water, which was dominated by snowmelt. Water that originated as snowmelt contributed to transpiration, unless other sources, such as recent rain, became available. In cases where xylem δ¹⁸O and δ²H signatures matched those of deeper saturated zone water, ³H data showed that xylem water was distinctly younger than the deep saturated zone water. During 2016, which experienced relatively normal snowpack in winter and seasonally dry summer conditions, mean summer saturated zone water and vegetation water were similar in δ¹⁸O, −12.4 ± 0.04 ‰ and − 12.5 ± 0.3 ‰, respectively, but were distinctly different in ³H, 5.5 ± 0.2 pCi/L and 13.7 ± 1.1 pCi/L, respectively. While δ¹⁸O shows that vegetation and meadow saturated zone water have similar origins, ³H shows they have dissimilar ages.     

An improved hemispherical photography model for stream surface shortwave radiation estimations in a central U.S. hardwood forest

Hemispherical photographs of forest canopies can be used to develop sophisticated models that predict incident below canopy shortwave radiation on the surface of interest (i.e. soil and water). Hemispherical photographs were collected on eight dates over the course of a growing season to estimate leaf area index and to quantify solar radiation incident on the surface of two stream reaches based on output from Gap Light Analyser and Hemisfer software. Stream reaches were shaded by a mixed‐deciduous Ozark border forested riparian canopy. Hemispherical photo model results were compared to observed solar radiation sensed at climate stations adjacent to each stream reach for the entire 2010 water year. Modeled stream‐incident shortwave radiation was validated with above‐stream pyranometers for the month of September. On average, the best hemispherical photo models underestimated daily averages of solar radiation by approximately 14% and 12% for E–W and N–S flowing stream reaches, respectively (44.7 W/m² measured vs 38.4 W/m² modeled E–W, 46.8 W/m² vs. 41.1 W/m²N–S). The best hemispherical photo models overestimated solar radiation relative to in–Stream pyranometers placed in the center of each stream reach by approximately 7% and 17% for E–W and N–S stream reaches respectively (31.3 W/m² measured vs 33.5 W/m² modeled E–W, 31.5 W/m² vs. 37.1 W/m²N–S). The model provides a geographically transferable means for quantifying changes in the solar radiation regime at a stream surface due to changes in canopy density through a growing season, thus providing a relatively simple method for estimating surface and water heating in canopy altered environments (e.g. forest harvest). Copyright © 2012 John Wiley & Sons, Ltd.     

Energy fluxes and melt rate of a seasonal snow cover in the Moroccan High Atlas

ABSTRACT

In this study, we used an energy balance model and two simple methods based on readily available data to identify the processes driving the point-scale energy and mass balance of the snowpack. Data were provided from an experimental site located at 3200 m. All models were evaluated by comparing observed and modelled snow water equivalents. Performances are variable from one season to the next and the energy balance model gives better results (mean of root mean square error, RMSE = 25 mm and r² = 0.90) than the two simplified approaches (mean of RMSE = 54 mm and r² = 0.70). There are significant amounts of snow sublimation but they are highly variable from season to season, depending on wind conditions (between 7 and 20% of the total). While the main source of energy for melting is net radiation, the amount of heat brought by sensible heat flux is significant for two of the most windy snow seasons.

Editor Z.W. Kundzewicz Associate editor not assigned     

Development and application of a spatially‐distributed Arctic hydrological and thermal process model (ARHYTHM)

A process‐based, spatially distributed hydrological model was developed to quantitatively simulate the energy and mass transfer processes and their interactions within arctic regions (arctic hydrological and thermal model, ARHYTHM). The model first determines the flow direction in each element, the channel drainage network and the drainage area based upon the digital elevation data. Then it simulates various physical processes: including snow ablation, subsurface flow, overland flow and channel flow routing, soil thawing and evapotranspiration. The kinematic wave method is used for conducting overland flow and channel flow routing. The subsurface flow is simulated using the Darcian approach. The energy balance scheme was the primary approach used in energy‐related process simulations (snowmelt and evapotranspiration), although there are options to model snowmelt by the degree‐day method and evapotranspiration by the Priestley–Taylor equation. This hydrological model simulates the dynamic interactions of each of these processes and can predict spatially distributed snowmelt, soil moisture and evapotranspiration over a watershed at each time step as well as discharge in any specified channel(s). The model was applied to Imnavait watershed (about 2·2 km²) and the Upper Kuparuk River basin (about 146 km²) in northern Alaska. Simulated results of spatially distributed soil moisture content, discharge at gauging stations, snowpack ablations curves and other results yield reasonable agreement, both spatially and temporally, with available data sets such as SAR imagery‐generated soil moisture data and field measurements of snowpack ablation, and discharge data at selected points. The initial timing of simulated discharge does not compare well with the measured data during snowmelt periods mainly because the effect of snow damming on runoff was not considered in the model. Results from the application of this model demonstrate that spatially distributed models have the potential for improving our understanding of hydrology for certain settings. Finally, a critical component that led to the performance of this modelling is the coupling of the mass and energy processes. Copyright © 2000 John Wiley & Sons, Ltd.     

Charred forests accelerate snow albedo decay: parameterizing the post‐fire radiative forcing on snow for three years following fire

As large, high‐severity forest fires increase and snowpacks become more vulnerable to climate change across the western USA, it is important to understand post‐fire disturbance impacts on snow hydrology. Here, we examine, quantify, parameterize, model, and assess the post‐fire radiative forcing effects on snow to improve hydrologic modelling of snow‐dominated watersheds having experienced severe forest fires. Following a 2011 high‐severity forest fire in the Oregon Cascades, we measured snow albedo, monitored snow, and micrometeorological conditions, sampled snow surface debris, and modelled snowpack energy and mass balance in adjacent burned forest (BF) and unburned forest sites. For three winters following the fire, charred debris in the BF reduced snow albedo, accelerated snow albedo decay, and increased snowmelt rates thereby advancing the date of snow disappearance compared with the unburned forest. We demonstrate a new parameterization of post‐fire snow albedo as a function of days‐since‐snowfall and net snowpack energy balance using an empirically based exponential decay function. Incorporating our new post‐fire snow albedo decay parameterization in a spatially distributed energy and mass balance snow model, we show significantly improved predictions of snow cover duration and spatial variability of snow water equivalent across the BF, particularly during the late snowmelt period. Field measurements, snow model results, and remote sensing data demonstrate that charred forests increase the radiative forcing to snow and advance the timing of snow disappearance for several years following fire. Copyright © 2016 John Wiley & Sons, Ltd.     

A detailed energy budget analysis of river supercooling and the importance of accurately quantifying net radiation to predict ice formation

River supercooling and ice formation is a regular occurrence throughout the winter in northern countries. The resulting frazil ice production can obstruct the flow through intakes along the river, causing major problems for hydropower and water treatment facilities, among others. Therefore, river ice modellers attempt to calculate the river energy budget and predict when supercooling will occur in order to anticipate and mitigate the effects of potential intake blockages. Despite this, very few energy budget studies have taken place during freeze-up, and none have specifically analysed individual supercooling events. To improve our understanding of the freeze-up energy budget detailed measurements of air temperature, relative humidity, barometric pressure, wind speed and direction, short- and longwave radiation, and water temperature were made on the Dauphin River in Manitoba. During the river freeze-up period of late October to early November 2019, a total of six supercooling events were recorded. Analysis of the energy budget throughout the supercooling period revealed that the most significant heat source was net shortwave radiation, reaching up to 298 W/m², while the most significant heat loss was net longwave radiation, accounting for losses of up to 135 W/m². Longwave radiation was also the most significant heat flux overall during the individual supercooling events, accounting for up to 84% of the total heat flux irrespective of flux direction, highlighting the importance of properly quantifying this flux during energy budget calculations. Five different sensible (Q_h) and latent (Q_e) heat flux calculations were also compared, using the bulk aerodynamic method as the baseline. It was found that the Priestley and Taylor method most-closely matched the bulk aerodynamic method on a daily timescale with an average offset of 8.5 W/m² for Q_h and 10.1 W/m² for Q_e, while a Dalton-type equation provided by Webb and Zhang was the most similar on a sub-daily timescale with average offsets of 20.0 and 14.7 W/m² for Q_h and Q_e, respectively.     

Variations in snowmelt energy and energy balance characteristics with larch forest density on Mt Iwate,Japan: observations and energy balance analyses

The variation in snowmelt energy and energy components were evaluated with respect to forest density. Surface snowmelt rates, surface evaporation from snow cover and meteorological elements were measured in the open and under sparse (411 trees/ha) and dense (1433 trees/ha) larch canopies. The surface snowmelt rate decreased as the forest density increased. Based on the observations and energy balance analyses, we concluded the following. (1) Albedo decreased while the bulk coefficient for latent heat increased with forest density. (2) The duration of snowmelt increased with forest density because the energy for nocturnal cooling of the snow cover decreased. (3) When comparing the open and forested sites, the changes in snowmelt energy with forest density were caused by sensible heat flux. However, the contribution of net radiation was highest in the forested sites. Therefore, the effects of forest cover on the snowmelt energy were different when comparing both the open and forested sites and the sparse and densely forested sites. (4) The ratio of net radiation to snowmelt energy increased with forest density; although both snowmelt energy and net radiation decreased with increased forest density, the snowmelt energy decreased more rapidly. Sensible heat also decreased as forest density increased. Both albedo and downward long‐wave radiation influenced net radiation. Copyright © 1999 John Wiley & Sons, Ltd.     

Slopewash,surface runoff and fine-litter transport in forest and landslide scars in humid-tropical steeplands,luquillo experimental forest,Puerto Rico

Rainfall, slopewash (the erosion of soil particles), surface runoff and fine-litter transport steepland sites in the Luquillo Experimental Forest, Puerto Rico (18° 20’ N, 65° 45’ W) were measured from 1991 to 1995. Hillslopes underlain by (1) Cretaceous tuffaceous sandstone and silstone in subtropical rain (tanonuco) forest with vegetation recovering from Hurricane Hugo (1989), and (2) Tertiary quartz diorite in subtropical lower mantone wet (colorado and dwarf) forest with undisturbed forest canopy were compared to recent landslide scars. Monthly surface runoff on these very steep hillslopes (24° to 43°) was only 0·2 to 0·5 per cent of monthly rainfall. Slopewash was higher in sandy loam soils whose parent material is quartz diorite (averaging 46 g m⁻² a⁻¹) than in silty clay loam soils derived from tuffaceous sandstone and siltstone where the average was 9 g m⁻² a⁻¹. Annual slopewash of 100 to 349 g m⁻² on the surfaces of two recent, small landslide scars was measured initially but slopewash decreased to only 3 to 4 g m⁻² a⁻¹ by the end of the study. The mean annual mass of fine litter (mainly leaves and twigs) transported downslope at the forested sites ranged from 5 to 8 g m⁻² and was lower at the tabonuco forest site, where post-Hurricane Hugo recovery is still in progress. Mean annual fine-litter transport was 2·5 g m⁻² on the two landslide scars. Copyright © 1999 John Wiley & Sons, Ltd.     

Seasonal snowpack dynamics and runoff in a cool temperate forest: lysimeter experiment in Niigata,Japan

Seasonal snowpack dynamics are described through field measurements under contrasting canopy conditions for a mountainous catchment in the Japan Sea region. Microclimatic data, snow accumulation, albedo and lysimeter runoff are given through the complete winter season 2002–03 in (1) a mature cedar stand, (2) a larch stand, and (3) a regenerating cedar stand or opening. The accumulation and melt of seasonal snowpack strongly influences streamflow runoff during December to May, including winter baseflow, mid‐winter melt, rain on snow, and diurnal peaks driven by radiation melt in spring. Lysimeter runoff at all sites is characterized by constant ground melt of 0·8–1·0 mm day⁻¹. Rapid response to mid‐winter melt or rainfall shows that the snowpack remains in a ripe or near‐ripe condition throughout the snow‐cover season. Hourly and daily lysimeter discharge was greatest during rain on snow (e.g. 7 mm h⁻¹ and 53 mm day⁻¹ on 17 December) with the majority of runoff due to rainfall passing through the snowpack as opposed to snowmelt. For both rain‐on‐snow and radiation melt events lysimeter discharge was generally greatest at the open site, although there were exceptions such as during interception melt events. During radiation melt instantaneous discharge was up to 4·0 times greater in the opening compared with the mature cedar, and 48 h discharge was up to 2·5 times greater. Perhaps characteristic of maritime climates, forest interception melt is shown to be important in addition to sublimation in reducing snow accumulation beneath dense canopies. While sublimation represents a loss from the catchment water balance, interception melt percolates through the snowpack and contributes to soil moisture during the winter season. Strong differences in microclimate and snowpack albedo persisted between cedar, larch and open sites, and it is suggested further work is needed to account for this in hydrological simulation models. Copyright © 2005 John Wiley & Sons, Ltd.     

WINTER RADIATION EXTINCTION AND REFLECTION IN A BOREAL PINE CANOPY: MEASUREMENTS AND MODELLING

Predicting the rate of snowmelt and intercepted snow sublimation in boreal forests requires an understanding of the effects of snow-covered conifers on the exchange of radiant energy. This study examined the amount of intercepted snow on a jack pine canopy in the boreal forest of central Saskatchewan and the shortwave and net radiation exchange with this canopy, to determine the effect of intercepted snow and canopy structure on shortwave radiation reflection and extinction and net radiation attenuation in a boreal forest. The study focused on clear sky conditions, which are common during winter in the continental boreal forest. Intercepted snow was found to have no influence on the clear-sky albedo of the canopy, the extinction of short wave radiation by the canopy or ratio of net radiation at the canopy top to that at the surface snow cover. Because of the low albedo of the snow-covered canopy, net radiation at the canopy top remains positive and a large potential source of energy for sublimation. The canopy albedo declines somewhat as the extinction efficiency of the underlying canopy increases. The extinction efficiency of short wave radiation in the canopy depends on solar angle because of the approximately horizontal orientation of pine branches. For low solar angles above the horizon, the extinction efficiency is quite low and short wave transmissivity through the canopy is relatively high. As the solar angle increases, extinction increases up to angles of about 50°, and then declines. Extinction of short wave radiation in the canopy strongly influences the attenuation of net radiation by the canopy. Short wave radiation that is extinguished by branches is radiated as long wave, partly downwards to the snow cover. The ratio of net radiation at the canopy top to that at the snow cover surface increases with the extinction of short wave radiation and is negative for low extinction efficiencies. For the pine canopy examined, the daily mean net radiation at the snow cover surface became positive when daily mean solar angles exceeded 22° in late March. Hence, canopy structure and solar angle control the net radiation at the snow cover surface during clear sky conditions and will govern the timing and rate of snowmelt. Models of intercepted snow sublimation and forest snowmelt could beneficially incorporate the canopy radiation balance, which can be extrapolated to stands of various canopy densities, coverage and heights in a physically based manner. Such models could hence avoid ‘empirical’ temperature index measures that cannot be extrapolated with confidence.     

Modelling snowpack bulk density using snow depth,cumulative degree-days,and climatological predictor variables

Snowpack water equivalent (SWE) is a key variable for water resource management in snow-dominated catchments. While it is not feasible to quantify SWE at the catchment scale using either field surveys or remotely sensed data, technologies such as airborne LiDAR (light detection and ranging) support the mapping of snow depth at scales relevant to operational water management. To convert snow depth to water equivalent, models have been developed to predict SWE or snowpack density based on snow depth and additional predictor variables. This study builds upon previous models that relate snowpack density to snow depth by including additional predictor variables to account for (1) long-term climatologies that describe the prevailing conditions influencing regional snowpack properties, and (2) the effect of intra- and inter-year variability in meteorological conditions on densification through a cumulative degree-day index derived from North American Regional Reanalysis products. A non-linear model was fit to 114 506 snow survey measurements spanning 41 years from 1166 snow courses across western North America. Under spatial cross-validation, the predicted densities had a root-mean-square error of 47.1 kg m⁻³, a mean bias of −0.039 kg m⁻³, and a Nash-Sutcliffe Efficiency of 0.70. The model developed in this study had similar overall performance compared to a similar regression-based model reported in the literature, but had reduced seasonal biases. When applied to predict SWE from simulated depths with random errors consistent with those obtained from LiDAR or Structure-from-Motion, 50% of the SWE estimates for April and May fell within −45 to 49 mm of the observed SWE, representing prediction errors of −15% to 20%.     

One‐year‐long runoff forecast by a single snowpack evaluation

A degree‐day‐based model is presented for a 1 year ahead runoff forecast, with 1 day time steps. The input information is a single snowpack evaluation collected at the beginning of the snowmelt season. The snow‐cover dynamics, the key information for long‐term snowmelt forecast, are described by the snow‐line dynamics, i.e. by the movements of the downhill snowpack limit. The snowmelt volume, estimated by the snow‐line dynamics, is the exogenous input of an autoregressive transformation model. The model is calibrated by a least‐squares procedure on the basis of observed daily runoff data and the corresponding measurements of the snowpack volume (one measurement per year). A real‐world case study on the Alto Tunuyan River basin (2380 km², Argentinean Andes) is presented. The 1 year ahead Alto Tunuyan River runoff patterns, computed for both calibration and validation periods, reveal high agreement with observed streamflows. Copyright © 2004 John Wiley & Sons, Ltd.     

Inertial oscillations in the Irish Sea

An inertial current with a maximum speed of about 0.20 m s¹ was observed in June 1981, by a current meter 17 m below the surface of the thermally stratified waters of the western Irish Sea. The dynamics of the surface waters have been modelled using contemporary wind data recorded at the site to drive a horizontally unbounded surface slab 20 m deep with a velocity decay time of about three inertial periods (1.86 days). From thermistor chain data, inertial frequency potential energy was observed to radiate beneath the thermocline, descending at a rate of about 0.06 cm s¹. The energy available for internal mixing from the observed oscillation may have been as much as 80 J m², compared with 150 J m² from wind-induced surface mixing. The mean surface mixing power in summer is about7.1 × 10⁴W m², while the mean shear-induced mixing power could be as much as2.0 × 10⁴W m². The likely occurrence of strong inertial events has been estimated using the wind record fcr 10 successive summers at Dublin airport. Large oscillations with speeds 0.20m s¹, or mixed layer kinetic energy(MLKE) > 400 J m², have an average duration of about 28 h. They occur about once every 2 months between mid-May and mid-August but seem to be more frequent from then until mid-September. Moderate oscillations (MLKE > 200 J m²) are more evenly distributed, occurring about once every 2 weeks. Small oscillations (MLKE > 40 J m²) occur nearly 50% of the time.     

ROAD SALT ACCUMULATION IN HIGHWAY SNOW BANKS AND TRANSPORT THROUGH THE UNSATURATED ZONE OF THE OAK RIDGES MORAINE,SOUTHERN ONTARIO

Pathways and fate of road deicing salt (NaCl) applied during the 1994–1995 winter were studied for a 14-km section of a major highway that crosses the Oak Ridges Moraine in southern Ontario. Total salt applications over the winter ranged from 29 to 74 kg m⁻¹ of highway, and NaCl concentrations in snow banks adjacent to the roadway reached 9400 mg l⁻¹ during the later stages of snow cover development. This salt was released to the ground surface during snowmelt. Sodium chloride (NaCl) loadings to soil from snow cover during the final melt phase were relatively uniform along the study section (3–5 kg NaCl m⁻¹ of highway). However, the snowpack at all transects retained <50% of applied NaCl, and this shortfall probably reflected direct runoff and infiltration of saline meltwater from the road surface into the adjacent shoulder and right-of-way. Cation exchange with Ca²⁺ in near-surface soils most likely resulted in preferential retention of Na⁺ relative to Cl⁻, although total storage of NaCl in upper soil horizons by winter's end was <15% of deicing salt applications. An environmental tracer (¹⁸O) was used to trace movement of saline meltwater through the unsaturated zone underlying the highway. Average meltwater particle velocities at a site underlain by loam soils were 0.02 m d⁻¹, and ca. 280 mm of water was displaced below a depth of 1.86 m over a 78-day period in the spring and summer of 1995. Sodium ion and chloride ion concentrations in water sampled in late summer 1995 at depths >2 m exceeded 500 mg l⁻¹ and 1000 mg l⁻¹, respectively. Approximately 75% of the net flux of NaCl below the upper soil was retained in the 0–2.8 m depth interval at this site, and results from more permeable soils traversed by the highway indicate an even greater penetration of the annual NaCl application into the unsaturated zone along the moraine. This saline water likely recharges groundwater in this portion of the Oak Ridges Moraine.     

CO2 flux evaluation over the evergreen coniferous and broad-leaved mixed forest in Dinghushan,China

The Dinghushan flux observation site, as one of the four forest sites of ChinaFLUX, aims to acquire long-term measurements of CO₂ flux over a typical southern subtropical evergreen coniferous and broad-leaved mixed forest ecosystem using the open path eddy covariance method. Based on two years of data from 2003 to 2004, the characteristics of temporal variation in CO₂ flux and its response to environmental factors in the forest ecosystem are analyzed. Provided two-dimensional coordinate rotation, WPL correction and quality control, poor energy-balance and underestimation of ecosystem respiration during nighttime implied that there could be a CO₂ leak during the nighttime at the site. Using daytime (PAR > 1.0 μmol⁻¹·m⁻²·s⁻¹) flux data during windy conditions (u* > 0.2 m·s⁻¹), monthly ecosystem respiration (Reco) was derived through the Michaelis-Menten equation modeling the relationship between net ecosystem C0₂ exchange (NEE) and photosynthetically active radiation (PAR). Exponential function was employed to describe the relationship between Reco and soil temperature at 5 cm depth (Ts05), then Reco of both daytime and nighttime was calculated respectively by the function. The major results are: (i) Derived from the Michaelis-Menten equation, the apparent quantum yield (α) was 0.0027±0.0011 mgCO₂·μmol⁻¹ photons, and the maximum photosynthetic assimilation rate (Amax) was 1.102±0.288 mgCO₂·m⁻²·s⁻¹. Indistinctive seasonal variation of α or Amax was consistent with weak seasonal dynamics of leaf area index (LAf) in such a lower subtropical evergreen mixed forest, (ii) Monthly accumulated Reco was estimated as 95.3±21.1 gC·m⁻²mon⁻¹, accounting for about 68% of the gross primary product (GPP). Monthly accumulated WEE was estimated as −43.2±29.6 gC·m⁻²·mon⁻¹. The forest ecosystem acted as carbon sink all year round without any seasonal carbon efflux period. Annual NEE of 2003 and 2004 was estimated as −563.0 and −441.2 gC·m⁻²·a⁻¹ respectively, accounting for about 32% of GPP.

     

A temperature‐precipitation‐based model of thirty‐year mean snowpack accumulation and melt in Oregon,USA

High‐resolution, spatially extensive climate grids can be useful in regional hydrologic applications. However, in regions where precipitation is dominated by snow, snowmelt models are often used to account for timing and magnitude of water delivery. We developed an empirical, nonlinear model to estimate 30‐year means of monthly snowpack and snowmelt throughout Oregon. Precipitation and temperature for the period 1971–2000, derived from 400‐m resolution PRISM data, and potential evapotranspiration (estimated from temperature and day length) drive the model. The model was calibrated using mean monthly data from 45 SNOTEL sites and accurately estimated snowpack at 25 validation sites: R² = 0·76, Nash‐Sutcliffe Efficiency (NSE) = 0·80. Calibrating it with data from all 70 SNOTEL sites gave somewhat better results (R² = 0·84, NSE = 0·85). We separately applied the model to SNOTEL stations located < 200 and ≥ 200 km from the Oregon coast, since they have different climatic conditions. The model performed equally well for both areas. We used the model to modify moisture surplus (precipitation minus potential evapotranspiration) to account for snowpack accumulation and snowmelt. The resulting values accurately reflect the shape and magnitude of runoff at a snow‐dominated basin, with low winter values and a June peak. Our findings suggest that the model is robust with respect to different climatic conditions, and that it can be used to estimate potential runoff in snow‐dominated basins. The model may allow high‐resolution, regional hydrologic comparisons to be made across basins that are differentially affected by snowpack, and may prove useful for investigating regional hydrologic response to climate change. Published in 2011 by John Wiley & Sons, Ltd.