Recently, solar-induced chlorophyll fluorescence (SIF) from satellites has shown potential for evaluating vegetation status and stress responses. Fluorescence quantum yield (ΦF) is essentially linked to vegetation stress. However, the complex physiological and structural responses of SIF and ΦFto drought need further study. This study normalized SIF as SIFnto account for angular variations and fluctuations in photosynthetically active radiation (PAR), aiming for more accurate drought monitoring. SIFnanomalies were compared to historical baselines (2019–2021 averages) of vegetation indices (VIs), raw SIF, and ΦFduring a 2019 drought in the North China Plain (NCP). The results show SIFnprovides an effective method for drought monitoring, showing the earliest decline compared to raw SIF, VIs, and ΦF. In the first two weeks of drought, SIFndecreased by 8.2%, 7.0%, 12.5%, and 8.2% across the four NCP subdivisions. SIFnoutperformed other indicators, proving sensitive to early drought detection. SIFnwas also examined for tracking drought alleviation by rainfall. The uncertainty under different viewing geometries was quantified. SIFnanomalies showed a strong correlation with rainfall anomalies (R: 0.45 ~ 0.52) and meteorological factors like PAR (R: 0.80 ~ 0.84) and relative humidity (R:0.52 ~ 0.54). The correlation of near-infrared reflectance (NIRv) and ΦFanomalies with SIF was weak during drought onset (R: 0.16 ~ 0.32) but strong at the end (R: 0.83 ~ 0.87). These suggest both canopy structure (mainly characterized by NIRv) and vegetation chlorophyll (ΦF) are impacted by drought and influence SIF at different stages.

Solar-induced chlorophyll fluorescence (SIF) is an effective indicator to track the gross primary productivity (GPP). However, there is still a lack of a clear understanding for the contribution of physiological and structural factors to the SIF-GPP relationship at the canopy scale. To quantify the influence of different SIF components, particularly the photon escape ratio (${f}_{\text {esc}}$) and fluorescence yield (${\Phi }_{F}$), on the SIF-GPP relationship, this study evaluated the performance of various vegetation indices (VIs), SIF components, light use efficiency (LUE), and GPP over two typical biomes (crops and rainforest), using a range of satellite remote sensing products. In August of each year from 2018 to 2020, both SIF and GPP over United States (U.S.) Corn Belt are higher than those over the Amazon rainforest, attributed to the consistent pattern of higher$f_{\text {esc}}$and LUE over crops than over rainforest, as well as${\Phi }_{F}$. Furthermore, the structural signals represented by${f}_{\text {esc}}$($R =0.41$–0.64) can better capture the LUE variations than$\boldsymbol {\Phi }_{F}$($R =0.10$–0.30) for each biome. This study highlights that$f_{\text {esc}}$, determined by canopy structure, has great potential to capture LUE and GPP changes within and across biomes.

The diurnal variation of land surface albedo (DVLSA) is crucial for understanding energy budgets and climate change. As topography complicates the radiative transfer processes, the estimation of DVLSA over rugged terrain becomes challenging. In this study, the topography-coupled DVLSA model (DVLSA_T) is developed to estimate DVLSA over rugged terrain. DVLSA_T represents DVLSA as a multiplication between the basic albedo and a diurnal variation factor. The basic albedo is the albedo at local noon with topographic effects removed, while the diurnal variation factor extends the albedo from local noon to different times of the day, accounting for topographic effects. Specifically, the diurnal variation factor of black-sky albedo (BSA) changes with the illumination geometry, integrating the topographic effects and U-shaped pattern of DVLSA. In contrast, the diurnal variation factor of white-sky albedo (WSA) is independent of illumination geometry and is solely influenced by topography. DVLSA_T shows good performance when compared with the 3-D radiative transfer simulations by the large-scale remote sensing data and image simulation framework (LESS) (BSA: coefficient of determination (${R}^{2}$) = 0.977; root-mean-square (RMSE) = 0.013; WSA:${R}^{2} =0.982$; and RMSE = 0.012) and sandbox measurements (blue-sky albedo:${R}^{2} = 0.904$and RMSE = 0.012). DVLSA_T also has a good agreement with in situ measurements, with an RMSE of 0.024 and an${R}^{2}$of 0.738. Our results demonstrate that DVLSA_T can effectively characterize DVLSA over rugged terrain.

Downward shortwave radiation (DSR) is critical to many surface processes, and many satellite-derived DSR products have been released. Few studies have validated DSR over mountains where it is highly heterogeneous and so the shortwave flux measured at ground stations does not match kilometer-scale DSR products. To tackle this challenge, we used a high spatial resolution (30 m) daily DSR over Sierra Nevada, Spain for 2008–2015, and a mountainous radiative transfer model to explore how topographic effects impacted the performances of DSR products. Four widely-used satellite products were selected as proxies for our evaluation: (i) MCD18A1 V6.1 (with a spatial resolution of 1 km); (ii) MSG DSR (~ 3.3 km); (iii) GLASS DSR V42 (0.05°); and (iv) BESS DSR (0.05°). There are three main findings under clear skies. Firstly, the product accuracies were slope-dependent, decreasing by 59.8–134.6% with slope ≥ 25° compared to areas with slope < 10°. Secondly, the product accuracies were aspect-dependent, exhibiting a higher degree of overestimation (i.e., average of 27.6 W/m²) on the north side and underestimation (i.e., average of -1.3 W/m²) on the south side. Thirdly, and finally, the product accuracies were time-dependent, exhibiting seasonal variations and pronounced overestimation in summer (i.e., 8.8 to 18.2 W/m²). Moreover, the impact of topography decreased with increasing cloud cover. Our findings can be applied to various mountainous areas due to the same mechanism of how topography influences the DSR estimation. This study corroborates the substantial uncertainties of the current DSR products in mountains and the necessity of incorporating topographic information into DSR estimations.

Semi-empirical bidirectional reflectance distribution function (BRDF) models are developed based on various spatial-resolution pixels. Because of its simplicity and physical significance, it is widely used in medium- and low-spatial-resolution quantitative remote sensing. With the emergence of high-spatial-resolution remote sensing data and the lack of high-spatial-resolution BRDF models, semi-empirical BRDF models have also been directly applied to high-spatial-resolution qualitative and quantitative remote sensing research. However, whether semi-empirical BRDF models can be directly applied to pixels with high resolution remains unclear. To answer this question, this letter quantitatively evaluates the applicability of semi-empirical BRDF models for remote sensing data with 0.5–30 m spatial resolution based on the WIDAS multi-angular observation dataset obtained during the HiWATER experiment in 2012. The results demonstrate that the semi-empirical BRDF models are not applicable at the 0.5 m pixel scale but are applicable at the 10 m pixel scale. There is a transitional pixel scale from not applicable to applicable between 0.5 and 10 m. We define this scale as the optimal minimum pixel scale (OMS) of semi-empirical BRDF models. The OMS is related to the spatial structure of the vegetation scene, and it is highly consistent with the canopy characteristic scale calculated based on the semivariogram method ($R^{2}=0.901$). Therefore, the range of the semivariogram can be used to estimate the OMS to answer the question of which scale semi-empirical BRDF models are applicable to high-spatial-resolution images.

The estimation of satellite-based albedo highly depends on the surface reflectance (SR). In mountainous areas, three types of SRs [i.e., the virtual SR (VSR) that is retrieved from the atmospheric correction model, the topographically corrected SR (TCSR) that is retrieved from the atmospheric and topographic correction model, and the sloping SR (SSR) that is retrieved from the physically bidirectional reflectance distribution function (BRDF)-based mountain-radiative-transfer (MRT) model] are commonly used to retrieve land surface albedo (SA). However, which type of SR is the best option for SA retrieval has not yet been quantitatively addressed. This letter assessed the performance of these three types of SRs on driving SA by comparison within situalbedo measurements over field sites in the Heihe River Basin, China. Our results show that these three types of albedos have consistent accuracy over flat sites with a root mean squared error (RMSE) smaller than 0.0320. Moreover, the sloping SA (SSA) retrieved from SSR shows the best agreement within situalbedo measurements over rugged sites with a bias of 0.0008, RMSE of 0.0338, relative RMSE (RMSER) of 12.92%, and correlation coefficient ($r$) of 0.89, followed by the topographically corrected SA (TCSA) from TCSR with a lager bias of 0.0208, RMSE of 0.0470, RMSERof 20.24%, and$r$of 0.69. The virtual SA (VSA) retrieved from VSR shows the largest uncertainty than the other two types of albedos, with an RMSE of 0.0516. These results illustrate that SSR is the best option of reflectance for satellite-based albedo retrieval over mountainous areas.

