Stochastic analysis of bed-thickness distribution of sediments

On formation of a bed and distribution of bed thickness, A. N. Kolmogorov presented a mathematical explanation that if repetitive alternations of material accumulation and erosion form a sequence of beds, the resultant bed-thickness distribution curve takes a shape truncated by the ordinate at zero thickness. In this truncated distribution curve, its continuation and extension from positive to negative thickness represents the distribution of beds with negative thickness, that is, the depth of erosion. When a distribution curve, including both positive and negative parts, is expressed by a function f(x),the ratio \(\int_0^\infty {f(x)dx to} \int_{ - \infty }^\infty {f(x)dx} \) ,called Kolmogorov's coefficient and designated as p,is a parameter representing the degree of accumulation in the depositional environment. On the assumption that f(x)is described by the Gaussian distribution function, the coefficient pfor Permian and Pliocene sequences in central Japan was calculated. The coefficients also were obtained from published data for different types of sediments from other areas. It was determined that they are more or less different depending on their depositional environments. The calculated results are summarized as follows: $$\begin{gathered} p = 0.80 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ p = 0.65 - 0.95for{\text{ }}nearshore{\text{ }}sediments \hfill \\ p = 0.55 - 0.95for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ p = 0.90 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ In addition, a ratio \(q = \int_0^\infty {xf(x)dx/} \int_{ - \infty }^\infty {|x|f(x)dx} \) ,called Kolmogorov's ratio in this paper, is introduced for estimating a degree of total thickness actually observed in the field relative to total thickness once present in a basin. The calculated results of Kolmogorov's ratio are as follows: $$\begin{gathered} q = 0.88 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ q = 0.68 - 0.98for{\text{ }}nearshore{\text{ }}sediments \hfill \\ q = 0.55 - 0.96for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ q = 0.92 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ The sedimentological significance of these values is discussed.     

Experimental determination of the reaction paragonite=jadeite+kyanite+H2O,and internally consistent thermodynamic data for part of the system Na2O−Al2O3−SiO2−H2O,with applications to eclogites and blueschists

Hydrothermal reversal experiments have been performed on the upper pressure stability of paragonite in the temperature range 550–740 ° C. The reaction $$\begin{gathered} {\text{NaAl}}_{\text{3}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{1 0}}} ({\text{OH)}}_{\text{2}} \hfill \\ {\text{ paragonite}} \hfill \\ {\text{ = NaAlSi}}_{\text{2}} {\text{O}}_{\text{6}} + {\text{Al}}_{\text{2}} {\text{SiO}}_{\text{5}} + {\text{H}}_{\text{2}} {\text{O}} \hfill \\ {\text{ jadeite kyanite vapour}} \hfill \\ \end{gathered}$$ has been bracketed at 550 ° C, 600 ° C, 650 ° C, and 700 ° C, at pressures 24–26 kb, 24–25.5 kb, 24–25 kb, and 23–24.5 kb respectively. The reaction has a shallow negative slope (? 10 bar °C^?1) and is of geobarometric significance to the stability of the eclogite assemblage, omphacite+kyanite. The experimental brackets are thermodynamically consistent with the lower pressure reversals of Chatterjee (1970, 1972), and a set of thermodynamic data is presented which satisfies all the reversal brackets for six reactions in the system Na₂O-Al₂O₃-SiO₂-H₂O. The Modified Redlich Kwong equation for H₂O (Holloway, 1977) predicts fugacities which are too high to satisfy the reversals of this study. The P-T stabilities of important eclogite and blueschist assemblages involving omphacite, kyanite, lawsonite, Jadeite, albite, chloritoid, and almandine with paragonite have been calculated using thermodynamic data derived from this study.     

The reversed alumina contents of orthopyroxene in equilibrium with spinel and forsterite in the system MgO-Al2O3-SiO2

