The spatial distribution of plant diversity and biomass informs management decisions to maintain biodiversity and carbon stocks in tropical forests. Optical remotely sensed data is often used for supporting such activities; however, it is difficult to estimate these variables in areas of high biomass. New technologies, such as airborne LiDAR, have been used to overcome such limitations. LiDAR has been increasingly used to map carbon stocks in tropical forests, but has rarely been used to estimate plant species diversity. In this study, we first evaluated the effect of using different plot sizes and plot designs on improving the prediction accuracy of species richness and biomass from LiDAR metrics using multiple linear regression. Second, we developed a general model to predict species richness and biomass from LiDAR metrics for two different types of tropical dry forest using regression analysis. Third, we evaluated the relative roles of vegetation structure and habitat heterogeneity in explaining the observed patterns of biodiversity and biomass, using variation partition analysis and LiDAR metrics. The results showed that with increasing plot size, there is an increase of the accuracy of biomass estimations. In contrast, for species richness, the inclusion of different habitat conditions (cluster of four plots over an area of 1.0 ha) provides better estimations. We also show that models of plant diversity and biomass can be derived from small footprint LiDAR at both local and regional scales. Finally, we found that a large portion of the variation in species richness can be exclusively attributed to habitat heterogeneity, while biomass was mainly explained by vegetation structure.

Accurate estimates of above ground biomass (AGB) are needed for monitoring carbon in tropical forests. LiDAR data can provide precise AGB estimations because it can capture the horizontal and vertical structure of vegetation. However, the accuracy of AGB estimations from LiDAR is affected by a co-registration error between LiDAR data and field plots resulting in spatial discrepancies between LiDAR and field plot data. Here, we evaluated the impacts of plot location error and plot size on the accuracy of AGB estimations predicted from LiDAR data in two types of tropical dry forests in Yucatán, México. We sampled woody plants of three size classes in 29 nested plots (80 m², 400 m² and 1000 m²) in a semi-deciduous forest (Kiuic) and 28 plots in a semi-evergreen forest (FCP) and estimated AGB using local allometric equations. We calculated several LiDAR metrics from airborne data and used a Monte Carlo simulation approach to assess the influence of plot location errors (2 to 10 m) and plot size on ABG estimations from LiDAR using regression analysis. Our results showed that the precision of AGB estimations improved as plot size increased from 80 m² to 1000 m² (R² = 0.33 to 0.75 and 0.23 to 0.67 for Kiuic and FCP respectively). We also found that increasing GPS location errors resulted in higher AGB estimation errors, especially in the smallest sample plots. In contrast, the largest plots showed consistently lower estimation errors that varied little with plot location error. We conclude that larger plots are less affected by co-registration error and vegetation conditions, highlighting the importance of selecting an appropriate plot size for field forest inventories used for estimating biomass.

In tropical dry forests, deciduousness (i.e., leaf shedding during the dry season) is an important adaptation of plants to cope with water limitation, which helps trees adjust to seasonal drought. Deciduousness is also a critical factor determining the timing and duration of carbon fixation rates, and affecting energy, water, and carbon balance. Therefore, quantifying deciduousness is vital to understand important ecosystem processes in tropical dry forests. The aim of this study was to map tree species deciduousness in three types of tropical dry forests along a precipitation gradient in the Yucatan Peninsula using Sentinel-2 imagery. We propose an approach that combines reflectance of visible and near-infrared bands, normalized difference vegetation index (NDVI), spectral unmixing deciduous fraction, and several texture metrics to estimate the spatial distribution of tree species deciduousness. Deciduousness in the study area was highly variable and decreased along the precipitation gradient, while the spatial variation in deciduousness among sites followed an inverse pattern, ranging from 91.5 to 43.3% and from 3.4 to 9.4% respectively from the northwest to the southeast of the peninsula. Most of the variation in deciduousness was predicted jointly by spectral variables and texture metrics, but texture metrics had a higher exclusive contribution. Moreover, including texture metrics as independent variables increased the variance of deciduousness explained by the models from R² = 0.56 to R² = 0.60 and the root mean square error (RMSE) was reduced from 16.9% to 16.2%. We present the first spatially continuous deciduousness map of the three most important vegetation types in the Yucatan Peninsula using high-resolution imagery.

Integrating information about the spatial distribution of carbon stocks and species diversity in tropical forests over large areas is fundamental for climate change mitigation and biodiversity conservation. In this study, spatial models showing the distribution of carbon stocks and the number of species were produced in order to identify areas that maximize carbon storage and biodiversity in the tropical forests of the Yucatan Peninsula, Mexico. We mapped carbon density and species richness of trees using L-band radar backscatter data as well as radar texture metrics, climatic and field data with the random forest regression algorithm. We reduced sources of errors in plot data of the national forest inventory by using correction factors to account for carbon stocks of small trees (<7.5 cm DBH) and for the temporal difference between field data collection and imagery acquisition. We created bivariate maps to assess the spatial relationship between carbon stocks and diversity. Model validation of the regional maps obtained herein using an independent data set of plots resulted in a coefficient of determination (R²) of 0.28 and 0.31 and a relative mean square error of 38.5% and 33.0% for aboveground biomass and species richness, respectively, at pixel level. Estimates of carbon density were influenced mostly by radar backscatter and climatic data, while those of species richness were influenced mostly by radar texture and climatic variables. Correlation between carbon density and species richness was positive in 79.3% of the peninsula, while bivariate maps showed that 39.6% of the area in the peninsula had high carbon stocks and species richness. Our results highlight the importance of combining carbon and diversity maps to identify areas that are critical—both for maintaining carbon stocks and for conserving biodiversity.