The impact of Great Basin Desert dust on the summer monsoon system over southwestern North America

Introduction Conclusions References


Introduction
The North American Monsoon (NAM) system is a major climate system that significantly affects the regional hydrological cycle over large areas of the southwestern United States (US) and northwestern Mexico.maritime sources that together produce convection that contributes over 70 % of annual precipitation to the region.The NAM summer rainfall typically begins in July with a pronounced increase from an extremely dry condition and lasts until September when the dry regime reestablishes (Adams and Comrie, 1997).The Great Basin Desert (GBD), the largest US desert, covers an arid expanse of ∼190 000 square miles and is bordered by the Sierra Nevada Range on the west and the Rocky Mountains on the east, the Columbia Plateau to the north and the Mojave and Sonoran deserts to the south.The Great Basin Desert is north of the NAM region.
The Great Basin Desert is the main source of dust aerosol over the southwestern US.Dust aerosol has important climate impact through changes in solar and terrestrial radiation and radiative and microphysical properties of clouds (e.g.Sokolik et al., 1998;Ginoux et al., 2001;Lau et al., 2009;Zhao et al., 2011).The dust lofted from the deserts near the monsoon regions may affect monsoon circulation and precipitation worldwide.Recent years, many efforts have begun to study dust impacts on the West African Monsoon (WAM) and the Asian Monsoon (AM) circulation and precipitation, using general circulation models (GCMs) and regional models combined with observations (e.g.Miller et al., 2004;Yoshioka et al., 2007;Lau et al., 2006Lau et al., , 2009;;Kim et al., 2010;Zhao et al., 2011).These studies have provided important understanding of the fundamental processes associated with dust effect on the two monsoon systems.Over West Africa, Lau et al. (2009) and Kim et al. (2010) showed that the Saharan dust radiative effect could strengthen the WAM, which is manifested in a northward shift of the West African precipitation over land.Kim et al. (2010) and Zhao et al. (2011) also found that dust could affect not only the total amount but also diurnal cycle of WAM precipitation simulated by a GCM and a regional model, respectively.Over Asia, Lau et al. (2006) elucidated a plausible "heat-pump" mechanism for dust impact on the Asian summer monsoon using the NASA finite-volume GCM.As dust from the deserts of western China, Pakistan, and the Middle East are transported into and stacked up against the northern and southern slopes of the Tibet Plateau, the dust heated air over the slopes can draw in warm and moist air over the Indian subcontinent, suggesting Figures that the dust loading may lead to an advance of the rainy periods and subsequently an intensification of the Indian summer monsoon.
Although the dust emitted from the Great Basin Desert over the Southwest US is less than that emitted from the Saharan desert and the deserts over Asia, the NAM system is also weaker than the monsoon systems over West Africa and Asia.Therefore, dust could still have significant effect on the NAM system.The IPCC AR4 climate model simulations projected warmer and dryer conditions for the southwestern US and northwestern Mexico in the 21st century.The projected drought parallels the severity of the 1930s Dust Bowl (Seager et al., 2007).The mega-droughts may lead to widespread dust storms (e.g.Woodhouse and Overpeck, 1998) that interact with the NAM system, strengthening or weakening the monsoon system through positive or negative feedbacks.To the best of our knowledge, the potential dust impact on NAM is not fully understood.
To improve understanding of dust impact on NAM in current and future climate conditions, as the first step, we used the regional model WRF-Chem to investigate the dust impact on the NAM system in the current climate condition.The objective of this study is to investigate the radiative forcing of GBD dust and how it affects the NAM circulation and precipitation at a regional scale.This study (1) examines both SW and LW radiative forcings of GBD dust during the NAM season over the Southwest US, and ( 2 WRF-Chem, a version of WRF (Skamarock et al., 2008), simulates trace gases and particulates simultaneously with the meteorological fields (Grell et al., 2005).Grell et al. (2005) implemented the RADM2 (Regional Acid Deposition Model 2) photochemical mechanism (Stockwell et al., 1990) and the MADE/SORGAM (Modal Aerosol Dynamics Model for Europe (MADE) and Secondary Organic Aerosol Model (SORGAM)) aerosol model (Ackermann et al., 1998;Schell et al., 2001) into WRF-Chem.MADE/SORGAM in WRF-Chem uses the modal approach with three modes (Aitken, accumulation, and coarse modes, assuming a log-normal distribution within each mode) to represent the aerosol size distribution.Different aerosol species within each size mode are assumed to be internally mixed so that all particles within a size mode have the same chemical composition.In MADE/SORGAM, aerosol species include sulfate, nitrate, ammonium, black carbon (BC), organic matters (OM), sea salt, mineral dust, and water.Aerosol optical properties for scattering (e.g.single-scattering albedo, asymmetry factor, and extinction) are computed as a function of wavelength and three-dimensional position.A complex refraction index is calculated by volume averaging for each mode for each chemical constituent of the aerosol.To efficiently compute the extinction efficiency (Q e ) and the scattering efficiency (Q s ), WRF-Chem follows a methodology described by Ghan et al. (2001) that performs the full Mie calculations once first to obtain a table of seven sets of Chebyshev expansion coefficients, and later the full Mie calculations are skipped and Q e and Q s are calculated using bilinear interpolation over the Chebyshev coefficients stored in the table.A detailed description of the computation of aerosol optical properties in WRF-Chem can be found in Fast et al. (2006) and Barnard et al. (2010).
The version of WRF-Chem model with MADE/SORGAM coupled with GOCART dust scheme was first used by Zhao et al. (2010) to investigate the sensitivities in simulating the size distributions and optical properties of Saharan dust to emitted dust size distributions and aerosol size treatments during the dry season (from January to February)  Zhao et al. (2010Zhao et al. ( , 2011)).The refractive index of dust is set to 1.53 + 0.003i as Zhao et al. (2010Zhao et al. ( , 2011) ) over West Africa because of insufficient information to constrain the values over the Southwest US.The rapid radiative transfer model (RRTMG) (Mlawer et al., 1997;Iacono et al., 2000) for both SW and LW radiation is used to simulate aerosol direct radiative effect.The Morrison scheme and Grell scheme are used to represent cloud microphysics and convection processes.
Although the aerosol activation to cloud droplet is included to account for cloud chemistry and aerosol wet deposition (Gustafson et al., 2007), dust indirect radiative effect on the NAM system is not investigated in this study, because the simulated dust coincides minimally with the NAM precipitation region and dust effect by acting as ice nuclei (IN) is not treated in the current version of WRF-Chem.The Noah land surface model and Mellor-Yamada-Janjic (MYJ) planetary boundary scheme are used.Large-scale meteorological fields are assimilated with lateral boundary and initial conditions from the North American Regional Reanalysis (NARR) data, which also provide the prescribed sea surface temperature (SST) for the simulations.Chemical lateral boundary conditions are from the default profiles in WRF-Chem, which are based on the averages of mid-latitude aircraft profiles from several field studies over the eastern Pacific Ocean (McKeen et al., 2002).Anthropogenic emissions are obtained Introduction

Conclusions References
Tables Figures

Back Close
Full

AERONET surface observation network
The Aerosol Robotic Network (AERONET) (Holben et al., 1998) has ∼100 identical globally distributed sun-and sky-scanning ground-based automated radiometers, which provide measurements of aerosol optical properties throughout the world (Dubovik and King, 2000;Dubovik et al., 2002).In this study, AERONET measured aerosol optical depth (AOD) at 675 nm and 440 nm from one Arizona site (Maricopa, 111.97 • W, 33.07 • N, Fig. 1) near the dust source region are used to derive the AOD at 550 nm (using the Angstr öm exponent) for comparison with model results.All of the retrievals of AOD are quality level 2, and the uncertainty of AOD measurements is about ±0.01 (Holben et al., 2001).In this study, the available data from 2003 to 2009 are used to evaluate the modeling results during the same period.

MODIS
The Moderate Resolution Imaging Spectroradiometer (MODIS) instruments with wide spectral range, high spatial resolution, and near daily global coverage are designed on board the NASA Aqua platforms to observe and monitor the Earth changes including tropospheric aerosols (Kaufman et al., 1997).The standard MODIS aerosol product does not retrieve aerosol information over bright surfaces (e.g.desert) due to a strong surface spectral contribution in the visible range (Kaufman et al., 1997).However, a new algorithm, called the "Deep Blue algorithm" (Hsu et al., 2006), has been integrated Figures

Back Close
Full with the existing MODIS algorithm to retrieve AOD even over bright surfaces.Therefore, the retrieved "deep blue" AOD from MODIS (Collection 5.1) (only available over land so far) (Levy et al., 2007;Remer et al., 2005) is used in this study.The MODIS on board the Aqua platform passes over the equator at ∼13:30 LT during daytime (Kaufman et al., 1997).In this study, the available data from 2003 to 2009 are used to evaluate the modeling results during the same period.Model results are sampled in the same overpass time as Aqua when comparing with MODIS retrievals.

MISR
Since February 2000, the Multi-angle Imaging SpectroRadiometer (MISR) instrument on board the NASA Terra platform observes AOD continuously at nine distinct zenith angles, ranging from 70 • afterward to 70 • forward, and in four narrow spectral bands centered at 446, 558, 672, and 866 nm.MISR's unique blend of directional and spectral data allows aerosol retrieval algorithms to be used that do not depend on the explicit radiometric surface properties.Therefore, MISR can retrieve aerosol properties even over highly reflective surfaces (e.g.deserts) (Diner et al, 1998;Martonchik et al., 2004).
The MISR on board the Terra platform passes over the equator at ∼10:45 LT during daytime (Diner et al, 2001).In this study, the available data from 2003 to 2009 are used to evaluate the modeling results during the same period.Model results are sampled in the same overpass time as Terra when comparing with MISR retrievals.

IMPROVE
The Interagency Monitoring for Protected Visual Environments (IMPROVE) program was started in 1985 by the US federal agencies EPA, National Park Services, Department of Agriculture-Forest Service, and other land management agencies as a part of the EPA Regional Haze program (Malm et al., 1994).There are more than 100 IMPROVE sites in the US regularly measuring aerosol chemical compositions as well IMPROVE).The mass of dust is estimated from the measured elements associated predominantly with soil, including Ti, Al, Si, Ca, and Fe.The uncertainties of estimated values results from the range and variation of mineral composition from a variety of soil types (see Malm et al., 1994 for details).Since there is no IMPROVE sites over the GBD, fine dust (with diameter less than 2.5 µm) mass concentrations from the 7 sites in the region (110 • N, mostly locate in the Arizona state, Fig. 1) near the GBD during 1995-2009 are used to evaluate model results.

