The Modular Earth Submodel System (MESSy) is an open, multi-institutional project providing a strategy for developing comprehensive Earth System Models (ESMs) with highly flexible complexity. The first version of the MESSy infrastructure and process submodels, mainly focusing on atmospheric chemistry, has been successfully coupled to an atmospheric General Circulation Model (GCM) expanding it into an Atmospheric Chemistry GCM (AC-GCM) for nudged simulations and into a Chemistry Climate Model (CCM) for climate simulations. Here, we present the second development cycle of MESSy, which comprises (1) an improved and extended infrastructure for the basemodel independent coupling of process-submodels, (2) new, highly valuable diagnostic capabilities for the evaluation with observational data and (3) an improved atmospheric chemistry setup. With the infrastructural changes, we place the headstone for further model extensions from a CCM towards a comprehensive ESM. The new diagnostic submodels will be used for regular re-evaluations of the continuously further developing model system. The updates of the chemistry setup are briefly evaluated.

Three types of reference simulations, as recommended by the Chemistry–Climate Model Initiative (CCMI), have been performed with version 2.51 of the European Centre for Medium-Range Weather Forecasts – Hamburg (ECHAM)/Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) model: hindcast simulations (1950–2011), hindcast simulations with specified dynamics (1979–2013), i.e. nudged towards ERA-Interim reanalysis data, and combined hindcast and projection simulations (1950–2100). The manuscript summarizes the updates of the model system and details the different model set-ups used, including the on-line calculated diagnostics. Simulations have been performed with two different nudging set-ups, with and without interactive tropospheric aerosol, and with and without a coupled ocean model. Two different vertical resolutions have been applied. The on-line calculated sources and sinks of reactive species are quantified and a first evaluation of the simulation results from a global perspective is provided as a quality check of the data. The focus is on the intercomparison of the different model set-ups. The simulation data will become publicly available via CCMI and the Climate and Environmental Retrieval and Archive (CERA) database of the German Climate Computing Centre (DKRZ). This manuscript is intended to serve as an extensive reference for further analyses of the Earth System Chemistry integrated Modelling (ESCiMo) simulations.

Power plants and large industrial facilities contribute more than half of global anthropogenic CO₂ emissions. Quantifying the emissions of these point sources is therefore one of the main goals of the planned constellation of anthropogenic CO₂ monitoring satellites (CO2M) of the European Copernicus program. Atmospheric transport models may be used to study the capabilities of such satellites through observing system simulation experiments and to quantify emissions in an inverse modeling framework. How realistically the CO₂ plumes of power plants can be simulated and how strongly the results may depend on model type and resolution, however, is not well known due to a lack of observations available for benchmarking. Here, we use the unique data set of aircraft in situ and remote sensing observations collected during the CoMet (Carbon Dioxide and Methane Mission) measurement campaign downwind of the coal-fired power plants at Bełchatów in Poland and Jänschwalde in Germany in 2018 to evaluate the simulations of six different atmospheric transport models. The models include three large-eddy simulation (LES) models, two mesoscale numerical weather prediction (NWP) models extended for atmospheric tracer transport, and one Lagrangian particle dispersion model (LPDM) and cover a wide range of model resolutions from 200 m to 2 km horizontal grid spacing. At the time of the aircraft measurements between late morning and early afternoon, the simulated plumes were slightly (at Jänschwalde) to highly (at Bełchatów) turbulent, consistent with the observations, and extended over the whole depth of the atmospheric boundary layer (ABL; up to 1800 m a.s.l. (above sea level) in the case of Bełchatów). The stochastic nature of turbulent plumes puts fundamental limitations on a point-by-point comparison between simulations and observations. Therefore, the evaluation focused on statistical properties such as plume amplitude and width as a function of distance from the source. LES and NWP models showed similar performance and sometimes remarkable agreement with the observations when operated at a comparable resolution. The Lagrangian model, which was the only model driven by winds observed from the aircraft, quite accurately captured the location of the plumes but generally underestimated their width. A resolution of 1 km or better appears to be necessary to realistically capture turbulent plume structures. At a coarser resolution, the plumes disperse too quickly, especially in the near-field range (0–8 km from the source), and turbulent structures are increasingly smoothed out. Total vertical columns are easier to simulate accurately than the vertical distribution of CO₂, since the latter is critically affected by profiles of vertical stability, especially near the top of the ABL. Cross-sectional flux and integrated mass enhancement methods applied to synthetic CO2M data generated from the model simulations with a random noise of 0.5–1.0 ppm (parts per million) suggest that emissions from a power plant like Bełchatów can be estimated with an accuracy of about 20 % from single overpasses. Estimates of the effective wind speed are a critical input for these methods. Wind speeds in the middle of the ABL appear to be a good approximation for plumes in a well-mixed ABL, as encountered during CoMet.

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 present study investigates the initiation of precipitating deep convection in an ensemble of convection‐resolving mesoscale models. Results of eight different model runs from five non‐hydrostatic models are compared for a case of the Convective and Orographically‐induced Precipitation Study (COPS). An isolated convective cell initiated east of the Black Forest crest in southwest Germany, although convective available potential energy was only moderate and convective inhibition was high. Measurements revealed that, due to the absence of synoptic forcing, convection was initiated by local processes related to the orography. In particular, the lifting by low‐level convergence in the planetary boundary layer is assumed to be the dominant process on that day. The models used different configurations as well as different initial and boundary conditions. By comparing the different model performance with each other and with measurements, the processes which need to be well represented to initiate convection at the right place and time are discussed. Besides an accurate specification of the thermodynamic and kinematic fields, the results highlight the role of boundary‐layer convergence features for quantitative precipitation forecasts in mountainous terrain. Copyright © 2011 Royal Meteorological Society

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.