Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Trace metals in Northern New England streams: Evaluating the role of road salt across broad spatial scales with synoptic snapshots

  • Jessica F. Wilhelm ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing

    wilhelmjf1@gmail.com

    Affiliation Center for the Environment, Plymouth State University, Plymouth, New Hampshire, United States of America

  • Daniel J. Bain,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Geology and Environmental Science, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Mark B. Green,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Center for the Environment, Plymouth State University, Plymouth, New Hampshire, United States of America, Northern Research Station, United States Forest Service, Durham, New Hampshire, United States of America

  • Kathleen F. Bush,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Project administration, Supervision, Validation, Writing – review & editing

    Affiliations Center for the Environment, Plymouth State University, Plymouth, New Hampshire, United States of America, New Hampshire Department of Health and Human Services, Concord, New Hampshire, United States of America

  • William H. McDowell

    Roles Data curation, Formal analysis, Project administration, Resources, Supervision, Validation, Writing – review & editing

    Affiliation Natural Resources and the Environment, University of New Hampshire, Durham, New Hampshire, United States of America

Abstract

Mobilization of trace metals from soils to surface waters can impact both human and ecosystem health. This study resamples a water sample archive to explore the spatial pattern of streamwater total concentrations of arsenic, cadmium, copper, lead, and zinc and their associations with biogeochemical controls in northern New England. Road deicing appears to result in elevated trace metal concentrations, as trace metal concentrations are strongly related to sodium concentrations and are most elevated when the sodium: chloride ratio is near 1.0 (~halite). Our results are consistent with previous laboratory and field studies that indicate cation exchange as a metal mobilization mechanism when road salt is applied to soils containing metals. This study also documents associations among sodium, chloride, dissolved organic carbon, iron, and metal concentrations, suggesting cation exchange mechanisms related to road deicing are not the only mechanisms that increase trace metal concentrations in surface waters. In addition to cation exchange, this study considers dissolved organic carbon complexation and oxidation-reduction conditions affecting metal mobility from soils in a salt-rich environment. These observations demonstrate that road deicing has the potential to increase streamwater trace metal concentrations across broad spatial scales and increase risks to human and ecosystem health.

Introduction

Salinization of surface waters occurs across seasonally snow-covered areas of the world through the practice of road-deicing [1], and has important consequences for water quality [24] and aquatic ecosystems [1,5,6]. Chemical road deicers accumulate in groundwater [7,8], increasing the total concentration of major cations in groundwater throughout the year [9]. The higher concentrations of chloride (Cl) and sodium (Na) cause physiological stresses to aquatic organisms, and thus degrade aquatic ecosystems [4]. Salinization due to deicers also impairs drinking water [10], for example, some groundwater aquifers contain Na concentrations that are unsafe for consumers with high blood pressure [11]. Further, recent work suggests that elevated Cl in water increases corrosivity of drinking water infrastructure and may contribute to elevated trace metals in drinking water (e.g. lead, copper) [12].

Deicer-driven salinization also impacts soils. Many studies have demonstrated the impact of salinization on soil cation exchange dynamics [13,14], which can mobilize trace metals from soils [15,16]. Soils with sorbed trace metals–either natural concentrations or concentrations enhanced due to human activities–are susceptible to altered ion exchange processes and mobilization [17]. Metal mobility is also enhanced by interactions among deicers, dissolved organic matter, and colloidal materials [15,17] (Fig 1).

thumbnail
Fig 1. Associations among sodium, chloride, dissolved organic carbon, iron, and metal concentrations.

This study considers cation exchange, dissolved organic carbon complexation, and oxidation-reduction conditions affecting metal mobility from soils in a salt-rich environment.

https://doi.org/10.1371/journal.pone.0212011.g001

Soils serve as an important regulator of trace metal dynamics, storing a range of potentially toxic metals. High concentrations of cations, particularly Na introduced during deicing, can displace metals via cation exchange and interact with dissolved organic carbon (DOC) and colloids to make soil metals available for potential transport to waters further downstream. In particular, Amrhein et al. [15] assessed trace metal mobility in roadside soils from cold weather climates using soil columns leached with road salt chemicals. They found metal concentrations increased during flushing by fresh water following exposure to NaCl. This simulated snowmelt seemed to mobilize low molecular weight DOC and bound metals. The fresh water flush contained loads of trace metals that were greater than those in the experimental high salt flush.

