ALBERT

Description: The Nanga Parbat Massif (NPM), Pakistan Himalaya, is an exhumed tract of Indian continental crust and represents an area of active crustal thickening and exhumation. While the most effective way to study the NPM at depth is through seismic imaging, interpretation depends upon knowledge of the seismic properties of the rocks. Gneissic, ‘mylonitic’ and cataclastic rocks emplaced at the surface were sampled as proxies for lithologies and fabrics currently accommodating deformation at depth. Mineral crystallographic preferred orientations (CPO) were measured via scanning electron microscope (SEM)/electron backscatter diffraction (EBSD), from which three-dimensional (3D) elastic constants, seismic velocities and anisotropies were predicted. Micas make the main contribution to sample anisotropy. Background gneisses have highest anisotropy (up to 10.4% shear-wave splitting, AVs) compared with samples exhibiting localized deformations (e.g. ‘mylonite’, 4.7% AVs; cataclasite, 1% AVs). Thus, mylonitic shear zones may be characterized by regions of low anisotropy compared to their wall rocks. CPO-derived sample elastic constants were used to construct seismic models of NPM tectonics, through which P-, S- and converted waves were ray-traced. Foliation orientation has dramatic effects on these waves. The seismic models suggest dominantly pure-shear tectonics for the NPM involving horizontal compression and vertical stretching, modified by localized ductile and brittle (‘simple’) shear deformations.

Description: Improving the accuracy of subsurface imaging is commonly the main incentive for including the effects of anisotropy in seismic processing. However, the anisotropy itself holds valuable information about rock properties and, as such, can be viewed as a seismic attribute. Here we summarize results from an integrated project that explored the potential to use observations of seismic anisotropy to interpret lithological and fluid properties (the SAIL project). Our approach links detailed petrofabric analyses of reservoir rocks, laboratory based measurements of ultrasonic velocities in core samples, and reservoir-scale seismic observations. We present results for the Clair field, a CarboniferousDevonian reservoir offshore Scotland, west of the Shetland Islands. The reservoir rocks are sandstones that are variable in composition and exhibit anisotropy on three length-scales: the crystal, grain and fracture scale. We have developed a methodology for assessing crystal-preferred-orientation (CPO) using a combination of electron back-scattered diffraction (EBSD), X-ray texture goniometry (XRTG) and image analysis. Modal proportions of individual minerals are measured using quantitative X-ray diffraction (QXRD). These measurements are used to calculate the intrinsic anisotropy due to CPO via Voigt-Reuss-Hill averaging of individual crystal elasticities and their orientations. The intrinsic anisotropy of the rock is controlled by the phyllosilicate content and to a lesser degree the orientation of quartz and feldspar; the latter can serve as a palaeoflow indicator. Our results show remarkable consistency in CPO throughout the reservoir and allow us to construct a mathematical model of reservoir anisotropy. A comparison of CPO-predicted velocities and those derived from laboratory measurements of ultrasonic signals allows the estimation of additional elastic compliance terms due to grain-boundary interactions. The results show that the CPO estimates are good proxies for the intrinsic anisotropy of the clean sandstones. The more micaceous rocks exhibit enhanced anisotropy due to interactions between the phyllosilicate grains. We then compare the lab-scale predictions with reservoir-scale measurements of seismic anisotropy, based on amplitude variation with offset and azimuth (AVOA) analysis and non-hyperbolic moveout. Our mathematical model provides a foundation for interpreting the reservoir-scale seismic data and improving the geological modelling of complex reservoirs. The observed increases in AVOA signal with depth can only be explained with an increase in fracturing beneath the major unit boundaries, rather than a change in intrinsic CPO properties. In general, the style and magnitude of anisotropy in the Clair field appears to be indicative of reservoir quality.

Description: Progressive deformation of upper mantle rocks via dislocation creep causes their constituent crystals to take on a non-random orientation distribution (crystallographic preferred orientation or CPO) whose observable signatures include shear-wave splitting and azimuthal dependence of surface wave speeds. Comparison of these signatures with mantle flow models thus allows mantle dynamics to be unraveled on global and regional scales. However, existing self-consistent models of CPO evolution are computationally expensive when used with 3-D and/or time-dependent convection models. Here we propose a new method, called ANPAR, which is based on an analytical parametrization of the crystallographic spin predicted by the second-order (SO) self-consistent theory. Our parametrization runs 2–6  x  10 4 times faster than the SO model and fits its predictions for CPO and crystallographic spin with a variance reduction 〉99 per cent. We illustrate the ANPAR model predictions for the deformation of olivine with three dominant slip systems, (010)[100], (001)[100] and (010)[001], for three uniform deformations (uniaxial compression, pure shear and simple shear) and for a corner-flow model of a spreading mid-ocean ridge.

Description: Indonesia is arguably one of the tectonically most complex regions on Earth today due to its location at the junction of several major tectonic plates and its long history of collision and accretion. It is thus an ideal location to study the interaction between subducting plates and mantle convection. Seismic anisotropy can serve as a diagnostic tool for identifying various subsurface deformational processes, such as mantle flow, for example. Here, we present novel shear wave splitting results across the Indonesian region. Using three different shear phases (local S, SKS, and downgoing S) to improve spatial resolution of anisotropic fabrics allows us to distinguish several deformational features. For example, the block rotation history of Borneo is reflected in coast-parallel fast directions, which we attribute to fossil anisotropy. Furthermore, we are able to unravel the mantle flow pattern in the Sulawesi and Banda region: We detect toroidal flow around the Celebes Sea slab, oblique corner flow in the Banda wedge, and sub-slab mantle flow around the arcuate Banda slab. We present evidence for deep, sub-520 km anisotropy at the Java subduction zone. In the Sumatran backarc, we measure trench-perpendicular fast orientations, which we assume to be due to mantle flow beneath the overriding Eurasian plate. These observations will allow to test ideas of, for example, slab–mantle coupling in subduction regions.

