Estimating global aerodynamic parameters in 1982–2017 using remote-sensing data and a turbulent transfer model

The Global Land Surface Satellite (GLASS) leaf area index (LAI) product is one of the most widely used global LAI products in the scientific community. The latest (version 5) GLASS LAI product has been generated from Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data. The purpose of this paper is to evaluate the quality of the version 5 GLASS LAI product. The GLASS LAI product was compared with the latest MODIS LAI product (MCD15A2 H, Collection 6) and the second version of Geoland2 (GEOV2) LAI product to evaluate their temporal and spatial discrepancies. A direct validation was conducted to compare these LAI products to the LAI values derived from the high-resolution reference maps from the Validation of Land European Remote Sensing Instruments (VALERI) and Implementing Multi-Scale Agricultural Indicators Exploiting Sentinels (IMAGINES) sites. The results show that the GLASS and GEOV2 LAI products have great spatial integrity. However, the MODIS LAI product contains many missing pixels in tropical areas. These LAI products follow fairly consistent seasonal characteristics. The spatial discrepancies of these LAI products mainly exist in forest areas, especially evergreen broadleaf forests where the GLASS LAI values are generally lower than the GEOV2 LAI values by approximately 1.0 LAI units and lower than the MODIS LAI values by 0.5 to 1.0 LAI units. The spatial distribution of these LAI products has slight discrepancies in savannahs, broadleaf crops, grasses/cereal crops and shrubs. The GLASS and GEOV2 LAI products capture a complete and reasonable temporal profile, contrasting with the MODIS LAI product, which shows dramatic fluctuations, particularly during the growing seasons. These LAI products show similar temporal trajectories and interannual variations for all biome types except evergreen broadleaf forests. The direct validation shows that the accuracy of the GLASS LAI product is better than the accuracy of the MODIS and GEOV2 LAI products. The coefficient of determination (R2) of the GLASS, MODIS and GEOV2 LAI products versus the LAI values derived from the high-resolution reference maps are 0.68, 0.47 and 0.55, respectively, and the root mean square error (RMSE) of these products are 0.86, 1.22 and 1.21, respectively.

The leaf area index (LAI) is a critical variable for forest ecosystem processes. Passive optical and active LiDAR remote sensing have been used to retrieve LAI. LiDAR data have good penetration to provide vertical structure distribution and deliver the ability to estimate forest LAI, such as the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2). Segment size and beam type are important for ICESat-2 LAI estimation, as they affect the amount of signal photons returned. However, the current ICESat-2 LAI estimation only covered a limited number of sites, and the performance of LAI estimation with different segment sizes has not been clearly compared. Moreover, ICESat-2 LAIs derived from strong and weak beams lack a comparative analysis. This study derived and evaluated LAI from ICESat-2 data over the National Ecological Observatory Network (NEON) sites in North America. The LAI estimated from ICESat-2 for different segment sizes (20, 100, and 200 m) and beam types (strong beam and weak beam) were compared with those from the airborne laser scanning (ALS) and the Copernicus Global Land Service (CGLS). The results show that the LAI derived from strong beams performs better than that of weak beams because more photon signals are received. The LAI estimated from the strong beam at the 200 m segment size shows the highest consistency with those from the ALS data (R = 0.67). Weak beams also present the potential to estimate LAI and have moderate agreement with ALS (R = 0.52). The ICESat-2 LAI shows moderate consistency with ALS for most forest types, except for the evergreen forest. The ICESat-2 LAI shows satisfactory agreement with the CGLS 300 m LAI product (R = 0.67, RMSE = 1.94) and presents a higher upper boundary. Overall, the ICESat-2 can characterize canopy structural parameters and provides the ability to estimate LAI, which may promote the LAI product generated from the photon-counting LiDAR.

Development of the GLASS 250-m leaf area index product (version 6) from MODIS data using the bidirectional LSTM deep learning model

Large-scale retrieval of leaf area index and vertical foliage profile from the spaceborne waveform lidar (GLAS/ICESat)

