Mapping tree carbon with airborne remote sensing

Forests are a major component of the global carbon cycle, and accurate estimation of forest carbon stocks and fluxes is important in the context of anthropogenic global change. Airborne laser scanning (ALS) data sets are increasingly recognized as outstanding data sources for high-fidelity mapping of carbon stocks at regional scales.We develop a tree-centric approach to carbon mapping, based on identifying individual tree crowns (ITCs) and species from airborne remote sensing data, from which individual tree carbon stocks are calculated. We identify ITCs from the laser scanning point cloud using a region-growing algorithm and identifying species from airborne hyperspectral data by machine learning. For each detected tree, we predict stem diameter from its height and crown-width estimate. From that point on, we use well-established approaches developed for field-based inventories: above-ground biomasses of trees are estimated using published allometries and summed within plots to estimate carbon density.We show this approach is highly reliable: tests in the Italian Alps demonstrated a close relationship between field- and ALS-based estimates of carbon stocks (r2 = 0·98). Small trees are invisible from the air, and a correction factor is required to accommodate this effect.An advantage of the tree-centric approach over existing area-based methods is that it can produce maps at any scale and is fundamentally based on field-based inventory methods, making it intuitive and transparent. Airborne laser scanning, hyperspectral sensing and computational power are all advancing rapidly, making it increasingly feasible to use ITC approaches for effective mapping of forest carbon density also inside wider carbon mapping programs like REDD++.We thank Dr L. Frizzera for help with field-data collection. ALS data acquisition was supported by the European Commission (Alpine Space 2-3-2-FR NEWFOR). MD was supported by Trees4Future (European Union FP7 284181) and a NERC grant NE/K016377/1. DAC was also supported by a grant from BBSRC and DEFRA to study ash dieback.This is the final version of the article. It first appeared from Wiley via https://doi.org/10.1111/2041-210X.1257

Estimating grassland vegetation cover with remote sensing: a comparison between Landsat-8, Sentinel-2 and PlanetScope imagery

Grassland fractional vegetation cover (FVC) accurate mapping on a large scale is crucial, since degraded grasslands contribute less to provisioning services, carbon storage, water purification, erosion control and biodiversity conservation. The spatial and temporal resolution of Sentinel-2 (S2) and PlanetScope (PS) data has never been explored for grassland FVC estimation so far and will enable researchers and agencies to quantify and map timelier and more precisely grassland processes. In this paper we compare FVC estimation models developed from Landsat-8 (L8), S2 and PS imagery. The reference grassland FVC dataset was obtained on the Paganella ski runs (46.15°N, 11.01°E, Italy) applying unsupervised classification to nadir grassland RGB photographs taken from 1.35 m above the soil. Fractional Response Models between reference FVC and 18 vegetation indices (VIs) extracted from satellite imagery were fitted and analysed. Then, leave-one-out cross validation and spatiotemporal change analysis were also performed. Our study confirms the robustness of the commonly used VIs based on the difference between NIR and the red wavelength region (R2 = 0.91 for EVI using S2 imagery) and indicate that VIs based on the red-edge spectral region are the best performing for PS imagery (R2 = 0.89 for RECI). Only medium to high spatial resolution imagery (S2 and PS) precisely mapped spatial patterns at the study site, since grasslands FVC varies at a fine scale. Previously available imagery at medium to low spatial and temporal resolution (e.g., L8) may still be interesting for analysis requiring long time-series of dat

Generative feature extraction from sentinel 1 and 2 data for prediction of forest aboveground biomass in the Italian Alps

—Aboveground biomass (AGB) is an important forest attribute directly linked to the forest carbon pool. The use of satellite remote sensing (RS) data has increased for AGB prediction due to their large footprint and low-cost availability. However, they have been limited due to saturation effect that leads to low prediction precision. In this article, we propose an innovative and dynamic architecture based on generative neural network that extracts target oriented generative features for predicting forest AGB using satellite RS data. These features are more resilient to mixed forest types and geographical conditions as compared to the traditional features and models. The effectiveness of the proposed features was assessed by experiments performed using multispectral, synthetic aperture radar, and combined dual-source datasets. The proposed model achieved best performance in terms of precision, model agreement, and overfitting as compared to the other conventional models for all analyzed datasets. The t-distributed stochastic neighbor embedding scatterplots of the generative features clearly show one dimension of the feature space associated with the target AGB. Feature importance analysis indicated that the produced generative features were more significant than the conventional analytical features. Also, the model provided a robust framework for homogeneous fusion of multisensor features from satellite RS data for predicting AG

Detection of grassland mowing frequency using time series of vegetation indices from Sentinel-2 imagery

5openInternationalItalian coauthor/editorManagement intensity deeply influences meadow structure and functioning, therefore affecting grassland ecosystem services. Conservation and management measures, including European Common Agricultural Policy subsidies, should therefore be based on updated and publicly available data about management intensity. The mowing frequency is a crucial trait to describe meadows management intensity, but the potential of using vegetation indices from Sentinel-2 imagery for its retrieval has not been fully exploited. In this work we developed on the Google Earth Engine platform a four-phases algorithm to identify mowing frequency, including i) vegetation index time-series computing, ii) smoothing and resampling, iii) mowing detection, and iv) majority analysis. Mowing frequency during 2020 of 240 ha of grassland fields in the Italian Alps was used for algorithm optimization and evaluation. Six vegetation indexes (EVI, GVMI, MTCI, NDII, NDVI, RENDVI783.740) were tested as input to the proposed algorithm. The Normalized Difference Infrared Index (NDII) showed the best performance, resulting in mean absolute error of 0.07 and 93% overall accuracy on average at the four sites used for optimization, at pixel resolution. A slightly lower accuracy (mean absolute error = 0.10, overall accuracy = 90%) was obtained aggregating the maps to management parcels. The algorithm showed a good generalization ability, with a similar performance between global and local optimization and an average mean absolute error of 0.12 and an overall accuracy of 89% on average on the sites not used for parameters optimization. The lowest accuracies occurred in intensively managed grasslands surveyed by one satellite orbit only. This study demonstrates the suitability of the proposed algorithm to monitor very fragmented grasslands in complex mountain ecosystems. Google Earth Engine was used to develop the model and will enable researchers, agencies and practitioners to easily and quickly apply the code to map grassland mowing frequency for extensive grasslands protection and conservation, for mowing event verification, or for forage system characterization.openAndreatta, Davide; Gianelle, Damiano; Scotton, Michele; Vescovo, Loris; Dalponte, MicheleAndreatta, D.; Gianelle, D.; Scotton, M.; Vescovo, L.; Dalponte, M

Automated machine learning driven stacked ensemble modelling for forest aboveground biomass prediction using multitemporal sentinel-2 data

Modelling and large-scale mapping of forest aboveground biomass (AGB) is a complicated, challenging and expensive task. There are considerable variations in forest characteristics that creates functional disparity for different models and needs comprehensive evaluation. Moreover, the human-bias involved in the process of modelling and evaluation affects the generalization of models at larger scales. In this paper, we present an automated machine learning (AutoML) framework for modelling, evaluation and stacking of multiple base models for AGB prediction. We incorporate a hyperparameter optimization procedure for automatic extraction of targeted features from multitemporal Sentinel-2 data that minimizes human-bias in the proposed modelling pipeline. We integrate the two independent frameworks for automatic feature extraction and automatic model ensembling and evaluation. The results suggest that the extracted target-oriented features have excessive contribution of red-edge and short-wave infrared spectrum. The feature importance scale indicates a dominant role of summer based features as compared to other seasons. The automated ensembling and evaluation framework produced a stacked ensemble of base models that outperformed individual base models in accurately predicting forest AGB. The stacked ensemble model delivered the best scores of R2 cv = 0.71 and RMSE = 74.44 Mgha-1 . The other base models delivered R2 cv and RMSE ranging between 0.38–0.66 and 81.27– 109.44 Mg ha-1 respectively. The model evaluation metrics indicated that the stacked ensemble model was more resistant to outliers and achieved a better generalization. Thus, the proposed study demonstrated an effective automated modelling pipeline for predicting AGB by minimizing human-bias and deployable over large and diverse forest area

Soil properties zoning of agricultural fields based on a climate-driven spatial clustering of remote sensing time series data

The identification of zones within an agricultural field that respond differently to environmental factors and agronomic management is a key requirement for the adoption of more precise and sustainable agricultural practices. Several approaches based on spatial clustering methods applied to different data sources, e.g. yield maps, proximal sensors and soil surveys, have been proposed in the last decades. The current availability of a huge amount of free remote sensing data allows to apply these approaches to agricultural areas where ground or proximal data are not available. However, in order to provide useful agronomic management information, it is essential that the zoning obtained by clustering is linked to the underlying spatial variability of soil properties. In this work we explore the hypothesis that the response of crop vigor to temporal climate variability, assessed by remote sensing data time series, selected to correspond to specific growth phases and seasonal climate patterns, provides indications on the variability of soil properties within agricultural fields, for both herbaceous and tree crops. NDVI time-series for 38 years (1984–2021) were obtained for fourteen non-irrigated herbaceous and tree crop fields in Central Italy, from multispectral satellites data (Landsat 5/7/8, Sentinel 2). The Standardized Precipitation-Evapotranspiration Index (SPEI) was used to classify time series into three climatic classes (dry/ normal/wet) for five different periods of the growth season, covering the main phenological phases. K-means clustering was used to identify patterns of crop growth from climatically classified image sets, as well as for all the bulked images for comparison (bulk clustering). Clustering results were compared with soil maps obtained from spatialized ground data, for soil texture (clay, silt and sand), soil organic matter and available soil water (ASW). The agreement between the different clustering results and soil maps was assessed by the Adjusted Rand Index. Agreement with soil maps varied depending on the field, the phenological phase considered and the soil property considered. Climate driven clustering from long, late growth season periods best matched soil properties, both for herbaceous and tree crops, despite being based on a limited number of images. The clustering from images spanning a longer growth period for dry years systematically outpaced the bulk clustering for silt, sand and ASW, while the clustering for normal climatic conditions was the best for organic matter. The performance of the matching between clustering and soil maps increased with soil variability significantly more (P < 0.05) than in the bulk clustering (mean slopes respectively 0.468 ± 0.167; 0.113 ± 0.270). The integration of the SPEI climatic index into the clustering procedure systematically improved the identification of zones with homogeneous soil properties, highlighting that a greater attention should be posed to the climate-crop-field interactions when using remotely sensed image

A method for continuous sub-annual mapping of forest disturbances using optical time series

Forest disturbances have a major impact on ecosystem dynamics both at local and global scales. Accordingly, it is important to acquire objective information about the location, nature and timing of such events to improve the understanding of their impact, update forest management policies and disturbance mitigation strategies. To this date, remotely sensed data have been widely used for the detection of stand replacing disturbances (SRD) such as windthrows and wildfires. In contrast, less effort has been devoted to the detection of non-stand replacing disturbances (NSRD), typically characterized by slower and gradual temporal dynamics. To address this gap, we propose a method for the automated detection of both SRD and NSRD. The proposed method can detect both past and recent disturbances, with a monthly temporal resolution, in a near real-time fashion by processing new images as they are acquired. Differently from existing approaches that handle the time series as a one-dimensional (1D) temporal trajectory, the method analyzes the sequence of images by organizing them in a two-dimensional (2D) grid-like structure. This representation allows us to model both the intra- and inter-annual variations of the time series taking advantage of the annual cyclical nature of the plant phenology. The method has been tested on study areas attacked by bark beetles achieving a user’s accuracy and producer’s accuracy of 0.91±0.08 and 0.81±0.07 (with 95% confidence intervals) for the disturbed areas, respectively

