Micrometeorological processes driving snow ablation in an Alpine catchment

Mountain snow covers typically become patchy over the course of a melting season. The snow pattern during melt is mainly governed by the end of winter snow depth distribution and the local energy balance. The objective of this study is to investigate micro-meteorological processes driving snow ablation in an Alpine catchment. For this purpose we combine a meteorological boundary-layer model (Advanced Regional Prediction System) with a fully distributed energy balance model (Alpine3D). Turbulent fluxes above melting snow are further investigated by using data from eddy-correlation systems. We compare modeled snow ablation to measured ablation rates as obtained from a series of Terrestrial Laser Scanning campaigns covering a complete ablation season. The measured ablation rates indicate that the advection of sensible heat causes locally increased ablation rates at the upwind edges of the snow patches. The effect, however, appears to be active over rather short distances of about 4–6 m. Measurements suggest that mean wind velocities of about 5 m s−1 are required for advective heat transport to increase snow ablation over a long fetch distance of about 20 m. Neglecting this effect, the model is able to capture the mean ablation rates for early ablation periods but strongly overestimates snow ablation once the fraction of snow coverage is below a critical value of approximately 0.6. While radiation dominates snow ablation early in the season, the turbulent flux contribution becomes important late in the season. Simulation results indicate that the air temperatures appear to overestimate the local air temperature above snow patches once the snow coverage is low. Measured turbulent fluxes support these findings by suggesting a stable internal boundary layer close to the snow surface causing a strong decrease of the sensible heat flux towards the snow cover. Thus, the existence of a stable internal boundary layer above a patchy snow cover exerts a dominant control on the timing and magnitude of snow ablation for patchy snow covers.

Atmospheric Flow Development and Associated Changes in Turbulent Sensible Heat Flux over a Patchy Mountain Snow Cover

Michael Lehning

In this study, the small-scale boundary layer dynamics and the energy balance over a fractional snow cover are numerically investigated. The atmospheric boundary layer flows over a patchy snow cover were calculated with an atmospheric model (Advanced Regional Prediction System) on a very high spatial resolution of 5 m. The numerical results revealed that the development of local flow patterns and the relative importance of boundary layer processes depend on the snow patch size distribution and the synoptic wind forcing. Energy balance calculations for quiescent wind situations demonstrated that well-developed katabatic winds exerted a major control on the energy balance over the patchy snow cover, leading to a maximum in the mean downward sensible heat flux over snow for high snow-cover fractions. This implies that if katabatic winds develop, total melt of snow patches may decrease for low snow-cover fractions despite an increasing ambient air temperature, which would not be predicted by most hydrological models. In contrast, stronger synoptic winds increased the effect of heat advection on the catchment's melt behavior by enhancing the mean sensible heat flux over snow for lower snow-cover fractions. A sensitivity analysis to grid resolution suggested that the grid size is a critical factor for modeling the energy balance of a patchy snow cover. The comparison of simulation results from coarse (50 m) and fine (5 m) horizontal resolutions revealed a difference in the spatially averaged turbulent heat flux over snow of 40%-70% for synoptic cases and 95% for quiescent cases.

Amer Meteorological Soc2015

How do patchy snow covers affect turbulent sensible heat fluxes? – Numerical analysis and experimental findings

Michael Lehning, Sebastian Schlögl

The surface energy balance of a snow cover significantly changes once the snow cover gets patchy. The substantial progress in knowledge about the surface energy balance of patchy snow covers is a mandatory requirement to reduce biases in flux parameterizations in larger scale meteorological or climatological models. The aim of this project was to numerically improve energy balance calculations late in the melting season when the spatial variability of turbulent fluxes is especially high owing to the complex feedback between bare/snow-covered areas and the atmosphere above. In order to account for the feedback between the atmosphere and the patchy snow-cover we calculated three-dimensional air temperature and wind velocity fields with the non-hydrostatic atmospheric model ARPS for an idealized flat test site initialized with different snow distributions and atmospheric conditions. The physics-based surface process model Alpine3D has been forced with these atmospheric fields close to the snow surface in order to resolve the small-scale spatial variability. We further initialized the model with atmospheric fields above the blending height as a reference case. The numerical analysis shows that for simulations initialized with fully-resolved atmospheric fields below the blending height, turbulent sensible heat fluxes are up to 50 W m-2 larger than for calculations forced without resolved atmospheric fields. This difference in turbulent sensible heat fluxes over snow increase with increasing number of snow patches and decreasing snow-cover fraction. This is mainly attributed to an increase in the mean near-surface air temperature over snow due to horizontal and vertical exchange processes induced by the heterogeneous land-surface. We used flux footprint estimations to analyse turbulence data measured during three ablation periods in the Dischma valley (Switzerland). This fundamental theory was deployed for eddy-covariance measurements revealing the origin of the measured turbulence as a function of the measurement height, atmospheric stability and wind speed. For patchy snow covers with flux footprints larger than the fetch distance, the observed turbulence is typically recorded as a combination of the turbulence over snow and the turbulence over the adjacent bare ground. These conditions were observed in the upper Dischma valley at the end of the melting period for a small snow cover fraction (< 30 %). The strong influence of the bare ground on the heat exchange over snow leads to much smaller measured heat fluxes than calculated in the physics-based one-dimensional snow model SNOWPACK. We show that at least two eddy-covariance measurements in different vertical heights and a snow cover mask to detect the fetch distance are required to clearly separate the turbulence over snow and the turbulence over the adjacent bare ground. Using the footprint technique turbulent sensible heat fluxes towards the snow patch significantly increase by this separation. This approach will allow an improved parametrization of turbulent heat fluxes over patchy snow-covers in larger-scale energy balance models.

2017

How Are Turbulent Sensible Heat Fluxes and Snow Melt Rates Affected by a Changing Snow Cover Fraction?

Michael Lehning

The complex interaction between the atmospheric boundary layer and the heterogeneous land surface is typically not resolved in numerical models approximating the turbulent heat exchange processes. In this study, we consider the effect of the land surface heterogeneity on the spatial variability of near-surface air temperature fields and on snow melt processes. For this purpose we calculated turbulent sensible heat fluxes and daily snow depth depletion rates with the physics-based surface process model Alpine3D. To account for the effect of a heterogeneous land surface (such as patchy snow covers) on the local energy balance over snow, Alpine3D is driven by two-dimensional atmospheric fields of air temperature and wind velocity, generated with the non-hydrostatic atmospheric model Advanced Regional Prediction System. The atmospheric model is initialized with a set of snow distributions [snow cover fraction (SCF) and number of snow patches] and atmospheric conditions (wind velocities) for an idealized flat test site. Numerical results show that the feedback of the heterogeneity of the land surface (snow, no snow) on the near-surface variability of the atmospheric fields result in a significant increase in the mean air temperature Delta T-a = 1.8 K (3.7 and 4.9 K) as the SCF is decreased from a continuous snow cover to 55% (25 and 5%). Mean air temperatures over snow heavily increase with increasing initial wind velocities and weakly increase with an increasing number of snow patches. Surface turbulent sensible heat fluxes and daily snow depth depletion rates are strongly correlated to mean air temperatures, leading to 22-40% larger daily snow depth depletion rates for patchy snow covers. Numerical results from the idealized test site are compared with a test site in complex terrain. As slope-induced atmospheric processes (such as the development of katabatic flows) modify turbulent sensible heat fluxes, the variation of the surface energy balance is larger in complex terrain than for an idealized flat test site.

FRONTIERS MEDIA SA2018

How do patchy snow covers affect turbulent sensible heat fluxes? – Numerical analysis and experimental findings

Michael Lehning, Sebastian Schlögl

2017

How Are Turbulent Sensible Heat Fluxes and Snow Melt Rates Affected by a Changing Snow Cover Fraction?

Michael Lehning

FRONTIERS MEDIA SA2018