Sebastian Schubert

The Influence of green areas and roof albedos on air temperatures during Extreme Heat Events in Berlin, Germany

The mesoscale atmospheric model COSMO-CLM (CCLM) with the Double Canyon Effect Parametrization Scheme (DCEP) is applied to investigate possible adaption measures to extreme heat events (EHEs) for the city of Berlin, Germany. The emphasis is on the effects of a modified urban vegetation cover and roof albedo on near-surface air temperatures. Five EHEs with a duration of 5 days or more are identified for the period 2000 to 2009. A reference simulation is carried out for each EHE with current vegetation cover, roof albedo and urban canopy parameters (UCPs), and is evaluated with temperature observations from weather stations in Berlin and its surroundings. The derivation of the UCPs from an impervious surface map and a 3-D building data set is detailed. Characteristics of the simulated urban heat island for each EHE are analysed in terms of these UCPs. In addition, six sensitivity runs are examined with a modified vegetation cover of each urban grid cell by 25%, 5% and 15%, with a roof albedo increased to 0.40 and 0.65, and with a combination of the largest vegetation cover and roof albedo, respectively. At the weather stations’ grid cells, the results show a maximum of the average diurnal change in air temperature during each EHE of 0.82 K and 0.48 K for the 25% and 15% vegetation covers, 0.50 K for the roof albedos of 0.65, and 0.63 K for the combined vegetation and albedo case. The largest effects on the air temperature are detected during midday.

A Double-Canyon Radiation Scheme for Multi-Layer Urban Canopy Models

We develop a double-canyon radiation scheme (DCEP) for urban canopy models embedded in mesoscale numerical models based on the Building Effect Parametrization (BEP). The new scheme calculates the incoming and outgoing longwave and shortwave radiation for roof, wall and ground surfaces for an urban street canyon characterized by its street and building width, canyon length, and the building height distribution. The scheme introduces the radiative interaction of two neighbouring urban canyons allowing the full inclusion of roofs into the radiation exchange both inside the canyon and with the sky. In contrast to BEP, we also treat direct and diffuse shortwave radiation from the sky independently, thus allowing calculation of the effective parameters representing the urban diffuse and direct shortwave radiation budget inside the mesoscale model. Furthermore, we close the energy balance of incoming longwave and diffuse shortwave radiation from the sky, so that the new scheme is physically more consistent than the BEP scheme. Sensitivity tests show that these modifications are important for urban regions with a large variety of building heights. The evaluation against data from the Basel Urban Boundary Layer Experiment indicates a good performance of the DCEP when coupled with the regional weather and climate model COSMO-CLM.

Dynamical analysis of blocking events: spatial and temporal fluctuations of covariant Lyapunov vectors

One of the most relevant weather regimes in the midlatitude atmosphere is the persistent deviation from the approximately zonally symmetric jet stream leading to the emergence of so-called blocking patterns. Such configurations are usually connected to exceptional local stability properties of the flow which come along with an improved local forecast skills during the phenomenon. It is instead extremely hard to predict onset and decay of blockings. Covariant Lyapunov Vectors (CLVs) offer a suitable characterization of the linear stability of a chaotic flow, since they represent the full tangent linear dynamics by a covariant basis which explores linear perturbations at all time scales. Therefore, we assess whether CLVs feature a signature of the blockings. As a first step, we examine the CLVs for a quasi-geostrophic beta-plane two-layer model in a periodic channel baroclinically driven by a meridional temperature gradient ΔT. An orographic forcing enhances the emergence of localized blocked regimes. We detect the blocking events of the channel flow with a Tibaldi-Molteni scheme adapted to the periodic channel. When blocking occurs, the global growth rates of the fastest growing CLVs are significantly higher. Hence, against intuition, the circulation is globally more unstable in blocked phases. Such an increase in the finite time Lyapunov exponents with respect to the long term average is attributed to stronger barotropic and baroclinic conversion in the case of high temperature gradients, while for low values of ΔT, the effect is only due to stronger barotropic instability. In order to determine the localization of the CLVs we compare the meridionally averaged variance of the CLVs during blocked and unblocked phases. We find that on average the variance of the CLVs is clustered around the center of blocking. These results show that the blocked flow affects all time scales and processes described by the CLVs.

Covariant Lyapunov vectors of a quasi‐geostrophic baroclinic model: analysis of instabilities and feedbacks

The classical approach for studying atmospheric variability is based on defining a background state and studying the linear stability of the small fluctuations around such a state. Weakly nonlinear theories can be constructed using higher order expansion terms. While these methods undoubtedly have great value for elucidating the relevant physical processes, they are unable to follow the dynamics of a turbulent atmosphere. We provide a first example of the extension of classical stability analysis to a nonlinearly evolving quasi-geostrophic flow. The so-called covariant Lyapunov vectors (CLVs) provide a covariant basis describing the directions of exponential expansion and decay of perturbations to the nonlinear trajectory of the flow. We use such a formalism to re-examine the basic barotropic and baroclinic processes of the atmosphere with a quasi-geostrophic beta-plane two-layer model in a periodic channel driven by a forced meridional temperature gradient ΔT. We explore three settings of ΔT, representative of relatively weak turbulence, well-developed turbulence and intermediate conditions. We construct the Lorenz energy cycle for each CLV describing the energy exchanges with the background state. A positive baroclinic conversion rate is a necessary but not sufficient condition for instability. Barotropic instability is present only for a few very unstable CLVs for large values of ΔT. Slowly growing and decaying hydrodynamic Lyapunov modes closely mirror the properties of the background flow. Following the classical necessary conditions for barotropic/baroclinic instability, we find a clear relationship between the properties of the eddy fluxes of a CLV and its instability. CLVs with positive baroclinic conversion seem to form a set of modes for constructing a reduced model of the atmospheric dynamics.

Exploring the Lyapunov instability properties of high-dimensional atmospheric and climate models

The stability properties of intermediate-order climate models are investigated by computing their Lyapunov exponents (LEs). The two models considered are PUMA (Portable University Model of the Atmosphere), a primitive-equation simple general circulation model, and MAOOAM (Modular Arbitrary-Order Ocean-Atmosphere Model), a quasi-geostrophic coupled ocean-atmosphere model on a beta-plane. We wish to to investigate the effect of the different levels of filtering on the instabilities and dynamics of the atmospheric flows. Moreover, we assess the impact of the oceanic coupling, the dissipation scheme and the resolution on the spectra of LEs. The PUMA Lyapunov spectrum is computed for two different values of the meridional temperature gradient defining the Newtonian forcing. The increase of the gradient gives rise to a higher baroclinicity and stronger instabilities, corresponding to a larger dimension of the unstable manifold and a larger first LE. The convergence rate of the rate functional for the large deviation law of the finite-time Lyapunov exponents (FTLEs) is fast for all exponents, which can be interpreted as resulting from the absence of a clear-cut atmospheric time-scale separation in such a model. The MAOOAM spectra show that the dominant atmospheric instability is correctly represented even at low resolutions. However, the dynamics of the central manifold, which is mostly associated to the ocean dynamics, is not fully resolved because of its associated long time scales, even at intermediate orders. This paper highlights the need to investigate the natural variability of the atmosphere-ocean coupled dynamics by associating rate of growth and decay of perturbations to the physical modes described using the formalism of the covariant Lyapunov vectors and to consider long integrations in order to disentangle the dynamical processes occurring at all time scales.