Next Article in Journal
A New AS-PCR Method to Detect CSN201 Allele, Genotyping at Ca-Sensitive Caseins Loci and Milk Traits Association Studies in Autochthonous Lazio Goats
Previous Article in Journal
Comprehensive Analysis of Differentially Expressed CircRNAs in the Ovaries of Low- and High-Fertility Sheep
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial- and Temporal-Trajectory Analysis of the Crested Ibis (Nipponia nippon) by Fusing Multiple Sources of Data

1
School of Geography and Information Engineering, China University of Geosciences, Wuhan 430074, China
2
Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China
3
School of Computer Science, China University of Geosciences, Wuhan 430074, China
4
Key Laboratory of Geological Survey and Evaluation of Ministry of Education, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Animals 2023, 13(2), 237; https://doi.org/10.3390/ani13020237
Submission received: 22 November 2022 / Revised: 6 January 2023 / Accepted: 7 January 2023 / Published: 9 January 2023
(This article belongs to the Section Birds)

Abstract

:

Simple Summary

This paper provides a methodology for analyzing the trajectory of the Crested Ibis (Nipponia nippon): its clustering and mining habitats, trajectory-data supplementation, the analysis of behavioral patterns, and finally, potential habitat mining. It was found that the characteristics of the Crested Ibis habitat in nature reserves include proximity to water sources and cultivated land close to woodlands. Additionally, there is a certain aggregation and sharing of night roosting sites, which are easily disturbed by human activities, and there is a strong correlation between the spatial proximity of the Crested Ibis night. Although human behavior threatens the survival of the Crested Ibis, they still depend on human behavior, mainly for food and water. Moreover, farmers are encouraged to store water in paddy fields in winter and to avoid using harmful substances such as pesticides. Luoshan County provides a protected area for the Crested Ibis; it was analyzed and found to be appropriately located, with an excellent overall habitat-rating. The Chinese Government has also set up corresponding ecological and environmental protection zones, and it is essential to provide solutions for the reintroduction of the Crested Ibis.

Abstract

The Crested Ibis (Nipponia nippon) is an endangered animal with an extremely high ecological, humanistic, and scientific value. However, this species still faces survival challenges, due to rapidly shrinking foraging grounds, the serious interference of human behavior, and increased habitat requirements. Geographical environment is a significant factor affecting Crested Ibis behavior-pattern analysis and habitat protection. The spatial and temporal trajectory contains habitat location and period information, a vital record of the Crested Ibis’ habits, and the basis of all research. Nevertheless, there are only a handful of studies on the missing trajectory data and fusing multiple sources of environmental data-research methods. We studied the spatial and temporal habitat use of the tracked Crested lbis by fusing multiple data-sources in China. This paper adopts the LSTM (long short-term memory) model to supplement the missing trajectory data and perform cluster mining; and a random forest model is used to predict the habitat of the Crested Ibis with high fitting accuracy (R2 = 84.9%). The results show that the Crested Ibis distribution-pattern is characterized by high altitude and proximity to woodland and rivers. Additionally, the habitat dependence on the village implicates human agricultural activities in positively impacting its reproduction. This paper provides a complete method for analyzing Crested Ibis’ spatial and temporal trajectory by fusing multi-source data, which is crucial for protecting the survival and reproduction of the Crested Ibis.

1. Introduction

Whilst not critically endangered, many animals are facing a very high probability of extinction in the wild in the near future [1]. The protection of endangered species constitutes an important part of the ecosystem, and the first task in this regard is habitat protection and restoration [2].
The Crested Ibis (Nipponia nippon), a medium-sized wading bird that was once distributed in East Asia, across China, Japan, Russia, Korea, and South Korea, is one of the most endangered birds in the world, and is classified as a Class I national priority animal [3]. In 2019, it was registered as a vulnerable species (VU) on China’s Red List and as endangered by the International Union for Conservation of Nature (IUCN) [4], due to its extremely high ecological, humanistic, and scientific value [5]. Its population has been on the brink of extinction for approximately a century [6]. However, by 2021, its numbers reached more than 7000 [7]. This is mainly due to the increased efforts to protect the natural environment of the Crested Ibis.
Understanding environmental influence is practical for the reproduction of the Crested Ibis. High-quality habitats can improve the survival ability and ensure the reproduction of Crested Ibis populations [8]. Altitude and wetlands greatly impact the Crested Ibis’ survival [6,9]. Changes in vegetation cover will affect the distribution and reproduction of the Crested Ibis [10]. Human activities should also be considered [11]. The Crested Ibis relies greatly on local farmers [12] and tends to retain its previous nests [13]. Therefore, farmland is vital to maintaining their populations [14,15], as they prefer to forage in shallow rivers [16,17], where food resources have recovered due to a ban on pesticide use [7]. Therefore, based on the spatial and temporal trajectories of the Crested Ibis, it is vital to explore and evaluate their habitat environment to ensure a sufficient and suitable habitat during the conservation and migration process [18].
Some implicit behavioral and activity patterns, potential information on habitat distribution and selection, population dynamics, and group behavior can be discovered through the spatial and temporal trajectory-data analysis of the Crested Ibis [19,20]. Research methods such as comprehensive habitat evaluation [3], maximum pixel method [10], automatic identification method [21] and the Maxent model [22,23], in which Wang Qi et al. [10] explored the vegetation characteristics of Crested Ibis habitats. A manipulative experimental method was adopted to summarize the minimum approach area for tourists, to discover the Crested Ibis’ tolerance to human behavioral disturbance during their wandering period [24]. The Crested Ibis has a certain adaptability and some requirements for living conditions in the wild. LU Shaohui et al. found the breeding factors of the Crested Ibis in the Dongzhai Nature Reserve in Henan, including weather, altitude and other interference factors [25]; Hu et al. evaluated and estimated potential collisions on highways and railways by observing the behavioral trajectories of the Crested Ibis [26]. There are also some studies that discuss the survival status of the Crested Ibis after release, and, based on Crested Ibis activity data, explore the best release strategies for the Crested Ibis [27], reproductive differences in the reintroduction of the Crested Ibis [28], and the Crested Ibis population viability [29].
However, some potential problems need to be addressed. Firstly, the problem of the missing trajectory of the Crested Ibis has not been effectively resolved. Moreover, multiple environmental data sources should be integrated, which can help understand the interactions between individuals and their environment, providing decisions and a basis for the conservation and management of the Crested Ibis. Finally, the Maxent model, developed by Steven J. Phillips, does not include Asia in its training data [30]. Additionally, the Crested Ibis suitability for species endemic to Asia remains to be proven, with overfitting problems and a low transferability [31]. Therefore, supplementing missing data, fusing multi-source environmental data, and finding suitable models to predict potential habitats, are required.
This paper integrates the spatial and temporal trajectory-data of the Crested Ibis with multiple environmental data sources. The study was conducted in Luoshan County in Henan Province, Yang County, Hanzhong City, and Chenggu County in Shaanxi Province in China. To address the problem of uneven data distribution caused by the missing trajectory data, an LSTM (long short-term memory) model was built, based on latitude and longitude location information to supplement missing data and ensure its integrity. The biological patterns of the Crested Ibis were analyzed, and a random forest model suitable for the Crested Ibis was designed and implemented, to predict the potential habitat distribution of the Crested Ibis in Luoshan County, Henan Province, China. Our paper can provide a complete method for an analysis of the trajectory of the Crested Ibis and suggestions for protecting the survival and reproduction of the Crested Ibis.

