Reaction rim growth in the system MgO-Al2O3-SiO2 under uniaxial stress

We synthesize reaction rims between thermodynamically incompatible phases in the system MgO-Al₂O₃-SiO₂ applying uniaxial load using a creep apparatus. Synthesis experiments are done in the MgO-SiO₂ and in the MgO-Al₂O₃ subsystems at temperatures ranging from 1150 to 1350 °C imposing vertical stresses of 1.2 to 29 MPa at ambient pressure and under a constant flow of dry argon. Single crystals of synthetic and natural quartz and forsterite, synthetic periclase and synthetic corundum polycrystals are used as starting materials. We produce enstatite rims at forsterite-quartz contacts, enstatite-forsterite double rims at periclase-quartz contacts and spinel rims at periclase-corundum contacts. We find that rim growth under the “dry” conditions of our experiments is sluggish compared to what has been found previously in nominally “dry” piston cylinder experiments. We further observe that the nature of starting material, synthetic or natural, has a major influence on rim growth rates, where natural samples are more reactive than synthetic ones. At a given temperature the effect of stress variation is larger than what is anticipated from the modification of the thermodynamic driving force for reaction due to the storage of elastic strain energy in the reactant phases. We speculate that this may be due to modification of the physical properties of the polycrystals that constitute the reaction rims or by deformation under the imposed load. In our experiments rim growth is very sluggish at forsterite-quartz interfaces. Rim growth is more rapid at periclase-quartz contacts. The spinel rims that are produced at periclase-corundum interfaces show parabolic growth indicating that reaction rim growth is essentially diffusion controlled. From the analysis of time series done in the MgO-Al₂O₃ subsystem we derive effective diffusivities for the Al₂O₃ and the MgO components in a spinel polycrystal as ${\rm D}_{MgO} = 1.4 \pm 0.2 \cdot 10^{-15}$  m²/s and ${\rm D}_{Al_2O_3} = 3.7 \pm 0.6 \cdot 10^{-16}$  m²/s for T?=?1350 °C and a vertical stress of 2.9 MPa.     

Solid–liquid–air coupling in multiphase porous media

This paper addresses various issues concerning the modelling of solid–liquid–air coupling in multiphase porous media with an application to unsaturated soils. General considerations based on thermodynamics permit the derivation and discussion of the general form of field equations; two cases are considered: a three phase porous material with solid, liquid and gas, and a two phase porous material with solid, liquid and empty space. Emphasis is placed on the presentation of differences in the formulation and on the role of the gas phase. The finite element method is used for the discrete approximation of the partial differential equations governing the problem. The two formulations are then analysed with respect to a documented drainage experiment carried out by the authors. The merits and shortcomings of the two approaches are shown. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling karst denudation on a synthetic landscape

This contribution presents the results of a numerical study of karst denudation on limestone plateaux. The landscape evolution model used incorporates not only long-range fluvial and short-range hill-slope processes, but also large-scale chemical dissolution of limestone surfaces. The relative efficiencies of fluvial and chemical processes are of equal importance to the landscape evolution of a plateau dropping to sea level along an escarpment. While fluvial processes have an impact confined mostly to river channels, the karst denudation process is more uniform, removing material also from the plateau surface. The combined effect of both processes results in a landscape evolution almost twice as effective as the purely erosional evolution of an insoluble landscape.     

Lichenometry in the Cordillera Blanca, Peru: “Little Ice Age” moraine chronology

This paper is a comparison and compilation of lichenometric and geomorphic studies performed by two independent teams in the Cordillera Blanca, Peru, in 1996 and 2002 on 66 “Little Ice Age” moraines of 14 glaciers. Using eleven new control points, we recalibrated the initial rapid growth phase of the previously established Rhizocarpon subgenus Rhizocarpon growth curve. This curve was then used to estimate the age of “Little Ice Age” moraines. The time of deposition of the most prominent and numerous terminal and lateral moraines on the Pacific-facing side of the Cordillera Blanca (between AD 1590 and AD 1720) corresponds to the coldest and wettest phase in the tropical Andes as revealed by ice-core data. Less prominent advances occurred between AD 1780 and 1880.     

Continental back-arc basin origin of some ophiolites from the Eastern Desert of Egypt

Summary Geochemical and petrographical data of three ophiolitic pillow metavolcanic occurrences from the central Eastern Desert of Egypt are presented. The investigated rocks show a subalkaline, tholeiitic affinity. Chemical data indicate that the metavolcanics have transitional within-plate basalt to island-arc basalt features, which are characteristics of basalts formed in ensialic back-arc basins. The association of the investigated ophiolites with volcanoclastic metasedimentary rocks of marine to continental facies is a further confirmation of their ensialic evolution. This suggestion, along with the geochronologic, isotopic and crustal growth rate evidences, revives interest in models that involve contribution from a pre-Pan-African continental crust at least in the southern part of the Egyptian Shield. Mixing between a depleted mantle-derived magma and an enriched crustal melt, a process similar to AFC (assimilation and fractional crystallization), is suggested for the evolution of the investigated rocks. This study provides evidence for formation of some ophiolites in the Eastern Desert of Egypt in continental (ensialic) back arc basins.     

Validation of turbulence closure parameterisations for stably stratified flows using the PROVESS turbulence measurements in the North Sea

A number of parameterisations for the simulation of mixing processes in the thermocline are compared and tested against the microstructure data of the PROVESS campaigns, conducted in the northern part of the North Sea during the autumn of 1998. The transport term in the turbulent kinetic energy equation is parameterised via the introduction of a third stability function S_k for turbulent energy diffusion. The formulations are compared with a simpler scheme based upon limiting conditions for turbulence variables. Improved results are obtained with a new form of S_k. The best agreement is, however, found with the simpler limiting scheme. This is explained in terms of a turbulence length scale theory for stably stratified turbulence. In agreement with previous laboratory and ocean data it is found that the ratios of the Thorpe and Kolmogorov scales to the Ozmidov length scale approach critical limiting values in the thermocline. The first of these conditions is satisfied when limiting conditions are implemented into the scheme, providing the necessary minimum value for the dissipation rate, whereas the schemes without limiting conditions fail to produce this critical ratio. The basic reason for this failure is that the Thorpe scale is overestimated, which is shown to be connected to an even larger overprediction of the dissipation rate of temperature variance. To investigate the impact of non-resolved advective processes and salinity stratification on the turbulence predictions, additional numerical experiments were conducted using a simple scheme for data assimilation. The best agreement is found again with the limiting scheme, which is able to make reasonable predictions for the dissipation rate without knowing the detailed shape of the mean stratification profile. It is shown that advective transport due to tidally and wind-driven motions has a non-negligible impact on vertical mixing. This is seen in the data and the models by periodic enhancements of turbulent mixing inside the thermocline.     

Nanoscale porosity in SAFOD core samples (San Andreas Fault)

With transmission electron microscopy (TEM) we observed nanometer-sized pores in four ultracataclastic and fractured core samples recovered from different depths of the main bore hole of the San Andreas Fault Observatory at Depth (SAFOD). Cutting of foils with a focused ion beam technique (FIB) allowed identifying porosity down to the nm scale. Between 40 and 50% of all pores could be identified as in-situ pores without any damage related to sample preparation. The total porosity estimated from TEM micrographs (1–5%) is comparable to the connected fault rock porosity (2.8–6.7%) estimated by pressure-induced injection of mercury. Permeability estimates for cataclastic fault rocks are 10^? 21–10^? 19 m² and 10^? 17 m² for the fractured fault rock. Porosity and permeability are independent of sample depth. TEM images reveal that the porosity is intimately linked to fault rock composition and associated with deformation. The TEM-estimated porosity of the samples increases with increasing clay content. The highest porosity was estimated in the vicinity of an active fault trace. The largest pores with an equivalent radius > 200 nm occur around large quartz and feldspar grains or grain-fragments while the smallest pores (equivalent radius < 50 nm) are typically observed in the extremely fine-grained matrix (grain size < 1 μm). Based on pore morphology we distinguish different pore types varying with fault rock fabric and alteration. The pores were probably filled with formation water and/or hydrothermal fluids at elevated pore fluid pressure, preventing pore collapse. The pore geometry derived from TEM observations and BET (Brunauer, Emmett and Teller) gas adsorption/desorption hysteresis curves indicates pore blocking effects in the fine-grained matrix. Observations of isolated pores in TEM micrographs and high pore body to pore throat ratios inferred from mercury injection suggest elevated pore fluid pressure in the low permeability cataclasites, reducing shear strength of the fault.     

Late Quaternary lake-level changes of Lake El'gygytgyn, NE Siberia

Lake El'gygytgyn is situated in a 3.6 Ma old impact crater in northeastern Siberia. Presented here is a reconstruction of the Quaternary lake-level history as derived from sediment cores from the southern lake shelf. There, a cliff-like bench 10 m below the modern water level has been investigated. Deep-water sediments on the shelf indicate high lake levels during a warm Mid-Pleistocene period. One period with low lake level prior to Marine Oxygen Isotope Stage (MIS) 3 has been identified, followed by a period of high lake level (10 m above present). In the course of MIS 2 the lake level dropped to − 10 m. At the end of MIS 2 the bench was formed and coarse beach sedimentation occurred. Subsequently, the lake level rose rapidly to the Holocene level. Changes in water level are likely linked to climate variability. During relatively temperate periods the lake becomes free of ice in summer. Strong wave actions transport sediment parallel to the coast and towards the outlet, where the material tends to accumulate, resulting in lake level rise. During cold periods the perennial lake ice cover hampers any wave activity and pebble-transport, keeping the outlet open and causing the lake level to drop.     

The electric field gradient in natural iron-doped chrysoberyl Al<Subscript>2</Subscript>BeO<Subscript>4</Subscript> and sinhalite MgAlBO<Subscript>4</Subscript> single crystals

This work reports on the evaluation of the electric field gradient (EFG) in natural chrysoberyl Al₂BeO₄ and sinhalite MgAlBO₄ using two different procedures: (1) experimental, with single crystal Mössbauer spectroscopy (SCMBS) on the three principal sections of each sample and (2) a “fully quantitative” method with cluster molecular orbital calculations based on the density functional theory. Whereas the experimental and theoretical results for the EFG tensor are in quantitative agreement, the calculated isomer shifts and optical d–d-transitions exhibit systematic deviations from the measured values. These deviations indicate that the substitution of Al and Mg with iron should be accompanied by considerable local expansion of the coordination octahedra.