Atmospheric transport and chemistry of trace gases in LMDz5B: evaluation and implications for inverse modelling

Representation of atmospheric transport is a major source of error in the estimation of greenhouse gas sources and sinks by inverse modelling. Here we assess the impact on trace gas mole fractions of the new physical parameterizations recently implemented in the atmospheric global climate model LMDz to improve vertical diffusion, mesoscale mixing by thermal plumes in the planetary boundary layer (PBL), and deep convection in the troposphere. At the same time, the horizontal and vertical resolution of the model used in the inverse system has been increased. The aim of this paper is to evaluate the impact of these developments on the representation of trace gas transport and chemistry, and to anticipate the implications for inversions of greenhouse gas emissions using such an updated model. Comparison of a one-dimensional version of LMDz with large eddy simulations shows that the thermal scheme simulates shallow convective tracer transport in the PBL over land very efficiently, and much better than previous versions of the model. This result is confirmed in three-dimensional simulations, by a much improved reproduction of the radon-222 diurnal cycle. However, the enhanced dynamics of tracer concentrations induces a stronger sensitivity of the new LMDz configuration to external meteorological forcings. At larger scales, the inter-hemispheric exchange is slightly slower when using the new version of the model, bringing them closer to observations. The increase in the vertical resolution (from 19 to 39 layers) significantly improves the representation of stratosphere/troposphere exchange. Furthermore, changes in atmospheric thermodynamic variables, such as temperature, due to changes in the PBL mixing modify chemical reaction rates, which perturb chemical equilibriums of reactive trace gases. One implication of LMDz model developments for future inversions of greenhouse gas emissions is the ability of the updated system to assimilate a larger amount of high-frequency data sampled at high-variability stations. Others implications are discussed at the end of the paper

A GCSS model intercomparison for a tropical squall line observed during toga-coare. II: Intercomparison of single-column models and a cloud-resolving model

This paper presents single-column model (SCM) simulations of a tropical squall-line case observed during the Coupled Ocean-Atmosphere Response Experiment of the Tropical Ocean/Global Atmosphere Programme. This case-study was part of an international model intercomparison project organized by Working Group 4 ‘Precipitating Convective Cloud Systems’ of the GEWEX (Global Energy and Water-cycle Experiment) Cloud System Study. Eight SCM groups using different deep-convection parametrizations participated in this project. The SCMs were forced by temperature and moisture tendencies that had been computed from a reference cloud-resolving model (CRM) simulation using open boundary conditions. The comparison of the SCM results with the reference CRM simulation provided insight into the ability of current convection and cloud schemes to represent organized convection. The CRM results enabled a detailed evaluation of the SCMs in terms of the thermodynamic structure and the convective mass flux of the system, the latter being closely related to the surface convective precipitation. It is shown that the SCMs could reproduce reasonably well the time evolution of the surface convective and stratiform precipitation, the convective mass flux, and the thermodynamic structure of the squall-line system. The thermodynamic structure simulated by the SCMs depended on how the models partitioned the precipitation between convective and stratiform. However, structural differences persisted in the thermodynamic profiles simulated by the SCMs and the CRM. These differences could be attributed to the fact that the total mass flux used to compute the SCM forcing differed from the convective mass flux. The SCMs could not adequately represent these organized mesoscale circulations and the microphysicallradiative forcing associated with the stratiform region. This issue is generally known as the ‘scale-interaction’ problem that can only be properly addressed in fully three-dimensional simulations. Sensitivity simulations run by several groups showed that the time evolution of the surface convective precipitation was considerably smoothed when the convective closure was based on convective available potential energy instead of moisture convergence. Finally, additional SCM simulations without using a convection parametrization indicated that the impact of a convection parametrization in forced SCM runs was more visible in the moisture profiles than in the temperature profiles because convective transport was particularly important in the moisture budget

Does convective aggregation need to be represented in cumulus parameterizations?

Tropical deep convection exhibits a variety of levels of aggregation over a wide range of scales. Based on a multisatellite analysis, the present study shows at mesoscale that different levels of aggregation are statistically associated with differing large-scale atmospheric states, despite similar convective intensity and large-scale forcings. The more aggregated the convection, the dryer and less cloudy the atmosphere, the stronger the outgoing longwave radiation, and the lower the planetary albedo. This suggests that mesoscale convective aggregation has the potential to affect couplings between moisture and convection and between convection, radiation, and large-scale ascent. In so doing, aggregation may play a role in phenomena such as “hot spots” or the Madden-Julian Oscillation. These findings support the need for the representation of mesoscale organization in cumulus parameterizations; most parameterizations used in current climate models lack any such representation. The ability of a cloud system-resolving model to reproduce observed relationships suggests that such models may be useful to guide attempts at parameterizations of convective aggregation

Intercomparison of the northern hemisphere winter mid-latitude atmospheric variability of the IPCC models

