Trends In Ocean Heat Content Research Articles

IAPv4 ocean temperature and ocean heat content gridded dataset

Abstract. Ocean observational gridded products are vital for climate monitoring, ocean and climate research, model evaluation, and supporting climate mitigation and adaptation measures. This paper describes the 4th version of the Institute of Atmospheric Physics (IAPv4) ocean temperature and ocean heat content (OHC) objective analysis product. It accounts for recent developments in quality control (QC) procedures, climatology, bias correction, vertical and horizontal interpolation, and mapping and is available for the upper 6000 m (119 levels) since 1940 (more reliable after ∼ 1957) for monthly and 1°×1° temporal and spatial resolutions. IAPv4 is compared with the previous version, IAPv3, and with the other data products, sea surface temperatures (SSTs), and satellite observations. It has a slightly stronger long-term upper 2000 m OHC increase than IAPv3 for 1955–2023, mainly because of newly developed bias corrections. The IAPv4 0–2000 m OHC trend is also higher during 2005–2023 than IAPv3, mainly because of the QC process update. The uppermost level of IAPv4 is consistent with independent SST datasets. The month-to-month OHC variability for IAPv4 is desirably less than IAPv3 and the other OHC products investigated in this study, the trend of ocean warming rate (i.e., warming acceleration) is more consistent with the net energy imbalance at the top of the atmosphere than IAPv3, and the sea level budget can be closed within uncertainty. The gridded product is freely accessible at https://doi.org/10.12157/IOCAS.20240117.002 for temperature data (Cheng et al., 2024a) and at https://doi.org/10.12157/IOCAS.20240117.001 for ocean heat content data (Cheng et al., 2024b).

Skill assessment of seasonal forecasts of ocean variables

There is growing demand for seasonal forecast products for marine applications. The availability of consistent and sufficiently long observational records of ocean variables permits the assessment of the spatial distribution of the skill of ocean variables from seasonal forecasts. Here we use state-of-the-art temporal records of sea surface temperature (SST), sea surface height (SSH) and upper 300m ocean heat content (OHC) to quantify the distribution of skill, up to 2 seasons ahead, of two operational seasonal forecasting systems contributing to the seasonal multi-model of the Copernicus Climate Change Services (C3S). This study presents the spatial distribution of the skill of the seasonal forecast ensemble mean in terms of anomaly correlation and root mean square error and compares it to the persistence and climatological benchmarks. The comparative assessment of the skill among variables sheds light on sources/limits of predictability at seasonal time scales, as well as the nature of model errors. Beyond these standard verification metrics, we also evaluate the ability of the models to represent the observed long-term trends. Results show that long-term trends contribute to the skill of seasonal forecasts. Although the forecasts capture the long-term trends in general, some regional aspects remain challenging. Part of these errors can be attributed to specific aspects of the ocean initialization, but others, such as the overestimation of the warming in the Eastern Pacific are also influenced by model error. Skill gains can be obtained by improving the trend representation in future forecasting systems. In the meantime, a forecast calibration procedure that corrects the linear trends can produce substantial skill gains. The results show that calibrated seasonal forecasts beat both the climatological and persistence benchmark almost at every location for all initial dates and lead times. Results demonstrate the value of the seasonal forecasts for marine applications and highlight the importance of representing the decadal variability and trends in ocean heat content and sea level.

Interhemispheric Contrasts of Ocean Heat Content Change Reveals Distinct Fingerprints of Anthropogenic Climate Forcings

AbstractDuring recent decades, both greenhouse gases (GHGs) and anthropogenic aerosols (AAs) drove major changes in the Earth's energy imbalance. However, their respective fingerprints in changes to ocean heat content (OHC) have been difficult to isolate and detect when global or hemispheric averages are used. Based on a pattern recognition analysis, we show that AAs drive an interhemispheric asymmetry within the 20°‐35° latitude band in historical OHC change due to the southward shift of the atmospheric and ocean circulation system. This forced pattern is distinct from the GHG‐induced pattern, which dominates the asymmetry in higher latitudes. Moreover, it is found that this significant aerosol‐forced OHC trend pattern can only be captured in analyzed periods of 20 years or longer and including 1975–1990. Using these distinct spatiotemporal characteristics, we show that the fingerprint of aerosol climate forcing in ocean observations can be distinguished from both the stronger GHG‐induced signals and internal variability.

Rectified Ocean Heat Uptake from Oscillatory Surface Forcing

