ALBERT

Description: A new module has been implemented in the fifth generation of the ECMWF/Hamburg (ECHAM5)/Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) model that simulates cloud-related processes on a much smaller grid. This so-called superparameterisation acts as a replacement for the convection parameterisation and large-scale cloud scheme. The concept of embedding a cloud-resolving model (CRM) inside of each grid box of a general circulation model leads to an explicit representation of cloud dynamics. The new model component is evaluated against observations and the conventional usage of EMAC using a convection parameterisation. In particular, effects of applying different configurations of the superparameterisation are analysed in a systematical way. Consequences of changing the CRM's orientation, cell size and number of cells range from regional differences in cloud amount up to global impacts on precipitation distribution and its variability. For some edge case setups, the analysed climate state of superparameterised simulations even deteriorates from the mean observed energy budget. In the current model configuration, different climate regimes can be formed that are mainly driven by some of the parameters of the CRM. Presently, the simulated total cloud cover is at the lower edge of the CMIP5 model ensemble. However, certain “tuning” of the current model configuration could improve the slightly underestimated cloud cover, which will result in a shift of the simulated climate. The simulation results show that especially tropical precipitation is better represented with the superparameterisation in the EMAC model configuration. Furthermore, the diurnal cycle of precipitation is heavily affected by the choice of the CRM parameters. However, despite an improvement of the representation of the continental diurnal cycle in some configurations, other parameter choices result in a deterioration compared to the reference simulation using a conventional convection parameterisation. The ability of the superparameterisation to represent latent and sensible heat flux climatology is independent of the chosen CRM setup. Evaluation of in-atmosphere cloud amounts depending on the chosen CRM setup shows that cloud development can significantly be influenced on the large scale using a too-small CRM domain size. Therefore, a careful selection of the CRM setup is recommended using 32 or more CRM cells to compensate for computational expenses.

Description: Lightning represents one of the dominant emission sources for NOx in the troposphere. The direct release of oxidised nitrogen in the upper troposphere does not only affect ozone formation, but also chemical and microphysical properties of aerosol particles in this region. This study investigates the direct impact of LNOx emissions on upper-tropospheric nitrate using a global chemistry climate model. The simulation results show a substantial influence of the lightning emissions on the mixing ratios of nitrate aerosol in the upper troposphere of more than 50 %. In addition to the impact on nitrate, lightning substantially affects the oxidising capacity of the atmosphere with substantial implications for gas-phase sulfate formation and new particle formation in the upper troposphere. In conjunction with the condensation of nitrates, substantial differences in the aerosol size distribution occur in the upper troposphere as a consequence of lightning. This has implications for the extinction properties of the aerosol particles and for the cloud optical properties. While the extinction is generally slightly enhanced due to the LNOx emissions, the response of the clouds is ambiguous due to compensating effects in both liquid and ice clouds. Resulting shortwave flux perturbations are of   ∼ −100 mW m−2 as determined from several sensitivity scenarios, but an uncertainty range of almost 50 % has to be defined due to the large internal variability of the system and the uncertainties in the multitude of involved processes. Despite the clear statistical significance of the influence of lightning on the nitrate concentrations, the robustness of the findings gradually decreases towards the determination of the radiative flux perturbations.

