shear wave polarizations near the north anatolian fault ii interpretation in terms of crack induced anisotropy summary the polarizations of shear waves recorded by networks of digital thxee component seismometers immediately above small earthquakes near the north anatolian fault in turkey display shear wave splitting on almost all shear wave seismograms recorded within the shear wave window this splitting is incompatible with source radiation patterns propagating through simple isotropic structures but is compatible with effective anisotropy of the internal structure of the rock along the ray paths this paper interprets the phenomena in terms of widespread crack induced anisotropy distributions of stress induced cracks model many features of the observations and synthetic polarization diagrams calculated for propagation through simulated cracked rock are similar to the observed patterns this evidence for widespread crack induced anisotropy lends strong support to the hypothesis of extensive dilatancy anisotropy eda suggested by laboratory experiments in subcritical crack growth the crucial evidence confirming some form of eda would be observations of temporal changes in shear wave splitting as the stress field alters the crack density and crack geometry there is some weak evidence for such temporal changes at one site but further analysis of suitable digital three component seismometer networks in seismic areas is required to confirm eda

body wave traveltimes and a spherically symmetric p and s wave velocity model summary the derivation of a spherically symmetric reference earth model is a subject which has recently regained interest after a decade of studies on lateral heterogeneity we now have the tools to approach the problem an iaspei working group on traveltime tables was formed in 1988 with the goal of deriving a reliable and up to date radial model of seismic velocities following the guidelines developed by that working group we have investigated the body wave traveltimes for several branches frequently reported in the isc bulletins and propose here a model for p and s wave velocities consistent with these data approximately 16000 best recorded shallow focus earthquakes are selected from 24 years of isc bulletins 1964 1987 and relocated using among others corrections for lateral heterogeneity from this dataset we derive summary traveltime curves specified for each degree of epicentral distance for p pcp p pkpab pkpbc pkikp and s to constrain velocity at the top of the outer core we include data for sks and skks sks from hales roberts 1970 1971 our discrete time distance curves are obtained by evaluating the shape of the residual density distribution independently for each 1 interval to limit contamination due to the interference of other branches precursors or late readings in the upper mantle the resulting model sp6 differs from iasp91 kennett engdahl 1991 only by slightly higher p and s velocities between 410 and 660 km in the lower mantle new velocity profiles have lower velocity gradients this is also true in comparison with prem dziewonski anderson 1981 at the top of the outer core vp is smaller than in prem and the velocity jump at the inner core boundary is reduced to the value of 0 62 km s 1 the new inner core radius is 1215 km our model predicts traveltimes that differ from the iasp91 tables considerably although generally less than 1 s except for sks and skks for which the differences are greater

anisotropy and the splitting of ps waves the ps wave is radially polarized in an isotropic spherically symmetric earth owing to the p to sv conversion at the free surface if we choose stations with weak sks splitting to minimize the effect of the receiver we can measure the anisotropy beneath the bounce point of the ps wave we find that with a deep source and epicentral range of 90 125 the ps phase can be distinguished from other phases owing to the limitation of recent large deep events with epicentral range between 90 and 125 we analyze data from stations in north america and events beneath the bonin islands and the fiji islands for events in these two regions the bounce points of the ps wave are roughly located northwest of the kuril islands and near the nova canton trough in the central pacific respectively because the first fresnel zone for the ps wave at this distance range is roughly 700 km across each ps splitting observation must be interpreted as a local average of anisotropic properties the average direction of the fast axis of anisotropy below the nova canton trough region is n104 2 e with a standard deviation 12 7 and a typical delay time of 1 8 s landward of the kuril islands the average direction of the fast axis is n127 6 e clockwise from north and the standard deviation is 11 3 with a typical delay time of 1 2 s we have compared the fast axis with the direction of fossil seafloor spreading and with the present day absolute plate motion the fast axis of anisotropy at the nova canton trough agrees well with the direction of absolute plate motion and poorly with the complex fossil spreading pattern this suggests that significant upper mantle shear anisotropy exists in this part of the central pacific and is consistent with mineral alignment owing to strain associated with present day plate motion similarly the fast axis beneath the kuril islands parallels the convergent motion of the slab coupled mode synthetics constructed with the strong born approximation can be used to model the interaction of the ps waves with anisotropic structure that has a horizontal axis of symmetry we use all free oscillations up to 35 mhz in the calculation the coupled mode synthetics show shear wave splitting with delay time and fast axis consistent with the ray theory prediction

a search for split sks waveforms in the north atlantic summary to test a prediction of the anisotropic properties of the upper mantle beneath the central north atlantic we measure splitting of the waveforms of sks phases recorded at wwssn stations in the azores pda and on bermuda bec because the analysis of individual analogue records does not provide sufficient resolution we have developed a new technique for representing the confidence limits that permits a straightforward estimation of the average splitting beneath a station from a group of events if events from a variety of different azimuths are recorded the average splitting parameters can be tightly constrained we found five records for pda and eight for bec that have good signal to noise ratio and a prominent sks phase only one of the phases at pda shows evidence of splitting the average anisotropy must be small if it exists at all at bec the average anisotropy is significant with a fast direction between nw se and n s and a time delay between the two split components of about 1 s both of these observations are significantly smaller than predicted by a model of the upper mantle anisotropy based on ss s differential traveltimes the absence of systematic anisotropy at pda may be attributable to a complex flow pattern in the upper mantle at this site near the eurasian african plate boundary in the vicinity of an active hotspot at bec the observed direction is consistent with the prediction and both the direction and the magnitude of the average anisotropy are indicative of a sublithospheric source to check the reliability of these results from analogue records we have performed a similar analysis for station kip on oahu where digital records are available for comparison the results from digital and analogue records are similar at kip the splitting is more variable than expected for uniform anisotropy with the fast directions for individual records ranging from about n26 e to about n83 e and with variable time delays similar variability is observed at all three stations the average anisotropy at kip is well constrained however with a time delay of about 0 6 s and fast direction n75 e the fast direction coincides with the fossil sea floor spreading direction beneath the hawaiian islands suggesting that the dominant source of anisotropy beneath kip is lithospheric

shear wave splitting and subcontinental mantle deformation splitting observations are interpreted in terms of the strain induced lattice preferred orientation of mantle minerals especially olivine we consider three hypotheses concerning the origin of the continental anlsotropy 1 strain associated with absolute plate motion as in the oceanic upper mantle 2 crustal stress and 3 the past and present internal deformation of the subcontinental upper mantle by tectonic episodes it is found that the last hypothesis is the most successful namely that the most recent significant episode of internal deformation appears to be the best predictor of q for stable continental regions this is interpreted as fossil anisotropy whereas for presently active regions such as alaska the anisotropy reflects present day tectonic activity in the stable portion of north america there is a good correlation between delay time and lithospheric thickness this is consistent with the anisotropy being localized in the subcontinental lithosphere and suggests that intrinsic anisotropy is approximately constant the acceptance of this hypothesis has several implications for subcontinental mantle deformation first it argues for coherent deformation of the continental lithosphere crust and mantle during orogenies this implies that the anisotropic portion of the lithosphere was present since the deformational episode and rules out the addition of undeformed material to this layer by subsequent underplating or conductive growth of the thermal boundary layer one of the most important issues in the study of orogenies is the need to reconcile the formation of thickened lithosphere with the paradoxically high mantle temperatures often associated with orogenic episodes most efforts to date have focussed on modes of deformation whereby the cold lithospheric mantle is removed by convective instability or delamination and replaced by warm asthenosphere these models however are incompatible with the evidence for preserved coherent lithospheric deformation rather the deformed mantle appears to have been heated in place we suggest that the elevated mantle temperatures may be due to the strain heating accompanying the deformation

seismic anisotropy tracing plate dynamics in the mantle elastic anisotropy is present where the speed of a seismic wave depends on its direction in earth9s mantle elastic anisotropy is induced by minerals that are preferentially oriented in a directional flow or deformation earthquakes generate two seismic wave types compressional p and shear s waves whose coupling in anisotropic rocks leads to scattering birefringence and waves with hybrid polarizations this varied behavior is helping geophysicists explore rock textures within earth9s mantle and crust map present day upper mantle convection and study the formation of lithospheric plates and the accretion of continents in earth history

kinematic model of active deformation in central asia the velocity field of present day deformation in central asia is modelled using a set of four rotating blocks siberia tarim tibet india on a spherical earth a best fit is inverted on the basis of estimated shortening rates across the main thrust zones himalayas tien shan and measured slip rates along the principal strike slip faults altyn tagh and karakorum separating those blocks the fit to the data implies that nearly all the present convergence between india and asia can be accounted for by slip partitioning on these four zones with as much as 50 absorbed by northeastwards extrusion of tibet this suggests that localised deformation governs the present mechanical behaviour of the central asian lithosphere

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

degrees 2 4 6 inferred from seismic tomography lateral heterogneities of geophysical fields gravity magnetic seismic velocities are usually expended into spherical harmonics free oscillation and geoid data show that heterogeneity in the deep earth s mantle is dominated by degree 2 however its geographical pattern and its location at depth is still questionable recent tomographic models of seismic velocity anisotropy and anelasticity obtained from geoscope and gdsn data make it possible to gain more insight on this problem they show that at depth larger than 400km 100km the degree 2 and t a less extent the degree 6 arise as the most important features these new models also display a large degree 4 for radial anisotropy at large depths a simple flow pattern can explain the predominance of these different degrees 2 4 6 and it is shown that degrees 2 4 6 in the transition zone are not independent they are compared to the corresponding degrees of the hotspot distribution and geoid a possible relation between the degree 2 of velocity and the degree 4 of radial anisotropy is also displayed however the different location of degree 2 in the upper mantle and the lower mantle is still puzzling

