The Impact of Aerosols on Satellite Radiance Data Assimilation Using NCEP Global Data Assimilation System

: Aerosol radiative effects have been studied extensively by climate and weather research communities. However, aerosol impacts on radiance in the context of data assimilation (DA) have received little research attention. In this study, we investigated the aerosol impacts on the assimilation of satellite radiances by incorporating time-varying three-dimensional aerosol distributions into the radiance observation operator. A series of DA experiments was conducted for August 2017. We assessed the aerosol impacts on the simulated brightness temperatures (BTs), bias correction and quality control (QC) algorithms for the assimilated infrared sensors, and analyzed temperature ﬁelds. We found that taking the aerosols into account reduces simulated BT in thermal window channels (8 to 13 µm) by up to 4 K over dust-dominant regions. The cooler simulated BTs result in more positive ﬁrst-guess departures, produce more negative biases, and alter the QC checks about 20%/40% of total/assimilated observations at the wavelength of 10.39 µm. As a result, assimilating aerosol-affected BTs produces a warmer analyzed lower atmosphere and sea surface temperature which have better agreement with measurements over the trans-Atlantic region.


Introduction
Aerosols affect Earth's energy balance through the absorption and scattering of solar and thermal radiation [1]. Aerosols also affect the climate through their indirect effects on cloud microphysics, reflectance, and precipitation [2][3][4][5][6][7]. In recent decades, the impacts of aerosol direct and indirect radiative effects on weather forecasts have been studied extensively [8][9][10][11][12][13]. Moreover, aerosols can influence the satellite radiance observations assimilated in the form of brightness temperature (BT) in data assimilation (DA) systems.
Many studies have demonstrated pronounced reduction of BT at the thermal infrared (IR) window region (8 to 13 µm) due to the presence of aerosols. During the Saharan Dust Experiment (SHADE), 2 to 4 K cooling effects of BTs in the IR window region by the nadir view measurements of Airborne Research Interferometer Evaluation System (ARIES) was reported [14]. In the Saharan Mineral Dust Experiment-2 (SAMUM-2), 1 to 4 K cooling by dust and 1 K cooling by smoke in the IR window region were reported [15]. The cooling effects are also reported by studies that incorporate climatological or artificial aerosol profiles into BT simulations [16][17][18][19][20][21][22]. Moreover, the quality of IR observation retrievals can be improved by considering aerosol radiative effects. Merchant et al. [19] indicated that the bias of the sea surface temperature (SST) retrieval of Spinning Enhanced Visible and Infrared Imager (SEVIRI) can be reduced with proper estimation of dust aerosol impacts on IR observations. Divakarla et al. [23] reported that the temperature retrievals of Infrared Atmospheric Sounding Interferometer (IASI) have less deviation from reanalysis when the dusty observations were removed.
Despite the aerosol impacts on radiance observations, the incorporation of aerosols into DA systems has received little research attention. Letertre-Danczak [24] demonstrated that rejecting IR observations affected by aerosols can improve the temperature analysis in the lower atmosphere. However, this approach prohibits the use of aerosol-affected IR measurements and thus the information content of the spectra is not fully exploited. Geer et al. [25] reported that about 60% high-spectral IR observations over aerosol-laden regions are removed, such as the trans-Atlantic region and Sahara Desert. Instead of rejecting aerosol-affected observations, Weaver et al. [26] and Kim et al. [27] demonstrated that assimilating aerosol-affected observations introduces a warmer analyzed temperature in the lower atmosphere.
The aerosol effects on the Earth's system and radiance observations have been wellrecognized, but aerosol information is not considered in most DA systems. This study aims to incorporate aerosol information into the radiance simulations to investigate its impacts on meteorological analyses. The present study is outlined as follows. The DA system and experimental designs are described in Section 2. The sensitivity of the DA system to aerosol information is presented in Section 3. Discussion and conclusions are presented in Section 4.