Bidirectional reflectance distribution function (BRDF) models are used to correct surface bidirectional effects and estimate land surface albedo. Many operational BRDF/albedo algorithms adopt a Roujean linear kernel-driven BRDF (RLKB) model because of its simple form and good performance in fitting multidirectional surface reflectance values. However, this model does not explicitly consider topographic effects, resulting in errors when applied over rugged terrain. To address this issue, we proposed a hybrid algorithm suitable for both flat and rugged terrain, called topographical kernel-driven (Topo-KD). First, we constructed a linear kernel-driven BRDF model considering terrain (LKB_T) which describes the topographic effects with a mountain radiative transfer (MRT) model. Then, the Topo-KD algorithm adaptively selects the most suitable model (RLKB or LKB_T) according to the terrain conditions and fitting residuals. The performances of Topo-KD and RLKB using the RossThick–LiSparseReciprocal (RTLSR) kernel are compared using simulated data sets and moderate-resolution imaging spectroradiometer (MODIS) observations. The results show that the BRDF of the pixel is affected by topography. But the RTLSR model does not specifically account for it, resulting in larger biases over rugged terrain than the Topo-KD algorithm in both the red and near-infrared (NIR) bands. The experiment using MODIS data sets demonstrates that the Topo-KD algorithm reduces fitting residuals in the red and NIR bands by 21.5% and 27.4% compared with the RTLSR model. These results indicate that the Topo-KD algorithm can be a better choice for retrieving land surface parameters and describing the radiative transfer process in mountainous areas.

The semiempirical kernel-driven model is commonly used for global surface reflectance characterization because of its simplicity and underlying physical meaning. However, the current kernel-driven reflectance models assume that the terrain is flat and homogeneous, and can induce significant errors in the surface reflectance estimation and subsequent parameter retrievals over rugged terrain. In this study, an improved topography-coupled kernel-driven (TCKD) reflectance model with the correction of diffuse skylight effects was proposed based on the diffused-equivalent slope model (dESM) and RossThick-LiTransit (RTLT) kernel-driven model. The TCKD model's accuracy and effectiveness were evaluated using surface reflectance simulated by the radiosity approach and the Moderate Resolution Imaging Spectroradiometer (MODIS) data. Against simulated data, the results show that the TCKD model can accurately capture the distortion of the reflectance shape and hemispherical distribution caused by the topographic effects. Compared to MODIS data, the TCKD model has an overall better performance than the RTLT model across different spatial scales and land cover types. When the mean slope is larger than 35° at the 500-m resolution, the TCKD model's near-infrared (NIR) root-mean-square error (RMSE) and the regression slope of the fitting line are 0.037 and 0.752, respectively, whereas those of the RTLT model are 0.049 and 0.645. Neglecting the diffuse skylight in the TCKD model can also lead to great bias in the reflectance retrievals. When the mean slope is 31°, as the ratio of diffuse skylight varies from 0 to 1, the NIR RMSE of the TCKD model decreases from 0.012 to 0.005, whereas that increases from 0.012 to around 0.02 if the diffuse skylight effects are neglected. These preliminary results demonstrate that the TCKD model is capable of improving the fitting ability of the kernel-driven model over rugged terrain and provides potentials for better retrieving and interpreting land surface parameters such as land surface albedo in mountainous areas.

The availability of global high-resolution land cover maps provides promising a priori knowledge for characterizing subpixel heterogeneity and improving predictions of directional reflectance of coarse-resolution pixels. Due to mutual shadowing and sheltering effects between the adjacent forest and cropland patches, the spectral nonlinear mixing of patchy ecotones is significant, especially when the sun illuminates the ecotone from the forest side with high solar zenith angle. The spectral linear mixture (SLM) approach leads to overestimation of the bidirectional reflectance factor (BRF) in the red band in the principal plane (PP), with a maximum absolute error (MAE) of 0.0063 and a maximum relative error (MRE) of 52.5%, and to underestimation in the near-infrared band in PP with an MAE of 0.0940 and an MRE of 14.5%. In a scenario with randomly distributed boundary orientations, the overestimation of SLM increases with the degree of fragmentation and the view zenith angle. We propose a Radiative Transfer model for patchy ECotones (RTEC). which improves R2from 0.61 to 0.94 in the red band of Landsat-8 directional reflectance at the validation site. The RTEC model provides an efficient and analytical approach for directional reflectance predictions over heterogeneous patchy landscapes at coarse resolution and will be used for biophysical parameter retrievals [e.g., the leaf area index (LAI)] in future applications.

Rugged terrain significantly complicates the land surface albedo modeling and retrievals in remote sensing. Neglecting the topographic effects may lead to large uncertainties when estimating land surface albedo over rugged terrain. In this study' the sensitivities of snow-free vegetation shortwave albedo to topography are quantitatively investigated and analyzed based on the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface albedo data and the global sensitivity analysis of a mechanistically-based land surface albedo parameterized model over rugged terrain. The results reveal that the topography can account for more than 30% of the total variabilities of the land surface albedo and the topographic effects on land surface albedo cannot be neglected even for the terrain with the mean slope of 10°.

