Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Boron-Bearing Kornerupine from Fiskenaesset, West Greenland: A Reexamination of Specimens from the Type Locality

In 1884, Lorenzen proposed the formula MgAI2SiO6 for his new mineral kornerupine from Fiskenæsset and did not suspect it to contain boron. Lacroix and de Gramont (1919) reported boron in Fiskenæsset kornerupine, while Herd (1973) found none. New analyses (ion microprobe mass analyser and spectrophotometric) of kornerupine in three specimens from the type locality, including the specimens analysed by Lorenzen and Herd, indicate the presence of boron in all three, in amounts ranging from 0.50 to 1.44 wt.% B203, e.g. (Li0.04 Na0.01 Ca0.01) (Mg3.49 Mn0.01 Fe0.17 Ti0.01 Al5.64)Σ9.30 (Si3.67 Al1.02 B0.31)Σ5 O21 (OH0.99 F0.01) for Lorenzen's specimen. Textures and chemical compositions suggest that kornerupine crystallized in equilibrium in the following assemblages, all with anorthite (An 92-95) and phlogopite (XFe = atomic Fe/(Fe + Mg) = 0.028-0.035): (1) kornerupine (0.045)-gedrite (0.067); (2) kornerupine (0.038-0.050)-sapphirine (0.032-0.035); and (3) kornerupine (0.050)-hornblende. Fluorine contents of kornerupine range from 0.01 to 0.06%, of phlogopite, from 0.09 to 0.10%. In the first assemblage, sapphirine (0.040) and corundum are enclosed in radiating bundles of kornerupine; additionally sapphirine, corundum, and/or gedrite occur with chlorite and pinite (cordierite?) as breakdown products of kornerupine. Kornerupine may have formed by reactions such as: gedrite + sapphirine + corundum + B203 (in solution) + H20 = kornerupine + anorthite + Na-phlogopite under conditions of the granulite facies. Boron for kornerupine formation was most likely remobilized by hydrous fluids from metasedimentary rocks occurring along the upper contact of the Fiskenæsset gabbro-anorthosite complex with amphibolite.

Prismatine: Revalidation for Boron-Rich Compositions in the Kornerupine Group

Kornerupine and prismatine were introduced independently by Lorenzen in 1884 (but published in 1886 and 1893) and by Sauer in 1886, respectively. Ussing (1889) showed that the two minerals were sufficiently close crystallographically and chemically to be regarded as one species. However, recent analyses of boron using the ion microprobe and crystal structure refinement, indicate that the boron content of one tetrahedral site in kornerupine ranges from 0 to 1. Kornerupine and prismatine, from their respective type localities of Fiskenaesset, Greenland and Waldheim, Germany, are distinct minerals, members of an isomorphic series differing in boron content. For this reason, we re-introduce Sauer's name prismatine for kornerupines with B > 0.5 atoms per formula unit (p.f.u.) of 22(O,OH,F), and restrict the name kornerupine sensu stricto to kornerupines with B < 0.5 p.f.u. Kornerupine sensu lato is an appropriate group name for kornerupine of unknown boron content. Kornerupine sensu stricto and prismatine from the type localities differ also in Fe2+/Mg ratio, Si - (Mg + Fe2+ + Mn) content, Al content, F content, colour, density, cell parameters, and paragenesis. Both minerals formed under granulite-facies conditions with sapphirine and phlogopite, but kornerupine sensu stricto is associated with anorthite and homblende or gedrite, whereas prismatine is found with oligoclase (An9-13), sillimanite, garnet, and/or tourmaline. Occurrences at other localities suggest that increasing boron content extends the stability range of prismatine relative to that of kornerupine sensu stricto.

Petrology and geochemistry at DSDP Site 89-462

The 16 samples of Deep Sea Drilling Project (DSDP) Leg 89 basalts that we analyzed for whole rock major and trace elements and for mineralogic compositions are identical to some of the basalts recovered during Leg 61. Leg 89 samples are mostly olivine-plagioclase-clinopyroxene sparsely phyric basalts and exhibit a wide variety of textures. These basalts have lower TiO2 at a given Mg/(Mg+Fe2+)*100 than MORB (midocean ridge basalt). We recognize three major chemical types of basalts in the Nauru Basin. We believe that different degrees of partial melting, modified by fractional crystallization and possibly by magma mixing at shallow depths, can explain the chemical differences among the three groups. This petrogenetic model is consistent with the observed downhole chemical-chronostratigraphic relations of the samples. New 87Sr/86Sr and U3Nd/144Nd analyses of basalt samples from DSDP Site 462 indicate that the Nauru Basin igneous complex is within the Sr-Nd isotopic range of ocean island basalt. Thus the Nauru Basin igneous complex resembles MORB in many aspects of its chemistry, morphology, and secondary alteration patterns (Larson, Schlanger, et al., 1981), but not in its isotopic characteristics. If it were not for the unambiguous evidence that the Nauru Basin complex was erupted off-ridge, the complex could easily be interpreted as normal oceanic layer 2. For this reason, we speculate that the Nauru Basin igneous complex was produced in an oceanic riftlike environment when multiple, fast-propagating rifts were formed during the fast seafloor spreading episode in the Cretaceous.

Geochemistry of minerals in CRP sediment cores from the Ross Sea, Antarctica

Sixty-four volcanic chists, sandstones and tephras between 5.95 and 618.19 meters below sea floor (mbsf) in the Cape Roberts Project cores 2 and 2A cores (CRP-2/2A) were examined for Cenozoic and Mesozoic volcanic components, using optical and Scanning Electron Microscopy. Minerals and glass shards in a selection of samples were analysed by electron microprobe fined with an EDAX detector. Laser-Ablation ICP-Mass-Spectrometry (ICP-MS) was used to determine rare earth elements and 14 additional trace elements in glass shards, pyroxenes and feldspars in order to pin-point the onset of McMurdo Volcanic Group (MVG) activity in the stratigraphic column. Pumices in tephra layers of peralkaline phonolite composition in Unit 7.2 -between 108 and 114 mbsf - were also analysed for trace elements by ICP-MS. This tephra unit is not reworked and its isotopic age (21.44 ± 0.05 Ma) is the age of deposition. The height of the eruptive column responsible for the deposition of the tephra was probably less than 8 km; the source was local, probably within 30 km from the drill site. Phonolite of unit 7.2 of CRP-2/2A has no direct petrogenetic relation with the peralkaline trachyte in the tephra-enriched layer of CRP-1 at 116.55 mbsf. Volcanic clasts and sand grains (glass shards, aegirine-augite, anorthoclase) related to Cenozoic activity of MVG were observed only starting from Unit 9.8, where they are dated at 24.22 ± 0.06 Ma at c. 280 mbsf. In this unit the lowest- occurring basaltic glass shard is found at 297.54 mbsf. Sampled McMurdo volcanics are generally vesicular and vary in composition from alkali basalt to trachyte and peralkaline phonolite. By contrast, below 320 mbsf, aphyric or slightly-porphyritic volcanic clasts become more abundant but they are all non-vesiculated, pigeconite and ilmenite-bearing basalts and dolerite of tholeiitic affinity. These rocks are considered to be related to lava flows and associated intrusions of Jurassic age (Kirkpatrick basalts and Ferrar dolerite). As in CRP-1, McMurdo volcanics appear to derive from a variety of lithologics. Besides glaciers, a dominant role of wind transportation from exposed volcanic rocks may be inferred from the contemporary occurrence of glass shards of different compositions at depths above 297.54 mbsf. These data confirm that the onset of magmatic activity in southern Victoria Land is considerably delayed (by about 24 Ma) with respect to northern Victoria Land.

Major and trace element analyses of ODP Leg 120 volcanic rocks (Table 1)

Widespread Lower Cretaceous magmatism occurred along the Indian-Australian/Antarctic margins, and in the juvenile Indian Ocean, during the rifting of eastern Gondwana. The formation of this magmatic province probably began around 120-130 Ma with the eruption of basalts on the Naturaliste Plateau and at Bunbury, western Australia. On the northeast margin of India, activity began around 117 Ma with the Rajmahal continental basalts and associated lamprophyre intrusions. The formation of the Kerguelen Plateau in the Indian Ocean began no later than 114 Ma. Ultramafic lamprophyres (alnoites) were emplaced in the Prince Charles Mountains near the Antarctic continental margin at ~ 110 Ma. These events are considered to be related to a major mantle plume, the remnant of which is situated beneath the region of Kerguelen and Heard islands at the present day. Geochemical data are presented for each of these volcanic suites and are indicative of complex interactions between asthenosphere-derived magmas and the continental lithosphere. Kerguelen Plateau basalts have Sr and Nd isotopic compositions lying outside the field for Indian Ocean mid-ocean ridge basalts (MORB) but, with the exception of Site 738 at the southern end of the plateau, within the range of more recent hotspot basalts from Kerguelen and Heard Islands. However, a number of the plateau tholeiites are characterized by lower 206Pb/204Pb ratios than are basalts from Kerguelen Island, and many also have anomalously high La/Nb ratios. These features suggest that the source of the Kerguelen Plateau basalts suffered contamination by components derived from the Gondwana continental lithosphere. An extreme expression of this lithospheric signature is shown by a tholeiite from Site 738, suggesting that the southernmost part of the Kerguelen Plateau may be underlain by continental crust. The Rajmahal tholeiites mostly fall into two distinct geochemical groups. Some Group I tholeiites have Sr and Nd isotopic compositions and incompatible element abundances, similar to Kerguelen Plateau tholeiites from Sites 749 and 750, indicating that the Kerguelen-Heard mantle plume may have directly furnished Rajmahal volcanism. However, their elevated 207Pb/204Pb ratios indicate that these magmas did not totally escape contamination by continental lithosphere. In contrast to the Group I tholeiites, significant contamination is suggested for Group II Rajmahal tholeiites, on the basis of incompatible element abundances and isotopic compositions. The Naturaliste Plateau and the Bunbury Basalt samples show varying degrees of enrichment in incompatible elements over normal MORB. The Naturaliste Plateau samples (and Bunbury Basalt) have high La/Nb ratios, a feature not inconsistent with the notion that the plateau may consist of stretched continental lithosphere, near the ocean-continent divide.