Abstract Ocean heat uptake is asymmetric with respect to the sign of radiative forcing. It is already known that surface cooling anomalies penetrate into the ocean faster than surface warming anomalies. Because of this asymmetry, the time-variable component of radiative forcing can induce a long-term, rectified cooling trend in ocean heat content, which is this work’s primary focus. Here, we explore this asymmetry and rectification on global and interannual scales, its implications, and its possible dependence on model parameters. We do so using a full-complexity global ocean–sea ice general circulation model and an idealized one-dimensional vertical mixing model, both forced with idealized abrupt and oscillatory surface forcing anomalies. In both models, the ocean heat uptake response to an abrupt cooling perturbation is larger than an equal-magnitude warming perturbation. This asymmetry is shown to be larger when the background vertical diffusivity is smaller, and is therefore likely model dependent. Sinusoidal oscillatory forcing with zero time mean induces a rectified cooling trend in both models, whose magnitude depends on both the diffusivity and the frequency of the interannual oscillatory forcing. The net rectified cooling can reach a rate of approximately −0.11 W m−2, which is substantial relative to the estimated anthropogenic warming rate of 0.47 W m−2 (von Schuckmann et al.). We discuss this rectification effect in the context of volcanic forcing in climate models, whose time-variable component may cause model-dependent ocean cooling in CMIP6 historical simulations, reaching 5%–30% the size of the total impact of volcanic forcing in our one-dimensional model. Correcting for this cooling may help reduce uncertainties in modeled ocean heat content evolution.

Heat Storage in the Upper Indian Ocean: The Role of Wind-Driven Redistribution

Abstract The heat content in the Indian Ocean has been increasing owing to anthropogenic greenhouse warming. Yet, where and how the anthropogenic heat is stored in the Indian Ocean have not been comprehended. Analysis of various observational and model-based datasets since the 1950s reveals a robust spatial pattern of the 0–700 m ocean heat content trend (ΔOHC), with enhanced warming in the subtropical southern Indian Ocean (SIO) but weak to minimal warming in the tropical Indian Ocean (TIO). The meridional temperature gradient between the TIO and SIO declined by 16.4% ± 7.5% during 1958–2014. The heat redistribution driven by time-varying surface winds plays a crucial role in shaping this ΔOHC pattern. Sensitivity experiments using a simplified ocean dynamical model suggest that changes in surface winds over the Indian Ocean, particularly those of the SIO, caused a convergence trend in the upper SIO and a divergence trend in the upper TIO. These wind changes primarily include the enhancements of westerlies in the Southern Ocean and the subtropical anticyclone in the SIO. Albeit with systematic biases, the ΔOHC pattern and surface wind changes simulated by phase 6 of the Coupled Model Intercomparison Project (CMIP6) models broadly resemble the observation and highlight the essence of external forcing in causing these changes. This heat storage pattern is projected to persist in the model-projected future, potentially impacting future climate.

Hierarchical Bayesian modeling of ocean heat content and its uncertainty

The accurate quantification of changes in the heat content of the world’s oceans is crucial for our understanding of the effects of increasing greenhouse gas concentrations. The Argo program, consisting of Lagrangian floats that measure vertical temperature profiles throughout the global ocean, has provided a wealth of data from which to estimate ocean heat content. However, creating a globally consistent statistical model for ocean heat content remains challenging due to the need for a globally valid covariance model that can capture complex nonstationarity. In this paper, we develop a hierarchical Bayesian Gaussian process model that uses kernel convolutions with cylindrical distances to allow for spatial nonstationarity in all model parameters while using a Vecchia process to remain computationally feasible for large spatial datasets. Our approach can produce valid credible intervals for globally integrated quantities that would not be possible using previous approaches. These advantages are demonstrated through the application of the model to Argo data, yielding credible intervals for the spatially varying trend in ocean heat content that accounts for both the uncertainty induced from interpolation and from estimating the mean field and other parameters. Through cross-validation, we show that our model outperforms an out-of-the-box approach as well as other simpler models. The code for performing this analysis is provided as the R package BayesianOHC.

A Comparative Study of the Argo‐Era Ocean Heat Content Among Four Different Types of Data Sets