The vast majority of deicer salinization and soil studies have been conducted either in the laboratory [15,18] or at limited spatial scales [17,1921], and the impact of these salinization processes on regional streamwater quality remains poorly understood. This study documents total trace metal concentrations across a large spatial scale in Northern New England, using an archive of synoptically sampled surface waters, and compares the spatial variability of trace metal concentrations with indicators of chemical deicers (Fig 2). Given the archival nature of the samples, the impacts of storage conditions are assessed as part of this analysis. Strong co-variance between the amount of urbanization and chemical indicators of road salt has been documented in this region previously [9,22]. Therefore this study focuses on soil-streamwater processes in areas impacted by road salt.

thumbnail
Fig 2. Metal concentrations during periods of fluctuating halite levels.

The sites representing the tail denote metal concentrations when higher Na+ is released from exchange sites, relative to lower Cl- flushed following road salt inputs. Corresponding Cl- is quickly flushed through the system following input, and therefore absent in periods not dominated by halite. The sites near the 1:1 halite peak, i.e. equal Na+ and Cl-, suggest high dissolved metal flux associated with periods of high halite inputs.

https://doi.org/10.1371/journal.pone.0212011.g002

Methods

In 2013 approximately 70 sites in northern New England were sampled synoptically in May, June, and July (sampling and storage methods detailed below). All sampling sites are located at road crossings of flowing waters. By New Hampshire State law, all navigable waters (broadly defined) are public spaces (Concord Mfg. Co. v. Robertson, 25 A. 718, 720 (N.H. 1890)) and these public spaces are accessible via road right of ways (Hartford v. Town of Gilmanton, 101 N.H. 424 (N.H. 1959)). Therefore, special collection permissions were not required for any of the locations. Further, this study did not involve endangered or protected species. We evaluated concentrations of As, Cd, Cu, Pb, and Zn in the archived samples and compared these concentrations with road salt indicators (Na, Cl, and Na: Cl), water quality parameters (DOC), and redox indicators, iron (Fe).

Site distribution

Sites were located in New Hampshire and northern Massachusetts and associated with the citizen science Lotic Volunteer Temperature, Electrical Conductivity, and Stage (LoVoTECS) sensing network [23] (Fig 3). Sites spanned a range of conditions including basin size, elevation, and land cover. Catchments ranged from 1.5 km2 to 2646 km2. Sample collection sites included four major New England rivers and their tributaries: the upper and lower Androscoggin, the Saco, the Merrimack, and the upper Connecticut. Study area climate is temperate; average July highs ranging between 24°C and 28°C, and average January lows ranging between -20°C and -9°C. Mean annual precipitation ranges between 40 cm and 118 cm [24]. The regional elevation ranges from 0 (MSL) to the highest point, Mount Washington, at 1,917 meters.

thumbnail
Fig 3. Site map.

Sites in Massachusetts and New Hampshire evaluated for metal concentrations and biogeochemistry.

https://doi.org/10.1371/journal.pone.0212011.g003

The bedrock geology of the greater New Hampshire region can be broadly separated into three composite terranes: the Grenville, Central Maine, and Nashoba-Casco-Miramichi [25]. Metasedimentary materials, found particularly in the southwestern part of the state, dominate most of New Hampshire. To the west, calc-silicate metasediments, and volcanic rocks and granitic gneisses are also important bed rock types [25].

Water sampling

Three synoptic sampling events were performed May 14th (n = 61), June 11th (n = 66), and July 16th (n = 77) 2013. Although most sites were sampled on all three-sample collection dates, not all sites were sampled on all three dates. This limits our ability to interpret temporal patterns. All samples were collected within a 12-hour period on the sampling date. Grab samples were collected 1 to 2 meters from the stream or riverbank. One-liter polypropylene sample vials were rinsed three times with filtered stream water and 30 mL borosilicate glass vials were rinsed three times with stream water before sample collection. One-liter polypropylene samples were field filtered and used for major ion and DOC analyses. The samples collected in 30 mL borosilicate glass vials were intended for water isotope analysis and therefore were neither filtered nor acidified at the time of collection. Immediately following collection, the samples were placed on ice, and stored at 4°C upon return to the laboratory.

There was some hydrologic variability among sampling events (Fig 4). Mean daily runoff for four USGS gages spanning our sampling area was 3.6 mm/d in May (flows were elevated due to residual snow melt), 4.2 mm/d in June (collected between two storms), and 2.1 mm/d in July (low flow conditions).

thumbnail
Fig 4. 2013 hydrograph comparisons at four sites spanning our sampling domain.

River discharge measurements from the Ammonoosuc (north central New Hampshire; USGS gage 01137500), Ashuelot (southwestern New Hampshire; USGS gage 01157000), Pemigewassett (central New Hampshire; USGS gage 01076500), and Lamprey (southeast New Hampshire; USGS gage 01064500). Dashed vertical lines indicate synoptic sampling events.

https://doi.org/10.1371/journal.pone.0212011.g004

Analytical methods

Dissolved organic carbon was measured with a Shimadzu TOC analyzer, and major ions (Na, K, Mg, Ca, Cl) were measured on a Dionex X2 ion chromatograph at the University of New Hampshire in 2013. These water chemistry measurements were made using sub-samples from the polypropylene bottles that were filtered and held frozen prior to analysis. Metals and metalloid concentrations (hereafter named “metals”) were measured with an Inductively-Coupled Plasma-Mass Spectrometer (ICP-MS, Perkin Elmer NexION 300X) at the University of Pittsburgh in February 2014, 8–10 months after collection. Metal concentrations were determined from sub samples of the synoptic water isotope archive described above. Sub-samples for metals analysis from the glass vials were diluted with sub-boil distilled 2% nitric acid (HNO3) to standardize matrix composition and spiked with an internal standard of Be, Ge, and Tl. The ICP was calibrated with a five-point curve and blanks and drift checks measured roughly every 10 samples. For elements with strong polyatomic interferences (e.g., As) kinetic energy discrimination collision cell methods were used.

These data were produced from a sample archive stored under conditions that do not follow typical conventions for dissolved metal preservation (samples were not filtered, not acidified, and stored in borosilicate vials). These samples should not be compared directly with either total or dissolved metal concentrations. To evaluate the impact of these storage conditions, a set of six samples were collected in July 2017 from sample sites that spanned the range of Na concentration and watershed land cover across our sampling domain, and collected and stored using both methods (filtered, acidified, polypropylene container and not filtered, not acidified, borosilicate glass container) (Fig 5). These samples were allowed to sit for four months to allow equilibration. Then metal concentrations in the two sets of samples were evaluated on the ICP-MS as described above and compared.

thumbnail
Fig 5. Comparison of metal concentration measurements on unfiltered, unacidified samples.

Samples were stored in borosilicate glass (total; x-axes) and filtered, acidified samples stored in polypropylene vials (dissolved; y-axes). Zinc is not shown because it was not detectable in unfiltered and unacidified samples (e.g. bdc, hbk, hod, nwd, pbb, prp), and filtered and acid-preserved samples (e.g. hbk, hod, pbb, prp), and therefore comparisons could not be made between total and dissolved zinc concentrations.

https://doi.org/10.1371/journal.pone.0212011.g005

Statistical methods

Data did not uniformly meet normality criteria and therefore relationships between trace metals and other biogeochemical parameters (e.g., DOC and Fe concentrations) were evaluated using Spearman’s rank order correlation with alpha = 0.05. All statistics were performed in R Computing Software [26].

Results

Direct comparison of preservation methods demonstrates that the metal concentrations measured in the archived samples are not precisely comparable to the filtered and acidified samples (Fig 5). However, in general, the signal strengths are related positively, though substantial noise is introduced. Therefore, the measurements presented throughout this paper seem to reflect overarching trends in synoptic metal concentrations, albeit a relatively noisier version of that signal.

Quality assurance procedures (e.g., internal standard values) identified three of the synoptic samples that were clear outliers, and these samples were removed from all subsequent data analysis. In the remaining 191 samples, metal concentrations varied widely across the study area (Table 1). Despite the variability, median concentrations were relatively consistent across the three summer months. High concentrations of As and Pb (max As = 18.1 μg/L, max Pb = 2.77 μg/L) occurred in samples with elevated Cl and Na concentrations (Fig 6). Further, the highest concentrations of As, Cd, and Pb occurred in samples with molar Na: Cl ratios that approached 1.0 (Figs 6 and 7). Copper, Pb, and sometimes As concentrations were associated with samples containing elevated DOC concentrations (Figs 8 and 9). Samples with higher concentrations of As were associated with samples with higher concentrations of Fe.

thumbnail
Fig 6. Arsenic and lead relationships with Na and Cl.

Relationship between As and Pb concentrations and Na and Cl concentrations and Na: Cl molar ratios in synoptic samples collected May 14, June 11, and July 16, 2013. Na:Cl molar ratio is colored using a gradient from high Cl concentrations (light pink), to low Cl concentrations (dark purple).

https://doi.org/10.1371/journal.pone.0212011.g006

thumbnail
Fig 7. Cadmium, copper, and zinc relationships with Na and Cl.

Relationships among Cd, Cu, and Zn concentrations and Na and Cl concentrations and Na: Cl molar ratios in synoptic samples collected May 14, June 11, and July 16, 2013. Na:Cl molar ratio is colored using a gradient from high Cl concentrations (light pink), to low Cl concentrations (dark purple).

https://doi.org/10.1371/journal.pone.0212011.g007

thumbnail
Fig 8. Arsenic and lead relationships with DOC and Fe.

Relationships between As and Pb concentrations and dissolved organic carbon (DOC) and Fe concentrations during May 14, June 11, and July 16, 2013.

https://doi.org/10.1371/journal.pone.0212011.g008

thumbnail
Fig 9. Cadmium, copper, and zinc relationships with DOC and Fe.

Relationships among Cd, Cu, and Zn concentrations and dissolved organic carbon (DOC) and Fe concentrations during May 14, June 11, and July 16, 2013.

https://doi.org/10.1371/journal.pone.0212011.g009

Elemental relationships

Metal concentration associations with DOC, Na and Cl, and Fe are consistent across months (Table 2). Over the course of the summer, arsenic and lead have the strongest relationships with all three indicators of chemical processes, i.e., cation exchange (Na, Cl) (Fig 6), dissolved organic carbon complexation (DOC), and oxidation-reduction reactions (Fe) (Fig 8).

thumbnail
Table 2. Evaluation of elemental relationships using Spearman’s rank order correlation.

https://doi.org/10.1371/journal.pone.0212011.t002

Arsenic.

Arsenic relationships with other metals and DOC are relatively consistent across months, and the strongest associations are observed in July (Figs 6 and 8, Table 2). Arsenic is strongly related to Na and Cl throughout the summer, and the highest concentrations of As occur as the molar Na: Cl ratio approaches one, the characteristic value of halite (Fig 6).

Cadmium.

Cadmium has no apparent relationships with Fe, DOC (Fig 8), nor with Na and Cl concentrations (Fig 7). Elevated cadmium concentrations occur as the Na: Cl ratio approaches one (Fig 7).

Copper.

Copper relationships with DOC and Fe are similar across months (Fig 9), and Cu concentrations are associated with Na and Cl (Fig 7). Copper concentrations are elevated with Na: Cl ratio values near 1.

Lead.

Lead is associated with DOC and Fe throughout the summer months (Fig 8), and Na and Cl concentrations in May, June, and July (Fig 6). Lead concentrations are highest in samples with Na: Cl ratios that approach the 1.0 molar ratio.

Zinc.

Zinc does not seem to correlate with Fe or DOC (Fig 9). Zinc concentrations are not related to Na or Cl concentrations (Fig 7).

Dissolved organic carbon, iron, sodium chloride.

Dissolved organic carbon is associated with Na throughout the summer months, and with Cl (Fig 10). Slight increases in DOC concentrations occur near 1.0 Na: Cl. Dissolved organic carbon is strongly correlated with Fe throughout May, June, and July. Iron is strongly associated with Na and Cl.

thumbnail
Fig 10. Na and Cl relationships with DOC.

Concentrations of Na, Cl, and Na: Cl molar ratio plotted against dissolved organic carbon (DOC) concentrations for all locations throughout northern New England on May 14, June 11, and July 16, 2013. Na:Cl molar ratio is colored using a gradient from high Cl concentrations (light pink), to low Cl concentrations (dark purple).

https://doi.org/10.1371/journal.pone.0212011.g010

Discussion

It is difficult to directly compare our water chemistry measurements with dissolved concentration data commonly reported in the literature. In this study, trace metal concentrations were measured in subsamples of an unfiltered, unacidified archive that was stored in borosilicate containers. This raises the potential for 1) the measurement of metals bound to particles (i.e., larger than 0.45 micron) that would be removed by filtration; 2) sorption of dissolved metals to the borosilicate surfaces; and 3) chemical precipitation processes and the formation of particulates that sorb dissolved trace metals (e.g., iron oxidation and the formation of iron hydroxides). Our direct comparison of samples subject to both the conditions for the synoptic sample archive and a standard dissolved metal measurement suggest that a mixture of these potential processes likely occurred, with concentrations in the synoptic archive measuring both above and below the concentrations measured in “dissolved metal” samples. Nonetheless, the measurement of the archived samples captures interesting patterns across broad spatial scales; although these patterns are likely more noisy than those that would be observed under optimal preservation methods, the limitations of archive conditions do not preclude assessment of spatial trends in metals.

Our data are consistent with road salt mobilization of trace metals and the resultant elevated streamwater concentrations of these metals. Observed Na and Cl concentrations range from highly dilute (min Na = 1800 μg/L; min Cl = 1320 μg/L) to highly concentrated (max Na = 160000 μg/L; max Cl = 420000 μg/L). The strong, positive correlation between trace metal concentrations (Pb, As, Cu, and Zn) and both Na and Cl indicate trace metal concentrations are elevated in salty waters. Further, the elevated concentrations of Pb, As, and Cd in samples where dissolved Na: Cl ratios are near 1:1 (i.e., the ratio in halite) provide a parallel and consistent line of evidence that road salt mobilizes trace metals. The Na:Cl ratio can vary in surface waters both above 1 (during periods after salting when lagged Na transport occurs) and below 1 (during the early salt flux where Na is sorbed and Cl is transported rapidly). In addition, other common sources of Na and Cl are not 1:1. For example, rainfall in experimental catchments is ~0.8 [27,28]. We do not have the data to claim causality, but the clear tendency warrants continued attention. While trace metals are mobilized, the concentrations are relatively low, and streamwater metal concentrations were always low relative to US EPA drinking water standards (Table 3). Some sites with a strong road salt signature had metal concentrations approaching levels of concern with concentrations above EPA freshwater chronic criteria for aquatic organisms measured at several stations. Comparisons of our measured trace metal values with water quality standards should be made cautiously due to the archival nature of the samples.

Previous studies of trace metal mobilization to surface waters by road salt have focused on limited spatial scales. Bäckström et al [16] studied two transects located at most 30m from the roadway, Löfgren [20] studied five small forested catchments (drainage areas 57–168 hectares), and Cunningham et al. [30] collected samples within the limits of a single upstate NY city (76 hectare college campus and one additional off campus site). Our sites span tens of thousands of km2, and allow evaluation of the road salt signal at much broader scales. Across these sites, and during each month, trace metals were consistently related to both Na and Cl. We interpret the strength of these relationships across a large scale, and their persistence across months, to indicate a general mechanism or set of mechanisms related to road salt driven salinization. Our data cannot reveal the patterns of hydrologic transport that would move these trace metals from soil pores to stream water. However, these mechanisms seem to alter trace metal concentrations regardless of heterogeneity in hydrologic transport processes across wide spatial scales.

Further, our data indicate that trace metals are mobilized by processes that parallel or interact with road salt chemical fluxes, processes that are consistent across broad spatial scales. Most importantly, As and Pb are related to both Fe and DOC concentrations across all months sampled. These associations suggest a number of potential mechanisms. First, Mn and Fe oxides in soils and sediments strongly sorb many cations [31]. In particular, As has strong adsorption affinities for reduced Fe species, Fe (III) oxyhydroxides, and can be released to groundwater during Fe-reducing biodegradation [32]. When these oxides are reduced, any sorbed cations are mobilized [33,34]. Therefore, any flushing of reducing zones by waters rich in road salt could mobilize dissolved metals in those waters. Alternatively, the association between Pb/As and DOC suggests that road salt chemical fluxes can enhance the DOC flux and therefore metals associated with the DOC. For example, Amrhein et al. [15,18] observed NaCl mobilizes metals associated with organic matter and colloids. However, ultimately, the covariance of Fe and DOC suggests increases in Fe could also simply result from increased amounts of DOC-sorbed Fe during increased DOC flux that arises from salting. Increase in soil salinity can cause DOC to flocculate with cations (Ca2+ and Mg2+) [35], enhancing metal mobility from soils to streams. Therefore, the NaCl interactions with DOC can potentially explain both associations (Pb/As with Fe and Pb/As with DOC) and reduction/oxidation dynamics may not be an important control. Our results do not allow us to partition the relative importance of salt pulses, redox conditions, and organic carbon to metal mobilization, though the consistent associations among salt, Pb/As, DOC, and Fe suggests these interactions should be more closely examined as inquiry into road salt effects continues.

