Anisotropic attenuation of stratified viscoelastic media

An important cause of seismic anisotropic attenuation is the interbedding of thin viscoelastic layers. However, much less attention has been devoted to layer‐induced anisotropic attenuation. Here, we derive a group of unified weighted average forms for effective attenuation from a binary isotropic, transversely isotropic‐ and orthorhombic‐layered medium in the zero‐frequency limit by using the Backus averaging/upscaling method and analyse the influence of interval parameters on effective attenuation. Besides the corresponding interval attenuation and the real part of stiffness, the contrast in the real part of the complex stiffness is also a key factor influencing effective attenuation. A simple linear approximation can be obtained to calculate effective attenuation if the contrast in the real part of stiffness is very small. In a viscoelastic medium, attenuation anisotropy and velocity anisotropy may have different orientations of symmetry planes, and the symmetry class of the former is not lower than that of the latter. We define a group of more general attenuation‐anisotropy parameters to characterize not only the anisotropic attenuation with different symmetry classes from the anisotropic velocity but also the elastic case. Numerical tests reveal the influence of interval attenuation anisotropy, interval velocity anisotropy and the contrast in the real part of stiffness on effective attenuation anisotropy. Types of effective attenuation anisotropy for interval orthorhombic attenuation and interval transversely isotropic attenuation with a vertical symmetry (vertical transversely isotropic attenuation) are controlled only by the interval attenuation anisotropy. A type of effective attenuation anisotropy for interval TI attenuation with a horizontal symmetry (horizontal transversely isotropic attenuation) is controlled by the interval attenuation anisotropy and the contrast in the real part of stiffness. The type of effective attenuation anisotropy for interval isotropic attenuation is controlled by all three factors. The magnitude of effective attenuation anisotropy is positively correlated with the contrast in the real part of the stiffness. Effective attenuation even in isotropic layers with identical isotropic attenuation is anisotropic if the contrast in the real part of stiffness is non‐zero. In addition, if the contrast in the real part of stiffness is very small, a simple linear approximation also can be performed to calculate effective attenuation‐anisotropy parameters for interval anisotropic attenuation.     

Geometrical characteristics of phase and group velocity surfaces in anisotropic media

The phase and group velocity surfaces are essential for wave propagation in anisotropic media. These surfaces have certain features that, especially, for shear waves result in complications for modelling and inversion of recorded wavefields. To analyse wave propagation in an anisotropic model, it is important to identify these features in both the phase and group domains. We propose few characteristics for this analysis: the energy flux angle, decomposed in the polar and azimuth angle correction angles and enhancement factor, which is able to characterize both singularity points and triplication zones. The very simple equation that controls the triplications is derived in the phase domain. The proposed characteristics are illustrated for elastic and acoustic anisotropic models of different symmetry classes.     

Discovery of an albite gneiss from the Ile de Groix (Armorican Massif, France): geochemistry and LA-ICP-MS U–Pb geochronology of its Ordovician protolith

For the first time, an albite orthogneiss has been recognised and dated within the HP–LT blueschist facies metabasites and metapelites of the Ile de Groix. It is characterised by a peraluminous composition, high LILE, Th and U contents, MORB-like HREE abundances and moderate Nb and Y values. A U–Pb age of 480.8?±?4.8?Ma was obtained by LA-ICP-MS dating of zircon and titanite. It is interpreted as the age of the magmatic emplacement during the Early Ordovician. Morphologically different zircon grains yield late Neoproterozoic ages of 546.6–647.4?Ma. Zircon and titanite U–Pb ages indicate that the felsic magmatism from the Ile de Groix is contemporaneous with the acid, pre-orogenic magmatism widely recognised in the internal zones of the Variscan belt, related to the Cambro-Ordovician continental rifting. The magmatic protolith probably inherited a specific chemical composition from a combination of orogenic, back-arc and anorogenic signatures because of partial melting of the Cadomian basement during magma emplacement. Besides, the late Devonian U–Pb age of 366?±?33?Ma obtained for titanite from a blueschist facies metapelite corresponds to the age of the HP–LT peak metamorphism.     

