Assessing the contribution of the ENSO and MJO to Australian dust activity based on satellite-and ground-based observations

Abstract. Despite Australian dust's critical role in the regional climate and
surrounding marine ecosystems, the controlling factors of the spatiotemporal
variations of Australian dust are not fully understood. Here we assess the
connections between observed spatiotemporal variations of Australian dust
with key modes of large-scale climate variability, namely the El
Niño–Southern Oscillation (ENSO) and Madden–Julian Oscillation
(MJO). Multiple dust observations from the Aerosol Robotic Network (AERONET),
weather stations, and satellite instruments, namely the Moderate Resolution
Imaging Spectroradiometer (MODIS) and Multi-angle Imaging SpectroRadiometer
(MISR), are examined. The assessed multiple dust observations consistently
identify the natural and agricultural dust hotspots in Australia, including
the Lake Eyre basin, Lake Torrens basin, Lake Frome basin, Simpson Desert,
Barwon–Darling basin, Riverina, Barkly Tableland, and the lee side of the Great
Dividing Range, as well as a country-wide, austral spring-to-summer peak in dust
activity. Our regression analysis of observed dust optical depth (DOD) upon an
ocean Niño index confirms previous model-based findings on the enhanced
dust activity in southern and eastern Australia during the subsequent austral
spring and summer dust season following the strengthening of austral
wintertime El Niño. Our analysis further indicates the modulation of the
ENSO–dust relationship with the MJO phases. During sequential MJO phases, the
dust-active center moves from west to east, associated with the eastward
propagation of MJO, with the maximum enhancement in dust activity at about 120,
130, and 140∘ E, corresponding to MJO phases 1–2, 3–4, and 5–6,
respectively. MJO phases 3–6 are favorable for enhanced ENSO modulation of
dust activity, especially the occurrence of extreme dust events, in
southeastern Australia, currently hypothesized to be attributed to the
interaction between MJO-induced anomalies in convection and wind and
ENSO-induced anomalies in soil moisture and vegetation.



