Back to report · Publisher record

Hawkins et al. (2020). Observed Emergence of the Climate Change Signal: From the Familiar to the Unknown

Page 1

Figures and page graphics, page 1

 

 

Observed emergence of the climate change signal: from the familiar to the unknown 

E. Hawkins

1

, D. Frame

2

, L. Harrington

3

, M. Joshi

4

, A. King

5

, M. Rojas

6

, and R. Sutton

1

  

1

 National Centre for Atmospheric Science, Department of Meteorology, University of Reading. 

2

 New Zealand Climate Change Research Institute, Victoria University of Wellington, 

Wellington 6012, New Zealand. 

3

 

4

 Climatic Research Unit, School of Environmental Sciences, University of East Anglia, UK 

5

 School of Earth Sciences and ARC Centre of Excellence for Climate Extremes, University of 

Melbourne. 

6

 Department of Geophysics and Centre for Climate and Resilience Research, CR2, University of 

Chile 

Corresponding author: Ed Hawkins (e.hawkins@reading.ac.uk) 

Key Points: 

 

The signal of changes in observed temperature and rainfall due to global warming has 
clearly emerged in many regions and at meso-scales 

 

Tropical regions have experienced the largest changes in temperature relative to the 
amplitude of internal variability 

 

Signals of increasing extreme rainfall are emerging more quickly than signals in mean 
rainfall over many parts of the UK 

 

 

 

Page 2

Figures and page graphics, page 2

 

 

Abstract 

Changes in climate are usually considered in terms of trends or differences over time. However, 
for many impacts requiring adaptation, it is the amplitude of the change relative to the local 
amplitude of climate variability which is more relevant. Here, we explore the concept of ‘signal-
to-noise’ in observations of local temperature, highlighting that many regions are already 
experiencing a climate which would be ‘unknown’ by late 19

th

 century standards. The emergence 

of temperature changes over both land and ocean is clearest in tropical regions, in contrast to the 
regions of largest change which are in the northern extra-tropics – broadly consistent with 
climate model simulations. Significant increases and decreases in rainfall have also already 
emerged in different regions with the UK experiencing a shift towards more extreme rainfall 
events, a signal which is emerging more clearly in some places than the changes in mean rainfall. 

Plain Language Summary 

Changes in climate are translated into impacts on society not just though the amount of change, 
but how this change compares to the variations in climate that society is used to. Here we 
demonstrate that significant changes, when compared to the size of past variations, are present in 
both temperature and rainfall observations over many parts of the world. 

1 Introduction 

It was first noted that surface air temperatures were increasing at both local and global scales 
more than 80 years ago [

Kincer

 1933, 

Callendar 

1938]. At the time it was unclear whether the 

observed changes were part of a longer term trend or a natural fluctuation – the ‘signal’ had not 
yet clearly emerged from the ‘noise’ of variability – although 

Callendar

 [1938] did suggest that 

the increase in atmospheric carbon dioxide concentrations was partly to blame. 

The concept of the emergence of a climate change signal has since been discussed extensively, 
often linked with the detection & attribution of climatic changes. For example, 

Madden & 

Ramanathan

 [1980] and 

Wigley & Jones

 [1981] could not robustly detect the carbon dioxide 

warming signal, but 

Hansen et al. 

[1988] predicted that the ratio of temperature change and the 

magnitude of interannual variability – the signal-to-noise ratio – would be above 3 in large parts 
of the tropics by the 2010s, with smaller values over high latitude land regions. 

Mahlstein et al.

 

[2011, 2012] subsequently demonstrated that the signal had indeed emerged in the observations, 
especially in the tropics in boreal summer, and with a similar pattern to that expected from 
climate model simulations. 

Lehner et al. 

[2017] subsequently highlighted emergence of observed 

temperature changes in both winter and summer in the northern extra-tropics. Significant 
changes in precipitation are often harder to detect because both thermodynamic and dynamic 
factors are crucial [e.g. 

Zappa & Shepherd

, 2017] and because internal variability in precipitation 

is larger. However, precipitation changes are apparent in some regions [e.g. 

Zhang et al.

 2007] 

including in extremes [e.g. 

Min et al.

 2011]. 

Many studies have also considered when further changes in climate will emerge, for both mean 
temperature [

Mahlstein et al.

 2011, 

Hawkins & Sutton

 2012] and precipitation [

Giorgi & Bi

 

2009, 

Fischer et al.

 2014]. Other studies have considered when changes in climate extremes 

should have emerged in the past [

King et al.

 2015] or future [

Diffenbaugh & Scherer

 2011, 

Page 3

Figures and page graphics, page 3

 

 

Fischer et al.

 2014]. However, rather than examine the timing of any climate emergence, we 

focus here on the related quantity – signal-to-noise.  

The clearest emergence of warming – and largest signal-to-noise values – tend to be found in the 
tropics, which are regions with large and vulnerable populations [

Frame et al.

 2017, 

Harrington 

et al. 

2017]. Signal-to-noise (S/N) is important for climate impacts, especially for ecosystems 

which have a limited ability to adapt and so large changes outside past experience could be 
particularly harmful [

Deutsch et al. 

2008; 

Beaumont et al.

 2011]. Crop growing areas also face 

unprecedented heat [

Battisti & Naylor

 2009] and changes in rainfall which may move outside 

past experiences [

Rojas et al.

 2019]. The impacts of shifts in snowfall [

Diffenbaugh et al.

 2012] 

and Köppen–Geiger zones [

Mahlstein et al.

 2013] have also been discussed in terms related to 

the natural variability of the local conditions. Quantifying the changes that have already occurred 
may help determine which regions are suffering the largest adverse consequences of a warming 
world. 

Here, we revisit the question of where and how the signal of climatic changes is emerging from 
the background noise of internal variability. In contrast to most previous studies we focus our 
analysis on observational datasets of temperature and precipitation, with model simulations used 
only to test the methodology. 

2. 

Observed emergence and signal-to-noise 

2.1 

Methodology 

Our aim is to produce estimates of signal-to-noise (S/N) for changes in observed climate 
variables without utilising data from any climate model simulations. The simple approach 
adopted is to linearly regress local variations in climate onto annual global mean surface 
temperature change (GMST), i.e. 

L

(t) = 

G

(t) + 

, 

where 

L

(t) is the local change (in temperature or precipitation) over time, 

G

(t) is a smoothed 

version of GMST change over the same period, 

 defines the linear scaling between 

L

 and 

G

, 

and 

 is a constant. 

Sutton et al.

 [2015] highlighted that a large fraction of variance in local 

climate changes can be represented by GMST changes, and 

Fischer et al.

 [2014] demonstrated 

that a similar regression approach provided robust estimates of S/N when examining future 
changes in precipitation in climate model simulations. 

For 

G

(t) we use GMST from the Berkeley Earth temperature dataset for 1850-2018 (

Rohde et al.

 

[2013], combined with HadSST3 from 

Kennedy et al.

 [2011]), relative to the mean of 1850-

1900, and smoothed with a lowess filter of 41-years to highlight the long-term variations (Figure 
1a). The conclusions are insensitive to whether the smoothing parameter is slightly larger or 
smaller. The ‘signal’ of global temperature change is defined as the value of the smoothed 

Page 4

Figures and page graphics, page 4

 

 

GMST in 2018 (

G

2018

 = 1.19K), the ‘signal’ of local climate change explained by GMST is 

G

 

and the ‘noise’ is defined as the standard deviation of the residuals (

L

 – 

G

).  

Although we do not formally attribute the observed change in GMST, and hence local changes, 
to particular radiative forcings or feedbacks, applying the method of 

Haustein et al.

 [2017] to 

derive a GMST change that is attributable to human activity gives 1.22K, similar to 

G

2018

. 

Although 1850-1900 is often considered as a proxy for ‘pre-industrial’ GMST, the 

Haustein et 

al. 

[2017] approach also suggests an additional anthropogenic warming of around 0.05K 

occurred between 1750 and 1850-1900, based on radiative forcing estimates back to 1750. 
Although this plausible pre-1850 attributable warming is not included in our analysis, we refer to 
the 1850-1900 period as the early-industrial era, rather than pre-industrial. 

2.2 

Example for annual mean temperatures in Oxford 

To demonstrate our approach we consider a case study of temperature change in Oxford, UK. 

Burt & Burt

 [2019] produced an extended temperature record for the Oxford Radcliffe 

Observatory with annual means available for 1814-2018. The temporal evolution of GMST and 
temperatures in Oxford are similar, showing that the ‘fingerprint’ of GMST change is clearly 
visible at the spatial scale of a single continuous weather station, although with more noise at the 
local scale (Figure 1b, also see 

Sutton et al.

 [2015]). We note that there is likely an urban heat 

island influence on temperatures in Oxford of around 0.1-0.2K [

Burt & Burt

 2019].  

We regress this local temperature dataset onto smoothed GMST and obtain 

 = 1.45 ± 0.25 (95% 

confidence interval). The ‘signal’ for Oxford is 

G

2018

 = 1.72 ± 0.30K and the ‘noise’, i.e. the 

local variations that are not explained by GMST variations, is 0.54K. Oxford therefore exhibits a 
S/N ratio of 3.2 ± 0.5 (Figure 1b). 

We adopt the language of 

Frame et al.

 [2017] to describe how the climate has changed from 

being familiar, to being ‘

unusual

’ relative to lived experience (S/N > 2), ‘

unknown

’ (S/N > 3), 

and here we introduce ‘

inconceivable

’ for S/N values above 5 (Fig. S1). Using this terminology, 

temperatures in Oxford have become unknown relative to the early-industrial era. Two other 
regional examples are illustrated in Fig. S2. 

2.3 

Local climate data and methodological tests 

We perform a similar S/N analysis for each land and ocean gridpoint in the Berkeley Earth 
temperature dataset (1850-2018) and in the GPCCv2018 land precipitation dataset (1891-2016, 

Schneider et al.

 [2017]). We use the 1º x 1º datasets for both Berkeley Earth and GPCC. We also 

use the HadUK-Grid dataset for the UK [

Hollis et al.

 2019] at 25km spatial resolution for 

monthly (1862-2017) and daily (1891-2017) precipitation data to examine changes in mean 
rainfall and extremes. Note that smoothed GMST (1850-2018) is used as 

G

 for both local 

temperature and precipitation analyses. 

As the local data is not necessarily available for all years back to 1850 we perform the regression 
only over the period where local temperatures or precipitation are defined. The signal relative to 
the early-industrial era can still be calculated assuming that the estimated regression parameter 
(

), is representative for the whole period, i.e. the signal is always 

G

2018

, irrespective of the 

Page 5

Figures and page graphics, page 5

 

 

time period used to calculate 

. However, we require that there must be at least 100 years of 

local climate data available. 

We test our methodology using a large ensemble of climate simulations for the historical period 
[

Maher et al.

 2019], specifically to examine the uncertainty due to internal variability in derived 

S/N values for temperature and precipitation. Figs. S3 and S4 demonstrate that the methodology 
produces S/N values with small uncertainties (typically <0.4 over land regions) and robust 
patterns. 

3. 

Emergence of unknown temperatures 

The map of the current observed signal of annual temperature change, relative to the early-
industrial era, is shown in Figure 2a. It shows the familiar pattern of more warming over land 
than over the oceans, more warming at high northern latitudes, and less warming in the tropical 
regions and the southern hemisphere. Virtually all locations have experienced more than 1K 
change since the early-industrial era, and many regions have exceeded 2K. The estimated noise 
shows a similar pattern with larger variability at higher northern latitudes, but the differences 
between the tropics and extra-tropics are more pronounced than for the signal (Figure 2b). 

The ratio of these two patterns results in a signal-to-noise (S/N) map with the largest values in 
the tropical regions (Figure 2c). Although these areas generally have smaller signals than higher 
latitude regions, they have experienced a larger amplitude change relative to the (smaller) 
background variations in temperature than other regions. This is important as societies, 
infrastructure and ecosystems are often adapted for the range of local climate experienced. S/N 
measures how far the climate is being shifted from that past range; the climate in large parts of 
the tropics has shifted such that the mean climate would have been inconceivable in the early-
industrial era. More than half of the land area has experienced S/N above 3, and so has moved 
into a climate that is unknown by early-industrial standards (Fig. S5). 

Over the oceans the largest S/N values are found in the tropical Atlantic and tropical Indian 
Oceans. Fish species such as tuna have already been seen to be moving away from the tropics to 
the sub-tropics, likely to avoid these warmer waters [

Monllor-Hurtado et al. 

2017]. Large parts 

of the North Atlantic have seen little warming overall, likely due to changes in ocean circulation 
providing a local cooling influence to offset global warming [e.g. 

Dima & Lohmann

 2010].  

Although there are variations in magnitude, the estimated S/N pattern is relatively robust to the 
choice of temperature dataset [

Morice et al.

 2012, 

Cowtan & Way

 2014, 

Lenssen et al. 

2019, 

Zhang et al.

 2019]. However, there are notable local differences between datasets over south-east 

USA and parts of South America (Fig. S6). The overall observed emergence pattern is broadly 
similar to that found in models under future climate change scenarios [

Frame et al.

, 2017] 

though there are regional-scale differences; especially in the oceans but over some land areas 
too. 

When considering how changes in climate may be experienced, it may in many cases be more 
relevant to examine seasonal or monthly timescales, depending on the impact being considered. 
For example, Figure 3 shows that S/N values can still be significant for monthly average 
temperatures. Again, the largest S/N values are found in the tropics and tend to be larger for the 

Page 6

Figures and page graphics, page 6

 

 

climatologically warmest month than the climatologically coldest month for each location. This 
is because weather variability tends to be larger in the colder months. Around 40% of land areas 
have moved into an unusual climate in their warmest months, and 20% in the coldest months 
(Fig. S5). This suggests a comparatively large increase in likelihood of heat-related extreme 
events in already warm months of already hot countries. One example is south-east Asia where 
the S/N values are large and the combined effects of El Nino events and climate change on 
extreme heat in the warmest months of the year has previously been noted [

Thirumalai et al. 

2017]. 

4. 

Emergence of unusual precipitation amounts 

The S/N analysis is repeated for annual mean precipitation using the GPCC dataset. In this case, 
some regions are getting significantly wetter and others are getting significantly drier (Figure 4) 
but, unsurprisingly, the signals are less clear than for temperature. Notable emergence of 
‘

unfamiliar

’ (S/N > 1) or unusual precipitation changes are observed in west Africa, Brazil, Chile 

and south-west Australia (drier), and the northern high latitudes and Argentina (wetter). The 
seasonal values of S/N are shown in Fig. S7. The changes in several of these regions have been 
discussed as being consistent with the expected response to increased greenhouse gas forcing, 
e.g. for south-west Australia [

Delworth & Zeng

 2014], for Chile [

Boisier et al.

 2016] and the 

northern extra-tropics [

Zhang et al.

 2007]. 

To demonstrate that this framework can be applied to a range of gridded datasets and spatial 
scales, we consider one small region in more detail. The UK has a gridded rainfall dataset 
available, covering 1891-2017 (daily) and 1862-2017 (monthly), which is suitable for examining 
changes in mean and extreme rainfall [

Hollis et al.

 2019].  

Figure 5 shows the signal and S/N for annual mean rainfall, highlighting a tendency for 
increasing rainfall in large parts of the northern UK and the western coasts of up to 20% per K of 
GMST change. The corresponding S/N values exceed 1 in several areas, and these tend to be 
mountainous regions. Fig. S8 shows the seasonal mean S/N values. 

When considering the wettest day of the year (RX1day) as 

L

(t), there is a clear signal of 

increasing extreme rainfall, but the pattern is strikingly different to the mean. This signal is 
visible across large parts of the UK, even in regions where there are only small changes in mean 
rainfall. The signal has only clearly emerged in a few locations (Fig. 5) but the spatial average of 
RX1day across the UK suggests an increase in extreme rainfall amounts of around 4mm (or 
11%) per K of GMST change (Fig. S9), which is around 8% per K of UK temperature change, 
approximately consistent with Clausius-Clapeyron expectations [

Pall et al.

 2007].  

