• Canada
• Publicationsclose
• Other research productsclose
• Open Access close
• CL close
• English close
• Hydrology and Earth System Sciences (HESS) close

Information about the raindrop size distribution (RSD) is vital for comprehending the precipitation microphysics, improving the rainfall estimation algorithms, and appraising the rainfall erosivity. Previous research has revealed that the RSD exhibits diversity with geographical location and weather type, which leads to the assessment of the region and weather-specific RSDs. Based on long-term (2004 to 2016) disdrometer measurements in northern Taiwan, this study attempts to demonstrate the RSD aspects of summer seasons that were bifurcated into two weather conditions, namely typhoon (TY) and non-typhoon (NTY) rainfall. The results show a higher concentration of small drops and a lower concentration of large-sized drops in TY compared to NTY rainfall, and this behavior persisted even after characterizing the RSDs into different rainfall rate classes. RSDs expressed in gamma parameters show higher mass-weighted mean diameter (Dm) and lower normalized intercept parameter (Nw) values in NTY than TY rainfall. Moreover, sorting these two weather conditions (TY and NTY rainfall) into stratiform and convective regimes revealed a larger Dm in NTY than in TY rainfall. The RSD empirical relations used in the valuation of rainfall rate (Z–R, Dm–R, and Nw–R) and rainfall kinetic energy (KE–R and KE–Dm) were enumerated for TY and NTY rainfall, and they exhibited profound diversity between these two weather conditions. Attributions of RSD variability between the TY and NTY rainfall to the thermodynamical and microphysical processes are elucidated with the aid of reanalysis, remote sensing, and ground-based data sets.

Abstract. Thresholds and hydrologic connectivity associated with runoff processes are a critical concept for understanding catchment hydrologic response at the event timescale. To date, most attention has focused on single runoff response types, and the role of multiple thresholds and flow path connectivities has not been made explicit. Here we first summarise existing knowledge on the interplay between thresholds, connectivity and runoff processes at the hillslope–small catchment scale into a single figure and use it in examining how runoff response and the catchment threshold response to rainfall affect a suite of runoff generation mechanisms in a small agricultural catchment. A 1.37 ha catchment in the Lang Lang River catchment, Victoria, Australia, was instrumented and hourly data of rainfall, runoff, shallow groundwater level and isotope water samples were collected. The rainfall, runoff and antecedent soil moisture data together with water levels at several shallow piezometers are used to identify runoff processes in the study site. We use isotope and major ion results to further support the findings of the hydrometric data. We analyse 60 rainfall events that produced 38 runoff events over two runoff seasons. Our results show that the catchment hydrologic response was typically controlled by the Antecedent Soil Moisture Index and rainfall characteristics. There was a strong seasonal effect in the antecedent moisture conditions that led to marked seasonal-scale changes in runoff response. Analysis of shallow well data revealed that streamflows early in the runoff season were dominated primarily by saturation excess overland flow from the riparian area. As the runoff season progressed, the catchment soil water storage increased and the hillslopes connected to the riparian area. The hillslopes transferred a significant amount of water to the riparian zone during and following events. Then, during a particularly wet period, this connectivity to the riparian zone, and ultimately to the stream, persisted between events for a period of 1 month. These findings are supported by isotope results which showed the dominance of pre-event water, together with significant contributions of event water early (rising limb and peak) in the event hydrograph. Based on a combination of various hydrometric analyses and some isotope and major ion data, we conclude that event runoff at this site is typically a combination of subsurface event flow and saturation excess overland flow. However, during high intensity rainfall events, flashy catchment flow was observed even though the soil moisture threshold for activation of subsurface flow was not exceeded. We hypothesise that this was due to the activation of infiltration excess overland flow and/or fast lateral flow through preferential pathways on the hillslope and saturation overland flow from the riparian zone.

