two dimensional modelling of subduction zone anisotropy with application to southwestern japan summary we present a series of 2 d numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate we use a kinematically defined slab geometry approximating the subduction of the philippine sea plate beneath eurasia through finite element modelling we explore the effects of different rheological and thermal constraints e g a low viscosity region in the wedge corner power law versus newtonian rheology the inclusion of thermal buoyancy forces and a temperature dependent viscosity law on the velocity and finite strain field in the mantle wedge from the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material we then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory and compare the predictions to 500 previously published shear wave splitting measurements from seventeen stations of the broad band f net array located in southwestern japan although the use of different model parameters can have a substantial effect on the character of the finite strain field the effect on the average predicted splitting parameters is small however the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge the splitting observed in southwestern japan agrees well with the predictions of trench perpendicular flow in the mantle wedge along with b type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench

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

alaska megathrust 1 seismicity 43 years after the great 1964 alaska megathrust earthquake 1 the largest moment release during the 1964 mw 9 2 alaska earthquake was on the portion of the megathrust under the eastern kenai peninsula and the prince william sound the area is currently locked geodetically and corresponds to where the yakutat terrane is subducting in 2006 2009 a seismic array consisting of 34 broadband seismometers was deployed in the region an automatic algorithm was used to detect 12 563 local earthquakes using 13 months of data from the experiment of these 9427 good quality earthquakes could be relocated in a joint inversion for hypocenters and velocity structure they were then relocated by double difference to generate a final catalog of 8308 hypocenters these microearthquakes delineate a deeper steeply dipping wadati benioff zone contiguous with an 11 km wide seismic zone dipping at 3 between 20 km and 40 km depth along the 1964 earthquake s rupture zone focal mechanisms do not show interplate thrust faulting events but mostly normal faulting events with t axes generally parallel the slab dip direction indicating them to be intraslab seismicity the shallow narrow band of seismicity lies within the subducting yakutat terrane and right below the thrust zone possibly thrust faulting has not yet resumed yet in this early phase of the earthquake cycle thick subducted sediments overlying the yakutat terrane could also form a large strong contact zone on a relatively smooth plate boundary which does not favor seismic sliding on small patches but ruptures homogenously in great earthquakes

model for the origin of the yakutat block an accreting terrane in the northern gulf of alaska a composite oceanic and continental terrane the yakutat block is currently colliding with and accreting to north america in the northern gulf of alaska marine geophysical data indicate that the block has moved with the pacific plate during pliocene and quaternary time but are indeterminant about its previous movement history a model that explains some observed geologic and structural features of the block is that it originated as a result of subduction of the kula farallon spreading center beneath north america in this model the spreading center that formed the basaltic basement of the yakutat block was offset from a spreading center to the south by a left lateral transform fault 400 500 km long subduction of the northern spreading center about 45 m y ago sliced off part of the north america continental margin attached it to the kula plate and initiated northward movement of the composite terrane as the yakutat block during eocene and oligocene time migration of the southern spreading center toward north america truncated the yakutat block on the south and left a fossil fracture zone that is now seen as the transition fault the southern spreading center subducted about 25 m y ago and the yakutat block has since traveled with the pacific plate to southern alaska

shear wave polarization anisotropy and subduction induced flow in the mantle wedge of northeastern japan abstract shear wave polarization anisotropy in the central part of northeastern ne japan was investigated we analyzed s phases of 293 intermediate depth earthquakes recorded at 77 stations and obtained 1286 splitting parameters the obtained results show the systematic variation in splitting parameters across the arc the leading shear wave polarization directions fast direction obtained at stations in the western side of the study area are oriented nearly e w which is sub parallel to the direction of relative plate motion between the pacific plate and the north american okhotsk plate in contrast fast directions obtained at stations in the eastern side of the study area are oriented approximately n s stations in the western side have a larger average delay time compared to those in the eastern side these variations suggest that the nature of anisotropy is quite different between the eastern and the western sides both experimental studies on the deformation of olivine and numerical simulations of plate driven flow predict that the lattice preferred orientation of olivine causes the fast direction sub parallel to the flow direction for the conditions in the mantle wedge of ne japan we thus infer that the lattice preferred orientation of olivine generated by flow induced strain is the most likely candidate for anisotropy that produces splitting with e w fast direction three possible causes of anisotropy in the eastern side are considered deformation of water rich olivine in the mantle wedge trench parallel flow in the mantle wedge due to along strike dip variations in the slab and anisotropy in the crust and the slab we could not however determine which causes are dominant further systematic study is required to reveal the cause of anisotropy in the eastern side

