Seasonal Variations and Long-term Trend of Mineral Dust Aerosols over The Taiwan Region

Atmospheric dust aerosols are known to affect the air quality and public health as well as climate and weather systems. An increasing number of modeling studies have related ice nucleation with the number concentrations of dust particles with a diameter larger than 500 nm (ND,d>500nm). In this paper, the seasonal variation, vertical properties, and long-term trend of ND,d>500nm over the Taiwan region are analyzed, using simulations from a global chemical transport model with size-resolved particle microphysics. Over Taiwan, ND,d>500nm shows a bimodal seasonal variation distribution with two peaks in spring (March–May) and fall-early winter (October–December). In the different seasons, ND,d>500nm varies by about one order of magnitude from summer to spring (0.06–1.23 cm in the boundary layer, 0.03–0.55 cm in the middle and lower troposphere, and 0.006–0.03 cm in the upper level). Vertically, ND,d>500nm profiles show the unimodal distribution, with the highest ND,d>500nm appears at ~1 km and decreasing with altitude. From surface to high levels, the frequencies of intense dust events decrease in fall (September–November) and increase in summer months (June–August). The long-term model results suggest a decreasing trend of the strong dust event frequencies and annual mean ND,d>500nm over Taiwan in the last two decades. From 1999 to 2018, the number of strong dust event days and ND,d>500nm decreased by 40–43% and 37–54%, respectively, under 4 km, and the decline is weaker at higher altitudes. The analysis suggests that these decrease trends are caused by the declining Asian dust emissions.


INTRODUCTION
Mineral dust is indicated to make a significant contribution to global atmospheric aerosol loading, for its large emission rate (up to 5000 Tg yr -1 ) (Engelstaedter et al., 2006) and long-range transportability Uno et al., 2009). Dust storms are observed to influence air quality and public health worldwide (Ganor et al., 2009). Mineral dust particles have significant impacts on climate and weather systems by altering the radiation forcing (Tegen et al., 1997;Liao and Seinfeld, 1998), cloud nucleation (Levin et al., 1996;DeMott et al., 2003), and convective invigoration (Andreae et al., 2004). The indirect influences of dust on clouds have also been revealed by a series of observations and numerical studies (Tao et al., 2012;Liu et al., 2012a, b;Fan et al., 2016;Kanji et al., 2017). Dust particles are recognized as one of the most important ice-nucleating particles (INPs), by influencing the heterogeneous freezing in the mixed-phase clouds (T > −37°C) Murray et al., 2012). Previous observation studies reveal that the long-range transported dust particles play a critical role in the development of cloud and precipitation (Ault et al., 2011;Creamean et al., 2013). Modeling studies indicate that the presence of mineral dust leads to the initiation of mixed-phase clouds and increases precipitation efficiency (Muhlbauer and Lohmann, 2009;Fan et al., 2014). Our recent study based on a long-term analysis of multiple datasets suggests a significant correlation between atmospheric dust loading and the summertime precipitation over the mountainous region in Taiwan (Zhang et al., 2020). Thus, an increasing number of numerical models have explicitly considered the dependence of ice nucleation rate on dust number concentrations Thompson and Eidhammer, 2013;DeMott et al., 2015).
Over the Taiwan region, dust storms are indicated by numerous studies to have significant impacts on air quality, visibility, and public health over the Taiwan area (Chen et al., 2005;Cheng et al., 2005;Liu et al., 2006;Chiu et al., 2008). Located in East Asia, Taiwan is highly influenced by the long-range transport of mineral dust aerosols from mainland China, the Middle East, and the Sahara (Chen et al., 2004;Lin et al., 2007;Liu et al., 2009;Chang et al., 2010;Lin et al., 2012;Shahsavani et al., 2012;Hsu et al., 2012). Recent research also suggests the wintertime river-dust as a local source of dust aerosols . Series of studies on atmospheric dust aerosols have been conducted in Taiwan, providing a profound understanding of mineral dust properties over this region (Liu and Shiu, 2001;Hsu et al., 2008;Lin et al., 2012;Tsai et al., 2014;Provençal et al., 2017). However, restricted by the observation technology, there are few studies on the long-term 3-dimensional dust properties in the area.
In this context, we use a global chemical transport model with size-resolved particle microphysics to study the mineral dust aerosols over the Taiwan region from 1999 to 2018. The main objective of the present work is to study the temporal and spatial characteristics of dust aerosols, including the typical value of dust number concentrations, seasonal variation, vertical distribution, and the long-term trends, that could provide a profound understanding of mineral dust aerosols and help to improve the dust-ice nucleation parameterizations over Taiwan.