2. Materials and Methods

The main methodological framework of the study is shown in Figure 1. Firstly, the spatial and temporal trajectory-data of the Crested Ibis were supplemented by the LSTM (long short-term memory) model. Secondly, according to the anniversary activity, the parameter-tuned DBSCAN (density-based spatial clustering of applications with noise) method was combined with the spatial exploration analysis to obtain the results of the analysis of biological and behavioral patterns of the Crested Ibis. Finally, a random forest model was trained to predict the Crested Ibis habitat. The habitat prediction results were normalized into ranks 1~10 for evaluation, and the final habitat-rank evaluation map and the contribution rate of impact factors were calculated to analyze the future potential habitat-distribution-range of the Crested Ibis.

2.1. Study Area and Data

We used BeiDou/GPS-enabled tracking backpacks (HQBG3621L) to collect spatial and temporal trajectory-data. The backpacks can automatically transmit the longitude, latitude, and flight times of the Crested Ibis once per hour. In this study, we collected the trajectory data of 7 individuals of the Crested Ibis, a total of 31,825 trajectory data (Table 1), from July 2014 to December 2021, and the study area included Luoshan County in Henan Province, Yang County, Hanzhong City, and Chenggu County in Shaanxi Province (Figure 2). The Dongzhai National Nature Reserve and Hanzhong National Nature Reserve were also in the study area.
Multi-source environmental data were collected to study the interaction between the Crested Ibis and its environment and to analyze its biological patterns and behavior (Table 2). They were DEM (digital elevation model) data, NDVI (normalized difference vegetation index) vegetation-index data, OSM (open-street map) road data, river-distribution vector data, and GlobeLand30 land-cover data (Table 2). The data set of DEM in Geospatial Data Cloud was ASTER GDEM 30M and the location was Luoshan County, Xinyang City, Henan Province. The data set of NDVI in Geospatial Data Cloud was MOD13Q1 250M and the location was Luoshan County, Xinyang City, Henan Province. The download map sheet in GlobeLand30 was N48_25, N48_30, N48_35, N49_25, N49_30, N49_35, N50_25, N50_30, and N50_35. The location of the vector data was Luoshan County, Xinyang City, Henan Province.

2.2. Trajectory Data Supplement

An LSTM model was constructed to supplement missing trajectories, to complete the data preparation work (Figure 1, left). Clearing the failed trajectory data and the different data sampling-time intervals resulted in missing trajectory data. Therefore, the LSTM complementary-data model was constructed. The latitude and longitude sequence-data was input. In this study, a six-layer sequence network model was built, with only adjacency between layers and no cross-layer connections, consisting of two layers of LSTM and two layers of dropout, and finally, an added dense layer and an activation layer. The input sequence data were entered into the LSTM for calculation in the first step, and the neurons were randomly removed by the dropout layer with a probability of 0.3 in the second step. After repeating the above two steps, the fully connected layer and activation layer output the final predicted latitude and longitude sequences.
The LSTM model adopted the loss function of MSE (mean squared error) (Supplementary Figure S1).
M S E = 1 m i = 1 m y i y i ^ 2
The distance error of the test set also needed to be evaluated, using the haversine method [36], where l a t ^ is the predicted latitude, lat is the true latitude, l o n ^ is the predicted longitude, and lon is the true longitude:
D = 2 r arcsin sin 2 l a t ^ l a t 2 + cos ( l a t ^ ) cos l a t sin 2 l o n ^ l o n 2

2.3. Night-Roosting-Points Identification and Analysis

