Long‐term investigations of the snow cover in a subalpine semi‐forested catchment

To improve spring runoff forecasts from subalpine catchments, detailed spatial simulations of the snow cover in this landscape is obligatory. For more than 30 years, the Swiss Federal Research Institute WSL has been conducting extensive snow cover observations in the subalpine watershed Alptal (central Switzerland). This paper summarizes the conclusions from past snow studies in the Alptal valley and presents an analysis of 14 snow courses located at different exposures and altitudes, partly in open areas and partly in forest. The long‐term performance of a physically based numerical snow–vegetation–atmosphere model (COUP) was tested with these snow‐course measurements. One single parameter set with meteorological input variables corrected to the prevailing local conditions resulted in a convincing snow water equivalent (SWE) simulation at most sites and for various winters with a wide range of snow conditions. The snow interception approach used in this study was able to explain the forest effect on the SWE as observed on paired snow courses. Finally, we demonstrated for a meadow and a forest site that a successful simulation of the snowpack yields appropriate melt rates. Copyright © 2006 John Wiley & Sons, Ltd.     

Determination of snow water equivalent over the eastern part of Turkey using passive microwave data

Information on regional snow water equivalent (SWE) is required for the management of water generated from snowmelt. Modeling of SWE in the mountainous regions of eastern Turkey, one of the major headwaters of Euphrates–Tigris basin, has significant importance in forecasting snowmelt discharge, especially for optimum water usage. An assimilation process to produce daily SWE maps is developed based on Helsinki University of Technology (HUT) model and AMSR‐E passive microwave data. The characteristics of the HUT emission model are analyzed in depth and discussed with respect to the extinction coefficient function. A new extinction coefficient function for the HUT model is proposed to suit models for snow over mountainous areas. Performance of the modified model is checked against the original, other modified cases and ground truth data covering the 2003–2007 winter periods. A new approach to calculate grain size and density is integrated inside the developed data assimilation process. An extensive validation was successfully performed by means of snow data measured at ground stations during the 2008–2010 winter periods. The root mean square error of the data set for snow depth and SWE between January and March of the 2008–2010 periods compared with the respective AMSR‐E footprints indicated that errors for estimated snow depth and predicted SWE values were 16.92 cm and 40.91 mm, respectively, for the 3‐year period. Validation results were less satisfactory for SWE less than 75.0 mm and greater than 150.0 mm. An underestimation for SWE greater than 150 mm could not be resolved owing to the microwave signal saturation that is observed for dense snowpack. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaluation of snow models in terrestrial biosphere models using ground observation and satellite data: impact on terrestrial ecosystem processes

Snow is important for water management, and an important component of the terrestrial biosphere and climate system. In this study, the snow models included in the Biome‐BGC and Terrestrial Observation and Prediction System (TOPS) terrestrial biosphere models are compared against ground and satellite observations over the Columbia River Basin in the US and Canada and the impacts of differences in snow models on simulated terrestrial ecosystem processes are analysed. First, a point‐based comparison of ground observations against model and satellite estimates of snow dynamics are conducted. Next, model and satellite snow estimates for the entire Columbia River Basin are compared. Then, using two different TOPS simulations, the default TOPS model (TOPS with TOPS snow model) and the TOPS model with the Biome‐BGC snow model, the impacts of snow model selection on runoff and gross primary production (GPP) are investigated. TOPS snow model predictions were consistent with ground and satellite estimates of seasonal and interannual variations in snow cover, snow water equivalent, and snow season length; however, in the Biome‐BGC snow model, the snow pack melted too early, leading to extensive underpredictions of snow season length and snow covered area. These biases led to earlier simulated peak runoff and reductions in summer GPP, underscoring the need for accurate snow models within terrestrial ecosystem models. Copyright © 2007 John Wiley & Sons, Ltd.     

Quantifying the role of the snowpack in generating water available for run‐off during rain‐on‐snow events from snow pillow records

Rain‐on‐snow events have generated major floods around the world, particularly in coastal, mountainous regions. Most previous studies focused on a limited number of major rain‐on‐snow events or were based primarily on model results, largely due to a lack of long‐term records from lysimeters or other instrumentation for quantifying event water balances. In this analysis, we used records from five automated snow pillow sites in south coastal British Columbia, Canada, to reconstruct event water balances for 286 rain‐on‐snow events over a 10‐year period. For large rain‐on‐snow events (event rainfall >40 mm), snowmelt enhanced the production of water available for run‐off (WAR) by approximately 25% over rainfall alone. For smaller events, a range of antecedent and meteorological factors influenced WAR generation, particularly the antecedent liquid water content of the snowpack. Most large events were associated with atmospheric rivers. Rainfall dominated WAR generation during autumn and winter events, whereas snowmelt dominated during spring and summer events. In the majority of events, the sensible heat of rain contributed less than 10% of the total energy consumed by snowmelt. This analysis illustrated the importance of understanding the amount of rainfall occurring at high elevations during rain‐on‐snow events in mountainous regions.     

Estimating the distribution of snow water equivalent and snow extent beneath cloud cover in the Salt–Verde River basin,Arizona

The temporal and spatial continuity of spatially distributed estimates of snow‐covered area (SCA) are limited by the availability of cloud‐free satellite imagery; this also affects spatial estimates of snow water equivalent (SWE), as SCA can be used to define the extent of snow telemetry (SNOTEL) point SWE interpolation. In order to extend the continuity of these estimates in time and space to areas beneath the cloud cover, gridded temperature data were used to define the spatial domain of SWE interpolation in the Salt–Verde watershed of Arizona. Gridded positive accumulated degree‐days (ADD) and binary SCA (derived from the Advanced Very High Resolution Radiometer (AVHRR)) were used to define a threshold ADD to define the area of snow cover. The optimized threshold ADD increased during snow accumulation periods, reaching a peak at maximum snow extent. The threshold then decreased dramatically during the first time period after peak snow extent owing to the low amount of energy required to melt the thin snow cover at lower elevations. The area having snow cover at this later time was then used to define the area for which SWE interpolation was done. The area simulated to have snow was compared with observed SCA from AVHRR to assess the simulated snow map accuracy. During periods without precipitation, the average commission and omission errors of the optimal technique were 7% and 11% respectively, with a map accuracy of 82%. Average map accuracy decreased to 75% during storm periods, with commission and omission errors equal to 11% and 12% respectively. The analysis shows that temperature data can be used to help estimate the snow extent beneath clouds and therefore improve the spatial and temporal continuity of SCA and SWE products. Copyright © 2004 John Wiley & Sons, Ltd.     

Uncertainty of snow water equivalent retrieved from AMSR‐E brightness temperature in northeast Asia

Accurate estimation of snow water equivalent (SWE) has been significantly recognized to improve management and analyses of water resource in specific regions. Although several studies have focused on developing SWE values based on remotely sensed brightness temperatures obtained by microwave sensor systems, it is known that there are still a number of uncertainties in SWE values retrieved from microwave radiometers. Therefore, further research for improving remotely sensed SWE values including global validation should be conducted in unexplored regions such as Northeast Asia. In this regard, we evaluated SWE through comparison of values produced by the Advanced Microwave Scanning Radiometer Earth Observing System (AMSR‐E) from December 2002 to February 2011 with in situ SWE values converted from snow‐depth observation data from four regions in the South Korea. The results from three areas showed similarities which indicated that the AMSR‐E SWE values were overestimated when compared with in situ SWE values, and their Mean Absolute Errors (MAE) by month were relatively small (1.1 to 6.5 mm). Contrariwise, the AMSR‐E SWE values of one area were significantly underestimated when compared with in situ SWE values and the MAE were much greater (4.9 to 35.2 mm). These results were closely related to AMSR‐E algorithm‐related error sources, which we analyzed with respect to topographic characteristics and snow properties. In particular, we found that snow density data used in the AMSR‐E SWE algorithm should be based on reliable in situ data as the current AMSR‐E SWE algorithm cannot reflect the spatio‐temporal variability of snow density values. Additionally, we derived better results considering saturation effect of AMSR‐E SWE. Despite the demise of AMSR‐E, this study's analysis is significant for providing a baseline for the new sensor and suggests parameters important for obtaining more reliable SWE. Copyright © 2013 John Wiley & Sons, Ltd.     

A conceptual model of water yield effects from beetle‐induced tree death in snow‐dominated lodgepole pine forests

In regions of western North America with snow‐dominated hydrology, the presence of forested watersheds can significantly influence streamflow compared to areas with other vegetation cover types. Widespread tree death in these watersheds can thus dramatically alter many ecohydrologic processes including transpiration, canopy solar transmission and snow interception, subcanopy wind regimes, soil infiltration, forest energy storage and snow surface albedo. One of the more important causes of conifer tree death is bark beetle infestation, which in some instances will kill nearly all of the canopy trees within forest stands. Since 1996, an ongoing outbreak of bark beetles (Coleoptera: Scolytidae) has caused widespread mortality across more than 600,000 km² of coniferous forests in western North America, including numerous Rocky Mountain headwaters catchments with high rates of lodgepole pine (Pinus contorta) mortality from mountain pin beetle (Dendroctonous ponderosae) infestations. Few empirical studies have documented the effects of MPB infestations on hydrologic processes, and little is known about the direction and magnitude of changes in water yield and timing of runoff due to insect‐induced tree death. Here, we review and synthesize existing research and provide new results quantifying the effects of beetle infestations on canopy structure, snow interception and transmission to create a conceptual model of the hydrologic effects of MPB‐induced lodgepole pine death during different stages of mortality. We identify the primary hydrologic processes operating in living forest stands, stands in multiple stages of death and long‐dead stands undergoing regeneration and estimate the direction of change in new water yield. This conceptual model is intended to identify avenues for future research efforts. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaluating the functionality and streamflow impacts of explicitly modelling forest–snow interactions and canopy gaps in a distributed hydrologic model

Many plot‐scale studies have shown that snow‐cover dynamics in forest gaps are distinctly different from those in open and continuously forested areas, and forest gaps have the potential to alter the magnitude and timing of snowmelt. However, the watershed‐level impacts of canopy gap treatment on streamflows are largely unknown. Here, we present the first research that explicitly assesses the impact of canopy gaps on seasonal streamflows and particularly late‐season low flows at the watershed scale. To explicitly model forest–snow interactions in canopy gaps, we made major enhancements to a widely used distributed hydrologic model, distributed hydrology soil vegetation model, with a canopy gap component that represents physical processes of snowpack evolution in the forest gap separately from the surrounding forest on the subgrid scale (within a grid typically 10–150 m). The model predicted snow water equivalent using the enhanced distributed hydrology soil vegetation model showed good agreement (R² > 0.9) with subhourly snow water equivalent measurements collected from open, forested, and canopy gap sites in Idaho, USA. Compared with the original model that does not account for interactions between gaps and surrounding forest, the enhanced model predicted notably later melt in small‐ to medium‐size canopy gaps (the ratio of gap radius (r) to canopy height (h) ≤ 1.2), and snow melt rates exhibited great sensitivity to changing gap size in medium‐size gaps (0.5 ≤ r/h ≤ 1.2). We demonstrated the watershed‐scale implications of canopy gaps on streamflow in the snow‐dominated Chiwawa watershed, WA, USA. With 24% of the watershed drainage area (about 446 km²) converted to gaps of 60 m diameter, the mean annual 7‐day low flow was increased by 19.4% (i.e., 0.37 m³/s), and the mean monthly 7‐day low flows were increased by 13.5% (i.e., 0.26 m³/s) to 40% (i.e., 1.76 m³/s) from late summer through fall. Lastly, in practical implementation of canopy gaps with the same total gap areas, a greater number of distributed small gaps can have greater potential for longer snow retention than a smaller number of large gaps.     

Estimating the distribution of snow water equivalent using remotely sensed snow cover data and a spatially distributed snowmelt model: A multi-resolution,multi-sensor comparison