It also should be noted that elevated As has been documented in groundwaters across the southern part of our sampling domain. These concentrations are attributed partially to naturally occurring As associated with weathering of bedrock minerals [36,37]. While our data demonstrate a relationship with road salt indicators, we cannot quantify the potential contribution of naturally occurring As. This is particularly challenging as areas rich in natural As sources are also the most urban. Deicer impacts covary spatially with bedrock As hotspots and observed associations could result from this covariance. However, the other trace metals we measured (Pb, Cu, Zn, and Cd) do not have major, documented natural sources in this region. Given the consistency in the observed relationships between these metals and deicer impacts, this covariance driving the As signal alone seems less likely. Potential anthropogenic sources of Pb, Cu, Zn, Cd, and As include road materials and mobile source emissions [38] and general trace metal pollution sources [37,38]. In particular, As and Pb enrichment associated with historical use of lead arsenate as a pesticide in orchards has been observed in study area streamwaters [39]. The bottom line is that there are many potential sources of trace metal contamination, however, road deicers and associated chemistry seem to mobilize the metals regardless of source.

Conclusions

The patterns in our data highlight several important aspects of surface water salinization. Persistent patterns in stream water chemistry across broad spatial scales suggest New England flowing waters are influenced by pools with elevated Na and Cl concentrations (e.g., soil and ground water). Soils, where most of the mechanistic work on metals mobilization has been conducted, are only ephemerally connected hydrologically to streams—primarily during storm events. In contrast, groundwater contributes a substantial fraction of streamwater flow, especially during low flows (i.e., our sampled conditions). Detailed hillslope studies to trace the movement of mobilized trace metals in soils [e.g., 21], and to determine whether they are laterally transported to drainage systems or streams, or vertically percolated into groundwater can clarify the relative importance of the pool. Second, the observed broad spatial interactions among NaCl, trace metals, Fe, and DOC indicate that documentation of road salt impacts should examine a broader range of chemical parameters. Without a broader perspective, mechanism can be misattributed, and therefore, mitigation based on potential misunderstandings will be less effective.

Supporting information

S1 File.

Table A. Streamwater and river collected ion site-specific information. Water sample data collected across northern New England and analyzed for metals and water quality; sites marked with (*) include those samples removed from subsequent statistical analysis. Samples indicated as ‘bdl’ are due to concentrations below detection limits. Table B. Land use and landscape site-specific information. Site characteristics across northern New England from the New Hampshire Land Cover archive (NLCD 2011); sites marked with (*) include those samples removed from subsequent statistical analysis. Table C. Contrast dissolved and total sample site-specific information. Resampled and analyzed ion data comprising a subset of our sites in northern New England. Samples indicated as ‘bdl’ are due to concentrations below detection limits.

https://doi.org/10.1371/journal.pone.0212011.s001

(DOC)

Acknowledgments

Support for this study was provided by the National Science Foundation, Plymouth State University Center for the Environment, University of Pittsburgh Department of Geology and Environmental Science, New Hampshire Established Program to Stimulate Competitive Research (NSF EPS-1101245), and the U.S. Forest Service. We especially thank the citizen scientist research group, Lotic Volunteer Temperature, Electric Conductivity, and Stage Sensing Network (LoVoTECS), and field technicians Ashley Inserillo and Errin Volitis, who coordinated the collection of the considerable data across Northern New England rivers and streams.

