Abstract
Predicting land susceptibility to wind erosion is necessary to mitigate the negative impacts of erosion on soil fertility, ecosystems, and human health. This study is the first attempt to model wind erosion hazards through the application of a novel approach, the graph convolutional networks (GCNs), as deep learning models with Monte Carlo dropout. This approach is applied to Semnan Province in arid central Iran, an area vulnerable to dust storms and climate change. We mapped 15 potential factors controlling wind erosion, including climatic variables, soil characteristics, lithology, vegetation cover, land use, and a digital elevation model (DEM), and then applied the least absolute shrinkage and selection operator (LASSO) regression to discriminate the most important factors. We constructed a predictive model by randomly selecting 70% and 30% of the pixels, as training and validation datasets, respectively, focusing on locations with severe wind erosion on the inventory map. The current LASSO regression identified eight out of the 15 features (four soil property categories, vegetation cover, land use, wind speed, and evaporation) as the most important factors controlling wind erosion in Semnan Province. These factors were adopted into the GCN model, which estimated that 15.5%, 19.8%, 33.2%, and 31.4% of the total area is characterized by low, moderate, high, and very high susceptibility to wind erosion, respectively. The area under curve (AUC) and SHapley Additive exPlanations (SHAP) of game theory were applied to assess the performance and interpretability of GCN output, respectively. The AUC values for training and validation datasets were estimated at 97.2% and 97.25%, respectively, indicating excellent model prediction. SHAP values ranged between −0.3 and 0.4, while SHAP analyses revealed that the coarse clastic component, vegetation cover, and land use were the most effective features of the GCN output. Our results suggest that this novel suite of methods is highly recommended for future spatial prediction of wind erosion hazards in other arid environments around the globe.
1 Introduction
Soil erosion by water and wind are the two main processes leading to its degradation (Webb et al., 2016; ; Soloki and Noori Nahad, 2019; Liu et al., 2021; Jia et al., 2022). Soil erosion by wind is an environmental hazard, which threatens land productivity and ecosystems, reduces atmospheric visibility during dust storms, and affects ecosystems, animal, and human health (Gregory et al., 2004). Wind erosion affects ∼28% of degraded land worldwide (Oldeman, 1994; ). However, predicting susceptibility to wind erosion remains a difficult task to be accurately quantified (Webb and Strong, 2011; ). On the other hand, predicting land susceptibility to wind erosion is highly beneficial for mitigating its negative consequences, such as depletion of soil organic carbon, decreased soil depth, reduced soil structure stability in situ, increased dust aerosol loading, degradation of air quality, and increased respiratory diseases off-site (Xu et al., 2018; ; ; Tian et al., 2021; ).
Wind erosion is a serious environmental problem in large areas of Iran, especially in the central, eastern, and southern parts of the country (; ). It is estimated that an area of 23.8 million ha is affected by wind erosion, of which 51% is susceptible to aeolian entrainment, 31% experiences significant aeolian transport, and 18% represents the area dominated by sedimentation (). The frequency of dust storms has increased substantially across Iran over the last 20 years, exacerbated by prolonged dry conditions during 2000–2002 and 2008–2010 (; ; Hamzeh et al., 2021; Rashki et al., 2021). In order to mitigate the effects of wind erosion and associated dust storms on the atmosphere, geosphere, and biosphere, the prediction and mapping of land susceptibility to wind erosion are essential (; Jafari et al., 2022; Pourhashemi et al., 2022).
Wind erosion depends on a range of variables such as vegetation cover, land surface characteristics, and climatic parameters (; Parajuli et al., 2014; Padarian et al., 2020). The contribution of each variable to erosion susceptibility remains poorly understood and is highly region-specific, while it is necessary to be defined in order to effectively model the spatial distribution of wind erosion hazards (; Kannangara et al., 2022). Feature selection (FS) is a critical step in the application of data mining (DM) models (Kursa and Rudnicki, 2010). Various FS algorithms, such as Boruta, leave one feature out (LOFO), and DM models (e.g., random forest, tree-based genetic algorithms, multivariate adaptive regression spline (MARS), CForest, and Cubist), have been widely used to select key controlling factors for environmental hazards such as flood, land subsidence, landslides, and dust emissions (Hamza et al., 2019; ).
Data mining (DM) is a component of data science (DS), which has frequently been widely used for spatial predictions in environmental and hydrological sciences (; ; ; Saha et al., 2021; ; ; Saha et al., 2022). Over the past few years, different types of DM models (e.g., ANFIS, frequency ratio, weights of evidence, Cubist, CForest, and extreme gradient boosting) have been applied to map land susceptibility to wind erosion and dust emissions in various parts of Iran, vulnerable to wind erosion, desertification, and climate change (; ; ). Furthermore, recurrent neural network (RNN), restricted Boltzmann machine (RBM), gcForest, bidirectional long short-term memory (Bi-LSTM), and Gaussian copula-based multivariate model, as deep learning (DL) models, have been used for the spatial modeling of dust provenance over the globe (; ). DL modeling is specifically interesting because this approach overcomes many of the restrictions placed on the modeling process by traditional machine learning (ML) or shallow models. In comparison with the traditional or shallow models, different architectures of DL transform the data presentation at one level into a representation at a higher level (LeCun et al., 2015). DL is increasingly being used as a tool to investigate a range of scientific problems, including spatial modeling of environmental issues (; Mohammadifar et al., 2021b; ; Mohammadifar et al., 2022). Typical DL models can be divided into four classes, namely, stacked auto-encoder (SAE), deep belief network (DBF), recurrent neural network (RNN), and convolutional neural network (CNN) (Zhang et al., 2018), of which CNN is one of the most popular for prediction goals (Reiss and Bogenberger, 2016; Zhang et al., 2016). One branch of CNN, the convolutional graph neural networks (ConvGNNs), has been applied to a wide range of research topics including computer vision (Monti et al., 2017), natural language processing, traffic prediction (Li et al., 2018), social analysis (), bioinformatics, and fraud detection (; Zhang et al., 2017; Zhang et al., 2019), but has not yet been investigated for environmental issues. The ConvGNN approach shows promising results for environmental hazard prediction due to the advantage of graphs to capture structural relations between data over analysis of data in isolation. Overall, ConvGNNs (Zhang et al., 2019) can be classified into the spectral-based techniques (; Kipf and Welling, 2016) and spatial-based techniques (Monti et al., 2017; Veličković et al., 2017). In addition, aspects of game theory may be used to assess the interpretability of predictive DM and DL models. Previous environmental studies have successfully made use of SHapley Additive exPlanations (SHAP) from game theory to interpret predictive DM models for digital soil mapping and soil erosion by water (Padarian et al., 2020; Mohammadifar et al., 2021a). These early results indicate that it may be beneficial to explore game theory through SHAP, in order to assess the interpretability of our GCN model outputs for wind erosion.
This study aimed to investigate novel tools for the prediction of susceptibility to wind erosion on the central Iranian Plateau (Semnan Province). To the best of our knowledge, this work is the first attempt at mapping wind erosion hazards by combining LASSO regression, graph convolutional networks (GCNs), and game theory. Game theory (SHAP) was applied to assess the interpretability of the DL results. The main objectives and goals of this study were as follows: 1) the identification of the most important factors controlling wind erosion for the study area using LASSO regression, 2) the spatial prediction of land susceptibility to wind erosion by applying a GCN deep learning model with Monte Carlo dropout, 3) the evaluation of the GCN model using the receiver-operating characteristic area under curve (ROC-AUC) approach, and 4) the examination of the interpretability of our predictive GCN model using game theory. We assess, on the basis of our case study, whether this approach could be applied for spatial mapping of dryland hazards such as wind erosion and dust plume transport worldwide.
2 Materials and methods
2.1 Study area
Semnan Province is a region covering 97,322 km2 in the northern part of the central Iranian Plateau, lying on the southern margins of the Alborz mountain range (Figure 1). The region is hyper-arid to semi-arid, with rough terrain and low population density, covering the northern part of the Kavir Desert. The altitude of the study area ranges substantially from 658 m in the south and central parts to 3,885 m in the north. Mean annual precipitation follows a steep gradient from 75 mm in the hyper-arid south to 699 mm in the more mountainous and vegetated northern part. The topography and climatic conditions contribute to significant variability in soil organic carbon stock (OCS) contents across the province, with values ranging from near zero in the central and southern parts to 78% in the north.
FIGURE 1
2.2 Mapping factors controlling wind erosion
Initially, maps for 15 potential factors controlling wind erosion in Semnan Province were generated (S1 and S2). The potential factors were selected based on previous studies investigating land susceptibility to wind erosion in Iran (; ; ). These factors include topography, which is expressed through a digital elevation model (DEM), and spatial distribution of precipitation, evaporation, surface wind speed, normalized difference vegetation index (NDVI), soil bulk density, nitrogen and organic carbon stock (OCS), surface cation exchange capacity (CEC), proportion of surface sediment of clay, silt, sand, and percentage of the coarse fragment, lithology, and land use. The topographical data for the DEM were downloaded from https://earthexplorer.usgs.gov with a spatial resolution of 30 m × 30 m. Climatic/meteorological data (precipitation, evaporation, and wind speed) were derived and interpolated from measurements at 27 synoptic meteorological stations in Semnan Province and surrounding areas averaged over the period 2012–2019, and mapped in ArcGIS. The NDVI data were retrieved for June 2020 using the MODIS (Moderate Resolution Imaging Spectroradiometer) satellite sensor on the assumption that June is representative for average annual vegetation cover in Semnan Province. We assume that since vegetation in the province is sparse and comprises bushes, shrubs, and minimal tree cover, variability in NDVI values is minimal over the spring and summer months in Semnan Province. Soil characteristics were obtained from publicly available data (https://soilgrids.org) and mapped with a spatial resolution of 50 m × 50 m. The geology of Semnan Province is characterized by a number of lithological units ranging in age from Paleozoic through to Quaternary. In this study, we used existing geological maps generated by the Iranian Geological Survey and Mineral Exploration department. Furthermore, we incorporated land use maps for the year 2020 from the Iran Forest, Rangeland, and Watershed Management Organization (IFRWMO). Since the spatial resolution of potential factors in various constructed maps is not uniform, all input layers used for modeling were converted to a consistent spatial resolution of 50 × 50 m. The details of input data used in the modeling process are presented in Table 1. It should be noted that most of the variables here are used as controls for the topographical, land use, and soil characteristics of the examined area that are fixed and not changed along season.
TABLE 1
| Data type | Source/spatial resolution/scale |
|---|---|
| DEM | USGS/30 m × 30 m |
| Climatic data | 27 synoptic meteorological stations/annual |
| Soil data | https://soilgrids.org//250 m × 250 m |
| NDVI | MODIS satellite product (MOD13Q1)/250 m × 250 m |
| Lithology map | Geological Survey and Mineral Exploration of Iran/1:200,000 |
| Land use map | Iran Forest, Rangeland, and Watershed Management Organization/1:200,000 |
Details of input data used in the modeling process.
2.3 Inventory map of wind erosion
We then developed an inventory map for predicting land susceptibility to wind erosion () over the study region, using a previously generated map of the regions for which active wind erosion has been monitored from the IFRWMO. Furthermore, 70% and 30% of the pixels (80 pixels and 33 pixels, respectively) in these regions (Figure 1) were randomly selected for training and test (validation) datasets, to build a new predictive model. We additionally selected 80% (training) and 20% (validation) of the training dataset for our predictive models.
3 Mapping wind erosion hazards using LASSO regression, deep learning, and game theory
In Figure 2, we summarize our workflow for mapping of land susceptibility and wind erosion hazards in Semnan Province, including feature selection by LASSO regression, the GCN deep learning model, and game theory-based interpretation. This workflow may be applied in other dryland environments and for spatial mapping of various environmental hazards around the globe. Each step is detailed in the following subsections.
FIGURE 2
3.1 Feature selection by LASSO regression
Feature selection (FS) describes the process of narrowing down the potential factors (explanatory features or variables), which can be used to describe a response feature (). In this case, the response feature is the susceptibility to wind erosion and the explanatory features are the 15 parameters (potential factors) mapped in Supplementary Figures S1, S2.
The most influential factors controlling wind erosion in our study area were selected using the LASSO regression, a type of FS algorithm proposed by Tibshirani (1996). LASSO is considered to be the optimal approach for FS in our study since it combines the qualities of ridge regression and subset selection, thus improving the model’s accuracy and interpretability. LASSO is a regression analysis that performs both variable selection and regularization in order to enhance the prediction accuracy and interpretability of the resulting model, which is a major reason for its selection against other regression approaches. Furthermore, LASSO is a forward-looking variable selection method for regression, which is a type of penalized least-squares approach incorporating L1-penalty function (Muthukrishnan and Rohini, 2016). The LASSO approach transforms each coefficient by a constant component λ, truncating at 0, with the objective of minimizing prediction errors (). Therefore, LASSO decreases the residual sum of squares, subject to the sum of the absolute value of the coefficient being less than a constant. During the application of FS, the features (variables) that still had a non-zero coefficient after the shrinking process were selected to be part of our model. The LASSO can be written as follows:where are the constant coefficient and the coefficient vector, respectively, indicates the outcome, and = ( is the covariate vector for the ith case.
3.2 Graph convolutional networks
Here, we undertake the first known study using convolutional graph neural networks (ConvGNNs) to investigate wind erosion susceptibility. The first adoption of neural network approaches to directed acyclic graphs was reported by Sperduti and Starita (1997). Neural networks have since been more widely applied in the form of graph neural networks (GNNs) of various kinds for a wide range of applications (Gori et al., 2005; Scarselli et al., 2009; ). Among these, graph convolutional networks (GCNs; Figure 3), first proposed by Kipf and Welling (2016), are a type of spectral-based convolutional GNN (Wu et al., 2020), which introduce a first-order approximation of Chebyshev spectral CNNs (hereafter ChebNet) ().
FIGURE 3
In this study, we train GCNs using a Monte Carlo dropout (). Assuming K = 1 and λmax = 2, the equation of is simplified as follows:where, K, G,Θ, and X indicate the layer index, a graph, the learnable model variable, and the feature vector, respectively.
To restrain the number of variables and avoid over-fitting, we assumed θ = θ0 = −θ1 for our GCNs, leading to the following definition of graph convolution:
To allow multiple channels of inputs and outputs, our GCN modifies Eq. 3 into a compositional layer as follows:
where and f (.) indicates an activation function. Applying empirically causes numerical instability to GCNs (Wu et al., 2020). To solve this problem, GCNs use a normalization trick to replace with and .Here, is the learnable model variable, while A and D indicate the graph adjacency matrix and the degree matrix of .
The predicted susceptibility values for all the pixels ranged between 0 and 1, with higher values indicating higher susceptibility to wind erosion. We divided susceptibility predictions into four classes: low (0–0.25), moderate (0.25–0.50), high (0.50–0.75), and very high (0.75–1).
3.3 Assessing GCN performance to spatial modeling of wind erosion
Previous studies have used several statistical indicators to assess the reliability and feasibility of predictive DM and DL models in environmental challenges (; ), including root-mean-square error (RMSE), mean absolute error (MAE), mean bias error (MBE), the Nash–Sutcliffe efficiency coefficient (NSC), Taylor diagram, and the receiver-operating characteristic–area under curve (ROC-AUC). In this study, we applied the ROC-AUC (Hong et al., 2017) to evaluate the developed GCN model for spatial prediction of wind erosion susceptibility in Semnan Province. Using the area under curve (AUC) value, we classified the accuracy of the model prediction into five categories: poor (50%–60%), moderate (60%–70%), good (70%–80%), very good (80%–90%), and excellent (90%–100%) (Tibshirani, 1996).
3.4 Interpretability of GCNs using game theory
In addition, we interpreted our GCN model for predicting wind erosion hazards using the game theory approach of SHapley Additive exPlanations (SHAP) (Lundberg and Lee, 2017; Kannangara et al., 2022). SHAP values indicate the contribution of a given factor to the final predictions of a model. SHAP specifies the explanation as follows:where g, ∈ (), and M are the explanation model, the coalition vector, and the maximum coalition size, respectively.
For calculating SHapley values, we fit a model fSU(i), involving factor i, and another additional model, fs. All possible subsets S ⊆ F comprise the empty set and the set F of all factors. The final factor contribution, ∈ R, is the weighted average of all marginal contributions (Padarian et al., 2020), which is given as
4 Results and discussion
4.1 Discriminating important from non-important features influencing wind erosion
The results of the LASSO regression, which was used to select the important factors responsible for wind erosion, are presented in Figure 4; Table 2. Figure 4 shows that AUC values decrease with the increasing number of the potential factors controlling wind erosion, which indicates that the model’s accuracy increases when more potential factors are considered within the LASSO regression.
FIGURE 4
TABLE 2
| Potential factor | Value | Decision |
|---|---|---|
| Soil nitrogen | 0 | Reject |
| OCS | −0.04 | Accept |
| Precipitation | 0 | Reject |
| Sand content | 0 | Reject |
| Silt content | 0 | Reject |
| Wind speed | −0.17 | Accept |
| Bulk density | 0.01 | Accept |
| CEC | 0.01 | Accept |
| Clay content | 0 | Reject |
| Coarse fragment | −0.01 | Accept |
| DEM | 0 | Reject |
| Evaporation | 0.04 | Accept |
| Lithology | 0 | Reject |
| Land use | 0.01 | Accept |
| NDVI | −5.03 | Accept |
Values of the potential factors controlling wind erosion, as estimated by the least absolute shrinkage and selection operator (LASSO) regression.
Out of the 15 features examined as potential factors controlling wind erosion, seven were identified as non-important in our study area and were excluded from further analysis. These non-important factors comprise soil nitrogen, precipitation, sand, silt, and clay content; topography (DEM); and lithology.
The LASSO regression identifies precipitation as a non-important factor for land susceptibility to wind erosion hazards in Semnan Province. Precipitation is generally considered an important factor controlling aeolian erosion and dust storms since low precipitation has a positive feedback effect for reducing soil moisture and vegetation growth, thereby enhancing erosion by wind (Sissakian et al., 2013; Parajuli et al., 2014). In the case of Semnan, the mountainous areas within the province receive relatively high rainfall and have a correspondingly lower vulnerability to aeolian erosion. This context reduces the relative importance of precipitation for dust emission in the study region, highlighting the sensitivity of local factors in the LASSO regression process. This suggests that in other arid regions, the current non-important factors may be highly sensitive to land susceptibility to wind erosion (; ; ).
Although percentage fractions of unconsolidated silt, clay, and sand often influence dust emission and transport (), our LASSO regression classified these factors as non-important for Semnan Province. Clearly, then, the fact that our LASSO regression finds typically influential factors to be non-important for our study region indicates that the classification of the importance of potential factors is highly region-specific and depends on local geographical, land cover, and meteorological conditions. Important factors are likely to depend on the type of soil, the dimension of desert areas, and topography.
LASSO criteria identified eight features as important variables for wind erosion in Semnan Province: wind speed, evaporation, OCS, bulk density, cation exchange capacity (CEC), coarse fragment, land use, and NDVI. Only these eight features entered into the spatial modeling process using GCNs to predict land susceptibility to wind erosion hazards in Semnan Province.
4.2 Map of wind erosion hazards
Figure 5 shows the map of susceptibility to wind erosion in Semnan Province, as produced by the GCN deep learning model using the eight important factors (Figure 4; Table 2).
FIGURE 5
According to the results, a substantial proportion of the region (c. 65%) exhibits high or very high susceptibility to wind erosion, which is 33.2% (= 32,300 km2) and 31.4% (= 30,581 km2), respectively. Areas classified as having low and moderate susceptibility to wind erosion represent 15.5% (= 15,136 km2) and 19.8% (= 19,305 km2) of the total area in Semnan Province, respectively.
The southern and central parts of the study area are classified as having very high susceptibility to wind erosion (Figure 5) and, therefore, suffer from extreme soil erosion hazards. These areas lie within the Kavir Desert, which is known as a generation point (dust source) for dust storms impacting the whole central Iranian Plateau (; Marjovvi et al., 2022). The Kavir Desert experiences high evaporation rates, low surface soil moisture, and high wind speeds; it contains only sparse vegetation cover and has low soil OCS and low CEC. Furthermore, these regions contain substantial unconsolidated fine-grained sand at the surface, sourced by alluvial input from inland streams. Dry river beds and ephemeral lakes have long been recognized as areas with high potential for wind erosion and dust emissions across the Middle East (Rashki et al., 2017). Such areas, with high sediment availability and frequent high wind speeds, often become the major sources of erodible sand for dust/sand storms in Semnan Province.
The part of Semnan Province with the lowest susceptibility to wind erosion is located in the north within the Alborz mountain range. This area (c. 35.5% of the province) is characterized by steep slopes, relatively high vegetation cover, and low evaporation rates, as well as by high proportions of silt-sized sediment and coarse fragments on the surface (Supplementary Figures S1, S2). Each of these factors represents negative feedback for wind erosion (Parajuli et al., 2019) and is likely to contribute to the GNS assessment of low susceptibility in our study. Although this part of Semnan Province was found to have the lowest susceptibility to wind erosion, it is nevertheless the part of the Iranian Plateau that is more susceptible to wind erosion than other areas (; ).
Figure 6 summarizes the results of sensitivity analysis for our GCN deep learning model. We found that the most sensitive features for the model are surface coarse fragment cover, evaporation, and land use; these factors exhibit the highest slopes on the graph and the largest variation in the normalized predictor change. The high importance of these features for soil susceptibility to wind erosion indicates strong regional dependence specific for Semnan Province. Using the same deep learning methods, we would likely find the deciding factors to be different in other arid environments (; ; ).
FIGURE 6
Overall, wind erosion and its contribution to desertification and land degradation is controlled by the arid climate, high evaporation rates, low precipitation, sparse vegetation cover, and intense winds (Keramat et al., 2011) in Semnan Province and in central Iran, more broadly. However, the balance and drivers of the determining factors—demonstrated here to be regionally specific—may be altered by future climatic change, including extensive droughts and increasing average temperatures, and by anthropogenic degradation of ecosystems and mismanagement of water, soil, and plant resources over the whole Middle East region (; Middleton, 2017; Middleton, 2019; ; Papi et al., 2022). On this basis, continued monitoring of the relative contributions of various controlling factors to dynamic phenomena such as desertification and dust storms in arid environments is recommended.
4.3 Assessment of GCN model performance
The performance of our GCN modeling results was assessed using the area under curve (AUC) approach on the training and test datasets (Figure 7). The calculated AUC values for the training and validation datasets were found to be the same (97.2% and 97.25%, respectively), thus indicating an excellent prediction accuracy (Yesilnacar, 2005). The results show that GCNs significantly improve the performance of the predictions for land susceptibility to wind erosion and support the findings of previous studies, which suggested that GCNs are an efficient graph model for semi-supervised learning () and constitute a useful tool for completing missing data (). The DL hybrid models (e.g., convolutional neural network-gated recurrent unit (CNN-GRU) and dense layer deep learning–random forest (DLDL-RF)) were also recently shown to be efficient methods for classifying dust sources in the Middle East ().
FIGURE 7
4.4 Relative contributions of the most important factors for wind erosion hazards
The SHAP values of the game theory were used to determine the contribution of each feature to our GCN model output for land susceptibility to wind erosion. We found that the contributions of the variables followed the order: coarse fragments > NDVI > land use > SOM > wind speed > CEC > bulk density > evaporation. Figure 8 illustrates the relative contributions of these parameters based on SHAP values and combines the feature importance with feature effects. The value, and therefore the relative importance of each parameter, is given by the color scheme. Overlapping points are projected in the y-axis to give a sense of the distribution of the SHAP values per feature, as per the published approach ().
FIGURE 8
Land susceptibility to wind erosion has previously been shown to be driven by a wide range of topographic and bio-climatic variables including surface roughness (which equates to our coarse fragment parameter), mineralogy, vegetation cover (comparable with NDVI), soil characteristics (e.g., soil texture, bulk density, and organic matter), wind speed, precipitation, and evaporation (Goudie and Middleton, 2006; Shao, 2008). Depending on the meteorological, topographical, and soil conditions in the study area, factors controlling the wind erosion may be highly variable from place to place. For example, reported that DEM-extracted variables are the most important factors that control wind erosion in the Jazmurian Basin, while in Central Asia, other parameters like precipitation, bulk density, slope, and lithology (among others) mostly controlled land susceptibility to wind erosion using the LOFO algorithm (). Precipitation can act as a roughness factor for controlling wind erosion in arid and semi-arid regions, while lithology can be another effective variable. Geological formations such as Qsd (sand dune unit) or Qal (Quaternary alluvial deposits) are the most susceptible units against wind erosion because they have single-grain structures with loose particles (Li et al., 2022; Li et al., 2023). The current results largely support these arguments while additionally identifying the main contributing factors and their relative contributions using novel techniques. Our novel approach also predicts the spatial variability of wind erosion hazards and enables the assessment of model performance.
Another important finding here is that factor contributions are strongly regionally controlled but are scale-dependent. A previous study focusing on Isfahan Province, immediately southwest of our study area (), based on the monotone multi-layer perception neural network (MMLPNN), found that topography (DEM), precipitation, and NDVI were the three most important factors controlling wind erosion and, therefore, contrasts with current results. Therefore, the characteristics and relative importance of potential controlling factors are highly region-specific, even for neighboring provinces on the central Iranian Plateau. Over larger scales, however, we start to see similarities in the major factors controlling wind erosion. Analyses using SHAP values for Central Asia—a similar approach to the methods used in this study—revealed that surface coarse fragment concentration, soil bulk density, and precipitation were the most important factors controlling land susceptibility to wind erosion over a much larger area than that investigated here (). Overall, the comparison with previous works that examined land susceptibility to wind erosion over arid areas in Iran and Central Asia revealed significant similarities in the factors that control wind erosion, despite some regional differences, and indicated the importance of defining scale in these sorts of studies.
We found that neural network techniques are a useful and accurate tool for investigating the spatial variability and intensity of land susceptibility to wind erosion (Luo et al., 2020). Our approach may be applied not only to mapping wind erosion hazards but also to a range of other environmental problems requiring spatial characterization, such as soil characteristics and the interaction between atmospheric circulation and the land surface. Our results provide an efficient means to measure soil erosion hazards, which can then be used for assessments of weather systems, ecosystem restoration, human health (Xiao et al., 2017; Wu et al., 2020; ; Tian et al., 2021), land management, and hazard mitigation in arid areas (Sharifikia, 2013; Xu et al., 2019). Policymakers, as end users of results regarding mapping of land susceptibility and wind erosion, can use maps that are produced in this study, to alert residents in high-risk areas for devastating phenomena like severe dust storms and apply appropriate mitigating strategies against natural hazards.
4.5 Advantages and limitations to our DL approach
Data science (DS) is an emerging research field that helps us better understand the complex mechanisms behind environmental phenomena (). DS has several advantages and can support a) systematic and objective exploration and visualization of data, b) improving data quality, c) discovery of patterns contained in large time series, and d) producing new validated and transferable knowledge (). DL models—a fully data-driven approach—are especially interesting in environmental issues because they overcome many of the restrictions placed on the modeling process by traditional machine learning models (LeCun et al., 2015; Van Dao et al., 2020).
The accuracy of the input data to the modeling process is specifically important for generating more accurate and reliable predictive maps for environmental applications. Additional accuracy may be achieved with higher spatial resolution of the imagery such as in NDVI, digital elevation, and soil mapping that would lead to improvements in spatial mapping. Due to the desert and remote characteristics of Semnan Province, meteorological stations are rather scarce with unsuitable spatial coverage. Therefore, interpolation of the meteorological data was associated with higher uncertainty. Meteorological variables (wind speed and rainfall) are also sensitive to the length of instrumental data records (Li et al., 2018). Detailed analysis of various natural and anthropogenic factors affecting land susceptibility to wind erosion and dust storms is crucial for efficient management and mitigation strategies of these phenomena in central Iran.
5 Conclusion
This study is the first attempt to apply GCNs to aeolian geomorphology with the aim of predicting land susceptibility to wind erosion. Our novel methodology comprises feature selection using LASSO regression, the prediction of wind erosion hazards by GCNs, and the assessment of the accuracy and elucidation of major contributing parameters identified in our GCN model output using game theory.
The case study region is Semnan Province in the central Iranian plateau. The applied methodology showed that ∼65% of the total area in Semnan Province is prone to high and very high wind erosion hazards, rendering the area as one of the most susceptible to wind erosion parts of the arid central Iranian Plateau. We also found that the main factors contributing to erosion susceptibility are not necessarily the most typical since wind erosion is largely controlled by the proportion of surface coarse fragments, vegetation cover, and land use, as well as by evaporation, wind speed, bulk density, and cation exchange capacity.
The game theory was used for the first time to test the performance and interpretability of the GCN model, using the area under curve (AUC) and SHapley Additive exPlanations (SHAP). The assessment showed that our GCN approach reliably characterized land susceptibility to wind erosion. The SHAP values showed that proportions of coarse fragments, NDVI, land use, and wind speed contributed mostly to land degradation and wind erosion in Semnan Province.
Our methodology can be used to predict the hazards of wind erosion in other arid environments around the globe and may be extended to address several environmental issues from small to large scales.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
HG: formal analysis, investigation, visualization, writing—original draft, supervision, project administration, and writing—review and editing. AM: software and formal analysis. KF: visualization, writing—original draft, and review and editing. YL: formal analysis and investigation. DK: visualization, writing—original draft, and review and editing.
Acknowledgments
The authors would like to thank the Faculty of Agriculture and Natural Resources, University of Hormozgan, Iran for supporting this joint research project.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenvs.2023.1187658/full#supplementary-material
References
1
AhmadiH. (1999). “Applied geomorphology,” in Water erosion (Tehran: Tehran University Press).
2
AkogluL.TongH.KoutraD. (2015). Graph based anomaly detection and description: A survey. Data Min. Knowl. Discov.29, 626–688. 10.1007/s10618-014-0365-y
3
AmiraslaniF.DragovichD. (2011). Combating desertification in Iran over the last 50 years: An overview of changing approaches. J. Environ. Manag.92 (1), 1–13. 10.1016/j.jenvman.2010.08.012
4
ArabameriA.Seyed DaneshA.SantoshM.CerdaA.Chandra PalS.GhorbanzadehO.et al (2022). Flood susceptibility mapping using meta-heuristic algorithms. Geomatics, Nat. Hazards Risk13 (1), 949–974. 10.1080/19475705.2022.2060138
5
BackstromL.LeskovecJ. (2011). “Supervised random walks: Predicting and recommending links in social networks,” in Proceedings of the fourth ACM international conference on Web search and data mining.
6
Bagheri-BodaghabadiM.JafariM. (2022). The dust deposition model (DDM): An empirical model for monitoring dust deposition using meteorological data over the Isfahan province in central Iran. CATENA211, 105952. 10.1016/j.catena.2021.105952
7
BelghaddarY.ChahinianN.SeriaiA.BegdouriA.AbdouR.DelenneC. (2021). Graph convolutional networks: Application to database completion of wastewater networks. Water13 (12), 1681. 10.3390/w13121681
8
BoroughaniM.PourhashemiS.HashemiH.SalehiM.AmirahmadiA.AsadiM. A. Z.et al (2020). Application of remote sensing techniques and machine learning algorithms in dust source detection and dust source susceptibility mapping. Ecol. Inf.56, 101059. 10.1016/j.ecoinf.2020.101059
9
BoroughaniN.PourhashemiS.GholamiH.KaskaoutisD. G. (2022). Predicting of dust storm source by combining remote sensing, statistic-based predictive models and game theory in the Sistan watershed, southwestern Asia. J. Arid Land13 (11), 1103–1121. 10.1007/s40333-021-0023-3
10
BorrelliP.PanagosP.BallabioC.LugatoE.WeynantsM.MontanarellaL. (2014). Towards a pan-European assessment of land susceptibility to wind erosion. Land Degrad. Dev.27, 1093–1105. 10.1002/ldr.2318
11
ChenJ.MaT.XiaoC. (2018). Fastgcn: Fast learning with graph convolutional networks via importance sampling. arXiv preprint arXiv:1801.10247.
12
ChenW.LeiX.ChakraborttyR.PalS. C.SahanaM.JanizadehS. (2021). Evaluation of different boosting ensemble machine learning models and novel deep learning and boosting framework for head-cut gully erosion susceptibility. J. Environ. Manag.284, 112015. 10.1016/j.jenvman.2021.112015
13
ChowdhuriI.PalS. C.ChakraborttyR.MalikS.DasB.RoyP.et al (2021a). Spatial prediction of landslide susceptibility using projected storm rainfall and land use in Himalayan region. Bull. Eng. Geol. Environ.80, 5237–5258. 10.1007/s10064-021-02252-z
14
ChowdhuriI.PalS. C.SahaA.ChakraborttyR.RoyP. (2021b). Evaluation of different DEMs for gully erosion susceptibility mapping using in-situ field measurement and validation. Ecol. Inf.65, 101425. 10.1016/j.ecoinf.2021.101425
15
ChowdhuriI.PalS. C.SahaA.ChakraborttyR.RoyP. (2022). Mapping of earthquake hotspot and coldspot zones for identifying potential landslide hotspot areas in the Himalayan region. Bull. Eng. Geol. Environ.81 (7), 257. 10.1007/s10064-022-02761-5
16
DefferrardM.BressonX.VandergheynstP. (2016). Convolutional neural networks on graphs with fast localized spectral filtering. Adv. neural Inf. Process. Syst.29.
17
DuniwayM. C.PfennigwerthA. A.FickS. E.NaumanT. W.BelnapJ.BargerN. N. (2019). Wind erosion and dust from US drylands: A review of causes, consequences, and solutions in a changing world. Ecosphere10 (3), e02650. 10.1002/ecs2.2650
18
EffrosynidisD.ArampatzisA. (2021). An evaluation of feature selection methods for environmental data. Ecol. Inf.61, 101224. 10.1016/j.ecoinf.2021.101224
19
EmamianA.RashkiA.KaskaoutisD. G.GholamiA.OppC.MiddletonN. (2021). Assessing vegetation restoration potential under different land uses and climatic classes in northeast Iran. Ecol. Indic.122, 107325. 10.1016/j.ecolind.2020.107325
20
FontiV.BelitserE. (2017). Feature selection using lasso. VU Amsterdam Res. Pap. Bus. Anal.30, 1–25.
21
ForootanE. (2022). Erosion susceptibility assessment using fuzzy logic and multi-influencing factors combination approach. Arabian J. Geosciences15 (5), 444. 10.1007/s12517-022-09598-y
22
GalY.GhahramaniZ. (2016). “Dropout as a bayesian approximation: Representing model uncertainty in deep learning,” in In: International conference on machine learning (PMLR).
23
GallicchioC.MicheliA. (2010). “Graph echo state networks,” in The 2010 international joint conference on neural networks (IJCNN) (IEEE).
24
GholamiH.MohamadifarA.CollinsA. L. (2020b). Spatial mapping of the provenance of storm dust: Application of data mining and ensemble modelling. Atmos. Res.223, 104716. 10.1016/j.atmosres.2019.104716
25
GholamiH.MohamadifarA.GolzariS.KaskaoutisD.CollinsA. L. (2021a). Using the Boruta algorithm and deep learning models for mapping land susceptibility to atmospheric dust emissions in Iran. Aeolian Res.50, 100682. 10.1016/j.aeolia.2021.100682
26
GholamiH.MohamadifarA.RahimiS.KaskaoutisD. G.CollinsA. L. (2021b). Predicting land susceptibility to atmospheric dust emissions in central Iran by combining integrated data mining and a regional climate model. Atmos. Pollut. Res.12 (4), 172–187. 10.1016/j.apr.2021.03.005
27
GholamiH.MohamadifarA.SorooshianA.JansenJ. D. (2020c). Machine-learning algorithms for predicting land susceptibility to dust emissions: The case of the Jazmurian Basin, Iran. Atmos. Pollut. Res.11 (8), 1303–1315. 10.1016/j.apr.2020.05.009
28
GholamiH.MohammadifarA.BuiD. T.CollinsA. L. (2020a). Mapping wind erosion hazard with regression-based machine learning algorithms. Sci. Rep.10 (1), 20494–20516. 10.1038/s41598-020-77567-0
29
GholamiH.MohammadifarA.MalakootiH.EsmaeilpourY.GolzariS.MohammadiF.et al (2021c). Integrated modelling for mapping spatial sources of dust in central Asia - an important dust source in the global atmospheric system. Atmos. Pollut. Res.12, 101173. 10.1016/j.apr.2021.101173
30
GholamiH.MohammadifarA. (2022). Novel deep learning hybrid models (CNN-gru and DLDL-RF) for the susceptibility classification of dust sources in the Middle East: A global source. Sci. Rep.12 (1), 19342. 10.1038/s41598-022-24036-5
31
GibertK.IzquierdoJ.Sànchez-MarrèM.HamiltonS. H.Rodríguez-RodaI.HolmesG. (2018). Which method to use? An assessment of data mining methods in environmental data science. Environ. Model. Softw.110, 3–27. 10.1016/j.envsoft.2018.09.021
32
GilletteD. A. (1979) Environmental factors affecting dust emissions by wind erosion in saharan dust. New York: John Wiley. Morales, C.
33
GoriM.MonfardiniG.ScarselliF. (2005). “A new model for learning in graph domains,” in Proceedings. 2005 IEEE international joint conference on neural networks (IEEE). 2005.
34
GoudieA. S.MiddletonN. J. (2006). Desert dust in the global system. Springer Science and Business Media.
35
GregoryJ. M.WilsonG. R.SinghU. B.DarwishM. M. (2004). Team: Integrated, process-based wind-erosion model. Environ. Model. Softw.19 (2), 205–215.
36
HamzaK.SabineG.WillieG. H. (2019). Digital mapping of soil carbon fractions with machine learning. Geoderma339, 40–58. 10.1016/j.geoderma.2018.12.037
37
HamzehN. H.KaskaoutisD. G.RashkiA.MohammadpourK. (2021). Long-term variability of dust events in southwestern Iran and its relationship with the drought. Atmosphere12, 1350. 10.3390/atmos12101350
38
HongH.PradhanB.BuiD. T.XuC.YoussefA. M.ChenW. (2017). Comparison of four kernel functions used in support vector machines for landslide susceptibility mapping: A case study at suichuan area (China). Geomatics, Nat. Hazards Risk8 (2), 544–569. 10.1080/19475705.2016.1250112
39
JafariR.AmiriM.AsgariF.TarkeshM. (2022). Dust source susceptibility mapping based on remote sensing and machine learning techniques. Ecol. Inf.72, 101872. 10.1016/j.ecoinf.2022.101872
40
JiaL.YuK.-x.LiZ.-b.LiP.ZhangJ.-z.WangA.-n.et al (2022). Temporal and spatial variation of rainfall erosivity in the Loess Plateau of China and its impact on sediment load. CATENA210, 105931. 10.1016/j.catena.2021.105931
41
KannangaraK.ZhouW.DingZ.HongZ. (2022). Investigation of feature contribution to shield tunneling-induced settlement using Shapley additive explanations method. J. Rock Mech. Geotechnical Eng.14 (4), 1052–1063. 10.1016/j.jrmge.2022.01.002
42
KeramatA.MarivaniB.SamsamiM. (2011). Climatic change, drought and dust crisis in Iran. Int. J. Geol. Environ. Eng.5 (9), 472–475.
43
KipfT. N.WellingM. (2016). Semi-supervised classification with graph convolutional networks. arXiv preprint arXiv:1609.02907.
44
KursaM. B.RudnickiW. R. (2010). Feature selection with the Boruta package. J. Stat. Softw.36 (11), 1–13. 10.18637/jss.v036.i11
45
LeCunY.BengioY.HintonG. (2015). Deep learning. nature521 (7553), 436–444. 10.1038/nature14539
46
LiQ.HanZ.WuX.-M. (2018). “Deeper insights into graph convolutional networks for semi-supervised learning,” in Proceedings of the AAAI conference on artificial intelligence.
47
LiY.SongY.ChenX.ShiZ.KaskaoutisD. G.GholamiH.et al (2023). Late Pleistocene dynamics of dust emissions related to westerlies revealed by quantifying loess provenance changes in North Tian Shan, Central Asia. Catena227, 107101. 10.1016/j.catena.2023.107101
48
LiY.SongY.KaskaoutisD. G.ZhangX.ChenX.ShukurovN.et al (2022). Atmospheric dust dynamics over central Asia: A perspective view from loess deposits. Godwana Res.109, 150–165. 10.1016/j.gr.2022.04.019
49
LiuX.SongH.LeiT.LiuP.XuC.WangD.et al (2021). Effects of natural and anthropogenic factors and their interactions on dust events in Northern China. CATENA196, 104919. 10.1016/j.catena.2020.104919
50
LundbergS. M.LeeS.-I. (2017). “A unified approach to interpreting model predictions,” in Advances in neural information processing systems, 30.
51
LuoQ.ZhenL.XiaoY.WangH. (2020). The effects of different types of vegetation restoration on wind erosion prevention: A case study in yanchi. Environ. Res. Lett.15 (11), 115001. 10.1088/1748-9326/abbaff
52
MarjovviA.SoleimaniM.MirghaffariN.KarimzadehH.YuanY.FangL. (2022). Monitoring, source identification and environmental risk of potentially toxic elements of dust in Isfahan Province, Central Iran. Bull. Environ. Contam. Toxicol.108 (5), 901–908. 10.1007/s00128-021-03446-7
53
MiddletonN. J. (2017). Desert dust hazards: A global review. Aeolian Res.24, 53–63. 10.1016/j.aeolia.2016.12.001
54
MiddletonN. (2019). Variability and trends in dust storm frequency on decadal timescales: Climatic drivers and human impacts. Geosci. Switz.9 (6), 261. 10.3390/geosciences9060261
55
MohammadifarA.GholamiH.CominoJ. R.CollinsA. L. (2021a). Assessment of the interpretability of data mining for the spatial modelling of water erosion using game theory. Catena200, 105178. 10.1016/j.catena.2021.105178
56
MohammadifarA.GholamiH.GolzariS. (2022). Assessment of the uncertainty and interpretability of deep learning models for mapping soil salinity using DeepQuantreg and game theory. Sci. Rep.12 (1), 15167. 10.1038/s41598-022-19357-4
57
MohammadifarA.GholamiH.GolzariS.CollinsA. L. (2021b). Spatial modelling of soil salinity: Deep or shallow learning models?Environ. Sci. Pollut. Res.28, 39432–39450. 10.1007/s11356-021-13503-7
58
MontiF.BronsteinM.BressonX. (2017). Geometric matrix completion with recurrent multi-graph neural networks. Adv. neural Inf. Process. Syst.30.
59
MuthukrishnanR.RohiniR. (2016). “Lasso: A feature selection technique in predictive modeling for machine learning,” in 2016 IEEE international conference on advances in computer applications (ICACA). IEEE.
60
OldemanL. R. (1994). “The global extent of land degradation,” in Land resilience and sustainable land use. Editors GreenlandD.SzabolcsI. (Wallingford: CAB International).
61
PadarianJ.McBratneyA. B.MinasnyB. (2020). Game theory interpretation of digital soil mapping convolutional neural networks. Soil6 (2), 389–397. 10.5194/soil-6-389-2020
62
PapiR.AttarchiS.AttarchiS.Darvishi BolooraniA.SamaniN. N. (2022). Characterization of hydrologic sand and dust storm sources in the Middle East. Sustainability14, 15352. 10.3390/su142215352
63
ParajuliS. P.StenchikovG. L.UkhovA.KimH. (2019). Dust emission modeling using a new high-resolution dust source function in WRF-chem with implications for air quality. J. Geophys. Res. D. Atmos. JGR124, 10109–10133. 10.1029/2019jd030248
64
ParajuliS. P.YangZ. L.KocurekG. (2014). Mapping erodibility in dust source regions based on geomorphology, meteorology, and remote sensing. J. Geophys. Res. Earth Surf.119 (9), 1977–1994. 10.1002/2014jf003095
65
PourhashemiS.AsadiM. A. Z.BoroughaniM.AzadiH. (2022). Mapping of dust source susceptibility by remote sensing and machine learning techniques (case study: Iran-Iraq border). Environ. Sci. Poll. Res.30, 27965–27979. 10.1007/s11356-022-23982-x
66
RashkiA.ArjmandM.KaskaoutisD. (2017). Assessment of dust activity and dust-plume pathways over Jazmurian Basin, southeast Iran. Aeolian Res.24, 145–160. 10.1016/j.aeolia.2017.01.002
67
RashkiA.MiddletonN.GoudieA. (2021). Dust storms in Iran – distribution, causes, frequencies and impacts. Aeolian Res.48, 100655. 10.1016/j.aeolia.2020.100655
68
ReissS.BogenbergerK. (2016). Validation of a relocation strategy for Munich's bike sharing system. Transp. Res. Procedia19, 341–349. 10.1016/j.trpro.2016.12.093
69
SahaA.PalS. C.ArabameriA.ChowdhuriI.RezaieF.ChakraborttyR.et al (2021). Optimization modelling to establish false measures implemented with ex-situ plant species to control gully erosion in a monsoon-dominated region with novel in-situ measurements. J. Environ. Manag.287, 112284. 10.1016/j.jenvman.2021.112284
70
SahaA.PalS. C.ChowdhuriI.IslamA. R. M. T.ChakraborttyR.RoyP. (2022). Application of neural network model-based framework approach to identify gully erosion potential hotspot zones in sub-tropical environment. Geocarto Int.37, 14758–14784. 10.1080/10106049.2022.2091042
71
ScarselliF.GoriM.TsoiA. C.HagenbuchnerM.MonfardiniG. (2009). The graph neural network model. IEEE Trans. neural Netw.20 (1), 61–80. 10.1109/tnn.2008.2005605
72
ShaoY. (2008). Physics and modelling of wind erosion. Springer Science and Business Media.
73
SharifikiaM. (2013). Environmental challenges and drought hazard assessment of Hamoun Desert Lake in Sistan region, Iran, based on the time series of satellite imagery. Nat. hazards65, 201–217.
74
SissakianV.Al-AnsariN.KnutssonS. (2013). Sand and dust storm events in Iraq. J. Nat. Sci.5 (10), 1084–1094. 10.4236/ns.2013.510133
75
SolokiH. R.Noori NahadF. (2019). Effect of soil characteristics and wind erosion on hydraulic structures of Sistan River. J. Nat. Environ. Hazards8 (19), 229–242.
76
SperdutiA.StaritaA. (1997). Supervised neural networks for the classification of structures. IEEE Trans. Neural Netw.8 (3), 714–735. 10.1109/72.572108
77
TianM.GaoJ.ZhangL.ZhangH.FengC.JiaX. (2021). Effects of dust emissions from wind erosion of soil on ambient air quality. Atmos. Pollut. Res.12 (7), 101108. 10.1016/j.apr.2021.101108
78
TibshiraniR. (1996). Regression shrinkage and selection via the lasso. J. R. Stat. Soc. Ser. B Methodol.58 (1), 267–288. 10.1111/j.2517-6161.1996.tb02080.x
79
Van DaoD.JaafariA.BayatM.Mafi-GholamiD.QiC.MoayediH.et al (2020). A spatially explicit deep learning neural network model for the prediction of landslide susceptibility. Catena188, 104451. 10.1016/j.catena.2019.104451
80
VeličkovićP.CucurullG.CasanovaA.RomeroA.LioP.BengioY. (2017). Graph attention networks. arXiv preprint arXiv:1710.10903.
81
WebbN. P.HerrickJ. E.Van ZeeJ. W.CourtrightE. M.HugenholtzC. H.ZobeckT. M.et al (2016). The National Wind Erosion Research Network: Building a standardized long-term data resource for aeolian research, modeling and land management. Aeolian Res.22, 23–36. 10.1016/j.aeolia.2016.05.005
82
WebbN. P.StrongC. L. (2011). Soil erodibility dynamics and its representation for wind erosion and dust emission models. Aeol. Res.3, 165–179. 10.1016/j.aeolia.2011.03.002
83
WuZ.PanS.ChenF.LongG.ZhangC.PhilipS. Y. (2020). A comprehensive survey on graph neural networks. IEEE Trans. neural Netw. Learn. Syst.32 (1), 4–24. 10.1109/tnnls.2020.2978386
84
XiaoY.XieG.ZhenL.LuC.XuJ. (2017). Identifying the areas benefitting from the prevention of wind erosion by the key ecological function area for the protection of desertification in Hunshandake, China. Sustainability9 (10), 1820. 10.3390/su9101820
85
XuJ.XiaoY.XieG.WangY.JiangY. (2019). Computing payments for wind erosion prevention service incorporating ecosystem services flow and regional disparity in Yanchi County. Sci. Total Environ.674, 563–579. 10.1016/j.scitotenv.2019.03.361
86
XuJ.XiaoY.XieG.ZhenL.WangY.JiangY. (2018). The spatio-temporal disparities of areas benefitting from the wind erosion prevention service. Int. J. Environ. Res. Public Health15 (7), 1510. 10.3390/ijerph15071510
87
YesilnacarE. K. (2005). The application of computational intelligence to landslide susceptibility mapping in Turkey. University of Melbourne.
88
ZhangQ.YangL. T.ChenZ.LiP. (2018). A survey on deep learning for big data. Inf. Fusion42, 146–157. 10.1016/j.inffus.2017.10.006
89
ZhangS.TongH.XuJ.MaciejewskiR. (2019). Graph convolutional networks: A comprehensive review. Comput. Soc. Netw.6 (1), 11–23. 10.1186/s40649-019-0069-y
90
ZhangS.ZhouD.YildirimM. Y.AlcornS.HeJ.DavulcuH.et al (2017). “Hidden: Hierarchical dense subgraph detection with application to financial fraud detection,” in Proceedings of the 2017 SIAM international conference on data mining (SIAM).
91
ZhangY.ThomasT.BrusselM.Van MaarseveenM. (2016). Expanding bicycle-sharing systems: Lessons learnt from an analysis of usage. PLoS one11 (12), e0168604. 10.1371/journal.pone.0168604
Summary
Keywords
deep learning, erodible susceptibility, feature selection, game theory, potential impact
Citation
Gholami H, Mohammadifar A, Fitzsimmons KE, Li Y and Kaskaoutis DG (2023) Modeling land susceptibility to wind erosion hazards using LASSO regression and graph convolutional networks. Front. Environ. Sci. 11:1187658. doi: 10.3389/fenvs.2023.1187658
Received
16 March 2023
Accepted
11 April 2023
Published
09 May 2023
Volume
11 - 2023
Edited by
Giulio Castelli, University of Florence, Italy
Reviewed by
Juan Cruz Colazo, Instituto Nacional de Tecnología Agropecuaria, Argentina
Jianhua Xiao, Northwest Institute of Eco-Environment and Resources (CAS), China
Updates
Copyright
© 2023 Gholami, Mohammadifar, Fitzsimmons, Li and Kaskaoutis.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hamid Gholami, hgholami@hormozgan.ac.ir
† Present address: Dimitris G. Kaskaoutis, Department of Chemical Engineering, University of Western Macedonia, Kozani, Greece
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.