European Centre Hamburg Research Articles

Formaldehyde and hydroperoxide distribution around the Arabian Peninsula – evaluation of EMAC model results with ship-based measurements

Abstract. Formaldehyde (HCHO), hydrogen peroxide (H2O2) and organic hydroperoxides (ROOH) play a key role in atmospheric oxidation processes. They act as sources and sinks for HOx radicals (OH + HO2), with OH as the primary oxidant that governs the atmospheric self-cleaning capacity. Measurements of these species allow for evaluation of chemistry-transport models which need to account for multifarious source distributions, transport, complex photochemical reaction pathways and deposition processes of these species. HCHO is an intermediate during the oxidation of volatile organic compounds (VOCs) and is an indicator of photochemical activity and combustion-related emissions. In this study, we use in situ observations of HCHO, H2O2 and ROOH in the marine boundary layer (MBL) to evaluate results of the general circulation model EMAC (ECHAM5/MESSy2 Atmospheric Chemistry; European Center HAMburg, Modular Earth Submodel System). The dataset was obtained during the Air Quality and Climate Change in the Arabian Basin (AQABA) ship campaign around the Arabian Peninsula in summer 2017. This region is characterized by high levels of photochemical air pollution, humidity and solar irradiation, especially in the areas around the Suez Canal and the Arabian Gulf. High levels of air pollution with up to 12 ppbv HCHO, 2.3 ppbv ROOH and relatively low levels of H2O2 (≤0.5 ppbv) were detected over the Arabian Gulf. We find that EMAC failed to predict absolute mixing ratios of HCHO and ROOH during high-pollution events over the Arabian Gulf, while it reproduced HCHO on average within a factor of 2. Dry deposition velocities were determined for HCHO and H2O2 at night with 0.77±0.29 cm s−1 for HCHO and 1.03±0.52 cm s−1 for H2O2 over the Arabian Sea, which were matched by EMAC. The photochemical budget of H2O2 revealed elevated HOx radical concentrations in EMAC, which resulted in an overestimation of H2O2 by more than a factor of 5 for the AQABA dataset. The underestimated air pollution over the Arabian Gulf was related to EMAC's coarse spatial resolution and missing anthropogenic emissions in the model.

Tropospheric ozone production and chemical regime analysis during the COVID-19 lockdown over Europe

Abstract. The COVID-19 (coronavirus disease 2019) European lockdowns have led to a significant reduction in the emissions of primary pollutants such as NO (nitric oxide) and NO2 (nitrogen dioxide). As most photochemical processes are related to nitrogen oxide (NOx≡ NO + NO2) chemistry, this event has presented an exceptional opportunity to investigate its effects on air quality and secondary pollutants, such as tropospheric ozone (O3). In this study, we present the effects of the COVID-19 lockdown on atmospheric trace gas concentrations, net ozone production rates (NOPRs) and the dominant chemical regime throughout the troposphere based on three different research aircraft campaigns across Europe. These are the UTOPIHAN (Upper Tropospheric Ozone: Processes Involving HOx and NOx) campaigns in 2003 and 2004, the HOOVER (HOx over Europe) campaigns in 2006 and 2007, and the BLUESKY campaign in 2020, the latter performed during the COVID-19 lockdown. We present in situ observations and simulation results from the ECHAM5 (fifth-generation European Centre Hamburg general circulation model, version 5.3.02)/MESSy2 (second-generation Modular Earth Submodel System, version 2.54.0) Atmospheric Chemistry (EMAC), model which allows for scenario calculations with business-as-usual emissions during the BLUESKY campaign, referred to as the “no-lockdown scenario”. We show that the COVID-19 lockdown reduced NO and NO2 mixing ratios in the upper troposphere by around 55 % compared to the no-lockdown scenario due to reduced air traffic. O3 production and loss terms reflected this reduction with a deceleration in O3 cycling due to reduced mixing ratios of NOx, while NOPRs were largely unaffected. We also study the role of methyl peroxyradicals forming HCHO (αCH3O2) to show that the COVID-19 lockdown shifted the chemistry in the upper-troposphere–tropopause region to a NOx-limited regime during BLUESKY. In comparison, we find a volatile organic compound (VOC)-limited regime to be dominant during UTOPIHAN.

Challenge of modelling GLORIA observations of upper troposphere–lowermost stratosphere trace gas and cloud distributions at high latitudes: a case study with state-of-the-art models

