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.

Summertime increases in upper-ocean stratification and mixed-layer depth

Subjects

An Author Correction to this article was published on 02 August 2021

Abstract

The surface mixed layer of the world ocean regulates global climate by controlling heat and carbon exchange between the atmosphere and the oceanic interior1,2,3. The mixed layer also shapes marine ecosystems by hosting most of the ocean’s primary production4 and providing the conduit for oxygenation of deep oceanic layers. Despite these important climatic and life-supporting roles, possible changes in the mixed layer during an era of global climate change remain uncertain. Here we use oceanographic observations to show that from 1970 to 2018 the density contrast across the base of the mixed layer increased and that the mixed layer itself became deeper. Using a physically based definition of upper-ocean stability that follows different dynamical regimes across the global ocean, we find that the summertime density contrast increased by 8.9 ± 2.7 per cent per decade (10−5–10−4 per second squared per decade, depending on region), more than six times greater than previous estimates. Whereas prior work has suggested that a thinner mixed layer should accompany a more stratified upper ocean5,6,7, we find instead that the summertime mixed layer deepened by 2.9 ± 0.5 per cent per decade, or several metres per decade (typically 5–10 metres per decade, depending on region). A detailed mechanistic interpretation is challenging, but the concurrent stratification and deepening of the mixed layer are related to an increase in stability associated with surface warming and high-latitude surface freshening8,9, accompanied by a wind-driven intensification of upper-ocean turbulence10,11. Our findings are based on a complex dataset with incomplete coverage of a vast area. Although our results are robust within a wide range of sensitivity analyses, important uncertainties remain, such as those related to sparse coverage in the early years of the 1970–2018 period. Nonetheless, our work calls for reconsideration of the drivers of ongoing shifts in marine primary production, and reveals stark changes in the world’s upper ocean over the past five decades.

This is a preview of subscription content

Access options

\$32.00

All prices are NET prices.

Data availability

All information about the source database used in the paper is available at https://github.com/jbsallee-ocean/GlobalMLDchange/tree/main/Databases. The resulting global maps of trends and climatological fields presented here are available at https://zenodo.org/record/4073174#.YA_jsC2S3XQ (https://doi.org/10.5281/zenodo.4073174.

Code availability

The code used to generate the analysis presented in the paper and its Supplementary Information is available at https://github.com/jbsallee-ocean/GlobalMLDchange.

References

1. 1.

Bindoff, N. L. et al. in IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (eds Pörtner, H.-O. et al.) Ch. 5 (IPCC, 2019).

2. 2.

Frankignoul, C. & Hasselmann, K. Stochastic climate models, part II application to sea-surface temperature anomalies and thermocline variability. Tellus 29, 289–305 (1977).

3. 3.

Bopp, L., Lévy, M., Resplandy, L. & Sallée, J. Pathways of anthropogenic carbon subduction in the global ocean. Geophys. Res. Lett. 42, (2015).

4. 4.

Sverdrup, H. U. On conditions for the vernal blooming of phytoplankton. ICES J. Mar. Sci. 18, 287–295 (1953).

5. 5.

Behrenfeld, M. J. et al. Climate-driven trends in contemporary ocean productivity. Nature 444, 752–755 (2006).

6. 6.

Boyce, D. G., Lewis, M. R. & Worm, B. Global phytoplankton decline over the past century. Nature 466, 591–596 (2010).

7. 7.

Polovina, J. J., Howell, E. A. & Abecassis, M. Ocean’s least productive waters are expanding. Geophys. Res. Lett. 35, L03618 (2008).

8. 8.

Li, G. et al. Increasing ocean stratification over the past half-century. Nat. Clim. Chang. 10, 1116–1123 (2020).

9. 9.

Yamaguchi, R. & Suga, T. Trend and variability in global upper-ocean stratification since the 1960s. J. Geophys. Res. Oceans 124, 8933–8948 (2019).

10. 10.

Young, I. R. & Ribal, A. Multiplatform evaluation of global trends in wind speed and wave height. Science 364, 548–552 (2019).

11. 11.

Young, I. R., Zieger, S. & Babanin, A. V. Global trends in wind speed and wave height. Science 332, 451–455 (2011).

