Trace Element Partitioning and Accessory Phase Saturation during H2O-Saturated Melting of Basalt with Implications for Subduction Zone Chemical Fluxes

Experimental phase equilibrium and trace element partitioningdata are reported for H₂O-saturated mid-ocean ridge basalt at2·5 GPa, 750–900°C and oxygen fugacities atthe nickel–nickel oxide buffer. Garnet, omphacite andrutile are present at all temperatures. Amphibole and epidotedisappear as residual phases above 800°C; allanite appearsabove 750°C. The Na–Al-rich silicate glass presentin all run products is likely to have quenched from a supercriticalliquid. Trace element analyses of glasses demonstrate the importantcontrol exerted by residual minerals on liquid chemistry. Inaddition to garnet, which controls heavy rare earth elements(HREE) and Sc, and rutile, which controls Ti, Nb and Ta, allanitebuffers the light REE (LREE; La–Sm) contents of liquidsto relatively low levels and preferentially holds back Th relativeto U. In agreement with previous experimental and metamorphicstudies we propose that residual allanite plays a key role inselectively retaining trace elements in the slab during subduction.Experimental data and analyses of allanite-bearing volcanicrocks are used to derive a model for allanite solubility inliquids as a function of pressure, temperature, anhydrous liquidcomposition and LREE content. The large temperature dependenceof allanite solubility is very similar to that previously determinedfor monazite. Our model, fitted to 48 datapoints, retrievesLREE solubility (in ppm) to within a factor of 1· 40over a pressure range of 0–4 GPa, temperature range of700–1200°C and for liquids with anhydrous SiO₂ contentsof 50–84 wt %. This uncertainty in LREE content is equivalentto a temperature uncertainty of only ± 27°C at 1000K, indicating the potential of allanite as a geothermometer.Silicic liquids from either basaltic or sedimentary protolithswill be saturated in allanite except for Ca-poor protolithsor at very high temperatures. For conventional subduction geothermsthe low solubility of LREE (+ Th) in liquids raises questionsabout the mechanism of LREE + Th transport from slab to wedge.It is suggested either that, locally, temperatures experiencedby the slab are high enough to eliminate allanite in the residueor that substantial volumes of H₂O-rich fluids must pass throughthe mantle wedge prior to melting. The solubility of accessoryphases in fluids derived from subducted rocks can provide importantconstraints on subduction zone thermal structure. KEY WORDS: subduction; experimental petrology; allanite; solubility; supercritical liquid; eclogite     

Origin of the Oceanic Lithosphere

In a global examination of mid-ocean ridge basalt (MORB) glasscompositions, we find that Na₈–Fe₈–depth variationsdo not support modeling of MORBs as aggregates of melt compositionsgenerated over a large range of temperature and pressure. However,the Na₈–Fe₈ variations are consistent with the compositionalsystematics of solidus melts in the plagioclase–spinellherzolite transition in the CaO–MgO–Al₂O₃–SiO₂–Na₂O–FeO(CMASNF) system. For natural compositions, the P–T rangefor melt extraction is estimated to be 1·2–1·5GPa and 1250–1280°C. This P–T range is a closematch with the maximum P–T conditions for explosive pressure-releasevaporization of carbonate-bearing melts. It is proposed thatfracturing of the lithosphere induces explosive formation andescape of CO₂ vapor. This provides the vehicle for extractionof MORBs at a relatively uniform T and P. The upper portionof the CO₂-bearing and slightly melted seismic low-velocityzone flows toward the ridge, rises at the ridge axis to lower-lithospheredepths, melts much more extensively during this rise, and releasesMORB melts to the surface driven by explosively escaping CO₂vapor. The residue and overlying crust produced by this meltingthen migrate away from the ridge axis as new oceanic lithosphere.The entire process of oceanic lithosphere creation involvesonly the upper 140 km. When lithospheric stresses shift fractureformation to other localities, escape of CO₂ ceases, the vehiclefor transporting melt to the surface disappears, and ridgesdie. Inverse correlations of Na₈ vs Fe₈ for MORB glasses areexplained by mantle heterogeneity, and positive variations superimposedon the inverse variations are consistent with progressive extractionof melts from short, ascending melting columns. The uniformlylow temperatures of MORB extraction are not consistent withthe existence of hot plumes on or close to ocean ridges. Inthis modeling, the southern Atlantic mantle from Bouvet to about26°N is relatively homogeneous, whereas the Atlantic mantlenorth of about 26°N shows significant long-range heterogeneity.The mantle between the Charlie Gibbs and Jan Mayen fracturezones is strongly enriched in FeO/MgO, perhaps by a trappedfragment of basaltic crust. Iceland is explained as the productof this enrichment, not a hot plume. The East Pacific Rise,Galapagos Ridge, Gorda Ridge, and Juan de Fuca Ridge samplemantle that is heterogeneous over short distances. The mantlebeneath the Red Sea is enriched in FeO/MgO relative to the mantlebeneath the northern Indian Ocean.     