AbstractWe conducted a comparative study of ocean heat content (OHC) in the top 2,000 m during the Argo‐era using 12 latest and representative global ocean data sets. The differences in the global and basins‐wide OHC trends were minor among the observation‐based data sets, and remarkable among the ocean reanalyzes (RAs). Some RAs might exhibit much higher or lower basins‐wide warming rates than the observation‐based data sets. In the top 700 m, RAs suggested similar large‐scale warming and cooling patterns, in agreement with the observation‐based data sets. Below 700 m, the major warming and cooling features were significantly different between RAs and the observation‐based data sets. All data sets suffered from relatively larger uncertainties in the highly dynamic regions. Special caution is suggested when estimating the OHC using only a single data set, especially a RA. Differences of RAs' OHC from observation‐based data sets were significantly reduced when considering their ensemble mean, to be further confirmed with a larger sample of data sets.

Imprint of intrinsic ocean variability on decadal trends of regional sea level and ocean heat content using synthetic profiles

The global ocean is warming and has absorbed 90% of the Earth Energy Imbalance over 2010–2018 leading to global mean sea level rise. Both ocean heat content (OHC) and sea level trends show large regional deviations from their global means. Both quantities have been estimated from in-situ observations for years. However, in-situ profile coverage is spatially uneven, leading to uncertainties when assessing both OHC and sea level trends, especially at regional scale. Recently, a new possible driver of regional sea level and OHC trends has been highlighted using eddy-permitting ensemble ocean simulations over multiple decades: non-linear ocean processes produce chaotic fluctuations, which yield random contributions to regional decadal OHC and sea level trends. In-situ measurements capture a combination of the atmospherically-forced response and this intrinsic ocean variability. It is therefore important to understand the imprint of the chaotic ocean variability recorded by the in-situ measurement sampling in order to assess its impact and associated uncertainty on regional budgets. A possible approach to investigate this problem is to use a set of synthetic in-situ-like profiles extracted from an ensemble of forced ocean simulations started from different states and integrated with the same atmospheric forcing. Comparisons between the original ensemble outputs and the remapped, subsampled, in-situ-like profiles elucidate the contribution of chaotic ocean variability to OHC and regional sea level trends. Our results show that intrinsic variability may be large in eddy-active regions in the gridded model outputs, and remains substantial when using the in-situ sampling-based estimates. Using the latter, the same result is also found on large scales, for which atmospheric forcing has been identified as the main driver. Our results suggest accounting for this intrinsic ocean variability when assessing regional OHC and sea level trend budgets on decadal time scales.

Variability and Drivers of Ocean Temperature Extremes in a Warming Western Boundary Current

Abstract Western boundary current (WBC) extensions such as the East Australian Current (EAC) southern extension are warming 2–3 times faster than the global average. However, there are nuances in the spatial and temporal variability of the warming that are not well resolved in climate models. In addition, the physical drivers of ocean heat content (OHC) extremes are not well understood. Here, using a high-resolution ocean model run for multiple decades, we show nonuniform warming trends in OHC in the EAC, with strong positive trends in the southern extension region (∼36°–38°S) but negative OHC trends equatorward of 33°S. The OHC variability in the EAC is associated with the formation of anticyclonic eddies, which is modulated by transport ∼880 km upstream (EAC mode) and the westward propagation of Rossby waves (eddy mode). Diagnosing the drivers of temperature extremes has implications for predictability both in the EAC and in WBCs more broadly, where ocean warming is already having considerable ecological impacts.

Trends, variability and predictive skill of the ocean heat content in North Atlantic: an analysis with the EC-Earth3 model

This study investigates linear trends, variability and predictive skill of the upper ocean heat content (OHC) in the North Atlantic basin. This is a region where strong decadal variability superimposes the externally forced trends, introducing important differences in the local warming rates and leading in the case of the Central Subpolar North Atlantic to an overall long-term cooling. Our analysis aims to better understand these regional differences, by investigating how internal and forced variability contribute to local trends, exploring also their role on the local prediction skill. The analysis combines the study of three ocean reanalyses to document the uncertainties related to observations with two sets of CMIP6 experiments performed with the global coupled climate model EC-Earth3: a historical ensemble to characterise the forced signals, and a retrospective decadal prediction system to additionally characterise the contributions from internal climate variability. Our results show that internal variability is essential to understand the spatial pattern of North Atlantic OHC trends, contributing decisively to the local trends and providing high levels of predictive skill in the Eastern Subpolar North Atlantic and the Irminger and Iceland Seas, and to a lesser extent in the Labrador Sea. Skill and trends in other areas like the Subtropical North Atlantic, or the Gulf Stream Extension are mostly externally forced. Large observational and modeling uncertainties affect the trends and interannual variability in the Central Subpolar North Atlantic, the only region exhibiting a cooling during the study period, uncertainties that might explain the very poor local predictive skill.

Observational Constraint on Greenhouse Gas and Aerosol Contributions to Global Ocean Heat Content Changes

AbstractObservations and climate models are combined to identify an anthropogenic warming signature in the upper ocean heat content (OHC) changes since 1971. We apply a new detection and attribution analysis developed by Ribes et al. that uses a symmetric treatment of the magnitude and the pattern of the climate response to each radiative forcing. A first estimate of the OHC response to natural, anthropogenic, greenhouse gas, and other forcings is derived from a large ensemble of CMIP5 simulations. Observational datasets from historical reconstructions are then used to constrain this estimate. A spatiotemporal observational mask is applied to compare simulations with actual observations and to overcome reconstruction biases. Results on the 0–700-m layer from 1971 to 2005 show that the global OHC would have increased since 1971 by 2.12 ± 0.21 × 107 J m−2 yr−1 in response to GHG emissions alone. But this has been compensated for by other anthropogenic influences (mainly aerosol), which induced an OHC decrease of 0.84 ± 0.18 × 107 J m−2 yr−1. The natural forcing has induced a slight global OHC decrease since 1971 of 0.13 ± 0.09 × 107 J m−2 yr−1. Compared to previous studies we have separated the effect of the GHG forcing from the effect of the other anthropogenic forcing on OHC changes. This has been possible by using a new detection and attribution (D&A) method and by analyzing simultaneously the global OHC trends over 1957–80 and over 1971–2005. This bivariate method takes advantage of the different time variation of the GHG forcing and the aerosol forcing since 1957 to separate both effects and reduce the uncertainty in their estimates.

Towards quantifying uncertainty in ocean heat content changes using synthetic profiles

Observational estimates of global ocean heat content (OHC) change are used to assess Earth’s energy imbalance over the 20th Century. However, intercomparison studies show that the mapping methods used to interpolate sparse ocean temperature profile data are a key source of uncertainty. We present a new approach to assessing OHC mapping methods using ‘synthetic profiles’ generated from a state-of-the-art global climate model simulation. Synthetic profiles have the same sampling characteristics as the historical ocean temperature profile data but are based on model simulation data. Mapping methods ingest these data in the same way as they would real observations, but the resultant mapped fields can be compared to a model simulation ‘truth’. We use this approach to assess two mapping methods that are used routinely for climate monitoring and initialisation of decadal forecasts. The introduction of the Argo network of autonomous profiling floats during the 2000s drives clear improvements in the ability of these methods to reconstruct the variability and spatial structure of OHC changes. At depths below 2000 m, both methods underestimate the magnitude of the simulated ocean warming signal. Temporal variability and trends in OHC are better captured in the better-observed northern hemisphere than in the southern hemisphere. At all depths, the sampling characteristics of the historical data introduces some spurious variability in the estimates of global OHC on sub-annual to multi-annual timescales. However, many of the large scale spatial anomalies, especially in the upper ocean, are successfully reconstructed even with sparse observations from the 1960s, demonstrating the potential to construct historical ocean analyses for assessing decadal predictions. The value of using accurate global covariances for data-poor periods is clearly seen. The results of this ‘proof-of-concept’ study are encouraging for gaining further insights into the capabilities and limitations of different mapping methods and for quantifying uncertainty in global OHC estimates.

Mechanisms underlying recent decadal changes in subpolar North Atlantic Ocean heat content

AbstractThe subpolar North Atlantic (SPNA) is subject to strong decadal variability, with implications for surface climate and its predictability. In 2004–2005, SPNA decadal upper ocean and sea‐surface temperature trends reversed from warming during 1994–2004 to cooling over 2005–2015. This recent decadal trend reversal in SPNA ocean heat content (OHC) is studied using a physically consistent, observationally constrained global ocean state estimate covering 1992–2015. The estimate's physical consistency facilitates quantitative causal attribution of ocean variations. Closed heat budget diagnostics reveal that the SPNA OHC trend reversal is the result of heat advection by midlatitude ocean circulation. Kinematic decompositions reveal that changes in the deep and intermediate vertical overturning circulation cannot account for the trend reversal, but rather ocean heat transports by horizontal gyre circulations render the primary contributions. The shift in horizontal gyre advection reflects anomalous circulation acting on the mean temperature gradients. Maximum covariance analysis (MCA) reveals strong covariation between the anomalous horizontal gyre circulation and variations in the local wind stress curl, suggestive of a Sverdrup response. Results have implications for decadal predictability.

A global probabilistic study of the ocean heat content low‐frequency variability: Atmospheric forcing versus oceanic chaos

