Scale effect, which is caused by a combination of model nonlinearity and surface heterogeneity, has been of interest to the remote sensing community for decades. However, there is no current analysis of scale effect in the ground-based indirect measurement of leaf area index (LAI), where model nonlinearity and surface heterogeneity also exist. This paper examines the scale effect on the indirect measurement of LAI. We built multiscale data sets based on realistic scenes and field measurements. We then implemented five representative methods of indirect LAI measurement at scales (segment lengths) that range from meters to hundreds of meters. The results show varying degrees of deviation and fluctuation that exist in all five methods when the segment length is shorter than 20 m. The retrieved LAI from either Beer's law or the gap-size distribution method shows a decreasing trend with increasing segment lengths. The length at which the LAI values begin to stabilize is about a full period of row in row crops and 100 m in broadleaf or coniferous forests. The impacts of segment length on the finite-length averaging method, the combination of gap-size distribution and finite-length methods, and the path-length distribution method are relatively small. These three methods stabilize at the segment scale longer than 20 m in all scenes. We also find that computing the average LAI of all of the short segment lengths, which is commonly done, is not as good as merging these short segments into a longer one and computing the LAI value of the merged one.

The fraction of absorbed photosynthetically active radiation (FAPAR) is a critical biophysical parameter in eco-environmental studies. Scaling of FAPAR from the field observation to the satellite pixel is essential for validating remote sensing FAPAR product and for further modeling applications. However, compared to spatial mismatches, few studies have considered temporal mismatches between in-situ and satellite observations in the scaling. This paper proposed a general methodology for scaling FAPAR from the field to the satellite pixel considering the temporal variation. Firstly, a temporal normalization method was proposed to normalize the in-situ data measured at different times to the time of satellite overpass. The method was derived from the integration of an atmospheric radiative transfer model (6S) and a FAPAR analytical model (FAPAR-P), which can characterize the diurnal variations of FAPAR comprehensively. Secondly, the logistic model, which derives smooth and consistent temporal profile for vegetation growth, was used to interpolate the in-situ data to match the dates of satellite acquisitions. Thirdly, fine-resolution FAPAR products at different dates were estimated from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data using the temporally corrected in-situ data. Finally, fine-resolution FAPAR were taken as reference datasets and aggregated to coarse resolution, which were further compared to the Moderate Resolution Imaging Spectroradiometer (MODIS) FAPAR product. The methodology is validated for scaling FAPAR from the field to the satellite pixel temporally and spatially. The MODIS FAPAR manifested a good consistency with the aggregated FAPAR with R² of 0.922 and the root mean squared error of 0.054.

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.

Accurate estimation of the fine-resolution fraction of absorbed photosynthetically active radiation (FPAR) across broad spatial extents and long time periods requires efficient and applicable methods. The existing methods can hardly provide a balance between accuracy, simplicity, and transferability through space and time. Within the remote-sensing trend-surface conceptual framework, this article proposes a scaling-based method to efficiently retrieve FPAR from fine-resolution satellite data using coarse-resolution FPAR products as a reference. The method was particularly developed and applied to Moderate Resolution Imaging Spectroradiometer (MODIS) FPAR product and Landsat imagery. First, necessary prior knowledge related to FPAR retrieval and scaling theories was used to explicitly linearize the complex relationship between MODIS FPAR and Landsat surface reflectance. Second, the explicit linear model for FPAR estimation was trained through one-pair image learning for each date to estimate FPAR from Landsat imagery in real time. Both homogeneous and heterogeneous cases were considered. The method was validated at ten selected worldwide sites from the Validation of Land European Remote Sensing Instruments (VALERI) program and derived an overall root mean squared error (RMSE) of 0.133. A long time series of FPAR data set at the 30-m resolution was generated at the regional scale (approximately 2000 km²) for 13 years (2000-2012). The results were accurate (RMSE = 0.072) and MODIS-consistent, which were significantly better than those of the normalized difference vegetation index (NDVI) downscaling-based and regression tree methods. The scaling-based method provides accurate, MODIS-consistent and spatially consistent FPAR estimates in real time, is highly transferrable through space and time, and allows for future extension of FPAR estimates to the era of the Landsat series satellites.

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.