After the data supplement was completed, it was necessary to identify and analyze the night roosting points of the Crested Ibis (Figure 1, middle). This paper classified the spatial and temporal trajectory-data of the Crested Ibis by period. According to the definitions in the International Crested Ibis Conservation Workshop [5], location points were categorized as night roost points, foraging points, and outing points, according to different periods in different months. The distance between adjacent points greater than 1 km represented different night roosting/foraging/outing points (Table 3) [37]. In this paper, the annual activity of the Crested Ibis was divided into wintering, reproductive and colonial periods. The reproductive period lasts from February to June, the colonial period from July to October, and the wintering period from November to January. Finally, the DBSCAN (density-based spatial-clustering-of-applications-with-noise) clustering algorithm was adopted for use on the Crested Ibis data.
Compared to the K-means method, the DBSCAN [38] method has good adaptability to non-convex clusters, its clustering benefits are unaffected by noise, and the clustering results achieve global optimality [39,40]. The DBSCAN algorithm selects night roosting points and marks noise (Supplementary Figure S2, left). Therefore, compared to the K-means method (Supplementary Figure S2), the clustering algorithm of DBSCAN, which is suitable for the spatial and temporal trajectory-data of the Crested Ibis, was used to identify its habitat and noise points that may be disturbed by human factors, after tuning it according to the characteristics of the Crested Ibis. To evaluate the suitability of the clustering parameters for the Crested Ibis’ spatial and temporal data, the silhouette coefficient proposed by Peter J. Rousseeuw (1986) was used [41]. The closer the silhouette coefficient is to 1 and the noise ratio is smaller, the better the suitability of the clustering parameters for the spatial and temporal data of the Crested Ibis (Supplementary Figure S3). Based on the algorithm’s results, the behavioral patterns and biological patterns of the Crested Ibis were initially analyzed by fusing multiple environmental-data sources.

2.4. Habitat Prediction Based on Random Forest

Finally, the trajectory data were further used for habitat prediction to study and evaluate the habitat status (Figure 1, right). We researched in Luoshan County only. Random forest is an integrated algorithm [42] based on bootstrapping in statistics, evolving the bagging method, which combines CART (the Classification and Regression Tree) and random subspace methods to construct a specific number of unrelated trees. The random forest method can prevent the overfitting problem caused by too much information and too high a correlation between trees. It achieves this by combining multiple decision trees and voting or taking the mean value of the decisions’ results to obtain the model’s final result.
The input training-set sample Dsample = {X1, X2, X3, X4, y}; X1 is the DEM data, X2 is the Euclidean-distance measure for roads, X3 is the Euclidean-distance measure for rivers, X4 is the vegetation index, and y is the number of Crested Ibis points in each grid. Training samples were obtained by repeated multiple sampling with put-back using the bagging method, and this sampling result was used as the training set to generate the decision tree. Using the four input feature-variables, n feature variables (n < 4) were randomly selected for each node in the tree. The best splitting point for the decision tree was determined with the selected n feature-variables. Each decision tree (classifier) was grown to its maximum capacity but not pruned, and the results of all decision trees were finally averaged, to construct a regression model.
The metric used to generate CART trees for this study is MSE [42], which determines the divided features, where D is the sample dataset, a is the input features, Ta is the mean set of the decision tree, y t v is the subset labels, and y t v ^ is the mean value of the subset labels:
M S E D , a = min t T a v , + y t v y t v ^ 2
The smaller the value obtained above, the better the sample-set in the tree fit to find the optimal partitioning attribute.
a b e s t = a r g m i n a M S E D , a
The optimal parameters to construct the model were selected (Supplementary Figure S4, Tables S1 and S2). Additionally, the evaluation metrics of the stochastic forest-habitat prediction model, R2, RMSE, MAE, and explained variance were chosen. The mean squared error RMSE (root-mean-square error) and MAE (mean absolute error) were adopted, to measure the difference between predicted and true values. For the random-forest-regression method, it is not sufficient to only explore the accuracy of the data prediction, but it is also necessary to observe whether the model learns the distribution and pattern of the data, and thus the explained variance and R2 are needed to evaluate the fitting ability of the model.
R 2 = 1 i = 1 n y i ^ y i 2 i = 1 n y ¯ y i 2
M A E = 1 n i = 1 n y i y i ^
R M S E = 1 n i = 1 n y i y i ^ 2
e x p l a i n e d _ v a r i a n c e = 1 V a r y i y i ^ V a r y i
Var denotes the variance of the model.

3. Results

3.1. Trajectory-Data Supplementation

After supplementation, the Crested Ibis trajectory data were mainly located in the central and western regions of Luoshan County (Figure 3).
With further study for additional details, we were able to find missing parts of the trajectory data (Figure 4a) and the supplementation of the data (Figure 4b). In the end, the trajectory data of the Crested Ibis were well supplemented.

3.2. Spatial Patterns of Trajectory Points

The trajectory data were clustered and noise was removed by DBSCAN for further analysis. The distribution of night roosting points for the Crested Ibis in Hantai District, Chenggu County, and Yang County tended to follow the rivers, with cultivated fields and ponds close to the rivers (Figure 5a). There were also a few more dispersed night roosting points in the north and south, and as they were in their wandering periods, they tended to select a wider range of night roosting points, with a tendency to choose areas at higher altitudes with taller and bushier trees.
The overall night-roosting-points range of the Crested Ibis was dense, with areas of night roosting points generally characterized by proximity to paddy fields, ponds, and surrounding woodland, providing both night roosting points and foraging for the Crested Ibis (Figure 5b). During the reproductive period, the choice of night roosting points fluctuated, with the birds roaming to more densely wooded areas.
The rest of the night roosting points in the three activity periods were similar in distribution, with Hantai District, Chenggu County, and Yang County in Shaanxi Province as the representative areas. Among the three activity periods of the Crested Ibis, the colonial period existed in the widest range, with a tendency to spread to higher altitudes on the north and south sides (Figure 5b). The reproductive period was more concentrated in the central parts of Hantai District, Chenggu County, and Yang County and was closer to the rivers than the colonial and wintering periods (Figure 5b). Additionally, there was a strong correlation between the spatial proximity of the Crested Ibis night roosting points and activity periods.
Taking Hantai District, Chenggu County, and Yang County in Shaanxi Province as representative regions, the distance-analysis pattern between the night roosting and foraging points was similar for the remaining two locations (Figure 6). Additionally, the foraging points that were very close together were selected very close together, and concentrated around the night roosting points.
Further studies showed a concentration of distances between 60 m and 595.5 m selected for the two kinds of points, with an interquartile spacing of 535.5 m and even a zero distance between them, indicating that the two kinds of points almost overlapped. The standard deviation of 555.9 m indicated a concentrated distance relationship between them and a weak dispersion of data. The mean value of 411 m further indicated the proximity of the night roosting points to the foraging points. With a maximum distance of 5796 m, this represented only 0.1% of the total data, and was most likely a result of the Crested Ibis finding other foraging points to feed at.

