Orographic Factors as a Predictor of the Spread of the Siberian Silk Moth Outbreak in the Mountainous Southern Taiga Forests of Siberia

: This research is dedicated to solving an urgent problem associated with the large-scale destruction of taiga forests by Siberian silk moth ( Dendrolimus sibiricus ) outbreaks. The dynamics of the damage to dark coniferous forest stands induced by the Siberian silk moth outbreaks in mid-altitude mountains were studied. A hypothesis was formulated based on the fundamental inﬂuence of the orography on the phytophage’s dispersal within the landscape, along with the climate, which acts as a secondary predictor—a catalyst for outbreaks. The study was carried out using Landsat − 8 satellite imagery time-series (from 2018 to 2020). The data were veriﬁed using a ﬁeld forest pathological survey of the territory. An assessment of the defoliated forest area and damage association with the landscape was carried out using an Advanced Spaceborne Thermal Emission and Reﬂection Radiometer (ASTER) digital elevation model. The assessment was aimed to detail the forecast parameters for an outbreak development in mid-altitude mountains using the orographic features—altitude, terrain slope, and slope aspect. Early warnings of phytophagous insect outbreaks in mountain southern taiga should be focused on the permanent monitoring of dark coniferous stands of the mossy group of forest types, covering altitude levels from 400 to 600 m, located on gentle terrains and slopes of up to 15 degrees. The greatest vulnerability to phytophage impacts was characterized as areas located at altitudes from 400 to 600 m. The upper limit of D. sibiricus distribution was 900 m above sea level. The results obtained provide comprehensive information on the Siberian silk moth potential reserves within the study area with the possibility of extrapolation to similar territories. The data will make it possible to model pest outbreaks based on orography and improve the forest pathological monitoring methods at the regional level.

Most studies base outbreak forecasting on climatic factors. Earlier observations have repeatedly proved the confinement of insect pest outbreaks to periods of drought and prolonged exposure to high temperatures [4,16,25]. From year to year, climatic indicators vary, and, only in the case of a certain temperature and humidity regime, do they act as predictors of an outbreak onset. Climatic factors facilitate the appropriate conditions for an outbreak but do not determine its dynamics. At the same time, a specific orographic base, a platform stable in time and space, characterizes any natural landscape.
The course of an outbreak from its beginning and further chronology is directly related to the site's location relative to the terrain within the natural complex. Therefore, according to the hypothesis put forward, the orographic factor is one of the fundamental predictors of pests living and spreading. The orographic factor is expressed through orographic features-altitude, terrain slope, and slope aspect.
The D. sibiricus phytophage's permanent habitats (reserves) are located within the host-tree range, directly confined to certain relief elements. Thus, the location of reserves is not accidental due to the species' bioecology specificity. Knowledge of a stand's entomoresistance level (identification of stands with low and high resistance to the phytophagous insect) within a specific relief will enable researchers to observe the current state in detail. This is why it is essential to determine parameters for predicting the onset and further development of an outbreak based on the orography. The advantages of a differentiated approach to the organization of background monitoring include the possibility of the early detection of mass defoliation and taking timely measures to prevent an outbreak.
The aim of the present study focused on the parameterization of D. sibiricus outbreaks to solve the following issues: to reveal the association between an insect outbreak occurrence and certain elements of the relief and analyze the spatiotemporal dynamics of the Siberian silk moth outbreak pattern based on the orographic features of the territory; and to assess the prospects for the early detection of Siberian silk moth outbreaks using the known patterns in species ecology and remote sensing methods in mountainous southern taiga forests conditions.

Study Area
The object of the study presents the area (266 thousand hectares) of the recent D. sibiricus outbreak within dark coniferous southern mountain taiga forests of the East-Sayan mountains (Krasnoyarsk Krai) ( Figure 1). The study area is located in the taiga zone of the Sayan region for forest conservation and belongs to the zone of medium forest pathology threat.
Medium-altitude mountains represent the landscape of the study area. The study site belongs to the South Siberian mountain forest zone and the Altai-Sayan mountains coniferous forest region. The region's climate is sharply continental and influenced by humid western air masses in the summer and the Siberian anticyclone in the winter. The growing season lasts for 149-151 days. The coldest month is January, with the lowest temperature standing at minus 50 • C. The warmest month is July, with the highest temperature of 38.6 • C. The mean annual precipitation is about 527 mm. The prevailing winds blow east-west with a mean wind speed of 2.8-4.7 m/s [26].
In the study region, the D. sibiricus mostly feeds on Siberian fir (Abies sibirica L.) and Siberian pine (Pinus sibirica Du Tour). The regional dark coniferous formations are dominated by the mossy group, represented by a whole set of forest types [27]. About half (41.3%) of the plantations belong to bilberry and bergenia forest types, mainly represented by cedar-fir stands (74.5% of the total area of bilberry and bergenia forest types). The forb-mossy group of plantations accounts for 28.1%. Bilberry-mossy cedar and fir type stands are confined to podzols soils. Medium-altitude mountains represent the landscape of the study area. The study site belongs to the South Siberian mountain forest zone and the Altai-Sayan mountains coniferous forest region. The region's climate is sharply continental and influenced by humid western air masses in the summer and the Siberian anticyclone in the winter. The growing season lasts for 149-151 days. The coldest month is January, with the lowest temperature standing at minus 50 °C. The warmest month is July, with the highest temperature of 38.6 °C. The mean annual precipitation is about 527 mm. The prevailing winds blow east-west with a mean wind speed of 2.8-4.7 m/s [26].
In the study region, the D. sibiricus mostly feeds on Siberian fir (Abies sibirica L.) and Siberian pine (Pinus sibirica Du Tour). The regional dark coniferous formations are dominated by the mossy group, represented by a whole set of forest types [27]. About half (41.3%) of the plantations belong to bilberry and bergenia forest types, mainly represented by cedar-fir stands (74.5% of the total area of bilberry and bergenia forest types). The forb-mossy group of plantations accounts for 28.1%. Bilberry-mossy cedar and fir type stands are confined to podzols soils.
The Bergenia cedar type of forests grows on poorly developed stony soils and mainly on steep slopes. Fir forests are represented by reed grassy and mossy-shrub types of forests and confined to habitats with a favorable combination of heat and humidity. Mossy pine forests occur fragmentarily and locally. Mountain-taiga larch forests are not widespread and are located, like spruce forests, exclusively in hollows and river valleys. Mature and overmature forest stands prevail.
The mean relative density is 0.5-0.7. In general, the share of such stands of the total area is more than 54.6%. High-density stands are represented mainly by cedar. ΙΙΙ-V bonitet classes (94.6%) characterize the conditions of forest growth on the territory. High-bonitet plantations are represented by birch and aspen stands. More than 85.0% of the low-bonitet forests are cedar-fir stands of long-mossy and green-mossy types. High-bonitet birch and aspen stands of the forb group of forest types occupy some areas of the plantations. The Bergenia cedar type of forests grows on poorly developed stony soils and mainly on steep slopes. Fir forests are represented by reed grassy and mossy-shrub types of forests and confined to habitats with a favorable combination of heat and humidity. Mossy pine forests occur fragmentarily and locally. Mountain-taiga larch forests are not widespread and are located, like spruce forests, exclusively in hollows and river valleys. Mature and overmature forest stands prevail.
The mean relative density is 0.5-0.7. In general, the share of such stands of the total area is more than 54.6%. High-density stands are represented mainly by cedar. III-V bonitet classes (94.6%) characterize the conditions of forest growth on the territory. High-bonitet plantations are represented by birch and aspen stands. More than 85.0% of the low-bonitet forests are cedar-fir stands of long-mossy and green-mossy types. High-bonitet birch and aspen stands of the forb group of forest types occupy some areas of the plantations.

Satellite Data Collection and Analysis
The study was carried out by a retrospective analysis based on the geographic information system (GIS) combined with the Earth Observation (EO) data. A Landsat-8 30 m Ground Sample Distance (GSD) time-series obtained from the Earth Explorer data portal (https://earthexplorer.usgs.gov/) was used for the analysis. There are well-proven high-precision methods for classifying areas damaged by phytophagous insects using such images in modern studies [13,21,[28][29][30].
Observations on pest outbreak development were conducted from June 2018 to May 2020. Table 1 shows the number of scenes according to satellite images for a specific date, demonstrating the state of the forests before the outbreak and at different damage stages.
The research presents the scenes that most clearly demonstrate the D. sibiricus distribution spatiotemporal dynamics within the study area. The period from 2018 to 2019 was taken for the analysis. The period was limited to 2019, as the D. sibiricus was not functioning during overwintering. Therefore, the situation in the image from 13.05.2020 entirely coincides with the image from 9 September 2019. Figure 2 shows a fragment of the satellite image, including the area of damaged tree stands.  EO data processing was carried out in the R program (v4.0.2) [31]. This is a free software environment for statistical computing and graphics. Visual analysis of the EO data was carried out using the QGIS program (v3.10) (https://qgis.org).
At the preparatory stage, a visual analysis of the areas damaged by the Siberian silk moth was carried out. Landsat-8 satellite images for 29 September 2018, 30 July 2019, and 9 September 2019 were the primary data source. Processing level L2 images were obtained from the U.S (Geological Survey website: (https://earthexplorer.usgs.gov/)). Clouds and cloud shadows were removed from the imagery during data preparation. The same date's images were glued together and re-projected into the WGS84 geographic coordinate system for complete coverage of the study area.
For better visibility of the damage, the following combination of channels was used: red-channel 7; green-channel 5; and blue-channel 3. With this channel combination, areas damaged by the Siberian silk moth were displayed in crimson, dark-coniferous stands in dark green, and deciduous and light-coniferous stands in light green. This remote sensing method made it possible to identify stands damaged by the phytophage. A sharply dissected erosional medium-altitude mountain relief with steep and prolonged mountain slopes, which often end as a river-cut cliff, dominates the territory. The mean height above sea level ranges from 470 to 1200 m. The degree of the slopes varies from 10 to 25 degrees, often reaching 40 degrees [27]. In the northern part of the study area, deciduous and light coniferous stands occupy relief elements. Wetlands are formed in depressions in the surface, which excludes this territory from the area of potential D. sibiricus distribution. Cedar and fir stands occur at higher elevations and are concentrated in the site's central and southern parts. EO data processing was carried out in the R program (v4.0.2) [31]. This is a free software environment for statistical computing and graphics. Visual analysis of the EO data was carried out using the QGIS program (v3.10) (https://qgis.org).
At the preparatory stage, a visual analysis of the areas damaged by the Siberian silk moth was carried out. Landsat-8 satellite images for 29 September 2018, 30 July 2019, and 9 September 2019 were the primary data source. Processing level L2 images were obtained from the U.S (Geological Survey website: (https://earthexplorer.usgs.gov/)). Clouds and cloud shadows were removed from the imagery during data preparation. The same date's images were glued together and re-projected into the WGS84 geographic coordinate system for complete coverage of the study area.
For better visibility of the damage, the following combination of channels was used: red-channel 7; green-channel 5; and blue-channel 3. With this channel combination, areas damaged by the Siberian silk moth were displayed in crimson, dark-coniferous stands in dark green, and deciduous and light-coniferous stands in light green. This remote sensing method made it possible to identify stands damaged by the phytophage. Later, this correspondence was established during field surveys of the territory (see Section 2.3).
Two popular classification algorithms were used to create a mask for damaged stands: random forest [32] and xgBoost [33], implemented in the caret package [34] for the R language. The classification aimed to divide each image's pixels into two classes-stands damaged by D. sibiricus, and the rest. Areas of both types were indicated on the images based on EO data visual assessment. A sample was formed from the selected polygons, which were divided into training and testing (a proportion of 75/25, respectively). EO data channels 1 to 7 with a spatial resolution of 30 m were included in the classification. The parameters of the classifier were selected separately for the setting by maximizing each Kappa measure. Cross-validation (cv) dividing the training sample seven-fold was used to search for hyper parameters. The result was two xgBoost models and one random forest model. The classification accuracy is presented in the form of standard statistical indicators in Table 2. As a result of tuning the models, the following parameters were obtained: 29 September 2018 (xgboost: nrounds = 600, max_depth = 20, eta = 0.1, gamma = 0, colsample_bytree = 0.9, min_child_weight = 1, and subsample = 1); 30 July 2019 (xgboost: nrounds = 900, max_depth = 5, eta = 0.1, gamma = 0, colsample_bytree = 0.9, min_child_weight = 1, and subsample = 1); 9 September 2019 (random forest: mtry = 4 and ntree = 500).
For XGBoost: nrounds-number of trees to grow; max_depth-max depth of the tree; eta-control of the learning rate; gamma-control regularization (gamma = 0 means no regularization); colsample_bytree-control the number of features (variables) supplied to a tree; min_child_weight-if the leaf node has a minimum sum of instance weight (calculated by second-order partial derivative) lower than min_child_weight, the tree splitting stops; subsample-controls the number of samples (observations) supplied to a tree. For random forest: mtry-the number of variables is randomly collected to be sampled at each split time; ntree-the number of branches to grow after each time split.
The resulting masks of damaged stands were used to extract the orographic characteristics of the area. The ASTER GDEM (Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model) v3 product (https://lpdaac.usgs. gov/news/nasa-and-meti-release-aster-global-dem-version-3/) was used as the basis for obtaining the characteristics of the relief. The spatial resolution of the digital elevation Land 2021, 10, 115 7 of 16 model was 30 m; the coordinate system was WGS84. The slopes' exposure and degree were calculated using standard methods in the raster package [35] and the terrain function was included in it. Eight neighboring cells of the raster were used to calculate these parameters.
For slopes of more than 5 degrees, the exposure from the azimuth was transformed into the names of the cardinal points in the following intervals: The terrain characteristics, which included the altitude, terrain slope, and slope aspect, were extracted from the damage mask for each date and converted into a data table that was analyzed using standard descriptive statistics methods. Statistical indicators of the distribution density were measured by the number of pixels.
Before starting the analysis, pixels related to damage at an altitude of more than 1000 m were regarded as random classification errors and removed from the data. Such a methodological approach was based on the available scientific information and our conclusions obtained from EO data processing. According to the scientific data, in the Altai-Sayan ecoregion forests, the highest pest density is observed at lower elevations. As a population moves up, its density decreases. The D. sibiricus population in the study region reaches 800-900 m above sea level and is rarely encountered at higher altitudes [4]. According to other researchers [13,18,36,37], global warming changes pest's ecology, facilitating elevation range expansion.
Potentially, this phenomenon was also considered in the analysis. However, when processing EO data, a clear boundary of the outbreak distribution at an altitude of up to 1000 m was revealed. Damages to forest stands by D. sibiricus above this altitude were removed from the set of analyzed data as classification errors. In total, 93 pixels were removed, which was 5.08 ha (0.82 ha for 30 July 2019 and 4.26 ha for 9 September 2019). When determining the area of damaged forest, stands subjected to cutting and forest fires were excluded. The classification algorithm somewhat underestimates damaged areas because it does not measure small damage to plantations located mainly at the outbreak's edges.
The analysis considered only the increase in forest territory damaged by the D. sibiricus. For example, areas damaged on 30 July 2019, did not include the damage detected earlier.
The same is true for damage to 9 September 2019. Such an approach allows analyzing only the forest canopy damage identified between the current and previous dates.
The determination of specific features of the plant damage locations, considering the peculiarities of orography, was carried out using the digital elevation model ASTER Global DEM Version 3 (https://lpdaac.usgs.gov/news/nasa-and-meti-release-aster-global-demversion-3/). The following layers of relief were used in the classification: height above mean sea level, slope exposure, and degree of slopes.
For analysis of the slope exposure, each exposure's damaged areas are expressed as a percentage of the total area of coniferous stands potentially included in the damaged zone, growing at this exposure. This ensures consideration of the factor of uneven distribution of food supply across slope exposures. The host-plant range was identified based on classified MODIS (Moderate-resolution Imaging Spectroradiometer) images [38] (MCD12Q1 v006) for 2017 with a spatial resolution of 500 m. The classification was based on the IGBP (International Geosphere-Biosphere Program) algorithm [39]. Pixels classified as an area covered by trees and more than 60% dominated by evergreen conifers or larch (I and III classes) were considered the host-plant range. The total area occupied by conifer plantations was 144,353.7 hectares. In all cases, exposures were calculated for slopes with an angle more than or equal to 5 degrees.

On-Ground Studies
According to forest pathologists of the Krasnoyarsk Center of Forest protection, the population's transition to the outbreak stage took place in the spring of 2019. The outbreak was facilitated by insufficient water supply during the growing season of 2015 and June, July, and September 2016. Notably, the daily mean air temperature was high during this period. As a result, favorable climatic conditions provoked the appearance of females with The EO data was verified during field surveys to visually assess the extent of the damage to trees caused by the D. sibiricus larvae. Field studies were conducted during May-June. The assessment was carried out using the method of laying out sample plots. The experimental sites were located at the following coordinates: No. Caterpillars were collected by striking the trunks of randomly selected model trees (five pieces on each sample plot), and their number was counted. In general, defoliation of the surveyed plantations ranged from 5-100%. According to remote sensing, the highest tree crown defoliation level was observed in forest areas where on-ground studies indicate lower indices of the pest density (from a few individuals to 129 individuals per MT). At the time of observation, this fact indicates a possible pest migration to nearby territories.

Results
In mountainous conditions, the distribution of the D. sibiricus depends on the height of a slope. Figure 3 presents the spatial characteristics of the damaged plantations' location, depending on the height above sea level within the study area. May-June. The assessment was carried out using the method of laying out sample plots. The experimental sites were located at the following coordinates: No. Caterpillars were collected by striking the trunks of randomly selected model trees (five pieces on each sample plot), and their number was counted. In general, defoliation of the surveyed plantations ranged from 5-100%. According to remote sensing, the highest tree crown defoliation level was observed in forest areas where on-ground studies indicate lower indices of the pest density (from a few individuals to 129 individuals per MT). At the time of observation, this fact indicates a possible pest migration to nearby territories.

Results
In mountainous conditions, the distribution of the D. sibiricus depends on the height of a slope. Figure 3 presents the spatial characteristics of the damaged plantations' location, depending on the height above sea level within the study area.   The graph shows that the median values of the altitude of tree damage in 2018 were lower than in subsequent periods. The range of variation of damage depending on altitude increased over time. The distribution of random variables appears to be approximately normal. Table 3 gives descriptive statistics of the distributions. Initial areas of sharp pest population growth were located at an altitude from 480 to 600 m above sea level. In July of the next growing season, the area of damaged forest stands increased by 29 times. Outbreaks were observed in the range from 430 m and reached an altitude of 883 m above sea level in 99.9% of cases. By September, the area of damage increased by another 4.4 times. In the mountainous part, the D. sibiricus outbreak distribution was limited by foothills of the ridge's northeastern slope and formed a long line with individual areas of local damage in the flattened terrain ( Figure 2).
The maximum elevation where stands defoliated by the D. sibiricus were found was 884 m above sea level, making a proportion of less than 0.5% of the total damaged area. Consequently, when predicting the development of Siberian silk moth outbreak for a given landscape condition, this factor should cover extreme points of the range from 400 to 900 m. To detect the pest population's activation early, one should focus on monitoring areas of dark coniferous forests at altitudes from 400 to 600 m.
Next, at the level of a given terrain, the possibility of identifying and understanding the natural structure of outbreaks onset places (the D. sibiricus reserves) was considered. The graph shows that the median values of the altitude of tree damage in 2018 were lower than in subsequent periods. The range of variation of damage depending on altitude increased over time. The distribution of random variables appears to be approximately normal. Table 3 gives descriptive statistics of the distributions. Initial areas of sharp pest population growth were located at an altitude from 480 to 600 m above sea level. In July of the next growing season, the area of damaged forest stands increased by 29 times. Outbreaks were observed in the range from 430 m and reached an altitude of 883 m above sea level in 99.9% of cases. By September, the area of damage increased by another 4.4 times. In the mountainous part, the D. sibiricus outbreak distribution was limited by foothills of the ridge's northeastern slope and formed a long line with individual areas of local damage in the flattened terrain ( Figure 2).
The maximum elevation where stands defoliated by the D. sibiricus were found was 884 m above sea level, making a proportion of less than 0.5% of the total damaged area. Consequently, when predicting the development of Siberian silk moth outbreak for a given landscape condition, this factor should cover extreme points of the range from 400 to 900 m. To detect the pest population's activation early, one should focus on monitoring areas of dark coniferous forests at altitudes from 400 to 600 m.
Next, at the level of a given terrain, the possibility of identifying and understanding the natural structure of outbreaks onset places (the D. sibiricus reserves) was considered. Figure 5 presents the dynamics of the distribution of areas damaged by D. sibiricus, depending on the terrain slope.  At the initial stage of outbreak development, the graph shows the concentration of damage sites on slopes from 8 to 14 degrees. In subsequent periods, the graphic lines shift to the left, to less steep slopes, and the concentration shifted 4 degrees, which characterizes the outbreak development process. From 30 July 2019 to 9 September 2019, damage zones associated with the slope steepness retained their general dynamics. During this period, damaged plantations were increasing, and insects were moving from reserve zones.
Insect outbreaks drive competitive relations, which increase the migratory activity of individuals. Therefore, if there is a food supply in nearby areas, pests move to more steep slopes. The pest's concentration and migration areas were mainly concentrated on gentle terrains and slopes up to 15 degrees. During field surveys of the damaged territory, the pattern of larval feeding was revealed. In mid-altitude mountains, larvae fed on needles from trees on the tops and middle parts of the slopes, avoiding lower areas due to higher humidity and cold microclimatic conditions. The graphically obtained results were confirmed by statistical characteristics of the distribution of areas regarding this feature, presented in Table 4. Damaged stands are mostly located on terrains with a slope angle of about 10 degrees, both in terms of the mean and median values. In 75% of cases during all observa- At the initial stage of outbreak development, the graph shows the concentration of damage sites on slopes from 8 to 14 degrees. In subsequent periods, the graphic lines shift to the left, to less steep slopes, and the concentration shifted 4 degrees, which characterizes the outbreak development process. From 30 July 2019 to 9 September 2019, damage zones associated with the slope steepness retained their general dynamics. During this period, damaged plantations were increasing, and insects were moving from reserve zones.
Insect outbreaks drive competitive relations, which increase the migratory activity of individuals. Therefore, if there is a food supply in nearby areas, pests move to more steep slopes. The pest's concentration and migration areas were mainly concentrated on gentle terrains and slopes up to 15 degrees. During field surveys of the damaged territory, the pattern of larval feeding was revealed. In mid-altitude mountains, larvae fed on needles from trees on the tops and middle parts of the slopes, avoiding lower areas due to higher humidity and cold microclimatic conditions. The graphically obtained results were confirmed by statistical characteristics of the distribution of areas regarding this feature, presented in Table 4. Damaged stands are mostly located on terrains with a slope angle of about 10 degrees, both in terms of the mean and median values. In 75% of cases during all observation periods, the outbreak spread predominantly in stands located on slopes of up to 13 degrees, sporadically reaching steep slopes of 30 degrees. Figure 6 shows the damaged area dynamics according to slope aspect differences at the specified date. The area increase value is expressed as a percentage of the coniferous stands total area stands on a particular slope aspect. Table 5 represents the statistical indicators of these distributions. tion periods, the outbreak spread predominantly in stands located on slopes of up to 13 degrees, sporadically reaching steep slopes of 30 degrees. Figure 6 shows the damaged area dynamics according to slope aspect differences at the specified date. The area increase value is expressed as a percentage of the coniferous stands total area stands on a particular slope aspect. Table 5 represents the statistical indicators of these distributions. The results show that the largest damaged areas recorded at the beginning of the outbreak cycle were located at the northern, northeastern, northwestern, and eastern slopes. This trend continued in subsequent periods. Empirically, this fact is related to the experimental site's location on the Eastern Sayan northeastern macroslope. The silvicultural characteristics of the forest area, where remote sensing indicated early signs of tree crown defoliation by the phytophage, were analyzed ( Table 6). The analysis clarified the confinement pattern of various forest types to certain relief elements in places of the D. sibiricus outbreak onset.  Figure 6. Dynamics of the area increase (%) of damaged dark coniferous stands depending on the slope aspect at a specific date. The results show that the largest damaged areas recorded at the beginning of the outbreak cycle were located at the northern, northeastern, northwestern, and eastern slopes. This trend continued in subsequent periods. Empirically, this fact is related to the experimental site's location on the Eastern Sayan northeastern macroslope.
The silvicultural characteristics of the forest area, where remote sensing indicated early signs of tree crown defoliation by the phytophage, were analyzed ( Table 6). The analysis clarified the confinement pattern of various forest types to certain relief elements in places of the D. sibiricus outbreak onset. A detailed analysis of the initial outbreak showed that damaged stands were represented mainly by mossy fir forests, occupying the entire spectrum of exposures on gentle terrains and slopes up to 15 degrees. Bilberry cedar forests were widely represented on gentle slopes with northerly, westerly, and easterly aspects. Bilberry fir forests were mainly concentrated on south-facing slopes. We concluded that the listed forest areas are recommended to be attributed to territories of low entomoresistance, and thus they should be considered as reserves of the Siberian silk moth when predicting an outbreak.
Considering the outbreak development in general ( Figure 5), the north-facing slope dominance smoothed over time since, due to larvae migration, the damaged area included nearby dark coniferous forests of the forb group of forest types located on slopes with various aspects. Mostly, these were bergenia cedar forests on steep slopes and reed grassy fir forests growing on flattened terrains and south-facing slopes. Hypothetically, within the study area, these plantations are not areas of the outbreak onset. However, when predicting an outbreak's spatial development, these stands should be attributed to the zone of subsequent (secondary) damage, which will be caused by the larvae migration from the reserves.

Discussion
Theoretically, the mountain taiga landscape's orography features determine the D. sibiricus outbreak development's spatiotemporal dynamics mainly by two factors-the phytophagous insect host-plant range and temperature and humidity regime of specific areas. Long-term observations have proven that all D. sibiricus outbreaks are concentrated on podzolic soils of mountainous regions. This coincidence expresses the silk moth's confinement to elevated parts of the relief, mainly to mountain ranges and plateaus. On the other hand, this pattern is associated with the pests' food plant distribution on these soils [14]. At the same time, D. sibiricus prefers warm habitats [3,4,6], so, as a rule, the pests' concentration is confined to well-warmed places. A favorable temperature regime is a fundamental criterion for determining the potential range of this pest's spread [3,5,19].
Dark coniferous formations within the study area were concentrated at the median height of 687 m above sea level. The onset of outbreak development was also confined to this altitude range. The share of 99% of the stands was located at an altitude of 397-1289 m. Thus, the upper limit of the Siberian silk moth's food supply was much higher than 900 m (this is the maximum altitude at which damage was detected). Therefore, the limit of damage depending on altitude was not associated with the food supply. At the same time, according to Yu. P. Kondakov [4], the upper limit of the D. sibiricus outbreak distribution at absolute altitudes of 800-900 m in mountainous southern taiga forests was associated with limiting climatic factors-long-term indicators of the sum of active temperatures 1200 • C and moisture coefficient of 2.0 (according to V. S. Mezentsev) [40].
Analysis of the D. sibiricus outbreak spatial distribution should consider the patterns of its ecological association with particular landscape elements [15,16,18]. The natural structure of places where an outbreak begins is essential. Complex indicators of the relief that allow for the prediction of pest reserve locations are the slope angle and exposure, combining to form the most favorable microclimatic conditions for a pest outbreak.
There are two factors determining the ground surface temperature in the mountains. The first is the aspect of the slope, which determines the energy resources of the habitat. The second factor is the position on the slope, which regulates the soil moisture reserves. Most of the previous studies mentioned that south-facing slopes were favorable for the D. sibiricus reserves due to greater insolation and higher ground surface temperature [3,5,6,13,19,20]. Such areas of dark coniferous forests are characterized by the most drained and warm conditions, where trees are growing on loose, relatively fresh, and porous soddy-podzolic soils with no signs of gleying [27]. Earlier melting of snow cover and later autumn frosts than on other forms of relief are observed at these areas, which prolong the larvae foliagefeeding period.
The soil aeration and hydrothermal regime, as well as the thickness of the litter layer, determine the best conditions for the wintering of D. sibiricus larvae in the study area [3,19]. However, in the given orographic landscape conditions, the opposite situation was revealed. The onset of damage and further development of the outbreak took place at slopes with northerly aspects. Regarding the location of the study area, this fact can be explained by the humidity conditions. In such conditions, the slopes with different aspects were less contrasting in terms of the heat and moisture supply.
The vegetation type or groups of formations remain unchanged, and differences appeared at the level of productivity and groups of forest types [26,27]. Therefore, under a given area's conditions, the slope exposure cannot act as a direct reference when predicting Land 2021, 10, 115 13 of 16 potential phytophage habitats. It is recommended to take the relief, including a particular type of forest and slope steepness, as the basic criterion. The initial hosting objects were cedar and fir forest stands of the mossy group of forest types located on gentle terrains and slopes of up to 15 degrees.
The mossy forest floor's temperature and humidity regime created the most favorable conditions for the Siberian silk moth development. Siberian silk moth larvae prefer to overwinter under the moss litter at a depth of 20 cm from the surface [3,5]. The clarification and adjustment of these area boundaries, considering the slope exposure, can be carried out in frames of field surveys. D. sibiricus avoids forest areas on steep erosional slopes and thin soils as well as wet and cold habitats, which typically form in lower-relief areas [3,5,19].
In the study area, the bergenia type of forest grows on poorly developed stony soils and mainly on steep slopes, so these cannot be areas of an outbreak onset. A specific pattern of a lesser proportion of damage to stands located in depressions was revealed during field surveys. Fir forests of the forb group of forest types were also characterized by late damage. The obtained results support the previous information on the pest outbreak's occurrence and spread in mid-altitude mountains. Thus, the present study results can act as parameters in the development of forecast models for the study area conditions.
The probabilistic forecasting of the D. sibiricus outbreak distribution in a specific terrain depends on certain other factors, more often associated with the silvicultural characteristics of forest stands. Defoliated forest areas are represented as irregular peripheral circles girding the outbreak epicenter with vast quantities of dead trees. Insect pest outbreaks spread from the center to the periphery as a damped wave-the farther from the core, the less the intensity of damage to forest stands [4]. An admixture of deciduous species negatively affects the dispersal of D. sibiricus.
Large areas of birch stands near the outbreak (with a proportion of six or more units) stop the spread of the pests completely [5]. As for the study area, such territories within the range of host plants were located at the upper limit of the outbreak spread, where little damage was observed, mainly limited by areas of birch-aspen stands. In other cases, deciduous forests growing in dark coniferous forests occupied an insignificant territory. Therefore, they could not provide a mosaic structure of the host-tree range, which could be a significant obstacle to the pest's migration to nearby forests.
Stand density does not play a significant role in the pest's dispersal, except for lowdensity stands, since the distance between trees is vital for larvae crawling. The most intense damage is often observed in forest stands with a relative density of 0.5-0.7 [5]. As for the study area, damaged sites included stands with a relative density of up to 0.8.
Researchers have indicated a nutritional selectivity of the pests [5,14]. For instance, the Siberian spruce needles are less attractive for D. sibiricus [12]. The presence of spruce in the forest composition is an indicator of unfavorable conditions for pest development. According to A.S. Isaev, in heavily damaged stands, the share of spruce in the forest composition does not exceed one unit (10%); at high values, the damage is reduced [18]. P.P. Okunev revealed, as a characteristic feature of the study area, the least degree of damage to Scots pine compared to other conifers [14].
The influence of the factors mentioned above on the larvae migratory ability in mountainous terrain requires additional studies aimed to provide an integrated approach in forecasting and monitoring the outbreak development within a particular area.

Conclusions
D. sibiricus permanently inhabits forest ecosystems in mountainous southern taiga, posing a constant threat of damage. Identifying the behavioral ecology of the phytophagous insect population in mountainous terrain opens up opportunities for improving monitoring systems and applying a digital elevation model to predict an outbreak's spread. This type of monitoring allows for the timely implementation of active forest protection measures. This is also important for the remaining fir and cedar forests near the studied Land 2021, 10, 115 14 of 16 defoliated areas, which continue to retain valuable biodiversity, being of great ecological and economic importance.
The study revealed specific Siberian silk moth distribution patterns aimed at detailing the preliminary observation zones in mountainous terrain. Background monitoring locations should be concentrated on dark coniferous forests covering altitude levels from 400 to 600 m with a potential distribution area of up to 900 m. Outbreak areas (reserves) may be located on gentle terrains and slopes of 10 to 15 degrees, excluding wet habitats and depressions. The pests prefer mossy forest floor, which is equally spread on slopes of both northern and southern exposures in mountain southern taiga forests.
Thus, the observation spectrum should include all areas of this forest type considering the above parameters. Based on the previous observations, south-facing slopes are assumed to be priority areas for reserves. The study results can serve as the basis for the parameterization of discrete data for mid-altitude mountain conditions when modeling the forecast of D. sibiricus outbreaks with the further possibility of extrapolating the revealed patterns to similar territories.

Funding:
The research was carried out within the projects "Fundamentals of forest protection from entomo-and fittings pests in Siberia" (№ FEFE-2020-0014) within the framework of the state assignment, set out by the Ministry of Education and Science of the Russian Federation, for the implementation by the Scientific Laboratory of Forest Health.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.