The monthly global 2° × 2° Extended Reconstructed Sea Surface Temperature (ERSST) has been revised and updated from version 4 to version 5. This update incorporates a new release of ICOADS release 3.0 (R3.0), a decade of near-surface data from Argo floats, and a new estimate of centennial sea ice from HadISST2. A number of choices in aspects of quality control, bias adjustment, and interpolation have been substantively revised. The resulting ERSST estimates have more realistic spatiotemporal variations, better representation of high-latitude SSTs, and ship SST biases are now calculated relative to more accurate buoy measurements, while the global long-term trend remains about the same. Progressive experiments have been undertaken to highlight the effects of each change in data source and analysis technique upon the final product. The reconstructed SST is systematically decreased by 0.077°C, as the reference data source is switched from ship SST in ERSSTv4 to modern buoy SST in ERSSTv5. Furthermore, high-latitude SSTs are decreased by 0.1°–0.2°C by using sea ice concentration from HadISST2 over HadISST1. Changes arising from remaining innovations are mostly important at small space and time scales, primarily having an impact where and when input observations are sparse. Cross validations and verifications with independent modern observations show that the updates incorporated in ERSSTv5 have improved the representation of spatial variability over the global oceans, the magnitude of El Niño and La Niña events, and the decadal nature of SST changes over 1930s–40s when observation instruments changed rapidly. Both long- (1900–2015) and short-term (2000–15) SST trends in ERSSTv5 remain significant as in ERSSTv4.

The NOAA/NESDIS/NCEI Daily Optimum Interpolation Sea Surface Temperature (SST), version 2.0, dataset (DOISST v2.0) is a blend of in situ ship and buoy SSTs with satellite SSTs derived from the Advanced Very High Resolution Radiometer (AVHRR). DOISST v2.0 exhibited a cold bias in the Indian, South Pacific, and South Atlantic Oceans that is due to a lack of ingested drifting-buoy SSTs in the system, which resulted from a gradual data format change from the traditional alphanumeric codes (TAC) to the binary universal form for the representation of meteorological data (BUFR). The cold bias against Argo was about −0.14°C on global average and −0.28°C in the Indian Ocean from January 2016 to August 2019. We explored the reasons for these cold biases through six progressive experiments. These experiments showed that the cold biases can be effectively reduced by adjusting ship SSTs with available buoy SSTs, using the latest available ICOADS R3.0.2 derived from merging BUFR and TAC, as well as by including Argo observations above 5-m depth. The impact of using the satellite MetOp-B instead of NOAA-19 was notable for high-latitude oceans but small on global average, since their biases are adjusted using in situ SSTs. In addition, the warm SSTs in the Arctic were improved by applying a freezing point instead of regressed ice-SST proxy. This paper describes an upgraded version, DOISST v2.1, which addresses biases in v2.0. Overall, by updating v2.0 to v2.1, the biases are reduced to −0.07° and −0.14°C in the global ocean and Indian Ocean, respectively, when compared with independent Argo observations and are reduced to −0.04° and −0.08°C in the global ocean and Indian Ocean, respectively, when compared with dependent Argo observations. The difference against the Group for High Resolution SST (GHRSST) Multiproduct Ensemble (GMPE) product is reduced from −0.09° to −0.01°C in the global oceans and from −0.20° to −0.04°C in the Indian Ocean.

The uncertainty in Extended Reconstructed SST (ERSST) version 4 (v4) is reassessed based upon 1) reconstruction uncertainties and 2) an extended exploration of parametric uncertainties. The reconstruction uncertainty (Ur) results from using a truncated (130) set of empirical orthogonal teleconnection functions (EOTs), which yields an inevitable loss of information content, primarily at a local level. The Ur is assessed based upon 32 ensemble ERSST.v4 analyses with the spatially complete monthly Optimum Interpolation SST product. The parametric uncertainty (Up) results from using different parameter values in quality control, bias adjustments, and EOT definition etc. The Up is assessed using a 1000-member ensemble ERSST.v4 analysis with different combinations of plausible settings of 24 identified internal parameter values. At the scale of an individual grid box, the SST uncertainty varies between 0.3° and 0.7°C and arises from both Ur and Up. On the global scale, the SST uncertainty is substantially smaller (0.03°–0.14°C) and predominantly arises from Up. The SST uncertainties are greatest in periods and locales of data sparseness in the nineteenth century and relatively small after the 1950s. The global uncertainty estimates in ERSST.v4 are broadly consistent with independent estimates arising from the Hadley Centre SST dataset version 3 (HadSST3) and Centennial Observation-Based Estimates of SST version 2 (COBE-SST2). The uncertainty in the internal parameter values in quality control and bias adjustments can impact the SST trends in both the long-term (1901–2014) and “hiatus” (2000–14) periods.

This analysis estimates uncertainty in the NOAA global surface temperature (GST) version 5 (NOAAGlobalTemp v5) product, which consists of sea surface temperature (SST) from the Extended Reconstructed SST version 5 (ERSSTv5) and land surface air temperature (LSAT) from the Global Historical Climatology Network monthly version 4 (GHCNm v4). Total uncertainty in SST and LSAT consists of parametric and reconstruction uncertainties. The parametric uncertainty represents the dependence of SST/LSAT reconstructions on selecting 28 (6) internal parameters of SST (LSAT), and is estimated by a 1000-member ensemble from 1854 to 2016. The reconstruction uncertainty represents the residual error of using a limited number of 140 (65) modes for SST (LSAT). Uncertainty is quantified at the global scale as well as the local grid scale. Uncertainties in SST and LSAT at the local grid scale are larger in the earlier period (1880s–1910s) and during the two world wars due to sparse observations, then decrease in the modern period (1950s–2010s) due to increased data coverage. Uncertainties in SST and LSAT at the global scale are much smaller than those at the local grid scale due to error cancellations by averaging. Uncertainties are smaller in SST than in LSAT due to smaller SST variabilities. Comparisons show that GST and its uncertainty in NOAAGlobalTemp v5 are comparable to those in other internationally recognized GST products. The differences between NOAAGlobalTemp v5 and other GST products are within their uncertainties at the 95% confidence level.

The mixed layer heat budget in the tropical Pacific is diagnosed using pentad (5 day) averaged outputs from the Global Ocean Data Assimilation System (GODAS), which is operational at the National Centers for Environmental Prediction (NCEP). The GODAS is currently used by the NCEP Climate Prediction Center (CPC) to monitor and to understand El Niño and La Niña in near real time. The purpose of this study is to assess the feasibility of using an operational ocean data assimilation system to understand SST variability. The climatological mean and seasonal cycle of mixed layer heat budgets derived from GODAS agree reasonably well with previous observational and model-based estimates. However, significant differences and biases were noticed. Large biases were found in GODAS zonal and meridional currents, which contributed to biases in the annual cycle of zonal and meridional advective heat fluxes. The warming due to tropical instability waves in boreal fall is severely underestimated owing to use of a 4-week data assimilation window. On interannual time scales, the GODAS heat budget closure is good for weak-to-moderate El Niños. A composite for weak-to-moderate El Niños suggests that zonal and meridional temperature advection and vertical entrainment/diffusion all contributed to the onset of the event and that zonal advection played the dominant role during decay of the event and the transition to La Niña. The net surface heat flux acts as a damping during the development stage, but plays a critical role in the decay of El Niño and the transition to the following La Niña. The GODAS heat budget closure is generally poor for strong La Niñas. Despite the biases, the GODAS heat budget analysis tool is useful in monitoring and understanding the physical processes controlling SST variability associated with ENSO. Therefore, it has been implemented operationally at CPC in support of NOAA’s ENSO forecasting.