1. Introduction
An unprecedented heatwave across British Columbia (BC), Canada, and the Pacific Northwest (PNW) region of the United States in June–July 2021 broke many daily maximum surface temperature records (Philip et al. 2021). This heatwave was dominated by an atmospheric blocking event associated with a strong upper-level ridge and anomalously high geopotential heights that persisted for longer than a week. Coincident with the event, hundreds of fires started across the region, prompting local and regional emergency declarations and the evacuation of thousands of people (BC Gov News 2021).
Persistent positive anomalies (PPAs) in midtropospheric geopotential heights (e.g., 500-hPa geopotential height) characterize atmospheric blocking highs and ridges (Elliott and Smith 1949; Dole and Gordon 1983). Synoptically, the PPAs represent an obstruction of the flow of planetary waves and include the blocking high, open ridges, and occasionally the poleward intrusion of subtropical ridges (Elliott and Smith 1949; Sousa et al. 2021). Although the term PPA refers to a variety of atmospheric blocking patterns, in general, it is well established that atmospheric blocking ridges have common attributes, including a persistent large anticyclonic anomaly and weakened or reversed zonal flow (Tibaldi and Molteni 1990; Liu 1994; Schwierz et al. 2004; Renwick 2005; Barriopedro et al. 2006; Woollings et al. 2018).
PPAs dynamically affect midlatitude weather phenomena through the inertia of zonal flow and atmospheric subsidence that can lead to extreme heatwaves (Horton et al. 2015; Mann et al. 2018; Lupo et al. 2019; Steinfeld and Pfahl 2019; Francis et al. 2020). This occurs through land–atmosphere feedbacks as evident in the mega-heatwaves of 2003 and 2010 in Europe (Dole et al. 2011; Miralles et al. 2014) and summertime heatwaves over North America (Lau and Nath 2012). Given the robust influence of PPAs on midlatitude surface weather, there exist ample studies that investigate the high pressure blocking ridges and their effects on surface variables such as temperature extremes (Woollings et al. 2014; Whan et al. 2016), precipitation (Antokhina et al. 2016), and heatwaves (Rodrigues and Woolings 2017; Schaller et al. 2018).
Synoptic-scale circulation patterns can similarly determine whether conditions are conducive for wildfire activity through their direct or indirect influence on fuel availability and moisture, receptivity to ignitions, and potential rate of fire spread. As such, several studies have identified a strong relationship between the occurrence of large wildfire events and the stagnation of upper airflow (Flannigan and Harrington 1988; Johnson and Wowchuk 1993; Skinner et al. 2002; Gedalof et al. 2005; Fauria and Johnson 2008; Rasilla et al. 2010; Petoukhov et al. 2018; Zhao and Liu 2019; Tan et al. 2019; Zhong et al. 2020; Abatzoglou et al. 2021b; Yasunari et al. 2021; Jain and Flannigan 2021). For example, Jain and Flannigan (2021) observed a close relationship between upper-air ridges and fire spread days across North America. Likewise, extreme fire events such as the 2016 Horse River fire coincided with a large amplitude ridge (Tan et al. 2019). These examples illustrate the role of PPA in specific fire events or fire behavior such as fire spread days; however, there exist gaps in better understanding the climatological relation between PPA and fire-weather conditions, fire ignition, and spread across the broad geography of western North America. Additionally, we do not understand how PPAs that favor fire weather conditions are changing and what mechanisms drive these changes.
While some studies report trends in atmospheric blocking ridges (Davini et al. 2012; Barnes et al. 2014; Gibson et al. 2020; Wazneh et al. 2021), there exist inconsistencies in these results, likely due to differences in working definitions, datasets, time, location, and season considered. Similarly, the dynamical processes resulting in blocking ridges over the Northern Hemisphere are still an active research area (Rodrigues and Woolings 2017; Riboldi et al. 2020). Some studies hypothesize the linkages between Arctic amplification (AA)—an enhanced rise in surface temperature in the Arctic compared to the remainder of the Northern Hemisphere (Francis and Vavrus 2012)—and increased weather extremes due to dynamical changes in atmospheric circulation (Francis and Vavrus 2015; Francis et al. 2018). AA has been suggested to increase the sinuosity of the Northern Hemisphere jet stream in the upper troposphere, likely due to weakening of the poleward temperature gradient (Francis and Vavrus 2015) that may favor persistent weather patterns and enhance extreme weather events in both summer and winter (Francis et al. 2018). However, other models and reanalysis-based studies have not shown strong relationships between AA and either waviness in upper airflow or midlatitude surface weather extremes (Screen 2014; Blackport and Screen 2020; Dai and Song 2020).
In the context of the aforementioned literature and gaps, the overall goal of this work is to investigate the link between PPAs and wildfire events and explore changes in PPA over 1979–2020 in western North America. Specifically, through this work, we address the following questions: 1) What is the climatology of PPAs during the main western North American fire season? 2) How do PPA events relate to wildfire activity and surface fire weather? 3) What are the trends of various PPA characteristics and what is the cause of any observed trends?
2. Data and methods
a. Study area
Our study domain is western North America covering 35°–75°N and 100°–160°W, an area that includes the western part of Canada and the contiguous United States, as well as Alaska (Fig. 1). One of the most fire-active areas globally, this region exhibits increasing trends in wildfire number and area burned (Schoennagel et al. 2017), which are expected to continue in the future, under warmer climate scenarios (Gillett et al. 2004; Flannigan et al. 2005; Balshi et al. 2009; Balch et al. 2017; Abatzoglou et al. 2021a; Kharuk et al. 2021). The 2016 Horse River Wildfire (in Fort McMurray, Alberta) (Tymstra et al. 2020), the September 2020 large fires in western Oregon (Abatzoglou et al. 2021b), the Camp Fire in California (Brewer and Clements 2020), and the 2017 record-breaking BC fire season (Kirchmeier-Young et al. 2019) are a few examples of the costliest and deadliest wildfire disasters that have occurred in this region in the past decade. The recent rise in large wildfires in this region is partially attributed to the warmer and drier climates associated with increased surface air temperature, earlier snowmelt, and a surge in lightning activity, especially in the northern latitudes (Abatzoglou and Williams 2016; Holzworth et al. 2021).

Map of the study domain. The blue rectangle indicates the area used for the identification of the PPA events. The map also shows the topography and indicates names of places mentioned in the main text.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Map of the study domain. The blue rectangle indicates the area used for the identification of the PPA events. The map also shows the topography and indicates names of places mentioned in the main text.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Map of the study domain. The blue rectangle indicates the area used for the identification of the PPA events. The map also shows the topography and indicates names of places mentioned in the main text.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
b. Data
We use the 500-hPa geopotential heights [Z500 (gpm)] and surface variables from the European Centre for Medium-Range Weather Forecasts (ERA5) global reanalysis data (Hersbach et al. 2020) for the extended summer season (May–September) of 1979–2020. We limited the study to the boreal summer months as the vast majority of wildfire events (92% of total fires with burned area >100 ha) within our study domain occur during these months (Fig. 2) despite recent increases in wildfire events outside of these months (Balch et al. 2017).

