The Kiryu Experimental Watershed: 50-years of rainfall-runoff data for a forest catchment in central Japan

The Kiryu Experimental Catchment (KEW) is a small (5.99 ha) forest catchment located in Shiga Prefecture, central Japan (34°58′ N, 136°00′ E; www.bluemoon.kais.kyoto-u.ac.jp/kiryu/contents.html ). Around this area, forest devastation occurred from ca. 1250 to ca. 150 years ago because of overuse of forest and timbers. Then, hillside forestation was carried out for more than 100 years to prevent soil erosion and support the timber industry, and consequently, most of this area is now covered with plantation forests mainly by Chamaecyparis obtusa Sieb. et Zucc. (Japanese cypress) planted around 1960's. This plantation forest is not actively managed. The KEW is one of the leading experimental forests with long-term monitoring data in Japan. Research in the KEW began in 1967 to elucidate the hydrological and biogeochemical processes in the forested catchment in relation to climate, geology, soil, and vegetation growth. Since then, the long-term hydrological data of precipitation, runoff and sediment transport are continuously monitoring. In this study, we provide the data and preliminarily discuss the rainfall–runoff patterns and sediment transport through 50 years in the KEW. The annual precipitation and the maximum daily rainfall have been greater than the average over the last decade. In response to the rainfall patterns, the ratio of annual direct runoff to precipitation was also larger in the last decade. The sediment transport in this decade was consequently larger than the preceding decades. Our data presented here suggest that a close relationship exists between the climate condition, rainfall–runoff response, sediment dynamics, as well as a slowly progressing change of forest condition.     

Using stable isotopes to identify major flow pathways in a permafrost influenced alpine meadow hillslope during summer rainfall period

Global warming has leaded to permafrost degradation, with potential impacts on the runoff generation processes of permafrost influenced alpine meadow hillslope. Stable isotopes have the potential to trace the complex runoff generation processes. In this study, precipitation, hillslope surface and subsurface runoff, stream water, and mobile soil water (MSW) at different hillslope positions and depths were collected during the summer rainfall period to analyse the major flow pathway based on stable isotopic signatures. The results indicated that (a) compared with precipitation, the δ²H values of MSW showed little temporal variation but strong heterogeneity with enriched isotopic ratios at lower hillslope positions and in deeper soil layers. (b) The δ²H values of middle-slope surface runoff and shallow subsurface flow were similar to those of precipitation and MSW of the same soil layer, respectively. (c) Middle-slope shallow subsurface flow was the major flow pathway of the permafrost influenced alpine meadow hillslope, which turned into surface runoff at the riparian zone before contributing to the streamflow. (d) The slight variation of δ²H values in stream water was shown to be related to mixing processes of new water (precipitation, 2%) and old water (middle-slope shallow subsurface flow, 98%) in the highly transmissive shallow thawed soil layers. It was inferred that supra-permafrost water levels would be lowered to a less conductive, deeper soil layer under further warming and thawing permafrost, which would result in a declined streamflow and delayed runoff peak. This study explained the “rapid mobilization of old water” paradox in permafrost influenced alpine meadow hillslope and improved our understanding of permafrost hillslope hydrology in alpine regions.     

Hydrological regimes in a tropical valley of New Caledonia (SW Pacific): Impacts of wildfires and invasive fauna

In New Caledonia wildfires and invasive mammals (deer and wild pigs) constitute the major agents of land surface degradation. Our study reveals the linkage between land cover and water balance on the northeast coast of New Caledonia (2400 mm annual rainfall) located on a micaschist basement. The hydrological regime of characteristic and representative land surfaces is assessed using a 1-year record from three 100 m² plots each, located in a forest area degraded by an invasive fauna, in a woody savannah which is regularly burned, and in a healthy forest area. The three plots present highly contrasting hydrological regimes, with annual and maximum runoff/rain ratios during a rain event of, respectively, 0.82, 0.16, 0.03, and 2.7, 0.7, 0.2, for the degraded forest, the savannah and the healthy forest. Such results suggest that subsurface flow originating from the contributing area above the degraded forest plot should exfiltrate inside the plot. A conceptual model for the degraded forest plot shows that water exfiltrating inside the plot represents 61% of the observed runoff. In savannahs, water should mainly be transferred downstream by subsurface flow within a thick organic soil layer limited by an impervious clay layer at a 20–30 cm depth. Savannahs are generally located above forests and generate the transfer of rainwater to downslope forests. Exfiltration into the forests can be the result of this transfer and depends on the thickness and permeability of the forest topsoils and on topographic gradients. Water exfiltration in forest areas highly degraded by pigs and deer enhances erosion and increases further degradation. It probably also limits percolation in the areas located downstream by increasing the amount of superficial runoff concentrated in gullies.     

