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.     

Stochastic gradient boosting classification trees for forest fuel types mapping through airborne laser scanning and IRS LISS-III imagery

This paper presents an application of Airborne Laser Scanning (ALS) data in conjunction with an IRS LISS-III image for mapping forest fuel types. For two study areas of 165 km² and 487 km² in Sicily (Italy), 16,761 plots of size 30-m × 30-m were distributed using a tessellation-based stratified sampling scheme. ALS metrics and spectral signatures from IRS extracted for each plot were used as predictors to classify forest fuel types observed and identified by photointerpretation and fieldwork. Following use of traditional parametric methods that produced unsatisfactory results, three non-parametric classification approaches were tested: (i) classification and regression tree (CART), (ii) the CART bagging method called Random Forests, and (iii) the CART bagging/boosting stochastic gradient boosting (SGB) approach. This contribution summarizes previous experiences using ALS data for estimating forest variables useful for fire management in general and for fuel type mapping, in particular. It summarizes characteristics of classification and regression trees, presents the pre-processing operation, the classification algorithms, and the achieved results. The results demonstrated superiority of the SGB method with overall accuracy of 84%. The most relevant ALS metric was canopy cover, defined as the percent of non-ground returns. Other relevant metrics included the spectral information from IRS and several other ALS metrics such as percentiles of the height distribution, the mean height of all returns, and the number of returns.     

Modeling biophysical properties of broad-leaved stands in the hyrcanian forests of Iran using fused airborne laser scanner data and ultraCam-D images

Inventories of mixed broad-leaved forests of Iran mainly rely on terrestrial measurements. Due to rapid changes and disturbances and great complexity of the silvicultural systems of these multilayer forests, frequent repetition of conventional ground-based plot surveys is often cost prohibitive. Airborne laser scanning (ALS) and multispectral data offer an alternative or supplement to conventional inventories in the Hyrcanian forests of Iran. In this study, the capability of a combination of ALS and UltraCam-D data to model stand volume, tree density, and basal area using random forest (RF) algorithm was evaluated. Systematic sampling was applied to collect field plot data on a 150 m × 200 m sampling grid within a 1100 ha study area located at 36°38′- 36°42′N and 54°24′–54°25′E. A total of 308 circular plots (0.1 ha) were measured for calculation of stand volume, tree density, and basal area per hectare. For each plot, a set of variables was extracted from both ALS and multispectral data. The RF algorithm was used for modeling of the biophysical properties using ALS and UltraCam-D data separately and combined. The results showed that combining the ALS data and UltraCam-D images provided a slight increase in prediction accuracy compared to separate modeling. The RMSE as percentage of the mean, the mean difference between observed and predicted values, and standard deviation of the differences using a combination of ALS data and UltraCam-D images in an independent validation at 0.1-ha plot level were 31.7%, 1.1%, and 84 m³ ha⁻¹ for stand volume; 27.2%, 0.86%, and 6.5 m² ha⁻¹ for basal area, and 35.8%, −4.6%, and 77.9 n ha⁻¹ for tree density, respectively. Based on the results, we conclude that fusion of ALS and UltraCam-D data may be useful for modeling of stand volume, basal area, and tree density and thus gain insights into structural characteristics in the complex Hyrcanian forests.     

Modeling Mediterranean forest structure using airborne laser scanning data

The conservation of biological diversity is recognized as a fundamental component of sustainable development, and forests contribute greatly to its preservation. Structural complexity increases the potential biological diversity of a forest by creating multiple niches that can host a wide variety of species. To facilitate greater understanding of the contributions of forest structure to forest biological diversity, we modeled relationships between 14 forest structure variables and airborne laser scanning (ALS) data for two Italian study areas representing two common Mediterranean forests, conifer plantations and coppice oaks subjected to irregular intervals of unplanned and non-standard silvicultural interventions. The objectives were twofold: (i) to compare model prediction accuracies when using two types of ALS metrics, echo-based metrics and canopy height model (CHM)-based metrics, and (ii) to construct inferences in the form of confidence intervals for large area structural complexity parameters.Our results showed that the effects of the two study areas on accuracies were greater than the effects of the two types of ALS metrics. In particular, accuracies were less for the more complex study area in terms of species composition and forest structure. However, accuracies achieved using the echo-based metrics were only slightly greater than when using the CHM-based metrics, thus demonstrating that both options yield reliable and comparable results. Accuracies were greatest for dominant height (Hd) (R² = 0.91; RMSE% = 8.2%) and mean height weighted by basal area (R² = 0.83; RMSE% = 10.5%) when using the echo-based metrics, 99th percentile of the echo height distribution and interquantile distance. For the forested area, the generalized regression (GREG) estimate of mean Hd was similar to the simple random sampling (SRS) estimate, 15.5 m for GREG and 16.2 m SRS. Further, the GREG estimator with standard error of 0.10 m was considerable more precise than the SRS estimator with standard error of 0.69 m.     

