shear wave splitting in northeast venezuela trinidad and the eastern caribbean abstract we present new measurements of splitting of core shear wave phases sks skks pks recorded at seven stations in the eastern caribbean region six stations in northeastern venezuela and trinidad constituted the secasa92 temporary broadband array the seventh station sjg is installed on puerto rico we also attempted to measure source side splitting of s waves from earthquakes in the lesser antilles subducted slab but with little success mean splitting parameters fast polarization direction and delay time t at the stations are consistently east west trending and over 1 s respectively delay times at the secasa92 array are amongst the highest yet measured reaching 2 1 s and generally over 1 5 s we note an approximately linear decrease in delay times here from maxima at the northernmost stations to a minimum at the southernmost the secasa92 array is situated in the wide south america caribbean sa ca plate boundary zone and sjg is within the north america caribbean na ca plate boundary zone in both cases trends are parallel to the plate boundary zones locally we interpret these results to indicate that the sa ca plate boundary deformation extends to depths of at least 200 km and that south american crust lithospheric mantle and asthenosphere are coherently deformed the extreme delay times we measure are probably the result of transpression in the sa ca plate boundary and overriding of the lesser antilles slab by south america compressional flattening results in extremes of foliation development and dextral shearing yields a strong lineation in upper mantle olivine by contrast the delay time at sjg 1 2 s is near the global average and reflects the fact that the na ca plate boundary zone is purely sinistral or includes minor transtension we note that recent counterclockwise rotation of puerto rico does not apparently affect upper mantle deformation beneath the island we suggest that this rotation is crustal driven by far field sinistral shear and partly decoupled from the mantle shear our results are consistent with the trench parallel mantle flow hypothesis russo and silver 1994 science 263 1105 1111 and indicate that asthenospheric mantle flow beneath the caribbean plate is eastward and probably drives the plate s eastward motion

lateral heterogeneity and azimuthal anisotropy in the north atlantic determined from ss s differential travel times to explore upper mantle heterogeneity beneath the north atlantic we have measured 70 ss s differential travel times using the waveform cross correlation method both the two way preliminary reference earth model prem and j b jeffreys bullen residuals exhibit a maximum variation of 17 s ranging from 10 to 7 and 9 to 8 s respectively the lack of correlation between ss s and s residuals along with the strong correlation between ss s and ss residuals suggests that ss s differential travel times eliminate source and receiver effects and sample the upper mantle anomalies near the ss bounce points without being strongly affected by lower mantle heterogeneities we find evidence for a strong local anomaly in the north atlantic residuals for bounce points immediately to the north of the azores gibraltar plate boundary average about 4 s faster than for the bounce points to the south of the boundary we infer that this velocity contrast exists in the lithosphere and continues below a depth of 200 km a series of linear multiple regression experiments demonstrate that variation in velocity with increasing age of the seafloor and azimuthal anisotropy are both significant contributors to the prem and j b residuals assuming that travel time delays are a linear function of the square root of the seafloor age we find the coefficient for age to be about 1 0 s m y this is smaller than but not significantly different from the corresponding slope for the s delays of intraplate earthquakes in the atlantic reported by j d duschenes and s c solomon 1977 and favors the presence of small amounts of water e g 0 1 in the upper mantle in modeling anisotropy we carry out an experiment using synthetic seismograms to examine the shear wave splitting phenomena in the anisotropic medium and investigate three possible forms for sh wave delay as a function of azimuth sinusoidal variations as a function of 2 4 and a combination of 4 predominant and 2 that mimics the delay expected in olivine crystals when the axis is aligned horizontally our regression models strongly prefer the predominance of a 4 to a 2 variation and the last type of anisotropy explains our data best from the preferred anisotropy model the olivine axis or the upper mantle shearing flow is inferred to be oriented n s the magnitude of the azimuthal variation suggests that anisotropy extends over a depth range of several hundred kilometers in the upper mantle

seismic azimuthal anisotropy beneath the pakistan himalayas teleseismic s sks and skks data collected from a temporary broadband array across the himalayan front in pakistan are analyzed for shear wave splitting parameters the sks and skks phases have ray paths originating from both the south pacific and colombia which have azimuths approximately 40 apart with respect to the pakistan array if significant seismic azimuthal anisotropy is present we should observe splitting associated with one of these ray paths no evidence was seen for any shear wave splitting beneath any of the stations in the array teleseismic s waves were also used in order to provide better azimuthal coverage for the shear wave splitting measurements we were able to correct for any source side anisotropy when needed no receiver side splitting was observed in any of the s wave data the lack of shear wave splitting beneath the pakistan array indicates that there is no appreciable large scale azimuthal anisotropy beneath this part of the himalayas therefore if there is any significant strain in the upper mantle beneath this area it must either be vertically oriented or if horizontal vertically vary in such a way that the integrated effect on s wave splitting is null

mantle deformation and tectonics constraints from seismic anisotropy in the western united states abstract we have examined shear wave splitting in teleseismic shear waves sks skks s from 15 stations in the western united states based on analysis of 123 records from 67 events the varied past and present tectonic styles in this region subduction transform faulting extension and stable domains are expected to cause a wide variety of anisotropic behavior and therefore make it an excellent natural laboratory fast polarization azimuths vary from e w to nearly n s and time delays t range from being barely detectable less than 0 6 s to 1 6 s most stations yielded consistent measurements independent of station event geometry the exceptions were stations situated very close to the san andreas fault which yielded well constrained but inconsistent splitting parameters these have been successfully modeled by two anisotropic layers with different horizontal symmetry axes the upper layer has a fast direction parallel to the fault and t of about 1 s the lower layer with t from 0 6 to 0 9 s is oriented e w near the san francisco bay area and ne sw in the mojave desert other measurements of e w fast are observed as far east as western nevada with large delay times of 1 3 1 5 s but disappear to the north stations in the northwestern basin and range have values of oriented at about 70 with delay times ranging from 0 7 to 1 2 s varies in other regions of the basin and range from 70 in eastern nevada to 20 in the transition zone between the colorado plateau and the basin and range with t of 1 0 s measurements are geographically coherent and in many cases changes correspond to geologic or tectonic boundaries the e w fast feature extends across several tectonic provinces but disappears near the southern boundary of the gorda plate and may therefore be related to the dynamics of plate interactions regional variations in the fast directions of anisotropy within the basin and range suggest that the present extension has not been strong enough to align the fast directions over a great depth extent fast polarization directions within the transition zone between the basin and range and colorado plateau are nearly parallel to the geologic and geophysical boundaries between the two regions and to other geophysical boundaries the size of the delay times suggests that in many of the areas studied either the mantle is almost twice as anisotropic as we estimate or the anisotropy and therefore mantle deformation extends into the sublithospheric mantle we used the two layer method to examine the effect of a postulated lower layer with fast direction parallel to the absolute motion of the north american plate and found the best estimate of contribution from a uniform layer is 0 1 s of t although a contribution of up to 0 7 s cannot be ruled out a value as low as 0 1 s suggests either that olivine is not strongly aligned in a horizontal direction in this layer or that the layer thickness is of the order of 10 km

