Dust Detection and Intensity Estimation Using Himawari-8/AHI Observation

: In this study, simple dust detection and intensity estimation methods using Himawari-8 Advanced Himawari Imager (AHI) data are developed. Based on the differences of thermal radiation characteristics between dust and other typical objects, brightness temperature difference (BTD) among four channels (BT 11 –BT 12 , BT 8 –BT 11 , and BT 3 –BT 11 ) are used together for dust detection. When considering the thermal radiation variation of dust particles over different land cover types, a dynamic threshold scheme for dust detection is adopted. An enhanced dust intensity index (EDII) is developed based on the reﬂectance of visible/near-infrared bands, BT of thermal-infrared bands, and aerosol optical depth (AOD), and is applied to the detected dust area. The AOD is retrieved using multiple temporal AHI observations by assuming little surface change in a short time period (i.e., 1–2 days) and proved with high accuracy using the Aerosol Robotic Network (AERONET) and cross-compared with MODIS AOD products. The dust detection results agree qualitatively with the dust locations that were revealed by AHI true color images. The results were also compared quantitatively with dust identiﬁcation results from the AERONET AOD and Ångström exponent, achieving a total dust detection accuracy of 84%. A good agreement is obtained between EDII and the visibility data from National Climatic Data Center ground measurements, with a correlation coefﬁcient of 0.81, indicating the effectiveness of EDII in dust monitoring.


Introduction
Mineral dust aerosol is tiny soil particle that is suspended in the atmosphere [1]. It is one of the most important aerosols as it affects the Earth's radiation balance, hydrological cycle, and biogeochemical cycles [1,2]. Previous studies confirm that mineral dust has significant impacts on regional and global climate via a series of complex processes [1][2][3][4][5][6][7]. Dust particles can produce strong backscattering to solar radiation at shortwave band, which significantly changes the albedo of the Earth-atmosphere system. They can also absorb incoming solar radiation and outgoing radiation, and thus alter the distribution of long-wave radiation [8,9]. Besides, dust aerosols can change the microphysical properties of cloud Figure 1 shows the land cover map for the study area, which is adopted from MODIS Land cover product (MCD12C1) [50]. Northwestern China (Taklimakan Desert and Gobi Desert) and Mongolia emit huge amounts of mineral dust every year, causing abundant dust storms, especially during spring seasons.
Remote Sens. 2018, 10, x FOR PEER REVIEW 3 of 19 Algorithms based on the combination of VIS/NIR and IR channels to detect dust are also developed. Miller [15] introduced an index using multi-channels reflectance, and BT at 11 µm and 12 µm to separate dust from other objects. Di et al. [49] developed an dust index, called the Enhanced Dust Index (EDI), which is based on aerosol optical depth (AOD), BTD (3)(4)(5)(6)(7)(8)(9)(10)(11), and the reflectance in 0.63 and 1.6 µm channels, and the dust index shows high consistency with visibility data from ground-based measurements.
In this study, a simple dust detection method using BT over both bright and dark surface is proposed using Himawari-8/AHI observations. An enhanced dust intensity index (EDII), which couples the AOD and spectral reflectance/BT, was also proposed. The AOD was retrieved using multiple temporal AHI VIS channels and validated using AERONET and cross-compared with MODIS AOD. A regression analysis between the visibility and EDII was performed for validation. The algorithm is applied to investigate AHI data over Northern China acquired in March-June 2017. This paper is organized as follows. Section 2 introduces the study area and the data used. Section 3 presents the algorithm, including the theory, dust detection methods, AOD retrieval, and dust index calculation. Section 4 describes the results and validation. A summary is given in Section 5. Figure 1 shows the land cover map for the study area, which is adopted from MODIS Land cover product (MCD12C1) [50]. Northwestern China (Taklimakan Desert and Gobi Desert) and Mongolia emit huge amounts of mineral dust every year, causing abundant dust storms, especially during spring seasons.