Equilibrium alumina contents of orthopyroxene coexisting with spinel and forsterite in the system MgO-Al₂O₃-SiO₂ have been reversed at 15 different P-T conditions, in the range 1,030–1,600° C and 10–28 kbar. The present data and three reversals of Danckwerth and Newton (1978) have been modeled assuming an ideal pyroxene solid solution with components Mg₂Si₂O₆ (En) and MgAl₂SiO₆ (MgTs), to yield the following equilibrium condition (J, bar, K): $$\begin{gathered} RT{\text{ln(}}X_{{\text{MgTs}}} {\text{/}}X_{{\text{En}}} {\text{) + 29,190}} - {\text{13}}{\text{.42 }}T + 0.18{\text{ }}T + 0.18{\text{ }}T^{1.5} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [0.013 + 3.34 \times 10^{ - 5} (T - 298) - 6.6 \times 10^{ - 7} P]P. \hfill \\ \end{gathered} $$ The data of Perkins et al. (1981) for the equilibrium of orthopyroxene with pyrope have been similarly fitted with the result: $$\begin{gathered} - RT{\text{ln(}}X_{{\text{MgTs}}} \cdot X_{{\text{En}}} {\text{) + 5,510}} - 88.91{\text{ }}T + 19{\text{ }}T^{1.2} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [ - 0.832 - 8.78{\text{ }} \times {\text{ 10}}^{ - {\text{5}}} (T - 298) + 16.6{\text{ }} \times {\text{ 10}}^{ - 7} P]{\text{ }}P. \hfill \\ \end{gathered} $$ The new parameters are in excellent agreement with measured thermochemical data and give the following properties of the Mg-Tschermak endmember: $$H_{f,970}^0 = - 4.77{\text{ kJ/mol, }}S_{298}^0 = 129.44{\text{ J/mol}} \cdot {\text{K,}}$$ and $$V_{298,1}^0 = 58.88{\text{ cm}}^{\text{3}} .$$ The assemblage orthopyroxene+spinel+olivine can be used as a geothermometer for spinel lherzolites, subject to a choice of thermodynamic mixing models for multicomponent orthopyroxene and spinel. An ideal two-site mixing model for pyroxene and Sack's (1982) expressions for spinel activities provide, with the present experimental calibration, a geothermometer which yields temperatures of 800° C to 1,350° C for various alpine peridotites and 850° C to 1,130° C for various volcanic inclusions of upper mantle origin.     

The Reaction Titanite+Kyanite=Anorthite+Rutile and Titanite-Rutile Barometry in Eclogites

Titanite and rutile are a common mineral pair in eclogites, and many equilibria involving these phases are potentially useful in estimating pressures of metamorphism. We have reversed one such reaction,

H<Subscript>2</Subscript>O diffusion in peralkaline to peraluminous rhyolitic melts

The diffusion of water in a peralkaline and a peraluminous rhyolitic melt was investigated at temperatures of 714–1,493 K and pressures of 100 and 500 MPa. At temperatures below 923 K dehydration experiments were performed on glasses containing about 2 wt% H₂O _t in cold seal pressure vessels. At high temperatures diffusion couples of water-poor (<0.5 wt% H₂O _t ) and water-rich (~2 wt% H₂O _t ) melts were run in an internally heated gas pressure vessel. Argon was the pressure medium in both cases. Concentration profiles of hydrous species (OH groups and H₂O molecules) were measured along the diffusion direction using near-infrared (NIR) microspectroscopy. The bulk water diffusivity () was derived from profiles of total water () using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between and Both methods consistently indicate that is proportional to in this range of water contents for both bulk compositions, in agreement with previous work on metaluminous rhyolite. The water diffusivity in the peraluminous melts agrees very well with data for metaluminous rhyolites implying that an excess of Al₂O₃ with respect to alkalis does not affect water diffusion. On the other hand, water diffusion is faster by roughly a factor of two in the peralkaline melt compared to the metaluminous melt. The following expression for the water diffusivity in the peralkaline rhyolite as a function of temperature and pressure was obtained by least-squares fitting:

