Changes in the structures of mixed-layer illite-smectite during flooding of terrigenous oil reservoirs

A difference-spectrum method is proposed for the qualitative assessment of changes of illite-smectite structures accompanying the flooding of oil reservoirs. The method permits one to get an open system and reduce the application of procedures based on Markov’s chain formalism. A computer simulation is made to obtain spectra by subtracting the spectrum of an ethylene glycol-saturated sample from the spectrum of an air-dried preparation throughout the entire range of concentrations of illite and smectite components with a short-range order factor R = 0 or R = 3.It has been established that only in the presence of filtration are the maximum and minimum of the spectra in the range of 12.5-9.4 A complicated by a number of local extrema, whose position is specified by the structure of intermediate phases. The flooding process first involves mixed-layer phases with R = 0, leading to a partial segregation of the structures into phase with one and two networks of interlayer H₂O. When the secondary mica particles break, phases with R = 3 appear along the boundaries of nanoblocks, first only with 1 H₂O and then only with 2 H₂O in labile interspaces. Their coexistence with the phases R = 0 in the sample proves the existence of percolation effects due to two-phase filtration in the porous medium. The fully flooded reservoir is always dominated by a mechanical mixture of illite-smectite phases of different nature with R = 0 and with different ratios of components. Transformation of mica that can drastically reduce oil production begins long before the appearance of flooding zones, which are revealed by standard logging methods.     

An experimental study of the partitioning of a rare earth element (Gd) in the system diopside—aqueous vapour

The partitioning of Gd in the experimental system diopside-aqueous vapor as a function of temperature, pressure, composition of the phases, time, grain size, solid-liquid ratio and Gd concentration has been investigated. A radioactive tracer measurement was used to determine Gd concentration in the separated phases. Diposides were reacted with aqueous vapor containing tracer Gd and reversibility was tested by reacting Gd-doped diopsides with pure aqueous vapor. Equilibration of Gd between the bulk of the diopside and the liquid was found to be limited by the slow rate of Gd diffusion in diopside, maximum value of D = 2 × 10^?15cm²sec^?1 at 800°C and 1 kb. Depending on whether the diopside was previously synthesized or synthesized from an oxide mix during the experiment, Gd concentrations were zoned in the crystal such that higher concentrations existed at the edges or center, respectively. Equilibrium is difficult to achieve in these experiments, but at the optimum experimental conditions for equilibration, the Gd diopside-aqueous vapor distribution coefficient is 20 ± 6 (800°C, 1 kb) in approximate agreement with previous results of 55 ± 23. Changing the composition of the aqueous vapor indicated that possible mechanisms for Gd substitution included coupling of Gd³⁺ with H⁺ or Na⁺ replacing 2Ca²⁺, or substitution of 2Gd³⁺ for 3Ca²⁺ with formation of a cation vacancy.     

蒙皂石—伊利石演化及有关问题的理论探讨

摘要：作者认为蒙皂石—伊利石的转化是多因素、多机制作用的结果．不同环境、不同机制形成的伊／蒙混层矿物,有序化出现时的伊利石层含量不同。伊／蒙混层由无序到有序是一个连续的渐变过程,存在过渡型混层矿物,多数伊／蒙混层矿物都是以某一类型为主的复杂集合体。作者指出分散蒙皂石是影响无序型伊／蒙混层矿物混层比计算精度的主要因素。     

Rare earth element fractionation in authigenic illite–smectite from Late Permian clastic rocks,Bowen Basin,Australia: implications for physico-chemical environments of fluids during illitization

REE analyses were performed on authigenic illitic clay minerals from Late Permian mudrocks, sandstones and bentonites from the Bowen Basin (Australia). The mixed-layer illite–smectite exhibit REE patterns with an obvious fractionation of the HREE from the LREE and MREE, which is an apparent function of degree of illitization reaction. The highly illitic (R≥3) illite–smectite from the northern Bowen Basin show a depletion of LREE relative to the less illitic (R=0 and 1) clays. In contrast, an enrichment of HREE for the illite-rich clays relative to less illitic clays is evident for the southern Bowen Basin samples. The North American Shale Composite-normalized (La/Lu)_sn ratios show negative correlations with the illite content in illite–smectite and positive correlations with the δ¹⁸O values of the clays for both the northern and southern Bowen Basin samples. These correlations indicate that the increasing depletion of LREE in hydrothermal fluids is a function of increasing water/rock ratios in the northern Bowen Basin. Good negative correlations between (La/Lu)_sn ratios and illite content in illite–smectite from the southern Bowen Basin suggest the involvement of fluids with higher alkalinity and higher pH in low water/rock ratio conditions. Increasing HREE enrichment with δ¹⁸O decrease indicates the effect of increasing temperature at low water/rock ratios in the southern Bowen Basin.Results of the present study confirm the conclusions of some earlier studies suggesting that REE in illitic clay minerals are mobile and fractionated during illitization and that this fact should be considered in studies of sedimentary processes and in identifying provenance. Moreover, our results show that REE systematic of illitic clay minerals can be applied as an useful technique to gain information about physico-chemical conditions during thermal and fluid flow events in certain sedimentary basins.     

