The most frequently used boundary-layer turbulence parameterization in numerical weather predicti... more The most frequently used boundary-layer turbulence parameterization in numerical weather prediction (NWP) models are turbulence kinetic energy (TKE) based-based schemes. However, these parameterizations suffer from a potential weakness, namely the strong dependence on an ad-hoc quantity, the so-called turbulence length scale. The physical interpretation of the turbulence length scale is difficult and hence it cannot be directly related to measurements or large eddy simulation (LES) data. Consequently, formulations for the turbulence length scale in basically all TKE schemes are based on simplified assumptions and are model-dependent. A good reference for the independent evaluation of the turbulence length scale expression for NWP modeling is missing. Here we propose a new turbulence length scale diagnostic which can be used in the gray zone of turbulence without modifying the underlying TKE turbulence scheme. The new diagnostic is based on the TKE budget: The core idea is to encapsulate the sum of the molecular dissipation and the cross-scale TKE transfer into an effective dissipation, and associate it with the new turbulence length scale. This effective dissipation can then be calculated as a residuum in the TKE budget equation (for horizontal sub-domains of different sizes) using LES data. Estimation of the scale dependence of the diagnosed turbulence length scale using this novel method is presented for several idealized cases.
Measurements of vertical velocity from vertically pointing Doppler lidars are used to derive the ... more Measurements of vertical velocity from vertically pointing Doppler lidars are used to derive the profiles of normalized vertical velocity variance. Observations were taken during the FESSTVaL (Field Experiment on Submesoscale Spatio-Temporal Variability in Lindenberg) campaign during the warm seasons of 2020 and 2021. Normalized by the square of the convective velocity scale, the average vertical velocity variance profile follows the universal profile of Lenschow et al. (1980). However, daily profiles still show a high day-today variability. We found that moisture transport and the content of moisture in the boundary layer could explain the remaining variability of the normalized vertical velocity variance. The magnitude of the normalized vertical velocity variance is highest on clear-sky days and decreases as the absolute humidity increases and surface latent heat flux decreases on cloud-topped days. This suggests that moisture content and moisture transport are limiting factors for the intensity of turbulence in the convective boundary layer. We also found that the intensity of turbulence decreases with an increase in the boundary layer cloud fraction during FESSTVaL, while the latent heating in the cloud layer was not a relevant source of turbulence in this case. We conclude that a new vertical velocity scale has to be defined that would take into account the moist processes in the convective boundary layer.
We evaluate the influence of a forest parametrization on the simulation of the boundary layer flo... more We evaluate the influence of a forest parametrization on the simulation of the boundary layer flow over moderate complex terrain in the context of the Perdigão 2017 field campaign. The numerical simulations are performed using the Weather Research and Forecasting model in large eddy simulation mode (WRF-LES). The short-term, high-resolution (40 m horizontal grid spacing) and long-term (200 m horizontal grid spacing) WRF-LES are evaluated for an integration time of 12 h and 1.5 months, respectively, with and without forest parameterization. The short-term simulations focus on low-level jet events over the valley, while the long-term simulations cover the whole intensive observation period (IOP) of the field campaign. The results are validated using lidar and meteorological tower observations. The mean diurnal cycle during the IOP shows a significant improvement of the along-valley wind speed and the wind direction when using the forest parametrization. However, the drag imposed by the parametrization results in an underestimation of the cross-valley wind speed, which can be attributed to a poor representation of the land surface characteristics. The evaluation of the high-resolution WRF-LES shows a positive influence of the forest parametrization on the simulated winds in the first 500 m above the surface.
Thermally driven local winds are ubiquitous in deep Alpine valleys during fair weather conditions... more Thermally driven local winds are ubiquitous in deep Alpine valleys during fair weather conditions resulting in a unique wind climatology for any given valley. The accurate forecasting of these local wind systems is challenging, as they are the result of complex and multi-scale interactions. Even more so, if the aim is an accurate forecast of the winds from the near-surface to the free atmosphere, which can be considered a prerequisite for the accurate prediction of mountain weather. This study combines the evaluation of the simulated surface winds in several Alpine valleys with a more detailed evaluation of the wind evolution for a particular location in the Swiss Rhone valley, at the town of Sion during the month of September 2016. Four numerical simulations using the COSMO model are evaluated, two using a grid spacing of 1.1 km and two with a grid spacing of 550 m. For each resolution, one simulation is initialised with the soil moisture from the COSMO analysis and one with an inc...
Turbulence kinetic energy (TKE) schemes are routinely used for turbulence parameterization in num... more Turbulence kinetic energy (TKE) schemes are routinely used for turbulence parameterization in numerical weather prediction models. A key component of these schemes is the so-called turbulence length scale. Novel scale-aware, budget-based diagnostics that account for the cross-scale transfer of variances are used to evaluate the performance of selected turbulence length scale formulations in the gray zone of turbulence. The diagnostics are computed using the coarse-graining method on high resolution large eddy simulation data for selected idealized cases. The vertical profiles and the temporal evolution of the turbulence length scales are analyzed. Additionally, the local normalized root mean square error and a non-local three-component technique tailored specifically to the turbulence length scale profiles are used for the evaluation. Based on our analyses, we recommend using turbulence length-scale formulations that depend not only on the boundary layer height, but also on the TKE ...
The single-column mode (SCM) of the ICON (ICOsahedral Nonhydrostatic) modeling framework is prese... more The single-column mode (SCM) of the ICON (ICOsahedral Nonhydrostatic) modeling framework is presented. The primary purpose of the ICON SCM is to use it as a tool for research, model evaluation and development. Thanks to the simplified geometry of the ICON SCM, various aspects of the ICON model, in particular the model physics, can be studied in a well-controlled environment. Additionally, the ICON SCM has a reduced computational cost and a low data storage demand. The ICON SCM can be utilized for idealized cases—several well-established cases are already included—or for semi-realistic cases based on analyses or model forecasts. As the case setup is defined by a single NetCDF file, new cases can be prepared easily by the modification of this file. We demonstrate the usage of the ICON SCM for different idealized cases such as shallow convection, stratocumulus clouds, and radiative transfer. Additionally, the ICON SCM is tested for a semi-realistic case together with an equivalent thre...
We thank the reviewers for their careful reading of our manuscript and constructive comments. The... more We thank the reviewers for their careful reading of our manuscript and constructive comments. Their comments helped us to reline our manuscript by clarifying several sections and by improving the presentation of the results. Point-by-point responses to individual questions and remarks for each reviewer are listed below. Note: The referee's comments are reproduced in extenso in bold characters. The line numbers mentioned in our answers are those of the revised version of the manuscript.
An intercomparison between 10 single-column (SCM) and 5 large-eddy simulation (LES) models is pre... more An intercomparison between 10 single-column (SCM) and 5 large-eddy simulation (LES) models is presented for a radiation fog case study inspired by the LANFEX field campaign. 7 of the SCMs represent single-column equivalents of operational numerical weather prediction (NWP) models, whilst 3 are research-grade SCMs designed for fog simulation, and the LES are designed to reproduce in the best manner currently possible the underlying physical processes governing fog formation. The LES model results are of variable quality, and do not provide a consistent baseline against which to compare the NWP models, particularly under high aerosol or cloud droplet number (CDNC) conditions. The main SCM bias appears to be toward over-development of fog, i.e. fog which is too thick, although the inter-model variability is large. In reality there is a subtle balance between water lost to the surface and water condensed into fog, and the ability of a model to accurately simulate this process strongly determines the quality of its forecast. Some NWP-SCMs do not represent fundamental components of this process (e.g. cloud droplet sedimentation) and therefore are naturally hampered in their ability to deliver accurate simulations. Finally, we show that modelled fog development is as sensitive to the shape of the cloud droplet size distribution, a rarely studied or modified part of the microphysical parametrization, as it is to the underlying aerosol or CDNC. 1 Introduction Most operational numerical weather prediction (NWP) centres will list errors in fog forecasting amongst their top model problems, with the requirement for improvement considered high-priority (Hewson, 2019). The key customer driving this is 1
The most frequently used boundary-layer turbulence parameterization in numerical weather predicti... more The most frequently used boundary-layer turbulence parameterization in numerical weather prediction (NWP) models are turbulence kinetic energy (TKE) based-based schemes. However, these parameterizations suffer from a potential weakness, namely the strong dependence on an ad-hoc quantity, the so-called turbulence length scale. The physical interpretation of the turbulence length scale is difficult and hence it cannot be directly related to measurements or large eddy simulation (LES) data. Consequently, formulations for the turbulence length scale in basically all TKE schemes are based on simplified assumptions and are model-dependent. A good reference for the independent evaluation of the turbulence length scale expression for NWP modeling is missing. Here we propose a new turbulence length scale diagnostic which can be used in the gray zone of turbulence without modifying the underlying TKE turbulence scheme. The new diagnostic is based on the TKE budget: The core idea is to encapsu...
Quarterly Journal of the Royal Meteorological Society, 2020
A realistic simulation of the atmospheric boundary layer (ABL) depends on an accurate representat... more A realistic simulation of the atmospheric boundary layer (ABL) depends on an accurate representation of the land–atmosphere coupling. Land surface temperature (LST) plays an important role in this context and the assimilation of LST can lead to improved estimates of the boundary layer and its processes. We assimilated synthetic satellite LST retrievals derived from a nature run as truth into a fully coupled, state‐of‐the‐art land–atmosphere numeric weather prediction model. As assimilation system a local ensemble transform Kalman filter was used and the control vector was augmented by the soil temperature and humidity. To evaluate the concept of the augmented control vector, two‐day case‐studies with different control vector settings were conducted for clear‐sky periods in March and August 2017. These experiments with hourly LST assimilation were validated against the nature run and overall, the RMSE of atmospheric and soil temperature of the first‐guess (and analysis) were reduced....
Diurnal valley winds frequently form over complex topography, particularly under fair weather con... more Diurnal valley winds frequently form over complex topography, particularly under fair weather conditions, and have a significant impact on the local weather and climate. Since diurnal valley winds result from complex and multi-scale interactions, their representation in numerical weather prediction models is challenging. Better understanding of these local winds based on observations is crucial to improve the accuracy of the forecasts. This study investigates the diurnal evolution of the three-dimensional mean wind structure in a deep Alpine valley, the Rhone valley at Sion, using data from a radar wind profiler and a surface weather station operated continuously from 1 September 2016 to 17 July 2017. In particular, the wind profiler data was analyzed for a subset of days on which fair weather conditions allowed for the full development of thermally driven winds. A pronounced diurnal cycle of the wind speed, as well as a reversal of the wind direction twice per day is documented for...
Journal of Geophysical Research: Atmospheres, 2019
The impact of precipitation in shallow cumulus convection on the moisture variance and third‐orde... more The impact of precipitation in shallow cumulus convection on the moisture variance and third‐order moments of moisture is investigated with the help of large‐eddy simulations. Three idealized simulations based on the Rain in Cumulus over the Ocean field experiment are analyzed: one nonprecipitating, on a smaller domain, and two precipitating cases, on a larger domain with different initial profiles of moisture. Results show that precipitation and the associated cloud organization lead to increased generation of higher‐order moments (HOM) of moisture compared to the nonprecipitating case. To understand the physical mechanism and the role of individual processes in this increase, budgets of HOM of moisture are studied. Microphysics directly decreases the generation of HOM of moisture, but this effect is not dominant. The gradient production term is identified as the main source term in the HOM budgets. The influence of the gradient production term on moisture variance is further exami...
We evaluate the near-surface representation of thermally driven winds in the Swiss Alps in a nume... more We evaluate the near-surface representation of thermally driven winds in the Swiss Alps in a numerical weather prediction model at km-scale resolution. In addition, the influence of grid resolution (2.2 km and 1.1 km), topography filtering, and land surface datasets on the accuracy of the simulated valley winds is investigated. The simulations are evaluated against a comprehensive set of surface observations for an 18-day fair-weather summer period in July 2006. The episode is characterized by strong diurnal wind systems and the formation of shallow convection over the mountains, which transitions to precipitating convection in some areas. The near-surface winds (10 m above ground level) follow a typical diurnal pattern with strong daytime up-valley flow and weaker nighttime down-valley flow. At a 2.2 km resolution the valley winds are poorly simulated for most stations, while at a 1.1 km resolution the diurnal cycle of the valley winds is well represented in most large (e.g., Rhein valley at Chur and Rhone valley at Visp) and medium-sized valleys (e.g., Linth valley at Glarus). In the smaller valleys (e.g., Maggia valley at Cevio), the amplitude of the valley wind is still significantly underestimated, even at a 1.1 km resolution. Detailed sensitivity experiments show that the use of high-resolution land surface datasets, for both the soil characteristics as well as for the land cover, and reduced filtering of the topography are essential to achieve good performance at a 1.1 km resolution.
The exchange of heat, momentum, and mass in the atmosphere over mountainous terrain is controlled... more 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 improvem...
Diurnal moist convection is an important element of summer precipitation over Central Europe and ... more Diurnal moist convection is an important element of summer precipitation over Central Europe and the Alps. It is poorly represented in models using parameterized convection. In this study, we investigate the diurnal cycle of convection during 11 days in June 2007 using the COSMO model. The numerical simulations are compared with satellite measurements of GERB and SEVIRI, AVHRR satellite-based cloud properties and ground-based precipitation and temperature measurements. The simulations use horizontal resolutions of 12 km (convection-parameterizing model, CPM) and 2 km (convection-resolving model, CRM) and either a one-moment microphysics scheme (1M) or a two-moment microphysics scheme (2M). They are conducted for a computational domain that covers an extended Alpine area from Northern Italy to Northern Germany. The CPM with 1M exhibits a significant overestimation of high cloud cover. This results in a compensation effect in the top of the atmosphere energy budget due to an underestimation of outgoing longwave radiation (OLR) and an overestimation of reflected solar radiation (RSR). The CRM reduces high cloud cover and improves the OLR bias from a domain mean of −20.1 to −2.6 W/m 2. When using 2M with ice sedimentation in the CRM, high cloud cover is further reduced. The stronger diurnal cycle of high cloud cover and associated convection over the Alps, compared to less mountainous regions, is well represented by the CRM but underestimated by the CPM. Despite substantial differences in high cloud cover, the use of a 2M has no significant impact on the diurnal cycle of precipitation. Furthermore, a negative mid-level cloud bias is found for all simulations.
Uploads
Papers by Juerg Schmidli