The origin of garnet in the anorthosite-charnockite suite of the Adirondacks

Detailed analysis of textural and chemical criteria in rocks of the anorthosite-charnockite suite of the Adirondack Highlands suggests that development of garnet in silica-saturated rocks of the suite occurs according to the reaction: $$\begin{gathered} {\text{Anorthite}} {\text{Orthopyroxene}} {\text{Quartz}} \hfill \\ {\text{2CaAl}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} + (6 - \alpha )({\text{Fe,Mg}}){\text{SiO}}_{\text{3}} + \alpha {\text{Fe - Oxide + (}}\alpha {\text{ - 2)SiO}}_{\text{2}} \hfill \\ {\text{Garnet}} {\text{Clinopyroxene}} \hfill \\ = {\text{Ca(Fe,Mg)}}_{\text{5}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{6}} {\text{O}}_{{\text{24}}} + {\text{Ca(Fe,Mg)Si}}_{\text{2}} {\text{O}}_{\text{6}} \hfill \\ \end{gathered} $$ , where α is a function of the distribution of Fe and Mg between the several coexisting ferromagnesian phases. Depending upon the relative amounts of Fe and Mg present, quartz may be either a reactant or a product. Using an aluminum-fixed reference frame, this reaction can be restated in terms of a set of balanced partial reactions describing the processes occurring in spatially separated domains within the rock. The fact that garnet invariably replaces plagioclase as opposed to the other reactant phases indicates that the aluminum-fixed model is valid as a first approximation. This reaction is univariant and produces unzoned garnet. It differs from a similar equation proposed by de Waard (1965) for the origin of garnet in Adirondack metabasic rocks, i.e. 6 Orthopyroxene+2 Anorthite = Clinopyroxene+Garnet+2 Quartz, the principle difference being that iron oxides (ilmenite and/or magnetite) are essential reactant phases in the present reactions. The product assemblage (garnet+clinopyroxene+plagioclase ± orthopyroxene ± quartz) is characteristic of the clinopyroxene-almandine subfacies of the granulite facies.     

Reaction spaces and P-T paths: from amphibole eclogite to greenschist facies in the Austroalpine domain (Oetztal Complex)

The transition from feldspar amphibolite to eclogite is a very wide P-T field that extends from some-where close to 5 kbar where the garnet-amphibole pair starts to appear, to 10–20 kbar at albite-out reaction, then up to 25–30 kbar where an hydrated phase such as amphibole can be stable with pyroxene and garnet. Thus the assemblage garnet (py)+ amphibole (tr)+epidote (cz)±plagioclase (ab)±clinopyroxene (di)±quartz (qz)±fluid is commonly reported in a large number of metamorphic terrains. These mineral phases are complex solid-solutions which adapt to variations in environmental conditions mainly by means of continuous reactions. The reaction space, introduced by. Thompson in 1982a, provides a very elegant and powerful tool to approach these high-variance assemblages. The reactions:

Kinetics of heterogeneous bubble nucleation in rhyolitic melts: implications for the number density of bubbles in volcanic conduits and for pumice textures

