Cementation of kerogen-rich marls by alkaline fluids released during weathering of thermally metamorphosed marly sediments. Part I: Isotopic (C,O) study of the Khushaym Matruk natural analogue (central Jordan)

The Khushaym Matruk site in central Jordan may represent a natural analogue depicting the interaction of alkaline solutions with a clayey sedimentary formation or with clay-rich confining barriers at the interface with concrete structures in waste disposal sites. In this locality, past spontaneous combustion of organic matter in a clayey biomicritic formation produced a ca. 60 m-thick layer of cement-marble containing some of the high-temperature phases usually found in industrial cements (e.g., spurrite, brucite, and Ca-aluminate). A vertical cross-section of the underlying sediments was used in order to study the interaction between cement-marbles and neighbouring clayey limestones under weathering conditions. A thermodynamic approach of the alteration parageneses (calcite–jennite–afwillite–brucite and CSH phases) in the cement-marbles constrains the interacting solutions to have had pH-values between 10.5 and 12. Over 3 m, the sediments located beneath the metamorphic unit were compacted and underwent carbonation. They display large C and O isotopic variations with respect to “pristine” sediments from the bottom of the section. Low δ¹³C-values down to −31.4‰/PDB show the contribution of CO₂ derived from the oxidization of organic matter and from the atmosphere to the intense carbonation process affecting that particular sedimentary level. The size of the C isotopic anomalies, their geometrical extent and their coincidence with the variations of other markers like the Zn content, the structure of organic matter, the mineralogical composition, all argue that the carbonation process was induced by the percolation of high pH solutions which derived from the alteration of cement-marbles. The temperature of the carbonation process remains conjectural and some post-formation O isotopic reequilibration likely affected the newly-formed carbonate. Carbonation induced a considerable porosity reduction, both in fractures and matrixes. The Khushaym Matruk site may have some bearing to the early life of a repository site, when water saturation of the geological formations hosting the concrete structures is incomplete, enabling simultaneous diffusion of alkaline waters and gaseous CO₂ in the near field.     

Methyl mercury distributions in relation to the presence of nano- and picophytoplankton in an oceanic water column (Ligurian Sea, North-western Mediterranean)

Recent findings on the distribution of methylated mercury (MeHg_T) in waters have highlighted the importance of organic carbon remineralization on the production of these compounds in the open ocean. Here, we present the first time-series (20 monthly samplings between July 2007 and May 2009) of high-resolution vertical profiles (10-12 depths in a 2350 m water column) of MeHg_T distributions in an open ocean environment, the Ligurian Sea (North-western Mediterranean Sea). Concentrations varied within the sub-picomolar range (general mean: 0.30 ± 0.17 pmol L⁻¹, n = 214) with the lowest values at the surface, increasing with depth up to the oxygen minimum zone, and decreasing slowly at greater depth. Concentrations in the surface waters never exceeded 0.15 pmol L⁻¹, while the highest concentrations (up to 0.82 pmol L⁻¹) were associated to the hypoxycline during the autumn bloom. A detailed vertical MeHg_T profile reveals a “double-peak” pattern, coincidental with the two microbial layers described by Tanaka and Rassoulzadegan (2002), the so-called “microbial food web” in the euphotic zone (<100 m) and the “microbial loop” in the aphotic zone (>100 m). Temporal variations in the MeHg_T abundance and distribution in the water column were linked to seasonality. The highest MeHg_T concentrations were found in the oxygen minimum zone during the period of stratification, and coincide with the greatest abundance of nano- and picophytoplankton (cyanobacteria, nanoflagellates, etc.) in the euphotic layer. None of our deep MeHg_T measurements (∼100 m above the sea bottom) revealed a significant sedimentary source of MeHg_T. We explored the correlation between MeHg_T concentrations and the apparent oxygen utilization, a proxy of organic matter remineralization, over the study period. Results of this study strengthen the hypothesis that net mercury methylation in the open ocean occurs in the water column, is linked to organic matter regeneration, and is promoted by the presence of small-sized nano- and picophytoplankton, that dominate under oligotrophic conditions.     

An absolute paleotemperature record from 10 to 6 Ka inferred from fluid inclusion D/H ratios of a stalagmite from Vancouver Island, British Columbia, Canada

By using continuous helium flow during the crushing of calcite speleothem samples, we are able to recover liberated inclusion waters without isotopic fractionation. A paleotemperature record for the Jacklah Jill Cave locality, Vancouver Island, BC, was obtained from a 30-cm tall stalagmite that grew 10.3-6.3 Ka ago, using δ¹⁸O values of the crushed calcite and of the inclusion water as inferred from its δD. It is found that the locality experienced mean annual temperature variations up to 11 °C over a 4-Ka period in the early Holocene. At the beginning of the period, local temperature quickly increased from a minimum of ∼1 °C to around 10 °C, but this early climate optimum, about 3 °C warmer than today, only lasted for ∼1200 years. About 8.6 Ka ago, temperature had declined to ∼7 °C, approximately the same as the modern cave temperature. Since then, the study area has experienced only minor temperature fluctuations, but there was a brief fall to ∼4 °C at around 7 Ka ago, which might be caused by a short lived expansion of local alpine glaciers. The long-term T-dependence of δD was 1.47‰/°C, identical to the value in modern precipitation.     

Hydrous mineralogy of CM and CI chondrites from infrared spectroscopy and their relationship with low albedo asteroids

IR spectroscopy is one of the few techniques that can directly probe water molecules in rocks. This method has been used to characterize the mineralogy of hydrated/hydrous carbonaceous chondrites, and to link known meteorite families with spectroscopic observations of low albedo asteroids. In this paper, we present measurements of the infrared transmission spectra of matrix chunks from 3 CI and 9 CM chondrites. Spectra were measured at ambient conditions and then at different temperatures along a dehydration path toward high-T (∼300 °C) under primary vacuum. At ambient conditions, the 3-μm spectral range is always dominated by adsorbed atmospheric water molecules. Upon moderate (∼100 °C) and high (∼300 °C) heating under low pressure (P < 10⁻⁴ mbar), adsorbed water and then phyllosilicates interlayer water are removed, revealing a residual absorption band around 3 μm. This band is a characteristic IR feature of the phyllosilicate phases which dominate the mineralogical assemblage of hydrated carbonaceous chondrites. Among the CM chondrites, the high-T spectra reveal a strong variability that appears correlated with the alteration classification scheme of Rubin et al. (2007) and Howard et al. (2009a). The 3-μm band continuously evolves from a broad feature peaking at 3550-3600 cm⁻¹ for the weakly altered CMs (Murchison-type) to a sharp asymmetric peak at ∼3675 cm⁻¹ for the more extensively altered samples (Cold Bokkeveld-type). We attribute this spectral evolution to variations in the chemistry of the phyllosilicate phases from Fe-rich to Mg-rich. On the other hand, the 10-μm spectral region shows a single broad peak which does not compare with known terrestrial serpentine spectra, probably due to high structural disorder of the chondrite phyllosilicate phases. The present work clearly shows that previously published reflectance spectra of chondrites are biased by the presence of adsorbed terrestrial water molecules. Laboratory data collected under dry conditions are needed to reinterpret the chondrite-asteroid connection from the comparison of their 3-μm absorption features.     

Cr diffusion in orthopyroxene: Experimental determination, Mn-Cr thermochronology, and planetary applications

We have determined Cr diffusion coefficients (D) in orthopyroxene parallel to the a-, b-, and c-axial directions as a function temperature at f(O₂) corresponding to those of the wüstite-iron (WI) buffer. Diffusion is found to be significantly anisotropic with D(//c) > D(//b) > D(//a), conforming to an earlier theoretical prediction. Increase of f(O₂) from WI buffer conditions to 4.5 log unit above the buffer at 950 and 1050 °C leads to decrease of D(Cr) by a factor of two to three, possibly suggesting significant contribution from an interstitial diffusion mechanism. We have used the diffusion data to calculate the closure temperatures (T_c) of the Mn-Cr decay system in orthopyroxene as a function of initial temperature (T₀), grain size (a) and cooling rate for spherical and plane sheet geometries. We also present graphical relations that permit retrieval of cooling rates from knowledge of the resetting of Mn-Cr ages in orthopyroxene during cooling, T₀ and a. Application of these relations to the Mn-Cr age data of the cumulate eucrite Serra de Magé yields a T_c of 830-980 °C, and cooling rates of 2-27 °C/Myr at T_c and ∼1-13 °C/Myr at 500 °C. It is shown that the cooling of Serra de Magé to the closure temperature of the Mn-Cr system took place at its original site in the parent body, and thus implies a thickness for the eucrite crust in the commonly accepted HED parent body, Vesta, of greater than 30 km. This thickness of the eucrite crust is compatible only with a model of relatively olivine-poor bulk mineralogy in which olivine constitutes 19.7% of the total asteroidal mass.     

Determination of sulfide ligands and association with natural organic matter

Group B metals, such as Hg, Cu, Ag, Pb and Cd bind strongly to reduced inorganic and organic S(II−) ligands. These S(II−) ligands, stable in oxic waters for significant periods of time, occur at the <1–100 s nM concentrations. It is hypothesized that S(II−) ligands are stabilized as Cu–S molecules associated with organic matter by multi-ligand binding or in nano-pore encapsulations in organic matter. S(II−) ligands are estimated by two methods: purge/trap analysis as Cr-reducible sulfide (CRS), and strong ligand (SL_T) from a competitive ligand titration with Ag(I). The CRS/SL_T ratio is nearly one for selected samples. CRS correlates reasonably well (r² ∼ 0.5) with organic C with a slope of 14.6 nM per mg C. The conditional binding constant of Ag–SL is 11.3 for effluent associated with waste-water and decreases for river waters from about 12–8.8 as the strong sites are occupied with Ag(I).     

Determining accurate temperature-time paths from U-Pb thermochronology: An example from the Kaapvaal craton, southern Africa

Thermochronology has revolutionized our understanding of the establishment and evolution of lithospheric thermal structure. However, many potential benefits provided by the application of diffusion theory to thermochronology have yet to be fully exploited. This study uses apatite (T_c = 450-550 °C) and titanite (T_c = 550-650 °C) U-Pb ID-TIMS thermochronology at the single- to sub-grain scale to separate the variable effects of volume diffusion of Pb from metamorphic (over)growth above and below the T_c of a mineral. Data are presented from two ca. 3227 Ma tonalite samples from north and south of the Barberton Greenstone Belt (BGB), southern Africa. Two distinct populations of apatite from a sample north of the BGB record fast cooling followed by metamorphic growth ∼10 Myr later. Both apatite and titanite dates from south of the BGB show a strong correlation with the grain size and record 100 Myr of post-emplacement cooling. Complex core-rim zoning observed in cathodoluminescence images of apatite is interpreted to reflect metamorphic overgrowth above the T_c. The age and topology of grain size versus date curves from titanite and apatite are used in combination with a finite-difference numerical model to show that slow, non-linear, cooling and not thermal resetting is responsible for the observed distribution. The thermal histories from either side of the BGB are very different and provide unique insight into the BGB’s tectonic evolution: a ∼70 Myr period of apparent stability after ca. 3.2 Ga terrane assembly was followed by fast exhumation south of the BGB that led to lower-crustal melting and intrusion of granitic batholiths ca. 3.14-3.10 Ga.     

How depositional conditions control input, composition, and degradation of organic matter in sediments from the Chilean coastal upwelling region

In a comprehensive study, we compared depositional conditions, organic matter (OM) composition, and organic carbon turnover in sediments from two different depositional systems along the Chilean continental margin: at ∼23° S off Antofagasta and at ∼36° S off Concepción. Both sites lie within the Chilean coastal upwelling system and have an extended oxygen minimum zone in the water column. However, the northern site (23° S) borders the Atacama Desert, while the southern site (36° S) has a humid hinterland. Eight surface sediment cores (up to 30 cm long) from water depths of 126-1350 m were investigated for excess ²¹⁰Pb (²¹⁰Pb_xs) activity, total organic and total inorganic carbon concentrations (TOC and TIC, respectively), C/N-ratios, organic carbon isotopic compositions (δ¹³C), chlorin concentrations, Chlorin Indices (CI), and sulfate reduction rates (SRR). Sediment accumulation rates obtained from ²¹⁰Pb-analysis were similar in both regions (0.04-0.15 cm yr⁻¹ at 23° S, 0.10-0.19 cm yr⁻¹ at 36° S), although total ²¹⁰Pb_xs fluxes indicated that the vertical particle flux was higher at 36° S than at 23° S. We propose that sediment focusing in isolated deposition centers led to high sediment accumulation rates at 23° S. Furthermore, there were no indications for sediment mixing at 23° S, while bioturbation was intense at 36° S. δ¹³C-values (−24.5‰ to −20.1‰ vs. VPDB) and C/N-ratios (molar, 8.6-12.8) were characteristic of a predominantly marine origin of the sedimentary OM in both investigated areas. The extent of OM alteration in the water column was partly reflected in the surface sediments as chlorin concentrations decreased and C/N-ratios and CI increased with increasing water depth of the sampling site. SRR were lower at 23° S (areal SRR 0.12-0.60 mmol m⁻² d⁻¹) than at 36° S (areal SRR 0.82-1.18 mmol m⁻² d⁻¹), which was partly due to the greater water depth of most of the sediments investigated in the northern region and consistent with a lower quality of the sedimentary OM at 23° S. Reaction rate constants for TOC degradation that were obtained from measured SRR (k_SRR; 0.0004-0.0022 yr⁻¹) showed a good correspondence to k_TOC that were derived from the depth profiles of TOC (0.0003-0.0014 yr⁻¹). Both, k_SRR and k_TOC, reflect differences in OM composition. At 36° S they were related to the degradation state of bulk OM (represented by C/N-ratios), whereas near 23° S they were related to the freshness of a small fraction of labile OM (represented by CI). Our study shows that although rates of organic carbon accumulation were similar in both investigated sites, the extent and kinetics of organic carbon degradation were closely linked to differing depositional conditions.     

The origin of the fumaroles of La Solfatara (Campi Flegrei, South Italy)

The analysis of gaseous compositions from Solfatara (Campi Flegrei, South Italy) fumaroles since the early 1980s, clearly reveals a double thermobarometric signature. A first signature at temperatures of about 360 °C was inferred by methane-based chemical-isotopic geoindicators and by the H₂/Ar geothermometer. These high temperatures, close to the critical point of water, are representative of a deep zone where magmatic gases flash the hydrothermal liquid, forming a gas plume. A second signature was found to be at around 200-240 °C. At these temperatures, the kinetically fast reactive species (H₂ and CO) re-equilibrate in a pure vapor phase during the rise of the plume. A combination of these observations with an original interpretation of the oxygen isotopic composition of the two dominant species, i.e. H₂O and CO₂, shed light on the origin of fumarolic fluids by showing that effluents are mixture between fluids degassed from a magma body and the vapor generated at about 360 °C by the vaporization of hydrothermal liquids. A typical ‘andesitic’ water type (δD ∼ −20‰, δ¹⁸O ∼10‰) and a CO₂-rich composition (X_CO2∼0.4) has been inferred for the magmatic fluids, while for the hydrothermal component a meteoric origin and a CO₂ fugacity fixed by fluid-rock reaction at high temperatures have been estimated. In the time the fraction of magmatic fluids in the fumaroles increased (up to ∼0.5) at each seismic and ground uplift crisis (bradyseism) which occurred at Campi Flegrei, suggesting that bradyseismic crises are triggered by periodic injections of CO₂-rich magmatic fluids at the bottom of the hydrothermal system.     

Thermal history of the unsaturated zone at Yucca Mountain,Nevada, USA

Secondary calcite, silica and minor amounts of fluorite deposited in fractures and cavities record the chemistry, temperatures, and timing of past fluid movement in the unsaturated zone at Yucca Mountain, Nevada, the proposed site of a high-level radioactive waste repository. The distribution and geochemistry of these deposits are consistent with low-temperature precipitation from meteoric waters that infiltrated at the surface and percolated down through the unsaturated zone. However, the discovery of fluid inclusions in calcite with homogenization temperatures (T_h) up to ∼80 °C was construed by some scientists as strong evidence for hydrothermal deposition. This paper reports the results of investigations to test the hypothesis of hydrothermal deposition and to determine the temperature and timing of secondary mineral deposition. Mineral precipitation temperatures in the unsaturated zone are estimated from calcite- and fluorite-hosted fluid inclusions and calcite δ¹⁸O values, and depositional timing is constrained by the ²⁰⁷Pb/²³⁵U ages of chalcedony or opal in the deposits. Fluid inclusion T_h from 50 samples of calcite and four samples of fluorite range from ∼35 to ∼90 °C. Calcite δ¹⁸O values range from ∼0 to ∼22‰ (SMOW) but most fall between 12 and 20‰. The highest T_h and the lowest δ¹⁸O values are found in the older calcite. Calcite T_h and δ¹⁸O values indicate that most calcite precipitated from water with δ¹⁸O values between −13 and −7‰, similar to modern meteoric waters.     

Late Quaternary glaciation and equilibrium-line altitudes of the Mayan Ice Cap, Guatemala, Central America

The Sierra los Cuchumatanes (3837 m), Guatemala, supported a plateau ice cap and valley glaciers around Montaña San Juan (3784 m) that totaled ∼ 43 km² in area during the last local glacial maximum. Former ice limits are defined by sharp-crested lateral and terminal moraines that extend to elevations of ∼ 3450 m along the ice cap margin, and to ca. 3000-3300 m for the valley glaciers. Equilibrium-line altitudes (ELAs) estimated using the area-altitude balance ratio method for the maximum late Quaternary glaciation reached as low as 3470 m for the valley glaciers and 3670 m for the Mayan Ice Cap. Relative to the modern altitude of the 0°C isotherm of ∼ 4840 m, we determined ELA depressions of 1110-1436 m. If interpreted in terms of a depression of the freezing level during maximal glaciation along the modern lapse rate of − 5.3°C km^− 1, this ΔELA indicates tropical highland cooling of ∼ 5.9 to 7.6 ± 1.2°C. Our data support greater glacial highland cooling than at sea level, implying a high tropical sensitivity to global climate changes. The large magnitude of ELA depression in Guatemala may have been partially forced by enhanced wetness associated with southward excursions of the boreal winter polar air mass.     

An in situ XAS study of copper(I) transport as hydrosulfide complexes in hydrothermal solutions (25-592 °C, 180-600 bar): Speciation and solubility in vapor and liquid phases

Chloride and hydrosulfide are the principal ligands assumed to govern transport of copper in hydrothermal fluids. Existing solubility experiments suggest that Cu(I)-hydrosulfide complexes are dominant compared to chloride complexes at low salinities in alkaline solutions (H₂S_(aq)/HS⁻ pH buffer), and may be important in transporting Cu in low density magmatic vapors, potentially controlling the liquid-vapor partitioning of Cu. This study provides the first in situ evidence of the solubility of copper sulfides and the nature and structure of the predominant Cu species in sulfur-containing fluids at temperatures up to 592 °C and pressures of 180-600 bar. XANES and EXAFS data show that at elevated T (?200 °C), Cu solubility occurs via a linear Cu complex. At 428 °C in alkaline solutions, Cu is coordinated by two sulfur atoms in a distorted linear coordination (angle ∼150-160°). This geometry is consistent with the species predicted by earlier solubility studies. In addition, in situ measurements of the solubility of chalcocite in 2 m NaHS solutions performed in this study are in remarkably good agreement with the solubilities calculated using available thermodynamic data for Cu(I)-hydrosulfide complexes, also supporting the interpretation of speciation in these studies and validating the extrapolation of low-T thermodynamic properties for to high P-T. Data on phase separation for the 2 m NaHS solution show that while significant amounts of copper can be partitioned into the vapor phase, there is no indication for preferential partitioning of Cu into the vapor. This is consistent with recent partitioning experiments conducted in autoclaves by Pokrovski et al. (2008a) and Simon et al. (2006). XANES data suggest that the species present in the low density phase is very similar to that present in the high density liquid, i.e., , although Cu(HS)(H₂S)⁰ cannot be excluded on the basis of XAS data.     

O/O and D/H ratios of pedogenic kaolinite in a North American Cenomanian laterite: Paleoclimatic implications

Kaolinite, gibbsite and quartz are the dominant minerals in samples collected from two outcrops of a Cenomanian (∼95 Ma) laterite in southwestern Minnesota. A combination of measured yields and isotope ratios permitted mass balance calculations of the δD and δ¹⁸O values of the kaolinite in these samples. These calculations yielded kaolinite δD values of about −73‰ and δ¹⁸O values of about +18.7‰. The δD and δ¹⁸O values appear to preserve information on the ancient weathering system.If formed in hydrogen and oxygen isotope equilibrium with water characterized by the global meteoric water line (GMWL), the kaolinite δD and δ¹⁸O values indicate a crystallization temperature of 22 (±5) °C. A nominal paleotemperature of 22 °C implies a δ¹⁸O value for the corresponding water of −6.3‰. The combination of temperature and meteoric water δ¹⁸O values is consistent with relatively intense rainfall at that mid-paleolatitude location (∼40°N) on the eastern shore of the North American Western Interior Seaway. The inferred Cenomanian paleosol temperature of ∼22 °C is in general accord with published mid-Cretaceous continental mean annual temperatures (MAT) estimated from leaf margin analyses of fossil plants.When compared with results from a published GCM-based Cenomanian climate simulation which specifies a latitudinal sea surface temperature (SST) gradient that was either near modern or smaller-than-modern, the kaolinite paleotemperature of 22 °C is closer to the GCM-predicted MAT for a smaller equator-to-pole temperature difference in the mid-Cretaceous. Moreover, the warm, kaolinite-derived, mid-paleolatitude temperature of 22 °C is associated with proxy estimates of high concentrations of atmospheric CO₂ in the Cenomanian. The overall similarity of proxy and model results suggests that the general features of Cenomanian continental climate in that North American locale are probably being revealed.     

Diffusion of Ar in muscovite

Hydrothermal treatment of closely sized muscovite aggregates in a piston-cylinder apparatus induced ⁴⁰Ar^∗ loss that is revealed in ⁴⁰Ar/³⁹Ar step heating spectra. Age spectra and Arrhenius data, however, differ from that expected from a single diffusion length scale. A numerical model of episodic loss assuming the presence of multiple diffusion domains yields excellent fits between synthetic and actual degassing spectra. We used this model to isolate ⁴⁰Ar^∗ loss from the grains that remained intact during hydrothermal treatment at 10 kbar permitting calculation of diffusion coefficients in the temperature range 730-600 °C. Diffusion data generated in this manner yield an activation energy (E) of 63 ± 7 kcal/mol and frequency factor (D_o) of 2.3  cm²/s. Experiments at 20 kbar yield diffusivities lower by about an order of magnitude and correspond to an activation volume of ∼14 cm³/mol. Together, these parameters predict substantially greater retentivity of Ar in muscovite than previously assumed and correspond to a closure temperature (T_c) of 425 °C for a 100 μm radius grain cooling at 10 °C/Ma at 10 kbar (T_c = 405 °C at 5 kbar. Age and log (r/r_o) spectra for the run products show strong correlations indicating that muscovites can retain Ar diffusion boundaries and mechanisms that define their natural retentivity during vacuum step heating. This may permit the application of high resolution, continuous ⁴⁰Ar/³⁹Ar thermochronology to low grade, regionally metamorphosed terranes.     

Paleoenvironmental implications of concentration and C/C ratios of Fe(CO3)OH in goethite from a mid-latitude Cenomanian laterite in southwestern Minnesota

A mid-Cretaceous (∼95 Ma) laterite in southwestern Minnesota contains pisolites that consist primarily of gibbsite, quartz, and kaolinite with smaller amounts of goethite (α-FeOOH) and hematite. The presence of minor berthierine (an Fe(II) sheet silicate) suggests that this Cenomanian laterite experienced some degree of low temperature reductive diagenesis during its burial history. The prospects for extracting useful paleoenvironmental information from the pisolitic goethite were explored by studying measured mole fraction (X_m) and δ¹³C_m values of the Fe(CO₃)OH component in solid solution in the goethite using the method of incremental vacuum dehydration-decarbonation.Data arrays that occupy distinctly different domains in plots of δ¹³C_m vs. 1/X_m suggest the existence of two generations of goethite in the pisolites. The apparently younger generation of goethite (“generation 2”) evolves CO₂ at 170 °C, while the older generation (“generation 1”) evolves CO₂ at 220 °C. The distribution of the data suggests that generation 2 goethite is a proxy for mixing of CO₂ from three distinct CO₂ sources in an oxidative environment which post-dated the reductive diagenesis. The small amount of generation 1 goethite seems to have persisted through the reductive diagenesis, and nine of the generation 1 goethite data appear to define a proxy, two-endmember, soil CO₂ mixing line. Such two-component mixing is consistent with expectations for a highly evolved, carbonate-free laterite (i.e., the pre-diagenetic Cenomanian weathering system). The δ¹³C_m values of these nine data points range from −23.1‰ to −13.7‰, whereas X_m values range from 0.0007 to 0.0222. Linear regression of these nine data yields a slope of 0.0064, which corresponds to an ancient tropospheric CO₂ concentration of about 1900 ppmV.Isotopic data from pisolitic kaolinite indicate a paleotemperature of about 24 °C at a paleolatitude of ∼40°N. This is substantially warmer than modern continental temperatures at such latitudes and is consistent with published indications of a generally warmer mid-Cretaceous climate. Moreover, the correspondence of a warmer mid-Cretaceous climate with the inferred, relatively high concentration of Cenomanian tropospheric CO₂ (∼1900 ppmV) is consistent with the idea that variations of atmospheric CO₂ have a relation to climate change. The results of this study emphasize the importance of careful evaluation of incremental dehydration-decarbonation data from natural goethites to assess the possibility that more than one generation of goethite is present in a sample. However, the results also indicate that the carbon isotope information recorded in admixed goethite generations may be sorted out and used in paleoenvironmental interpretations.     

Organic geochemistry of Miocene source rocks from the Banat Depression (SE Pannonian Basin, Serbia)

The origin, depositional environment and maturity of petroleum source rocks were determined via conventional whole rock and biomarker analysis of samples from wells in the Banat Depression, where the most important Serbian oil and gas fields are located. The organic matter (OM) in organic-rich upper Tertiary siltstones and marls consists predominantly of Type II kerogen. Numerous biomarker parameters indicated mixed algal-terrestrial OM, related to a brackish or freshwater environment, whose salinity decreased from Middle to Upper Miocene. The OM was deposited under variable redox conditions, reducing to sub-oxic.The wells in the Banat Depression experienced variable high rates of rapid heating, providing an opportunity for examining the applicability of different thermal indicators in a hyperthermal basin. Rock-Eval and numerous biomarker parameters indicate that the main stage of oil generation begins at ca.130 °C and vitrinite reflectance (Rc) ca. 0.63% and reaches a maximum at ca.145-150 °C and Rc ca. 0.72-0.75%, while the late stage of oil generation starts at ca.155 °C and ca. Rc 0.78%, which corresponds, depending on geothermal gradient, to relative depths of 2100-2300 m, 2600-2900 m and 3050-3100 m, respectively. The naphthalene and phenanthrene maturity parameters proved to be less applicable than the biomarker ratios, particularly in the early to moderate maturation range. The newly proposed parameter C(14a)-homo-26-nor-17α(H)-hopane/C₃₀hopane (C₃₀HH/C₃₀H) proved applicable to a wide range of maturity.     

Low-temperature anisotropic diffusion of helium in zircon: Implications for zircon (U-Th)/He thermochronometry

In the last decade the zircon (U-Th)/He (ZHe) thermochronometer has been applied to a variety of geologic problems. Although bulk diffusion coefficients for He in zircon are available from laboratory step-heating experiments, little is known about the diffusion mechanism(s) and their dependence on the crystallographic structure of zircon. Here, we investigate the diffusion of He in perfectly crystalline zircon using atomistic simulation methods that provide insights into the structural pathways of He migration in zircon. Empirical force fields and quantum-mechanical calculations reveal that the energy barriers for He diffusion are strongly dependent on structure. The most favorable pathway for He diffusion is the [0 0 1] direction through the open channels parallel to the c-axis (, activation energy for tracer diffusion of a He atom along [0 0 1]). In contrast, energy barriers are higher in other directions where narrower channels for He diffusion are identified, such as [1 0 0], [1 0 1], and [1 1 0] (ΔE^∗ of 44.8, 101.7, and 421.3 kJ mol⁻¹, respectively). Molecular dynamics simulations are in agreement with these results and provide additional insight in the diffusion mechanisms along different crystallographic directions, as well as the temperature dependence. Below the closure temperature of He in zircon [T_c ∼ 180 °C, Reiners P. W., Spell T. L., Nicolescu S., and Zanetti K. A. (2004) Zircon (U-Th)/He thermochronometry: He diffusion and comparisons with Ar-40/Ar-39 dating. Geochim. Cosmochim. Acta68, 1857-1887], diffusion is anisotropic as He moves preferentially along the [0 0 1] direction, and calculated tracer diffusivities along the two most favorable directions differ by approximately five orders of magnitude (D_[001]/D_[100] ∼ 10⁵, at T = 25 °C). Above this temperature, He atoms start to hop between adjacent [0 0 1] channels, along [1 0 0] and [0 1 0] directions (perpendicular to the c-axis). The diffusion along [1 0 0] and [0 1 0] is thermally activated, such that at higher temperatures, He diffusion in zircon becomes nearly isotropic (D_[001]/D_[100] ∼ 10, at T = 580 °C). These results suggest that the anisotropic nature of He diffusion at temperatures near the closure temperature should be considered in future diffusivity experiments. Furthermore, care should be taken when making geologic interpretations (e.g., exhumation rates, timing of cooling, etc.) from this thermochronometer until the effects of anisotropic diffusion on bulk ages and closure temperature estimates are better quantified.     

An XAS study of the structure and thermodynamics of Cu(I) chloride complexes in brines up to high temperature (400 °C, 600 bar)

The transport and deposition of copper in saline hydrothermal fluids are controlled by the stability of copper(I) complexes with ligands such as chloride. Despite their role in the formation of most hydrothermal copper deposits, the nature and stability of Cu(I) chloride complexes in highly saline brines remains controversial. We present new X-ray absorption data (P = 600 bar, T = 25-400 °C, salinity up to 17.2 m Cl), which indicate that the linear (x = 1, 2) complexes are stable up to supercritical conditions. Distorted trigonal planar complexes predominate at room temperature and at high salinity (>3 m LiCl): subtle changes in the XANES spectrum with increasing salinity may reflect geometric distortions of this complex. Similar changes were observed in UV-Vis data [Liu, W., Brugger, J., McPhail, D.C., Spiccia, L., 2002. A spectrophotometric study of aqueous copper(I) chloride complexes in LiCl solutions between 100 °C and 250 °C. Geochim. Cosmochim. Acta66, 3615-3633], and were erroneously interpreted as a new species, . Our XAS data and ab-initio XANES calculations show that this tetrahedral species is not present to any significant degree in our solutions. The stability of the complexe decreases with increasing temperature; under supercritical conditions and in brines under magmatic-hydrothermal conditions (e.g., 15.58 m Cl, 400 °C, 600 bar), only the linear Cu(I) chloride complexes were observed. This result and the instability of the complex are also consistent with the recent ab-initio molecular dynamic calculations of Sherman [Sherman D. M.(2007) Complexation of Cu⁺ in hydrothermal NaCl brines: ab-initio molecular dynamics and energetics. Geochim. Cosmochim. Acta71, 714-722]. This study illustrates the power of the quantitative nature of XANES and EXAFS measurements for deciphering the speciation of weak transition metal complexes up to magmatic-hydrothermal conditions.The systematic XANES data are used to retrieve the formation constant for at 150 °C, which is in good agreement with the reinterpretation of the UV-Vis data of Liu et al. (Liu et al., 2002). At high temperatures (?400 °C), the solubility of chalcopyrite in equilibrium with hematite-magnetite-pyrite and K-feldspar-muscovite-quartz calculated with the new properties is lower than that calculated using the previous model, and the calculated solubilities are at the lower end of the range of values measured in brine inclusions from porphyry copper systems.     

Characterization of elemental release during microbe-basalt interactions at T = 28 °C

This study used batch reactors to characterize the rates and mechanisms of elemental release during the interaction of a single bacterial species (Burkholderia fungorum) with Columbia River Flood Basalt at T = 28 °C for 36 days. We primarily examined the release of Ca, Mg, P, Si, and Sr under a variety of biotic and abiotic conditions with the aim of evaluating how actively metabolizing bacteria might influence basalt weathering on the continents. Four days after inoculating P-limited reactors (those lacking P in the growth medium), the concentration of viable planktonic cells increased from ∼10⁴ to 10⁸ CFU (Colony Forming Units)/mL, pH decreased from ∼7 to 4, and glucose decreased from ∼1200 to 0 μmol/L. Mass-balance and acid-base equilibria calculations suggest that the lowered pH resulted from either respired CO₂, organic acids released during biomass synthesis, or H⁺ extrusion during uptake. Between days 4 and 36, cell numbers remained constant at ∼10⁸ CFU/mL and pH increased to ∼5. Purely abiotic control reactors as well as control reactors containing inert cells (∼10⁸ CFU/mL) showed constant glucose concentrations, thus confirming the absence of biological activity in these experiments. The pH of all control reactors remained near-neutral, except for one experiment where the pH was initially adjusted to 4 but rapidly rose to 7 within 2 days. Over the entire 36 day period, P-limited reactors containing viable bacteria yielded the highest Ca, Mg, Si, and Sr release rates. Release rates inversely correlate with pH, indicating that proton-promoted dissolution was the dominant reaction mechanism. Both biotic and abiotic P-limited reactors displayed low P concentrations. Chemical analyses of bacteria collected at the end of the experiments, combined with mass-balances between the biological and fluid phases, demonstrate that the absence of dissolved P in the biotic reactors resulted from microbial P uptake. The only P source in the basalt is a small amount of apatite (∼1.2%), which occurs as needles within feldspar grains and glass. We therefore conclude that B. fungorum utilized apatite as a P source for biomass synthesis, which stimulated elemental release from coexisting mineral phases via pH lowering. The results of this study suggest that actively metabolizing bacteria have the potential to influence elemental release from basalt in continental settings.     

Reduced sulfur species in a stratified seawater lake (Rogoznica Lake, Croatia); seasonal variations and argument for organic carriers of reactive sulfur

Over a period of a year, Hg⁰-reactive, total reduced sulfur species (RSS_T), as well as a non-volatile fraction that cannot be gas-stripped at pH ∼2 (RSS_NV), have been measured by voltammetry in a stratified, saline lake. In the hypolimnion, RSS_T is dominated by unusually high (up to 5 mM) dissolved divalent sulfur (S^−II), present as H₂S + HS⁻ and as inorganic polysulfides (H_xS_n^x−2). Less abundant RSS_NV is attributed to dissolved zero-valent sulfur (S⁰) in inorganic polysulfides. Assuming negligible contribution of organic S⁰ species in the hypolimnion, the equilibrium distribution of polysulfide ions is calculated; S₅²⁻ is found to predominate. In the epilimnion, all RSS_T consists of RSS_NV within analytical uncertainty. Through spring and summer, RSS_T and RSS_NV display little vertical or seasonal variation, but they increase dramatically when stratification breaks down in autumn. Based on decay rate, RSS during mixing events is attributed to dissolved S₈ from oxidation of sulfide and decomposition of inorganic polysulfides. This hypothesis quantitatively predicts precipitation of elemental sulfur in a year when colloidal sulfur was observed and predicts no precipitation in a year when it was not observed. Except during mixing events, the entire water column is undersaturated with respect to both rhombic sulfur and biologic sulfur, and the limited variations of RSS exclude hydrophobic and volatile aqueous S₈ as a major species. During such periods, RSS (typically 8 nM) may be associated with organic carbon, perhaps as adsorbed S₈ or as covalently bound polysulfanes or polysulfides. The hypolimnion is viewed as a zero-valent sulfur reactor that creates S⁰-containing, dissolved organic macromolecules during stable stratification periods. Some are sufficiently degradation-resistant and hydrophilic to be dispersed throughout the lake during mixing events, subsequently giving rise to ∼10⁻⁸ M RSS in the oxic water column. Voltammetrically determined RSS in oxic natural waters has often been described as “sulfide” or “metal complexed sulfide”, implying an oxidation state of S^−II; we argue that RSS in oxic Rogoznica Lake waters is mainly S⁰.