Structure and Evolution of Non-Lake-Effect Snow Producing Alberta Clippers

: Alberta Clippers (clippers) have long been associated with lake-effect snow (LES) events due to their frequent passage over the Great Lakes basin. However, not all clippers produce LES, and no research has inquired into which synoptic ﬁelds most inﬂuence LES formation. This study analyzes clippers during non-LES situations to further knowledge on which atmospheric variables most regulate LES development on the synoptic scale. As no such database currently exists, a clipper repository is developed using National Centers for Environmental Prediction Reanalysis data. The repository is then cross referenced with a previously developed LES repository to identify clippers responsible for LES. Composite synoptic-scale patterns were then constructed on the remaining non-LES clippers to identify synoptic conditions that ultimately inhibited LES formation. This analysis is supplemented by an assessment of lake surface conditions in each composite to evaluate how inﬂuential the lake characteristics were in the suppression of LES activity. In total, 51 non-LES clippers were identiﬁed, tracked, and separated into three composite map types that exhibited unique storm track and spatial characteristics. Permutation testing revealed that lake surface conditions were not signiﬁcantly ( p ≤ 0.05) different between LES and non-LES associated clippers implying the main LES inhibition factors were meteorological.


Introduction
Mid-latitude cyclones (hereafter referred to as 'cyclones') have been studied frequently, as knowledge of the dynamics, structure, and tracks comprising these storms is critical for mitigation of their impacts. Cyclones form in areas of enhanced thermal contrasts as baroclinic instability leads to the transformation of potential to kinetic energy [1] (p. 98). These cyclones have a dominant role in influencing atmospheric processes on various spatial and temporal scales. For example, prior research has linked cyclones with severe tornado outbreaks [2,3], global poleward heat transport [4,5], regional seasonal precipitation distribution [6], and blizzard conditions [7]. Though cyclones can form in any region of enhanced baroclinicity, the complex topography of North America (elevated landmasses adjacent to major ocean currents, etc.) render distinct regions favorable for cyclogenesis [8], resulting in frequent cyclone development and propagation over the continental U.S. This is most pronounced during the North American cool season when the polar jet stream is strongest. Previous studies have identified three major cyclogenesis regions and associated storm tracks that have unique characteristics and impacts: the Colorado cyclone, East Coast cyclone, and Alberta Clipper [8][9][10][11][12][13][14].
One of the less studied aspects of cyclones is their indirect contribution to the formation of lake-effect snow (LES) off the North American Great Lakes (hereafter referred to as 'Great Lakes'). While LES occurs across the globe (e.g., Great Salt Lake (North America), Great Slave Lake (North America), Lake Baikal (Asia), Baltic Sea (Europe)) [15] (p. 17), Typically, clippers first evolve as an antecedent west coast cyclone makes landfall west/northwest of the Canadian Rocky Mountains [37]. While the cyclone propagates east towards Alberta, Canada, cross mountain flow increases, resulting in leeward adiabatic warming and the development of a lee trough [38,39]. As the cyclone continues propagating and traverses the Rockies, synoptic vertical ascent associated with quasigeostrophic (QG) omega forcing becomes prevalent atop of the low-level lee trough, resulting in contrasting vertical circulations characterized by low-level descent and mid-and upper-level ascent. Vortex stretching ensues which, when combined with enhanced baroclinicity via Typically, clippers first evolve as an antecedent west coast cyclone makes landfall west/northwest of the Canadian Rocky Mountains [37]. While the cyclone propagates east towards Alberta, Canada, cross mountain flow increases, resulting in leeward adiabatic warming and the development of a lee trough [38,39]. As the cyclone continues propagating and traverses the Rockies, synoptic vertical ascent associated with quasigeostrophic (QG) omega forcing becomes prevalent atop of the low-level lee trough, resulting in contrasting vertical circulations characterized by low-level descent and mid-and upper-level ascent. Vortex stretching ensues which, when combined with enhanced baroclinicity via the cross mountain zonal thermal contrast as well as differential positive vorticity advection (PVA), culminates in the development of a clipper on the lee side of the Canadian Rockies. Once developed, clippers track east/southeast ( Figure 1) where thermal advection and upperlevel forcing patterns are strongest [40].
Clippers are generally considered weak cyclones that feature high propagation speeds and rarely produce substantial precipitation [14,41] owing to a deficiency of moisture availability and lower thermodynamic energy availability. Clippers are also not typically associated with severe convective weather due to their relatively weak thermal gradients and QG forcing compared to other storm tracks [14]. However, clippers have been linked with strong wind events with gusts up to 38 m s −1 which can cause power outages, infrastructure damage, wildfires, and blizzard conditions [42,43]. Fall winds (also known as frontal chinooks) are frequently observed leeward of the Rockies as continental polar air associated with the clipper's cold front is advected across the mountain range, resulting in higher density air to flow and accelerate across the eastern slopes. Additionally, as the clipper propagates and matures, an anticyclone will typically build in its wake which creates a pressure dipole and tight isobaric gradient between the two pressure systems. As [35,36] note, in addition to strong flow, this dipole has also been linked with LES off the Great Lakes.
Previous studies have outlined the following conditions as typical for a synoptic environment comprising Great Lakes LES [35,36,[44][45][46]: • mid-tropospheric low-pressure anomaly situated near the Hudson Bay; • cyclone east/northeast of the Great Lakes (typically a clipper); • anticyclone west/southwest of the Great Lakes.
This surface dipole structure ( Figure 2) results in synoptic scale cold air advection (CAA) (linked with the cyclone and associated cold front) and LES conducive PBL wind regimes characterized by west/northwesterly flow. This pattern promotes air parcel trajectories over maximum amounts of open lake fetch. Synoptic scale forcing plays a secondary role in LES formation by modifying the low-level wind field which plays an instrumental role in LES morphology. That is, the low-level flow regime, which is dictated by the overlying synoptic environment, directly governs which types of snow bands (long-lake axis parallel, widespread coverage, etc. [20]) form during LES events. This was observed in [35] and [36] who found that, for a PBL thermodynamic profile conducive to LES formation, subtle variations in the surface dipole structure driven by deviations in the synoptic scale pattern led to varying LES snow band types and impacts.  From [35] highlighting the evolution of the common dipole structure present during a LES event off the eastern Great Lakes (Lakes Erie and Ontario).

Alberta Clipper Repository
The research objectives require a robust database of historic clippers. The work in [37] developed such a repository using European Centre for Medium-Range Weather Forecasts (ECMWF) Tropical Ocean Global Atmosphere (TOGA) data during 15 boreal cold seasons (October-March) spanning 1986/87 to 2000/01. They classified a cyclone as a clipper based on five propagation and intensity criteria based on mean sea level pressure (MSLP) and vorticity fields (see [37] for a full description). Our study utilized a similar approach, adapting these criteria to National Centers for Environmental Prediction (NCEP) National Center for Atmospheric Research (NCAR) Reanalysis data [47] that spans 11 meteorological winters (December-February) including 1997/98, 1998/99, and Of course, not all clippers result in this dipole structure, and even if the dipole is formed, LES still may not develop. The exact reasons for this are not known, despite advances in understanding the synoptic-scale's role in driving LES. Understanding the unique synoptic-scale dynamics of non-LES associated cyclones will supplement our current knowledge of LES-conducive synoptic-scale patterns and assist local forecasters with predicting LES with greater lead time. As such, the purpose of our study is to examine the synoptic scale characteristics of wintertime clippers that do not result in eventual Great Lakes LES production and differentiate those patterns from traditional LES synoptic-scale configurations. The resulting patterns will be compared with their underlying mesoscale and lake-surface conditions to ensure differences between LES and non-LES environments are primarily meteorological. Section 2 describes methods employed for developing a clipper repository, reanalysis datasets used, and statistical analysis techniques for creating the synoptic composites. Section 3 presents a preliminary analysis of the clipper repository and a comprehensive examination of the synoptic composites and provides a discussion on the influence of the synoptic environment on LES suppression by analyzing the mesoscale environment. Finally, Section 4 summarizes the results and discusses future work.

Alberta Clipper Repository
The research objectives require a robust database of historic clippers. The work in [37] developed such a repository using European Centre for Medium-Range Weather Forecasts (ECMWF) Tropical Ocean Global Atmosphere (TOGA) data during 15 boreal cold seasons (October-March) spanning 1986/87 to 2000/01. They classified a cyclone as a clipper based on five propagation and intensity criteria based on mean sea level pressure (MSLP) and vorticity fields (see [37] for a full description). Our study utilized a similar approach, adapting these criteria to National Centers for Environmental Prediction (NCEP) National Center for Atmospheric Research (NCAR) Reanalysis data [47] that spans 11 meteorological winters (December-February) including 1997/98, 1998/99, and 2005/06 to 2013/2014. These specific time periods were selected to align with LES repositories in [35] and [36] which were used to categorize clippers as associated/disassociated with LES. Meteorological winter (December-February) was selected as the temporal domain for this study owing to the climatological frequency of clippers [37], and only 48 (16.5%) LES cases occurred outside the winter months. The NCEP/NCAR reanalysis dataset was selected to develop the clipper repository as it features an equal spatial (2.5 • ) and higher temporal (6 h) resolution to that of the ECMWF TOGA data used by [37] and has been used in previous work [14] characterizing clipper environments. All cyclone identification and tracking was performed on a spatial domain spanning 25 • N-75 • N and 150 • W-30 • W, a domain seen in previous studies [8,14,37] when characterizing clipper climatology. For a cyclone/local MSLP minimum of interest to be classified a clipper, the following conditions must be met (note these were based heavily on the work in [37]):

1.
Cyclogenesis in the lee of the Canadian Rocky Mountains within the Canadian provinces of Yukon, Northwest Territories, British Columbia, and/or Alberta as evidenced by a closed circulation and/or evident local MSLP minimum (based on 2 mb intervals); 2.
Closed circulation and/or evident local MSLP minimum present for 60 h once the system makes major progress east of the Canadian Rocky Mountains (hereafter referred to as 'time of departure'); 3.
System propagates southeast towards U.S-Canadian border immediately after time of departure whereafter it tracks east and/or southeast; 4.
System is located east of 90 • W within 60 h after time of departure (a measure of cyclone propagation speed).
Based on these criteria, 78 total clippers were identified over the 11 meteorological winters previously mentioned (Table 1). On average, 7.1 clippers occurred every winter, consistent with [37] who found that approximately 7.2 clippers occur per winter. Interestingly, the intraseasonal distribution differed slightly with [37], as most clippers identified in our study occurred in January, while their work found most clippers occurred in December. These differences were attributed to the small temporal span in both studies. Once identified, clippers were visually tracked using the Grid Analysis and Display System [48] and ArcGIS Pro version 2.5.0 [49]. These mapping and visualization tools were used in tandem to locate and record the position of the lowest MSLP from the NCEP/NCAR data, noted as the clipper's center, an approach consistent with numerous prior cyclone-tracking studies with similar objectives [37,[50][51][52]. Finally, the 'end time' for each clipper was identified as the NCEP/NCAR timestep where the MSLP minima was no longer present in the domain as it crossed 50 • W or 60 • N. This methodology resulted in a clipper repository that contained the following cyclone characteristics:

•
Time of departure date and location; • End time date and location; • Complete track of system throughout its duration as marked by NCEP/NCAR Reanalysis grid points at six-hour timesteps.
Once finalized, each clipper system in the repository was categorized a LES or non-LES producing clipper by cross-referencing the clipper repository with the LES repository from [35,36]. For a clipper to be linked to LES formation, it must have influenced the surface flow regime over the associated lake during the LES event start time. This influence was determined visually by examining all LES events that started within the life span of a given clipper and comparing the MSLP setup with the surface wind field over the lake. If that wind field was predominantly a result of the clipper system under investigation, the clipper was logged as being 'LES associated'. This approach yielded 21 (26.9%) of the 78 total clippers being LES associated, with the remaining 57 cases non-LES associated.

Data and Statistical Methods
Once established, the non-LES clipper database (N = 57 cases) was used to construct a synoptic climatology using a T-mode principal component analysis (PCA) combined with a k-means cluster analysis (CA) to obtain climatological map types. These methods have been shown to be useful in numerous studies [3,35,36,44] as they identify temporal modes of variability among atmospheric data and use that information to group individual members (i.e., clippers for this study) into distinct clusters from which composites can be constructed.
Though the NCEP/NCAR reanalysis was optimal for developing the clipper repository (owing to its reliance on the synoptic scale), its spatial and temporal resolution were not ideal for representing the mesoscale conditions associated with the clippers. Instead, the North American Regional Reanalysis (NARR) dataset [53] was used to construct the synoptic climatology of the non-LES clippers as it features a higher spatial (32 km) and temporal (3 h) resolution than the NCEP/NCAR reanalysis while encompassing the study period and spatial domain. In total, 41 meteorological variables were retained from the NARR to characterize the synoptic and mesoscale conditions within each clipper, including five surface fields (mean sea-level pressure (MSLP), 10-m zonal and meridional wind components, skin temperature, and specific humidity) and five three-dimensional fields (zonal and meridional (u and v) wind components, geopotential height, temperature (T), and specific humidity (q)) captured at seven isobaric levels (1000 mb, 925 mb, 850 mb, 700 mb, 500 mb, 300 mb, and 250 mb). All fields were retained for 96 h, beginning with the respective clipper's time of departure over a spatial domain that extended from 25 • N-65 • N and 130 • W-40 • W and included 30,352 NARR gridpoints (a near match to the domain used for the clipper repository).
Importantly, in this study, it was unknown which timestep(s) of the NARR would best characterize the distinctions between LES and non-LES environments, as each clipper had different evolutionary characteristics (cyclogenesis location, propagation speed, etc.). As no LES occurred in the non-LES clippers, we estimated the geographic position of the clippers where LES would most likely occur as the locations where these distinctions should be found. This was achieved by computing the mode longitude of the 21 LES associated clippers (as discussed above) at the time LES began. This longitude (75 • W) marked the location where clipper associated LES was most likely to form, though this longitudinal distribution was clearly bimodal (Figure 3), an important limitation of this approach. We addressed this issue by selecting NARR fields from two longitudinal positions in the clipper's lifespan (75 • W and 90 • W) to serve as an objective spatial and temporal framework from which non-LES associated clippers could be compared with their LES counterparts. Note that six non-LES and two LES clippers (roughly 10%) did not track east of 75 • W and fell outside of the bimodal distribution. These systems were excluded from further analysis, as they did not characterize the overlying patterns of the other clippers (resulting in 51 and 19 total non-LES and LES clippers, respectively).
Once the non-LES fields were retained, synoptic-scale composites of the dominant patterns in these environments were constructed following [35,36]. The composite methodology required several tuning/optimization steps before final map types could be constructed. These steps included determining which NARR fields best distinguish the clipper environments, the number of principal components (PCs) to retain from the PCA, and number of clusters/map types from the CA. These steps were completed by conducting numerous PCAs and CAs on 11 combinations of NARR fields representing the 51 non-LES cases ( Table 2) and quantifying variance explained by the resulting PCs, as well as silhouette coefficients [54] associated with each CA. through the progression of the experiments. This suggests that MSLP was the most distinguishing atmospheric field among the NARR fields tested, which parallels the findings of [35,36] who found that MSLP was the least correlated atmospheric field among synoptic composites of LES events. It should be noted that MSLP was the only NARR field to be tested individually, owing to its characterization of the dipole structure thought to be of critical importance for LES formation in previous work [35,36,[44][45][46]. With the optimal PCA-CA configuration identified, a nonhierarchical k-means CA was used to separate the 51 non-LES clippers into three distinct clusters ( Figure 4) based on variability structures identified by the PCA. Clippers in each cluster were averaged to construct to three sets of synoptic composites that depicted atmospheric conditions for all clippers in each group (map types) at each reference longitude (75° W and 90° W). Finally, a set of mean composites for the 19 LES clippers were constructed as a reference to compare against the non-LES patterns derived from the composite analysis described above.  We employed methods from [36] to determine the optimal PCA/CA configuration and tested all permutations of retaining 2-10 PCs and 2-10 clusters for each NARR configuration (i.e., 'Experiment'; Table 2). The PC-cluster combination selected to construct the composites was based on the amount of variance explained by the PCs, the average silhouette coefficient, and the frequency of misclustered cases (negative silhouettes). Among the 11 experiments, 'Experiment 1', which consisted of only retaining MSLP, two PCs, and three clusters, was deemed best, as its PCs explained the most variance (25.9%) and featured the second highest silhouette coefficient (0.442) with only one misclustered clipper (Table 3).  Figure 4 shows the distinctiveness of the PC loadings, revealing that Cluster 1 is generally associated with PC1 loading values that are positive, while Cluster 2 had positive PC2 and negative PC1 values and Cluster 3 had negative PC1 and PC2 values. Interestingly, the addition of more NARR variables appeared to deteriorate analysis performance as metrics used to evaluate the PCA and CA generally decreased through the progression of the experiments. This suggests that MSLP was the most distinguishing atmospheric field among the NARR fields tested, which parallels the findings of [35,36] who found that MSLP was the least correlated atmospheric field among synoptic composites of LES events. It should be noted that MSLP was the only NARR field to be tested individually, owing to its characterization of the dipole structure thought to be of critical importance for LES formation in previous work [35,36,[44][45][46]. Table 3. Statistical results from experiments used to identify the optimal PCA and CA configuration. Each of the listed values from each experiment represent the most optimal PC-cluster configuration associated with the NARR field(s) retained for that specific experiment. Starred fields represent the configuration used for the statistical analysis.

Diagnostic Variables
Following [35,36], MSLP and upper-level geopotential height fields were used to diagnose synoptic-scale structure and forcing patterns. Additionally, a derived quasi-geo- With the optimal PCA-CA configuration identified, a nonhierarchical k-means CA was used to separate the 51 non-LES clippers into three distinct clusters (Figure 4) based on variability structures identified by the PCA. Clippers in each cluster were averaged to construct to three sets of synoptic composites that depicted atmospheric conditions for all clippers in each group (map types) at each reference longitude (75 • W and 90 • W). Finally, a set of mean composites for the 19 LES clippers were constructed as a reference to compare against the non-LES patterns derived from the composite analysis described above.

Diagnostic Variables
Following [35,36], MSLP and upper-level geopotential height fields were used to diagnose synoptic-scale structure and forcing patterns. Additionally, a derived quasigeostrophic (QG) variable was calculated to assess synoptic-scale vertical motion. When assessing synoptic-scale vertical motion, using the traditional QG omega diagnostic approach can prove difficult in situations when differential geostrophic vorticity advection and temperature advection counter one another, yielding indeterminate vertical motion insight even though such motion may be present. This issue was present in our analysis (not shown), so we elected to utilize a derived QG diagnostic that blends both terms in the QG omega equation by coupling geostrophic horizontal shear with the horizontal temperature gradient on an isobaric surface, a quantity known as the Q-vector [55]. Q is directly related to QG omega through: where Q is defined as: This relationship shows that areas with Q-vector convergence (divergence) are colocated with synoptic-scale ascent (descent). Following the methods of [14], static stability (σ) was excluded from the Q calculations as it can be divided out as a scalar without altering the direction of Q (as σ is almost always positive for large-scale synoptic analysis).
In addition to the synoptic-scale analysis, a mesoscale analysis was completed which characterized the role of surface-atmosphere stability and lapse rates in LES suppression. Low-level (1000-850 mb) lapse rates were calculated over a NARR grid point ( Figure 5) centered over each lake (resulting in five lapse rates for five lakes) to evaluate stability. These lake-centric grid points were selected as they feature the highest lake surface temperatures due to the lakes' bathymetry patterns and are co-located the location of where LES associated convection would be most likely to develop initially. Finally, surface specific humidity (q) fields were evaluated to assess atmospheric moisture content.
To ensure the LES suppression mechanisms were meteorological, lake surface conditions were also analyzed separately given their importance on LES development. Specifically, if stark differences in the lake surface temperatures and lake ice cover arose between LES and non-LES clippers, this would suggest lake conditions were the primary factors differentiating LES and non-LES cases. Lake temperature data were retained from the daily Great Lakes Surface Environmental Analysis (GLSEA) Surface Water Temperature Data archive [56], while lake ice cover was based on the GLSEA Great Lakes Average Ice Cover Data [56] which features daily lake average ice cover. It should be noted that the ice cover dataset began December 2008 which resulted in a portion of the clippers (42.9%) being excluded from the ice cover analysis.
LES and non-LES clippers, this would suggest lake conditions were the primary factors differentiating LES and non-LES cases. Lake temperature data were retained from the daily Great Lakes Surface Environmental Analysis (GLSEA) Surface Water Temperature Data archive [56], while lake ice cover was based on the GLSEA Great Lakes Average Ice Cover Data [56] which features daily lake average ice cover. It should be noted that the ice cover dataset began December 2008 which resulted in a portion of the clippers (42.9%) being excluded from the ice cover analysis.

Influence of Lake Surface Characteristics on LES Suppression
To quantify the role of lake surface conditions on LES suppression, Table 4 presents average lake surface temperatures and ice covers of all LES and non-LES clippers across all Great Lakes. Permutation tests of the mean [57] (p. 182) were used to establish statistical significance to the differences in LES and non-LES lake temperature and ice covers at each reference longitude.

Influence of Lake Surface Characteristics on LES Suppression
To quantify the role of lake surface conditions on LES suppression, Table 4 presents average lake surface temperatures and ice covers of all LES and non-LES clippers across all Great Lakes. Permutation tests of the mean [57] (p. 182) were used to establish statistical significance to the differences in LES and non-LES lake temperature and ice covers at each reference longitude. Though warmer surface temperatures (roughly 0.25 • C) and lower ice covers (roughly 3.5%) were observed with LES clippers, which are indictive of a more LES conducive environment, these differences were not statistically significant for any lake (Table 4). Across all lakes and reference longitudes, ice cover differed more than lake surface temperatures, as evidenced by smaller p-values (not shown). The lack of statistical contrast implies that the primary forcing mechanisms suppressing convective activity associated with non-LES clippers were not primarily based on lake conditions, but instead on the mesoscale and synoptic-scale environment.

Synoptic Analysis
Clippers were mostly evenly distributed among the three resulting clusters (N1 = 18, N2 = 19, and N3 = 14). While storm characteristics varied among all clippers analyzed, each composite featured distinct underlying structures ( Table 5) in terms of average storm speed, duration, and intensity. Cluster 1's composites were characterized by typical clipper propagation characteristics [37] highlighted by a closed low-pressure circulation in the Canadian Rockies lee (53 • N, 112.5 • W) that tracked southeast through the first 24 h before taking on a more westerly track clipping the U.S.-Canadian border ( Figure 5). Cluster 1 featured the strongest non-LES systems characterized by the lowest central pressure values, slowest propagation speeds, longest average duration, and were most comparable to LES systems of all clusters (Table 5). However, Cluster 1 systems generally tracked further north than LES clippers without as much meridional variation altering the surface flow regime (not shown). Composite NARR fields exemplified these characteristics with a 500 mb trough present at both reference longitudes collocated with a surface low-pressure system (Figures 6a and 7a). The initial westward displacement of the upper-level feature from the surface cyclone at 97.5 • W (not shown) suggests the system was in the 'open wave' stage of development based on the Norwegian Cyclone Model [58].
Upper-level forcing in Cluster 1's composite (Figure 6a) was weaker relative to LES systems (Figure 6d) as only a small area of Q convergence was observed over the southcentral coast of Lake Superior at 90 • W (Figure 6a). This was likely due to weak (0-5 m s −1 ) southerly flow resulting in weak warm air advection (WAA) that enhanced rising motion across the western Great Lakes region (Figure 7a). Though the flow characteristics were overall similar between Cluster 1 and the LES composite (not shown), the cross-basin near-surface temperature gradient was exceptionally higher for LES systems, resulting in heightened WAA and Q-vector convergence ( Figure 7d). As the upper-level feature continued to strengthen and propagate east, rising motion was observed over the eastern Great Lakes basin as evidenced by a pocket of large Q-vector convergence over the eastern lee of Lake Ontario (Figure 8a). However, the magnitude of this forcing decreased due to the increase in CAA as the clipper propagated east and surface winds acquired a northerly component (not shown). Additionally, Cluster 1 clipper's westward displacement from the upper-level feature decreased, resulting in eventual weakening of the system, and marking its transition into the mature phase. This was evidenced by MSLP tendencies as the Cluster 1 clipper's central pressure was 8 mb higher at 75 • W than 97.5 • W (Figure 8a). In contrast, the LES clipper's central pressure was 4 mb lower at 75 • W than 97.5 • W. Through Cluster 2's composite clipper progression, an upper-level vorticity maximum developed north of the Great Lakes basin as the trough-ridge pattern damped, resulting in minimal Q-vector convergence at the location when LES was most likely to form (Figure 8b). This pattern contrasted LES systems that strengthened throughout their progression. As the clipper exited the Great Lakes basin (Figure 9b), an anticyclone originating from western Canada propagated southeastward, roughly following the Cluster 2 composite clipper. This resulted in the traditional high-low pressure dipole structure coupled with large-scale CAA over the north central U.S, a pattern commonly seen in previous studies [35,36] during LES episodes (as well as in the LES composites). However, the absence of upper-level forcing and the relatively stable environment over the lakes (further discussed below) suppressed convective activity. Note that the strength of the gradient between the dipole structure was higher for LES systems as well, featuring stronger high-(1030 mb) and low-pressure (1008 mb) systems which produced faster winds (5-10 m s −1 ). This suggests that the intensity of the dipole structure may indirectly be a differentiating factor between LES and non-LES clippers. The synoptic structure and propagation of Cluster 3 notably differed from the first two clusters and most matched the LES composite, even though its intensity characteristics most differed. Similar to the LES composite, Cluster 3's storm track featured meridional variation absent from Clusters 1 and 2 as it originated at the northernmost location (54.6° N) and followed the southernmost track ( Figure 5). Cluster 3 clippers propagated tario) LES conducive environment as the southwest-northeast pressure gradient resulted in southwesterly flow across a large fetch across the two lakes. This contrasts the LES dipole that featured a purely zonal pressure gradient leading to westerly winds (not shown) across most of the Great Lakes. However, upper-level forcing was minimalized through Cluster 3′s progression due to strong CAA (Figure 9c) and, as in Cluster 2, the strength of the dipole was weaker than the LES composites which generated weaker flow (0-5 m s −1 ) (not shown).  Cluster 2 composites followed a similar storm track to Cluster 1, though the overall track position was further north than LES clippers ( Figure 5). Cluster 2 clippers were on average much less intense (6.3 mb higher central MSLP) than LES systems and Cluster 1 and featured shorter lifespans and faster propagation speeds (Table 5). This was partially attributed to Cluster 2 featuring the highest percentage (47.4%) of December cases when the polar jet stream is weaker compared to January due to the lessened meridional equatorpole temperature gradient [1]. Composites showed an initial upper-level trough and associated surface cyclone that advanced zonally, strengthening through its progression. Synoptic-scale ascent was maximized over the same general region as Cluster 1 (south central Lake Superior coast) based on Q-vector fields while the clipper was located at 90 • W ( Figure 6b). As in Cluster 1, this strengthening was attributed to WAA associated with southerly flow across the Great Lakes basin while the cyclone propagated into the study area (Figure 7b). However, the southerly flow and associated WAA inhibited convective development, making LES formation unlikely contrasted to the westerly flow observed with the LES composite (not shown).
Through Cluster 2's composite clipper progression, an upper-level vorticity maximum developed north of the Great Lakes basin as the trough-ridge pattern damped, resulting in minimal Q-vector convergence at the location when LES was most likely to form (Figure 8b). This pattern contrasted LES systems that strengthened throughout their progression. As the clipper exited the Great Lakes basin (Figure 9b), an anticyclone originating from western Canada propagated southeastward, roughly following the Cluster 2 composite clipper. This resulted in the traditional high-low pressure dipole structure coupled with large-scale CAA over the north central U.S, a pattern commonly seen in previous studies [35,36] during LES episodes (as well as in the LES composites). However, the absence of upper-level forcing and the relatively stable environment over the lakes (further discussed below) suppressed convective activity. Note that the strength of the gradient between the dipole structure was higher for LES systems as well, featuring stronger high-(1030 mb) and low-pressure (1008 mb) systems which produced faster winds (5-10 m s −1 ). This suggests that the intensity of the dipole structure may indirectly be a differentiating factor between LES and non-LES clippers.

Mesoscale and Stability Analysis
As most of the differences in the LES and non-LES composites were in the structure and magnitude of temperature advection, a mesoscale analysis focusing on stability and temperature advection was completed to complement the synoptic analysis. For Cluster 1, WAA patterns enhanced atmospheric stability at 90° W, which suppressed LES across the study region (Figure 7a) by modifying static stability. Relatively high inland surface temperatures combined with the low-level WAA produced the lowest low-level lapse rates (3.5 °C km −1 ) from Cluster 1's environment of all composites considered (and much lower than the LES composite of 6.7 °C km −1 ). The flow pattern from the WAA was also unsupportive of LES, as numerous studies [18][19][20][21][22] have observed that north/northwesterly flow, not southerly flow, accommodates LES formation by ensuring maximum fetch across each lake (except Lake Michigan, due to its meridional orientation). Later in the clipper's track, the surface wind pattern veered significantly, resulting in westerly flow as the dominant regime when the clipper was centered 75° W (Figure 9a). This pattern would be most similar to a typical LES setup as more of the necessary ingredients were observed (e.g., CAA near the backend of the cyclone's life cycle and westerly flow). However, 1000 mb temperatures over the southern half of the Great Lakes basin were above freezing throughout the clippers progression which, combined with enhanced stability and minimal upper-level forcing, may be why Cluster 1 clippers did not result in LES. These high temperatures also resulted in relatively high atmospheric moisture content with specific humidity values ranging from 3.5-4.5 g kg −1 (Figures 7a and 9a). Lastly, the horizontal The synoptic structure and propagation of Cluster 3 notably differed from the first two clusters and most matched the LES composite, even though its intensity characteristics most differed. Similar to the LES composite, Cluster 3's storm track featured meridional variation absent from Clusters 1 and 2 as it originated at the northernmost location (54.6 • N) and followed the southernmost track ( Figure 5). Cluster 3 clippers propagated faster than Clusters 1 and 2, and were weakest based on central MSLP ( Table 5). The initial composite pattern featured an upper-level local vorticity maximum over northern North Dakota that was coupled with a synoptic scale 500 mb low pressure anomaly centered over the eastern Hudson Bay, and an associated surface low-pressure system over north Minnesota (not shown). As the surface clipper system progressed southeast, the vorticity maximum strengthened initially, owing to height falls, resulting in the development of a large-scale positively tilted upper-level trough at 90 • W (Figure 6c). Simultaneously, strong Q-vector convergence was present over Michigan's lower peninsula via strong low-level WAA (Figure 7c), shifting the upper-level trough to a neutral phase as the system propagated east. Q-vector convergence was maximized over western Lake Superior adjacent to the backend of the cyclone where WAA was strongest. Overall, synoptical vertical forcing was strongest of all clusters and closely matched the LES composite.
As the clipper tracked eastward, the lessened displacement of the upper-level vorticity maximum from the surface cyclone reduced baroclinicity and resulted in the clipper weakening while still tracking 1 • latitude south of the Cluster 1 and 2 tracks. Simultaneously, an anticyclone originating from the Rocky Mountains began building in the clipper's wake from considerable CAA and anticyclonic vorticity advection (AVA) (not shown), resulting in the surface dipole structure observed in Cluster 2 and the LES composite. The orientation of the dipole suggests an eastern Great Lakes (Lakes Erie and Ontario) LES conducive environment as the southwest-northeast pressure gradient resulted in southwesterly flow across a large fetch across the two lakes. This contrasts the LES dipole that featured a purely zonal pressure gradient leading to westerly winds (not shown) across most of the Great Lakes. However, upper-level forcing was minimalized through Cluster 3 s progression due to strong CAA (Figure 9c) and, as in Cluster 2, the strength of the dipole was weaker than the LES composites which generated weaker flow (0-5 m s −1 ) (not shown).

Mesoscale and Stability Analysis
As most of the differences in the LES and non-LES composites were in the structure and magnitude of temperature advection, a mesoscale analysis focusing on stability and temperature advection was completed to complement the synoptic analysis. For Cluster 1, WAA patterns enhanced atmospheric stability at 90 • W, which suppressed LES across the study region ( Figure 7a) by modifying static stability. Relatively high inland surface temperatures combined with the low-level WAA produced the lowest low-level lapse rates (3.5 • C km −1 ) from Cluster 1's environment of all composites considered (and much lower than the LES composite of 6.7 • C km −1 ). The flow pattern from the WAA was also unsupportive of LES, as numerous studies [18][19][20][21][22] have observed that north/northwesterly flow, not southerly flow, accommodates LES formation by ensuring maximum fetch across each lake (except Lake Michigan, due to its meridional orientation). Later in the clipper's track, the surface wind pattern veered significantly, resulting in westerly flow as the dominant regime when the clipper was centered 75 • W (Figure 9a). This pattern would be most similar to a typical LES setup as more of the necessary ingredients were observed (e.g., CAA near the backend of the cyclone's life cycle and westerly flow). However, 1000 mb temperatures over the southern half of the Great Lakes basin were above freezing throughout the clippers progression which, combined with enhanced stability and minimal upper-level forcing, may be why Cluster 1 clippers did not result in LES. These high temperatures also resulted in relatively high atmospheric moisture content with specific humidity values ranging from 3.5-4.5 g kg −1 (Figures 7a and 9a). Lastly, the horizontal temperature gradient over the region was minimal compared to the LES composite, resulting in weaker CAA even though the flow regimes were similar.
Most mesoscale structures evident in Cluster 1 were also present in the Cluster 2 composites. Weak WAA owing to a minimal temperature gradient was present over the study region as the clipper entered the basin, and this advection quickly transitioned to CAA as winds acquired a westerly component in the wake of the cyclone (Figures 7b and 9b). Inland 1000 mb temperatures initially ranged from 0 • C to −4 • C across most the basin before decreasing in the northern basin via synoptic scale CAA, supporting snowfall formation (Figure 9b). This CAA occurred while the composite clipper tracked through 75 • W, as an anticyclone propagated southeast setting up the previously mentioned dipole structure. Once again, the flow regime and CAA present at 75 • W denotes where the environment was most conducive to LES generation. The Cluster 2 environment was less stable compared to Cluster 1, as indicated by higher low-level lapse rates (4.37 • C km −1 on average). However, this was still notably less than that of the LES composite, again suggesting high atmospheric stability was a major factor in LES suppression for Cluster 2. Additionally, similar to Cluster 1, near surface temperatures co-located with the strongest upper-level synoptic-scale support were again at or above freezing, inhibiting snow formation despite ample atmospheric moisture ( Figure 9b).
As observed with the synoptic fields, the mesoscale and surface characteristics of Cluster 3 differed the most from the other clusters. Near surface temperatures were exceptionally low over the Great Lakes basin, fluctuating between 0 • C off the southern Lake Erie and Ontario coasts to −10 • C over northern Lake Superior (Figure 7c), matching much more closely with what was observed with the LES composite. Southwesterly flow (not shown) occurred as the clipper entered the Great Lakes basin with considerable WAA observed east of Lake Superior. This continued as the clipper tracked east of 90 • W, resulting in inland temperatures warming 2-3 • C (Figure 7c). However, once the clipper crossed the Great Lakes basin, surface winds veered to the north, resulting in strong CAA due to the strong temperature gradient observed across the basin (Figure 9c). Upon the clipper's exit from the domain, surface temperatures fell dramatically, with some areas south of Lake Superior reaching −14 • C (Figure 9c). In general, the Cluster 3 mesoscale environment was most conducive to LES formation. Cluster 3 featured the steepest average lapse rates (5.91 • C km −1 ) and lowest inland surface temperatures with the entire Great Lakes basin featuring below freezing temperatures through the clippers progression.
To quantify statistical significance among the static stability differences in the LES and non-LES composites, permutation tests were run on lapse rates calculated for each storm in the LES and non-LES clipper repositories. To avoid an asymmetry of case distribution between these two clipper sets (19 LES clippers vs. 51 non-LES clippers), non-LES clippers were separated into their previously assigned clusters, leading to three sets of permutation tests run across each reference longitude. p-values calculated from the permutation tests showed that stability was a significant (p ≤ 0.05) differentiating factor between Cluster 1 clippers and LES clippers across all lakes (Table 6). This was also observed over Lakes Superior and Michigan for Cluster 2. This result is not surprising given the large inland near surface temperature contrast between these clusters and the LES composite. Similarly, the results between the Cluster 3 and LES composite lapse rates were not significantly different, which was expected based on results presented above. Interestingly, this suggests that the overlying mesoscale moisture profile was not important with regard to LES development, consistent with [21], which noted high levels of atmospheric moisture are not necessarily a mandatory requirement for LES formation and even small amounts of water vapor can result in substantial snowfall.

Summary and Future Work
This study's research objective was to analyze the spatial and temporal characteristics linked with non-LES associated clippers through the development of synoptic composites from a newly defined repository of clipper systems. Once clippers were identified and tracked, each was associated/disassociated with LES formation by cross referencing a previously developed LES repository. In total, 78 clippers were identified from the tracking methods, of which 57 were non-LES and 21 were LES systems. Composites were then constructed by blending a PCA with a CA that grouped non-LES clipper systems with similar atmospheric and propagation characteristics at two reference longitudes (75 • W and 90 • W) where LES was most likely to occur. Three non-LES clusters were selected based on PCA and CA evaluation metrics. Composites of the mean LES clippers were also constructed to evaluate potential separating atmospheric characteristics crucial to LES development that differed from the non-LES composites. Additionally, daily lake surface temperatures and ice cover data were retained to characterize the lake surface for both LES and non-LES clippers. Permutation tests revealed that the overlying lake environment did not significantly (p ≤ 0.05) contrast during LES and non-LES clippers, indicating the primary differences were meteorological.
Cluster 1 featured the strongest clippers (lowest average MSLP) and propagated close to parallel with the climatological clipper track as defined by [37]. Cluster 2 exhibited a similar track to Cluster 1, though its composite average strength was weaker (5.9 mb higher central MSLP) which can be attributed to their contrast in intraseasonal distribution to the rest of the clipper repository (most Cluster 2 cases occurred in December). The Cluster 3 composite mirrored Cluster 1 in terms of intraseasonal distribution (most clippers in January), but its track had the greatest latitudinal variation and featured the weakest systems on average. This southward shift in these clippers likely contributed to their weaker magnitude (relative to Cluster 1) as weaker planetary vorticity was present along most of the Cluster 3 track.
Analyses revealed several differentiating characteristics among the clusters and LES composite. Low-level stability appeared to be a primary distinguishing factor for Cluster 1 as low-level lapse rates were significantly (p ≤ 0.05) lower than LES clippers across all the Great Lakes. Furthermore, above freezing inland surface temperatures and weak CAA also factored into why LES activity was suppressed. As [21,59] notes, ample CAA is needed to create large vertical temperature gradients which generate turbulent heat and moisture fluxes and destabilize the overlying polar air mass. Cluster 1 mesoscale patterns were observed, albeit to a lesser degree, with Cluster 2. Additionally, Lakes Superior and Michigan had significantly lower lapse rates in their non-LES environments, while the other three lakes, while not statistically significant, did possess large differences in static stability when compared to the LES cases. Upper-level forcing appeared to be another distinguishing factor between these clusters (1 and 2) and LES systems. Analysis of Q-vector divergence fields indicated large amounts of synoptic ascent across most of the Great Lakes basin early in the LES composite clipper's life cycle that was absent for Clusters 1 and 2.
These trends were not observed with Cluster 3, as the mesoscale and thermodynamic environment was similar to the LES composite. The notable differences with Cluster 3's composite were synoptic in nature and largely based on the strength of the surface dipole. LES clippers were much stronger (21.8 mb lower central MSLP) than Cluster 3 systems on average, as was the anticyclone that built in the clippers wake. The subsequently strong winds in the LES composite ultimately enhanced low-level vertical heat and moisture fluxes via turbulent mixing between the lake surface and low-level atmosphere, distinguishing these patterns from the Cluster 3 composite.
As previous research has focused on the synoptic environment during LES events, the purpose of this research was to provide a baseline diagnosis of the synoptic conditions during non-LES situations associated with cyclonic systems that most frequently result in LES (i.e., clippers). These differences primarily included the presence and magnitude of synoptic forcing present, low-level stability, and the strength of the surface dipole. Future research will further investigate these meteorological traits through the development of a diagnostic objective classification model that categorizes LES and non-LES clippers based on results from this study. Reference [59] demonstrated that the climatological spatial snowfall patterns over Lake Michigan contain enough of a synoptic signal to objectively classify LES from synoptically driven snowfall. The authors plan to further this work by developing a machine learning based classifier using the results of this work. Optimizing the classifier will provide insight into which spatial scales and atmospheric fields are most important regarding LES development/suppression related to clippers. An analysis of surface temperature fields of all 19 LES and 51 non-LES cases revealed that the differentiating atmospheric fields separating these two systems goes beyond whether temperatures were above freezing. Knowledge of these physical traits will aid local forecasters and provide the foundation for future prognostic efforts.