We performed decompression experiments to simulate the ascent of a phenocryst-bearing rhyolitic magma in a volcanic conduit. The starting materials were bubble-free rhyolites water-saturated at 200 MPa–800°C under oxidizing conditions: they contained 6.0 wt% dissolved H₂O and a dense population of hematite crystals (8.7 ± 2 × 10⁵ mm⁻³). Pressure was decreased from the saturation value to a final value ranging from 99 to 20 MPa, at constant temperature (800°C); the rate of decompression was either 1,000 or 27.8 kPa/s. In all experiments, we observed a single event of heterogeneous bubble nucleation beginning at a pressure P _N equal to 63 ± 3 MPa in the 1,000 kPa/s series, and to 69 ± 1 MPa in the 27.8 kPa/s series. Below P _N, the degree of water supersaturation in the liquid rapidly decreased to a few 0.1 wt%, the nucleation rate dropped, and the bubble number density (BND) stabilized to a value strongly sensitive to decompression rate: 80 mm⁻³ at 27.8 kPa/s vs. 5,900 mm⁻³ at 1,000 kPa/s. This behaviour is like the behavior formerly described in the case of homogeneous bubble nucleation in the rhyolite-H₂O system and in numerical simulations of vesiculation in ascending magmas. Similar degrees of water supersaturation were measured at 27.8 and 1,000 kPa/s, implying that a faster decompression rate does not result in a larger departure from equilibrium. Our experimental results imply that BNDs in acid to intermediate magmas ascending in volcanic conduits will depend on both the decompression rate and the number density of phenocrysts, especially the number density of magnetite microphenocrysts (1–100 mm⁻³), which is the only mineral species able to reduce significantly the degree of water supersaturation required for bubble nucleation. Very low BNDs (≈1 mm⁻³) are predicted in the case of effusive eruptions ( ≈ 0.1 kPa/s). High BNDs (up to 10⁷ mm⁻³) and bimodal bubble size distributions are expected in the case of explosive eruptions: (1) a relatively small number density of bubbles (1–100 mm⁻³) will first nucleate in the lower part of the conduit ( ≈ 10 kPa/s), either at high pressure on magnetite or at lower pressure on quartz and feldspar (or by homogeneous nucleation in the liquid) and (2) then, extreme decompression rates near the fragmentation level ( ≈ 10³ kPa/s) will trigger a major nucleation event leading to the multitude of small bubbles, typically a few micrometers to a few tens of micrometers in diameter, which characterizes most silicic pumices.     

Reversed experimental calibration of the garnet-clinopyroxene Fe — Mg exchange thermometer

The temperature-sensitive Fe,Mg exchange equilibrium,

(Ca,Mg)2Si2O6 clinopyroxenes: A solution model based on nonconvergent site-disorder

Experiments at high pressure and temperature indicate that excess Ca may be dissolved in diopside. If the (Ca, Mg)₂Si₂O₆ clinopyroxene solution extends to more Ca-rich compositions than CaMgSi₂O₆, macroscopic regular solution models cannot strictly be applied to this system. A nonconvergent site-disorder model, such as that proposed by Thompson (1969, 1970), may be more appropriate. We have modified Thompson's model to include asymmetric excess parameters and have used a linear least-squares technique to fit the available experimental data for Ca-Mg orthopyroxene-clinopyroxene equilibria and Fe-free pigeonite stability to this model. The model expressions for equilibrium conditions \(\mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction A) and \(\mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction B) are given by: 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Mg}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ W_{21} [2(X_{{\text{Ca}}}^{{\text{M2}}} )^3 - (X_{{\text{Ca}}}^{{\text{M2}}} ] \hfill \\ {\text{ + 2W}}_{{\text{22}}} [X_{{\text{Ca}}}^{{\text{M2}}} )^2 - (X_{{\text{Ca}}}^{{\text{M2}}} )^3 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{Wo}}}^{{\text{opx}}} )^2 \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Ca}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ 2W_{21} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^2 - (X_{{\text{Mg}}}^{{\text{M2}}} )^3 ] \hfill \\ {\text{ + W}}_{{\text{22}}} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^3 - (X_{{\text{Mg}}}^{{\text{M2}}} )^2 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{En}}}^{{\text{opx}}} )^2 \hfill \\ \hfill \\ \end{gathered} $$ where 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = 2.953 + 0.0602{\text{P}} - 0.00179{\text{T}} \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = 24.64 + 0.958{\text{P}} - (0.0286){\text{T}} \hfill \\ {\text{W}}_{{\text{21}}} = 47.12 + 0.273{\text{P}} \hfill \\ {\text{W}}_{{\text{22}}} = 66.11 + ( - 0.249){\text{P}} \hfill \\ {\text{W}}^{{\text{opx}}} = 40 \hfill \\ \Delta {\text{G}}_*^0 = 155{\text{ (all values are in kJ/gfw)}}{\text{.}} \hfill \\ \end{gathered} $$ . Site occupancies in clinopyroxene were determined from the internal equilibrium condition 1 $$\begin{gathered} \Delta G_{\text{E}}^{\text{O}} = - {\text{RT 1n}}\left[ {\frac{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}{{X_{{\text{Ca}}}^{{\text{M2}}} \cdot X_{{\text{Mg}}}^{{\text{M1}}} }}} \right] + \tfrac{1}{2}[(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} )(2{\text{X}}_{{\text{Ca}}}^{{\text{M2}}} - 1) \hfill \\ {\text{ + }}\Delta G_*^0 (X_{{\text{Ca}}}^{{\text{M1}}} - X_{{\text{Ca}}}^{{\text{M2}}} ) + \tfrac{3}{2}(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} ) \hfill \\ {\text{ (1}} - 2X_{{\text{Ca}}}^{{\text{M1}}} )(X_{{\text{Ca}}}^{{\text{M1}}} + \tfrac{1}{2})] \hfill \\ \end{gathered} $$ where δG _E ⁰ =153+0.023T+1.2P. The predicted concentrations of Ca on the clinopyroxene Ml site are low enough to be compatible with crystallographic studies. Temperatures calculated from the model for coexisting ortho- and clinopyroxene pairs fit the experimental data to within 10° in most cases; the worst discrepancy is 30°. Phase relations for clinopyroxene, orthopyroxene and pigeonite are successfully described by this model at temperatures up to 1,600° C and pressures from 0.001 to 40 kbar. Predicted enthalpies of solution agree well with the calorimetric measurements of Newton et al. (1979). The nonconvergent site disorder model affords good approximations to both the free energy and enthalpy of clinopyroxenes, and, therefore, the configurational entropy as well. This approach may provide an example for Febearing pyroxenes in which cation site exchange has an even more profound effect on the thermodynamic properties.     

