Next Article in Journal
Effect of Corporate Social Responsibility on Community Resilience: Empirical Evidence in the Nickel Mining Industry in Southeast Sulawesi, Indonesia
Next Article in Special Issue
Exploring the Optimal Cost-Benefit Solution for a Low Impact Development Layout by Zoning, as Well as Considering the Inundation Duration and Inundation Depth
Previous Article in Journal
Effects of SOP Lagoon Additive on Gaseous Emissions from Stored Liquid Dairy Manure
Previous Article in Special Issue
Watershed-Based Evaluation of Automatic Sensor Data: Water Quality and Hydroclimatic Relationships
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Disproportionate Water Quality Impacts from the Century-Old Nautanen Copper Mines, Northern Sweden

1
Department of Physical Geography, Stockholm University, 10691 Stockholm, Sweden
2
Faculty of Geography, M.V. Lomonosov Moscow State University Leninskie Gory 1, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(4), 1394; https://doi.org/10.3390/su12041394
Submission received: 23 January 2020 / Revised: 7 February 2020 / Accepted: 11 February 2020 / Published: 13 February 2020
(This article belongs to the Special Issue Watershed Modelling and Management for Sustainability)

Abstract

:
Pollution from small historical mining sites is usually overlooked, in contrast to larger ones. Especially in the Arctic, knowledge gaps remain regarding the long-term mine waste impacts, such as metal leakage, on water quality. We study the small copper (Cu) mines of Nautanen, northern Sweden, which had been in operation for only six years when abandoned approximately 110 years ago in 1908. Measurements from field campaigns in 2017 are compared to synthesized historical measurement data from 1993 to 2014, and our results show that concentrations of Cu, Zn, and Cd on-site as well as downstream from the mining site are order(s) of magnitude higher than the local background values. This is despite the small scale of the Nautanen mining site, the short duration of operation, and the long time since closure. Considering the small amount of waste produced at Nautanen, the metal loads from Nautanen are still surprisingly high compared to the metal loads from larger mines. We argue that disproportionately large amounts of metals may be added to surface water systems from the numerous small abandoned mining sites. Such pollution loads need to be accounted for in sustainable assessments of total pollutant pressures in the relatively vulnerable Arctic environment.

1. Introduction

Historical mining operations commonly deposited waste rocks, tailings, and metallurgical slags on pristine ground or in nearby lakes or streams [1]. As a consequence, historical mine waste has direct contact with terrestrial and aquatic environments and ecosystems, which implies a high risk for the development of acid mine drainage if the waste contains sulfidic material. Compared to present mining operations in most countries, where mine effluents are becoming increasingly controlled, abandoned historical mine waste represents an uncontrolled source of metal pollution at the landscape level [2], which hampers long-term sustainable resource management.
Long-term contamination from historical mining is usually reconstructed from bed sediments sampled from downstream surface water systems. For instance, sediments from the Big River, central USA, have elevated metal concentrations along a 170-km stretch as a result of 120 years of mining [3]. Mass flows of metals through surface water, on the other hand, directly reflect continued loads of mobile metal phases into freshwater ecosystems. Considering riverine mass flows of copper (Cu), the historical mining districts in the Odiel drainage basin in Spain have, for example, been estimated to contribute as much as 1300 tons of Cu per year to the surface water system [4]. Within the drainage basin, which covers 2300 km2, mining operations began in prehistoric times [5], which could explain the magnitude of the Cu mass flows. A release of up to 160 tons of Cu per year was estimated for the Britannia mine, in south-western Canada, to coastal areas along an 80 km underground tunnel system [6]. The more than 1000-year-old Falun mine, in central Sweden, released (pre-remediation) about 30 tons Cu per year to downstream surface waters [7]. Analogously, much of the existing literature on historical mining regards large scale and relatively well-defined “worst-case” pollution sites where a massive metal pressure on aquatic and terrestrial ecosystems is obvious.
Research on long-term (> 100 years) mine waste degradation and spreading processes are essentially lacking from the Arctic, in contrast to the relatively well documented impacts of historic mines, e.g., in the Mediterranean or from the temperate climate zones of Europe or North America. Follow-up long-term monitoring and groundwater data are most often not available from northern latitudes [8]. This knowledge gap should be of increasing concern, as the Arctic is likely to face high pressures from mineral extraction in the near future [9]. In addition, the Arctic is already facing multiple pressures on the biotic and abiotic environment due to climate change and other societal activities [10].
Current knowledge of historical mine waste impacts in the Arctic has mostly been gained from sites that were abandoned during the last few decades. For example, the Black Angel mine in western Greenland closed in 1990 and its subsea tailing deposit was found to contribute to metal accumulation in local fish tissue [11]. Studies on older mine waste are rare, which hampers the understanding of long-term impacts of Arctic mining operations. Existing (70 years old) examples include the Sherritt-Gordon mine in central Canada and the Laver mine in northern Sweden, which, currently, continue to release metals to their respective local environment [12,13]. Furthermore, the impact of the numerous, small mining sites (e.g., < 1 million tons of produced tailings) are poorly investigated, which inhibits, e.g., the upscaling of total historical loads to regional scales.
The Nautanen Cu mining site in northern Sweden is a clear example of a small historical site comprising mine waste from only six years of active production more than 100 years ago (1902–1908). Existing monitoring data from the Nautanen site provide us with a rare opportunity to study the long-term environmental impacts from a well-defined abandoned Arctic mining site. This may increase the understanding of historical mining impacts in other regions of the Arctic, including Russia [14,15] and Finland [16]. For example, soil Cu contamination has previously been found impact biota in the vicinity of the mine tailings, e.g., lowering the cold tolerance of the oligochaete worm [17]. The water-borne spreading of Cu and its potential consequences for ecosystem functions needs to be further studied, especially in colder climates. The hereto unexploited Cu ore at Nautanen is currently under consideration for extraction, which highlights the need for a thorough characterization of the site’s current metal contamination. The main objectives of this study were therefore to examine (i) the present situation and potential temporal trends in how the Nautanen mines impact surface water metal concentrations (relative to regional backgrounds), (ii) present and cumulative metal riverine mass flows since the mine closed, and (iii) the downstream mass flows in relation to the amount of deposited mining waste within the catchment, facilitating a comparison of impacts with other Arctic mines (of various sizes).

2. Materials and Methods

2.1. Site Description

2.1.1. Historical Background

The Nautanen copper mines were owned by the company Nautanen Kopparfält AB. Mining production begun in 1902 and was in operation only for six years before shutting down in 1908 due to unexpectedly limited ore reserves [18]. About 72,000 tons of unprocessed ore were mined during this short period, which resulted in 5700 tons of Cu concentrate and 4600 tons of iron (Fe) concentrate [19]. Assuming that all the mined Cu concentrate was chalcopyrite (CuFeS2 with 35% Cu [20]), which is the most common Cu mineral at Nautanen, the total Cu production was approximately 2000 tons. The mining site included five larger mines named: Dagny, Hoppet, Maria, Max, and Fredrik (Figure 1), a small residential area and a cable railway system. The crusher, smelter, and power station were located in the central parts of the mining site, hereafter referred to as the industrial area (Figure 1) [21]. As of 2017, there was about 10,000 tons (5500 m3) of unconfined tailings and up to 10,000 tons (5000 m3) of metallurgical slag products in the industrial area [22]. After the mining site closed, the buildings and infrastructure components were removed and Nautanen was abandoned a few years later. While in operation, the mines at Nautanen were estimated to have produced a total of 68,000 tons of waste rock, which mostly were dumped outside the Max and Maria mines [18]. As the waste rock was found to contain a relatively high Cu grade, most of the dry-lying heaps were removed in 2005–2008 for re-processing at the nearby, currently active Cu mine, Aitik, which is owned by Boliden AB [23]. Waste rocks were still present in 2017 in the Maria Lake (Figure 1) and in smaller heaps spread across the mining site. So far, no thorough remediation has been implemented at Nautanen. The area is primarily used for recreational purposes and winter-spring reindeer herding by the Baste čearru Sami community [24]. An exploration permit has allowed Boliden AB to investigate the Nautanen area for possible new mining since 2009 [25]. Some relatively recent publicly available documents therefore describe the ore deposits and potential mining operations [26].

2.1.2. Climate and Environmental Setting

The annual average temperature is −1.6°C with summer (June–August) and winter (December–February) temperatures averaging around 11 °C and −13 °C respectively (based on data of the 25-year period 1993–2017 from the Climatic Research Unit (CRU) at the University of East Anglia described by Harris et al. [27] (see Supplementary Materials S1 for details). The same dataset showed that annual average temperatures have increased by 1.0 °C during the last 110 years (1908–2017). The annual average precipitation is 560 mm/yr (for 1993–2017) and compared to the first quarter of the century (1901–1925), total annual precipitation has increased by 23% (based on CRU data). Nautanen is located within a larger active mining district where the Malmberget mine (Fe) and Aitik mine (Cu) are the two largest mines in operation within a 10–15 km distance (Figure 1). Nautanen and Aitik are connected through the Nautanen Deformation Zone (north-north-west direction), which consists of metamorphosed volcanoclastic sedimentary rocks with Cu-Au ± Fe mineralization [28]. The formation hosts mainly sulfide deposits of chalcopyrite (CuFeS2), pyrite (FeS2), and magnetite (Fe3O4) [20]. The landscape is characterized by Veiki moraine (ice walled lake plains) where the geomorphology of the Nautanen area is characterized by exposed bedrock or till [29]. The headwaters of the Imetjoki stream system originate from two smaller lakes situated about 1.5 km north of the Nautanen mining site (Figure 1). The Imetjoki stream flows along the northern side of the mining site partly passing through wetlands. The stream drains the Northern Lake and Imetjärvi, the latter being the outlet for the Maria Lake. In the industrial area, Imetjoki flows across the tailings where the groundwater level has been measured to be located at only 0.5–1.5 m depth [18]. Downstream of the Max mine the groundwater level was estimated to be located below 3 m depth [22]. Imetjoki joins the Nietsajoki stream at 4 km downstream of the mining site, which gives the Imetjoki a drainage area of 6.6 km2. Nietsajoki joins the Lina River another 8 km downstream and then eventually drains into the Kalix River and the Baltic Sea.

2.2. Methods

2.2.1. Field Measurements and Data Synthesis

Two field campaigns were conducted at Nautanen in 2017; one during May (snowmelt conditions) and one in August (after ground thawing; Figure 1). Water quality sampling and discharge measurements were performed simultaneously at nine sampling points, complemented with seven additional sampling points for water quality only. The sampling points were chosen to represent four zones: (i) local background conditions at locations judged to be unaffected by the past mining activities (i.e., in the Imetjoki stream, upstream of the main mining zone, and in two reference rivers in the close vicinity; Figure 2a), (ii) at the main mining zone (Figure 2b), (iii) directly downstream of the main mining zone (Figure 2c; D1 in Figure 1), and (iv) 4 km downstream of the main mining zone (Figure 2d; D2 in Figure 1).
The two reference rivers were the Nietsajoki stream and the Lina River (Figure 1). Water samples for analyses of metal concentrations and base chemistry were taken in the streams in triplicates for total concentration and dissolved concentration (having passed a 0.22 µm Filtropur filter), and acidified according to the sampling protocol in the Supplementary Materials S2. pH, electric conductivity, alkalinity, water temperature, and turbidity were also measured in situ (see Section S2 in the Supplementary Materials for instruments and sampling details). The stream discharge at each point was estimated from measuring water velocity with an A. Ott Kempten propeller current meter (at 60% of the total depth from water surface at each width increment) by wading along a cross section of the stream.
In addition to the conducted local Nautanen measurements, we considered regional background concentrations reported by Herbert et al. [30], reflecting element concentration ranges observed in northern Sweden between 1996 and 2007 in streams judged to be unaffected by point pollutant sources. Existing historical field measurement data on water quality in the Nautanen region were also synthesized from publicly available documents both online and in state archives, mostly in the form of technical consultancy reports ordered by governmental agencies. More specifically, water quality measurements have been conducted at Nautanen on several occasions and open data were available from 1993 to 1994 [19,31], 2001 [18], 2005 to 2009, and 2014 [23]. Laboratory analyses and details of all measured and historical samples are described in Section S2 in the Supplementary Materials. Specifically, the elements copper (Cu), zinc (Zn), cobalt (Co), cadmium (Cd), nickel (Ni), molybdenum (Mo), barium (Ba), manganese (Mn), aluminum (Al), and iron (Fe), together with dissolved organic carbon (DOC) and sulfate (SO4), were analyzed in detail due to their common association with mining activities and related geochemical processes [32,33,34].

2.2.2. Mass Flow Analysis

Mass flows of elements at Nautanen were calculated in two ways considering the Imetjoki stream and its tributaries; as snapshot mass flows, and as long-term average mass flows. In addition, the cumulative element mass carried by Imetjoki was estimated over the approximately 110 years that has passed since the mine closed in 1908. The quantifications considered conditions in four zones (or locations): (i) upstream of the main mining zone, (ii) at the main mining zone, and downstream of the main mining zone at locations (iii) D1 and (iv) D2 (Figure 1 and Figure 2).
The snapshot element concentration within each sampling zone was calculated by averaging the inflowing tributary water with the main Imetjoki stream water. Specifically, for sampling points within a considered zone j, with measured discharge q1 and element concentration c1, where Imetjoki was joined by a tributary with measured discharge q2 and measured total element concentration c2, the snapshot element concentration in Imetjoki immediately after the confluence point was calculated as a flow weighted average concentration, c ¯ i . This was done for each of the four sampling zones (j; upstream, main mining zone, D1 and D2) along the main Imetjoki stream according to:
c ¯ j = c 1 · q 1 + c 2 · q 2 q 1 + q 2 .
The corresponding snapshot element mass flow for sampling zone j, mfj, was calculated as:
m f j = c ¯ j · ( q 1 + q 2 ) .
Concentration enrichment factors and relative mass flow increases were calculated by normalizing the respective zonal concentration and mass flow values by the corresponding values for the upstream zone. Since upstream sampling could not be performed in the May field campaign (due to snow conditions), we only used results from the August field campaign (where sampling in all zones were performed) in the mass flow analysis.
The long-term average mass flows ( M F ¯ L , j ) for a 25-year period (1993–2017) in each of the sampling zones j were calculated by:
M F ¯ L , j = C ¯ L , j · Q ¯ L , j ,
where C ¯ L , j is the long-term average concentration for each of the sampling zones, and Q ¯ L , j is a corresponding long-term average discharge. Regarding C ¯ L , j , we used the assumption of (approximate) geochemical stationarity, i.e., that the inter-annual metal concentration variation is typically low compared to the much higher variability in discharge [35]. This allowed us to estimate C ¯ L , j   as   approximately   equal   to   C ¯ j , i.e., that long-term average concentration would be equivalent to our average zonal snapshot concentration ( C ¯ j ). Concerning the long-term average discharge, Q ¯ L , j , we used a water balance approach based on monthly precipitation and temperature measurements provided by the Climatic Research Unit (CRU) at the University of East Anglia [27] to estimate the average discharge at the downstream location D2 ( Q ¯ L , D 2 ; see Section S3 in the Supplementary Materials for water balance calculations), for the period 1993–2017. Our resulting evapotranspiration (estimated from equations from Turc [36] and Langbein [37]) was furthermore compared to independent evapotranspiration quantifications based on available, region-specific precipitation and stream flow data from the Swedish Meteorological and Hydrological Institute (SMHI). Q ¯ L , D 2 was further scaled to quantify the discharge in the upstream sampling zones, based on the relative catchment area of each sampling zone (Aj), yielding the following equation for calculating long-term average discharge of sampling zone j:
Q ¯ L , j =   Q ¯ L , D 2 · A j A D 2 .
The average mass flows at the downstream location D2 ( M F ¯ L , D 2 ) were furthermore compared to the total mass of deposited tailings and slag products at the industrial area to yield a so-called “stream load-to-tailings” ratio. The industrial area was considered as the main pollutant source within the catchment [22] and the ratio was, therefore, calculated by dividing the amount (mass) of tailings and slag products with M F ¯ L , D 2 . This stream load-to-tailings ratio can be seen as a first-order indicator, facilitating order-of-magnitude comparisons of environmental impacts between mining sites of different sizes. The ratio is mainly relevant when particulate flows are negligible (as in Nautanen), in which case, the geochemical impacts from land disturbance can be approximated by the mass of deposited wastes. For comparison, under conditions where particulate flows are dominant, one can alternatively consider surface erosion rates per unit of disturbed lands (e.g., Walling and Fang [38]).
Last, the 110-year (1908–2017) cumulative mass for each sampling zone ( M 110 , j ) was estimated by:
M 110 , j = C ¯ 100 , j · Q ¯ 110 , j · 110 .
Similar to Equation (3), C ¯ 100 , j is here used as equivalent to the snapshot C ¯ j . Q ¯ 110 , j is the average discharge for the 110-year period (1908–2017; see Section S3 in the Supplementary Materials).

3. Results

3.1. Water Quality

The stream discharge measured at Nautanen in May 2017 was relatively low for the season. Much of the winter snow still remained on the ground, without contributing to melt water-related flow peaks that frequently occur in May. The August 2017 discharge levels were consistent with previously recorded measurements (August–September in 2001–2002 [18]). Both in May and August, most elements existed predominantly in the dissolved phase (a t-test showed no significant difference at p = 0.05 between dissolved and total concentrations). The pH ranged between 3.1 and 6.4 at the industrial area (i.e., where the tailings are located within the main mining zone). The same pH-range (3.8–6.3) was measured at the lake outlets in the upstream zone. Other sampling points within the main mining zone had a pH around 6.1, while the downstream locations D1 and D2, as well as the reference Nietsajoki and Lina rivers had pH closer to 6.8 (see Tables S2 and S3 in the Supplementary Materials for a complete list of all water quality parameters and discharge).
Both of the field campaigns conducted in 2017 (May and August) are represented in Figure 3, where the total (unfiltered) concentrations for selected elements are grouped by sampling zone. Each element is shown together with its corresponding regional background concentration range (shown in grey; no regional values were available for Mo, Ba, and SO4). Our results show that, except for Fe, all metal concentrations in the water samples collected from the main mining zone (red boxes; Figure 3) were significantly higher compared to local background concentrations (green boxes) and concentrations at the downstream locations D1 and D2 (blue boxes). For instance, the median Cu, Zn, Co, and Ni concentrations at the main mining zone exceeded median local background concentrations by as much as one to three orders of magnitude. Many metals (Cu, Zn, Co, Cd, and Mn) additionally showed considerably higher mining zone concentrations than the regional background concentrations. Even at the downstream locations D1 and D2, the metal concentrations were in most cases elevated compared to background values. Notably, the trend of Fe differed from other metals, following instead the same trend as DOC in that the downstream concentrations in general exceed those of the mining zone.

3.2. Concentration Enrichment and Element Mass Flow Increases

Concentration enrichment factors and snapshot mass flow increases are presented in Figure 4 for the selected elements considering the main mining zone (red bars) and downstream locations (blue bars), relative to upstream conditions (green area), where the enrichment factors and relative increases of mass flows per definition are of equal unity. The highest concentration enrichment factor of all investigated elements was found for Cu at the main mining zone, with a value of about 140 (Figure 4a). The mass flow of Cu increased even more than that, by 290 times over the main mining zone (Figure 4b), which is due to a higher stream water discharge at the main mining zone than upstream of it. Many elements (Cu, Zn, Co, Cd, Ni, Mn, and SO4) showed considerable concentration enrichment (factors >20; average 65) as well as relative mass flow increased (>50; average 140) over the main mining zone. At the downstream location D1, these elements displayed distinctly lower relative mass flow increases, although the increases were still greater than or equal to 8 (average 26). The downstream location D2 displayed similar mass flow increased relative to the upstream conditions at the D1 location. In contrast, due to dilution effects, the concentration enrichments were much less pronounced at the D2 location than at D1. Other elements, like Fe and DOC displayed no or low concentration enrichments (factors near unity) over the mining site, whereas mass flows still increased all the way along the main stretch of the stream, due to the element supply with inflowing water from tributary streams.
Considering the surveyed 25-year period (1993–2017), the average mass flow of Cu leaving the upstream zone was estimated at 1.2 kg/yr, whereas the mass flow of Cu leaving the main mining zone was estimated at 450 kg/yr. At both the downstream locations D1 and D2, the mass flow of Cu was lower, although still considerable at about 70 kg/yr. The identical average mass flow values of Cu at these two locations (D1 and D2) implies that no major net addition or net removal of Cu occur along the 4 km pathway that separates the points. These mass flows were based on long-term average discharges of 450,000 m3/yr at the D1 location and of 2,000,000 m3/yr at the D2 location for the period 1993–2017 (i.e., Q ¯ L , D 1 and Q ¯ L , D 2 , respectively). During this 25-year period, the annual discharges showed a coefficient of variation of 0.23. Q ¯ L , D 2   was in turn derived from a calculated actual evapotranspiration (ETa) of 260 mm/yr, which is within the reported range of 260–310 mm/yr for 1981–2010 for watersheds close to Gällivare City, as calculated by the Swedish Meteorological and Hydrological Institute [39]. Previous estimates of annual mass flows of Cu ranged between 78 and 240 kg/yr for the main mining zone. The total Cu mass that has been transported through each sampling zone during the approximately 110 years that has passed since the mine closed (1908–2017) was calculated to 120 kg from upstream of the main mining zone, 43 tons for the main mining zone, and 6.7 tons for the downstream locations D1 and D2. During these 110 years the discharge was estimated to increase with approximately 35% from the first quarter of the century (1901–1925) to the last (1993–2017), in response to changing hydroclimatic conditions. Finally, the stream load-to-tailings ratio for Nautanen was estimated at 1:105 per year based on the average Cu load in surface waters 4 km downstream of the main mining zone (70 kg/yr) and the amount of deposited tailings and slag products in the industrial area (20,000 tons).

3.3. Historical Cu Measurements

The historical measurements of total Cu concentrations in Nautanen are synthesized in Figure 5, which shows the results from 1993 to 2014, together with the present field campaign results from 2017. The scattering in local background levels reflects different sampling locations between the various historical measurement campaigns. This is also the case for the wide spread Cu concentrations at the mining zone. The highest on-site Cu value was measured in 2017, pertaining to a tributary stream flowing directly over the tailings that had not been monitored before. A general spread in concentrations across the main mining zone has been seen throughout the different measurement campaigns conducted between 1993 and 2017. However, the Cu concentration values at the downstream locations D1 and D2 has been similar over this time period. The lowest D1 concentration in Figure 5 (at 20 µg/L) pertains to a small stream downstream of the Fredrik mine (Figure 1), and unlike the other D1 measurement locations, is therefore not representative of the main Imetjoki channel. The Cu concentrations after Imetjoki and Nietsajoki joins (smaller dark blue points in Figure 5) dropped by one order of magnitude, due to dilution from the much higher discharge of Nietsajoki. These concentrations (3.9–7.9 µg/L) are still higher than most of the local background values (see Table S4 in the Supplementary Materials for a full list of synthesized historical data). For other elements than Cu the historically available records show little information. Zn and Cd were measured at a few locations in 1993–1994, 2001, and 2014, and the results show similar concentration ranges in their corresponding sampling zone as were found in 2017 (Figure 3).

4. Discussion

The historical activities at the Nautanen mining site cover a short time period (1902–1908) and are well documented. From the results of this study, we inferred that Cu concentrations were relatively constant during the surveyed 25-year period (1993–2017). For example, average Cu concentrations just downstream of the Nautanen mining site, and 4 km downstream of it, were 140 and 43 µg/L, respectively, according to measurements during the period of 1993–2014. These concentrations were similar to the average values of the conducted 2017-campaigns (130 and 38 µg/L respectively). Furthermore, our results showed that the average Cu loads in surface waters 4 km downstream of the main mining zone (70 kg/yr during 1993–2017) relative to the amount of tailings and slag products produced (20,000 tons), i.e., the stream load-to-tailings ratio, was to 1:105 per year.
The much larger abandoned Laver Cu mine in Sweden (1.2 million tons tailings) and the Sherritt-Gordon Cu mine in Canada (7.4 million tons tailings) both have stream load-to-tailings ratios of 1:107 per year, based on the stream mass flows of Cu at similar downstream distances as for Nautanen [12,13]. This means that these mining sites, compared to Nautanen, release significantly less Cu to surface waters relative to their total waste production. The tailings of the Sherritt-Gordin mine are, for instance, roughly 370 times larger than at Nautanen, whereas the Cu release rate is only about six times larger. Both the Laver and Sherritt-Gordon sites share similar characteristics as Nautanen; e.g., climatic conditions, a long time (70 years) since abandonment, and (mostly) un-remediated waste.
Our results therefore demonstrate that the combined impact of the remaining factors is large, including, e.g., impacts of varying abundance and availability of Cu in the host rock, weathering potential of the waste, and hydrological flow paths, to mention a few [32]. Notably, the existence of systematic differences in stream load-to-tailings ratios between small and large mining sites cannot be excluded. Some governing factors may, for example, be scale dependent, such as the surface area to volume ratio of the waste heaps, which, in turn, can impact weathering rates. The fact that the small Nautanen mining site had disproportionately high environmental impact, compared to the few existing quantifications of larger Arctic mining sites, illustrates a need to increase the general understanding of complex hydrogeochemical processes across different Arctic mining sites.
In Sweden, closed and abandoned sulfide mines have been estimated to load surface waters with a total of about 12 tons/yr of Cu [40]. This estimate is, however, highly uncertain and is based on a survey where the relevant County Administrative Boards provided data for monitored sulfide mines, which only represent a small fraction of the roughly 1000 known abandoned sulfide mines in Sweden [41]. Most of the smaller mining sites are unmonitored and could, therefore, not be included in those data sets. A similar regional-scale synthesis of historical mining pollution in England and Wales estimated the Cu loading of surface waters to be 19 tons/yr [2]. Thus, pollution from abandoned mines is associated with a substantial metal mass flux at the landscape-level and unaccounted-for smaller mines (similar to Nautanen) may increase this estimate even more.
Whereas the above-discussed metal concentrations and loads regard conditions in the surface water systems downstream of the main mining zone of Nautanen, our results additionally showed that the on-site stream water, at the middle of the mining site, has much higher aqueous concentrations and hence carries considerably higher metal loads (with, for example, a flow weighted average Cu concentration of 990 µg/L and a load of 450 kg/yr). Considering the observed absence of temporal trends during the last quarter of a century, together with the assumption that aqueous concentrations and metal loads should not have been lower than that in the beginning of the 20th century, when the mine was new, we estimate the total, on-site Cu loading of the stream corresponds to at least 43 tons since the closure of the mine (1908–2017).
For comparison, the total amount of produced Cu until 1908 was 2000 tons. Thus, for every ton of produced Cu, about 20 kg has been dissolved into the on-site stream water. Taken together, apart from observed spreading of metals to downstream waters, the above results reflect an on-going and considerable spatial re-distribution of metals within the local mining site. Specifically, (I) a net dissolution of metals in the upper parts of the main mining zone, leading to the observed order-of-magnitude in-stream mass flow increases (the difference between the green and red bars in Figure 4b), and (II) a net retention of metals in the lower parts leading to order-of-magnitude in-stream mass flow decreases (the difference between the red and blue bars in Figure 4b). In particular, there was a small stream that carried very high mass flows of Cu, Zn, Co, Cd, and Ni in the middle of the Nautanen mining site, just upstream of the junction with the main stream Imetjoki.
The observed decrease in Cu mass flows of about 380 kg/yr (450 − 70 = 380) as the small stream joins Imetjoki could potentially be explained by precipitation from the aqueous phase, e.g., as a result of diffuse inflows of groundwater with a different chemical composition (as previously seen in, e.g., Kimball et al. [42] and Palumbo-Roe and Dearden [43]). Even if the precipitated metals (Cu) would have accumulated primarily in the stream sediments, they are at risk of re-mobilizing in the future due to changing biogeochemical conditions or due to extreme weather events. In any case, the vegetation at and near the main mining zone of Nautanen has already been shown to be negatively impacted by the mining activities [23].
Projected future warming (+4 and +5.6 °C for RCP 4.5 and RCP 8.5, respectively) by the year 2100 for Norrbotten County is expected to shorten the period during which streams and lakes are frozen [44]. Specifically, winter (December-February) temperatures are predicted to increase the most, which will increase the likelihood of precipitation occurring as rain instead of snow, thus enabling an increased movement of water and metals through the water system over the year. This also implies that the water equivalent of the accumulated winter snow pack may decrease, which, in turn, would reduce the amplitude of the spring snowmelt runoff [10]. Even if high meltwater flushes were becoming less prevalent, peak precipitation events have the potential to increase erosion and the transport of material from the freely exposed tailings at Nautanen. Furthermore, a potential effect from a projected increased and more variable precipitation (by +20% and +30% for RCP 4.5 and RCP 8.5, respectively) is a higher mean groundwater table and larger level fluctuations around the mean value. In such conditions, the washing-out of metals from tailings to the mobile aqueous phase at Nautanen would increase [12].

5. Conclusions

From the results of the field measurements at the Nautanen mining site (1993–2017), we conclude that:
  • Despite the small spatial scale of the Nautanen mining site, the short duration of operation, and the long time (approximately 110 years) since closure, the average concentrations of Cu, Zn, and Cd on-site (990, 280, and 1.0 µg/L, respectively) and downstream of the mining site (150, 50, and 0.2 µg/L, respectively) were generally considerably above local (7.2, 3.0, and 0.010 µg/L, respectively) and regional (1.6, 11, and 0.10 µg/L, respectively) background values. In particular, downstream Cu concentrations were consistently high throughout the surveyed 25-year period (1993–2017).
  • The on-site mass flows of Cu, Zn, Co, and Cd were estimated to be between 100 and 300 times higher than upstream of the mining site. For Cu, the average mass flows were 1.2 kg/yr upstream of the main mining zone, 450 kg/yr on the main mining zone, and about 70 kg/yr 100 m to 4 km downstream of the site. Many metals exhibited a similar spatial pattern, indicating an on-going considerable retention of metals at or near the main mining zone, corresponding to 380 kg/yr for Cu.
  • Compared to other major abandoned mines in the Arctic, the metal loads from the Nautanen mining site were found to be unexpectedly high relative to the (small) amount of tailings and slag products produced. More generally, our approach based on stream load-to-tailing ratios show that small abandoned mining sites, which are numerous, could add disproportionately large amounts of metals to the surface water systems. Such effects need to be accounted for in assessments of total pollutant pressures in the relatively sensitive Arctic environment.

Supplementary Materials

The following are available online at https://www.mdpi.com/2071-1050/12/4/1394/s1. This supplementary material provides a description of the field measurement sampling protocols and water balance calculations. It also tabulates water quality parameters from the field measurements in 2017, and synthesizes historical data (period: 1993–2014). S1: Temperature and precipitation data, S2: Water sampling procedure and laboratory analyses, Table S1: Water sampling protocol for the Nautanen field campaign 2017, S3: Water balance calculations, Table S2: Results in base chemistry of surface water samples from the 2017 measurement campaign at Nautanen, Table S3: Results of total element concentrations in surface water samples from the 2017 measurement campaign at Nautanen, Table S4: Synthesis of historical field measurement campaigns at Nautanen.

Author Contributions

All authors participated in the field work; formal analysis: S.F.; data curation: S.F.; Supervision: J.J. and G.R.; writing-original draft: S.F.; writing-review and editing: J.J., G.R. and S.R.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Nordforsk Centre of Excellence for Resource Extraction and Sustainable Arctic Communities (REXSAC, project no. 76938) with additional funding from Göran Gustafsson Foundation for field work. S.C. was supported by Russian Science Foundation project no. 18-17-00086.

Acknowledgments

We acknowledge Carl-Magnus Mörth (Dept. of Geological Sciences), Marcus Sundbom, and Pär Hjelmquist (Dept. of Environmental Science) at the laboratories at Stockholm University for analyzing water quality data. We also thank Camilla Winqvist for providing historical archive data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lottermoser, B. Mine Wastes: Characterization, Treatment and Environmental Impacts, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2010; ISBN 978-3-642-12418-1. [Google Scholar]
  2. Mayes, W.M.; Potter, H.A.B.; Jarvis, A.P. Inventory of aquatic contaminant flux arising from historical metal mining in England and Wales. Sci. Total Environ. 2010, 408, 3576–3583. [Google Scholar] [CrossRef] [PubMed]
  3. Pavlowsky, R.T.; Lecce, S.A.; Owen, M.R.; Martin, D.J. Legacy sediment, lead, and zinc storage in channel and floodplain deposits of the Big River, Old Lead Belt Mining District, Missouri, USA. Geomorphology 2017, 299, 54–75. [Google Scholar] [CrossRef]
  4. Olías, M.; Cánovas, C.R.; Nieto, J.M.; Sarmiento, A.M. Evaluation of the dissolved contaminant load transported by the Tinto and Odiel rivers (South West Spain). Appl. Geochem. 2006, 21, 1733–1749. [Google Scholar] [CrossRef]
  5. Olías, M.; Nieto, J.M.; Sarmiento, A.M.; Cerón, J.C.; Cánovas, C.R. Seasonal water quality variations in a river affected by acid mine drainage: The Odiel River (South West Spain). Sci. Total Environ. 2004, 333, 267–281. [Google Scholar] [CrossRef] [PubMed]
  6. Barry, K.L.; Grout, J.A.; Levings, C.D.; Nidle, B.H.; Piercey, G.E. Impacts of acid mine drainage on juvenile salmonids in an estuary near Britannia Beach in Howe Sound, British Columbia. Can. J. Fish. Aquat. Sci. 2000, 57, 2032–2043. [Google Scholar] [CrossRef]
  7. Lindeström, L.; Tröjbom, M. Konsekvenser för Faluån, Runn och Dalälven av Åtgärder på Gruvavfall i Falun [English: Consequences for the Faluån River, Runn and Dalälven River from Remediation of Mine Waste in Falun]; Faluprojektet: Report 6403; Swedish Environmental Protection Agency: Stockholm, Sweden, 2010; ISBN 978-91-620-6403-7. [Google Scholar]
  8. Widerlund, A.; Öhlander, B.; Ecke, F. Mining and Sustainable Development: Environmental Aspects of Mining; Luleå University of Technology: Luleå, Sweden, 2014. [Google Scholar]
  9. Government Offices of Sweden. Sweden’s Mineral Strategy: For Sustainable Use of Sweden’s Mineral Resources that Creates Growth throughout the Country. Available online: https://www.government.se/49b757/contentassets/78bb6c6324bf43158d7c153ebf2a4611/swedens-minerals-strategy.-for-sustainable-use-of-swedens-mineral-resources-that-creates-growth-throughout-the-country-complete-version (accessed on 20 January 2020).
  10. Callaghan, T.V.; Johansson, M.; Brown, R.D.; Groisman, P.Y.; Labba, N.; Radionov, V.; Bradley, R.S.; Blangy, S.; Bulygina, O.N.; Christensen, T.R.; et al. Multiple Effects of Changes in Arctic Snow Cover. AMBIO 2011, 40, 32–45. [Google Scholar] [CrossRef] [Green Version]
  11. Søndergaard, J.; Halden, N.; Bach, L.; Gustavson, K.; Sonne, C.; Mosbech, A. Otolith Chemistry of Common Sculpins (Myoxocephalus scorpius) in a Mining Polluted Greenlandic Fiord (Black Angel Lead-Zinc Mine, West Greenland). Water. Air. Soil Pollut. 2015, 226, 336. [Google Scholar] [CrossRef]
  12. Alakangas, L.; Öhlander, B.; Lundberg, A. Estimation of temporal changes in oxidation rates of sulphides in copper mine tailings at Laver, Northern Sweden. Sci. Total Environ. 2010, 408, 1386–1392. [Google Scholar] [CrossRef] [PubMed]
  13. Moncur, M.C.; Ptacek, C.J.; Lindsay, M.B.J.; Blowes, D.W.; Jambor, J.L. Long-term mineralogical and geochemical evolution of sulfide mine tailings under a shallow water cover. Appl. Geochem. 2015, 57, 178–193. [Google Scholar] [CrossRef]
  14. Malinovsky, D.; Rodushkin, I.; Moiseenko, T.; Öhlander, B. Aqueous transport and fate of pollutants in mining area: A case study of Khibiny apatite–nepheline mines, the Kola Peninsula, Russia. Environ. Geol. 2002, 43, 172–187. [Google Scholar]
  15. Efimov, V.A.; Chalov, S.R.; Efimova, L.E.; Ivanov, V.A.; Jarsjö, J.; Fischer, S. Impact of mining activities on the surface water quality (case study of Khibiny mountains, Russia). IOP Conf. Ser. Earth Environ. Sci. 2019, 263, 012008. [Google Scholar] [CrossRef]
  16. Leppänen, J.J.; Luoto, T.P.; Weckström, J. Spatio-temporal impact of salinated mine water on Lake Jormasjärvi, Finland. Environ. Pollut. 2019, 247, 1078–1088. [Google Scholar] [CrossRef] [PubMed]
  17. Boas, S.W.; Slotsbo, S.; Patricio Silva, A.L.; Larsen, M.M.; Damgaard, C.; Holmstrup, M. Increased frequency of freeze-thaw events in a future climate can significantly increase negative effects of copper on enchytraeids. Appl. Soil Ecol. 2016, 107, 272–278. [Google Scholar] [CrossRef]
  18. Bergman, H. Gällivare Kommun: Huvudstudierapport Nautanen 2002 [English: Gällivare Municipality: Main Study Report Nautanen 2002]; Technical Report. No. 0223; Gällivare Municipality: Gällivare, Sweden, 2002. (In Swedish) [Google Scholar]
  19. Larborn, L. Inventering av Äldre Gruvavfall i Norrbotten [English: Inventory of Older Mine Waste in Norrbotten]; Report 4; County Administrative Board Norrbotten: Luleå, Sweden, 1993. (In Swedish) [Google Scholar]
  20. Geological Survey of Sweden (SGU) Map Viewer, Ore and Minerals. Available online: https://www.sgu.se/produkter/kartor/kartvisaren/bergkartvisare/malm-och-mineral/ (accessed on 4 February 2019).
  21. Johansson, M. Arkeologisk Dokumentation av Gruvområdet RAÄ 748, Nautanens Kopparfält, Fornvårdsobjekt nr 81 [English: Archeological Documentation of the Mining Area RAÄ 748, Nautanen Copper ore Field, Object nr 81]; Norrbottens Museum: Luleå, Sweden, 2004; Dnr. 429-2004; p. 45. (In Swedish) [Google Scholar]
  22. Hifab, (Envipro). Teknisk Beskrivning Nautanen [English: Technical Description Nautanen]; Technical Report. No. 313 064; Gällivare Municipality: Gällivare, Sweden, 2009. (In Swedish) [Google Scholar]
  23. Jonasson, C.; Grahn, S. Uppföljande Miljökontroll Efterbehandling av Nautanens Gruvområde [English: Follow-Up Monitoring Remediation of Nautanen Mining Area]; Technical Report; Gällivare Municipality: Gällivare, Sweden, 2015. (In Swedish) [Google Scholar]
  24. Sami Parliament. Baste. Available online: https://www.sametinget.se/baste (accessed on 28 January 2019).
  25. Årebäck, H.; Lindblom, M. Samrådsunderlag Angående Ansökan om Bearbetningskoncession Nautanen K nr: 1: Inför Samrådsmöten den 20 Februari 2018 [English: Consultation Material Concerning Application for Exploitation Concession Nautanen K nr: 1: For Consultation February 20th 2018]; Boliden Mineral AB: Stockholm, Sweden, 2018. (In Swedish) [Google Scholar]
  26. Årebäck, H.; Dean, B. Boliden Summary Report: Resources and Reserves 2018: Nautanen; Technical Report; Boliden AB: Stockholm, Sweden, 2018. [Google Scholar]
  27. Harris, I.; Jones, P.D.; Osborn, T.J.; Lister, D.H. Updated high-resolution grids of monthly climatic observations – the CRU TS3.10 Dataset. Int. J. Climatol. 2014, 34, 623–642. [Google Scholar] [CrossRef] [Green Version]
  28. Lynch, E.P.; Jönberger, J.; Bauer, T.E.; Sarlus, Z.; Martinsson, O. Barents Project 2014: Meta-Volcanosedimentary Rocks in the Nautanen Area, Norrbotten: Preliminary Lithological and Deformation Characteristics; Geological Survey of Sweden: Uppsala, Sweden, 2015. (In Swedish) [Google Scholar]
  29. Geological Survey of Sweden (SGU) Map Viewer, Soil Type 1:25000–1:100000. Available online: https://www.sgu.se/produkter/kartor/kartvisaren/jordkartvisare/ (accessed on 4 February 2019).
  30. Herbert, R.; Björkvald, L.; Wällstedt, T.; Johansson, K. Bakgrundshalter av Metaller i Svenska Inlands- och Kustvatten [English: Background Concentrations of Metals in Swedish Inland and Coastal Water]; Report 12; Swedish University of Agricultural Sciences: Uppsala, Sweden, 2009. (In Swedish) [Google Scholar]
  31. Eriksson, L.G. Utläckage av Lakvatten från Gruvområdet - Nautanen [English: Effluents from the Mining Area - Nautanen]; Technical Report. MRAP 94061; County Administrative Board Norrbotten: Luleå, Sweden, 1994. (In Swedish) [Google Scholar]
  32. Nordstrom, D.K. Hydrogeochemical processes governing the origin, transport and fate of major and trace elements from mine wastes and mineralized rock to surface waters. Appl. Geochem. 2011, 26, 1777–1791. [Google Scholar] [CrossRef]
  33. Van der Perk, M. Part II: Sources, role, and behaviour of substances in soil and water. In Soil and Water Contamination: From Molecular to Catchment Scale; Taylor & Francis: Leiden, The Netherlands, 2006; ISBN 0-415-40943-8. [Google Scholar]
  34. Levitan, D.M.; Schreiber, M.E.; Seal, R.R.; Bodnar, R.J.; Aylor, J.G. Developing protocols for geochemical baseline studies: An example from the Coles Hill uranium deposit, Virginia, USA. Appl. Geochem. 2014, 43, 88–100. [Google Scholar] [CrossRef]
  35. Basu, N.B.; Destouni, G.; Jawitz, J.W.; Thompson, S.E.; Loukinova, N.V.; Darracq, A.; Zanardo, S.; Yaeger, M.; Sivapalan, M.; Rinaldo, A.; et al. Nutrient loads exported from managed catchments reveal emergent biogeochemical stationarity. Geophys. Res. Lett. 2010, 37, L23404. [Google Scholar] [CrossRef]
  36. Turc, L. The Water Balance of Soils - Relation between Precipitation, Evaporation and Flow; Annales Agronoiques; Institut National de la Recherche Agronomique: Paris, France, 1954. (In French) [Google Scholar]
  37. Langbein, W.B. Annual Runoff in the United States; US Geol. Surv: Washington, DC, USA, 1949. [Google Scholar]
  38. Walling, D.E.; Fang, D. Recent trends in the suspended sediment loads of the world’s rivers. Glob. Planet. Chang. 2003, 39, 111–126. [Google Scholar] [CrossRef]
  39. Swedish Meteorological and Hydrological Institute (SMHI) Evapotranspiration. Available online: https://www.smhi.se/kunskapsbanken/hydrologi/avdunstning-1.30720 (accessed on 10 July 2019).
  40. Ejhed, H.; Liljeberg, M.; Olshammar, M.; Wallin, M.; Rönnback, P.; Stenström, A. Bruttobelastning på Vatten av Metaller från Punktkällor och Diffusa Källor - Slutrapport [English: Gross Load of Metals in Water from Diffuse and Point Sources]; Report 41; Svenska MiljöEmissionsData: Norrköping, Sweden, 2010; (In Swedish, English Summary). [Google Scholar]
  41. Swedish Environmental Protection Agency. Uppföljning av Efterbehandlingsprojekt Inom Gruvsektorn: Åtgärder, Kostnader och Resultat [English: Survey of Remediation Projects in the Mining Sector: Measures, Costs and results]; Report 5190; Swedish EPA: Stockholm, Sweden, 2002; ISBN 91-620-5190-3. (In Swedish, English Summary). [Google Scholar]
  42. Kimball, B.A.; Runkel, R.L.; Walton-Day, K.; Bencala, K.E. Assessment of metal loads in watersheds affected by acid mine drainage by using tracer injection and synoptic sampling: Cement Creek, Colorado, USA. Appl. Geochem. 2002, 17, 1183–1207. [Google Scholar] [CrossRef]
  43. Palumbo-Roe, B.; Dearden, R. The hyporheic zone composition of a mining-impacted stream: Evidence by multilevel sampling and DGT measurements. Appl. Geochem. 2013, 33, 330–345. [Google Scholar] [CrossRef] [Green Version]
  44. Berglöv, G.; Asp, M.; Berggreen-Clausen, S.; Björck, E.; Axén Mårtensson, J.; Nylén, L.; Ohlsson, A.; Persson, H.; Sjökvist, E. Framtidsklimat i Norrbottens Län - Enligt RCP-Scenarier [English: Future climate in Norrbotten County - Based on RCP Scenarios]; Swedish Meteorological and Hydrological Institute: Norrköping, Sweden, 2015; Klimatologi Nr 32; p. 91. ISBN 1654-2258. (In Swedish) [Google Scholar]
Figure 1. Map of the Nautanen mining site in the Imetjoki drainage basin, where sampling zones for the measurement campaigns in May and August 2017 are marked. Local background samples were taken within the green rectangles (in the upstream part of the Imetjoki drainage basin and in the reference rivers Nietsajoki and Lina River). Samples of the main mining site were taken within the red rectangles, and downstream samples were taken within the blue rectangles.
Figure 1. Map of the Nautanen mining site in the Imetjoki drainage basin, where sampling zones for the measurement campaigns in May and August 2017 are marked. Local background samples were taken within the green rectangles (in the upstream part of the Imetjoki drainage basin and in the reference rivers Nietsajoki and Lina River). Samples of the main mining site were taken within the red rectangles, and downstream samples were taken within the blue rectangles.
Sustainability 12 01394 g001
Figure 2. Photographs showing conditions at the sampling zones: (a) local background (here specifically the upstream lake of the Imetjoki stream); (b) main mining; (c) directly downstream of the main mining zone, D1; and (d) 4 km downstream of Nautanen mining site, D2.
Figure 2. Photographs showing conditions at the sampling zones: (a) local background (here specifically the upstream lake of the Imetjoki stream); (b) main mining; (c) directly downstream of the main mining zone, D1; and (d) 4 km downstream of Nautanen mining site, D2.
Sustainability 12 01394 g002
Figure 3. Boxplots of total log concentrations (µg/L) of selected parameters from Nautanen in May and August 2017 per sampling zone: local background (incl. upstream and reference rivers), mining, and downstream. Whiskers represent 1.5 times the inter-quartile length and black points show values outside that range. Grey areas indicate regional background ranges for northern Sweden (except for Mo, Ba, and SO4 where no regional background values were found). All concentrations are given in Table S3 in the Supplementary Materials.
Figure 3. Boxplots of total log concentrations (µg/L) of selected parameters from Nautanen in May and August 2017 per sampling zone: local background (incl. upstream and reference rivers), mining, and downstream. Whiskers represent 1.5 times the inter-quartile length and black points show values outside that range. Grey areas indicate regional background ranges for northern Sweden (except for Mo, Ba, and SO4 where no regional background values were found). All concentrations are given in Table S3 in the Supplementary Materials.
Sustainability 12 01394 g003
Figure 4. Relative change (in log scale) from upstream to the main mining zone (red color) and between the upstream and downstream sampling zones (blue colors) shown as (a) concentration enrichment factors and (b) as relative mass flow increases at Nautanen in August 2017. Values >1 represent an increase, <1 represents a decrease, and values equal to 1 indicates no change.
Figure 4. Relative change (in log scale) from upstream to the main mining zone (red color) and between the upstream and downstream sampling zones (blue colors) shown as (a) concentration enrichment factors and (b) as relative mass flow increases at Nautanen in August 2017. Values >1 represent an increase, <1 represents a decrease, and values equal to 1 indicates no change.
Sustainability 12 01394 g004
Figure 5. Total Cu concentrations (in log scale µg/L) from water samples from Nautanen taken in 1993–1994, 2001, 2005–2009, 2014, and 2017, grouped by sampling location (local background, mining zone, and downstream locations). All concentrations and measurement locations are outlined in Tables S3 and S4 in the Supplementary Materials.
Figure 5. Total Cu concentrations (in log scale µg/L) from water samples from Nautanen taken in 1993–1994, 2001, 2005–2009, 2014, and 2017, grouped by sampling location (local background, mining zone, and downstream locations). All concentrations and measurement locations are outlined in Tables S3 and S4 in the Supplementary Materials.
Sustainability 12 01394 g005

Share and Cite

MDPI and ACS Style

Fischer, S.; Rosqvist, G.; Chalov, S.R.; Jarsjö, J. Disproportionate Water Quality Impacts from the Century-Old Nautanen Copper Mines, Northern Sweden. Sustainability 2020, 12, 1394. https://doi.org/10.3390/su12041394

AMA Style

Fischer S, Rosqvist G, Chalov SR, Jarsjö J. Disproportionate Water Quality Impacts from the Century-Old Nautanen Copper Mines, Northern Sweden. Sustainability. 2020; 12(4):1394. https://doi.org/10.3390/su12041394

Chicago/Turabian Style

Fischer, Sandra, Gunhild Rosqvist, Sergey R. Chalov, and Jerker Jarsjö. 2020. "Disproportionate Water Quality Impacts from the Century-Old Nautanen Copper Mines, Northern Sweden" Sustainability 12, no. 4: 1394. https://doi.org/10.3390/su12041394

APA Style

Fischer, S., Rosqvist, G., Chalov, S. R., & Jarsjö, J. (2020). Disproportionate Water Quality Impacts from the Century-Old Nautanen Copper Mines, Northern Sweden. Sustainability, 12(4), 1394. https://doi.org/10.3390/su12041394

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop