Newly improved formulations of the bare soil evaporation and the surface temperature are presented, using the multilayer land surface scheme TERRA of the Consortium for Small-scale Modeling (COSMO) atmospheric model. The simulations were carried out in offline mode with atmospheric forcing data from the Meteorological Observatory Lindenberg–Richard-Aßmann-Observatory of the German Meteorological Service. The results show that the bare soil evaporation simulated by the reference version of TERRA is substantially overestimated under wet conditions, and underestimated under dry conditions. Furthermore, the amplitude of the diurnal cycle of the surface temperature is systematically underestimated. In contrast, the diurnal cycles of the temperatures in the soil are overestimated instead. The new description of the bare soil evaporation in TERRA is based on a resistance formulation analogue to Ohm’s law, while the surface temperature is now based on the skin temperature formulation by Viterbo and Beljaars. The new formulation improves the simulated bare soil evaporation substantially. In particular, the overestimation under wet conditions is reduced, also acting against an extensive drying of the soil during the annual cycle. Additionally, the underestimation under dry conditions is reduced as well. Furthermore, the simulated amplitude of the diurnal cycle of the surface temperature is substantially increased. In particular, a nocturnal warm bias is systematically reduced. In addition to this, the new formulations were also applied in coupled mode in the COSMO model, resulting in improved diurnal cycles of near-surface temperature and dew point.

A realistic representation of processes that are not resolved by the model grid is one of the key challenges in Earth-system modelling. In particular, the non-linear nature of processes involved makes a representation of the link between the atmosphere and the land surface difficult. This is especially so when the land surface is horizontally strongly heterogeneous. In the majority of present day Earth system models two strategies are pursued to couple the land surface and the atmosphere. In the first approach, surface heterogeneity is not explicitly accounted for, instead effective parameters are used to represent the entirety of the land surface in a model’s grid box (parameter-aggregation). In the second approach, subgrid-scale variability at the surface is explicitly represented, but it is assumed that the blending height is located below the lowest atmospheric model level (simple flux-aggregation). Thus, in both approaches the state of the atmosphere is treated as being horizontally homogeneous within a given grid box. Based upon the blending height concept, an approach is proposed that allows for a land-surface–atmosphere coupling in which horizontal heterogeneity is considered not only at the surface, but also within the lowest layers of the atmosphere (the VERTEX scheme). Below the blending height, the scheme refines the turbulent mixing process with respect to atmospheric subgrid fractions, which correspond to different surface features. These subgrid fractions are not treated independently of each other, since an explicit horizontal component is integrated into the turbulent mixing process. The scheme was implemented into the JSBACH model, the land component of the Max Planck Institute for Meteorology’s Earth-system model, when coupled to the atmospheric general circulation model ECHAM. The single-column version of the Earth system model is used in two example cases in order to demonstrate how the effects of surface heterogeneity are transferred into the atmosphere, influencing local stability and the turbulent mixing process. Furthermore, a simple flux-aggregation scheme was implemented into the JSBACH model. By comparing single-column simulations utilizing the VERTEX scheme and the simple flux-aggregation scheme, it can be shown that the horizontal disaggregation of the turbulent mixing process has a substantial impact on the simulated mean state of a grid box. Here, the differences between simulations with the two schemes may, in certain cases, be even larger than the differences between simulations with the simple flux-aggregation scheme and simulations in which surface heterogeneity is not explicitly accounted for (i.e., a parameter-aggregation scheme).

The increase in built surfaces constitutes the main reason for the formation of the Urban Heat Island (UHI), that is a metropolitan area significantly warmer than its surrounding rural areas. The urban heat islands and other urban-induced climate feedbacks may amplify heat stress and urban flooding under climate change and therefore to predict them correctly has become essential. Currently in the COSMO model, cities are represented by natural land surfaces with an increased surface roughness length and a reduced vegetation cover, but this approach is unable to correctly reproduce the UHI effect. By increasing the model resolution, a representation of the main physical processes that characterize the urban local meteorology should be addressed, in order to better forecast temperature, moisture and precipitation in urban environments. Within the COSMO Consortium a bulk parameterization scheme (TERRA_URB or TU) has been developed. It parametrizes the effects of buildings, streets and other man-made impervious surfaces on energy, moist and momentum exchanges between the surface and atmosphere, and additionally accounts for the anthropogenic heat flux as a heat source from the surface to the atmosphere. TU implements an impervious water-storage parameterization, and the Semi-empirical Urban canopy parametrization (SURY) that translates 3D urban canopy into bulk parameters. This paper presents evaluation results of the TU scheme in high-resolution simulations with a recent COSMO model version for selected European cities, namely Turin, Naples and Moscow. The key conclusion of the work is that the TU scheme in the COSMO model reasonably reproduces UHI effect and improves air temperature forecasts for all the investigated urban areas, despite each city has very different morphological characteristics. Our results highlight potential benefits of a new turbulence scheme and the representation of skin-layer temperature (for vegetation) in the model performance. Our model framework provides perspectives for enhancing urban climate modelling, although further investigations in improving model parametrizations, calibration and the use of more realistic urban canopy parameters are needed.