Model and Data
(1) GEOS-Chem-APM model and dust simulations: Limited by the lack of long-term 3-D quantitative observations, dust mass, and number concentration simulations are used in our research. The dust simulation is based on a global 3-D chemical transport model (GEOS-Chem) with size-resolved Advanced Particle Microphysics (APM) model incorporated (Yu and Luo, 2009). GEOS-Chem-APM model is driven by Global Modeling and Assimilation Office (GMAO) Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorology fields. The GEOS-Chem-APM model was run globally at 2° × 2.5° horizontal resolution and 47 vertical layers (with 14 layers from the surface to ~2 km), from 1999 to 2018. In the GEOS-Chem-APM, dust emission is calculated with the Mineral Dust Entrainment and Deposition (DEAD) scheme (Zender et al., 2003), and dust aerosols with diameters from 30 nm to 25 µm are represented by 15 bins. In previous studies, GEOS-Chem has been widely validated and used for dust simulations (Prenni et al., 2009;Xu et al., 2017;Zhang et al., 2019). Hourly Dust simulations by GEOS-Chem-APM are outputted at two grids of LLN (Lulin) site (23.47°N, 120°E) and TP (Taipei) site (25.09°N, 121.56°E) over the Taiwan region. Previous measurements indicate that the heterogeneous freezing of cloud ice particles is closely related to the dust number concentration with a diameter larger than 500 nm (ND,d>500nm) Creamean et al., 2013). Thus ND,d>500nm (procced into local time) is used to represent the strength of daily dust aerosol loading at each level over Taiwan. Mass concentration of dust smaller than 10 µm (PM10) is also output to compare to the observations to validate the model simulation.
(2) CALIPSO dust observation: Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) Lidar Level 2 Vertical Feature Mask Data (https://www-calipso.larc.nasa.gov) provided by the NASA Langley Research Center's (LaRC) ASDC DAAC and NASA Earth Science Data and Information System (ESDIS) project are used in this study (Winker, 2018). CALIPSO observations from 2006 to 2018 with satellite track passing through Taiwan (21.7°-25.5°N, 119.8°-122.2°E) are used to represent the atmospheric dust loading in this area. CALIPSO level 2 classified aerosol data provides information on the vertical properties of dust, polluted dust, biomass burning, polluted continental, clean continental, and clean marine aerosols. In this study, the classified "dust" and "polluted dust" pixels are treated as dust signals. As there is no quantitative dust concentration in the CALIPSO data, the dust ratio is used as the qualitative reference, which is defined as the ratios of dust pixels to all pixels within 23°-25°N under 10 km. Dust ratios are used to compare with and validate the GEOS-Chem ND,d>500nm simulations. The CALIPSO observations with more than 20% missing data in this region are excluded.
(3) PM10. Aerosol Speciation: Long-term in-situ aerosol sampling (2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017) was conducted at the Cape Fuguei Research Station (25.30°N, 121.54°E, 10 m) at the northern tip of Taiwan Island (Chou et al., 2017). Mass concentrations of PM10 aerosols, organic carbon (OC), elemental carbon (EC), Na + , NH4 + , K + , Mg 2+ , Ca 2+ , Cl -, NO3 -, SO4 2were observed at this station. As there is no direct observation of mineral dust, Calcium ion concentration is used as a proxy. Ca is considered as one of the basic elements composing mineral dust (Perry et al., 1997;Formenti et al., 2008;Remoundaki et al., 2011), by following the method of Song and Carmichael (2001), we assume a calcium/dust ratio of 6.8% to estimate the mass loading of dust particles (Song and Carmichael, 2001). The derived PM10 dust mass concentrations are compared with the PM10 dust concentrations simulated by the GEOS-Chem-APM model at the surface level in the same region.  Examples of (b) dusty and (c) clean days observed by the CALIPSO satellite (ND = "not determined", CM = "clean marine", D = "Dust", PC = "Polluted continental", CC = "clean continental", PD = "Polluted dust", S = "smoke" O = "other"); and (d) CALIPSO observed dust pixel ratios (red spots) and model ND,d>500nm simulations (blue columns) in the selected days, and grey shades are those days CALIPSO has over 20% missing data (not used in the comparison). obvious difference, indicating that CALISPO observations can capture the dust signals be used to identify the occurrence of dust events. Fig. 1(d) shows the comparisons between CALIPSO dust ratios and GEOS-Chem ND,d>500nm simulations for all 187 days. The result suggests that ND,d>500nm simulations are generally consistent with the CALIPSO observations (r = 0.55) and the GEOS-Chem model can simulate the temporal variation of atmospheric dust loading and observed strong dust events over the Taiwan area. However, as the quantitative dust concentration is not available in the CALIPSO dataset, dust ratios may not be proportional to the dust number concentration. In Table 1, 26 strong dust events observed by CALIPSO (with dust ratio > 10%) are selected to qualitatively compare with the GEOS-Chem dust simulations on the corresponding date. The result shows that, in the 26 days with CALIPSO dust signals, GEOS-Chem simulated 25 dust events (top 25%, ND,d>500nm > 0.1 cm -3 ) (96%) and 21 strong dust events (top 10%, ND,d>500nm > 0.2 cm -3 ) (81%). The result suggests that the GEOS-Chem model can simulate most of the intermittent strong dust events, which is one of the major sources of atmospheric dust over the Taiwan region (Chen et al., 2004;Hsu et al., 2012;Lin et al., 2012;Chou et al., 2017).

Validation of Model Dust Simulations
To validate the long-term GEOS-Chem-APM dust simulation, we compare the surface PM10 dust mass concentration simulations with the in-situ dust observations. To filter out some background noises, observation samplings with the lower 20% values are eliminated, by following the approach of Hand et al. (2016) and Achakulwisut et al. (2017). Fig. 2 shows the comparison between the time series of the observed and simulated dust mass concentrations that the temporal variation of model dust simulation generally matches with the observation, and the model result can simulate the observed strong dust events. In the 1928 observation days (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017), the model dust simulation shows a high correlation with the observation (r = 0.69). The mean model dust concentration is 8.73 µg m -3 , close to the mean value of dust observation (7.32 µg m -3 ). Fig. 2 also reveals that, to some extent, the GEOS-Chem model underestimated the dust concentration during weak dust days and overestimated the strong dust events. This  underestimation in weak dust event days may be caused by (1) the Calcium ion mass concentration of non-mineral dust is involved in observation; (2) in the GEOS-Chem model, the local dust emission over Taiwan is underestimated, such as the wintertime river-dust emission in Taiwan ; and (3) the uncertainties associated to the mass mixing ratio of calcium in dust particles. Nevertheless, the overall well agreement between dust observation and simulation suggests that the GEOS-Chem-APM model can simulate the mineral dust temporal variation and the occurrence of strong dust events over the Taiwan region.

The Temporal and Spatial Properties of Dust Aerosols over Taiwan
As mentioned in the Model and data Section, daily (local time, converted from hourly data) ND,d>500nm simulations are output at two sites of LLN (23.47°N, 120°E) and Taipei (25.09°N, 121.56°E), covering the whole of Taiwan region. ND,d>500nm simulations at the two grids are averaged to represent the atmospheric dust loading over Taiwan. Fig. 3 shows the 20-years (1999-2018) ND,d>500nm simulations at different altitudes over the Taiwan region. Fig. 3 exhibits that dust transport events have a significant influence on atmospheric dust loading in Taiwan, with ND,d>500nm varying by 2-3 orders of magnitude from clean to adjacent dust events days in the same seasons. The result also indicates significant annual and seasonal variations of ND,d>500nm: annually, the maximum ND,d >500nm in different years can vary by about ten times (from 38.4 cm -3 in 2010 to 6.8 cm -3 in 2009); seasonally, the frequencies of strong dust event in spring and fall are much higher than that in summer and winter. At most times, ND,d >500nm at the surface (the black line), 1 km (blue), and 2 km (orange) are much stronger than that at higher levels, indicating that dust transporting events over Taiwan primarily occur at low latitudes. Below we focus on the typical value, seasonal variation, and long-term trend of dust aerosols over the Taiwan area.  0.03 cm -3 in the upper level). In a year, the monthly mean ND,d>500nm has the highest value in April (~1.77 cm -3 ), and ND,d >500nm is the lowest in July (~0.03 cm -3 ) which is two orders of magnitude smaller than that in April. Fig. 4 also suggests a unimodal vertical distribution of atmospheric dust aerosols over Taiwan, with monthly mean ND,d>500nm having the peak value at lower levels (surface~1 km).