TRMM
The 3B-43 product of the Tropical Rainfall Measuring Mission (TRMM) data is used.
The purpose of algorithm 3B-43 is to produce TRMM merged high quality (HQ)/infrared (IR) precipitation and root-mean-square (RMS) precipitation-error estimates.These gridded estimates are on a monthly temporal resolution and a 0.25 × 0.25 degree spatial resolution in a global belt extending from 50 • S to 50 • N latitude.Detailed description of the dataset can be found in Huffman et al. (2001Huffman et al. ( , 2007)).In this study, the available data from 2000 to 2009 are used to evaluate the modeling results during the same period.

PERSIANN
Precipitation Estimation from Remote Sensing Information using Artificial Neural Network (PERSIANN) provides hourly precipitation data on a 0.25 × 0.25 degree spatial resolution in a global belt extending from 60 • S to 60 • N latitude (Hsu et al., 1999;Sorooshian et al., 2000).In the PERSIANN system, a neural network trained by precipitation from the TRMM Microwave Imager (TMI) and other microwave measurements (Hsu et al., 1997(Hsu et al., , 1999) ) was used to estimate 30-min precipitation rates from infrared (IR) and visible imagery from geostationary satellites.In this study, the available data from 2000 to 2009 are used to evaluate the modeling results during the same period.Introduction

