Sentinel-2数据的冬小麦地上干生物量估算及评价

作物生物量快速精确的监测对于农业资源的合理利用与农田的精准管理具有重要意义。近年来,遥感技术因其独特的优势已被广泛用于作物生物量的估算中。本文主要针对不同宽波段植被指数在冬小麦生物量(文中的生物量均是指地上干生物量)估算方面的表现进行探索。首先利用欧洲空间局最新的Sentinel-2A卫星数据提取出17种常见的植被指数,之后分别构建其与相应时期内采集的冬小麦地上生物量间的最优估算模型,通过分析两者间的相关性与敏感性,获取适宜进行生物量估算的指数。最后,绘制了研究区的生物量空间分布图。结果表明,所选的植被指数均与生物量显著相关。其中,红边叶绿素指数(CI_(re))与生物量的估算精度最高(决定性系数R~2为0.83;均方根误差RMSE为180.29 g·m~(–2))。虽然相关性较高,但部分指数,如归一化差值植被指数(NDVI)等在生物量较高时会出现饱和现象,从而导致生物量的低估。而加入红边波段的指数不仅能够延缓指数的饱和趋势,而且能够提高反演精度。此外,通过敏感性分析发现,归一化差值指数和比值指数分别在作物生长的早期和中后期对生物量的变化保持较高的敏感性。由于红边比值指数(SR_(re))和MERIS叶绿素敏感指数(MTCI)在冬小麦全生长季内一直对生物量的变化保持高灵敏性,二者是生物量估算中最为稳定的指数。     

Monitoring of water stress in wheat using multispectral indices derived from Landsat-TM

Detection of crop water stress is crucial for efficient irrigation water management. Potential of Satellite data to provide spatial and temporal dynamics of crop growth conditions makes it possible to monitor crop water stress at regional level. This study was conducted in parts of western Uttar Pradesh and Haryana. Multi-temporal Landsat data were used for detecting wheat crop water stress using vegetation indices (VIs), viz. vegetation water stress index (VWSI) and land surface wetness index water stress factor (Ws_LSWI). The estimated water stress from satellite data-based VIs was validated by water stress factor (Ws) derived from flux-tower data. The study observed Ws_LSWI to be better index for water stress detection. The results indicated that Ws_LSWI was superior over other index showing RMSE = 0.12, R² = 0.65, whereas VWSI showed overestimated values with mean RD 4%.     

WorldView-2纹理的森林地上生物量反演

使用高空间分辨率卫星WorldView-2的多光谱遥感影像,构建植被指数和纹理因子等遥感因子与森林地上生物量的关系方程,并计算模型估测精度和均方根误差,探索高分辨率数据的光谱与纹理信息在温带森林地上生物量估测应用中的潜力。以黑龙江省凉水自然保护区温带天然林及天然次生林为研究对象,通过灰度共生矩阵(GLCM)、灰度差分向量(GLDV)及和差直方图(SADH)对高分辨率遥感影像进行纹理信息提取,并利用外业调查的74个样地地上生物量与遥感因子建立参数估计模型。提取的遥感因子包括6种植被指数(比值植被指数RVI、差值植被指数DVI、规一化植被指数NDVI、增强植被指数EVI、土壤调节植被指数SAVI和修正的土壤调节植被指数MSAVI)以及3类纹理因子(GLCM、GLDV和SADH)。为避免特征变量个数较多对估测模型造成过拟合,利用随机森林算法对提取的遥感因子进行特征选择,将最优的特征变量输入模型参与建模估测。采用支持向量回归(SVR)进行生物量建模及验证,结果显示选入模型的和差直方图均值(sadh_mean)、灰度共生矩阵方差(glcm_var)和差值植被指数(DVI)等遥感因子对森林地上生物量有较好的解释效果;植被指数+纹理因子组合的模型获得较精确的AGB估算结果(R2=0.85,RMSE=42.30 t/ha),单独使用植被指数的模型精度则较低(R~2=0.69,RMSE=61.13 t/ha)。     

Estimation of Crown Closure and Tree Density Using Landsat TM Satellite Images in Mixed Forest Stands

The objective of this study was to investigate the relationship between crown closure and tree density in mixed forest stands using Landsat Thematic Mapper (TM) reflectance values (TM 1- TM 5 and TM 7) and six vegetation indices (SR, DVI, SAVI, NDVI, TVI and NLI). In this study, multiple regression analysis was used to estimate the relationships between the crown closure and tree density (number of tree stems per hectare) using reflectance values and vegetation indices (VIs). The results demonstrated that the model that used SR and DVI had the best performances in terms of crown closure (R²?=?0.674) and the model that used the DVI and SAVI had the best performances in terms of tree density (R²?=?0.702). The regression model that used TM 1, TM 3 together with TM 4 showed the performances of the crown closure (R²?=?0.610) and the regression model that used TM 1 showed the performances of the tree density (0.613). Results obtained from this research show that vegetation indices (VIs) were a better predictor of crown closure and tree density than other TM bands.     

Estimating salinity stress in sugarcane fields with spaceborne hyperspectral vegetation indices

The presence of salt in the soil profile negatively affects the growth and development of vegetation. As a result, the spectral reflectance of vegetation canopies varies for different salinity levels. This research was conducted to (1) investigate the capability of satellite-based hyperspectral vegetation indices (VIs) for estimating soil salinity in agricultural fields, (2) evaluate the performance of 21 existing VIs and (3) develop new VIs based on a combination of wavelengths sensitive for multiple stresses and find the best one for estimating soil salinity. For this purpose a Hyperion image of September 2, 2010, and data on soil salinity at 108 locations in sugarcane (Saccharum officina L.) fields were used. Results show that soil salinity could well be estimated by some of these VIs. Indices related to chlorophyll absorption bands or based on a combination of chlorophyll and water absorption bands had the highest correlation with soil salinity. In contrast, indices that are only based on water absorption bands had low to medium correlations, while indices that use only visible bands did not perform well. From the investigated indices the optimized soil-adjusted vegetation index (OSAVI) had the strongest relationship (R² = 0.69) with soil salinity for the training data, but it did not perform well in the validation phase. The validation procedure showed that the new salinity and water stress indices (SWSI) implemented in this study (SWSI-1, SWSI-2, SWSI-3) and the Vogelmann red edge index yielded the best results for estimating soil salinity for independent fields with root mean square errors of 1.14, 1.15, 1.17 and 1.15 dS/m, respectively. Our results show that soil salinity could be estimated by satellite-based hyperspectral VIs, but validation of obtained models for independent data is essential for selecting the best model.     

Comprehensive study of the biophysical parameters of agricultural crops based on assessing Landsat 8 OLI and Landsat 7 ETM+ vegetation indices

For three agricultural crop types, winter wheat (Triticum aestivum L.), barley (Hordeum vulgare L.), and canola (Brassica napus L.), we estimated biophysical parameters including fresh and dry biomass, leaf area index (LAI), and vegetation water content, for which we found the equivalent water thickness (EWT), fuel moisture content per fresh weight (FMC_FW), and fuel moisture content per dry weight (FMC_DW). We performed these estimations using data from the newly launched Landsat 8 Operational Land Imager (OLI) sensor, as well as its predecessor the Landsat 7 Enhanced Thematic Mapper Plus (ETM+). Progress in the design of the new sensor (i.e., Landsat 8), including narrower near-infrared (NIR) wavebands, higher signal-to-noise ratio (SNR), and greater radiometric resolution highlights the necessity to investigate the biophysical parameters of agricultural crops, especially compared to data from its predecessor. This study aims to evaluate vegetation indices (VIs) derived from the Landsat 8 OLI and the Landsat 7 ETM+. Both the Landsat 8 OLI and Landsat 7 ETM+ VIs agreed well with in-situ data measurements. However, the Landsat 8 OLI-derived VIs were generally more consistent with in situ data than the Landsat 7 ETM+ VIs. We also note that the Landsat 8 OLI is better able to capture the small variability of the VIs because of its higher SNR and wider radiometric range; in addition, the saturation phenomenon occurred earlier for the Landsat 7 ETM+ than for the Landsat 8 OLI. This indicates that the new sensor is better able to estimate the biophysical parameters of crops.     

Delineation of management zones in agricultural fields using cover–crop biomass estimates from PlanetScope data

Several methods have been proposed to delineate management zones in agricultural fields, which can guide interventions of the farmers to increase crop yield. In this study, we propose a new approach using remote sensing data to delineate management zones at three farm sites located in southern Brazil. The approach is based on the hypothesis that the measured aboveground biomass (AGB) of the cover crops is correlated with the measured cash-crop yield and can be estimated from surface reflectance and/or vegetation indices (VIs). Therefore, we used seven different statistical models to estimate AGB of three cover crops (forage turnip, white oats, and rye) in the season prior to cash-crop planting. Surface reflectance and VIs were used as predictors to test the performance of the models. They were obtained from high spatial and temporal resolution data of the PlanetScope (PS) constellation of satellites. From the time series of 30 images acquired in 2017, we used the PS data that matched the dates of the field campaigns to build the models. The results showed that the satellite AGB estimates of the cover crops at the date of maximum VI response at the beginning of the flowering stage were useful to delineate the management zones. The cover-crop AGB models that presented the highest coefficient of determination (R²) and the lowest root mean square (RMSE) in the validation and test datasets were Support Vector Machine (SVM), Cubist (CUB) and Stochastic Gradient Boosting (SGB). For most models and cover crops, the Enhanced Vegetation Index (EVI) and the Normalized Difference Vegetation Index (NDVI) were the two most important AGB predictors. At the date of maximum VI at the beginning of the flowering stage, the correlation coefficients (r) between the cover-crop AGB and the cash-crop yield (soybean and maize) ranged from +0.70 for forage turnip to +0.78 for rye. The fuzzy unsupervised classification of the cover-crop AGB estimates delineated two management zones, which were spatially consistent with those obtained from cash-crop yield. The comparison between both maps produced overall accuracies that ranged from 61.20% to 68.25% with zone 2 having higher cover-crop AGB and cash-crop yield than zone 1 over the three sites. We conclude that satellite AGB estimates of cover crops can be used as a proxy for generating management zone maps in agricultural fields. These maps can be further refined in the field with any other type of method and data, whenever necessary.     

Combining UAV-based plant height from crop surface models,visible, and near infrared vegetation indices for biomass monitoring in barley

In this study we combined selected vegetation indices (VIs) and plant height information to estimate biomass in a summer barley experiment. The VIs were calculated from ground-based hyperspectral data and unmanned aerial vehicle (UAV)-based red green blue (RGB) imaging. In addition, the plant height information was obtained from UAV-based multi-temporal crop surface models (CSMs). The test site is a summer barley experiment comprising 18 cultivars and two nitrogen treatments located in Western Germany. We calculated five VIs from hyperspectral data. The normalised ratio index (NRI)-based index GnyLi (Gnyp et al., 2014) showed the highest correlation (R² = 0.83) with dry biomass. In addition, we calculated three visible band VIs: the green red vegetation index (GRVI), the modified GRVI (MGRVI) and the red green blue VI (RGBVI), where the MGRVI and the RGBVI are newly developed VI. We found that the visible band VIs have potential for biomass prediction prior to heading stage. A robust estimate for biomass was obtained from the plant height models (R² = 0.80–0.82). In a cross validation test, we compared plant height, selected VIs and their combination with plant height information. Combining VIs and plant height information by using multiple linear regression or multiple non-linear regression models performed better than the VIs alone. The visible band GRVI and the newly developed RGBVI are promising but need further investigation. However, the relationship between plant height and biomass produced the most robust results. In summary, the results indicate that plant height is competitive with VIs for biomass estimation in summer barley. Moreover, visible band VIs might be a useful addition to biomass estimation. The main limitation is that the visible band VIs work for early growing stages only.     

Analysis of common canopy vegetation indices for indicating leaf nitrogen accumulations in wheat and rice

The common spectra wavebands and vegetation indices (VI) were identified for indicating leaf nitrogen accumulation (LNA), and the quantitative relationships of LNA to canopy reflectance spectra were determined in both wheat (Triticum aestivum L.) and rice (Oryza sativa L.). The 810 and 870 nm are two common spectral wavebands indicating LNA in both wheat and rice. Among all ratio vegetation indices (RVI), difference vegetation indices (DVI) and normalized difference vegetation indices (NDVI) of 16 wavebands from the MSR16 radiometer, RVI (870, 660) and RVI (810, 660) were most highly correlated to LNA in both wheat and rice. In addition, the relations between VIs and LNA gave better results than relations between single wavebands and LNA in both wheat and rice. Thus LNA in both wheat and rice could be indicated with common VIs, but separate regression equations are better for LNA monitoring.     

Using low-density discrete Airborne Laser Scanning data to assess the potential carbon dioxide emission in case of a fire event in a Mediterranean pine forest

The aim of study is to map the carbon dioxide (CO₂) emission of the aboveground tree biomass (AGB) in case of a fire event. The suitability of low point density, discrete, multiple-return, Airborne Laser Scanning (ALS) data and the influence of several characteristics of these data and the study area on the results obtained have been evaluated. A sample of 45 circular plots representative of Pinus halepensis Miller stands were used to fit and validate the model of AGB. The ALS point clouds were processed to obtain the independent variables and a multivariate linear regression analysis between field data and ALS-derived variables allowed estimation of AGB. Then, the influence of several characteristics on the residuals of the model was analyzed. Finally, conversion factors were applied to obtain the CO₂ values. The AGB model presented a R² value of 0.84 with a relative root-mean-square error of 27.35%. This model included ALS variables related to vegetation height variability and to canopy density. Terrain slope, aspect, canopy cover, scan angle and the number of laser returns did not influence AGB estimations at plot level.     

Sensitivity of narrowband vegetation indices to boreal forest LAI,reflectance seasonality and species composition

There is growing evidence that imaging spectroscopy could improve the accuracy of satellite-based retrievals of vegetation attributes, such as leaf area index (LAI) and biomass. In this study, we evaluated narrowband vegetation indices (VIs) for estimating overstory effective LAI (LAI_eff) in a southern boreal forest area for the period between the end of snowmelt and maximum LAI using three Hyperion images and concurrent field measurements. We compared the performance of narrowband VIs with two SPOT HRVIR images, which closely corresponded to the imaging dates of the Hyperion data, and with synthetic broadband VIs computed from Hyperion images. According to the results, narrowband VIs based on near infrared (NIR) bands, and NIR and shortwave infrared (SWIR) bands showed the strongest linear relationships with LAI_eff over its typical range of variation and for the studied period of the snow-free season. The relationships were not dependent on dominant tree species (coniferous vs. broadleaved), which is an advantage in heterogeneous boreal forest landscapes. The best VIs, particularly those based on NIR spectral bands close to the 1200 nm liquid water absorption feature, provided a clear improvement over the best broadband VIs.     

Investigation of terrain illumination effects on vegetation indices and VI-derived phenological metrics in subtropical deciduous forests

We used RapidEye and Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra data to study terrain illumination effects on 3 vegetation indices (VIs) and 11 phenological metrics over seasonal deciduous forests in southern Brazil. We applied TIMESAT for the analysis of the Enhanced Vegetation Index (EVI) and the Normalized Difference Vegetation Index (NDVI) derived from the MOD13Q1 product to calculate phenological metrics. We related the VIs with the cosine of the incidence angle i (Cos i) and inspected percentage changes in VIs before and after topographic C-correction. The results showed that the EVI was more sensitive to seasonal changes in canopy biophysical attributes than the NDVI and Red-Edge NDVI, as indicated by analysis of non-topographically corrected RapidEye images from the summer and winter. On the other hand, the EVI was more sensitive to terrain illumination, presenting higher correlation coefficients with Cos i that decreased with reduction in the canopy background L factor. After C-correction, the RapidEye Red-Edge NDVI, NDVI, and EVI decreased 2%, 1%, and 13% over sunlit surfaces and increased up to 5%, 14%, and 89% over shaded surfaces, respectively. The EVI-related phenological metrics were also much more affected by topographic effects than the NDVI-derived metrics. From the set of 11 metrics, the 2 that described the period of lower photosynthetic activity and seasonal VI amplitude presented the largest correlation coefficients with Cos i. The results showed that terrain illumination is a factor of spectral variability in the seasonal analysis of phenological metrics, especially for VIs that are not spectrally normalized.     

Investigating the impact of the temporal resolution of MODIS data on measured phenology in the prairie grasslands

ABSTRACT

The temporal resolution of vegetation indices (VIs) determines the details of seasonal variation in vegetation dynamics observed by remote sensing, but little has been known about how the temporal resolution of VIs affects the retrieval of land surface phenology (LSP) of grasslands. This study evaluated the impact of temporal resolution of MODIS NDVI, EVI, and per-pixel green chromatic coordinate (GCCpp) on the quality and accuracy of the estimated LSP metrics of prairie grasslands. The near-surface PheonoCam phenology data for grasslands centered over Lethbridge PhenoCam grassland site were used as the validation datasets due to the lack of in situ observations for grasslands in the Prairie Ecozone. MODIS Nadir Bidirectional Reflectance Distribution Function (BRDF)-Adjusted Reflectance (NBAR) data from 2001 to 2017 were used to compute the time series of daily reference and to simulate 2–32 day MODIS VIs. The daily reference and simulated multi-day time series were fitted with the double logistic model, and the LSP metrics were then retrieved from the modeled daily time series separately. Comparison within satellite-based estimates showed no significant difference in the phenological metrics derived from daily reference and multi-day VIs resampled at a time step less than 18 days. Moreover, a significant decline in the ability of multi-day VIs to predict detailed temporal dynamics of daily reference VIs was revealed as the temporal resolution increased. Besides, there were a variety of trends for the onset of phenological transitions as the temporal resolution of VIs changed from 1 to 32 days. Comparison with PhenoCam phenology data presented small and insignificant differences in the mean bias error (MBE) and the mean absolute error (MAE) of grassland phenological metrics derived from daily, 8-, 10-, 14-, and 16-day MODIS VIs. Overall, this study suggested that the MODIS VIs resampled at a time step less than 18 days are favorable for the detection of grassland phenological transitions and detailed seasonal dynamics in the Prairie Ecozone.     

Comparing methods for mapping canopy chlorophyll content in a mixed mountain forest using Sentinel-2 data

The Sentinel-2 Multi-Spectral Imager (MSI) has three spectral bands centered at 705, 740, and 783 nm wavelengths that exploit the red-edge information useful for quantifying plant biochemical traits. This sensor configuration is expected to improve the prediction accuracy of vegetation chlorophyll content. In this work, we assessed the performance of several statistical and physical-based methods in retrieving canopy chlorophyll content (CCC) from Sentinel-2 in a heterogeneous mixed mountain forest. Amongst the algorithms presented in the literature, 13 different vegetation indices (VIs), a non-parametric statistical approach, and two radiative transfer models (RTM) were used to assess the CCC prediction accuracy. A field campaign was conducted in July 2017 to collect in situ measurements of CCC in Bavarian forest national park, and the cloud-free Sentinel-2 image was acquired on 13 July 2017. The leave-one-out cross-validation technique was used to compare the VIs and the non-parametric approach. Whereas physical-based methods were calibrated using simulated data and validated using the in situ reference dataset. The statistical-based approaches, such as the modified simple ratio (mSR) vegetation index and the partial least square regression (PLSR) outperformed all other techniques. As such the modified simple ratio (mSR3) (665, 865) gave the lowest cross-validated RMSE of 0.21 g/m² (R² = 0.75). The PLSR resulted in the highest R² of 0.78, and slightly higher RMSE =0.22 g/m² than mSR3. The physical-based approach-INFORM inversion using look-up table resulted in an RMSE =0.31 g/m², and R² = 0.67. Although mapping CCC using these methods revealed similar spatial distribution patterns, over and underestimation of low and high CCC values were observed mainly in the statistical approaches. Further validation using in situ data from different terrestrial ecosystems is imperative for both the statistical and physical-based approaches' effectiveness to quantify CCC before selecting the best operational algorithm to map CCC from Sentinel-2 for long-term terrestrial ecosystems monitoring across the globe.     

A hyperspectral experiment over tropical forests based on the EO-1 orbit change and PROSAIL simulation

ABSTRACT

We designed a unique hyperspectral experiment from the Earth Observing One (EO-1) orbit change to evaluate solar illumination effects over tropical forests in Brazil. Ten nadir-viewing Hyperion images collected over a fixed site and period of the year (July to August) were selected for analysis. We evaluated variations in reflectance and in 16 narrowband vegetation indices (VIs) with increasing solar zenith angle (SZA) from the pre-drift (2004–2008) to the EO-1 drift period (2011–2016). To detect changes in reflectance and shadows, we applied spectral mixture analysis (SMA) and principal component analysis (PCA) and calculated the similarity spectral angle (θ) between the vegetation spectra measured with variable SZA. The magnitude of the illumination effects was also evaluated from change-point analysis and nonparametric Mann-Whitney U tests applied over the time series. Finally, we complemented our experiment using the PROSAIL model to simulate the VIs variation with increasing SZA resultant from satellite drift. The results showed significant changes in Hyperion reflectance and VIs, especially when the EO-1 crossed the study area at earlier times and larger SZA in 2015 (9:05 a.m.; SZA = 59°) and 2016 (8:30 a.m.; SZA = 67°). Compared to the pre-drift period (10:30 a.m.; SZA = 45°), the SZA differences of 14° (2015) and 22° (2016) increased the shade fractions and decreased the vegetation brightness. PCA separated the pre-drift and drift reflectance datasets, showing shifts in scores due to changes in brightness. θ increased with SZA, indicating changes in the shape of the vegetation spectra with drift. For most VIs, the change-point analysis indicated 2015 (SZA = 59°) as the predominant year of detected changes. Compared to the EO-1 original orbit, the Plant Senescence Reflectance Index (PSRI), Anthocyanin Reflectance Index (ARI) and Structure Insensitive Pigment Index (SIPI) presented the largest positive changes during drift, while the Photochemical Reflectance Index (PRI), Visible Atmospherically Resistant Index (VARI) and Enhanced Vegetation Index (EVI) had the largest negative changes. The effect size of the illumination geometry on these VIs was large, as indicated by increasing values of the Cohen’s r metric toward 2016. The anisotropy of the Hyperion VIs was generally consistent with that from PROSAIL in the simulated pre-drift and drift periods. Focusing on structural indices, it affected the relationships between VIs and simulated leaf area index (LAI) at large SZA.     

Estimation of vegetation LAI from hyperspectral reflectance data: Effects of soil type and plant architecture

The retrieval of canopy biophysical variables is known to be affected by confounding factors such as plant type and background reflectance. The effects of soil type and plant architecture on the retrieval of vegetation leaf area index (LAI) from hyperspectral data were assessed in this study. In situ measurements of LAI were related to reflectances in the red and near-infrared and also to five widely used spectral vegetation indices (VIs). The study confirmed that the spectral contrast between leaves and soil background determines the strength of the LAI–reflectance relationship. It was shown that within a given vegetation species, the optimum spectral regions for LAI estimation were similar across the investigated VIs, indicating that the various VIs are basically summarizing the same spectral information for a given vegetation species. Cross-validated results revealed that, narrow-band PVI was less influenced by soil background effects (0.15 ≤ RMSE_cv ≤ 0.56). The results suggest that, when using remote sensing VIs for LAI estimation, not only is the choice of VI of importance but also prior knowledge of plant architecture and soil background. Hence, some kind of landscape stratification is required before using hyperspectral imagery for large-scale mapping of vegetation biophysical variables.     

Integrating temperature vegetation dryness index (TVDI) and regional water stress index (RWSI) for drought assessment with the aid of LANDSAT TM/ETM+ images

This paper presents a new drought assessment method by spatially and temporally integrating temperature vegetation dryness index (TVDI) with regional water stress index (RWSI) based on a synergistic approach. With the aid of LANDSAT TM/ETM data, we were able to retrieve the land-use and land-cover (LULC), vegetation indices (VIs), and land surface temperature (LST), leading to the derivation of three types of modified TVDI, including TVDI_SAVI, TVDI_ANDVI and TVDI_MSAVI, for drought assessment in a fast growing coastal area, Northern China. The categorical classification of four drought impact levels associated with the RWSI values enables us to refine the spatiotemporal relationship between the LST and the VIs. Holistic drought impact assessment between 1987 and 2000 was carried out by linking RWSI with TVDIs group wise. Research findings indicate that: (1) LST and VIs were negatively correlated in most cases of low, medium, and high vegetation cover except the case of high density vegetation cover in 2000 due to the effect of urban heat island (UHI) effect; (2) the shortage of water in 1987 was more salient than that that in 2000 based on all indices of TVDI and RWSI; and (3) TVDIs are more suitable for monitoring mild drought, normal and wet conditions when RWSI is smaller than 0.752; but they are not suitable for monitoring moderate and severe drought conditions.     

Historical forest biomass dynamics modelled with Landsat spectral trajectories

Estimation of forest aboveground biomass (AGB) is informative of the role of forest ecosystems in local and global carbon budgets. There is a need to retrospectively estimate biomass in order to establish a historical baseline and enable reporting of change. In this research, we used temporal spectral trajectories to inform on forest successional development status in support of modelling and mapping of historic AGB for Mediterranean pines in central Spain. AGB generated with ground plot data from the Spanish National Forest Inventory (NFI), representing two collection periods (1990 and 2000), are linked with static and dynamic spectral data as captured by Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) sensors over a 25 year period (1984–2009). The importance of forest structural complexity on the relationship between AGB and spectral vegetation indices is revealed by the analysis of wavelet transforms. Two-dimensional (2D) wavelet transforms support the identification of spectral trajectory patterns of forest stands that in turn, are associated with traits of individual NFI plots, using a flexible algorithm sensitive to capturing time series similarity. Single-date spectral indices, temporal trajectories, and temporal derivatives associated with succession are used as input variables to non-parametric decision trees for modelling, estimation, and mapping of AGB and carbon sinks over the entire study area. Results indicate that patterns of change found in Normalized Difference Vegetation Index (NDVI) values are associated and relate well to classes of forest AGB. The Tasseled Cap Angle (TCA) index was found to be strongly related with forest density, although the related patterns of change had little relation with variability in historic AGB. By scaling biomass models through small (∼2.5 ha) spatial objects defined by spectral homogeneity, the AGB dynamics in the period 1990–2000 are mapped (70% accuracy when validated with plot values of change), revealing an increase of 18% in AGB irregularly distributed over 814 km² of pines. The accumulation of C calculated in AGB was on average 0.65 t ha⁻¹ y⁻¹, equivalent to a fixation of 2.38 t ha⁻¹ y⁻¹ of carbon dioxide.     

机载激光雷达及高光谱的森林乔木物种多样性遥感监测

利用机载LiDAR和高光谱数据并结合37个地面调查样本数据,基于结构差异与光谱变异理论,通过相关分析法分别筛选了3个最优林冠结构参数和6个最优光谱指数,在单木尺度上利用自适应C均值模糊聚类算法,在神农架国家自然保护区开展森林乔木物种多样性监测,实现了森林乔木物种多样性的区域成图。研究结果表明,(1)基于结合形态学冠层控制的分水岭算法可以获得较高精度的单木分割结果(R~2=0.88,RMSE=13.17,P0.001);(2)基于LiDAR数据提取的9个结构参数中,95%百分位高度、冠层盖度和植被穿透率为最优结构参数,与Shannon-Wiener指数的相关性达到R~2=0.39—0.42(P0.01);(3)基于机载高光谱数据筛选的16个常用的植被指数中,CRI、OSAVI、Narrow band NDVI、SR、Vogelmann index1、PRI与Shannon-Wiener指数的相关性最高(R~2=0.37—0.45,P0.01);(4)在研究区,利用以30 m×30 m为窗口的自适应模糊C均值聚类算法可预测的最大森林乔木物种数为20,物种丰富度的预测精度为R~2=0.69,RMSE=3.11,Shannon-Wiener指数的预测精度为R~2=0.70,RMSE=0.32。该研究在亚热带森林开展乔木物种多样性监测,是在区域尺度上进行物种多样性成图的重要实践,可有效补充森林生物多样性本底数据的调查手段,有助于实现生物多样性的长期动态监测及科学分析森林物种多样性的现状和变化趋势。