Accurate and updated finer resolution maps of rubber plantations and stand ages are needed to understand and assess the impacts of rubber plantations on regional ecosystem processes. This study presented a simple method for mapping rubber plantation areas and their stand ages by integration of PALSAR 50-m mosaic images and multi-temporal Landsat TM/ETM+ images. The L-band PALSAR 50-m mosaic images were used to map forests (including both natural forests and rubber trees) and non-forests. For those PALSAR-based forest pixels, we analyzed the multi-temporal Landsat TM/ETM+ images from 2000 to 2009. We first studied phenological signatures of deciduous rubber plantations (defoliation and foliation) and natural forests through analysis of surface reflectance, Normal Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), and Land Surface Water Index (LSWI) and generated a map of rubber plantations in 2009. We then analyzed phenological signatures of rubber plantations with different stand ages and generated a map, in 2009, of rubber plantation stand ages (≤5, 6–10, >10 years-old) based on multi-temporal Landsat images. The resultant maps clearly illustrated how rubber plantations have expanded into the mountains in the study area over the years. The results in this study demonstrate the potential of integrating microwave (e.g., PALSAR) and optical remote sensing in the characterization of rubber plantations and their expansion over time.

Knowing the stand age of rubber tree (Hevea brasiliensis) plantations is vitally important for best management practices, estimations of rubber latex yields, and carbon cycle studies (e.g., biomass, carbon pools, and fluxes). However, the stand age (as estimated from the establishment year of rubber plantation) is not available across large regions. In this study, we analyzed Landsat time series images from 1987–2015 and developed algorithms to identify (1) the establishment year of rubber plantations; and (2) the pre-conversion land cover types, such as old rubber plantations, evergreen forests, and cropland. Exposed soil during plantation establishment and linear increases in canopy closure during non-production periods (rubber seedling to mature plantation) were used to identify the establishment year of rubber plantations. Based on the rubber plantation map for 2015 (overall accuracy = 97%), and 1981 Landsat images since 1987, we mapped the establishment year of rubber plantations on Hainan Island (R² = 0.85/0.99, and RMSE = 2.34/0.54 years at pixel/plantation scale). The results show that: (1) significant conversion of croplands and old rubber plantations to new rubber plantations has occurred substantially in the northwest and northern regions of Hainan Island since 2000, while old rubber plantations were mainly distributed in the southeastern inland strip; (2) the pattern of rubber plantation expansion since 1987 consisted of fragmented plantations from smallholders, and there was no tendency to expand towards a higher altitude and steep slope regions; (3) the largest land source for new rubber plantations since 1988 was old rubber plantations (1.26 × 10⁵ ha), followed by cropland (0.95 × 10⁵ ha), and evergreen forests (0.68 × 10⁵ ha). The resultant algorithms and maps of establishment year and pre-conversion land cover types are likely to be useful in plantation management, and ecological assessments of rubber plantation expansion in China.

The increasing expansion of rubber plantations throughout East and Southeast Asia urgently requires improved methods for effective mapping and monitoring. The phenological information from rubber plantations was found effective in rubber mapping. Previous studies have mostly applied rule-pixel-based phenology approaches for rubber plantations mapping, which might result in broken patches in fragmented landscapes. This study introduces a new paradigm by combining phenology information with object-based classification to map fragmented patches of rubber plantations in Xishuangbanna. This research first delineated the time windows of the defoliation and foliation phases of rubber plantations by acquiring the temporal profile and phenological features of rubber plantations and natural forests through the Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) data. To investigate the ability of finer resolution images at capturing the temporal profile or phenological information, 30 m resolution Landsat image data were used to capture the temporal profile, and a phenology algorithm to separate rubber plantations and natural forests was then defined. The derived phenology algorithm was used by both the object-based and pixel-based classification to investigate whether the object-based approach could improve the mapping accuracy. Whether adding the phenology information to the object-based classification could improve rubber plantation mapping accuracy in mountainous Xishuangbanna was also investigated. This resulted in three approaches: rule-pixel-based phenology, rule-object-based phenology, and nearest-neighbor-object-based phenology. The results showed that the rule-object-based phenology approaches (with overall accuracy 77.5% and Kappa Coefficients of 0.66) and nearest-neighbor-object-based phenology approach (91.0% and 0.86) achieved a higher accuracy than that of the rule-pixel-based phenology approach (72.7% and 0.59). The results proved that (1) object-based approaches could improve the accuracy of rubber plantation mapping compared to the pixel-based approach and (2) incorporating the phenological information from vegetation improved the overall accuracy of the thematic map.