why fast polarization directions of sks seismic waves are parallel to mountain belts abstract the splitting of sks waves shows that the fast velocity is commonly parallel to the trend of mountain belts this has been related to the anisotropy induced by frozen mantle flow at lithospheric depth when the flow direction is parallel to mountain belts this flow reflects a transcurrent motion parallel to the boundaries of the plates whose collision is responsible for mountain building considering the width of major shear zones transcurrent flow should not reorient the mantle fabric on a transverse width larger than 20 40 km unless the flow is partitioned between several major transcurrent shear zones therefore such a scale seems to be too small to explain the sks results another contribution to the observed mantle anisotropy is suggested to be the fabric induced by mantle flow during a prior event of rifting at the colliding plate boundaries which are derived from passive margins

models of seismic anisotropy in the deep continental lithosphere seismological observations sks wave polarizations systematic p residual variations azimuthal dependence of pn and surface wave velocities or a dispersion of surface waves are not consistent with isotropic if laterally heterogeneous upper mantle structure therefore an anisotropy should be considered as an a priori aspect of future large scale studies of mantle structure most studies of anisotropy however have assumed horizontal or vertical axes of symmetry but such orientations cannot explain bipolar patterns of spatial variations of p residuals which we have observed at many seismological stations on the basis of the petrophysical properties of real upper mantle rocks we consider anisotropy formed either by hexagonal or by orthorhombic aggregates composed of olivine orthopyroxene and clinopyroxene rotations of the aggregates about vertical and horizontal axes allow us to find the three dimensional orientations of symmetry axes that fit combinations of both p and s seismological observations in central europe and in western north america the orientations with plunging symmetry axes velocity extremes seem to be consistent across large spatially uniform tectonic units and change abruptly at important suture zones

shear wave splitting in the upper mantle wedge above the tonga subduction zone summary intermediate period seismograms of shear waves from deep earthquakes are examined for shear wave splitting diagnostic of seismic velocity anisotropy in the upper mantle wedge between the subducting tonga slab and stations in the fiji islands earthquakes with near vertical incidence angles are selected to minimize complications of converted phases and the free surface horizontal component seismograms for 25 deep earthquakes recorded at station lak on the lau ridge show splitting of up to 0 9 s with an average of 0 45 0 25 s between fast and slow s wave components anisotropy in the upper mantle is the most likely cause of the splitting the ese orientation of the fast velocity direction at station lak corresponds to an axis of symmetry of the anisotropic mantle and is approximately perpendicular to magnetic lineations in the lau basin and the northeastern portion of the south iiji basin the observations are consistent with a model of flow alignment of olivine 100 and orthopyroxene 001 axes parallel to the direction of basin extension seismograms for eight deep earthquakes at station sva on viti levu show up to 1 0 s splitting with an average of 0 65 0 29 s and a weak nne preferred orientation of fast velocity polarizations of shear wave particle motions corrected for anisotropy are within 36 5 26 of polarixations predicted for nine moment tensor solutions the magnitude of shearwave splitting and estimates of the thickness of the lithosphere in the fiji region imply that the anisotropy cannot be contained in the lithosphere alone either the anisotropy is distributed in both the lithosphere and asthenosphere or the asthenosphere alone is anisotropic

seismic anisotropy and mantle creep in young orogens summary seismic anisotropy provides evidence for the physical state and tectonic evolution of the lithosphere we discuss the origin of anisotropy at various depths and relate it to tectonic stress geotherms and rheology the anisotropy of the uppermost mantle is controlled by the orthorhombic mineral olivine and may result from ductile deformation dynamic recrystallization or annealing anisotropy beneath young orogens has been measured for the seismic phase pn that propagates in the uppermost mantle this anisotropy is interpreted as being caused by deformation during the most recent thermotectonic event and thus provides information on the process of mountain building whereas tectonic stress and many structural features in the upper crust are usually orientated perpendicular to the structural axis of mountain belts pn anisotropy is aligned parallel to the structural axis we interpret this to indicate mountainparallel ductile i e creeping deformation in the uppermost mantle that is a consequence of mountain perpendicular compressive stresses the preferred orientation of the fast axes of some anisotropic minerals such as olivine is known to be in the creep direction a consequence of the anisotropy of strength and viscosity of orientated minerals in order to explain the anisotropy of the mantle beneath young orogens we extend the concept of crustal escape or extrusion tectonics to the uppermost mantle we present rheological model calculations to support this hypothesis mountain perpendicular horizontal stress determined in the upper crust and mountain parallel seismic anisotropy in the uppermost mantle require a zone of ductile decoupling in the middle or lower crust of young mountain belts examples for stress and mountain parallel pn anisotropy are given for tibet the alpine chains and young mountain ranges in the americas finally we suggest a simple model for initiating mountain parallel creep

laboratory and seismological observations of lower mantle isotropy we have carried out seismological and mineral physics investigations to identify the extent and origin of elastic anisotropy in the lower mantle based on observations of shear wave splitting we conclude that the lower mantle is effectively isotropic this result is surprising since the lower mantle is composed of elastically anisotropic minerals silicate perovskite mgo sio 2 that have been extensively strained by geologic processes to reconcile the seismic observations we have measured texture development during deformation and high temperature recrystallization of mantle silicates we show that these experiments provide a direct explanation for the shear wave splitting measurements the lower mantle appears isotropic because silicate perovskite maintains an isotropic texture during deformation and recrystallization

seismic anisotropy of the pacific ocean inferred from long period surface waves dispersion abstract the dispersion of surface rayleigh and love waves in the period range 40 300 s along a large number of paths allows the estimation of both the azimuthal anisotropy and the shear wave polarization anisotropy the regional dispersion is determined taking into account simultaneously its dependence with age and an azimuthal factor the pacific ocean has been divided into 5 regions for rayleigh waves and into 3 regions for love waves this partition discriminates the regions of extreme age which show a fast variation of dispersion with age from the regions of intermediate age where the variation is weak a variation of 2 of rayleigh wave group velocity with the azimuth of the path measured with respect to the direction of spreading is displayed up to very long period on the contrary the azimuthal anisotropy for love waves is difficult to resolve for rayleigh waves the present day direction of plate motion seems to agree best with the direction of maximum velocity on the other hand the isotropic inversion of the regional dispersion curves indicates except for young regions a discrepancy between rayleigh wave and love wave models with this hypothesis sh velocities are higher than sv velocities for the regions older than 23 ma down to a depth of 300 km which is indicative of the presence of polarization anisotropy the latter very weak for the young part of the ocean increases with age and reaches 7 for the oldest region

global anisotropic phase velocity maps for higher mode love and rayleigh waves summary it is well established that the earth s uppermost mantle is anisotropic but there are no clear observations of anisotropy in the deeper parts of the mantle surface waves are well suited to observe anisotropy since they carry information about both radial and azimuthal anisotropy fundamental mode surface waves for commonly used periods up to 200 s are sensitive to structure in the first few hundred kilometres and therefore do not provide information on anisotropy below higher mode surface waves have sensitivities that extend to and beyond the transition zone and should thus give insight about azimuthal anisotropy at greater depths we have measured higher mode love and rayleigh phase velocities using a model space search approach which provides us with consistent relative uncertainties from measurement to measurement and from mode to mode from these phase velocity measurements we constructed global anisotropic phase velocity maps prior to inversion we determine the optimum relative weighting for anisotropy we present global azimuthal phase velocity maps for higher mode rayleigh waves up to the sixth higher mode and love waves up to the fifth higher mode with corresponding average model uncertainties the anisotropy we derive is robust within the uncertainties for all modes given the ray theoretical sensitivity kernels of rayleigh and love wave modes the source of anisotropy is complex but mainly located in the asthenosphere and deeper our models show a good correspondence with other studies for the fundamental mode but we have been able to achieve higher resolution

the depth distribution of mantle anisotropy beneath the tonga subduction zone abstract shear wave splitting recorded by a temporary deployment of broadband seismometers located above the tonga subduction zone yields strong constraints on the depth distribution of anisotropy in the mantle beneath the tonga back arc region splitting parameters obtained for local s and teleseismic sks phases were modeled using a ray based method that predicts shear wave splitting on individual phase paths and determines the best fitting values of anisotropic strength orientation and depth extent splitting in teleseismic sks phases is identical to that in s phases from local earthquakes that occur at the base of the transition zone demonstrating that there is no significant splitting due to anisotropy in the lower mantle splitting times for local s phases do not vary significantly with depth and for models in which the strength of anisotropy is laterally uniform they rule out anisotropy in the transition zone finally the observed splitting indicates that anisotropy of more than 0 8 exists in the upper mantle with a fast symmetry axis at 60 w roughly parallel to the absolute plate motion vector of the subducting pacific lithosphere