In practice, the catchment hydrologist is often confronted with the task of predicting discharge without having the needed records for calibration. Here, we report the discharge predictions of 10 modellers – using the model of their choice – for the man-made Chicken Creek catchment (6 ha, northeast Germany, Gerwin et al., 2009b) and we analyse how well they improved their prediction in three steps based on adding information prior to each following step. The modellers predicted the catchment's hydrological response in its initial phase without having access to the observed records. They used conceptually different physically based models and their modelling experience differed largely. Hence, they encountered two problems: (i) to simulate discharge for an ungauged catchment and (ii) using models that were developed for catchments, which are not in a state of landscape transformation. The prediction exercise was organized in three steps: (1) for the first prediction the modellers received a basic data set describing the catchment to a degree somewhat more complete than usually available for a priori predictions of ungauged catchments; they did not obtain information on stream flow, soil moisture, nor groundwater response and had therefore to guess the initial conditions; (2) before the second prediction they inspected the catchment on-site and discussed their first prediction attempt; (3) for their third prediction they were offered additional data by charging them pro forma with the costs for obtaining this additional information. Holländer et al. (2009) discussed the range of predictions obtained in step (1). Here, we detail the modeller's assumptions and decisions in accounting for the various processes. We document the prediction progress as well as the learning process resulting from the availability of added information. For the second and third steps, the progress in prediction quality is evaluated in relation to individual modelling experience and costs of added information. In this qualitative analysis of a statistically small number of predictions we learned (i) that soft information such as the modeller's system understanding is as important as the model itself (hard information), (ii) that the sequence of modelling steps matters (field visit, interactions between differently experienced experts, choice of model, selection of available data, and methods for parameter guessing), and (iii) that added process understanding can be as efficient as adding data for improving parameters needed to satisfy model requirements. Hydrology and Earth System Sciences, 18 (6) ISSN:1027-5606 ISSN:1607-7938

Assessing the impacts of both natural (e.g., tidal forcing from the ocean) and human-induced changes (e.g., dredging for navigation, land reclamation) on estuarine morphology is particularly important for the protection and management of the estuarine environment. In this study, a novel analytical approach is proposed for the assessment of estuarine morphological evolution in terms of tidally averaged depth on the basis of the observed water levels along the estuary. The key lies in deriving a relationship between wave celerity and tidal damping or amplification. For given observed water levels at two gauging stations, it is possible to have a first estimation of both wave celerity (distance divided by tidal travelling time) and tidal damping or amplification rate (tidal range difference divided by distance), which can then be used to predict the morphological changes via an inverse analytical model for tidal hydrodynamics. The proposed method is applied to the Lingdingyang Bay of the Pearl River Estuary, located on the southern coast of China, to analyse the historical development of the tidal hydrodynamics and morphological evolution. The analytical results show surprisingly good correspondence with observed water depth and volume in this system. The merit of the proposed method is that it provides a simple approach for understanding the decadal evolution of the estuarine morphology through the use of observed water levels, which are usually available and can be easily measured. National Key R&D of China (Grant No. 2016YFC0402601), National Natural Science Foundation of China (Grant No. 51979296, 51709287, 41706088, 41476073), Fundamental Research Funds for the Central Universities (No.18lgpy29) and from the Water Resource Science and Technology Innovation Program of Guangdong Province (Grant No. 2016-20, 2016-21). The work of the second author was supported by FCT research contracts IF/00661/2014/CP1234. info:eu-repo/semantics/submittedVersion

