Predicting forest carbon stocks from high resolution satellite data in dry forests of Zimbabwe: exploring the effect of the red-edge band in forest carbon stocks estimation

In this study, we tested whether the inclusion of the red-edge band as a covariate to vegetation indices improves the predictive accuracy in forest carbon estimation and mapping in savanna dry forests of Zimbabwe. Initially, we tested whether and to what extent vegetation indices (simple ratio SR, soil-adjusted vegetation index and normalized difference vegetation index) derived from high spatial resolution satellite imagery (WorldView-2) predict forest carbon stocks. Next, we tested whether inclusion of reflectance in the red-edge band as a covariate to vegetation indices improve the model's accuracy in forest carbon prediction. We used simple regression analysis to determine the nature and the strength of the relationship between forest carbon stocks and remotely sensed vegetation indices. We then used multiple regression analysis to determine whether integrating vegetation indices and reflection in the red-edge band improve forest carbon prediction. Next, we mapped the spatial variation in forest carbon stocks using the best regression model relating forest carbon stocks to remotely sensed vegetation indices and reflection in the red-edge band. Our results showed that vegetation indices alone as an explanatory variable significantly (p < 0.05) predicted forest carbon stocks with R² ranging between 45 and 63% and RMSE ranging from 10.3 to 12.9%. However, when the reflectance in the red-edge band was included in the regression models the explained variance increased to between 68 and 70% with the RMSE ranging between 9.56 and 10.1%. A combination of SR and reflectance in the red edge produced the best predictor of forest carbon stocks. We concluded that integrating vegetation indices and reflectance in the red-edge band derived from high spatial resolution can be successfully used to estimate forest carbon in dry forests with minimal error.     

Validation of MODIS Vegetation Continuous Fields for monitoring deforestation and forest degradation: two cases in Mexico

This study assesses whether MODIS Vegetation Continuous Fields percent tree cover (PTC) data can detect deforestation and forest degradation. To assess the usefulness of PTC for detecting deforestation, we used a data set consisting of eight forest and seven non-forest categories. To evaluate forest degradation, we used data from two temperate forest types in three conservation states: primary (dense), secondary (moderately degraded) and open (heavily degraded) forest. Our results show that PTC can differentiate temperate forest from non-forest categories (p = 0.05) and thus suggests PTC can adequately detect deforestation in temperate forests. In contrast, single-date PTC data does not appear to be adequate to detect forest degradation in temperate forests. As for tropical forest, PTC can partially discriminate between forest and non-forest categories.     

Spectral anisotropy of subtropical deciduous forest using MISR and MODIS data acquired under large seasonal variation in solar zenith angle

Recent studies in Amazonian tropical evergreen forests using the Multi-angle Imaging SpectroRadiometer (MISR) and the Moderate Resolution Imaging Spectroradiometer (MODIS) have highlighted the importance of considering the view-illumination geometry in satellite data analysis. However, contrary to the observed for evergreen forests, bidirectional effects have not been evaluated in Brazilian subtropical deciduous forests. In this study, we used MISR data to characterize the reflectance and vegetation index anisotropies in subtropical deciduous forest from south Brazil under large seasonal solar zenith angle (SZA) variation and decreasing leaf area index (LAI) from the summer to winter. MODIS data were used to observe seasonal changes in the normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI). Topographic effects on their determination were inspected by dividing data from the summer to winter and projecting results over a digital elevation model (DEM). By using the PROSAIL, we investigated the relative contribution of LAI and SZA to vegetation indices (VI) of deciduous forest. We also simulated and compared the MISR NDVI and EVI response of subtropical deciduous and tropical evergreen forests as a function of the large seasonal SZA amplitude of 33°. Results showed that the MODIS-MISR NDVI and EVI presented higher values in the summer and lower ones in the winter with decreasing LAI and increasing SZA or greater amounts of canopy shadows viewed by the sensors. In the winter, NDVI reduced local topographic effects due to the red-near infrared (NIR) band normalization. However, the contrary was observed for the three-band EVI that enhanced local variations in shaded and sunlit surfaces due to its strong dependence on the NIR band response. The reflectance anisotropy of the MISR bands increased from the summer to winter and was stronger in the backscattering direction at large view zenith angles (VZA). EVI was much more anisotropic than NDVI and the anisotropy increased from the summer to winter. It also increased from the forward scatter to the backscattering direction with the predominance of sunlit canopy components viewed by MISR, especially at large VZA. Modeling PROSAIL results confirmed the stronger anisotropy of EVI than NDVI for the subtropical deciduous and tropical evergreen forests. PROSAIL showed that LAI and SZA are coupled factors to decrease seasonally the VIs of deciduous forest with the first one having greater importance than the latter. However, PROSAIL seasonal variations in VIs were much smaller than those observed with MODIS data probably because the effects of shadows in heterogeneous canopy structures or/and cast by emergent trees and from local topography were not modeled.     

