Texture augmented detection of macrophyte species using decision trees

Image classification using multispectral sensors has shown good performance in detecting macrophytes at the species level. However, species level classification often does not utilize the texture information provided by high resolution images. This study investigated whether image texture provides useful vector(s) for the discrimination of monospecific stands of three floating macrophyte species in Quickbird imagery of the South Nation River. Semivariograms indicated that window sizes of 5 × 5 and 13 × 13 pixels were the most appropriate spatial scales for calculation of the grey level co-occurrence matrix and subsequent texture attributes from the multispectral and panchromatic bands. Of the 214 investigated vectors (13 Haralick texture attributes * 15 bands + 9 spectral bands + 10 transformations/indices), feature selection determined which combination of spectral and textural vectors had the greatest class separability based on the Mann–Whitney U-test and Jefferies–Matusita distance. While multispectral red and near infrared (NIR) performed satisfactorily, the addition of panchromatic-dissimilarity slightly improved class separability and the accuracy of a decision tree classifier (Kappa: red/NIR/panchromatic-dissimilarity – 93.2% versus red/NIR – 90.4%). Class separability improved by incorporating a second texture attribute, but resulted in a decrease in classification accuracy. The results suggest that incorporating image texture may be beneficial for separating stands with high spatial heterogeneity. However, the benefits may be limited and must be weighed against the increased complexity of the classifier.     

日光温室的纹理特征分层提取

针对高分辨率影像上日光温室的信息提取问题,该文提出了利用支持向量机、最近邻算法结合纹理特征在不同层上分别提取连片日光温室和独栋日光温室的方法。实验表明:纹理特征能提高分类精度,在大尺度的层上,分类精度提升幅度较大,但在小尺度的层上,分类精度提升幅度会比较小;并不是参与运算特征数越多,分类精度越高,多数情况下光谱+纹理组合的分类精度最高;提取连片日光温室的最优方案是支持向量机和光谱+形状+纹理(7像素×7像素),总精度为92.86%,Kappa系数为0.90,而提取独栋日光温室最优方案为SVM和光谱+纹理(11像素×11像素),总精度为88.39%,Kappa系数为0.86。     

Estimating soil salinity in Pingluo County of China using QuickBird data and soil reflectance spectra

Soil salinization is a worldwide environmental problem with severe economic and social consequences. In this paper, estimating the soil salinity of Pingluo County, China by a partial least squares regression (PLSR) predictive model was carried out using QuickBird data and soil reflectance spectra. At first, a relationship between the sensitive bands of soil salinity acquired from measured reflectance spectra and the spectral coverage of seven commonly used optical sensors was analyzed. Secondly, the potentiality of QuickBird data in estimating soil salinity by analyzing the correlations between the measured reflectance spectra and reflectance spectra derived from QuickBird data and analyzing the contributions of each band of QuickBird data to soil salinity estimation Finally, a PLSR predictive model of soil salinity was developed using reflectance spectra from QuickBird data and eight spectral indices derived from QuickBird data. The results indicated that the sensitive bands covered several bands of each optical sensor and these sensors can be used for soil salinity estimation. The result of estimation model showed that an accurate prediction of soil salinity can be made based on the PLSR method (R² = 0.992, RMSE = 0.195). The PLSR model's performance was better than that of the stepwise multiple regression (SMR) method. The results also indicated that using spectral indices such as intensity within spectral bands (Int1, Int2), soil salinity indices (SI1, SI2, SI3), the brightness index (BI), the normalized difference vegetation index (NDVI) and the ratio vegetation index (RVI) as independent model variables can help to increase the accuracy of soil salinity mapping. The NDVI and RVI can help to reduce the influences of vegetation cover and soil moisture on prediction accuracy. The method developed in this paper can be applied in other arid and semi-arid areas, such as western China.     

Assessing the accuracy of hyperspectral and multispectral satellite imagery for categorical and Quantitative mapping of salinity stress in sugarcane fields