constraints on north atlantic upper mantle anisotropy from s and ss phases several analyses of ssh sh differential travel times disagree on whether these data contain evidence for azimuthal anisotropy in the north atlantic upper mantle to explore possible explanations for these discrepant results we calculated the azimuthal variation in ssh sh travel times predicted by an anisotropic upper mantle model theoretical travel times were obtained by cross correlating s and ss phases on transverse component reflectivity synthetic seismograms we determined that predicted azimuthal travel time patterns are sensitive to the relative amplitudes of radial and transverse particle motions incident on the anisotropic medium if a large variation in the relative incident amplitudes exists in a dataset the overlap of different azimuthal travel time patterns may obscure anisotropic signatures we analyzed shear wave splitting in a dataset of 12 ss and 10 s phase observations to constrain the presence of anisotropy in the north atlantic upper mantle a comparison of ss splitting parameters with s and sks particle motions suggests that shear wave splitting apparent in the ss phases is at least partially due to anisotropy in the vicinity of their bounce points although contributions from anisotropy in the source and station regions cannot be ruled out observed ss splitting is consistent with an olivine model in which horizontal a axes are aligned at an average of n24 w in the central north atlantic 15 40 n and at an average of n23 e in the northern north atlantic 40 50 n

global anisotropy in the upper mantle inferred from the regionalization of phase velocities a data set of 2600 paths for rayleigh waves and 2170 paths for love waves was selected in order to resolve azimuthal anisotropy of surface waves in the period range 70 250 s the epicenters were chosen so that the paths provide uniform spatial and azimuthal coverage of earth a synthetic seismogram is computed for each path and compared to the observed seismogram in order to infer the phase velocity along each path to process such a large data set a new tomographic technique is designed a good correlation is found with surface tectonics in the whole period range for love waves and at shorter periods for rayleigh waves the directions of maximum velocities are also well correlated with plate tectonics motions for the 2 azimuthal term of rayleigh waves when the power spectrum of heterogeneities is computed it is found that for love waves the even orders are usually larger than the odd orders for rayleigh waves the degrees 2 and 6 dominate only at long periods suggesting a deep origin for these components

lateral heterogeneity and azimuthal anisotropy of the upper mantle love and rayleigh waves 100 250 s the lateral heterogeneity and apparent anisotropy of the upper mantle are studied by measuring rayleigh and love wave phase velocities in the period range 100 250 s spherical harmonic descriptions of the lateral heterogeneity are obtained for order and degree up to 1 m 10 slow regions are evident at the east pacific rise northeast africa tibet tasman sea southwestern north america and triple junctions in the northern atlantic and indian oceans fast regions occur in australia western pacific and the southern atlantic these features are also found by a completely different analysis based on the backus gilbert method the backus gilbert method also shows that the obtained phase velocities are averaged values within an area of about 2000 km radius and the errors are about 1 of the phase velocity in the zeroth order spherically symmetric earth inversion for azimuthal dependence shows that for low angular order the fast phase velocity directions seem to correlate well with the plate motion vectors however resolution analysis by the backus gilbert method shows that the current data do not have enough resolution for everywhere on the globe only a few regions currently have adequate azimuthal coverage thus confirmation of the above correlation requires a more complete data set

anisotropy of the african continent inferred from surface waves abstract a no block regionalization is performed over the african continent to retrieve simultaneously local surface wave velocities and the 2 azimuthal variation of rayleigh waves the geographical distribution of phase and group velocities shows general agreement with global models but with enhanced lateral resolution a correlation with surface tectonics is found even at long periods low velocities occur in the tectonically active regions and high velocities in the shield areas the azimuthal anisotropy analysis results in a general north south fastest direction for rayleigh waves which could be correlated with the african plate s absolute motion the amplitude is of the order of 2 an isotropic inversion with data corrected for effects of surficial layers is performed very slight differences appear between the corrected and uncorrected results lateral heterogeneities of rayleigh and love waves are then inverted simultaneously with the 2 azimuthal terms of rayleigh waves by vectorial tomography to retrieve the 3d structures of s wave velocity and anisotropy throughout the african continent good correlation between low velocities and small anisotropies is found and their location at depth suggests an origin for african hotspots

analysis of seismic sv waves in the core s penumbra vertically polarized s waves sv attenuate rapidly as they diffract around the earth s core whereas horizontally polarized s waves sh propagate to large distances in the core shadow zone the amplitude decay of diffracted sv signals is so strong that few studies have been made of these phases despite their acute sensitivity to velocity structure just above the core mantle boundary we analyze sv signals in the penumbra of the core s shadow finding systematic waveform complexities indicative of local stratification and apparent anisotropy at the base of the mantle with implications for the dynamic processes in this internal boundary layer

mineralogy and composition of the upper mantle seismic velocities are calculated for two petrological models of the upper mantle an olivine rich assemblage pyrolite and a garnet clinopyroxene rich olivine bearing assemblage piclogite these are compared with recent seismic profiles for various tectonic provinces the shield data is most consistent with a cold olivine and orthopyroxene rich lid the seismic lithosphere extending to 150 km followed by a high temperature gradient and or a change in mineralogy that serves to decrease the velocity from 200 to 400 km the velocities follow a 1400 c adiabat the rise tectonic mantle is much slower presumably hotter and is likely to be above the solidus to depths of at least 300 km the high v p v s ratio of the lower oceanic lithosphere in the western pacific is most consistent with eclogite

petrological constraints on seismic anisotropy the purpose of this paper is to investigate the correlations between anisotropic parameters for different orientations and mineralogies for realistic mineralogical and petrological models of the upper mantle such correlations make it possible to reduce the number of independent parameters required in seismic modelling of the upper mantle the variation with depth and with degree of random orientation is also investigated it is shown that the anisotropic combinations involved in an equivalent transversely isotropic medium azumithally averaged velocities with a vertical symmetry axis are strongly correlated but the combinations involved in azimuthal anisotropy are more weakly correlated comparison with available seismic data does not yet allow discrimination based on anisotropy between different petrological models it is shown that even if a proportion of random orientation is introduced it should be possible to achieve petrological discrimination from seismic data for some orientation configurations if the inverse problem is well posed i e takes account of the correlations between anisotropic parameters these correlations are incorporated in an inversion process to show that it is possible to explain the dataset used to derive the prem by alternative anisotropic models for the uppermost 220 km

pyroxene orientation within the upper mantle petrofabric analyses of orthopyroxene and olivine for a suite of rocks from the bay of islands ophiolite newfoundland show systematic relationships between the minerals9 fabrics the plane of the mohorovicic discontinuity determined from the contact of the gabbro and ultramafic sections and the spreading direction inferred from the attitude of the sheeted dikes the orthopyroxene c crystallographic axes maxima which are aligned subparallel to the spreading direction coincide with olivine a axes maxima the orthopyroxene b axes and olivine c axes also lie approximately within the plane of the mohorovicic discontinuity but parallel the ridge whereas the orthopyroxene a axes maxima and olivine b axes maxima are approximately perpendicular to the mohorovicic discontinuity using olivine and orthopyroxene elastic constants and their pressure derivatives velocities have been calculated from composite petrofabric data and contoured for varying proportions of olivine and pyroxene the results of these calculations show that within a vertical direction compressional wave velocities change only slightly with increasing pyroxene content whereas parallel to the mohorovicic discontinuity velocities decrease with increasing pyroxene content anisotropy decreases from 4 for a rock with 100 olivine to 2 6 for a rock with 50 olivine and 50 orthopyroxene also within the plane of mohorovicic discontinuity the maximum and minumum velocities occur at approximately the same azimuth for differing percentages of orthopyroxene a comparison of the calculated anisotropy and the observed seismic anisotropy off the coast of california shows a striking agreement from this comparison it is estimated that the percentage of orthopyroxene within the upper mantle of this region varies between 0 and 30

mid wisconsinan eolian deposits of the kittigazuit formation tuktoyaktuk coastlands northwest territories canada the kittigazuit formation is a late quaternary sand unit commonly observed throughout the tuktoyaktuk coastlands and in subbottom sediments of the southern beaufort shelf stratigraphic and sedimentology data including sedimentary structures grain surface characteristics and heavy and light mineralogy assist in characterizing the deposit and indicate that it is eolian in origin plant and arthropod macrofossils suggest that although the summer climate during deposition was as warm or slightly warmer than today conditions were likely more arid permafrost is interpreted as being widespread during deposition accelerator mass spectrometry radiocarbon dates indicate that the kittigazuit formation was deposited between 37 and 33 ka bp the unit is therefore interpreted as a mid wisconsinan eolian dune deposit formed by reworking of underlying alluvial sediments of the kidluit formation under nonglacial conditions glaciogenic sediments overlying the kittigazuit formation indicate that glacial ice cover

aspherical earth structure from fundamental spheroidal mode data free oscillation data reveal heterogeneity in the earth s mantle whose geographical pattern is dominated by spherical harmonics of angular degree two and correlates well with the hydrostatically referenced geoid the heterogeneity can be modelled as localized in the transition zone 420 670 km depth and may be related to a large scale component of mantle convection

a quantitative evaluation of the contribution of crustal rocks to the shear wave splitting of teleseismic sks waves abstract the seismic elastic properties of the crust have been modelled using hypothetical polycrystals with typical crustal compositions and commonly measured petrofabrics this modelling allows us to better evaluate and quantify the crustal contribution to splitting of teleseismic sks waves mafic rocks such as anorthosite or pyroxene bearing gabbros have complex s waves properties i e with small shear wave anisotropies

the three dimensional seismological model a priori constrained confrontation with seismic data we compare the predictions of an a priori model of the upper mantle with seismic observations of surface waves and eigenmodes the 3 dimensional seismological model a priori constrained 3smac has been developed by natal and ricard 1996 it is based on the interpretation by geodynamicists of the near surface layers of the earth on distributions of temperature pressure and composition as a function of depth and then on estimates of seismic parameters density velocities attenuations from solid state laboratory measurements as a function of temperature and pressure the 3smac predictions are confronted with observations consisting of phase velocities for love and rayleigh waves in the period range of 70 250 s montagner and taniraoto 1990 we first show that tomographic inversions applied to 3smac synthetics induce a strong smoothing of the heterogeneities this ca sts doubt on the meaning of the spectra of mantle heterogeneities revealed by tomography we then show that most of the love and rayleigh fundamental mode observations for periods less than 200 s are satisfactorily predicted by 3smac the major differences come from the seismic velocities under the red sea and southeast china which are much slower than what is estimated from 3smac as well as those under greenland which are not as fast as the other cratonic areas because the lithosphere is thinner than 100 km under oceans and thinner than 300 km under continents in 3smac we suggest that the existence of deeper lithospheric anomalies as proposed in many tomographic models is mostly due to a spurious effect of the inversion rather than implied by surface wave data half of the variance of the degree 2 anomaly mapped by low degree eigenmode observations can be explained by lithospheric velocity structures the other half is highly correlated with the distribution of deep slabs but its amplitude is a factor of 3 or 4 larger than that predicted by 3smac the lithospheric anomalies present a degree 6 pattern well correlated with the distribution of hotspots even when the thermal anomalies that could be associated with plumes are not included in 3smac our results emphasize the importance of giving very close attention to surface corrections in tomographic models

anisotropy in kuril kamchatka subduction zone structure to constrain the orientation and extent of seismic velocity anisotropy in the kuril kamchatka subduction zone we have developed a technique that isolates anisotropy at shallow depths by differencing splitting observed in teleseismic shear wave phases olivine a axis orientations inferred from 29 ss s phase pairs and two local s phases vary between the major tectonic components of the trench back arc system although the ss paths do not uniquely constrain a single model of anisotropy a axes are consistent with finite strain patterns predicted by plate convergence and back arc spreading and shearing

uppermantle anisotropy and the oceanic lithosphere rayleigh wave and love wave dispersion data for oceanic paths cannot be satisfied by an isotropic uppermantle and it is incorrect to invert these two datasets independently for separate isotropic structures available phase and group velocity data are inverted for oceanic structure as a function of age taking into account anelastic dispersion sphericity and anisotropy the resulting models are quite different from previous results which ignore the above effects and the extra parameters involved in proper anisotropic inversion the models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age and an anisotropic low velocity zone with an age dependent anisotropy the lid or seismic lithosphere reaches a maximum thickness of 50 km which is about one half the thickness obtained in previous isotropic or pseudo isotropic inversions the seismic lithosphere therefore may be comparable in thickness to the elastic or flexural lithosphere raising the possibility that both are controlled by effects other than temperature stress and time such as mineralogy crystal orientation or partial melting the velocities in the low velocity zone are higher than inferred by previous surface wave studies

the nature of the lehmann discontinuity from its seismological clapeyron slopes abstract while the existence of the lehmann discontinuity at 220 km in earth s mantle has been known for over 40 years from seismology it is still debated what causes this discontinuity here we report measurements of seismological clapeyron slopes for the lehmann discontinuity which are key to explaining the seismic discontinuity as either a manifestation of phase transitions or of other mineral physical processes the clapeyron slopes are measured by correlating discontinuity depths with local velocity perturbations from a tomographic model assuming that the velocity perturbations are solely due to temperature variations we find that in most regions the lehmann discontinuity is characterised by a regionally varying negative seismological clapeyron slope reflections from greater depths 250 350 km which we associate with the x discontinuity appear in clusters without a clear clapeyron slope in seeking mineralogical explanations of the lehmann and x discontinuities we can reject hypotheses in disagreement with these observations known phase transitions in the upper mantle above 400 km depth all have positive clapeyron slopes in the case of the lehmann discontinuity the only remaining hypothesis for a negative clapeyron slope is that it represents the transition in deformation mechanism from dislocation to diffusion creep the only exception is the middle east region where a positive claperyon slope is observed which could be explained by a phase transition from coesite to stishovite in the case of the x discontinuity we are not aware of mineral physical mechanisms possessing the required behaviour

anisotropy and shear velocity heterogeneities in the upper mantle long period surface waves are used to map lateral heterogeneities of velocity and anisotropy in the upper mantle the dispersion curves are expanded in spherical harmonics up to degree 6 and inverted to find the depth structure the data are corrected for the effect of surface layers and both love and rayleigh waves are used shear wave velocity and shear polarization anisotropy can be resolved down to a depth of about 450 km the shear wave velocity distribution to 200 km depth correlates with surface tectonics except in a few anomalous regions below that depth the correlation vanishes cold subducted material shows up weakly at 350 km as fast s wave anomalies in the transition region a large scale pattern appears with fast mantle in the south atlantic s anisotropy at 200 km can resolve uprising or downwelling currents under some ridges and subduction zones the pacific shows a nw se fabric

lehmann discontinuity as the base of an anisotropic layer beneath continents long period surface wave r 1 g 1 body wave s ss sss and scs reverberation data have been inverted to obtain anisotropic structures along seismic corridors that sample australia and the western pacific these models support the proposal that the lehmann discontinuity beneath stable continents represents a transition from an anisotropic lithosphere to a more isotropic material in the lower part of the continental tectosphere

3 dimensional structure of the indian ocean inferred from long period surface waves to improve the lateral resolution of the first global 3 dimensional models of seismic wave velocities regional studies have to be undertaken the dispersion of rayleigh waves along 86 paths across the indian ocean and surrounding regions is investigated in the period range 40 300 s the regionalization of group velocity according to the age of the sea floor shows an increase of velocity with age up to 150 s only similar to the results in the pacific ocean but here this relationship vanishes more quickly at long period therefore the correlation of the deep structure with surface tectonics seems to be shallower in the indian ocean than in the pacific ocean a tomographic method is applied to compute the geographical distributions of group velocity and azimuthal anisotropy and then the 3 d structure of s wave velocity horizontal wavelengths of 2000 km for velocity and 3000 km for azimuthal anisotropy distribution can be resolved except for the central part of the south east indian ridge which displays high velocities at all depths the inversion corroborates a good correlation between lithospheric structure down to 120 km and surface tectonics low velocities along the central and southeast indian ridges velocity increasing with the age of the sea floor high velocities under african indian and australian shields at greater depths the low velocity zones under the gulf of aden and the western part of the southeast indian ridges hold but the low velocity anomaly of the central indian ridge is offset eastward the low velocity anomalies suggest uprising material and complex plate boundary

joint interpretation of upper mantle anisotropy based on teleseismic p travel time delays and inversion of shear wave splitting parameters abstract we inverted shear wave splitting parameters and simultaneously analysed delays of teleseismic longitudinal waves to obtain a self consistent three dimensional 3 d image of the anisotropic upper mantle beneath the continents efficiency of the simultaneous 3 d analyses of p residual spheres and shear wave polarizations is demonstrated on data from two regions southern sweden and the western usa the anisotropic inferences of the subcrustal lithosphere are to a first approximation represented by homogeneous hexagonal or orthorhombic media with plunging symmetry axes

upper mantle anisotropy a preliminary model seismic anisotropy in the upper mantle can be explained by crystallographic mineral alignment achieved through dislocation motion the physical mechanism of mineral alignment requires upper mantle shear flow which reorients and aligns minerals by dislocation glide and climb governed by the dominant glide system of each mineral the dominant glide systems are assumed to be 100 0l0 for olivine and 001 100 for the pyroxenes these yield a predominantly orthorhombic fabric with the olivine 100 and the pyroxene 001 axes aligned in the upper mantle flow direction and the olivine 010 and the pyroxene 100 axes aligned normal to the upper mantle flow plane these glide systems have a threshold temperature of enhanced mobility of 1100 1200 k which yields a solid state thermally defined lithosphere asthenosphere boundary in an olivine pyroxene mantle consistent with recent seismic determinations of the thickness of the lithosphere mantle anisotropy due to mineral alignment is then actively maintained below this boundary asthenosphere and mesosphere and is a fossil state above this boundary lithosphere we use the dominant glide systems to establish the crystallographic orientation of olivine and pyroxene in calculating the maximum seismic anisotropy of two petrologic models pyrolite and piclogite for the upper mantle by an extrapolation of single crystal anisotropic mineral elastic properties to a depth of 400 km the real earth seismic anisotropy will be bound by the limits of maximum anisotropy from perfect mineral alignment and minimum anisotropy isotropy from random mineral alignment the seismic anisotropy of the upper 220 km is best represented by the pyrolite model which reduces to quasi hexagonal symmetry with the unique axis in the direction of mantle flow the lehmann discontinuity is conjectured to be due to a change in composition from pyrolite to piclogite and therefore may represent a change in anisotropy the piclogite model has orthorhombic symmetry distinctly different from that of the pyrolite the piclogite anisotropy model can appear quasi isotropic when measured by transverse isotropy parameterization if the mantle flow is mainly horizontal with a horizontal shear plane that is the shear flow gradient duh dr dominates