effect of three dimensional slab geometry on deformation in the mantle wedge implications for shear wave anisotropy shear wave splitting observations from many subduction zones show complex patterns of seismic anisotropy that commonly have trench parallel fast directions three dimensional flow may give rise to trench parallel stretching and provide an explanation for these patterns of seismic anisotropy along strike variations in slab geometry produce trench parallel pressure gradients and are therefore a possible mechanism for three dimensional flow in this study we quantify the effects of variable slab dip curved slabs oblique subduction and slab edges on flow geometry and finite strain in the mantle wedge of subduction zones temperature dynamic pressure velocity and strain are calculated with high resolution three dimensional finite element models these models include temperature and stress dependent rheology and parameterized slab and trench geometry thick layers 20 60 km with strong trench parallel stretching are observed in the mantle wedge when slab geometry involves a transition to slab dip less than 15 or strong curvature in the slab in these cases strong trench parallel stretching develops when flow lines have an oblique to trench normal orientation this suggests that trench parallel seismically fast directions may not indicate trench parallel flow lines in systems with large along strike variations an oblique component of stretching is confined to a 20 30 km layer above the slab in systems with oblique subduction the effects of slab edges include strong toroidal flow and focusing in the mantle near slab edges and trench parallel flow that extends 50 100 km into the core of the mantle wedge

natural hazards triggered mud eruption the indonesian government ruled that the lusi mud eruption was triggered by drilling and held an oil company responsible instead a curved rock layer capping the mud reservoir may have amplified passing seismic waves and the trigger may have been natural

seismicity of the wrangell and aleutian wadati benioff zones and the north american plate along the trans alaska crustal transect chugach mountains and copper river basin southern alaska the configuration of the aleutian and wrangell wadati benioff zones in southern alaska between longitudes 142 w and 148 w is determined from regionally recorded earthquakes from 1971 to 1986 the earthquakes are located with a simple velocity model calibrated with seismic refraction measurements conducted for the trans alaska crustal transect tact west of the transect meridian approximately 145 20 w and beneath the chugach mountains shocks in the depth range 25 60 km compose the shallow part of the northwest to north dipping aleutian zone this zone can be traced updip and seaward beneath prince william sound east of the transect meridian and north of 61 n latitude the weakly active north northeast dipping wrangell zone extends to a depth of at least 100 km beneath the wrangell volcanoes and along strike for at least 150 km the two zones appear to be continuous at least in the depth range 20 45 km and may define adjacent limbs of a buckle in the subducted oceanic lithosphere north of 61 n focal mechanisms in the shallow aleutian zone are compatible with a uniform stress field characterized by in plane least compressive stress directed toward the steeply dipping deeper part of the aleutian zone north of cook inlet stress orientation in the wrangell zone is unresolved because of the low level of seismicity no focal mechanisms corresponding to slip on a gently northward dipping megathrust have been observed along the transect meridian a change of focal mechanism with depth near the top of the dipping seismic zone is inferred to reflect the interplate boundary beneath the chugach mountains this boundary coincides with a band of strong reflectors imaged by tact vertical seismic profiling the overriding plate appears to be aseismic beneath the northern chugach mountains but minor shallow seismicity occurs beneath the copper river basin first motion patterns for these shallow shocks suggest the most compressive stress within the overriding plate is oriented subhorizontally and approximately north south

active tectonics of the st elias orogen alaska observed with gps measurements 1 we use data from campaign and continuous gps sites in southeast and south central alaska to constrain a regional tectonic block model for the st elias orogen active tectonic deformation in the orogen is dominated by the effects of the collision of the yakutat block with southern alaska our results indicate that 37 mm yr of convergence is accommodated along a relatively narrow belt of n nw dipping thrust faults in the eastern half of the orogen with the present day deformation front running through icy bay and beneath the malaspina glacier near the bering glacier the collisional thrust fault regime transitions into a broad northwest dipping decollement as the yakutat block basement begins to subduct beneath the counterclockwise rotating elias block the location of this transition aligns with the gulf of alaska shear zone implying that the pacific plate is fragmenting in response to the yakutat collision our model indicates that the bering glacier region is undergoing internal deformation and could correspond to the final stage of offscraping and accretion of sediments from the yakutat block prior to subduction predicted block motions at the western edge of the orogen suggest that the crust is laterally escaping along the aleutian fore arc