This study evaluates the feasibility of hyperspectral and multispectral satellite imagery for categorical and quantitative mapping of salinity stress in sugarcane fields located in the southwest of Iran. For this purpose a Hyperion image acquired on September 2, 2010 and a Landsat7 ETM+ image acquired on September 7, 2010 were used as hyperspectral and multispectral satellite imagery. Field data including soil salinity in the sugarcane root zone was collected at 191 locations in 25 fields during September 2010. In the first section of the paper, based on the yield potential of sugarcane as influenced by different soil salinity levels provided by FAO, soil salinity was classified into three classes, low salinity (1.7–3.4 dS/m), moderate salinity (3.5–5.9 dS/m) and high salinity (6–9.5) by applying different classification methods including Support Vector Machine (SVM), Spectral Angle Mapper (SAM), Minimum Distance (MD) and Maximum Likelihood (ML) on Hyperion and Landsat images. In the second part of the paper the performance of nine vegetation indices (eight indices from literature and a new developed index in this study) extracted from Hyperion and Landsat data was evaluated for quantitative mapping of salinity stress. The experimental results indicated that for categorical classification of salinity stress, Landsat data resulted in a higher overall accuracy (OA) and Kappa coefficient (KC) than Hyperion, of which the MD classifier using all bands or PCA (1–5) as an input performed best with an overall accuracy and kappa coefficient of 84.84% and 0.77 respectively. Vice versa for the quantitative estimation of salinity stress, Hyperion outperformed Landsat. In this case, the salinity and water stress index (SWSI) has the best prediction of salinity stress with an R² of 0.68 and RMSE of 1.15 dS/m for Hyperion followed by Landsat data with an R² and RMSE of 0.56 and 1.75 dS/m respectively. It was concluded that categorical mapping of salinity stress is the best option for monitoring agricultural fields and for this purpose Landsat data are most suitable.     

黄河三角洲地区GF-3雷达数据与Sentinel-2多光谱数据湿地协同分类研究

黄河三角洲湿地的动态变化监测对湿地资源合理利用、开发保护具有重要意义。采用C波段全极化高分三号（GF-3）合成孔径雷达数据与欧洲空间局哨兵二号（Sentinel-2B）多光谱数据，分析了黄河三角洲湿地7类地物的光谱、指数、极化散射以及纹理等特征信息，分别基于最大似然法（maximum likelihood，ML）、决策树（decision tree，DT）、支持向量机（support vector machine，SVM）方法实现了有监督分类，评估了两者协同与单独应用于湿地地物分类与识别的能力，结果表明，两者协同分类时，其总体精度分别可达90.4%、95.4%、95.7%，均明显高于两者单独分类的结果，证明了GF-3雷达数据与多光谱数据在湿地协同分类方面的可靠性和应用潜力。     

Fusion of RADARSAT-2 and multispectral optical remote sensing data for LULC extraction in a tropical agricultural area

In this study, we investigated the performance of different fusion and classification techniques for land cover mapping in Hilir Perak, Peninsula Malaysia using RADAR and Landsat-8 images in a predominantly agricultural area. The fusion methods used are Brovey Transform, Wavelet Transform, Ehlers and Layer Stacking and their results classified into seven different land cover classes which include (1) pixel-based classifiers (spectral angle mapper (SAM), maximum likelihood (ML), support vector machine (SVM)) and (2) Object-based (rule-based and standard nearest neighbour (NN)) classifiers. The result shows that pixel-based classification achieved maximum accuracy of the optical data classification using SVM in Landsat-8 with 74.96% accuracy compared to SAM and ML. For multisource data classification, the highest overall accuracy recorded for layer stacking (SVM) was 79.78%, Ehlers fusion (SVM) with 45.57%, Brovey fusion (SVM) with 63.70% and Wavelet fusion (SVM) 61.16%. And for object-based classifiers, the overall classification accuracy is 95.35% for rule-based and 76.33% for NN classifier, respectively. Based on the analysis of their performances, object-based and the rule-based classifiers produced the best classification accuracy from the fused images.     

Enhancing the performance of regional land cover mapping

Different pixel-based, object-based and subpixel-based methods such as time-series analysis, decision-tree, and different supervised approaches have been proposed to conduct land use/cover classification. However, despite their proven advantages in small dataset tests, their performance is variable and less satisfactory while dealing with large datasets, particularly, for regional-scale mapping with high resolution data due to the complexity and diversity in landscapes and land cover patterns, and the unacceptably long processing time. The objective of this paper is to demonstrate the comparatively highest performance of an operational approach based on integration of multisource information ensuring high mapping accuracy in large areas with acceptable processing time. The information used includes phenologically contrasted multiseasonal and multispectral bands, vegetation index, land surface temperature, and topographic features. The performance of different conventional and machine learning classifiers namely Malahanobis Distance (MD), Maximum Likelihood (ML), Artificial Neural Networks (ANNs), Support Vector Machines (SVMs) and Random Forests (RFs) was compared using the same datasets in the same IDL (Interactive Data Language) environment. An Eastern Mediterranean area with complex landscape and steep climate gradients was selected to test and develop the operational approach. The results showed that SVMs and RFs classifiers produced most accurate mapping at local-scale (up to 96.85% in Overall Accuracy), but were very time-consuming in whole-scene classification (more than five days per scene) whereas ML fulfilled the task rapidly (about 10 min per scene) with satisfying accuracy (94.2–96.4%). Thus, the approach composed of integration of seasonally contrasted multisource data and sampling at subclass level followed by a ML classification is a suitable candidate to become an operational and effective regional land cover mapping method.     

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.     

Detection and mapping the spatial distribution of bracken fern weeds using the Landsat 8 OLI new generation sensor

Bracken fern is an invasive plant that presents serious environmental, ecological and economic problems around the world. An understanding of the spatial distribution of bracken fern weeds is therefore essential for providing appropriate management strategies at both local and regional scales. The aim of this study was to assess the utility of the freely available medium resolution Landsat 8 OLI sensor in the detection and mapping of bracken fern at the Cathedral Peak, South Africa. To achieve this objective, the results obtained from Landsat 8 OLI were compared with those derived using the costly, high spatial resolution WorldView-2 imagery. Since previous studies have already successfully mapped bracken fern using high spatial resolution WorldView-2 image, the comparison was done to investigate the magnitude of difference in accuracy between the two sensors in relation to their acquisition costs. To evaluate the performance of Landsat 8 OLI in discriminating bracken fern compared to that of Worldview-2, we tested the utility of (i) spectral bands; (ii) derived vegetation indices as well as (iii) the combination of spectral bands and vegetation indices based on discriminant analysis classification algorithm. After resampling the training and testing data and reclassifying several times (n = 100) based on the combined data sets, the overall accuracies for both Landsat 8 and WorldView-2 were tested for significant differences based on Mann-Whitney U test. The results showed that the integration of the spectral bands and derived vegetation indices yielded the best overall classification accuracy (80.08% and 87.80% for Landsat 8 OLI and WorldView-2 respectively). Additionally, the use of derived vegetation indices as a standalone data set produced the weakest overall accuracy results of 62.14% and 82.11% for both the Landsat 8 OLI and WorldView-2 images. There were significant differences {U (100) = 569.5, z = −10.8242, p < 0.01} between the classification accuracies derived based on Landsat OLI 8 and those derived using WorldView-2 sensor. Although there were significant differences between Landsat and WorldView-2 accuracies, the magnitude of variation (9%) between the two sensors was within an acceptable range. Therefore, the findings of this study demonstrated that the recently launched Landsat 8 OLI multispectral sensor provides valuable information that could aid in the long term continuous monitoring and formulation of effective bracken fern management with acceptable accuracies that are comparable to those obtained from the high resolution WorldView-2 commercial sensor.     

Texture-based classification of sub-Antarctic vegetation communities on Heard Island

This study was the first to use high-resolution IKONOS imagery to classify vegetation communities on sub-Antarctic Heard Island. We focused on the use of texture measures, in addition to standard multispectral information, to improve the classification of sub-Antarctic vegetation communities. Heard Island’s pristine and rapidly changing environment makes it a relevant and exciting location to study the regional effects of climate change. This study uses IKONOS imagery to provide automated, up-to-date, and non-invasive means to map vegetation as an important indicator for environmental change. Three classification techniques were compared: multispectral classification, texture based classification, and a combination of both. Texture features were calculated using the Grey Level Co-occurrence Matrix (GLCM). We investigated the effect of the texture window size on classification accuracy. The combined approach produced a higher accuracy than using multispectral bands alone. It was also found that the selection of GLCM texture features is critical. The highest accuracy (85%) was produced using all original spectral bands and three uncorrelated texture features. Incorporating texture improved classification accuracy by 6%.     

Assessing the potential of multi-seasonal WorldView-2 imagery for mapping West African agroforestry tree species

High resolution satellite systems enable efficient and detailed mapping of tree cover, with high potential to support both natural resource monitoring and ecological research. This study investigates the capability of multi-seasonal WorldView-2 imagery to map five dominant tree species at the individual tree crown level in a parkland landscape in central Burkina Faso. The Random Forest algorithm is used for object based tree species classification and for assessing the relative importance of WorldView-2 predictors. The classification accuracies from using wet season, dry season and multi-seasonal datasets are compared to gain insights about the optimal timing for image acquisition. The multi-seasonal dataset produced the most accurate classifications, with an overall accuracy (OA) of 83.4%. For classifications based on single date imagery, the dry season (OA = 78.4%) proved to be more suitable than the wet season (OA = 68.1%). The predictors that contributed most to the classification success were based on the red edge band and visible wavelengths, in particular green and yellow. It was therefore concluded that WorldView-2, with its unique band configuration, represents a suitable data source for tree species mapping in West African parklands. These results are particularly promising when considering the recently launched WorldView-3, which provides data both at higher spatial and spectral resolution, including shortwave infrared bands.     

Multi range spectral feature fitting for hyperspectral imagery in extracting oilseed rape planting area

Spectral feature fitting (SFF) is a commonly used strategy for hyperspectral imagery analysis to discriminate ground targets. Compared to other image analysis techniques, SFF does not secure higher accuracy in extracting image information in all circumstances. Multi range spectral feature fitting (MRSFF) from ENVI software allows user to focus on those interesting spectral features to yield better performance. Thus spectral wavelength ranges and their corresponding weights must be determined. The purpose of this article is to demonstrate the performance of MRSFF in oilseed rape planting area extraction. A practical method for defining the weighted values, the variance coefficient weight method, was proposed to set up criterion. Oilseed rape field canopy spectra from the whole growth stage were collected prior to investigating its phenological varieties; oilseed rape endmember spectra were extracted from the Hyperion image as identifying samples to be used in analyzing the oilseed rape field. Wavelength range divisions were determined by the difference between field-measured spectra and image spectra, and image spectral variance coefficient weights for each wavelength range were calculated corresponding to field-measured spectra from the closest date. By using MRSFF, wavelength ranges were classified to characterize the target's spectral features without compromising spectral profile's entirety. The analysis was substantially successful in extracting oilseed rape planting areas (RMSE ≤ 0.06), and the RMSE histogram indicated a superior result compared to a conventional SFF. Accuracy assessment was based on the mapping result compared with spectral angle mapping (SAM) and the normalized difference vegetation index (NDVI). The MRSFF yielded a robust, convincible result and, therefore, may further the use of hyperspectral imagery in precision agriculture.     

Evaluation of forest fire on Madeira Island using Sentinel-2A MSI imagery

A forest fire started on August 8th, 2016 in several places on Madeira Island causing damage and casualties. As of August 10th the local media had reported the death of three people, over 200 people injured, over 950 habitants evacuated, and 50 houses damaged. This study presents the preliminary results of the assessment of several spectral indices to evaluate the burn severity of Madeira fires during August 2016. These spectral indices were calculated using the new European satellite Sentinel-2A launched in June 2015. The study confirmed the advantages of several spectral indices such as Normalized Difference Vegetation Index (NDVI), Green Normalized Difference Vegetation Index (GNDVI), Normalized Burn Ratio (NBR) and Normalized Difference Vegetation Index (NDVIreXn) using red-edge spectral bands to assess the post-fire conditions. Results showed high correlation between NDVI, GNDVI, NBR and NDVIre1n spectral indices and the analysis performed by Copernicus Emergency Management Service (EMSR175), considered as the reference truth. Regarding the red-edge spectral indices, the NDVIre1n (using band B5, 705 nm) presented better results compared with B6 (740 nm) and B7 (783 nm) bands. These preliminary results allow us to assume that Sentinel-2 will be a valuable tool for post-fire monitoring. In the future, the two twin Sentinel-2 satellites will offer global coverage of the Madeira Archipelago every five days, therefore allowing the simultaneous study of the evolution of the burnt area and reforestation information with high spatial (up to 10 m) and temporal resolution (5 days).     

基于CART决策树的高分二号洞庭湖区湿地分类方法

高分卫星遥感湿地分类的关键在于解决“同物异谱、异物同谱”难题。本文将当前应用前景广泛的亚米级国产高分二号（GF-2）影像和CART决策树面向对象分类算法相结合,以湖南沅江为例进行洞庭湖区典型湿地的分类提取工作,选取包括光谱信息、几何特征、地形特征和纹理特征等多维对象特征对分类器进行训练,构建了多维特征湿地分类方法。试验区总体分类精度优于传统方法,可为基于GF-2影像的洞庭湖湿地分类提供技术参考。     

Improved coastal wetland mapping using very-high 2-meter spatial resolution imagery

Accurate wetland maps are a fundamental requirement for land use management and for wetland restoration planning. Several wetland map products are available today; most of them based on remote sensing images, but their different data sources and mapping methods lead to substantially different estimations of wetland location and extent. We used two very high-resolution (2 m) WorldView-2 satellite images and one (30 m) Landsat 8 Operational Land Imager (OLI) image to assess wetland coverage in two coastal areas of Tampa Bay (Florida): Fort De Soto State Park and Weedon Island Preserve. An initial unsupervised classification derived from WorldView-2 was more accurate at identifying wetlands based on ground truth data collected in the field than the classification derived from Landsat 8 OLI (82% vs. 46% accuracy). The WorldView-2 data was then used to define the parameters of a simple and efficient decision tree with four nodes for a more exacting classification. The criteria for the decision tree were derived by extracting radiance spectra at 1500 separate pixels from the WorldView-2 data within field-validated regions. Results for both study areas showed high accuracy in both wetland (82% at Fort De Soto State Park, and 94% at Weedon Island Preserve) and non-wetland vegetation classes (90% and 83%, respectively). Historical, published land-use maps overestimate wetland surface cover by factors of 2–10 in the study areas. The proposed methods improve speed and efficiency of wetland map production, allow semi-annual monitoring through repeat satellite passes, and improve the accuracy and precision with which wetlands are identified.     

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.     

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.     

Remotely sensed characterization of forest fuel types by using satellite ASTER data

The characterization of fuel types is very important for computing spatial fire hazard and risk and simulating fire growth and intensity across a landscape. However, due to the complex nature of fuel characteristic a fuel map is considered one of the most difficult thematic layers to build up. The advent of sensors with increased spatial resolution may improve the accuracy and reduce the cost of fuels mapping. The objective of this research is to evaluate the accuracy and utility of imagery from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite imagery. In order to ascertain how well ASTER data can provide an exhaustive classification of fuel properties a sample area characterized by mixed vegetation covers was analysed. The selected sample areas has an extension at around 60 km² and is located inside the Sila plateau in the Calabria Region (South of Italy). Fieldwork fuel type recognitions, performed before, after and during the acquisition of remote sensing ASTER data, were used as ground-truth dataset to assess the results obtained for the considered test area. The method comprised the following three steps: (I) adaptation of Prometheus fuel types for obtaining a standardization system useful for remotely sensed classification of fuel types and properties in the considered Mediterranean ecosystems; (II) model construction for the spectral characterization and mapping of fuel types based on a maximum likelihood (ML) classification algorithm; (III) accuracy assessment for the performance evaluation based on the comparison of ASTER-based results with ground-truth. Results from our analysis showed that the use ASTER data provided a valuable characterization and mapping of fuel types being that the achieved classification accuracy was higher than 90%.     

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.     

Understanding the optical responses of leaf nitrogen in Mediterranean Holm oak (Quercus ilex) using field spectroscopy

The direct estimation of nitrogen (N) in fresh vegetation is challenging due to its weak influence on leaf reflectance and the overlaps with absorption features of other compounds. Different empirical models relate in this work leaf nitrogen concentration ([N]_Leaf) on Holm oak to leaf reflectance as well as derived spectral indices such as normalized difference indices (NDIs), the three bands indices (TBIs) and indices previously used to predict leaf N and chlorophyll. The models were calibrated and assessed their accuracy, robustness and the strength of relationship when other biochemicals were considered. Red edge was the spectral region most strongly correlated with [N]_Leaf, whereas most of the published spectral indexes did not provide accurate estimations. NDIs and TBIs based models could achieve robust and acceptable accuracies (TBI_{1310,1720,730}: R² = 0.76, [0.64,0.86]; RMSE (%) = 9.36, [7.04,12.83]). These models sometimes included indices with bands close to absorption features of N bonds or nitrogenous compounds, but also of other biochemicals. Models were independently and inter-annually validated using the bootstrap method, which allowed discarding those models non-robust across different years. Partial correlation analysis revealed that spectral estimators did not strongly respond to [N]_Leaf but to other leaf variables such as chlorophyll and water, even if bands close to absorption features of N bonds or compounds were present in the models.