GLASS叶面积指数产品验证

在国家高技术研究发展计划(863计划)重点项目的支持下,已利用MODIS和AVHRR地表反射率数据生成了1981年—2012年的GLASS(Global LAnd Surface Satellite)叶面积指数(LAI)产品。本文从两个方面对GLASS LAI产品的质量进行分析和评价:(1)与现有的全球LAI产品进行比较,分析GLASS LAI产品的时空变化特征;(2)利用LAI的地面测量数据,对GLASS LAI的精度进行评价。研究结果表明:GLASS LAI与CCRS LAI在高纬度和赤道附近区域的差异较大;相对而言,GLASS LAI与MODIS(主算法反演)和CYCLOPES LAI在空间分布上具有更好的一致性;GLASS和CYCLOPES LAI的时间序列曲线连续平滑,MODIS LAI在一些区域的植被生长季节存在剧烈的跳跃;与LAI的地面测量数据进行比较,GLASS LAI产品的R2为0.76,RMSE为0.51,结果明显优于MODIS和C YCLOPES LAI产品。     

中国区域MuSyQ叶面积指数产品验证与分析

对多源遥感数据协同生产的2010年—2015年中国区域1 km空间分辨率5天合成的MuSyQ（Multi-source data Synergized Quantitative remote sensing production system）叶面积指数LAI产品进行验证。参考现有的LAI产品（MODIS c5,GLASS LAI）和中国生态系统研究网络部分农田和森林站点可用的LAI地面测量数据,从时空连续性、时空一致性、精度和准确性等方面对中国区域的MuSyQ LAI产品进行定性和定量分析与评价。结果表明：（1） MuSyQ LAI产品在保证精度优于MODIS产品的情况下,时间分辨率和时空连续性均有提高。MuSyQ LAI与其他LAI产品（MODIS c5,GLASS LAI）在整体上有很好的一致性（RMSE=1.0,RMSE=0.81）,但对常绿阔叶林高值处的描述不稳定;（2） 与LAI地面测量数据相比,MuSyQ LAI产品与地面参考图对比结果较好（最高相关性（R²=0.54）和较低总体误差（RMSE=0.96））,其在阔叶作物生长季高值处有些许低估且在某些阔叶林站点有些高估。整体上,MuSyQ LAI产品呈现出较高的精度,可靠的空间分布和连续稳定的时间分布,且对森林LAI的描述具有更可靠的动态范围。     

Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data

A sufficient number of satellite acquisitions in a growing season are essential for deriving agronomic indicators, such as green leaf area index (GLAI), to be assimilated into crop models for crop productivity estimation. However, for most high resolution orbital optical satellites, it is often difficult to obtain images frequently due to their long revisit cycles and unfavorable weather conditions. Data fusion algorithms, such as the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), have been developed to generate synthetic data with high spatial and temporal resolution to address this issue. In this study, we evaluated the approach of assimilating GLAI into the Simple Algorithm for Yield Estimation model (SAFY) for winter wheat biomass estimation. GLAI was estimated using the two-band Enhanced Vegetation Index (EVI2) derived from data acquired by the Operational Land Imager (OLI) onboard the Landsat-8 and a fusion dataset generated by blending the Moderate-Resolution Imaging Spectroradiometer (MODIS) data and the OLI data using the STARFM and ESTARFM models. The fusion dataset had the temporal resolution of the MODIS data and the spatial resolution of the OLI data. Key parameters of the SAFY model were optimised through assimilation of the estimated GLAI into the crop model using the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm. A good agreement was achieved between the estimated and field measured biomass by assimilating the GLAI derived from the OLI data (GLAI_L) alone (R² = 0.77 and RMSE = 231 g m⁻²). Assimilation of GLAI derived from the fusion dataset (GLAI_F) resulted in a R² of 0.71 and RMSE of 193 g m⁻² while assimilating the combination of GLAI_L and GLAI_F led to further improvements (R² = 0.76 and RMSE = 176 g m⁻²). Our results demonstrated the potential of using the fusion algorithms to improve crop growth monitoring and crop productivity estimation when the number of high resolution remote sensing data acquisitions is limited.     

