Boundary layer processes over a patchy mountain snow cover and associated changes in the energy balance

Heterogeneous land-covers induce the formation of secondary flows and small-scale boundary layer processes. While turbulent fluxes over a continuous snow cover are mostly governed by the spatial variability of wind velocity and solar radiation, local advection of sensible heat and boundary layer decoupling become additionally important for the surface energy exchange over patchy snow covers. Furthermore, evolving katabatic and anabatic winds strongly change the mean flow patterns over snow free and snow covered areas causing positive and negative buoyancy production. In this study we use an atmospheric model (Advanced Regional Prediction System) to investigate boundary layer processes over patchy snow covers in an alpine catchment. The model is initialized with low to high snow cover ratios based on field measurements of snow distribution. We particularly analyse the effect of heat advection, boundary layer decoupling and changing patterns of secondary flows on the energy balance at the local and catchment scale. The numerical results show that the changing land-cover controls the length scale of internal boundary layers, while secondary flows drive the strength of boundary layer decoupling and heat advection. The relative importance of boundary layer processes depend on the snow patch size distribution and the synoptic wind. Stronger synoptic winds and low snow cover ratios increase the effect of heat advection and decrease the impact of boundary layer decoupling on the catchment`s melt behaviour. The typical lengths of stable internal boundary layers coincide with snow patch size distribution, which is consistent with typical length scales of snow accumulation patterns. As those length scales are in the order of tens of meters for the studied alpine catchment, very small grid sizes (below 10 m) are required to adequately model boundary layer processes which are shown to significantly alter the energy balance over 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

Large-Eddy Simulation of the Convective Atmospheric Boundary Layer over Heterogeneous Land Surfaces

Chad Higgins, Daniel Nadeau, Eric Richard Pardyjak, Marc Parlange

The parameterization of the atmospheric boundary layer is essential for accurate numerical weather predictions. The near-surface values of air temperature or wind speed for instance are highly dependent on the complex land – atmosphere interactions over heterogeneous terrain. Over such surfaces, several open challenges remain regarding the growth of internal boundary layers, the determination of mixing layer heights and the spatial distribution of heat and momentum fluxes. The large-eddy simulation (LES) code we apply is based on the robust Lagrangian scale-dependent dynamic subgrid-scale model (Bou-Zeid et al., 2005, Physics of Fluids). The flow is driven by a mean pressure gradient expressed in terms of horizontal geostrophic winds. The code is pseudo-spectral with spectral decomposition in the horizontal dimensions. Land surface heterogeneities take the form of spatially distributed patches of surface temperature and roughness derived from satellite imagery and land use analysis. The simulation of the complete diurnal cycle is possible but we restrict our discussion to fully convective conditions. Data from LITFASS – 2003 (Lindenberg Inhomogeneous Terrain Fluxes between Atmosphere and Surface: a long-term Study) are used to validate the simulation results. The LES domain covers a 99-m meteorological tower with turbulent measurements and five surface micrometeorological stations. We investigate the presence of large-scale turbulent structures that are typical for daytime conditions with a particular interest in the blending height. Overall, this study illustrates how important is the accurate parameterization of heterogeneous terrain in studies of the atmospheric boundary layer.

2009

Large-Eddy Simulation of the Convective Atmospheric Boundary Layer over Heterogeneous Land Surfaces

Chad Higgins, Daniel Nadeau, Marc Parlange

The parameterization of the atmospheric boundary layer is crucial for accurate numerical weather predictions. Over heterogeneous terrain, several open challenges remain regarding the growth of internal boundary layers, the determination of mixing layer heights and the spatial distribution of heat and momentum fluxes. The large-eddy simulation (LES) code we apply is based on the very robust Lagrangian scale-dependent dynamic subgrid-scale model. The flow is driven by a mean pressure gradient expressed in terms of horizontal geostrophic winds. The code is pseudo-spectral with spectral decomposition in the horizontal dimensions. Land surface heterogeneities take the form of distributed patches of surface temperature and momentum roughness derived from satellite imagery and land use analysis. Although recent developments in the code have enabled the simulation of the complete diurnal cycle, the focus here is restricted to convective conditions. Data from LITFASS 2003 (standing for Lindenberg Inhomogeneous Terrain Fluxes between Atmosphere and Surface: a long-term Study) are utilized to validate the simulation results. The LES domain covers a 99-m meteorological tower with turbulent measurements and a few surface micrometeorological stations. We investigate the presence of large-scale turbulent structures that are typical for daytime conditions with a particular interest in the blending height. Overall, this study illustrates how important is the accurate parameterization of heterogeneous terrain in studies of the atmospheric boundary layer.

2009