Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

Bed form-induced hyporheic exchange and geochemical hotspots

Small-scale bed form topographies control hyporheic exchange and biogeochemical processes within aquatic sediments, which ultimately affect water quality and nutrient cycling at the watershed scale. The impact of three-dimensional and small-scale bed form topographies on hyporheic exchange and solute mixing is investigated in the present work. The effect of bed form morphologies on the development of zones of enhanced reaction rates (i.e., hotspots) is also studied. A computational fluid dynamics model to simulate river flow over bed forms is combined with a subsurface flow and multicomponent reactive solute transport model. A wide variety of bed form topographies are generated using geometric models by varying parameters controlling curvature as well as bed form wavelength and amplitude. The results in this research suggest that out-of-phase bed forms generate more complex hyporheic flow patterns which reduce the efficiency of solute transformations. Higher phase shifts in bed form shapes result in overall higher average velocity, larger zones of enhanced pressure and reaction rates, and higher amounts of solute exchange. Moreover, the bed form shapes control the reaction process for a wide range of sediment conductivities. Here this study advances the understanding of the effects of complex and small scale morphological features on hyporheic exchange processes including the rate and spatio-temporal distribution of reaction hotspots.

54 ENVIRONMENTAL SCIENCES↗

Mineralogic and petrologic implications of Viking geochemical results from Mars - Interim report

Chemical results from four samples of Martian fines delivered to Viking landers 1 and 2 are remarkably similar in that they all have high iron; moderate magnesium, calcium, and sulfur; low aluminum; and apparently very low alkalies and trace elements. This composition is best interpreted as representing the weathering products of mafic igneous rocks. A mineralogic model, derived from computer mixing studies and laboratory analog preparations, suggests that Mars fines could be an intimate mixture of about 80% iron-rich clay, about 10% magnesium sulfate (kieserite), about 5% carbonate (calcite), and about 5% iron oxides (hematite, magnetite, maghemite, goethite). The mafic nature of the present fines (distributed globally) and their probable source rocks seems to preclude large-scale planetary differentiation of a terrestrial nature.

Baird, A. K.↗

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Data and scripts associated with a manuscript investigating dissolved organic matter and microbial community linkages across seven globally distributed rivers

This data package is associated with the publication “Meta-metabolome ecology reveals that geochemistry and microbial functional potential are linked to organic matter development across seven rivers” submitted to Science of the Total Environment. This data package includes the data necessary to replicate the analyses presented within the manuscript to investigate dissolved organic matter (DOM) development across broad spatial distances and within divergent biomes. Specifically, we included the Fourier transform ion cyclotron mass spectrometry (FTICR-MS) data, geochemistry data, annotated metagenomic data, and results from ecological null modeling analyses in this data package. Additionally, we included the scripts necessary to generate the figures from the manuscript. Complete metagenomic data associated with this data package can be found at the National Center for Biotechnology (NCBI) under Bioproject PRJNA946291. This dataset consists of (1) four folders; (2) a file-level metadata (flmd) file; (3) a data dictionary (dd) file; (4) a factor sheet describing samples; and (5) a readme. The FTICR Data folder contains (1) the processed Fourier transform ion cyclotron mass spectrometry (FTICR-MS) data; (2) a transformation-weighted characteristics dendrogram generated from the FTICR-MS data; and (3) the script used to generate all FTICR-MS related figures. The Geochemical Data folder contains (1) the single geochemistry data file and (2) the R script responsible for generating associated figures. The Metagenomic Data folder contains (1) annotation information across different levels; (2) carbohydrate active enzyme (CAZyme) information from the dbCAN database (Yin et al., 2012); (3) phylogenetic tree data (FASTAs, alignments, and tree file); and (4) the scripts necessary to analyze all of these data and generate figures. The Null Modeling Data folder contains (1) data generated during null modeling for each river and all rivers combined and (2) the R scripts necessary to process the data. All files are .csv, .pdf, .tsv, .tre, .faa, .afa, .tree, or .R.

54 ENVIRONMENTAL SCIENCES↗

Roughrider Carbon Storage Hub (Final Report)