Seasonal Variations and Vertical Properties of Mineral Dust
To further study the typical values and vertical distributions of mineral dust over Taiwan, we choose the lower 10% and lower quartile of ND,d>500nm at each level to represent the dust number concentrations in weak dust days in specific seasons, and the upper quartile and upper 10% to represent those in strong dust days. The statistical results are essential to providing a profound understanding of atmospheric dust aerosol and can be used in numerical studies for heterogeneous ice nucleation calculation. Fig. 5 shows that, in the four seasons, dust concentrations in both clean and dusty days have similar distribution properties, with peak ND,d>500nm appearing at low altitude, and decreasing with height. In spring, fall, and winter, ND,d>500nm decreases by about two orders of magnitude from low to high elevations; in summer, ND,d>500nm varies by roughly an order of magnitude, less dramatically than in other seasons. Fig. 5 shows the significant distinctions between clean and dusty days, with ND,d>500nm varying by 2-3 orders of magnitude at each level in the same season. On clean days, dust number concentrations are lower than 0.1 cm -3 , which have limited influences on air quality and cloud development. In dust days, ND,d >500nm can reach higher than 3 cm -3 in the boundary layer, which could influence that air quality and public health (Uematsu et al., 1983;Li et al., 1996;Cheng et al., 2005;Liu et al., 2006;Chiu et al., 2008); and the ND,d>500nm can reach about 0.2-2 cm -3 in the middle and high levels. Previous laboratory experiments suggest that the activation ratio of mineral dust is about 0.5-3% at ~−20°C (Zimmermann et al., 2008;Niemand et al., 2012), the dust-contributed INP number concentration can reach about 1-60 L −1 over Taiwan, high enough to substantially influence the development of cloud and precipitation (Creamean et al., 2013;Fan et al., 2014).
Previous studies indicate that the atmospheric dust loading is highly influenced by the intermittent long-range transport of dust from mainland China, the Middle East, and Sahara (Chen et al., 2003;Hsu et al., 2012;Lin et al., 2012;Chou et al., 2017). We choose the mean values of ND,d>500nm during the top 25%, 10%, and 1% dust events at the six levels for four seasons during 1999-2018, to represent the strength of dust events, strong dust events, and extreme dust events, respectively. The mean and maximum ND,d>500nm at the six levels in four seasons are also listed. The statistical results are given in Tables 2-5. Table 2 shows that in spring (MAM), the mean value varies from 0.012 cm -3 (8 km) to 1.281 cm -3 (1 km); from the mean state, ND,d>500nm can increase by ~3, 5, 17 times during dust events, strong and extreme events, respectively. In summer (JJA) (Table 3), the dust has the lowest number concentration over the year, with the mean ND,d>500nm of about 0.005 (8 km) ~0.062 cm -3 (1 km) although ND,d>500nm during extreme dust events can be a factor of ~30 higher. The result suggests that although the mean ND,d>500nm is lower than in other seasons, the transport of dust in summer has more significant contributions to the atmospheric dust loading over Taiwan. In fall (SON) (Table 3), vertically, the mean values at different heights vary from 0.471 (surface) to 0.002 cm -3 (8 km), and from the mean value, ND,d>500nm increases by ~3, 6, 20 times when dust events, strong and extreme events occur. In winter (Table 5), ND,d>500nm at low altitudes (surface-2 km) is much higher than that in summer, but lower at higher altitudes (6-8 km). In Fig. 5 and Tables 2-5, the results suggest the shifting vertical properties of dust particles, that in spring and summer, the highest ND,d>500nm appears at Mean ND,d>500nm during dust events (at the upper 25% of overall 7305 days), strong dust events (at the upper 10%), and extreme events (at the upper 1%) as well as maximum ND,d>500nm at 6 height model levels in spring (MAM). The peak value in each dust event category is marked with shadow.   a higher altitude than in fall and winter. The shift of the dust sources and transport modes could be important factors of this seasonal variation. It is also notable that, although ND,d>500nm at 4-8 km is much lower than that at lower levels, middle-high level ND,d>500nm can reach more than 0.1 cm -3 and 0.4 cm -3 when dust events and extreme events happen. Influenced by the dynamic uplifting of terrain and convection systems, ND,d>500nm at the middle-high level could be even higher. Due to the low atmospheric temperatures at the middle and upper troposphere, dust particles could impact the cloud development through heterogeneous freezing and convective invigoration (DeMott et al., 2003;Andreae et al., 2004;Creamean et al., 2013;Fan et al., 2014).
Considering the importance of dust events over Taiwan, the strong dust events are studied in aspects of both frequency and strength, which are factors influencing the atmospheric dust aerosol loading. Strong dust events (top 10%) usually last for a couple of days, in this research a single day is defined as "an event day", and the number of days is used to represent the frequency of dust events. Fig. 6 shows the frequencies (columns) and ND,d>500nm (dashed lines) of strong dust events at six selected altitudes. The GEOS-Chem-APM model results (Figs. 6(a), 6(b), and 6(c)) show the bimodal seasonal distributions of strong dust events at low altitude (surface, 1 km, and 2 km), with high frequencies and concentrations mainly occur in spring (March, April, and May), and fallwinter (October, November, and December). At surface and 1 km levels, there are more strong dust events in spring, about twice that in fall and winter. However, the dust event strengths are close in spring (with highest ND,d>500nm = 4.29 cm -3 at the surface and 4.76 cm -3 at 1 km occur in April) and fall (with highest ND,d>500nm = 4.10 cm -3 at the surface and 3.89 cm -3 at 1 km occur in October). The result reveals that, although the frequency of dust events in fall and winter is lower than in spring, when dust storm happens, the atmospheric dust loading can be very high, influencing the environment and air quality over Taiwan. Fig. 6 shows that from the surface to high altitudes, the frequencies of strong dust storm events gradually increase in summer and early fall (JJAS) and decrease in late fall and winter (OND). Strong dust events at the middle troposphere (4, 6, and 8 km) show the unimodal seasonal variation; this pattern could be clearly distinguished from the boundary layer. The GEOS-Chem-APM model results exhibit different seasonal patterns of strong dust events at low (surface, 1 and 2 km) and middle-upper troposphere (4, 6, and 8 km), which may be caused by the various sources and heights of dust transport. The long-range transport of dust from east Asia at low altitudes during late winter and spring has been investigated by previous observation and modeling studies (Chen et al., 2004;Lin et al., 2012). In the mid-to-upper troposphere over Taiwan, dust might have a strong contribution from the Middle East and Sahara . A recent study also indicates the river-dust event as a local source of dust during the wintertime .

