Modelling the effects of topographic heterogeneity on distribution of Nitraria tangutorum Bobr. species in deserts using LiDAR-data

Microclimate ecology is attracting renewed attention because of its fundamental importance in understanding how organisms respond to climate change. Many hot issues can be investigated in desert ecosystems, including the relationship between species distribution and environmental gradients (e.g., elevation, slope, topographic convergence index, and solar insolation). Species Distribution Models (SDMs) can be used to understand these relationships. We used data acquired from the important desert plant Nitraria tangutorum Bobr. communities and desert topographic factors extracted from LiDAR (Light Detection and Ranging) data of one square kilometer in the inner Mongolia region of China to develop SDMs. We evaluated the performance of SDMs developed with a variety of both the parametric and nonparametric algorithms (Bioclimatic Modelling (BIOCLIM), Domain, Mahalanobi, Generalized Linear Model, Generalized Additive Model, Random Forest (RF), and Support Vector Machine). The area under the receiver operating characteristic curve was used to evaluate these algorithms. The SDMs developed with RF showed the best performance based on the area under curve (0.7733). We also produced the Nitraria tangutorum Bobr. distribution maps with the best SDM and suitable habitat area of the Domain model. Based on the suitability map, we conclude that Nitraria tangutorum Bobr. is more suited to southern part with 0–20 degree slopes at an elevation of approximately 1010 m. This is the first attempt of modelling the effects of topographic heterogeneity on the desert species distribution on a small scale. The presented SDMs can have important applications for predicting species distribution and will be useful for preparing conservation and management strategies for desert ecosystems on a small scale.