NCEP Global Data Assimilation System
In this study, we use version 14 of the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS v14). The forecast model is National Oceanic and Atmospheric Administration (NOAA) Environmental Modeling System (NEMS) spectral dynamical core with GFS physics. The analysis system, Global Data Assimilation System (GDAS), is a hybrid four-dimensional ensemble-variational (4DEnVar) [28,29] analysis system based on Gridpoint Statistical Interpolation (GSI). NCEP GDAS generates meteorological analysis every 6 h (i.e., 00Z, 06Z, 12Z, and 18Z) with the following information: 1.
Prescribed static background error covariance and observation error covariance The First Guess at Appropriate Time (FGAT) [30] method is applied to derive the first-guess fields at the corresponding time of observations from −3, 0, +3 h of the actual analysis time.
To assimilate radiance observations in GSI, the Community Radiative Transfer Model (CRTM) [31,32] is used as the radiance observation operator. It simulates the BT at the top of the atmosphere and derives the Jacobian matrix (i.e., the first derivative) of BT with respect to model state variables. More information of the CRTM is described in Section 2.2.
In a variational DA system, the magnitude and shape of analysis increments are determined by the background errors, observation errors, and the Jacobian for the corresponding state variables. In addition to the impacts on analysis increments, the Jacobian for surface emissivity and temperature derived by CRTM are used in the variational bias correction (VarBC) scheme [33,34] and quality control (QC) in GSI. VarBC estimates the biases within the first-guess departures as parts of control variables during the minimization process in the GSI. Four types of BC predictors are actively applied to IR sensors: (1) global offset, (2) temperature lapse rate and square of temperature lapse rate, (3) surface emissivity sensitivity, and (4) fourth-order polynomial formula of scan angle [34].
Regarding the QC procedures for IR observations, a clear-sky approach including the checks for clouds, sensitivity of surface temperature, and gross error, etc. is applied [35].
The QC check of clouds is based on the bias-corrected first-guess departures because cloud fields in first-guess were not ingested in CRTM simulations for IR sensors. Another QC check is based on the increment of SST derived through the physical retrieval algorithm [36]. When the increment of SST is larger than a channel dependent threshold, the observations at the corresponding channels are rejected. As shown here, aerosol information is not explicitly considered in the QC procedures for IR observations.

CRTM
The CRTM was developed at the Joint Center for Satellite Data Assimilation (JCSDA) with contributions from NOAA scientists. It uses the advanced double-adding (ADA) method [37] for fast and accurate calculations of multiple-scattering in radiative transfer. The adjoint method is applied in the CRTM to obtain the Jacobians. Ding et al. [38] reported that the CRTM provides accurate BT simulations under both clear-sky and cloudy-sky conditions against the simulations from the line-by-line radiative transfer model. Li et al. [39] showed that the CRTM Jacobian calculation is in good agreement with the finite differential method under cloudy conditions.
The CRTM considers the effects of aerosols from ultraviolet to infrared regions, but not for the microwave region since the wavelength is much larger than the particle size of aerosols. Spherical particles were assumed with the lognormal distribution for aerosols. The aerosol refractive indices are adopted from the Optical Properties of Aerosols and Clouds (OPAC) [40]. More information of the CRTM aerosol module is described in Liu and Lu [41]. Although CRTM has the capability to consider the aerosol scattering effect when solving the radiative transfer equation, aerosol information is not ingested into GSI for aerosol-aware BT simulations.