Long-term high-quality global leaf area index (LAI) and fraction of absorbed photosynthetically active radiation (FAPAR) products are urgently needed for the study of global change, climate modeling, and many other problems. As the successor of the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor, the Visible Infrared Imaging Radiometer Suite (VIIRS) will continue to provide global environmental measurements. This paper aims to generate longer time series Global LAnd Surface Satellite (GLASS) LAI and FAPAR products after the era of the MODIS sensor. To ensure spatial and temporal consistencies between GLASS LAI/FAPAR values retrieved from different satellite observations, the GLASS LAI/FAPAR retrieval algorithms were adapted in this study to retrieve LAI and FAPAR values from VIIRS surface reflectance time-series data. After reprocessing of the VIIRS surface reflectance to remove remaining effects of cloud contamination and other factors, a database generated from the GLASS LAI product and the reprocessed VIIRS surface reflectance for all Benchmark Land Multisite Analysis and Intercomparison of Products (BELMANIP) sites was used to train general regression neural networks (GRNNs). The reprocessed VIIRS surface reflectance data from an entire year were entered into the trained GRNNs to estimate the one-year LAI values, which were then used to calculate FAPAR values. A cross-comparison indicates that the LAI and FAPAR values retrieved from VIIRS surface reflectance were generally consistent with the GLASS, MODIS and Geoland2/BioPar version 1 (GEOV1) LAI/FAPAR values in their spatial patterns. The LAI/FAPAR values retrieved from VIIRS surface reflectance achieved good agreement with the GLASS LAI/FAPAR values (R2 = 0.8972 and RMSE = 0.3054; and R2 = 0.9067 and RMSE = 0.0529, respectively). However, validation of the LAI and FAPAR values derived from VIIRS reflectance data is now limited by the scarcity of LAI/FAPAR ground measurements.

Leaf area index (LAI) is a key input for many land surface models, ecological models, and yield prediction models. In order to make the model simulation and/or prediction more reliable and applicable, it is crucial to know the characteristics and uncertainties of remotely sensed LAI products before they are input into models. In this study, we conducted a comparison of four global remotely sensed LAI products—Global Land Surface Satellite (GLASS), Global LAI Product of Beijing Normal University (GLOBALBNU), Global LAI Map of Chinese Academy of Sciences (GLOBMAP), and Moderate-resolution Imaging Spectrometer (MODIS) LAI, while the former three products are newly developed by three Chinese research groups on the basis of the MODIS land reflectance product over China between 2001 and 2011. Direct validation by comparing the four products to ground LAI observations both globally and over China demonstrates that GLASS LAI shows the best performance, with R2 = 0.70 and RMSE = 0.96 globally and R2 = 0.94 and RMSE = 0.61 over China; MODIS performs worst (R2 = 0.55, RMSE = 1.23 globally and R2 = 0.03, RMSE = 2.12 over China), and GLOBALBNU and GLOBMAP performs moderately. Comparison of the four products shows that they are generally consistent with each other, giving the smallest spatial correlation coefficient of 0.7 and the relative standard deviation around the order of 0.3. Compared with MODIS LAI, GLOBALBNU LAI is the most similar, followed by GLASS LAI and GLOBMAP. Large differences mainly occur in southern regions of China. LAI difference analysis indicates that evergreen needleleaf forest (ENF), woody savannas (SAV) biome types and temperate dry hot summer, temperate warm summer dry winter and temperate hot summer no dry season climate types correspond to high standard deviation, while ENF and grassland (GRA) biome types and temperate warm summer dry winter and cold dry winter warm summer climate types are responsible for the large relative standard deviation of the four products. Our results indicate that although the three newly developed products have improved the accuracy of LAI estimates, much work remains to improve the LAI products especially in ENF, SAV, and GRA regions and temperate climate zones. Findings from our study can provide guidance to communities regarding the performance of different LAI products over mainland China.