Application of 3D triangulations of airborne laser scanning data to estimate boreal forest leaf area index

We propose 3D triangulations of airborne Laser Scanning (ALS) point clouds as a new approach to derive 3D canopy structures and to estimate forest canopy effective LAI (LAI_e). Computational geometry and topological connectivity were employed to filter the triangulations to yield a quasi-optimal relationship with the field measured LAI_e. The optimal filtering parameters were predicted based on ALS height metrics, emulating the production of maps of LAI_e and canopy volume for large areas. The LAI_e from triangulations was validated with field measured LAI_e and compared with a reference LAI_e calculated from ALS data using logarithmic model based on Beer’s law. Canopy transmittance was estimated using All Echo Cover Index (ACI), and the mean projection of unit foliage area (β) was obtained using no-intercept regression with field measured LAI_e. We investigated the influence species and season on the triangulated LAI_e and demonstrated the relationship between triangulated LAI_e and canopy volume. Our data is from 115 forest plots located at the southern boreal forest area in Finland and for each plot three different ALS datasets were available to apply the triangulations. The triangulation approach was found applicable for both leaf-on and leaf-off datasets after initial calibration. Results showed the Root Mean Square Errors (RMSEs) between LAI_e from triangulations and field measured values agreed the most using the highest pulse density data (RMSE = 0.63, the coefficient of determination (R²) = 0.53). Yet, the LAI_e calculated using ACI-index agreed better with the field measured LAI_e (RMSE = 0.53 and R² = 0.70). The best models to predict the optimal alpha value contained the ACI-index, which indicates that within-crown transmittance is accounted by the triangulation approach. The cover indices may be recommended for retrieving LAI_e only, but for applications which require more sophisticated information on canopy shape and volume, such as radiative transfer models, the triangulation approach may be preferred.     

Effect of flying altitude,scanning angle and scanning mode on the accuracy of ALS based forest inventory

Airborne laser scanning (ALS) is a widely used technology in the mapping of environment and forests. Data acquisition costs and the accuracy of the forest inventory are closely dependent on some extrinsic parameters of the ALS survey. These parameters have been assessed in numerous studies about a decade ago, but since then ALS devices have developed and it is possible that previous findings do not hold true with newer technology. That is why, the effect of flying altitudes (2000, 2500 or 3000 m), scanning angles (±15° and ±20° off nadir) and scanning modes (single- and multiple pulses in air) with the area-based approach using a Leica ALS70HA-laser scanner was studied here. The study was conducted in a managed pine-dominated forest area in Finland, where eight separate discrete-return ALS data were acquired. The comparison of datasets was based on the bootstrap approach with 5-fold cross validation. Results indicated that the narrower scanning angle (±15° i.e. 30°) led to slightly more accurate estimates of plot volume (RMSE%: 21–24 vs. 22.5–25) and mean height (RMSE%: 8.5–11 vs. 9–12). We also tested the use case where the models are constructed using one data and then applied to other data gathered with different parameters. The most accurate models were identified using the bootstrap approach and applied to different datasets with and without refitting. The bias increased without refitting the models (bias%: volume 0 ± 10, mean height 0 ± 3), but in most cases the results did not differ much in terms of RMSE%. This confirms previous observations that models should only be used for datasets collected under similar data acquisition conditions. We also calculated the proportions of echoes as a function of height for different echo categories. This indicated that the accuracy of the inventory is affected more by the height distribution than the proportions of echo categories.     

PTrees: A point-based approach to forest tree extraction from lidar data

This paper introduces PTrees, a multi-scale dynamic point cloud segmentation dedicated to forest tree extraction from lidar point clouds. The method process the point data using the raw elevation values (Z) and compute height (H = Z − ground elevation) during post-processing using an innovative procedure allowing to preserve the geometry of crown points. Multiple segmentations are done at different scales. Segmentation criteria are then applied to dynamically select the best set of apices from the tree segments extracted at the various scales. The selected set of apices is then used to generate a final segmentation. PTrees has been tested in 3 different forest types, allowing to detect 82% of the trees with under 10% of false detection rate. Future development will integrate crown profile estimation during the segmentation process in order to both maximize the detection of suppressed trees and minimize false detections.     

Scaling estimates of vegetation structure in Amazonian tropical forests using multi-angle MODIS observations

Detailed knowledge of vegetation structure is required for accurate modelling of terrestrial ecosystems, but direct measurements of the three dimensional distribution of canopy elements, for instance from LiDAR, are not widely available. We investigate the potential for modelling vegetation roughness, a key parameter for climatological models, from directional scattering of visible and near-infrared (NIR) reflectance acquired from NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS). We compare our estimates across different tropical forest types to independent measures obtained from: (1) airborne laser scanning (ALS), (2) spaceborne Geoscience Laser Altimeter System (GLAS)/ICESat, and (3) the spaceborne SeaWinds/QSCAT. Our results showed linear correlation between MODIS-derived anisotropy to ALS-derived entropy (r² = 0.54, RMSE = 0.11), even in high biomass regions. Significant relationships were also obtained between MODIS-derived anisotropy and GLAS-derived entropy (0.52 ≤ r² ≤ 0.61; p < 0.05), with similar slopes and offsets found throughout the season, and RMSE between 0.26 and 0.30 (units of entropy). The relationships between the MODIS-derived anisotropy and backscattering measurements (σ⁰) from SeaWinds/QuikSCAT presented an r² of 0.59 and a RMSE of 0.11. We conclude that multi-angular MODIS observations are suitable to extrapolate measures of canopy entropy across different forest types, providing additional estimates of vegetation structure in the Amazon.     

Characterization of aboveground biomass in an unmanaged boreal forest using Landsat temporal segmentation metrics

Canada is dominated by forested ecosystems which are subject to various inventory and management practices, with more northern boreal forests subject to neither. Our objectives were to measure the capacity of temporal trajectory metrics for estimating selected forest attributes in a northern Canadian boreal forest context using Landsat imagery and investigate the importance of different types of temporal trajectory metrics. Results indicated that Wetness was the best Tasseled Cap (TC) component for aboveground biomass estimation (R² = 50%, RMSE% = 56%), and the combination of simple and complex metrics from all TC components produced the highest R² (62%) and lowest RMSE% (49%). Using a similar combination of variables, other forest attributes were estimated equally reliably with lower RMSE% values. The most important temporal trajectory metrics were simple and described TC component values at each point of change in the temporal trajectory, however the most important variables overall were environmental variables.     

Timber production assessment of a plantation forest: An integrated framework with field-based inventory,multi-source remote sensing data and forest management history

Timber production is the purpose for managing plantation forests, and its spatial and quantitative information is critical for advising management strategies. Previous studies have focused on growing stock volume (GSV), which represents the current potential of timber production, yet few studies have investigated historical process-harvested timber. This resulted in a gap in a synthetical ecosystem service assessment of timber production. In this paper, we established a Management Process–based Timber production (MPT) framework to integrate the current GSV and the harvested timber derived from historical logging regimes, trying to synthetically assess timber production for a historical period. In the MPT framework, age-class and current GSV determine the times of historical thinning and the corresponding harvested timber, by using a “space-for-time” substitution. The total timber production can be estimated by the historical harvested timber in each thinning and the current GSV. To test this MPT framework, an empirical study on a larch plantation (LP) with area of 43,946 ha was conducted in North China for a period from 1962 to 2010. Field-based inventory data was integrated with ALOS PALSAR (Advanced Land-Observing Satellite Phased Array L-band Synthetic Aperture Radar) and Landsat-8 OLI (Operational Land Imager) data for estimating the age-class and current GSV of LP. The random forest model with PALSAR backscatter intensity channels and OLI bands as input predictive variables yielded an accuracy of 67.9% with a Kappa coefficient of 0.59 for age-class classification. The regression model using PALSAR data produced a root mean square error (RMSE) of 36.5 m³ ha⁻¹. The total timber production of LP was estimated to be 7.27 × 10⁶ m³, with 4.87 × 10⁶ m³ in current GSV and 2.40 × 10⁶ m³ in harvested timber through historical thinning. The historical process-harvested timber accounts to 33.0% of the total timber production, which component has been neglected in the assessments for current status of plantation forests. Synthetically considering the RMSE for predictive GSV and misclassification of age-class, the error in timber production were supposed to range from −55.2 to 56.3 m³ ha⁻¹. The MPT framework can be used to assess timber production of other tree species at a larger spatial scale, providing crucial information for a better understanding of forest ecosystem service.     

Assessing the utility WorldView-2 imagery for tree species mapping in South African subtropical humid forest and the conservation implications: Dukuduku forest patch as case study

Indigenous forest biome in South Africa is highly fragmented into patches of various sizes (most patches < 1 km²). The utilization of timber and non-timber resources by poor rural communities living around protected forest patches produce subtle changes in the forest canopy which can be hardly detected on a timely manner using traditional field surveys. The aims of this study were to assess: (i) the utility of very high resolution (VHR) remote sensing imagery (WorldView-2, 0.5–2 m spatial resolution) for mapping tree species and canopy gaps in one of the protected subtropical coastal forests in South Africa (the Dukuduku forest patch (ca.3200 ha) located in the province of KwaZulu-Natal) and (ii) the implications of the map products to forest conservation. Three dominant canopy tree species namely, Albizia adianthifolia, Strychnos spp. and Acacia spp., and canopy gap types including bushes (grass/shrubby), bare soil and burnt patches were accurately mapped (overall accuracy = 89.3 ± 2.1%) using WorldView-2 image and support vector machine classifier. The maps revealed subtle forest disturbances such as bush encroachment and edge effects resulting from forest fragmentation by roads and a power-line. In two stakeholders’ workshops organised to assess the implications of the map products to conservation, participants generally agreed amongst others implications that the VHR maps provide valuable information that could be used for implementing and monitoring the effects of rehabilitation measures. The use of VHR imagery is recommended for timely inventorying and monitoring of the small and fragile patches of subtropical forests in Southern Africa.     

An integrated remote sensing and GIS approach for monitoring areas affected by selective logging: A case study in northern Mato Grosso,Brazilian Amazon

Forest cover disturbances due to processes such as logging and forest fires are a widespread issue especially in the tropics, and have heavily affected forest biomass and functioning in the Brazilian Amazon in the past decades. Satellite remote sensing has played a key role for assessing logging activities in this region; however, there are still remaining challenges regarding the quantification and monitoring of these processes affecting forested lands. In this study, we propose a new method for monitoring areas affected by selective logging in one of the hotspots of Mato Grosso state in the Brazilian Amazon, based on a combination of object-based and pixel-based classification approaches applied on remote sensing data. Logging intensity and changes over time are assessed within grid cells of 300 m × 300 m spatial resolution. Our method encompassed three main steps: (1) mapping forest/non-forest areas through an object-based classification approach applied to a temporal series of Landsat images during the period 2000–2015, (2) mapping yearly logging activities from soil fraction images on the same Landsat data series, and (3) integrating information from previous steps within a regular grid-cell of 300 m × 300 m in order to monitor disturbance intensities over this 15-years period. The overall accuracy of the baseline forest/non-forest mask (year 2000) and of the undisturbed vs disturbed forest (for selected years) were 93% and 84% respectively. Our results indicate that annual forest disturbance rates, mainly due to logging activities, were higher than annual deforestation rates during the whole period of study. The deforested areas correspond to circa 25% of the areas affected by forest disturbances. Deforestation rates were highest from 2001 to 2005 and then decreased considerably after 2006. In contrast, the annual forest disturbance rates show high temporal variability with a slow decrease over the 15-year period, resulting in a significant increase of the ratio between disturbed and deforested areas. Although the majority of the areas, which have been affected by selective logging during the period 2000–2014, were not deforested by 2015, more than 70% of the deforested areas in 2015 had been at least once identified as disturbed forest during that period.     

Object-based class modelling for multi-scale riparian forest habitat mapping

Object-based class modelling allows for mapping complex, hierarchical habitat systems. The riparian zone, including forests, represents such a complex ecosystem. Forests within riparian zones are biologically high productive and characterized by a rich biodiversity; thus considered of high community interest with an imperative to be protected and regularly monitored. Satellite earth observation (EO) provides tools for capturing the current state of forest habitats such as forest composition including intermixture of non-native tree species. Here we present a semi-automated object based image analysis (OBIA) approach for the mapping of riparian forests by applying class modelling of habitats based on the European Nature Information System (EUNIS) habitat classifications and the European Habitats Directive (HabDir) Annex 1. A very high resolution (VHR) WorldView-2 satellite image provided the required spatial and spectral details for a multi-scale image segmentation and rule-base composition to generate a six-level hierarchical representation of riparian forest habitats. Thereby habitats were hierarchically represented within an image object hierarchy as forest stands, stands of homogenous tree species and single trees represented by sunlit tree crowns. 522 EUNIS level 3 (EUNIS-3) habitat patches with a mean patch size (MPS) of 12,349.64 m² were modelled from 938 forest stand patches (MPS = 6868.20 m²) and 43,742 tree stand patches (MPS = 140.79 m²). The delineation quality of the modelled EUNIS-3 habitats (focal level) was quantitatively assessed to an expert-based visual interpretation showing a mean deviation of 11.71%.     

Testing the discrimination and detection limits of WorldView-2 imagery on a challenging invasive plant target

