Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Opposite response of strong and moderate positive Indian Ocean Dipole to global warming

Abstract

A strong positive Indian Ocean Dipole (pIOD) induces weather extremes such as the 2019 Australian bushfires and African floods. The impact is influenced by sea surface temperature (SST), yet models disagree on how pIOD SST may respond to greenhouse warming. Here we find increased SST variability of strong pIOD events, with strong equatorial eastern Indian Ocean cool anomalies, but decreased variability of moderate pIOD events, dominated by western warm anomalies. This opposite response is detected in the Coupled Model Inter-comparison Project (CMIP5 and CMIP6) climate models that simulate the two pIOD regimes. Under greenhouse warming, the lower troposphere warms faster than the surface, limiting Ekman pumping that drives the moderate pIOD warm anomalies; however, faster surface warming in the equatorial western region favours atmospheric convection in the west, strengthening equatorial nonlinear advection that forces the strong pIOD cool anomalies. Climate extremes seen in 2019 are therefore likely to occur more frequently under greenhouse warming.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Identification of observed strong and moderate pIOD events.
Fig. 2: Selection of models based on pIOD nonlinearity.
Fig. 3: Projected opposite response of moderate and strong pIOD SST.
Fig. 4: Mechanisms for the projected opposite change in moderate and strong pIOD.

Similar content being viewed by others

Data availability

Data relevant to the paper can be downloaded from websites listed below:

OISST v2 at https://www.esrl.noaa.gov/psd/data/gridded/data.noaa.oisst.v2.html;

ORA-S5 at https://icdc.cen.uni-hamburg.de/daten/reanalysis-ocean/easy-init-ocean/ecmwf-oras5.html;

SODA3.3.1 at https://www2.atmos.umd.edu/~ocean/index_files/soda3.3.1_mn_download.htm;

GODAS at https://www.esrl.noaa.gov/psd/data/gridded/data.godas.html;

CMIP5 database at http://www.ipcc-data.org/sim/gcm_monthly/AR5/;

CMIP6 database at https://esgf-node.llnl.gov/projects/cmip6/;

19 CMIP6 models are used in this study, including: MRI-ESM2–0 (refs. 40,41); CNRM-ESM2–1 (refs. 42,43); CNRM-CM6–1 (refs. 44,45); BCC-CSM2-MR (refs. 46,47); EC-Earth3-Veg (refs. 48,49); NorESM2-LM (refs. 50,51); MIROC6 (refs. 52,53); CESM2 (refs. 54,55); CAMS-CSM1-0 (refs. 56,57); CESM2-WACCM (refs. 58,59); MIROC-ES2L (refs. 60,61); IPSL-CM6A-LR (refs. 62,63); INM-CM4-8 (refs. 64,65); INM-CM5-0 (refs. 66,67); CanESM5 (refs. 68,69); MPI-ESM1-2-LR (refs. 70,71); UKESM1-0-LL (refs. 72,73); MCM-UA-1-0 (refs. 74,75); GFDL-CM4 (refs. 76,77).

Code availability

Codes for calculating EOF, the parameter |α| can be downloaded from https://drive.google.com/open?id=1d2R8wKpFNW-vMIfoJsbqIGPIBd9Z_8rj.