Description: We present a sensitivity study on transatlantic dust transport, a process which has many implications for the atmosphere, the ocean and the climate. We investigate the impact of key processes that control the dust outflow, i.e., the emission flux, convection schemes and the chemical aging of mineral dust, by using the EMAC model following Abdelkader et al. (2015). To characterize the dust outflow over the Atlantic Ocean, we distinguish two geographic zones: (i) dust interactions within the Intertropical Convergence Zone (ITCZ), or the dust–ITCZ interaction zone (DIZ), and (ii) the adjacent dust transport over the Atlantic Ocean (DTA) zone. In the latter zone, the dust loading shows a steep and linear gradient westward over the Atlantic Ocean since particle sedimentation is the dominant removal process, whereas in the DIZ zone aerosol–cloud interactions, wet deposition and scavenging processes determine the extent of the dust outflow. Generally, the EMAC simulated dust compares well with CALIPSO observations; however, our reference model configuration tends to overestimate the dust extinction at a lower elevation and underestimates it at a higher elevation. The aerosol optical depth (AOD) over the Caribbean responds to the dust emission flux only when the emitted dust mass is significantly increased over the source region in Africa by a factor of 10. These findings point to the dominant role of dust removal (especially wet deposition) in transatlantic dust transport. Experiments with different convection schemes have indeed revealed that the transatlantic dust transport is more sensitive to the convection scheme than to the dust emission flux parameterization. To study the impact of dust chemical aging, we focus on a major dust outflow in July 2009. We use the calcium cation as a proxy for the overall chemical reactive dust fraction and consider the uptake of major inorganic acids (i.e., H2SO4, HNO3 and HCl) and their anions, i.e., sulfate (SO42−), bisulfate (HSO4−), nitrate (NO3−) and chloride (Cl−), on the surface of mineral particles. The subsequent neutralization reactions with the calcium cation form various salt compounds that cause the uptake of water vapor from the atmosphere, i.e., through the chemical aging of dust particles leading to an increase of 0.15 in the AOD under subsaturated conditions (July 2009 monthly mean). As a result of the radiative feedback on surface winds, dust emissions increased regionally. On the other hand, the aged dust particles, compared to the non-aged particles, are more efficiently removed by both wet and dry deposition due to the increased hygroscopicity and particle size (mainly due to water uptake). The enhanced removal of aged particles decreases the dust burden and lifetime, which indirectly reduces the dust AOD by 0.05 (monthly mean). Both processes can be significant (major dust outflow, July 2009), but the net effect depends on the region and level of dust chemical aging.

Description: The new submodel SVOC for the Modular Earth Submodel System (MESSy) was developed and applied within the ECHAM5/MESSy Atmospheric Chemistry (EMAC) model to simulate the atmospheric cycling and air–surface exchange processes of semivolatile organic pollutants. Our focus is on four polycyclic aromatic hydrocarbons (PAHs) of largely varying properties. Some new features in input and physics parameterizations of tracers were tested: emission seasonality, the size discretization of particulate-phase tracers, the application of poly-parameter linear free-energy relationships in gas–particle partitioning, and re-volatilization from land and sea surfaces. The results indicate that the predicted global distribution of the 3-ring PAH phenanthrene is sensitive to the seasonality of its emissions, followed by the effects of considering re-volatilization from surfaces. The predicted distributions of the 4-ring PAHs fluoranthene and pyrene and the 5-ring PAH benzo(a)pyrene are found to be sensitive to the combinations of factors with their synergistic effects being stronger than the direct effects of the individual factors. The model was validated against observations of PAH concentrations and aerosol particulate mass fraction. The annual mean concentrations are simulated to the right order of magnitude for most cases and the model well captures the species and regional variations. However, large underestimation is found over the ocean. It is found that the particulate mass fraction of the benzo(a)pyrene is well simulated, whereas those of other species are lower than observed.

Description: Climate sensitive urban planning involves the implementation of green infrastructure as one measure to mitigate excessive heat in urban areas. Depending on thermal conditions, certain trees tend to emit more biogenic volatile organic compounds, which act as precursors for ozone formation, thus hampering air quality. Combining a theoretical approach from a box model analysis and microscale modeling from the microclimate model ENVI-met, we analyze this relationship for a selected region in Germany and provide the link to air quality prediction and climate sensitive urban planning. A box model study was conducted, indicating higher ozone levels with higher isoprene concentration, especially in NO-saturated atmospheres. ENVI-met sensitivity studies showed that different urban layouts strongly determine local isoprene emissions of vegetation, with leaf temperature, rather than photosynthetic active radiation, being the dominant factor. The impact of isoprene emission on the ozone in complex urban environments was simulated for an urban area for a hot summer day with and without isoprene. A large isoprene-induced relative ozone increase was found over the whole model area. On selected hot spots we find a clear relationship between urban layout, proximity to NOx emitters, tree-species-dependent isoprene emission capacity, and increases in ozone concentration, rising up to 500% locally.