12. 12.

Sprintall, J. & Tomczak, M. Evidence of the barrier layer in the surface layer of the tropics. J. Geophys. Res. Oceans 97, 7305–7316 (1992).

13. 13.

Kataoka, T., Kimoto, M., Watanabe, M. & Tatebe, H. Wind–mixed layer–SST feedbacks in a tropical air–sea coupled system: application to the Atlantic. J. Clim. 32, 3865–3881 (2019).

14. 14.

Sallée, J. B., Matear, R., Rintoul, S. R. & Lenton, A. Localised subduction of anthropogenic carbon dioxide in the Southern Hemisphere oceans. Nat. Geosci. 5, 579–584 (2012).

15. 15.

Llort, J., Lévy, M., Sallée, J. B. & Tagliabue, A. Nonmonotonic response of primary production and export to changes in mixed-layer depth in the Southern Ocean. Geophys. Res. Lett. 46, 3368–3377 (2019).

16. 16.

Pörtner, H.-O. et al. Summary for policymakers. IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (eds Pörtner, H.-O. et al.) (IPCC, 2019).

17. 17.

Somavilla, R., González-Pola, C. & Fernández-Diaz, J. The warmer the ocean surface, the shallower the mixed layer. How much of this is true? J. Geophys. Res. Oceans 122, 7698–7716 (2017).

18. 18.

Capotondi, A., Alexander, M. A., Bond, N. A., Curchitser, E. N. & Scott, J. D. Enhanced upper ocean stratification with climate change in the CMIP3 models. J. Geophys. Res. Oceans 117, C04031 (2012).

19. 19.

Helm, K. P., Bindoff, N. L. & Church, J. A. Observed decreases in oxygen content of the global ocean. Geophys. Res. Lett. 38, L23602 (2011).

20. 20.

Kwiatkowski, L. et al. Twenty-first century ocean warming, acidification, deoxygenation, and upper-ocean nutrient and primary production decline from CMIP6 model projections. Biogeosciences 17, 3439–3470 (2020).

21. 21.

Pellichero, V., Sallée, J.-B., Schmidtko, S., Roquet, F. & Charrassin, J.-B. The ocean mixed layer under Southern Ocean sea-ice: seasonal cycle and forcing. J. Geophys. Res. Oceans 122, 1608–1633 (2017).

22. 22.

Treasure, A. M. et al. Marine mammals exploring the oceans pole to pole: a review of the MEOP consortium. Oceanography 30, 132–138 (2017).

23. 23.

Schmidtko, S., Johnson, G. C. & Lyman, J. M. MIMOC: a global monthly isopycnal upper-ocean climatology with mixed layers. J. Geophys. Res. Oceans 118, 1658–1672 (2013).

24. 24.

Schmidtko, S., Heywood, K. J., Thompson, A. F. & Aoki, S. Multidecadal warming of Antarctic waters. Science 346, 1227–1231 (2014).

25. 25.

Schmidtko, S., Stramma, L. & Visbeck, M. Decline in global oceanic oxygen content during the past five decades. Nature 542, 335–339 (2017).

26. 26.

Kuusela, M. & Stein, M. L. Locally stationary spatio-temporal interpolation of argo profiling float data. Proc. R. Soc. Lond. A 474, 20180400 (2018).

27. 27.

Buckingham, C. E. et al. The contribution of surface and submesoscale processes to turbulence in the open ocean surface boundary layer. J. Adv. Model. Earth Syst. 11, 4066–4094 (2019).

28. 28.

de Boyer Montégut, C., Madec, G., Fischer, A. S., Lazar, A. & Iudicone, D. Mixed layer depth over the global ocean: an examination of profile data and a profile-based climatology. J. Geophys. Res. 109, C12003 (2004).

29. 29.

Sallée, J. B., Speer, K. G. & Rintoul, S. R. Zonally asymmetric response of the Southern Ocean mixed-layer depth to the Southern Annular Mode. Nat. Geosci. 3, 273–279 (2010).

30. 30.