References

  1. Saji, N. H., Goswami, B. N., Vinayachandran, P. N. & Yamagata, T. A dipole mode in the tropical Indian Ocean. Nature 401, 360–363 (1999).

    CAS  Google Scholar 

  2. Webster, P. J., Moore, A. M., Loschnigg, J. P. & Leben, R. R. Coupled ocean–atmosphere dynamics in the Indian Ocean during 1997–98. Nature 401, 356–360 (1999).

    Article  CAS  Google Scholar 

  3. Schott, F. A., Xie, S.-P. & McCreary Jr, J. P. Indian Ocean circulation and climate variability. Rev. Geophys. https://doi.org/10.1029/2007RG000245 (2009).

  4. Ummenhofer, C. C. et al. What causes southeast Australia’s worst droughts? Geophys. Res. Lett. https://doi.org/10.1029/2008GL036801 (2009).

  5. Ashok, K., Guan, Z. & Yamagata, T. Influence of the Indian Ocean Dipole on the Australian winter rainfall. Geophys. Res. Lett. https://doi.org/10.1029/2003GL017926 (2003).

  6. Cai, W., Cowan, T. & Raupach, M. Positive Indian Ocean Dipole events precondition southeast Australia bushfires. Geophys. Res. Lett. https://doi.org/10.1029/2009GL039902 (2009).

  7. Abram, N. J. et al. Seasonal characteristics of the Indian Ocean Dipole during the Holocene epoch. Nature 445, 299–302 (2007).

    Article  CAS  Google Scholar 

  8. Abram, N. J., Gagan, M. K., Cole, J. E., Hantoro, W. S. & Mudelsee, M. Recent intensification of tropical climate variability in the Indian Ocean. Nat. Geosci. 1, 849–853 (2008).

    Article  CAS  Google Scholar 

  9. Behera, S. K. et al. Paramount impact of the Indian Ocean Dipole on the East African short rains: A CGCM study. J. Clim. 18, 4514–4530 (2005).

    Article  Google Scholar 

  10. Hashizume, M., Chaves, L. F. & Minakawa, N. Indian Ocean Dipole drives malaria resurgence in East African highlands. Sci. Rep. 2, 269 (2012).

    Article  Google Scholar 

  11. Abram, N. J., Gagan, M. K., McCulloch, M. T., Chappell, J. & Hantoro, W. S. Coral Reef Death during the 1997 Indian Ocean Dipole Linked to Indonesian Wildfires. Science 301, 952–955 (2003).

    Article  CAS  Google Scholar 

  12. Page, S. E. et al. The amount of carbon released from peat and forest fires in Indonesia during 1997. Nature 420, 61–65 (2002).

    Article  CAS  Google Scholar 

  13. Wang, G., Cai, W., Yang, K., Santoso, A. & Yamagata, T. A unique feature of the 2019 extreme positive Indian Ocean Dipole event. Geophys. Res. Lett. https://doi.org/10.1029/2020GL088615 (2020).

  14. Richards, L. Brew, N. & Smith, L. 2019–20 Australian Bushfires—Frequently Asked Questions: a Quick Guide (Parliament of Australia, 2020); https://www.aph.gov.au/About_Parliament/Parliamentary_Departments/Parliamentary_Library/pubs/rp/rp1920/Quick_Guides/AustralianBushfires

  15. Cai, W. et al. Projected response of the Indian Ocean Dipole to greenhouse warming. Nat. Geosci. 6, 999–1007 (2013).

    Article  CAS  Google Scholar 

  16. Zheng, X.-T. et al. Indian Ocean dipole response to global warming in the CMIP5 multimodel ensemble. J. Clim. 26, 6067–6080 (2013).

    Article  Google Scholar 

  17. Hui, C. & Zheng, X.-T. Uncertainty in Indian Ocean Dipole response to global warming: the role of internal variability. Clim. Dyn. 51, 3597–3611 (2018).

    Article  Google Scholar 

  18. Vecchi, G. A. & Soden, B. J. Global warming and the weakening of the tropical circulation. J. Clim. 20, 4316–4340 (2007).

    Article  Google Scholar 

  19. Xie, S.-P. et al. Global warming pattern formation: sea surface temperature and rainfall. J. Clim. 23, 966–986 (2010).

    Article  Google Scholar 

  20. Cai, W. et al. Increased frequency of extreme Indian Ocean Dipole events due to greenhouse warming. Nature 510, 254–258 (2014).

    Article  CAS  Google Scholar 

  21. Weller, E. et al. More-frequent extreme northward shifts of eastern Indian Ocean tropical convergence under greenhouse warming. Sci. Rep. 4, 6087 (2014).

    Article  CAS  Google Scholar 

  22. Du, Y., Cai, W. & Wu, Y. A new type of the Indian Ocean Dipole since the mid-1970s. J. Clim. 26, 959–972 (2013).

    Article  Google Scholar 

  23. Anil, N., Kumar, R., R., M., Sajeev, R. & Saji, P. K. Role of distinct flavours of IOD events on Indian summer monsoon. Nat. Hazards 82, 1317–1326 (2016).

    Article  Google Scholar 

  24. Verdon-Kidd, D. C. On the classification of different flavours of Indian Ocean Dipole events. Int. J. Climatol. 38, 4924–4937 (2018).

    Article  Google Scholar 

  25. Endo, S. & Tozuka, T. Two flavors of the Indian Ocean Dipole. Clim. Dyn. 46, 3371–3385 (2016).

    Article  Google Scholar 

  26. Du, Y., Xie, S.-P., Huang, G. & Hu, K. Role of air–sea interaction in the long persistence of El Niño–induced north Indian Ocean warming. J. Clim. 22, 2023–2038 (2009).

    Article  Google Scholar 

  27. Yang, K. et al. Oceanic processes in ocean temperature products key to a realistic presentation of positive Indian Ocean Dipole nonlinearity. Geophys. Res. Lett. 46, e2020GL089396 (2020).

    Google Scholar 

  28. Karamperidou, C., Jin, F.-F. & Conroy, J. L. The importance of ENSO nonlinearities in tropical pacific response to external forcing. Clim. Dyn. 49, 2695–2704 (2017).

    Article  Google Scholar 

  29. Cai, W. et al. Increased variability of eastern Pacific El Niño under greenhouse warming. Nature 564, 201–206 (2018).

    Article  CAS  Google Scholar 

  30. Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment design. Bull. Am. Meteorol. Soc. 93, 485–498 (2012).

    Article  Google Scholar 

  31. 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).

    Article  Google Scholar 

  32. Abram, N. J. et al. Palaeoclimate perspectives on the Indian Ocean Dipole. Quat. Sci. Rev. 237, 106302 (2020).

    Article  Google Scholar 

  33. Doi, T., Behera, S. K. & Yamagata, T. Predictability of the super IOD Event in 2019 and its link with El Niño Modoki. Geophys. Res. Lett. 47, e2019GL086713 (2020).

    Google Scholar 

  34. Behringer, D. & Xue, Y. Evaluation of the global ocean data assimilation system at NCEP: the Pacific Ocean. In Eighth Symposium on Integrated Observing and Assimilation Systems for Atmosphere, Oceans, and Land Surface (Americal Meteorological Society, 2004).

  35. Zuo, H., Balmaseda, M. A., Tietsche, S., Mogensen, K. & Mayer, M. The ECMWF operational ensemble reanalysis–analysis system for ocean and sea ice: a description of the system and assessment. Ocean Sci. 15, 779–808 (2019).

    Article  Google Scholar 

  36. Reynolds, R. W., Rayner, N. A., Smith, T. M., Stokes, D. C. & Wang, W. An improved in situ and satellite SST analysis for climate. J. Clim. 15, 1609–1625 (2002).

    Article  Google Scholar 

  37. Carton, J. A., Chepurin, G. A. & Chen, L. SODA3: A new ocean climate reanalysis. J. Clim. 31, 6967–6983 (2018).

    Article  Google Scholar 

  38. Lorenz, E. N. Empirical Orthogonal Functions and Statistical Weather Prediction Statistical Forecast Project Report 1 (Department of Meteorology, Massachusetts Institute of Technology, 1956).

  39. Austin, P. C. & Tu, J. V. Bootstrap methods for developing predictive models. The American Statistician 58, 131–137 (2004).

    Article  Google Scholar 

  40. Yukimoto, S. et al. MRI MRI-ESM2.0 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.6842

  41. Yukimoto, S. et al. MRI MRI-ESM2.0 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.6929

  42. Seferian, R. CNRM-CERFACS CNRM-ESM2-1 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2018); https://doi.org/10.22033/ESGF/CMIP6.4068

  43. Voldoire, A. CNRM-CERFACS CNRM-ESM2-1 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.4226

  44. Voldoire, A. CMIP6 Simulations of the CNRM-CERFACS Based on CNRM-CM6-1 Model for CMIP Experiment Historical Version 20200102 (Earth System Grid Federation, 2018); https://doi.org/10.22033/ESGF/CMIP6.4066

  45. Voldoire, A. CNRM-CERFACS CNRM-CM6-1 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.4224

  46. Wu, T. et al. BCC BCC-CSM2MR Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2018); https://doi.org/10.22033/ESGF/CMIP6.2948

  47. Xin, X. et al. BCC BCC-CSM2MR Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.3050

  48. EC-Earth Consortium (EC-Earth) EC-Earth-Consortium EC-Earth3-Veg Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.4706

  49. EC-Earth Consortium (EC-Earth) EC-Earth-Consortium EC-Earth3-Veg Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.4914

  50. Seland, Ø. et al. NCC NorESM2-LM Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.8036

  51. Seland, Ø. et al. NCC NorESM2-LM Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.8319

  52. Tatebe, H. & Watanabe, M. MIROC MIROC6 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2018); https://doi.org/10.22033/ESGF/CMIP6.5603

  53. Shiogama, H., Abe, M. & Tatebe, H. MIROC MIROC6 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.5771

  54. Danabasoglu, G. NCAR CESM2 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.7627

  55. Danabasoglu, G. NCAR CESM2 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.7768

  56. Rong, X. C. A. M. S. CAMS_CSM1.0 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.9754

  57. Rong, X. CAMS CAMS-CSM1.0 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.11052

  58. Danabasoglu, G. NCAR CESM2-WACCM Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.10071

  59. Danabasoglu, G. NCAR CESM2-WACCM Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.10115

  60. Hajima, T. et al. MIROC MIROC-ES2L Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.5602

  61. Tachiiri, K. et al. MIROC MIROC-ES2L Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.5770

  62. Boucher, O. et al. IPSL IPSL-CM6A-LR Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2018); https://doi.org/10.22033/ESGF/CMIP6.5195

  63. Boucher, O. et al. IPSL IPSL-CM6A-LR Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.5271

  64. Volodin, E. et al. INM INM-CM4-8 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.5069

  65. Volodin, E. et al. INd INM-CM4-8 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.12337

  66. Volodin, E. et al. INM INM-CM5-0 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.5070

  67. Volodin, E. et al. INM INM-CM5-0 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.12338

  68. Swart, N. C. et al. CCCma CanESM5 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.3610

  69. Swart, N. C. et al. CCCma CanESM5 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.3696

  70. Wieners, K.-H. et al. MPI-M MPI-ESM1.2-LR Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.6595

  71. Wieners, K.-H. et al. MPI-M MPI-ESM1.2-LR Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.6705

  72. Tang, Y. et al. MOHC UKESM1.0-LL Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.6113

  73. Good, P. et al. MOHC UKESM1.0-LL Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.6405

  74. Stouffer, R. UA MCM-UA-1-0 Model Output Prepared for CMIP6 CMIP Historical Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.8888

  75. Stouffer, R. UA MCM-UA-1-0 Model Output Prepared for CMIP6 ScenarioMIP ssp585 Version 20200102 (Earth System Grid Federation, 2019); https://doi.org/10.22033/ESGF/CMIP6.13901

  76. Guo, H. et al. NOAA-GFDL GFDL-CM4 Model Output Version 20200102 (Earth System Grid Federation, 2018); https://doi.org/10.22033/ESGF/CMIP6.1402

  77. Guo, H. et al. NOAA-GFDL GFDL-CM4 Model Output Prepared for CMIP6 ScenarioMIP Version 20200102 (Earth System Grid Federation, 2018); https://doi.org/10.22033/ESGF/CMIP6.9242

Download references

Acknowledgements

This work is supported by National Key R&D Program of China (grant no. 2018YFA0605700). L.W. is supported by the National Natural Science Foundation of China (NSFC) projects (grant nos. 41490643, 41490640, U1606402 and 41521091). W.C., G.W., A.S. and B.N. are supported by CSHOR and the Earth System and Climate Change Hub of the Australian Government’s National Environment Science Program. CSHOR is a joint research Center for Southern Hemisphere Oceans Research between QNLM and CSIRO. G.H. is supported by NSFC (grant nos. 41831175 and 91937302), CAS XDA20060501 and COMS2019Q03. K.Y. is supported by China Postdoctoral Science Foundation (grant no. 2018M640168) and a scholarship from China Scholarship Council. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups for producing and making available their model output. For CMIP, the US Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. We are grateful to the reanalysis groups for making the datasets publicly available.

Author information

Authors and Affiliations

Authors

Contributions

W.C. conceived the study and wrote the initial manuscript in discussion with K.Y. K.Y. performed all analysis and generated final figures. All authors contributed to interpreting results, discussion of the associated dynamics and improvement of this paper.

Corresponding authors

Correspondence to Wenju Cai or Lixin Wu.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Climate Change thanks Jasti Chowdary and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Characteristics of strong pIOD events.

a–c, SON SST anomalies (°C, color shading) and wind stress anomalies (N m−2, vectors) for the 1994, 1997, and 2006 strong pIOD event, respectively. The pattern is dominated by strong cooling in the equatorial eastern Indian Ocean. Strong equatorial easterly anomalies extend to the central-western Indian Ocean. The anomalies are referenced to the average over the 1982–2015 period. d–f, SON SST anomalies (in °C) relative to the SST threshold for tropical convection, defined as the SON SST average over 20°S to 20°N, for the 1994, 1997, and 2006 strong pIOD event, respectively. Areas with SST lower than the convection threshold have been masked; there is a large area of suppressed convection in the equatorial eastern Indian Ocean. g–i, Equatorial vertical potential temperature anomalies (in °C) over SON of 1994, 1997, and 2006, respectively. The equatorial vertical potential temperature is calculated as the average between 2.5°S to 2.5°N. Strong equatorial cooling is seen in the eastern Indian Ocean. Vertical potential temperature and wind stress data are from ORA-S5.

Extended Data Fig. 2 Characteristics of moderate pIOD events.

ac, SON SST anomalies (°C, color shading) and wind stress anomalies (N m−2, vectors) for the 1982, 1987, and 2015 moderate pIOD event, respectively. The pattern is dominated by a broad-scale warming over the western Indian Ocean. Easterly anomalies are confined to the equatorial central-eastern Indian Ocean. The anomalies are referenced to the average over the 1982–2015 period. df, SON SST anomalies (in °C) relative to SST threshold for tropical convection, defined as SON SST average over 20°S to 20°N for 1982, 1987, and 2015 moderate pIOD event, respectively. Areas with SST lower than the convection threshold have been masked; there is a small area of suppressed convection in the equatorial southeastern Indian Ocean. gi, Equatorial vertical potential temperature anomalies (in °C) over SON of 1982, 1987, and 2015, respectively. The equatorial vertical potential temperature is calculated as the average between 2.5°S to 2.5°N. There is a general lack of cooling in the equatorial eastern Indian Ocean as cooling occurs only off Sumatra-Java. Vertical potential temperature and wind stress data are from ORA-S5.

