Yinggang Zhan

Accurate estimation of fine-resolution FPAR is urgently needed for modeling land surface processes at finer scales. While traditional methods can hardly balance universality, efficiency, and accuracy, methods using coarse-resolution products as a reference are promising for operational fine-resolution FPAR estimation. However, current methods confront major problems of underrepresentation of FPAR-reflectance relations within coarse-resolution FPAR products, particularly for densely vegetated areas. To overcome this limitation, this article has developed an enhanced scaling method that proposes an outlier removal procedure and a method weighting the selected samples and models FPAR through weighted multiple linear regression (MLR) between the coarse-resolution FPAR product and the aggregated fine-resolution surface reflectance. Meanwhile, a random forest regression (RFR) method has also been implemented for comparison. Both methods were particularly applied to Landsat 8 OLI and MODIS FPAR data on the Google Earth Engine. Their performance was tested on a regional scale for an entire year. The results of the enhanced scaling method were closer to the in situ measurements (RMSE = 0.058 and R² = 0.768) and were more consistent with the MODIS FPAR (RMSE = 0.091 and R² =0.894) than those of the RFR, particularly over densely vegetated pixels. This indicates that a well-designed simple MLR-based method can outperform the more sophisticated RFR method. The enhanced scaling method is also less sensitive to the number of training samples than the RFR method. Moreover, both methods are insensitive to land cover maps, and their computation efficiency depends on the number of images to be estimated.

As a direct indicator of vegetation photosynthesis, the fraction of absorbed photosynthetically active radiation (FPAR) serves as a critical input in a series of land surface models. While existing satellite FPAR products are generally at coarse resolutions ranging from 250 m to 1 km, operational FPAR products at fine resolution are urgently needed in studying land surface processes at the plot scale. However, existing methods for estimating fine-resolution FPAR were mainly designed for Landsat data, and few studies have attempted to develop algorithms for Sentinel-2 data. In particular, the operational estimation of decameter FPAR has a higher requirement for the algorithms in terms of generalizability, efficiency, accuracy, and adaptability to Sentinel-2 data. In this paper, we developed a retrieval chain on the Google Earth Engine (GEE) platform to estimate FPAR by learning the relationship between MODIS FPAR and Sentinel-2 surface reflectance. Scale-consistent multilinear models were used to model the relationship between MODIS FPAR and Sentinel-2 surface reflectance, and the model coefficients were regressed from the selected training samples. To account for the spectral and spatial characteristics of the Sentinel-2 data, we designed criteria for selecting training samples and compared different band combinations. Three strategies for band combination were used: (1) green, red, and near infrared (NIR) bands at 10 m resolution (i.e., three bands); (2) green, red, NIR, and red edge (RE) 1, RE2, and RE3 bands at 20 m resolution (i.e., five bands); and (3) green, red, NIR, RE1, RE2, RE3, shortwave infrared1 (SWIR1) and SWIR2 bands at 20 m resolution (i.e., eight bands). Meanwhile, the official Sentinel Application Platform (SNAP) method has also been implemented to estimate the Sentinel FPAR at 10 m and 20 m resolutions for comparison. Both methods were applied to the western Guanzhong area, Shaanxi Province, China, for FPAR estimation of all cloud-free Sentinel-2 images in 2021. The results show that the scaling-based method using five bands at 20 m resolution was the most accurate compared to the in situ measurements (RMSE = 0.076 and R² = 0.626), which outperformed the SNAP method at 10 m and 20 m resolutions and the scaling-based method using other strategies. The results of the scaling-based method using all three strategies were highly consistent with the MODIS FPAR product, while the SNAP method systematically underestimated FPAR values compared to the MODIS FPAR products. The proposed method is more ready-to-use and more efficient than SNAP software. Considering that the service of the MODIS sensor is overdue, the proposed method can be extended to alternatives to MODIS products, such as VIIRS and Sentinel-3 data.

The accurate estimation of phenological metrics from satellite data, especially the start of season (SOS), is of great significance to enhance our understanding of trends in vegetation phenology under climate change at regional or global scales. However, for regions with winter snow cover, such as the alpine grasslands on the Tibetan Plateau, the presence of snow inevitably contaminates satellite signals and introduces bias into the detection of the SOS. Despite recent progress in eliminating the effect of snow cover on SOS detection, the mechanism of how snow cover affects the satellite-derived vegetation index (VI) and the detected SOS remains unclear. This study investigated the effect of snow cover on both VI and SOS detection by combining simulation experiments and real satellite data. Five different VIs were used and compared in this study, including four structure-based (i.e., NDVI, EVI2, NDPI, NDGI) VIs and one physiological-based (i.e., NIRv) VI. Both simulation experiments and satellite data analysis revealed that the presence of snow can significantly reduce the VI values and increase the local gradient of the growth curve, allowing the SOS to be detected. The bias in the detected SOS caused by snow cover depends on the end of the snow season (ESS), snow duration parameters, and the snow-free SOS. An earlier ESS results in an earlier estimate of the SOS, a later ESS results in a later estimate of the SOS, and an ESS close to the snow-free SOS results in small bias in the detected SOS. The sensitivity of the five VIs to snow cover in SOS detection is NDPI/NDGI < NIRv < EVI2 < NDVI, which has been verified in both simulation experiments and satellite data analysis. These findings will significantly advance our research on the feedback mechanisms between vegetation, snow, and climate change for alpine ecosystems.

Despite its wide applications, the spatiotemporal fusion of coarse- and fine-resolution satellite images is limited primarily to the availability of clear-sky fine-resolution images, which are commonly scarce due to unfavorable weather, and such a limitation might cause errors in spatiotemporal fusion. Thus, the effective use of limited fine-resolution images, while critical, remains challenging. To address this issue, in this paper we propose a new phenological similarity strategy (PSS) to select the optimal combination of image pairs for a prediction date. The PSS considers the temporal proximity and phenological similarity between the base and prediction images and computes a weight for identifying the optimal combination of image pairs. Using the PSS, we further evaluate the influence of input data on the fusion accuracy by varying the number and temporal distribution of input images. The results show that the PSS (mean R = 0.827 and 0.760) outperforms the nearest date (mean R = 0.786 and 0.742) and highest correlation (mean R = 0.821 and 0.727) strategies in both the enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM) and the linear mixing growth model (LMGM), respectively, for fusing Landsat 8 OLI and MODIS NDVI datasets. Furthermore, base images adequately covering different growth stages yield better predictability than simply increasing the number of base images.