On a global expansion of the disturbing function in the planar elliptic restricted three-body problem

Starting with a simple Taylor-based expansion of the inverse of the distance between two bodies, we are able to obtain a series expansion of the disturbing function of the three-body problem (planar elliptic case) which is valid for all points of the phase space outside the immediate vicinity of the collision points. In particular, the expansion is valid for very high values of the eccentricity of the perturbed body. Furthermore, in the case of an interior mean-motion resonant configuration, the above-mentioned expression is easily averaged with respect to the synodic period, yielding once again a global expansion of (R) valid for very high eccentricities.Comparisons between these results and the numerically computed exact function are presented for various resonances and values of the eccentricity. Maximum errors are determined in each case and their origin is established. Lastly, we discuss the applicability of the present expansion to practical problems.     

Dust supply from river floodplains: the case of the lower Huang He (Yellow River) recorded in a loess–palaeosol sequence from the Mangshan Plateau

The Mangshan Plateau is located on the south bank of the Huang He (Yellow River) just west of the city of Zhengzhou, well outside the Loess Plateau in central China. Mixing models of the grain‐size data indicate that the loess deposits are mixtures of three loess components. Comparison of the mixing model with existing models established for a series of loess–palaeosol sequences from the Loess Plateau indicates that the Mangshan loess has been supplied from a proximal dust source, the Huang He floodplain, during major dust outbreaks. The high accumulation rates, the composition of the loess components, and especially the high proportions of a sandy loess component support this. Owing to the exceptionally high accumulation rates, the Mangshan grain size, magnetic susceptibility and carbonate records provide a high‐resolution archive of environmental and climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Tephrochronology and tephrostratigraphy of two Pleistocene continental fossiliferous successions from central Italy

The stratigraphical context of two Middle Pleistocene fossiliferous palaeosols from Central Italy (Abruzzo and Tuscany) have been studied. Small mammals and molluscs occur in both palaeosols, which are covered by tephra layers that were analysed using an interdisciplinary approach. Application of fission‐track dating to apatites separated from the Case Picconetto tephra (Pescara, Abruzzo), yielded an age of 0.48 ± 0.04 Ma, indistinguishable from those previously determined for the Campani Quarry (Lower Valdarno, Tuscany) (0.46 ± 0.05 Ma and 0.48 ± 0.05 Ma). Geochemical and petrographic investigations indicate that these tephra originated from different volcanoes, the Alban Hills Volcanic Complex and the Vico Volcano (Latium) respectively. Small mammal and mollusc assemblages indicate different palaeoclimatic and palaeoenvironmental conditions for the Case Picconetto and Campani Quarry palaeosols. Warm and humid conditions can be inferred for the Campani Quarry site, whereas open and cold conditions can be inferred for Case Picconetto. On the basis of faunal data, fission‐track dates and attribution of tephra to specific volcanic eruptions, we suggest a correlation of these faunas with marine oxygen isotope stage 14 (Case Picconetto) and with marine oxygen isotope stage 11 (Campani Quarry), respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.