Spaceborne LiDAR provides a promising method for large-scale characterizing LAI. However, the quality of point cloud data from spaceborne LiDAR, especially ICESat-2, is susceptible to atmosphere and background noise, introducing considerable uncertainty in LAI retrieval. Thus, efficiently screening out the high-quality point cloud is a significant guarantee for high-quality LAI retrieval. In this study, we proposed a quality control (QC) method that employed the number of 10 m windows without ground points in the ICESat-2 100 m segment as the QC flag. This method divided segments into 11 QC flags from 0 to 10 and was applied to LAI retrieval across Chinese forests from 2019 to 2020. The field measurements at locations identical to ICESat-2 ground tracks were used to validate the ICESat-2 LAI at different QC flags. The results showed that the proposed method effectively improved point cloud quality recognition and LAI accuracy, with ICESat-2 LAI (QC < 3) reducing RMSE by 26.36% compared to all ICESat-2 LAIs. It also showed good agreement with MODIS and GLASS LAI and mitigated saturation issues in passive optical imagery. The ICESat-2 LAI with QC < 3 performed better in deciduous broadleaved, evergreen needle-leaved, deciduous needle-leaved, and mixed forests, but not in evergreen broadleaved forests. ICESat-2 LAI was particularly adept at capturing high LAI values, which had the highest proportion of LAI values over 6.0 compared to MODIS and GLASS LAI. The proposed method has the potential for large-scale and high-quality LAI retrieval using ICESat-2 data on a global scale.

In this paper, we develop the Slant Leaf Area Index based Method (SLAIM) to estimate the leaf area of individual trees from single-scan Terrestrial Laser Scanner (TLS) data by introducing the concept of Slant Leaf Area Index (SLAI). SLAI quantifies the amount of leaves along the view direction and can be retrieved at given view zeniths using gap probability. Subsequently, leaf area can be accumulated by SLAI across the whole crown. The innovative SLAIM offers several advantages. Firstly, it operates with single-scan point clouds, which are more accessible than multiple-scan data. Secondly, it effectively corrects the clumping effect resulting from non-uniform leaf distribution. Both simulated and field-measured TLS point clouds of trees are used to test the method. The results show that the error of SLAIM is less than 10% in most cases.

Leaf area parameters are crucial in ecosystem studies. As ecophysiological models advance toward finer detail, accurately estimating LA at various scales becomes essential, particularly for diverse units like urban individual trees. Several algorithms based on terrestrial laser scanning (TLS) data have been developed to obtain the LA of individual trees. However, their use at the stand level needs further research. In this study, the comparative shortest-path algorithm (CSP) is introduced for the automatic individual tree segmentation, thereby facilitating the application of the path length distribution model (PATH) for leaf area estimation at the stand level. Using high-density TLS data, we presented a bottom-up estimation of stand leaf area index (LAI) from 50 individual tree measurements and validated the results at different scales. At the tree scale, the LA derived from TLS and allometric model were highly correlated, with an R-value of 0.83. At the stand scale, the proposed method provides consistent results with the allometric and TRAC instrument measurements, performing better than vertical upward photography. Generally, 23 shared stations under the forest are enough to accurately obtain the LA of 50 trees and the LAI in an urban forest stand. Sensitivity analysis shows that the method is not sensitive to TLS scan resolution and parameters used in tree crown envelope reconstruction. The proposed bottom-up approach provides a new way of estimating the LAI at stand level using TLS and has the advantage of providing multi-level leaf area information and avoiding the scale effect.

Individual trees are fundamental to urban ecosystems as they play an important role in energy transfer, pollutant removal, and habitat formation. Leaf area (LA) is an important factor to quantify the effect of individual trees on urban ecosystems. Terrestrial laser scanners (TLSs) are widely recognized as the most accurate devices for tree structural measurements. However, they face challenges in estimating LA from LA index (LAI) for individual trees primarily due to arbitrary and confusing horizontal projection areas. Occlusion and clumping effects further hinder the objective and accurate LA measurements of individual trees. Therefore, we developed the slant leaf area index-based method (SLAIM) to estimate the LA of individual trees from single-scan TLS data by introducing the concept of slant leaf area index (SLAI). SLAI quantifies the amount of leaves along the view direction, and it can be retrieved at given view zeniths using gap probability. Subsequently, LA can be accumulated by SLAI across the whole crown. Tests with simulated and field-measured TLS point clouds demonstrate SLAIM’s accuracy, with the relative errors (REs) in LA below 10% in most cases. Stratified LA validation reveals an$R^{2}$exceeding 0.77 across all realistic crowns, along with a root-mean-square error (RMSE) under 2 m2. SLAIM’s advantages include compatibility with single-scan point clouds, effective correction of clumping effects, and consideration of variations in leaf projection coefficients at different zeniths. SLAIM proves more efficient and practical for actual LA measurements, showcasing its potential for advanced urban ecosystem research.

The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) employs a unique multibeam photon counting approach to acquire a near-continuously sampled profile and provides more precise technology for mapping the leaf area index (LAI) at the global scale. The inversion accuracy of LAI is affected by the clumping effect, which has been an open question for spaceborne laser scanning (SLS). Here, we present a segmented method based on the path length distribution model to calculate the clumping-corrected LAI independently using ICESat-2 data. The results showed that the LAI derived by the proposed method with a 200 m segment was consistent with the airborne laser scanning (ALS)-derived LAI, with a root mean squared error (RMSE) of 0.37. A satisfactory agreement (RMSE$=1.03$) was also shown between moderate resolution imaging spectroradiometer (MODIS) LAI and ICESat-2 LAI. Moreover, the LAI derived by the proposed method was on average 31.72% higher than the LAIe derived by Beer’s law, which indicated that the proposed method achieved the purpose of correcting the clumping effect. The gap probability was calculated by the 200 m moving window and the path length distribution was obtained by the 1 m moving window as the model input had the highest accuracy. In addition, the limitation of the point cloud data and the time lag of ICESat-2 acquisitions and ALS observations may affect the inversion accuracy of LAI. This study proposed a feasible way to correct the clumping effect and invert LAI independently using ICESat-2 data, which has the potential to characterize vegetation structure precisely at regional and global scales.

The Qinghai–Tibet Plateau is rich in water resources with numerous lakes, rivers, and glaciers, and, as a source of many rivers in Central Asia, it is known as the Asian Water Tower. Under global climate change, it is critical to understand the current influencing factors on surface water area in this region. Although there are numerous studies on surface water mapping, they are still limited by temporal/spatial resolution and record length. Moreover, the complicated topographic condition makes it challenging to map the surface water accurately. Here, we proposed an automatic two-step annual surface water classification framework using long time-series Landsat images and topographic information based on the Google Earth Engine (GEE) platform. The results showed that the producer accuracy (PA) and user accuracy (UA) of the surface water map in the Qinghai–Tibet Plateau in 2020 were 99% and 90%, respectively, and the Kappa coefficient reached 0.87. Our dataset showed high consistency with high-resolution images, indicating that the proposed large-scale water mapping method has great application potential. Furthermore, a new annual surface water area dataset on the Qinghai–Tibet Plateau from 2000 to 2020 was generated, and its relationship with climate, vegetation, permafrost, and glacier factors was explored. We found that the mean surface water area was about 59 481 km2, and there was a significant increasing trend (=322 km2/year,$p < 0.01$) during 2000–2020 in the plateau. Greening, warming, and wetting climate conditions contributed to the increase of surface water area. Active layer thickness and permafrost types may be the most related to the decrease of surface water area. This study provides important information for ecological assessment and protection of the plateau and promotes the implementation of sustainable development goals related to surface water resources.

Soil salinity plays a key role in influencing the soil dielectric constant and soil backscatter coefficient. However, soil moisture (SM) retrieval models constructed based on active microwave data hardly consider soil salinity. Thus, obtaining the SM datasets with various salinity on regional and local scales is difficult. This study aimed to employ theoretical model simulation to investigate the errors of SM retrieval due to not considering the impact of soil salinity. Then, three typical saline soil dielectric constant models were validated and compared based on the experimental measurement datasets. Results show that the WYR saline soil dielectric constant model has excellent performance. The soil salinity mainly affects the imaginary part of the dielectric constant and the effect of salinity on the soil dielectric constant is more significant when the SM has larger values. In addition, in retrieving SM with soil salinity more than 10 g/kg, the retrieval result of SM has an absolute error of 0.04$\text{m}^{3}/\text{m}^{3}$and a relative error of 5% when not considering the soil salinity impact. In retrieving SM with soil salinity less than 10 g/kg, the retrieved SM error increased by 2%, and the absolute error increased by 0.01$\text{m}^{3}/\text{m}^{3}$as soil salinity increased by 3 g/kg. We believe that The study will give a theoretical reference for establishing the SM retrieval model in saline soil areas using microwave data.

Canopy clumping index (CI) characterizes the extent of the nonrandom spatial distribution of foliage elements within a canopy and is critical for determining the radiative transfer, photosynthesis, and transpiration processes in the canopy. It is widely perceived that CI increases with zenith angle (θ), because between-crown gaps decrease in size and number with increasing θ. In this study, we demonstrate that this is not always true. Analytical equations between CI and θ are first developed based on widely-used forest canopy gap fraction theories. The results show that the zenith angular variation of CI is closely related to crown projected area or crown shapes (i.e., the ratio of the crown height to its diameter, RHD): CI increases with θ for canopies with “tower” crowns (RHD > 1), but decreases with θ for “umbrella” crowns (RHDin-situmeasurements and multi-angular remote sensing.

Clumping effect denotes the nonrandomness of foliage. It deviates from the random distribution assumption of Beer’s law which is usually applied to leaf area index (LAI) retrieval from large-footprint full-waveform light detection and ranging (LiDAR). Some studies correct for large gaps-induced between-crown clumping, yet ignore the within-crown clumping. The error of LAI caused by these clumping effects and the influence of the forest structure parameters on them have not been quantitatively studied. This study quantified the between-crown, within-crown, and total clumping indices through a theoretical derivation, clarifying the mechanism of clumping; we used airborne LiDAR point clouds data in 11 290 footprints (diameter = 25 m) to estimate these indices in real forests. We found that: 1) the underestimation of LAI caused by directly applying Beer’s law could be up to 93%, and it decreases with fractional crown coverage but increases with crown length and leaf area density; 2) the method of correcting between-crown clumping improves LAI retrieval for cylindrical canopies effectively; however, 3) considerable underestimation (up to 58%) exists if we neglect the within-crown clumping for other canopies, which has not been realized before; and 4) both the between-crown and the within-crown clumping can be the dominant contributor, and the within-crown clumping was greater than the between-crown clumping in 47% of the studied footprints. In the two physically based LAI retrieval methods, Beer’s law has been commonly used due to its simplicity. Pathways to improve future LAI retrieval would be instrument improvement to capture the between-crown gaps and method study to correct the within-crown clumping further.

To improve our capacity to map long-term vegetation dynamics in heterogeneous landscapes, this study proposed a new prior knowledge-based spatiotemporal enhancement method, namely, PK-STEM, to fuse MODIS and Landsat FPAR products following the remote sensing trend surface framework. PK-STEM uses historical Landsat FPAR images as prior knowledge and fuses them with new satellite-derived FPAR data. PK-STEM can work in three modes: 1) using only MODIS data; 2) using only Landsat data; and 3) using both MODIS and Landsat data. This study retrieved FPAR from Landsat images using a scaling-based method and tested the performance of PK-STEM in a regional application. For the entire year of 2012, we compared the performance of PK-STEM in different modes and with that of two typical spatiotemporal fusion methods, the enhanced spatial and temporal adaptive reflectance model (ESTARFM) and unmixing-based linear mixing growth model (LMGM). Then, a long time series FPAR data set at 30-m resolution and eight-day intervals was generated for 13 years (2000–2012). Our results show that PK-STEM in mode III is the most robust and accurate (root mean squared error (RMSE) = 0.062; mean$R = 0.851$) among the three modes and more accurate than ESTARFM (mean RMSE = 0.065; mean$R = 0.776$) and LMGM (mean RMSE = 0.074; mean$R = 0.734$). For the 12 years (2000–2011), PK-STEM also achieves high accuracies with mean RMSE = 0.066 and$R = 0.938$. PK-STEM is very flexible with a continual update mechanism and is efficient for long time series applications.

Leaf inclination angle distribution plays an important role in indirect leaf area measurement methods. Path length distribution method (PATH) is an indirect leaf area measurement method based on Beer's Law, which has been applied to isolated trees with Terrestrial Laser Scanning (TLS). However, it can only set the leaf projection G to a constant value. In this paper, the PATH method considering leaf inclination models is introduced, which allows the G function to be a variable corresponding to leaf inclination. On the basis of this method, this paper explores the effect of different leaf inclination angle distribution assumptions on the inversion of the leaf area of isolated trees based on TLS. The results show that compared with the measured leaf inclination, the relative errors of the inversion results based on the six typical leaf inclination assumptions are between −14.3 % to +41.2%, which indicates that leaf inclination has a significant effect on the inversion of the leaf area. Further, experiments in this paper show that the mean value of the G function is a relatively accurate representation of it.

Occlusion effect, an inherent problem of terrestrial laser scanning (TLS) measurements, limits the potential of TLS data in tree attribute estimation. Multiple scans seek to mitigate this effect to provide enhanced scan completeness. However, the numbers and locations of the scans (i.e., the scan design) are usually determined via a subjective assessment of the tree density, spatial patterns of trees, and attributes to be derived. These could cause suboptimal scan completeness and limit tree attribute estimation. This study proposed an iterative-mode scan design to minimize the occlusion effect. First, we introduced a PoTo index based on visibility analysis to evaluate how many trees can be scanned from a location and to select effective candidates for the optimal TLS location. Second, we introduced a cumulative degree of ring closure (CDRC) to quantify the scan completeness for each candidate and determine the optimal TLS location. The TLS data sets of virtual forests with field-measured and synthetic plot parameter settings were simulated according to iterative- and regular-mode designs by using a Heidelberg light detection and ranging (LiDAR) Operations Simulator (HELIOS). The results demonstrated that an iterative-mode design can improve the scan completeness of trees compared to the regular-mode design. The tree attribute (diameter at breast height (DBH), tree height, stem curve, and crown volume) estimates of the iterative-mode design were less erroneous than those of the regular-mode design (e.g., the root-mean-square error (RMSE) could decrease the stem curve estimation by 38% and the crown volume estimation by 15%). This study suggests that the iterative-mode design can obtain an improved quality of the TLS data, especially for dense stands.

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 km2) 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 airborne laser scanner (ALS) provides great potential for mapping the leaf area index (LAI) at the landscape scale using grid cell statistics, while its application is restricted by the lack of clumping information, which has been an unsolved issue highlighted for a long time. ALS generally provides an effective LAI because its footprint is too large to capture small gaps to apply traditional ground-based clumping correction methods. Here, we present a grid cell method based on path length distribution model to calculate the clumping-corrected LAI using ALS data without the requirement of additional field measurements. We separated the within- and between-crown areas to consider between-crown clumping, and used the path length distribution as estimated by local canopy height distribution to consider 3-D foliage profile and within-crown clumping. The path length distribution model takes advantage of the 3-D information rather than the gap size distribution, thus avoiding the limitation of large ALS footprint. With the 0.4-m-footprint ALS data, the results are generally promising and a multilevel clumping analysis is consistent with landscape flown. The ALS LAIs of different resolutions are consistent, with a difference of less than 5% from 5- to 250-m resolutions. Due to its consistency and simple configuration, the method provides an opportunity to map the clumping-corrected LAI operationally and strengthens the ability of airborne lidar to monitor vegetation change and validate the satellite product. This grid cell method based on path length distribution is worth further testing and application using more recent laser technology.

Spatial scale effects, defined as the phenomenon that the estimates at multiple resolution are inconsistent, have aroused wide concerns in remote sensing studies. But very few studies have paid attention to the effects of scale on phenological studies. This paper investigated the scale effects in estimating phenological transitional dates from remote sensing data. A prior-knowledge vegetation index (VI) time series at 30 m resolution was composed, based on which the time series at 240 m, 480 m and 960 m were derived. The green-up onset and dormancy onset dates were then estimated from the VI time series using a double-logistic plant growth model. The derived estimates at multiple resolutions were compared and the effects of spatial scales were verified. Landscape heterogeneities were found to be related to spatial scale effects. The changes in the estimated green-up onset dates and dormancy onset dates exhibited different patterns with the coarsening of spatial resolution.

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.

Non-randomness within canopies and woody component are two factors limiting the accuracy of indirect leaf area index (LAI) measurement. Here we combine the path length distribution model and Multispectral Canopy Imager (MCI) together for the first time to improve the accuracy. The results show that non-randomness within canopies underestimates 17.1%-28.2% LAI, while woody component overestimates 14.6%-27.8% LAI in four forest sites. Although these two factors were sometimes offset, the degree of non-randomness within canopies and the proportion of woody component vary in different forests. More attention should be paid to the impact of the non-randomness within canopies and the woody component, especially in coniferous forest dominated by tree trunks and branches.

The ground-measured emissivity is always affected by many kinds of noises, which lead the retrieval accuracy to be out of expectation. This paper investigates the influence of three major noises (formula simplification, surface temperature measurement, and temperature emissivity separation algorithm) on the spectral emissivity by using simulation data based on radiative transfer model and field measured data from portable 102F infrared spectrometer. The findings of this paper can provide some suggestions for the further emissivity measurement.

Accurate spectral calibration of the in-flight sensors is crucial for processing and exploration of remotely sensed data. This paper developed a strategy to make spectral recalibration (i.e., spectral response function, central wavelength, and bandwidth) for in-flight broadband sensor using a device-responsivity-decomposition model with a priori knowledge and an optimization algorithm. Sensitivity analysis indicates that an accurate result requires the targets to be observed under a dry and clear atmospheric condition (column water vapor2and visibility > 23 km) and no more than 5% error is included in the measured data. The new strategy was used to retrieve the spectral parameters along with radiometric calibration coefficients for a multichannel camera onboard an unmanned aerial vehicle from simultaneously remotely sensed and ground measured data sets over 19 (15 color-scaled and four gray-scaled) man-made surface targets, and the retrieved results were validated with a similar data set over another four man-made targets. It demonstrated that the camera's spectral parameters were accurately retrieved and an error less than 3.5 W/m2/μm/sr was brought to the channel radiance.

This paper presents a portable Multi-Angle Observation System (MAOS) to quickly collect bi-directional reflectance factor (BRF) and directional thermal radiance of land surface along with the spectroradiometer and thermal radiometer. The new system is able to make more than 13 zenith measurements in six minutes at an arbitrary azimuth direction, with the angle-controlling accuracy better than 2°. More observations are sampled in the hot-spot direction. All operations of the MAOS and data-processing are automatically controlled by the computer. Field campaign of winter wheat canopy shows that the MAOS had captured the angular variations of the BRF.

Gap fraction is a very important parameter to the indirect estimation of the true Leaf Area Index. In this paper, we combined the multispectral digital imageries (RGB color imagery and Near-Infrared imagery), which were obtained from a new device called Multispectral Canopy Imager (MCI), to estimate gap fraction. A new method incorporated with CIE L*a*b* color space has also been proposed to segment the multispectral digital imagery. The preliminary results of the estimated gap fraction have been showed in the conclusions section and been proved to be very well.