Time series of fractional snow covered area (SCA) estimates from Landsat Enhanced Thematic Mapper (ETM+), Moderate Resolution Imaging Spectroradiometer (MODIS), and Advanced Very High Resolution Radiometer (AVHRR) data were combined with a spatially distributed snowmelt model to reconstruct snow water equivalent (SWE) in the Rio Grande headwaters (3419 km²). In this reconstruction approach, modeled snowmelt over each pixel is integrated during the period of satellite-observed snow cover to estimate SWE. Due to underestimates in snow cover detection, maximum basin-wide mean SWE using MODIS and AVHRR were, respectively, 45% and 68% lower than SWE estimates obtained using ETM+ data. The mean absolute error (MAE) of SWE estimated at 100-m resolution using ETM+ data was 23% relative to observed SWE from intensive field campaigns. Model performance deteriorated when MODIS (MAE = 50%) and AVHRR (MAE = 89%) SCA data were used. Relative to differences in the SCA products, model output was less sensitive to spatial resolution (MAE = 39% and 73% for ETM+ and MODIS simulations run at 1 km resolution, respectively), indicating that SWE reconstructions at the scale of MODIS acquisitions may be tractable provided the SCA product is improved. When considering tradeoffs between spatial and temporal resolution of different sensors, our results indicate that higher spatial resolution products such as ETM+ remain more accurate despite the lower frequency of acquisition. This motivates continued efforts to improve MODIS snow cover products.     

Detecting the flow relationships between deep and shallow aquifers in an exploited groundwater system,using long‐term monitoring data and quantitative hydrogeology: the Acque Albule basin case (Rome,Italy)

Five years of hydrogeological monitoring and field activities performed in the complex hydrogeological system of the Acque Albule basin (AAB) were conducted to define the hydrogeological setting, the relationship between deep and shallow aquifers and a conceptual groundwater flow model of this exploited area using conventional quantitative techniques. The basin, which is located close to Rome (Italy) on the west side of the Apennine chain and just north of the Colli Albani volcano, subsided after development of a north–south fault system (about 115 000 y bp). The AAB experiences intense hydrothermal activity, which has produced a large travertine deposit (80‐m thick). The travertine deposit constitutes a fractured aquifer that is the final destination of more than 5 m³ s^‐1 of water and is strongly dewatered by quarry activities. The complex hydrogeology of this basin was investigated, revealing two main hydraulically connected aquifers, one thermalised and partly confined into the limestone bedrock and one unconfined in the travertine. The two aquifers are separated by a non‐continuous clayey aquiclude. The hydrogeological survey and geological characterisation contributed to the development of the groundwater flow conceptual model. Analysis and comparison of the monitored levels highlighted the pattern of flow between the deep and shallow parts of the flow system. Copyright © 2012 John Wiley & Sons, Ltd.     

Measurement of evaporation from the forest floor in a deciduous forest throughout the year using microlysimeter and closed‐chamber systems

Evaporation from the forest floor (E_FF) in a deciduous broadleaf forest was measured using microlysimeter and closed‐chamber systems. The microlysimeter was used at six points in the experimental basin, and measurements gave different E_FF values at different points. This could be attributed to the local photoenvironment of each sampling point, rather than to litter conditions, if the spatial variation in air temperature (Ta) or vapour pressure deficit (VPD) at the forest floor was small within this basin. A detachable microlysimeter measured condensation in the litter layer during the night, indicating that the litter layer, as well as the mulch layer, played a role in preventing evaporation from the soil layer. The closed‐chamber system made it possible to continuously measure long‐term E_FF. E_FF was closely related to VPD; even during the night, when solar radiation was zero, E_FF amounted to 14·0% of the daily E_FF. The daily E_FF was 0·20 ± 0·13 mm day^?1 during the study period, with two seasonal peaks: in late spring (0·31 mm day^?1 in April) and early fall (0·22 mm day^?1 in September). The former peak has been reported from two deciduous forests in Japan and is strongly related to the solar radiation reaching the forest floor when the trees are dormant. Copyright © 2008 John Wiley & Sons, Ltd.     

Long‐term monitoring of a mercury contaminated estuary (Ria de Aveiro,Portugal): the effect of weather events and management in mercury transport

The main aim of this research was to assess the mercury transport from an estuarine basin with a background of anthropogenic contamination during a spring tidal cycle (year 2009) and compare it with two previous tidal cycles (years 1994 and 1999), as part of a long‐term monitoring program. Results showed that effective mercury transport occurs both in the dissolved and particulate fractions (0.18 and 0.20 kg per tidal cycle, respectively), and despite an overall decrease in environmental contamination, results more than double previous findings on particulate transport in the system. These findings result essentially from changes in the tidal prism (net export of 2 million m³ of water), given that both dissolved and particulate concentrations did not increase over time. Hydrodynamic simulations were performed to evaluate the effect of physical disturbance (dredging) and weather events (increased freshwater flow) in these processes, and results suggest the increased freshwater flow into the system as the main forcing function for the mercury transport increment. These results highlight the importance of long‐term monitoring programs, since despite an overall improvement in local contamination levels, the enhancement of transport processes through hydrological changes increases environmental pressure away from the contamination source. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating river discharge from very high‐resolution satellite data: a case study in the Yangtze River,China

