Cited by 150Climate variability and modelsOceanographic and Atmospheric ProcessesMarine Plankton and Biogeochemistry

The DOI matches its Crossref record; 74 of the 77 references it cites were found in a registry or an index; 36 of the 101 sources its claims cite were checked for support.

Tendencies, variability and persistence of sea surface temperature anomalies

Open access
Bulgin, Claire E., Merchant, Christopher J., Ferreira, David·2020·Scientific Reports

Abstract

Quantifying global trends and variability in sea surface temperature (SST) is of fundamental importance to understanding changes in the Earth’s climate. One approach to observing SST is via remote sensing. Here we use a 37-year gap-filled, daily-mean analysis of satellite SSTs to quantify SST trends, variability and persistence between 1981–2018. The global mean warming trend is 0.09 K per decade globally, with 95% of local trends being between −0.1 K and + 0.35 K. Excluding perennial sea-ice regions, the mean warming trend is 0.11 K per decade. After removing the long-term trend we calculate the SST power spectra over different time periods. The maximum variance in the SST power spectra in the equatorial Pacific is 1.9 K 2 on 1–5 year timescales, dominated by ENSO processes. In western boundary currents characterised by an intense mesoscale activity, SST power on sub-annual timescales dominates, with a maximum variance of 4.9 K 2. Persistence timescales tend to be shorter in the summer hemisphere due to the shallower mixed layer. The median short-term persistence length is 11–14 days, found over 71–79% of the global ocean area, with seasonal variations. The mean global correlation between monthly SST anomalies with a three-month time-lag is 0.35, with statistically significant correlations over 54.0% of the global oceans, and notably in the northern and equatorial Pacific, and the sub-polar gyre south of Greenland. At six months, the mean global SST anomaly correlation falls to 0.18. The satellite data record enables the detailed characterisation of temporal changes in SST over almost four decades.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

Using a 37-year gap-filled daily satellite SST analysis (1981–2018), the paper quantifies global SST trends, variability, and persistence. It finds a global mean warming trend of 0.09 K per decade, ENSO-dominated 1–5 year variability in the equatorial Pacific, and short-term anomaly persistence typically 11–14 days over most of the ocean.

Key findings

  • Global mean SST warming trend is 0.09 ± 0.006 K per decade including perennial sea-ice regions and 0.11 ± 0.008 K per decade excluding them; 95% of local trends lie between −0.1 K and +0.35 K.
  • After detrending, maximum SST power-spectrum variance in the equatorial Pacific is 1.9 K² on 1–5 year timescales, dominated by ENSO; in western boundary currents, sub-annual power reaches a maximum variance of 4.9 K² in the Brazilian-Falkland confluence.
  • Median short-term SSTA persistence length is 11–14 days, found over 71–79% of the global ocean area, with seasonal variations.
  • Mean global SSTA correlation at a three-month lag is 0.35, with statistically significant correlations over 54.0% of the global oceans; at six months the mean global correlation falls to 0.18.

Methods

  • Uses the ESA CCI SST analysis: gap-filled, daily-mean satellite SST from Sept 1981–Dec 2018, originally 0.05° grid and spatially averaged to 0.25°, combining AVHRR and ATSR-series retrievals.
  • Trends are calculated with the Theil-Sen estimator; variability is assessed from decad (10-day) detrended SSTA power spectra over sub-annual, 1–5 year, and >5 year timescales.
  • Short-term persistence is derived from daily detrended SSTA autocorrelations, fitted with a Savitzky-Golay filter and exponential decay, using a correlation threshold of 0.6; longer persistence/re-emergence is assessed via Pearson correlations at 3-, 6-, and 12-month lags.

Limitations

  • Stated: the 37-year record is relatively short for estimating >5-year SST power; 12-month lag correlations use 36 years; persistence estimates are sensitive to temporal resolution, persistence definition, and fit window; re-emergence can skew exponential fits; high spatial resolution may make co-located re-emergence harder to identify.
  • Stated: uncertainty values exclude instrument instability over time, estimated at about 0.05 K per decade; point-by-point significance over a spatially correlated field means some local significant areas are expected by chance; dust events can bias retrieved SST cold and create artificial short-term persistence in the equatorial Atlantic.
  • Evident from the text: the analysis is satellite-era only and depends on the ESA CCI product’s retrieval and gap-filling characteristics, with limited ability to resolve long-term or re-emergence signals in noisy or advective regions.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s41598-020-64785-9

Cites

74 of 75 references verified0 citation errors · 1 not checkable

Statistics

Nothing to recompute.

Figures

Figure 1 (a) Global sea-surface temperature anomaly (SSTA) trends between 1981–2018 relative to a climatology between 1982–2010 at 0.25 degree resolution. (b) Timeseries of SSTA between 1981–2018 for the global data (raw data in grey, smoothed data in black, multiplied by a factor of 10), the Atlantic Meridional Variability Region [0–65 N, 75–7.5 W30 (blue, multiplied by a factor of 3), the Nino 3.4 region [5N-5S, 170–120 W31 (red) and the NE Pacific warm waters [32–54 N, 165–138 W32 (cyan). Decadal tendencies (degrees kelvin per decade) calculated from data for (c) 1981–1989, (d) 1990–1999, (e) 2000–2009 and (f) 2010–2018 are at one degree resolution. Note that the colour scale is different (+/−4 K) for decadal panels (c–f) than for the full-period figure (a,+/−0.8 K).
Figure 1 (a) Global sea-surface temperature anomaly (SSTA) trends between 1981–2018 relative to a climatology between 1982–2010 at 0.25 degree resolution. (b) Timeseries of SSTA between 1981–2018 for the global data (raw data in grey, smoothed data in black, multiplied by a factor of 10), the Atlantic Meridional Variability Region [0–65 N, 75–7.5 W30 (blue, multiplied by a factor of 3), the Nino 3.4 region [5N-5S, 170–120 W31 (red) and the NE Pacific warm waters [32–54 N, 165–138 W32 (cyan). Decadal tendencies (degrees kelvin per decade) calculated from data for (c) 1981–1989, (d) 1990–1999, (e) 2000–2009 and (f) 2010–2018 are at one degree resolution. Note that the colour scale is different (+/−4 K) for decadal panels (c–f) than for the full-period figure (a,+/−0.8 K).CC BY · PubMed Central
Figure 2 (a) SST variability in the monthly climatology from 1982–2010. (b–d) SST integrated amplitude attributable to different parts of the SST power spectrum: (b) sub-annual timescales, (c) timescales of 1–5 years, (d) timescales longer than 5 years.
Figure 2 (a) SST variability in the monthly climatology from 1982–2010. (b–d) SST integrated amplitude attributable to different parts of the SST power spectrum: (b) sub-annual timescales, (c) timescales of 1–5 years, (d) timescales longer than 5 years.CC BY · PubMed Central
Figure 3 Seasonally-dependent short-term persistence of SST anomalies starting in different months ((a) December, (b) March, (c) June and (d) September). Persistence is calculated by correlating SST anomalies from the 1st, 15th and final day of each month with the SST anomalies on subsequent days, over a 30-day period across all years. These daily correlations are fitted using a 15-day Savitzky-Golay filter with the time series truncated where the smoothed data first drop below a correlation of 0.6. An exponential curve fitted to the remaining time series is used to define short-term persistence, using a correlation threshold of 0.658. Ocean regions with no short-term SSTA persistence calculable by this method are coloured white. Sea-ice pixels are shown in light grey, and land in dark grey.
Figure 3 Seasonally-dependent short-term persistence of SST anomalies starting in different months ((a) December, (b) March, (c) June and (d) September). Persistence is calculated by correlating SST anomalies from the 1st, 15th and final day of each month with the SST anomalies on subsequent days, over a 30-day period across all years. These daily correlations are fitted using a 15-day Savitzky-Golay filter with the time series truncated where the smoothed data first drop below a correlation of 0.6. An exponential curve fitted to the remaining time series is used to define short-term persistence, using a correlation threshold of 0.658. Ocean regions with no short-term SSTA persistence calculable by this method are coloured white. Sea-ice pixels are shown in light grey, and land in dark grey.CC BY · PubMed Central
Figure 4 SSTA correlations with a three-month time-lag from (a) December, (b) March, (c) June and (d) September. Correlations are calculated for 36 complete years of data (1982–2017). Sea-ice pixels are shown in light grey and land in dark grey. Statistically significant correlations are indicated by the overlaid hatching.
Figure 4 SSTA correlations with a three-month time-lag from (a) December, (b) March, (c) June and (d) September. Correlations are calculated for 36 complete years of data (1982–2017). Sea-ice pixels are shown in light grey and land in dark grey. Statistically significant correlations are indicated by the overlaid hatching.CC BY · PubMed Central
Figure 5 SSTA correlations with a six-month time-lag from (a) December, (b) March, (c) June and (d) September. Correlations are calculated for 36 complete years of data (1982–2017). Sea-ice affected pixels are shown in light grey and land in dark grey. Statistically significant correlations are indicated by the overlaid hatching.
Figure 5 SSTA correlations with a six-month time-lag from (a) December, (b) March, (c) June and (d) September. Correlations are calculated for 36 complete years of data (1982–2017). Sea-ice affected pixels are shown in light grey and land in dark grey. Statistically significant correlations are indicated by the overlaid hatching.CC BY · PubMed Central
Figure 6 SSTA correlations with a twelve-month time-lag from (a) December, (b) March, (c) June and (d) September. Correlations are calculated for 36 complete years of data (1982–2017). Sea-ice affected pixels are shown in light grey and land in dark grey. Statistically significant correlations are indicated by the overlaid hatching.
Figure 6 SSTA correlations with a twelve-month time-lag from (a) December, (b) March, (c) June and (d) September. Correlations are calculated for 36 complete years of data (1982–2017). Sea-ice affected pixels are shown in light grey and land in dark grey. Statistically significant correlations are indicated by the overlaid hatching.CC BY · PubMed Central
Figure 7 (a) Uncertainty in the Theil-Sen slope fit82 of the sea-surface temperature anomaly (SSTA) trend between 1981–2018 at 0.25 degree resolution. (b) Timeseries of SSTA uncertainties between 1981–2018 for global data (black), the Atlantic Meridional Variability Region [0–65 N, 75–7.5 W30 (blue), the Nino 3.4 region [5N-5S, 170–120 W31 (red) and the NE Pacific warm waters [32–54 N, 165–138 W32 (cyan). Uncertainties in decadal tendencies (degrees kelvin per decade) calculated from data for (c) 1981–1989, (d) 1990–1999, (e) 2000–2009 and (f) 2010–2018 are at one degree resolution. Note that the colour scale is different (0–1.5 K) for decadal panels (c–f) than for the full-period figure (a, 0–0.4 K).
Figure 7 (a) Uncertainty in the Theil-Sen slope fit82 of the sea-surface temperature anomaly (SSTA) trend between 1981–2018 at 0.25 degree resolution. (b) Timeseries of SSTA uncertainties between 1981–2018 for global data (black), the Atlantic Meridional Variability Region [0–65 N, 75–7.5 W30 (blue), the Nino 3.4 region [5N-5S, 170–120 W31 (red) and the NE Pacific warm waters [32–54 N, 165–138 W32 (cyan). Uncertainties in decadal tendencies (degrees kelvin per decade) calculated from data for (c) 1981–1989, (d) 1990–1999, (e) 2000–2009 and (f) 2010–2018 are at one degree resolution. Note that the colour scale is different (0–1.5 K) for decadal panels (c–f) than for the full-period figure (a, 0–0.4 K).CC BY · PubMed Central

Verification

Claims the paper makes

69 Own result · 21 No full text of the sources · 21 Own interpretation · 19 Partly supported · 14 No basis given · 7 Sources not found · 5 Supported · 4 Supported so far · 4 Own method · 2 Couldn’t confirm · accuracy not yet measured

Show all 105 claims

  • It may state two separate claims as one. 49
  • Probably not a substantial claim. 13
  • The quote only partly states this claim. 1
61 of 61
  • The authors state this as an assumption.

    The daily data are fully correlated over the spatial scale when regridding from 0.05° to 0.25°.

    “assuming that these data are fully correlated over this spatial scale”
  • The authors state this as an assumption.

    The temporal correlation length scale is 3 days and the spatial correlation is 100 km for uncertainty propagation.

    “assuming a temporal correlation length scale of 3 days and a spatial correlation of 100 km”
2 of 2

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

2 of 2

Citations

Cites

Show all 77 references

Research works printed but not linked: 6

  1. Hartmann, D. L., Albert M. G, Klein T. & Rusticucci, M. Intergovernmental Panel on Climate Change, Climate Change 2013, Physical Science Basis, Chapter 2: Observations: Atmosphere and Surface, https://www.ipcc.ch/report/ar5/wg1/ (2013).
  2. Casey, K. S. & Cornillon, P. A Comparision of Satellite and In Situ-Based Sea Surface Temperature Climatologies. Journal of Climate 12, 1848–1863 (1998). · DOI
  3. Hanley, D. E., Bourassa, M. A., O’Brien, J. J., Smith, S. R. & Spade, E. R. A Quantitative Evaluation of ENSO Indices. Journal of Climate 16, 1249–1258 (2003). · DOI
  4. Rayner, N. A. & Kennedy, J. J. Comparison of satellite and in situ measured estimates of sea surface temperature focusing on the ESA Climate Change Initiative v2.1 products. Earth and Space Science (In review, 2020).
  5. Reynolds, R. W., Rayner, N. A., Smith, T. M., Stokes, D. C. & Wang, W. An Improved In Situ and Satellite SST Analysis for Climate. Journal of Climate 15, 1610–1625 (2002). · DOI
  6. Armour, K. C. & Bitz, C. M. Observed and projected trends in Antarctic sea ice. US CLIVAR Variations 13, 12–19 (2015).

Cited by

Sorted by citations

Metadata

Authors
Bulgin, C. E., Merchant, C. J., & Ferreira, D.
Venue
Scientific Reports
Published
2020-05-14
Volume
10
Issue
1
Pages
7986-7986
ISSN
2045-2322 (ISSN-L)
OpenAlex
W3022238508
Cited by
150
References
77
Citations by year
2026: 15 · 2025: 24 · 2024: 39 · 2023: 21 · 2022: 27 · 2021: 24
FWCI
6.55
Type
journal article
Language
en
Open access
Locations
  • Scientific Reports (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • CentAUR (University of Reading) (acceptedVersion)
  • PubMed Central (submittedVersion, cc-by)
Funders
  • European Space Agency
  • Sight Research UK (NE/R016518/1, nceo020006)
  • Natural Environment Research Council (nceo020006, NE/R016518/1)
Topics
  • Physical Sciences › Environmental Science › Global and Planetary Change › Climate variability and models
  • Physical Sciences › Earth and Planetary Sciences › Oceanography › Oceanographic and Atmospheric Processes
  • Physical Sciences › Earth and Planetary Sciences › Oceanography › Marine Plankton and Biogeochemistry
Keywords
sea surface temperature anomalies, ocean variability, power spectrum, ENSO, western boundary currents, mesoscale variability, mixed layer depth, seasonal variation, global ocean warming

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-16, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-16 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review