3.3. Potential Habitat Mining

The accuracy evaluation indicators of the random forest model: The R2 was 84.9%, and the explained variance was 85.7% (Table 4), indicating that the random forest model designed in this experiment performed well in fitting the data from the Crested Ibis. The RMSE was 21.0, and the MAE was 12.6 (Table 4), indicating that the difference between the predicted and true values was small, and the prediction model was relatively accurate.
In the southern part of Luoshan County, areas with high DEM (digital elevation model) and high vegetation cover provide foraging grounds for the Crested Ibis. In the southern part of the county, there is a certain amount of woodland cover, which is also more suitable for the Crested Ibis to roost and breed (Figure 7). In Luoshan County, the anthropogenic impact on the Crested Ibis is both positive and negative, as the Crested Ibis does not like to be disturbed by humans, so the habitat level of the Crested Ibis is lower in areas closer to roads but slightly higher in areas closer to rivers, where artificial paddy fields and ponds provide foraging grounds for the Crested Ibis.
In addition to the vegetation cover and elevation data, the location of the habitat with respect to roads and rivers should also be considered when selecting habitats for the Crested Ibis. When selecting future conservation sites for this species, one may prefer areas with a higher vegetation-cover and elevation that are located a certain distance from rivers and roads, to ensure a better foraging area for the Crested Ibis while providing a night roosting environment.
Further quantitative analyses of each of the drivers of Crested Ibis habitat-class evaluation were carried out, with larger values indicating a greater influence of this variable on Crested Ibis habitats (Figure 8). The most influential variable in predicting Crested Ibis habitats based on the random forest model is the vegetation index NDVI (normalized difference vegetation index) (0.359), which reflects the degree of variation in the quality of the region’s ecological environment. Additionally, the NDVI has the greatest influence on the selection of Crested Ibis habitat, indicating a strong correlation between Crested Ibis habitat and the degree of vegetation cover. Next is DEM (0.277), which showed that the Crested Ibis is better adapted to higher altitudes and prefers higher-altitude areas for night roosting points. The remaining two variables are River_Dist (0.206), an assessment of Euclidean-distance to rivers, and Road_Dist (0.145), an assessment of Euclidean-distance to roads (Figure 8). Areas closer to rivers and roads, where there is more water and arable land, often meet the dietary needs of the Crested Ibis.
We also found that a very large part of the wetlands is a highly suitable habitat. In addition, the southern part of Luoshan County, where the natural forest is located, corresponds to a wider range of habitats, with a high suitability in the northern part, and the southern part is more suitable for a Crested Ibis habitat than the northern part. At the same time, we also found that Dongzhai Reserve is located between wetlands and natural forests in Luoshan County (Figure 9a). It is a very advantageous geographical location, and an excellent place for the Crested Ibis to inhabit.
The extensive coverage of the ecological function of reserves protects most habitats of high suitability. The Dongzhai Reserve is in the southern part of Luoshan County and is covered by a biodiversity reserve, and a water conservation area also covers the area where the wetlands are located (Figure 9b). The water conservation area in Luoshan County will not only protect the local production and living environment, but also better guarantee the habitat suitability for the survival of the Crested Ibis. It shows that Luoshan County is relatively good in terms of environmental protection, ensuring the survival of the Crested Ibis.

4. Discussion

We collected a large number of the Crested Ibis trajectory data in the Henan and Shaanxi provinces of China from 2014 to 2021 with the tracking backpacks. By analyzing the data, we find that the Crested Ibis prefer cultivated lands close to forest land, paddy fields, and ponds, for the night. This is associated with the food provided by paddy fields and ponds for the Crested Ibis [6,7,9,12,15]. The areas with a high concentration of foraging points can provide more night roosting points, so the night roosting points of the Crested Ibis are close to areas with a high concentration of paddy fields and ponds. Another finding is that the activity areas of the Crested Ibis are close to the north and south in this study. The principal reason is that the Crested Ibis need tall trees for egg-laying and nesting during the breeding season [18], and there are many fruit groves and economic forests in the north and south of the Yang-xian Reserve in Shaanxi Province providing a large amount of forest land for roosting.
There are numerous factors influencing Crested Ibis’ habitat selection. The analysis shows that NDVI makes the largest contribution to the habitat assessment model and largely reflects habitat quality. We believe the NDVI helps interpret the impact of the habitat of the Crested Ibis, and it also validates the fact that vegetation changes have an important effect on the survival of the Crested Ibis [10,14]. It can also explain the conflicts with the findings of Zheng, L.’s research suggesting that altitude is the main factor affecting the habitat of the Crested Ibis, without considering the vegetation [22]. Besides the vegetation, the altitude and distance from rivers also influence the Crested Ibis’ habitat selection. The high-altitude areas are consistent with the Crested Ibis’ habit of choosing nocturnal roosting [6], and rivers can provide more fish, shrimps, and crabs [9,12]. Taking Luoshan County as an example, there are large areas (including many protected areas) suitable for the Crested Ibis to survive, especially in the southern areas. Considering the fact that the Crested Ibis usually reuse the same breeding site without frequent human disturbance [13], we advocate the protection of paddy fields and reservoirs for local farmers where the Crested Ibis once stayed, encourage active water-storage in winter, and prohibit the extensive use of pesticides for maintaining a stable ecological environment.
Human behavior also has a considerable impact on the habit selection of the Crested Ibis [11,12,43], and we believe this kind of impact should be evaluated more comprehensively. On the one hand, some artificial features could provide adequate food for the Crested Ibis; on the other hand, human activities might disturb the life of the Crested Ibis [22]. There are developed paddy fields and ponds along the roads, which are conducive to the Crested Ibis’ survival. However, human activities probably disturb the life of the Crested Ibis, and the continued expansion of the road network could also threaten the survival of the Crested Ibis [44]. Thus, revealing the relationship between Crested Ibis survival and human activities needs further study, and it might depend on specific cases.
There is a potential problem that might affect the accuracy of the analysis, which is the uneven data distribution caused by the missing trajectory data [45,46,47]. To deal with the problem, we built an LSTM module to supplement missing data and ensure its integrity in this study. Compared to the Skip-gram model [45], the LSTM module is more suitable for predicting the trajectories of individual Crested Ibis. The LSTM model is able to preserve information during long periods for learning-order dependence in sequence prediction problems [48]. This result also validates Wijeyakulasuriya, D. A.’s study, which adopted LSTM to predict the trajectory of the Larus fuscus [47]. Although the LSTM module can supplement the trajectory data, this supplementary process also generates errors for habitat prediction. Thus, in the future study, we will evaluate the prediction and obtain more Crested Ibis trajectory-data for broader analysis.