The Ulan Buh Desert is the 8th largest desert in China, and is located in the mid-latitudinal zone.It is characterized by a typical continental climate with little precipitation.The desert is located in a transitional zone with a temperate climate ranging from semiarid to arid conditions 1 .It is one of the driest regions at similar latitudes in the world and has a very sensitive ecological environment, which creates a fragile desert ecotone.The desert ecotone is the key ecoregion and is suitable for researching ecosystem degradation and recovery mechanisms.Desert ecotones may have high ecological significance, which quantitatively reveals the interactions between plants and environmental factors.The spatial distribution pattern of the desert plants is therefore the result of the interactions between the plant communities and the harsh environmental factors 2 .
Spatial distributions are fundamental to species ecology 3 and help identify threats to the conservation of plant species, manage the impacts of biological invasion 4 , and plan conservation strategies.Climate, edaphic, and topographic factors are the main factors affecting species distributions 5 .These factors are affected by a variety of biophysical processes including competitive interactions between plant species and their dispersal history 6,7 .SDMs, like habitat suitability or ecological models, are developed with the observed species occurrence as a response variable and environmental factors as predictor variables 8 .Many modelling techniques are available to developed SDMs [9][10][11] .Although previous SDMs studies have compared different modelling techniques for their precision [12][13][14][15][16][17] , most of them describe species distributions over a large scale, and therefore have low predictive performance in small-localized condition 6,[18][19][20] .Most SDMs are calibrated only with climate data at a low spatial resolution, such as the climate data acquired from WorldClim(https:// www.world clim), with a 1 km 2 grid resolution.However, some studies have indicated that the predictive performance of SDMs can be substantially improved by incorporating other environmental factors in addition to climate factors 21 or by considering finer-scale processes, such as topo-climatic processes 22 .Topographic heterogeneity is the most important factor affecting the spatial distribution of desert plants 23 .In the mountainous regions, topography may precisely describe local temperature, light, and humidity, Topography significantly affects light and water availability and soil development.Furthermore, topography affects species composition, structure, appearance, and dynamics of the plant communities [24][25][26] .
Building on the previous work done in this field, the current study intends to answer a few important questions including 'What distribution patterns are followed by the desert plant species?' and 'Which topography factors play a major role in affecting the species distribution?' Intense solar radiation, high evapotranspiration, extreme temperatures, degraded soil and vegetation cover, low atmospheric moisture, and harsh topographic features are important factors in the desert environment.The desert is extremely arid with < 60-100 mm mean annual precipitation, and therefore, its moisture regime is the most critical factor for regulating biological processes and species distribution.The rainfall spatially varies at different scales, not only at the regional scale ranging from 0.1 to 10 km, but also at a local scale of few meters.Rainfall is the main source of moisture for desert plants, which triggers species' survival and growth 27 .Many studies have explored the appropriate methods for understanding the relationships between vegetation dynamics, soil moisture, and temperature in the desert ecosystem [28][29][30] .However, under similar climate and desert environmental conditions, the influence of topography on species distribution has rarely been investigated in a small scale.
Additionally, several studies have been conducted on the topographic factors affecting animal communities in the desert.In recent decades, some studies have been conducted to model the environmental variables influencing kit fox distribution.For example, Maxent fixed-effects and generalized linear mixed-effects models have been developed to describe kit fox distribution 30 .These modelling studies show the large effects of elevation on kit fox survival and growth, followed by that of the slope and canopy height of the vegetation at the landscape scale.The Maxent algorithm was used to predict the habitat in the Mojave desert and parts of the Sonoran desert considering the geographical and topographical factors at the spatial grid of one square kilometer 31 .The field-validated site-and landscape-level SDMs were developed to identify the potentially rare and endemic plant habitats in the Great Basin of Western North America 32 , in which the potential habitat combined with the elevation, slope, aspect, rock type, and geologic processes are mapped with the resolution of 100 square meters.Similar to climatic factors, non-climate predictors (e.g., topography and habitat) are also important at a finer scale 33 .Thus, the analysis of these factors on a finer scale could be the most appropriate for developing SDMs and increasing their predictive accuracy 34 .As mentioned above, topography is one of the most important predictors of SDMs.It controls the habitat structure, biota, and all growing conditions locally 35 , especially in the desert, where topographic features are more unfavorable than other land features.Thus, a high-resolution topographic map is necessary to develop SDMs.
For dominant species, such as Nitraria tangutorum Bobr. in the desert, a small geographic range and narrow habitat may lead to acute habitat limitation, which is strongly correlated with other environmental factors.The microclimate, which influences ecosystem dynamics and processes, is often disregarded in ecology and evolution 36 .It is meaningful to investigate the geographical variations of Nitraria tangutorum Bobr.species in the desert on a small scale, such as within the scale of 100 m, where topographic factors significantly affect species distribution.Existing studies have shown that geographic and topographic predictors can substantially improve the fitting performance of SDMs 37 .
High-resolution terrain data are rarely available.Remote sensing technology, which offers alternative solutions to various data sets, can supply accurate high-resolution terrain data using the light detection and ranging (LiDAR) technology.Remote sensing data is important for determining species distribution 38,39 .In recent years, LiDAR datasets have become useful for modelling biodiversity and habitat analyses.For example, LiDAR data could improve the prediction of composition and changes in the plant communities 40,41 , improve the descriptions of animal behaviors and their spatial distributions 42 , and enhance the species richness prediction 17 .
In this study, we developed SDMs for revealing the relationships between desert plants and topographic factors affecting species distributions.Topographic variables used in the SDMs were derived from LiDAR data (resolution: 0.15 m).The SDMs could be used to predict the precise distribution of species covering an area of one square kilometer with a resolution of 10 m.Different modelling algorithms proposed in this study were evaluated based on the dataset obtained from 67 Nitraria tanguotrum species in the Ulan Buh Desert, an arid land in northern China.The objectives of this study were to: (1) compare various modelling algorithms including parametric and non-parametric modeling based on statistical indexes, and determine the best to develop SDMs; (2) evaluate the contribution of topographic factors to the precise prediction of distribution of the plant species of interest in the desert; (3) determine the effects of topographic heterogeneity on the distribution of Nitraria tangutorum Bobr; and (4) present the habitat suitability map for Nitraria tangutorum Bobr.and show the corresponding management implications of the presented SDMs.
Habitat suitability and prediction maps for Nitraria tangutorum Bobr.distribution.We produced the maps (Sup.Figures 1-7) using R 3.4.3and ARCGIS 10.5, which show the habitat suitability for Nitraria tangutorum Bobr.The potentially suitable habitat distribution maps of Nitraria tangutorum Bobr. in the study area were obtained from different algorithms.Furthermore, the corresponding species' presence/absence with each algorithm was mapped.The Nitraria tangutorum Bobr.species distributed mainly in the eastern part of the study area.The spatial distributions of Nitraria tangutorum Bobr.simulated by the different models were significantly different.Concerning the most suitable area for Nitraria tangutorum Bobr, the most suitable habitat area simulated by the Domain model was the largest, and the most suitable habitat area simulated by BIOCLIM model was the smallest.The suitable habitat areas simulated by GLM and GAM models were identical; the suitable habitat simulated by GLM and GAM was the largest, while that simulated by BIOCLIM was the smallest.
The driven factors of topography on the Nitraria tangutorum Bobr.distribution.The responses of the driven factor-topography with different models were compared for Nitraria tangutorum Bobr.Different models delineated the topography characteristics of the potential distribution areas.With BIOCLIM, the existence probability followed a trend of increase at first and then decreasing trend, and aspect of the south part was on the top.For DEM, the existence probability trend of the species was steeper, reaching a peak at an altitude of approximately 1010 m.For the profile factor, the species Nitraria tangutorum Bobr.appeared more suitable at the predicted value of zero, and was also suitable for slope of 0°-20°.With the Domain model, the existence probability trend of the species was similar to that of the BIOCLIM model, but it seemed steeper.In the Malanhanobis model, at the aspect of the north by east 10°, the species Nitraria tangutorum Bobr.exhibited minimum fitness, and at the altitude range of 1007-1013 m, exhibited maximum fitness.Regarding the topography factor of the profile, there were two crests at the predicted value near zero.With the GLM and GAM models, the existence probability trend of the species was nearly the same as the linear trend line, but in the RF, the trend was nonlinear.Regarding the factor of aspect, predicted value fluctuated between 0.4 and 0.8, but for DEM factor, www.nature.com/scientificreports/predicted value was the highest at an altitude of approximately 1010 m.With the factors of profile and slope, at the degree of zero, predicted value reached the top, at a slope of 20°-60°, exhibiting a rising tendency.Analysis of the maps showed the distribution of the species Nitraria tangutorum Bobr.was more suitable in the south, with an altitude of approximately 1010 m, from a slope between 0° and 20°, and general trend was the same with different models.The GLM or GAM was found inadequate for the distribution modelling of Nitraria tangutorum Bobr.

Discussion
Importance of LiDAR data and its ecological implications with selected predictor variables.How does climate change affect organisms?In the past, our limited ability to map and monitor the microclimatic variations using modelling through the incorporation of the microclimate.The microclimate sensors set in the networks provided point-based approaches and weather stations provided macroclimate data.However, means of interpolating and downscaling data were lacking 43 .LiDAR is valuable for mapping microclimates, as it provides spatially continuous and sub-meter-scale information about ground topography 44 .Aerial photography provides an approach for assessing topography.However, these are less accurate than LiDAR at deriving terrain elevation 45 .An important advantage of UAVs is that they are flexible, enabling the collection of time series of aerial imagery over a period at a reasonable cost.In our study, LiDAR data were used to obtain DEM for SDMs, thereby supporting more effective conservation monitoring, management, and policy decision for a sustainable future 46 .
Our results showed that Nitraria tangutorum Bobr.species are specialists in the range of their habitat resources.In other words, Nitraria tangutorum Bobr.has a relatively narrow ecological niche.
It is limited to making the distribution by considering only topographic variables, the role of habitat heterogeneity in the desert ecosystem may be underestimated.If other environmental factors, such as light, soil, and water are considered, niche theory may be able to better explain the existence mechanism of species.Therefore, future studies need to consider more bio-variables to simulate the distribution of the species, and a permanent plot should be designed to contain the entire spatial gradient of the microclimate conditions in the study system.Long-term microclimate data series are required to complement the dynamics of the microclimate around the species.Gathering georeferenced microclimate data from different habitat types across the globe would  www.nature.com/scientificreports/significantly promote progress in microclimate ecology 47 .UAVs equipped with multispectral or hyperspectral sensors and laser scanner systems are very promising, providing high-resolution data in the context of microhabitats, vegetation structures, and topography 48 .
Suggestion for improvement of SDMs.In our study, not only AUC of the models but also the responses of the variables with the models were different.The parametric regression-based models (GLM and GAM) had lower mean performance, whereas nonparametric or machine learning-based models (RF and SVM) were more accurate.GLMs are widely used in SDM studies 49 , but have poor nonlinear response performance.At the same time, the results confirmed the effectiveness and robustness of the machine-learning techniques.Similarly, Marmion et al. found a single algorithm with lower accuracy than RF in predicting the plant distributions 50 .The parameters of the model affect its accuracy; for example, the performance of RF depends on two important parameters, ntree and mty.As for the ntree, the optimal value using only the default may not be the best; therefore, future research should test and adjust the model parameters to the optimum state.In this study, various parametric and nonparametric algorithms were evaluated for SDMs.In the next step, semi-parametric models 51 shall be used for better results.
A stable and reliable mathematical algorithm is required to develop more accurate SDMs.In previous studies of community-level predictions of the species distributions, the choice of how species and site-level occurrence probabilities are combined into species assemblage predictions is more important than the choice of the model type used 52 .Following research shall use ensemble models to predict species distribution.As many researchers have expressed, no evidence has proven that one model is steadily accurate, but it is proven to be better than other models.
The number of environmental factors that accurately represent the habitat characteristics of a species was sufficient SDM has some limitations, for example, over-prediction of species richness per site and sensitivity to methodological biases.In addition, sample size, sample prevalence, sample design, model techniques, imperfect detection of species, or the choice of environmental variables could affect the uncertainty of the predictions 53 .
The relationships between species occurrence and sets of predictor variables explored by the models produce two kinds of useful outputs: estimates of the probability that species might occur at a given unrecorded location 50 , and estimates of an area's suitability for species 16 .Climate determines the distribution of desert species 54 .Previous research has shown that the roots of desert plant species develop and depend on the underwater 55 .As the study area was only one square kilometer, the variable of climate change could not be obvious, and species distribution could have been disturbed by the level of underground water.However, further studies are necessary in this regard.It is well known that biotic interactions play an important role in creating accurate SDMs for many species, particularly at small geographic scales, such as host requirements and competition 6,56 .In this study, the model did not include the interaction between species and the influence of predators or human activities.More efforts should be dedicated to include fine-scale environmental measures, such as micro-temperature, underwater, moisture, and rainfall, for fine-scale species mapping and management.
Management implication of SDMs.Species conservation and management are key to maintaining regional ecological balance, especially in ecologically fragile areas such as the Ulan Buh Desert.This study could distinguish the habitat requirements of Nitraria tangutorum Bobr., especially for the topography, to provide a reference for protection.Nitraria tangutorum Bobr.was found to be distributed mainly in the eastern part of the study area.This finding is consistent with the results of the actual survey.The suitable habitat area is larger than the real area, which indicates that the area of Nitraria tangutorum Bobr. is more widely distributed under current environmental conditions.Some of the suitable habitat characteristics shown in our study were consistent with those of other studies 57 .For example, the population size of Nitraria tangutorum Bobr.increases with sand depth, making it suitable for dunes.To expand the area of Nitraria tangutorum Bobr., the land should be dune, and will be a good choice at an altitude of approximately 1010 m, better to face the aspect of the south part, within the slope of 20°.Because of this characteristic, Nitraria tangutorum Bobr.can prevent quicksand from moving forward and land desertification.Nitraria tangutorum Bobr. is an excellent species for sand fixation.

Study workflow.
The study process shown in the flowchart (Fig. 2) consists of four parts: (1) preparing data, including species data and factors of topography; (2) fitting different models; (3) evaluating models by Area Under the Curve (AUC) and comparing the driving factors of different models; and (4) estimating areas suitable for habitat using the developed SDMs.

Study area and data analysis. Study area. The Ulan Buh Desert is located in the northern part of
Inner Mongolia Province, 106°38′42″E-106°57′00″ E, 40°17′24″ N-40°28′36″ N (Fig. 3), covering an area of 14,905.13km 2 .The elevation ranges from 1009 to 1016 m, and the topographical profile on the surface varies from 1211.58 to 1339.63 m 58 .The climate is characterized by semi-arid to arid conditions, with a mean annual precipitation of 147.4 mm, mean annual evaporation of 2458.4 mm, mean annual relative humidity of 47.3%, and mean annual temperature of 6.8 °C in the desert 59 .The eastern edge of the desert is important for dividing the desert and grassland in Central Asia.Our study was based on control plots, such as desertification processes and control treatments 60 .(1) Species data collection The ground surveys were conducted between August and September 2019.100 sample plots of Nitraria tangutorum Bobr.plant community was established in one square kilometer at the eastern edge of the study area.Each 100 × 100 m sample plot was divided into 100 small quadrats with an area of 10 × 10 m (Fig. 3).The center points GPS location is southeast corner of the sample F5, and the accuracy of the point is within a half meter.The grids that were created in the same direction as we sampled in field, all are facing north.335 occurrence records for Nitraria tangutorum Bobr.were collected in the area (Fig. 5).The formal identification of the plant material used in my study has complied with the Chinese Virtual Herbarium (CVH): (https:// www.cvh.ac.cn/ spms/ detail.php? id= f76ac eff), and a specimen of this material has been deposited in a publicly available herbarium.And we were employed by the Chinese Academy of Forestry, which is also the State Forestry Administration Dengkou Desert ecosystem positioning observation station.To ensure that we have permission www.nature.com/scientificreports/ to collect Nitraria tangutorum Bobr.We confirm that all methods were performed following current guidelines and regulations in China, Huoyan zhou, Xiao Zhou, and Xiaodi Zhao undertook the formal plant identification, and field investigation were carried out.

Species data processing. Nitraria tangutorum
(2) Processing species data The data were cleaned for use in our modelling.The longitudes and latitudes of the sample plot center were recorded as coordinates, and if the record was out of the boundary of the study area, it was checked with the UAV-RGB (Unmanned Aerial Vehicle-Red Green Blue) data and overlaid on the map.If the canopy was outside the plot boundary, the root inside was preserved.
The coordinates were cross-checked by a visual or 'overlay' function using ArcGIS 10.5.This study applied the coordinate function from the 'SP' (spatial point) package to create a Spatial Point Data Frame, and then the over function from 'SP' to do a point-in-polygon query with the study boundary.Moreover, the occurrence of each grid would be only once, without duplicate records; the function ' duplicate' could be used to remove duplicates in the software R3.4.3.
The species data included the presence and absence data, the resolution of the grid was 10 m, and each grid retained only presence or absence data, reducing the bias of the model predictions by space sampling 61 .

Environment data processing. LiDAR data
The UAV technique with the SPAN_IGM_STIM300 IMU (Fig. 4a) was used to gather airborne LiDAR data during the study conducted on October 24, 2019.The LiDAR parameters were a pulse rate of 125 kHz, a scan angle (FOV(field of view)) of ± 22.5°, and a laser beam divergence (IFOV) of 0.5 mrad.An average point density of 41.62 points per m2 and a footprint diameter of 25 cm was obtained.www.nature.com/scientificreports/

Topographic variables
The digital elevation model (DEM) was derived from the UAV LiDAR point cloud (Fig. 4b) with a resolution of 15 cm (Fig. 5a) and then resampled to 10 m (Fig. 5b) with image analysis tools.Other topographic variables, including slope, aspect, and profile (Fig. 6), were obtained using the spatial analysis model in ArcGIS 10.5.Predictor variables were organized as raster (grid) files for species distribution modelling.After generating the set of predictor variables (rasters) and the occurrence points, our next step was to extract the values of the predictors at the locations of the points using the raster tools of 'extract' in ArcGIS 10.5.To select the significant predictor variables, the variance inflation factor (VIF) was calculated.A VIF of 10 was acceptable as it showed no significant collinearity among the predictor variables 62 .All variables were resampled using ArcGIS 10.5 to unify the space dimension of 10 m, which was the same value as the resolution of the species samples.
Forms for SDMs.To predict species distribution, a number of statistical and machine-learning approaches have been used 63,64 .In conjunction with geographic information systems and remote sensing, the methods used in species distribution modelling can be classified as profile, regression, and machine learning methods.The profile method only considers the presence (occurrence) data, whereas the regression and machine learning methods use both the presence and absence data.The profile method includes BIOCLIM 65 , Domain 66 , and Mahalanobis 67 .In addition to these methods, we also analyzed and compared the predictive performance of SDMs developed using regression methods, such as the GLM 68 and GAM 69,70 , and some machine learning methods, such as RF 71 and SVM.www.nature.com/scientificreports/Sampling bias may have been present in the occurrence data used 72 .Random sampling was performed to record the presence or absence data with the random points function in the 'dismo' package.An 'extent' tool was used to further restrict the area from where random locations were drawn 8,73 .The 'dismo' package containing the k-fold function facilitates the data partitioning and creates a vector that assigns each row in the data matrix to a group (between 1 and k with k = 5).
We used seven modelling algorithms, including both nonparametric and parametric algorithms, to develop the SDMs, seven modelling algorithms are briefly described below.

Bioclimatic modelling (BIOCLIM).
BIOCLIM is a classic climate envelope model 65 .Although it generally does not perform as well as some other modelling methods 74 , particularly in the context of climate change 72 , it is still being used, because the algorithm is easy to understand. 66has been extensively used for species distribution modelling.The Domain algorithm computes Gower's metric 75 between the environmental variables at any location and those at any of the known locations of occurrence, which were defined as training sites in this study.

Domain. The Domain algorithm
A suitable means of quantifying the similarity between the two sites is provided by Gower's metric.In Euclidean Q-dimensional space, the distance d between two points A and B is defined as To equalize the contribution from each climatic attribute, Gower's metric uses range standardization.In this application, the standardization method is preferred over variance standardization because it is less susceptible to bias arising from dense clusters of sample points.The complementary similarity measure ( R AB ) is defined as R AB is constrained between 0 and 1 for points within the ranges used in Eq. 1 but may yield negative values for points outside this range.S A is defined as the maximum similarity between candidate point A and the com- bination of known record site Tm: S A is evaluated for all grid points in a target area, where m denotes the number of known records, and a matrix of continuously varying similarity values is generated, which can be displayed as a grayscale, thematic, or contour map.As with all models discussed here, the values generated are not probability estimates but the degrees of classification confidence 66 .

Mahalanobis.
The Mahal function implements SDMs based on Mahalanobis distance 67 .The Mahalanobis distance considers the correlations of the variables in the dataset, and it is not dependent on the scale of measurements.
(1) where α µv is the dispersion, (α) µ is the mean value, P is the independent variable, and D 2 is the sample value of the 2 statistic.Both the dispersion and mean values are subject to sampling fluctuations.

Generalized linear models (GLM).
In GLM, X k (k = 1, • • •, r) is the vector of k predictor variables, which is related to the expected value µ = E(Y ) of the response variable Y (denote the response variable) through a link function g(•) (Eq.5).The use of GLM in species distribution modelling is provided in 68 .
where α is a constant called the intercept, and β is the vector of k regression coefficients (one for each predictor).(Eq.( 6)).
It is assumed that the topographic factor is normally distributed and is defined as X, and the species distribution probability as Y; Therefore in the response distribution of GLM, which uses the "family = gaussian".

Generalized additive models (GAM).
GAM 69,70 are extensions of GLMs.In GAMs, the linear predictor is the sum of the smoothing functions.This makes GAMs very flexible and they can fit complex functions.This  www.nature.com/scientificreports/also makes them very similar to machine-learning methods.The GAMs were implemented using the 'mgcv' package R 3.4.3.
where µ = E(Y X 1 , X 2 , • • •X p ) , n is for the linear prediction, and s i (•) is the nonparametric smooth function.This model does not require any assumptions.It consists of a random component, an additive component, and a link function connecting the two parts.The distribution of the response variables(Y) belonging to the family of exponential distributions can be binomial.

Random forest (RF).
The RF method is an extension of the regression and classification trees (CART) 71 .
This method can be implemented using the "random forest" function in the package with the same package name in R 3.4.3.The function RF can take a formula, or in two separate arguments, a data frame with the predictor variables and a vector with the response.If the response variable is a factor (categorical) type, RF performs classification; otherwise, it performs regression.In species distribution modelling, because RF showed remarkable modelling results, it is considered a viable option.
(1) Suppose there is a dataset ) ( X im is the vector of m predictor variables, y i are the response variables) with feature number n, sampling with replacement can generate a sampling space (m * n) m * n .(2) Build a basic learner (a decision tree): ) generates a decision tree for each sample (where K < < M) and records the results of each decision tree h j (x).
(3) Training times to T H(x) = max T t=1 φ(h(x) = y) , where φ(x) is an algorithm (absolute majority voting, majority voting, weighted voting method, etc.).RF requires two parameters: (1) mtry, the number of predictor variables performing data partitioning at each node, and (2) ntree, the total number of trees to be grown in the model run.In this study, ntree (number of trees to grow) was set to 500 and mtry to 10 after some initial tuning experiments. 76applies a simple linear method to data; however, in a highdimensional feature space, nonlinearly exists in the input space.However, in practice, this does not involve any computations in a high-dimensional space.This simplicity, combined with state-of-the-art performance on many learning problems (classification, regression, and novelty detection), has contributed to the popularity of SVM 77 .They were first used in species distribution modelling by Guo et al. 78 .There are several implementations of SVM in R 3.4.3.The most useful implementations in our context are the function 'ksvm' in the package 'kernlab' and the 'svm' function in the package ' e1071' .The 'ksvm' includes many different SVM formulations and kernels, provides useful options, and features a method for plotting, but it lacks a proper model selection tool.The 'svm' function in package ' e1071' includes a model selection tool: the 'tune' function 77 .