NGAC v2
In this study, we incorporated the 3-hourly three-dimensional aerosol distributions from NEMS GFS Aerosol Component (NGAC) into the BT simulation of CRTM. NGAC is a global aerosol forecast system operated by NCEP since September 2012 [42,43]; it has been replaced by Global Ensemble Forecast System-Aerosols (GEFS-Aerosols) [44] in September 2020. NGAC coupled NEMS GFS with the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model [45,46]. NGAC v2 provided multi-species forecasts out to 120 h at T126 (∼1 • ) resolution on 00Z and 12Z every day. The aerosol species include the 5 bins for both dust and sea salt, sulfate aerosols, and hydrophobic and hydrophilic carbonaceous aerosols (black carbon and organic carbon). The NCEP archived NGAC v2 data is retrieved to provide aerosol information. Details about the verification and evaluation of NGAC v2 can be found in Bhattacharjee et al. [47]. Figure 1 sketches the procedures of assimilating radiance observations for the two GDAS experiments during August 2017, initialized from NCEP's archived GDAS analysis. Due to resource constraints, GDAS experiments are conducted at the lower resolution of T670 (∼30 km) for GFS forecasts and T254 (∼80 km) for ensemble members and analysis. The first experiment, CTL, is an aerosol-blind cycling experiment. In the GSI, CRTM simulates the first-guess BTs based on the short-range forecasts from the previous cycle. Then, the first-guess departures are obtained by subtracting the first-guess BTs from the observed BTs of satellite sensors. The first-guess departures are also called observation minus first-guess (OMFs). After the BC and QC processes, the OMFs involved in the minimization of cost function are determined. The atmospheric and surface analyses are obtained after the minimization is completed and thus provided the initial conditions for short-range GFS forecasts as the first-guesses to the next cycle. The second experiment, AER, is an aerosol-aware offline analysis experiment. It has the identical procedures and the same first-guesses and observations as CTL, but it incorporates the NGAC v2 aerosols and only performs the analysis step. Like other meteorological fields in first-guesses, NGAC v2 aerosol fields are interpolated by the FGAT and ingested into the BT simulation of CRTM in the AER experiment. Under this offline non-cycled configuration, AER analyses do not feed back to the cycling GDAS. Hence, the sensitivity of GSI analysis to aerosol-aware CRTM calculations (i.e., CTL versus AER in Section 3) can be clearly demonstrated. Figure 2 shows the 550 nm aerosol optical depth (AOD) distributions for different aerosol species in NGAC v2, averaged over August 2017. It shows the trans-Atlantic transport of Saharan dust aerosols as well as the major dust source regions, e.g., the Sahara Desert, the Middle East, and northwestern China. Carbonaceous aerosols are seen over the areas with active biomass burning, e.g., Canada, Russia, South America, and tropical Africa. Sea salt aerosols are mainly distributed over the tropical region of the Northwest Pacific Ocean and the Indian Ocean, and the Antarctic Ocean. China is the major source region of the sulfate aerosols. Overall, dust and carbonaceous aerosols are the dominant species during the study period.

Observational Dataset
Both experiments used the non-restricted copy of the observational dataset in NCEP operational system. During the study period, the level 1 product of 29 satellite sensors (15 infrared and 14 microwave sensors) are ingested into GDAS. All satellite data is retrieved and pre-processed by NOAA National Environmental Satellite, Data, and Information Service (NESDIS). Table 1 lists the satellite radiance observations of IR sensors ingested into the GDAS framework. Please note that not all channels of each sensor are assimilated in GDAS. For instance, 117 channels of the 281 channel subset of Atmospheric Infrared Sounder (AIRS) on Aqua are assimilated. Both microwave and IR radiance observations are thinned to 145 km grid points in the GSI. Other in situ measurements (e.g., rawinsonde, synop, buoy, ship, etc.) were retrieved from Global Telecommunications System (GTS) and processed by PREPBUFR (https://www.emc.ncep.noaa.gov/mmb/data_processing/ prepbufr.doc/document.htm, accessed on 10 February 2021) at NCEP. The observational dataset used in GDAS is publicly available at the National Climatic Data Center (NCDC) website (https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/globaldata-assimilation-system-gdas, accessed on 10 February 2021). It should be noted that the dataset retrieved through NCDC website could be different to the dataset used in this study due to the data policy. Most discussions in the results section (Section 3) focus on the aerosol impacts on observations of IASI, which are representative of other IR sensors. Also, assimilating IASI observations provides significant deduction of forecast errors in DA systems [48][49][50][51]. Briefly, IASI is a cross-track IR scanner that covers the spectral domain from 645 to 2760 cm . As a result, the observations from the orbits used in GDAS cycles at 00Z and 12Z cover the trans-Atlantic region which is primarily affected by Saharan dust in our experiments. More details of IASI and its application can be found in Hilton et al. [52] and the references therein.

Statistical Analysis
To quantify the aerosol impacts on the GSI analysis, the root-mean-square (RMS) and mean are applied to first-guess departures (OMFs) and analysis departures (OMAs). The associated information of IR observations in both experiments from 1-28 August 2017 is regridded to 2.5 • by 2.5 • for further diagnosis. Please note that the information for IR diagnosis includes (1) OMFs, (2) OMAs, and (3) bias estimates from BC. The mean of this information from the data points within a gridbox is derived. Additionally, the significance of the difference between experiments is calculated by the T-test method in each gridbox. Besides the RMS and mean, the linear regression and determination coefficients are derived to evaluate analyses against in situ measurements. Figure 3 displays the first-guess BT differences between CTL and AER (AER minus CTL), averaged from 1-28 August 2017, over aerosol-affected regions (total AOD at 550 nm > 0.3) and dominated by different aerosol species. It should be noted that IASI on MetOp-A is overlapped by IASI MetOp-B due to similar values. The dominant species is determined by the fraction of species AOD at 550 nm over 0.65. For observations dominated by dust aerosols, both high-spectral IR sensors ( Figure 3) and low-spectral IR sensors (not shown) show aerosol cooling effects, which is about 1 to 2 K near 4 µm and 3 to 4 K near 10 µm. These two spectral regions with aerosol cooling effects are located at thermal infrared windows, 3.6 to 4.0 µm and 8 to 14 µm, respectively. Dust aerosols (Figure 3a) generate stronger cooling effects on BT simulations than sea salt, carbonaceous and sulfate aerosols (Figure 3b-d). The cooling effects of BT in the present study have shown similar responses at the atmospheric window channels as previous studies based on aerosol profiles specified from models [17,18,22,27] and the field campaigns [14,15], especially the 8.33-12.5 µm region. For dust aerosols, the channel at 10.39 µm produces the strongest cooling effects among assimilated channels of IASI (Figure 3a), but the effects are spatially dependent.   Figure 4d). The 12Z 5 August 2017, was selected due to a strong dust plume in NGAC v2 which is about 0.6 to 1.5 of AOD at 550 nm over the East Atlantic. It should be noted that for both observations and experiments in Figure 4a-c, the small land-sea difference of BTs near the Nigeria region is attributed to the presence of clouds. It also implies that GFS forecasts in CTL capture the cloud signals and thus reduce the land-sea differences of surface temperature at this region. For the aerosol-blind CTL experiment (Figure 4b), simulated BTs are warmer than observed BTs near Africa and the west part of South America (Venezuela and Columbia). When aerosols are considered in the BT derivation (i.e., the AER experiment), simulated BTs become cooler than observations in the trans-Atlantic region. The simulated BTs in CTL qualitatively agree with observations better than AER over the mid-Atlantic (around 40

Brightness Temperature
• W and 10 • to 20 • N). However, the root-mean-square (RMS) deviation between observations and simulated BTs is 8.96 K and 9.25 K for AER and CTL in this domain, respectively. While the RMSE of AOD at 550 nm from NGAC v2 forecasts is about 0.1 over this region [46], it leads to about 0.4 K uncertainty of BT simulations at this channel (not shown) which does not significantly affect the result presented here. Differences in simulated BTs (Figure 4d) show that the largest differences occur within the region enclosed by the higher AOD at 550 nm contours (brown lines), which spans across North Africa and the tropical Atlantic. This indicates that the differences are due to the cooling effect associated with higher aerosol loading. It should be noted that the magnitude of BT differences are larger over North Africa than the Gulf of Guinea, which is due to the different aerosol types. Over North Africa, dust aerosols are the dominant species, while carbonaceous aerosols dominate near the Gulf of Guinea.

First-Guess Departure
The BT cooling effect shown in Figures 3 and 4 subsequently affects the first-guess departures (OMFs) of BT which determine the analysis increments. Figure 5 illustrates the gridded monthly OMFs before QC and BC of IASI on MetOp-A at 10.39 µm for (a) CTL, (b) AER, and (c) the differences (i.e., AER minus CTL). Both CTL (Figure 5a) and AER (Figure 5b) have negative OMFs over most regions. This feature is mainly attributed to the cloud information not incorporated into the BT simulation under the clear-sky framework. Due to the cooler BTs simulated in AER, the OMF differences are globally positive in Figure 5c. It is statistically significant over the dust-impacted areas, such as the trans-Atlantic region, Sahara Desert, and Middle East. This indicates that the incorporation of aerosols in BT simulations shifts the OMFs toward more positive. As a result, AER will have larger positive analysis increments than CTL. Stippling in (c) indicates the difference is 99% significant (p ≤ 0.01) Figure 6 shows monthly mean OMFs (upper panel) and RMS OMFs (lower panel) before QC and BC from AIRS, Cross-track Infrared Sounder (CrIS), and IASI. Across most channels, all four high-spectral IR sensors have less negative mean OMFs in AER than CTL which is the same response as Figure 5 shown above. The smaller RMS OMFs indicate that considering aerosol radiance effect (i.e., AER) can improve the OMFs statistics through the entire thermal IR spectrum, especially the atmospheric window regions. Due to the identical first-guesses in CTL and AER, the standard deviation of mean and RMS OMFs along analysis cycles does not change significantly (not shown). It should be noted that cloud-affected observations were included in the statistics. Since clouds have stronger attenuation than aerosols, the impact of considering aerosol information in BT simulations on the statistics of OMFs is limited. The gap near 11.1 µm (900 cm −1 ) of AIRS on Aqua in Figure 6a is due to no available observations after screening out the unreasonable observed BT (around zero degree) at that particular channel.