Extended Data Fig. 3 Heat budget analysis based on observations.

a, Composite of heat budget terms over the equatorial Indian Ocean (10°S–10°N, 40°E–100°E) in August-September-October (ASO). The uncertainty bar on each composite represents the range over the three strong or three moderate pIOD events. b, c, Composite of ASO temperature tendency during strong and moderate pIOD events. The green box marks the area where the main warm temperature tendency is located over the southwestern Indian Ocean (10°S–0°, 50°E–80°E). d, Relationship between normalized temperature tendency and the Ekman pumping term in ASO averaged over the main warming area in the southwestern Indian ocean (10°S–0°, 50°E–80°E, green box in c), with statistical properties shown. e, Map of ASO correlation coefficients between the Ekman pumping term averaged over 10°S–0°, 50°E–80°E and zonal wind stress (TauX) anomalies. Correlation coefficients of ±0.27, ±0.32, ±0.41 are statistically significant at the 90%, 95%, and 99% confidence levels, respectively. The area indicated by the red box is used to calculate the average zonal wind stress anomalies. f, Time series of normalized ASO Ekman pumping term averaged over the western Indian Ocean (10°S–0°, 50°E–80°E) (red curve) and zonal wind stress (TauX) anomalies averaged over the equatorial Indian Ocean (5°S–5°N, 60°E–100°E, red box in e) (green curve). The Ekman pumping term has been multiplied by −1 for plotting. The data used in this figure are from ORA-S5 for the period of 1979–2018.

