Satellite-Based Aerosol Classification for Capital Cities in Asia Using a Random Forest Model

Aerosol types in Asian capital cities were classified using a random forest (RF) satellite-based aerosol classification model during 2018–2020 in an investigation of the contributions of aerosol types, with or without Aerosol Robotic Network (AERONET) observations. In this study, we used the recently developed RF aerosol classification model to detect and classify aerosols into four types: pure dust, dust-dominated aerosols, strongly absorbing aerosols, and non-absorbing aerosols. Aerosol optical and microphysical properties for each aerosol type detected by the RF model were found to be reasonably consistent with those for typical aerosol types. In Asian capital cities, pollution-sourced aerosols, especially non-absorbing aerosols, were found to predominate, although Asian cities also tend to be seasonally affected by natural dust aerosols, particularly in East Asia (March–May) and South Asia (March–August). No specific seasonal effects on aerosol type were detected in Southeast Asia, where there was a predominance of non-absorbing aerosols. The aerosol types detected by the RF model were compared with those identified by other aerosol classification models. This study indicates that the satellite-based RF model may be used as an alternative in the absence of AERONET sites or observations.


Introduction
Aerosol characteristics such as particle size and absorption properties vary with their emission sources [1], and it is difficult to define categories of aerosol type due to their diverse properties and temporospatial variability [2][3][4]. In climate research, aerosol type is an important input parameter in calculating radiative forcing, and can cause uncertainties in outputs [5][6][7][8]. Aerosol type is also a useful input in satellite aerosol retrieval algorithms [9,10] and transport models [11,12].
In the investigation of predominant aerosol types on regional scales, several classification methods have been proposed using ground-based Aerosol Robotic Network (AERONET) observations, which provide aerosol optical (e.g., aerosol optical depth, single scattering albedo, and phase function) and microphysical (e.g., size distribution) properties. Gobbi et al. [13] classified aerosol types for unsupervised classification based on the Ångström exponent (AE) and its variations, as well as aerosol optical depth (AOD). Other studies introduced cluster analysis with inputs for aerosol properties such as refractive index, AE, mean radius, asymmetry factors, and single scattering albedo (SSA) [14][15][16]. In threshold-based aerosol classification, AOD, SSA, and aerosol particle size (AE and fine mode fraction; FMF) are commonly involved [17][18][19][20][21]. These threshold-based aerosol classification methods were applied to investigate aerosol types on a global scale [17], Southeast Asia [18], Beijing [19], Europe, the Middle East, North Africa, the Arabian Peninsula [20], and Saudi Arabia [21]. The depolarization ratio, which is sensitive to aerosol particle shape,

Method
The overall flow chart of this study is shown in Figure 1. Aerosol types at 13:30 local time (overpass time of satellite inputs) were classified and seasonally and annually investigated in several capital cities in Asia for the period from 2018 to 2020. Locations of the regions of interest (Asian capital cities) are specified in Section 2.1. To classify aerosol types in Asian capital cities, we used the RF model that had been globally trained using a dataset for January 2018 to June 2020 [31], as described in Section 2.2. To test the feasibility of the RF model in Asian capital cities, we first investigated AERONET aerosol optical (spectral dependence of SSA, depolarization ratio, and FMF) and microphysical (volume size distribution and effective radius) properties for aerosol types detected by the RF model. The RF-based aerosol types of Asian capital cities were compared with those of the AERONET-based and other satellite aerosol classification models ( Figure 1). The descriptions of AERONET and satellite data used in this study are in Section 2.3. Remote Sens. 2021, 13, x FOR PEER REVIEW 3 of 14