Abstract. Water vapour and ozone are important for the thermal and radiative balance of the upper troposphere (UT) and lowermost stratosphere (LMS). Both species are modulated by transport processes. Chemical and microphysical processes affect them differently. Thus, representing the different processes and their interactions is a challenging task for dynamical cores, chemical modules and microphysical parameterisations of state-of-the-art atmospheric model components. To test and improve the models, high-resolution measurements of the UT–LMS are required. Here, we use measurements taken in a flight of the GLORIA (Gimballed Limb Observer for Radiance Imaging of the Atmosphere) instrument on HALO (High Altitude and LOng Range Research Aircraft). The German research aircraft HALO performed a research flight on 26 February 2016 that covered deeply subsided air masses of the aged 2015/16 Arctic vortex, high-latitude LMS air masses, a highly textured region affected by troposphere-to-stratosphere exchange and high-altitude cirrus clouds. Therefore, it provides a challenging multifaceted case study for comparing GLORIA observations with state-of-the-art atmospheric model simulations in a complex UT–LMS region at a late stage of the Arctic winter 2015/16. Using GLORIA observations in this manifold scenario, we test the ability of the numerical weather prediction (NWP) model ICON (ICOsahedral Nonhydrostatic) with the extension ART (Aerosols and Reactive Trace gases) and the chemistry–climate model (CCM) EMAC (ECHAM5/MESSy Atmospheric Chemistry – fifth-generation European Centre Hamburg general circulation model/Modular Earth Submodel System) to model the UT–LMS composition of water vapour (H2O), ozone (O3), nitric acid (HNO3) and clouds. Within the scales resolved by the respective model, we find good overall agreement of both models with GLORIA. The applied high-resolution ICON-ART set-up involving an R2B7 nest (local grid refinement with a horizontal resolution of about 20 km), covering the HALO flight region, reproduces mesoscale dynamical structures well. Narrow moist filaments in the LMS observed by GLORIA at tropopause gradients in the context of a Rossby wave breaking event and in the vicinity of an occluded Icelandic low are clearly reproduced by the model. Using ICON-ART, we show that a larger filament in the west was transported horizontally into the Arctic LMS in connection with a jet stream split associated with poleward breaking of a cyclonically sheared Rossby wave. Further weaker filaments are associated with an older tropopause fold in the east. Given the lower resolution (T106) of the nudged simulation of the EMAC model, we find that this model also reproduces these features well. Overall, trace gas mixing ratios simulated by both models are in a realistic range, and major cloud systems observed by GLORIA are mostly reproduced. However, we find both models to be affected by a well-known systematic moist bias in the LMS. Further biases are diagnosed in the ICON-ART O3, EMAC H2O and EMAC HNO3 distributions. Finally, we use sensitivity simulations to investigate (i) short-term cirrus cloud impacts on the H2O distribution (ICON-ART), (ii) the overall impact of polar winter chemistry and microphysical processing on O3 and HNO3 (ICON-ART and EMAC), (iii) the impact of the model resolution on simulated parameters (EMAC), and (iv) consequences of scavenging processes by cloud particles (EMAC). We find that changing the horizontal model resolution results in notable systematic changes for all species in the LMS, while scavenging processes play a role only in the case of HNO3. We discuss the model biases and deficits found in this case study that potentially affect forecasts and projections (adversely) and provide suggestions for further model improvements.

A parameterization of long-continuing-current (LCC) lightning in the lightning submodel LNOX (version 3.0) of the Modular Earth Submodel System (MESSy, version 2.54)

Abstract. Lightning flashes can produce a discharge in which a continuing electrical current flows for more than 40 ms. Such flashes are proposed to be the main precursors of lightning-ignited wildfires and also to trigger sprite discharges in the mesosphere. However, lightning parameterizations implemented in global atmospheric models do not include information about the continuing electrical current of flashes. The continuing current of lightning flashes cannot be detected by conventional lightning location systems. Instead, these so-called long-continuing-current (LCC) flashes are commonly observed by extremely low-frequency (ELF) sensors and by optical instruments located in space. Reports of LCC lightning flashes tend to occur in winter and oceanic thunderstorms, which suggests a connection between weak convection and the occurrence of this type of discharge. In this study, we develop a parameterization of LCC lightning flashes based on a climatology derived from optical lightning measurements reported by the Lightning Imaging Sensor (LIS) on board the International Space Station (ISS) between March 2017 and March 2020. We use meteorological data from reanalyses to develop a global parameterization that uses vertical velocity at the 450 hPa pressure level as a proxy for the ratio of LCC to total lightning in thunderstorms. We implement this parameterization into the LNOX submodel of the Modular Earth Submodel System (MESSy) for usage within the European Center HAMburg general circulation model (ECHAM)/MESSy Atmospheric Chemistry (EMAC) model and compare the observed and simulated climatologies of LCC lightning flashes using six different lightning parameterizations. We find that the best agreement between the simulated and observed spatial distribution is obtained when using a novel combined lightning parameterization based on the cloud-top height over land and on the convective precipitation over ocean.

