Satellite retrieval of aerosol combined with assimilated forecast
- 1Japan Aerospace Exploration Agency, Tsukuba, 305-8505, Japan
- 2Research Institute for Applied Mechanics, Kyushu University, Fukuoka, 816-8580, Japan
- 3Atmosphere and Ocean Research Institute, The University of Tokyo, Chiba, 277-8568, Japan
- 4Meteorological Research Institute, Tsukuba, 305-0052, Japan
- apresent address: Remote Sensing Technology Center of Japan, Tsukuba, 305-8505, Japan
Correspondence: Mayumi Yoshida (email@example.com)
We developed a new aerosol satellite retrieval algorithm combining a numerical aerosol forecast. In the retrieval algorithm, the short-term forecast from an aerosol data assimilation system was used as an a priori estimate instead of spatially and temporally constant values. This method was demonstrated using observation of the Advanced Himawari Imager onboard the Japan Meteorological Agency's geostationary satellite Himawari-8. Overall, the retrieval results incorporated strengths of the observation and the model and complemented their respective weaknesses, showing spatially finer distributions than the model forecast and less noisy distributions than the original algorithm. We validated the new algorithm using ground observation data and found that the aerosol parameters detectable by satellite sensors were retrieved more accurately than an a priori model forecast by adding satellite information. Further, the satellite retrieval accuracy was improved by introducing the model forecast instead of the constant a priori estimates. By using the assimilated forecast for an a priori estimate, information from previous observations can be propagated to future retrievals, leading to better retrieval accuracy. Observational information from the satellite and aerosol transport by the model are incorporated cyclically to effectively estimate the optimum field of aerosol.
Aerosols have a fundamental influence on the energy budget of the earth's climate system through the scattering and absorption of solar radiation. The fifth assessment report of the Intergovernmental Panel on Climate Change (IPCC, 2014) stated that radiative forcing of the total aerosol effect in the atmosphere, including cloud adjustments due to aerosols, is −0.9 W m−2. The report also highlighted that the range of uncertainties in these radiative forcing estimations remains large (−1.9 W m−2 to −0.1 W m−2). Identifying the frequency and properties of aerosols over the globe by satellite measurements is essential in estimating the radiation budget and the impacts of aerosols on climate systems.
In satellite aerosol remote sensing, not all aerosol properties can be accurately detected by satellite sensors, as there are more unknown aerosol parameters (e.g., particle size distributions, vertical density distribution, shape, refractive index) than the actual information obtained by the sensors. Most studies use assumptions or information about aerosol parameters and limit the number of parameters retrieved. For example, Higurashi and Nakajima (1999) and Fukuda et al. (2013) assumed fixed complex refractive indices (1.5–0.005i in Higurashi and Nakajima (1999), and 1.503– for small-mode particles and 1.445– for coarse-mode particles based on sulfate and sea spray models, respectively, in Fukuda et al., 2013) and retrieved the aerosol optical thickness and Ångström exponent. Some studies assumed aerosol particle models according to the location and season. For example, Kaufman et al. (1997), Remer et al. (2005), and Levy et al. (2007) estimated aerosol optical thickness and fine-mode fraction over a dark target using the Moderate Resolution Imaging Spectroradiometer (MODIS), by selecting the fine-dominated particle models as a function of geography and season. Hsu et al. (2004) retrieved the optical thickness and type of aerosols over desert regions using blue channels (<500 nm), where the surface reflectance was relatively low, assuming dust or a mixture of dust and smoke depending on the region and season. Jeong et al. (2016) used a priori information according to the location and season. However, these studies did not take temporal changes into account. Because it is impossible to completely fix aerosol type as a function of geographical location and season, the unrealistic assumptions hence implemented lead to one of the major causes of retrieval error.
Aerosol data assimilation methods using satellite data have also been developed to obtain better initial conditions for the aerosol transport model. The aerosol data assimilation study was first developed with low earth orbit (LEO) satellites (Benedetti et al., 2009; Saide et al., 2013; Dai et al., 2014; Rubin et al., 2016; Yumimoto et al., 2015). In recent years, assimilation studies have been extended to using geostationary satellites with large spatial coverage and fine observation frequencies (Saide et al., 2014; Lee et al., 2016; Yumimoto et al., 2016, 2018; Dai et al., 2019; Jin et al., 2019).
Due to the development of such assimilation studies, the satellite data have contributed to improving aerosol forecast simulations. However, no studies have utilized assimilated model forecast as an a priori estimate of the retrieval. Since satellite sensors cannot accurately detect all parameters and unrealistic assumptions of aerosol parameters are a major cause of retrieval errors as mentioned above, adding the model information is expected to improve the retrieval accuracy. Therefore, in this study, we utilize the forecast of an aerosol transport model for a priori estimates of the retrieval. This allows the aerosol information in the aerosol transport model to be used for retrieval. By using the assimilated forecast, information from previous satellite observations can be propagated to future satellite retrievals through the aerosol transport model.
The sections in this study are organized as follows: Sect. 2 explains the retrieval methodology in detail. Section 3.1 presents the results of application to the Advanced Himawari Imager (AHI) onboard Himawari-8. Section 3.2 describes the validation of the estimations using ground observations, and Sect. 3.3 tests the worst-case scenario. Finally, Sect. 4 summarizes our findings.
The aerosol retrieval algorithm in this study is based on Yoshida et al. (2018). As an a priori estimate of the retrieval, the algorithm introduces aerosol forecast from a transport model that has assimilated previous satellite observations. Given the general applicability of the retrieval algorithm by Yoshida et al. (2018), the methodology explained in this section can also be applied to various sensors. Here, we demonstrate the algorithm using the Himawari-8/AHI whose assimilation system is operationally available. The AHI has six observation bands from visible to near-infrared wavelength ranges and observes the top-of-atmosphere (TOA) radiance at a resolution of 0.5–2.0 km over Asia and Oceania at 10 min intervals (Bessho et al., 2016).
Figure 1 depicts an overview of the algorithm, showing the process of using forecast data for a priori estimates of the retrieval. In the original retrieval process, the Level-2 (L2) aerosol optical thickness at 500 nm (τ), the Ångström exponent between 400 and 600 nm (α), and the single-scattering albedo at 500 nm (ω) are retrieved using Level-1 (L1) AHI-observed radiance every 10 min around time T0 as per Yoshida et al. (2018). The Level-3 (L3) τ and α at T0 are then estimated using L2 products in 1 h by an hourly-combined algorithm (Kikuchi et al., 2018). The hourly-combined algorithm is a method that (1) minimizes cloud contamination using the difference between aerosol and cloud spatiotemporal variability characteristics and (2) interpolates the aerosol retrievals using 1 h of data and the movement of clouds within the hour (see Kikuchi et al., 2018, for more details).
The L3 τ at T0 is then assimilated into a global aerosol transport model by the 2D-Var assimilation system (Yumimoto et al., 2018). For the aerosol transport model, we use MASINGAR (Model of Aerosol Species IN the Global AtmospheRe; Tanaka et al., 2003; Tanaka and Chiba, 2005) developed at the Meteorological Research Institute (MRI) of the Japan Meteorological Agency (JMA). MASINGAR covers the major tropospheric aerosol components (i.e., black and organic carbon, mineral dust (10 size bins), sea salt (10 size bins), sulfate aerosols) and their precursors (e.g., sulfur dioxide (SO2), dimethyl sulfide, terpenes), and is coupled online with an atmospheric general circulation model (MRI-AGCM3; Yukimoto et al., 2012). The model's grid resolution is set to horizontal Gaussian TL479 (960×480 grids, about 0.375∘) and 40 vertical layers in hybrid sigma–pressure coordinates from the ground to 0.4 hPa. The integration time step is set to 600 s. Anthropogenic emissions of SO2 and black and organic carbon are taken from the MACCity emission inventory (Granier et al., 2011). Daily biomass burning emission flux is taken from the Global Fire Assimilation System (GFAS, Kaiser et al., 2012) version 1.2 provided by the European Centre of Medium Range Forecast (ECMWF). The horizontal wind components and sea surface temperature are nudged toward the global analyses and forecasts of JMA (GANAL). The forecast from the assimilation system serves as the operational sand and dust forecasting by JMA, the aerosol property model product in the JAXA Himawari Monitor (https://www.eorc.jaxa.jp/ptree/index.html, last access: 20 January 2021) and a member of the ICAP multi-model ensemble (MME) (Xian et al., 2019). The volume concentration (then τ) of each aerosol component at the next time (T1) is then forecasted using the assimilated aerosol transport model.
In the new retrieval process, we retrieve the L2 aerosol properties (τ, α, and ω) from AHI-observed radiance at T1 using these L4 forecasts for a priori estimates of the retrieval. In this way, the information from previous observations at T0 is used for the next aerosol retrievals at T1 through the aerosol transport model. The retrieval obtained at T1 is further used in the same way to derive the retrieval at the following time step (T=T2, not shown) by using L4 forecasts for an a priori estimate. Figure 6 compares the improved retrieval with the original retrieval at T1 as later described in Sect. 3.2. The methodology for using the forecast as a priori estimates of the retrieval is detailed as follows: in the retrieval process, the final retrieval parameters (τ, α, and ω) are calculated from the set of aerosol parameters (τ, external mixing ratio of dry-volume concentration of fine particles ηf, and external mixing ratio of the dry-volume concentration of dust particles for the coarse model ) defined by Yoshida et al. (2018). Here, the imaginary part of the refractive index (mi) for the fine-aerosol model varies with change in such that the fine and coarse models exhibit the same ω at 500 nm (see Yoshida et al., 2018, for more details). The α and ω are calculated from the retrieved ηf and (i.e., mi for fine-aerosol model) using the tables previously calculated by radiative transfer code called the System for the Transfer of Atmospheric Radiation, whose development was initially led by the University of Tokyo (STAR; Nakajima and Tanaka, 1986, 1988; Stamnes et al., 1988). The detailed aerosol setting is explained in Yoshida et al. (2018) and is outlined in Appendix A. Appendix B shows the relationship of the final retrieval parameters (α and ω) with the set of aerosol parameters (ηf and ). We retrieve the aerosol parameters (τ,ηf, and ) using an optimal estimation method (Rodgers, 2000). The state vector of a set of aerosol parameters is derived by minimizing object function J (Eq. 1). It uses the measurement vector of a gas-corrected observed reflectance set , and simulated TOA reflectance , where N is the channel number.
where is the vector of a prior estimate of x and Se and Sa are the covariance matrices of R and xa, respectively. The calculations of R, F(x), and Se are the same as those of Yoshida et al. (2018), but we apply canonical correlation analysis to find the optimal coordinate system and converted R, F(x), and Se to the coordinate system whose dimension is reduced to the number of retrieved parameters (i.e., three). In the original retrieval process, we used spatially and temporally constant values of xa and Sa that were derived from climatology analysis and assumed that the non-diagonal component of covariance matrices was set to 0 (Yoshida et al., 2018).
To introduce a more realistic a priori estimate and covariances into the retrieval process, we employ the forecast from the aerosol assimilation system instead of the constants. The model forecast includes the total aerosol optical thickness at 500 and 870 nm and the absorption aerosol optical thickness at 500 nm derived from the modeled volume concentration and extinction cross section of each aerosol component (Yumimoto et al., 2017). We assign an a priori estimate xa as follows: the model's total aerosol optical thickness at 500 nm is used for τa. is derived from the ratio of total aerosol optical thickness between 500 and 870 nm. As the selection of , we use the model's ω as calculated from the total and absorption aerosol optical thickness at 500 nm.
The assimilation system uses an ensemble method to calculate the background error covariance matrix (Yumimoto et al., 2018). In the method, the ensemble was collected from forecast values within ±2 h of the targeted hour of the five previous forecasts (Fig. 2). We employ this method to define Sa. The model ensemble enables Sa to include the non-diagonal component and express the error of aerosol transport. However, Sa from model ensemble may become too small when the model does not predict the aerosol event itself. For that reason, in order to estimate total Sa, we add a model absolute error () to the error estimated from the ensemble () as follows:
where , , and are the standard deviations of τa, , and , respectively, estimated from the ensemble. , , and are the same as those of the model absolute error.
is set to or (whichever is larger) as follows:
where is the root mean square error (RMSE) of the model's τ from ground observation (0.399 in Fig. 6c) and is the standard deviation of τ for 5 years as calculated by the free-run model without assimilation. The free-run model's spatial resolution is around 1.2∘, and the standard deviation is calculated for MAM (March, April, and May), JJA (June, July, and August), SON (September, October, and November), and DJF (December, January, and February) (Fig. 3). (0.093) is calculated from RMSE of the model's α (0.223 in Fig. 6f) at α of 1.2 and of 0.5. is set to 0.5 because takes a value from 0. to 1, and ω (which is uniquely determined by ) has little correlation with the ground observation data (Fig. 6i). As the non-diagonal component of cannot currently be calculated from our limited database, we use the non-diagonal components of as those of Sa.
3.1 Results of application to Himawari-8
We applied the methodology described in Sect. 2 to the Himawari-8/AHI. We retrieved τ, ηf, and and then derived ω and α at 10 min intervals from the calibrated L1 data subsampled at 0.05∘ using the method described in Sect. 2. The channels used for the retrieval are the same as those used by Yoshida et al. (2018), which are channels 1 (0.46 µm), 2 (0.51 µm), 3 (0.64 µm), 4 (0.86 µm), and 5 (1.6 µm) over land and channels 4 and 5 over the ocean. As the number of satellite channels (two) used over the ocean is less than the number of retrieval parameters (three), not all parameters are stably retrieved by satellite data. Therefore, over the ocean, which is the least sensitive to satellite observation, is set to 0 (i.e., non-absorbing aerosol) at this time because the aerosol over the ocean is generally less absorbing than that over land and using the model's as over the ocean leads to a worse estimation of τ (not shown). However, using non-absorbing aerosol over the ocean causes a big problem in the case of dust or smoke transported over the ocean, so we will use the model's as over the ocean after obtaining a better model's in the future. Note that over land is properly retrieved from satellite data (i.e., not set to 0) using the model's as an a priori estimate, since the number of satellite channels (five) used over land is greater than the number of retrieval parameters (three).
Figures 4 and 5 compare the retrieval results from the new algorithm using for Sa with the original algorithm (Figs. 4, 5a and b). The retrieval results from the new algorithms using only for Sa are also shown to evaluate the effect of on the retrieval result (Figs. 4, 5c). Figure 4 depicts the retrieved τ, ηf, and at 02:00 UTC on 19 May 2016 when aerosols originating from wildfires near Lake Baikal in Russia reached Japan. The model's xa used for retrieval in Fig. 4b and c is indicated in Fig. 4d. The , , and used for retrieval in Fig. 4b are shown in Fig. 4e. The white regions indicate the area where retrieval is not executed due to the presence of clouds, etc. The 2 h forecasts starting from 00:00 UTC on 19 May were assimilated with L3 merged τ at 03:00, 06:00, and 09:00 UTC on 18 May and at 00:00 UTC on 19 May and then used for an a priori estimate (Fig. 4d). Figure 5 is the same as Fig. 4 except for another case at 05:00 UTC on 7 May 2017, when Asian dust was observed in Japan. The 5 h forecast starting from 00:00 UTC on 7 May (and assimilated at 03:00, 06:00, and 09:00 UTC on 6 May and at 00:00 UTC on 7 May) is used for an a priori estimate (Fig. 5d). These short-term forecasts with data assimilation are considered relatively realistic compared to long-term forecasts or a free run without assimilation (Yumimoto et al., 2018). If only model ensemble error is used for Sa (Figs. 4c, 5c), that is, the absolute error is not included in Sa, all retrieved parameters (especially ηf and over land) are highly dependent on an a priori estimate (Figs. 4d, 5d). However, when using an appropriate Sa containing absolute error, the retrieved τ, ηf, and are updated by satellite data or remain close to an a priori estimate depending on the location (Figs. 4b, 5b). Specifically, spatially finer τ distributions than the model forecast are retrieved for cases of both wildfire aerosol (Fig. 4b) and Asian dust (Fig. 5b) due to the relatively coarser model horizontal resolution. Similar ηf is retrieved over open ocean in Figs. 4b and Fig 5b, and the large ηf (i.e., small particle) and small ηf (i.e., large particle) are successfully retrieved in areas corresponding to wildfire aerosol (Fig. 4b) and Asian dust (Fig. 5b), respectively. This distribution is also expressed in the forecast model in Fig. 5d but cannot be expressed sufficiently in the forecast model in Fig. 4d because information about the aerosol particle size (e.g., α, ηf) is not assimilated into the model. That is, by using an appropriate Sa, both the model and satellite data are used for estimating the aerosol parameters. In addition, the local noise in τ and ηf is apparently reduced for this algorithm (Figs. 4b, 5b) as compared with the original algorithm (Figs. 4a, 5a). This will be discussed in Sect. 3.2.
We conducted a preliminary validation of our method by comparing the retrieved τ, α, and ω from the Himawari-8/AHI with those from ground observation known as the Aerosol Robotic Network (AERONET). AERONET's τ and α were derived from Level 2.0 quality-assured Version 3 direct sun algorithm data (Giles et al., 2019; O'Neill et al., 2003), and ω was derived from Version 3 AERONET inversion products (Dubovik and King, 2000; Dubovik et al., 2000; Dubovik et al., 2002a, b, 2006; Sinyuk et al., 2007). AERONET's ω at 500 nm was calculated from linear interpolation of ω at 440 and 675 nm. In this study, the 60 AERONET sites on the full disk of Himawari-8 were used for the validation. We used the AERONET data averaged over 10 min of AHI observation time. For our retrieval data, we used τ, α, and ω estimated from AHI L1 radiance data subsampled at 0.05∘ nearest to the AERONET sites. Initial validation was conducted for 6 months (March, April, May, June, and July 2018 and February 2019). Long-term validation will be required in future studies.
Figure 6 compares the τ, α, and ω retrieved from the AHI with those from AERONET. The validation of α and ω is limited to cases where the simultaneously retrieved τ are greater than 0.3 because there is little information of α and ω from satellite observation for thin aerosol layers. The total number of validation points (14 711, 14 031, and 521) from this algorithm (Fig. 6b, e, and h) is about 6 %–7 % higher than those (13 714, 13 137, and 493) from the original algorithm (Fig. 6a, d, and g), which means that the new algorithm successfully retrieved the aerosol in more cases than the original algorithm. Here, the total number of validation points for ω is less than those for τ and α because the number of ω data from AERONET inversion products is less than those of τ and α from the direct sun measurements.
For the τ estimations (Fig. 6a, b, c), the root mean square error (RMSE), mean bias (MB), and correlation (0.290, −0.099, and 0.758) from this algorithm (Fig. 6b) are all better than those (0.399, −0.224, and 0.572) from the model forecast (i.e., an a priori estimate) in Fig. 6c, which means that satellite information is very effective for the retrieval of τ. In addition, the RMSE (0.290) in Fig. 6b is better than that (0.307) without the forecast model (Fig. 6a), which means that the model information is also effective and the improved algorithm shows better performance than the original algorithm. The MB (−0.099) in Fig. 6b is worse than that (−0.023) in Fig. 6a, probably because the large outlier in Fig. 6a is improved in Fig. 6b. Figure 7 shows an example of the retrieval results of the outlier (red asterisks in Fig. 6a, b, and c) and shows that the outlier of the original algorithm (Fig. 7a) is improved in the new algorithm by constraining τ to the model's τa. In addition, the retrieval results around the red circles show that the new algorithm successfully retrieved the τ, ηf, and even where the original algorithm failed to retrieve them. Thus, integrating the model and satellite information resulted in an improvement of the τ estimations.
For the α estimations (Fig. 6d, e, f), large variance in the original method is considerably reduced by this method. The RMSE and correlation (0.271 and 0.581) from this algorithm (Fig. 6e) are much better than those (0.429 and 0.353) from the original algorithm without the forecast model (Fig. 6d), which indicates that the new algorithm could improve the precision of α estimations by adding more accurate α (RMSE of 0.223) information from the model. In addition, the MB (−0.052) from this algorithm (Fig. 6e) is better than that (−0.057) from the model forecast (Fig. 6f), due to the improvement of negative bias in the large α in the model forecast. Thus, the new α can be retrieved with good accuracy by utilizing the relatively accurate model's α and correcting the bias by using the satellite data.
For the ω estimations (Fig. 6g, h, i), the RMSE, MB, and correlation (0.035, 0.000, and 0.550) from this algorithm (Fig. 6h) are better than those (0.048, −0.002, and 0.176) from the model forecast (Fig. 6i), which indicates the effectiveness of satellite information for ω retrieval. In addition, this algorithm improved RMSE, MB, and correlation by introducing the model forecast. Note that the current system assimilates only total τ (total amount of aerosols), and information about the fraction of fine and absorbing aerosols from the retrieval is not fed back to the model forecast. The improved retrieval accuracy of ω can be expected if the model's ω becomes more realistic in the future, such as by assimilating the satellite's ω to the model. Considering the validation results of τ, α, and ω, this new algorithm effectively improved the retrieval accuracy using information from both the model and the satellite by setting appropriate Sa and Se.
We also investigated the cause of the possible large deviation between the retrieved parameters from the new algorithm and the ground observation. Figures 8, 9, and 10 show the validation results of τ, α, and ω, respectively, when the chi-square value (χ2) and the uncertainties of the retrieved three parameters (τ, ηf, and ) are smaller than a threshold. The chi-square value (χ2) is calculated as follows:
It shows the closeness of the retrieved value to the observed value. The covariance matrix of the uncertainties of the retrieved parameters is calculated using the law of error propagation, as follows:
where A is the Jacobian matrix. Se is the covariance matrix of R and calculated from sum of sensor noise and the uncertainty in TOA reflectance that results from surface reflectance uncertainty (Yoshida et al., 2018). In reality, the Se is almost determined by the uncertainty in TOA reflectance that results from surface reflectance uncertainty because sensor noise is much smaller. Therefore, the is mostly caused by the surface reflectance uncertainty. Figure 8 shows that RMSE for τ decreases as the threshold of χ2 or Sτ becomes strict (i.e., decreases). On the other hand, RMSE for α (in Fig. 9) is not dependent on the threshold of but decreases as the χ2 threshold decreases. RMSE for ω (in Fig. 10) is almost independent on the threshold of and χ2. Next, in Fig. 11 we investigated how the retrieved accuracy (difference between aerosol parameters retrieved from AHI and those of AERONET) depends on the model's (i.e., a priori) accuracy. The retrieved accuracy of α and ω has strong linear relationships (a correlation of 0.801 and 0.739, respectively) to the model's accuracy, while that of τ has a moderate linear relationship (a correlation of 0.622). Summarizing these results, the retrieved accuracy of τ depends on all of the closeness to the observed value, accuracy of the surface reflectance estimation, and accuracy of an a priori estimate, while the accuracy of an a priori estimate is critical for the retrieved accuracy of α and ω. Thus, introducing more realistic a priori estimates into the new retrieval algorithm instead of the constant values in the original algorithm led to the improvement of RMSE. It is also shown that the improvement of a numerical aerosol forecast by improving the aerosol transport model and the assimilation method and increasing the assimilation frequency may further improve the retrieval accuracy in the future.
3.3 Worst-case scenario
We have shown that the new retrieval algorithm using the forecast of an aerosol transport model improves the retrieval accuracy. However, in order to use this algorithm constantly (such as in an operational system), the effects of the model forecast (a priori estimate) that deviate from reality must be examined because the model forecast may miss an aerosol event. Therefore, we conducted a sensitivity test to investigate the impact on the retrieval results of using unrealistic forecast as an a priori estimate. Figure 12 shows the retrieval results on the same day as in Fig. 4, except for using the forecast on another day (27 April 2018) as an a priori estimate of the retrieval (Fig. 12d). If only is used as Sa (Fig. 12c), all parameters (especially ηf and ) are retrieved unrealistically by being dependent on the unrealistic a priori estimate. However, when using an appropriate Sa (Eq. 2), the retrieved parameters are updated well by satellite data with less dependence on an unrealistic a priori estimate (Fig. 12b). Even in such an extremely worst-case scenario, this new algorithm is apparently not significantly worse than the current algorithm, especially where the model forecast is missing an aerosol event, which may occur in the model forecast for natural aerosols (e.g., mineral dust and smoke from biomass burning).
We developed a new satellite aerosol retrieval algorithm combining a numerical aerosol forecast. In the retrieval algorithm, the short-term forecast from an aerosol data assimilation system was used for an a priori estimate instead of spatially and temporally constant values. This is the first study that utilizes the assimilated model forecast of aerosol as an a priori estimate of the satellite retrieval. We applied this new algorithm to the Himawari-8/AHI and confirmed that the aerosol parameters detectable by satellite sensors were retrieved more accurately (RMSE of 0.290 for τ and 0.035 for ω) than an a priori model forecast (RMSE of 0.399 for τ and 0.048 for ω) by adding satellite information. Moreover, the satellite retrieval accuracy was improved (RMSE of 0.290 for τ, 0.271 for α, and 0.035 for ω) by using the model forecast as compared with those using constant a priori estimates (RMSE of 0.307 for τ and 0.429 for α and 0.039 for ω). As a result, aerosol retrievals were improved by effectively incorporating both model and satellite information, depending on each covariance. By using the assimilated forecast as an a priori estimate, information from previous observations can be propagated to future retrievals, thereby leading to better retrieval accuracy. In this way, satellite observation and model simulation are used synergistically to continuously estimate the optimum field of aerosol. Future work would include applying the methodology proposed in this study to polar-orbiting satellites and combining them with geostationary satellite measurements, in order to offer consistent geostationary and polar-orbiting estimates and thereby improve aerosol properties over the globe.
We assume that the aerosol model is an external mixture of fine and coarse particles (ηf is the external mixing ratio of the dry-volume concentration of fine particles). We set the fine-aerosol model based on the average properties of the fine mode for categories 1–6 by Omar et al. (2005). For the coarse-aerosol model, we set the external mixture of the pure marine aerosol on the basis of the model illustrated by Sayer et al. (2012) and the dust model based on the coarse model of category 1 (dust) as illustrated by Omar et al. (2005). is the external mixing ratio of the dry-volume concentration of dust particles for the coarse model.
Regarding each aerosol size, we use a monomodal lognormal volume size (rd) distribution, which is defined as follows:
where Cv is the particle volume concentration, rv is the volume median radius, and σ is the standard deviation. rv is set to 0.143, 2.59, and 2.834 (σ is 1.537, 2.054, and 1.908) for fine, coarse marine, and coarse dust, respectively, based on the observations by Omar et al. (2005) and Sayer et al. (2012). Regarding the aerosol shape, we assume a spherical model for the fine and coarse marine models and a non-spherical model for the coarse dust model (Nakajima et al., 1989). The aerosol vertical distribution is set to the same distribution that was used for rural (dominant at 0–2 km), sea spray (below 2 km), and yellow sand (4–8 km) for fine, coarse marine, and coarse dust in the STAR code, respectively. The real part of the refractive index is set to 1.439, 1.362, and 1.452 for fine, coarse marine, and coarse dust, respectively, and the imaginary part of the refractive index (mi) is set to and 0.0036 at all wavelengths for coarse marine, and coarse dust, respectively, based on Sayer et al. (2012) and Omar et al. (2005). The mi for the fine-aerosol model is perturbed to represent non-absorbing and absorbing aerosols. To decrease the number of derived parameters, the mi for the fine-aerosol model varies with change in such that the fine and coarse models exhibit the same ω at 500 nm.
Figure B1 shows the relations of the final retrieval parameters α and ω with the external mixing ratio of dry-volume concentration of fine particles (ηf) and external mixing ratio of the dry-volume concentration of dust particles for the coarse model (). The ω at 500 nm can be uniquely determined by the (Fig. B1a), since for the coarse aerosol changes in conjunction with mi for the fine aerosol so that the ω at 500 nm has the same value without depending on the ηf. Note that the ω at wavelengths other than 500 nm are dependent on not only but also ηf. The α is mainly determined by ηf but also depends slightly on (Fig. B1b).
Himawari-8/AHI aerosol data are available from the JAXA Himawari Monitor (2020) site: https://www.eorc.jaxa.jp/ptree/index.html (last access: 20 January 2021). The retrieved data in this study can be requested directly from the lead author (firstname.lastname@example.org).
MY developed the retrieval code and analyzed the data with significant conceptual input from KY and critical technical support from TMN, MK, and HM. KY and TYT prepared the assimilated forecast data for test cases and long time validations, respectively. MY prepared the paper with contributions from all co-authors.
The authors declare that they have no conflict of interest.
The authors are grateful to the Open CLASTER project for allowing us to use the RSTAR package for this research. We would like to thank the AERONET project and its staff for establishing and maintaining the AERONET sites considered in this investigation. We also thank Haruma Ishida for providing the cloud detection algorithm (CLAUDIA). Finally, we appreciate the valuable discussions and support provided by Takashi Maki, Tsuyoshi Sekiyama, Makiko Hashimoto, and Teruyuki Nakajima.
This research has been supported by the JSPS Grants-in-Aid for Scientific Research (grant no. JP16H02946) and the JAXA Earth Observation Priority Research (grant no. PI.ER2GCF212).
This paper was edited by Stelios Kazadzis and reviewed by Alexander Kokhanovsky and two anonymous referees.
Benedetti, A., Morcrette, J.-J., Boucher, O., Dethof, A., Engelen, R. J., Fisher, M., Flentje, H., Huneeus, N., Jones, L., Kaiser, J. W., Kinne, S., Mangold A., Razinger, M. , Simmons, A. J., and Suttie, M.: Aerosol analysis and forecast in the European Centre for Medium-Range Weather Forecasts Integrated Forecast System: 2. Data assimilation, J. Geophys. Res., 114, D13205, https://doi.org/10.1029/2008JD011115, 2009.
Bessho, K., Date, K., Hayashi, M., Ikeda, A., Imai, T., Inoue, H., Kumagai, Y., Miyakawa, T., Murata, H., Ohno, T., Okuyama, A., Oyama, R., Sasaki, Y., Shimazu, Y., Shimoji, K., Sumida, Y., Suzuki, M., Taniguchi, H., Tsuchiyama, H., Uesawa, D., Yokota, H., and Yoshida, R.: An introduction to Himawari-8/9 – Japan's new-generation geostationary meteorological satellites, J. Meteor. Soc. Japan, 94, 151–183, 2016.
Dai, T., Schutgens, N. A. J., Goto, D., Shi, G., and Nakajima, T.: Improvement of aerosol optical properties modeling over Eastern Asia with MODIS AOD assimilation in a global non-hydrostatic icosahedral aerosol transport model, Environ. Pollut., 195, 319–329, https://doi.org/10.1016/j.envpol.2014.06.021, 2014.
Dai, T., Cheng, Y., Suzuki, K, Goto, D., Kikuchi, M., Schutgens, N., Yoshida, M., Zhang, P., Husi, L., Guangyu, S., and Nakajima, T.: Hourly aerosol assimilation of Himawari-8 AOT using the four-dimensional local ensemble transform Kalman filter, J. Adv. Model. Earth Syst., 11, 680–711, https://doi.org/10.1029/2018MS001475, 2019.
Dubovik, O. and King, M. D.: A flexible inversion algorithm for retrieval of aerosol optical properties from Sun and sky radiance measurements, J. Geophys. Res., 105, 20673–20696, 2000.
Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kaufman, Y. J., Eck, T. F., and Slutsker, I.: Accuracy assessment of aerosol optical properties retrieval from AERONET sun and sky radiance measurements, J. Geophys. Res., 105, 9791–9806, 2000.
Dubovik, O., Holben, B. N., Eck, T. F., Smirnov, A., Kaufman, Y. J., King, M. D., Tanré, D., and Slutsker, I.: Variability of absorption and optical properties of key aerosol types observed in worldwide locations, J. Atmos. Sci., 59, 590–608, 2002a.
Dubovik, O., Holben, B. N., Lapyonok, T., Sinyuk, A., Mishchenko, M. I., Yang, P., and Slutsker, I.: Non-spherical aerosol retrieval method employing light scattering by spheroids, Geophys. Res. Lett., 10, 54-1–54-4, https://doi.org/10.1029/2001GL014506, 2002b.
Dubovik, O., Sinyuk, A., Lapyonok, T., Holben, B. N., Mishchenko, M., Yang, P., Eck, T. F., Volten, H., Munoz, O., Veihelmann, B., van der Zander, Sorokin, M., and Slutsker, I.: Application of light scattering by spheroids for accounting for particle non-sphericity in remote sensing of desert dust, J. Geophys. Res., 111, D11208, https://doi.org/10.1029/2005JD006619, 2006.
Fukuda, S., Nakajima, T., Takenaka, H., Higurashi, A., Kikuchi, N., Nakajima, T. Y., and Ishida, H.: New approaches to removing cloud shadows and evaluating the 380 nm surface reflectance for improved aerosol optical thickness retrievals from the GOSAT/TANSO-Cloud and Aerosol Imager, J. Geophys. Res.-Atmos., 118, 13520–13531, 2013.
Giles, D. M., Sinyuk, A., Sorokin, M. G., Schafer, J. S., Smirnov, A., Slutsker, I., Eck, T. F., Holben, B. N., Lewis, J. R., Campbell, J. R., Welton, E. J., Korkin, S. V., and Lyapustin, A. I.: Advancements in the Aerosol Robotic Network (AERONET) Version 3 database – automated near-real-time quality control algorithm with improved cloud screening for Sun photometer aerosol optical depth (AOD) measurements, Atmos. Meas. Tech., 12, 169–209, https://doi.org/10.5194/amt-12-169-2019, 2019.
Granier, C., Bessagnet, B., Bond, T., D'Angiola, A., Denier van der Gon, H., Frost, G. J., Heil, A., Kaiser, J. W., Kinne, S., Klimont, Z., Kloster, S., Lamarque, J.-F., Liousse, C., Masui, T., Meleux, F., Mieville, A., Ohara, T., Raut, J.-C., Riahi, K., Schultz, M. G., Smith, S. J., Thompson, A., van Aardenne, J., van der Werf, G. R., and van Vuuren, D. P.: Evolution of anthropogenic and biomass burning emissions of air pollutants at global and regional scales during the 1980–2010 period, Clim. Change, 109, 163–190, https://doi.org/10.1007/s10584-011-0154-1, 2011.
Higurashi, A. and Nakajima T.: Development of a two-channel aerosol retrieval algorithm on a global scale using NOAA AVHRR, J. Atmos. Sci., 56, 924–941, 1999.
Hsu, N. C., Tsay, S.-C., King, M. D., and Herman J. R.: Aerosol properties over bright-reflecting source regions, IEEE T. Geosci. Remote, 42, 557–569, 2004.
IPCC, Climate Change: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Core Writing Team, Pachauri, R. K., and Meyer, L. A., Geneva, Switzerland, 151 pp., 2014.
JAXA Himawari Monitor: P-Tree System, available at: https://www.eorc.jaxa.jp/ptree/index.html (last access: 20 January 2021), 2020.
Jeong, U., Kim, J., Ahn, C., Torres, O., Liu, X., Bhartia, P. K., Spurr, R. J. D., Haffner, D., Chance, K., and Holben, B. N.: An optimal-estimation-based aerosol retrieval algorithm using OMI near-UV observations, Atmos. Chem. Phys., 16, 177–193, https://doi.org/10.5194/acp-16-177-2016, 2016.
Jin, J., Segers, A., Heemink, A., Yoshida, M., Han, W., and Lin, H.-X: Dust emission inversion using Himawari-8 AODs over East Asia: An extreme dust event in May 2017, J. Adv. Model. Earth Syst., 11, 446–467, https://doi.org/10.1029/2018MS001491, 2019.
Kaiser, J. W., Heil, A., Andreae, M. O., Benedetti, A., Chubarova, N., Jones, L., Morcrette, J.-J., Razinger, M., Schultz, M. G., Suttie, M., and van der Werf, G. R.: Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power, Biogeosciences, 9, 527–554, https://doi.org/10.5194/bg-9-527-2012, 2012.
Kaufman, Y. J., Tanré, D., Remer, L. A., Vermote, E. F., Chu, A., and Holben, B. N.: Operational remote sensing of tropospheric aerosol over land from EOS moderate resolution imaging spectroradiometer, J. Geophys. Res.-Atmos., 102, 17051–17067, 1997.
Kikuchi M., Murakami, H., Suzuki, K., Nagao, T. M., and Higurashi, A.: Improved Hourly Estimates of Aerosol Optical Thickness using Spatiotemporal Variability Derived from Himawari-8 Geostationary Satellite, IEEE T. Geosci. Remote, 2018, 3442–3455, https://doi.org/10.1109/TGRS.2018.2800060, 2018.
Lee, S., Song, C. H., Park, R. S., Park, M. E., Han, K. M., Kim, J., Choi, M., Ghim, Y. S., and Woo, J.-H.: GIST-PM-Asia v1: development of a numerical system to improve particulate matter forecasts in South Korea using geostationary satellite-retrieved aerosol optical data over Northeast Asia, Geosci. Model Dev., 9, 17–39, https://doi.org/10.5194/gmd-9-17-2016, 2016.
Levy, R. C., Remer, L. A., Mattoo, S., Vermote, E. F., and Kaufman Y. J.: Second-generation operational algorithm: Retrieval of aerosol properties over land from inversion of Moderate Resolution Imaging Spectroradiometer spectral reflectance, J. Geophys. Res., 112, D13211, https://doi.org/10.1029/2006JD007811, 2007.
Nakajima, T. and Tanaka, M.: Matrix formulation for the transfer of solar radiation in a plane-parallel scattering atmosphere, J. Quant. Spectrosc. Ra., 35, 13–21, 1986.
Nakajima, T. and Tanaka, M.: Algorithms for radiative intensity calculations in moderately thick atmospheres using a truncation approximation, J. Quant. Spectrosc. Ra., 40, 51–69, 1988.
Nakajima, T., Tanaka, M., Yamano, M., Shiobara, M., Arao, K., and Nakanishi, Y.: Aerosol optical characteristics in the yellow sand events observed in May, 1982 in Nagasaki – part II model, J. Meteorol. Soc. Jpn., 67, 279–291, 1989.
Omar, A. H., Won, J.-G., Winker, D. M., Yoon, S.-C., Dubovik, O., and McCormick, M. P.: Development of global aerosol models using cluster analysis of Aerosol Robotic Network (AERONET) measurements, J. Geophys. Res., 110, D10S14, https://doi.org/10.1029/2004JD004874, 2005.
O'Neill, N. T., Eck, T. F., Smirnov, A., Holben, B. N., and Thulasiraman, S.: Spectral discrimination of coarse and fine mode optical depth, J. Geophys. Res., 108, 4559–4573, https://doi.org/10.1029/2002JD002975, 2003.
Remer, L. A., Kaufman, Y. J., Tanré, D., Mattoo, S., Chu, D. A., Martins, J. V., R.-R. Li, Ichoku, C., Levy, R. C., Kleidman, R. G., Eck, T. F., Vermote, E., and Holben, B. N.: The MODIS aerosol algorithm, products, and validation, J. Atmos. Sci., 62, 947–973, 2005.
Rodgers, C. D.: Inverse Method for Atmospheric Sounding: Theory and Practice, World Scientific, Singapore, 238 pp., 2000.
Rubin, J. I., Reid, J. S., Hansen, J. A., Anderson, J. L., Collins, N., Hoar, T. J., Hogan, T., Lynch, P., McLay, J., Reynolds, C. A., Sessions, W. R., Westphal, D. L., and Zhang, J.: Development of the Ensemble Navy Aerosol Analysis Prediction System (ENAAPS) and its application of the Data Assimilation Research Testbed (DART) in support of aerosol forecasting, Atmos. Chem. Phys., 16, 3927–3951, https://doi.org/10.5194/acp-16-3927-2016, 2016.
Saide, P. E., Carmichael, G. R., Liu, Z., Schwartz, C. S., Lin, H. C., da Silva, A. M., and Hyer, E.: Aerosol optical depth assimilation for a size-resolved sectional model: impacts of observationally constrained, multi-wavelength and fine mode retrievals on regional scale analyses and forecasts, Atmos. Chem. Phys., 13, 10425–10444, https://doi.org/10.5194/acp-13-10425-2013, 2013.
Saide, P. E., Kim, J., Song, C. H., Choi, M., Cheng, Y. and Carmichael, G. R.: Assimilation of next generation geostationary aerosol optical depth retrievals to improve air quality simulations, Geophys. Res. Lett., 41, 9188–9196, https://doi.org/10.1002/2014GL062089, 2014.
Sayer, A. M., Smirnov, A., Hsu, N. C., and Holben, B. N.: A pure marine aerosol model, for use in remote sensing applications, J. Geophys. Res., 117, D05213, https://doi.org/10.1029/2011JD016689, 2012.
Sinyuk, A., Dubovik, O., Holben, B., Eck, T. F., Breon, F.-M., Martonchik, J., Kahn, R., Diner, D. J., Vermote, E. F., Roger, J.-C., Lapyonok, T., and Slutsker, I.: Simultaneous retrieval of aerosol and surface properties from a combination of AERONET and satellite, Rem. Sens. of Env., 107, 90–108, https://doi.org/10.1016/j.rse.2006.07.022, 2007.
Stamnes, K., Tsay, S. C., Wiscombe, W., and Jayaweera, K.: Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media, Appl. Optics, 27, 2502–2509, 1988.
Tanaka, T. Y. and Chiba, M.: Global simulation of dust aerosol with a chemical transport model, MASINGAR, J. Meteorol. Soc. Jpn., 83A, 255–278, https://doi.org/10.2151/jmsj.83A.255, 2005.
Tanaka, T. Y., Orito, K., Sekiyama, T. T., Shibata, K., Chiba, M., and Tanaka, H.: MASINGAR, a global tropospheric aerosol chemical transport model coupled with MRI/JMA98 GCM: Model description, Pap. Meteorol. Geophys., 53, 119–138, https://doi.org/10.2467/mripapers.53.119, 2003.
Xian, P., Reid, J. S., Hyer, E. J., Sampson, C. R., Rubin, J. I., Ades, M., Asencio, N., Basart, S., Benedetti, A., Bhattacharjee, P., Brooks, M. E., Colarco, P. R., Da Silva, A., Eck, T. F., Guth, J., Jorba, O., Kouznetsov, R., Kipling, Z., Sofiev, M., Carcia-Pando, C. P., Pradhan, Y., Tanaka, T., Wang, J., Westphal, D. L., Yumimoto, K., and Zhang, J.: Current State of the global operational aerosol multi-model ensemble: an update from the International Cooperative for Aerosol Prediction (ICAP), Q. J. Roy. Meteor. Soc., 145, qj.3497, https://doi.org/10.1002/qj.3497, 2019.
Yoshida, M., Kikuchi, M., Nagao, T. M., Murakami, H., Nomaki, T., and Higurashi, A.: Common retrieval of atmospheric aerosol properties for imaging satellite sensors, J. Meteorol. Soc. Jpn., 96B, 193–209, https://doi.org/10.2151/jmsj.2018-039, 2018.
Yukimoto, S., Adachi, Y., Hosaka, M., Sakami, T., Yoshimura, H., Hirabara, M., Tanaka, T. Y., Shindo, E., Tsujino, H., Deushi, M., Mizuta, R., Yabu, S., Obata, A., Nakano, H., Koshiro, T., Ose, T., and Kitoh, A.: A new global climate model of the Meteorological Research Institute: MRI-CGCM3 – model description and basic performance, J. Meteorol. Soc. Jpn., 90A, 23–64, https://doi.org/10.2151/jmsj.2012-A02, 2012.
Yumimoto, K., Murakami, H., Yanaka, T. Y., Sekiyama, T. T., Ogi, A., and Maki, T.: Forecasting of Asian dust storm that occurred on May 10–13, 2011, using an ensemble-based data assimilation system, Particuology, 28, 121–130, https://doi.org/10.1016/j.partic.2015.09.001, 2015.
Yumimoto, K., Nagao, T. M., Kikuchi, M., Sekiyama, T. T., Murakami, H., Tanaka, T. Y., and Maki, T.: Aerosol data assimilation using data from Himawari-8, a next-generation geostationary meteorological satellite, Geophys. Res. Lett., 43, 5886–5894, https://doi.org/10.1002/2016GL069298, 2016.
Yumimoto, K., Tanaka, T. Y., Oshima, N., and Maki, T.: JRAero: the Japanese Reanalysis for Aerosol v1.0, Geosci. Model Dev., 10, 3225–3253, https://doi.org/10.5194/gmd-10-3225-2017, 2017.
Yumimoto, K., Tanaka, T., Yoshida, M, Kikuchi, M., Nagao, T. M., Murakami, H., and Maki, T.: Assimilation and Forecasting Experiment for Heavy Siberian Wildfire Smoke in May 2016 with Himawari-8 Aerosol Optical Thickness, J. Meteorol. Soc. Jpn., 96B, 133–149, https://doi.org/10.2151/jmsj.2018-035, 2018.