Introduction
Australia represents a major contributor of dust to the Southern Hemisphere (Tanaka and Chiba, 2006), influencing regional climate and the marine ecosystems of the surrounding ocean basins.The substantial aerosol loading in the atmosphere from dust storms in Australia exerts a direct effect on the radiation budget through the absorption and scattering of incoming shortwave radiation and the absorption and emission of outgoing longwave radiation (Choobari et al., 2013;Miller et al., 2004;Tegen and Lacis, 1996).Dust aerosols also produce an indirect effect by influencing the nucleation, microphysics, development, and optical properties of clouds, thereby altering rainfall patterns (DeMott et al., 2010).The deposition of transported dust over ocean affects ocean biogeochemistry through changes to the iron supply (Gabric et al., 2010;Jickells et al., 2005).Australia's arid and semiarid regions provide a key supply of iron to the Southern Ocean and Antarctica by dust, during the present (Li et al., 2010), glacial (Lamy et al., 2014), and interglacial (Revel-Rolland et al., 2006) periods.Since most of the Southern Ocean is iron-limited (Sunda and Huntsman, 1997), the transport and deposition of Australian dust affect its productivity and carbon uptake (Boyd et al., 2004;Gabric et al., 2002).Therefore, deeper understanding of the spatiotemporal variations in Australian dust emissions and Published by Copernicus Publications on behalf of the European Geosciences Union.
Y. Yu and P. Ginoux: Assessing the contribution of ENSO and MJO to Australian dust activity their driving mechanisms will have broad implications on the regional and global climate.
The current understanding of Australian dust activity largely extends from interpretation of ground-based observations and satellite aerosol products, which have not been thoroughly intercompared over this region.Using a dust storm index (DSI) derived from dust storm and visibility observations made at Australian Bureau of Meteorology (BoM) stations (McTainsh et al., 1998), O'Loingsigh et al. (2014) assessed the spatial distribution of frequency and intensity of dust activity at 160 stations across Australia.Based on dust loading derived from the satellite aerosol optical depth (AOD), namely the Deep Blue algorithm applied to Moderate Resolution Imaging Spectroradiometer (MODIS) data from the polar-orbiting Terra and Aqua satellites, Ginoux et al. (2012) identified the Lake Eyre basin as the leading natural dust source in Australia, consistent with previous ground-based (McTainsh, 1989) and satellite-based (Bullard et al., 2008;Prospero et al., 2002) identification of dust sources.Ginoux et al. (2012) further identified agricultural dust sources in the Murray-Darling Basin in southeastern Australia, including the Victorian Big Desert, Riverina, and the Barwon-Darling basin, consistent with an earlier satellite-based dust source identification (Prospero et al., 2002) and model-based wind erodibility during dry years (Webb et al., 2006); however, these agricultural dust sources generated minimal dust storm frequency at nearby weather stations (McTainsh, 1989;Mc-Tainsh et al., 1989, 1998, 2007;O'Loingsigh et al., 2014).Ginoux et al. (2012) attributed the potential inconsistency in dust source maps among datasets to the various temporal coverage of each dataset.The apparent inconsistency between satellite-and ground-based spatial distribution of dust frequency and intensity could also be a result of the differential spatial coverage of these datasets.Therefore, it is critical to rigorously cross-validate these observations of dustiness in Australia.
Observations and general circulation models (GCMs) have shown substantial variability in the occurrence and intensity of dust emissions across Australia on interannual to decadal timescales, primarily driven by persistent anomaly in rainfall associated with Pacific sea surface temperature (SST) fluctuations, particularly El Niño-Southern Oscillation (ENSO) events (Bullard and Mctainsh, 2003;Evans et al., 2016;Lamb et al., 2009;Risbey et al., 2009;Strong et al., 2011;Webb et al., 2006).Corresponding to the ENSOinduced rainfall anomalies, during El Niño conditions, there is increased wind erosion in central and southeastern Australia, while during La Niña years the sources are shifted to the southwestern regions (Webb et al., 2006).Based on the Geophysical Fluid Dynamics Laboratory (GFDL) climate model (CM3), Evans et al. (2016) further uncovered the role of climate-vegetation interactions in amplifying and persisting ENSO's modulation on dust emission in southeastern Australia.The modulation on dust emission and loads was simulated to further amplify the ENSO-related rainfall vari-ability across eastern Australia (Rotstayn et al., 2011).On longer timescales, Lamb et al. (2009) revealed a pronounced and consistent dust maximum during 1959 to 1973 and a much more dust-free period after 1977 across central eastern Australia, based on the frequency of dust events reported at weather stations.This regime shift in Australian dustiness was attributed to wind field changes associated with oscillations in the Pacific climate system, including the latitudinal displacement of the South Pacific Convergence Zone, and SST changes of the Pacific Decadal Oscillation and North Pacific Oscillation (Lamb et al., 2009).From the paleoclimate perspective, the geochemical characteristics of Australian dust deposition in New Zealand over the last 8000 years have been used to identify corresponding dust sources, thereby inferring patterns of aridity in eastern Australia and climate variability associated with ENSO during the Holocene (Marx et al., 2009).
On intraseasonal to interseasonal timescales, the variability in Australia's atmospheric and terrestrial states features the dominant signal from the Madden-Julian Oscillation (MJO) (Notaro, 2018;Risbey et al., 2009;Wheeler et al., 2009;Yu and Notaro, 2020); yet, the potential role of MJO on Australia's dust emission and transport remains understudied.The MJO is characterized by eastward-propagating, large-scale wave-like disturbances in the tropics, particularly across the tropical Indo-Pacific region, with a typical cycle of 30-60 d (Madden andJulian, 1971, 1972).During an MJO event, anomalous convection acts as a tropical heat source, emitting stationary Rossby waves that propagate into the extratropics and significantly modulate the extratropical circulation (Matthews et al., 2004;Seo and Son, 2012).Previous modeling and observational analyses have identified significant influence of MJO on rainfall and circulation across Australia (Marshall et al., 2013;Risbey et al., 2009;Wheeler et al., 2009).An observational study by Wheeler et al. (2009) identified rainfall responses across extratropical Australia to MJO-induced changes in vertical motion and meridional moisture transport during austral summer and other seasons, respectively.Marshall et al. (2013) uncovered that the observed probability of an upper decile heat event varies according to MJO phase and time of year, with the greatest impact of the MJO on extreme heat occurring over southern Australia in austral spring during phases 2-3.The convectively active MJO phases 5-6 are also responsible for anomalous vegetation growth in the northern Australian monsoon region (Notaro, 2018), which further supports circulation changes over a broader region across the continent (Yu and Notaro, 2020).
Furthermore, MJO may interact with ENSO on the modulation of Australian dust emission.First, surface westerly anomalies introduced by the MJO can force downwelling oceanic Kelvin waves and influence ENSO evolution (Kessler et al., 1995;McPhaden and Taft, 1988), thereby triggering ENSO's modulation on Australian dust emission.Second, ENSO modulates MJO's propagation, resulting in dif-ferentiated spatiotemporal evolution of MJO (Wei and Ren, 2019) and its influence on Australia's regional climate.Furthermore, the convection and circulation anomalies introduced by MJO can enhance or weaken the basic response of regional climate to ENSO events, depending on the phase of MJO, as demonstrated by an observational study on extreme precipitation over northern South America (Shimizu et al., 2017).Despite MJO's critical influence on the regional climate, its direct or indirect role in modulating dust emission or concentration in Australia has, to our knowledge, never been explicitly investigated in either observations or models.
Motivated by the current knowledge gap in the large-scale climate drivers of Australian dust activity, the present study first assesses the multiple satellite and ground-based observations of dustiness in Australia, and then establishes the connection between the observed spatiotemporal variations in Australian dust activity with ENSO and MJO.We further provide hypotheses regarding ENSO and MJO's modulation on Australian dust activity, through ENSO's cumulative influence on vegetation and soil properties and MJO's shortterm perturbation on convection and extreme precipitation.Sections 2, 3, and 4 present the data and methods, results, and conclusions and discussion, respectively.
2 Data and methods

Dust optical depth proxies
Dust optical depth (DOD) is a column integration of extinction coefficient by mineral particles.The current study examines DOD proxies during 2000-2019 from MODIS onboard the polar-orbiting Terra and Aqua satellites, the Multiangle Imaging SpectroRadiometer (MISR) instrument (Diner et al., 1998) on Terra, and the ground-based AErosol RObotic NETwork (AERONET) (Holben et al., 1998) sun photometers, with distinct retrieval algorithms.

