The Road Map to Classify the Potential Risk of Wind Erosion

: Environmental degradation, for example, by wind erosion, is a serious global problem. Despite the enormous research on this topic, complex methods considering all relevant factors remain unpublished. The main intent of our paper is to develop a methodological road map to identify key soil–climatic conditions that make soil vulnerable to wind and demonstrate the road map in a case study using a relevant data source. Potential wind erosion (PWE) results from soil erosivity and climate erosivity. Soil erosivity directly reﬂects the wind-erodible fraction and indirectly reﬂects the soil-crust factor, vegetation-cover factor and surface-roughness factor. The climatic erosivity directly reﬂects the drought in the surface layer, erosive wind occurrence and clay soil-speciﬁc winter regime, making these soils vulnerable to wind erosion. The novelty of our method lies in the following: (1) all relevant soil–climatic data of wind erosion are combined; (2) different soil types “sand” and “clay” are evaluated simultaneously with respect to the different mechanisms of wind erosion; and (3) a methodological road map enables its application for various conditions. Based on our method, it is possible to set threshold values that, when exceeded, trigger landscape adjustments, more detailed in situ measurements or indicate the need for speciﬁc management.


Introduction
Soil erosion, along with soil compaction, acidification, organic matter decline, soil biodiversity decline, soil pollution and others, is one of the crucial negative externalities of intense agricultural management that reduce both production and ecological soil functions. Mitigation measures to these degradation processes are frequently sought immediately after significant public interest and social demand are raised. Aroused legislative pressure in terms of soil degradation was, for instance, put on the filtering function of soil for groundwater and surface water protection [1]; nevertheless, the problems seemingly relevant only to agriculture have often been omitted so far. Soil degradation processes involving erosion globally threaten food security, human health, and ecosystem services and stability [2,3]. Various anthropogenic and natural disturbances of the landscape can increase the risks of wind erosion, which can be intensified with drought. It can be assumed that wind erosion will be even more serious in the future due to lower soil moisture content [4][5][6][7][8]. Dry and bare soil surfaces are highly susceptible to wind erosion, while soil particles driven by wind endanger emerging plants of agricultural crops. Longterm phenological data as a response to climate change and detected changes in their onset with a link to wind erosion in temperate climate were investigated by [9][10][11][12][13]. In terms of dry and semidry areas, [14] pointed to the important role of vegetation cover in wind erosion reduction.
However, according to [15], the majority of northern European countries claim antidegradation of soil (or at least point pollution and water erosion and landslides) in their policy documents. Today a paradigm shift is needed in this matter: to move from environmental protection to sustainable use and management. Though it cannot be agreed upon more, specific tools for wind-erosion control have been neglected, partly due to a lack of appropriate methods of its evaluation. As the physics of wind erosion is complex, soil and atmospheric and land-surface processes must be taken into account to assess the wind-erosion susceptibility of soils. Webb and Pierre [16] critically reviewed the drivers of anthropogenic dust emissions and current evaluation approaches. The existence of gaps in the research into the soil-erosion processes in Europe, especially the wind-erosion processes was emphasized by [17]. The common indicators and measurements of wind erosion and air quality used for monitoring various ecosystems by scientists and land managers were summarized by [18] while also taking into account data availability and the technical difficulty of understanding how the measurements indicate the erosion risk or outcomes. Most of these models require detailed input data, which are not always easily available and cannot be simply adapted to conditions or climates different from those in which they were developed [19]. Studies on a regional scale produce varied percentages of the affected area, but the comparability of these data is limited by the differences in methodology and definitions.
Wind-erosion modelling started in the early 1960s for the semiquantitative estimation of soil losses, with the wind erosion equation (WEQ) [20] being considered a standard. This method was globally used to calculate the potential risk of wind erosion [21][22][23]. In contrast, some precise geographic information system (GIS)-based wind-erosion models, such as the Wind Erosion on European Light Soils (WEELS) assessment, are highly data demanding, and thus their practical application is limited.
Wind erosion is not only a concern for agricultural land degeneration but also air quality due to the emission of fine particulate matter into the atmosphere. Aeolian sediment transport and dust emission pose serious hazards to human health. Increasing pressures on land use leads to dust emission and atmospheric pollution in terms of particulate matter (PM) dispersion. The main source of PM in the countryside is traffic on unpaved roads, and even though wind erosion does not work by the same mechanism it can still contribute, especially when the soil surface is bare. Ref. [24] quantified anthropogenic dust emissions from agricultural land under arid conditions when all the dust schemes showed a dependence of anthropogenic dust on land cover over agricultural lands. Yulevitch et al. (2020) [25] identified wind-velocity-triggering dust emissions from loess soils as 4 m s −1 . Dust emission processes have major implications for soil loss and human exposure to air pollution. [26] studied the impacts of soil disturbances by human activities on soil aggregation and dust fluxes, and recorded the substantial loss of PM10 under most experimental conditions. These researchers highlight the significant implications for soil nutrient resources in annual balance and management strategies, as well as for PM loading in the atmosphere and the risk of air pollution. Moreover, the atmosphere presents a new vehicle for microplastics, which can be transported by wind together with soil particles and then spread to a wider environment. The dependency of fibrous and nonfibrous airborne microplastics on wind while emphasizing the need to include airborne pathways when consolidating the impact of microplastics on the broader environment and human health indicated [27]. Fine suspended particulates in the atmosphere have also caused disruption to the operation of aircraft and vehicles, and present a health hazard to humans and animals.
The current state of erosion research lacks the knowledge of where and when wind erosion occurs in Europe, as well as the erosion intensity, which poses a threat to agricultural productivity [17].
For this reason, we present a general methodology to identify and spatially express the risk of wind erosion (as risk categories from 1 to 6) while emphasizing all phenomena that need to be taken into account. We do not intend to provide the results in the absolute values of soil loss but rather a relative comparison of areas with different exposure risks to wind erosion. Our study thus focuses on the potential wind erosion (PWE) that represents the soil-climatic-driven susceptibility to erosive wind. Alternatively, by means of the actual wind erosion (AWE), we understand real soil loss under the field conditions caused by a high wind velocity with an erosive effect. Additionally, we employ the case study for central European conditions to practically demonstrate the methodological protocol. The method might help to identify and prioritize areas with the highest risk of wind erosion and thus the highest need for anti-degradation measures, and management strategies and action. This approach also works as a supportive tool for financial and management decision making.