Estimating forest standing biomass in savanna woodlands as an indicator of forest productivity using the new generation WorldView-2 sensor

Accurate and up-to-date information on forest dendrometric traits, such as above ground biomass is important in understanding the contribution of terrestrial ecosystems to the regulation of atmsopheric carbon, especially in the face of global environmental change. Besides, dendrometric traits information is critical in assessing the healthy and the spatial planning of fragile ecosystems, such as the savanna dry forests. The aim of this work was to test whether red-edge spectral data derived from WorldView-2 multispectral imagery improve biomass estimation in savanna dry forests. The results of this study have shown that biomass estimation using all Worldview-2 raw spectral bands without the red-edge band yielded low estimation accuracies (R² of 0.67 and a RMSE-CV of 2.2 t ha^?1) when compared to when the red-edge band was included as a co-variate (R² of 0.73 and a RMSE-CV of 2.04 t ha^?1). Also, similar results were obseved when all WorldView-2 vegetation indices (without the red-edge computed ones), producing slightly low accuracies (R² of about 0.67 and a RMSE-CV of 2.20 t ha^?1), when compared to those obtained using all indices and RE-computed indices(R² of 0.76 and a RMSE-CV of 1.88 t ha^?1). Overall, the findings of this work have demontrated the potential and importance of strategically positioned bands, such as the red-edge band in the optimal estimation of indigeonus forest biomass. These results underscores the need to shift towards embracing sensors with unique and strategeically positioned bands, such as the forthcoming Sentinel 2 MSI and HysPIRI which have a global footprint.     

Elephants move faster in small fragments of low productivity in Amboseli ecosystems: Kenya

Understanding factors affecting the behaviour and movement patterns of the African elephant is important for wildlife conservation, especially in increasingly human-dominated savanna landscapes. Currently, knowledge on how landscape fragmentation and vegetation productivity affect elephant speed of movement remains poorly understood. In this study, we tested whether landscape fragmentation and vegetation productivity explains elephant speed of movement in the Amboseli ecosystem in Kenya. We used GPS collar data from five elephants to quantify elephant speed of movement for three seasons (wet, dry and transitional). We then used multiple regression to model the relationship between speed of movement and landscape fragmentation, as well as vegetation productivity for each season. Results of this study demonstrate that landscape fragmentation and vegetation productivity predicted elephant speed of movement poorly (R² < 0.4) when used as solitary covariates. However, a combination of the covariates significantly (p < 0.05) explained variance in elephant speed of movement with improved R² values of 0.69, 0.45, 0.47 for wet, transition and dry seasons, respectively.     

Estimation of aboveground biomass and carbon in a tropical rain forest in Gabon using remote sensing and GPS data

The knowledge of biomass stocks in tropical forests is critical for climate change and ecosystem services studies. This research was conducted in a tropical rain forest located near the city of Libreville (the capital of Gabon), in the Akanda Peninsula. The forest cover was stratified in terms of mature, secondary and mangrove forests using Landsat-ETM data. A field inventory was conducted to measure the required basic forest parameters and estimate the aboveground biomass (AGB) and carbon over the different forest classes. The Shuttle Radar Topography Mission (SRTM) data were used in combination with ground-based GPS measurements to derive forest heights. Finally, the relationships between the estimated heights and AGB were established and validated. Highest biomass stocks were found in the mature stands (223 ± 37 MgC/ha), followed by the secondary forests (116 ± 17 MgC/ha) and finally the mangrove forests (36 ± 19 MgC/ha). Strong relationships were found between AGB and forest heights (R² > 0.85).     