We compare, for the overlapping time frame 1962-2000, the estimate of the northern hemisphere (NH) mid-latitude winter atmospheric variability within the XX century simulations of 17 global climate models (GCMs) included in the IPCC-4AR with the NCEP and ECMWF reanalyses. We compute the Hayashi spectra of the 500hPa geopotential height fields and introduce an integral measure of the variability observed in the NH on different spectral sub-domains. Only two high-resolution GCMs have a good agreement with reanalyses. Large biases, in most cases larger than 20%, are found between the wave climatologies of most GCMs and the reanalyses, with a relative span of around 50%. The travelling baroclinic waves are usually overestimated, while the planetary waves are usually underestimated, in agreement with previous studies performed on global weather forecasting models. When comparing the results of various versions of similar GCMs, it is clear that in some cases the vertical resolution of the atmosphere and, somewhat unexpectedly, of the adopted ocean model seem to be critical in determining the agreement with the reanalyses. The GCMs ensemble is biased with respect to the reanalyses but is comparable to the best 5 GCMs. This study suggests serious caveats with respect to the ability of most of the presently available GCMs in representing the statistics of the global scale atmospheric dynamics of the present climate and, a fortiori, in the perspective of modelling climate change.Comment: 39 pages, 8 figures, 2 table

An object-based model for convective cold pool dynamics

A simple model of the organization of atmospheric moist convection by cold outflows is presented. The model consists of two layers: a lower layer where instability gradually builds up, and an upper layer where instability is rapidly released. Its formulation is inspired by Abelian sandpile models: instability and convection are both represented in terms of particles that are coupled to a lattice grid. An excess of particles in the lower layer triggers a particle release into the upper (cloud) layer. Particles in the upper layer also induce particle movement in the lower layer: this reverse coupling represents the effect of precipitation and the associated cold outflows. The model shows two behavioral regimes. Activity is scattered when the reverse coupling is weak, but when it is strong, convection forms cellular patterns. Though this model does not contain a detailed representation of physical processes in convection, it captures some key dynamical features of precipitating convection seen in satellite observations and LES studies. These include the formation of open cells, temporal oscillations in convective intensity, hysteresis, and the effect of precipitation on the scale of convection. We argue that an object-based representation of convection may be able to capture properties of convective organization that are missing in traditional parameterizations

Departures from convective equilibrium with a rapidly-varying surface forcing

Convective equilibrium is a long-standing and useful concept for understanding many aspects of the behaviour of deep moist convection. For example, it is often invoked in developing parameterizations for large-scale models. However, the equilibrium assumption may begin to break down as models are increasingly used with shorter timesteps and finer resolutions. Here we perform idealized cloud-system resolving model simulations of deep convection with imposed time variations in the surface forcing. A range of rapid forcing timescales from 1 − 36hr are used, in order to induce systematic departures from equilibrium. For the longer forcing timescales, the equilibrium assumption remains valid, in at least the limited sense that cycle-integrated measures of convective activity are very similar from cycle to cycle. For shorter forcing timescales, cycle-integrated convection becomes more variable, with enhanced activity on one cycle being correlated with reduced activity on the next, suggesting a role for convective memory. Further investigation shows that the memory does not appear to be carried by the domain-mean thermodynamic fields but rather by structures on horizontal scales of 5 − 20km. Such structures are produced by the convective clouds and can persist beyond the lifetime of the cloud, even through to the next forcing cycle

Modeling Haboob Dust Storms in Large-Scale Weather and Climate Models

Recent eld campaigns have shown that haboob dust storms, formed by convective cold pool outflows, contribute a significant fraction of dust uplift over the Sahara and Sahel in summer. However, in-situ observations are sparse and haboobs are frequently concealed by clouds in satellite imagery. Furthermore, most large-scale weather and climate models lack haboobs, because they do not explicitly represent convection. Here a one-year long model run with explicit representation of convection delivers the first full seasonal cycle of haboobs over northern Africa. Using conservative estimates, the model suggests that haboobs contribute one fifth of the annual dust-generating winds over northern Africa, one fourth between May and October, and one third over the western Sahel during this season. A simple parameterization of haboobs has recently been developed for models with parameterized convection, based on the downdraft mass flux of convection schemes. It is applied here to two model runs with different horizontal resolutions, and assessed against the explicit run. The parameterization succeeds in capturing the geographical distribution of haboobs and their seasonal cycle over the Sahara and Sahel. It can be tuned to the different horizontal resolutions, and different formulations are discussed with respect to the frequency of extreme events. The results show that the parameterization is reliable and may solve a major and long-standing issue in simulating dust storms in large-scale weather and climate models

A seamless assessment of the role of convection in the water cycle of the West African Monsoon

A suite of 40 day UK Met Office Unified Model simulations over West Africa during summer 2006 are analyzed to investigate the causes of biases in the position of the rainbelt and to understand the role of convection in the regional water budget. The simulations include climate, global operational, and limited area runs (grid spacings from 1.5 to 40 km), including two 12 km runs, one with parameterized and one with explicit convection. The most significant errors in the water cycle terms occur in the simulations with parameterized convection, associated with the diurnal cycle and the location of the convection. Errors in the diurnal cycle increase the northward advection of moisture out of the Sahel toward the Sahara but decrease the advection of moisture into the Sahel from further south, which limits the availability of moisture for Sahelian rainfall. These biases occur within the first 24 h, showing that they originate from the representation of fast physical processes, specifically, the convection scheme. Once these rainfall regimes have been established, the terms of the water budgets act to reinforce the biases, effectively locking the rainbelt's latitude. One of the simulations with parameterized convection does, however, produce a better latitudinal distribution of rainfall because on the first day it is better able to trigger convection in the Sahel. Accurate representation of the diurnal cycle of convection and the ability to trigger convection in a high convective inhibition environment is key to capturing the water cycle of the region and will improve the representation of the West African Monsoon