Description: The Modular Earth Submodel System (MESSy) provides an interface to couple submodels to a base model via a highly flexible data management facility (Jöckel et al., 2010). In the present paper we present the four new radiation related submodels RAD, AEROPT, CLOUDOPT, and ORBIT. The submodel RAD (including the shortwave radiation scheme RAD_FUBRAD) simulates the radiative transfer, the submodel AEROPT calculates the aerosol optical properties, the submodel CLOUDOPT calculates the cloud optical properties, and the submodel ORBIT is responsible for Earth orbit calculations. These submodels are coupled via the standard MESSy infrastructure and are largely based on the original radiation scheme of the general circulation model ECHAM5, however, expanded with additional features. These features comprise, among others, user-friendly and flexibly controllable (by namelists) online radiative forcing calculations by multiple diagnostic calls of the radiation routines. With this, it is now possible to calculate radiative forcing (instantaneous as well as stratosphere adjusted) of various greenhouse gases simultaneously in only one simulation, as well as the radiative forcing of cloud perturbations. Examples of online radiative forcing calculations in the ECHAM/MESSy Atmospheric Chemistry (EMAC) model are presented.

Description: Enhanced aerosol abundance in the upper troposphere and lower stratosphere (UTLS) associated with the Asian summer monsoon (ASM) is referred to as the Asian Tropopause Aerosol Layer (ATAL). The chemical composition, microphysical properties, and climate effects of aerosols in the ATAL have been the subject of discussion over the past decade. In this work, we use the ECHAM/MESSy Atmospheric Chemistry (EMAC) general circulation model at a relatively fine grid resolution (about 1.1×1.1∘) to numerically simulate the emissions, chemistry, and transport of aerosols and their precursors in the UTLS within the ASM anticyclone during the years 2010–2012. We find a pronounced maximum of aerosol extinction in the UTLS over the Tibetan Plateau, which to a large extent is caused by mineral dust emitted from the northern Tibetan Plateau and slope areas, lofted to an altitude of at least 10 km, and accumulating within the anticyclonic circulation. We also find that the emissions and convection of ammonia in the central main body of the Tibetan Plateau make a great contribution to the enhancement of gas-phase NH3 in the UTLS over the Tibetan Plateau and ASM anticyclone region. Our simulations show that mineral dust, water-soluble compounds, such as nitrate and sulfate, and associated liquid water dominate aerosol extinction in the UTLS within the ASM anticyclone. Due to shielding of high background sulfate concentrations outside the anticyclone from volcanoes, a relative minimum of aerosol extinction within the anticyclone in the lower stratosphere is simulated, being most pronounced in 2011, when the Nabro eruption occurred. In contrast to mineral dust and nitrate concentrations, sulfate increases with increasing altitude due to the larger volcano effects in the lower stratosphere compared to the upper troposphere. Our study indicates that the UTLS over the Tibetan Plateau can act as a well-defined conduit for natural and anthropogenic gases and aerosols into the stratosphere.

Description: Hydroxyl radical (OH) plays critical roles within the troposphere, such as determining the lifetime of methane (CH4), yet is challenging to model due to its fast cycling and dependence on a multitude of sources and sinks. As a result, the reasons for variations in OH and the resulting CH4 lifetime (τCH4), both between models and in time, are difficult to diagnose. We apply a neural network (NN) approach to address this issue within a group of models that participated in the Chemistry-Climate Model Initiative (CCMI). Analysis of the historical specified dynamics simulations performed for CCMI indicates that the primary drivers of τCH4 differences among ten models are the flux of UV light to the troposphere (indicated by the photolysis frequency JO1D) due mostly to clouds, mixing ratio of tropospheric ozone (O3), the abundance of nitrogen oxides (NOx≡NO+NO2), and details of the various chemical mechanisms that drive OH. Water vapor, carbon monoxide (CO), the ratio of NO:NOx, and formaldehyde (HCHO) explain moderate differences in τCH4, while isoprene, CH4, the photolysis frequency of NO2 by visible light (JNO2), overhead O3 column, and temperature account for little-to-no model variation in τCH4. We also apply the NNs to analysis of temporal trends in OH from 1980 to 2015. All models that participated in the specified dynamics historical simulation for CCMI demonstrate a decline in τCH4 during the analysed timeframe. The significant contributors to this trend, in order of importance, are tropospheric O3, JO1D, NOx, and H2O, with CO also causing substantial interannual variability in OH burden. Finally, the identified trends in τCH4 are compared to calculated trends in the tropospheric mean OH concentration from previous work, based on analysis of observations. The comparison reveals a robust result for the effect of rising water vapor on OH and τCH4, imparting an increasing and decreasing trend of about 0.5 % decade−1, respectively. The responses due to NOx, O3 column, and temperature are also in reasonably good agreement between the two studies, though a discrepancy in the CH4 response highlights a need for further examination of the CH4 feedback on the abundance of OH.