Materials and Methods
The aim of this study is to generate a methodology to identify and spatially express the risk of wind erosion by taking into account all necessary phenomena/factors (not only soil). The study is focused on PWE, and the methodology is applied and explained in a case study of the Czech Republic located in central Europe.

Description of the Case Study Area
The area of interest is located in central Europe (12 • -19 • E, 48 • -51 • N) and represents part of the Czech Republic, namely, the area where the risk of wind erosion is expected due to the relatively dry and windy climate (i.e., agricultural land up to an altitude of 500 masl)-see Figure 1. Ref. [28] reports a long-term (period 1990-2019) average annual precipitation of 614.4 mm, an average annual temperature of 9 • C, an average wind velocity of 2.5 m s −1 in the spring season (March-May) and an average wind velocity of 2.2 m s −1 in the autumn season (September-November). The lowest wind velocities occur during summer, while spring and autumn are windier [29]. For this reason, we present a general methodology to identify and spatially express the risk of wind erosion (as risk categories from 1 to 6) while emphasizing all phenomena that need to be taken into account. We do not intend to provide the results in the absolute values of soil loss but rather a relative comparison of areas with different exposure risks to wind erosion. Our study thus focuses on the potential wind erosion (PWE) that represents the soil-climatic-driven susceptibility to erosive wind. Alternatively, by means of the actual wind erosion (AWE), we understand real soil loss under the field conditions caused by a high wind velocity with an erosive effect. Additionally, we employ the case study for central European conditions to practically demonstrate the methodological protocol. The method might help to identify and prioritize areas with the highest risk of wind erosion and thus the highest need for anti-degradation measures, and management strategies and action. This approach also works as a supportive tool for financial and management decision making.

Materials and Methods
The aim of this study is to generate a methodology to identify and spatially express the risk of wind erosion by taking into account all necessary phenomena/factors (not only soil). The study is focused on PWE, and the methodology is applied and explained in a case study of the Czech Republic located in central Europe.

Description of the Case Study Area
The area of interest is located in central Europe (12°-19° E, 48°-51° N) and represents part of the Czech Republic, namely, the area where the risk of wind erosion is expected due to the relatively dry and windy climate (i.e., agricultural land up to an altitude of 500 masl)-see Figure 1. Ref. [28] reports a long-term (period 1990-2019) average annual precipitation of 614.4 mm, an average annual temperature of 9 °C, an average wind velocity of 2.5 m s −1 in the spring season (March-May) and an average wind velocity of 2.2 m s −1 in the autumn season (September-November). The lowest wind velocities occur during summer, while spring and autumn are windier [29].

Soil and Climate Factors of PWE
Wind erosion and dust emissions are controlled by climate, soil and vegetation properties [16]. The simplified PWE is given by the soil erosivity (SE) and climatic erosivity (CE).