Holte, J., Talley, L. D., Gilson, J. & Roemmich, D. An Argo mixed layer climatology and database. Geophys. Res. Lett. 44, 5618–5626 (2017).

31. 31.

Rhein, M. et al. in Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the IPCC (eds Stocker, T. F. et al.) 255–316 (Cambridge Univ. Press, 2013).

32. 32.

Durack, P. J. & Wijffels, S. E. Fifty-year trends in global ocean salinities and their relationship to broad-scale warming. J. Clim. 23, 4342–4362 (2010).

33. 33.

Durack, P. J., Wijffels, S. E. & Matear, R. J. Ocean salinities reveal strong global water cycle intensification during 1950 to 2000. Science 336, 455–458 (2012).

34. 34.

Bindoff, N. L. et al. in Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the IPCC (eds Stocker, T. F. et al.) 867–952 (Cambridge Univ. Press, 2013).

35. 35.

Armour, K. C., Marshall, J., Scott, J. R., Donohoe, A. & Newsom, E. R. Southern Ocean warming delayed by circumpolar upwelling and equatorward transport. Nat. Geosci. 9, 549–554 (2016).

36. 36.

Marshall, J. et al. The ocean’s role in the transient response of climate to abrupt greenhouse gas forcing. Clim. Dyn. 44, 2287–2299 (2015).

37. 37.

Belcher, S. E. et al. A global perspective on Langmuir turbulence in the ocean surface boundary layer. Geophys. Res. Lett. 39, L18605 (2012).

38. 38.

Barkan, R., Winters, K. B. & McWilliams, J. C. Stimulated imbalance and the enhancement of eddy kinetic energy dissipation by internal waves. J. Phys. Oceanogr. 47, 181–198 (2017).

39. 39.

Hu, S. et al. Deep-reaching acceleration of global mean ocean circulation over the past two decades. Sci. Adv. 6, eaax7727 (2020).

40. 40.

Siegelman, L. et al. Enhanced upward heat transport at deep submesoscale ocean fronts. Nat. Geosci. 13, 50–55 (2020).

41. 41.

Su, Z., Wang, J., Klein, P., Thompson, A. F. & Menemenlis, D. Ocean submesoscales as a key component of the global heat budget. Nat. Commun. 9, 775 (2018).

42. 42.

Gill, A. E. & Niiler, P. P. The theory of the seasonal variability in the ocean. Deep-Sea Res. Oceanogr. Abstr. 20, 141–178 (1973).

43. 43.

Oschlies, A., Brandt, P., Stramma, L. & Schmidtko, S. Drivers and mechanisms of ocean deoxygenation. Nat. Geosci. 11, 467–473 (2018).

44. 44.

Killworth, P. D., Chelton, D. B. & de Szoeke, R. A. The speed of observed and theoretical long extratropical planetary waves. J. Phys. Oceanogr. 27, 1946–1966 (1997).

45. 45.

Durack, P. J., Wijffels, S. E. & Boyer, T. P. in Ocean Circulation and Climate Vol. 103 (eds. Siedler, G. et al.) Ch. 28, 727–757 (Academic Press, 2013).

46. 46.

Mantyla, A. W. The treatment of inconsistencies in Atlantic deep water salinity data. Deep Sea Res. Part I 41, 1387–1405 (1994).

47. 47.

Sloyan, B. M. et al. The global ocean ship-based hydrographic investigations program (GO-SHIP): a platform for integrated multidisciplinary ocean science. Front. Mar. Sci. 6, 445 (2019).

48. 48.

Riser, S. C. et al. Fifteen years of ocean observations with the global Argo array. Nat. Clim. Chang. 6, 145–153 (2016).

49. 49.

Roemmich, D. et al. On the future of argo: a global, full-depth, multi-disciplinary array. Front. Mar. Sci. 6, 439 (2019).

50. 50.

Roquet, F. et al. A Southern Indian Ocean database of hydrographic profiles obtained with instrumented elephant seals. Sci. Data 1, 140028 (2014).

51. 51.

Durack, P. J., Gleckler, P. J., Landerer, F. W. & Taylor, K. E. Quantifying underestimates of long-term upper-ocean warming. Nat. Clim. Chang. 4, 999–1005 (2014).