Feedbacks of plant phenology to the regional climate system affect fluxes of energy, water, CO2, biogenic volatile organic compounds as well as canopy conductance, surface roughness length, and are influencing the seasonality of albedo. We performed simulations with the regional climate model COSMO-CLM (CCLM) at three locations in Germany covering the period 1999 to 2015 in order to study the sensitivity of grass phenology to different environmental conditions by implementing a new phenology module. We provide new evidence that the annually-recurring standard phenology of CCLM is improved by the new calculation of leaf area index (LAI) dependent upon surface temperature, day length, and water availability. Results with the new phenology implemented in the model show a significantly higher correlation with observations than simulations with the standard phenology. The interannual variability of LAI improves the representation of vegetation in years with extremely warm winter/spring (e.g., 2007) or extremely dry summer (e.g., 2003) and shows a more realistic growth period. The effect of the newly implemented phenology on atmospheric variables is small but tends to be positive. It should be used in future applications with an extension on more plant functional types.

A new urban canopy scheme for the ICON atmospheric model is presented. Increasing the resolution of atmospheric models for numerical weather prediction (NWP) or climate simulations allows, among others, for a more realistic description of the processes at the land surface. Here, one field of growing interest are the processes in urban areas. Beside their relevance for the meteorological modelling, there is a general trend in most countries that the number of people living in towns is significantly increasing. During the recent years, an urban canopy parameterization was developed for the multi-layer land surface scheme TERRA of the Consortium for Small-scale Modeling (COSMO) mesoscale atmospheric model. This parameterisation, TERRA_URB, originally developed for the climate version of COSMO and then ported to the NWP version, was shown to be able to reproduce the key urban meteorological features for different European cities. In the framework of the transition of the COSMO Consortium to the ICON model, TERRA_URB needs to be implemented in ICON. Furthermore, an updated set of urban canopy parameters needs to be provided, for describing the urban characteristics down to a mesh size of 1 km, and below. For these purposes, the COSMO Consortium organises the dedicated Priority Project CITTA’. First results are presented for TERRA_URB in the ICON limited-area model ICON-LAM for different cities of interest of the CITTA’ partners. The preliminary results indicate already that urban features like the urban heat island effect are well represented. This is in agreement with the experiences with TERRA_URB in the COSMO model, both the climate as well as the NWP version.

Land surface processes have a significant impact on near-surface atmospheric phenomena. They determine, among others, near-surface sensible and latent heat fluxes and the radiation budget, and thus influence atmosphere and land characteristics, such as temperature and humidity, the structure of the planetary boundary layer, and even cloud formation processes. It is therefore important to simulate the land surface processes in atmospheric models as realistically as possible.Verifications have shown that the amplitude of the diurnal cycle of the surface temperature simulated by the land surface scheme TERRA of the COSMO atmospheric model is systematically underestimated. In contrast, the diurnal cycles of the temperatures in the soil are overestimated, instead. This means that the other components of the surface energy balance are biased as well, for instance, the surface turbulent heat fluxes or the ground heat flux.Data from the Meteorological Observatory Lindenberg of the German Meteorological Service (DWD) were used to analyse this model behaviour. In the standard model configuration of TERRA, there is no representation of the vegetation in the surface energy balance. This means, there is no energy budget including a temperature for the vegetation layer. Furthermore, the insulating effects by the vegetation at the sub-canopy level are missing as well. In this work, a scheme providing both of these missing model characteristics was implemented in TERRA. As a result, the simulated diurnal amplitude of the surface temperature is increased and the one of the soil temperature is reduced, both leading to better agreements with the measurements. These improvements are found in TERRA in offline mode, using Lindenberg observations, as well as in coupled mode in the atmospheric models of DWD, i.e. the limited-area COSMO model and the global ICON model.