The extension of NASA’s Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission beyond the cryosphere to include the study of vegetation and the land surface has significantly advanced global Earth observation. This review synthesizes findings from 293 peer-reviewed articles and provides a comprehensive assessment of the mission’s contributions to terrestrial ecosystem science. We begin by outlining the mission’s objectives, instrumentation, and core data products—particularly ATL03 (geolocated photons) and ATL08 (terrain and canopy heights)—which have supported a wide array of applications in forestry, terrain analysis, and environmental monitoring. Our analysis shows that most studies focus on temperate and tropical broadleaf forests, with a strong emphasis on estimating canopy height, terrain elevation, and forest structure. Terrain metrics derived from ICESat-2 products typically achieve high accuracy, while vegetation variables, such as canopy height, cover, leaf area index (LAI), and aboveground biomass, demonstrate moderate to strong accuracies across biomes. While parametric models remain the most used approach, machine learning methods are expanding, and more than one third of the reviewed literature incorporates synergistic analyses with other satellite missions. Despite its versatility, ICESat-2 remains underutilized in several key domains—including nonforest vegetation, fire ecology, wildlife habitat assessment, urban monitoring, and disturbance detection—mainly due to the scarcity of standardized algorithms and validated reference datasets. However, the growing availability of open source processing tools presents a significant opportunity for expanding the user base and fostering innovation in these emerging areas. Future advancements, such as the forthcoming ATL18 gridded canopy height and terrain products, alongside synergies with missions like the European Space Agency (ESA)’s BIOMASS and the NASA–Indian Space Research Organisation (ISRO) Synthetic Aperture Radar (NISAR), promise to enhance multisensor integration and ecosystem monitoring. Overall, ICESat-2 continues to evolve as a powerful resource for characterizing vegetation structure, ecosystem dynamics, and land surface processes on a global scale.

The Leaf Area Index (LAI) strongly influences crop biomass production and yields. The variation characteristic of LAI and the development of crop growth models can provide a theoretical basis for predicting crops’ water consumption, fruit quality and yields. This paper analyzes the relationship between measurements of aboveground grape biomass and trends in LAI and dry biomass production in grapes grown in the Turpan area. The LAI changes in grapes were estimated using the modified logistic model, the modified Gaussian model, the log-normal model, the cubic polynomial model, and the Gaussian model. Universal models of LAI were established in which the applied irrigation quota was applied to calculate the maximum LAI. The relationship between the irrigation quota and biomass production, yields, and the harvest index was investigated. The developed models could accurately predict the LAI of grapevines grown in an extremely arid area. However, the Gaussian and cubic polynomial models produced less accurate results than the other models tested. The Michaelis–Menten model analyzed the relationship between biomass and LAI, providing a numerical method for predicting dynamic changes in grapevine LAI. Moreover, the crop biomass increased linearly with the irrigation quota for quotas between 6375 and 13,200 m3/hm. This made it possible to describe the grape yield and harvest index with a quadratic polynomial function, which increases during the early stages of the growing season and then decreases. The analyses of the relationship between yield and harvest index provide important theoretical insights that can be used to improve water use efficiency in grape cultivation and to identify optimal irrigation quotas.

The leaf area index (LAI) is a critical parameter for characterizing the structure and function of vegetation in ecosystems. Currently, operational LAI products always have coarse spatial resolution, and fine-resolution LAI maps are urgently needed for ecological environment assessment and the precise monitoring of cropland growth. LAI downscaling methods are efficient at improving the spatial resolution of LAI products but often ignore the scaling effect of the model. In this study, a novel model-downscaling method is proposed for fine-resolution LAI estimation. It uses scaling equations of model parameters (SEMPs) to describe the scaling relations of models at different spatial resolutions and construct a downscaled model from a coarse-resolution model. Landsat Normalized Difference Vegetation Index (NDVI) at 30 m and Global LAnd Surface Satellite (GLASS) LAI at 1 km spatial resolutions are used because they are readily available. The downscaled model is evaluated by a fine-resolution model directly constructed with fine-resolution data. The fine-resolution LAI values estimated by this model-downscaling method are evaluated with field LAI measurements. The validation results show that the proposed method can generate highly accurate LAIs, with an RMSE of 0.821 at the Pshenichne cropland site in Ukraine and an RMSE of 0.515 at the Camerons forest site in Australia when compared with field LAI measurements. The results are also better than those of Ovakoglou’s downscaling method. These results demonstrate that the model-downscaling method for fine-resolution LAI estimation is viable and referable for related studies.

Variations in leaf area index (LAI) are critical to research on forest ecosystem structure and function, especially carbon and water cycle, and their responses to climate change. Using the ensemble empirical mode decomposition (EEMD) method and global inventory modeling and mapping studies (GIMMS) LAI3g dataset from 1982 to 2010, we analyzed the nonlinear feature and spatial difference of forest LAI variability over China for the past 29years in this paper. Results indicated that the national-averaged forest LAI was characterized by quasi-3- and quasi-7-year oscillations, which generally exhibited a rising trend with an increasing rate. When compared with 1982, forest LAI change by 2010 was more evident than that by 1990 and 2000. The largest increment of forest LAI occurred in Central and South China, while along the southeastern coastal areas LAI increased at the fastest pace. During the study period, forest LAI experienced from decrease to increase or vice versa across much of China and varied monotonically for only a few areas. Focusing on regional-averaged trend processes, almost all eco-geographical regions showed continuously increasing trends in forest LAI with different magnitudes and speeds, other than tropical humid region and temperate humid/subhumid region, where LAI decreased initially and increased afterwards.

