the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Weather persistence on subseasonal to seasonal timescales: a methodological review
Alexandre Tuel
Olivia Martius
Persistence is an important concept in meteorology. It refers to surface weather or the atmospheric circulation either remaining in approximately the same state (quasistationarity) or repeatedly occupying the same state (recurrence) over some prolonged period of time. Persistence can be found at many different timescales; however, subseasonal to seasonal (S2S) timescales are especially relevant in terms of impacts and atmospheric predictability. For these reasons, S2S persistence has been attracting increasing attention from the scientific community. The dynamics responsible for persistence and their potential evolution under climate change are a notable focus of active research. However, one important challenge facing the community is how to define persistence from both a qualitative and quantitative perspective. Despite a general agreement on the concept, many different definitions and perspectives have been proposed over the years, among which it is not always easy to find one's way. The purpose of this review is to present and discuss existing concepts of weather persistence, associated methodologies and physical interpretations. In particular, we call attention to the fact that persistence can be defined as a global or as a local property of a system, with important implications in terms of methods and impacts. We also highlight the importance of timescale and similarity metric selection and illustrate some of the concepts using the example of summertime atmospheric circulation over western Europe.
Surface weather persistence at subseasonal to seasonal (S2S) timescales can have severe impacts on human and natural systems. Longlasting dry conditions, for instance, can lead to droughts and wildfires and can affect agriculture and energy production. Longlasting wet spells may cause severe flooding and crop loss. Persistent surface weather can result either from quasistationary, longlived atmospheric circulation conditions (quasistationarity) or from repeated, shorterlived circulation features (recurrence). Recurrence refers to the repeated occurrence of similar largescale circulation types or weather systems within some (S2S) time interval, usually with brief interruptions. Many recent highimpact weather and climate events were linked to persistent quasistationary or recurrent weather conditions. An example for recurrence are the western European floods in July 2021 that occurred at the end of an extreme wet spell in western Europe. The wet spell resulted from repeated atmospheric blocks and Rossby wave breaking episodes (Tuel et al., 2022b). Other examples of recurrence include the floods in the UK during winter 2013–2014 and in Queensland (Australia) in February–April 2022 both of which were caused by sequences of cyclones (Huntingford et al., 2014; Wikipedia, 2022; Floodlist, 2022). An example for quasistationarity is the catastrophic flooding in Pakistan in summer 2022 that resulted (in part) from longlasting and particularly heavy monsoon rains (Mallapaty, 2022). Intense heat waves and associated atmospheric circulations also tend to be persistent (often arising from a combination of recurrence and quasistationarity) (Lorenz et al., 2010), as in western Europe in 2003 (Black et al., 2004; GarcíaHerrera et al., 2010), western Russia in 2010 (Drouard and Woollings, 2018; Di Capua et al., 2021), or China (World Meteorological Organization, 2022) and India (Bloomberg, 2022) in 2022.
S2S weather persistence offers the potential for improved predictability at the S2S timescale (Franzke et al., 2011), which is highly relevant for risk preparedness and is attracting increased attention from the research community (e.g., Vitart et al., 2017; Meehl et al., 2021; Domeisen et al., 2022). However, a persistent state is not necessarily highly predictable, and persistent states with low predictability can cause large errors in subseasonal weather forecasts (Quandt et al., 2017).
Persistence is also an important aspect of climate model evaluation and climate projections. Whether global climate models are able to correctly simulate persistence is key to the robustness of longterm projections, especially of highimpact weather events – all the more so as climate projections suggest enhanced persistence (Li and Thompson, 2021; Hoffmann et al., 2021; Tuel and Martius, 2021a) in the future.
Characterizing weather persistence is therefore key to our understanding of the atmospheric circulation and its predictability and the associated hazards. Previous studies have focused on weather persistence from varied perspectives. Some assessed the persistence of specific weather systems or features, like atmospheric blocking (e.g., Liu, 1994), Rossby waves (Röthlisberger et al., 2019) or teleconnection patterns (Barnes and Hartmann, 2010). Others analyzed specific episodes of particularly persistent weather conditions (e.g., Black et al., 2004; Di Capua et al., 2021; Tuel et al., 2022b) or characterized the overall tendency of the atmospheric circulation and surface weather to exhibit persistence (e.g., MacDonald, 1992; Li and Thompson, 2021; Hoffmann et al., 2021). However, while previous studies generally agree on what persistence means conceptually, past work on this topic has involved many different definitions, often causing confusion and leading to different interpretations of persistence. Many case studies have also described observed situations as persistent based on subjective analyses rather than quantitative metrics. A further source of confusion is that “recurrence” is used in the literature to refer not only to successive occurrences of the same weather pattern at close intervals – what we will focus on in this review – but also to the states of the atmosphere with the highest probability of occurrence (Michelangeli et al., 1995).
It is difficult to give a unique definition of weather persistence. Besides, it may not even be desirable as different interpretations are possible and useful, depending on the system and timescale of analysis and on the motivations and goals of the study. Our goal here is therefore to review existing concepts of weather persistence, associated methodologies and physical interpretations. We present and structure a wide variety of approaches, definitions, techniques and metrics that have been used to analyze these concepts and that allow for answering one or several of the following questions.

Is there persistence in the data?

What are the persistent timescales in the data?

In which specific periods does persistence occur?

What are the persistent locations in the state space?
While persistence occurs at many different timescales, we specifically focus on the S2S timescale – often highly impact relevant and certainly important for predictability – but most methods and arguments in principle also apply to longer timescales. An overview of the methods is given in Fig. 1. We begin in Sect. 2 by introducing the two aspects of persistence: quasistationarity and recurrence. Section 3 then discusses several different perspectives on persistence to consider when choosing an analysis methodology. Finally, we present a detailed list of methods to detect or quantify persistence in Sects. 4 and 5. We illustrate most of the methods with examples taken from the literature or with our own analyses of summertime atmospheric circulation over Europe (details about the data we use are given in the data availability section). We keep the interpretation of the results to a minimum, the point being to illustrate the methods and not to analyze European summer circulation persistence in detail.
Persistence in a dynamical system (climate variable, atmospheric circulation field, etc.) arises from the repeated occurrence of the same value(s) or pattern(s) over a period of time. Successive occurrences can follow each other continuously – a situation we refer to as “quasistationary” – or in an interrupted sequence – in which case we speak of recurrence. Persistence therefore comes in two flavors, quasistationarity and recurrence, which we illustrate with the example of extreme warm conditions in a region north of the Black Sea in Fig. 2. Here, extreme warmth is defined as a daily mean temperature exceeding 1.5 standard deviations above its annual cycle. In summer 2010, the region experienced persistent (quasistationary) extreme warmth, with temperatures remaining continuously above the threshold for 41 d straight (Fig. 2a). By contrast, in spring 1979, extreme warm temperatures occurred frequently throughout May, but the three extreme warm episodes within 1 month were separated by at least 5 d with nearaverage temperatures (Fig. 2c). The corresponding evolution of the atmospheric circulation can be conceptually described as a slowmoving trajectory for quasistationary conditions (Fig. 2b) and a more rapidly evolving trajectory that revisits the same point repetitively (Fig. 2d).
The literature often uses persistence as a synonym for quasistationarity (e.g., Dole and Gordon, 1983; Barnston and Livezey, 1987; Franzke et al., 2011; Di Lorenzo and Mantua, 2016; Fereday, 2017; Liu et al., 2018; Du et al., 2019; Francis et al., 2020; Hoffmann et al., 2021; Li and Thompson, 2021), i.e., persistence is associated with longlived flow anomalies and little change in atmospheric circulation and surface weather. Recurrence has by contrast attracted less scientific attention. Most studies on recurrence have focused on extreme or impactful events (e.g., Mailier et al., 2006; Barton et al., 2016; Dacre and Pinto, 2020; Tuel and Martius, 2022a) and Rossby waves (e.g., Röthlisberger et al., 2019; Ali et al., 2021). Yet, from a physical perspective, quasistationarity and recurrence are intimately related. Recurring weather systems typically result from quasistationary favorable largescale conditions, like sea surface temperature (SST) anomalies or the location of extratropical jets (e.g., Tuel and Martius, 2022b). Temporal dependence, or memory, in a system can thus translate in practice into both quasistationarity and recurrence (Franzke, 2013). Additionally, from the impacts perspective, it makes sense to look at quasistationarity and recurrence together, since both can cause prolonged, impactful surface weather. Recurrent Rossby waves, for instance, modulate the persistence of surface temperature and precipitation anomalies (Röthlisberger et al., 2019; Ali et al., 2021) while longlived SST anomaly patterns like El Niño–Southern Oscillation (ENSO) can remotely trigger recurrent extreme weather (Gershunov and Barnett, 1998).
Note that recurrence, as we define it here, is sometimes referred to as “temporal” or “serial clustering” (Franzke, 2013), for instance in the case of recurrent cyclones (Mailier et al., 2006) or heavy precipitation events (Barton et al., 2016). Note also that recurrence is also frequently used in meteorology to refer to preferred patterns that repeatedly occur in a time series but not necessarily at close intervals (e.g., Vigaud et al., 2018; Kornhuber et al., 2019; Son et al., 2021) – in the case of atmospheric circulation patterns, one speaks of circulation or weather “regimes” (Michelangeli et al., 1995). This is different from our definition, in which recurrence specifically relates to the repeated occurrence of the same patterns over S2S timescales. Such patterns may be rare in the full dataset and would therefore not be considered “recurrent” in the regime approach but can be highly relevant from an impacts perspective. Hannachi et al. (2017) gave a comprehensive overview of the weather regime approach and associated methodologies. We also discuss the relevance of the weather regime perspective for quasistationarity analysis in Sect. 4.2.1.
Because quasistationarity and recurrence are two faces of the same coin, distinguishing one from the other may not be evident nor necessarily relevant.
First, the distinction often depends on the variable of interest. Recurrent weather systems can indeed result in quasistationary surface weather anomalies and vice versa. For example, the longlived heat waves of the 2010 and 2021 summers in western Russia (Fig. 2a) and the Baltic were linked to recurrent atmospheric blocks (Drouard and Woollings, 2018; Tuel et al., 2022b). Likewise, prolonged droughts or wet spells may result from recurrent Rossby wave activity (Röthlisberger et al., 2019; Ali et al., 2021, 2022). Quasistationary surface warm and humid conditions can also trigger recurrent thunderstorm activity (Mohr et al., 2020). Conversely, recurrent extreme precipitation events (e.g., Barton et al., 2022) or extratropical cyclones (e.g., Dacre and Pinto, 2020) are often linked to quasistationary jet states that last for much longer than the lifetime of individual weather systems.
Second, the longer the timescale of analysis, the less obvious the difference between quasistationarity and recurrence becomes. Impacts, for example, often depend on anomalies of surface temperature or precipitation averaged or accumulated over several weeks to months (e.g., droughts). Weekly or monthly values may thus sometimes be preferred to daily ones, in which case synopticscale quasistationarity and recurrence would both result in large weekly or monthly anomalies that would result from simple “persistence”.
Nevertheless, the distinction between quasistationarity and recurrence remains highly relevant for several reasons: from a methodological perspective (Sects. 4 and 5) but also for forecasting and process understanding at the synoptic timescale and for some stakeholders like insurers (for whom it matters whether impacts resulted from a single or multiple events). For process understanding, taking a weather systems perspective is often very relevant. In such a case, distinct weather systems (like cyclones) can in principle be separated from one another, and longlived single systems be distinguished from multiple shortlived systems occurring in close succession. It therefore matters whether persistence is driven by recurrence or quasistationarity. The distinction is also important to assess whether numerical models simulate persistence for the right physical reasons.
Before reviewing how persistence can be apprehended and quantified, let us begin with some basic notation and definitions. In the following, we denote the dynamical system under analysis with x(t)∈ℝ^{m}.
x(t) could, for example, be the weather evolution in a region, the time evolution of one variable at one location or the general circulation. (x(t))_{t} evolves within a state space 𝒮 that consists of the set of all possible system states ${\left\{{s}_{i}\right\}}_{i}$ that the system can occupy. A single system value and a single system state corresponds to each time t, x(t). The same system state can be attained at multiple different time steps. The series of successive system values constitutes the system trajectory. In practice, only a finite number of observations are available $\left(\mathit{x}\right({t}_{k}){)}_{\mathrm{1}\le k\le N}$ to characterize the persistence of a system.
Characterizing a system as “persistent” can mean different things. Consequently, it is important to always specify the perspective that is taken to avoid confusion. First, there are different flavors to persistence (global, state or episodic persistence; see Sect. 3.1). Second, persistence can be studied from a Lagrangian or an Eulerian perspective (Sect. 3.2). Finally, persistence is linked to a similarity metric (Sect. 3.3) and a notion of timescale (Sect. 3.4), for which different approaches are possible.
3.1 Global, state and episodic persistence
Persistence (whether quasistationarity or recurrence) is a broad concept that covers different kinds of behavior in dynamical systems. We might say, for instance, that temperature is more persistent than precipitation because temperature evolves, on average, over longer timescales than precipitation. In this sense, persistence quantifies the system's inertia. However, if we qualify last summer's weather as particularly persistent, we mean something different: namely, last summer's weather varied much less than what one could have reasonably expected in a normal summer. Here, persistence refers to some unusual behavior of the system over a particular period. Saying that zonal jets or atmospheric blocks are persistent again means something else, i.e., that these particular states of the circulation tend to be more longlived or recurrent than other states.
This leads us to make the distinction between three types of persistence: global, state and episodic persistence.

Global persistence. This characterizes the tendency to quasistationarity or recurrence across the whole trajectory of the dynamical system. We choose this term because persistence in this sense is a “global” property of the system, meaning that it is not restricted to any particular system state or time period. Global quasistationarity translates into the tendency for the system to change little at small timescales (successive values being close to each other). In mathematical terms, this translates into $<\left\frac{\mathrm{d}\mathit{x}\left(t\right)}{\mathrm{d}t}\right>\ll \frac{{\mathit{\sigma}}_{x}}{T}$, where $<\cdot >$ is a time average, σ_{x} is the standard deviation of the series x(t) and T its typical timescale of evolution. Global persistence can sometimes be impact relevant (e.g., trends in global persistence under climate change can be important for impacts). However, it does not make distinctions between the different system states. Consequently, global persistence is not suited to characterize persistent system states or specific time intervals with persistent system behavior and is generally not the bestsuited approach for risk assessment. Global persistence is, however, strongly related to intrinsic system predictability, since present values of the system largely constrain its future values. It can thus yield important information for numerical forecasting, including at the S2S timescale. Most global persistence methods focus on quasistationarity – like the autocorrelation coefficient (Sect. 4.1.1) or the Hurst exponent (Sect. 4.1.2) – but some exist for recurrence as well (Sect. 5.2.2). Many studies look at global persistence, and we may cite, for instance, MacDonald (1992); Pfleiderer and Coumou (2018); Pfleiderer et al. (2019) and Li and Thompson (2021), who analyze the persistence of temperature series, and Hoffmann et al. (2021), who consider 10 d atmospheric flow persistence.

State persistence. This relates to the persistent behavior of specific system states. State space persistence analysis consists of either identifying persistent system states – e.g., with the quasistationary (Sect. 4.2.2) or optimally persistent pattern (Sect. 4.2.3) methods – or characterizing the persistence of given states with methods like residence times (Sect. 4.2.5) for quasistationarity or Ripley's K for recurrence (Sect. 5.1.3). State persistence is highly relevant from the impact, forecasting and process understanding perspectives. Knowing which states are persistent and to what extent is useful (i) to make the link to surface impacts and (ii) to know which weather patterns or sequences may be more predictable than others. It also helps to determine the physical processes that support or are responsible for the persistence. The state approach to persistence has been used to shed light on quasistationary states of the North Atlantic circulation and their predictability (Faranda et al., 2017b), to characterize the quasistationarity of continentalscale weather patterns (Francis et al., 2018, 2020) or drought (Ford and Labosier, 2014), or to quantify recurrence in extratropical cyclones (Mailier et al., 2006) and extreme precipitation (Tuel and Martius, 2021a).