Soil Erosivity (SE)
SE generally means the ability of a soil to resist wind forces. With a certain amount of generalization, we can follow the pattern of [30,31] and define SE as a combination of the wind-erodible fraction (EF) factor, soil-crust (SC) factor, vegetation-cover (VC) factor and surface-roughness (SR) factor.
The EF factor defined by [32] is based on the soil's texture and chemical properties, such as the content of sand, silt and clay and the content of organic-matter calcium carbonate, while [33] simultaneously employed the clay and organic-matter content to obtain the SC factor.
We are fully aware that the mechanism of wind erosion differs with soil type. Across many studies (for instance [34]), clay soils are not considered to be vulnerable to wind erosion. That is why they are not generally linked to any risk category in terms of the PWE. Some authors [35][36][37][38][39][40] claim that specific clay soils are highly vulnerable to the erosive effects of wind after winter under specific weather conditions. Notably, mainly clay soils show substantial annual changes in soil cohesion. Periodic freezing and thawing of the soil profile together with an enormous water content cause aggregate disintegration. The combination of particular soil properties with a specific winter regime makes clay soil vulnerable to wind erosion [41].
In terms of the PWE estimation, we approached soil grouping with a certain measure of generalization. Soils with clay contents up to 45% were assigned as sandy-to-loam soils (throughout the whole paper, they were indexed as "SAND"). This category thus includes soils that have light texture, lighter to medium-heavy texture, and medium-heavy texture. That is, we include soils of sandy, loamy sand, sandy loam and loamy textures. The second group includes soils with a clay content above 45%, i.e., clay loam soils, clay soils and clay (hereafter indexed as "CLAY"). While speaking of sandy or clay soils, we indicate only that part of the soil profile that is directly exposed to climatic influence.
Nonetheless, the sandy soils generally evince a higher vulnerability to wind erosion than clay soil, so in terms of the SE degree, we assign them degrees of 2, 4 or 6, whereas clay soil is assigned degrees of 1, 3 or 5.
Soil erosivity of sandy soil (SE SAND ) Lighter soils can generally be seen as having a higher proportion of macropores (lower overall porosity, larger mineral grains). These soils are typically characterized by higher aeration (hence, they are often more overheated than heavier soils) and better rainwater permeability. As they are not classified as medium-heavy soils, due to their low sorption capacity, they are often quite dry, showing a higher acidity and typically a lower content of mineral nutrients (and organic substances in general). Due to their predisposition, the lightest soils are very vulnerable to wind and water erosion. This susceptibility is because particles with a diameter of 0.05-0.08 mm (the boundary between dust and fine sand) are most easily transported by wind. According to the soil vulnerability factor, we differentiate the whole group according to Table 1 in order to obtain Degrees of SE SAND (Table 2).
Case study on SE SAND The data source for soil properties was the national soil database "Estimated Pedological Ecological Unit" (EPEU) in the form of a five-digit code developed in the Czech Republic in the 1970s to evaluate the main pedological and climatological conditions of agricultural land from its production ability viewpoint. The EPEU could be successfully used for the estimation of wind-erosion risk in terms of the SE SAND degree ( Figure 2). Table 1. Description of sandy soil potentially vulnerable to wind erosion.

Soil Group Detail Description
A The most vulnerable soils, i.e., those whose upper part of the soil profile shows a strongly arenic character. These sandy soils (or partly sandy loam soils) of different genesis are developed on light soil-forming substrates.
B Endangered soils are formed by the sand-loam granularity and can be described as lighter to medium-heavy soils.
They are developed on similar parent substrates as the most erosion-endangered soils, and even on the highest levels of the soil profile we can identify finer material, most often of eolic origin (i.e., loess) or solifluction sediments with a high proportion of marl detritus. In addition to the aforementioned soil types, which often acquire modal, stagnic or gleyic subtype, we can also include lighter Haplic Luvisols.
C Soils susceptible to wind erosion represent an aggregate group, which includes a varied mixture of soil types with the variable genesis of formation and soil-forming substrate and is mainly associated with the medium-heavy (loam) grain size of the upper part of the soil profile.
They are the most fertile Chernozems and grey Phaeozems of the modal or luvic subtype. Furthermore, quality modal Rendzic Leptosols, Cambisols and lastly modal, gleyic and stagnic Fluvisols.  Soil erosivity of clay soil (SECLAY) Heavier soils are characterized by a lower proportion of macropores, thus showing lower water permeability for rainwater and overall aeration. If these soils accumulate rainwater at all, then due to high sorption capacity they lead poorly across the profile. As a result of all these phenomena, in heavier soils we encounter a double extreme; that is, in excess, the water concentrates on the soil surface and causes the site to become wet, and in the absence of water, hardening and cracking occur.
Some clay soils became vulnerable to wind erosion after they underwent the specific winter regime, which led to the creation of specific soil aggregates. The effects of moisture and temperature can cause significant changes in the aggregates and thus the erodibility of soil. Standardized methods allow aggregate-size distribution and aggregate stability, which affect wind erosion [42]. Bullock (1988) [43] measured the stability of soil aggregates Soil erosivity of clay soil (SE CLAY ) Heavier soils are characterized by a lower proportion of macropores, thus showing lower water permeability for rainwater and overall aeration. If these soils accumulate rainwater at all, then due to high sorption capacity they lead poorly across the profile. As a result of all these phenomena, in heavier soils we encounter a double extreme; that is, in excess, the water concentrates on the soil surface and causes the site to become wet, and in the absence of water, hardening and cracking occur. Some clay soils became vulnerable to wind erosion after they underwent the specific winter regime, which led to the creation of specific soil aggregates. The effects of moisture and temperature can cause significant changes in the aggregates and thus the erodibility of soil. Standardized methods allow aggregate-size distribution and aggregate stability, which affect wind erosion [42]. Bullock (1988) [43] measured the stability of soil aggregates throughout the year using wet sieving techniques. Major decreases in cohesion were found when the minimum daily air temperature fell to or below 0 • C during the winter and early spring months (caused by soil freezing at high water contents). Refs. [41,44] tried to identify the effects of freeze/thaw cycles and soil moisture on wind erosion with the help of wind-tunnel experiments. The surface-soil moisture content decreased in thawed soil and increased in frozen soil after wind erosion. The wind-driven sediment flow in the thawed soil was always greater than the flow of frozen soil with the same moisture content, but this difference became negligible at a moisture content above 3.38% (gravimetric moisture content) [44].
According to [41], the porosity of black soil with initial moisture contents on the surface of soil samples of 5-8% can increase by 1.4-1.5 times after the freeze/thaw process. The soil structure reaches a new stable state after 6 freeze/thaw cycles. The wind erosion increases with the number of freeze/thaw cycles. After 6-9 freeze/thaw cycles, the wind erosion of soil with moisture content of 5% and 7% increases by 1.2 and 2.0 times, respectively. The influence of freeze/thaw cycles on wind erosion is not significant when moisture exceeds 8%.
Therefore, in order to get SE CLAY we defined soil potentially vulnerable to wind erosion (Table 3) and following criteria of the specific winter regime, which are combined according to Matrix 1 in Table 4): (a) soil properties SE CLAY -SP (Table 1) number of freeze-thaw and thaw-freeze episodes during the days without snow cover (ii) number of days with wet soil surface, allowing soil aggregate dispersion Table 3. Description of clay soil potentially vulnerable to wind erosion.

