Introduction

In recent decades, episodes of warm ocean temperature extremes have been associated with more intense and frequent impacts on marine organisms, ecosystems and reliant human industries around the world1,2,3. By analogy with their atmospheric counterpart, these extreme ocean temperature events have been termed “marine heatwaves” (MHWs)4,5. Some of the most prominent events, together with the unprecedented warming during the boreal summer of 2023 are presented in Box 1. MHWs influence regional climate phenomena and often drive substantial impacts on the marine environment. For example, MHWs in the Indian Ocean have been found to modulate the monsoon winds and rains over the Indian subcontinent, impacting water and food security over the region6. MHWs interact with and intensify tropical cyclones, making them more destructive7,8,9,10. Biological MHW impacts include mass mortality events in invertebrates, fish, birds and marine mammals1,11,12,13, coral bleaching14,15, declines in foundation species3,16,17 and entire ecosystem restructuring18,19, with far-reaching socioeconomic impacts20.

Recent reviews and perspectives13,21,22 have outlined major steps forward in understanding MHW characteristics, drivers, and predictability, along with the economic impacts they cause. However, in this rapidly evolving field, more recent research has provided new insights into MHWs, while generating important new questions and research avenues. Although MHW research has primarily considered temperature extremes at the ocean surface, subsurface temperature extremes may be more intense and longer-lasting than their surface counterparts23,24,25,26,27. Given the prevalence of life throughout the water column, subsurface MHWs need to be closely observed, mechanistically understood, and skillfully predicted. In addition, while the physical characterization of MHWs has mainly focused on large-scale events (Box 1), MHWs are now also studied in more localized coastal areas, marginal seas, and fjords28,29,30,31, where they are negatively impacting the local ecology and coastal communities3,16. MHWs are also increasingly being examined along with other extreme conditions, like high acidity or low-oxygen32,33, sea level extremes34, floods35, droughts36, severe weather events37 or even terrestrial heat waves over the adjacent land38. These “compound events” act as multiple stressors for marine life and societies.

The ability to predict MHWs and compound events from days to seasons in advance is key for stakeholder preparation and mitigation efforts39. Skillful forecasts require enhanced understanding of MHW drivers to assess their predictability, and prediction systems that realistically capture the processes underpinning that predictability21,40,41. While progress has been made in prediction activities42,43, additional improvements could be achieved through a deepened understanding of the relative roles of different MHW drivers, and dynamical model improvements, which include an assessment of the sensitivity of MHW forecasts to model resolution44,45,46,47.

As the oceans continue to warm with anthropogenic climate change48,49, defining MHWs under non-stationary conditions becomes increasingly challenging, as commonly used definitions will lead to a permanent MHW state in areas experiencing sufficient warming50 (Fig. 1). In addition, separating the processes internal to the climate system from those of anthropogenic origin51,52 is key to the mechanistic understanding of the nature of MHWs and the assessment of their predictability and their future changes.

Fig. 1: Influence of trends on marine heatwave definition.
figure 1

a SST anomaly (SSTA, °C) relative to the 1982–2011 climatology over the Mediterranean during July 2022 (dashed line in b) from NOAA-OISSTv2.176. Anomalies include the trend signal. b Monthly SST anomalies (seasonal cycle removed) averaged over the western Mediterranean (the region bounded by the black lines and the coast, i.e., north of 40°N, and west of 12°E). c As in (a), but with SST anomalies linearly detrended. d as in (b), but with anomalies linearly detrended. In (b) and (d), the thick blue horizontal line indicates the baseline period used to compute the climatology and to define the 90th percentile (thin horizontal blue lines), which for simplicity is chosen to be seasonally independent. In the presence of a trend, the western Mediterranean tends to be in a quasi-permanent MHW state toward the end of the record, while removal of the trend highlights isolated events since 2015, and also reveals more pronounced extreme events at the beginning of the record. The MHW in the boreal summer of 2003 emerges as an extremely intense event, irrespective of the presence of a trend signal.

This article extends previous reviews13,21,22 by highlighting the new emerging areas in MHW research outlined above, including: a critical re-evaluation of MHW definitions and their detection, both at the surface and in the subsurface, in the presence of climate change; observational needs and new emerging “observing” strategies; advances in the understanding of both surface and subsurface MHW drivers to aid prediction efforts; compound events and their prediction; and investigations to assess future MHW projections using empirical approaches and state-of-the-art modeling systems. This review also provides a perspective on new and promising avenues for advancing our understanding and prediction capabilities of ocean extremes in the context of our changing climate.

Defining a marine heatwave

Defining a MHW involves multiple choices, each leading to outcomes with distinct implications. These choices may be motivated by the need to understand the physical drivers or impacts of a MHW, or they can be constrained by the characteristics of the available data, like record length or temporal resolution. For simplicity, MHWs have typically been analyzed using local definitions53. However, since MHWs have a three-dimensional structure that evolves over time, other approaches are emerging25,54,55 to facilitate the tracking of extended surface or subsurface events over time.

