Extreme temperature and precipitation events occurring in Australia in recent decades have caused significant socio-economic and environmental impacts, and thus determining the factors contributing to these extremes is an active area of research. Many recently occurring record-breaking temperature and rainfall events have now been examined from an extreme event attribution (EEA) perspective. This paper describes a set of studies that have examined the causes of extreme climate events using various general circulation models (GCMs), presenting a comprehensive methodology for GCM-based attribution of extremes of temperature and precipitation observed on large spatial and temporal scales in Australia. First, we review how Coupled Model Intercomparison Project Phase 5 (CMIP5) models have been used to examine the changing odds of observed extremes. Second, we review how a large perturbed initial condition ensemble of a single climate model (CESM) has been used to quantitatively examine the changing characteristics of Australian heat extremes. For each approach, methodological details and applications are provided and limitations highlighted. The conclusions of this methodological review discuss the limitations and uncertainties associated with this approach and identify key unexplored applications of GCM-based attribution of extremes. Ideally, this information will be useful for the application of the described extreme event attribution approaches elsewhere.
Since around the year 2000, many high-impact extreme and record-breaking temperature and precipitation events have occurred over Australian regions (Bureau of Meteorology, 2013a, b, c, 2016a). During 2001–2010, the Murray–Darling Basin (MDB) region, in south-eastern Australia, experienced an extreme dry period that had severe environmental and socio-economic consequences (Van Dijk et al., 2013). The Millennium Drought in the MDB (defined here as the period from 2001 to 2010) was the longest uninterrupted series of years since at least 1900 with below median observed rainfall. The persistent heavy rainfall of 2010–2012 (Bureau of Meteorology, 2012) ended the drought of the preceding decade but resulted in severe flooding in northern and eastern Australia. This led to billions of dollars of flood damage, the evacuation of thousands and 33 deaths (Queensland Floods Commission of Inquiry, 2012). Furthermore, following the consecutive strong La Niña events of 2010 to 2012, a period of well above-average temperatures occurred across Australia (Bureau of Meteorology, 2014). In the period since, record-breaking extremes were observed on sub-daily to multi-year timescales and on spatial scales ranging from single stations to Australia-wide (Bureau of Meteorology, 2016b). Many of these events were the most extreme since high-quality national records began in 1910.
Australia has a highly variable climate that is influenced by several large-scale climate modes such as El Niño–Southern Oscillation (ENSO), the Interdecadal Pacific Oscillation (IPO), the Indian Ocean Dipole (IOD), the Southern Annular Mode (SAM), the Subtropical Ridge (STR), and the Madden Julian Oscillation (MJO) (King et al., 2014; Maher and Sherwood, 2014; Min et al., 2013). These large-scale modes may also interact with local- and regional-scale processes, such as soil moisture feedbacks, to impact the severity, duration or likelihood of extreme events such as heatwaves (Perkins et al., 2015). Anthropogenic climate change is also a key influence in the observed characteristics of some extreme weather and climate events (Herring et al., 2018). Understanding the influences of both these anthropogenic forcings and/or natural climatic variability on recent temperature and rainfall extremes affecting Australia, as well as elsewhere, has become an active research avenue.
Extreme event attribution (EEA) studies focus on understanding a particular observed extreme weather or climate event. These studies typically combine observational and model data to determine whether various factors (e.g. anthropogenic greenhouse gas composition changes) contributed to a specific observed aspect of an observed extreme event, such as its intensity, magnitude or frequency. Likelihood-based EEA approaches compare the probability of occurrence of an extreme event in the current climate with the occurrence in a counterfactual climate without anthropogenic climate change. There are a number of different methods possible for applying the risk-based EEA approach, including using observations, regional climate models, or coupled or atmosphere-only general circulation models (GCMs).
Numerous previous studies have made quantitative assessments of the
influences of anthropogenic forcings on the likelihood of an observed event
using fraction of attributable risk (FAR) values. A FAR value is usually
defined as
This paper contributes to a series that outlines different methodological approaches to examining and applying event attribution approaches to observed extremes. Many extreme weather and climate events occurring in Australia have now been examined from an event attribution perspective using a cognate methodology. Broadly, this methodology is a quantitative event attribution approach using GCMs and some Earth system models (ESMs), which is outlined in detail here. As a large number of studies have now been conducted focused on Australian weather and climate extremes (see Lewis et al., 2017b), we detail the methodological approaches used in these analyses. In this paper, the climate models and attribution approaches that have been applied to Australia's observed extremes are described and compared with a level of methodological detail not provided previously. The specific purpose of this contribution to this EEA-focused series is (i) to provide detailed information about model-based EEA studies applied to Australia and hence (ii) to provide guidance for others aiming to interpret these previously published results or ideally apply the described EEA approaches elsewhere.
The Australian-focused studies detailed in our paper use differing EEA model frameworks that we call either multi-model or single-model ensembles, with the applications, advantages and limitations of each discussed throughout. The multi-model ensemble is constituted by many different climate models, each with a small number of contributing ensemble members or realizations. The single-model ensemble is alternatively constituted by a large number of ensemble members of one model only that differ in terms of model physics or initial conditions in order to sample internal climate variability. Each model framework provides a different level of conditioning, which provides a methodological focus on specific aspects of weather and climate such as sea surface conditions. This degree of conditioning is important for interpreting attribution statements (National Academies of Sciences, Engineering, and Medicine, 2016) and is discussed throughout.
In the first examples, attribution statements on observed extremes were made by analysing models participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) (Taylor et al., 2012). CMIP5 provides a framework for coordinated and standardized climate change experiments, including detection and attribution experiments. More than 50 models, from over 20 international groups, participated in CMIP5 by contributing data, although the number of models and realizations varies depending on experiments. Experiments are tiered by importance and participation for each experiment is voluntary. All the models are fully coupled (atmosphere–ocean GCMs and ESMs) and of varying horizontal and vertical resolutions, and are run under standardized time-evolving forcings (e.g. time-varying concentrations of various atmospheric constituents such as greenhouse gases).
Multiple CMIP5 experiments can be used together to analyse changes in the
characteristics of observed weather and climate extremes that are associated
with various forcings. These experiments are
historical, 1850–2005, using time-evolving atmospheric compositions due to
observed anthropogenic and volcanic influences, solar forcings and emissions
of short-lived species from natural and anthropogenic aerosols. historicalNat, 1850–2005, using time-evolving volcanic influences and solar
forcings. piControl. Control runs provided for each model which are long, freely
evolving climate simulations with greenhouse gas concentrations appropriate
for circa 1850 that permit the analysis of a large number of model years. rcp8.5. Representative Concentration Pathway (RCP) simulations of
2006–2100. RCPs are scenarios that assume particular policies are
implemented to achieve greenhouse gas emissions targets. These scenarios
include a range of future projections for future populations and
technological responses (Moss et al.,
2010). For instance, in the rcp8.5 future scenario, the radiative forcing
increases before reaching a level of about 8.5 W m
We note that the use of multiple CMIP5 models does not condition to aspects
of the observed state of the climate (such as sea surface temperatures) that
prevailed during the event in question. Simulated SSTs in a given year in
CMIP5 realizations do not relate directly to those observed. This
unconditional approach using CMIP5 is described as the most comprehensive
and easy to interpret (National Academies of Sciences,
Engineering, and Medicine, 2016, p. 51), as all aspects of weather and
climate are considered and analysis is not limited to particular climatic
situations. Coupled model ensembles may be used with a degree of
conditioning on climate modes by subsetting climate model years
characterized by a given mode of variability (e.g. the NINO3.4 index
could be calculated in models to discriminate between phases of the El
Niño–Southern Oscillation). This approach has been used in EEA studies,
but the conditioning is still, by necessity, less restricted than for
atmosphere-only simulations (discussed further in Sect. 2.3).
We next describe a single-model ensemble for EEA, which samples the
influence of internal variability in the attribution of climate extremes to
human influence. In recent years, this approach has been demonstrated using
two versions of the Community Earth System Model (CESM). The first version
is a 21-member ensemble of CESM-CAM4, run at approximately 2
The specifics of the multi-member CESM ensemble differ slightly compared to CMIP5. Firstly, instead of a historicalNat simulation, a long control is used to represent the counterfactual world. Here, greenhouse gas forcings are held to pre-industrial levels, and the model freely evolves under these constraints for the length of the simulation. Secondly, all realizations forced by historical and anthropogenic emissions are identical, except in a tiny atmospheric perturbation in the atmospheric temperature initial conditions of each simulation. This is enough to trigger completely different realizations of internal variability over the length of the simulation. The start of each realization is forced by historical conditions, and the rcp8.5 scenario from the start of 2006. The older version of CESM employed in attribution studies has a control run of 982 years and 21 anthropogenically forced realizations that commence in 1950. This compares to the newer version with a control run of 1800 years and 35 anthropogenically forced realizations that commence in 1920. The CESM ensemble provides an unconditional attribution framing.
We refer later to CMIP5 comparisons with results that are also derived from other model frameworks, which we describe here. First, we describe the Attribution of extreme weather and Climate Events (ACE) framework (Stott et al., 2016). Unlike the coupled ocean–atmosphere CMIP5 experiments, these are atmosphere-only model simulations driven by prescribed sea surface temperatures (SSTs). This model approach uses two large ensembles of simulations over the extreme event period: one ensemble simulates the actual climate, including all known external forcings (anthropogenic forcings are long-lived greenhouse gases, aerosols, tropospheric and stratospheric ozone and land-use changes, and natural forcings are from volcanic aerosols and solar irradiance). The second ensemble provides various possible representations of a “natural” climate without or with only minor human influences. In these experiments, an estimate of the anthropogenic contribution to observed SST conditions is calculated and removed. Both ensembles are comprised of members that are run with slightly different (perturbed) initial conditions. In these ACE experiments, the prescribed SSTs imposed mean that the attribution statement is conditioned on the current climate state.
We also refer to HAPPI (Half a Degree Additional warming, Prognosis, and
Projected Impacts) simulations (Mitchell et al., 2017).
These are decade long simulations of various scenarios, including the
present day (2006–20015) driven by observed SSTs and sea ice, and a
historicalNat simulation and stabilized at higher warming thresholds
(1.5 and 2
We next discuss the specifics of how these model datasets have been used in previous studies to determine attribution statements of observed Australian extremes, though note that other methodological approaches are also suitable for EEA using these models. We first detail a set of CMIP5-based attribution analyses focused on quantifying changes in the risk of Australian seasonal-scale events that were the most extreme in the observational record. Second, we discuss single-model attribution approaches. The suitability and application of these methods to other events and regions are discussed in Sect. 5.
Different CMIP5 models were included in ensembles used to analyse observed Australian extremes, based on data availability at the time of publication and the ability of models to capture critical aspects of observed climate. In the CMIP5-based attribution examples discussed in Sect. 4, a subset of available models was used that ideally facilitates the determination of meaningful attribution statements about causes of changes in the odds of an event. That is, models ideally capture the right mechanisms, with realistic frequency and characteristics, necessary to reproduce the event under consideration (Christidis et al., 2013a), and a sufficiently large ensemble of numerous models was available for analysis. Several approaches have been applied for assessing the validity of model simulations for Australian seasonal-scale extremes, although note that there are numerous evaluation procedures we do not discuss here.
In one example examining Australia-wide summer temperatures
(Lewis and Karoly, 2013a), models were selected
based on their skill in capturing interannual variability in Australian
summer temperatures. First, a time series of Australia areal-average observed
summer temperatures was calculated from the Australian Water Availability
Project (AWAP) gridded dataset from 1910 to the present (Jones et
al., 2009). Using a bootstrap resampling procedure with replacement, 2000
time series were synthesized from 10-year blocks of observed data. For each
synthetic time series, a standard deviation was calculated and a spread of
standard deviations determined. This resampling procedure was then applied
to all available CMIP5 historical simulations to calculate simulated
standard deviations of Australian summer temperatures for model years
1911–2005 (the period of modelled–observed overlap). In instances where models
contributed more than one realization to a CMIP5 experiment (e.g. r1i1p1,
r2i1p1), data were assessed collectively for all available
realizations. Next, this study compared observed and simulated data using a
Perkins skill score (Perkins et al., 2007).
The Perkins score is defined as
Conceptual representation of approaches to evaluating CMIP5
models, showing a visual comparison of a high
Further Australia EEA studies used a two-sided Kolmogorov–Smirnov test (KS)
to compare modelled distributions of variables with observed ones. For example, a
study of the consecutive record-breaking spring temperatures experienced in
Australian in 2013/2014 used KS tests to compare temperatures in CMIP5
historical simulations from various models to observed Australian annual
average temperatures. Only models where the distributions (estimated using a
kernel density function) were statistically indistinguishable (
These CMIP5 model selection steps were modified in subsequent Australia EEA studies, based on the event analysed. Analysis of the 2010–2012 eastern Australian extreme rainfall required that models capture several elements of observed climate variability (Lewis and Karoly, 2014a). Here, a larger set of evaluation criteria was applied to CMIP5 models to investigate this extreme event with a more complex climatological context than previous studies focused on large-scale long-duration temperature extremes (see Fig. 2). In this case, CMIP5 models were selected for analysis based on (1) their representation of monthly surface air temperature variability in the NINO3.4 region, (2) their representation of Australian rainfall variability and (3) rainfall amount associated with ENSO conditions (e.g. teleconnected relationship). This study considered simulated values to be physically realistic where values were in the 5th–95th percentile window of observed values determined using a bootstrap resampling method, whereby 10 000 time series of 80-year length were generated from observed anomalies. Models were only included where they are realistic compared to observation for all three criteria.
Schematic of the three-step CMIP5 model evaluation used for examining influences on heavy rainfall of 2010–2012. Left-hand side shows flowchart evaluation steps (for the NINO3.4 index and south-eastern Australian (SEA) rainfall in December–February (DJF)) and the number of models suitable, and right-hand side demonstrates evaluation of model realizations against observations for the metrics considered.
Australian EEA studies using CMIP5 models have employed FAR values as a means to quantify the change in probability of an observed event that can be assigned to anthropogenic climate forcings (e.g. greenhouse gases) (see Lewis and Karoly, 2013, 2014b). The calculation of a FAR value first requires that an event be precisely defined and a threshold of exceedance determined. Various studies have defined extreme events differently from the examples highlighted here and have also demonstrated that attribution statements can be sensitive to the spatio-temporal scales used to define the event (Angélil et al., 2014; Cattiaux and Ribes, 2018; Uhe et al., 2016). Hence, each application of this EEA approach must consider event definitions and the potential sensitivity of FAR values to definitions.
Seasonal-scale Australia temperature analyses defined a threshold of the
second most extreme entry in the observed record
(Lewis and Karoly, 2013a, 2014b). This
temperature anomaly (
FAR values compare the probability of an event exceeding the pre-defined
threshold (e.g.
The probabilities in Eq. (1) are determined through several possible approaches. The count-based estimate simply determines the number of times the defined threshold was exceeded, relative to the total sample size. Alternatively, probabilities can be calculated as the area under the fitted distribution exceeding the event threshold, compared to the entire distribution (Lewis and Karoly, 2013b). Distributions can be fitted using a kernel density function or a fitted generalized extreme value (GEV) model which is suited to examining rare events.
Studies also attempt to assess uncertainty in FAR calculations. For each experiment, only a single FAR value is obtained for each realization of each model. A bootstrap resampling procedure was applied in this set of CMIP5-based studies to evaluate the uncertainty associated with FAR estimated. In determining the FAR values associated with the 2012/2013 Australian summer temperatures, each distribution of temperature was bootstrap resampled 10 000 times (using in each iteration sub-samples of all years from only 50 % of available model simulations) and a distribution of FAR values was then calculated (Lewis and Karoly, 2013b). This calculated distribution of 10 000 FAR values represents the uncertainty associated with using different models and provides a basis for communicating FAR ranges. This 2012/2013 study, for example, reported both the median and 10th percentile FAR values, meaning they are exceeded by 90 % of the values in the bootstrapped FAR distributions. These are described respectively as “best estimate” and “very likely” values. By providing both the best estimate (median) and very likely (10th percentile) values, an estimate of the uncertainty in FAR values is conveyed.
Perkins et al. (2014) used the 21-member version of CESM to determine the anthropogenic signal behind heatwave intensity and frequency during the 2012/2013 Australian summer. The total number of heatwave days (heatwave frequency), as well as the hottest heatwave day (heatwave intensity) (see Perkins and Alexander, 2013), were computed using the excess heat factor heatwave metric (Nairn and Fawcett, 2013). This metric was used since it underpins operational heatwave forecasts in Australia. Areally averaged heatwave frequency and intensity for the 2012/2013 austral summer were extracted from AWAP (Jones et al., 2009). The frequency of each metric was compared via the forced and controlled simulations, as outlined above. Over the period 1984–2012, heatwave intensity increased by 2-fold, and frequency by 3-fold, relative to the counterfactual climate.
The calculations of FAR using a multi-member single-model ensemble are similar to those described for CMIP5. That is, evaluations based on the similarity between the observed and model distributions, such as the Kolmogorov–Smirnov test, are employed, however for each individual realization against the observations. Since Australian climate extremes occurring after 2005 have been analysed by this approach, the historical and rcp8.5 realizations were spliced together. Previous attribution studies using CESM have found no significant difference between each CESM realization and the observational dataset that includes the extreme event of interest (Perkins-Kirkpatrick et al., 2019a; Perkins et al., 2014; Perkins and Gibson, 2015). A more detailed evaluation of the weather systems that drive Australian heatwaves is provided by Perkins and Gibson (2015), with the model performing adequately across all 21 members.
After evaluation, bootstrapping with replacement is performed 10 000 times to create new samples of the counterfactual and factual simulations, each of which are 50 % the size of the original. This allows for uncertainties in model sampling. For each bootstrapped iteration, the observed extreme of interest is isolated and its frequency compared in the control and forced samples, resulting in a distribution of 10 000 FAR values. As a conservative measure, the 10th percentile FAR value is reported. This allows for 90 % confidence that the true FAR value based on CESM is at least as big as this value.
It is notable that small atmospheric perturbations initiate substantial
differences in the resolving climate variability (see Fig. 3), and,
therefore, within the extremes themselves. For example, Fischer et al. (2013) demonstrated that the spatially aggregated distribution of changes in
numerous extremes in CESM is comparable to that in CMIP5. Ranges in
regional heatwave trends also significantly vary across CESM realizations
(Perkins and Fischer, 2013). The representation of
synoptic systems, which ultimately drive extreme events over short
timescales, also demonstrates striking differences among CESM realizations
(Schaller et al., 2018). Thus, using just one
or a few realizations from a given model may not be an accurate estimate of
Example of the different realizations of ENSO, as defined by the NINO3.4 index, in a version of CESM. For comparison, a subset of only 5 of 35 available realizations is shown using different colours which demonstrate vastly different ENSO conditions in realizations differing based on perturbed initial atmospheric conditions.
We next detail examples that have employed the methods described above to examine observed weather and climate extremes in Australia.
The first study employing CMIP5 detection and attribution experiments to
quantitatively assess the relative influence of different climatic forcings
on an observed extreme was focused on the record hot Australia-wide
2012/2013 summer (Lewis and Karoly, 2013a). The
Australia-wide mean temperatures from December 2012 to February 2013 were the
hottest in the observed high-quality record (extending to 1910), with
average conditions exceeding the observed 1911–1940 mean by 1.32
By using FAR values to compare the likelihood of extreme Australia-wide
summer temperatures between the CMIP5 experiments, this analysis showed that
it was very likely (
Quantitative attribution using CMIP5 models also contributed to a later
study analysing factors affecting Great Barrier Reef (GBR) bleaching
(Lewis and Mallela, 2018). In 2016, the GBR endured a
significant bleaching event: 93 % of the northern, 700 km stretch of coral
was bleached, and by June,
For each metric, probabilities of exceedance were compared in CMIP5
detection and attribution experiments. Results varied depending on the
metric. There was a significant increase in the likelihood of extreme
January–March temperature anomalies when anthropogenic forcings are included
(median FAR
The extreme heavy precipitation of 2010–2012 occurring over Australia was also examined from an attribution perspective (Lewis and Karoly, 2014a). Over this period, eastern Australia experienced its heaviest ever 2-year accumulated rainfall. Rainfall records were broken on daily through to seasonal timescales for large spatial regions. The heavy rainfall and disastrous accompanying flooding were coincident with two strong, consecutive La Niña events, which are typically associated with enhanced rainfall in the eastern Australian rainfall. Using five CMIP5 climate models, King et al. (2013) examined the anthropogenic influence on rainfall totals and maximum consecutive 5 d rainfall (Rx5day) across south-eastern Australia in 2012. This analysis determined that there was little robust change in the risk of extreme rainfall events between the 1861–1890 and 1976–2005 periods. Overall, the magnitude of cool sea surface temperatures in the NINO3.4 region was found to have a greater effect on Rx5day compared to the magnitude of anomalously warm local SSTs.
A further study by Lewis and Karoly (2014a) examined this period of record heavy precipitation using CMIP5 models in combination with two sets of simulations conducted as part of the ACE initiative, in which ensembles of simulations are produced representing the recent climate with, and then without, the effects of human influences. This study explored the heavy rainfall in several defined Australia regions on multiple seasonal to multi-year timescales, and results were broadly in agreement with King et al. (2013). The approach used aimed to determine the robustness of FAR values for the 2010–2012 heavy Australia rainfall to changes in the attribution framework. Results showed that attribution statements for the rainfall events were variable, with FAR values sensitive to the attribution parameters considered, including event thresholds, regions and seasons. Furthermore, estimates of the attributable change in rainfall risk depended on the model datasets considered. This study argued that consideration of model outputs from several datasets (e.g. CMIP5 and ACE) was useful for establishing robust attribution statements for extreme rainfall events.
Analysis of the 2014 Australian May heatwave demonstrated a 23-fold increase
in experiencing an event of similar concurrent length and magnitude over the
period 1995–2020 (Perkins and Gibson, 2015). This study was
focused on a 19 d heatwave during 8–26 May 2014,
areally averaged for Australia, that had a magnitude of
In analysing the 2017/2018 Tasman Sea marine heatwave, the 35-member ensemble
of CESM was used in conjunction with CMIP5 to compare the two approaches
(Perkins-Kirkpatrick et al., 2019a). This was
to address whether different amounts of anthropogenic warming were detected
via the two ensembles. During November 2017–January 2018, sea surface
temperatures (SSTs) were 1.7
When the Paris Agreement was drafted in December 2015, with the aim of
limiting global warming to well below 2
Here we briefly discuss the “time-slicing” techniques used to analyse
extremes in future warmer worlds using a multi-model ensemble such as CMIP5.
The CMIP5 ensemble, and in particular the simulations run under different
representative concentration pathways, are designed to investigate transient
climate change projections under plausible greenhouse gas emissions
scenarios. Careful application of a “time-slicing method”
(James et al., 2017) allows for analysis of
1.5 and 2
In the CMIP5 models, a pre-industrial climate baseline may be defined from a
choice of experiments including the control simulations or historical
simulations with anthropogenic influences removed. Subsequently, time slices
from individual climate model projections can be extracted (see Fig. 4 for
a schematic). For example, King et al. (2017)
extracted all model years within decades where the global-average
temperature was 1.3–1.7 the use of all RCPs or just a subset of scenarios, the definition of a model pre-industrial baseline, and the temperature range and length of the averaging window.
For (1), if the difference that is being inferred between the 1.5
and 2
Schematic illustrating the time-slice method. Global-average annual
temperature anomalies (relative to 1901–2005 historicalNat average) are
shown for ACCESS1.3 rcp8.5 and rcp4.5 for 2006–2100. The boxes show the
years which are extracted from each RCP to construct the 1.5
This framework for analysis was then applied in King et al. (2017) to examine Australian extremes under the
Paris Agreement global warming limits. They extended existing attribution
analyses for events like the record hot “Angry summer” of 2012/2013
(Lewis and Karoly, 2013) to examine their
likelihoods under future warming scenarios. Consistent with Lewis and Karoly (2013), King et al. (2017) found a large increase in the likelihood of hot
Australian summers exceeding the previous record, from
This framework illustrates the utility of the CMIP5 ensemble and event
attribution methodology beyond examining extreme events only in the current
climate context. Analysis of the 1.5 and 2
The CMIP5-based event attribution methodology described above is useful for estimating the roles of anthropogenic and natural climate influences in extreme events. By making comparisons between model simulations representing the current climate and that of a counterfactual world without anthropogenic influences, the effect humans have had on extreme events may be quantified. Overall, the presented method is limited in its applications for several key reasons.
First, the CMIP5-based approaches are best suited to large spatio-temporal events (sub-continental to continental spatial-scale extreme events on daily to multi-year timescales), rather than short-duration or small-scale events that likely require large ensemble sizes and/or finer horizontal resolution. Second, the CMIP5 method is useful in sampling the anthropogenic attributable signal across varying climate model configurations. This means that differences in model physics, parameterization schemes, the simulation/exclusion of certain physical mechanisms that drive extreme events, as well as climate model sensitivity, are accounted for when quantifying FAR values. However, most models employed using this method are represented by just one or at most a few realizations of the same experiment (e.g. historical, historicalNat, rcp8.5). Therefore, the influence of internal variability is unlikely fully accounted for, especially at the individual model level.
Using a large ensemble of a single model in attribution assessments should
not currently be seen as a substitute to other approaches, particularly when
this type of ensemble is limited to, at best, a handful of individual
models. Rather, this approach is complementary to existing methods, by
including a robust anthropogenic estimate from plausible realizations of
fully coupled internal climate variability. For this to be effective,
evaluation of the model
This paper has outlined a cogent examination of recent extreme temperature and precipitation events in Australia from an extreme event attribution perspective using GCMs (and some ESMs). The models used in these examples are not conditioned to specific atmospheric or oceanic states and hence data are available and accessible for use. The methodologies can be applied to developing attribution statements for extreme events occurring elsewhere, provided that event-specific model evaluation is performed first. For precipitation-related extremes in particular, we recommend that rigorous model evaluation is necessary prior to attribution and that multiple model-based approaches are combined (for example, the combination of GCM- and AMIP-based results). In some extreme event cases, conditional EEA model frameworks that are conditioned by specific boundaries (SSTs or seasonal forecasts) may be more robust (National Academies of Sciences, Engineering, and Medicine, 2016).
There are some important overall limitations of the methods we have presented here. Firstly, while these studies provide robust evaluations of models used against observed climates, they typically employ singular analytical approaches. That is, one set of models (CMIP5 or CESM) was used to examine the characteristics of observed extremes for most of the reviewed studies, which may affect the attribution statement provided. Independent re-examination of other extreme event attribution studies found that attribution statements are broadly similar for extreme temperature events upon re-evaluation, but disagree for many extreme precipitation events (Angélil et al., 2017). Similarly, the Lewis and Karoly (2014a) study discussed here found that aspects of attributing the heavy Australian rainfall over 2010–2012 to specific causes were model dependent. While the set of studies we have reviewed provide robust evidence of anthropogenic influences on observed Australia temperature-related extremes, in combination with further results (e.g. Knutson et al., 2014; Christidis et al., 2013a), a comprehensive analysis of uncertainty has not yet been undertaken. We suggest a combination of approaches should be used in future attribution studies.
We also note that the approaches are present-day focused and consider only
anthropogenic influences on climate events already observed. As the world is
projected to further warm under growing greenhouse emissions, it is also
helpful to extend the extreme event attribution framework to understand
extremes in future, warmer climates. This was future-focused analysis
demonstrated in the time-slicing approach. For EEA studies quantifying the
human influence in temperature extremes using an attribution framework,
there is a general trend towards a greater degree of confidence in the
results of temperature-related studies. For some temperature extremes,
studies have identified events that are virtually impossible without
anthropogenic influences (FAR values of one)
(Imada
et al., 2018; Knutson et al., 2018; Lewis and Karoly, 2014b; Walsh et al.,
2018; Perkins-Kirkpatrick et al., 2019). More recent studies are
increasingly presenting risk ratio (RR) rather than FAR values, as a clearer
communication of attribution statements (see discussion in
National Academies of Sciences, Engineering, and Medicine, 2016) and in light
of the rapid saturation of FAR values for some event types
(Harrington et al., 2018). Reflecting on the results of the
studies presented here and these FAR issues, we suggest that further
applications of the methods we have outlined present risk ratio
(
While these studies have provided comprehensive evidence of the factors
(anthropogenic and natural) contributing to specific observed extremes,
several unexplored research directions should be noted. As the validity of
all event attribution statements depends on the representation of observed
climates and variability in models, a systematic analysis of uncertainty in
attribution statements in the Australian context would be useful for
providing recommendations around robust attribution statements. While a
larger sample of independent models could make up for
The temperature and precipitation-based studies highlighted here effectively demonstrate how anthropogenic climate change is already affecting extremes and is likely to do so further under the Paris Agreement warming levels. However, these studies only broach this key research area, focusing for the most part on large spatial- and temporal-scale extremes such as seasonal heat. Further studies exploring the contributing factors to complex climatological events (such as the emerging 2018 drought in New South Wales), compound events (such as recent heatwaves co-occurring with catastrophic fire weather in eastern Australia) or event impacts (such as on the Great Barrier Reef ecosystem or on human health) would be scientifically valuable as well as useful to a range of stakeholders and policy-makers. Future work attributing Australia's weather and climate extremes will target these complex events and impacts using CMIP6, in addition to other GCM-based datasets.
We have presented no new analyses here. AWAP, CMIP5 (
All the authors participated in the writing and editing of the manuscript. SCL developed Figs. 1 and 2, SEP Fig. 3 and ADK Fig. 4.
The authors declare that they have no conflict of interest.
This article is part of the special issue “Methods for detection and attribution of anthropogenic climate change”. It is not associated with a conference.
We acknowledge the support of the NCI facility in Australia and we acknowledge the World Climate Research Programme's Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modelling groups for producing and making available their model output. We also acknowledge the CESM Large Ensemble Community Project and supercomputing resources provided by NSF/CISL/Yellowstone. We thank the Bureau of Meteorology, the Bureau of Rural Sciences and CSIRO for providing the Australian Water Availability Project (AWAP) data.
Sophie C. Lewis is funded through the Australian Research Council (DE160100092). Sarah E. Perkins-Kirkpatrick is funded through the Australian Research Council (FT170100106). Andrew D. King is funded through the Australian Research Council (DE180100638).
This paper was edited by Francis Zwiers and reviewed by two anonymous referees.