Regions of Interest
The regions of interest were selected from capital cities in East Asia, South Asia, and Southeast Asia with populations of >1 million. About 40% of the 19 capitals lack an AER-ONET observation site (Table 1). To simultaneously compare AERONET and satellite observation data, the location of the AERONET site was used as the location of the corresponding city, with satellite data for near that location being collected. For cities with no AERONET site, coordinates were obtained from the GeoHack system of Wikitech (https://geohack.toolforge.org/, accessed on 20 April 2021). Table 1. Regions of interest containing capital cities in East Asia, South Asia, and Southeast Asia with populations of >1,000,000 [47]. The coordinates of cities are those of an AERONET site; coordinates of cities without AERONET sites were obtained from the GeoHack system of Wikitech (https://geohack.toolforge.org/, accessed on 20 April 2021).

Regions of Interest
The regions of interest were selected from capital cities in East Asia, South Asia, and Southeast Asia with populations of >1 million. About 40% of the 19 capitals lack an AERONET observation site (Table 1). To simultaneously compare AERONET and satellite observation data, the location of the AERONET site was used as the location of the corresponding city, with satellite data for near that location being collected. For cities with no AERONET site, coordinates were obtained from the GeoHack system of Wikitech (https://geohack.toolforge.org/, accessed on 20 April 2021).

Random Forest Model for Satellite Aerosol Classification
We used the satellite classification model developed previously [31] and based on the RF approach, which builds an ensemble of decision trees by making predictions determined by a majority vote of each tree [48]. Choi et al. [31] trained using a set of observational data including target (AERONET aerosol type) and input (satellite measurement) variables to detect aerosol types without input of AERONET data. Satellite input parameters of the RF model include 6 satellite-based products: aerosol index, column densities of trace gases (CO and tropospheric NO 2 ), solar zenith angle, land-cover type, and percentage urban area. For example, column amounts of CO and tropospheric NO 2 were selected for identifying the presence of carbonaceous particles and scattering-dominant aerosols, respectively [30,31]. In the case of sulfur oxides (SOx), SO 2 is the only one that can be measured by satellite observations, but it was excluded due to the greater uncertainty in satellite SO 2 column data [49] compared with that of column CO [50] and tropospheric NO 2 [51] data. The aerosol index is the parameter for identifying absorbing and non-absorbing aerosols, but it is not the only parameter useful for identifying dust aerosols [21]. The RF model was found to have sensitivity for distinguishing dust aerosols by training with satellite input variables (e.g., aerosol index, column amount of trace gases) and the AERONET-based aerosol type dataset [31]. The training dataset of the model included 60% of data collected globally from January 2018 to June 2020 [31]. Optimal hyperparameters were selected based on 5-fold validation while training the RF model [31]. The previously developed RF model allowed classification of 7 aerosol types: pure dust, dust-dominant mixed, pollution-dominant mixed, and pollution aerosols (strongly, moderately, weakly, and non-absorbing) with an accuracy of 56% [31]. The accuracy of the RF model was improved to 72% by merging aerosol classes into 4 types: pure dust (PD), dust-dominated mixed (DDM), and strongly absorbing (SA) and non-absorbing (NA) aerosols [31]. Therefore, in this study, an RF model able to classify 4 aerosol types was utilized to take advantage of its reasonable accuracy. Choi et al. [31] reported that aerosol types classified by the RF model had consistent aerosol optical properties (FMF and spectral dependence of SSA) compared with the typical values for individual aerosol types. More detailed descriptions for development and validation of the model are in Choi et al. [31].

Satellite and AERONET Observation Data
As referred above, the RF model uses 6 satellite input variables including aerosol index, column densities of CO and tropospheric NO 2 , solar zenith angle, land-cover type, and percentage urban area [31]. In this study, we additionally tried to compare the results using other satellite aerosol classification models [26,27]. In Lee et al. [26], aerosol index, AOD, and AE are used as input variables. In Torres et al. [27], aerosol index and column densities of CO are input variables. To collect those satellite data, aerosol index, column densities of CO and tropospheric NO 2 , and solar zenith angle were obtained from TROPOMI measurements onboard the Sentinel-5 Precursor mission. Aerosol index and solar zenith angle were from Sentinel-5P TROPOMI Aerosol Index 1-Orbit L2. Column densities of CO and Tropospheric NO 2 were obtained from Sentinel-5P TROPOMI Carbon Monoxide CO Column 1-Orbit L2 and Sentinel-5P TROPOMI Tropospheric NO2 1-Orbit L2, respectively. AOD was from MODIS measurements onboard Aqua (MYD04_L2). Landcover type and percentage urban area were from a MODIS level-3 product (MCD12C1; Terra+Aqua). In this study, we used the AERONET version 3 level-1.5 data when AOD was >0.4 to minimize uncertainties in AERONET Level 1.5 data. The aerosol optical and microphysical properties were collected, including FMF, SSA, depolarization ratio, effective radius, and volume size distribution. Both satellite input data and AERONET observation data in this Section were collected for the period from 2018 to 2020. Because the overpass time of TROPOMI and MODIS satellites is 13:30 local time, AERONET data for 13:00-14:00 were collected.

Aerosol Properties of the Classified Aerosol Types in Asian Capital Cities
We investigated aerosol properties of the aerosol types detected by the RF model, which included the four types (PD, DDM, NA, and SA) during 2018-2020 in the 11 cities where there is an AERONET site (Table 1). Satellite aerosol types were classified based on the nearest pixel of each AERONET site. The AERONET-derived aerosol properties were collected near overpass time, from 13:00 to 14:00. Figure 2 and Table 2 indicate the AERONET-derived aerosol properties for aerosol types identified by the RF model, including spectral variation of SSA, mean volume size distribution, effective radius, SSA difference (dSSA 440/1020 = SSA at 1020 nm-SSA at 440 nm), and depolarization ratio. The depolarization ratio at 1020 nm was utilized because the AERONET-derived ratio at 1020 nm is reported to have the highest correlation with lidar data [52]. As shown in Figure 2b, volume size distributions of the NA and SA types were similar, with high FMF (0.88 for NA and 0.89 for SA) and low effective radius (0.42 µm for NA and 0.38 µm for SA). Depolarization ratios of the NA and SA types were 0.06 and 0.04, respectively, showing the predominance of spherical fine-mode aerosols detected by the RF model. However, SSA spectral variation for the NA type differs from that of the SA type ( Figure 2a). SSA spectral variability is reported as an indicator of aerosol particle size, type, and source [53][54][55].
In the case of absorbing fine mode aerosols (SA), SSA tends to decrease with increasing wavelength [53][54][55]. Non-absorbing fine mode aerosols (NA) were reported to have no spectral dependence on SSA [53]. In this study, the dSSA 440/1020 value for the NA type was found to be zero, whereas that for the SA type was negative (−0.03). These spectral tendencies and aerosol particle sizes are consistent with previous studies [53][54][55]. Thus, the RF model detects pollution aerosols (both NA and SA) well. However, standard deviations of dSSA 440/1020 for NA and SA were 0.04 and 0.07, respectively. Thus, there was a presence of misclassification of NA and SA. Notably, a higher standard deviation (0.07) of dSSA 440/1020 of SA appeared, showing more misclassification of SA than NA. For PD aerosols, the volume size distribution (Figure 2b) showed a predominance of coarse-mode particles with an FMF of 0.30 and effective radius of 0.70 µm. The fraction of fine-mode particles of PD was higher than that for the DDM, which is aerosol mixed with dust with FMF of 0.58 and effective radius of 0.49 µm. Depolarization ratios, indicators of non-spherical aerosols, were 0.19 for DDM and 0.32 for PD, indicating that PD aerosol detected by the RF model reflects well the characteristics of that type even though only nine such types were detected. As reported in previous studies [53][54][55], the SSA of dustcontaining aerosols (PD and DDM) tend to increase with increasing wavelength with positive dSSA 440/1020 values usually being associated with the coarse mode particle (dust or mixtures of dust). The RF model also detects dust aerosols and aerosols mixed with dust well.

Aerosol Properties of the Classified Aerosol Types in Asian Capital Cities
We investigated aerosol properties of the aerosol types detected by the RF model, which included the four types (PD, DDM, NA, and SA) during 2018-2020 in the 11 cities where there is an AERONET site (Table 1). Satellite aerosol types were classified based on the nearest pixel of each AERONET site. The AERONET-derived aerosol properties were collected near overpass time, from 13:00 to 14:00. Figure 2 and Table 2 indicate the AERONET-derived aerosol properties for aerosol types identified by the RF model, including spectral variation of SSA, mean volume size distribution, effective radius, SSA difference (dSSA440/1020 = SSA at 1020 nm-SSA at 440 nm), and depolarization ratio. The depolarization ratio at 1020 nm was utilized because the AERONET-derived ratio at 1020 nm is reported to have the highest correlation with lidar data [52]. As shown in Figure 2b, volume size distributions of the NA and SA types were similar, with high FMF (0.88 for NA and 0.89 for SA) and low effective radius (0.42 µm for NA and 0.38 µm for SA). Depolarization ratios of the NA and SA types were 0.06 and 0.04, respectively, showing the predominance of spherical fine-mode aerosols detected by the RF model. However, SSA spectral variation for the NA type differs from that of the SA type (Figure 2a). SSA spectral variability is reported as an indicator of aerosol particle size, type, and source [53][54][55]. In the case of absorbing fine mode aerosols (SA), SSA tends to decrease with increasing wavelength [53][54][55]. Non-absorbing fine mode aerosols (NA) were reported to have no spectral dependence on SSA [53]. In this study, the dSSA440/1020 value for the NA type was found to be zero, whereas that for the SA type was negative (−0.03). These spectral tendencies and aerosol particle sizes are consistent with previous studies [53][54][55]. Thus, the RF model detects pollution aerosols (both NA and SA) well. However, standard deviations of dSSA440/1020 for NA and SA were 0.04 and 0.07, respectively. Thus, there was a presence of misclassification of NA and SA. Notably, a higher standard deviation (0.07) of dSSA440/1020 of SA appeared, showing more misclassification of SA than NA. For PD aerosols, the volume size distribution (Figure 2b) showed a predominance of coarse-mode particles with an FMF of 0.30 and effective radius of 0.70 µm. The fraction of fine-mode particles of PD was higher than that for the DDM, which is aerosol mixed with dust with FMF of 0.58 and effective radius of 0.49 µm. Depolarization ratios, indicators of nonspherical aerosols, were 0.19 for DDM and 0.32 for PD, indicating that PD aerosol detected by the RF model reflects well the characteristics of that type even though only nine such types were detected. As reported in previous studies [53][54][55], the SSA of dust-containing aerosols (PD and DDM) tend to increase with increasing wavelength with positive dSSA440/1020 values usually being associated with the coarse mode particle (dust or mixtures of dust). The RF model also detects dust aerosols and aerosols mixed with dust well.

Frequency of Occurrence of Aerosol Types for 3 Years
Over 19 capital cities of Asia, aerosol types were classified for a three-year period using the RF model and compared with those classified using the AERONET-based aerosol classification model [22] and two satellite threshold-based aerosol classification models [27,28]. Figures 3 and 4 indicate the frequency distributions of aerosol-type detections in the capital cities, with the pie charts applying to cases with >200 detections. In the AERONET-based model, aerosol types were classified using AERONET data collected from near real time (13:00-14:00) of the TROPOMI overpass. Of the 11 capital cities with AERONET measurement sites, only Beijing (China; Beijing) and Seoul (Korea; Yonsei_University) had >200 detections. The time interval was extended to 12:00-15:00 to capture data for Dhaka (Bangladesh; Dhaka_University) and Taipei (Taiwan; Taipei_CWB). Figure 3a shows the aerosol-type distribution for the 19 capital cities with populations of >1,000,000. In general, spherical aerosols (NA and SA) had a higher incidence than DDM and PD, except for Kabul, Afghanistan. The NA proportion was high in Singapore (98%) and that of SA was high in Phnom Penh, Cambodia (38%). The highest proportions of PD (31%) and DDM (31%) occurred in Kabul, Afghanistan, which is close to sources of dust aerosols. In East Asia, PD was rarely detected over capital cities except for Ulaanbaatar, Mongolia (12%). However, spherical aerosols (NA and SA) were also detected in Ulaanbaatar as reported in a previous study [56]. DDM was partly detected over East Asia with a range from 3 to 10%, which was affected by the transport of desert dust [57][58][59]. In Southeast Asia, proportions of SA (29%) in capital cities on the Indo-China peninsula (Hanoi, Bangkok, and Phnom Penh) were higher than those (5%) in Jakarta, Kuala Lumpur, Singapore, and Manila. It may be affected by biomass burning events in the Indo-China peninsula [60,61]. South Asia is reported to be affected by dust aerosols transported from West Asia [62] and deserts in South Asia [62,63]. The average proportion of dust-related aerosols (PD and DDM) in South Asia was calculated to be 25%, which is the higher proportion of dust-related aerosols in Southeast Asia (5%) and East Asia (7%).
Thus, the RF model reflects the influence of adjacent sources well. Figure 3b shows the aerosol-type distribution based on AERONET data. For Beijing and Seoul, which were classified using observation data for 13:00-14:00 local time, aerosol-type distributions were similar. Dhaka and Taipei involved AERONET data for 12:00-15:00 local time, and the RF model classified significantly fewer SA-type aerosols. This underestimation of the SA type may be caused by internal uncertainties in the RF model, as reported previously [31], and could be overcome by training the RF model with more data. As more AERONET observation data become available in the satellite observation period, it may be possible to compare and verify aerosol-type classifications at more AERONET stations in the future. However, the RF-based classification can be an alternative tool to classify aerosols in the regions where AERONET-based classification is not possible due to the lack of an AERONET observatory such as Tokyo, Ulaanbaatar, Pyongyang (East Asia), Kabul, Kathmandu, Islamabad (South Asia), Phnom Penh, or Kuala Lumpur (Southeast Asia). Remote Sens. 2021, 13, x FOR PEER REVIEW 7 of 14 Aerosol-type classifications for Asian capital cities from threshold-based satellite models [27,28] are shown in Figure 4 (with >200 classifications). Aerosol types were classified into three types (sulfate, carbonaceous, and desert dust) by Torres et al. [27] in 18 capital cities (Figure 4a) and classified into six types (smoke, sulfate, dust, sea salt and sulfate, dust and smoke, and uncertain) by Lee et al. [26] for Beijing, Seoul, Islamabad, Kabul, and New Delhi (Figure 4b). In general, both satellite models also detected nonabsorbing aerosols with high occurrence rates in all capital cities. There were high contributions of carbonaceous aerosol (strongly absorbing fine mode aerosol) detected by Torres et al. [27] in Pyongyang (49%), Dhaka (44%), and Beijing (46%). Moreover, the contribution of carbonaceous aerosol detected by Torres et al. [27] was found to be similar to the AERONET-based aerosol type (Figure 3b) in Dhaka and Taipei. Torres et al. [27] also Aerosol-type classifications for Asian capital cities from threshold-based satellite models [27,28] are shown in Figure 4 (with >200 classifications). Aerosol types were classified into three types (sulfate, carbonaceous, and desert dust) by Torres et al. [27] in 18 capital cities (Figure 4a) and classified into six types (smoke, sulfate, dust, sea salt and sulfate, dust and smoke, and uncertain) by Lee et al. [26] for Beijing, Seoul, Islamabad, Kabul, and New Delhi (Figure 4b). In general, both satellite models also detected non-absorbing aerosols with high occurrence rates in all capital cities. There were high contributions of carbonaceous aerosol (strongly absorbing fine mode aerosol) detected by Torres et al. [27] in Pyongyang (49%), Dhaka (44%), and Beijing (46%). Moreover, the contribution of carbonaceous aerosol detected by Torres et al. [27] was found to be similar to the AERONET-based aerosol type (Figure 3b) in Dhaka and Taipei. Torres et al. [27] also detected dust aerosol with a maximum occurrence rate (20%) in Kabul among other capital cities while Lee et al. [26] classified smoke (strongly absorbing fine mode aerosol) aerosol in Kabul with a high occurrence rate of 29%.
detected dust aerosol with a maximum occurrence rate (20%) in Kabul among other capital cities while Lee et al. [26] classified smoke (strongly absorbing fine mode aerosol) aerosol in Kabul with a high occurrence rate of 29%.

Seasonal Variations in Aerosol Type
The aerosol-type occurrence rate was calculated seasonally in the capital cities of East Asia, South Asia, and Southeast Asia ( Figure 5). For the seasonal variation, four seasons were labeled as season 1, season 2, season 3, and season 4, respectively, because there are different ways to divide the four seasons: South Asia and Southeast Asia (pre-monsoon, monsoon, post-monsoon, and winter), East Asia (spring, summer, autumn, and winter). As reported in earlier study [22], pollution particles (NA and SA) are predominant in East Asia with the NA type occurring frequently throughout all seasons, and NA being the most prevalent in Hong Kong and Taipei, possibly due to large contributions of maritime

Seasonal Variations in Aerosol Type
The aerosol-type occurrence rate was calculated seasonally in the capital cities of East Asia, South Asia, and Southeast Asia ( Figure 5). For the seasonal variation, four seasons were labeled as season 1, season 2, season 3, and season 4, respectively, because there are different ways to divide the four seasons: South Asia and Southeast Asia (pre-monsoon, monsoon, post-monsoon, and winter), East Asia (spring, summer, autumn, and winter). As reported in earlier study [22], pollution particles (NA and SA) are predominant in East Asia with the NA type occurring frequently throughout all seasons, and NA being the most prevalent in Hong Kong and Taipei, possibly due to large contributions of maritime aerosol, which is strongly light-scattering and highly hygroscopic (Figure 5a). Contributions of the SA type usually peak in season 4 (winter), likely due to increased coal combustion for domestic heating [22]. Over East Asia, the DDM and PD types are usually more prevalent in season 1 (March-May, spring), especially in Ulaanbaatar where the occurrence rate of pure dust and dust mixed with pollution aerosols (PD and DDM) peaked at 32%. Contributions of PD were <3% in capital cities except for Ulaanbaatar (12%) near desert regions. Besides Ulaanbaatar, a higher DDM (20%) type was detected in Beijing in season 1 (March-May, spring) due to the Yellow wind phenomenon in the area [57]. Seasonal aerosol-type occurrence rates in South Asia are shown in Figure 5b. Overall, there were fewer SA and more DDM and PD detections than in East Asia. Dust-containing (DDM and PD) aerosols were usually detected from season 1 to season 2 (March-August, spring-summer). As reported in previous studies, dust aerosols are usually transported during the pre-monsoon (season 1) and South Asian monsoon (season 2) [64,65]. Dust-dominated aerosols were more frequently identified in Kabul than in other capital cities during season 1-season 3 (March-November), with a PD contribution of >50% in season 2 (June-August, South Asian monsoon), likely due to the effect of desert regions [64,65]. From season 3 to season 4 (September-February), contributions of dust aerosol were significantly lower than those of pollution aerosol (SA and NA). In Southeast Asia (Figure 5c) overall, the contribution of NA was predominant, possibly due to maritime aerosols, as these cities are located near the coast. The NA type includes both anthropogenic and marine aerosols [20]. In capital cities of the Indo-china Peninsula (Phnom Penh, Bangkok, and Hanoi), contributions of SA were generally at their maximum in season 1 (March-May, pre-monsoon) and minimum in season 2 (June-August, Southeast Asian monsoon) due to the seasonality of biomass burning activity as reported in previous studies [66,67]. There was negligible seasonality in Jakarta, Kuala Lumpur, Manila, and Singapore with an average proportion of NA of 98%.

Discussion
Aerosol types in 19 Asian capital cities were identified using a previously developed satellite aerosol-classification algorithm [31] and aerosol optical and microphysical prop-

Discussion
Aerosol types in 19 Asian capital cities were identified using a previously developed satellite aerosol-classification algorithm [31] and aerosol optical and microphysical properties investigated for aerosol types detected using an RF model for Asian cities where AERONET sites are located. For cities (e.g., Kabul and Ulaanbaatar) near desert regions, the predominant type was dust aerosol. Dust aerosols were generally mixed with anthropogenic aerosols in many capital cities and detected as the DDM type. In Kabul, the contribution of PD was the highest (>50%) in South Asian monsoon, but the NA type was predominant in most cities. Aerosol types classified by the RF model were compared with those of AERONET and satellite classification models, with those of the RF method generally exhibiting trends similar to those of AERONET-based aerosol types. The other satellite-based classification algorithm [27] seems to classify SA aerosols similarly to AERONET. Seasonal variations in aerosol type detected using the RF model were investigated for capital cities in East Asia, South Asia, and Southeast Asia. NA was the predominant type in most cities and seasons. In East Asia, the dust aerosol (PD and DDM) incidence increased in spring, and the SA type was highest in winter, likely due to increased coal combustion for domestic heating [22]. In particular, a high proportion of DDM was detected in Beijing due to the Yellow wind phenomenon in the area [57]. In South Asia, dust aerosols (PD and DDM) increased between the pre-monsoon and monsoon seasons (March-August) [64,65], while in Southeast Asia NA was generally much more predominant than any seasonal changes.

Conclusions
In this work, aerosol types were classified using the RF model over Asian capital cities for a three-year period. Aerosol properties for each aerosol type were found to be reasonable, showing similar aerosol properties for the typical aerosol type. The globally trained RF model seems to be appropriate for use in Asian cities. This study demonstrated that the RF model can be used to provide aerosol type information with spatially continuous coverage even where there are no AERONET sites. Long-term results of satellite-based aerosol classification methods may be useful in providing climatological input data in the future. Aerosol detections here were relevant to near satellite overpass time (13:30), and if satellite observation data at a range of times were applied to the RF algorithm, it may be possible to extend its temporal coverage. For example, geostationary observations can provide information for trace gases and aerosols with high temporal resolution. Moreover, a combination of satellite-based measurement data, numerical weather prediction models, and chemical transport models should improve the accuracy of the RF model.  Data Availability Statement: For the results and data generated during the study, please contact the first author.