The HANTS-fitted RSEI constructed in the vegetation growing season reveals the spatiotemporal patterns of ecological quality.

Yuxi, located in China's central plateau of Yunnan, is grappling with ecological and environmental challenges as it continues to develop its economy. While ecological quality assessment serves as the foundation for ecological protection, it is pivotal to have reliable and long-term methods for assessing the ecological status to support informed decision-making in ecological protection. Reliable and long-term methods for assessing ecological status in order to facilitate informed decision-making in ecological protection are applied. This study utilized Landsat data to reconstruct four indices (greenness, wetness, dryness, and heat) during the vegetation growth in Yuxi from 2000 to 2020 that employs Harmonic Analysis of Time Series (HANTS) method. Subsequently, the annual Remote Sensing Ecological Index (RSEI) was computed by using the reconstructed indices to evaluate ecological quality in Yuxi. Additionally, spatiotemporal patterns and determinants of Yuxi's ecological quality are unveiled through Sen's slope estimator and Mann-Kendall test (Sen + MK) trend analysis, spatial auto-correlation analysis, and geographical detectors applied to year-by-year RSEI data. The findings in the paper indicate that the accuracy of the RSEI is significantly influenced by the vegetation season, suggesting that constructing the RSEI model with data from the vegetation growth season is crucial. Moreover, the HANTS optimization method effectively enhances the ecological indices used in the RSEI model, leading to smoother and more continuous filling of missing data. The difference between the reconstructed RSEI and the original RSEI falls within the range of - 0.15 to 0.15. Yuxi has an average RSEI of 0.54 to emphasis a moderate level of comprehensive ecological quality. Compared with river valley plains, the ecological quality of mountainous areas is higher, and the ecological quality of Yuxi presents a distinct center-edge pattern. From 2000 to 2020, Yuxi's ecological quality exhibited fluctuations, with a slight overall improvement. Land use patterns, particularly in forestry land and impervious surfaces, are identified as the main drivers of these changes. The research offers valuable insights for scientific decision-making related to sustainable development and ecological protection.