AbstractA global 1/4° ocean/sea ice 50‐member ensemble simulation is used to disentangle the low‐frequency imprints of the atmospherically forced oceanic variability and of the chaotic intrinsic oceanic variability (IOV) on the large‐scale (10° × 10°) ocean heat content (OHC) between 1980 and 2010. The IOV explains most of the interannual‐to‐decadal large‐scale OHC variance over substantial fractions of the global ocean area that increase with depth: 9%, 22%, and 31% in the 0–700 m, 700–2000 m and 2000 m bottom layers, respectively. Such areas concern principally eddy‐active regions, mostly found in the Southern Ocean and in western boundary current extensions, and also concern the subtropical gyres at intermediate and deep levels. The oceanic chaos may also induce random multidecadal fluctuations so that large‐scale regional OHC trends computed on the 1980–2010 period cannot be unambiguously attributed to the atmospheric forcing in several oceanic basins at various depths. These results are likely to raise detection and attribution issues from real observations.

Accuracy of Global Upper Ocean Heat Content Estimation Expected from Present Observational Data Sets

The simplest global mapping method and dense data coverage for the global oceans by the latest observation network ensure an estimate of global ocean heat content (OHC) within a satisfactory uncertainty for the last 60 years. The observational database conditionally presented a level high enough for practical use for the global OHC estimation when applying bias corrections of expendable bathythermograph, assuming that the other severe observational biases are not included in the database. Uncertainties in annual global mean temperatures averaged vertically from the surface to 1,500 m are within 0.01 K for the period from 1955 onward, when only sampling errors are taken into account. Those in annual mean global OHC of an improved objective analysis for 0-1,500 m depth is 16ZJ on average throughout the period. Compared to previous studies, the new objective analysis provides a higher estimation of the global 0-1,500 m OHC trend for a longer period from 1955 to 2015, which is an increase of 350 ± 57ZJ with a 95% confidence interval.

Insights into Earth’s Energy Imbalance from Multiple Sources

Abstract The current Earth’s energy imbalance (EEI) can best be estimated from changes in ocean heat content (OHC), complemented by top-of-atmosphere (TOA) radiation measurements and an assessment of the small non-ocean components. Sustained observations from the Argo array of autonomous profiling floats enable near-global estimates of OHC since 2005, which reveal considerable cancellation of variations in the upper 300 m. An analysis of the monthly contributions to EEI from non-ocean components (land and ice) using the Community Earth System Model (CESM) Large Ensemble reveals standard deviations of 0.3–0.4 W m−2 (global); largest values occur in August, but values are below 0.75 W m−2 greater than 95% of the time. Global standard deviations of EEI of 0.64 W m−2 based on top-of-atmosphere observations therefore substantially constrain ocean contributions, given by the tendencies of OHC. Instead, monthly standard deviations of many Argo-based OHC tendencies are 6–13 W m−2, and nonphysical fluctuations are clearly evident. It is shown that an ocean reanalysis with multivariate dynamical data assimilation features much better agreement with TOA radiation, and 44% of the vertically integrated short-term OHC trend for 2005–14 of 0.8 ± 0.2 W m−2 (globally) occurs below 700-m depth. Largest warming occurs from 20° to 50°S, especially over the southern oceans, and near 40°N in all ocean analyses. The EEI is estimated to be 0.9 ± 0.3 W m−2 for 2005–14.

Ocean heat content variability and change in an ensemble of ocean reanalyses

Accurate knowledge of the location and magnitude of ocean heat content (OHC) variability and change is essential for understanding the processes that govern decadal variations in surface temperature, quantifying changes in the planetary energy budget, and developing constraints on the transient climate response to external forcings. We present an overview of the temporal and spatial characteristics of OHC variability and change as represented by an ensemble of dynamical and statistical ocean reanalyses (ORAs). Spatial maps of the 0–300 m layer show large regions of the Pacific and Indian Oceans where the interannual variability of the ensemble mean exceeds ensemble spread, indicating that OHC variations are well-constrained by the available observations over the period 1993–2009. At deeper levels, the ORAs are less well-constrained by observations with the largest differences across the ensemble mostly associated with areas of high eddy kinetic energy, such as the Southern Ocean and boundary current regions. Spatial patterns of OHC change for the period 1997–2009 show good agreement in the upper 300 m and are characterized by a strong dipole pattern in the Pacific Ocean. There is less agreement in the patterns of change at deeper levels, potentially linked to differences in the representation of ocean dynamics, such as water mass formation processes. However, the Atlantic and Southern Oceans are regions in which many ORAs show widespread warming below 700 m over the period 1997–2009. Annual time series of global and hemispheric OHC change for 0–700 m show the largest spread for the data sparse Southern Hemisphere and a number of ORAs seem to be subject to large initialization ‘shock’ over the first few years. In agreement with previous studies, a number of ORAs exhibit enhanced ocean heat uptake below 300 and 700 m during the mid-1990s or early 2000s. The ORA ensemble mean (±1 standard deviation) of rolling 5-year trends in full-depth OHC shows a relatively steady heat uptake of approximately 0.9 ± 0.8 W m−2 (expressed relative to Earth’s surface area) between 1995 and 2002, which reduces to about 0.2 ± 0.6 W m−2 between 2004 and 2006, in qualitative agreement with recent analysis of Earth’s energy imbalance. There is a marked reduction in the ensemble spread of OHC trends below 300 m as the Argo profiling float observations become available in the early 2000s. In general, we suggest that ORAs should be treated with caution when employed to understand past ocean warming trends—especially when considering the deeper ocean where there is little in the way of observational constraints. The current work emphasizes the need to better observe the deep ocean, both for providing observational constraints for future ocean state estimation efforts and also to develop improved models and data assimilation methods.

Heat content variability in the North Atlantic Ocean in ocean reanalyses.

Warming of the North Atlantic Ocean from the 1950s to 2012 is analyzed on neutral density surfaces and vertical levels in the upper 2000 m. Three reanalyses and two observational data sets are compared. The net gain of 5 × 1022 J in the upper 2000 m is roughly 30% of the global ocean warming over this period. Upper ocean heat content (OHC) is dominated in most regions by heat transport convergence without widespread changes in the potential temperature/salinity relation. The heat convergence is associated with sinking of midthermocline isopycnals, with maximum sinking occurring at potential densities σ0 = 26.4−27.3, which contain subtropical mode waters. Water masses lighter than σ0 = 27.3 accumulate heat by increasing their volume, while heavier waters lose heat by decreasing their volume. Spatially, the OHC trend is nonuniform: the low latitudes, 0–30°N are warming steadily while large multidecadal variability occurs at latitudes 30–65°N.Key Points Heat content change dominated by heat transport convergence Due to widespread sinking trend of midthermocline isopycnals over 50+ years

Earth's energy imbalance since 1960 in observations and CMIP5 models.

Observational analyses of running 5 year ocean heat content trends (Ht) and net downward top of atmosphere radiation (N) are significantly correlated (r ∼ 0.6) from 1960 to 1999, but a spike in Ht in the early 2000s is likely spurious since it is inconsistent with estimates of N from both satellite observations and climate model simulations. Variations in N between 1960 and 2000 were dominated by volcanic eruptions and are well simulated by the ensemble mean of coupled models from the Fifth Coupled Model Intercomparison Project (CMIP5). We find an observation-based reduction in N of − 0.31 ± 0.21 W m−2 between 1999 and 2005 that potentially contributed to the recent warming slowdown, but the relative roles of external forcing and internal variability remain unclear. While present-day anomalies of N in the CMIP5 ensemble mean and observations agree, this may be due to a cancelation of errors in outgoing longwave and absorbed solar radiation.Key PointsObserved maximum in ocean heat content trend in early 2000s is likely spuriousNet incoming radiation (N) reduced by 0.31 ± 0.21 W m−2 during the warming pausePresent-day estimates of N may contain opposing errors in radiative components

Influences of the Choice of Climatology on Ocean Heat Content Estimation

AbstractThe choice of climatology is an essential step in calculating the key climate indicators, such as historical ocean heat content (OHC) change. The anomaly field is required during the calculation and is obtained by subtracting the climatology from the absolute field. The climatology represents the ocean spatial variability and seasonal circle. This study found a considerable weaker long-term trend when historical climatologies (constructed by using historical observations within a long time period, i.e., 45 yr) were used rather than Argo-period climatologies (i.e., constructed by using observations during the Argo period, i.e., since 2004). The change of the locations of the observations (horizontal sampling) during the past 50 yr is responsible for this divergence, because the ship-based system pre-2000 has insufficient sampling of the global ocean, for instance, in the Southern Hemisphere, whereas this area began to achieve full sampling in this century by the Argo system. The horizontal sampling change leads to the change of the reference time (and reference OHC) when the historical-period climatology is used, which weakens the long-term OHC trend. Therefore, Argo-period climatologies should be used to accurately assess the long-term trend of the climate indicators, such as OHC.