Impact of climate change on runoff from a mid‐latitude mountainous catchment in central Japan

Hydrologic balance in high‐altitude, mid‐latitude mountain areas is important in terms of the water resources available to associated lowlands. This study examined how current and historical shifts in precipitation (P) patterns and concurrent increases in temperature (T) affected runoff (Q) and other hydrologic components in a mid‐latitude mountain catchment of central Japan, using a combination of long‐term data and a simplified hydrologic model, along with their stochastic treatment. The availability of intensive meteorological and hydrological data from the period 1997–2001 allowed the derivation of key relationships for the current climate that tie the forcing term to the parameters or state variables. By using the data recorded in the period 1965–2001, the force for driving the historical simulation was generated. Based on this model and historical shifts in P and T, the probability density functions of Q (pdf(Q)) was computed. A main novelty in this study is that such a stochastic representation, which is useful for considering the influence of projected shifts in environmental factors on the hydrologic budget, was provided. Despite the large increase in the rate of T in winter and spring, pdf(Q) in spring and summer varied appreciably during the time studied mainly because of an increase in snowmelt. An interannual change in whole‐year Q was robust to shifts in T because while Q in spring increased, in summer it decreased, implying a crucial effect of global warming on mountain hydrologic regimes is change in the timing of Q. Copyright © 2009 John Wiley & Sons, Ltd.     

Snowfall interception in a deciduous Nothofagus forest and implications for spatial snowpack distribution

Native Nothofagus forests in the midlatitude region of the Andes Cordillera are notorious biodiversity hot spots, uniquely situated in the Southern Hemisphere such that they develop in snow‐dominated reaches of this mountain range. Spanning a smaller surface area than similar ecosystems, where forests and snow coexist in the Northern Hemisphere, the interaction between vegetation and snow processes in this ecotone has received lesser attention. We present the first systematic study of snow–vegetation interactions in the Nothofagus forests of the Southern Andes, focusing on how the interplay between interception and climate determines patterns of snow water equivalent (SWE) variability. The Valle Hermoso experimental catchment, located in the Nevados de Chillán vicinity, was fitted with eight snow depth sensors that provided continuous measurements at varying elevations, aspect, and forest cover. Also, manual measurements of snow properties were obtained during snow surveys conducted during end of winter and spring seasons for 3 years, between 2015 and 2017. Each year was characterized by distinct climatological conditions, with 2016 representing one of the driest winters on record in this region. Distance to canopy, leaf area index, and total gap area were measured at each observational site. A regression model was built on the basis of statistical analysis of local parameters to model snow interception in this kind of forest. We find that interception implied a 23.2% reduction in snow accumulation in forested sites compared with clearings. The interception in these deciduous trees represents, on average, 23.6% of total annual snowfall, reaching a maximum measured interception value of 13.8‐mm SWE for all snowfall events analysed in this research.     

Spatio‐temporal controls of snowmelt and runoff generation during rain‐on‐snow events in a mid‐latitude mountain catchment

A network of 30 standalone snow monitoring stations was used to investigate the snow cover distribution, snowmelt dynamics, and runoff generation during two rain‐on‐snow (ROS) events in a 40 km² montane catchment in the Black Forest region of southwestern Germany. A multiple linear regression analysis using elevation, aspect, and land cover as predictors for the snow water equivalent (SWE) distribution within the catchment was applied on an hourly basis for two significant ROS flood events that occurred in December 2012. The available snowmelt water, liquid precipitation, as well as the total retention storage of the snow cover were considered in order to estimate the amount of water potentially available for the runoff generation. The study provides a spatially and temporally distributed picture of how the two observed ROS floods developed in the catchment. It became evident that the retention capacity of the snow cover is a crucial mechanism during ROS. It took several hours before water was released from the snowpack during the first ROS event, while retention storage was exceeded within 1 h from the start of the second event. Elevation was the most important terrain feature. South‐facing terrain contributed more water for runoff than north‐facing slopes, and only slightly more runoff was generated at open compared to forested areas. The results highlight the importance of snowmelt together with liquid precipitation for the generation of flood runoff during ROS and the large temporal and spatial variability of the relevant processes. Copyright © 2015 John Wiley & Sons, Ltd.     

Evapotranspiration and runoff in a forest watershed,western Japan

A Note has been published for this article in Hydrological Processes 18(4) 2004, 825. Both water and heat balances were studied in a conifer plantation watershed in south‐west Japan, within the warm‐temperate East Asia monsoon area. Forest cover in the watershed consists mainly of Japanese cedar (Cryptomeria japonica) and Japanese cypress (Chamaecyparis obtusa) plantations. Precipitation and runoff have been observed since 1991, so evapotranspiration can be compared with the water balance. Two meteorological observation towers were built to monitor evapotranspiration in the watershed. The annual average precipitation, amount of runoff and losses were 2166, 1243 and 923 mm, respectively. The evapotranspiration (latent heat flux) agreed well with the water balance losses. The average annual evapotranspiration at the tower built in the centre of the watershed was 902 mm; evapotranspiration at the other tower, further upslope, was 875 mm. The observed evapotranspiration was 39% to 40% of the average precipitation (2166 mm). The mean net radiation was c. 2·6 GJ m^?2 year^?1, and is considered a representative value of the net radiation (Rn) in coniferous plantations in this region. This region is classified in the humid zone based on the ratio of net radiation (Rn) to the energy required to evaporate the rainfall (λR). The mean annual evaporation of canopy‐intercepted water was 356 mm or about 15% of the average precipitation. Copyright © 2003 John Wiley & Sons, Ltd.     

Soil C and N response to changes in winter precipitation in a subalpine forest ecosystem,NW Italy

Among the potential effects of climate change on subalpine forest ecosystems during the winter season, the shift in snowline towards higher altitudes and the increase in frequency of rain events on the snowpack are of particular interest. Here, we present the results of a 2‐year field experiment conducted in a forest stand (Larix decidua) in NW Italy at 2020 m a.s.l. From 2009 to 2011, we monitored soil physical characteristics (temperature and moisture), and soil and soil solution chemistry, in particular carbon (C) and nitrogen (N) forms and their change in time, as affected by simulated late snowpack accumulation and rain on snow events. Late snowpack accumulation determined a stronger effect on soil thermal and moisture regimes than rain on snow events. Also soil chemistry was significantly affected by late snowfall simulation. Although microbial biomass C and N were not reduced by soil freezing, soil contents of the more labile dissolved organic carbon and inorganic N increased when the soil was affected by mild/hard freezing. Variations in the soil solution were shifted with respect to those observed in soil, with an increase in N‐NO₃^? concentrations occurring during spring and summer. This study highlights the potential N loss in subalpine soils under changing environmental conditions driven by a changing climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of snowpack removal on energy balance,melt and runoff in a small supraglacial catchment

Modelling melt and runoff from snow‐ and ice‐covered catchments is important for water resource and hazard management and for the scientific study of glacier hydrology, dynamics and hydrochemistry. In this paper, a distributed, physically based model is used to determine the effects of the up‐glacier retreat of the snowline on spatial and temporal patterns of melt and water routing across a small (0·11 km²) supraglacial catchment on Haut Glacier d'Arolla, Switzerland. The melt model uses energy‐balance theory and accounts for the effects of slope angle, slope aspect and shading on the net radiation fluxes, and the effects of atmospheric stability on the turbulent fluxes. The water routing model uses simplified snow and open‐channel hydrology theory and accounts for the delaying effects of vertical and horizontal water flow through snow and across ice. The performance of the melt model is tested against hourly measurements of ablation in the catchment. Calculated and measured ablation rates show a high correlation (r² = 0·74) but some minor systematic discrepancies in the short term (hours). These probably result from the freezing of surface water at night, the melting of the frozen layer in the morning, and subsurface melting during the afternoon. The performance of the coupled melt/routing model is tested against hourly discharge variations measured in the supraglacial stream at the catchment outlet. Calculated and measured runoff variations show a high correlation (r² = 0·62). Five periods of anomalously high measured discharge that were not predicted by the model were associated with moulin overflow events. The radiation and turbulent fluxes contribute c. 86% and c. 14% of the total melt energy respectively. These proportions do not change significantly as the surface turns from snow to ice, because increases in the outgoing shortwave radiation flux (owing to lower albedo) happen to be accompanied by decreases in the incoming shortwave radiation flux (owing to lower solar incidence angles) and increases in the turbulent fluxes (owing to higher air temperatures and vapour pressures). Model sensitivity experiments reveal that the net effect of snow pack removal is to increase daily mean discharges by c. 50%, increase daily maximum discharges by >300%, decrease daily minimum discharges by c. 100%, increase daily discharge amplitudes by >1000%, and decrease the lag between peak melt rates and peak discharges from c. 3 h to c. 50 min. These changes have important implications for the development of subglacial drainage systems. Copyright © 2002 John Wiley & Sons, Ltd.     

Distribution of fossil seawater and its role in Neogene sedimentary rock landslides in Niigata,eastern marginal region of the Japan Sea

The Neogene marine sedimentary rock area in the eastern marginal region of the Japan Sea is an area with some of the highest landslide densities in Japan. Some of the landslides in this area have been known to involve saline groundwater, which can be the cause of these landslides. In order to demonstrate the relationships between landslides and saline water, topographic, geological, groundwater, and electromagnetic surveys were performed in the eastern marginal region of the Japan Sea. Many landslides and gravitational slope deformations with linear depressions and small scarps were recognized in the study area. The resistivity profile obtained by an electromagnetic survey suggests that there is a wide zonal distribution of saline water with salt concentrations equivalent to seawater at depths of 50–100 m or more and that the groundwater shallower than 50 m has an electrical conductivity of less than 100 mS/m. The shallow resistive groundwater is inferred to be meteoric water that replaced the saline groundwater, which likely weakened the bedrock, resulting in landslides. A ridge of competent tuff overlying mudstone has many linear depressions from gravitational slope deformation and low‐resistivity water to a depth of 600 m, which suggests that the mudstone was weakened by water replacement and deformed under the tuff caprock. The saline groundwater is inferred to be fossil seawater trapped in pores during sediment deposition, which is brought near the ground surface along with rocks by tectonic movement in the hills. Thus, the saline water and its fresh water replacement are among the important basic causes of the landslides. The oil well data obtained in the eastern marginal region of the Japan Sea suggest that such saline water replacement has occurred widely and that replacement is likely one of the predispositions for the frequent landslides there.     

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.     

Application and evaluation of a snowmelt runoff model in the Tamor River basin,Eastern Himalaya using a Markov Chain Monte Carlo (MCMC) data assimilation approach

Previous studies have drawn attention to substantial hydrological changes taking place in mountainous watersheds where hydrology is dominated by cryospheric processes. Modelling is an important tool for understanding these changes but is particularly challenging in mountainous terrain owing to scarcity of ground observations and uncertainty of model parameters across space and time. This study utilizes a Markov Chain Monte Carlo data assimilation approach to examine and evaluate the performance of a conceptual, degree‐day snowmelt runoff model applied in the Tamor River basin in the eastern Nepalese Himalaya. The snowmelt runoff model is calibrated using daily streamflow from 2002 to 2006 with fairly high accuracy (average Nash–Sutcliffe metric ~0.84, annual volume bias < 3%). The Markov Chain Monte Carlo approach constrains the parameters to which the model is most sensitive (e.g. lapse rate and recession coefficient) and maximizes model fit and performance. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall compared with simulations using observed station precipitation. The average snowmelt contribution to total runoff in the Tamor River basin for the 2002–2006 period is estimated to be 29.7 ± 2.9% (which includes 4.2 ± 0.9% from snowfall that promptly melts), whereas 70.3 ± 2.6% is attributed to contributions from rainfall. On average, the elevation zone in the 4000–5500 m range contributes the most to basin runoff, averaging 56.9 ± 3.6% of all snowmelt input and 28.9 ± 1.1% of all rainfall input to runoff. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall versus snowmelt compared with simulations using observed station precipitation. Model experiments indicate that the hydrograph itself does not constrain estimates of snowmelt versus rainfall contributions to total outflow but that this derives from the degree‐day melting model. Lastly, we demonstrate that the data assimilation approach is useful for quantifying and reducing uncertainty related to model parameters and thus provides uncertainty bounds on snowmelt and rainfall contributions in such mountainous watersheds. Copyright © 2013 John Wiley & Sons, Ltd.     

A weighing lysimeter in a regenerating eucalypt forest: Design,construction, and performance

A 36 tonne monolithic weighing lysimeter (3.7 m diameter and 1.5 m deep) was installed in Kioloa State Forest near Batemans Bay, New South Wales, Australia, to provide a continuous record of water use from a regenerating natural eucaplyt community, five years old at installation. The resolution is equivalent to 0.05 mm of evaporation and provides hourly information for diurnal variation of forest evaporation. The sensitivity of 14 ppm compares favourably with those reported for other lysimeters, similar in magnitude, in the U.S.A. Design criteria for selecting, isolating and suspending a sample community ‘in situ’ with an electro-mechanical balance are outlined. Procedures adopted during installation are described with emphasis placed on precautions taken to minimize disturbance and to avert possible changes in character of the isolated natural community. Lysimeter performance was assessed both in terms of its accuracy and its operation as a remote facility with infrequent attention.     

Changes to flood peaks of a mountain river: implications for analysis of the 2013 flood in the Upper Bow River,Canada

The mountain headwater Bow River at Banff, Alberta, Canada, was subject to a large flood in June 2013, over which considerable debate has ensued regarding its probability of occurrence. It is therefore instructive to consider what information long‐term streamflow discharge records provide about environmental change in the Upper Bow River basin above Banff. Though protected as part of Banff National Park, since 1885, the basin has experienced considerable climate and land cover changes, each of which has the potential to impact observations, and hence the interpretations of flood probability. The Bow River at Banff hydrometric station is one of Canada's longest‐operating reference hydrological basin network stations and so has great value for assessing changes in flow regime over time. Furthermore, the station measures a river that provides an extremely important water supply for Calgary and irrigation district downstream and so is of great interest for assessing regional water security. These records were examined for changes in several flood attributes and to determine whether flow changes may have been related to landscape change within the basin as caused by forest fires, conversion from grasslands to forest with fire suppression, and regional climate variations and/or trends. Floods in the Upper Bow River are generated by both snowmelt and rain‐on‐snow (ROS) events, the latter type which include flood events generated by spatially and temporally large storms such as occurred in 2013. The two types of floods also have different frequency characteristics. Snowmelt and ROS flood attributes were not correlated significantly with any climate index or with burned area except that snowmelt event duration correlated negatively to the Pacific Decadal Oscillation. While there is a significant negative trend in all floods over the past 100 years, when separated based on generating process, neither snowmelt floods nor large ROS floods associated with mesoscale storms show any trends over time. Despite extensive changes to the landscape of the basin and in within the climate system, the flood regime remains unchanged, something identified at smaller scales in the region but never at larger scales. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of changes in canopy interception on stream runoff response and recovery following clear-cutting of a Japanese coniferous forest in Fukuroyamasawa Experimental Watershed in Japan

Understanding changes in evapotranspiration during forest regrowth is essential to predict changes of stream runoff and recovery after forest cutting. Canopy interception (Ic) is an important component of evapotranspiration, however Ic changes and the impact on stream runoff during regrowth after cutting remains unclear due to limited observations. The objective of this study was to examine the effects of Ic changes on long-term stream runoff in a regrowth Japanese cedar and Japanese cypress forest following clear-cutting. This study was conducted in two 1-ha paired headwater catchments at Fukuroyamasawa Experimental Watershed in Japan. The catchments were 100% covered by Japanese coniferous plantation forest, one of which was 100% clear-cut in 1999 when the forest was 70 years old. In the treated catchment, annual runoff increased by 301 mm/year (14% of precipitation) the year following clear-cutting, and remained 185 mm/year (7.9% of precipitation) higher in the young regrowth forest for 12–14 years compared to the estimated runoff assuming no clear-cutting. The Ic change was −358 mm/year (17% of precipitation) after cutting and was −168 mm/year (6.7% of precipitation) in the 12–14 years old regrowth forest compared to the observed Ic during the pre-cutting period. Stream runoff increased in all seasons, and the Ic change was the main fraction of evapotranspiration change in all seasons throughout the observation period. These results suggest that the change in Ic accounted for most of the runoff response following forest cutting and the subsequent runoff recovery in this coniferous forest.     

Soil partitioning and surface store controls on spring runoff from a boreal forest peatland basin in north‐central Manitoba,Canada

The mid‐ to high‐boreal forest in Canada occupies the discontinuous permafrost zone, and is often underlain by glaciolacustrine sediments mantled by a highly porous organic mat. The result is a poorly drained landscape dominated by wetlands. Frost‐table dynamics and surface storage conditions help to control runoff contributions from various landscape elements, hydrological linkages between these elements, and basin streamflow during spring snowmelt. Runoff components and pathways in a forested peatland basin were assessed during two spring snowmelts with contrasting input and basin conditions. Runoff from relatively intense melt (up to 16 mm day^?1) on slopes with limited soil thawing combined with large pre‐melt storage in surface depressions to produce high flows composed primarily of meltwater (78% of the 0·29 m³ s^?1 peak discharge) routed over wetland surfaces and through permeable upper peat layers. Melt intensity was less in the subsequent year (maximum of 10 mm day^?1) and active layer development was relatively greater (0·2 m deeper at the end of spring melt), resulting in less slope runoff. Coupling of reduced slope contributions with lower storage levels in basin wetlands led to relatively subdued streamflows dominated by older water (73% of the 0·09 m³ s^?1 peak discharge) routed through less‐permeable deeper peat layers and mineral soil. Interannual differences in runoff conditions provide important insight for the development of distributed hydrological models for boreal forest basins and into potential influences on biogeochemical cycling in this landscape under a warming climate. Copyright © 2001 John Wiley & Sons, Ltd.     

Seasonal soil water dynamics in the jarrah forest,Western Australia. II: Results from a site with fine-textured soil profiles

Seasonal soil water dynamics were measured at a fine-textured, upslope site within the jarrah forest of southwest Western Australia and compared to the results from a coarse-textured hillslope transect. Gravity drainage dominated during winter and early spring. This reversed in early summer and an upward potential gradient was observed to 7 m depth. A shallow ephemeral saturation zone was observed above a clay pan at 1.5 m depth. This saturation zone persisted through late winter and early spring, contrasting with the short-lived saturation in the duricrust on the hillslope transect. The annual maximum to minimum unsaturated soil water storage was about 530 mm, 50 mm greater than the hillslope transect and higher than most values reported elsewhere in Australia. Significant soil water content changes following winter rain were generally restricted to 6 m but at one site occurred to 9 m. These depths were significantly less than the coarser-textured hillslope transect. Soil water drying rates averaged 5 mm day^?1 during extended dry periods compared to 3.5 mm day^?1 on the hillslope transect. The drying rate occurred uniformly through the profile until late summer when a significant decrease in the upper 3 m was observed.     

Seasonal soil water dynamics in the jarrah forest,Western Australia. I: Results from a hillslope transect with coarse-textured soil profiles

Seasonal soil water dynamics were measured on a hillslope transect in the jarrah forest of southwest Western Australia over the period 1984-86 using mercury manometer tensiometers, gypsum blocks, and a neutron moisture meter. The soil water potential gradients indicated downward vertical drainage flux through winter and spring. There was generally a change to an upwards flux in early summer which was sustained through to autumn. A shallow ephemeral saturation zone was identified in and above a duricrust layer, lasting up to three days after heavy, late winter rainfall. The annual maximum to minimum unsaturated soil water storage on the hillslope was approximately 400 mm to 6 m depth and 480 mm to 15 m depth. This did not change significantly in years of substantially different winter rainfall. The magnitude of seasonal soil water storage was similar to other forested areas with deep soil profiles. The depth of observable infiltration was dependent on annual rainfall. This was consistent with the observation that groundwater levels responded to rainfall over the whole hillslope in wet years but only responded on the lower slopes in dry years. The average summer drying rate of the soil profile to 6 m depth of 3.5 mm day^?1 was within the range of values reported for forests elsewhere. In late summer, following an extended drought period, the drying rate decreased downslope but increased midslope.     

Long-term precipitation and stream discharge records at seven forested experimental watersheds along a latitudinal transect in Japan: Jozankei,Kamabuchi, Takaragawa,Tsukuba, Tatsunokuchi-yama,Kahoku and Sarukawa

This data note introduces a database of long-term daily total precipitation and stream discharge data for seven forested watersheds in Japan that have been continuously monitored by the Forestry and Forest Products Research Institute. Three of the watersheds started data collection in the 1930s. Forest cover across the sites ranges from cool to warm temperate regions with the latitude spanning from 31 to 44° N and annual precipitation ranging from 1200 to 3000 mm yr⁻¹. The effects of vegetation change via clearcutting, thinning and forest fire (among other stressors) on stream discharge can be analysed from the long-term observation sites. Moreover, this multi-site dataset allows for inter- and intra-site comparisons of annual water loss (difference of annual precipitation and stream discharge). These long-term datasets can provide comprehensive insights into the effects of climate change and other stressors on forested ecosystems, not only in Japan but across a spectrum of forest types, if combined with other long-term records from other forested watersheds across the world.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.