Fine-resolution Leaf Area Index Research Articles

Three-dimensional estimation of deciduous forest canopy structure and leaf area using multi-directional, leaf-on and leaf-off airborne lidar data

Airborne laser scanning (ALS) has been widely used to map gap probability and leaf area index (LAI) distribution at plot and landscape scales. As an indirect measurement, most ALS methods to estimate LAI combine waveform or point density information with supporting field measurements such as the leaf angle distribution, gap probability, or direct LAI measures. The development of a more independent estimation approach would facilitate more widespread use of existing ALS data to investigate patterns of forest structure and build realistic 3-D vegetation scenes to simulate remote sensing imagery and energy balance. Here, we develop a data processing workflow (named PVlad) using ALS point cloud apparent reflectance to estimate LAI and voxel-based leaf area density (LAD), aiming to reduce the need for associated field measurements such as the gap probability. The adaptation of the path volume (PV) concept derived from apparent reflectance integrates information from multi-directional ALS pulses, and quantifies the percentage exploration of each voxel for classification and occlusion correction, such that rigorous volumetric sampling approaches can be developed to derive LAI and LAD. The PVlad workflow was applied to discrete-return lidar data (Riegl VQ480i) acquired by NASA Goddard's LiDAR, Hyperspectral and Thermal Imager (G-LiHT) Airborne Imager during leaf-on (summer) and leaf-off (spring) conditions at the Smithsonian Environmental Research Center (SERC). The estimates of LAI and LAD captured structural differences between mature, logged, and intermediate-aged stands over eight deciduous forest plots. The derived LAI values were compared to field litter collection measurements, and the derived LAD vertical distribution was compared to the output of the VoxLAD model using terrestrial laser scan (TLS) field survey data. Using voxel sizes ranging from 0.5 m to 5 m, overall LAI estimation showed linear fitting coefficient bias <0.035 and RMSE<0.5m2/m2 for 1 and 2 m voxel sizes, and vertical LAD distribution showed strong correlation with R≈0.9 and RMSE≈0.028m2/m3 for 0.5 and 1m voxel sizes. For every forest stand, upper-canopy LAD had a low variance for voxel sizes of ≤ 2m. Application of PVlad to the G-LiHT and other similar ALS data archives enables the development of fine-resolution LAI map products, including voxelization of LAD for ecosystem science and radiative transfer simulations of remote sensing imagery or surface energy balance.

ValLAI_Crop, a validation dataset for coarse-resolution satellite LAI products over Chinese cropland

Numerous validation efforts have been conducted over the last decade to assess the accuracy of global leaf area index (LAI) products. However, such efforts continue to face obstacles due to the lack of sufficient high-quality field measurements. In this study, a fine-resolution LAI dataset consisting of 80 reference maps was generated during 2003–2017. The direct destructive method was used to measure the field LAI, and fine-resolution LAI images were derived from Landsat images using semiempirical inversion models. Eighty reference LAI maps, each with an area of 3 km × 3 km and a percentage of cropland larger than 75%, were selected as the fine-resolution validation dataset. The uncertainty associated with the spatial scale effect was also provided. Ultimately, the fine-resolution reference LAI dataset was used to validate the Moderate Resolution Imaging Spectroradiometer (MODIS) LAI product. The results indicate that the fine-resolution reference LAI dataset builds a bridge to link small sampling plots and coarse-resolution pixels, which is extremely important in validating coarse-resolution LAI products.

A Model-Downscaling Method for Fine-Resolution LAI Estimation

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.

Spatially and Temporally Continuous Leaf Area Index Mapping for Crops through Assimilation of Multi-resolution Satellite Data

