Hydrogeochemistry and strontium isotopes in the Arno River Basin (Tuscany, Italy): Constraints on natural controls by statistical modeling

• Published in: Journal of hydrology (Amsterdam) Vol. 360; no. 1; pp. 166 - 183
• Main Authors: Nisi, B., Buccianti, A., Vaselli, O., Perini, G., Tassi, F., Minissale, A., Montegrossi, G.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 15-10-2008; [Amsterdam; New York]: Elsevier; Elsevier Science
• Subjects: Arno River Basin; Central Italy; Earth sciences; Earth, ocean, space; Exact sciences and technology; Freshwater; Geochemistry; Hydrology; Hydrology. Hydrogeology; Inverse modeling; Isotope geochemistry; Isotope geochemistry. Geochronology; Mineralogy; Running waters; Silicates; Strontium isotopes; Water geochemistry; Southern Europe; Italy; Arno River Basin; Europe; Tuscany Italy; Italy, Arno R; Italy, Tuscany; Water geochemistry; Arno River Basin; Running waters; Central Italy; Strontium isotopes; Inverse modeling; tributaries; springs; mixing; bromine; inverse modeling; sulfates; nitrates; drainage basins; trace elements; hydrochemistry; rubidium; thermal waters; Sr-87/Sr-86; fluorine; stable isotopes; dissolution; least-squares; concentration; Rb/Sr; chlorine; water-rock interaction; statistical analysis; chemical composition
• ISSN: 0022-1694; 1879-2707
• DOI: 10.1016/j.jhydrol.2008.07.030

In this paper the chemistry of major and trace elements and, for the first time, strontium isotopic ratios measured in running waters from the Arno River Basin (Tuscany, central-northern Italy) and thermal springs discharging in the same hydrographic system are presented and discussed. Classical graphical methods (e.g. mixing diagrams) have here been improved to identify, in a correct statistical sample space, extreme chemical compositions attributable to the action of geochemical processes and/or inherited from specific lithologies (namely contributions or components 1, 2 and 3) to be used in inverse modeling procedures, due to the absence of clear end-members. A linear least squares problem, with non-negativity constraints and distances, as required for compositional data (convex linear mixing) was solved by considering the contribution of the most discriminant variables given by the
87Sr/
86Sr isotopic ratios and the concentration of Ca
2+, Mg
2+, Sr and Rb. Following such approach, it can be assumed (
p
<
0.05) that component 1, characterized by a Ca
2+(Mg
2+)–
HCO
3
-
radiogenic-rich (
87Sr/
86Sr
=
0.71274; Rb/Sr
=
0.039), represents a dominant feature at basin scale, with a weight ranging from 69% to 100%. Much lower percentages are related to component 2, represented by a Ca
2+(Mg
2+)–
SO
4
2
-
facies with intermediate
87Sr/
86Sr (0.70874) and low Rb/Sr (2.8
×
10
−4) ratios and component 3, identified by Ca
2+–
HCO
3
-
facies with less radiogenic
87Sr/
86Sr (0.70827) and low Mg
2+/Ca
2+ (0.011) ratios. These components are mainly dominated by the dissolution of evaporitic rocks and/or mixing with thermal waters in the southern part of the Arno River Basin and by dissolution of the carbonatic fraction, kinetically favored with respect to that of the silicatic minerals, in the upper reaches of the main course and its tributaries, respectively.