References

  1. 1. Kelly VR, Lovett GM, Weathers KC, Findlay SEG, Strayer DL, Burns DJ, et al. Long-Term Sodium Chloride Retention in a Rural Watershed: Legacy Effects of Road Salt on Streamwater Concentration. Environ Sci Technol. 2008;42: 410–415. pmid:18284139
  2. 2. Corsi SR, De Cicco LA, Lutz MA, Hirsch RM. River chloride trends in snow-affected urban watersheds: increasing concentrations outpace urban growth rate and are common among all seasons. Sci Total Environ. 2015;508: 488–497. pmid:25514764
  3. 3. Kaushal S, McDowell W, Wollheim W, Johnson T, Mayer P, Belt K, et al. Urban Evolution: The Role of Water. Water. 2015;7: 4063–4087.
  4. 4. Williams DD, Williams NE, Cao Y. Road salt contamination of groundwater in a major metropolitan area and development of a biological index to monitor its impact. Water Res. 2000;34: 127–138.
  5. 5. Hanes RE, Zelazny LW, Verghese KG, Bosshart RP, Carson EW, Blaser RE, et al. EFFECTS OF DEICING SALTS ON PLANT BIOTA AND SOIL—EXPERIMENTAL PHASE. NCHRP Rep. 1976; Available: https://trid.trb.org/view.aspx?id=66581
  6. 6. Corsi SR, Graczyk DJ, Geis SW, Booth NL, Richards KD. A Fresh Look at Road Salt: Aquatic Toxicity and Water-Quality Impacts on Local, Regional, and National Scales. Environ Sci Technol. 2010;44: 7376–7382. pmid:20806974
  7. 7. Findlay SEG, Kelly VR. Emerging indirect and long-term road salt effects on ecosystems: Findlay & Kelly. Ann N Y Acad Sci. 2011;1223: 58–68. pmid:21449965
  8. 8. Ostendorf DW, Palmer RN, Hinlein ES. Seasonally varying highway de-icing agent contamination in a groundwater plume from an infiltration basin. Hydrol Res. 2009;40: 520–532.
  9. 9. Daley ML, Potter JD, McDowell WH. Salinization of urbanizing New Hampshire streams and groundwater: effects of road salt and hydrologic variability. J North Am Benthol Soc. 2009;28: 929–940.
  10. 10. Jackson RB, Jobbágy EG. From icy roads to salty streams. Proc Natl Acad Sci U S A. 2005;102: 14487–14488. pmid:16203970
  11. 11. Talukder MRR, Rutherford S, Huang C, Phung D, Islam MZ, Chu C. Drinking water salinity and risk of hypertension: A systematic review and meta-analysis. Arch Environ Occup Health. 2016; 1–13.
  12. 12. Stets EG, Lee CJ, Lytle DA, Schock MR. Increasing chloride in rivers of the conterminous U.S. and linkages to potential corrosivity and lead action level exceedances in drinking water. Sci Total Environ. 2018;613–614: 1498–1509. pmid:28797521
  13. 13. Mason CF, Norton SA, Fernandez IJ, Katz LE. Deconstruction of the chemical effects of road salt on stream water chemistry. J Environ Qual. 1999;28: 82–91.
  14. 14. Shanley JB. Effects of Ion Exchange on Stream Solute Fluxes in a Basin Receiving Highway Deicing Salts. Journal Environmental Quality. 1994;23: 977–986.
  15. 15. Amrhein C, Strong JE, Mosher PA. Effect of deicing salts on metal and organic matter mobilization in roadside soils. Environ Sci Technol. 1992;26: 703–709.
  16. 16. Bäckström M, Karlsson S, Bäckman L, Folkeson L, Lind B. Mobilisation of heavy metals by deicing salts in a roadside environment. Water Res. 2004;38: 720–732. pmid:14723942
  17. 17. Norrström A-C, Bergstedt E. The impact of road de-icing salts (NaCl) on colloid dispersion and base cation pools in roadside soils. Water Air Soil Pollut. 2001;127: 281–299.
  18. 18. Amrhein C, Strong JE. The effect of deicing salts on trace metal mobility in roadside soils. J Environ Qual. 1990;19: 765–772.
  19. 19. Bäckström M, Karlsson S, Bäckman L, Folkeson L, Lind B. Mobilisation of heavy metals by deicing salts in a roadside environment. Water Res. 2004;38: 720–732. pmid:14723942
  20. 20. Löfgren S. The Chemical Effects of Deicing Salt on Soil and Stream Water of Five Catchments in Southeast Sweden. Water Air Soil Pollut. 2001;130: 863–868.
  21. 21. Rossi RJ, Bain DJ, Elliott EM, Divers M, O’Neill B. Hillslope soil water flowpaths and the dynamics of roadside soil cation pools influenced by road deicers: Soil Water Flowpaths and Interactions of Roadside Soils and Road Salt. Hydrol Process. 2017;31: 177–190.
  22. 22. Zuidema S, Wollheim WM, Mineau MM, Green MB, Stewart RJ. Controls of Chloride Loading and Impairment at the River Network Scale in New England. J Environ Qual. 2018;47: 839. pmid:30025050
  23. 23. Inserillo EA, Green MB, Shanley JB, Boyer JN. Comparing catchment hydrologic response to a regional storm using specific conductivity sensors. Hydrol Process. 2017;31: 1074–1085.
  24. 24. New Hampshire State Climate Office WebSite [Internet]. 2014 [cited 17 Apr 2018]. Available: http://unh.edu/stateclimatologist/
  25. 25. Peters SC, Blum JD, Klaue B, Karagas MR. Arsenic Occurrence in New Hampshire Drinking Water. Environ Sci Technol. 1999;33: 1328–1333.
  26. 26. Team RC. R: A language and environment for statistical computing [Internet]. Vienna, Austria; 2014. 2017.
  27. 27. Lovett GM, Likens GE, Buso DC, Driscoll CT, Bailey SW. The biogeochemistry of chlorine at Hubbard Brook, New Hampshire, USA. Biogeochemistry. 2005;72: 191–232.
  28. 28. Neal C, Kirchner JW, others. Sodium and chloride levels in rainfall, mist, streamwater and groundwater at the Plynlimon catchments, mid-Wales: inferences on hydrological and chemical controls. Hydrol Earth Syst Sci Discuss. 2000;4: 295–310.
  29. 29. US EPA O. National Recommended Water Quality Criteria—Aquatic Life Criteria Table. In: US EPA [Internet]. 3 Sep 2015 [cited 17 Apr 2018]. Available: https://www.epa.gov/wqc/national-recommended-water-quality-criteria-aquatic-life-criteria-table
  30. 30. Cunningham MA, Snyder E, Yonkin D, Ross M, Elsen T. Accumulation of deicing salts in soils in an urban environment. Urban Ecosyst. 2008;11: 17–31.
  31. 31. Chao TT. Selective Dissolution of Manganese Oxides from Soils and Sediments with Acidified Hydroxylamine Hydrochloride. Soil Sci Soc Am J. 1972;36: 764–768.
  32. 32. Ziegler BA, McGuire JT, Cozzarelli IM. Rates of As and Trace-Element Mobilization Caused by Fe Reduction in Mixed BTEX–Ethanol Experimental Plumes. Environ Sci Technol. 2015;49: 13179–13189. pmid:26486694
  33. 33. Lovely DR, Phillips EJ, Lonergan DJ. Enzymatic versus nonenzymatic mechanisms for Fe (III) reduction in aquatic sediments. Environ Sci Technol States. 1991;25. Available: http://www.osti.gov/scitech/biblio/5118160
  34. 34. Smedley PL, Kinniburgh DG. A review of the source, behaviour and distribution of arsenic in natural waters. Appl Geochem. 2002;17: 517–568.
  35. 35. Fram MS, Fujii R, Weishaar JL, Bergamaschi BA, Aiken GR. How DOC composition may explain the poor correlation between specific trihalomethane formation potential and specific UV absorbance. Water Resour Investig Rep WRI. 1999; 99–4018.
  36. 36. Boudette EL, Canney FC, Cotton JE, Davis RI, Ficklin H, Motooka JM. High levels of arsenic in the groundwaters of southeastern New Hampshire; a geochemical reconnaissance [Internet]. U.S. Geological Survey,; 1985. Report No.: 85–202. Available: http://pubs.er.usgs.gov/publication/ofr85202
  37. 37. Peters SC, Blum JD. The source and transport of arsenic in a bedrock aquifer, New Hampshire, USA. Appl Geochem. 2003;18: 1773–1787.
  38. 38. Zehetner F, Rosenfellner U, Mentler A, Gerzabek MH. Distribution of Road Salt Residues, Heavy Metals and Polycyclic Aromatic Hydrocarbons across a Highway-Forest Interface. Water Air Soil Pollut. 2009;198: 125–132.
  39. 39. Peryea FJ, Kammereck R. Phosphate-Enhanced Movement of Arsenic Out of Lead Arsenate-Contaminated Topsoil and Through Uncontaminated Subsoil. Water Air Soil Pollut. 1997;93: 243–254.