Snowfall Limit Forecasts and Hydrological Modeling

Hydrological flood forecasting in mountainous areas requires accurate partitioning between rain and snowfall to properly estimate the extent of runoff contributing areas. Here a method to make use of snowfall limit information-a standard output of limited-area models (LAMs)-for catchment-scale hydrological modeling is proposed. LAMs consider the vertical, humid, atmospheric structure in their snowfall limit calculations. The proposed approach is thus more physically based than inferring snowfall limit estimates based on (dry) ground temperature measurements, which is the standard procedure in most hydrological models. The presented case study uses forecast reanalyses from the Consortium for Small-Scale Modeling (COSMO) limited-area model as input for discharge simulation in a topographically complex catchment in the Swiss Alps. Results suggest that the use of COSMO snowfall limits during spring snowmelt periods can provide more accurate runoff simulations than routine procedures, with practical implications for operational hydrology in Alpine regions.

Improving Alpine Flood Prediction through Hydrological Process Characterization and Uncertainty Analysis

Cara Christine Tobin

Among the many challenges of Alpine flood prediction is describing complex, meteo-hydrological processes in a simplified, robust manner that can be easily integrated into operational forecasting. In this dissertation, improved methods to characterize these processes are developed and integrated into the hydrological modeling component of an operational flood forecasting system used in the Swiss Alps. Detailed studies are conducted to improve hydrological model inputs, processes and outputs. Improvements, detailed in four chapters of this thesis, address the overarching goal of this work – the reduction of flood forecasting uncertainty. The accuracy of flood predictions in Alpine areas is contingent upon adequate interpolation of meteorological forcings, which has significant impacts on discharge volumes and flood peaks. This thesis demonstrates an improvement in the interpolation of temperature and precipitation inputs using a robust variogram which considers anisotropy and using a geostatistical interpolation method to distribute inputs in space and time. Results show that using elevation as the external drift factor better describes orographically-induced precipitation and temperature gradients. Also, the consideration of anisotropy is integral in detailing spatial patterns of precipitation induced by storm advection. Hydrological flood forecasting in mountainous areas also requires accurate partitioning between rain and snowfall to properly estimate the extent of runoff contributing areas. Partitioning is improved in this work by using a new method to integrate Limited Area Model output. Unlike standard hydrological procedures in inferring snowfall limit estimates based on dry, ground temperature measurements, Limited Area Model output considers the vertical, humid, atmospheric structure in its snowfall limit calculations. In effect, this method provides good estimates of runoff contributing areas in the spring as evidenced by validation on discharge measurements and satellite images of snow coverage. Accurately describing snowmelt processes on a sub-daily scale is also of critical importance in Alpine flood forecasting. However, the complex topography of the study region has limited observations available for validation. This thesis presents the development of a new physically-based snowmelt method applicable to regions with limited data. This method uses only daily minimum and maximum temperatures to mimic the effects of radiation. A comparative analysis of snowmelt methods is validated with snow lysimeter data and with a unique, distributed meteorological dataset collected by a wireless weather station network. Results demonstrate that the new method is competitive with more complex snowmelt methods as shown by accurately reproducing diurnal snowmelt cycles. Conveying limits of certainty on flood prediction outputs to users is critical because of epistemic and aleatory errors inherent to environmental modeling. Due to the presence of these errors, the GLUE methodology and multi-criteria performance ideas have been adapted with a fit-for-purpose uncertainty estimation technique in the final part of this thesis. With this method, hydrological model parameters are constrained based on hydrograph behavior, with a particular focus on flood peak response. A key component of the technique is a visualization tool which shows acceptable ensembles of discharge with respect to individual and combined criteria. By integrating the aforementioned input and process improvements into the hydrological model, calibration achieves model outputs that capture observed river discharge. Also, the uncertainty associated with hydrological modeling output error is reduced. Findings of this thesis are applicable to operational flood forecasting in general and have proven utility in improving hydrological model predictions in mountainous regions. Due to the novelty of the developments in terms of new methods or the use of tools and data sources previously unexploited in flood forecasting, further testing of the improvements is recommended. Future research in quantifying the chain of uncertainty produced by combining probabilistic forecast inputs with the hydrological output ensembles is also critical when the improved flood forecasting model becomes fully operational.

EPFL2012

Modelling Drifting Snow Sublimation and Investigating Small-Scale Snow Transport

Christine Dorothea Groot Zwaaftink

Drifting snow has a large influence on the mountain snow cover and the overlying atmospheric boundary layer. Consequently, drifting snow also affects avalanches, snow hydrology, vegetation and climate. Considerable effort is devoted to understanding the process in order to model the effects on the mountain snow cover or the overlying atmosphere. This thesis work aims to I) provide estimates of drifting snow sublimation within a small catchment in the Swiss Alps and II) test an alternative method to describe snow transport on short time and small spatial scales. Drifting snow sublimation moistens and cools the surrounding air and causes a loss of snow mass. Based on the feedback mechanisms integrating air temperature, humidity and snow concentration, it should be seen as a self-limiting process. Drifting snow sublimation (including these feedback mechanisms) has been implemented in Alpine3D, which is a high-resolution snow transport model for alpine terrain. This model has been applied to a small area in the Swiss Alps for a case study of approximately two days. Results show that, in general, the negative feedbacks of sublimation on the snow mass concentration, air temperature, and humidity are small. They are, however, relevant on the slope scale, as could be retrieved from analyzing the deposition on a leeward slope. Compared to a reference simulation without sublimation, drifting snow sublimation reduced the deposition by approximately 12% on this leeward slope when omitting all feedback mechanisms. In a simulation where all feedbacks were included, drifting snow sublimation reduced the snow deposition on this slope by 10%. Drifting snow sublimation thus appears to be a significant process for a leeward slope influenced by drifting snow. On the contrary, the spatially averaged reduction of deposition due to drifting snow sublimation in the complete modelled domain was 2.3%. During this case study, this value was comparable to surface sublimation. At the same field site a simulation of a complete season was performed and in this case, sublimation of drifting snow was accounted for only in the absence of concurrent precipitation. As for the case study, the slope affected the most by drifting snow sublimation was shown to be a leeward slope for SE wind. This could be explained by generally warmer and dryer conditions during periods with SE wind. In this particular slope, the snow amount was reduced by about 1.8% due to drifting snow sublimation. For the complete domain however, the model results indicate that drifting snow sublimation is much smaller than surface sublimation and negligible for the total water equivalent. Thus, at the studied field site drifting snow sublimation seems only significant locally or on short time scales. Limitations of drifting snow simulations based on mean wind fields are seen as the variability of the snow cover is underestimated on scales of few meters. Moreover, these models are also not able to represent snow transport on very short time scales. High-frequency observations of drifting snow and sand saltation have shown that turbulent structures have an influence on the snow or sand mass flux. Improvements in modelling drifting snow on such time scales can thus be reached by considering turbulence. Here, a combination of wind fields from large eddy simulations and a Lagrangian stochastic model was used to describe snow particle trajectories over flat terrain. This model reproduced a time series of the snow mass flux that qualitatively corresponds to observations from the Weissfluhjoch Versuchsfeld. Furthermore, results show that the use of the surface shear stress rather than the friction velocity, as in conventional models, resulted in a spatial variation of drifting snow despite homogeneous snow properties and the flatness of the considered terrain. The distribution of the surface shear stress in space and time was found reasonable yet narrow based on the ratio between the standard deviation and the mean shear stress, which was 0.17 compared to a value of 0.4 reported for measurements. Nonetheless, this variable surface shear stress played a key role in the modelled spatial variation of drifting snow. This also points at a relevant influence of the particle aerodynamic entrainment process for drifting snow. The results thus show that this model seems to overcome limitations of conventional three-dimensional snow transport models as it gives a more realistic description of drifting snow on short time scales and it may be used to gain further insight in drifting snow.

EPFL2013

Flow and Temperature Dynamics in the Hydrologic Response of Alpine Catchments: Travel Time Formulation and Geomorphologic Signatures

Francesco Comola, Michael Lehning, Andrea Rinaldo, Bettina Schaefli

We present a travel time formulation of water and energy transport at sub-catchment scale. The derived equations are implemented in Alpine3D, a physically-based model of snow processes, which provides the necessary boundary conditions to perform hydro-thermal response simulations of Alpine catchments. The model set-up accounts for advective and non- advective energy fluxes to perform spatially distributed simulations of streamflow and temperature in river networks having an arbitrary degree of geomorphological complexity. The model gives reliable predictions of streamflow and tem- perature, as shown by comparing modeled and measured hydrographs and thermographs at the outlet of the Dischma catchment (45 km2) in the Swiss Alps. Our model setup is applied to investigate the role of hillslope aspects, representing the main control on radiation and snowmelt patterns, in the flow regime of the study catchment. The distributed simu- lation results show that snowmelt-induced discharge exhibits a visible geomorphologic signature of aspects at sub-catch- ment scale, but this progressively fades out going from headwater streams to the outlet. Accordingly, the geomorphologic signature is scale-dependent: it is significant at small scales where the high aspect correlation generates predominant orientations but is lost at larger scales where aspects are de-correlated and different orientations are averaged out. We further apply the model to investigate the geomorphologic signature of drainage density in the thermal regime of the study catchment. The results show that the contribution of the advective energy fluxes becomes progressively smaller when the drainage density increases, while the one of the non-advective energy fluxes becomes larger. Moreover, such variations balance out at the catchment outlet, where the temperature signal is not sensitive to the increasing drainage density. The relevance of the performed invetigations stems from the increasing scientific interest concerning the impacts of the warming climate on water resources management and temperature-influenced ecological processes.

2014