Bias Correction and Quality Control
The BC algorithm in the GSI estimates and corrects the biases in the OMFs of radiance observations. Figure 7 illustrates the bias estimates gridded at 2.  Figure 7a displays that the CTL has globally positive values which means OMFs will decrease after the bias correction. In contrast, negative biases are estimated in AER (Figure 7b) over the Sahara Desert, South Asia, and the west Pacific, where aerosol loadings are higher (as shown in Figure 2). Together with OMFs statistics shown in Figure 5c, the OMFs of AER will become more positive than CTL over those regions. It should be noted that the land-sea difference is less significant in AER. Figure 7c indicates that the differences mainly distribute over land and the high latitude regions in both hemispheres.   Figure 8 shows the differences (i.e., AER minus CTL) of bias estimates contributed by temperature lapse rate and surface emissivity sensitivity predictors. Here only these two predictors are displayed because the change of bias estimates from other predictors are small. Figure 8a displays the result of lapse rate term which is estimated by the column integral of the product of layer differences of transmittance and temperature (∑ ∆τ × ∆T). It shows that AER has smaller bias estimates over Africa and the trans-Atlantic regions. Since the first-guesses are identical in CTL and AER, the change of bias estimates of the lapse rate term is subjected to the change of the transmittance profile from CRTM due to the presence of aerosols. Figure 8b displays the results of surface emissivity sensitivity term which is estimated by the BT Jacobian for surface emissivity from CRTM. Please note that this term is only applied to the observations over the non-water surface, including the ice and snow surface at polar regions, to account for the uncertainties in surface types used in CRTM simulations. The significant reduction of this term in AER indicates that the BT Jacobian for surface emissivity is changed considerably after including aerosols into CRTM simulations. It also explains the less land-sea difference in AER (Figure 7b,c). In contrast to the lapse rate term, further study is needed to address how aerosols affect the calculation of BT Jacobian of surface emissivity in the CRTM.  Figure 7, but for bias estimates differences (AER minus CTL) from the predictors of (a) lapse rate term and (b) emissivity term. Stippling indicates the difference is 99% significant (p ≤ 0.01). In the GSI, the QC procedure (described in Section 2.1) checks the bias-corrected OMFs and determines the dataset involved in the minimization process. Table 2 lists the number of observations passed or rejected due to various QC checks for the CTL (top row) and AER (bottom row) at 10.39 µm for IASI on MetOp-A through 1-28 August 2017. Here we list the results of four QC checks used for this IASI channel. Only the observations that passed QC checks (the first column) were assimilated in the GSI. For this channel on IASI, the AER run assimilated fewer observations (cf. 235,997 vs. 262,758). It rejected more observations due to surface emissivity ( s f c ) and temperature (T s f c ) sensitivity (cf. 259,595 vs. 166,291), but rejected less observations due to presence of cloud layers in model (cf. 572,345 vs. 638,676) than the CTL run. The larger OMF in AER also resulted in a small number of rejected observations due to the gross error check. Table 3 lists the changes in the QC results at 10.39 µm for IASI on MetOp-A in AER with respect to CTL. Considering aerosol information could change the QC results by about 21% of observations. For the assimilated observations in CTL, around 39% observations' QC statistics were changed in AER. About 13% and 25% observations rejected in CTL due to clouds and s f c and T s f c sensitivity were changed, respectively. Moreover, most rejected observations in CTL due to physical temperature retrievals were changed in AER (89.8% in Table 3). These changes indicated that QC checks should be carefully addressed for assimilating aerosol-affected observations. Figure 9 displays the horizontal distribution of OMFs rejected by the QC check of surface emissivity and temperature for the same cycle, channel, sensor, and region as Figure 4. Compared to CTL (Figure 9a), AER (Figure 9b) rejected more observations due to this QC check over mid-Atlantic, which contributes to the changes of QC statistics shown in the fifth column of Table 2. As documented in Liang and Weng [35], the QC check of surface emissivity and temperature for IR sensors is determined by a "skin temperature sensitivity related parameter". The observation at a particular channel is rejected when the parameter is larger than 5% of channel observation error. Briefly, the parameter is proportional to the bias-corrected OMF which is larger in AER due to the cooling effect of BT simulations and the negative bias estimates over this region. Therefore, given the same observation error prescribed in AER and CTL, GSI tends to reject the observations with larger OMF (i.e., AER in this case).  The feature of more rejected observations over water at 10.39 µm in AER is shown in other sensors as well. Table 4 lists the count of the assimilated observations over different surface types (water, land, and mixed) for IASI, AIRS, and CrIS in CTL and AER near 10.4 µm during the period. IASI onboard MetOp-B and AIRS have a similar response as IASI onboard MetOp-A, less observations over water and more over land were assimilated in AER. However, the response is different for CrIS, where AER assimilated more observations over water than CTL at this channel. CrIS observations over land are not assimilated because the surface temperature sensitivity at this channel is higher than 0.2 (K/K) which is rejected in the QC procedure of CrIS. Overall, the total number of the assimilated data counts at this channel is reduced in AER. The Gaussianity of probability density function (PDF) of OMFs is a common metric to qualitatively assess whether the biases in the assimilated dataset were removed properly through the BC algorithm. Figure 10 shows the PDFs between CTL and AER before BC (histograms) and after BC (dashed lines) near 10.39 µm for IASI, CrIS, and AIRS during 1-28 August 2017; all PDFs are computed based on the OMFs after QC. The histograms reflect that AER has more positive OMFs in all the four sensors. Comparing the BC influences on CTL and AER, the corrected PDFs of CTL peak at zero line which indicates that the BC algorithm properly removes the biases in OMFs. On the other hand, the pattern of the corrected PDFs skews toward positive slightly in AER. This indicates that BC coefficients from CTL may be unsuitable to estimate biases for aerosol-affected OMFs. Among the four sensors in AER, the OMFs of CrIS on Suomi-NPP (Figure 10b) are better corrected. This feature could be attributed to the rejection of the observations over land for CrIS by QC. It suggests that the BC algorithm may need some adjustment for using aerosol-affected observations over land.

