ALBERT

Description: A new version of the Global Model of Aerosol Processes (GLOMAP) is described, which uses a two-moment pseudo-modal aerosol dynamics approach rather than the original two-moment bin scheme. GLOMAP-mode simulates the multi-component global aerosol, resolving sulfate, sea-salt, dust, black carbon (BC) and particulate organic matter (POM), the latter including primary and biogenic secondary POM. Aerosol processes are simulated in a size-resolved manner including primary emissions, secondary particle formation by binary homogeneous nucleation of sulfuric acid and water, particle growth by coagulation, condensation and cloud-processing and removal by dry deposition, in-cloud and below-cloud scavenging. A series of benchmark observational datasets are assembled against which the skill of the model is assessed in terms of normalised mean bias (b) and correlation coefficient (R). Overall, the model performs well against the datasets in simulating concentrations of aerosol precursor gases, chemically speciated particle mass, condensation nuclei (CN) and cloud condensation nuclei (CCN). Surface sulfate, sea-salt and dust mass concentrations are all captured well, while BC and POM are biased low (but correlate well). Surface CN concentrations compare reasonably well in free troposphere and marine sites, but are underestimated at continental and coastal sites related to underestimation of either primary particle emissions or new particle formation. The model compares well against a compilation of CCN observations covering a range of environments and against vertical profiles of size-resolved particle concentrations over Europe. The simulated global burden, lifetime and wet removal of each of the simulated aerosol components is also examined and each lies close to multi-model medians from the AEROCOM model intercomparison exercise.

Description: In 2008 the temporal focus of the Palaeoclimate Modelling Intercomparison Project was expanded to include a model intercomparison for the mid-Pliocene warm period (3.29–2.97 million years ago). This project is referred to as PlioMIP (Pliocene Model Intercomparison Project). Two experiments have been agreed upon and comprise phase 1 of PlioMIP. The first (Experiment 1) will be performed with atmosphere-only climate models. The second (Experiment 2) will utilise fully coupled ocean-atmosphere climate models. The aim of this paper is to provide a detailed model intercomparison project description which documents the experimental design in a more detailed way than has previously been done in the literature. Specifically, this paper describes the experimental design and boundary conditions that will be utilised for Experiment 1 of PlioMIP.