The measurement of river discharge is necessary for understanding many water‐related issues. Traditionally, river discharge is estimated by measuring water stage and converting the measurement to discharge by using a stage–discharge rating curve. Our proposed method for the first time couples the measurement of water‐surface width with river width–stage and stage–discharge rating curves by using very high‐resolution satellite data. We used it to estimate the discharge in the Yangtze (Changjiang) River as a case study. The discharges estimated at four stations from five QuickBird‐2 images matched the ground observation data very well, demonstrating that the proposed approach can be regarded as ancillary to traditional field measurement methods or other remote methods to estimate river discharge. Copyright © 2004 John Wiley & Sons, Ltd.     

Impact of an inter‐basin water transfer and reservoir operation on a karst open streamflow hydrological regime: an example from the Dinaric karst (Croatia)

This paper describes the hydrological changes caused by inter‐basin water transfer and the reservoir development on the hydrological regimes of two rivers. The Sabljaki Reservoir in the Zagorska Mre?nica River and the Bukovik Reservoir in the upper Dobra River began operation in 1959. Both are part of the hydroelectric power plant (HEPP) Gojak, whose installed capacity is 50 m³/s. Their water volumes at the spillway altitudes of 320·10 and 320·15 m a. s. l. are 3·3 × 10⁶ and 0·24 × 10⁶ m³ respectively. Both the Dobra and Mre?nica Rivers are losing, sinking and underground karst rivers. A 9376‐m‐long tunnel provides water from the Sabljaki Reservoir to the HEPP Gojak, which was constructed in the Lower Dobra River. The Sabljaki Reservoir is located in the Pla?ki karst polje, while the Bukovik Reservoir is located in the neighbouring Ogulin karst polje. The consequences of the inter‐basin water transfer are strong and have caused abrupt changes in the hydrological regimes of the downstream sections of both rivers. At the same time, the construction and development of both the reservoirs have also caused hydrological changes to the upstream section of the Upper Dobra River. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of forest on annual water yield in the mountains of an arid inland river basin: a case study in the Pailugou catchment on northwestern China's Qilian Mountains

The effect of forests on annual water yield is an unresolved central issue in forest hydrology despite years of study. There has been a particular shortage of research in the mountains of arid inland river basins. In the present study, we examined the effects of forests on hydrology using data on precipitation, evaporation, canopy interception, transpiration, and runoff from 1994 to 2008 for the Pailugou catchment of northwestern China's Qilian Mountains. We modelled the water balance to assess the contribution of different vegetation types to annual water yield. In our study area, Picea crassifolia forest covered 38·5% of the catchment area, but contributed little to annual water yield. For an annual average precipitation of 407·1 mm (from 2003 to 2008) at an elevation of 2700 m, the runoff depth from the forest was 11·6 mm, accounting for only 3·5% of total annual water yield of the catchment. For an annual average precipitation of 374·1 mm (from 1994 to 2002), the runoff depth from the forest was ? 14·3 mm (i.e. 5·9% of total annual water yield of the catchment was consumed to sustain tree growth). This has significant implications, because forests are increasingly being planted in the Qilian Mountains, and this may decrease the downstream water supply. Thus, the relationship between the forest and water yield must be better understood to permit the establishment of an appropriate regional level of forest cover. Copyright © 2011 John Wiley & Sons, Ltd.     

Investigation of runoff generation in a pristine,poorly gauged catchment in the Chilean Andes I: A multi‐method experimental study