Statistical mechanics of coupled solid solutions in the dilute limit

A statistical mechanical analysis of the limiting laws for coupled solid solutions shows that the random model, in which the configurational entropy is calculated as if atoms mix randomly on each crystallographic site, is correct as a first approximation. In coupled solid solutions, since atoms of different valence substitute on the same sites, significant short-range order which reduces the entropy can be expected. A first-order correction is rigorously obtained for the entropy in dilute binary short-range ordered coupled solid solutions: $$\bar S^{{\text{XS}}} {\text{/R = }}Q\left( {{\text{e}}^{--H_{\text{A}} /{\text{R}}T} \left( {\frac{{H_{\text{A}} }}{{{\text{R}}T}} + 1} \right) - 1} \right)N_2^a N_4^b ,$$ where Q is the number of positions an associated cation pair can assume per formula unit, H _A is the association energy per formula unit, and N ₂ ^a and N ₄ ^b are the site occupancy fractions for atoms 2 and 4 that are dilute on sites a and b. S ^XS is the configurational entropy minus the random model entropy. Aluminous pyroxenes on the joints diopside-jadeite and diopside-CaTs are examined as examples. A generalization for dilute multiple component solutions, including possible long-range ordering variations is given by: $$\frac{{\bar S^{{\text{XS}}} }}{{\text{R}}}{\text{ = }}\sum\limits_i {\sum\limits_j {\sum\limits_k {Q_i } } \left( {{\text{e}}^{--H_{\text{A}}^{j{\text{ }}k{\text{, }}i} /{\text{R}}T} \left( {\frac{{H_{\text{A}}^{j{\text{ }}k{\text{, }}i} }}{{{\text{R}}T}} + 1} \right) - 1} \right)N_j^l N_k^m ,} $$ where i labels each crystallographically distinct pair, j and k label atomic species, l and m label crystallographic sites, and the N's are site occupancy fractions for the solute atoms. A total association model is examined as well as the partial association and random models. Real solution behavior must lie between the total association model and the random model. Molecular models in which the ideal activity is proportional to a mole fraction, which in itself is not always unambiguously defined, do not lie in this range and furthermore have no physical justification.     

Vibrational spectroscopy of beryllium aluminosilicates: Heat capacity calculations from band assignments

Far-infrared, mid-IR, and Raman powder spectra were measured on six phases (bromellite, chrysoberyl, phenakite, bertrandite, beryl, and euclase) in the system BeO-Al₂O₃-SiO₂-H₂O. A single-crystal absorption spectrum of IR fundamentals in beryl is also presented, which more closely resembles the powder absorption spectrum than it does absorption spectra calculated from single-crystal reflection data. Assignments of the SiO₄ and BeO₄ internal vibrations are made in accordance with each mineral's symmetry and composition and by comparison to structural analogs. Heat capacities C _v calculated for these partial band assignments agree with C _v derived from experimental C _p for all six phases, provided that Kieffer's (1979c) model is slightly modified to correctly enumerate both Si-O and Be-O stretching modes in the high frequency region (>750 cm^?1). Si-O stretching bands were found to out-number Be-O stretching modes in the high-energy region of the vibrational spectra with two exceptions: (1) For those phases containing oxygen ions not coordinated to silicon, vibrations occurring at v>1,080 cm^?1 that are attributable to Be-O (H) stretching must be treated separately in the model in order to calculate C _v accurately. (2) Minerals consisting entirely of interlocking Si and Be tetrahedra (i.e., phases without Al or OH) can be modeled by one optic continuum representing all optical modes. These results, along with the occurrence of very low energy lattice vibrations for Be-silicates within Al, suggests that although Be-O bonds are generally weaker than neighboring Si-O bonds, Be mimics the network-forming characteristic of Si to a limited extent.     