Conclusions References
Tables Figures

Back Close
Full Figure 4 shows the WRF-Chem simulated 24-h averaged spatial distribution of clearsky dust SW, LW, and net (SW + LW) radiative forcing at the top and bottom of the atmosphere and in the atmosphere during the summer period over North America.The dust radiative forcing is derived from the WRF-Chem simulations with and without dust.At the top of atmosphere (TOA), the dust SW radiative forcing is negative throughout the domain.The pattern is consistent with the spatial distribution of dust mass concentration (Fig. 2) with the highest forcing of −1.25 W m −2 over the desert region.The dust LW radiative forcing at TOA is mostly close to zero with some slightly positive values over the domain.Overall, the simulated dust net (SW+LW) radiative forcing at TOA has a cooling effect of −0.1 W m −2 on domain average with a maximum of −1.15 W m −2 over the GBD.In the atmosphere, dust absorbs SW and warms the atmosphere.Dust LW effect cools the atmosphere.The dust LW forcing is higher over the desert because of the hotter desert surface and stronger interaction between LW and coarse mode dust particles that are spatially confined to the desert due to their shorter lifetimes and distance being carried.In the atmosphere, dust produces a net America.Both TRMM and PERSIANN show a typical monthly monsoonal evolution.Monsoon onset occurs in June in most of western Mexico south of ∼30 • N and reaches the international border and then the southwestern US in July and August.Rainfall intensifies strongly and rapidly in July, with maximum rainfall rates along the western foothills of the Sierra Madre Occidental in northwestern Mexico, consistent with previous studies (e.g.Gutzler, 2009).Although, TRMM observes the spatial variability of the NAM precipitation, it underestimates the NAM precipitation compared to the PER-SIANN data.It has been found that TRMM has low biases for NAM precipitation (Gutzler, 2009).The model successfully simulates the observed spatial distribution of the monthly monsoonal evolution with the northward migration of the rain-belt.The model reproduces the monsoonal precipitation maximum across northwest Mexico.To further investigate the NAM feature, the NAM region is defined in the black box (20 following previous studies (Gutzler, 2004(Gutzler, , 2009)).
Figure 5b shows the monthly and seasonal mean spatial distribution of precipitation over the NAM region from the WRF-Chem simulations.It has been stated that the NA monsoon is most visible in Arizona and New Mexico states of the US and the northwest of Mexico (Adams and Comire, 1997).Four sub-regions (black boxes) are further defined for time series analysis within the NAM region: DESERT, CORE, AZNM, and TXNM.The DESERT region is defined to cover the dust source region with a very low precipitation (<10 mm month −1 ).The CORE region is defined to cover the tongue of the NAM rain-belt that receives the highest precipitation among the four sub-regions.The AZNM region is defined as the region over Arizona and New Mexico, which are the two states that are most affected by the NAM in the southwestern US and receive ∼50 % of their annual precipitation during the NAM season (Adams and Comrie, 1997).The TXNM region is defined adjacent to the east of the CORE and AZNM regions.
To examine lower tropospheric moisture transports in the NAM region, longitudeheight cross-section of zonal and meridional moisture fluxes are investigated.Figure 6 shows the WRF-Chem simulated monthly mean cross-section of zonal circulation and meridional moisture fluxes at 29 • N from 115 • W to 105 • W (Fig. 5a) during the NAM Figures

Back Close
Full The low-level southerly moisture fluxes increase from 10 g kg −1 m s −1 to 50 g kg −1 m s −1 from June to August.On the other side, the eastern slopes of the Sierra Madre Occidental, the model also simulates the near-surface southerly moisture fluxes from Gulf of Mexico.Corresponding to the two southerly jets confined to the lower atmosphere, the model also simulates a northerly moist jet at 700-800 hPa at ∼109 • W near the top of the Sierra Madre Occidental.Besides the meridional jets, the zonal moisture fluxes show a convergence at the western slope of the Sierra Madre Occidental, which is consistent with the maximum NAM precipitation over the region.The convergence is stronger in July-August than June.increases precipitation by 6 mm month −1 (12 %) in July and 3 mm month −1 (6 %) in August.Over the TXNM region, dust increases precipitation by 12 mm month −1 (25 %) in July and 15 mm month −1 (40 %) in August.Although the monthly variation of precipitation over the TXNM region corresponds least with the NAM migration, the dustinduced precipitation change increases most significantly along with the NAM migration.From the spatial distribution of dust-induced precipitation, statistically significant impact of dust is mainly found within the AZNM and TXNM regions (not shown).The dust-induced precipitation change over the two regions exceeds the 90 % statistical significance level based on the Student's t-test.

Dust impact on precipitation and circulation
Dust affects precipitation mainly through changes in the surface energy budget and atmospheric diabatic heating profiles via direct radiative forcing and changes in cloud property via semi-direct and indirect radiative forcing.Locally, dust aerosol cools the surface and warms the atmosphere during daytime, hence increases atmospheric stability (e.g.Zhao et al., 2011).However, over the NAM region, since most of the dust concentrates over the DESERT region, the dust direct radiative forcing over the three sub-regions (CORE, AZNM, and TXNM) is relatively small (Fig. 4).In addition, dust semi-direct and indirect effect is also limited over the three sub-regions because dust is mostly confined to its source region and there is no treatment of dust as ice nuclei (IN) in our simulations.However, the heated air by dust over the DESERT region could remotely affect the overall monsoonal circulations (e.g.Lau et al, 2006Lau et al, , 2009)).GBD dust-induced heating, particularly in August, resulting in a shift of moisture convergence to the east of the Sierra Madre Occidental.Overall, the GBD dust helps to draw in warm and moist air from the Gulfs of California and Mexico, resulting in a general increase in NAM precipitation over all three sub-regions.
The TXNM region is taken as an example to further illustrate the correlation between dust-induced diabatic heating and the NAM precipitation change.Figure 9 shows the monthly daytime dust-induced diabatic heating change over the DESERT region along with the dust-induced precipitation change over the TXNM region for May-August, 1995-2009 from the WRF-Chem simulations.In May, dust-induced atmospheric diabatic heating over the DESERT region is small.It increases significantly from 0.05 K day −1 to 0.5 K day −1 from May to July-August.The dust heating effect reaches a maximum in July.Corresponding to the dust heating in the lower atmosphere (below 800 hPa), the dust-induced precipitation change over the TXNM region gradually increases from June to July when dust heating in the lower atmosphere over the DESERT region is largest.