Average percentage of monthly fires (area burned > 100 ha) within the study domain.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Average percentage of monthly fires (area burned > 100 ha) within the study domain.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Average percentage of monthly fires (area burned > 100 ha) within the study domain.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Although the ERA5 data are available at higher spatial (0.25° × 0.25°) and temporal (hourly) resolution, we aggregated the data to 1° × 1° spatial resolution and to daily averages because the PPAs are reasonably well represented at this coarser resolution (Barriopedro et al. 2006; Barnes et al. 2012) and this provides greater computational efficiency. Fire-related surface variables used in this analysis included the 2-m air temperature, precipitation, 2-m dewpoint temperature, and 10-m meridional and zonal wind components extracted for the study domain. We calculate relative humidity (RH) and vapor pressure deficit (VPD) from daily air and dewpoint temperature using vapor pressures following Tsonis (2013).
In addition, we use the Canadian Fire Weather Index System (Van Wagner 1987) fire weather indices that were calculated using ERA5 surface data (McElhinny et al. 2020). This system accounts for the effects of interseasonal drought and reflects the effect of weather on fuel moisture and potential fire behavior. The fire weather indices used include the fine fuel moisture code (FFMC), Duff moisture code (DMC), drought code (DC), initial spread index (ISI), and the fire weather index (FWI). The FFMC and DMC represent moisture conditions for shaded litter fuels and decomposed organic matter underneath litter with ∼16-h and 2-week equilibrium drying times, respectively. The DC represents deeper soil moisture with about a 53-day equilibrium drying time, and it is a useful indicator of seasonal effects of weather on fuel moisture. The ISI depends on FFMC and wind, and it represents the potential fire spread rate. Finally, FWI is a general index of fire danger and represents a numeric rating of potential fire intensity (Van Wagner 1987; McElhinny et al. 2020).
We used observed fire data from two data products, the Global Fire Atlas (GFA) and the fourth-generation global fire emission database (GFED4.1), as they provide different fire-related variables of interest. The GFA provides fire ignition timing and location, fire duration, size (minimum ∼21 ha), and perimeter of individual wildfires for 2003–18 (Andela et al. 2019), whereas GFED4.1 provides the fraction of monthly burned area at 0.25° × 0.25° spatial resolution from 1997 to 2016 (Giglio et al. 2013). The GFA data are based on the Moderate Resolution Imaging Spectroradiometer (MODIS) burned-area dataset at 500-m resolution (2019). Giglio et al. (2013) derived the GFED4.1 data through a combination of 500-m MODIS burned area and active fire data from the Visible Infrared Imaging Radiometer Suite (VIIRS). The GFA and GFED4.1 fire parameters are in good agreement with observed independent fire parameters, especially in forest and shrublands (Giglio et al. 2013; Andela et al. 2019), and provide a consistent fire database for our study domain.
3. Methods
a. Identification of PPA
Given the class of synoptic patterns that blocking ridges represent, various definitions and identification procedures exist, perhaps owing partly to the research implications of the blocking events of choice (Woollings et al. 2018). Broadly, midlatitude atmospheric blocking ridges are identified based on three approaches: a determination of the persistent positive departure above time-averaged geopotential height field (Dole and Gordon 1983; Renwick 2005; Gibson et al. 2020; Miller et al. 2020), use of a dynamical index based on integrated potential vorticity (Pelly and Hoskins 2003; Schwierz et al. 2004; Small et al. 2014), and measurement of the reversal of meridional flow (Tibaldi and Molteni 1990; Pinheiro et al. 2019; Wazneh et al. 2021). The dynamic potential vorticity–based approach is tied to strong polar dynamics that are not common during summer months within the study domain (Schwierz et al. 2004). The traditional blocking climatology based on measurement of reversal of meridional flow shows no or minimal blocking episodes over our study domain during the spring and/or summer months (e.g., Fig. S1 in the online supplemental material) consistent with other studies (Barriopedro et al. 2006; Drouard and Woollings 2018; Pinheiro et al. 2019; Wazneh et al. 2021). We therefore follow the first approach to identify PPAs that highlights the anomalous two-dimensional area of the blocking and/or ridges and it is similar to previous studies such as Dole and Gordon (1983), Renwick (2005), Parsons et al. (2016), and Miller et al. (2020). However, we modify this approach to facilitate the identification of PPA events during boreal summer by considering the persistence of anomalies relative to the climatological mean of geopotential heights.
We develop and implement a feature tracking algorithm that follows the spatial and temporal extent of PPA grid cells within the study domain that satisfies the above criteria and allows us to label individual PPA events. Specifically, we identify Z′ grid cells that are of at least M magnitude, label them as potential PPAs, and only retain the resulting grid cells that persist for at least D days. Then, the algorithm identifies the geometric centroid of each cluster of contiguous PPA grid cells for each day. It tracks the movement of the centroid of the PPA until the contiguous PPA area is greater than ∼80 000 km2. We choose this minimum area threshold to ensure that the small events are not included. When multiple PPAs are identified within the study domain in a single day, individual clusters based on their centroid are tracked as individual PPA events. To include all PPA events that occurred within May–September, we allow backward (forward) tracking of the PPA events if they start (end) on 1 May (30 September). For example, Fig. S3 shows the evolution of one of the PPA events identified following this approach. The PPA identified here are sensitive to the threshold used to define them. A sensitivity analysis with more relaxed M and D criteria (M = 0.5 × σ and D = 3 days) incorporates weaker PPA events, thus increasing PPA numbers by ∼25%. However, it does not show a substantial difference in PPA relationships to fire weather and trends in PPA characteristics.
b. Statistical analysis
We implement multiple descriptive statistics to analyze PPA–wildfire associations and PPA event characteristics. Several definitions were considered: 1) The percentage of days with PPA (“PPA %-days”) were defined as the fraction of days (for May–September) that satisfy both M and D criteria; 2) the frequency of PPA was defined by counts of individual (spatially contiguous) PPA events; 3) PPA duration is the number of days for which each PPA event is detected; 4) the spatial extent of PPA grid cells for the duration of each PPA event provides the average area of that PPA event; and 5) the PPA strength is the product of PPA magnitude, size, and duration.
To calculate the percentage of wildfires during each PPA event, we determine whether the fire perimeter of individual large fire events (burned area > 500 ha) overlaps with at least one PPA grid cell within the fire duration or up to three days prior to the fire ignition. If there is an overlap, we assign the fire event as a fire associated with a PPA event. Thus, the number of fires related to PPAs out of the total number of fires provides the percentage of PPA-related fire events. Similarly, the sum of PPA days that overlap with a fire duration and have at least one common grid cell between fire perimeter and PPA grid cells provides the percentage of active fire days linked to PPAs.
To further understand the synoptic conditions associated with PPA days having at least one fire ignition, we used self-organizing maps (SOMs) to determine the dominant observed patterns. The SOM with four nodes (2 × 2 grids) provides visualization of the circulation patterns of PPA days with at least one fire ignition within PPA extent. SOMs are widely used unsupervised artificial neural networks well adapted for pattern recognition and classification; details on SOMs can be found in Kohonen (2001). The composite of circulation patterns for each node discerns the synoptic pattern associated with PPA days. In addition, we quantify the percentage of PPA days associated with each SOM node.
The interannual Spearman’s rank correlation coefficient (ρ) correlates PPA %-days and the fraction of averaged burned area over 1997–2016. Similarly, it relates the monthly PPA %-days and monthly fire weather indices over 1979–2020; ρ assesses the monotonic relation between the variables without the assumption of normality. The differences between daily values and climatological means provide daily anomalies of fire-related variables and fire weather indices. The nonparametric Mann–Whitney–Wilcoxon (MWW) rank-sum test (Wilks 2011) shows the difference in the median of the distribution of fire-related surface variables during PPA days and non-PPA days. The MWW test is used here to determine whether the median of the anomalies in surface variables are significantly greater or less during PPA days than during non-PPA days. Moreover, we explore the temporal evolution of the fire-related surface variables before, during, and after the PPA events. The lead–lag relationship between maximum PPA strength and daily surface variable anomalies presents the surface variable response to PPA events. We calculate the daily PPA strength from 15 days before to 15 days after the day with maximum PPA strength. We mask the grid cells of fire-related variables using the PPA area calculated on the day with maximum PPA strength. Finally, we calculate daily lead–lag average values of the fire-related variables for the masked area.
We perform the Mann–Kendall test (Mann 1945; Kendall 1975) and the Theil–Sen (Sen 1968) trend estimate over 1979–2020, which yield monotonic trends and trend magnitude of the PPA characteristics (Sen 1968). This nonparametric approach is robust to outliers and does not require data normality. We use the R programming language (R Core Team 2019) and the packages raster (Hijmans and Etten 2012), igraph (Csardi and Nepusz 2006), zyp (Bronaugh and Werner 2019), kohonen (Wehrens and Kruisselbrink 2018), and tidyverse (Wickham et al. 2019) to perform all analysis.
4. Results
a. PPA climatology and fire weather
On average, 17% (26 days) of May–September are PPA days with the greatest frequency occurring in the northeastern Pacific (Figs. 1 and 3a). PPA days are relatively frequent west of the Rocky Mountains compared to the prairies and boreal plains of central North America. Higher values of PPA days exist over the northeast Pacific Ocean and Beaufort Sea expanding inland to Alaska, BC, and the PNW region of the United States and the boreal plains respectively (Figs. 1 and 3a). The spatially averaged monthly percentage of PPA days is 19% in May, 17% in June, July, and August, and 16% in September. Moreover, the Hovmöller diagram of average PPA days’ Z500 anomalies within the study domain indicates higher PPA amplitudes during the early and late periods in the fire season relative to July and August (Fig. S4).

(a)–(f) Average percentage of PPA days, 1979–2020. The blue rectangle in (a) shows the domain considered for the identification and tracking of the PPA events.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

(a)–(f) Average percentage of PPA days, 1979–2020. The blue rectangle in (a) shows the domain considered for the identification and tracking of the PPA events.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
(a)–(f) Average percentage of PPA days, 1979–2020. The blue rectangle in (a) shows the domain considered for the identification and tracking of the PPA events.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
We identify a total of 517 PPA events during 1979–2020 (May–September) within our study domain, with an average of 12 (SD of 3) PPA events per year. The average duration of a PPA event is 12 days, with 24% and 4% of events having lengths of ≥15 and ≥31 days, respectively (Fig. 4). The average magnitude and area of the PPA events are 99 m and 2.2 × 106 km2, respectively.

Distribution of the duration of PPA events, 1979–2020.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Distribution of the duration of PPA events, 1979–2020.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Distribution of the duration of PPA events, 1979–2020.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
On average, 40% of fires with burned area > 500 ha were ignited within PPAs during 2003–18; this association increased to 60% for fires with area burned > 20 000 ha. Similarly, 29% of fire days with burned area > 500 ha were also PPA days and that number increased to 38% for fires with burned area > 20 000 ha.
The SOM classification of the circulation patterns associated with PPA-days having at least one fire ignition indicates that the majority of fire ignitions co-occurred with a 500-hPa open-wave ridge over western North America (Fig. 5). Among the four nodes, nodes 1 and 3, with a small difference in the location and PPA amplitude represent 47% and 23% of fire-related PPA days where a strong ridge and positive anomalies build over northern high-latitude regions including Alaska, northern British Columbia, Yukon, and Northwest Territories (Fig. 5). Nodes 2 and 4 represent the circulation for 16% and 14% of PPA days with at least one fire ignition. These two nodes exhibit the strongest PPA amplitude and display Rex- and omega-like blocking patterns, respectively. Moreover, a monthly breakdown of SOM analysis indicates small differences in the percentage of PPA-day circulation among the four nodes (Fig. S5). The difference in the percentage of PPA-day circulation patterns among all four nodes is lower (∼10%) early in the fire season, and increases to ∼17% later in the fire season (Fig. S5).

The 2 × 2 SOM patterns of 500-hPa geopotential height and anomalies (gpm) corresponding to fire ignition PPA days. The number at the top left of each panel indicates the percentage of the days of the occurrence of each pattern.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

The 2 × 2 SOM patterns of 500-hPa geopotential height and anomalies (gpm) corresponding to fire ignition PPA days. The number at the top left of each panel indicates the percentage of the days of the occurrence of each pattern.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
The 2 × 2 SOM patterns of 500-hPa geopotential height and anomalies (gpm) corresponding to fire ignition PPA days. The number at the top left of each panel indicates the percentage of the days of the occurrence of each pattern.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
b. PPA and wildfire relation
The OR is >1 for 89% of grid cells (49% significant) within the study area (Fig. 6a). The interpretation of OR > 1 is that wildfire ignitions are more likely to occur during or following (up to three days) a PPA event. On average, the odds of wildfire ignition concurrent with a PPA increases by a factor of 7 across the study domain; this number is even greater for higher latitude and altitude regions (Figs. 1 and 6a). For example, the odds of fire ignition with the presence of PPA increases by a factor of 8 above 50°N latitude whereas it is only a factor of 5 below 50°N.

(a) The OR of occurrence of wildfire ignition during the presence of PPA, 2003–18. The white grid cells indicate insufficient data for OR calculation. (b) Spearman’s correlation coefficients between the average May–September burned area and percentage of PPA days, 1997–2016. The black dots indicate grid cells with a significant value (p < 0.1). The density plots on the left side indicate the distribution of mean values along the latitudes.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

(a) The OR of occurrence of wildfire ignition during the presence of PPA, 2003–18. The white grid cells indicate insufficient data for OR calculation. (b) Spearman’s correlation coefficients between the average May–September burned area and percentage of PPA days, 1997–2016. The black dots indicate grid cells with a significant value (p < 0.1). The density plots on the left side indicate the distribution of mean values along the latitudes.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
(a) The OR of occurrence of wildfire ignition during the presence of PPA, 2003–18. The white grid cells indicate insufficient data for OR calculation. (b) Spearman’s correlation coefficients between the average May–September burned area and percentage of PPA days, 1997–2016. The black dots indicate grid cells with a significant value (p < 0.1). The density plots on the left side indicate the distribution of mean values along the latitudes.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Most of the study domain exhibits a positive correlation between the average May–September burned area and PPA %-days during 1997–2016 (Fig. 6b). Higher-latitude regions, such as Alaska, Yukon, and Northwest Territories, experience a greater interannual positive correlation between the burned area and PPA %-days (Figs. 1 and 6). This is the same region that exhibits formation of a strong 500-hPa open-wave ridge during PPA days with fire ignition (Fig. 5). Similarly, there exists a significant positive interannual correlation between May–September PPA %-days and fire weather indices (Fig. 7).

The Spearman’s correlation coefficients between May–September average fire weather indices and PPA-% days, 1979–2020. The black dots indicate grid cells with a significant correlation ( p < 0.1). The density plots on the right side show the distribution of mean positive correlation coefficients as a function of latitude.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

The Spearman’s correlation coefficients between May–September average fire weather indices and PPA-% days, 1979–2020. The black dots indicate grid cells with a significant correlation ( p < 0.1). The density plots on the right side show the distribution of mean positive correlation coefficients as a function of latitude.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
The Spearman’s correlation coefficients between May–September average fire weather indices and PPA-% days, 1979–2020. The black dots indicate grid cells with a significant correlation ( p < 0.1). The density plots on the right side show the distribution of mean positive correlation coefficients as a function of latitude.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
The surface air temperature, VPD, and fire weather indices (FWI, FFMC, DMC, DC, ISI) show larger positive anomalies during PPA days than during non-PPA days, while the atmospheric moisture–related variables (daily RH and precipitation) show lower anomalies (Fig. 8). Spatially, the magnitudes of PPA-day surface variable anomalies are greater inland except for precipitation and wind speed. The average wind speed anomalies during PPA days are significantly lower than non-PPA days wind speed (Table S1); however, few grid cells along the coast and farther inland exhibit positive wind speed anomalies (Fig. 8i). Further, the MWW test results indicate that the distribution of mean air temperature, VPD, and FWI output anomalies during PPA days are significantly higher than during non-PPA days (Table S1). Similarly, the average RH and precipitation are significantly lower during PPA days than non-PPA days (Table S1).

Fire-related surface variables (temperature, precipitation, RH, VPD, and wind) and fire weather indices (FWI, FFMC, DMC, DC, and ISI) anomalies and the distribution of anomalies averaged over 1979–2020 during PPA days. The box-and-whisker plots show the distribution of anomalies during PPA days. The blue horizontal lines indicate median anomalies; red dots show mean anomalies; vertical lines present the range of PPA and non-PPA day anomalies within the 1.5 interquartile range. The black dots indicate outliers of the anomalies.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Fire-related surface variables (temperature, precipitation, RH, VPD, and wind) and fire weather indices (FWI, FFMC, DMC, DC, and ISI) anomalies and the distribution of anomalies averaged over 1979–2020 during PPA days. The box-and-whisker plots show the distribution of anomalies during PPA days. The blue horizontal lines indicate median anomalies; red dots show mean anomalies; vertical lines present the range of PPA and non-PPA day anomalies within the 1.5 interquartile range. The black dots indicate outliers of the anomalies.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Fire-related surface variables (temperature, precipitation, RH, VPD, and wind) and fire weather indices (FWI, FFMC, DMC, DC, and ISI) anomalies and the distribution of anomalies averaged over 1979–2020 during PPA days. The box-and-whisker plots show the distribution of anomalies during PPA days. The blue horizontal lines indicate median anomalies; red dots show mean anomalies; vertical lines present the range of PPA and non-PPA day anomalies within the 1.5 interquartile range. The black dots indicate outliers of the anomalies.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
The lead–lag relation (15 days before and after from the day of maximum PPA strength) between PPA strength and surface variables indicate that the temperature and VPD anomalies increase with PPA growth while the precipitation and RH anomalies decrease (Fig. 9). The peak increase (decrease) in temperature and VPD (RH) anomalies coincides with maximum PPA strength, whereas the precipitation deficit peaks two days earlier than the PPA peak (Fig. 9). For all these variables, the anomalies weaken with the weakening of PPA. The fire weather indices increase with PPA strength and decrease with its decay except for DC that has persistent high values even after the weakening of PPA events consistent with the continued drying of deeper soils after the maximum PPA strength occurs (Fig. 9).

Lead–lag relation between PPA strength and fire-related surface variables and fire weather indices. Here 0-day represents the day of maximum PPA strength for each event; the red lines indicate PPA strength with maximum strength on 0-day. The blue lines indicate averaged PPA day anomalies of each parameter.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Lead–lag relation between PPA strength and fire-related surface variables and fire weather indices. Here 0-day represents the day of maximum PPA strength for each event; the red lines indicate PPA strength with maximum strength on 0-day. The blue lines indicate averaged PPA day anomalies of each parameter.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Lead–lag relation between PPA strength and fire-related surface variables and fire weather indices. Here 0-day represents the day of maximum PPA strength for each event; the red lines indicate PPA strength with maximum strength on 0-day. The blue lines indicate averaged PPA day anomalies of each parameter.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
c. PPA trends
Having examined the relationship between PPA occurrence and wildfire ignitions, we examine trends in PPA characteristics from 1979 to 2020. Increases in the percentage of days during May–September with PPAs were observed for most of the study domain, with significant increases over the Pacific and Alaska (Fig. 10). During 1979–2020 May–September PPA days increased by eight days on average over our study area. A significant rise in PPA days occurred over Alaska, BC, and the PNW region, whereas no significant increase was observed farther inland within our study domain (Fig. 10).

The trend in May–September PPA-%days, 1979–2020. The black dots indicate grid cells with a significant trend magnitude (p < 0.1).
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

The trend in May–September PPA-%days, 1979–2020. The black dots indicate grid cells with a significant trend magnitude (p < 0.1).
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
The trend in May–September PPA-%days, 1979–2020. The black dots indicate grid cells with a significant trend magnitude (p < 0.1).
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
The average number of PPA events shows a small statistically nonsignificant rise (1.6 events, p = 0.12) over the period of 1979–2020 (Fig. 11a). PPA events’ average duration and magnitude exhibit small (not statistically significant) increasing trends, while the average area of the PPA events shows a statistically significant expansion of 1.06 × 104 km2 yr−1 (p = 0.08) over 1979–2020, representing a 20% increase over the period of record (Fig. 11c).

Trends in average PPA (a) number, (b) duration, (c) area, and (d) magnitude, 1979–2020. The red lines represent trends from the Theil–Sen estimator while the inset on each plot gives the trend magnitude and corresponding p value; p < 0.1 indicates statistically significant trends.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Trends in average PPA (a) number, (b) duration, (c) area, and (d) magnitude, 1979–2020. The red lines represent trends from the Theil–Sen estimator while the inset on each plot gives the trend magnitude and corresponding p value; p < 0.1 indicates statistically significant trends.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Trends in average PPA (a) number, (b) duration, (c) area, and (d) magnitude, 1979–2020. The red lines represent trends from the Theil–Sen estimator while the inset on each plot gives the trend magnitude and corresponding p value; p < 0.1 indicates statistically significant trends.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
The breakdown of PPA characteristics into different latitudinal bands reveals variations in the trends by latitude despite the latitude-dependent normalization of Z500. The number of PPA events increases significantly only along the lower latitudinal band (30°–45°N; 0.04 events yr−1) while a significant increase in the PPA duration (0.16 days yr−1) occurs along the highest latitude band (60°–75°N) (Fig. S6). However, the PPA area increases along all latitudinal bands, with a significant increase along the higher latitudes ( p = 0.05) (Fig. S6).
In addition, the monthly breakdown of PPAs indicates substantial growth in the PPA area early in the fire season (May–July) (Fig. S7); other PPA characteristics, namely number, magnitude, and duration, do not exhibit significant monthly changes (not shown). In addition, we explore trends on only those PPA events that intersect with the land area of western North America as they directly influence surface fire weather conditions. They exhibit trends similar to those of all PPA events (Fig. S8).
5. Discussion
a. Comparison with previous studies
Differing approaches for defining PPA events (Barnes et al. 2012; Woollings et al. 2018) impede direct comparisons with previous studies. The climatology of summertime PPAs observed here is spatially consistent with previous studies based on Z500 anomaly measurement algorithms (Dunn-Sigouin et al. 2013; Miller et al. 2020) or dynamic potential vorticity–based approaches (Pelly and Hoskins 2003; Small et al. 2014); however, there is a small difference in magnitude compared to our approach. For example, we find the June, July, and August average PPA %-days as 17% (16 days) over the same region as our study domain. This is higher than Dunn-Sigouin et al. (2013) and Miller et al. (2020), who report ∼10%–12% days and <5 days summer PPA %-days with M of 1.5 × σ and 2 × σ, respectively. While previous studies have used higher thresholds, our approach is less conservative in that it includes lower-magnitude persistent events, especially during summer. Similarly, Small et al. (2014), using a dynamic potential vorticity–based approach, observed ∼8%–10% of days in summer months blocked over the Gulf of Alaska and western North America whereas we find 17% of days blocked for the same months (Fig. 3).
b. PPA and wildfires
Previous studies have established links between summertime blocking ridges, heatwaves (Gedalof et al. 2005; Dole et al. 2011; Miralles et al. 2014; Mann et al. 2018; Schaller et al. 2018), and wildfires (Flannigan et al. 2005; Zhong et al. 2020; Jain and Flannigan 2021; Yasunari et al. 2021) but mostly for specific fire events and/or small regions. This study extends previous efforts relating PPAs and fire weather, fire ignition, and burned area and quantifies the role of PPAs on fire ignitions across the broad geography of western North America. There exists a substantial association between PPA and fire ignitions and fire burned area; this relation is greater in western North America’s higher-latitude regions than in lower latitudes (Fig. 6). In addition, the significant increase in PPA days, PPA area, and duration that occurs in the higher-latitude band (60°–75°N) in recent decades (Fig. S6) may result in enhanced surface fire weather in this region. This, in turn, leads to disproportionately higher PPA–fire associations over northern latitudes where greater surface weather anomalies are required to reach extreme fire weather conditions.
The presence of a PPA amplifies fire weather through cascading effects on surface weather conditions as a result of land–atmosphere interactions. These interactions result in enhanced air subsidence, warm air advection, and moisture diversion that leads to higher temperatures, clear skies, lower relative humidity, and low evaporative fraction (Horton et al. 2016). As a result, a positive feedback loop is created with anomalously high surface temperatures. We observe strikingly high positive anomalies in temperature during PPAs (Fig. 8a) that acts as a major predictor of wildfire and area burned (Flannigan and Harrington 1988; Gillett et al. 2004). Similarly, the presence of a PPA leads to the forced diversion of jet streams and associated baroclinic disturbances leading to negative precipitation anomalies within PPA-influenced areas, raising air temperature and reducing RH (Gedalof et al. 2005; Gochis et al. 2015; Milrad et al. 2015; Antokhina et al. 2016; Yang et al. 2021). All these processes combined during PPA further enhance the FWI and its components; thus, a PPA generates surface conditions that can sustain and contribute to the rapid spread of wildfires, whether they are caused by lightning or humans.
c. PPA expansion and mechanisms
The increase in PPA size observed here is consistent with the findings of Lyon et al. (2019) and Nabizadeh et al. (2019); these studies indicate an increase in the size of the Northern Hemisphere blocking events and heatwaves over the contiguous United States under future climate change scenarios. The increase in PPA size affects the spatial and temporal characteristics of extreme fire weather conditions and exposes larger areas to increased fire danger and for a longer duration. Furthermore, from a fire management perspective, the expansion of the PPA area makes it harder to suppress or control fires because there could be increased fire activity in a larger domain within the same time period.
Changes in PPA characteristics may result from either dynamic changes in the atmosphere (i.e., frequency of anticyclonic episodes) or from thermodynamic changes related to global warming and/or other land-use changes (Horton et al. 2015). We hypothesize that the significant expansion of the PPA observed in this study is likely due to the thermal expansion of the atmosphere associated with background warming. To verify this hypothesis, we first calculated the trend on Z500 anomalies and Z thickness field (difference in Z at 500 hPa and 1000 hPa) and found that they exhibit very similar trends implying that the changes observed in the Z500 field are largely due to thermodynamic changes (Fig. 12). Then, we detrended the Z500 field and recalculated PPA characteristics and their trends. We did not find significant trends in PPA area or other characteristics when the detrended Z500 field was used (Fig. 13), which further confirms that the trends observed in PPA characteristics are due to thermodynamic changes in the atmosphere. Thus, the expansion is likely associated with an increase in lower tropospheric temperature in recent decades (Vinnikov and Grody 2003).