Role of convection\u2013circulation coupling in the propagation mechanism of the Madden\u2013Julian Oscillation over the Maritime Continent in a climate model

This study investigates the role of convection–circulation coupling on the simulated eastward propagation of the Madden–Julian Oscillation (MJO) over the Maritime Continent (MC). Experiments are conducted with the European Centre Hamburg Model Version 5 (ECHAM5) coupled with the one-column ocean model—Snow-Ice-Thermocline (SIT) and two different cumulus schemes, Nordeng-Tiedtke (E5SIT-Nord) and Tiedtke (E5SIT-Tied). During the early phase of MJO composites, the E5SIT-Nord simulation reveals stronger intraseasonal anomalies in the apparent heat source (Q1) over the convective center, however, the E5SIT-Tied produces a stronger background Q1, suggesting that deep convection prevails over the MC but does not couple with the MJO circulation. Similarly, in the E5SIT-Tied simulation, in-column moisture is kept mostly by local deep convection over the MC, which is in contrast to the well-correlated relationship between moisture anomaly and MJO circulation in E5SIT-Nord. A case study based on an observational MJO reveals similar biases concerning of convection–circulation coupling emerges within a few days of simulations. The E5SIT-Tied simulation produces weaker heating at the convective center of the MJO than the E5SIT-Nord a few days after model initiation, resulting weaker subsidence to the east and less favorable for propagation. The present findings highlight the instantaneous responses of cumulus parameterization schemes to MJO-related environmental changes can further affect intraseasonal variability through altering convection–circulation coupling over the MC. Physical schemes of moist convection are essential to realistically represent this coupling and thereby improve the simulation of the eastward propagation of the MJO.

Analysis of Aircraft Routing Strategies for North Atlantic Flights by Using AirTraf 2.0

Climate-optimized routing is an operational measure to effectively reduce the climate impact of aviation with a slight increase in aircraft operating costs. This study examined variations in the flight characteristics among five aircraft routing strategies and discusses several characteristics of those routing strategies concerning typical weather conditions over the North Atlantic. The daily variability in the North Atlantic weather patterns was analyzed by using the European Center Hamburg general circulation model (ECHAM) and the Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) model in the specified dynamics mode from December 2008 to August 2018. All days of the ten complete winters and summers in the simulations were classified into five weather types for winter and into three types for summer. The obtained frequency for each of the weather types was in good agreement with the literature data; and then representative days for each weather type were selected. Moreover, a total of 103 North Atlantic flights of an Airbus A330 aircraft were simulated with five aircraft routing strategies for each representative day by using the EMAC model with the air traffic simulation submodel AirTraf. For every weather type, climate-optimized routing shows the lowest climate impact, at which a trade-off exists between the operating costs and the climate impact. Cost-optimized routing lies between the time- and fuel-optimized routings and achieves the lowest operating costs by taking the best compromise between flight time and fuel use. The aircraft routing for contrail avoidance shows the second lowest climate impact; however, this routing causes extra operating costs. Our methodology could be extended to statistical analysis based on long-term simulations to clarify the relationship between the aircraft routing characteristics and weather conditions.

A revised dry deposition scheme for land–atmosphere exchange of trace gases in ECHAM/MESSy v2.54

