Calcite precipitation instability under laminar, open-channel flow

We present a 2D numerical model for the growth of calcite from supersaturated aqueous solutions under laminar, open-channel flow conditions. The model couples solution chemistry, precipitation at solution/calcite interfaces, hydrodynamics, diffusion and degassing. The model output is compared with experimental results obtained using an oversaturated calcite solution produced by mixing CaCl₂ and Na₂CO₃. The precipitation rate is observed to increase when the supersaturated solution flows over an obstruction, leading to a growth instability that causes the formation of terraces. At relatively high flow rates, the most important mechanism for this behaviour seems to be hydrodynamic advection of dissolved species either towards or away from the calcite surface, depending on location relative to the obstruction, which deforms the concentration gradients. At lower flow rates, steepening of diffusion gradients around protrusions becomes important. Enhanced degassing over the obstruction due to shallowing and pressure drop is not important on small scales. Diffusion controlled transport close to the calcite surface can lead to a fingering-type growth instability, which generates porous textures. Our results are consistent with existing diffusive boundary layer theory, but for flow over non-smooth surfaces, simple calcite precipitation models that include empirical correlations between fluid flow rate and calcite precipitation rate are inaccurate.     

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.     

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.     

Modeling households’ decisions on reconstruction of houses damaged by earthquakes––Japanese case study

In this study, households’ decisions on reconstruction of damaged houses were modeled, using questionnaire data in Japan. Characteristics of households’ decisions were investigated using parameter estimation results. The effects of subsidizing policy were evaluated as follows. First, subsidy for rebuilding may be effective for the households whose houses were heavily damaged. Second, there is a possibility that subsidy accelerated rebuilding of houses by the households having children. Third, subsidy for rebuilding may not be effective for elderly people’s households and households in depopulated areas.