A Modeling Study on the Downslope Wind of “Katevatos” in Greece and Implications for the Battle of Arachova in 1826

Downslope winds and lee gravity waves are common features of mountainous environments. A similar weather type at Mt. Parnassos in Arachova, Greece is known as “Katevatos” and has devastating results for the population and visitors at the local touristic resorts. In this study, we analyze three incidents of this atmospheric pattern at local scale resolution (1 × 1 km) with WRF model. This is the first study of this local weather hazard, and the following key factors are identified. (I) The main synoptic forcing is the propagation of an upper-level trough from central Europe towards the Balkans. (II) The associated generation of a surface low-pressure system over the Aegean Sea results in a northeast flow in the lower troposphere that is perpendicular to the main topographic ridge of Mt. Parnassos. (III) Generation of gravity waves and downward reflection of wave energy at the critical level between the upper level flow and the undercutting northeast current result in the formation of “Katevatos” downslope wind at the lee side of the mountain. This hurricane-scale wind is accompanied with horizontal transport of frozen rain and snow from the mountain tops towards the village of Arachova. This wind pattern appeared also during the battle of Arachova in November 1826 between the Greek and Ottoman forces resulting in enormous casualties due to the adverse weather conditions.


Introduction
The year 2021 marks the anniversary for the 200 years since the Greek revolution of 1821 against the Ottoman Empire that gave birth to the modern-era Greek state. One of the most important conflicts during this revolution took place in November 1826 at the mountainous village of Arachova that is located at 1000 m elevation at the south slope of Mt. Parnassos. According to the historic sources, a severe weather type that was prevailing during these days included very low temperatures, gale-scale winds, and adverse snowstorm [1][2][3][4][5]. The Greek barracks during the battle were at the lower elevations, allowing the troops to take cover and protection inside the village houses, while the Ottoman forces remained at higher elevations at the mountain slopes and thus were exposed to the extreme weather. This condition resulted in an almost complete elimination of the Ottoman forces due to the cold, the wind, and the severe snowstorm. This specific pattern is well known to the locals as "Katevatos", which obviously reflects the katabatic nature of the wind and the associated phenomena. However, we need to distinguish here that the term "katabatic" is commonly reserved in meteorology for weak to moderate flows that occur from the cooling of the air masses adjacent to a mountain slope, while "katevatos" is actually a hurricane-force downslope windstorm of great local importance with wind speeds reaching up to 100 m s −1 .
In general, the complex topographic features of Greece including a 15,000 km long coastline, over 6000 islands and a mountainous inland with steep orography, canyons, valleys, etc., result in the generation of numerous local weather patterns. Nowadays, these unique sceneries are the attraction of millions of visitors from around the world. However, the formation of local scale weather types played an important role and affected the fate of several battles and wars throughout the long Greek history. An early example is already evident for the battle of Salamis in 480 B.C., when the Greek admirals took advantage of a local sea-breeze to defeat the Persian fleet in the straits of Salamis [6].
Several local katabatic winds are known worldwide (e.g., "Bora" in the Adriatic Sea [7], "Mistral" in south France (Rhone valley) [8], etc.) associated with the formation of mountain gravity waves. Foehn-type winds are also typical in mountainous regions (e.g., "Chinook" in the US Rocky Mountains [9], "Livas" in Crete [10], etc.). Due to their local nature, such flows are not commonly included in mesoscale forecasts and analysis fields [11]. For example, [12] describe how the formation of a Foehn flow over the leeward slopes of Crete can only be described at a modeling grid-space of 1 × 1 km or finer. As a result, the "Katevatos" wind that is discussed here has not yet been extensively analyzed, despite its great importance for civil protection at the very popular touristic destination of Arachova and its' historic role. The specific area is one of the most important touristic destinations in Greece, attracting thousands of visitors per year due to winter sports activities but also due to the very important archaeological site of Delfi that is located 12 km west of Arachova. For this reason, apart from the profound scientific importance, a deeper understanding of this phenomenon is also needed in terms of civil security and protection against extreme weather events. In this study, we investigate the various physical processes that take place along several atmospheric scales (from synoptic to mesoscale and local) leading to the formation of this unique weather type. A conceptual approach of katabatic winds formation as well as the methodology used in this study are described in Section 2, followed by the analysis of our results in Section 3. Section 4 includes the conclusions as well as a discussion on the historical implications of "Katevatos".