Extended Data Fig. 4 Anomaly patterns of strong and moderate pIOD in some selected models.

A composite of SON SST anomalies (°C) for strong pIOD with an S-index>1.5 s.d. and moderate pIOD with an M-index>1.25 s.d. in a CESM1-BGC, b CNRM-CM5, c GFDL-ESM2M, d MPI-ESM-LR. Left panels show strong pIOD and right panels show moderate pIOD. These show that anomaly patterns can be rather different from one model to another.

Extended Data Fig. 5 Uncertainty in projected pIOD SST using the DMI.

a, Comparison of the standard deviation of the DMI in the present-day (1900–1999) and future (2000–2099) climate in 20 models that simulate a nonlinear coefficient α greater than 50% of the observed in the quadratic fit between the first two principal components (PC1 and PC2) from EOF analysis on SON SST anomalies, PC2(t) = α [PC1(t)]2+β PC1(t) + γ. Although a total of 15 out of the 20 selected models simulate lower DMI variability in the future period (red bars) than in the present-day climate (blue bars), the multi-model mean decrease is not statistically significant. The error bar represents the 95% confidence level. Models that simulate an increase in variability are indicated by green circles. Models from CMIP6 are shown in purple. b, Same as a, but for 34 models that simulate an α smaller than 50% of the observed value. A total of 17 out of the 34 models (50%) simulate smaller variability in the DMI in the future period than in the present-day period. There is no inter-model consensus, and the multi-model mean change is not statistically significant.

Extended Data Fig. 6 Projected change in pIOD frequency.

a, Moderate and b, strong pIOD using a threshold of M-index>1.25 s.d. and S-index>1.5 s.d., respectively. Shown is the comparison of frequency (events per 100 years) in the present-day (1900–1999) and future (2000–2099) climate in the 20 selected models, that is, with an α greater than 50% of the observed value. A total of 17 out of the 20 models (85%) simulate a decreased frequency in moderate pIOD events, with a statistically significant decrease of 52% in the multi-model mean, from 9.60 events per 100 years in the present-day (blue bars) to 4.65 events per 100 years (red bars) in the future climate. By contrast, a total of 16 out of the 20 models (80%) simulate an increased frequency in strong pIOD events, with a statistically significant increase of 66% in the multi-model average, from 7.5 events per 100 years in the present-day (blue bars) to 12.45 events per 100 years (red bars) in the future climate. Models that simulate an opposite change to the multi-model mean are indicated by green circles. c, Same as a, but using a threshold of M-index>1.5 s.d. A total of 18 out of the 20 selected models (90%) simulate a decreased frequency in moderate pIOD events, with a statistically significant decrease of 60% in the multi-model ensemble mean. d, Same as b, but using a threshold of S-index>1 s.d. A total of 18 out of the 20 models (90%) simulate an increased frequency in strong pIOD, with a statistically significant increase of 44% in the multi-model ensemble mean. Models from CMIP6 are indicated in purple.

Extended Data Fig. 7 Sensitivity of projected changes in pIOD variability to model selections.