Description: A new module has been implemented in the ECHAM5/MESSy Atmospheric Chemistry (EMAC) Model that simulates cloud related processes on a much smaller grid. This so called superparametrisation acts as a replacement for the convection parametrisation and large-scale cloud scheme. The concept of embedding an ensemble of cloud resolving models (CRMs) inside of each grid box of a general circulation model leads to an explicit representation of cloud dynamics. The new model component is evaluated against observations and the conventional usage of EMAC using a convection parametrisation. In particular, effects of applying different configurations of the superparametrisation are analyzed in a systematical way. Consequences of changing the CRMs orientation, cell size and number of cells range from regional differences in cloud amount up to global impacts on precipitation distribution and its variability. For some edge case setups the analysed climate state of superparametrised simulations even deteriorates from the mean observed energy budget. In the current model configuration different climate regimes can be formed that are mainly driven by some of the parameters of the CRM. Presently, the simulated cloud cover is at the lower edge of the CMIP5 model ensemble indicating that the hydrological overturning is too efficient. However, certain "tuning" of the current model configuration could improve the currently underestimated cloud cover, which will result in a shift of the climate. The simulation results show that especially tropical precipitation is better represented with the superparamerisation in the EMAC model configuration. Furthermore, the diurnal cycle of precipitation is heavily affected by the choice of the CRM parameters. However, despite an improvement of the representation of the continental diurnal cycle in some configurations, other parameter choices result in a deterioration compared to the reference simulation using a conventional convection parameterisation. The ability of the superparametrisation to represent latent and sensible heat flux climatology is dependent on the chosen CRM setup. Further interactions of the planetary boundary layer and the free troposphere can significantly influence cloud development on the large-scale. Therefore a careful selection of the CRM setup is recommended to compensate for computational expenses.

Description: The hydroxyl radical (OH) plays critical roles within the troposphere, such as determining the lifetime of methane (CH4), yet is challenging to model due to its fast cycling and dependence on a multitude of sources and sinks. As a result, the reasons for variations in OH and the resulting methane lifetime (τCH4), both between models and in time, are difficult to diagnose. We apply a neural network (NN) approach to address this issue within a group of models that participated in the Chemistry-Climate Model Initiative (CCMI). Analysis of the historical specified dynamics simulations performed for CCMI indicates that the primary drivers of τCH4 differences among 10 models are the flux of UV light to the troposphere (indicated by the photolysis frequency JO1D), the mixing ratio of tropospheric ozone (O3), the abundance of nitrogen oxides (NOx≡NO+NO2), and details of the various chemical mechanisms that drive OH. Water vapour, carbon monoxide (CO), the ratio of NO:NOx, and formaldehyde (HCHO) explain moderate differences in τCH4, while isoprene, methane, the photolysis frequency of NO2 by visible light (JNO2), overhead ozone column, and temperature account for little to no model variation in τCH4. We also apply the NNs to analysis of temporal trends in OH from 1980 to 2015. All models that participated in the specified dynamics historical simulation for CCMI demonstrate a decline in τCH4 during the analysed timeframe. The significant contributors to this trend, in order of importance, are tropospheric O3, JO1D, NOx, and H2O, with CO also causing substantial interannual variability in OH burden. Finally, the identified trends in τCH4 are compared to calculated trends in the tropospheric mean OH concentration from previous work, based on analysis of observations. The comparison reveals a robust result for the effect of rising water vapour on OH and τCH4, imparting an increasing and decreasing trend of about 0.5 % decade−1, respectively. The responses due to NOx, ozone column, and temperature are also in reasonably good agreement between the two studies.