Search

The ocean takes up about 93% of the global warming heat entering Earth’s climate system. In addition, the associated thermal expansion contributes substantially to sea-level rise. Hence, quantifying the oceanic heat uptake rate and its statistical significance has been a research focus. Here we use gridded ocean heat content maps to examine regional trends in ocean warming for 0–700 m depth from 1993–2019 and 1968–2019, periods based on sampling distributions. The maps are from four research groups, three based on ocean temperature alone and one combining ocean temperature with satellite altimeter sea-level anomalies. We show that use of longer periods results in larger percentages of ocean area with statistically significant warming trends and less ocean area covered by statistically significant cooling trends. We discuss relations of these patterns to climate phenomena, including the Pacific Decadal Oscillation, the Atlantic Meridional Overturning Circulation and global warming. A large proportion of anthropogenic heat energy is being taken up by ocean warming. Analysis of yearly 0–700 m ocean heat content maps from four different estimates shows that the longer the period over which regional trends are estimated, the larger the area of statistically significant warming.

The Earth system is accumulating energy due to human-induced activities. More than 90% of this energy has been stored in the ocean as heat since 1970, with ∼60% of that in the upper 700 m. Differences in upper-ocean heat content anomaly (OHCA) estimates, however, exist. Here, we use a dataset protocol for 1970–2008—with six instrumental bias adjustments applied to expendable bathythermograph (XBT) data, and mapped by six research groups—to evaluate the spatiotemporal spread in upper OHCA estimates arising from two choices: 1) those arising from instrumental bias adjustments and 2) those arising from mathematical (i.e., mapping) techniques to interpolate and extrapolate data in space and time. We also examined the effect of a common ocean mask, which reveals that exclusion of shallow seas can reduce global OHCA estimates up to 13%. Spread due to mapping method is largest in the Indian Ocean and in the eddy-rich and frontal regions of all basins. Spread due to XBT bias adjustment is largest in the Pacific Ocean within 30°N–30°S. In both mapping and XBT cases, spread is higher for 1990–2004. Statistically different trends among mapping methods are found not only in the poorly observed Southern Ocean but also in the well-observed northwest Atlantic. Our results cannot determine the best mapping or bias adjustment schemes, but they identify where important sensitivities exist, and thus where further understanding will help to refine OHCA estimates. These results highlight the need for further coordinated OHCA studies to evaluate the performance of existing mapping methods along with comprehensive assessment of uncertainty estimates.

The energy radiated by the Earth toward space does not compensate the incoming radiation from the Sun leading to a small positive energy imbalance at the top of the atmosphere (0.4-1 Wm(-2)). This imbalance is coined Earth's Energy Imbalance (EEI). It is mostly caused by anthropogenic greenhouse gas emissions and is driving the current warming of the planet. Precise monitoring of EEI is critical to assess the current status of climate change and the future evolution of climate. But the monitoring of EEI is challenging as EEI is two orders of magnitude smaller than the radiation fluxes in and out of the Earth system. Over 93% of the excess energy that is gained by the Earth in response to the positive EEI accumulates into the ocean in the form of heat. This accumulation of heat can be tracked with the ocean observing system such that today, the monitoring of Ocean Heat Content (OHC) and its long-term change provide the most efficient approach to estimate EEI. In this community paper we review the current four state-of-the-art methods to estimate global OHC changes and evaluate their relevance to derive EEI estimates on different time scales. These four methods make use of: (1) direct observations of in situ temperature; (2) satellite-based measurements of the ocean surface net heat fluxes; (3) satellite-based estimates of the thermal expansion of the ocean and (4) ocean reanalyses that assimilate observations from both satellite and in situ instruments. For each method we review the potential and the uncertainty of the method to estimate global OHC changes. We also analyze gaps in the current capability of each method and identify ways of progress for the future to fulfill the requirements of EEI monitoring. Achieving the observation of EEI with sufficient accuracy will depend on merging the remote sensing techniques with in situ measurements of key variables as an integral part of the Ocean Observing System.