Phase Relations and Melting of Anhydrous K-bearing Eclogite from 1200 to 1600{degrees}C and 3 to 5 GPa

To investigate eclogite melting under mantle conditions, wehave performed a series of piston-cylinder experiments usinga homogeneous synthetic starting material (GA2) that is representativeof altered mid-ocean ridge basalt. Experiments were conductedat pressures of 3·0, 4·0 and 5·0 GPa andover a temperature range of 1200–1600°C. The subsolidusmineralogy of GA2 consists of garnet and clinopyroxene withminor quartz–coesite, rutile and feldspar. Solidus temperaturesare located at 1230°C at 3·0 GPa and 1300°C at5·0 GPa, giving a steep solidus slope of 30–40°C/GPa.Melting intervals are in excess of 200°C and increase withpressure up to 5·0 GPa. At 3·0 GPa feldspar, rutileand quartz are residual phases up to 40°C above the solidus,whereas at higher pressures feldspar and rutile are rapidlymelted out above the solidus. Garnet and clinopyroxene are theonly residual phases once melt fractions exceed 20% and garnetis the sole liquidus phase over the investigated pressure range.With increasing melt fraction garnet and clinopyroxene becomeprogressively more Mg-rich, whereas coexisting melts vary fromK-rich dacites at low degrees of melting to basaltic andesitesat high melt fractions. Increasing pressure tends to increasethe jadeite and Ca-eskolaite components in clinopyroxene andenhance the modal proportion of garnet at low melt fractions,which effects a marked reduction in the Al₂O₃ and Na₂O contentof the melt with pressure. In contrast, the TiO₂ and K₂O contentsof the low-degree melts increase with increasing pressure; thusNa₂O and K₂O behave in a contrasted manner as a function ofpressure. Altered oceanic basalt is an important component ofcrust returned to the mantle via plate subduction, so GA2 maybe representative of one of many different mafic lithologiespresent in the upper mantle. During upwelling of heterogeneousmantle domains, these mafic rock-types may undergo extensivemelting at great depths, because of their low solidus temperaturescompared with mantle peridotite. Melt batches may be highlyvariable in composition depending on the composition and degreeof melting of the source, the depth of melting, and the degreeof magma mixing. Some of the eclogite-derived melts may alsoreact with and refertilize surrounding peridotite, which itselfmay partially melt with further upwelling. Such complex magma-genesisconditions may partly explain the wide spectrum of primitivemagma compositions found within oceanic basalt suites. KEY WORDS: eclogite; experimental petrology; mafic magmatism; mantle melting; oceanic basalts     

Comparison of arsenic concentrations in simultaneously-collected groundwater and aquifer particles from Bangladesh,India, Vietnam,and Nepal

One of the reasons the processes resulting in As release to groundwater in southern Asia remain poorly understood is the high degree of spatial variability of physical and chemical properties in shallow aquifers. In an attempt to overcome this difficulty, a simple device that collects groundwater and sediment as a slurry from precisely the same interval was developed in Bangladesh. Recently published results from Bangladesh and India relying on the needle-sampler are augmented here with new data from 37 intervals of grey aquifer material of likely Holocene age in Vietnam and Nepal. A total of 145 samples of filtered groundwater ranging in depth from 3 to 36 m that were analyzed for As (1–1000 μg/L), Fe (0.01–40 mg/L), Mn (0.2–4 mg/L) and S (0.04–14 mg/L) are compared. The P-extractable (0.01–36 mg/kg) and HCl-extractable As (0.04–36 mg/kg) content of the particulate phase was determined in the same suite of samples, in addition to Fe(II)/Fe ratios (0.2–1.0) in the acid-leachable fraction of the particulate phase. Needle-sampler data from Bangladesh indicated a relationship between dissolved As in groundwater and P-extractable As in the particulate phase that was interpreted as an indication of adsorptive equilibrium, under sufficiently reducing conditions, across 3 orders of magnitude in concentrations according to a distribution coefficient of 4 mL/g. The more recent observations from India, Vietnam and Nepal show groundwater As concentrations that are often an order of magnitude lower at a given level of P-extractable As compared to Bangladesh, even if only the subset of particularly reducing intervals characterized by leachable Fe(II)/Fe >0.5 and dissolved Fe >0.2 mg/L are considered. Without attempting to explain why As appears to be particularly mobile in reducing aquifers of Bangladesh compared to the other regions, the consequences of increasing the distribution coefficient for As between the particulate and dissolved phase to 40 mL/g for the flushing of shallow aquifers of their initial As content are explored.     

Atmospheric mercury near Salmon Falls Creek Reservoir in southern Idaho

Gaseous elemental mercury (GEM) and reactive gaseous mercury (RGM) were measured over 2-week seasonal field campaigns near Salmon Falls Creek Reservoir in south-central Idaho from the summer of 2005 through the fall of 2006 and over the entire summer of 2006 using automated Tekran Hg analyzers. GEM, RGM, and particulate Hg (HgP) were also measured at a secondary site 90 km to the west in southwestern Idaho during the summer of 2006. The study was performed to characterize Hg air concentrations in the southern Idaho area for the first time, estimate Hg dry deposition rates, and investigate the source of observed elevated concentrations. High seasonal variability was observed with the highest GEM (1.91 ± 0.9 ng m⁻³) and RGM (8.1 ± 5.6 pg m⁻³) concentrations occurring in the summer and lower values in the winter (1.32 ± 0.3 ng m⁻³, 3.2 ± 2.9 pg m⁻³ for GEM, RGM, respectively). The summer-average HgP concentrations were generally below detection limit (0.6 ± 1 pg m⁻³). Seasonally averaged deposition velocities calculated using a resistance model were 0.034 ± 0.032, 0.043 ± 0.040, 0.00084 ± 0.0017 and 0.00036 ± 0.0011 cm s⁻¹ for GEM (spring, summer, fall and winter, respectively) and 0.50 ± 0.39, 0.40 ± 0.31, 0.51 ± 0.43 and 0.76 ± 0.57 cm s⁻¹ for RGM. The total annual RGM + GEM dry deposition estimate was calculated to be 11.9 ± 3.3 μg m⁻², or about 2/3 of the total (wet + dry) deposition estimate for the area. Periodic elevated short-term GEM (2.2–12 ng m⁻³) and RGM (50–150 pg m⁻³) events were observed primarily during the warm seasons. Back-trajectory modeling and PSCF analysis indicate predominant source directions to the SE (western Utah, northeastern Nevada) and SW (north-central Nevada) with fewer inputs from the NW (southeastern Oregon and southwestern Idaho).     

The passivation of calcite by acid mine water. Column experiments with ferric sulfate and ferric chloride solutions at pH 2

Column experiments, simulating the behavior of passive treatment systems for acid mine drainage, have been performed. Acid solutions (HCl or H₂SO₄, pH 2), with initial concentrations of Fe(III) ranging from 250 to 1500 mg L⁻¹, were injected into column reactors packed with calcite grains at a constant flow rate. The composition of the solutions was monitored during the experiments. At the end of the experiments (passivation of the columns), the composition and structure of the solids were measured. The dissolution of calcite in the columns caused an increase in pH and the release of Ca into the solution, leading to the precipitation of gypsum and Fe–oxyhydroxysulfates (Fe(III)–SO₄–H⁺ solutions) or Fe–oxyhydroxychlorides (Fe(III)–Cl–H⁺ solutions). The columns worked as an efficient barrier for some time, increasing the pH of the circulating solutions from 2 to 6–7 and removing its metal content. However, after some time (several weeks, depending on the conditions), the columns became chemically inert. The results showed that passivation time increased with decreasing anion and metal content of the solutions. Gypsum was the phase responsible for the passivation of calcite in the experiments with Fe(III)–SO₄–H⁺ solutions. Schwertmannite and goethite appeared as the Fe(III) secondary phases in those experiments. Akaganeite was the phase responsible for the passivation of the system in the experiments with Fe(III)–Cl–H⁺ solutions.     

Geochemical characterization of surface waters and groundwater resources in the Managua area (Nicaragua,Central America)

This paper reports new geochemical data on dissolved major and minor constituents in surface waters and ground waters collected in the Managua region (Nicaragua), and provides a preliminary characterization of the hydrogeochemical processes governing the natural water evolution in this area. The peculiar geological features of the study site, an active tectonic region (Nicaragua Depression) characterized by active volcanism and thermalism, combined with significant anthropogenic pressure, contribute to a complex evolution of water chemistry, which results from the simultaneous action of several geochemical processes such as evaporation, rock leaching, mixing with saline brines of natural or anthropogenic origin. The effect of active thermalism on both surface waters (e.g., saline volcanic lakes) and groundwaters, as a result of mixing with variable proportions of hyper-saline geothermal Na–Cl brines (e.g., Momotombo geothermal plant), accounts for the high salinities and high concentrations of many environmentally-relevant trace elements (As, B, Fe and Mn) in the waters. At the same time the active extensional tectonics of the Managua area favour the interaction with acidic, reduced thermal fluids, followed by extensive leaching of the host rock and the groundwater release of toxic metals (e.g., Ni, Cu). The significant pollution in the area, deriving principally from urban and industrial waste-water, probably also contributes to the aquatic cycling of many trace elements, which attain concentrations above the WHO recommended limits for the elements Ni (∼40 μg/l) and Cu (∼10 μg/l) limiting the potential utilisation of Lake Xolotlan for nearby Managua.     

Phosphorus retention by phlogopite-rich mine tailings

Sewage water rich in P poses a threat to natural water systems. In households that are not connected to a sewer network, an easily maintained purification technique is needed to remove P from domestic waste waters. In the present study, the P retention properties of phlogopite-rich mine tailings were studied to determine whether the untapped tailings can be utilised in wastewater treatment for this purpose. On the basis of its chemical and mineralogical properties, this material, produced in the enrichment process of apatite ore in Siilinjärvi, Finland, can be assumed to act as a P sorbent. The P retention capacity and the P sorption mechanisms of phlogopite-rich tailings were studied with desorption–sorption isotherms and various chemical extractions. The impact of the artificial weathering of the material as well as the reaction time and particle size on P retention were also investigated. The fluctuating shape of the isotherm determined for the untreated material suggests the operation of several sorption mechanisms. The results imply that phlogopite-rich tailings retain P efficiently through specific sorption by amorphous Al and Fe oxides and possibly also by calcite present as an accessory mineral. Acidification of the material, extended reaction time and decrease in particle-size markedly increased the P retention.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.