Summary and conclusions
In this study, the radiative forcing of GBD dust and its feedback to the monsoon precipitation and circulation over the NAM region are investigated using the WRF-Chem model.over the source region.In the atmosphere, the GBD dust absorbs SW and warms the atmosphere by 0.05 W m −2 on domain average and up to 0.6 W m −2 over the source region.At the surface, the GBD dust has a cooling effect of −0.15 W m −2 on domain average and −1.38 W m −2 on maximum over the source region.Overall, the net (LW + SW) radiative effect of GBD dust is atmospheric warming, surface cooling, and negative TOA forcing over North America.
The WRF-Chem successfully simulates the observed spatial distribution of the monthly monsoonal evolution with the northward migration of the rain-belt and properly resolves the low level southerly jet along the western foothills of the Sierra Madre Occidental, extending from the surface to 900 hPa.The GBD dust affects the NAM precipitation mostly over three NAM sub-regions: CORE, AZNM, and TXNM in July and August.Dust increases precipitation by 3 mm month −1 (1 %), 6 mm month −1 (12 %), and 12 mm month −1 (25 %) in July and 12 mm month −1 (6 %), 3 mm month −1 (6 %), and 15 mm month −1 (40 %) in August over the CORE, AZNM, and TXNM regions, respectively.The dust-induced precipitation change over the AZNM and TXNM regions exceeds the 90 % statistical significance level based on the Student's t-test.The dust-induced precipitation change is facilitated by the dust-heated air over the DESERT region that helps to draw warm and moist air from the Gulfs of California and Mexico.In July and August, the southerly moisture fluxes at lower atmosphere are strengthened by ∼10 % on both western and eastern slopes of the Sierra Madre Occidental.In the meantime, the easterly moisture fluxes are also strengthened, resulting in an eastward shift of moisture convergence, leading to the largest precipitation increase in the TXNM region.The dust-induced precipitation changes are highly correlated with the dust-induced lower atmospheric (below 800 hPa) diabatic heating (up to 0.5 K day −1 in July and August) over the DESERT region.This mechanism of remote impact of dust on the overall monsoonal circulations is similar as that proposed by Lau et al. (2006Lau et al. ( , 2009) ) over Asia and West Africa.
Compared to the previous studies on the Saharan dust impact on the West Africa Monsoon (WAM) system (e.g.Lau et al. 2009), the GBD dust plays a similar role as a Figures "heat pump" to strengthen the monsoon circulation and re-distributes the precipitation.
Although the GBD dust forcing over the southwestern US is much smaller than the Saharan dust over West Africa, it still has a significant impact on the NAM circulation and precipitation, because the diabatic heating that drives the NAM is also much smaller than that of WAM.For example, residual diagnosis of NAM diabatic heating based on two global reanalysis data suggests a maximum averaged low level heating no more than 1.5 K day −1 in the Sierra Madre Occidental (Barlow et al., 1998) (1-2 K day −1 at 700 hPa simulated in this study).Therefore the dust heating of 0.5 K day −1 estimated from our simulations could produce significant effects on the NAM circulation depending on its location and vertical distribution.However, the diurnal cycle of NAM precipitation is not changed significantly by the GBD dust, which is different from the role of Saharan dust found by previous studies (e.g.Kim et al., 2010;Zhao et al., 2011).The negligible GBD dust impact on the diurnal cycle of NAM precipitation is related to the low dust concentration and therefore negligible impact on atmospheric stability over the precipitation region.Overall, this study highlights the interaction between the GBD dust and NAM system and motivates further investigation of possible dust feedback on the monsoon under climate change and the projected mega-drought conditions over the southwestern US in the 21st century.Finally, several caveats should be noted in this study.First, in the absence of measurements of GBD dust optical property, this study uses the same refractive index of dust as our previous studies over West Africa (Zhao et al., 2010(Zhao et al., , 2011)).The imaginary part of dust refractive index is uncertain (e.g.Balkanski et al., 2007;McConnell et al., 2010;Zhao et al., 2011) and may change the dust-induced atmospheric diabatic heating, which triggers the NAM circulation and precipitation change.In order to constrain the GBD dust absorption, measurements of the GBD dust optical properties are required in the future.Second, although the dust impact on diabatic heating and hence the NAM circulation is evident, the dust-induced changes of NAM precipitation may be sensitive to the cloud microphysics and convective schemes used in the model, which is beyond the scope of this study.We note, however, that the simulations reported Introduction