lowermost mantle anisotropy beneath the north pacific from differential s scs splitting summary seismic anisotropy is an important tool for studying the nature origin and dynamics of the lowermost mantle d we introduce differential s scs splitting as a tool for removing the effect of near source and near receiver anisotropy to estimate splitting accrued in the d region this is applicable to events recorded at epicentral distances between 60 and 85 near source anisotropy has often been ignored in previous studies of lowermost mantle anisotropy we apply differential s scs splitting to records from canadian national seismic network stations of western pacific earthquakes these sample the lowermost mantle beneath the north pacific the residual splitting in scs which we attribute to d shows lag times between 1 0 and 3 9 s given the near horizontal ray path of scs in d we interpret the recovered fast directions as the orientation of the fast shear wave in the plane defined by the vertical and transverse directions and observe a clearly non vti transverse isotropy with a vertical axis of symmetry style of anisotropy the largest population of results shows an approximately southeasterly dipping symmetry axis which we speculate might be explained by descending palaeoslab material being swept horizontally across the core mantle boundary towards an upwelling region beneath the central pacific non vti symmetry and the many possible contributions to d anisotropy from lower mantle minerals melt and subducted materials suggest that our understanding of the lowermost mantle could be greatly improved by trying to resolve a more general style of anisotropy

seismic structure of the upper mantle in a central pacific corridor the seismic structure of the tonga hawaii corridor has been investigated by combining two data sets revenaugh and jordan s reflectivity profile from scs reverberations which provides travel times to and impedance contrasts across the major mantle discontinuities and 1500 new observations of frequency dependent phase delays for the three component s ss and sss body waves and the r1 and g1 surface waves which constrain the velocity structure within this layered framework the shear waves turning in the upper mantle showed significant splitting of the sh and sv components indicative of shallow polarization anisotropy the data set was inverted in conjunction with attenuation and mineralogical constraints to obtain a complete spherically symmetric radially anisotropic structure the final model pa5 is characterized by a high velocity anisotropic lid bounded at 68 km depth by a large negative g discontinuity a low velocity anisotropic layer below g extending to a small l discontinuity at 166 km an isotropic steep gradient region between 166 km and 415 km and transition zone discontinuities at 415 507 and 651 km the depth of the radial anisotropy in pa5 is shallower than in most previous studies based on surface waves and higher modes the average value of radial shear anisotropy in the lid 3 7 is consistent with the magnitude expected from the spreading controlled models of olivine orientation while anisotropy in the low velocity zone which is required by our data set could be induced either by paleostrains that took place near the ridge crest or by shearing in the asthenosphere as a result of present day plate motions on the basis of recent work by hirth and kohlstedt we suggest that the g discontinuity is caused by a rapid increase in the water content of mantle minerals with depth marking the fossilized lower boundary of the melt separation zone active during crust formation the high gradient zone between 200 and 400 km is a characteristic feature of convecting oceanic upper mantle and is probably controlled by a steady decrease in the homologous temperature over this depth interval the average shear velocity gradient in the transition zone is lower than in most previous seismic models in better agreement with the predictions for a pyrolitic composition

elastic thickness structure of the andaman subduction zone implications for convergence of the ninetyeast ridge abstract we use the bouguer coherence morlet isostatic response function technique to compute the spatial variation of effective elastic thickness te of the andaman subduction zone the recovered te map resolves regional scale features that correlate well with known surface structures of the subducting indian plate and the overriding burma plate the major structure on the india plate the ninetyeast ridge ner exhibits a weak mechanical strength which is consistent with the expected signature of an oceanic ridge of hotspot origin however a markedly low strength 0

inversion of multimode surface wave data evidence for sub lithospheric anisotropy summary multiparameter inversions of multimode dispersion data are performed for two large regions the pacific ocean and north america anisotropy is taken into account by considering transversely isotropic structures with a vertical axis of symmetry two fundamental questions are studied in detail 1 how to make the inverted models consistent when using different sets of parameters 2 what is the significance of transversely isotropic inversion for the actual earth s structure it is proved that full consistency of the inverted models can be achieved by properly taking into account some a priori informations on the model and it is shown that the use of transversely isotropic models with vertical axis of symmetry does not cause severe limitations when interpreting the data the models we have obtained are discussed in the light of these investigations considering an olivine rich upper mantle we make a tentative interpretation of these models in terms of preferred orientation of the a axis of the crystals in one fixed horizontal direction

shear wave splitting in the mantle of the pacific summary the bulk of previously known data on seismic anisotropy in the suboceanic mantle was obtained from observations of azimuthal variations of the surface wave and p velocities our study of azimuthal anisotropy in the mantle of the pacific is based on observations of shear wave splitting the idea of our method of measurement vinnik farra 1992 is to analyse evolution of the waveform of the long period multiple reflected scs and sscs phases as a function of the order of reflection the estimates of anisotropy correspond to circular regions with a radius of several hundred kilometres around the surface bounce points the technique was applied to the digital recordings of deep and intermediate events from the circum pacific seismic belt at seismograph stations in and around the pacific the estimates of polarization direction of the fast split wave and time lag of the slow wave are obtained for a few regions of the pacific basin marginal basins and the andes in the back arc extensional basins whose spreading has continued to the present time or ended later than about 15ma the inferred direction of extension in the mantle is close to that imprinted in the crust lau basin japan sea in the back arc basins whose active development ended 15 20 ma or earlier south fiji and parece vela basins the inferred directions of extension in the mantle are in strong disagreement with the crustal data in the central andes the region of plate collision and low angle subduction the inferred direction of flow in the mantle is perpendicular to the plate boundary contrary to the regularity previously found in several other collisional belts the fast direction of anisotropy in the mantle beneath the southwestern part of the 100ma old pacific plate is very close to the absolute plate motion direction the magnitude of anisotropy corresponds to the lowest values found in the continental mantle and is much lower than in most of the marginal basins and the andes the weakness of this presumably asthenospheric anisotropy in comparison with asthenospheric anisotropy under some old continental plates is probably related to a difference in upper mantle viscosities

finite strain calculations of continental deformation 2 comparison with the india asia collision zone numerical experiments on a thin viscous sheet model for deformation of continental lithosphere subjected to an indenting boundary condition yield distributions of crustal thickness of stress and strain rate and of latitudinal displacements that may be compared with observations in the india asia collision zone a simple indenting boundary condition applied to initially laterally homogeneous sheets obeying a power law rheology produces results that are in broad agreement with the observations provided that the power law exponent is three or greater and the sheet can support vertically integrated stress differences of 2 1013 5 1012 n m 1 in regions in front of the indenter under these conditions the calculated deformation shows accommodation of convergence primarily by crustal thickening to produce a plateau in front of the indenter palaeomagnetic data from india and tibet and the observed distribution of topography suggest that much of the post eocene convergence of india with asia has been taken up by deformation within asia that involved crustal thickening the principal difference between calculation and observation is the absence from the calculated strain rate fields of east west extension of the plateau in front of the indenting boundary the calculations show that once such a plateau is formed the buoyancy force associated with the crustal thickness contrast inhibits further thickening and the plateau strains at less than half the rate of its immediate surroundings seismically determined regional strain rates exhibit a similar distribution with the tibetan plateau straining at about one quarter the rate of the tien shan and ningxia gansu regions calculated principal compressive stress orientations and regional strain rates agree with the seismically determined quantities in the mongolia baikal tien shan tibet and ningxia gansu regions of asia to within the uncertainty of the latter the vertically integrated stresses that are calculated for the viscous sheet are comparable with those that can be supported by a theologically stratified continental lithosphere obeying laboratory determined flow laws we suggest that the thin viscous sheet model described in this paper and its companion gives a simple and physically plausible description of the observed deformation in central asia in this description the predominant mechanism of accommodation of continental convergence is diffuse crustal thickening with shear on vertical planes playing a subsidiary role once large crustal thickness contrasts have been established

a search for source side mantle anisotropy shear wave splitting due to source side mantle anisotropy is estimated from broadband teleseismic s waves using a suite of 19 subduction zone events more than 50km depth and recorded at several rstn stations we assume that there are at most two anisotropic regions along the ray path one on the source side and one on the receiver side and that receiver side splitting is known from analysis of sks it is found that splitting is not detectable below the slab in most cases there are however two specific regions one beneath peru and one beneath kamchatka that indicate significant source side anisotropy based on two nearby intermediate focus events in each region splitting delay times are found to be 1 0 1 5s with the fast polarization direction perpendicular to the trench assuming 4 intrinsic anisotropy this would correspond to path lengths of 100 170 km below the slab both regions are characterized by abrupt changes in along strike slab properties a change in slab dip in the case of peru and a cessation of deep focus seismicity in the case of kamchatka we speculate that the anisotropy is caused by an associated local disturbance in asthenospheric flow that the shear wave splitting parameters of s and sks are consistent with each other in most cases suggests that large scale anisotropy in the lower mantle is very small probably no more than a fraction of a percent

a simple model of the 1982 83 californian el nino steric height steric height anomaly sea level and sea level anomaly show that the 1940 41 and the 1982 83 californian el nino events were resonant responses of the california current to the expansion and intensification of the aleutian low and decrease in strength of the north pacific high in both cases the subsurface anomalies were dynamically produced by onshore transport of subarctic water and a downward tilt of the inshore thermocline sea surface temperature anomalies however resulted from thermodynamic as well as dynamic processes most previous californian el nino or anti el nino events probably were resonant or anti resonant responses of the seasonal cycle of the california current to anomalous large scale atmospheric forcing

anisotropy of the asthenosphere the higher mode data of the pacific revisited surface wave data previously used by the authors to constrain upper mantle anisotropy in the pacific ocean are re examined after an error was found in the calculations of the love wave phase velocities sub lithospheric anisotropy is still required by the corrected data but it extends to shallower depth than we originally proposed using different parameterizations of the upper mantle we present strong evidence that the bottom of the anisotropic layer has been reached around 250 km in the pacific beneath north america no revision of our data is required and anisotropy could extend to around 400 km depth our results are in agreement with an upper mantle anisotropy due to preferred orientation of olivine crystals in an overall horizontal shear strain parallel to the present day plate motion we suggest that this preferred orientation of the crystals mainly occurs within the upper mantle low velocity layer

rayleigh wave phase velocities in the pacific with implications for azimuthal anisotropy and lateral heterogeneities summary the lateral distribution of fundamental mode rayleigh wave phase velocities in the pacific has been calculated in order to determine the variation in velocity as a function of the age of the oceanic plate the importance of azimuthal anisotropy and the presence of secondary lateral heterogeneities the data set used in this study consists of phase velocity measurements in the period range 20 125 s for 178 paths traversing the pacific both the pure path and spherical harmonic inversion techniques are used in this investigation with an emphasis placed on determining the effectiveness and resolving power of these methods in calculating the velocity distribution in an oceanic regime the first order velocity age relation which dominantly characterizes surface wave dispersion in an oceanic environment can be adequately modelled by either the pure path or spherical harmonic techniques azimuthal anisotropy is the dominant second order effect and is best developed in regions less than 80 myr in age in this age zone the faster direction of wave propagation is either in the direction of fossil seafloor spreading or in the direction of present day absolute plate motion in the western pacific there is a correlation with paleo relative motion fossil seafloor spreading for periods less than 50 s lateral velocity heterogeneities which are not related to the age of the oceanic plate are found by pure path inversions and are correlated with regions associated with shallow residual depth anomalies this effect slow velocities with respect to the velocity age model monotonically decreases as a function of period and is indicative of a shallow origin age independent anomalies were also modelled by the application of the pure path and spherical harmonic methods in a sequential inversion

observational evidence for diffracted sv in the shadow of the earth s core on the basis of observational evidence it has been assumed for a long time in the seismological literature that the diffracted s wave in the deep shadow of the earth s core is polarized as sh in good agreement with models that do not contain a well pronounced low velocity zone at the base of the mantle while this seems often to be the case broad band horizontal records from the geoscope network show that for some wavepaths the diffracted sv is strong enough to be clearly observed at distances between 107 and 117 moreover in some cases the diffracted sv is shifted with respect to sh by approximately a quarter period such observations should put additional constraints on the structure near the core mantle boundary among other possibilities lateral variations of the gradient at the base of the mantle allowing for regions with negative gradient and possibly azimuthal anisotropy in d should be considered in order to explain these data

mountain building strike parallel motion and mantle anisotropy abstract geophysical and geological evidence indicates that tectonic movements parallel to the strike of orogenic belts and suture zones are the most important mode of displacement during continental collision and oblique subduction several hundred to a few thousand kilometers of displacement may have occurred this way resulting in a pervasive ductile deformation of the lithosphere consequently both crustal and mantle fabrics may reflect these movements and their study using both geological and geophysical methods could supply important clues to lithosphere behavior during mountain building on the basis of kinematic analysis of ductilely deformed crustal rocks in zones of movement a variety of tectonic models have been proposed which differ in 1 1 the timing of strike parallel motions early before or during continental accretion or late after accretion was completed 2 2 the geometric pattern of strike parallel displacement zones low angle or high angle and 3 3 the location of orogen parallel movement zones in outboard terranes on the subducting plate along the suture zone or within the continental hangingwall large scale displacements parallel to mountain belts also suggest that they occurred along lithospheric faults rooted in the upper mantle although few studies have been carried out from this view point thermal structural and petrological data support the conclusion that lithospheric faults involving the upper mantle may be discriminated from crustal faults rooted in some crustal decoupling zone where available data on mantle flow beneath orogenic belts are in agreement with mountain parallel displacements inferred from crustal studies this has an important bearing on the seismic anisotropy pattern which can be predicted in the lithospheric mantle below mountain belts the fast velocity direction which is parallel to the flow direction should therefore be recorded parallel to the strike of the belt

p wave anisotropic tomography of the alps the first tomographic images of p wave azimuthal and radial anisotropies in the crust and upper mantle beneath the alps are determined by joint inversions of arrival time data of local earthquakes and teleseismic events our results show the south dipping european plate with a high velocity high v anomaly beneath the western central alps and the north dipping adriatic plate with a high v anomaly beneath the eastern alps indicating that the subduction polarity changes along the strike of the alps the p wave azimuthal anisotropy is characterized by mountain chain parallel fast velocity directions fvds in the western central alps and ne sw fvds in the eastern alps which may be caused by mantle flow induced by the slab subductions our results reveal a negative radial anisotropy i e vph vpv in the low velocity mantle wedge which may reflect the subvertical plate subduction and its induced mantle flow the results of anisotropic tomography provide important new information on the complex mantle structure and dynamics of the alps and adjacent regions

shear wave polarization anisotropy in the upper mantle beneath honshu japan shear wave polarization anisotropy in the wedge portion of the upper mantle between a subducting plate and the earth s surface is investigated using three component seismograms of intermediate depth and deep earthquakes recorded at 14 local stations in honshu japan eighty nine high quality seismograms were selected from a period of 3 years the data used in this study are restricted such that incidence angles are smaller than the critical angle of 30 to the earth s surface in order to avoid phase shifts in the shear wave train to find directions of the maximum and minimum velocities in split shear waves where shear waves are resolved into two phases with the maximum time separation each set of the two horizontal component seismograms is rotated in the horizontal plane the split shear waves thus obtained are again recombined after the correction of anisotropy and the anisotropy corrected particle motion is compared with the focal mechanism for a cross check of the observed anisotropy directions of the maximum axes are plotted on azimuthincidence angle stereograms at each station the stereograms and the cross sections of seismic ray paths show that 1 the anisotropic material is distributed at intermediate location between earthquake sources and receiving stations and 2 the anisotropic region is separated into two parts one in the north of the present study area with the polarization of the maximum velocity shear wave trending 0 to 30 from the north north anisotropy and the other in the south with it trending 90 to 120 south anisotropy the maximum time delays between the two shear waves along a vertical seismic ray is about l s for both the anisotropic regions the horizontal extent of the anisotropic area in the north is 50 km at depths of 50 to 150 km perhaps prevalent in west honshu however the depth of anisotropic material is not well constrained because of insufficient coverage of seismic ray paths and angles if we assume the vertical extent of the anisotropic material to be 100 km the maximum velocity contrast would be 4 if we adopt a crack alignment model for the observed anisotropy cracks are inferred to be distributed vertically trending 0 to 30 in the north and 90 to 120 in the south if we assume an olivine alignment model olivine crystals are inferred to be aligned in the north with the a axis of the slip direction trending 0 to 30 and the slip plane normal to the b axis being vertical and in the south anisotropy with the axis trending 90 to 120 and the slip plane being horizontal

anisotropy beneath 9 stations of the geoscope broadband network as deduced from shear wave splitting polarization of sks scs and s waves has been analyzed on digital 3 components seismograms recorded at 9 stations of the recently installed geoscope broadband network splitting was observed for about 15 records and is interpreted as being due to seismic anisotropy beneath the stations from every record it is possible to retrieve the fast direction at the station and the time difference between the fast and slow polarization directions some 30 observations of quasi linear s pulses also help constrain the direction of fast velocity synthetic seismograms were built to test the validity of the method used to retrieve the parameters of anisotropy typical time differences are about 1 sec corresponding for example to a 80 km thick zone of perfectly oriented pyrolite the fast directions obtained at 3 continental stations in america westford santa cruz heredia could be related to plate tectonic processes on the other hand the fast directions found in 3 oceanic islands kipapa papete la reunion differ strongly from the oceanic plate velocity direction there

3 d upper mantle shear velocity and attenuation from fundamental mode free oscillation data summary we present the results of global measurements of spheroidal and toroidal fundamental mode eigenfrequencies as well as fundamental spheroidal mode attenuation made from a data set consisting primarily of records from the geoscope network data processing includes variable filtering and account is taken of higher order effects along the dispersion branch local frequency and attenuation maps are first obtained as a function of angular order by a regionalization scheme based on the introduction of a correlation length in a second step after correction for shallow layer effects shear velocity maps at different depths in the mantle are obtained by simultaneous inversion of spheroidal and toroidal local frequency maps the results confirm the correlation of shear velocity with large scale tectonic features in the first 200 km of the mantle below that depth the maps are dominated by a degree 2 pattern whose power is peaked around a depth of 300 km shallower than inferred in the first 3 d global studies of the upper mantle and correlated with a degree 2 in attenuation thus confirming the results of our earlier study of degree 2 in the upper mantle we also report estimates of spherically averaged eigenfrequencies for fundamental spheroidal and toroidal mode branches

shear wave splitting in the mantle ps phases distribution with depth of mantle anisotropy is uncertain to better constrain this distribution we analyse records of the mantle ps phases at the grf array in germany two components of these phases sv and t are detected by stacking teleseismic records of events well distributed in azimuth and distance the detected t component signals are generated not so much by splitting of sv as by conversion from p to sh at the mantle discontinuities to reconcile these observations with previously known observations of splitting in sks at least three anisotropic layers are required the data favour a model with nne oriented fast velocity at the top of the mantle this is consistent with the results of other seismic studies beneath this layer the fast direction of anisotropy changes by around 90 our analysis suggests as well that a layer at the base of the mantle transition zone mtz is anisotropic this anisotropy might arise in the convective boundary layer between the upper and the lower mantle

seismic anisotropy in the d layer we present observations of diffracted sv svd for a path between the fiji tonga islands and the eastern coast of north america at distances greater than 110 observed features of s diffracted suggest that coupling between svd and shd can be ruled out as a first order effect for this path arrivals of shd are late relative to iasp91 travel times by about 10 s and those of svd are late relative to shd by 3 s for most records the slope of the log svd shd spectral ratio is around 3hz 1 in the range 0 06 0 15 hz a transversely isotropic low velocity layer in the lowermost mantle with a thickness of 200 300 km may account for most of the observed properties of svd

composition of the upper mantle geophysical tests of two petrological models the elastic properties of candidate mantle phases are used to test the viability of olivine rich pyrolitic and cao al 2o 3 rich eclogitic assemblages for the mantle high temperature adiabats for each phase of interest are constructed and compared to mantle seismic properties both pyrolitic and eclogitic assemblages satisfy the seismic properties between 200 and 400 km between 400 and 670 km depth an eclogitic assemblage yields a superior match to velocities and velocity gradients the 400 km seismic discontinuity may represent a chemical boundary between pyrolite and picritic eclogite piclogite or phase transformations in the olivine orthopyroxene components of a piclogitic assemblage containing about 16 olivine high velocity gradients in the transition zone may be explained by the transformation of ca rich cpx to majorite garnet seismic properties at the top of the lower mantle are consistent with pyrolite piclogite or perovskite implying that the 670 km discontinuity may be a chemical boundary

three dimensional structure of the upper mantle beneath the atlantic ocean inferred from long period rayleigh waves 1 group and phase velocity distributions group and phase velocity distributions of rayleigh wave fundamental modes are determined in the atlantic ocean using a method of regionalization without a priori constraints for a period range from 50 to 200 s eighty six direct source station paths are studied and 38 out of the best data from weidner s 1974 phase velocities are added to our data set for periods lower than 100 s the importance of shallow layer corrections is described in detail the general features of the distributions agree with global models and smaller scale heterogeneities are revealed a spatial resolution of 1500 km is needed to detect fine lateral heterogeneities in the atlantic the small lateral extent of young oceanic regions and strong lateral velocity gradients do not allow to quantify velocity versus age relationships using our data set on the other hand the regionalization without a priori constraints of phase velocities confirms the presence of a high velocity anomaly beneath the central part of the atlantic for periods greater than 100 s

evidence from core reflected shear waves for anisotropy in the earth s mantle seismologists have begun to recognize that the upper mantle at depth is anisotropic and that seismological data are not fully explicable if anisotropy is ignored1 5 on entering an anisotropic zone a plane shear wave splits into two orthogonally polarized phases6 an effect observed by several researchers7 notably by ando et al 5 8 9 the near vertical scs wave an s wave reflected once at the core mantle boundary is one of the most suitable phases to study the phenomenon of shear wave birefringence because it is isolated from other phases and because its particle motion is nearly horizontal and thus free from s to p conversion the split scs shear waves travel with different velocities in the anisotropic zone and consequently there is a time delay between their arrival at stations my aim was to determine this time delay and the azimuth of polarization using the scs records of two horizontal components from the japanese islands for a deep shock beneath the sea of okhotsk from the total of 41 stations 34 show a polarization of the maximum velocity phase in the nnw sse direction with an average time advance of 0 8s from the minimum velocity phase the relevant anisotropic zone seems to be deep seated

seismic anisotropy in the eastern united states deep structure of a complex continental plate we have analyzed shear wave splitting recorded by portable and permanent broadband and long period stations located in the eastern united states teleseismic shear waves sks skks and pks were used to retrieve the splitting parameters the orientation of the fast wave polarization plane and the delay time t in total 120 seismic events were processed allowing for more than 600 splitting measurements within the appalachians stations located in the western external part are characterized by t 1s and trending n50 70 e in the south and central regions and n30 40 e in the north closely following the trend of the orogenic belt in these areas the transition region between north and central is characterized by t 1 1 3 s and by e w trending that are at a high angle to the regional geologic trend measurements at two stations located in the eastern internal part of the belt indicate very weak anisotropy the large scale pattern of anisotropy is not consistent with that predicted for simple asthenospheric flow beneath the plate splitting along the southern and eastern margins of the continent is consistent with that expected for grenvillian deformation an alternative model of asthenospheric flow around the cratonic keel cannot be ruled out within the cratonic core the correlation between t and lithospheric thickness suggests a lithospheric anisotropy smaller length scale variations also argue for a significant contribution of lithospheric structures the fabric responsible for shear wave splitting may have formed during tectonic episodes that affected the eastern united states i e the grenville and appalachian orogenies and the subsequent rifting of the north atlantic ocean our observations in the western appalachians suggest that the anisotropy may be preserved since the grenvillian orogeny the absence of detectable splitting in the two stations in the eastern appalachians is attributed to the igneous intrusions related to the atlantic rifting the measurements in the transition between the northern and central southern appalachians constitute an intriguing anomaly whose e w have little obvious relation to the regional surface geology we suggest two possible causes 1 the local dominance of asthenospheric flow motivated by the proximity of a pervasive low velocity anomaly and 2 lithospheric deformation in a transcontinental strike slip fault zone active during the appalachian collision

measurements of mantle wave velocities and inversion for lateral heterogeneities and anisotropy 3 inversion lateral heterogeneity in the earth s upper mantle is investigated by inverting dispersion curves of long period surface waves 100 330 s models for seven different tectonic regions are derived by inversion of regionalized great circle phase velocity measurements from our previous studies we also obtain a representation of upper mantle heterogeneities with no a priori regionalization from the inversion of the degree 6 spherical harmonic expansion of phase and group velocities the data are from the observation of about 200 paths for love waves and 250 paths for rayleigh waves for both the regionalized and the spherical harmonic inversions corrections are applied to take into account lateral variations in crustal thickness and other shallow parameters these corrections are found to be important especially at low spherical harmonic order the trench region and fast velocities down to 250 km under shields below 200 km under the oceans both s velocity and s anisotropy support a model of small scale convection in which cold blobs detach from the bottom of the lithosphere when its age is large enough the spherical harmonic models clearly demonstrate a posteriori the relation between surface tectonics and s velocity heterogeneities in the first 250 km all shields are fast most ridges are slow below 300 km a belt of fast mantle follows the pacific subduction zones however at greater depths large scale heterogeneities that seem to bear no relationship to surface tectonics are observed the most prominent feature at 450 km is a fast velocity region under the south atlantic ocean smaller scale heterogeneities that are not related to surface tectonics are also mapped at shallower depths an anomalously slow region centered in the south central pacific is possibly linked to intense hot spot activity a very fast region southeast of south america may be related to subduction of old pacific plate between 200 and 400 km a belt of sv sh anisotropy follows part of the ridge and subduction systems indicating vertical mantle flow in these regions the spherical harmonic results open new horizons for the understanding of convection in the mantle perspectives for the improvement of the models presented are discussed

global patterns of azimuthal anisotropy and deformations in the continental mantle summary we present a summary of measurements of azimuthal anisotropy in the continental mantle based on the sks technique and performed mostly with the active participation of the authors the directions of polarization of the fast quasi shear wave and the time delays between the quasi shear waves are obtained at nearly 70 locations in all continents except antarctica these data are interpreted in terms of lattice preferred orientation of olivine which is caused by deformations in the mantle the depth interval responsible for anisotropy is unknown but the data suggest that it may reach at least 300 km the fast directions in sks do not show clear correlation with the fast directions of the teleseismic p at the same seismograph stations in the regions of present day convergence the fast direction of anisotropy usually aligns with the plate boundary this correlation implies that the direction of shortening is the same in the crust and the upper mantle in the regions of rifting the inferred direction of mantle flow usually aligns with the direction of extension in the crust outside the regions of recent tectonic activity we most likely observe a combined effect of frozen anisotropy in the subcrustal lithosphere and of recently formed anisotropy in the asthenosphere on a global scale in these regions there is a positive correlation between the absolute plate velocity directions and the fast directions of anisotropy the correlation is especially strong in central and eastern parts of north america a clear absence of any evidence of large scale azimuthal anisotropy in the data of long range refraction profiling for the upper 100 km of the mantle of that region implies that the effect in sks is generated mainly at greater depths in the asthenosphere orientation of olivine at these depths reflects recent and present day flow in the mantle rather than processes of a distant geologic past

azimuthal anisotropy in the lithosphere from observations of long period s waves summary the paper gives a detailed description of the technique for measuring azimuthal anisotropy which was suggested earlier by vinnik kosarev makeyeva 1984 and used by kind et ai 1985 the technique is based on the observations of long period converted phases like sks additionally we describe a generalization of this approach which makes use of long period s waves of arbitrary polarization both modifications were applied to records of the grf array in southern germany the results of this analysis if combined with the data on p velocities bamford 1977 suggest that the direction of the fast velocity in the lithosphere of the region varies with depth

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

regional and azimuthal dependence of phase velocities of mantle rayleigh waves in the pacific ocean abstract the lateral heterogeneity and azimuthal anisotropy of the upper mantle in the pacific ocean are investigated by analysing rayleigh wave phase velocities in the period range 60 200 s rayleigh waves at periods below 150 s travel fast in the direction parallel to the present day plate motion rayleigh wave velocities at longer periods show a similar azimuthal dependence in most parts of the pacific the average magnitude of the velocity anisotropy in the pacific ocean changes from 1 at 60 s to 0 4 at 200 s the azimuth independent phase velocities increase with the increase in age of the sea floor in the period range 60 120 s but the correlation is weak at longer periods the trade off between the lateral variation and the azimuthal variation is weak except for marginal places of the studied area

shear wave splitting as a diagnostic of variable anisotropic structure of the upper mantle beneath central fennoscandia abstract we analyze splitting of shear waves recorded during the svekalapko passive seismic experiment in south central finland to study fabrics of the mantle lithosphere of the precambrian region and thus to bring information into a debate on existence of plate tectonics or its forms in the early stage of continent formation geographical variations of the splitting parameters and their distinct dependence on direction of wave propagation through the upper mantle allow us to identify six domains of the central fennoscandian mantle lithosphere including the proterozoic archean transition and to model their fabrics by joint inversion of body wave anisotropic parameters fabrics of the archean mantle lithosphere can be approximated by a peridotite aggregate with lineation a dipping to the ne on the other hand anisotropy of the proterozoic mantle lithosphere is weaker and we model its fabric by the a c foliations dipping to the se we present a 3d self consistent anisotropic model of the proterozoic and archean upper mantle along the sw ne profile in the south central finland boundaries of inter growing wedges of the proterozoic and archean mantle lithospheres explain the longitudinal and shear wave propagation and polarization mantle xenolith ages surface wave tomography and location of the upper mantle reflectors we interpret the six anisotropic domains as fragments of mantle lithosphere retaining an old fossil olivine fabric which was created before these micro continents assembled

where can seismic anisotropy be detected in the earth s mantle in boundary layers during the last 30 years considerable evidence of seismic anisotropy has accumulated demonstrating that it is present at all scales but not in all depth ranges we detail which conditions are necessary to detect large scale seismic anisotropy firstly minerals must display a strong anisotropy at the microscopic scale and or the medium must be finely layered secondly the relative orientations of symmetry axes in the different crystals must not counteract in destroying the intrinsic anisotropy of each mineral and there must be efficient mechanisms of orientation of minerals and aggregates finally the strain field must be coherent at large scale in order to preserve long wavelength anisotropy part of shallow anisotropy can be related to the past strain field frozen in anisotropy however the deep anisotropy is due to the present strain field all these conditions are fulfilled only in boundary layers of convective mantle we review in this paper the seismic data sets which provide insight into the location at depth of large scale anisotropy from the d layer up to the lithosphere in addition to the well documented seismic anisotropy in the lithosphere and asthenosphere there is new evidence of seismic anisotropy in the upper 400 660 km and lower 660 900 km transition zones and in the d layer nonetheless the bulk of the lower mantle seems close to isotropy if we assume the hypothesis that seismic anisotropy is associated with boundary layers in convective systems these observations strongly suggest that the transition zone is a boundary layer which makes the pasage of matter between the upper and the lower mantle difficult however this general statement does not rule out flow circulation between the upper and lower mantles finally the geophysical mineral physics and geological applications are briefly reviewed an intercomparison between surface wave anisotropy and body wave anisotropy data sets is presented we discuss the scientific potential of seismic anisotropy and how it makes it possible to gain more insight into continental root deformation and geodynamics processes

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

seismic anisotropy in ocean basins resistive drag of the sublithospheric mantle we use finite element models to evaluate the effect of asthenospheric strain associated with the motion of an oceanic plate over a presumed stationary mantle on shear wave splitting observations modeled velocity profiles display a clear strain localization within a horizontal shear zone several tens of kilometers wide developed between an almost rigid mechanical lithosphere and a mildly deformed upper mantle for young oceanic lithosphere cooling results in migration of the maximum shear strain rate towards deeper levels and progressive widening of the shear zone shear strain accumulates with plate displacement the deformed layer thickness depends on plate age and velocity seismic anisotropy depends on the preferred orientations of olivine developed and thus on the strain field for a ray propagating vertically e g sks the fast wave will be polarized parallel to the flow direction i e the apm the delay time will depend on the thickness of the sheared layer and on its intrinsic anisotropy it will therefore increase away from the ridge at a progressively decreasing rate and depend on the absolute plate velocity the estimated anisotropies agree with seismic anisotropy measurements in old domains of ocean basins suggesting that away from mid ocean ridges this simple model may provide a first approximation of the mechanism responsible for the observed seismic anisotropy

seismic anisotropy in the deep mantle boundary layers and the geometry of mantle convection an attempt is made to explore the geodynamical significance of seismic anisotropy in the deep mantle on the basis of mineral physics the mineral physics observations used include the effects of deformation mechanisms on lattice and shape preferred orientation the effects of pressure on elastic anisotropy and the nature of lattice preferred orientation in deep mantle minerals in dislocation creep regime many of these issues are still poorly constrained but a review of recent results shows that it is possible to interpret deep mantle seismic anisotropy in a unified fashion based on the solid state processes without invoking partial melting the key notions are i the likely regional variation in the magnitude of anisotropy as deformation mechanisms change from dislocation to diffusion creep or superplasticity associated with a change in the stress level and or grain size in the convecting mantle with a high rayleigh number and ii the change in elastic anisotropy with pressure in major mantle minerals particularly in mg fe o the results provide the following constraints on the style of mantle convection i the sh sv anisotropy in the bottom transition zone and the sv sh anisotropy in the top lower mantle can be attributed to anisotropy structures lattice preferred orientation and or laminated structures caused by the horizontal flow in this depth range suggesting the presence of a mid mantle boundary layer due to partially layered convection ii the observed no significant seismic anisotropy in the deep mantle near subduction zones implies that deformation associated with subducting slabs is due mostly to diffusion creep or superplasticity and therefore slabs are weak in the deep mantle and hence easily deformed when encountered with resistance forces and iii the sh sv anisotropy in the cold thick portions of the d layer is likely to be due to horizontally aligned shape preferred orientation in perovskite plus magnesiowustite aggregates formed by strong horizontal shear motion in the recent past

seismic anisotropy beneath the continents probing the depths of geology seismic anisotropy beneath continents is analyzed from shear wave splitting recorded at more than 300 continental seismic stations anisotropy is found to be a ubiquitous property that is due to mantle deformation from past and present orogenic activity the observed coherence with crustal deformation implies that the mantle plays a major if not dominant role in orogenies no evidence is found for a continental asthenospheric decoupling zone suggesting that continents are coupled to general mantle circulation

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

constraints from seismic anisotropy on the nature of the lowermost mantle the d layer lies at the bottom of the earth s rocky mantle and separates it from the liquid metal alloy core this region extending from the core mantle boundary to a few hundred kilometres above fig 1 is geodynamically analogous to the more easily studied lithosphere at the top of the mantle the structure of d may reflect the style of lower mantle convection the nature of core mantle interaction and perhaps even the fate of subducting lithosphere1 observations of lithospheric seismic anisotropy have provided valuable insight into the nature of the upper mantle boundary layer but discussion of lower mantle seismic anisotropy has been somewhat contentious2 5 here we present evidence from seismic waves that have traversed the lowermost mantle beneath the caribbean region for a zone of seismic anisotropy below the d discontinuity which in this region lies 250 km above the core mantle boundary the anisotropy is most probably due to horizontal layering or aligned inclusions of a material with differing shear wave velocity if d is a graveyard for subducted lithosphere a plausible explanation of the anisotropy may be the contrast between cold lithospheric mantle and material that formerly constituted the oceanic crust which may have lower shear wave velocity owing to the presence of melt

evidence for anisotropy in the deep mantle beneath alaska we consider the possibility that the velocity structure of d is anisotropic the data we examined consist of seismograms from 9 deep japanese earthquakes recorded at wwssn receiver stations in north america the source receiver combinations span distances of 70 106 with associated s waves passing through d beneath alaska differential travel times of the s scd scs and sks phases are used to constrain the velocity structure in d shear waves refracted by d are observed beyond 72 2 and provide a sensitive measurement of the velocity structure in d beyond 93 the vertically polarized sv and horizontally polarized sh shear waves often appear distinctly split although at distances less than 89 the components are more nearly synchronous near 94 sh occurs as a double arrival sv in this range however remains a single arrival roughly synchronous with the second sh arrival we have been unable to reproduce these effects in isotropic model synthetics synthetics for transversely isotropic models have been computed that do match these waveforms the anisotropy was constrained to be only within d with a vertical symmetry axis we conclude that these observations may be explained by an anisotropic d layer the d discontinuity may be due to a transition to anisotropic mantle a few hundred kilometers above the core mantle boundary