The effects of land cover change on carbon stock dynamics in a dry Afromontane forest in northern Ethiopia

Background

Forests play an important role in mitigating global climate change by capturing and sequestering atmospheric carbon. Quantitative estimation of the temporal and spatial pattern of carbon storage in forest ecosystems is critical for formulating forest management policies to combat climate change. This study explored the effects of land cover change on carbon stock dynamics in the Wujig Mahgo Waren forest, a dry Afromontane forest that covers an area of 17,000 ha in northern Ethiopia.

Results

The total carbon stocks of the Wujig Mahgo Waren forest ecosystems estimated using a multi-disciplinary approach that combined remote sensing with a ground survey were 1951, 1999, and 1955 GgC in 1985, 2000 and 2016 years respectively. The mean carbon stocks in the dense forests, open forests, grasslands, cultivated lands and bare lands were estimated at 181.78?±?27.06, 104.83?±?12.35, 108.77?±?6.77, 76.54?±?7.84 and 83.11?±?8.53 MgC ha^?1 respectively. The aboveground vegetation parameters (tree density, DBH and height) explain 59% of the variance in soil organic carbon.

Conclusions

The obtained estimates of mean carbon stocks in ecosystems representing the major land cover types are of importance in the development of forest management plan aimed at enhancing mitigation potential of dry Afromontane forests in northern Ethiopia.

     

Warming trends in Patagonian subantartic forest

The forests in the Aysén region (ca. 43–49 °S, Chile) have a high degree of wilderness and cover more than 4.8 million hectares, making it one of the largest areas of subantarctic forest in the Southern Hemisphere. The impact of global warming on this region is poorly documented. The main objective of this work was to analyze the normalized difference vegetation index (NDVI), land surface temperature (LST) and precipitation over Aysen forests in the context of ongoing global warming. We used average monthly images of LST and NDVI derived from the MODIS sensor covering the period 2001–2016 and precipitation from gridded datasets. The Aysén region was divided into three nested spatial scales: i) regional, ii) regional considering only forests, iii) local scale considering an evergreen subantarctic forest area covering around 5 × 5 km and a local deciduous forest area (dominated by Nothofagus pumilio). Trend analysis showed a warming rate of +0.78 K/decade (p ≤ 0.05) over the subantarctic forest zone, greening of +0.01/decade for NDVI (p ≤ 0.05) over the western zone, and a drying trend (p ≤ 0.05) over the eastern zone. The minimum temperature anomalies showed an increase of about 4.5 K during the period under analysis. LST, NDVI and precipitation were also analyzed here. The recent trends in temperature, greening and precipitation over the forests of Aysén detected in this research contribute to a better understanding of global warming impacts on subantarctic forests in the southern tip of South America. Nevertheless, to get a better estimation of the impact of global warming at multiple scales is needed to have better quality and quantity of data in situ.     

View-illumination effects on hyperspectral vegetation indices in the Amazonian tropical forest

Because of the pointing capability of the Hyperion/Earth Observing-One (EO-1) to improve the revisit time of the scene, temporal series of narrowband vegetation indices (VIs) can be generated to study the phenology of the Amazonian tropical forests. In this study, 10 selected narrowband VIs calculated from Hyperion nadir and off-nadir data and from different view directions (forward scattering and backscattering) were analyzed for their sensitivity to view-illumination effects along the dry season on the Seasonal Semi-deciduous Forest. Data analysis was also supported by PROSAIL modeling to simulate the spectral response of this forest type in both directions. Hyperion and PROSAIL results showed that the Enhanced Vegetation Index (EVI) and Photochemical Reflectance Index (PRI) were the two more anisotropic VIs, whereas the Normalized Difference Vegetation Index (NDVI), Structure Insensitive Pigment Index (SIPI) and the Vogelmann Red Edge Index (VOG) were comparatively less sensitive to view-illumination effects. When compared to the other VIs and because of the greater dependence on the near-infrared (NIR) reflectance, EVI showed a different spectral behavior. EVI increased from forward scattering to backscattering and with decreasing solar zenith angle (SZA) towards the end of the local dry season, due to reduction in shading and enhancement of the illumination effects. On the other hand, PRI was higher with increasing shading in the forward scattering direction, as deduced from the PROSAIL simulation. Results emphasized the importance of taking into account bidirectional effects when analyzing temporal series of VIs collected over tropical forests by imaging spectrometers with pointing capability or even by multispectral sensors with large field-of-view (FOV).     

Seasonal trend analysis (STA) of MODIS vegetation index time series for the mangrove canopy of the Teacapan-Agua Brava lagoon system,Mexico

Monthly time series, from 2001 to 2016, of the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) from MOD13Q1 products were analyzed with Seasonal Trend Analysis (STA), assessing seasonal and long-term changes in the mangrove canopy of the Teacapan-Agua Brava lagoon system, the largest mangrove ecosystem in the Mexican Pacific coast. Profiles from both vegetation indices described similar phenological trends, but the EVI was more sensitive in detecting intra-annual changes. We identified a seasonal cycle dominated by Laguncularia racemosa and Rhizophora mangle mixed patches, with the more closed canopy occurring in the early autumn, and the maximum opening in the dry season. Mangrove patches dominated by Avicennia germinans displayed seasonal peaks in the winter. Curves fitted for the seasonal vegetation indices were better correlated with accumulated precipitation and solar radiation among the assessed climate variables (Pearson’s correlation coefficients, estimated for most of the variables, were r ≥ 0.58 p < 0.0001), driving seasonality for tidal basins with mangroves dominated by L. racemosa and R. mangle. For tidal basins dominated by A. germinans, the maximum and minimum temperatures and monthly precipitation fit better seasonally with the vegetation indices (r ≥ 0.58, p < 0.0001). Significant mangrove canopy reductions were identified in all the analyzed tidal basins (z values for the Mann-Kendall test ≤ ?1.96), but positive change trends were recorded in four of the basins, while most of the mangrove canopy (approximately 87%) displayed only seasonal canopy changes or canopy recovery (z > ?1.96). The most resilient mangrove forests were distributed in tidal basins dominated by L. racemosa and R. mangle (Mann-Kendal Tau t ≥ 0.4, p ≤ 0.03), while basins dominated by A. germinans showed the most evidence of disturbance.     

The improvement of an object-oriented classification using multi-temporal MODIS EVI satellite data

Abstract

This paper investigates the contribution of multi-temporal enhanced vegetation index (EVI) data to the improvement of object-based classification accuracy using multi-spectral moderate resolution imaging spectral-radiometer (MODIS) imagery. In object-oriented classification, similar pixels are firstly grouped together and then classified; the produced result does not suffer the speckled appearance and closer to human vision. EVI data are from the MODIS sensor aboard Terra spacecraft. 69 EVI data (scenes) were collected during the period of three years (2001–2003) in a mountainous vegetated area. These data sets were used to study the phenology of the land cover types. Different land cover types show distinct fluctuations over time in EVI values and this information might be used to improve object-oriented land cover classification. Two experiments were carried out: one was only with single date MODIS multispectral data, and the other one including also the 69 EVI images. Eight classes were distinguished: temperate forest, tropical dry forest, grassland, irrigated agriculture, rain-fed agriculture, orchards, lava flows and human settlement. The two classifications were evaluated with independent verification data, and the results showed that with multi-temporal EVI data, the classification accuracy was improved 5.2%. Evaluated by McNemar's test, this improved was significant, with significance level p=0.01.     

Effect of flying altitude and pulse repetition frequency on laser scanner penetration rate for digital elevation model generation in a tropical forest

In tropical forests, the penetration ability of airborne laser scanning (ALS) may be limited because of highly dense vegetation cover. However, in the typical planning of ALS surveys, the ability of laser pulses to penetrate forests is not considered. Nine round-trip flight lines covering the area of a tropical forest on the northeast side of the Tsengwen Reservoir in Taiwan were designed in this study. Five flight lines flew at altitudes of 1.525, 1.830, 2.135, 2.440, and 2.745 km, and the other four had pulse repetition frequencies (PRFs) of 100, 150, 200, and 250 kHz. The laser penetration index (LPI) is a quantitative index measuring the penetration ability of the ALS and consists of the ratio of the number of laser pulses reaching the forest floor to the total number of laser pulses. The LPI was used to represent the laser penetration rate and investigate the influence of flying altitude and PRF on the LPI. The results showed that as the flying altitude decreased by 1 km, the average LPI increased by 10%, and as the PRF decreased by 50 kHz, the average LPI increased by 2%. The effect of the LPI on digital elevation models (DEMs) was confirmed by visual images obtained by DEMs at five altitudes. The DEM obtained at an altitude of 2.745 km was coarsely textured, whereas that obtained at an altitude of 1.525 km was finely textured. The in-situ height data obtained from the electronic Global Navigation Satellite System (eGNSS) were compared with the data of the ALS-generated DEMs. The results indicated that when the LPI ≥60%, the height difference between the in situ data and DEM data was not prominent. However, when the LPI <60%, the ALS-derived DEM could be higher or lower than the in-situ height; the largest difference between the two was 1.7 m. The LPI of a forest should be considered for ALS survey planning, especially when consistent DEM precision for large tropical forest areas is paramount.     

Biomass allometric equation and expansion factor for a mountain moist evergreen forest in Mozambique

Background

Worldwide, forests are an important carbon sink and thus are key to mitigate the effects of climate change. Mountain moist evergreen forests in Mozambique are threatened by agricultural expansion, uncontrolled logging, and firewood collection, thus compromising their role in carbon sequestration. There is lack of local tools for above-ground biomass (AGB) estimation of mountain moist evergreen forest, hence carbon emissions from deforestation and forest degradation are not adequately known. This study aimed to develop biomass allometric equations (BAE) and biomass expansion factor (BEF) for the estimation of total above-ground carbon stock in mountain moist evergreen forest.

Methods

The destructive method was used, whereby 39 trees were felled and measured for diameter at breast height (DBH), total height and the commercial height. We determined the wood basic density, the total dry weight and merchantable timber volume by Smalian’s formula. Six biomass allometric models were fitted using non-linear least square regression. The BEF was determined based on the relationship between bole stem dry weight and total dry weight of the tree. To estimate the mean AGB of the forest, a forest inventory was conducted using 27 temporary square plots. The applicability of Marzoli’s volume equation was compared with Smalian’s volume equation in order to check whether Marzoli’s volume from national forest inventory can be used to predict AGB using BEF.

Results

The best model was the power model with only DBH as predictor variable, which provided an estimated mean AGB of 291?±?141 Mg ha^?1 (mean?±?95% confidence level). The mean wood basic density of sampled trees was 0.715?±?0.182 g cm^?3. The average BEF was of 2.05?±?0.15 and the estimated mean AGB of 387?±?126 Mg ha^?1. The BAE from miombo woodland within the vicinity of the study area underestimates the AGB for all sampled trees. Chave et al.’s pantropical equation of moist forest did not fit to the Moribane Forest Reserve, while Brown’s equation of moist forest had a good fit to the Moribane Forest Reserve, having generated 1.2% of bias, very close to that generated by the selected model of this study. BEF showed to be reliable when combined with stand mean volume from Marzoli’s National Forestry Inventory equation.

Conclusion

The BAE and the BEF function developed in this study can be used to estimate the AGB of the mountain moist evergreen forests at Moribane Forest Reserve in Mozambique. However, the use of the biomass allometric model should be preferable when DBH information is available.

     

激光雷达和点云切片算法结合的森林有效叶面积指数估算

叶面积指数(LAI)是定量模拟和刻画植被冠层结构、生理过程以及研究森林生态系统碳水循环中物质和能量交换的重要生物物理参数之一。定量分析林冠元素的3维空间分布是准确估算森林叶面积指数的重要基础和关键步骤。本文利用地面激光雷达扫描系统(TLS),获取了不同树种、密度、年龄和空间分布的森林3维点云数据,利用"径向半球点云切片"和"点云法向量重建"算法分别计算森林的角度孔隙率和消光系数,进而得到森林冠层的有效叶面积指数。通过与利用传统光学仪器得到的结果比较发现:单站半球式激光雷达扫描得到森林样方尺度的结果与LAI-2200和数字半球摄影观测所得结果的相关性分别为R2=0.7084(N=9,p0.01)和0.7409(N=14,p0.01)。通过角度分辨率(LBA)和角度孔隙率(AGF)的敏感性分析,建议径向半球切片算法的角度分辨率应参考中心单站TLS半球式扫描所设置的采样间距。本文所采用的径向半球切片算法可以有效地利用单站半球式扫描所获取的3维点云进行森林冠层有效叶面积指数的定量估算。     