Conclusions References
Tables Figures

Back Close
Full in this study reproduced the seasonal cycle and spatial distribution of NAM precipitation rather well, hence lending some confidence in simulated response to changes in dust diabatic heating.Third, although the GBD dust-induced atmospheric heating influences the NAM circulation and precipitation, the effect cannot be fully simulated using a regional modeling framework because the same SST lower boundary conditions and atmospheric lateral boundary conditions are used in the simulations with and without dust.However, this limitation may not have large impacts because the GBD dust is rather spatially confined so its impacts on large-scale circulation may not extend well beyond the boundaries of our simulation domain.On the other hand, NAM precipitation is generally not well simulated by GCM's because higher spatial resolution is needed to represent the complex landscape and its interactions with regional and large-scale circulation over the NAM region.Last, dust indirect effects on convective cloud and precipitation are not represented in the current convective microphysical parameterization, although convective precipitation dominates during the NAM season (generating over 90 % of precipitation in the simulations).In addition, the dust aerosol cannot be  B., Clark, D., Frouin, R., Halthore, R., Karneli, A., O'Neill, N., Pietras, C., Pinker, R., Voss, K., and Zibordi, G.: An emerging ground-based aerosol climatology: Aerosol optical depth from AERONET, J. Geophys.Res., 106, 12067-12097, 2001.Hsu, K., Gao, X., Sorooshian, S., and Gupta, H. V.: Precipitation estimation from remotely sensed information using artificial neural networks.,J. Appl. Meteor., 36, 1176-1190, 1997.Introduction