Updated extent, area, and spatial distribution of tropical evergreen forests from inventory data provides valuable knowledge for research of the carbon cycle, biodiversity, and ecosystem services in tropical regions. However, acquiring these data in mountainous regions requires labor-intensive, often cost-prohibitive field protocols. Here, we report about validated methods to rapidly identify the spatial distribution of tropical forests, and obtain accurate extent estimates using phenology-based procedures that integrate the Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat imagery. Firstly, an analysis of temporal profiles of annual time-series MODIS Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), and Land Surface Water Index (LSWI) was developed to identify the key phenology phase for extraction of tropical evergreen forests in five typical lands cover types. Secondly, identification signatures of tropical evergreen forests were selected and their related thresholds were calculated based on Landsat NDVI, EVI, and LSWI extracted from ground true samples of different land cover types during the key phenology phase. Finally, a map of tropical evergreen forests was created by a pixel-based thresholding. The developed methods were tested in Xishuangbanna, China, and the results show: (1) Integration of Landsat and MODIS images performs well in extracting evergreen forests in tropical complex mountainous regions. The overall accuracy of the resulting map of the case study was 92%; (2) Annual time series of high-temporal-resolution remote sensing images (MODIS) can effectively be used for identification of the key phenology phase (between Julian Date 20 and 120) to extract tropical evergreen forested areas through analysis of NDVI, EVI, and LSWI of different land cover types; (3) NDVI and LSWI are two effective metrics (NDVI ≥ 0.670 and 0.447 ≥ LSWI ≥ 0.222) to depict evergreen forests from other land cover types during the key phenology phase in tropical complex mountainous regions. This method can make full use of the Landsat and MODIS archives as well as their advantages for tropical evergreen forests geospatial inventories, and is simple and easy to use. This method is suggested for use with other similar regions.

Rubber (Hevea brasiliensis Muell.) plantations constitute one of the most important agro-ecosystems in the tropical region of China and Southeast Asia, playing an important role in the carbon budget there. Accurately obtaining their biomass over a large area is challenging because of difficulties in acquiring the Diameter at Breast Height (DBH) through remote sensing and the problem of biomass saturation. The stand age, which is closely related to the forest biomass, was proposed for biomass estimation in this study. A stand age map at an annual scale for Hainan Island, which is the second largest natural rubber production base in China, was generated using all Landsat and Sentinel-2 (LS2) data (1987–2017). Scatter plots and the correlation coefficient method were used to explore the relationship (e.g., biomass saturation) between rubber biomass and different LS2-based variables. Subsequently, a regression model fitted with the stand age (R² = 0.96) and a Random Forest (RF) model parameterizing with LS2-based variables and/or the stand age were respectively employed to estimate rubber biomass for Hainan Island. The results show that rubber biomass was saturated around 65 Mg/ha with all LS2-based variables. The regression model estimated biomass accurately (R² = 0.79 and Root Mean Square Error (RMSE) = 14.00 Mg/ha) and eliminated the saturation problem significantly. In addition to LS2-based variables, adding a stand age parameter to the RF models was found to significantly improve the prediction accuracy (R² = 0.82–0.96 and RMSE = 4.08–10.59 Mg/ha, modeling using samples of different biomass sizes). However, all RF models overestimated the biomass of young plantations and underestimated the biomass of old plantations. A hybrid model integrating the optimal results of RF and regression models reduced estimation bias and generated the best performance (R² = 0.83 and RMSE = 12.48 Mg/ha). The total rubber biomass of Hainan Island in 2017 was about 5.40 × 10⁷ Mg. The northward and westward expansions after 2000 had great impact on the biomass distribution, leading to a higher biomass density for the inland coastal strip from south to northeast and a lower biomass density in the northern and western regions.

The ongoing global warming trajectory poses extensive challenges to plant ecosystems, with rubber plantations particularly vulnerable due to their influence on not only the longevity of the growth cycle and rubber yield, but also the complex interplay of carbon, water, and energy exchanges between the forest canopy and atmosphere. However, the response mechanism of phenology in rubber plantations to climate change remains unclear. This study concentrates on sub-optimal environment rubber plantations in Yunnan province, Southwest China. Utilizing the Google Earth Engine (GEE) cloud platform, multi-source remote sensing images were synthesized at 8-day intervals with a spatial resolution of 30-meters. The Normalized Difference Vegetation Index (NDVI) time series was reconstructed using the Savitzky-Golay (S-G) filter, coupled with the application of the seasonal amplitude method to extract three crucial phenological indicators, namely the start of the growing season (SOS), the end of the growing season (EOS), and the length of the growing season (LOS). Linear regression method, Pearson correlation coefficient, multiple stepwise regression analysis were used to extract of the phenology trend and find the relationship between SOS, EOS and climate factors. The findings demonstrated that 1) the phenology of rubber plantations has undergone dynamic changes over the past two decades. Specifically, the SOS advanced by 9.4 days per decade (R² = 0.42, p< 0.01), whereas the EOS was delayed by 3.8 days per decade (R² = 0.35, p< 0.01). Additionally, the LOS was extended by 13.2 days per decade (R² = 0.55, p< 0.01); 2) rubber phenology demonstrated a notable sensitivity to temperature fluctuations during the dry season and precipitation patterns during the rainy season. The SOS advanced 2.0 days (r =−0.19, p< 0.01) and the EOS advanced 2.8 days (r =−0.35, p< 0.01) for every 1°C increase in the cool-dry season. Whereas a 100 mm increase in rainy season precipitation caused the SOS to be delayed by 2.0 days (r = 0.24, p< 0.01), a 100 mm increase in hot-dry season precipitation caused the EOS to be advanced by 7.0 days (r =-0.28, p< 0.01); 3) rubber phenology displayed a legacy effect of preseason climate variations. Changes in temperature during the fourth preseason month and precipitation during the fourth and eleventh preseason months are predominantly responsible for the variation in SOS. Meanwhile, temperature changes during the second, fourth, and ninth preseason months are primarily responsible for the variation in EOS. The study aims to enhance our understanding of how rubber plantations respond to climate change in sub-optimal environments and provide valuable insights for sustainable rubber production management in the face of changing environmental conditions.