1. Introduction
Aerosol exists widely in the atmosphere and is one of the important factors affecting air quality and climate change [
1,
2]. Aerosol Optical Depth (AOD) is one of the important parameters to describe aerosol optical properties, and accurate AOD information can be used to assess global climate change and atmospheric pollution [
3].
At present, aerosol optical characteristics monitoring research mainly includes ground observation and satellite remote sensing observation. Ground-based observations, such as the NASA-led Aerosol Robotic Network (AERONET), provide high-precision measurements of aerosol optical properties free of charge to researchers around the world [
4]. However, the available data from ground observation stations in China are still limited, and ground-based observations cannot provide detailed monitoring results of aerosol optical properties in a large area. The development of aerosol satellite remote sensing observation has a history of over 40 years [
5,
6]. The first satellite sensor to provide aerosol optical depth products is AVHRR (Advanced Very High-Resolution Radiometer), which is a scanning radiometer with five bands and has been mounted on NOAA (National Ocean and Atmosphere Administration) since 1978 [
2,
7,
8,
9]. A large number of aerosol retrieval algorithms have emerged [
1,
6,
10,
11] after years of development of aerosol satellite remote sensing.
So far, the most widely used algorithms are Dark Target (DT) and DeepBlue (DB) algorithms officially adopted by MODIS [
12,
13,
14]. Aerosol inversion is based on multiangle sensors, represented by AATSR [
15,
16] and MISR sensors [
17]. AATSR has dual-angle observation capabilities. When establishing the relationship between the reflectance of different observation angles, it is usually assumed that the surface bidirectional reflection characteristics are negligible compared with the scattering characteristics of the surface and the atmosphere [
10,
18,
19,
20]. MISR is a true multiangle observation sensor the inversion algorithm was assumed to have the same BRF shape function as different bands, depending on the sun observation geometry [
21]. Xue et al. [
22] proposed the Synergetic Retrieval of Aerosol Properties (SRAP), which used Terra/Aqua data to produce long-term aerosol optical depth datasets in the Asian region. Lyapustin et al. [
23] used the short-wave infrared 2.1 μm band and BRDF model and proposed the Multiangle implementation of the atmospheric correction (MAIAC) algorithm.
In comparison to polar-orbiting satellites, geostationary satellites can give a wealth of high coverage, high frequency, and multispectral images to monitor aerosol pollution, allowing for the study of aerosol changes on an hourly scale. For aerosol inversion using geostationary satellite data, the most widely used algorithm is the synthetic background method, which assumes that the surface reflectance is stable over time and selected the sub-dark pixel top of atmosphere reflectance (TOA) to obtain surface reflectance within this time window. Knapp et al. [
24] used a synthetic background method; 4 km spatial resolution and 30-min temporal resolution were obtained from GOES. Choi et al. [
25] used the dark blue band (412 nm) to synthesize the background reflectance at each time based on the monthly observation data for GOCI and then realized aerosol inversion. However, due to the existence of clouds, shadows and the complex relationship between the aerosol, and changes in surface and zenith angle, the true absolute surface reflectance is selected that is uncertain. Currently, some studies have moved to the strategy of obtaining the surface reflectance ratio in advance for AOD inversion rather than directly employing absolute reflectance. She et al. [
26] used the stability of BRDF kernel coefficient ratio in a period, combined with the optimal estimation (OE) algorithm, and the synchronous solution of real-time surface reflectance and AOD was realized from the AHI sensor. Su et al. [
27] also proposed the HiPARA algorithm for the AHI sensor; the sensor band was made the atmospheric correction by MOD09 CMA data, then the surface reflectance ratio database was established and achieved high-precision inversion results based on AHI data. Xie et al. [
28] constructed a multi-channel inversion method based on FY-4A data using MODIS AOD data and believed that the ratio of surface reflectance remained unchanged within two weeks and successfully achieved aerosol optical depth inversion in South Asia. The above ratio library methods have achieved very good inversion results. However, the FY-4A currently has no official AOD products available. In this paper, we draw on these ideas, then a new aerosol optical depth retrieval algorithm was proposed for FY-4A/AGRI, which is called the Band Ratio Library (BRL) algorithm.
3. Research Methodology
3.1. Cloud/Water/Snow Pixel Mask Methods
The BRL algorithm can only invert the aerosol optical depth of clear sky pixels over land, so cloud/water pixels need to be eliminated. We used the 4 km cloud products officially released by FY4A for cloud masks based on visible observation and IR observation (
http://satellite.nsmc.org.cn/, accessed on 1 April 2022). In the cloud mask (CLM) product datasets, the cloud states are described as follows, 0: Cloud, 1: Probably Cloud, 2: Probably Clear, and 3: Clear. In this study, we only selected Clear Cloud to ensure the accuracy of the cloud mask and used the MCD12C1 land cover product to remove water pixels. For snow pixels and ephemeral water pixels, we adopt the mask method proposed in this study [
32]. An example diagram of cloud mask results at 04:00 (UTC) on 24 June 2020 is shown in
Figure 2.
3.2. Aerosol Types
The BRL algorithm used five aerosol models, including the continental model (C1), Desert model (C2), Urban model (C3), Maritime model (C4), and Biomass burning model (C5) in the 6S model. Here, C1, C3, and C4 come from a mixture of four basic components (dust-like component, oceanic component, water-soluble component, and soot component), and C2 and C5 come from AERONET measurements [
33,
34]. The five aerosol models are shown in
Table 3.
In general, the aerosol type is unknown before inversion. Because FY-4A/AGRI provides L1 datasets every 15 min, the surface reflectance changes very little at two adjacent moments. Therefore, for optimal aerosol type selection, the following cost function by minimizing the
was adopted [
35]:
In Equation (2), i is the number of observations, λ is the FY-4A spectral channel used for the inversion, and is surface reflectance.
3.3. Inverse Method
In general, for Lambertian surfaces that do not consider surface anisotropy, the equation characterizing the process of the atmospheric radiative transfer model is as follows [
36]:
In Equation (1), , and represent the solar zenith angle, view zenith angle, and relative azimuth angle between them, respectively. and are the total amount of light transmitted downward and upward, respectively. represents the Lambertian surface reflectance, and is the atmospheric path reflectance, including aerosol and Rayleigh scattering. S is the atmosphere spherical albedo.
Before AOD retrieval, an important step is to determine the surface reflectance to separate from the satellite TOA signal. The correction procedure for water vapor and ozone is similar to this study [
37]. We used the method of surface reflectance band ratio library to realize the inversion of AOD. This method has been proven to be effective in retrieving AOD with high accuracy. Firstly, we adopted the MODIS “combined” AOD to do atmosphere correction. In general, the uncertainty of atmospheric corrections increases significantly with increasing aerosol loading [
14,
27,
38]. According to prior studies, the absolute inaccuracy of the MODIS combined AOD rises as the AOD rises; however, there is no discernible shift between 0 and 0.5 [
28]. To guarantee the accuracy of atmospheric correction, these datasets under the condition of AOD < 0.5 was selected for atmospheric correction. Secondly, we selected the “Continental model” aerosol type, then calculated the surface reflectance at 03:00 and 05:00 UTC by running the 6S model. In many previous studies, it is assumed that the two-channel surface reflectance ratio remains unchanged within one month, which can be used to retrieve accurate AOD [
27,
39]. According to Zhang et al. [
40], AOD can be retrieved on bright surfaces with reasonable accuracy if we have good estimates of surface reflectance and can use a single surface ratio value to represent an area of approximately 10 km. This hypothesis was adopted in this study. The monthly average surface reflectance ratio library of Red/NIR obtained is shown in
Figure 3. The band ratio library method has good applicability. Compared with the lack of the DT algorithm in the bright surface area, our algorithm can be applied to the bright surface area. Compared with the DB algorithm for arid and semi-arid regions, our algorithm can dynamically measure the surface reflectance and use the reflectance ratio instead of the absolute reflectance to reduce the error, and variation in scene brightness brought on by topography is one of the reasons why the ratios are more stable than the absolute reflectance [
40].
Next, we used Atmospheric Radiative Transfer Model to estimate the theoretical TOA reflectivity. However, running the model is time-consuming. Therefore, the 6S model is used to make a look-up table (LUT); 6S is a high-precision radiation transmission model; when the zenith angle of solar and view (SZA/VZA) is within 75 degrees, the simulation accuracy of TOA is less than 0.4% [
41]. Therefore, the effective condition for retrieval is that the SZA/VZA is within 75 degrees. The detailed parameters and flowchart are shown in
Table 4 and
Figure 4, respectively.
3.4. Sensitivity Analysis
Figure 5 shows the deviation of the AOD results due to uncertainty of ±5% in the surface reflectance ratio. The sample data were from the location of the Beijing site on June 21, 2020. It can be seen that when the AOD value is higher, the AOD deviation caused by the uncertainty of the ratio is smaller. In the case of low AOD (AOD < 0.5), the ratio error will lead to a significant deviation of the AOD inversion results, which may be related to the dominance of surface signals under low AOD. Note that here the ratio uncertainties are the same for different AOD conditions, but the atmospheric correction uncertainty is lower for low AOD conditions, as mentioned earlier. The accuracy of the atmospheric correction is difficult to guarantee under high AOD conditions. Therefore, this sensitivity analysis does not affect the correctness of atmospheric correction experiments under low AOD conditions.
Figure 6 shows the deviation of the TOA results of different aerosol models in the forward modeling process, where (SZA = 30, VZA = 30, RAA = 30, and SR = 0.25). We can see that when AOD ≤ 0.5, the TOA simulated by C1, C2, and C5 are very close, while C3 and C5 are quite different. C5 is a pure scattering aerosol (SSA = 0.99), while C3 is a strongly absorbing aerosol (SSA = 0.65, of which the soot component accounts for 22%). Of course, if the AOD is very small, the TOA simulated by several aerosol models has almost no difference because the signal of the aerosol is already very weak at this time.
5. Conclusions
There is little research on the retrieval of FY-4A AOD, and the official AOD product of the FY-4A has not yet been made public. In this article, the hour-level AOD over land from FY-4A/AGRI was extracted using a unique, straightforward approach called the BRL algorithm. This algorithm created a library of monthly average band ratios in pixels that may be used to dynamically calculate surface reflectance.
Both MODIS AOD (DT and DB) and AGRI AOD (DB) showed a similar spatial distribution pattern. In regional aerosol studies, AGRI-retrieved AOD datasets exhibit superior temporal resolution compared to MODIS AOD. The spatiotemporal variation of local pollution can be captured using AOD inversion based on geostationary satellite data, and atmospheric pollution can be monitored with high temporal resolution.
However, we still need to consider additional work to make up for some deficiencies. First, this research only uses one month’s data, and the stability of the inversion algorithm needs more time to test, for example, with the characteristics of seasonal and hourly changes. Second, the continental aerosol was used for atmospheric correction in this research since the retrieval results indicated that it was appropriate. Moreover, future research will continue to examine how various aerosol models affect atmospheric correction in the study region.