High-resolution Leaf Area Index Research Articles

Use of Satellite-based Leaf Area Index Data for Monitoring Green Gram Crop Growth in Ikombe-Katangi Area, Machakos County, Kenya

Monitoring the conditions under which crops grow and estimating their yield is essential to the process of economic development in any nation. The traditional approaches to crop monitoring are labor- and resource-intensive, as well as limited in their ability to cover expansive geographic regions. Since remote sensing data is integrated with ground measurements, it has been used as a potentially useful tool to extract biophysical variables such as leaf area index (LAI), biomass, and phenology. Leaf Area Index (LAI) is a key variable that bridges remote sensing observations to the quantification of agroecosystem processes such as yield. The potential of satellite-based LAI for monitoring crop growth was investigated and compared to field measurements. This was the objective of his study. High resolution leaf area index was retrieved from Landsat 8 OLI imageries, and ground leaf area index measurements were taken in Ikombe-Katangi Machakos in Machakos county. An equation based on regression was developed to estimate the leaf area index using the normalized difference vegetation index that was derived from Landsat-8 OLI (NDVI). According to the findings, the derived LAI exhibited strong linear relationships with the leaf area index that was measured on the ground for green gram crops, with RMSE values being 0.09846, and R2 values of 0.9249. The overall findings shed light on the viability of using multispectral data, in estimating leaf area index in a very fragmented agricultural landscape, such as that found in the Ikombe and Katangi areas of Machakos, Kenya. Therefore, accurate crop monitoring on a large scale can be accomplished through remote sensing data in conjunction with LAI measurements taken from the ground. The implementation of the leaf area index is subject to the environmental condition which requires to be investigated as the use of remote sensing data is advocated for. The advancement in the use of satellite data for extraction LAI data is key towards transforming crop production within Kenya ad Africa at large. Scarce resources for crop growth monitoring can now be augmented with open-source satellite datasets such as Landsat 8 in LAI estimation.

A High Spatiotemporal Enhancement Method of Forest Vegetation Leaf Area Index Based on Landsat8 OLI and GF-1 WFV Data

The leaf area index (LAI) is a crucial parameter for analyzing terrestrial ecosystem carbon cycles and global climate change. Obtaining high spatiotemporal resolution forest stand vegetation LAI products over large areas is essential for an accurate understanding of forest ecosystems. This study takes the northwestern part of the Inner Mongolia Autonomous Region (the northern section of the Greater Khingan Mountains) in northern China as the research area. It also generates the LAI time series product of the 8-day and 30 m forest stand vegetation growth period from 2013 to 2017 (from the 121st to the 305th day of each year). The Simulated Annealing-Back Propagation Neural Network (SA-BPNN) model was used to estimate LAI from Landsat8 OLI, and the multi-period GaoFen-1 WideField-View satellite images (GF-1 WFV) and the spatiotemporal adaptive reflectance fusion mode (STARFM) was used to predict high spatiotemporal resolution LAI by combining inversion LAI and Global LAnd Surface Satellite-derived vegetation LAI (GLASS LAI) products. The results showed the following: (1) The SA-BPNN estimation model has relatively high accuracy, with R2 = 0.75 and RMSE = 0.38 for the 2013 LAI estimation model, and R2 = 0.74 and RMSE = 0.17 for the 2016 LAI estimation model. (2) The fused 30 m LAI product has a good correlation with the LAI verification of the measured sample site (R2 = 0.8775) and a high similarity with the GLASS LAI product. (3) The fused 30 m LAI product has a high similarity with the GLASS LAI product, and compared with the GLASS LAI interannual trend line, it accords with the growth trend of plants in the seasons. This study provides a theoretical and technical reference for forest stand vegetation growth period LAI spatiotemporal fusion research based on high-score data, and has an important role in exploring vegetation primary productivity and carbon cycle changes in the future.

Evaluation of the High-Resolution MuSyQ LAI Product over Heterogeneous Land Surfaces

In recent years, the retrieval and validation of remotely-sensed leaf area index (LAI) products over complex land surfaces have received much attention due to the high-precision land surface model simulations and applications in global climate change. However, most of these related researches mainly focus on coarse resolution products. This is because few products have been specifically designed for solving the problems derived from complex land surfaces in mountain areas until now. MuSyQ LAI is a new product derived from Gaofen-1 (GF-1) satellite data. This product is characterized with a temporal resolution of 10 days and a spatial resolution of 16 m. As is well known, high-resolution products have less uncertainties because of the homogeneities of sub-pixel. Therefore, to evaluate the precision and uncertainty of MuSyQ LAI, an up-scaling strategy was employed here to validate MuSyQ LAI for three mountain regions in Southwest China. The validation strategy can be divided into three parts. First, a regression model was built by in situ LAI measured by LAI-2200 and the normalized difference vegetation index (NDVI) from unmanned aerial vehicle (UAV) images to obtain a 0.5 m resolution LAI map. Second, an up-scaled LAI map with a spatial resolution consistent with MuSyQ LAI was calculated by the pixel-averaging method from the UAV-based LAI map. Third, the MuSyQ LAI was validated by the up-scaled UAV-based LAI in pixel scale. Simultaneously, the sources of uncertainty were analyzed and compared from the view of data source, retrieval model, and scale effects. The results suggested that MuSyQ LAI in the study areas are significantly underestimated by 53.69% due to the complex terrain and heterogeneous land cover. There are three main reasons for the underestimation. The differences between GF-1 reflectance and UAV-based reflectance employed to estimate LAI are the largest factors for the validation results, even accounting for 61.47% of the total bias. Subsequently, the scale effects led to about 28.44% bias. Last but not least, the models employed to retrieve LAI contributed merely 10.09% uncertainties to the total bias. In conclusion, the accuracy of MuSyQ LAI still has a large space to be improved from the view of reflectance over complex terrain. This study is quite important for applications of MuSyQ LAI products and also provides a reference for the improvement and application of other high-resolution remotely sensed LAI products.

A Method for Retrieving Coarse-Resolution Leaf Area Index for Mixed Biomes Using a Mixed-Pixel Correction Factor

The leaf area index (LAI) is a key structural parameter of vegetation canopies. Accordingly, several moderate-resolution global LAI products have been produced and widely used in the field of remote sensing. However, the accuracy of the current moderate-resolution global LAI products cannot satisfy the requirements recommended by the LAI application communities, especially in heterogeneous areas composed of mixed land cover types. In this study, we propose a mixed-pixel correction (MPC) method to improve the accuracy of LAI retrievals over heterogeneous areas by considering the influence of heterogeneity caused by the mixture of different biome types with the help of high-resolution land cover maps. The DART-simulated LAI, the aggregated Landsat LAI, and the site-based high-resolution LAI reference maps are used to evaluate the performance of the MPC method. The results indicate that the MPC method can reduce the influences of spatial heterogeneity and biome misclassification to obtain the LAI with much better accuracy than the Moderate Resolution Imaging Spectroradiometer (MODIS) main algorithm, given that the high-resolution land cover map is accurate. The root mean square error (RMSE) (bias) decreases from 0.749 (0.486) to 0.414 (0.087), while the R2 increases from 0.084 to 0.524, and the proportion of pixels that fulfill the uncertainty requirement of the GCOS increases from 38.2% to 84.6% for the results of site-based high-resolution LAI reference maps. Spatially explicit information about vegetation fractional cover can further reduce uncertainties induced by variations in canopy density for the results of DART simulated data. The proposed method shows potential for improving global moderate-resolution LAI products.

Moderate resolution LAI prediction using Sentinel-2 satellite data and indirect field measurements in Sikkim Himalaya.

The leaf area index (LAI) has been traditionally used as a photosynthetic variable. LAI plays an essential role in forest cover monitoring and has been identified as one of the important climate variables. However, due to challenges in field sampling, complex topography, and availability of cloud-free optical satellite data, LAI assessment on larger scale is still unexplored in the Sikkim Himalayan area. We used two optical instruments, digital hemispherical photography (DHP) and LAI-2200C, to assess the LAI across four different forests following 20 × 20m2 elementary sampling units (ESUs) in the Himalayan state of Sikkim, India. The use of Sentinel-2 derived vegetation indices (VIs) demonstrated a better correlation with the DHP based LAI estimates than using LAI-2200C. Further, the combination of both reflectance bands and VIs were integrated to predict the LAI maps using random forest model. The temperate evergreen forests demonstrated the highest LAI value, while the predicted maps exhibited LAI maxima of 3.4. The estimated vs predicted LAI for DHP and LAI-2200C based estimation demonstrated reasonably good (R2 = 0.63 and R2 = 0.68, respectively) agreement. Further, improvements on the LAI prediction can be attempted by minimizing errors from the inherent field protocols, optimizing the density of field measurements, and representing heterogeneity. The recent rise of frequent forest fires in Sikkim Himalaya prompts for better understanding of fuel load in terms of surface fuel or canopy fuel that can be linked to LAI. The high-resolution LAI map could serve as input to forest fuel bed characterization, especially in seasonal forests with significant variations in green leaves and litter, thereby offering inputs for forest management in changing climate.

CubeSat constellations provide enhanced crop phenology and digital agricultural insights using daily leaf area index retrievals

Satellite remote sensing has great potential to deliver on the promise of a data-driven agricultural revolution, with emerging space-based platforms providing spatiotemporal insights into precision-level attributes such as crop water use, vegetation health and condition and crop response to management practices. Using a harmonized collection of high-resolution Planet CubeSat, Sentinel-2, Landsat-8 and additional coarser resolution imagery from MODIS and VIIRS, we exploit a multi-satellite data fusion and machine learning approach to deliver a radiometrically calibrated and gap-filled time-series of daily leaf area index (LAI) at an unprecedented spatial resolution of 3 m. The insights available from such high-resolution CubeSat-based LAI data are demonstrated through tracking the growth cycle of a maize crop and identifying observable within-field spatial and temporal variations across key phenological stages. Daily LAI retrievals peaked at the tasseling stage, demonstrating their value for fertilizer and irrigation scheduling. An evaluation of satellite-based retrievals against field-measured LAI data collected from both rain-fed and irrigated fields shows high correlation and captures the spatiotemporal development of intra- and inter-field variations. Novel agricultural insights related to individual vegetative and reproductive growth stages were obtained, showcasing the capacity for new high-resolution CubeSat platforms to deliver actionable intelligence for precision agricultural and related applications.

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

Leaf area index (LAI) is a terrestrial essential climate variable that is required in a variety of ecosystem and climate models. The Global LAnd Surface Satellite (GLASS) LAI product has been widely used, but its current version (V5) from Moderate Resolution Imaging Spectroradiometer (MODIS) data has several limitations, such as frequent temporal fluctuation, large data gaps, high dependence on the quality of surface reflectance, and low computational efficiency. To address these issues, this paper presents a deep learning model to generate a new version of the LAI product (V6) at 250-m resolution from MODIS data from 2000 onward. Unlike most existing algorithms that estimate one LAI value at one time for each pixel, this model estimates LAI for 2 years simultaneously. Three widely used LAI products (MODIS C6, GLASS V5, and PROBA-V V1) are used to generate global representative time-series LAI training samples using K-means clustering analysis and least difference criteria. We explore four machine learning models, the general regression neural network (GRNN), long short-term memory (LSTM), gated recurrent unit (GRU), and Bidirectional LSTM (Bi-LSTM), and identify Bi-LSTM as the best model for product generation. This new product is directly validated using 79 high-resolution LAI reference maps from three in situ observation networks. The results show that GLASS V6 LAI achieves higher accuracy, with a root mean square (RMSE) of 0.92 at 250 m and 0.86 at 500 m, while the RMSE is 0.98 for PROBA-V at 300 m, 1.08 for GLASS V5, and 0.95 for MODIS C6 both at 500 m. Spatial and temporal consistency analyses also demonstrate that the GLASS V6 LAI product is more spatiotemporally continuous and has higher quality in terms of presenting more realistic temporal LAI dynamics when the surface reflectance is absent for a long period owing to persistent cloud/aerosol contaminations. The results indicate that the new Bi-LSTM deep learning model runs significantly faster than the GLASS V5 algorithm, avoids the reconstruction of surface reflectance data, and is resistant to the noises (cloud and snow contamination) or missing values contained in surface reflectance than other methods, as the Bi-LSTM can effectively extract information across the entire time series of surface reflectance rather than a single time point. To our knowledge, this is the first global time-series LAI product at the 250-m spatial resolution that is freely available to the public (www.geodata.cn and www.glass.umd.edu).

Multiscale leaf area index assimilation for Moso bamboo forest based on Sentinel-2 and MODIS data

