Spatial Analysis as a Tool for Plant Population Conservation: A Case Study of Tamarix chinensis in the Yellow River Delta, China

: Saltcedar ( Tamarix chinensis ) is undergoing population declination and fragmentation due to climate change and human disturbance. The existing restoration strategies usually focus on improving the environmental conditions based on the environment–saltcedar relationship, while they ignore the role of spatial autocorrelation resulting from biological interaction and ecological processes. This oversight limits the efﬁciency and sustainability of the restoration. Here, we explored the spatial pattern of the saltcedar population in the Yellow River Delta, China, and its relationship with environmental factors, incorporating spatial autocorrelation. The plant and soil parameters were extracted by an airborne LiDAR system integrated with ﬁxed soil environment measurements. The environment–saltcedar relationship incorporating spatial autocorrelation was evaluated with different regression models. Results showed that saltcedars aggregated at small scales (2–6 m), resulting from intraspeciﬁc facilitation and wind dispersal of seeds, while intraspeciﬁc competition was responsible for the random distribution at large scales (>10 m). The long-distance dispersal of seeds through water explained the signiﬁcant positive spatial autocorrelation of saltcedars at distances up to 125 m. Consequently, resulting from intraspeciﬁc facilitation and seed dispersal, aggregation distribution and positive spatial autocorrelation within the saltcedar population improved the adaptability of saltcedar to environmental stress and thereby reduced the impact of environmental factors on the abundance of saltcedar.


Introduction
As the dominant shrubs in saltmarsh wetlands, saltcedar (Tamarix chinensis Lour.) populations make a valuable contribution in preventing seawater intrusion [1]. The spatial distribution and ecological processes of saltcedar populations are strongly influenced by soil moisture and salinity gradients in saltmarsh wetlands [2]. Due to the changeable hydrological regimes resulting from climate change and human disturbance, the soil physicochemical properties and their spatial distribution in coastal saltmarsh wetland have changed significantly [3]. The intensified soil salinization and changes in the soil salinity gradient cause serious degradation of saltcedars, including a reduction in abundance, population declination, as well as fragmentation of spatial distribution patterns [4,5]. The common methods to restore the saltcedar population include hydrological regulation, seedling transplantation, and propagule cuttage. These measures aim to provide suitable physicochemical environmental conditions for individual growth and population restoration, while ignoring the intrinsic biological interactions and ecological processes. As a tial autocorrelation analysis were both applied to characterize the spatial pattern of the saltcedar population. Then, combined with soil environmental factors, spatial autoregressive models were used to analyze the environment-saltcedar relationship incorporating the spatial autocorrelation of saltcedars. These results were jointly used to analyze the generation process and mechanisms driving the spatial patterns of the saltcedar population at different spatial scales. By understanding how saltcedars distribute at different spatial scales, it would be expected to guide the restoration of plant populations and biodiversity conservation in coastal saltmarsh wetlands.

Study Area Description
Saltcedars are mainly distributed in the Yellow River Delta, China (117 • 31 -119 • 18 E, 36 • 55 -38 • 16 N). The delta is flat (<15 m above sea level) and covers a large area (>5400 km 2 ). Saltcedars, together with other typical halophytes (e.g., Phragmites australis, Suaeda heteroptera, and Imperata cylindrica), form various community types. These communities are important in wind break, dike consolidation, climate regulation, biodiversity maintenance, and community succession in the delta areas.
Aerial surveys for areas where saltcedars were typically distributed were carried out to obtain DEM, orthophoto, geographic coordinates, and information about plant parameters (abundance, density, tree height, and crown diameter) during October 2018 using an UAV-LiDAR system (Figure 1). The grid dividing method was used to divide the entire area into 35 quadrats. Topsoil samples (~−20 cm) were collected using a five-point sampling method in every quadrat to measure soil water and salt content. Soil samples were dried at 60°C until a constant mass was reached in an oven. The samples were weighed before and after being dried. The difference in these values was used to calculate the soil moisture. Dried soil samples were ground and passed through a 35-mesh screen (Φ = 20 cm). The samples were then mixed within deionized water (mass ratio, soil:water = 1:5) in a 50 mL centrifugal tube. After shaking for 0.5 h, the supernatant was left to stand for 24 h. Salinity was then measured using an electronic meter (Bante 540-DH, Shanghai Bante Instruments Co., Ltd., Shanghai, China) [3].

Photogrammetry Workflow
A six-rotor-wing unmanned aerial vehicle (UAV, DJI Matrice 600 pro) equipped with a mini-LiDAR system (Genius V+Rfans, SureStar Laser Technology Co.Ltd., Beijing, China) and a Sony (A7R) compact digital camera with a 48-megapixel image sensor was used to collect point clouds and aerial photographs of the study area ( Figure 1). The drone was flown in automated survey mode, whereby it followed a pre-programmed global positioning system (GPS) path in DJI GO application software, to avoid human piloting error and to achieve a consistent forward and side overlap of ≥80% between aerial images. The UAV was flown at 50 m altitude with a velocity of 4 m/s. In the meantime, the camera triggered a photo every 2 s. The aerial survey generated 1177 photographs and point clouds for the entire targeted area. We then recorded the centroid coordinates of each photograph as x, y, and z using a novel real-time kinematic-global positioning system (RTK-GPS) with centimeter-level relative positioning accuracy. We checked for system error by conducting post-processing of data of position and orientation system (POS, including differential global positioning system (DGPS) and inertial measurement unit (IMU)) with POSPac UAV 8.2.1 software. The flight trajectories of UAV and exterior orientation elements of photographs were calculated as well. By combining trajectories and laser ranging data, the original point cloud data (.Las) were obtained with UI_v3.9.3_Rfans software. Throughout the entire process, the different coordinate systems (IMU coordinate system, laser scanning coordinate system, UAV coordinate system, navigation coordinate system, and geocentric coordinate system) were unified to the World Geodetic System-1984 (WGS-1984) coordinate system. The original point clouds were classified automatically by TerraSolid software (TerraScan and TerraModeler modules) to establish DEM and contour lines of the study area by separating ground points. This software was also used to separate different vegetation types based on tree height. Individual plant segmentation and identification was processed by PointCloudCatalyst software. By combining the original photos with the exterior orientation element of photos, a digital orthophoto map (DOM) of the study area was integrated with Pix4Dmapper 2.0 software after manually removing images during takeoff and landing.

Photogrammetry Workflow
A six-rotor-wing unmanned aerial vehicle (UAV, DJI Matrice 600 pro) equipped with a mini-LiDAR system (Genius V+Rfans, SureStar Laser Technology Co.Ltd., Beijing, China) and a Sony (A7R) compact digital camera with a 48-megapixel image sensor was used to collect point clouds and aerial photographs of the study area ( Figure 1). The drone was flown in automated survey mode, whereby it followed a pre-programmed global positioning system (GPS) path in DJI GO application software, to avoid human piloting error and to achieve a consistent forward and side overlap of ≥80% between aerial images. The UAV was flown at 50 m altitude with a velocity of 4 m/s. In the meantime, the camera triggered a photo every 2 s. The aerial survey generated 1177 photographs and point clouds for the entire targeted area. We then recorded the centroid coordinates of each photograph as x, y, and z using a novel real-time kinematic-global positioning system (RTK-GPS) with centimeter-level relative positioning accuracy. We checked for system error by conducting post-processing of data of position and orientation system (POS, including differential global positioning system (DGPS) and inertial measurement unit (IMU)) with POSPac UAV 8.2.1 software. The flight trajectories of UAV and exterior orientation elements of photographs were calculated as well. By combining trajectories and laser ranging data, the original point cloud data (.Las) were obtained with UI_v3.9.3_Rfans software. Throughout the entire process, the different coordinate systems (IMU coordinate system, laser scanning coordinate system, UAV coordinate system, navigation coordinate system, and geocentric coordinate system) were unified to the World Geodetic System-1984 (WGS-1984) coordinate system. The original point clouds were classified automatically by TerraSolid software (TerraScan and TerraModeler The geographic coordinate information of saltcedars was extracted from individual plant segmentation and identification, and was manually checked to eliminate invalid points. Invalid points included obvious mismatches between coordinate points and plant position in DOM, or non-target plants. Finally, 3083 effective saltcedar trees were confirmed. To further analyze the relationship between environmental factors (soil moisture (%), soil salinity (ppt), and ground elevation (m)) and the spatial distribution pattern of the saltcedar population, the entire study area was split into 35 grids using ArcGIS 10.6 (see Section 2.1). For each grid, information on saltcedars was extracted, including abundance, density (trees per 100 m 2 ), average tree height (m), and average crown diameter (m). The geographical coordinates of each grid had a one-to-one correlation with the coordinates of the soil sampling quadrats. After clipping the edge of the study area regularly, data of soil and saltcedar features obtained in the 33 quadrats were used in the following spatial analysis.
A single-sample Kolmogorov-Smirnov test was used to test the normality of raw data in SPSS Statistics 23.0. Logarithmic transformation was used to ensure that the data followed a normal distribution when raw data were not normally distributed. All data were expressed as mean ± standard error.

. Analysis of Spatial Distribution Patterns
Second-order statistics methods (Ripley's K function and the pair-correlation function g(r)) were applied to analyze the spatial distribution pattern of the saltcedar population, based on the spatial geographical coordinates of every saltcedar tree (3083 trees in total). Ripley's K and g(r) function were defined as: A heterogeneous Poisson process (HPP) was adopted as a null hypothesis model when calculating the g(r) function. Gaussian kernel function (bandwidth R = 30.00 m) was used to estimate the intensity, and to compare patterns. We performed 199 Monte Carlo simulations, and selected an estimate of 95% simulation envelopes. The maximum distance was half of the shortest side of the plot (400 m), and the step length was 1 m. If the function value calculated from the realistic population distribution data was above the envelope, it represented an aggregated distribution of saltcedar population; if the real data fell within the envelope, they represented a random distribution of saltcedar population; otherwise, they represented a regular distribution. Spatial point pattern analysis was conducted in Programita 2014 software [27,34].

Spatial Autocorrelation Analysis
The spatial autocorrelation index (Moran's I) was used to evaluate the degree of spatial correlation of the saltcedar population and environmental factors. Moran's I reflects the similarity of the attribute values of adjacent positions. The spatial weight matrix was constructed based on queen contiguity (order of contiguity = 1). Global spatial autocorrelation was used to determine whether aggregation features existed for plants in a targeted area; however, the type and location of the spatial pattern needed to be further analyzed by local indicators of spatial association (LISA). Global and local Moran's I were measured and defined, respectively, as: I generally ranges between [−1 and 1]; I < 0 represents a negative autocorrelation, I = 0 represents no correlation, and I > 0 represents a positive autocorrelation. Usually, a significance test of Moran's I is performed with a Z score to judge the significance of spatial autocorrelation. When the Z score is greater than 1.96 or less than −1.96 (a = 0.05), it indicates a significant spatial autocorrelation; when the Z score is in the interval [−1.96, 1.96], it indicates an independent random distribution. Global and local spatial autocorrelation of soil and saltcedar features in the study area were calculated with Geoda 1.16 software.
Autocorrelation coefficients were calculated separately against different distance classes and plotted as spatial correlograms [35]. These nondirectional correlograms indicated the degree of association between values of a variable at different spatial scales. The significance of Moran's I coefficients in correlograms was calculated by comparing them with those obtained from the 199 Monte Carlo simulations. The significance of an entire correlogram was usually calculated using a Bonferroni correction [36]. Spatial correlograms of soil and saltcedar features were calculated in PASSaGE 2.0 (PASSaGE: Pattern Analysis, Spatial Statistics and Geographic Exegesis) [37].

Spatial Regression Analysis
The abundance of a plant is an important quantitative index to reflect the relationship between plant distribution and environmental factors [38]. We used the abundance of saltcedar as a dependent variable, while other variables (including soil moisture, soil salinity, elevation, tree density, crown diameter, and tree height of saltcedars) were used as explanatory variables. To eliminate the influence of dimensions and magnitude between variables, all variables were standardized as zero-mean normalization in SPSS Statistics 23.0 before the regression analysis. Data standardization was calculated as x is the mean value of all sample data, and σ is the standard deviation of all sample data.
To compare the effectiveness in predicting the spatial distribution of the saltcedar population, two spatial autoregressive models (spatial lag model (SLM) and spatial error model (SEM)) and one classical multiple linear regression model (estimated by ordinary least square (OLS)) were selected. SLM assumes that the spatial autocorrelation of the response variable is caused by an inherent property, whereas SEM assumes that the autoregressive process is found only in the error term because of inherent or induced spatial autocorrelation [19]. The general forms of these regression models were defined as: where w is the weighted matrices according to the contiguity of spatial elements, ρ and λ are the regression parameters to be estimated, β is the regression coefficient of explanatory variables, ε is the spatially dependent error term, η~N (0, σ 2 I n ), η is the random error term, σ 2 is the variance of η, and I is the unit matrix. Spatial autocorrelation of the residuals of OLS, SLM, and SEM was also tested. The goodness of fit of the regression models was examined by R 2 , maximum likelihood logarithm (LIK), Akaike's Information Criterion (AIC), and Schwartz index (SC). Higher R 2 and LIK (or lower AIC and SC) indicated the strong explanatory ability of these models. It indicated a significant difference in the goodness of fit between two models if the AICs between the two models differed by more than three. Spatial regression analyses were calculated in Geoda 1.16. 40% of all saltcedars, respectively. For all 3083 trees, the average tree height and average crown diameter of saltcedars were 1.84 ± 0.01 m and 1.18 ± 0.01 m, respectively. Average soil salinity, average soil moisture, and average ground elevation were 1.87 ± 0.12 ppt, 21.19 ± 0.26%, and 1.23 m, respectively. Furthermore, the result of the Kolmogorov-Smirnov test showed that the original data (soil properties and saltcedar features) followed a normal distribution (p > 0.05). The coefficient of variation (CV) reflected the spatial variability of soil properties and saltcedar features. The CVs of soil salinity, density, and abundance of saltcedar (36%, 47%, and 36%, respectively) showed strong variability (CV > 30%). Maps of soil and saltcedar features based on the inverse distance weighting interpolation method also demonstrated a heterogeneous spatial distribution ( Figure 3). and average crown diameter of saltcedars were 1.84 ± 0.01 m and 1.18 ± 0.01 m, respectively. Average soil salinity, average soil moisture, and average ground elevation were 1.87 ± 0.12 ppt, 21.19 ± 0.26%, and 1.23 m, respectively. Furthermore, the result of the Kolmogorov-Smirnov test showed that the original data (soil properties and saltcedar features) followed a normal distribution (p > 0.05). The coefficient of variation (CV) reflected the spatial variability of soil properties and saltcedar features. The CVs of soil salinity, density, and abundance of saltcedar (36%, 47%, and 36%, respectively) showed strong variability (CV > 30%). Maps of soil and saltcedar features based on the inverse distance weighting interpolation method also demonstrated a heterogeneous spatial distribution ( Figure 3).

Spatial Distribution Patterns of Saltcedars
Second-order statistics showed that Ripley's K function described a successive random-clustered-regular-random distribution with scale increasing. In comparison, the pair-correlation function g(r) showed a successive regular-clustered-random distribution (Figure 4). The spatial distribution pattern of the saltcedar population was clustered at a distance of 2-6 m and random above 10 m. The saltcedar population presented the strongest aggregation effect at the distance of 6 m.

Spatial Distribution Patterns of Saltcedars
Second-order statistics showed that Ripley's K function described a successive randomclustered-regular-random distribution with scale increasing. In comparison, the paircorrelation function g(r) showed a successive regular-clustered-random distribution (Figure 4). The spatial distribution pattern of the saltcedar population was clustered at a distance of 2-6 m and random above 10 m. The saltcedar population presented the strongest aggregation effect at the distance of 6 m.

Spatial Distribution Patterns of Saltcedars
Second-order statistics showed that Ripley's K function described a successive random-clustered-regular-random distribution with scale increasing. In comparison, the pair-correlation function g(r) showed a successive regular-clustered-random distribution (Figure 4). The spatial distribution pattern of the saltcedar population was clustered at a distance of 2-6 m and random above 10 m. The saltcedar population presented the strongest aggregation effect at the distance of 6 m.

Spatial Autocorrelation Analysis of Soil and Saltcedar Variables
Both soil environmental factors and saltcedar features were significantly and positively spatially autocorrelated (Table 1). In this study, the longest boundary of the study

Spatial Autocorrelation Analysis of Soil and Saltcedar Variables
Both soil environmental factors and saltcedar features were significantly and positively spatially autocorrelated (Table 1). In this study, the longest boundary of the study area was 685 m and was divided into ten distance classes. The upper bound of each distance class was 125, 178, 221, 256, 305, 338, 391, 443, 502, and 685 m, successively. However, due to the large breadth of the class, the coefficients in the largest distance class were often unreliable [35]. After removing unreliable coefficients in the largest distance class, the correlograms of soil properties and saltcedar features at different distance scales showed that spatial autocorrelation decreased with the increasing distance ( Figure 5). The spatial autocorrelation of elevation, soil salinity, and moisture all peaked in the 0-125 m distance class, indicating a significant positive spatial autocorrelation between closer points at 0-125 m. The spatial autocorrelation of the saltcedar population changed noticeably, especially for the abundance of saltcedar. Overall, soil properties and the saltcedar population showed a positive significant spatial autocorrelation up to 125 m.
Sustainability 2021, 13, x FOR PEER REVIEW 9 of 16 area was 685 m and was divided into ten distance classes. The upper bound of each distance class was 125, 178, 221, 256, 305, 338, 391, 443, 502, and 685 m, successively. However, due to the large breadth of the class, the coefficients in the largest distance class were often unreliable [35]. After removing unreliable coefficients in the largest distance class, the correlograms of soil properties and saltcedar features at different distance scales showed that spatial autocorrelation decreased with the increasing distance ( Figure 5). The spatial autocorrelation of elevation, soil salinity, and moisture all peaked in the 0-125 m distance class, indicating a significant positive spatial autocorrelation between closer points at 0-125 m. The spatial autocorrelation of the saltcedar population changed noticeably, especially for the abundance of saltcedar. Overall, soil properties and the saltcedar population showed a positive significant spatial autocorrelation up to 125 m.  The coefficient in the largest distance class was often unreliable, due to the large breadth of the class and/or paucity of relevant pairs of points (depending on how the classes were designed); thus, this value was often not plotted in correlogram analyses [35]. All correlograms were globally significant (p < 0.05).

Quantification of Factors Influencing the Distribution of Saltcedars
Stepwise regression was used to remove explanatory variables with high VIF values (>7.5) from the regression model. The OLS model between the abundance of saltcedars and explanatory variables was delineated as: YZabundance = −0.01 + 0.37 XZmoisture −0.26 XZsalinity Figure 5. Nondirectional spatial correlogram on the features of (A) soil and (B) saltcedar (Tamarix chinensis). Solid circles were significant at p < 0.05, while asterisks were nonsignificant. The upper bound of each distance class was 125, 178, 221, 256, 305, 338, 391, 443, and 502 m, successively. The coefficient in the largest distance class was often unreliable, due to the large breadth of the class and/or paucity of relevant pairs of points (depending on how the classes were designed); thus, this value was often not plotted in correlogram analyses [35]. All correlograms were globally significant (p < 0.05).

Quantification of Factors Influencing the Distribution of Saltcedars
Stepwise regression was used to remove explanatory variables with high VIF values (>7.5) from the regression model. The OLS model between the abundance of saltcedars and explanatory variables was delineated as: Y Zabundance = −0.01 + 0.37 X Zmoisture −0.26 X Zsalinity + 0.57 X Zcrown , (R 2 = 0.52, p < 0.01). Using the same variables in the OLS, the environment-saltcedar relationship was also delineated through SLM and SEM ( Table 2). The contribution of soil moisture and soil salinity to the abundance of saltcedar was reduced by 11.67% and 10.04% in the SLM and by 13.57% and 4.77% in the SEM, respectively. Pairs of predicted and observed values for the abundance of saltcedars in the three regression models were all discretely distributed around the 1:1 reference line ( Figure 6). Spatial autocorrelation and correlograms of residuals showed no spatial autocorrelation in the residuals of OLS, SLM, and SEM ( Figure 6). Thus, these models were suitable for evaluating environment-saltcedar relationships. Overall, there was no significant difference in the prediction performance of the three regression models for the abundance of saltcedar. However, comparison of the goodness of fit among three regression models indicated that R 2 in SEM was greater than that in SLM and OLS. Furthermore, the AIC values of SEM were the lowest among the three regression models, in which a value of SEM differed by over three from the values in OLS and SLM ( Table 2). The goodness of fit showed that, compared with OLS and SLM, SEM could better explain the variation in the abundance of saltcedars. Overall, SEM had more advantages in explaining the relationships between environmental factors and saltcedar population distribution in saltmarsh. was greater than that in SLM and OLS. Furthermore, the AIC values of SEM were the lowest among the three regression models, in which a value of SEM differed by over three from the values in OLS and SLM ( Table 2). The goodness of fit showed that, compared with OLS and SLM, SEM could better explain the variation in the abundance of saltcedars. Overall, SEM had more advantages in explaining the relationships between environmental factors and saltcedar population distribution in saltmarsh. Solid circles indicate the significant spatial autocorrelation (p < 0.05) and asterisks indicate the nonsignificant spatial autocorrelation. All correlograms were globally significant (p < 0.05). Solid circles indicate the significant spatial autocorrelation (p < 0.05) and asterisks indicate the nonsignificant spatial autocorrelation. All correlograms were globally significant (p < 0.05).

Spatial Pattern of Saltcedar Population and Mechanisms Driving Its Formation
Our results showed that the saltcedar population exhibited an aggregated distribution at 2-6 m scales, and a random distribution at large scale (>10 m). Previous studies also reported a similar spatial distribution pattern of the saltcedar population [12,14,39,40]. The aggregation patterns of plant populations at small scales are often the result of intraspecific facilitation and local seed dispersal. Seeds of saltcedar are light and small, with hairs on the epidermis, making them highly suited to dispersal to different distances by wind or water. A previous study reported that saltcedar populations could disperse around 2.5 km yr −1 (approximately 6.85 m d −1 ) by wind and 11 km yr −1 (approximately 34.25 m d −1 ) along rivers and other waterways [41]. For wind-dispersed seeds, they usually germinate within 24 h in moist soil [42]. The intensive dispersal and germination in the short term resulted in aggregated distribution in small areas. Therefore, the aggregated distribution at 2-6 m might correspond to the distance at which the wind dispersal of seeds influenced the spatial distribution of saltcedar populations.
Biotic interaction is another important factor driving the spatial patterns of plant populations at different scales. Generally, the increasing abundance of plants at small scales promotes a balance between aboveground and underground biomass [43] or branch and leaf growth [44]. These balances contribute towards regulating the net photosynthetic and transpiration rates of plants, and benefiting their adaptation in the heterogeneous soil environments. Furthermore, the densely packed distribution of halophytes results in habitat amelioration (i.e., decline in soil salinity and increase in soil water content), enhancing seed germination and seedling growth [45,46]. As the intensity of environmental stress declines, the intraspecies facilitation weakens while intraspecific competition increases [47,48]. This phenomenon is the result of increased demand for resources by individual plants, which intensifies competition among individuals [49]. Moreover, with the increasing resistance of individual to soil environmental stresses, the interdependence among saltcedars declines. As a result of the reducing interdependence and increasing biotic competition, the plant population presents random or uniform distribution at a large scale [50]. Due to the transformation of these biological interactions, the spatial distribution pattern of the saltcedar population gradually changed from the aggregating pattern at small scales (2-6 m) to the random distribution at large scales (>10 m).

Origin of Spatial Autocorrelation in Saltcedar Population
Spatial autocorrelation is considered to be a general attribute of ecological variables and can be found at any scale in the ecosystem [24]. Spatial autocorrelation of plant populations is caused by the dependence of spatially autocorrelated environmental variables (induced or exogenous spatial autocorrelation) or by the inherent biological processes of the variable itself (inherent or endogenous spatial autocorrelation) [23,51]. In this study, results showed that soil environmental factors and saltcedar features both had positively significant spatial autocorrelation ( Table 1). The saltcedar population had a positive spatial autocorrelation at distances up to 125 m ( Figure 5). This indicated that saltcedar trees within a 125 m distance tended to have similar occupancy. For saltcedar populations, their spatial autocorrelation is not only caused by the spatially autocorrelated soil salinity and moisture, but also related to the ecological processes [23]. Dispersal is a major factor driving endogenous spatial autocorrelation [52]. A previous study reported that saltcedar seeds could germinate in water, even when afloat, and that seedlings could survive submerged for a few weeks [42]. Thus, the viability of the water-dispersed seeds of saltcedar could be maintained for long periods and at great distances. Even when germinating in water, seedlings could continue dispersing along rivers or other waterways until reaching suitable conditions for colonization, such as open sunny ground with low competition [42]. Therefore, the largest scale of saltcedars presenting a positively significant spatial autocorrelation (125 m) might correspond to the distance at which long-distance seed dispersal by water influenced population distributions.
The regression analysis in this study showed that the abundance of saltcedars was influenced significantly by soil salinity and moisture rather than elevation, which was consistent with previous studies [1,53]. However, only considering soil conditions could not fully explain the abundance of saltcedar. Crown diameter has a closely positive correlation with the abundance of saltcedars. This may be because the increasing crown diameter is conducive to maintaining soil moisture and providing a shading condition for the growth of seedlings, which can be regarded as intraspecific promotion [17]. These results reflected the existence of intraspecific facilitation within the saltcedar population. Moreover, it is noteworthy that the classical multiple linear regression model (OLS) in this study could not describe the environment-saltcedar relationship very well (R 2 = 0.52). Two potential factors might have led to this result. One factor was that the OLS model ignored the spatial autocorrelation, which was an important factor in the spatial distribution of the saltcedar population. By comparison, SEM, which incorporated the spatial autocorrelation of saltcedars, presented a better explanatory ability in evaluating the relationship between soil environmental factors and plant spatial distribution. This was consistent with the result of other studies [25,26,31,54]. Another factor might be due to the accuracy of individual plant segmentation and identification. In this research, different vegetation types were separated based on tree height. Some young trees with low height might have been ignored or classified into other types. This might also have resulted in the unsatisfactory correlation coefficient. In the future, an aerial survey of vegetation types, and a combination of a UAV-LiDAR system with multispectral imaging technology or simultaneous localization and mapping (SLAM) technology, would promote the quality of point clouds and the accuracy of individual plant segmentation and identification [55]. Other effective measures include increasing the density of fixed monitoring sites, prolonging the monitoring time, and enriching the monitoring soil elements.
Overall, as a result of intraspecific facilitation and seed dispersal, the saltcedar population was aggregated at 2-6 m scales and presented a positive spatial autocorrelation at distances up to 125 m. The existence of positive spatial autocorrelation in the saltcedar population led to overestimation of environmental factors when analyzing the environmentsaltcedar relationship by the traditional OLS method. This study highlighted the importance of spatial autocorrelation when analyzing the environment-saltcedar relationships. Future study should cautiously estimate the influence of soil environmental factors on the spatial distribution of saltcedar populations.

Implications for Ecological Management and Restoration of Saltcedar Population
Currently, saltcedar populations, as well as other halophytes in coastal saltmarsh, are undergoing increasingly serious degradation and fragmentation because of the multiple impacts of global climate change and human disturbance [4,5]. Generally, spatial pattern analysis plays an instructional role for the restoration and management of plant populations. Based on life history strategies and distribution patterns of species, different implementation approaches are adopted for population restoration. Previous studies suggest that restoration efforts would be more effective for trees in forests and corals with an r-selected life history if organisms were added to landscapes in uniform and gridded patterns [56,57]. Conversely, for the plants with a competitive disadvantage or inhabiting harsh environments, aggregated arrangement could be more effective in increasing flowering individuals and reproductive biomass, which would contribute to maintaining the plant population in the following generations [58].
In this study, the results suggested an aggregation distribution of the saltcedar population at small scales and a positive significant relationship between the abundance and crown diameter of saltcedars. Consistent with previous studies, these results implied that saltcedar individuals could be transplanted as patches with scopes of 2-6 m. Aggregation of saltcedar individuals would promote the survival rates and colonization ability after transplantation into a new habitat [39]. Some engineering measures in favor of hydrological connectivity are also recommended to restore the saltcedar population through assisting seed dispersal. For example, ditch excavation provides channels for seed spreading with water [59]. Modifying microtopography contributes to seed interception and provides a surface runoff channel for seed dispersal [60,61].
In addition, the UAV-LiDAR system shows good application prospects in topographic mapping, plant population monitoring, and pattern analysis in saltmarsh. As an emerging remote sensing technology, the UAV-LiDAR system allows the efficient survey of target areas at different spatial scales, especially local-scale monitoring [62]. UAV-LiDAR systems overcome the disadvantages of traditional remote sensing approaches, which are not able to meet the time frequency requirements in analyzing ecological processes and the required high spatial resolution in analyzing fine-scale information [63]. Future research on the spatial pattern analysis and restoration of plant populations in saltmarsh should be conducted over the long term and at broad scales, combining UAV-LiDAR systems with diverse monitoring methods, such as remote sensing technologies and multispectral or hyperspectral imaging techniques. Such information would assist in understanding the distribution patterns of plant populations at different spatiotemporal scales and guide theoretically their ecological management and restoration in saltmarsh.

Conclusions
This study extracted plant and soil parameters by integrating the UAV-LiDAR system with fixed soil environment measurements, and further analyzed the spatial pattern of the saltcedar population and the environment-saltcedar relationship, incorporating the fine-scale spatial autocorrelation. The results showed that saltcedars presented clustered distribution at small scales (2-6 m) due to intraspecific facilitation and wind dispersal of seeds. In comparison, intraspecific competition was responsible for the random distribution of saltcedars at large scale (>10 m). The saltcedar population had a positive significant spatial autocorrelation, with Moran's I peaking at 0-125 m, which was related to water-dispersed seeds. After incorporating the spatial autocorrelation of saltcedars when analyzing the environment-saltcedar relationship, the interpretation ability of soil moisture and salinity for the abundance of saltcedars decreased by 11.67% and 10.04% in SLM, and by 13.57% and 4.77% in SEM, respectively. Meanwhile, there was a significant positive correlation between the abundance and crown diameter of saltcedars, which had a stronger effect on the abundance of saltcedar (>55%) than soil factors.
Consequently, these results suggested that due to intraspecific facilitation and seed dispersal, the saltcedar population was clustered and autocorrelated spatially at different scales. This spatial pattern improved the adaptability of saltcedars to environmental stress and thus reduced the impact of environmental factors on the abundance of saltcedar. Therefore, during the restoration of the saltcedar population in the future, biotic interactions and seed dispersal should be given more attention, in addition to the existing measures (e.g., improving environmental factors). Transplanting saltcedars as patches with appropriate size and constructing channels to facilitate seed dispersal are both recommended measures to restore and manage the saltcedar population in saltmarsh. These measures will profit the restoration and management of the saltcedar population as well as other halophytes in coastal saltmarsh.