Global Ocean Heat Content Anomalies based on Argo data

NOTE for users: please use the latest version of the product at https://zenodo.org/doi/10.5281/zenodo.10182972. Ocean Heat Content Anomalies (OHCA) are calculated subtracting the mean over the period 2005-2020 from the monthly time series. Yearly OHCA time series are then calculated. OHC fields are mapped using locally stationary Gaussian processes with data-driven decorrelation scales (Kuusela and Stein, 2018). A linear time trend was included in the estimate of the mean field (along with spatial terms and harmonics for the annual cycle). In the present version, mapping is done in latitude and longitude with monthly subsets of data (a future version will add time to the mapping). Mapping is done separately for different vertical sections. Different vertical sections are combined to estimate: 1. Global OHC timeseries (e.g., for level 0-2000m: GCOS_0000_2000_OHCA_J_m2_oc, for OHC in J/m2; GCOS_0000_2000_OHCA_ZJ, for OHC in ZJ; the attribute “GCOS_area” is included for both variable types in the netcdf file and it tells the corresponding surface area); 2. Volume averaged temperature anomaly (global) timeseries (e.g., for level 0-2000m: GCOS_0000_2000_vol_ave_temp_anom, in degC; the attribute “GCOS_volume'' is included for this variable type in the netcdf file and it tells the corresponding volume). Regions of the ocean that are shallower than 300 m or are not sufficiently well sampled by the Argo array are not included.