The Asian summer monsoon (ASM) traps convectively lifted boundary layer pollutants inside its upper-tropospheric lower-stratospheric Asian monsoon anticyclone (AMA). It is associated with a seasonal and spatially confined enhanced aerosol layer, called the Asian Tropopause Aerosol Layer (ATAL). Due to the dynamical variability of the AMA, the dearth of in situ observations in this region, the complexity of the emission sources and of transport pathways, knowledge of the ATAL properties in terms of aerosol budget, chemical composition, as well as its variability and temporal trend is still largely uncertain. In this work, we use the Community Earth System Model (CESM 1.2 version) based on the coupling of the Community Atmosphere Model (CAM5) and the MAM7 (Modal Aerosol Model) aerosol module to simulate the composition of the ATAL and its decadal trends. Our simulations cover a long-term period of 16 years from 2000 to 2015. We identify a typical “double-peak” vertical profile of aerosols for the ATAL. We attribute the upper peak (around 100 hPa, predominant during early ATAL, e.g., in June) to dry aerosols, possibly from nucleation processes, and the lower peak (around 250 hPa, predominant for a well-developed and late ATAL, e.g., in July and August) to cloud-borne aerosols associated with convective clouds. We find that mineral dust (present in both peaks) is the dominant aerosol by mass in the ATAL, showing a large interannual variability but no long-term trend, due to its natural variability. The results between 120 and 80 hPa (dry aerosol peak) suggest that for aerosols other than dust the ATAL is composed of around 40 % of sulfate, 30 % of secondary and 15 % of primary organic aerosols, 14 % of ammonium aerosols and less than 3 % of black carbon. Nitrate aerosols are not considered in MAM7. The analysis of the anthropogenic and biomass burning aerosols shows a positive trend for all aerosols simulated by CESM-MAM7.
During boreal summer, major convective activity is driven by the Asian
summer monsoon (ASM). The ASM-related convection combines both land
convection over mainland Asia and maritime convection over surrounding seas.
This dynamical mechanism acts as a pathway for the transport of trace gases
and pollutants from the boundary layer to the UTLS (Upper Troposphere Lower
Stratosphere) (Randel and Park, 2006; Park et al., 2007; Pan et al., 2016;
Gottschaldt et al., 2017). The upper atmospheric circulation is dominated by
the related Asian monsoon anticyclone (AMA), which is known to contain enhanced concentration of tropospheric trace gases and aerosols (Randel and
Park, 2006; Park et al., 2007, 2008), due to rapid lifting from
the boundary layer by deep convection and subsequent horizontal confinement.
The AMA is confined by the subtropical westerly jet stream in the north
(
The sources, chemical composition and spatial and temporal variability of the ATAL are not yet well understood. Recent observations from the StratoClim (Stratospheric and upper tropospheric processes for better climate predictions) aircraft campaign in 2017 and a few recent balloon measurements from the BATAL (Balloon measurement campaigns of the Asian Tropopause Aerosol Layer) 2015 campaign suggest that aerosol particles in the ATAL may contain large amounts of sulfate as well as organics, nitrates (including ammonium nitrate), black carbon and dust (Vernier et al., 2015, 2018; Höpfner et al., 2019). Different indications on the ATAL composition have been given by a number of modeling studies. Fadnavis et al. (2013), using the aerosol–chemistry–climate model ECHAM5-HAMMOZ, studied the transport of aerosols to the UTLS and showed persistent maxima in black carbon, organic carbon, sulfate, and mineral dust aerosols within the anticyclone throughout the ASM (from July to September). Yu et al. (2015), using the CESM1 (Community Earth System Model) global Earth system model coupled with the CARMA (Community Aerosol and Radiation Model for Atmospheres) aerosol model, suggested that the ATAL might be principally composed of secondary organic and sulfate aerosols as well as of primary organic aerosols. Fadnavis et al. (2017) performed model simulations with the ECHAM6-HAM (European Centre Hamburg Model 6.3-Hamburg Aerosol Model) global aerosol–climate model, and their simulations showed a persistent maximum of carbonaceous aerosols in the ATAL region. Ma et al. (2019), using the ECHAM/MESSy (Modular Earth Submodel System) for the Atmospheric Chemistry (EMAC) general circulation model coupled with the Global Modal-aerosol eXtension (GMXe) aerosol module, found that mineral dust and water-soluble compounds, like nitrate and sulfate, are the principal typology of aerosols over the Tibetan Plateau, within the AMA. Using the GEOS-Chem (Goddard Earth Observing System with Chemistry) chemical transport model, Fairlie et al. (2020) found significant amounts of sulfate, ammonium, organic aerosols and nitrate in the ATAL, with a predominant contribution of nitrate, as was identified previously by Gu et al. (2016) using an earlier version of the model. Therefore, existing modeling studies have proven able to simulate the enhanced concentration of aerosols in the AMA region, even if a very large uncertainty in the composition of the ATAL remains.
In several studies, dust has been shown to be a major contributor to the aerosol burden in the Asian upper troposphere during summer. Xu et al. (2015), using CALIOP and MISR (Multi-angle Imaging SpectroRadiometer) satellite data, found that dust is one of the predominant aerosols over the Tibetan Plateau, most probably originating from the Taklamakan desert and lofted from the surface to an altitude of about 10 km. Ma et al. (2019) simulated a broad maximum of dust surface concentration at the northern edge of the Tibetan Plateau up to 10 km. Their model results showed that the enhancement of dust aerosols is still visible up to 16 km above the Tibetan Plateau, with a maximum shifted to the east and south as a consequence of the influence of anticyclonic circulation. Large amounts of dust have also been reported by Lau et al. (2018) in the mid and upper troposphere over India and China from May to June transported from the Middle East desert and then from July to August trapped and accumulated within the AMA and contributing to the ATAL formation.
A rising temporal trend of the ATAL optical signature in the AMA region has been observed (Vernier et al., 2015). The recent rising trends of sulfur dioxide and volatile organic compounds emissions in India have been proposed as a candidate for explaining the appearance of the ATAL and its evolution. Continental convective regions have also been shown to be the main contributors to the air trapped within the AMA with northern India and south of the Tibetan Plateau as specific source areas (e.g., Tissier and Legras, 2016; Legras and Bucci, 2020). Bergman et al. (2013), using Lagrangian backward trajectories, showed that the anticyclone is connected to the boundary layer through a vertical conduit centred over northeastern India, Nepal, and southern Tibet. In the recent BATAL campaign, Vernier et al. (2018) used back-trajectory calculations to point at the north of India as a principal region source for the ATAL. Lau et al. (2018), based on MERRA-2 reanalysis, reported that the Himalayas Gangetic Plain (HGP) region and the Sichuan Basin (SB) of southwestern China are two important regions with strong vertical transport of CO, carbonaceous aerosols and dust from the surface to the UTLS. On the other hand, the simulations of Fairlie et al. (2020) suggested that the anthropogenic sources from India contribute up to 40 % of sulfate and up to 65 % of organic and ammonium aerosols in the western ATAL region, whereas China contributes up to 60 % (both sulfate and organic aerosols) in the eastern ATAL region.
It is also important to note that the ATAL formation and possible spatial and temporal variability are closely related to the dynamical variability of the AMA. For example, Basha et al. (2020) suggested that the spatial extent and strength of the AMA are greater during July and August compared to June and September and that the decadal variability is bigger at the edges of the anticyclone. As a consequence of the variability of atmospheric dynamics, some years show a stronger monsoon activity than others (Lau et al., 2018; Basha et al., 2020; Yuan et al., 2019), and this affects the ATAL formation, location and composition. Several studies have shown that the AMA exhibits intraseasonal variability between the Iranian Plateau and the Tibetan Plateau with a quasi-biweekly oscillation (e.g., Zhang et al., 2002; Yan et al., 2011; Nützel et al., 2016; Pan et al. 2016; Wei et al., 2019).
This study provides further insight into the chemical composition of the ATAL and assesses its decadal variability composition and aerosol trends for the first time. To assess this, we have carried out long-term modeling of the ATAL using the Community Earth System Model (CESM 1.2), which embeds the Community Atmosphere Model (CAM5) coupled with the MAM7 (Modal Aerosol Model) aerosol module. Our simulations cover an overall extended period of 16 years, from 15 January 2000 to 15 December 2015. Yuan et al. (2019) derived decadal trends for carbonaceous aerosols and dust in the ATAL using only meteorological reanalysis data, while in the present study a detailed chemistry and microphysical modeling is used to estimate trends for a more comprehensive set of aerosol compositions.
The present paper is structured as follows. In Sect. 2, we describe the model and correlative data used for its validation. The validation is discussed in Sect. 3. Results are presented and discussed in Sect. 4. Conclusions are drawn in Sect. 5.
Model simulations were performed using the global Community Earth System Model (CESM1.2), based on the Community Atmospheric Model (CAMS 5.1) with its full chemical core for both troposphere and stratosphere, coupled with the Modal Aerosol Model (MAM7). The MAM7 module treats the aerosol microphysics, size distribution and both internal and external mixing using seven modes. The seven modes are, specifically, accumulation (a1), Aitken (a2), primary carbon (a3), fine dust and sea salt (a5 and a4), and coarse dust and sea salt (a7 and a6) (Liu et al., 2012). Extraterrestrial aerosols are neglected in our model. Table 1 lists the aerosols and dry diameter size ranges of each mode. The size distributions of each mode are assumed to be log-normal.
Predicted species for interstitial and cloud-borne components (see text) of each aerosol mode in MAM7 and dry diameter size ranges.
The total number of transported aerosol tracers by the seven log-normal modes in MAM7 is 31. The transported precursor gas species are SO
Wet removal of soluble gas-phase species combines two processes: in-cloud,
or nucleation scavenging (rainout), which is the local uptake of soluble
gases and aerosols by the formation of initial cloud droplets and their
conversion to precipitation, and below-cloud, or impaction scavenging
(washout), which is the collection of soluble species from the interstitial
air by falling droplets or from the liquid phase via accretion processes.
The transfer of soluble gases into liquid condensate is calculated using
Henry's law, assuming equilibrium between the gas and liquid phases. This is the standard scheme used in CAM5.1 (Lamarque et al., 2012), although as noted by Fairlie et al. (2020), a more physically based treatment of wet scavenging of SO
The MAM7 module explicitly treats the microphysics of sulfate (SO
In MAM7 the aerosol particles (AP) can exist in the “interstitial” state (AP that are suspended in clear or cloudy air) and “cloud-borne” state (AP attached to or contained within different hydrometeors, such as cloud droplets and/or ice crystals). MAM7 distinguishes between cloud-borne aerosols that are within stratiform clouds and the interstitial aerosols which include both clear-sky AP and AP contained within convective clouds. This means that the AP in convective cloud droplets are lumped with the interstitial AP in the model, and the interstitial aerosol mixing ratios include the truly interstitial (i.e., “clear-sky/dry”) AP and the “convective” cloud-borne AP.
As has been detailed in Wang et al. (2013), in CAM5-MAM7 cloud-borne aerosols in stratiform clouds are treated in a prognostic way in CAM5: their mixing ratios are saved between model time steps and evolve as a result of source, sink, and transport processes. Their activation is parameterized using vertical velocity (resolved and sub-grid turbulent) and aerosol properties of all the modes, following Abdul-Razzak and Ghan (2000). The stratiform-cloud-borne AP are assumed to not interact with convective clouds. AP in convective clouds are treated diagnostically: their mixing ratios are diagnosed each model time step (with no “memory”) from the interstitial aerosol mixing ratios.
Both interstitial and cloud-borne aerosol particles are subject to wet and dry removal deposition. CESM-MAM7 distinguishes between “in-cloud” and “below-cloud” wet removal. In-cloud wet removal involves activation of interstitial aerosol to become cloud-borne, followed by conversion of cloud droplets (and the cloud-borne aerosol particles) to precipitation. Below-cloud wet removal involves direct capture of interstitial aerosols by precipitation particles through a number of processes (e.g., inertial impaction, Brownian diffusion) and is relatively inefficient for aerosol in the accumulation mode size range. For a complete description of the CAM5-MAM7 model, see Liu et al. (2012).
In our configuration, land, sea ice, and rivers are interactive processes in CESM, whereas oceans are prescribed. The model horizontal grid resolution is
The following emissions are used in our simulations. Biogenic emissions for
CO, isoprene, C
The biomass burning emissions for CMIP6 are based on merged satellite observation and fire models (van Marle et al., 2017), using GFED4 (Global Fire Emissions Database version 4), which includes small-magnitude fires (available from 1997 to 2015).
The emission of sea salt aerosols from the ocean follows the
parameterization of Märtensson et al. (2003), for aerosols with
a geometric diameter
The meteorology in the model has been nudged using MERRA2 (Modern-Era
Retrospective analysis for Research and Applications, Version 2,
In the standard configuration of CESM-MAM7 the vertical transport of interstitial aerosols and trace gases by deep convective clouds uses updraft and downdraft mass fluxes from the Zhang–McFarlane parameterization (Zhang and McFarlane, 1995). Currently this vertical transport is calculated separately from wet removal. Cloud-borne aerosols associated with large-scale stratiform clouds are assumed to not interact with the convective clouds.
Vertical transport by shallow convective clouds is treated similarly, using mass fluxes from the Park and Bretherton (2009) shallow convection parameterization.
We run our simulations for 16 years, from 15 January 2000 to 15 December 2015, using the CESM1.2 (CAM5) initial atmosphere state file at that date.
Our simulations have been compared to satellite data from the Microwave Limb Sounder (MLS) and the Atmospheric Chemistry Experiment–Fourier Transform Spectrometer (ACE-FTS).
The MLS was launched in July 2004 onboard the NASA Aura satellite. Measurements in the millimeter and submillimeter wavelength ranges are continuously made during both night and day every 165 km along the
suborbital track, covering latitudes from 82
The ACE-FTS instrument is an infrared solar occultation spectrometer, providing profiles of the Earth since February 2004 from the Canadian
satellite SCISAT-1 (Bernath et al., 2005). It operates in the wavelength
range from 2.2 to 13.3
We compare CO measurements from MLS and ACE-FTS with our simulations. While a direct comparison of aerosol extinction observations from various satellite instruments with CESM-MAM7 is not easy, e.g., due to the interference of clouds, using a trace gas (like CO) is a more straightforward approach for a comparison. In fact, three-dimensional summer distributions of CO show a distinct enhancement in the AMA and have proven an ideal tracer to identify the AMA's location and to track the transport processes to the AMA (e.g., Park et al., 2008; Barret et al., 2016; Santee et al., 2017). The CO comparison enables a test of the model's capacity to reproduce the large-scale dynamical and morphological features, which is related to the aerosol distribution.
Figure 1a and b show the average summer (June–July–August) CO distribution, for the year 2008, observed by MLS in the UTLS (Fig. 1a) and produced by CESM-MAM7 (Fig. 1b), at the pressure level of 147 hPa, for MLS, and 150 hPa (average between 160 and 140 hPa, three levels), for CESM-MAM7. The locations of the general enhancement of CO mixing ratios in the AMA and of the absolute maximum above India are well reproduced by the model (i.e., they are consistent with MLS observations). It should be noted that the pressure levels used in this comparison, for CESM-MAM7 and MLS, are not exactly identical. In addition, the vertical resolutions differ as well (about 5.1 km, for MLS, and about 1 km, for CESM-MAM7). Furthermore, the temporal samplings of satellite and model data also differ: for CESM-MAM7 the temporal sampling is twice a day (noon and midnight), whereas MLS samples the Earth on distinct orbits, with a full global coverage every 3 d. Even though it is therefore possible that intensive short-time events are missed by either CESM-MAM7 or MLS, the sampling bias is not expected to present a significant source of discrepancies for 3-month averages, as shown in Fig. 1a and b.
Compared to MLS observations, the model underestimates the CO mixing ratio by about 30 %. One possible reason for this underestimation could be an underestimation of biomass burning emissions in the model (obtained from GFED4), which are a significant source of CO. We have also compared CESM-MAM7 HCN mixing ratios (a strong biomass burning tracer) with ACE-FTS HCN observations (comparison not shown here). This latter comparison shows a marked underestimation of modeled HCN amounts, which supports the hypothesis of an underestimation of biomass burning emissions. Stroppiana et al. (2010) compared different biomass burning inventories for CO. For 2003, they found that the CO emissions range from 365 Tg (GFED3) to 1422 Tg (VGT – Vegetation Emission Inventory, CNRS-LA) (Tansey et al., 2008), with GFED at the low end of this variability. Unlike GFED3, GFED4 includes upgrades like the inclusion of small fire-burned areas and a revised fuel consumption parameterization that causes global emissions to increase in comparison with the previous version. However, the effects of these adjustments vary spatially and, in particular regions like the southeast of Asia or the north and south of Africa, the CO biomass burning emissions are lower (see Van de Werf et al., 2017). This could explain the low bias in CO mixing ratios for our comparisons with satellite measurements. On the other hand, as mentioned in Sect. 2.1, the CEDS anthropogenic inventory uses a scaling process to match the CEDS emissions estimates with available inventories. In the case of anthropogenic emissions for CO, the last year from local inventories available is 2008 in Asia (from REAS) and 2010 in China (from the MEIC-Multi-resolution Emission Inventory for China). As a result, the extrapolation during 2010–2015 may be an additional source of uncertainties for comparisons with observations over this period.
While reproducing monthly average features is a probing test for our simulations, catching shorter-term processes and variability is even more challenging towards the description of a complex phenomenon as the ATAL. Thus, we have also tested the model's ability to reproduce observed daily specific features. Figure 1c shows a 3 d average from 4 to 6 July 2005. During this short time period, a multi-centric AMA is observed by MLS, with rather multiple maxima in eastern Asia, instead of a classical individual maximum above the Himalayan region. Our CESM-MAM7 simulations reasonably reproduce this pattern. They show a distributed pattern with maxima above eastern Asia but also above western Asia (Fig. 1d), which is very consistent with MLS observations (Fig. 1c). For 3 d averages the sampling bias can play a significant role in explaining the different patterns observed for MLS and CESM-MAM7. Therefore, some short-term features might not have been captured by the MLS instrument. Nevertheless, our simulations are very consistent with MLS observations for this short-term configuration.
We have also tested the vertical structures of CESM-MAM7 simulations, using
an ACE-FTS CO mixing ratio profile in the UTLS (Fig. 1e). Observations with
ACE-FTS have been chosen because of their better vertical resolution with
respect to MLS. It has to be noted that the location and time of the ACE-FTS
measurement profile and the model output are not exactly the same but agree within
Because of the sparse sampling of ACE-FTS data in the AMA, we have provided an additional comparison of the monthly vertical distribution of CO for the whole year of 2008 between MLS data and modeled CO (Fig. S1). The comparison, while showing an underestimation of the modeled vertical amounts of CO, especially below the 150 hPa level, presents spatial distributions of CO which are in good agreement.
The comparison of simulated CO with observed MLS and ACE-FTS CO in the UTLS allows us to conclude that, except for a possible underestimation of CO emissions, the model is able to reproduce the position and spatial extent of the Asian monsoon anticyclone in our simulations.
Figure 2 shows the CESM-MAM7 regional distribution, over an extended area centered around the AMA region, of different aerosol types: sulfate, SOA, POM, BC, ammonium and mineral dust. The accumulation mode (a1) is shown here for all aerosol types, except for mineral dust (for which fine soil dust mode (a5) is shown). These maps represent average aerosol concentrations, for July–August 2014, at three different pressure levels: 120, 100 and 80 hPa, respectively (approximately 15.0, 16.5 and 18.0 km). Concentrations of sea salt particles, also modeled in our study, are negligible and therefore are not shown in Fig. 2. The model reproduces the horizontal distribution of the ATAL, i.e., an increase in aerosol concentration in the AMA region with elevated aerosol concentration at 120–100 hPa (upper troposphere) and noticeably decreasing for pressures lower than 80 hPa (altitudes higher than 18.0 km, lower stratosphere).
Figure 2 shows that dust is the principal aerosol species in the ATAL, in
terms of mass concentration, in our simulations. These results agree with
some previous modeling studies (e.g., Fadnavis et al., 2013; Ma et al., 2019). Our results show an aerosol dust concentration at 100 hPa of about 100 ng m
In addition, Wu et al. (2019), using CESM1-CAM5 with the default scheme for the dust emissions (Zender et al., 2003), the same scheme used in the present study, showed that the model overestimates dust extinction over the Taklamakan and Gobi deserts during the summer period. Such high biases in dust extinction have been attributed to excessive convective transport, lack of secondary activation of aerosol entrained into convective updrafts and strong dust transport in the upper troposphere from Africa and the Middle East. These hypotheses, together with differences in the model resolution, could explain the higher dust amounts in our CESM-MAM7 simulations, which use the same default scheme for the generation of dust.
The discrepancies observed between different models could also result from the different schemes used for the dust lifting as well as the sensitivity of dust release to surface conditions, particularly to surface winds and soil properties.
Other main aerosol components contributing to the ATAL in our simulations are sulfates, followed by SOA, POM, ammonium and to a lesser extent BC. Yu
et al. (2015), using the CESM1/CARMA model, suggested that the ATAL (at altitude levels between 230 and 100 hPa) is principally composed of organics (
As discussed in previous studies, the spatial extent, strength and position
of the AMA are highly variable due to the dynamical seasonal variability of the ASM (e.g., Randel and Park, 2006; Garny and Randel, 2013; Lau et al., 2018; Basha et al., 2020). Due to this dynamical variability, the tracer concentrations are strongly controlled by the oscillations and shedding of
the AMA that therefore affect the ATAL extent and composition. In order to determine the aerosol burden within the ATAL, we have defined a simple criterion to isolate the ATAL horizontal extent, i.e., where there is a high probability of finding AMA air masses, based on a threshold on the geopotential
height (GPH) values. Similar empirical selections of high GPH values to
represent anticyclone boundaries have been used in a number of previous
works, e.g., Highwood and Hoskins (1998), Bergman et al. (2013), Barret et al. (2016), and Pan et al. (2016). For the subsequent analysis, we identify the AMA region based on GPH values higher than 16.7 km at 100 hPa. Based on this
criterion, a wide region from around 20–130
Spatial distribution of the aerosol mass concentration, averaged over July–August 2014, from CESM-MAM7 simulations, for six
different aerosol types. From top to bottom row: sulfate, SOA, POM, BC,
ammonium (in the accumulation mode) and mineral dust (in the fine dust soil
mode). From left to right column: 120, 100 and 80 hPa pressure levels. Note
the different color-scale ranges. The black lines in the map represent the geopotential height
In Fig. 3 we show CESM-MAM7 vertical aerosol mass mixing ratio profiles for the accumulation mode, averaged from June to August within the blue box of Fig. 2, for 2 selected years, 2000 and 2014. Our focus on the accumulation mode is justified by the fact that it is the principal mode that contributes to the ATAL (see Fig. S2 in Supplement), with mostly anthropogenic origin. In this first analysis, we have excluded dust. Dust is still the most important ATAL component, in our simulations, in terms of mass, but its burden and variability are mostly subject to natural factors and their variability.
A vertical region with marked localized increase in the concentrations of all the aerosol types is observed between 300 and 80 hPa. This is what is
expected as a manifestation of the ATAL, as it is broadly the vertical
region where the AMA is located. The vertical structure of the AMA-related
dynamics has been investigated by several authors (e.g., Park et al., 2009; Bergman et al., 2013; Garny and Randel, 2013; Brunamonti et al., 2018; Bian
et al., 2020), showing evidence of deep convection and confinement extending
up to 1.5–2.0 km above the cold-point tropopause. Enhanced aerosol
backscatter also reveals the signature of the ATAL over the same altitude
range (Vernier et al., 2015; Brunamonti et al., 2018). This location
suggests that the existence of the layer is tied to a large-scale vertical
transport in the anticyclone, i.e., around 200 to 80 hPa (
Our simulations show a characteristic “double-peak” vertical configuration
with two relative maxima, one at higher altitudes (
This “double-peak” vertical structure can be found in some observations
from recent aircraft and balloon campaigns but is not discussed. For instance, particle-counting observations during the 2015 BATAL campaign (Vernier et
al., 2018) have shown two maxima in the aerosol concentration profile, at
We have also tried to separate the overall in-cloud and purely dry aerosols (these latter ones likely coming from nucleation/condensation
processes). In order to analyze the contribution of dry aerosols to the ATAL, we have carried out an analysis to reduce the contribution from convective
cloud-borne aerosols. For this purpose, we have filtered out the profiles,
in our blue box, for which the extinction coefficient is larger than an
arbitrary threshold (
It is worth noticing that, for these 2 selected years (2000 and 2014), the aerosol profiles can differ from one aerosol type to another but are quite similar for a given month/year, and a double- or single-peak structure is observed. This variability observed in the ATAL's vertical profiles also reflects the aspect of the dynamical variability of the AMA, which can be put into a relationship with both the long-range transport and convection, as was shown in previous studies (e.g., Qie et al., 2014; Pan et al., 2016; Santee et al., 2017).
Modeling vertical profiles of aerosol mass concentration of sulfate, SOA, POM, ammonium and BC in the accumulation mode (a1) in ng/kg
averaged between 20–35
Figure 4a–d show the annual average aerosol total mass concentrations for all the aerosol types simulated by CESM-MAM7, in the period 2000–2015, for all modes (Fig. 4a, c) and the isolated accumulation mode (Fig. 4b, d). To account for the whole double-peak phenomenology and to isolate the single dry AP peak (see discussion in Sect. 4.2), the concentrations are averaged between 200–80 hPa (Fig. 4a, b) and 120–80 hPa (Fig. 4c, d). These two vertical ranges allow the differentiation of the ATAL composition based on in-cloud processes or, from another point of view, description of how the composition changes depending on the altitude. No filter has been applied to show the contributions of all aerosols.
The aerosol type that dominates the ATAL, for both altitude ranges, is dust, followed by sulfates and organic particles (secondary and primary). The comparison between Fig. 4a and c shows that at higher altitudes the amount of sulfates increases slightly and, more markedly, dust amount decreases. Figure 4e shows the percent contribution of aerosols types to the ATAL, between 120 and 80 hPa. It is evident that although less dust reaches higher altitudes, this aerosol type is still the mass-dominant aerosol type in the ATAL, contributing around 60 %. Even if there still is large disagreement among reported studies about the exact amount of dust present in the ATAL, it is clear that in our study this natural component contributes significantly to the ATAL seasonal build-up due to its transport from the nearby desert regions, like the Taklamakan and Thar deserts, and the northern slope areas of the Tibetan Plateau (Lau et al., 2018; Ma et al., 2019). As was mentioned in Sect. 4.1, the difference in the amount of dust reported by the different authors may be related to the different schemes used for the generation of dust, e.g., how the topography is represented in the model, the resolution of the model and the parameterization of the convection processes.
Wu et al. (2018) evaluated dust emissions in East Asia simulated by 15 climate models participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) during 1961–2005. They found discrepancies with the observations for all the models, because climate models may not sufficiently represent the trends of surface wind speeds and precipitation. This indicates that there is still a need to improve the representation of the dust cycle in climate models to simulate long-term dust changes.
With the intention of analyzing the composition of the ATAL in terms of anthropogenic and biomass burning emissions, we discuss more in detail the
contribution of the non-dust aerosols, for which the accumulation mode at
two different altitude ranges is shown in Fig. 4b, d. Excluding dust
particles, the accumulation mode (a1, size range: 0.056–0.26
Sulfate aerosols from moderate-to-strong volcanic eruptions, with injection
in the UTLS, can also interact with the dynamical features of the AMA (e.g., Sellitto at al., 2017) and, under certain conditions, can impact the ATAL aerosol population. Larger sulfate concentrations in 2009 and 2011 are
linked to the volcanic eruptions of Sarychev (June 2009) and Nabro (June
2011). These eruptions injected large quantities of SO
Evolution of the total aerosol mass concentration of all
the aerosol types present in CESM-MAM7 in all the modes averaged at
20–35
Excluding dust and focusing on the mostly anthropogenic accumulation mode,
Fig. 4f suggests that the fraction of the ATAL of anthropogenic origin is composed of about 40 % sulfate, 30 % SOA, 15 % POM, 14 % ammonium
and less than 3 % BC. Compared to the results reported by Yu et al. (2015), our results show about the same percentage of sulfate in the ATAL
but less organics, i.e.,
In the following, we evaluate the decadal trends of the different aerosol types in the ATAL. In particular, we have estimated the trends for the dust in the fine soil dust mode (Fig. 5a) and all other aerosol types in the accumulation mode (Fig. 5b, c). The concentrations for each year are averaged between 120 and 80 hPa pressure levels and over the domain defined by the blue box of Fig. 2, excluding the years with volcanic eruptions impacting the UTLS, i.e., 2005–2009 and 2011–2012 (Manam: April 2005; Soufrière Hills: August 2006; Tavurvur: October 2006; Okmok: August 2008; Kasatochi: August 2008; Sarychev: June 2009; Nabro: June 2011, taken from Khaykin et al., 2017; see Table 3 in their paper).
As can been seen from Fig. 5a, dust does not display any clear trend. The
Figure 5b, c show the trends for all the aerosols in the accumulation mode
averaged in our box over the 120–80 hPa vertical level range, respectively, without and with the extinction filter applied so as to isolate dry from
in-cloud (including from convective clouds) aerosols. All the aerosol types
show an increase over the simulated 16 years. This mirrors the increase in the emissions in Asia. From Fig. 5b, it can be seen that sulfate aerosol trends in the ATAL roughly double their concentration from
Figure 5c shows the trends of dry aerosols, i.e., with the extinction filter applied, in the ATAL between 120 and 80 hPa. The comparison between Fig. 5c with b, together with the values reported in Table 2, shows that the increasing trends and correlation values are slightly smaller than values reported without applying the filter. This reflects the fact that at 120–80 hPa the dry aerosols contribute a larger fraction of the ATAL than convective cloud-borne aerosols.
The analysis of differences without and with the application of the
extinction filter (i.e., (dry
We have also carried out the same analysis for the larger altitude interval
of the ATAL, i.e., between 200 and 80 hPa (Fig. 5d and e). More convective cloud-borne aerosols are present in this case. Thus, the differences for the
cases without versus with the extinction filter (calculated from Table 2) are
larger than the previous case (
Averaged aerosol mass concentration and percentage of the
increase from 2000 to 2015 for SO
Aerosol mass concentration trends simulated by CESM-MAM7
averaged between 20–35
Figure 6a and b show the aerosol optical depth (AOD) at 550 nm averaged for July–August and between 20 and 35
AOD values from about 0.0019, in 2000, to about 0.004, in 2015, are found over the 120–80 hPa range where dry aerosols dominate (Fig. 6b). Between 200 and 80 hPa, higher AOD values are obtained as a result of a large contribution of convective cloud-borne aerosols at this altitude range. The difference between the AOD values obtained for the two altitude ranges in Fig. 6a and b points to the importance of what we have identified as convective in-cloud aerosols.
Figure 6c and d show the temporal evolution of the yearly ATAL AOD from 2000 to 2015 for our selected box (20–35
The attribution of the possible causes to the increase in the aerosol content and optical depth in the ATAL between 2000 and 2015 (e.g., increase in Asian emissions, more efficient vertical transport or different chemical/microphysical processes) requires further investigations and the continuous monitoring of ATAL burden and properties in the future.
In this paper, we have presented the results for our long-term simulation, i.e., 16 years (15 January 2000–15 December 2015), to investigate the composition and trends of the specific ATAL aerosols using the CESM-MAM7 model. The model was driven by the CMIP6 emission inventory for the anthropogenic and biomass burning emissions of the principal trace gases and aerosols, while the biogenic emissions were taken from the MEGAN-MACC inventory. During summer, a confinement of polluted air masses has been found within the AMA region, which is tied to the ATAL position. The model results show overall good agreement with the space–time behavior of CO in the UTLS region observed by the MLS and ACE-FTS space-borne instruments, despite a possible underestimation in the CO burden due to the underestimation of surface emissions. In particular, the horizontal distribution of modeled CO is in good agreement with MLS data, and the vertical structure in the AMA shows a maximum near 150 hPa in agreement with the available ACE-FTS observations.
Our model results indicate that dust is a dominating aerosol type in terms of mass in the ATAL, in agreement with other studies (e.g., Lau et al., 2018; Ma et al., 2019). However, the lack of in situ or satellite measurements of dust in the AMA region makes the validation of this result difficult. Our modeled burdens of dust in the ATAL are larger than what has been reported in the past (e.g., Fadnavis et al., 2013; Yu et al., 2015; Fairlie et al., 2020). The higher amount of dust found in our model could be due to excessive convective transport, a lack of secondary activation of aerosols entrained into convective updrafts, a too strong dust transport in the upper troposphere from Africa and the Middle East (Wu et al., 2019), as well as the sensitivity of dust emissions to the resolution of the model (Brühl et al., 2018; Wu et al., 2019).
The differences between the simulated dust burdens between different models can be linked to the different physical processes computed for dust emissions (e.g., wind speed, hydrological parameters and soil properties).
Apart from dust, the average partitioning for other aerosol types contained in the ATAL (from anthropogenic and biomass burning emissions) is the following: 40 % sulfate, 30 % secondary organic aerosols, 15 % primary organic matter, 14 % ammonium and less than 3 % black carbon. Nitrate aerosols are expected to be an important aerosol component in Asia (e.g., Höpfner et al., 2019) due to the increase in nitrogen oxides and ammonia emissions but are not simulated in our work.
For non-dust aerosols the accumulation mode dominates the anthropogenic and
biomass burning ATAL aerosols. A marked positive trend of anthropogenic and
biomass burning aerosol concentrations is found, with up to a factor of 2 increase in mass concentrations between 2000 and 2015. It is important to
note that the simulated aerosol trends depend on the emission inventory used. For example, Zheng et al. (2018) showed that after 2013 the SO
Our simulations reveal a double-peak structure in the vertical profile of aerosols of the ATAL, highlighting the contribution of two types of aerosols, i.e., “cloud-borne” aerosols, including those from convective clouds and “clear-sky/dry” aerosols. The CESM-MAM7 simulations have allowed us to analyze separately the contributions of these two types of aerosols. Dry aerosols contribute to the higher peak (peaking around 80–120 hPa) and convective cloud-borne aerosol to the lower peak (peaking around 200–250 hPa). We show that the contribution of the convective cloud-borne aerosols to the ATAL generally increases during the phases of the mature and late ATAL, in July–August, shifting the maximum of aerosol concentrations to lower altitudes. The dry aerosols are generally dominating in the early phases of the ATAL. This “double-peak” vertical structure has been observed in recent balloon and aircraft campaigns (e.g., Vernier et al., 2018; Höpfner et al., 2019) but has not been discussed in detail so far. These observations support our simulation results, which in turn provide a possible explanation for the observations. Given the uncertainties discussed throughout the paper, the ability of our simulations to represent the reality of the convective transport in the ASM is not entirely clear, but the model results provide hypotheses for follow-up studies.
The obtained AOD values show an enhancement by a factor
The release version 1.2.2 of CESM can be downloaded from
Atmospheric Forcing data, regridded from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA2) data and prepared as meteorological variables to run CESM and WRF simulations
are available at
The supplement related to this article is available online at:
AB, PS, GB and FJ designed the research and the analysis and interpretation of the model results. AB performed the model simulation with the support of FJ. CK performed the satellite analysis from MLS and ACE-FTS data. BL was involved in the discussion and results interpretation. AB prepared the manuscript with contributions and discussions from all the co-authors.
The authors declare that they have no conflict of interest.
This article is part of the special issue “StratoClim stratospheric and upper tropospheric processes for better climate predictions (ACP/AMT inter-journal SI)”. It is not associated with a conference.
The authors wish to thank the CaSciModOT structure (Calcul Scientifique et Modélisation Orléans-Tours), part of the French national network of complex systems (RNSC – Réseau National des Systèmes Complexes), thanks to which the simulations could be completed.
Adriana Bossolasco would also like to thank the NCAR/CESM online discussion board for many helpful technical discussions that helped throughout this study, and thanks especially to Louisa Emmons and Simone Tilmes.
Furthermore, the authors thank the ACE-FTS and MLS teams.
This research has been supported by the Agence Nationale de La Recherche (grant no. ANR-17-CE01-0015, TTL-Xing), the VOLTAIRE project (ANR-10-LABX-100-01) funded by ANR through the PIA (Programme d'Investissement d’Avenir), and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation, grant no. 409585735).
This paper was edited by Gabriele Stiller and reviewed by three anonymous referees.