fabric related seismic anisotropy in upper mantle xenoliths evidence from measurements and calculations the seismic properties and their directional dependence have been measured in four mantle peridotites from picritic tuffs vitim baikal region russia geothermobarometric data suggest a source depth for the xenoliths corresponding to the spinel garnet transition zone 40 90km the experimental data include the measurement of p and s wave velocities and densities at confining pressures to 600 mpa and temperatures to 600 c and the determination of velocity anisotropy and shear wave splitting the measured velocities are compared with velocity calculations based on lattice preferred orientation lpo of olivine orthopyroxene the modal composition of the rocks and the single crystal stiffness coefficients macroscopic velocity anisotropy is strongly controlled by the lpo of olivine measured 600 mpa room temperature and calculated p wave anisotropies compare fairly well although the experimentally determined velocities are markedly lower than the calculated velocities the correspondence between measured and calculated shear wave splitting data is less pronounced this is attributed to grain boundary effects retrogressive reactions interstitial glass incomplete crack closure and to the fact that the lpo of the cpx fraction has not been included in the calculation the magnitude of p wave anisotropy corresponds to that observed in mantle peridotites exposed at the surface ivrea zone northern italy ronda southern spain

a comment on the form of the geometrical spreading equations with some numerical examples of seismic ray tracing in inhomogeneous anisotropic media summary we point out that while the equations of seismic ray geometrical spreading given recently by norris do differ as stated from equations given by cerveny it does not imply that the latter equations are wrong the two sets of equations differ only in form and in a way which in part at least can be ascribed to different choices of hamiltonian by the two authors quantitatively though the two sets of equations are entirely equivalent we also present some numerical results of ray tracing in anisotropic models simulating a continential rift a spreading ridge and a subduction zone these three structures span a range of geological mechanisms for seismic anisotropy though definitive conclusions cannot be easily drawn when there is both anisotropy and inhomogeneity the results do indicate the magnitude of ray path travel time and amplitude variations to be expected for p waves when anisotropy is introduced

sks splitting and the seismic anisotropy of the mantle beneath the hikurangi subduction zone new zealand abstract the analysis of 12 sks phases recorded on broadband stations above the hikurangi subduction zone in new zealand shows clear evidence of mantle anisotropy with the fast direction 28 5 almost parallel to the strike of subduction and the dominant geology of the region the slow shear wave delay times show a systematic change with the azimuth of the arrivals which if hexagonal symmetry is assumed indicates that either the axis of symmetry of the anisotropic volume beneath the subduction zone is not horizontal or that more than one anisotropic layer is present the magnitude of the delays 1 5 0 4s suggests that the anisotropy is most probably confined to the top 300 km of the mantle

modelling the anisotropic seismic properties of partially molten rocks found at mid ocean ridges abstract the problem of modelling the seismic properties of mid ocean ridge rocks of the axial magma chamber amc and the low velocity triangle lvt immediately below it has been addressed using two samples from oman ophiolite as examples the specimens are a layered gabbro from the lower oceanic crustal sequence and a harzburgite from the upper most mantle section at the palaeo ridge axis the seismic properties have been simulated at a temperature of 1200 c and pressure of 200 mpa so that the basalt melt is above its solidus various effective medium methods are discussed in the perspective of modelling rocks that have a strong background elastic anisotropy due to crystal preferred orientation cpo and an introduced anisotropy due to oriented melt filled inclusions calculations using various methods show a wide range of predicted melt fractions for a given seismic velocity the self consistent scheme scs and differential effective medium dem are compared in some detail a tensorial model was developed using a poro elastic method of gassman at low seismic frequency and a standard dem with isolated basalt inclusions at high frequency at low frequency the basalt pore fluid is considered to be everywhere connected and the pressure uniform whereas at high frequency basalt inclusions are isolated assuming the medium to be a standard linear solid the low and high frequency velocities were used to calculate the anisotropy of attenuation calculations with spherical basalt inclusions show that the seismic velocities decrease and attenuation increase with increasing melt fraction the symmetry of the background anisotropy due to cpo is preserved but gradually reduced with increasing melt fraction for the seismic velocities observed for the axial magma chamber amc with v p of 3 5 to 4 7 km s the model predicts 50 70 melt for the low velocity triangle lvt below the amc with v p of 4 8 5 6 km s the model predicts 35 50 melt predictions from the observed attenuations q 1 at the amc of 0 05 0 02 and in the lvt of 0 02 0 01 are around 60 and 45 melt respectively seismic anisotropy has been modelled with ellipsoidal basalt pancake shaped inclusions with their circular sections in the xy foliation plane only a 2 3 of ellipsoidal basalt inclusions are required to over print the anisotropy of the background medium there is a rapid decrease in v p and increase of q 1 for propagation normal to the foliation z with increasing axial ratio of the inclusions with increasing axial ratio the model predicts decreasing amounts of melt for given velocity v p in the z direction for the lvt only 15 25 melt is predicted for a 10 1 inclusion instead of the 35 50 of the spherical inclusion

trench parallel flow beneath the nazca plate from seismic anisotropy shear wave splitting of s and sks phases reveals the anisotropy and strain field of the mantle beneath the subducting nazca plate cocos plate and the caribbean region these observations can be used to test models of mantle flow two dimensional entrained mantle flow beneath the subducting nazca slab is not consistent with the data rather there is evidence for horizontal trench parallel flow in the mantle beneath the nazca plate along much of the andean subduction zone trench parallel flow is attributale utable to retrograde motion of the slab the decoupling of the slab and underlying mantle and a partial barrier to flow at depth resulting in lateral mantle flow beneath the slab such flow facilitates the transfer of material from the shrinking mantle reservoir beneath the pacific basin to the growing mantle reservoir beneath the atlantic basin trenchparallel flow may explain the eastward motions of the caribbean and scotia sea plates the anomalously shallow bathymetry of the eastern nazca plate the long wavelength geoid high over western south america and it may contribute to the high elevation and intense deformation of the central andes

geodetic vertical velocities affected by recent rapid changes in polar motion summary secular motion of earth s rotation pole results in large scale secular deformation of earth here we investigate the magnitude of the deformation that has resulted from the rapid motion of the rotation pole to the east since 2005 we show that geodetic gnss doris vlbi and slr estimates of vertical velocity since 2005 have been biased by up to 0 38 mm yr 1 relative to the longer term deformation pattern the largest signals occur within regions that include the u s pacific coast europe and south pacific islands where geodetic measurements provide essential measurements of tide gauge vertical movement and important constraints on models of glacial isostatic adjustment consequently geodetic vertical velocities based on recent data should not be interpreted as being identical to centennial or longer term vertical land movement since 2010 the effect is further amplified by the overprediction of the iers polar motion model relative to the ongoing secular change in pole position during this time geodetic vertical velocities based on the iers pole tide model are not just biased relative to the long termratesbutalsofromactualpost 2010earthdeformation forgeophysicalorreference frame studies seeking geodetic vertical velocities that are representative of decadal timescales where interannual variation is considered noise the correction for this non linear effect is straightforward requiring an elastic computation using a reference rate of polar motion that is linear over the timescales of interest

a study of azimuthal p residuals and shear wave splitting across the kunlun range northern tibetan plateau abstract the kunlun fault system is an east west transcurrent fault bordering the tibetan plateau that has experienced left lateral motion since 3 ma a linear array of one component and three component seismological stations was deployed across the kunlun range for 4 months p travel time residuals define three basic units along the profile the north kunlun unit the south kunlun unit and the bayan har unit p residuals are a function of azimuth p residuals are fast in the southern stations for southern azimuths and in the northern stations for northern azimuths measurements of shear wave splitting of the s scs and sks phases have been performed on deep events delay times and fast polarization directions have been computed by finding the maximum of the cross correlation between the two horizontal components rotated at different angles delay times from 0 40 s to 2 52 s remain in the range of those found by mcnamara et al 1994 j geophys res 99 13655 13665 at two sites along the same profile p residual travel times and the splitting of s waves are discussed in relation to models of the geological evolution of this region the northern tibetan plateau appears as a juxtaposition of different terranes the azimuthal p residual pattern implies that p wave anisotropy is about 3 in the lithosphere on both sides of the kunlun fault and that fast velocities are oriented outwards from the kunlun fault plane s wave splitting is observed far from the fault and implies that the lithosphere beneath northern tibet has undergone an important deformation

interpretation of sks waves using samples from the subcontinental lithosphere abstract the seismic properties of five south african kimberlite nodules have been calculated from petrofabric measurements these garnet lherzolite and harzburgite composition nodules 56 80 olivine are known from previous studies to have originated at depths of 120 170 km in the subcontinental lithosphere their in situ seismic properties have been calculated by extrapolating the elastic constants to the appropriate temperature 900 1050 c and pressure 3 0 3 5 gpa conditions the average of the five samples has a maximum s wave anisotropy v s of 3 7 which is sufficiently high to explain previously reported teleseismic sks delay times t in continental shield areas of up to 1 7 s the biggest delay time is observed when the foliation is vertical the lineation horizontal and the fastest s wave polarized parallel to the lineation i e an orientation for transcurrent motion delay times of over 1 s can only be generated by this orientation a strong mantle crust mechanical coupling is suggested for such situations the potential use of seismic anisotropy as an indicator of strain in the lithosphere has been investigated using experimental and simulated fabric data for olivine the most abundant phase 70 in the upper mantle the v p and v s seismic anisotropy increases with fabric strength and finite strain recrystallization tends to reduce fabric strength and seismic anisotropy resulting in saturation values for experimental and simulated fabrics which correspond to approximately 50 60 strain or 8 9 v s for pure olivine aggregates a survey of 25 naturally deformed peridotites of oceanic origin suggests an average maximum s wave anisotropy for the olivine component of 9 or about 8 for a lherzolite or harzburgite rock composition when the orthopyroxene component is taken into account the ophiolite samples are twice as anisotropic as the kimberlite nodules if an average anisotropy value is representative of the lithosphere then the sks delay times represent variations in anisotropic layer thickness for delay times over 1 s as the orientation is constrained to be foliation vertical and lineation horizontal the magnitude of s wave anisotropy is less sensitive than v p to variations in olivine volume fraction in range 50 100 and to deformation or fabric intensity further suggesting that s wave anisotropy is particularly apt for the determination of the anisotropic lithosphere thickness for smaller delay times there is some trade off between structural orientation and thickness

quantitative measurements of shear wave polarizations at the anza seismic network southern california implications for shear wave splitting and earthquake prediction we analyze shear wave polarizations from local earthquakes recorded by the anza network in southern california using an automated method which provides unbiased and quantitative measurements of the polarization and the duration of linear motion following the shear wave arrival the linearity interval initial shear wave particle motions are strongly aligned at four stations a feature that is not predicted by focal mechanisms the particle motion alignment is most likely caused by shear wave splitting due to anisotropy beneath these stations a result supported by the clear shear wave splitting seen in a borehole recording near one of the anza stations these results are consistent with an earlier analysis of these data by peacock et al 1988 however our analysis does not support claims by crampin et al 1990 that shear wave splitting delay times at station knw exhibit temporal variations which can be correlated with the occurrence of the north palm springs earthquake ml 5 6 of july 8 1986 automatically determined linearity intervals scatter widely from 0 02 to 0 15 s and exhibit no clear temporal trends we find a correlation between earthquake moment and the linearity interval possibly a result of longer effective source time functions for the larger events the inability to identify a distinct slow quasi shear wave pulse for the vast majority of these events indicates that scattering strongly affects the particle motion even in the very early shear wave coda analysis of earthquake clusters with similar waveforms recorded at knw shows that seismic green functions are stable throughout the observational period and that most linearity interval variation is due to source and ray path differences between events if shear wave splitting is causing the observed delay times between horizontal components the waveform stability for events in these clusters restricts any temporal changes in shear wave splitting delay times to less than 5 10

p and s wave anisotropy and shear wave splitting at pressure and temperature in possible mantle rocks and their relation to the rock fabric the seismic properties and their directional dependence have been measured in three potentially mantle derived rocks peridotite ivrea zone italy lherzolitic xenolith massif central france serpentinite the experimental data include the measurement of elastic wave velocities and densities at confining pressures up to 600 mpa and temperatures up to 600 c 600 mpa and the determination of velocity anisotropy and shear wave splitting the petrophysical data can be related to mineral composition and to microstructural characteristics of the rocks macroscopic anisotropy of p and s wave velocities is strongly controlled in the peridodites by the olivine and in the serpentinite by the serpentine antigorite lattice preferred orientation lpo respectively shear wave splitting was found to be an important characteristic of the anisotropic rocks the study shows that for the interpretation of shear wave splitting the direction of wave propagation in relation to the lpos of the major minerals and to the symmetry of their lattice fabric is of great importance in the ivrea peridotite exhibiting rhombic symmetry there is practically no shear wave splitting to be observed for s waves propagating parallel to the olivine 100 axes and the maximum shear wave splitting is not coupled with the direction of highest p wave velocity in the transversely isotropic massif central xenolith shear waves do not split when propagating parallel to the olivine 010 maximum coinciding with the direction of the radial symmetry axis of the olivine fabric the experimental data are principally in accordance with velocity calculations based on the measured olivine lpo and the olivine single crystal stiffness coefficients implications for the interpretation of seismic data are discussed

seismic studies on the trans hudson orogen of western canada abstract the discovery of diamondiferous kimberlites in the late 1980 s within the central saskatchewan segment of the trans hudson orogen prompted an aggressive diamond exploration program in the region the presence of this apparently uncharacteristic mineralization within the proterozoic orogen was initially a puzzle recent seismic reflection and geochronology evidence however suggests that the glennie domain of the orogen may be cored by a microcontinental archean block while geochemical results imply that segments of the precambrian basement of southern saskatchewan mainly comprise archean crust which has been reworked by thermotectonism we present the initial results of two seismic programs designed to elucidate the lithospheric structure associated with the economically targeted segment of the orogen preliminary interpretation of wide angle reflection data along a 730 km north south profile in north central saskatchewan reveals a very complex structure with the crustal thickness varying from a maximum of 50 km to a minimum of 37 km over a distance of 200 km beneath the glennie domain a teleseismic feasibility study with an 8 station array with individual seismographs operating for 4 to 6 months supports and extends this image of complexity relative travel time differentials across a 500 km array are as large as 700 ms while key features of these delay times correlate with the crustal structure a significant portion of the variations must be assigned to deeper levels of the lithosphere receiver function analysis for crustal structure is shown to be feasible in this environment where the phanerozoic sediment thickness varies from 0 to 2 1 km and provides evidence for large variations 7 km in crustal thickness in the southwestern area of the province sks analysis at a single station also adds to the evidence which shows rapid variations in anisotropy within the orogen a teleseismic program with seventeen 3 component stations is now underway to clarify the nature of this complex crustal and upper mantle structure

seismic anisotropy beneath the shumagin islands segment of the aleutian alaska subduction zone lithospheric strain and mantle flow in subduction zone back arcs reflect the fundamental dynamics of lithospheric subduction in this study we inferred possible strain geometries in the mantle beneath the eastern aleutians by mapping seismic anisotropy with shear wave splitting we analyzed s phases of local events recorded by the shumagin seismic network and obtained splitting parameters fast direction and splitting time t and their 95 confidence regions to avoid complexity from converted phases data were restricted to paths at incidence angles less than 32 and to events deeper than 50 km the s phases sample the forearc and back arc of the shumagin segment from 159 7 w to 163 8 w at depths up to 216 km in our data set of 84 local events 44 of 127 s phases yielded constrained splitting parameters with consistent fast directions that are approximately parallel to the volcanic arc roughly 60 the remaining 83 s phases produced null measurements and the fast directions inferred from these phase polarizations agree with the fast axis alignment in the constrained measurements splitting times vary from 0 1 to 0 35 s and show an overall increase with hypocentral depth observed splitting parameters are well matched by predicted values for olivine rich mantle wedge models with 1 sv anisotropy where the olivine b axis is arc orthogonal and the a axis is vertical no arc parallel strain or arcparallel no vertical strain these results are consistent with mantle wedge strain models in which arc normal compression is accompanied by arc parallel or vertical shearing or extension

seismic anisotropy in the fore arc region of the hikurangi subduction zone new zealand abstract shear waves from local earthquakes recorded on an array of nine broad band seismograph stations have been used to investigate the anisotropic structure within the crust of the overlying australian plate and within and below the subducting slab underlying the southern north island new zealand earthquakes near the plate interface were used to image the crustal anisotropy whereas deeper down dip subduction zone events were used to sample the subducting slab and the mantle below apart from a small amount of variable near surface anisotropy top few kilometres of the crust the measured fast shear wave directions of both the crust of the overlying plate 51 18 and the mantle 41 15 are similar and are nearly parallel to the strike of the subduction zone and the dominant geological features of the region sks measurements which mainly sample the mantle beneath the subducting plate give a value of 28 11 and measurements made elsewhere in the new zealand region using sks and scs2 produce similar results the measurements of crustal anisotropy are similar to those found further south on the wellington peninsula and confirm the presence of near surface anisotropy the maximum delay times for crustal earthquakes indicate a pervasive shear wave velocity anisotropy of 4 in the crust of the overlying plate and the increase in the delay times with depth for the subduction zone earthquakes translates to a shear wave velocity anisotropy of 1 4 in the mantle below the subducting slab these results when combined with the sks measurements suggest trench parallel flow in the mantle below the subducted slab and a coherence of the deformation processes in the crust of the overlying plate and the mantle beneath the subducting slab

on the relation between seismic anisotropy and finite strain because seismic anisotropy in the upper mantle is due primarily to lattice preferred orientation lpo of olivine crystals induced by deformation observations of anisotropy can provide constraints on mantle deformation and flow in this paper we use a previously published theory for deformation induced lpo ribe and yu 1991 to quantify the relation between seismic anisotropy and deformation for two model compositions a pure olivine dunite aggregate and a harzburgite comprising 70 olivine and 30 enstatite analytical and numerical solutions of the governing equations show that the predicted lpo and hence also the seismic anisotropy depends only on the finite strain and not on the deformation path by which that strain is produced the magnitude of the anisotropy depends only on the ratios c1 c2 and c2 c3 where c1 c2 c3 are the lengths of the principle axes of the finite strain ellipsoid and the directions of maximum and minimum compressional wave velocity vp coincide with the axes c1 and c3 respectively we calculate the percent anisotropies for compresional and shear waves propagating along the three principle axes of the finite strain ellipsoid as functions of the axial ratios c1 c2 and c2 c3 for our two model compositions however inferences of the magnitude of finite strain obtained using these results are likely to be minimum estimates because our theory neglects the effects of dynamic recrystallization finally we present as a byproduct of this study a new direct method for calculating the finite strain produced by an arbitiary deformation history

seismic constraints on dynamics of the mantle of the kaapvaal craton abstract shear wave splitting in the upper mantle of the kaapvaal craton in southern africa is measured by applying the sks technique to the records of a portable array of digital broad band seismographs the data are interpreted in terms of lattice preferred orientation in olivine fast wave polarization directions at all stations of the array are close to the absolute plate motion apm direction of southern africa since the end of the jurassic this alignment is similar to that previously reported for the north american craton and suggests that the large scale component of mantle anisotropy on both sides of the atlantic can be related to resistive drag at the base of the plates a likely depth range of the corresponding deformation is 150 400 km analysis of the mantle converted phases suggests that the high velocity mantle root of the kaapvaal craton which translates coherently with the plate resides in the same depth range the bottom of the root is found at a depth close to 380 km we conclude that the root is deformed by the recent plate motion but the deformations are not strong enough to be seen in the available tomographic models

earthquake strain rates and instantaneous relative motions within central and eastern asia summary we use the spatial distribution of moment tensors of earthquakes in this century to estimate the velocity field in asia within a eurasian reference frame in a least squares inversion strain rates on the surface of the earth are matched with continuous spline functions in order to recover the velocity gradient tensor associated with the seismic moment release in asia earthquakes account for 40 60 per cent of the expected motion of india relative to eurasia with the missing component of strain rate equivalent to about 20 mm yr 1 of n s shortening between siberia and india in this solution south china rotates counterclockwise and moves eastwards relative to siberia using rigid plate constraints we next investigate the characteristics of the complete horizontal strain field in asia that accommodates plate motions our strain rate solutions are analogous to the response of a newtonian thin viscous sheet in which the rate of work done by the straining medium in accommodating the velocity boundary conditions is a minimum in these solutions the euler pole for india relative to eurasia is constrained nuvel 1a demets et al 1994 but in the process of fitting the vlbi velocity at shanghai china ward 1994 heki et al 1995 the euler pole for south china is determined in the inversion a solution that both fits the velocity at shanghai china and yields a strain rate field consistent with the earthquake mechanisms is one where the south china block has a motion relative to siberia described by the pole at 51 n 131 e 0 3 deg myr 1 comparison of the complete strain field that accommodates plate motion with the seismic strains indicates that earthquake moment release rates in this century within mongolia are about a factor of 4 larger than the long term rate within gansu ningxia the earthquake moment rates have been about a factor of 2 higher than the long term rate the strike slip faulting within mongolia gansu ningxia western sichuan and yunnan is possibly a direct result of velocity boundary conditions imposed on the south china block by forces unrelated to continental collision of india and eurasia such as forces associated with subduction along the margins of south eastern asia verification of this requires a better understanding of the role of pre existing zones of weakness within the asian continental lithosphere

the strength of mantle silicates at high pressures and room temperature implications for the viscosity of the mantle the rheological properties of the upper mantle important in convective and tectonic processes have long been studied through low pressure creep experiments on olivine1 5 in comparison our understanding of deformation in the deep earth is incomplete because there are few experimental constraints on the rheology of the transition zone 400 670 km depth and the lower mantle estimates of the viscosity at these depths have had to rely on geophysical models of for example post glacial uplift the geoid and changes in the earth s rate of rotation6 13 recently we developed techniques for measuring solid state creep at mantle pressures14 and here we use these methods to make direct laboratory measurements of the strength of key mantle minerals we study mg fe 2sio4 olivine representative of about 60 of the upper mantle mg fe 2sio4 spinel which constitutes a large fraction of the transition zone and mg fe sio3 perovskite and perovskite mg fe o magnesiowiistite assemblages which comprise most of the lower mantle all at room temperature and pressures of up to 60 gpa our results suggest that the transition zone may form a layer of relatively high viscosity or strength between the upper and lower mantle a simple theological model shows that this behaviour may be compatible with current geophysical observations

a student s garden of anisotropy anisotropy is seismology s most venerable skeleton in the closet of all elastic propagation problems none eludes our intuition more than ani sotropy even the formidable problem of propagation in a three dimen sional heterogeneous medium seems less intimidating because we have an intuitive sense and some computational tools for the behavior of classical rays and even of diffractions in isotropic media this intuition is precious to us it is a product of centuries of study and use of wave propagation and is for seismologists the most valuable heritage from mathematical physics many significant advances in seismology especially the inference of gross earth structure from travel time analysis depended on geometric ray theory in isotropic media no comparable progress has been made in the study of anisotropy in the earth except where simple travel time differences could be predicted the simplest of these observations is the azimuthal dependence of the apparent velocities of pn refractions in the upper mantle under oceans e g raitt et al 1969 bibee shor 1976 and under continents e g vetter minster 1981 the orientation of the pn anisotropy in the oceanic upper mantle clearly indicates a tectonic origin using multicomponent seismic data we can observe shear wave birefrin gence in body wave arrivals crampin 1984 1985 and surface wave propagation nataf et al 1984 tanimoto anderson 1985 lerner lam 1985 all of which are travel time phase differences recent work in marine seismology gives convincing evidence for anisotropy in the basaltic upper oceanic crust and places some constraint on its variability with depth stephen 1985 used the time dependence of the polarization of shear wave arrivals at a borehole receiver deep sea drilling project dsdp site 504b in 6 myr old crust to show that oriented cracking in

role of melt during deformation in the deep crust deformation in the deep crust is strongly influenced by the presence of melt injected melt or magma weakens the crust because strain will tend to localize where melt is present the amount of strain a pluton may accommodate is dependent on the length of time it takes for a pluton to crystallize and the strain rate for plutons that intrude into rocks which are near the solidus temperature of the melt crystallization times can be quite long 1myr partial melting of deep crustal rocks can lead to melt enhanced embrittlement this occurs because the volume change for most melting reactions is positive therefore when the rate of melt production outpaces the rate at which melt can leave the system the melt pressure increases eventually the melt pressure may become sufficiently high that the melting rocks behave in a brittle fashion and fracture conjugate sets of dilatant shear fractures filled with melt occur in migmatite from the central gneiss belt canada this suggests that melt enhanced embrittlement occurred in these rocks an expression which relates the magnitude of differential stress to the angle between conjugate dilatant shear fractures is derived assuming that migmatite has a small tensile strength differential stresses are 20 mpa in migmatitic rocks at the time melt enhanced embrittlement occurs the occurrence of melt enhanced embrittlement shows that a switch in deformation mechanism from plastic flow to cataclasis is possible in the deep crust during melting furthermore repeated episodes of melt enhanced embrittlement in migmatitic rocks may be an efficient mechanism for extracting melt from partially melted terrains

variations in sks splitting across western canada measurements of sks splitting parameters have been made at stations of the canadian national seismograph network for events over the period 1993 mid1994 here we report on results obtained from western stations in the northern cordillera the alignment of polarization directions with major transcurrent fault systems suggests that the crust and upper mantle are being deformed coherently splitting parameters in the southern cordillera however are oriented approximately perpendicular to the strike of major tectonic features notably the cascadia subduction zone this contrast in the character of fast polarization directions between the southern and northern cordillera is thought to reflect the change in the tectonic style of the adjacent plate margin from convergent along the southern b c coast to transform north of the triple junction one plausible explanation is that upper mantle anisotropy in the southern region is due primarily to mantle flow perpendicular to the margin involving upwelling at the juan de fuca ridge subduction along the coast and associated back arc asthenospheric flow beneath the b c interior

shear wave splitting around the northern atlantic frozen pangaean lithospheric anisotropy abstract although the number of teleseismic shear wave splitting measurements has considerably increased world wide during the last decade there is no consensus on a unique and universal interpretation of these data splitting might primarily originate in the lithospheric and or asthenospheric mantle but the technique has no vertical resolution and therefore leaves this question of anisotropy location unresolved generally the consequences of the two possibilities are used to infer which is more likely if splitting occurs in the lithosphere the splitting parameters the orientation of the fast split shear wave o and the delay time t should be related to frozen or active deformation on the other hand if the splitting occurs in the asthenosphere the observed splitting parameters should be related to the present day flow induced by the differential motion between the tectonic plates and the underlying mantle the aim of this paper is to test these hypotheses against shear wave splitting observations from western europe and the eastern u s indeed except in a few places splitting on either side of the atlantic on hercynian caledonian or grenvillian structures do not appear to have been deeply and pervasively affected by the opening of this ocean this provides us the opportunity to place these results in the pangaean reference frame of the hercynian caledonian belt before the opening of the atlantic in many places the close parallelism between lithospheric structure and the absolute plate motion apm impedes a reliable identification of the origin of the anisotropy especially an asthenospheric flow channelled around the deeply rooted north american craton may be compatible with the large scale anisotropy pattern in the eastern u s despite this ambiguity fast polarization directions o observed on or near some major lithospheric structures correlate better with the trend of the local tectonic structures than with the expected asthenospheric flow therefore we suggest that these lithospheric structures developed before the opening of the atlantic remain frozen since hercynian caledonian or grenvillian times and are responsible for the observed splitting at these places

anisotropic upper mantle stratigraphy and architecture of the slave craton the earth s physical properties show a dominantly radial structure which is the result of compositional differentiation isochemical phase changes1 and rheological layering2 rheological layering is perhaps the most difficult to investigate using conventional seismological techniques because the seismic manifestation of this property elastic anisotropy may closely mimic the effects of isotropic heterogeneity3 nonetheless an improved characterization of earth s rheological structure promises important insights into such processes as plate dynamics and continental evolution here i present a methodology for effectively characterizing sharp transitions in anisotropic upper mantle structure using the coda of teleseismic p waves application to seismograms from the slave craton reveals a well developed stratigraphy at least in part anisotropic with major boundaries occurring at nominal depths of 75 135 and 195 km the geometry and sharpness of these discontinuities suggest a structural origin perhaps involving shallow subduction

sks splitting beneath continental rift zones we present measurements of sks splitting at 28 digital seismic stations and 35 analog stations in the baikal rift zone siberia and adjacent areas and at 17 stations in the east african rift in kenya and compare them with previous measurements from the rio grande rift of north america fast directions in the inner region of the baikal rift zone are distributed in two orthogonal directions ne and nw approximately parallel and perpendicular to the ne strike of the rift in the adjacent siberian platform and northern mongolian fold belt only the rift orthogonal fast direction is observed in southcentral mongolia the dominant fast direction changes to rift parallel again although a small number of measurements are still rift orthogonal for the axial zones of the east african and rio grande rifts fast directions are oriented on average nne that is rotated clockwise from the n s trending rift all three rifts are underlain by low velocity upper mantle as determined from teleseismic tomography rift related mantle flow provides a plausible interpretation for the rift orthogonal fast directions the rift parallel fast directions near the rift axes can be interpreted by oriented magmatic cracks in the mantle or small scale mantle convection with rift parallel flow the agreement between stress estimates and corresponding crack orientations lends some weight to the suggestion that the rift parallel fast directions are caused by oriented magmatic cracks

seismic anisotropy in the mantle beneath an oceanic spreading centre beneath an active mid ocean ridge the mantle upwells in response to the divergence of the newly formed plates leading to high temperatures and pressure release melting below the ridge axis the width of the upwelling region and the amount of melting depend on mantle rheology1 5 but all models predict a maximum decrease in seismic velocity at the ridge axis it has also been suggested however that the alignment of anisotropic minerals by shear in the upwelling mantle will increase seismic velocity for rays travelling subvertically through the upwelling zone6 7 here we report the observation of a consistent pattern of anomalously early p wave arrival times at an array of ocean bottom seismographs deployed across the axis of the southern mid atlantic ridge p waves from distant earthquakes arrive earlier at stations near the axis than at those further away our results are consistent with a model of anisotropy in which the degree of mineral alignment is greatest directly beneath the ridge axis and significant anisotropy extends tens of kilometres from the axis

seismological evidence for three dimensional melt migration beneath the east pacific rise the extent to which crustal processes along mid ocean ridges are controlled by either the pattern of mantle upwelling or the mode of magma injection into the crust is not known models of mantle upwelling vary from two dimensional passive flow1 to three dimensional diapiric flow2 3 4 similarly beneath a ridge segment bounded by tectonic offsets crustal magma chambers may be replenished continuously along the ridge5 6 7 or at a central injection zone2 3 4 from which magma migrates towards the segment s ends here we present tomographic images that reveal the seismic structure and anisotropy of the uppermost mantle beneath the east pacific rise the anisotropy is consistent with two dimensional mantle flow diverging from the rise whereas the anomalous isotropic structure requires a three dimensional but continuous distribution of melt near the crust mantle interface our results indicate that crustal magma chambers are replenished at closely spaced intervals along axis and that crustal systems inherit characteristics of scale from melt transport processes originating in the mantle

an automated analytical method to determine shear wave splitting abstract we introduce and discuss an automated analytical method for determination of shear wave splitting to determine the polarization direction of the first arriving shear wave fast wave we evaluate the linearity of the particle motion of shear waves by maximizing as a function of azimuth the aspect ratio i e the ratio of the projections of the particle displacement on orthogonal axes in the horizontal plane in a time window that contains only the first shear wave arrival the azimuth at which the maximum aspect ratio occurs is taken as the polarization direction of the fast shear wave then the time separation between the split shear waves is determined by time shifting the component with the slow shear wave arrival toward the fast shear wave arrival until maximum linearity in particle motion is observed the polarization direction of the fast shear wave can be accurately predicted for seismograms having good signal to noise ratios as are typical of local earthquake recordings similar procedures are used to evaluate particle motion in a vertical plane in the time interval including the first shear wave arrival to examine the apparent incident angle in order to reveal the likelihood of interfering s p conversions

azimuthal anisotropy in the earth from observations of sks at geoscope and nars broadband stations abstract we have analyzed observations of low frequency sks waves recorded on horizontal components at several broadband stations of the geoscope and nars arrays splitting of sks is consistently observed and in good agreement with a hypothesis of azimuthal anisotropy due to transverse anisotropy with a horizontal axis in the upper mantle beneath the receivers we have designed an inversion method which takes advantage of observations made simultaneously at different azimuths and which permits us to obtain estimates of the direction of fast axis of anisotropy as well as the time delay between the orthogonal quasi shear waves in general the directions of fast velocity obtained both for stations around the atlantic and around the pacific appear to be related to directions of past and present plate motions

seismic anisotropy and mantle flow seismic anisotropy in the mantle is primarily due to deformation induced lattice preferred orientation lpo of olivine crystals the aim of this paper is to better understand how such lpo is produced by flow and to learn whether observations of seismic anisotropy can be used to infer the local direction of flow within the mantle the basis of the work is a two dimensional theory which describes how lpo evolves during an arbitrary time dependent deformation ribe 1989 i investigate the evolution of lpo beneath ocean ridges above subducted slabs and in thermal convective flows the ocean ridge model predicts that the seismic anisotropy produced at a fast spreading ridge should be greater in magnitude than that produced at a slow spreading ridge and that the fast axis should be more nearly horizontal the subduction zone model predicts that strong seismic anisotropy should be produced above subducted slabs with the fast axis aligned with the slab regardless of the subduction angle the evolution of lpo in typical convective flows is extremely complex and depends strongly on such factors as the mode of heating and the surface boundary conditions in general there is no simple relation between the lpo and the local flow direction in any of the models considered the attempt to use seismic anisotropy to map flow patterns on a global scale is therefore probably unwarranted however the lpo does become aligned with the flow direction when the flow is a progressive simple shear such a flow occurs in limited portions of the mantle beneath the lithosphere at a fast spreading ridge above subducted slabs and in boundary layers beneath rigid surfaces such as continents observations of seismic anisotropy may therefore indicate the sense of shear at shallow depths in the mantle

seismic anisotropy beneath the lower half of the north island new zealand teleseismic scs and sks events recorded on nine broadband seismograph stations have been used to investigate seismic anisotropy beneath the lower half of the north island new zealand this area lies above the hikurangi subduction zone and the array provides ray paths which sample the mantle both above and below the slab shear wave splitting measurements give similar fast polarizations and delay times at each station the average sks fast polarization is approximately ne sw subparallel to the strike of subduction and the major geological features with an average sks delay time of 1 6 0 1 s this lack of variation in splitting parameters suggests that similar fast polarizations are found in both the mantle wedge and the subslab mantle the anisotropy in the lithospheric portion of the mantle wedge is most likely caused by the preferred orientation of olivine due to the shear deformation associated with oblique convergence any anisotropy in the slab is probably due to fossil mineral alignment anisotropy in the asthenosphere is most likely caused by the preferred orientation of olivine due to asthenospheric flow the similar ne sw fast polarizations found in the asthenosphere both above and below the slab suggest that the mantle flow is in a trench parallel direction in both regions

shear wave splitting in the appalachians and the pyrenees importance of the inherited tectonic fabric of the lithosphere splitting of teleseismic shear waves has been measured in the appalachians eastern usa and the pyrenees western europe using data recorded by permanent and portable stations from a comparison of the results it appears that an interpretation of the recorded seismic anisotropy in terms of geodynamics is not straightforward successive geodynamic events have generated structures that may have resulted in a similar pattern of mantle flow and that therefore may have contributed in the development of the recorded anisotropy combining geological and geophysical arguments it appears that the mantle anisotropy measured across the appalachians and the pyrenees may not be systematically appalachian or pyrenean in age but may be mainly due to a lithospheric structure formed during earlier major tectonic events i e the grenvillian and the hercynian orogenies respectively we suggest that during major episodes of continent assembly a pervasive tectonic fabric is developed in the lithospheric mantle in the subsequent evolution of the continent this fabric may induce a significant mechanical anisotropy that will drastically influence the mechanical behaviour of the iithosphere when submitted to new tectonic events

mantle anisotropy beneath northwest pacific subduction zones to assess the location strength and orientation of seismic anisotropy in the southern kuril japan and izu bonin subduction zones we analyzed shear wave splitting in local s phases and teleseismic sks phases fast directions from these phases are roughly parallel to the direction of absolute pacific plate motion wnw in izu bonin roughly parallel to the strike of the trench nne in central honshu and roughly n in the southern kuril back arc in the vicinity of sakhalin island assuming that the orientation and strength of anisotropy does not vary with depth modeling of splitting times from local s and teleseismic phases recorded at sakhalin island requires that the maximum depth of anisotropy lies between 480 km and 950 km in contrast splitting times for local s phases from events in izu bonin rule out anisotropy in the transition zone all the data are consistent with a model in which the lower transition zone 520 660 km and lower mantle are largely isotropic and in which anisotropy occurs intermittently in the upper transition zone 410 520 km possibly due to preferred orientation of spinel assuming that splitting in the upper mantle is produced by preferred orientation of olivine the observed fast directions indicate that the geometry of back arc strain varies systematically between subduction zones the relationship of fast directions and plate motions suggests that the subducting slab exerts significant control on back arc flow but that flow correlated with upper plate deformation must also exist