Description: The Model for Ozone and Related chemical Tracers, version 4 (MOZART-4) is an offline global chemical transport model particularly suited for studies of the troposphere. The updates of the model from its previous version MOZART-2 are described, including an expansion of the chemical mechanism to include more detailed hydrocarbon chemistry and bulk aerosols. Online calculations of a number of processes, such as dry deposition, emissions of isoprene and monoterpenes and photolysis frequencies, are now included. Results from an eight-year simulation (2000–2007) are presented and evaluated. The MOZART-4 source code and standard input files are available for download from the NCAR Community Data Portal (http://cdp.ucar.edu).

Description: The use of global three-dimensional (3-D) models with satellite observations of CO2 in inverse modeling studies is an area of growing importance for understanding Earth's carbon cycle. Here we use the GEOS-Chem model (version 8-02-01) CO2 mode with multiple modifications in order to assess their impact on CO2 forward simulations. Modifications include CO2 surface emissions from shipping (~0.19 Pg C yr−1), 3-D spatially-distributed emissions from aviation (~0.16 Pg C yr−1), and 3-D chemical production of CO2 (~1.05 Pg C yr−1). Although CO2 chemical production from the oxidation of CO, CH4 and other carbon gases is recognized as an important contribution to global CO2, it is typically accounted for by conversion from its precursors at the surface rather than in the free troposphere. We base our model 3-D spatial distribution of CO2 chemical production on monthly-averaged loss rates of CO (a key precursor and intermediate in the oxidation of organic carbon) and apply an associated surface correction for inventories that have counted emissions of CO2 precursors as CO2. We also explore the benefit of assimilating satellite observations of CO into GEOS-Chem to obtain an observation-based estimate of the CO2 chemical source. The CO assimilation corrects for an underestimate of atmospheric CO abundances in the model, resulting in increases of as much as 24% in the chemical source during May–June 2006, and increasing the global annual estimate of CO2 chemical production from 1.05 to 1.18 Pg C. Comparisons of model CO2 with measurements are carried out in order to investigate the spatial and temporal distributions that result when these new sources are added. Inclusion of CO2 emissions from shipping and aviation are shown to increase the global CO2 latitudinal gradient by just over 0.10 ppm (~3%), while the inclusion of CO2 chemical production (and the surface correction) is shown to decrease the latitudinal gradient by about 0.40 ppm (~10%) with a complex spatial structure generally resulting in decreased CO2 over land and increased CO2 over the oceans. Since these CO2 emissions are omitted or misrepresented in most inverse modeling work to date, their implementation in forward simulations should lead to improved inverse modeling estimates of terrestrial biospheric fluxes.

Description: Successful simulation of cloud-aerosol interactions (indirect aerosol effects) in climate models requires relating grid-scale aerosol, dynamic, and thermodynamic fields to small-scale processes like aerosol activation. A turbulence and cloud parameterization, based on multi-variate probability density functions of sub-grid vertical velocity, temperature, and moisture, has been extended to treat aerosol activation. Multi-variate probability density functions with dynamics (MVD PDFs) offer a solution to the problem of the gap between the resolution of climate models and the scales relevant for aerosol activation and a means to overcome the limitations of diagnostic estimates of cloud droplet number concentration based only on aerosol concentration. Incorporated into the single-column version of GFDL AM3, the MVD PDFs successfully simulate cloud properties including precipitation for cumulus, stratocumulus, and cumulus-under-stratocumulus. The extension to treat aerosol activation predicts droplet number concentrations in good agreement with large eddy simulations (LES). The droplet number concentrations from the MVD PDFs match LES results more closely than diagnostic relationships between aerosol concentration and droplet concentration. In the single-column model simulations, as aerosol concentration increases, droplet concentration increases, precipitation decreases, but liquid water path can increase or decrease.

Description: The current version of the Dutch Atmospheric Large-Eddy Simulation (DALES) is presented. DALES is a large-eddy simulation code designed for studies of the physics of the atmospheric boundary layer, including convective and stable boundary layers as well as cloudy boundary layers. In addition, DALES can be used for studies of more specific cases, such as flow over sloping or heterogeneous terrain, and dispersion of inert and chemically active species. This paper contains an extensive description of the physical and numerical formulation of the code, and gives an overview of its applications and accomplishments in recent years.

Description: We present a new aerosol microphysics and gas aerosol partitioning submodel (Global Modal-aerosol eXtension, GMXe) implemented within the ECHAM/MESSy Atmospheric Chemistry model (EMAC, version 1.8). The submodel is computationally efficient and is suitable for medium to long term simulations with global and regional models. The aerosol size distribution is treated using 7 log-normal modes and has the same microphysical core as the M7 submodel (Vignati et al., 2004). The main developments in this work are: (i) the extension of the aerosol emission routines and the M7 microphysics, so that an increased (and variable) number of aerosol species can be treated (new species include sodium and chloride, and potentially magnesium, calcium, and potassium), (ii) the coupling of the aerosol microphysics to a choice of treatments of gas/aerosol partitioning to allow the treatment of semi-volatile aerosol, and, (iii) the implementation and evaluation of the developed submodel within the EMAC model of atmospheric chemistry. Simulated concentrations of black carbon, particulate organic matter, dust, sea spray, sulfate and ammonium aerosol are shown to be in good agreement with observations (for all species at least 40% of modeled values are within a factor of 2 of the observations). The distribution of nitrate aerosol is compared to observations in both clean and polluted regions. Concentrations in polluted continental regions are simulated quite well, but there is a general tendency to overestimate nitrate, particularly in coastal regions (geometric mean of modelled values/geometric mean of observed data ≈2). In all regions considered more than 40% of nitrate concentrations are within a factor of two of the observations. Marine nitrate concentrations are well captured with 96% of modeled values within a factor of 2 of the observations.

Description: Ideally, a validation and assimilation scheme should maintain the physical principles embodied in the model and be able to evaluate and assimilate lower dimensional features (e.g., discontinuities) contained within a bulk simulation, even when these features are not directly observed or represented by model variables. We present such a scheme and suggest its potential to resolve or alleviate some outstanding problems that stem from making and applying required, yet often non-physical, assumptions and procedures in common operational data assimilation. As proof of concept, we use a sea-ice model with remotely sensed observations of leads in a one-step assimilation cycle. Using the new scheme in a sixteen day simulation experiment introduces model skill (against persistence) several days earlier than in the control run, improves the overall model skill and delays its drop off at later stages of the simulation. The potential and requirements to extend this scheme to different applications, and to both empirical and statistical multivariate and full cycle data assimilation schemes, are discussed.

Description: Calculating the equilibrium composition of atmospheric aerosol particles, using all variations of Köhler theory, has largely assumed that the total solute concentrations define both the water activity and surface tension. Recently however, bulk to surface phase partitioning has been postulated as a process which significantly alters the predicted point of activation. In this paper, an analytical solution to calculate the removal of material from a bulk to a surface layer in aerosol particles has been derived using a well established and validated surface tension framework. The applicability to an unlimited number of components is possible via reliance on data from each binary system. Whilst assumptions regarding behaviour at the surface layer have been made to facilitate derivation, it is proposed that the framework presented can capture the overall impact of bulk-surface partitioning. Demonstrations of the equations for two and five component mixtures are given while comparisons are made with more detailed frameworks capable at modelling ternary systems at higher levels of complexity. Predictions made by the model across a range of surface active properties should be tested against measurements. Indeed, reccomendations are given for experimental validation and to assess sensitivities to accuracy and required level of complexity within large scale frameworks. Importantly, the computational efficiency of using the solution presented in this paper is roughly a factor of 20 less than a similar iterative approach, a comparison with highly coupled approaches not available beyond a 3 component system.

Description: The main characteristics of the new version 1.2 of the three-dimensional Earth system model of intermediate complexity LOVECLIM are briefly described. LOVECLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. The atmospheric component is ECBilt2, a T21, 3-level quasi-geostrophic model. The ocean component is CLIO3, which consists of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is of 3° by 3°, and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetation model that simulates the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the ocean carbon cycle is represented by LOCH, a comprehensive model that takes into account both the solubility and biological pumps. The ice sheet component AGISM is made up of a three-dimensional thermomechanical model of the ice sheet flow, a visco-elastic bedrock model and a model of the mass balance at the ice-atmosphere and ice-ocean interfaces. For both the Greenland and Antarctic ice sheets, calculations are made on a 10 km by 10 km resolution grid with 31 sigma levels. LOVECLIM1.2 reproduces well the major characteristics of the observed climate both for present-day conditions and for key past periods such as the last millennium, the mid-Holocene and the Last Glacial Maximum. However, despite some improvements compared to earlier versions, some biases are still present in the model. The most serious ones are mainly located at low latitudes with an overestimation of the temperature there, a too symmetric distribution of precipitation between the two hemispheres, and an overestimation of precipitation and vegetation cover in the subtropics. In addition, the atmospheric circulation is too weak. The model also tends to underestimate the surface temperature changes (mainly at low latitudes) and to overestimate the ocean heat uptake observed over the last decades.