52. 52.

Brainerd, K. E. & Gregg, M. C. Surface mixed and mixing layer depths. Deep Sea Res. Part I 42, 1521–1543 (1995).

53. 53.

Thomson, R. E. & Fine, I. V. Estimating mixed layer depth from oceanic profile data. J. Atmos. Ocean. Technol. 20, 319–329 (2003).

54. 54.

Holte, J. & Talley, L. A new algorithm for finding mixed layer depths with applications to Argo data and subantarctic mode water formation. J. Atmos. Ocean. Technol. 26, 1920–1939 (2009).

55. 55.

Sallée, J. B., Wienders, N., Speer, K. & Morrow, R. Formation of subantarctic mode water in the southeastern Indian Ocean. Ocean Dyn. 56, 525–542 (2006).

56. 56.

Pedlosky, J. Geophysical Fluid Dynamics (Springer, 1987).

57. 57.

Jakobsen, P. K., Ribergaard, M. H., Quadfasel, D., Schmith, T. & Hughes, C. W. Near-surface circulation in the northern North Atlantic as inferred from Lagrangian drifters: variability from the mesoscale to interannual. J. Geophys. Res. Oceans 108, 3251 (2003).

58. 58.

Lacasce, J. H. Floats and f/H. J. Mar. Res. 58, 61–95 (2000).

59. 59.

Dijkstra, E. W. et al. A note on two problems in connexion with graphs. Numer. Math. 1, 269–271 (1959).

60. 60.

Roquet, F. et al. Delayed-mode calibration of hydrographic data obtained from animal-borne satellite relay data loggers. J. Atmos. Ocean. Technol. 28, 787–801 (2011).

61. 61.

Boehme, L. et al. Animal-borne CTD-satellite relay data loggers for real-time oceanographic data collection. Ocean Sci. 5, 685–695 (2009).

62. 62.

Kennedy, J., Rayner, N., Atkinson, C. & Killick, R. An ensemble data set of sea surface temperature change from 1850: the Met Office Hadley Centre HadSST. 4.0. 0.0 data set. J. Geophys. Res. Atmos. 124, 7719–7763 (2019).

63. 63.

Banzon, V., Smith, T. M., Chin, T. M., Liu, C. & Hankins, W. A long-term record of blended satellite and in situ sea-surface temperature for climate monitoring, modeling and environmental studies. Earth Syst. Sci. Data 8, 165–176 (2016).

64. 64.

Sutherland, P. & Melville, W. K. Field measurements of surface and near-surface turbulence in the presence of breaking waves. J. Phys. Oceanogr. 45, 943–965 (2015).

65. 65.

Sullivan, P. P., McWilliams, J. C. & Melville, W. K. Surface gravity wave effects in the oceanic boundary layer: large-eddy simulation with vortex force and stochastic breakers. J. Fluid Mech. 593, 405–452 (2007).

66. 66.

Li, M., Zahariev, K. & Garrett, C. Role of Langmuir circulation in the deepening of the ocean surface mixed layer. Science 270, 1955–1957 (1995).

67. 67.

Li, M. & Garrett, C. Mixed layer deepening due to Langmuir circulation. J. Phys. Oceanogr. 27, 121–132 (1997).

68. 68.

Edson, J. B. et al. On the exchange of momentum over the open ocean. J. Phys. Oceanogr. 43, 1589–1610 (2013).

69. 69.

Tokinaga, H. & Xie, S.-P. Wave- and anemometer-based sea surface wind (waswind) for climate change analysis. J. Clim. 24, 267–285 (2011).

70. 70.

Ozmidov, R. On the turbulent exchange in a stably stratified ocean. Atmos. Oceanic Phys. 1, 493–497 (1965).

71. 71.

Riley, J. J. & Lelong, M.-P. Fluid motions in the presence of strong stable stratification. Annu. Rev. Fluid Mech. 32, 613–657 (2000).

72. 72.

Klein, P. & Lapeyre, G. The oceanic vertical pump induced by mesoscale and submesoscale turbulence. Annu. Rev. Mar. Sci. 1, 351–375 (2009).