Soil Group Detail Description
A Clay loam to clay Chernozems and Phaeozems, as deep-humid soils with a diagnostic phaeozem horizon, occurring in a relatively humid climate.
Both soil types are developed on unconsolidated heavy substrates containing carbonates in the soil profile (in the Czech Republic especially on flysch rocks, limnic Tertiary, tuffs and tuffites or chalk marlstone) and usually without a significant skeleton. Phaeozems are bound to the depressed position of Chernozems, having a pronounced hydric regime, often prone to waterlogging and are characterized by an even higher carbon content than Chernozems. Typical representatives in these groups are Chernozems pelic, Phaeozems pelic, and Phaeozems pelic-gleyic. Vertisols exhibiting similar properties were not considered for their limited area or margin of interest in this study.

B
Heavy Fluvisols, i.e., soils of flat areas along watercourses that are geologically built by fluvial or deluviofluvial sediments, respectively.
They are mostly skeleton-free, rarely slightly skeleton-like and tend to temporarily get wet. The diagnostic feature is the stratification of the soil profile with an irregular distribution of organic substances. Typical representatives are Fluvisols pelic and Fluvisols pelic-gleyic, which then exhibit marked reductomorphic features less than 60 cm from the soil mineral surface. Heavy colluvisols (not defined in the World Reference Base for Soil Resources) were not considered for their non-significant incidence.
C Soils developed on solifluction deposits and regoliths of diagenetic carbonate-silicate rocks occurring mainly in flysch areas (typically marine lower Cenozoic claystones and marlstones). Less often, they are tied to the greywacke and slate of the culm, basic effusives, chalk marlstones or rocks of the limnic tertiary.
They are modal pelozems (not defined in the World Reference Base for Soil Resources, Cambisols modal), Cambisols pelic, Regosols pelic and calcic Leptosols pelic. They are skeletonless to medium skeletal, and often take the form of a slightly stagnic character. If they are heavily stagnic with a tendency to temporarily get wet, then they are marked as stagnic or as Stagnosols pelic. Individual degrees of SE CLAY brings Table 5. Case study on SE CLAY The data source for soil properties was the national soil database EPEU. The data source for specific winter regime determination was the Czech Hydrometeorological Institute (CHMI) database, particularly the data on the "condition of bare soil surface" (CBSS) in combination with temperature data. The methodological approach to SE CLAY is described in detail in [45]. The SE CLAY degree for the area of the case study is presented in Figure 3. tendency to temporarily get wet, then they are marked as stagnic or as Stagnosols pelic.  Case study on SECLAY The data source for soil properties was the national soil database EPEU. The data source for specific winter regime determination was the Czech Hydrometeorological Institute (CHMI) database, particularly the data on the "condition of bare soil surface" (CBSS) in combination with temperature data. The methodological approach to SECLAY is described in detail in [45]. The SECLAY degree for the area of the case study is presented in Figure 3.  The SE VC factor comprises two different dimensions. First, the protective effect of soil coverage (SE VC-C )-mainly vegetation-and second, the protective effect of windbreaks (SE VC-W )-mainly tree shelters. Borrelli (2016) [46] and others expressed SE VC-C by using the percentage of the surface covered with nonerodible plant material. SE VC-W has often been based on optical porosity of windbreaks. The SE SC factor affects wind erosion by dissipating the wind erosivity, while [46] used a land cover map to derive the roughness length.
We integrated both SE VC (i.e., SE VC-C and SE VC-W ), and SE SR factors indirectly by choosing the relevant periods for CE, i.e., spring and autumn, which are characterized by: (i) Little to no crop-vegetation cover and a flat and smooth soil surface with limited roughness (short before or after sowing). (ii) Low protective effect of tree windbreaks. The optimal efficiency of the windbreak is mainly given by the optimal porosity (39%) and structure of the windbreak-multiplerow design [47]. Ref. [48] suggested that porosity between 40% and 50% is the most efficient. The highest values of optical porosity (up to 50%) are achieved in the nonvegetation period. Almost full foliage of windbreaks occurs in approximately 10 days in the spring. Autumn leaf fall lasts longer, and a gradual change in the optical porosity is longer. (iii) With a tendency to be windier, according to [29], the average wind velocity is increasing by 1-2 m s −1 in the Czech Republic during these periods. It can be concluded that wind erosion is very intense at the beginning of spring, especially after the dry winter, and has a major impact on the regional agriculture.
SE SC is relevant mainly to heavy textured soil "CLAY" and is reflected in their classification into the different risk categories in Tables 3 and 5.

