Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Canonical theory of perturbations using eccentric anomaly as independent variable

It is shown in this paper how to build a canonical transformation of variables, so that the eccentric anomaly becomes the new independent variable. In the case of eccentric elliptical orbits it changes the equations of motion so, that they can be integrated analytically to any order of approximation comparatively easy.     

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.     

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.     

An explicit representation of steady state response of a beam on an elastic foundation to moving harmonic line loads

A closed‐form deflection response of a beam rest is presented in this paper using the integral transform method. The theory of linear partial differential equations is used to represent the deflection of beam subjected to a moving harmonic line load in integration form. The solution is finally carried out using the inverse Fourier transform. To evaluate the integration analytically, poles of the integrand are identified with the help of algebraic equation theory. Residue theorem is then utilized to represent the integration as a contour integral in the complex plane. Closed‐form deflections and numerical results are provided for different combinations of load frequency and velocity. Copyright © 2002 John Wiley & Sons, Ltd.     

Solution to the rotation of the elastic earth by method of rigid dynamics

Hamiltonian mechanics is applied to the problem of the rotation of the elastic Earth. We first show the process for the formulation of the Hamiltonian for rotation of a deformable body and the derivation of the equations of motion from it. Then, based on a simple model of deformation, the solution is given for the period of Euler motion, UT1 and the nutation of the elastic Earth. In particular it is shown that the elasticity of the Earth acts on the nutation so as to decrease the Oppolzer terms of the nutation of the rigid Earth by about 30 per cent. The solution is in good agreement with results which have been obtained by other, different approaches.     

Porphyroblast inclusion-trail orientation data: eppure non son girate*!

Extensive examination of large numbers of spatially orientated thin sections of orientated samples from orogens of all ages around the world has demonstrated that porphyroblasts do not rotate relative to geographical coordinates during highly non-coaxial ductile deformation of the matrix subsequent to their growth. This has been demonstrated for all tectonic environments so far investigated. The work also has provided new insights and data on metamorphic, structural and tectonic processes including: (1) the intimate control of deformation partitioning on metamorphic reactions; (2) solutions to the lack of correlation between lineations that indicate the direction of movement within thrusts and shear zones, and relative plate motion; and (3) a possible technique for determining the direction of relative plate motion that caused orogenesis in ancient orogens.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.