Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Axial convergence in a well-mixed estuary

Transverse secondary circulations involving surface convergence, observed in a well-mixed estuary in North Wales, are made visible by the collection of surface material along an axial line which extends continuously for many kilometres through the estuary. The circulation and axial convergence, however, are seen only during the flood phase of the tide and no similar behaviour has been observed during the ebb phase.Convergent circulations in the estuary are associated with small but steady transverse density gradients in the cross-section, produced by non-uniform advection of the longitudinal gradient through the channel. A diagnostic model, using measured mean distributions of cross-sectional density, indicates surface transverse velocities (～0.1 ms^?1) similar to those observed in the estuary. The model further predicts appreciable transverse divergent currents at a fractional depth of 0.75: a prediction which has been tested in the estuary using a vertical array of accurately resolving current direction indicators.     

Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite

Two distinct age estimates for eclogite-facies metamorphism in the Sanbagawa belt have been proposed: (i) c.  120–110 Ma based on a zircon SHRIMP age for the Western Iratsu unit and (ii) c.  88–89 Ma based on a garnet–omphacite Lu–Hf isochron age from the Seba and Kotsu eclogite units. Despite the contrasting estimates of formation ages, petrological studies suggest the formation conditions of the Western Iratsu unit are indistinguishable from those of the other two units—all ∼20 kbar and 600–650 °C. Studies of the associated geological structures suggest the Seba and Western Iratsu units are parts of a larger semi-continuous eclogite unit. A combination of geochronological and petrological studies for the Western Iratsu eclogite offers a resolution to this discrepancy in age estimates. New Lu–Hf dating for the Western Iratsu eclogite yields an age of 115.9 ± 0.5 Ma that is compatible with the zircon SHRIMP age. However, petrological studies show that there was significant garnet growth in the Western Iratsu eclogite before eclogite facies metamorphism, and the early core growth is associated with a strong concentration of Lu. Pre-eclogite facies garnet (Grt1) includes epidote–amphibolite facies parageneses equilibrated at 550–650 °C and ∼10 kbar, and this is overgrown by prograde eclogite facies garnet (Grt2). The Lu–Hf age of c.  116 Ma is strongly skewed to the isotopic composition of Grt1 and is interpreted to reflect the age of the pre-eclogite phase. The considerable time gap ( c.  27 Myr) between the two Lu–Hf ages suggests they may be related to separate tectonic events or distinct phases in the evolution of the Sanbagawa subduction zone.     

Pressure-induced incipient amorphization of α-quartz and transition to coesite in an eclogite from Antarctica: a first record and some consequences

Monocrystalline quartz inclusions in garnet and omphacite from various eclogite samples from the Lanterman Range (Northern Victoria Land, Antarctica) have been investigated by cathodoluminescence (CL), Raman spectroscopy and imaging, and in situ X‐ray (XR) microdiffraction using the synchrotron. A few inclusions, with a clear‐to‐opalescent lustre, show ‘anomalous’ Raman spectra characterized by weak α‐quartz modes, the broadening of the main α‐quartz peak at 465 cm^?1, and additional vibrations at 480–485, 520–523 and 608 cm^?1. CL and Raman imaging indicate that this ‘anomalous’α‐quartz occurs as relicts within ordinary α‐quartz, and that it was preserved in the internal parts of small quartz inclusions. XR diffraction circular patterns display irregular and broad α‐quartz spots, some of which show an anomalous d‐spacing tightening of ～2%. They also show some very weak, hazy clouds that have d‐spacing compatible with coesite but not with α‐quartz. Raman spectrometry and XR microdiffraction characterize the anomalies with respect to α‐quartz as (i) a pressure‐induced disordering and incipient amorphization, mainly revealed by the 480–485 and 608‐cm^?1 Raman bands, together with (ii) a lattice densification, evidenced by d‐spacing tightening; (iii) the cryptic development of coesite, 520–523 cm^?1 being the main Raman peak of coesite and (iv) Brazil micro‐twinning. This ‘anomalous’α‐quartz represents the first example of pressure‐induced incipient amorphization of a metastable phase in a crustal rock. This issue is really surprising because pressure‐induced amorphization of metastable α‐quartz, observed in impactites and known to occur between 15 and 32 GPa during ultrahigh‐pressure (UHP) experiments at room temperature, is in principle irrelevant under normal geological P–T conditions. A shock (due to a seism?) or a local overpressure at the inclusion scale (due to expansion mismatch between quartz and its host mineral) seem the only geological mechanisms that can produce such incipient amorphization in crustal rocks. This discovery throws new light on the modality of the quartz‐coesite transition and on the pressure regimes (non‐lithostatic v. lithostatic) during high‐pressure/UHP metamorphism. In particular, incipient amorphization of quartz could favour the quartz‐coesite transition, or allow the growth of metastable coesite, as already experimentally observed.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Multiple longshore sand bars in the upper Chesapeake Bay

In the upper Chesapeake Bay (Maryland, U.S.A.) field surveys were conducted at 18 multiple longshore sand bar sites. The multiple bar systems were found in water depths less than approximately 2 m (mean sea level), and exhibited mild bottom slopes of 0·0052 or less. The number of bars composing each system ranged from four to 17 and the spacing between the crests typically increased in the offshore direction, ranging from 12 to 70 m. Bar height also typically increased with distance offshore and ranged from 0·03 to 0·61 m. A grain size analysis of crest and trough sediment did not reveal any significant differences and the sediment was categorized as ‘fine sand’. A review of the literature data indicated that the Chesapeake Bay multiple bars possessed similar characteristics to those found in Gelding Bay (Baltic Sea); similarities in fetch, wave height and tidal range between the two bays may account for this finding. The surf-scaling parameter indicated that the multiple bar systems were extremely dissipative with regard to wave energy, and wave height appeared to be an important factor in controlling bar spacing and bar height. A multiple wave break point hypothesis was discussed as a possible mechanism for the formation of Chesapeake Bay multiple longshore bars, and limited observational evidence appeared to support such a mechanism.     

Ultraviolet reflectance properties of asteroids

An analysis of the UV spectra of 28 asteroids obtained with the Internal Ultraviolet Explorer (IUE) satellite is presented. The spectra lie within the range 2100–3200 Å. Our results are examined in terms of both asteroid classification and of current ideas concerning the surface mineralogy of asteroids. For all the asteroids examined, UV reflectivity declines approximately linearly toward shorter wavelengths. In general, the same taxonomic groups are seen in the UV as in the visible and IR, although there is some evidence for asteroids with anomalous UV properties and for UV subclasses within the S class. No mineral absorption features are reported of strength similar to the strongest features in the visible and IR regions, but a number of shallow absorptions do occur and may provide valuable information on the surface composition of many asteroids.     

Organotemplate silica deposition in Neoproterozoic deep-marine environments: evidence from the Penganga Group, Adilabad, India

Silica in bedded cherts interstratified with manganese carbonates in a deep‐water carbonate ramp succession of the Neoproterozoic Penganga Group, India, displays film‐like, spherical to rod‐shaped, and tubular branching microstructures. The microstructures resemble mineralized extracellular polysaccharides, biofilms and bacterial morphologies. The microstructures suggest silicification by nucleation of silica on organic‐templates or indirectly as sorbed species accumulating on organic templates. Given that similar microstructures have also been documented in Archean cherts it is proposed that organotemplates might have been an important sink for the deposition of silica in Precambrian deep‐water marine environments.     

Fluid production and isograd reactions at contacts of carbonate-rich and carbonate-poor layers during progressive metamorphism

ABSTRACT With increasing temperature during prograde metamorphism reactions will occur first at the lithological contacts of mixed pelite and calcsilicate terranes. At these interfaces, a fluid of lower chemical potential of H₂O and CO₂ than that required to produce a fluid in either layer can be produced whether reaction is caused by fluid infiltration or is initially fluid absent. If the interface region does not allow fluid transport then as temperature increases, a fluid pressure greater than lithostatic can develop. At some degree of over-pressure relative to rock pressure, the fluid hydraulically fractures the rock and a gradient in fluid composition away from the contact can be produced. These phenomena occur at the compositional interfaces whenever univariant reactions in the differing layers cross on a temperature vs. mole fraction of CO₂ diagram with slopes of opposite sign. The first occurrence of these reaction products at lithological contacts delineates an isograd that defines temperature as well as the mole fraction of CO₂ at constant pressure in systems open to fluid transport. These isograds can be contrasted with fluid-producing isograds in closed systems. As an illustration of possible effects, the reactions quartz + clinozoisite + muscovite = anorthite + K-feldspar + H₂O and phlogopite + quartz + calcite = tremolite + K-feldspar + H₂O + CO₂ at 4 kbar are analysed and equations for fluid production and transport are developed.     

SHRIMP U–Pb ages of diagenetic and hydrothermal xenotime from the Archaean Witwatersrand Supergroup of South Africa

SHRIMP dating of xenotime overgrowths on detrital zircon grains can constrain maximum durations since diagenesis and therefore provide minimum dates of sediment deposition. Thus, xenotime dating has significant economic application to Precambrian sediment-hosted ore deposits, such as Witwatersrand Au–U, for which there are no precise depositional ages. The growth history of xenotime in the Witwatersrand Supergroup is texturally complex, with several phases evident. The oldest authigenic xenotime ²⁰⁷Pb/²⁰⁶Pb age obtained in sandstone underlying the Vaal Reef is 2764 ± 5 Myr (1 σ), and most likely represents a mixture of diagenetic and hydrothermal growth. Nevertheless, this represents the oldest authigenic mineral age yet recorded in the sequence and provides a minimum age of deposition. Other xenotime data record a spread of ages that correspond to numerous post-diagenetic thermotectonic events (including a Ventersdorp event at ≈ 2720 Ma) up to the ≈2020 Ma Vredefort event.