A physically based approach to model LAI from MODIS 250 m data in a tropical region

A time series of leaf area index (LAI) has been developed based on 16-day normalized difference vegetation index (NDVI) data from the Moderate Resolution Imaging Spectroradiometer (MODIS) at 250 m resolution (MOD250_LAI). The MOD250_LAI product uses a physical radiative transfer model which establishes a relationship between LAI, fraction of vegetation cover (FVC) and given patterns of surface reflectance, view-illumination conditions and optical properties of vegetation. In situ measurements of LAI and FVC made at 166 plots using hemispherical photography served for calibration of model parameters and validation of modelling results. Optical properties of vegetation cover, summarized by the light extinction coefficient, were computed at the local (pixel) level based on empirical models between ground-measured tree crown architecture at 85 sampling plots and spectral values in Landsat ETM+ bands. Influence of view-illumination conditions on optical properties of canopy was simulated by a view angle geometry model incorporating the solar zenith angle and the sensor viewing angle. The results revealed high compatibility of the produced MOD250_LAI data set with ground truth information and the 30 m resolution Landsat ETM+ LAI estimated using the similar algorithm. The produced MOD250_LAI was also compared with the global MODIS 1000-m LAI product (MOD15A2 LAI). Results show good consistency of the spatial distribution and temporal dynamics between the two LAI products. However, the results also showed that the annual LAI amplitude by the MOD15A2 product is significantly higher than by the MOD250_LAI. This higher amplitude is caused by a considerable underestimation of the tropical rainforest LAI by the MOD15A2 during the seasonal phases of low leaf production.     

Directional effects on NDVI and LAI retrievals from MODIS: A case study in Brazil with soybean

The Moderate Resolution Imaging Spectroradiometer (MODIS) is largely used to estimate Leaf Area Index (LAI) using radiative transfer modeling (the “main” algorithm). When this algorithm fails for a pixel, which frequently occurs over Brazilian soybean areas, an empirical model (the “backup” algorithm) based on the relationship between the Normalized Difference Vegetation Index (NDVI) and LAI is utilized. The objective of this study is to evaluate directional effects on NDVI and subsequent LAI estimates using global (biome 3) and local empirical models, as a function of the soybean development in two growing seasons (2004–2005 and 2005–2006). The local model was derived from the pixels that had LAI values retrieved from the main algorithm. In order to keep the reproductive stage for a given cultivar as a constant factor while varying the viewing geometry, pairs of MODIS images acquired in close dates from opposite directions (backscattering and forward scattering) were selected. Linear regression relationships between the NDVI values calculated from these two directions were evaluated for different view angles (0–25°; 25–45°; 45–60°) and development stages (<45; 45–90; >90 days after planting). Impacts on LAI retrievals were analyzed. Results showed higher reflectance values in backscattering direction due to the predominance of sunlit soybean canopy components towards the sensor and higher NDVI values in forward scattering direction due to stronger shadow effects in the red waveband. NDVI differences between the two directions were statistically significant for view angles larger than 25°. The main algorithm for LAI estimation failed in the two growing seasons with gradual crop development. As a result, up to 94% of the pixels had LAI values calculated from the backup algorithm at the peak of canopy closure. Most of the pixels selected to compose the 8-day MODIS LAI product came from the forward scattering view because it displayed larger LAI values than the backscattering. Directional effects on the subsequent LAI retrievals were stronger at the peak of the soybean development (NDVI values between 0.70 and 0.85). When the global empirical model was used, LAI differences up to 3.2 for consecutive days and opposite viewing directions were observed. Such differences were reduced to values up to 1.5 with the local model. Because of the predominance of LAI retrievals from the MODIS backup algorithm during the Brazilian soybean development, care is necessary if one considers using these data in agronomic growing/yield models.     

Retrieval of effective leaf area index (LAIe) and leaf area density (LAD) profile at individual tree level using high density multi-return airborne LiDAR

As an important canopy structure indicator, leaf area index (LAI) proved to be of considerable implications for forest ecosystem and ecological studies, and efficient techniques for accurate LAI acquisitions have long been highlighted. Airborne light detection and ranging (LiDAR), often termed as airborne laser scanning (ALS), once was extensively investigated for this task but showed limited performance due to its low sampling density. Now, ALS systems exhibit more competing capacities such as high density and multi-return sampling, and hence, people began to ask the questions like—“can ALS now work better on the task of LAI prediction?” As a re-examination, this study investigated the feasibility of LAI retrievals at the individual tree level based on high density and multi-return ALS, by directly considering the vertical distributions of laser points lying within each tree crown instead of by proposing feature variables such as quantiles involving laser point distribution modes at the plot level. The examination was operated in the case of four tree species (i.e. Picea abies, Pinus sylvestris, Populus tremula and Quercus robur) in a mixed forest, with their LAI-related reference data collected by using static terrestrial laser scanning (TLS). In light of the differences between ALS- and TLS-based LAI characterizations, the methods of voxelization of 3D scattered laser points, effective LAI (LAIe) that does not distinguish branches from canopies and unified cumulative LAI (ucLAI) that is often used to characterize the vertical profiles of crown leaf area densities (LADs) was used; then, the relationships between the ALS- and TLS-derived LAIes were determined, and so did ucLAIs. Tests indicated that the tree-level LAIes for the four tree species can be estimated based on the used airborne LiDAR (R² = 0.07, 0.26, 0.43 and 0.21, respectively) and their ucLAIs can also be derived. Overall, this study has validated the usage of the contemporary high density multi-return airborne LiDARs for LAIe and LAD profile retrievals at the individual tree level, and the contribution are of high potential for advancing forest ecosystem modeling and ecological understanding.     

基于HJ-1B影像的内蒙古草地叶面积指数反演

针对HJ-1B星CCD数据,以内蒙古草地为研究对象,应用植被指数模型、PROSAIL模型以及热暗点指数模型对研究区草地进行叶面积指数定量估算研究,同时利用LAI-2000获取的实测数据进行反演结果检验。结果显示:复合型植被指数统计模型的R2=0.7,物理模型为0.67;而传统植被指数模型R2仅为0.61。实验表明热暗点植被指数模型比传统植被指数模型更适于进行内蒙古草地LAI定量反演。     

Estimation of leaf area index using PROSAIL based LUT inversion,MLRA-GPR and empirical models: Case study of tropical deciduous forest plantation,North India

Forests play a vital role in biological cycles and environmental regulation. To understand the key processes of forest canopies (e.g., photosynthesis, respiration and transpiration), reliable and accurate information on spatial variability of Leaf Area Index (LAI), and its seasonal dynamics is essential. In the present study, we assessed the performance of biophysical parameter (LAI) retrieval methods viz. Look-Up Table (LUT)-inversion, MLRA-GPR (Machine Learning Regression Algorithm- Gaussian Processes Regression) and empirical models, for estimating the LAI of tropical deciduous plantation using ARTMO (Automated Radiative Transfer Models Operator) tool and Sentinel-2 satellite images. The study was conducted in Central Tarai Forest Division, Haldwani, located in the Uttarakhand state, India. A total of 49 ESUs (Elementary Sampling Unit) of 30 m × 30 m size were established based on variability in composition and age of plantation stands. In-situ LAI was recorded using plant canopy imager during the leaf growing, peak and senescence seasons. The PROSAIL model was calibrated with site-specific biophysical and biochemical parameters before used to the predicted LAI. The plantation LAI was also predicted by an empirical approach using optimally chosen Sentinel-2 vegetation indices. In addition, Sentinel-2 and MODIS LAI products were evaluated with respect to LAI measurements. MLRA-GPR offered best results for predicting LAI of leaf growing (R² = 0.9, RMSE = 0.14), peak (R² = 0.87, RMSE = 0.21) and senescence (R² = 0.86, RMSE = 0.31) seasons while LUT inverted model outperformed VI’s based parametric regression model. Vegetation indices (VIs) derived from 740 nm, 783 nm and 2190 nm band combinations of Sentinel-2 offered the best prediction of LAI.     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.     

Comparison of UAV and WorldView-2 imagery for mapping leaf area index of mangrove forest

Unmanned Aerial Vehicle (UAV) remote sensing has opened the door to new sources of data to effectively characterize vegetation metrics at very high spatial resolution and at flexible revisit frequencies. Successful estimation of the leaf area index (LAI) in precision agriculture with a UAV image has been reported in several studies. However, in most forests, the challenges associated with the interference from a complex background and a variety of vegetation species have hindered research using UAV images. To the best of our knowledge, very few studies have mapped the forest LAI with a UAV image. In addition, the drawbacks and advantages of estimating the forest LAI with UAV and satellite images at high spatial resolution remain a knowledge gap in existing literature. Therefore, this paper aims to map LAI in a mangrove forest with a complex background and a variety of vegetation species using a UAV image and compare it with a WorldView-2 image (WV2).In this study, three representative NDVIs, average NDVI (AvNDVI), vegetated specific NDVI (VsNDVI), and scaled NDVI (ScNDVI), were acquired with UAV and WV2 to predict the plot level (10 × 10 m) LAI. The results showed that AvNDVI achieved the highest accuracy for WV2 (R² = 0.778, RMSE = 0.424), whereas ScNDVI obtained the optimal accuracy for UAV (R² = 0.817, RMSE = 0.423). In addition, an overall comparison results of the WV2 and UAV derived LAIs indicated that UAV obtained a better accuracy than WV2 in the plots that were covered with homogeneous mangrove species or in the low LAI plots, which was because UAV can effectively eliminate the influence from the background and the vegetation species owing to its high spatial resolution. However, WV2 obtained a slightly higher accuracy than UAV in the plots covered with a variety of mangrove species, which was because the UAV sensor provides a negative spectral response function(SRF) than WV2 in terms of the mangrove LAI estimation.     

Mapping paddy rice planting area in rice-wetland coexistent areas through analysis of Landsat 8 OLI and MODIS images

Accurate and up-to-date information on the spatial distribution of paddy rice fields is necessary for the studies of trace gas emissions, water source management, and food security. The phenology-based paddy rice mapping algorithm, which identifies the unique flooding stage of paddy rice, has been widely used. However, identification and mapping of paddy rice in rice-wetland coexistent areas is still a challenging task. In this study, we found that the flooding/transplanting periods of paddy rice and natural wetlands were different. The natural wetlands flood earlier and have a shorter duration than paddy rice in the Panjin Plain, a temperate region in China. We used this asynchronous flooding stage to extract the paddy rice planting area from the rice-wetland coexistent area. MODIS Land Surface Temperature (LST) data was used to derive the temperature-defined plant growing season. Landsat 8 OLI imagery was used to detect the flooding signal and then paddy rice was extracted using the difference in flooding stages between paddy rice and natural wetlands. The resultant paddy rice map was evaluated with in-situ ground-truth data and Google Earth images. The estimated overall accuracy and Kappa coefficient were 95% and 0.90, respectively. The spatial pattern of OLI-derived paddy rice map agrees well with the paddy rice layer from the National Land Cover Dataset from 2010 (NLCD-2010). The differences between Rice_Landsat and Rice_NLCD are in the range of ±20% for most 1-km grid cell. The results of this study demonstrate the potential of the phenology-based paddy rice mapping algorithm, via integrating MODIS and Landsat 8 OLI images, to map paddy rice fields in complex landscapes of paddy rice and natural wetland in the temperate region.     

植被叶面积指数估算的红边指数建模分析

针对在路域环境监测中,如何精确估算叶面积指数问题,该文提出以长韶娄高速路域为研究区,筛选出4种常用植被指数和4种红边指数两类指数,分别构建了经验模型和机器学习的反演模型,利用Sentinel-2影像数据和同步的LAI-2000地面实测数据完成路域植被叶面积指数反演。结果表明,红边波段参与运算的植被指数与植被叶面积指数敏感性是显著相关,红边指数在反演精度上更优。由此可知,相较于常见植被指数,红边指数增强了其与叶面积指数的敏感性,提高了叶面积指数估算模型精度。     

Assessment of leaf carotenoids content with a new carotenoid index: Development and validation on experimental and model data

Leaf carotenoids content (LCar) is an important indicator of plant physiological status. Accurate estimation of LCar provides valuable insight into early detection of stress in vegetation. With spectroscopy techniques, a semi-empirical approach based on spectral indices was extensively used for carotenoids content estimation. However, established spectral indices for carotenoids that generally rely on limited measured data, might lack predictive accuracy for carotenoids estimation in various species and at different growth stages. In this study, we propose a new carotenoid index (CARI) for LCar assessment based on a large synthetic dataset simulated from the leaf radiative transfer model PROSPECT-5, and evaluate its capability with both simulated data from PROSPECT-5 and 4SAIL and extensive experimental datasets: the ANGERS dataset and experimental data acquired in field experiments in China in 2004. Results show that CARI was the index most linearly correlated with carotenoids content at the leaf level using a synthetic dataset (R² = 0.943, RMSE = 1.196 μg/cm²), compared with published spectral indices. Cross-validation results with CARI using ANGERS data achieved quite an accurate estimation (R² = 0.545, RMSE = 3.413 μg/cm²), though the RBRI performed as the best index (R² = 0.727, RMSE = 2.640 μg/cm²). CARI also showed good accuracy (R² = 0.639, RMSE = 1.520 μg/cm²) for LCar assessment with leaf level field survey data, though PRI performed better (R² = 0.710, RMSE = 1.369 μg/cm²). Whereas RBRI, PRI and other assessed spectral indices showed a good performance for a given dataset, overall their estimation accuracy was not consistent across all datasets used in this study. Conversely CARI was more robust showing good results in all datasets. Further assessment of LCar with simulated and measured canopy reflectance data indicated that CARI might not be very sensitive to LCar changes at low leaf area index (LAI) value, and in these conditions soil moisture influenced the LCar retrieval accuracy.     

Mapping forest leaf area index using reflectance and textural information derived from WorldView-2 imagery in a mixed natural forest area in Florida,US

The leaf area index (LAI) of plant canopies is an important structural parameter that controls energy, water, and gas exchanges of plant ecosystems. Remote sensing techniques may offer an alternative for measuring and mapping forest LAI at a landscape scale. Given the characteristics of high spatial/spectral resolution of the WorldView-2 (WV2) sensor, it is of significance that the textural information extracted from WV2 multispectral (MS) bands will be first time used in estimating and mapping forest LAI. In this study, LAI mapping accuracies would be compared from (a) spatial resolutions between 2-m WV2 MS data and 30-m Landsat TM imagery, (b) the nature of variables between spectrum-based features and texture-based features, and (c) sensors between TM and WV2. Therefore spectral/textural features (SFs) were first selected and tested; then a canonical correlation analysis was performed with different data sets of SFs and LAI measurement; and finally linear regression models were used to predict and map forest LAI with canonical variables calculated from image data. The experimental results demonstrate that for estimating and mapping forest LAI, (i) using high resolution data (WV2) is better than using relatively low resolution data (TM); (ii) extracted from the same WV2 data, texture-based features have higher capability than that of spectrum-based features; (iii) a combination of spectrum-based features with texture-based features could lead to even higher accuracy of mapping forest LAI than their either one separately; and (iv) WV2 sensor outperforms TM sensor significantly. However, we need to address the possible overfitting phenomenon that might be brought in by using more input variables to develop models. In addition, the experimental results also indicate that the red-edge band in WV2 was the worst on estimating LAI among WV2 MS bands and the WV2 MS bands in the visible range had a much higher correlation with ground measured LAI than that red-edge and NIR bands did.     

Diurnal and seasonal impacts of urbanization on the urban thermal environment: A case study of Beijing using MODIS data

Beijing has experienced rapid urbanization and associated urban heat island effects and air pollution. In this study, a contribution index was proposed to explore the effect of urbanization on land surface temperature (LST) using Moderate-Resolution Imaging Spectroradiometer (MODIS)-derived data with high temporal resolution. The analysis indicated that different zones and landscapes make diurnally and seasonally different contributions to the regional thermal environment. The differences in contributions by the three main functional zones resulted from differences in their landscape compositions. The roles of landscapes in this process varied diurnally and seasonally. Urban land was the most important contributor to increases in regional LSTs. The contributions of cropland and forest varied distinctly between daytime and nighttime owing to differences in their thermal inertias. Vegetation had a notable cooling effect as the normalized vegetation difference index (NDVI) increased during summer. However, when the NDVI reached a certain value, the nighttime LST shifted markedly in other seasons. The results suggest that urban design based on vegetation partitions would be effective for regulating the thermal environment.     

Retrieval of leaf area index in alpine wetlands using a two-layer canopy reflectance model

In this paper, we focused on the retrieval of the LAI in an alpine wetland located in western part of China in late August and early July 2011. A two-layer canopy reflectance model (ACRM) was used to establish the relationships between the LAI and the reflectance of near-infrared (NIR) and red (RED) wavebands. The reflectance data were derived from Landsat TM L1T product and the Terra and Aqua MODIS 16-day and 8-day composite reflectance products (MOD/MYD09) at 250 m resolution. Due to the lack of the information about some major input parameters for ACRM, which are sensitive to model outputs in the reflectance of NIR and RED wavebands, the inverse problem was ill-posed. To overcome this problem, a method of increasing the sensitivity of the LAI while reducing the influence of other model free parameters based on the study of free parameters’ sensitivity to the ACRM outputs and the region’s features was studied. The area of interest was divided into two parts using the approximately statistic normalized difference vegetation index (NDVI) value around 0.5. One part was sparse vegetation (0.1 < NDVI < 0.5), which is more sensitive to soil background effects and less sensitive to the canopy biophysical and biochemical variables. The other part was dense vegetation (0.5 ≤ NDVI < 1.0), which is less sensitive to soil background effects and more sensitive to plant canopies and leaf parameters. Then, the relationships of ρ_nir–LAI and ρ_red–LAI were established using a look-up table algorithm for the two parts. Furthermore, a regularization technique for fast pixel-wise retrieval was introduced to reduce the elements of LUT sets while maintaining a relatively high accuracy. The results were very promising compared to the field measured LAI values that the correlation (R²) of the measured LAI values and retrieved LAI values reached 0.95, and the root-mean-square deviation (RMSD) was 0.33 for late August, 2011, while the R² reached 0.82 and RMSD was 0.25 for early July 2011.     

地面站点叶面积指数观测的空间代表性评价——以CERN站网观测为例

在叶面积指数LAI(Leaf Area Index)产品真实性检验中,地面站点的多时相连续观测LAI数据是重要的验证数据来源。当站点观测范围与产品像元尺度不一致时,站点观测LAI直接用于产品验证可能为验证结果带来误差。因此,在验证之前需要分析站点观测对像元尺度的空间代表性,选择空间代表性好的观测来验证产品,从而减小尺度效应带来的验证误差。以往的研究只是简单的定性说明研究区域,并直接用站点测量数据对产品进行验证,缺少一套系统的站点观测在产品像元尺度内空间代表性评价的方法体系。本文提出了站点LAI观测的空间代表性评价方法,建立了评价指标DVTP(Dominant Vegetation Type Percent)、RSSE(Relative Spatial Sampling Error)和CS(Coefficient of Sill),构建了空间代表性评价分级体系。以中国生态系统研究网络CERN(Chinese Ecosystem Research Network)农田站和森林站LAI观测为例,对站点观测在1 km产品像元尺度内空间代表性进行评价,并分析评价前后站点观测对MODIS LAI产品验证精度的影响。结果显示,本文提出的方法能够有效地对不同站点LAI观测在产品像元尺度内空间代表性进行质量分级,且年际间的站点观测空间代表性较为一致。评价方法能够去掉在特定产品像元尺度下空间代表性不好的观测数据,一定程度上提高验证数据集对产品验证精度的可靠性。     

Estimation of LAI and above-ground biomass in deciduous forests: Western Ghats of Karnataka, India

This study demonstrates the potentials of IRS P6 LISS-IV high-resolution multispectral sensor (IGFOV  6 m)-based estimation of biomass in the deciduous forests in the Western Ghats of Karnataka, India. Regression equations describing the relationship between IRS P6 LISS-IV data-based vegetation index (NDVI) and field measured leaf area index (ELAI) and estimated above-ground biomass (EAGB) were derived. Remote sensing (RS) data-based leaf area index (PLAI) image is generated using regression equation based on NDVI and ELAI (r² = 0.68, p ≤ 0.05). RS-based above-ground biomass (PAGB) image was generated based on regression equation developed between PLAI and EAGB (r² = 0.63, p ≤ 0.05). The mean value of estimated above-ground biomass and RS-based above-ground biomass in the study area are 280(±72.5) and 297.6(±55.2) Mg ha⁻¹, respectively. The regression models generated in the study between NDVI and LAI; LAI and biomass can also help in generating spatial biomass map using RS data alone. LISS-IV-based estimation of biophysical parameters can also be used for the validation of various coarse resolution satellite products derived from the ground-based measurements alone.     

Satellite-based automated burned area detection: A performance assessment of the MODIS MCD45A1 in the Brazilian savanna

Burnings, which cause major changes to the environment, can be effectively monitored via satellite data, regarding both the identification of active fires and the estimation of burned areas. Among the many orbital sensors suitable for mapping burned areas on global and regional scales, the moderate resolution imaging spectroradiometer (MODIS), on board the Terra and Aqua platforms, has been the most widely utilized. In this study, the performance of the MODIS MCD45A1 burned area product was thoroughly evaluated in the Brazilian savanna, the second largest biome in South America and a global biodiversity hotspot, characterized by a conspicuous climatic seasonality and the systematic occurrence of natural and anthropogenic fires. Overall, September MCD45A1 polygons (2000–2012) compared well to the Landsat-based reference mapping (r² = 0.92) and were closely accompanied, on a monthly basis, by MOD14 and MYD14 hotspots (r² = 0.89), although large omissions errors, linked to landscape patterns, structures, and overall conditions depicted in each reference image, were observed. In spite of its spatial and temporal limitations, the MCD45A1 product proved instrumental for mapping and understanding fire behavior and impacts on the Cerrado landscapes.     

Crop area mapping in West Africa using landscape stratification of MODIS time series and comparison with existing global land products

In Africa, food security early warning systems use satellite-derived data concerning crop conditions and agricultural production. Such systems can be improved if they are provided with a more reliable estimation of the cultivated area at national scale. This paper evaluates the potential of using time series from the MODerate resolution Imaging Spectroradiometer MOD13Q1 (16-day composite of normalized difference vegetation index at 250 m resolution) to extract cultivated areas in the fragmented rural landscapes of Mali. To this end, we first stratified Southern Mali into 13 rural landscapes based on the spatio-temporal variability of NDVI and textural indices, using an object-oriented classification scheme.The accuracy of the resulting map (MODIS_crop) and how it compares with existing coarse-resolution global land products (GLC2000 Africa, GLOBCOVER, MODIS V05 and ECOCLIMAP-II), was then assessed against six crop/non-crop maps derived from SPOT 2.5 m resolution images used as references. For crop areal coverage, the MODIS_crop cultivated map was successful in assessing the overall cultivated area at five out of the six validation sites (less than 6% of the absolute difference), while in terms of crop spatial distribution, the producer accuracy was between 33.1% and 80.8%. This accuracy was linearly correlated with the mean patch size index calculated on the SPOT crop maps (r² = 0.8). Using the Pareto boundary as an accuracy assessment method at the study sites, we showed that (i) 20-40% of the classification crop error was due to the spatial resolution of the MODIS sensor (250 m), and that (ii) compared to MODIS V05, which otherwise performed better than the other existing products, MODIS_crop generally minimized omission-commission errors. A spatial validation of the different products was carried out using SPOT image classifications as reference. In the corresponding error matrices, the fraction of correctly classified pixels for our product was 70%, compared to 58% for MODIS V05, while it ranged between 40% and 51% for the GLC2000, the ECOCLIMAP-II and the GLOBCOVER.