Range of validity of seismic ray and beam methods in general inhomogeneous media – I. General theory

Summary. The limitations of asymptotic wave theory and its geometrical manifestations are newly formalized and scrutinized. Necessary and sufficient conditions for the existence of acoustic and seismic rays and beams in general inhomogeneous media are expressed in terms of new physical parameters: the threshold frequency ω₀ associated with the P/S decoupling condition, the cut-off frequency ω_c associated with the radiation-zone condition, the total curvature of the wavefront and the Fresnel-zone radius.
The analysis is facilitated with the introduction of a new ancillary functional – the hypereikonal which is capable of representing ordinary as well as evanescent waves. The hypereikonal is the natural extension of the eikonal theory.
With the aid of the above new parameters, simple conditions are obtained for the decoupled far field, the decoupled near field, two point dynamic ray tracing, paraxial wavefields and Gaussian beams.     

An approximate method for evaluating the shear band thickness in saturated soils

A formula for the thickness of a shear band formed in saturated soils under a simple shear or a combined stress state has been proposed. It is shown that the shear band thickness is dependent on the pore pressure properties of the material and the dilatancy rate, but is independent of the details of the combined stress state. This is in accordance with some separate experimental observations. Copyright © 2004 John Wiley & Sons, Ltd.     

A WKBJ spectral method for computation of SV synthetic seismograms in a cylindrically symmetric medium

Summary. A method of synthetic seismogram computation for teleseismic SV -waves is developed in order to treat quantitatively SV -waves in problems of body wave source inversion and source—receiver structure studies. The method employs WKBJ theory for a generalized ray in a vertically inhomogeneous half-space and the propagator matrix technique for waves in near-surface homogeneous layers. Wavenumber integration is done along the real axis of the wavenumber plane and anelasticity is included by using complex velocity in all regions of the earth model. The near-surface source structure is taken into account in the computation for the case of the shallow source by allowing a point source to be located in the homogeneous layers. Source and receiver area structures are also allowed to differ. A general moment tensor point source is considered.     

The geochemistry of lead in rivers,estuaries and the continental shelf of the Southeastern United States

Transport by southeastern rivers has insignificant influence on the lead concentration of southeastern shelf waters. If the rate of mobilization and transport of lead by these rivers is representative of uncontaminated fluvial transport during the Pleistocene, only about 5% of the prehistoric output of lead to the North Atlantic can be accounted for by river input.Lead concentrations in southeastern shelf waters are similar to those of North Atlantic Water from the upper 1000 m, which is probably similar to the intrusion source water for the shelf. Atmospheric inputs to the shelf of the same magnitude as observed for the Western North Atlantic are difficult to reconcile given the residence time of shelf waters and their lead concentration unless the rate of loss of lead to shelf sediments is about the same as the atmospheric flux.     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Zircon U–Pb age, trace element and Hf isotope composition of Kongling terrane in the Yangtze Craton: refining the timing of Palaeoproterozoic high-grade metamorphism

U–Pb age, trace element and Hf isotope compositions of zircon were analysed for a metasedimentary rock and two amphibolites from the Kongling terrane in the northern part of the Yangtze Craton. The zircon shows distinct morphological and chemical characteristics. Most zircon in an amphibolite shows oscillatory zoning, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high formation temperature, high trace element contents, clear negative Eu anomaly, as well as HREE-enriched patterns, suggesting that it is igneous. The zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2857 ± 8 Ma, representing the age of the magmatic protolith. The zircon in the other two samples is metamorphic. It has low Th/U ratios, low trace element concentrations, variable HREE contents (33.8 ≥ Lu_N≥2213; 14.7 ≤ Lu_N/Sm_N ≤ 354) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.000030–0.001168). The data indicate that the zircon formed in the presence of garnet and under upper amphibolite facies conditions. The metamorphic zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2010 ± 13 Ma. These results combined with previously obtained Palaeoproterozoic metamorphic ages suggest a c. 2.0 Ga Palaeoproterozoic collisional event in the Yangtze Craton, which may result from the assembly of the supercontinent Columbia. The zircon in two samples yields weighted mean two-stage Hf model ( T _DM2) ages of 3217 ± 110 and 2943 ± 50 Ma, respectively, indicating that their protoliths were mainly derived from Archean crust.     

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.