Accurate information regarding forest tree species composition is useful for a wide range of applications, both for forest management and scientific research. Remote sensing is an efficient tool for collecting spatially explicit information on forest attributes. With the launch of the Sentinel-2 mission, new opportunities have arisen for mapping tree species owing to its spatial, spectral, and temporal resolution. The short revisit cycle (five days) is crucial in vegetation mapping because of the reflectance changes caused by phenological phases. In our study, we evaluated the utility of the Sentinel-2 time series for mapping tree species in the complex, mixed forests of the Polish Carpathian Mountains. We mapped the following nine tree species: common beech, silver birch, common hornbeam, silver fir, sycamore maple, European larch, grey alder, Scots pine, and Norway spruce. We used the Sentinel-2 time series from 2018, with 18 images included in the study. Different combinations of Sentinel-2 imagery were selected based on mean decrease accuracy (MDA) and mean decrease Gini (MDG) measures, in addition to temporal phonological pattern analysis. Tree species discrimination was performed using the Random Forest classification algorithm. Our results showed that the use of the Sentinel-2 time series instead of single date imagery significantly improved forest tree species mapping, by approximately 5–10% of overall accuracy. In particular, combining images from spring and autumn resulted in better species discrimination.

Accurate information on forest tree species composition is vital for various scientific applications, as well as for forest inventory and management purposes. Country-wide, detailed species maps are a valuable resource for environmental management, conservation, research, and planning. Here, we performed the classification of 16 dominant tree species and genera in Poland using time series of Sentinel-2 imagery. To generate comprehensive spectral–temporal information, we created Sentinel-2 seasonal aggregations known as spectral–temporal metrics (STMs) within the Google Earth Engine (GEE). STMs were computed for short periods of 15–30 d during spring, summer, and autumn, covering multi-annual observations from 2018 to 2021. The Polish Forest Data Bank served as reference data, and, to obtain robust samples with pure stands only, the data were validated through automated and visual inspection based on very-high-resolution orthoimagery, resulting in 4500 polygons serving as training and test data. The forest mask was derived from available land cover datasets in GEE, namely the ESA WorldCover and Dynamic World dataset. Additionally, we incorporated various topographic and climatic variables from GEE to enhance classification accuracy. The random forest algorithm was employed for the classification process, and an area-adjusted accuracy assessment was conducted through cross-validation and test datasets. The results demonstrate that the country-wide forest stand species mapping achieved an accuracy exceeding 80 %; however, this varies greatly depending on species, region, and observation frequency. We provide freely accessible resources, including the forest tree species map and training and test data: https://doi.org/10.5281/zenodo.10180469 (Grabska-Szwagrzyk, 2023a).

Site index is the most commonly used measure of potential site productivity, providing key information for forest management practices. It is determined using species-specific site index models that take into account climatic and edaphic factors. To reliably estimate the site index, appropriate models are necessary. In Poland, however, outdated guidelines, i.e., yield tables, are used to determine site classes, which result in the inappropriate estimation of height growth and increments of stands. Therefore, the aim of this study is to develop new site index models for the main forest forming tree species in Poland, in a total of eight species. For the development of site index models, we used growth trajectories of 3052 sample trees, representing the whole range of geographic locations and site conditions. Five dynamic models were selected and parametrized to develop the site index models. The models were evaluated using quantitative measures of goodness of fit (MAE, R2, and AIC), the analysis of residuals, and the assessment of how the model reflects the biological phenomena of height growth. Results showed that depending on the species, different models have the highest predictive ability. There are significant differences in results using traditional yield tables and developed site index models. For most of the species, the largest differences characterized either the youngest or the oldest age classes. These differences can be attributed to the changes in growth conditions from the time when yield tables were developed. Growth dynamics of forest stands may also show spatial variability, thus, in future research additional site variables and, regional variability should be taken into account.

The changes in vegetation phenology observed in recent decades may have significant implications for carbon, water, and energy fluxes as well as ecosystem productivity. Phenology is an important indicator of a changing environment, and studying species-specific phenology can help understand their responses to varying climate conditions. While satellite-based phenology provides long-term information for large areas, it describes different traits than ground-based phenology. Nonetheless, emerging research on remote sensing-based phenology has shown its usefulness in determining inter- and intra-specific differences, as well as interannual changes. In this study, we examined the potential of the Sentinel-2 time series from 2018 to 2022 in assessing the spring phenology of seven temperate tree species. Three vegetation indices were calculated, namely, the Medium Resolution Imaging Spectrometer (MERIS) Terrestrial Chlorophyll Index, the Normalized Difference Vegetation Index, and the Enhanced Vegetation Index. We examined interspecific variations in the start of the season (SOS) and interannual variability. Additionally, we explored intraspecific variations in spring phenology among common beech stands across an elevation gradient. Our findings reveal substantial differences in Sentinel-2-derived phenology across species, with larch, birch, and hornbeam displaying an earlier SOS compared with oak and alder. The latest SOS dates were observed in 2021, whereas the earliest occurred in 2018 and the difference exceeded 20 days. It corresponded to the diverse spring meteorological conditions in the examined years, with the highest heat accumulation in 2018 and the lowest in 2021. Additionally, the three examined indices showed different phenological curves, with MERIS Terrestrial Chlorophyll Index exhibiting distinct patterns with a delayed start of the spring increase compared with the other two indices. The relationship between elevation and SOS dates was generally positive, but there was high variability in some years. The earliest SOS dates occurred at ~600 m a.s.l. and larger interannual variations were observed for stands located on the highest elevations. These results underscore the potential of Sentinel-2 in monitoring the spring phenology of temperate species, offering insights into inter- and intra-specific variations, as well as year-to-year changes. In the future, longer Sentinel-2 time series will improve estimates and facilitate the tracking of long-term changes. Such advancements hold relevance for forest management, carbon budget predictions, and understanding ecosystem dynamics.

Climate change and severe extreme events, i.e., changes in precipitation and higher drought frequency, have a large impact on forests. In Poland, particularly Norway spruce and Scots pine forest stands are exposed to disturbances and have, thus experienced changes in recent years. Considering that Scots pine stands cover approximately 58% of forests in Poland, mapping these areas with an early and timely detection of forest cover changes is important, e.g., for forest management decisions. A cost-efficient way of monitoring forest changes is the use of remote sensing data from the Sentinel-2 satellites. They monitor the Earth’s surface with a high temporal (2–3 days), spatial (10–20 m), and spectral resolution, and thus, enable effective monitoring of vegetation. In this study, we used the dense time series of Sentinel-2 data from the years 2015–2019, (49 images in total), to detect changes in coniferous forest stands dominated by Scots pine. The simple approach was developed to analyze the spectral trajectories of all pixels, which were previously assigned to the probable forest change mask between 2015 and 2019. The spectral trajectories were calculated using the selected Sentinel-2 bands (visible red, red-edge 1–3, near-infrared 1, and short-wave infrared 1–2) and selected vegetation indices (Normalized Difference Moisture Index, Tasseled Cap Wetness, Moisture Stress Index, and Normalized Burn Ratio). Based on these, we calculated the breakpoints to determine when the forest change occurred. Then, a map of forest changes was created, based on the breakpoint dates. An accuracy assessment was performed for each detected date class using 861 points for 46 classes (45 dates and one class representing no changes detected). The results of our study showed that the short-wave infrared 1 band was the most useful for discriminating Scots pine forest stand changes, with the best overall accuracy of 75%. The evaluated vegetation indices underperformed single bands in detecting forest change dates. The presented approach is straightforward and might be useful in operational forest monitoring.

Forest stand reflectance at the canopy level results from various factors, such as vegetation chemical properties, leaf morphology, canopy structure, and tree sizes. These factors are dependent on the species, age, and health statuses of trees, as well as the site conditions. Sentinel-2 imagery with the high spatial, spectral, and temporal resolution, has enabled analysis of the relationships between vegetation properties and their spectral responses at large spatial scales. A comprehensive study of these relationships is needed to understand the drivers of vegetation spectral patterns and is essential from the point of view of remote sensing data interpretation. Our study aimed to quantify the site and forest parameters affecting forest stands reflectance. The analysis was conducted for common beech-, silver fir- and Scots pine-dominated stands in a mountainous area of the Polish Carpathians. The effect of stands and site properties on reflectance in different parts of the growing season was captured using the dense time series provided by Sentinel-2 from 2018–2019. The results indicate that the reflectance of common beech stands is mainly influenced by elevation, particularly during spring and autumn. Other factors influencing beech stand reflectance include the share of the broadleaved understory, aspect, and, during summer, the age of stands. The reflectance of coniferous species, i.e., Scots pine and silver fir, is mainly influenced by the age and stand properties, namely the crown closure and stand density. The age is a primary driver for silver fir stands reflectance changes, while the stand properties have a large impact on Scots pine stands reflectance. Also, the understory influences Scots pine stands reflectance, while there appears to be no impact on silver fir stands. The influence of the abovementioned factors is highly diverse, depending on the used band and time of the season.