The role of sulfur in chemical weathering and atmospheric CO2 fluxes: Evidence from major ions, δCDIC, and δSSO4 in rivers of the Canadian Cordillera

Water samples from the Fraser, Skeena and Nass River basins of the Canadian Cordillera were analyzed for dissolved major element concentrations (HCO₃⁻, SO₄²⁻, Cl⁻, Ca²⁺, Mg²⁺, K⁺, Na⁺), δ¹³C of dissolved inorganic carbon (δ¹³C_DIC), and δ³⁴S of dissolved sulfate (δ³⁴S_SO4) to quantify chemical weathering rates and exchanges of CO₂ between the atmosphere, hydrosphere, and lithosphere. Weathering rates of silicates and carbonates were determined from major element mass balance. Combining the major element mass balance with δ³⁴S_SO4 (−8.9 to 14.1‰_CDT) indicates sulfide oxidation (sulfuric acid production) and subsequent weathering of carbonate and to a lesser degree silicate minerals are important processes in the study area. We determine that on average, 81% of the riverine sulfate can be attributed to sulfide oxidation in the Cordilleran rivers, and that 25% of the total weathering cation flux can be attributed to carbonate and silicate dissolution by sulfuric acid. This result is validated by δ¹³C_DIC values (−9.8 to −3.7‰ _VPDB) which represents a mixture of DIC produced by the following weathering pathways: (i) carbonate dissolution by carbonic acid (−8.25‰) > (ii) silicate dissolution by carbonic acid (−17‰) ≈ (iii) carbonate dissolution by sulfuric acid derived from the oxidation of sulfides (coupled sulfide-carbonate weathering) (+0.5‰).δ³⁴S_SO4 is negatively correlated with δ¹³C_DIC in the Cordilleran rivers, which further supports the hypothesis that sulfuric acid produced by sulfide oxidation is primarily neutralized by carbonates, and that sulfide-carbonate weathering impacts the δ¹³C_DIC of rivers. The negative correlation between δ³⁴S_SO4 and δ¹³C_DIC is not observed in the Ottawa and St. Lawrence River basins. This suggests other factors such as landscape age (governed by tectonic uplift) and bedrock geology are important controls on regional sulfide oxidation rates, and therefore also on the magnitude of sulfide-carbonate weathering—i.e., it is more significant in tectonically active areas.Calculated DIC fluxes due to Ca and Mg silicate weathering by carbonic acid (38.3 × 10³ mol C · km⁻² · yr⁻¹) are similar in magnitude to DIC fluxes due to sulfide-carbonate weathering (18.5 × 10³ mol C · km⁻² · yr⁻¹). While Ca and Mg silicate weathering facilitates a transfer of atmospheric CO₂ to carbonate rocks, sulfide-carbonate weathering can liberate CO₂ from carbonate rocks to the atmosphere when sulfide oxidation exceeds sulfide deposition. This implies that in the Canadian Cordillera, sulfide-carbonate weathering can offset up to 48% of the current CO₂ drawdown by silicate weathering in the region.     

Geochemistry of volcanic and hydrothermal fluids and volatile budget of the Kamchatka-Kuril subduction zone

A data base for the composition and emission rates of more than 100 thermal manifestations including boiling geothermal systems and 23 volcanoes along the 1900 km long Kamchatka-Kuril (KK) arc is presented. These results were used to estimate mean fluxes of volatiles from the KK arc. The fluxes from the KK arc are compared with the fluxes from the best studied Central American (CA) arc and with the compiled literature data on global fluxes. The error ranges and the OUT/IN (in)balance calculations are also discussed. The estimated fluxes of volatiles from volcanic fumaroles and the observed, normalized to the Cl content, fluxes from hydrothermal systems are very close, with the higher hydrothermal flux from Kuril Islands due to a larger number of the acidic Cl-SO₄ springs on the Islands and their outflow rates. The total volcanic SO₂ flux from the whole KK arc is estimated to be higher than 3000 t/d. The measured S and C fluxes from hydrothermal systems are much lower than the volcanic output due to the loss of these components in the upper crust (mineral precipitation). The Cl/³He ratio is inferred to be a stable indicator of the arc setting for hydrothermal and volcanic fluids with a mean value of (2 ± 4) × 10⁹. Comparison of the obtained volcano-hydrothermal fluxes with fluxes calculated from the erupted solid volcanic products at Kamchatka and Kurils during Holocene time reveals that the total estimated volatile output from the KK arc is compatible with the total magmatic output if the intruded to erupted ratio is close to 7, i.e. almost the same as assumed for the Central American arc. Calculated fluxes as well as the ratios for OUT/IN fluxes (volcanic + hydrothermal output/slab + mantle input) for CO₂, S, H₂O, Cl, N₂, ⁴He and ³He from the KK arc normalized to the arc length are in general close to the global estimates. The fractions of CO₂ and S in the total volatile output at KK arc derived directly from the mantle wedge are 18% and 16% (mole basis), respectively. Fractions of mantle derived H₂O, N₂ and Cl are much lower, less that 5% of their output.     

Impact of nitrogenous fertilizers on carbonate dissolution in small agricultural catchments: Implications for weathering CO2 uptake at regional and global scales

The goal of this study was to highlight the occurrence of an additional proton-promoted weathering pathway of carbonate rocks in agricultural areas where N-fertilizers are extensively spread, and to estimate its consequences on riverine alkalinity and uptake of CO₂ by weathering. We surveyed 25 small streams in the calcareous molassic Gascogne area located in the Garonne river basin (south-western France) that drain cultivated or forested catchments for their major element compositions during different hydrologic periods. Among these catchments, the Hay and the Montoussé, two experimental catchments, were monitored on a weekly basis. Studies in the literature from other small carbonate catchments in Europe were dissected in the same way. In areas of intensive agriculture, the molar ratio (Ca + Mg)/HCO₃ in surface waters is significantly higher (0.7 on average) than in areas of low anthropogenic pressure (0.5). This corresponds to a decrease in riverine alkalinity, which can reach 80% during storm events. This relative loss of alkalinity correlates well with the content in surface waters. In cultivated areas, the contribution of atmospheric/soil CO₂ to the total riverine alkalinity (CO_{2 ATM-SOIL}/HCO₃) is less than 50% (expected value for carbonate basins), and it decreases when the nitrate concentration increases. This loss of alkalinity can be attributed to the substitution of carbonic acid (natural weathering pathway) by protons produced by nitrification of N-fertilizers (anthropogenic weathering pathway) occurring in soils during carbonate dissolution. As a consequence of these processes, the alkalinity over the last 30 years shows a decreasing trend in the Save river (one of the main Garonne river tributaries, draining an agricultural catchment), while the nitrate and calcium plus magnesium contents are increasing.We estimated that the contribution of atmospheric/soil CO₂ to riverine alkalinity decreased by about 7-17% on average for all the studied catchments. Using these values, the deficit of CO₂ uptake can be estimated as up to 0.22-0.53 and 12-29 Tg¹ yr⁻¹ CO₂ on a country scale (France) and a global scale, respectively. These losses represent up to 5.7-13.4% and only 1.6-3.8% of the total CO₂ flux naturally consumed by carbonate dissolution, for France and on a global scale, respectively. Nevertheless, this loss of alkalinity relative to the Ca + Mg content relates to carbonate weathering by protons from N-fertilizers nitrification, which is a net source of CO₂ for the atmosphere. This anthropogenic CO₂ source is not negligible since it could reach 6-15% of CO₂ uptake by natural silicate weathering and could consequently partly counterbalance this natural CO₂ sink.     

Chemical weathering, river geochemistry and atmospheric carbon fluxes from volcanic and ultramafic regions on Luzon Island, the Philippines

We investigated rates of chemical weathering of volcanic and ophiolitic rocks on Luzon Island, the Philippines. Luzon has a tropical climate and is volcanically and tectonically very active, all factors that should enhance chemical weathering. Seventy-five rivers and streams (10 draining ophiolites, 65 draining volcanic bedrock) and two volcanic hot springs were sampled and analyzed for major elements, alkalinity and ⁸⁷Sr/⁸⁶Sr. Cationic fluxes from the volcanic basins are dominated by Ca²⁺ and Mg²⁺ and dissolved silica concentrations are high (500-1900 μM). Silica concentrations in streams draining ophiolites are lower (400-900 μM), and the cationic charge is mostly Mg²⁺. The areally weighted average CO₂ export flux from our study area is 3.89 ± 0.21 × 10⁶ mol/km²/yr, or 5.99 ± 0.64 × 10⁶ mol/km²/yr from ophiolites and 3.58 ± 0.23 × 10⁶ mol/km²/yr from volcanic areas (uncertainty given as ±1 standard error, s.e.). This is ∼6-10 times higher than the current best estimate of areally averaged global CO₂ export by basalt chemical weathering and ∼2-3 times higher than the current best estimate of CO₂ export by basalt chemical weathering in the tropics. Extrapolating our findings to all tropical arcs, we estimate that around one tenth of all atmospheric carbon exported via silicate weathering to the oceans annually is processed in these environments, which amount to ∼1% of the global exorheic drainage area. Chemical weathering of volcanic terranes in the tropics appears to make a disproportionately large impact on the long-term carbon cycle.     

Chemical weathering processes in the Yalong River draining the eastern Tibetan Plateau,China

To better understand chemical weathering and controlling processes in the Yalong River of the eastern Tibetan Plateau, this study presents major ion concentrations and stable isotopes of the dissolved loads. The isotopic compositions (δ¹³C-DIC, δ³⁴S and δ¹⁸O-SO₄) of the dissolved loads are very useful to quantify solute sources and define the carbon budget related with chemical weathering in riverine systems. The isotopic composition of sulphate demonstrates that most of the sulphate is derived from sulphide oxidation, particularly in the upper reach of the Yalong River. The correlations between δ¹³C-DIC, water chemistry and isotopes of sulphate, suggest that the carbon dynamics are mainly affected by carbonate weathering by sulphuric acid and equilibration processes. Approximately 13% of the dissolved inorganic carbon in the Yalong River originates from carbonate weathering by strong acid. The CO₂ consumption rates are estimated to be 2.8 × 10⁵ mol/km²/yr and 0.9 × 10⁵ mol/km²/yr via carbonate and silicate weathering in the Yalong River, respectively. In this study, the influence of sulphide oxidation and metamorphic CO₂ on the carbon budget is estimated for the Yalong River draining the eastern Tibetan Plateau.     

Silicate weathering in anoxic marine sediments

Two sediment cores retrieved at the northern slope of Sakhalin Island, Sea of Okhotsk, were analyzed for biogenic opal, organic carbon, carbonate, sulfur, major element concentrations, mineral contents, and dissolved substances including nutrients, sulfate, methane, major cations, humic substances, and total alkalinity. Down-core trends in mineral abundance suggest that plagioclase feldspars and other reactive silicate phases (olivine, pyroxene, volcanic ash) are transformed into smectite in the methanogenic sediment sections. The element ratios Na/Al, Mg/Al, and Ca/Al in the solid phase decrease with sediment depth indicating a loss of mobile cations with depth and producing a significant down-core increase in the chemical index of alteration. Pore waters separated from the sediment cores are highly enriched in dissolved magnesium, total alkalinity, humic substances, and boron. The high contents of dissolved organic carbon in the deeper methanogenic sediment sections (50-150 mg dm⁻³) may promote the dissolution of silicate phases through complexation of Al³⁺ and other structure-building cations. A non-steady state transport-reaction model was developed and applied to evaluate the down-core trends observed in the solid and dissolved phases. Dissolved Mg and total alkalinity were used to track the in-situ rates of marine silicate weathering since thermodynamic equilibrium calculations showed that these tracers are not affected by ion exchange processes with sediment surfaces. The modeling showed that silicate weathering is limited to the deeper methanogenic sediment section whereas reverse weathering was the dominant process in the overlying surface sediments. Depth-integrated rates of marine silicate weathering in methanogenic sediments derived from the model (81.4-99.2 mmol CO₂ m⁻² year⁻¹) are lower than the marine weathering rates calculated from the solid phase data (198-245 mmol CO₂ m⁻² year⁻¹) suggesting a decrease in marine weathering over time. The production of CO₂ through reverse weathering in surface sediments (4.22-15.0 mmol CO₂ m⁻² year⁻¹) is about one order of magnitude smaller than the weathering-induced CO₂ consumption in the underlying sediments. The evaluation of pore water data from other continental margin sites shows that silicate weathering is a common process in methanogenic sediments. The global rate of CO₂ consumption through marine silicate weathering estimated here as 5-20 Tmol CO₂ year⁻¹ is as high as the global rate of continental silicate weathering.     

Impact of nitrogen fertilizers on the natural weathering-erosion processes and fluvial transport in the Garonne basin

Knowledge of the impact of N-fertilizers on the weathering-erosion processes of soils in intensively cultivated regions is of prime importance. Nitrification of NH₄⁻ fertilizers produces HNO₃ in the basin of the Garonne river, enhancing soil degradation. Their influence on the weathering rates was determined by calculating the consumption rate of atmospheric/soil CO₂ by soil weathering and erosion, and its contribution to the total dissolved riverine HCO₃⁻. This contribution was found to be less than 50% which corresponds normally to a complete carbonate dissolution by carbonic acid, suggesting that part of the alkalinity in the river waters is due to carbonate dissolution by an acid other than carbonic acid, probably HNO₃.     

River dissolved and solid loads in the Lesser Antilles: New insight into basalt weathering processes

We present here the first available estimations of chemical weathering and associated atmospheric CO₂ consumption rates as well as mechanical erosion rate for the Lesser Antilles. The chemical weathering (100–120 t/km²/year) and CO₂ consumption (1.1–1.4 × 10⁶ mol/km²/year) rates are calculated after subtraction of the atmospheric and hydrothermal inputs in the chemical composition of the river dissolved loads. These rates thus reflect only the low-temperature basalt weathering. Mechanical erosion rates (approx. 800–4000 t/km²/year) are estimated by a geochemical mass balance between the dissolved and solid loads and mean unaltered rock. The calculated chemical weathering rates and associated atmospheric CO₂ consumption rates are among the highest values worldwide but are still lower than those of other tropical volcanic islands and do not fit with the HCO₃⁻ concentration vs. 1/T correlation proposed by Dessert et al. (2001). The thick soils and explosive volcanism context of the Lesser Antilles are the two possible keys to this different weathering behaviour; the development of thick soils limits the chemical weathering and the presence of very porous pyroclastic flows allows an important water infiltration and thus subsurface weathering mechanisms, which are less effective for atmospheric CO₂ consumption.     

Volcanic and anthropogenic contributions to global weathering budgets

We evaluate whether the global weathering budget is near steady state for the pre-anthropogenic modern environment by assessing the magnitude of acidity-generating volcanic exhalations. The weathering rate induced by volcanic acid fluxes, of which the CO₂ flux is the most important, can be expressed as an average release rate of dissolved silica, based on a model feldspar-weathering scheme, and the ratio of carbonate-to-silicate rock weathering. The theoretically predicted flux of silica from chemical weathering is slightly smaller than the estimated global riverine silica flux. After adjustment for carbonate weathering, the riverine dissolved bicarbonate flux is larger than the volcanic carbon degassing rate by a factor of about three. There are substantial uncertainties associated with the calculated and observed flux values, but the modern system may either not be in steady state, or additional, “unknown” carbon sources may exist. The closure errors in the predicted budgets and observed riverine fluxes suggest that continental weathering rates might have had an impact on atmospheric CO₂ levels at a time scale of 10³-10⁴ years, and that enhanced weathering rates during glacial periods might have been a factor in the reduced glacial atmospheric CO₂ levels. Recent anthropogenic emissions of carbon and sulfur have a much larger acid-generating capacity than the natural fluxes. Estimated potential weathering budgets to neutralize these fluxes are far in excess of observed values. A theoretical scenario for a return to steady state at the current anthropogenic acidity emissions (disregarding the temporary buffering action of the ocean reservoir) requires either significantly lower pH values in continental surface waters as a result of storage of strong acids, and/or higher temperatures as a result of enhanced atmospheric CO₂ levels in order to create weathering rates that can neutralize the total flux of anthropogenic and natural background acidity.     

The CO2 system in rivers of the Australian Victorian Alps: CO2 evasion in relation to system metabolism and rock weathering on multi-annual time scales