new evidence for segmentation of the alaska subduction zone the purpose of this article is to provide additional evidence for segmentation of the subducted plate in alaska and to introduce a catalog of the relocated earthquakes for future studies we used the joint hypocenter determination method to relocate 14 099 subduction zone earthquakes that occurred from july 1988 to july 1998 and were located between 58 n and 65 n latitude the earthquake data were taken from the alaska earthquake information center catalog the selected earthquakes were divided into 16 blocks on the basis of their hypocentral locations and each block was relocated separately average epicenter shift was 3 8 km and average upward and downward depth shifts were 4 1 and 4 4 km respectively roughly the same number of earthquakes shifted upward 47 and downward 53 the overall change with respect to the initial locations is that the seismicity became more compact revealing details about the fine structure of the wadati benioff zone in particular we were able to identify more precisely the boundary between the kenai and mckinley blocks of the subducting plate in addition there is evidence for plate segmentation within the mckinley block manuscript received 15 november 2000

reconciling surface plate motions with rapid three dimensional mantle flow around a slab edge a central axiom of the theory of plate tectonics is that the surface plates are coupled to and carried along by convection in the earth s mantle the situation is rather different at subduction zones however margarete jadamec and magali billen now present calculations that quantify this variance they use three dimensional numerical models of buoyancy driven deformation with realistic slab geometry for the alaska subduction transform system to determine the surrounding mantle flow field they show that near a subduction zone edge mantle flow velocities can have magnitudes of greater than ten times the surface plate motions and that the complex mantle flow field is consistent with observations from seismic anisotropy these results imply local rapid transport of geochemical signatures through subduction zones the direction of tectonic plate motion at the earth s surface and the flow field of the mantle inferred from seismic anisotropy are well correlated globally implying large scale coupling between the mantle and surface plates but the fit is poor at subduction zones here three dimensional numerical models of buoyancy driven deformation for the alaska subduction transform system are used to find the origin of the regional decoupling of flow the results imply local rapid transport of geochemical signatures through subduction zones the direction of tectonic plate motion at the earth s surface and the flow field of the mantle inferred from seismic anisotropy are well correlated globally suggesting large scale coupling between the mantle and the surface plates1 2 the fit is typically poor at subduction zones however where regional observations of seismic anisotropy suggest that the direction of mantle flow is not parallel to3 4 5 6 7 and may be several times faster than6 plate motions here we present three dimensional numerical models of buoyancy driven deformation with realistic slab geometry for the alaska subduction transform system and use them to determine the origin of this regional decoupling of flow we find that near a subduction zone edge mantle flow velocities can have magnitudes of more than ten times the surface plate motions whereas surface plate velocities are consistent with plate motions8 and the complex mantle flow field is consistent with observations from seismic anisotropy5 the seismic anisotropy observations constrain the shape of the eastern slab edge and require non newtonian mantle rheology the incorporation of the non newtonian viscosity9 10 results in mantle viscosities of 1017 to 1018 pa s in regions of high strain rate 10 12 s 1 and this low viscosity enables the mantle flow field to decouple partially from the motion of the surface plates these results imply local rapid transport of geochemical signatures through subduction zones and that the internal deformation of slabs decreases the slab pull force available to drive subducting plates

oxygen isotope evidence for major fluid flow along the contact zone of the rum ultrabasic intrusion recent studies indicate that the rum tertiary ultrabasic intrusion formed in situ and was not emplaced as a fault bounded plug the suggestion that the main ring fault was the primary pathway for the flow of meteoric hydrothermal fluids on rum is therefore seriously flawed oxygen isotope evidence is presented indicating that the contact zone of the intrusion was the major pathway for meteoric fluids during cooling of the pluton 18 o depletions of over 12 correlate with hydrothermal alteration assemblages indicating that the bulk of the interaction with meteoric fluids took place at low temperatures 200 450 c

subduction of oceanic asthenosphere evidence from sub slab seismic anisotropy 1 the oceanic asthenosphere is characterized as a low viscosity channel down to 200 300 km depth separating the cold lithosphere from above and it is intimately linked to a layer of low seismic velocity and prominent seismic anisotropy observed globally beneath ocean basins while subduction of tectonic plates in convergent margins is well recognized the fate of oceanic asthenosphere remains enigmatic we demonstrate that subduction of the oceanic asthenosphere characterized by weak azimuthal anisotropy and strong radial anisotropy explains the essence of sub slab shear wave splitting patterns where the fast splitting direction changes from predominantly trench parallel or sub parallel under relatively steep subduction zones to frequently trench normal under shallow subduction zones to explain the observed splitting time the thickness of the subducted asthenosphere is estimated to be 100 50 km

