scholarly journals A new radiation infrastructure for the Modular Earth Submodel System (MESSy, based on version 2.51)

2016 ◽  
Vol 9 (6) ◽  
pp. 2209-2222 ◽  
Author(s):  
Simone Dietmüller ◽  
Patrick Jöckel ◽  
Holger Tost ◽  
Markus Kunze ◽  
Catrin Gellhorn ◽  
...  

Abstract. 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.

2016 ◽  
Author(s):  
Simone Dietmüller ◽  
Patrick Jöckel ◽  
Holger Tost ◽  
Markus Kunze ◽  
Cathrin Gellhorn ◽  
...  

Abstract. The Modular Earth Submodel System (MESSy) provides an interface to couple submodels to a basemodel 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 (with 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) on-line 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 on-line radiative forcing calculations in the ECHAM/MESSy Atmospheric Chemistry (EMAC) model are presented.


2014 ◽  
Vol 7 (5) ◽  
pp. 2503-2516 ◽  
Author(s):  
K. Klingmüller ◽  
B. Steil ◽  
C. Brühl ◽  
H. Tost ◽  
J. Lelieveld

Abstract. The modelling of aerosol radiative forcing is a major cause of uncertainty in the assessment of global and regional atmospheric energy budgets and climate change. One reason is the strong dependence of the aerosol optical properties on the mixing state of aerosol components, such as absorbing black carbon and, predominantly scattering sulfates. Using a new column version of the aerosol optical properties and radiative-transfer code of the ECHAM/MESSy atmospheric-chemistry–climate model (EMAC), we study the radiative transfer applying various mixing states. The aerosol optics code builds on the AEROPT (AERosol OPTical properties) submodel, which assumes homogeneous internal mixing utilising the volume average refractive index mixing rule. We have extended the submodel to additionally account for external mixing, partial external mixing and multilayered particles. Furthermore, we have implemented the volume average dielectric constant and Maxwell Garnett mixing rule. We performed regional case studies considering columns over China, India and Africa, corroborating much stronger absorption by internal than external mixtures. Well-mixed aerosol is a good approximation for particles with a black-carbon core, whereas particles with black carbon at the surface absorb significantly less. Based on a model simulation for the year 2005, we calculate that the global aerosol direct radiative forcing for homogeneous internal mixing differs from that for external mixing by about 0.5 W m−2.


1998 ◽  
Vol 11 (8) ◽  
pp. 1883-1905 ◽  
Author(s):  
O. P. Sharma ◽  
H. Le Treut ◽  
G. Sèze ◽  
L. Fairhead ◽  
R. Sadourny

Abstract The sensitivity of the interannual variations of the summer monsoons to imposed cloudiness has been studied with a general circulation model using the initial conditions prepared from the European Centre for Medium-Range Forecasts analyses of 1 May 1987 and 1988. The cloud optical properties in this global model are calculated from prognostically computed cloud liquid water. The model successfully simulates the contrasting behavior of these two successive monsoons. However, when the optical properties of the observed clouds are specified in the model runs, the simulations show some degradation over India and its vicinity. The main cause of this degradation is the reduced land–sea temperature contrast resulting from the radiative effects of the observed clouds imposed in such simulations. It is argued that the high concentration of condensed water content of clouds over the Indian land areas will serve to limit heating of the land, thereby reducing the thermal contrast that gives rise to a weak Somali jet. A countermonsoon circulation is, therefore, simulated in the vector difference field of 850-hPa winds from the model runs with externally specified clouds. This countermonsoon circulation is associated with an equatorial heat source that is the response of the model to the radiative effects of the imposed clouds. Indeed, there are at least two clear points that can be made: 1) the cloud–SST patterns, together, affect the interannual variability; and 2) with both clouds and SST imposed, the model simulation is less sensitive to initial conditions. Additionally, the study emphasizes the importance of dynamically consistent clouds developing in response to the dynamical, thermal, and moist state of the atmosphere during model integrations.


2007 ◽  
Vol 7 (10) ◽  
pp. 2503-2515 ◽  
Author(s):  
C. Cagnazzo ◽  
E. Manzini ◽  
M. A. Giorgetta ◽  
P. M. De F. Forster ◽  
J. J. Morcrette

Abstract. In order to improve the representation of ozone absorption in the stratosphere of the MAECHAM5 general circulation model, the spectral resolution of the shortwave radiation parameterization used in the model has been increased from 4 to 6 bands. Two 20-years simulations with the general circulation model have been performed, one with the standard and the other with the newly introduced parameterization respectively, to evaluate the temperature and dynamical changes arising from the two different representations of the shortwave radiative transfer. In the simulation with the increased spectral resolution in the radiation parameterization, a significant warming of almost the entire model domain is reported. At the summer stratopause the temperature increase is about 6 K and alleviates the cold bias present in the model when the standard radiation scheme is used. These general circulation model results are consistent both with previous validation of the radiation scheme and with the offline clear-sky comparison performed in the current work with a discrete ordinate 4 stream scattering line by line radiative transfer model. The offline validation shows a substantial reduction of the daily averaged shortwave heating rate bias (1–2 K/day cooling) that occurs for the standard radiation parameterization in the upper stratosphere, present under a range of atmospheric conditions. Therefore, the 6 band shortwave radiation parameterization is considered to be better suited for the representation of the ozone absorption in the stratosphere than the 4 band parameterization. Concerning the dynamical response in the general circulation model, it is found that the reported warming at the summer stratopause induces stronger zonal mean zonal winds in the middle atmosphere. These stronger zonal mean zonal winds thereafter appear to produce a dynamical feedback that results in a dynamical warming (cooling) of the polar winter (summer) mesosphere, caused by an increased downward (upward) circulation in the winter (summer) hemisphere. In addition, the comparison of the two simulations performed with the general circulation model shows that the increase in the spectral resolution of the shortwave radiation and the associated changes in the cloud optical properties result in a warming (0.5–1 K) and moistening (3%–12%) of the upper tropical troposphere. By comparing these modeled differences with previous works, it appears that the reported changes in the solar radiation scheme contribute to improve the model mean temperature also in the troposphere.


2011 ◽  
Vol 29 (5) ◽  
pp. 955-963 ◽  
Author(s):  
D. Goto ◽  
K. V. S. Badarinath ◽  
T. Takemura ◽  
T. Nakajima

Abstract. Aerosols have great impacts on atmospheric environment, human health, and earth's climate. Therefore, information on their spatial and temporal distribution is of paramount importance. Despite numerous studies have examined the variation and trends of BC and AOD over India, only very few have focused on their spatial distribution or even correlating the observations with model simulations. In the present study, a three-dimensional aerosol transport-radiation model coupled with a general circulation model. SPRINTARS, simulated atmospheric aerosol distributions including BC and aerosol optical properties, i.e., aerosol optical thickness (AOT), Ångström Exponent (AE), and single scattering albedo (SSA). The simulated results are compared with both BC measurements by aethalometer and aerosol optical properties measured by ground-based skyradiometer and by satellite sensor, MODIS/Terra over Hyderabad, which is a tropical urban area of India, for the year 2008. The simulated AOT and AE in Hyderabad are found to be comparable to ground-based measured ones. The simulated SSA tends to be higher than the ground-based measurements. Both these comparisons of aerosol optical properties between the simulations with different emission inventories and the measurements indicate that, firstly the model uncertainties derived from aerosol emission inventory cannot explain the gaps between the simulations and the measurements and secondly the vertical transport of BC and the treatment of BC-containing particles can be the main issue in the global model to solve the gap.


2006 ◽  
Vol 6 (6) ◽  
pp. 11067-11092 ◽  
Author(s):  
C. Cagnazzo ◽  
E. Manzini ◽  
M. A. Giorgetta ◽  
P. M. De F. Forster

Abstract. In order to improve the representation of the shortwave radiative transfer in the MAECHAM5 general circulation model, the spectral resolution of the shortwave radiation parameterization used in the model has been increased and extended in the UV-B and UV-C bands. The upgraded shortwave parameterization is first validated offline with a 4 stream discrete-ordinate line-by-line model. Thereafter, two 20-years simulations with the MAECHAM5 middle atmosphere general circulation model are performed to evaluate the temperature changes and the dynamical feedbacks arising from the newly introduced parameterization. The offline clear-sky comparison of the standard and upgraded parameterizations with the discrete ordinate model shows considerable improvement for the upgraded parameterization in terms of shortwave fluxes and heating rates. In the simulation with the upgraded ratiation parameterization, we report a significant warming of almost the entire atmosphere, largest at 1 hPa at the stratopause, and stronger zonal mean zonal winds in the middle atmosphere. The warming at the summer stratopause alleviates the cold bias present in the model when the standard radiation scheme is used. The stronger zonal mean zonal winds induce a dynamical feedback that results in a dynamical warming (cooling) of the polar winter (summer) mesosphere, caused by an increased downward (upward) circulation in the winter (summer) hemisphere. In the troposphere, the changes in the spectral resolution and the associated changes in the cloud optical parameters introduce a relatively small warming and, consistenly, a moisteneing. The warming occurs mostly in the upper troposphere and can contribute to a possible improvement of the model temperature climatology.


2013 ◽  
Vol 13 (4) ◽  
pp. 11395-11425 ◽  
Author(s):  
C. Brühl ◽  
J. Lelieveld ◽  
M. Höpfner ◽  
H. Tost

Abstract. A multiyear study with the atmospheric chemistry general circulation model EMAC with the aerosol module GMXe at high altitude resolution demonstrates that the sulfur gases COS and SO2, the latter from low-latitude volcanic eruptions, predominantly control the formation of stratospheric aerosol. The model consistently uses the same parameters in the troposphere and stratosphere for 7 aerosol modes applied. Lower boundary conditions for COS and other long-lived trace gases are taken from measurement networks, while estimates of volcanic SO2 emissions are based on satellite observations. We show comparisons with satellite data for aerosol extinction (e.g. SAGE) and SO2 in the middle atmosphere (MIPAS on ENVISAT). This corroborates the interannual variability induced by the Quasi-Biennial Oscillation, which is internally generated by the model. The model also realistically simulates the radiative effects of stratospheric and tropospheric aerosol including the effects on the model dynamics. The medium strength volcanic eruptions of 2005 and 2006 exerted a nonnegligible radiative forcing of up to −0.6 W m−2 in the tropics, while the large Pinatubo eruption caused a maximum though short term tropical forcing of about −10 W m−2. The study also shows that observed upper stratospheric SO2 can be simulated accurately only when a sulphur sink on meteoritic dust is included and the photolysis of gaseous H2SO4 in the near infrared is higher than assumed previously.


2014 ◽  
Vol 7 (3) ◽  
pp. 3367-3402
Author(s):  
K. Klingmüller ◽  
B. Steil ◽  
C. Brühl ◽  
H. Tost ◽  
J. Lelieveld

Abstract. The modelling of aerosol radiative forcing is a major cause of uncertainty in the assessment of global and regional atmospheric energy budgets and climate change. One reason is the strong dependence of the aerosol optical properties on the mixing state of aerosol components like black carbon and sulphates. Using a new column version of the aerosol optical properties and radiative transfer code of the atmospheric chemistry-climate model EMAC, we study the radiative transfer applying various mixing states. The aerosol optics code builds on the AEROPT submodel which assumes homogeneous internal mixing utilising the volume average refractive index mixing rule. We have extended the submodel to additionally account for external mixing, partial external mixing and multi-layered particles. Furthermore, we have implemented the volume average dielectric-constant and Maxwell Garnett Mixing rule. We performed regional case studies considering columns over China, India and Africa, corroborating much stronger absorption by internal than external mixtures. Well mixed aerosol is a good approximation for particles with a black carbon core, whereas particles with black carbon at the surface absorb significantly less. Based on a model simulation for the year 2005 we calculate that the global aerosol direct radiative-forcing for homogeneous internal mixing differs from that for external mixing by about 0.5 W m−2.


2018 ◽  
Vol 11 (3) ◽  
pp. 1059-1076 ◽  
Author(s):  
Astrid Kerkweg ◽  
Christiane Hofmann ◽  
Patrick Jöckel ◽  
Mariano Mertens ◽  
Gregor Pante