These findings are consistent with 

Min et al.

 [2011] who showed that the signal of changes in 

extreme rainfall were detectable and attributable to human activity over large parts of the 
northern hemisphere land areas, and with 

Fischer et al.

 [2014] used climate model simulations to 

suggest that emergence of changes in extreme rainfall can occur earlier than changes in mean 
rainfall. Continued recovery of millions of undigitized weather observations, including for daily 

Page 7

Figures and page graphics, page 7

 

 

rainfall, will improve and lengthen these gridded datasets [e.g. 

Ashcroft et al.

 2018; 

Hawkins et 

al.

 2019]. 

5. 

Summary and discussion 

We have estimated the signal-to-noise ratio (S/N) of observed temperature and precipitation 
changes since the early-industrial era (1850-1900). Although we do not formally attribute these 
local changes to specific radiative forcings or feedbacks, the emergence of significantly different 
climates is related to increases in GMST, which itself is largely due to anthropogenic factors 
[e.g. 

IPCC

 2018]. 

Consistent with previous studies and expectations from climate model simulations, the largest 
S/N values for historical temperature changes are seen in the tropical regions, over both land and 
ocean. Large regions have already experienced a shift to a climate state that is unknown, and 
even inconceivable, compared to that in the late 19th century. These signals of change are also 
clear in monthly average temperatures, with warmer months showing more significant changes. 

Precipitation signals are emerging in several regions when considering observed rainfall changes, 
particularly West Africa, parts of South America and northern Eurasia. Some regions in South 
America and central Africa exhibit simultaneously high S/N for temperature (S/N>4) and 
significantly drier precipitation (S/N<-1) which may compound impacts. 

As a demonstration of the methods in a data-rich region, and over a range of spatial scales, our 
analysis shows there are clear shifts towards more annual rainfall over the UK, focussed over 
northern and western areas. Significant increases in extreme heavy rainfall are emerging over 
large parts of the UK and are emerging more quickly than changes in mean rainfall in some 
places. The magnitude of the increase in extreme rainfall (~8% per K of local temperature 
change) is approximately consistent with expectations from the Clausius-Clapeyron relationship. 

Many of the largest global shifts in climate, relative to the background variability, are found in 
countries with large, vulnerable populations, and this will be exacerbated if policy targets such as 
those in the Paris Agreement are not met [

Frame et al. 

2017, 

King & Harrington

 2018]. There 

are also implications for ecosystems in these regions, which may not be able to adapt to such an 
unknown climate, especially given the rates of change. The rates of change of signal-to-noise to 
which societies and ecosystems can adapt is an important topic for future analyses. 

Acknowledgments 

EH and RS were supported by the National Centre for Atmospheric Science, and EH was funded 
by the INDECIS project. ADK was funded by the Australian Research Council (DE180100638). 
MJ is funded by UK Natural Environment Research Council grants NE/S004645/1 and 
NE/N018486/1. MR acknowledges funding from FONDAP-CONICYT 15110009. 

References 

Ashcroft, L., J. R. Coll, A. Gilabert, P. Domonkos, M. Brunet, E. Aguilar, M. Castella, J. Sigro, 

I. Harris, P. Unden and P. Jones, A rescued dataset of sub-daily meteorological 

Page 8

Figures and page graphics, page 8

 

 

observations for Europe and the southern Mediterranean region, 1877–2012, Earth Syst. 
Sci. Data, 10, 1613–1635doi: 10.5194/essd-10-1613-2018, 2018 

Battisti, D. S., and R. L. Naylor, Historical warnings of future food insecurity with 

unprecedented seasonal heat, Science, 323(5911), 240–244, doi:10. 
1126/science.1164363, 2009. 

Beaumont, L. J., A. Pitman, S. Perkins, N. E. Zimmermann, N. G. Yoccoz, and W. Thuiller, 

Impacts of climate change on the worlds most exceptional ecoregions, Proceedings of the 
National Academy of Sciences, 108 (6), 2306– 2311, doi:10.1073/pnas.1007217108, 
2011. 

Boisier, J. P., R. Rondanelli, R. D. Garreaud, and F. Munoz, Anthropogenic and natural 

contributions to the southeast pacific precipitation decline and recent megadrought in 
central Chile, Geophysical Research Letters, 43(1), 413–421, doi:10.1002/2015gl067265, 
2016. 

Burt, S. and Burt, T., Oxford Weather and Climate since 1767, Oxford University Press, 2019 

Callendar, G. S., The artificial production of carbon dioxide and its influence on temperature, 

Quarterly Journal of the Royal Meteorological Society, 64(275), 223–240, 
doi:10.1002/qj.49706427503, 1938. 

Cowtan, K., and R. G. Way, Coverage bias in the HadCRUT4 temperature series and its impact 

on recent temperature trends, Quarterly Journal of the Royal Meteorological Society, 
140(683), 1935–1944, doi: 10.1002/qj.2297, 2014. 

Delworth, T. L., and F. Zeng, Regional rainfall decline in Australia attributed to anthropogenic 

greenhouse gases and ozone levels, Nature Geoscience, 7(8), 583–587, 
doi:10.1038/ngeo2201, 2014. 

Deutsch, C. A., J. J. Tewksbury, R. B. Huey, K. S. Sheldon, C. K. Ghalambor, D. C. Haak, and 

P. R. Martin, Impacts of climate warming on terrestrial ectotherms across latitude, 
Proceedings of the National Academy of Sciences, 105(18), 6668–6672, 
doi:10.1073/pnas.0709472105, 2008. 

Diffenbaugh, N. S., and M. Scherer, Observational and model evidence of global emergence of 

permanent, unprecedented heat in the 20th and 21st centuries, Climatic Change, 107(3-4), 
615–624, doi:10.1007/ s10584-011-0112-y, 2011. 

Diffenbaugh, N. S., M. Scherer and M. Ashfaq, Response of snow-dependent hydrologic 

extremes to continued global warming, Nature Climate Change, 3, 379, doi: 
10.1038/nclimate1732, 2013 

Dima, M. and G. Lohmann, Evidence for Two Distinct Modes of Large-Scale Ocean Circulation 

Changes over the Last Century, Journal of Climate, 23, 5, doi: 10.1175/2009JCLI2867.1, 
2010 

Fischer, E. M., J. Sedlacek, E. Hawkins, and R. Knutti, Models agree on forced response pattern 

of precipitation and temperature extremes, Geophysical Research Letters, 41(23), 8554–
8562, doi:10.1002/2014gl062018, 2014. 

Page 9

Figures and page graphics, page 9

 

 

Frame, D., M. Joshi, E. Hawkins, L. J. Harrington, and M. de Roiste, Population-based 

emergence of unfamiliar climates, Nature Climate Change, 7(6), 407–411, 
doi:10.1038/nclimate3297, 2017. 

Giorgi, F., and X. Bi, Time of emergence (TOE) of GHG-forced pre- cipitation change hot-spots, 

Geophysical Research Letters, 36(6), doi: 10.1029/2009gl037593, 2009. 

Hansen, J., I. Fung, A. Lacis, D. Rind, S. Lebedeff, R. Ruedy, G. Russell, and P. Stone, Global 

climate changes as forecast by Goddard institute for space studies three-dimensional 
model, Journal of Geophysical Research, 93(D8), 9341, doi:10.1029/jd093id08p09341, 
1988. 

Harrington, L., D. Frame, E. Hawkins and M. Joshi, Seasonal cycles enhance disparities between 

low- and high-income countries in exposure to monthly temperature emergence with 
future warming, Environ. Res. Lett., 12, 114039, doi: 10.1088/1748-9326/aa95ae, 2017 

Haustein, K., M. R. Allen, P. M. Forster, F. E. L. Otto, D. M. Mitchell, H. D. Matthews, and D. 

J. Frame, A real-time global warming index, Scientific Reports, 7(1), 
doi:10.1038/s41598-017-14828-5, 2017. 

Hawkins, E., and R. Sutton, Time of emergence of climate signals, Geophysical Research 

Letters, 39 (1), doi:10.1029/2011gl050087, 2012. 

Hawkins, E., S. Burt, P. Brohan, M. Lockwood, H. Richardson, M. Roy, and S. Thomas, Hourly 

weather observations from the Scottish highlands (1883–1904) rescued by volunteer 
citizen scientists, Geoscience Data Journal, doi:10.1002/gdj3.79, 2019. 

Hollis, D., M. McCarthy, M. Kendon, T. Legg, and I. Simpson, HadUK-grid—a new UK dataset 

of gridded climate observations, Geoscience Data Journal, doi:10.1002/gdj3.78, 2019. 

IPCC, Summary for Policymakers. In: Global Warming of 1.5°C. An IPCC Special Report on 

the impacts of global warming of 1.5°C above pre-industrial levels and related global 
greenhouse gas emission pathways, in the context of strengthening the global response to 
the threat of climate change, sustainable development, and efforts to eradicate poverty 
[Masson-Delmotte, V., P. Zhai, H.-O. Pörtner, D. Roberts, J. Skea, P.R. Shukla, A. 
Pirani, W. Moufouma-Okia, C. Péan, R. Pidcock, S. Connors, J.B.R. Matthews, Y. Chen, 
X. Zhou, M.I. Gomis, E. Lonnoy, T. Maycock, M. Tignor, and T. Waterfield (eds.)]. 
World Meteorological Organization, Geneva, Switzerland, 32 pp., 2018 

Kennedy, J. J., N. A. Rayner, R. O. Smith, D. E. Parker, and M. Saunby, Reassessing biases and 

other uncertainties in sea surface temperature observations measured in situ since 1850: 
1. measurement and sampling uncertainties, Journal of Geophysical Research, 116(D14), 
doi:10.1029/ 2010jd015218, 2011. 

Kincer, J. B., Is our climate changing? A study of long-time temperature trends, Monthly 

Weather Review, 61 (9), 251– 259, doi:10.1175/1520-0493(1933)61

⟨

251:ioccas

⟩

2.0.co;2, 

1933. 

King, A. D., and L. J. Harrington, The inequality of climate change from 1.5 to 2c of global 

warming, Geophysical Research Letters, 45 (10), 5030–5033, 
doi:10.1029/2018gl078430, 2018. 

Page 10

Figures and page graphics, page 10

 

 

King, A. D., M. G. Donat, E. M. Fischer, E. Hawkins, L. V. Alexander, D. J. Karoly, A. J. 

Dittus, S. C. Lewis, and S. E. Perkins, The timing of anthropogenic emergence in 
simulated climate extremes, Environmental Research Letters, 10(9), 094,015, 
doi:10.1088/1748-9326/10/9/094015, 2015. 

Lehner, F., C. Deser, and L. Terray, Toward a new estimate of “time of emergence” of 

anthropogenic warming: Insights from dynamical adjustment and a large initial-condition 
model ensemble, Journal of Climate, 30(19), 7739–7756, doi:10.1175/jcli-d-16-0792.1, 
2017. 

Lenssen, N. J. L., G. A. Schmidt, J. E. Hansen, M. J. Menne, A. Persin, R. Ruedy, and D. Zyss, 

Improvements in the GISTEMP uncertainty model, Journal of Geophysical Research: 
Atmospheres, 124(12), 6307–6326, doi: 10.1029/2018jd029522, 2019. 

Madden, R. A., and V. Ramanathan, Detecting climate change due to in- creasing carbon 

dioxide, Science, 209 (4458), 763–768, doi:10.1126/science. 209.4458.763, 1980. 

Maher, N., et al., The Max Planck Institute grand ensemble: Enabling the exploration of climate 

system variability, Journal of Advances in Modeling Earth Systems, 11(7), 2050–2069, 
doi:10.1029/2019ms001639, 2019. 

Mahlstein, I., R. Knutti, S. Solomon, and R. W. Portmann, Early onset of significant local 

warming in low latitude countries, Environmental Research Letters, 6(3), 034,009, 
doi:10.1088/1748-9326/6/3/034009, 2011. 

Mahlstein, I., G. Hegerl, and S. Solomon, Emerging local warming signals in observational data, 

Geophysical Research Letters, 39 (21), doi: 10.1029/2012gl053952, 2012. 

Mahlstein, I., J. S. Daniel, and S. Solomon, Pace of shifts in climate regions increases with 

global temperature, Nature Climate Change, 3, 739, doi: 10.1038/nclimate1876, 2013 

Min, S.-K., X. Zhang, F. W. Zwiers, and G. C. Hegerl, Human contribution to more-intense 

precipitation extremes, Nature, 470(7334), 378–381, doi: 10.1038/nature09763, 2011. 

Monllor-Hurtado, A., M. G. Pennino, and J. L. Sanchez-Lizaso, Shift in tuna catches due to 

ocean warming, PLOS ONE, 12(6), e0178,196, doi: 10.1371/journal.pone.0178196, 
2017. 

Morice, C. P., J. J. Kennedy, N. A. Rayner, and P. D. Jones, Quantifying uncertainties in global 

and regional temperature change using an ensemble of observational estimates: The 
HadCRUT4 data set, Journal of Geophysical Research: Atmospheres, 117(D8), n/a–n/a, 
doi:10.1029/2011jd017187, 2012. 

Pall, P., M. R. Allen & D. A. Stone, Testing the Clausius–Clapeyron constraint on changes in 

extreme precipitation under CO

2

 warming, Climate Dynamics, 28, 351, doi: 

10.1007/s00382-006-0180-2, 2007 

Rohde, R., R. A. Muller, R. Jacobsen, E. Muller, S. Perlmutter, A. Rosenfeld, J. Wurtele, D. 

Groom and C. Wickham, A New Estimate of the Average Earth Surface Land 
Temperature Spanning 1753 to 2011, Geoinfor Geostat, 1, 1, doi:10.4172/2327-
4581.1000101, 2013 

Rojas, M., F. Lambert, J. Ramirez-Villegas, and A. J. Challinor, Emergence of robust 

precipitation changes across crop production areas in the 21st century, Proceedings of the 

Page 11

Figures and page graphics, page 11

 

 

National Academy of Sciences, 116 (14), 6673– 6678, doi:10.1073/pnas.1811463116, 
2019. 

Schneider, U., P. Finger, A. Meyer-Christoffer, E. Rustemeier, M. Ziese, and A. Becker, 

Evaluating the hydrological cycle over land using the newly- corrected precipitation 
climatology from the global precipitation climatology centre (GPCC), Atmosphere, 
8(12), 52, doi:10.3390/atmos8030052, 2017. 

Sutton, R., E. Suckling, and E. Hawkins, What does global mean tempera- ture tell us about local 

climate?, Philosophical Transactions of the Royal Society A: Mathematical, Physical and 
Engineering Sciences, 373(2054), 20140,426, doi:10.1098/rsta.2014.0426, 2015. 

Thirumalai, K., P. N. DiNezio, Y. Okumura, and C. Deser, Extreme temperatures in southeast 

Asia caused by El Nino and worsened by global warming, Nature Communications, 8(1), 
doi:10.1038/ncomms15531, 2017. 

Wigley, T. M. L. and P. D. Jones, Detecting CO

2

-induced climatic change, 

Nature

, 292, 205, 

1981 

Zappa, G. and T.G. Shepherd, Storylines of Atmospheric Circulation Change for European 

Regional Climate Impact Assessment. J. Climate, 30, 6561–6577, doi: 10.1175/JCLI-D-
16-0807.1, 2017 

Zhang, H.-M., et al., Updated temperature data give a sharper view of climate trends, Eos, 100, 

doi:10.1029/2019eo128229, 2019. 

Zhang, X., F. W. Zwiers, G. C. Hegerl, F. H. Lambert, N. P. Gillett, S. Solomon, P. A. Stott, and 

T. Nozawa, Detection of human influence on twentieth-century precipitation trends, 
Nature, 448(7152), 461–465, doi: 10.1038/nature06025, 2007. 

 

 

 

 

 

Page 12

Figures and page graphics, page 12

 

 

 

 

 

Figure 1: 

Emergence of global and local temperature change from 1850-2018. (top) GMST 

(grey), smoothed with a 41-year lowess filter (black). (bottom) Oxford annual temperature (grey) 
and scaled smoothed GMST (black). The correlation between Oxford temperatures and 
smoothed GMST is 0.67, and if the Oxford data is also smoothed with a 41-year lowess filter the 
correlation increases to 0.98. The shaded bands indicate 1 and 2 standard deviations of the noise. 

 

 

Page 13

Figures and page graphics, page 13

 

 

 

 

 

 

Figure 2: 

Signal, noise (both in K) and S/N for observed annual mean temperature change in the 

Berkeley Earth dataset. Many tropical regions show the smallest signal, but also the smallest 
noise and largest S/N. Grey regions denote lack of sufficient data. The S/N values in stippled 
areas are not significantly different from zero. 

 

 

Page 14

Figures and page graphics, page 14

 

 

 

Figure 3: 

Signal-to-noise ratio for monthly average temperatures, for the climatologically 

warmest (left) and coldest (right) months at each grid point. Grey regions denote lack of 
sufficient data. The S/N values in stippled areas are not significantly different from zero.

  

 

Figure 4: 

Signal-to-noise ratio for annual mean precipitation over land using the GPCC dataset. 

Blue colours denote regions becoming wetter, and red colours denote regions that are becoming 
drier. Grey regions are either unobserved (oceans) or deserts (<250mm/year). Stippling indicates 
where the regression parameter is not statistically significant from zero. 

 

 

Page 15

Figures and page graphics, page 15

 

 

 

Figure 5: 

Signal (left) and signal-to-noise ratio (right) for annual mean precipitation over the UK (top row, 1862-

2017) and extreme daily rainfall (RX1day, bottom row, 1891-2017) using the HadUK-Grid dataset. The signal is 
presented in units of % per K of GMST change. Blue colours denote regions becoming wetter, and red colours 
denote regions that are becoming drier. Stippling in the S/N panels indicates where the regression parameter is not 
statistically significant from zero. 
 

 

Page 16

Figures and page graphics, page 16

 

 

 

Geophysical Research Letters 

Supporting Information for 

Observed emergence of the climate change signal: from the familiar to the unknown 

E. Hawkins

1

, D. Frame

2

, L. Harrington

3

, M. Joshi

4

, A. King

5

, M. Rojas

6

, and R. Sutton

1

  

1

 National Centre for Atmospheric Science, Department of Meteorology, University of Reading. 

2

 New Zealand Climate Change Research Institute, Victoria University of Wellington, 

Wellington 6012, New Zealand. 

3

 

4

 Climatic Research Unit, School of Environmental Sciences, University of East Anglia, UK 

5

 School of Earth Sciences and ARC Centre of Excellence for Climate Extremes, University of 

Melbourne. 

6

 Department of Geophysics and Centre for Climate and Resilience Research, CR2, University of 

Chile 

 

Contents of this file  

 

Text S1 to S3 
Figures S1 to S9 
 

 

 

Page 17

Figures and page graphics, page 17

 

 

S.1 Shifting distributions 
 

The emergence of a signal can be visualised using shifting normal distributions (Fig. S1). 

Frame et al.

 (2017) 

described S/N>1 as a shift to an ‘

unfamiliar

’ climate, S/N>2 as an ‘

unusual

’ climate and S/N>3 as an ‘

unknown

’ 

climate, in terms of an individual’s lifetime. We add the term ‘

inconceivable

’ for S/N>5, as the new mean climate 

would be experienced once every 3 million years in the old climate.  
 
Two regional average examples are shown in Fig. S2, for tropical America and northern America, highlighting the 
differences in signal and noise characteristics. Even though northern America has a larger signal, the change is more 
apparent in tropical America. 
 

 

Figure S1: 

Shifting a normal distribution by 0 (black) to 6 (dark red) standard deviations.

   

 
 

 

 
Figure S2: 

Two regional examples of how observed temperature changes have become apparent, using the Berkeley 

Earth temperature dataset. The red shaded bands represent 1 and 2 standard deviations of the noise.

 

 
 

 

Page 18

Figures and page graphics, page 18

 

 

S.2 Using model simulations to test the emergence methodology 
 

We can test the robustness of the methodology to estimate the S/N using a large ensemble of model simulations. 

Maher et al.

 (2019) describe the 100-member ensemble of the MPI GCM, from which we use the simulated SAT for 

the historical period (1850-2005), extended to 2018 with the RCP4.5 scenario. First, we apply the same 
methodology used for the observations to each ensemble member individually. The ensemble mean S/N, which is 
expected to be smoother than the observed S/N due to averaging, is shown in Fig. S3a, and the spread in S/N across 
the ensemble is shown in Fig. S3c. The uncertainty in S/N is generally between 0.2-0.4 over land, which is typically 
far smaller than the mean S/N. The maritime continent, North Atlantic and Southern Ocean are regions with largest 
uncertainty in this GCM. The percentage uncertainty in S/N is less than 30% over most land areas (Fig. S3d). A 
simpler approach, which is not possible using observations, is to calculate the S/N by averaging the simulated 
temperature anomaly patterns in 2018, relative to the mean of 1850-1900, from all ensemble members, and dividing 
by the standard deviation of the 2018 anomalies (Fig. S3b). This pattern is virtually identical to Fig. S3a, 
highlighting that the regression approach produces S/N estimates that are robust. These results also demonstrate that 
the uncertainty in S/N due to simulated internal variability is relatively small. 
 
Note that the patterns of simulated S/N in this ensemble are noticeably different from the observed patterns. One 
important example is in parts of west Africa where the MPI ensemble S/N is close to zero but is larger than 5 in the 
observations. India also has a low S/N in the ensemble, but significant values in the observations. This finding 
highlights the benefit of using the observations alone, as in the current study. 
 
Fig. S4 shows the same maps for simulated precipitation change in the MPI ensemble. Again, the two methods 
produce similar patterns (Fig. S4a, b), with the ensemble method showing slightly larger values. The simulated 
uncertainty in S/N due to internal variability is typically 0.3-0.4 over land regions. The patterns are again different 
from that derived from the observations, especially in west Africa which is significantly wetter in the simulations 
but drier in the observations. 

 

Page 19

Figures and page graphics, page 19

 

 

 

 

Figure S3: 

Testing the S/N methodology using the MPI Large Ensemble (

Maher et al.

 2019). (top left) S/N 

calculated as for the observations in each individual ensemble member, averaged across the 100-members. (top 
right) Mean simulated temperature in 2018 minus the average of 1850-1900 across all ensemble members, divided 
by the standard deviation of simulated temperature in 2018. (bottom left) Standard deviation in the S/N estimated 
using the observational method across the 100-members. (bottom right) The percentage uncertainty in S/N, i.e. 
bottom left panel divided by top left. 

 

Page 20

Figures and page graphics, page 20

 

 

 

Figure S4: 

as Fig. S3 for precipitation

. 

 
 

 

Page 21

Figures and page graphics, page 21

 

 

 
S.3 Additional metrics 
 

Figure S5 shows the fraction of land area which has a S/N for temperature exceeding the value indicated, using the 
Berkeley Earth dataset. For the annual mean, around 15% of the land area has a S/N larger than 5, and 40% shows a 
S/N larger than 2 for the warmest climatological month of the year. The warmest months tend to show larger S/N 
values than the coldest months. 
 
Figure S6 repeats the S/N temperature analysis using other datasets: HadCRUT4 (

Morice et al.

 2012), 

Cowtan & 

Way

 (2014, hereafter CW14) infilled version of HadCRUT4, GISTEMP (

Lenssen et al.

 2019) and NOAA 

GlobTemp (

Zhang et al. 

2019). For this sensitivity test we have used the same smoothed GMST from Berkeley 

Earth in all cases. These datasets generally produce similar patterns to that from Berkeley Earth (Fig. 2c), but with 
varying amplitudes. NOAA GlobTemp has larger S/N values in the tropics than the other datasets and Berkeley 
Earth has larger S/N for the south-east USA. There are other notable differences for west Africa and parts of south 
America, mainly due to different estimates for the signal, rather than the noise (not shown). There is consistent 
agreement that the tropical Atlantic and Indian Oceans exhibit the highest S/N for the ocean areas, and that there has 
been very little warming overall in the central North Atlantic. 
 
Figure S7 shows the S/N patterns for precipitation in different seasons, highlighting that the west Africa signals are 
present in all seasons except DJF, and the south-west Australia drying signal is mainly present in JJA. The wetter 
northern latitude signal is mainly present in DJF and MAM. 
 
Figure S8 shows the S/N patterns for UK mean precipitation in different seasons. There are tendencies towards 
wetter seasons, except for JJA where the S/N is rarely significant. Note that the observed signal in southern UK is 
for drier summers but it has not yet emerged. 
 
Figure S9 shows the UK mean RX1day time-series with maps for two example years. 
 

 

Page 22

Figures and page graphics, page 22

 

 

 

 

Figure S5: 

The fraction of land area with an observed temperature S/N larger than the ratio shown, for different 

seasons, the annual average, and warmest and coldest months (using the Berkeley Earth dataset). 
 
 

 

Page 23

Figures and page graphics, page 23

 

 

 

 

 

 

 
Figure S6: 

Observed S/N for temperature using the CW14 dataset (top left), HadCRUT4 (top right), GISTEMP 

(bottom left) and NOAA GlobTemp (bottom right). Stippled cells indicate that the regression coefficient is not 
statistically significant. Grey regions are where there is less than 100 years of data in that location for that dataset. 

 
 
 
 
 

Page 24

Figures and page graphics, page 24

 

 

 

Figure S7: 

Signal-to-noise for precipitation in different seasons. Grey regions are either unobserved (oceans), have 

a seasonal precipitation of less than 62.5mm or annual precipitation less than 250mm. Stippled regions denote areas 
where the regression parameter is not statistically significant. 

 

Page 25

Figures and page graphics, page 25

 

 

 

 
Figure S8: 

Signal-to-noise for UK mean precipitation in different seasons. Stippled regions denote areas where the 

regression parameter is not statistically significant. 

 

Page 26

Figures and page graphics, page 26

 

 

 

 
Figure S9:

 UK extreme rainfall (RX1day, mm): average across the UK (1891-2017, black line) and regression on 

GMST (red dashed line), and maps for two example years (1968 and 2003). 1968 shows the effect of three 
significant storm events, in contrast to 2003 which mainly shows larger rainfall over higher orographic features.