an olivine fabric database an overview of upper mantle fabrics and seismic anisotropy we present a unique database of 110 olivine petrofabrics and their calculated seismic properties the samples come from a variety of the upper mantle geodynamic environments ophiolites subduction zones and kimberlites with a wide range of micro structures a phenomenological relationship is established between p and s wave seismic anisotropy and the degree of crystal alignment fabric strength seismic anisotropy increases rapidly at low fabric strength before reaching a near saturation level of 15 to 20 for p waves and 10 to 15 for s waves despite a large variation in the symmetry of fabric patterns the average seismic properties of the different fabric micro structural and geodynamic settings have similar anisotropies in both magnitude and symmetry hence it would seem possible to determine some measure of fabric strength from seismic anisotropy if the dimensions of the anisotropic region are known but not geodynamic environment or details of the petrofabric pattern a simple pattern of seismic anisotropy characterises the average sample of the database which has the following features the polarisation plane of the fastest s wave is parallel or sub parallel to the foliation plane the maximum shear wave splitting is parallel to the y structural direction in foliation plane and normal to the lineation the maximum of the p wave velocity is parallel to the high concentration of 100 axes which is sub parallel to the lineation the 100 orientation distribution has the greatest influence on the p wave seismic anisotropy the 100 and 001 orientation distributions have the greatest influence on the symmetry of s wave anisotropy although the magnitude of anisotropy is influenced by the distribution of all three principal axes although the database only contains olivine petrofabrics this statistical study clearly shows that seismic anisotropy can be used to deduce the orientation of the structural frame in the upper mantle 1998 elsevier science b v all rights reserved

natural examples of olivine lattice preferred orientation patterns with a flow normal a axis maximum tectonic plate motion is thought to cause solid state plastic flow within the underlying upper mantle and accordingly lead to the development of a lattice preferred orientation of the constituent olivine crystals the mechanical anisotropy that results from such preferred orientation typically produces a direction of maximum seismic wave velocity parallel to the plate motion direction this has been explained by the existence of an olivine preferred orientation with an a axis maximum parallel to the induced mantle flow direction in subduction zones however the olivine a axes have been inferred to be arranged roughly perpendicular to plate motion which has usually been ascribed to localized complex mantle flow patterns recent experimental work suggests an alternative explanation under conditions of high water activity a b type olivine preferred orientation may form with the a axis maximum perpendicular to the flow direction natural examples of such b type preferred orientation are however almost entirely unknown here we document widespread b type olivine preferred orientation patterns from a subduction type metamorphic belt in southwest japan and show that these patterns developed in the presence of water our discovery implies that mantle flow above subduction zones may be much simpler than has generally been thought

shear wave splitting continental keels and patterns of mantle flow in this study we investigated the origin of seismic anisotropy in the mantle beneath north america in particular we evaluated whether shear wave splitting patterns in eastern north america are better explained by anisotropy caused by lithospheric deformation anisotropy due to mantle flow beneath the lithosphere or a combination of both we examined new measurements of shear wave splitting from the missouri to massachusetts broadband seismometer array moma the north american mantle anisotropy and discontinuity experiment nomad as well as splitting parameters from several previous studies we developed a simple finite difference model that approximates mantle flow around a complex three dimensional continental lithospheric keel to evaluate potential anisotropy from mantle flow beneath the lithosphere in eastern north america we compared shear wave splitting observations to predicted splitting parameters calculated using this mantle flow model our results indicate that a significant portion of observed shear wave splitting in eastern north america can be explained by mantle flow around the continental keel however shear wave splitting patterns in a few regions of eastern north america indicate that a component of lithospheric anisotropy must exist particularly in regions containing the largest keel thicknesses for eastern north america as well as for splitting observations in australia europe and south america we favor a model in which anisotropy is controlled by a combination of both lithospheric deformation and subcontinental mantle flow

seismic anisotropy under central alaska from sks splitting observations 1 seismic anisotropy under central alaska is studied using shear wave splitting observations of sks waves recorded on the broadband experiment across the alaska range beaar 1999 2001 splitting results can be divided into two distinct regions separated by the 70 km contour of the subducting pacific plate waves that travel through the thicker mantle wedge show fast directions that are parallel to the strike of the slab these slab parallel directions appear to indicate along strike flow in the mantle wedge and splitting delay times increase with path length in the mantle wedge suggesting anisotropy of 7 9 0 9 the region of along strike flow corresponds to high seismic attenuation and hence high temperatures along strike flow here may be driven by secular shallowing of the slab driven by subduction of buoyant yakutat terrane crust or by toroidal flow around the east end of the aleutian slab waves traveling southeast of the 70 km contour sample the pacific plate and the nose of the mantle wedge they show fast directions that parallel the direction of plate motion these fast directions are most likely due to flow under the subducting pacific plate and or anisotropy within the subducting pacific lithosphere the high splitting delay times 0 8 1 7 s associated with these convergence parallel directions cannot be produced in the mantle wedge which is 10 30 km thick here thus anisotropy shows a sharp 90 change in fabric associated with the onset of high temperature wedge flow