Over the past decade, the majority of studies have adopted a common framework for defining MHWs. Following the widely used Hobday et al.53 framework, a MHW occurs at a given location when daily sea surface temperature anomalies exceed the seasonally-varying 90th percentile climatology for five days or more (with dips below this threshold for two days or less ignored). The 90th percentile climatology is typically based on a fixed reference period, or “baseline”.

These threshold criteria were chosen in analogy with atmospheric heatwaves56, and were not necessarily dictated by specific impacts in the marine environment. As such, other definitions have also been employed22, including, for example, definitions using the 99th percentile57, approaches using monthly data41,42,51,52,58,59 instead of daily data (Fig. 1), annual maximum temperatures60, or cumulative temperatures exceeding fixed thresholds, a criterion commonly used for coral bleaching monitoring and prediction61,62,63. Attempts to incorporate information on biological impacts has led to the creation of MHW hazard indices, where species-tailored metrics were co-developed with stakeholders using absolute temperatures64. With a fixed baseline, MHW conditions will become increasingly common as the ocean warms57,65, potentially leading to a “permanent” MHW state in regions experiencing a high level of warming10 (Fig. 1). These changing characteristics may reflect the risk these events pose to some marine organisms, particularly those with slow adaptation rates13. However, considering a fixed baseline limits our ability to distinguish the slow climate change-related processes from the faster processes associated with internal modes of climate variability or synoptic weather conditions40, with implications for understanding events’ predictability and assessing their prediction skill66. Thus, there has been a recent call to remove the effects of mean warming when defining MHWs by detrending temperature time series67 or using a shifting baseline period68, especially for future projections. To define MHW characteristics, the decision to use a temperature threshold that remains fixed or changes over time will ultimately depend on the application being studied, the importance of maintaining consistency with past studies and the characteristics of the data record.

Given the availability of satellite-derived sea surface temperature (SST), many MHW studies have relied on daily gridded satellite data, starting in the early 1980s. However, different datasets will have varying temporal and spatial resolutions, different interpolation accuracies, and may be sporadic and contain data gaps. Modified MHW definitions may be appropriate for different datasets and specific applications. For instance, monthly means can be used in regions characterized by long ocean memory (e.g., the tropical Pacific), or when the focus is on long-lasting MHWs20,69. In all definitions, the temporal and spatial scales of the dynamics at play need to be considered.

Spatially, MHWs can cover large horizontal areas and extend deep into the water column. MHW structure is linked to their drivers and needs to be included in their characterization. While horizontal extent has been considered in studies assessing MHW projections57,70, a number of new techniques have been developed to track connected MHW regions at the surface54,55,71 or in three-dimensional space25,72, accounting for the splitting and merging of MHW regions. These techniques treat MHWs as objects that evolve in space and time providing an illustration of their areas of influence. Similar algorithms have also been applied to ocean acidification extremes in the Northeast Pacific73, allowing an assessment of the severity of their impacts. To date, these tracking algorithms are purely statistical and do not incorporate information about event dynamics, but the use of tracking will provide new opportunities for understanding the extent of MHW systems as they evolve through time, and facilitate the identification of their underlying dynamics.

Observations for characterizing marine heatwaves

Observations are the foundation for characterizing and understanding MHWs. For more than a century, a diversity of ways to measure ocean temperature have been developed from in situ stationary and moving platforms (both passive and active) to remotely sensed methods (Box 2; Table 1 of Oliver et al.22).

The challenge for observing MHWs is to measure ocean temperature at high temporal resolution and over a long period (decades) to define a threshold for extremes, while accounting for the inherent variability of temperature at different timescales. Advances in understanding MHWs globally have relied largely on satellite derived sea surface temperature products blended with near-surface in situ data provided from surface drifting buoys and ship underway systems (e.g., products such as the Operational SST and Ice Analysis (OSTIA) system74,75 and NOAA Daily Optimum Interpolation SST v2.1 dataset76). On the other hand, long-term in situ temperature measurements from water samples and moorings have been crucial for characterizing temperature extremes at the daily timescale, although not representative of large areas. There are coastal locations distributed worldwide where ocean temperatures have been recorded since before the satellite era65,77,78,79, providing insight into long-term trends of local ocean temperatures and changes in the frequency of temperature extremes. Only a few sites include sustained measurements extending through the water column, which have been crucial to understand subsurface MHW characteristics and drivers23,24,80. Globally, the network of Argo floats has provided a transformative capability to study subsurface events by sampling ocean temperatures27,81 over the upper 2000 m for more than two decades82. However, using raw Argo profiles poses major analytical challenges, and Argo-derived gridded datasets are primarily available at monthly time resolution and lack coverage over continental shelves and marginal seas. Nevertheless, the combination of different observational platforms, for example, Argo and coastal moorings29, is proving extremely valuable to achieve a comprehensive view of MHWs (Box 2).

