Pollution-Induced Changes in the Composition of Atmospheric Deposition and Soil Waters in Coniferous Forests at the Northern Tree Line

: This study examines the dynamics of the composition of atmospheric precipitation and soil water in coniferous forests under the inﬂuence of atmospheric emissions from the Severonickel Copper–Nickel Smelter in Russia’s Murmansk region. We studied dwarf shrub-green moss spruce forests and lichen-shrub pine forests, the most common in the boreal zone. Our results showed a signiﬁcant intra-(below and between the crowns) and inter-biogeocenotic (spruce and pine forests) variation in the composition of atmospheric precipitation and soil water in forests exposed to air pollution. The concentrations of main pollutants in atmospheric fallout and soil water are tens (sulfates) and hundreds (heavy metals) times higher than in the background areas and typically higher below the crowns. The long-term dynamics (between 1999 and 2020) of the composition of atmospheric fallout and soil water in coniferous forests in the background areas and defoliating forests demonstrates a signiﬁcant increase in nickel concentrations in recent years. This may be due to an increase in nickel concentrations in aerosols propagating over considerable distances. In pollution-induced sparse forests, a trend was found toward a decrease in the concentration of pollutants, which may indicate a decrease in the fallout of pollutants in the composition of larger particles close to the smelter.


Introduction
Intensive industrial development has led to widespread environmental pollution and the relevance of studying and preserving forest ecosystems has increased. Forest ecosystems growing at their northern limit of distribution are exposed to long-term impacts from mining and smelting operations. One of the least studied components of forest biogeocenoses is atmospheric precipitation and soil water, the study of the composition of that makes it possible to obtain more information about the biogeochemical cycles of elements in forest ecosystems. Numerous studies of the chemical composition of atmospheric and soil waters [1][2][3][4][5][6][7][8][9][10][11] showed that the chemical composition of atmospheric precipitation and soil solution changes significantly after contact with the forest canopy, depending on the species composition of woody plants and the impact of air pollution. In turn, the amount and quality of precipitation and soil water has a significant impact on the state of forest ecosystems [1,12]. Between the tree crowns at the northern limit of forest distribution, the precipitation is predominantly atmospheric, while below the crowns autogenic (biogenic) precipitation dominates [13]. When exposed to industrial air pollution, a significant transformation of the composition of atmospheric precipitation and soil water is observed, which consists of consistent changes in acidity and concentrations of heavy metals, sulfates, and nutrients in forests at different stages of degradation.
The Murmansk region is the northernmost region of the European part of Russia and its entire area belongs to the Russian Arctic. For more than eighty years, the forest ecosystems of Murmansk region have been exposed to intense industrial air pollution, which spreads over long distances. The main source of pollution is the Severonickel Smelter (Kola Mining and Metallurgical Company JSC (Murmansk, Russia)), located in the central part of Murmansk region in the city of Monchegorsk. The pollutants that affect the functioning of forest biogeocenoses include acid-forming compounds of sulfur and nitrogen, as well as compounds of heavy metals (nickel, copper, cobalt, lead, cadmium, chromium) [1,14]. The decrease in atmospheric emissions by the smelter and the decrease in the air pollution load on ecosystems causes it to be possible to consider the dynamic trends in the composition of individual components of forest ecosystems [14,15]. To predict the dynamics of biogeochemical cycles of elements in forests and propose measures aimed at increasing their stability under the combined action of natural and anthropogenic factors, it is necessary to estimate the spatial variation in the composition and properties of atmospheric precipitation and soil water, taking into account long-term changes.
The purpose of this study is (1) to assess the chemical composition of atmospheric precipitation and soil water in the northern taiga forests of the Murmansk region at different stages of technogenic digression due to the influence of a copper-nickel smelter, taking into account intra-and inter-biogeocenotic variation based on long-term observations and (2) to assess the long-term dynamic chemical composition of atmospheric fallout and soil water.

Materials and Methods
The objects of this study were atmospheric fallout and soil water in dwarf shrub-green moss spruce forests and lichen-shrub pine forests. The field studies were carried out on 6 permanent monitoring sites (PMSs) of the Institute of Industrial Ecology Problems of the North at the KSC RAS in the Murmansk region, Russia (Figure 1). The PMSs were established in typical forest ecosystems of the region in automorphic landscape positions. The PMSs are within transport accessibility for regular surveys. The spread of pollutants from the smelter site is controlled by the height of the pipe stacks (110-160 m), topography, and the strength and direction of the prevailing winds [16]. The predominance of southern and southwestern winds in the region in winter and northern and northwestern winds in summer [17] contribute to the formation of a plume of pollution from the smelter's emissions. The field studies were carried out on 6 permanent monitoring sites representing various stages of pollution-induced degradation: background areas, defoliating forests, and pollution-induced sparse forests ( Figure 1). The PMSs are located along the pollution gradient of the Severonickel Smelter at various distances from the source of pollution: in pollution-induced sparse forests at 7-10 km, in defoliating forests at 28-31 km, and in the background (control) areas at more than 150 km [1]. The background site represents the regional background conditions in the study area and meets all the globally recommended criteria for control sites [18]. The soils in pine forests are illuvial-ferruginous podzols (Carbic Podzols, WRB) and in spruce forests are illuvial-humus podzols (Rustic Podzols, WRB) [19,20].
The snow sampling was carried out annually from 1991/95 to 2020 under and between the crowns in triplicate. The snow cores were extracted using a plastic collector for the entire snow depth to the soil surface annually before the start of snowmelt (usually in the first week of April).
The rainwater sampling was carried out annually from 1999 to 2020 using sediment collectors (5 under the crowns, 5 between the crowns), which are plastic pipes with a funnel. A plastic bag is placed inside the pipe and secured with a cap. To prevent the ingress of plant litter, insects, and other particles, the surface of the pipe is covered with a removable fine mesh created with a synthetic material before it is secured with a cap. The snow sampling was carried out annually from 1991/95 to 2020 under and between the crowns in triplicate. The snow cores were extracted using a plastic collector for the entire snow depth to the soil surface annually before the start of snowmelt (usually in the first week of April).
The rainwater sampling was carried out annually from 1999 to 2020 using sediment collectors (5 under the crowns, 5 between the crowns), which are plastic pipes with a funnel. A plastic bag is placed inside the pipe and secured with a cap. To prevent the ingress of plant litter, insects, and other particles, the surface of the pipe is covered with a removable fine mesh created with a synthetic material before it is secured with a cap.
The soil water sampling was carried out annually from 1991/95 to 2020 using gravitytype lysimeters designed by D. Derom [21] (6 per PMS), located at different depths in accordance with the soil horizons (organogenic A0, eluvial E + B, illuvial ВС/С), and taking into account the mosaic structure of the biogeocenosis (under the crowns and between the crowns) [17].
During the field sampling campaign, the monthly volume of rain and soil water accumulated in the sediment collectors and lysimeters was measured using plastic measuring utensils. The sampling of lysimetric waters and rainwaters was carried out monthly from May to October.
The water samples were filtered through a blue-ribbon paper filter. The pH was measured potentiometrically, metals (Ca, Mg, Ni, Cu) by atomic absorption spectrophotometry, SO4 2− by ion-exchange chromatography, C by chromatometry or permanganatometry, depending on the concentration. The acid-neutralizing capacity (ANC) was calculated as the difference between the sum of basic cations and anions of mineral acids in mmol/L.
The descriptive statistics and trends to assess the long-term dynamics of the composition of atmospheric precipitation and soil water were prepared in Microsoft Excel 2019. To compare the composition of atmospheric precipitation and soil water at different stages of pollution-induced degradation of Northern taiga forests and to compare the composition of under-crown and inter-crown precipitation, the Mann-Whitney test and the software Statistica 13.3 were used. The soil water sampling was carried out annually from 1991/95 to 2020 using gravitytype lysimeters designed by D. Derom [21] (6 per PMS), located at different depths in accordance with the soil horizons (organogenic A0, eluvial E + B, illuvial BC/C), and taking into account the mosaic structure of the biogeocenosis (under the crowns and between the crowns) [17].
During the field sampling campaign, the monthly volume of rain and soil water accumulated in the sediment collectors and lysimeters was measured using plastic measuring utensils. The sampling of lysimetric waters and rainwaters was carried out monthly from May to October.
The water samples were filtered through a blue-ribbon paper filter. The pH was measured potentiometrically, metals (Ca, Mg, Ni, Cu) by atomic absorption spectrophotometry, SO 4 2− by ion-exchange chromatography, C by chromatometry or permanganatometry, depending on the concentration. The acid-neutralizing capacity (ANC) was calculated as the difference between the sum of basic cations and anions of mineral acids in mmol/L.
The descriptive statistics and trends to assess the long-term dynamics of the composition of atmospheric precipitation and soil water were prepared in Microsoft Excel 2019. To compare the composition of atmospheric precipitation and soil water at different stages of pollution-induced degradation of Northern taiga forests and to compare the composition of under-crown and inter-crown precipitation, the Mann-Whitney test and the software Statistica 13.3 were used.

Composition of Atmospheric Precipitation in the Form of Snow
The winter atmospheric fallout coincides with the period of biological dormancy. In boreal forests, the snow cover is present between one hundred and two hundred days a year, which determines the significant role of precipitation in the form of snow in biogeochemical cycles. Currently, industrial air pollution has a great influence on the formation of the composition of precipitation in the Russian North. A general discussion of the composition of snow fallout in the forests exposed to industrial air pollution, as well as of the effects of the tree canopy on the spatial distribution of elements in atmospheric fallout, is provided in a number of studies [1,7, [22][23][24][25]. The study revealed a significant influence of forest canopies Intra-biogeocenotic variability. In pine and spruce forests at the stages of defoliating forest and pollution-induced sparse forest, the concentrations of Ca, Mg, C, Ni, Cu, and SO 4 2− in the snow water under the crowns are significantly, up to 3 times, higher (p < 0.05)) compared to between the crowns. In the background area, a higher content was found under the crowns for C, Ca, and SO 4 2− (p > 0.05) ( Table 1). The revealed pattern is associated with the leaching of these elements from tree tissues, which is due to frequent winter thaws on the Kola Peninsula [24,26]. The pH value of the snow under the crowns is lower than between the crowns at all surveyed PMSs. The acid-neutralizing capacity never dropped to negative values in the background areas and in defoliating pine and spruce forests. The fallout values of compounds of elements did not have pronounced intra-biogeocenotic differences, which is associated, on the one hand, with increased concentrations of elements under the crowns, on the other hand, with a higher amount of precipitation between the crowns. Effects of air pollution. Industrial air pollution has a significant impact on the formation of snow water composition. In defoliating pine and spruce forests, the concentrations of many elements differ from the background area. Under the crowns and between the crowns, there was a significant (p < 0.05) increase in the concentrations of SO 4 2− compounds of up to three times, Cu compounds up to four times, and Ni compounds up to seven times compared with the background area. The concentrations of the main nutrients (Ca and Mg) were higher, concentrations C and the pH lower compared to the background area, which is most pronounced under the crowns. In pollution-induced sparse spruce and pine forests, compared to the background area, the concentrations of the main pollutants were significantly higher: SO 4 2− up to 4 times, Cu up to 100 times, and Ni up to 130 times under and between the crowns. In addition, there was a decrease in C concentrations and an increase in the actual acidity and concentrations of basic cations in pollution-induced sparse spruce and pine forests compared to the background area, which can be explained both by the influence of air emissions [1] and by the leaching of these elements from the tree canopy. The ANC index under the crowns in spruce and pine forests, as well as between the crowns in spruce forests was significantly lower than in the background area. As in defoliating forests, an increase in pollutant concentrations compared to the background area indicates a significant impact of industrial pollution on forest ecosystems as the source of emissions becomes closer.
Fallout of elements. The fallout of all elements, except for C, in the spruce and pine defoliating forests and pollution-induced sparse forests was significantly higher compared to the background areas both between and under the crowns, which, under pollution conditions, can be associated with (1) an increase in the concentrations of the elements that constitute the emissions, (2) an increase in the concentrations of elements as a result of leaching from damaged needles, and (3) an increase in fallout due to the deterioration of the crowns (defoliation) and an increase in their permeability.
Inter-biogeocenotic variation. When comparing the concentrations and fallout of elements in the snow water in spruce and pine forests at different stages of degradation, significant differences were found only under the crowns of pollution-induced sparse forests, where the concentrations of Cu, Ni, and SO 4 2− in the snow water are higher in pollution-induced sparse forests of spruce than of pine.
Multi-year dynamics. In the period from 1991/99 to 2020, the concentrations of elements in the snow of the Northern taiga forests at different stages of pollution-induced degradation are characterized by significant variability.
Under and between the crowns in the background Northern taiga forests, a trend of a major increase in the concentrations of copper and nickel in 2016 and 2017 was found. No other distinct long-term trends were found.
In defoliating forests, the long-term dynamics of the concentrations of elements in snow demonstrates high variability both under and between the crowns. Between the crowns in defoliating forests, a decrease in sulfate concentrations was observed, from 1.1 to 0.87 mg/L (R 2 = 0.37) in spruce forests and from 0.87 to 0.62 mg/L (R 2 = 0.29) in pine forests. There was a decrease in copper concentrations under the crowns in spruce forests from 0.08 to 0.004 mg/L (R 2 = 0.36) and between the crowns from 0.02 to 0.003 mg/L (R 2 = 0.49). This may be due to the reduction in pollutant emissions into the atmosphere from the smelter.
In pollution-induced sparse forests, no clear trends of the concentrations of Ni, Cu, and SO 4 2− in winter precipitation were identified, with the exception of a decrease in the Cu content in snow between the crowns in pine-from 0.13 to 0.0004 mg/L (R 2 = 0.34). This may be due to the ongoing deposition of large particles near the smelter.

Composition of Atmospheric Fallout in the Form of Rain
Atmospheric precipitation in the form of rain plays an important role in the cycling of chemical elements and the functioning of forest ecosystems. The chemical composition of rainfall changes significantly after passing through the forest canopy. During this interaction, physicochemical reactions occur, leading to a change in water acidity and the concentrations of most of the elements contained in it. Currently, atmospheric precipitation chemistry is significantly influenced by air pollution and boreal forests act as a modifier of pollution flows [3,4,7,9,14,25,[27][28][29], etc. The studies show that in the rainwater near the smelter, the concentrations of SO 4 2− and the concentrations of the main cations (Ca and Mg) were several times higher than the background values. The content of heavy metals in the rain was higher by an order of magnitude and more than in the background area. It is noted that filtered gas and dust emissions and chimneys are the dominant source of rainwater pollution, the reduction in emissions by the plant itself has significantly reduced the amount of pollutants coming from pipes into rainwater.
Intra-biogeocenotic variability. The composition of rainwater, as well as meltwater, at all stages of degradation in spruce and pine forests is characterized by pronounced intra-biogeocenotic variability: the concentrations of Ca, Mg, Ni, C, Cu, and SO 4 2− in the rainwater under the crowns are significantly higher (p < 0.05) than between the crowns ( Table 2). In addition, in the rainwater of spruce and pine forests, the ANC index under the crowns are significantly higher than between the crowns and the acidity under the tree crowns also significantly increases. The elevated concentrations of elements under the crowns indicate leaching of element compounds from woody plants [1,5,17]. Effects of air pollution. Rainwater in defoliating spruce and pine forests, compared with the background areas, is characterized by a significant increase in pollutant concentrations: copper up to 17 times, Ni up to 43 times, and sulfates up to 5 times. This pattern is most pronounced under the crowns. The concentrations of Ca and Mg in defoliating spruce and pine forests, compared to the background areas, are significantly higher, both under the crowns (with the exception of magnesium concentrations in spruce forest, p > 0.05) and between the crowns, which is explained by the leaching of these elements from tree crowns due to acid rain. When exposed to air pollution, the tissues of woody plants are depleted in inactive elements, such as calcium, which accumulates in aging and often damaged needles [17,30,31]. The concentration of C is significantly reduced compared to the background area, which is due to crown defoliation.
The rainwater in pollution-induced sparse pine and spruce forests, under and between the crowns, differs compared to the background by a sharp increase in the concentrations of pollutants that constitute emissions: Cu up to 75 times, Ni up to 200 times, and SO 4 2− up to 6 times. Between the crowns in spruce and pine forests, the concentrations of calcium and magnesium in rainwater are higher compared to the background, which may be due to their high content in industrial emissions and dusting of areas not covered with vegetation (pollution-induced wastelands). In the rainwaters of sparse pine and spruce forests, a decrease in C concentrations is observed compared to the background area due to the degradation of the crowns of woody plants.
An increase in the acidity of rainwater between and under the crowns was found, which indicates a high level of industrial air pollution. Compared to the background areas, the indicator of the acid-neutralizing capacity was significantly (p < 0.05) lower in the rainwater of defoliating spruce forests and pollution-induced sparse pine and spruce forests under and between the crowns. This pattern is explained by a lower content of organic acids due to thinning of crowns, due to the premature fall of needles and a significant increase in the concentration of sulfates under the influence of technogenic air pollution [1].
Fallout of elements. In pine and spruce forests at all stages of degradation, the fallout of compounds of elements has clear intra-biogeocenotic differences: the fallout of elements with rain under the crowns is in most cases significantly higher than between the crowns. At the same time, between the crowns, the amount of atmospheric precipitation in spruce and pine forests in the background area, as well as in defoliating spruce forests, is up to two times higher than under the crowns. In defoliating pine forests and pollutioninduced sparse spruce and pine forests, the differences in precipitation levels are leveled, which is associated with the thinning of tree crowns due to the impact of pollution. In defoliating forests and pollution-induced sparse forests, the deposition of SO 4 2− , Cu, and Ni is much higher, both under and between the crowns compared to the background area. This can be explained by an increase in their concentrations and an increase in the volume of precipitation due to an increase in the transmission capacity of the tree crowns. Fallouts, as well as carbon concentrations in rain in defoliating forests and pollution-induced sparse forests, are significantly lower compared to the background area, with the exception of defoliating pine forests, here the differences are not significant (p > 0.05). Similar results were obtained in earlier studies, where it was found that in pollution-induced sparse forests, the concentration of carbon in atmospheric fallout was lower compared to the background values, which is associated with a sharp decrease in the mass of needles and a significant influx of sulfates, nitrates, and chlorides with emissions from the copper-nickel smelter [1].
Inter-biogeocenotic variability. In the background area, the concentrations Ca, Mg, C, Ni, Cu, and SO 4 2− between and under the crowns are significantly higher in spruce forests compared to pine forests. In the background area and in defoliating forests, the pH and ANC values are higher in pine forests than in spruce forests. In defoliating spruce forests, the concentrations of elements in rain under tree crowns are significantly higher than in pine forests, with the exception of C concentrations. In pollution-induced sparse forests, the pH and ANC values are higher in the pine forest compared to the spruce forest, the concentrations of Ca, Mg, Ni, Cu, and SO 4 2− , on the contrary, are higher in the spruce forest, both under the crowns and between the crowns of the trees.
In the background area and in defoliating forests, the volume of precipitation penetrating through the tree canopy is significantly higher in pine forests. There were no differences in the amount of precipitation in the pollution-induced sparse forests. The fallout of Ca, Mg, Ni, Cu, and SO 4 2− is significantly higher under the crowns in spruce forests, which is associated with higher concentrations in the rain. The elevated concentrations of element compounds in the rainfall under the crowns of spruce forests are associated with a larger surface area and more pronounced barrier functions of spruce tree crowns [1, 11,13].
Multi-year dynamics. The atmospheric fallout of element compounds with rain during the period from 1999 to 2020 at different succession stages of coniferous forests is characterized by high variability. In the background areas, during the period from 1999 to 2020, under and between the crowns of Northern taiga forests, a trend of increasing Ni concentrations was found, most pronounced in pine forests (Figure 2). This can be explained by an increase in the level of pollutants in industrial emissions in the form of aerosols propagating over considerable distances, which indicates an increase in pollution loads [32,33].
Inter-biogeocenotic variability. In the background area, the concentrations Ca, Mg, C, Ni, Cu, and SO4 2− between and under the crowns are significantly higher in spruce forests compared to pine forests. In the background area and in defoliating forests, the pH and ANC values are higher in pine forests than in spruce forests. In defoliating spruce forests, the concentrations of elements in rain under tree crowns are significantly higher than in pine forests, with the exception of С concentrations. In pollution-induced sparse forests, the pH and ANC values are higher in the pine forest compared to the spruce forest, the concentrations of Ca, Mg, Ni, Cu, and SO4 2− , on the contrary, are higher in the spruce forest, both under the crowns and between the crowns of the trees.
In the background area and in defoliating forests, the volume of precipitation penetrating through the tree canopy is significantly higher in pine forests. There were no differences in the amount of precipitation in the pollution-induced sparse forests. The fallout of Ca, Mg, Ni, Cu, and SO4 2− is significantly higher under the crowns in spruce forests, which is associated with higher concentrations in the rain. The elevated concentrations of element compounds in the rainfall under the crowns of spruce forests are associated with a larger surface area and more pronounced barrier functions of spruce tree crowns [1, 11,13].
Multi-year dynamics. The atmospheric fallout of element compounds with rain during the period from 1999 to 2020 at different succession stages of coniferous forests is characterized by high variability. In the background areas, during the period from 1999 to 2020, under and between the crowns of Northern taiga forests, a trend of increasing Ni concentrations was found, most pronounced in pine forests (Figure 2). This can be explained by an increase in the level of pollutants in industrial emissions in the form of aerosols propagating over considerable distances, which indicates an increase in pollution loads [32,33]. Under the crowns in defoliating pine forests, a trend was observed toward decreasing copper concentrations in the rainwater from 0.14 to 0.04 mg/L (R 2 = 0.32). No other significant trends were found in the long-term dynamics of the concentrations of main pollutants in the rainwater in defoliating forests and pollution-induced sparse forests.

Soil Water Composition
The soil is the most important component of the forest ecosystem and serves as an accumulator for most of the elements present in the biosphere. The soil solution is the most active phase of the soil. Since the soil water is in direct contact with plant roots, information on soil water composition is of high indicative value for assessing the possible Under the crowns in defoliating pine forests, a trend was observed toward decreasing copper concentrations in the rainwater from 0.14 to 0.04 mg/L (R 2 = 0.32). No other significant trends were found in the long-term dynamics of the concentrations of main pollutants in the rainwater in defoliating forests and pollution-induced sparse forests.

Soil Water Composition
The soil is the most important component of the forest ecosystem and serves as an accumulator for most of the elements present in the biosphere. The soil solution is the most active phase of the soil. Since the soil water is in direct contact with plant roots, information on soil water composition is of high indicative value for assessing the possible negative impacts on forests. The chemical composition of soil water enables it to be possible to obtain information on the availability of nutrients and compounds of elements potentially toxic to biota in soils [34][35][36]. According to studies, significant intraprofile variability can be observed in the soil waters of boreal forests, which is expressed in a decrease in the concentrations of all elements with depth and is explained by the presence of biogeochemical barriers. It is also possible to observe pronounced intra-biogeocentric variability; soil waters in under-crown spaces are more concentrated and acidic than in intercrown spaces, which is due to the formation of intense flows of acid-forming substances with crown and stem waters. The variability of the composition of soil waters depending on the type of biogeocenosis is noted: it is shown that in the soil waters of spruce forests, the concentrations of elements, as a rule, are higher than in pine forests.
When approaching the source of pollution, there is a violation of the functioning of all components of the biogeocenosis, in connection with this, the concentrations of nutrients in soil water and the ANC index decrease, and an increase in the concentrations of the compounds of the elements included in the emissions is observed. These changes are most noticeable in woody, especially in spruce, parcels due to the high sorption capacity of the crowns of coniferous trees.
Intra-biogeocenotic variability. In the background area, in lysimetric water samples from all soil horizons in spruce and pine forests, the concentrations of Ca, Mg, C, Ni, Cu, and SO 4 2− were significantly (p < 0.05) higher under the crowns than between the crowns (Tables 3-5), which is associated with the leaching of element compounds from plant tissues by atmospheric precipitation. An exception is the concentration of Cu in the water from the A0 soil horizon in pine forests, as well as in the water from the BC soil horizon in spruce and pine forests (p > 0.05).
The soil water acidity under spruce crowns in the background area is significantly lower (p < 0.05) than between the crowns. This is explained by the higher Ca concentration in crown rainwater and in the soils under spruce. In addition, the dense and low spruce crowns retain a significant amount of precipitation, which prevents the intensive removal of bases from the soil. The organogenic and mineral horizons of soils under spruce trees are richer in available calcium compounds, the content of which in the falling needles of spruce is much higher than in pine [36]. In pine forests, significant differences in water acidity in the under-crown and inter-crown spaces were not found.
The ANC index in the water samples from all soil horizons under and between the crowns in the background spruce and pine forests decreases with depth and has a positive value. A decrease in the ANC index with depth can be observed in spruce and pine forests under and between the crowns, which is associated with a decrease in the concentrations of low molecular weight aliphatic organic acids [37]. In spruce and pine forests, the ANC index in the soil water under the crowns is significantly higher than between the crowns, with the exception of soil waters from the illuvial horizons of pine forests.
The intra-biogeocenotic differences in defoliating pine and spruce forests show that the concentrations of Ca, Mg, Cu, Ni, and SO 4 2− in the soil water from all soil horizons under the crowns are significantly higher than between the crowns. It should be noted that in pine forests, the concentration of Cu in the soil water from the soil horizon A0 between the crowns is significantly (p < 0.05) higher than under the crowns. In contrast to the background area, the acidity of soil water in defoliating spruce forests is significantly lower (p < 0.05) between the crowns compared to under the crowns. Under conditions of industrial air pollution, this can be explained by an increase in the supply of acid-forming substances from the atmosphere, especially under the crowns, and a disruption in the functioning of the phytocenosis. In defoliating spruce forests, the ANC index in the soil water from all soil horizons between the crowns is significantly higher than under the crowns, which may be due to the intense flow of acid-forming substances with precipitation under spruce crowns. In defoliating forests, the acidity and composition of atmospheric fallout is controlled by two additional factors-the level of atmospheric pollution and the degree of damage to the tree canopy. Under the conditions of aerotechnogenic pollution, due to the presence of the tree canopy, the lower tiers of vegetation, the soils of the undercrown spaces, and the trees themselves are exposed to a more powerful technogenic impact than the vegetation and soils of the inter-crown spaces [27]. The ANC indicator has a positive value in the soil water of spruce and pine forests under the crowns and between the crowns. Table 3. Concentration of element compounds in the soil water from organic soil horizons (A0) in spruce and pine forests in 1995-2020, mg/L, ANC in mmol/L. In pollution-induced sparse forests, the concentrations of Ca, Mg, Cu, C, Ni, and SO 4 2− and the acidity are significantly higher in the soil water under the crowns compared between the crowns. In spruce forests, the ANC indicator has a negative value in the soil water both under and between the crowns.

Stages of
Effect of air pollution. In the soil water under the crowns of defoliating spruce and pine forests, the concentrations of the main pollutants are much higher in comparison with the background areas: sulfates up to 4 times, copper up to 30 times, and nickel up to 40 times, which is most pronounced under the crowns. In the soil waters from all soil horizons under and between the crowns of spruce forests, the acidity of lysimetric waters increases in comparison with the background areas. In pine forests, in comparison with the background areas, the acidity also increases in the soil water from A0 horizons under-the crowns and in the soil water from BC horizons, both under and between the crowns. However, in the soil water from the E + B horizons in pine forests, on the contrary, a decrease in actual acidity was observed compared to the background. This may be due to an increase in the concentration of bases associated with the composition of parent rocks. In soil under crown waters in defoliating spruce and pine forests, the value of the ANC index is significantly lower than in the background area, except for waters from the A0 horizons of spruce forest soils in inter-crown spaces [1,36]. Table 4. Concentration of element compounds in the soil water from eluvial soil horizons (E + B) in spruce and pine forests in 1995-2020, mg/L, ANC in mmol/L. In technogenic light forests, the concentrations of the main pollutants in soil waters are significantly higher than in the background area: up to 10 times for SO 4 2− , up to 100 times for Cu, and up to 500 times for Ni, which is more pronounced in the under-crown space. The pH value of soil water in the under-crown spaces of pine woodlands, under-crown and inter-crown spaces of spruce woodlands is lower compared to the background. In the spruce and pine light forests, in soil waters, with the exception of waters from the soil horizon A0 of the under-crown spaces of the spruce forest, an increase in the concentrations of Ca and Mg is observed compared to the background territory, which is explained by parent rocks rich in bases. In sparse forests, the value of the ANC index in the soil water under and between the crowns is lower than in the background area. There is a decrease in C concentrations in waters from all soil horizons in defoliating pine and spruce forests and technogenic light forests under crowns compared to the background, which is due to the degradation of the state of the forest stand and crown defoliation. Table 5. Concentration of element compounds in the soil water from illuvial soil horizons (IS) in spruce and pine forests in 1995-2020, mg/L, ANC in mmol/L. The increase in water acidity and concentrations of the main pollutants in defoliating forests and pollution-induced sparse forests indicates a significant negative impact of industrial air pollution on forest biogeocenoses. The decrease in the ANC index in defoliating forests and pollution-induced sparse forests compared to the background areas is associated with a decrease in the concentrations of low molecular weight aliphatic acids and an increase in the concentrations of mineral acid anions.

Stages of
Removal of elements. In the background area in the Northern taiga forests, the removal of elements with lysimetric waters has high variability. Under the crowns, in the soil water from all soil horizons in spruce forests and from A0 and E + B horizons in pine forests, the removal of the majority of elements is significantly higher than between the crowns. This is due to the higher concentration of elements in the soil water under the crowns.
In defoliating spruce forests, the removal of elements with soil water is higher under the crowns than between the crowns. In defoliating pine forests in the soil water from the A0 horizon, the removal of Ca, Mg, C, Cu, and Ni increases in the water between the crowns. This can be explained by the large amount of soil water between the crowns compared to under the crowns. The removal of sulfates, copper, and nickel in the soil water in defoliating Northern taiga forests increases (p < 0.05) compared to the background area.
In pollution-induced sparse forests, the removal of elements with soil water under the crowns is higher than between the crowns. The removal of the main pollutants with water from soil horizons in sparse spruce and pine forests is significantly higher than in the background area. It should be noted that the amount of soil water from all soil horizons under the crowns in sparse forests is higher than in the background area, with the exception of soil water from the E + B horizon in pine forests. This can be explained both by an increase in the volume of precipitation due to thinning of the crowns and an increase in their transmission capacity when exposed to pollution and with a smaller volume of water absorption from the soil by the remaining living plants.
Inter-biogeocenotic variability. In the soil water from all soil horizons under the crowns of background spruce forests, the concentrations of Ca, Mg, SO 4 2− , and Ni in the soil water from the A0 horizon are significantly higher than in pine forests. The precipitation passing through the forest canopy becomes more concentrated and spruce undergoes a deeper transformation than pine, due to a longer and denser crown. Between the crowns, in the soil water from the A0 horizon, the concentrations of elements are either comparable or higher in pine forests. In the soil water from mineral horizons, generally, inter-biogeocenotic differences similar to the soil water from organogenic soil horizons were observed. The acidity of the soil water from all soil horizons under and between the crowns in pine forests is significantly higher than in spruce forests (with the exception of the soil water from the BC horizons between the crowns. This can be explained by a higher calcium content in the soils of spruce forests [1]. In defoliating spruce forests, in the soil water from all soil horizons under and between the crowns, the concentrations of most elements are significantly higher than in pine forests. In spruce forests, the ANC index in the soil water under the crowns is lower and, between the crowns, higher compared to pine forests. In the context of industrial air pollution, this can be explained by a higher flux of acid-forming substances under the crowns of spruce, known to intercept more acid-forming substances that the more open pine crowns. In pollution-induced sparce spruce forests, in the soil water from all soil horizons, the concentrations of sulfates, calcium, and nickel are higher than in pine forests. The water acidity is higher under the crowns of pine forests. The ANC index is lower in the soil water of spruce forests, both under and between the crowns, which is due to a more intense flux of acid-forming substances under the crowns of spruce forests.
Multi-year dynamics. In the Northern taiga forests, the long-term dynamics of element concentrations in soil waters are characterized by high variability.
In background spruce forests in the soil water from organogenic soil horizons, the long-term dynamics (in the period from 1999 to 2020) shows a trend of increasing Ni concentrations: under the crowns from 0.005 to 0.008 mg/L (R 2 = 0.47) and between the crowns from 0.002 to 0.009 mg/L (R 2 = 0.39). In the background pine forests, this trend was not found; on the contrary, there is a trend under the crowns toward decreasing Cu concentrations from 0.01 to 0.001 mg/L (R 2 = 0.46).
In defoliating forests in the soil water from A0 soil horizons, trends toward increasing concentrations of Ni, Cu, and SO 4 2− in recent years were found under the crowns (Figure 3). At this stage of degradation, a decrease in sulfate concentrations in the soil water from organic soil horizons between the spruce crowns from 10 to 5 mg/L (R 2 = 0.46) and from 10 to 2 mg/L (R 2 = 0.61) in pine was found. In spruce forests, there was a decrease in sulfate concentrations in soil water from 43 to 18 mg/L (R 2 = 0.33) (horizons E + B) and from 20 to 13 mg/L (R 2 = 0.61) (horizon BC).
In defoliating forests in the soil water from A0 soil horizons, trends toward increasing concentrations of Ni, Cu, and SO4 2− in recent years were found under the crowns ( Figure  3). At this stage of degradation, a decrease in sulfate concentrations in the soil water from organic soil horizons between the spruce crowns from 10 to 5 mg/L (R 2 = 0.46) and from 10 to 2 mg/L (R 2 = 0.61) in pine was found. In spruce forests, there was a decrease in sulfate concentrations in soil water from 43 to 18 mg/L (R 2 = 0.33) (horizons E + B) and from 20 to 13 mg/L (R 2 = 0.61) (horizon BC). In the soil water in pollution-induced sparse forests, a trend was found of decreasing concentrations of the main pollutants (Cu, Ni, SO 4 2− ). The most significant trends toward lower concentrations of copper, nickel, and sulfates were observed in the soil water from the organic soil horizon in sparse spruce forests, both under and between the crowns (Figure 4). The observed trends toward lower concentrations of the main pollutants at the pollution-induced sparse forest stage may indicate a decrease in the fallout of pollutants in the composition of large particles near the smelter.

Spruce Pine
Thus, despite the decrease in the level of pollutants in the emissions of the Severonickel Smelter, in the Northern taiga forests, the long-term dynamics of the concentrations of the main pollutants in atmospheric fallout and soil water in the background area and in defoliating forests demonstrate significant trends toward their increase, which is most pronounced for nickel concentrations. This may be due to an increase in nickel concentrations in aerosols propagating over considerable distances. However, it is worth noting the decrease in the concentrations of the main pollutants in atmospheric fallout and soil water in pollution-induced sparse spruce and pine forests, which may indicate a decrease in pollutant fallout as part of larger particles near the smelter.
In the soil water in pollution-induced sparse forests, a trend was found of decreasing concentrations of the main pollutants (Cu, Ni, SO4 2− ). The most significant trends toward lower concentrations of copper, nickel, and sulfates were observed in the soil water from the organic soil horizon in sparse spruce forests, both under and between the crowns (Figure 4). The observed trends toward lower concentrations of the main pollutants at the pollution-induced sparse forest stage may indicate a decrease in the fallout of pollutants in the composition of large particles near the smelter.
Under the crowns Between the crowns