Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/ doi:10.5194/esd-8-177-2017 Author(s) 2017. CC Attribution 3.0 License.
Roger N. Jones and James H. Ricketts
Victoria Institute of Strategic Economic Studies, Victoria University, Victoria 8001, Melbourne, Australia
Correspondence to: Roger N. Jones ([email protected])
Received: 13 August 2016 Discussion started: 23 August 2016Revised: 20 February 2017 Accepted: 21 February 2017 Published: 16 March 2017
Abstract. Interactions between externally forced and internally generated climate variations on decadal timescales is a major determinant of changing climate risk. Severe testing is applied to observed global and regional surface and satellite temperatures and modelled surface temperatures to determine whether these interactions are independent, as in the traditional signal-to-noise model, or whether they interact, resulting in step-like warming. The multistep bivariate test is used to detect step changes in temperature data. The resulting data are then subject to six tests designed to distinguish between the two statistical hypotheses, hstep and htrend. Test 1: since the mid-20th century, most observed warming has taken place in four events: in 1979/80 and 1997/98 at the global scale, 1988/89 in the Northern Hemisphere and 196870 in the Southern Hemisphere. Temperature is more step-like than trend-like on a regional basis. Satellite temperature is more step-like than surface temperature. Warming from internal trends is less than 40 % of the total for four of ve global records tested (18802013/14). Test 2: correlations between step-change frequency in observations and models (18802005) are 0.32 (CMIP3) and 0.34 (CMIP5). For the period 19502005, grouping selected events (1963/64, 196870, 1976/77, 1979/80, 1987/88 and 199698), the correlation increases to 0.78. Test 3: steps and shifts (steps minus internal trends) from a 107-member climate model ensemble (20062095) explain total warming and equilibrium climate sensitivity better than internal trends. Test 4: in three regions tested, the change between stationary and non-stationary temperatures is step-like and attributable to external forcing. Test 5: step-like changes are also present in tide gauge observations, rainfall, ocean heat content and related variables. Test 6: across a selection of tests, a simple stepladder model better represents the internal structures of warming than a simple trend, providing strong evidence that the climate system is exhibiting complex system behaviour on decadal timescales. This model indicates that in situ warming of the atmosphere does not occur; instead, a store-and-release mechanism from the ocean to the atmosphere is proposed. It is physically plausible and theoretically sound. The presence of step-like rather than gradual warming is important information for characterising and managing future climate risk.
1 Introduction
The dominant paradigm for how the climate changes over decadal timescales is based on the standard signal-to-noise model, where the externally driven signal of climate change forms a trend surrounded by the internally generated noise of climate variability. Here, the external driver of interest is radiative forcing produced by anthropogenic greenhouse gas emissions, mediated by other anthropogenic emissions such as sulfate aerosols and black carbon. This paradigm is
Reconciling the signal and noise of atmospheric warming on decadal timescales
widely represented by trend analysis, which extracts a monotonic signal from a noisy time series (e.g. North et al., 1995;Hegerl and Zwiers, 2011; Santer et al., 2011). The resulting methodology dominates climate practice, forming the basis for detection and attribution, projection, prediction and char-acterisation of climate risk.
However, it is not the only theoretically plausible representation of a changing climate (Palmer, 1999; Branstator and Selten, 2009; Solomon et al., 2011; Kirtman et al., 2013).The two main hypotheses that describe how externally driven
Published by Copernicus Publications on behalf of the European Geosciences Union.
178 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
and internally generated climate may be related over decadal timescales are (Corti et al., 1999; Hasselmann, 2002)
H1. Externally forced climate change and internally generated natural variability change independently of each other.
H2. They interact, for example, where patterns of the response project principally onto modes of climate variability (Corti et al., 1999) or form a two-way relationship (Branstator and Selten, 2009).
These interactions can lead to a range of different outcomes. For global mean surface temperature, the signal is generally portrayed as following a linear pathway that conforms to the relationship T = F , where T is temperature,
F is forcing and is a constant related to feedback processes (Ramaswamy et al., 2001; Andrews et al., 2015). This is widely accepted for both H1 and H2 over longer timescales (e.g. > 50 years), but how boundary-limited and initial conditions uncertainties combine over shorter timescales remains unclear.
For H1, if the response to external forcing is considered to be independent of variability over shorter timescales (< 50 years), the trend model will hold, despite often being obscured by variability. Such variability is generally represented as stochastic behaviour in annual to decadal phenomena, where teleconnections, lagged effects and regime changes all potentially interact (Solomon et al., 2011; Kirtman et al., 2013). Alternatively, instead of a gradual line or curve, a segmented trend is sometimes proposed, where the signal of atmospheric warming is modied by varying decadal regimes governing oceanic sources and sinks of heat (Meehl et al., 2013; Cahill et al., 2015; Trenberth, 2015). All these statistical models are linked by the representation of warming as a gradual process, leading to the gradualistic narrative of change (Jones et al., 2013).
The potential behaviour of warming under H2 has many possible permutations because the signal may project onto the regime-like structures of decadal climate variability, or it may dynamically modify those structures. Although a number of nonlinear and often abrupt changes in climate are recognised as part of decadal change, these are overwhelmingly attributed to changes in climate variability. Here, we deal with one such type of response, manifesting as step changes. Step changes have been detected in warming and related climatic variables by several different methods (Jones, 2010; Reid and Beaugrand, 2012; Jones et al., 2013; Belolipetsky, 2014; Belolipetsky et al., 2015; Bartsev et al., 2016; Reid et al., 2016); in one case, step-like warming over SE Australia was attributed to anthropogenic forcing (Jones, 2012). The purpose of this paper is to detect step changes in a range of temperature records and to apply severe testing to steps and trends to determine which carries the greater part of the warming signal. The results are used to determine whether H1 or H2 is the more viable hypothesis and, if the signal is shown to be non-gradual, to explore the nature of the interaction between external forcing and internal variability.
We apply a methodology combining theoretical-mechanistic and statistical-inductive reasoning to test which statistical model, step or trend, better represents the warming signal on decadal timescales. It is applied to the substantive null of model adequacy approach described by Mayo and Cox (2010) as part of severe testing principles articulated by Mayo and Spanos (2010). Although a test may provide a small p value for the null hypothesis, other tests may do so as well, in which case the hypothesis that test represents is provisional. Support for both H1 and H2 in the literature shows this to be the case. The presence of several statistical models with similar p values also shows that there are viable alternatives to the simple trend model (Seidel and Lanzante, 2004).
A substantive null of model adequacy is where a test closely supports a hypothesis and where a rival test has a high probability of detecting a specic discrepancy from that hypothesis, if that rival hypothesis is correct (Mayo and Cox, 2010). The testing model can be adapted for a single hypothesis or rival hypotheses. If the rival test fails, then the original hypothesis succeeds; if the rival test succeeds, then the original test should also have a low probability of detecting a specic discrepancy from the rival hypothesis. When rival hypotheses are being tested, conrmation and falsication provide two different views of the same issue.
The theoretical-mechanistic component describes plausible, alternative physical processes in the climate system required to sustain steps and trends respectively. Step changes are measured using an objective rule-based multistep adaptation of the bivariate test of Maronna and Yohai (1978) to analyse regional and global surface air temperature, global satellite temperature of the lower troposphere and global mean temperature from the CMIP3 and CMIP5 climate model archives. The data produced by those analyses are then subject to six tests designed to distinguish between steps and trends as the main drivers of the anthropogenic climate signal over decadal timescales.
2 Methodology
The process of theoretical-mechanistic and statistical-inductive reasoning requires matching scientic hypotheses (H) with statistical hypotheses (h) in order to distinguish between alternative hypotheses. The next few sections detail how this has been carried out. This employs a hierarchy of models between theory and data, as suggested by Suppes (1962) and articulated by Haig (2016). Underlying theory is used to inform plausible mechanisms for alternative types of change (steps and trends), experimental analyses test those mechanisms and statistical models that detect those alternative types of change are used to prepare climate data for testing. By and large, statistical models are used to undertake error testing, whereas the experimental analyses under-
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 179
take probative testing designed to provide evidence for the hypotheses being tested.
Here, linearity of response is dened by the T = F re
lationship, where forcing produces a continuous response in temperature that can be masked by climate variability. Even if the function increases over time (e.g. Rypdal and Rypdal, 2014; Andrews et al., 2015), the response will be gradual but will accelerate with increasing forcing. This relationship is also used to dene the concept of model equilibrium climate sensitivity (ECS), measured as the atmospheric warming caused by a forcing of 2 CO2 in the atmospheric com
ponent of a climate model. The relationship between step-like and trend-like behaviour in climate model output and ECS can be used to test how strongly each responds to radiative forcing. The results will show whether forcing produces gradual or episodic warming over decadal timescales.
2.1 Development of physical mechanisms for probative testing
Application of a theoretical-mechanistic process starts with well-agreed theoretical positions (core theory) and then builds on those theories to explore alternative mechanisms required to support competing hypotheses. The exploration of plausible mechanisms produces probative criteria for severe testing. This paper cannot undertake a full survey of the theory behind anthropogenic global warming, but the trapping of heat by added greenhouse gases, creating an imbalance between the surface and the top of the atmosphere and between the equator and the poles, is widely agreed on as the foundational theory, i.e. radiative transfer theory and global warming resulting from the enhanced greenhouse effect (IPCC, 2013). However, between the time when heat is trapped in the atmosphere and when it is measured as a change in temperature there is a gap in understanding, which has competing explanations. These explanations focus on where that trapped heat is stored within the climate system and how it is subsequently distributed. Because H1 implies a gradual signal and H2 a discontinuous or episodic signal, represented here as step-like change, these pathways will be distinctly different.
For H1, close adherence to a warming trend implies that the atmosphere warms gradually. If so, this must occur via either or both of the following processes:
1. A measurable proportion of radiatively forced anthropogenic warming trapped in the atmosphere is retained in situ, as represented by models of radiative convective transfer (Ramanathan and Coakley, 1978), gradually warming the air mass, especially over land. Such warming would also be expected to produce a trend in lower troposphere satellite temperatures since the air mass warms gradually from the surface.
2. Most of the heat trapped by anthropogenic greenhouse gas forcing is absorbed by the ocean, with the ocean re-
taining an estimated 93 % of historically trapped heat (Levitus et al., 2012; Roemmich et al., 2015). Models of upwelling diffusion assume a constant release of heat into the atmosphere (Raper et al., 2001, 2002) and the assumption of gradual release follows through into much of the literature. Recent papers discuss the role of decadal variability within the oceans mediating trends in atmospheric warming (England et al., 2014;Watanabe et al., 2014; Dai et al., 2015; Meehl, 2015;Trenberth, 2015; Meehl et al., 2016) through variations in ocean surface temperatures and/or overturning processes.
This combination of processes forms the dominant paradigm, where the anthropogenic warming signal is widely considered largely as forming a monotonic trend (Swanson et al., 2009; Zhou and Tung, 2013; Ji et al., 2014). However, mental (conceptual) models held by individual scientists vary widely (Benestad, 2016). Under a scenario of changing decadal regimes, it is also possible that internally driven step changes could be detected in temperature time series, forming a stepladder, as suggested by Trenberth (2015). However, if H1 were to hold, these would have to be unrelated to forcing.
Non-gradual warming (H2) requires mechanisms such as regime change combining with storage and release processes.On decadal timescales, oceanatmosphere interaction is the only realistic source for such changes. If warming is mediated by the hydrothermal oceanatmosphere system, it could be entrained by the nonlinear processes involved in the distribution of energy skywards and polewards from the equator through quasi-oscillatory systems (Ozawa et al., 2003;Lucarini and Ragone, 2011). Lucarini and Ragone (2011) describe the overall process of distribution of heat energy within the climate system as the generation of entropy, where moist static energy is transformed into mechanical energy like a heat engine. This could ip between different states, modulated by Lorenzian strange attractors as described by Palmer (1993). One important distinguishing characteristic for nonlinear behaviour in a changing climate is whether it is internally generated and essentially random, whereas if it is forced, the response will be related to changing boundary conditions (Lorenz, 1975; Hasselmann, 2002). Distinguishing between these possibilities is the focus of the testing regime: whether gradual or step-like changes provide the better explanation for the response to external forcing.
2.2 Development of severe testing
The aim of severe testing is to produce highly probed (evidential) rather than highly probable results (Mayo, 2005). A hypothesis H passes a severe test T with data x if (Mayo and Spanos, 2010)
1. x agrees with H.
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
Test 4 Can step-like change be identied using attribution methods?
Test 5 Do other climate variables also undergo step changes?
Test 6 Are temperature time series more step-like or trend-like?
The rst four tests can be considered largely probative, where hstep and htrend are tested to determine whether H1 or
H2 provides the better explanation for the relationship between external forcing and internal variability. The last two focus mainly on error testing to see how well hstep and htrend
explain the climate data. The combination of different tests means that deriving a single probability through an objective process is not possible. The procedure we follow here uses a two-sided test between hstep and htrend as representatives of H1 and H2. Paraphrasing Mayo and Spanos (2010) to address the results, with very high probability, Tests 16 would produce a result that accords less well with H2 than does H1, if H2 were false or incorrect (and conversely).
2.3 Statistical testing
2.3.1 The multistep MaronnaYohai bivariate test
The MaronnaYohai bivariate test (MYBT, Maronna and Yohai, 1978) is used to detect step changes in temperature data. This test has been widely used to detect inhomogeneities in climate variables (Potter, 1981; Bcher and Dessens, 1991; Kirono and Jones, 2007; Sahin and Cigizoglu, 2010), decadal regime shifts in climate-related data and step changes in a wide range of climatic time series (Buishand, 1984; Vivs and Jones, 2005; Boucharel et al., 2011; Jones, 2012; Jones et al., 2013). One of us (Jones) has been using it for 25 years, both for adjusting inhomogeneous data (Jones, 1995; Kirono and Jones, 2007) and for detecting abrupt changes in climate variables. Surprisingly, the MYBT is rarely included in reviews of change point analysis techniques (Rodionov, 2005; Reeves et al., 2007) despite being on par or better than other techniques (Vivs and Jones, 2005). For example, it performed similarly to the STARS test in Jones et al. (2013) but has the advantage of not needing tuning and being able to accommodate a reference data set, providing a degree of exibility that few other tests have.That made it our testing model of choice, especially because all six tests used here compare step changes in time series to a null reference, and Test 4 assesses step changes between correlated variables.
The test was adapted from being able to only assess single change points by developing an objective set of rules that would detect a minimal and stable conguration of multiple step changes. Previously, this involved a trial-and-error process of constructing a robust set of step changes one at a time.A multistep, rule-based application of the MYBT was devel-
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
180 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
2. with very high probability, test T produces a result that accords less well with H than does x, if H is false or incorrect.
Two sets of data are produced, representing competing statistical hypotheses hstep and htrend. These are linked to rival hypotheses H1 and H2. Previous statistical testing of alternative structures for warming has been inconclusive. For example, when Seidel and Lanzante (2004) tested trends, steps, segmented trends, and step-and-trend statistical models, no single model stood out. They concluded that detection and attribution studies should consider abrupt changes.Studies that extract short-term components of climate variability from time series producing a more trend-like result (Foster and Rahmstorf, 2011; Werner et al., 2015) or decompose temperature time series into separate signal and noise components (Wu et al., 2011; Yao et al., 2016) all implicitly assume H1. Consequently, the exact nature of change on decadal timescales remains an open question (Trenberth, 2015). If warming conforms to a long-term complex trend and is additive (Marvel et al., 2015), such studies will only produce a trend-like output because they are not congured to detect alternative structures. However, because they are framed on H1, these tests do not show that such structures do not exist.
Therefore, htrend has never been severely tested to the point where its alternatives have been eliminated. The usual null hypothesis for htrend is no trend has emerged from background variability. Accordingly, the null hypothesis testing of trends is usually carried out assuming H1. Where step changes are detected, they are generally attributed to internal variability. However, non-gradual change on decadal timescales has become part of the climate wars, being used to challenge global warming theory on the basis that if ob-served change is not gradual, climate change is either disproved or overstated (e.g. Legates et al., 2015). Evidence of nonlinear change, such as step change, is therefore widely associated with challenges to global warming theory (e.g. see Skeptical Science, 2015). This asymmetry in null hypotheses means that severe testing needs to cover both H1 and H2, testing hstep against htrend.
The following six tests are used to test the relationship between gradual and step-like change and their responses to external forcing:
Test 1 What patterns of step changes can be detected in temperature observations? Do particular dates and locations line up with known events or processes?
Test 2 Do models forced by historical emissions reproduce the patterns of steps changes shown in observations?
Test 3 What is the relationship between different components of change steps, internal trends and shifts to each other and to total warming and ECS?
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 181
oped to carry this out (Ricketts, 2015; see Supplement for details).
The test adapts the formulation of Bcher and Dessens (1991) by testing a single serially independent variate (xi) against a reference variate (yi) using a random time series following Vivs and Jones (2005). The important outputs of the test in a time series of length N are (1) the Ti statistic, which is dened for times i < N, (2) the Ti0 value, which is the maximum Ti value, (3) i0, the time associated with Ti0, (4) shift at that time, and (5) p, the probability of zero shift. Note that i0 is the last year prior to the change. In this paper, we routinely give the year of change.
A single time series analysis consists of a screening pass, followed by a convergent pass. In both passes, we apply a resampling test to each segment being examined, where the test is repeated 100 times, resampling the random number reference series. The screening pass starts from the most signicant shift in a time series, determined using the resampling test and if p < 0.01, the series is divided into shorter time series either side of the step and these are tested until all steps have been detected. This is a recursive procedure whereby the rst steps detected may be inuenced by asyet-unlocated steps. The convergent pass then serially renes these segments to provide a causal sequence. The convergent process is repeated until a stable set of step changes is produced.
The analysis above is run 100 times. This procedure may produce several different but related solutions (sets of change dates); the most common solution is returned as the best estimate. Alternatives often indicate the presence of localised events embedded in larger-scale areally averaged data. Most historical temperature records analysed contain one or two stable congurations for surface temperature and zero or one for satellite temperature. Climate model data may produce a larger number of stable solutions, especially for the higher forcing scenarios.
Mean annual data for observations are considered serially independent, and in most cases applied in the paper, the MYBT is reliable. Deseasonalised quarterly and monthly data can be used to locate a shift within 1 year, but are not serially independent and are thus used here in combination with the t test either side of the change date to assess signicance. A resampling test that shufes data either side of a shift will also indicate whether a change point is abrupt or the time series is trend-like. As discussed in Sect. 4.3, 21st century model data are not serially independent under high rates of forcing.
For error testing, we routinely use thresholds of p < 0.01 for the bivariate test (exceptions are noted), and p < 0.01, p < 0.05 and non-signicant (NS, p > 0.05) for trend analysis and the t test.
2.3.2 Regional attribution
Regional attribution of step changes (Test 4) uses a technique detailed in Jones (2012). The basic methodology is suitable for continental mid-latitude areas where annual average maximum temperature (Tmax) is correlated with total rainfall (P ), and minimum temperature (Tmin) is correlated with Tmax (Power et al., 1998; Nicholls et al., 2004;Karoly and Braganza, 2005). For Central England Temperature, a largely maritime climate, diurnal temperature is assessed against precipitation instead of Tmax. The method uses the following steps:
1. Homogenous regional average data are obtained for Tmax, Tmin and P .
2. A period of stationary climate is calculated by testing when the relationship between Tmin and Tmax undergoes a statistically signicant step change. The relationship between Tmax and P will change at the same time or a later date.
3. Linear regressions are calculated between each pair (Tmax/P and Tmin/Tmax) for the stationary period.
4. Externally forced warming is estimated for the non-stationary period using these regressions.
5. The results are tested for step changes.
2.3.3 Observed data
Time series tested here are mean annual global air temperature anomalies from ve groups (NCDC, Peterson and Vose, 1997; GISS, Hansen et al., 2010; HadCRU, Morice et al., 2012; BEST, Rohde et al., 2012; C&W, Cowtan and Way, 2014), hemispheric temperatures from three groups (Had-CRU, NCDC and GISS) and zonal temperatures from two groups (NCDC and GISS) to see how prevalent step changes are, whether they coincide across different records and to investigate the relationship between step changes and trends.Lower tropospheric satellite temperatures from two groups (UAH, Christy et al., 2003; Christy et al., 2007; RSS, Mears and Wentz, 2009) are also tested.
For the regional data, Australian data were sourced from the Australian Bureau of Meteorology, Texas data from the National Climate Data Center and central England temperatures from the Met Ofce Hadley Climate Centre. Tide gauge records were sourced from the Permanent Service for Mean Sea Level and the ocean heat content records from the KNMI Climate Explorer. The specic records used are described in the Supplement.
2.3.4 Model data
Simulated mean global surface temperature from the CMIP3 and CMIP5 climate model archives is also tested. The analy-
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
182 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
sis is carried out in two parts. The rst part investigates simulated 20th century temperatures to determine how well the models reproduce the pattern of step changes in the observed data. The second part analyses how step changes evolve over the 21st century under the different radiative concentration pathways (RCPs). The output data are provided in the Supplement.
2.3.5 Metrics
Measurement of change where nonlinear behaviour is present is not an exact process, and there is no established terminology that carries commonly understood technical meanings; thus, we dene here a limited number of terms used in the paper. The MYBT measures total change between segments of a time series, ignoring any trend that may be present. We refer to these as steps. Internal trends are calculated between steps and the distance between the end of one trend and the start of the next is referred to as a shift. We call the process of calculating steps then trends the step-and-trend model. Steps, internal trends and shifts all provide data for severe testing.
Shifts and internal trends are not strictly additive; summed over a number of steps, they can add up to more or less than the change in temperature measured between the beginning and end of a series. These differences are largest in records containing reversals and negative trends.
The main phenomena analysed are (Fig. 1)
steps, which are the measurements of the whole change across a discontinuity, assuming stationarity produced by the bivariate test. This assumes no trend either side of the step.
internal trends, which are the measurements of trends between steps using ordinary least square trend analysis.
shifts, which are the measurements of the internal step between the end of a preceding trend and the beginning of the next trend.
trend / step ratio, which is the ratio between total internal trends and total steps in a multistep time series.Because shifts and internal trends are not additive, this measure gives a slight preference to trends over shifts as a ratio.
trend / shift ratio, which is the ratio between total internal trends and internal shifts (steps minus trends).
3 Results observations
3.1 Global and zonal temperatures
This section undertakes global, hemispheric and zonal analyses to determine temporal and spatial patterns of step changes
in observed temperature, consistent with Test 1. All series were tested from their earliest recorded date (1850 and 1880), and results from 1880 to 2014 are shown. Step changes meeting the p < 0.01 threshold in global and zonal temperatures show a great deal of structure. Downward steps occur in the late 19th and early 20th centuries, upward steps between 1912 and 1938, with one downward step in 1964.From 1968, upward steps dominate, with one exception in the high Southern Hemisphere (SH) latitudes in 2007 (Fig. 2).
The 1997 step change was global, with some regional steps occurring in 1996 and 1998. A global step change occurred in 1979/80, also registering in many regions, except for the Northern Hemisphere mid- and high latitudes. All other step changes occurred across more limited regions, with some being conned solely to land or to ocean. The 1997 step is the largest at 0.31 0.01 C. The 1979/80 step is the next largest
at 0.22 0.03 C. The greater variation in size of 1979/80 is
affected by the timing and size of previous steps and trends.
In the rst half of the 20th century, three global records show positive steps in 1920/21 and in 1937, and two in 1930 (Fig. 2). The GISS record also shows a downward step in 1902, coinciding with the Northern Hemisphere (NH) ocean, the tropics and the Southern Hemisphere. The two groups are based on the early 20th century differences: GISS, BEST and C&W in one group and HadCRU and NCDC in the other.The anomaly averaged from all ve records shows upward step changes in 1930, 1979 and 1997, coinciding with the HadCRU and NCDC records.
Differences emerge between ocean and land records. The global HadSST (HadCRU) records shifts in 1937, 1979 and 1997, whereas the ERSST (NCDC) records shifts in 1890, 1930, 1977, 1987 and 1997. Global land records from both CRU and NCDC shifted in 1920/21, 1980 and 1997. Northern hemispheric land and ocean step changes are consistent across three records: in 1924/25, 1987 and 1997. The NH ocean shows a downward step in 1902/03 and is less consistent between the two records tested for subsequent upward steps. The SH is consistent across 1937, 1979 and 1997, with two records showing a downward step in 1890 and an upward step in 1969.
The tropics show a downward step in 1902/03 and upward steps in 1926, 1979 and 1997. Three NH mid-latitude records show upward steps in 1920, 1921 or 1930, in 1987/88 and 1997/98. One zonal record also shows a downward step in 1964. The two NH high-latitude records show a single downward step in 1902 and in 2005, both step upwards in 1921 and 1994 and a single step upwards in 2005. The three SH mid-latitude records show a downward step in 1887 and one in 1902 and upward steps in 1933 or 1937, 1968 or 1970, 1977/1978 or 1984, and 1997 or 1998. SH high-latitude data are not very reliable, being absent for NCDC 6090 S. The
GISS 6490 S average anomaly steps downward in 1912 and upward in 1955.
Figure 3 shows the internal trends and their error significance for the ve global mean temperature records. Steps
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 183
0.4
1.0
1880 1900 1920 1940 1960 1980 2000
0.8
0.6
Anomaly (C)
0.4
0.2
0.0
0.2
Year
Figure 1. Record of mean annual surface temperature anomalies 18802014 from the Hadley Centre and Climate Research Unit (HadCRU), showing step changes (p < 0.01) and internal trends and shifts taken from the end of one internal trend to the start of the next across a step.
and trends are consistent for the last two periods 1979/80 to 1996 and 1997 to 2013/14 but diverge in the middle of the record due to differences in the timing and magnitude of steps and accompanying internal trends. Data quality may be an issue in the earlier parts of the record. For example, the versions of GISS data used here show ve steps in 1902, 1920, 1937, 1980 and 1997, whereas a version previous to 2013 stabilised at steps in 1930, 1979 and 1997, consistent with the average anomaly of all ve records. This indicates that the timing and magnitude of steps in the early 20th century can be inuenced by adjustments made to improve data quality. However, all global step change dates coincide with regional steps, showing that while the relative importance of dates associated with step changes may be different, the dates themselves are quite stable. This gives us added condence that we are not detecting false positives.
Internal trends are mainly p > 0.05 in the early record, the exception being the GISS 192037 period. The 1979/80 to 1996 trend is at p < 0.01 in two records (HadCRU and NCDC) and p < 0.05 in the other three records. The NH step change in 1987 seen in all three records tested strongly inuences this trend, which is examined further in the next section. The post-1997 period is p > 0.05 in two records and p < 0.05 in three records.
3.1.1 Step / trend and shift / trend ratios
There is no objective way to partition shifts and internal trends. Giving the rst preference to internal trends in calculating ratios gives a slight preference to gradual change in contrast to episodic change, preferring the methodological status quo. Expressed as a ratio between internal trends and steps, four global records range between 0.32 and 0.38, with the GISS record yielding a ratio of 0.62 due to the cool reversal in the early 20th century. For trends and shifts, the ratio
ranges between 0.44 and 0.58, with the GISS record being an outlier at 1.38.
Test 2 aims to determine whether trends or steps are more prominent at the regional level than at the global scale. The global trend / step ratio for the HadCRU record, for example, is 0.55 (0.30 / 0.55 C), 0.31 for the NH, 0.28 for the
SH and 0.33 for the tropics (30 N30 S); this is close to the average of the two hemispheres. When divided into land and ocean, the HadCRU and NCDC records show 0.90 and1.15 for land and 0.16 and 0.26 for ocean respectively, which shows the oceans to be more step-like and the land to have roughly equal measure. The SH ocean is very step-like (0.16) and SH land is less so (0.39). The mid-latitudes are also very step-like, as is the tropical ocean. High ratios (> 1) often involve a temporary cool reversal around the early 20th century.
This also holds for single steps on a regional basis. In 199798 the global shift was 0.16 0.01 C, a ratio of about
50 % compared to the step change of 0.32 C. For the NH, this ratio varied between 57 and 68 % for three land and three ocean data sets. For the NH mid-latitudes, land and ocean from two data sets (NCDC 3060 N, GISS 2444 N)
show a step / shift ratio that measures 0.43 / 0.44 C, close to a 1 : 1 ratio, which indicates no trend.
The more step-like character of both the oceans and the mid-latitudes is consistent with those areas being the loci of change in terms of decadal regimes and nonlinear equator-to-pole transport. This is inconsistent with the hypothesis of gradual warming. Varying shift dates and rates of change at regional scales contribute to the global record being more trend-like than individual regions.
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
184 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
1880
1890
1900
1910
1920
1930
1940
1950
1960
1970
1980
1990
2000
2010
NH land
Southern Hemisphere
SH ocean
SH land
Tropics
NH tropics
NH mid- lat
NH high lat
SH tropics
GISSHadCRUNCDCBESTC&WHadSSTERSSTCRUtemNCDCGISSHadCRUNCDCHadSST NHERSST NHCRUtem NHNCDC land NHGISSHadCRUNCDCHadSST SHERSST SHCRUtem SHNCDC land SH HadCRUNCDC 20 S20 N GISS 24 S24 N Tropic ocean HadSST 30 S30 N
GISS 024 N NCDC 030 N GISS 24 N44 N NCDC 30 N60 N
GISS 44 N64 N NCDC 60 N90 N GISS 64 N90 N GISS 024 S NCDC 030 S GISS 24 S44 S NCDC 30 S60 S GISS 44 S64 S SH high lats GISS 64 S90 S
NCDC 60 N60 S NCDC 20 N90 N GISS 24 N90 N NCDC 20 S90 S GISS 24 S90 S
Global
Global ocean
Global land
Northern Hemisphere
NH ocean
SH mid- lats
Other latitudinal zones
quarters of 1997 dominate the UAH global record, which is less so for the RSS record.
Quarterly anomalies for the RSS and UAH satellites and HadCRU and GISS surface mean global temperatures were compared to provide more precision on dates of step changes.Quarterly time series are affected by autocorrelation due to the El NioSouthern Oscillation (ENSO) for the bivariate test, making results robust for timing but not for probabilities for false positive (Type I) errors. The Students t test (two sided, unequal variance), which is insensitive to serial correlation, was used as a backup.
For the quarterly results, RSS shifted in December, January and February (DJF) 1987/88 by 0.11 C (p < 0.05 for
MYBT and p < 0.1 Students t test) and UAH shifts in DJF 1987/88 by 0.09 C (p > 0.05 for MYBT and p < 0.05 for Students t test). For surface temperature, HadCRU and GISS shifted in JJA 1987 by 0.14 and 0.15 C respectively
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
Figure 2. Dates of statistically signicant step changes (p < 0.01) in 18802014 for a range of mean annual temperature records. Downward steps are blue and upward steps are red. Records are sourced from the Goddard Institute of Space Studies (GISS); Hadley Centre and Climate Research Unit: HadCRU (land and ocean), HadSST (ocean), CRUtem (land); National Climatic Data Center: NCDC (land, land and ocean); ERSST (ocean); Berkeley Earth Surface Temperature (BEST); and Cowtan and Way (C&W). See Supplement for details.
3.2 Satellite-era records
A comparison of surface and lower tropospheric satellite temperatures straties records according to altitude and source of measurement, which is also consistent with Test 2.Satellite records of annual and seasonal lower-troposphere anomalies sourced from the RSS and UAH records beginning in December 1978 were analysed for step changes (1979 2014). Mean annual global and zonal temperatures show 1995 and 1998 as the two main step dates, with 1995 being more prominent at the global scale (Table 1). Seasonal temperatures were assessed to distinguish between these dates.For individual seasons, steps in 1995 are dominated by the NH June, July and August (JJA) and September, October and November (SON) periods, especially on land. This can be traced back to warm El Nio conditions in 1994/95. For the quarterly time series (4 seasons 36 years), the JJA and SON
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 185
1.0
1880 1900 1920 1940 1960 1980 2000
0.8
0.6
0.4
Anomaly (C)
0.2
0.0
-0.2
-0.4
Year
HadCRU
C&W
BEST
NCDC
GISS
Period Change Trend sig.
Period Change Trend sig.
Period Change Trend sig.
Period Change Trend sig.
Period Change Trend sig. 18801901 0.03 NS
18501929 0.01 NS
18501919 0.06 NS
18501920 0.05 NS
18801929 0.07 NS
1902 0.13
1930 0.26
1920 0.16
1921 0.15
1930 0.26
19021919 0.06 NS1920 0.18
19201936 0.12 NS
19201936 0.09 NS
19201936 0.15 **
1937 0.16
1937 0.19
1937 0.18
19301979 0.01 NS
19371979 0.10 *
19371979 0.06 NS
19301978 0.09 NS
19371979 0.01 NS
1980 0.20
1980 0.19
1980 0.20
1979 0.26
1980 0.26
19801996 0.19 **
19801996 0.17 *
19801996 0.16 *
19791996 0.18 **
19801996 0.17 *
1997 0.32
1997 0.33
1997 0.31
1997 0.30
1997 0.31
19972014 0.11 NS
19972013 0.17 *
19972013 0.14 *
19972014 0.10 NS
19972014 0.14 *
Trend 0.30
0.31
0.27
0.30
0.49
Step 0.78
0.85
0.86
0.82
0.80
Shift 0.55
0.59
0.63
0.54
0.36
Trend/step 0.38
0.37
0.32
0.36
0.62
Trend/shift 0.55
0.52
0.44
0.56
1.38
one global (1979/80) and one largely northern hemispheric (1987/88).
When all four records are plotted on a common baseline of 19791998, the surface and satellite temperatures display similar shifts but different internal trends (Fig. 4). Shown this way, the supposed differences between surface and satellite trends are largely removed. The satellite data contain signi-cant negative internal trends over 19791986 (RSS p < 0.01, UAH p < 0.05), but otherwise they show p > 0.05. The surface data show signicant positive internal trends over 1997 2014 (GISS p < 0.01, HadCRU p < 0.05), but otherwise they show p > 0.05. The decline post-1981 and lower trends in
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
Figure 3. Mean global anomalies of surface temperature with internal trends. The annual anomalies (dotted lines) from ve records (Had-CRU, C&W, BEST, NCDC, GISS) are taken from a 18801899 baseline. Internal trends (dashed lines) are separated by step changes detected by the bivariate test at the p < 0.01 error level. The size of each step (in red) and change in temperature of each internal trend (in black) is shown in the gure table along with its signicance, where NS is p > 0.05, is p > 0.01< 0.05 and is p < 0.01. Totals of trends, steps,
shifts (change from one trend to the next) and ratios are also shown.
(p < 0.01 for both tests). On an annual basis, the bivariate test registers 1987/88 at the p < 0.05 level. The lower error probabilities in the satellite records are due to the slightly lower shift size and higher variance. RSS shifted in JJA 1997 by 0.23 C, UAH shifted in DJF 1997/98 by 0.26 C, Had-
CRU in JJA 1997 by 0.26 C and GISS in SON 1997 by0.25 C (all with p < 0.01 for both tests). These four data sets show consistent shift dates in 1997 and similar shift dates in 1986/87, showing that the signicant step change in the NH is present at the global scale. This suggests that the period of accelerated trend noted by many for 19761998 (e.g. Tren-berth, 2015) is actually a period containing two step changes,
186 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
-0.4
0.8
Jan-79 Dec-83 Dec-88 Dec-93 Dec-98 Dec-03 Dec-08 Dec-13
0.6
Anomaly (C)
RSS
RSS trend
UAH
UAH trend
HadCRU
HadCRU trend
GISS
GISS trend
0.4
0.2
0.0
-0.2
Year
Figure 4. Quarterly mean lower tropospheric satellite (RSS, UAH) and surface (HadCRU, GISS) temperature anomalies on a common baseline 19792014. Annual anomalies (dotted lines) and internal trends (dashed lines) are separated by step changes.
Table 1. Dates of step changes for lower tropospheric satellite temperature anomalies, with annual time series and quarterly breakdowns in parentheses (DJF, MAM, JJA, SON), and quarterly time series. Data sources include remote sensing systems (RSS) and the University of Alabama, Huntsville (UAH).
Region Annual time series (quarterly breakdown) Quarterly time series
RSS UAH RSS UAH
Global land and ocean 1995 (98, 98, 95, 95) 1995 (97, 98, 94, 95) JJA 1997 SON 1997 Global land 1995 (95, 98, 95, 95) 1998 (98, 98, 94, 95) SON 1994 SON 1997 Global ocean 1998 (98, , 97, 95) 1995 (97, , , 95) JJA 1997 SON 1997 NH land and ocean 1995 (98, 98, 94, 94) 1998 (98, 98, 94, 94) JJA1997 SON 1997 NH land NA 1998 (98, 98, 98, 98) NA JJA 1997 NH ocean NA 1994 (, , , 94) NA JJA 1997 SH land and ocean 1995 (98, , , 95) 1995 (97, , 87, 95) SON 1997 SON 1997 SH land NA 1995 (95, , 91, 95) NA MAM 2002 SH ocean NA 1995 (97, , , 95) NA DJF 1998 Tropics land and ocean 1995 (, , , 93) (, , , 95) JJA 1997 JJA 1997 Tropics land 1995 (, , , 87) 1995 (98, , 95, 95) SON 1997 JJA1997 Tropics ocean 1995 (, , , 95) (, , , ) JJA 1997 NH ex-trop land and ocean 1998 (95, 98, 98, 94) 1998 (98, 98, 98, 94) SON 1997 DJF 1998 NH ex-trop land 1998 (, 98, 94, 94) 1998 (, 98, 98, 98) MAM 1994 DJF 1998 NH ex-trop ocean 1998 (99, 98, 98, 94) 1994 (02, 98, , 94) SON 1997 MAM 1998 SH ex-trop land and ocean 1998 (96, , , 95) 1996 (97, , , 95) DJF 1998 DJF 2001 SH ex-trop land 1995 (, , , ) 2001 (03, , , 02) JJA 1995 MAM 2002 SH ex-trop ocean 1998 (96, , , ) 1996 (97, , , 95) DJF 1998 DJF 1998 N polar land and ocean 1995 (03, 95, 98, 95) 1995 (05, 95, 98, 95) DJF 2000 MAM 1998 N polar land 1995 (, 94, 98, 95) 1995 (, 89, 98, ) DJF 2005 MAM 2000 N polar ocean 1995 (03, 05, 98, 95) 1995 (05, 95, 98, 95) MAM 2002 MAM 1998 S polar land and ocean S polar land
NA: not available.
the early 1990s in the satellite data are likely due to volcanic eruptions, which amplify cooling at altitude (Free and Lanzante, 2009). The differences in internal trends post-1996 may be due to orbital decay that has not been fully allowed for in the satellite record, with cooling from above affecting the satellite data and heating from below affecting the surface data, or a combination of these.
Unless substantially contaminated by artefacts, these changes do not represent gradual warming in the atmosphere, but they may represent regime-like change controlled from the surface. The capacity of the oceans to emit sufcient heat during El Nio events and absorb it during La Nia events to cause large warming anomalies at the global scale suggests
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 187
that available heat energy is not a limiting factor for abrupt changes.
In Fig. 4, both surface and satellite temperature records are very step-like. The trend / shift ratios for the HadCRU and GISS records are 0.19 and 0.27 respectively and for the RSS and UAH records they are 0.55 and 0.40 respectively,
showing the effect of the negative internal trends. Shifts are consequently higher than steps in the satellite data. These are clearly due to the presence of the ENSO cycle within the data, where La Nia events precede shifts and El Nio events accompany them. If they are not assumed to be a contaminating inuence of noise affecting the signal, there is no clear way to allow for them; therefore, the data are analysed and presented as they are. As we discuss later in the paper, it appears that El Nio has an active role in step-like warming.
3.3 Regional attribution
This section on regional attribution covers the issue of stationarity and the character of change over regional areas, and it addresses Test 4. Regional attribution of step changes in annual temperature has previously been carried out for southeastern Australia (SEA, Jones, 2012) and is repeated here for Texas and central England. The methodology is suitable for continental mid-latitude areas where annual average minimum temperature (Tmin) is correlated with maximum temperature (Tmin/Tmax), and Tmax is correlated with total annual rainfall (Tmax/P ) (Power et al., 1998; Nicholls et al., 2004; Karoly and Braganza, 2005). For maritime areas such as central England, diurnal temperature range (DTR) is used (DTR / P ) instead of Tmax/P . The method uses the bivariate method to test the dependent variable against the reference variable. A shift in the dependent variable denotes a regime change.
SEA climate was stationary until 1967 when a step change increased Tmin by 0.6 C with respect to Tmax (Jones, 2012).
Six independent climate model simulations for the same region become non-stationary by the same means between 1964 and 2003, showing steps of 0.4 to 0.7 C (Jones, 2012).
Texas becomes non-stationary in 1990, with an increase in Tmin/Tmax of 0.5 C. Tmax increased by 0.8 C against P in 1998. For central England, Tmin increased against DTR by0.3 C and Tmax against P by 0.9 C in 1989. Tmax also increased against P in 1911 by 0.5 C (Table 2).
The stationary period is used to established regression relationships that calculate Tmax and Tmin from P and Tmax respectively. These regressions are used to estimate how Tmax and Tmin would have evolved during the non-stationary period. The residual is then attributed to anthropogenic regional warming and is tested using the bivariate test. Here, the residuals for Tmax and Tmin are averaged to estimate externally forced warming (TavARW).
In SEA, TavARW shifted upward by 0.5 C in 1973 (Fig. 5).
Similar patterns were found for 11 climate model simulations for SEA, undergoing a series of step changes until 2100
(a)
(b)
(c)
2.0
Central England
1875 1895 1915 1935 1955 1975 1995
1.5
1.0
Anomaly (C)
0.5
0.0
0.5
1.0
1.5
2.0
Year
2.0
Texas
1895 1915 1935 1955 1975 1995
1.5
1.0
Anomaly (C)
0.5
0.0
0.5
1.0
1.5
Year
2.0
South-eastern Australia
1910 1930 1950 1970 1990 2010
1.5
1.0
Anomaly (C)
0.5
0.0
0.5
1.0
Year
Figure 5. Anomalies of annual mean temperature attributed to nonlinear changes where the inuences of interannual variability have been removed for (a) central England, (b) Texas and (c) southeastern Australia. Internal trends (dashed lines) are separated by step changes (p < 0.01).
(Jones, 2012). For Texas, TavARW shifted by 0.8 C in 1990.
Central England temperature shifted upward by 0.7 C in 1989 and by 0.5 C in 1911. Using the full record for central England average temperature from 1659, a signicant step change was found in 1920, whereas using a starting date of 1878 identies 1911. Given that the second mode identied in the longer test is 1911, we conclude that the 1911 date is an artefact of the starting date in 1878, and a step change in 1920, consistent with NH data, would register if earlier data were available.
None of the internal trends in Fig. 5 achieve p < 0.05. The trend / shift ratios for Tav (not shown in Fig. 5) and attributed to external forcing (TavARW) are 0.23 and 0.88 respectively for
SEA, 0.45 and 0.53 for Texas, and 0.01 and 0.33 for cen
tral England (18782014). The lower ratio in SEA TavARW is because reduced rainfall post-1997 produces lower attributed TmaxARW, but if that rainfall reduction is also a response to external forcing (Timbal et al., 2010), TmaxARW will be un-
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
188 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
Table 2. Year of non-stationarity in regional temperature for south-eastern Australia, Texas and central England. Data source, year of rst change greater than 1 standard deviation for Tmax against P and Tmin against Tmax, or DTR /P using the bivariate test. The stationary period is also shown.
Data source Tmax/P Tmin/Tmax DTR /P Stationary period
Year Change Year Change Year Change (SEA)
SE Australia 1999 0.7 1968 0.6 19101967 Texas 1998 0.8 1990 0.5 18951990 Central UK 1989 0.9 N/S 1989 0.3 18781988
1911 0.5
derestimated. The negative ratio for Texas is because TavARW contains negative internal trends, mostly after 1990 (largely a rainfall effect on Tmax). For central England, the ratio for Tav has been calculated from the long-term record from 1659, which shows no step changes or trends between 1701 and 1920. Late 20th century warming in both central England and the continental US elsewhere has also been analysed as nonlinear (Franzke, 2012; Capparelli et al., 2013).
These results show that the transition from stationarity to non-stationarity is abrupt for regional temperature at three locations on three continents and for six independent climate model simulations for one of those locations (SE Australia). The close association of the observed transition in SEA in 1968 with the widespread shift date over the SH mid-latitudes indicates that the onset of the warming signal in these broader regions is abrupt (Jones, 2012). The changes in central England in 1989 and Texas in 1990 may also be associated with a widespread step change in the NH mid-latitudes in 1987/88 (Overland et al., 2008; Boucharel et al., 2009; Lo and Hsu, 2010; Reid and Beaugrand, 2012; North et al., 2013; Menberg et al., 2014; Reid et al., 2016).
The low trend / shift ratios shown for ocean and some zonal areas also occur over the three land areas analysed.This suggests that shifts may be more distinct at regional scales, integrating into a more trend-like global average. This is the case for sea level rise data, where individual tide gauge records exhibit stepladder-like behaviour at individual locations and global mean sea level follows a curve (Jones et al., 2013).
3.4 Other climate variables
If climate changes in a stepwise manner, it would be expected that other variables would show signs of this (Test 5). Instances of step changes in the literature are widespread, and are mentioned elsewhere in this paper (e.g. Table 6). For rainfall, notable examples are a step change in the Sahel in 1970 (LHte et al., 2002; Mah and Paturel, 2009), south-west Western Australia (WA) in the late 1960s and early 1970s (Li et al., 2005; Power et al., 2005; Hope et al., 2010) and the western US in the 1930s (Narisma et al., 2007). Similar changes have been detected in streamow records world-
wide, showing that regime changes in moisture have been a long-standing aspect of climate variability (Whetton et al., 1990). A few more recent changes have been directly attributed to increasing gases, although south-west WA is an exception (Cai and Cowan, 2006; Timbal et al., 2006; Del-worth and Zeng, 2014), with large-scale shifts in synoptic types accompanying a rapid decrease in rainfall (Hope et al., 2006). The bivariate test identies a step change in southwest WA winter rainfall in 1969 (shown in Fig. 6a), with an upward step in summer rainfall in northern Australia 1 year later.
Ocean heat content of the upper ocean also shows step changes occurring in 1977, 1996 and 2003 (Fig. 6b). Changes in long-run tide gauge records also show a stepladder-like process of sea level rise, with the San Francisco record, quality controlled and dating back to 1855, being a good example; it shows step changes in 1866, 1935, 1957 and 1982 (Fig. 6c). Step changes in the Fremantle tide gauge data records, one of the longest in the Southern Hemisphere, shows that most of the decline in the average return intervals of extreme events noted by Church et al. (2006) before and after 1950 occurred in two events (Fig. 6d) in the late 1940s and the late 1990s. This variation in rise was noted by White et al. (2014). None of the internal trends in Fig. 6ad attain p < 0.05, showing the dynamic nature of change and limited trend-like behaviour in these examples.
4 Results models
4.1 20th century simulations (18612014)
These sections report on the multistep analysis of 102 simulations of global mean surface warming from the CMIP3 archive and 295 simulations from the CMIP5 archive. Further information on the archives can be found in the Supplement. The relevant test for models is to identify phenomena similar to observations. Here we describe analyses of the timing of change points and their relationship with known regime changes and the measurement of the relative contributions of steps, shifts and internal trends in the temperature record (covering Tests 1, 2 and 3).
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 189
SW WA w inter N Aust summer
10
(a) (b)
(c) (d)
1000
200 1900 1920 1940 1960 1980 2000
15
1955 1965 1975 1985 1995 2005 2015
900
800
Rainfall (mm)
600
Ocean Heat Content (J 1022)
10
700
5
500
0
400
5
300
Year
Year
7250
1855 1875 1895 1915 1935 1955 1975 1995 2015
6950
1910 1930 1950 1970 1990 2010
7200
6900
7150
6850
Tide guage (mm)
Tide guage (mm)
7100
6800
7050
6750
7000
6700
6950
6650
6900
6600
6850
6550
6800
6500
6750
6450
Year
Year
Figure 6. Records showing internal trends separated by step changes of (a) total rainfall for south-west Western Australia (winter) and northern Australia (summer, 19002015); (b) global ocean heat content of the top 700 m (19552014); (c) tide gauge data for San Francisco, US (18552015) and (d) Fremantle, Australia (19121925, 19272015). Step changes (p < 0.01) identied by the bivariate test.
Starting with observations, the percentages of annual steps (p < 0.01) in the 45 time series of mean annual surface temperature from Fig. 2 are shown in Fig. 7a. Two-thirds of all historical records shift in 1997 and one-third shifts in 1980 and 1937. Lesser peaks of 1015 % occur in 1920, 1921, 1926, 1930, 1968/69, 1987 and 1988. The three shifts in 1979/80, 1987/88 and 1997/98 are the main contributors to the higher rate of trend noted from around 1970. Because these peaks measure how strongly steps occur globally and regionally, percentages denote how pervasive a step is. The models only register a signicant step at the global scale, meaning they will only pick up the most extensive step changes. Any steps occurring below the assigned level of probability (p < 0.01) will show up as part of a trend, as is the case for 1987/88 in the observations.
Figure 7b shows step changes from the CMIP3-combined SRES A1B and A2 simulations for the 20th and 21st centuries: 84 are independent and 18 are ensemble averages. The CMIP3 models were driven by observed forcing, including sulfate aerosols, until 19992000 and not all contain natural forcings (see Table S2 in the Supplement). They do a reasonable job of capturing the three main post-1950 peaks.Figure 7cf show the CMIP5 RCP2.6, RCP4.5, RCP6.0 and RCP8.5 ensemble results respectively. The models were driven by observed forcing, including natural volcanic and solar forcing, until 2005. Visually, the CMIP5 results illustrate the observed peaks and troughs better than CMIP3. This is presumably due to the improved representation of forcing factors and physical processes and to improved model resolution (Table S3).
The RCP4.5 result (Fig. 7d), with 107 independent members, is the largest multi-model ensemble (MME). The three major post-1950 step changes are reproduced as follows: 55 % (58 of 107) of the runs undergo a step change in 1996 98 (17 % step in 1996, 16 % in 1997 and 22 % in 1998), 40 % of the runs peak in 197678, just missing the observed peak in 1979/80, and 19 % peak in 198688. In the mid-1970s, the models may have picked up the observed regime shift in 1976/77 in the Pacic Ocean (Ebbesmeyer et al., 1991;Miller et al., 1994; Mantua et al., 1997; Hare and Mantua, 2000) as a contemporaneous increase in warming. With weak El Nios affecting observations during 19771980 (Wolter and Timlin, 2011), this step change may have been delayed in the observed temperature record until 1979/80.
Of the pre-1950 peaks, the models peak around 1916, rather than 1920, and 1936/37 forms a minor peak, less prominent than in the observations. The volcanic eruptions of Krakatoa (1883) and Mount Agung (1963) both feature in the model simulations but less so in the observations. The mid-20th century period of little change is also reasonably well reproduced.
Correlations over the full period 18802005 between observations and the CMIP3 and CMIP5 models are 0.32 and0.34 respectively (p < 0.01). For the period 19502005, the correlations rise to 0.45 and 0.40 respectively. If specic events from 1963/64, 196870, 1976/77, 1979/80, 1987/88 and 199698 are grouped, and all other years are analysed individually, then the correlation increases to 0.78 for both CMIP3 and CMIP5 records (note that this treats the simulated and observed peaks in the 1970s separately). We con-
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
190 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
0 %
(a)
80 %
1880 1900 1920 1940 1960 1980 2000
(b)
25 %
1861 1881 1901 1921 1941 1961 1981 2001 2021 2041 2061 2081
70 %
Number of observations
20 %
60 %
40 %
Number of models
50 %
15 %
30 %
10 %
20 %
5 %
10 %
0 %
Year
Year
(c)
25 %
1861 1881 1901 1921 1941 1961 1981 2001 2021 2041 2061 2081
(d)
25 %
1861 1881 1901 1921 1941 1961 1981 2001 2021 2041 2061 2081
20 %
20 %
Number of models
Number of models
15 %
15 %
10 %
10 %
5 %
5 %
0 %
0 %
Year
Year
(e)
25 %
1861 1881 1901 1921 1941 1961 1981 2001 2021 2041 2061 2081
(f)
25 %
1861 1881 1901 1921 1941 1961 1981 2001 2021 2041 2061 2081
20 %
20 %
Number of models
Number of models
15 %
15 %
10 %
10 %
5 %
5 %
0 %
0 %
Year
Year
Figure 7. Step changes in observed and simulated surface temperatures. Frequency in percent of statistically signicant step changes from (a) global, hemispheric and zonal averages (45, 18802014); (b) global mean warming from 102 model simulations from the CMIP3 archive for SRES A1b and A2 emission scenarios; (cf) global mean warming 19612100 from the CMIP5 archive for the (c) RCP2.6 pathway (61), (d) RCP4.5 pathway (107), (e) RCP6.0 pathway (47) and (f) RCP8.5 pathway (80).
sider this a reasonable test because all these dates have been linked to regime changes or break points in temperature in the literature. Finessing the exact years involved around these events makes little difference to the result; thus, the correlation is robust.
Although collectively the model ensembles reproduce the observed peaks, single models do not fare as well. We experimented with a skill score that matched steps between models and observations, but the resulting scores did not correlate with any other factor. The only event reproduced widely by the models was the 199698 step change, peaking in 1997, when 58 of the 107 MME (55 %) underwent a step change, although 40 % of the MME produces a step in 197678.
4.2 Relationship between steps and trends over time
Here, we report on the relationships between steps, shifts, and trends; the magnitude of warming; and ECS to estimate the proportion of signal in each warming component, addressing Test 3. Total warming over time can be represented by straightforward differencing, or change measured from a simple trend and the sum of various components, such as the sum of steps and of shifts and trends. All come up with slightly different answers but describe a process that over many decades largely conforms to a trend.
Warming components measured here are steps, the internal trends between steps and the shifts from one trend to the next. Counting shifts as the remainder between internal trends preferences trends over shifts (by about 5 % in the hindcast period). When each is contrasted with an indepen-
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 191
(a)
1.75
(b)
1.8
R = 0.34
R = 0.87
R = 0.13
R = 0.00
R = 0.01
R = 0.00
Warmingcomponent18612005 (C)
Warming component 18612005 (C)
1.50
1.6
1.4
1.25
1.2
1.00
1.0
0.75
0.8
0.50
0.6
0.4
0.25
0.2
0.00
0.0
0.25
0.2
0.25 0.50 0.75 1.00 1.25 1.50
2 2.5 3 3.5 4 4.5 5
Total warming 18612005 (C)
Steps Trends Shifts
Equilibrium climate sensitivity (C)
Steps Trends Shifts
(c)
1.25
(d)
3.0
R = 0.96
1.00
2.5
0.75
Total trends (C)
2.0
0.50
18612005 Observations
Warmingcomponent20062095 (C)
1.5
0.25
1.0
0.00
0.5
0.0
0.25
0.5 1.0 1.5 2.0 2.5
0.00 0.25 0.50 0.75 1.00 1.25
Total shifts (C)
Total warming 20062095 (C)
Steps Trends Shifts
(e)
3.0
(f)
1.6
Warmingcomponent20052095 (C)
CSIRO GISS-E2-H GISS-E2-R EC-EARTH IPSL Canada CCSM Hadley MIROC MPI
R = 0.65
R = 0.18
R = 0.52
2.5
1.4
1.2
2.0
Shifts 20062095
1.0
1.5
0.8
1.0
0.6
0.4
0.5
0.2
0.0
0.0
2 2.5 3 3.5 4 4.5 5
0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75
Equilibrium climate sensitivity (C)
Steps Trends Shifts
Trends 20062095
and total internal trends correlate poorly with ECS (0.01, 0.01, 0.07 and 0.09, all NS; Table 4, Fig. 8b).
The ratio of total internal trends to total steps slightly favours shifts (mean 0.44), ranging between 0.09 and 1.22.
A low ratio means that trends either cancel each other out or are negligible. A high ratio usually indicates that the time series contains one or more negative shifts and/or a number of positive trends. Observations t comfortably within this distribution, with ratios of 0.32 to 0.38, except for the GISS time series, which has a ratio of 0.62 because of a downward shift and upward trend in the early part of the record (Fig. 8c).The MME ratios are slightly negative with respect to total warming (0.14, NS), suggesting that the mix of shifts and
trends is largely unrelated to the amount of hindcast warming (18612005).
For the historical period, total warming and its various components steps, shifts or trends are unrelated to ECS.
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
Figure 8. Multi-model ensemble (RCP4.5, 107 members) characteristics of hindcast (18612005) and projection (20062095) periods. (a) Relationship between total warming and steps, trends and shifts (18612005); (b) relationship between ECS and steps, trends and shifts (18612005); (c) total shifts and total trends for 19612005 with observed points from ve warming records; (d) relationship between total warming and steps, trends and shifts (20062095); (e) relationship between ECS and steps, trends and shifts (20062095); (f) total shifts and total trends for 20052095 from individual climate models.
dent variable such as ECS, this poses a strong test for shifts because internal trends estimate Hstep in each time series.
The hindcast (18612005) and projection (20062095) components of the RCP4.5 107-member ensemble were analysed separately.
For the hindcasts (18612005), total warming (the 2000 05 average minus the 186199 average) is positively correlated with total steps (0.93, p < 0.01). Their means are 0.97 and 0.94 C. The correlation between total warming and internal trends is 0.36 (p < 0.01) and is 0.58 between total warming and shifts (p < 0.01). Shifts therefore explain 2.5 times the variance explained by internal trends in estimating total warming (Fig. 8a). A simple linear trend measured over the entire period has the same correlation with steps (0.93, p < 0.01) but averages 0.76 C, thus underestimating total warming by 0.18 C. Total warming, total steps, total shifts
192 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
Table 3. Steps collated for each decade from 1876 to 2195 from the RCP4.5 MME, showing total steps up and down and the correlation between step size and ECS. The second part of the table shows the correlations between total warming, steps and trends over the observed and simulated periods and ECS. Correlations are classied as not signicant (NS, p > 0.05), p < 0.05 ( ) and p < 0.01 ( ). Total correlations with the MME are n = 107 and with ECS are n = 92.
Change and period Steps up Steps down Correlation with ECS Signicance
Steps 18761885 0 26 0.40
Steps 18861895 13 1 0.32 NS
Steps 18961905 7 1 0.09 NS
Steps 19061915 31 0 0.27 NS Steps 19161925 65 0 0.27
Steps 19261935 17 1 0.09 NS Steps 19361945 33 0 0.20 NS Steps 19461955 6 1 0.85
Steps 19561965 4 12 0.52
Steps 19661975 29 0 0.33 NS Steps 19761985 56 0 0.41
Steps 19861995 34 0 0.49
Steps 19962005 101 0 0.19 NS
Steps 20062015 83 0 0.68
Steps 20162025 82 0 0.65
Steps 20262035 70 0 0.74
Steps 20362045 82 0 0.66
Steps 20452055 75 0 0.57
Steps 20562065 65 0 0.67
Steps 20662075 61 0 0.60
Steps 20762085 51 0 0.66
Steps 20862095 27 0 0.82
Mean ( C) Range ( C)
Warming 18612005 0.9 0.41.4 0.01 NS
Warming 20062095 1.5 0.72.4 0.81
Steps 18612005 1.0 0.31.5 0.01 NS
Steps 20062095 1.6 0.72.5 0.81
Shifts 18612005 0.6 0.01.2 0.07 NS
Shifts 20062095 0.8 0.31.5 0.72
Trends 18612005 0.4 0.01.0 0.09 NS
Trends 20062095 0.8 0.11.6 0.43
The relationship between total shifts and total internal trends is negative (0.47, p < 0.01), which is to be expected, but
the lack of a relationship between the shift / trend ratios and warming or ECS suggests that this uncertainty is stochastic.
For the projection period, total warming over 200695 is based on the difference between 5-year averages centred on 2006 and 2095. Total warming averages 1.55 C, total steps average 1.57 C and they are highly correlated (0.98, p < 0.01). The correlation between shifts and internal trends with total warming is 0.70 and 0.74 respectively, with trends having a slightly higher correlation (Fig. 8d). However, correlations between ECS and total steps, shifts and trends are0.81, 0.72 and 0.43 respectively (all p < 0.01, Fig. 8e). This shows that the time series are becoming more trend-like at higher rates of forcing, when compared to the hindcast period. Shifts have 2.9 times more explanatory power than trends with respect to ECS, but 0.9 times the explanatory
power with respect to total warming over 20062095. We take this to mean that shifts (steps minus internal trends) carry most of the signal and that trends are more random since they are affected by short-term (interannual) stochastic behaviour. Some of the signal embedded in trends could also be due to shifts occurring at regional scales, which are too small to register statistically as steps at the global scale.
The ratio of trends to steps is 0.51, ranging from 0.14 to0.88. The ratio of trends to shifts favours trend (1.22) but has a large range (3.25 to 0.15). The correlations of both ratios with warming are very low (0.07 and 0.03 respectively, NS). This paradox, where there is no correlation with the amount of warming but there is with ECS, when both ECS and warming are correlated, can be viewed by plotting the different modelling groups according to the relationship between shifts and trends. Individual models plot along linear pathways, as was the case for the hindcast ensemble (Fig. 8f).
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 193
(a)
Figure 9. Correlations between ECS and linear trends, total step changes, warming to date and quadratic trends (a) from 1861 to the current decade (warming to date: 186199 average subtracted from current decadal average) and (b) from 1961. Dotted lines mark p < 0.01.
The high-sensitivity models plot towards the upper right and lower sensitivity models plot towards the lower left. The trend / step ratios for these individual groups vary widely; the CSIRO eight-model ensemble has ratios from 0.25 to 0.56 and the GISS-E2-R 17-member ensemble ranges from 0.17 to 0.72. The potential for the same model to produce very different shift / trend ratios shows high stochastic uncertainty, probably generated by oceanatmosphere interactions. The timing of these interactions appears to be largely unrelated to climate sensitivity, although the warming response to steps when they do occur is related to sensitivity.
Interestingly, the GISS models form two groups, the main difference being the ocean conguration (see Schmidt et al., 2014a), where the Russell ocean model produces more step-like outcomes and the HYCOM ocean model produces more trend-like outcomes.
For each individual decade from 18761875 to 2086 2095, correlations were performed between step size and ECS (Table 3). The late 19th century produces downward steps in response to the Krakatoa eruption in 1883 and is negatively correlated with ECS. Positive steps dominate from 1886 through to 1945 and are positively correlated at levels of low or no signicance. The period 1946 to 1965 is negatively correlated with ECS; in 195665, corresponding with the 1963 Mount Agung eruption, downward steps result in a negative correlation of 0.52 (p < 0.05). Correlations
between ECS and step size become positive after 1965, being 0.41 for 197685 and 0.49 for 198695 (both p < 0.01). For the decade 19962005, 101 of the 107 members of the MME underwent an upward step, but the correlation with ECS is only 0.19 (NS). This low correlation may partly be due to a rebound from the negative forcing of the 1991 Mount Pinatubo eruption in the models, which has been overestimated by about one-third (Schmidt et al., 2014b). Correlations for the forcing period (20062095) rose to 0.68 in 200615 and vary between 0.57 and 0.82 for subsequent decades to 2095.
The lack of predictability in the hindcasts is a result of negative aerosol forcing due to volcanic eruptions and anthropogenic sources occurring after 1950. The more-sensitive
models produce strong positive and negative responses depending on the direction of forcing, whereas in the less-sensitive models this effect is reduced. This effect cancels out any consistent relationship between ECS and step size over the historical period. The implication of this nding is that the magnitude of 20th century warming in the models has little predictive skill and is not a reliable guide to potential future risk.
The hindcast results are also uncorrelated with the 21st century projections. Total warming (18612005) is negatively correlated with 21st century warming (200695,
0.25, p 0.01) and is uncorrelated with respect to ECS
(0.01). Total steps from the hindcast and forecast peri
ods show similar negative correlations. Internal trends 1861 2005 are also uncorrelated with future total warming, steps or trends. This strongly indicates that 20th century warming may not be a good guide to future warming, if observations are being affected in a similar way.
A nal analysis looks at the explanatory power of different change models with respect to ECS over time. Linear and quadratic trends, steps and warming to date are calculated for successive decades for each ensemble member and the results correlated with ECS. Both trends and warming difference respond to negative forcing in the rst part of the record. Step changes are less volatile, remaining close to zero until increasing from 1995 and remaining higher than the other models until the end of the century (Fig. 9a). The standard error measured from total accrued warming was also least out of the three statistical models. Although it would be possible to derive a closer t for some of those models with a greater number of factors, step changes clearly carry the greatest signal with respect to ECS over time. The analysis repeated from 1965 produces a similar result (Fig. 9b).
This result is further evidence that step changes carry the signal. Warming to date assesses any warming irrespective of its cause, whereas if step changes are part of a direct response to forcing, they would be a better predictor. This is the case for climate models, and may therefore apply to observations as well. The advantage for using warming to date as a measure is that it has roughly a decades advantage over sta-
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
1.00
0.75
0.50
0.25
0.00
0.25
0.50
(b)
1.00
1965 1985 2005 2025 2045 2065 2085
0.75
Correlation coefficient
Correlation coefficient
0.50
0.25
0.00
0.25
0.50
1875 1895 1915 1935 1955 1975 1995 2015 2035 2055 2075 2095
Linear trend Step changes Warming to date Quadratic trend
Step changes Warming to date Quadratic trend
194 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
Figure 10. Global mean surface temperature as analysed by the multistep bivariate test. (a) Step and trend breakdown of global mean surface temperature in the RCP2.6, 4.5, 6.5 and 8.0 simulations from the HadGEM-ES model, run 3; (be) Ti0 results from a 40-year moving window for the RCP2.6, 4.5, 6.5 and 8.0 simulations respectively.
tistical tests, which require hindsight. Therefore, unless the physical mechanism(s) for steps become known, both have roughly equivalent predictive skill at the present time.
4.3 21st century forcing proles
If increased forcing raises the rate of entropy production, we would expect to see step-like behaviour becoming more trend-like over time. Such behaviour would involve either
an increase in the frequency and distribution of regional step changes that would integrate to become more trend-like at the global scale, or
we would see an increase in the rate of diffuse warming, producing widespread trend-like behaviour.
If either is the case, then simulations for the four different emission pathways, RCP2.6, 4.5, 6.0 and 8.5, should show this.
Figure 7cf show the percentage of step changes in any given year for the multi-model ensemble for each of these pathways. For RCP2.6, peaks occur until about 2050, after which the ensemble stabilises. Some models step downward, the earliest of which occurs in 2051. Individual members stabilise between 2018 and 2092, with 48 of the nal shifts being positive and 13 negative. This timing is weakly correlated
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
haviour under increasing emissions, RCP2.6 shows a return to step-like changes, stabilising around 2050; RCP4.5 shows a return to step-like change late century; RCP6.0 shows increasing trend-like behaviour over the century and RCP8.5 shows a consistent trend until the end of the century, with few steps.
An indication of change at the regional scale and how it may relate to global change is illustrated by using selected CMIP3 models for south-east Australia, as described in Jones (2012). For example, for the CSIRO Mark3.5 A1B simulation, for global mean warming, internal trends comprise 52 % of total warming from 2006 to 2095, whereas for SEA Tmax the ratio is 13 % and Tmin 47 %. These were consistent for A1B- and A2-forced simulations, which are roughly equivalent to RCP4.5 and 6.0. The number of step changes is also notable: four and ve at the local scale and 12 at the global scale (Fig. 11). The higher ratio for Tmin compared to Tmax may be due to Tmin being related to large-scale sea surface temperature patterns and Tmax being related to more local soil moisture patterns, as is the case for the central and western US (Alfaro et al., 2006). Jones et al. (2013) showed that such changes at the local scale produce signi-cant increases in impact risks.
These analyses do not support increasing trend-like behaviour at the local scale, and therefore favour the rst alternative above, but further work across more regions is required to conrm this.
5 Testing of steps versus trends
Earlier sections have identied steps and trends in temperature and tested how trends, steps and trendshift relationships relate to total warming and the independent variable ECS. This section examines how well trend, step and step trend models reproduce the temperature records examined throughout the paper. This tests htrend against hstep. The error value assigning p < h0 is not the principal measure being sought. Instead, the statistical model that combines low error with unstructured residuals while sustaining physically plausible assumptions is preferred. Another aim is, if possible, to provide likelihoods for severe testing.
Four statistical models are tested: ordinary least square trend, LOWESS, step, and step and trend. The LOWESS model (locally weighted regression; Cleveland and Devlin, 1988) was applied with a bandwidth of 0.5 to assess sensitivity to uctuations in the data, contrasting those with both the trend and step models. It is not considered a valid statistical rival because it is tted without regard to physical process.Likewise, although the step-and-trend model will t well to the data, the step model is the one used for severe testing, being a straightforward measure of hstep. The trend model represents htrend.
With the data produced, we look at goodness of t (r2), the residual sum of squares (ResSS), cumulative residuals (
PR)
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 195
with ECS (0.18, NS). ECS is uncorrelated with the size of the nal shift or with the gradient of the following trend. The RCP4.5 ensemble produces frequent steps that peak around 2025 and decline towards the end of the century. RCP6 produces a fairly constant rate of steps and RCP8.5 produces sustained steps throughout the century, peaking in the 2080s at a higher rate than 199698.
This evolution shows a stepladder-like process in the 20th century that changes into an elevator-like process in the 21st century, becoming more trend-like with increasing forcing.Depending on the subsequent rate of forcing, trend-like processes can either recede back to a step-like process or even stabilise. The HadGEM2-ES single model ensemble is used to illustrate this (Fig. 10a).
This ensemble shares the same historical forcing until 2005. It warms by less than observations until 2010, with a reversal in 19641980; then it warms substantially in a series of steps over the next few decades. It undergoes a step change of 0.37 C and shift of 0.18 C in 1998, 1 year after the ob-served shift. The next steps occur in 2012, 2013, 2014 and 2015 in the four simulations, ranging from 0.40 to 0.49 C in absolute terms and 0.19 to 0.27 C as the shift from the pre-step trend to the post-step trend. The rst half of the 21st century shows the inuence of decadal variability on mediating step changes. In 2021, the RCP2.6 simulation undergoes a step change and is higher than the others for most of that decade. The RCP6.0 simulation is lower than the others from 2025 to 2045 before accelerating under a sustained step-and-trend process. The relative proportion of internal trends to total warming under the four scenarios is 0.34, 0.60, 0.57 and 0.79 for warming of 1.9, 2.9, 3.7 and 5.3 C respectively.
The RCP4.5 has a higher trend ratio, showing the stochastic uncertainty inherent in the simulations.
Like most statistical tests that detect change points, the bivariate test is considerably weakened under autocorrelated data, where its timing is fairly robust but p(H0) becomes increasingly sensitive. Such autocorrelations may be caused by simple trends, with lag-1 or longer lag processes inuencing the complex nature of warming. Removing these without assuming an underlying process is difcult. Thus, one way of assessing its inuence is to pass a moving window through a time series. If the data are step-like and largely free of autocorrelation, a distinct step will produce a line of horizontal Ti0 statistics on a single date as it passes through the window.
If there are no steps within a window period and autocorrelation is low, background Ti0 values will return to low values (single digits). With autocorrelation, background Ti0 values remain above the p < 0.01 threshold and form a cloud, rather than steps producing horizontal lines.
In Fig. 10be, successive horizontal lines extending right from low Ti0 values indicate stepladder-like behaviour in the 20th century. Horizontal lines that stay on the right without returning to low Ti0 values indicate both step-like and trending behaviour. A cloud to the far right, as in Fig. 10e, shows a trend-dominated process. Summarising 21st century be-
196 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
2.5
(a) (b)
(c)
5.5
1870 1920 1970 2020 2070
4.5
1870 1920 1970 2020 2070
5.0
4.0
4.5
4.0
3.5
3.5
3.0
Anomaly (C)
Anomaly (C)
3.0
2.5
2.5
2.0
2.0
1.5
1.0
1.5
0.5
1.0
0.0
0.5
0.5
1.0
0.0
1.5
2.0
0.5
1.0
Year
Year
4.0
1870 1920 1970 2020 2070
3.5
3.0
2.5
Anomaly (C)
2.0
1.5
1.0
0.5
0.0
0.5
Year
Figure 11. Anomalies of annual mean temperature showing internal trends separated by step changes from the CSIRO Mk3.5 A1B simulation: (a) Maximum temperature of south-eastern Australia; (b) minimum temperature of south-eastern Australia and (c) global mean surface temperature. Internal trends (dashed lines) are separated by step changes (p < 0.01).
and cumulative residuals squared (
PR2 the trend model behaves more poorly than the other three.
The LOWESS test performs less well than the autocorrelation and heteroscedasticity tests for the 40-year windows. Although the LOWESS model performs well over the whole record, it is subject to deviations within the record that cancel each other out akin to cutting corners. The step and trend model performs worst for F -hetero over the whole record, but the best over 40-year windows. This is due to high variance within the early part of the record and is an issue of precision, as standard error of this relationship is almost half that of the trend model (not shown, but is similar to the
PR2). Residuals (R) show how much variance is explained by the model, cumulative residuals will show whether residuals are showing structure not explained by the model and cumulative residuals squared show accumulating error, including rapid changes not accounted for. Four more tests have been added to these: F tests for autocorrelation (F -auto) and heteroscedasticity (F -hetero) of the residuals over the whole record and percentage of exceedance over moving 40-year windows. Whites test (White, 1980) is used for heteroscedasticity. The rst four of these tests use absolute error, or the amount of a time series not explained by the statistical models, and the second four show patterns, working on accuracy and precision. The statistical models that fail to combine both are therefore the weakest.
Results are shown in Fig. 12 and Tables 4 and 5. The data and statistical models for the HadCRU record for 18802014 are shown in Fig. 12a. Cumulative residuals that track close to zero (Fig. 12b) show the model mimicking the data closely and sustained departures show signicant deviation. Here, the trend model deviates substantially and the LOWESS model less so, while the step and step-and-trend models deviate least. This follows through to the cumulative residuals squared. The less change the better, whereas upward kinks show rapid changes or large outliers (positive or negative) not incorporated into the model (Fig. 12c). Trend analysis produces an r2 value of 0.76 and residual sum of squares of0.87, and the other three statistical models have an r2 of 0.87
and ResSS of 0.8. For
PR2 relationship). The step model is clearly superior to the trend model for the moving window tests. The results for the other four long-term global warming records, BEST, C&W, GISS and NCDC, are not shown but have similar results.
These tests, omitting LOWESS, were carried out for Had-CRU 19652014, a period with a sustained radiative forcing signal (Fig. 12d). The results for the different statistical models are similar, with r2 values of 0.85, 0.86 and 0.89 respectively. The step-and-trend model is still the best performed, but the step model is only slightly better than the trend model.This is due to the NH shift in 1987/88 being incorporated into the global mean trend. Dividing this time series into quarters will bring 1987/88 into the picture but will also make both the MYBT and Students t tests more sensitive.
Also shown in Table 4 are the zonal temperatures from NCDC 3060 N (18802014), where total internal trends
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 197
0.4
1.0
1880 1900 1920 1940 1960 1980 2000
1.0
1965 1975 1985 1995 2005
0.8
0.8
0.6
Anomaly (C)
Anomaly (C)
0.4
0.6
0.2
0.4
0.0
0.2
0.2
0.0
(a)
Anomaly Linear LOWESS Step Step trend
(d)
Anomaly Linear Step Step trend
4
1880 1900 1920 1940 1960 1980 2000
0.6
1965 1975 1985 1995 2005
3
0.4
Cumulative residuals
Cumulative residuals
2
0.2
1
0.0
0
0.2
1
0.4
2
0.6
3
0.8
(b)
Linear LOWESS Step Step trend
(e)
Linear Step Step trend
3.0
1880 1900 1920 1940 1960 1980 2000
0.5
1965 1975 1985 1995 2005
Cumulative residuals squared
Cumulative residuals squared
2.5
0.4
2.0
0.3
1.5
0.2
1.0
0.5
0.1
0.0
0.0
(c)
Linear LOWESS Step Step trend
(f)
Linear Step Step trend
Figure 12. Testing three models to mean global anomalies of surface temperature from the HadCRU record for 18802014 (ac) and 1965 2014 (df); (a) and (d) mean annual anomalies and linear step change and shift-and-trend models; (b) and (e) show cumulative residuals for each model, where success is measured as tracking close to zero; (c) and (f) show the cumulative sum of residuals squared, where upward steps show non-linearity not explained by each model.
are slightly negative (0.04 C) and shifts are positive
(1.13 C or 106 % of steps). The pattern of results is similar to those for the global HadCRU record but the residuals are slightly more than double and the cumulative residuals are almost double, showing the step-like structure of this record.Here, the step model is clearly superior to the trend model, which fails Whites test for the whole record, fails the 40-year F -auto at a level of 51 % and has a ResSS double that for steps. This record is entirely made up of steps, showing the lack of trend occurring within some regions.
The quarterly record of HadCRU from Fig. 4 (19652014) is more ne-grained, incorporating the 1987/88 shift (Table 4). If warming is gradual, the results for trends should be
scalable; however, they perform less well at this timescale.
The respective r2 results are 0.69, 0.72, 0.75 and 0.76, whereas the differences in the cumulative residuals are 2.0,0.5, 0.7 and 0.2, where zero is a perfect score. Here, the LOWESS model performs similarly to the step model because it closely follows the data. The step model performs better than the trend model for HadCRU quarterly data, and almost as well as the step-and-trend model. For the GISS quarterly data, the results are similar.
The satellite records are more step-like than surface temperature when measured using cumulative residuals. The step-and-trend model for the 40-step window heteroscedasticity tests for satellite data fails for both RSS and UAH. This
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
198 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
Table 4. Results of eight tests on four statistical models for selected observed global temperature data (except where noted). The statistical models tested are trends (power shown), LOWESS (0.5 total series smoothing), steps, and steps and trends. Results include the adjusted r2 value, the residual sum of squares (SS), cumulative residuals and squared cumulative residuals. F tests for the whole series are shown, with p < 0.05 and p < 0.01 noted if registered, otherwise p > 0.05. F test failure for 40-year period autocorrelation and heteroscedasticity is measured at p < 0.01.
Cumulative Cumulative F test auto- F test hetero- 40-year periods 40-year periods
residuals residuals2 correlation scedasticity fail F test fail F test
Model r2 Residual SS (PR yr1) ( PR2 yr1) (F , pH0) (F , pH0) autocorrelation heteroscedasticity
HadCRU 18612014
Trend 0.76 2.6 1.2 1.3 0.0 3.7 58 % 13 %
LOWESS 0.87 1.4 0.7 0.8 0.3 1.0 28 % 13 %
Step 0.87 1.4 0.5 0.8 0.7 3.2 0 % 0 %
Steptrend 0.87 1.3 0.1 0.8 0.2 5.8, 0.05 0 % 0 % HadCRU 19652014Trend 0.85 0.43 0.20 0.24 0.0 1.2 0 % 0 %
Step 0.86 0.40 0.20 0.21 0.4 0.7 0 % 0 %
Steptrend 0.89 0.31 0.06 0.18 0.0 1.4 0 % 0 % NCDC 3060 N 18802014
Trend 0.64 6.3 1.8 2.3 0.0 10.2, 0.01 51 % 9 %
LOWESS 0.79 3.7 0.9 1.6 0.2 3.0 19 % 0 %
Step 0.83 2.9 0.3 1.4 0.0 3.0 0 % 1 %
Steptrend 0.83 2.9 0.2 1.4 0.0 3.2, 0.05 1 % 0 % HadCRU quarterly 19792014Trend 0.69 1.7 2.0 3.5 0.0 1.1 20 % 3 %
LOWESS 0.72 1.6 0.5 3.3 0.2 2.8 3 % 5 %
Step 0.75 1.4 0.7 2.8 0.0 0.2 0 % 0 %
Steptrend 0.76 1.3 0.2 2.7 0.0 0.4 0 % 4 % GISS quarterly 19792014Trend 0.67 1.9 1.6 4.1 0.0 1.1 20 % 0 %
LOWESS 0.69 1.8 0.5 3.9 0.1 2.2 6 % 2 %
Step 0.71 1.6 0.9 3.4 0.0 0.0 4 % 0 %
Steptrend 0.72 1.6 0.3 3.3 0.0 0.6 0 % 0 % RSS quarterly 19792014Trend 0.40 3.4 4.4 6.9 0.0 1.2 11 % 6 %
LOWESS 0.46 3.1 1.1 6.4 0.3 2.3 4 % 14 %
Step 0.52 2.7 0.9 5.5 0.0 0.3 4 % 8 %
Steptrend 0.53 2.6 0.7 5.1 0.0 1.3 0 % 37 % UAH quarterly 19792014Trend 0.35 3.6 3.1 7.4 0.0 1.8 6 % 9 %
LOWESS 0.39 3.4 1.0 7.2 0.1 3.3, 0.05 4 % 20 %
Step 0.46 3.0 1.5 6.1 0.0 0.7 7 % 12 %
Steptrend 0.46 2.9 0.8 5.8 0.0 1.5 4 % 42 %
is due to two instances of short-term departures on an otherwise stable background that measures heteroscedasticity as signicant with the F test: (1) a warm period during 1998, which is represented as a single step but lasts four quarters, and (2) a small warming event associated with an El Nio event in 2010 lasting two quarters. Removing this short-term warming from these sequences removes the heteroscedastic-
ity. Therefore, although not all deviations are removed by representing the satellite record as stepwise, it still provides a better explanation of change than the trend model.
Simulated global annual mean surface temperatures from climate models show results consistent with observations (Table 5). The data from Fig. 10 were analysed in the same way, except that quadratic (RCP4.5, RCP6.0), cubic
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 199
Table 5. Results of eight tests on four statistical models for representing global mean warming from HadGEM-ES climate model run 3 RCP2.6, 4.5, 6.0 and 8.5, showing the amount of warming for different measures. The statistical models tested are trends (power shown), LOWESS (0.5 total series smoothing), steps, and steps and trends. Results include the adjusted r2 value, the residual sum of squares (SS), cumulative residuals and squared cumulative residuals. F tests for the whole series are shown, with p < 0.05 and p < 0.01 noted if registered, otherwise p > 0.05. F test failure for 40-year period autocorrelation and heteroscedasticity is measured at p < 0.01.
Warming Steps Trends Shifts Pathway ( C) ( C) ( C) ( C)
RCP2.6 1.93 2.29 0.65 1.24 RCP4.5 2.93 3.30 1.76 1.07 RCP6.0 3.65 3.86 2.09 1.75 RCP8.5 5.34 5.35 4.24 1.41
Cumulative Cumulative F test auto- F test hetero- 40-year periods 40-year periods residual residual2 autocorrelation heteroscedasticity fail F test fail F test
Model r2 Residual SS (
PR yr1) ( PR2 yr1) (F , pH0) (F , pH0) correlation scedasticity
RCP2.6
Trend (x4) 0.95 3.9 4.7 3.6 0.4 8.9, 0.01 75 % 18 % LOWESS 0.96 4.7 7.7 2.8 6.9, 0.01 0.4 64 % 31 % Step 0.98 1.1 0.04 1.2 0.1 10.7, 0.01 1 % 3 % Steptrend 0.98 0.9 0.01 1.1 0.0 12.1, 0.01 0 % 4 %
RCP4.5
Trend (x2) 0.95 8.8 16.6 4.8 0.8 2.1 77 % 73 % LOWESS 0.99 3.9 13.3 2.5 2.3 4.1, 0.05 61 % 45 % Step 0.98 2.4 0.5 1.4 0.0 5.7, 0.05 19 % 14 % Steptrend 0.99 1.0 0.02 1.1 0.0 13.4, 0.01 0 % 2 %
RCP6.0
Trend (x2) 0.97 4.5 51.1 5.2 3.7 23.5, 0.01 63 % 56 % LOWESS 0.98 2.9 24.6 2.4 0.9 8.3, 0.01 52 % 31 % Step 0.99 1.2 0.06 1.2 0.1 9.7, 0.01 2 % 5 % Steptrend 0.99 0.6 0.01 1.1 0.0 17.9, 0.01 0 % 20 %
RCP8.5
Trend (x3) 0.99 4.3 4.5 3.1 0.0 11.8, 0.01 62 % 39 % LOWESS 0.992 3.1 66.6 2.8 2.0 4.5, 0.05 45 % 22 % Step 0.99 8.1 2.0 1.7 0.2 106.7, 0.01 13 % 18 % Steptrend 0.997 0.7 0.01 1.1 0.0 12.0, 0.01 0 % 3 %
(RCP8.5) and quartic (RCP2.6) polynomial functions were used instead of a linear trend. The LOWESS model used here at 0.5 record length is relatively low resolution, providing 120-year smoothing. The step model outperforms both the trend and the LOWESS model in all simulations, with the exception of the ResSS in the RCP8.5 simulation. The RCP2.6 simulation is the most step-like. In the RCP4.5 simulation, the step model does slightly worse than in the RCP6.0 simulation, which is actually more step-like. This shows the role of stochastic uncertainty in the warming process, as portrayed in Fig. 8f. The RCP8.5 simulation is the most trend-like; the step model fails in the nal decades of the 21st century because the bivariate test detects no steps, but the climate continues to warm. This is what we would expect if shifts became more local and more frequent, integrating into a curve at the global level, much like sea level rise does today.
6 Severe testing summary
A range of statistical tests have been used to examine hstep and htrend as representatives of scientic hypotheses H1 and H2. The focus is on whether atmospheric warming is gradual, forming a monotonic or even segmented trend or is stepwise and periodic, forming a complex trend over time.
As stated in the introductory sections, no single test can undertake that task. We rely on the multistep MaronnaYohai bivariate test to identify step changes in the input data, but we make as few assumptions as possible. A total of six tests with links to the two substantive hypotheses were proposed earlier in the paper. These are designed to pinpoint discrepancies between H1 and H2 by analysing the temperature data they seek to explain. The data generated consist of steps, trends and shifts calculated using the multistep MYBT model and least square trend analysis. The use of statistical models such as LOWESS is for sensitivity testing and not part of the probative assessment.
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
200 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
The test results are summarised through the following ndings:
Test 1 What patterns of step changes can be detected in temperature observations?
Global and regional analyses of steps show a highly coherent pattern of change points, where warming in the second half of the 20th century aligns with known regime changes associated with changes in decadal variability (Table 6). These events comprise the major proportion of historical warming until 2014.
Analysis of steps, internal trends and shifts in observations attributes higher proportions of warming to shifts at the zonal scale (up to 100 %), moving to lower proportions at the global scale. Three regional assessments also contain high shift / step ratios, with trends playing a lesser role.
This effect is larger in mid-latitude regions and with SST, indicating the role of equator-to-pole hydrothermal transport of energy in the ocean atmosphere system. Their timing shows that a strong role is being played by decadal variability.
Surface and satellite temperatures undergo contemporaneous shifts at the global scale, largely removing the discrepancy between trends within the two data sets. Both surface and satellite temperature records are very step-like, with surface trend / shift ratios of 0.19 and 0.27 and satellite ratios of 0.55
and 0.40 showing the effect of downward internal
trends. Shifts are consequently higher than steps in the satellite data.
Test 2 Do models reproduce the patterns of steps changes shown in observations?
Correlations between step change frequency in the observed 44-member group of global and regional data and the CMIP3 and CMIP5 MMEs analysed (18802005) are 0.32 and 0.34 respectively (p < 0.01). For the period 19502005, correlations rise to 0.45 and 0.40 respectively. Grouping specic events (1963/64, 196870, 1976/77, 1979/80, 1987/88 and 199698) and analysing other years individually, correlation increases to 0.78 for both CMIP3 and CMIP5 records. Variations in forcing, especially from volcanoes, may affect the timing and direction of step changes, but they are not the sole cause, given that 21st century simulations produce step changes from smoothly varying changes in forcing.
Of a 107-member MME (CMIP5 RCP4.5), 58 members show a step change in 199698, reproducing the observed change in 1997 within 1 year.
Test 3 What is the relationship between different components of change?
For simulated historical warming during 1861 2005, the r2 values for steps, shifts and trends in explaining total warming are 0.87, 0.43 and 0.13 respectively. Simulated warming for this period is not correlated with ECS.
For the 21st century (20062095), the r2 values for steps, shifts and trends in explaining total warming are 0.96, 0.54 and 0.49 respectively. The r2 values for steps, shifts and trends in explaining ECS are0.65, 0.52 and 0.18 respectively.
Test 4 Can step-like change be identied using attribution methods?
In all three locations on three continents tested, and for six independent climate model simulations for south-eastern Australia, warming commenced with a step change in Tmin and sometimes Tmax. Warming is not slowly emergent in any of these data as would be expected if it were gradual. The coincident timing of shifts in south-eastern Australia with southern hemispheric step changes and those in the UK and US with northern hemispheric changes, suggest that warming has commenced abruptly in different areas of the globe at different times and that the separation between stationarity and nonstationarity in the temperature record is abrupt.
Test 5 Do other climate variables also undergo step changes?
Step changes exhibiting similar timing have been shown for tide gauge observations, rainfall, ocean heat content, forest re danger index and a range of other climate variables, in addition to many impact variables (Jones et al., 2013). These are overwhelmingly attributed to random climate variability, including abrupt changes identied as part of decadal regime change.
Test 6 Are temperature time series more step-like or trend-like?
For observations and selected model data, the simple stepladder model performs better than the monotonic trend model for goodness of t (r2), the residual sum of squares (ResSS), cumulative (
PR) residuals and cumulative residuals squared (
PR2), Whites test for heteroscedasticity, a moving 40-year window regression of the residuals and a moving 40-year window of Whites test.
Table 6 summarises the major tests undertaken with expected outcomes for htrend and hstep. While objections
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 201
Table 6. Selected test results that distinguish between htrend and hstep. The null positions for each are generally not considered diametric. There is no generally accepted null with respect to htrend that references nonlinear change, whereas for Hstep the null has no signicant stepwise change points, or if it does, they are completely random and do not contain an external forcing signal.
Test results Evidence htrend hstep Supporting literature
Global warming 1895 2014
Trend / step ratio 0.320.38 (4 records),0.58 (1 record)
Trend / shift ratio 0.440.58 (4 records), 1.38 (1 record)
Gradual change, uctuations but no steps
Substantial fraction of record contains steps
Varotsos et al. (2014); Belolipetsky et al. (2015); Bartsev et al. (2016)
Overland et al. (2008); Chikamoto et al. (2012a, b); Reid and Beau-grand (2012); Menberg et al. (2014)
Regime changes 1997
29 in 1997, 37 in 199698 of 45 global and regional records
Extreme El Nio 1997/98, stochastic event
Stepwise change points identied in temp. and physically related records
Overland et al. (2008); Boucharel et al. (2009); Lo and Hsu (2010); Reid and Beaugrand (2012); North et al. (2013); Menberg et al. (2014); Reid et al. (2016)
1987/886 in 1987, 4 in 1988 of 44 regional records. Global ocean NH, NH mid-lat
El Nio, stochastic event Stepwise change points identied in temp. and physically related records
197915 in 1979, 7 in 1980, 5 in 1977, 1 in 1976 of 44 global and regional records. Global, tropics, SH
N Pacic regime shift 1976/77, El Nio 1978/79
Stepwise change points identied in temp. and physically related records
Hare and Mantua (2000); Overland et al. (2008); Meehl et al. (2009); Fischer et al. (2012); Reid and Beaugrand (2012); Menberg et al. (2014)
19694 in 1969, 8 in 196870, Southern Hemisphere
El Nio, stochastic event Stepwise change points identied in temp. and physically related records
Li et al. (2005); Hope et al. (2010); Jones (2012)
Scalability of regional records
Records more step-like at zonal and regional scales and over the oceans.
Regional records would be trend-like if warming were diffuse and gradual
Regional records more step-like, large-scale records more trend-like.
None located
Attribution Stepwise attribution forSE Australia (obs. and models), Texas (obs.),central England (obs.)
Gradual emergence of signal Abrupt emergence of signal Jones (2012)
None located
Quarterly surface and satellite temperature 19792014
Surface and satellite records share similar shifts but not trends
Signicant trend for periods > 30 years
Contemporaneous stepwise change points in independently measured records
Simulated temperature patterns 18612005
Clustering on runs test highly non-random (p 0.0 run test)
Signicant correlations between timing of steps in models and obs. CMIP30.32, CMIP5 0.34 18802005.
No matching patterns, randomicity
Matching stepwise changes between models and observations
None located
Simulated temperature quantities 18612005
Trend / step ratio 0.44 0.22 Gradual change, deviations but
no steps
Substantial fraction of record contains shifts
None located
None located
Simulated temperature relationships with independent variable ECS RCP4.5 20062095
Correlation and r2 between ECS and total warming 0.81 and 0.65, steps 0.81 and 0.65, shifts 0.72 and 0.52, and internal trends 0.43 and 0.18
Shifts random with respect to forcing
Shifts and steps more highly correlated with ECS and warming than trends
None located
Autocorrelation and heteroscedasticity observations18802014
Steps perform better than simple trends (failure rate trends 58 1 % autoc.,
10 4 % heterosc.; steps 2 4 % au
toc., 0 % heterosc. 40-year window)
Trends serially independent data, variations due to independent processes
Steps perform better than trends to explain autocorrelation and heteroscedasticity
None located
Autocorrelation and heteroscedasticity observations19652014
Trends and steps pass all tests for annual data, steps slightly better correlated than trends (0.86, 0.85 HadCRU)
Trends serially independent data, variations due to independent processes
Steps perform better than trends to explain autocorrelation and heteroscedasticity
None located
Autocorrelation and heteroscedasticity quarterly observations surface temp 19792014
Trends fail 40-year autoc. 20 %, steps 0 %, accumulated error trend / step 2.9 Little difference in heterosc.
Trends serially independent data, variations due to independent processes
Steps perform better than trends to explain autocorrelation and heteroscedasticity
None located
Autocorrelation and heteroscedasticity quarterly observations satellite temp 19792014
Accumulated error trend / step 4.4, 0.9 and 3.1, 2.1 RSS and UAHTrends and steps show little difference in autoc. and heterosc. (except for steps with 24 % v 8 % heterosc.)
Trends serially independent data, variations due to independent processes
Steps perform better than trends to explain cumulative error, little difference between autocorrelation and heteroscedasticity
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
202 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
could be made to each of these on an individual basis, collectively they show that for externally forced warming on decadal scales, hstep is better supported than htrend.
In summary, these tests show that hstep is a close approximation of the data when analysing decadal-scale warming.Over the long term, this warming conforms to a complex trend that can be simplied as a monotonic curve, but the actual pathway is step-like. As outlined in Sect. 3.3, this rules out gradual warming, either in situ in the atmosphere or as a gradual release from the ocean, in favour of a more abrupt process of storage and release. This conclusion supports the substantive hypothesis H2 over H1, where the climate change and variability interact, rather than varying independently.
7 Proposed mechanisms for step-like warming
The correlation between step-like warming and ECS in the models, between the timing of steps in model hindcasts and observations and between steps and known regime changes in observations (Table 6), provides strong evidence that warming is non-gradual on decadal timescales. The high correlations of steps and shifts with model ECS indicate that atmospheric feedback processes respond to abrupt releases of heat into the atmosphere. The presence of negligible internal trends occurring over some oceanic regions, the region 30 60 N, and in tropospheric satellite temperatures, suggests that little of the heat being trapped in the atmosphere by anthropogenic greenhouse gases actually remains there.
One justication given for rejecting externally driven step-like warming is that it is presumed that there is no plausible physical mechanism for this (Cahill et al., 2015; Foster and Abraham, 2015). However, to suggest that the stepwise release of heat energy is physically implausible overlooks the energetics of the oceanatmosphere system. Hydrodynamic processes are quite capable of supplying the energy required (Ozawa et al., 2003; Lucarini and Ragone, 2011; Ghil, 2012).The atmosphere contains as much heat energy as the top3.2 m of ocean (Bureau of Meteorology, 2003). About 93 % of historically added heat currently resides in the ocean (Levitus et al., 2012; Roemmich et al., 2015), whereas the atmosphere contains about 3 % of the total. A similar amount of the heat has been stored within the land mass (Balmaseda et al., 2013) and on an annual basis a similar ux is absorbed in melting ice (Hansen et al., 2011). A physical reorganisation of the oceanatmosphere system, as part of a regime change, is therefore large enough to provide the relatively small amount of energy required to cause abrupt sea surface and atmospheric warming (Roemmich et al., 2015;Reid et al., 2016), as shown by rapid changes in shallow ocean heat content (Fig. 6b; Roemmich and Gilson, 2011;Reid, 2016).
For example, Reid et al. (2016) in describing the late 1980s regime change, show it was associated with large-
scale shifts in temperature and multiple impacts across terrestrial and marine systems, mainly in the Northern Hemisphere. Changes in the North Pacic in 1977 were considered even more extensive (Hare and Mantua, 2000), as were those in 1997/98 involving both the Pacic and Atlantic oceans (Chikamoto et al., 2012a, b). In developing tests for detection and attribution, Jones (2012) noted two types of regime change over land: one where codependent variables such as maximum temperature and rainfall undergo a step change but remain in a stationary relationship, and the other, non-stationary change, where warming undergoes a step change independent of rainfall change. This suggests that although regime changes are a normal part of internal climate variability, they can be enhanced, releasing extra heat. The step changes summarised in Table 6 coincide with El Nio events but the heat emitted by other El Nio events dissipates and is absorbed back into the ocean within months; thus, an added mechanism is required. We propose that there is negligible in situ atmospheric warming and that almost all of the added heat trapped by anthropogenic greenhouse gases is absorbed by and stored in the ocean. It is subsequently released through the action of oscillatory mechanisms associated with regime shifts.
Most heat (long-wave radiation) is trapped near the ground or ocean surface and much of that is radiated downwards (Trenberth, 2011). The atmosphere as a whole has little intrinsic heat memory and does not warm independently of the surface. This is supported by observations on land where the overpassing air mass takes on the characteristics of the underlying surface, achieving energy balance within a 300 m distance (Morton, 1983). When passing from land to water, this will see all of the available heat energy taken up by water if the temperature of the air mass exceeds that of water (Morton, 1983, 1986), with the temperature of the overpassing air mass reaching equilibrium with the water beneath within a very short time. Very little of the heat trapped over land can be absorbed by the land surface, but will be transported from land to ocean within a few days to a few weeks, where it can be absorbed (the high latitudes being an exception). Given that the atmosphere interacts with the top 70 m of ocean over an annual cycle (Hartmann, 1994), there is ample opportunity for the majority of available heat trapped over land that is not absorbed by land, lakes and ice to be absorbed by the ocean.
In terms of energy budgets, the additional direct forcing from anthropogenic greenhouse gases is roughly 1.5 %(2.3 Wm2, IPCC 2013) of the estimated total annual budget of 155 Wm2 trapped mainly by water vapour and CO2 (Schmidt et al., 2010). Since > 90 % of that 1.5 % is already accepted as being absorbed by the ocean, it is not clear why the roughly 3 % of that 2 % (0.07 Wm2) not absorbed by land, snow and ice would remain in the atmosphere if its absorption by the ocean is not energy limited,i.e. in the low to mid-latitudes. Negligible internal trends in lower tropospheric satellite temperatures also indicate that
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 203
the air column is not warming in situ but exhibits stable temperatures punctuated by step changes (Fig. 4). This suggests that climate forms a series of oscillating steady-state regimes, with the temperature of the atmosphere being controlled by oceanatmosphere interactions.
Step-like warming requires a trigger and release mechanism. Recently, Peyser et al. (2016) linked dynamic sea level in the Pacic Ocean, measured using an eastwest see-saw index, to rapid changes in global mean surface temperature.In 19961997, that index underwent a west-to-east see-saw movement of 149 mm. This would mark the release of a large tongue of warm water from the western Pacic warm pool to the east, making heat available for discharge into the atmosphere. Based on a linear regression between the see-saw index and surface temperature calculated from control runs of 38 CMIP5 climate models, they estimate a jump in surface temperature of 0.29 0.10 C in 19971998, close to
our estimate of 0.32, or 0.25 C if 1987/88 is taken into account. They estimated another see-saw change of 111 mm in 2014/15 as contributing to a rapid warming of 0.21 0.07 C
in 2016. We interpret their observations of rapid sea level rise in the western Pacic region as representing the sustained storage of heat in the Indo-Pacic warm pool. Heat absorbed in the tropical Pacic is blown westward into the warm pool, where it accumulates, maintaining the tropical Pacic as a region of generally low warming (Power et al., 2016). As the warm pool reaches critical limits, it becomes unstable, releasing surplus heat as a tongue of warm water from the western to eastern Pacic during an El Nio event.
Meehl et al. (2016) have also suggested that the negative phase of the Interdecadal Pacic Oscillation that commenced in 1997/98 (Overland et al., 2008; Meehl et al., 2013) could change to positive during 20152019 as part of oscillatory mechanisms associated with the build-up of heat in the western Pacic. OKane et al. (2014) provide evidence that such changes may be identied years in advance. An accompanying regime change emplacing large areas of warmer water required to sustain higher temperatures after the initial outburst is consistent with widespread coral bleaching in 2014 2016 (Normile, 2016), rivalling that of 1998. Note that both Peyser et al. (2016) and Meehl et al. (2013) interpret their results as variability acting on a long-term trend; however, we reinterpret their ndings as supporting a heat pulse and regime change, producing step-like warming.
In storing heat for redistribution, the Indo-Pacic warm pool acts as a global heat engine (Bosc et al., 2009), a function it has fullled for millions of years over a wide range of climatic changes (Gagan et al., 2004; de Garidel-Thoron et al., 2005; Abram et al., 2009). The storage and release mechanism identied by Peyser et al. (2016) may therefore be an additional response to a build-up of heat over and above oscillations associated with ongoing decadal regime change.Storage and release mechanisms may exist in other ocean basins but would need to be identied.
8 Discussion
There are many reasons as to why H1 where climate change and variability are considered to be independent of each other has dominated climate research despite the lack of a conclusive theoretical or statistical case. They include historical, social, theoretical and political considerations too broad to cover here.
Benestad (2016) reviews models used to build a mental picture of the greenhouse effect, nominating radiative convective and heat balance models as two types historically used for this purpose. He describes the basic processes of radiative transfer as being well understood but insufcient to explain the warming process. Radiative transfer theory constitutes core greenhouse theory. However, the subsequent process of heat diffusion through the climate system is less well understood, although the understanding that if greenhouse gases are increased, the atmosphere will warm until the radiative balance at the top of the atmosphere is achieved also constitutes core theory.
Our conclusion that the atmosphere does not warm in situ will challenge many who consider that to be a basic part of the greenhouse effect. However, an exhaustive search of the literature failed to nd any direct evidence that this actually takes place. We nd it hard to perceive how an additional increment of long-wave radiation on the order of 0.2 Wm2
(direct forcing and feedback derived from Schmidt et al. 2010) can behave differently to the 155 Wm2 produced
in the atmosphere year to year without being absorbed by the wider climate system. Given that climate models exhibit step-like warming, where the abrupt component carries the greater part of the signal than internal trends, they produce emergent behaviour that is not identied by mainstream analytic approaches.
Overwhelmingly, model- and statistically based studies represent the global warming signal as changing gradually.Some are prescriptive because of their structure or because they apply simplied assumptions about a more complex climate system, other models examine a small part of the system, and some have a historical legacy bestowing familiarity and reliability. Modern climate models are almost as complex as the climate, and thus need to be understood through simpler models (Held, 2005; Benestad, 2016), forming a nested modelling approach from simple through to complex (Schneider and Dickinson, 1974; Ghil, 2015). The linking of trend analysis methods with gradual change may overlook the distinction between process-based and diagnostic models. A diagnostic model may identify a trend without necessarily indicating a gradual process. A large part of the climate wars has been fought over this very point.
Nonlinear responses in climate are being investigated by researchers, with an interest in complex system behaviour via dynamical systems and related theory. Our conclusions suggest that the processes of radiative transfer and subsequent warming take place in two separate domains of the climate
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
elevator-like record of change. Whether the oscillatory systems themselves change under greater forcing (e.g. RCP8.5) or whether warming itself becomes more diffuse has yet to be investigated. Note that these step changes are quite different to those catalogued by Drijfhout et al. (2015), who used a different method to screen the CMIP5 model ensemble for abrupt shifts that could be considered as singularities, locating 37 ocean, sea ice, snow cover, permafrost and terrestrial biosphere changes.
Statistical characterisations of changing climate variables are becoming more probabilistic, with probability distribution functions increasingly being produced from climate model ensembles. However, the presence of non-gradual change suggests that statistics developed from the path-wise analysis of individual simulations (as was carried out in this paper and as suggested by Ghil, 2015) are required, especially higher-order statistics that represent extreme events potentially subject to step changes. For example, re risk in Victoria, Australia, increased abruptly by 38 % between 197297 and 19982010, driven by a step change in climate (Jones et al., 2013). Because methods for detection, attribution, climate forecasting and characterisation of future climate risk are almost totally dependent on being scaled to gradual change in mean variables, a stepwise process will require a substantial rethink as to how these activities can be conceptualised.
For example, seamless links between weather and climate forecasting over a range of timescales are proposed as a key scientic target (Palmer et al., 2008; Hoskins, 2013). The Global Framework for Climate Services (World Meteorological Organization, 2011) reects that
Weather and climate research are closely intertwined; progress in our understanding of climate processes and their numerical representation is common to both. Seamless prediction (on timescales from a few hours to centuries) needs to be further developed and extended to aspects across multiple disciplines relevant to climate processes (World Meteorological Organization, 2010).
Solomon et al. (2011) state that Long experience in weather and climate forecasting has shown that forecasts are of little utility without a priori assessment of forecast skill and reliability. The assumption that the processes involved are timescale invariant indicates that what seamless prediction means in a decision-support context has not been fully thought through. For the moment, decadal prediction concentrates on ensemble mean change in variables that show skill in climate models, whereas the prospect of non-gradual change carries the greater risk. Under this type of framing, climate services remain supply driven rather than demand driven (Gunasekera et al., 2014; Street, 2016). Projections of mean change also overlook the considerable literature on scenarios that have arisen because of the failure of multi-year
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
204 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
system, separated by a delay. The absorption of radiation is a linear process that is quite separate from the behaviour of turbulent dissipation of heat energy within the climate system, which is fundamentally nonlinear (Ozawa et al., 2003).Developments based on deterministic nonlinear and stochastic linear behaviour originating from work by Lorenz (1963) and Hasselmann (1976) respectively explore a range of interrelated phenomena, such as non-equilibrium stable states, oscillators, strange attractors, bifurcations and entropy production, in order to develop a unied theory of climate (Ozawa et al., 2003; Lucarini et al., 2014; Franzke et al., 2015; Ghil, 2015). Studying how the free and forced aspects of change combine to alter the statistical properties of climate is a specic goal (Lucarini and Sarno, 2011; Ghil, 2012, 2015).
Our focus is on understanding the role of linear and nonlinear behaviour in changing climate risk over decadal timescales, specically how initial condition and boundary-limited uncertainties (as described by Lorenz, 1975 and Hasselmann, 2002) combine. Initial-condition uncertainty is boundary limited, varying within a certain amplitude, with the outcome depending on the pathway taken within those limits (Lorenz, 1975). There is also a time-dependent window that serves as a predictability barrier. Changing boundary conditions are intransitive, with the outcome being insensitive to initial conditions. The nested nature of climate phenomena over different timescales results in decadal-scale climates being both an initial-condition and intransitive process, combining to produce stochastically driven step changes in warming that integrate into a long-term complex trend. The coincident timing of step changes in both observations and models (Fig. 7) suggests that other factors, such as short-term volcanic forcing, can also inuence the timing of step changes.
Lorenz (1968) referred to the outcome of forced climate change on century timescales as almost intransitive. The almost is due to initial condition uncertainties operating within the boundary limitations of decadal variability. The almost-intransitive model (Lorenz, 1968) is described via linear response theory (Lucarini et al., 2014; Ragone et al., 2016) and shown to be robust for concepts such as effective radiative forcing (Hansen et al., 2005) and effective climate sensitivity (Andrews et al., 2015), although these phenomena would be sensitive to bifurcations if they were to occur (Hasselmann, 2002).
If the ocean takes up the additional available heat from anthropogenic greenhouse gases while maintaining steady-state conditions within an oscillatory system of climate regimes, it can be considered as acting homeostatically with respect to the atmosphere (e.g. Kleidon, 2004). Heat will accumulate in the shallow ocean until such a time as it becomes unstable and is released as part of a stepwise regime change. The new regime, being warmer, enhances vertical and horizontal heat uxes, which is consistent with a more energetic system.Sustained forcing would produce a series of regime changes becoming successively warmer and forming a stepladder or
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 205
predictions of mean change in systems that exhibit considerable non-linearity (Wack, 1985a, b; Brjeson et al., 2006).
9 Conclusions
Here, we have adapted and applied severe testing principles proposed by Mayo and Spanos (2010) to determine the role that step changes play in decadal-scale warming.This involves the linking of scientic hypotheses H1 and H2 with statistical hypotheses htrend and hstep and subjecting them to severe testing. Paraphrasing the severity principle of Mayo (2010), the results of Tests 16 provide evidence for hypothesis H2 if and only if hstep passes a severe test with very high probability, where htrend would have uncovered the falsity of H2, and yet no such error is detected.Error and probative testing of steps against trends lends little support for the proposition that the climate warms gradually.If trend-like behaviour were dominating warming or were on an even footing with step-like change, these tests would have identied it. H1 is only suitable for intransitive estimates of change, where the initial conditions, pathway and non-linear components of forcing are unimportant.
Surface and tropospheric warming on decadal timescales is dominated by stepwise changes in temperature (Reid and Beaugrand, 2012; Jones et al., 2013; Belolipetsky et al., 2015; Bartsev et al., 2016; Reid et al., 2016). The basic physical mechanism for moving from H1 to H2 is deceptively simple: instead of warming occurring in situ in the atmosphere and/or being released gradually from the ocean, all available heat from additional greenhouse gases not absorbed by the land surface, snow and ice and in lakes is absorbed by the ocean. There, it is entrained into the nonlinear processes of climate variability, where the added forcing interacts with those processes. The most plausible explanation for step-like behaviour is that steady-state decadal regimes are punctuated by step-like bursts of warming that are subsequently maintained by higher sea surface temperature emplaced by ocean atmosphere regime changes.
This conclusion does not invalidate the considerable literature that assesses long-term (> 50 years) climate change as a relatively linear process and the warming response as being broadly additive with respect to forcing (e.g. Lucarini et al., 2010; Marvel et al., 2015). However, the signal-to-noise model of a gradually changing mean surrounded by random climate variability poorly represents warming on decadal timescales. The separation of signal and noise into good and bad is likewise poor framing for the purposes of understanding and managing risk in fundamentally nonlinear systems (Koutsoyiannis, 2010). As we show, the presence of such changes within climate models does not indicate a need to fundamentally change how climate modelling is carried out. It does, however, indicate a need to change how the results are analysed.
Climate conceptualised as a mechanistic system and described using classical statistical methods is substantially different from climate conceptualised as a complex system. With record atmospheric and surface ocean temperatures in 2015/16 variously being described as a singular event, a rein-vigoration of trend-like warming or a wholesale shift to a new climate regime, this issue is too important to be left unresolved.
Code availability. Included with the Supplement as a zip le (Python and R modules).
Data availability. Included with the Supplement as Excel les.
The Supplement related to this article is available online at http://dx.doi.org/10.5194/esd-8-177-2017-supplement
Web End =doi:10.5194/esd-8-177-2017-supplement .
Author contributions. R. N. Jones conceived the study,J. H. Ricketts coded and tested the multistep model, R. N. Jones developed the severe testing regime for the results and withJ. H. Ricketts undertook analyses, J. H. Ricketts put together the SI and maintained quality control and R. N. Jones led the paper with contributions from J. H. Ricketts.
Competing interests. The authors declare that they have no conict of interest.
Acknowledgements. J. H. Ricketts is the holder of a Victoria University postgraduate research scholarship. Data sources include the Met Ofce Hadley Centre, National Aeronautics and Space Administration Goddard Institute for Space Studies, National Oceanographic Data Center and United States National Climatic Data Center, Berkeley Earth, Cowtan and Way, Permanent Service for Mean Sea Level, and the Australian Bureau of Meteorology. CMIP3 and CMIP5 archives are made available by the modelling groups, the Program for Climate Model Diagnosis and Intercom-parison (PCMDI) and the WCRPs Working Group on Coupled Modelling (WGCM). The U.S. Department of Energys Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. D. Kelly ODay provided the macro templates, which have been adapted to provide the step and trend charts. Two anonymous reviewers and the editor made valuable comments that helped focus the paper.
Edited by: M. CrucixReviewed by: two anonymous referees
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
206 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
References
Abram, N. J., McGregor, H. V., Gagan, M. K., Hantoro, W. S., and Suwargadi, B. W.: Oscillations in the southern extent of the Indo-Pacic Warm Pool during the mid-Holocene, Quaternary Sci.Rev., 28, 27942803, 2009.
Alfaro, E. J., Gershunov, A., and Cayan, D.: Prediction of Summer Maximum and Minimum Temperature over the Central and Western United States: The Roles of Soil Moisture and Sea Surface Temperature, J. Climate, 19, 14071421, doi:http://dx.doi.org/10.1175/JCLI3665.1
Web End =10.1175/JCLI3665.1 http://dx.doi.org/10.1175/JCLI3665.1
Web End = , 2006.
Andrews, T., Gregory, J. M., and Webb, M. J.: The dependence of radiative forcing and feedback on evolving patterns of surface temperature change in climate models, J. Climate, 28, 1630 1648, 2015.
Balmaseda, M. A., Trenberth, K. E., and Klln, E.: Distinctive climate signals in reanalysis of global ocean heat content, Geophys.Res. Lett., 40, 17541759, doi:http://dx.doi.org/10.1002/grl.50382
Web End =10.1002/grl.50382 http://dx.doi.org/10.1002/grl.50382
Web End = , 2013.Bartsev, S. I., Belolipetskii, P. V., Degermendzhi, A. G., Ivanova,Y. D., Pochekutov, A. A., and Saltykov, M. Y.: Refocusing on the dynamics of the Earths climate, Her. Russ. Acad. Sci., 86, 135142, doi:http://dx.doi.org/10.1134/s1019331616020015
Web End =10.1134/s1019331616020015 http://dx.doi.org/10.1134/s1019331616020015
Web End = , 2016.Belolipetsky, P. V.: The Shifts Hypothesis an alternative view of global climate change, arXiv preprint arXiv:1406.5805, arXiv:1406.5805, 2014.
Belolipetsky, P. V., Bartsev, S., Ivanova, Y., and Saltykov, M.: Hidden staircase signal in recent climate dynamic, Asia-Pac. J. Atmos. Sci., 51, 323330, doi:http://dx.doi.org/10.1007/s13143-015-0081-6
Web End =10.1007/s13143-015-0081-6 http://dx.doi.org/10.1007/s13143-015-0081-6
Web End = , 2015.Benestad, R. E.: A mental picture of the greenhouse effect, Theor.
Appl. Climatol., 110, doi:http://dx.doi.org/10.1007/s00704-016-1732-y
Web End =10.1007/s00704-016-1732-y http://dx.doi.org/10.1007/s00704-016-1732-y
Web End = , 2016.
Brjeson, L., Hjer, M., Dreborg, K.-H., Ekvall, T., and Finnveden,G.: Scenario types and techniques: towards a users guide, Futures, 38, 723739, 2006.
Bosc, C., Delcroix, T., and Maes, C.: Barrier layer variability in the western Pacic warm pool from 2000 to 2007, J. Geophys. Res.-Oceans, 114, C06023, doi:http://dx.doi.org/10.1029/2008JC005187
Web End =10.1029/2008JC005187 http://dx.doi.org/10.1029/2008JC005187
Web End = , 2009.Boucharel, J., Dewitte, B., Garel, B., and du Penhoat, Y.:
ENSOs non-stationary and non-Gaussian character: the role of climate shifts, Nonlin. Processes Geophys., 16, 453473, doi:http://dx.doi.org/10.5194/npg-16-453-2009
Web End =10.5194/npg-16-453-2009 http://dx.doi.org/10.5194/npg-16-453-2009
Web End = , 2009.
Boucharel, J., Dewitte, B., Penhoat, Y., Garel, B., Yeh, S.-W., and
Kug, J.-S.: ENSO nonlinearity in a warming climate, Clim. Dynam., 37, 20452065, doi:http://dx.doi.org/10.1007/s00382-011-1119-9
Web End =10.1007/s00382-011-1119-9 http://dx.doi.org/10.1007/s00382-011-1119-9
Web End = , 2011.Branstator, G. and Selten, F.: Modes of Variability and Climate Change, J. Climate, 22, 26392658, doi:http://dx.doi.org/10.1175/2008JCLI2517.1
Web End =10.1175/2008JCLI2517.1 http://dx.doi.org/10.1175/2008JCLI2517.1
Web End = , 2009.
Bcher, A. and Dessens, J.: Secular trend of surface temperature at an elevated observatory in the Pyrenees, J. Climate, 4, 859868, 1991.
Buishand, T.: Tests for detecting a shift in the mean of hydrological time series, J. Hydrol., 73, 5169, 1984.
Bureau of Meteorology: The greenhouse effect and climate change,Bureau of Meteorology, Melbourne, 2003.
Cahill, N., Rahmstorf, S., and Parnell, A. C.: Change points of global temperature, Environ. Res. Lett., 10, 084002, doi:http://dx.doi.org/1748-9326/10/8/084002
Web End =1748- http://dx.doi.org/1748-9326/10/8/084002
Web End =9326/10/8/084002 , 2015.
Cai, W. and Cowan, T.: SAM and regional rainfall in IPCC AR4 models: Can anthropogenic forcing account for southwest West-
ern Australian winter rainfall reduction?, Geophys. Res. Lett., 33, L24708, doi:http://dx.doi.org/10.1029/2006GL028037
Web End =10.1029/2006GL028037 http://dx.doi.org/10.1029/2006GL028037
Web End = , 2006.
Capparelli, V., Franzke, C., Vecchio, A., Freeman, M. P., Watkins,N. W., and Carbone, V.: A spatiotemporal analysis of US station temperature trends over the last century, J. Geophys. Res.-Atmos., 118, 74277434, 2013.
Chikamoto, Y., Kimoto, M., Ishii, M., Watanabe, M., Nozawa, T., Mochizuki, T., Tatebe, H., Sakamoto, T. T., Komuro, Y., and Shiogama, H.: Predictability of a stepwise shift in Pacic climate during the late 1990s in hindcast experiments using MIROC, J. Meteorol. Soc. Jpn, 90, 121, 2012a.
Chikamoto, Y., Kimoto, M., Watanabe, M., Ishii, M., and Mochizuki, T.: Relationship between the Pacic and Atlantic stepwise climate change during the 1990s, Geophys. Res. Lett., 39, L21710, 2012b.
Christy, J. R., Spencer, R. W., Norris, W. B., Braswell, W. D., and Parker, D. E.: Error estimates of version 5.0 of MSU-AMSU bulk atmospheric temperatures, J. Atmos. Ocean. Tech., 20, 613629, 2003.
Christy, J. R., Norris, W. B., Spencer, R. W., and Hnilo, J. J.: Tropospheric temperature change since 1979 from tropical radiosonde and satellite measurements, J. Geophys. Res.-Atmos., 112, D06102, doi:http://dx.doi.org/10.1029/2005JD006881
Web End =10.1029/2005JD006881 http://dx.doi.org/10.1029/2005JD006881
Web End = , 2007.
Church, J. A., Hunter, J. R., McInnes, K. L., and White, N. J.: Sea-level rise around the Australian coastline and the changing frequency of extreme sea-level events, Aust. Meteorol. Mag., 55, 253260, 2006.
Cleveland, W. S. and Devlin, S. J.: Locally weighted regression: an approach to regression analysis by local tting, J. Am. Stat. Assoc., 83, 596610, 1988.
Corti, S., Molteni, F., and Palmer, T. N.: Signature of recent climate change in frequencies of natural atmospheric circulation regimes, Nature, 398, 799802, 1999.
Cowtan, K. and Way, R. G.: Coverage bias in the HadCRUT4 temperature series and its impact on recent temperature trends, Q. J. Roy. Meteor. Soc., 140, 19351944, doi:http://dx.doi.org/10.1002/qj.2297
Web End =10.1002/qj.2297 http://dx.doi.org/10.1002/qj.2297
Web End = , 2014. Dai, A., Fyfe, J. C., Xie, S.-P., and Dai, X.: Decadal modulation of global surface temperature by internal climate variability, Nature Climate Change, 5, 555559, 2015.de Garidel-Thoron, T., Rosenthal, Y., Bassinot, F., and Beaufort, L.:
Stable sea surface temperatures in the western Pacic warm pool over the past 1.75 million years, Nature, 433, 294298, 2005. Delworth, T. L. and Zeng, F.: Regional rainfall decline in Australia attributed to anthropogenic greenhouse gases and ozone levels, Nature Geosci., 7, 583587, doi:http://dx.doi.org/10.1038/ngeo2201
Web End =10.1038/ngeo2201 http://dx.doi.org/10.1038/ngeo2201
Web End = , 2014. Drijfhout, S., Bathiany, S., Beaulieu, C., Brovkin, V., Claussen, M.,
Huntingford, C., Scheffer, M., Sgubin, G., and Swingedouw, D.: Catalogue of abrupt shifts in Intergovernmental Panel on Climate Change climate models, P. Natl. Acad. Sci., 112, E5777E5786, 2015.
Ebbesmeyer, C. C., Cayan, D. R., McLain, D. R., Nichols, F. H.,
Peterson, D. H., and Redmond, K. T.: 1976 step in the Pacic climate: forty environmental changes between 19681975 and 19771984, Proceedings of the Seventh annual Pacic Climate (PACLIM) Workshop, Asilomar, California, April, 1990, 1991. England, M. H., McGregor, S., Spence, P., Meehl, G. A., Timmer-mann, A., Cai, W., Gupta, A. S., McPhaden, M. J., Purich, A., and Santoso, A.: Recent intensication of wind-driven circula-
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 207
tion in the Pacic and the ongoing warming hiatus, Nature Climate Change, 4, 222227, 2014.
Fischer, T., Gemmer, M., Liu, L., and Su, B.: Change-points in climate extremes in the Zhujiang River Basin, South China, 1961 2007, Climatic Change, 110, 783799, doi:http://dx.doi.org/10.1007/s10584-011-0123-8
Web End =10.1007/s10584-011- http://dx.doi.org/10.1007/s10584-011-0123-8
Web End =0123-8 , 2012.
Foster, G. and Abraham, J.: Lack of evidence for a slowdown in global temperature, US CLIVAR, 13, 69, 2015.
Foster, G. and Rahmstorf, S.: Global temperature evolution 19792010, Environ. Res. Lett., 6, 044022, 2011.
Franzke, C.: Nonlinear trends, long-range dependence, and climate noise properties of surface temperature, J. Climate, 25, 4172 4183, 2012.
Franzke, C. L., OKane, T. J., Berner, J., Williams, P. D., and Lucarini, V.: Stochastic climate theory and modeling, Climate Change, 6, 6378, 2015.
Free, M. and Lanzante, J.: Effect of volcanic eruptions on the vertical temperature prole in radiosonde data and climate models, J.Climate, 22, 29252939, 2009.
Gagan, M. K., Hendy, E. J., Haberle, S. G., and Hantoro, W. S.:
Post-glacial evolution of the Indo-Pacic warm pool and El Nino-Southern Oscillation, Quatern. Int., 118, 127143, 2004.Ghil, M.: Climate variability: nonlinear and random effects, Encyclopedia of Atmospheric Sciences, Elsevier, 16, 2012.
Ghil, M.: A mathematical theory of climate sensitivity or, How to deal with both anthropogenic forcing and natural variability?, in: Climate Change: Multidecadal and Beyond, edited by: Chang,C.-P., Ghil, M., Latif, M., and Wallace, J. M., World Scientic Publishing Company, London, Singapore, 2015.
Gunasekera, D., Troccoli, A., and Boulahya, M. S.: Energy and meteorology: partnership for the future, in: Weather Matters for Energy, edited by: Troccoli, A., Dubus, L., and Haupt, S. E., Springer, 497511, 2014.
Haig, B. D.: Tests of statistical signicance made sound, Educ. Psychol. Meas., doi:http://dx.doi.org/10.1177/0013164416667981
Web End =10.1177/0013164416667981 http://dx.doi.org/10.1177/0013164416667981
Web End = , 2016.
Hansen, J., Sato, M., Ruedy, R., Nazarenko, L., Lacis, A., Schmidt,G. A., Russell, G., Aleinov, I., Bauer, M., Bauer, S., Bell, N., Cairns, B., Canuto, V., Chandler, M., Cheng, Y., Del Genio, A., Faluvegi, G., Fleming, E., Friend, A., Hall, T., Jackman, C., Kelley, M., Kiang, N., Koch, D., Lean, J., Lerner, J., Lo, K., Menon,S., Miller, R., Minnis, P., Novakov, T., Oinas, V., Perlwitz, J., Perlwitz, J., Rind, D., Romanou, A., Shindell, D., Stone, P., Sun,S., Tausnev, N., Thresher, D., Wielicki, B., Wong, T., Yao, M., and Zhang, S.: Efcacy of climate forcings, J. Geophys. Res.-Atmos., 110, D18104, doi:http://dx.doi.org/10.1029/2005JD005776
Web End =10.1029/2005JD005776 http://dx.doi.org/10.1029/2005JD005776
Web End = , 2005.Hansen, J., Ruedy, R., Sato, M., and Lo, K.: Global surface temperature change, Rev. Geophys., 48, RG4004, 2010.
Hansen, J., Sato, M., Kharecha, P., and von Schuckmann, K.:
Earths energy imbalance and implications, Atmos. Chem. Phys., 11, 1342113449, doi:http://dx.doi.org/10.5194/acp-11-13421-2011
Web End =10.5194/acp-11-13421-2011 http://dx.doi.org/10.5194/acp-11-13421-2011
Web End = , 2011.Hare, S. R. and Mantua, N. J.: Empirical evidence for North Pacic regime shifts in 1977 and 1989, Prog. Oceanogr., 47, 103145, 2000.
Hartmann, D. L.: Global Physical Climatology, International Geophysics, Academic Press, San Diego USA and London UK, 411 pp., 1994.
Hasselmann, K.: Stochastic climate models part I. Theory, Tellus,28, 473485, 1976.
Hasselmann, K.: Is Climate Predictable?, in: The Science of Disasters: Climate Disruptions, Heart Attacks, and Market Crashes, edited by: Bunde, A., Kropp, J., and Schellnhuber, H. J., Springer, Berlin Heidelberg, 141188, 2002.
Hegerl, G. and Zwiers, F.: Use of models in detection and attribution of climate change, Climate Change, 2, 570591, 2011.
Held, I. M.: The Gap between Simulation and Understanding in Climate Modeling, B. Am. Meteorol. Soc., 86, 16091614, doi:http://dx.doi.org/10.1175/BAMS-86-11-1609
Web End =10.1175/BAMS-86-11-1609 http://dx.doi.org/10.1175/BAMS-86-11-1609
Web End = , 2005.
Hope, P. K., Drosdowsky, W., and Nicholls, N.: Shifts in the synoptic systems inuencing southwest Western Australia, Clim. Dynam., 26, 751764, doi:http://dx.doi.org/10.1007/s00382-006-0115-y
Web End =10.1007/s00382-006-0115-y http://dx.doi.org/10.1007/s00382-006-0115-y
Web End = , 2006.
Hope, P., Timbal, B., and Fawcett, R.: Associations between rainfall variability in the southwest and southeast of Australia and their evolution through time, Int. J. Climatol., 30, 13601372, doi:http://dx.doi.org/10.1002/joc.1964
Web End =10.1002/joc.1964 http://dx.doi.org/10.1002/joc.1964
Web End = , 2010.
Hoskins, B.: The potential for skill across the range of the seamless weather-climate prediction problem: a stimulus for our science,Q. J. Roy. Meteor. Soc., 139, 573584, 2013.
IPCC: Climate Change 2013: The Physical Science Basis, in: Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P.M., Cambridge University Press, Cambridge, UK and New York, USA, 1535 pp., 2013.
Ji, F., Wu, Z., Huang, J., and Chassignet, E. P.: Evolution of land surface air temperature trend, Nature Climate Change, 4, 462 466, doi:http://dx.doi.org/10.1038/nclimate2223
Web End =10.1038/nclimate2223 http://dx.doi.org/10.1038/nclimate2223
Web End = , 2014.
Jones, R. N.: Modelling the Hydrologic and Climatic Controls of
Closed Lakes, Western Victoria, PhD, Civil and Environmental Engineering, University of Melbourne, 283 pp., 1995.
Jones, R. N.: North central Victorian climate: past, present and future, Proceedings of the Royal Society of Victoria, 122, 147160, 2010.
Jones, R. N.: Detecting and attributing nonlinear anthropogenic regional warming in southeastern Australia, J. Geophys. Res., 117, D04105, doi:http://dx.doi.org/10.1029/2011jd016328
Web End =10.1029/2011jd016328 http://dx.doi.org/10.1029/2011jd016328
Web End = , 2012.
Jones, R. N., Young, C. K., Handmer, J., Keating, A., Mekala, G. D., and Sheehan, P.: Valuing Adaptation under Rapid Change, National Climate Change Adaptation Research Facility, Gold Coast, Australia, 182 pp., 2013.
Karoly, D. J. and Braganza, K.: A new approach to detection of anthropogenic temperature changes in the Australian region, Mete-orol. Atmos. Phys., 89, 5767, doi:http://dx.doi.org/10.1007/s00703-005-0121-3
Web End =10.1007/s00703-005-0121-3 http://dx.doi.org/10.1007/s00703-005-0121-3
Web End = , 2005.
Kirono, D. and Jones, R.: A bivariate test for detecting inhomogeneities in pan evaporation time series, Aust. Meteorol. Mag., 56, 93103, 2007.
Kirtman, B., Power, S., Adedoyin, A. J., Boer, G., Bojariu, R.,
Camilloni, I., Doblas-Reyes, F., Fiore, A., Kimoto, M., Meehl,G., Prather, M., Sarr, A., Schr, C., Sutton, R., Oldenborgh, G. J.v., Vecchi, G., and Wang, H.-J.: Near-term Climate Change: Projections and Predictability, in: Climate Change 2013: The Physical Science Basis. Working Group I contribution to the IPCC 5th Assessment Report, edited by: Stocker, T. F., Qin, D., Plattner,G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge and New York, 121, 2013.
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
208 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
Kleidon, A.: Beyond Gaia: Thermodynamics of Life and
Earth System Functioning, Climatic Change, 66, 271319, doi:http://dx.doi.org/10.1023/B:CLIM.0000044616.34867.ec
Web End =10.1023/B:CLIM.0000044616.34867.ec http://dx.doi.org/10.1023/B:CLIM.0000044616.34867.ec
Web End = , 2004.Koutsoyiannis, D.: HESS Opinions A random walk on water,
Hydrol. Earth Syst. Sci., 14, 585601, doi:http://dx.doi.org/10.5194/hess-14-585-2010
Web End =10.5194/hess-14-585- http://dx.doi.org/10.5194/hess-14-585-2010
Web End =2010 , 2010.
LHte, Y., Mah, G. I. L., Som, B., and Triboulet, J. P.: Analysis of a Sahelian annual rainfall index from 1896 to 2000, the drought continues, Hydrolog. Sci. J., 47, 563572, doi:http://dx.doi.org/10.1080/02626660209492960
Web End =10.1080/02626660209492960 http://dx.doi.org/10.1080/02626660209492960
Web End = , 2002.
Legates, D. R., Soon, W., and Briggs, W. M.: Climate consensus and misinformation: A rejoinder to Agnotology, scientic consensus, and the teaching and learning of climate change, Sci. Educ.-Netherlands, 24, 299318, 2015.
Levitus, S., Antonov, J. I., Boyer, T. P., Baranova, O. K., Garcia, H.E., Locarnini, R. A., Mishonov, A. V., Reagan, J. R., Seidov, D., Yarosh, E. S., and Zweng, M. M.: World ocean heat content and thermosteric sea level change (02000 m), 19552010, Geophys.Res. Lett., 39, L10603, doi:http://dx.doi.org/10.1029/2012GL051106
Web End =10.1029/2012GL051106 http://dx.doi.org/10.1029/2012GL051106
Web End = , 2012.
Li, F., Chambers, L., and Nicholls, N.: Relationships between rainfall in the southwest of Western Australia and near global patterns of sea-surface temperature and mean sea-level pressure variability, Aust. Meteorol. Mag., 54, 2333, 2005.
Lo, T. T. and Hsu, H. H.: Change in the dominant decadal patterns and the late 1980s abrupt warming in the extratropical Northern Hemisphere, Atmos. Sci. Lett., 11, 210215, 2010.
Lorenz, E. N.: Deterministic nonperiodic ow, J. Atmos. Sci., 20,130141, 1963.
Lorenz, E. N.: Climatic determinism, Meteor. Mon., 8, 13, 1968. Lorenz, E. N.: Climate Predictability, in: The Physical Basis of Climate and Climate Modelling, edited by: Bolin, B., Doos, B., Godson, W., Hasselmann, K., Kutzbach, J., and Sawyer, J., World Meteorological Organisation, Geneva, 132136, 1975.Lucarini, V. and Ragone, F.: Energetics of climate models: net energy balance and meridional enthalpy transport, Rev. Geophys., 49, RG1001, doi:http://dx.doi.org/10.1029/2009RG000323
Web End =10.1029/2009RG000323 http://dx.doi.org/10.1029/2009RG000323
Web End = , 2011.
Lucarini, V. and Sarno, S.: A statistical mechanical approach for the computation of the climatic response to general forcings, Nonlin. Processes Geophys., 18, 728, doi:http://dx.doi.org/10.5194/npg-18-7-2011
Web End =10.5194/npg-18-7-2011 http://dx.doi.org/10.5194/npg-18-7-2011
Web End = , 2011.
Lucarini, V., Fraedrich, K., and Lunkeit, F.: Thermodynamics of climate change: generalized sensitivities, Atmos. Chem. Phys., 10, 97299737, doi:http://dx.doi.org/10.5194/acp-10-9729-2010
Web End =10.5194/acp-10-9729-2010 http://dx.doi.org/10.5194/acp-10-9729-2010
Web End = , 2010.
Lucarini, V., Blender, R., Herbert, C., Ragone, F., Pascale, S., and
Wouters, J.: Mathematical and physical ideas for climate science, Rev. Geophys., 52, 809859, doi:http://dx.doi.org/10.1002/2013RG000446
Web End =10.1002/2013RG000446 http://dx.doi.org/10.1002/2013RG000446
Web End = , 2014.Mah, G. and Paturel, J.-E.: 18962006 Sahelian annual rainfall variability and runoff increase of Sahelian Rivers, C. R. Geosci., 341, 538546, doi:http://dx.doi.org/10.1016/j.crte.2009.05.002
Web End =10.1016/j.crte.2009.05.002 http://dx.doi.org/10.1016/j.crte.2009.05.002
Web End = , 2009.
Mantua, N. J., Hare, S. R., Zhang, Y., Wallace, J. M., and Francis, R. C.: A Pacic interdecadal climate oscillation with impacts on salmon production, B. Am. Meteorol. Soc., 78, 10691079, 1997.
Maronna, R. and Yohai, V. J.: A bivariate test for the detection of a systematic change in mean, J. Am. Stat. Assoc., 73, 640645, 1978.
Marvel, K., Schmidt, G. A., Shindell, D., Bonls, C., LeGrande, A.N., Nazarenko, L., and Tsigaridis, K.: Do responses to different
anthropogenic forcings add linearly in climate models?, Environ. Res. Lett., 10, 104010, 2015.
Mayo, D. G.: Evidence as passing severe tests: highly probable versus highly probed hypotheses, in: Scientic Evidence: Philosophical Theories and Applications, edited by: Achinstein, P., John Hopkins University Press, Baltimore and London, 95127, 2005.
Mayo, D. G.: Learning from error, severe testing, and the growth of theoretical knowledge, in: Error and Inference: Recent Exchanges on Experimental Reasoning, Reliability, and the Objectivity and Rationality of Science, edited by: Mayo, D. G. and Spanos, A., Cambridge University Press, Cambridge UK and New York USA, 2857, 2010.
Mayo, D. G. and Cox, D.: Frequentist Statistics as a Theory of Inductive Inference, in: Error and Inference: Recent Exchanges on Experimental Reasoning, Reliability, and the Objectivity and Rationality of Science, edited by: Mayo, D. G. and Spanos,A., Cambridge University Press, Cambridge UK and New York USA, 247275, 2010.
Mayo, D. G. and Spanos, A.: Error and Inference: Recent Exchanges on Experimental Reasoning, Reliability, and the Objectivity and Rationality of science, Cambridge University Press, Cambridge UK, 419 pp., 2010.
Mears, C. A. and Wentz, F. J.: Construction of the RSS V3. 2 lower-tropospheric temperature dataset from the MSU and AMSU microwave sounders, J. Atmos. Ocean. Tech., 26, 14931509, 2009. Meehl, G. A.: Decadal climate variability and the early-2000s hiatus, US CLIVAR, 13, 16, 2015.
Meehl, G. A., Hu, A., and Santer, B. D.: The Mid-1970s Climate Shift in the Pacic and the Relative Roles of Forced versus Inherent Decadal Variability, J. Climate, 22, 780792, doi:http://dx.doi.org/10.1175/2008JCLI2552.1
Web End =10.1175/2008JCLI2552.1 http://dx.doi.org/10.1175/2008JCLI2552.1
Web End = , 2009.
Meehl, G. A., Hu, A., Arblaster, J. M., Fasullo, J., and Trenberth,K. E.: Externally Forced and Internally Generated Decadal Climate Variability Associated with the Interdecadal Pacic Oscillation, J. Climate, 26, 72987310, doi:http://dx.doi.org/10.1175/JCLI-D-12-00548.1
Web End =10.1175/JCLI-D-12- http://dx.doi.org/10.1175/JCLI-D-12-00548.1
Web End =00548.1 , 2013.
Meehl, G. A., Hu, A., and Teng, H.: Initialized decadal prediction for transition to positive phase of the Inter-decadal Pacic Oscillation, Nature Communications, 7, 11718, doi:http://dx.doi.org/10.1038/ncomms11718
Web End =10.1038/ncomms11718 http://dx.doi.org/10.1038/ncomms11718
Web End = , 2016.
Menberg, K., Blum, P., Kurylyk, B. L., and Bayer, P.: Observed groundwater temperature response to recent climate change, Hydrol. Earth Syst. Sci., 18, 44534466, doi:http://dx.doi.org/10.5194/hess-18-4453-2014
Web End =10.5194/hess-18-4453- http://dx.doi.org/10.5194/hess-18-4453-2014
Web End =2014 , 2014.
Miller, A. J., Cayan, D. R., Barnett, T. P., Graham, N. E., and Ober-huber, J. M.: The 197677 climate shift of the Pacic Ocean, Oceanography, 7, 2126, 1994.
Morice, C. P., Kennedy, J. J., Rayner, N. A., and Jones, P.D.: Quantifying uncertainties in global and regional temperature change using an ensemble of observational estimates: The HadCRUT4 data set, J. Geophys. Res.-Atmos., 117, D08101, doi:http://dx.doi.org/10.1029/2011JD017187
Web End =10.1029/2011JD017187 http://dx.doi.org/10.1029/2011JD017187
Web End = , 2012.
Morton, F. I.: Operational estimates of areal evapotranspiration and their signicance to the science and practice of hydrology, J. Hydrol., 66, 176, 1983.
Morton, F. I.: Practical estimates of lake evaporation, J. Clim. Appl.Meteorol., 25, 371387, 1986.
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming 209
Narisma, G. T., Foley, J. A., Licker, R., and Ramankutty, N.: Abrupt changes in rainfall during the twentieth century, Geophys. Res.Lett., 34, L06710, doi:http://dx.doi.org/10.1029/2006GL028628
Web End =10.1029/2006GL028628 http://dx.doi.org/10.1029/2006GL028628
Web End = , 2007.
Nicholls, N., Dellamarta, P., and Collins, D.: 20th century changes in temperature and rainfall in New South Wales, Aust. Meteorol.Mag., 53, 263268, 2004.
Normile, D.: El Nios warmth devastating reefs worldwide, Science, 352, 1516, doi:http://dx.doi.org/10.1126/science.352.6281.15
Web End =10.1126/science.352.6281.15 http://dx.doi.org/10.1126/science.352.6281.15
Web End = , 2016.
North, G. R., Kim, K.-Y., Shen, S. S. P., and Hardin, J.W.: Detection of Forced Climate Signals. Part 1: Filter Theory, J. Climate, 8, 401408, doi:http://dx.doi.org/10.1175/1520-0442(1995)008<0401:DOFCSP>2.0.CO; 2
Web End =10.1175/1520- http://dx.doi.org/10.1175/1520-0442(1995)008<0401:DOFCSP>2.0.CO; 2
Web End =0442(1995)008<0401:DOFCSP>2.0.CO; 2 , 1995.
North, R. P., Livingstone, D. M., Hari, R. E., Kster, O., Nieder-hauser, P., and Kipfer, R.: The physical impact of the late 1980s climate regime shift on Swiss rivers and lakes, Inland Waters, 3, 341350, 2013.
OKane, T. J., Matear, R. J., Chamberlain, M. A., and Oke, P. R.: ENSO regimes and the late 1970s climate shift: The role of synoptic weather and South Pacic ocean spiciness, J. Comput.Phys., 271, 1938, 2014.
Overland, J., Rodionov, S., Minobe, S., and Bond, N.: North Pacic regime shifts: Denitions, issues and recent transitions, Prog.Oceanogr., 77, 92102, 2008.
Ozawa, H., Ohmura, A., Lorenz, R. D., and Pujol, T.: The second law of thermodynamics and the global climate system: A review of the maximum entropy production principle, Rev. Geophys., 41, 1018, doi:http://dx.doi.org/10.1029/2002RG000113
Web End =10.1029/2002RG000113 http://dx.doi.org/10.1029/2002RG000113
Web End = , 2003.
Palmer, T. N.: A nonlinear dynamical perspective on climate change, Weather, 48, 314326, doi:http://dx.doi.org/10.1002/j.1477-8696.1993.tb05802.x
Web End =10.1002/j.1477- http://dx.doi.org/10.1002/j.1477-8696.1993.tb05802.x
Web End =8696.1993.tb05802.x , 1993.
Palmer, T. N.: A nonlinear dynamical perspective on climate prediction, J. Climate, 12, 575591, 1999.
Palmer, T. N., Doblas-Reyes, F. J., Weisheimer, A., and Rodwell, M.J.: Toward Seamless Prediction: Calibration of Climate Change Projections Using Seasonal Forecasts, B. Am. Meteorol. Soc., 89, 459470, doi:http://dx.doi.org/10.1175/bams-89-4-459
Web End =10.1175/bams-89-4-459 http://dx.doi.org/10.1175/bams-89-4-459
Web End = , 2008.
Peterson, T. C. and Vose, R. S.: An overview of the Global Historical Climatology Network temperature database, B. Am. Meteorol. Soc., 78, 28372849, 1997.
Peyser, C. E., Yin, J., Landerer, F. W., and Cole, J. E.: Pacic sea level rise patterns and global surface temperature variability, Geophys. Res. Lett., 43, 86628669, doi:http://dx.doi.org/10.1002/2016GL069401
Web End =10.1002/2016GL069401 http://dx.doi.org/10.1002/2016GL069401
Web End = , 2016.
Potter, K.: Illustration of a new test for detecting a shift in mean in precipitation series, Mon. Weather Rev., 109, 20402045, 1981.Power, S., Tseitkin, F., Torok, S., Lavery, B., and McAvaney, B.:
Australian temperature, Australian rainfall, and the Southern Oscillation, 19101996: coherent variability and recent changes, Aust. Meteorol. Mag., 47, 85101, 1998.
Power, S., Sadler, B., and Nicholls, N.: The inuence of climate science on water management in Western Australia Lessons for climate scientists, B. Am. Meteorol. Soc., 86, 839844, doi:http://dx.doi.org/10.1175/bams-86-6-839
Web End =10.1175/bams-86-6-839 http://dx.doi.org/10.1175/bams-86-6-839
Web End = , 2005.
Power, S., Delage, F., Wang, G., Smith, I., and Kociuba, G.: Apparent limitations in the ability of CMIP5 climate models to simulate recent multi-decadal change in surface temperature: implications for global temperature projections, Clim. Dynam., 117, doi:http://dx.doi.org/10.1007/s00382-016-3326-x
Web End =10.1007/s00382-016-3326-x http://dx.doi.org/10.1007/s00382-016-3326-x
Web End = , 2016.
Ragone, F., Lucarini, V., and Lunkeit, F.: A new framework for climate sensitivity and prediction: a modelling perspective, Clim. Dynam., 46, 14591471, 2016.
Ramanathan, V. and Coakley, J.: Climate modeling through radiative-convective models, Rev. Geophys. Space GE, 16, 465 489, 1978.
Ramaswamy, V., Boucher, O., Haigh, J., Hauglustaine, D., Haywood, J., Myhre, G., Nakajima, T., Shi, G. Y., and Solomon, S.: Radiative forcing of climate change, in: Climate Change 2001: The Scientic Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Houghton, J. T., Ding, Y., Griggs, D.J., Noguer, M., van der Linden, P. J., Dai, X., Maskell, K., and Johnson, C. A., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 349416, 2001.
Raper, S. C. B., Gregory, J. M., and Osborn, T. J.: Use of an upwelling-diffusion energy balance climate model to simulate and diagnose A/OGCM results, Clim. Dynam., 17, 601613, doi:http://dx.doi.org/10.1007/pl00007931
Web End =10.1007/pl00007931 http://dx.doi.org/10.1007/pl00007931
Web End = , 2001.
Raper, S. C. B., Gregory, J. M., and Stouffer, R. J.: The role of climate sensitivity and ocean heat uptake on AOGCM transient temperature response, J. Climate, 15, 124130, 2002.
Reeves, J., Chen, J., Wang, X. L., Lund, R., and Lu, Q. Q.: A Review and Comparison of Changepoint Detection Techniques for Climate Data, J. Appl. Meteorol. Clim., 46, 900915, doi:http://dx.doi.org/10.1175/JAM2493.1
Web End =10.1175/JAM2493.1 http://dx.doi.org/10.1175/JAM2493.1
Web End = , 2007.
Reid, P. C.: Ocean warming: setting the scene, in: Explaining ocean warming:causes, scale, effects and consequences, edited by: Laf-foley, D. and Baxter, J. M., IUCN, Gland, Switzerland, 2016. Reid, P. C. and Beaugrand, G.: Global synchrony of an accelerating rise in sea surface temperature, J. Mar. Biol. Assoc. UK, 92, 14351450, doi:http://dx.doi.org/10.1017/S0025315412000549
Web End =10.1017/S0025315412000549 http://dx.doi.org/10.1017/S0025315412000549
Web End = , 2012.
Reid, P. C., Hari, R. E., Beaugrand, G., Livingstone, D. M., Marty,C., Straile, D., Barichivich, J., Goberville, E., Adrian, R., and Aono, Y.: Global impacts of the 1980s regime shift, Glob. Change Biol., 22, 682703, doi:http://dx.doi.org/10.1111/gcb.13106
Web End =10.1111/gcb.13106 http://dx.doi.org/10.1111/gcb.13106
Web End = , 2016. Ricketts, J. H.: A probabilistic approach to climate regime shift detection based on Maronnas bivariate test, The 21st International Congress on Modelling and Simulation (MODSIM2015), Gold Coast, Queensland, Australia, 2015.
Rodionov, S. N.: A brief overview of the regime shift detection methods, Large-Scale Disturbances (Regime Shifts) and Recovery in Aquatic Ecosystems: Challenges for Management Toward Sustainability, UNESCO-ROSTE/BAS Workshop on Regime Shifts, Varna, Bulgaria, 1416 June 2005.
Roemmich, D. and Gilson, J.: The global ocean imprint of ENSO, Geophys. Res. Lett., 38, L13606, doi:http://dx.doi.org/10.1029/2011GL047992
Web End =10.1029/2011GL047992 http://dx.doi.org/10.1029/2011GL047992
Web End = , 2011.
Roemmich, D., Church, J., Gilson, J., Monselesan, D., Sutton, P., and Wijffels, S.: Unabated planetary warming and its ocean structure since 2006, Nature Climate Change, 5, 240245, doi:http://dx.doi.org/10.1038/nclimate2513
Web End =10.1038/nclimate2513 http://dx.doi.org/10.1038/nclimate2513
Web End = , 2015.
Rohde, R., Muller, R. A., Jacobsen, R., Muller, E., Perlmutter, S., Rosenfeld, A., Wurtele, J., Groom, D., and Wickham, C.: A new estimate of the average earth surface land temperature spanning 1753 to 2011, Geoinformatics & Geostatistics: An Overview, 1, 1000101, doi:http://dx.doi.org/10.4172/2327-4581.1000101
Web End =10.4172/2327-4581.1000101 http://dx.doi.org/10.4172/2327-4581.1000101
Web End = , 2012.
www.earth-syst-dynam.net/8/177/2017/ Earth Syst. Dynam., 8, 177210, 2017
210 R. N. Jones and J. H. Ricketts: Reconciling the signal and noise of atmospheric warming
Rypdal, M. and Rypdal, K.: Long-memory effects in linear response models of Earths temperature and implications for future global warming, J. Climate, 27, 52405258, 2014.
Sahin, S. and Cigizoglu, H. K.: Homogeneity analysis of Turkish meteorological data set, Hydrol. Process., 24, 981992, 2010.Santer, B. D., Mears, C., Doutriaux, C., Caldwell, P., Gleckler, P.J., Wigley, T. M. L., Solomon, S., Gillett, N. P., Ivanova, D.,
Karl, T. R., Lanzante, J. R., Meehl, G. A., Stott, P. A., Taylor,K. E., Thorne, P. W., Wehner, M. F., and Wentz, F. J.: Separating signal and noise in atmospheric temperature changes: The importance of timescale, J. Geophys. Res., 116, D22105, doi:http://dx.doi.org/10.1029/2011jd016263
Web End =10.1029/2011jd016263 http://dx.doi.org/10.1029/2011jd016263
Web End = , 2011.
Schmidt, G. A., Ruedy, R. A., Miller, R. L., and Lacis, A. A.: Attribution of the present-day total greenhouse effect, J. Geophys.Res.-Atmos., 115, D20106, doi:http://dx.doi.org/10.1029/2010JD014287
Web End =10.1029/2010JD014287 http://dx.doi.org/10.1029/2010JD014287
Web End = , 2010.Schmidt, G. A., Kelley, M., Nazarenko, L., Ruedy, R., Russell, G.L., Aleinov, I., Bauer, M., Bauer, S. E., Bhat, M. K., and Bleck,R.: Conguration and assessment of the GISS ModelE2 contributions to the CMIP5 archive, Journal of Advances in Modeling Earth Systems, 6, 141184, 2014a.
Schmidt, G. A., Shindell, D. T., and Tsigaridis, K.: Reconciling warming trends, Nat. Geosci., 7, 158160, doi:http://dx.doi.org/10.1038/ngeo2105
Web End =10.1038/ngeo2105 http://dx.doi.org/10.1038/ngeo2105
Web End = , 2014b.
Schneider, S. H. and Dickinson, R. E.: Climate modeling, Rev. Geophys., 12, 447493, 1974.
Seidel, D. J. and Lanzante, J. R.: An assessment of three alternatives to linear trends for characterizing global atmospheric temperature changes, J. Geophys. Res., 109, D14108, doi:http://dx.doi.org/10.1029/2003jd004414
Web End =10.1029/2003jd004414 http://dx.doi.org/10.1029/2003jd004414
Web End = , 2004.
Skeptical Science: The Escalator, available at: http://www.skepticalscience.com/graphics.php?g=47
Web End =http://www. http://www.skepticalscience.com/graphics.php?g=47
Web End =skepticalscience.com/graphics.php?g=47 (last access: 23 February 2016), John Cook, Brisbane, 2015.
Solomon, A., Goddard, L., Kumar, A., Carton, J., Deser, C., Fukumori, I., Greene, A. M., Hegerl, G., Kirtman, B., Kushnir, Y., Newman, M., Smith, D., Vimont, D., Delworth, T., Meehl, G. A., and Stockdale, T.: Distinguishing the Roles of Natural and Anthropogenically Forced Decadal Climate Variability, B. Am. Meteorol. Soc., 92, 141156, doi:http://dx.doi.org/10.1175/2010bams2962.1
Web End =10.1175/2010bams2962.1 http://dx.doi.org/10.1175/2010bams2962.1
Web End = , 2011.Street, R. B.: Towards a leading role on climate services in Europe:
A research and innovation roadmap, Climate Services, 1, 25, 2016.
Suppes, P.: Models of data, in: Studies in the Methodology and Foundations of Science, edited by: Nagel, E., Suppes, P., and Tarski, A., Stanford University Press, Stanford CA, 252261, 1962.
Swanson, K. L., Sugihara, G., and Tsonis, A. A.: Long-term natural variability and 20th century climate change, P. Natl. Acad. Sci., 106, 1612016123, doi:http://dx.doi.org/10.1073/pnas.0908699106
Web End =10.1073/pnas.0908699106 http://dx.doi.org/10.1073/pnas.0908699106
Web End = , 2009.Timbal, B., Arblaster, J. M., and Power, S.: Attribution of the late-twentieth-century rainfall decline in southwest Australia, J. Climate, 19, 20462062, 2006.
Timbal, B., Arblaster, J., Braganza, K., Fernandez, E., Hendon, H., Murphy, B., Raupach, M., Rakich, C., Smith, I., Whan, K., and Wheeler, M.: Understanding the anthropogenic nature of the ob-served rainfall decline across south-eastern Australia, The Centre for Australian Weather and Climate Research, Melbourne, 180 pp., 2010.
Trenberth, K. E.: Changes in precipitation with climate change, Climate Research, 47, 123138, 2011.
Trenberth, K. E.: Has there been a hiatus?, Science, 349, 691692, doi:http://dx.doi.org/10.1126/science.aac9225
Web End =10.1126/science.aac9225 http://dx.doi.org/10.1126/science.aac9225
Web End = , 2015.
Varotsos, C. A., Franzke, C. L., Efstathiou, M. N., and Degermendzhi, A. G.: Evidence for two abrupt warming events of SST in the last century, Theor. Appl. Climatol., 116, 5160, 2014.
Vivs, B. and Jones, R. N.: Detection of Abrupt Changes in Australian Decadal Rainfall (18901989), CSIRO Atmospheric Research, Melbourne, 54, 2005.
Wack, P.: The Gentle Art of Reperceiving Scenarios: Shooting the Rapids (part 2 of a two-part article), Harvard Business Review, NovemberDecember, 214, 1985a.
Wack, P.: The Gentle Art of Reperceiving Scenarios: Uncharted Waters Ahead (part 1 of a two-part article), Harvard Business Review, SeptemberOctober, 7389, 1985b.
Watanabe, M., Shiogama, H., Tatebe, H., Hayashi, M., Ishii, M., and Kimoto, M.: Contribution of natural decadal variability to global warming acceleration and hiatus, Nature Climate Change, 4, 893897, doi:http://dx.doi.org/10.1038/nclimate2355
Web End =10.1038/nclimate2355 http://dx.doi.org/10.1038/nclimate2355
Web End = , 2014.
Werner, R., Valev, D., Danov, D., Guineva, V., and Kirillov, A.: Analysis of global and hemispheric temperature records and prognosis, Adv. Space Res., 55, 29612973, doi:http://dx.doi.org/10.1016/j.asr.2015.03.005
Web End =10.1016/j.asr.2015.03.005 http://dx.doi.org/10.1016/j.asr.2015.03.005
Web End = , 2015.
Whetton, P., Adamson, D., and Williams, M.: Rainfall and river ow variability in Africa, Australia and East Asia linked to El Nio-Southern Oscillation events, Geological Society of Australia Symposium Proceedings, 25 August 1988, 7182, 1990. White, H.: A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity, Econometrica: Journal of the Econometric Society, 48, 817838, 1980.
White, N. J., Haigh, I. D., Church, J. A., Koen, T., Watson, C. S.,
Pritchard, T. R., Watson, P. J., Burgette, R. J., McInnes, K. L., You, Z.-J., Zhang, X., and Tregoning, P.: Australian sea levels Trends, regional variability and inuencing factors, Earth-Sci. Rev., 136, 155174, doi:http://dx.doi.org/10.1016/j.earscirev.2014.05.011
Web End =10.1016/j.earscirev.2014.05.011 http://dx.doi.org/10.1016/j.earscirev.2014.05.011
Web End = , 2014. Wolter, K. and Timlin, M. S.: El Nio/Southern Oscillation behaviour since 1871 as diagnosed in an extended multivariate ENSO index (MEI.ext), Int. J. Climatol., 31, 10741087, doi:http://dx.doi.org/10.1002/joc.2336
Web End =10.1002/joc.2336 http://dx.doi.org/10.1002/joc.2336
Web End = , 2011.
World Meteorological Organization: Position Paper on Global Framework for Climate Services, World Meteorological Organisation, Geneva, 2010.
World Meteorological Organization: Climate Knowledge for Action: A Global Framework for Climate Services Empowering the Most Vulnerable, World Meteorological Organization, Geneva, 247 pp., 2011.
Wu, Z., Huang, N., Wallace, J., Smoliak, B., and Chen, X.: On the time-varying trend in global-mean surface temperature, Clim. Dynam., 37, 759773, doi:http://dx.doi.org/10.1007/s00382-011-1128-8
Web End =10.1007/s00382-011-1128-8 http://dx.doi.org/10.1007/s00382-011-1128-8
Web End = , 2011. Yao, S.-L., Huang, G., Wu, R.-G., and Qu, X.: The global warming hiatus a natural product of interactions of a secular warming trend and a multi-decadal oscillation, Theor. Appl. Climatol., 123, 349360, doi:http://dx.doi.org/10.1007/s00704-014-1358-x
Web End =10.1007/s00704-014-1358-x http://dx.doi.org/10.1007/s00704-014-1358-x
Web End = , 2016.
Zhou, J. and Tung, K.-K.: Deducing Multidecadal Anthropogenic Global Warming Trends Using Multiple Regression Analysis, J. Atmos. Sci., 70, 38, doi:http://dx.doi.org/10.1175/JAS-D-12-0208.1
Web End =10.1175/JAS-D-12-0208.1 http://dx.doi.org/10.1175/JAS-D-12-0208.1
Web End = , 2013.
Earth Syst. Dynam., 8, 177210, 2017 www.earth-syst-dynam.net/8/177/2017/
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright Copernicus GmbH 2017
Abstract
Interactions between externally forced and internally generated climate variations on decadal timescales is a major determinant of changing climate risk. Severe testing is applied to observed global and regional surface and satellite temperatures and modelled surface temperatures to determine whether these interactions are independent, as in the traditional signal-to-noise model, or whether they interact, resulting in step-like warming. The multistep bivariate test is used to detect step changes in temperature data. The resulting data are then subject to six tests designed to distinguish between the two statistical hypotheses, h<sub>step</sub> and h<sub>trend</sub>. Test 1: since the mid-20th century, most observed warming has taken place in four events: in 1979/80 and 1997/98 at the global scale, 1988/89 in the Northern Hemisphere and 1968-70 in the Southern Hemisphere. Temperature is more step-like than trend-like on a regional basis. Satellite temperature is more step-like than surface temperature. Warming from internal trends is less than 40% of the total for four of five global records tested (1880-2013/14). Test 2: correlations between step-change frequency in observations and models (1880-2005) are 0.32 (CMIP3) and 0.34 (CMIP5). For the period 1950-2005, grouping selected events (1963/64, 1968-70, 1976/77, 1979/80, 1987/88 and 1996-98), the correlation increases to 0.78. Test 3: steps and shifts (steps minus internal trends) from a 107-member climate model ensemble (2006-2095) explain total warming and equilibrium climate sensitivity better than internal trends. Test 4: in three regions tested, the change between stationary and non-stationary temperatures is step-like and attributable to external forcing. Test 5: step-like changes are also present in tide gauge observations, rainfall, ocean heat content and related variables. Test 6: across a selection of tests, a simple stepladder model better represents the internal structures of warming than a simple trend, providing strong evidence that the climate system is exhibiting complex system behaviour on decadal timescales. This model indicates that in situ warming of the atmosphere does not occur; instead, a store-and-release mechanism from the ocean to the atmosphere is proposed. It is physically plausible and theoretically sound. The presence of step-like - rather than gradual - warming is important information for characterising and managing future climate risk.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer