The Aerosol Module in the Community Radiative Transfer Model (v2.2 and v2.3): accounting for aerosol transmittance effects on the radiance observation operator

The Community Radiative Transfer Model (CRTM), a sensor-based radiative transfer model, has been used within the Gridpoint Statistical Interpolation (GSI) system for directly assimilating radiances from infrared and microwave sensors. We conducted numerical experiments to illustrate how including aerosol radiative effects in CRTM calculations changes the GSI analysis. Compared to the default aerosol-blind calculations, the aerosol influences reduced simulated brightness temperature (BT) in thermal window channels, particularly over dust-dominant regions. A case study is presented, which illustrates how failing to correct for aerosol transmittance effects leads to errors in meteorological analyses that assimilate radiances from satellite infrared sensors. In particular, the case study shows that assimilating aerosol-affected BTs significantly affects analyzed temperatures in the lower atmosphere across several regions of the globe. Consequently, a fully cycled aerosol-aware experiment improves 1–5 d forecasts of wind, temperature, and geopotential height in the tropical troposphere and Northern Hemisphere stratosphere. Whilst both GSI and CRTM are well documented with online user guides, tutorials, and code repositories, this article is intended to provide a joined-up documentation for aerosol absorption and scattering calculations in the CRTM and GSI. It also provides guidance for prospective users of the CRTM aerosol option and GSI aerosol-aware radiance assimilation. Scientific aspects of aerosol-affected BT in atmospheric data assimilation are briefly discussed.

1 Introduction

An accurate and computationally efficient radiative transfer model is essential in radiance assimilation for supporting weather prediction, physical retrievals for satellite environmental data records, and inter-comparison between different remote sensing instruments. The Community Radiative Transfer Model (CRTM) is a radiative transfer model used extensively within satellite and remote sensing systems (Weng, 2007; Han et al., 2007). It was primarily designed for computing satellite radiances and has been widely used within the Gridpoint Statistical Interpolation (GSI, Wu et al., 2002; Kleist et al., 2009) system for directly assimilating radiances from infrared (IR) and microwave (MW) sensors. Specifically, clear-sky radiance calculations are carried out within the CRTM given the atmospheric scattering and absorption profile, surface emissivity and reflectivity, and source functions. For cloudy radiance simulations (Stegmann et al., 2018), vertical profiles of hydrometeor variables (e.g., cloud liquid water path and ice water path) are also required. Note that the CRTM was not designed to enact composition–radiation interaction effects within spectral longwave and shortwave radiative transfer calculations in general circulation models. Instead, the CRTM was developed to support monochromatic satellite radiance assimilation from longwave and microwave sensors, and for satellite retrieval algorithm development.

Past studies have demonstrated that aerosols significantly impact the simulation of brightness temperature (BT) in the IR channels. BT is “a descriptive measure of radiation in terms of the temperature of a hypothetical blackbody emitting an identical amount of radiation at the same wavelength” (American Meteorological Society, 2012). A reduction in retrieved BT of 2–4 K in the atmospheric window region due to a strong dust outbreak was reported during the Saharan Dust Experiment (SHADE) campaign (Highwood et al., 2003). Pierangelo et al. (2004) and Peyridieu et al. (2009) showed that the dust cooling effects may reach 3 K in tropical atmospheric conditions depending on the dust burden. Diaz et al. (2001) found that there is a significant increase in the errors of sea surface temperature (SST) retrievals in the presence of enhanced aerosol loading in the atmosphere. The dust effects on satellite-derived SST are constrained by accounting for dust absorption (Weaver et al., 2003), applying a dust correction scheme (Nalli and Stowe, 2002; Merchant et al., 2006), or removing dust-contaminated observations (Divakarla et al., 2012).