The Multi-Angle Imaging Spectroradiometer (MISR) is a powerful sensor for leaf area index (LAI) mapping with its simultaneous multi-angle observations. However, the LAI product derived from MISR observations has low temporal resolution, which is unsatisfactory for many applications. This paper presents an algorithm that expands the MISR LAI product to high temporal resolution with the aid of Moderate Resolution Imaging Spectroradiometer (MODIS) data. The algorithm establishes relationships between the MISR LAI and the MODIS red/near-infrared band ratio (simple ratio (SR)) pixel by pixel using coincident data of these two sensors for the past nine years. Using these pixel-based SR-LAI relationships, a new LAI product with the merits of the original MISR product and high temporal resolution is obtained from MODIS surface reflectance. The expanded LAI series was compared with the original MISR and MODIS LAI products, as well as field LAI measurements made at the Baohe and Maoershan forest sites and the Hulunbeier grassland site, to assess the algorithm's performance. The results show that the temporal coverage of the MISR LAI improved from 15.5% to 65.2% in an 8-day composite, and the mean root-mean-square error is 0.74 for the vegetated pixels. This LAI product has similar temporal consistency and seasonal dynamics to the existing MODIS LAI product generated from the main algorithm, but is more robust against the low quality of reflectance inputs. The expanded LAI product differs with field measurements by about 11.5%, with agreement to field observations at all three sites within an accuracy of 0.8 LAI.

Abstract. Systematic evaluation of the carbon cycle physical and biological variables simulated in Earth System Model (ESM) participating in the Coupled Model Intercomparison Project Phase 6 (CMIP 6) is fundamental to the understanding of terrestrial ecosystems, as well as to future projections. Leaf Area Index (LAI), Gross Primary Productivity (GPP), Net Primary Productivity (NPP), Net Ecosystem Productivity (NEP) and Land Surface Temperature (LST) as key indicators of carbon cycle performance in ESM outputs, play a critical role in evaluating ecosystem functions. Assessing these metrics can provide valuable insights into the biases in model-simulated ecosystems and offer guidance for model improvement. In this study, we assessed the interannual trends performance of LAI, GPP, NPP, NEP and LST simulated by 12 CMIP6 ESMs during the historical period by using satellite LAI, NPP, NEP, LST and CSIF data as observations. The findings indicate that: (1) There are significant uncertainties in the overall trends and interannual variability in LAI, NPP, and LST captured by the CMIP6 ESM. Meanwhile, simulated GPP and NEP trends were lower than observations with discrepancies reaching 0.03 yr−1 for GPP and 2.46 gCm-2yr-1 for NEP. (2) Spatially, CMIP6 ESMs exhibited widespread underestimation of trends in LAI, GPP, NPP, and NEP across China. The MME underestimated these variables in 46.29 % (LAI), 43.47 % (GPP), 49.81 % (NPP), and 61.34 % (NEP) of the study area. Meanwhile, the simulated LST trend is underestimated in northern China, while its overestimations in western and southern China. (3) ESMs inadequate responsiveness to anthropogenic and environmental forcing and incomplete mechanistic representation of plant respiration pathways struggled accurate simulation of trends in LAI, GPP, NPP, NEP and LST.

