Research Article |
Corresponding author: EI Kotova ( ecopp@yandex.ru ) Academic editor: Aleksandr I. Malov
© 2019 EI Kotova, IA Kuznetsova, VV Kriauciunas, SA Iglovsky, NS Larionov.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Kotova EI, Kuznetsova IA, Kriauciunas VV, Iglovsky SA, Larionov NS (2019) Contribution of the atmospheric channel to lead contamination of soils in the Arctic territories. Arctic Environmental Research 19(2): 56-64. https://doi.org/10.3897/issn2541-8416.2019.19.2.56
|
Although the Arctic territories have long been considered “the gold standard of purity”, global climate change and environmental pollution are having a significant impact on the state of Arctic ecosystems. In particular, industrial complexes combined with transboundary transport are having a negative impact. The aim of this work is to determine the contribution of atmospheric Pb to the contamination of soils of the Arctic territories using the trajectory statistics method which combines correlation and factor statistical processing approaches while taking the actual lead content in the upper soil horizon into account. In order to assess the transfer of pollutants, an analysis of multiannual air mass tracts and impurity transfer trajectories to certain points was carried out. The objects of study are the soil of the Khorey-Ver settlement, the basin of the Shapkin and Sula rivers, northeastern and southern Svalbard and the islands of the Southern Archipelago of Novaya Zemlya. It was established that the primary contribution to airborne pollution of the Nenets Autonomous District (NAO) during the summer period comes from sources located within a radius of about 400 km. Aerogenic influence of the mainland on the territory of Svalbard is practically absent, with the main atmospheric transport of substances coming from nearby marine areas. The transfer of impurities from the mainland to the Novaya Zemlya District is also insignificant. During the summer months, there is practically no airborne lead contamination of the island territories. Some influences on the territory of Novaya Zemlya and the NEO can be traced to sources on the Kola Peninsula, the coasts of the Barents and Kara seas, as well as long-distance transport from the southern regions. On the territory of the NAO, the number of heavy metal (HM) precipitations increases from west to east, possibly indicating a stronger impact of emissions from the Norilsk industrial hub than from the enterprises of the Kola Peninsula on this territory. In terms of seasonal dynamics, the island points are characterised by an increase in HM fluxes during the winter period; conversely, for the NAO territory, such an increase takes place during the summer. This discrepancy is primarily due to the changing nature of the atmospheric circulation.
lead, soil, Arctic Region, atmospheric transport
The impact of global climate change and environmental pollution on the state of the oceanic and marine ecosystems, as well as on land-based bodies of water, is especially significant in the Arctic territories. Until recently, the state of the Arctic territories was considered as a “gold standard of purity”. At the present time, the northern regions of Eurasia can be very conditionally considered as background in terms of anthropogenic impact on their environment. Both the long-range transport of anthropogenic pollution and the presence of powerful industrial complexes in the Arctic itself are having an impact. Despite the fairly large number of individual publications and reviews devoted to this issue (
The aim of this work is to determine the contribution of atmospheric Pb to the contamination of soils of the Arctic territories using the trajectory statistics method.
To date, there are a large number of methods for mathematical modelling of direct and inverse problems in mathematical physics, including problems of contaminant transfer.
A number of models have already been developed that describe processes of emission, distribution and neutralisation of polluting aerosol substances in atmospheric air. All of them, in one form or another, rely on the law of conservation of mass and use grid methods to solve the corresponding equation (
When using mathematical modelling methods for the formulation of a direct problem, the investigated physical phenomenon is described by partial differential equations with the corresponding initial and boundary conditions.
Three-dimensional hydrodynamic models are typically used for assessing the transfer of contaminants over a complex underlying surface (
The acronym LDPM (Lagrangian particle dispersion models) is used in the literature to describe the method of calculating turbulent diffusion on the basis of random trajectories of elements of atmospheric turbulence. The best-known models for calculating deterministic trajectories are implemented in the FLEXTRA (NILU) and HYSPLIT (NOAA) programs accessible via the Internet (
In order to calculate the return trajectories, the HYSPLIT4 model was used in combination with data from the reanalysis of NCEP / NCAR meteorological elements (
For six selected points (Fig.
Location of calculated points: T. 1 – Khorey-Ver, T. 2 – Shapkino, T. 3 – Sula, T. 4 – NW Svalbard, T. 5 – South Svalbard, T. 6 – Novaya Zemlya
The calculation was carried out for the winter and summer periods between 2010 and 2015.
Our calculations of long-range pollution transport in the Arctic are based on the statistical air mass trajectory method, which has been widely used in similar studies (
In recent years, the method of analysing the trajectories of air masses has been used not only to study the composition and genesis of specific air samples taken during experimental work, but also to study the average regularities in the processes of long-range air transport and its components in the atmosphere (
The basis of the approach used in this work is the analysis of multiannual arrays of air-mass and pollutant transfer trajectories to an arbitrary geographic location, the methodology of which was originally developed at the IFA RAS and improved over recent years (
The modelling of the distribution of contaminants during their long-range atmospheric transport is performed similarly to that detailed in the works (
On the grid (1° × 1°), the function Zij(X) was calculated for each cell with coordinates (ij) – a function of potential sources characterising all atmospheric processes determining the transfer of impurity X from a source m to point n: probability, duration, range, vertical deposition on the surface and horizontal spreading.
The value Zij is calculated according to (
where nij – the number of trajectory points in the cell; tij – the average time of air movement from a given cell to the point of observation; hij – the average height of the layer of mixing, calculated along the trajectories of air from the cell to the point of observation; Lij – the average length of these trajectories; N – the total number of points of daily trajectories for the billing period; b – the coefficient of expansion of the air flux in the horizontal direction (b ≈ 0.15); Кij – the average sedimentation velocity of the impurity on the surface as it moves from the cell to the observation point, which is chosen as the average of two values – for the cell (ij) and for the observation point cell.
It is assumed that the admixture is distributed uniformly in the vertical dimension within the surface mixing layer, whose height is calculated
simultaneously with the trajectories using the HYSPLIT4 model and then averaged over hij.
The parameters of the rate of removal of contaminants from the airstream are presented in (Vinogradova 2015). These correspond approximately to the classical climatic zoning of the territory.
As indicators of environmental pollution, Pb fluxes to the surface were calculated, consisting of chemically passive (along the transfer path) contaminants transported over long distances predominantly on aerosol particles of submicron size.
Data on emissions of lead into the atmosphere from Europe were selected from the EMEP database, which is accessible for interactive use via the Internet (
Within the framework of field studies, soil samples were selected from 75 sample areas in the territories shown in Fig.
The average lead content (in terms of mg Pb per kg of soil absolute dry weight.) In the upper soil layer was obtained from the results of X-ray fluorescence analysis of soil samples using a specialised technique (
Statistical data processing was carried out using software from StatSoft Inc. (2011). STATISTICA (data analysis software system), version 10. www.statsoft.com.
Statistical data processing consisted in establishing significant correlations between the considered parameters of the objects of study. Factor analysis using study.
The spatial distribution of the function Zij(X), reflecting the degree of influence of anthropogenic sources of pollution on the territories under consideration, is presented in Figures
The average heavy metal (HM) fluxes to the surface, which were calculated taking the emission of metal sources into account, formed the basis for calculating pair correlations, factor analysis and plotting the average lead content in the soil and fluxes given in Table
Territory | Pb fluxes, mkg·km-2·season-1 | Mean Pb concentration for surface soils, mg·kg-1 | |
---|---|---|---|
Winter | Summer | ||
Khorey-ver | 30.54 | 51.41 | 11.88 |
Shapkino | 32.82 | 31.20 | 9.57 |
Sula | 21.68 | 11.45 | 5.41 |
Northeast Svalbard | 18.78 | 5.55·10-5 | 2.47 |
South Svalbard | 13.59 | 2.03·10-6 | 2.47 |
Novaya Zemlya | 3.96 | 95.19·10-3 | 6.15 |
Paired correlations between the average lead content in soils and fluxes are shown in Table
Pb fluxes | ||
---|---|---|
Winter | Summer Pb | |
Average Pb content in soils | 0.6453, p = 0.166 | 0.9328, p = 0.007 |
Factor 1 | Factor 2 | Factor 3 | |
---|---|---|---|
Pb flux winter | 0.75* | 0.00 | 0.66 |
Pb flux summer | 0.93 | 0.32 | 0.19 |
Average Pb content in soils | 0.78 | 0.60 | 0.11 |
Explained Variance | 6.88 | 4.87 | 2.19 |
Proportion of Total Variance | 0.49 | 0.35 | 0.16 |
Significant correlations of lead content in soils with p <0.01 with fluxes in the summer period should be noted.
The graph of average lead content in soils versus flux during the summer period is shown in Figure
The spatial distribution of the function Zij(X) (Figures
As seen in Figures
The obtained calculations are consistent with the data of Roshydromet on the burden of atmospheric precipitation of lead in the area (
According to the calculations, the minimum delivery of lead pollutants, both in the summer and in the winter periods, is expected to enter the territory of the Novaya Zemlya archipelago. During the summer months, there is practically no airborne lead contamination of the island territories.
On the territory of the NAO, the quantity of lead precipitation increases from west to east, possibly indicating a stronger impact of emissions from the Norilsk industrial hub than from the Kola Peninsula on this territory.
In terms of seasonal dynamics, the island points are characterised by an increase in HM fluxes during the winter period; conversely, for the NAO territory, such an increase takes place during the summer. This may be due to the summer contribution of the wind uplift of lead in the composition of fine fractions of eroded soils (terrigenous component), as well as a change in the character of the atmospheric circulation.
Correlation analysis of data on lead fluxes and data on its average content in the soil of the studied territories showed that the summer flux makes a significant contribution (Table
Factor analysis using the varimax method allowed three factors to be identified, of which only the first (49%) is closely related to lead fluxes in the summer and winter months, playing a significant role in the formation of the soil composition of the studied territories (Table
It is established that the primary contribution to airborne pollution of the NAO during the summer period comes from sources located within a radius of about 400 km. During the summer, aerogenic influence from the mainland onto the territory of Svalbard is practically absent, with the main atmospheric transport of substances coming from nearby marine areas. The transfer of impurities from the mainland to the Novaya Zemlya District is also insignificant. Therefore, there is practically no airborne lead contamination of the island territories during the summer months. At the same time, some influence of sources from the Kola Peninsula, the coast of the Barents and Kara Seas on the territory of Novaya Zemlya and NAO, can also be traced. A contribution of long-range transport from the southern regions is also made to the airborne pollution of the NAO and Novaya Zemlya territories. On the territory of the NAO, the quantity of lead precipitation increases from west to east, possibly indicating a stronger impact of emissions from the Norilsk industrial hub than from the Kola Peninsula on this territory. In terms of seasonal dynamics, island points are characterised by an increase in lead fluxes in winter due to the addition of remote sources. In the case of the NAO, this factor is present during the summer season due to dust and soil contributions.
The present work was carried out with the financial support of a grant for the implementation of the state task “Complex isotope-geochemical studies of the quality of the environment and identification of its transformation processes in the coastal territories of the European North both present and past” (No. AAAA-A19-119011890018-3).