Using UAV LiDAR to Extract Vegetation Parameters of Inner Mongolian Grassland

: The accurate estimation of grassland vegetation parameters at a high spatial resolution is important for the sustainable management of grassland areas. Unmanned aerial vehicle (UAV) light detection and ranging (LiDAR) sensors with a single laser beam emission capability can rapidly detect grassland vegetation parameters, such as canopy height, fractional vegetation coverage (FVC) and aboveground biomass (AGB). However, there have been few reports on the ability to detect grassland vegetation parameters based on RIEGL VUX ‐ 1 UAV LiDAR ( Riegl VUX-1) systems. In this paper, we investigated the ability of Riegl VUX ‐ 1 to model the AGB at a 0.1 m pixel resolution in the Hulun Buir grazing platform under different grazing intensities. The LiDAR ‐ derived minimum, mean, and maximum canopy heights and FVC were used to estimate the AGB across the entire grazing platform. The flight height of the LiDAR ‐ derived vegetation parameters was also analyzed. The following results were determined: (1) The Riegl VUX ‐ 1 ‐ derived AGB was predicted to range from 29 g/m 2 to 563 g/m 2 under different grazing conditions. (2) The LiDAR ‐ derived maximum canopy height and FVC were the best predictors of grassland AGB (R 2 = 0.54, root ‐ mean ‐ square error (RMSE) = 64.76 g/m 2 ). (3) For different UAV flight altitudes from 40 m to 110 m, different flight heights showed no major effect on the derived canopy height. The LiDAR ‐ derived canopy height decreased from 9.19 cm to 8.17 cm, and the standard deviation of the LiDAR ‐ derived canopy height decreased from 3.31 cm to 2.35 cm with increasing UAV flight altitudes. These conclusions could be useful for estimating grasslands in smaller areas and serving as references for other remote sensing datasets for estimating grasslands in larger areas.


Introduction
Grasslands are an important part of the land surface ecosystem and account for approximately 27% of the global surface [1][2][3][4]. Grassland aboveground biomass (AGB) is a significant grassland productivity indicator informative for both engineers and grassland managers [5,6]. Moreover, AGB drives ecosystem processes via vegetation properties and biodiversity across multiple trophic levels in northern Sweden [7]. As one of the most important indices of grasslands, AGB can be used to evaluate carbon cycling and the net primary productivity for grasslands [8][9][10]. AGB can be effectively estimated using vegetation structure data, including canopy height and fractional vegetation coverage (FVC). FVC includes the ground layer (grasses, leaf litter, moss and lichen), mid layer (mid-sized trees, mid-storey cover) and canopy layer (hollow-bearing trees, canopy depth) as well as attributes capturing vegetation conditions (tree dieback, mistletoe) [11].
Traditional AGB field sampling efficiency is time consuming and costly. Biomass is traditionally measured by quadrats that were not defoliated by cutting or grazing during the growing season [12]. Although the data are accurate, data collection is time consuming, costly and laborious [13], and accurately acquiring grassland AGB data is difficult in larger areas. Passive optical remote sensing is an advanced grassland ecosystem monitoring tool in the management of grassland structure research on parameters such as vegetation coverage, AGB and vegetation growth [14][15][16]. L-band passive microwaves (1-2 GHz) have a strong ability to penetrate vegetation due to their long wavelength spectral range (15-30 cm) [17]. However, passive optical remote sensing suffers from saturation problems when estimating grassland biomass, especially in areas with high vegetation cover. The saturation problem of the vegetation index is characterized by lower accuracy in the vegetation index because the spectral bands (6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18) weakly penetrate through the vegetation layer compared with AGB accumulation when the grassland vegetation density reaches a certain threshold. In addition, UAV optical cameras have been used to improve the estimation accuracy of vegetation biomassbased structures from motion (SfM) algorithms and canopy height models (CHMs) [18]. In recent years, plant structure parameters (such as plant height) have become the focus of UAV-based remote sensing approaches for crop monitoring. The canopy surface model has a robust estimator of crop biomass. However, monitored arable crops typically suffer differences in floristic composition, and the co-occurrence of phenology leads to high spatial-temporal heterogeneity [19]. The emergence of this heterogeneity presents a great challenge for remote sensing to detect biomass. Electromagnetic radiation information received by remote sensing cannot accurately reflect a change in biomass, and remote sensing models cannot accurately estimate the distribution area of high biomass accumulation.
Light detection and ranging (LiDAR) acquires three-dimensional information (x, y, z) and emits a beam of light pulses, hits the object, and reflects it back; this information is ultimately received by the LiDAR sensor. Combined with a specified laser height and laser scanning angle, the laser position is obtained from global positioning systems and inertial navigation systems. As a method of 3D data acquisition, LiDAR technology has been widely used in agriculture [15,20], forestry [21][22][23] and many other industries [24,25]. LiDAR data offer a unique opportunity to measure three-dimensional vertical structures. For example, canopy height measurements rely on the laser pulse return from the vegetation and ground in different physical objects and different elevations. In recent years, unmanned aerial vehicle (UAV) LiDAR has also been used to monitor plant structures [26]. Owing to their low weight and cost, UAV LiDAR has been widely used in vegetation monitoring. Xiangqian Wu et al. [27] used a Velodyne LiDAR sensor to detect individual trees and estimate forests in China. Wang et al. [28] used Velodyne's HDL-32E UAV LiDAR system to study the Hulun Buir grassland ecosystem. In summary, Velodyne's LiDAR sensor can effectively extract vegetation parameters in forests and grasslands and rapidly evaluate AGB.
From the perspective of soil and ecology science, some researchers have carried out experiments on the Hulun Buir grazing platform to study the spatial distribution characteristics of soil carbon, nitrogen and grassland species under different grazing intensities [29][30][31]. However, there is still a lack of research on remote sensing detection of grassland vegetation on the Hulun Buir grazing platform. In this paper, we obtain RIEGL VUX-1 UAV LiDAR (Riegl VUX-1)-derived grassland vegetation structures to model the grassland AGB of the Hulun Buir grazing platform. This paper evaluates the ability of the Riegl VUX-1 sensor to extract grassland parameters. Using LiDAR data to obtain the canopy height, FVC and AGB prediction maps with a 0.1 m resolution under different grazing intensities are evaluated. We also analyze different flight heights of LiDAR-derived vegetation structures in the Hulun Buir grazing platform. The relevant conclusions can provide a theoretical reference for steppe grasslands under different grazing conditions in Inner Mongolia.

Study Area
The Hulun Buir grazing platform was established in 2009. The stocking rates were set to 0, 0.23, 0.34, 0.46, 0.69, and 0.92 AU ha −1 , which are referred to as G0, G1, G2, G3, G4, and G5, respectively [32] (where 1 AU = 500 kg of adult cattle) ( Figure 1). The area is located 8 km east of Hailaer, Inner Mongolia, China (49°20′24″ N, 119°59′44″ E). The area of each plot was 300 m × 160 m, and the total area was 90 hm 2 . In the Hulun Buir grazing platform area, the altitude is approximately 666-680 m, and the terrain slope is less than 3°. The average annual rainfall is between 350 mm and 400 mm, and approximately 80% of rainfall occurs during the growing season [28]. Grazing lasts approximately 120 days (from 1 June to 1 October) annually [33]. The vegetation characteristics in this area are those of typical steppe. The dominant vegetation species are Leymus chinensis, Stipa capillata, Carex pediformis, Filifolium sibiricum, Bupleurum scorzonerifolium, Galium verum and Pulsatilla turczaninovii. First, we used TerraScan (Terrasolid Ltd., Helsinki, Finland) for module classification in the Terrasolid software package to preprocess the LiDAR data. LiDAR noise points were eliminated, and the LiDAR point cloud data were filtered and normalized. Second, we used LiDAR data divided into vegetation points, ground points and other points. The ground points were then used to generate a digital terrain model (DTM), and the vegetation points were then used to generate a digital surface model (DSM). The DSM minus DTM was gridded without interpolation to establish a canopy height model (CHM) with a 0.1 m spatial resolution. Third, we used LiDAR data to obtain the LiDAR-derived minimum canopy height, mean canopy height, maximum canopy height and FVC from CHM. Fourth, we analyzed a simple regression model between the LiDAR-derived vegetation structures (the minimum, mean and maximum canopy height and FVC) and ground-measured data with all training samples. Finally, we chose the best model to obtain grassland prediction maps at the plot scale in the grazing platform. The workflow is shown in Figure 2. The LiDAR workflow was similar to that in other studies [34,35] and is described below.

Ground Measurement Data Collection
Ground measurements started on 11 September 2018, with observations of the vegetation parameters (canopy height, FVC and AGB) of this steppe. As shown in Figure  2, six sample points were selected under each plot. Among them, we selected four groundmeasured samples as model training samples, and two ground-measured samples were used as model validation samples. Ninety-six ground-measured 50 cm x 50 cm quadrats (in the W2 and W4 plots, the number of ground-measured samples was not sufficient; thus, these plots were not included in the model calculation) were established following a line transecting each plot [36].
The ground-measured grassland canopy height and FVC and AGB measurements were carried out in a 50 cm × 50 cm field quadrat. Nine measurement points were evenly selected throughout the quadrat. We recorded the mean, maximum and minimum values of the natural height of the grassland with a ruler. The FVC was ground-measured using an ordinary digital camera (Canon PowerShot G7 X, Tokyo, Japan). The manual threshold (10-30%) method was used to obtain vegetation pixels from camera red-and green-band images perpendicular to the ground at a height of 1 m. After identifying and counting the green vegetation pixels using these values [37], the AGB was then clipped at the ground level and dried at 65 °C for 48 h to a constant weight [30]. We also used the real-time kinematic product (RTK-G975, Beijing Hezhongsizhuang Technology Co., Ltd., Beijing, China). A reference station and mobile station were used to form the carrier phase, and the fixed RTK sites were used to obtain the geographic center coordinates of the 96 quadrats. The measurement accuracy was at the centimeter level, which was used to determine the central coordinate position of each quadrat.

LiDAR Data Collection
On 11 September 2018, we designed three UAV LiDAR flights for use in the study area. The UAV LiDAR data collected on the first and second flights covered 18 plots at a flight altitude of 119 m. With increasing flight altitude, the standard deviation of the LiDAR-derived canopy height decreased sharply [28]. The third UAV LiDAR flight was conducted over the W0 plot at heights of 40-110 m (at intervals of 10 m), and the flight speed was 8 m/s. Figure 3a shows the UAV flight lines. The HS-600 UAV system adopts a DJI Matrice 600, and the weight of the HS-600 UAV system is 4.5 kg on the platform, which integrates a Riegl VUX-1 sensor with an AP15(X) inertial navigation system. As shown in Figure 3b, the Riegl VUX-1 sensor can provide a measurement rate of up to 500,000 pulses/s and a large field angle of view of 330°. The LiDAR data adopted the pulse ranging principle of the Riegl VUX-1 sensor through nearinfrared laser beam scanning. The UAV flight parameters of the Riegl VUX-1 sensor are shown in Table 1.

LiDAR Data Analysis
LiDAR point data were mainly classified as vegetation, ground and other noise by the triangulated irregular network (TIN) filtering method [38], which is embedded in TerraScan (Terrasolid Ltd., Helsinki, Finland) for module classification in the Terrasolid software package to create a digital terrain model (DTM) and a digital surface model (DSM) from the LiDAR data. Classification points were assessed to ensure that most LiDAR points were correctly classified. The relative height LiDAR data were gridded without interpolation to establish a canopy height model (CHM) with a 0.1 m spatial resolution. To determine the CHM, some researchers compared the DSM at different vegetation growing season times with the initial DTM for the first measurement [39]. The CHM was defined as the difference between the DSM and DTM, which is expressed as follows: To obtain a better match between the LiDAR-derived parameters and quadrat measurements, the side length of the window was set to 1.8 m in the interactive data language (IDL), which was slightly larger than the side length of the quadrat. The minimum canopy height derived from the LiDAR data corresponds to the quadrat with a centre coordinate position (i,j) [28] and is expressed as follows: where h u, v indicates the LiDAR-derived minimum canopy height image pixel value u, v , and CHM i, j is derived from the LiDAR image pixel value at centre coordinates i, j of the CHM. Using the same method, , and , can be obtained as follows: The FVC corresponding to the LiDAR data window with centre coordinates of (i,j) was calculated as the ratio of the LiDAR-derived canopy height pixel value (≥2 cm) to all CHM(i,j) pixel values and is expressed as follows [28]: where , indicates the number of canopy returns (we set the threshold to ≥2 cm) within , , and indicates the number of all returns within , . The grassland model fitting degree R 2 represents the correlation between the independent and dependent variables. The RMSE reflects the regression model error between the calculated and predicted values. A smaller RMSE value corresponds to a more precise regression model. The RMSE of the regression model was calculated from the 64 training samples and 32 validation samples. The R 2 and RMSE values can be expressed as follows: where n represents all validation samples, y represents the ith ground-measured sample value, and y represents the ith predicted sample value.
In multiple linear regression analysis, the results do not improve when additional independent variables are introduced because some values have a certain approximately linear or collinear relationship in remote-sensing information, resulting in strong correlations between multiple variables and easily causing data redundancy resulting in complex model calculations. To avoid the adverse effects of collinearity and correlation among multiple independent variables, this paper used Pearson correlation to filter the independent variables according to the needs of the data.
We used SPSS 19.0 software (IBM SPSS Statistics for Windows, Version 19.0. Armonk, NY: IBM Corp, USA), and Pearson correlation testing of the LiDAR-derived minimum, mean and maximum canopy heights and FVC of 64 training samples was carried out.
We validated the standard deviation (σ) of the LiDAR-derived canopy height in the W0 plot to analyze the influence of different flight heights ranging from 40 m to 110 m.
The symbol σ represents the standard deviation of the LiDAR-derived canopy height data and is defined as follows: where represents LiDAR-derived canopy height pixel values, indicates the average LiDAR-derived canopy height value at the W0 plot, and is the number of all sampling points at the W0 plot.

Ground-Matured Data Analysis
We used box plots to analyze the ground-measured data. Figure 4 shows the 96 ground-measurement quadrat data points. Ground-measured vegetation parameters showed a decreasing trend with increasing grazing intensity. Under the G0 grazing intensity, the ground-measured canopy height reached a maximum of 56.7 cm, the median canopy height was 45.61 cm, and the outlier value was 23.1 cm, while under the G5 grazing intensity, the outlier value was 6 cm, the median canopy height was 3.95 cm, and the minimum of the ground-measured canopy height was 2.5 cm. Under the G2 grazing intensity, the canopy height had a large box range. This conclusion indicates that the samples were scattered with large fluctuations ranging from 11.33 cm to 32.8 cm in canopy height. The maximum ground-measured FVC was 87.3%, and a large fluctuation ranged from 63.8% to 87.3% of FVC in the G3 subplot. The minimum ground-measured FVC was 48.7% in the G5 subplot. The maximum AGB was 443.4 g/m 2 under the G0 grazing intensity. The median AGB in the G5 subplot was 96.18 g/m 2 , with a small fluctuation range from 62.6 g/m 2 to 166.16 g/m 2 of AGB.

LiDAR-Derived Vegetation Parameters
Using the LiDAR-derived results specified in Figure 5a,b, spatial distribution maps of the canopy height and FVC at the plot level using pixel data under different grazing intensities in the grazing platform were generated. Figure 5a,b shows that the LiDARderived canopy height and FVC at each subplot followed the expected order, with low grazing intensity observed for high vegetation. At the quantitative level, reduced grassland parameter values were associated with increased grazing intensity in each plot. Figure. 5c,d shows the LiDAR-derived grassland structures under different grazing intensities at the plot level. Under the G0 grazing intensity, the LiDAR-derived canopy height was 26.25 cm and the LiDAR-derived FVC was 78.69%. Under the G5 grazing intensity, the LiDAR-derived canopy height was 2.76 cm and the LiDAR-derived FVC was 62.27%. Additionally, the ground-measured canopy height and FVC based on 96 quadrats exhibited decreasing trends with increasing grazing intensity. The LiDAR-derived grassland parameter trends in all quadrats were associated with the ground-measured data.

Relationship between AGB and Canopy Height/FVC
We analyzed all ground-measured grassland parameters for each quadrat, as displayed in Figure 6. We analyzed the relationship between the ground-measured canopy height, FVC and AGB. The correlation of all ground-measured data was statistically significant (p < 0001). Analysis of the ground-measured mean canopy height indicated a correlation with the AGB (R 2 = 0.59, RMSE = 60.39 g/m 2 ). The ground-measured maximum canopy height was correlated with the AGB (R 2 = 0.53, RMSE = 57.33 g/m 2 ). The minimum ground-measured canopy height showed the lowest correlation (R 2 = 0.48, RMSE = 60.46 g/m 2 ). The ground-measured FVC showed a correlation with R 2 = 0.50 and RMSE = 63.93 g/m 2 . Thus, the results showed that grassland AGB can be better estimated using the ground-measured mean and maximum canopy heights and FVC. The multiple linear regression results were also promising. Additionally, when the ground-measured maximum canopy height and FVC were combined, the AGB regression was slightly improved from R 2 = 0.53 to R 2 = 0.54 and RMSE = 57.33 g/m 2 to RMSE = 55.79 g/m 2 . The results of the ground-measured analysis are similar to previous relevant conclusions [5,28]. This conclusion shows that there is a high correlation between canopy height and FVC, and the fitting degree estimating AGB from LiDAR data was limited because the ground-measured vegetation parameter data accuracy was higher than that of LiDAR-derived data. This conclusion also shows a weaker correlation between groundmeasured AGB and LiDAR-derived vegetation parameters.
As shown in Figure 7, all LiDAR-derived canopy heights and FVC data points from the 64 training samples exhibited statistically significant correlations (p < 0.001) with the ground-measured parameters. However, the LiDAR data always underestimated the ground-measured canopy height and FVC. The LiDAR-derived canopy heights were approximately 1:3.26-1:1.13 as large as the ground-measured canopy heights, and the LiDAR-derived FVC was approximately 1:1.3 as large as the ground-measured FVC. Further research showed the fitting effects of the LiDAR-derived mean canopy height (R 2 = 0.81, RMSE = 4.61 cm) and maximum canopy height (R 2 = 0.92, RMSE = 3.61 cm). In terms of the FVC, the fitting effect was mainly between 75% and 85% (R 2 = 0.51, RMSE = 4.89%), indicating that in most cases, the LiDAR-derived parameters underestimated the groundmeasured vegetation parameters. Because of this correlation, the ground-measured data could be used to correct the LiDAR-derived canopy heights and FVC. Due to the LiDAR signal penetrating the canopy, overlapping leaves appeared as a cluster in the image. If the vegetation was sparse and open, some LiDAR signals could pass through leaves and bounce back to the LiDAR. These openings appeared as holes in the image [40]. Analysis of LiDAR and ground-measured data revealed that the LiDARderived canopy height data were lower than the ground-measured data in the study area. There was a strong correlation between LiDAR and ground-measured data. In grassland ecosystems, ground-measured data can be used for calibration to improve the estimation accuracy of ground-measured canopy height and FVC. According to the analysis in Table  2, the LiDAR-derived vegetation parameters had significantly high levels of correlation (r > 0.6), and the lowest correlation coefficient (r = 0.62) was found between the LiDARderived mean canopy height and FVC. Therefore, we did not need to combine all LiDARderived canopy heights to calculate the grassland AGB.

LiDAR-Derived AGB
For the regression model analysis, we compared the LiDAR-derived canopy height and FVC to those of the grassland AGB model, and the results showed a statistically significant difference (p < 0.001) from the grassland AGB. A simple power-law regression model to predict the grassland AGB is suggested in Table 3, and the LiDAR-derived maximum canopy height and FVC were identified as prediction factors (R 2 = 0.54, RMSE = 64.76 g/m 2 ). Previous studies have established the power-law regression model as one of the most appropriate models. Because each grassland vegetation parameter can establish a relative growth model and exhibit heteroscedasticity phenomena in biomass models, it is better to build a power function model of independent variables or combined variables [41,42]. The regression model equation is as follows: However, the AGB regression model resulted in almost no significant improvement in the grassland AGB model with only the LiDAR-derived canopy height. As shown in Table 3, the fitting degree of the LiDAR-derived maximum canopy height was better than that of the LiDAR-derived mean canopy height (R 2 = 0.41-0.53, RMSE = 71.98-65.26 g/m 2 ). After the FVC factor was added, the fitting effect of the power function model remained basically unchanged, and the RMSE generally declined, showing that adding the FVC effectively improved the model accuracy. Overall, the AGB model could be more effectively established by combining the power-law models of the canopy height and the FVC derived from the LiDAR point cloud data.
Using the results in Figure 5 and Table 3, the different LiDAR-derived regression models were evaluated using spatial distribution maps of the AGB with a 0.1 m spatial resolution. As shown in Figure 8a-c, low grazing intensity was observed for high AGB, while high grazing intensity was observed for low AGB. Figure 8d-f shows the LiDAR-derived AGB under different grazing intensities at the plot level. The LiDAR-derived AGB showed a decreasing trend due to high grazing conditions. The LiDAR-derived grassland parameter trends in all quadrats were associated with the ground-measured data. However, overestimated or underestimated deviations may have been caused by topographic errors in the Riegl VUX-1 sensor. Overall, the Riegl VUX-1 sensor discrete-return LiDAR could distinguish grassland vegetation structures under different grazing conditions.

Analysis of the Flight Altitude
We also analyzed different flight altitudes (40-110 m) in the study area at the W0 subplot. In Figure 9, the analysis of flight height results shows that the LiDAR-derived canopy height exhibits a decreasing trend. The LiDAR-derived canopy height decreased from 9.19 cm to 8.17 cm, and the standard deviation of the LiDAR-derived canopy height decreased from 3.31 cm to 2.35 cm. However, the LiDAR-derived FVC showed an increasing trend with increasing flight height, and the LiDAR-derived FVC increased from 30.05% to 68.46%. The slowly decreasing LiDAR-derived canopy height results also showed that Riegl VUX-1 has high accuracy at a certain flight height. Therefore, more accurate data must be obtained, and we need to set a low flight altitude in the study area.

Discussion
We found that when estimating AGB based on the relevant vegetation index of ordinary optical images, saturation phenomenon easily occurs in the high vegetation coverage areas because ordinary optical images have difficulty detecting canopy height [43]. Based on the backscatter SAR sensor used to detect vegetation parameters, the change in signal wavelength was relatively weak, so it was difficult to accurately detect grassland canopy structure parameters and AGB [44,45]. Due to the characteristics of LiDAR data, three-dimensional canopy structure information can be directly obtained. The emergence of GEDI, ICESAT-2 and NISAR multisource sensor data fusion methods will contribute to the evaluation of AGB maps. Compared with other prediction maps, these maps have higher accuracy, more spatial coverage cover and higher spatial resolution [46,47].
All results show that UAV LiDAR data can estimate grassland AGB under different grazing intensities. The results also suggest that LiDAR-derived vegetation structures were well correlated with the ground-measured relative vegetation parameters ( Figure 6). The LiDAR-derived maximum canopy height combined with the FVC showed the strongest correlation with the AGB and the highest accuracy (R 2 = 0.54, RMSE = 64.76 g/m 2 ) among the examined canopy structure indices. The analysis of LiDAR-derived parameters also showed that the improvement in the R 2 values was limited after adding the LiDARderived FVC values, as shown in Table 3. The RMSE values decreased to 61.48-71.23 g/m 2 when the LiDAR-derived canopy height and FVC were combined to calculate the grassland AGB. In recent years, some researchers have used LiDAR and optical imagebased structure from motion (SfM) alternative 3D point cloud methods to improve model accuracy [48,49]. If we increase the amount of ground-measured sample data, we can use random forest and SfM algorithms to build a more advanced AGB model [50]. The accuracy and fitting effect of the model require further exploration.
We analyzed all LiDAR-derived vegetation parameters and found that all parameters showed strong correlations, as shown in Table 3 (r > 0.6). Other studies focusing on different vegetation species, including conifer/softwood or deciduous/hardwood forests, have also found similar relationships and conclusions [51][52][53]. Such a similar finding can be attributed to the capability of UAV LiDAR sensors to penetrate the grassland canopy at different flight heights. The main difference is underestimation under different UAV LiDAR densities and wavelengths [54,55]. Another influence on LiDAR-derived vegetation structures is the difference in flight heights [28]. Figure 9 shows that different flight altitudes influence the LiDAR-derived grassland parameters from 40 to 110 m. Different flight heights showed no major effect on the derived canopy height. These results are consistent with previous conclusions [56]. The slowly decreasing LiDARderived canopy height results also showed that Riegl VUX-1 exhibited a high accuracy for flight heights of 40 to 110 m. Therefore, more accurate data are needed at low flight altitudes (<40 m) in the study area. Further research is needed using the optimal flight height of 90-100 m for the Hulun Buir grassland with an FVC of ~65%. This result is similar to the results of ground-measured FVC in Hulun Buir grassland. This conclusion can be used to set the flight altitude of UAVs at 90-100 m to accurately and quickly calculate the Hulun Buir grassland FVC.
Different grazing intensities were also analyzed between LiDAR-derived vegetation structures and ground-measured data. Figure 6 shows that the LiDAR-derived canopy height and FVC distribution maps and values were consistent with the grazing intensities shown in Figure 2. Low grazing intensity was observed for high vegetation, and high grazing intensity was observed for low vegetation. The discrete-return Riegl VUX-1 is an advanced three-dimensional sensor that can estimate distribution maps under different grazing intensities. Previous scientists have been limited by severe saturation problems in optical images in grasslands with high vegetation cover when detecting canopy height [2,57]. We found that UAV LiDAR has the advantage of addressing optical sensors, which typically suffer from saturation problems [34,58].
Different UAV LiDAR sensors have varying estimation precisions for the canopy height, grassland AGB and standard deviations. For example, the UAV Velodyne HDL-32E LiDAR sensor has the potential to provide highly accurate grassland vegetation parameters on a large scale [5,59]. Wang et al. [28] used Velodyne's HDL-32E LiDAR sampling densities of approximately 26 points/m 2 at a flight height of 40 m. The LiDARderived mean canopy height obtained from this sensor exhibited a strong linear correlation (R 2 = 0.583, RMSE = 4.9 cm) with the ground-measured mean height. In this paper, the Riegl VUX-1 sensor received point clouds at a higher density (102 points/m 2 ) at a flight height of 119 m. We found that the accuracy of Riegl VUX-1 was better than that of Velodyne's HDL-32E in estimating the maximum canopy height and vegetation AGB. The LiDAR-derived maximum canopy height showed a strong correlation (R 2 = 0.92, RMSE = 3.61 cm) with the ground-measured maximum canopy height. Therefore, when the accuracy requirements are higher, it is recommended to use the Riegl VUX-1 sensor. The reason for the Riegl VUX-1 sensor having higher accuracy is the use of a single-line laser beam with a positioning and orientation system. Different raster resolutions have a great influence on segmentation vegetation accuracy. The results show that the CHM raster resolution was 0.1 m × 0.1 m, and the LiDAR-derived maximum canopy height combined with the FVC was best for estimating the AGB using a power-law regression model (R 2 = 0.54, RMSE = 64.76 g/m 2 ). Similar to the results of Koch et al. [60] and Pu et al. [61]. When the raster resolution increased, image pixels were divided into several parts, and the p value increased. In addition, the segmentation distance of the LiDAR point cloud also affects the p value [27]. Next, we will use more advanced LiDAR sensors to further improve the LiDAR image resolution and obtain a more accurate LiDAR-derived grassland AGB perdition map.
The Riegl VUX-1 sensor provides three-dimensional vertical structures due to its single-line laser beam, which can effectively penetrate vegetation. Therefore, LiDAR data can effectively represent grassland vegetation parameters. According to the results of this study, with increasing grazing intensity, the LiDAR-derived canopy height, FVC and AGB decreased gradually. This important conclusion can be effectively extended to the entire Inner Mongolia grassland grazing area and provide a theoretical reference for sustainable grazing of the entire Inner Mongolia steppe grassland.
UAV LiDAR has been widely used to extract vegetation parameters and detect AGB and other characteristics [62][63][64]. The Riegl VUX-1 sensor has great potential to provide high-precision grassland vegetation parameters over large areas. However, the FVC is not accurate; we should combine LiDAR with high-precision remote sensing images to retrieve the FVC, increase the number of ground-measured samples, and use the random forest method to improve the AGB modelling accuracy in the future. LiDAR data contain three-dimensional spatial information and have a high sampling frequency, although the point cloud data lack the spectral information and texture information of the target objects. Therefore, our future research will aim to combine LiDAR data and multispectral data to determine a grassland AGB model under different grazing intensities.

Conclusions
This study showed different regression methods used to estimate grassland AGB under different grazing intensities from UAV LiDAR data in the Hulun Buir grazing platform in Northeast China. These conclusions can be extended to the whole Inner Mongolian grassland ecosystem. Under different grazing intensities, LiDAR data can effectively obtain the prediction map of canopy height, FVC and AGB. The primary conclusions are as follows. (1) The Riegl VUX-1 sensor could effectively estimate grassland vegetation parameters, such as canopy height and FVC. Among the LiDAR-derived grassland vegetation structures, the LiDAR-derived maximum canopy height combined with the FVC was best for estimating the AGB using a power-law regression model (R 2 = 0.54, RMSE = 64.76 g/m 2 ). (2) The Riegl VUX-1 sensor underestimated the groundmeasured grassland canopy heights; the LiDAR-derived canopy heights were approximately 1:3.26-1:1.13 as high as the ground-measured canopy heights. Therefore, it is necessary to calibrate the LiDAR-derived canopy heights using ground-measured data. (3) Different flight heights had a major effect on the LiDAR-derived FVC. Therefore, it is necessary to use the optimal flight height of 90-100 m for the Hulun Buir grassland with an FVC of ~65%.