The phase of precipitation and its distribution at the surface can affect water resources and the regional water cycle of a region. A field project was held in March–April 2015 on the eastern slope of the Canadian Rockies to document precipitation characteristics and associated atmospheric conditions. During the project, 60 % of the particles documented were rimed in relatively warm and dry conditions. Rain–snow transitions also occurred aloft and at the surface in sub-saturated conditions. Ice-phase precipitation falling through a saturated atmospheric layer with temperatures > 0 ∘C will start melting. In contrast, if the melting layer is sub-saturated, the ice-phase precipitation undergoes sublimation, which increases the depth of the rain–snow transition. In this context, this study investigates the role of sublimation and riming in precipitation intensity and type reaching the surface in the Kananaskis Valley, Alberta, during March–April 2015. To address this, a set of numerical simulations of an event of mixed precipitation observed at the surface was conducted. This event on 31 March 2015 was documented with a set of devices at the main observation site (Kananaskis Emergency Services, KES), including a precipitation gauge, disdrometer, and micro rain radar. Sensitivity experiments were performed to assess the impacts of temperature changes from sublimation and the role of the production of graupel (riming) aloft in the surface precipitation evolution. A warmer environment associated with no temperature changes from sublimation leads to a peak in the intensity of graupel at the surface. When the formation of graupel is not considered, the maximum snowfall rate occurred at later times. Results suggest that unrimed snow reaching the surface is formed on the western flank and is advected eastward. In contrast, graupel would form aloft in the Kananaskis Valley. The cooling from sublimation and melting by rimed particles increases the vertical shear near KES. Overall, this study illustrated that the presence of graupel influenced the surface evolution of precipitation type in the valley due to the horizontal transport of precipitation particles.

Understanding hydrological responses to reforestation is an important subject in watershed management, particularly in large forested watersheds ( > 1000 km2). In this study, we selected two large forested watersheds (Pingjiang and Xiangshui) located in the upper reach of the Poyang Lake watershed, southeastern China (with an area of 3261.4 and 1458 km2, respectively), along with long-term data on climate and hydrology (1954–2006) to assess the effects of large-scale reforestation on streamflow. Both watersheds have similar climate and experienced comparable and dramatic forest changes during the past decades, but with different watershed properties (e.g., the topography is much steeper in Xiangshui than in Pingjiang), which provides us with a unique opportunity to compare the differences in hydrological recovery in two contrasted watersheds. Streamflow at different percentiles (e.g., 5, 10, 50 and 95 %) were compared using a combination of statistical analysis with a year-wise method for each watershed. The results showed that forest recovery had no significant effects on median flows (Q50%) in both watersheds. However, reforestation significantly reduced high flows in Pingjiang, but had limited influence in Xiangshui. Similarly, reforestation had significant and positive effects on low flows (Q95%) in Pingjiang, while it did not significantly change low flows in Xiangshui. Thus, hydrological recovery is limited and slower in the steeper Xiangshui watershed, highlighting that watershed properties are also important for determining hydrological responses to reforestation. This finding has important implications for designing reforestation and watershed management strategies in the context of hydrological recovery.

Abstract. Measurements of snowmelt and turbulent heat fluxes were made during the snowmelt periods of two years at two neighbouring tundra sites in the Yukon, one in a sheltered location with tall shrubs exposed above deep snow and the other in an exposed location with dwarf shrubs covered by shallow snow. The snow was about twice as deep in the valley as on the plateau at the end of each winter and melted out about 10 days later. The site with buried vegetation showed a transition from air-to-surface heat transfers to surface-to-air heat transfers as bare ground became exposed during snowmelt, but there were daytime transfers of heat from the surface to the air at the site with exposed vegetation even while snow remained on the ground. A model calculating separate energy balances for snow and exposed vegetation, driven with meteorological data from the sites, is found to be able to reproduce these behaviours. Averaged over 30-day periods the model gives about 8 Wm−2 more sensible heat flux to the atmosphere for the valley site than for the plateau site. Sensitivity of simulated fluxes to model parameters describing vegetation cover and density is investigated.

Permafrost strongly controls hydrological processes in cold regions. Our understanding of how changes in seasonal and perennial frozen ground disposition and linked storage dynamics affect runoff generation processes remains limited. Storage dynamics and water redistribution are influenced by the seasonal variability and spatial heterogeneity of frozen ground, snow accumulation and melt. Stable isotopes are potentially useful for quantifying the dynamics of water sources, flow paths and ages, yet few studies have employed isotope data in permafrost-influenced catchments. Here, we applied the conceptual model STARR (the Spatially distributed Tracer-Aided Rainfall–Runoff model), which facilitates fully distributed simulations of hydrological storage dynamics and runoff processes, isotopic composition and water ages. We adapted this model for a subarctic catchment in Yukon Territory, Canada, with a time-variable implementation of field capacity to include the influence of thaw dynamics. A multi-criteria calibration based on stream flow, snow water equivalent and isotopes was applied to 3 years of data. The integration of isotope data in the spatially distributed model provided the basis for quantifying spatio-temporal dynamics of water storage and ages, emphasizing the importance of thaw layer dynamics in mixing and damping the melt signal. By using the model conceptualization of spatially and temporally variable storage, this study demonstrates the ability of tracer-aided modelling to capture thaw layer dynamics that cause mixing and damping of the isotopic melt signal.

Abstract. Understanding and modeling snowmelt-runoff generation in seasonally-frozen regions is a major challenge in hydrology. Partly, this is because the controls on hillslope-scale snowmelt-runoff generation are potentially extensive and their hierarchy is poorly understood. Understanding the relative importance of controls (e.g. topography, vegetation, land use, soil characteristics, and precipitation dynamics) on runoff response is necessary for model development, spatial extrapolation, and runoff classification schemes. Multiple interacting process controls, the nonlinearities between them, and the resultant threshold-like activation of runoff, typically are not observable in short-term experiments or single-season field studies. Therefore, long-term datasets and analyses are needed. Here, we use a 52-year dataset of runoff, precipitation, soil water content, snow cover, and meteorological data from three monitored c.5 ha hillslopes on the Canadian Prairies to determine the controls on snowmelt-runoff, their time-varying hierarchy, and the interactions between the controls. We use decision tree learning to extract information from the dataset on the controls on runoff ratio. Our analysis shows that there was a variable relationship between total spring runoff amount and either winter snowfall amount or snow cover water equivalent. Other factors came into play to control the fraction of precipitated water that infiltrated into the frozen ground. In descending order of importance, these were: total snowfall, snow cover, fall soil surface water content, melt rate, melt season length, and fall soil profile water content. While mid-winter warm periods in some years likely increased soil water content and/or led to development of impermeable ice lenses that affected the runoff response, hillslope memory of fall soil moisture conditions played a strong role in the spring runoff response. The hierarchy of these controls was condition-dependent, with the biggest differences between high and low snow cover seasons, and wet and dry fall soil moisture conditions. For example, when snow cover was high, the top three controls on runoff ratio matched the overall hierarchy of controls, with fall soil surface water content being the most important of these. By comparison, when snow cover was low, fall soil surface content was relatively unimportant and superseded by four other controls. Existing empirical methods for predicting infiltration into frozen ground failed to adequately predict runoff response at our site. Our analysis of the hierarchy of controls on meltwater runoff will aid in focusing new model approaches and understanding what to focus future measurement campaigns on in snowmelt-dominated, seasonally-frozen regions.

The 0 °C temperature threshold is critical to many meteorological and hydrological processes driven by melting and freezing in the atmosphere, surface and sub-surface and by the associated precipitation varying between rain, freezing rain, wet snow and snow. This threshold, linked with freeze-thaw, is especially important in cold regions such as Canada. This study develops a Canada-wide perspective on near 0 °C conditions with a particular focus on the occurrence of its associated precipitation. Since this analysis requires hourly values of surface temperature and precipitation type observations, it was limited to 92 stations over the 1981–2011 period. In addition, nine stations representative of various climatic regions are selected for further analysis. Near 0 °C conditions are defined as periods when the surface temperature is between −2 °C and 2 °C. Near 0 °C conditions occur often across all regions of the country although the annual number of days and hours and the duration of these events varies dramatically. Various forms of precipitation (including rain, freezing rain, wet snow and ice pellets) are sometimes linked with these temperatures with highest fractions tending to occur in Atlantic Canada. Trends of most temperature-based and precipitation-based indicators show little or no change despite a systematic warming in annual temperatures. Over the annual cycle, near 0 °C temperatures and precipitation often exhibit a pattern with short durations near summer driven by the diurnal cycle, while longer durations tend to occur more towards winter associated with storms. There is also a tendency for near 0 °C temperatures to occur more often than expected relative to other temperature windows; due at least in part to diabatic cooling and heating occurring with melting and freezing, respectively, in the atmosphere and at the surface.