Abstract
The southwestern United States is currently in a multi-decade drought that has developed since a precipitation maximum in the 1980s. While anthropogenic warming has made the drought more severe, it is the decline in winterâspring precipitation that has had a more profound effect on water resources and ecosystems. This precipitation decline is not well understood beyond its attribution to the post-1980 La Niña-like cooling trend in tropical sea surface temperatures, which caused a North Pacific anti-cyclonic atmospheric circulation trend conducive to declining precipitation in the southwestern United States. Using a hierarchy of model simulations, we show that, even under El Niño-like sea surface temperature trends, there is a tendency towards a North Pacific anti-cyclonic circulation trend and declining precipitation in the southwestern United States, counter to the canonical El Niño teleconnection. This unintuitive yet robust circulation change arises from non-additive responses to tropical mean sea surface temperature warming and radiative effects from anthropogenic aerosols. The post-1980 period exhibits the fastest southwestern US soil moisture drying among past and future periods of similar length due to the combination of this forced precipitation decline and anthropogenic warming. While the precipitation trend might reverse due to future projected El Niño-like warming and aerosol emissions reduction, it is unlikely to substantially alleviate the currently projected future drought risk.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
27,99 ⬠/ 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
269,00 ⬠per year
only 22,42 ⬠per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout






Similar content being viewed by others