To overcome issues associated with sparse and inconsistent observations, and extend analyses of MHWs back to the pre-satellite era, many studies have leveraged ocean reanalyses—models constrained by observations through data assimilation—to understand MHW drivers and dynamical processes41,83,84, and analyze MHW characteristics at both the ocean surface and in the subsurface25,26. Ocean reanalyses offer uniform data coverage in time and space, in some cases at high-resolution (e.g., GLORYS12v185), thus also facilitating the characterization of MHWs on continental shelves86,87. Reanalysis products are subject to model errors and biases and may differ in their representation of MHWs over past decades88. Thus, they should be used with care to study extremes, especially in areas where limited observations were assimilated (e.g., the deep ocean or shelf regions). Some variables, like biogeochemical properties, are much less constrained by observations than physical quantities like temperature.

At longer timescales, a useful approach for increasing the sample size of extreme events is to use empirical models trained on observations, like Linear Inverse Models (LIMs89), to produce multi-millennia synthetic time series. These synthetic data share similar statistical properties (covariances, autocorrelation, event evolution) with observations35,51,58 and allow exploration of the full range of possible MHW realizations that are consistent with the dynamics and noise structure of the training data.

Drivers of surface and subsurface marine heatwaves

The processes driving MHWs affect their characteristics, including duration, intensity and vertical structure, and are key for predicting their evolution. MHWs are driven by local heat fluxes associated with synoptic atmospheric conditions or ocean advective and mixing processes, and are sensitive to the ocean state (e.g., mixed layer depth). These local drivers may themselves be modulated by large-scale modes of climate variability or anthropogenic warming, and vary regionally and seasonally. In the extra-tropics, intense MHWs are commonly associated with persistent atmospheric highs, resulting in increased insolation and decreased wind speeds that reduce turbulent heat losses and vertical ocean mixing40,90,91,92,93,94,95. Associated shallower mixed layers can further amplify the warming from surface heat fluxes59,83,96,97. More broadly, heat budget analyses indicate that increased insolation and reduced evaporative cooling typically dominate the build-up of MHWs while decay is commonly driven by increased turbulent heat losses98. In boundary current regions, anomalous warm oceanic advection is often important. Key examples include the 2011 Ningaloo Niño5,99,100 and the long-lived 2015/16 Tasman Sea MHW101,102. Advection-driven MHWs typically have a smaller surface area, but last longer21 and may reach deeper103 than atmospherically-driven events21. In the tropical Pacific, MHWs associated with El Niño Southern Oscillation (ENSO) are dynamically driven104, with surface heat fluxes damping temperatures during both the onset and decay phases98,105. In high-energy regions, like western boundary currents, oceanic mesoscale eddies and meanders, which often cross onto the shelf, can contribute to the onset, intensity, and longevity of MHWs at small temporal and spatial scales8,47,81,106. Changes in atmosphere-ocean interactions, including positive cloud feedbacks (reduction of low-cloud cover leading to enhanced insolation) in response to the initial SST anomalies, may contribute to the maintenance and intensification of those anomalies, and increase their persistence, as documented for long-lasting events in the northeast Pacific107,108.

Large-scale modes of climate variability can affect the likelihood of MHW occurrences regionally6,40, by modulating the local drivers and the initial upper-ocean stratification. For example, the 2013/14 MHW in the southwest Atlantic was forced by atmospheric conditions associated with a wave train triggered by the Madden-Julian Oscillation in the tropical Indian Ocean36. ENSO events are associated with a significant increase in the frequency, intensity and duration of MHWs across many parts of the global oceans65,90, enhancing the forecast skill of MHWs and ocean acidity extremes42,109, and influencing MHW projections52. For example, stronger equatorial Pacific easterly winds during La Niña events lead to an enhancement of the Indonesian Throughflow and Leeuwin Current, thereby transporting warm tropical waters to the Western Australian coast, creating favorable conditions for the development of MHWs99. ENSO also affects the Northeast Pacific Ocean through both oceanic and atmospheric pathways110. However, ENSO’s influence on MHWs may depend on the location of ENSO-related SST anomalies84, and may be mediated by other modes of variability at interannual or decadal timescales40,41,94,111,112. For example, a pre-existing positive Indian Ocean Dipole can increase the likelihood and predictability of MHWs off Western Australia up to 20 months in advance112, while in the Northeast Pacific, MHW onset is influenced by low-frequency variability related to the Pacific Decadal Oscillation (PDO)41,111. Interactions between tropical basins113,114 can also contribute to MHW development, as exemplified by the 2020 MHWs in the Northwest Pacific and South China Sea115,116 and the unprecedented Northwest Pacific event of 2022117. Finally, MHWs in the far-eastern tropical Pacific (“coastal El Niño events”), result primarily from the constructive interference of the North and South Pacific Meridional Modes118,119, and are not necessarily related to basin-wide ENSO conditions35. Assessing the relative contributions and links between large-scale drivers is critical to fully understand and exploit the inherent system predictability and improve predictions120.

