Why do strong winter winds over Turkey play an important role in the formation of precipitation and snow in relation to synoptic patterns? Can it be shown using cluster analysis? What is the role of blocking in snowfall in mountainous regions? Does the Scandinavian block over the Black Sea play a role in heavy precipitation in Turkey?

We investigate the patterns of strong winds and the correlation between synoptic patterns in Turkey during the winter season using cluster analysis technique. We used the k-means algorithm to detect surface patterns of strong winds. In addition, we used the self-organizing map (SOM) technique to identify clusters of geopotential height at the level of 500 hPa, mean temperature at the level of 850 hPa, and mean sea level pressure. We used the dataset from the New European Wind Atlas (NEWA) project to analyze surface weather conditions and the ERA5 dataset to study upper-level weather conditions. The k-means algorithm, when applied to the land surface data in Turkey, identifies six distinct clusters. These clusters are mainly located around the Taurus Mountains, which extend in an east-west and northeast direction along the Black Sea coast. The formation of these clusters is controlled by the terrain features and their physical characteristics. The higher-level clusters, consisting of nine SOM nodes, are not affected by the terrain and weather systems that are characteristic of the Mediterranean macroclimate. These clusters are identified in the Eastern Mediterranean, Black Sea and inner Aegean regions, emphasizing the influence of topography on surface patterns.The assets that modern societies have developed are severely threatened by the strong winds and violent storms of today [1]. Studies on climatology and mid-latitude cyclones have grown in recent years to evaluate severe storms and strong winds. In order to represent this link, global or regional climate models are used. These models attempt to capture the relationship between the frequency and intensity of cyclones and extreme winds. According to research [1,2], severe winds are occurring more frequently in Europe now than they did at the end of the twenty-first century. Numerous studies have been carried out in the existing body of literature to develop models for climatological variables, including temperature, precipitation, and strong or extreme wind patterns. Most of the scholarly work in the field of study employs statistical methodologies, including but not limited to Extreme Value Theory, Generalized Pareto Distribution (GPD), Annual Maxima, and Peak Over Threshold. In addition, statistical decomposition techniques, such as Principal Component Analysis (PCA) and Empirical Orthogonal Function (EOF), are frequently employed in climatological analyses. Furthermore, in the field of meteorological analysis, clustering methods like Ward’s and k-means, as well as artificial network approaches like Self-Organizing Map (SOM) are also commonly utilized machine learning methods. The main objective of all of these techniques is to classify data (observations or variables) into defined subgroups with similar attributes. Using meteorological stations or gridded data and considering characteristics like temperature, wind, and precipitation, this division can yield equally distributed observations. However, an accurate evaluation requires precise data aggregation so that the aggregation of ensemble members helps assess uncertainty and model performance. As a prerequisite for statistical models to give accurate results, long-term data helps classical statistical models to yield better results. Extreme Value Theory and GPD analysis are proven statistical models utilized in climatological studies [3–7]. Advanced statistical approaches, such as cluster analysis, PCA, and EOF, are employed for the purposes of exploring data and reducing dimensionality. Their primary focus lies on comprehending the structure and patterns within large quantities of data, often employed for analyzing severe phenomena such as wind events. EOF analysis [8] is an additional statistical approach employed in climate studies. It aims to identify structures within a two-dimensional dataset that account for the greatest amount of variation, like PCA [9]. It is also employed to decrease the size of a dataset containing multiple variables into a smaller dataset with fewer variables, hence lowering variance. One dimension reflects the structure to be discovered, whereas the other does not represent its realization size. For instance, when studying spatial structures that vary over time, time is the structure size and space is the sample dimension. EOFs are first-dimensional structures from the analysis. The variability of wind direction statistics for both average and extreme wind occurrences [10], as well as surface wind patterns [11,12], has been evaluated by the application of EOF analysis in several studies. Aside from traditional statistical techniques, clustering method k-means, as well as artificial network approaches like SOM, are commonly employed in the study of meteorological data. Cluster analysis utilizes a similarity matrix to group climatic time-series datasets that have similar characteristics. These categories can be further analyzed in detail to get insight from the shared properties of the datasets in each climatic sub-region group. Cluster analysis in climatology or meteorology aims to establish distinct categories of synoptic types or climate regimes, as well as to identify climatic subregions by grouping stations and/or grid points [13]. The examination of clustering methods is typically categorized into two groups: hierarchical and non-hierarchical procedures. As a non-hierarchical clustering method, the k-means algorithm is widely acknowledged as a prominent clustering approach and is extensively employed in academic research. Several extensions of the k-means algorithm have been presented in academic literature. The k-means clustering algorithm, developed by Wong and Hartigan in 1979 [14] is a method of vector quantization that facilitates the partitioning of data into a predefined number of clusters. The k-means approach and its extensions, commonly used in unsupervised learning for clustering in pattern recognition and machine learning, are consistently impacted by the initializations of clusters that are predetermined. According to Sinaga and Yang [15], it may be inferred that the k-means algorithm does not fall under the category of unsupervised clustering methodologies. The implementation of the k-means approach is quite straightforward, and it is capable of handling large-scale data collections. The proposed algorithm ensures convergence and exhibits generalizability to clusters of varying forms and sizes, including elliptical clusters. There exist numerous studies that employ the k-means clustering on wind event analysis due to its simplicity and linear nature. The study conducted by Pena et al. [16] examines the correlation between strong wind events and atmospheric circulation patterns in Catalonia through the utilization of k-means clustering. The subject of Bernardino and his friends’ discussion was the categorization of synoptic and local-scale wind patterns in a coastal area of the Tyrrhenian Sea (Italy) using the k-means clustering method [17]. Another study implemented the k-means algorithm to cluster the multivariate wind speed data from the 75 Turkish provinces. The clustering was based on monthly mean wind speed datasets [18]. Burlando utilized k-means cluster analysis to evaluate the synoptic-scale surface wind climate patterns in the Mediterranean Sea [19]. The study employed k-means and EOF analysis to evaluate the occurrence of severe windstorms over Europe in both current and future climate conditions [20]. On the other hand, the SOM [21] is a categorization technique commonly employed in climate research. Unlike k-means clustering, SOMs can effectively reveal complex and non-linear relationships within high-dimensional data spaces. They utilize a neural network algorithm and an unsupervised learning approach to discern patterns within datasets. The method generates a two-dimensional array of categories that represent the probability density function of the input data. The array contains a greater number of categories in regions where the data density is higher. The ultimate allocation organizes comparable attributes in proximity, disparate ones at a further distance, and gradual transitional categories in the intervening space [22]. Cavazos employed the SOM technique to uncover significant intra-seasonal changes in the North American monsoon. This is achieved by analyzing various atmospheric variables daily, including 850 hPa meridional winds, 700 hPa specific humidity, 500 hPa geopotential heights, and 850–500 hPa thickness [23]. The SOM is employed for the purpose of categorizing the variables of the mid-troposphere, namely the air temperature, geopotential height, and specific humidity at the 700 hPa level. This classification is conducted to identify generalized atmospheric patterns. Additionally, the SOM is utilized to recreate the synoptic patterns of climate in the Antarctic region, based on data obtained from ice cores [24]. The classification of meteorological station data using the SOM technique is discussed in another study to evaluate many factors, including temperature, humidity, wind, daylight hours, and solar radiation, among others [25]. In their study, Khedairia and Khadir conducted a classification analysis on meteorological data from the Annaba region, located in the North-East of Algeria. They employed SOM and k-means clustering methods for their research [26]. The origin of wind gusts in Australia was investigated using the SOM approach by Spassiani et al., focusing on meteorological factors [27]. In another study, Zhao et al. utilized SOM to identify the specific characteristics of largescale atmospheric circulation patterns that contribute to extreme wind events in Beijing, China during the spring and winter seasons [28]. In a study conducted by Kim et al., the researchers examined the categorization of strong wind events and gusty winds in Korea. They utilized the SOM technique to identify and classify locations with similar wind characteristics [29].In Türkiye, there has been specialized research conducted on the spatial dispersion of temperature, humidity, and precipitation, as well as their correlation with synoptic and large-scale wave patterns. In a study done by Türke¸s et al. [30], k-means and hierarchical clustering were used to review Türkiye’s climatic zones and compare historical periods. The dataset contained 1950–2010 air temperature and precipitation. Western, central, and southern regions that have a dry summer subtropical Mediterranean climate experienced less rain than northern and eastern regions. Kömü¸scü et al. [31] investigated the spatiotemporal variability of precipitation clusters in Türkiye through the application of k-means clustering. During the time span from 1980 to 2020, the analysis identified five distinct groups, each consisting of four to eight sub-periods lasting either 5 or 10 years. The Arctic Oscillation has a stronger association with the Aegean and Mediterranean regions due to the presence of clustering structures and teleconnection patterns. The research indicates that precipitation patterns changed throughout time. Another study examined unusual air movements at the 500 hPa geopotential level during periods of severe North Atlantic Oscillation (NAO) Index in Türkiye. The findings indicate that the average yearly, winter, spring, autumn, and partially summer composite precipitation averages are generally more humid compared to the long-term average circumstances during the negative phase of the NAO Index. On the other hand, positive NAO Index responses are associated with drier conditions throughout the year and in all seasons except for summer. Noticeable variations in precipitation levels were particularly evident in the western and central regions of Türkiye. The NAO pattern exerted influence on the weather conditions in the North Atlantic and Europe, resulting in less precipitation during the positive phase of the NAO Index [32]. Studies on extreme winds and associated synoptic patterns, aside from variables like temperature and precipitation, are quite rare in research done especially for Türkiye. In one study, mesoscale models and climatic forecasts were used to investigate how climate change may affect wind energy resources in Türkiye, a Mediterranean hotspot. The findings indicate that Türkiye’s wind resources are expected to be impacted by climate change, with summertime increases and wintertime declines. During the winter, wind speeds in the Aegean Sea, Aegean Region, Marmara Region, and Marmara Sea are expected to decrease, but may increase at other times. The energy production predictions for two wind farms were also computed using mesoscale model data [33]. Tornadoes, which are abrupt natural calamities, have resulted in substantial damage in Türkiye in recent times, as indicated by another study. A synoptic pattern classification was devised to categorize tornado occurrences from 2000 to 2020. This classification is based on factors such as the circulation type at the 500 hPa level, the placement of troughs and ridge axes, and surface characteristics. The study analyzed a total of 520 tornado occurrences and 408 days, categorizing them into three distinct types: “Sharp Trough”, “Spanning Trough”, and “Closed System” [34]. Türkiye possesses a highly complex geographical terrain. Depending on their dynamic and physical characteristics as well as the trajectories they follow, synoptic and largescale atmospheric systems that influence Türkiye interact with this complex topographic structure in different ways, resulting in a structure that changes both seasonally and regionally. The topography and atmospheric systems have a significant impact on winds, which is an important meteorological parameter. Türkiye has seasonal and local variations in its wind patterns as a result of its coastline being bordered by the sea on three sides and the existence of elevated plateaus in its interior regions. This study examines the relationship between the wind regime in Türkiye and synoptic and large-scale wave patterns, with a focus on the winter season. The study specifically investigates how the systems that are effective on the synoptic scale, particularly during winter, impact the wind regime in Türkiye. The wind regime in Türkiye is classified based on this consideration. This study aims to present the initial efforts to analyze extreme winds over Türkiye using k-means clustering on ground-level data and SOM clustering on upper-level data, respectively. The k-means clustering approach was chosen due to its simplicity and ability to differentiate similar patterns. [35]. Moreover, the SOM algorithm provides an efficient means of objective classification by classifying and visualizing synoptic states in the time series of atmospheric state variables, utilizing input data to discover archetypical states [22]. Accordingly, this study analyzed the association between the clustering results at the lower and upper levels and investigated the impact of upper-level patterns on extreme wind speeds at the lower level. The paper is organized into three main sections, excluding this introductory section. In Section 2, the two datasets are introduced: 100 m above ground level for wind speed cluster analysis and a set parameter for upper-level clustering. The data selection criteria are also explained. Valuable insights regarding the wind speed data are given before proceeding to the clustering process. Next, the clustering process at 100 m above ground level and synoptic level are described in detail, and cluster validation metrics are discussed. In Section 3, both k-means and SOM clustering results are evaluated in synoptic terms.The winter subtropical high-pressure systems in the Northern Hemisphere are centered over the oceans because of the spatial expansion and strengthening of the thermic high pressures in Siberia and North America (Canada/Greenland). These high-pressure systems exhibit a regional weakening because of strong cooling in large continental areas in the north. However, dynamic semicontinuous and cold low-pressure systems at about 60 degrees north latitude in the North Atlantic (Iceland) and North Pacific (Aleutian) successfully extend their influence on general flows over the ocean between the southern 30–35 degrees north latitudes. These systems are located on the North Pacific and North Atlantic. Because of this, there is a high likelihood of intense precipitation and storms occurring over regions where these pressure systems are effective [35]. In winter, extratropical cyclones that originate from Iceland and their related frontal systems move towards southern latitudes, impacting the mid-latitudes and playing a crucial role in the development of strong storms. Mediterranean frontal systems and their associated cyclones, which originate in the Mediterranean and travel eastward, exert an influence on Türkiye by following different trajectories. The initial orbit follows a south-west to north-east trajectory, while the second orbit runs from west to east along the southern Aegean and Mediterranean coasts of Türkiye. The cyclone’s travel along both circles can result in powerful southwesterly and southerly winds on its front (east side), as well as strong northerly and northeasterly winds and severe weather events on its back (west side) [35,40]. The upper-level variables were filtered so that the winter days for the 30-year period were included in the final dataset. The final dimensions for the reduced dataset for each variable are 2707 × 121 × 201, where 2707 is the number of days within the selected period, and 121 × 201 is the total number of grid points for the selected region. The final dataset has temporal alignment, so that each row corresponds to a winter day (2707 samples) and each column (grid point) corresponds to a feature (24,321 features). Unit transformations have also been applied to obtain the actual values for each variable. Data normalization or standardization is a common preprocessing step for SOM that guarantees scale homogeneity [58]. Methods such as Z-score normalization or Min-Max scaling can be incorporated into this procedure. In order to normalize the input, Min-Max scaling has been applied to the reshaped dataset. Prior to SOM clustering, the EOF analysis was applied to the min-max normalized datasets separately for dimensional reduction purposes. As EOF per variable is applied independently, the aim is to avoid the potential issue of one variable’s magnitude dominating the combined EOF analysis due to differences in units or variability. The threshold for cumulative variance for the mean sea-level pressure, geopotential height, and temperature is set to 90%, 90%, and 85%, respectively. Based on the EOF analysis 8, 8, and 9, a total of 25 PCs have been determined as the optimal number of PCs that explained the desired level of cumulative variance as 91.19%, 91.54%, and 86% per upper-level variable, mean sea-level pressure, geopotential height, and temperature, respectively (Figure 6).Results 3.1. k-Means Clustering Results The k-means clustering successfully characterized the wind dispersion patterns in the Aegean Region, considering the influence of the terrain. Nevertheless, the impact of the Black Sea and Mediterranean Region Mountain ranges on Central Anatolia was adequately represented by the k-means clustering algorithm with a 6-cluster test (Figure 7a). Also, the average values of all clusters were illustrated in Figure 7b. The resulting graph was used to highlight the cluster grid points in terms of contour plots.The area along the Black Sea coast. The Taurus Mountains and the Black Sea Mountains are both situated in parallel to the sea. A contiguous arrangement of clusters may be observed along a linear trajectory spanning from the Eastern Mediterranean to the Eastern Black Sea. A distinct pattern has emerged along a linear trajectory extending from the Black Sea Region to the eastern part of Marmara. An alternative cluster arrangement is evident, originating from the Central Black Sea Region and extending to the southern part of Central Anatolia. Comparable characteristics are also observed in the inner Aegean region. Another recurring pattern is observed throughout the Aegean coast. The Mediterranean and Black Sea exhibit distinct structural characteristics. It is clear from looking at the various clustering setups that the terrain affects the surface patterns. 3.2. EOF Analysis and SOM Clustering Results The experiments conducted for upper-level cluster analysis were examined with a different approach. The EOF analysis was independently applied to the min-max normalized datasets for the purpose of dimensional reduction. The application of EOF to each variable is done separately to prevent the possible problem of one variable’s magnitude overpowering the combined EOF analysis as a result of the variations in units or variability. Additionally, Figure 8 displays the findings of the first three EOFs, illustrating the presence of different patterns. Atmosphere 2024, 15, x FOR PEER REVIEW 14 of 26 m8. The first three Empirical Orthogonal Functions (EOFs) during DJF for a period of 30 years (1989–2018): (a) MSLP; (b) 500 hPa geopotential height; (c) 850 hPa temperature. The third EOFs, which consist of GPH with an explained variance (EV) of 12%, MSLP with an EV of 17%, and T with an EV of 11%, are most likely to represent the less persistent, more localized, or transitional characteristics of the atmospheric state. A strong gradient is observed in MSLP over Türkiye, which can be attributed to distinct atmospheric patterns present across Europe, East Africa, and Arabia. Similarly, in the EOF3 GPH analysis, a strong gradient was observed along the west-east axis across Türkiye. Figure 8. The first three Empirical Orthogonal Functions (EOFs) during DJF for a period of 30 years (1989–2018): (a) MSLP; (b) 500 hPa geopotential height; (c) 850 hPa temperature. On the other hand, in typical circumstances, EOF results are used to discern the main patterns of variability in upper-level data, particularly in the identification of synoptic patterns. GPH, a measure of the actual height of a pressure level in the atmosphere, is essential for comprehending synoptic weather patterns. An EOF examination of GPH at the 500 hPa level may identify prevailing wave patterns, such as ridges and troughs, which govern the movement of weather systems [59]. The following EOF modes could show secondary patterns, such as localized weather occurrences or smaller-scale characteristics. Less dominant but still significant patterns, such as a cutoff low, which can be associated with strong winds and stormy conditions, can be detected via the third EOF mode. Furthermore, MSLP, which plays a crucial role in explaining surface meteorological conditions, displays influence over wind patterns, storm courses, and frontal systems [60]. The first EOF mode often aligns with the principal synoptic pattern, including factors such as the location and intensity of high and low-pressure systems. Successive EOF modes might potentially detect additional characteristics, such as frontal systems, pressure gradients, or localized pressure anomalies. More localized pressure anomalies that drive mesoscale wind events, which can be particularly extreme, can be captured by the third EOF mode. In addition, EOF analysis can detect significant temperature anomaly patterns linked to occurrences like heatwaves, cold spells, or the location of air masses in the temperature study [61]. The initial EOF mode can uncover large-scale temperature variations, such as the distribution of warm and cold air masses or the existence of temperature anomalies. Sequential EOF modes have the potential to capture smaller-scale temperature fluctuations, such as boundaries between air masses or differences in temperature inside specific areas. The illustration of the first EOFs of each synoptic variable explains the largest fraction of the variance in the dataset. The three leading EOFs of the data, as shown in Figure 8, account for 70.5%, 63.95%, and 58.94% of the total variance for the MSLP, Z500, and T850 variables, respectively. In Figure 8, the initial EOFs of GPH and T, which collectively account for 27% of the variance, along with the first EOF of MSLP, which accounts for 35% of the variance, are likely the most significant synoptic-scale patterns for the region. The majority of the EOF1 MSLP visuals exhibit a dominant pattern, which is consistently observed across the whole region encompassing Türkiye, the Mediterranean, and western Italy. The areas indicated by positive values in the figure, which show an opposite pattern, are also visible in the northern part of the Caspian Sea and the Gulf of Basra. The EOF1 GPH exhibits a pronounced gradient in Türkiye, which arises from differing values in the East Anatolia region compared to the northwestern part of Europe. This scenario depicts strong winds, particularly in the western areas. In Türkiye, the EOF1 T exhibits a prominent pattern at its center. The second EOFs, which consist of 24% GPH, 19% MSLP, and 21% T, represent secondary yet significant patterns. The EOF2 MSLP exhibits pronounced gradients in Türkiye, with opposite values in northern and northeastern Europe and in Mediterranean and northern Africa. The EOF2 GPH is characterized by a unique pattern centered over Greece and Italy, and a contrasting pattern in Arabia. Türkiye is located within a strong gradient that originates from the west. A significant temperature gradient is developing over Türkiye in EOF2 T due to the positive influence of different patterns in Europe and Arabia. These variations may be attributed to the irregularities in the placement of the ridges and troughs (GPH), the distinct pressure centers (MSLP), and the related thermal gradients (T). The interaction of these patterns can result in fluctuations in both the speed and direction of wind, which can give rise to rare cases of extreme wind events. The third EOFs, which consist of GPH with an explained variance (EV) of 12%, MSLP with an EV of 17%, and T with an EV of 11%, are most likely to represent the less persistent, more localized, or transitional characteristics of the atmospheric state. A strong gradient is observed in MSLP over Türkiye, which can be attributed to distinct atmospheric patterns present across Europe, East Africa, and Arabia. Similarly, in the EOF3 GPH analysis, a strong gradient was observed along the west-east axis across Türkiye. Figure 8 illustrates those different patterns at the Z500, MSLP, and T850 levels in the north-west, south-east, north-east, and south-west regions of Türkiye, which contribute to the formation of a strong gradient area across the country. In synoptic climatology, EOFs decompose a dataset into orthogonal (independent) modes of variability, which helps in understanding the spatial patterns and their temporal evolution. However, SOMs are considered as a complementary tool to classify largescale atmospheric patterns, like those represented by the selected PCs from EOF analysis.The SOM arranges many synoptic scenarios into a grid, such as a 3 × 3 grid in our empirical setup. Each node in the grid indicates either a typical pattern or a group of similar patterns. The utilization of this grid arrangement enables a straightforward and precise representation and examination of intricate atmospheric dynamics. EOFs/PCs and SOMs together offer a potent toolkit for recognizing, categorizing, and displaying atmospheric patterns in synoptic climatology. EOFs provide a linear and orthogonal breakdown of the data, but SOMs offer a non-linear and more adaptable approach to classify and analyze the data. SOMs can capture intricate and nuanced correlations in the data, beyond the limitations of EOF which impose orthogonality constraints on the patterns. The Principal Components (PCs) generated from the EOF analysis are used as input for the SOM. This ensures that the SOM training is concentrated on the most important patterns of variability discovered by the EOF analysis so that SOM helps in interpreting large-scale drivers of local climate variability and is a part of statistical downscaling methodologies [62]. Once the SOM cluster analysis was completed, the days that belong to each SOM node, referred as SOM neurons, were extracted. Each SOM node was visualized based on the mean GPH, MSLP, and T as contour plots together (Figures 9–13). A comprehensive analysis was conducted to see how these nodes may accurately represent the cumulative variability observed in the EOFs. This section primarily examines how the SOM nodes can exhibit different synoptic-scale characteristics that may impact extreme surface winds.The significant pressure gradients appear to create regions of strong winds throughout Türkiye. The occurrence of low-pressure and high-pressure areas, namely in the Mediterranean and the Black Sea, coupled with the corresponding frontal systems, greatly impact the frequency of days marked by extreme winds. Furthermore, the Low-Pressure Centre, primarily located over the Black Sea, exerts a significant impact on wind intensification due to its frontal effect in the area. Lastly, the Azores and Siberia, places known for their high-pressure systems, play a significant role in creating substantial gradients where southerly and southwesterly winds dominate over Türkiye. As a result of the strong pressure gradient, areas with extreme winds are formed, particularly in the western Mediterranean, Aegean, Marmara, and Western Karadeniz regions. Within the Central and Eastern Black Sea region, the intensification of the gradient results in the formation of strong wind zones caused by westerly flows along these coastal areas. Corroborating this, there exists a significant temperature gradient at the 850 hPa level. The upper-level wave pattern shows that, the upper-level trough, linked to the low-pressure center over Italy, extends its influence on the western coasts of Türkiye. These regions are situated within the ascending flow area of the trough, resulting in a stronger impact on increasing the winds in this area. The ridge area exerts its effect over the eastern and southeastern areas. This refers to the decreased frequency of windy days in this area. The southerly airflow originating from the Mediterranean exerts an important accelerating influence, particularly on the southwestern Aegean and western Mediterranean. This dominant southerly flow pattern aligns with the high frequency of strong winds experienced, particularly in the western and northwestern regions of Türkiye. Generally, there exists a significant relationship between the frequency of days with high wind speeds and the strong pressure and temperature gradients generated by various ground-level pressure systems (Figure 9a). • SOM Node (C2): Observations indicate that a dynamic low-pressure system, originating from Iceland, has extended its area of influence towards southern Europe. There is a thermal high-pressure system located in the eastern region of Türkiye. The prevailing wind direction in Türkiye is generally from the southwest, while westerly winds are noted around the Central Black Sea coast and the Eastern Black Sea region. The southeastern region of Türkiye experiences easterly winds because of the influence of the thermal high-pressure system. The strong pressure gradient observed in the western and northwestern areas, as well as along the Black Sea coast, correlates with the high number of days characterized by extreme winds in these specific locations. The strong temperature gradient at the 850 hPa level along the Eastern Black Sea coast also influences the number of days with high winds in this area. Similarly, when analyzing the wave pattern at higher altitudes, a strong gradient in geopotential height is observed over the western and northwestern parts of Türkiye. This, in turn, leads to a strong zonal flow at the 500 hPa level, which has a direct impact on the increase of surface winds. (Figure 9b)SOM Node (C3): A high-pressure center is located over the cold central European continent. Furthermore, there exists a cold-origin Siberian high-pressure center located on the east of the Caspian Sea. The low-pressure system centered on the Red Sea affects Türkiye. The relationship between the high-pressure system situated over Europe and the low-pressure system over the Red Sea is linked to the number of days with extreme winds over the Aegean region. The temperature field at 850 hPa shows that the temperature gradient in the eastern Mediterranean region aligns with the pressure gradient at the surface level. The impact of the trough on Türkiye is observed at the 500 hPa level. The central and eastern Mediterranean Region, characterized by significant pressure and temperature gradients, is within the region of ascending flow. This accounts for the increased number of days characterized by high winds in the area (Figure 10a). • SOM Node (C4): It is observed in this pattern that Türkiye is being affected by a low-pressure center originating from the Mediterranean, with its center located in the eastern Mediterranean region. The pressure gradient is strong in the northwestern region of Türkiye and the eastern Mediterranean. Over Türkiye, a strong temperature gradient is seen at the 850 hPa level. The upper-level wave trough is situated to the west of the ground-based low-pressure center in this pattern, which illustrates the highly baroclinic unstable structure of a frontal system over the area. The number of days with high winds along the central and eastern Mediterranean coast is influenced by the rising flow impact of the upper-level trough moving in a south-west to northeast direction (Figure 10b). • SOM Node (C5): The Azores-originating high-pressure center exerts its influence over a broad region encompassing Türkiye, as well as western and central Europe. A low-pressure center is situated above Cyprus. The impact of a cold low-pressure system originating from Iceland extends into the northeastern region of the Black Sea in Eastern Europe. It is seen that a strong pressure and temperature gradient formed along Türkiye’s Mediterranean coast in a west-east direction when combined with an examination of the temperature at 850 hPa and the pressure at ground level. Hence, it may be asserted that a frontal system is effective in the area. Also, it is seen that there is a pressure gradient throughout the eastern Black Sea. Examining the 500 hPa wave pattern reveals that the trough connected to a cold-originating low-pressure system over Eastern Europe has extended its impact to Türkiye. The Black Sea region is within the expanding zone of the trough’s flow. The convergence of various ground and upper-level patterns is linked to areas characterized by a significant frequency of days with extreme winds (Figure 11a). • SOM Node (C6): This node shows similar patterns with the C5 node. However, at this point, it is observed that the Azores high-pressure center has a more stable form and extends its range of influence to northern Europe. The flow pattern over Türkiye, which was predominantly northwestern and western in the preceding pattern, is therefore classified as northern in this instance. Both the ground-level pressure field and the 850 hPa level temperature field show strong pressure and temperature gradients in areas where the number of extreme windy days in this pattern is high. When looking at the geopotential height values, it becomes evident that the highpressure area over Europe corresponds to a ridge that impacts the western part of Türkiye. Meanwhile, the middle and eastern portions of Türkiye are influenced by a trough. The trough is mostly affecting the Mediterranean and Southeastern Anatolia regions, since they are located inside its ascending flow area (Figure 11b). • SOM Node (C7): As the low-pressure center affects northern Europe and the Western Mediterranean, the high-pressure center originating from Siberia is observed to extend its sphere of effect towards eastern and central Europe. While Türkiye has predominantly high atmospheric pressure, there is a localized region of modest low pressure over the Eastern Black Sea. The MSLP and 850 hPa temperature indicate the presence of a strong gradient aligned in the west-east direction, running parallel to the shoreline in the Mediterranean area. Consistent with this arrangement, the number of days with strong winds is greater along this trajectory. The occurrence of extreme winds, particularly in the eastern part of the Marmara Sea and the Inner Aegean Region, may be attributed to two factors: the presence of a low-pressure system near the surface, which creates a restricted track into the Black Sea, and the temperature gradient at the 850 hpa level. The upper-level wave pattern typically exhibits a ridge effect in regions that encompass Türkiye (Figure 12a). • SOM Node (C8): This pattern involves the influence of a high-pressure center that originates from the Azores and its accompanying ridge at the 500 hPa level. It impacts a large area encompassing central Europe, Türkiye, and northern Africa. The western areas of Türkiye are situated in a col, whereas the Anatolian plateau experiences a high-pressure center. A low-pressure system, which originated from Iceland, together with its associated 500 hPa trough, is now observable across eastern and northern Europe. The locations that have a high number of days with strong winds are found to be consistent with areas with strong pressure and temperature gradients, particularly those that stretch from the central Mediterranean to the southwest-northeast direction. Similarly, there exists a gradient of pressure over the central region, the Black Sea, and Georgia. (Figure 12b).• SOM Node (C9): The Mediterranean origin low-pressure center, which is centered over Italy, is shown to have strong southern flows that affect a large portion of Türkiye. A regional high-pressure center is present in the southeastern part of Türkiye, while a low-pressure center is located over the Caspian Sea. On the low-pressure center at ground level, there is a trough at 500 hPa, but it gives way to a ridge toward Türkiye’s east. There exists a strong pressure gradient at ground level in Türkiye, apart from the southeastern area. Particularly during periods when a low-pressure system dominates Türkiye, with prevailing winds originating from the southwest and south Mediterranean, there are strong upward movements associated with the low-pressure system. Additionally, southerly winds over the open sea intensify the flow, resulting in generally strong winds observed throughout Türkiye when this weather pattern is in effect. As seen by the observed pattern, there is a rising trend in the number of days characterized by extreme winds. Nevertheless, it is observed that the southeastern area has a limited number of days characterized by extreme winds due to the dominance of high pressure (Figure 13). The significant pressure gradients appear to create regions of strong winds throughout Türkiye. The occurrence of low-pressure and high-pressure areas, namely in the Mediterranean and the Black Sea, coupled with the corresponding frontal systems, greatly impact the frequency of days marked by extreme winds. Furthermore, the Low-Pressure Centre, primarily located over the Black Sea, exerts a significant impact on wind intensification due to its frontal effect in the area. Lastly, the Azores and Siberia, places known for their high-pressure systems, play a significant role in creating substantial gradients over Türkiye. These gradients are essential for the development of extreme winds. To summarize, it has been observed that most of the SOM nodes were representative enough to identify the common pressure systems which are effective throughout Türkiye during the winter period. Moreover, the vast majority of the extreme wind events were explained by the SOM nodes. Conclusions : In this study, the relationship between the ground-level wind speeds and the upperlevel parameters was investigated. The main objective was to accurately characterize the extreme wind conditions during the wintertime in Türkiye and examine their correlation with synoptic patterns. Cluster analysis techniques were employed to identify the correlation between them. Two distinct atmospheric levels—the upper level and the ground level—were considered during the calculations. The lower-level data used in the first section of the analysis was taken from the NEWA dataset, which includes 48 samples daily at 6732 grid points over a 30-year period at 100 m. The 98th percentile of the data was obtained for each grid point, and these extracted values were considered as the extreme values for that grid point. The k-means clustering algorithm was utilized to classify the given datasets. After doing a k-means clustering analysis, six clusters were identified. These clusters show key patterns along the Taurus and Black Sea mountains, demonstrating the importance of topography in meteorology. The complicated connection between terrain and wind patterns is shown by the linear trajectories of clusters from the Eastern Mediterranean to the Eastern Black Sea and from the Black Sea Region to Eastern Marmara, as well as comparable patterns in the Central Black Sea and inner Aegean. These findings give a sophisticated picture of regional wind dynamics, enabling meteorological and environmental investigations in the area. Also, it was discovered that the clustering results match the widespread presence of high or strong winds throughout Türkiye. Moreover, the analysis of the dataset revealed that the winter season in the study area exhibits the most frequent occurrence of windy conditions characterized by high wind velocities.The second phase of the study consisted of using upper-level datasets to investigate synoptic patterns. The analysis was performed using the MSLP, Z500, and T850 datasets. The data source for this analysis was the ERA5 Reanalysis data, covering DJF of a 30-year period (1998–2018). The combined application of EOF and SOM provided a detailed and multifaceted view of the atmospheric dynamics affecting Türkiye, particularly in terms of synoptic-scale wind patterns. The EOF technique was individually applied to the datasets to reduce their dimensionality while preserving the variance. The scree test indicated that the eight Z500, MSLP, and nine T850 EOFs gathered were considered adequate. Furthermore, the initial three findings of the EOF analysis were presented just to demonstrate the presence of different patterns. The SOM technique was implemented on the MSLP, T850, and Z500 datasets, which were obtained from the EOF results. A total of nine SOM nodes were generated, following a structured grid arrangement with dimensions of 3 × 3. Each node is assigned a unique set of days for a given atmospheric state. Three primary patterns were identified for the winter season because of these analyses, despite some variations among them. The first pattern is seen to be particularly efficient in patterns with low-pressure centers located over the Mediterranean region or Italy and extending towards Türkiye. These patterns have a significant impact on Türkiye, particularly in the southern regions. Over the open sea, southern flows have the tendency to accelerate. The existence of a temperature gradient at the 850 hPa level in the relevant region, along with the ascending flow zone of the 500 hPa trough, indicates the effect of the frontal system in this area. Furthermore, the presence of these patterns in the low-pressure center in the west and the high-pressure center in the east results in strong pressure gradients throughout a broad region of Türkiye. Consequently, extreme winds are experienced during the period when these patterns are active. Examples of this pattern include nodes C1, C2, C7, and C9. A secondary dominant pattern is evident over Türkiye and the southern Mediterranean region, encompassing Cyprus. This pattern is characterized by the presence of a low-pressure center originating from the Mediterranean and an upper-level 500 hPa wave trough. Additionally, a frontal system with a temperature gradient at the 850 hPa level affects the region. The flow pattern described has a stimulating impact on the western areas of Türkiye, particularly when combined with the north-eastern winds flowing over the Black Sea. Similarly, in the Mediterranean region, the movement is intensified by winds blowing from the south. When these patterns are successful, there is an increase in the number of days with extreme winds in Türkiye. This pattern was detected at node C4. In the last pattern, the occurrence of windy days in regions with a high-pressure center over central Europe and a low-pressure center over East Africa, Asia, or Russia might be attributed to the significant pressure and temperature gradients in these locations. Illustrative nodes are C3, C6, C5, and C8. Additionally, it has been observed that the initial three EOF outcomes effectively illustrate the variations in patterns that were gathered from the SOM calculations. The rise in elevation observed when moving from the coastal regions to the inland sections of Türkiye, particularly during winter, is a contributing element to the disparity in heating between the sea and the land. Examining temperature gradients across all SOM nodes reveals that synoptic and large-scale patterns also reflect disparities in land and sea warming, particularly at the 850 hPa level over the Mediterranean. These variations are attributed to synoptic and large-scale patterns. This phenomenon is also observed in SOM nodes C2, C5, C6, C7, and C8. Temperature differences between the sea and land on the Black Sea coast are found in SOM nodes C1 and C2. Pattern C4 is characterized by a low-pressure center and its related frontal system, which exerts influence throughout the entirety of Türkiye. The pattern demonstrates the impact of significant temperature variations throughout the whole region of Türkiye. Temperature gradients, particularly those running parallel to coasts, are a result of the significant temperature difference between the sea and land. As a follow-up study, the days for each SOM node are identified, and the NEWA dataset is used to calculate extreme wind event frequency for each NEWA grid point. Frequency analysis reveals the association between SOM node climate conditions and extreme wind events. It has been found that some SOM clusters effectively represent the known patterns of the research region. However, several clusters within the SOM did not accurately depict the expected trends during the DJF period. In addition, significant pressure gradients seem to produce extreme wind zones across Türkiye during the DJF period. The presence of low-pressure and high-pressure zones, namely in the Mediterranean and the Black Sea, along with the associated frontal systems, significantly influences the frequency of days characterized by extreme winds. Moreover, the Low-Pressure Center, predominantly situated above the Black Sea, exerts a substantial influence on the strengthening of winds as a result of its frontal effect in the region. Finally, the Azores and Siberia, renowned for their high-pressure systems, have a considerable impact on the formation of huge gradients over Türkiye. These gradients are crucial for the formation of extreme winds over Türkiye. Furthermore, the local geography, characterized by Türkiye’s three-sided coastal location and the existence of elevated plateaus and mountains in the interior, also contributes to the amplification of these systems. The occurrence of mountain and sea winds is a localized phenomenon. The channel and valley impact of these systems, which influence the movement in hilly regions, also contribute to local acceleration or reduction at a small or regional level. The existence of mountain ranges running from west to east in the northern and southern parts of Türkiye, along with the rugged mountainous terrain in the Aegean region, leads to an acceleration of the flow of air. This acceleration is caused by the influence of the mountains, as well as the reinforcement of the airflow towards the continental interior due to the formation of channels in the mountainous areas. In conclusion, this study also examined the correlation between the k-means and SOM clustering results. In the scope of this analysis, the relevant statistics were extracted using both clustering outcomes in terms of event-day counts at both levels and mean wind speed at lower-levels. Even SOM clustering efforts were focused on the winter period, the SOM results show similarities with the all-time k-means clustering results in terms of lower-level extreme wind events. To gain a deeper comprehension of the correlation between k-means clustering and SOM, one might employ several statistical techniques to do computations. Furthermore,there are plans to carry out a comparable investigation for the remaining three sets of seasonal data. In addition, it would be worthy to conduct future study utilizing more precise grids or alternative data reduction approaches to enhance our findings. Lastly, future study may utilize Z500, MSLP, and T850 gradient computations to do more precise quantitative calculations. This will furthermore be beneficial for substantiating study findings and validating the correctness of the model.