Forest Communities Spatial Modeling as a Basis for Assessing the Sequestration Potential of Ecosystems (Republic of Tatarstan, Russia)

1. Introduction

Terrestrial and water ecosystems play a significant role in carbon sequestration through the creation and retention of organic carbon mass within plant biomass and soil. To date, there exists a plethora of publications dedicated to various aspects of this issue. A comprehensive review of numerous methods and practical cases is provided in the work by Kurbatova [1]. A common premise across these studies is the identification of desired parameters for individual objects (for example, carbon content in the wood of a specific tree species), which are then extrapolated to groups of objects (such as all trees of this species in similar ecological conditions) and multiplied by the number of trees to estimate the overall carbon stock in a forest stand. In this case, the commonality among all objects in the group is the plant species concept, which doesn't pose difficulties, since the taxonomy of higher vascular plants is well-established for the most part.

When referring to a volumetric area object, such as a soil layer, a series of samples, where carbon content is determined in the laboratory, is extrapolated to the entire area occupied by this type of soil, assuming equal thickness (for instance) of the humus horizon. Here, two problematic aspects emerge: object classification (in our case, soil classification), and the determination of their area distribution (direct mapping or mapping through decryption) [2,3].

When addressing the carbon exchange between the entire ecosystem as a whole and the atmosphere, similar to the previous example, it is necessary to classify ecosystems in one way or another, and then determine the distribution of each class for area calculation.

Thus, one of the "cornerstones" in creating carbon balance models is the classification of objects for extrapolation of field observations. Vegetation changes over time and space as a result of ecological processes acting on plant populations and communities at different temporal and spatial scales, making it one of the integral indicators of the ecosystem as a whole [4]. This is because the vegetation component of ecosystems is formed depending on both abiotic factors (moisture, light, soil fertility) and anthropogenic influences. Hence, we propose utilizing the spatial model of vegetation cover as a basis for evaluating the sequestration potential of ecosystems.

Recently, renewed interest in vegetation classification has emerged worldwide, and efforts have been undertaken both nationally and internationally to develop new classification systems using standardized procedures [5,6]. Moreover, there's a growing interest in harmonizing approaches and standardizing the informational content of classifications [7]. This interest is driven by the practical necessity of using vegetation typologies and the increasing recognition of their scientific bases [8]. The classification problem is closely tied to assigning a geobotanical description of a particular community to a syntaxonomic unit of an existing vegetation classification. Experience exists in solving this task using expert systems [9].

Very recently, ideas have emerged regarding the amalgamation of vegetation classification procedures and the solutions provided by expert systems, employing the apparatus of fuzzy sets [10]. The proposed algorithms align with the notions of continuous vegetation cover and the individual behavior of plant species concerning environmental factors.

Therefore, all these research directions aim to minimize the subjectivity of researchers at all stages of creating vegetation cover maps from the development of classification to assigning a specific geographical point to one type of ecosystem or another.

The analysis of terrestrial ecosystems, particularly the classification and mapping of vegetation, is a pivotal task for managing natural resources [11]. Vegetation serves as the foundation for all living beings, playing an indispensable role in various ecological and environmental processes. In recent years, remote sensing technologies have substantially evolved, offering a robust tool for accurate and efficient vegetation mapping. A critical aspect of leveraging remote sensing data for vegetation mapping is the choice of classification approach, which largely dictates the accuracy and reliability of the resulting maps. Among the various classification approaches, unsupervised classification has emerged as a vital technique in minimizing researcher subjectivity at all stages of vegetation mapping, from developing a classification scheme to assigning a particular geographic point to a specific ecosystem type [12,13,14,15,16]. Unlike supervised classification, which relies on prior knowledge and extensive field samples, unsupervised classification operates without any a priori knowledge, solely depending on spectrally pixel-based statistics [13]. This inherent characteristic of unsupervised classification significantly reduces the potential bias that may arise from researchers' prior knowledge or assumptions, thereby enhancing the objectivity of the classification process. Unsupervised classification encompasses a variety of clustering algorithms, including rule-based, neural network, and SVM-based approaches, each with its unique advantages and applications in different scenarios. The utilization of unsupervised classification is not confined to preparing for supervised classifiers but extends to visualizing and monitoring similar areas in a scene, capturing a wide range of land use signatures over areas with extensive and unknown ranges of vegetation types. This broad capture capability underscores the versatility and efficacy of unsupervised classification in handling diverse ecosystems and varied vegetation types. Furthermore, the cost-effectiveness of unsupervised classification, in terms of both time and financial resources, is another significant advantage, especially when compared to supervised classification that necessitates a large number of field samples. The reduced need for extensive field sampling is particularly beneficial in scenarios such as agricultural production management, where timely access to spatial distribution information of crops is of great importance [17,18].

In light of the above, optimized unsupervised classification techniques have been proposed to further enhance the accuracy and efficiency of vegetation mapping. For instance, a two-step approach combining Jenks classification and Agglomerative Hierarchical Clustering (AHC) has been implemented to derive vegetation cover types from SAVI data [19]. Such advancements in unsupervised classification techniques are indicative of the ongoing efforts to minimize researcher subjectivity, thereby contributing to more accurate and reliable vegetation mapping.

The burgeoning interest and continuous advancements in unsupervised classification underscore its importance in the realm of vegetation mapping [17]. By minimizing researcher subjectivity, unsupervised classification paves the way for more objective, accurate, and reliable vegetation maps, which are indispensable for a myriad of ecological, environmental, and agricultural applications. The exploration and adoption of unsupervised classification techniques, coupled with the integration of other emerging technologies, hold promise for further advancing the field of vegetation mapping, fostering a deeper understanding of terrestrial ecosystems, and facilitating more informed decision-making in natural resource management [20].

The main goal of this study is to create a spatial model of forest vegetation in the classes of the Braun-Blanquet system. This classification approach is based on information about the composition and abundance of plant species in a community. Each species is characterized by its own functional characteristics, which ultimately can give a knowledge about a particular community productivity, which in turn can be used to judge the sequestration capacity of ecosystems.

2. Materials and Methods

2.1. Study Area

The study area is located within the Republic of Tatarstan (Russian Federation) situated in the eastern part of European Russia (Figure 1). The territory is divided by river valleys of the Volga, Kama and Vyatka into natural areas [21]including landscapes related to the Eastern European boreal subtaiga and Eastern semi-humid broadleaved forests, the boundary between them being determined by the isoline of hydrothermal coefficient of 1 [22]. In zonal terms, the northern part of the territory belongs to the ecoregion of Sarmatian mixed forests, the southern part to the ecoregion of Eastern European forest-steppes [23].

The area includes various landscape types and vegetation units ranging from the formation of the southern taiga to the fragments of typical steppe [24]. In the river valleys on sandy or sandy-clayey substrate various types of pine forests are formed. Meadow formations have developed as a result of forest clearing and mowing. Steppe has remained mostly on the steep rocky slopes that are not suitable for agricultural development. The predominant types of wetlands are fens.

The flora and vegetation of the Republic of Tatarstan have been studied mainly by researchers from Kazan University. This university has a long tradition of floristic and vegetation studies since the 19th century in connection with S. Korshinsky, A. Gordyagin and M. Markov. Much phytosociological data has been collected during more than 120 years by various authors [25].

2.2. Data Acquisition

The data acquisition phase of this study primarily revolved around the utilization of multi-temporal satellite imagery from the Landsat 8 and Landsat 9 platforms. These satellites are part of the established Landsat program, which has been a critical resource for earth observation, aiding in various environmental monitoring [26] and resource management projects [27].

Landsat 8, launched by NASA in 2013, is equipped with two sensors: the Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS). The OLI captures imagery across nine spectral bands, covering visible, near-infrared, and shortwave infrared wavelengths, with a spatial resolution of 30 meters for most bands and 15 meters for the panchromatic band. Conversely, the TIRS captures thermal infrared imagery with a resolution of 100 meters, which is resampled to 30 meters for product delivery.

Landsat 9, akin to Landsat 8, continues the mission of providing high-quality imagery for earth observation. With similar sensors, Landsat 9 ensures a consistent data stream, which is indispensable for long-term environmental monitoring [28].

Several advantages underpin the choice of Landsat imagery for this study. The temporal resolution, with a revisit time of 16 days for both Landsat 8 and Landsat 9, provides an 8-day revisit at the equator when combined. This frequency allows for the capture of temporal dynamics within the forest communities over the study period from 2018 to 2022. Additionally, the 30-meter spatial resolution of Landsat imagery offers a balanced trade-off between detail and coverage, enabling the analysis of forest structure and composition at a landscape scale. Furthermore, the wide range of spectral bands available facilitates the computation of various vegetation indices crucial for assessing different aspects of vegetation health and structure.

In the image processing phase, the raw Landsat imagery underwent a series of pre-processing steps to correct atmospheric and geometric distortions. This was achieved through the preprocess method in the eemont [29] extension of Google Earth Engine (GEE) [30], which provides functionalities for atmospheric correction and cloud masking, thereby improving the quality and consistency of the imagery. Temporal filtering was also performed to include data from April to November for the years 2018 to 2022, to account for seasonal variations in vegetation cover and exclude winter months where snow cover could interfere with the analysis. Subsequently, a comprehensive suite of vegetation indices was derived from the pre-processed imagery using the spectralIndices method in eemont. These indices serve as proxies for various vegetation attributes and are critical for the unsupervised classification task. The Landsat 8 and Landsat 9 imagery were then merged to form a continuous time-series dataset. This dataset was utilized to compute statistical composites like median, mean, maximum, and standard deviation for each of the vegetation indices across the temporal dimension. These composites encapsulate the central tendency and variability of vegetation conditions over the study period, providing a robust dataset for the unsupervised classification.

2.3. Google Earth Engine Platform

The GEE platform played a quintessential role in the data processing and analysis stages of this study. GEE is a cloud-based platform designed for the planetary-scale environmental data analysis. It provides a robust infrastructure necessary for processing large datasets, which is essential given the volume of multi-temporal satellite imagery utilized in this study. GEE houses an extensive archive of satellite imagery and geospatial datasets, which is continuously updated, thereby providing a rich repository of data for environmental studies. For this investigation, the Landsat 8 and Landsat 9 imagery were directly accessed and processed within the GEE environment, which significantly streamlined the data handling process.

One of the significant advantages of using GEE [31] is its ability to handle and process large volumes of data in a parallelized and efficient manner. The platform employs a server-side processing paradigm, which means that the data does not need to be downloaded to local machines for analysis, thus circumventing the traditional bottlenecks associated with large-scale data processing. This is particularly advantageous for studies like ours, where the analysis spans a large geographical area over an extended period.

The GEE platform supports a variety of data processing and analysis tasks essential for this study. These include image collection filtering, image reduction, computation of vegetation indices, and unsupervised classification. The platform also provides a range of pre-built algorithms and machine learning models, alongside the ability to implement custom algorithms, which was pivotal for the unsupervised classification task undertaken.

Furthermore, GEE provides a user-friendly interface and a powerful scripting environment through its Code Editor, which supports both JavaScript and Python. For this study, the Python programming language was employed for scripting, facilitated by the ee and eemont libraries. These libraries provided the necessary tools and extended functionalities for image processing, feature extraction, and statistical analysis, all of which were crucial for the successful execution of the unsupervised classification of forest communities.

Moreover, the GEE platform also enables easy sharing and exporting of results, which is crucial for collaboration and further analysis. The results from the unsupervised classification were exported to Google Drive, which allowed for a seamless transition to the post-classification analysis stage of the study.

2.4. Vegetation Indices Computation

The computation of vegetation indices (VIs) is pivotal in harnessing the spectral information embodied in the satellite imagery. VIs serves as proxies for various biophysical parameters and are instrumental in delineating vegetation characteristics. In this study, a gamut of VIs was derived from the pre-processed Landsat 8 and Landsat 9 imagery, leveraging the spectralIndices method within the eemont library in Google Earth Engine.

In our comprehensive analysis of vegetation indices, we adopted a meticulous approach to categorize and scrutinize various indices based on their underlying characteristics and applications. The vegetation indices were classified into distinct categories, each representing a specific aspect of vegetation characteristics captured via remote sensing (Table 1).

The array of VIs computed in this study furnishes a robust dataset encapsulating a myriad of vegetation attributes. These indices, each with its unique domain of application, collectively contribute to a comprehensive feature set forming the basis for the unsupervised classification of forest communities within the study area. Through the judicious computation and utilization of these VIs, a rich tapestry of information regarding vegetation dynamics over time and across the study area was woven, paving the avenue for an in-depth analysis of forest communities through unsupervised classification.

2.5. Unsupervised Classification

The mapping of forest communities in an unsupervised manner was a focal point in this analysis, aimed at uncovering the intrinsic structure and spatial distribution of vegetation types within the study area. The technique employed was rooted in the Weka X-Means clustering algorithm [69] facilitated by the GEE platform, extending the functionality of the traditional K-means algorithm by adaptively ascertaining the optimal number of clusters within a designated range.

2.5.1. Data Pre-Processing

A cornerstone to the ensuing classification was the meticulous pre-processing of data. The ensemble of vegetation indices computed from prior steps was amalgamated to form a cohesive dataset. Subsequently, statistical reducers were employed to distill the temporal essence of the data by generating median, mean, maximum, and standard deviation composites for each index. This was pivotal in encapsulating the temporal dynamics of vegetation attributes over the study timeframe. To refine the dataset further, a mask was applied utilizing a dataset from the Global Land Analysis and Discovery (GLAD) laboratory at the University of Maryland [86], effectively sieving out non-forested areas and thus honing the focus on forested regions.

2.5.2. Training Dataset Preparation

The basis of the clustering process was the preparation of a robust training dataset. A total of 20,000 pixels were randomly sampled across the study area, each pixel representing a data point with a suite of features encapsulated in the vegetation indices and statistical composites. This dataset, rendered at a spatial resolution of 30 meters, served as the crucible for training the clusterer, providing a diverse yet representative glimpse into the spectral and temporal heterogeneity of the vegetation within the study area.

2.5.3. Clustering Process

With a well-curated training dataset at hand, the stage was set for the clustering process. The Weka X-Means clustering algorithm was trained, iterating through a specified range of clusters, specifically from 36 to 50, to dynamically ascertain the optimal number of clusters that best encapsulated the inherent structure of the data. Once trained, the clustering model was applied in a stratified manner to each grid within the study area, engendering a separate classified image for each grid. This stratified approach ensured a localized analysis, catering to the unique vegetative composition and spatial heterogeneity inherent within each grid.

2.5.4. Mosaic Creation and Export

Post-classification, the separate classified images for each grid were melded into a single coherent image through a mosaicking process. This process, facilitated by the mosaic method in GEE, seamlessly stitched together the classified images based on the region boundaries, yielding a comprehensive classified image of the entire study area.

2.6. Field Data Collection and Processing

The primary source of data for this study validation and verification was procured from the extensive "Flora" database [25], which houses a vast collection of over 17,000 relevés. A selective extraction process executed to acquire relevés documented post the year 2000, which also contained geographical references, ensuring relevance and contemporaneity of the data. For maintaining a standardized taxonomy across the dataset, the nomenclature of plant species was meticulously aligned with the globally recognized "The Plant List" database [70]. A deliberate exclusion of bryophyte species from the selection carried out to maintain a focused approach towards vascular plant communities. The classification of the relevés executed at the level of vegetation classes as delineated in the Braun-Blanquet system [87,88]. Automation played a crucial role in the classification process. The expert systems EuroVeg Checklist [9] employed within the JUICE 7.1 software package [89] to facilitate a precise and efficient automated classification. The integration of this advanced expert system streamlined the classification process, ensuring an objective and methodical analysis of the vegetation data.

This rigorous methodology in data extraction, nomenclature standardization, and classification, bolstered by automated expert system, constituted a robust framework for field data collection and processing. This framework was instrumental in ensuring the accuracy and reliability of the classification, which in turn significantly contributed to a comprehensive understanding of vegetation and habitat distributions, aligning with the broader objectives of this study.

3. Results and discussion

The conducted cluster analysis facilitated the delineation of 44 unique statistically distinguishable classes, which formed the foundation of the raster model of forest communities in the Republic of Tatarstan (Figure 2). In the derived model, pixel values are present, some of which exhibit minimal occurrence, introducing a level of 'noise' into the model, while others are frequently occurring values (Figure 3). Such 'noisy' values either need to be eliminated or associated with another class. To address this issue, an approach based on "tail cutting" was adopted – at the 1 and 99 percentiles, which enabled the filtering out of all minimally occurring classes (values 3, 10, 34, 40).

To correlate the classes of the raster model with the existing Braun-Blanquet classification system, two approaches were employed – automated and expert. For the automated search for correspondences, an approach based on the analysis of frequency tables of values from the raster model and corresponding coordinates of geobotanical descriptions was utilized.

The implementation of automated identification, to determine which values of the raster model correspond to a specific type according to the classification of plant communities, was carried out in several stages. Initially, classes of field data were identified, where a vast variety of raster values were observed, and classes where they were more homogeneous (Table 2, Figure 4). In the subsequent stage, grouping and aggregation were performed. The grouping of raster values was conducted based on the community type according to the Braun-Blanquet classification to unveil common patterns (Figure 2) and explore the possibility that certain raster values are more prevalent in specific classes of field data. Following this, statistical analysis was employed to uncover significant associations between raster values and classes of field data. This allowed for the recognition of several patterns that may indicate a correspondence between groups of raster values and specific classes of field data. The results of this analysis can be reviewed in the Table S1 of the Supplementary Materials.

The meticulous procedure of aggregating the raster model values has been instrumental in unequivocally identifying the principal community types, which are denoted as FAG (Carpino-Fagetea sylvaticae), BRA (Brachypodio pinnati-Betuletea pendulae), PIC (Vaccinio-Piceetea), ALN (Alnetea glutinosae), POP (Alno glutinosae-Populetea albae), and PUB (Quercetea pubescentis). It’s pertinent to note that the class QUE (Quercetea robori-petraeae) entirely aligns with the FAG group based on the raster model values. The remaining classes either showcase low statistical significance or have been included erroneously due to inaccuracies in the spatial referencing of the geobotanical descriptions. However, these raster model values were not discarded; instead, they were expertly allocated.

Post aggregation, a Z-score assessment was meticulously conducted to measure the deviation of a particular observation from the mean value in a distribution. This is calculated as the ratio of the difference between the observation and the mean value to the standard deviation. The Z-score is pivotal for identifying anomalous values, discerning data percentiles, and determining the position of an observation within the mean value's confidence interval. Additionally, it facilitates data normalization, thus simplifying the analysis and comparison of data. The application of this assessment has enabled a clear correlation of at least one (main) raster model class with a forest community as per the Braun-Blanquet classification. In certain instances, where a high degree of correspondence between multiple raster model classes and one type according to Braun-Blanquet was observed, both raster values were associated with the corresponding type. The remaining values were expertly distributed based on the geobotanical descriptions of areas, aiming to align them with one of the 6 types of forest communities as per the Braun-Blanquet classification found within the study area. In cases of derivative and deadwood forests, the raster model values were categorized under the class OTHERS.

The resultant distribution of pixels across classes of forest communities presented in Table 3. A comprehensive cartographic representation of the raster model of forest communities within the Republic of Tatarstan has been illustrated in Figure 5. This representation is pivotal as it offers a spatial visualization of the forest community distribution, encapsulating the diverse forest communities within the specified region. It provides a visual aid that significantly contributes to a deeper understanding and interpretation of the forest community distributions and their classifications, aligning with the broader objectives of this study.

For the validation of the obtained classification, the finalized raster model encompassing seven classes was integrated with a comprehensive stack of vegetation indices and their respective statistical metrics, amounting to a total of 292 indicators. Numerous statistical metrics and tests are available for the assessment of the statistical significance pertaining to the differences observed between groups. In conducting a rigorous and thorough analysis to evaluate the statistical significance of various vegetation indices across distinct vegetation classes, we employed a variety of statistical techniques to ascertain the robustness and validity of our findings (Table S2-S8, Supplementary Materials)

The dataset, characterized by a diverse range of vegetation classes, was subjected to a stringent preprocessing stage. The filtration of outliers was deemed paramount to ensure the accuracy and reliability of the results. A meticulous categorization of the vegetation indices was carried out into six specific groups: Reflectance-based indices, Water content indices, Soil-adjusted indices, Chlorophyll indices, Non-linear indices, and Structural indices. This categorization significantly facilitated a nuanced and contextually relevant analysis, thereby enabling a more precise evaluation of the indices.

Utilizing the Kruskal-Wallis test, a non-parametric method, significant variances in the median values of the vegetation indices across different vegetation classes were identified. This preliminary step proved instrumental in pinpointing specific indices that merited further investigation. For a more detailed examination, a selection of statistical tests was employed, namely the t-test for comparing the mean values of two groups; the Kruskal-Wallis Test as a non-parametric alternative to ANOVA for three or more groups; Cohen's d effect size to gauge the magnitude of differences between two groups; and the Dunn Test as a post-hoc analysis following the Kruskal-Wallis Test. To account for Type I error during the tests, significance levels were adjusted using the Bonferroni correction. The findings of the statistical analysis are documented in the Table S9 at Supplementary Materials.

The meticulous approach towards statistical analysis entailed in this study was pivotal in delving deeper into the dynamics of vegetation indices across distinctive classes. The robust preprocessing of the dataset ensured a solid foundation for the subsequent statistical assessments. The strategic categorization of the vegetation indices provided a structured framework for a focused analysis, enhancing the interpretability and relevance of the findings.

Furthermore, the employment of both parametric (t-test) and non-parametric (Kruskal-Wallis Test and Dunn Test) statistical tests facilitated a well-rounded comprehension of the distribution and relationships among the vegetation indices across different vegetation classes. The utilization of Cohen's d effect size furnished a clearer understanding of the practical significance of the observed differences, complementing the p-value analysis from the t-test and Kruskal-Wallis Test.

Lastly, the implementation of the Bonferroni correction in the statistical analysis underscores a meticulous approach towards controlling the false positive rate, thereby fortifying the reliability of the findings. The detailed exposition of the statistical analysis outcomes in the supplementary materials ensures transparency and provides a solid foundation for further evaluation and interpretation of the results, significantly contributing to the broader understanding of vegetation indices dynamics across varying vegetation classes.

Indices that demonstrated significant disparities were then subjected to post-hoc analysis using the Mann-Whitney U tests for pairwise comparisons between vegetation classes. This granular approach enabled us to identify the precise pairs of vegetation classes that exhibited significant differences, providing invaluable insights. Furthermore, we calculated the median differences for each vegetation index across the vegetation classes. This provided a clear, straightforward metric to gauge the extent of the differences, adding another layer of understanding to our analysis.

By considering the p-values from the Kruskal-Wallis test, the outcomes of the post-hoc tests, and the median differences, we could deduce the most representative vegetation indices for each category. This comprehensive and holistic methodology ensured that our conclusions were not just statistically robust but also practically relevant.

However, it is imperative to note that statistical significance does not always equate to practical or ecological significance. The results must be interpreted within the specific research context, considering the practical applications of the vegetation indices. The outlier filtration process was crucial in enhancing data quality, subsequently bolstering the reliability of our findings. Additionally, when interpreting the results from the post-hoc tests, it was essential to adjust for multiple comparisons to mitigate the risk of false positives.

Incorporating a discussion perspective, this analysis sheds light on the intricate dynamics of vegetation indices across different vegetation classes. The identification of the most representative indices for each category stands as a testament to the importance of a thorough and contextually aware statistical approach. Similar methodologies have found applications in various ecological and remote sensing research endeavors, underscoring the universality and relevance of our approach.

Evaluating holistically, it is noteworthy that for all community groups, across all vegetation indices and statistical metrics, the highest level of significance was achieved, registering at 0.0 or approaching it. This signifies the accuracy in the distribution of results from the unsupervised classification across community groups (Figure 6, Figures S1–S17 in supplementary materials).

Within the selected groups of vegetation indices, one most significant vegetation index was chosen from each group post-testing, and the variations in value distribution within community groups were analyzed, as shown in Figure 7.

The detailed analysis of vegetation indices across various forest formations unveils distinct patterns, reflecting the intricate interplay between vegetation characteristics and environmental conditions. Alnus (ALN), Fagus (FAG), and Pubescent Oak (PUB) formations consistently exhibit traits of robust and healthy vegetation, as evidenced by higher values in indices related to chlorophyll content, vegetation structure, and vitality. This can be attributed to their inherent biological characteristics and adaptation to their respective habitats, which might provide optimal conditions for growth and development.

In contrast, Brassica (BRA) and Picea (PIC) formations generally display lower values across the indices, suggesting less favorable conditions for vegetation. This could be a result of environmental stressors, less optimal soil conditions, or specific physiological adaptations of these species. The lower chlorophyll content and sparse vegetation structure might be indicative of these challenges.

Populus (POP) and OTHERS categories present a diverse range of conditions, as seen in the broad distribution across the indices. This diversity could be a reflection of the heterogeneity in species composition, age structures, and health conditions within these categories. The wide array of vegetation states captured in these categories underscores the complexity of forest ecosystems and the multitude of factors that influence vegetation dynamics.

Populus (POP), characterized by rapid growth and adaptability, tends to exhibit a broad range of values across the vegetation indices, indicating varying degrees of health and vitality within this formation. This variability could be attributed to the diverse age structures and health conditions of Populus species, as well as their ability to thrive in a variety of environmental conditions.

The OTHERS category, encompassing a mix of different species and vegetation types, showcases an even wider variability in the vegetation indices. This heterogeneity reflects the complexity of these mixed formations, where multiple factors including species diversity, interspecies interactions, and environmental conditions come into play, influencing the overall vegetation health and structure.

The observed patterns in the vegetation indices across different forest formations provide valuable insights into the health and vitality of the vegetation, as well as the environmental conditions prevailing in these habitats. These patterns are closely aligned with the known biological characteristics and ecological preferences of the dominant species in each formation, offering a nuanced understanding of forest ecosystem dynamics.

The robustness of Alnus, Fagus, and Pubescent Oak formations can be linked to their ecological adaptability and resilience, which enable them to maintain high levels of chlorophyll content and structural integrity even in challenging conditions. On the other hand, the lower index values observed in Brassica and Picea formations highlight the potential challenges faced by these species, possibly due to environmental stressors or less optimal growth conditions.

These insights underscore the importance of continuous monitoring and analysis of vegetation indices for effective forest management and conservation efforts. By leveraging the nuanced information provided by these indices, we can gain a deeper understanding of forest ecosystem dynamics, identify areas in need of conservation or restoration, and implement targeted interventions to promote forest health and resilience.

4. Conclusions

The results of the study showed that through unsupervised classification of multispectral satellite images, areas of vegetation cover are identified that correspond with high confidence to vegetation classes in the Braun-Blanquet system. Based on the obtained model and data on the sequestration potential of various types of forest communities, future studies are planned to determine the carbon balance of ecosystems of the Republic of Tatarstan. These data can also serve as a cartographic base for planning and implementing climate projects.