Abstract. Dry deposition to vegetation is a major sink of ground-level ozone and is responsible for about 20 % of the total tropospheric ozone loss. Its parameterization in atmospheric chemistry models represents a significant source of uncertainty for the global tropospheric ozone budget and might account for the mismatch with observations. The model used in this study, the Modular Earth Submodel System version 2 (MESSy2) linked to the fifth-generation European Centre Hamburg general circulation model (ECHAM5) as an atmospheric circulation model (EMAC), is no exception. Like many global models, EMAC employs a “resistance in series” scheme with the major surface deposition via plant stomata which is hardly sensitive to meteorology, depending only on solar radiation. Unlike many global models, however, EMAC uses a simplified high resistance for non-stomatal deposition which makes this pathway negligible in the model. However, several studies have shown this process to be comparable in magnitude to the stomatal uptake, especially during the night over moist surfaces. Hence, we present here a revised dry deposition in EMAC including meteorological adjustment factors for stomatal closure and an explicit cuticular pathway. These modifications for the three stomatal stress functions have been included in the newly developed MESSy VERTEX submodel, i.e. a process model describing the vertical exchange in the atmospheric boundary layer, which will be evaluated for the first time here. The scheme is limited by a small number of different surface types and generalized parameters. The MESSy submodel describing the dry deposition of trace gases and aerosols (DDEP) has been revised accordingly. The comparison of the simulation results with measurement data at four sites shows that the new scheme enables a more realistic representation of dry deposition. However, the representation is strongly limited by the local meteorology. In total, the changes increase the dry deposition velocity of ozone up to a factor of 2 globally, whereby the highest impact arises from the inclusion of cuticular uptake, especially over moist surfaces. This corresponds to a 6 % increase of global annual dry deposition loss of ozone resulting globally in a slight decrease of ground-level ozone but a regional decrease of up to 25 %. The change of ozone dry deposition is also reasoned by the altered loss of ozone precursors. Thus, the revision of the process parameterization as documented here has, among others, the potential to significantly reduce the overestimation of tropospheric ozone in global models.

Weaker cooling by aerosols due to dust–pollution interactions

Abstract. The interactions between aeolian dust and anthropogenic air pollution, notably chemical ageing of mineral dust and coagulation of dust and pollution particles, modify the atmospheric aerosol composition and burden. Since the aerosol particles can act as cloud condensation nuclei, this affects the radiative transfer not only directly via aerosol–radiation interactions, but also indirectly through cloud adjustments. We study both radiative effects using the global ECHAM/MESSy atmospheric chemistry-climate model (EMAC) which combines the Modular Earth Submodel System (MESSy) with the European Centre/Hamburg (ECHAM) climate model. Our simulations show that dust–pollution–cloud interactions reduce the condensed water path and hence the reflection of solar radiation. The associated climate warming outweighs the cooling that the dust–pollution interactions exert through the direct radiative effect. In total, this results in a net warming by dust–pollution interactions which moderates the negative global anthropogenic aerosol forcing at the top of the atmosphere by (0.2 ± 0.1) W m−2.

Newly developed aircraft routing options for air traffic simulation in the chemistry–climate model EMAC 2.53: AirTraf 2.0

Abstract. Aviation contributes to climate change, and the climate impact of aviation is expected to increase further. Adaptations of aircraft routings in order to reduce the climate impact are an important climate change mitigation measure. The air traffic simulator AirTraf, as a submodel of the European Center HAMburg general circulation model (ECHAM) and Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) model, enables the evaluation of such measures. For the first version of the submodel AirTraf, we concentrated on the general setup of the model, including departure and arrival, performance and emissions, and technical aspects such as the parallelization of the aircraft trajectory calculation with only a limited set of optimization possibilities (time and distance). Here, in the second version of AirTraf, we focus on enlarging the objective functions by seven new options to enable assessing operational improvements in many more aspects including economic costs, contrail occurrence, and climate impact. We verify that the AirTraf setup, e.g., in terms of number and choice of design variables for the genetic algorithm, allows us to find solutions even with highly structured fields such as contrail occurrence. This is shown by example simulations of the new routing options, including around 100 North Atlantic flights of an Airbus A330 aircraft for a typical winter day. The results clearly show that AirTraf 2.0 can find the different families of optimum flight trajectories (three-dimensional) for specific routing options; those trajectories minimize the corresponding objective functions successfully. The minimum cost option lies between the minimum time and the minimum fuel options. Thus, aircraft operating costs are minimized by taking the best compromise between flight time and fuel use. The aircraft routings for contrail avoidance and minimum climate impact reduce the potential climate impact which is estimated by using algorithmic climate change functions, whereas these two routings increase the aircraft operating costs. A trade-off between the aircraft operating costs and the climate impact is confirmed. The simulation results are compared with literature data, and the consistency of the submodel AirTraf 2.0 is verified.

The Flexible Ocean and Climate Infrastructure version 1 (FOCI1): mean state and variability

