Global mean surface temperatures (GMST) exhibited a smaller rate of warming during 1998–2013, compared to the warming in the latter half of the 20th Century. Although, not a “true” hiatus in the strict definition of the word, this has been termed the “global warming hiatus” by IPCC (2013). There have been other periods that have also been defined as the “hiatus” depending on the analysis. There are a number of uncertainties and knowledge gaps regarding the “hiatus.” This report reviews these issues and also posits insights from a collective set of diverse information that helps us understand what we do and do not know. One salient insight is that the GMST phenomenon is a surface characteristic that does not represent a slowdown in warming of the climate system but rather is an energy redistribution within the oceans. Improved understanding of the ocean distribution and redistribution of heat will help better monitor Earth's energy budget and its consequences. A review of recent scientific publications on the “hiatus” shows the difficulty and complexities in pinpointing the oceanic sink of the “missing heat” from the atmosphere and the upper layer of the oceans, which defines the “hiatus.” Advances in “hiatus” research and outlooks (recommendations) are given in this report.

A self-organizing map (SOM) neural network was developed from Argo gridded datasets in order to estimate a subsurface temperature anomaly (STA) from remote sensing data. The SOM maps were trained using anomalies of sea surface temperature (SST), height (SSH), and salinity (SSS) data from Argo gridded monthly anomaly datasets, labeled with Argo STA data from 2005 through 2010, which were then used to estimate the STAs at different depths in the North Atlantic from the sea surface data. The estimated STA maps and time series were compared with Argo STAs including independent datasets for validation. In the Gulf Stream path areas, the STA estimations from the SOM algorithm show good agreement with in situ measurements taken from the surface down to 700-m depth, with a correlation coefficient larger than 0.8. Sensitivity of the SOM, when including salinity, shows that with SSS anomaly data in the SOM training process reveal the importance of SSS information, which can improve the estimation of STA in the subtropical ocean by up to 30%. In subpolar basins, the monthly climatology SST and SSH can also help to improve the estimation by as much as 40%. The STA time series for 1993–2004 in the midlatitude North Atlantic were estimated from remote sensing SST and altimetry time series using the SOM algorithm. Limitations for the SOM algorithm and possible error sources in the estimation are briefly discussed.

Retrieving multi-temporal and large-scale thermohaline structure information of the interior of the global ocean based on surface satellite observations is important for understanding the complex and multidimensional dynamic processes within the ocean. This study proposes a new ensemble learning algorithm, extreme gradient boosting (XGBoost), for retrieving subsurface thermohaline anomalies, including the subsurface temperature anomaly (STA) and the subsurface salinity anomaly (SSA), in the upper 2000 m of the global ocean. The model combines surface satellite observations and in situ Argo data for estimation, and uses root-mean-square error (RMSE), normalized root-mean-square error (NRMSE), and R² as accuracy evaluations. The results show that the proposed XGBoost model can easily retrieve subsurface thermohaline anomalies and outperforms the gradient boosting decision tree (GBDT) model. The XGBoost model had good performance with average R² values of 0.69 and 0.54, and average NRMSE values of 0.035 and 0.042, for STA and SSA estimations, respectively. The thermohaline anomaly patterns presented obvious seasonal variation signals in the upper layers (the upper 500 m); however, these signals became weaker as the depth increased. The model performance fluctuated, with the best performance in October (autumn) for both STA and SSA, and the lowest accuracy occurred in January (winter) for STA and April (spring) for SSA. The STA estimation error mainly occurred in the El Niño-Southern Oscillation (ENSO) region in the upper ocean and the boundary of the ocean basins in the deeper ocean; meanwhile, the SSA estimation error presented a relatively even distribution. The wind speed anomalies, including the u and v components, contributed more to the XGBoost model for both STA and SSA estimations than the other surface parameters; however, its importance at deeper layers decreased and the contributions of the other parameters increased. This study provides an effective remote sensing technique for subsurface thermohaline estimations and further promotes long-term remote sensing reconstructions of internal ocean parameters.