Impacts of the North Atlantic Oscillation on winter precipitations and storm track variability in southeast Canada and the northeast United States

The North Atlantic Oscillation (NAO) affects atmospheric variability from eastern North America to Europe. Although the link between the NAO and winter precipitations in eastern North America has been the focus of previous work, only few studies have considered extreme precipitation and hitherto provided clear physical explanations on these relationships. In this study we revisit and extend the analysis of the effect of the NAO on mean and heavy winter precipitations over a large domain covering southeast Canada and the northeastern United States. Furthermore, we use the recent ERA5 reanalysis dataset (1979–2018), which currently has the highest available horizontal resolution for a global reanalysis (0.25), to track extratropical cyclones to delve into the physical processes behind the relationship between NAO and precipitation, snowfall, snowfall-to-precipitation ratio (S/P ), and snow cover depth anomalies in the region. In particular, our results show that positive NAO phases are associated with less snowfall over a wide region covering Nova Scotia, New England and the Mid-Atlantic of the United States relative to negative NAO phases. Over the same area, the analysis of heavy snowfall revealed that there are up to twice as many heavy snowfall events during negative phases compared to positive phases. Therefore, a significant negative correlation is also seen between S/P and the NAO over this region. This is due to a decrease (increase) in cyclogenesis of coastal storms near the United States east coast during positive (negative) NAO phases, as well as a northward (southward) displacement of the mean storm track over North America.


Introduction
The North Atlantic Oscillation (NAO) is the dominant mode of atmospheric variability in the North Atlantic (Hurrell, 1995). The NAO refers to swings in the atmospheric pressure difference between the Icelandic low and Azores high and is a key factor in the cool-season climate variability from the eastern coast of the United States to Siberia and from the Arctic to the subtropical Atlantic. A common measure of the NAO phase is the so-called NAO index (NAOI) that is commonly defined as the difference in the normalized sea level pressure (SLP) anomalies between Stykkishólmur/Reykjavik, Iceland, and either Lisbon, Portugal, or Ponta Delgada, Azores (Hurrell, 1995). Positive NAO phases are associated with a deepening of the Icelandic low and a strengthening of the Azores high. The increased SLP gradient consequently leads to an enhanced westerly flow and a northward shift of the midlatitude storm track (Rogers, 1990;Hurrell and Van Loon, 1997). The negative phase of the NAO is associated with weaker westerlies, an increase in meandering of the jet stream and higher-than-average pressure over Greenland and Iceland, which often form blocking patterns over the North Atlantic (Shabbar et al., 2001).
While the link between the NAO and winter precipitations in northern Europe and in the Mediterranean region is moderately strong (Hurrell, 1995), the effects of the NAO phase on precipitation and snowfall amounts in eastern North America are known to be weaker (e.g., Bradbury et al., 2003). Various studies have shown that a slight negative correlation exists between the winter NAO phase and total winter precipitation in New England (Bradbury et al., 2002a;Ning and Bradley, 2015). A significant negative correlation be-J. Chartrand and F. S. R. Pausata: Impacts of the NAO on precipitation and storm track variability tween NAO and seasonal snowfall in New England (Hartley and Keables, 1998) and in other parts of northeastern US has also been reported (Kocin and Uccellini, 2004;Notaro et al., 2006;Morin et al., 2008). Moreover, Huntington et al. (2004) found a significant negative correlation between the NAO and the snowfall-to-precipitation ratio (S/P ) in the northeastern US. Kocin and Uccellini (2004) have studied the snowfall climatology of the northeast urban corridor in depth and have noticed that winters with a low NAO almost always coincide with the high-snowfall years. While the low-frequency variations in the NAO proved to be quite well correlated with seasonal snowfall in that region, their study also demonstrated that the occurrence of high-impact winter storms is strongly linked to the daily NAO value. However, Archambault et al. (2008) and Notaro et al. (2006) have found that precipitation event frequency in the northeastern US is slightly higher in positive NAO conditions. Several studies (Stone et al., 2000;Bonsal and Shabbar, 2008;Whan and Zwiers, 2017) have investigated the relationship between precipitation and teleconnection patterns in Canada but have found that winter precipitation is generally poorly correlated with the NAO in eastern Canada.
As extratropical cyclones account for almost all of the winter precipitation in eastern North America (Pfahl and Wernli, 2012), the variability in winter precipitation and snowfall can be explained by the shift in the storm track over North America and the North Atlantic. Several studies have focused specifically on the storm track variability in the North Atlantic and showed that the positive (negative) phase of the NAO favors a more northern (southern) storm track over the North Atlantic (Rogers, 1990;Serreze et al., 1997;Riviere and Orlanski, 2007;Pinto et al., 2009), as well as an eastward displacement of the mean trough axis over North America (Bradbury et al., 2002b). In a broader region,  have found that a strong positive NAO is always associated simultaneously with more frequent cyclone activity in the high Canadian Arctic and less frequent activity on the east coast.
In this study we investigate the effect of the NAO on precipitation and snowfall over a wide domain covering southeastern Canada as well as the northeastern US. Although previous studies have delved into the link between the NAO and winter precipitation in eastern North America, they have not considered extreme precipitation and snowfall. Knowing that extreme winter storms are some of the most damaging meteorological phenomena that occur in these regions, it is highly relevant to look at the changes in heavy precipitation event occurrence as well. Furthermore, we use the recent highresolution ERA5 reanalysis dataset  to track extratropical cyclones and to calculate the blocking frequency over the North Atlantic Basin. In doing so, we aim at providing a physical explanation for the precipitation and snowfall variability based on the storm track variability associated with the phases of the NAO.
The paper is organized as follows. Data and methodology are described and discussed in Sect. 2. Section 3 presents the results obtained from statistical analysis of the winter climate and storm track variability with the NAO. Section 4 presents a discussion of the results. Finally, the main findings are presented in Sect. 5, along with the conclusion.
2 Data and methods

Datasets
The NAOI used in this analysis is the monthly stationbased index (Hurrell, 1995) based on the difference of normalized SLP between Lisbon, Portugal, and Stykkishólmur/Reykjavik in Iceland (taken from the National Center for Atmospheric Research (NCAR) website: https://climatedataguide.ucar.edu/climate-data/ hurrell-north-atlantic-oscillation-nao-index-station-based, last access: 24 November 2020; Hurrell and National Center for Atmospheric Research Staff, 2020). Positive values of the NAOI are given by higher-than-average SLP over Portugal and lower-than-average SLP over Iceland. As in previous studies (Ning and Bradley, 2015;Bradbury et al., 2002a;Huntington et al., 2004), we consider the December to March period (DJFM), because consistent snowfall and winter conditions usually occur during these 4 months over most of the study area. The monthly NAOI has been used to calculate the link between winter precipitations and the NAO.
The daily snowfall, total precipitation and snowpack depth in snow water equivalent (SWE) data come from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 high-resolution (0.25 • ) reanalysis dataset for the period 1979period -2018period (Hersbach et al., 2020. The precipitation and snowfall totals used for calculating correlations are the summed total precipitation and snowfall for each month at every grid point of the domain of interest. The change in heavy precipitation and snowfall events has been analyzed by using various extreme climate indices (R95N, R10mm, Rx1day, Rx5days, R20mm). However, the results are very similar among these indices within our domain. Therefore, for clarity in this paper we only show the results for the relationship between the NAO and the number of days with precipitation (and snowfall) above the 95th percentile of climatological daily precipitation (snowfall) during the cold season.
The S/P was calculated by dividing snowfall by total precipitation, as in Huntington et al. (2004). To analyze how the annual maximum of snowpack depth varies depending on the winter (DJFM) NAO phase, the seasonal maximum of SWE was taken over the course of the snow year (1 August to 31 July), as the maximum of snow depth is sometimes achieved later than March in some parts of the domain. The winter climatology of total precipitation, snowfall, S/P and maximum seasonal SWE using the ERA5 dataset are shown in Fig. 1, as well as the 95th percentile of daily precipitation and snowfall.
Linear correlation coefficients with the monthly NAOI were first calculated for mean and extreme precipitation, snowfall, and S/P . For the snowpack SWE, correlation coefficients were calculated between maximum SWE and seasonal (DJFM) NAO. Correlations were calculated using the Pearson R-value method at every grid point. The statistical significance of these correlations using the p value was determined at the 0.05 level. For statistical comparison between high-NAO and low-NAO months, the variables were averaged for the 40 months with the highest NAO (first quartile) and for the 40 months with the lowest NAO (fourth quartile). For the snowpack maximum SWE analysis, the 10 winters with the highest NAOI were compared to the 10 years with the lowest NAOI. To measure the statistical significance of these comparisons, the Student t test was used with a significance level of 0.05.
Daily precipitation and snowfall data from the Global Historical Climatology Network (GHCN-Daily) served as validation data for the same period as the ERA5 dataset . The data from 29 stations well distributed in the domain have been analyzed. Of these stations, 14 are located in eastern Canada, 6 are located in New England and 9 are located in the Mid-Atlantic region of the United States (Fig. 2). These stations were chosen based on the availability of data in the 1979-2018 period. Each station has near-complete precipitation data in this time period, which is essential for the comparison between results. Many stations that had data covering the whole time period of interest were rejected because of the high percentage of missing data during some years. Station inclusion criteria were intended to prevent stations with short, nonrepresentative data records from biasing the analysis. The same calculations on monthly precipitation, snowfall and S/P that have been done on ERA5 gridded data have been performed using the station-based data. Comparison between station-based and reanalysis data also helped in validating the ERA5 data on a regional scale.

Storm tracking algorithm
In order to track individual cyclone trajectories on the ERA5 dataset, we adapted the existing storm tracking algorithms by Murray and Simmonds (1991) and Hanley and Caballero (2012). We modified these methods rather than directly use one of them to improve the algorithm performance in tracking the extratropical cyclones over our domain of interest (including over mountains), when using the highresolution and 3-hourly output of the ERA5 reanalysis. The main features of the algorithm are described here below.
Because of the high resolution of ERA5, we first use a Cressman smoother (400 km influence radius) to filter out small-scale stationary local minima in the SLP field in the intermountain regions and retain only synoptic-scale low-pressure systems. Local minima in the SLP field (SLP value is smaller than the eight neighboring grid points) are then taken as low-pressure centers for the tracking process.
Cyclone trajectories are created through linear projection: to forecast the second position in a cyclone trajectory, the 500 hPa wind field is used; afterward, the locations of future low-pressure centers are extrapolated (linearly projected) by using the difference between the previous and current center's locations. Matches between low-pressure centers at each time step are attempted by looking for low-pressure centers within a 350 km search radius from the predicted trajectories. If several matches are possible for a single trajectory, absolute departure from predicted position and pressure with predicted value is used for matching low-pressure centers. If no low-pressure centers are found within 350 km of a predicted cyclone position, this cyclone trajectory is completed. A low-pressure center that is not matched to any previous storm track is used as a first center for a new storm trajectory.
After each complete storm trajectory, two elimination criteria are applied in order to remove the stationary orographic features and short-lived features from the list of tracks found; a minimum total displacement of more than 10 • and a minimum lifetime of 48 h.
The cyclone tracking has been performed in a limited area of the ERA5 grid (180 • W-0 • E, 10-80 • N) to make calculations faster. As our study area is much smaller and far enough from the borders of the tracking area, this did not bias the analysis of storm tracks within our domain of interest. An average of 381 cyclone trajectories per winter were found (a total of 15 248 for the 1979-2018 period).
The track density is defined in this study as the number of cyclone trajectories per area of 1000 km 2 . Each cyclone trajectory is only counted once for each grid point near the trajectory. For the purpose of the study, we equate changes in the cyclone track density to changes in storm tracks. Cyclogenesis density is defined as the number of cyclogenesis events (first points in trajectories) per area of 1000 km 2 . The spatial distribution of both track density and cyclogenesis density climatology (Fig. 3) are well in agreement with previous work that used different tracking methods (e.g., Neu et al., 2013). In particular, the main region of lee cyclogenesis east of the Canadian and Colorado Rockies is very well represented, as well as the areas of maximum coastal cyclogenesis due to the land-sea temperature contrast just offshore of Cape Hatteras in North Carolina (Fig. 3a).

Blocking frequency
Blocks near Greenland (i.e., Greenland block) are known to be responsible for deeper troughs over eastern North America (e.g., Resio and Hayden, 1975), which consequently have an influence on cyclogenesis and cyclone trajectories. For that reason, we investigated the changes in blocking frequency during positive and negative NAO months over the North Atlantic Basin (30-70 • N, 85 • W-0 • E). The analysis of the blocking frequency using the ERA5 dataset was performed using a bidimensional index that identifies reversals in the meridional gradient of the 500 hPa geopotential height (Pausata et al., 2015;Anstey et al., 2013 ;Tibaldi and Molteni, 1990). At each grid point of latitude φ and longitude λ, the northward and southward meridional gradients of geopotential heights are respectively estimated following Eqs. (1) and (2): A blocking event at a grid point is diagnosed when two conditions are verified: (1) N (φ, λ) > 0, indicating a reversal of the climatological conditions with easterlies equatorward of the grid point, and (2) S (φ, λ) < −10 m •−1 , indicating westerlies poleward of the grid point. The blocking fre-  quency (%) is then obtained by the number of time steps that a grid point is defined as "blocked" divided by the total number of time steps.

Results
In this section we first present the relationship between the NAO and total precipitation, snowfall and S/P (Sect. 3.1): for each variable, the spatial pattern of correlations with the NAO index is discussed as well as the difference in monthly averages between high-NAO and low-NAO months. In Sect. 3.2, heavy precipitation and snowfall events are considered, as well as the link between seasonal NAO and maximum snowpack depth. Finally, in Sect. 3.3 we present the effect of the NAO on storm tracks and cyclogenesis in North America.
3.1 Changes in total precipitation, snowfall and snowfall-to-total-precipitation ratio For total precipitation, the results show that only a weak relationship with the NAO exists over most of the domain of interest. Nonetheless, the coastal regions of New England and the Mid-Atlantic show slightly negative correlations between −0.10 and −0.20, with limited areas showing significant values (Fig. 4b). In relative terms, these coastal areas received up to 15 %-20 % less precipitation during positive NAO compared to negative NAO months (Fig. 4a). When moving northward or westward from the northeastern US coast, the results show a gradual reversal of the correlation coefficient, which becomes slightly positive in southern Newfoundland and around the Great Lakes. In Atlantic Canada, the correlation coefficient is as high as 0.3 on the southern coast of Newfoundland. The correlation is instead moderately negative over Labrador, being the only area that shows a clear relationship between the NAO and mean precipitation (Fig. 4b). When looking at snowfall ( Fig. 4c and d), the relationship between snowfall and the NAO is much stronger than for total precipitation, especially for the northeastern US ( Fig. 4c  and d). In a wide corridor from southern Nova Scotia to the southeast US, a significant negative correlation is shown together with notable decrease in snowfall during positive NAO. In the Mid-Atlantic region, the snowfall that is received during months with positive NAO is on average half the snowfall during negative NAO months. The negative correlation between snowfall and the NAO is statistically significant all over the northeastern US except at the northernmost parts of the region, which implies that high snowfall is associated with the negative phase of the NAO. Moving northward into southeastern Canada, the ERA5 results show correlations fading toward zero and becoming weakly positive in most of the region between 45 and 50 • N. In Newfound-land and Labrador, the snowfall relationship pattern with the NAO is similar to that between precipitation and NAO, which is not surprising as almost all precipitation falls as snow during the cold season in this region (Fig. 1c).
The relationship between winter snowfall-to-precipitation ratio and the NAO (Fig. 4e and f) shows a very similar pattern in the northeastern US as the link between snowfall and NAO ( Fig. 4c and d). As for snowfall, southern New England and the Mid-Atlantic are the areas with the strongest correlation between NAO and S/P . Everywhere within the northeastern US, positive NAO tends to bring a lower S/P ratio, the only exception being northern Maine. Southern parts of eastern Canada also shows a decrease in S/P during positive NAO months. In southern Nova Scotia, this decrease is as high as 0.05 (i.e., from 45 % to 40 % of S/P ), and the negative correlation is significant. In Newfoundland and northeastern Quebec, the S/P ratio is instead higher during positive NAO.
All over the domain of interest, the ERA5 and stationbased results are consistent with each other. Only small differences exist between the results. Notably, for the correlation between snowfall and the NAO over southern Quebec, a slight negative correlation also exists with station results, although not with ERA5 results. Overall, these results indicate that the NAO has a much larger impact on the variability of winter precipitation and snowfall in the northeastern US compared to southeastern Canada.

Changes in heavy precipitation, snowfall and maximum snow depth
The relationship between the monthly NAO and heavy precipitation is fairly similar to the relationship between the NAO and total precipitation. However, the results for heavy precipitation are noisier, and there are more inconsistencies between reanalysis-based and station-based results, which is explained by the much larger variability in the amount of days per months with heavy precipitation compared to the variability of the mean precipitation. In the coastal regions of the northeastern US, there is a slight negative correlation that reaches statistical significance in some areas. There are on average 20 % to 30 % less heavy precipitation days during positive NAO months compared to negative NAO months in that region. The inverse results are seen over the southern half of Atlantic Canada and over the Great Lakes region, where positive NAO months see on average slightly more heavy precipitation days than negative NAO months. The relationship is again much stronger when we look at heavy snowfall (days with snowfall above 95th percentile) rather than heavy precipitation ( Fig. 5c and d). A significant negative correlation exists in a large area covering the entire Mid-Atlantic region, southern New England and southern Nova Scotia. The densely populated regions between Boston, MA, and Washington, DC, are found to be affected up to twice as much by heavy snowfall days during negative NAO months compared to positive months (Fig. 5c). Heavy snowfall events over eastern Canada are less affected by the NAO phase, with very limited areas that show significant correlations, except over eastern Labrador, where a negative correlation is seen between the two variables.
While the spatial distribution of the correlation coefficient between the NAO and SWE (Fig. 5e and f) is less homogeneous, it is coherent with the results for mean and heavy snowfall. This is not surprising as high snowfall is well correlated with maximum SWE in winter. Significant negative correlation is again shown in the Mid-Atlantic coastal region. As south of 40 • N snow that falls during storms usually stays on the ground only for few days, this also suggests that there is an increase in heavy snowstorms affecting the Mid-Atlantic during negative NAO winters. Although not statistically significant, negative NAO winters tend to be linked to a 15 %-30 % increase in snow cover depth in New England compared to positive NAO. In eastern Quebec, there is a significant positive correlation between SWE and NAO, which is also consistent with the positive correlation with snowfall and S/P .
As the snowpack is strongly affected by temperature, negative correlations and increase in maximum SWE in the northeastern US during negative NAO could also be partly explained by the colder conditions prevailing in that region during the negative phase (Fig. 6), which was also pointed out in previous studies (Notaro et al., 2006;Ning and Bradley, 2015). In a similar manner, colder temperature witnessed during a positive NAO phase in eastern Quebec ( Fig. 6; also in Wettstein and Mearns, 2002) can explain the positive correlation between SWE and the NAO found in that region.

Changes in storm tracks
In order to better understand the causes of the precipitation and snowfall anomalies discussed in the previous sections, we analyze the changes in track density and cyclogenesis across North America associated with the two phases of the NAO (Fig. 7).
During positive NAO conditions, the results show a considerable positive anomaly on lee cyclogenesis in Western Canada compared to neutral conditions (Fig. 7a), giving way to an increase in cyclonic activity in Canada and the Arctic during a positive NAO phase (Fig. 7b). However, over the US, there is no significant change in storm tracks, except for a slight decrease in lee cyclogenesis over Colorado (Colorado lows) and Montana.
The most notable result is the increase in cyclone occurrence in the western North Atlantic just offshore of the northeastern US coast during a negative phase of the NAO (Fig. 7d). While a decrease in cyclone occurrence is seen in that region during a positive phase (Fig. 7b), the anomaly is not as strong as during a negative phase. This implies that negative NAO conditions particularly favor an increase in Shaded circles represent station-based correlation coefficients that are statistically significant at the 0.05 level, and diamonds represent values that are not statistically significant. (c, d) As in panels (a) and (b) but for the relationship between monthly NAO and the number of heavy snowfall days (above 95th percentile of daily snowfall). (e, f) As in panels (a) and (b) but for the relationship between cold season averaged NAO and annual maximum snowpack depth. frequency of nor'easters (coastal storms). Near both Cape Hatteras and Cape Cod, which are two typical regions of coastal storm cyclogenesis on the US east coast (Davis et al., 1993), cyclogenesis is favored during negative NAO (Fig. 7c). Another important aspect of the changes in the storm track is the southwest-northeast-oriented dipole in the cyclone occurrence anomaly over the eastern half of the con-tinent during negative NAO months (Fig. 7d). This dipole suggests a southward shift in the storm track as the storm track density decreases over southeast Canada and the northeast US but increases over the southeast US.
The spatial patterns of correlation between the NAO and track density (Fig. 8) are consistent with the anomalies of cyclone occurrence for both positive and negative NAO. Large  areas of significant negative correlation are obtained over the southeast US and in the western North Atlantic between 35 and 45 • N, as well as over Labrador. Over Canada's Western and Arctic regions, the correlation is instead positive, with statistical significance that is reached over a large area.
The results also show that more cyclones tend to follow the typical North Atlantic storm track during positive phases, with more cyclones passing directly over Newfoundland and over the North Atlantic, south and east of Greenland during positive NAO phases.  Cyclogenesis and the storms displacement are directly related to the upper-air circulation. For that reason, we have also investigated the relationship between NAO and the mean 500 hPa geopotential height. The link between these two variables is more evident over the North Atlantic (Fig. 9), as it is closely related to the definition of the NAO. A dipole pattern of strong positive correlation in the subtropical and midlatitude North Atlantic and a negative correlation in the subpolar North Atlantic are clearly shown. This pattern extends over eastern North America. Over southeastern Canada and the northeastern US, a positive correlation is seen, indicating higher 500 hPa geopotential heights during positive NAO conditions. The difference in the average position of the 5400 m isohypse between positive and negative phases could indicate that troughs are more likely to be present over the eastern part of the continent during negative NAO.

Discussion
The results presented in the previous section provide insight on the mechanism of winter climate variability in relationship with the NAO. Changes in mean circulation patterns associated to the NAO over North America and the North Atlantic produce significant shifts in cyclogenesis and associated storm tracks, directly impacting regional precipitation and snowfall in the region.
The increase in coastal cyclogenesis near the US east coast during negative phases (Sect. 3.3) is likely caused by more frequent or deeper troughs over eastern North America, as seen in the difference in the average position of the 5400 m isohypse between positive and negative phases (Fig. 9). A region of divergence aloft and positive vorticity advection then more often overlaps the high temperature contrast region near the northern edge of the Gulf Stream offshore of the US east coast, leading to stronger surface cyclogenesis. While their cyclogenesis is favored during a negative NAO phase, the spatial pattern of anomalies shows that cyclones tend to take a more zonal and southern path in the Atlantic (Fig. 7d). Rogers (1990), Serreze et al. (1997) and Pinto et al. (2009) also observed that during negative NAO phases coastal storms often diverge from the continent near the 45th parallel and follow a track almost directly east, as opposed to continuing on a more typical northeastward trajectory following the coastline. This result suggests that they tend to affect the region in the vicinity of their formation area, with strong onshore flow and heaviest snowfall north of their path. But as they rapidly move offshore, their effects are not felt as much in areas further up the east coast. This could explain why the Mid-Atlantic region and New England have slight negative correlations between NAO and precipitation (mean and heavy) but not Atlantic Canada (Figs. 4b and 5b). The more zonal storm track of coastal storms during negative phases is likely caused by the much higher frequency of blockings near Iceland and Greenland, as presented in Fig. 10 (also in Shabbar et al., 2001;Woollings et al., 2008). The presence of a blocking high over Greenland forces the storms to move south of it, which is visible in the significant decrease in number of storms over Newfoundland and south of Greenland (Fig. 7d).
As the position of the rain/snow boundary during winter storms is strongly linked to the location relative to the storm center (e.g., Donaldson and Stewart, 1989), the southward shift in the storm track over the eastern half of the US explains the more frequent heavy snowfall event, higher snowfall and S/P witnessed all over the northeastern US and Nova Scotia during negative NAO months. This relationship found between snowfall, S/P and the NAO is consistent with previous studies that found that major snowstorms that affected the northeastern US are often associated with a negative NAO (Kocin and Uccellini, 2004;Hartley and Keables, 1998). Over Labrador, a lower (higher) cyclone occurrence during positive (negative) NAO relative to neutral (Fig. 7d) is consistent with the negative correlations found between precipitation and the NAO and between snowfall and the NAO. The higher occurrence of cyclones over Labrador and the western Labrador Sea during negative NAO may be also caused by the more frequent development of a Greenland blocking high during these NAO phases (Fig. 10). While most of the coastal storms will take a zonal path well south of the blocking high, some coastal storms, along with systems coming from the west, will be steered north, west of Greenland because of the anticyclonic circulation associated with the blocking high. Therefore, the results show a higher occurrence of cyclones in the western part of the Labrador Sea, bringing higher-than-normal precipitation to Labrador and southwestern Greenland (Fig. 4a), which was also pointed out in Auger et al. (2017). As they take a more northern track, these low-pressure systems advect warm air in the region, lowering the S/P .
In western Canada, a positive anomaly on lee cyclogenesis during positive phases (Fig. 8a) gives way to an increase in Alberta clippers forming over the plains and traveling eastward. Alberta clippers (or Canadian clippers) are cyclones that form frequently east (leeward) of the Canadian Rockies, over Alberta or the Canadian territories during the cold season. They are called clippers as they typically travel very rapidly over the continent (clippers were some of the fastest moving ships of the 19th century). Although they are not usually associated with very heavy precipitation, there is an increase in their occurrence during positive NAO over Canada and the northern US (Fig. 7a). This response could have a significant impact on the occurrence of storms in the Prairies region of Canada, such as an increase in prairie blizzard frequency. Furthermore, it explains the slight positive correlation between heavy precipitation and the NAO over the Great Lakes region (Fig. 5b), as these storms tend to move eastward and pass over the Great Lakes. However, these events are mostly in the form of rain as no correlation between heavy snowfall and the NAO is present in that region.

Summary and conclusion
In this study we examined the variability of the winter climate in relation with the NAO during the winter season over North America, with a particular focus on heavy precipitation and snowfall in the northeastern US and southeastern Canada. In order to better understand these changes in precipitation and snowfall, we also investigated the storm track and cyclogenesis variability associated with the NAO phases, as well as the changes in blocking frequency across the North Atlantic Basin.
First, we show that while there is only a slight negative correlation between precipitation and the NAO in the northeastern US, average and heavy snowfall are instead significantly affected by the NAO phase in that region. Both reanalysis and station-based results show that positive NAO phases tend to bring considerably less snowfall compared to negative NAO phases months over a wide region covering Nova Scotia, New England and the Mid-Atlantic of the United States. Henceforth, a significant negative correlation is also seen between the snowfall-to-precipitation ratio and the NAO in the same region. These results are explained by the increase in cyclogenesis of coastal storms near the US east coast during negative NAO phases, as well as a southward shift in the mean storm track over the United States. This increase in cyclogenesis is likely caused by more frequent or deeper troughs over the eastern US. During negative NAO phases, the more frequent blocking over Greenland also forces cyclones to follow a more southerly and zonal storm track in the North Atlantic, leading to fewer cyclones tracking directly over New England and Atlantic Canada. As a result, negative correlations between the NAO and precipitation, snowfall, and S/P in the coastal northeastern US slowly fade when moving northward 45 • N. Moreover, positive NAO phases are linked to a considerable increase in lee cyclogenesis east of the Canadian Rockies compared to neutral NAO conditions. As cyclones that form due to lee cyclogenesis usually travel eastward across the continent to affect southeastern Canada, a slight positive correlation is found between cyclone occurrence and the NAO over most of eastern Canada. This causes more frequent heavy precipitation events to occur over the Great Lakes region during positive NAO months. Even with the increase (decrease) of cyclone occurrence during positive (negative) NAO over southeastern Canada, the only regions where total precipitation is significantly positively correlated to the NAO phase are southern Ontario and limited parts of eastern Quebec and southern Newfoundland.
To conclude, using a combination of station-based and high-resolution ERA5 reanalysis data over a longer and more recent period  compared to previous studies, our study provides additional results on the mechanisms responsible for winter precipitations and storm track variability in the northeastern US and southeastern Canada. In particular, very few studies have hitherto investigated in detail the NAO impacts in Canada, and in particular on heavy precipitation and snowfall in this region. Given that global forecast models can predict reasonably well large-scale regimes such as the NAO at a lead time of several weeks (e.g., Johansson, 2007;Scaife et al., 2014;Black et al., 2017), the results presented in this study could prove valuable for improving the subseasonal forecasting of the probability of the occurrence of precipitation events, in particular heavy snowfall events.
Author contributions. JC performed the analysis and wrote the manuscript. FSRP contributed to the analysis design, interpretation of the results and writing of the manuscript.