Continuous large-scale mapping of forest canopy height is crucial for estimating and reporting forest carbon content, analyzing forest degradation and restoration, or to model ecosystem variables such as aboveground biomass. Over the last years, the spaceborne Light Detection and Ranging (LiDAR) sensor specifically designed to acquire forest structure information, Global Ecosystem Dynamics Investigation (GEDI), has been used to extract forest canopy height information over large areas. Yet, GEDI has no spatial coverage for most forested areas in Canada and other high latitude regions. On the other hand, the spaceborne LiDAR called Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) provides a global coverage but was not specially developed to study forested ecosystems. Nonetheless, both spaceborne LiDAR sensors obtain point-based information, making spatially continuous forest canopy height estimation very challenging. This study compared the performance of both spaceborne LiDAR, GEDI and ICESat-2, combined with ALOS-2/PALSAR-2 and Sentinel-1 and -2 data to produce continuous canopy height maps in Canada for the year 2020. A set-aside dataset and airborne LiDAR (ALS) from a national LiDAR campaign were used for accuracy assessment. Both maps overestimated canopy height in relation to ALS data, but GEDI had a better performance than ICESat-2 with a mean difference (MD) of 0.9 m and 2.9 m, and a root mean square error (RMSE) of 4.2 m and 5.2 m, respectively. However, as both GEDI and ALS have no coverage in most of the hemi-boreal forests, ICESat-2 captures the tall canopy heights expected for these forests better than GEDI. PALSAR-2 HV polarization was the most important covariate to predict canopy height, showing the great potential of L-band in comparison to C-band from Sentinel-1 or optical data from Sentinel-2. The approach proposed here can be used operationally to produce annual canopy height maps for areas that lack GEDI and ICESat-2 coverage.

Leaf area index (LAI) products retrieved from observations acquired on one occasion have obvious discontinuity in the time series owing to cloud coverage and other factors, and the accuracy may not meet the needs of many applications. Effectively utilizing data assimilation techniques to retrieve biophysical parameters from the time series of remote-sensing data has attracted special interest. The data assimilation technique is based on a reasonable consideration of dynamic change rules of biophysical parameters and time series observational quantities, thereby improving the quality of retrieved profiles. In this article, a variational assimilation procedure for retrieving LAI from the time series of remote-sensing data is developed. The procedure is based on the formulation of an objective function. A dynamic model is constructed based on the climatology from multi-year Moderate Resolution Imaging Spectroradiometer (MODIS) LAI data to evolve LAI in time, and a radiative transfer model is coupled with the dynamic model to simulate a time series of surface reflectances. A shuffled complex evolution method (developed at the University of Arizona; SCE-UA) optimization algorithm is then used to minimize the objective function and estimate the dynamic model states and the parameters of the coupled model from the MODIS reflectance data with a higher quality in a given time window. The variational assimilation method is applied to the MODIS surface reflectance data for the whole of 2008 at the Heihe river basin to produce regional LAI mapping results. The ground LAI data measured in situ are used to develop algorithms to estimate LAI from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) surface reflectance, and ASTER LAI maps are produced for each ASTER scene using the algorithms developed. Then the ASTER LAI maps are aggregated to compare with the new LAI products. It is found that the variational assimilation method is able to produce temporal continuous LAI products and that accuracy has been improved over the MODIS LAI standard product.

The Leaf Area Index (LAI) is a crucial vegetation parameter for climate and ecological models. Reflectance anisotropy contains valuable supplementary information for the retrieval of properties of an observed target surface. Previous studies have utilized multi-angular reflectance data and physically based Bidirectional Reflectance Distribution Function (BRDF) models with detailed vegetation structure descriptions for LAI estimation. However, the optimal selection of viewing angles for improved inversion results has received limited attention. By optimizing directional observations and integrating the PROSAIL and Ross–Li models, this study aims to enhance LAI estimation from MODIS BRDF data. A dataset of 20,000 vegetation parameter combinations was utilized to identify the directions in which the PROSAIL model exhibits higher sensitivity to LAI changes and better consistency with the Ross–Li BRDF models. The results reveal significant variations in the sensitivity of the PROSAIL model to LAI changes and its consistency with the Ross–Li model over the viewing hemisphere. In the red band, directions with high sensitivity to LAI changes and strong model consistency are mainly found at smaller solar and viewing zenith angles. In the near-infrared band, these directions are distributed at positions with larger solar and viewing zenith angles. Validation using field measurements and LAI maps demonstrates that the proposed method achieves comparable accuracy to an algorithm utilizing 397 viewing angles by utilizing reflectance data from only 30 directions. Moreover, there is a significant improvement in computational efficiency. The accuracy of LAI estimation obtained from simulated multi-angle data is relatively high for LAI values below 3.5 when compared with the MODIS LAI product from two tiles. Additionally, there is also a slight improvement in the results when the LAI exceeds 4.5. Overall, our results highlight the potential of utilizing multi-angular reflectance in specific directions for vegetation parameter inversion, showcasing the promise of this method for large-scale LAI estimation.