The Roughrider Carbon Storage Hub was a 2-year project (October 2023 – September 2025) conducted by the Energy & Environmental Research Center (EERC) focused on advancing the feasibility of a commercial-scale carbon dioxide (CO 2 ) geologic storage hub in McKenzie County, North Dakota. The project’s objective was to investigate the potential that stacked storage complexes (multiple deep saline formations) can safely and economically store at least 50 million tonnes of CO 2 within 30 years. The captured CO 2 would be sourced from industrial emitters including project partner ONEOK, Inc.’s gas-processing plants and a planned gas-to-liquids facility. Drilling of the Roughrider 1 stratigraphic test well (14,979-ft total depth) was completed in November 2024. The wellbore intersected four candidate storage formations: Inyan Kara, Broom Creek, Mission Canyon, and Black Island–Deadwood. Operational challenges, including a stuck drill string, were resolved without long-term impact. A comprehensive logging and coring program was conducted, followed by successful well abandonment and site reclamation. Over 660 ft of 4-in. whole core was retrieved. Core plug samples were processed and analyzed for petrophysical and geochemical properties. Results confirmed promising porosity and permeability in the Inyan Kara and Broom Creek Formations and removal of the Mission Canyon and Black Island–Deadwood horizons from further investigation. Data derived from the logging and coring program were used to improve initial geologic models built from legacy data. CO 2 injection simulations showed that the Inyan Kara alone can feasibly store the target mass of CO 2 . Because of subtle differences in geologic structure and porosity trends between the formations, a stacked storage scenario using the Broom Creek and Inyan Kara Formations resulted in a larger overall plume area than using the Inyan Kara alone. Preliminary CO 2 pipeline routes from the industrial sources were mapped utilizing existing rights of way and evaluated for capacity and cost using U.S. Department of Energy Office of Fossil Energy and Carbon Management/National Energy Technology Laboratory models and U.S. Environmental Protection Agency emissions data. Integrating capture, transport, and storage cost estimates with policy incentives (e.g., 45Q credits) provided a total cost-per-ton analysis. Results indicate that the small scale of the volumes to be transported over the cumulative large distances does not support the project’s financial viability. However, the groundwork laid during this project from geological, regulatory, and social perspectives positions the Roughrider hub site as a promising candidate for commercial carbon storage in North Dakota, especially if the economy of scale is introduced for CO 2 transportation to the hub site.

01 COAL, LIGNITE, AND PEAT↗

Impact of Mineral Reactive Surface Area on Forecasting Geological Carbon Sequestration in a CO2-EOR Field

Mineral reactive surface area (RSA) is one of the key factors that control mineral reactions, as it describes how much mineral is accessible and can participate in reactions. This work aims to evaluate the impact of mineral RSA on numerical simulations for CO2 storage at depleted oil fields. The Farnsworth Unit (FWU) in northern Texas was chosen as a case study. A simplified model was used to screen representative cases from 87 RSA combinations to reduce the computational cost. Three selected cases with low, mid, and high RSA values were used for the FWU model. Results suggest that the impact of RSA values on CO2 mineral trapping is more complex than it is on individual reactions. While the low RSA case predicted negligible porosity change and an insignificant amount of CO2 mineral trapping for the FWU model, the mid and high RSA cases forecasted up to 1.19% and 5.04% of porosity reduction due to mineral reactions, and 2.46% and 9.44% of total CO2 trapped in minerals by the end of the 600-year simulation, respectively. The presence of hydrocarbons affects geochemical reactions and can lead to net CO2 mineral trapping, whereas mineral dissolution is forecasted when hydrocarbons are removed from the system.

58 GEOSCIENCES↗

NEWTS USGS Produced Waters Database

Data from the USGS Produced Waters Database on produced water compositions from across the US. Data has been filtered for completeness of data attributes. Data has been transposed and converted for ease of input into aqueous chemistry modeling software including OLI Systems and Geochemist's Workbench. Original data from: Blondes, M. S., Gans, K. D., Engle, M. A., Kharaka, Y. K., Reidy, M. E., Saraswathula, V., Thordsen, J. J., Rowan, E. L., & Morrissey, E. A. (2019). U.S. Geological Survey National Produced Waters Geochemical Database v2.3 [Data set]. U.S. Geological Survey. https://doi.org/10.5066/F7J964W8

Geochemists Workbench↗

Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation

Shallow subsurface microbial and geochemical processes in watersheds are dynamic and responsive to seasonal and long-term environmental change. At high elevation or latitude, such changes may occur during winter months when normal sampling is impossible, impractical, unsafe, or may actually alter the processes being studied. Yet accurate modeling of biogeochemical parameters in these environments requires sampling that can capture events and persistent cold climate trends. Gaps in current models exist where data have not been collected during extended snow- or ice-covered periods. The goals of our exploratory study were to examine the biogeochemical processes occurring in a well-studied high altitude watershed during winter when the system is largely inaccessible. We tested the hypothesis that during snow cover in the East River (ER) watershed, episodic excursions of microbial community structure and biogeochemical processes and concentrations fluctuated from values extrapolated from pre- and post-snow time periods. Through our research we identified previously unresolved subsurface-surface transport phenomena and a heretofore under-appreciated methane cycle occurring during wintertime months, thus providing valuable insight into watershed processes. Our research showed that autonomous sampling techniques, first used in marine systems, is extremely versatile and can be adapted for use in harsh continental settings to examine temporal changes in freshwater chemistry and microbial community assemblage. Further, we found that by sampling during periods when systems are hard to access, will lend insight and better accountability for global budgets for methane, an important greenhouse gas. By strengthening our knowledge of wintertime biogeochemical relationships in high altitude watersheds our research leads to more complete descriptions of cryptic processes and ultimately will inform and improve the reactive transport models that simulate processes in these systems.

54 ENVIRONMENTAL SCIENCES↗

Cogenetic Rock Fragments from a Lunar Soil: Evidence of a Ferroan Noritic-Anorthosite Pluton on the Moon

The impact that produced North Ray Crater, Apollo 16 landing site, exhumed rocks that include relatively mafic members of the lunar ferroan anorthositic suite. Bulk and mineral compositions indicate that a majority of 2-4 mm lithic fragments from sample 67513, including impact breccias and monomict igneous rocks, are related to a common noritic-anorthosite precursor. Compositions and geochemical trends of these lithic fragments and of related samples collected along the rim of North Ray Crater suggest that these rocks derived from a single igneous body. This body developed as an orthocumulate from a mixture of cumulus plagioclase and mafic intercumulus melt, after the plagioclase had separated from any cogenetic mafic minerals and had become concentrated into a crystal mush (approximately 70 wt% plagioclase, 30 wt% intercumulus melt). We present a model for the crystallization of the igneous system wherein "system" is defined as cumulus plagioclase and intercumulus melt. The initial accumulation of plagioclase is analogous to the formation of thick anorthosites of the terrestrial Stillwater Complex; however, a second stage of formation is indicated, involving migration of the cumulus-plagioclase-intercumulus-melt system to a higher crustal level, analogous to the emplacement of terrestrial massif anorthosites. Compositional variations of the lithic fragments from sample 67513 are consistent with dominantly equilibrium crystallization of intercumulus melt. The highly calcic nature of orthocumulus pyroxene and plagioclase suggests some reaction between the intercumulus melt and cumulus plagioclase, perhaps facilitated by some recrystallization of cumulus plagioclase. Bulk compositions and mineral assemblages of individual rock fragments also require that most of the mafic minerals fortned in close contact with cumulus plagioclase, not as separate layers. The distribution of compositions (and by inference, modes) has a narrow peak at anorthosite and a broader, larger peak at noritic anoilhosite. Characteristics of the samples and their geochemical trends imply an origin in a system that was large relative to the (unknown) size of the impact that produced the breccias of ferroan noritic-anorthosite composition that were excavated later by the formation of North Ray Crater, and they appear to be consistent with an origin of the suite within a perched plagioclase cumulate. If the Moon's crust formed by accumulation of plagioclase in a magma ocean, ferroan noritic anothosite, formed as an orthocumulate, is an alternative to extensive adcumulus formation of ferroan anorthosite (greater than 90 vol% plagiclase). This provides a relatively mafic ferroan anorthositic component (approximately 15 vol% mafics), which is required by mass-balance models of compositions of polymict lunar-crustal materials. The inferred bulk composition of the system of cumulus plagioclase and intercumulus melt is similar to that of ferroan regolith breccia MacAlpine Hills 88104/5, a lunar-highland meteorite, and may represent a common and widespread component of the Moon's early highland crust.

Jolliff, B. L.↗

Influence of geochemical reactions on the creep behavior of Mt. Simon sandstone

The capture and subsurface storage of carbon dioxide is a sustainable option that is currently pursued worldwide to mitigate greenhouse gas effect. However, predicting the long-term mechanical integrity of CO 2 underground storage systems remains a challenge. To address that question, it is essential to understand the influence of fluid-rock chemo-mechanical interactions on the long-term and on the time-dependent mechanical properties. In turn, the long-term mechanical response and the time-dependent mechanical behavior can be represented by the creep response. We investigate the impact of CO 2 -induced geochemical reaction on the creep response of Mt. Simon sandstone with a 50–400 micron grain size. We perform static and dynamic flow experiments on Mt. Simon sandstone specimens under geological conditions, at a temperature of T = 50 - 53 °C and for CO 2 pressures of P = 8:62; 17:2 MPa under both static flow and dynamic flow-through conditions. After aging, we employ creep indentation testing, high-resolution SEM-EDS, computer vision, machine learning, and micromechanics modeling to probe changes on the microstructure and mechanical properties. Following both static and dynamic flow-through experiments, we observe a 10–22% decrease in quartz volume fraction and an increase in both the microporosity (7–28%) and nanoporosity (60–65%). Additional CO 2 -induced microstructural changes include an enlargement of pore throats and the formation of channels. These observations point to the presence of CO 2 -induced K-feldspar dissolution and clay dissolution reactions. Here, the macroscopic creep behavior is logarithmic and the macroscopic creep modulus varies depending on the microporosity and the relative quartz and feldspar content. As a result of these geochemical reactions and of the related microstructural changes, a 55-60% decrease in the macroscopic logarithmic creep modulus is predicted.

58 GEOSCIENCES↗

Endothermic physiology of extinct megatooth sharks

The evolution of the extinct megatooth shark, Otodus megalodon , and its close phylogenetic relatives remains enigmatic. A central question persists regarding the thermophysiological origins of these large predatory sharks through geologic time, including whether O. megalodon was ectothermic or endothermic (including regional endothermy), and whether its thermophysiology could help to explain the iconic shark’s gigantism and eventual demise during the Pliocene. To address these uncertainties, we present unique geochemical evidence for thermoregulation in O. megalodon from both clumped isotope paleothermometry and phosphate oxygen isotopes. Our results show that O. megalodon had an overall warmer body temperature compared with its ambient environment and other coexisting shark species, providing quantitative and experimental support for recent biophysical modeling studies that suggest endothermy was one of the key drivers for gigantism in O. megalodon and other lamniform sharks. The gigantic body size with high metabolic costs of having high body temperatures may have contributed to the vulnerability of Otodus species to extinction when compared to other sympatric sharks that survived the Pliocene epoch.

60 APPLIED LIFE SCIENCES↗

On the composition of the lunar interior

There is now abundant geophysical and geochemical evidence suggesting that the moon has a thick plagioclase rich outer shell. This is most easily explained by early and extensive melting of a CaO and Al2O3 rich moon followed by fractional crystallization involving plagioclase flotation. Melilite is probably an important constituent of the interior. This model explains the seismic velocities, the mean density, and the moment of inertia of the moon. The moon is 73-88% high-temperature condensate.

Anderson, D. L.↗

Radar-anomalous, high-altitude features on Venus

Over nearly all of the surface of Venus the reflectivity and emissivity at centimeter wavelengths are about 0.15 and 0.85 respectively. These values are consistent with moderately dense soils and rock populations, but the mean reflectivity is about a factor of 2 greater than that for the Moon and other terrestrial planets. Pettingill and Ford, using Pioneer Venus reflectivities and emissivities, found a number of anomalous features on Venus that showed much higher reflectivities and much lower emissivities with both values approaching 0.5. These include Maxwell Montes, a number of high regions in Aphrodite Terra and Beta Regio, and several isolated mountain peaks. Most of the features are at altitudes above the mean radius by 2 to 3 km or more. However, such features have been found in the Magellan data at low altitudes and the anomalies do not exist on all high structures, Maat Mons being the most outstanding example. A number of papers have been written that attempt to explain the phenomena in terms of the geochemistry balance of weathering effects on likely surface minerals. The geochemists have shown that the fundamentally basaltic surface would be stable at the temperatures and pressures of the mean radius in the form of magnetite, but would evolve to pyrite and/or pyrrhotite in the presence of sulfur-bearing compounds such as SO2. Pyrite will be stable at altitudes above 4 or 5 km on Venus. Although the geochemical arguments are rather compelling, it is vitally important to rationally look at other explanations for radar and radio emission measurements such as that presented by Tryka and Muhleman. The radar reflectivity values are retrieved from the raw Magellan backscatter measurements by fitting the Hagfors' radar scattering model in which a surface roughness parameters and a normal incidence electrical reflectivity are estimated. The assumptions of the theory behind the model must be considered carefully before the results can be believed. These include that the surface roughness exists only at horizontal scales large compared to the wavelength, the vertical deviations are gaussianly distributed, there is no shadowing, and that the reflection occurs at the interface of two homogeneous dielectric half-spaces. Probably all these conditions are violated at the anomalous features under discussion. The most important of these is the homogeneity of the near surface of Venus, particularly in highlands. Under the assumptions of the theory, all of the radio energy is reflected by the impedance jump at the very boundary. However, in heterogeneous soil some fraction of the illuminating energy is propagated into the soil and then scattered back out by impedance discontinuities such as rock, voids, and cracks. In light soils, the latter effect can overwhelm the scattering effects of the true surface and greatly enhance the backscatter power, suggesting a much higher value of an effective dielectric constant that would be estimated from Hagfors' model.

Muhleman, Duane O.↗

Differentiation of Asteroid 4 Vesta: Core Formation by Iron Rain in a Silicate Magma Ocean

Geochemical observations of the eucrite and diogenite meteorites, together with observations made by NASA's Dawn spacecraft while orbiting asteroid 4 Vesta, suggest that Vesta resembles H chondrites in bulk chemical composition, possible with about 25 percent of a CM-chondrite like composition added in. For this model, the core is 15 percent by mass (or 8 percent by volume) of the asteroid, with a composition of 73.7 percent by weight Fe, 16.0 percent by weight S, and 10.3 percent by weight Ni. The abundances of moderately siderophile elements (Ni, Co, Mo, W, and P) in eucrites require that essentially all of the metallic phase in Vesta segregated to form a core prior to eucrite solidification. The combination of the melting phase relationships for the silicate and metal phases, together with the moderately siderophile element concentrations together require that complete melting of the metal phase occurred (temperature is greater than1350 degrees Centigrade), along with substantial (greater than 40 percent) melting of the silicate material. Thus, core formation on Vesta occurs as iron rain sinking through a silicate magma ocean.

Kiefer, Walter S.↗

A Case for Nebula Scale Mixing Between Non-Carbonaceous and Carbonaceous Chondrite Reservoirs: Testing the Grand Tack Model with Chromium Isotopic Composition of Almahata Sitta Stone 91A

There is an increasing number of Cr-O-Ti isotope studies that show solar system materials are divided into two main populations, one carbonaceous chondrite (CC)-like and the other is non-carbonaceous (NC)-like, with minimal mixing attributed to a gap opened in the protoplanetary disk due to Jupiter's formation. The Grand Tack model suggests there should be large-scale mixing between S- and C-type asteroids, an idea supported by our recent work on chondrule (Delta)17O-ε54Cr isotope systematics. The Almahata Sitta (AhS) meteorite provides a unique opportunity to test the Grand Tack model. The meteorite fell to Earth in October 2008 and has been linked to the asteroid 2008 TC3 which was discovered just prior to the fall of the AhS stones. The AhS meteorite is composed of up to 700 individual pieces with approx.140 of those pieces having some geochemical and/or petrologic studies. Almahata Sitta is an anomalous polymict ureilite with other meteorite components, including enstatite, ordinary, and carbonaceous chondrites with an approximate abundance of 70% ureilites and 30% chondrites. This observation has lead to the suggestion that TC3 2008 was a loosely aggregated rubble pile-like asteroid with the non-ureilite sample clasts within the rubble-pile. Due to the loosely-aggregated nature of AhS, the object disintegrated during atmospheric entry resulting in the weakly held clasts falling predominantly as individual stones in the AhS collection area. However, recent work has identified one sample of AhS, sample 91A, which may represent two different lithologies coexisting within a single stone. The predominate lithology type in 91A appears to be that of a C2 chondrite based on mineralogy but also contains olivine, pyroxene, and albite that have ureilite-like compositions. Previous Cr isotope investigations into AhS stones are sparse and what data is available show nearly uniform isotopic composition similar to that of typical ureilites with negative ε54Cr values.

Sanborn, M. E.↗

Diamondites: Evidence for a Distinct Tectono-Thermal Diamond-Forming Event Beneath the Kaapvaal Craton

The petrogenesis and relationship of diamondite to well-studied monocrystalline and fibrous diamonds are poorly understood yet would potentially reveal new aspects of how diamond-forming fluids are transported through the lithosphere and equilibrate with surrounding silicates. Of twenty-two silicate- and oxide-bearing diamondites investigated, most yielded garnet intergrowths (n = 15) with major element geochemistry (i.e. Ca-Cr) classifying these samples as low-Ca websteritic or eclogitic. The garnet REE patterns fit an equilibrium model suggesting the diamond-forming fluid shares an affinity with high-density fluids (HDF) observed in fibrous diamonds, specifically on the join between the saline–carbonate end- members. The δ13C values for the diamonds range from -5.27 to -22.48 ‰ (V-PDB) with δ18O values for websteritic garnets ranging from +7.6 to +5.9 ‰ (V-SMOW). The combined C-O stable isotope data support a model for a hydrothermally altered and organic carbon-bearing subducted crustal source(s) for the diamond- and garnet-forming media. The nitrogen aggregation states of the diamonds require that diamondite-formation event(s) pre-dates fibrous diamond-formation and post-dates most of the gem monocrystalline diamond-formation events at Orapa. The modelled fluid compositions responsible for the precipitation of diamondites match the fluid-poor and fluid-rich (fibrous) monocrystalline diamonds, where all grow from HDFs within the saline-silicic-carbonatitic ternary system. However, while the nature of the parental fluid(s) share a common lithophile element geochemical affinity, the origin(s) of the saline, silicic, and/or carbonatitic components of these HDFs do not always share a common origin. Therefore, it is wholly conceivable that the diamondites are evidence of a distinct and temporally unconstrained tectono-thermal diamond-forming event beneath the Kaapvaal craton.

S Mikhail↗

Martian volatiles - Their degassing history and geochemical fate

Implications of the Ar-40 content of the Martian atmosphere for Martian volatile history are considered. It is assumed that the potassium content of both Mars and earth is similar to that of normal chondrites, that Mars possesses a water content similar to that of normal chondrites, and that earth's water content is equal to or somewhat less than Mars'. A simple earth-analogous Mars degassing model is outlined in which the Martian surface volatile inventory is presumed to be identical to earth's, but scaled to the smaller mass and surface area of Mars. Physical storage of volatiles in the Martian regolith is examined along with the occurrence of volatile-containing mineral phases in the regolith, exospheric escape as a massive volatile sink, and the effect of the time history of degassing on the surface volatile inventory. Implications of earth's degassing history for Mars are assessed, and models describing the overall intensity of Mars degassing are presented. It is concluded that the maximum amounts of water and other volatiles that could be stored in the Martian regolith are marginally compatible with those required by the earth-analogous model, although a lower atmospheric Ar-40 concentration and regolith volatile inventory would be easier to reconcile with observational constraints.

Fanale, F. P.↗

Climatic consequences of very high CO2 levels in Earth's early atmosphere

Earth has approximately 60 bars of carbon dioxide tied up in carbonate rocks, or roughly 2/3 the amount of CO2 of the atmosphere of Venus. Two different lines of evidence, one based on thermodynamics and the other on geochemical cycles, indicate that a substantial fraction of thes CO2 may have resided in the atmosphere during the first few hundred million years of the Earth's history. A natural question which arises is whether this much CO2 would have resulted in a runaway greenhouse effect. One dimensional radiative/convective model calculations presented showed that the surface temperature of a hypothetical primitive atmosphere containing 20 bars of CO2 was less than 100 C; thus no runaway greenhouse effect would have occurred. The climatic stability of the early atmosphere is a consequence of three factors: reduced solar luminosity at that time, an increase in planetary albedo caused by Rayleigh scattering by CO2, and the stabilizing effects of a moist convection. The latter two factors are sufficient to prevent a CO2 induced runaway greenhouse effect on the present Earth as well, for CO2 levels up to 100 bars. Further studies are being undertaken to determine whether a runaway greenhouse effect could have occurred during the latter stages of the accretion process and, if so, whether it would have collapsed one the influx of material slowed down.

Katsing, J. F.↗