5. Conclusions

Based on the spatial and temporal trajectory-data and multi-source geographic and environmental data of the Crested Ibis, this paper uses the clustering method suitable for Crested Ibis to determine its habitat. Moreover, we analyze the biological patterns and behavioral patterns of the Crested Ibis, including the characteristics of night roosting points, the types of roosting vegetation, and the distance analysis between night roosting points and foraging points. Since the trajectories of the Crested Ibis are missing, to different degrees, an LSTM network was constructed to supplement the missing trajectory data. Additionally, this paper predicts the habitat distribution of the Crested Ibis in Luoshan County using a random forest method adjusted for its suitability, and discusses the potential habitat conditions, the impact of human activities, and the current habitat status and conservation status of the surrounding environment.
This paper provides a complete analysis tool and process for the analysis of spatial and temporal trajectory-data of the Crested Ibis, and provides methods and suggestions for experts to make future decisions.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ani13020237/s1, Figure S1: Loss graph. The horizontal coordinate is the epoch in the training process, and the vertical coordinate is the output value of the loss function; Figure S2: Clustering schematic—DBSCAN (left) Kmeans (right); Figure S3: optimal parameters: silhouette score (left) signal-to-noise ratio (right). The horizontal coordinates are Minpts, ranging from 1 to 20; the vertical coordinates are the silhouette score and the noise ratio. The silhouette score was close to 1 and the noise ratio was relatively small when Minpts was 3; Figure S4: Error bar graph. The horizontal coordinate is the group ID, and the vertical coordinate is the MSE score. Error bars represent the standard deviation of scores for each group. Since the error of group 10 is too large, it is not shown in the figure; Table S1: Random Forest Parameter Tuning. n_estimators is the number of subtrees built for the whole model; min_samples_split is the minimum number of samples required for node division; min_samples_leaf is the minimum number of leaf nodes; max_features is the maximum number of features selected; max_depth is the maximum depth of the tree; bootstrap is whether the sampling method selects bootstrap; Table S2: Random forest optimal parameters.

Author Contributions

Conceptualization, Z.C., X.J. and Y.Z.; methodology, Z.C., X.J. and Y.Z.; software, Z.C. and Y.Z.; validation, X.J., Z.C. and Y.Z.; formal analysis, Y.Z.; writing—original draft preparation, Z.C., X.J. and Y.Z.; writing—review and editing, X.J., Z.C. and Y.Z.; visualization, Z.C., X.J. and Y.Z.; funding acquisition, Z.C. and X.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (41871305, 42122025), the Opening Fund of Key Laboratory of Geological Survey and Evaluation of Ministry of Education (GLAB2022ZR06) and Fundamental Research Funds for the Central Public Welfare, Visual Analysis of Animal Trajectory Monitoring Data-Taking Crested Ibis as an Example (CAFYBB2021SY008).

Informed Consent Statement

Not applicable.

Data Availability Statement

The trajectory data are available upon request from the corresponding author. Additional environmental data can be obtained from the data sources in Table 2.

Acknowledgments

Special thanks go to the editor and reviewers for their insightful comments and constructive suggestions that substantially improved the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mace, G.M.; Collar, N.; Cooke, J.; Gaston, K.; Ginsberg, J.; Leader-Williams, N.; Maunder, M.; Milner-Gulland, E. The Development of New Criteria for Listing Species on the IUCN Red List. Verlag Nicht Ermittelbar. 1992. Available online: https://www.researchgate.net/profile/Georgina-Mace/publication/236679197_The_development_of_new_criteria_for_listing_species_on_the_IUCN_Red_List/links/0046352cece3a00c66000000/The-development-of-new-criteria-for-listing-species-on-the-IUCN-Red-List.pdf (accessed on 6 December 2022).
  2. Fahrig, L. How much habitat is enough? Biol. Conserv. 2001, 100, 65–74. [Google Scholar] [CrossRef]
  3. Li, X.H.; Li, D.M.; Ding, C.Q.; Cao, Y.H.; Lu, X.R.; Fu, W.K.; Ma, Z.J.; Lu, B.Z.; Zhai, J. A preliminary evaluation of the habitat quality of the Crested ibis (Nipponia nippon). Biodiv. Sci. 1999, 7, 161. [Google Scholar]
  4. IUCN. The IUCN Red List of Threatened Species. Version 2019-3. IUCN Global Species Programme Red List Unit United Kingdom. 2019. Available online: https://www.iucnredlist.org/ (accessed on 10 December 2019).
  5. Li., J.; Zhu, X. The value of Crested Ibis and its status and role in natural museums. In Proceedings of the International Crested Ibis Conservation Symposium, Shaanxi, China, 6 September 1999. [Google Scholar]
  6. Li, X.; Tian, H.; Li, D.J.B. Why the crested ibis declined in the middle twentieth century. Biodivers. Conserv. 2009, 18, 2165–2172. [Google Scholar] [CrossRef]
  7. Ma, L.; Li, X.; Zhai, T.; Zhang, Y.; Song, K.; Holyoak, M.; Sun, Y.J.A. Changes in the Habitat Preference of Crested Ibis (Nipponia nippon) during a Period of Rapid Population Increase. Animals 2021, 11, 2626. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Wang, F.; Cui, Z.; Li, M.; Li, X.; Ye, X.; Yu, X. Can we reestablish a self-sustaining population? A case study on reintroduced Crested Ibis with population viability analysis. Avian Res. 2021, 12, 1–10. [Google Scholar] [CrossRef]
  9. Sun, Y.; Skidmore, A.K.; Wang, T.; van Gils, H.A.; Wang, Q.; Qing, B.; Ding, C. Reduced dependence of Crested Ibis on winter-flooded rice fields: Implications for their conservation. PLoS ONE 2014, 9, e98690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Wang, Q.; Zhang, J.; Yan, W.B.; Wu, J. Spatial-temporal characteristics of vegetation cover in typical habitat of Crested Ibis. Ecol. Sci. 2019, 38, 193. [Google Scholar]
  11. Yu, F.; Sun, Y.; Wang, T.; Skidmore, A.K.; Ding, C.; Ye, X. Linking the past and present to predict the distribution of Asian crested ibis (Nipponia nippon) under global changes. Integr. Zool. 2022, 17, 1095–1105. [Google Scholar] [CrossRef]
  12. Xinhai, L.; Zhijun, M.; Changqing, D.; Tianqing, Z.; Dianmo, L. Relationship between the distribution of Crested Ibis and local farmers. Acta Zool. Sin 2002, 48, 725–732. [Google Scholar]
  13. Dong, R.; Ye, X.; Zhong, L.; Li, X.; Li, M.; Wang, H.; Yu, X. Effects of breeding success, age and sex on breeding dispersal of a reintroduced population of the Crested Ibis (Nipponia nippon) in Ningshan County, China. Avian Res. 2018, 9, 1–7. [Google Scholar] [CrossRef] [Green Version]
  14. Endo, C.; Nagata, H. Seasonal changes of foraging habitats and prey species in the Japanese Crested Ibis Nipponia nippon reintroduced on Sado Island, Japan. Bird Conserv. Int. 2013, 23, 445–453. [Google Scholar] [CrossRef] [Green Version]
  15. Li, L.; Hu, R.; Huang, J.; Bürgi, M.; Zhu, Z.; Zhong, J.; Lü, Z. A farmland biodiversity strategy is needed for China. Ecol. Evol. 2020, 4, 772–774. [Google Scholar] [CrossRef] [PubMed]
  16. Li, X.; Fu, J.; Zhai, T.; Zhang, Y.; Bruford, M.W.; Sun, Y.; Zhan, X. Understanding recovery is as important as understanding decline: The case of the crested ibis in China. Land 2022, 11, 1817. [Google Scholar] [CrossRef]
  17. Huo, Z.; Guo, J.; Li, X.; Yu, X. Post-fledging dispersal and habitat use of a reintroduced population of the Crested Ibis (Nipponia nippon). Avian Res. 2014, 5, 1–9. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, Y.; Ye, Y.; Zhang, Y.; Barras, A.G.; Wang, C.; Qing, B.; Ding, C. Tall trees drive the nest-site selection of wild Crested Ibis Nipponia nippon. Bird Conserv. Int. 2022, 32, 1–12. [Google Scholar] [CrossRef]
  19. Rew, J.; Park, S.; Cho, Y.; Jung, S.; Hwang, E. Animal movement prediction based on predictive recurrent neural network. Sensors 2019, 19, 4411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Ji, G.L.; Zhao, J. Research Progress in Pattern Mining for Big Spatio-temporal Trajectories. J. Data Acquis. Process. 2015, 30, 47–58. [Google Scholar]
  21. Jiang, X.; Yang, T.; Liu, D.; Zheng, Y.; Chen, Y.; Li, F. An Automatic Identification Method of Crested Ibis (Nipponia nippon) Habitat Based on Spatiotemporal Density Detection. Animals 2022, 12, 2220. [Google Scholar] [CrossRef] [PubMed]
  22. Zheng, L.; Wang, Y.; Zhu, J.; Wang, K.; Cai, D.; Qin, Y.; Guo, Y.; Niu, H.; Bu, Y. Habitat evaluation for reintroduced crested Ibis (Nipponia nippon) in Dongzhai national nature reserve, China, based on a maximum entropy model. Pak. J. Zool. 2018, 50, 1319–1327. [Google Scholar] [CrossRef]
  23. Sun, Y.; Wang, T.; Skidmore, A.K.; Palmer, S.C.; Ye, X.; Ding, C.; Wang, Q. Predicting and understanding spatio-temporal dynamics of species recovery: Implications for Asian crested ibis Nipponia nippon conservation in China. Divers. Distrib. 2016, 22, 893–904. [Google Scholar] [CrossRef] [Green Version]
  24. Zhou, X.H.; Jiang, L.; Wang, Q.; Zhang, W.; Cheng, K.; Ding, J. The tolerance of human disturbance on Crested Ibis (Nipponia nippon) of wardering stage. Acta Ecol. Sin. 2009, 29, 5176–5184. [Google Scholar]
  25. Lu, S.H.; Huang, Z.X.; Yuan, G.J.; Cai, D.J.; Zhao, H.; Mei, X.X.; Wang, J.; Zhu, J. Influencing factors the reproduction of wild of Nipponia nippon in Dongzhai nature reserve, Henan Province. Ecol. Sci. 2021, 40, 184. [Google Scholar]
  26. Hu, H.; Tang, J.; Wang, Y.; Zhang, H.; Lin, Y.; Su, L.; Liu, Y.; Zhang, W.; Wang, C.; Wu, D. Evaluating bird collision risk of a high-speed railway for the crested ibis. Transp. Res. Part D Transp. Environ. 2020, 87, 102533. [Google Scholar] [CrossRef]
  27. Li, M.; Tuohetahong, Y.; Lin, F.; Dong, R.; Wang, H.; Wu, X.; Ye, X.; Yu, X. Factors affecting post-release survival and dispersal of reintroduced Crested Ibis (Nipponia nippon) in Tongchuan City, China. Avian Res. 2022, 13, 100054. [Google Scholar] [CrossRef]
  28. Liu, D.; Zhang, G.; Wang, C.; Qing, B.; Lu, J. Breeding variation in a reintroduced crested ibis nipponia nippon population in Central China. Pak. J. Zool. 2020, 52, 1595–1598. [Google Scholar] [CrossRef]
  29. Okahisa, Y.; Nagata, H. Evaluation of ongoing Crested Ibis (Nipponia nippon) reintroduction using an integrated population model and Bayesian population viability analysis. IBIS 2022, 164, 1104–1122. [Google Scholar] [CrossRef]
  30. Phillips, S.J.; Dudík, M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography 2008, 31, 161–175. [Google Scholar] [CrossRef]
  31. Zhu, G.; Qiao, H. Effect of the Maxent model’s complexity on the prediction of species potential distributions. Biodivers. Sci. 2016, 24, 1189. [Google Scholar] [CrossRef] [Green Version]
  32. Geospatial Data Cloud. Available online: https://www.gscloud.cn/search (accessed on 26 December 2022).
  33. GlobeLand30. Available online: http://www.globeland30.org/defaults_en.html?src=/Scripts/map/defaults/En/download_en.html&head=download&type=data (accessed on 26 December 2022).
  34. National Catalogue Service for Geographic Information. Available online: https://www.webmap.cn/commres.do?method=result25W (accessed on 26 December 2022).
  35. OpenStreetMap Database. Available online: https://www.openstreetmap.org/export#map=12/32.1715/114.6279 (accessed on 26 December 2022).
  36. Robusto, C.C. The cosine-haversine formula. Am. Math. Mon. 1957, 64, 38–40. [Google Scholar] [CrossRef]
  37. Wang, K.F.; Shi, D.C. Observation of Anniversary and Daily Activity of the Created Ibis. In Proceedings of the International Crested Ibis Conservation Symposium, Shaanxi, China, 6 September 1999. [Google Scholar]
  38. Ester, M.; Kriegel, H.-P.; Sander, J.; Xu, X. A Density-Based Algorithm for Discovering Clusters in Large Spatial Databases with Noise. In Proceedings of the KDD’96: Proceedings of the Second International Conference on Knowledge Discovery and Data Mining, Portland, OR, USA, 2–4 August 1996; AAAI Press: Palo Alto, CA, USA, 1996; pp. 226–231. [Google Scholar]
  39. Zaki, M.J.; Meira, W., Jr. Data Mining and Machine Learning: Fundamental Concepts and Algorithms; Cambridge University Press: Cambridge, UK, 2020; pp. 375–380. [Google Scholar]
  40. Har-Peled, S.; Sadri, B. How fast is the k-means method? Algorithmica 2005, 41, 185–202. [Google Scholar] [CrossRef] [Green Version]
  41. Rousseeuw, P.J. Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. JCAM 1987, 20, 53–65. [Google Scholar] [CrossRef] [Green Version]
  42. Breiman, L. Random forests. Mach Learn 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, C.; Liu, D.; Qing, B.; Ding, H.; Cui, Y.; Ye, Y.; Lu, J.; Yan, L.; Ke, L.; Ding, C. The current population and distribution of wild Crested Ibis Nipponia nippon. Chin. J. Zool. 2014, 49, 666–671. [Google Scholar]
  44. Tao, A.M.Y.W.S.; Lü, Y.K.H.W.Z. Research status and challenges of road impacts on wildlife in China. Biodivers. Sci. 2022, 22209. [Google Scholar]
  45. Ardakani, I.; Hashimoto, K.; Yoda, K. Context-based semantical vector representations for animal trajectories. Adv. Robot. 2019, 33, 118–133. [Google Scholar] [CrossRef]
  46. Hirakawa, T.; Yamashita, T.; Tamaki, T.; Fujiyoshi, H.; Umezu, Y.; Takeuchi, I.; Matsumoto, S.; Yoda, K. Can AI predict animal movements? Filling gaps in animal trajectories using inverse reinforcement learning. Ecosphere 2018, 9, e02447. [Google Scholar] [CrossRef]
  47. Wijeyakulasuriya, D.A.; Eisenhauer, E.W.; Shaby, B.A.; Hanks, E.M. Machine learning for modeling animal movement. PLoS ONE 2020, 15, e0235750. [Google Scholar] [CrossRef]
  48. Gers, F.A.; Schmidhuber, J.; Cummins, F. Learning to forget: Continual prediction with LSTM. Neural Comput. 2000, 12, 2451–2471. [Google Scholar] [CrossRef]