The dissolution kinetics of shallow marine carbonates in seawater: A laboratory study

The dissolution kinetics of shallow water marine carbonates (low-Mg calcite, aragonite and Mg-calcites) were investigated in seawater (S = 35) at 25°C and a P_CO2 of 10^?2.5 atm. using the pH-stat method. Carbonate dissoluton rates (μmoles g^?1 hr^?1) fit the empirical kinetic expression, R = k(1 - Ω)ⁿ, where R = dissolution rate, k = rate constant, Ω = saturation state, and n = order of reaction. Reaction orders were near 2.9 for low-Mg calcites, 2.5 for aragonites and 3.4 for Mg-calcites.The rate constant, k, expressed as μmoles g^?1 hr^?1, varied by nearly a factor of ten for the different samples, reflecting differences in amount of reactive surface area. Reactive surface area of the biogenic phases ranged from 0.3% to 66% of the total surface area determined by the BET gas adsorption method. The discrepancy between reactive and total surface area was greatest for samples with high BET surface areas (> 1 m² g^?1) and delicate microstructures.Relative dissolution rates of the various biogenic carbonates as a function of seawater calcium carbonate ion molal product (IMP) were related to both mineral stability and grain microstructure. In seawater undersaturated with respect to aragonite, finely crystalline aragonites dissolved more rapidly than thermodynamically less stable high Mg-calcites (15–18 mole% MgCO₃) with lower reactive surface areas. Therefore, under certain conditions, differences in grain microstructural complexity can override thermodynamic constraints and lead to selective dissolution of a thermodynamically more stable mineral phase.     

Thermochemistry of glasses along joins of pyroxene stoichiometry in the system Ca2Si2O6-Mg2Si2O6-Al4O6

Enthalpies of solution in 2PbO · B₂O₃ at 974 K have been measured for glasses along the joins Ca₂Si₂O₆ (Wo)-Mg₂Si₂O₆ (En) and Mg₂Si₂O₆-MgAl₂SiO₆ (MgTs). Heats of mixing are symmetric and negative for Wo-En with W_H = ?31.0 ± 3.6 kJ mol^?. Negative heats of mixing were also found for the En-MgTs glasses (W_H = ?33.4 ± 3.7 kJ mol^?).Enthalpies of vitrification of pyroxenes and pyroxenoids generally increase with decreasing alumina content and with decreasing basicity of the divalent cation.Heats of mixing along several glassy joins show systematic trends. When only non-tetrahedral cations mix (outside the aluminosilicate framework), small exothermic heats of mixing are seen. When both nontetrahedral and framework cations mix (on separate sublattices, presumably), the enthalpies of mixing are substantially more negative. Maximum enthalpy stabilization near compositions with Al/Si ≈ 1 is suggested.     

Surface destabilization and laboratory-induced non-stoichiometry in San Carlos olivine

Annealing experiments on natural olivine (Mg_1-x Fe _x )₂SiO₄ (with x≈0.11) crystals (San Carlos, Arizona, spinel-lherzolite context) have been performed between T=1,100° C and 1,500° C for oxygen partial pressures pO₂=10^?3 to 10^?13 bar and times of 1 to 140 h in CO/CO₂ or H₂/H₂O gas mixtures. Even specimens annealed within the T-pO₂ theoretical stability field (TSF) calculated for stoichiometric olivine (Nitsan 1974) show systematic alterations developed within the first few microns of the surface of the crystals. Pyroxene crystals or melt form on the original olivine surface even at T=1,100° C, with preference of pyroxene when T<1,350° C and SiO₂-rich glass if T>1,350° C. This glass (rhyolite-like) can concentrate calcium from the starting olivine, and aluminum when Cr-Al spinels are present as inclusions. These observations are in contradiction with the TSF. They are obviously due to the presence of platinum used as a container of our samples, even if the contact between olivine and platinum is very weak. Rapid surficial diffusion of iron toward platinum (or via the gas phase) induces a Fe-depleted surface. According to the TSF, this more forsteritic surface should have a wider pO₂ range of stability. This is not the case, just because this situation is largely out of equilibrium. This iron loss induces a departure from cationic stoichiometry: (Mg, Fe)_2(1?δ), SiO₄ with δ small and positive. We extended the model that Nakamura and Schmalzried (1983) (N.S.) developed for fayalite (x=1) to our natural olivine composition, under the assumption that the majority defects are magnesium vacancies, Fe³⁺ occupying octahedral and tetrahedral sites, and the more complex neutral defect corresponding to Coulombic attraction between neighboring Fe³⁺ ions. We have recalculated the olivine stability field in pO₂ vs. δ space at T=1,300° C using this model for x≈0.1 (its extreme limit of validity) and conclude that olivine is stable only in a very narrow range in pO₂ which depends on δ. The calculation shows also that when olivine has nearly cationic stoichiometry (δ=0) as we believe for our starting material, the pO₂ range of stability is narrower than indicated by the TSF. In particular, it explains why Fe precipitates from the olivine (δ=0) (in absence of any other precipitation of SiO₂-rich phases) between 10^?11 and 10^?13 bar at 1,300° C; this was not predicted by the TSF. Magnetite or iron precipitates also coexist with SiO₂-rich exsolutions or pyroxene when pO₂ is close to the upper or lower boundaries of the TSF, respectively. The N.S. model may have important implications for the interpretation of the existence of partial melting and/or the low-viscosity/low velocity zone in the upper mantle.     

Template plotting

The template transformation allows rapid straight line plotting of fluid-bearing reactions for conditions exceeding the two phase field in water-carbon dioxide. The system CaO-MgO-SiO₂-H₂O-CO₂ with the six solid phase tale, tremolite, quartz, magnesite, calcite, and dolomite is treated at 2Kb. These phases, in T-NH₂O space at constant pressure, participate in two invariant points and eight reactions. Published T-NH₂O data for three of these reactions at 2Kb are transformed using the template into LnK-1/T space to yield the parameters of the best (by eye) straight lines. The parameters of the other five reaction lines are derived by linear combination of these parameters. The reactions belonging to one invariant point in T-NH₂O space will be straight lines that intersect at a single point in ln K-1/T space. This point is not the transformation of the invariant point, and it occurs because of the linear relationship between the reactions. Errors need to be considered in transforming the system back into T-NH₂O space, especially when reactions lines with shallow slopes are involved.     

Some chemographic relationships in n-component systems

A petrologic problem of fundamental importance is to determine whether 2 or more mineral assemblages can be related to one another by continuous or discontinuous facies changes, or whether their bulk compositions occupy non-overlapping regions of composition space. A general method is developed by which 2 regions of n-dimensional space whose vertices are defined by the phases present are tested for compositional overlap. This is accomplished by generating mass balance equations of the type:

$\text{∑}\text{i = 1}\text{m}\text{a}{}_{\mathrm{i}}\text{A}{}_{\mathrm{i}}\text{=}\text{∑}\text{j = 1}\text{k}\text{b}{}_{\mathrm{j}}\text{B}{}_{\mathrm{j}}$

where A_i is the ith phase in one region and B_j is the th phase in the other. If any such equation satisfies the requirement that the sign of each a_i is the same, and that the sign of each b_j is the opposite for all i, j such that: k + m = n + 1 then the 2 regions overlap in phase space.By eliminating all overlapping assemblages in a given set, the bivariant fields bounded by univariant equilibria in n-dimensional systems are completely specified. All bulk compositions are considered within the space defined by the phases that participate in the bounding reactions. An extension of the method generates in sequence all bivariant fields and associated reactions about any invariant point. A further extension is applied to multi-system analysis.     

The heat-capacity of ilmenite and phase equilibria in the system Fe-T-O

Low temperature adiabatic calorimetry and high temperature differential scanning calorimetry have been used to measure the heat-capacity of ilmenite (FeTiO₃) from 5 to 1000 K. These measurements yield S₂₉₈⁰ = 108.9 J/(mol · K). Calculations from published experimental data on the reduction of ilmenite yield Δ₂₉₈⁰(I1) = ?1153.9 kJ/(mol · K). These new data, combined with available experimental and thermodynamic data for other phases, have been used to calculate phase equilibria in the system Fe-Ti-O. Calculations for the subsystem Ti-O show that extremely low values of $\text{?O}{}_2$ are necessary to stabilize TiO, the mineral hongquiite reported from the Tao district in China. This mineral may not be TiO, and it should be re-examined for substitution of other elements such as N or C. Consideration of solid-solution models for phases in the system Fe-Ti-O allows derivation of a new thermometer/oxybarometer for assemblages of ferropseudobrookite-pseudobrookite_ss and hematite-ilmenite_ss. Preliminary application of this new thermometer/oxybarometer to lunar and terrestrial lavas gives reasonable estimates of oxygen fugacities, but generally yields subsolidus temperatures, suggesting re-equilibration of one or more phases during cooling.     