a, b, Comparison of variability of the moderate pIOD (M-index) and strong pIOD (S-index), respectively, in the present-day (1900–1999) and future (2000–2099) climate in 30 models that simulate a nonlinear coefficient α greater than 33.3% of the observed value. In a, a total of 27 out of the 30 models (90%) simulate a decrease in M-index variability in the future period (red bars) from that in the present-day period (blue bars). This leads to a multi-model mean decrease of 17%, statistically significant above the 95% confidence level, as indicated by the error bars. In b, a total of 23 out of the 30 selected models (77%) simulate an increase in S-index variability in the future climate, with a multi-model mean increase of 19%, statistically significant above the 95% confidence level. Models from CMIP6 are indicated in purple. Models that simulate an opposite change to the multi-model mean are indicated by green circles.

Extended Data Fig. 8 Sensitivity of projected changes in pIOD variability to emission scenario.

Shown are results for emission scenario RCP 4.5. a, b, Comparison of variability of the moderate pIOD (M-index) and strong pIOD (S-index), respectively, in the present-day (1900–1999) and future (2000–2099) climate in 15 out of 33 CMIP5 models that simulate a nonlinear coefficient α greater than 50% of the observed value. In a, a total of 14 out of the 15 models (93%) simulate a decrease in M-index variability in the future period (red bars) from that in the present-day period (blue bars). This leads to a multi-model mean decrease of 13%, statistically significant above the 95% confidence level, as indicated by the error bars. In b, a total of 13 out of the 15 selected models (87%) simulate an increase in S-index variability in the future climate, with a multi-model mean increase of 16%, statistically significant above the 95% confidence level. Models that simulate an opposite change to the multi-model mean are indicated by green circles. A total of 33 CMIP5 models are used here as outputs from CMCC-CESM and FGOALS-g2 are unavailable.

Extended Data Fig. 9 Relationship between projected changes in ENSO and in pIOD.

a, Scatter plot of changes in September, October, and November (SON) S-index variability versus changes in SON Ni\({\tilde{\mathrm n}}\)o3.4 index variability. The change is defined as the difference between the present-day (1900–1999) and future (2000–2099) period scaled by the corresponding increase in global mean SST in each model. b, Same as a, but between changes in SON M-index variability and in SON Ni\({\tilde{\mathrm n}}\)o3.4 index variability. c, d, Same as a, b, but for changes in strong pIOD (S-index>1.5 s.d.) or moderate pIOD (M-index>1.25 s.d.) frequency (events per 100 years) vs frequency of El Ni\({\tilde{\mathrm n}}\)o (Ni\({\tilde{\mathrm n}}\)o3.4>1.0 s.d.). The associated correlation coefficient and P-value are plotted. There is no relationship between changes in El Niño and S-index in a and c. Although there is a positive correlation in b and d, the changes in M-index and in El Niño are in the opposite direction, contravening the relationship between El Ni\({\tilde{\mathrm n}}\)o and M-index.

Extended Data Fig. 10 Projected mean state change in the tropical Indian Ocean.

Shown are averages over the 20 models that simulate a nonlinear coefficient α greater than 50% of the observed value. a, Multi-model ensemble averaged mean state changes of SON SST (in °C) between the future (2000–2099) and the present-day (1900–1999) climate. b, Same as a, but for the mean state change in surface wind stress (in N m−2); changes in wind stress magnitude (in N m−2) are indicated by colors. c, Same as a, but for the mean state change in latent heat flux (in W m−2). d, Same as a, but for the mean state change in total heat flux (in W m−2); positive values indicate an upward transfer of energy. Stippled areas indicate where changes are statistically significant above the 95% confidence level according to a two-tailed Student’s t-tests.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cai, W., Yang, K., Wu, L. et al. Opposite response of strong and moderate positive Indian Ocean Dipole to global warming. Nat. Clim. Chang. 11, 27–32 (2021). https://doi.org/10.1038/s41558-020-00943-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41558-020-00943-1

This article is cited by

Search

Quick links

Nature Briefing Anthropocene

Sign up for the Nature Briefing: Anthropocene newsletter — what matters in anthropocene research, free to your inbox weekly.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing: Anthropocene