The patterns of dissolved inorganic C (DIC) and aqueous CO₂ in rivers and estuaries sampled during summer and winter in the Australian Victorian Alps were examined. Together with historical (1978–1990) geochemical data, this study provides, for the first time, a multi-annual coverage of the linkage between CO₂ release via wetland evasion and CO₂ consumption via combined carbonate and aluminosilicate weathering. δ¹³C values imply that carbonate weathering contributes ∼36% of the DIC in the rivers although carbonates comprise less than 5% of the study area. Baseflow/interflow flushing of respired C3 plant detritus accounts for ∼50% and atmospheric precipitation accounts for ∼14% of the DIC. The influence of in river respiration and photosynthesis on the DIC concentrations is negligible. River waters are supersaturated with CO₂ and evade ∼27.7 × 10⁶ mol/km²/a to ∼70.9 × 10⁶ mol/km²/a CO₂ to the atmosphere with the highest values in the low runoff rivers. This is slightly higher than the global average reflecting higher gas transfer velocities due to high wind speeds. Evaded CO₂ is not balanced by CO₂ consumption via combined carbonate and aluminosilicate weathering which implies that chemical weathering does not significantly neutralize respiration derived H₂CO₃. The results of this study have implications for global assessments of chemical weathering yields in river systems draining passive margin terrains as high respiration derived DIC concentrations are not directly connected to high carbonate and aluminosilicate weathering rates.     

Understanding the Cyclicity of Chemical Weathering and Associated CO<Subscript>2</Subscript> Consumption in the Brahmaputra River Basin (India): The Role of Major Rivers in Climate Change Mitigation Perspective

Weathering of rocks that regulate the water chemistry of the river has been used to evaluate the CO₂ consumption rate which exerts a strong influence on the global climate. The foremost objective of the present research is to estimate the chemical weathering rate (CWR) of the continental water in the entire stretch of Brahmaputra River from upstream to downstream and their associated CO₂ consumption rate. To establish the link between the rapid chemical weathering and thereby enhance CO₂ drawdown from the atmosphere, the major ion composition of the Brahmaputra River that drains the Himalaya has been obtained. Major ion chemistry of the Brahmaputra River was resolved on samples collected from nine locations in pre-monsoon, monsoon and post-monsoon seasons for two cycles: cycle I (2011–2012) and cycle II (2013–2014). The physico-chemical parameters of water samples were analysed by employing standard methods. The Brahmaputra River was characterized by alkalinity, high concentration of Ca²⁺ and HCO₃ ^? along with significant temporal variation in major ion composition. In general, it was found that water chemistry of the river was mainly controlled by rock weathering with minor contributions from atmospheric and anthropogenic sources. The effective CO₂ pressure (log\({{\text{P}}_{{\text{C}}{{\text{O}}_{\text{2}}}}}\)) for pre-monsoon, monsoon and post-monsoon has been estimated. The question of rates of chemical weathering (carbonate and silicate) was addressed by using TDS and run-off (mm year^?1). It has been found that the extent of CWR is directly dependent on the CO₂ consumption rate which may be further evaluated from the perspective of climate change mitigation The average annual CO₂ consumption rate of the Brahmaputra River due to silicate and carbonate weathering was found to be 0.52 (×10⁶ mol Km^?2 year^?1) and 0.55 (×10⁶ mol Km^?2 year^?1) for cycle I and 0.49 (×10⁶ mol Km^?2 year^?1) and 0.52 (×10⁶ mol Km^?2 year^?1) for cycle II, respectively, which were significantly higher than that of other Himalayan rivers. Estimation of CWR of the Brahmaputra River indicates that carbonate weathering largely dominates the water chemistry of the Brahmaputra River.     

Geochemistry of the dissolved load of the Changjiang Basin rivers: Anthropogenic impacts and chemical weathering

This study focuses on the chemical and Sr isotopic compositions of the dissolved load of the rivers of the Changjiang Basin, one of the largest riverine systems in the world. Water samples were collected in August 2006 from the main tributaries and the main Changjiang channel. The chemical and isotopic analyses indicated that four major reservoirs (carbonates, silicates, evaporites and agriculture/urban effluents) contribute to the total dissolved solutes. The overall chemical weathering (carbonate and silicate) rate for the Changjiang is approximately 40 ton/km²/year or 19 mm/kyr, similar to that of the Ganges-Brahmaputra system, and the basin is characterized by carbonate and silicate weathering rates ranging from 17 to 56 ton/km²/year and from 0.7 to 7.1 ton/km²/year, respectively. In the lower reach of the Changjiang main channel, the weathering rates are estimated to be 36 and 2.2 ton/km²/year for carbonates and silicates, respectively. It appears that sulphuric acid may dominate chemical weathering reactions for some sub-basins. The budgets of CO₂ consumption are estimated to be 646 × 10⁹ and 191 × 10⁹ mol/year by carbonate and silicate weathering, respectively. The contribution of the anthropogenic inputs to the cationic TDS of the Changjiang is estimated to be 15-20% for the most downstream stations. Our study suggested that the Changjiang is strongly impacted by human activities and is very sensitive to the change of land use.     

Redox potential (Eh) and anion effects of pyrite (FeS2) leaching at pH 1

Pyrite plays the central role in the environmental issue of acid rock drainage. Natural weathering of pyrite results in the release of sulphuric acid which can lead to further leaching of heavy and toxic metals from other associated minerals. Understanding how pyrite reacts in aqueous solution is critical to understanding the natural weathering processes undergone by this mineral. To this end an investigation of the effect of solution redox potential (Eh) and various anions on the rate of pyrite leaching under carefully controlled conditions has been undertaken.Leaching of pyrite has been shown to proceed significantly faster at solution Eh of 900 mV (SHE) than at 700 mV, at pH 1, for the leach media of HCl, H₂SO₄ and HClO₄. The predominant effect of Eh suggests electrochemical control of pyrite leaching with similar mechanism(s) at Eh of 700 and 900 mV albeit with different kinetics. Leach rates at 700 mV were found to decrease according to HClO₄ > HCl > H₂SO₄ while at 900 mV the leach rate order was HCl > HClO₄ > H₂SO₄. Solution Fe³⁺ activity is found to continually increase during all leaches; however, this is not accompanied by an increase in leach rate.Synchrotron based photoemission electron microscopy (PEEM) measurements showed a localised distribution of adsorbed and oxidised surface species highlighting that pyrite oxidation and leaching is a highly site specific process mediated by adsorption of oxidants onto specific surface sites. It appears that rates may be controlled, in part, by the propensity of acidic anions to bind to the surface, which varies according to , thus reducing the reactive or effective surface area. However, anions may also be involved in specific reactions with surface leach products. Stoichiometric dissolution data (Fe/S ratio), XPS and XRD data indicate that the highest leach rates (in HCl media at 900 mV Eh) correlate with relatively lower surface S abundance. Furthermore, there are indications that solution Cl⁻ assists oxidation especially at higher Eh through the prevention of surface S⁰ buildup at reactive surface sites.     

Stream geochemistry, chemical weathering and CO2 consumption potential of andesitic terrains, Dominica, Lesser Antilles

Recent studies of chemical weathering of andesitic-dacitic material on high-standing islands (HSIs) have shown these terrains have some of the highest observed rates of chemical weathering and associated CO₂ consumption yet reported. However, the paucity of stream gauge data in many of these terrains has limited determination of chemical weathering product fluxes. In July 2006 and March 2008, stream water samples were collected and manual stream gauging was performed in watersheds throughout the volcanic island of Dominica in the Lesser Antilles. Distinct wet and dry season solute concentrations reveal the importance of seasonal variations on the weathering signal. A cluster analysis of the stream geochemical data shows the importance of parent material age on the overall delivery of solutes. Observed Ca:Na, HCO₃:Na and Mg:Na ratios suggest crystallinity of the parent material may also play an important role in determining weathering fluxes. From total dissolved solids concentrations and mean annual discharge calculations we calculate chemical weathering yields of (6-106 t km⁻² a⁻¹), which are similar to those previously determined for basalt terrains. Silicate fluxes (3.1-55.4 t km⁻² a⁻¹) and associated CO₂ consumption (190-1575 × 10³ mol km⁻² a⁻¹) determined from our study are among the highest determined to date. The calculated chemical fluxes from our study confirm the weathering potential of andesitic-dacitic terrains and that additional studies of these terrains are warranted.     

GEOCARBSULF: A combined model for Phanerozoic atmospheric O2 and CO2

A model for the combined long-term cycles of carbon and sulfur has been constructed which combines all the factors modifying weathering and degassing of the GEOCARB III model [Berner R.A., Kothavala Z., 2001. GEOCARB III: a revised model of atmospheric CO₂ over Phanerozoic time. Am. J. Sci. 301, 182-204] for CO₂ with rapid recycling and oxygen dependent carbon and sulfur isotope fractionation of an isotope mass balance model for O₂ [Berner R.A., 2001. Modeling atmospheric O₂ over Phanerozoic time. Geochim. Cosmochim. Acta65, 685-694]. New isotopic data for both carbon and sulfur are used and new feedbacks are created by combining the models. Sensitivity analysis is done by determining (1) the effect on weathering rates of using rapid recycling (rapid recycling treats carbon and sulfur weathering in terms of young rapidly weathering rocks and older more slowly weathering rocks); (2) the effect on O₂ of using different initial starting conditions; (3) the effect on O₂ of using different data for carbon isotope fractionation during photosynthesis and alternative values of oceanic δ¹³C for the past 200 million years; (4) the effect on sulfur isotope fractionation and on O₂ of varying the size of O₂ feedback during sedimentary pyrite formation; (5) the effect on O₂ of varying the dependence of organic matter and pyrite weathering on tectonic uplift plus erosion, and the degree of exposure of coastal lands by sea level change; (6) the effect on CO₂ of adding the variability of volcanic rock weathering over time [Berner, R.A., 2006. Inclusion of the weathering of volcanic rocks in the GEOCARBSULF model. Am. J. Sci.306 (in press)]. Results show a similar trend of atmospheric CO₂ over the Phanerozoic to the results of GEOCARB III, but with some differences during the early Paleozoic and, for variable volcanic rock weathering, lower CO₂ values during the Mesozoic. Atmospheric oxygen shows a major broad late Paleozoic peak with a maximum value of about 30% O₂ in the Permian, a secondary less-broad peak centered near the Silurian/Devonian boundary, variation between 15% and 20% O₂ during the Cambrian and Ordovician, a very sharp drop from 30% to 15% O₂ at the Permo-Triassic boundary, and a more-or less continuous rise in O₂ from the late Triassic to the present.     

Contribution of {\text{P}}_{{{\text{CO}}_{ 2} {\text{eq}}}} and 13CTDIC Evaluation to the Identification of CO2 Sources in Volcanic Groundwater Systems: Influence of Hydrometeorological Conditions and Lava Flow Morphologies—Application to the Argnat Basin (Chaîne des Puys, Massif Central, France)

Mineralization of groundwater in volcanic aquifers is partly acquired through silicates weathering. This alteration depends on the dissolution of atmospheric, biogenic, or mantellic gaseous CO₂ whose contributions may depend on substratum geology, surface features, and lava flow hydrological functionings. Investigations of $ {\text{P}}_{{{\text{CO}}_{ 2} {\text{eq}}}} $ and δ¹³C_TDIC (total dissolved inorganic carbon) on various spatiotemporal scales in the unsaturated and saturated zones of volcanic flows of the Argnat basin (French Massif Central) have been carried out to identify the carbon sources in the system. Mantellic sources are related to faults promoting CO₂ uplift from the mantle to the saturated zone. The contribution of this source is counterbalanced by infiltration of water through the unsaturated zone, accompanied by dissolution of soil CO₂ or even atmospheric CO₂ during cold periods. Monitoring and modeling of δ¹³C_TDIC in the unsaturated zone shows that both $ {\text{P}}_{{{\text{CO}}_{ 2} {\text{eq}}}} $ and δ¹³C_TDIC are controlled by air temperature which influences soil respiration and soil-atmosphere CO₂ exchanges. The internal geometry of volcanic lava flows controls water patterns from the unsaturated zone to saturated zone and thus may explain δ¹³C heterogeneity in the saturated zone at the basin scale.