Description: Seismic anisotropy in the Earth's lowermost mantle (D '' ) is often attributed to the alignment of MgSiO 3 post-perovskite (ppv) by the movement of dislocations in response to mantle flow. However, ppv's plastic yield surface is not known; nor do we know if this is the main deformation mechanism. We make use of a heterogeneous, generally anisotropic model of elasticity in D '' derived from a 3-D model of mantle flow, which is obtained by inversion of geophysical observables. Unlike previous approaches, completely general, 3-D flow and full anisotropy are permitted, yielding more information to compare with observations than has been possible before. We model observations of anisotropy in D '' by calculating the shear wave splitting predicted in ScS waves for a series of models of ppv plasticity. We find that observations in regions of the lowermost mantle beneath subduction zones are best fit by a model which accommodates slip on (010). Our results show that, within one standard deviation, slip on (010)—or a mechanism giving the same style of anisotropy—explains D '' anisotropy beneath these regions.

Description: Hudson Bay Lithospheric Experiment (HuBLE) was designed to understand the processes that formed Laurentia and the Hudson Bay basin within it. Receiver function analysis shows that Archaean terranes display structurally simple, uniform thickness, felsic crust. Beneath the Palaeoproterozoic Trans-Hudson Orogen (THO), thicker, more complex crust is interpreted as evidence for a secular evolution in crustal formation from non-plate-tectonic in the Palaeoarchaean to fully developed plate tectonics by the Palaeoproterozoic. Corroborating this hypothesis, anisotropy studies reveal 1.8 Ga plate-scale THO-age fabrics. Seismic tomography shows that the Proterozoic mantle has lower wavespeeds than surrounding Archaean blocks; the Laurentian keel thus formed partly in post-Archaean times. A mantle transition zone study indicates ‘normal’ temperatures beneath the Laurentian keel, so any cold mantle down-welling associated with the regional free-air gravity anomaly is probably confined to the upper mantle. Focal mechanisms from earthquakes indicate that present-day crustal stresses are influenced by glacial rebound and pre-existing faults. Ambient-noise tomography reveals a low-velocity anomaly, coincident with a previously inferred zone of crustal stretching, eliminating eclogitization of lower crustal rocks as a basin formation mechanism. Hudson Bay is an ephemeral feature, caused principally by incomplete glacial rebound. Plate stretching is the primary mechanism responsible for the formation of the basin itself.

Description: Cap rock integrity is an essential characteristic of any reservoir to be used for long-term $${\mathrm{CO}}_{2}$$ storage. Seismic AVOA (amplitude variation with offset and azimuth) techniques have been applied to map HTI anisotropy near the cap rock of the Weyburn field in southeast Saskatchewan, Canada, with the purpose of identifying potential fracture zones that may compromise seal integrity. This analysis, supported by modeling, observes the top of the regional seal (Watrous Formation) to have low levels of HTI anisotropy, whereas the reservoir cap rock (composite Midale Evaporite and Ratcliffe Beds) contains isolated areas of high intensity anisotropy, which may be fracture-related. Properties of the fracture fill and hydraulic conductivity within the inferred fracture zones are not constrained using this technique. The predominant orientations of the observed anisotropy are parallel and normal to the direction of maximum horizontal stress (northeast–southwest) and agree closely with previous fracture studies on core samples from the reservoir. Anisotropy anomalies are observed to correlate spatially with salt dissolution structures in the cap rock and overlying horizons as interpreted from 3D seismic cross sections.

Description: Observations of shear wave splitting provide unambiguous evidence of the presence of anisotropy in the Earth's lowermost mantle, a region known as D '' . Much recent work has attempted to use these observations to place constraints on strain above the core–mantle boundary (CMB), as this may help map flow throughout the mantle. Previously, this interpretation has relied on the assumption that waves can be modelled as infinite-frequency rays, or that the Earth is radially symmetric. Due to computational constraints it has not been possible to test these approximations until now. We use fully 3-D, generally anisotropic simulations of ScS waves at the frequencies of the observations to show that ray methods are sometimes inadequate to interpret the signals seen. We test simple, uniform models, and for a D '' layer as thin as 50 km, significant splitting may be produced, and we find that recovered fast orientations usually reflect the imposed fast orientation above the CMB. Ray theory in these cases provides useful results, though there are occasional, notable differences between forward methods. Isotropic models do not generate apparent splitting. We also test more complex models, including ones based on our current understanding of mineral plasticity and elasticity in D '' . The results show that variations of anisotropy over even several hundred kilometres cause the ray-theoretical and finite-frequency calculations to differ greatly. Importantly, models with extreme mineral alignment in D '' yield splitting times not dissimilar to observations ( t ≤ 3 s), suggesting that anisotropy in the lowermost mantle is probably much stronger than previously thought—potentially ~10 per cent shear wave anisotropy or more. We show that if the base of the mantle is as complicated as we believe, future studies of lowermost mantle anisotropy will have to incorporate finite-frequency effects to fully interpret observations of shear wave splitting.