Use of Observation
While the previous section focused on the response at a single channel, we would like to understand how the aerosol impacts on data usage vary with the assimilated channels and with the analysis cycles. For the response at different channels, Figure 11a shows the distribution of the data count differences (AER minus CTL) of assimilated observations from the 164 channels of IASI at each analysis cycle (i.e., every 6 h) during August 2017. When a box locates at the positive range in Figure 11a, it means that AER assimilated more observations from that channel in most analysis cycles. In the wavelengths from 9 to 13 µm, AER assimilates fewer observations than CTL in most analysis cycles. However, for the channels with longer wavelength, e.g., 13.5 to 14.5 µm, AER assimilates more data than CTL. As a result of compensation in data count difference between 9 to 13 µm and 13.5 to 14.5 µm, AER assimilates slightly more observations for IASI onboard MetOp-A. This is shown for each analysis cycle in Figure 11b, which displays the distribution of data count differences from all assimilated channels. When a box locates at the positive range in Figure 11b, it means AER assimilated more observations than CTL for the majority of assimilated channels at that analysis cycle. The mean of assimilated channels (red line) in Figure 11b indicates that the average number of assimilated observations in AER is slightly more than CTL at each analysis cycle. Figure 11. Boxplots for assimilated observation number difference (AER minus CTL) from IASI on MetOp-A in (a) the distribution of analysis cycles along assimilated channels and (b) the distribution of assimilated channels along analysis cycles. The 5th, 25th, 50th, 75th, and 95th percentile are displayed. Red line in (b) is the average value over all the assimilated channels. Table 5 lists the average number of assimilated observations for the IR sensors. It also lists the average of differences (AER minus CTL), and their corresponding temporal (STD T ) and spectral (STD S ) standard deviations of differences. Larger value of STD T represents that the averaged difference of data count over assimilated channels varies considerably at each analysis cycle. Larger value of STD S represents that the averaged difference of data count over each analysis cycle varies considerably among each assimilated channel. For monitoring sensors, all available channels were included in the statistics. Table 5 shows that AER assimilates more observations than CTL from most IR sensors, except the Advanced Very High Resolution Radiometer (AVHRR) on MetOp-A and NOAA-18 and SEVIRI on Meteosat-10. The STD T differences among sensors also reflect the orbits of satellites. The sensors on polar orbiting satellites have larger STD T than geostationary satellites due to the heterogeneous distribution of aerosols. The larger STD S in IASI, CrIS, and AIRS could be attributed to the large set of assimilated channels which include aerosol sensitive and insensitive channels. Overall, the statistics indicate that more information contents from IR sensors were exploited in the AER experiment. This feature is primarily attributed to the channels of wavelength from 13.5 to 14.5 µm (Figure 11a) which are not measured in AVHRR and not the assimilated channel of SEVIRI in our experiments. Table 5. Average assimilated observations number among infrared sensors over 1- 28 August 2017. Rows below the dash line are monitored-only sensors. STD T and STD S stand for the temporal and spectral standard deviation of difference between AER and CTL, respectively. The relative changes of each sensor in AER are indicated in the parenthese. 3.5. Temperature Analysis Figure 12 presents the mean analysis difference between CTL and AER of (a) sea surface temperature and (b) temperature at 850 hPa. Also shown is the vertical cross-section of mean temperature differences (10 • E-80 • W) meridionally averaged over the northern (Figure 12c). The differences indicate that analyzed SSTs in AER are about 0.3 to 0.5 K warmer than in CTL in the Atlantic Ocean (Figure 12a). This warming effect also extends to 850 hPa (Figure 12b) over Africa to the tropical Atlantic Ocean. However, there is a region with cooler 850 hPa temperatures near the Africa coast. Figure 12c indicates a warmer low level of atmosphere (∼0.15 K) crossing Africa to the Atlantic Ocean and a cooler mid-level over the Atlantic Ocean near Africa coast. This change of vertical temperature structure may induce more instability over this region. Consequently, the feature could affect the evolution of Africa easterly wave and the subsequent tropical cyclone activity. The OMAs from in situ datasets (e.g., radiosonde, buoy, etc.) can be used to evaluate the analysis performance because the measurements are more accurate than remote sensing data. Figure 13 displays the scatterplots of SSTs between measurements and analysis over four regions with larger SST differences during August 2017. As shown in Figure 12a, the West Atlantic has the strongest warming effects (∼0.5 K); the East Atlantic region has the warming effects from 0.2-0.4 K; Arabian Sea region ranges around 0.1 to 0.2 K warming; West Pacific is about 0.2 K cooling in analyzed SSTs. The RMS and mean OMAs are labeled in the title of each panel.