As a key parameter that represents the structural characteristics and biophysical changes of crop canopy, the leaf area index (LAI) plays a significant role in monitoring crop growth and mapping yield. A considerable amount of farmland is dispersed with strong spatial heterogeneity. The existing time series satellite LAI products fail to capture spatial distributions and growth changes of crops due to coarse spatial resolutions and spatio-temporal discontinuities. Therefore, it becomes crucial for fine resolution LAI mapping in time series over crop areas. A two-stage data assimilation scheme was developed for dense time series LAI mapping in this study. A LAI dynamic model was first constructed using multi-year MODIS LAI data. This model coupled with the PROSAIL radiative transfer model, and MOD09A1 reflectance data were used to retrieve temporal LAI profiles at the 500 m resolution with the assistance of the very fast simulated annealing (VFSA) algorithm. Then, the LAI dynamics at the 500 m scale were incorporated as prior information into the Landsat 8 OLI reflectance data for time series LAI mapping at the 30 m resolution. Finally, the spatio-temporal continuities and retrieval accuracies of assimilated LAI values were assessed at the 500 m and 30 m resolutions respectively, using the MODIS LAI product, fine resolution LAI reference map and field measurements. The results indicated that the assimilated the LAI estimations at the 500 m scale effectively eliminated the spatio-temporal discontinuities of the MODIS LAI product and displayed reasonable temporal profiles and spatial integrity of LAI. Moreover, the 30 m resolution LAI retrievals showed more abundant spatial details and reasonable temporal profiles than the counterparts at the 500 m scale. The determination coefficient R2 between the estimated and field LAI values was 0.76 with a root mean square error (RMSE) value of 0.71 at the 30 m scale. The developed method not only improves the spatio-temporal continuities of the LAI at the 500 m scale, but also obtains 30 m resolution LAI maps with fine spatial and temporal consistencies, which can be expected to meet the needs of analysis on crop dynamic changes and yield mapping in fragmented and highly heterogeneous areas.

A Multiscale Assimilation Approach to Improve Fine-Resolution Leaf Area Index Dynamics

Fine spatial details of vegetation growth are usually lost in leaf area index (LAI) products obtained from coarse spatial resolution satellite sensors. This may bring uncertainties in ecosystem process models, which usually require LAI products with fine spatiotemporal resolutions. Successful downscaling of LAI dynamics to fine spatial resolution is very important for meeting the demands of these models. Hence, a multiscale multisensor approach using the ensemble Kalman smoother (EnKS) technique is proposed in this paper. The LAI dynamics at a coarser spatial resolution are incorporated as prior information into the remotely sensed observations for time series LAI estimation at a finer spatial resolution. Downscaled LAI dynamics are evaluated based on spatial distribution and temporal trajectory. The results indicate the assimilated LAI to be in good agreement with the reference values at the different spatial scales. For example, the coefficient of determination ( $R^{2}$ ) between the reference values and fine-resolution LAI results retrieved by the proposed approach is 0.71 with a root-mean-square-error (RMSE) value of 0.65 on Julian day 185 at the Agro site. The method has proved to be effective for downscaling LAI dynamics, which improves the spatiotemporal patterns of fine-resolution LAI retrievals with respect to earlier methods.

Mapping daily leaf area index at 30 m resolution over a meadow steppe area by fusing Landsat, Sentinel-2A and MODIS data

ABSTRACTThe leaf area index (LAI) is a key vegetation canopy structure parameter and is closely associated with vegetation photosynthesis, transpiration, and energy balance. Developing a landscape-scale LAI dataset with a high temporal resolution (daily) is essential for capturing rapidly changing vegetation structure at field scales and supporting regional biophysical modeling efforts. In this study, two daily 30 m LAI time series from 2014 to 2016 over a meadow steppe site in northern China were generated using a spatial and temporal adaptive reflectance fusion model (STARFM) combined with an LAI retrieval radiative transfer model (PROSAIL). Gap-filled Landsat 7, Landsat 8 and Sentinel-2A surface reflectance (SR) images were used to generate fine-resolution LAI maps with the PROSAIL look-up table method. Two daily 500 m moderate-resolution imaging spectroradiometer (MODIS) LAI product-the existing MCD15A3H LAI product and one was generated from the MCD43A4 SR product and the PROSAIL model, were used to provide temporally continuous LAI variations. The STARFM model was then used to fuse the fine-resolution LAI maps with the two 500 m LAI products separately to generate two daily 30 m LAI time series. Both results were assessed for three types of pasture (mowed pasture, grazing pasture, and fenced pasture) using ground measurements from 2014–2015. The results showed that the PROSAIL-generated LAI maps all exhibited a high accuracy, and the root mean squared errors (RMSEs) for the Landsat 7 LAI and Landsat 8 LAI compared to the ground-measured LAI were 0.33 and 0.28 respectively. The Landsat LAI maps also showed good agreement and similar spatial patterns with the Sentinel-2A LAI with mean differences between ± 0.5. The MCD43A4_PROSPECT LAI product exhibited similar seasonal variability to the ground measurements and to the Landsat and Sentinel-2A LAIs, and these data are also smoother and contain fewer noisy points than the gap-filled MCD15A3H LAI product. Compared to the ground measurements, the daily 30 m LAI time series fused from the fine-resolution LAI maps and PROSPECT generated MODIS LAI product demonstrated better performance with an RMSE of 0.44 and a mean absolute error (MAE) of 0.34, which is an improvement from the LAI time series fused from the fine-resolution LAI maps and the existing MCD15A3H LAI product (RMSE of 0.56 and MAE of 0.42). The latter dataset also exhibited abnormal temporal fluctuations, which may have been caused by the interpolation method. The results also demonstrated the very good performance of the STARFM model in grazing and mowed pasture with homogeneous surfaces compared to fenced pasture with smaller patch sizes. The Sentinel-2A data offers increased landscape vegetation observation frequency and provides temporal information about canopy changes that occur between Landsat overpass dates. The scheme developed in this study can be used as a reference for regional vegetation dynamic studies and can be applied to larger areas to improve grassland modeling efforts.

Integration of Landsat ETM+ Data with Field Measurements for Mapping Leaf Area Index in the Grasslands of Central Kazakhstan

A fine-resolution leaf area index (LAI) data set over a 150 km × 150 km region in central Kazakhstan is retrieved using Landsat ETM+ imagery and ground-based LAI inferred from hemispherical photography and direct measurements. Regression analysis and geostatistics are applied for developing empirical models of LAI from Landsat ETM+ data. The best accuracy is achieved using a model employing a canonical index that combines all the contributions of individual Landsat ETM+ bands into a single index (R 2 = 0.67; RMSE = 0.21). This model is then applied for mapping LAI at a regional scale.

Mapping Spatial and Temporal Variations of Leaf Area Index for Winter Wheat in North China

Leaf area index (LAI) is an important parameter in a number of models related to ecosystem functioning, carbon budgets, climate, hydrology, and crop growth simulation. Mapping and monitoring the spatial and temporal variations of LAI are necessary for understanding crop growth and development at regional level. In this study, the relationships between LAI of winter wheat and Landsat TM spectral vegetation indices (SVIs) were analyzed by using the curve estimation procedure in North China Plain. The series of LAI maps retrieved by the best regression model were used to assess the spatial and temporal variations of winter wheat LAI. The results indicated that the general relationships between LAI and SVIs were curvilinear, and that the exponential model gave a better fit than the linear model or other nonlinear models for most SVIs. The best regression model was constructed using an exponential model between surface-reflectance-derived difference vegetation index (DVI) and LAI, with the adjusted R 2 (0.82) and the RMSE (0.77). The TM LAI maps retrieved from DVI-LAI model showed the significant spatial and temporal variations. The mean TM LAI value (30 m) for winter wheat of the study area increased from 1.29 (March 7, 2004) to 3.43 (April 8, 2004), with standard deviations of 0.22 and 1.17, respectively. In conclusion, spectral vegetation indices from multi-temporal Landsat TM images can be used to produce fine-resolution LAI maps for winter wheat in North China Plain.