Episodic persistence. This is tied to specific time intervals (or “episodes”) during which the system exhibits quasistationary or recurrent behavior. It is in that sense a purely local property that characterizes the anomalous behavior of the system over a limited time period. Importantly, episodic persistence can occur simply by chance in any dynamical system, even in systems that exhibit no global persistence. Episodic persistence is therefore not necessarily relevant for predictability and process understanding. Persistent periods can still always be analyzed to look for potential drivers that discriminate pure “statistical flukes” from possibly predictable events. Episodic persistence is, however, well suited for impacts analysis because it can make a direct link between periods of persistent weather and impacts. Relevant methods include running window techniques (Sect. 4.3) for quasistationarity and window counts (Sect. 5.1.1) for recurrence. Hoffmann et al. (2021) investigated for example quasistationarity in 10 d sequences of atmospheric circulation, while Bevacqua et al. (2020) and Kopp et al. (2021) looked at subseasonal periods with recurrent cyclones and extreme precipitation events.
While it helps to capture the various interpretations of persistence, this classification is not perfect, and there is some overlap between categories. In practice, the state or episodic perspectives can also be used to quantify global persistence (by averaging persistence metrics across systems states or time intervals) and global persistence metrics can be computed on subsets of the data to identify persistent periods. Some methods, like recurrence plots (Sect. 5.2.2), can even deal with all three types of persistence.
The three types of persistence are also not independent from one another. Global persistence, for instance, can emerge from repeated occurrences in one or a handful of system states only, while the rest of the trajectory, if analyzed separately, may not be qualified as persistent. There are also strong relationships between state and episodic persistence. The most common system states to occur during persistent periods are indeed likely to be persistent states. Correspondingly, persistent states, when they occur, are likely to be associated with persistent periods. Persistent states can thus be uncovered from the knowledge of persistent periods, for instance with pattern recognition or clustering algorithms applied to system values during persistent periods, or simply by averaging system values during persistent periods (e.g., Faranda et al., 2017b; Hoffmann et al., 2021).
However, state persistence only characterizes the average behavior of system states – it does not imply that all occurrences of a persistent state will necessarily be persistent. Consequently, there is no onetoone relationship between persistent states and persistent periods. Some system states can behave in a persistent way under certain conditions but not under others. The lifetime and travel speed of atmospheric blocks, for instance, is affected by land–atmosphere feedbacks or upstream latent heating (Steinfeld et al., 2020). Similarly, extratropical cyclones may occur in sequences but also as single events (Dacre and Pinto, 2020). Additionally, persistent periods may be associated with a variety of system states. While we expect persistent states to be the most frequent ones during persistent periods, some states which on average are not persistent can still at times behave persistently by pure chance.
However, this classification is useful to illustrate the different methodological ways that weather persistence can be tackled, and we rely on it to structure the description of methods in Sects. 4 and 5.
3.2 Lagrangian and Eulerian perspectives
Weather persistence is most often analyzed from an Eulerian perspective, i.e., persistence of the same weather pattern or conditions at a fixed location in space. By contrast, in the Lagrangian perspective the focus is on the persistence of a given weather pattern in time. In the Eulerian framework, the system x(t) represents a time series over a domain fixed with time, whereas in the Lagrangian one x(t) follows individual weather patterns or the background atmospheric flow. The Eulerian stationarity of a quantity ψ(x,t) translates into $\partial \mathit{\psi}/\partial t=\mathrm{0}$, while the Lagrangian stationarity implies $\partial \mathit{\psi}/\partial t+\mathit{u}\cdot \mathrm{\nabla}\mathit{\psi}=\mathrm{0}$, where u is the background flow. Similarly, quasistationarity translates into $\partial \mathit{\psi}/\partial t\approx \mathrm{0}$ in the Eulerian framework and $\partial \mathit{\psi}/\partial t+\mathit{u}\cdot \mathrm{\nabla}\mathit{\psi}\approx \mathrm{0}$ in the Lagrangian one. For example, temperature anomalies can be tracked over time (e.g., Kornhuber and TamarinBrodsky, 2021) or analyzed at a fixed location (e.g., Pfleiderer and Coumou, 2018; Li and Thompson, 2021). Similarly, weather systems such as blocking, cyclones or vortices (Bray and Cavallo, 2022) can be analyzed at a fixed location or following the weather systems. For example, Økland and Lejenäs (1987) contrast the persistence of blocking at fixed longitudes against the persistence of individual blocking episodes. Another example is Kossin (2018), who characterize the persistence in tropical cyclones by their translation speed.
Both the Eulerian and the Lagrangian perspectives are relevant for impact and risk assessment. The former links persistence to impacts at a given location, and the latter highlights impacts along the trajectory of a weather pattern (system) during its lifetime. Indeed, the same weather system can produce hazardous weather over large areas, putting strain on the resources of insurance companies or governments. The two perspectives can be brought together by considering the translation speed and lifetime of the tracked weather systems. Systems with long lifetimes and low translation speed lead to both Eulerian and Lagrangian persistence. By contrast, longlived systems that travel quickly are persistent from a Lagrangian perspective only. Likewise, slowmoving but shortlived systems are not Lagrangian persistent, but Eulerian persistence can still be detected if several such systems occur over the same area in close succession.
3.3 Quantifying similarity
Assessing persistence typically requires quantifying the selfsimilarity of system values x(t) with a metric (e.g., Wharton et al., 2008; Zerzucha and Walczak, 2012; Ali et al., 2020; Ontañón, 2020). By selfsimilarity, we refer to the tendency for successive values of x(t) to be similar to each other according to some metric. This is not to be confused with the concept of geometric selfsimilarity in fractal geometry. Metric selection is an important step that should be done with care because it conditions how persistence is quantified and some metrics are not suited to certain kinds of data (e.g., heavily skewed). There are two main classes of similarity metrics.

Categorical metrics. With these metrics, system values are either similar (if they belong to the same category) or not (if they do not).