Figure 1. The main methodological framework of the study.
Figure 1. The main methodological framework of the study.
Animals 13 00237 g001
Figure 2. Location of the study area; ① in Shaanxi Province and ② in Henan Province.
Figure 2. Location of the study area; ① in Shaanxi Province and ② in Henan Province.
Animals 13 00237 g002
Figure 3. Trajectory supplementation diagram. The red rectangle in (a) is enlarged to (b). The value of each pixel unit in this map is the number of trajectory points with a resolution of 200 m; the higher the number, the darker the color.
Figure 3. Trajectory supplementation diagram. The red rectangle in (a) is enlarged to (b). The value of each pixel unit in this map is the number of trajectory points with a resolution of 200 m; the higher the number, the darker the color.
Animals 13 00237 g003
Figure 4. Trajectory supplementation diagram (detail). The picture shows the trajectory data for one day: (a) denotes the missing data, (b) denotes the supplemented data. The value of each pixel unit is the number of trajectory points.
Figure 4. Trajectory supplementation diagram (detail). The picture shows the trajectory data for one day: (a) denotes the missing data, (b) denotes the supplemented data. The value of each pixel unit is the number of trajectory points.
Animals 13 00237 g004
Figure 5. (a) Night roosting points in Hantai District, Chenggu County, and Yang County. (b) Distribution of activity periods at night roosting points. (a) The gray areas indicate the extent of their night roosting points. The darker the red, the greater the number of night roosting points. (b) The reproductive period lasts from February to June, the colonial period from July to October, and the wintering period from November to January.
Figure 5. (a) Night roosting points in Hantai District, Chenggu County, and Yang County. (b) Distribution of activity periods at night roosting points. (a) The gray areas indicate the extent of their night roosting points. The darker the red, the greater the number of night roosting points. (b) The reproductive period lasts from February to June, the colonial period from July to October, and the wintering period from November to January.
Animals 13 00237 g005
Figure 6. Analysis of the distance between the night roosting points and the foraging points. The Euclidean-distance analysis of the night roosting points in this region shows that the darker the color of the foraging points, the closer the foraging point is to the night roosting points.
Figure 6. Analysis of the distance between the night roosting points and the foraging points. The Euclidean-distance analysis of the night roosting points in this region shows that the darker the color of the foraging points, the closer the foraging point is to the night roosting points.
Animals 13 00237 g006
Figure 7. Crested Ibis habitat-rating map. The red indicates the most suitable places for a high habitat-rating; the blue indicates the least suitable places, for a low habitat-rating.
Figure 7. Crested Ibis habitat-rating map. The red indicates the most suitable places for a high habitat-rating; the blue indicates the least suitable places, for a low habitat-rating.
Animals 13 00237 g007
Figure 8. Impact factor contributions. The vertical coordinates River_Dist and Road_Dist indicate the Euclidean-distance assessment for roads and rivers.
Figure 8. Impact factor contributions. The vertical coordinates River_Dist and Road_Dist indicate the Euclidean-distance assessment for roads and rivers.
Animals 13 00237 g008
Figure 9. (a) High-suitability habitat; (b) Ecological function reserve.
Figure 9. (a) High-suitability habitat; (b) Ecological function reserve.
Animals 13 00237 g009
Table 1. Crested Ibis Data Sheet. ID is the number of each individual; date is the time range of each individual’s data recording; numbers of trajectory is the number of trajectory points for each individual.
Table 1. Crested Ibis Data Sheet. ID is the number of each individual; date is the time range of each individual’s data recording; numbers of trajectory is the number of trajectory points for each individual.
IDDateLocationNumbers of Trajectory
4B04A016 July 2014–30 December 2021Luoshan County, Henan Province6541
4B03DB16 July 2014–30 December 2021Luoshan County, Henan Province1878
CAFL00424 July 2015–30 December 2021Hantai District, Chenggu County, Yang County in Shaanxi Province18,600
801324 July 2015–30 December 2021Yang County in Shaanxi Province220
801424 July 2015–30 December 2021Yang County in Shaanxi Province2769
802224 July 2015–30 December 2021Chenggu County, Yang County in Shaanxi Province947
802624 July 2015–30 December 2021Chenggu County, Yang County in Shaanxi Province870
Table 2. Environmental data sets used in the analyses, and their sources and resolution. Resources can be downloaded from the website in the citation following the data name.
Table 2. Environmental data sets used in the analyses, and their sources and resolution. Resources can be downloaded from the website in the citation following the data name.
NameSourcesResolution(m)
Google MapsSatellite LANDSAT-8, QuickBird, etc.0.5 m
DEM (digital elevation model) [32]Geospatial Data Cloud30 m
GlobeLand30 land cover data [33]GlobeLand3030 m
NDVI (normalized difference vegetation index) [32]Geospatial Data Cloud250 m
River vector data [34]National Catalogue Service for Geographic Information/
Road vector data [35]OpenStreetMap database/
Table 3. Selection rules of study periods for location points.
Table 3. Selection rules of study periods for location points.
Location PointsMarch–OctoberNovember–February
Night roost points21:00–5:0020:00–6:00
Foraging points8:00–16:009:00–16:00
Outing points5:00–8:00
16:00–21:00
6:00–7:00
16:00–20:00
Table 4. Evaluation indicators. R2 is the degree of model fit; RMSE is the root-mean-square error; MAE is the mean absolute error; explained variance is the degree of explanation of the model.
Table 4. Evaluation indicators. R2 is the degree of model fit; RMSE is the root-mean-square error; MAE is the mean absolute error; explained variance is the degree of explanation of the model.
Precision Evaluation IndicatorsR2RMSEMAEExplained_Variance
Value84.9%21.012.685.7%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, Y.; Jiang, X.; Chen, Z. Spatial- and Temporal-Trajectory Analysis of the Crested Ibis (Nipponia nippon) by Fusing Multiple Sources of Data. Animals 2023, 13, 237. https://doi.org/10.3390/ani13020237

AMA Style

Zhou Y, Jiang X, Chen Z. Spatial- and Temporal-Trajectory Analysis of the Crested Ibis (Nipponia nippon) by Fusing Multiple Sources of Data. Animals. 2023; 13(2):237. https://doi.org/10.3390/ani13020237

Chicago/Turabian Style

Zhou, Yulong, Xian Jiang, and Zhanlong Chen. 2023. "Spatial- and Temporal-Trajectory Analysis of the Crested Ibis (Nipponia nippon) by Fusing Multiple Sources of Data" Animals 13, no. 2: 237. https://doi.org/10.3390/ani13020237

APA Style

Zhou, Y., Jiang, X., & Chen, Z. (2023). Spatial- and Temporal-Trajectory Analysis of the Crested Ibis (Nipponia nippon) by Fusing Multiple Sources of Data. Animals, 13(2), 237. https://doi.org/10.3390/ani13020237

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop