Spatially and temporally continuous LAI data sets based on an integrated filtering method: Examples from North America

Leaf Area Index (LAI) is an important biophysical variable for characterizing the land surface vegetation. Global LAI product has been routinely produced from the MODerate resolution Imaging Spectroradiometer (MODIS) aboard the Terra and Aqua satellite platforms. However, the MODIS standard LAI product is not continuous both spatially and temporally. To fill the gaps and improve the quality, we have developed a data filtering algorithm. This filter, called the temporal spatial filter (TSF), integrates both spatial and temporal characteristics for different plant functional types. The spatial gaps are first filled with the multi-year averages of the same day. If the values are missing over all years, the pixel is filled with a new estimate using the vegetation continuous field-ecosystem curve fitting method. The TSF integrates both the multi-seasonal average trend (background) and the seasonal observation. We implement this algorithm using the MODIS Collection 4 LAI product over North America. Comparison of the TSF results with the Savitzky-Golay filter indicates that the TSF performs much better in restoring the spatial and temporal distribution of seasonal LAI trends. The new LAI product has been validated by comparing with field measurements and the derived LAI maps from ETM+ data at a broadleaf forest site and an agricultural site. The validation results indicate that the new LAI product agrees better with both the field measurements and LAI values obtained from the ETM+ than does the MODIS LAI standard product, which usually shows higher LAI values.     

Estimating time-series leaf area index based on recurrent nonlinear autoregressive neural networks with exogenous inputs

The leaf area index (LAI) is a key parameter in many meteorological, environmental and agricultural models. At present, global LAI products from several sensors have been released. These single sensor-based LAI products are generally discontinuous in time and cannot characterize the status of natural vegetation growth very well. In this study, by fusing Moderate Resolution Imaging Spectroradiometer (MODIS) and Satellite Pour l'Observation de la Terre (SPOT) VEGETATION products, time-series LAIs were used to train recurrent nonlinear autoregressive neural networks with exogenous inputs (NARXNNs) for six typical vegetation types. The exogenous inputs included time-series reflectances in the red, near-infrared and shortwave infrared bands as well as the corresponding sun-viewing angles. These NARXNNs subsequently served to predict the time-series LAI. The validation results show that the predicted LAI of the NARXNN is not only more continuous and stable than the MODIS LAI as a function of time but is also much closer to the ground truth. Thus, the proposed method may be helpful for improving the quality of the MODIS LAI.     

Comparison and validation of MODIS and VEGETATION global LAI products over four BigFoot sites in North America

A new set of recently developed leaf area index (LAI) algorithms has been employed for producing a global LAI dataset at 1 km resolution and in time-steps of 10 days, using data from the Satellite pour l'observation de la terre (SPOT) VEGETATION (VGT) sensor. In this paper, this new LAI product is compared with the global MODIS Collection 4 LAI product over four validation sites in North America. The accuracy of both LAI products is assessed against seven high resolution ETM+ LAI maps derived from field measurements in 2000, 2001, and 2003. Both products were closely matched outside growing season. The MODIS product tended to be more variable than the VGT product during the summer period when the LAI was maximum. VGT and ETM+ LAI maps agreed well at three out of the four sites. The median relative absolute error of the VGT LAI product varied from 24% to 75% at 1 km scale and it ranged from 34% to 88% for the MODIS LAI product. The importance of correcting field measurements for the clumping effect is illustrated at the deciduous broadleaf forest site (HARV). Inclusion of the sub-pixel land cover information improved the quality of LAI estimates for the prairie grassland KONZ site. Further improvement of the global VGT LAI product is suggested by production and inclusion of pixel-specific global foliage clumping index and forest background reflectance maps that would serve as an input into the VGT LAI algorithms.     

Comparisons of land cover and LAI estimates derived from ETM+ and MODIS for four sites in North America: a quality assessment of 2000/2001 provisional MODIS products

The MODIS land science team produces a number of standard products, including land cover and leaf area index (LAI). Critical to the success of MODIS and other sensor products is an independent evaluation of product quality. In that context, we describe a study using field data and Landsat ETM+ to map land cover and LAI at four 49-km² sites in North America containing agricultural cropland (AGRO), prairie grassland (KONZ), boreal needleleaf forest, and temperate mixed forest. The purpose was to: (1) develop accurate maps of land cover, based on the MODIS IGBP (International Geosphere-Biosphere Programme) land cover classification scheme; (2) derive continuous surfaces of LAI that capture the mean and variability of the LAI field measurements; and (3) conduct initial MODIS validation exercises to assess the quality of early (i.e., provisional) MODIS products. ETM+ land cover maps varied in overall accuracy from 81% to 95%. The boreal forest was the most spatially complex, had the greatest number of classes, and the lowest accuracy. The intensive agricultural cropland had the simplest spatial structure, the least number of classes, and the highest overall accuracy. At each site, mapped LAI patterns generally followed patterns of land cover across the site. Predicted versus observed LAI indicated a high degree of correspondence between field-based measures and ETM+ predictions of LAI. Direct comparisons of ETM+ land cover maps with Collection 3 MODIS cover maps revealed several important distinctions and similarities. One obvious difference was associated with image/map resolution. ETM+ captured much of the spatial complexity of land cover at the sites. In contrast, the relatively coarse resolution of MODIS did not allow for that level of spatial detail. Over the extent of all sites, the greatest difference was an overprediction by MODIS of evergreen needleleaf forest cover at the boreal forest site, which consisted largely of open shrubland, woody savanna, and savanna. At the agricultural, temperate mixed forest, and prairie grassland sites, ETM+ and MODIS cover estimates were similar. Collection 3 MODIS-based LAI estimates were considerably higher (up to 4 m² m⁻²) than those based on ETM+ LAI at each site. There are numerous probable reasons for this, the most important being the algorithms' sensitivity to MODIS reflectance calibration, its use of a prelaunch AVHRR-based land cover map, and its apparent reliance on mainly red and near-IR reflectance. Samples of Collection 4 LAI products were examined and found to consist of significantly improved LAI predictions for KONZ, and to some extent for AGRO, but not for the other two sites. In this study, we demonstrate that MODIS reflectance data are highly correlated with LAI across three study sites, with relationships increasing in strength from 500 to 1000 m spatial resolution, when shortwave-infrared bands are included.     

Variational retrieval of leaf area index from MODIS time series data: examples from the Heihe river basin,north-west China

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.     

中国西南山区GEOV１、GLASS和MODISLAI产品的对比分析

叶面积指数(LeafAreaIndex,LAI)是表征植被生物物理变化和冠层结构特征的关键参数,目前存在多个全球范围、长时间序列LAI产品,对其进行验证是LAI产品应用的重要前提,然而目前山区的验证工作尤其少见.在我国西南山区选取６个典型样区,考虑山区复杂地形特征,从产品时空完整性以及对山区植被时空特征表征能力等方面对GEOV１、GLASS和MODISLAI产品进行对比分析.研究结果表明:①相比于地形平坦地区,在山区随海拔和地形起伏度的增加,LAI产品时空完整性呈递减的趋势,其中,GEOV１LAI表现最差,MODISLAI次之,GLASSLAI表现最好;②GLASSLAI和GEOV１LAI的空间分布合理且具有较好的一致性,MODISLAI的空间分布和二者存在差异,３种LAI产品均难以准确反映山区植被垂直带谱的变化特征;③草地类型LAI产品间差值较小,林地和农作物GLASSLAI和GEOV１LAI产品一致性较好,MODISLAI产品和二者存在较大的差异;④GLASSLAI时间序列曲线平滑且连续,GEOV１LAI存在时间不连续现象,MODISLAI季相变化中的波动现象比较严重;各产品不仅难以准确反映冬季的常绿针叶林LAI,而且难以准确表征样区内农田作物轮作的物候信息.对比分析有助于发现LAI产品在山区存在的问题,并为今后LAI产品的算法改进提供帮助和参考.     

黑河及汉江流域MODIS叶面积指数产品质量评价

叶面积指数（LAI）是MODIS地面队伍生产的一系列标准产品之一，对其进行独立的质量评价有助于用户了解数据的适用性。本文用近同时相的Landsat影像及野外实测LAI数据获得了黑河及汉江两个研究区高分辨率的Landsat LAl分布图。基于此，对MODIS LAI数据进行了质量评价，评价指标包括统计特征和空间特征。分析结果表明，就统计特征而言，MODIS LAI数据值一般低于Landsat的LAI值，在植被覆盖较好的汉江区低估约10％，在植被覆盖稀疏的黑河区，LAI值低估达58％；就空间特征而言，两个研究区的结果都表明MODIS LAI数据无法很好地体现植被空间分布信息，在黑河区存在大量低槽被覆盖像元被归类为非植被覆盖区的情况。     

呼伦贝尔草甸草原MODIS/LAI产品验证

以内蒙古呼伦贝尔草甸草原为研究区域,利用2013年6期地面实测数据,结合HJ-1A/B CCD高分辨率影像,经过辐射校正与模型建立,对研究区域草原生长季的MODIS/LAI产品进行验证。结果表明:在时间上,MODIS/LAI产品能够较好地反映草原的长势与物候变化。在空间上,由于MODIS/LAI产品输入数据的不确定性,MODIS/LAI产品与地面情况存在一定偏差(ΔLAI=0.59m2/m2),在呼伦贝尔草甸草原草场整个生长季都存在高估现象,平均相对误差为40%。在生长初期和末期,较大的地表异质性使MODIS/LAI产品高估现象较严重;生长中期高估现象减小,相对误差在30%以内。研究结果对了解该地区的MODIS/LAI产品精度与使用该地区MODIS/LAI产品具有重要指导意义。     

A multisensor fusion approach to improve LAI time series

High-quality and gap-free satellite time series are required for reliable terrestrial monitoring. Moderate resolution sensors provide continuous observations at global scale for monitoring spatial and temporal variations of land surface characteristics. However, the full potential of remote sensing systems is often hampered by poor quality or missing data caused by clouds, aerosols, snow cover, algorithms and instrumentation problems. A multisensor fusion approach is here proposed to improve the spatio-temporal continuity, consistency and accuracy of current satellite products. It is based on the use of neural networks, gap filling and temporal smoothing techniques. It is applicable to any optical sensor and satellite product. In this study, the potential of this technique was demonstrated for leaf area index (LAI) product based on MODIS and VEGETATION reflectance data. The FUSION product showed an overall good agreement with the original MODIS LAI product but exhibited a reduction of 90% of the missing LAI values with an improved monitoring of vegetation dynamics, temporal smoothness, and better agreement with ground measurements.     

Assessment of the SeaWinds scatterometer for vegetation phenology monitoring across China

Vegetation phenology tracks plants' lifecycle events, revealing the response of vegetation to global climate changes. Changes in vegetation phenology also influence fluxes of carbon, water, and energy at local and global scales. In this study, we analysed a time series of Ku-band radar backscatter measurements from the SeaWinds scatterometer on board the Quick Scatterometer (QuickSCAT) to examine canopy phenology from 2003 to 2005 across China. The thaw season SeaWinds backscatter and Moderate Resolution Imaging Spectroradiometer (MODIS) leaf area index (LAI) time series were significantly correlated in 20 of the 22 sites (p < 0.05). A weighted curve-fitting method was applied to detect the start of season and end of season from both data sets. The SeaWinds scatterometer generally detected earlier timing of spring leaf-out and later fall senescence than the MODIS LAI data sets. The SeaWinds backscatter detected phenological metrics in 75.85% of mainland China. Similar spatial patterns were observed from the SeaWinds backscatter and MODIS LAI time series; however, the average standard deviation of the scatterometer-detected metrics was lower than that of MODIS LAI products. Overall, the phenological information from the SeaWinds scatterometer could provide an alternative view on the growth dynamics of land-surface vegetation.     

Comparison of spatial sampling strategies for ground sampling and validation of MODIS LAI products

The development of an efficient ground sampling strategy is critical to assess uncertainties associated with moderate- or coarse-resolution remote-sensing products. This work presents a comparison of estimating spatial means from fine spatial resolution images using spatial random sampling (SRS), Block Kriging (BK), and Means of Surface with Nonhomogeneity (MSN) at 1 km² spatial scale. Towards this goal, we focus on the sampling strategies for ground data measurements and provide an assessment of the MODIS LAI product validated by the spatial means estimated by the above-mentioned three methods. The results of this study indicate that: (1) for its effective stratification strategies and its criteria of minimum mean square estimation error, MSN demonstrates the lowest mean squared estimation error for estimating the means of stratified nonhomogeneous surface; (2) BK is efficient in estimating the means of homogeneous surfaces without bias and with minimum mean squared estimation errors. The MODIS LAI product is assessed using the means estimated by SRS, BK, and MSN based on Landsat 8 OLI and SPOT HRV fine-resolution LAI maps. For heterogeneous surfaces, MSN results in low RMSE and high accuracy of MODIS LAI product compared with BK and SRS, whereas for homogeneous surfaces, the statistical parameters outputted by these three methods are similar. These results reveal that MSN is an effective method for estimating the spatial means for heterogeneous surfaces. There are differences in the accuracies of MODIS LAI product assessed by these three methods.     

Generating vegetation leaf area index Earth system data record from multiple sensors. Part 2: Implementation,analysis and validation

The evaluation of a new global monthly leaf area index (LAI) data set for the period July 1981 to December 2006 derived from AVHRR Normalized Difference Vegetation Index (NDVI) data is described. The physically based algorithm is detailed in the first of the two part series. Here, the implementation, production and evaluation of the data set are described. The data set is evaluated both by direct comparisons to ground data and indirectly through inter-comparisons with similar data sets. This indirect validation showed satisfactory agreement with existing LAI products, importantly MODIS, at a range of spatial scales, and significant correlations with key climate variables in areas where temperature and precipitation limit plant growth. The data set successfully reproduced well-documented spatio-temporal trends and inter-annual variations in vegetation activity in the northern latitudes and semi-arid tropics. Comparison with plot scale field measurements over homogeneous vegetation patches indicated a 7% underestimation when all major vegetation types are taken into account. The error in mean values obtained from distributions of AVHRR LAI and high-resolution field LAI maps for different biomes is within 0.5 LAI for six out of the ten selected sites. These validation exercises though limited by the amount of field data, and thus less than comprehensive, indicated satisfactory agreement between the LAI product and field measurements. Overall, the inter-comparison with short-term LAI data sets, evaluation of long term trends with known variations in climate variables, and validation with field measurements together build confidence in the utility of this new 26 year LAI record for long term vegetation monitoring and modeling studies.     

Cross-scalar satellite phenology from ground, Landsat, and MODIS data

Phenological records constructed from global mapping satellite platforms (e.g. AVHRR and MODIS) hold the potential to be valuable tools for monitoring vegetation response to global climate change. However, most satellite phenology products are not validated, and field checking coarse scale (≥ 500 m) data with confidence is a difficult endeavor. In this research, we compare phenology from Landsat (field scale, 30 m) to MODIS (500 m), and compare datasets derived from each instrument. Landsat and MODIS yield similar estimates of the start of greenness (r² = 0.60), although we find that a high degree of spatial phenological variability within coarser-scale MODIS pixels may be the cause of the remaining uncertainty. In addition, spatial variability is smoothed in MODIS, a potential source of error when comparing in situ or climate data to satellite phenology. We show that our method for deriving phenology from satellite data generates spatially coherent interannual phenology departures in MODIS data. We test these estimates from 2000 to 2005 against long-term records from Harvard Forest (Massachusetts) and Hubbard Brook (New Hampshire) Experimental Forests. MODIS successfully predicts 86% of the variance at Harvard forest and 70% of the variance at Hubbard Brook; the more extreme topography of the later is inferred to be a significant source of error. In both analyses, the satellite estimate is significantly dampened from the ground-based observations, suggesting systematic error (slopes of 0.56 and 0.63, respectively). The satellite data effectively estimates interannual phenology at two relatively simple deciduous forest sites and is internally consistent, even with changing spatial scale. We propose that continued analyses of interannual phenology will be an effective tool for monitoring native forest responses to global-scale climate variability.     

森林叶面积指数遥感反演与空间尺度转换研究

以贵州省黎平县为研究区,着重研究森林叶面积指数(LAI)的ETM遥感信息反演和向1km空间尺度转换算法.通过LAI-2000的针叶林和阔叶林等植被类型的LAI实地观测,建立实测LAI与ETM影像归一化植被指数(NDVI)的相关关系并进行LAI遥感制图,并在陆地覆盖类型遥感分类信息提取的基础上,发展了针叶林、混交林和空旷地三种地表类型LAI的向上空间尺度转换算法,以对粗分辨MODIS遥感数据的LAI产品实现LAI算法的转换与校正,并通过示例应用显示了本研究空间尺度转换算法的有效性.     

Improvements of the MODIS terrestrial gross and net primary production global data set

MODIS primary production products (MOD17) are the first regular, near-real-time data sets for repeated monitoring of vegetation primary production on vegetated land at 1-km resolution at an 8-day interval. But both the inconsistent spatial resolution between the gridded meteorological data and MODIS pixels, and the cloud-contaminated MODIS FPAR/LAI (MOD15A2) retrievals can introduce considerable errors to Collection4 primary production (denoted as C4 MOD17) results. Here, we aim to rectify these problems through reprocessing key inputs to MODIS primary vegetation productivity algorithm, resulting in improved Collection5 MOD17 (here denoted as C5 MOD17) estimates. This was accomplished by spatial interpolation of the coarse resolution meteorological data input and with temporal filling of cloud-contaminated MOD15A2 data. Furthermore, we modified the Biome Parameter Look-Up Table (BPLUT) based on recent synthesized NPP data and some observed GPP derived from some flux tower measurements to keep up with the improvements in upstream inputs. Because MOD17 is one of the down-stream MODIS land products, the performance of the algorithm can be largely influenced by the uncertainties from upstream inputs, such as land cover, FPAR/LAI, the meteorological data, and algorithm itself. MODIS GPP fits well with GPP derived from 12 flux towers over North America. Globally, the 3-year MOD17 NPP is comparable to the Ecosystem Model-Data Intercomparison (EMDI) NPP data set, and global total MODIS GPP and NPP are inversely related to the observed atmospheric CO₂ growth rates, and MEI index, indicating MOD17 are reliable products. From 2001 to 2003, mean global total GPP and NPP estimated by MODIS are 109.29 Pg C/year and 56.02 Pg C/year, respectively. Based on this research, the improved global MODIS primary production data set is now ready for monitoring ecological conditions, natural resources and environmental changes.     

Assessment of the MODIS LAI product for Australian ecosystems

The leaf area index (LAI) product from the Moderate Resolution Imaging Spectroradiometer (MODIS) is important for monitoring and modelling global change and terrestrial dynamics at many scales. The algorithm relies on spectral reflectances and a six biome land cover classification. Evaluation of the specific behaviour and performance of the product for regions of the globe such as Australia are needed to assist with product refinement and validation. We made an assessment of Collection 4 of the MODIS LAI product using four approaches: (a) assessment against a continental scale Structural Classification of Australian Vegetation (SCAV); (b) assessment against a continental scale land use classification (LUC); (c) assessment against historical field-based measurement of LAI collected prior to the Terra Mission; and (d) direct comparison of MODIS LAI with coincident field measurements of LAI, mostly from hemispherical photography. The MODIS LAI product produced a wide variety of geographically and structurally specific temporal response profiles between different classes and even for sub-groups within classes of the SCAV. Historical and concurrent field measurements indicated that MODIS LAI was giving reasonable estimates for LAI for most cover types and land use types, but that major overestimation of LAI occurs in some eastern Australian open forests and woodlands. The six biome structural land cover classification showed some significant deviations in class allocation compared to the SCAV particularly where grasslands are allocated to shrubland, savanna woodlands are allocated to shrubland, savanna and broadleaf forest, and open forests are allocated to savanna and broadleaf forest. The land cover and LAI products could benefit from some additional examination of Australian data addressing the structural representation of Eucalypt canopies in the “space of canopy realisation” for savanna and broadleaf forest classes.     

An assessment of the MODIS collection 5 leaf area index product for a region of mixed coniferous forest

Canopy leaf area index (LAI), defined as the single-sided leaf area per unit ground area, is a quantitative measure of canopy foliar area. LAI is a controlling biophysical property of vegetation function, and quantifying LAI is thus vital for understanding energy, carbon and water fluxes between the land surface and the atmosphere. LAI is routinely available from Earth Observation (EO) instruments such as MODIS. However EO-derived estimates of LAI require validation before they are utilised by the ecosystem modelling community. Previous validation work on the MODIS collection 4 (c4) product suggested considerable error especially in forested biomes, and as a result significant modification of the MODIS LAI algorithm has been made for the most recent collection 5 (c5). As a result of these changes the current MODIS LAI product has not been widely validated. We present a validation of the MODIS c5 LAI product over a 121 km² area of mixed coniferous forest in Oregon, USA, based on detailed ground measurements which we have upscaled using high resolution EO data. Our analysis suggests that c5 shows a much more realistic temporal LAI dynamic over c4 values for the site we examined. We find improved spatial consistency between the MODIS c5 LAI product and upscaled in situ measurements. However results also suggest that the c5 LAI product underestimates the upper range of upscaled in situ LAI measurements.     

Spectral mixture analysis for bi-sensor wetland mapping using Landsat TM and Terra MODIS data

Spatial and temporal resolution is essential for understanding the spatial and temporal characteristics and dynamics of wetland ecosystems. However, single satellite imagery with both high spatial resolution and high temporal frequency is currently unavailable. Instead, the development of a bi-sensor monitoring technique utilizing spatial details of middle-to-high resolution data and temporal details of coarse spatial resolution data is highly desirable. For the initial work on our time-series bi-sensor wetland mapping, the applicability of multiple endmember spectral mixture analysis (MESMA) using single-date bi-sensor imagery with different orbiting periods was investigated. Landsat-5 Thematic Mapper (TM) and Terra Moderate Resolution Image Spectrometer (MODIS) data were utilized in the Poyang Lake area in China and the Great Salt Lake area in the USA to examine three decisive elements in utilizing MESMA: (1) the method of optimal endmember selection; (2) the threshold between two- and three-endmember models; and (3) the treatment of shade fractions. As a result, we found that (1) the number of spectra for an endmember spectrum similar to other endmember spectra meeting the modelling restrictions of maximum and minimum land-cover fractions and root mean square error (RMSE) within a class (In_CoB), the number of spectra for an endmember spectrum similar to other endmember spectra meeting the modelling restrictions outside of a class (Out_CoB), the ratio of In_CoB to Out_CoB multiplied by the inverse number of spectra within the class (CoBI) and the endmember average RMSE (EAR) were optimal endmember selection methods for the TM maps, whereas CoBI, EAR and minimum average spectral angle (MASA) were optimal endmember selection methods for the MODIS maps; (2) the MODIS maps were more sensitive to change in the two- and three-endmember modelling thresholds than the TM maps; and (3) the addition of shade fractions to dark water fractions were an appropriate shade treatment. This research demonstrated how MESMA can be applied for multi-scale mapping of wetland ecosystems, how the difference in observation dates between the TM and MODIS data affects the agreement in land-cover fractions and how spectral similarity between dark water and shade affects the agreement in land-cover fractions.     

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

ABSTRACT

The 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 re?ectance 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.     

Prototyping of Land-SAF leaf area index algorithm with VEGETATION and MODIS data over Europe

The Satellite Application Facility on Land Surface Analysis (Land-SAF) aims to provide land surface variables for the meteorological and environmental science communities from EUMETSAT satellites. This study assesses the performance of a simplified (i.e. random distribution of vegetation is assumed) version of the Land-SAF algorithm for the estimation of Leaf Area Index (LAI) when prototyped with VEGETATION (processed in CYCLOPES program) and MODIS reflectances. The prototype estimates of LAI are evaluated both by comparison with validated CYCLOPES and MODIS LAI products derived from the same sensors and directly through comparison with ground-based estimates. Emphasis is given on evaluating the impact of the algorithm and input data on LAI retrieval discrepancies. Analysis is achieved over Europe for the 2000-2003 period. The results demonstrate the capacity of the Land-SAF algorithm to retrieve consistent LAI estimates from multiple optical sensors even when their reflectances present systematic differences. High spatial and temporal consistencies between Land-SAF prototype estimates and existing LAI products are found. The differences between Land-SAF and CYCLOPES LAI are lower than their uncertainties (RMSE (relative RMSE) within 0.4 (30%)). Land-SAF prototype estimates and MODIS LAI show larger discrepancies mainly due to differences in the vegetation structure representation and algorithm assumptions (RMSE ranging from 0.2 (30%) up to 0.8 (40%)). Land-SAF prototype provides higher LAI values than MODIS for herbaceous canopies (i.e. shrubs, grasses and crops) and lower values for woody biomes (i.e. savannas and forests). Direct validation indicates that LAI estimates from prototyping of the Land-SAF algorithm with CYCLOPES and MODIS reflectances achieve similar performances (differences with ground measurements are lower than 0.5 LAI units in 60% and 50% of the cases, respectively) as CYCLOPES and MODIS LAI products. Results from this prototyping exercise appear useful for improved retrieval of LAI and constitute a step forward for refinement, validation and consolidation of the Land-SAF algorithm.