Dry season temperature and rainy season precipitation significantly affect the spatio-temporal pattern of rubber plantation phenology in Yunnan province.

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 (R2 = 0.42, p< 0.01), whereas the EOS was delayed by 3.8 days per decade (R2 = 0.35, p< 0.01). Additionally, the LOS was extended by 13.2 days per decade (R2 = 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.

Diversified land conversion deepens understanding of impacts of rapid rubber plantation expansion on plant diversity in the tropics.

Understanding the status and changes of plant diversity in rubber (Hevea brasiliensis) plantations is essential for sustainable plantation management in the context of rapid rubber expansion in the tropics, but remains very limited at the continental scale. In this study, we investigated plant diversity from 10-meter quadrats in 240 different rubber plantations in the six countries of the Great Mekong Subregion (GMS)-where nearly half of the world's rubber plantations are located-and analyzed the influence of original land cover types and stand age on plant diversity using Landsat and Sentinel-2 satellite imagery since the late 1980s. The results indicate that the average plant species richness of rubber plantations is 28.69 ± 7.35 (1061 species in total, of which 11.22 % are invasive), approximating half the species richness of tropical forests but roughly double that of the intensively managed croplands. Time-series satellite imagery analysis revealed that rubber plantations were primarily established in place of cropland (RPC, 37.72 %), old rubber plantations (RPORP, 27.63 %), and tropical forests (RPTF, 24.12 %). Plant species richness in RPTF (34.02 ± 7.62) was significantly (p < 0.001) higher than that in RPORP (26.41 ± 7.02) and RPC (26.34 ± 5.37). More importantly, species richness can be maintained for the duration of the 30-year economic cycle, and the number of invasive species decreases as the stand ages. Given diverse land conversions and changes in stand age, the total loss of species richness due to rapid rubber expansion in the GMS was 7.29 %, which is far below the traditional estimates that only consider tropical forest conversion. In general, maintaining higher species richness at the earliest stages of cultivation has significant implications for biodiversity conservation in rubber plantations.

Mapping paddy rice planting area in rice-wetland coexistent areas through analysis of Landsat 8 OLI and MODIS images.

Accurate and up-to-date information on the spatial distribution of paddy rice fields is necessary for the studies of trace gas emissions, water source management, and food security. The phenology-based paddy rice mapping algorithm, which identifies the unique flooding stage of paddy rice, has been widely used. However, identification and mapping of paddy rice in rice-wetland coexistent areas is still a challenging task. In this study, we found that the flooding/transplanting periods of paddy rice and natural wetlands were different. The natural wetlands flood earlier and have a shorter duration than paddy rice in the Panjin Plain, a temperate region in China. We used this asynchronous flooding stage to extract the paddy rice planting area from the rice-wetland coexistent area. MODIS Land Surface Temperature (LST) data was used to derive the temperature-defined plant growing season. Landsat 8 OLI imagery was used to detect the flooding signal and then paddy rice was extracted using the difference in flooding stages between paddy rice and natural wetlands. The resultant paddy rice map was evaluated with in-situ ground-truth data and Google Earth images. The estimated overall accuracy and Kappa coefficient were 95% and 0.90, respectively. The spatial pattern of OLI-derived paddy rice map agrees well with the paddy rice layer from the National Land Cover Dataset from 2010 (NLCD-2010). The differences between RiceLandsat and RiceNLCD are in the range of ±20% for most 1-km grid cell. The results of this study demonstrate the potential of the phenology-based paddy rice mapping algorithm, via integrating MODIS and Landsat 8 OLI images, to map paddy rice fields in complex landscapes of paddy rice and natural wetland in the temperate region.

Mapping paddy rice planting area in wheat-rice double-cropped areas through integration of Landsat-8 OLI, MODIS, and PALSAR images.

As farmland systems vary over space and time (season and year), accurate and updated maps of paddy rice are needed for studies of food security and environmental problems. We selected a wheat-rice double-cropped area from fragmented landscapes along the rural-urban complex (Jiangsu Province, China) and explored the potential utility of integrating time series optical images (Landsat-8, MODIS) and radar images (PALSAR) in mapping paddy rice planting areas. We first identified several main types of non-cropland land cover and then identified paddy rice fields by selecting pixels that were inundated only during paddy rice flooding periods. These key temporal windows were determined based on MODIS Land Surface Temperature and vegetation indices. The resultant paddy rice map was evaluated using regions of interest (ROIs) drawn from multiple high-resolution images, Google Earth, and in-situ cropland photos. The estimated overall accuracy and Kappa coefficient were 89.8% and 0.79, respectively. In comparison with the National Land Cover Data (China) from 2010, the resultant map better detected changes in the paddy rice fields and revealed more details about their distribution. These results demonstrate the efficacy of using images from multiple sources to generate paddy rice maps for two-crop rotation systems.

Mapping paddy rice planting areas through time series analysis of MODIS land surface temperature and vegetation index data.

Knowledge of the area and spatial distribution of paddy rice is important for assessment of food security, management of water resources, and estimation of greenhouse gas (methane) emissions. Paddy rice agriculture has expanded rapidly in northeastern China in the last decade, but there are no updated maps of paddy rice fields in the region. Existing algorithms for identifying paddy rice fields are based on the unique physical features of paddy rice during the flooding and transplanting phases and use vegetation indices that are sensitive to the dynamics of the canopy and surface water content. However, the flooding phenomena in high latitude area could also be from spring snowmelt flooding. We used land surface temperature (LST) data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor to determine the temporal window of flooding and rice transplantation over a year to improve the existing phenology-based approach. Other land cover types (e.g., evergreen vegetation, permanent water bodies, and sparse vegetation) with potential influences on paddy rice identification were removed (masked out) due to their different temporal profiles. The accuracy assessment using high-resolution images showed that the resultant MODIS-derived paddy rice map of northeastern China in 2010 had a high accuracy (producer and user accuracies of 92% and 96%, respectively). The MODIS-based map also had a comparable accuracy to the 2010 Landsat-based National Land Cover Dataset (NLCD) of China in terms of both area and spatial pattern. This study demonstrated that our improved algorithm by using both thermal and optical MODIS data, provides a robust, simple and automated approach to identify and map paddy rice fields in temperate and cold temperate zones, the northern frontier of rice planting.

Mapping paddy rice planting area in cold temperate climate region through analysis of time series Landsat 8 (OLI), Landsat 7 (ETM+) and MODIS imagery.

Accurate and timely rice paddy field maps with a fine spatial resolution would greatly improve our understanding of the effects of paddy rice agriculture on greenhouse gases emissions, food and water security, and human health. Rice paddy field maps were developed using optical images with high temporal resolution and coarse spatial resolution (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)) or low temporal resolution and high spatial resolution (e.g., Landsat TM/ETM+). In the past, the accuracy and efficiency for rice paddy field mapping at fine spatial resolutions were limited by the poor data availability and image-based algorithms. In this paper, time series MODIS and Landsat ETM+/OLI images, and the pixel- and phenology-based algorithm are used to map paddy rice planting area. The unique physical features of rice paddy fields during the flooding/open-canopy period are captured with the dynamics of vegetation indices, which are then used to identify rice paddy fields. The algorithm is tested in the Sanjiang Plain (path/row 114/27) in China in 2013. The overall accuracy of the resulted map of paddy rice planting area generated by both Landsat ETM+ and OLI is 97.3%, when evaluated with areas of interest (AOIs) derived from geo-referenced field photos. The paddy rice planting area map also agrees reasonably well with the official statistics at the level of state farms (R2 = 0.94). These results demonstrate that the combination of fine spatial resolution images and the phenology-based algorithm can provide a simple, robust, and automated approach to map the distribution of paddy rice agriculture in a year.