effect of temperature on the b to c type olivine fabric transition and implication for flow pattern in subduction zones abstract the effect of temperature on the b to c type fabric transition in olivine was investigated by large strain shear deformation experiments a new series of deformation experiments was conducted at 1273 1373 k 2 0 gpa and strain rates ranging from 3 8 10 4 to 3 2 10 5 s 1 under water saturated conditions and combined with previous results at 1473 k the influence of temperature on deformation fabrics was investigated two types of olivine fabrics were identified one in which the olivine 0 0 1 axis is subparallel to the shear direction and the 0 1 0 plane is parallel to the shear plane b type and another one in which the 0 0 1 axis is parallel to the shear direction and the 1 0 0 plane is parallel to the shear plane c type the b type olivine fabric occurs at higher stress than the c type but the stress magnitude at which the b to c type fabric transition occurs decreases significantly with decreasing temperature the temperature dependence of the transition stress is parameterized using a range of exponential flow law that is constrained by the experimental study when the experimental results are extrapolated to lower temperatures the fabric transition between b and c type is predicted to occur at 1 mpa and 1000 k this prediction is consistent with the olivine fabrics in naturally deformed peridotites found in plate convergent regions the present study indicates that when some amounts of water 200 ppm h si are present the b type fabric will dominate in low temperature regions of the upper mantle whereas the c type or e type fabric will dominate at high temperatures the complex seismic anisotropy in subduction zones can be explained by changes in the dominant type of olivine fabrics caused by spatial variations in physical and chemical variables such as temperature and water content rather than by complexities in mantle wedge flow pattern

anisotropy and flow in pacific subduction zone back arcs we have obtained constraints on the strength and orientation of anisotropy in the mantle beneath the tonga southern kuril japan and izu bonin subduction zones using shear wave splitting in s phases from local earthquakes and in teleseismic core phases such as sks the observed splitting in all four subduction zones is consistent with a model in which the lower transition zone 520 660 km and lower mantle are isotropic and in which significant anisotropy occurs in the back arc upper mantle the upper transition zone 410 520 km beneath the southern kurils appears to contain weak anisotropy the observed fast directions indicate that the geometry of back arc strain in the upper mantle varies systematically across the western pacific rim beneath izu bonin and tonga fast directions are aligned with the azimuth of subducting pacific plate motion and are parallel or sub parallel to overriding plate extension however fast directions beneath the japan sea western honshu and sakhalin island are highly oblique to subducting plate motion and parallel to present or past overriding plate shearing models of back arc mantle flow that are driven by viscous coupling to local plate motions can reproduce the splitting observed in tonga and izu bonin but further three dimensional flow modeling is required to ascertain whether viscous plate coupling can explain the splitting observed in the southern kurils and japan the fast directions in the southern kurils and japan may require strain in the back arc mantle that is driven by regional or global patterns of mantle flow

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

olivine fabric transitions and shear wave anisotropy in the ryukyu subduction system recent shear wave splitting measurements from the fore arc region of the ryukyu subduction system show large magnitude 0 3 1 6 s trenchparallel splitting in both local and teleseismic phases the similarity of splitting parameters associated with shallow local s and teleseismic phases suggests that the source of anisotropy is located in the fore arc mantle one explanation for this pattern of shear wave splitting involves a transition from commonly observed high temperature olivine fabrics with flow parallel seismically fast directions to a flow normal b type olivine fabric in the cold fore arc mantle of the ryukyu wedge we test the b type fabric hypothesis by comparing observed splitting parameters to those predicted from geodynamic models that incorporate olivine fabric development the distribution of olivine fabric is calculated with high resolution thermomechanical models of the ryukyu subduction zone that include realistic slab geometry and an experimentally based wet olivine rheology we conclude that b type fabric can explain the magnitude and trench parallel orientation of deep local s phases that sample the core of the forearc mantle however our calculations show that b type fabric alone cannot account for large magnitude trench parallel splitting associated with teleseismic phases that sample the shallow tip of the fore arc mantle alternative models for trench parallel teleseismic splitting in the shallow tip of the fore arc mantle involve the addition of crustal or slab anisotropy and highly anisotropic foliated antigorite serpentinite 2008 elsevier b v all rights reserved