Abstract. A new Earth system model, the Flexible Ocean and Climate Infrastructure (FOCI), is introduced. A first version of FOCI consists of a global high-top atmosphere (European Centre Hamburg general circulation model; ECHAM6.3) and an ocean model (Nucleus for European Modelling of the Ocean v3.6; NEMO3.6) as well as sea-ice (Louvain-la-Neuve sea Ice Model version 2; LIM2) and land surface model components (Jena Scheme for Biosphere Atmosphere Coupling in Hamburg; JSBACH), which are coupled through the OASIS3-MCT software package. FOCI includes a number of optional modules which can be activated depending on the scientific question of interest. In the atmosphere, interactive stratospheric chemistry can be used (ECHAM6-HAMMOZ) to study, for example, the effects of the ozone hole on the climate system. In the ocean, a biogeochemistry model (Model of Oceanic Pelagic Stoichiometry; MOPS) is available to study the global carbon cycle. A unique feature of FOCI is the ability to explicitly resolve mesoscale ocean eddies in specific regions. This is realized in the ocean through nesting; first examples for the Agulhas Current and the Gulf Stream systems are described here. FOCI therefore bridges the gap between coarse-resolution climate models and global high-resolution weather prediction and ocean-only models. It allows to study the evolution of the climate system on regional and seasonal to (multi)decadal scales. The development of FOCI resulted from a combination of the long-standing expertise in ocean and climate modeling in several research units and divisions at the Helmholtz Centre for Ocean Research Kiel (GEOMAR). FOCI will thus be used to complement and interpret long-term observations in the Atlantic, enhance the process understanding of the role of mesoscale oceanic eddies for large-scale oceanic and atmospheric circulation patterns, study feedback mechanisms with stratospheric processes, estimate future ocean acidification, and improve the simulation of the Atlantic Meridional Overturning Circulation changes and their influence on climate, ocean chemistry and biology. In this paper, we present both the scientific vision for the development of FOCI as well as some technical details. This includes a first validation of the different model components using several configurations of FOCI. Results show that the model in its basic configuration runs stably under pre-industrial control as well as under historical forcing and produces a mean climate and variability which compares well with observations, reanalysis products and other climate models. The nested configurations reduce some long-standing biases in climate models and are an important step forward to include the atmospheric response in multidecadal eddy-rich configurations.

Potential Influence of a Developing La Niña on the Sea-Ice Reduction in the Barents–Kara Seas

The potential influence of a developing La Niña on Arctic sea-ice annual variability is investigated using both observational data and an atmospheric general circulation model. It is found that during the developing phase of an eastern Pacific (EP) La Niña event in June, July, and August (JJA) and September, October, and November (SON), the sea-ice concentration (SIC) over the Barents–Kara Seas declines more than 15%. The local atmospheric circulation pattern associated with the EP La Niña is characterized as a weak decrease in geopotential height over the Barents–Kara Seas, combined with an anticyclone in the North Atlantic. The corresponding southerly winds push warm waters northward into the key sea-ice reduction region and directly accelerate sea-ice melt. Meanwhile, the abundant moisture contained in the lower troposphere is transported into the Arctic region by winds resulting from the local barotropic structure. The humid atmosphere contributes to both net shortwave and longwave radiation and thus indirectly accelerates the decline in sea ice. Simulations by the European Centre Hamburg Model, version 5.4, are forced by observed sea surface temperature anomalies associated with EP La Niña events. The results of the simulations capture the North Atlantic anticyclone and reproduce the moisture transport, which supports the premise that an EP La Niña plays a crucial role in sea-ice reduction over the Barents–Kara sector from the perspective of atmospheric circulation and net surface heat flux.

Water stable isotope spatio-temporal variability in Antarctica in 1960–2013: observations and simulations from the ECHAM5-wiso atmospheric general circulation model

Abstract. Polar ice core water isotope records are commonly used to infer past changes in Antarctic temperature, motivating an improved understanding and quantification of the temporal relationship between δ18O and temperature. This can be achieved using simulations performed by atmospheric general circulation models equipped with water stable isotopes. Here, we evaluate the skills of the high-resolution water-isotope-enabled atmospheric general circulation model ECHAM5-wiso (the European Centre Hamburg Model) nudged to European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis using simulations covering the period 1960–2013 over the Antarctic continent. We compare model outputs with field data, first with a focus on regional climate variables and second on water stable isotopes, using our updated dataset of water stable isotope measurements from precipitation, snow, and firn–ice core samples. ECHAM5-wiso simulates a large increase in temperature from 1978 to 1979, possibly caused by a discontinuity in the European Reanalyses (ERA) linked to the assimilation of remote sensing data starting in 1979. Although some model–data mismatches are observed, the (precipitation minus evaporation) outputs are found to be realistic products for surface mass balance. A warm model bias over central East Antarctica and a cold model bias over coastal regions explain first-order δ18O model biases by too-strong isotopic depletion on coastal areas and underestimated depletion inland. At the second order, despite these biases, ECHAM5-wiso correctly captures the observed spatial patterns of deuterium excess. The results of model–data comparisons for the inter-annual δ18O standard deviation differ when using precipitation or ice core data. Further studies should explore the importance of deposition and post-deposition processes affecting ice core signals and not resolved in the model. These results build trust in the use of ECHAM5-wiso outputs to investigate the spatial, seasonal, and inter-annual δ18O–temperature relationships. We thus make the first Antarctica-wide synthesis of prior results. First, we show that local spatial or seasonal slopes are not a correct surrogate for inter-annual temporal slopes, leading to the conclusion that the same isotope–temperature slope cannot be applied for the climatic interpretation of Antarctic ice core for all timescales. Finally, we explore the phasing between the seasonal cycles of deuterium excess and δ18O as a source of information on changes in moisture sources affecting the δ18O–temperature relationship. The few available records and ECHAM5-wiso show different phase relationships in coastal, intermediate, and central regions. This work evaluates the use of the ECHAM5-wiso model as a tool for the investigation of water stable isotopes in Antarctic precipitation and calls for extended studies to improve our understanding of such proxies.