MHWs extend into the subsurface ocean, with depth structures that may vary considerably depending on the region27, the leading driving mechanisms and the local bathymetry, whether in open ocean or on the shelf. MHWs can be confined to the mixed layer (“shallow MHWs”), driven by enhanced air-sea heat fluxes, ocean advection or reduced wind-induced turbulent mixing (Fig. 2a, e), or they can penetrate well below the mixed layer27. In shallow coastal regions, they can even extend to the ocean bottom (“Extended” events; Fig. 2c, d)24,86 due to the intrusion of warm eddies and western boundary current meanders onto the shelf (Fig. 2c) or through warm advection by alongshore currents (Fig. 2d), as shown by data from a near-shore mooring site in eastern Australia24. Near the shelf, deep warm anomalies can result from downwelling processes and may coexist with surface cooling24 (Fig. 2b). More generally, subsurface intensification through the dynamical movement of the thermocline may result from local Ekman pumping or from the passage of large-scale planetary waves (Fig. 2f), processes that occur ubiquitously throughout the ocean26. These subsurface anomalies are often larger than surface anomalies, due to the movement of the strong vertical temperature gradients around the thermocline23,24,26,27,81,86,121. As they evolve over time121 (Fig. 2e–h), subsurface MHWs can extend below the mixed layer during its seasonal shoaling and persist at depth even though the surface layer cools (Fig. 2g). They may sometimes be entrained back into a deepening mixed layer, and produce a delayed surface warming, a process known as “re-emergence”122 (Fig. 2h). Such evolutions were found in a simulation of the eastern tropical and North Pacific121 and in Argo observations of Northeast Pacific MHWs during 2004–2020123, and are likely to occur in other regions124.

Fig. 2: Vertical structures of MHWs.
figure 2

ad Possible vertical structures of MHWs near the shelf, including: “shallow” MHWs which do not penetrate below the mixed layer (a); “Bottom” intensified events due to a downwelling thermocline near the bottom, resulting, for example, from alongshore winds, as illustrated for the Southern Hemisphere (b); “Extended” profiles from the surface to the bottom due to intrusion of warm eddies or western boundary meanders into the shelf (c) or due to warm alongshore advection (d). eh Temporal evolution of subsurface MHWs associated with: changes in upper-ocean mixing for shallow events (e); propagation of oceanic Rossby waves causing variations in thermocline depth (f); persistence of deep anomalies with no surface signature due to mixed layer shoaling (g); and re-emergence of deep anomalies at the surface when the mixed layer deepens (h). The subsurface structure of MHWs depends on the processes involved in their formation, as well as the region’s stratification and circulation.

Compound and cascading events

Compound events are generally defined as a combination of extreme conditions and/or hazards that contribute to societal or environmental risk125. As such, they stress both natural and human systems, causing socioeconomic impacts such as loss of essential ecosystem services and income125. Understanding their underlying physical processes is thus critical for predictability assessments. During a compound event, extreme conditions can occur simultaneously (e.g., high ocean temperatures and low oxygen concentrations in the ocean) or in close sequence, where one event can increase the system vulnerability to a successive event. For instance, droughts and heatwaves can lead to a higher risk of flash floods over land. Events can also occur concurrently over different regions with large-scale consequences, as exemplified by the widespread impacts on fisheries caused by MHWs and low upper ocean nutrient levels during El Niño events.

Marine heatwaves and terrestrial extremes

Our understanding is more advanced for extremes and compound events over land126. However, work into ocean-land compound events is growing. For instance, atmospheric blocking over eastern South America and the western South Atlantic is associated with persistent high-pressure centers that can reduce cloud cover and latent heat loss, leading to simultaneous drought conditions over land and heatwaves in the adjacent ocean36. Similarly, synoptic conditions driving terrestrial heatwaves in some locations around Australia are found to be conducive to the warming of the ocean, increasing the likelihood of a concurrent MHW38. More generally, as many extreme extra-tropical MHWs are associated with persistent high-pressure centers90, such systems, straddling the land and ocean, might plausibly lead to compound marine-terrestrial temperature extremes in coastal regions. MHWs can also be related to enhanced evaporation and transport of humidity, inducing heavy rainfall along coastal regions, such as during the Tasman Sea MHW in 2015–16101 or the coastal MHW off Peru in 2017127.

Marine heatwaves and ocean biogeochemical extremes

Given their potential impacts on marine organisms, there is growing interest in ocean biogeochemical extremes that can co-occur with temperature extremes (Fig. 3), including high acidity (OAX)73,128,129,130,131, low oxygen (LOX)132 and low chlorophyll extremes (LChl)133,134. These stressors may act additively or synergistically135. For example, compound MHW-LOX events can have detrimental effects on aerobic metabolic rates, especially in ectotherms, i.e., cold-blooded organisms136,137,138. Additionally, compound MHW-OAX events can adversely affect molluscs139 or warm-water corals140, while MHW-LChl events are often associated with extremely low fish biomass conditions141. Moreover, it is plausible that the concurrent extreme ocean acidity conditions amplified the devastating effects of the Northeast Pacific Marine Heatwave of 2014–2015142 (Box 1, top panel), also known as the “Blob”83.

Fig. 3: Near-surface biogeochemical anomalies and compound conditions during some impactful MHWs.
figure 3