MODIS
Following Pu et al. (2020), daily DOD is retrieved from collection 6.1, level 2, MODIS Deep Blue aerosol products (Hsu et al., 2013;Sayer et al., 2013), including aerosol optical depth (AOD), single-scattering albedo (ω), and the Ångström exponent (α).All the daily variables are first interpolated to a 0.1 • × 0.1 • grid using the algorithm described by Ginoux et al. (2010).To account for dust's absorption of solar radiation and separate dust from scattering aerosols, such as sea salt, we require the single-scattering albedo at 470 nm to be less than 0.99 for the retrieval of DOD.Based on the size distribution of dust towards the coarse range and to separate it from fine particles, DOD is retrieved as a continuous function of AOD and the Ångström exponent: DOD = AOD × (0.98 − 0.5089α + 0.051α 2 ). ( This retrieval of DOD is on the basis of the Ångström exponent's sensitivity to particle size, with smaller values of Ångström exponent indicating larger particles (Eck et al., 1999), and the previously established relationship between Ångström exponent and fine-mode AOD (Anderson et al., 2005).In short, MODIS DOD represents the optical depth of absorbing, coarse-mode aerosols that are often dust over bare ground or sparsely vegetated regions.Details about the retrieval process and estimated errors are summarized by Pu and Ginoux (2018a).MODIS DOD products have been widely used for the identification and characterization of dust sources (Baddock et al., 2009(Baddock et al., , 2016;;Ginoux et al., 2012), as well as examination of variations in regional and global dustiness (Pu et al., 2019(Pu et al., , 2020;;Pu andGinoux, 2017, 2018a).Following the recommendation from Baddock et al. (2016) and previous applications of MODIS DOD (Pu et al., 2019(Pu et al., , 2020;;Pu andGinoux, 2017, 2018a), here we use DOD with a low-quality flag of QA = 1, under the assumptions that (1) dust sources are better detected using DOD with a low-quality flag and that (2) retrieved aerosol products are poorly flagged over dust source regions.

MISR
Benefiting from its multiangle observations, MISR data can be used to directly retrieve AOD and particle properties (Diner et al., 1998).In the current study, version 23, level 2, daily MISR 550 nm coarse-mode AOD (cmAOD) and nonspherical AOD (nsAOD) at 4.4 km resolution (Garay et al., 2020) are compared with MODIS DOD.The MISR nonspherical AOD fraction is often referred to as "fraction of total AOD due to dust", as dust is the primary nonspherical aerosol particle in the atmosphere, especially over desert regions such as those found in the arid and semiarid regions in Australia (Kalashnikova et al., 2005).The MISR nonspherical AOD has been used to examine variations in dustiness in North Africa and the Middle East (Yu et al., 2013(Yu et al., , 2016(Yu et al., , 2018a(Yu et al., , 2020)).Similar to our use of MODIS DOD with a low-quality flag, here we analyze the raw MISR cmAOD and nsAOD retrieval without quality filtering.MISR cmAOD and nsAOD are also interpolated to a 0.1 • × 0.1 • grid using the algorithm described by Ginoux et al. (2010).Due to its relatively narrow swath of ∼ 380 km, MISR samples the study region about every 10 d.

AERONET
The version 3, level 2 (cloud-screened and quality-assured), sub-daily AERONET coarse-mode AOD (cmAOD) at 500 nm obtained from 18 sun photometers across Australia (Giles et al., 2019) and retrieved by the Spectral Deconvolution Algorithm (SDA) (O'Neill et al., 2003) is analyzed here along with DOD from MODIS and cmAOD from MISR.In the analysis of annual mean and seasonal cycle, AERONET cmAOD monthly data are first screened by removing those https://doi.org/10.Annual mean and seasonal cycle are only analyzed for the 15 AERONET stations with at least 5 months of data for at least 3 years.Colocated cmAOD from AERONET and satellite products are further compared.Here a "collocated observation" is identified when there is available MODIS DOD or MISR cmAOD over the 0.1 • grid covering the AERONET site within ±0.5 h of the corresponding AERONET site observation.Although further spatial smoothing may improve the consistency between AERONET and satellite measurements (Yu et al., 2013), here we keep the fine satellite pixels to evaluate the accuracy of satellite products at their original spatial resolution.At each AERONET site, one satellite observation is often associated with multiple AERONET observations in time.In this case, AERONET observations are temporally averaged, resulting in only one pair of colocated and averaged satellite-AERONET DOD observations for a given colocated incident at each AERONET site.Larger temporal averaging windows, such as ±1 h, do not improve the consistency between satellite and AERONET measurements, likely due to the fine spatial scale considered in the current study.

DSI from weather stations
The present study analyzes meteorological records of dust activity, based on nine weather codes that are related to dust events as defined by the World Meteorological Organization (WMO).
where severe dust storm (SDS) = 1 if a decreased (code 33), stable (code 34), or begun/increasing (code 35) severe dust storm with visibility < 200 m is reported at least once and 0 otherwise; moderate dust storm (MDS) = 1 if a decreased (code 30), stable (code 31), or begun or increasing (code 32) slight or moderate sand or dust storm with visibility < 1000 m is reported at least once and 0 otherwise; and local dust event (LDE) = 1 if raised dust or sand (code 07), well-developed dust whirls (i.e., "dust devils", 08), or distant or past dust storm (code 09) is reported at least once and 0 otherwise.The credibility and temporal stability of DSI was evaluated in detail by O'Loingsigh et al. (2014).Similar to the AERONET data availability screening, annual mean and seasonal cycle are only analyzed for 182 weather stations with at least 5 months of effective data, namely with at least 5 d of DSI available during these months for at least 3 years during 2000-2019.

Seasonal cycle of dustiness
To achieve statistically meaningful analysis of the dustiness annual cycle, the peak month of each dustiness measure, namely DOD from MODIS, cmAOD and nsAOD from MISR, cmAOD from AERONET, and DSI from weather stations, is obtained via a two-step approach.First, a sinusoid function of month is fitted for each dustiness measure, where i stands for the calendar month (1 for January, 2 for February, . . ., and 12 for December); D(i) is the 20-year average dustiness in month i; and α, β, and γ are estimated by minimizing the square error between the predicted and observed D(i)'s (i = 1 to 12).Following this, the peak month of dustiness is obtained from the predicted dustiness among 12 months.The peak month is regarded as statistically meaningful only if (1) the predicted and observed seasonal cycle of dustiness are significantly correlated with correlation exceeding 0.58 (n = 12), based on the Student's t test at a significance level of 0.05, (2) the root-mean-square error between the predicted and observed dustiness is below a quarter of the annual mean dustiness, and (3) the amplitude of the predicted dustiness seasonal cycle (maximum minus minimum) exceeds half of the maximum value among 12 months.

Large-scale climate indices, environmental variables, and associated analysis
The connection between dust emission and large-scale climate indices is established here using MODIS DOD, which has both morning and afternoon coverage for almost all pixels in Australia, and station DSI during 2003-2019.Regression and composite analyses are conducted to assess the role of ENSO and MJO in regulating Australian dust emission by modulating soil and vegetation properties and perturbing convection and extreme precipitation activities, respectively.

Ocean Niño Index and regression analysis
To assess ENSO's modulation on Australian dustiness, an Ocean Niño Index (ONI) is analyzed.ONI is calculated as the 3-month running mean of Extended Reconstructed Sea Surface Temperature, version 5 (ERSSTv5) (Huang et al., 2017), SST anomalies in the Niño 3.4 region (5 • N-5 • S, 120-170 • W), based on centered 30-year base periods updated every 5 years (Climate Prediction Center, 2020).The influence of ENSO on DOD and DSI is quantified based on regression of seasonal average of daily DOD and occurrence of extremely high daily DOD during December-February (DJF) and September-November (SON) upon antecedent 3-month-averaged Niño 3.4 (sample size = 17 based on 17 years of data).The regression analysis is performed with stations that have more than 2 weeks of daily DSI during the focal season (DJF or SON) of at least 12 out of the 17 years.Although regression analysis does not directly infer causality, the resultant identification of covariability between Australian dust and antecedent ENSO state indicates higher likelihood of the latter driving the former than the opposite.To account for the non-Gaussian distribution of DOD and DSI, here the significance of regression coefficient is obtained through a Monte Carlo permutation test with 1000 iterations, following Yu and Notaro (2020).In each iteration, the time series of DOD or DSI is randomly scrambled, leading to a random estimate of the regression coefficient on ONI.The probability distribution function (PDF) of the random regression coefficients is used to test if the regressions in the original, non-permutated data are statistically significant.In the current study, a significance level of 0.05 is used to indicate statistically significant results.Given the insignificant autocorrelation at a 1-year lag with all the dust and environmental variables across the major dusty regions in central and southeastern Australia (Fig. S1 in the Supplement), the current statistical significance test does not account for the potential problem with random scrambling caused by autocorrelation.

Real-time multivariate MJO index and composite and regression analysis
To assess the potential influence of MJO and its interaction with ENSO on Australian dust activity, the real-time multivariate MJO index (RMM) (Wheeler and Hendon, 2004) is examined.RMM is derived as the principal components (PCs) corresponding to the leading two empirical orthogonal functions (EOF) of the combined fields of near-equatorially averaged 850 hPa zonal wind, 200 hPa zonal wind, and satellite-observed outgoing longwave radiation (OLR) data.Longer-timescale variability resulting from ENSO and other interannual variations with periods longer than about 200 d is removed prior to the EOF analysis.The combination of PC1 (RMM1) and PC2 (RMM2) of RMM reflects the magnitude and phase of the MJO.When the amplitude is greater than 1, eight MJO phases are determined by the sign of RMM1 and RMM2.Phases 1 and 2 mark the time when the MJO's convective envelope is centered near the western Indian Ocean, and phases 5-6 mark the time when the envelope is near northern Australia (Wheeler et al., 2009).Composite analysis is conducted for DOD, frequency of extremely high DOD, and DSI in each of the consecutive two MJO phases (phases 1-2, 3-4, 5-6, and 7-8) during the dust season of September-February.The results are expressed as the differences between the phase-specific DOD or DSI and the all-phase seasonal averages.The composite analysis is applied to stations that have more than 7 d of daily DSI in each MJO phase group during the dust season (September to February) of at least 12 out of the 17 years.A Monte Carlo bootstrap test with 1000 iterations is used to determine the significance of anomalies in dustiness during various MJO phases.In each iteration, daily dustiness measures are randomly sampled with the same size as a particular group of MJO phases.These randomly sampled dustiness measures are used to construct a PDF of sample mean dustiness and test if the mean dustiness during specific MJO phases is lower than the 2.5th or higher than the 97.5th percentile of the PDF.
Further, regression of dustiness upon ONI is performed for each MJO phase group to evaluate potential role of https://doi.org/10.5194/acp-21-8511-2021Atmos.Chem.Phys., 21, 8511-8530, 2021   In each iteration, daily dustiness measures are randomly sampled from the entire dust season with the same size as a particular group of MJO phases and averaged to obtain a random-phase mean dustiness measure for each year.The time series of these random-phase mean dustiness measures is regressed on the antecedent ONI, resulting in a PDF of the regression coefficients to test if the regression coefficient from the realistic, phase-specific dustiness is lower than the 2.5th or higher than the 97.5th percentile of the PDF.

Other environmental variables
To examine the potential mechanisms underlying the modulation of ENSO and MJO on Australian dustiness, we assess the connection between these large-scale climate drivers and various environmental factors such as surface wind speed, precipitation, soil moisture, and leaf area index (LAI) across Australia.The data sources of these environmental variables are outlined in Table 1.Regression and composite analyses are applied to these environmental variables, similar to those applied to the dustiness observations.To account for the non-Gaussian distribution of these environmental variables, the statistical significance of the regression and composite sig-nals are evaluated by the aforementioned non-parametric approaches.

Comparison of multiple observations of dustiness
MODIS DOD from both Terra and Aqua and station-based DSI consistently identify the natural and agricultural dust hotspots in Australia, including the Lake Eyre basin, Lake Torrens basin, Lake Frome basin, Simpson Desert, Barwon-Darling basin, Riverina, Barkly Tableland, and the lee side of the Great Dividing Range (Figs. 1 and 2).The annual mean MODIS DOD reaches 0.2 over Lake Eyre, Lake Torrens, and Lake Frome, where over 30 % of days observe a DOD exceeding 0.2, the 98th percentile of all MODIS DODs across Australia.MISR cmAOD and nsAOD are generally lower than MODIS DOD in magnitude and exhibit minimal spatial variation.Moreover, MISR only captures the margin of Lake Eyre basin and Barkly Tableland and shows relatively low dustiness over the Barwon-Darling basin and mostly invalid retrievals over the Lake Torrens and Lake Frome basins.The failure of MISR in capturing the agricultural dust hotspots in Australia is likely attributed to a lack of dust-smoke mixture in the look-up table involved in the operational MISR retrieving algorithm (Kahn et al., 2010).This shortage leads to less reliable retrievals of both cmAOD and nsAOD from MISR over the fire-prone dust hotspots in southeastern Australia (Garay et al., 2020).The spatial distribution of mean cmAOD from AERONET is largely consistent with the satellite observations.The apparently high DOD from MODIS and cmAOD from AERONET over the coastal region is https://doi.org/10.5194/acp-21-8511-2021Atmos.Chem.Phys., 21, 8511-8530, 2021 likely caused by the abundance of sea salt aerosol and its mixture with dust and biomass burning aerosols.DOD from MODIS, cmAOD from AERONET and MISR, and nsAOD MISR exhibit a generally consistent seasonal cycle in dust activity, which peaks in austral spring to summer, namely November, December, and January, across most of the country (Figs. 3 and 4).In particular, the seasonal cycle in DOD proxies is generally consistent between all satellite instruments and AERONET at most sites in Australia.The seasonal cycle in DOD proxies and DSI are highly consistent in Birdsville and Tinga Tingana, located near the dust hotspots in Simpson Desert and Lake Eyre basin, respectively (Fig. 4).The largest disagreement between satellite and station-based observations of dustiness occurs over the Barwon-Darling basin and its northern downwind regions in eastern Australia, where MODIS DOD and station DSI indicate peak dustiness in November to December, while cmAOD and nsAOD from MISR indicates undetectable seasonal cycle peak.The inconsistency between MISR and other sources of dustiness measure is likely attributed to the uncertainty in MISR's retrieval of dust-smoke mixtures (Kahn et al., 2010).Over the Lake Eyre-Torrens-Frome basin, the morning satellite observations, namely MODIS-Terra and MISR, display a summertime peak in dustiness, while the afternoon satellite observation from MODIS-Aqua and all-day station observations indicate a springtime peak.This contrast between the seasonal cycles in morning and afternoon dustiness suggests a seasonally varying diurnal evolution of dust emission in southern central Australia.
The general comparison between colocated satellite DOD or cmAOD and AERONET cmAOD exhibits reasonable quality of satellite retrievals over the majority of Australia, but wider spreads of DOD from both MODIS-Terra and MODIS-Aqua, and cmAOD from MISR, especially corresponding to colocated high cmAOD from AERONET (Figs. 5 and 6).The wide spread of MISR cmAOD, compared with colocated AERONET cmAOD, is partly attributed to the limited spectral range of MISR.Very few MODIS DOD retrievals reach lower than 0.005, likely due to the numeri-cal limits of retrieving algorithm.Furthermore, both MODIS and MISR display a wider spread at higher DOD or cmAOD and an overall underestimation, especially when AERONET DOD exceeds 0.1 (Fig. 5).This underestimation of high optical depth has been reported by previous global validations of total AOD from MODIS (Sayer et al., 2019;Wei et al., 2019) and MISR (Garay et al., 2020), as well as MODIS DOD (Pu and Ginoux, 2018b).The underestimation of high DOD potentially leads to the deteriorated correlations between colocated satellite DOD or cmAOD and AERONET DOD over the dustiest region near the Lake Eyre basin, compared with less dusty regions in Australia (Fig. 6).Given the distinct retrieval algorithms involved in the satellite DOD, cmAOD, and AERONET cmAOD, the moderate but significant correlations (p < 0.001) between colocated, thousands-satellite DOD or cmAOD and AERONET cmAOD (Fig. 5   According to the comparison between colocated MODIS DOD and MISR cmAOD and nsAOD, as well as satellite DOD proxies versus station dust storm observations, coarsemode AOD is an acceptable approximation of dustiness over the dust hotspots in central and southeastern Australia (Figs. 7 and 8).While the correlation between MODIS DOD and MISR cmAOD exceeds 0.4 across the majority of the country, even exceeding 0.6 over central eastern Australia, the correlation between MODIS DOD and MISR nsAOD reaches 0.6 only near the major dust source regions, e.g., the Lake Eyre-Torrens-Frome basin, the Simpson Desert, and northern downwind area of the Barwon-Darling basin.Fur-thermore, higher MODIS DOD is observed on reported dusty days at most stations in Australia, especially during countrywide local dust and moderate dust events, as well as severe dust events in the south and east with MODIS-Aqua (Fig. 8).Insufficient number of collocations between ground observations and MISR overpasses likely leads to the minimal signal in MISR cmAOD or nsAOD, as previously shown in the annual mean (Fig. 2).
Overall, the general consistency between MODIS DOD and colocated AERONET cmAOD, MISR cmAOD, and MISR nsAOD and qualitative consistency between MODIS DOD and station dust observations provides confidence in

Modulation from ENSO and MJO on Australian dustiness
According to regression analysis applied to multiple dust observation data sets and Niño 3.4 at various antecedent time, austral wintertime El Niño supports enhanced dust activity in southern and eastern Australia during the subsequent austral summer dust season (Fig. 9).An SST anomaly of +1 • C in the Niño 3.4 region during July to September (JAS) leads to an increase in daily mean DOD of about 0.05 over the lee side of the Great Dividing Range, 0.04 over the Barwon-Darling basin, and 0.03 over Riverina during the subsequent December to February (DJF).The +1 • C warming in the Niño 3.4 region during JAS causes an increase in the frequency of extreme DOD of about 5 % over the lee side of the Great Dividing Range and 2 % over Riverina, and an increase in DSI of about 1 % over the Barwon-Darling basin and Riverina during the subsequent DJF.The El Niño condition in the austral autumn and winter also leads to enhanced dust emissions across the Simpson Desert and Barwon-Darling basin during the subsequent September-November (SON).The lagged response in Australian dust activity to ENSO is attributed to ENSO's persistent and cumulative influence on the regional soil moisture and LAI (Fig. 10).The El Niñoinduced inhibited rainfall across eastern and central Australia (e.g., Risbey et al., 2009)   causes a reduction in surface roughness and strengthened surface wind that further enhances dust emission.The response magnitude in soil moisture and LAI in austral summer and spring peaks after 3-6 months of the ENSO signal, supporting the 3-6 month lag in the dustiness response to ENSO.The currently identified importance of vegetation in the modulation of ENSO on dust emissions in Australia confirms the model-based finding about the role of climate-vegetation interactions in amplifying and persisting ENSO's modulation on dust emission in southeastern Australia by Evans et al. (2016).
According to the composite of DOD, frequency of extreme DOD, and station-based DSI, the dust-active center moves from west to east associated with the eastward propagation of MJO, with a maximum enhancement in dust activity at about 120, 130, and 140 • E, corresponding to MJO phases 1-2, 3-4, and 5-6, respectively (Fig. 11).During MJO phases 5-6, i.e., the convection-active phases for Australia, the increased surface wind speed over the majority of the continent, especially over the dust hotspots in the Lake Eyre-Torrens-Frome basin and Riverina, appears responsible for the enhanced dustiness (Fig. 12).Surprisingly, the enhanced dustiness over the central and eastern Australian dust hotspots seems to be associated with anomalously wet conditions during all MJO phases.Given that central and southern Australia generally receive less than 1 mm of rainfall on an average day, we hypothesize that over these arid or semiarid regions, enhanced rainfall during the MJO phases 3-6 in austral spring and summer associated with enhanced convection and occurrence of thunderstorms supports higher occurrence of haboob type of dust events.Several case studies have reported haboob dust events in the central and eastern Australia  9a-h.(McTainsh et al., 2005;Shao et al., 2007).Strong et al. (2011) found that about 24 % of dust storms in the lower Lake Eyre basin during 2005-2006 are associated with thunderstorms.Our alternative hypothesis relies on the supply of fine particles by occasional flooding from MJO-induced storms.For supply-limited and/or transport-limited dust sources such as those in southeastern Australia, lack of occasional storms under drier conditions usually leads to the failure of sediment replenishment, thereby leading to anomalously inactive dust emission (Arcusa et al., 2020;Bullard and Mctainsh, 2003).
ENSO's regulation of dust emission varies in magnitude by MJO phases, with MJO phases 3-6 favorable for enhanced ENSO regulation on dust activity, especially the occurrence of extreme dust events, in southeastern Australia (Fig. 13).An SST anomaly of +1 • C in the Niño 3.4 region in austral winter is associated with an increased DOD by over 0.05, an increased frequency of extreme DOD by over 5 %, and an increased DSI by 2 % over the Barwon-Darling basin and Riverina during MJO phases 5-6 in austral spring and summer.MJO phases 3-4 feature a moder-ately enhanced dustiness over the Lake Eyre-Torrens-Frome basin in response to antecedent El Niño.We hypothesize that the enhanced response in dustiness across southeastern Australia to ENSO during MJO phases 3-6 can be attributed to the interplay between MJO-induced anomalies in convection, rainfall, and wind and the ENSO-induced anomalies in soil moisture and vegetation.While the dry soils and diminished vegetation caused by El Niño provide favorable conditions for dust emission (Fig. 10), the active convections and elevated occurrence of extreme precipitation during MJO phases 3-6, as well as strengthened surface wind during MJO phases 5-6, likely trigger more dust emission and extreme dust events across southeastern Australia through either haboob type of dust events or additional sediment supply by occasional flooding (Fig. 12).

Discussion and conclusions
The current study investigates the contribution of large-    the modulation of Australian dust activity on the intraseasonal to interannual time scales.Multiple sources of dustiness measurements, namely DOD from MODIS, cmAOD and nsAOD from MISR, cmAOD from AERONET, and DSI from weather stations, are intercompared in terms of their annual mean, seasonal cycle, and day-to-day variations over a 20-year period from 2000 to 2019.These assessed dust observations consistently identify the natural and agricultural dust hotspots in Australia, including the Lake Eyre-Torrens-Frome basin, Simpson Desert, Barwon-Darling basin, Riverina, Barkly Tableland, and the lee side the Great Dividing Range, and a country-wide dust peak during austral spring-to-summer, confirming the previous ground-based (McTainsh, 1989) and satellite-based (Bullard et al., 2008;Ginoux et al., 2012;Prospero et al., 2002) identification of dust sources.Furthermore, the intercomparison between the multiple dust observations demonstrates the credibility of MODIS DOD -a widely analyzed satellite dust observation with optimal temporal and spatial coverage -over the arid to semiarid regions in central and southeastern Australia.Regression analysis of MODIS DOD upon Niño 3.4 SST confirms the previous model-based findings by Evans et al. (2016) on the enhanced dust activity in southern and eastern Australia during the subsequent austral summer dust season following the strengthening of austral wintertime El Niño.Composites of dustiness during sequential MJO phases demonstrates the propagation of the dust-active center from west to east associated with the eastward movement of MJO, with maximum enhancement in dust activity at about 120, 130, and 140 • E corresponding to MJO phases 1-2, 3-4, and 5-6, respectively.Our analysis further indicates the modulation of the ENSO-dust relationship with the MJO phases; MJO phases 3-6 are favorable for amplifying ENSO's modulation on dust activity, especially the occurrence of extreme dust events in southeastern Australia.
Although the current study demonstrates the general reliability of MODIS DOD over the arid and semiarid regions in Australia, uncertainties of this product should be noted.For example, the retrieval of MODIS DOD relies on the light-absorbing and coarse-mode nature of dust and is unable to distinguish between dust and the coarse-mode part of biomass burning aerosols (e.g., Noyes et al., 2020) et al., 2017).Given the potential contamination from biomass burning aerosols, our interpretation of the currently examined connection between dust and climatic drivers mainly focuses on central and southeastern Australia.In addition, for haboob dust events, which often occur with the presence of convective clouds, MODIS and MISR algorithms are unlikely to perform aerosol retrievals.Given the single assumption on dust particle shape involved in nsAOD, the MISR nsAOD is often regarded as a better proxy of DOD than coarse-mode AOD.But the limited temporal coverage of MISR makes it less useful for studying the day-to-day variations and extreme events of dust activity, especially corresponding to MJO.Typically, MISR only samples about 5 d during each MJO phase group (phases 1-2, 3-4, 5-6, and 7-8) per dust season (September to February) over most pixels in Australia.Furthermore, the retrieval of the dust-smoke mixtures, typically present over the southeastern shrublands and grasslands in Australia, is subject to huge uncertainty in the operational MISR aerosol product (Garay et al., 2020;Kahn et al., 2010).Therefore, MISR cmAOD and nsAOD are analyzed here only to support the reliability of MODIS DOD in representing dust activity.Overall, the optimal spatial and temporal coverage of MODIS aerosol products with an over 20-year record warrant its application for studying the spatiotemporal variations and environmental drivers of global aerosol loads.The current analysis on the connection between environmental factors, such as LAI, soil moisture, wind, and precipitation, and ENSO and MJO leads to the hypothetical mechanisms underlying the identified modulation of ENSO and MJO on Australian dustiness.We hypothesize that the dry soils and diminished vegetation resulting from the El Niñoinduced rainfall reduction provide favorable conditions for dust emission during the subsequent season; the enhanced convective activity and strengthened surface wind during MJO phases 3-6 likely triggers more dust emission and extreme dust events across southeastern Australia during the El Niño-associated dry years, thereby amplifying ENSO's modulation on dust emission.Under the hypothesized mechanism, we expect more pronounced MJO enhancement of ENSO's modulation on dust following El Niño than La Niña conditions.One explanatory hypothesis for this relationship builds partly on the occurrence of haboob dust storms and their connection with MJO-induced anomalies in deep convection over southeastern Australia.An alternative hypothesis relies on the supply of sediments by MJO-induced storms and their resulting occasional flooding.Our results shed light on a potential linkage between extreme precipitation and enhanced dust emission in Australia.These alternatives motivate further evaluation of these hypothesized mechanisms underlying the modulation of ENSO and MJO on dust emissions across Australia in an Earth system model.Furthermore, the current hypotheses regarding the influence of ENSO and MJO are established upon regression and composite analyses, which do not directly infer causality.Advanced statistical approaches, such as the Stepwise Generalized Equilibrium Feedback Assessment (SGEFA) (Yu et al., 2017(Yu et al., , 2018b)), will be useful to evaluate the role of large-scale climate modes and local environmental changes in the emission and transport of Australian dust.In addition, the present study focuses on the natural drivers of Australian dust activity, while anthropogenic dust emission from land use change is a key contributor to total dust emission in Australia (Ginoux et al., 2012;Tegen et al., 2004;Webb and Pierre, 2018).Indeed, disturbed soil and vegetation from land use, such as pastoral and agricultural activity in eastern Australia, have caused substantial increase in dust emission and deposition during the 20th century (Brahney et al., 2019;Cattle, 2016).The modulation of land use on dust emission and transport from Australia may also be quantified and compared with natural drivers through future Earth system modeling.
Author contributions.YY conceived the study, analyzed the data, and wrote the manuscript with contributions from PG. PG retrieved the MODIS DOD data from MODIS Deep Blue aerosol products.

Figure 1 .
Figure 1.Dominant land cover type across Australia in 2012, marked with key dust source regions.The land cover data is from MODIS, based on the International Geosphere-Biosphere Programme (IGBP) classification scheme (Sulla-Menashe and Friedl, 2018).The land cover fraction of 18 classes is interpolated from 0.05 to 0.1 • .

Figure 2 .
Figure 2. Annual mean dust activity over Australia during 2000-2019.In (a-d), the map represents average DOD from (a) MODIS-Terra and (b) MODIS-Aqua and (c) coarse-mode AOD (cmAOD) and (d) non-sphere AOD (nsAOD) from MISR.Filled circles represent cmAOD from 15 AERONET sites, which are identical in (a-d).In (e-h), the map represents frequency of DOD proxies exceeding the 98th percentile of all observations from each instrument, namely 0.205 for DOD from MODIS-Terra, 0.198 for DOD from MODIS-Aqua, 0.095 for cmAOD from MISR, and 0.062 for nsAOD from MISR.Filled circles represent the Dust Storm Index (DSI, %) at 1489 weather stations, which is identical in (e-h).Dots in (e-h) represent stations with annual mean DSI smaller than 2 %.Grey indicates pixels with satellite sampling less than 10 d during the analyzed 20 years.

Figure 3 .
Figure 3. Month of peak dust activity over Australia during 2000-2019.In (a-d), the map represents average DOD from MODIS-Terra and MODIS-Aqua and cmAOD and nsAOD from MISR.Filled circles represent DOD from 15 AERONET sites, which are identical in (a-d).In (e-h), the map represents frequency of DOD proxies exceeding the 98th percentile of all observations from each instrument.Filled circles represent DSI from 182 weather stations.Only statistically meaning seasonal cycle peak month is shown for each dustiness measure, as outlined in Sect.2.3.

Figure 4 .
Figure 4. Seasonal cycle of dust activity at 15 AERONET sites in Australia.Lines represent mean DOD from MODIS-Terra (blue) and MODIS-Aqua (green), cmAOD from MISR (solid orange) and nsAOD from MISR (dashed orange), and cmAOD from AERONET (black) by month, referring to the left y axis.Bars, referring to the right y axis, represent mean DSI over weather stations located within 100 km of each AERONET site.The seasonal cycle of dust activity from satellite and in-site observations is obtained from the active years of each AERONET site.The AERONET sites are presented by longitude from west to east, matching the letters in Fig. 3a.

Figure 5 .
Figure 5.Comparison of MODIS DOD and MISR cmAOD against colocated AERONET cmAOD.(a-c) Joint probability density (%) of colocated cmAOD from AERONET and (a) DOD from MODIS-Terra, (b) DOD from MODIS-Aqua, and (c) cmAOD from MISR.(d-f) Boxplot of the difference in colocated measurements between (d) DOD from MODIS-Terra, (b) DOD from MODIS-Aqua, and (c) cmAOD from MISR and cmAOD from AERONET, as a function of AERONET cmAOD.The boxplots show the 5th, 25th, 50th, 75th, and 95th percentiles of the DOD difference.The sample size (N), correlation coefficients (R), and root-mean-square error (RMSE) with colocated AERONET cmAOD are provided in (d-f) for the corresponding satellite product.

Figure 6 .
Figure 6.Correlation between colocated hourly cmAOD from AERONET and satellite measurements at 18 AERONET sites in Australia.The satellite measurements include (a) DOD from MODIS-Terra, (b) DOD from MODIS-Aqua, and (c) cmAOD from MISR.The sizes of the dots indicate the number of colocated observations.A missing circle in (c) indicates no collocation between MISR and AERONET observations.

Figure 7 .
Figure 7. Correlation between colocated, daily MODIS-Terra DOD and MISR (a) cmAOD and (b) nsAOD.The number of collocations per pixel between MODIS and MISR during 2000-2019 generally varies from 100 to 800 for DOD and 20 to 200 for nsAOD.Grey indicates areas with fewer than 10 collocations per pixel.

Figure 8 .
Figure 8. Difference in satellite DOD proxies between dusty days and clear days at weather stations in Australia.The color of the filled dots represents the difference in (a, e, i) DOD from MODIS-Terra, (b, f, j) DOD from MODIS-Aqua, (c, g, k) cmAOD from MISR, and (d, h, i) nsAOD from MISR between days with no reported dust observation and days with reported (a-d) "local dust", (e-h) "moderate dust", and (i-l) "severe dust".The size of the dots indicates number of days with both weather observations and valid satellite retrievals within 25 km at each station.A plus sign indicates a significant positive difference in DOD between dusty and clear days based on a Monte Carlo bootstrap test (p < 0.05).

Figure 9 .
Figure 9. Regression of anomalies in seasonal dust activity in (a, b, e, f, i, j) December-February (DJF) and (c, d, g, h, k, l) September-November (SON) upon antecedent Niño 3.4.Analyzed dust variables include seasonal (a-d) DOD averaged from MODIS-Terra and MODIS-Aqua, (e-h) frequency of daily DOD anomaly exceeding 3 times the interannual standard deviation, and (i-l) DSI.(a, c, e, g, i, k) Regression coefficient between (a, e, i) DJF dust and antecedent July-September (JAS) Niño 3.4 (ENSO leading dust for 5 months), and (c, g, k) SON dust and antecedent May-July (MJJ) Niño 3.4 (ENSO leading dust for 4 months).In (a, c, e, g), the stitches indicate regions with significant regression coefficient (p < 0.05), based on Monte Carlo bootstrap test; and the slashes further denotes regions with significant positive correlation between MODIS DOD and MISR nsAOD.In (i, k), only statistically significant regression coefficients (p < 0.05) are shown.(b, d, f, h, j, l) Boxplot of the regression coefficient of seasonal dustiness on Niño 3.4 at different antecedent times as a function of the leading time of Niño 3.4, showing the 5th, 25th, 50th, 75th, and 95th percentiles of regression coefficient at all pixels within the dust source region (25-35 • S, 135-155 • E).

Figure 10 .
Figure 10.Regression of anomalies in seasonal LAI, soil moisture, and surface wind speed in (a, b, e, f, i, j) December-February (DJF) and (c, d, g, h, k, l) September-November (SON) upon Niño 3.4.(a, c, e, g, i, k) Regression coefficient between (a) LAI, (e) soil moisture, and (i) surface wind speed during DJF and Niño 3.4 during antecedent July-September (JAS) (ENSO leading environmental conditions for 5 months), and SON (c) LAI, (g) soil moisture, and (k) surface wind speed during SON and Niño 3.4 during antecedent May-July (MJJ) (ENSO leading environmental conditions for 4 months).Figure elements are the same as Fig. 9a-h.

Figure 11 .
Figure 11.Mean anomalies (% climatology) in dust activity by MJO phase in September-February.Analyzed dust variables include (a-d) DOD averaged from MODIS-Terra and MODIS-Aqua, (e-h) frequency of daily DOD anomaly exceeding 3 times the interannual standard deviation, and (i-l) DSI.The anomalies for each consecutive two MJO phases, namely phases 1-2, 3-4, 5-6, and 7-8, are calculated as the percentage differences between these two phases and the long-term average during September-February of 2003-2019.The composites consist of 376, 492, 613, and 450 d with RMM > 1 for phases 1-2, 3-4, 5-6, 7-8, respectively.In (a-h), the stitches indicate regions with significant percentage difference with the climatology (p < 0.05), based on a Monte Carlo bootstrap test, and the slashes further denote regions with significant positive correlation between MODIS DOD and MISR nsAOD.In (i-l), only statistically significant percentage differences (p < 0.05) are shown.

Figure 12 .
Figure 12.Mean anomalies in daily precipitation and surface wind speed by MJO phase in September-February.Figure elements are the same as in Fig. 11a-h.Grey dashes indicates areas with seasonal mean rainfall less than 1 mm d −1 .

Figure 13 .
Figure 13.Regression of anomalies in dust activity during each consecutive two MJO phases in September-February upon the antecedent June-August (JJA) Niño 3.4 by MJO phase.Analyzed dust variables include (a-d) DOD averaged from MODIS-Terra and MODIS-Aqua, (eh) frequency of daily DOD anomaly exceeding 3 times the interannual standard deviation, and (i-l) DSI.The anomalies for each consecutive two MJO phases, namely phases 1-2, 3-4, 5-6, and 7-8, are calculated as the differences between these two phases and the long-term average during September-February of 2003-2019.In (a-h), the stitches indicate regions with significant regression coefficient (p < 0.05), based on the Monte Carlo bootstrap test for the null hypothesis of a phase-invariant regression coefficient, and the slashes further denote regions with significant positive correlation between MODIS DOD and MISR nsAOD.In (i-l), only statistically significant regression coefficients (p < 0.05) are shown.

Table 1 .
Summary of the observational datasets for environmental variables applied in the current analyses across Australia.
) demonstrate the reliability of MODIS DOD and MISR cmAOD in representing coarse-mode aerosol loads.