Leaf area index (LAI) is an important driving factor in forest ecosystems carbon cycle, the acquisition of multiscale LAI can illustrate the carbon cycle at different scales and its dynamic response to environmental changes. Moso bamboo forest (MBF) has strong carbon sequestration capability and unique phenological characteristics of on-years and off-years. In this study, based on the MODIS LAI and Sentinel-2 reflectance products in 2018–2019, we coupled the Hierarchical Bayesian Network (HBN) algorithm, LAI dynamic model and PROSAIL model to acquire multiscale spatiotemporal assimilated LAI (500-m, 100-m, 20-m) of MBF. First, the results showed that band7 (B7) and band8a (B8a) of Sentinel-2 were highly sensitive to LAI; the simulated leaf and canopy reflectance by PROSAIL demonstrated higher accuracies and lower errors in these two bands. Second, after applying the Savitzky-Golay smoothing, compared with the observed LAI (Obs_eLAI), the accuracy of MODIS LAI was improved by 121%, and the error was reduced by 24%. Third, there was a significant correlation between the multiscale assimilated LAI (HBN_LAIs) and Obs_eLAI (R500m_HBN_LAI2=0.80,R100m_HBN_LAI2=0.82,R20m_HBN_LAI2=0.80). Finally, the assimilated results in 2018 (off-year) were better than those in 2019 (on-year); the variations in multiscale spatiotemporal assimilated LAI were consistent with the actual growth trends of MBF. This study provides a feasible way to accurately obtain multiscale high-resolution LAI for the carbon cycle simulation of bamboo forests.

A data-driven approach to estimate leaf area index for Landsat images over the contiguous US

Leaf Area Index (LAI) is a fundamental vegetation biophysical variable serving as an essential input to many land surface and atmospheric models. Long-term LAI maps are typically generated with satellite images at moderate spatial resolution (0.25 to 1 km), such as those from the Moderate Resolution Imaging Spectroradiometer (MODIS). While useful for regional-scale land surface modeling, these moderate resolution products often cannot resolve spatial heterogeneity important for many agricultural and hydrological applications. This paper proposes an approach to map LAI at 30-m resolution based on Landsat images for the Contiguous US (CONUS) consistent with the MODIS product, aimed at multi-scale modeling applications. The algorithm was driven by 1.6 million spatially homogeneous samples derived from MODIS LAI and Landsat surface reflectance products from 2006 to 2018. Based on these samples, we trained separate random forest models to estimate LAI from Landsat surface reflectance for eight biomes of the National Land Cover Database (NLCD). A balanced sample design regarding the saturation status of MODIS LAI and a machine-learning-based noise detection technique were introduced to mitigate the trade-off in estimation accuracy between medium LAI (e.g., 3 to 4, unsaturated) and high LAI (e.g., 4–6, saturated).This approach was evaluated using ground measurements from 19 National Ecological Observatory Network (NEON) sites and eight independent sites from other sources. These sites comprise a representative sample of forests, grasslands, shrublands, and croplands across the US. For NEON sites, the LAI estimates show an overall Root Mean Squared Error (RMSE) of 0.8 with r2 of 0.88. For the eight independent sites, the Landsat LAI algorithm achieves RMSE between 0.52 and 0.91. The uncertainty in Landsat estimated LAI varies across biomes and locations. The proposed algorithm was implemented on the Google Earth Engine platform, allowing for the rapid generation of long-term high-resolution LAI records from the 1980s using Landsat images (code is available at https://github.com/yanghuikang/Landsat-LAI). Our findings also highlight the importance of sample balance on regression-based modeling in remote sensing applications.

Time series high-resolution leaf area index estimation and change monitoring in the Saihanba area

Time series high-resolution leaf area index estimation and change monitoring in the Saihanba area

Quantifying Uncertainty and Bridging the Scaling Gap in the Retrieval of Leaf Area Index by Coupling Sentinel-2 and UAV Observations

Leaf area index (LAI) estimates can inform decision-making in crop management. The European Space Agency’s Sentinel-2 satellite, with observations in the red-edge spectral region, can monitor crops globally at sub-field spatial resolutions (10–20 m). However, satellite LAI estimates require calibration with ground measurements. Calibration is challenged by spatial heterogeneity and scale mismatches between field and satellite measurements. Unmanned Aerial Vehicles (UAVs), generating high-resolution (cm-scale) LAI estimates, provide intermediary observations that we use here to characterise uncertainty and reduce spatial scaling discrepancies between Sentinel-2 observations and field surveys. We use a novel UAV multispectral sensor that matches Sentinel-2 spectral bands, flown in conjunction with LAI ground measurements. UAV and field surveys were conducted on multiple dates—coinciding with different wheat growth stages—that corresponded to Sentinel-2 overpasses. We compared chlorophyll red-edge index (CIred-edge) maps, derived from the Sentinel-2 and UAV platforms. We used Gaussian processes regression machine learning to calibrate a UAV model for LAI, based on ground data. Using the UAV LAI, we evaluated a two-stage calibration approach for generating robust LAI estimates from Sentinel-2. The agreement between Sentinel-2 and UAV CIred-edge values increased with growth stage—R2 ranged from 0.32 (stem elongation) to 0.75 (milk development). The CIred-edge variance between the two platforms was more comparable later in the growing season due to a more homogeneous and closed wheat canopy. The single-stage Sentinel-2 LAI calibration (i.e., direct calibration from ground measurements) performed poorly (mean R2 = 0.29, mean NRMSE = 17%) when compared to the two-stage calibration using the UAV data (mean R2 = 0.88, mean NRMSE = 8%). The two-stage approach reduced both errors and biases by >50%. By upscaling ground measurements and providing more representative model training samples, UAV observations provide an effective and viable means of enhancing Sentinel-2 wheat LAI retrievals. We anticipate that our UAV calibration approach to resolving spatial heterogeneity would enhance the retrieval accuracy of LAI and additional biophysical variables for other arable crop types and a broader range of vegetation cover types.

Generating High Resolution LAI Based on a Modified FSDAF Model

Leaf area index (LAI) is an important parameter for monitoring the physical and biological processes of vegetation canopy. Due to the constraints of cloud contamination, snowfall, and instrument conditions, most of the current satellite remote sensing LAI products have lower resolution that cannot satisfy the needs of vegetation remote sensing application in areas of high heterogeneity. We proposed a new model to generate high resolution LAI, by combining linear pixel unmixing and the Flexible Spatiotemporal Data Fusion (FSDAF) method. This method derived the input data of FSDAF by downscaling the MODIS (Moderate Resolution Imaging Spectroradiometer) data with a linear spectral mixture model. Through the improved input parameters of the algorithm, the fusion of MODIS LAI and LAI at Landsat spatial resolution estimated by Support Vector Regression model was realized. The fusion accuracy of generated LAI data was validated based on Sentinel-2 LAI products. The results showed that strong correlation between predicted LAI and Sentinel-2 LAI in sample sites was observed with higher correlation coefficients and lower Root Mean Square Error. Compared to the simulation results of FSDAF, the modified FSDAF model showed higher accuracy and reflected more spatial details in the boundary areas of different land cover types.

Deriving high-spatiotemporal-resolution leaf area index for agroecosystems in the U.S. Corn Belt using Planet Labs CubeSat and STAIR fusion data

Leaf area index (LAI) is a key variable for characterizing crop growth conditions and estimating crop productivity. Despite continuing efforts to develop LAI estimation algorithms, LAI datasets still need improvement at spatial and temporal resolutions to meet the requirements of agricultural applications. Advancements in data fusion technique and the emergence of new satellite data provide opportunities for LAI data at higher resolutions in both space and time. In this study, we derived new LAI estimations by leveraging novel satellite remote sensing datasets, STAIR fusion (MODIS-Landsat fusion) and Planet Labs' CubeSat data (through a reprocessed pipeline) for a typical agricultural landscape in the U.S. Corn Belt. The STAIR fused data and our reprocessed CubeSat data have both fine spatial resolutions (30 m and 3.125 m, respectively) and high frequencies (daily for both). To reliably estimate LAI from these advanced satellite datasets, we used two methods: inversion of a radiative transfer model (RTM), and empirical relationship with vegetation index (VI) calibrated from field measured LAI. Compared to the ground-truth LAI collected at 36 sites across the study region, reliable approximations were achieved by both LAI estimations based on PROSAIL RTM (STAIR: R2 = 0.69 and root mean squared error (RMSE) = 1.12 (m2 m−2), CubeSat: R2 = 0.76 and RMSE = 1.09 (m2 m−2)), and LAI estimations based on Green Wide Dynamic Range Vegetation Index (GrWDRVI) (STAIR: R2 = 0.75, RMSE = 1.10 (m2 m−2), CubeSat: R2 = 0.76, RMSE = 1.08 (m2 m−2), where validation ground-truth is independent from calibration data). Newly estimated high-resolution LAI data were aggregated at 500 m resolution and compared with MODIS and VIIRS LAI products, revealing substantial uncertainties and biases in these two products. We also demonstrated phenology stage estimation at fine spatial resolutions based on our high-frequency LAI data. The proposed LAI estimation methods at both high spatial resolution and temporal frequency can be applied to the entire U.S. Corn Belt and provide significant advancement to crop monitoring and precision agriculture.

A Data Assimilation Method for Simultaneously Estimating the Multiscale Leaf Area Index From Time-Series Multi-Resolution Satellite Observations

Current global leaf area index (LAI) products are generally produced from single-temporal satellite observations acquired by a single sensor. These LAI products are usually spatiotemporally discontinuous and inaccurate for some vegetation types in many areas, which limit the applications of these LAI products to the understanding of land dynamics. In this paper, a new data assimilation method was proposed to estimate multiscale and temporally continuous LAI values from multi-sensor time-series satellite observations with different spatial resolutions. An ensemble multiscale tree (EnMsT) was used to establish the conversion relationships between different spatial resolution LAI values, and dynamic models of the LAI at different spatial scales were constructed to evolve LAI at the corresponding spatial scales over time. At each time step, a multiscale Kalman filter (MKF) was introduced to fuse the predicted LAI values from the dynamic models at different spatial scales and to construct a forecasted EnMsT. When satellite observations were available, an ensemble multiscale filter (EnMsF) technique was applied to update the LAI values at each node of the EnMsT. The method was applied to estimate temporally continuous multiscale LAI values from the time series of Thematic Mapper (TM) or Enhanced Thematic Mapper Plus (ETM+) surface reflectance data and Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data at several sites with different vegetation types. The estimated multiscale LAI values were compared with the MODIS and GEOV2 LAI products, and the reference LAI values at the corresponding scales aggregated from the high-resolution LAI surface images. The estimated LAI values with the finest spatial resolution were also validated by ground measurements from the selected sites. The results show that the new method is able to simultaneously estimate temporally continuous multiscale LAI values by assimilating satellite observations with different spatial resolutions, and the estimated multiscale LAI values are well consistent with the reference LAI values at the corresponding scales over the selected sites. The root-mean-square error (RMSE) and coefficient of determination of the retrieved LAI values at the finest spatial scale against the ground measurements over the selected sites are 0.539 and 0.788, respectively.

High resolution mapping of agricultural water productivity using SEBAL in a cultivated African catchment, Tanzania

The application of remote sensing techniques for WPET mapping in data scarce regions is gaining more recognition since it can cover large areas with minimal field observations. Important concerns are the generation of high-resolution WPET maps and addressing the question on how accurate the results are. This study aims at high resolution (10 m) mapping and evaluation of the spatial variability of biomass, yield, ET and WPET in the Makanya river catchment using the automated Surface Energy Balance Algorithm for Land (pySEBAL) with SENTINEL-2 and LANDSAT-8 images, local land use map and locally calibrated leaf area index (LAI) inputs. A coupled phenological variability and supervised classification approach on high resolution images generated a high accuracy LULC layer which was used to map the WPET in the agricultural lands. The pySEBAL results were evaluated in view of local information on crop yields, water allocation and agricultural management practices in the different agro-ecological zones within the catchment. Calibration of high-resolution satellite LAI generated products with error estimates within acceptable levels of uncertainty. The simulated crop yields were in agreement with reported crop yields. The results showed relatively high WPET in the highlands and low WPET in the midland and lowland areas of the catchment. The latter was attributed to high transmission losses, low irrigation efficiencies, poor agricultural practices and pest/disease attack. When applying SEBAL in African cultivated catchments, it is highly recommended to use SENTINEL-2 data in addition to LANDSAT-8, and to use local information, especially for the ground truthing of land use maps, phenology, crop practices and crop yields.

A Recursive Update Model for Estimating High-Resolution LAI Based on the NARX Neural Network and MODIS Times Series

Leaf area index (LAI) remote sensing data products with a high resolution (HR) and long time series are in demand in a wide variety of applications. Compared with long time series LAI products with 1 km resolution, LAI products with high spatial resolution are difficult to acquire because of the lack of remote sensing observations in long-term sequences and the lack of estimation methods applicable to highly variable land-cover types. To address these problems, we proposed a recursive update model to estimate 30 m resolution LAI based on the updated Nonlinear Auto-Regressive with Exogenous Inputs (NARX) neural network and MODIS time series. First, we used a variety of HR satellite remote sensing observations to produce HR datasets for recent years. Historical low spatial resolution MODIS products were employed as background information and used to calculate the initial parameters of the NARX neural network for each pixel. Subsequently, one year’s reflectance from the HR dataset was used as the new observation that was input into the NARX model to estimate the HR LAI of that year, and the background and HR data were then used for remodeling to update the NARX model parameters. This procedure was recursively repeated year by year until both MODIS background data and all HR data were involved in the modeling. Finally, we obtained an LAI time series with 30 m resolution. In the cropland study area in Hebei Province, China, the results were compared with LAI measurements from ground sites in 2013 and 2014. A high degree of similarity existed between the results for the two study years ( RMSE 2013 = 0.288 and RMSE 2014 = 0.296 ). The HR LAI estimates showed favorable spatiotemporal continuity and were in good agreement with the multisample ground survey LAI measurements. The results indicated that for data with a rapid revisit cycle and high spatial resolution, the recursive update model based on the NARX neural network has excellent LAI estimation performance and fairly strong fault-tolerance capability.

The Retrieval of 30-m Resolution LAI from Landsat Data by Combining MODIS Products

Leaf area index (LAI) is a critical vegetation structural parameter in biogeochemical and biophysical ecosystems. High-resolution LAI products play an essential role in regional studies. Empirical methods, which normally use field measurements as their training samples and have been identified as the most commonly used approaches to retrieve structural parameters of vegetation from high-resolution remote-sensing data, are limited by the quality of training samples. Few efforts have been made to generate training samples from existing global LAI products. In this study, two methods (a homogeneous and pure pixel filter method (method A) and a pixel unmixing method (method B)) were developed to extract training samples from moderate-resolution imaging spectroradiometer (MODIS) surface reflectance and LAI products, and a support vector regression (SVR) algorithm trained by the samples was used to retrieve the high-resolution LAI from Landsat data at Baoding, situated in the Hebei Province in China, and Des Moines, situated in Iowa, United States. For the homogeneous and pure pixel filter method, two different sets of training samples were designed. One was composed of upscaled Landsat reflectance at the 500-m resolution and MODIS LAI products (dataset A1); the other was composed of MODIS reflectance and LAI products (dataset A2). With them, two inversion models were developed using SVR. For the pixel unmixing method, the training samples (dataset B) were extracted from unmixed MODIS surface reflectance and LAI products at 30-m resolution, and the third inversion model was obtained with them. LAI inversion results showed that good agreement with field measurements was achieved using these three inversion models. The R2 (coefficient of determination) value and the root mean square error (RMSE) value were computed to assess the results. For all tests, the R2 values are higher than 0.74 and RMSE values are less than 0.73. These tests showed that three models for the two methods combined with MODIS products can retrieve 30-m resolution LAI from Landsat data. The results of the pixel unmixing method was slightly better than that of the homogeneous and pure pixel filter method.

A Critical Comparison of Remote Sensing Leaf Area Index Estimates over Rice-Cultivated Areas: From Sentinel-2 and Landsat-7/8 to MODIS, GEOV1 and EUMETSAT Polar System

Leaf area index (LAI) is a key biophysical variable fundamental in natural vegetation and agricultural land monitoring and modelling studies. This paper is aimed at comparing, validating and discussing different LAI satellite products from operational services and customized solution based on innovative Earth Observation (EO) data such as Landsat-7/8 and Sentinel-2A. The comparison was performed to assess overall quality of LAI estimates for rice, as a fundamental input of different scale (regional to local) operational crop monitoring systems such as the ones developed during the “An Earth obseRvation Model based RicE information Service” (ERMES) project. We adopted a multiscale approach following international recognized protocols of the Committee on Earth Observation Satellites (CEOS) Land Product Validation (LPV) guidelines in different steps: (1) acquisition of representative field sample measurements, (2) validation of decametric satellite product (10–30 m spatial resolution), and (3) exploitation of such data to assess quality of medium-resolution operational products (~1000 m). The study areas were located in the main European rice areas in Spain, Italy and Greece. Field campaigns were conducted during three entire rice seasons (2014, 2015 and 2016—from sowing to full-flowering) to acquire multi-temporal ground LAI measurements and to assess Landsat-7/8 LAI estimates. Results highlighted good correspondence between Landsat-7/8 LAI estimates and ground measurements revealing high correlations (R2 ≥ 0.89) and low root mean squared errors (RMSE ≤ 0.75) in all seasons. Landsat-7/8 as well as Sentinel-2A high-resolution LAI retrievals, were compared with satellite LAI products operationally derived from MODIS (MOD15A2), Copernicus PROBA-V (GEOV1), and the recent EUMETSAT Polar System (EPS) LAI product. Good agreement was observed between high- and medium-resolution LAI estimates. In particular, the EPS LAI product was the most correlated product with both Landsat/7-8 and Sentinel-2A estimates, revealing R2 ≥ 0.93 and RMSE ≤ 0.53 m2/m2. In addition, a comparison exercise of EPS, GEOV1 and MODIS revealed high correlations (R2 ≥ 0.90) and RMSE ≤ 0.80 m2/m2 in all cases and years. The temporal assessment shows that the three satellite products capture well the seasonality during the crop phenological cycle. Discrepancies are observed mainly in absolute values retrieved for the peak of rice season. This is the first study that provides a quantitative assessment on the quality of available operational LAI product for rice monitoring to both the scientific community and users of agro-monitoring operational services.

A Method to Upscale the Leaf Area Index (LAI) Using GF-1 Data with the Assistance of MODIS Products in the Poyang Lake Watershed

A leaf area index is a key parameter reflecting the growth changes of vegetation and one of the most important canopy structural parameters for performing quantitative analyses of many ecological and climate models. Although using high-resolution satellite data and the radiative transfer model (RTM) can be used to generate high resolution LAI products, the RTM method has some problems because its temporal resolution is low, the input parameters are more appropriate for a physics model, and some parameters are difficult to obtain. Problems that urgently need to be solved include improving the temporal-spatial resolution for LAI products and localizing LAI products. To explore an applicable method for the high-resolution LAI products in a small basin and to improve the inversion accuracy, we propose an approach for GF-1 WFV LAI retrieval using MOD15A2 data and the measured LAI of the Poyang Lake watershed. Empirical models were used to retrieve high resolution LAI values, and the results show that these models are well designed for analyzing time-series satellite data. Good correlations were obtained between the NDVI of the GF-1 WFV data, the retrieved LAI values and the MODIS LAI data from samples acquired in both summer and winter. The exponential NDVI model obtained the best LAI value estimation results from the GF-1 WFV data (R2 = 0.697, RMSE = 1.100); the best synthetic validation of the RMSE is 0.883, close to the optimum model. Therefore, the retrieval results more fully reflect the growth process of the different features. This study proposed an upscale method for developing a high spatial resolution GF-1 satellite standard LAI products retrieval model using MODIS data. The proposed method will be helpful for efficiently improving the temporal-spatial resolution of LAI products to benefit the extraction of vegetation parameter information and dynamic land use monitoring.

Leaf Area Index Estimation Using Chinese GF-1 Wide Field View Data in an Agriculture Region.

Leaf area index (LAI) is an important vegetation parameter that characterizes leaf density and canopy structure, and plays an important role in global change study, land surface process simulation and agriculture monitoring. The wide field view (WFV) sensor on board the Chinese GF-1 satellite can acquire multi-spectral data with decametric spatial resolution, high temporal resolution and wide coverage, which are valuable data sources for dynamic monitoring of LAI. Therefore, an automatic LAI estimation algorithm for GF-1 WFV data was developed based on the radiative transfer model and LAI estimation accuracy of the developed algorithm was assessed in an agriculture region with maize as the dominated crop type. The radiative transfer model was firstly used to simulate the physical relationship between canopy reflectance and LAI under different soil and vegetation conditions, and then the training sample dataset was formed. Then, neural networks (NNs) were used to develop the LAI estimation algorithm using the training sample dataset. Green, red and near-infrared band reflectances of GF-1 WFV data were used as the input variables of the NNs, as well as the corresponding LAI was the output variable. The validation results using field LAI measurements in the agriculture region indicated that the LAI estimation algorithm could achieve satisfactory results (such as R2 = 0.818, RMSE = 0.50). In addition, the developed LAI estimation algorithm had potential to operationally generate LAI datasets using GF-1 WFV land surface reflectance data, which could provide high spatial and temporal resolution LAI data for agriculture, ecosystem and environmental management researches.