seismic anisotropy of oceanic upper mantle shear wave splitting methodologies and observations we develop methodologies to obtain accurate measurements of shear wave splitting and apply these techniques to examine the pattern of oceanic upper mantle anisotropy to obtain high quality estimates of receiver splitting at island stations we devise a stacking method that finds the optimum splitting parameters and 95 error region for teleseismic shear wave phases from a suite of earthquake events to obtain additional measurements in oceanic regions we develop techniques to measure splitting parameters and errors for ss phases that sample upper mantle anisotropy at their bounce points however we find the data are often of low resolution and anomalous characteristics are sometimes found that make splitting difficult to interpret ten second splitting of ss is observed across the banjo seismic array but we cannot unambiguously attribute this signal to mantle anisotropy at the bounce the receiver splitting methods are used to assess the adequacy of a two anisotropic layer model for the pacific region with fast polarization azimuth in the lithosphere oriented in the fossil spreading direction and in the asthenosphere oriented in the absolute plate motion direction this model has been proposed to explain surface wave data in the pacific ocean but our splitting results demonstrate that oceanic anisotropy patterns are more heterogeneous than would be predicted while the island splitting measurements could reflect the influence of individual hotspot upwellings hotspot effects do not appear to be universally dominant we propose that splitting observations alternatively indicate broad scale differences in the underlying character of oceanic upper mantle anisotropy associated with coherant patters of lithospheric structure and asthenospheric flow in particular splitting surface wave models and regional studies all support a model where lithospheric anisotropy throughout the south pacific has been erased or reoriented toward the absolute plate motion direction whereas more limited observations in the north pacific indicate that the fossil lithospheric signature appears to be preserved

sks splitting beneath alaska regional variability and implications for subduction processes at a slab edge abstract many aspects of subduction zone geodynamics remain poorly understood including the pattern of flow in the mantle surrounding subducting slabs because of the relationship between deformation and seismic anisotropy measurements of shear wave splitting in subduction zones can shed light on this pattern of flow while observations of splitting in subduction zones are numerous and robust the source and geodynamical implications of the anisotropy remain imperfectly constrained the alaska subduction zone associated with the subduction of the pacific plate beneath north america is a particularly complex tectonic setting dramatic lateral variations in slab morphology seismicity and volcanic character are present beneath the region and previous splitting observations and modeling studies have suggested the possibility of significant along strike mantle flow in the region here we present sks splitting measurements from more than 50 stations of the permanent broadband ak network located throughout alaska splitting patterns show significant regional variability but proximate stations often exhibit similar splitting behavior stations in southern alaska exhibit nearly trench normal fast directions while in southeasternmost alaska fast directions trend parallel to the plate boundary a group of stations in the eastern part of the array exhibit predominantly null splitting while stations to the northwest of this group tend to exhibit 1 s or more of splitting with fast directions that are roughly ne sw stations east of the kenai peninsula region exhibit complicated splitting patterns that suggest complex anisotropy we discuss a number of potential factors that may contribute to the complexity in splitting patterns observed in alaska including frozen anisotropic structure in the overriding plate lithosphere asthenospheric shear due to absolute plate motion trench parallel flow in the mantle wedge two dimensional entrained flow beneath the slab and mantle upwelling at the slab edge

new gps constraints on the motion of the yakutat block global positioning system gps measurements were made at yakutat on the yakutat terrane of southern alaska to investigate the motion of the yakutat block with respect to the north american plate and to help constrain motion along the fairweather fault the velocity of yakutat derived from the gps data is 44 1 1 9 mm yr toward n37 4 w relative to stable north america the magnitude of this velocity is similar to that of the pacific plate predicted by nuvel1a although there is a significant difference in the azimuth of these two vectors the motion of yakutat relative to north america is almost exactly parallel to the strike of the fairweather fault suggesting that most deformation inboard of yakutat is right lateral strike slip on the fairweather fault or faults parallel to it and that significant motion normal to the fairweather fault occurs offshore of yakutat the gps velocity at yakutat is also used to help constrain the slip rate and locking depth of the fairweather fault

mapping the megathrust beneath the northern gulf of alaska using wide angle seismic data in the northern gulf of alaska and prince william sound we have used wide angle seismic reflection refraction profiles earthquake studies and laboratory measurements of physical properties to determine the geometry of the prince william and yakutat terranes the aleutian megathrust and the subducting pacific plate in this complex region the yakutat terrane is underthrust beneath the prince william terrane and both terranes are interpreted to be underlain by the pacific plate wide angle seismic reflection refraction profiles recorded along five seismic lines are used to unravel this terrane geometry modeled velocities in the upper crust of the prince william terrane to 18 km depth agree closely with laboratory velocity measurements of orca group phyllites and quartzofeldspathic graywackes the chief components of the prince william terrane to hydrostatic pressures as high as 600 mpa 6 kbar a landward dipping reflector at depths of 16 24 km is interpreted as the base of the prince william terrane this reflector corresponds to the top of the wadati benioff zone seismicity and is interpreted as the megathrust immediately beneath the megathrust is a 4 km thick 6 9 km s refractor which we infer to be the source of a prominent magnetic anomaly and which is interpreted by us and previous workers to be gabbro in eocene age oceanic crust of the underthrust yakutat terrane wide angle seismic data magnetic anomaly data and tectonic reconstructions indicate that the yakutat terrane has been underthrust beneath the prince william terrane for at least a few hundred kilometers wide angle seismic data are consistent with a 9 to 10 landward dip of the subducting pacific plate beneath the outer shelf and slope distinctly different from the inferred average 3 to 4 dip of the overlying 6 9 km s refractor and the wadati benioff seismic zone beneath the inner shelf our preferred interpretation of the geophysical data is that one composite plate composed of the pacific plate of a fairly uniform thickness and the yakutat plate of varying thickness is subducting beneath southern alaska

plate tectonics and net lithosphere rotation over the past 150 my we have developed an improved model of global digital palaeo plate boundaries and plate motion to describe the distribution and history of plates since the late jurassic from this history we computed net lithospheric rotation nr through time confirming the so called westward drift but only for the past 30 myr the nr has significantly smaller magnitudes 0 13 my past 5 my than for some other plate models it averages to 0 11 0 03 my for the past 50 my with a small but systematic increase toward the present the westward drift seen only for the past 30 my is attributed to the increased dominance of a steadily growing and accelerating pacific plate nr shows peaks with time but only an early tertiary peak of 0 33 my when the indian plate was undergoing the largest known acceleration deceleration can be interpreted with some confidence we find a linear decreasing trend in net rotation over the past 150 my but attribute this trend to increasing reconstruction uncertainties back in time as subduction consumed more than half of the oceanic crust since the jurassic after removing a linear time trend we find a nr average of about 0 12 my for the past 150 my

trench parallel flow and seismic anisotropy in the mariana and andean subduction systems shear wave splitting measurements above the mantle wedge of the mariana1 and southern andean2 3 subduction zones show trench parallel seismically fast directions close to the trench and abrupt rotations to trench perpendicular anisotropy in the back arc these patterns of seismic anisotropy may be caused by three dimensional flow associated with along strike variations in slab geometry1 2 3 4 5 the mariana and andean subduction systems are associated with the largest along strike variations of slab geometry observed on earth6 7 and are ideal for testing the link between slab geometry and solid state creep processes in the mantle here we show with fully three dimensional non newtonian subduction zone models that the strong curvature of the mariana slab and the transition to shallow slab dip in the southern andes give rise to strong trench parallel stretching in the warm arc and warm back arc mantle and to abrupt rotations in stretching directions that are accompanied by strong trench parallel stretching these models show that the patterns of shear wave splitting observed in the mariana and southern andean systems may be caused by significant three dimensional flow induced by along strike variations in slab geometry

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

unusual mantle poisson s ratio subduction and crustal structure in central alaska 1 subduction zones alter the upper plate in a variety of ways including metasomatism via slab derived fluids and accretion of exotic terranes these processes should produce distinctive seismic signatures as seismic velocities are sensitive to compositional variations a new method is developed to estimates depths to interfaces and poisson s ratio vp vs ratio for dipping layered structures by stacking teleseismic receiver functions it also estimates crustal vp vs and thickness using this method we analyze a broadband data set from the central alaska subduction zone to seek evidence for mantle wedge composition crustal thickness varies from 27 to 45 km and crustal vp vs indicates an intermediate to mafic composition expected for these accreted island arc terranes inversions for mantle structure confirm the presence of a 15 20 km thick low velocity zone atop of the downgoing plate to 130 km depth perhaps subducted crust of an exotic terrane vp vs in the mantle wedge allows for 15 15 serpentinization where the slab is 80 km deep and the wedge is cold where the slab is deeper vp vs for the mantle wedge is unusually low 1 7 this value is lower than predicted for any common mantle mineral and may indicate significant quantities of quartz provided that past physical properties measurements are not grossly in error these results suggest that the mantle wedge in subduction zones can occasionally differ substantially in bulk composition from normal subcrustal mantle perhaps because collisions of exotic terranes cause tectonic mixing of crustal material into the mantle

mantle dynamics and seismic anisotropy abstract observations of seismic anisotropy yield some of the most direct constraints available on both past and present day deformation in the earth s mantle insight into the character of mantle flow can also be gained from the geodynamical modeling of mantle processes on both global and regional scales we highlight recent progress toward understanding mantle flow from both observations and modeling and discuss outstanding problems and avenues for progress particularly in the integration of seismological and geodynamical constraints to understand seismic anisotropy and the deformation that produces it to first order the predictions of upper mantle anisotropy made by global mantle circulation models match seismological observations well beneath the ocean basins but the fit is poorer in regions of greater tectonic complexity such as beneath continental interiors and within subduction systems in many regions of the upper mantle models of anisotropy derived from surface waves are seemingly inconsistent with shear wave splitting observations which suggests that our understanding of complex anisotropic regions remains incomplete observations of anisotropy in the d layer hold promise for improving our understanding of dynamic processes in the deep earth but much progress remains to be made in characterizing anisotropic structure and relating it to the geometry of flow geochemical heterogeneity or phase transitions major outstanding problems related to understanding mantle anisotropy remain particularly regarding the deformation and evolution of continents the nature of the asthenosphere subduction zone geodynamics and the thermo chemical state of the lowermost mantle however we expect that new seismological deployments and closer integration of observations with geodynamical models will yield rapid progress in these areas

shear wave splitting to test mantle deformation models around hawaii seismic anisotropy allows us to study mantle deformation and it can thus help to constrain mantle flow in the vicinity of hotspots hypotheses for the cause of seismic anisotropy in this environment include the parabolic asthenospheric flow paf model radial flow from a mantle plume impinging on a moving lithosphere is dragged by the plate in the direction of absolute plate motion in map view this gives a parabolic pattern of flow opening in the direction of plate motion we present new shear wave splitting observations from land and ocean stations around the hawaiian islands that can be explained by the parabolic flow model the observations suggest asthenospheric anisotropy under the hawaiian islands which may be explained if dislocation creep persists to deeper depths there than in other regions perhaps due to the higher temperatures near hotspots

insights into mantle structure and flow beneath alaska based on a decade of observations of shear wave splitting sks shear wave splitting measurements from three program for array seismic studies of the continental lithosphere experiments broadband experiment across the alaska range alaska receiving cross transect of the inner core and multidisciplinary observations of subduction which form a north south transect across alaska show a remarkably simple pattern of two large anisotropy domains in the northern domain extending from the 70 km contour of the subducting pacific plate north to the arctic ocean fast directions are consistently in the ne sw direction these directions are essentially parallel to the absolute plate motion direction in northern alaska and parallel to the strike of the subducting plate above the mantle wedge suggesting that they represent some combination of plate scale asthenospheric flow in the upper mantle and flow along the subducting plate in the mantle wedge a strong wedge component beneath the alaska range is required to explain systematics of splitting delay times in the southern domain which extends south from the 70 km depth contour to the subducting plate fast directions are in the nw se direction a 90 rotation from the northern domain these fast directions are parallel to the dip of the subducting plate in the direction of convergence and represent entrained flow beneath the subducting slab the pacific plate absolute motion approximately parallels local convergence two major factors seem to control flow in these regions absolute plate motion in the north and the subduction of the pacific plate in the south although both subduction driven wedge flow and absolute plate motion contribute to the southern part of the northern regime

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

subduction of oceanic asthenosphere a critical appraisal in central alaska song and kawakatsu 2012 have shown that the sub slab fast splitting pattern observed in most subduction zones can be a direct consequence of subduction of the oceanic asthenosphere that has strong radial anisotropy this model not only explains the non intuitive trench parallel splitting pattern in most of subduction zones but also predicts the trench normal behavior fast polarization direction sub parallel to the absolute plate motion of the incoming plate observed in several shallow subduction zones the general validity of such a scenario is crucial to fundamental understandings of the development of mantle anisotropy in sub lithosphere or and sub slab conditions the nature and formation of oceanic asthenosphere as well as the flow pattern and mass transport near subduction zones to validate this scenario we examine sks splitting patterns observed across the fore arc in central alaska here the fast splitting direction varies from plate motion sub parallel near the trench to mostly trench parallel beyond the 100 km slab isodepth contour while being strongly variable in between after taking into account the rotation of anisotropy symmetry in the oceanic asthenosphere with respect to the local plate motion obliquity and down dip variations in the slab dip we reproduce a general 90 degree switch in fast splitting direction as well as the back azimuth dependent splitting pattern across the entire fore arc the current validation further augments the idea that apart from anisotropy in the mantle wedge and the subducting slab subduction of the oceanic asthenosphere is likely to be the dominant source of seismic anisotropy in central alaska and potentially in many subduction zones furthermore this result also provides alternative views to models of seismic anisotropy in the mantle wedge and on the length scale in which the 3d mantle flow may be important