Crystal field and covalency of octahedral chromium in natural [8](Mg1?xCax) 3[6] (Al0.67Cr0.33)2Si3O12 garnets from upper mantle rocks

In the course of a thorough study of the influences of the second coordination sphere on the crystal field parameters of the 3d ^N -ions and the character of 3d ^N –O bonds in oxygen based minerals, 19 natural Cr³⁺-bearing (Mg,Ca)-garnets from upper mantle rocks were analysed and studied by electronic absorption spectroscopy, EAS. The garnets had compositions with populations of the ^[8] X-sites by 0.881 ± 0.053 (Ca + Mg) and changing Ca-fractions in the range 0.020 ≤ w _Ca[8] ≤ 0.745, while the ^[6] Y-site fraction was constant with x _Cr³⁺ _[6] = 0.335 ± 0.023. The garnets had colours from deeply violet-red for low Ca-contents (up to x _Ca = 0.28), grey with 0.28 ≤ x _Ca ≤ 0.4 and green with 0.4 ≤ x _Ca. The crystal field parameter of octahedral Cr³⁺ 10Dq decreases strongly on increasing Ca-fraction from 17,850 cm⁻¹ at x _Ca[8] = 0.020 to 16,580 cm⁻¹ at x _Ca[8] = 0.745. The data could be fit with two model which do statistically not differ: (1) two linear functions with a discontinuity close to x _Ca[8] ≈ 0.3,

Re-evaluation of lead isotopic data,southern Massif Central,France

Three independent Pb isotope homogenizing processes operating on large volumes of rock material during limited intervals in the Phanerozoic have been used to define a unique evolutionary curve for rock and ore lead isotopic compositions of the southern Massif Central, France. The model is

The Hydrolysis of Al(III) in NaCl solutions: A Model for M(II), M(III), and M(IV) Ions

Understanding the identity and stability of the hydrolysis products of metals is required in order to predict their behavior in natural aquatic systems. Despite this need, the hydrolysis constants of many metals are only known over a limited range of temperature and ionic strengths. In this paper, we show that the hydrolysis constants of 31 metals [i.e. Mn(II), Cr(III), U(IV), Pu(IV)] are nearly linearly related to the values for Al(III) over a wide range of temperatures and ionic strengths. These linear correlations allow one to make reasonable estimates for the hydrolysis constants of +2, +3, and +4 metals from 0 to 300°C in dilute solutions and 0 to 100°C to 5 m in NaCl solutions. These correlations in pure water are related to the differences between the free energies of the free ion and complexes being almost equal $$ \Updelta {\text{G}}^\circ \left( {{\text{Al}}^{3 + } } \right) - \Updelta {\text{G}}^\circ \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{{\left( {3 - j} \right)}} } \right) \cong \Updelta {\text{G}}^\circ \left( {{\text{M}}^{n + } } \right) - \Updelta {\text{G}}^\circ \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{{\left( {n - j} \right)}} } \right) $$ The correlation at higher temperatures is a result of a similar relationship between the enthalpies of the free ions and complexes $$ \Updelta {\text{H}}^\circ \left( {{\text{Al}}^{3 + } } \right) - \Updelta {\text{H}}^\circ \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{3 - j} } \right) \cong \Updelta {\text{H}}^\circ \left( {{\text{M}}^{n + } } \right) - \Updelta {\text{H}}^\circ \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{n - j} } \right) $$ The correlations at higher ionic strengths are the result of the ratio of the activity coefficients for Al(III) being almost equal to that of the metal. $$ \gamma \left( {{\text{M}}^{n + } } \right)/\gamma \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{n - j} } \right) \cong \gamma \left( {{\text{Al}}^{3 + } } \right)/\gamma \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{3 - j} } \right) $$ The results of this study should be useful in examining the speciation of metals as a function of pH in natural waters (e.g. hydrothermal fresh waters and NaCl brines).     

Temperature coefficients of elastic constants of single crystal MgO between 80 and 1,300 K

Elastic constants of single crystal MgO have been measured by the rectangular parallelepiped resonance (RPR) method at temperatures between 80 and 1,300 K. Elastic constants C _ij (Mbar=10³ kbar) and their temperature coefficients (kbar/K) are: $$\begin{gathered} {\text{ }}C_{{\text{11}}} {\text{ }}C_{{\text{12}}} {\text{ }}C_{{\text{44}}} {\text{ }}K_s {\text{ }}C_s \hfill \\ C_{ij} {\text{ 300 K 2}}{\text{.966 0}}{\text{.959 1}}{\text{.562 1}}{\text{.628 1}}{\text{.004}} \hfill \\ \partial C_{ij} {\text{/}}\partial T{\text{100 K }} - {\text{0}}{\text{.259 0}}{\text{.013 }} - {\text{0}}{\text{.072 }} - {\text{0}}{\text{.078 }} - {\text{0}}{\text{.136}} \hfill \\ {\text{ 300K }} - {\text{0}}{\text{.596 0}}{\text{.068 }} - {\text{0}}{\text{.122 }} - {\text{0}}{\text{.153 }} - {\text{0}}{\text{.332}} \hfill \\ {\text{ 800 K }} - {\text{0}}{\text{.619 0}}{\text{.009 }} - {\text{0}}{\text{.152 }} - {\text{0}}{\text{.200 }} - {\text{0}}{\text{.314}} \hfill \\ {\text{ 1,300 K }} - {\text{0}}{\text{.598 0}}{\text{.036 }} - {\text{0}}{\text{.130 }} - {\text{0}}{\text{.223 }} - {\text{0}}{\text{.218}} \hfill \\ \end{gathered} $$ By combining the present results with the previous data on the thermal expansivity and specific heat, the thermodynamic properties of magnesium oxide are presented and discussed. The elastic parameters of MgO at very high temperatures in the earth's lower mantle are also clarified.     

Opening and resetting temperatures in heating geochronological systems

We present a theoretical model for diffusive daughter isotope loss in radiochronological systems with increasing temperature. It complements previous thermochronological models, which focused on cooling, and allows for testing opening and resetting of radiochronometers during heating. The opening and resetting temperatures are, respectively,

Experimental study of Fe-Mg distribution between biotite and orthopyroxene at P=490 MPa

Partitioning of Mg and Fe between coexisting biotite and orthopyroxene has been experimentally determined at temperatures 700, 750 and 800° C and 490 MPa total pressure in the system KAlO₂-MgO-FeO-SiO₂-H₂O. Oxygen fugacity was controlled by the QFM buffer. Starting materials were synthetic minerals of differing Fe/(Fe+Mg) values. Run products were analyzed for partitioning of components by a microprobe. Orthopyroxene was established to be notably inhomogeneous, whereas biotite was essentially homogeneous. To establish equilibrium relations, statistical treatment of the results of each experiment in addition to the whole complex of experimental data was applied. The regression equations for isotherms of the Fe-Mg partitioning between the minerals studied have been obtained. As a result, the equation for a two-dimensional regression may be written as: $$\begin{gathered} Y = (A + A_1 t + A_2 t^2 )(X - X^4 ) + (B + B_1 t + B_1 t^2 )(X^2 - X^4 ) + \hfill \\ (C + C_1 t + C_1 t^2 )(X^3 - X^4 ) + X^4 {\text{ where }}Y = X_{{\text{Opx}}}^{{\text{Fe}}} ;{\text{ X}} = {\text{X}}_{{\text{Bi}}}^{{\text{Fe}}} ; \hfill \\ t = 1000/T,K, \hfill \\ \begin{array}{*{20}c} {A = {\text{ }}4.59398,} & {A_1 = - {\text{ }}8.29838,} & {A_2 = {\text{ }}4.97316,} \\ {B = - 11.13731,} & {B_1 = {\text{ }}28.19304,} & {B_2 = - 20.98240,} \\ {A = {\text{ }}8.25072,} & {C_1 = - 20.80485,} & {C_2 = {\text{ }}15.35967} \\ \end{array} \hfill \\ {\text{ }}\sigma = 0.0143{\text{ }} \hfill \\ \end{gathered}$$ . This equation enables extrapolation of partitioning isotherms over a wide range of temperatures.     

Thermodynamic Properties of Brucite Determined by Solubility Studies and Their Significance to Nuclear Waste Isolation

Solubility experiments were conducted for the dissolution reaction of brucite, Mg(OH)₂ (cr): Experiments were conducted from undersaturation in deionized (DI) water and 0.010–4.4 m NaCl solutions at 22.5°C. In addition, brucite solubility was measured from supersaturation in an experiment in which brucite was precipitated via dropwise addition of 0.10 m NaOH into a 0.10 m MgCl₂ solution also at 22.5°C. The attainment of the reversal in equilibrium was demonstrated in this study. The solubility constant at 22.5°C at infinite dilution calculated from the experimental results from the direction of supersaturation by using the specific interaction theory (SIT) is: with a corresponding value of 17.0 ± 0.2 (2σ) when extrapolated to 25°C. The dimensionless standard chemical potential (μ°/RT) of brucite derived from the solubility data in 0.010 m to 4.4 m NaCl solutions from undersaturation extrapolated to 25°C is −335.76 ± 0.45 (2σ), with the corresponding Gibbs free energy of formation of brucite, , being −832.3 ± 1.1 (2σ) kJ mol⁻¹. In combination with the auxiliary thermodynamic data, the is calculated to be 17.1 ± 0.2 (2σ), based on the above Gibbs free energy of formation for brucite. This study recommends an average value of 17.05 ± 0.2 in logarithmic unit as solubility constant of brucite at 25°C, according to the values from both supersaturation and undersaturation. Sandia National Laboratories is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy’s National Nuclear Security Administration under Contract DE-AC04-94AL85000.     

The Equation of State for Caspian Sea Waters

The density ρ of Caspian Sea waters was measured as a function of temperature (273.15–343.15) K at conductivity salinities of 7.8 and 11.3 using the Anton-Paar Densitometer. Measurements were also made on one of the samples (S = 11.38) diluted with water as a function of temperature (T = 273.15–338.15 K) and salinity (2.5–11.3). These latter results have been used to develop an equation of state for the Caspian Sea (σ = ±0.007 kg m⁻³)

The combined effect of F and H2O on interdiffusion between peralkaline dacitic and rhyolitic melts

The effective binary diffusion coefficient (EBDC) of silicon has been measured during the interdiffusion of peralkaline, fluorine-bearing (1.3 wt% F), hydrous (3.3 and 6 wt% H₂O), dacitic and rhyolitic melts at 1.0 GPa and temperatures between 1100°C and 1400°C. From Boltzmann-Matano analysis of diffusion profiles the diffusivity of silicon at 68 wt% SiO₂ can be described by the following Arrhenius equations (with standard errors): $$\begin{gathered} {\text{with 1}}{\text{.3 wt\% F and 3}}{\text{.3\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.66}} \times {\text{10}}^{ - {\text{9}}} } \\ { - {\text{1}}{\text{.86}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ {\text{with 1}}{\text{.3 wt\% F and 6}}{\text{.0\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.51}} \times {\text{10}}^{ - {\text{8}}} } \\ { - {\text{1}}{\text{.77}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ \end{gathered} $$ where D is in m²s^?1 and activation energies are in kJ/mol. Diffusivities measured at 64 and 72 wt% SiO₂ are only slightly different from those at 68 wt% SiO₂ and frequently all measurements are within error of each other. Silicon, aluminum, iron, magnesium, and calcium EBDCs were also calculated from diffusion profiles by error function inversion techniques assuming constant diffusivity. With one exception, silicon EBDCs calculated by error function techniques are within error of Boltzmann-Matano EBDCs. Average diffusivities of Fe, Mg, and Ca were within a factor of 2.5 of silicon diffusivities whereas Al diffusivities were approximately half those of silicon. Alkalies diffused much more rapidly than silicon and non-alkalies, however their diffusivities were not quantitatively determined. Low activation energies for silicon EBDCs result in rapid diffusion at magmatic temperatures. Assuming that water and fluorine exert similar effects on melt viscosity at high temperatures, the viscosity can be calculated and used in the Eyring equation used to determine diffusivities, typically to within a factor of three of those measured in this study. This correlation between viscosity and diffusivity can be inverted to calculate viscosities of fluorine- and water-bearing granitic melts at magmatic temperatures; these viscosities are orders of magnitude below those of hydrous granitic melts and result in more rapid and effective separation of granitic magmas from partially molten source rocks. Comparison of Arrhenius parameters for diffusion measured in this study with Arrhenius parameters determined for diffusion in similar compositions at the same pressure demonstrates simple relationships between Arrhenius parameters, activation energy-E_a, kJ/mol, pre-exponential factor-D_o, m²s^?1, and the volatile, X=F or OH^?, to oxygen, O, ratio of the melt {(X/X+O)}: $$\begin{gathered} {\text{E}}a = - {\text{1533\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }} + {\text{213}}{\text{.3}} \hfill \\ {\text{D}}_{\text{O}} = {\text{2}}{\text{.13}} \times {\text{10}}^{ - {\text{6}}} {\text{exp}}\left[ { - {\text{6}}{\text{.5\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }}} \right] \hfill \\ \end{gathered} $$ These relationships can be used to estimate diffusion in various melts of dacitic to rhyolitic composition containing both fluorine and water. Calculations for the contamination of rhyolitic melts by dacitic enclaves at 800°C and 700°C provide evidence for the virtual inevitability of diffusive contamination in hydrous and fluorine-bearing magmas if they undergo magma mixing of any form.