Retrieval of tropical forest biomass information from ALOS PALSAR data

Advanced Land Observing Satellite Phased Array L-band Synthetic Aperture Radar (ALOS PALSAR) data from different observation modes were analysed to determine (1) which observation mode most accurately retrieves tropical forest biomass information and (2) whether different modes, when considered together, yield improved results in comparison to identical data-sets analysed independently. We performed regression analysis to estimate above-ground forest biomass using PALSAR backscatter data for natural and planted forests in south-eastern Bangladesh. The coefficient of determination (r ²) was lower or equal to 0.499 (n = 70) when PALSAR data from different observation modes were separately considered, but increased sharply when one class (rubber) is dropped and average backscatter of fine beam single (FBS) and polrimetric (PLR) modes are used in the analysis. The results presented in this article are useful for both regional and global forest biomass inventories and fixing acquisition modes for planned L-band SAR missions.     

Vegetation response to intensive commercial horticulture and environmental changes within watersheds in central highlands,Kenya, using AVHRR NDVI data

Expansion and heterogeneous clustering of commercial horticulture within the central highlands of Kenya after the mid-1990s impact watersheds and the sustainable resource management. This is distressing since climate conditions for world horticultural regions are projected to change, making such farming extremely difficult and costly to the environment. To understand the scope of impact on vegetation, the study evaluated (1) interannual variability in averaged normalized difference vegetation index (NDVI); (2) trends in average annual NDVI before and after 1990 – the presumed onset of rapid horticulture; and (3) relationship between the average annual NDVI and large-scale commercial farms, population density, and mean annual rainfall in subwatersheds. Time-series analysis of long-term Global Inventory Modeling and Mapping Studies NDVI data were analyzed as indicator of vegetation condition. NDVI trends before 1990s (1982–1989) and after 1990s (1990–2006) were evaluated to determine the slope (sign), and the Spearman’s correlation coefficient (strength). Overall, results show considerable variations in vegetation condition due largely to mixed factors including intensive farming activities, drought, and rainfall variation. Statistical analysis shows significant differences in slopes before 1990 and after 1990 (p < 0.05 and p < 0.1 respectively). Negative (decline) trends were common after 1990, linked to increased commercial horticulture and related anthropogenic disturbances on land cover. There was decline in vegetation over densely populated subwatersheds, though low NDVI values in 1984 and 2000 were the effect of severe droughts. Understanding the linkage between vegetation responses to the effects of human-induced pressure at the subwatershed scale can help natural resource managers approach conservation measures more effectively.     

The influence correction of boreal forest vegetation on SRTM data

The digital elevation model based on SRTM is very convenient for a wide range of studies but requires correction due to the influence of forest vegetation. The present study was conducted to analyse the effect of boreal forests on altitudes, aspects and slopes calculated from the SRTM. A DEM based on topographic maps at 1:100 000 scale was used as a reference. The linear regression analysis showed low data correlation in forested areas. The presence of different types of forests and felling in the woods leads to a complex distribution of deviations from the SRTM. A simple correction method was proposed, using a forest mask, built according to Landsat, and forest heights indicated on the topographic maps. After correction, the correlation coefficient between the altitudes increased by 0.05–0.14, the share of matching aspects by 1–4% and the share of matching slopes by 2–8%.     

Sentinel-2卫星落叶松林龄信息反演

林龄结构信息能够有效反映区域森林群落不同生长阶段的固碳能力,对于评估森林生态系统的健康状况具有重要意义。本研究以中国温带典型优势树种落叶松林为研究对象,分别选择其芽萌动期、展叶期和落叶期时段的Sentinel-2影像,采用多元线性回归（MLR）、随机森林（RF）、支持向量机回归（SVR）、前馈反向传播神经网络（BP）以及多元自适应回归样条（MARS）等5种方法依次构建落叶松林龄反演模型。通过相关性分析首先确定最佳遥感反演物候期,并在此基础上根据相关性差异筛选出5个最优特征变量用于模型反演,分别为冠层含水量（CWC）,归一化水体指数（NDWI）,叶面积指数（LAI）,光合有效辐射吸收率（FAPAR）和植被覆盖度（FVC）。研究结果表明,展叶期为落叶松林最佳遥感反演物候期。除植被衰减指数（PSRI）以及落叶期的NDVI、RVI外,落叶松林龄与各指标之间均呈负相关关系,其中与冠层含水量（CWC）的相关性最高,pearson相关系数达到-0.74（p<0.01）。此外,不同模型反演结果表明,随机森林模型（RF）为最佳落叶松林龄估测模型,其平均决定系数R²和平均均方根误差RMSE分别为0.89和2.91 a;多元线性回归模型（MLR）的林龄估测结果最差,其平均决定系数R²和平均均方根误差RMSE仅为0.57和5.69 a,非线性模型能更好的解释林龄与建模变量之间的关系。     

Augmenting forest inventory attributes with geometric optical modelling in support of regional susceptibility assessments to bark beetle infestations

Assessment of the susceptibility of forests to mountain pine beetle (Dendroctonus ponderosae Hopkins) infestation is based upon an understanding of the characteristics that predispose the stands to attack. These assessments are typically derived from conventional forest inventory data; however, this information often represents only managed forest areas. It does not cover areas such as forest parks or conservation regions and is often not regularly updated resulting in an inability to assess forest susceptibility. To address these shortcomings, we demonstrate how a geometric optical model (GOM) can be applied to Landsat-5 Thematic Mapper (TM) imagery (30 m spatial resolution) to estimate stand-level susceptibility to mountain pine beetle attack. Spectral mixture analysis was used to determine the proportion of sunlit canopy and background, and shadow of each Landsat pixel enabling per pixel estimates of attributes required for model inversion. Stand structural attributes were then derived from inversion of the geometric optical model and used as basis for susceptibility mapping. Mean stand density estimated by the geometric optical model was 2753 (standard deviation ± 308) stems per hectare and mean horizontal crown radius was 2.09 (standard deviation ± 0.11) metres. When compared to equivalent forest inventory attributes, model predictions of stems per hectare and crown radius were shown to be reasonably estimated using a Kruskal–Wallis ANOVA (p < 0.001). These predictions were then used to create a large area map that provided an assessment of the forest area susceptible to mountain pine beetle damage.     

An assessment of chlorophyll-a concentration spatio-temporal variation using Landsat satellite data,in a small tropical reservoir

Traditional approaches to monitoring aquatic systems are often limited by the need for data collection which often is time-consuming, expensive and non-continuous. The aim of the study was to map the spatio-temporal chlorophyll-a concentration changes in Malilangwe Reservoir, Zimbabwe as an indicator of phytoplankton biomass and trophic state when the reservoir was full (year 2000) and at its lowest capacity (year 2011), using readily available Landsat multispectral images. Medium-spatial resolution (30 m) Landsat multispectral Thematic Mapper TM 5 and ETM+ images for May to December 1999–2000 and 2010–2011 were used to derive chlorophyll-a concentrations. In situ measured chlorophyll-a and total suspended solids (TSS) concentrations for 2011 were employed to validate the Landsat chlorophyll-a and TSS estimates. The study results indicate that Landsat-derived chlorophyll-a and TSS estimates were comparable with field measurements. There was a considerable wet vs. dry season differences in total chlorophyll-a concentration, Secchi disc depth, TSS and turbidity within the reservoir. Using Permutational multivariate analyses of variance (PERMANOVA) analysis, there were significant differences (p < 0.0001) for chlorophyll-a concentration among sites, months and years whereas TSS was significant during the study months (p < 0.05). A strong positive significant correlation among both predicted TSS vs. chlorophyll-a and measured vs. predicted chlorophyll-a and TSS concentrations as well as an inverse relationship between reservoir chlorophyll-a concentrations and water level were found (p < 0.001 in all cases). In conclusion, total chlorophyll-a concentration in Malilangwe Reservoir was successfully derived from Landsat remote sensing data suggesting that the Landsat sensor is suitable for real-time monitoring over relatively short timescales and for small reservoirs. Satellite data can allow for surveying of chlorophyll-a concentration in aquatic ecosystems, thus, providing invaluable data in data scarce (limited on site ground measurements) environments.