Himawari-8 (H8) Advanced Himawari Imager (AHI) Data
Himawari-8 (H8) satellite was launched on 7 October 2014 by the Japan Meteorological Agency (JMA), and started operation on 7 July 2015. H8 carries a new payload, called Advanced Himawari Imager (AHI), which can provide multi-spectral image. As shown in Table 1, AHI has 16 channels, including three visible (VIS) channels (0.47, 0.51, 0.67 µm), three near infrared (NIR) channels, and 10 mid and thermal infrared channels. AHI Level-2 calibrated data from JMA have a spatial Remote Sens. 2018, 10,490 4 of 20 coverage of 120 • by 120 • , and centered at 0 • N, 140 • E with a spatial resolution of 2-10 km for different products. AHI sensor can scan the full-disk every 10 min. It offers unique advantages for dust detection, as well as tracking dust movement. In this paper, AHI Level-2 calibrated data with a spatial resolution of 2 km were used, which have been converted to reflectance for VIS and NIR channels, and to BT for thermal infrared channels. Data are obtained from Japan Aerospace Exploration Agency (JAXA) Earth Observation Research Center (EORC), which are provided by JMA (download from http://www.eorc.jaxa.jp/ptree/terms.html).

MODIS/Aqua and Terra
MODIS aerosol products are used to validate the retrieved AODs. The current Collection 6 (C6) aerosol product provides AOD from both dark-target and deep-blue algorithms and a combination of them for a better spatial coverage [51]. MODIS C6 aerosol product has been intensively evaluated at both the regional and global scale [52,53]. According to Levy et al. [51] the MODIS AOD has expected errors (EE) of ±(0.05 + 0.15τ) for dark-target algorithm and ±(0.05 + 0.2τ) for deep-blue algorithm [54]. In this study, dark-target and deep-blue combined AOD was used for a better spatial coverage.

AERONET Data
AERONET AOD and Ångström exponent were used to validate the retrieved AODs and to define the dust area for AHI dust detection validation. AERONET is a ground-based aerosol network [55], which provides spectral AOD and inversion product derived from the direct and diffuse radiation measured by Cimel sun/sky-radiometers [56]. The Level-1.5 data are cloud screened, and Level-2 data are cloud screened and quality assured. Both Level-1.5 and Level-2 AOD from ten AERONET sites in the study area during March-June in 2017 were used to validate the retrieved AODs from AHI data. AOD and Ångström exponent from both Level-1.5 and Level-2 were also used to define dust area for AHI dust detection results validation.

Visibility Data
Visibility data are used to show the effectiveness of EDII. Visibility data from global surface summary of day product produced by the National Climatic Data Center (NCDC) were used to compare with the dust index. These daily data are built by using the Integrated Surface Data, which includes global data obtained from the USAF Climatology Center, located in the Federal Climate Complex within NCDC. Visibility data from this product are mean visibility for the day. There are hundreds of stations located in China, which provide data since 1970s. In this study, daily mean visibility data during March-June in 2017 were used.

Therory Base of Dust Detection and Intensity Indication
As dust aerosols show distinguished radiative properties over a wide wavelength range from VIS to IR over different underlying backgrounds, it is possible to effectively identify them through Remote Sens. 2018, 10, 490 5 of 20 the distinct spectral reflective and thermal radiation characteristics. Figure 2 shows the comparison of spectral reflective ( Figure 2a) and thermal radiation (Figure 2b) characteristics of several typical objects, including thin dust, thick dust, cloud, bare/desert ground, and dark ground. The spectral reflective and thermal radiation data are derived from AHI sensor and the samples were collected from different dust events in mainland China during 2016 and 2017. The mean reflectance and BT of all the samples were calculated. Dusts show a higher reflectance in the NIR wavelength region and lower reflectance in blue wavelength region when compared with other objects (Figure 2a). The reflectance of dust increases with the wavelength, reaches maximum value at 1.61 µm, and drops slightly at 2.25 µm. The big difference between reflectance at 0.47 and 1.61 µm suggests the potential of utilizing VIS/NIR reflectance to indicate information about dust intensity.
In the thermal-infrared channels, cloud shows low BT, while ground shows high BT at most of channels (Figure 2b). BT at 11.2 µm (BT 11 ) is greater than that at 12.4 µm (BT 12 ) for ground and cloud, while it is lower for dust. BTD between 8.6 and 11 µm (BTD (8-11)) is larger for dust than for the ground and cloud. In addition, the big differences between BT at 3.9 (BT 3 ) and 11 µm for dust particles were evident. Consequently, the dust can be identified by using the BTDs between the different channels [35], e.g., the difference between 3.9 and 11 µm. Besides, Ackerman [36] found that the increasing of dust optical depth increases the BTD (8-11) and BTD (11)(12). Thus, the combination of 3.9, 8.6, 11.2, and 12.4 µm BT is very promising to identify the dust. Furthermore, as a major contribution to aerosol loading, dust aerosols are strongly related to high AODs.

Therory Base of Dust Detection and Intensity Indication
As dust aerosols show distinguished radiative properties over a wide wavelength range from VIS to IR over different underlying backgrounds, it is possible to effectively identify them through the distinct spectral reflective and thermal radiation characteristics. Figure 2 shows the comparison of spectral reflective ( Figure 2a) and thermal radiation (Figure 2b) characteristics of several typical objects, including thin dust, thick dust, cloud, bare/desert ground, and dark ground. The spectral reflective and thermal radiation data are derived from AHI sensor and the samples were collected from different dust events in mainland China during 2016 and 2017. The mean reflectance and BT of all the samples were calculated. Dusts show a higher reflectance in the NIR wavelength region and lower reflectance in blue wavelength region when compared with other objects (Figure 2a). The reflectance of dust increases with the wavelength, reaches maximum value at 1.61 µm, and drops slightly at 2.25 µm. The big difference between reflectance at 0.47 and 1.61 µm suggests the potential of utilizing VIS/NIR reflectance to indicate information about dust intensity.
In the thermal-infrared channels, cloud shows low BT, while ground shows high BT at most of channels ( Figure 2b). BT at 11.2 µm (BT11) is greater than that at 12.4 µm (BT12) for ground and cloud, while it is lower for dust. BTD between 8.6 and 11 µm (BTD (8-11)) is larger for dust than for the ground and cloud. In addition, the big differences between BT at 3.9 (BT3) and 11 µm for dust particles were evident. Consequently, the dust can be identified by using the BTDs between the different channels [35], e.g., the difference between 3.9 and 11 µm. Besides, Ackerman [36] found that the increasing of dust optical depth increases the BTD (8-11) and BTD (11)(12). Thus, the combination of 3.9, 8.6, 11.2, and 12.4 µm BT is very promising to identify the dust. Furthermore, as a major contribution to aerosol loading, dust aerosols are strongly related to high AODs.

Dust Detection
According to Seemann et al. [57], the emissivity in IR channels varied with different surface types, and the thresholds of BTD are strongly related to the source regions [26,58]. It is necessary to consider the regional difference of surface types in dust detection. In order to address this issue, efforts have been made to analyze BT characteristics over different regions. Figure 3 depicts the scatterplots of BTD Remote Sens. 2018, 10, 490 6 of 20 as a function of different surface types, including Gobi (Ground A), relatively dark ground (Ground B), and desert (Ground C). As the plots show, there is no clearly gap between dust-free samples and dusty samples. BTDs under clear sky vary from place to place. The BTD characteristics of dust also vary greatly over typical surface types. It is difficult to distinguish the dust area using a single fixed BTD threshold over different surface types. Instead, if we focus on one surface type at a time, relatively clear boundaries can be found.
To determine a suitable dust identification method for different regions, the spatial distribution of the background BTDs over China mainland was analyzed. Figure 4 shows the clear sky BTDs at 05:00 UTC in May 2017. The clear sky reference BT was obtained by selecting the maximum BT at 11 µm within 30-day window for every hour. The BTD distributions show geographic and surface-type dependency by comparing the background BTDs with the normalized difference vegetation Index (NDVI) and altitude ( Figure 5). On the one hand, the BTDs were strongly related to NDVI. The plots of BTD (3-11) verse BTD (8)(9)(10)(11) show two clusters: (i) when NDVI is greater than 0.3, BTD (8-11) is generally greater than −5 K and BTD (3-11) ranged 1-9; (ii) when NDVI is less than 0.3, the BTD (3-11) tends to be large than 5 K. On the other hand, the BTD (11-12) scatters were also separated into two clusters in the altitude space: (i) when the altitude is less than 3000, the BTD (11-12) is larger than 2.5 K, (ii) when altitude is greater than 3000, the BTD (11-12) is generally greater than 1 and less than 4 K. Thus, a NDVI value of 0.3 and an altitude value of 3000 were combined to classify land surface into three types: (1) arid and semiarid region, (2) relatively dark region, and (3) high altitude region. For these three regions, different thresholds were tested to optimize dust detection. Training data were manually selected based on the AHI true color images. A total of more than 100 scenes of dusty days from March-June in 2016 and 2017 were used for training, specific dates of AHI data used for this training are listed in Table 2. The following dust identification thresholds were determined: For dust detection over arid and semiarid region: For dust detection over relatively dark region: For dust detection over high altitude region: Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 19

Dust Detection
According to Seemann et al. [57], the emissivity in IR channels varied with different surface types, and the thresholds of BTD are strongly related to the source regions [26,58]. It is necessary to consider the regional difference of surface types in dust detection. In order to address this issue, efforts have been made to analyze BT characteristics over different regions. Figure 3 depicts the scatterplots of BTD as a function of different surface types, including Gobi (Ground A), relatively dark ground (Ground B), and desert (Ground C). As the plots show, there is no clearly gap between dust-free samples and dusty samples. BTDs under clear sky vary from place to place. The BTD characteristics of dust also vary greatly over typical surface types. It is difficult to distinguish the dust area using a single fixed BTD threshold over different surface types. Instead, if we focus on one surface type at a time, relatively clear boundaries can be found.
To determine a suitable dust identification method for different regions, the spatial distribution of the background BTDs over China mainland was analyzed. Figure 4 shows the clear sky BTDs at 05:00 UTC in May 2017. The clear sky reference BT was obtained by selecting the maximum BT at 11 µm within 30-day window for every hour. The BTD distributions show geographic and surface-type dependency by comparing the background BTDs with the normalized difference vegetation Index (NDVI) and altitude ( Figure 5). On the one hand, the BTDs were strongly related to NDVI. The plots of BTD (3-11) verse BTD (8)(9)(10)(11) show two clusters: (i) when NDVI is greater than 0.3, BTD (8-11) is generally greater than −5 K and BTD (3-11) ranged 1-9; (ii) when NDVI is less than 0.3, the BTD (3-11) tends to be large than 5 K. On the other hand, the BTD (11-12) scatters were also separated into two clusters in the altitude space: (i) when the altitude is less than 3000, the BTD (11-12) is larger than 2.5 K, (ii) when altitude is greater than 3000, the BTD (11-12) is generally greater than 1 and less than 4 K. Thus, a NDVI value of 0.3 and an altitude value of 3000 were combined to classify land surface into three types: (1) arid and semiarid region, (2) relatively dark region, and (3) high altitude region. For these three regions, different thresholds were tested to optimize dust detection. Training data were manually selected based on the AHI true color images. A total of more than 100 scenes of dusty days from March-June in 2016 and 2017 were used for training, specific dates of AHI data used for this training are listed in Table 2. The following dust identification thresholds were determined: For dust detection over arid and semiarid region: For dust detection over relatively dark region: For dust detection over high altitude region:

Dust Index
Di et al. [49] proposed an Enhanced Dust Index for dust detection over the Xinjiang region using Indian geostationary satellite-INSAT 3A/D IMGER data. The EDI is calculated using observations at VIS, NIR, and TIR channels, and AOD: In which, . and . denote the TOA reflectance at 1.6 and 0.6 µm, while BT3.7 and BT11 refer to the brightness temperature at 3.7 and 11 µm. a = 0.1, b = 10, and c = 0.1 are three coefficients to balance different component contributions. It is proved that this EDI can indicate dust intensity effectively over the Xinjiang region. However, EDI does not work well for northern China in our tests using AHI data, as a single EDI without dust mask might lead to incorrect dust identification due to high AOD caused by fine mode aerosol (e.g., frequently haze weather in northern China).
A new dust index based on the EDI, but adaptive to the AHI sensor data, is designed. Firstly, dust has the minimum and maximum reflectance values at the AHI 0.47 µm and 1.61 µm bands, respectively, based on the analysis in Section 3.1 and in previous studies [59][60][61]. Thus, the AHI 0.47 and 1.61 µm reflectance difference is effective to describe the dust intensity. Similar to the NDDI (( . − . )/ ( . + . ) ) [32], ( . − . )/( . + . ) was employed as a component to describe the dust intensity. Secondly, dust strongly increases the AHI BT at 3.9 µm, while it decreases the AHI BT at 11 µm. The BTD between AHI 3.7 and 11 µm shows a positive relationship with dust intensity. Thirdly, the AOD is also strongly related to the dust intensity, which has been reported by previous studies [53,54]. Consequently, an enhanced dust intensity index (EDII) was proposed for AHI sensor data:

Dust Index
Di et al. [49] proposed an Enhanced Dust Index for dust detection over the Xinjiang region using Indian geostationary satellite-INSAT 3A/D IMGER data. The EDI is calculated using observations at VIS, NIR, and TIR channels, and AOD: In which, R 1.6 and R 0.6 denote the TOA reflectance at 1.6 and 0.6 µm, while BT 3.7 and BT 11 refer to the brightness temperature at 3.7 and 11 µm. a = 0.1, b = 10, and c = 0.1 are three coefficients to balance different component contributions. It is proved that this EDI can indicate dust intensity effectively over the Xinjiang region. However, EDI does not work well for northern China in our tests using AHI data, as a single EDI without dust mask might lead to incorrect dust identification due to high AOD caused by fine mode aerosol (e.g., frequently haze weather in northern China).
A new dust index based on the EDI, but adaptive to the AHI sensor data, is designed. Firstly, dust has the minimum and maximum reflectance values at the AHI 0.47 µm and 1.61 µm bands, respectively, based on the analysis in Section 3.1 and in previous studies [59][60][61]. Thus, the AHI 0.47 and 1.61 µm reflectance difference is effective to describe the dust intensity. Similar to the NDDI ((R 2.1 − R 0.47 )/(R 2.1 + R 0.47 )) [32], (R 1.6 − R 0.47 )/(R 1.6 + R 0.47 ) was employed as a component to describe the dust intensity. Secondly, dust strongly increases the AHI BT at 3.9 µm, while it decreases the AHI BT at 11 µm. The BTD between AHI 3.7 and 11 µm shows a positive relationship with dust intensity. Thirdly, the AOD is also strongly related to the dust intensity, which has been reported by previous studies [53,54]. Consequently, an enhanced dust intensity index (EDII) was proposed for AHI sensor data: where A = 1, B = 10, and C = 0.1 and the AHI AOD retrieval is detailed in the following. These coefficients are determined based on Di et al. [49] and adapted to AHI data by analysis of AHI dust and non-dust observations. Typical areas with dust, deserts, and non-desert land surface are selected as training samples. A bunch of coefficient combinations are tested to make sure the EDII values of dust pixels are greater than 0, and the EDII values of pixels for land surface without dust are less than 0. In this way, the balance among the three component contributions is achieved and the EDII values can also be used for dust detection with 0 values as a boundary value. An aerosol retrieval algorithm assuming a stable surface Bidirectional Reflectance Distribution Function (BRDF) in a short time period (i.e., 1-2 days) was designed to retrieve hourly AOD and surface BRDF parameters from multi-temporal AHI data [62,63]. Unlike the dark target method, this algorithm can retrieve AOD over bright surface and does not make any assumptions on the VIS and NIR reflectance relationships. Unlike the deep blue method, this algorithm does not require historical observation data to construct land surface reflectance. This algorithm is suitable for satellite data with high frequency observations (e.g., 10 min for AHI) so that there are enough observations for a robust model inversion.
The radiative transfer function with coupled BRDF are expressed as [64]: where θ s /θ v is the solar/view zenith angle and ϕ is the relative azimuth angle. ρ TOA (λ, θ s , θ v , ϕ) represents the top of atmosphere (TOA) reflectance, ρ Atm (λ, θ s , θ v , ϕ) is the contribution of the atmospheric path reflectance to the TOA reflectance. T(λ, θ s ) and T(λ, θ v ) are the total transmission of light on the sun-land path and the sensor-land path, respectively. S(λ) is the atmosphere hemispherical albedo. ρ sur (λ) is the surface spectral albedo, while ρ sur (λ, θ s , θ v , ϕ) is the surface bidirectional reflectance. The Ross-Thick-Li-Sparse (RTLS) model was employed to describe the surface BRDF, which expressed as [65] ρ sur (λ, θ s , The three kernel coefficients f i (i = iso, vol, geo) depend on the wavelength, while the kernels K iso , K vol , and K geo depend only on the sun-view geometry. The expressions for all three kernels can be found in Wanner et al. [65]. The RTLS model is also used in MODIS BRDF products [66].
This aerosol retrieval algorithm assumes that the surface bidirectional reflective property (the three BRDF kernel coefficients) change little within short time period (e.g., 1-2 days). Using multi-spectral and multi-temporal AHI observations, the AOD and BRDF kernel coefficients are retrieved through an optimization method, i.e., to minimize the following cost function [67,68]: where x denotes the unknown parameters to be retrieved (i.e., the AOD and three kernel parameters for each band), x 0 is the initial state of the variables, which are derived from the MODIS BRDF product for the kernel coefficients. F(x) is the calculated TOA reflectance based on the radiative transfer forward model (Equation (6)) and variables x. y 0 is the observed AHI TOA reflectance. S y and S a are the error covariance matrices for the AHI TOA reflectance and the initial state of the variables (x 0 ), respectively. The uncertainty of AHI TOA reflectance is band-independent. Consequently, the diagonal components of S y are calculated by multiplying the spectral reflectance with the radiometric calibration accuracy for each band. The diagonal components of S a are obtained from the uncertainty of MODIS BRDF product. The Levenberg-Marquardt method is used for the minimization to obtain the AOD and BRDF kernel parameters. This joint retrieval of AOD and BRDF kernel parameters were proven to be effective [63,68,69].

Results and Validation
The dust detection and dust index algorithms were applied to AHI data acquired from March to June in 2017. Dust detection results are qualitatively compared with the RGB images, and quantitatively compared with the dust detection results from the AERONET measurements. EDII is evaluated using visibility data from the ground-based observations. In addition, we compared the retrieved AHI AOD with MODIS aerosol product and AERONET AOD observations.

Results of Dust Detection and Dust Index
For qualitative evaluation, AHI Red, Green, and Blue (RGB) composite images and EDII images from 02:00 to 08:00 UTC at 2 h interval on 3 May and 4 May 2017 were shown in Figures 6 and 7, respectively. These two days are when a typical Asian dust event over northern China. Note that the EDII values are shown only for the detected dust area and for the area with an EDII greater than 0. The detected dust areas agree qualitatively with the dust locations that are revealed by AHI RGB images. As shown in Figures 6 and 7 Figures 6 and 7. Besides, the dust plume over both bright land surface (e.g., Gobi desert) and 'Dark land' (e.g., North China), and even the small dust plumes over Taklimakan desert, were detected. However, a small part of dust plume is missed, such as the dust storm over Inner Mongolia in Figure 7. This may be because it is mixed with cloud. quantitatively compared with the dust detection results from the AERONET measurements. EDII is evaluated using visibility data from the ground-based observations. In addition, we compared the retrieved AHI AOD with MODIS aerosol product and AERONET AOD observations.

Results of Dust Detection and Dust Index
For qualitative evaluation, AHI Red, Green, and Blue (RGB) composite images and EDII images from 02:00 to 08:00 UTC at 2 h interval on 3 May and 4 May 2017 were shown in Figures 6 and 7, respectively. These two days are when a typical Asian dust event over northern China. Note that the EDII values are shown only for the detected dust area and for the area with an EDII greater than 0. The detected dust areas agree qualitatively with the dust locations that are revealed by AHI RGB images. As shown in Figures 6 and 7 Figures 6 and 7. Besides, the dust plume over both bright land surface (e.g., Gobi desert) and 'Dark land' (e.g., North China), and even the small dust plumes over Taklimakan desert, were detected. However, a small part of dust plume is missed, such as the dust storm over Inner Mongolia in Figure 7. This may be because it is mixed with cloud.    The visual comparison between EDII and AHI true colors image also showed the consistency between the EDII value and dust intensity. This is evident by high EDII over central dust area (e.g., Inner Mongolia) and relatively small EDII over marginal part of dust areas (e.g., Xinjiang). The changes of the high EDII values are in agreement with the movement of the dust storms.
For quantitative evaluation, the dust detection results are compared with the dust identification results from the AERONET AOD and Ångström exponent. Ångström exponent is an indicator of aerosol particle size, i.e., a smaller value indicates larger particles, with value near 0 indicating a coarse mode dominant aerosol [70]. A dust criterion was adapted from Ciren et al. [25], i.e., any AERONET observation with AOD at 1020 nm greater than 0.3 and Ångström exponent between 440 and 870 nm less than 0.6 was classified as dusty observation. In this study, four stations (see Table 2) that are located in northern China covering the dust regions were used.  The visual comparison between EDII and AHI true colors image also showed the consistency between the EDII value and dust intensity. This is evident by high EDII over central dust area (e.g., Inner Mongolia) and relatively small EDII over marginal part of dust areas (e.g., Xinjiang). The changes of the high EDII values are in agreement with the movement of the dust storms.
For quantitative evaluation, the dust detection results are compared with the dust identification results from the AERONET AOD and Ångström exponent. Ångström exponent is an indicator of aerosol particle size, i.e., a smaller value indicates larger particles, with value near 0 indicating a coarse mode dominant aerosol [70]. A dust criterion was adapted from Ciren et al. [25], i.e., any AERONET observation with AOD at 1020 nm greater than 0.3 and Ångström exponent between 440 and 870 nm less than 0.6 was classified as dusty observation. In this study, four stations (see Table 2) that are located in northern China covering the dust regions were used. The visual comparison between EDII and AHI true colors image also showed the consistency between the EDII value and dust intensity. This is evident by high EDII over central dust area (e.g., Inner Mongolia) and relatively small EDII over marginal part of dust areas (e.g., Xinjiang). The changes of the high EDII values are in agreement with the movement of the dust storms.
For quantitative evaluation, the dust detection results are compared with the dust identification results from the AERONET AOD and Ångström exponent. Ångström exponent is an indicator of aerosol particle size, i.e., a smaller value indicates larger particles, with value near 0 indicating a coarse mode dominant aerosol [70]. A dust criterion was adapted from Ciren et al. [25], i.e., any AERONET observation with AOD at 1020 nm greater than 0.3 and Ångström exponent between 440 and 870 nm less than 0.6 was classified as dusty observation. In this study, four stations (see Table 2) that are located in northern China covering the dust regions were used.
Dust detection results using hourly AHI observations from 01:00 UTC to 10:00 UTC from March to June in 2017 are evaluated. The AERONET observations are time-space matched with the dust detection results from AHI data. AERONET measurements were required within 15 min of satellite scan time. Only the satellite results within 25 km distance from AERONET station are used. For satellite dust detections, if more than one half of the cloud-free pixels were identified as dusty, the area is considered as dusty.
The performance of the dust detection was evaluated using three statistical indices [25], including the Accuracy, Probability of correct detection (PCD), and Probability of false detection (PFD):  Table 3 gives the evaluation results. A total of 172 matched points were obtained, with an Accuracy of 84%, a PCD of 77%, and a low PFD of 9%. Note that only the AERONET observations with AOD greater than 0.3 were used in matching points to balance the dust and non-dust matched points, i.e., to avoid an extremely large number of non-dust observations. Station-based statistics shows large differences. The best performance occurred at AOE_Baotou station, which has accuracy, PCD, and PFD values of 81%, 86%, and 13%, respectively. The worst detection was found at Beijing station, with a relatively low PCD of 69%. This is because Beijing is affected by both the transported dust particle and the haze that is caused by industrial pollution and traffic emission. The mixed particles exhibit more complex properties making dust detection difficult. Xianghe station shows similar behavior with Beijing. The Dalanzadgad station has low Accuracy and high PFD due to the relatively low number (25) of matchups and high 'DN' number. Table 3. Statistics metrics for comparison of AHI-based dust identification and AERONET observations. Many studies have shown that the visibility is a good indicator of dust concentration [71,72]. In order to demonstrate that the EDII values are able to describe the intensity of dust, a regression analysis between the obtained EDII and ground measurements of visibility was performed. When considering that the low visibility might result from the high aerosol loading caused by fine mode particles, only data from stations located in northwest China were used to reduce the influence of the fine mode aerosols. AHI detection results within 25 km from the station are used. Daily average of the AHI EDII results was calculated and only AHI detection results with at least three cloud-free observations in one day are used. The data used for this comparison is from March-June 2017. A total of 116 matchups from 31 stations were obtained. Figure 9 depicts the scatter plot and the regression curve. A good agreement was obtained between the visibility and EDII, a quantitative relationship was obtained as: y = 1.67· exp(−0. 25

Results of Aerosol Optical Depth Retrieval
Retrieved AODs were compared with the MODIS C6 and AERONET AOD data. Figure 10 shows the comparison between the MODIS C6 AOD with our retrieved AOD at 03:00 UTC and 06:00 UTC. AODs from JMA L2 aerosol product were also employed to compare with our retrieval, as shown in the right column. Their corresponding RGB composite images were shown in Figure 7. Our retrieved AODs show similar spatial distribution with that of the MODIS C6 AOD, while the JMA L2 AODs missed some retrieval over strong dust area. Both MODIS and our retrieved AODs reveal the sharp contrast between dusty and clear areas.
To assess the reliability of AHI AOD retrievals, both AERONET and MODIS AODs are used to quantitatively compare with our retrieved AODs. The AHI AOD and AERONET AOD were compared for a period from March to June 2017 over the study area for the ten AERONET sites. For validation, AERONET AODs were converted to 550 nm using Angstrom exponent. Figure 11a shows the comparison of AHI AOD and AERONET AOD. These two datasets show good consistency, with R 2 = 0.79. Figure 11b is the cross comparison between AODs from MODIS and our AHI retrieval on 4 May 2017. It has a R 2 = 0.81, showing good agreement.

Results of Aerosol Optical Depth Retrieval
Retrieved AODs were compared with the MODIS C6 and AERONET AOD data. Figure 10 shows the comparison between the MODIS C6 AOD with our retrieved AOD at 03:00 UTC and 06:00 UTC. AODs from JMA L2 aerosol product were also employed to compare with our retrieval, as shown in the right column. Their corresponding RGB composite images were shown in Figure 7. Our retrieved AODs show similar spatial distribution with that of the MODIS C6 AOD, while the JMA L2 AODs missed some retrieval over strong dust area. Both MODIS and our retrieved AODs reveal the sharp contrast between dusty and clear areas.

Results of Aerosol Optical Depth Retrieval
Retrieved AODs were compared with the MODIS C6 and AERONET AOD data. Figure 10 shows the comparison between the MODIS C6 AOD with our retrieved AOD at 03:00 UTC and 06:00 UTC. AODs from JMA L2 aerosol product were also employed to compare with our retrieval, as shown in the right column. Their corresponding RGB composite images were shown in Figure 7. Our retrieved AODs show similar spatial distribution with that of the MODIS C6 AOD, while the JMA L2 AODs missed some retrieval over strong dust area. Both MODIS and our retrieved AODs reveal the sharp contrast between dusty and clear areas.
To assess the reliability of AHI AOD retrievals, both AERONET and MODIS AODs are used to quantitatively compare with our retrieved AODs. The AHI AOD and AERONET AOD were compared for a period from March to June 2017 over the study area for the ten AERONET sites. For validation, AERONET AODs were converted to 550 nm using Angstrom exponent. Figure 11a shows the comparison of AHI AOD and AERONET AOD. These two datasets show good consistency, with R 2 = 0.79. Figure 11b is the cross comparison between AODs from MODIS and our AHI retrieval on 4 May 2017. It has a R 2 = 0.81, showing good agreement.

Conclusions
In this study, we presented a dust detection method and a dust index to estimate the dust intensity over China mainland using AHI data. Based on a dynamic threshold of BTD among four channels: BTD (11)(12), BTD (3)(4)(5)(6)(7)(8)(9)(10)(11), and BTD , the dust detection method can achieve fast dust identification for the different surface types over northern China. By combining the spectral reflectance at VIS/NIR, BT at TIR, and AOD, the enhanced dust intensity index (EDII) was developed with AOD that was estimated from multiple temporal AHI observations.
The dust detection and density estimation (EDII) algorithms were applied on the AHI data acquired in March-June 2017 in mainland China. The AHI dust detection results and intensity map were examined by the analysis of a strong dust event over northern China during 3-4 May 2017 using the AHI true color images and the air mass trajectories data. It was shown that the dust intensity map well captured the dust area and its movement. This dust detection method is simple, fast, and ideally suited for monitoring dust in real time using geostationary satellite H8 data. Unlike other dust detection methods using a universal threshold value on BTDs, the dynamic threshold scheme with threshold values adaptive to different land surfaces allowed for more accurate dust identification under various surface conditions. The AERONET measurement derived dust detection result was used to validate the AHI dust detection result, showing a dust detection agreement of 84%. In addition to identifying the dust area, the EDII can be used as a semi-quantitative parameter to describe the dust intensity and showed good agreement with the visibility data (a correlation coefficient of 0.81). In addition, to gain more confidence on the EDII, the retrieved AODs used for To assess the reliability of AHI AOD retrievals, both AERONET and MODIS AODs are used to quantitatively compare with our retrieved AODs. The AHI AOD and AERONET AOD were compared for a period from March to June 2017 over the study area for the ten AERONET sites. For validation, AERONET AODs were converted to 550 nm using Angstrom exponent. Figure 11a shows the comparison of AHI AOD and AERONET AOD. These two datasets show good consistency, with R 2 = 0.79. Figure 11b is the cross comparison between AODs from MODIS and our AHI retrieval on 4 May 2017. It has a R 2 = 0.81, showing good agreement.

Conclusions
In this study, we presented a dust detection method and a dust index to estimate the dust intensity over China mainland using AHI data. Based on a dynamic threshold of BTD among four channels: BTD (11)(12), BTD (3)(4)(5)(6)(7)(8)(9)(10)(11), and BTD , the dust detection method can achieve fast dust identification for the different surface types over northern China. By combining the spectral reflectance at VIS/NIR, BT at TIR, and AOD, the enhanced dust intensity index (EDII) was developed with AOD that was estimated from multiple temporal AHI observations.
The dust detection and density estimation (EDII) algorithms were applied on the AHI data acquired in March-June 2017 in mainland China. The AHI dust detection results and intensity map were examined by the analysis of a strong dust event over northern China during 3-4 May 2017 using the AHI true color images and the air mass trajectories data. It was shown that the dust intensity map well captured the dust area and its movement. This dust detection method is simple, fast, and ideally suited for monitoring dust in real time using geostationary satellite H8 data. Unlike other dust detection methods using a universal threshold value on BTDs, the dynamic threshold scheme with threshold values adaptive to different land surfaces allowed for more accurate dust identification

Conclusions
In this study, we presented a dust detection method and a dust index to estimate the dust intensity over China mainland using AHI data. Based on a dynamic threshold of BTD among four channels: BTD (11)(12), BTD (3)(4)(5)(6)(7)(8)(9)(10)(11), and BTD , the dust detection method can achieve fast dust identification for the different surface types over northern China. By combining the spectral reflectance at VIS/NIR, BT at TIR, and AOD, the enhanced dust intensity index (EDII) was developed with AOD that was estimated from multiple temporal AHI observations.
The dust detection and density estimation (EDII) algorithms were applied on the AHI data acquired in March-June 2017 in mainland China. The AHI dust detection results and intensity map were examined by the analysis of a strong dust event over northern China during 3-4 May 2017 using the AHI true color images and the air mass trajectories data. It was shown that the dust intensity map well captured the dust area and its movement. This dust detection method is simple, fast, and ideally suited for monitoring dust in real time using geostationary satellite H8 data. Unlike other dust detection methods using a universal threshold value on BTDs, the dynamic threshold scheme with threshold values adaptive to different land surfaces allowed for more accurate dust identification under various surface conditions. The AERONET measurement derived dust detection result was used to validate the AHI dust detection result, showing a dust detection agreement of 84%. In addition to identifying the dust area, the EDII can be used as a semi-quantitative parameter to describe the dust intensity and showed good agreement with the visibility data (a correlation coefficient of 0.81). In addition, to gain more confidence on the EDII, the retrieved AODs used for EDII calculation were validated using AERONET and MODIS AODs, both showing high correlation coefficient values.
Future work will focus on the dust detection at night time to achieve a continuous monitoring of dust. Other than the algorithm difference for night time dust detection, we noted that the dust validation is also challenging especially considering that there is no 'ground truth' data for dust quantification at either day or night time. Visibility data from ground measurement used in this study can be related to satellite derived dust intensity, assuming that there is no horizontal and vertical dust distribution difference. Instead, ground aerosol measurements from active remote sensing with optical aerosol characterization could be utilized for dust detection validation. Lidar and ceilometer can provide observations during both day and night time, can also provide vertical information of dust layer, and can distinguish dust layers with other layers [73][74][75]. Future study would consider using Lidar and ceilometer data for validation subject to the ground station availability. The new generation of geostationary satellite, like FY-4 equipped with a Multi-channel scanning radiometer, which is similar to AHI, can provide better sight for the China and Asia area. This presented method can be also adapted for FY-4 dust detection, to achieve dust detection for China and West Asia with more complete coverage.