Sensors
The slope of four regression lines in each region implies that the analyzed SSTs in AER are a slightly better fit to the measurements except the West Atlantic. The coefficients of determination, R 2 , are similar over the four regions, which means that the analyzed SSTs are not significantly improved by considering aerosols in BT simulations. While there are no significant differences over these regions, the smaller RMS OMAs indicate that the SST analyses in AER are in better agreement with measurements. Also, the smaller absolute mean OMAs implies the AER experiment has reduced the warm biases in the West Pacific and the cold biases in other three regions. It should be noted that, compared to the large sampling size in the West Atlantic (>30,000 measurements), the other three regions have less than 100 measurements which could be less statistically meaningful, especially over the West Pacific (only 24 points during the entire period). O'Carroll et al. [53] reported that IASI SSTs have stronger cold biases than AVHRR SSTs against measurements from drifting buoys. Compare the SST retrieval methods from both, the aerosol corrections (e.g., the Saharan Dust Index from SEVIRI [19]) are applied in the retrieval method of AVHRR [54] but not in the method of IASI [55]. Similarly, cold biases of SST analyses in the CTL experiment can be reduced by considering aerosol information in BT simulations (i.e., AER).

Discussion and Conclusions
This study investigates the aerosol impacts on meteorological analysis by incorporating aerosol information into the radiance observation operator. Two GDAS experiments were conducted for 1-28 August 2017. These include an aerosol-blind cycling analysis experiment (CTL) and an aerosol-aware offline analysis experiment (AER). Overall, the aerosol-aware BT simulations lead to the following features in the AER run: (1) lower simulated BTs at the top of the atmosphere, (2) more positive OMFs, (3) more negative estimated biases, (4) considerable changes in the outcome of QC checks, (5) more IR observations assimilated, (6) warmer analyzed SST over the Atlantic Ocean and lower atmosphere over Africa to the central Atlantic Ocean, and (7) a slightly better agreement with SST measurements. In addition to the experiments presented in this study, a companion aerosol-aware cycling experiment was conducted. The results can be found at Grogan and Lu [56] and Wei et al. [57].
In terms of the aerosol information, there can be significant seasonal and annual variation of the loading, species, and vertical distribution in regions around the world. Also, the forecast errors of the aerosol model contribute uncertainties to the BT simulations. Our results suggest that cases with lower dust loading or dominated by a non-dust species would produce smaller cooling effects on the BT simulations. This would lead to smaller warming effects on the analyzed temperatures in the lower atmosphere. Future research should be conducted to comprehensively evaluate the effects of aerosols on the analysis system under a wider range of aerosol loading.
As the radiance observation operator in the GSI, the CRTM is the key component for direct assimilation of radiance observations. In this study, aerosol induced changes are driven by the aerosol-aware CRTM simulations, which uses the spherical assumption for aerosols. It has been recognized that non-spherical aerosols can affect the accuracy of radiance calculations in the IR window [58,59]. Hence, the spherical assumption could introduce uncertainties in our results. This issue has been addressed by the latest public release of the CRTM v2.4 (as of October 2020), which incorporates the non-spherical optical properties into the CRTM aerosol module. Nonetheless, the transmittance and Jacobians from the CRTM are incorporated into the BC, QC, and minimization processes for assimilating radiance observations in the GSI. It is still unclear how aerosols affect the CRTM simulation, the Jacobians in particular. Therefore, further investigation is needed to understand the response of the CRTM simulation to the presence of aerosols.
Our results showed that considering aerosol information in the CRTM reduced the discrepancies between the simulated and observed BTs. This in turn produced considerable changes in the outcome of QC checks and a more positive skewness in the PDFs of biascorrected OMFs in the IR window region. Therefore, future studies are needed to explore the aerosol-aware QC and the aerosol-based BC predictor for assimilating aerosol-affected IR observations in the system. For aerosol-aware QC, an observation-based aerosol detection method can be used to identify the aerosol-affected IR observations. An aerosol-based BC predictor could be incorporated into the bias estimation to account for the uncertainties of the modeling aerosols in OMFs. These capabilities can constrain the aerosol effects on IR observations, and thus the aerosol-affected IR observations can be assimilated in the analysis system such as the all-sky DA approach for cloudy observations [60].

Funding:
The present work is funded by NOAA NWS NGGPS program (award number NA15NWS4680008).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data that supports the findings of this study is archived at NESDISfunded Supercomputer for Satellite Simulations and Data Assimilation Studies (S4) cluster and can be made available to the readers upon request. The observations used in the GDAS experiment are available in the NCDC website (see text).