Abstract Low seismic velocities and/or strong seismic anisotropy are often interpreted as being caused by partial melt. To better understand this, we used numerical modelling, varying the shape and amount of melt, to show how seismic phases are affected by melt. We observed that seismic waves are more sensitive to the shape than to the amount of melt. Rayleigh wave velocities were almost always reduced in the presence of melt, while Pn/wide-angle P-wave refraction and Love wave velocities showed low velocity anomalies for vertically aligned melt, but little anomaly for horizontally aligned melt. These data can therefore be used to determine the alignment of melt. Shear wave splitting/receiver functions showed strong anisotropy and can be used to constrain the strike of vertically aligned partial melt. We showed that melt in the mantle beneath Ethiopia is probably stored in low aspect ratio disc-like inclusions, suggesting that the melt is not in textural equilibrium. We estimated that 2–7% of the vertically aligned melt is stored beneath the Main Ethiopian Rift, >6% of the horizontally and vertically aligned melt is stored beneath the Red Sea Rift and 1–6% of the horizontally aligned melt is stored beneath the Danakil microplate. This supports the idea of strong shear-derived segregation of melt in the narrow Main Ethiopian Rift compared with that observed beneath Afar.

Abstract 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.

Abstract 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.

Abstract 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 Carboniferous–Devonian 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.

Abstract The rifting of continents and eventual formation of ocean basins is a fundamental component of plate tectonics, yet the mechanism for break-up is poorly understood. The East African Rift System (EARS) is an ideal place to study this process as it captures the initiation of a rift in the south through to incipient oceanic spreading in north-eastern Ethiopia. Measurements of seismic anisotropy can be used to test models of rifting. Here we summarize observations of anisotropy beneath the EARS from local and teleseismic body-waves and azimuthal variations in surface-wave velocities. Special attention is given to the Ethiopian part of the rift where the recent EAGLE project has provided a detailed image of anisotropy in the portion of the Ethiopian Rift that spans the transition from continental rifting to incipient oceanic spreading. Analyses of regional surface-waves show sub-lithospheric fast shear-waves coherently oriented in a northeastward direction from southern Kenya to the Red Sea. This parallels the trend of the deeper African superplume, which originates at the core–mantle boundary beneath southern Africa and rises towards the base of the lithosphere beneath Afar. The pattern of shear-wave anisotropy is more variable above depths of 150 km. Analyses of splitting in teleseismic phases (SKS) and local shear-waves within the rift valley consistently show rift-parallel orientations. The magnitude of the splitting correlates with the degree of magmatism and the polarizations of the shear-waves align with magmatic segmentation along the rift valley. Analysis of surf ace-wave propagation across the rift valley confirms that anisotropy in the uppermost 75 km is primarily due to melt alignment. Away from the rift valley, the anisotropy agrees reasonably well with the pre-existing Pan-African lithospheric fabric. An exception is the region beneath the Ethiopian plateau, where the anisotropy is variable and may correspond to pre-existing fabric and ongoing melt-migration processes. These observations support models of magma-assisted rifting, rather than those of simple mechanical stretching. Upwellings, which most probably originate from the larger super-plume, thermally erode the lithosphere along sites of pre-existing weaknesses or topographic highs. Decompression leads to magmatism and dyke injection that weakens the lithosphere enough for rifting and the strain appears to be localized to plate boundaries, rather than wider zones of deformation.

Abstract The link between petrofabric (LPO) and seismic properties of an amphibolite-facies quartzo-feldspathic shear zone is explored using SEM/EBSD. The shear-zone LPO develops by a combination of slip systems and dauphine twinning, with a -maximum parallel to lineation ( X ) and c -maximum normal to mylonitic foliation ( XY ). The LPO are used to predict elastic parameters, from which the three-dimensional seismic properties of different shear-zone regions are derived. Results suggest that LPO evolution is reflected in the seismic properties but the precise impact is not simple. In general, the P-wave velocity ( V P ) minimum is parallel to the a -axis maximum; the direction of maximum shear-wave splitting ( AV S ) is parallel to mylonitic foliation; and the V P maximum and AV S minimum are parallel to the c -axis maximum. The seismic anisotropy predicted is significant and increases from shear zone wallrock to mature mylonite. The P-wave anisotropy ranges from 11 to 13%, fast and slow shear waves’ anisotropies range from 6 to 15% and the magnitude of shear-wave splitting ranges from 9 to 16%. Nevertheless, such anisotropy requires a considerable thickness of rock with this LPO before it becomes seismically visible (i.e. 100s of m for local earthquakes; 10s of m for controlled source experiments). However, reflections and mode conversions provide much better resolution, particularly across tectonic boundaries. The low symmetry and strong anisotropy due to the LPO suggest that multi-azimuth wide-angle reflection data may be useful in the determination of the deformation characteristics of deep shear zones.

