Search Results

Now showing 1 - 10 of 10
  • Item
    The global aerosol-climate model ECHAM6.3-HAM2.3-Part 2: Cloud evaluation, aerosol radiative forcing, and climate sensitivity
    (Katlenburg-Lindau : Copernicus, 2019) Neubauer, David; Ferrachat, Sylvaine; Siegenthaler-Le Drian, Colombe; Stier, Philip; Partridge, Daniel G.; Tegen, Ina; Bey, Isabelle; Stanelle, Tanja; Kokkola, Harri; Lohmann, Ulrike
    The global aerosol–climate model ECHAM6.3–HAM2.3 (E63H23) as well as the previous model versions ECHAM5.5–HAM2.0 (E55H20) and ECHAM6.1–HAM2.2 (E61H22) are evaluated using global observational datasets for clouds and precipitation. In E63H23, the amount of low clouds, the liquid and ice water path, and cloud radiative effects are more realistic than in previous model versions. E63H23 has a more physically based aerosol activation scheme, improvements in the cloud cover scheme, changes in the detrainment of convective clouds, changes in the sticking efficiency for the accretion of ice crystals by snow, consistent ice crystal shapes throughout the model, and changes in mixed-phase freezing; an inconsistency in ice crystal number concentration (ICNC) in cirrus clouds was also removed. Common biases in ECHAM and in E63H23 (and in previous ECHAM–HAM versions) are a cloud amount in stratocumulus regions that is too low and deep convective clouds over the Atlantic and Pacific oceans that form too close to the continents (while tropical land precipitation is underestimated). There are indications that ICNCs are overestimated in E63H23. Since clouds are important for effective radiative forcing due to aerosol–radiation and aerosol–cloud interactions (ERFari+aci) and equilibrium climate sensitivity (ECS), differences in ERFari+aci and ECS between the model versions were also analyzed. ERFari+aci is weaker in E63H23 (−1.0 W m−2) than in E61H22 (−1.2 W m−2) (or E55H20; −1.1 W m−2). This is caused by the weaker shortwave ERFari+aci (a new aerosol activation scheme and sea salt emission parameterization in E63H23, more realistic simulation of cloud water) overcompensating for the weaker longwave ERFari+aci (removal of an inconsistency in ICNC in cirrus clouds in E61H22). The decrease in ECS in E63H23 (2.5 K) compared to E61H22 (2.8 K) is due to changes in the entrainment rate for shallow convection (affecting the cloud amount feedback) and a stronger cloud phase feedback. Experiments with minimum cloud droplet number concentrations (CDNCmin) of 40 cm−3 or 10 cm−3 show that a higher value of CDNCmin reduces ERFari+aci as well as ECS in E63H23.
  • Item
    Modelling mineral dust emissions and atmospheric dispersion with MADE3 in EMAC v2.54
    (Katlenburg-Lindau : Copernicus, 2020) Beer, Christof G.; Hendricks, Johannes; Righi, Mattia; Heinold, Bernd; Tegen, Ina; Groß, Silke; Sauer, Daniel; Walser, Adrian; Weinzierl, Bernadett
    It was hypothesized that using mineral dust emission climatologies in global chemistry climate models (GCCMs), i.e. prescribed monthly-mean dust emissions representative of a specific year, may lead to misrepresentations of strong dust burst events. This could result in a negative bias of model dust concentrations compared to observations for these episodes. Here, we apply the aerosol microphysics submodel MADE3 (Modal Aerosol Dynamics model for Europe, adapted for global applications, third generation) as part of the ECHAM/MESSy Atmospheric Chemistry (EMAC) general circulation model. We employ two different representations of mineral dust emissions for our model simulations: (i) a prescribed monthly-mean climatology of dust emissions representative of the year 2000 and (ii) an online dust parametrization which calculates wind-driven mineral dust emissions at every model time step. We evaluate model results for these two dust representations by comparison with observations of aerosol optical depth from ground-based station data. The model results show a better agreement with the observations for strong dust burst events when using the online dust representation compared to the prescribed dust emissions setup. Furthermore, we analyse the effect of increasing the vertical and horizontal model resolution on the mineral dust properties in our model. We compare results from simulations with T42L31 and T63L31 model resolution (2.8∘×2.8∘ and 1.9∘×1.9∘ in latitude and longitude, respectively; 31 vertical levels) with the reference setup (T42L19). The different model versions are evaluated against airborne in situ measurements performed during the SALTRACE mineral dust campaign (Saharan Aerosol Long-range Transport and Aerosol-Cloud Interaction Experiment, June–July 2013), i.e. observations of dust transported from the Sahara to the Caribbean. Results show that an increased horizontal and vertical model resolution is able to better represent the spatial distribution of airborne mineral dust, especially in the upper troposphere (above 400 hPa). Additionally, we analyse the effect of varying assumptions for the size distribution of emitted dust but find only a weak sensitivity concerning these changes. The results of this study will help to identify the model setup best suited for future studies and to further improve the representation of mineral dust particles in EMAC-MADE3.
  • Item
    EC-Earth3-AerChem: a global climate model with interactive aerosols and atmospheric chemistry participating in CMIP6
    (Katlenburg-Lindau : Copernicus, 2021-9-13) van Noije, Twan; Bergman, Tommi; Le Sager, Philippe; O'Donnell, Declan; Makkonen, Risto; Gonçalves-Ageitos, María; Döscher, Ralf; Fladrich, Uwe; von Hardenberg, Jost; Keskinen, Jukka-Pekka; Korhonen, Hannele; Laakso, Anton; Myriokefalitakis, Stelios; Ollinaho, Pirkka; Pérez García-Pando, Carlos; Reerink, Thomas; Schrödner, Roland; Wyser, Klaus; Yang, Shuting
    This paper documents the global climate model EC-Earth3-AerChem, one of the members of the EC-Earth3 family of models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6). EC-Earth3-AerChem has interactive aerosols and atmospheric chemistry and contributes to the Aerosols and Chemistry Model Intercomparison Project (AerChemMIP). In this paper, we give an overview of the model, describe in detail how it differs from the other EC-Earth3 configurations, and outline the new features compared with the previously documented version of the model (EC-Earth 2.4). We explain how the model was tuned and spun up under preindustrial conditions and characterize the model's general performance on the basis of a selection of coupled simulations conducted for CMIP6. The net energy imbalance at the top of the atmosphere in the preindustrial control simulation is on average −0.09 W m−2 with a standard deviation due to interannual variability of 0.25 W m−2, showing no significant drift. The global surface air temperature in the simulation is on average 14.08 ∘C with an interannual standard deviation of 0.17 ∘C, exhibiting a small drift of 0.015 ± 0.005 ∘C per century. The model's effective equilibrium climate sensitivity is estimated at 3.9 ∘C, and its transient climate response is estimated at 2.1 ∘C. The CMIP6 historical simulation displays spurious interdecadal variability in Northern Hemisphere temperatures, resulting in a large spread across ensemble members and a tendency to underestimate observed annual surface temperature anomalies from the early 20th century onwards. The observed warming of the Southern Hemisphere is well reproduced by the model. Compared with the ECMWF (European Centre for Medium-Range Weather Forecasts) Reanalysis version 5 (ERA5), the surface air temperature climatology for 1995–2014 has an average bias of −0.86 ± 0.05 ∘C with a standard deviation across ensemble members of 0.35 ∘C in the Northern Hemisphere and 1.29 ± 0.02 ∘C with a corresponding standard deviation of 0.05 ∘C in the Southern Hemisphere. The Southern Hemisphere warm bias is largely caused by errors in shortwave cloud radiative effects over the Southern Ocean, a deficiency of many climate models. Changes in the emissions of near-term climate forcers (NTCFs) have significant effects on the global climate from the second half of the 20th century onwards. For the SSP3-7.0 Shared Socioeconomic Pathway, the model gives a global warming at the end of the 21st century (2091–2100) of 4.9 ∘C above the preindustrial mean. A 0.5 ∘C stronger warming is obtained for the AerChemMIP scenario with reduced emissions of NTCFs. With concurrent reductions of future methane concentrations, the warming is projected to be reduced by 0.5 ∘C.
  • Item
    CAPRAM reduction towards an operational multiphase halogen and dimethyl sulfide chemistry treatment in the chemistry transport model COSMO-Muscat(5.04e)
    (Katlenburg-Lindau : Copernicus, 2020) Hoffmann, Erik H.; Schrödner, Roland; Tilgner, Andreas; Wolke, Ralf; Herrmann, Hartmut
    A condensed multiphase halogen and dimethyl sulfide (DMS) chemistry mechanism for application in chemistry transport models is developed by reducing the CAPRAM DMS module 1.0 (CAPRAM-DM1.0) and the CAPRAM halogen module 3.0 (CAPRAM-HM3.0). The reduction is achieved by determining the main oxidation pathways from analysing the mass fluxes of complex multiphase chemistry simulations with the air parcel model SPACCIM (SPectral Aerosol Cloud Chemistry Interaction Model). These simulations are designed to cover both pristine and polluted marine boundary layer conditions. Overall, the reduced CAPRAM-DM1.0 contains 32 gas-phase reactions, 5 phase transfers, and 12 aqueous-phase reactions, of which two processes are described as equilibrium reactions. The reduced CAPRAM-HM3.0 contains 199 gas-phase reactions, 23 phase transfers, and 87 aqueous-phase reactions. For the aqueous-phase chemistry, 39 processes are described as chemical equilibrium reactions. A comparison of simulations using the complete CAPRAM-DM1.0 and CAPRAM-HM3.0 mechanisms against the reduced ones indicates that the relative deviations are below 5 % for important inorganic and organic air pollutants and key reactive species under pristine ocean and polluted conditions. The reduced mechanism has been implemented into the chemical transport model COSMO-MUSCAT and tested by performing 2D simulations under prescribed meteorological conditions that investigate the effect of stable (stratiform cloud) and more unstable meteorological conditions (convective clouds) on marine multiphase chemistry. The simulated maximum concentration of HCl is of the order of 109 molecules cm−3 and that of BrO is around 1×107 molecules cm−3, reproducing the range of ambient measurements. Afterwards, the oxidation pathways of DMS in a cloudy marine atmosphere have been investigated in detail. The simulations demonstrate that clouds have both a direct and an indirect photochemical effect on the multiphase processing of DMS and its oxidation products. The direct photochemical effect is related to in-cloud chemistry that leads to high dimethyl sulfoxide (DMSO) oxidation rates and a subsequently enhanced formation of methane sulfonic acid compared to aerosol chemistry. The indirect photochemical effect is characterized by cloud shading, which occurs particularly in the case of stratiform clouds. The lower photolysis rate affects the activation of Br atoms and consequently lowers the formation of BrO radicals. The corresponding DMS oxidation flux is lowered by up to 30 % under thick optical clouds. Moreover, high updraught velocities lead to a strong vertical mixing of DMS into the free troposphere predominately under cloudy conditions. The photolysis of hypohalous acids (HOX, X = Cl, Br, or I) is reduced as well, resulting in higher HOX-driven sulfite-to-sulfate oxidation in aerosol particles below stratiform clouds. Altogether, the present model simulations have demonstrated the ability of the reduced mechanism to be applied in studying marine aerosol–cloud processing effects in regional models such as COSMO-MUSCAT. The reduced mechanism can be used also by other regional models for more adequate interpretations of complex marine field measurement data.
  • Item
    A new Lagrangian in-time particle simulation module (Itpas v1) for atmospheric particle dispersion
    (Katlenburg-Lindau : Copernicus, 2021) Faust, Matthias; Wolke, Ralf; Münch, Steffen; Funk, Roger; Schepanski, Kerstin
    Trajectory models are intuitive tools for airflow studies. But in general, they are limited to non-turbulent, i.e. laminar flow, conditions. Therefore, trajectory models are not particularly suitable for investigating airflow within the turbulent atmospheric boundary layer. To overcome this, a common approach is handling the turbulent uncertainty as a random deviation from a mean path in order to create a statistic of possible solutions which envelops the mean path. This is well known as the Lagrangian particle dispersion model (LPDM). However, the decisive factor is the representation of turbulence in the model, for which widely used models such as FLEXPART and HYSPLIT use an approximation. A conceivable improvement could be the use of a turbulence parameterisation approach based on the turbulent kinetic energy (TKE) at high temporal resolution. Here, we elaborated this approach and developed the LPDM Itpas, which is coupled online to the German Weather Service's mesoscale weather forecast model COSMO. It benefits from the prognostically calculated TKE as well as from the high-frequency wind information. We demonstrate the model's applicability for a case study on agricultural particle emission in eastern Germany. The results obtained are discussed with regard to the model's ability to describe particle transport within a turbulent boundary layer. Ultimately, the simulations performed suggest that the newly introduced method based on prognostic TKE sufficiently represents the particle transport.
  • Item
    Implementation of aerosol-cloud interactions in the regional atmosphere-aerosol model COSMO-Muscat(5.0) and evaluation using satellite data
    (Katlenburg-Lindau : Copernicus, 2017) Dipu, Sudhakar; Quaas, Johannes; Wolke, Ralf; Stoll, Jens; Mühlbauer, Andreas; Sourdeval, Odran; Salzmann, Marc; Heinold, Bernd; Tegen, Ina
    The regional atmospheric model Consortium for Small-scale Modeling (COSMO) coupled to the Multi-Scale Chemistry Aerosol Transport model (Muscat) is extended in this work to represent aerosol-cloud interactions. Previously, only one-way interactions (scavenging of aerosol and in-cloud chemistry) and aerosol-radiation interactions were included in this model. The new version allows for a microphysical aerosol effect on clouds. For this, we use the optional two-moment cloud microphysical scheme in COSMO and the online-computed aerosol information for cloud condensation nuclei concentrations (Cccn), replacing the constant Cccn profile. In the radiation scheme, we have implemented a droplet-size-dependent cloud optical depth, allowing now for aerosol-cloud-radiation interactions. To evaluate the models with satellite data, the Cloud Feedback Model Intercomparison Project Observation Simulator Package (COSP) has been implemented. A case study has been carried out to understand the effects of the modifications, where the modified modeling system is applied over the European domain with a horizontal resolution of 0.25°g × g0.25°. To reduce the complexity in aerosol-cloud interactions, only warm-phase clouds are considered. We found that the online-coupled aerosol introduces significant changes for some cloud microphysical properties. The cloud effective radius shows an increase of 9.5g%, and the cloud droplet number concentration is reduced by 21.5g%.
  • Item
    On the application and grid-size sensitivity of the urban dispersion model CAIRDIO v2.0 under real city weather conditions
    (Katlenburg-Lindau : Copernicus, 2022) Weger, Michael; Baars, Holger; Gebauer, Henriette; Merkel, Maik; Wiedensohler, Alfred; Heinold, Bernd
    There is a gap between the need for city-wide air-quality simulations considering the intra-urban variability and mircoscale dispersion features and the computational capacities that conventional urban microscale models require. This gap can be bridged by targeting model applications on the gray zone situated between the mesoscale and large-eddy scale. The urban dispersion model CAIRDIO is a new contribution to the class of computational-fluid dynamics models operating in this scale range. It uses a diffuse-obstacle boundary method to represent buildings as physical obstacles at gray-zone resolutions in the order of tens of meters. The main objective of this approach is to find an acceptable compromise between computationally inexpensive grid sizes for spatially comprehensive applications and the required accuracy in the description of building and boundary-layer effects. In this paper, CAIRDIO is applied on the simulation of black carbon and particulate matter dispersion for an entire mid-size city using a uniform horizontal grid spacing of 40gm. For model evaluation, measurements from five operational air monitoring stations representative for the urban background and high-traffic roads are used. The comparison also includes the mesoscale host simulation, which provides the boundary conditions. The measurements show a dominant influence of the mixing layer evolution at background sites, and therefore both the mesoscale and large-eddy simulation (LES) results are in good agreement with the observed air pollution levels. In contrast, at the high-traffic sites the proximity to emissions and the interactions with the building environment lead to a significantly amplified diurnal variability in pollutant concentrations. These urban road conditions can only be reasonably well represented by CAIRDIO while the meosocale simulation indiscriminately reproduces a typical urban-background profile, resulting in a large positive model bias. Remaining model discrepancies are further addressed by a grid-spacing sensitivity study using offline-nested refined domains. The results show that modeled peak concentrations within street canyons can be further improved by decreasing the horizontal grid spacing down to 10gm, but not beyond. Obviously, the default grid spacing of 40gm is too coarse to represent the specific environment within narrow street canyons. The accuracy gains from the grid refinements are still only modest compared to the remaining model error, which to a large extent can be attributed to uncertainties in the emissions. Finally, the study shows that the proposed gray-scale modeling is a promising downscaling approach for urban air-quality applications. The results, however, also show that aspects other than the actual resolution of flow patterns and numerical effects can determine the simulations at the urban microscale.
  • Item
    Coupling aerosols to (cirrus) clouds in the global EMAC-MADE3 aerosol–climate model
    (Katlenburg-Lindau : Copernicus, 2020) Righi, Mattia; Hendricks, Johannes; Lohmann, Ulrike; Beer, Christof Gerhard; Hahn, Valerian; Heinold, Bernd; Heller, Romy; Krämer, Martina; Ponater, Michael; Rolf, Christian; Tegen, Ina; Voigt, Christiane
    A new cloud microphysical scheme including a detailed parameterization for aerosol-driven ice formation in cirrus clouds is implemented in the global ECHAM/MESSy Atmospheric Chemistry (EMAC) chemistry–climate model and coupled to the third generation of the Modal Aerosol Dynamics model for Europe adapted for global applications (MADE3) aerosol submodel. The new scheme is able to consistently simulate three regimes of stratiform clouds – liquid, mixed-, and ice-phase (cirrus) clouds – considering the activation of aerosol particles to form cloud droplets and the nucleation of ice crystals. In the cirrus regime, it allows for the competition between homogeneous and heterogeneous freezing for the available supersaturated water vapor, taking into account different types of ice-nucleating particles, whose specific ice-nucleating properties can be flexibly varied in the model setup. The new model configuration is tuned to find the optimal set of parameters that minimizes the model deviations with respect to observations. A detailed evaluation is also performed comparing the model results for standard cloud and radiation variables with a comprehensive set of observations from satellite retrievals and in situ measurements. The performance of EMAC-MADE3 in this new coupled configuration is in line with similar global coupled models and with other global aerosol models featuring ice cloud parameterizations. Some remaining discrepancies, namely a high positive bias in liquid water path in the Northern Hemisphere and overestimated (underestimated) cloud droplet number concentrations over the tropical oceans (in the extratropical regions), which are both a common problem in these kinds of models, need to be taken into account in future applications of the model. To further demonstrate the readiness of the new model system for application studies, an estimate of the anthropogenic aerosol effective radiative forcing (ERF) is provided, showing that EMAC-MADE3 simulates a relatively strong aerosol-induced cooling but within the range reported in the Intergovernmental Panel on Climate Change (IPCC) assessments.
  • Item
    SALSA2.0: The sectional aerosol module of the aerosol-chemistry-climate model ECHAM6.3.0-HAM2.3-MOZ1.0
    (Katlenburg-Lindau : Copernicus, 2018) Kokkola, Harri; Kühn, Thomas; Laakso, Anton; Bergman, Tommi; Lehtinen, Kari E. J.; Mielonen, Tero; Arola, Antti; Stadtler, Scarlet; Korhonen, Hannele; Ferrachat, Sylvaine; Lohmann, Ulrike; Neubauer, David; Tegen, Ina; Siegenthaler-Le Drian, Colombe; Schultz, Martin G.; Bey, Isabelle; Stier, Philip; Daskalakis, Nikos; Heald, Colette L.; Romakkaniemi, Sami
    In this paper, we present the implementation and evaluation of the aerosol microphysics module SALSA2.0 in the framework of the aerosol-chemistry-climate model ECHAM-HAMMOZ. It is an alternative microphysics module to the default modal microphysics scheme M7 in ECHAM-HAMMOZ. The SALSA2.0 implementation within ECHAM-HAMMOZ is evaluated against observations of aerosol optical properties, aerosol mass, and size distributions, comparing also to the skill of the M7 implementation. The largest differences between the implementation of SALSA2.0 and M7 are in the methods used for calculating microphysical processes, i.e., nucleation, condensation, coagulation, and hydration. These differences in the microphysics are reflected in the results so that the largest differences between SALSA2.0 and M7 are evident over regions where the aerosol size distribution is heavily modified by the microphysical processing of aerosol particles. Such regions are, for example, highly polluted regions and regions strongly affected by biomass burning. In addition, in a simulation of the 1991 Mt. Pinatubo eruption in which a stratospheric sulfate plume was formed, the global burden and the effective radii of the stratospheric aerosol are very different in SALSA2.0 and M7. While SALSA2.0 was able to reproduce the observed time evolution of the global burden of sulfate and the effective radii of stratospheric aerosol, M7 strongly overestimates the removal of coarse stratospheric particles and thus underestimates the effective radius of stratospheric aerosol. As the mode widths of M7 have been optimized for the troposphere and were not designed to represent stratospheric aerosol, the ability of M7 to simulate the volcano plume was improved by modifying the mode widths, decreasing the standard deviations of the accumulation and coarse modes from 1.59 and 2.0, respectively, to 1.2 similar to what was observed after the Mt. Pinatubo eruption. Overall, SALSA2.0 shows promise in improving the aerosol description of ECHAM-HAMMOZ and can be further improved by implementing methods for aerosol processes that are more suitable for the sectional method, e.g., size-dependent emissions for aerosol species and size-resolved wet deposition.
  • Item
    ICONGETM v1.0 – flexible NUOPC-driven two-way coupling via ESMF exchange grids between the unstructured-grid atmosphere model ICON and the structured-grid coastal ocean model GETM
    (Katlenburg-Lindau : Copernicus, 2021) Bauer, Tobias Peter; Holtermann, Peter; Heinold, Bernd; Radtke, Hagen; Knoth, Oswald; Klingbeil, Knut
    Two-way model coupling is important for representing the mutual interactions and feedbacks between atmosphere and ocean dynamics. This work presents the development of the two-way coupled model system ICONGETM, consisting of the atmosphere model ICON and the ocean model GETM. ICONGETM is built on the latest NUOPC coupling software with flexible data exchange and conservative interpolation via ESMF exchange grids. With ICON providing a state-of-the-art kernel for numerical weather prediction on an unstructured mesh and GETM being an established coastal ocean model, ICONGETM is especially suited for high-resolution studies. For demonstration purposes the newly developed model system has been applied to a coastal upwelling scenario in the central Baltic Sea.