Conclusions References
Tables Figures

Back Close
Full The NAM system is characterized by a combination of seasonally warm land surfaces and atmospheric moisture from nearby Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | ) investigates the potential dust effects on circulation and precipitation during the NAM season.Dust effects on cloud and precipitation through aerosol indirect effects are not considered in this study as the primary dust source region is generally removed from the monsoon core region where clouds are prominent (Sect.4.1).The paper is organized as follows.Sections 2 and 3 detail the WRF-Chem model and the measurements used in this study.The SW and LW radiative forcings of the GBD dust and their impact on the NAM circulation and precipitation are analyzed in Sect. 4. The paper concludes in Sect Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | as fine and coarse mode aerosol mass concentrations (http://vista.cira.colostate.edu/Discussion Paper | Discussion Paper | Discussion Paper | Figure1shows the monthly mean spatial distribution of 10-m wind circulation from NARR reanalysis data and WRF-Chem simulations in June, July, August, and summer averaged for 1995-2009.The dust emissions from WRF-Chem simulations are also shown.In general, the model well captures the surface wind circulation compared to the NARR reanalysis data.Dust emission occurs mostly over the southwestern US.The monthly variation of dust emission is small during summer.WRF-Chem simulates a total dust emission of 2.4 Tg for the summer season.Figure2shows the monthly mean spatial distribution of wind circulation at 800 hPa from NARR reanalysis data and WRF-Chem simulations in June, July, August, and summer averaged for 1995-2009.The dust mass concentration at 800 hPa from WRF-Chem simulations is also shown.The model simulated wind circulation is consistent with the reanalysis data.In general, dust concentrates near the GBD region in all months.The model simulates higher dust concentration over the source region in July and August than that in June.The dust in

(
Figure5ashows the monthly and seasonal mean spatial distribution of precipitation from TRMM and PERSIANN observations and WRF-Chem simulations over North Discussion Paper | Discussion Paper | Discussion Paper | season, extending across the Gulf of California onto the high topography of the Sierra Madre Occidental.At 29 • N, the Gulf of California extends zonally between 113 • W and 111 • W. The model correctly resolves the low level southerly jet along the western foothills of the Sierra Madre Occidental, extending from surface to 900 hPa.The model clearly exhibits a monthly-mean southerly moisture fluxes maximum near 111 • W in the vicinity of Gulf of California.

Figure 7
Figure 7 shows the monthly precipitation and dust-induced change over the three NAM sub-regions (CORE, AZNM, and TXNM) from WRF-Chem simulations during May-August averaged for 1995-2009.The dust-induced precipitation change is calculated by subtracting the results from the WRF-Chem simulations without dust from the results from the WRF-Chem simulations with dust.In general, the WRF-Chem simulations without dust show that precipitation increases from ∼40 mm month −1 to ∼200 mm month −1 over the CORE region and from ∼24 mm month −1 to ∼65 mm month −1 over the AZNM region during May-August corresponding to the NAM northward migration, while the precipitation over the TXNM region reduces from ∼56 mm month −1 to ∼40 mm month −1 .The dust-induced precipitation change mainly occurs in July and August over all three regions.Over the CORE region, dust increases precipitation by 3 mm month −1 (1 %) in July and 12 mm month −1 (6 %) in August.Over the AZNM region, dust

Figure 8
Figure8shows the monthly mean cross sections of dust-induced change of zonal and meridional moisture fluxes at 29• N from 115• W to 105• W during the NAM season averaged for 1995-2009.Same as Fig.7, the dust-induced circulation change is calculated by subtracting the results from the WRF-Chem simulations without dust from the results from the WRF-Chem simulations with dust.In general, the dust impact on the moisture transport is negligible in June, but it becomes significant in July and August.In July and August, the southerly moisture fluxes are strengthened by 1-6 g kg −1 m s −1 (∼10 %) on both western and eastern slopes of the Sierra Madre Occidental.In the meantime, the easterly moisture fluxes are also strengthened by the To establish statistical significance of dust effects, simulation results for May-August of 1995-2009 are analyzed to capture potential dominant dust effects above interannual variability.Both dust SW and LW direct radiative effect are included.Generally, the model simulates reasonable aerosol mass concentrations and optical depth over the NAM region during summer.Compared to the observations, the model successfully reproduces the magnitude and spatial distribution of NAM precipitation.The WRF-Chem simulates a total GBD dust emission of 2.4 Tg during the NAM season.The GBD dust concentrates near the source region in the lower atmosphere (below 800 hPa) during summer, resulting in negative TOA forcing by up to −1.25 W m Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | activated as ice nuclei (IN) to influence clouds in the Morrison microphysics scheme in this study.Despite the low occurrence of dust over the NAM cloud/precipitation region, episodic intrusions of dust to the CORE region due to anomalous transport could result in important aerosol indirect effects that deserve further study in the future.Discussion Paper | Discussion Paper | Discussion Paper |

Figure 1 10 Fig. 1 .Figure 2 Fig. 2 .Figure 3 Fig. 3 .Figure 4 Fig. 4 .Figure 5a
Figure 1 Spatial distributions of 10-m wind from NARR reanalysis data and WRF-Chem 4 simulations and dust emissions from WRF-Chem simulations in June, July, August, and 5 summer averaged for 1995-2009.The seven upward black triangles (rightmost of the top 6 panel) represent the IMPROVE sites; the one downward black triangle (rightmost of the 7 bottom panel) represents the AERONET site Maricopa in Arizona.8 9 10

Fig. 5a .Figure 5b
Fig. 5a.Spatial distributions of precipitation from TRMM and PERSIANN data and WRF-Chem simulations in June, July, August, and summer averaged for 2000-2009.The region in black box is defined as the North America Monsoon (NAM) region (20 • N-40 • N, 120 • W-95 • W).The black line within the black box represents the cross section at 29 • N from 115 • W to 105 • W.

Fig. 6 .Fig. 7 .Figure 8
Fig. 6.Monthly mean cross sections of zonal circulation and meridional moisture fluxes at 29 • N from • W to 105 • W (Fig. 5a) during the NAM season.The color contour shows the meridional moisture fluxes.The positive values represent southerlies and negative values represent northerlies.

Fig. 8 .Figure 9 Fig. 9 .
Fig. 8. Monthly mean cross sections of dust-induced change of zonal and meridional moisture fluxes at 29 • N from 115 • W to 105 • W during the NAM season.The color contour shows the change in meridional moisture fluxes.The positive values represent increase in southerlies and negative values represent increase in northerlies.