Variance Reduction in Analytical Chemistry : New Numerical Methods in Chemometrics and Molecular Simulation

This thesis is based on five papers addressing variance reduction in different ways. The papers have in common that they all present new numerical methods. Paper I investigates quantitative structure-retention relationships from an image processing perspective, using an artificial neural network to preprocess three-dimensional structural descriptions of the studied steroid molecules. Paper II presents a new method for computing free energies. Free energy is the quantity that determines chemical equilibria and partition coefficients. The proposed method may be used for estimating, e.g., chromatographic retention without performing experiments. Two papers (III and IV) deal with correcting deviations from bilinearity by so-called peak alignment. Bilinearity is a theoretical assumption about the distribution of instrumental data that is often violated by measured data. Deviations from bilinearity lead to increased variance, both in the data and in inferences from the data, unless invariance to the deviations is built into the model, e.g., by the use of the method proposed in paper III and extended in paper IV. Paper V addresses a generic problem in classification; namely, how to measure the goodness of different data representations, so that the best classifier may be constructed. Variance reduction is one of the pillars on which analytical chemistry rests. This thesis considers two aspects on variance reduction: before and after experiments are performed. Before experimenting, theoretical predictions of experimental outcomes may be used to direct which experiments to perform, and how to perform them (papers I and II). After experiments are performed, the variance of inferences from the measured data are affected by the method of data analysis (papers III-V).

Constraints on mantle melting from major and trace element systematics in residual abyssal peridotites

ZusammenfassungDie Bildung von mittelozeanischen Rückenbasalten (MORB) ist einer der wichtigsten Stoffflüsse der Erde. Jährlich wird entlang der 75.000 km langen mittelozeanischen Rücken mehr als 20 km3 neue magmatische Kruste gebildet, das sind etwa 90 Prozent der globalen Magmenproduktion. Obwohl ozeanische Rücken und MORB zu den am meisten untersuchten geologischen Themenbereichen gehören, existieren weiterhin einige Streit-fragen. Zu den wichtigsten zählt die Rolle von geodynamischen Rahmenbedingungen, wie etwa Divergenzrate oder die Nähe zu Hotspots oder Transformstörungen, sowie der absolute Aufschmelzgrad, oder die Tiefe, in der die Aufschmelzung unter den Rücken beginnt. Diese Dissertation widmet sich diesen Themen auf der Basis von Haupt- und Spurenelementzusammensetzungen in Mineralen ozeanischer Mantelgesteine.Geochemische Charakteristika von MORB deuten darauf hin, dass der ozeanische Mantel im Stabilitätsfeld von Granatperidotit zu schmelzen beginnt. Neuere Experimente zeigen jedoch, dass die schweren Seltenerdelemente (SEE) kompatibel im Klinopyroxen (Cpx) sind. Aufgrund dieser granatähnlichen Eigenschaft von Cpx wird Granat nicht mehr zur Erklärung der MORB Daten benötigt, wodurch sich der Beginn der Aufschmelzung zu geringeren Drucken verschiebt. Aus diesem Grund ist es wichtig zu überprüfen, ob diese Hypothese mit Daten von abyssalen Peridotiten in Einklang zu bringen ist. Diese am Ozeanboden aufgeschlossenen Mantelfragmente stellen die Residuen des Aufschmelz-prozesses dar, und ihr Mineralchemismus enthält Information über die Bildungs-bedingungen der Magmen. Haupt- und Spurenelementzusammensetzungen von Peridotit-proben des Zentralindischen Rückens (CIR) wurden mit Mikrosonde und Ionensonde bestimmt, und mit veröffentlichten Daten verglichen. Cpx der CIR Peridotite weisen niedrige Verhältnisse von mittleren zu schweren SEE und hohe absolute Konzentrationen der schweren SEE auf. Aufschmelzmodelle eines Spinellperidotits unter Anwendung von üblichen, inkompatiblen Verteilungskoeffizienten (Kd's) können die gemessenen Fraktionierungen von mittleren zu schweren SEE nicht reproduzieren. Die Anwendung der neuen Kd's, die kompatibles Verhalten der schweren SEE im Cpx vorhersagen, ergibt zwar bessere Resultate, kann jedoch nicht die am stärksten fraktionierten Proben erklären. Darüber hinaus werden sehr hohe Aufschmelzgrade benötigt, was nicht mit Hauptelementdaten in Einklang zu bringen ist. Niedrige (~3-5%) Aufschmelzgrade im Stabilitätsfeld von Granatperidotit, gefolgt von weiterer Aufschmelzung von Spinellperidotit kann jedoch die Beobachtungen weitgehend erklären. Aus diesem Grund muss Granat weiterhin als wichtige Phase bei der Genese von MORB betrachtet werden (Kapitel 1).Eine weitere Hürde zum quantitativen Verständnis von Aufschmelzprozessen unter mittelozeanischen Rücken ist die fehlende Korrelation zwischen Haupt- und Spuren-elementen in residuellen abyssalen Peridotiten. Das Cr/(Cr+Al) Verhältnis (Cr#) in Spinell wird im Allgemeinen als guter qualitativer Indikator für den Aufschmelzgrad betrachtet. Die Mineralchemie der CIR Peridotite und publizierte Daten von anderen abyssalen Peridotiten zeigen, dass die schweren SEE sehr gut (r2 ~ 0.9) mit Cr# der koexistierenden Spinelle korreliert. Die Auswertung dieser Korrelation ergibt einen quantitativen Aufschmelz-indikator für Residuen, welcher auf dem Spinellchemismus basiert. Damit kann der Schmelzgrad als Funktion von Cr# in Spinell ausgedrückt werden: F = 0.10×ln(Cr#) + 0.24 (Hellebrand et al., Nature, in review; Kapitel 2). Die Anwendung dieses Indikators auf Mantelproben, für die keine Ionensondendaten verfügbar sind, ermöglicht es, geochemische und geophysikalischen Daten zu verbinden. Aus geodynamischer Perspektive ist der Gakkel Rücken im Arktischen Ozean von großer Bedeutung für das Verständnis von Aufschmelzprozessen, da er weltweit die niedrigste Divergenzrate aufweist und große Transformstörungen fehlen. Publizierte Basaltdaten deuten auf einen extrem niedrigen Aufschmelzgrad hin, was mit globalen Korrelationen im Einklang steht. Stark alterierte Mantelperidotite einer Lokalität entlang des kaum beprobten Gakkel Rückens wurden deshalb auf Primärminerale untersucht. Nur in einer Probe sind oxidierte Spinellpseudomorphosen mit Spuren primärer Spinelle erhalten geblieben. Ihre Cr# ist signifikant höher als die einiger Peridotite von schneller divergierenden Rücken und ihr Schmelzgrad ist damit höher als aufgrund der Basaltzusammensetzungen vermutet. Der unter Anwendung des oben erwähnten Indikators ermittelte Schmelzgrad ermöglicht die Berechnung der Krustenmächtigkeit am Gakkel Rücken. Diese ist wesentlich größer als die aus Schweredaten ermittelte Mächtigkeit, oder die aus der globalen Korrelation zwischen Divergenzrate und mittels Seismik erhaltene Krustendicke. Dieses unerwartete Ergebnis kann möglicherweise auf kompositionelle Heterogenitäten bei niedrigen Schmelzgraden, oder auf eine insgesamt größere Verarmung des Mantels unter dem Gakkel Rücken zurückgeführt werden (Hellebrand et al., Chem.Geol., in review; Kapitel 3).Zusätzliche Informationen zur Modellierung und Analytik sind im Anhang A-C aufgeführt

Partial melting within the lower crust – Constraints from bulk rock geochemical data and trace element mineral analyses of granulites and migmatites from Central Finland

This PhD thesis concerns geochemical constraints on recycling and partial melting of Archean continental crust. A natural example of such processes was found in the Iisalmi area of Central Finland. The rocks from this area are Middle to Late Archean in age and experienced metamorphism and partial melting between 2.7-2.63 Ga. The work is based on extensive field work. It is furthermore founded on bulk rock geochemical data as well as in-situ analyses of minerals. All geochemical data were obtained at the Institute of Geosciences, University of Mainz using X-ray fluorescence, solution ICP-MS and laser ablation-ICP-MS for bulk rock geochemical analyses. Mineral analyses were accomplished by electron microprobe and laser ablation ICP-MS. Fluid inclusions were studied by microscope on a heating-freezing-stage at the Geoscience Center, University Göttingen. Part I focuses on the development of a new analytical method for bulk rock trace element determination by laser ablation-ICP-MS using homogeneous glasses fused from rock powder on an Iridium strip heater. This method is applicable for mafic rock samples whose melts have low viscosities and homogenize quickly at temperatures of ~1200°C. Highly viscous melts of felsic samples prevent melting and homogenization at comparable temperatures. Fusion of felsic samples can be enabled by addition of MgO to the rock powder and adjustment of melting temperature and melting duration to the rock composition. Advantages of the fusion method are low detection limits compared to XRF analyses and avoidance of wet-chemical processing and use of strong acids as in solution ICP-MS as well as smaller sample volumes compared to the other methods. Part II of the thesis uses bulk rock geochemical data and results from fluid inclusion studies for discrimination of melting processes observed in different rock types. Fluid inclusion studies demonstrate a major change in fluid composition from CO2-dominated fluids in granulites to aqueous fluids in TTG gneisses and amphibolites. Partial melts were generated in the dry, CO2-rich environment by dehydration melting reactions of amphibole which in addition to tonalitic melts produced the anhydrous mineral assemblages of granulites (grt + cpx + pl ± amph or opx + cpx + pl + amph). Trace element modeling showed that mafic granulites are residues of 10-30 % melt extraction from amphibolitic precursor rocks. The maximum degree of melting in intermediate granulites was ~10 % as inferred from modal abundances of amphibole, clinopyroxene and orthopyroxene. Carbonic inclusions are absent in upper-amphibolite facies migmatites whereas aqueous inclusion with up to 20 wt% NaCl are abundant. This suggests that melting within TTG gneisses and amphibolites took place in the presence of an aqueous fluid phase that enabled melting at the wet solidus at temperatures of 700-750°C. The strong disruption of pre-metamorphic structures in some outcrops suggests that the maximum amount of melt in TTG gneisses was ~25 vol%. The presence of leucosomes in all rock types is taken as the principle evidence for melt formation. However, mineralogical appearance as well as major and trace element composition of many leucosomes imply that leucosomes seldom represent frozen in-situ melts. They are better considered as remnants of the melt channel network, e.g. ways on which melts escaped from the system. Part III of the thesis describes how analyses of minerals from a specific rock type (granulite) can be used to determine partition coefficients between different minerals and between minerals and melt suitable for lower crustal conditions. The trace element analyses by laser ablation-ICP-MS show coherent distribution among the principal mineral phases independent of rock composition. REE contents in amphibole are about 3 times higher than REE contents in clinopyroxene from the same sample. This consistency has to be taken into consideration in models of lower crustal melting where amphibole is replaced by clinopyroxene in the course of melting. A lack of equilibrium is observed between matrix clinopyroxene / amphibole and garnet porphyroblasts which suggests a late stage growth of garnet and slow diffusion and equilibration of the REE during metamorphism. The data provide a first set of distribution coefficients of the transition metals (Sc, V, Cr, Ni) in the lower crust. In addition, analyses of ilmenite and apatite demonstrate the strong influence of accessory phases on trace element distribution. Apatite contains high amounts of REE and Sr while ilmenite incorporates about 20-30 times higher amounts of Nb and Ta than amphibole. Furthermore, trace element mineral analyses provide evidence for magmatic processes such as melt depletion, melt segregation, accumulation and fractionation as well as metasomatism having operated in this high-grade anatectic area.

Protein and lipid binding parameters in rainbow trout (Oncorhynchus mykiss) blood and liver fractions to extrapolate from an in vitro metabolic degradation assay to in vivo bioaccumulation potential of hydrophobic organic chemicals

Binding of hydrophobic chemicals to colloids such as proteins or lipids is difficult to measure using classical microdialysis methods due to low aqueous concentrations, adsorption to dialysis membranes and test vessels, and slow kinetics of equilibration. Here, we employed a three-phase partitioning system where silicone (polydimethylsiloxane, PDMS) serves as a third phase to determine partitioning between water and colloids and acts at the same time as a dosing device for hydrophobic chemicals. The applicability of this method was demonstrated with bovine serum albumin (BSA). Measured binding constants (K(BSAw)) for chlorpyrifos, methoxychlor, nonylphenol, and pyrene were in good agreement with an established quantitative structure-activity relationship (QSAR). A fifth compound, fluoxypyr-methyl-heptyl ester, was excluded from the analysis because of apparent abiotic degradation. The PDMS depletion method was then used to determine partition coefficients for test chemicals in rainbow trout (Oncorhynchus mykiss) liver S9 fractions (K(S9w)) and blood plasma (K(bloodw)). Measured K(S9w) and K(bloodw) values were consistent with predictions obtained using a mass-balance model that employs the octanol-water partition coefficient (K(ow)) as a surrogate for lipid partitioning and K(BSAw) to represent protein binding. For each compound, K(bloodw) was substantially greater than K(S9w), primarily because blood contains more lipid than liver S9 fractions (1.84% of wet weight vs 0.051%). Measured liver S9 and blood plasma binding parameters were subsequently implemented in an in vitro to in vivo extrapolation model to link the in vitro liver S9 metabolic degradation assay to in vivo metabolism in fish. Apparent volumes of distribution (V(d)) calculated from the experimental data were similar to literature estimates. However, the calculated binding ratios (f(u)) used to relate in vitro metabolic clearance to clearance by the intact liver were 10 to 100 times lower than values used in previous modeling efforts. Bioconcentration factors (BCF) predicted using the experimental binding data were substantially higher than the predicted values obtained in earlier studies and correlated poorly with measured BCF values in fish. One possible explanation for this finding is that chemicals bound to proteins can desorb rapidly and thus contribute to metabolic turnover of the chemicals. This hypothesis remains to be investigated in future studies, ideally with chemicals of higher hydrophobicity.

Magma differentiation fractionates Mo isotope ratios: Evidence from the Kos Plateau Tuff (Aegean Arc)

We investigated high temperature Mo isotope fractionation in a hydrous supra-subduction volcano-plutonic system (Kos, Aegean Arc, Greece) in order to address the debate on the δ98/95Mo variability of the continental crust. In this igneous system, where differentiation is interpreted to be dominated by fractional crystallization, bulk rock data from olivine basalt to dacite show δ98/95Mo ratios increasing from +0.3 to +0.6‰ along with Mo concentrations increasing from 0.8 to 4.1 μg g−1. Data for hornblende and biotite mineral separates reveal the extraction of light Mo into crystallizing silicates, with minimum partition coefficients between hornblende- silicate melt and biotite-silicate melt of 0.6 and 0.4 δ98/95Mo, respectively. Our data document significant Mo isotope fractionation at magmatic temperatures, hence, the igneous contribution to continental runoff is variable, besides probable source-related variability. Based on these results and published data an average continental δ98/95Mo of +0.3 to +0.4‰ can be derived. This signature corresponds more closely to the average of published data of dissolved Mo loads of large rivers than previous estimates and is consistent with an upper limit of δ98/95Mo = 0.4‰ of the Earth's upper crust as derived from the analysis of molybdenites.

The behaviour of incompatible elements during hydrous melting of metasomatized peridotite at 4–6 GPa and 1000 °C–1200 °C

Trace element behavior during hydrous melting of a metasomatized garnet–peridotite was examined at pressures of 4–6 GPa and temperatures of 1000 °C–1200 °C, conditions appropriate for fluid penetrating the mantle wedge atop the subducting slab. Experiments were performed in a rocking multi-anvil apparatus using a diamond-trap setup. The compositions of the fluid and melt phases were measured using the cryogenic LA-ICP-MS technique. The water-saturated solidus of the K-lherzolite composition is located between 900 °C and 1000 °C at 4 GPa and between 1000 °C and 1100 °C at 5 and 6 GPa. The partition coefficients between fluid or melt and clinopyroxene reveal an asymmetric MREE trough with a minimum at Dy. The clinopyroxene in equilibrium with aqueous fluids is characterized by DUfluid–cpx > DThfluid–cpx while DUmelt–cpx tends to be similar to DThmelt–cpx. The partition coefficients between fluid or melt and garnet reveal very strong light to heavy REE fractionation, DLa/DLu from 95 (hydrous melt) to 1600 (aqueous fluid). The LILE are highly incompatible with partition coefficients > 50. The behavior of HFSE are decoupled, with DZr,Hf close to 1 while DNb,Ta > 10. Garnet is characterized by DUmelt/fluid–garnet < DThmelt/fluid–garnet. A comparison of our experimental partitioning results for trivalent cations as well as the results from the literature and the calculations carried out using the lattice strain model adapted to the presence of water in the bulk system indicates that H2O in the fluid or melt phase has a prominent effect on trace element partitioning. Garnet in mantle rocks in equilibrium with an aqueous fluid is characterized by significantly higher Do(3 +) for REE in the X site of the garnet compared with the partitioning values of the optimal cation in garnet in equilibrium with hydrous melts. Our data show for the first time that the change in the nature of the mobile phase (fluid vs. melt) does affect the affinities of trace elements into the garnet crystal at conditions below the second critical endpoint of the system. The same also applies for clinopyroxene, although this is less clear. Consequently, our new data allow for refinements in predictive modeling of element transfer from the slab to the mantle wedge and of possible compositions of metasomatized mantle that sources OIB magmatism.

Variations of clinopyroxene/melt element partitioning during assimilation of olivine/peridotite by low-Mg diorite magma

The effects of crystal chemistry and melt composition on the control of clinopyroxene/melt element partitioning (D) during the assimilation of olivine/peridotite by felsic magma have been investigated in Mesozoic high-Mg diorites from North China. The assimilation resulted in significant increase of Mg, Cr and Ni and only slight (< 30%) decrease of incompatible elements of the magma, and the compositional variations have been mirrored by the normally and reversely zoned clinopyroxene microphenocrysts formed at the early stage of the magma evolution. The Mg# [100 × Mg / (Mg + Fe)] values of the reversely zoned clinopyroxenes increase from 65 to 75 in the core to 85–90 in the high-Mg midsection, and reduce back to 73–79 at the rim. Trace element profiles across all these clinopyroxene domains have been measured by LA-ICP-MS. The melt trace element composition has been constrained from bulk rock analyses of the fine-grained low- and high-Mg diorites. Clinopyroxene/melt partition coefficients for rare earth elements (REE) and Y in the high-Mg group zonings (Mg# > 73–79, DDy < 1.2) are positively correlated with tetrahedral IVAl and increase by a factor of 3–4 as tetrahedral IVAl increases from 0.01 to 0.1 per formula unit (pfu). These systematic variations are interpreted to be controlled by the clinopyroxene composition. In contrast, partition coefficients for low-Mg group zonings (Mg# < 75–79, DDy > 1.2) are elevated by up to an order of magnitude (for REE and Y) or more (for Zr and Hf) at similar IVAl, indicating dominant control of melt composition/structure. DZr and DHf show a larger sensitivity to the compositional change of crystal and melt than DREE. DTi values for the low- and high-Mg zonings show a uniform dependence on IVAl. DSr and DLi are insensitive to the compositional change of clinopyroxene and melt, resulting in Sr depletions in the clinopyroxene zonings with elevated REE without crystallization of plagioclase. Our observations show that crystal chemistry and melt composition/structure may alternatively control clinopyroxene/melt partitioning during the assimilation of peridotite by felsic magma, and may be useful for deciphering clinopyroxene compositions and related crust–mantle processes.

Apatite as an indicator of fluid salinity: An experimental study of chlorine and fluorine partitioning in subducted sediments

In order to constrain the salinity of subduction zone fluids, piston-cylinder experiments have been conducted to investigate the partitioning behaviour of Cl and F in subducted sediments. These experiments were performed at H2O-undersaturated conditions with a synthetic pelite starting composition containing 800 ppm Cl, over a pressure and temperature range of 2.5–4.5 GPa and 630–900 °C. Repetitive experiments were conducted with 1900 ppm Cl + 1000 ppm F, and 2100 ppm Cl. Apatite represents the most Cl-abundant mineral phase, with Cl concentration varying in the range 0.1–2.82 wt%. Affinity for Cl decreases over the following sequence: aqueous fluid > apatite ⩾ melt > other hydrous minerals (phengite, biotite and amphibole). It was found that addition of F to the Cl-bearing starting composition significantly lowers the Cl partition coefficients between apatite and melt (DClAp–melt) and apatite and aqueous fluid (DClAp–aq). Cl–OH exchange coefficients between apatite and melt (KdCl–OHAp–melt) and apatite and aqueous fluid (KdCl–OHAp–aq) were subsequently calculated. KdCl–OHAp–melt was found to vary from 1 to 58, showing an increase with temperature and a decrease with pressure and displaying a regular decrease with increasing H2O content in melt. Mole fractions of Cl and OH in melt were calculated based on an ideal mixing model for H2O, OH, O, Cl and F. The Cl contents of other hydrous minerals (phengite, biotite and amphibole) fall between 200 and 800 ppm, with resultant Cl partition coefficients from 0.02 to 0.49, appearing independent of the bulk Cl and F content. Preliminary data from this study show that the partitioning behaviour of F is strongly in favour of apatite relative to melt and phengite, with DFAp–melt = 15–51. Apatites from representative eclogite facies metasediments were examined and found to have low Cl contents close to ∼100 ppm. Calculations using our experimentally determined KdCl–OHAp–aq of 0.004 at 2.5 GPa, 630 °C indicate a low salinity character (0.5–2 wt% NaCleq) for the fluid formed during dehydration of subducted oceanic sediment at ∼80 km depth.

Geochemistry of ODP Leg 135 vitrophyric rhyolite samples

In-situ proton-microprobe analyses are presented for glasses, plagioclases, pyroxenes, olivines, and spinels in eleven samples from Sites 834-836, 839, and 841 (vitrophyric rhyolite), plus a Tongan dacite. Elements analyzed are Mn, Ni, Cu, Zn, Ga, Rb, Sr, Y, Zr, Pb, and Sn (in spinels only). The data are used to calculate two sets of partition coefficients, one set based on the ratio of element in mineral/element in coexisting glass. The second set of coefficients, thought to be more robust, is corrected by application of the Rayleigh fractionation equations, which requires additional use of modal data. Data are presented for phenocryst core-rim phases and microphenocryst-groundmass phases from a few samples. Comparison with published coefficients reveals an overall consistency with those presented here, but with some notable anomalies. Examples are relatively high Zr values for pyroxenes and abnormally low Mn values in olivines and clinopyroxenes from Site 839 lavas. Some anomalies may reflect kinetic effects, but interpretation of the coefficients is complicated, especially in olivines from Sites 836 and 839, by possible crystal-liquid disequilibrium resulting from mixing processes.

Geochemistry of Fe-Mn crusts in sediments of the Central Indian Ridge

Layered Fe-Mn crusts from the off-axis region of the first segment of the Central Indian Ridge north of the Rodrigues Triple Junction were studied geochemically and mineralogically. Vernadite (delta-MnO2) is the main mineral oxide phase. 230Thxs and Co concentrations suggest high growth rates of up to 29 mm/Myr and a maximum age of the basal crust layer of 1 Ma. Whereas most of the major and minor elements show concentrations which are typical of hydrogenetic formation, Co, Pb, Ni and Ti concentrations are strikingly lower. Concentrations and distribution of the strictly trivalent rare-earths and yttrium (REY) are typical of hydrogenetic ferromanganese oxide precipitates, but in marked contrast, the crusts are characterized by negative CeSN (shale normalized) anomalies and (Ce/Pr)SN ratios less than unity. Profiles through the crusts reveal only minor variations of the REY distribution and (Ce/Pr)SN ratios range from 0.45 to 0.68 (compared to ratios of up to 2 for typical hydrogenetic crusts from the Central Indian Basin). The apparent bulk partition coefficients between the crusts and seawater suggest that for the strictly trivalent REY the adsorption-desorption equilibrium has been reached. Positive Ce anomalies in the partition coefficient patterns reveal preferential uptake of Ce, but to a lesser extent than in normal hydrogenetic crusts. A new parameter (excess Ce, Cexs) to quantify the degree of decoupling of Ce from REY(III) is established on the basis of partition coefficients. Cexs/Cebulk ratios suggest that the CIR crusts formed by precipitation of Fe-Mn oxides from a hydrothermal plume and that in hydrothermal plumes and normal seawater the enrichment of Ce results from the same oxidative sorption process. The growth rates, calculated with 230Thxs data as well as with the Co formula, are inversely related to Cexs.

Isotopic and chemical compositions of rocks and minerals from the TAG hydrothermal mound, ODP Site 957

Strontium- and oxygen-isotopic measurements of samples recovered from the Trans-Atlantic Geotraverse (TAG) hydrothermal mound during Leg 158 of the Ocean Drilling Program provide important constraints on the nature of fluid-rock interactions during basalt alteration and mineralization within an active hydrothermal deposit. Fresh Mid-Ocean Ridge Basalt (MORB), with a 87Sr/86Sr of 0.7026, from the basement beneath the TAG mound was altered at both low and high temperatures by seawater and altered at high temperature by near end-member black smoker fluids. Pillow breccias occurring beneath the margins of the mound are locally recrystallized to chlorite by interaction with large volumes of conductively heated seawater (>200°C). The development of a silicified, sulfide-mineralized stockwork within the basaltic basement follows a simple paragenetic sequence of chloritization followed by mineralization and the development of a quartz+pyrite+paragonite stockwork cut by quartz-pyrite veins. Initial alteration involved the development of chloritic alteration halos around basalt clasts by reaction with a Mg-bearing mixture of upwelling, high-temperature (>300°C), black smoker-type fluid with a minor (<10%) proportion of seawater. Continued high-temperature (>300°C) interaction between the wallrock and these Mg-bearing fluids results in the complete recrystallization of the wallrock to chlorite+quartz+pyrite. The quartz+pyrite+paragonite assemblage replaces the chloritized basalts, and developed by reaction at 250-360°C with end-member hydrothermal fluids having 87Sr/86Sr ~0.7038, similar to present-day vent fluids. The uniformity of the 87Sr/86Sr ratios of hydrothermal assemblages throughout the mound and stockwork requires that the 87Sr/86Sr ratio of end-member hydrothermal fluids has remained relatively constant for a time period longer than that required to change the interior thermal structure and plumbing network of the mound and underlying stockwork. Precipitation of anhydrite in breccias and as late-stage veins throughout most of the mound and stockwork, down to at least 125 mbsf, records extensive entrainment of seawater into the hydrothermal deposit. 87Sr/86Sr ratios indicate that most of the anhydrite formed from ~2:1 mixture of seawater and black smoker fluids (65%±15% seawater). Oxygen-isotopic compositions imply that anhydrite precipitated at temperatures between 147°C and 270°C and require that seawater was conductively heated to between 100°C and 180°C before mixing and precipitation occurred. Anhydrite from the TAG mound has a Sr-Ca partition coefficient Kd ~0.60±0.28 (2 sigma). This value is in agreement with the range of experimentally determined partition coefficients (Kd ~0.27-0.73) and is similar to those calculated for anhydrite from active black smoker chimneys from 21°N on the East Pacific Rise. The d18O (for SO4) of TAG anhydrite brackets the value of seawater sulfate oxygen (~9.5?). Dissolution of anhydrite back into the oceans during episodes of hydrothermal quiescence provides a mechanism of buffering seawater sulfate oxygen to an isotopically light composition, in addition to the precipitation and dissolution of anhydrite within the oceanic basement during hydrothermal recharge at the mid-ocean ridges.

Geochemistry of clinopyroxene megacrysts, xenocrysts, and phenocrysts in diabase dikes from ODP Hole 140-504B

We report the major, rare earth, and other trace element compositions of clinopyroxenes from two Leg 140, Hole 504B diabase dikes. These pyroxenes reflect a complex history of crystal growth and magma evolution. The large ranges of composition found reflect incorporation of exotic phenocrysts into the melt, the early formation of crystal clots before dike intrusion during an undercooling event, and in-situ fractionation of melt during and following dike emplacement. Some of the pyroxenes occur in coarse two- and three-phase glomerocrysts, which may be ôprotogabbrosö representing early stages of melt crystallization in the lower crust. Large variations in trace element composition are found. These likely reflect heterogeneous nucleation and growth of plagioclase and pyroxene in the melt, as well as complex interface kinetics that may affect partition coefficients during rapid crystal growth expected during undercooling. This can explain the formation of irregular chemical sector zoning in some equant anhedral phenocrysts. Undercooling of magmas in the lower crust most likely reflects input of fresh hot melt into a stagnating melt-storage zone. Dikes intruded upward from an inflated melt-storage zone during such a cycle are likely to be larger than those intruded from the storage zone between such cycles, when it would be deflated, consistent with the greater overall thickness of the phyric dikes in the Leg 140 section of Hole 504B.

Sr/Ca ratios of Cenozoic benthic foraminifera

A Cenozoic multi-species record of benthic foraminiferal calcite Sr/Ca has been produced and is corrected for interspecific offsets (typically less than 0.3 mmol/mol) and for the linear relationship between decreasing benthic foraminiferal Sr/Ca and increasing water depth. The water depth correction, determined from Holocene, Late Glacial Maximum and Eocene paleowater-depth transects, is ~0.1 mmol/mol/km. The corrected Cenozoic benthic foraminiferal Sr/Ca record ranges from 1.2 to 2.0 mmol/mol, and has been interpreted in terms of long-term changes in seawater Sr/Ca, enabling issues related to higher-resolution variability in Sr/Ca to be ignored. We estimate that seawater Sr/Ca was ~1.5 times modern values in the late Cretaceous, but declined rapidly into the Paleogene. Following a minimum in the Eocene, seawater Sr/Ca increased gradually through to the present day with a minimum superimposed on this trend centered in the late Miocene. By assuming scenarios for changing seawater calcium concentration, and using published carbonate accumulation rate data combined with suitable values for Sr partition coefficients into carbonates, the seawater Sr/Ca record is used to estimate global average river Sr fluxes. These fluxes are used in conjunction with the seawater strontium isotope curve and estimates of hydrothermal activity/tectonic outgassing to calculate changes in global average river 87Sr/86Sr through the Cenozoic. The absolute magnitude of Sr fluxes and isotopic compositions calculated in this way are subject to relatively large uncertainties. Nevertheless, our results suggest that river Sr flux increased from 35 Ma to the present day (roughly two-fold) accompanied by an overall increase in 87Sr/86Sr (by ~0 to 0.001). Between 75 and 35 Ma, river 87Sr/86Sr also increased (by ~0.001 to 0.002) but was accompanied by a decrease (two- to three-fold) in river Sr flux.

Chemical profiles in sediment and basalt at DSDP Hole 74-525A

Forty sediment and four basement basalt samples from DSDP Hole 525A, Leg 74, as well as nine basalt samples from southern and offshore Brazil, were subjected to instrumental neutron activation analysis. Thirty-two major, minor, and trace elements were determined. The downcore element concentration profiles and regression analyses show that the rare earth elements (REE) are present in significant amounts in both the carbonate and noncarbonate phases in sediments; Sr is concentrated in the carbonate phase, and most of the other elements determined exist mainly in the noncarbonate phase. The calculated partition coefficients of the REE between the carbonate phase and the free ion concentrations in seawater are high and increase with decreasing REE ionic radii from 3.9 x 10**6 for La to 15 x 10**6 for Lu. Calculations show that the lanthanide concentrations in South Atlantic seawater have not been changed significantly over the past 70 Ma. The Ce anomaly observed in the carbonate phase is a redox indicator of ancient seawater. Study of the Ce anomaly reveals that seawater was anoxic over the Walvis Ridge during the late Campanian. As the gap between South America and West Africa widened and the Walvis Ridge subsided from late Campanian to late Paleocene times, the water circulation of the South Atlantic improved and achieved oxidation conditions about 54 Ma that are similar to present seawater redox conditions in the world oceans. The chemical compositions of the basement rocks correspond to alkalic basalts, not mid-ocean ridge basalts (MORBs). The results add more evidence to support the hypothesis that the Walvis Ridge was formed by a series of volcanos moving over a "hot spot" near the Mid-Atlantic Ridge. From the chemical composition and REE pattern, one 112 Ma old basalt on the Brazilian continental shelf has been identified as an early stage MORB. To date, this is the oldest oceanic tholeiite recovered from the South Atlantic. This direct evidence indicates that the continental split between South America and Africa commenced > 112 Ma.

Geochemistry of sediments beneath the Tonga-Kermadec arc

The fate of subducted sediment and the extent to which it is dehydrated and/or melted before incorporation into arc lavas has profound implications for the thermo-mechanical nature of the mantle wedge and models for crustal evolution. In order to address these issues, we have undertaken the first measurements of 10Be and light elements in lavas from the Tonga-Kermadec arc and the sediment profile at DSDP site 204 outboard of the trench. The 10Be/9Be ratios in the Tonga lavas are lower than predicted from flux models but can be explained if (a) previously estimated sediment contributions are too high by a factor of 2-10, (b) the top 1-22 m of the incoming sediment is accreted, (c) large amounts of sediment erosion are proposed, or (d) the sediment component takes several Myr longer than the subducting plate to reach the magma source region beneath Tonga. The lavas form negative Th/Be-Li/Be arrays that extend from a depleted mantle source composition to lower Th/Be and Li/Be ratios than that of the bulk sediment. Thus, these arrays are not easily explained by bulk sediment addition and, using partition coefficients derived from experiments on the in-coming sediment, we show that they are also unlikely to result from fluid released during dehydration of the sediment (or altered oceanic crust). However, partial melts of the dehydrated sediment residue formed at ~800 °C during the breakdown of amphibole +/- plagioclase and in the absence of cordierite have significantly lowered Th/Be ratios. The lava arrays can be successfully modelled as 10-15% partial melts of depleted mantle after it has been enriched by the addition of 0.2-2% of these partial melts. Phase relations suggest that this requires that the top of the subducting crust reaches temperatures of ~800 °C by the time it attains ~ 80 km depth which is in excellent agreement with the results of recent numerical models incorporating a temperature-dependent mantle viscosity. Under these conditions the wet basalt solidus is also crossed yet there is no recognisable eclogitic signal in the lavas suggesting that on-going dehydration or strong thermal gradients in the upper part of the subducting plate inhibit partialmelting of the altered oceanic crust.