The on-line coupled atmospheric chemistry model system MECO(n) – Part 5: Expanding the Multi-Model-Driver (MMD v2.0) for 2-way data exchange including data interpolation via GRID (v1.0)

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.

Ozone impacts of gas–aerosol uptake in global chemistry transport models

Abstract. The impact of six heterogeneous gas–aerosol uptake reactions on tropospheric ozone and nitrogen species was studied using two chemical transport models, the Meteorological Synthesizing Centre-West of the European Monitoring and Evaluation Programme (EMEP MSC-W) and the European Centre Hamburg general circulation model combined with versions of the Hamburg Aerosol Model and Model for Ozone and Related chemical Tracers (ECHAM-HAMMOZ). Species undergoing heterogeneous reactions in both models include N2O5, NO3, NO2, O3, HNO3, and HO2. Since heterogeneous reactions take place at the aerosol surface area, the modelled surface area density (Sa) of both models was compared to a satellite product retrieving the surface area. This comparison shows a good agreement in global pattern and especially the capability of both models to capture the extreme aerosol loadings in east Asia. The impact of the heterogeneous reactions was evaluated by the simulation of a reference run containing all heterogeneous reactions and several sensitivity runs. One reaction was turned off in each sensitivity run to compare it with the reference run. The analysis of the sensitivity runs confirms that the globally most important heterogeneous reaction is the one of N2O5. Nevertheless, NO2, HNO3, and HO2 heterogeneous reactions gain relevance particularly in east Asia due to the presence of high NOx concentrations and high Sa in the same region. The heterogeneous reaction of O3 itself on dust is of minor relevance compared to the other heterogeneous reactions. The impacts of the N2O5 reactions show strong seasonal variations, with the biggest impacts on O3 in springtime when photochemical reactions are active and N2O5 levels still high. Evaluation of the models with northern hemispheric ozone surface observations yields a better agreement of the models with observations in terms of concentration levels, variability, and temporal correlations at most sites when the heterogeneous reactions are incorporated. Our results are loosely consistent with results from earlier studies, although the magnitude of changes induced by N2O5 reaction is at the low end of estimates, which seems to fit a trend, whereby the more recent the study the lower the impacts of these reactions.

Lagrangian process attribution of isotopic variations in near-surface water vapour in a 30-year regional climate simulation over Europe

Abstract. Stable water isotopes are naturally available tracers of moisture in the atmosphere. Due to isotopic fractionation, they record information about condensation and evaporation processes during the transport of air parcels, and therefore present a valuable means for studying the global water cycle. However, the meteorological processes driving isotopic variations are complex and not very well understood so far, in particular on short (hourly to daily) timescales. This study presents a Lagrangian method for attributing the isotopic composition of air parcels to meteorological processes, which provides new insight into the isotopic history of air parcels. It is based on the temporal evolution of the isotope ratios, the humidity, the temperature, and the location of the air parcels. Here these values are extracted along 7-day backward trajectories started every 6 hours from near the surface in a 30-year regional climate simulation over Europe with the isotope-enabled version of the model of the Consortium for Small-Scale Modelling (COSMOiso). The COSMOiso simulation has a horizontal resolution of 0.25∘ and is driven at the lateral boundaries by a T106 global climate simulation with the isotope-enabled version of the European Centre Hamburg model (ECHAMwiso). Both simulations are validated against measurements from the Global Network of Isotopes in Precipitation (GNIP), which shows that nesting COSMOiso within ECHAMwiso improves the representation of δ2H and deuterium excess in monthly accumulated precipitation. The method considers all isotopic changes that occur inside the COSMOiso model domain, which, on average, correspond to more than half of the mean and variability in both δ2H and deuterium excess at the air parcels' arrival points. Along every trajectory, the variations in the isotope values are quantitatively decomposed into eight process categories (evaporation from the ocean, evapotranspiration from land, mixing with moister air, mixing with drier air, liquid cloud formation, mixed phase cloud formation, ice cloud formation, and no process). The results show that for air parcels arriving over the ocean, evaporation from the ocean is the primary factor controlling δ2H and deuterium excess. Over land, evapotranspiration from land and mixing with moister air are similarly important. Liquid and mixed phase cloud formation contribute to the variability of δ2H and deuterium excess, especially over continental Europe. In summary, the presented method helps to better understand the linkage between the meteorological history of air parcels and their isotopic composition, and may support the interpretation of stable water isotope measurements in future.

Synoptic Control of Contrail Cirrus Life Cycles and Their Modification Due to Reduced Soot Number Emissions

AbstractThe atmospheric state, aircraft emissions, and engine properties determine formation and initial properties of contrails. The synoptic situation controls microphysical and dynamical processes and causes a wide variability of contrail cirrus life cycles. A reduction of soot particle number emissions, resulting, for example, from the use of alternative fuels, strongly impacts initial ice crystal numbers and microphysical process rates of contrail cirrus. We use the European Centre/Hamburg (ECHAM) climate model version 5 including a contrail cirrus modul, studying process rates, properties, and life cycles of contrail cirrus clusters within different synoptic situations. The impact of reduced soot number emissions is approximated by a reduction in the initial ice crystal number, exemplarily studied for 80%. Contrail cirrus microphysical and macrophysical properties can depend much more strongly on the synoptic situation than on the initial ice crystal number. They can attain a large cover, optical depth, and ice water content in long‐lived and large‐scale ice‐supersaturated areas, making them particularly climate‐relevant. In those synoptic situations, the accumulated ice crystal loss due to sedimentation is increased by around 15% and the volume of contrail cirrus, exceeding an optical depth of 0.02, and their short‐wave radiative impact are strongly decreased due to reduced soot emissions. These reductions are of little consequence in short‐lived and small‐scale ice‐supersaturated areas, where contrail cirrus stay optically very thin and attain a low cover. The synoptic situations in which long‐lived and climate‐relevant contrail cirrus clusters can be found over the eastern U.S. occur in around 25% of cases.

A 60-year ice-core record of regional climate from Adélie Land, coastal Antarctica

Abstract. A 22.4 m-long shallow firn core was extracted during the 2006/2007 field season from coastal Adélie Land. Annual layer counting based on subannual analyses of δ18O and major chemical components was combined with 5 reference years associated with nuclear tests and non-retreat of summer sea ice to build the initial ice-core chronology (1946–2006), stressing uncertain counting for 8 years. We focus here on the resulting δ18O and accumulation records. With an average value of 21.8 ± 6.9 cm w.e. yr−1, local accumulation shows multi-decadal variations peaking in the 1980s, but no long-term trend. Similar results are obtained for δ18O, also characterised by a remarkably low and variable amplitude of the seasonal cycle. The ice-core records are compared with regional records of temperature, stake area accumulation measurements and variations in sea-ice extent, and outputs from two models nudged to ERA (European Reanalysis) atmospheric reanalyses: the high-resolution atmospheric general circulation model (AGCM), including stable water isotopes ECHAM5-wiso (European Centre Hamburg model), and the regional atmospheric model Modèle Atmosphérique Régional (AR). A significant linear correlation is identified between decadal variations in δ18O and regional temperature. No significant relationship appears with regional sea-ice extent. A weak and significant correlation appears with Dumont d'Urville wind speed, increasing after 1979. The model-data comparison highlights the inadequacy of ECHAM5-wiso simulations prior to 1979, possibly due to the lack of data assimilation to constrain atmospheric reanalyses. Systematic biases are identified in the ECHAM5-wiso simulation, such as an overestimation of the mean accumulation rate and its interannual variability, a strong cold bias and an underestimation of the mean δ18O value and its interannual variability. As a result, relationships between simulated δ18O and temperature are weaker than observed. Such systematic precipitation and temperature biases are not displayed by MAR, suggesting that the model resolution plays a key role along the Antarctic ice sheet coastal topography. Interannual variations in ECHAM5-wiso temperature and precipitation accurately capture signals from meteorological data and stake observations and are used to refine the initial ice-core chronology within 2 years. After this adjustment, remarkable positive (negative) δ18O anomalies are identified in the ice-core record and the ECHAM5-wiso simulation in 1986 and 2002 (1998–1999), respectively. Despite uncertainties associated with post-deposition processes and signal-to-noise issues, in one single coastal ice-core record, we conclude that the S1C1 core can correctly capture major annual anomalies in δ18O as well as multi-decadal variations. These findings highlight the importance of improving the network of coastal high-resolution ice-core records, and stress the skills and limitations of atmospheric models for accumulation and δ18O in coastal Antarctic areas. This is particularly important for the overall East Antarctic ice sheet mass balance.