Catchment scale hydrological process studies in southern Chile are of special interest as little research at this scale has been carried out in this region. In particular, the young volcanic ash soils, which are typical for this area, are not well understood in their hydrological behaviour. In addition, extensive land use changes require detailed knowledge of hydrological processes in disturbed as well as undisturbed catchments in order to estimate resulting risks of erosion, eutrophication, floods and droughts. This study focuses on data collection and experimental determination of relevant processes in an undisturbed forested catchment in the Andes of southern Chile. The here gained understanding of runoff generation can serve as a reference for comparison with sites subject to human intervention, improving estimation of the effects of land use change. Owing to the lack of long‐term data for this catchment it was necessary to replace long time series by a multitude of experimental methods covering as many aspects of the runoff generation process as possible. The methods used in this investigation include: measurements of streamflow, rainfall, throughfall, water chemistry, soil water dynamics, groundwater dynamics, soil physics, soil mineralogy, geo‐electrical sounding, and tracer techniques. Methods and equipment used during field campaigns are described and evaluated for usefulness versus expenditure (labour and financial costs). Selected results and the hypotheses developed from these findings are presented. The results suggest the importance of fast processes for rainfall runoff response on the one hand as well as considerable dampening effects of a large subsurface storage on the other hand. Copyright © 2008 John Wiley & Sons, Ltd.     

Short‐term variability in isotopic composition of precipitation: A case study from the Midwestern United States

Oxygen and deuterium isotopes in precipitation were analysed to define local isotopic trends in Iowa, US. The area is far inland from an oceanic source and the observed averages of δ¹⁸O and δ D are ? 6·43‰ and ? 41·35‰ for Ames, ? 7·53‰ and ? 51·33‰ for Cedar Falls, and ? 6·01‰ and ? 38·19‰ for Iowa City, respectively. Although these data generally follow global trends, they are different when compared to a semi‐arid mid‐continental location in North Platt, Nebraska. The local meteoric water lines of Iowa are δ D = 7·68 δ¹⁸O + 8·0 for Ames, δ D = 7·62 δ¹⁸O + 6·07 for Cedar Falls, and δ D = 7·78 δ¹⁸O + 8·61 for Iowa City. The current Iowa study compares well with a study conducted in Ames, Iowa, 10 years earlier. The differences between Iowa and Nebraska studies are attributed to a variable climate across the northern Great Plains ranging from sub‐humid in the east to semi‐arid in the west. Iowa being further east in the region is more strongly influenced by a moist sub‐humid to humid climate fed by the tropical air stream from the Gulf of Mexico. The average d‐excess values are 10·06‰ for Ames, 8·92‰ for Cedar Falls and 9·92‰ for Iowa City. Eighty seven percent of the samples are within the global d‐excess range of 0‰ and 20‰. The results are similar to previous studies, including those by National Atmospheric Deposition Programs and International Atomic Energy Agency. It appears that the impact of recycled water or secondary evaporation on δ¹⁸O values of area precipitation is minimal. Copyright © 2008 John Wiley & Sons, Ltd.     

Multi‐site calibration using a grid‐based event rainfall–runoff model: a case study of the upstream areas of the Nakdong River basin in Korea

Most runoff analyses using a grid‐based distributed model use one parameter group calibrated at the outlet of a watershed, instead of dividing the watershed into subwatersheds. Significant differences between the observed value and the simulation result of the subwatersheds can occur if just one parameter group is used in all subwatersheds that have different hydrological characteristics from each other. Therefore, to improve the simulation results of the subwatersheds within a watershed, a model calibrated at every subwatershed needs to be used to reflect the characteristics of each subwatershed. In this study, different parameter groups were set up for one or two sites using a distributed model, the GRM (Grid based Rainfall‐runoff Model), and the evaluations were based on the results of rainfall–runoff analysis, which uses a multi‐site calibration (MSC) technique to calibrate the model at the outlet of each site. The Hyangseok watershed in Naeseong River, which is a tributary of Nakdong River in Korea, was chosen as the study area. The watershed was divided into five subwatersheds each with a subwatershed outlet that was applied to the calibration sites . The MSC was applied for five cases. When a site was added for calibration in a watershed, the runoff simulation showed better results than the calibration of only one site at the most downstream area of the watershed. The MSC approach could improve the simulation results on the calibrated sites and even on the non‐calibrated sites, and the effect of MSC was improved when the calibrated site was closer to the runoff site. Copyright © 2014 John Wiley & Sons, Ltd.