Theoretical Background
Formation of lee gravity waves is a common feature in atmospheric fluid dynamics. Lee waves form due to the vertical oscillation of the air mass, after the flow has been vertically perturbed to overcome a mountain ridge. The existence of such waves is often evident by the formation of lenticular clouds at the lee side of the mountain and by the generation of severe downslope winds ( Figure 1). Considering Newton's second law, the balance at the vertical direction takes the form of Equation (1), where z' is the mean elevation of the air mass parcel, z is the vertical elevation, g is the gravitational constant and ρ is the density of the air mass. Assuming a vertically stabilized density field with denser fluid laying below the lighter ones and a constant rate of density change in the vertical ∂ρ ∂z < 0, constant , Equation (1) obtains a solution in the wave form of Equation (2). This solution shows that the air mass parcel is oscillating in the vertical while traveling horizontally. Parameters A and B in Equation (2) are constants, and N is the classical Brunt-Vaisala frequency (Equation (3)) where θ is the potential temperature (e.g., [13]). In a stable air mass where N 2 > 0 a vertically perturbed air parcel will have the tendency to come back to its equilibrium level, and the frequency of the oscillation is given by N.
In an unstable air mass, N 2 < 0 and the vertically perturbed air parcel will move away from the equilibrium level (i.e., free convection). Favorable tuning between atmospheric and topographic properties can lead to the formation of downslope windstorms analogous to a hydraulic jump when the flow over a mountain transitions from subcritical to supercritical, or when vertically-propagating mountain waves undergo partial reflection at wave-induced critical layers (wave breaking) or atmospheric critical layers (e.g., wind shear regimes) [14][15][16][17][18].
Atmosphere 2021, 12, 993 3 of 16 back to its equilibrium level, and the frequency of the oscillation is given by N. In an un stable air mass, N 2 < 0 and the vertically perturbed air parcel will move away from the equilibrium level (i.e., free convection). Favorable tuning between atmospheric and topo graphic properties can lead to the formation of downslope windstorms analogous to a hydraulic jump when the flow over a mountain transitions from subcritical to supercriti cal, or when vertically-propagating mountain waves undergo partial reflection at wave induced critical layers (wave breaking) or atmospheric critical layers (e.g., wind shear re gimes) [14][15][16][17][18]. (3

Model Configuration
The modeling study is performed with WRF-ARW v4.1.3 model [19,20] set up in a 3 grid structure (12 × 12 km, 3 × 3 km, 1 × 1 km grid space) and 41 vertical sigma levels up to 50 mb. Initial and boundary fields are from the National Centers for Environmenta Prediction (NCEP) final analysis dataset (FNL) at 1° × 1° resolution. The daily updated Sea Surface Temperature (SST) is from the NCEP 1° × 1° analysis. A 6 h spin-up time is as sumed in all runs. The model microphysics are computed with the Morrison 2-momen scheme [21], shortwave and longwave radiation is computed with RRTMG scheme [22] the revised MM5 scheme and the NOAH Unified model are used for surface layer and soil physics, respectively [23,24], the YSU scheme is used for boundary layer physics [25 and the Grell 3D ensemble scheme [26,27] is used for cumulus parameterization.
Three separate runs are performed for the incidents of 13 February 2019, 6 January 2020 and 15 February 2021. The above days are selected due to their unusually high speeds and due to their characterization by the local community of Arachova as typical "Kateva tos" events. For all three days, the recorded wind speed at the station of the National Ob servatory of Athens in Arachova [28] is the maximum monthly record. The station is lo

Model Configuration
The modeling study is performed with WRF-ARW v4.1.3 model [19,20] set up in a 3-grid structure (12 × 12 km, 3 × 3 km, 1 × 1 km grid space) and 41 vertical sigma levels up to 50 mb. Initial and boundary fields are from the National Centers for Environmental Prediction (NCEP) final analysis dataset (FNL) at 1 • × 1 • resolution. The daily updated Sea Surface Temperature (SST) is from the NCEP 1 • × 1 • analysis. A 6 h spin-up time is assumed in all runs. The model microphysics are computed with the Morrison 2-moment scheme [21], shortwave and longwave radiation is computed with RRTMG scheme [22], the revised MM5 scheme and the NOAH Unified model are used for surface layer and soil physics, respectively [23,24], the YSU scheme is used for boundary layer physics [25] and the Grell 3D ensemble scheme [26,27] is used for cumulus parameterization.
Three separate runs are performed for the incidents of 13 February 2019, 6 January 2020 and 15 February 2021. The above days are selected due to their unusually high speeds and due to their characterization by the local community of Arachova as typical "Katevatos" events. For all three days, the recorded wind speed at the station of the National Observatory of Athens in Arachova [28] is the maximum monthly record. The station is located at elevation of 1069 m at 38 •

Modeling Results and Analysis of the Atmospheric Processes Forming "Katevatos"
The analysis is performed separately for the three cases of "Katevatos" in Arachova to identify the main atmospheric features that drive the phenomenon. All three events present similar findings in terms of the general atmospheric structure. For this reason, in the following paragraphs we present plots only from the February 2019 case, except when certain features are found to vary significantly between the runs. In such a case, separate plots are shown for each run to discuss similarities and differences.
At synoptic scale, the main atmospheric disturbance is found to be an eastward propagating trough that is evident by the geopotential height and temperature field at 500 mb in Figure 2a, on 13 February 2019, 18:00. The trough is deepening over Greece with geopotential heights reaching 5400 m and temperatures falling below −36 • C. The wind flow at 500 mb is dominated by the trough circulation, and as seen by the wind field on 14 February 2019, 00:00 UTC in Figure 2b, the wind direction over central Greece at this tropospheric level is SSE. Cyclogenesis is evident at the lower tropospheric levels on 14 February 2019, 00:00 UTC, as seen by the sea level pressure contours in Figure 2c. The center of the surface depression is found over the southern parts of Aegean Sea, accompanied by the "packing" of isobars along the SW-NE axis over the Aegean. The combination of this low-pressure system with the higher sea-level pressure over central Europe results in the generation of a strong northeast wind flow channeling through the Aegean with maximum speeds on 14 February 2019, 00:00 UTC exceeding 30 m s −1 at 850 mb ( Figure 2d). The acceleration of wind speed that is evident in this plot over mainland Greece, the Aegean islands and the Dalmatian Alps, indicates the interaction between the wind flow and the orography. As will be explained in the following sections, backing of the wind from SSE at 500 mb to NNE at 850 mb over the area of interest plays an important role for the generation of the "Katevatos" downslope flow. Furthermore, the advected NE air masses are enriched with vapor from the Aegean Sea, and thus, they become lighter and more unstable. As seen by the total modeled accumulated snowfall from all runs at the intermediate model domain (Figure 3), the northeast current is associated with heavy snowfall over central Greece. This is mostly evident for the upwind slopes that are perpendicular to the incoming flow, due to mechanical convection when the advected air masses reach the inland mountain ridges. As seen also in Figure 3, the greater area of Mt. Parnassos is one of the most affected by snowfall areas in mainland Greece in our simulations.
A more detailed representation of the local atmospheric conditions for all three cases is given in Figure 4 from the inner modeling domain at 1 × 1 km grid space. In all cases, the highest wind speeds (reaching 33 m s −1 in Figure 4a,b) are evident at the lee (south) slope of Mt. Parnassos and particularly at the greater area of Arachova. It is worth mentioning that wind speeds above 33 m s −1 are classified as hurricane in the Saffir-Simpson scale. The 2 m-temperature is also very low in the greater area of interest, reaching down to -8 • C in Figure 4c. Although freezing conditions are evident overall in the domain, the 2 m temperatures at the lee side of the mountain are still much higher than the windward ones, indicating a Foehn type adiabatic heating of the air masses as they descend the lee slopes ( Figure 5). This combination of extreme wind speed and very cold temperatures marks the appearance of "Katevatos" as the most extreme downslope flow over the entire modeling domain. Downslope windstorms are determined by the vertical structure of the atmosphere above the mountain top and by the destabilization of the relatively laminar windward flow as the air masses climb over the mountain top. The resulting atmospheric structure that is shown by the isentropic lines in Figure 6a is typical of vertically propagating gravity waves and downslope windstorm formation as described already in early theoretical studies, e.g., [15,16]. The modeled flow in Figure 6a corresponds to the high drag regime cases that are shown in Figure 1b,c in [15] and (Figures 4b and 6b) in [16]. The cold air of the northeast current is forced to climb the windward slope in order to overpass the ridge of Mt. Parnassos. When the flow crosses the top of the mountain, the denser air parcels fall towards the ground due to buoyancy and the oscillation is advected several kilometers downwind. The generation of vertically propagating gravity waves is evident by the upwind tilting isentropes above the mountain top. Increased stability is also evident near the mountain top in Figure 6a by the "packing" of isentropic lines between 2 and 4 km height above the upwind slope. The tropopause is also evident by the very stable layer above 9 km. A critical level is defined at about 6 km height in the upwind flow where the cross-barrier flow (i.e., the meridional wind component) is zero (Figure 6a). This critical level is produced at the interface between the SE winds associated with the trough circulation in the middle troposphere and the undercutting NE  (Figure 6b). According to [15], downward reflection of the gravity waves occurs at critical layers with Richardson number less than 1 /4 . Such areas are evident in Figure 6c at about 6-7 km height above the windward slope and also at the tropopause (9-10 km). The downward reflection of gravity wave energy that occurs at this critical level, results in the generation and intensification of the downslope wind. The constant potential temperature surfaces (isentropic lines) in Figure 6a are actually the traveling paths of the air parcels.
Atmosphere 2021, 12, 993 A more detailed representation of the local atmospheric conditions for is given in Figure 4 from the inner modeling domain at 1 × 1 km grid spac the highest wind speeds (reaching 33 m s −1 in Figure 4a,b) are evident at slope of Mt. Parnassos and particularly at the greater area of Arachova. It tioning that wind speeds above 33 m s −1 are classified as hurricane in the S scale. The 2 m-temperature is also very low in the greater area of interest, r to -8 °C in Figure 4c. Although freezing conditions are evident overall in th 2 m temperatures at the lee side of the mountain are still much higher than ones, indicating a Foehn type adiabatic heating of the air masses as they d slopes ( Figure 5). This combination of extreme wind speed and very cold marks the appearance of "Katevatos" as the most extreme downslope flow modeling domain. Earlier studies on the effects of asymmetric mountain slopes to the generation of downslope windstorms have shown that the mountain shape has a peculiar influence on gravity wave structure and that the wave drags are significantly enhanced for mountain ridges with gentle windward slopes and steeper lee slopes [29][30][31][32]. A steeper lee orography is also the case for Mt. Parnassos as can be seen by the topographic contour lines in Figure 4 and by the cross-section plots (e.g., Figure 6). The steeper lee slope acts as an additional enhancing factor to the windstorm force. Obviously, the flow domain is in fact more complex, especially at areas of cloud formation; however, our results are in line with the idealized theoretical concepts for the formation of lee gravity waves and downslopes. Favorable tuning between the critical layer height, mountain top heights and flow properties (wind speed and shear, stability) can result in partial reflection of waves towards the ground and generation of downslope winds [14][15][16]. Such conditions are evident in Figures 6a and 7a where the abrupt wind shift at about 6-7 km acts as critical surface for wave reflection.   Downslope windstorms are determined by the vertical structure of the atmosphere above the mountain top and by the destabilization of the relatively laminar windward flow as the air masses climb over the mountain top. The resulting atmospheric structure that is shown by the isentropic lines in Figure 6a is typical of vertically propagating gravity waves and downslope windstorm formation as described already in early theoretical studies, e.g., [15,16]. The modeled flow in Figure 6a corresponds to the high drag regime cases that are shown in Figure 1b,c in [15] and Figure 4.6b in [16]. The cold air of the northeast current is forced to climb the windward slope in order to overpass the ridge of Mt. Parnassos. When the flow crosses the top of the mountain, the denser air parcels fall towards the ground due to buoyancy and the oscillation is advected several kilometers downwind. The generation of vertically propagating gravity waves is evident by the upwind tilting isentropes above the mountain top. Increased stability is also evident near the mountain top in Figure 6a by the "packing" of isentropic lines between 2 and 4 km height above the upwind slope. The tropopause is also evident by the very stable layer above 9 km. A critical level is defined at about 6 km height in the upwind flow where the crossbarrier flow (i.e., the meridional wind component) is zero (Figure 6a). This critical level is produced at the interface between the SE winds associated with the trough circulation in the middle troposphere and the undercutting NE current arriving below (Figure 6b). According to [15], downward reflection of the gravity waves occurs at critical layers with Richardson number less than ¼. Such areas are evident in Figure 6c at about 6-7 km height above the windward slope and also at the tropopause (9-10 km). The downward reflection of gravity wave energy that occurs at this critical level, results in the generation and intensification of the downslope wind. The constant potential temperature surfaces (isentropic lines) in Figure 6a are actually the traveling paths of the air parcels. Earlier studies on the effects of asymmetric mountain slopes to the generation of downslope windstorms have shown that the mountain shape has a peculiar influence on gravity wave structure and that the wave drags are significantly enhanced for mountain ridges with gentle windward slopes and steeper lee slopes [29][30][31][32]. A steeper lee orography is also the case for Mt. Parnassos as can be seen by the topographic contour lines in Figure 4 and by the cross-section plots (e.g., Figure 6). The steeper lee slope acts additional enhancing factor to the windstorm force. Obviously, the flow domain is in fact more complex, especially at areas of cloud formation; however, our results are in line with the idealized theoretical concepts for the formation of lee gravity waves and downslopes. Favorable tuning between the critical layer height, mountain top heights and flow properties (wind speed and shear, stability) can result in partial reflection of waves towards the ground and generation of downslope winds [14][15][16]. Such conditions are evident in Figures 6a and 7a where the abrupt wind shift at about 6-7 km acts as critical surface for wave reflection. The downslope flow is shown in Figure 8 at the time of maximum horizontal speed occurrence for each simulation. It is defined by a shallow layer of hurricane-scale wind speeds (>40 m s −1 ) rolling down from the mountain top towards the village of Arachova. As seen previously in Figure 6b, the horizontal wind during the onset stages of the phenomenon presents a significant vertical shear and is actually rotating in the vertical between 0 • and 360 • downwind of the lee slope. After the windstorm is fully established, northeast directions (0-90 • ) prevail overall the domain from the surface up to 6 km height (not shown).
Cloud development and generation of precipitation during "Katevatos" events are mostly related to the windward slope of Mt. Parnassos. Mechanical uplifting of the moist air masses coming from the Aegean Sea result in vertical cloud development and increased mixing ratios of liquid (cloud and rain droplets) and frozen (ice, snow, hail) condensates at the windward slope as seen in Figure 9a. Most of the precipitation falls in the form of snow and supercooled rain droplets before the air parcels overcome the mountain top. However, during the most intense phases of the downslope windstorm, large amounts of rain and snow are transported horizontally (parallel to the slope) by the wind, from the mountain top towards lower elevations and Arachova. This process can be seen in  form of snow and supercooled rain droplets before the air parcels overcome the mountain top. However, during the most intense phases of the downslope windstorm, large amounts of rain and snow are transported horizontally (parallel to the slope) by the wind, from the mountain top towards lower elevations and Arachova. This process can be seen in Figure 9b by comparing the vertical structure of the precipitating clouds over the windward slope to the very shallow layer of condensates rolling down the lee slope of the mountain.

Conclusions and Discussion
Our study presents an analysis of the main atmospheric processes that interplay for the generation of a unique and understudied winter weather type that appears in Mt. Parnassos. The specific "Katevatos" pattern results from a combination between synoptic scale forcing, mesoscale wind flow and local-scale perturbations caused by the steep mountain topography and the specific landscape of Arachova (mountain passages, canyons, slopes etc.). The main features that are identified in our analysis are as follows: (I)

Conclusions and Discussion
Our study presents an analysis of the main atmospheric processes that interplay for the generation of a unique and understudied winter weather type that appears in Mt. Parnassos. The specific "Katevatos" pattern results from a combination between synoptic scale forcing, mesoscale wind flow and local-scale perturbations caused by the steep mountain topography and the specific landscape of Arachova (mountain passages, canyons, slopes etc.). The main features that are identified in our analysis are as follows: (I) The vertical displacement of the NE flow to overcome Mt. Parnassos results in vertical oscillation of the air mass at the lee side of the mountain due to buoyancy; (II) Generation of vertically propagating mountain waves is evident over the mountain top and up to the tropopause; (III) Reflection of wave energy towards the surface occurs at the critical level that appears along the interface between the upper level SE flow and the undercutting NE wind resulting in the generation of the downslope wind; (IV) Although precipitation occurs mostly over the windward slope, an avalanche of ice condensates and supercooled rain droplets are also carried downstream by the "Katevatos" current, finally resulting in a hazardous event that combines extreme wind, very low temperatures and snowfall. In fact, "Katevatos" is even more hazardous than described in our modeling analysis because even at this dense grid space of 1 × 1 km, the simulations are not capable of reproducing the highest recorded wind gusts. Very local topographic properties such as mountain openings and canyons increase the severity of the phenomenon due to channeling effects and local accelerations of the downslope flow. However, the similar mesoscale fields that are found for all three analyzed cases suggest two forecasting guidelines for the appearance of Katevatos: (1) Formation of a low-pressure system centered east of the Parnassos mountain (i.e., in the Aegean Sea) and (2) Formation of a NE wind current arriving at the windward slope from the surface up to about 6 km accompanied by abrupt wind shear above this level.
After analyzing the current events, we can postulate that the weather patterns shown in this study are responsible for the "Katevatos" event that was recorded during the battle of Arachova in November of 1826 between Greek and Ottoman forces. As described in the historical sources [1][2][3][4][5], this critical battle took place under weather conditions similar to the ones found in this study. The prevalence of the Greek forces in the battle of Arachova was crucial for the fate of the Greek revolution, and this conflict had a major contribution to the final stages of the Greek-Ottoman war. The battle lasted for 5-6 days and ended on 24 November 1826. The Ottoman forces, leaded by general Moustambeis, were trapped at the elevated slopes of Mt. Parnassos, while the Greek forces of general Karaiskakis remained protected from the extreme weather at lower elevations inside Arachova (see Figure 10). When the Ottomans tried to escape through the mountain passages, the appearance of "Katevatos" made their attempt impossible and resulted in heavy casualties, practically eliminating their forces [1][2][3][4]. The impact of "Katevatos" was so important that it was attributed to divine intervention by the Greeks [5].

Arachova
Greek troops Ottoman troops