Evaluation of Total Precipitable Water from CRCM4 using the NVAP-MEaSUREs Dataset and ERA-Interim Reanalysis Data

The fourth-generation Canadian Regional Climate Model’s (CRCM4) precipitable water is evaluated and compared with observational data and ERA-Interim reanalysis data over five Canadian basins with simulations driven by ERA-Interim (two) and global climate models (two). Considering the 22 years of data available in the observations, we analyze precipitable water’s behaviour through its annual cycle, its daily distribution, and its annual daily maxima. For the simulations driven by reanalyses, differences in annual daily maximum values and their correlations with observations are examined. In general, the values for precipitable water simulated by CRCM4 are similar to those observed, and the model reproduces both the interannual and inter-basin variabilities. The simulation at 15 km resolution produces higher extreme values than simulations performed at 45 km resolution and higher than the observations taken at coarser resolution (1°), without much influence on the mean behaviour. Some underestimation is found with the simulation driven by the Canadian Centre for Climate Modelling and Analysis Model, version 3, a sign of a cold and dry bias, whereas the run driven by the European Centre Hamburg Model, version 5, is much closer to the observations, pointing to the importance of closely considering the regional–global model combination. Overall, CRCM4's ability to reproduce the major characteristics of observed precipitable water makes it a possible tool for providing precipitable water data that could serve as a basis for probable maximum precipitation and probable maximum flood studies at the basin scale.

Comparison of aerosol optical properties above clouds between POLDER and AeroCom models over the South East Atlantic Ocean during the fire season

AbstractAerosol properties above clouds have been retrieved over the South East Atlantic Ocean during the fire season 2006 using satellite observations from POLDER (Polarization and Directionality of Earth Reflectances). From June to October, POLDER has observed a mean Above‐Cloud Aerosol Optical Thickness (ACAOT) of 0.28 and a mean Above‐Clouds Single Scattering Albedo (ACSSA) of 0.87 at 550 nm. These results have been used to evaluate the simulation of aerosols above clouds in five Aerosol Comparisons between Observations and Models (Goddard Chemistry Aerosol Radiation and Transport (GOCART), Hadley Centre Global Environmental Model 3 (HadGEM3), European Centre Hamburg Model 5‐Hamburg Aerosol Module 2 (ECHAM5‐HAM2), Oslo‐Chemical Transport Model 2 (OsloCTM2), and Spectral Radiation‐Transport Model for Aerosol Species (SPRINTARS)). Most models do not reproduce the observed large aerosol load episodes. The comparison highlights the importance of the injection height and the vertical transport parameterizations to simulate the large ACAOT observed by POLDER. Furthermore, POLDER ACSSA is best reproduced by models with a high imaginary part of black carbon refractive index, in accordance with recent recommendations.

Response of Northern Hemispheric Air Temperature to Arctic Sea Ice Decline

Two numerical experiments were performed by using the Community Atmosphere Model version 3 (CAM3) with different sea ice datasets to assess the contribution of the decline of Arctic sea ice to warming in the Northern Hemisphere. One experiment was driven by observed sea ice cover data; for the other one, the authors used the sea ice data of the 4×CO2 scenario simulated by the fourth-generation European Centre Hamburg atmospheric general circulation Model of Istituto Nazionale di Geofisica e Vulcanologia, Italy (INGV_ECHAM4). The comparison of the two experiments indicates that the decline of the Arctic sea ice leads to a dramatic warming over the high latitudes of the Northern Hemisphere, characterized by a maximum warming of more than 26℃ over the Arctic region. The significant warming is closely related to the enhanced atmospheric heat source. A 40–60 Wm-2 increase in the apparent heat source was simulated in winter due to the decline of Arctic sea ice. In contrast, no significant change was found in the atmospheric apparent heat source in summer. As a result, the summer temperature change induced by the decline of Arctic sea ice appears to be weak. This study suggests that accurate sea ice cover data is crucial for future climate projection of air temperature in high latitudes.