- 1Key Laboratory of Earth Observation of Hainan Province, Hainan Aerospace Information Research Institute, Wenchang, China
- 2Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing, China
- 3College of Resources and Environment, University of Chinese Academy of Sciences, Beijing, China
- 4Hainan Jianfengling National Key Field Station for Forest Ecosystem, Research Institute of Tropical Forestry, Chinese Academy of Forestry, Guangzhou, China
Satellite remote sensing data is essential for large-scale, timely, and repeatable monitoring of forest species diversity. While various methods have been applied to satellite-based diversity estimation at regional scales, selecting suitable sensor and monitoring period remains challenging, especially in tropical forests. This study aims to identify the optimal time window, spatial resolution, and metrics for species diversity estimation in the Jianfengling tropical forest in southern China. We constructed stepwise linear regression models for estimating Richness, Simpson, and Shannon-Wiener indices using in-situ species diversity and heterogeneity metrics of spectra and structure. For analyzing phenology influence, we utilized six Sentinel-2 images acquired bimonthly from January to November. For evaluating scale dependency, we resampled the GF2 image to five spatial resolutions ranging from 0.8 to 10 m. The results indicated that the suitable phenological periods for species diversity estimation were at the beginning and end of the growing season, especially September performing the best for all diversity indices. Among four types of heterogeneity metrics, spectral information consistently explained most variance in species diversity indices across all periods. The optimal spatial resolution for estimating Richness and Shannon-Wiener index was 4–5 m, which corresponded to the average tree crown size. The texture features made a significant contribution compared to other metrics. Our study highlights that species diversity monitoring is highly dependent on the spatiotemporal scales of remote sensing data. It may offer practical guidance for selecting appropriate data and methods for species diversity monitoring in tropical forests.
1 Introduction
Forest species diversity is an important indicator of forest ecosystem health, which plays a vital role in maintaining ecosystem services, functions, and stability (Pennekamp et al., 2018). Tropical forests, as one of the most biologically diverse ecosystems on Earth, are crucial for mitigating climate change (Csillik et al., 2019). However, they are suffering a rapid loss of biodiversity due to increasing pressures from human disturbances, climate change, and biological invasions (Montràs-Janer et al., 2024). Moreover, the dense evergreen vegetation and complex structure characteristics make it particularly challenging to estimate species diversity accurately and understand the dynamics timely (Wiegand et al., 2017). Therefore, effective methods to monitor forest species diversity across different spatial and temporal scales are urgently required for assessing biodiversity status and guiding sustainable forest management.
Traditional forest surveys could provide accurate species diversity measurements at local scale. However, limited accessibility and high costs of large-scale in-situ data collection lead to calls for remotely sensed monitoring as a complement (Kerr and Ostrovsky, 2003). Remote sensing approaches for species diversity estimation can be broadly divided into two categories: direct classification of species or functional types and indirect estimation based on the relationships between species diversity indices and multivariate heterogeneity indicators (Turner, 2014). Many studies have explored the relationships between spectral heterogeneity metrics and species diversity indices at landscape scale utilizing airborne hyperspectral imagery, but its limited coverage and low repeatability preclude widespread usage for monitoring species diversity across larger spatial and temporal scales (Fassnacht et al., 2022; Féret and Asner, 2014). In this context, multispectral satellite data that effectively balance spatiotemporal issues provide a unique opportunity for mapping species diversity at large spatial scales, which have been successfully applied to various ecosystems (Hernández-Stefanoni et al., 2012; Lopes et al., 2017). However, their potential in tropical forest ecosystems still requires further exploration.
Some studies have estimated tropical forest diversity using spaceborne data (Ganivet and Bloomberg, 2019). For instance, Kumar et al. (2022) calculated forest species diversity using information theory-based indices derived from Sentinel-2 imagery. Njomaba et al. (2024) explored the potential of PlanetScope to predict tropical species diversity using stepwise linear regression analysis. However, it has been observed that the relationship exhibits phenological sensitivity and scale dependency, and is strongly influenced by factors such as community composition, topography, and selected heterogeneity metrics, leading to high variability with positive, negative, or no correlation (Badourdine et al., 2022; Pangtey et al., 2023; Wang et al., 2018a). As a result, understanding and assessing the influence of phenology, scale, and various metrics on species diversity estimation is crucial for achieving more reliable monitoring in tropical forests.
Vegetation phenology refers to the periodic rhythms for growth and development of plants to adapt to interannual or seasonal changes in the environment (Richardson, 2018). An increasing number of studies have used multi-temporal images covering different phenological periods to improve tree species and diversity mapping accuracy (Blickensdörfer et al., 2024; Liu et al., 2023b; Persson et al., 2018). However, the necessity of combining multi-temporal images and the determination of optimal phenological periods for species diversity estimation is still challenged. Some studies showed improved performance with multi-temporal images, especially during the transition of the growing season (Immitzer and Atzberger, 2023; Xi et al., 2021), while others reported no advantage over single-date images (Pouteau et al., 2018; Torresani et al., 2019). Additionally, the subtle phenological dynamics of evergreen vegetation and limited cloud-free imagery in tropical regions result in gaps in understanding the impact of phenology on forest species diversity estimation. Thus, it is meaningful to identify suitable time windows for monitoring species diversity in tropical forests, further guiding the acquisition of temporally matched scenes.
Spatial scale is a central topic in both ecology and remote sensing (Gamon et al., 2020). An increasing number of studies have adopted the multi-scale, multi-source framework to explore how sensor characteristics and resolution affect species diversity estimation, using data from Sentinel-2, Landsat, RapidEye, IKONOS, and WorldView-2 (Mallinis et al., 2020; Nagendra et al., 2010; Wang et al., 2022). Some studies have shown that higher spatial resolution does not necessarily improve species diversity prediction, and medium resolution (10/15 m) may be optimal (Liu et al., 2024). High-resolution imagery can introduce intra-species spectral variance due to canopy shadows, while low-resolution images may obscure inter-species spectral variance under highly mixed conditions (Rocchini et al., 2010). Therefore, finding the optimal spatial resolution is critical. Given that satellite data with very high spatial and temporal resolution (e.g. GF2 and PlanetScope) allows for better matching of pixel size with crown size, it is essential to explore its applicability and determine the suitable spatial resolution for forest species diversity monitoring.
Variations in species resource-use and growth strategies shape environmental complexity across multiple dimensions (e.g., horizontal and vertical), which can be effectively characterized through heterogeneity metrics that integrate multi-source remote sensing information. Remote sensing-based heterogeneity metrics can be generally categorized into spectral and structural heterogeneity. The spectral heterogeneity metric is a bridge between spectral diversity and species diversity, which has been proven to affect species diversity predictions (Torresani et al., 2021). Nowadays, various metrics such as coefficient of variation (CV), convex hull area (CHA), convex hull volume (CHV), and spectral angle mapper (SAM), have been proposed to capture multidimensional spectral heterogeneity (Gholizadeh et al., 2018; Kruse et al., 1993). In addition, forest structural characteristics and their heterogeneity are considered to be proxies of forest species diversity (Ma et al., 2022; Zhao et al., 2018). Texture features, as expressions of spatial structure, also play a significant role in species diversity estimation (Liu et al., 2023a). However, no single heterogeneity metric was found to be universally applicable across all species diversity estimation scenarios. Therefore, it is essential to assess the importance of different metrics in estimating species diversity, particularly considering the variations under various phenological phases and resolutions.
In this study, we evaluated the impact of phenology and spatial resolution on species diversity estimation in a tropical rainforest using multi-temporal Sentinel-2 images and very-high-resolution GF2 images, then compared the performance of various metrics across different periods and resolutions. Specifically, we aimed to address the following questions: (i) Which is the crucial phenological stage for canopy species diversity estimation in tropical forests? (ii) What is the best spatial resolution for species diversity estimation? (iii) How do heterogeneity metrics quantitatively contribute to species diversity monitoring?
2 Materials and methods
2.1 Study area
The study area is located in the Jianfengling National Nature Reserve in the southwest of Hainan province, China (JFL, 18°20′-18°57′N, 108°41′-109°12′E), covering an area of approximately 0.54 km2 (Figure 1). It is a typical representative of China’s tropical rainforest and belongs to the northern margin type of Asian tropical rainforest (Zhai et al., 2013). The complex and specialized structure of the rainforest provides conditions for a high diversity of plants. The dominant vegetation type is tropical montane rainforest, varying along the elevation gradient from 600 to 1200 m above sea level. This region is characterized by a tropical island monsoon climate, with an average annual temperature of 19.8°C and annual precipitation ranging from 1300 to 3700 mm. It shows distinct dry and wet seasons, with the wet season from May to October and the dry season from November to April (Xu et al., 2015; Zhu, 2017). The forest canopy across this study area comprises more than 30 dominant evergreen tree species, including Gironniera subaequalis, Quercus patelliformis, Alniphyllum fortune, Castanopsis chinensis, Symplocos anomala, Dacrycarpus imbricatus, and Castanopsis fissa, with varied flowering and fruiting periods (Supplementary Figure S1).

Figure 1. Overview of the study area and field data. (a) Location of the JFL National Nature Reserve with ChinaCover land cover data (Wu et al., 2017); (b) Location of the study area and sample plots with GF2 image; (c) LiDAR-derived canopy height model (CHM); (d) in-situ photos.
2.2 Field measurements and diversity metrics
Field measurements were conducted from December 12 to 24 in 2023, across 20 sample plots (20 × 20 m) within the study area (Figure 1). The coordinates of the four corners of each plot were determined by integrating the Real Time Kinematic (RTK) GPS/GLONASS System with Total Station, with errors under 10 cm. In each plot, tree parameters including species name, diameters at breast height (DBH), crown classes (dominant, co-dominant, intermediate, and suppressed trees), tree height, and crown diameters in two directions (east-west and south-north) were measured to identify upper canopy trees and calculate in-situ species diversity. All individual trees with DBH ≥ 5 cm were recorded and upper canopies located on the first or second layers (dominant and co-dominant trees) were used to calculate species diversity. The average crown diameter was 4–5 m and some plots contained large trees (Figure 2).
Based on the field measurements, three commonly used species diversity metrics, namely Richness (S) (Gaston, 2000), Shannon-Wiener index (H) (Shannon, 1948), and Simpson index (D) (Simpson, 1949), were calculated within each plot to quantify different aspects of species diversity (Supplementary Table S1). Among them, Richness represents the number of different species. The Shannon-Wiener and Simpson indexes share certain similarities as both take species evenness into account, but the former focuses on rare species, whereas the latter emphasizes dominant species. The calculation formula of the Shannon-Wiener (H) and Simpson (D) were as Equations 1, 2:
where S is the total number of species in plot, and is the proportional abundance of the species , represented by the ratio of individuals of species within a plot.
2.3 Remote sensing data
2.3.1 Satellite data
Sentinel-2 can effectively capture spectral changes of vegetation across different periods due to its high revisit frequency of 5–10 days, enabling us to analyze the impact of phenology on species diversity estimation models. Considering the cloudy and rainy weather in Southern China, we downloaded all nearly cloud-free Sentinel-2 surface reflectance imagery over the study area from 2022 to 2024 to cover the entire phenological cycle. These level-2A images, consist of 13 wavebands covering visible, near- and shortwave-infrared spectra, were atmospherically corrected by ESA using the Sen2Cor algorithm (Louis et al., 2016). Finally, we obtained six Sentinel-2 images representing different phenological periods (Table 1). Due to cloud shading obscuring one plot in the July image, only 19 plots were used for this period. We also conducted a visual matching inspection between Sentinel-2 and UAV imagery to ensure that geolocation shift is acceptable. For analysis, we excluded the 60 m atmospheric bands since they do not contain surface information, and resampled the 20 m bands to 10 m using bilinear interpolation based on SNAP software. Then, after analyzing the NDVI distribution histogram and manually inspecting the NDVI values of non-vegetation and shadow pixels, a NDVI threshold (0.3) was applied to mask most of the non-forested pixels while preserving forest pixels from scenes during non-growing seasons.
GF2 is an optical satellite with sub-meter spatial resolution and revisits within 5 days, which can capture vegetation information effectively. GF2 images feature one panchromatic band (0.8 m) and four VIS-NIR multispectral bands (blue, green, red, and near-infrared) with a 4 m spatial resolution, offering valuable data for our analysis of spatial scales (Zhou et al., 2020). We selected an available scene of cloud-free GF2 image (29/01/2023, download from Natural Resources Satellite Remote Sensing Cloud Service Platform) that was closest to the field sampling period and performed image pre-processing, including radiometric calibration, atmospheric correction (FLAASH module based on MODTRAN radiative transfer model), geometric correction, image fusion, and cropping, using ENVI 5.3 software (Xia et al., 2023). Finally, the NDVI threshold of 0.5 was applied to mask shadows. According to the average crown diameter, we resampled the fused image from the original (0.8 m) to coarser spatial resolutions (3 m, 4 m, 5 m, and 10 m) utilizing the nearest neighbor resampling algorithm (Liu et al., 2024). It is worth noting that there were only two scenes of GF2 data available from 2022 to 2024 due to persistent cloud and rain cover, which were insufficient to support phenological analysis. Therefore, GF2 data was only used for multi-resolution analysis.
2.3.2 LiDAR data
Given the absence of spaceborne LiDAR coverage in the study area, we employed unmanned aerial vehicle (UAV) LiDAR as an alternative to analyze the importance of structural heterogeneity metrics for species diversity estimation. The UAV LiDAR data was simultaneously collected with field surveys using an FT-800H laser scanner (LuoJiaYiYun Optoelectronic Technology Co., Ltd., Wuhan, China) mounted on a DJI M300 UAV platform. The scanner provided a wide field of view of 330° and a measurement accuracy of 1 cm. The average point density was 1,058 points/m2. The UAV LiDAR data were denoised and filtered to generate a digital elevation model (DEM) based on classified ground points, and a digital surface model (DSM) based on the first pulse reflections. The canopy height model (CHM) was derived by subtracting the DEM from the DSM at 0.5 m resolution. UAV LiDAR data processing was performed using LiDAR360 V7.2 (GreenValley International Inc., Beijing, China). Given the subtle structural changes, a single period of structural features was sufficient for analysis.
2.4 Spectral and structural heterogeneity metrics
After synthesizing the rationale and applicability of commonly used metrics, four types of heterogeneity metrics were selected: spectral information metrics, spectral diversity metrics, texture features, and structural diversity metrics (Table 2). Spectral information metrics derived from spectral bands or vegetation indices (VIs) directly represent spectral reflectance features and highlight specific properties of vegetation. After comprehensively considering variable redundancy and spatial resolution of different bands, six representative VIs from Sentinel-2 and the corresponding three from GF2 were calculated to characterize pigment content, specific leaf area, and water content. We adopted the coefficients of these commonly used vegetation indices, as they have been validated for both sensors in previous studies (Liu et al., 2024; Xi et al., 2023; Zhou et al., 2020). For normalized difference spectral index (NDSI), we modified the band combination to reflect photosynthetic pigments and vegetation growth (Patil et al., 2007). The average values of these spectral variables (i.e., bands and VIs) were then computed for each plot. Spectral diversity metrics included CV, CHA, and SAM, where CVs were calculated based on VIs or multi-bands and CHA and SAM were calculated using all spectral bands. Eight texture features were derived from the gray level co-occurrence matrix (GLCM) based on the first principal component of Sentinel-2 and GF2 multispectral imagery. Considering the plot size, moving window sizes of GLCM were set as 27 × 27, 7 × 7, 5 × 5, 3 × 3, and 3 × 3 for the 0.8 m, 3 m, 4 m, 5 m, and 10 m resolution images respectively. For the phenological analysis based on Sentinel-2 images, we excluded the texture parameter of Cor because it was almost the same among the 20 plots. Rao’s Q index calculated based on CHM, which incorporates horizontal variation in canopy vertical structure, was selected as the representative of structural diversity (Torresani et al., 2020). We directly retrieved Rao’s Q at plot-scale using the 0.5 m CHM for temporal analysis. But in spatial scale analysis, we resampled CHM (0.8, 3, 4, 5, and 10 m) to match the optical satellite data, and then calculated Rao’s Q at plot-scale.
2.5 Species diversity estimation
To evaluate the impact of phenological stages and spatial resolution on forest species diversity estimation, a series of multiple stepwise regression analyses were conducted based on multi-temporal Sentinel-2 data and resampled GF2 data. To avoid overfitting caused by multicollinearity among variables, we calculated correlation coefficients (r) among all initial variables and removed highly correlated variables (r > 0.8) that had lower correlation with in-situ species diversity indices (Supplementary Figure S3-S13). Remaining variables were ranked based on their correlation with species diversity indices, and predictor variables were selected using forward stepwise regression. Additionally, variables with a large variance inflation factor (VIF > 10) were removed. Finally, model accuracy was assessed via leave-one-out cross-validation, using metrics of the coefficient of determination (R2), root mean square error (RMSE), and the mean absolute error (MAE). In terms of variable importance, hierarchical partitioning was used to separate the amount of variation explained by each predictor (Groemping, 2006).
3 Results
3.1 Impact of phenology
Based on best-fit models using multi-temporal Sentinel-2 data, the optimal monitoring period for Richness (R2 = 0.70, RMSE = 3.39, MAE = 2.49), Shannon-Wiener (R2 = 0.71, RMSE = 0.32, MAE = 0.26), and Simpson (R2 = 0.48, RMSE = 0.11, MAE = 0.09) consistently occurred in September (Table 3; Figure 3). Moreover, the model performance of Richness and Shannon-Wiener was notably superior to that of Simpson (Supplementary Figure S2).

Table 3. Model accuracies from multi-temporal Sentinel-2 images for Richness, Shannon-Wiener, and Simpson.

Figure 3. The temporal variation of coefficient of determination (R2) for the stepwise regression models for (a) Richness, (b) Shannon-Wiener, and (c) Simpson (*, 0.05< p-value <0.1, significant; **, 0.01< p-value <0.05, very significant; ***, p-value <0.01, extremely significant).
3.2 Impact of spatial resolution
The cross-validation accuracies of the models for estimating species diversity at different spatial resolutions are shown in Table 4 and Figure 4. GF2 images were more advantageous for accurate estimation of Richness compared to Simpson and Shannon-Wiener. The species richness estimation model based on 5 m resolution GF2 data performed the best among all spatial resolution scenarios (R2 = 0.62, RMSE = 4.11, MAE = 3.06). The accuracy for Shannon-Wiener estimation was highest at the resolution of 4 m (R2 = 0.24, RMSE = 0.53, MAE = 0.39), while the Simpson estimation model performed best at 0.8 m (R2 = 0.26, RMSE = 0.12, MAE = 0.10).

Table 4. Model accuracies from multi-spatial resolution GF2 images for Richness, Shannon-Wiener, and Simpson.

Figure 4. The spatial variation of coefficient of determination (R2) for the stepwise regression models for (a) Richness, (b) Shannon-Wiener, and (c) Simpson (*, 0.05< p-value <0.1, significant; **, 0.01< p-value <0.05, very significant; ***, p-value <0.01, extremely significant).
Specifically, the model at 0.8 m for Richness showed a considerable coefficient of determination (R2 = 0.50) and outperformed Shannon-Wiener and Simpson. Models at fine scale (0.8 m) and approaching crown size scales (4–5 m) had similar estimation accuracy for Shannon-Wiener, as same as for Simpson. Compared to the estimation based on Sentinel-2 in January, GF2 data with 10 m resolution provided a higher R2 for Simpson (R2 = 0.22), while lower R2 values for Richness and Shannon-Wiener.
3.3 Importance of heterogeneity metrics
In phenological effect analysis, spectral information metrics exhibited the highest explanatory power, accounting for 48% of variance in Richness and over 60% for Shannon-Wiener and Simpson. Texture features explained more variance in Richness (39%) than in Shannon-Wiener (12%) and Simpson (9%). Spectral diversity metrics explained 11%, 13%, and 18% of the variance for Richness, Shannon-Wiener, and Simpson, respectively. Structural diversity metrics had the lowest explanatory power, contributing 10% for Shannon-Wiener, 4% for Simpson, and 2% for Richness (Figure 5a–c). Additionally, the temporal variability of variable contributions (Figure 5d) indicated that the explanatory power of spectral information remained relatively stable across all phenological periods. However, the metric with the greatest variation in importance differed depending on diversity indices. For Richness, the structural diversity metric showed the greatest temporal variation in importance, while it was the spectral diversity metric for Shannon-Wiener, and the texture feature for Simpson.

Figure 5. The overall explanatory power of four types of heterogeneity metrics on species diversity indices in phenological effect analysis (a-c); and the explained variance’s CV of four types of heterogeneity metrics in different temporal models (d).
The hierarchical partitioning of the variation explained by each variable for species diversity models at different phenological stages is shown in Figure 6. Specifically, B2 and SR played important roles in the January models for all species diversity indices (33% and 38% for Richness, 23% and 77% for Shannon-Wiener, 30% and 70% for Simpson). B11 explained variation in Shannon-Wiener and Simpson models from May to September, and contributed to Richness models in July (26%) and September (18%). Additionally, NDWI and NDSI were involved in Richness estimation model at November, while acted on Shannon-Wiener models in May and July, and contributed to Simpson models across these three periods (May, July, and November). Moreover, the texture feature of Mean independently explained all variance in Richness in the May model, CV_NDWI accounted for all variance in Shannon-Wiener in the November model, and Hom explained all variance in Simpson in the March model.

Figure 6. Hierarchical partitioning of the variation explained for (a) Richness, (b) Shannon-Wiener, and (c) Simpson in different phenological models.
Figure 7a–c showed significant differences in the importance of four types of heterogeneity metrics for species diversity estimation at different spatial scales. Overall, texture features were the most important, explaining 40%, 48%, and 57% of the variance in Richness, Shannon-Wiener, and Simpson index, respectively. Spectral information contributed 21% to Richness, 32% to Shannon-Wiener, and 20% to Simpson. Spectral diversity performed better for Richness (25%) and Simpson (21%) than Shannon-Wiener (9%). Structural diversity showed stronger explanatory power for Richness (14%) and Shannon-Wiener (11%) than Simpson (2%). Moreover, the contribution of texture features remained stable across scales for all three diversity indices (Figure 7d). Spectral diversity varied most in contribution for Shannon-Wiener estimation across different spatial resolutions, while spectral information and structural diversity varied significantly in importance for Simpson estimation across scales.

Figure 7. The overall explanatory power of four types of heterogeneity metrics on species diversity indices in spatial scale effect analysis (a-c); and the explained variance’s CV of four types of heterogeneity metrics in different spatial models (d).
As shown in Figure 8, texture features Dis, Cont, and Mean completely explained Richness at 10 m and 4 m, and Simpson at 3 m. NDVI was the most important spectral information variable for Shannon-Wiener, contributing significantly at 0.8 m, 3 m and 4 m resolutions (34%, 46%, and 27%). CHA contributed to Richness estimation at three scales (0.8 m: 14%, 3 m: 27%, and 5 m: 33%), while only one scale for Shannon-Wiener (24% at 5 m) and Simpson (10% at 0.8 m). Rao’s Q based on CHM was involved in three spatial models for Richness (0.8 m, 3 m, and 5 m), and Shannon-Wiener (0.8 m, 4 m, and 5 m), but only one model for Simpson (0.8 m).

Figure 8. Hierarchical partitioning of the variation explained for (a) Richness, (b) Shannon-Wiener, and (c) Simpson in different spatial models.
4 Discussion
4.1 Phenological sensitivity of species diversity estimation
We identified September as the optimal phenological period, which aligns with the consensus that the growing season (May–September) is the most informative period and a primary choice for remote sensing-based diversity monitoring. However, inconsistencies exist in literature regarding the most efficient single-date time window for estimating forest species diversity. While many studies underscored the importance of transitional seasons, such as green-up (May) and senescence (October) (Liu et al., 2023a), others suggested that peak growing seasons like July can obtain higher accuracy (Arekhi et al., 2017; Chrysafis et al., 2020; Yang et al., 2022). These differences stem from the biochemical and morphological changes among species caused by crucial phenological events, such as leaf unfolding and flowering in spring, and leaf discoloration and falling in autumn (Kollert et al., 2021; Xi et al., 2023). In our tropical site, some tree species begin fruiting in September, leading to the transformation of dominant pigments in leaves and driving color changes from green to yellow and red, which enhances tree species differentiation (Kapoor et al., 2022). Although the effect of shrubs and grasses was relatively low under dense canopy cover in tropical forests, seasonal variations in non-tree vegetation or background signals may also affect the relationship between spectral heterogeneity and species diversity (Madonsela et al., 2021), potentially explaining the reduced accuracy in November. The phenological asynchrony in different forest layers responding to warming temperatures in May may increase the spectral contribution of non-canopy vegetation, further reducing accuracy. These findings underscore the importance of understory vegetation in species diversity assessments, which may be better captured in the future with advancements in sensor technologies, such as hyperspectral LiDAR (Bai et al., 2024) and multi-angle observations and satellite thermal infrared sensing (Adams et al., 2025).
Our results highlighted that late spring and autumn are critical temporal windows for forest species diversity monitoring. However, frequent cloud cover and rainfall in tropical regions increase the uncertainty in acquiring images during a single optimal period. Moreover, effective phenological time for diversity monitoring varies by vegetation type. For example, broadleaved tree species are typically more distinguishable in autumn, while conifers exhibit greater spectral changes in spring (Grabska et al., 2019). Hence, if September imagery is unavailable, May or November images can also be considered based on vegetation characteristics. We acknowledge that bimonthly temporal selection may be insufficient to fully capture subtle phenological changes in tropical evergreen vegetation. New generation satellites with higher spatiotemporal resolution are expected to address this limitation. For instance, the PlanetScope dataset characterized by daily revisit frequency and 3 m spatial resolution exhibits considerable potential to comprehensively assess forest diversity. Besides, continuous advancements in the fusion of multi-source remote sensing data offer promising opportunities for more effective monitoring of tropical forests at both high temporal frequency and fine spatial resolution. Notably, Multi-Temporal and Spectral-Temporal-Metric methods are increasingly applied to forest species diversity mapping (Rahmanian et al., 2023; Vanguri et al., 2024). We recommend performing a similar phenological sensitivity analysis before image composition to reduce data redundancy and improve accuracy.
4.2 Scale dependency of species diversity estimation
Our results indicated that the optimal spatial resolution varied among species diversity indices. The model accuracy of Richness and Shannon-Wiener was highest when the spatial resolution approaches the average crown size (4–5 m). This is because, at this resolution, the spectral properties effectively capture the chemical, morphological, and structural characteristics of individual trees (Zheng et al., 2022). However, Simpson performed better at finer resolutions (0.8 m), likely due to its sensitivity to dominant species. Considering the overall low abundance of dominant species caused by the high diversity in this study area, at resolutions of 5 m or alike, the spectral information of the coniferous dominant species is easily concealed by that of broadleaved species, while it could be preserved at fine resolution. Although finer resolutions may introduce noise, appropriate data preprocessing effectively mitigates this risk, ensuring reliable results. Another point worth noting is the number of in-situ plots, which may not fully represent the community composition and distribution in the study area. Despite our best efforts to sample as evenly as possible, the limited in-situ plots with high abundance of dominant species may impact the results.
Some previous studies have compared the impact of spatial scale on species diversity estimation using different satellite images with multiple resolutions (Gyamfi-Ampadu et al., 2021). However, differences in spectral properties of satellite data sources may cause bias in evaluating scale effects. To address this, we employed a single type of satellite data and resampling method, focusing on identifying suitable monitoring scales to provide informed recommendations. Moreover, evidence shows that resampling high-resolution imagery to coarser resolutions may improve the estimation accuracy up to a threshold, beyond which spectral mixing reduces performance (Liu et al., 2024). Additionally, a coarser resolution resulted in fewer pixels available within the plot for calculating the diversity indices, partly explaining the lower accuracy observed at 10 m resolution in our study. Beyond the critical spatial scales tested in our study, we suggest future research paying more attention to the integrated utilization of nested plots and UAV data to better understand the effects of grain and extent on forest species diversity estimation. A critical next step is to assess how spatial resolution degradation influence spectral integrity and metric performance, with UAV imagery providing a reliable ground reference. UAV provides very high-resolution data, which can match well with the nested plots at different scales, thereby enabling spatial extrapolation and selecting the optimal monitoring scale. Adjusting the field sampling design to be remotely driven may be also a viable strategy to probe deeply into scale effects. Besides, while this study was targeted at tropical forests, the methodology for scale detection is broadly applicable to other ecosystems. Notably, in areas with large terrain undulations, such as mountainous and canyon areas, terrain correction and shadow factors need to be carefully taken into consideration.
4.3 Variable importance in species diversity estimation
Our results indicated that the importance of spectral heterogeneity metrics in forest species diversity estimation has temporal and spatial effects. In our study, the blue band was prominent in January, probably associated with carotenoid uptake during leaf senescence (Hennessy et al., 2020). The near-infrared and red-edge bands were significant in late spring, attributed to their association with strong chlorophyll uptake and internal leaf structure during early growth stages. The SWIR band contributed to monitoring species diversity across several phenological periods, potentially due to its sensitivity to canopy leaf water content (Ming et al., 2024). However, for GF2, we only found a weak role of Green, Red, and NDVI, which may be related to the sensor band range and the limited number of bands. But interestingly, when resampling to 10 m resolution, the spectral band mattered more significantly to Simpson and Shannon-Wiener estimation compared to spectral diversity indices that integrate multidimensional information. This phenomenon, where more complex metrics using full-range spectral data perform worse on coarse-resolution images, has also been reported (Rossi et al., 2021; Wang et al., 2018b). In addition, spectral diversity metrics calculated by Sentinel-2 showed limited explanatory power for species diversity, while CHA based on GF2 explained a larger proportion of the variance. This discrepancy may be attributed to the coarse spatial resolution of Sentinel-2 and its matching issue with the plot size. Previous studies also reported that the relationship between species and spectral diversity is more significant at larger plot size (Marselis et al., 2019). More refined spatial matching strategies, such as increasing field plot size, applying sampling buffers or comparing the performance of satellite data extracted from different corner coordinates (Liu et al., 2023b), could contribute to mitigate the effects of positional offset. In addition, different atmospheric correction methods may affect metric consistency across time and space (Chraibi et al., 2022). While the use of normalized vegetation indices (VIs) can mitigate some atmospheric effects, band reflectance values used in spectral diversity metrics remain sensitive to correction accuracy. Given the complexity of atmospheric conditions in tropical regions, future research could investigate the consistency and uncertainty of different atmospheric correction methods, and its potential influence on spectral diversity estimation. Also, the scale aggregation led to various degrees of spectral mixing, which mitigated the importance of spectral heterogeneity metrics. Therefore, the availability of more satellite-based data with high temporal and spatial resolution or multi-temporal UAV data is called for further exploring the suitability of spectral heterogeneity metrics at different spatial and temporal scales.
Our analysis also emphasized the critical role of structural heterogeneity metrics, especially texture features, in forest species diversity estimation. These features derived from the spatial variability in image pixels can effectively reflect variations in canopy structure, leaf arrangement, and key aspects of ecosystem heterogeneity that influence forest diversity patterns. Compared to other metrics, the superior performance of texture features likely benefits from the fact that they can be computed through diverse algorithms, enabling the quantification of multidimensional image heterogeneity (Farwell et al., 2021). With the emergence of more high-precision data, texture features offer promising potential to enhance species diversity prediction models. The limited structural diversity metrics (Rao’s Q and texture) in our study may underestimate the importance of structural metrics for species diversity estimation. Thus, we suggest fully leveraging the advantages of LiDAR point cloud data and incorporating innovative structural heterogeneity indices, such as canopy entropy, to reduce saturation effects in structurally complex forests, and improve the accuracy of diversity estimation. Moreover, with the recent availability of high-resolution, wall-to-wall structural parameter datasets, such as the canopy height map at 1 m resolution (Tolan et al., 2024), the influence of structural heterogeneity can be integrated in future large-scale applications.
5 Conclusion
Overall, we assessed the impacts of phenology and spatial resolution on species diversity estimation in a typical tropical forest using Sentinel-2 and GF2 satellite data and emphasized the importance of various heterogeneity metrics across temporal and spatial scales. Our results indicate that the optimal phenological periods for estimating species diversity are at the beginning and end of the growing season, while the ideal spatial resolution aligns with the tree crown size. These findings provide guidance on selecting appropriate phenological periods and spatial scales, thereby improving species diversity monitoring in complex evergreen tropical forests to achieve more accurate and efficient estimations. Despite focusing on a specific tropical forest, our methodology could offer concrete spatiotemporal prior knowledge for integrating suitable remote sensing data and optimizing fieldwork designs in various ecosystems. Broader exploration across globally diverse tropical forests would enhance the understanding of tropical-specific phenological windows and spatial scales. Additionally, our study provides valuable insights into promoting the development of dynamic heterogeneity metrics that incorporate spatiotemporal variation and calling for in-depth analysis of the impact of scale changes on metrics themselves. Future research could explore the synergistic application of multi-source satellite imagery, UAV data, and machine learning algorithms, which may help bridge the scale gap between in-situ measurements and large-scale satellite observations, and further enhance the accuracy and scalability of species diversity monitoring.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: Multi-temporal Sentinel-2 data were retrieved from https://dataspace.copernicus.eu/. GF2 data was retrieved from http://www.sasclouds.com/. UAV LiDAR data was private and not publicly available, please contact the corresponding author via email if interested. Field plots data was available in the Supplementary Material.
Author contributions
PZ: Conceptualization, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing, Data curation, Investigation, Software, Validation. YZ: Conceptualization, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing, Funding acquisition, Project administration, Resources, Supervision. ZJZ: Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing. CX: Investigation, Methodology, Software, Visualization, Writing – review & editing. JW: Investigation, Writing – review & editing. XM: Investigation, Writing – review & editing. ZFZ: Investigation, Writing – review & editing. JC: Investigation, Writing – review & editing. TZ: Investigation, Writing – review & editing. DZ: Investigation, Resources, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the Hainan Provincial Natural Science Foundation of China (Grant No. 423CXTD389), the National Natural Science Foundation of China (Grant No. 42301410), and the China Postdoctoral Science Foundation (Grant No. 2023M743587).
Acknowledgments
We extend our heartfelt thanks to Hainan Jianfengling National Key Field Station for Forest Ecosystem for providing a data acquisition platform. We thank all colleagues for their support with the fieldwork.
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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.2025.1582910/full#supplementary-material
References
Adams, J. S., Damm, A., Werfeli, M., Gröbner, J., and Naegeli, K. (2025). Across-scale thermal infrared anisotropy in forests: Insights from a multi-angular laboratory-based approach. Remote Sens. Environ. 326, 114766. doi: 10.1016/j.rse.2025.114766
Arekhi, M., Yılmaz, O. Y., Yılmaz, H., and Akyüz, Y. F. (2017). Can tree species diversity be assessed with Landsat data in a temperate forest? Environ. Monit. Assess. 189, 586. doi: 10.1007/s10661-017-6295-6
Badourdine, C., Féret, J. B., Pélissier, R., Vincent, G., Sankey, T., and Ling, F. (2022). Exploring the link between spectral variance and upper canopy taxonomic diversity in a tropical forest: influence of spectral processing and feature selection. Remote Sens. Ecol. Conserv. 9, 235–250. doi: 10.1002/rse2.306
Bai, J., Niu, Z., and Wang, L. (2024). A theoretical demonstration on the independence of distance and incidence angle effects for small-footprint hyperspectral LiDAR: Basic physical concepts. Remote Sens. Environ. 315, 114452. doi: 10.1016/j.rse.2024.114452
Blickensdörfer, L., Oehmichen, K., Pflugmacher, D., Kleinschmit, B., and Hostert, P. (2024). National tree species mapping using Sentinel-1/2 time series and German National Forest Inventory data. Remote Sens. Environ. 304, 114069. doi: 10.1016/j.rse.2024.114069
Chraibi, E., de Boissieu, F., Barbier, N., Luque, S., and Féret, J.-B. (2022). Stability in time and consistency between atmospheric corrections: Assessing the reliability of Sentinel-2 products for biodiversity monitoring in tropical forests. Int. J. Appl. Earth Observ. Geoinform. 112:102884. doi: 10.1016/j.jag.2022.102884
Chrysafis, I., Korakis, G., Kyriazopoulos, A. P., and Mallinis, G. (2020). Predicting tree species diversity using geodiversity and sentinel-2 multi-seasonal spectral information. Sustainability 12, 9250. doi: 10.3390/su12219250
Csillik, O., Kumar, P., Mascaro, J., O’Shea, T., and Asner, G. P. (2019). Monitoring tropical forest carbon stocks and emissions using Planet satellite data. Sci. Rep. 9, 17831. doi: 10.1038/s41598-019-54386-6
El-Shikha, D. M., Waller, P., Hunsaker, D., Clarke, T., and Barnes, E. (2007). Ground-based remote sensing for assessing water and nitrogen status of broccoli. Agric. Water Manage. 92, 183–193. doi: 10.1016/j.agwat.2007.05.020
Farwell, L. S., Gudex-Cross, D., Anise, I. E., Bosch, M. J., Olah, A. M., Radeloff, V. C., et al. (2021). Satellite image texture captures vegetation heterogeneity and explains patterns of bird richness. Remote Sens. Environ. 253, 112175. doi: 10.1016/j.rse.2020.112175
Fassnacht, F. E., Müllerová, J., Conti, L., Malavasi, M., Schmidtlein, S., and Rocchini, D. (2022). About the link between biodiversity and spectral variation. Appl. Vegetation Sci. 25, e12643. doi: 10.1111/avsc.12643
Féret, J.-B. and Asner, G. P. (2014). Mapping tropical forest canopy diversity using high-fidelity imaging spectroscopy. Ecol. Appl. 24, 1289–1296. doi: 10.1890/13-1824.1
Gamon, J. A., Wang, R., Gholizadeh, H., Zutta, B., Townsend, P. A., and Cavender-Bares, J. (2020). Consideration of scale in remote sensing of biodiversity. Remote Sens. Plant Biodivers., 425–447. doi: 10.1007/978-3-030-33157-3_16
Ganivet, E. and Bloomberg, M. (2019). Towards rapid assessments of tree species diversity and structure in fragmented tropical forests: A review of perspectives offered by remotely-sensed and field-based data. For. Ecol. Manage. 432, 40–53. doi: 10.1016/j.foreco.2018.09.003
Gao, B.-C. (1996). NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 58, 257–266. doi: 10.1016/S0034-4257(96)00067-3
Gholizadeh, H., Gamon, J. A., Zygielbaum, A. I., Wang, R., Schweiger, A. K., and Cavender-Bares, J. (2018). Remote sensing of biodiversity: Soil correction and data dimension reduction methods improve assessment of α-diversity (species richness) in prairie ecosystems. Remote Sens. Environ. 206, 240–253. doi: 10.1016/j.rse.2017.12.014
Gitelson, A. A., Gritz, Y., and Merzlyak, M. N. (2003). Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 160, 271–282. doi: 10.1078/0176-1617-00887
Grabska, E., Hostert, P., Pflugmacher, D., and Ostapowicz, K. (2019). Forest stand species mapping using the sentinel-2 time series. Remote Sens. 11, 1197. doi: 10.3390/rs11101197
Groemping, U. (2006). Relative importance for linear regression in R: the package relaimpo. J. Stat. Softw. 17, 1–27. doi: 10.18637/jss.v017.i01
Gyamfi-Ampadu, E., Gebreslasie, M., and Mendoza-Ponce, A. (2021). Evaluating multi-sensors spectral and spatial resolutions for tree species diversity prediction. Remote Sens. 13, 1033. doi: 10.3390/rs13051033
Hennessy, A., Clarke, K., and Lewis, M. (2020). Hyperspectral classification of plants: A review of waveband selection generalisability. Remote Sens. 12, 113. doi: 10.3390/rs12010113
Hernández-Stefanoni, J. L., Gallardo-Cruz, J. A., Meave, J. A., Rocchini, D., Bello-Pineda, J., and López-Martínez, J. O. (2012). Modeling α- and β-diversity in a tropical forest from remotely sensed and spatial data. Int. J. Appl. Earth Observ. Geoinform. 19, 359–368. doi: 10.1016/j.jag.2012.04.002
Huete, A., Didan, K., Miura, T., Rodriguez, E. P., Gao, X., and Ferreira, L. G. (2002). Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 83, 195–213. doi: 10.1016/S0034-4257(02)00096-2
Immitzer, M. and Atzberger, C. (2023). Tree species diversity mapping—Success stories and possible ways forward. Remote Sens. 15, 3074. doi: 10.3390/rs15123074
Jordan, C. F. (1969). Derivation of leaf-area index from quality of light on the forest floor. Ecology 50, 663–666. doi: 10.2307/1936256
Kapoor, L., Simkin, A. J., George Priya Doss, C., and Siva, R. (2022). Fruit ripening: dynamics and integrated analysis of carotenoids and anthocyanins. BMC Plant Biol. 22, 27. doi: 10.1186/s12870-021-03411-w
Kerr, J. T. and Ostrovsky, M. (2003). From space to species: ecological applications for remote sensing. Trends Ecol. Evol. 18, 299–305. doi: 10.1016/s0169-5347(03)00071-5
Kollert, A., Bremer, M., Löw, M., and Rutzinger, M. (2021). Exploring the potential of land surface phenology and seasonal cloud free composites of one year of Sentinel-2 imagery for tree species mapping in a mountainous region. Int. J. Appl. Earth Observ. Geoinform. 94, 102208. doi: 10.1016/j.jag.2020.102208
Kruse, F. A., Lefkoff, A. B., Boardman, Y. J. W., Heidebrecht, K. B., Shapiro, A. T., Barloon, P. J., et al. (1993). The spectral image processing system (SIPS)—interactive visualization and analysis of imaging spectrometer data. Remote Sens. Environ. 44, 145–163. doi: 10.1016/0034-4257(93)90013-N
Kumar, P., Dobriyal, M., Kale, A., Pandey, A. K., Tomar, R. S., and Thounaojam, E. (2022). Calculating forest species diversity with information-theory based indices using sentinel-2A sensor’s of Mahavir Swami Wildlife Sanctuary. PloS One 17, e0268018. doi: 10.1371/journal.pone.0268018
Liu, X., Frey, J., Munteanu, C., Denter, M., and Koch, B. (2024). Tree species diversity mapping from spaceborne optical images: The effects of spectral and spatial resolution. Remote Sens. Ecol. Conserv. 10, 463–479. doi: 10.1002/rse2.383
Liu, X., Frey, J., Munteanu, C., Still, N., and Koch, B. (2023a). Mapping tree species diversity in temperate montane forests using Sentinel-1 and Sentinel-2 imagery and topography data. Remote Sens. Environ. 292, 113576. doi: 10.1016/j.rse.2023.113576
Liu, Y., Zhang, R., Lin, C.-F., Zhang, Z., Zhang, R., Shang, K., et al. (2023b). Remote sensing of subtropical tree diversity: The underappreciated roles of the practical definition of forest canopy and phenological variation. For. Ecosyst. 10, 100122. doi: 10.1016/j.fecs.2023.100122
Lopes, M., Fauvel, M., Ouin, A., and Girard, S. (2017). Spectro-temporal heterogeneity measures from dense high spatial resolution satellite image time series: application to grassland species diversity estimation. Remote Sens. 9, 993. doi: 10.3390/rs9100993
Louis, J., Debaecker, V., Pflug, B., Main-Korn, M., Bieniarz, J., Mueller-Wilm, U., et al. (2016). “Sentinel-2 sen2Cor: L2A processor for users,” in Proceedings Living Planet Symposium 2016 SP-740, 1–8.
Ma, Q., Su, Y., Hu, T., Jiang, L., Mi, X., Lin, L., et al. (2022). The coordinated impact of forest internal structural complexity and tree species diversity on forest productivity across forest biomes. Fundam. Res. 4, 1185–1195. doi: 10.1016/j.fmre.2022.10.005
Madonsela, S., Cho, M., Ramoelo, A., and Mutanga, O. (2021). Investigating the relationship between tree species diversity and landsat-8 spectral heterogeneity across multiple phenological stages. Remote Sens. 13, 2467. doi: 10.3390/rs13132467
Mallinis, G., Chrysafis, I., Korakis, G., Pana, E., and Kyriazopoulos, A. P. (2020). A random forest modelling procedure for a multi-sensor assessment of tree species diversity. Remote Sens. 12, 1210. doi: 10.3390/rs12071210
Marselis, S. M., Tang, H., Armston, J., Abernethy, K., Alonso, A., Barbier, N., et al. (2019). Exploring the relation between remotely sensed vertical canopy structure and tree species diversity in Gabon. Environ. Res. Lett. 14, 094013. doi: 10.1088/1748-9326/ab2dcd
Ming, L., Liu, J., Quan, Y., Li, M., Wang, B., and Wei, G. (2024). Mapping tree species diversity in a typical natural secondary forest by combining multispectral and LiDAR data. Ecol. Indic. 159, 111711. doi: 10.1016/j.ecolind.2024.111711
Montràs-Janer, T., Suggitt, A. J., Fox, R., Jönsson, M., Martay, B., Roy, D. B., et al. (2024). Anthropogenic climate and land-use change drive short- and long-term biodiversity shifts across taxa. Nat. Ecol. Evol. 8, 739–751. doi: 10.1038/s41559-024-02326-7
Nagendra, H., Rocchini, D., Ghate, R., Sharma, B., and Pareeth, S. (2010). Assessing plant diversity in a dry tropical forest: comparing the utility of landsat and ikonos satellite images. Remote Sens. 2, 478–496. doi: 10.3390/rs2020478
Njomaba, E., Ofori, J. N., Guuroh, R. T., Aikins, B. E., Nagbija, R. K., and Surový, P. (2024). Assessing forest species diversity in Ghana’s tropical forest using planetScope data. Remote Sens. 16, 463. doi: 10.3390/rs16030463
Pangtey, D., Padalia, H., Bodh, R., Rai, I. D., and Nandy, S. (2023). Application of remote sensing-based spectral variability hypothesis to improve tree diversity estimation of seasonal tropical forest considering phenological variations. Geocarto Int. 38, 2178525. doi: 10.1080/10106049,.2023.2178525
Patil, V. D., Adsul, P. B., and Deshmukh, L. S. (2007). Studies on spectral reflectance under normal and nitrogen, phosphorus and pest and disease stress condition in soybean (Glycine max L.). J. Indian Soc. Remote Sens. 35, 351–359. doi: 10.1007/BF02990790
Pennekamp, F., Pontarp, M., Tabi, A., Altermatt, F., Alther, R., Choffat, Y., et al. (2018). Biodiversity increases and decreases ecosystem stability. Nature 563, 109–112. doi: 10.1038/s41586-018-0627-8
Persson, M., Lindberg, E., and Reese, H. (2018). Tree species classification with multi-temporal sentinel-2 data. Remote Sens. 10, 1794. doi: 10.3390/rs10111794
Pouteau, R., Gillespie, T. W., and Birnbaum, P. (2018). Predicting tropical tree species richness from normalized difference vegetation index time series: the devil is perhaps not in the detail. Remote Sens. 10, 698. doi: 10.3390/rs10050698
Rahmanian, S., Nasiri, V., Amindin, A., Karami, S., Maleki, S., Pouyan, S., et al. (2023). Prediction of plant diversity using multi-seasonal remotely sensed and geodiversity data in a mountainous area. Remote Sens. 15, 387. doi: 10.3390/rs15020387
Richardson, A. D. (2018). Tracking seasonal rhythms of plants in diverse ecosystems with digital camera imagery. New Phytol. 222, 1742–1750. doi: 10.1111/nph.15591
Rocchini, D., Balkenhol, N., Carter, G. A., Foody, G. M., Gillespie, T. W., He, K. S., et al. (2010). Remotely sensed spectral heterogeneity as a proxy of species diversity: Recent advances and open challenges. Ecol. Inf. 5, 318–329. doi: 10.1016/j.ecoinf.2010.06.001
Rossi, C., Kneubühler, M., Schütz, M., Schaepman, M. E., Haller, R. M., Risch, A. C., et al. (2021). Spatial resolution, spectral metrics and biomass are key aspects in estimating plant species richness from spectral diversity in species-rich grasslands. Remote Sens. Ecol. Conserv. 8, 297–314. doi: 10.1002/rse2.244
Shannon, C. E. (1948). A mathematical theory of communication. Bell System Tech. J. 27, 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x
Tolan, J., Yang, H.-I., Nosarzewski, B., Couairon, G., Vo, H. V., Brandt, J., et al. (2024). Very high resolution canopy height maps from RGB imagery using self-supervised vision transformer and convolutional decoder trained on aerial lidar. Remote Sens. Environ. 300, 113888. doi: 10.1016/j.rse.2023.113888
Torresani, M., Feilhauer, H., Rocchini, D., Féret, J. B., Zebisch, M., and Tonon, G. (2021). Which optical traits enable an estimation of tree species diversity based on the Spectral Variation Hypothesis? Appl. Vegetation Sci. 24, e12586. doi: 10.1111/avsc.12586
Torresani, M., Rocchini, D., Sonnenschein, R., Zebisch, M., Hauffe, H. C., Heym, M., et al. (2020). Height variation hypothesis: A new approach for estimating forest species diversity with CHM LiDAR data. Ecol. Indic. 117, 106520. doi: 10.1016/j.ecolind.2020.106520
Torresani, M., Rocchini, D., Sonnenschein, R., Zebisch, M., Marcantonio, M., Ricotta, C., et al. (2019). Estimating tree species diversity from space in an alpine conifer forest: The Rao’s Q diversity index meets the spectral variation hypothesis. Ecol. Inf. 52, 26–34. doi: 10.1016/j.ecoinf.2019.04.001
Tucker, C. J. (1979). Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 8, 127–150. doi: 10.1016/0034-4257(79)90013-0
Vanguri, R., Laneve, G., and Hościło, A. (2024). Mapping forest tree species and its biodiversity using EnMAP hyperspectral data along with Sentinel-2 temporal data: An approach of tree species classification and diversity indices. Ecol. Indic. 167, 112671. doi: 10.1016/j.ecolind.2024.112671
Wang, R., Gamon, J. A., Cavender-Bares, J., Townsend, P. A., and Zygielbaum, A. I. (2018a). The spatial sensitivity of the spectral diversity–biodiversity relationship: an experimental test in a prairie grassland. Ecol. Appl. 28, 541–556. doi: 10.1002/eap.1669
Wang, R., Gamon, J. A., Schweiger, A. K., Cavender-Bares, J., Townsend, P. A., Zygielbaum, A. I., et al. (2018b). Influence of species richness, evenness, and composition on optical diversity: A simulation study. Remote Sens. Environ. 211, 218–228. doi: 10.1016/j.rse.2018.04.010
Wang, D., Qiu, P., Wan, B., Cao, Z., and Zhang, Q. (2022). Mapping α- and β-diversity of mangrove forests with multispectral and hyperspectral images. Remote Sens. Environ. 275, 113021. doi: 10.1016/j.rse.2022.113021
Wiegand, T., May, F., Kazmierczak, M., and Huth, A. (2017). What drives the spatial distribution and dynamics of local species richness in tropical forest? Proc. R. Soc. B.: Biol. Sci. 284, 20171503. doi: 10.1098/rspb.2017.1503
Wu, B., Qian, J., and Zeng, Y. (2017). Land cover atlas of the people’s republic of China (1:1,000,000) (Beijing: Sinomaps Press).
Xi, Y., Ren, C., Tian, Q., Ren, Y., Dong, X., and Zhang, Z. (2021). Exploitation of time series sentinel-2 data and different machine learning algorithms for detailed tree species classification. IEEE J. Selected Topics Appl. Earth Observ. Remote Sens. 14, 7589–7603. doi: 10.1109/jstars.2021.3098817
Xi, Y., Zhang, W., Brandt, M., Tian, Q., and Fensholt, R. (2023). Mapping tree species diversity of temperate forests using multi-temporal Sentinel-1 and -2 imagery. Sci. Remote Sens. 8, 100094. doi: 10.1016/j.srs.2023.100094
Xia, Q., Li, J., Dai, S., Zhang, H., and Xing, X. (2023). Mapping high-resolution mangrove forests in China using GF-2 imagery under the tide. Natl. Remote Sens. Bull. 27, 1320–1333. doi: 10.11834/jrs.20221848
Xu, H., Li, Y., Lin, M., Wu, J., Luo, T., Zhou, Z., et al. (2015). Community characteristics of a 60 ha dynamics plot in the tropical montane rain forest in Jianfengling, Hainan Island. Biodivers. Sci. 23, 192–201. doi: 10.17520/biods.2014157
Xu, C., Zeng, Y., Zheng, Z., Zhao, D., Liu, W., Ma, Z., et al. (2022). Assessing the impact of soil on species diversity estimation based on UAV imaging spectroscopy in a natural alpine steppe. Remote Sens. 14, 671. doi: 10.3390/rs14030671
Yang, Q., Wang, L., Huang, J., Lu, L., Li, Y., Du, Y., et al. (2022). Mapping plant diversity based on combined SENTINEL-1/2 data—Opportunities for subtropical mountainous forests. Remote Sens. 14, 492. doi: 10.3390/rs14030492
Zhai, D., Xu, J., Dai, Z., Cannon, C. H., and Grumbine, R. E. (2013). Increasing tree cover while losing diverse natural forests in tropical Hainan, China. Regional Environ. Change 14, 611–621. doi: 10.1007/s10113-013-0512-9
Zhao, Y., Zeng, Y., Zheng, Z., Dong, W., Zhao, D., Wu, B., et al. (2018). Forest species diversity mapping using airborne LiDAR and hyperspectral data in a subtropical forest in China. Remote Sens. Environ. 213, 104–114. doi: 10.1016/j.rse.2018.05.014
Zheng, Z., Zeng, Y., Schuman, M. C., Jiang, H., Schmid, B., Schaepman, M. E., et al. (2022). Individual tree-based vs pixel-based approaches to mapping forest functional traits and diversity by remote sensing. Int. J. Appl. Earth Observ. Geoinform. 114, 103074. doi: 10.1016/j.jag.2022.103074
Zhou, J., Dian, Y., Wang, X., Yao, C., Jian, Y., Li, Y., et al. (2020). Comparison of GF2 and SPOT6 imagery on canopy cover estimating in northern subtropics forest in China. Forests 11, 407. doi: 10.3390/f11040407
Keywords: tropical forest, species diversity, phenology, spatial scale, sentinel-2, GF2
Citation: Zhao P, Zeng Y, Zheng Z, Xu C, Wu J, Mu X, Zhou Z, Chen J, Zhang T and Zhao D (2025) Species diversity estimation in a typical tropical forest: which phenological stage and spatial resolution are suitable? Front. Plant Sci. 16:1582910. doi: 10.3389/fpls.2025.1582910
Received: 25 February 2025; Accepted: 25 July 2025;
Published: 18 August 2025.
Edited by:
Ram P. Sharma, Tribhuvan University, NepalReviewed by:
Juan Carlos Jimenez, University of Valencia, SpainShihua Li, University of Electronic Science and Technology of China, China
Zhihui Wang, Guangzhou Institute of Geography, China
Smrithy Vijayan, Agharkar Research Institute, India
Copyright © 2025 Zhao, Zeng, Zheng, Xu, Wu, Mu, Zhou, Chen, Zhang and Zhao. 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: Zhaoju Zheng, emhlbmd6aGFvanVAYWlyY2FzLmFjLmNu
†These authors have contributed equally to this work and share first authorship