Abstract The Superior Province, which forms the nucleus of North America, is the largest preserved Archaean crustal block in the world and may have originated as the result of a widespread crustal accretion event ( c. 2.7 Ga) manifested in Archaean cratons worldwide. An understanding of the accretionary evolution of this craton is the objective of the continuing Lithoprobe transect in the Western Superior Province, Canada. The geophysical components of the transect include seismic reflection and refraction, magnetotellurics and teleseismic experiments. The Teleseismic Western Superior Transect (TWiST) was designed to explore the structural and physical properties of the subcrustal lithosphere and their implications for proposed accretionary models. A north-south-trending array of 17 broadband three-component seismometers was deployed between May and November 1997. Surface-wave analyses, SKS-splitting studies and travel-time tomography show variations in the velocity structure and anisotropy between the southern end of the transect, a region affected by Keweenawan rifting, and the northern part, which lies in the Proterozoic Trans-Hudson shear zone. Surface waves reveal evidence for a thin high-velocity layer, 5–20 km thick, beneath a 37–43 km thick crust and above c. 250 km of high-velocity continental root. This thin layer is also visible in wide-angle refraction data from the southern end of the line and may be evidence of underplating during terrane accretion. Discrepancies in the Love and Rayleigh waves and surface-wave particle motions show evidence for an anisotropic mantle. SKS analysis shows large amounts (up to 2 s) of shear-wave splitting with a roughly eastwest trend in the fast-shear-wave polarization direction for most stations. This conforms with crustal deformation trends. Stations in the younger Trans-Hudson orogen show much less splitting. Detailed analysis at a permanent station in the Western Superior shows evidence for two layers of anisotropy. A thinner upper layer is aligned with the surface geology, indicating crust-mantle coupling during craton formation, whereas a thicker lower layer is aligned with the direction of absolute plate motion. Tomographic results show a featureless mantle beneath the Sachigo proto-craton and more heterogeneity towards the south end of the line. A steeply dipping slab-like feature in the lithosphere correlates with wide-angle refraction and deep-reflection seismic profiles. A similar high-velocity feature continues well into the transition zone, but its origin remains to be understood. Towards the southern end of the line there is a deep-seated low-velocity anomaly, which may be associated with Keweenawan plume activity. As a whole, the seismic results show many features that support ideas of subduction-related accretion of a thick stable Archaean tectosphere. There are, however, interesting details that are to date unique to the Western Superior Province. These include thicker than normal Archaean crust, a slab-like velocity anomaly in the mantle transition zone, and large SKS splitting in the Archaean Superior Province but little splitting in the surrounding Trans-Hudson Proterozoic shear zone.

Abstract Lithosphere that formed in Archaean and possibly early Proterozoic time is thicker, more buoyant, and geochemically distinct from lithosphere that formed after about 2.3 Ga. Mantle xenolith and seismic data indicate that some cratonic roots, or ‘keels’, extend to depths of c. 250 km, compared with normal continental lithosphere of thickness 150 km or less; yet many cratons have experienced uplift, dyking and kimberlite emplacement, suggesting interactions with hot, rising asthenosphere referred to as mantle plumes. Plumes supply additional heat to the base of the lithospheric plates, whose base can be heated and entrained in the flow (thermal erosion). How have these cratonic keels persisted despite their interactions with mantle plumes? The geometry of cratonic keels during their interactions with mantle plumes is a critical factor controlling keel preservation. To a laterally spreading plume head, cratonic keels appear as major obstacles, and the hot, buoyant plume material ponds beneath thinner lithosphere. Our model simulations show that deep keels deflect mantle plume material and that steep gradients at the lithosphere-asthenosphere boundary between Archaean keels and ‘normal’ lithosphere will focus flow, leading to localized adiabatic decompression melting. Plume processes can lead to a reduction in the breadth of a cratonic root where the plume rises beneath the craton, regardless of the initial breadth of the craton. Where the plume rises beneath a craton the hot plume material will spread laterally beneath the keel and attain thicknesses of tens of kilometres. This transfers heat to the base of the lithosphere and could generate small volumes of melt at considerable depth, depending on the composition of the lower lithosphere. We have used model simulations of plumes beneath Africa to predict the magnitude and direction of seismic anisotropy caused by lateral flow of hot plume material beneath and around a cratonic keel. The shear-wave splitting in our models is greatest at the edge of the cratonic keel, and its azimuth is parallel to the plume flow direction.