73. 73.

Lapeyre, G., Klein, P. & Hua, B. L. Oceanic restratification forced by surface frontogenesis. J. Phys. Oceanogr. 36, 1577–1590 (2006).

74. 74.

Thomas, L. N. Destruction of potential vorticity by winds. J. Phys. Oceanogr. 35, 2457–2466 (2005).

75. 75.

Thomas, L. N. & Lee, C. M. Intensification of ocean fronts by down-front winds. J. Phys. Oceanogr. 35, 1086–1102 (2005).

76. 76.

D’Asaro, E., Lee, C., Rainville, L., Harcourt, R. & Thomas, L. Enhanced turbulence and energy dissipation at ocean fronts. Science 332, 318–322 (2011).

77. 77.

Tandon, A. & Garrett, C. Mixed layer restratification due to a horizontal density gradient. J. Phys. Oceanogr. 24, 1419–1424 (1994).

Acknowledgements

This project received funding from the European Union’s Horizon 2020 Research and Innovation programme under grant agreement number 821001. V.P., C.A., E.P. and L.V. received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation programme (grant agreement 637770). A.N.G. acknowledges the support of the Royal Society and the Wolfson Foundation. P.S. received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation programme (grant agreement No 805186). We thank I. Young for providing the percentage rate change of 10-m wind speed from recently published analysis. We thank L. Sigelman and G. Madec for comments and discussions that greatly helped us to refine our study.

Author information

Authors

Contributions

J.-B.S. designed the experiment and performed the computations and data analyses; V.P., C.A., E.P. and L.V. helped with the development of the global database and its analysis, and evaluated the analysis; S.S. developed the mapping method; and A.N.G. and P.S. provided expertise on surface ocean turbulence and associated scaling arguments. M.K. provided expertise on the statistical methods used in this study. All authors discussed the results and wrote the manuscript.

Corresponding author

Correspondence to Jean-Baptiste Sallée.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review information Nature thanks Paul Durack, Stephen Riser 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 figures and tables

Extended Data Fig. 1 Trends associated with the vertical structure of the upper ocean.

a, b, Schematics showing idealized density profiles in the upper ocean for the cases in which the mixed layer and pycnocline are shallower (a) and deeper (b) than 200 m. The black line shows the typical shape of the density profile with a total density increase of ∂ρ across the pycnocline (thickness h) and the mixed layer (thickness H). The dashed red lines show the density profiles after lightening of the mixed layer with no change of mixed-layer depth, and the dotted red lines show the density profiles after lightening of the mixed layer concomitant with a deepening of the mixed layer.

Extended Data Fig. 2 Geographical distribution of available observations.

ad, Number of mixed layer estimates in 1° × 1° longitude–latitude bins: from all available observation sources (a), from ship profiles (b), from Argo profiles (c) and from instrumented marine mammal observations (d). e, f, Maximum time span (in years) covered by the combined dataset in 1° × 1° longitude–latitude bins in summer (e) and in winter (f).

Extended Data Fig. 3 Impact of linear-regression choices on mean mixed layer.

ad, Winter (a, b) and summer (c, d) mean mixed-layer depth computed using slightly different linear-regression models: Choice 1 (a, c; covariance between observations) and Choice 2 (b, d; no covariance between observations). See Methods for more details.

Extended Data Fig. 4 Impact of linear-regression choices on summer mixed-layer depth, stratification trends and their associated standard errors.

ah, 1970–2018 summer trend for mixed-layer depth (a, c) and summer pycnocline stratification (e, g) and their associated standard error: standard error of mixed-layer depth trend (b, d) and standard error of summer pycnocline stratification trend (f, h). a, b, e, f show the solution computed with the linear-regression model Choice 1 (covariance between observations); c, d, g, h show the solution computed with the linear-regression model Choice 2 (no covariance between observations).

Extended Data Fig. 5 1970–2018 trends in winter upper-ocean stratification and mixed-layer depth.

