<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJE</journal-id><journal-title-group><journal-title>Open Journal of Ecology</journal-title></journal-title-group><issn pub-type="epub">2162-1985</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/oje.2022.128028</article-id><article-id pub-id-type="publisher-id">OJE-119140</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Assessment of the State of Forests Based on Joint Statistical Processing of Sentinel-2B Remote Sensing Data and the Data from Network of Ground-Based ICP-Forests Sample Plots
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Alexander</surname><given-names>S. Alekseev</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Dmitry</surname><given-names>M. Chernikhovskii</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Forest Inventory, Management and GIS, Saint-Petersburg State Forest Technical University, Saint-Petersburg, Russian Federation</addr-line></aff><pub-date pub-type="epub"><day>09</day><month>08</month><year>2022</year></pub-date><volume>12</volume><issue>08</issue><fpage>513</fpage><lpage>528</lpage><history><date date-type="received"><day>22,</day>	<month>June</month>	<year>2022</year></date><date date-type="rev-recd"><day>9,</day>	<month>August</month>	<year>2022</year>	</date><date date-type="accepted"><day>12,</day>	<month>August</month>	<year>2022</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The research was carried out on the territory of the Karelian Isthmus of the Leningrad Region using Sentinel-2B images and data from a network of ground sample plots. The ground sample plots are located in the studied territory mainly in a regular manner, laid and surveyed according to the ICP-Forests methodology with some additions. The total area of the sample plots is a small part of the entire study area. One of the objectives of the study was to determine the possibility of using the k-NN (nearest neighbor method) to assess the state of forests throughout the whole studied territory by joint statistical processing of data from ground sample plots and Sentinel-2B imagery. The data of the ground-based sample plots were divided into 2 equal parts, one for the application of the k-NN method, the second for checking the results of the method application. The systematic error in determining the mean damage class of the tree stands on sample plots by the k-NN method turned out to be zero, the random error is equal to one point. These results offer a possibility to determine the state of the forest in the entire study area. The second objective of the study was to examine the possibility of using the short-wave vegetation index (SWVI) to assess the state of forests. As a result, a close statistically reliable dependence of the average score of the state of plantations and the value of the SWVI index was established, which makes it possible to use the established relationship to determine the state of forests throughout the studied territory. The joint use and statistical processing of remotely sensed data and ground-based test areas by the two studied methods make it possible to assess the state of forests throughout the large studied area within the image. The results obtained can be used to monitor the state of forests in large areas and design appropriate forestry protective measures.
 
</p></abstract><kwd-group><kwd>Remote Sensing</kwd><kwd> Sentinel-2B Imagery</kwd><kwd> ICP-Forest Sample Plot</kwd><kwd> Tree Stand</kwd><kwd> Damage Class</kwd><kwd> k-NN (Nearest Neighbor Method)</kwd><kwd> Vegetation Index SWVI</kwd><kwd> Nonlinear Regression</kwd><kwd> Systematic Error</kwd><kwd> Random Error</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The state of forests has always been of great importance from the point of view of ecology, nature conservation, forestry and forestry industry due to their unique economic, environmental and social functions. In recent years, according to the ICP-Forests monitoring system [<xref ref-type="bibr" rid="scirp.119140-ref1">1</xref>], there has been a deterioration in the state of forests in Europe, for example, almost a third of the test areas of the system are moderately or severely damaged, and the dynamics over time of mean defoliation of tree crowns shows its increase, on the one hand, and a decrease in the number of completely healthy trees, on the other.</p><p>The main damaging forest factors are forest fires, insect attacks, pests and diseases, unfavorable weather conditions and different kinds of anthropogenic influences. It is shown that in Europe visible signs of tree damage were presented on almost half of the trees monitored in 2019. The most frequently reported biotic factors causing damage to trees were insects and drought became the main abiotic factor.</p><p>In Russia, forest fires are the main factor damaging forests, the areas of forests that they damage have increased substantially in recent years and amount to many millions of hectares [<xref ref-type="bibr" rid="scirp.119140-ref2">2</xref>].</p><p>Thus, it is confirmed that forest damage occurs over vast areas of millions of hectares. To counteract this, it is necessary to organize continuous monitoring of forests. The condition of forests is assessed as a result of their survey on networks of test areas, for example, according to the ICP-Forests program or by special expeditionary research. In both cases, the share of the surveyed area is a small part of the entire controlled territory, which makes it necessary to transfer information from the surveyed test areas or routes to the entire territory. Such a task must be solved with the help of remotely sensed data of the Earth, data from ground test areas and routes, as well as special statistical methods for processing combined information from these two sources.</p><p>The close to border forests health status is of large issue of both countries Russian Federation and Finland first due to the fact of possible chances for сross-border transport of pollutants, forest fires, pests and diseases. This territory is relatively dense populated. There is a well-developed road network from the both sides of the border. This area intensively used for different kinds of economic and social activities. Some of them such as forestry, recreation, tourism, sport, hunting, fishing, ecological farming, schooling and training are sensitive to forest health status. That is why authorities at regional and municipal levels, including forest authorities and nongovernmental groups of interests representing local people, are interested in the information on study area forests status. It is practical to use for the tree stands health status estimations on study area the level 1 regular grid of PSP established according ICP-Forests methodologies. At least two motives may also be referred to in want of such a decision: first, ICP-Forests program used to be applied for a long time in both countries, so about the identical grids of sample plots are already existed in each country, second, the methodology of field assessment of sample plots is exactly the identical and the results of all estimations and assessments will be comparable, which is very essential for between countries comparisons and conclusions.</p><p>The main objectives of this study were to test the utility of the k-NN method to assess the state of forests throughout the whole studied territory by joint statistical processing of data from ground sample plots and Sentinel-2B imagery, as well as examine the possibility of using the short-wave vegetation index (SWVI) calibrated on ground-based data to assess the state of forests of study area. In the course of the study, the idea of mandatory joint use of Earth remote sensing materials and ground-based sample areas was considered. To obtain reliable characteristics of the state of forests, it is necessary to calibrate and verify remotely sensed data based on field sample plots. Otherwise, the picture of the state of forests may be distorted.</p></sec><sec id="s2"><title>2. Background</title><p>The k-NN method of remotely sensed data classification occupies a special place among other methods for researchers in the field of forest management and forest monitoring. This method makes it possible to combine in one procedure the use of Earth remotely sensed data and data obtained from ground-based polygons. This possibility may be of interest for assessing the state of forests. The k-NN method for studying forests was first proposed and then put into practice by Finnish researchers [<xref ref-type="bibr" rid="scirp.119140-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref6">6</xref>]. In the early 2000s, many studies were published aimed at a comprehensive analysis of the possibilities and limitations of the k-NN method [<xref ref-type="bibr" rid="scirp.119140-ref7">7</xref>] - [<xref ref-type="bibr" rid="scirp.119140-ref12">12</xref>].</p><p>Part of the publications related to k-NN is devoted to the results of the experimental use of the method to solve the problems of the national forest inventory (NFI) in individual countries that differ in natural and economic conditions, forest management features, as well as tasks solved within the framework of the NFI [<xref ref-type="bibr" rid="scirp.119140-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref15">15</xref>]. Experiments on the use of the k-NN method were carried out in the NFI system of Sweden, USA, Norway, China, Germany and some other countries [<xref ref-type="bibr" rid="scirp.119140-ref16">16</xref>].</p><p>In recent decades, the study of the use of various vegetation indices for classifying and assessing the characteristics and health of vegetation, soils, and crops has become one of the most popular areas of remote methods. Additional reasons for the development of this direction are substantial technical progress associated with an increase in the number of satellites, improvement in the characteristics of remotely sensed materials, increasing their availability for numerous users, and continuous improvement of software and processing algorithms.</p><p>The basic foundations of the definition and interpretation of vegetation indices are set out in a number of textbooks and monographs [<xref ref-type="bibr" rid="scirp.119140-ref17">17</xref>] - [<xref ref-type="bibr" rid="scirp.119140-ref23">23</xref>]. Remotely sensed data processing programs, geoinformation systems and special applications can serve as tools for determining and studying vegetation indices based on remotely sensed materials [<xref ref-type="bibr" rid="scirp.119140-ref24">24</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref25">25</xref>].</p><p>Brief overviews of spectral characteristics of vegetation as properly as 27 vegetation indices which might also be determined by software ENVI are performed in [<xref ref-type="bibr" rid="scirp.119140-ref26">26</xref>]. It is mentioned that the predominant advantage of vegetation indices is the ease of obtaining them and a broad vary of tasks solved with their help. It is also noted that for particular natural conditions and number tasks, some indexes can provide more accurate predictions than others. The determination of an index that displays the investigated property as precisely as possible is carried out via comparing the results of index calculations with field data. It is additionally noted that any vegetation indices provide only relative estimates of the properties of vegetation cover, which can be interpreted and recalculated into absolute ones with the involvement of field data. The index values are influenced by using the characteristics of the sensor and shooting conditions [<xref ref-type="bibr" rid="scirp.119140-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref27">27</xref>].</p><p>A large number of studies are devoted to solving the problems of identifying and analyzing the relationship between vegetation indices and forest characteristics, as well as developing and improving forest monitoring methodology based on satellite data [<xref ref-type="bibr" rid="scirp.119140-ref28">28</xref>] - [<xref ref-type="bibr" rid="scirp.119140-ref33">33</xref>].</p></sec><sec id="s3"><title>3. Materials and Methods</title><p>The studies became carried out on the territory of the Vyborg district of the Leningrad region. The trial regions are organized as a regular network with a step of 4 - 8 km. Every trial area was implemented as permanent sample plot (PSP) and is a 4-element cluster North-South and West-East orientated, including 4 subplots on each of which became evaluated the state of 6 specially selected sample trees. The methodology of undertaking field work is defined in detail in [<xref ref-type="bibr" rid="scirp.119140-ref34">34</xref>] and effectively was extensively used in some of different special researches [<xref ref-type="bibr" rid="scirp.119140-ref35">35</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref36">36</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref37">37</xref>] [<xref ref-type="bibr" rid="scirp.119140-ref38">38</xref>].</p><p>A total of 99 sample plots were surveyed, during the field work and then included in the database for further processing. In <xref ref-type="table" rid="table1">Table 1</xref> data on the predominant tree species in the surveyed sample areas are given. Six species of woody plant became tested on PSP: Scots pine (Pinus sylvestris), Norway spruce (Picea abies), Birch (Betula pendula), Aspen (Populus tremula), Black alder (Alnus glutinosa) and Gray alder (Alnus incana). Four of them Scots pine, Norway spruce, Birch and Aspen occupied a dominant position on PSP even as two species Black and Grey alder have been presented through single trees.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Dominating tree species on sample plots located in Vyborg district of Leningrad region</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Tree species</th><th align="center" valign="middle" >Number of sample plots</th><th align="center" valign="middle" >Percent, %</th></tr></thead><tr><td align="center" valign="middle" >Scots pine</td><td align="center" valign="middle" >59</td><td align="center" valign="middle" >59.6</td></tr><tr><td align="center" valign="middle" >Norway spruce</td><td align="center" valign="middle" >28</td><td align="center" valign="middle" >28.3</td></tr><tr><td align="center" valign="middle" >Birch</td><td align="center" valign="middle" >11</td><td align="center" valign="middle" >11.1</td></tr><tr><td align="center" valign="middle" >Aspen</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1.0</td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" >99</td><td align="center" valign="middle" >100.0</td></tr></tbody></table></table-wrap><p>From the data in <xref ref-type="table" rid="table1">Table 1</xref> and <xref ref-type="fig" rid="fig1">Figure 1</xref>, it can be seen that Scots pine and Norway spruce predominate on the territory—they occupy 87.9% of the total number of sample plots. The health status of 2416 trees was assessed using the ICP-Forests methodology based on estimates of defoliation and discoloration of trees crowns assessments during field work. As a result, every sample tree become attributed via health status grade: 0—complete health tree, 1—slightly damaged, 2—moderately damaged, 3—severely damaged and 4—dead tree. Species distribution of the trees is provided in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p><xref ref-type="table" rid="table2">Table 2</xref> data demonstrate that in the total number of model trees evaluated in the study area in the Vyborgsky district of the Leningrad region, Scots pine and Norway spruce predominate.</p><p>The developed method of field observations at each PSP included the registration in a circle with a radius of 100 m from the central tree of a set of additional characteristics related to the tree stand health and vitality: occurrence of infrastructural elements, garbage dumps, footprints of cuttings, thinning, forest fires, wind withdrawals, snow and frost damage to trees, recreation footprints (fire places, campfire, walking trail) and other signs of forest damage. With the help of maps and satellite data, the proximity of the PSP to quarries and settlements was determined.</p><p>To characterize the forest health status, mean damage class was calculated for each sample plots (see <xref ref-type="table" rid="table3">Table 3</xref>). This variable is considered as a quantitative measure of forest health status at tree stands level and as initial data for spatial analysis in geographic information systems (GIS).</p><p>To be able to interpret <xref ref-type="table" rid="table3">Table 3</xref> data statistics the subsequent scale was used: full healthy tree stands were those with a mean damage class in the range of 0 - 0.5, slightly damaged 0.51 - 1.50, moderately 1.51 - 2.50, severely damaged 2.51 - 3.50, dying and dead 3.51 - 4.00. <xref ref-type="table" rid="table3">Table 3</xref> information says that stands of all tree species are either healthy (spruce, birch) or slightly damaged (pine), broadleaves are in better condition than coniferous ones, and pine stands are most damaged. In general, the state of stands in the border region of Russian Federation and Finland on the part of Russia, within the Vyborg district of the Leningrad region, may be assessed as good: 56 sample plots were assessed as completely healthy and 43 as slightly damaged. To increase the statistical representation of the tree stands with negative health status, forest allotments with severe</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Species distribution of the trees estimated on the grid of sample plots</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Tree species</th><th align="center" valign="middle" >Number of trees</th><th align="center" valign="middle" >Percent, %</th></tr></thead><tr><td align="center" valign="middle" >Scots pine</td><td align="center" valign="middle" >1329</td><td align="center" valign="middle" >55.0</td></tr><tr><td align="center" valign="middle" >Norway spruce</td><td align="center" valign="middle" >668</td><td align="center" valign="middle" >27.6</td></tr><tr><td align="center" valign="middle" >Birch</td><td align="center" valign="middle" >383</td><td align="center" valign="middle" >15.9</td></tr><tr><td align="center" valign="middle" >Aspen</td><td align="center" valign="middle" >31</td><td align="center" valign="middle" >1.3</td></tr><tr><td align="center" valign="middle" >Alder gray</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.1</td></tr><tr><td align="center" valign="middle" >Alder black</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.1</td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" >2416</td><td align="center" valign="middle" >100</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Mean damage class of the main forest species and groups of species on the surveyed sample areas</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Tree species or group of species</th><th align="center" valign="middle" >Mean damage class</th></tr></thead><tr><td align="center" valign="middle" >Scots pine</td><td align="center" valign="middle" >0.69</td></tr><tr><td align="center" valign="middle" >Norway spruce</td><td align="center" valign="middle" >0.40</td></tr><tr><td align="center" valign="middle" >Coniferous</td><td align="center" valign="middle" >0.59</td></tr><tr><td align="center" valign="middle" >Birch</td><td align="center" valign="middle" >0.28</td></tr><tr><td align="center" valign="middle" >Broadleaves</td><td align="center" valign="middle" >0.28</td></tr><tr><td align="center" valign="middle" >All species</td><td align="center" valign="middle" >0.54</td></tr></tbody></table></table-wrap><p>damaged, dying and dead tree stands were added to the ground based data. To solve this problem, forest management data from standard geoinformation databases of the North-Western forestry enterprise were used. With the help of WinGIS and PLP-2015 programs, a set of forest parcels was selected from databases, for which the degree of damage was marked as severely damaged, drying out, as well as dead stands. Then, a set of polygons with the boundaries of forest parcels with severely damaged, drying and dead stands was exported from the forest management databases to the GIS. The suitability of selected plots for further analysis was assessed using fine spatial resolution images from the public domain (Bing, Google and Yandex). In total, 153 ground sample areas were used during processing, which were divided into two equal groups. Training trial plots were used to classify the territory according to the degree of vegetation damage, and control plots were used to assess the quality of the classification.</p><sec id="s3_1"><title>3.1. k-NN Classification Method</title><p>A general description of the kNN method is as follows [<xref ref-type="bibr" rid="scirp.119140-ref7">7</xref>]. The spectral distance, d<sub>pi,p</sub> is computed in the feature space from the pixel p to be classified to each pixel p<sub>i</sub> for which the ground measurement or class is known. For each pixel p, take k-nearest field plot pixels (in the feature space) and denote the distances from the pixel p to the nearest field plot pixels by d p i , p , ⋯ , d p k , p ( d p i , p ≤ ⋯ ≤ d p k , p ). The estimate of the variable value for the pixel p is then expressed as a function</p><p>of the closest units. In case of numerical variable all of the nearest neighbor values are taken in consideration, but each of them has weighed proportionally to the distance from classifiable pixel. A commonly used function for weighting distances is:</p><p>w ( p i ) p = 1 d ( p i ) p t / ∑ j = 1 k 1 d ( p i ) p t , (1)</p><p>where: w ( p i ) p —vector of weighting values, d ( p i ) p t —feature space distance from classifiable (p) and sample data (p<sub>i</sub>), t = 2.</p><p>The estimate of the variable m for pixel p is then:</p><p>m ^ p = ∑ i = 1 k w p i , p m p i , (2)</p><p>where:</p><p>m ^ p —numerical value of the parameter to be classified for the classifiable pixel;</p><p>m p i —value of the parameter to be classified for nearest neighbor pixel p<sub>i</sub>.</p><p>Classification of the Sentinel image by the k-NN method according to the damage classes of tree stands was performed on the basis of the training part of the sample plots. To determine the quality of the classification, the results of the classification were compared with the data of the control sample areas. The following formulas were used to calculate systematic and random (standard deviation) errors:</p><p>Δ x = ∑ i = 1 n ( x r i − x g i ) n , (3)</p><p>here ∆x is a systematic error in determining the damage class of tree stands, points; x r i and x g i are the damage classes of tree stands obtained as a result of Sentinel image classification by the k-NN method and ground-based methods, respectively; n is the number of control sample plots.</p><p>SD = ∑ i = 1 n ( Δ x i − Δ x ) 2 n − 1 , (4)</p><p>here, SDis a random error in determining the state class of stands, points;</p><p>Δ x i = x r i − x g i is an error in determining the tree stand damage class by ground and remote methods for the control test sample plot i, points.</p><p>The Sentinel-2B satellite image (date of shooting July 31, 2018) was used to classify the state of the forest. A single multi-spectral image was formed from 4 spectral channels with a spatial resolution of 10 m: 8—NIR, 4—Red, 3—Green, 2—Blue. On <xref ref-type="fig" rid="fig2">Figure 2</xref> is shown the location of the ground sample plots in the study area.</p><p>Pre-processing of survey materials (atmospheric correction, multi-spectral image formation) was performed using GIS QGIS 3.10.11-A Coruna and the Semi-Automatic classification plugin [<xref ref-type="bibr" rid="scirp.119140-ref39">39</xref>]. To perform classification by the k-NN method was carried out using the QGIS Desktop 2.14.12 GIS with the k-NN classifier plugin [<xref ref-type="bibr" rid="scirp.119140-ref40">40</xref>]. The initial settings of the k-NN classifier plugin used for classification were: the number of nearest neighbor k = 5, the number of spectral channels—four, and the classification mode—for continuous values.</p></sec><sec id="s3_2"><title>3.2. Short-Wave Vegetation Index SWVI</title><p>In the process of the study, a set of known vegetation indices [<xref ref-type="bibr" rid="scirp.119140-ref26">26</xref>] was considered to assess the state of forest stands based on the Sentinel image. But substantial statistical dependencies between the values of the index and the damage classes of forest stands were obtained only for the short-wave vegetation index SWVI. The short-wave vegetation index SWVI is calculated from the NIR and SWIR bands:</p><p>SWVI = NIR − SWIR NIR + SWIR , (5)</p><p>here, NIR is the reflection of vegetation cover in the near-infrared part of the spectrum with a range of 0.76 - 0.9 microns; SWIR is the reflection of vegetation cover in the mid-infrared part of the spectrum with a range of 1.55 - 1.75 microns.</p><p>The Statgraphics 18-X64 software was used to find and study the relationship between the value of the short-wave vegetation index SWVI and the damage classes of tree stands. The possibilities of variance and regression analysis were used.</p></sec></sec><sec id="s4"><title>4. Results and Discussion</title><sec id="s4_1"><title>4.1. k-NN Classification Method</title><p>After classification of the Sentinel-2B image by k-NN method made on the base of training PSP calculations of systematic and random errors were done for testing PSP to estimate the classification quality. On <xref ref-type="fig" rid="fig3">Figure 3</xref> is shown one of the results of the k-NN classification—a thematic map, coloured by tree stand damage classes.</p><p>After calculating the errors, the following values were obtained: the systematic error was equal to zero (∆x = 0 points), and the random error was 1 point (SD = 1 point). These error values mean that there is no systematic error in the resulting classification, and as a result of a random error, the health status of forest stands can be mistakenly determined by 1 point better or worse compared to the actual value. Taking into account that the obtained error values can be assessed as not very substantial, and the use of the k-NN method can be considered suitable for classifying forest stands by damage classes.</p><p>increase the reliability of the results. This method can be used to solve many The proposed method for classifying the health state of forests is based on the joint use of remote sensing data and field sample plots, which can substantially</p><p>pressing problems of modern sustainable forest management, such as:</p><p>- planning activities for the use, reproduction, protection and conservation of forests based on their condition;</p><p>- identification and mapping of forests according to the degree of damage;</p><p>- development of strategies for forest management planning for damaged forests, taking into account preserving their stability and sustainability;</p><p>- identification and study of the causes of forest damage and planning of measures to mitigate and eliminate them.</p></sec><sec id="s4_2"><title>4.2. Short-Wave Vegetation Index SWVI</title><p>The results of the analysis of variance showed a significant relationship between the damage classes of forest stands and the vegetation index SWVI (see <xref ref-type="table" rid="table4">Table 4</xref>).</p><p>The F-ratio (Fisher criteria), reflecting the level a ratio of the between-group estimate to the within-group estimate is 44.53. Since the P-value of the F-test is less than 0.05, there is a statistically significant difference between the mean vegetation indexes SWVI from one level of damage class to another at the 95% significance level. Based on the fact that the P-value of the F-test is less than 0.05, it can be argued that there is a statistically significant difference in the mean values of the SWVI vegetation indices between damage classes at a 95% significance level. Graphical analysis of the average values of damage classes shows a significant (more than 4 times) difference between the average index value from 0.08 for dead stands to 0.34 for healthy stands, which indicates the sensitivity of the SWVI index to the state of forest stands. But the graph also shows that the values of the SWVI index for some classes practically do not differ from each other. For</p><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> ANOVA table for vegetation index SWVI by 5 damage classes</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Variability source</th><th align="center" valign="middle" >Sum of Squares</th><th align="center" valign="middle" >Degree of freedom</th><th align="center" valign="middle" >Mean square</th><th align="center" valign="middle" >F-ratio</th><th align="center" valign="middle" >P-value</th></tr></thead><tr><td align="center" valign="middle" >Between damage classes</td><td align="center" valign="middle" >1.63300</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.408251</td><td align="center" valign="middle" >44.53</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >Within damage classes</td><td align="center" valign="middle" >1.35677</td><td align="center" valign="middle" >148</td><td align="center" valign="middle" >0.009167</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" >2.98977</td><td align="center" valign="middle" >152</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr></tbody></table></table-wrap><p>example, healthy and slightly damaged stands (see <xref ref-type="fig" rid="fig4">Figure 4</xref>).</p><p>The relationship between mean trees stands damage classes and vegetation index SWVI may be approximate by invert logistic curve:</p><p>Mean damage class = 4 / ( 1 + a ∗ exp ( − b ∗ SWVI ) )</p><p>where a and b are curve parameters. The shape of the curve statistically reliable describing the relationship between vegetation index SWVI and tree stands mean damage classes is shown on <xref ref-type="fig" rid="fig5">Figure 5</xref> and illustrated the uncertainty in determining neighboring damage classes by the value of index.</p><p>The large coefficient of determination of the curve equal to 83.5% once again confirms the close relationship between the vegetation index SWVI and the average class of damage to stands on PSP.</p><p>In order to increase the resolution of damage class determination it was decided to build new aggregated classes: the first aggregated class included healthy and slightly damaged stands, the second—moderately damaged and severely damaged stands, and the third—dead stands. The results of the analysis of variance for the aggregated classes of the vegetation index SWVI are presented in <xref ref-type="table" rid="table5">Table 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref>.</p><p>The results of the analysis of variance for the aggregated classes demonstrated a significant increase in the statistical relationship between the SWVI index and the aggregated damage classes. According to the performed analysis of variance for vegetation index SWVI by 3 aggregated damage classes, the F-ratio (Fisher criteria) equals 82.56, that almost twice more than in <xref ref-type="table" rid="table4">Table 4</xref> (44.53). <xref ref-type="fig" rid="fig6">Figure 6</xref> shows that the 95% confidence intervals for mean values of index SWVI calculated for 3 aggregated damage classes haven’t interceptions and can be used for damage class determination on higher stage of confidence.</p><p>The relationship for damage class calculation using the vegetation index SWVI as an explaining variable in this case of aggregated classes is linear with large coefficient of determination:</p><p>I ^ A = 3.753 − 8.453 ∗ SWVI , R 2 = 89.4 % ,</p><p>where I<sub>A</sub>—aggregated damage class. The resulting equation can be used to calculate the health state of forests and develop thematic maps based on remotely sensed data.</p><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> ANOVA table for vegetation index SWVI by 3 aggregated damage classes</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Variability source</th><th align="center" valign="middle" >Sum of Squares</th><th align="center" valign="middle" >Degree of freedom</th><th align="center" valign="middle" >Mean square</th><th align="center" valign="middle" >F-ratio</th><th align="center" valign="middle" >P-value</th></tr></thead><tr><td align="center" valign="middle" >Between damage classes</td><td align="center" valign="middle" >1.56663</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.7833</td><td align="center" valign="middle" >82.56</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >Within damage classes</td><td align="center" valign="middle" >1.42314</td><td align="center" valign="middle" >148</td><td align="center" valign="middle" >0.0095</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" >2.98977</td><td align="center" valign="middle" >152</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr></tbody></table></table-wrap></sec></sec><sec id="s5"><title>5. Conclusion</title><p>The k-NN approach of forest classification based totally on the joint use of remotely sensed data and ground-based sample areas has been broadly used in many countries for national forest inventories, mainly in Finland. With its help, objective data for large areas have been acquired during the forest inventory and the method demonstrated its practical effectiveness. In this study the k-NN approach was used to predict the condition and extent of damage to tree stands, which are also of great significance for solving some specific problems of sustainable forest management. The results of the application of developed methodology to assess the state of forests over large areas based on the k-NN approach were positive. A very important result is that the approach does not give systematic errors, and a random error creates uncertainty only when determining neighboring classes of damage to tree stands, which is not critical for its application. This method makes it possible to extrapolate data on the health status of tree stands on PSP, which in the geographical sense are points, to the entire territory represented through the satellite image. Applicability of the k-NN approach for estimation of forest health in assembles with ICP-Forest ground PSP is of significance due to the fact this network covers all the Europe and similar monitoring grids exist in the number of countries.</p><p>Some other end results of this study consist in revealed high sensitivity of short-wave vegetation index SWVI on tree stand health state as reflected by using its mean damage class. The value of SWVI varies as much as 4 times for healthy and dead tree stands. The higher the index value, the better the health status of plantings. If in analysis use the 3 aggregated damage classes instead of 5 preliminaries usually used the resolution ability of vegetation index SWVI elevates as much as twice and relationship between index value and aggregated damage classes became linear with large determination coefficient. This feature may be successively used for thematic mapping by means of GIS technology for visualization and evaluation of forest health state at the regions of interest. The practical application of this approach may be related to monitoring the condition of forests in remote areas. Thematic maps that replicate the state of plantings may be interesting to various fascinated organizations and agencies of citizens, such as state forestry management bodies, environmental organizations, hunters, fishermen and different recreational workers.</p><p>Generally, the methods primarily based on k-NN classification and use of vegetation index SWVI can be recommended for solving many issues of forestry related to monitoring the state of forests. Mainly, it could be used to organize nonstop monitoring of the health status of forests on the basis of acquiring updated materials of remote sensing for constantly informing involved companies of the population about the consequences.</p></sec><sec id="s6"><title>Acknowledgements</title><p>This study was carried out with the financial support of the project No 21-16-00065: “The role of insects and pathogens in the weakening and death of coniferous stands in the North-West of the Russian Federation: quantification and monitoring” funded by the Russian Science Foundation.</p></sec><sec id="s7"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s8"><title>Cite this paper</title><p>Alekseev, A.S. and Chernikhovskii, D.M. (2022) Assessment of the State of Forests Based on Joint Statistical Processing of Sentinel-2B Remote Sensing Data and the Data from Network of Ground-Based ICP-Forests Sample Plots. Open Journal of Ecology, 12, 513-528. https://doi.org/10.4236/oje.2022.128028</p></sec></body><back><ref-list><title>References</title><ref id="scirp.119140-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">ICP-Forests (2021) Tree Health Is Deteriorating in the European Forests. https://icp-forests.org/pdf/ICPForestsBriefNo5.pdf</mixed-citation></ref><ref id="scirp.119140-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Filipchuk, A.N., Malysheva, N.V., Zolina, T.A., Yugov, A.N. and Mironov, R.Y. (2021) Assessing the Possible Risks of Including the Reserved Forests in National Reporting under the UN Convention on Climate Change. Forestry Information, 2, 90-105.</mixed-citation></ref><ref id="scirp.119140-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Haapanen, R. and Ek, A. (2001) Software and Instructions for kNN Applications in Forest Resources. Description and Estimation. Staff Paper Series 152, University of Minnesota, Minneapolis. https://conservancy.umn.edu/handle/11299/37198</mixed-citation></ref><ref id="scirp.119140-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Tomppo, E., Czaplewski, R. and M&amp;#228;kisara K. (2002) The Role of Remote Sensing in Global Forest Assessment. Working Paper No. 61, Forest Resources Assessment, Rome. http://www.fao.org/docrep/006/ad650e/AD650E00.htm#TopOfPage</mixed-citation></ref><ref id="scirp.119140-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Kangas, A. and Maltamo, M. (2006) Forest Inventory: Methodology and Applications. Springer, Dordrecht. https://doi.org/10.1007/1-4020-4381-3</mixed-citation></ref><ref id="scirp.119140-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Tomppo, E., Haakana, M., Katila, M. and Per&amp;#228;saari, J. (2008) Multi-Source National Forest Inventory: Methods and Applications. Springer Science &amp; Business Media, Berlin.</mixed-citation></ref><ref id="scirp.119140-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Franco-Lopez, H., Ek, A.R. and Bauer, M.E. (2001) Estimation and Mapping of Forest Stand Density, Volume, and Cover Type Using the K-Nearest Neighbor’s Method. Remote Sensing of Environment, 77, 251-274. https://doi.org/10.1016/S0034-4257(01)00209-7</mixed-citation></ref><ref id="scirp.119140-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Haapanen, R., Ek, A.R., Bauer, M.E. and Finley, A.O. (2004) Delineation of Forest/Non Forest Land Use Classes Using Nearest Neighbor Methods. Remote Sensing of Environment, 89, 265-271. https://doi.org/10.1016/j.rse.2003.10.002</mixed-citation></ref><ref id="scirp.119140-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Koukal, T., Suppan, F. and Schneider, W. (2005) The Impact of Relative Radiometric Calibration on the Accuracy of kNN-Predictions of Forest Attributes. FORESTSAT 2005: Operational Tools in Forestry Using Remote Sensing Techniques, Bor&amp;#229;s, 31 May-3 June 2005, 17-21. https://doi.org/10.1016/j.rse.2006.08.016</mixed-citation></ref><ref id="scirp.119140-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Koukal, T., Suppan, F. and Schneider, W. (2007) The Impact of Relative Radiometric Calibration on the Accuracy of kNN-Predictions of Forest Attributes. Remote Sensing of Environment, 110, 431-437.</mixed-citation></ref><ref id="scirp.119140-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Gjertsen, A. (2007) Accuracy of Forest Mapping Based on Landsat TM Data and a kNN-Based Method. Remote Sensing of Environment, 110, 420-430. https://doi.org/10.1016/j.rse.2006.08.018</mixed-citation></ref><ref id="scirp.119140-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Meng, Q., Cieszewski, C.J., Madden, M. and Borders, B. (2007) K Nearest Neighbor Method for Forest Inventory Using Remote Sensing Data. GIScience &amp; Remote Sensing, 44, 149-165. https://doi.org/10.2747/1548-1603.44.2.149</mixed-citation></ref><ref id="scirp.119140-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">McInerney, D., Pekkarinen, A. and Haakana M. (2005) Combining Landsat ETM+ with Field Data for Ireland’s National Forest Inventory—A Pilot Study for Co. Clare. FORESTSAT 2005: Operational Tools in Forestry Using Remote Sensing Techniques, Bor&amp;#229;s, 31 May-3 June 2005, 12-16.</mixed-citation></ref><ref id="scirp.119140-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Beaudoin, A., Bernier, P.Y., Guindon, L., Villemaire, P., Guo, X.J., Stinson, G., Bergeron, T., Magnussen, S. and Hall, R.J. (2014) Mapping Attributes of Canada’s Forests at Moderate Resolution through kNN and MODIS Imagery. Canadian Journal of Forest Research, 44, 521-532. https://doi.org/10.1139/cjfr-2013-0401</mixed-citation></ref><ref id="scirp.119140-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Chernikhovskii, D.M. and Alekseev, A.S. (2019). The Method for Determining Forest Characteristics Based on Earth Remote Sensing Materials, Forest Management Data and the k-NN Algorithm (Case Study of Lodeynopol’skoe Forest District of Leningrad Region). Lesnoy Zhurnal, 4, 45-65. https://doi.org/10.17238/issn0536-1036.2019.4.45</mixed-citation></ref><ref id="scirp.119140-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">McRoberts, R.E. and Tomppo, E.O. (2007) Remote Sensing Support for National forest Inventories. Remote Sensing of Environment, 110, 412-419. https://doi.org/10.1016/j.rse.2006.09.034</mixed-citation></ref><ref id="scirp.119140-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Franklin, S. (2001) Remote Sensing for Sustainable Forest Management. CRC Press, Boca Raton. https://doi.org/10.1201/9781420032857</mixed-citation></ref><ref id="scirp.119140-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Schowengerdt, R.A. (2007) Remote Sensing: Models and Methods for Image Processing. 3rd Edition, Elsevier, Amsterdam.</mixed-citation></ref><ref id="scirp.119140-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Malysheva, N. (2012) Automated Interpretation of Aerospace Imageries of Forest Stands. Moscow State Forest University Press, Moscow.</mixed-citation></ref><ref id="scirp.119140-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Lillesand, T. (2015) Remote Sensing and Image Interpretation. 7th Edition, Wiley, Hoboken.</mixed-citation></ref><ref id="scirp.119140-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Yifang, B. (2016) Multitemporal Remote Sensing: Methods and Applications, Springer International Publishing, Cham.</mixed-citation></ref><ref id="scirp.119140-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Kurbanov, E., Vorobiev, O. and Lezhnin, S. (2015) Thematic Mapping of Vegetation Cover from Satellite Imagery: Validation and Accuracy Assessment. Volga State University of Technology, Yoshkar-Ola.</mixed-citation></ref><ref id="scirp.119140-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Bartalev, S., Egorov, V., Zharko, V., Loupian, E., Plotnikov, D., Khvostikov, S. and Shabanov, N. (2016) Land Cover Mapping over Russia Using Earth Observation Data. Russian Academy of Sciences’ Space Research Institute, Moscow.</mixed-citation></ref><ref id="scirp.119140-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Fisher, R., Hobgen, S., Mandaya, I., Kaho, N. R. and Zulkarnain, R. (2017) Satellite Image Analysis and Terrain Modeling—A Practical Manual for Natural Resource Management, Disaster Risk and Development Planning Using Free Geospatial Data and Software. Version 2. SAGA GIS 4, 150. https://sagatutorials.wordpress.com/</mixed-citation></ref><ref id="scirp.119140-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Eklunhd, L. and J&amp;#246;nsson, P. (2017) Timesat 3.3 with Seasonal Trend Decomposition and Parallel Processing: Software Manual. http://web.nateko.lu.se/timesat/docs/TIMESAT33_SoftwareManual.pdf</mixed-citation></ref><ref id="scirp.119140-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Cherepanov, A. and Druzhinina, E. (2009) Spectral Properties of Vegetation and Vegetation Indexes. Geomatics, 3, 28-32. https://sovzond.ru/upload/uf/082/0821bcc05ada9cb61aa16b9128bb09e9.pdf</mixed-citation></ref><ref id="scirp.119140-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Tokareva, O. (2010) Processing and Interpretation of Remote Sensing Data. Tomsk Polytechnic University Press, Tomsk.</mixed-citation></ref><ref id="scirp.119140-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Waring, R., Milner, K., Jolly, W., Phillips, L. and McWethy, D. (2006) Assessment of Site Index and Forest Growth Capacity across the Pacific and Inland Northwest U.S.A. with a MODIS Satellite-Derived Vegetation Index Forest. Ecology and Management, 228, 285-291. https://doi.org/10.1016/j.foreco.2006.03.019</mixed-citation></ref><ref id="scirp.119140-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Vorobyev, O., Kurbanov, E., Polevshikova, Y. and Leznin, S. (2016) Assessment of Dynamics and Disturbance of Forest cover in the Middle Povolzhje by Landsat Images. Current Problems in Remote Sensing of the Earth from Space, 13, 124-134. https://doi.org/10.21046/2070-7401-2016-13-3-124-134</mixed-citation></ref><ref id="scirp.119140-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Bartalev, S., Stytsenko, F., Khvostikov, S. and Loupian, E. (2017) Methodology of Post-Fire Tree Mortality Monitoring and Prediction Using Remote Sensing Data. Current Problems in Remote Sensing of the Earth from Space, 14, 176-193. https://doi.org/10.21046/2070-7401-2017-14-6-176-193</mixed-citation></ref><ref id="scirp.119140-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Vorobyev, O. and Kurbanov, E. (2017) Remote Monitoring of Vegetation Regeneration Dynamics on Burnt Areas of Mari Zavolzhje Forests. Current Problems in Remote Sensing of the Earth from Space, 14, 84-97. https://doi.org/10.21046/2070-7401-2017-14-2-84-97</mixed-citation></ref><ref id="scirp.119140-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Loranty, M., Davydov, S., Kropp, H., Alexander, H., Mack, M., Natali, S. and Zimov, N. (2018) Vegetation Indices Do Not Capture Forest Cover Variation in Upland Siberian Larch Forests. Remote Sensing, 10, Article No. 1686. https://doi.org/10.3390/rs10111686</mixed-citation></ref><ref id="scirp.119140-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Belova, E. and Ershov, D. (2019) Using Landsat Time Series for Assessing Reforestation on Clear Cuts in Bryansk Region. Forest Science Issues, 2, 1-20. https://doi.org/10.31509/2658-607x-2019-2-4-1-20</mixed-citation></ref><ref id="scirp.119140-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">ICP-Forests Manual (2020) http://icp-forests.net/page/icp-forests-manual</mixed-citation></ref><ref id="scirp.119140-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Alekseev, A.S. (2003) Forest Ecosystems Monitoring. Text Book. 2nd Edition, Saint-Petersburg State Forest Technical Academy Publishing, Saint-Petersburg.</mixed-citation></ref><ref id="scirp.119140-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Alekseev, A.S., Vetrov, L.S., Gurjanov, M.O., Nikiforchin, I.V., Chernikhovsky, D.M. and Chernov, I.M. (2020) Analysis of the Tree Stands Health Status in the Near Border Area of Russia and Finland Based on the Regular Grid of Sample Plots and GIS Technologies. IOP Conference Series: Earth and Environmental Science, 507, Article ID: 012001. https://doi.org/10.1088/1755-1315/507/1/012001</mixed-citation></ref><ref id="scirp.119140-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Alekseev, A.S. and Chernikhovsky, D.M. (2021) Assessment of the Health Status of Tree Stands Based on Sentinel-2B Remote Sensing Materials and the Short-Wave Vegetation Index SWVI. IOP Conference Series: Earth and Environmental Science, 876, Article ID: 012003. https://doi.org/10.1088/1755-1315/876/1/012003</mixed-citation></ref><ref id="scirp.119140-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">Alekseev, A.S., Chernikhovsky, D.M., Vetrov, L.S., Gurjanov, M.O. and Nikiforchin, I.V. (2021) Determination of the State of Forests Based on a Regular Grid of Ground-Based Sample Plots and Sentinel-2B Satellite Imagery Using the k-NN (“Nearest Neighbor”) Method. IOP Conference Series: Earth and Environmental Science, 876, Article ID: 012002. https://doi.org/10.1088/1755-1315/876/1/012002</mixed-citation></ref><ref id="scirp.119140-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">Congedo, L. (2020) Semi-Automatic Classification Plugin. Documentation Release 7.0.0.1.</mixed-citation></ref><ref id="scirp.119140-ref40"><label>40</label><mixed-citation publication-type="other" xlink:type="simple">Supervised kNN Classifier (2017) https://github.com/m6ev/kNN/find/master</mixed-citation></ref></ref-list></back></article>