Climatic Erosivity (CE)
Kouchami-Sardoo (2019) [49] developed a quantification of wind-erosion risk by incorporating the assessments of consequence and vulnerability. The model used was evaluated by scenario testing, sensitivity analysis, and wind tunnel measurements, which provided reasonable estimates of the soil vulnerability, consequence and risk to/of wind erosion. The results showed that weather factors were among the most important parameters affecting wind erosion risk.
CE is basically defined as a measure of the climatic tendency to produce conditions conducive to wind erosion. Ref. [30] expressed the CE using the equation developed by [50] employing the mean monthly wind speed (m s −1 ) at a 2 m height and the potential evapotranspiration (mm) and the precipitation (mm). However, these are macroclimatic parameters that do not take into account the threshold wind speed and soil surface.
To obtain the CE degrees, we employed a complex approach that takes into account relevant climate variables, i.e., the drought in the surface layer (CE-DSL) and risk of erosive wind occurrence (CE-REW).
CE-DSL Soils identified as vulnerable (see SE chapter) increase their susceptibility to wind erosion with the drought occurrence in their surface layer. In terms of wind erosion, the moisture of the soil surface is crucial mainly in the periods when the soil is uncovered by vegetation. Since the 1970s, remote sensing data have been of limited use for the purposes of assessing soil surface moisture [51,52]. For example, the successful use of the spatial resolution and frequency of imaging was not possible for the evaluation of soil-erosion hazards at the regional scale, especially when the land structure significantly differs. In those areas, it is worth applying data from sufficiently dense ground measurement networks of soil moisture if available.
CE-REW Wind velocity is important for assessing wind erosion mainly in terms of the AWE. In this case, it is necessary to evaluate the actual wind velocity during a specific wind erosion event and relate it to real soil loss. In the case of the PWE, actual wind velocity is not taken into account, but average long-term wind conditions along with other parameters of SE and CE define the complex vulnerability of soil to wind erosion.
Wind direction and velocity measurements are normally performed 10 m above the surface. The wind speed required to initiate soil transport is referred to as the threshold wind speed (TWS). The TWS is a particularly important parameter but is difficult to determine it in wide-scale field studies. Most studies on the relationship between the particle size distribution of sediment and wind speed have focused on wind-tunnel experiments [53][54][55][56].
However, many factors that affect wind erosion and dust emission, e.g., clod structure, crusting, and soil shear strength, are difficult to simulate concurrently in a wind tunnel; thus, the results of wind-tunnel experiments cannot completely reflect the change in the particle-size distribution of sediment with wind speed [57]. It is preferable to study the effect of wind speed on the particle-size distribution of sediment through field observations. However, it is difficult to completely collect and accurately determine the particle-size distribution of sediments at different wind speeds in the field [58,59]. The divergence in the erodibility for different sizes of particles at different wind speeds plays a critical role in changing the regularity of the aggregate-size distribution of windblown sediment with increasing wind speed.
The TWS varies with soil type. Based on wind-tunnel experiments, [60] determined the TWC for various soil types (sandy, loamy-sandy, sandy-loamy and loamy) in dry and wet conditions (3.3 m s −1 for sandy and 22 m s −1 for loamy-sandy), with the lowest/critical value of 3.3 m s −1 . According to [61], the TWS ranged from 8.0 m s −1 to 22.0 m s −1 . For Hungarian soils (central Europe) on the basis of wind-tunnel experiments, [62] stated the TWSs for soils with sand texture (5.7-8.1 m s −1 ), loamy sand texture (6.1-9.1 m s −1 ), sandy loam texture (7.4-9.8 m s −1 ), loam texture (9.3-11.6 m s −1 ) and silt loam texture (10.0-11.6 m s −1 ). The wind data measured by weather stations at a height of 10 m can be easily interpolated to any height of 0 m to 10 m using models with the estimated parameters mentioned in [63]. This study provides results that were used for the derivation of the wind speed at ground level. It can be derived that wind speeds of 10 m s −1 and 10 m correspond to 3.3 m s −1 at ground level.
Considering two principally different mechanisms of wind erosion due to different soil types, we again distinguish between CE of sandy soils (CE SAND ) and CE of clay soils (CE CLAY ).
Climatic erosivity of sandy soil (CE SAND ) Sandy soils, identified as vulnerable (see the SE SAND chapter), increase their susceptibility to wind erosion with drought in their surface layer. Contrary to clay soils, the SE SAND of sandy soils is not affected by a specific winter regime but is simply given by their texture, i.e., the content of the EF in both seasons with low SE VC and SE SC (spring and autumn) factors are taken into account in terms of the drought in the surface layer. Considering the general comments on the CE CLAY , only the spring period was taken into account in terms of the CE SAND -REW.
The CE SAND -DSL and CE SAND -REW, when joined together according to Matrix 2 ( Table 6), reveal the degrees of CE SAND . Case study on CE SAND CE SAND -DSL To identify drought in the soil-surface layer, we employed CBSS data. This indicator is specified in [45]. With respect to the spring and autumn vulnerability of sandy soils described above, CBSS data (code 0) from spring (March to May) and autumn (September to October) were employed for CE SAND identification. CBSS code 0 is defined as: "soil surface is completely dry but not frozen; soil particles crumble to dust when squeezed, wind picks up dust from soil". The long-term (30yr) average number of days with CBSS code 0 was calculated for 70 stations all over the area of interest. The area of interest was divided accordingly into 5 categories of the CE SAND -DSL.
CE SAND -REW The risk of erosive winds in the spring and autumn was based on 10-min or 15-min wind speed data from 103 stations all over the area of interest for the longest possible (9 years) period in terms of data availability. The percentage of 10/15 min with a wind gust above 10 m s −1 of the total number of measurements (approx. 136,000 values for each station) was determined. The threshold value of the wind speed 10 m s −1 (measured at a height of 10 m above the surface) results from the conclusions of field measurements, when it was found that the ground wind speed is approximately one-third to one-half the wind speed at 10 m. A speed of 3.3 m s −1 is then the minimum towing speed for soil particles on dry light soil. The area of interest was divided accordingly into 5 categories of CE SAND -REW.
The CE SAND degree for the area of the case study is presented in Figure 4.  Case study on CESAND CESAND-DSL To identify drought in the soil-surface layer, we employed CBSS data. This indicator is specified in [45]. With respect to the spring and autumn vulnerability of sandy soils described above, CBSS data (code 0) from spring (March to May) and autumn (September to October) were employed for CESAND identification. CBSS code 0 is defined as: "soil surface is completely dry but not frozen; soil particles crumble to dust when squeezed, wind picks up dust from soil". The long-term (30yr) average number of days with CBSS code 0 was calculated for 70 stations all over the area of interest. The area of interest was divided accordingly into 5 categories of the CESAND-DSL.
CESAND-REW The risk of erosive winds in the spring and autumn was based on 10-min or 15-min wind speed data from 103 stations all over the area of interest for the longest possible (9 years) period in terms of data availability. The percentage of 10/15 min with a wind gust above 10 m s −1 of the total number of measurements (approx. 136,000 values for each station) was determined. The threshold value of the wind speed 10 m s −1 (measured at a height of 10 m above the surface) results from the conclusions of field measurements, when it was found that the ground wind speed is approximately one-third to one-half the wind speed at 10 m. A speed of 3.3 m s −1 is then the minimum towing speed for soil particles on dry light soil. The area of interest was divided accordingly into 5 categories of CESAND-REW.
The CESAND degree for the area of the case study is presented in Figure 4. Climatic erosivity of clay soil (CE CLAY ) In terms of the mechanism, CE CLAY is analogous to CE SAND with one modification. CE CLAY deals with critical period determination. Due to the specific winter regime, which invokes the creation of aggregates removable by wind, the critical period for clay soils in terms of the drought occurrence is mainly early spring. In the spring months, these soils show considerable potential erodibility due to the disintegration of aggregates in the winter, while in the autumn they tend to be mostly compact, and unless they are damaged by inappropriate cultivation they are almost invulnerable to wind erosion. Therefore, only spring droughts were taken into account. Considering the general comments on CE CLAY , only the spring period was taken into account in terms of the CE CLAY -REW.
"Drought in surface layer" and "Risk of erosive wind occurrence" when joined together according to Matrix 2 (Table 6) reveal the degrees of CE CLAY (Table 7). Case study on CE CLAY CE CLAY -DSL and CE CLAY -REW were estimated analogically to CE SAND -DSL and CE SAND -REW. However, with respect to the general definition of CE CLAY , only data from the spring period were used.
The CE CLAY degrees for the area of the case study are presented in Figure 5.