The biogeochemical quantities are shown for the month and the area of the MHWs displayed in the top panel of Box1’s figure. The footprint of those MHWs is indicated by gray lines. a Percentile associated with the mean chlorophyll anomaly during the MHWs, compared to the local empirical distribution of chlorophyll monthly anomalies from 1998 to 2018. b Percentile associated with the mean [H+] anomaly during the MHWs, compared to the local empirical distribution of [H+] monthly anomalies from 1982 to 2019, based on observationally-derived data143. c Extent of the MHWs co-occurring with a low chlorophyll extreme event (MHW-LChl, in blue), a high acidity event (MHW-OAX, in red), and both (MHW-LChl-OAX, in yellow). LChl events are defined as events with chlorophyll anomaly percentiles on panel (a) lower than their 10th percentile, and OAX events as events with [H + ] anomaly percentiles exceeding their 90th percentile. The chlorophyll data, corresponding to the mean chlorophyll concentration within the mixed layer, are obtained from the NASA Ocean Biogeochemical Model reconstruction176, and are publicly available for 1998–2021 (https://gmao.gsfc.nasa.gov/gmaoftp/rousseaux/Carlos/NOBM/).

Compound MHW-OAX events are more likely to occur in the subtropics than in the equatorial Pacific and mid-to-high latitudes, as high temperatures in the subtropics strongly increase the hydrogen ion [H+] concentration (i.e., acidity)143. At higher-latitudes, lower background temperatures limit this effect, while in the equatorial Pacific, reduced Dissolved Inorganic Carbon due to weaker upwelling leads to a decreased [H+] concentration, counteracting the effect of temperature. Conversely, hotspots of compound MHW-LChl events are found in the equatorial Pacific, along the boundaries of the subtropical gyres and in the northern Indian Ocean, often associated with El Niño events133 and enhanced nutrient limitation on phytoplankton growth134,144. Notably, the North Pacific MHW in 2014–2016 was identified as a quadruple compound event during some phases of its development, involving high temperature, low oxygen, high acidity and low chlorophyll levels32,133,145. For example, in January 2014, the extreme warming of the Blob (Box 1, top panel) co-occurred with low chlorophyll over part of the MHW area (Fig. 3c).

Climate model projections indicate that long-term trends in acidification, deoxygenation, and nutrient decline in the low-latitude upper ocean will persist for decades146,147, amplifying the frequency, intensity and scale of compound MHWs and biogeochemical extremes32,143. Notably, even when using a shifting baseline, whereby the effect of long-term warming and OAX are removed, OAX events and compound MHW-OAX events are expected to increase due to projected increases in the seasonal and diurnal variations in [H+]130,148,149.

Despite initial studies, understanding ocean compound extreme events is still in its infancy32. A global perspective on the temporal and spatial characteristics of these events, especially at depth, and a mechanistic understanding of relevant processes and their cascading impacts on ecosystems, are currently missing, mainly due to the lack of available subsurface data.

Climate models representation of marine heatwaves

Given the sparsity of long-term observational records, especially at depth, numerical models can help us to better understand MHW characteristics and their drivers. Global Earth system models (ESMs), which include both physical and biogeochemical components, are essential to provide future projections of both MHWs and biogeochemical extremes. But how well do climate models represent MHWs? Global coupled climate models vary in their degree of fidelity in representing climatological characteristics of basic MHW metrics (frequency, intensity and duration) at both daily44,60, and monthly51,52 timescales. Generally, CMIP-type ESMs tend to overestimate the duration of MHWs47,52,150 (Fig. 4). In addition, regions of strong ocean currents are especially problematic in models without eddy-permitting resolution44,45,47, due to the models’ inability to capture the influence of mesoscale eddies on MHW development in those regions106. Observed changes in MHW characteristics over the historical period are also challenging for models to simulate, although those stemming from mean state changes are better represented than those due to changes in internal variability51. However, the observational record of surface MHWs is relatively short, consisting of approximately 40 years for daily SSTs derived from satellite remote sensing, and approximately 100 years for monthly SSTs measured by ships of opportunity, with subsurface MHW records being even shorter. Such observational records only provide a limited sample of all the possible realizations that are consistent with the dynamics and noise of the climate system.

Fig. 4: Fidelity in climate models’ representation of MHW statistics.
figure 4

af Composite MHW intensity (°C) and (gl) composite MHW duration (months) during 1950–2020 from the 100-member of the Community Earth System Model version 2 (CESM2) SMILE and observations (ERSSTv5)175. a, g Ensemble average; (b, h) Observations; (c, i) Ensemble average minus Observations; (d, j) Ensemble maximum; (e, k) Ensemble minimum; (f, l) Ensemble maximum minus minimum. Gray shading in (c, i) indicates that observations lie within the 5th–95th percentile range of the CESM2 Large Ensemble. Adapted from Deser et al.52 © American Meteorological Society, used with permission.

Unlike observations, coupled climate models offer the potential for multiple realizations of the past and future, thereby enhancing sample sizes of extreme events. In particular, so-called “Single Model Initial-condition Large Ensembles” (SMILEs) have become a powerful tool in climate research for studying the simulated characteristics of internal variability and forced responses on local and regional scales151. SMILEs consist of many historical and future scenario simulations (generally 30–100) for a particular model, each starting from slightly different initial conditions, and allow a clean separation between the forced signal (the ensemble mean) and internal variability/extremes (departures from the ensemble mean). The power of SMILEs is only beginning to be exploited for the study of MHWs and their projected changes52,57,143,152.

Figure 4 illustrates the effect of sampling uncertainty on MHW characteristics during the historical period 1950–2020, obtained from the 100-member CESM2 SMILE52, after the forced signal is removed. The composite MHW intensity (Fig. 4a) and duration (Fig. 4g) metrics based on all ensemble members mask the considerable range found across individual realizations (Fig. 4d–f, j–l), underscoring the need for SMILEs to guard against sampling uncertainty. Since the single observational record (Fig. 4b, h) provides a limited sample size of extreme events, it may be challenging to separate true model biases from apparent biases stemming from inadequate sampling. One approach is to assess whether the characteristics of the single observed composite MHW lie outside the plausible (5th–95th percentile) range across SMILE members, in which case the model shows a likely bias (Fig. 4c, i). ESMs are also used for seasonal predictions of MHWs and other biogeochemical extremes. Thus, the fidelity of models in accurately simulating such extremes is critical for assessing the reliability of their predictions.

Prediction of marine heatwaves and associated biogeochemical extremes

Understanding MHW predictability and building effective prediction systems can greatly benefit marine management21. For example, accurate seasonal predictions of MHWs have the capacity to transform resource management practices that affect ecosystem services such as fisheries, aquaculture, and tourism21,43. Motivated by many potential benefits, recent research has quantified subseasonal-to-seasonal MHW predictability and forecast skill using dynamical and statistical approaches41,42,153,154,155.

Forecast systems based on global climate models have been used to estimate MHW probabilistic forecast skill and errors by comparing initialized hindcasts (i.e., retrospective forecasts) with the actual evolution of historical temperature anomalies. Results indicate that, for many open-ocean regions, these dynamical forecast systems are capable of skillfully predicting MHW onset, intensity, and duration several months in advance in both the surface and subsurface ocean42,153,154. Forecast skill, quantified by the correlation of ensemble-mean SST with that of the observations, is generally higher in the tropical and northeast extratropical Pacific, beating the skill associated with statistical (damped persistence) forecasts42,109,154. MHW forecast skill is also higher in the subsurface (0-40 m) than the surface when compared to a reanalysis product, though subsurface skill outside the tropics is primarily due to persistence153.

While in some cases dynamical forecast systems can produce skillful predictions of MHWs multiple months in advance, this is not always true. For example, 8.5-month lead forecasts initialized in July 1997 predicted an elevated likelihood of surface ocean MHW occurrence in the eastern Tropical Pacific, Gulf of Alaska, California Current, subtropical Atlantic and Indian Oceans, and the Pacific sector of the Southern Ocean during March 1998 (Fig. 5a). However, dynamical 8.5-month lead forecasts initialized in March 2013 predicted low probability of surface ocean MHWs nearly everywhere in the global ocean for November 2013 (Fig. 5b). The observed SST anomalies in March 1998 (Fig. 5c) and November 2013 (Fig. 5d) indicate that the forecasts generated in July 1997 were more accurate than the March 2013 forecasts. The accuracy of the July 1997 initialized forecast is primarily due to the development of the 1997/1998 El Niño event, as ENSO predictability imparts prediction skill to initialized forecasts of MHWs21,42. The March 2013 initialized forecast provides little indication of the development of the Blob in late-201383,156. This comparison highlights the difficulties in forecasting MHWs that are driven by stochastic atmospheric processes, like the Blob155, or energized by modes of variability not accurately captured by the models. On the other hand, surface and subsurface MHWs that are associated with ENSO variability and/or oceanic teleconnections may be predictable several months in advance21,42,109,153.

Fig. 5: Dependence of MHW forecast skill on El Niño.
figure 5

Forecasted MHW probabilities for two periods: (a) March 1998, and (b) November 2013, based on probabilistic forecasts of linearly detrended anomalies from the North American Multimodel Ensemble177 initialized 8.5 months prior (i.e., July 1997 and March 2013, respectively). White contour indicates 30% probability of occurring MHW conditions. cd Observed monthly SST anomalies (°C) for the two forecasted periods. Black contours indicate observed MHW conditions.

Statistical MHW forecasts may have similar skill as forecasts from dynamical models, while requiring substantially less computational resources. McAdam et al.153 showed that a simple statistical persistence forecast can skillfully predict the number of subsurface MHW days one season in advance in approximately half of the ocean, but it underestimates the number of events compared to the reanalysis product used as validation. More complex statistical models, including empirical-dynamical models such as Linear Inverse Models (LIMs), can be used to probe sources of predictive skill for particular regions or events. For example, LIM-based studies showed that a decadal mode of variability was a precursor for MHW growth in the Northeast Pacific Blob region41, and that predictability of MHWs off Western Australia was enhanced up to 20 months in advance by the presence of a positive Indian Ocean Dipole157.

ESM dynamical forecast systems display promising levels of forecast skill for surface and subsurface biogeochemical properties affected by MHWs, such as oxygen, acidity, or productivity158,159,160. Recent studies have explored dynamical forecast skill of ocean biogeochemical extremes. For example, Mogen et al.109 showed that a coupled model produces skillful forecasts of OAX events associated with aragonite saturation state anomalies, at lead times of up to twelve months in some regions, and further identify ENSO events as playing a key role in predicting this type of extremes. Such findings inspire efforts to include biogeochemical predictions in operational forecasting systems for MHWs.

Marine heatwaves in a changing climate

The ocean has stored more than 90% of the excess heat161 that has accumulated in the Earth System due to human-induced increases in radiative forcing agents, resulting in ocean warming that is projected by climate models to become large and widespread by the end of the century (Fig. 6a). Such slow background warming exacerbates naturally-occurring temperature excursions, resulting in increased frequency, intensity and duration of extreme SST events. Indeed, attribution studies have shown that the majority of the most impactful MHWs worldwide over recent decades could not have occurred without the influence of global warming50,57,70. In the presence of the global warming trend, climate models project large increases in the frequency, intensity, duration and spatial extent of warm temperature extremes, with the magnitude of the increase becoming progressively larger at higher warming levels57.

Fig. 6: Projected changes in MHWs in one climate model Large Ensemble.
figure 6

a Changes in mean SST (2070–2100 minus 1970–2000) based on the ensemble mean of the CESM2 large ensemble according to the SSP370 scenario. b Changes in composite MHW intensity (2070–2100 minus 1970–2000) due to internal variability divided by the intensity changes due to both changes in variability-plus-mean state. ce Histograms of area-averaged SST (°C) from the CESM2 large ensemble for (c) Arctic (poleward of 67°N), (d) western tropical Pacific (8°S–6°N, 155°E–175°W), and (e) northeast Atlantic (35°–62°N, 30°–0°W) based on all months from all 100 ensemble members during 1970–2000 (gray) and 2070–2100 (blue) after removing the ensemble-mean climatological seasonal cycle for each period. The regions considered for computing the histograms are shown by the boxes in (b). Purple histograms are the same as the blue histograms but with the mean state change (2070–2100 minus 1970–2000) added back in. The 10th and 90th percentiles of each distribution are shown as vertical solid lines, and the 50th percentile is shown as a vertical dashed line. The number in the upper right of (a) indicates the global mean ocean temperature difference (°C), while the number in the upper right of (b) indicates the fractional area (%) of values in the range −0.1 and +0.1. Adapted from Deser et al.52 © American Meteorological Society, used with permission.

In addition to the long-term ocean warming trend, climate change can also affect ocean extremes through changes in variability. An increase in mean SSTs, relative to, for example, pre-industrial levels or early historical periods, will result in a shift of the probability density function (PDF) toward larger values, enhancing the likelihood of more severe events (Fig. 6c–e). Changes in variability, however, alter the PDF’s width, which can also affect the probability of SST extremes (Fig. 6c–e). Moreover, changes in internal SST variability may be asymmetric, and lead to increased probabilities for either warm or cold extremes (Fig. 6d). The relative influence of the warming trend vs. anthropogenically-induced changes in internal variability on MHW statistics varies geographically51,52, with the long-term warming trend often accounting for more than 90% of the total changes46,50,51,52,57,143,162, as illustrated for CESM2 in Fig. 6b. Exceptions are the Arctic, where internal variability can account for 30–40% of MHW intensity changes, and the Northeast Atlantic, with values up to 80% (Fig. 6b). While separating the effects of the temperature trend and internal variability on MHW characteristics is critical for the mechanistic understanding of MHWs, and for assessing events’ predictability, this separation is challenging. The climate change trend may be nonlinear51,79, and failure to accurately account for such nonlinearity may result in an apparent change in internal variability51. Approaches used to estimate the forced trend in observations for MHW studies include the use of univariate50 or multivariate51,163,164 statistical approaches, while large ensembles can be used in the modeling context.

There are several ways by which anthropogenic forcing can alter internal climate variability. Mixed layer shoaling may occur with global warming97,152, resulting in increased mixed-layer temperatures for the same level of atmosphere-to-ocean heat exchange. The projected increase in upper-ocean stratification165 and ocean heat content162 can alter the characteristics of key large-scale drivers of MHWs. For example, increased stratification in the equatorial Pacific has been related to future enhancements of ENSO amplitude in several climate models166, while in the extra-tropics, stronger stratification will result in faster oceanic Rossby waves and shorter adjustment processes, potentially leading to reduced growth and predictability of decadal modes of variability like the PDO167. In addition, changes in extra-tropical atmospheric circulation variability driven by mean state changes and by teleconnections from changing ENSO behavior, could alter MHW characteristics through impacts on air-sea heat and momentum exchange168. Dramatic changes in Arctic sea-ice coverage and amplified Arctic warming may have been responsible for the changes in atmospheric circulation and Northeast Pacific surface heat fluxes that led to the unprecedented MHWs in that region in recent decades169. The reduction in sea-ice will also result in an increase in MHW activity near the marginal ice zone52,124.

Changes in ENSO characteristics are particularly critical for future MHWs. Consistent with the observed association between El Niño events and the enhanced likelihood of MHW occurrence40,65, multi-model large ensembles project a significant reduction in MHW areal coverage, intensity and duration during ENSO-neutral periods relative to all periods, when the mean warming component is removed52. Thus, changes in ENSO variability might significantly influence the statistics of MHWs in the future, highlighting the critical need of constraining the spread in expected ENSO changes170, and achieving more reliable future projections.

Summary and future perspectives

MHWs are an active and fast evolving area of research, and significant progress has been made in the last few years. Definitions of MHWs have been critically re-evaluated to best characterize these events and their drivers in the presence of the climate change trend, and approaches have been developed that incorporate spatial dimensions and time evolution25,54,55. On the observational side, multi-platform systems capable of providing the three-dimensional structure of MHWs in near real-time are now emerging (Box 2). Additional advances include a deepened understanding of local and remote drivers of MHWs41,112,120,171,172, explorations of subsurface MHWs and their possible structures, investigations of land-ocean and physical-biogeochemical compound events, future projections of MHWs and related uncertainty52, and evolving efforts in MHW prediction42,109,153,154. Yet, to achieve a more robust assessment of MHW predictability, additional investigations are needed to better understand large-scale drivers of MHWs in different regions and during different seasons, and their interplay in altering local processes responsible for MHW growth, evolution and persistence. MHW definitions should also be extended to reflect MHW mechanisms, and allow event characterization based on their primary drivers.

A key question for MHW prediction and projection, is whether climate models currently used for seasonal predictions and for future projections can realistically simulate MHW mechanisms beyond basic, local surface statistics (like frequency, intensity and duration). For example, can models simulate events similar to the most prominent and impactful MHWs in the historical record? Are these events driven by the same local and remote influences as in nature? Do they have similar subsurface characteristics? Assessing models’ fidelity in simulating modes of variability that can influence MHW development is also critical. Given the strong association between ENSO events and MHW occurrences40,65,90, the reliability of simulated MHWs in both present and future scenarios strongly depend on the models’ ability to realistically simulate ENSO. However, ENSO representation in climate models still shows significant biases, and its future projections vary significantly across models170, calling for an in-depth understanding of model differences and biases, and concerted efforts toward model improvement.

In order to provide forecasts that are useful for stakeholders, greater focus is needed on higher resolution global and regional models, that are able to resolve processes occurring on the shelf or at scales relevant for coastal topography (e.g., embayments, fjords, coral atolls, etc.). Regional models, which are currently under development for some regions173, should include biogeochemistry, and be used for prediction and projection applications. The availability of observations at these scales is also critical for model development and validation.

While many studies have discussed MHW impacts, this area of research is still evolving. For example, some long-lasting MHW impacts are just emerging, like the decline of the humpback whales in the North Pacific since 2014, attributed to loss of prey after the 2014-16 MHW174. Conversely, other research suggests that MHWs are not a dominant driver of change in demersal fishes over the recent decades87. Aspects in need of further research include: (1) influence of MHWs on local atmospheric extremes, like atmospheric rivers and heatwaves; (2) connections between temperature extremes and oceanic biogeochemical extremes32,133; and (3) long-term and cumulative consequences of MHWs on marine life across trophic levels, as well as assessment of recovery times in different regions. Also, given the reported impact of MHWs on air-sea CO2 fluxes69 and their documented association with cloud feedback in some regions107,108, a deeper exploration of possible MHW influences on the Earth’s carbon and radiation budgets may be important.

Observations of physical, biogeochemical and ecological quantities, at the surface and especially in the ocean subsurface, are key to all of the above aspects of MHW research. Given that MHWs can occur anytime, anywhere, concerted efforts to improve our global capacity to observe the state of the ocean both at the surface and in the subsurface are needed to properly monitor MHWs and their cascading impacts, as well as to constrain ocean reanalyses and assess model performance. Sustained observations, both globally (e.g., satellite, (deep-) Argo) or regionally (e.g., moorings), in conjunction with ocean reanalyses and long time series from empirical models, are necessary to examine long-term changes in ocean properties and robustly assess the statistics of extreme warm events relative to those long-term changes. On the other hand, systems that can be rapidly deployed for real-time monitoring (Box 2), provide not only a comprehensive characterization of individual events, but also immediate guidance to decision-makers. Such multi-platform systems, however, may not be feasible everywhere. Assessing which observations can most effectively monitor MHWs in different regions is a critical issue that the MHW community must address.

Enhanced understanding of MHWs and their impacts is essential to guide and support adaptation and mitigation strategies. The tremendous level of ecological, economical and societal losses resulting from these ocean extremes calls for urgent actions to drastically reduce greenhouse gas emissions in order to limit the devastating consequences of climate change.