- 1Institute of Fragile Eco-environment, School of Geographic Science, Nantong University, Nantong, China
- 2State Key Laboratory of Grassland Agro-ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, China
- 3School of Agriculture and Environment, College of Science, Massey University, Palmerston North, New Zealand
- 4Forestry Station of Huangnan Prefecture of Qinghai Province, Tongren, China
Grassland health assessment (GHA) is a bridge of study and management of grassland ecosystem. However, there is no standardized quantitative indicators and long-term monitor methods for GHA at a large scale, which may hinder theoretical study and practical application of GHA. In this study, along with previous concept and practices (i.e., CVOR, the integrated indexes of condition, vigor, organization and resilience), we proposed an assessment system based on the indicators monitored by unmanned aerial vehicles (UAVs)-UAVCVOR, and tested the feasibility of UAVCVOR at typical household pastures on the Qinghai-Tibetan Plateau, China. Our findings show that: (1) the key indicators of GHA could be measured directly or represented by the relative counterpart indicators that monitored by UAVs, (2) there was a significantly linear relationship between CVOR estimated by field- and UAV-based data, and (3) the CVOR decreased along with the increasing grazing intensity nonlinearly, and there are similar tendencies of CVOR that estimated by the two methods. These findings suggest that UAVs is suitable for GHA efficiently and correctly, which will be useful for the protection and sustainable management of grasslands.
1 Introduction
Grassland covers about one third of the world’s land area, it plays an irreplaceable role in maintaining ecological functions (e.g., carbon storage and water conservation), suppling forage for livestock and wild animals (O’Mara, 2012; Wang et al., 2019). However, grasslands are degraded or degrading continually because of anthropogenic activities (especially grazing) and climate changes (Wang et al., 2019; Wei et al., 2022a). Therefore, it is urgent to estimate grassland health status accurately and efficiently, and furthermore, adapt reclamation activities to make sure the sustainability of the grassland ecosystem.
Grassland ecosystem is a complex system, and it is necessary to establish an integrated grassland health assessment system (GHA, Xu and Guo, 2015). Generally, GHA patterns could be categorized as qualitative or quantitative in nature (Lepak et al., 2022). Qualitative methods estimate or judge grassland conditions based on systematic observations. To date, the most widely adopted qualitative assessment approach is “Interpreting Indicators of Rangeland Health” (IIRH) (Pyke et al., 2002). It estimates three ecosystem attributes, i.e., soil/site stability, biotic integrity, and hydrologic function, and a suite of qualitative indicators collected by interdisciplinary specialists. Despite it was increasing used in USA and Australia, few researchers have applied the IIRH to assess grassland ecosystems health status across large spatial extents (de Soyza et al., 2000a), and a lack of specialists and comparability resulting in inconsistent standards is the major cause limiting the use of IIRH. Quantitative methods were established on numerical indicators, and scientists are always wrestled with how to characterize and keep track of grassland changes using quantitative monitor (Hou et al., 2004; Lepak et al., 2022). The widely adopted quantitative assessment approach is mainly derived from the concept of ecosystem health, e.g., the CVOR integrated index (i.e., condition, vigor, organization and resilience, hereafter COVR) which is proposed based on the VOR index (Hou et al., 2004). CVOR was widely used to assess GHA of various types of grassland in Inner Mongolia (Wang et al., 2008), Ningxia (Yu et al., 2018), and Xijiang Province (Lu et al., 2017). Usually, the quantitative assessing frameworks are developed by process-based studies (Pyke et al., 2002), however, the available resources are often limited and hard to support the work, except for a small number of grassland ecological sites (Miller, 2008), and only focused on the field-based indicators in majority of previous work, which featured with inefficient and labor-intensive (Whitford et al., 1998; de Soyza et al., 2000b; Pyke et al., 2002).
The need for affordable and effective methods across expansive grassland led to lots of efforts focused on the development of indicators that could be reliably detected (de Soyza et al., 2000a; Ludwig et al., 2007; Sun et al., 2018). To date, however, the widely recognized assessment systems are lack, and how to select effective evaluation indexes and establish appropriate evaluation methods with new tools is still largely unknown (Wang et al., 2018). In the last decade, along with unmanned aerial vehicles (UAVs) was applied in studies of species composition (Sun et al., 2018; Sun et al., 2022), vegetation coverage (Chen et al., 2016) and above-ground biomass (Zhang et al., 2018), new opportunities are emerging to enable novel GHA systems that are suitable for multiple scale and featured with high comparability and efficiency.
The Qinghai-Tibetan Plateau (QTP) is the highest plateau in the world, and is often referred to as the Third pole (Yao et al., 2022). Alpine grassland is the most dominant vegetation type, and grazing is the most extensive management mode on the QTP (Harris, 2010). Reasonable grazing management is regarded as the most effective way to utilize and protect alpine grassland (Sun et al., 2020), and GHA is the first step in practical application. Actually, GHA is gradually recognized as a vital component for reasonable management of grassland ecosystems (Herrick et al., 2006; Miller, 2008). In this study, along with previous concept and practices (mainly the CVOR), we proposed an assessment method based on the indicators collected by unmanned aerial vehicles (UAVs) - UAVCVOR, and tested the feasibility of UAVCVOR at the typical household pastures on the QTP. Our specific objectives were to: (1) make sure whether UAVCVOR is suitable for GHA, especially for alpine meadow, (2) compare the differences of CVOR that were estimated by UAV- and field-based sampling methods, and (3) reveal the GHA changes along with the gradients of grazing intensity. Results of this study could have significant implications on reasonable monitor and management of grassland and socio-ecological sustainable development.
2 Materials and methods
2.1 Study sites
In 2017, This study was conducted at Azi Research Station in Gansu Province, China (101°52′07.9″E, 33°24′24.1″N; 3547m a.s.l.) (Figure 1). The study area is located in a humid region, mean annual precipitation > 600 mm and mean annual air temperature is approximately 1.1°C (Sun et al., 2015). The soil is alpine meadow soil (Gao and Li, 1995). The plant community is dominated by Poaceae and Cyperaceae, and meanwhile, various dicotyledonous species are also common, e.g., the Polygonaceae, Saxifragaceae, Ranunculaceae, Asteraceae (Ma et al., 2010).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1150859/fpls-14-1150859-HTML-r1/image_m/fpls-14-1150859-g001.jpg)
Figure 1 Study area in the Eastern Qinghai-Tibetan Plateau (the triangle in the insert). Gray areas indicate sampling household pastures, lines indicate pasture borders, black circles represent campsites and back square represent the control plot.
Three typical household pastures (the aera is 48.53–113.64 ha) that are primarily used for grazing in warm season were selected, the pastures are gentle topography (slope< 5°) in case the slope and aspect affect the plant and herbivores (Figure 1). As conventional grazing pattern, yaks are free grazing during daytime and were penned at night. The pastures we selected have been continually grazed for more than 30 years in this way, and it generates a radial gradient of grazing intensity from the campsites (Chillo et al., 2015; Sun et al., 2018). Therefore, it provides suitable research areas to study the relationships between GHA and grazing intensity, and meanwhile, test the feasibility of the proposed UAVCVOR.
Nine sample points (the definition is the centers of Belt and Quadrat routes in this study, Figure 2) were set up randomly between campsites and the margin of the pastures. The distances of the sample points from the corresponding campsite were measured based on georeferenced orthomosaic of the pastures, which were generated using overlapping aerial photographs acquired from the Mosaic flight mode (e.g., Figures 2D, E). It was recognized that grazing intensities are the proxy parameters of inverse distance from the water point or campsite (Wesuls et al., 2013; Chillo et al., 2015; Sun et al., 2018). Three fencing paddocks (used as control treatments of multi-stocking rates and grazing system experiments) set up since 2010 were selected as the control plots (Sun et al., 2015; Sun et al., 2018). All the monitor missions were conducted with permissions from both the landowner and airspace authority.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1150859/fpls-14-1150859-HTML-r1/image_m/fpls-14-1150859-g002.jpg)
Figure 2 Indicators of grassland health assessment sampling and monitor by unmanned aerial vehicle and field-based verification (Using the No.1 pasture as an example). (A) the real-time control interface for Rectangle mode, (B) Belt mode, (C) Quadrat mode, (D) the georeferenced orthomosaic of the studied pasture (dotted black line indicates borders of the pasture, and yellow circle represents campsite area), (E) Mosaic mode, (F) species composition and biomass field-based verification, and (G) soil sampling.
2.2 Sampling in the field and collecting data
Based on the empirical and experimental studies (Hou et al., 2004; Shan et al., 2012; Yu et al., 2018), the VOR integrated index was used to describe biological characteristics of grassland ecosystems, which included vigor, organization and resilience and were described by aboveground biomass, richness, and the ratio of the species that indicate restoration and degraded degree (RRD) (Shan et al., 2012; Yu et al., 2018). To describe the integrated environmental factor, the soil organic carbon (SOC) was sampled to establish a new practical integrated index - CVOR (Hou et al., 2004). In this study, we followed the previous work and calculated the CVOR using the same methods at each sample point (Hou et al., 2004; Shan et al., 2012; Yu et al., 2018). To explore relationships between the CVOR values that were estimated by field-and UAV-based methods, three quadrats were sampled within the Quadrat routes randomly (detailed as below). The aboveground biomass and species composition were measured and average values were used in statistical analyses.
UAV-based sampling was also conducted. At each sample point, a Belt (the ground sampling distance was about 0.09 cm and covered 2.6 m × 3.5 m on ground), Quadrat and Rectangle route were set up using the FragMAP system (Yi, 2017) (Figure 2). Given that fractional vegetation cover (FVC) is one of the key controlling factors in transpiration, photosynthesis and other terrestrial processes, and meanwhile, it is also an important variable that could be used to describe vegetation quality and reflect changes of ecosystem (Graetz et al., 1988; Hirano et al., 2004). In this study, we selected FVC as the Condition index, and Rectangle routes were set to collect the FVC of sampling points. In brief, 12 way points were set within an area of 20,000 m2 (i.e., 100 m × 200 m), the height was set at 20 m and the UAV took one aerial photograph vertically at each way point (to avoiding disturbance factor, e.g., the aerial photographs taken on the rooftop, road, campsite or river were excluded from this study, and finally 3–5 aerial photographs around each sampling point were selected to estimate the FVC (Figure S1). The ground resolution was about 1 cm and the FVC could be calculated by the threshold method, the flight mode and data collection were introduced in Chen et al. (2016) and Yi (2017).
The biomass (regarded as Vigor index) was estimated by Quadrat mode. Briefly, 28 aerial photographs were taken vertically downward and 8 azimuth angles (2 m height), then the aerial photographs were matched and extracted point cloud information to estimate the biomass (Zhang et al., 2018). The richness (regarded as the Organization index) and RRD (regarded as the Resilience index) were estimated based on the aerial photographs that were collected by Belt mode, i.e., all species were identified visually and recorded within each aerial photograph, and the indexes were then calculated for each sample point (there are 74 plant species and 3 ones were not identified on the aerial photographs, see Sun et al., 2018 for details). DJI Mavic Pro UAV was used to take aerial photographs (DJI Innovation Company, China), which equipped with a built-in 20-megapixel RGB camera, and featured with terrain following function.
2.3 Calculation of CVOR
The integrated index CVOR (both field- and UAV-based methods, field CVOR and UAVCVOR) were calculated followed the previous studies (Hou et al., 2004; Shan et al., 2012; Yu et al., 2018):
where C is the integrated environmental index, C ∈ [0, 1] and it will be valued as 1 if C > 1; VOR is the integrated biological index, VOR = (WV × V + WO × O + WR × R), where WV, WO and WR are the weighting coefficients that indicate the importance of each single index, avoiding the errors that could be caused by spatial heterogeneity and temporal variability. Given that the sampling background was consistent in this study, the values were calculated as WV = WO = WR= 1/3, VOR ∈ [0, 1] and it will be valued as 1 if VOR > 1; CVOR ∈ [0, 1] and it will be valued as 1 if CVOR > 1.
For field CVOR, the soil organic carbon concentration (SOC) was detected by the K2Cr2O7-H2SO4 oxidation method (Nelson and Sommers, 1996), and the Condition index (CSOC) for SOC was calculated as:
where SOCxand SOCck are the concentrations of soil sampled in the sampling points and control plot, respectively.
The aboveground biomass was clipped from the quadrats and oven-dried to constant weight at 65 °C, then the Vigor index (VfB) for aboveground biomass was calculated as:
where Bxand Bck are the aboveground biomass sampled in the sampling points and control plot, respectively.
Different plant functional types’ relative abundance could be regarded as an important factor that affects ecosystem responses to herbivores grazing (de Soyza et al., 2000a; Lunt et al., 2007). The richness was measured by counting all the species within the quadrats, then the Organization index (OfR) for richness was calculated as:
where Rx and Rck are the richness within quadrats sampled in the sampling points and control plot, respectively.
The RRD was measured by counting the functional groups within the quadrats, then the Resilience index (RfD) for RRD was calculated as:
where Pg and Pd are the proportion of restored and degraded species in the sampling points, and Pgck and Pdck are the proportion of restoration and degraded species in control plot, respectively.
For UAVCVOR, the FVC was retrieved by a threshold method which is conducted a Java-based software (Chen et al., 2016; Yi, 2017), and the Condition index (CFVC) for the FVC was calculated as:
where FVCxand FVCck are the vegetation coverage in the sampling points and control plot, respectively.
The biomass was measured based on the point cloud from the aerial photographs, and the Vigor index (VUB) for the UAV-based biomass was calculated as:
where Bxand Bck are the aboveground biomass sampled in the sampling points and control plot, respectively.
The richness was measured based on species frequencies occurred in the aerial photographs collected along Belt routes, and the Organization index (OUR) for the UAV-based richness was calculated as:
where Rxand Rck are the species frequencies occurred within Belt routes sampled in the sampling points and control plot, respectively.
The RRD was measured by counting the functional groups within the Belt routes, and the Resilience index (RUD) for the UAV-based richness was calculated as:
where Pg and Pd are the proportion of restoration and degraded species in sampling points, and Pgck and Pdck are the proportion of restoration and degraded species in control plot, respectively.
2.4 Data analysis
To test the normality of data, a goodness-of-fit test (based on Shapiro-Wilk test and univariate procedure) was used firstly for all the collected data. For the regression analysis, the R2 and P values were used to compare the performance of sampling methods (i.e., UAVCVOR vs. field CVOR). To select the final regression models which indicated the effects of grazing intensity on GHA, the likelihood ratio tests were used to compare the regression models (n = 27), and select the optimal regression models to reveal the relationships between grazing intensity and GHA. The statistical analyses were performed with R version 4.1.3.
3 Results
3.1 The characteristics of UAVCVOR
The SOC increased significantly with the increase of FVC, but the rate of increase was significantly faster when FVC≥ 0.98 compared to that when FVC< 0.98, showing a segmented function between the SOC and FVC (Figure 3A). After transformation based on the segmented function, the FVC was positively corelated to the transitional indicator (FVCSOC, Figure 3B), reflecting the condition of the sample points. RRD could be extracted from aerial photographs that were collected by FragMAP system and analyzed easily and efficiently, and it was significant corelated to the field-based RRD (Figure 3C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1150859/fpls-14-1150859-HTML-r1/image_m/fpls-14-1150859-g003.jpg)
Figure 3 Relationships among indicators that estimated by field- and UAV-based methods. (A) relationship between vegetation coverage (FVC) and soil organic carbon concentration (SOC) (SOC = 228.29 FVC -137.71, R2 = 0.612, P< 0.001 when FVC< 0.980, and SOC = 3938.6 FVC - 3786.9, R2 = 0.725, P< 0.001 when FVC ≥ 0.980); (B) relationship between the transitional indicator (FVCSOC) and actual measured FVC, (C) relationship between ratio of the species that indicate restoration and degraded degree (RDD) that estimated by field- and UAV-based methods (RRDfield = 1.87 RRDUAV - 0.84, R2 = 0.713, P< 0.001).
Compare to the nearest sample sites from campsites, the farthest sample sites exhibited higher estimated values for aboveground biomass, richness and RDD (Figure 4, P< 0.05). The UAV-based method estimated higher values of the richness both at the nearest and farthest sample points, and the difference was greater at the nearest sample points (Figure 4B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1150859/fpls-14-1150859-HTML-r1/image_m/fpls-14-1150859-g004.jpg)
Figure 4 Comparation of the biological index of the nearest and farthest sample points of the tree household pastures that estimated by field- and UAV-based method. (A) aboveground biomass, (B) richness, (C) ratio of the species that indicate restoration and degraded degree (RRD).
3.2 Comparation of the CVORs that were assessed by field- and UAV-based methods
There were significant linear relationships between CVOR estimated by the field- and UAV-based methods (P< 0.001, Figure 5). The variation range (0.44 – 0.98) of values estimated by UAV-based method was smaller than that (0.15 – 0.93) estimated by the field-based method, especially around areas with higher grazing intensities.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1150859/fpls-14-1150859-HTML-r1/image_m/fpls-14-1150859-g005.jpg)
Figure 5 Relationships between grassland health assessment integrated index that were estimated by the traditional field- and UAV-based methods (fieldCVOR = 1.87 UAVCVOR - 0.84, R2 = 0.713, P< 0.001).
3.3 Relationships between GHA and grazing intensity
CVOR index was initially decreased nonlinearly, and then kept stable with increasing grazing intensity (Figure 6). Although the tendency was same, the UAVCVOR exhibited higher values than that estimated by the field-based method for the same sampling point, and furthermore, the differences increased with the increasing grazing intensity (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1150859/fpls-14-1150859-HTML-r1/image_m/fpls-14-1150859-g006.jpg)
Figure 6 The grassland health assessment integrated indicator (COVR) along the grazing intensity (fieldCVOR = 0.0000268 GI -1.45 + 0.176, R2 = 0.866, P< 0.001; UAVCVOR = 0.0000926 GI -1.22 + 0.474, R2 = 0.884, P< 0.001).
4 Discussion
4.1 UAVCVOR establishment for GHA
It is challenged for grassland managers to make management decisions across large landscapes based on limited time and financial resource (Lepak et al., 2022). GHA is a bridge of grassland study and management (Xu and Guo, 2015). Followed the methodology of Hou et al. (2004); Yu et al. (2018) and etc., we established the integrated index-CVOR based on the indicators represent the environmental conditions and biological characteristics (i.e., C + VOR). Similar to the SOC concentration, the vegetation coverage is regarded as a sensitive indicator of land degradation and desertification (Hirano et al., 2004; Jiapaer and Bao, 2011). In this study, the relationship between FVC and SOC was a segmented function (Figure 5), but the normalized transformation made it possible to indicate the SOC level by FVCSOC. In addition, in this study, the vegetation and soil type were similar, thus the FVC could be an appropriate indicator for the Condition index.
Regarding the biological indexes, it has been found that the UAV-based method could estimate the aboveground biomass and richness of alpine grassland efficiently and correctly (Sun et al., 2018; Zhang et al., 2018). Similarly, the UAV-based method could also be used to estimate the RRD (the significant linear relationship between the values that measured by UAV- and field-based methods, Figure 3C). Therefore, the CVOR estimated by UAV-based method is reliable and efficient.
4.2 Changes of GHA along grazing intensities
Quantitative data describing condition gradients of specific ecosystem will be of utility in assessment, monitoring, and sustainable management of grassland ecosystems (Bosch and Kellner, 1991; Wei et al., 2022b), meanwhile, it is also an efficient study method to the scientists engaging in grassland ecosystem related research activities (Herrick et al., 2006). Specifically, resource inhomogeneous distribution or management mode could result in phenomena of radial gradient of grazing intensity in specific areas such as household pastures, which is regarded as a useful way to understand and predict the relationships between grasslands and herbivores (Chillo et al., 2015; Sun et al., 2018). The superior advantage is that it includes multi-grazing intensities, rather than the only grazing & enclosure or several specific grazing intensities treatments (Lunt et al., 2007), which is helpful to reveal the potential relationships between grazing intensity and the response indictors (Sun et al., 2018; Teague and Kreuter, 2020). In this study, we demonstrate a nonlinear relationship between the CVOR (estimated by field- and UAV-based methods) and grazing intensities (Figure 6). Meanwhile, compared to traditional field-based method, the UAV-based method could sample within larger ranges both for single monitor unit and total monitor areas, which meaning the higher representativeness and lower errors for the specific monitor area (Sun et al., 2018), especially for the species richness (Figure 4B). Furthermore, the higher representativeness improves estimation accuracy (Fuhlendorf and Engle, 2001), especially for the grassland ecosystem featured with heterogeneity. In this study, the CVOR values estimated by the UAV-based method were higher than those measured by the field-based method could be resulted from the higher representativeness (Figures 5 and 6, Sun et al., 2018). Hence, we conclude that the UAV-based method is appropriate for GHA, and it could be an optional way for scientific research and practical guidance, such as provides dynamic optimum grazing intensities at various scales according to the GHA, and fills the gap between botanical characteristics that were studied in plot experimental and natural communities (Genung et al., 2020; Sun et al., 2022).
4.3 Advantages of UAV-based GHA
In general, different biophysical parameters of grasslands require scale-matched monitoring method prior to capture their spatial or temporal variation (Bradley and Millington, 2006; He et al., 2006; He et al., 2007). Thus, it is better if the methodology of GHA could include analyses at multiple spatio-temporal scales. However, up to now, it is still a challenge to find suitable scale (both temporal and spatial scales) thresholds for GHA modeling (Xu and Guo, 2015). In recent years, the UAVs had been used to fill the spatial gap between the field-based sampling and satellite images (Chen et al., 2016; Zhang et al., 2023), which overcomes the scale issue for ecological studies. Similarly, the UAVs controlled by FragMAP system could be used for GHA estimate within a similar phenological period at a large scale. The results would be more precise because of the larger sampling area and higher representativeness, and which could reduce sampling bias in the heterogeneous landscapes (Sun et al., 2018). Besides, the UAVs can reach the areas where are difficult for human walking (e.g., wetland), and data are collected efficiently (Floreano and Wood, 2015). Therefore, the UAV-based method makes it possible to carry out GHA in situ at a large scale, which is especially useful on the QTP featured with high elevation and low oxygen (Sun et al., 2022).
Generally, GHA requires long-term repeated monitoring to realize comparability (Arif et al., 2016; Lepak et al., 2022). The UAV-based method (controlled by FragMAP system) proposed in this study realized flights repeatedly and efficiently at different times with fixed height (Yi, 2017; Sun et al., 2018). Meanwhile, the large numbers of sampling locations that represent general patterns of GHA among ecological sites (de Soyza et al., 2000a). Moreover, the UAV-based noninvasive sample method could keep the sample sites at ideal state which makes sure that the results are more believable and comparable, and it is crucial for the control or reference sites that could be destroyed by sampling repeatedly (Lepak et al., 2022).
4.4 Limitations of UAVCVOR and suggestions to improve UAV-based GHA
We proposed and tested the UAV-based GHA method at typical household pastures on the QTP. The method exhibited more efficiency and representativeness of larger sample aeras (e.g., Sun et al., 2018). However, we do acknowledge the limitations of UAV-based GHA method. First, some species such as the creeping or low-growing plants in UAV photographs could not be identified (Sun et al., 2018). Fortunately, higher resolution of UAVs could allow us to identify detail information, meanwhile, taking additional aerial photographs manually at close distance (e.g., 0.5 m) will make it possible to find detail under the canopy as the airflow made by the UAVs’ propellors (Sun et al., 2018). Second, to date, the visual identification process requires substantial time and professional knowledge, the target identification automatically by increasing the machine learning algorithm could further promote the efficiency of UAV-based method (Lu and He, 2017; Menshchikov et al., 2021).
Though we demonstrated the feasibility of UAVCVOR at the household pastures in this study, the feasibility of the UAV-based integrated indexes to different types of grassland at multiple scales is still necessary. Meanwhile, the specific indicators included in CVOR could make adjustments according to study need, e.g., taking the bioclimatic dataset (download from www.worldclim.org/current or other datasets) as Condition index at regional scale. In addition, with filling special gap between the traditional field-based samples and satellite images by UAV technology, the standardized space-air-ground integrated observation system for GHA could be established, and it will be easier to select the control or reference sites scientifically around the large areas. Finally, we may integrate the single monitor points and thus improve the datasets, which will be beneficial to vertical and horizontal comparison and applied in site-specific management strategies and environmental protection (Miller, 2008; Ye et al., 2011).
It is usually hard to separate the effects caused by grazing and succession in ecosystems featured with the dynamic succession processes. Thus, the long-term fixed points monitor system is urgent to develop to solve the confusing issues. The UAV-based sampling processes of GHA was separated into two parts: cooperative UAV-field sampling and net-cooperation data extraction (Sun et al., 2022). UAV-based sampling is more efficiency which could reduce the field work time, and net-cooperation data analysis could be carried out synergistically in laboratories, which overcome the temporal and spatial limitations of data analysis of GHA (e.g., the species composition, Sun et al., 2018); and more important is that the datasets dased on aerial photographs could be the critical point in integrating the qualitative or quantitative protocols (Lepak et al., 2022). Therefore, the UAVCVOR will contribute to study real-world conditions at a large scale, which could be helpful to separate the effects caused by disturbance and succession, and will contribute to the sustainable development of grassland ecosystem.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
Conceptualization, SY and YS; methodology, YS, QM, and SY; software, YS and YLuo; validation, YY, GD, QM, and YLuo; investigation, YLiao, WJ, WW and YS; data curation, YS, YY and XW; writing-original draft preparation, YS, YLuo, WJ, WW and XW; writing-review and editing, SY and YS; visualization, YLiao, GD, and YY; supervision, SY; project administration, YS. All authors contributed to the article and approved the submitted version.
Funding
This study was jointly supported by the grants from the National Natural Science Foundation of China (42071056, 31901393) and Natural Science Foundation of Jiangsu Province (BK20201439), a grant from State Key Laboratory of Grassland Agro-Ecosystems (Lanzhou University) (SKLGAE201906) and National University’ Innovation and Entrepreneurship Training Program (202320304054Z, 202310304145Y).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1150859/full#supplementary-material
References
Arif, M. S., Gülch, E., Tuhtan, J., Thumser, P., Haas, C. (2016). An investigation of image processing techniques for substrate classification based on dominant grain size using RGB images from UAV. Int. J. Remote Sens. 38, 1–23. doi: 10.1080/01431161.2016.1249309
Bosch, O. J. H., Kellner, K. (1991). The use of a degradation gradient for the ecological interpretation of condition assessments in the western grassland biome of southern Africa. J. Arid Environ. 21, 21–29. doi: 10.1016/S0140-1963(18)30724-9
Bradley, A., Millington, A. (2006). Spatial and temporal issues in determining biomass burning regimes in Bolivia and Peru. Int. J. Remote Sens. 27, 2221–2253. doi: 10.1080/01431160500396550
Chen, J. J., Yi, S. H., Qin, Y., Wang, X. Y. (2016). Improving estimates of fractional vegetation cover based on UAV in alpine grassland on the Qinghai-Tibetan Plateau. Int. J. Remote Sens. 37, 1922–1936. doi: 10.1080/01431161.2016.1165884
Chillo, V., Ojeda, R., Anand, M. (2015). A novel approach to assess livestock management effects on biodiversity of drylands. Ecol. Indic. 50, 69–78. doi: 10.1016/j.ecolind.2014.10.009
de Soyza, A., Whitford, W., Turner, S., Zee, J., Johnson., A. (2000a). Assessing and monitoring the health of western rangeland watersheds. Environ. Monit. Assess. 64, 153–166. doi: 10.1007/978-94-011-4343-1_14
de Soyza, A., Zee, J., Whitford, W., Neale, A., Tallent, N., Herrick, J., et al. (2000b). Indicators of Great Basin rangeland health. J. Arid Environ. 45, 289–304. doi: 10.1006/jare.2000.0645
Floreano, D., Wood, R. J. (2015). Science, technology and the future of small autonomous drones. Nature 521, 460–466. doi: 10.1038/nature14542
Fuhlendorf, S. D., Engle, D. M. (2001). Restoring heterogeneity on rangelands: ecosystem management based on evolutionary grazing patterns. BioScience 51, 625–632. doi: 10.1641/0006-3568(2001)051[0625:RHOREM]2.0.CO;2
Gao, Y. X., Li, M. S. (1995). Soil regionalization of the qinghai-tibetan plateau. Mountain Res. 13, 203–211. doi: 10.16089/j.cnki.1008-2786.1995.04.001
Genung, M. A., Fox, J., Winfree, R. (2020). Species loss drives ecosystem function in experiments, but in nature the importance of species loss depends on dominance. Global Ecol. Biogeogr. 29, 1531–1541. doi: 10.1111/geb.13137
Graetz, R. D., Pech, R. P., Davis, A. W. (1988). The assessment and monitoring of sparsely vegetated rangelands using calibrated landsat data. Int. J. Remote Sens. 9, 1201–1222. doi: 10.1080/01431168808954929
Harris, R. B. (2010). Rangeland degradation on the Qinghai-Tibetan plateau: A review of the evidence of its magnitude and causes. J. Arid Environ. 74, 1–12. doi: 10.1016/j.jaridenv.2009.06.014
He, Y., Guo, X., Si, B. (2007). Detecting grassland spatial variation by a wavelet approach. Int. J. Remote Sens. 28, 1527–1545. doi: 10.1080/01431160600794621
He, Y., Guo, X., Wilmshurst, J., Si, B. (2006). Studying mixed grassland ecosystems II: Optimum pixel size. Can. J. Remote. Sens. 32, 108–115. doi: 10.5589/m06-018
Herrick, J. E., Bestelmeyer, B. T., Archer, S., Tugel, A. J., Brown, J. R. (2006). An integrated framework for science-based arid land management. J. Arid Environ. 65, 319–335. doi: 10.1016/j.jaridenv.2005.09.003
Hirano, Y., Yasuoka, Y., Ichinose, T. (2004). Urban climate simulation by incorporating satellite-derived vegetation cover distribution into a mesoscale meteorological model. Theor. Appl. Climatol. 79, 175–184. doi: 10.1007/s00704-004-0069-0
Hou, F. J., Yu, Y. W., Fu, H., Zhu, Z. Y., Liu, Z. L. (2004). CVOR index for health evaluation of Alashan grazing land. Acta Pratac. Sin. 13, 117–126. doi: 10.3321/j.issn:1004-5759.2004.04.021
Jiapaer, G., Bao, A. (2011). A comparison of methods for estimating fractional vegetation cover in arid regions. Agr. For. Meteorol. 151, 1698–1710. doi: 10.1016/j.agrformet.2011.07.004
Lepak, N., Newingham, B., Kachergis, E., Toledo, D., Moffitt, J. (2022). Where do qualitative assessments fit in an era of increasingly quantitative monitoring? Perspectives from interpreting indicators of rangeland health. Rangelands 44, 39–49. doi: 10.1016/j.rala.2021.07.008
Lu, B., He, Y. (2017). Species classification using unmanned aerial vehicle (UAV)-acquired high spatial resolution imagery in a heterogeneous grassland. ISPRS J. Photogramm. Remote Sens. 128, 73–85. doi: 10.1016/j.isprsjprs.2017.03.011
Lu, J., Hu, Y. K., Yue, P., Li, K. H. (2017). Assessment on the health of alpine steppe in Bayinbuluk based on CVOR index. Arid Zone Res. 34, 862–869. doi: 10.13866/j.azr.2017.04.20
Ludwig, J., Bastin, G., Chewings, V., Eager, R., Liedloff, A. (2007). Leakiness: A new index for monitoring the health of arid and semiarid landscapes using remotely sensed vegetation cover and elevation data. Ecol. Indic. 7, 442–454. doi: 10.1016/j.ecolind.2006.05.001
Lunt, I. D., Eldridge, D. J., Morgan, J. W., Witt, G. B. (2007). A framework to predict the effects of livestock grazing and grazing exclusion on conservation values in natural ecosystems in Australia. Aust. J. Bot. 55, 401–415. doi: 10.1071/BT06178
Ma, M. J., Zhou, X. H., Wang, G., Ma, Z., Du, G. Z. (2010). Seasonal dynamics in alpine meadow seed banks along an altitudinal gradient on the Tibetan Plateau. Plant Soil 336, 291–302. doi: 10.1007/s11104-010-0480-5
Menshchikov, A., Shadrin, D., Prutyanov, V., Lopatkin, D., Sosnin, S., Tsykunov, E., et al. (2021). Real-time detection of hogweed: UAV platform empowered by deep learning. IEEE Trans. Comput. 70, 1175–1188. doi: 10.1109/TC.2021.3059819
Miller, M. (2008). Broad-scale assessment of rangeland health, Grand Staircase-Escalante National Monument, USA. Rangeland Ecol. Manage. 61, 249–262. doi: 10.2111/07-107.1
Nelson, D., Sommers, L. (2008). Total carbon, organic carbon and organic matter, in: Methods of soil analysis. Part 3: chemical methods, edited by: Nelson, D., Sparks, D. L., Page, A. L., Helmke, P. A., Loeppert, R. H., Soltanpour, P. N., Tabatabai, M. A., Johnston, C. T., Sumner, M. E.Soil Science Society of America Book Series, Number 5, Soil Science Society of America, Wisconsin, 1996, 961–1010. doi: 10.2111/07-107.1
O’Mara, F. (2012). The role of grasslands in food security and climate change. Ann. Bot. 110, 1263–1270. doi: 10.1093/aob/mcs209
Pyke, D. A., Herrick, J. E., Shaver, P., Pellant, M. (2002). Rangeland health attributes and indicators for qualitative assessment. J. Range Manage. 55, 584–597. doi: 10.2307/4004002
Shan, G. L., Chen, G., Liu, Z. L., Yan, Z. J., Chu, X. H. (2012). VOR and CVOR index for health evaluation of typical steppe in Inner Mongolia. Acta Agrestia Sin. 20, 401–406. doi: 10.11733/j.issn.1007-0435.2012.03.003
Sun, Y., Angerer, P. J., Hou, F. J. (2015). Effects of grazing systems on herbage mass and liveweight gain of Tibetan sheep in Eastern Qinghai-Tibetan Plateau, China. Rangeland J. 37, 181–190. doi: 10.1071/RJ14062
Sun, J., Liu, M., Fu, B. J., Kemp, D., Zhao, W., Liu, G., et al. (2020). Reconsidering the efficiency of grazing exclusion using fences on the Tibetan Plateau. Sci. Bull. 65, 1405–1414. doi: 10.1016/j.scib.2020.04.035
Sun, Y., Yi, S. H., Hou, F. J. (2018). Unmanned aerial vehicle methods makes species composition monitoring easier in grasslands. Ecol. Indic. 95, 825–830. doi: 10.1016/j.ecolind.2018.08.042
Sun, Y., Yuan, Y. X., Luo, Y. F., Ji, W. X., Bian, Q. Y., Zhu, Z. Q., et al. (2022). An improved method for monitoring multiscale plant species diversity of alpine grassland using UAV: A case study in the Source Region of the Yellow River, China. Front. Plant Sci. 13, 905715. doi: 10.3389/fpls.2022.905715
Teague, R., Kreuter, U. (2020). Managing grazing to restore soil health, ecosystem function, and ecosystem services. Front. Sustain. Food Syst. 4, 534187. doi: 10.3389/fsufs.2020.534187
Wang, L., Delgado-Baquerizo, M., Wang, D. L., Isbell, F., Liu, J., Feng, C., et al. (2019). Diversifying livestock promotes multidiversity and multifunctionality in managed grasslands. Proc. Natl. Acad. Sci. 116, 201807354. doi: 10.1073/pnas.1807354116
Wang, L. X., Liu, Z. L., Liu, H. M., Wang, W., Liang, C. Z., Qiao, J., et al. (2008). Assessment of typical steppe ecosystem health in Inner Mongolia. Acta Ecologica Sin. 28, 544–550. doi: 10.3321/j.issn:1000-0933.2008.02.012
Wang, D. L., Wang, L., Liu, J. S., Zhu, H., Zhong, Z. W. (2018). Grassland ecology in China: Perspectives and challenges. Front. Agric. Sci. Eng. 5, 24–43. doi: 10.15302/J-FASE-2018205
Wei, W. R., An, S. H., Zheng, Q. Y., Qin, M. S., Chen, T. (2022b). Structural changes in vegetation coincident with reseeding Elymus nutans can increase perceived predation risk of plateau pikas (Ochotona curzoniae). Appl. Anim. Behav. Sci. 255, 105745. doi: 10.1016/j.applanim.2022.105745
Wei, W. R., Yao, X., Zhang, Y., Zhen, Q. Y., Qin, M. S., Tang, Z. M., et al. (2022a). Vegetation restoration measures: Increasing plant height suppresses population densities of plateau pikas. Land Degrad Dev. 34, 2201–2213. doi: 10.1002/ldr.4600
Wesuls, D., Pellowski, M., Suchrow, S., Oldeland, J., Jansen, F., Dengler, J. (2013). The grazing fingerprint: Modelling species responses and trait patterns along grazing gradients in semi-arid Namibian rangelands. Ecol. Indic. 27, 61–70. doi: 10.1016/j.ecolind.2012.11.008
Whitford, W., Soyza, D. A., Zee, V. J., Herrick, J. E., Havstad, K. M. (1998). Vegetation, soil, and animal indicators of rangeland health. Environ. Monit. Assess. 51, 179–200. doi: 10.1023/A:1005987219356
Xu, D. D., Guo, X. L. (2015). Some insights on grassland health assessment based on remote sensing. Sensors 15, 3070–3089. doi: 10.3390/s150203070
Yao, T. D., Bolch, T., Chen, D. L., Gao, J., Immerzeel, W., Piao, S. L., et al. (2022). The imbalance of the Asian water tower. Nat. Rev. Earth Environ. 3, 618–632. doi: 10.1038/s43017-022-00299-4
Ye, X., Zhou, H. K., Zhao, X. Q., Wen, J., Chen, Z., Duan, J. C. (2011). Review on grassland ecosystem health. Pratac. Sci. 28, 549–560. doi: 10.3969/j.issn.1001-0629.2011.04.005
Yi, S. H. (2017). FragMAP: a tool for long-term and cooperative monitoring and analysis of small-scale habitat fragmentation using an unmanned aerial vehicle. Int. J. Remote Sens. 38, 2686–2697. doi: 10.1080/01431161.2016.1253898
Yu, H. Q., Jiang, Q., Wang, Z. J., He, J. L., Chen, He. (2018). Application of VOR and CVOR index for health assessment of desert steppe in Ningxia-Taking Yanchi County for example. Acta Agrestia Sin. 26, 584–590. doi: 10.11733/j.issn.1007-0435.2018.03.009
Zhang, H. F., Sun, Y., Chang, L., Qin, Y., Chen, J. J., Qin, Y., et al. (2018). Estimation of grassland canopy height and aboveground biomass at the quadrat scale using unmanned aerial vehicle. Remote Sens. 10, 851. doi: 10.3390/rs10060851
Keywords: grassland ecosystem, grassland health monitor, biological index, FragMAP, environmental condition
Citation: Luo Y, Ji W, Wu W, Liao Y, Wei X, Yang Y, Dong G, Ma Q, Yi S and Sun Y (2023) Grassland health assessment based on indicators monitored by UAVs: a case study at a household scale. Front. Plant Sci. 14:1150859. doi: 10.3389/fpls.2023.1150859
Received: 25 January 2023; Accepted: 04 September 2023;
Published: 20 September 2023.
Edited by:
Juanle Wang, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Huakun Zhou, Chinese Academy of Sciences (CAS), ChinaWei Wanrong, China West Normal University, China
Copyright © 2023 Luo, Ji, Wu, Liao, Wei, Yang, Dong, Ma, Yi and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yi Sun, sunyi@ntu.edu.cn