Mapping a European spruce bark beetle outbreak using sentinel-2 remote sensing data

Insect outbreaks affect forests, causing the deaths of trees and high economic loss. In this study, we explored the detection of European spruce bark beetle (Ips typographus, L.) outbreaks at the individual tree crown level using multispectral satellite images. Moreover, we explored the possibility of tracking the progression of the outbreak over time using multitemporal data. Sentinel-2 data acquired during the summer of 2020 over a bark beetle–infested area in the Italian Alps were used for the mapping and tracking over time, while airborne lidar data were used to automatically detect the individual tree crowns and to classify tree species. Mapping and tracking of the outbreak were carried out using a support vector machine classifier with input vegetation indices extracted from the multispectral data. The results showed that it was possible to detect two stages of the outbreak (i.e., early, and late) with an overall accuracy of 83.4%. Moreover, we showed how it is technically possible to track the evolution of the outbreak in an almost bi-weekly period at the level of the individual tree crowns. The outcomes of this paper are useful from both a management and ecological perspective: it allows forest managers to map a bark beetle outbreak at different stages with a high spatial accuracy, and the maps describing the evolution of the outbreak could be used in further studies related to the behavior of bark beetle

Extracting flowering phenology from grassland species mixtures using time-lapse cameras

Understanding the impacts of climate change on plant phenology is crucial for predicting ecosystem responses. However, accurately tracking the flowering phenology of individual plant species in grassland species mixtures is challenging, hindering our ability to study the impacts of biotic and abiotic factors on plant reproduction and plant-pollinator interactions. Here, we present a workflow for extracting flowering phenology from grassland species mixtures using near-surface time-lapse cameras. We used 89 image series acquired in plots with known species composition at the Jena trait-based experiment (Germany) to develop random forest classifiers, which were used to classify images and compute time series of flower cover for each species. The high temporal resolution of time-lapse cameras allowed to select images in proper light conditions, and to extract vegetation indices and texture metrics to improve discrimination among flowering species. The random forest classifiers showed a high accuracy in predicting the cover of Leucanthemum vulgare, Ranunculus acris, and Knautia arvensis flowers, whereas graminoid flowers were harder to predict due to their green-to-brownish colours. The proposed workflow can be applied in climate change studies, ecosystem functioning, plant community ecology, and biodiversity change research, including the investigation of effects of species richness on individual species' flowering phenology. Our method could be a valuable tool for understanding the impacts of climate change on plant reproduction and ecosystem dynamic

Large-area inventory of species composition using airborne laser scanning and hyperspectral data

5openInternationalInternational coauthor/editorTree species composition is an essential attribute in stand-level forest management inventories and remotely sensed data might be useful for its estimation. Previous studies on this topic have had several operational drawbacks, e.g., performance studied at a small scale and at a single tree-level with large fieldwork costs. The current study presents the results from a large-area inventory providing species composition following an operational area-based approach. The study utilizes a combination of airborne laser scanning and hyperspectral data and 97 field sample plots of 250 m2 collected over 350 km2 of productive forest in Norway. The results show that, with the availability of hyperspectral data, species-specific volume proportions can be provided in operational forest management inventories with acceptable results in 90% of the cases at the plot level. Dominant species were classified with an overall accuracy of 91% and a kappa-value of 0.73. Species-specific volumes were estimated with relative root mean square differences of 34%, 87%, and 102% for Norway spruce (Picea abies (L.) Karst.), Scots pine (Pinus sylvestris L.), and deciduous species, respectively. A novel tree-based approach for selecting pixels improved the results compared to a traditional approach based on the normalized difference vegetation index.openØrka, Hans Ole; Hansen, Endre Hofstad; Dalponte, Michele; Gobakken, Terje; Næsset, ErikØrka, H.O.; Hansen, E.H.; Dalponte, M.; Gobakken, T.; Næsset, E