contrasting lithospheric structure between the colorado plateau and great basin initial results from colorado plateau great basin passcal experiment the colorado plateau great basin cpgb passcal experiment provides broadband seismic data that allow us to determine crustal and upper mantle structure in the colorado plateau and eastern great basin twelve portable broadband seismometers throughout utah and nevada are supplemented by permanent broadband stations in the region receiver functions indicate an average crustal thickness of 35 5 km in the great basin gb with a strong p to s ps conversion from the moho the crust in the transition between the great basin and colorado plateau cp produces a strong intracrustal conversion the cp has somewhat thicker crust 42 5 km and weaker moho conversion br surface wave phase velocities are 3 4 higher than cp phase velocities in the 20 35 s range consistent with a thinner crust under the gb the phase velocities at periods greater than 40 s are indistinguishable between the two provinces indicating that differences between upper mantle velocities are minor shear wave splitting yields fast polarizations that parallel the western and northern boundaries of the cp with delay times decreasing towards the center of the plateau

compression along the northern apennines evidence from the mw 5 3 monghidoro earthquake in this paper we present the seismological data recorded during the deployment of a dense three component seismic network installed a few hours after the 2003 mw 5 3 monghidoro earthquake in northern apennines the main shock focal solutions derived from polarities distribution and body wave modelling of regional broadband data show a ne sw striking reverse mechanism accurate relative locations of aftershocks and the inversion of focal mechanisms show that earthquakes occurred on a nw dipping backthrust within the adria lithosphere under a nw trending horizontal compression the observed compression is a secondary process possibly explained by differential motion within the adriatic lithosphere fault geometry and kinematics is controlled by pre existing structures

shear wave splitting beneath western united states in relation to plate tectonics we have examined shear wave splitting in teleseismic shear waves from 26 broadband stations in the western united states fast polarization directions and delay times t show spatial variations that are coherent within geologic provinces stations located near the san andreas fault show clear evidence for fault parallel anisotropy in the crust and upper mantle 115 125 km thickness this can be explained by the finite strain associated with the relative plate motion between the north american and pacific plates the lateral extent of this strain field is probably narrow to the west because stations 55 km west of the san andreas fault do not show fault parallel anisotropy in southern california station lac located 80 km east of the san andreas fault shows large fault parallel anisotropy this suggests that the pacific north american plate boundary in the mantle might be displaced to the east in southern california a deeper e w oriented fast direction of anisotropy underlies the fault parallel anisotropic layer in the vicinity of the san andreas fault an e w fast feature is also present beneath the western basin and range and the foothills of the sierra nevada although local variations are present the magnitude of delay times suggests that this feature resides in the asthenosphere we interpret this feature as the asthenospheric flow in the slabless window left behind the farallon plate the flow induced anisotropy may partially be frozen in at shallow depths station orv is located near the southern edge of the gorda slab where no anisotropy is detected the absence of anisotropy at this location could therefore mark a boundary between farallon associated flow and regions where e w oriented asthenospheric flow did not occur the lack of evidence for ne sw fast orientation within the walker lane shear belt of western nevada suggests that this crustal feature does not extend into the mantle or that is not as well developed as that beneath the san andreas fault stations located over the young subducting gorda plate mark a change in the fast direction to nearly ne sw this direction aligns well with the maximum compressive stress direction in the overlying north american plate and the ne sw directed internal shearing of the gorda plate the anisotropic thicknesses calculated from delay times suggest roughly double that expected for purely lithospheric contributions this implies that the anisotropic thickness may include some of the asthenosphere alternatively using a higher anisotropy of 8 can bring thicknesses in line with other measures of lithospheric thicknesses the correspondence between the fast directions and the present plate tectonic deformations suggest that mapping upper mantle deformation through seismic anisotropy is a viable method and that asthenospheric flow may be a significant contributor to seismic anisotropy

lateral variations in lowermost mantle shear wave anisotropy beneath the north pacific and alaska s waves recorded by long period world wide standardized seismograph network and broadband stations in north america for deep northwest pacific subduction zone earthquakes provide evidence for anisotropy in the lowermost mantle shear velocity structure beneath the north pacific and alaska systematic delays of up to 4 s are observed between longitudinal components sv and transverse components sh of motion for core reflected scs waves as well as for core grazing and diffracted s waves the absence of significant splitting for s waves that have turning points more than a few hundred kilometers above the core mantle boundary indicates that anisotropy is localized within the d region the lowermost portion of the mantle sv sh differential arrival times for both scs and sdiff along with path length estimates assuming a 250 km thick d region indicate spatial variations in the strength of shear wave anisotropy the strongest anisotropy 1 1 5 is found in the eastern part of the study area with systematic reduction in magnitude toward the west a transverse isotropy model can explain the data with the velocity structure for horizontally polarized waves vsh having a 2 3 discontinuous shear velocity increase at the top of d as proposed in earlier studies of the region and a similar structure for s wave particle motion in the direction normal to the core mantle boundary vsv but with the velocity jump at the top of d and the velocity within d being reduced from that for vsh by 0 5 1 5 km s large uncertainties exist for velocity gradients above and below the velocity jump but the requirement of a reduced vsv relative to vsh in d is clear synthetic waveforms calculated by using separate isotropic structures for sh and sv match the observations well and constrain the basic anisotropic structure because the shear waves all traverse the region with near grazing geometries the study area exhibits strong lateral variations in lower mantle shear velocity structure and variable thickness of the d layer topography of the d layer is not well resolved because of trade off with volumetric heterogeneity

seismic anisotropy above a subducting plate shear wave splitting observed in northeastern colombia has provided evidence of seismic anisotropy in a shear zone immediately above a subducting plate in an upper mantle composed mainly of olivine 57 and orthopyroxene 17 the splitting can be interpreted by wave propagation in an anisotropic medium of orthorhombic symmetry that results from alignment of these intrinsically anisotropic minerals the mechanism of alignment is most likely the shearing associated with the subduction aided by fluids migrating from the subducting plate when the plate exceeds 100 km in depth