Support vector machines (SVM). SVM
Given the training examples, x i (i = 1, 2, • • •m) indicates the samples to be classified, y ∈ {+1, −1} indicates the marked label value for the samples, + 1 denotes the point above the separating hyperplane, and -1 denotes the point under the separating hyperplane, where K(x, x i ) is the Gaussian kernel function, K(x, x i ) is nonlinearity, and * i is the optimal solution to the primal problem obtained by solving the dual problem.Note that if a sample x i is not a support vector, to maximize the Lagrangian function, there must be a λ* = 0 non-support vector cor- responding to * equal to zero, which theoretically shows that the decision boundary of the SVM is only related to the support vector.
Evaluation and comparison of SDMs.Accuracy assessment.Different measures may be used to evaluate the prediction ability of each model [79][80][81] .Many measures for evaluating models based on presence-absence or presence-only data are threshold dependent.This implies that a threshold must first be set.However, many statistics are threshold independent, such as the correlation coefficients and the area under the receiver operator curve (AUROC, generally abbreviated as AUC).The AUC value is the area of ROC and the X-axis under the ROC curve; the larger the area, the higher the relationship, and the greater the creditability.When the AUC value is greater than 0.7, the result is fair, and when the AUC is higher than 0.9, the result is excellent 82 .We applied this measure to evaluate the SDMs.
In this framework, an algorithm is supposed to predict "positive" or "negative".However, some concepts are confusing, as summarized in Table 2.

Suitable habitat index (SHI).
The probability of species occurrence was extracted using ArcGIS.The lowest probability value was selected to distinguish the suitable zone from the non-suitable zone 83 .The habitats of the entire study area were divided into three categories using the frequency statistical method 84 : unsuitable habitat (SHI ≤ 0.3), low suitable habitat (0.3 < SHI ≤ 0.5), and suitable habitat (SHI > 0.5).

Figure 1 .
Figure 1.The boxplots of different models (BIOCLIM: classic climate envelope model; Domain computes the Gower distance between environmental variables and locations of occurrence; Mahalanobis model based on the mahalanobis distance; GLM: Generalized Linear Model; GAM: Generalized Additive Model; RF: Random Forest; SVM: Support Vector Machine).

Figure 2 .
Figure 2. Workflow of the study.ARCGIS 10.5 is the software of geographic information system, R3.4.3 is the data processing platform, AUC is the area under the receiver operator curve, which evaluates the fitted models.Species distribution models include Bioclimatic Modelling (BIOCLIM), Domain, Mahalanobis, Generalized Linear Model (GLM), Generalized Additive Model (GAM), Random Forest (RF), Support Vector Machine (SVM).SHI denotes suitable area index.

Figure 3 .
Figure 3. Location of study site in the Ulan Buh desert (the top of the figure is the northwest arid area in China; the bottom left is the study area located in Ulan Buh desert; and the bottom right is the plot of Ulan Buh desert.The number of rows are 0,1,2,3,…and the columns are A,B,C,…in the proper sequence).

Figure 4 .
Figure 4.The Unmanned Aerial Vehicle (UAV) Light Detection and Ranging (LiDAR); (a) is the UAV LiDAR; (b) is the digital elevation model (DEM), and legend describes the elevation of the point.

Figure 5 .
Figure 5.The original and rescaled LiDAR data with different resolution, (a) is LiDAR DEM with a resolution of 0.15 m, (b) is LiDAR DEM with a resolution of 10 m, and the purple points indicate the occurrence of Nitraria tangutorum Bobr.

Figure 6 .
Figure 6.Maps of the predictor variables.The horizontal and vertical coordinates represent latitude and longitude, respectively.The (a) indicates different aspect of the study area, (b) indicates different elevation of the study area, (c) indicates different profile of the study area, and (d) indicates different slope of the study area

( 7 )Table 2 .
g(µ) = s 0 + s 1 (X 1 ) + s 2 (X 2 ) + • • •s P (X P )n = s 0 + K(x,x i )≤−1ConclusionThis is the first attempt at modelling the relationships between the spatial distribution of desert plant species and topographical factors extracted from the LiDAR data.The results indicated a general relationship between model performance and modelled species distribution on a small scale.Topographical factors should be contingent on the modelled species distribution.The LiDAR technique has certain applications in species distribution simulations.The RF produced the highest AUC when AUC was integrated into TNR + TPR (false positive rate plus true positive rate); Manhalanobis showed the best SDMs to predict the distribution of Nitraria tangutorum Bobr. in the desert on a small scale.The suitable habitat distribution area of Nitraria tangutorum Bobr.predicted in the study was larger than the actual area.A flat land with an aspect of the south at an altitude near 1010 m was found to be the driving topography of the Nitraria tangutorum Bobr.distribution.Thus, this study provides a reliable reference for the restoration and management of desert vegetation.Receiver operation characteristics.A: true positive (TP) ; A + B: positive (P); B: false negative (FN); C + D: negative; C: false positive (FP); D: true negative (TN); TPR: True positive rate = A/(A + B) = TP/ (TP + FN) = TP/P; TNR: False positive rate = C/(C + D) = FP/(FP + TN); and Accuracy (ACC) = (A + D)/ (A + B + C + D). https://doi.org/10.1038/s41598-023-40678-5