Lithium and boron isotopes in illite-smectite: The importance of crystal size

Clay minerals record chemical data about the past, acting like natural computer memory chips. To retrieve the data we must understand how they are stored. To achieve this we have examined the isotopic information revealed by two trace elements, lithium and boron, that are incorporated into the common clay minerals illite-smectite (I-S) during diagenesis. We used hydrothermal experiments at 300°C, 100 MPa, to speed up the reaction of smectite to illite that normally occurs during slow (10-100 Ma) sediment burial. During illitization, Li substitutes into the octahedral sites and B enters the tetrahedral sites of the silicate framework. Both Li and B are also adsorbed in the interlayer of smectite, but Li is preferred over B in the exchange sites. To determine the equilibrium isotope fractionation of the two trace elements it is important to remove these adsorbed interlayer species. By measuring the isotopic composition of Li and B in the silicate framework during reaction, we can address the relative timing of element exchange in the different crystallographic sites. Furthermore, because illitization of smectite is a crystal growth process (not an isomorphous replacement) we have examined the effect of crystal size on the isotope fractionation.The results show that Li and B approach an isotopic steady state when R1 ordering occurs, long before oxygen isotopes equilibrate with the fluid. The isotopic fractionation (α_{mineral-water}) for Li (0.989) is similar to that for B (0.984) at 300°C. However, when separated into <0.2, 0.2-2.0, and >2.0 μm fractions, there are significant differences in measured isotope ratios by as much as 9‰. Crystal growth mechanisms and surface energy effects of nanoscale crystals may explain the observed isotopic differences. The fact that different crystals equilibrate at different rates (based on size) may be applied to natural samples to reveal the changing paleofluid history, provided we understand the conditions of equilibrium. This has very important implications for the interpretation of diagenetic environments, fluid flow, and surficial geochemical cycling.     

Sorption and distribution of adsorbed metals in three soils of India

The mobility and bioavailability of heavy metals depends on the metal retention capacity of soil and also on the geochemical phases with which metals are associated. Laboratory batch experiments were carried out to study the sorption and distribution of Cd, Ni and Pb in 3 soils differing in their physicochemical properties from India: Oxyaquic Haplustalf (SL1), Typic Haplustalf (SL2) and Typic Haplustert (SL3). The heavy metal adsorption was studied by isotherms and the distribution coefficient (K_D) for each metal was obtained from the linear regressions of the concentration of metal remaining in equilibrium solution and the amount adsorbed. In general, the sorption capacity for all the metals decreased in the order: SL3>SL2>SL1. Among metals, the sorption capacity in all the soils decreased in the order: Pb>>Ni>Cd. Distribution of sorbed metals at various equilibrating concentrations was studied by sequential extraction. Results showed significant differences in the distribution of metals in these soils. At higher additions (such as 200 μM l⁻¹) most of the metals were extracted in their more mobile fractions, exchangeable and/or inorganic in contrast to their original partitioning in soils, where they were preferentially associated with the less mobile residual fraction. Largest percentages of metals extracted in the exchangeable fraction corresponded to those soil–metal systems with smaller K_D values, e.g. Cd, Ni and Pb in SL1 and Cd and Ni in SL2. In neutral and alkaline soils (SL2, pH=7.1, and SL3, pH=8.6) Pb was predominantly extracted from the inorganic fractions and this corresponded to higher K_D values for Pb in these soils. The predominance of metals associated with the exchangeable fraction together with low K_D values indicates higher mobility of metals retained in the acidic soil (SL1, pH=5.2) compared with the others.     

A Bonded Particle Model Simulation of Shear Strength and Asperity Degradation for Rough Rock Fractures

Different failure modes during fracture shearing have been introduced including normal dilation or sliding, asperity cut-off and degradation. Attempts have been made to study these mechanisms using analytical, experimental and numerical methods. However, the majority of the existing models simplify the problem, which leads to unrealistic results. With this in mind, the aim of this paper is to simulate the mechanical behaviour of synthetic and rock fracture profiles during direct shear tests by using the two-dimensional particle flow computer code PFC2D. Correlations between the simulated peak shear strength and the fracture roughness parameter D _R1 recently proposed by Rasouli and Harrison (2010) are developed. Shear test simulations are carried out with PFC2D and the effects of the geometrical features as well as the model micro-properties on the fracture shear behaviour are studied. The shear strength and asperity degradation processes of synthetic profiles including triangular, sinusoidal and randomly generated profiles are analysed. Different failure modes including asperity sliding, cut-off, and asperity degradation are explicitly observed and compared with the available models. The D _R1 parameter is applied to the analysis of synthetic and rock fracture profiles. Accordingly, correlations are developed between D _R1 and the peak shear strength obtained from simulations and by using analytical solutions. The results are shown to be in good agreement with the basic understanding of rock fracture shear behaviour and asperity contact degradation.     

Natrolite,part I: Refinement of high-order data,separation of internal and external vibrational amplitudes from displacement parameters

A single crystal of natrolite, Na₂Al₂Si₃O₁₀·2H₂O, was studied by X-ray diffraction methods at room temperature. The intensities were measured with MoKα radiation (λ = 0.7107 Å) in a complete sphere of reflection up to sin θ/λ = 0.903 Å^?1. The structure was refined in the orthorhombic space group Fdd2 with a = 18.2929 (7) Å, b = 18.6407(9) Å, c = 6.5871(6) Å, V = 2246 ų, Z = 8. A refinement of high-order diffraction data yielded reliability factors of R(F) = 0.9%, R _w(F) = 0.8%, GoF = 1.40 for 1856 high-angle reflections (0.7 ?in θ/λ <0.903 Å^?1) and R(F) = 1.0%, R _W(F) = 1.2%, GoF = 3.07 for all 3471 independent reflections in the complete sphere of reflection. The T-O distances as well as the T-O-T angles were found to be strongly influenced by the different bond strengths received by the individual oxygen atoms. The T O distances calculated using Baur's extended valence rule agree on average within 0.003 Å with the observed values. An analysis of the mean square displacement amplitudes allowed a separation of the external and internal vibrational amplitudes along the T-O bonds as well as along the Na O and H₂O-O bond directions and the calculation of force constants. The internal vibrational amplitudes (ΔU) of the T-O vibrations are in the range of 5 to 11 × 10^-4 Ų, that is about one order of magnitude smaller than the mean square displacement amplitudes of the external vibrations. The corresponding force constants are F = 354 to 824 Nm^?1. The values of the force constants of the motion of the Na-ion and the water molecule against the framework oxygen atoms lie in the range between F = 57 and 293 Nm^?1. This is the first instance where displacement amplitudes from a zeolite structure refinement could be apportioned between contributions from internal and external vibrations for individual bonds.     

Oxygen isotope measurements of individual unmelted Antarctic micrometeorites

We present oxygen isotope measurements of 28 unmelted Antarctic micrometeorites measuring 150-250 μm (long axis) collected in the South Pole water well. The micrometeorites were all unmelted and classified as either fine-grained, scoriaceous, coarse-grained or composite (a mix of two other classes). Spot analyses were made of each micrometeorite type using an ion microprobe. The oxygen isotope values were measured relative to standard mean ocean water (SMOW) and range from δ¹⁸O = 3‰ to 60‰ and δ¹⁷O = −1‰ to 32‰, falling along the terrestrial fractionation line (TFL) within 2σ errors. Several analytical spots (comprising multiple phases) were made on each particle. Variability in the oxygen isotope ratios was observed among micrometeorite types, between micrometeorites of the same type and between analytical spots on a single micrometeorite indicating that micrometeorites are isotopically heterogeneous. In general, the lowest isotope values are associated with the coarse-grained micrometeorites whereas most of the fine-grained and scoriaceous micrometeorites have an average δ¹⁸O ? 22‰, suggesting that the matrix in micrometeorites is isotopically heavier than the anhydrous silicate phases. The oxygen isotope values for the coarse-grained micrometeorites, composed mainly of anhydrous phases, do not lie along the carbonaceous chondrite anhydrous mineral (CCAM) line, as observed for olivines, pyroxenes and some kinds of chondrules in carbonaceous chondrites, suggesting that coarse-grained MMs are not related to chondrules, as previously thought. Our measurements span the same range as values found for melted micrometeorites in other studies. Although four of the micrometeorites have oxygen isotope values lying along the TFL, close to the region where the bulk CI carbonaceous chondrites are found, 21 particles have very enriched ¹⁷O and ¹⁸O values that have not been reported in previous analyses of chondrite matrix material, suggesting that they could be a new type of Solar System object. The parent bodies of the micrometeorites with higher ¹⁸O values may be thermal metamorphosed carbonaceous asteroids that have not been found as meteorites either because they are friable asteroids that produce small particles rather than rocks upon collision with other bodies, or because the rocks they produce are too friable to survive atmospheric entry.     

An internally consistent set of thermodynamic data for twentyone CaO-Al2O3-SiO2- H2O phases by linear parametric programming

The technique of linear parametric programming has been applied to derive sets of internally consistent thermodynamic data for 21 condensed phases of the quaternary system CaO-Al₂O₃-SiO₂-H₂O (CASH) (Table 4). This was achieved by simultaneously processing:

1. calorimetric data for 16 of these phases (Table 1), and
2. experimental phase equilibria reversal brackets for 27 reactions (Table 3) involving these phases.

Calculation of equilibrium P-T curves of several arbitrarily picked reactions employing the preferred set of internally consistent thermodynamic data from Table 4 shows that the input brackets are invariably satisfied by the calculations (Fig. 2a). By contrast, the same equilibria calculated on the basis of a set of thermodynamic data derived by applying statistical methods to a large body of comparable input data (Haas et al. 1981; Hemingway et al. 1982) do not necessarily agree with the experimental reversal brackets. Prediction of some experimentally investigated phase relations not included into the linear programming input database also appears to be remarkably successful. Indications are, therefore, that the thermodynamic data listed in Table 4 may be used with confidence to predict geologic phase relations in the CASH system with considerable accuracy. For such calculated phase diagrams and their petrological implications, the reader's attention is drawn to the paper by Chatterjee et al. (1984).     

Laboratory simulation of meteoritic noble gases. I. Sorption of xenon on carbon: Trapping experiments

We have studied trapping of radioactive ¹²⁷Xe in three types of carbon: carbon black (lamp black  LB), pyrolyzed polyvinylidene chloride (PVDC), and pyrolyzed acridine (C₁₃H₉N). A total of 86 samples were exposed to Xe at T between 100 and 1000°C, for times between 5 min and 240 hours, at p_xe ~ 5 × 10^?7 atm. Excess gas phase and loosely sorbed Xe were pumped away and the remaining, tightly bound Xe was measured by γ-spectrometry.At 100°C,× >90% of the Xe desorbs within a few minutes' pumping but a small amount remains even after 4000 min. Distribution coefficients for this tightly bound Xe are ~1 × 10^?2, 1 and 10 ccSTP/g atm for LB, acridine and PVDC carbons. The tightly bound Xe consists of two components. One occurs over the entire range 100–1000°C, becoming less abundant at high T; it appears to be physisorbed. The other occurs only at T > 500°C and is probably due to volume diffusion. The adsorbed component in LB has an apparent ΔH between ?2.3 and ?5.7 kcal/mole. The diffused component, which occurs in LB and possibly in acridine carbon, has an activation energy Q = 27 ± 8 kcal/mole and a diffusion coefficient D = 1.3 × 10^?17 cm²/sec at 1000°C. These values are comparable to those found for other types of amorphous carbon (Morrisonet al., 1963; Nakai et al., 1960).The low-T component displays two paradoxical features: low ΔH_ads, in the range for Xe physisorbed on carbon, but exceedingly long adsorption or desorption times (~10³ min at 100–400 or 1000°C). Although these long times seem to suggest a high energy process such as chemisorption, our results are best explained by a model that invokes physisorption within a labyrinth of micropores—of atomic dimensions—known to exist in amorphous carbons. The long adsorption/desorption times reflect either the long distances (~5 cm) Xe atoms must migrate by random walk to enter or leave the labyrinth, or the long times needed for Xe atoms to traverse tight spots or constricted pores that connect interior and exterior surfaces of the carbon (activated entry). Both variants of this model predict long equilibration times for the observed ΔH_ads of ?2 to ?6 kcal/mole. Apparently, xenon can be tightly trapped in carbon without resorting to high-energy bonding or to exotic mechanisms.These results suggest that “planetary” type noble gases in meteorites, located at or near grain surfaces of amorphous carbon, may be trapped by adsorption in micropores, whereas components such as CCFXe, which are uniformly distributed in their carrier phases, may be trapped by mechanisms such as volume diffusion or ion implantation.