Accurately estimating the spatial–temporal distribution of downward surface shortwave radiation (DSSR) is essential for terrestrial ecological modeling and climate change research. The accurate georegistration of digital elevation model (DEM) has become one of the significant bottlenecks for improving the DSSR accuracy over rugged terrain. To clearly understand and quantitatively evaluate the impact of geolocation bias on the DSSR estimation under clear sky, this letter conducts a systematical simulation research in Dayekou Basin of China based on a developed remote sensing satellite-based DSSR estimation scheme over rugged terrain. The results demonstrate that the proposed approach can accurately capture the high temporal and spatial heterogeneities of DSSR, and the DSSR estimations are sensitive to geolocation bias. When the horizontal bias is lower than half a pixel, the deviations of the direct radiation could lead to above 600 W/m2due to the illumination angle effects and shadow effects. The consequence of the bias on the diffuse and reflected radiation from adjacent terrains is little because of their relatively small values and low-spatial heterogeneities under clear sky in general except for the deep valley areas. The trends of the total radiation errors with the geolocation bias are identical in different days (scenes), and the error is related to the solar zenith angle. In addition, the more rugged the terrain, the greater the influence of geolocation bias on the radiation accuracy.

Topographic effects on canopy reflectance play a pivotal role in the retrieval of surface biophysical variables over rugged terrain. In this paper, we proposed a new canopy anisotropic reflectance model for discrete forests, Geometric Optical and Mutual Shadowing and Scattering-from-Arbitrarily-Inclined-Leaves model coupled with Topography (GOSAILT), which considers the effects of slope, aspect, geotropic nature of tree growth, multiple scattering, and diffuse skylight. GOSAILT-simulated areal proportions of four scene components (i.e., sunlit crown, shaded crown, sunlit background, and shaded background) were evaluated using the Geometric Optical model for Sloping Terrains (GOST) model. The canopy reflectances simulated by GOSAILT were validated against two reflectance data sets: Discrete anisotropic radiative transfer (DART) simulations and wide-angle infrared dual-model line/area array scanner (WIDAS) observations. Compared with a horizontal surface, the forest canopy reflectance over a steep slope (60°) is significantly distorted with absolute (relative) bias values of 0.048 (79.60%) and 0.056 (12.02%) for the red and near-infrared (NIR) bands, respectively. The GOSAILT-simulated component areal proportions show close agreements with GOST. Moreover, GOSAILT simulations have high overall accuracy (red band: coefficient of determination (R2) = 0.96; root-mean-square error (RMSE) = 0.003; and mean absolute percentage error (MAPE) = 3.91%; and NIR band: R2= 0.78, RMSE = 0.019; MAPE = 3.94%) when compared with the DART simulations. These extensive validations indicate good performances of GOSAILT in canopy reflectance simulations over sloped surfaces.

Rugged terrain complicates the BRDF modeling mainly by the modulation of sun-target-sensor geometry and shadowing effects. An improved kernel-driven BRDF model coupled with topography (KDCT) is put forward by combining the RTLSR model used in the algorithm for MODIS bidirectional reflectance anisotropies of land surface (AMBRALS) and the anisotropic reflectance model for rugged terrain (dESM). The improved model was compared with the original RTLSR model by using the simulated data based on the radiosity approach and the MODIS reflectance data. The validation results revealed that the improved KDCT model outperforms the RTLSR model without topographic consideration and can significantly improve the ability of the kernel-driven model to process the multi-angular reflectance measurements over rugged terrain.

Heterogeneous terrain significantly complicates signals received by airborne or satellite sensors. It has been demonstrated that both solar direct beam and diffuse skylight illumination conditions are significant factors influencing the anisotropy of reflectance over mountainous areas. Several models and methods have been developed to account for topographic effects on surface reflectance at the pixel level in remote sensing. However, subtopographic effects are generally neglected for low-spatial-resolution pixels due to the complex law of radiative transfer and the limitations of higher spatial resolution digital elevation models, which can lead to deviations in reflectance estimation. Accurately estimating the subtopographic effects on anisotropic reflectance over composite sloping terrain under different illumination conditions presents a challenge for remote sensing models and applications. In this paper, the diffused equivalent slope model (dESM) was developed, which is an anisotropic reflectance simulation model coupled with diffuse skylight over composite sloping terrain. The corresponding subtopographic impact factor was also proposed to exhibit how microslope topography affects reflectance over composite sloping terrain under different illumination conditions. Simulated reflectance data sets simulated by the radiosity method and Moderate Resolution Imaging Spectroradiometer reflectance data were used to evaluate the performance of the dESM model. The results reveal that the dESM model can accurately capture the reflectance anisotropy over composite sloping terrain under different illumination conditions, and the subtopographic impact factor can account for the effects of microslope topography, shadow, and illumination conditions.