Continuous metrics. These metrics measure the degree of similarity between two system values in a continuous way.
Categorical metrics focus on specific features of the system, such as the presence of a given weather pattern or the occurrence of a specific event. They require the set of all x(t) values to be classified into distinct categories (usually from two to a few dozen). Categories can represent predefined system states of interest (e.g., a warm or cold anomaly, a given phase of a teleconnection mode, or the occurrence of a specific weather pattern, like a block) (e.g., Pinto et al., 2014; Drouard and Woollings, 2018; Pfleiderer and Coumou, 2018; Ali et al., 2021; Kopp et al., 2021) but can also be obtained objectively with dimension reduction or pattern recognition methods. Examples include principal component analysis (or empirical orthogonal functions, EOFs) (Fereday, 2017, e.g.,), optimally interpolated patterns (OPPs) (Hannachi, 2008), selforganizing maps (SOMs) (e.g., Francis et al., 2018; Weiland et al., 2021; Rousi et al., 2022b) or clustering algorithms (probabilistic, hierarchical or nonhierarchical) (e.g., Demuzere et al., 2011; Hannachi et al., 2012; Grams et al., 2017). When x(t) represents 2D circulation data (like sea level pressure or geopotential height), the ensemble of system states is often referred to as “weather regimes” (Michelangeli et al., 1995; Grams et al., 2017; Francis et al., 2018) (see Sect. 4.2.1). In EOF analysis, the distance metric is the L^{2} norm (Euclidean distance), but most clustering methods can work with any distance metric. The challenge with dimension reduction methods is choosing the number of categories to retain (EOFs, clusters, SOM nodes, etc.). A high number can capture rare states of the system, but this comes at the cost of making persistence more difficult to assess (since sequences when the system remains in the exact same category will become less frequent).
Continuous metrics measure the degree of similarity between any pair of system values in a continuous way. Common examples of continuous metrics include the Euclidean distance (e.g., Faranda et al., 2017b), pattern correlation(e.g., Mo and Ghil, 1987), or more complex similarity indices like the Teweles–Wobus score (e.g., Horton et al., 2017; Blanchet et al., 2018) or the image structural similarity index (SSIM) (e.g., Hoffmann et al., 2021).
In comparison to categorical metrics, continuous metrics offer the advantage that they do not require specifying features or events of interest beforehand. They are also more flexible insofar as similarity can be quantified with respect to any system value as reference and not just representative values for each category. Since they do not require simplifying the state space, continuous metrics may also be able to pick up rare persistent patterns that are missed by dimension reduction methods that focus on the most common patterns in a series. These advantages come at a cost: working with continuous metrics, especially complex ones, can be more computationally intensive and sometimes more difficult to interpret physically.
3.4 Persistence timescales
Persistence is linked to a notion of timescale during which the system continuously remains in the same state (for quasistationarity) or occupies that same state repeatedly (for recurrence). There are three common ways to approach persistence timescales.
The first option is to choose a single, fixed timescale for analysis. This choice can be guided by impact and forecasting considerations, by physical knowledge of the underlying system, or by observations of persistent events (e.g., Huntingford et al., 2014; Lawrence et al., 2020; Overland and Wang, 2021; Rakovec et al., 2022; Rousi et al., 2022a). This is the most common approach to analyze episodic persistence, but it is also applicable to global persistence. Quasistationarity can be quantified by the average similarity between n successive states (for continuous similarity metrics see, e.g., Kolstad et al., 2017; Hoffmann et al., 2021) or by the variety of system states during an nstep window (for categorical metrics see, e.g., Fereday, 2017; Richardson et al., 2019). Quasistationarity can also be inferred from extreme anomalies of circulation, temperature or precipitation during nstep windows. For instance, Gálfi et al. (2019) and Tuel and Martius (2023) identify persistent warm and cold spells by averaging temperature anomalies over 1–3 weeks.
Recurrence can similarly be assessed by calculating the number of times that a particular system state or event occurs during nstep windows (e.g., Mailier et al., 2006; Pinto et al., 2014; Kopp et al., 2021; Tuel and Martius, 2022a). A fixed timescale can also be used as a threshold to separate quasistationary from nonstationary events, by requiring quasistationary events to last at least n steps. For instance, Francis et al. (2018) and Francis et al. (2020) define persistent periods by requiring the circulation pattern in a region to remain in the same state for at least 4 consecutive days. Mann et al. (2018) similarly define persistent resonant wave events as those lasting at least 10 d. Dole and Gordon (1983) also take this approach for the persistence of pointwise geopotential anomalies. Note that this fixed timescale can also consist of a single time step. In this case, persistence characterizes how much past system values determine future ones. Li and Thompson (2021), for instance, characterize persistence with the lag1 autocorrelation coefficient, Röthlisberger and Martius (2019) look at 1 d transition probabilities between different system states, and Vautard (1990); Michelangeli et al. (1995) and Hannachi et al. (2017) calculate time derivatives of geopotential fields to identify quasistationary states.
The second option is to select an analysis method that explores a range of timescales and pinpoints the relevant persistence timescales, sometimes accompanied by some notion of statistical significance. This approach only works for global persistence. Autocorrelation analysis, for instance, detects the timescales at which the system exhibits significant lagged memory (Sect. 4.1.1). Spectral analysis can likewise highlight important timescales of variability that can be linked to quasistationarity (Sect. 4.1.2). For recurrence, methods like Ripley's K function indicate the timescales at which recurrence is statistically significant (Sect. 5.1.3).
Finally, the third option is to characterize persistence not by a single timescale but by a distribution of timescales. To assess quasistationarity, one can typically work with the distribution of persistent event durations. Persistent events are periods during which the system satisfies a persistence criterion: for quasistationarity, successive system values must be similar, and for recurrence the same event must occur multiple times, with each occurrence being separated by at most n time steps from the previous one. For recurrence, it is also possible to consider the distribution of interevent times (e.g., Altmann and Kantz, 2005). The distribution of event lengths can be further summarized by considering the average or maximum event length (Faranda et al., 2017b; Rousi et al., 2022b) or by modeling it with an exponential or power law distribution (Sect. 4.2.5). This approach has been applied to numerous cases: heat waves (Lorenz et al., 2010), droughts (Meng et al., 2017; Moon et al., 2018), wet spells (Ali et al., 2021), atmospheric blocks (Liu, 1994), circulation patterns (Huguenin et al., 2020), and midlatitude cyclone clustering (Bevacqua et al., 2020).
The diversity of perspectives on persistence translates into a wide range of methods, of which we give an overview in the following two sections dedicated, respectively, to quasistationarity and recurrence. Following the distinctions introduced in Sect. 3, we separate methods that quantify global, state, and episodic persistence (though some methods can be used for more than one). We also specify whether methods can only be used with a single timescale or whether they quantify persistence across timescales. An overview of methods is shown in Table 1.
4.1 Global quasistationarity
We begin with several methods that quantify global quasistationarity in onedimensional time series. They characterize quasistationarity in the series as a whole but are generally unable to identify quasistationary states or periods. However, they often allow the user to characterize the timescales of variability and persistence in the data and are hence relevant for system predictability and process understanding.
4.1.1 Autocorrelation
Autocorrelation is a frequently used measure of quasistationarity in weather and climate science. If X_{t} is a continuous, onedimensional process of mean μ and variance σ^{2}, its Pearson autocorrelation coefficient at lag k is defined as follows:
where 𝔼 denotes the expectancy with respect to the distribution of X_{t}. Potential cycles and longterm trends should be removed from the data before analysis (Weiss and Weiss, 1999). In Eq. (1), X_{t} can also be replaced by its rank, which yields the alternative Spearman autocorrelation that is more robust to nonlinear behavior in the data. The definition of Eq. (1) can also be extended to thirdorder statistics to capture interactions coming from nonlinear correlation, yielding the bicorrelation $\mathit{\rho}({k}_{\mathrm{1}},{k}_{\mathrm{2}})=\frac{\mathbb{E}\left[\left({X}_{t}\mathit{\mu}\right)\left({X}_{t+{k}_{\mathrm{1}}}\mathit{\mu}\right)\left({X}_{t+{k}_{\mathrm{2}}}\mathit{\mu}\right)\right]}{{\mathit{\sigma}}^{\mathrm{3}}}$. (Pires and Hannachi, 2021).
Several summary metrics for autocorrelation exist, like the autocorrelation timescale:
the decorrelation time (Hannachi, 2021)
where N is the number of time steps (if ρ is integrable); or the characteristic time ${T}_{\mathit{\rho}}^{\mathrm{c}}$ (Trenberth, 1984) (Fig. 3a):
The absolute value on ρ(k) in Eqs. (3) and (4) is not required but can help avoid underestimated timescales in case ρ exhibits oscillations (Franzke et al., 2005).
Both ${T}_{\mathit{\rho}}^{\mathrm{d}}$ and ${T}_{\mathit{\rho}}^{\mathrm{c}}$ roughly approximate the average time between successive independent values. The lag1 autocorrelation ρ(1) can also be used as a quasistationarity metric (Li and Thompson, 2021), which is equivalent to fitting a linear regression model between successive system values: $\mathit{x}(t+\mathrm{1})=\mathit{\beta}\mathit{x}\left(t\right)+\mathit{\u03f5}\left(t\right)$, where ϵ(t) is a whitenoise process. If x(t) is normalized, the regression slope β is equal to ρ(1). Note that more complex models are possible: linear models for continuous series can be extended by the large class of autoregressivemoving average (ARMA) models.
Autocorrelation has many advantages: it has a simple definition, the ease of interpretation (ρ(k) of which is related to the linear regression coefficient of X_{t+k} against X_{t}); it has high flexibility; it is already implemented in common programming languages; it requires no subjective threshold choice; and its results can be easily reproduced. By varying the lag k, it can measure quasistationarity at all timescales. It also comes with a notion of statistical significance: given a confidence level, it is possible to say whether the obtained autocorrelation is significant (indicating a link at lag k) or not, pointing to the relevant quasistationarity timescales in the series. On the downside, autocorrelation only works for onedimensional data and requires a large number of values as input. Therefore, it is best suited to measuring quasistationarity globally. One can compute autocorrelation on a subset of the data only (summer values or a specific time interval, for instance) if the number of data points is large enough, in which case autocorrelation may be used to characterize quasistationarity locally in time. Additionally, autocorrelation only highlights linear relationships (though bicorrelation can help capture nonlinearity). Finally, autocorrelation measures the strength of the connection between lagged system values but not how far apart they might be in the state space.
Example applications include Horel (1985a) and Barnston and Livezey (1987), who calculated lag autocorrelation on principal component time series of monthly Northern Hemisphere geopotential fields. MacDonald (1992) used autocorrelation to detect quasistationarity in monthly temperature series, as did Weiss and Weiss (1999) to assess quasistationarity in ENSO. Weatherhead et al. (2010) characterized daily temperature quasistationarity in the Arctic based on lag1 autocorrelation. Degenhardt and Ólafsson (2019) and Kolstad et al. (2015) calculated lag1 autocorrelation to highlight intraseasonal quasistationarity of monthly mean temperatures in Iceland and Europe, respectively, and Li and Thompson (2021) applied autocorrelation to daily temperature series and found it was strongly related to the average length of warm and cold episodes across the world. Kolstad et al. (2017) even used autocorrelation analysis in a causal discovery framework by regressing temperature values against previous ones and including potential covariates.
We show the characteristic time for daily rescaled 500 hPa geopotential height over Europe during summer in Fig. 3a. While the results say nothing about the role of individual weather systems, the larger values over the British Isles and western Russia seem consistent with more frequent persistent atmospheric blocks over these regions.
4.1.2 Asymptotic methods to characterize variability across timescales
Persistence in a system is associated with memory effects that lead to variability being concentrated at long rather than short timescales. How variability in the series is distributed across timescales is therefore an important indicator of global persistence and can point to relevant persistence timescales. Specifically, persistence often translates into scaling laws: in the time series variability as a function of frequency or timescale (for quasistationarity) (Bunde et al., 2002) but also in the distribution of interevent times (for recurrence, see Sect. 5.1.4). The slope of these laws indicates the degree of persistence. Time series also often exhibit different scaling laws over different frequency intervals, highlighting how persistence may differ across timescales in the data.
A first scaling law can be obtained by considering how autocorrelation decreases as a function of time lag. The more quasistationary a series, the less rapidly its autocorrelation should decrease. Conceptually, time series can be broadly divided between shortrange and longrange dependence series. For a shortrange dependent series (like autoregressive models), the autocorrelation decreases rapidly with the time lag, eventually reaching 0 after a certain lag or decaying exponentially.
as k→∞ and with $\mathrm{0}<\left\mathit{\alpha}\right<\mathrm{1}$ (with α<0 corresponding to antipersistence). By contrast, in a longrange dependent series, the autocorrelation decays following a power law:
where d is called the dependence parameter ($\mathrm{0}<d<\frac{\mathrm{1}}{\mathrm{2}}$) (Beran, 2017; Franzke, 2013). A lower d is associated with a slower decay of the autocorrelation function and hence with more quasistationarity systems. White noise has d=0. Witt and Malamud (2013) and Franzke (2013) discuss several ways how d can be estimated.
Like autocorrelation, α and d (the dependence parameter) measure global quasistationarity in continuous time series because they are based on asymptotic relationships. However, they cannot highlight specific quasistationarity timescales in a series as autocorrelation does. Applications to atmospheric time series show that temperature exhibits longterm dependence (e.g., Yuan et al., 2010; KoscielnyBunde et al., 1998; Eichner et al., 2003), while precipitation exhibits either short or longterm dependence (Potter, 1979; Hannachi, 2014; Yang and Fu, 2019). Franzke (2013) computed d and found longrange dependence in North Atlantic winds (Fig. 3b).
The Hurst coefficient (or exponent) H (Hurst, 1951) is a common measure of memory in a time series which is obtained from a second scaling law. Specifically, H characterizes how a time series (X_{t})_{t} fluctuates relative to its mean. Noting
the cumulative fluctuations in (X_{t})_{t} relative to its mean, calculated over intervals of size n, Hurst argued, empirically from geophysical time series, that these fluctuations could exhibit scaleinvariant properties. In other words, cumulative fluctuations over two different timescales n and m are related through
where stands for equality in distribution. H ranges from 0 to 1. $H<\mathrm{1}/\mathrm{2}$ indicates antipersistent behavior, such that successive increments of X_{t} relative to its mean (X_{t}−𝔼[X_{t}]) tend to be negatively correlated, and the time series fluctuates substantially at short timescales. By contrast, $H>\mathrm{1}/\mathrm{2}$ points to persistent behavior, in which successive increments are positively correlated, and variability is concentrated at long timescales. $H=\mathrm{1}/\mathrm{2}$ corresponds to white noise (no temporal correlation). Thus, the higher H is, the smoother the time series. H is also theoretically related to the dependency parameter (as $H=\mathrm{2}d+\mathrm{1}$) and the power spectrum exponent (see below) (Franzke et al., 2020).
H can be estimated in many ways (see, e.g., Koutsoyiannis, 2003; De las Nieves López García and Requena, 2019; Franzke et al., 2020). In the literature, H has mainly been used to characterize longterm dependence or (quasi)stationarity (Mandelbrot and Wallis, 1969), for instance in series of monthly or annual mean temperature (MacDonald, 1992; Kumar et al., 2013), precipitation (Bunde et al., 2013), or drought indices (Tatli, 2015). However, some studies also calculated it for daily time series (e.g., Rehman and Siddiqi, 2009; Velásquez Valle et al., 2013).
Like the Hurst exponent, spectral analysis also characterizes how a time series' variability is distributed across timescales. The power spectrum is commonly defined as the Fourier transform of the autocorrelation function ρ:
where f is the frequency and k the time lag. In a pure whitenoise series, the variability is distributed equally across frequencies. Hence, the power spectral density is a constant. If temporal dependence is present, however, the power spectrum typically exhibits a power law decrease with frequency f:
where β, called the power spectrum exponent, indicates the degree of quasistationarity in the time series ($\mathrm{0}<\mathit{\beta}<\mathrm{1}$). For statistically stationary series, one can also show that $\mathit{\beta}=\mathrm{2}H\mathrm{1}$, where H is the Hurst coefficient (Parzen, 1986).
The larger β is, the more variance is concentrated at low frequencies. This implies a memory effect at low frequencies that relates to quasistationarity in the time series. Such scaling behavior is very common in climatological series, and the spectrum is often divided into distinct scaling regimes corresponding to specific frequency intervals. Fraedrich and Larnder (1993) discuss the example of precipitation in continental Europe and relate the different regimes to specific physical processes and timescales (from individual storms at high frequencies to climate fluctuations at low frequencies) (Fig. 4). Yang and Fu (2019) obtain similar results using hourly and daily precipitation data for the United States. Pelletier and Turcotte (1997) also used power spectra to quantify quasistationarity in various monthly climatic series, as did Ault et al. (2014) for drought quasistationarity and MacDonald (1992) for temperature quasistationarity. Telesca et al. (2016) analyzed quasistationary regimes in 10 min wind series across Switzerland.
Note that, as with the autocorrelation, the definition of the power spectrum can be extended to thirdorder statistics to yield the bispectrum (the Fourier transform of the bicorrelation function) (Pires and Hannachi, 2021). While the power spectrum characterizes how the signal's variance is distributed across timescales, the bispectrum provides the contribution of each pair of frequency to the signal's skewness. Pires and Hannachi (2021), for instance, calculated the bispectrum of the 3monthly ENSO series to better understand its predictability.
4.2 State quasistationarity
We now turn to methods that focus on the quasistationarity of specific system states. Unlike global methods, which take onedimensional time series as input, several of the following methods are directly applicable to multidimensional data. We begin with methods that identify quasistationary states from the system trajectory (Sects. 4.2.1–4.2.4), before discussing methods that quantify the average quasistationarity of a system state (Sects. 4.2.5–4.2.6).
4.2.1 Weather regimes
We begin with an identification method for quasistationary states based on the concept of weather regimes (Michelangeli et al., 1995). This concept emerges from the realization that the extratropical atmospheric circulation evolves mainly as a succession of a handful of largescale circulation patterns (Hannachi et al., 2017). These preferred subseasonal flow patterns, or weather regimes, tend to be quasistationary over timescales of a few days to a few weeks and are therefore strongly related to weather persistence. They account for much of the lowfrequency atmospheric variability at intraseasonal timescales (Pandolfo, 1993; Hannachi et al., 2017). The existence of weather regimes in the midlatitudes has long been recognized, such as the concept of Grosswetterlagen (Baur, 1951) or atmospheric blocking (e.g., Namias, 1964). They have attracted considerable attention, in particular because of the potential longrange predictability they offer (e.g., Ghil and Robertson, 2002; Büeler et al., 2021) and their link to surface impacts (e.g., Grams et al., 2017).
There are many ways to calculate weather regimes (see Huth et al., 2008, and Hannachi et al., 2017, for detailed overviews). The most common methods include pattern recognition and dimensionality reduction techniques, applied to a proxy variable for the atmospheric circulation like 500 hPa geopotential height or sea level pressure. Examples include orthogonal pattern decomposition techniques like principal component analysis (in the temporal or spectral domain) (Fereday, 2017; Grams et al., 2017) or optimally interpolated patterns (OPPs) (Hannachi, 2008), selforganizing maps (Huth et al., 2008; Francis et al., 2020; Weiland et al., 2021), archetypal analysis (Hannachi and Trendafilov, 2017; Chapman et al., 2022), and clustering algorithms (Fig. 5). The latter can be probabilistic (e.g., Gaussian mixture models; Woollings et al., 2010), hierarchical (e.g., Ward clustering; Hannachi et al., 2012) and nonhierarchical (e.g., k means or partitioning around medoids; Grams et al., 2017). Various statistics (gap statistic, silhouette coefficient, etc.) can help objectively select an optimal number of clusters, many of which are available from the R
package clusterCrit
(Desgraupes, 2018). Note that, in practice, the input data should be normalized to remove longterm and seasonal trends to focus exclusively on intraseasonal variability (Grams et al., 2017).
Regimes can also be identified as local maxima in the (multidimensional) probability distribution function (PDF) of the target field, obtained empirically through e.g., kernel smoothing (Kimoto and Ghil, 1993; Woollings et al., 2010), or with more complex tools of network theory (Mukhin et al., 2022) and topology (Strommen et al., 2022). Finally, Franzke et al. (2011) identify regimes in the North Atlantic jet position with a hidden Markov model (HMM). HMMs are a powerful tool that brings together Markov models and Gaussian mixture models. Given N unknown (hidden) states, the HMM models the distribution of the observed series conditionally on each hidden state, with the sequence of hidden states following a firstorder Markov process (Franzke et al., 2008). The transition matrix, hidden states and conditional distributions can be estimated simultaneously.
The main drawback of the weather regime approach is that it is biased toward preferred or frequent states. Quasistationary but rare flow patterns that fall outside the range of the major regimes may therefore be missed. Additionally, certain patterns may fall in between regimes and are not classified or are misclassified. Despite this, weather regimes are very useful because they transform complex multidimensional systems into categorized, onedimensional series (according to which regime the system is closest to at each time step).
4.2.2 Quasistationary states
One major disadvantage of weather regimes is that they strongly simplify the state space. They also focus on the patterns that account for most of the variability in the data, regardless of their persistence. They may thus miss rare yet impactrelevant persistent patterns. Several other techniques exist to directly extract quasistationary patterns from the system trajectory.
Quasistationary circulation patterns can be seen as mathematically quasistationary solutions of the atmosphere's equations of evolution (Mo and Ghil, 1987). Such solutions are characterized by average system time derivatives close to zero, meaning that the system tends to remain in their vicinity longer than elsewhere in the state space – hence their link to quasistationary states. Quasistationary states can therefore be directly identified from the system's dynamics by looking for states whose time derivative is close to zero. If we know the system's exact evolution equations (in the case of simplified models, for example), strictly stationary solutions can be directly computed (e.g., Charney and DeVore, 1979; Legras and Ghil, 1985; Mo and Ghil, 1987). In practice, however, this is rarely the case, and quasistationary (or “metastable”) states are defined in a statistical sense only as those for which system tendencies (i.e., time derivatives) approach zero (Vautard, 1990):
𝒯(x^{*}) is the composite tendency at x^{*}, defined as the average (or areaaverage for multidimensional fields) of instantaneous tendencies at all occurrences of x^{*} – in practice at all times when the system is in the neighborhood of x^{*}:
where $\mathrm{d}(\cdot ,\cdot )$ is a similarity metric and d_{0} some small threshold. It is also possible to weigh the terms in Eq. (12) according to $\mathrm{d}({\mathit{x}}^{*},\mathit{x}(t\left)\right)$ (Vautard, 1990; Michelangeli et al., 1995). When dealing with complex multidimensional systems, it is simpler to calculate tendencies on the time series of the system's leading principal components. Equation (12) can then more easily be solved by minimizing $\left\mathcal{T}\right(\mathit{x}{)}^{\mathrm{2}}$ in x, with distinct solutions corresponding to different quasistationary states. In practice, instantaneous tendencies at any x^{*} can exhibit a large variance, and time series are commonly lowpass filtered to remove the shortterm noise that would complicate solving for Eq. (12). Additionally, $\left\mathcal{T}\right(\mathit{x}{)}^{\mathrm{2}}$ can be difficult to minimize as it is piecewise constant (due to the finite number of observed x values). Only approximate solutions are achievable. It may thus be difficult to know precisely how many zeros of $\left\mathcal{T}\right(\mathit{x}{)}^{\mathrm{2}}$ exist and whether two approximate solutions correspond to the same minimum. Vautard (1990) presents a method to select relevant solutions.
It is important to note that even if ${T}_{i}\left({\mathit{x}}^{*}\right)\approx \mathrm{0}$, instantaneous tendencies $\frac{\mathrm{d}\mathit{x}}{\mathrm{d}t}{}_{\mathit{x}\left(t\right)}$ in Eq. (12) may be far from zero (Vautard, 1990). Indeed, the time derivative of x(t) usually depends not only on x but also on other timedependent variables that may evolve independently from x.
Haines and Hannachi (1995) and Hannachi (1997) estimated quasistationary states over the North Pacific in the output from a global climate forced by perpetual January conditions, by projecting simplified dynamics (e.g., quasigeostrophy) onto the leading modes of variability in the global climate model simulation. Vautard (1990) found four main quasistationary patterns in wintertime daily 700 hPa geopotential height fields over the Atlantic and analyzed their quasistationarity and onset and break characteristics. Michelangeli et al. (1995) also looked at 700 hPa geopotential fields and compared quasistationary states with the leading EOFs over the Atlantic and Pacific oceans during winter. Mo and Ghil (1987) did a similar comparison for the Southern Hemisphere circulation during the austral winter.
Figure 6 shows five quasistationary Z500 anomaly patterns over Europe in summer. We obtained them from the 10 leading EOFs of the 5 d averaged Z500 fields, following Vautard (1990). The patterns were obtained by clustering the resulting approximate solutions with the highest number of clusters for which the resulting patterns were subjectively different. The solutions include blockingtype patterns over Scandinavia and western Russia (Fig. 6a, b) and zonal patterns (Fig. 6c, d) that are similar to the results of Figs. 5 and 11 and to the canonical patterns of variability over the European–Atlantic sector (Grams et al., 2017).
4.2.3 Optimally persistent patterns
The technique of optimally persistent patterns (OPPs) was introduced by DelSole (2001). OPP analysis is conceptually similar to principal component analysis; however, instead of patterns that maximize the variance in the observed series, OPPs are defined as the patterns whose time components (i.e., the projection of the observed series onto the OPP) are the most persistent. Persistence is here measured by the decorrelation time (Eq. 3) or, alternatively, by the squared decorrelation time ${T}_{\mathrm{2}}=\mathrm{1}+\mathrm{2}{\sum}_{k=\mathrm{1}}^{N}\mathit{\rho}(k{)}^{\mathrm{2}}$ DelSole (2001) (with ρ(k) the autocorrelation function and N the number of observations). If the dynamical system x(t) is mdimensional, then for an mdimensional pattern u, the time component of u is defined as y(t)=u^{T}x(t). If ρ_{y}(k) is its associated autocorrelation function, finding OPPs then consists of maximizing $\mathrm{1}+\mathrm{2}{\sum}_{k=\mathrm{1}}^{N}{\mathit{\rho}}_{y}\left(k\right)$ or $\mathrm{1}+\mathrm{2}{\sum}_{k=\mathrm{1}}^{N}{\mathit{\rho}}_{y}^{\mathrm{2}}\left(k\right)$. In the first case (maximizing the decorrelation time), the optimization reduces to an eigenvalue problem, details of which can be found in Hannachi (2021). The leading eigenvector maximizes the decorrelation time. OPPs can then be obtained by projecting the observed data x(t) onto the time series associated with each eigenvector. The second case (maximizing the squared decorrelation time) leads to a more complex nonlinear optimization problem that can be solved iteratively (see details in DelSole, 2001). Note that, as in Sect. 4.2.2, the system should first be embedded into a lowerdimensional space, e.g., by projecting it onto the set of leading EOFs. Figure 7 shows the leading OPP obtained by DelSole (2001) by maximizing the squared decorrelation time for daily Northern Hemisphere Z500 fields.
OPPs are especially relevant for forecasting since they correspond to the patterns with the most lowfrequency variability. However, like EOFs and quasistationary states, OPPs are mathematical objects that do not necessarily correspond to “real” patterns ever attained by the system.
Other methods have been proposed to identify persistent patterns in a spatiotemporal field. One consists of minimizing the onestepahead forecast error, where the forecast is obtained by projecting the observed field x(t) onto a set of temporal patterns. These patterns, named “predictive oscillation patterns” by Kooperberg and O'sullivan (1996), can thus be considered the “most predictable” (i.e., persistent) patterns of a spatiotemporal field. Another possibility is to minimize the interpolation error variance between x(t) and its estimate $\widehat{\mathit{x}}\left(t\right)$ obtained from all $\left(\mathit{x}\right({t}^{\prime}){)}_{{t}^{\prime}\ne t}$ values. In this case one speaks of optimally interpolated patterns (Hannachi, 2008). Further details for both methods can be found in Hannachi (2021).
4.2.4 Extreme values and dynamical systems theory
The quasistationary state method only focuses on the most quasistationary states of a system, and OPPs characterize quasistationarity for a handful of possible patterns only. By contrast, dynamical systems theory provides a convenient framework to describe the quasistationarity of any system state (Lucarini et al., 2016). In this framework, quasistationarity is defined for any point x_{0} of the state space as the inverse of the average residence time of trajectories around x_{0}. The residence time is calculated based on the distance between successive system values, with two successive values deemed similar if their distance is below some small threshold. For any state x_{0}, the probability that x(t) will remain in a close neighborhood around x_{0} (a ball of radius ϵ) can be approximated by Faranda et al. (2017a) and Faranda et al. (2017b)
where $\mathrm{d}(\cdot ,\cdot )$ is a distance function (Euclidean distance in Faranda et al., 2017a, b) and θ is called the extremal index. ϵ is usually chosen to be the second percentile of the d(x(t),x_{0}) values. θ can be estimated in various ways (Hamidieh et al., 2009; Holešovský and Fusek, 2022), with two common ones being the intervals estimator of Ferro and Segers (2003) and the gap estimator of Süveges (Faranda et al., 2017a).
In extreme value statistics, θ is the inverse average duration of consecutive sequences of extreme events and is used to cluster events (Ferro and Segers, 2003). A large θ therefore indicates that event occurrences tend to be isolated, while a low θ indicates that events occur as part of a larger group. In dynamical systems theory, quasistationarity is then defined as $\mathrm{1}/\mathit{\theta}$ (some studies directly use θ as a measure of quasistationarity e.g., Franzke, 2013). A perfectly stationary point of the system (where $\frac{\mathrm{d}\mathit{x}}{\mathrm{d}t}=\mathrm{0}$) has infinite stationarity. By contrast, if trajectories immediately leave the neighborhood of x_{0}, the stationarity is equal to 1. In practice, long trajectories of x(t) are required to explore all possible states of the state space (also called the set of attractor points or simply the attractor) and to best approximates sequences of states on the attractor. Note that the choice of ϵ constrains the values that θ can reach. In particular, very small values (necessary for Eq. 13 to hold) usually lead to $\mathrm{1}/\mathit{\theta}$ being around 1–2 d (e.g., Faranda et al., 2017b, a; Holmberg et al., 2023), meaning that the resulting persistence metric is highly local in time and not necessarily relevant for S2S timescales.
The advantage of the method is that it is grounded in mathematical theory and, unlike approaches based on, e.g., transition probabilities, it does not require categorizing the data. It also provides an easily interpretable index of quasistationarity for each system state (the average residence time of the system's trajectory around that state). Holmberg et al. (2023) have, for instance, taken this approach to analyze the link between atmospheric circulation persistence and warm spells in Europe. Faranda et al. (2019) averaged θ values across all time steps to quantify global persistence in the North Atlantic flow. Going back to our example of the European summer circulation, the three major persistent Z500 anomaly patterns are obtained by clustering the top 3 % (115) most persistent daily patterns (Fig. 8). We chose five cluster centers empirically, after testing for 2–15 centers and selecting the highest number of clusters for which the resulting patterns were subjectively different. As in Fig. 6, persistent states include blocking regimes (over Scandinavia, western Russia and the North Sea; Fig. 8a, b, c), a zonal regime (Fig. 8d) and a European ridge regime (Fig. 8e).
4.2.5 Residence times
A common way to quantify quasistationarity in time series relies on the concept of “residence time”. Residence times extend the dynamical systems approach of the previous section. For a continuous series, the residence time R at time t is defined as the time during which the system remains similar to its value in t:
where $\mathrm{d}(\cdot ,\cdot )$ is a similarity metric and ϵ is a small threshold. For a categorical series, R is similarly defined as the time during which the system remains in its state x(t) before transitioning to another state:
The definition can be relaxed to allow for brief interruptions in a sequence of similar states (one “average” day between two sequences of warm or wet days, for instance) (e.g., Ali et al., 2021). This definition can be extended to sets $\stackrel{\mathrm{\u0303}}{\mathcal{S}}\subset \mathcal{S}$ of several system states: $R\left(t\right)={max}_{k\ge \mathrm{0}}\left\{\mathit{x}(t+{k}^{\prime})\in \stackrel{\mathrm{\u0303}}{\mathcal{S}}\phantom{\rule{0.125em}{0ex}}\forall \phantom{\rule{0.125em}{0ex}}\mathrm{0}\le {k}^{\prime}\le k\right\}$ (Richardson et al., 2019). Note that Eqs. (14)–(15) take an Eulerian perspective; their parallel in the Lagrangian perspective is the concept of “survival time”, i.e., the duration of a specific weather system or pattern along its trajectory (Liu et al., 2018; von Lindheim et al., 2021).
The residence time approach can characterize episodic, state and global quasistationarity and is particularly useful for predictability and risk assessment (De Luca et al., 2019; Francis et al., 2020; Berkovic and RavehRubin, 2022). From the time perspective, large values of R(t) indicate the most quasistationary periods. A minimum threshold is often defined to separate quasistationary from nonstationary periods: for instance, 2 d for weather regimes (De Luca et al., 2019; Francis et al., 2020) and extreme precipitation (Du et al., 2022), 5–6 d for warm spells (Berkovic and RavehRubin, 2022; Rousi et al., 2022b), 5–25 d for geopotential anomalies (Dole and Gordon, 1983), or two seasons for droughts (Ford and Labosier, 2014).
In the state space, the quasistationarity of any state x_{0} can be described from the distribution of its residence times ${R}_{{\mathit{x}}_{\mathrm{0}}}\sim \mathbb{P}\left\{R\left(t\right)\phantom{\rule{0.125em}{0ex}}\mathrm{s}.\mathrm{t}.\phantom{\rule{0.33em}{0ex}},\mathit{x}\left(t\right)={\mathit{x}}_{\mathrm{0}}\phantom{\rule{0.125em}{0ex}}\mathrm{and}\phantom{\rule{0.125em}{0ex}}\mathit{x}(t\mathrm{1})\ne {\mathit{x}}_{\mathrm{0}}\right\}$ (Fig. 9). The easiest is to calculate the mean (Kyselý and Domonkos, 2006; Kučerová et al., 2017; Richardson et al., 2019), maximum (Rousi et al., 2022b) or some extreme percentile (Zolina et al., 2013) of ${R}_{{\mathit{x}}_{\mathrm{0}}}$. Like autocorrelation, it is also possible to characterize quasistationarity by the dependence of $\mathbb{P}\left({R}_{{\mathit{x}}_{\mathrm{0}}}=n\right)$ on n (Sharma and Panu, 2014). Light distribution tails indicate an exponential decay and shortterm memory (see Sect. 4.2.6), while heavy tails (e.g., power law (Bunde et al., 2013) or qexponential (Weber et al., 2019) distributions) point to longrange dependence in the system. For example, Pfleiderer and Coumou (2018) fitted exponential models to the distribution of consecutive warm and cold spells to quantify temperature persistence in the Northern Hemisphere, and Huguenin et al. (2020) did the same for weather types over Central Europe. Many other types of distribution can also be fitted (e.g., Deni et al., 2010; Zolina et al., 2013). Residence times can also be modeled as functions of covariates. For instance, Röthlisberger et al. (2019) and Ali et al. (2021) apply a Weibull regression model to relate the duration of dry and wet spells to Rossby wave activity in the midlatitudes.
Finally, averaging R(t) over time t, or ${R}_{{\mathit{x}}_{\mathrm{0}}}$ across all possible system states, provides a convenient and easily interpretable global quasistationarity metric (e.g., Pfleiderer et al., 2019).
4.2.6 Transition probabilities
Residence times characterize how long the system remains in the same state x_{0}. By contrast, transition probabilities P describe the likelihood that the state of the system does not change at the next time step:
Transition probabilities can be further divided by marginal probabilities to yield probability ratios (e.g., Kolstad et al., 2015):
A high transition probability P(x_{0}) or PR(x_{0}) means that the system likely remains in x_{0}, which indicates persistence (Fig. 10a). Like residence times, transition probabilities can be computed on specific system states only or averaged across all system states (according to their frequency) to yield a measure of global quasistationarity. Their statistical significance can be estimated by comparing observed transition probabilities to those obtained from a large set of randomly shuffled series. However, unlike residence times, transition probabilities cannot identify quasistationary periods: they only work in the state space.
Økland and Lejenäs (1987) applied this method to blocking quasistationarity; Ford and Labosier (2014), Wilby et al. (2016), and Moon et al. (2018) to droughts; Fereday (2017) to largescale sealevel pressure patterns; Guilbert et al. (2015) to daily precipitation; and Röthlisberger and Martius (2019) to dry and warm periods (Fig. 10a).
It is often convenient to assume that the distribution of system states at t+1 only depends on the system state at t, in other words to approximate the system with a firstorder Markov chain model (Sericola, 2013). Under this assumption, transition probabilities are directly related to the distribution of residence times. Indeed, if the probability of remaining in state x_{0} does not depend on how long the system has previously been in state x_{0}: $\mathbb{P}\left(\mathit{x}(t+\mathrm{1})={\mathit{x}}_{\mathrm{0}}\left\mathit{x}\right(t)={\mathit{x}}_{\mathrm{0}}\right)=\mathit{\alpha}$, then the probability that the residence time of x_{0}, ${R}_{{\mathit{x}}_{\mathrm{0}}}$ exceeds n is equal to the probability of finding a sequence of at least n consecutive x_{0} values, which specifically leads to the following equation:
where the subscript in ℙ_{t} specifies that the probability is taken with respect to the time variable, whereas in $\mathbb{P}\left({T}_{{\mathit{x}}_{\mathrm{0}}}\ge n\right)$ it is taken with respect to the ensemble of residence times for x_{0}. Thanks to the firstorder Markov assumption, Eq. (18) factors as
Thus,
where $\mathit{\beta}={\mathbb{P}}_{t}\left(\mathit{x}(t\mathrm{1})\ne {\mathit{x}}_{\mathrm{0}};\mathit{x}\left(t\right)={\mathit{x}}_{\mathrm{0}}\right)$ does not depend on t for a (statistically) stationary system, and
The logarithm of the distribution of residence times is therefore a linear function of n, with the slope (log (α)) equal to the transition probability. The Markov assumption is thus consistent with a shortterm dependent system in which residence time probabilities decay exponentially with n (see Sect. 4.2.5 and Bunde et al., 2013). For instance, Huguenin et al. (2020) fit such exponential laws and use their slope as measures of weather regime quasistationarity over Central Europe in current and future climates. This approach is well suited to comparing two different series (e.g., two different periods or locations) because it separates changes in the marginal frequency of state x_{0} (only affecting β) from changes in quasistationarity (only affecting α).
4.2.7 Time series modeling
Another possibility to investigate state quasistationarity is to fit a statistical model to the data that links successive system values with each other:
where stands for equality in distribution and {θ(t)}_{k=1…t} are covariates. Transition probabilities (Sect. 4.2.6), for example, are equivalent to modeling a system's evolution by a transition matrix between all possible system states. Hidden Markov models (Sect. 4.2.1) likewise estimate transition matrices between hidden states. Logistic regression can be used to assess quasistationarity in a given state. Several studies used it to quantify quasistationarity in droughts, using previous precipitation or soil moisture anomalies as predictor variables (θ(t) in Eq. 22) (Ford and Labosier, 2014; Meng et al., 2017) (Fig. 10b). The link to drought occurrence at the previous time step is implicitly taken into account with previous precipitation or soil moisture anomalies. If previous system values are explicitly included as covariates to the model, one speaks of autologistic regression (Wolters, 2017).
4.3 Running window methods for episodic persistence
The running window approach identifies persistent periods of a given fixed length in continuous or categorical data. It assesses the degree of quasistationarity over fixed time intervals by means of a “similarity index”. Given a time interval $\mathcal{T}=\left\{{t}_{\mathrm{1}},\mathrm{\dots},{t}_{n}\right\}$ of length n and a similarity metric $d(\cdot ,\cdot )$, the similarity index S_{n}(𝒯) is defined as the average similarity between system values in 𝒯. This similarity can be computed with reference to one time step x(t_{k}) in particular (typically the first one):
or it can measure the average similarity across all pairs of values:
Different time intervals can then be ranked according to their degree of quasistationarity. Mo and Ghil (1987), for instance, use pattern correlation as similarity metric and set a threshold of 0.5 to separate quasistationary from nonstationary periods. Another possibility is to keep the N periods with the largest S_{n}(𝒯) (i.e., to use a percentilebased threshold). Quasistationary states can then be detected with, e.g., pattern recognition techniques applied to the most quasistationary periods (Mo and Ghil, 1987). We show an example for summertime European circulation in Fig. 11. The most common are blockingtype patterns over northern Europe and western Russia (Fig. 11a, b, c), but we also find two pronounced zonal patterns with a southward or northwardshifted jet (Fig. 11d, e).
Averaging S_{n}(𝒯) values across multiple periods of length n can additionally provide for a measure of global quasistationarity, as Hoffmann et al. (2021) did to analyze 10 d persistence in summer atmospheric circulation. When working with categorical data, similarity indices can be averaged for each system state to highlight the most quasistationary ones (Horel, 1985b).
The running window approach is well suited to compute temporal trends in quasistationarity, since the similarity index can be defined at all time steps (Hoffmann et al., 2021). One limitation, however, is that results for the same system but calculated for different period lengths are not directly comparable since the marginal distribution of S_{n}(𝒯) depends on n. Furthermore, S_{n}(𝒯) only measures the average similarity between successive values. A high S_{n}(𝒯) does not guarantee a high degree of quasistationarity because there could be breaks in between sequences of similar values. This is especially relevant for long periods during which strict persistence is unlikely to occur.
While we try to give a comprehensive view of common methods used in quasistationarity analysis, we could not include all the methods that exist in the literature. Kornhuber and TamarinBrodsky (2021), for instance, define quasistationarity in a Lagrangian context as the zonal velocity of individual weather systems (in their case, localized temperature anomalies), while Hoskins and Woollings (2015) use Rossby wave phase speed as a proxy for weather quasistationarity in the midlatitudes. Finally, in the framework of dynamical systems theory, quasistationarity can also be equated to the system remaining for some period of time in a small subset of the state space. During that period, the system may explore different configurations, but fewer configurations are explored than if it had been able to evolve freely across the whole state space (Fereday, 2017). In this sense, persistent states are states from which the system takes a (statistically) long time to reach the rest of its state space. This perspective relates to the intransitivity theory of Lorenz (1990), which postulates that the evolution of the atmospheric circulation on S2S timescales can, under certain conditions, be governed by a few wellseparated attractors, each with their own preferred states (Weiland et al., 2021).
We now discuss methods to capture recurrence in atmospheric data (see Table 2 for an overview). In contrast to quasistationarity, most methods here only quantify recurrence in a specific system state. They typically require binary time series as input (where “1” represents an occurrence of the system state or event of interest) and count event occurrences over time. Recurrence is therefore always defined for a given state or event, usually with the potential for large impacts like tropical (Mumby et al., 2011) and extratropical cyclones (Mailier et al., 2006; Vitolo et al., 2009), windstorms (Khare et al., 2015), or extreme precipitation episodes (Kopp et al., 2021). For this reason, we did not explicitly divide this section into global, state and episodic methods. Common methods are available to characterize episodic persistence (Sect. 5.1.1) and the persistence of given states (Sect. 5.1.2–5.2.1). However, few definitions of “global recurrence” exist. Additionally, no study has, to our knowledge, tried to analyze S2S recurrence in a “nonsupervised” way, i.e., that would allow us to objectively identify recurrent states. In Sect. 5.2.2, we discuss the potential of recurrence plots to fill that gap.
Because most recurrence methods work directly with binary time series, point processes are a convenient theoretical tool in recurrence analysis. Point processes are a class of probability models for the random occurrence of points in a space (one or multidimensional space). The most simple hypothesis that can be made is that of complete serial randomness, i.e., points occur completely independently of one another. In this case, recurrence occurs by chance only. The homogeneous Poisson point process is a simple model without memory for binary series in which events occur without correlation and with a constant intensity. In a homogeneous Poisson process, the number of points in individual time intervals are independent, and the number of points in an interval of length τ is Poisson distributed with the rate λτ, where λ is the process intensity rate. This model is useful to simulate binary series with complete serial randomness, and therefore provides an empirical basis to assess the statistical significance of various recurrence metrics in observations.
Separating distinct occurrences of a state, pattern or event is critical to distinguish recurrence from quasistationarity. Over short periods, like a season, it is in principle possible to visually separate multiple systems occurring in close succession, for instance with Hovmöller diagrams (e.g., Tuel et al., 2022a; Rousi et al., 2022a). Most of the time, however, separation cannot be done by hand, and objective criteria are required. One option is to use physically based detection and tracking algorithms – e.g., for cyclones (Hodges, 1995; Wernli and Schwierz, 2006) or atmospheric blocks (Schwierz et al., 2004; Steinfeld, 2021) – that can separate distinct weather systems. A second option is to use a minimum duration between events to classify them as distinct, with the choice depending on the underlying physical system or on impact considerations. Barton et al. (2016), Kopp et al. (2021) and Tuel and Martius (2021a), for instance, require a minimum of 2 d with nonextreme precipitation between extreme precipitation episodes (a typical timescale for extratropical cyclones). In Pfleiderer et al. (2019), a single day is used to separate distinct dry or warm periods.
5.1 Diagnostic methods
5.1.1 Window counts
A simple way to identify recurrent periods is to look for periods with high event counts. One possibility is to set a time window and require the number of event counts during this window to exceed some threshold (two or higher; Fig. 12a). Kopp et al. (2021) and Tuel and Martius (2022a) take this approach to analyze recurrent extreme precipitation episodes at the 2 and 3week timescales. Pinto et al. (2014) likewise identify cyclone clusters over the British Isles as periods with four or more consecutive cyclones within 7 d. Another possibility is to look for clusters of events in which each event is separated by at most n time steps from the preceding one. For instance, Bevacqua et al. (2020) define cyclone clusters as sequences of N≥2 consecutive cyclones transiting within a given area and separated by 24 h or less. This method is largely impact driven and does not aim at modeling the binary series, assessing the significance of the clustering or calculating return periods for large event counts (Dacre and Pinto, 2020). For this, further steps must be taken.
5.1.2 Dispersion metrics
Dispersion metrics characterize the distribution of event occurrences in a time series. In a series with no tendency to clustering, high event counts over a given duration should be much less frequent than in a series where clustering is prevalent. If events occur with the same average probability in the two series, then in the clustered series, sequences with no events should mechanically be more frequent, and the variance in event counts will be higher. In a homogeneous Poisson process (see above) N_{k}(τ), event counts in disjoint intervals of length τ (indexed by k) follow a Poisson distribution with parameter λτ. Hence, the expected value μ(τ) and variance σ(τ)^{2} of N_{k}(τ) are both equal to λτ. Two common statistical dispersion metrics quantify deviations from this homogeneous behavior. The dispersion statistic is a measure of temporal correlation over different timescales. It is defined as Mailier et al. (2006)
ψ is related to the Fano factor ($\mathrm{FF}\left(\mathit{\tau}\right)=\mathit{\psi}\left(\mathit{\tau}\right)+\mathrm{1}$) (Telesca, 2007), also called index of dispersion (Mailier et al., 2006). A positive ψ(τ) (σ(τ)^{2}>μ(τ)) indicates overdispersion, meaning that events tend to occur in clusters. A negative ψ(τ) points to underdispersion, meaning that events occur more regularly than in a random process. Nearzero values are consistent with a homogeneous Poisson process (Fig. 12b).
Similarly, the Allan factor (AF) is defined as the variance of successive event counts over an interval of length τ divided by twice the average event count in τ steps (Telesca, 2007):
For a homogeneous Poisson process, AF(τ)=1, while AF(τ)>1 indicates clustering behavior. ψ and AF can be evaluated for a single window τ^{*}, but assessing how they evolve with τ helps better detect recurrent behavior in the series (Serinaldi and Kilsby, 2013). In the absence of clustering, ψ(τ) and AF(τ) remain constant with τ. In clustered processes, however, both metrics scale with τ, often with a power law behavior over some range [τ_{1},τ_{2}], which can be described by an exponent (Telesca, 2007) similar to the spectral analysis of persistence (Sect. 4.1.2). One advantage of ψ and AF is that they are insensitive to the underlying event frequency, which makes them convenient to compare different locations or event identification algorithms with each other (Dacre and Pinto, 2020). However, the corresponding null hypothesis – that the absence of recurrence is equivalent to exponentially distributed event return times under the Poisson assumption – is not always relevant, and other approaches may be better suited to detect deviations from memoryless processes (Blender et al., 2015).
The statistical significance of ψ(τ) and AF(τ) can be assessed empirically through Monte Carlo sampling of random series without memory or analytically. Mailier et al. (2006) discuss a chisquared test for ψ(τ) (for a fixed window τ), while Serinaldi and Kilsby (2013) introduce the sampling distribution of the empirical AF estimator. Note that seasonality in event counts can artificially inflate the dispersion statistic, meaning that care must be taken either to remove the seasonality (e.g., Tuel and Martius, 2021a) or to simulate random processes with the same seasonality as the observed signal (λ function of the month, for instance). Note also that ψ(τ) and AF(τ) both involve μ(τ) in the denominator; in other words, these metrics are normalized by the marginal event frequency. The same ψ(τ) can be obtained for series with different event frequencies, which makes dispersion metrics less relevant for impacts than other methods (e.g., window counts; Sect. 5.1.1).
Since its introduction by Mailier et al. (2006), the dispersion statistic has been frequently applied in recurrence analysis, namely to detect temporal or serial clustering in cyclones, both tropical (Mumby et al., 2011; Wolff et al., 2016) and extratropical (Vitolo et al., 2009; Pinto et al., 2013, 2014; Economou et al., 2015; Pinto et al., 2016). The Allan factor has been applied to daily precipitation series (Serinaldi and Kilsby, 2013), wave storms (Besio et al., 2017) and to extreme wind speeds in Switzerland (Telesca et al., 2020).
Alternative metrics that characterize the distribution of event counts have been proposed in the literature. For example, Kopp et al. (2021), working with extreme precipitation counts over 2 to 4week timescales, proposed a new metric, S_{cl}(τ) – which they named the “clustering metric” – based on a weighted sum of window event counts. Like the dispersion statistic and Allan factor, S_{cl}(τ) requires choosing a window size τ. Clustering periods of length τ are then identified and ranked by applying a moving window sum to the binary series. The period with the most events is selected first; the corresponding events are removed from the series (to avoid any intersection between clustering periods); and the procedure continues until a minimum number of events is obtained or until no more clustering periods (with at least two events) are found. This yields a set of K periods with decreasing event counts n_{k}(τ): ${n}_{\mathrm{1}}\left(\mathit{\tau}\right)\ge {n}_{\mathrm{2}}\left(\mathit{\tau}\right)\ge \mathrm{\dots}\ge {n}_{K}\left(\mathit{\tau}\right)$. Kopp et al. (2021) then define their clustering metric by
where q_{k} are weights that can be defined in various ways, as long as they verify (i) that q_{k} decreases with k (periods with fewer events are given less weight) and (ii) that $q\left(i\right)q(i+\mathrm{1})>q(i+\mathrm{1})q(i+\mathrm{2})$. In the end, this approach is equivalent to a scoring system for clustering periods, and the resulting metric S_{cl}(τ) correlates well with the dispersion statistic.
5.1.3 Ripley's K function
Given a system state x_{0}, Ripley's K function measures the average number of occurrences of x_{0} within a given neighborhood around a random occurrence of x_{0}. It was originally developed for spatial data, but it can be applied in a temporal setting to characterize recurrence in x_{0}. Let the binary variable Y_{t}=𝟙{X_{t} is in state x_{0}}. Ripley's K function for a neighborhood of size n is then equal to Ripley (1981) and Dixon (2014):
where λ=𝔼[Y_{t}] (frequency of x_{0} in X_{t}). The statistical significance of K(n) can then be assessed by comparing it to values obtained for a homogeneous Poisson process of same length and with the same average rate of occurrence as Y_{t} (Fig. 13). In Eq. (28), λ^{−1} can also be discarded, meaning that K(n) represents an average event count and is more directly related to impacts. Ripley's K was used in meteorology to characterize S2S recurrence in extreme precipitation events by Barton et al. (2016), Tuel and Martius (2021a), and Tuel and Martius (2021b). Stephenson et al. (2004) and Hannachi (2010) also used Ripley's K to characterize the clustering of climate modes in the state space.
5.1.4 Distribution of interevent times
Quasistationarity in a given system state x_{0} can be described by the corresponding distribution of residence times (Sect. 4.2.5). Similarly, the distribution of interevent times (also called “recurrence times” or “waiting times”) characterizes the recurrence in x_{0} from a global perspective (Altmann and Kantz, 2005). The recurrence time T_{i} is the duration between the ith and (i+1)th occurrences of state x_{0} in a time series. We can easily see that both short and long recurrence times will be more frequent in a series that exhibits clustering than in a series where events occur independently of one another. The presence of a heavy tail in the distribution of T_{i} therefore points to a tendency for temporal clustering in the series. If x_{0} is rare enough, the distribution of T_{i} may converge to wellknown distributions like the exponential (a sign of shortterm dependence, as with homogeneous Poisson processes), stretched exponential or power law distribution (Eichner et al., 2007; Corral, 2015). When events correspond to extremes above a given threshold, the exponents that characterize these distribution are even related to the autocorrelation coefficient of the original series (Santhanam and Kantz, 2008; Kalra and Santhanam, 2021).
5.2 Stochastic modeling of recurrence
5.2.1 Event series
We have so far discussed diagnostic methods that characterize recurrence in binary time series. A different approach consists of fitting point processes to the series. A point process is a random process representing the occurrence times of specific events. It can be characterized by event times T(n) (time of the nth event in the series) or by the distribution of event counts N(t_{1},t_{2}) for any interval [t_{1},t_{2}]. We can define its intensity or hazard function λ(t) by
and, for any t_{1}<t_{2}, its intensity measure Λ(t_{1},t_{2}) as
Many different point processes have been described in the statistics literature (e.g., Cox and Isham, 1980), but clustering analyses in weather and climate studies have mainly relied on Poisson processes. In a Poisson process, λ(t) is a deterministic function of time. If λ(t) is constant, the process is said to be homogeneous; otherwise, the process is said to be inhomogeneous. For any t_{1}<t_{2}, N(t_{1},t_{2}) follows a Poisson distribution with rate Λ(t_{1},t_{2}):
and event numbers in disjoint intervals are independent from each other.
The point process approach has mainly been used to relate recurrence in extreme events to atmospheric or climate variability. To that end, λ(t) is assumed to depend on given time covariates X(t). A linear dependence corresponds to Poisson generalized linear models (GLMs):
For instance, Villarini et al. (2011), Tuel and Martius (2022a), and Tuel and Martius (2022b) used Poisson GLMs to relate temporal clustering in heavy precipitation to largescale modes of climate and atmospheric variability. Mailier et al. (2006) did the same for extratropical cyclones in the North Atlantic. The linear constraint can be lifted by using general additive models (GAMs) in which the relationship between λ(t) and the covariates is specified by flexible smooth functions (Hastie, 1992). Villarini et al. (2013) used GAMs to relate flood recurrence in the US Midwest to modes of climate variability like ENSO, and Barton et al. (2022) used GAMs to link temporal clustering in extreme precipitation across Europe to weather regimes in the North Atlantic. Another extension of the Poisson GLM is the Cox regression model (Smith and Karr, 1986), in which
where λ_{0}(t) is a baseline intensity or hazard function and exp (β⋅X(t)) accounts for the timevarying effects of the selected covariates. Examples include Villarini et al. (2013), Mallakpour et al. (2017), and Yang and Villarini (2019), who applied Cox regression to recurrence in floods and extreme precipitation. More complex models are also possible; for instance, Khare et al. (2015) introduced a Poisson mixture model for windstorm clustering, in which λ(t) is expressed as a stochastic function of time modulated by a gamma distribution (Khare et al., 2015). Note that models can be fitted directly by maximizing the likelihood function (Yang and Villarini, 2019) or indirectly by working with window event counts (Mailier et al., 2006; Tuel and Martius, 2022a; Barton et al., 2022).
5.2.2 Recurrence plots
All previously discussed methods share the same limitation: they require binary time series as input, meaning that they are designed to characterize recurrence only in a given system state. It is not possible with such methods to explore recurrence in an “unsupervised” way, i.e., to automatically detect recurrent system states or sequences of states and associated recurrent periods. Recurrence plots (RPs) may offer a solution to this problem.
An RP is a graphical representation of a system's selfsimilarity with time (Marwan et al., 2007). Mathematically, it is defined as a twodimensional binary matrix $\mathbf{R}\in \mathit{\{}\mathrm{0},\mathrm{1}{\mathit{\}}}^{N\times N}$, with time on both axes, that consists of pairwise comparisons of system values:
with $d(\cdot ,\cdot )$ a distance or similarity metric (by convention, $d\left(\mathit{x}\right({t}_{i}),\mathit{x}({t}_{i})=\mathrm{0}$), Θ the Heaviside step function, and ϵ a similarity threshold. If $d\left(\mathit{x}\right({t}_{i}),\mathit{x}({t}_{j}\left)\right)<\mathit{\u03f5}$, then x(t_{j}) and x(t_{j}) are similar, and $\mathbf{R}(i,j)=\mathrm{1}$ represents a recurrence of the system. By definition, R is symmetric and $\forall \phantom{\rule{0.33em}{0ex}}\mathrm{1}\le i\le N,\phantom{\rule{0.33em}{0ex}}\mathbf{R}(i,i)=\mathrm{1}$. A conceptual example is shown in Fig. 15a.
Before calculating R, the system x can first be reduced to a lowerdimensional space (as in Sect. 4.2.2), for example by projecting x onto a set of principal components (Marwan et al., 2007; Mukhin et al., 2022). Note that R depends on the selected similarity metric and on ϵ. For a given system x, R is thus not uniquely defined. Choosing a value for ϵ should be done carefully. ϵ should be small enough to account for real recurrences in the data but not too small to have enough recurrences to analyze. Several rules of thumb have been proposed in the literature (see Marwan et al., 2007, and references therein). For instance, ϵ can be selected by requiring an average recurrence rate of 1 %–10 %. For weather and climate data, ϵ could also be selected based on physical considerations and expert judgment.
RPs are a convenient way to visualize a system's trajectory, especially for multidimensional systems (like circulation fields). RPs capture persistent behavior: vertical lines indicate quasistationarity, while diagonal lines (outside of the main diagonal) indicate recurrence. RPs can thus highlight periods when the trajectory of a system roughly visits the same sequence of states or parts of the state space.
Several measures have been proposed to quantify the presence of specific patterns in RPs. They are known collectively as “recurrence quantification analysis” (Marwan, 2008). Common ones include the recurrence rate:
the determinism, which quantifies the fraction of recurrence points forming diagonal lines of minimum length l_{min}:
where D(l) is the number of diagonal lines of length l; and the laminarity, which quantifies the fraction of recurrence points forming vertical lines of minimum length l_{min}:
where V(l) is the number of vertical lines of length l. All three can be used as global measures of persistence (e.g., RamirezAmaro and FigueroaNazuno, 2006). However, these indices do not a priori discriminate between quasistationarity and persistence (high values of RR, for instance, could be due to either). To do so would require making sure that similar time steps are separated by “0” in the RP (Fig. 15b).
RPs can also detect local recurrence in time. With RPs, one can extend the counting approach introduced in Sect. 5.1.1 to multipleday sequences of complex multidimensional data. Specifically, to search for sequences of p steps that recur over periods of q steps ($\mathrm{0}<p<q$), we can simply define a recurrence index RI_{p,q}(t) as the number of diagonal lines of length p in the submatrix of R indexed by $[t,t+p]\times [t,t+q]$ (Fig. 15b). To make sure that distinct occurrences of the same sequence are separated, we can require diagonals to be separated by a minimum number of days (2 in Fig. 15b). It is also possible to allow for short breaks (e.g., one step) in the diagonals.
RPs have recently gained attention in a range of disciplines that deal with complex systems (Goswami, 2019) but only rarely in environmental sciences. RamirezAmaro and FigueroaNazuno (2006) investigated the recurrence properties of major teleconnection patterns with RPs. Yiou et al. (2018) indirectly relied on RPs to identify the most recurrent states of North Atlantic circulation at the intraseasonal timescale. Adeniji et al. (2018) used RPs to analyze the recurrence characteristics of hourly wind speed in Nigeria, and Ray et al. (2019) did so to investigate daily temperature and humidity data across India. Recently, Mukhin et al. (2022) used RPs to detect weather regimes in the Northern Hemisphere. In Fig. 16, we show results obtained for 3 d sequences of similar upperlevel circulation patterns over Europe in summer. We calculate RI_{p,q}(t) with p=3 d and q=21 d and identify the two most common sequences among the ≈30 21 d windows with ${\mathrm{RI}}_{p,q}\left(t\right)\ge \mathrm{2}$. The method identifies recurring blocks and cyclones over Scandinavia as the most frequent recurring sequences of circulation patterns. The results appear physically meaningful since counts of individual blocks and cyclones during the corresponding windows range from 1–2 (for blocks) and 4–11 (for cyclones). More accurate patterns could possibly be obtained by choosing more than two clusters, but the results of Fig. 16 are only meant to illustrate the potential of the method.
Weather persistence on S2S timescales has been a topic of research since the early days of meteorology. Quasistationary or recurrent behavior are common features of weather dynamics and are strongly related to fundamental physical processes, weather predictability and surface weather impacts. Studying weather persistence is therefore important for theoretical and practical reasons. One challenge is that persistence remains a very broad concept that relates to different behaviors in dynamical systems. We propose a typology/structure for the broad concepts related to persistence. Namely that persistence is used to describe the average behavior of a system across its whole trajectory (global persistence), to refer to specific segments of this trajectory (episodic persistence) or to qualify the behavior of particular system states (state persistence).
A wide range of methods has been introduced in the literature to describe persistence in weather and climate series, and several exist for each type of persistence. They offer many distinct and often complementary perspectives on persistence. Some methods quantify persistence in a statistical framework, which can be very relevant for weather forecasting. Others focus instead on persistent periods, including statistical flukes; this is a useful approach for risk assessment, but it is one that says nothing about the overall behavior of the series. Other methods still aim to identify which weather patterns tend to be more persistent than others.
The diversity of existing methods presented in this review reflects the fact that persistence is a multifaceted concept. While we can agree on a general definition of the concept, many options exist when it comes to actually quantify persistence in realworld data. Though different methods may be related, each sheds light on a particular aspect of persistence. What is meant by persistence cannot be dissociated from the metric used to quantify it. The choice of method should be guided by the end goal, whether it is process understanding, risk and impact assessment, or predictability. Future research should nevertheless perhaps consider testing the robustness of results to the choice of persistence metric more systematically. This could be particularly important for better characterizing potential trends in persistence under climate change and their associated impacts. As a final note, while we centered our review on S2S persistence, most of the methods we discussed apply in principle to other timescales, be they subdaily or multidecadal. What really differs is how to define the system to analyze (daily time steps may not be relevant for interannual variability, for instance) and interpret persistence.
R code implementing several of the methods presented in this paper is available at https://github.com/Quriosity129/persistence (Tuel, 2023). Code to calculate the dynamical system metrics (Sect. 4.2.4) is available for R at https://github.com/thaos/dtheta (Thao, 2021) or from the R package CSTools
(PerezZanon et al., 2022) and for Python at https://github.com/yrobink/CDSK (Robin, 2021). The structural similarity index (SSIM) can be calculated with the Python package scikitimage
(van der Walt et al., 2014).
We illustrate some of the methods (Figs. 3, 6, 8, 11 and 16) with data from the ERA5 reanalysis (Hersbach et al., 2020; available from https://doi.org/10.24381/cds.bd0915c6). We use geopotential height at 500 hPa and 2 m temperature for the months of June, July and August and over the 1979–2020 period. We remove their seasonality and longterm trends by normalizing the data with a moving 30 d, 7year window (as in Pfleiderer and Coumou, 2018). The calculation of the blocking and cyclone indices in Fig. 16 is described in Tuel and Martius (2022a).
AT: conceptualization, formal analysis, visualization, software, writing – original draft. OM: conceptualization, funding acquisition, supervision, writing – review and editing.
At least one of the (co)authors is a member of the editorial board of Earth System Dynamics. The peerreview process was guided by an independent editor, and the authors also have no other competing interests to declare.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This research has been supported by the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. 178751).
This paper was edited by Christian Franzke and reviewed by Abdel Hannachi and two anonymous referees.
Adeniji, A. E., Olusola, O. I., and Njah, A. N.: Comparative study of chaotic features in hourly wind speed using recurrence quantification analysis, Aip. Adv., 8, 025102, https://doi.org/10.1063/1.4998674, 2018. a
Ali, S. K., Aydam, Z. M., and Rashed, B. M.: Similarity metrics for classification: A Review, Iop. Conf. Ser.Mat. Sci., 928, 032052, https://doi.org/10.1088/1757899X/928/3/032052, 2020. a
Ali, S. M., Martius, O., and Röthlisberger, M.: Recurrent Rossby Wave Packets Modulate the Persistence of Dry and Wet Spells Across the Globe, Geophys. Res. Lett., 48, e2020GL091452, https://doi.org/10.1029/2020GL091452, 2021. a, b, c, d, e, f, g
Ali, S. M., Röthlisberger, M., Parker, T., Kornhuber, K., and Martius, O.: Recurrent Rossby waves and southeastern Australian heatwaves, Weather Clim. Dynam., 3, 1139–1156, https://doi.org/10.5194/wcd311392022, 2022. a
Altmann, E. G. and Kantz, H.: Recurrence time analysis, longterm correlations, and extreme events, Phys. Rev. E, 71, 056106, https://doi.org/10.1103/PhysRevE.71.056106, 2005. a, b
Ault, T. R., Cole, J. E., Overpeck, J. T., Pederson, G. T., and Meko, D. M.: Assessing the Risk of Persistent Drought Using Climate Model Simulations and Paleoclimate Data, J. Clim., 27, 7529–7549, https://doi.org/10.1175/JCLID1200282.1, 2014. a
Barnes, E. A. and Hartmann, D. L.: Dynamical Feedbacks and the Persistence of the NAO, J. Atmos. Sci., 67, 851–865, https://doi.org/10.1175/2009JAS3193.1, 2010. a
Barnston, A. G. and Livezey, R. E.: Classification, Seasonality and Persistence of LowFrequency Atmospheric Circulation Patterns, Mon. Weather Rev., 115, 1083–1126, https://doi.org/10.1175/15200493(1987)115<1083:CSAPOL>2.0.CO;2, 1987. a, b
Barton, Y., Giannakaki, P., von Waldow, H., Chevalier, C., Pfahl, S., and Martius, O.: Clustering of RegionalScale Extreme Precipitation Events in Southern Switzerland, Mon. Weather Rev., 144, 347–369, https://doi.org/10.1175/MWRD150205.1, 2016. a, b, c, d
Barton, Y., Rivoire, P., Koh, J., Ali, M. S., Kopp, J., and Martius, O.: On the temporal clustering of European extreme precipitation events and its relationship to persistent and transient largescale atmospheric drivers, Weather Clim. Extrem., 38, 100518, https://doi.org/10.1016/j.wace.2022.100518, 2022. a, b, c, d
Baur, F.: ExtendedRange Weather Forecasting, American Meteorological Society, Boston, MA, 814–833, https://doi.org/10.1007/9781940033709_66, 1951. a
Beran, J.: Statistics for LongMemory Processes, Routledge, 315 pp., https://doi.org/10.1201/9780203738481, 2017. a
Berkovic, S. and RavehRubin, S.: Persistent warm and dry extremes over the eastern Mediterranean during winter: The role of North Atlantic blocking and central Mediterranean cyclones, Q. J. Roy. Meteor. Soc., 148, 2384–2409, https://doi.org/10.1002/qj.4308, 2022. a, b
Besio, G., Briganti, R., Romano, A., Mentaschi, L., and De Girolamo, P.: Time clustering of wave storms in the Mediterranean Sea, Nat. Hazards Earth Syst. Sci., 17, 505–514, https://doi.org/10.5194/nhess175052017, 2017. a
Bevacqua, E., Zappa, G., and Shepherd, T. G.: Shorter cyclone clusters modulate changes in European wintertime precipitation extremes, Environ. Res. Lett., 15, 124005, https://doi.org/10.1088/17489326/abbde7, 2020. a, b, c, d
Black, E., Blackburn, M., Harrison, G., Hoskins, B., and Methven, J.: Factors contributing to the summer 2003 European heatwave, Weather, 59, 217–223, https://doi.org/10.1256/wea.74.04, 2004. a, b
Blanchet, J., Stalla, S., and Creutin, J.D.: Analogy of multiday sequences of atmospheric circulation favoring large rainfall accumulation over the French Alps, Atmos. Sci. Lett., 19, e809, https://doi.org/10.1002/asl.809, 2018. a
Blender, R., Raible, C. C., and Lunkeit, F.: Nonexponential return time distributions for vorticity extremes explained by fractional Poisson processes, Q. J. Roy. Meteor. Soc., 141, 249–257, https://doi.org/10.1002/qj.2354, 2015. a
Bloomberg: Extreme Heat to Persist in India for Third Straight Month, https://www.bloomberg.com/ (last access: 8 September 2023), 2022. a
Bray, M. T. and Cavallo, S. M.: Characteristics of longtrack tropopause polar vortices, Weather Clim. Dynam., 3, 251–278, https://doi.org/10.5194/wcd32512022, 2022. a
Bunde, A., Havlin, S., KoscielnyBunde, E., and Schellnhuber, H. J.: Atmospheric Persistence Analysis: Novel Approaches and Applications, pp. 170–191, Springer Berlin Heidelberg, Berlin, Heidelberg, https://doi.org/10.1007/9783642562570_5, 2002. a
Bunde, A., Büntgen, U., Ludescher, J., Luterbacher, J., and von Storch, H.: Is there memory in precipitation?, Nat. Clim. Change, 3, 174–175, https://doi.org/10.1038/nclimate1830, 2013. a, b, c
Büeler, D., Ferranti, L., Magnusson, L., Quinting, J. F., and Grams, C. M.: Yearround subseasonal forecast skill for Atlantic–European weather regimes, Q. J. Roy. Meteor. Soc., 147, 4283–4309, https://doi.org/10.1002/qj.4178, 2021. a
Cassou, C., Terray, L., and Phillips, A. S.: Tropical Atlantic Influence on European Heat Waves, J. Clim., 18, 2805–2811, https://doi.org/10.1175/JCLI3506.1, 2005. a
Chapman, C. C., Monselesan, D. P., Risbey, J. S., Feng, M., and Sloyan, B. M.: A largescale view of marine heatwaves revealed by archetype analysis, Nat. Commun., 13, 7843, https://doi.org/10.1038/s4146702235493x, 2022. a
Charney, J. G. and DeVore, J. G.: Multiple Flow Equilibria in the Atmosphere and Blocking, J. Atmos. Sci, 36, 1205–1216, https://doi.org/10.1175/15200469(1979)036<1205:MFEITA>2.0.CO;2, 1979. a
Corral, A.: Scaling in the timing of extreme events, Chaos Soliton Fract., 74, 99–112, https://doi.org/10.1016/j.chaos.2015.01.011, 2015. a
Cox, D. and Isham, V.: Point Processes, Routledge, 188 pp., https://doi.org/10.1201/9780203743034, 1980. a
Dacre, H. F. and Pinto, J. G.: Serial clustering of extratropical cyclones: a review of where, when and why it occurs, npj Clim. Atmos. Sci., 3, 48, https://doi.org/10.1038/s41612020001529, 2020. a, b, c, d, e
De las Nieves López García, M. and Requena, J. P. R.: Different methodologies and uses of the Hurst exponent in econophysics, Estud. Econ., 37, 96–108, 2019. a
De Luca, P., Harpham, C., Wilby, R. L., Hillier, J. K., Franzke, C. L. E., and Leckebusch, G. C.: Past and Projected Weather Pattern Persistence with Associated MultiHazards in the British Isles, Atmosphere, 10, 577, https://doi.org/10.3390/atmos10100577, 2019. a, b
Degenhardt, L. and Ólafsson, H.: Persistence of observed air temperatures in Iceland, Int. J. Climatol., 39, 1262–1275, https://doi.org/10.1002/joc.5875, 2019. a
DelSole, T.: Optimally Persistent Patterns in TimeVarying Fields, J. Atmos. Sci., 58, 1341–1356, https://doi.org/10.1175/15200469(2001)058<1341:OPPITV>2.0.CO;2, 2001. a, b, c, d, e
Demuzere, M., Kassomenos, P., and Philipp, A.: The COST733 circulation type classification software: an example for surface ozone concentrations in Central Europe, Theor. Appl. Climatol., 105, 143–166, https://doi.org/10.1007/s0070401003784, 2011. a
Deni, S. M., Jemain, A. A., and Ibrahim, K.: The best probability models for dry and wet spells in Peninsular Malaysia during monsoon seasons, Int. J. Climatol., 30, 1194–1205, https://doi.org/10.1002/joc.1972, 2010. a
Desgraupes, B.: clusterCrit: An R Package for Computing Clustering Quality Indices, R package version 1.2.8, https://CRAN.Rproject.org/package=clusterCrit (last access: 8 September 2023), 2018. a
Di Capua, G., Sparrow, S., Kornhuber, K., Rousi, E., Osprey, S., Wallom, D., van den Hurk, B., and Coumou, D.: Drivers behind the summer 2010 wave train leading to Russian heatwave and Pakistan flooding, npj Clim. Atmos. Sci., 4, 55, https://doi.org/10.1038/s41612021002119, 2021. a, b
Di Lorenzo, E. and Mantua, N.: Multiyear persistence of the 2014/15 North Pacific marine heatwave, Nat. Clim. Change, 6, 1042–1047, https://doi.org/10.1038/nclimate3082, 2016. a
Dixon, P. M.: Ripley's K Function, in: Wiley StatsRef: Statistics Reference Online, John Wiley & Sons, Ltd, Chichester, UK, https://doi.org/10.1002/9781118445112.stat07751, 2014. a
Dole, R. M. and Gordon, N. D.: Persistent Anomalies of the Extratropical Northern Hemisphere Wintertime Circulation: Geographical Distribution and Regional Persistence Characteristics, Mon. Weather Rev., 111, 1567–1586, https://doi.org/10.1175/15200493(1983)111<1567:PAOTEN>2.0.CO;2, 1983. a, b, c
Domeisen, D. I., White, C. J., AfarganGerstman, H., Muñoz, Á. G., Janiga, M. A., Vitart, F., Wulff, C. O., Antoine, S., Ardilouze, C., Batté, L., Bloomfield, H. C., Brayshaw, D. J., Camargo, S. J., CharltonPérez, A., Collins, D., Cowan, T., del Mar Chaves, M., Ferranti, L., Gómez, R., González, P. L., González Romero, C., Infanti, J. M., Karozis, S., Kim, H., Kolstad, E. W., LaJoie, E., Lledó, L., Magnusson, L., Malguzzi, P., ManriqueSuñén, A., Mastrangelo, D., Materia, S., Medina, H., Palma, L., Pineda, L. E., Sfetsos, A., Son, S.W., Soret, A., Strazzo, S., and Tian, D.: Advances in the subseasonal prediction of extreme events: Relevant case studies across the globe, B. Am. Meteorol. Soc., E1473–E1501, https://doi.org/10.1175/BAMSD200221.1, 2022. a
Drouard, M. and Woollings, T.: Contrasting Mechanisms of Summer Blocking Over Western Eurasia, Geophys. Res. Lett., 45, 12,040–12,048, https://doi.org/10.1029/2018GL079894, 2018. a, b, c
Du, H., Alexander, L. V., Donat, M. G., Lippmann, T., Srivastava, A., Salinger, J., Kruger, A., Choi, G., He, H. S., Fujibe, F., Rusticucci, M., Nandintsetseg, B., Manzanas, R., Rehman, S., Abbas, F., Zhai, P., Yabi, I., Stambaugh, M. C., Wang, S., Batbold, A., Oliveira, P. T., Adrees, M., Hou, W., Zong, S., Santos e Silva, C. M., Lucio, P. S., and Wu, Z.: Precipitation From Persistent Extremes is Increasing in Most Regions and Globally, Geophys. Res. Lett., 46, 6041–6049, https://doi.org/10.1029/2019GL081898, 2019. a
Du, H., Donat, M. G., Zong, S., Alexander, L. V., Manzanas, R., Kruger, A., Choi, G., Salinger, J., He, H. S., Li, M.H., Fujibe, F., Nandintsetseg, B., Rehman, S., Abbas, F., Rusticucci, M., Srivastava, A., Zhai, P., Lippmann, T., Yabi, I., Stambaugh, M. C., Wang, S., Batbold, A., de Oliveira, P. T., Adrees, M., Hou, W., e Silva, C. M. S., Lucio, P. S., and Wu, Z.: Extreme Precipitation on Consecutive Days Occurs More Often in a Warming Climate, B. Am. Meteorol. Soc., 103, E1130–E1145, https://doi.org/10.1175/BAMSD210140.1, 2022. a
Economou, T., Stephenson, D. B., Pinto, J. G., Shaffrey, L. C., and Zappa, G.: Serial clustering of extratropical cyclones in a multimodel ensemble of historical and future simulations, Q. J. Roy. Meteor. Soc., 141, 3076–3087, https://doi.org/10.1002/qj.2591, 2015. a
Eichner, J. F., KoscielnyBunde, E., Bunde, A., Havlin, S., and Schellnhuber, H.J.: Powerlaw persistence and trends in the atmosphere: A detailed study of long temperature records, Phys. Rev. E, 68, 046133, https://doi.org/10.1103/PhysRevE.68.046133, 2003. a
Eichner, J. F., Kantelhardt, J. W., Bunde, A., and Havlin, S.: Statistics of return intervals in longterm correlated records, Phys. Rev. E, 75, 011128, https://doi.org/10.1103/PhysRevE.75.011128, 2007. a
Faranda, D., Messori, G., AlvarezCastro, M. C., and Yiou, P.: Dynamical properties and extremes of Northern Hemisphere climate fields over the past 60 years, Nonlinear Proc. Geophys., 24, 713–725, https://doi.org/10.5194/npg247132017, 2017a. a, b, c, d
Faranda, D., Messori, G., and Yiou, P.: Dynamical proxies of North Atlantic predictability and extremes, Sci. Rep.UK, 7, 41278, https://doi.org/10.1038/srep41278, 2017b. a, b, c, d, e, f, g
Faranda, D., AlvarezCastro, M. C., Messori, G., Rodrigues, D., and Yiou, P.: The hammam effect or how a warm ocean enhances large scale atmospheric predictability, Nat. Commun., 10, 1316, https://doi.org/10.1038/s41467019093058, 2019. a
Fereday, D.: How Persistent Are North AtlanticEuropean Sector Weather Regimes?, J. Clim., 30, 2381–2394, https://doi.org/10.1175/JCLID160328.1, 2017. a, b, c, d, e, f
Ferro, C. A. T. and Segers, J.: Inference for clusters of extreme values, J. Roy. Stat. Soc. Ser. B, 65, 545–556, https://doi.org/10.1111/14679868.00401, 2003. a, b
Floodlist: Australia? More Floods in Queensland After Widespread Heavy Rainfall, https://floodlist.com/australia/queenslandfloodsmay2022 (last access: 8 September 2023), 2022. a
Ford, T. and Labosier, C. F.: Spatial patterns of drought persistence in the Southeastern United States, Int. J. Climatol., 34, 2229–2240, https://doi.org/10.1002/joc.3833, 2014. a, b, c, d, e
Fraedrich, K. and Larnder, C.: Scaling regimes of composite rainfall time series, Tellus A, 45, 289–298, https://doi.org/10.3402/tellusa.v45i4.14893, 1993. a, b
Francis, J. A., Skific, N., and Vavrus, S. J.: North American Weather Regimes Are Becoming More Persistent: Is Arctic Amplification a Factor?, Geophys. Res. Lett., 45, 11414–11422, https://doi.org/10.1029/2018GL080252, 2018. a, b, c, d
Francis, J. A., Skific, N., and Vavrus, S. J.: Increased persistence of largescale circulation regimes over Asia in the era of amplified Arctic warming, past and future, Sci. Rep.UK, 10, 14953, https://doi.org/10.1038/s41598020719454, 2020. a, b, c, d, e, f
Franzke, C., Majda, A. J., and VandenEijnden, E.: LowOrder Stochastic Mode Reduction for a Realistic Barotropic Model Climate, J. Atmos. Sci., 62, 1722–1745, https://doi.org/10.1175/JAS3438.1, 2005. a
Franzke, C., Crommelin, D., Fischer, A., and Majda, A. J.: A Hidden Markov Model Perspective on Regimes and Metastability in Atmospheric Flows, J. Clim., 21, 1740–1757, https://doi.org/10.1175/2007JCLI1751.1, 2008. a
Franzke, C., Woollings, T., and Martius, O.: Persistent Circulation Regimes and Preferred Regime Transitions in the North Atlantic, J. Atmos. Sci., 68, 2809–2825, https://doi.org/10.1175/JASD11046.1, 2011. a, b, c, d
Franzke, C. L. E.: Persistent regimes and extreme events of the North Atlantic atmospheric circulation, Phil. T. R. Soc. A, 371, 20110471, https://doi.org/10.1098/rsta.2011.0471, 2013. a, b, c, d, e, f, g
Franzke, C. L. E., Barbosa, S., Blender, R., Fredriksen, H.B., Laepple, T., Lambert, F., Nilsen, T., Rypdal, K., Rypdal, M., Scotto, M. G., Vannitsem, S., Watkins, N. W., Yang, L., and Yuan, N.: The Structure of Climate Variability Across Scales, Rev. Geophys., 58, e2019RG000657, https://doi.org/10.1029/2019RG000657, 2020. a, b
Gálfi, V. M., Lucarini, V., and Wouters, J.: A large deviation theorybased analysis of heat waves and cold spells in a simplified model of the general circulation of the atmosphere, J. Stat. Mech.Theory E, 2019, 033404, https://doi.org/10.1088/17425468/ab02e8, 2019. a
GarcíaHerrera, R., Díaz, J., Trigo, R. M., Luterbacher, J., and Fischer, E. M.: A Review of the European Summer Heat Wave of 2003, Crit. Rev. Env. Sci. Tec., 40, 267–306, https://doi.org/10.1080/10643380802238137, 2010. a
Gershunov, A. and Barnett, T. P.: ENSO Influence on Intraseasonal Extreme Rainfall and Temperature Frequencies in the Contiguous United States: Observations and Model Results, J. Clim., 11, 1575–1586, https://doi.org/10.1175/15200442(1998)011<1575:EIOIER>2.0.CO;2, 1998. a
Ghil, M. and Robertson, A. W.: “Waves” vs. “particles” in the atmosphere's phase space: A pathway to longrange forecasting?, P. Natl. Acad. Sci. USA, 99, 2493–2500, https://doi.org/10.1073/pnas.012580899, 2002. a
Goswami, B.: A Brief Introduction to Nonlinear Time Series Analysis and Recurrence Plots, Vibration, 2, 332–368, https://doi.org/10.3390/vibration2040021, 2019. a
Grams, C. M., Beerli, R., Pfenninger, S., Staffell, I., and Wernli, H.: Balancing Europe's windpower output through spatial deployment informed by weather regimes, Nat. Clim. Change, 7, 557–562, https://doi.org/10.1038/nclimate3338, 2017. a, b, c, d, e, f, g
Guilbert, J., Betts, A. K., Rizzo, D. M., Beckage, B., and Bomblies, A.: Characterization of increased persistence and intensity of precipitation in the northeastern United States, Geophys. Res. Lett., 42, 1888–1893, https://doi.org/10.1002/2015GL063124, 2015. a
Haines, K. and Hannachi, A.: Weather Regimes in the Pacific from a GCM, J. Atmos. Sci., 52, 2444–2462, https://doi.org/10.1175/15200469(1995)052<2444:WRITPF>2.0.CO;2, 1995. a
Hamidieh, K., Stoev, S., and Michailidis, G.: On the Estimation of the Extremal Index Based on Scaling and Resampling, J. Comput. Graph. Stat., 18, 731–755, https://doi.org/10.1198/jcgs.2009.08065, 2009. a
Hannachi, A.: LowFrequency Variability in a GCM: ThreeDimensional Flow Regimes and Their Dynamics, J. Clim., 10, 1357–1379, https://doi.org/10.1175/15200442(1997)010<1357:LFVIAG>2.0.CO;2, 1997. a
Hannachi, A.: A New Set of Orthogonal Patterns in Weather and Climate: Optimally Interpolated Patterns, J. Clim., 21, 6724–6738, https://doi.org/10.1175/2008JCLI2328.1, 2008. a, b, c
Hannachi, A.: On the Origin of PlanetaryScale Extratropical Winter Circulation Regimes, J. Atmos. Sci., 67, 1382–1401, https://doi.org/10.1175/2009JAS3296.1, 2010. a
Hannachi, A.: Intermittency, autoregression and censoring: a firstorder AR model for daily precipitation, Meteorol. Appl., 21, 384–397, https://doi.org/10.1002/met.1353, 2014. a
Hannachi, A.: Patterns Identification and Data Mining in Weather and Climate, Springer Nature, https://doi.org/10.1007/9783030670733, 2021. a, b, c
Hannachi, A. and Trendafilov, N.: Archetypal Analysis: Mining Weather and Climate Extremes, J. Clim., 30, 6927–6944, https://doi.org/10.1175/JCLID160798.1, 2017. a
Hannachi, A., Woollings, T., and Fraedrich, K.: The North Atlantic jet stream: a look at preferred positions, paths and transitions, Q. J. Roy. Meteor. Soc., 138, 862–877, https://doi.org/10.1002/qj.959, 2012. a, b
Hannachi, A., Straus, D. M., Franzke, C. L. E., Corti, S., and Woollings, T.: Lowfrequency nonlinearity and regime behavior in the Northern Hemisphere extratropical atmosphere, Rev. Geophys., 55, 199–234, https://doi.org/10.1002/2015RG000509, 2017. a, b, c, d, e, f
Hastie, T. J.: Generalized Additive Models, in Statistical Models in S, Routledge, 59 pp., https://doi.org/10.1201/97802037385357, 1992. a
Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., MuñozSabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J. N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, https://doi.org/10.1002/qj.3803, 2020. a, b
Hodges, K. I.: Feature Tracking on the Unit Sphere, Mon. Weather Rev., 123, 3458–3465, https://doi.org/10.1175/15200493(1995)123<3458:FTOTUS>2.0.CO;2, 1995. a
Hoffmann, P., Lehmann, J., Fallah, B., and Hattermann, F. F.: Atmosphere similarity patterns in boreal summer show an increase of persistent weather conditions connected to hydroclimatic risks, Sci. Rep.UK, 11, 22893, https://doi.org/10.1038/s4159802101808z, 2021. a, b, c, d, e, f, g, h, i, j
Holešovský, J. and Fusek, M.: Improved interexceedancetimesbased estimator of the extremal index using truncated distribution, Extremes, 25, 695–720, https://doi.org/10.1007/s10687022004448, 2022. a
Holmberg, E., Messori, G., Caballero, R., and Faranda, D.: The link between European warmtemperature extremes and atmospheric persistence, Earth Syst. Dynam., 14, 737–765, https://doi.org/10.5194/esd147372023, 2023. a, b
Horel, J. D.: Persistence of the 500 mb Height Field during Northern Hemisphere Winter, Mon. Weather Rev., 113, 2030–2042, https://doi.org/10.1175/15200493(1985)113<2030:POTMHF>2.0.CO;2, 1985a. a
Horel, J. D.: Persistence of Wintertime 500 mb Height Anomalies over the Central Pacific, Mon. Weather Rev., 113, 2043–2048, https://doi.org/10.1175/15200493(1985)113<2043:POWMHA>2.0.CO;2, 1985b. a
Horton, P., Jaboyedoff, M., and Obled, C.: Global Optimization of an Analog Method by Means of Genetic Algorithms, Mon. Weather Rev., 145, 1275–1294, https://doi.org/10.1175/MWRD160093.1, 2017. a
Hoskins, B. and Woollings, T.: Persistent Extratropical Regimes and Climate Extremes, Curr. Clim. Change Rep., 1, 115–124, https://doi.org/10.1007/s4064101500208, 2015. a
Huguenin, M. F., Fischer, E. M., Kotlarski, S., Scherrer, S. C., Schwierz, C., and Knutti, R.: Lack of Change in the Projected Frequency and Persistence of Atmospheric Circulation Types Over Central Europe, Geophys. Res. Lett., 47, e2019GL086132, https://doi.org/10.1029/2019GL086132, 2020. a, b, c
Huntingford, C., Marsh, T., Scaife, A. A., Kendon, E. J., Hannaford, J., Kay, A. L., Lockwood, M., Prudhomme, C., Reynard, N. S., Parry, S., Lowe, J. A., Screen, J. A., Ward, H. C., Roberts, M., Stott, P. A., Bell, V. A., Bailey, M., Jenkins, A., Legg, T., Otto, F. E. L., Massey, N., Schaller, N., Slingo, J., and Allen, M. R.: Potential influences on the United Kingdom's floods of winter 2013/14, Nat. Clim. Change, 4, 769–777, https://doi.org/10.1038/nclimate2314, 2014. a, b
Hurst, H. E.: LongTerm Storage Capacity of Reservoirs, T. Am. Soc. Civ. Eng., 116, 770–799, https://doi.org/10.1061/TACEAT.0006518, 1951. a
Huth, R., Beck, C., Philipp, A., Demuzere, M., Ustrnul, Z., Cahynová, M., Kyselý, J., and Tveito, O. E.: Classifications of Atmospheric Circulation Patterns, Ann. Ny. Acad. Sci., 1146, 105–152, https://doi.org/10.1196/annals.1446.019, 2008. a, b
Kalra, D. S. and Santhanam, M. S.: Inferring long memory using extreme events, Chaos, 31, 113131, https://doi.org/10.1063/5.0064432, 2021. a
Khare, S., Bonazzi, A., Mitas, C., and Jewson, S.: Modelling clustering of natural hazard phenomena and the effect on re/insurance loss perspectives, Nat. Hazards Earth Syst. Sci., 15, 1357–1370, https://doi.org/10.5194/nhess1513572015, 2015. a, b, c
Kimoto, M. and Ghil, M.: Multiple Flow Regimes in the Northern Hemisphere Winter. Part I: Methodology and Hemispheric Regimes, J. Atmos. Sci., 50, 2625–2644, https://doi.org/10.1175/15200469(1993)050<2625:MFRITN>2.0.CO;2, 1993. a
Kolstad, E. W., Sobolowski, S. P., and Scaife, A. A.: Intraseasonal Persistence of European Surface Temperatures, J. Clim., 28, 5365–5374, https://doi.org/10.1175/JCLID150053.1, 2015. a, b
Kolstad, E. W., Barnes, E. A., and Sobolowski, S. P.: Quantifying the role of landatmosphere feedbacks in mediating nearsurface temperature persistence, Q. J. Roy. Meteor. Soc., 143, 1620–1631, https://doi.org/10.1002/qj.3033, 2017. a, b
Kooperberg, C. and O'sullivan, F.: Predictive Oscillation Patterns: A Synthesis of Methods for SpatialTemporal Decomposition of Random Fields, J. Am. Stat. Assoc., 91, 1485–1496, https://doi.org/10.1080/01621459.1996.10476716, 1996. a
Kopp, J., Rivoire, P., Ali, S. M., Barton, Y., and Martius, O.: A novel method to identify subseasonal clustering episodes of extreme precipitation events and their contributions to large accumulation periods, Hydrol. Earth Syst. Sci., 25, 5153–5174, https://doi.org/10.5194/hess2551532021, 2021. a, b, c, d, e, f, g, h, i
Kornhuber, K. and TamarinBrodsky, T.: Future Changes in Northern Hemisphere Summer Weather Persistence Linked to Projected Arctic Warming, Geophys. Res. Lett., 48, e2020GL091603, https://doi.org/10.1029/2020GL091603, 2021. a, b
Kornhuber, K., Osprey, S., Coumou, D., Petri, S., Petoukhov, V., Rahmstorf, S., and Gray, L.: Extreme weather events in early summer 2018 connected by a recurrent hemispheric wave7 pattern, Environ. Res. Lett., 14, 054002, https://doi.org/10.1088/17489326/ab13bf, 2019. a
KoscielnyBunde, E., Bunde, A., Havlin, S., Roman, H. E., Goldreich, Y., and Schellnhuber, H.J.: Indication of a Universal Persistence Law Governing Atmospheric Variability, Phys. Rev. Lett., 81, 729–732, https://doi.org/10.1103/PhysRevLett.81.729, 1998. a
Kossin, J. P.: A global slowdown of tropicalcyclone translation speed, Nature, 558, 104–107, https://doi.org/10.1038/s4158601801583, 2018. a
Koutsoyiannis, D.: Climate change, the Hurst phenomenon, and hydrological statistics, Hydrol. Sci. J., 48, 3–24, https://doi.org/10.1623/hysj.48.1.3.43481, 2003. a
Kučerová, M., Beck, C., Philipp, A., and Huth, R.: Trends in frequency and persistence of atmospheric circulation types over Europe derived from a multitude of classifications, Int. J. Climatol., 37, 2502–2521, https://doi.org/10.1002/joc.4861, 2017. a
Kumar, S., Merwade, V., Kinter, J. L., and Niyogi, D.: Evaluation of Temperature and Precipitation Trends and LongTerm Persistence in CMIP5 TwentiethCentury Climate Simulations, J. Clim., 26, 4168–4185, https://doi.org/10.1175/JCLID1200259.1, 2013. a
Kyselý, J. and Domonkos, P.: Recent increase in persistence of atmospheric circulation over Europe: comparison with longterm variations since 1881, Int. J. Climatol., 26, 461–483, https://doi.org/10.1002/joc.1265, 2006. a
Lawrence, Z. D., Perlwitz, J., Butler, A. H., Manney, G. L., Newman, P. A., Lee, S. H., and Nash, E. R.: The Remarkably Strong Arctic Stratospheric Polar Vortex of Winter 2020: Links to RecordBreaking Arctic Oscillation and Ozone Loss, J. Geophys. Res.Atmos, 125, e2020JD033271, https://doi.org/10.1029/2020JD033271, 2020. a
Legras, B. and Ghil, M.: Persistent Anomalies, Blocking and Variations in Atmospheric Predictability, J. Atmos. Sci., 42, 433–471, https://doi.org/10.1175/15200469(1985)042<0433:PABAVI>2.0.CO;2, 1985. a
Li, J. and Thompson, D. W. J.: Widespread changes in surface temperature persistence under climate change, Nature, 599, 425–430, https://doi.org/10.1038/s4158602103943z, 2021. a, b, c, d, e, f, g, h
Liu, P., Zhu, Y., Zhang, Q., Gottschalck, J., Zhang, M., Melhauser, C., Li, W., Guan, H., Zhou, X., Hou, D., Peña, M., Wu, G., Liu, Y., Zhou, L., He, B., Hu, W., and Sukhdeo, R.: Climatology of tracked persistent maxima of 500hPa geopotential height, Clim. Dynam., 51, 701–717, https://doi.org/10.1007/s0038201739500, 2018. a, b
Liu, Q.: On the definition and persistence of blocking, Tellus A, 46, 286–298, https://doi.org/10.1034/j.16000870.1994.t01200004.x, 1994. a, b
Lorenz, E. N.: Can chaos and intransitivity lead to interannual variability?, Tellus A, 42, 378–389, https://doi.org/10.3402/tellusa.v42i3.11884, 1990. a
Lorenz, R., Jaeger, E. B., and Seneviratne, S. I.: Persistence of heat waves and its link to soil moisture memory, Geophys. Res. Lett., 37, L09703, https://doi.org/10.1029/2010GL042764, 2010. a, b
Lucarini, V., Faranda, D., Freitas, A. C. G. M. M., de Freitas, J. M. M., Holland, M., Kuna, T., Nicol, M., Todd, M., and Vaienti, S.: Extremes and Recurrence in Dynamical Systems, John Wiley and Sons, Inc, Hoboken, NJ, USA, https://doi.org/10.1002/9781118632321, 2016. a
MacDonald, G.: Persistence in Climate, Tech. Rep., https://irp.fas.org/agency/dod/jason/persistence.pdf (last access: 7 September 2023), 1992. a, b, c, d, e
Mailier, P. J., Stephenson, D. B., Ferro, C. A. T., and Hodges, K. I.: Serial Clustering of Extratropical Cyclones, Mon. Weather Rev., 134, 2224–2240, https://doi.org/10.1175/MWR3160.1, 2006. a, b, c, d, e, f, g, h, i, j, k
Mallakpour, I., Villarini, G., Jones, M. P., and Smith, J. A.: On the use of Cox regression to examine the temporal clustering of flooding and heavy precipitation across the central United States, Global Planet. Change, 155, 98–108, https://doi.org/10.1016/j.gloplacha.2017.07.001, 2017. a
Mallapaty, S.: Why are Pakistan’s floods so extreme this year?, https://www.nature.com/articles/d41586022028136 (last access: 7 September 2023), 2022. a
Mandelbrot, B. B. and Wallis, J. R.: Some longrun properties of geophysical records, Water Resour. Res., 5, 321–340, https://doi.org/10.1029/WR005i002p00321, 1969. a
Mann, M. E., Rahmstorf, S., Kornhuber, K., Steinman, B. A., Miller, S. K., Petri, S., and Coumou, D.: Projected changes in persistent extreme summer weather events: The role of quasiresonant amplification, Sci. Adv., 4, eaat3272, https://doi.org/10.1126/sciadv.aat3272, 2018. a
Marwan, N.: A historical review of recurrence plots, Europ. Phys. J. Spec. Top., 164, 3–12, https://doi.org/10.1140/epjst/e2008008291, 2008. a
Marwan, N., Carmen Romano, M., Thiel, M., and Kurths, J.: Recurrence plots for the analysis of complex systems, Phys. Rep., 438, 237–329, https://doi.org/10.1016/j.physrep.2006.11.001, 2007. a, b, c
Meehl, G. A., Richter, J. H., Teng, H., Capotondi, A., Cobb, K., DoblasReyes, F., Donat, M. G., England, M. H., Fyfe, J. C., Han, W., Kim, H., Kirtman, B. P., Kushnir, Y., Lovenduski, N. S., Mann, M. E., Merryfield, W. J., Nieves, V., Pegion, K., Rosenbloom, N., Sanchez, S. C., Scaife, A. A., Smith, D., Subramanian, A. C., Sun, L., Thompson, D., Ummenhofer, C. C., and Xie, S.P.: Initialized Earth System prediction from subseasonal to decadal timescales, Nat. Rev. Earth Environ., 2, 340–357, https://doi.org/10.1038/s4301702100155x, 2021. a
Meng, L., Ford, T., and Guo, Y.: Logistic regression analysis of drought persistence in East China, Int. J. Climatol., 37, 1444–1455, https://doi.org/10.1002/joc.4789, 2017. a, b
Michelangeli, P.A., Vautard, R., and Legras, B.: Weather Regimes: Recurrence and Quasi Stationarity, J. Atmos. Sci., 52, 1237–1256, https://doi.org/10.1175/15200469(1995)052<1237:WRRAQS>2.0.CO;2, 1995. a, b, c, d, e, f, g
Mo, K. C. and Ghil, M.: Statistics and Dynamics of Persistent Anomalies, J. Atmos. Sci., 44, 877–902, https://doi.org/10.1175/15200469(1987)044<0877:SADOPA>2.0.CO;2, 1987. a, b, c, d, e, f
Mohr, S., Wilhelm, J., Wandel, J., Kunz, M., Portmann, R., Punge, H. J., Schmidberger, M., Quinting, J. F., and Grams, C. M.: The role of largescale dynamics in an exceptional sequence of severe thunderstorms in Europe May–June 2018, Weather Clim. Dynam., 1, 325–348, https://doi.org/10.5194/wcd13252020, 2020. a
Moon, H., Gudmundsson, L., and Seneviratne, S. I.: Drought Persistence Errors in Global Climate Models, J. Geophys. Res.Atmos., 123, 3483–3496, https://doi.org/10.1002/2017JD027577, 2018. a, b
Mukhin, D., Hannachi, A., Braun, T., and Marwan, N.: Revealing recurrent regimes of midlatitude atmospheric variability using novel machine learning method, Chaos, 32, 113105, https://doi.org/10.1063/5.0109889, 2022. a, b, c
Mumby, P. J., Vitolo, R., and Stephenson, D. B.: Temporal clustering of tropical cyclones and its ecosystem impacts, P. Natl. Acad. Sci. USA, 108, 17626–17630, https://doi.org/10.1073/pnas.1100436108, 2011. a, b
Namias, J.: Seasonal persistence and recurrence of European blocking during 1958–1960, Tellus, 16, 394–407, https://doi.org/10.1111/j.21533490.1964.tb00176.x, 1964. a
Økland, H. and Lejenäs, H.: Blocking and persistence, Tellus A, 39, 33–38, https://doi.org/10.1111/j.16000870.1987.tb00286.x, 1987. a, b
Ontañón, S.: An overview of distance and similarity functions for structured data, Artif. Intell. Rev., 53, 5309–5351, https://doi.org/10.1007/s1046202009821w, 2020. a
Overland, J. E. and Wang, M.: The 2020 Siberian heat wave, Int. J. Climatol., 41, E2341–E2346, https://doi.org/10.1002/joc.6850, 2021. a
Pandolfo, L.: Observational Aspects of the LowFrequency Intraseasonal Variability of the Atmosphere in Middle Latitudes, Adv. Geophys., 34, 93–174, https://doi.org/10.1016/S00652687(08)604355, 1993. a
Parzen, E.: Quantile spectral analysis and longmemory time series, J. Appl. Probab., 23, 41–54, https://doi.org/10.2307/3214341, 1986. a
Pelletier, J. D. and Turcotte, D. L.: Longrange persistence in climatological and hydrological time series: analysis, modeling and application to drought hazard assessment, J. Hydrol., 203, 198–208, https://doi.org/10.1016/S00221694(97)001029, 1997. a
PerezZanon, N., et al.: CSTools: Assessing Skill of Climate Forecasts on SeasonaltoDecadal Timescales, R package version 4.1.1, https://CRAN.Rproject.org/package=CSTools (last access: 7 September 2023), 2022. a
Pfleiderer, P. and Coumou, D.: Quantification of temperature persistence over the Northern Hemisphere landarea, Clim. Dynam., 51, 627–637, https://doi.org/10.1007/s003820173945x, 2018. a, b, c, d, e
Pfleiderer, P., Schleussner, C.F., Kornhuber, K., and Coumou, D.: Summer weather becomes more persistent in a 2 degree C world, Nat. Clim. Change, 9, 666–671, https://doi.org/10.1038/s4155801905550, 2019. a, b, c
Pinto, J. G., Bellenbaum, N., Karremann, M. K., and DellaMarta, P. M.: Serial clustering of extratropical cyclones over the North Atlantic and Europe under recent and future climate conditions, J. Geophys. Res.Atmos., 118, 12476–12485, https://doi.org/10.1002/2013JD020564, 2013. a
Pinto, J. G., Gómara, I., Masato, G., Dacre, H. F., Woollings, T., and Caballero, R.: Largescale dynamics associated with clustering of extratropical cyclones affecting Western Europe, J. Geophys. Res.Atmos., 119, 13704–13719, https://doi.org/10.1002/2014JD022305, 2014. a, b, c, d
Pinto, J. G., Ulbrich, S., Economou, T., Stephenson, D. B., Karremann, M. K., and Shaffrey, L. C.: Robustness of serial clustering of extratropical cyclones to the choice of tracking method, TELLUS A, 68, 32204, https://doi.org/10.3402/tellusa.v68.32204, 2016. a
Pires, C. A. L. and Hannachi, A.: Bispectral analysis of nonlinear interaction, predictability and stochastic modelling with application to ENSO, Tellus A, 73, 1–30, https://doi.org/10.1080/16000870.2020.1866393, 2021. a, b, c
Potter, K. W.: Annual precipitation in the northeast United States: Long memory, short memory, or no memory?, Water Resour. Res., 15, 340–346, https://doi.org/10.1029/WR015i002p00340, 1979. a
Quandt, L.A., Keller, J. H., Martius, O., and Jones, S. C.: Forecast Variability of the Blocking System over Russia in Summer 2010 and Its Impact on Surface Conditions, Weather Forecast., 32, 61–82, https://doi.org/10.1175/WAFD160065.1, 2017. a
Rakovec, O., Samaniego, L., Hari, V., Markonis, Y., Moravec, V., Thober, S., Hanel, M., and Kumar, R.: The 20182020 MultiYear Drought Sets a New Benchmark in Europe, Earth's Future, 10, e2021EF002394, https://doi.org/10.1029/2021EF002394, 2022. a
RamirezAmaro, K. and FigueroaNazuno, J.: Recurrence Plot Analysis and its Application to Teleconnection Patterns, in: 2006 15th International Conference on Computing, 65–72, https://doi.org/10.1109/CIC.2006.59, 2006. a, b
Ray, R., Khondekar, M. H., Ghosh, K., and Bhattacharjee, A. K.: Complexity and periodicity of daily mean temperature and dewpoint across India, J. Earth Syst. Sci., 128, 143, https://doi.org/10.1007/s120400191174x, 2019. a
Rehman, S. and Siddiqi, A.: Waveletbased Hurst exponent and fractal dimensional analysis of Saudi climatic dynamics, Chaos Soliton Fract., 40, 1081–1090, https://doi.org/10.1016/j.chaos.2007.08.063, 2009. a
Richardson, D., Kilsby, C. G., Fowler, H. J., and Bárdossy, A.: Weekly to multimonth persistence in sets of daily weather patterns over Europe and the North Atlantic Ocean, Int. J. Climatol., 39, 2041–2056, https://doi.org/10.1002/joc.5932, 2019. a, b, c
Ripley, B. D.: Spatial Statistics, Wiley Series in Probability and Statistics, John Wiley & Sons, Inc., Hoboken, NJ, USA, https://doi.org/10.1002/0471725218, 1981. a
Robin, Y.: CDSK (Chaotic Dynamical System Kit), MIT [code], https://github.com/yrobink/CDSK (last access: 8 July 2023), 2021. a
Röthlisberger, M. and Martius, O.: Quantifying the Local Effect of Northern Hemisphere Atmospheric Blocks on the Persistence of Summer Hot and Dry Spells, Geophys. Res. Lett., 46, 10101–10111, https://doi.org/10.1029/2019GL083745, 2019. a, b, c
Röthlisberger, M., Frossard, L., Bosart, L. F., Keyser, D., and Martius, O.: Recurrent SynopticScale Rossby Wave Patterns and Their Effect on the Persistence of Cold and Hot Spells, J. Clim., 32, 3207–3226, https://doi.org/10.1175/JCLID180664.1, 2019. a, b, c, d, e
Rousi, E., Fink, A. H., Andersen, L. S., Becker, F. N., BeobideArsuaga, G., Breil, M., Cozzi, G., Heinke, J., Jach, L., Niermann, D., Petrovic, D., Richling, A., Riebold, J., Steidl, S., SuarezGutierrez, L., Tradowsky, J., Coumou, D., Düsterhus, A., Ellsäßer, F., Fragkoulidis, G., Gliksman, D., Handorf, D., Haustein, K., Kornhuber, K., Kunstmann, H., Pinto, J. G., WarrachSagi, K., and Xoplaki, E.: The extremely hot and dry 2018 summer in central and northern Europe from a multifaceted weather and climate perspective, EGUsphere [preprint], https://doi.org/10.5194/egusphere2022813, 2022a. a, b
Rousi, E., Kornhuber, K., BeobideArsuaga, G., Luo, F., and Coumou, D.: Accelerated western European heatwave trends linked to morepersistent double jets over Eurasia, Nat. Commun., 13, 3851, https://doi.org/10.1038/s4146702231432y, 2022b. a, b, c, d
Santhanam, M. S. and Kantz, H.: Return interval distribution of extreme events and longterm memory, Phys. Rev. E, 78, 051113, https://doi.org/10.1103/PhysRevE.78.051113, 2008. a
Schwierz, C., CrociMaspoli, M., and Davies, H. C.: Perspicacious indicators of atmospheric blocking, Geophys. Res. Lett., 31, L06125, https://doi.org/10.1029/2003GL019341, 2004. a
Sericola, B.: Markov Chains, John Wiley & Sons, Inc., Hoboken, NJ USA, https://doi.org/10.1002/9781118731543, 2013. a
Serinaldi, F. and Kilsby, C. G.: On the sampling distribution of Allan factor estimator for a homogeneous Poisson process and its use to test inhomogeneities at multiple scales, Physica A, 392, 1080–1089, https://doi.org/10.1016/j.physa.2012.11.015, 2013. a, b, c
Sharma, T. C. and Panu, U. S.: Modeling of hydrological drought durations and magnitudes: Experiences on Canadian streamflows, J. Hydrol. Reg. Stud., 1, 92–106, https://doi.org/10.1016/j.ejrh.2014.06.006, 2014. a
Smith, J. A. and Karr, A. F.: Flood Frequency Analysis Using the Cox Regression Model, Water Resour. Res., 22, 890–896, https://doi.org/10.1029/WR022i006p00890, 1986. a
Son, R., Wang, S.Y. S., Kim, S. H., Kim, H., Jeong, J.H., and Yoon, J.H.: Recurrent pattern of extreme fire weather in California, Environ. Res. Lett., 16, 094031, https://doi.org/10.1088/17489326/ac1f44, 2021. a
Steinfeld, D.: ConTrack – Contour Tracking of circulation anomalies in weather and climate data, Zenodo, https://doi.org/10.5281/ZENODO.4765560, 2021. a
Steinfeld, D., Boettcher, M., Forbes, R., and Pfahl, S.: The sensitivity of atmospheric blocking to upstream latent heating – numerical experiments, Weather Clim. Dynam., 1, 405–426, https://doi.org/10.5194/wcd14052020, 2020. a
Stephenson, D. B., Hannachi, A., and O'Neill, A.: On the existence of multiple climate regimes, Q. J. Roy. Meteor. Soc., 130, 583–605, https://doi.org/10.1256/qj.02.146, 2004. a
Strommen, K., Chantry, M., Dorrington, J., and Otter, N.: A topological perspective on weather regimes, Clim. Dynam., 60, 1415–1445, https://doi.org/10.1007/s0038202206395x, 2022. a
Tatli, H.: Detecting persistence of meteorological drought via the Hurst exponent, Meteorol. Appl., 22, 763–769, https://doi.org/10.1002/met.1519, 2015. a
Telesca, L.: Timeclustering of NAT HAZARDS, Nat. Hazards, 6, 593–601, https://doi.org/10.1007/s110690069023z, 2007. a, b, c
Telesca, L., Lovallo, M., and Kanevski, M.: Power spectrum and multifractal detrended fluctuation analysis of highfrequency wind measurements in mountainous regions, Appl. Energ., 162, 1052–1061, https://doi.org/10.1016/j.apenergy.2015.10.187, 2016. a
Telesca, L., Guignard, F., Laib, M., and Kanevski, M.: Analysis of temporal properties of extremes of wind measurements from 132 stations over Switzerland, Renew. Energ., 145, 1091–1103, https://doi.org/10.1016/j.renene.2019.06.089, 2020. a
Thao, S.: dtheta, MIT [code], https://github.com/thaos/dtheta (last access: 8 September 2023), 2021. a
Trenberth, K. E.: Some Effects of Finite Sample Size and Persistence on Meteorological Statistics.Part II: Potential Predictability, Mon. Weather Rev., 112, 2369–2379, https://doi.org/10.1175/15200493(1984)112<2369:SEOFSS>2.0.CO;2, 1984. a
Tuel, A.: Quantifying persistence in weather data, MIT [code], https://doi.org/10.5281/ZENODO.8329531, 2023. a
Tuel, A. and Martius, O.: A global perspective on the subseasonal clustering of precipitation extremes, Weather Clim. Extrem., 33, 100348, https://doi.org/10.1016/j.wace.2021.100348, 2021a. a, b, c, d, e, f
Tuel, A. and Martius, O.: A climatology of subseasonal temporal clustering of extreme precipitation in Switzerland and its links to extreme discharge, Nat. Hazards Earth Syst. Sci., 21, 2949–2972, https://doi.org/10.5194/nhess2129492021, 2021b. a, b
Tuel, A. and Martius, O.: Subseasonal Temporal Clustering of Extreme Precipitation in the Northern Hemisphere: Regionalization and Physical Drivers, J. Clim., 35, 3537–3555, https://doi.org/10.1175/JCLID210562.1, 2022a. a, b, c, d, e, f
Tuel, A. and Martius, O.: The influence of modes of climate variability on the subseasonal temporal clustering of extreme precipitation, Science, 25, 103855, https://doi.org/10.1016/j.isci.2022.103855, 2022b. a, b
Tuel, A. and Martius, O.: On the persistence of warm and cold spells in the Northern Hemisphere extratropics: regionalisation, synopticscale dynamics, and temperature budget, EGUsphere [preprint], https://doi.org/10.5194/egusphere20221478, 2023. a
Tuel, A., Schaefli, B., Zscheischler, J., and Martius, O.: On the links between subseasonal clustering of extreme precipitation and high discharge in Switzerland and Europe, Hydrol. Earth Syst. Sci., 26, 2649–2669, https://doi.org/10.5194/hess2626492022, 2022a. a
Tuel, A., Steinfeld, D., Ali, S. M., Sprenger, M., and Martius, O.: Largescale drivers of persistent extreme weather during early summer 2021 in Europe, Geophys. Res. Lett., 49, e2022GL099624, https://doi.org/10.1029/2022GL099624, 2022b. a, b, c
van der Walt, S., Schönberger, J. L., NunezIglesias, J., Boulogne, F., Warner, J. D., Yager, N., Gouillart, E., Yu, T., and the scikitimage contributors: scikitimage: image processing in Python, PeerJ, 2, e453, https://doi.org/10.7717/peerj.453, 2014. a
Vautard, R.: Multiple Weather Regimes over the North Atlantic: Analysis of Precursors and Successors, Mon. Weather Rev., 118, 2056–2081, https://doi.org/10.1175/15200493(1990)118<2056:MWROTN>2.0.CO;2, 1990. a, b, c, d, e, f, g
Velásquez Valle, M. A., Medina García, G., Sánchez Cohen, I., Klaudia Oleschko, L., Ruiz Corral, J. A., and Korvin, G.: Spatial Variability of the Hurst Exponent for the Daily Scale Rainfall Series in the State of Zacatecas, Mexico, J. Appl. Meteorol. Clim., 52, 2771–2780, https://doi.org/10.1175/JAMCD130136.1, 2013. a
Vigaud, N., Robertson, A., and Tippett, M. K.: Predictability of Recurrent Weather Regimes over North America during Winter from Submonthly Reforecasts, Mon. Weather Rev., 146, 2559–2577, https://doi.org/10.1175/MWRD180058.1, 2018. a
Villarini, G., Smith, J. A., Baeck, M. L., Vitolo, R., Stephenson, D. B., and Krajewski, W. F.: On the frequency of heavy rainfall for the Midwest of the United States, J. Hydrol., 400, 103–120, https://doi.org/10.1016/j.jhydrol.2011.01.027, 2011. a
Villarini, G., Smith, J. A., Vitolo, R., and Stephenson, D. B.: On the temporal clustering of US floods and its relationship to climate teleconnection patterns, Int. J. Climatol., 33, 629–640, https://doi.org/10.1002/joc.3458, 2013. a, b
Vitart, F., Ardilouze, C., Bonet, A., Brookshaw, A., Chen, M., Codorean, C., Déqué, M., Ferranti, L., Fucile, E., Fuentes, M., Hendon, H., Hodgson, J., Kang, H.S., Kumar, A., Lin, H., Liu, G., Liu, X., Malguzzi, P., Mallas, I., Manoussakis, M., Mastrangelo, D., MacLachlan, C., McLean, P., Minami, A., Mladek, R., Nakazawa, T., Najm, S., Nie, Y., Rixen, M., Robertson, A. W., Ruti, P., Sun, C., Takaya, Y., Tolstykh, M., Venuti, F., Waliser, D., Woolnough, S., Wu, T., Won, D.J., Xiao, H., Zaripov, R., and Zhang, L.: The Subseasonal to Seasonal (S2S) Prediction Project Database, B. Am. Meteorol. Soc., 98, 163–173, https://doi.org/10.1175/BAMSD160017.1, 2017. a
Vitolo, R., Stephenson, D. B., Cook, I. M., and MitchellWallace, K.: Serial clustering of intense European storms, Meteorol. Z., 18, 411–424, https://doi.org/10.1127/09412948/2009/0393, 2009. a, b
von Lindheim, J., Harikrishnan, A., Dörffel, T., Klein, R., Koltai, P., Mikula, N., Müller, A., Névir, P., Pacey, G., Polzin, R., and Vercauteren, N.: Definition, detection, and tracking of persistent structures in atmospheric flows, https://doi.org/10.48550/ARXIV.2111.13645, 2021. a
Weatherhead, E., Gearheard, S., and Barry, R.: Changes in weather persistence: Insight from Inuit knowledge, Glob. Environ. Change, 20, 523–528, https://doi.org/10.1016/j.gloenvcha.2010.02.002, 2010. a
Weber, J., Reyers, M., Beck, C., Timme, M., Pinto, J. G., Witthaut, D., and Schäfer, B.: Wind Power Persistence Characterized by Superstatistics, Sci. Rep.UK, 9, 19971, https://doi.org/10.1038/s41598019562861, 2019. a
Weiland, R. S., van der Wiel, K., Selten, F., and Coumou, D.: Intransitive Atmosphere Dynamics Leading to Persistent HotDry or ColdWet European Summers, J. Clim., 34, 6303–6317, https://doi.org/10.1175/JCLID200943.1, 2021. a, b, c
Weiss, J. P. and Weiss, J. B.: Quantifying Persistence in ENSO, J. Atmos. Sci., 56, 2737–2760, https://doi.org/10.1175/15200469(1999)056<2737:QPIE>2.0.CO;2, 1999. a, b
Wernli, H. and Schwierz, C.: Surface Cyclones in the ERA40 Dataset (1958–2001). Part I: Novel Identification Method and Global Climatology, J. Atmos. Sci., 63, 2486–2507, https://doi.org/10.1175/JAS3766.1, 2006. a
Wharton, E., Panetta, K., and Agaian, S.: Human visual system based similarity metrics, in: 2008 IEEE International Conference on Systems, Man Cybernet., 685–690, https://doi.org/10.1109/ICSMC.2008.4811357, 2008. a
Wikipedia: 2022 Eastern Australia floods, https://en.wikipedia.org/wiki/2022_Eastern_Australia_floods (last access: 7 September 2023), 2022. a
Wilby, R. L., Noone, S., Murphy, C., Matthews, T., Harrigan, S., and Broderick, C.: An evaluation of persistent meteorological drought using a homogeneous Island of Ireland precipitation network, Int. J. Climatol., 36, 2854–2865, https://doi.org/10.1002/joc.4523, 2016. a
Witt, A. and Malamud, B. D.: Quantification of LongRange Persistence in Geophysical Time Series: Conventional and BenchmarkBased Improvement Techniques, Surv. Geophys., 34, 541–651, https://doi.org/10.1007/s1071201292178, 2013. a
Wolff, N. H., Wong, A., Vitolo, R., Stolberg, K., Anthony, K. R. N., and Mumby, P. J.: Temporal clustering of tropical cyclones on the Great Barrier Reef and its ecological importance, Coral Reefs, 35, 613–623, https://doi.org/10.1007/s0033801614009, 2016. a
Wolters, M. A.: Better Autologistic Regression, Front. Appl. Mathemat. Stat., 3, 24, https://doi.org/10.3389/fams.2017.00024, 2017. a
Woollings, T., Hannachi, A., and Hoskins, B.: Variability of the North Atlantic eddydriven jet stream, Q. J. Roy. Meteor. Soc., 136, 856–868, https://doi.org/10.1002/qj.625, 2010. a, b
World Meteorological Organization: Extreme weather in China highlights climate change impacts and need for early warnings, https://public.wmo.int/en/ (last access: 7 September 2023), 2022. a
Yang, L. and Fu, Z.: Processdependent persistence in precipitation records, Physica A, 527, 121459, https://doi.org/10.1016/j.physa.2019.121459, 2019. a, b
Yang, Z. and Villarini, G.: Examining the capability of reanalyses in capturing the temporal clustering of heavy precipitation across Europe, Clim. Dynam., 53, 1845–1857, https://doi.org/10.1007/s0038201904742z, 2019. a, b
Yiou, P., Cattiaux, J., Ribes, A., Vautard, R., and Vrac, M.: Recent Trends in the Recurrence of North Atlantic Atmospheric Circulation Patterns, Complexity, 2018, 3140915, https://doi.org/10.1155/2018/3140915, 2018. a
Yuan, N., Fu, Z., and Mao, J.: Different scaling behaviors in daily temperature records over China, Physica A, 389, 4087–4095, https://doi.org/10.1016/j.physa.2010.05.026, 2010. a
Zerzucha, P. and Walczak, B.: Concept of (dis)similarity in data analysis, Trac.Trend. Anal. Chem., 38, 116–128, https://doi.org/10.1016/j.trac.2012.05.005, 2012. a
Zolina, O., Simmer, C., Belyaev, K., Gulev, S. K., and Koltermann, P.: Changes in the Duration of European Wet and Dry Spells during the Last 60 Years, J. Clim., 26, 2022–2047, https://doi.org/10.1175/JCLID1100498.1, 2013. a, b
 Abstract
 Introduction
 Persistence: quasistationarity or recurrence?
 Considerations and methodological approaches for persistence analysis in weather and climate data
 Quasistationarity
 Recurrence
 Summary and outlook
 Code availability
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Persistence: quasistationarity or recurrence?
 Considerations and methodological approaches for persistence analysis in weather and climate data
 Quasistationarity
 Recurrence
 Summary and outlook
 Code availability
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Financial support
 Review statement
 References