Ocean warming accounts for the majority of the earth’s recent energy imbalance. Historic ocean heat content (OHC) changes are important for understanding changing climate. Calculations of OHC anomalies (OHCA) from in situ measurements provide estimates of these changes. Uncertainties in OHCA estimates arise from calculating global fields from temporally and spatially irregular data (mapping method), instrument bias corrections, and the definitions of a baseline climatology from which anomalies are calculated. To investigate sensitivity of OHCA estimates for the upper 700 m to these different factors, the same quality-controlled dataset is used by seven groups and comparisons are made. Two time periods (1970–2008 and 1993–2008) are examined. Uncertainty due to the mapping method is 16.5 ZJ for 1970–2008 and 17.1 ZJ for 1993–2008 (1 ZJ = 1 × 1021 J). Uncertainty due to instrument bias correction varied from 8.0 to 17.9 ZJ for 1970–2008 and from 10.9 to 22.4 ZJ for 1993–2008, depending on mapping method. Uncertainty due to baseline mean varied from 3.5 to 14.5 ZJ for 1970–2008 and from 2.7 to 9.8 ZJ for 1993–2008, depending on mapping method and offsets. On average mapping method is the largest source of uncertainty. The linear trend varied from 1.3 to 5.0 ZJ yr−1 (0.08–0.31 W m−2) for 1970–2008 and from 1.5 to 9.4 ZJ yr−1 (0.09–0.58 W m−2) for 1993–2008, depending on method, instrument bias correction, and baseline mean. Despite these complications, a statistically robust upper-ocean warming was found in all cases for the full time period.

Data from full-depth closely sampled hydrographic sections and Argo floats are analyzed to inform the design of a future Deep Argo array. Here standard errors of local decadal temperature trends and global decadal trends of ocean heat content and thermosteric sea level anomalies integrated from 2000 to 6000 dbar are estimated for a hypothetical 5° latitude × 5° longitude × 15-day cycle Deep Argo array. These estimates are made using temperature variances from closely spaced full-depth CTD profiles taken during hydrographic sections. The temperature data along each section are high passed laterally at a 500-km scale, and the resulting variances are averaged in 5° × 5° bins to assess temperature noise levels as a function of pressure and geographic location. A mean global decorrelation time scale of 62 days is estimated using temperature time series at 1800 dbar from Argo floats. The hypothetical Deep Argo array would be capable of resolving, at one standard error, local trends from −1 in the quiescent abyssal North Pacific to about 26 m °C decade−1 below 2000 dbar along 50°S in the energetic Southern Ocean. Larger decadal temperature trends have been reported previously in these regions using repeat hydrographic section data, but those very sparse data required substantial spatial averaging to obtain statistically significant results. Furthermore, the array would provide decadal global ocean heat content trend estimates from 2000 to 6000 dbar with a standard error of ±3 TW, compared to a trend standard error of ±17 TW from a previous analysis of repeat hydrographic data.

Characteristics of eddies in the Gulf of Alaska are assessed from January 2003 through April 2012. Ensemble statistics for eddy subsurface water properties on isopycnals are computed using temperature and salinity profiles from Argo profiling floats located within eddies, which are identified in sea-surface height using objective techniques. Ninety cyclonic and 154 anticyclonic eddies are identified during this period. The anticyclonic eddies are strongly nonlinear and exhibit significant warm subsurface temperature anomalies and associated salty anomalies on isopycnals while no clear distinguishing subsurface anomalies on isopycnals are detected in association with the cyclonic eddies. Heat and freshwater fluxes for the eddies are estimated from integrations in depth coordinates. The anticyclonic eddies transport heat both westward off the continental shelf into the Subarctic Gyre and westward within the Alaskan Stream. However, they transport salt into the Subarctic Gyre and freshwater within the Alaskan Stream. In both pathways eddy heat and freshwater transport show possible year-to-year fluctuations, varying from 0 to 50.4 × 1018 J a−1 and −16.8 to +7.4 km3 a−1, respectively. The anticyclonic eddies are capped by relatively fresh water year-round.

Ocean heat content anomalies are analyzed from 1950 to 2011 in five distinct depth layers (0–100, 100–300, 300–700, 700–900, and 900–1800 m). These layers correspond to historic increases in common maximum sampling depths of ocean temperature measurements with time, as different instruments—mechanical bathythermograph (MBT), shallow expendable bathythermograph (XBT), deep XBT, early sometimes shallower Argo profiling floats, and recent Argo floats capable of worldwide sampling to 2000 m—have come into widespread use. This vertical separation of maps allows computation of annual ocean heat content anomalies and their sampling uncertainties back to 1950 while taking account of in situ sampling advances and changing sampling patterns. The 0–100-m layer is measured over 50% of the globe annually starting in 1956, the 100–300-m layer starting in 1967, the 300–700-m layer starting in 1983, and the deepest two layers considered here starting in 2003 and 2004, during the implementation of Argo. Furthermore, global ocean heat uptake estimates since 1950 depend strongly on assumptions made concerning changes in undersampled or unsampled ocean regions. If unsampled areas are assumed to have zero anomalies and are included in the global integrals, the choice of climatological reference from which anomalies are estimated can strongly influence the global integral values and their trend: the sparser the sampling and the bigger the mean difference between climatological and actual values, the larger the influence.

Global climate change results from a small yet persistent imbalance between the amount of sunlight absorbed by the Earth and the thermal radiation emitted back to space. A revised analysis of measured changes in the net radiation imbalance at the top of the atmosphere, and the ocean heat content to a depth of 1,800 m, suggests that these two sets of observations are consistent within error margins.

The relative significance of short multi-year linear trends in the global integral of 0–700 m ocean heat content anomaly (OHCA) is investigated by examining the overlapping segments of the 16-year OHCA curve from Lyman et al. (Nature 465:334–337, 2010). Segments of 4 years and less are found not to be significantly different from each other or from 0 at the 90% confidence interval. Likewise, short 5- to 7-year segments are not statistically different from each other. Ten-year and longer trends are significant and provide a useful comparison for satellite observations of the radiation imbalance at the top of the atmosphere.

Two significant instrument biases have been identified in the in situ profile data used to estimate globally integrated upper-ocean heat content. A large cold bias was discovered in a small fraction of Argo floats along with a smaller but more prevalent warm bias in expendable bathythermograph (XBT) data. These biases appear to have caused the bulk of the upper-ocean cooling signal reported by Lyman et al. between 2003 and 2005. These systematic data errors are significantly larger than sampling errors in recent years and are the dominant sources of error in recent estimates of globally integrated upper-ocean heat content variability. The bias in the XBT data is found to be consistent with errors in the fall-rate equations, suggesting a physical explanation for that bias. With biased profiles discarded, no significant warming or cooling is observed in upper-ocean heat content between 2003 and 2006.

The effects of irregular in situ ocean sampling on estimates of annual globally integrated upper ocean heat content anomalies (OHCA) are investigated for sampling patterns from 1955 to 2006. An analytical method is presented for computing the effective area covered by an objective map for any given in situ sampling distribution. To evaluate the method, appropriately scaled sea surface height (SSH) anomaly maps from Archiving, Validation, and Interpretation of Satellite Oceanographic data (AVISO) are used as a proxy for OHCA from 1993 to 2006. Use of these proxy data demonstrates that the simple area integral (SI) of such an objective map for sparse datasets does not agree as well with the actual integral as the weighted integral (WI), defined as the simple integral weighted by the ratio of the total area over the “observed” area. From 1955 to 1966, in situ ocean sampling is inadequate to estimate accurately annual global integrals of the proxy upper OHCA. During this period, the SI for the sampling pattern of any given year underestimates the 13-yr trend in proxy OHCA from 1993 to 2006 by around 70%, and confidence limits for the WI are often very large. From 1967 to 2003 there appear to be sufficient data to estimate annual global integrals. Limited by the constraints of this analysis, the SI for any given year’s sampling pattern still underestimates the 1993–2006 13-yr trend in the proxy by around 30%, but the WI matches the trend well with small confidence limits. For 2004 through 2006 in situ sampling, with near-global in situ Argo data coverage, the 1993–2006 13-yr trend in the proxy is equally well represented by the SI or WI.

The evolution of ocean temperature measurement systems is presented with a focus on the development and accuracy of two critical devices in use today (expendable bathythermographs and conductivity-temperature-depth instruments used on Argo floats). A detailed discussion of the accuracy of these devices and a projection of the future of ocean temperature measurements are provided. The accuracy of ocean temperature measurements is discussed in detail in the context of ocean heat content, Earth's energy imbalance, and thermosteric sea level rise. Up-to-date estimates are provided for these three important quantities. The total energy imbalance at the top of atmosphere is best assessed by taking an inventory of changes in energy storage. The main storage is in the ocean, the latest values of which are presented. Furthermore, despite differences in measurement methods and analysis techniques, multiple studies show that there has been a multidecadal increase in the heat content of both the upper and deep ocean regions, which reflects the impact of anthropogenic warming. With respect to sea level rise, mutually reinforcing information from tide gauges and radar altimetry shows that presently, sea level is rising at approximately 3 mm yr-1 with contributions from both thermal expansion and mass accumulation from ice melt. The latest data for thermal expansion sea level rise are included here and analyzed. Key Points Oceanographic techniques and analysis have improved over many decadesThese improvements allow more accurate Earth-energy balance estimatesUnderstanding of ocean heat content and sea-level rise has also increased ©2013. American Geophysical Union. All Rights Reserved.

In recent decades, over nine-tenths of Earth's top-of-the-atmosphere energy imbalance has been stored in the ocean, which is rising as it warms. Combining satellite sea-level data with ocean mass data or model results allows insights into ocean warming.

This community white paper outlines the requirements of the future observing system necessary for measuring and advancing understanding of global ocean heat uptake and heat content variability, with an emphasis on the in situ observing system. We review the progress made in observation-based estimates of ocean heat uptake since Ocean Obs '99 and propose a future observational strategy. Some of the key scientific questions addressed are: 1. What future observations are required to monitor global ocean heat content? 2. How has new technology improved our ability to make estimates of ocean heat uptake? 3. What are the current estimates of global and regional ocean heat uptake and what are the uncertainties? 4. What is the impact of instrumental biases and gridding methodology on estimates of ocean heat uptake?

The upper ocean acts as a giant heat sink and has absorbed the majority of excess energy generated by anthropogenic greenhouse gasses. This makes ocean heat content, potentially, a key indicator of climate change. But to be useful for evaluating the global energy balance and as a constraint on climate models, the measurement uncertainties of such a key indicator need to be well understood. At present the magnitude of the oceanic heat uptake is highly uncertain, with patterns of inter-annual variability in particular differing among estimates. In a major international collaboration, Lyman et al. compare the available upper-ocean heat content anomaly curves and examine the sources of uncertainly attached to them — including the difficulties in correcting bias in expendable bathythermograph data. They find that, uncertainties notwithstanding, there is clear and robust evidence for a warming trend of 0.64 watts per square metre between 1993 and 2008. The upper 300 m of the world's oceans act as a giant heat sink and have absorbed the majority of the excess energy generated by anthropogenic greenhouse gases. But the magnitude of the oceanic heat uptake is uncertain, and differing estimates have led to questions regarding the closure of the global energy budget. Here, a comparison of ocean heat content estimates is presented; the conclusion is that a robust warming of 0.64 W m−2 occurred from 1993 to 2008. A large (∼1023 J) multi-decadal globally averaged warming signal in the upper 300 m of the world’s oceans was reported roughly a decade ago1 and is attributed to warming associated with anthropogenic greenhouse gases2,3. The majority of the Earth’s total energy uptake during recent decades has occurred in the upper ocean3, but the underlying uncertainties in ocean warming are unclear, limiting our ability to assess closure of sea-level budgets4,5,6,7, the global radiation imbalance8 and climate models5. For example, several teams have recently produced different multi-year estimates of the annually averaged global integral of upper-ocean heat content anomalies (hereafter OHCA curves) or, equivalently, the thermosteric sea-level rise5,9,10,11,12,13,14,15,16. Patterns of interannual variability, in particular, differ among methods. Here we examine several sources of uncertainty that contribute to differences among OHCA curves from 1993 to 2008, focusing on the difficulties of correcting biases in expendable bathythermograph (XBT) data. XBT data constitute the majority of the in situ measurements of upper-ocean heat content from 1967 to 2002, and we find that the uncertainty due to choice of XBT bias correction dominates among-method variability in OHCA curves during our 1993–2008 study period. Accounting for multiple sources of uncertainty, a composite of several OHCA curves using different XBT bias corrections still yields a statistically significant linear warming trend for 1993–2008 of 0.64 W m-2 (calculated for the Earth’s entire surface area), with a 90-per-cent confidence interval of 0.53–0.75 W m-2.