Invasive plants pose significant threats to biodiversity and ecosystem function globally, leading to costly monitoring and management effort. While remote sensing promises cost-effective, robust and repeatable monitoring tools to support intervention, it has been largely restricted to airborne platforms that have higher spatial and spectral resolutions, but which lack the coverage and versatility of satellite-based platforms. This study tests the ability of the WorldView-2 (WV2) eight-band satellite sensor for detecting the invasive shrub mesquite (Prosopis spp.) in the north-west Pilbara region of Australia. Detectability was challenged by the target taxa being largely defoliated by a leaf-tying biological control agent (Gelechiidae: Evippe sp. #1) and the presence of other shrubs and trees. Variable importance in the projection (VIP) scores identified bands offering greatest capacity for discrimination were those covering the near-infrared, red, and red-edge wavelengths. Wavelengths between 400 nm and 630 nm (coastal blue, blue, green, yellow) were not useful for species level discrimination in this case. Classification accuracy was tested on three band sets (simulated standard multispectral, all bands, and bands with VIP scores ≥1). Overall accuracies were comparable amongst all band-sets (Kappa = 0.71–0.77). However, mesquite omission rates were unacceptably high (21.3%) when using all eight bands relative to the simulated standard multispectral band-set (9.5%) and the band-set informed by VIP scores (11.9%). An incremental cover evaluation on the latter identified most omissions to be for objects <16 m². Mesquite omissions reduced to 2.6% and overall accuracy significantly improved (Kappa = 0.88) when these objects were left out of the confusion matrix calculations. Very high mapping accuracy of objects >16 m² allows application for mapping mesquite shrubs and coalesced stands, the former not previously possible, even with 3 m resolution hyperspectral imagery. WV2 imagery offers excellent portability potential for detecting other species where spectral/spatial resolution or coverage has been an impediment. New generation satellite sensors are removing barriers previously preventing widespread adoption of remote sensing technologies in natural resource management.     

Mapping and estimating land change between 2001 and 2013 in a heterogeneous landscape in West Africa: Loss of forestlands and capacity building opportunities

In West Africa, accurate classification of land cover and land change remains a big challenge due to the patchy and heterogeneous nature of the landscape. Limited data availability, human resources and technical capacities, further exacerbate the challenge. The result is a region that is among the more understudied areas in the world, which in turn has resulted in a lack of appropriate information required for sustainable natural resources management. The objective of this paper is to explore open source software and easy-to-implement approaches to mapping and estimation of land change that are transferrable to local institutions to increase capacity in the region, and to provide updated information on the regional land surface dynamics. To achieve these objectives, stable land cover and land change between 2001 and 2013 in the Kara River Basin in Togo and Benin were mapped by direct multitemporal classification of Landsat data by parameterization and evaluation of two machine-learning algorithms. Areas of land cover and change were estimated by application of an unbiased estimator to sample data following international guidelines. A prerequisite for all tools and methods was implementation in an open source environment, and adherence to international guidelines for reporting land surface activities. Findings include a recommendation of the Random Forests algorithm as implemented in Orfeo Toolbox, and a stratified estimation protocol − all executed in the QGIS graphical use interface. It was found that despite an estimated reforestation of 10,0727 ± 3480 ha (95% confidence interval), the combined rate of forest and savannah loss amounted to 56,271 ± 9405 ha (representing a 16% loss of the forestlands present in 2001), resulting in a rather sharp net loss of forestlands in the study area. These dynamics had not been estimated prior to this study, and the results will provide useful information for decision making pertaining to natural resources management, land management planning, and the implementation of the United Nations Collaborative Programme on Reducing Emissions from Deforestation and Forest Degradation in Developing Countries (UN-REDD).     

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.     

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.     

Estimating temporal changes in soil carbon stocks at ecoregional scale in Madagascar using remote-sensing

Soil organic carbon (SOC) plays an important role in climate change regulation notably through release of CO2 following land use change such a deforestation, but data on stock change levels are lacking. This study aims to empirically assess SOC stocks change between 1991 and 2011 at the landscape scale using easy-to-access spatially-explicit environmental factors. The study area was located in southeast Madagascar, in a region that exhibits very high rate of deforestation and which is characterized by both humid and dry climates. We estimated SOC stock on 0.1 ha plots for 95 different locations in a 43,000 ha reference area covering both dry and humid conditions and representing different land cover including natural forest, cropland, pasture and fallows. We used the Random Forest algorithm to find out the environmental factors explaining the spatial distribution of SOC. We then predicted SOC stocks for two soil layers at 30 cm and 100 cm over a wider area of 395,000 ha. By changing the soil and vegetation indices derived from remote sensing images we were able to produce SOC maps for 1991 and 2011. Those estimates and their related uncertainties where combined in a post-processing step to map estimates of significant SOC variations and we finally compared the SOC change map with published deforestation maps. Results show that the geologic variables, precipitation, temperature, and soil-vegetation status were strong predictors of SOC distribution at regional scale. We estimated an average net loss of 10.7% and 5.2% for the 30 cm and the 100 cm layers respectively for deforested areas in the humid area. Our results also suggest that these losses occur within the first five years following deforestation. No significant variations were observed for the dry region. This study provides new solutions and knowledge for a better integration of soil threats and opportunities in land management policies.     

Detecting Sirex noctilio grey-attacked and lightning-struck pine trees using airborne hyperspectral data,random forest and support vector machines classifiers

The visual progression of sirex (Sirex noctilio) infestation symptoms has been categorized into three distinct infestation phases, namely the green, red and grey stages. The grey stage is the final stage which leads to almost complete defoliation resulting in dead standing trees or snags. Dead standing pine trees however, could also be due to the lightning damage. Hence, the objective of the present study was to distinguish amongst healthy, sirex grey-attacked and lightning-damaged pine trees using AISA Eagle hyperspectral data, random forest (RF) and support vector machines (SVM) classifiers. Our study also presents an opportunity to look at the possibility of separating amongst the previously mentioned pine trees damage classes and other landscape classes on the study area. The results of the present study revealed the robustness of the two machine learning classifiers with an overall accuracy of 74.50% (total disagreement = 26%) for RF and 73.50% (total disagreement = 27%) for SVM using all the remaining AISA Eagle spectral bands after removing the noisy ones. When the most useful spectral bands as measured by RF were exploited, the overall accuracy was considerably improved; 78% (total disagreement = 22%) for RF and 76.50% (total disagreement = 24%) for SVM. There was no significant difference between the performances of the two classifiers as demonstrated by the results of McNemar’s test (chi-squared; χ² = 0.14, and 0.03 when all the remaining ASIA Eagle wavebands, after removing the noisy ones and the most important wavebands were used, respectively). This study concludes that AISA Eagle data classified using RF and SVM algorithms provide relatively accurate information that is important to the forest industry for making informed decision regarding pine plantations health protocols.