Phenology Modelling and Forest Disturbance Mapping with Sentinel-2 Time Series in Austria

Assessment of spongy moth infestation impacts on forest productivity and carbon loss using the Sentinel-2 satellite remote sensing and eddy covariance flux data

Background

Deciduous forests in eastern North America experienced a widespread and intense spongy moth (Lymantria dispar) infestation in 2021. This study quantified the impact of this spongy moth infestation on carbon (C) cycle in forests across the Great Lakes region in Canada, utilizing high-resolution (10 × 10 m2) Sentinel-2 satellite remote sensing images and eddy covariance (EC) flux data. Study results showed a significant reduction in leaf area index (LAI) and gross primary productivity (GPP) values in deciduous and mixed forests in the region in 2021.

Results

Remote sensing derived, growing season mean LAI values of deciduous (mixed) forests were 3.66 (3.18), 2.74 (2.64), and 3.53 (2.94) m2 m-2 in 2020, 2021 and 2022, respectively, indicating about 24 (14)% reduction in LAI, as compared to pre- and post-infestation years. Similarly, growing season GPP values in deciduous (mixed) forests were 1338 (1208), 868 (932), and 1367 (1175) g C m-2, respectively in 2020, 2021 and 2022, showing about 35 (22)% reduction in GPP in 2021 as compared to pre- and post-infestation years. This infestation induced reduction in GPP of deciduous and mixed forests, when upscaled to whole study area (178,000 km2), resulted in 21.1 (21.4) Mt of C loss as compared to 2020 (2022), respectively. It shows the large scale of C losses caused by this infestation in Canadian Great Lakes region.

Conclusions

The methods developed in this study offer valuable tools to assess and quantify natural disturbance impacts on the regional C balance of forest ecosystems by integrating field observations, high-resolution remote sensing data and models. Study results will also help in developing sustainable forest management practices to achieve net-zero C emission goals through nature-based climate change solutions.

Harmonised statistics and maps of forest biomass and increment in Europe

Forest biomass is an essential resource in relation to the green transition and its assessment is key for the sustainable management of forest resources. Here, we present a forest biomass dataset for Europe based on the best available inventory and satellite data, with a higher level of harmonisation and spatial resolution than other existing data. This database provides statistics and maps of the forest area, biomass stock and their share available for wood supply in the year 2020, and statistics on gross and net volume increment in 2010-2020, for 38 European countries. The statistics of most countries are available at a sub-national scale and are derived from National Forest Inventory data, harmonised using common reference definitions and estimation methodology, and updated to a common year using a modelling approach. For those counties without harmonised statistics, data were derived from the State of Europe's Forest 2020 Report at the national scale. The maps are coherent with the statistics and depict the spatial distribution of the forest variables at 100 m resolution.

Spatial and remote sensing monitoring shows the end of the bark beetle outbreak on Belgian and north-eastern France Norway spruce (Picea abies) stands

In 2022, Europe emerged from eight of the hottest years on record, leading to significant spruce mortality across Europe. The particularly dry weather conditions of 2018 triggered an outbreak of bark beetles (Ips typographus), causing the loss of thousands of hectares of Norway spruce stands, including in Wallonia and North-eastern France. A methodology for detecting the health status of spruce was developed based on a dense time series of satellite imagery (Sentinel-2). The time series of satellite images allowed the modelling of the spectral response of healthy spruce forests over the seasons: a decrease in photosynthetic activity of the forest canopy causes deviations from this normal seasonal vegetation index trajectory. These anomalies are caused by a bark beetle attack and are detected automatically. The method leads in the production of an annual spruce health map of Wallonia and Grand-Est. The goal of this paper is to assess the damage caused by bark beetle using the resulting spruce health maps. A second objective was to compare the influence of basic variables on the mortality of spruce trees in these two regions. Lasted 6 years (2017-2022), bark beetle has destroyed 12.2% (23,674 ha) of the spruce area in Wallonia and Grand-Est of France. This study area is composed of three bioclimatic areas: Plains, Ardennes and Vosges, which have not been equally affected by bark beetle attacks. The plains were the most affected, with 50% of spruce forests destroyed, followed by the Ardennes, which lost 11.3% of its spruce stands. The Vosges was the least affected bioclimatic area, with 5.6% of spruce stands lost. For the most problematic sites, Norway spruce forestry should no longer be considered.

Towards precision forestry: Drought response from remote sensing-based disturbance monitoring and fine-scale soil information in Central Europe

Prolonged drought and susceptibility to biotic stressors induced an extensive calamity in Norway spruce (Picea abies (L.) Karst.) and widespread crown defoliation in European beech (Fagus sylvatica L.) in Central Europe. For future management decisions, it is crucial to link changes in canopy cover to site conditions. However, current knowledge on the role of soil properties for drought-induced forest disturbance is limited due to the scarcity and low spatial resolution of soil information. We present a fine-scale assessment on the role of soil properties for forest disturbance in Norway spruce and European beech derived from optical remote sensing. A forest disturbance modeling framework based on Sentinel-2 time series was applied on 340 km² in low mountain ranges of Central Germany. Spatio-temporal information on forest disturbance was calculated at 10 m spatial resolution in the period 2019-2021 and intersected with high-resolution soil information (1:10,000) based on roughly 2850 soil profiles. We found distinct differences in disturbed area, depending on soil type, texture, stoniness, effective rooting depth and available water capacity (AWC). For spruce, we found a polynomial relationship between AWC (R² = 0.7) and disturbance, with highest disturbed area (65 %) for AWC between 90 and 160 mm. Interestingly, we found no evidence for generally higher disturbance on shallow soils, although stands on the deepest soils were significantly less affected. Noteworthy, sites affected first did not necessarily exhibit highest proportions of disturbed area post-drought, indicating recovery or adaptation. We conclude that site- and species-specific understanding of drought impacts benefits from a combination of remote sensing and fine-scale soil information. Since our approach revealed which sites were affected first and most, it qualifies for prioritizing in situ monitoring activities to most vulnerable stands in acute drought conditions as well as for developing long-term strategies for reforestation and site-specific risk assessment for precision forestry.

Grassland Use Intensity Classification Using Intra-Annual Sentinel-1 and -2 Time Series and Environmental Variables

Detailed spatial data on grassland use intensity is needed in several European policy areas for various applications, e.g., agricultural management, supporting nature conservation programs, improving biodiversity strategies, etc. Multisensory remote sensing is an efficient tool to collect information on grassland parameters. However, there is still a lack of studies on how to process, combine, and implement large radar and optical image datasets in a joint observation framework to map grassland types on large heterogeneous study areas. In our study, we assessed the usefulness of 2521 Sentinel-1 and 586 Sentinel-2 satellite images and topographic data for mapping grassland use intensity. We focused on the distinction between intensively and extensively managed permanent grassland in a large heterogeneous study area in Slovenia. We provided dense Satellite Image Time Series (SITS) for 2017, 2018 and 2019 to identify important differences, e.g., management practices, between the two grassland types analysed. We also investigated the effectiveness of combining two different remote-sensing products, the optical Normalised Difference Vegetation Index (NDVI) and radar coherence. Grassland types were distinguished using an object-based approach and the Random Forest classification. With the use of SITS only, the models achieved poor performance in the case of cloudy years (2018). However, the performance improved with additional features (environmental variables). The feature selection method based on Mean Decrease Accuracy (MDA) provided a deeper insight into the high-dimensional multisensory SITS. It helped select the most relevant features (acquisition dates, environmental variables) that distinguish between intensive and extensive grassland types. The addition of environmental variables improved the overall classification accuracy by 7–15%, while the feature selection additionally improved the final overall classification accuracy (using all available features) by 2–3%. Although the reference dataset was limited (1259 training samples), the final overall classification accuracy was above 88% in all years analysed. The results show that the proposed Random Forest classification using combined multisensor data and environmental variables can provide better and more stable information on grasslands than single optical or radar data SITS on large heterogeneous areas. Therefore, a combined approach is recommended to distinguish different grassland types.

Detection of Tree Decline (Pinus pinaster Aiton) in European Forests Using Sentinel-2 Data

Moderate-resolution satellite imagery is essential to detect conifer tree decline on a regional scale and address the threat caused by pinewood nematode (PWN), (Bursaphelenchus xylophilus. This is a quarantine organism responsible for pine wilt disease (PWD), which has caused substantial ecological and economic losses in the maritime pine (Pinus pinaster) forests of Portugal. This study describes the first instance of a pre-operational algorithm applied to Sentinel-2 imagery to detect PWD-compatible decline in maritime pine. The Random Forest model relied on a pre-wilting and an in-season image, calibrated with data from a 24-month long field campaign enhanced with Worldview-3 data and the analysis of biological samples (hyperspectral reflectance, pigment quantification in needles, and PWN identification). Independent validation results attested to the good performance of the model with an overall accuracy of 95%, particularly when decline affects more than 30% of the 100 m2 pixel of Sentinel-2. Spectral angle mapper applied to hyperspectral measurements suggested that PWN infection cannot be separated from other drivers of decline in the visible-near infrared domain. Our algorithm can be employed to detect regional decline trends and inform subsequent aerial and field surveys, to further investigate decline hotspots.

Extension of an Open GEOBIA Framework for Spatially Explicit Forest Stratification with Sentinel-2

Spatially explicit information about forest cover is fundamental for operational forest management and forest monitoring. Although open-satellite-based earth observation data in a spatially high resolution (i.e., Sentinel-2, ≤10 m) can cover some information needs, spatially very high-resolution imagery (i.e., aerial imagery, ≤2 m) is needed to generate maps at a scale suitable for regional and local applications. In this study, we present the development, implementation, and evaluation of a Geographic Object-Based Image Analysis (GEOBIA) framework to stratify forests (needleleaf, broadleaved, non-forest) in Luxembourg. The framework is exclusively based on open data and free and open-source geospatial software. Although aerial imagery is used to derive image objects with a 0.05 ha minimum size, Sentinel-2 scenes of 2020 are the basis for random forest classifications in different single-date and multi-temporal feature setups. These setups are compared with each other and used to evaluate the framework against classifications based on features derived from aerial imagery. The highest overall accuracies (89.3%) have been achieved with classification on a Sentinel-2-based vegetation index time series (n = 8). Similar accuracies have been achieved with classification based on two (88.9%) or three (89.1%) Sentinel-2 scenes in the greening phase of broadleaved forests. A classification based on color infrared aerial imagery and derived texture measures only achieved an accuracy of 74.5%. The integration of the texture measures into the Sentinel-2-based classification did not improve its accuracy. Our results indicate that high resolution image objects can successfully be stratified based on lower spatial resolution Sentinel-2 single-date and multi-temporal features, and that those setups outperform classifications based on aerial imagery only. The conceptual framework of spatially high-resolution image objects enriched with features from lower resolution imagery facilitates the delivery of frequent and reliable updates due to higher spectral and temporal resolution. The framework additionally holds the potential to derive additional information layers (i.e., forest disturbance) as derivatives of the features attached to the image objects, thus providing up-to-date information on the state of observed forests.

Earth Observation for Phenological Metrics (EO4PM): Temporal Discriminant to Characterize Forest Ecosystems

The study of vegetation phenology has great relevance in many fields since the importance of knowing timing and shifts in periodic plant life cycle events to face the consequences of global changes in issues such as crop production, forest management, ecosystem disturbances, and human health. The availability of high spatial resolution and dense revisit time satellite observations, such as Sentinel-2 satellites, allows high resolution phenological metrics to be estimated, able to provide key information from time series and to discriminate vegetation typologies. This paper presents an automated and transferable procedure that combines validated methodologies based on local curve fitting and local derivatives to exploit full satellite Earth observation time series to produce information about plant phenology. Multivariate statistical analysis is performed for the purpose of demonstrating the capacity of the generated smoothed vegetation curve, temporal statistics, and phenological metrics to serve as temporal discriminants to detect forest ecosystems processes responses to environmental gradients. The results show smoothed vegetation curve and temporal statistics able to highlight seasonal gradient and leaf type characteristics to discriminate forest types, with additional information about forest and leaf productivity provided by temporal statistics analysis. Furthermore, temporal, altitudinal, and latitudinal gradients are obtained from phenological metrics analysis, which also allows to associate temporal gradient with specific phenophases that support forest types distinction. This study highlights the importance of integrated data and methodologies to support the processes of vegetation recognition and monitoring activities.

A First Assessment of Canopy Cover Loss in Germany’s Forests after the 2018–2020 Drought Years

Central Europe was hit by several unusually strong periods of drought and heat between 2018 and 2020. These droughts affected forest ecosystems. Cascading effects with bark beetle infestations in spruce stands were fatal to vast forest areas in Germany. We present the first assessment of canopy cover loss in Germany for the period of January 2018–April 2021. Our approach makes use of dense Sentinel-2 and Landsat-8 time-series data. We computed the disturbance index (DI) from the tasseled cap components brightness, greenness, and wetness. Using quantiles, we generated monthly DI composites and calculated anomalies in a reference period (2017). From the resulting map, we calculated the canopy cover loss statistics for administrative entities. Our results show a canopy cover loss of 501,000 ha for Germany, with large regional differences. The losses were largest in central Germany and reached up to two-thirds of coniferous forest loss in some districts. Our map has high spatial (10 m) and temporal (monthly) resolution and can be updated at any time.

Mapping annual land disturbance and reclamation in rare-earth mining disturbance region using temporal trajectory segmentation

Rare-earth mining has caused extensive damage to soil, vegetation, and water, significantly threatening ecosystems. Monitoring environmental disturbance caused by rare-earth mining is necessary to protect the ecological environment. A spatiotemporal remote sensing monitoring method for mining to reclamation processes in a rare-earth mining area using multisource time-series satellite images is described. In this study, the normalized difference vegetation index (NDVI) is used to evaluate the mining impact. Regression analysis is conducted to relate the HJ-1B CCD and Landsat 5/8 data to reduce the NDVI error related to sensor differences between different datasets. The analysis method of NDVI trajectory data of ground objects is proposed, and areas of environmental disturbance caused by rare-earth mining are identified. Pixel-based trajectories were used to reconstruct the temporal evolution of vegetation, and a temporal trajectory segmentation method is established based on the vegetation changes in different disturbance stages. The temporal trajectory of the rare-earth disturbance points is segmented to extract features in each stage to obtain the disturbance year, recovery year, and recovery cycle and evaluate the vegetation recovery after rare-earth mining disturbance. We applied the method to a stack of 20 multitemporal images from 2000 to 2019 to analyze vegetation disturbance due to rare-earth mining and vegetation recovery in the upper reaches of the Guangdong-Hong Kong-Macao Greater Bay Area, China. The results show the following. (1) Mining industry in the study area experienced rapid expansion before 2008, but growth slowed since the policies implemented by the government since 2009 to restrict rare-earth mining. (2) The continuous influence to the land caused by rare-earth mining can last for decades; however, the reclamation activities shorten the recovery cycle of mining land from 5 to 3 years.

Applicability of Smoothing Techniques in Generation of Phenological Metrics of Tectona grandis L. Using NDVI Time Series Data

Information on phenological metrics of individual plant species is meager. Phenological metrics generation for a specific plant species can prove beneficial if the species is ecologically or economically important. Teak, a dominating tree in most regions of the world has been focused on in the present study due to its multiple benefits. Forecasts on such species can attain a substantial improvement in their productivity. MODIS NDVI time series when subjected to statistical smoothing techniques exhibited good output with Tukey’s smoothing (TS) with a low RMSE of 0.042 compared to single exponential (SE) and double exponential (DE). Phenological metrics, namely, the start of the season (SOS), end of the season (EOS), maximum of the season (MAX), and length of the season (LOS) were generated using Tukey-smoothed MODIS NDVI data for the years 2003–2004 and 2013–2014. Post shifts in SOS and EOS by 14 and 37 days respectively with a preshift of 28 days in MAX were observed in the year 2013–2014. Preshift in MAX was accompanied by an increase in greenness exhibiting increased NDVI value.LOS increased by 24 days in the year 2013–2014, showing an increase in the duration of the season of teak. Dates of these satellite-retrieved phenological occurrences were validated with ground phenological data calculated using crown cover assessment. The present study demonstrated the potential of a spatial approach in the generation of phenometrics for an individual plant species, which is significant in determining productivity or a crucial trophic link for a given region.