Traveltime parameters in a tilted elliptical anisotropic medium‡

In this paper, we derive analytical expressions for one‐way and two‐way kinematical parameters in elliptical tilted transverse isotropy media. We show that the homogeneous elliptical tilted transverse isotropy models result in hyperbolic moveout with a reflection point sideslip x₀, which can be considered as an additional traveltime parameter for one‐way wave propagation. For homogeneous elliptical tilted transverse isotropy models we show that the inversion of one‐way traveltime parameters suffers from the ambiguity for large tilts. It is shown that the accuracy of the inversion is sensitive to the error in x₀. We also derive and invert the traveltime parameters for a vertically heterogeneous elliptical tilted transverse isotropy model with a tilt gradually changing with depth. The a priori knowledge of parameter δ is very important for inversion. The wrong choise of this parameter results in significant errors in inverted model parameters.     

Noble gases in presolar diamonds III: Implications of ion implantation experiments with synthetic nanodiamonds

Abstract— A series of experiments carried out by Koscheev et al. (1998, 2001, 2004, 2005) showed that the bimodal release of heavy noble gases from meteoritic nanodiamonds can be reproduced by a single implanted component. This paper investigates the implications of this result for interpreting the noble gas compositions of meteoritic nanodiamonds and for their origin and history. If the bimodal release exhibited by meteorite diamonds reflects release of the P3 noble gas component, then the composition inferred for the pure Xe‐HL end member changes slightly, the excesses of heavy krypton isotopes that define Kr‐H become less extreme, evidence appears for a Kr‐L component, and the nucleosynthetic contribution to argon becomes much smaller. After correction for cosmogenic neon inherited from the host meteorites, the neon in presolar diamonds shows evidence for pre‐irradiation, perhaps in interstellar space, and a nucleosynthetic component perhaps consistent with a supernova source. After a similar correction, helium also shows evidence for presolar irradiation and/or a nucleosynthetic component. For the case of presolar irradiation, due to the small size of the diamonds, a large entity must have been irradiated and recoiling product nuclei collected by the nanodiamonds. The high ³He/²¹Ne ratio (?43) calls for a target with a (C + O)/heavier‐element ratio higher than in chondritic abundances. Bulk gas + dust (cosmic abundances) meet this criteria, as would solids enriched in carbonaceous material. The long recoil range of cosmogenic ³He argues against a specific phase. The excess ³He in presolar diamonds may represent trapped cosmic rays rather than cosmogenic ³He produced in the vicinity of the diamond crystals.     

Effect of Decoupling of Lithospheric Plates on the Observed Geoid

A joint effect of weak zones, dividing lithospheric plates, and lateral viscosity variations (LVV) in the whole mantle on the observed geoid is investigated by a new numerical approach. This technique is based on the substantially revised method introduced by Zhang and Christensen (Geophys J Int 114:531–547, 1993) for solving the Navier–Stokes–Poisson equations in the spectral domain with strong LVV. Weak plate boundaries (WPB) are introduced based on an integrated global model of plate boundary deformations GSRM (Kreemer et al. in Geophys J Int 154:8–34, 2003). The effect of WPB on the geoid is significant and reaches ?40 to 70 m with RMS ~20 m. The peaks are observed over large subduction zones in South America and the southwestern Pacific in agreement with previous studies. The positive geoid anomaly in South America could be explained largely by a dynamic effect of decoupling of the Nazca and South American plates. The negative changes of the geoid mostly relate to mid-oceanic ridges. The amplitude of the effect depends on the viscosity contrasts at WPB compared with the plate viscosity until its value reaches the limit of 2.5–3 orders of magnitude. This value might be considered as a level at which the plates are effectively decoupled. The effect of WPB exceeds the effect of LVV in the whole mantle and generally does not correlate with it. However, inclusion of LVV reduces the geoid perturbations due to WPB by about 10 m. Therefore, it is important to consider all factors together. The geoid changes mainly result from changes of the dynamic topography, which are about ?300 to +500 m. The obtained results show that including WPB may significantly improve the reliability of instantaneous global dynamic models.     