PWE
Overall PWE consists of both PWESAND (Table 8) and PWECLAY (Table 9) Figure 5. Map of the CE CLAY degrees for the area under the case study.

PWE CLAY
Join SE CLAY 1, 3, 5 (see Table 5) with CE CLAY (see Table 7-Matrix 3). PWE SAND and PWE CLAY , when combined, reveal aggregate PWE degrees. The grouping according to Table 10 belongs to one of the risk categories of the PWE. The highest exposure >30 Case study on PWE The resulting PWE combining both PWE SAND and PWE CLAY and respected categorization of the case study area is presented in Figure 6. 3 Moderate exposure (12, 18⟩ 4 High exposure (18, 24⟩ 5 Very high exposure (24, 30⟩ 6 The highest exposure >30 Case study on PWE The resulting PWE combining both PWESAND and PWECLAY and respected categorization of the case study area is presented in Figure 6.

Results
The main goal of the paper was to develop a complex methodological protocol (expressed as "road map"-see Figure 7) to identify and classify the risk of wind erosion given as PWE degrees and their subsequent distribution into six PWE categories (Table  10). Particular steps of the road map creation are gradually and described in detail in the methodology. Figure 7 thus summarizes the whole methodological process of PWE derivation. We also bring the "road map" application to demonstrate the entire method on a particular case study area in central Europe that represents part of the Czech Republic. The indicators employed in the case study to express key soil-climatic elements are summarized in Table 11. The "road map" being applied to the case -study area classifies it as 46% out of exposure and 54% under exposure, while the majority of the area falls into subtle and slight exposure categories of PWE ( Figure 8).

Results
The main goal of the paper was to develop a complex methodological protocol (expressed as "road map"-see Figure 7) to identify and classify the risk of wind erosion given as PWE degrees and their subsequent distribution into six PWE categories (Table 10). Particular steps of the road map creation are gradually and described in detail in the methodology. Figure 7 thus summarizes the whole methodological process of PWE derivation. We also bring the "road map" application to demonstrate the entire method on a particular case study area in central Europe that represents part of the Czech Republic. The indicators employed in the case study to express key soil-climatic elements are summarized in Table 11. The "road map" being applied to the case -study area classifies it as 46% out of exposure and 54% under exposure, while the majority of the area falls into subtle and slight exposure categories of PWE ( Figure 8).

Soil-Climatic Element
Indicator Employed in the Case Study

SESAND
Sandy soil groups based mainly on texture

Discussion
Gradual progress in natural risk assessment, quantification and modelling of water erosion has been connected with development of computer science such as computer simulation, machine learning, artificial intelligence, geospatial technology and remote sensing in recent decades. Enhanced approaches such as concept of connectivity [64], which allows building models with a holistic approach to simulate system dynamics to design better measuring and modelling schemes for water and sediment dynamics are yet to reveal their potential.
On the contrary, advanced machine learning and related techniques have not been sufficiently applied on the issue of wind erosion except for several surveys, for instance [65]. Since wind erosion is actually an area-wide phenomenon its spatial variability and regional specificity cause methodological obstacles of its monitoring and modelling. Physical process-based models, such as the WEQ, the wind erosion prediction system (WEPS), and the revised wind erosion equation (RWEQ), in combination with GIS methods, are frequently used for soil-risk quantification. However, these models have mainly been developed for the agro-ecological conditions of North America. These methods are primarily applicable to relatively small areas that are known to be susceptible to wind erosion [66,67]. Furthermore, compromising methods of wind-risk estimation have been employed in Germany. The standard in question is DIN 19706:2013-02 "Soil quality-Determination of the soil erosion risk caused by wind". The map "Potential risk of wind erosion on arable soils in Germany" is based on pedological (soil category, humus content) and climatic factors (mean annual wind speed at 10 m above ground level); the land-use information is derived from the CORINE land-cover data set. A similar concept

Discussion
Gradual progress in natural risk assessment, quantification and modelling of water erosion has been connected with development of computer science such as computer simulation, machine learning, artificial intelligence, geospatial technology and remote sensing in recent decades. Enhanced approaches such as concept of connectivity [64], which allows building models with a holistic approach to simulate system dynamics to design better measuring and modelling schemes for water and sediment dynamics are yet to reveal their potential.
On the contrary, advanced machine learning and related techniques have not been sufficiently applied on the issue of wind erosion except for several surveys, for instance [65]. Since wind erosion is actually an area-wide phenomenon its spatial variability and regional specificity cause methodological obstacles of its monitoring and modelling. Physical process-based models, such as the WEQ, the wind erosion prediction system (WEPS), and the revised wind erosion equation (RWEQ), in combination with GIS methods, are frequently used for soil-risk quantification. However, these models have mainly been developed for the agro-ecological conditions of North America. These methods are primarily applicable to relatively small areas that are known to be susceptible to wind erosion [66,67]. Furthermore, compromising methods of wind-risk estimation have been employed in Germany. The standard in question is DIN 19706:2013-02 "Soil quality-Determination of the soil erosion risk caused by wind". The map "Potential risk of wind erosion on arable soils in Germany" is based on pedological (soil category, humus content) and climatic factors (mean annual wind speed at 10 m above ground level); the land-use information is derived from the CORINE land-cover data set. A similar concept for wind soil assessment in Iran (vegetation-cover data based on remote sensing from the Moderate Resolution Imaging Spectroradiometer Normalized Difference Vegetation Index-MODIS NDVI, a basic soil texture map and wind speed data from meteorological stations) was used by [68].
The findings of our case study were compared with various surveys on wind erosion previously conducted.
According to [46], who developed unified methodology of the wind risk estimation based on the CE, soil erodibility, vegetation cover and landscape roughness (Index of Land Susceptibility to Wind Erosion-ILSWE) for 36 European countries, the Czech Republic belongs among the countries with low land susceptibility to wind erosion (just 1.8% of the area falls into moderate susceptibility but not into high land susceptibility). Similarly, the analysis of wind-erosion risk based on the number of erosive days by [69] classifies the majority of the Czech Republic as slightly endangered (0.03-1.5 erosive days per year). It well corresponds to our findings given in our results. However, very high values of the erosivity index (mainly determined by the wind velocity and modified by the moisture conditions) and small or medium values of the erodibility index (mainly attributed to the texture and organic matter content) in the climate of the Czech Republic were reported by [70]. Analogically, a map of the distribution of wind-erosion risk within NUTS-x for the EU-27 in 2025 [71] classifies the majority of the Czech area as regions with high or the highest wind-erosion risk (2-5.57 and more than 5.57 t ha −1 year −1 ). An evaluation of topsoil loss by wind was published by [69]. The most extensive and severe wind erosion is mapped in southeastern Europe, Romania, Ukraine and Russia. According to this study, the Czech Republic is the most endangered among all countries of western and central Europe.
Unlike the abovementioned studies, our approach to PWE also encompassed a module of clay soils, which are not generally considered to be vulnerable to wind erosion. That is why they are not being linked to any risk category in terms of the PWE. However, some authors [35,36,41] claim that the combination of particular soil properties with a specific winter regime makes specific clay soils highly vulnerable to wind erosion. The periodic freezing and thawing of soil combined with enormous moisture cause aggregate disintegration. Previously, published and generally used methods of the PWE estimation (even though they are more profound in terms of soil-crust factors, crop factors, etc.) by [30][31][32][33] do not take into account one of the crucial climatic elements of the PWE, i.e., PWE long-term wind conditions. The areas that tend to be windier are logically under a higher risk of wind erosion than calm areas. However, it is important to emphasize that long-term wind conditions, as a key factor of PWE, cannot be mistaken for actual wind and are a key factor of the AWE.

Conclusions
This study contributes to worldwide land anti-degradation efforts and provides new insights into the methodological approach of wind-erosion risk assessment. We identified and employed key soil-climatic elements that make soil vulnerable to wind erosion. Only when all of them are combined it enables wind-erosion risk across various soil and climate environments to be evaluated. We provided a parallel application of the method on the particular case study, which also suggests usable indicators for each element. An interregional comparation is possible, as long as all employed indicators are homogenous across regions. However, following the road map, it is also possible to classify any region in terms of the wind-erosion risk even with the application of different indicators according to the data available. The novelty of our method lies in: (i) all relevant soil-climatic data of wind erosion are combined, (ii) different soil types "sand" and "clay" are evaluated simultaneously with respect to the different mechanisms of wind erosion, and (iii) a methodological road map enables its application for various conditions. To be effective, wind-erosion research requires the explicit articulation of its objectives, which should be interpreted and translated into real actions. Our results provide a useful tool for management decisions. High-quality data and analyses that can support decision making are demanding. Based on our method, it is possible to set threshold values that, when exceeded, trigger landscape adjustments, more detailed in situ measurements, or indicate a need for specific management.
Ongoing and planned future research will attempt to merge the map of the PWE with anti-erosion measures (windbreaks) and revise the PWE accordingly. The resulting GIS-based software would then allow not only an evaluation of the efficiency of existing measures, but also a judgment of various plans of the building of the new ones to accomplish healthy, resilient and sustainable landscapes.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to licensing policy of data providers (Czech Hydrometeorological Institute and Research Institute for Soil and Water Conservation).