a, b, Map of the 1970–2018 winter 0–200 m trend (a; $${N}_{200}^{2}$$ trend in s−2 dec−1) and pycnocline stratification trend (b; N2 trend in s−2 dec−1); thick black lines show the zonal-median value and thin black lines show the 33th–66th percentile. Regions with no significant trend (see Methods) are shaded in grey on the map. c, As in a, b, but for winter mixed-layer trend, in m dec−1 (mixed-layer deepening is shown as a negative trend).

Extended Data Fig. 6 Regional time series of winter pycnocline stratification and mixed-layer depth anomaly.

a, Winter climatological mixed-layer depth (same as Fig. 2f) with three specific regions of interest outlined by red contours: North Atlantic subpolar convection region (A); Southern Ocean Indian sector convection region (B); and Southern Ocean Pacific sector convection region (C). b, d, f, Winter stratification anomaly time series and associated trends for regions A (b), B (d) and C (f). c, e, g, Winter mixed-layer depth anomaly times series and associated trends for regions A (c), B (e) and C (g). A negative depth anomaly refers to a deepening. Each time series shows: thin grey line, the annual median percentage anomaly (from the local climatological seasonal cycle), computed for each individual observation; error bars referring to the 33th–66th percentile range of the percentage anomaly (error bars are shown in black (grey) when more (fewer) than 50 data points are used in the annual statistics); the associated five-year smoothed median time series superimposed in blue; and a linear trend from 1970–2018, shown by the red line, if greater than twice its standard error.

Extended Data Fig. 7 Comparison between mixed-layer temperature trend and SST trends.

a, Summer mixed-layer mean temperature trend from 1970 to 2018, as estimated in this study. b, Summer SST trend from 1982 to 2018, as estimated from the satellite-based product GHRSSTv2. c, Box plot showing the median (red) and interquartile range (blue box) of local summer SST trend estimates from this study (mixed-layer mean temperature from 1970–2018), from the satellite-based product GHRSSTv2 (SST from 1982–2018) and from the in situ observation reconstruction product HadSSTv4 (SST from 1970–2018). The whiskers extend to the most extreme data points.

Extended Data Fig. 8 Difference between Argo- and ship-based derived mixed-layer depth.

a, b, Difference between mixed-layer depth (MLD) estimates coming from nearby (sampled within 330 km and 1.5 day) Argo and ship-based observation profiles (that is, co-located in time and space) for all instances for which we derived a smaller (a) or a greater (b) mixed-layer depth from the Argo profile than from the ship-based profile. c, Histogram of all differences. Because Argo started in the 2000s, the co-located profiles cover only the years 2000–2018.

Extended Data Fig. 9 1970–2018 trends in summer pycnocline stratification and mixed-layer depth when using only ship-based profiles (removing all Argo and MEOP programme observations).

a, Map of the 1970–2018 summer pycnocline stratification trend (N2 trend, in s−2 dec−1) along with the zonal-median value: median (thick black line) and 33th–66th percentile (thin black line). The red shading shows the global 33th–66th percentile range of the local trend estimates. Regions with no significant trend (see Methods) are shaded in grey on the map. b, As in a but for the summer mixed-layer trend, in m dec−1 (mixed-layer deepening is shown as a negative trend).

Supplementary information

Supplementary Information

This file contains Supplementary Text, Supplementary Figures 1 – 3 and Supplementary Table 1.

Rights and permissions

Reprints and Permissions

Sallée, JB., Pellichero, V., Akhoudas, C. et al. Summertime increases in upper-ocean stratification and mixed-layer depth. Nature 591, 592–598 (2021). https://doi.org/10.1038/s41586-021-03303-x

• Accepted:

• Published:

• Issue Date:

• DOI: https://doi.org/10.1038/s41586-021-03303-x

• Recent upper Arctic Ocean warming expedited by summertime atmospheric processes

• Zhe Li
• Qinghua Ding
• Axel Schweiger

Nature Communications (2022)

• Impact of intensifying nitrogen limitation on ocean net primary production is fingerprinted by nitrogen isotopes

• Pearse J. Buchanan
• Olivier Aumont
• Alessandro Tagliabue

Nature Communications (2021)

• Biogeochemical extremes and compound events in the ocean

• Nicolas Gruber
• Philip W. Boyd
• Meike Vogt

Nature (2021)