The Long-term Trend of Dust Aerosols over Taiwan
The above results show that the properties of dust events vary at different vertical levels. Accordingly, the dust event frequency and strength are used to study the long-term trend of atmospheric dust aerosols at each level over the Taiwan region. The number of dust event days (NOE) with the top 10% ND,d>500nm of overall 7305 days (1999-2018) is used to represent the frequency of strong dust events, and the annual mean ND,d>500nm represents the atmospheric dust loading.
The decline of ND,d>500nm over Taiwan is significant in the past two decades and it is worthwhile to analyze the underlying drivers for the decrease. Various reasons could contribute to this trend of dust loading over Taiwan, and we think one of the most important influencing factors is the Fig. 7. Long-term (1999Long-term ( -2018 trends of annual mean ND,d>500nm (solid line) and the dust numbers of strong dust event day (columns) at six levels. Broken lines are the linear trends of NOE (black) and ND,d>500nm (red).
decreasing of East Asia dust emission. Previous studies indicate that the dust emission and outbreaks over East Asia show a notable downtrend over the past decades (Wang et al., 2008;Kurosaki et al., 2011;Wu et al., 2018). The trend has been suggested to be caused by various reasons, such as the increased vegetation coverage, the reduced spring cyclone frequency, the rising surface air temperature, etc.
To figure out the relationship between East Asian dust emission and atmospheric dust over the Taiwan region, the annual mean and springtime (the dustiest season) Asian dust emissions and ND,d>500nm are given in Fig. 8. Fig. 8 shows that the annual and springtime Asian dust emissions have similar interannual characteristics as the Taiwan ND,d>500nm at lower levels, with r = ~0.8 and ~0.7 respectively, and dust emission also shows a similar decreasing trend as the Taiwan dust loading. Fig. 8(a) shows that the annual mean dust emission decreased by ~12% (decreasing rate of -0.01 kg m -2 yr -1 ) from 1999 to 2018. In Fig. 8(b), the spring dust emission shows a more significant decline, which decreased by ~17% (decreasing rate of -0.03 kg m -2 yr -1 ). The high correlation between that dust emission and Taiwan ND,d>500nm suggests that the decline of East Asian dust emission is one of the primary drivers of the ND,d>500nm decreasing trend over the Taiwan region. Besides the Asian dust, we also analyzed other dust sources, such as Sahara and the Middle East. However, the analysis shows no obvious correlations between dust in Taiwan and these regions.
As discussed above, that strong dust event frequency and atmospheric dust loading have distinct decreasing trends and interannual variations at low and high altitudes. At the same time, Fig. 8 shows that the East Asian dust emission has a high correlation with dust aerosols from the surface to 2 km, but with an insignificant correlation with dust at high altitudes. One possible reason is that the atmospheric dust over Taiwan has various transport modes at different heights.
Previous studies indicate that, over Taiwan, mineral dust outflowed from East Asia is mainly under 4 km (Lin et al., 2007;Liu et al., 2009;Chang et al., 2010;Lin et al., 2012), and at the midto-upper troposphere, dust loading is influenced by the transport from the Middle East and Sahara . At the same time, other possible reasons, such as the large-scale circulations, etc., may also influence the dust properties at different levels over Taiwan. Further studies will be processed to figure out possible influencing factors.

CONCLUSIONS
Based on the long-term (1999-2018) GEOS-Chem-APM model simulations, this study provides the seasonal variation, vertical properties, and long-term trend of the atmospheric dust aerosols over Taiwan. The 20-year model simulation is validated by comparing it with the in-situ dust measurement at the Cape Fuguei Research Station (25.30°N, 121.54°E) at the northern tip of Taiwan Island (Chou et al., 2017).
The result reveals that the atmospheric dust loading in Taiwan is highly influenced by dust events and that ND,d>500nm can vary by two orders of magnitude from clean to adjacent dusty days. During some extreme dust events, the ND,d>500nm can reach more than 20 cm -3 which could significantly influence the public health, air quality, radiation budget, and the heterogeneous freezing of cloud ice particles. It can be seen from the long-term analysis that, the mineral dust shows clear seasonal characteristics with a bimodal distribution: monthly mean ND,d>500nm and dust event frequencies are highest in spring (MAM) and fall-winter (ON). Monthly mean ND,d>500nm varies by two orders of magnitude from the lowest value in July (~0.03 cm -3 ) to the highest in April (~1.77 cm -3 ). Vertically, the dust loading has the peak value at low altitude (at ~1 km in spring-summer, and close to surface in fall-winter) and ND,d>500nm decreases significantly with height. It is worth noting that the monthly variations of mineral dust are different at low and high altitudes. The result shows that from the surface to the 8 km level, the frequency of strong dust events declines in winter and increases in summer. By analyzing the twenty-year GEOS-Chem simulations, we find a long-term decreasing trend of mineral dust over the Taiwan region, both for the frequency of strong dust events and the annual mean dust number concentration. The result indicates that, over Taiwan, the number of strong dust event days (NOE) decreases by ~40%, and the annual mean ND,d>500nm decrease by ~45% at low levels (under 4 km), and the decline tendency is less significant at upper levels (6-8 km), with NOE decreases by ~12% and the annual mean ND,d>500nm decrease by ~30%. At the same time, NOE and ND,d>500nm also show different interannual variation patterns at low (surface-4 km) and high altitudes (6-8 km). The twentyyear (1999-2018) analysis shows a significant correlation between Asian dust emission and the ND,d>500nm in this region. The long-range transported dust from Asia is the major source of the dust aerosols over Taiwan, and the decline of Asian dust emission is one of the primary drivers of the ND,d>500nm decreasing trend in this region. Over Taiwan, the distinct seasonal characteristics and long-term trends of mineral dust at low and high altitudes suggest that the dust transport modes and sources may shift vertically and seasonally. Previous studies indicate that, over Taiwan, the long-term transport of dust from East Asia arrives Taiwan mainly under 4 km (Lin et al., 2007;Liu et al., 2009;Chang et al., 2010;Lin et al., 2012), and at the mid-to-upper troposphere, the atmospheric dust loading is influenced by the dust from the Middle East and Sahara .
Numerous modeling and observational works indicate the impacts of mineral dust on climate and weather systems, and an increasing number of numerical models have taken account of dust aerosols. Our study reveals that, over Taiwan, the atmospheric dust loading is highly variable and can vary by about 1-3 orders of magnitude in different conditions. The significant space-time variations of mineral dust may result in significant impacts on cloud and precipitation, via direct and indirect effects, highlighting the necessity of more detailed and realistic dust number concentrations in the dust-aware weather forecasting models. We hope this study could provide some understandings of the properties and typical values of mineral dust over Taiwan and help to guide the assumptions of dust concentration in dust-related numerical studies.

ACKNOWLEDGMENTS
• This research was supported by the NSF Partnership for International Research and Education (PIRE) Program between the United States and Taiwan under grant OISE-1545917, and NASA under grants NNX17AG35G and 80NSSC19K1275. • The GEOS-Chem model is a community model managed by the Atmospheric Chemistry Modeling Group at Harvard University with support from NASA. The work described in this paper is based on GEOS-Chem version 10-01. GEOS-Chem is a freely accessible community model that can be downloaded from http://acmg.seas.harvard.edu/geos/. • The CALIPSO data used in this paper are obtained from the NASA Langley Research Center Atmospheric Science Data Center, and can be downloaded from https://eosweb.larc.nasa.go v/project/calipso/cal_lid_l2_vfm_standard_v4_20. • The long-term measurement of dust concentration at the Cape Fuguei Research Station is conducted by the Research Center of Environmental Changes (RCEC) with support from the Academia Sinica and the MOST, Taiwan.