Andrea Burke

Theodore M Present

Guillaume Paris

Emily C M Rae

Brodie H Sandilands

Jérôme Gaillardet

Bernhard Peucker-Ehrenbrink

Woodward W Fischer

James W Mcclelland

Robert G M Spencer

Britta M Voss

Jess F Adkins

Sulfur isotopes in rivers: Insights into global weathering budgets, pyrite oxidation, and the modern sulfur cycle

Keywords: sulfur, rivers, weathering, pyrite

come   L'archive ouverte pluridisciplinaire

Introduction

The biogeochemical sulfur cycle is intimately linked to the cycles of carbon and oxygen (e.g. [START_REF] Berner | Burial of Organic-Carbon and Pyrite Sulfur in Sediments Over Phanerozoic Time -a New Theory[END_REF]. Reconstructing the sources and sinks of sulfur to the marine environment in the past is thus important for understanding long-term changes in climate and the redox processes operating in Earth's surface environments. The sulfur isotope compositions of sources and sinks provide a sensitive tracer of the processes that drive the sulfur cycle because there are large isotope fractionations that occur associated with cycling sulfur between oxidized and reduced phases (e.g. [START_REF] Garrels | Coupling of the sedimentary sulfur and carbon cycles; an improved model[END_REF].

Microbial sulfate reduction, for instance, imparts a large S isotope fractionation (ε≈0 to -70‰ (e.g. [START_REF] Habicht | Isotope fractionation by sulfate-reducing natural populations and the isotopic composition of sulfide in marine sediments[END_REF][START_REF] Sim | Effect of electron donors on the fractionation of sulfur isotopes by a marine Desulfovibrio sp[END_REF])), leaving-on average-pyrite and other sulfide-bearing minerals with lower sulfur isotope ratios than seawater and sedimentary sulfate.

Reconstructions of sulfur isotope ratios through geologic time from marine sedimentary rocks have typically been used to infer past changes in the burial flux of reduced sulfur (pyrite) relative to the removal of oxidized sulfur in the form of sulfate (evaporite deposits) [START_REF] Kump | Modeling Atmospheric O-2 in the Global Sedimentary Redox Cycle[END_REF]. In a simple isotope box model of the marine sulfur reservoir, variations in the isotopic composition of marine sulfate are interpreted as being driven by relative changes in these outputs, while typically assuming that the input of sulfur to the ocean has remained constant through time. Recent work [START_REF] Halevy | Sulfate Burial Constraints on the Phanerozoic Sulfur Cycle[END_REF] , however, has highlighted the need to consider changes in the flux and the isotopic composition of sulfur to the ocean. Riverine sulfur is the major source of sulfate to the ocean, supplying approximately 4.7 Tmol/y today (including 1.3 Tmol/y from anthropogenic sources [START_REF] Meybeck | Global Occurrence of Major Elements in Rivers[END_REF]). Thus in order to fully understand the secular changes in δ 34 S of seawater, we need to better constrain both the modern values for, and controls on, the isotopic composition of riverine sulfate.

The modern riverine sulfur isotopic composition can also inform estimates of chemical weathering fluxes, with important implications for the carbon cycle. Sulfur isotopes in rivers can provide insight into how much riverine S is sourced from dissolution of sedimentary sulfate minerals versus oxidative weathering of pyrite (OWP) [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF]. OWP produces sulfuric acid, which is a source of acidity for chemical weathering and which, when it interacts with carbonate minerals, can lead to a net release of CO2, in contrast to the sink of CO2 associated with silicate weathering (e.g. [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF][START_REF] Torres | The acid and alkalinity budgets of weathering in the Andes-Amazon system: Insights into the erosional control of global biogeochemical cycles[END_REF]2015). Previous estimates of global OWP fluxes range from 0.5 to 0.65 Tmol/y [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF][START_REF] Francois | Modelling the Phanerozoic carbon cycle and climate; constraints from the 87 Sr/ 86 Sr isotopic ratio of seawater[END_REF][START_REF] Lerman | CO2 and H2SO4 consumption in weathering and material transport to the ocean, and their role in the global carbon balance[END_REF]. However, recent studies that use sulfur isotopes (and sulfate-oxygen isotopes) from individual catchments indicate that estimates of global OWP flux are potentially much too low. The sum of OWP fluxes (0.15 Tmol/y) from just three river basins (Mackenzie [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF], Kaoping [START_REF] Das | Disproportionately high rates of sulfide oxidation from mountainous river basins of Taiwan orogeny: Sulfur isotope evidence[END_REF], and Jialing [START_REF] Li | Identification of dissolved sulfate sources and the role of sulfuric acid in carbonate weathering using dualisotopic data from the Jialing River, Southwest China[END_REF]) can account for a third of previous global OWP flux estimates, despite covering less than 2% of global land area [START_REF] Das | Disproportionately high rates of sulfide oxidation from mountainous river basins of Taiwan orogeny: Sulfur isotope evidence[END_REF]. Underestimating global OWP by this magnitude may result in substantial overestimates of the modern-day sink of CO2 associated with chemical weathering.

Previous estimates of δ 34 S of river water

Previous estimates of the global sulfur isotopic composition of rivers come from either measurements of river water from a single geographical region [START_REF] Ivanov | The Sulphur Cycle in Continental Reservoirs[END_REF] or back-of-the-envelope calculations based on simple geochemical assumptions [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF]. The previous data-based study that included the largest amount of river data was limited to the Eurasian continent, and reported an average riverine δ 34 S of 9.2‰ [START_REF] Ivanov | The Sulphur Cycle in Continental Reservoirs[END_REF].

The rivers sampled represent only 7% of the total global riverine discharge and have a total sulfate flux of 0.4 Tmol/y, accounting for only 9% of the total riverine sulfate flux. The limited geographic extent of this estimate raises the question of how representative the value of 9.2‰ is for the global riverine δ 34 S input to the oceans, especially given that many of the rivers sampled are weathering large evaporitic deposits of Cambrian/Ordovician age that are exposed across the Siberian Platform (Huh et al., 1998b;[START_REF] Ivanov | The Sulphur Cycle in Continental Reservoirs[END_REF].

These deposits might bias the riverine δ 34 S to high values, since evaporites have δ 34 S values reflecting the seawater δ 34 S during the time of deposition, and range from between 10 to 30‰ (e.g. [START_REF] Kampschulte | The sulfur isotopic evolution of Phanerozoic seawater based on the analysis of structurally substituted sulfate in carbonates[END_REF].

Geochemical calculations tend to form the basis of the most commonly cited sulfur isotope compositions for modern riverine input estimates. Isotope mass balance models of the sulfur cycle have typically employed a riverine δ 34 S value of around 7-8‰ (e.g. [START_REF] Garrels | Coupling of the sedimentary sulfur and carbon cycles; an improved model[END_REF][START_REF] Halevy | Sulfate Burial Constraints on the Phanerozoic Sulfur Cycle[END_REF][START_REF] Kump | Modeling Atmospheric O-2 in the Global Sedimentary Redox Cycle[END_REF][START_REF] Kurtz | Early Cenozoic decoupling of the global carbon and sulfur cycles[END_REF]. These values can be traced back to assumptions about the relative contributions of sulfide and sulfate weathering to the riverine sulfate budget. Specifically, it was assumed that the abundance of sedimentary sulfate minerals is equal to the abundance of sedimentary sulfide minerals and that gypsum weathers twice as fast as pyrite [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF]. These two assumptions imply that sulfate mineral weathering should contribute twice as much sulfate to rivers as pyrite weathering. Thus, if a δ 34 S value of 17‰ is assumed for sulfate in evaporite minerals and a δ 34 S value of -12‰ is assumed for pyrite, then a simple river isotope mass balance predicts an average riverine δ 34 S of between 7 to 8‰ (ignoring anthropogenic and other minor sources of sulfate to rivers). It is important to note that because these calculations assumed a fixed ratio of riverine sulfur from sulfide weathering to sulfate weathering, this isotopic composition cannot then be used to calculate the relative proportion of sulfide weathering. Nonetheless, this approach and its origins are important and informative to the field.

Given the large uncertainties in these estimates of the relative fluxes in the modern biogeochemical sulfur cycle, and the resulting implications for weathering and the modern carbon cycle, the aims of this paper are to: (1) reevaluate the global modern sulfur isotopic composition of riverine sulfate, and

(2) estimate the modern flux of pyrite-derived sulfate supplied to the ocean from rivers using two different and complementary methods: a weathering endmember decomposition, and a simple sulfur isotope mass balance.

Methods

Measurement of river water sulfate and δ 34 S

River waters were sampled either opportunistically or as part of a number of field campaigns between years 1993 -2013. Details of all rivers measured in this study for sulfur isotopes can be found in Supplementary Table 1. Previously published sulfur isotope data from main stems of rivers were compiled from the literature, and can be found in Supplementary Table 2, along with the new main stem data from this study. The locations of the main stems of rivers included in this study can be seen in Figure 1.

The concentration of sulfate in river waters was determined by ion chromatography with a Dionex ICS-2000, using an AS-19 column and 20 mM KOH eluent at the Environmental Analysis Center at Caltech. River samples were then dried down and re-dissolved in 0.01 M HCl. The sulfate was purified from its matrix with an anion exchange column as described in [START_REF] Paris | Neoarchean carbonate-associated sulfate records positive 33S anomalies[END_REF].

New measurements of sulfur isotopes in rivers were made by MC-ICP-MS on a Neptune Plus at Caltech [START_REF] Paris | MC-ICP-MS measurement of δ 34 S and∆ 33 S in small amounts of dissolved sulfate[END_REF]. Measurement by MC-ICP-MS reduces sample size requirements by three orders of magnitude over traditional gas source mass spectrometric methods, and thus only 20 nmol of sulfate were needed for each sample. Typical rivers have micromolar concentrations of sulfate, thus sample sizes were in the range of 100 uL to a few mL of river water depending on concentration. An in-house sodium sulfate solution was used as a bracketing standard on the MC-ICP-MS to correct for instrumental mass bias.

Consistency in chemistry and measurement was monitored with multiple full replicates of a seawater standard (21.04 ± 0.17‰ V-CDT 2 s.d.; n = 20) and an in-house consistency standard from a filtered river water sample collected from the headwaters of the Arroyo Seco in Angeles National Forest, California near Switzer Falls (4.11 ± 0.24‰ V-CDT 2 s.d.; n = 10).

Complete chemistry blanks were monitored along with every set of 10 samples, and contained an average of 0.1 nmol of sulfate. As the smallest samples measured had 20 nmol of sulfate, the blank contamination contributes at most 0.5% of the total sulfate measured, and typically contributes closer to 0.1%, as most samples were analyzed with at least 100 nmol of sulfate. The δ 34 S of the blank is typically close to zero, with a long-term average for all blanks of 2.4 ± 6.5‰ (2 s.d.; n = 16). The measured δ 34 S of all river samples have been blank corrected with the average of all blanks, and the uncertainty on both the amount of blank and its isotopic ratio is propagated through to the final uncertainties on the riverine δ 34 S values.

Global average δ 34 S values were calculated based on the data in Supplementary Table 2. Sulfate fluxes were calculated for these main stem rivers by multiplying sulfate concentration by discharge, and these fluxes were used to calculate the global flux-weighted average and standard deviation of δ 34 S.

End-member weathering calculation

Previous studies have used coupled δ 34 S and δ 18 Ο in sulfate to estimate OWP fluxes in individual river catchments [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF]e.g. Karim and Veizer, 2000;[START_REF] Turchyn | Isotope evidence for secondary sulfide precipitation along the Marsyandi River, Nepal, Himalayas[END_REF]. These estimates do not rely on knowing the δ 34 S of the pyrite being weathered which can have a large range within an individual basin (e.g. [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF]. This is a very powerful technique, but measurement by gas source mass spectrometry requires larger sample sizes, and has only been used on a few rivers, in contrast to the numerous δ 34 S measurements on riverine sulfate that have been previously published. Thus, to maximize the number of rivers included in our global database, we used riverine δ 34 S and major anion and cation concentrations in an end-member weathering calculation following [START_REF] Gaillardet | Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[END_REF] to further constrain the global budget of sulfide weathering.

Rivers that were either measured or compiled for δ 34 S that also had Cl, Na, Mg, Ca, and Sr concentration data (42 rivers total, representing 41% of the global riverine discharge) were included in this end-member deconvolution. A correction for atmospheric deposition from rainwater was made following [START_REF] Gaillardet | Chemical and physical denudation in the Amazon River basin[END_REF] based on the Cl concentration of the river. Each river was assigned a critical value of Cl (ranging from 20 to 100 µM, Supplementary Table 2) based on the evapotranspiration factor of the river basin or nearby rivers calculated from the GEMS-GLORI database [START_REF] Meybeck | GEMS/Water contribution to the global register of river inputs (GLORI)[END_REF], and typical values of Cl concentration in rainwater (14 µM, [START_REF] Moller | The Na/Cl ratio in rainwater and the seasalt chloride cycle[END_REF]. If the concentration of Cl in the river was less than the critical value of Cl, then all of the Cl in that river was attributed to rainwater and typical rainwater ratios of Na to Cl (1.15) were used to determine the rainwater-derived fraction of Na for that river. The remaining Na in the river was assumed to come from the weathering of silicates, carbonates, and evaporites. If the concentration of Cl in the river was greater than the critical value, it was assumed that the contribution of Cl from rainwater was equal to the critical value for that river, with the remaining Cl sourced from evaporites. Following the approach of [START_REF] Gaillardet | Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[END_REF], we then used a series of linear equations to solve for the proportions of sodium in the river that are attributed to evaporite, carbonate, and silicate weathering. The linear equations for each river were solved 10,000 times using a random sampling of weathering end-member values from a normal distribution of the values and associated uncertainties listed in Table 1. The median and standard deviation of these Monte Carlo simulations were used as the end-member fraction and uncertainty for each river (see Supplementary Information).

We calculated the proportion of sulfate that could be attributed to each end-member using the S/Na in the end-members (Table 1), and the S/Na measured in the rivers following equation 1.

𝑆 ! = ! !" ! ! ! ! !" !"# Equation 1
where (S/Na) are molar ratios of sulfate to sodium, the subscript i refers to the silicate, evaporite, carbonate, and rainwater end-members (Table 1), the subscript riv refers to the riverine molar ratio, α is the fraction of sodium attributed to each of the end-members (i) from the end-member weathering calculation, and Si is the fraction of sulfate in the river attributed to each of the end-members (i). Any sulfate that cannot be attributed to one of the four endmembers, we call "excess sulfate". We used Monte Carlo simulations using the end-member values and uncertainties from Table 1 as described above. Results of this calculation are sensitive to the end-member ratios selected, especially for evaporites, which can vary widely. We use the Ca/Na, Mg/Na, and Sr/Na end-member ratios for silicates and carbonates from [START_REF] Gaillardet | Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[END_REF] and [START_REF] Negrel | Erosion Sources Determined by Inversion of Major and Trace-Element Ratios and Strontium Isotopic-Ratios in River Water -the Congo Basin Case[END_REF], and rainwater ratios from [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF] and [START_REF] Gaillardet | Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[END_REF]. For evaporites, we calculate molar ratios based on mineral abundances in evaporites given in [START_REF] Lerman | CO2 and H2SO4 consumption in weathering and material transport to the ocean, and their role in the global carbon balance[END_REF], assuming congruent weathering of evaporites. These evaporite based ratios are consistent with the median ratio calculated from rivers that drain evaporites [START_REF] Meybeck | Concentrations des eaux fluviales en éléments majeurs et apports en solution aux océans[END_REF][START_REF] Meybeck | Concentrations des eaux fluviales en éléments majeurs et apports en solution aux océans[END_REF], giving additional support for our choice of endmembers 1 . The S/Na value for carbonates is based on the S/Ca molar ratio of 0.001 in marine carbonates [START_REF] Volkov | The Sulfur Cycle in Oceans[END_REF]) and the Ca/Na value of carbonates of 60 [START_REF] Negrel | Erosion Sources Determined by Inversion of Major and Trace-Element Ratios and Strontium Isotopic-Ratios in River Water -the Congo Basin Case[END_REF]. Finally, the S/Na value of the silicate endmember is calculated from global estimates of the abundance of S and Na2O in continental igneous and metamorphic rocks [START_REF] Holser | Isotope geochemistry of sedimentary sulfates[END_REF][START_REF] Rudnick | Composition of the continental crust[END_REF]. This estimate of S/Na in silicate rocks deliberately does not include sedimentary sulfides (e.g. from shales). Consequently, any sulfate derived from oxidative weathering of sedimentary pyrite would be incorporated into the excess sulfate parameter described above.

Results

The river waters analyzed in this study yield a large range of δ 34 S values, between -13.4‰ and 21.7‰ (Figures 1 and2), and represent more than 46% of the global water discharge to the ocean (3.88 x 10 4 km 3 /y (Peucker-Ehrenbrink, 2009)). This dataset includes smaller tributaries, as well as time series measurements from individual major rivers (Fraser, Lena, Ob', Yenisei, Yukon, Mackenzie, Kolyma, and Congo rivers). Within an individual river basin, there can be a large range in the δ34 S of tributaries. For instance, δ 34 S values from the Congo River and its tributaries range from 3.7 to 12.5‰ (n = 7), those from the Ganges River system range from -1 to 13.3‰ (n = 3), and those in the Amazon River basin range from 4.5 to 13.4‰ (n= 19), though most of the Amazonian tributaries have a narrow range of 4.5 to 7‰ (n=16), consistent with the findings of [START_REF] Longinelli | Isotope geochemistry of the Amazon basin: a reconnaissance[END_REF].

The δ 34 S measured in river time series samples can also show a large range, highlighting temporal heterogeneity in the isotopic composition of riverine sulfate, that is likely driven by varying contributions from tributaries that drain different lithologies (Figure 3). For instance the δ 34 S of the Lena River has a range of 7‰ between May and November, the δ 34 S of the Mackenzie River has a range of 4.5‰ between March and September, and the δ 34 S of the Fraser River has a range of 5‰ over the course of a year (Supplementary Table 1;

Figure 3). In contrast, time series data from the Congo, Kolyma, Ob', Yenisei, and Yukon rivers have more constant (<2‰ variation) δ 34 S over the course of several months to a year (Supplementary Table 1).

The flux-weighted average of the rivers gives a δ 34 S value of 4.4 ± 4.5‰ (1 s.d.). If the polluted rivers (([Na] -[Cl])/[Cl] < 0.1) are excluded 2 from this estimate, the average δ 34 S of riverine sulfate is 4.8 ± 4.9‰, and accounts for more than 43% of the water discharge to the ocean. The similarity between these two values suggests that the most polluted rivers do not meaningfully change the estimate of the global riverine δ 34 S value. We therefore take the value of 4.8‰

as our best estimate of the modern pre-anthropogenic δ 34 S value of river water.

The flux-weighted distribution of the isotopic composition of riverine sulfate (Figure 4) shows that half of the rivers have δ 34 S values between 1.8 and 8.2‰ We can assess how representative our subset of rivers is of global riverine discharge by comparing the flux-weighted Mg/Na and Ca/Na ratios and strontium isotopes from our dataset with global estimates of these values from previously published studies. Our flux weighted Mg/Na ratio of 0.66 and Ca/Na ratio of 1.63 are within the range of previous global estimates of these values of 0.57 to 0.70 and 1.32 to 1.67, respectively [START_REF] Meybeck | GEMS/Water contribution to the global register of river inputs (GLORI)[END_REF][START_REF] Miller | Reassessing the surface cycling of molybdenum and rhenium[END_REF]. The same is true if we scale the river fluxes by the large scale drainage regions [START_REF] Graham | Five-minute, 1/2°, and 1° data sets of continental watersheds and river networks for use in regional and global hydrologic and climate system modeling studies[END_REF][START_REF] Peucker-Ehrenbrink | Quantitative bedrock geology of the continents and large-scale drainage regions[END_REF] following [START_REF] Miller | Reassessing the surface cycling of molybdenum and rhenium[END_REF]. In this case, our subset of rivers yields a Mg/Na ratio of 0.68 and a Ca/Na ratio of 1.70, compared to previously published values of 0.55 to 0.71 and 1.27 to 1.74, respectively [START_REF] Meybeck | GEMS/Water contribution to the global register of river inputs (GLORI)[END_REF][START_REF] Miller | Reassessing the surface cycling of molybdenum and rhenium[END_REF].

Thus, our elemental ratios support the use of this subset of rivers for global budget estimates. The flux weighted average 87 Sr/ 86 Sr ratio from our subset of rivers with published ratios (40% of global water discharge) is 0.7115, slightly higher than estimates of global riverine 87 Sr/ 86 Sr (0.7111 to 0.7114 [START_REF] Peucker-Ehrenbrink | Continental bedrock and riverine fluxes of strontium and neodymium isotopes to the oceans[END_REF][START_REF] Vance | Variable Quaternary chemical weathering fluxes and imbalances in marine geochemical budgets[END_REF]). The slightly more radiogenic volcanically dominated river basins (e.g. smaller rivers on volcanic islands and along active margins) in our subset of rivers. This might bias our estimate of global riverine δ 34 S to higher values, but the effect is likely small given that rivers draining volcanic islands (e.g. Reunion, Iceland) have relatively low concentrations of sulfate (34-70 µM) and have δ 34 S values only slightly (1-3‰) lower than our global average (Tables 1 and2).

Weathering end-member decomposition

End-member decomposition allows us to tease apart the different sources of sulfate to rivers. A plot of δ 34 S of riverine sulfate versus the percentage of Na derived from silicate weathering illustrates the importance of pyrite weathering in influencing the sulfur isotopic composition of rivers (Figure 5). For rivers with no Na derived from silicate weathering, the only other two sources of Na are evaporite and carbonate minerals. As both of these lithologies reflect the sulfur isotopic composition of seawater when they were formed, they should both have high δ 34 S values (~17‰, [START_REF] Holser | Isotope geochemistry of sedimentary sulfates[END_REF]). Similarly, for rivers that derive 100% of their Na from silicate weathering (none from evaporites or carbonates), the δ 34 S of river water should be similar to the silicate value of 4‰ [START_REF] Holser | Isotope geochemistry of sedimentary sulfates[END_REF]. These two end-members determine the red mixing lines in Figure 5. Most data plot below these mixing lines, implying a 34 S-depleted source of sulfate to the rivers that does not also add Na. The oxidative weathering of sedimentary pyrite could reasonably constitute such a source.

We further quantify the sources of sulfate to rivers by using the endmember values of S/Na from Table 1 and Equation 1to convert fractional contributions of Na from the different end-members into fractional contributions of S from the different end-members. Based on our end-member decomposition, evaporite weathering contributes 31 ± 6% (1.5 ± 0.2 Tmol/y) to the global riverine sulfate budget. We find that the contribution of rainwater sulfate to rivers in our database is minor, contributing about 1.6% (0.07 Tmol/y) to the total riverine sulfate flux, consistent with independent estimates of the contribution of cyclic sea salts (0.09 -0.13 Tmol/y; [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF][START_REF] Eriksson | The yearly circulation of sulfur in nature[END_REF]) to riverine sulfate. The contribution of carbonate and silicate weathering is also minor with only 0.01 and 0.09 Tmol/y derived from each of those end-members, respectively.

The remaining sulfate in each river that cannot be readily attributed to carbonate, silicate, and evaporite weathering, or to rainwater (sea-salt) constitutes a parameter we call "excess sulfate". This "excess sulfate" is a combination of the other major sources of sulfate to rivers: pollution, volcanic atmospheric deposition, biogenic emissions, and the oxidative weathering of sedimentary pyrite. This relative contribution of excess sulfate ranges from 0 to more than 97% for all rivers, with the exception of the Narmada and Tapti rivers, which have negative excess sulfate values. A negative excess sulfate value implies that the sulfate in those rivers can be more than accounted for by weathering of the three end-members and rainwater inputs, and thus could indicate an unidentified sink for sulfate along the river (e.g. in the river flood plain). For example, from sulfur isotope data it was suggested that removal of riverine sulfate might occur via sulfate reduction and sulfide precipitation in weakly developed soils in the Marsyandi River catchment in the Himalaya [START_REF] Turchyn | Isotope evidence for secondary sulfide precipitation along the Marsyandi River, Nepal, Himalayas[END_REF]. However, the extent to which this process occurs in other river systems remains to be established [START_REF] Torres | The acid and alkalinity budgets of weathering in the Andes-Amazon system: Insights into the erosional control of global biogeochemical cycles[END_REF]. The excess sulfate values thus represent minimum values, and for the remainder of this study we will ignore the two rivers in our compilation that have a negative excess sulfate.

The flux weighted global percentage of excess sulfate is 65 ± 6% (3.1 ± 0.2 Tmol/y), which is dominated by pyrite weathering and pollution. The contributions to riverine sulfate from volcanoes and biogenic emissions are both relatively small (0.34 and 0.14 Tmol/y, respectively [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF]).

Using an estimate of the pollutive sulfate in rivers of 1.3 Tmol/y [START_REF] Meybeck | Global Occurrence of Major Elements in Rivers[END_REF] ) and the Kolyma (-4.2 to -5.6‰) rivers, which also drain catchments rich in sedimentary rocks, including black shales [START_REF] Edmond | Fluvial geochemistry of the eastern slope of the northeastern Andes and its foredeep in the drainage of the Orinoco in Colombia and Venezuela[END_REF]Huh et al., 1998a). These examples are consistent with the notion that local lithology plays a key role in setting the sulfur isotopic composition of rivers, and the important role of black shale weathering in the sulfur budget.

Plots of the δ 34 S versus inverse concentration of sulfate and versus the riverine Na/SO4 ratio (Figure 2) show a triangular pattern, indicative of at least three end-members contributing to sulfate in rivers. The most sulfate-rich rivers span a wide range of δ 34 S values (-13 to 21‰), whereas the rivers with the lowest sulfate concentrations cluster tightly between ~5-10 ‰. Rivers plotting in the high concentration-high δ 34 S region, such as the Lena River, are likely dominated by sulfate minerals weathered from evaporites. Gypsum (CaSO4) dissolves easily, and rivers draining evaporitic rocks usually have high concentrations of total dissolved solids, including sulfate. These rivers have Na/S values ranging between 1.7 and 5, consistent with an evaporite endmember.

Rivers plotting in the high concentration-low δ 34 S region could have sulfate sources dominated by the oxidative weathering of pyrite, which has low δ 34 S values, or by pollution, which is often cited in the literature as having low or negative δ 34 S signatures (e.g. [START_REF] Ivanov | The Sulphur Cycle in Continental Reservoirs[END_REF]. We screened for polluted rivers using ancillary anion and cation data [START_REF] Fiege | Chemical Geology[END_REF][START_REF] Miller | Reassessing the surface cycling of molybdenum and rhenium[END_REF][START_REF] Peterson | Arctic Great Rivers Observatory 1 Biogeochemistry Data[END_REF][START_REF] Voss | Tracing river chemistry in space and time: Dissolved inorganic constituents of the Fraser River, Canada[END_REF] River, which is influenced by Cl-rich brines [START_REF] Gaillardet | Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[END_REF]. As shown in Figure 2, the polluted rivers tend to have high sulfate concentrations (>100 µM, with the exception of the Doubs whose main sulfate sources are atmospheric [START_REF] Calmels | Sensitivity of carbonate weathering to soil CO2 production by biological activity along a temperate climate transect[END_REF]). However, there is a large range of δ 34 S for these rivers; thus they do not simply populate the bottom left of the graph, as would be expected if pollution always had a low δ 34 S signature. Indeed, measurements of δ 34 S on oil, gas, coal and sulphide ores range from -25 to +30‰ [START_REF] Newman | Sulphur Isotope Variations in the Atmosphere[END_REF], and δ 34 S of fertilizers, another anthropogenic sulfate source, can range from -6.5 to 20.7‰ [START_REF] Szynkiewicz | Anthropogenic sulfate loads in the Rio Grande, New Mexico (USA)[END_REF][START_REF] Vitoria | Fertilizer Characterization: Isotopic Data (N, S, O, C, and Sr)[END_REF]. Since pollution cannot be the sole reason for the high sulfate concentration and low δ 34 S, it follows that the oxidation of pyrite is a major source of sulfate to these rivers. This is consistent with the abundance of shales and pyrites in the catchment areas of the rivers (e.g. the Kolyma (Huh et al., 1998a), Orinoco [START_REF] Edmond | Fluvial geochemistry of the eastern slope of the northeastern Andes and its foredeep in the drainage of the Orinoco in Colombia and Venezuela[END_REF], Buller [START_REF] Robinson | Discrimination of sulfur sources in pristine and polluted New Zealand river catchments using stable isotopes[END_REF], Kaoping [START_REF] Das | Disproportionately high rates of sulfide oxidation from mountainous river basins of Taiwan orogeny: Sulfur isotope evidence[END_REF], Skeena [START_REF] Spence | The role of sulfur in chemical weathering and atmospheric CO2 fluxes: Evidence from major ions, δ13CDIC, and δ34SSO4 in rivers of the Canadian Cordillera[END_REF], and Stikine rivers [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF]).

Rivers and tributaries with low concentrations of sulfate (e.g. Amazon,

Congo, and Niger rivers) tend to have δ 34 S values of 5 -10‰ (Figure 2). The low sulfate concentrations likely reflect the lack of any significant weathering of evaporites or pyrite-rich shales in these drainage basins. Their low concentration of sulfate could also be attributed in part to the large flow of water and the dilution of solutes, as the concentration of total dissolved solids for these rivers is low at 35-59 mg/L [START_REF] Meybeck | GEMS/Water contribution to the global register of river inputs (GLORI)[END_REF]. However, the small range of δ 34 S values indicates a similar source of sulfate to the rivers, leaving us to conclude that silicate weathering of igneous rock and atmospheric deposition of sulfate (sea salts) are the major sources of sulfate to these rivers and tributaries, consistent with previous studies [START_REF] Negrel | Erosion Sources Determined by Inversion of Major and Trace-Element Ratios and Strontium Isotopic-Ratios in River Water -the Congo Basin Case[END_REF][START_REF] Gaillardet | Chemical and physical denudation in the Amazon River basin[END_REF].

Rivers with a high ratio of Na/S also plot within this narrower range of δ 34 S, consistent with a silicate or atmospheric source.

The wide range in δ 34 S from tributaries from an individual river basin (e.g. Congo , Ganges, Mackenzie) highlights the spatial heterogeneity of sulfur sources to rivers, that is most likely controlled by heterogeneity of the bedrock geology-and commensurate differences in the isotopic composition of the sulfur-bearing phases-in their catchments. The differences in the temporal variability of δ 34 S between individual river basins indicate that some rivers (e.g.

Lena, Mackenzie, and Fraser) have weathering conditions that are more heterogeneous than other river basins with more constant δ 34 S, and that the varied lithologies in these river catchments contribute differently to the dissolved load at distinct times throughout the year. For instance, in the case of the Fraser River, the δ 34 S of sulfate over the course of a year shows a strong correlation (r 2 = 0.88) with the dissolved 87 Sr/ 86 Sr ratio in these waters [START_REF] Voss | Tracing river chemistry in space and time: Dissolved inorganic constituents of the Fraser River, Canada[END_REF] The importance of pyrite weathering in setting the δ 34 S value of rivers globally can be seen in Figure 6. This figure shows the relationship between the measured δ 34 S versus the δ 34 S values predicted from a sulfur isotope mass balance in each river based on the weathering end-member decomposition. This sulfur isotope mass balance uses the proportion of sulfate derived from the weathering and rain water end-members and assumes end-member δ 34 S values of 17‰ for evaporite and carbonate weathering, 4‰ for silicate weathering, and 20‰ for rainwater [START_REF] Herut | Sources of sulfur in rainwater as indicated by isotopic δ34S data and chemical composition, Israel[END_REF][START_REF] Holser | Isotope geochemistry of sedimentary sulfates[END_REF]. The 1:1 line in Figure 6 indicates rivers whose measured δ 34 S values can be entirely accounted for by sulfur input from sources besides pyrite weathering and pollution. Most rivers plot to the left of this line indicating that the measured values of δ 34 S are lower than predicted. Furthermore, the rivers that plot farthest from the line have large excess sulfate proportions (calculated from the end-member decomposition), implying that the sulfur isotopic composition of this excess sulfate has a low δ 34 S value, consistent with a pyrite source3 .

Global dissolved sulfur isotope mass balance for rivers

The riverine δ 34 S values reported here (4.4±4.5‰ and 4.8 ± 4.9‰

excluding the most polluted rivers) are lower than the previous compilation value from Eurasian rivers of 9‰ [START_REF] Ivanov | The Sulphur Cycle in Continental Reservoirs[END_REF] and the values of 7 to 8‰

typically used in models of the global sulfur cycle [START_REF] Halevy | Sulfate Burial Constraints on the Phanerozoic Sulfur Cycle[END_REF]e.g. Kump and Garrels, 1986;[START_REF] Kurtz | Early Cenozoic decoupling of the global carbon and sulfur cycles[END_REF]. At face value these lower values of the δ 34 S of rivers implies that a greater fraction of riverine sulfate derives from reduced sources (e.g. pyrite) than previously estimated. This finding is in line with recent studies illustrating the high relative contribution of oxidative weathering of pyrite to riverine sulfate based on detailed studies of individual river basins (e.g.

Mackenzie [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF], Kaoping [START_REF] Das | Disproportionately high rates of sulfide oxidation from mountainous river basins of Taiwan orogeny: Sulfur isotope evidence[END_REF] , Marysandi [START_REF] Turchyn | Isotope evidence for secondary sulfide precipitation along the Marsyandi River, Nepal, Himalayas[END_REF], Amazon [START_REF] Torres | The acid and alkalinity budgets of weathering in the Andes-Amazon system: Insights into the erosional control of global biogeochemical cycles[END_REF] and Ganges-Brahmaputra [START_REF] Galy | Weathering processes in the Ganges-Brahmaputra basin and the riverine alkalinity budget[END_REF])), and we can now expand those conclusions to the global budget.

Using a simple terrestrial weathering isotope mass balance akin to that used in [START_REF] Kump | Modeling Atmospheric O-2 in the Global Sedimentary Redox Cycle[END_REF], we can calculate the fraction of pyrite weathering contributing to riverine sulfate:

fowp = -(δevap -δriv)/(δpy -δevap) Equation 2
where fowp is the fraction of pyrite-derived sulfate, δevap is the δ 34 S of evaporite (17‰ [START_REF] Holser | Isotope geochemistry of sedimentary sulfates[END_REF])), δriv is the δ 34 S of riverine sulfate (4.8‰, above), and δpy is the δ 34 S of pyrite from sedimentary rocks (-12‰ [START_REF] Holser | Isotope geochemistry of sedimentary sulfates[END_REF]). This gives a fowp in rivers of 0.42, substantially higher than the estimate of fowp of 0.28 if the riverine δ 34 S value is taken to be 9‰ [START_REF] Ivanov | Major Fluxes of the Global Biogeochemical Cycle of Sulphur[END_REF].

These estimates are very sensitive to the assumption of the end-member δ 34 S values of the oxidized (sulfate salts) and reduced (sulfide and disulfide) sedimentary reservoirs, which are difficult to constrain from a global perspective. Estimates of the global δ 34 S values of oxidized S in sediments range from 12.1‰ to 19‰, and of reduced sediments from -7.4‰ to -16‰ (Holser et al., 1988 and references therein). These different combinations of δevap and δpy give a range of fowp between 0.35 to 0.42 (compared to fowp = 0.16 to 0.29 if the river δ 34 S value is taken to be 9‰), with a median fowp of 0.41. Despite the large uncertainties in the correct end-member isotope values, the estimate of fowp = 0.41 is broadly consistent with the independently derived fluxes from the weathering end-member decomposition that indicate that pyrite weathering contributes up to 46% of the sulfate flux in rivers due to weathering. The weathering end-member decomposition is based solely on major ion concentrations and not isotope values, and so is independent of the isotope mass balance calculation in Equation 2.

Implications for the modern sulfur cycle and marine isotope mass balance

Our flux estimate of riverine sulfate from evaporite weathering (1.5 Tmol/y) that we derived from weathering end-member decomposition is larger than the assumed value of 1 Tmol/y of sulfate weathering by [START_REF] Berner | Burial of Organic-Carbon and Pyrite Sulfur in Sediments Over Phanerozoic Time -a New Theory[END_REF] and [START_REF] Garrels | Phanerozoic Cycles of Sedimentary Carbon and Sulfur[END_REF], which is widely used in box models of the sulfur cycle and is smaller than the value estimated more recently by [START_REF] Lerman | CO2 and H2SO4 consumption in weathering and material transport to the ocean, and their role in the global carbon balance[END_REF] (1.98 Tmol/y). Our estimate of the sulfide weathering flux from the weathering end-member decomposition (1.3 Tmol/y) is about twice as large as the commonly cited estimates of the global OWP flux (0.5-0.64

Tmol/y [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF][START_REF] Lerman | CO2 and H2SO4 consumption in weathering and material transport to the ocean, and their role in the global carbon balance[END_REF]), and is consistent with our estimate from an isotope mass balance discussed in the previous section that is based on the revised global average riverine δ 34 S value of 4.8‰. Thus the isotope and river chemistry data independently imply a larger contribution of oxidative weathering of pyrite to the riverine sulfate flux than previously estimated 4 .

It can be informative to investigate the derivation of previous estimates of the sulfate flux from OWP in order to understand why they were significantly lower than our new value. Previous estimates of global OWP from [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF] and [START_REF] Lerman | CO2 and H2SO4 consumption in weathering and material transport to the ocean, and their role in the global carbon balance[END_REF] both relied on assumptions of the relative abundance of pyrite in sedimentary rocks, though they do this differently. The 0.5 Tmol/y estimate [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF] assumed that the riverine ratio of calcium from weathering of sedimentary rocks to pyrite-derived sulfate must be the same as the ratio of calcium to pyrite sulfur in common 4 In order to close the riverine sulfate budget, previous estimates attributed any remaining sulfate that did not derive from pyrite or evaporite weathering to pollution (e.g. [START_REF] Berner | Global environment: water, air, and geochemical cycles[END_REF], rather than using an independent estimate of sulfate pollution [START_REF] Meybeck | Global Occurrence of Major Elements in Rivers[END_REF] as we have done here.

sedimentary rocks (which was taken to be 8.5 and 0.3 weight % of calcium and pyrite-S, giving a molar pyrite-S/Ca ratio of 0.04 [START_REF] Garrels | Coupling of the sedimentary sulfur and carbon cycles; an improved model[END_REF]).

The molar ratios of excess sulfate/Ca from our end-member decomposition are more typically around 0.1 and can be as high as 0.6 for catchments with abundant shales (e.g. Kolyma (Huh et al., 1998a), Orinoco [START_REF] Edmond | Fluvial geochemistry of the eastern slope of the northeastern Andes and its foredeep in the drainage of the Orinoco in Colombia and Venezuela[END_REF],

and Kaoping [START_REF] Das | Disproportionately high rates of sulfide oxidation from mountainous river basins of Taiwan orogeny: Sulfur isotope evidence[END_REF]). Thus the assumed ratio of calcium to pyrite sulfur in common sedimentary rocks might not fully account for the presence of pyrite-rich shales. The OWP flux estimate of 0.64 Tmol/y [START_REF] Lerman | CO2 and H2SO4 consumption in weathering and material transport to the ocean, and their role in the global carbon balance[END_REF] assumes that 24% of sedimentary sulfur is in the reduced form of pyrite. This percentage is significantly lower than previous estimates, which are closer to 50% (e.g. Holser et al., 1988 and references therein). Thus both of these prior estimates potentially underestimated the global abundance of pyrite S in sedimentary rocks by a factor of two, which could help to explain the underestimation of global OWP fluxes.

Our revised estimate of the OWP flux is similar to the estimate of global pyrite burial in ocean sediments calculated from the estimates of global organic carbon burial and the C/S value in modern sediments (1.22 Tmol/y; [START_REF] Berner | Burial of organic carbon and pyrite sulfur in the modern ocean; its geochemical and environmental significance[END_REF]). It is also similar to the upper range of estimates of pyrite burial (36% of the total natural, pre-anthropogenic flux) based on marine isotope mass balance and constraints from marine δ 33 S [START_REF] Tostevin | Multiple sulfur isotope constraints on the modern sulfur cycle[END_REF]. Our upward revision of global pyrite weathering brings the oceanic input and output of reduced sulfur, which previously differed by a factor of 2, into close agreement.

The uncertainties on these modern fluxes from the weathering endmember decomposition of river concentration data alone are still too large (0.2 Tmol/y) to assess whether the modern sulfur cycle is in steady state with regard to the relative amounts of oxidized and reduced sulfur that are weathered and buried. However, if we assume steady-state, our new riverine δ 34 S value can be used in a marine isotope mass balance that solves for the fraction of the sulfate removed from the ocean through pyrite burial (fpyr): fpyr = (δrivδsw)/(Δ 34 Spyrite-seawater) Equation 3where δsw is the δ 34 S of seawater (21‰; [START_REF] Paris | MC-ICP-MS measurement of δ 34 S and∆ 33 S in small amounts of dissolved sulfate[END_REF]) and Δ 34 Spyrite-seawater is the isotope fractionation between seawater sulfate and sedimentary pyrite in the modern ocean (δpyδsw). A plot of the solutions to Equation 3 for our new value for δriv (4.8‰) and the previous value for δriv (9‰) is shown in Figure 7 for a range of Δ 34 Spyrite-seawater. Our lower value of δriv = 4.8‰ implies that there is a greater fraction of pyrite burial than previously predicted based on a δriv = 9‰.

Using this model to determine a specific value for the modern fraction of pyrite burial is not currently possible because the global average value of Δ 34 Spyriteseawater is difficult to constrain given the large range of δ 34 S values measured in modern pyrites.

If, however, we assume that the fraction of pyrite burial (fpyr) is equivalent to the fraction of pyrite weathered (fowp) as calculated either by the weathering end-member decomposition from major ion chemistry (0.46) or riverine isotope mass balance (0.41), then our data suggest that Δ 34 Spyrite-seawater should be -35‰

or -39‰, respectively. These estimates are roughly consistent with canonical assumptions of, Δ 34 Spyrite-seawater = -35‰ (e.g. [START_REF] Garrels | Phanerozoic Cycles of Sedimentary Carbon and Sulfur[END_REF], but are much more positive than flux weighted estimates of this value based on pyrite measurements from modern continental shelf, slope and rise sediments (Δ 34 Spyrite-seawater = -48‰ [START_REF] Yu Lein | The Sulphur Cycle in Oceans, in: The Global Biogeochemical Sulphur Cycle, The Mass-Isotopic Balance of Sulphur in Oceanic Sediments[END_REF]). Modelling studies that use constraints from δ 33 S also predict a greater difference between the δ 34 S values of seawater and bulk pyrite buried (Δ 34 Spyrite-seawater = -56 ± 5‰) [START_REF] Tostevin | Multiple sulfur isotope constraints on the modern sulfur cycle[END_REF]. These larger offsets in the isotopic composition of seawater and buried pyrite would imply that the fraction of pyrite buried is 0.29 and 0.34 for Δ 34 Spyriteseawater = -56 and Δ 34 Spyrite-seawater = -48, respectively (Figure 7), significantly less than our two independent estimates of the fraction of pyrite weathered. If correct, these estimates indicate that the biogeochemical sulfur cycle may not be balanced, and that more pyrite is being weathered than is being buried in modern marine sediments. This conclusion would be consistent with recent ice core evidence that shows a decrease in atmospheric O2 over the past 800,000 years [START_REF] Stolper | A Pleistocene ice core record of atmospheric O2 concentrations[END_REF]. The magnitude of the imbalance between weathering and burial fluxes of reduced carbon and sulfur that is required to explain the long term oxygen trend in ice cores is ~2% [START_REF] Stolper | A Pleistocene ice core record of atmospheric O2 concentrations[END_REF]; this number is well within uncertainties of the weathering flux estimates (~0.2 Tmol/y). However, given the large range of estimates of Δ 34 Spyrite-seawater more work constraining global Δ 34 Spyrite-seawater from modern sediments would help to resolve this issue.

Our increased estimate of global OWP means that the previously published OWP rates from three river basins (Mackenzie [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF], Kaoping [START_REF] Das | Disproportionately high rates of sulfide oxidation from mountainous river basins of Taiwan orogeny: Sulfur isotope evidence[END_REF], and Jialing [START_REF] Li | Identification of dissolved sulfate sources and the role of sulfuric acid in carbonate weathering using dualisotopic data from the Jialing River, Southwest China[END_REF]) only account for ~12%

(rather than a third) of the total global OWP flux estimate. Nonetheless, the fact that these rivers drain less than 2% of global surface area highlights the heterogeneity of pyrite weathering; simply put, some river basins will contribute disproportionately to the global pyrite weathering budget, based on their tectonic and geological characteristics. The spatial heterogeneity of pyrite weathering thus has important implications for the long-term variations in the δ 34 S of riverine input: changes in paleogeography, paleohydrology, eustasy, and the formation and uplift of sedimentary basins will change both the amounts and types of rock exposed to weathering, which in turn will likely change the value of riverine δ 34 S input to the oceans through time.

Finally, the increased estimate of the OWP flux has implications for the modern carbon cycle and carbon-weathering feedbacks. As pyrite oxidation coupled to carbonate weathering can provide a source of CO2 to the atmosphere (e.g. [START_REF] Calmels | Sustained sulfide oxidation by physical erosion processes in the Mackenzie River basin: Climatic perspectives[END_REF][START_REF] Das | Disproportionately high rates of sulfide oxidation from mountainous river basins of Taiwan orogeny: Sulfur isotope evidence[END_REF][START_REF] Lerman | CO2 and H2SO4 consumption in weathering and material transport to the ocean, and their role in the global carbon balance[END_REF][START_REF] Torres | The acid and alkalinity budgets of weathering in the Andes-Amazon system: Insights into the erosional control of global biogeochemical cycles[END_REF]2015), the greater flux of pyrite-derived sulfate indicates a more significant role for this process in the modern biogeochemical sulfur and carbon cycles. Modern weathering and carbon budgets should be reassessed on a global scale to account for an increased global contribution from OWP.

Conclusion

We determined the modern δ 34 S value of riverine sulfate to be 4. Jun606" Aug606" Oct606" Dec606" Feb607" Apr607" Jun607" Aug607" Oct607" dissolved" 87 Sr/ 86 Sr"

Jun""""Aug""""Oct""""Dec"""Feb""""Apr""""Jun""""Aug"""""Oct""" isotope mass balance, using the S/Na end-member ratios from Table 1, and assuming that the δ 34 S value of silicate rocks is 4‰, that the δ 34 S value of evaporites is 17‰, and that there is no contribution of Na from carbonates. The red dashed line is the predicted δ 34 S based on a riverine isotope mass balance, using the S/Na end-member ratios from Table 1, and assuming that the δ 34 S value of silicate rocks is 4‰, that the δ 34 S value of evaporites and carbonates is 17‰, and that carbonates contribute 12% of the Na (maximum contribution found for all rivers in this study), with the remainder coming from silicates and evaporites. oxidized sulfur for a given riverine sulfate δ 34 S value, and assuming steady state isotope mass balance in the ocean. The red line with stars represents our estimate of the modern riverine δ 34 S value from this study (4.8‰). Previous compilations of the riverine δ 34 S value were higher (~9‰, black circles; [START_REF] Ivanov | The Sulphur Cycle in Continental Reservoirs[END_REF]). Estimates of the relative contribution of weathering of pyrite (fowp) to modern riverine sulfate from our two independent approaches are shown for reference (fowp = 0.46 (solid line, based on our weathering end-member decomposition from major ion chemistry of the rivers) and fowp = 0.41 (dashed line, based on riverine isotope mass balance)).

(

  interquartile range), including the rivers with the largest flux of sulfate (Amazon and Yangtze rivers). Notable exceptions to this include the Mississippi and Lena rivers, which contribute a large flux of sulfate to the oceans, but have low (-5‰) and high δ 34 S (19‰) values, respectively.

  . Dissolved 87 Sr/ 86 Sr values decrease along the flow path of the Fraser River reflecting the underlying geology, as the river flows from Precambrian/Paleozoic sedimentary rocks of the Rocky Mountains, through the young, magmatic and metamorphic rocks of the Coast Range (Cameron and Hattori, 1997). The time series data indicate a greater relative contribution of the Coast Range to the dissolved ion budget during the winter. The sulfur isotope compositions imply that during the winter months, more sulfate is derived from the oxidative weathering of 34 S-poor sulfide-bearing phases in the Coast Range, whereas during the summer a greater proportion of the sulfate derives from sulfate minerals within Paleozoic strata.

FiguresFigure 1 .Figure 2 .

 12 Figures

Figure 3 .

 3 Figure 3. Time series data from Fraser River samples collected between June

Figure 4 .Figure 5 .

 45 Figure 4. Stacked bar graph showing the δ 34 S values (‰, VCDT) and the flux of

Figure 6 .

 6 Figure 6. Plot of predicted δ 34 S (based on assumptions of end-members from

Figure 7 .

 7 Figure 7. Sensitivity of the proportion of sulfur buried marine basins as pyrite

Table 1 . End-member molar ratios (±2 s.d.)

 1

	Molar Ratio	Rain	Silicate	Carbonate	Evaporite
	Ca/Na	0.023 ± 0.01	0.35 ± 0.25	60 ± 30	0.5 ± 0.5
	Mg/Na	0.11 ± 0.01	0.25 ± 0.2	30 ± 15	0.1 ± 0.08
	Sr/Na	0.0002 ± 0.0001	0.003 ± 0.001	0.04 ± 0.02	0.003 ± 0.002
	Cl/Na	1.15 ± 0.1	-	-	1 ± 0.2
	S/Na	0.06 ± 0.01	0.019 ± 0.005	0.06 ± 0.03	0.4 ± 0.2

  , we estimate that the global contribution of oxidative weathering of pyrite to riverine sulfate is 1.3 Tmol/y, or about 28% of the total sulfate flux from allThe highest δ 34 S values measured (14.4 to 21.7‰) come from the Lena and Yenisei rivers which drain the Siberian Platform-a region of large evaporite and carbonate platforms deposited from late Proterozoic through Paleozoic time(Huh et al., 1998b). Since these rocks are abundant in sulfate phases derived from seawater, they are characterized by relatively high δ 34 S values. The sulfate minerals weather readily, resulting in both high concentrations of sulfate and high δ 34 S values for these rivers. The lowest δ 34 S values (-8.5 to -13.4‰) come from tributaries of the Santa Clara River (e.g. Sespe) in Southern California, which weather sandstones, siltstones, and organic-rich shales of the Monterey Formation and equivalent units. Sedimentary evaporites are absent from the underlying strata, and thus sulfate in these rivers is likely dominated by the oxidation of pyrite and organic sulfur molecules, which, due to the biological isotope fractionations associated with sulfate reduction in organic-rich sediments during their diagenesis and lithification, results in lower overall δ 34 S values. The lowest δ 34 S values measured in major rivers come from the Orinoco(-5.1‰

	4 Discussion
	4.1 Controls on riverine sulfate δ 34 S

sources to the ocean (4.7 Tmol/y), and 46% of the flux of sulfate from weathering (i.e. 2.8 Tmol/y, excluding sulfate from pollution, volcanoes, biogenic emissions, and cyclic sea salts).

  . If the value of ([Na] -[Cl])/[Cl] is less than 0.1, then the river is considered polluted. The main sources of Cl to rivers are sea salt aerosols, evaporites, and pollution. As both sea salt aerosols and evaporites deliver Na and Cl in roughly equal proportions, an excess of Cl relative to Na may be indicative of a polluted river. The following rivers have high Cl relative to Na: Connecticut, Danube, Indus, Doubs, Lena, Mississippi, Neman, Rhine, Seine, and St. Lawrence. These rivers are also located near large human populations or near industrial centers. The one exception to this is the Lena

  Tmol S/y, respectively. This estimate of pyrite weathering is twice as high as previous estimates, and is consistent with isotope mass balance constraints on the δ 34 S of river water that imply that ~40% of non-anthropogenic riverine sulfate derives from weathered pyrite. Since pyrite oxidation can act as a source of CO2 to the atmosphere, this new estimate should reduce the global carbon sink attributed to chemical weathering. Basin lithology exerts a strong control on the δ 34 S values of riverine sulfate, questioning the validity of assuming constant δ 34 S of riverine input to the ocean over geological time. Stratigraphic trends observed in sulfur isotopes in ancient sedimentary rocks might therefore reflect secular changes in the isotopic composition of the riverine inputs, rather than changes in the relative burial of marine pyrite versus sulfate salt evaporites.

	1.3
	4 ±
	4.5‰, based on measurements of rivers that account for more than 46% of
	global freshwater discharge. Removing highly polluted rivers from the estimate
	gives a similar δ 34 S value of riverine sulfate of 4.8 ± 4.9‰. These δ 34 S values are
	lower than previous estimates of the isotopic composition of riverine sulfate, and
	imply a greater contribution of pyrite weathering to riverine sulfate than has
	previously been assumed. An end-member decomposition of weathering
	sources to these rivers indicates that the relative contributions of evaporite
	(sulfate) and pyrite (sulfide) weathering to riverine sulfate budgets are ~1.5 and

The one exception to this statement is evaporite Mg/Na. The Mg/Na ratio based on mineral abundances is 0.03, while the Mg/Na ratio from rivers draining evaporites is 0.1. However, the final results are not sensitive to the choice of this parameter, with only a 0.02 Tmol/y difference in the estimate of the global flux of sulfate from evaporites, well within the uncertainty of the calculation (0.2 Tmol/y).

The only river with sufficient data to estimate a pre-anthropogenic, natural δ

S and flux of sulfate is the Mississippi River.[START_REF] Killingsworth | Significant Human Impact on the Flux and δ 34S of Sulfate from the Largest River in North America[END_REF] used an isotope mixing model to estimate a pre-anthropogenic δ 34 S = -5‰, and a concentration of sulfate of 115 µM. This pre-anthropogenic estimate is lower than the measured δ 34 S = -0.72‰ and sulfate concentration of 462 µM from this study (which is similar to the isotope and concentration measurements from 2009-2013 from[START_REF] Killingsworth | Significant Human Impact on the Flux and δ 34S of Sulfate from the Largest River in North America[END_REF]).

Note that the most polluted rivers (see text) have been excluded from these calculations and figures.

Acknowledgements

This research was funded by a Foster and Coco Stanback postdoctoral fellowship and a Marie Curie Career Integration Grant (CIG14-631752) to AB.

JFA acknowledges the support of NSF-OCE grant 1340174 and NSF-EAR grant 1349858. WF acknowledges the support of a grant from the David and Lucile Packard Foundation.