Accelerated continental-scale snowmelt and ecohydrological impacts in the four largest Siberian river basins in response to spring warming

River runoff from the four largest Siberian river basins (the Ob, Yenisei, Lena, and Kolyma) considerably contributes to freshwater flux into the Arctic Ocean from the Eurasian continent. However, the effects of variation in snow cover fraction on the ecohydrological variations in these basins are not well understood. In this study, we analysed the spatiotemporal variability of the maximum snow cover fraction (SCF_max) in the four Siberian river basins. We compared the SCF_max from 2000 to 2016 with data in terms of monthly temperature and precipitation, night-time surface temperatures, the terrestrial water storage anomaly (TWSA), the normalised difference vegetation index (NDVI), and river runoff. Our results exhibit a decreasing trend in the April SCF_max values since 2000, largely in response to warming air temperatures in April. We identified snowmelt water as the dominant control on the observed increase in the runoff contribution in May across all four Siberian river basins. In addition, we detected that the interannual river runoff was predominantly controlled by interannual variations in the TWSA. The NDVI in June was strongly controlled by the timing of the snowmelt along with the surface air temperature and TWSA in June. The rate of increase in the freshwater flux from the four Siberian rivers decreased from 2000 to 2016, exhibiting large interannual variations corresponding to interannual variations in the TWSA. However, we identified a clear increase trend in the freshwater flux of ~4 km³/year when analysing the long-term 39-year historical record (1978–2016). Our results suggest that continued global warming will accelerate the transition towards the earlier timing of snowmelt and spring freshwater flux into the Arctic Ocean. Our findings also highlight the effects of earlier snowmelt on ecohydrological changes in the Northern Hemisphere.     

Stability of model performance and parameter values on two catchments facing changes in climatic conditions

Abstract

Hydrological models are often used for studying the hydrological effects of climate change; however, the stability of model performance and parameter values under changing climate conditions has seldom been evaluated and compared. In this study, three widely-used rainfall–runoff models, namely the SimHYD model, the HBV model and the Xin’anjiang model, are evaluated on two catchments subject to changing climate conditions. Evaluation is carried out with respect to the stability in their performance and parameter values in different calibration periods. The results show that (a) stability of model performance and parameter values depends on model structure as well as the climate of catchments, and the models with higher performance scores are more stable in changing conditions; (b) all the tested models perform better on a humid catchment than on an arid catchment; (c) parameter values are also more stable on a humid catchment than on an arid catchment; and (d) the differences in stability among models are somewhat larger in terms of model efficiency than in model parameter values.     

Evaluation of dynamically downscaled reanalysis precipitation data for hydrological application

Skilful and reliable precipitation data are essential for seasonal hydrologic forecasting and generation of hydrological data. Although output from dynamic downscaling methods is used for hydrological application, the existence of systematic errors in dynamically downscaled data adversely affects the skill of hydrologic forecasting. This study evaluates the precipitation data derived by dynamically downscaling the global atmospheric reanalysis data by propagating them through three hydrological models. Hydrological models are calibrated for 28 watersheds located across the southeastern United States that is minimally affected by human intervention. Calibrated hydrological models are forced with five different types of datasets: global atmospheric reanalysis (National Centers for Environmental Prediction/Department of Energy Global Reanalysis and European Centre for Medium‐Range Weather Forecasts 40‐year Reanalysis) at their native resolution; dynamically downscaled global atmospheric reanalysis at 10‐km grid resolution; stochastically generated data from weather generator; bias‐corrected dynamically downscaled; and bias‐corrected global reanalysis. The reanalysis products are considered as surrogates for large‐scale observations. Our study indicates that over the 28 watersheds in the southeastern United States, the simulated hydrological response to the bias‐corrected dynamically downscaled data is superior to the other four meteorological datasets. In comparison with synthetically generated meteorological forcing (from weather generator), the dynamically downscaled data from global atmospheric reanalysis result in more realistic hydrological simulations. Therefore, we conclude that dynamical downscaling of global reanalysis, which offers data for sufficient number of years (in this case 22 years), although resource intensive, is relatively more useful than other sources of meteorological data with comparable period in simulating realistic hydrological response at watershed scales. Copyright © 2013 John Wiley & Sons, Ltd.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

δ53Cr values of catchment runoff exhibit seasonality regardless of bedrock Cr concentrations: Role of periodically changing chromium export fluxes

One pre-requisite for the construction of a global chromium isotope mass balance is detailed understanding of Cr isotope systematics in the critical zone where redox-processes can modify the isotope signature of geogenic Cr input into the hydrosphere. A Cr isotope inventory of bedrock, soil, and runoff was performed in a Central European headwater catchment underlain by amphibolite, situated in the vicinity of two previously studied catchments underlain by different bedrock types (serpentinite and leucogranite). Fresh bedrock in the amphibolite catchment NAZ contained ~300 mg/kg Cr, serpentinite at PLB contained ~800 mg/kg Cr, and leucogranite at LYS contained ~2 mg/kg Cr. Monthly hydrochemical monitoring at all three sites revealed higher Cr(VI) export fluxes in winter than in summer. NAZ was characterized by a distinct seasonality in the δ⁵³Cr values, with minima during winter/spring snowmelts (−0.35‰) and maxima during dry summers (0.40‰). Similar seasonality in δ⁵³Cr values had been reported from PLB and LYS. Bedrock at all three sites had similar Cr isotope composition close to −0.10‰, a value indistinguishable from the δ⁵³Cr value of bulk silicate Earth (BSE). Positive mean δ⁵³Cr value of NAZ runoff indicated Cr-isotope fractionations during weathering of geogenic Cr(III), combined with adsorption of the resulting Cr(VI) on soil particles during pedogenesis. However, the mass-weighted mean δ⁵³Cr of NAZ runoff was lower (−0.08‰), indistinguishable from the Cr isotope signature of bedrock. The same pattern of lower mass-weighted mean δ⁵³Cr values of runoff, compared to arithmetic mean δ⁵³Cr values of runoff, were observed also at PLB and LYS. We suggest that elevated Cr runoff fluxes in winter remove some of the residual isotopically light Cr that accumulated in the soil during summer. Seasonality in runoff δ⁵³Cr values appears to be a relatively widespread phenomenon, de-coupled from Cr availability for chemical weathering.     

Urban stream temperature patterns: Spatial and temporal heterogeneity in the Philadelphia region,Pennsylvania, USA

Stream temperature is a critical water quality parameter that is not fully understood, particularly in urban areas. This study explores drivers contributing to stream temperature variability within an urban system, at 21 sites within the Philadelphia region, Pennsylvania, USA. A comprehensive set of temperature metrics were evaluated, including temperature sensitivity, daily maximum temperatures, time >20°C, and temperature surges during storms. Wastewater treatment plants (WWTPs) were the strongest driver of downstream temperature variability along 32 km in Wissahickon Creek. WWTP effluent temperature controlled local (1–3 km downstream) temperatures year-round, but the impacts varied seasonally: during winter, local warming of 2–7°C was consistently observed, while local cooling up to 1°C occurred during summer. Summer cooling and winter warming were detected up to 12 km downstream of a WWTP. Comparing effects from different WWTPs provided guidelines for mitigating their thermal impact; WWTPs that discharged into larger streams, had cooler effluent, or had lower discharge had less effect on stream temperatures. Comparing thermal regimes in four urban headwater streams, sites with more local riparian canopy had cooler maximum temperatures by up to 1.5°C, had lower temperature sensitivity, and spent less time at high temperatures, although mean temperatures were unaffected. Watershed-scale impervious area was associated with increased surge frequency and magnitude at headwater sites, but most storms did not result in a surge and most surges had a low magnitude. These results suggest that maintaining or restoring riparian canopy in urban settings will have a larger impact on stream temperatures than stormwater management that treats impervious area. Mitigation efforts may be most impactful at urban headwater sites, which are particularly vulnerable to stream temperature disruptions. It is vital that stream temperature impacts are considered when planning stormwater management or stream restoration projects, and the appropriate metrics need to be considered when assessing impacts.