The theoretical treatment of turbulence is largely based on the assumption of horizontally homogeneous and flat underlying surfaces. Correspondingly, approaches developed over the years to measure ...turbulence statistics in order to test this theoretical understanding or to provide model input, are also largely based on the same assumption of horizontally homogeneous and flat terrain. Here we discuss aspects of turbulence measurements that require special attention in mountainous terrain. We especially emphasize the importance of data quality (flux corrections, data quality assessment, uncertainty estimates) and address the issues of coordinate systems and different post-processing options in mountainous terrain. The appropriate choice of post-processing methods is then tested based on local scaling arguments. We demonstrate that conclusions drawn from turbulence measurements obtained in mountainous terrain are rather sensitive to these post-processing choices and give suggestions as to those that are most appropriate.
Coupling of the earth’s surface with the atmosphere is achieved through an exchange of momentum, energy, and mass in the atmospheric boundary layer. In mountainous terrain, this exchange results from ...a combination of multiple transport processes, which act and interact on different spatial and temporal scales, including, for example, orographic gravity waves, thermally driven circulations, moist convection, and turbulent motions. Incorporating these exchange processes and previous studies, a new definition of the atmospheric boundary layer in mountainous terrain, a mountain boundary layer (MBL), is defined. This paper summarizes some of the major current challenges in measuring, understanding, and eventually parameterizing the relevant transport processes and the overall exchange between the MBL and the free atmosphere. Further details on many aspects of the exchange in the MBL are discussed in several other papers in this issue.
Eddy-covariance data from five stations in the Inn Valley, Austria, are analyzed for stable conditions to determine the gap scale that separates turbulent from large-scale, non-turbulent motions. The ...gap scale is identified from (co)spectra calculated from different variables using both Fourier analysis and multi-resolution flux decomposition. A correlation is found between the gap scale and the mean wind speed and stability parameter
z
/
L
that is used to determine a time-varying filter time, whose performance in separating turbulent and non-turbulent motions is compared to the performance of constant filter times between 0.5 and 30 min. The impact of applying different filter times on the turbulence statistics depends on the parameter and location, with a comparatively smaller impact on the variance of the vertical wind component than on the horizontal components and the turbulent fluxes. Results indicate that a time-varying filter time based on a multi-variable fit taking both mean wind speed and stability into account and a constant filter time of 2–3 min perform best in that they remove most of the non-turbulent motions while at the same time capturing most of the turbulence. For the studied sites and conditions, a time-varying filter time does not outperform a well chosen constant filter time because of relatively small variations in the filter time predicted by the correlation with mean flow parameters.
Global warming will lead to earlier beginnings and prolongation of growing seasons in temperate regions and will have pronounced effects on phenology and life-history adaptation in many species. ...These changes were not easy to simulate for actual phenologies because of the rudimentary temporal (season) and spatial (regional) resolution of climate model projections. We investigate the effect of climate change on the regional incidence of a pest insect with nearly worldwide distribution and very high potential for adaptation to season length and temperature--the Codling Moth, Cydia pomonella. Seasonal and regional climate change signals were downscaled to the hourly temporal scale of a pest phenology model and the spatial scale of pest habitats using a stochastic weather generator operating at daily scale in combination with a re-sampling approach for simulation of hourly weather data. Under future conditions of increased temperatures (2045-2074), the present risk of below 20% for a pronounced second generation (peak larval emergence) in Switzerland will increase to 70-100%. The risk of an additional third generation will increase from presently 0-2% to 100%. We identified a significant two-week shift to earlier dates in phenological stages, such as overwintering adult flight. The relative extent (magnitude) of first generation pupae and all later stages will significantly increase. The presence of first generation pupae and later stages will be prolonged. A significant decrease in the length of overlap of first and second generation larval emergence was identified. Such shifts in phenology may induce changes in life-history traits regulating the life cycle. An accordingly life-history adaptation in photoperiodic diapause induction to shorter day-length is expected and would thereby even more increase the risk of an additional generation. With respect to Codling Moth management, the shifts in phenology and voltinism projected here will require adaptations of plant protection strategies to maintain their sustainability.
“Valley-wind days” are characterized by synoptically undisturbed, clear-sky conditions, which lead to the formation of thermally-driven slope- and valley-wind circulations in mountain regions. A ...simple method is presented to identify these conditions in the Inn Valley, Austria, using ERA-Interim geopotential height fields and a clear-sky index, which is calculated from measurements of longwave incoming radiation, air temperature, and humidity at a single site on the valley floor. As the method is based on identifying weak synoptic-scale flows and clear skies, it can also be applied to the identification of ideal conditions for other thermally-driven circulations. The mean diurnal cycle of the valley-wind circulation on these days is briefly discussed.
The exchange of heat, momentum, and mass in the atmosphere over mountainous terrain is controlled by synoptic-scale dynamics, thermally driven mesoscale circulations, and turbulence. This article ...reviews the key challenges relevant to the understanding of exchange processes in the mountain boundary layer and outlines possible research priorities for the future. The review describes the limitations of the experimental study of turbulent exchange over complex terrain, the impact of slope and valley breezes on the structure of the convective boundary layer, and the role of intermittent mixing and wave–turbulence interaction in the stable boundary layer. The interplay between exchange processes at different spatial scales is discussed in depth, emphasizing the role of elevated and ground-based stable layers in controlling multi-scale interactions in the atmosphere over and near mountains. Implications of the current understanding of exchange processes over mountains towards the improvement of numerical weather prediction and climate models are discussed, considering in particular the representation of surface boundary conditions, the parameterization of sub-grid-scale exchange, and the development of stochastic perturbation schemes.
Assessing the avalanche danger level requires snow stratigraphy and instability data. As such data are usually sparse, we investigated whether distributed snow cover modeling can be used to provide ...information on spatial instability patterns relevant for regional avalanche forecasting. Using Alpine3D, we performed spatially distributed simulations to evaluate snow instability for the winter season 2016–17 in the region of Davos, Switzerland. Meteorological data from automatic weather stations were interpolated to 100 m horizontal resolution and precipitation was scaled with snow depth measurements from airborne laser scanning. Modeled snow instability metrics assessed for two different weak layers suggested that the weak layer closer to the snow surface was more variable. Initially, it was less stable than the weak layer closer to the ground, yet it stabilized faster as the winter progressed. In spring, the simulated snowpack on south-facing slopes stabilized faster than on north-facing slopes, in line with the regional avalanche forecast. In the winter months January to March 2017, simulated instability metrics did not suggest that the snowpack on south-facing slopes was more stable, as reported in the regional avalanche forecast. Although a validation with field data is lacking, these model results still show the potential and challenges of distributed modeling for supporting operational avalanche forecasting.
The correct simulation of the atmospheric boundary layer (ABL) is crucial for reliable weather forecasts in truly complex terrain. However, common assumptions for model parametrizations are only ...valid for horizontally homogeneous and flat terrain. Here, we evaluate the turbulence parametrization of the numerical weather prediction model COSMO with a horizontal grid spacing of
Δ
x
=
1.1
km
for the Inn Valley, Austria. The long-term, high-resolution turbulence measurements of the i-Box measurement sites provide a useful data pool of the ABL structure in the valley and on slopes. We focus on days and nights when ABL processes dominate and a thermally-driven circulation is present. Simulations are performed for case studies with both a one-dimensional turbulence parametrization, which only considers the vertical turbulent exchange, and a hybrid turbulence parametrization, also including horizontal shear production and advection in the budget of turbulence kinetic energy (TKE). We find a general underestimation of TKE by the model with the one-dimensional turbulence parametrization. In the simulations with the hybrid turbulence parametrization, the modelled TKE has a more realistic structure, especially in situations when the TKE production is dominated by shear related to the afternoon up-valley flow, and during nights, when a stable ABL is present. The model performance also improves for stations on the slopes. An estimation of the horizontal shear production from the observation network suggests that three-dimensional effects are a relevant part of TKE production in the valley.
A scheme to represent the impact of urban buildings on airflow in mesoscale atmospheric models is presented. In the scheme, the buildings are not explicitly resolved, but their effects on the ...grid-averaged variables are parameterised. An urban quarter is characterised by a horizontal building size, a street canyon width and a building density as a function of height. The module computes the impact of the horizontal (roof and canyon floor) and vertical (walls) surfaces on the wind speed, temperature and turbulent kinetic energy. The computation of the shortwave and longwave radiation, needed to compute the temperature of the urban surfaces, takes into account the shadowing and radiation trapping effects induced by the urban canyons. The computation of the turbulent length scales in the TKE equation is also modified to take into account the presence of the buildings. The parameterisation is introduced into a mesoscale model and tested in a bidimensional case of a city over flat terrain. The new parameterisation is shown to be able to reproduce the most important features observed in urban areas better than the traditional approach which is based only on the modification of the roughness length, thereby retaining the Monin-Obukhov similarity theory. The new surface exchange parameterisation is furthermore shown to have a strong impact on the dispersion characteristics of air pollutants in urban areas.PUBLICATION ABSTRACT
•Local scaling approach is valid at a forested site with highly inhomogeneous fetch.•Flux-variance and flux-gradient relationships respond differently to inhomogeneity.•Roughness sublayer influence ...is observed only for the flux-variance relationships.•Above the RSL scaled TKE depends linearly on stability.•Surface inhomogeneity has insignificant effect on dimensionless wind shear.
The local scaling approach was examined based on the multi-level measurements of atmospheric turbulence in the wintertime (December 2008–February 2009) stable atmospheric boundary layer (SBL) established over a heterogeneous surface influenced by mixed agricultural, industrial and forest surfaces. The heterogeneity of the surface was characterized by spatial variability of both roughness and topography. Nieuwstadt’s local scaling approach was found to be suitable for the representation of all three wind velocity components. For neutral conditions, values of all three non-dimensional velocity variances were found to be smaller at the lowest measurement level and larger at higher levels in comparison to classical values found over flat terrain. Influence of surface heterogeneity was reflected in the ratio of observed dimensionless standard deviation of the vertical wind component and corresponding values of commonly used similarity formulas for flat and homogeneous terrain showing considerable variation with wind direction. The roughness sublayer influenced wind variances, and consequently the turbulent kinetic energy and correlation coefficients at the lowest measurement level, but not the wind shear profile. The observations support the classical linear expressions for the dimensionless wind shear (ϕm) even over inhomogeneous terrain after removing data points associated with the flux Richardson number (Rf) greater than 0.25. Leveling-off of ϕm at higher stabilities was found to be a result of the large number of data characterized by small-scale turbulence (Rf>0.25). Deviations from linear expressions were shown to be mainly due to this small-scale turbulence rather than due to the surface heterogeneities, supporting the universality of this relationship. Additionally, the flux-gradient dependence on stability did not show different behavior for different wind regimes, indicating that the stability parameter is sufficient predictor for flux-gradient relationship. Data followed local z-less scaling for ϕm when the prerequisite Rf≤0.25 was imposed.