The impact of aerosol-affected BTs on the meteorological analysis fields has also been investigated. Wei et al. (2021) used the Global Data Assimilation System (GDAS, Kleist et al., 2009) to assess the aerosol impact on the meteorological analysis. To do this, two GDAS experiments were conducted: a control cycled experiment, where aerosol transmittance effects are not considered, and an offline non-cycled experiment, where aerosol transmittance effects are considered in the BT calculations. The offline experiment uses identical observations and first guesses as the control experiment, and thus the response of atmospheric analysis to aerosol-aware radiance calculations can be clearly demonstrated. The experimental setup in Wei et al. (2021) followed the methodology presented in Kim et al. (2018), which is based on the Goddard Earth Observing System (GEOS) atmospheric data assimilation system (ADAS). Note that GEOS-ADAS and GDAS both used GSI and CRTM, although the version and configuration differed. The studies by Kim et al. (2018) and Wei et al. (2021) reported that (i) there is a considerable cooling effect on simulated BT when aerosols are considered, (ii) including aerosol transmittance effects in the BT calculation improves the fit to observations over the dust-laden regions, and (iii) the offline aerosol-aware experiment produces warmer analyzed SST (0.3–0.5 K) over the Atlantic Ocean. Wei et al. (2021) also reported a warmer analyzed lower atmosphere (0.15 K) over Africa and the central Atlantic Ocean in the offline aerosol-aware experiment.

The experiments conducted in Kim et al. (2018) and Wei et al. (2021) were based on the application of the CRTM aerosol absorption and scattering routines. While aerosol absorption and scattering options are available from CRTM version 2.2 onwards, to our knowledge the documentation of the CRTM aerosol module (Liu and Lu, 2016) has yet to be updated. Here we presented a joined-up documentation for aerosol absorption and scattering calculations in the CRTM and GSI. In addition, we provide guidance for prospective users of running aerosol-affected GSI analysis. Scientific aspects of aerosol-affected BT in atmospheric data assimilation are also briefly discussed.

2 GSI and CRTM

Below, we provide a brief introduction to the GSI in Sect. 2.1 and a description of the CRTM aerosol option in Sect. 2.2. In Sect. 2.3, a description of running aerosol-aware GSI analysis is given.

2.1 GSI

The multi-partner-developed GSI is an incremental three-dimensional variational (3D-Var) data assimilation system (Wu et al., 2002; Kleist et al., 2009). GSI, alone or combined with an ensemble system, has been used widely by modeling centers and the research community for a range of applications. For instance, it is used operationally by the National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Prediction (NCEP) for medium-range weather forecasting. It is also used by the National Aeronautics and Space Administration (NASA) Global Modeling and Assimilation Office (GMAO) for recent production of the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2; Gelaro et al., 2017). The community version of the GSI system has been supported and maintained by the Developmental Testbed Center (DTC; http://dtcenter.org, last access: 8 February 2022). Note that DTC is scheduled to cease all activities supporting the GSI user community by the end of December 2021. However, community GSI-related assets (website, forum, and repository) built by DTC will remain available to and usable by the community.

GSI can assimilate a wide range of observations, including conventional observations (such as radiosonde observations), radar data, satellite retrievals (for example global positioning system (GPS) radio occultation sounding data), satellite radiance data, etc. For IR satellite instruments, GSI has the capability to assimilate radiances from the Advanced Infrared Sounder (AIRS) on AQUA; Infrared Atmospheric Sounding Interferometer (IASI) on METOP-A and METOP-B; Cross-track Infrared Sounder (CrIS) on S-NPP; High resolution Infrared Radiation Sounder (HIRS) on METOP-A, METOP-B, and NOAA-19; Advanced Very High Resolution Radiometer (AVHRR) on NOAA-18 and METOP-A; Spinning Enhanced Visible and Infrared Imager (SEVIRI) on M08 and M10; and Geostationary Operational Environmental Satellite (GOES) Sounders (sndrD1, sndrD2, sndrD3, and sndrD4) on GOES-15. A comprehensive list of all observations assimilated and monitored by GDAS can be found at the web page for “Observational Data Processing at NCEP” (https://www.emc.ncep.noaa.gov/emc/pages/infrastructure/obs-data-processing.php, last access: 8 February 2022).

Despite the broad applications of GSI, the publicly released version handles only clear-sky radiances for IR sensors. Without correcting for aerosol transmittance effects, systematic biases may be introduced into the meteorological analysis fields when observations affected by aerosols are assimilated. The aerosol-aware option (discussed in Sect. 2.2) reduces such errors by enabling aerosols to influence GSI's radiance observation operator, CRTM, which calculates the BT and Jacobians (radiance first derivative). This option, however, may fluctuate the amount of observations assimilated in GSI because the quality control (QC) algorithm screens out observations based on measured BTs and aerosol-free simulated BTs. Thus, an improved QC algorithm is needed to fully exploit radiance measurements under all sky conditions. The technical issues regarding the QC procedure have been discussed in Kim et al. (2018) and Wei et al. (2021).

2.2 CRTM aerosol module

The CRTM, a one-dimensional radiative transfer model (Liu and Weng, 2006), is being developed at the US Joint Center for Satellite Data Assimilation (JCSDA) with algorithm and software input from JCSDA collaborating research institutions. The CRTM is composed of four modules, which include gaseous transmittance, surface emission and reflection, cloud and aerosol absorption and scattering, and a solver for radiative transfer (Han et al., 2006). Given an atmospheric profile of temperature, cloud and surface properties, and gaseous constituents and aerosol concentrations, the CRTM is called within the GSI to calculate BTs for satellite sensors from IR sounders to MW imagers. Here, we describe the aerosol scattering and absorption scheme in CRTM version 2. We refer the readers to Han et al. (2006) for the full details regarding CRTM version 1.

Absorption by atmospheric trace gases, such as water vapor and carbon dioxide, is parameterized using the Optical Depth in Absorber Space (ODAS) and the Optical Depth in Pressure Space (ODPS) algorithms (Chen et al., 2012), which are based on rigorous line-by-line calculations from the Line-By-Line Radiative Transfer Model (LBLRTM, Clough et al., 1992). For enacting aerosol attenuation effects, the CRTM uses pre-computed look-up tables, which calculate aerosol optical properties, specifically the extinction coefficient, single-scattering albedo, asymmetry factor, and phase function coefficients. The CRTM version 2.2 and 2.3 (Johnson et al., 2021) contain the optical look-up table based on the aerosol types of the mass-based Goddard Chemistry Aerosol Radiation and Transport (GOCART, Chin et al., 2002; Colarco et al., 2010) module, for their radiative effects from the ultraviolet to the infrared. Operationally, given aerosol types, radius, concentration, and ambient relative humidity, CRTM generates aerosol optical profiles that the radiative transfer solver requires for multi-scattering simulations and radiance calculations. The effect of aerosols on MW sensors is not considered yet because the impact of aerosols on MW radiance is usually very small, given that aerosol size is generally much smaller than MW wavelengths (Petty, 2006). There are ongoing and planned CRTM development efforts to incorporate more aerosol optical tables (such as the Community Multiscale Air Quality model, CMAQ). With the expansion of the aerosol schemes, a new releasing and versioning system for optical tables is essential and currently under discussion. This article, however, mainly discusses the GOCART model, which is the default aerosol scheme in the CRTM version 2.

The GOCART model (Chin et al., 2002, 2014), a bulk aerosol scheme, simulates major tropospheric aerosol components, including dust, sea salt, black carbon (BC), organic carbon (OC), and sulfate. It is one of the most widely used aerosol modules in the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem; see Ukhov et al., 2021, and references therein). It is used in the GEOS framework at GMAO for near-real-time aerosol forecasts (Colarco et al., 2010) as well as in MERRA reanalysis (Buchard et al., 2015) and MERRA-2 reanalysis (Randles et al., 2017). It is also implemented in the Global Forecast System (GFS) framework at NCEP for near-real-time global aerosol forecasts (Lu et al., 2016; Wang et al., 2018; Zhang et al., 2021).

When GOCART was selected as the aerosol module within WRF-Chem, it was configured with 14 GOCART aerosol species (Liu et al., 2011): sulfate, hydrophobic and hydrophilic OC and BC, sea salt in four particle size bins (with radii of 0.1–0.5, 0.5–1.5, 1.5–5, and 5–10 µm), and dust particles in five particle size bins (with radii of 0.1–1.0, 1.0–1.8, 1.8–3, 3–6, and 6–10 µm). A default CRTM look-up table has been used for pre-calculated aerosol optical property parameters for the 14 GOCART aerosol species (Liu et al., 2007; Liu and Lu, 2016). We assume that the particles are spherical and externally mixed. We also assume lognormal size distributions for sulfate and carbonaceous aerosols as well as for each sea salt and dust bin. The lognormal size distribution for N particles can be expressed as follows (d'Almeida et al., 1991):

$$\begin{array}{}\text{(1)}& n\left(\ln r\right)={\displaystyle \frac{N}{\sqrt{\mathrm{2}\mathit{\pi }}\ln \left({\mathit{\sigma }}_{\mathrm{g}}\right)}}\exp \left[-{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{\left({\displaystyle \frac{\ln r-\ln r_{\mathrm{g}}}{\ln \left({\mathit{\sigma }}_{\mathrm{g}}\right)}}\right)}^{\mathrm{2}}\right],\end{array}$$

where r is a radius, r_g is the geometric median radius, and σ_g is the geometric mean standard deviation. The kth moment of the distribution can be expressed as follows (Binkowski and Roselle, 2003):

$$\begin{array}{}\text{(2)}& M_k=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{r}^{k}n\left(\ln r\right)\mathrm{d}\ln \left(r\right)=r_{\mathrm{g}}^k\exp \left[{\displaystyle \frac{k^{\mathrm{2}}}{\mathrm{2}}}{\ln }^{\mathrm{2}}\left({\mathit{\sigma }}_{\mathrm{g}}\right)\right],\end{array}$$

where M₀ is the number N of aerosol particles, and M₂ and M₃ are proportional to the total particulate surface area and volume, respectively. Thus, the effective radius (r_eff) can be defined as follows:

$$\begin{array}{}\text{(3)}& r_{\mathrm{eff}}={\displaystyle \frac{M_{\mathrm{3}}}{M_{\mathrm{2}}}}=r_{\mathrm{g}}\exp \left[{\displaystyle \frac{\mathrm{5}}{\mathrm{2}}}{\ln }^{\mathrm{2}}\left({\mathit{\sigma }}_{\mathrm{g}}\right)\right].\end{array}$$

Table 1 lists the GOCART size parameters (particle density, effective radius, and geometric standard deviation) and refractive indices at 550 nm used in CRTM version 2. The optical properties of each aerosol species are computed based on Mie scattering theory. Hydrophilic aerosol particle size increases as relative humidity (RH) of the ambient atmosphere increases. Therefore, the water content in aerosol needs to be considered when calculating the refractive index. The effective radius growth factor for hygroscopic aerosols may be theoretically calculated or obtained from a pre-calculated look-up table (d'Almeida et al., 1991). In this study, the hygroscopic growth factor used for the GOCART model (Chin et al., 2002) is adopted and given in Table 2. Once the growth factor a_g is evaluated, the refractive index n_r for the hygroscopic aerosol can be calculated using a volume mixing method as follows:

$$\begin{array}{}\text{(4)}& n_{\mathrm{r}}=n_{\mathrm{w}}+(n_{\mathrm{o}}-n_{\mathrm{w}})\times a_{\mathrm{g}}^{\mathrm{3}},\end{array}$$

where n_o and n_w are the refractive indices for dry aerosols and water, respectively. We adopt the refractive index n_o from the Optical Properties of Aerosols and Clouds (OPAC) dataset (Hess et al., 1998), while the water refractive index is given by Hale and Querry (1973).

Table 1 Goddard Chemistry Aerosol Radiation and Transport (GOCART) size distribution parameters and refractive indices at 550 nm for dry aerosols.

Table 2 Hygroscopic aerosol growth factor a_g as a function of the ambient relative humidity (RH).

The GOCART model used by GMAO and NCEP for aerosol forecast and reanalysis has evolved to use five sea salt size bins (with radii of 0.03–0.1, 0.1–0.5, 0.5–1.5, 1.5–5, and 5–10 µm). The first sub-micron sea salt bin was added to facilitate optical properties and aerosol–cloud interaction studies (Colarco et al., 2010) but was excluded from the previous GOCART versions as well as the WRF-Chem GOCART model. While GMAO's GEOS and NCEP's GFS contain 15 GOCART aerosol species, the CRTM aerosol module has also not yet been modified to include the new added sub-micron sea salt bin (see Table 1). To overcome this discrepancy, the latest GSI/CRTM release (i.e., GSI 3.7 and CRTM 2.3) combines the mixing ratios from the two sub-micron sea salt bins in order to use the aerosol optical property parameters from the original GOCART model. This limitation is acknowledged in this article and will be addressed in a future CRTM release (see Sect. 4).

While the CRTM is primarily designed for computing satellite radiances, an additional module was added to CRTM by Liu and Lu (2016) to compute aerosol optical depth (AOD). This CRTM-AOD module enables the GSI system to assimilate AOD observations (Liu et al., 2011; Schwartz et al., 2014; Pagowski et al., 2014). This article, however, is focused on the observation operator for radiance, and we refer the reader to Pagowski et al. (2014) for the description of the AOD observation operator and GSI AOD data assimilation.

2.3 Running aerosol-aware GSI analysis

The operational version of GSI maintained by the NOAA/NCEP Environmental Modeling Center (EMC) is utilized in the present study. Its source code and associated static files are distributed through the GitHub repository (https://github.com/NOAA-EMC/GSI, last access: 8 February 2022). An open-access repository of GSI source code is archived on Zenodo (Lueken et al., 2021). To run the GSI analysis, the reader can refer to the user guide for GSI v3.7 (the latest version released as of April 2021), which is available at https://dtcenter.ucar.edu/com-GSI/users/docs/users_guide/html_v3.7/index.html (last access: 8 February 2022​​​​​​​). In addition, an online tutorial is available at https://dtcenter.ucar.edu/com-GSI/users/tutorial/online_tutorial/index_v3.7.php (last access: 8 February 2022​​​​​​​). For CRTM, the user guide and tutorials can be found at https://www.jcsda.org/jcsda-project-community-radiative-transfer-model (last access: 8 February 2022​​​​​​​). Thus, only a brief description of aerosol-affected BT calculations is given.

A regression test “global_C96_fv3aerorad” has been introduced into the NOAA/EMC GSI code repository (pull request no. 32) to assure the functionality of aerosol-aware BT derivations in GSI/CRTM works as expected. This regression test uses a sample background file taken from the aerosol member of the Global Ensemble Forecast System (GEFS-Aerosol; Zhang et al., 2021). All 15 GOCART aerosol species are passed along to the CRTM. In addition to the background file, a user needs to modify the configuration files, anavinfo and satinfo, in the “fix” directory. The anavinfo file is the information file to set control and analysis variables. The satinfo file is the information file to specify satellite channels to be assimilated and associated parameters. For an aerosol-aware experiment where aerosol absorption and scattering are included in BT calculations, aerosol species are specified in the “chem_guess” section of anavinfo, and sensors and channels are set to 1 in the “iaerosol” column of satinfo. The reader can refer to the fv3aerorad_satinfo.txt and anavinfo_fv3aerorad for the aerosol-aware configuration. The corresponding namelist (gsiparm.anl) can be found in the “global_C96_fv3aerorad” section (line 2931–3046) in regression_namelists.sh under the “regression” directory. It should be noted that the namelist variable, “lread_ext_aerosol”, determines how GSI ingests the aerosol information from background files or external files. An open-access repository of fixed files and sample data for the “fv3aerorad” regression test is archived at Zenodo (Lu et al., 2021).

3 Numerical results

3.1 Aerosol impacts on BT calculations

To illustrate how an aerosol transmittance correction is required within satellite radiances assimilated into meteorological data assimilation systems, we present a detailed analysis of a single-cycle GSI experiment (the AER experiment) using GOCART fields from MERRA-2 at 12:00 Z​​​​​​​ on 22 June 2020​​​​​​​. This time is chosen because it captures a strong Saharan dust event that covers the trans-Atlantic region. A baseline GSI experiment (the CTL experiment) with the anavinfo and satinfo resource files reverted back to the default aerosol-blind configuration was also conducted. Both experiments used the same first-guess fields and assimilated identical conventional and satellite observations within a ±3 h assimilation window. In AER, the aerosol transmittance effects were only considered in the CRTM simulation for IR sensors.

Figure 1 Aerosol column mass density (kg m⁻²) from MERRA-2 at 12:00 Z on 22 June 2020: (a) dust, (b) sea salt, (c) carbonaceous, and (d) sulfate.

Figure 1 shows the global aerosol column mass density distribution from MERRA-2 at 12:00 Z on 22 June 2020. The panels a, b, c, and d depict dust, sea salt, carbonaceous, and sulfate aerosols, respectively. Dust plumes spread over northern Africa, the tropical Atlantic Ocean, the Middle East, and northwestern China. Wind-driven sea salt aerosols are seen over tropical and Southern Hemisphere oceans. Carbonaceous and sulfate aerosols mainly appear in areas with extensive biomass burning and fuel combustion activities (note: one order smaller than dust and sea salt). The overall aerosol loading is dominated by mineral dust. Wu et al. (2020) evaluated the dust spatiotemporal variations of MERRA-2 against satellite observations and global model simulations. They found that MERRA-2 agrees well with satellite observations due to the assimilation of satellite AOD. But in North America and the Arctic, the dust burden in MERRA-2 is much larger than those in other models despite having similar dust emissions fluxes. The high dust burden over these regions is due to a higher mass fraction of fine dust and enhanced dust transport. Furthermore, Bullard et al. (2016) reported that large gaps exist in our understanding of basic characteristics of high-latitude dust sources. This highlights the importance of representing aerosol emissions, transport, removal, and size distribution in global models in correctly simulating aerosol spatiotemporal distributions.

Figure 2a shows the first-guess BT differences of IASI on board METOP-A between the two experiments (AER − CTL) in the IR atmospheric window channels over dust-, sea-salt-, carbonaceous-, and sulfate-dominant regions. The stratification criterion for each type is where the fraction of column mass density of the dominant species, from MERRA-2, is larger than 0.65 (shown in Fig. 2b). Figure 2a shows that dust aerosols generate the stronger cooling effects, about 0.7 K at the thermal IR window region (∼10 µm), than other species. The importance of correcting for aerosol transmittance effects within BT algorithms has been reported in previous studies (Sokolik, 2002; Weaver et al., 2003; Pierangelo et al., 2004; Matricardi, 2005; Merchant et al., 2006; Kim et al., 2018; Wei et al., 2021). Table 3 describes the range and the average of total aerosol column mass density over the regions with different dominant aerosol species. It shows that the total loading of aerosols is similar over the dust- and carbonaceous-aerosol-dominated regions. This indicates that the stronger cooling effects by dust aerosol on BT in the IR window region is not due to stronger loading. Note that in the Northern Hemisphere, the high-latitude region is characterized as dust-dominant except for the Russian Far East in MERRA-2 (Fig. 2b). While anomalous or erroneous modeled aerosol loading may bias the results, the finding that dust has the largest impact on the BT simulations, reported in this study and previous studies, remains unchanged. Therefore, we focus our remaining analysis on dust over tropical Africa and the mid-Atlantic.

Figure 2 (a) The differences (AER − CTL) of first-guess brightness temperatures in the IR window region of IASI on board METOP-A. (b) The corresponding regions dominated by different aerosol species at 12:00 Z on 22 June 2020. The data counts for each species are labeled in panel (b).

Table 3 The range of aerosol column mass density (kg m⁻²) from MERRA-2 at the regions dominated by different aerosol species (fraction over 0.65) of IASI on board METOP-A at the cycle of 12:00 Z on 22 June 2020.

Figure 3 displays the AER − CTL difference in the simulated BTs and their respective first-guess departures (observed minus first guess, OMF) calculated at the 10.39 µm channel from IASI on board METOP-A. The figure focuses on North Africa and the trans-Atlantic region, where a large dust plume spans the region. Significant aerosol cooling (∼4 K) in BT was found in the aerosol-aware experiment (Fig. 3a) due to the large plume. Comparing the first guess departures from CTL and AER experiments (Fig. 3b and c) shows that OMFs for AER are warmer than CTL (compare 0.27 K vs. −0.09 K). Note that some observations assimilated in CTL were rejected in AER (near 55^∘ W and 15^∘ N) and vice versa (near 65^∘ W and 15^∘ N, and over Africa). This feature suggests that the quality control has been influenced by including aerosol transmittance effects in CRTM. Over the trans-Atlantic region, the aerosol-aware experiment assimilated several observations with larger first-guess departures located in the strong dust plume (Fig. 3d). Figure 4 presents the scatter plot of dust column mass density versus OMF differences (AER − CTL) for these data points assimilated in AER at 12:00 Z on 22 June 2020. The data points with large OMF differences are corresponding to the areas with higher dust loading. Nevertheless, when considering aerosol information, the root-mean-square first-guess departures decreased 0.08 globally and 0.42 K over the trans-Atlantic region at this channel (not shown here). This implies that simulated BTs in the aerosol-aware run are in better agreement with the observations.

Figure 3 (a) Simulated BT differences (AER − CTL), (b) bias-corrected OMF from the CTL experiment, (c) bias-corrected OMF from the AER experiment, and (d) OMF differences (AER − CTL) for 10.39 µm channel of IASI on board METOP-A. All the data are from the analysis cycle at 12:00 Z on 22 June 2020. Contours of total column mass density from MERRA-2 are plotted in panel (d).

Figure 4 The scatter plot of dust column mass density from MERRA-2 against the first-guess departure differences (AER − CTL) assimilated in AER experiment (without bias correction) at 12:00 Z on 22 June 2020.

Figure 5 shows (a) the global differences in analyzed temperature at 900 hPa between the two experiments and (b) the total aerosol column mass density incorporated into the GSI/CRTM system. When aerosol transmittance effects are considered in the BT calculations, the air temperatures are not only adjusted over aerosol-laden regions but also across the globe. The impact is shown outside aerosol-active regions, which could be attributed to the change from the spatial correlation in the GSI background error covariance. Over the trans-Atlantic region where the dust loading is high (shown in Fig. 1a), the AER experiment produces 0.5 to 1 K of warming relative to CTL. As dust travels off the west coast of Africa into the Atlantic, the particles are lifted and carried by the Saharan air layer (SAL), around 800–600 hPa (Diaz et al., 1976; Karyampudi et al., 1999). In the case of 12:00 Z on 22 June 2020, MERRA-2 captured the dust transport within SAL, and air mass is increasingly composed of fine dust particles due to the gravitational settling of coarser particles (not shown here). Wei et al. (2022) conducted a series of CRTM v2.3 experiments using idealized dust profiles and reported that mass loading and the altitude of the dust layer are the primary and secondary factors affecting the BT simulations, respectively; changes in the fine versus coarse particle partition show little influence on the BT simulations. Based on these results we speculate that elevated dust plume retains unneglected influences on BT calculations (Fig. 3a). Experiments with robust estimated aerosol distributions over extended time periods are needed to quantify the sensitivity of GSI analysis to aerosol-aware CRTM calculations. This paper, however, is intended to provide a joined-up documentation for the CRTM aerosol option, and thus unraveling these questions is beyond the scope of this study.

Figure 5 (a) The differences (AER − CTL) of analyzed temperature (K) at 900 hPa and (b) the corresponding aerosol column mass density (kg m⁻²) from MERRA-2 at 12:00 Z on 22 June 2020.

3.2 Aerosol impacts on the analysis

The experiments reported in this section were produced with the NCEP GFS version 14 and the corresponding GDAS. Our experiments used a coarser resolution, T670 (∼30 km) for the model and T254 (∼80 km) for the analysis, different from the NCEP operational GFSv14 configuration at T1534 (∼13 km) and T574 (∼27 km). The experiments covered the August 2017 period, initialized from NCEP's archived GDAS analysis on 25 July at 00:00 Z. The analysis cycles every 6 h (at 00:00, 06:00, 12:00, and 18:00 Z), with a ±3 h assimilation window and continuous data utilization. The control experiment (CTL_cyc) was an aerosol-blind fully cycled experiment where aerosol effects on radiances are not considered (as they are by default). The aerosol experiment (AER_cyc) was an aerosol-aware fully cycled experiment where aerosol-affected satellite radiances are taken into account. Here, we used CRTM version 2.2.4. Time-varying three-dimensional GOCART aerosols were taken from NCEP's archived NEMS GFS Aerosol Component (NGAC) v2, which simulates the emission, transport, and removal of the GOCART aerosols (Wang et al., 2018).

Figure 6 displays the statistics of analysis departures (observation minus analysis, OMA) from CTL_cyc and AER_cyc to evaluate the performance of temperature analysis at the lower atmosphere over the tropical region (20^∘ S–20^∘ N). The positive value of mean OMAs indicates that both experiments have cold biases in the tropical region. It shows neutral impact on the root mean square (RMS) and slightly positive impact on the cold biases. The latter implies that the departure of temperature analysis becomes larger when considering aerosol transmittance effects during the data assimilation (i.e., AER_cyc).

Figure 6 The comparison of the RMS and mean analysis departures (observation minus analysis, OMA) against in situ measurements (e.g., radiosonde) of temperature with pressure over 1000 hPa at the tropical region (20^∘ S–20^∘ N) during 00:00 Z on 1 August to 18:00 Z on 28 August 2017.

Medium-range forecasts of AER_cyc are examined against CTL_cyc using the verification package from NOAA/NCEP EMC (https://www.emc.ncep.noaa.gov/gmb/STATS_vsdb, last access: 8 February 2022). Figure 7 displays the scorecard of anomaly correlation and root-mean-square error (RMSE) for the day-1, -3, and -5 forecasts over 1–28 August 2017. Anomaly correlation coefficients show neutral to positive impact on day-1 forecasts of wind and temperature fields when aerosol cooling effects in BTs are considered. The RMSE scorecards show the forecast improvements in the wind, temperature, and height fields throughout the troposphere over the tropics (20^∘ S–20^∘ N) and at the upper level over the Northern Hemisphere (20–80^∘ N). For the Southern Hemisphere (20–80^∘ S), however, there is neutral impact or degradation in the forecasts, which is likely due to cloud contamination and mixture of sea salt and aged smoke/sulfate aerosols. Compared to both hemispheres, the tropical forecasts show the most improved statistics in the aerosol-aware analysis, which may be attributed to larger aerosol loading in this region. While the RMSE scorecard focuses on background (i.e., time-averaged) fields, it should be noted that evaluation of the aerosol impacts on the analysis and forecasts of African easterly wave that developed Hurricane Harvey and Gert in 2017 is presented in Grogan et al. (2021).

Figure 7 Scorecard of anomaly correlation and RMSE of comparison between AER_cyc and CTL_cyc. Green colors mean AER_cyc is better than CTL_cyc at 95 % (filled box), 99 % (), and 99.9 % () significance level. Red colors mean AER_cyc is worse than CTL_cyc at 95 % (filled box), 99 % (), and 99.9 % () significance level. Grey boxes mean no statistically significant difference between AER_cyc and CTL_cyc. Blue boxes are not statistically relevant. The statistics are calculated between 20 to 80^∘ of latitude for both hemispheres. The data between 20^∘ S and 20^∘ N are used for the tropical region.

4 Conclusions and future outlook

This article described aerosol absorption and scattering calculations of the CRTM version 2 in the GSI analysis. We also conducted sensitivity experiments to investigate the aerosol-affected GSI analysis in both single-cycle and fully cycled runs. Both GSI and CRTM are well documented with user guides, tutorials, and code repositories available online. This article is primarily a joined-up documentation for aerosol absorption and scattering calculations in the CRTM version 2 and GSI. It also provides guidance for prospective users of the CRTM aerosol option. Scientific aspects of aerosol-affected BT in atmospheric data assimilation are briefly discussed. Specifically, numerical experiments were conducted to illustrate how including aerosol radiative effects in CRTM changes the GSI analysis. We found that taking the aerosols into account reduces simulated BT in thermal window channels over dust-dominant regions. Assimilating aerosol-affected BTs produces a warmer analyzed lower atmosphere. From the verification scorecard, neutral to positive results are found in the fully cycled, aerosol-aware experiment.

The CRTM team, in coordination with its partners and collaborators, is building a robust capability to accurately and consistently simulate the emission, absorption, and scattering properties of all (radiatively important) atmospheric constituents. There are several ongoing and planned efforts to enhance the CRTM aerosol module. For example, more aerosol optical look-up tables have been added and the calculations of aerosol optical properties are being evaluated. In addition, the CRTM is being refactored toward a more flexible aerosol interface to handle aerosol optical look-up tables as well as to support aerosol specifications from other operational aerosol models, such as CMAQ. Other aerosol-related efforts include but are not limited to improving the physical representation of aerosols and including active sensors such as aerosol lidar. These developments, once implemented and tested, will be reported in future paper.