Trend in May–September daily (a) Z500 anomalies and (b) the Z thickness (difference in 500-hPa and 1000-hPa Z field), 1979–2020. The black dots indicate grid cells with a significant trend magnitude ( p < 0.1).
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Trend in May–September daily (a) Z500 anomalies and (b) the Z thickness (difference in 500-hPa and 1000-hPa Z field), 1979–2020. The black dots indicate grid cells with a significant trend magnitude ( p < 0.1).
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Trend in May–September daily (a) Z500 anomalies and (b) the Z thickness (difference in 500-hPa and 1000-hPa Z field), 1979–2020. The black dots indicate grid cells with a significant trend magnitude ( p < 0.1).
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Trends in average PPA (a) number, (b) duration, (c) area, and (d) magnitude, 1979–2020 using the detrended Z500. The red lines represent trends while inset on each plot lists trend magnitude and p values; p < 0.1 indicates statistically significant trends.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1

Trends in average PPA (a) number, (b) duration, (c) area, and (d) magnitude, 1979–2020 using the detrended Z500. The red lines represent trends while inset on each plot lists trend magnitude and p values; p < 0.1 indicates statistically significant trends.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
Trends in average PPA (a) number, (b) duration, (c) area, and (d) magnitude, 1979–2020 using the detrended Z500. The red lines represent trends while inset on each plot lists trend magnitude and p values; p < 0.1 indicates statistically significant trends.
Citation: Journal of Climate 35, 19; 10.1175/JCLI-D-21-0926.1
In the context of inconsistent findings on what drives changes in persistent extreme weather events in midlatitude regions of the Northern Hemisphere (Francis and Vavrus 2012; Blackport and Screen 2020; Dai and Song 2020; Francis et al. 2020; Riboldi et al. 2020; Sun et al. 2022), we observe that the trends in PPA characteristics such as area and strength across western North America are primarily due to background warming. This agrees with previous work such as that of Riboldi et al. (2020), who found an increase in extreme temperatures due to slow progression of atmospheric flow; however, that study did not determine the role of AA in reducing the progression of atmospheric flow both in winter and summer months. Similarly, Blackport and Screen (2020) report that the changes in midlatitude temperature extremes are likely due to thermodynamic effects and not as a forced response of tropospheric waviness to AA.
6. Conclusions
While the presence of atmospheric blocking ridges and extreme temperatures that include heatwave events has been studied previously (Sillmann and Croci-Maspoli 2009; Trenberth et al. 2015; Grotjahn et al. 2016; Horton et al. 2016; Rodrigues and Woolings 2017; Mann et al. 2018; Parente et al. 2018; Sousa et al. 2018; Gibson et al. 2020; Abatzoglou et al. 2021b; Jain and Flannigan 2021), this study expands our understanding of extended summertime persistent positive midtropospheric height anomalies and its association to both wildfire and fire weather in western North America. Although our study is limited to western North America, we expect similar relationships between PPA and fire activity in areas where summertime PPA events frequently occur. This study’s findings may improve our understanding on the connections between PPA and wildfires, signal how it may change in future warmer climate scenarios, and enhance the predictability of extreme wildfire events in this region. In addition, the quantification of PPA-related fire ignition may help to improve fire weather predictability with improvements in models that are able to forecast PPAs.
To conclude, this analysis reveals that the May–September PPAs substantially impact fire weather conditions in western North America, with greater influence at higher latitudes. Furthermore, significant spatial expansion of PPA events over 1979–2020 exposes larger areas and populations of western North America to amplified fire weather conditions that are expected to continue expanding due to climate change. Specific key findings are as follows:
-
On average, 26 out of 153 days of May–September experience PPAs over western North America.
-
40% (60%) of wildfire events ignitions with burned area > 500 ha (>20 000 ha) are associated with PPAs during 2003–18.
-
There is a presence of an open-wave 500-hPa ridge over the Canadian Rockies during the majority of PPA days with fire ignitions.
-
The presence of PPAs significantly increases the odds of occurrence of fire ignition across most of western North America by a factor of 7; this factor is greater in higher latitudes than lower latitudes.
-
A strong interannual positive correlation exists between the percentage of PPA days and the average fire burned areas during 1997–2016.
-
The distribution of PPA-day anomalies of temperature, VPD, FFMC, DMC, DC, ISI, and FWI (precipitation, RH, wind speed) are significantly higher (lower) than that of non-PPA days anomalies.
-
The temperature and VPD (RH) reaches a maximum (minimum) on the same day as maximum PPA strength, whereas the precipitation anomalies are minimum two days earlier and the FFMC, DMC, DC, ISI, and FWI are maximum 2–6 days later than the maximum PPA strength.
-
The number of PPA events, their duration, and magnitude have not changed significantly between 1979 and 2020; however, the area of PPAs has increased significantly during this time period.
-
The spatial expansion of PPA events is largely due to thermodynamic changes in the lower atmosphere in recent decades.
This work establishes the relationship between PPA and wildfires for the historical period of 1979–2020 over western North America. However it should be noted that there exist many drivers for fire ignition and spread in addition to fire-related variables discussed in this paper such as prolonged drought, tree composition and mortality, orographic controls, and fire management strategies (Gedalof et al. 2005; Schoennagel et al. 2017). However, these variables are not directly influenced by the presence of a PPA within synoptic-scale timeframe, and are therefore not discussed in detail here. Further investigation of the influence of PPAs on extreme fire weather indices, the relationship between the occurrence of PPAs, teleconnections, and fire events, the combined influence of PPA and prolonged drought on fire weather, and PPA fire weather relations under future climate scenarios will help to improve the predictability of wildfire events in the context of anthropogenic global warming.
Acknowledgments.
Thanks to Niels Andela for providing additional two years of the Global Fire Atlas (GFA) data in addition to a publicly available database and to Dante Castellanos for help on extending the FWI data used here. Thanks to three anonymous referees who provided constructive comments on this work that led to a much improved paper.
Data availability statement.
All datasets used in this paper are publicly available through the references cited and can be accessed through the web portals https://cds.climate.copernicus.eu/cdsapp#!/search?type=dataset (ERA5 data, accessed May 2020), https://www.globalfiredata.org/ (GFA and GEFD4.1 data, accessed May 2020), and https://zenodo.org/record/3626193#.X9pTFdhKg4s (FWI data, accessed June 2020). The R algorithm to identify and track PPA events is available by request to corresponding author.
REFERENCES
Abatzoglou, J. T., and A. P. Williams, 2016: Impact of anthropogenic climate change on wildfire across western US forests. Proc. Natl. Acad. Sci. USA., 113, 11 770–11 775, https://doi.org/10.1073/pnas.1607171113.
Abatzoglou, J. T., D. S. Battisti, A. P. Williams, W. D. Hansen, B. J. Harvey, and C. A. Kolden, 2021a: Projected increases in western US forest fire despite growing fuel constraints. Commun. Earth Environ., 2, 227, https://doi.org/10.1038/s43247-021-00299-0.
Abatzoglou, J. T., D. E. Rupp, L. W. O’Neill, and M. Sadegh, 2021b: Compound extremes drive the western Oregon wildfires of September 2020. Geophys. Res. Lett., 48, e2021GL092520, https://doi.org/10.1029/2021GL092520.
Andela, N., D. C. Morton, L. Giglio, R. Paugam, Y. Chen, S. Hantson, G. R. van der Werf, and J. T. Randerson, 2019: The Global Fire Atlas of individual fire size, duration, speed, and direction. Earth Syst. Sci. Data, 11, 529–552, https://doi.org/10.5194/essd-11-529-2019.
Antokhina, O. Y., P. N. Antokhin, Y. V. Martynova, and V. I. Mordvinov, 2016: The impact of atmospheric blocking on spatial distributions of summertime precipitation over Eurasia. IOP Conf. Ser. Earth Environ. Sci., 48, 012035, https://doi.org/10.1088/1755-1315/48/1/012035.
Balch, J. K., B. A. Bradley, J. T. Abatzoglou, R. C. Nagy, E. J. Fusco, and A. L. Mahood, 2017: Human-started wildfires expand the fire niche across the United States. Proc. Natl. Acad. Sci. USA, 114, 2946–2951, https://doi.org/10.1073/pnas.1617394114.
Balshi, M. S., A. D. McGuire, P. Duffy, M. Flannigan, J. Walsh, and J. Melillo, 2009: Assessing the response of area burned to changing climate in western boreal North America using a Multivariate Adaptive Regression Splines (MARS) approach. Global Change Biol., 15, 578–600, https://doi.org/10.1111/j.1365-2486.2008.01679.x.
Barnes, E. A., J. Slingo, and T. Woollings, 2012: A methodology for the comparison of blocking climatologies across indices, models and climate scenarios. Climate Dyn., 38, 2467–2481, https://doi.org/10.1007/s00382-011-1243-6.
Barnes, E. A., E. Dunn-Sigouin, G. Masato, and T. Woollings, 2014: Exploring recent trends in Northern Hemisphere blocking. Geophys. Res. Lett., 41, 638–644, https://doi.org/10.1002/2013GL058745.
Barriopedro, D., R. García-Herrera, A. R. Lupo, and E. Hernández, 2006: A climatology of Northern Hemisphere blocking. J. Climate, 19, 1042–1063, https://doi.org/10.1175/JCLI3678.1.
BC Gov News, 2021: Provincial state of emergency declared. BC Gov News, Government of British Columbia, accessed 30 July 2021, https://news.gov.bc.ca/releases/2021PSSG0047-001404.
Blackport, R., and J. A. Screen, 2020: Insignificant effect of Arctic amplification on the amplitude of midlatitude atmospheric waves. Sci. Adv., 6, eaay2880, https://doi.org/10.1126/sciadv.aay2880.
Brewer, M. J., and C. B. Clements, 2020: The 2018 Camp Fire: Meteorological analysis using in situ observations and numerical simulations. Atmosphere, 11, 47, https://doi.org/10.3390/atmos11010047.
Bronaugh, D., and A. Werner, 2019: zyp: Zhang + Yue-Pilon Trends, version 0.10-1.1. R package, https://rdrr.io/cran/zyp/.
Csardi, G., and T. Nepusz, 2006: The igraph software package for complex network research. InterJournal Complex Syst., 1695, 1–9, https://www.bibsonomy.org/bibtex/2a3db966c7a0ce1715c56ce04c806afd9/yourwelcome.
Dai, A., and M. Song, 2020: Little influence of Arctic amplification on mid-latitude climate. Nat. Climate Change, 10, 231–237, https://doi.org/10.1038/s41558-020-0694-3.
Davini, P., C. Cagnazzo, S. Gualdi, and A. Navarra, 2012: Bidimensional diagnostics, variability, and trends of Northern Hemisphere blocking. J. Climate, 25, 6496–6509, https://doi.org/10.1175/JCLI-D-12-00032.1.
Dole, R. M., and N. D. Gordon, 1983: Persistent anomalies of the extratropical Northern Hemisphere wintertime circulation: Geographical distribution and regional persistence characteristics. Mon. Wea. Rev., 111, 1567–1586, https://doi.org/10.1175/1520-0493(1983)111<1567:PAOTEN>2.0.CO;2.
Dole, R. M., and Coauthors, 2011: Was there a basis for anticipating the 2010 Russian heat wave? Geophys. Res. Lett., 38, L06702, https://doi.org/10.1029/2010GL046582.
Drouard, M., and T. Woollings, 2018: Contrasting mechanisms of summer blocking over western Eurasia. Geophys. Res. Lett., 45, 12 040–12 048, https://doi.org/10.1029/2018GL079894.
Dunn-Sigouin, E.,-S. W. Son, and H. Lin, 2013: Evaluation of Northern Hemisphere blocking climatology in the global environment multiscale model. Mon. Wea. Rev., 141, 707–727, https://doi.org/10.1175/MWR-D-12-00134.1.
Elliott, R. D., and T. B. Smith, 1949: A study of the effect of large blocking highs on the general circulation in the Northern Hemisphere westerlies. J. Appl. Meteor., 6, 68–85, https://doi.org/10.1175/1520-0469(1949)006<0068:ASOTEO>2.0.CO;2.
Fauria, M. M., and E. A. Johnson, 2008: Climate and wildfires in the North American boreal forest. Philos. Trans. Roy. Soc., 363B, 2317–2329, https://doi.org/10.1098/rstb.2007.2202.
Flannigan, M. D., and J. B. Harrington, 1988: A study of the relation of meteorological variables to monthly provincial area burned by wildfires in Canada (1953–80). J. Appl. Meteor., 27, 441–452, https://doi.org/10.1175/1520-0450(1988)027<0441:ASOTRO>2.0.CO;2.
Flannigan, M. D., K. A. Logan, B. D. Amiro, W. R. Skinner, and B. J. Stocks, 2005: Future area burned in Canada. Climatic Change, 72, 1–16, https://doi.org/10.1007/s10584-005-5935-y.
Francis, J. A., and S. J. Vavrus, 2012: Evidence linking Arctic amplification to extreme weather in mid-latitudes. Geophys. Res. Lett., 39, L06801, https://doi.org/10.1029/2012GL051000.
Francis, J. A., and S. J. Vavrus, 2015: Evidence for a wavier jet stream in response to rapid Arctic warming. Environ. Res. Lett., 10, 014005, https://doi.org/10.1088/1748-9326/10/1/014005.
Francis, J. A., N. Skific, and S. J. Vavrus, 2018: North American weather regimes are becoming more persistent: Is Arctic amplification a factor? Geophys. Res. Lett., 45, 11 414–11 422, https://doi.org/10.1029/2018GL080252.
Francis, J. A., N. Skific, and S. J. Vavrus, 2020: Increased persistence of large-scale circulation regimes over Asia in the era of amplified Arctic warming, past and future. Sci. Rep., 10, 14953, https://doi.org/10.1038/s41598-020-71945-4.
Gedalof, Z., D. L. Peterson, and N. J. Mantua, 2005: Atmospheric, climatic, and ecological controls on extreme wildfire years in the northwestern United States. Ecol. Appl., 15, 154–174, https://doi.org/10.1890/03-5116.
Gibson, P. B., D. E. Waliser, B. Guan, M. J. DeFlorio, F. M. Ralph, and D. L. Swain, 2020: Ridging associated with drought across the western and southwestern United States: Characteristics, trends, and predictability sources. J. Climate, 33, 2485–2508, https://doi.org/10.1175/JCLI-D-19-0439.1.
Giglio, L., J. T. Randerson, and G. R. van der Werf, 2013: Analysis of daily, monthly, and annual burned area using the fourth-generation Global Fire Emissions Database (GFED4). J. Geophys. Res. Biogeosci., 118, 317–328, https://doi.org/10.1002/jgrg.20042.
Gillett, N. P., A. J. Weaver, F. W. Zwiers, and M. D. Flannigan, 2004: Detecting the effect of climate change on Canadian forest fires. Geophys. Res. Lett., 31, L18211, https://doi.org/10.1029/2004GL020876.
Gochis, D., and Coauthors, 2015: The great Colorado flood of September 2013. Bull. Amer. Meteor. Soc., 96, 1461–1487, https://doi.org/10.1175/BAMS-D-13-00241.1.
Grotjahn, R., and Coauthors, 2016: North American extreme temperature events and related large scale meteorological patterns: A review of statistical methods, dynamics, modeling, and trends. Climate Dyn., 46, 1151–1184, https://doi.org/10.1007/s00382-015-2638-6.
Hersbach, H., and Coauthors, 2020: The ERA5 global reanalysis. Quart. J. Roy. Meteor. Soc., 146, 1999–2049, https://doi.org/10.1002/qj.3803.
Hijmans, R. J., and J. van Etten, 2012: raster: Geographic Analysis and Modeling. with raster data version 3.5-21. R package, https://rdrr.io/cran/raster/.
Holzworth, R. H., J. B. Brundell, M. P. McCarthy, A. R. Jacobson, C. J. Rodger, and T. S. Anderson, 2021: Lightning in the Arctic. Geophys. Res. Lett., 48, e2020GL091366, https://doi.org/10.1029/2020GL091366.
Horton, D. E., N. C. Johnson, D. Singh, D. L. Swain, B. Rajaratnam, and N. S. Diffenbaugh, 2015: Contribution of changes in atmospheric circulation patterns to extreme temperature trends. Nature, 522, 465–469, https://doi.org/10.1038/nature14550.
Horton, R. M., J. S. Mankin, C. Lesk, E. Coffel, and C. Raymond, 2016: A review of recent advances in research on extreme heat events. Curr. Climate Change Rep., 2, 242–259, https://doi.org/10.1007/s40641-016-0042-x.
Jain, P., and M. Flannigan, 2021: The relationship between the polar jet stream and extreme wildfire events in North America. J. Climate, 34, 6247–6265, https://doi.org/10.1175/JCLI-D-20-0863.1.
Johnson, E. A., and D. R. Wowchuk, 1993: Wildfires in the southern Canadian Rocky Mountains and their relationship to mid-tropospheric anomalies. Can. J. For. Res., 23, 1213–1222, https://doi.org/10.1139/x93-153.
Kendall, M. G., 1975: Rank Correlation Methods. 4th ed. Charles Griffin, 210 pp.
Kharuk, V. I., E. I. Ponomarev, G. A. Ivanova, M. L. Dvinskaya, S. C. P. Coogan, and M. D. Flannigan, 2021: Wildfires in the Siberian taiga. Ambio, 50, 1953–1974, https://doi.org/10.1007/s13280-020-01490-x.
Kirchmeier-Young, M. C., N. P. Gillett, F. W. Zwiers, A. J. Cannon, and F. S. Anslow, 2019: Attribution of the influence of human-induced climate change on an extreme fire season. Earth’s Future, 7, 2–10, https://doi.org/10.1029/2018EF001050.
Kohonen, T., 2001: Self-Organazing Maps. 3rd ed. Springer, 502 pp.
Lau, N.-C., and M. J. Nath, 2012: A model study of heat waves over North America: Meteorological aspects and projections for the twenty-first century. J. Climate, 25, 4761–4784, https://doi.org/10.1175/JCLI-D-11-00575.1.
Liu, Q., 1994: On the definition and persistence of blocking. Tellus, 46, 286–298, https://doi.org/10.1034/j.1600-0870.1994.t01-2-00004.x.
Lupo, A. R., A. D. Jensen, I. I. Mokhov, A. V. Timazhev, T. Eichler, and B. Efe, 2019: Changes in global blocking character in recent decades. Atmosphere, 10, 92, https://doi.org/10.3390/atmos10020092.
Lyon, B., A. G. Barnston, E. Coffel, and R. M. Horton, 2019: Projected increase in the spatial extent of contiguous US summer heat waves and associated attributes. Environ. Res. Lett., 14, 114029, https://doi.org/10.1088/1748-9326/ab4b41.
Mann, H. B., 1945: Nonparametric test against trend. Econometrica, 13, 245–259, https://doi.org/10.2307/1907187.
Mann, M. E., S. Rahmstorf, K. Kornhuber, B. A. Steinman, S. K. Miller, S. Petri, and D. Coumou, 2018: Projected changes in persistent extreme summer weather events: The role of quasi-resonant amplification. Sci. Adv., 4, eaat3272, https://doi.org/10.1126/sciadv.aat3272.
McElhinny, M., J. F. Beckers, C. Hanes, M. Flannigan, and P. Jain, 2020: A high-resolution reanalysis of global fire weather from 1979 to 2018—Overwintering the drought code. Earth Syst. Sci. Data, 12, 1823–1833, https://doi.org/10.5194/essd-12-1823-2020.
Miller, R. L., G. M. Lackmann, and W. A. Robinson, 2020: A new variable-threshold persistent anomaly index: Northern Hemisphere anomalies in the ERA-Interim reanalysis. Mon. Wea. Rev., 148, 43–62, https://doi.org/10.1175/MWR-D-19-0144.1.
Milrad, S. M., J. R. Gyakum, and E. H. Atallah, 2015: A meteorological analysis of the 2013 Alberta flood: Antecedent large-scale flow pattern and synoptic-dynamic characteristics. Mon. Wea. Rev., 143, 2817–2841, https://doi.org/10.1175/MWR-D-14-00236.1.
Miralles, D. G., A. J. Teuling, C. C. van Heerwaarden, and J. V. G. de Arellano, 2014: Mega-heatwave temperatures due to combined soil desiccation and atmospheric heat accumulation. Nat. Geosci., 7, 345–349, https://doi.org/10.1038/ngeo2141.
Nabizadeh, E., P. Hassanzadeh, D. Yang, and E. A. Barnes, 2019: Size of the atmospheric blocking events: Scaling law and response to climate change. Geophys. Res. Lett., 46, 13 488–13 499, https://doi.org/10.1029/2019GL084863.
Parente, J., M. G. Pereira, M. Amraoui, and E. M. Fischer, 2018: Heat waves in Portugal: Current regime, changes in future climate and impacts on extreme wildfires. Sci. Total Environ., 631–632, 534–549, https://doi.org/10.1016/j.scitotenv.2018.03.044.
Parsons, S., J. A. Renwick, and A. J. McDonald, 2016: An assessment of future Southern Hemisphere blocking using CMIP5 projections from four GCMs. J. Climate, 29, 7599–7611, https://doi.org/10.1175/JCLI-D-15-0754.1.
Pelly, J. L., and B. J. Hoskins, 2003: A new perspective on blocking. J. Atmos. Sci., 60, 743–755, https://doi.org/10.1175/1520-0469(2003)060<0743:ANPOB>2.0.CO;2.
Petoukhov, V., S. Petri, K. Kornhuber, K. Thonicke, D. Coumou, and H. J. Schellnhuber, 2018: Alberta wildfire 2016: Apt contribution from anomalous planetary wave dynamics. Sci. Rep., 8, 12375, https://doi.org/10.1038/s41598-018-30812-z.
Philip, S. Y., and Coauthors, 2021: Rapid attribution analysis of the extraordinary heatwave on the Pacific Coast of the US and Canada June 2021. Earth Syst. Dyn., https://doi.org/10.5194/esd-2021-90, in press.
Pinheiro, M. C., P. A. Ullrich, and R. Grotjahn, 2019: Atmospheric blocking and intercomparison of objective detection methods: Flow field characteristics. Climate Dyn., 53, 4189–4216, https://doi.org/10.1007/s00382-019-04782-5.
Rasilla, D. F., J. C. García-Codron, V. Carracedo, and C. Diego, 2010: Circulation patterns, wildfire risk and wildfire occurrence at continental Spain. Phys. Chem. Earth, 35, 553–560, https://doi.org/10.1016/j.pce.2009.09.003.
R Core Team, 2019: A language and environment for statistical computing. R Foundation for Statistical Computing, https://www.R-project.org/.
Renwick, J. A., 2005: Persistent positive anomalies in the Southern Hemisphere circulation. Mon. Wea. Rev., 133, 977–988, https://doi.org/10.1175/MWR2900.1.
Riboldi, J., F. Lott, F. D’Andrea, and G. Rivière, 2020: On the linkage between Rossby wave phase speed, atmospheric blocking, and Arctic amplification. Geophys. Res. Lett., 47, e2020GL087796, https://doi.org/10.1029/2020GL087796.
Rodrigues, R. R., and T. Woollings, 2017: Impact of atmospheric blocking on South America in austral summer. J. Climate, 30, 1821–1837, https://doi.org/10.1175/JCLI-D-16-0493.1.
Schaller, N., J. Sillmann, J. Anstey, E. M. Fischer, C. M. Grams, and S. Russo, 2018: Influence of blocking on northern European and western Russian heatwaves in large climate model ensembles. Environ. Res. Lett., 13, 054015, https://doi.org/10.1088/1748-9326/aaba55.
Schoennagel, T., and Coauthors, 2017: Adapt to more wildfire in western North American forests as climate changes. Proc. Natl. Acad. Sci. USA., 114, 4582–4590, https://doi.org/10.1073/pnas.1617464114.
Schwierz, C., M. Croci-Maspoli, and H. C. Davies, 2004: Perspicacious indicators of atmospheric blocking. Geophys. Res. Lett., 31, L06125, https://doi.org/10.1029/2003GL019341.
Screen, J. A., 2014: Arctic amplification decreases temperature variance in northern mid- to high-latitudes. Nat. Climate Change, 4, 577–582, https://doi.org/10.1038/nclimate2268.
Sen, P. K., 1968: Estimates of the regression coefficient based on Kendall’s tau. J. Amer. Stat. Assoc., 63, 1379–1389, https://doi.org/10.1080/01621459.1968.10480934.
Sillmann, J., and M. Croci-Maspoli, 2009: Present and future atmospheric blocking and its impact on European mean and extreme climate. Geophys. Res. Lett., 36, L10702, https://doi.org/10.1029/2009GL038259.
Skinner, W. R., and Coauthors, 2002: A 500 hPa synoptic wildland fire climatology for large Canadian forest fires, 1959–1996. Theor. Appl. Climatol., 71, 157–169, https://doi.org/10.1007/s007040200002.
Small, D., E. Atallah, and J. R. Gyakum, 2014: An objectively determined blocking index and its Northern Hemisphere climatology. J. Climate, 27, 2948–2970, https://doi.org/10.1175/JCLI-D-13-00374.1.
Sousa, P. M., R. M. Trigo, D. Barriopedro, P. M. M. Soares, and J. A. Santos, 2018: European temperature responses to blocking and ridge regional patterns. Climate Dyn., 50, 457–477, https://doi.org/10.1007/s00382-017-3620-2.
Sousa, P. M., D. Barriopedro, R. García-Herrera, T. Woollings, and R. M. Trigo, 2021: A new combined detection algorithm for blocking and subtropical ridges. J. Climate, 34, 7735–7758, https://doi.org/10.1175/JCLI-D-20-0658.1.
Steinfeld, D., and S. Pfahl, 2019: The role of latent heating in atmospheric blocking dynamics: A global climatology. Climate Dyn., 53, 6159–6180, https://doi.org/10.1007/s00382-019-04919-6.
Stephenson, D. B., 2000: Use of the “odds ratio” for diagnosing forecast skill. Wea. Forecasting, 15, 221–232, https://doi.org/10.1175/1520-0434(2000)015<0221:UOTORF>2.0.CO;2.
Sun, X., and Coauthors, 2022: Enhanced jet stream waviness induced by suppressed tropical Pacific convection during boreal summer. Nat. Commun., 13, 1288, https://doi.org/10.1038/s41467-022-28911-7.
Tan, X., S. Chen, T. Y. Gan, B. Liu, and X. Chen, 2019: Dynamic and thermodynamic changes conducive to the increased occurrence of extreme spring fire weather over western Canada under possible anthropogenic climate change. Agric. For. Meteor., 265, 269–279, https://doi.org/10.1016/j.agrformet.2018.11.026.
Tibaldi, S., and F. Molteni, 1990: On the operational predictability of blocking. Tellus, 42A, 343–365, https://doi.org/10.1034/j.1600-0870.1990.t01-2-00003.x.
Trenberth, K. E., J. T. Fasullo, and T. G. Shepherd, 2015: Attribution of climate extreme events. Nat. Climate Change., 5, 725–730, https://doi.org/10.1038/nclimate2657.
Tsonis, A. A., 2013: An Introduction to Atmospheric Thermodynamics. 2nd ed. Cambridge University Press, 199 pp.
Tymstra, C., B. J. Stocks, X. Cai, and M. D. Flannigan, 2020: Wildfire management in Canada: Review, challenges, and opportunities. Prog. Disaster Sci., 5, 100045, https://doi.org/10.1016/j.pdisas.2019.100045.