Sensitivity of the southern annular mode to greenhouse gas emission scenarios

The leading mode of southern hemisphere (SH) climatic variability, the southern annular mode (SAM), has recently seen a shift towards its positive phase due to stratospheric ozone depletion and increasing greenhouse gas (GHG) concentrations. Here we examine how sensitive the SAM (defined as the leading empirical orthogonal function of SH sea level pressure anomalies) is to future GHG concentrations. We determine its likely evolution for three intergovernmental panel on climate change (IPCC) special report on emission scenarios (SRES) for austral summer and winter, using a multi-model ensemble of IPCC fourth assessment report models which resolve stratospheric ozone recovery. During the period of summer ozone recovery (2000–2050), the SAM index exhibits weakly negative, statistically insignificant trends due to stratospheric ozone recovery which offsets the positive forcing imposed by increasing GHG concentrations. Thereafter, positive SAM index trends occur with magnitudes that show sensitivity to the SRES scenario utilised, and thus future GHG emissions. Trends are determined to be strongest for SRES A2, followed by A1B and B1, respectively. The winter SAM maintains a similar dependency upon GHG as summer, but over the entire twenty-first century and to a greater extent. We also examine the influence of ozone recovery by comparing results to models that exclude stratospheric ozone recovery. Projections are shown to be statistically different from the aforementioned results, highlighting the importance of ozone recovery in governing SAM-evolution. We therefore demonstrate that the future SAM will depend both upon GHG emissions and stratospheric ozone recovery.     

Low-frequency anisotropy in fractured and layered media

Computing effective medium properties is very important when upscaling data measured at small scale. In the presence of stratigraphic layering, seismic velocities and anisotropy parameters are scale and frequency dependent. For a porous layer permeated by aligned fractures, wave-induced fluid flow between pores and fractures can also cause significant dispersion in velocities and anisotropy parameters. In this study, we compare the dispersion of anisotropy parameters due to fracturing and layering at low frequencies. We consider a two-layer model consisting of an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce anisotropy parameters dispersion due to fractures (meso-scale) in the sand layer. This intrinsic dispersion is added to anisotropy parameters dispersion induced by layering (macro-scale) at low frequencies. We derive the series coefficients that control the behaviour of anisotropy parameters at low frequencies. We investigate the influences of fracture length and fracture density on fracturing effect, layering effect and combined effect versus frequency and volume fraction of sand layer. Numerical modelling results indicate that the frequency dependence due to layering is not always the dominant effect of the effective properties of the medium. The intrinsic dispersion is not negligible compared with the layering effect while evaluating the frequency-dependent properties of the layered medium.     

Frequency-dependent PP and PS reflection coefficients in fractured media

When a porous layer is permeated by mesoscale fractures, wave-induced fluid flow between pores and fractures can cause significant attenuation and dispersion of velocities and anisotropy parameters in the seismic frequency band. This intrinsic dispersion due to fracturing can create frequency-dependent reflection coefficients in the layered medium. In this study, we derive the frequency-dependent PP and PS reflection coefficients versus incidence angle in the fractured medium. We consider a two-layer vertical transverse isotropy model constituted by an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce the intrinsic dispersion due to fracturing in the sand layer. Based on the series coefficients that control the behaviour of velocity and anisotropy parameters in the fractured medium at low frequencies, we extend the conventional amplitude-versus-offset equations into frequency domain and derive frequency-dependent amplitude-versus-offset equations at the elastic–anelastic surface. Increase in fracture length or fracture density can enlarge the frequency dependence of amplitude-versus-offset attributes of PP and PS waves. Also, the frequency dependence of magnitude and phase angle of PP and PS reflection coefficients increases as fracture length or fracture density increases. Amplitude-versus-offset type of PP and PS reflection varies with fracture parameters and frequency. What is more, fracture length shows little impact on the frequency-dependent critical phase angle, while the frequency dependence of the critical phase angle increases with fracture density.