Data availability
The information to access CESM2-LE and single-forcing simulations can be found at https://www.cesm.ucar.edu/community-projects/lens2/data-sets; Research Data Archive (RDA) at NCAR provides access to ERSSTv5 (ref. 50; https://www.ncei.noaa.gov/) and ERA5 (ref. 52; https://www.ecmwf.int/en/forecasts/datasets); GPCC is accessed through ref. 51.
Code availability
Figures are generated with Python, and maps are generated with the library Cartopy68, including the basemap from Natural Earth69. Codes to generate the results in this study are available via Zenodo at https://doi.org/10.5281/zenodo.14990892 (ref. 70).
Change history
21 October 2025
A Correction to this paper has been published: https://doi.org/10.1038/s41561-025-01774-5
References
Gleick, P. H. Roadmap for sustainable water resources in southwestern North America. Proc. Natl Acad. Sci. USA 107, 21300â21305 (2010).
MedellÃn-Azuara, J. et al. Economic analysis of the 2016 California drought on agriculture. California WaterBlog https://californiawaterblog.com/2016/08/15/economic-analysis-of-the-2016-california-drought-for-agriculture/ (2016).
Prein, A. F., Holland, G. J., Rasmussen, R. M., Clark, M. P. & Tye, M. R. Running dry: the US Southwestâs drift into a drier climate state. Geophys. Res. Lett. 43, 1272â1279 (2016).
Lehner, F., Deser, C., Simpson, I. R. & Terray, L. Attributing the US Southwestâs recent shift into drier conditions. Geophys. Res. Lett. 45, 6251â6261 (2018).
Diffenbaugh, N. S., Swain, D. L. & Touma, D. Anthropogenic warming has increased drought risk in California. Proc. Natl Acad. Sci. USA 112, 3931â3936 (2015).
Williams, A. P. et al. Large contribution from anthropogenic warming to an emerging North American megadrought. Science 368, 314â318 (2020).
Ault, T. R., Mankin, J. S., Cook, B. I. & Smerdon, J. E. Relative impacts of mitigation, temperature, and precipitation on 21st-century megadrought risk in the American Southwest. Sci. Adv. 2, e1600873 (2016).
Juang, C. S. et al. Rapid growth of large forest fires drives the exponential response of annual forestâfire area to aridity in the western United States. Geophys. Res. Lett. 49, e2021GL097131 (2022).
Jacobson, T. W. P. et al. An unexpected decline in spring atmospheric humidity in the interior southwestern United States and implications for forest fires. J. Hydrometeorol. 25, 373â390 (2024).
Lukas, J. & Payton, E. Colorado River Basin Climate and Hydrology: State of the Science (Western Water Assessment, University of Colorado Boulder, 2020).
Delworth, T. L., Zeng, F., Rosati, A., Vecchi, G. A. & Wittenberg, A. T. A link between the hiatus in global warming and North American drought. J. Clim. 28, 3834â3845 (2015).
Polade, S. D., Gershunov, A., Cayan, D. R., Dettinger, M. D. & Pierce, D. W. Precipitation in a warming world: assessing projected hydro-climate changes in California and other Mediterranean climate regions. Sci. Rep. 7, 10783 (2017).
Seager, R. et al. Climate variability and change of Mediterranean-type climates. J. Clim. 32, 2887â2915 (2019).
Shepherd, T. G. Atmospheric circulation as a source of uncertainty in climate change projections. Nat. Geosci. 7, 703â708 (2014).
Schmidt, D. F. & Grise, K. M. The response of local precipitation and sea level pressure to Hadley cell expansion. Geophys. Res. Lett. 44, 10573â10582 (2017).
Kuo, Y., Kim, H. & Lehner, F. Anthropogenic aerosols contribute to the recent decline in precipitation over the US Southwest. Geophys. Res. Lett. 50, e2023GL105389 (2023).
Carrillo, C. M. et al. Megadrought: a series of unfortunate La Niña events? JGR Atmos. 127, e2021JD036376 (2022).
Seager, R. & Hoerling, M. Atmosphere and ocean origins of North American droughts. J. Clim. 27, 4581â4606 (2014).
Allen, R. J. & Luptowitz, R. El Niño-like teleconnection increases California precipitation in response to warming. Nat. Commun. 8, 16055 (2017).
Seager, R. et al. Ocean-forcing of cool season precipitation drives ongoing and future decadal drought in southwestern North America. NPJ Clim. Atmos. Sci. 6, 141 (2023).
Wills, R. C. J., Dong, Y., Proistosecu, C., Armour, K. C. & Battisti, D. S. Systematic climate model biases in the largeâscale patterns of recent seaâsurface temperature and seaâlevel pressure change. Geophys. Res. Lett. 49, e2022GL100011 (2022).
Seager, R., Henderson, N. & Cane, M. Persistent discrepancies between observed and modeled trends in the tropical Pacific Ocean. J. Clim. 35, 4571â4584 (2022).
Coats, S. & Karnauskas, K. B. Are simulated and observed twentieth century tropical Pacific sea surface temperature trends significant relative to internal variability? Geophys. Res. Lett. 44, 9928â9937 (2017).
Dong, L. & Leung, L. R. Winter precipitation changes in California under global warming: contributions of CO2, uniform SST warming, and SST change patterns. Geophys. Res. Lett. 48, e2020GL091736 (2021).
Lehner, F. & Deser, C. Origin, importance, and predictive limits of internal climate variability. Environ. Res. Clim. 2, 023001 (2023).
Dow, W. J., Maycock, A. C., Lofverstrom, M. & Smith, C. J. The effect of anthropogenic aerosols on the Aleutian Low. J. Clim. 34, 1725â1741 (2021).
Kang, J. M., Shaw, T. A. & Sun, L. Anthropogenic aerosols have significantly weakened the regional summertime circulation in the Northern Hemisphere during the satellite era. AGU Adv. 5, e2024AV001318 (2024).
Allen, R. J., Lamarque, J., WatsonâParris, D. & Olivié, D. Assessing California wintertime precipitation responses to various climate drivers. J. Geophys. Res. Atmos. 125, e2019JD031736 (2020).
Wang, Y., Hu, K., Huang, G. & Tao, W. Asymmetric impacts of El Niño and La Niña on the PacificâNorth American teleconnection pattern: the role of subtropical jet stream. Environ. Res. Lett. 16, 114040 (2021).
Kushnir, Y., Seager, R., Ting, M., Naik, N. & Nakamura, J. Mechanisms of Tropical Atlantic SST influence on North American precipitation variability. J. Clim. 23, 5610â5628 (2010).
Xu, M., Zhan, R. & Zhao, J. Distinct responses of tropical cyclone activity to spatio-uniform and nonuniform SST warming patterns. Environ. Res. Lett. 19, 064020 (2024).
Lu, J., Vecchi, G. A. & Reichler, T. Expansion of the Hadley cell under global warming. Geophys. Res. Lett. 34, 2006GL028443 (2007).
Baek, S. H. et al. Precipitation, temperature, and teleconnection signals across the combined North American, Monsoon Asia, and Old World drought atlases. J. Clim. 30, 7141â7155 (2017).
Zeppetello, L. R. V., Zhang, L. N., Battisti, D. S. & Laguë, M. M. How much does landâatmosphere coupling influence summertime temperature variability in the western United States? J. Clim. 37, 3457â3478 (2024).
Alessi, M. J. & Rugenstein, M. Potential nearâterm wetting of the southwestern United States if the eastern and Central Pacific cooling trend reverses. Geophys. Res. Lett. 51, e2024GL108292 (2024).
Persad, G. G., Samset, B. H. & Wilcox, L. J. Aerosols must be part of climate risk assessments. Nature 611, 662â664 (2022).
Qiu, W., Collins, M., Scaife, A. A. & Santoso, A. Tropical Pacific trends explain the discrepancy between observed and modelled rainfall change over the Americas. NPJ Clim. Atmos. Sci. 7, 201 (2024).
Chung, E.-S. et al. Reconciling opposing Walker circulation trends in observations and model projections. Nat. Clim. Change 9, 405â412 (2019).
Seager, R. et al. Strengthening tropical Pacific zonal sea surface temperature gradient consistent with rising greenhouse gases. Nat. Clim. Change 9, 517â522 (2019).
Heede, U. K. & Fedorov, A. V. Eastern equatorial Pacific warming delayed by aerosols and thermostat response to CO2 increase. Nat. Clim. Change 11, 696â703 (2021).
Hwang, Y.-T., Xie, S.-P., Chen, P.-J., Tseng, H.-Y. & Deser, C. Contribution of anthropogenic aerosols to persistent La Niña-like conditions in the early 21st century. Proc. Natl Acad. Sci. USA 121, e2315124121 (2024).
Hartmann, D. L. The Antarctic ozone hole and the pattern effect on climate sensitivity. Proc. Natl Acad. Sci. USA 119, e2207889119 (2022).
Kim, H., Kang, S. M., Kay, J. E. & Xie, S.-P. Subtropical clouds key to Southern Ocean teleconnections to the tropical Pacific. Proc. Natl Acad. Sci. USA 119, e2200514119 (2022).
Dong, Y., Armour, K. C., Battisti, D. S. & Blanchard-Wrigglesworth, E. Two-way teleconnections between the Southern Ocean and the tropical Pacific via a dynamic feedback. J. Clim. 35, 6267â6282 (2022).
Kohyama, T., Hartmann, D. L. & Battisti, D. S. La Niñaâlike mean-state response to global warming and potential oceanic roles. J. Clim. 30, 4207â4225 (2017).
Shin, S.-I. & Sardeshmukh, P. D. Critical influence of the pattern of tropical ocean warming on remote climate trends. Clim. Dyn. 36, 1577â1591 (2011).
Watanabe, M., Iwakiri, T., Dong, Y. & Kang, S. M. Two competing drivers of the recent Walker circulation trend. Geophys. Res. Lett. 50, e2023GL105332 (2023).
Chen, M. et al. Why do DJF 2023/24 upperâlevel 200âhPa geopotential height forecasts look different from the expected El Niño response? Geophys. Res. Lett. 51, e2024GL108946 (2024).
Schmidt, D. F. & Grise, K. M. Impacts of subtropical highs on summertime precipitation in North America. J. Geophys. Res. Atmos. 124, 11188â11204 (2019).
Huang, B. et al. Extended reconstructed sea surface temperature, version 5 (ERSSTv5): upgrades, validations, and intercomparisons. J. Clim. 30, 8179â8205 (2017).
Schneider, U., Hänsel, S., Finger, P., Rustemeier, E. & Ziese, M. GPCC Full Data Monthly Product Version 2022 at 1.0°: Monthly Land-Surface Precipitation from Rain-Gauges Built on GTS-Based and Historical Data Global Precipitation Climatology Centre (GPCC, http://gpcc.dwd.de/) at Deutscher Wetterdienst https://doi.org/10.5676/DWD_GPCC/FD_M_V2022_100 (2022)
Hersbach, H. et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 146, 1999â2049 (2020).
Martens, B. et al. GLEAM v3: satellite-based land evaporation and root-zone soil moisture. Geosci. Model Dev. 10, 1903â1925 (2017).
Rohde, R. A. & Hausfather, Z. The Berkeley Earth Land/Ocean Temperature Record. Earth Syst. Sci. Data 12, 3469â3479 (2020).
Rodgers, K. B. et al. Ubiquity of human-induced changes in climate variability. Earth Syst. Dynam. 12, 1393â1411 (2021).
Simpson, I. R. et al. The CESM2 single-forcing large ensemble and comparison to CESM1: implications for experimental design. J. Clim. 36, 5687â5711 (2023).
Danabasoglu, G. et al. The Community Earth System Model Version 2 (CESM2). J. Adv. Model. Earth Syst. 12, e2019MS001916 (2020).
Cook, B. I. et al. Uncertainties, limits, and benefits of climate change mitigation for soil moisture drought in southwestern North America. Earths Future 9, e2021EF002014 (2021).
Ziehn, T. et al. The Australian Earth System Model: ACCESS-ESM1.5. J. South. Hemisph. Earth Syst. Sci. 70, 193â214 (2020).
Eyring, V. et al. Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geosci. Model Dev. 9, 1937â1958 (2016).
Shin, S.-I., Sardeshmukh, P. D., Newman, M., Penland, C. & Alexander, M. A. Impact of annual cycle on ENSO variability and predictability. J. Clim. 34, 171â193 (2021).
Balmaseda, M. A., Mogensen, K. & Weaver, A. T. Evaluation of the ECMWF ocean reanalysis system ORAS4. Q. J. R. Meteorol. Soc. 139, 1132â1161 (2013).
Penland, C. & Matrosova, L. Studies of El Niño and interdecadal variability in tropical sea surface temperatures using a nonnormal filter. J. Clim. 19, 5796â5815 (2006).
DeRepentigny, P. Enhanced simulated early 21st century Arctic sea ice loss due to CMIP6 biomass burning emissions. Sci. Adv. 8, eabo2405 (2022).
Fasullo, J. T. et al. Spurious late historicalâera warming in CESM2 driven by prescribed biomass burning emissions. Geophys. Res. Lett. 49, e2021GL097420 (2022).
Yang, W., Hsieh, T.-L. & Vecchi, G. A. Hurricane annual cycle controlled by both seeds and genesis probability. Proc. Natl Acad. Sci. USA 118, e2108397118 (2021).
Hsieh, T., Yang, W., Vecchi, G. A. & Zhao, M. Model spread in the tropical cyclone frequency and seed propensity index across global warming and ENSOâlike perturbations. Geophys. Res. Lett. 49, e2021GL097157 (2022).
Elson, P. et al. SciTools/cartopy v.0.21.1. Zenodo https://doi.org/10.5281/zenodo.7430317 (2022).
Made with Natural Earth. Free vector and raster map data @ naturalearthdata.com
Kuo, Y.-N. & Shin, S.-I. ynkuo/kuo24_lim-toga_swus: v2. Zenodo https://doi.org/10.5281/zenodo.14990892 (2025).
Acknowledgements
We thank W. Yang, A. Hoell and H. Kim for helpful discussions and comments. Y.-N.K. and F.L. were supported by NOAA MAPP award NA21OAR4310349. F.L. and J.M.A. acknowledge support from the US Department of Energy, Office of Science, Office of Biological & Environmental Research (BER), Regional and Global Model Analysis (RGMA) component of the Earth and Environmental System Modeling Program under award number DE-SC0022070. F.L. also acknowledges support from National Science Foundation (NSF) IA 1947282. I.R.S., C.D. and A.S.P. acknowledge funding from the NSF National Center for Atmospheric Research (NCAR), which is a major facility sponsored by the National Science Foundation under cooperative agreement no. 1852977. I.R.S. also acknowledges support from NOAA MAPP awards NA20OAR4310413 and NA23OAR4310634. M.N. and S.-I.S. acknowledge support from NOAA Cooperative Agreement NA22OAR4320151. S.W. and J.M.A. were supported by the Australian Research Council Centre of Excellence for Climate Extremes (grant CE170100023). The ACCESS simulations were undertaken with the assistance of resources and services from the National Computational Infrastructure, which is supported by the Australian Government. We also acknowledge the CESM Climate Variability and Change working group for making available the regular CESM TOGA, AAER and xAAER simulations used in this work. Simulations were conducted on UCARâs supercomputers Cheyenne (https://doi.org/10.5065/D6RX99HX) and Derecho (https://doi.org/10.5065/qx9a-pg09), operated by NCARâs Computational and Information Systems Laboratory.
Author information
Authors and Affiliations
Contributions
Y.-N.K., F.L., C.D., M.N. and I.R.S. conceptualized this study and the experimental set-ups. M.N. and S.-I.S. generated the synthetic sea surface temperature (SST) with the linear inverse model. Y.-N.K., A.S.P. and S.W. conducted the prescribed SST experiments (LIM-TOGA and F2000climo sensitivity experiments). Y.-N.K. performed the data analyses and visualizations. Y.-N.K., F.L., C.D., I.R.S., S.-I.S. and J.M.A. wrote and edited the initial versions of the manuscript. All the authors participated in discussions on interpreting the results and contributed to the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Geoscience thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Thomas Richardson and James Super, in collaboration with the Nature Geoscience team.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 DJFMAM sea level pressure (psl) and precipitation (pr) trends over 1980-2014.
Trends from a, global ocean global atmosphere (GOGA; that is, standard AMIP experiment) simulation (20; including both CAM6/CLM5 and UM7.3/CABLE), b, the difference between GOGA (20) and TOGA (20), c, GOGA from CAM6/CLM5 (10), d, GOGA from UM7.3/CABLE (10), e, TOGA from CAM6/CLM5 (10), f, TOGA from UM7.3/CABLE (10). The Hatching/Stippling in (a, c, d, e, f) indicates 67% of the ensemble members agree with the sign of ensemble mean trend. The Hatching/Stippling in b indicates the difference between GOGA and TOGA is significant at 95% confidence level (t test).
Extended Data Fig. 2 Same as Figs. 2c-2h but from ACCESS-ESM1.5.
a,c,e, El Niño-like SST trends and associated psl and pr responses in ACCESS-ESM1.5 from PiCtrl (22) (a), Hist (4) (c) and LIM-TOGA (10) (e). b,d,f, Same as a,c,e, but for La Niña-like SST trends. Hatching/stippling indicates psl/pr trends where 67% of the ensemble members agree with the sign of the ensemble mean. Contours represent sea-level pressure trends starting from ±0.1âhPaâdecadeâ1 (red solid lines for positive values; blue dashed lines for negative values), with contour spacing 0.3âhPaâdecadeâ1.
Extended Data Fig. 3 The residual DJFMAM psl, pr trends in fully coupled large ensemble and single forcing simulations and the sea surface temperature trends in single forcing simulations.
a, The residual trends of psl and pr from the Anthropogenic Aerosols (AAER) and everything-but-anthropogenic aerosols (xAAER) forcing simulations. The residual trends are calculated as the difference between the all forcing large ensemble (Fig. 1g) and the sum of the AAER (Fig. 4b) and xAAER (Fig. 4c). Hatching/stippling indicate psl/pr trends with 67% of the ensemble members agreeing with the sign of the ensemble mean for model simulations. Contours represent sea-level pressure trends starting from ±0.1âhPaâdecadeâ1 (red solid lines for positive values; blue dashed lines for negative values), with contour spacing 0.3âhPaâdecadeâ1. The tropical SST trends from b, AAER, and c, xAAER simulations.
Extended Data Fig. 4 DJFMAM trends of 200 hPa zonal wind (U200), 200 hPa geopotential height (Z200), and sea level pressure (psl) for 1980â2014.
(a, b) anthropogenic aerosol (AAER) simulation, (c, d) radiative forcing only (RF-only) simulation, (e, f) ElNiño-like-only simulation, and (g, h) TOGA: El Niño-like simulation. Black contours in (a, c, e, g) are the DJFMAM climatological U200 from the F2000climo control run, starting from 20âm/s with a contour spacing 10âm/s, and hatching indicates 67% of the ensemble members agree on the change in U200. Red/Blue contours in (b, d, f, h) are psl trends shown in the Main text, and hatching indicates 67% of the ensemble members agree on the change in Z200.
Extended Data Fig. 5 Continued DJFMAM trends of 200 hPa zonal wind (U200), 200 hPa geopotential height (Z200), and sea level pressure (psl) for 1980â2014.
(a, b) residual from ElNiño-like-only and RF-only, (c, d) residual from AAER and xAAER, (e, f) xAAER, and (g, h) ALL. Black contours in (a, c, e, g) are the DJFMAM climatological U200 from the F2000climo control run, starting from 20âm/s with a contour spacing 10âm/s, and hatching indicates 67% of the ensemble members agree on the change in U200. Red/Blue contours in (b, d, f, h) are psl trends shown in the Main text, and hatching indicates 67% of the ensemble members agree on the change in Z200.
Extended Data Fig. 6 DJFMAM lower tropospheric static stability change under 2âK uniform tropical warming.
The lower tropospheric static stability is defined as the difference of the potential temperatures at 700 hPa and at the surface. Stippling indicates 67% of the ensemble members agree with the sign of ensemble mean change.
Supplementary information
Supplementary Information
Supplementary Tables S1, S2.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Kuo, YN., Lehner, F., Simpson, I.R. et al. Recent southwestern US drought exacerbated by anthropogenic aerosols and tropical ocean warming. Nat. Geosci. 18, 578â585 (2025). https://doi.org/10.1038/s41561-025-01728-x
Received:
Accepted:
Published:
Version of record:
Issue date:
DOI: https://doi.org/10.1038/s41561-025-01728-x