Abstract. As part of the Modular Earth Submodel System (MESSy), the Multi-Model-Driver (MMD v1.0) was developed to couple online the regional Consortium for Small-scale Modeling (COSMO) model into a driving model, which can be either the regional COSMO model or the global European Centre Hamburg general circulation model (ECHAM) (see Part 2 of the model documentation). The coupled system is called MECO(n), i.e., MESSy-fied ECHAM and COSMO models nested n times. In this article, which is part of the model documentation of the MECO(n) system, the second generation of MMD is introduced. MMD comprises the message-passing infrastructure required for the parallel execution (multiple programme multiple data, MPMD) of different models and the communication of the individual model instances, i.e. between the driving and the driven models. Initially, the MMD library was developed for a one-way coupling between the global chemistry–climate ECHAM/MESSy atmospheric chemistry (EMAC) model and an arbitrary number of (optionally cascaded) instances of the regional chemistry–climate model COSMO/MESSy. Thus, MMD (v1.0) provided only functions for unidirectional data transfer, i.e. from the larger-scale to the smaller-scale models.Soon, extended applications requiring data transfer from the small-scale model back to the larger-scale model became of interest. For instance, the original fields of the larger-scale model can directly be compared to the upscaled small-scale fields to analyse the improvements gained through the small-scale calculations, after the results are upscaled. Moreover, the fields originating from the two different models might be fed into the same diagnostic tool, e.g. the online calculation of the radiative forcing calculated consistently with the same radiation scheme. Last but not least, enabling the two-way data transfer between two models is the first important step on the way to a fully dynamical and chemical two-way coupling of the various model instances.In MMD (v1.0), interpolation between the base model grids is performed via the COSMO preprocessing tool INT2LM, which was implemented into the MMD submodel for online interpolation, specifically for mapping onto the rotated COSMO grid. A more flexible algorithm is required for the backward mapping. Thus, MMD (v2.0) uses the new MESSy submodel GRID for the generalised definition of arbitrary grids and for the transformation of data between them.In this article, we explain the basics of the MMD expansion and the newly developed generic MESSy submodel GRID (v1.0) and show some examples of the abovementioned applications.


2020 ◽  
Vol 20 (1) ◽  
pp. 281-301 ◽  
Author(s):  
Le Kuai ◽  
Kevin W. Bowman ◽  
Kazuyuki Miyazaki ◽  
Makoto Deushi ◽  
Laura Revell ◽  
...  

Abstract. The top-of-atmosphere (TOA) outgoing longwave flux over the 9.6 µm ozone band is a fundamental quantity for understanding chemistry–climate coupling. However, observed TOA fluxes are hard to estimate as they exhibit considerable variability in space and time that depend on the distributions of clouds, ozone (O3), water vapor (H2O), air temperature (Ta), and surface temperature (Ts). Benchmarking present-day fluxes and quantifying the relative influence of their drivers is the first step for estimating climate feedbacks from ozone radiative forcing and predicting radiative forcing evolution. To that end, we constructed observational instantaneous radiative kernels (IRKs) under clear-sky conditions, representing the sensitivities of the TOA flux in the 9.6 µm ozone band to the vertical distribution of geophysical variables, including O3, H2O, Ta, and Ts based upon the Aura Tropospheric Emission Spectrometer (TES) measurements. Applying these kernels to present-day simulations from the Chemistry-Climate Model Initiative (CCMI) project as compared to a 2006 reanalysis assimilating satellite observations, we show that the models have large differences in TOA flux, attributable to different geophysical variables. In particular, model simulations continue to diverge from observations in the tropics, as reported in previous studies of the Atmospheric Chemistry Climate Model Intercomparison Project (ACCMIP) simulations. The principal culprits are tropical middle and upper tropospheric ozone followed by tropical lower tropospheric H2O. Five models out of the eight studied here have TOA flux biases exceeding 100 mW m−2 attributable to tropospheric ozone bias. Another set of five models have flux biases over 50 mW m−2 due to H2O. On the other hand, Ta radiative bias is negligible in all models (no more than 30 mW m−2). We found that the atmospheric component (AM3) of the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model and Canadian Middle Atmosphere Model (CMAM) have the lowest TOA flux biases globally but are a result of cancellation of opposite biases due to different processes. Overall, the multi-model ensemble mean bias is -133±98 mW m−2, indicating that they are too atmospherically opaque due to trapping too much radiation in the atmosphere by overestimated tropical tropospheric O3 and H2O. Having too much O3 and H2O in the troposphere would have different impacts on the sensitivity of TOA flux to O3 and these competing effects add more uncertainties on the ozone radiative forcing. We find that the inter-model TOA outgoing longwave radiation (OLR) difference is well anti-correlated with their ozone band flux bias. This suggests that there is significant radiative compensation in the calculation of model outgoing longwave radiation.


Sign in / Sign up

Export Citation Format

Share Document