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

a comparison between the tectonic stress measured in situ and stress parameters from induced seismicity at monticello reservoir south carolina during may and early june 1979 five three component digital seismographs were deployed near the monticello reservoir south carolina to determine the stress drops of reservoir induced earthquakes and compare them with earlier hydrofracture stress measurements of zoback and hickman this issue records for 327 events were obtained 10 of these events each recorded on four or more stations were analyzed to determine stress drops depths for the 10 events ranged from 0 5 to 1 4 km fault plane orientations were calculated from an inversion of the body wave displacements individual radiation pattern corrections were applied in calculating the moments at each station final estimates of seismic moment were obtained by averaging logarithmically the estimates from each component for each event an arithmetic moment was obtained from the individual estimates of the source radius and stress drop was then calculated once from the final estimates of seismic moment and source radius for each event moments of the 10 digitally recorded earthquakes ranged from 2 5 1017 to 2 7 1018 dyne cm for the p wave data and from 4 6 1017 to 4 6 1018 dyne cm for the s wave data the stress drops s wave data only generally ranged from about 0 2 to 4 bars four earthquakes were recorded by a strong motion instrument at the dam site magnitudes m for these events ranged from 2 8 to 3 0 moments taken as averages of the values from each component ranged from 1 7 to 3 4 1020 dyne cm and the stress drops from 13 to 92 bars large differences between the corner frequencies of the two components factors of 2 5 and 4 for two of the events correspond to a large relative uncertainty for these stress drop values nevertheless tens of bars of stress are apparently released for the largest events in the monticello sequence but only a few bars are released for the smaller magnitude events in comparison the in situ stress data measured in two holes give an estimate of shear stress of about 45 bars at 0 2 km depth and 30 to 70 bars at 1 km depth in the northern hole and from 15 20 bars near the surface to almost 70 bars at depth in the southern hole the variation in shear stress inferred from in situ measurements at the deeper levels indicates an inconsistency between the thrust faulting mechanisms of the earthquakes and the state of stress which is more consistent with thrust faulting at 0 2 km depth and strike slip or normal faulting below thus the estimates of shear stress from in situ measurements are of the same order as the stress drops of the m 2 8 3 0 events but about 1 order of magnitude larger than the stress drops obtained for the smaller events

a kinematic model for recrystallization and texture development in olivine polycrystals the interpretation of seismic anisotropy in the mantle requires a knowledge of the relationship between the lattice preferred orientation lpo of crystals and the convective flow field in order better to understand this link we present a model for the evolution of lpo in olivine aggregates that deform by both intracrystalline slip and dynamic recrystallization dynamic recrystallization depends on the dislocation density of the grains which is a function of the applied local stress grains with a large density of dislocations lower their bulk strain energy by nucleating strain free sub grains at a rate proportional to a dimensionless nucleation parameter grains with high energy are then invaded by grains with low energy by grain boundary migration at a rate proportional to a dimensionless grain boundary mobility m the value of is constrained by observed lpo patterns in experimentally deformed olivine aggregates and m is constrained by the temporal evolution of the strength of the lpo for m 125 75 and 3 the model predictions agree well with the experimental results numerical calculations of lpo using our model are significantly faster than those based on viscoplastic self consistent or equilibrium based theories making the model especially suitable for applications for complex convective flows

three dimensional instantaneous mantle flow induced by subduction 1 we conduct three dimensional subduction experiments by a finite element approach to study flow around slabs which are prescribed based on a transient stage of upper mantle subduction from a laboratory model instantaneous velocity field solutions are examined focusing on the toroidal vs poloidal components as a function of boundary conditions plate width and viscosity contrast between slab and mantle we show how the slab to mantle viscosity ratio determines the strength of toroidal flow and find that the toroidal flow component peaks for slab mantle viscosity ratios 100 independent of slab width or geometry

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

constraints on upper mantle viscosity from the flow induced pressure gradient across the australian continental keel 1 the thickness of continental lithosphere varies considerably from tectonically active to cratonic regions where it can be as thick as 250 300 km embedded in the upper mantle like a ship when driven to move by a velocity imposed at the surface a continental keel is expected to induce a pressure gradient in the mantle we hypothesize that the viscosity of the asthenosphere or the shear coupling between lower lithosphere and asthenosphere should control this pressure effect and thus the resulting dynamic topography we perform three dimensional finite element calculations to examine the effects of forcing a continental keel by an imposed surface velocity with the australian region as a case study when the upper mantle is strong but still weaker than the lower mantle positive dynamic topography is created around the leading edge and negative dynamic topography is created around the trailing edge of the keel which is measurable by positive and negative geoid anomalies respectively for a weak upper mantle the effect is much reduced we analyze geoidal and gravity anomalies in the australian region by spatiospectral localization using slepian functions the method allows us to remove a best fit estimate of the geographically localized low spherical harmonic degree contributions regional geoid anomalies thus filtered are on the order of 10 m across the australian continent with a spatial pattern similar to that predicted by the models the comparison of modeled and observed geoid anomalies places constraints on mantle viscosity structure models with a two layer mantle cannot sufficiently constrain the ratio of viscosity between the upper and lower mantle the addition of a third weak upper mantle layer an asthenosphere amplifies the effects of keels our three layer models with lower mantle viscosity of 3 1022 pa s suggest that the upper mantle asthenosphere is 300 times weaker than the lower mantle while the transition zone 400 670 km depths has a viscosity varying between 1021 and 1022 pa s

upper mantle seismic anisotropy beneath antarctica and the scotia sea region summary recent investigations on shear wave splitting from recordings of permanent and temporary antarctic seismological stations have lead to a greater understanding of the upper mantle dynamics of the scotia sea region and the continental margin in the eastern weddell sea in terms of their tectonic evolution the analysis of shear wave splitting from teleseismic core sks skks pks and direct s waves reveals the seismic anisotropy and the strain field of the upper mantle similar to the caribbean anisotropy structures in the antarctic peninsula and scotia sea regions are assumed to be influenced by mantle flows in easterly directions around the subducting nazca plate in general anisotropy polarization directions in the scotia sea do not contradict this hypothesis with polarizations oriented nearly e w and therefore aligning with the suggested mantle flow patterns anisotropy strengths decrease from delay times of dt 1 8 s pmsa palmer station in the west towards the east with delay times of dt 0 3 s beneath hope south georgia and cand candlemas south sandwich islands nevertheless a lithospheric and therefore fossil origin cannot be ruled out only the exceptionally high delay times at pmsa probably originate in part from recent asthenospheric flow around the subduction slab of the former phoenix plate beneath the northwestern margin of the antarctic peninsula the continental margin of western dronning maud and coats land plays a crucial role in understanding the early processes during the break up of gondwana upper mantle seismic anisotropy with delay times well over dt 1 s in this region gives new constraints on ancient deformation processes during break up and former episodes two layer modelling reveals archaean anisotropy in the upper layer corresponding well to polarization directions of the south african kaapvaal craton lower layer anisotropy is assumed to have been created during early gondwana rifting stages

comparison of azimuthal seismic anisotropy from surface waves and finite strain from global mantle circulation models summary we present global models of strain accumulation in mantle flow to compare the predicted finite strain ellipsoid fse orientations with observed seismic anisotropy the geographic focus is on oceanic and young continental regions where we expect our models to agree best with azimuthal anisotropy from surface waves finite strain derived models and alignment with the largest fse axes lead to better model fits than the hypothesis of alignment of fast propagation orientation with absolute plate motions our modelling approach is simplified in that we are using a linear viscosity for flow and assume a simple relationship between strain and anisotropy however results are encouraging and suggest that similar models can be used to assess the validity of assumptions inherent in the modelling of mantle convection and lithospheric deformation our results substantiate the hypothesis that seismic anisotropy can be used as an indicator for mantle flow circulation derived models can contribute to the establishment of a quantitative functional relationship between the two

the taiwan ryukyu subduction collision complex folding of a viscoelastic slab and the double seismic zone 1 the termination of the ryukyu trench against eurasia and the oblique subduction of the philippine sea plate create a subduction collision complex offshore taiwan which has not previously been elucidated in detail we combine traveltime data from the seismic networks in taiwan and japan to better illuminate how the subducting ryukyu slab deforms in this subduction collision zone more than 5000 events recorded by both networks were relocated with the double difference method using an optimal regional one dimensional velocity model the offshore seismicity indicates that the double seismic zone with a gap of 15 20 km exists in the subducting slab in the depth range of 40 80 km focal mechanisms suggest that the double seismic zone is caused by east west compression resulting from oblique convergence the improved hypocentral locations for the first time reveal folding of the slab into a horizontal curvature larger in magnitude than and opposite in sign to that of the ryukyu trench in the depth range 50 100 km the anomalous curvature together with the focal mechanisms suggests that the slab folds against the eurasian lithosphere and that this deformation cannot be fully elastic we model this deformation mode as the developing instability of a viscoelastic maxwell layer embedded in a viscous medium the characteristic wavelength of the instability i e 250 km is consistent with folding of a slab whose viscosity is 100 times higher than that of the surrounding mantle for an along strike elastic membrane strain as small as 0 01 or more than 3 orders of magnitude higher if 5 elastic strain is allowed

mantle flow deflected by interactions between subducted slabs and cratonic keels subducting slabs can influence mantle flow but the importance of neighbouring continental cratons is little understood geodynamical modelling constrained by seismic data that identify regions of mantle flow beneath the caribbean south american plate margin shows that the deep rooted south american craton acts to deflect and enhance mantle flow into a narrow channel

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

reservoir induced earthquakes at sainte marguerite 3 quebec canada located some 90 km north of the town of sept iles quebec the sainte marguerite 3 sm 3 hydroelectric project was one of the most ambitious in eastern north america in recent years the project included the construction of a 171 m high dam that created a 140 km long reservoir some 120 m deep near the dam toe although the reservoir is located in a weakly seismic area of the canadian shield reservoir induced earthquakes occurred a first event mn 1 3 was recorded in september 1999 17 months after the reservoir started filling between october 2000 and march 2001 four additional events were detected by the canadian national seismograph network the largest being mn 1 7 during a field survey that took place between july and october 2001 276 events of magnitude mn 1 2 were detected within 30 km of a three component seismograph single station locations confirmed by agreement between p phase and s phase azimuths were determined for 182 events most earthquakes fall in four areas with estimated f

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

mantle flow at a slab edge seismic anisotropy in the kamchatka region the junction of the aleutian island and the kamchatka peninsula defines a sharp turn in the boundary of the pacific and north american plates terminating the subduction zones of the northwest pacific the regional pattern of shear wave birefringence near the junction indicates that trench parallel strain follows the seismogenic benioff zone but rotates to trench normal beyond the slab edge asthenospheric mantle is inferred to flow around and beneath the disrupted slab edge and may influence the shallowing dip of the benioff zone at the aleutian junction

shear wave splitting at the edge of the ryukyu subduction zone abstract intraslab events recorded by ocean bottom seismometers in the okinawa trough provide an extended depiction of shear wave splitting in the southwest section of the ryukyu subduction zone at 100 200 km from the western edge of the subduction system we observed trench normal fast polarization direction in the back arc compatible with 2d slab or rifting driven corner flow towards the edge the fast directions are sub parallel to the trench in the arc back arc region and rotate to trench normal within 50 km of the edge splitting constrained by land stations with paths mostly in the mantle wedge exhibits similar trench normal fast directions in the subduction edge zone further inland the dominant component of fast directions becomes roughly parallel to the taiwan orogenic fabric splitting of p to s phases converted at the moho of the okinawa trough and of s phases from shallow events suggest that crustal anisotropy may affect the measured splitting but the observed pattern reflects predominantly mantle anisotropy the variation in splitting along the okinawa trough cannot be explained by a b type a type olivine fabric transition in the mantle wedge it may indicate the presence of an along arc flow in the mantle wedge towards the edge where it is blocked and deflected by the eurasian lithosphere this scenario bolsters previous studies suggesting a significant impact of eurasian lithosphere on the dynamics of the ryukyu subduction system

development of mantle seismic anisotropy during subduction induced 3 d flow 1 the dynamics of subduction can be indirectly constrained by studying the induced mantle flow however inferring the circulation of the mantle around subducting plates from the interpretation of shear wave splitting patterns remains elusive we calculated the strain induced lattice preferred orientation lpo developed in 3 d models of subduction where retreat motions are maximized and found that in the mantle layer entrained with the downgoing slab the seismic anisotropy is trench perpendicular and becomes trench parallel deeper where the toroidal flow accommodates slab retreat synthetic sks splitting shows that in the fore arc slab rollback favors trench parallel polarization of the fast shear wave component while plate advance enhances trench perpendicular seismic anisotropy the interference between these two competing mechanisms yields subslab delay times of 0 5 to 1 3 sec comparable with those observed at most natural subduction zones the magnitude of the subslab trench parallel splitting is independent of the rate at which the slab migrates instead it is proportional to the amount of retreat

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

structure of the upper mantle and transition zone beneath southeast asia from traveltime tomography 1 tomographic images of the mantle beneath east asia were obtained from the inversion of traveltime data from global and regional seismograph networks and from temporary arrays on and around the tibetan plateau our results are consistent with previous studies but the unprecedented resolution of mantle heterogeneity provides new insight into the large scale tectonic framework of the continental india asia collision in the western part of the study region and subduction of the oceanic lithosphere in the east in the realm of continental collision west of 100 e a relatively slow p wave speed characterizes the upper mantle beneath much of the tibetan plateau but the wave speed is high beneath cratonic india the southern and western part of the tibetan plateau hindu kush and the tian shan in the subduction realm east of 110 e the main structures are i pronounced low wave speed anomalies at a depth of between 100 and 400 km beneath asia s southeastern seaboard and the back arc regions of ongoing subduction ii narrow fast anomalies in the upper mantle beneath major subduction zones and iii widespread fast anomalies at a depth of 500 700 km beneath the sea of japan the northern part of the philippine sea plate and southeastern china if the latter anomalies represent stagnant slabs their fragmented nature and large lateral extent suggest that they are produced by different episodes of subduction beneath western pacific island arcs along the old se margin of asia or during the mesozoic collision of cratonic units in southeast asia attribution to ancient subduction systems implies that slab fragments can reside in the transition zone for at least several tens of millions of years shallow slow anomalies beneath the red river fault region connect to deep anomalies beneath the south china fold belt and south china sea suggesting a causal relationship between the evolution of the continental lithosphere of sw china and deeper mantle processes between the collision and the subduction realms tomography reveals high wave speed continental roots beneath the western part of the north china craton ordos block and the south china or yangtze craton sichuan basin to a depth of 300 km

hydrothermal sediments as a potential record of seawater nd isotope compositions the rainbow vent site 36 14 n mid atlantic ridge 1 geochemical compositions and sr and nd isotopes were measured in two cores collected 2 and 5 km from the rainbow hydrothermal vent site on the mid atlantic ridge overall the cores record enrichments in fe and other metals from hydrothermal fallout but sequential dissolution of the sediments allows discrimination between a leach phase easily leachable and a residue phase refractory the oxy anion and transition metal distribution combined with rare earth element ree patterns suggest that 1 the leach fraction is a mixture of biogenic carbonate and hydrothermal fe mn oxy hydroxide with no significant contribution from detrital material and 2 99 5 of the ree content of the leach fraction is of seawater origin in addition the leach fraction has an average 87sr 86sr ratio indistinguishable from modern seawater at 0 70916 although we lack the nd value of present day deep water at the rainbow vent site we believe that the ree budget of the leach fraction is predominantly of seawater origin we suggest therefore that the leach fraction provides a record of local seawater nd values nd isotope data from these cores span the period of 4 14 ka 14c ages and yield nd values for north east atlantic deep water neadw that are higher 9 3 to 11 1 than those observed in the nearby madeira abyssal plain from the same depth 12 4 0 9 this observation suggests that either the iceland scotland overflow water isow and lower deep water contributions to the formation of neadw are higher along the mid atlantic ridge than in the surrounding basins or that the relative proportion of isow was higher during this period than is observed today this study indicates that hydrothermal sediments have the potential to provide a higher resolution record of deep water nd values and hence deepwater circulation patterns in the oceans than is possible from other types of sediments

la seismicite provoquee au reservoir lg 2 the construction of a rock and earth filled dam in the la grande river valley quebec has created a hydroelectric reservoir with a maximum depth of 145 m and a volume of 61 7 109 m3 no seismicit

influence of surrounding plates on 3d subduction dynamics 1 our 3d modelling study shows that the presence of lithospheric plates around a subducting plate has a significant influence on subduction dynamics in particular on trench retreat rate slab dip and lateral shortening of the subducting plate neighbouring plates prevent unrealistic plate behaviour with no need for complex rheologies because at the earth s surface plates form a continuous shell they should not be neglected

role of faults nonlinear rheology and viscosity structure in generating plates from instantaneous mantle flow models concentrated strain within plate margins and a significant toroidal component in global plate motion are among the most fundamental features of plate tectonics a significant proportion of strain in plate margins is accommodated through motion on major tectonic faults the decoupling influence of faulted plate margins primarily results from history dependent lithospheric deformation rather than from instantaneous stress weakening rheologies for instantaneous mantle flow models we argue that faults should be treated as preexisting mechanical structures with models incorporating preexisting faults a power law rheology with an exponent of 3 and slab pull and ridge push forces we demonstrate that nonlinear interaction between weak faults and this power law rheology produces plate like motion our models show that in order to produce plate like motion the frictional stress on faults needs to be small and the asthenosphere viscosity should be much weaker than that of lithosphere while both plateness and the ratio of toroidal to poloidal velocities are reduced with increasing fault coupling the viscosity contrast between the lithosphere and asthenosphere only influences plateness this shows that both diagnostics plateness and the ratio of toroidal to poloidal velocities are necessary to characterize plate motion the models demonstrate that weak transform faults can guide plate motion this guiding property of transform faults and the decoupling of thrust faults result in oblique subduction where the strike of subducted slabs is oblique to transform faults subducted slabs beneath a dipping fault produce oceanic trench and fore bulge topography and principal stresses consistent with subduction zone observations

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

shear wave splitting from local events beneath the ryukyu arc trench parallel anisotropy in the mantle wedge we present shear wave splitting measurements from local slab earthquakes at eight seismic stations of the japanese f net array located in the ryukyu arc we obtained high quality splitting measurements for 70 event station pairs and found that the majority of the measured fast directions were parallel to the strike of the trench and perpendicular to the convergence direction splitting times for individual measurements ranged from 0 25 to 2 s most values were between 0 75 and 1 25 s both the fast directions and the split times were similar to results for teleseismic s k ks and s wave splitting at the same stations which suggests that the anisotropy is located in the mantle wedge above the slab we considered several mantle deformation scenarios that would result in predominantly trench parallel fast directions and concluded that for the ryukyu subduction system the most likely explanation for the observations is corner flow in the mantle wedge combined with b type olivine fabric in this model the flow direction in the wedge is perpendicular to the trench but the fast axes of olivine crystals tend to align perpendicular to the flow direction resulting in trench parallel shear wave splitting

trench parallel flow in the southern ryukyu subduction system effects of progressive rifting of the overriding plate 1 the okinawa trough in the ryukyu subduction system is one of a few actively extending back arc basins within the continental lithosphere recent shear wave splitting measurements show variable fast directions along the trough suggesting a complex three dimensional mantle flow field in this study we use numerical subduction models to demonstrate that the thickness variations of the continental lithosphere bounding the edge of a subduction zone can result in complicated mantle circulation and regional dynamics the model results for the southern ryukyu show a combination of two effects along the okinawa trough the increasing thickness of the overriding eurasian plate toward taiwan due to a gradually diminishing rifting induces pressure gradients that drive trench parallel flow to the edge in the shallow portion of the mantle wedge and the thick eurasian lithosphere to the west effectively blocks and diverts this along arc flow and limits the toroidal flow around the slab edge below it the toroidal flow thus enters the mantle wedge at depths of more than about 100 km opposite in direction with and largely below the along arc flow these combined geometry effects of the eurasian lithosphere create an intricate three dimensional flow structure at the southern edge of the ryukyu subduction zone model predictions for lattice preferred orientations of olivine aggregates show rotation patterns that agree with the observed shear wave splitting patterns this three dimensional scenario echoes the geochemical and seismological evidence worldwide that indicates complex depth varying mantle circulations in subduction systems

timescales for the evolution of seismic anisotropy in mantle flow 1 we study systematically the relationship between olivine lattice preferred orientation and the mantle flow field that produces it using the plastic flow recrystallization model of kaminski and ribe 2001 in this model a polycrystal responds to an imposed deformation rate tensor by simultaneous intracrystalline slip and dynamic recrystallization by nucleation and grain boundary migration numerical solutions for the mean orientation of the a axes of an initially isotropic aggregate deformed uniformly with a characteristic strain rate show that the lattice preferred orientation evolves in three stages 1 for small times recrystallization is not yet active and the average a axis follows the long axis of the finite strain ellipsoid 2 for intermediate times the fabric is controlled by grain boundary migration and the average a axis rotates toward the orientation corresponding to the maximum resolved shear stress on the softest slip system 3 for the fabric is controlled by plastic deformation and average a axis rotates toward the orientation of the long axis of the finite strain ellipsoid corresponding to an infinite deformation the infinite strain axis in more realistic nonuniform flows lattice preferred orientation evolution depends on a dimensionless grain orientation lag parameter x defined locally as the ratio of the intrinsic lattice preferred orientation adjustment timescale to the timescale for changes of the infinite strain axis along path lines in the flow explicit numerical calculation of the lattice preferred orientation evolution in simple fluid dynamical models for ridges and for plume ridge interaction shows that the average a axis aligns with the flow direction only in those parts of the flow field where 1 calculation of provides a simple way to evaluate the likely distribution of lattice preferred orientation in a candidate flow field at low numerical cost

arc parallel flow within the mantle wedge evidence from the accreted talkeetna arc south central alaska 1 residual mantle exposures in the accreted talkeetna arc alaska provide the first rock analog for the arc parallel flow that is inferred from seismic anisotropy at several modern arcs the peridotites exposed at the base of the jurassic talkeetna arc have a moho parallel foliation and indicate dislocation creep of olivine at temperatures of 1000 to 1100 c slip occurred chiefly on the 001 100 slip system which has only rarely been observed to be the dominant slip system in olivine stretching lineations and olivine 100 slip directions are subparallel to the long axis of the talkeetna arc for over 200 km indicating that mantle flow was parallel to the arc axis the alignment of the olivine 100 axes yields a calculated s wave anisotropy with the fast polarization direction parallel to the arc thus 1 the fast polarization directions observed parallel to some modern arcs now have an exposed geological analog 2 arc parallel fast polarization directions can be caused by anisotropic peridotites and do not require the presence of fracture zones fluid filled pockets or glide on the 010 001 h2o induced slip system 3 seismic anisotropy beneath modern arcs may be due to slip on 001 100 with a horizontal foliation rather than slip on 010 100 with a vertical foliation and 4 the observed dominance of the 001 100 slip system may be due to high h2o concentrations suggesting that strain in the oceanic upper mantle may be accommodated dominantly by 001 100 olivine slip

tomography of the westernmost ryukyu subduction zone and the serpentinization of the fore arc mantle 1 three dimensional tomographic images of the subducting philippine sea slab and the fore arc mantle beneath ne taiwan and the westernmost ryukyu were generated in this study more than 5600 events recorded simultaneously by seismic networks in taiwan and japan were relocated for the inversion for variations in vp vs and vp vs analysis of the trade off between the data variance reduction and model variance helps to determine the appropriate strictness for regularization to avoid either overinterpretatioin or underinterpretation of data the regularization parameters were also chosen to ensure suppression of artificial vp vs anomalies the subducting slab is characterized by high vp high vs and intermediate to low vp vs notable in the mantle wedge is the high vp vs anomalies that abut the surface of the subducting slab at depths of 30 80 km the previously identified positive vp vs channel connecting the slab and the arc volcano interpreted to be a melt pathway is not reproduced in this study we convert vs and vp vs to perturbations of temperature and serpentinization of the mantle the slab is cooled by 200 c 400 c relative to the mantle in accord with the estimates from theoretical modeling of subduction the serpentinization reaches 15 or 2 water content at 50 km depth in the fore arc mantle we interpret the peak serpentinization as hydrated by the water released from the basalt eclogite metamorphic reaction in the oceanic crust of the subducting philippine sea plate the spatial limitation of the present seismic networks in this region with respect to subduction zone events hinders a full description of the pattern of melting in much of the mantle wedge resolution tests of the tomographic inversion provide a basis to guide our interpretation to better resolved regions

dynamic interaction of cold anomalies with the mid ocean ridge flow field and its implications for the australian antarctic discordance abstract negative thermal anomalies beneath a mid ocean ridge are dynamically isolated from the ambient upwelling and diverging flow field in the asthenosphere whose viscosity is on the order of 5 10 19 pa s or less this study examines on what condition a near ridge cold anomaly ascends with the upwelling passive flow and spread off axis three dimensional numerical modeling demonstrates that for a given magnitude of the cold anomaly the viscosity of the asthenosphere the spreading rate and the interference from continental rifting are the dominant controlling factors to the ascent descent of the anomaly to overcome the weight of such an anomaly and couple it with the upwelling either the spreading rate or the asthenospheric viscosity has to be high in a low viscosity asthenosphere the cold anomaly also ascends during the early stage of continental rifting due to the enhanced upwelling induced by the thick continental lithosphere the dynamic interaction between the cold anomaly and the ambient flow renders a transient nature of the subsidence of the seafloor which may lead to exaggerated temperature variation estimated by using a conduction model alone the scenarios examined are employed to place a constraint on dynamic models recently proposed for the australian antarctic discordance in which the source of the negative anomaly is hypothesized to be deeply rooted in the upper mantle with the asthenospheric viscosity less than 10 20 pa s the upwelling of the cooler material from great depths which causes a significant topographic low at the discordance is made possible only by rifting of the australian continent off the gondwanaland

stress strain and b type olivine fabric in the fore arc mantle sensitivity tests using high resolution steady state subduction zone models 1 shear wave splitting observations at several subduction systems show trench parallel fast directions in the fore arc mantle the presence of b type olivine fabric in the mantle wedge may provide an explanation for this pattern of anisotropy under low temperature and hydrated conditions sensitivity tests are shown that provide insights into the distribution and magnitude of b type fabric using two dimensional high resolution kinematic dynamic subduction zone models these models include a wet olivine rheology and a simulated aseismic creep zone in the fore arc mantle with parameterized viscous coupling between the slab and mantle wedge the calculated thermal structure is in reasonable agreement with heat flow observations and q tomography for small amounts of coupling and moderate amounts of shear heating in the seismogenic and aseismic creep zones these small amounts of viscous coupling give rise to a slowly flowing fore arc mantle with thermal stress and strain conditions suitable for significant b type fabric development the minimum shear stress necessary to drive geologically significant flow in the predicted b type region is sensitive to the magnitude of shear heating along the aseismic creep zone and ranges from less than 100 to 300 mpa the minimum time required to generate sufficient finite strain for fabric development in the b type region is around 10 myr depending on the amount of coupling and shear heating the predicted distribution of olivine fabric is consistent with a pattern of shear wave anisotropy with trench parallel fast polarization directions above the fore arc and a rapid rotation to trench perpendicular in the arc and back arc

an integrated global model of present day plate motions and plate boundary deformation summary in this paper we present a global model gsrm 1 of both horizontal velocities on the earth s surface and horizontal strain rates for almost all deforming plate boundary zones a model strain rate field is obtained jointly with a global velocity field in the process of solving for a global velocity gradient tensor field in our model we perform a least squares fit between model velocities and observed geodetic velocities as well as between model strain rates and observed geological strain rates model velocities and strain rates are interpolated over a spherical earth using bi cubic bessel splines we include 3000 geodetic velocities from 50 different mostly published studies geological strain rates are obtained for central asia only and they are inferred from quaternary fault slip rates for all areas where no geological information is included a priori constraints are placed on the style and direction but not magnitude of the model strain rate field these constraints are taken from a seismic strain rate field inferred from normalized focal mechanisms of shallow earthquakes we present a global solution of the second invariant of the model strain rate field as well as strain rate solutions for a few selected plate boundary zones generally the strain rate tensor field is consistent with geological and seismological data with the assumption of plate rigidity for all areas other than the plate boundary zones we also present relative angular velocities for most plate pairs we find that in general there is a good agreement between the present day plate motions we obtain and longterm plate motions but a few significant differences exist the rotation rates for the indian arabian and nubian plates relative to eurasia are 30 13 and 50 per cent slower than the nuvel1a estimate respectively and the rotation rate for the nazca plate relative to south america is 17 per cent slower on the other hand caribbean north america motion is 76 per cent faster than the nuvel 1a estimate while crustal blocks in the india eurasia collision zone move significantly and self consistently with respect to bounding plates only a very small motion is predicted between the nubian and somalian plates by integrating plate boundary zone deformation with the traditional modelling of angular velocities of rigid plates we have obtained a model that has already been proven valuable in for instance redefining a no net rotation model of surface motions and by confirming a global correlation between seismicity rates and tectonic moment rates along subduction zones and within zones of continental deformation

reservoir triggered seismicity at the toulnustouc hydroelectric project quebec north shore canada hydro quebec s toulnustouc hydroelectric complex is located approximately 120 km north of baie comeau on the quebec north shore region the project involved the construction of a 77 m high dam and a 45 m high dike that allowed the expansion of an existing 235 km2 reservoir by some 22 km2 when completely filled the new part of the reservoir has a maximum depth of 75 m starting on 10 february 2005 and over a period of four days the water of the original reservoir was used to start filling the new part the water level of the entire reservoir increased by another 10 m after the spring snowmelt the project also included a 9 8 km long headrace tunnel constructed in precambrian bedrock that creates a 170 m water head between the reservoir and the powerhouse electricity production started in may 2005 the complex is located in an area of low level earthquake activity on the canadian shield four stations of the canada national seismograph network cnsn located within 100 km of the reservoir provide a location threshold with a magnitude of about 1 2 m n to complement the cnsn two portable seismographs were installed near the flooded area and the headrace tunnel about one week before flooding started between february and mid may 2005 a total of four triggered earthquakes were recorded by the field instruments and the cnsn stations the first event a magnitude 1 4 m n occurred 16 days after flooding of the extension of the reservoir started three events were located beneath the reservoir and one was near the tunnel between the reservoir and the powerhouse where the water pressure is most intense numerous smaller but unlocatable events m n 1 0 were detected by one of the field stations about two years later two events separated by one minute and recognized as multiplets occurred near the reservoir the small stress perturbations that led to the reservoir triggered seismicity at toulnustouc support the hypothesis that the mid to upper crust of intraplate areas is in a state of stress close to failure

mantle flow in the south sandwich subduction environment from source side shear wave splitting 1 the occurrence of seismic anisotropy in the earth s upper mantle is a global phenomenon related to subcrustal deformation and flow processes the shear wave splitting analysis method has led to a global set of anisotropy maps mainly derived from receiver side analysis remote places with few seismometers deployed remain unexamined source side splitting analysis allows mapping of mantle fabrics in these regions here we investigate seismic anisotropy in the south sandwich islands subduction environment core reflected scs waves recorded at the neumayer seismograph network are corrected for well constrained receiver anisotropy and then analysed for source anisotropy sub slab mantle minerals are aligned horizontally almost parallel to the trench indicating a westward flow around the subducting slab this is consistent with a model of horizontal mantle flow due to slab rollback that was previously inferred from marine and geochemical studies in the back arc region

an eulerian technique for thermomechanical modeling of lithospheric extension a two dimensional eulerian spline function finite element method for modeling the coupled thermomechanical processes which control lithospheric extension is presented the model domain comprises a cross section through the lithosphere and part of the asthenosphere material may enter or leave through the bottom and side boundaries on the side boundary a fixed velocity or a fixed force condition can be applied highly nonlinear and thermally activated creep laws are used for the brittle and ductile rheologies of the lithosphere compositional differences between crust and mantle are described by a continuous function which is advected and which affects density rheology and heat production an approximate topography is calculated from the normal stress at a free slip surface boundary a set of sample case studies is presented which highlights the capabilities and possible applications of the method it also serves to study the influence of numerical parameters and of conceptual assumptions in various nontrivial applications the examples comprise narrow rifting continuing into the stage of seafloor spreading formation of a wide basin formation of a metamorphic core complex and impact of a mantle plume onto the bottom of the lithosphere in the last case the entire upper mantle is modeled and the instability of a thermal bottom boundary layer gives rise to a diapiric mantle plume the main weakness of the eulerian approach as compared to the lagrangian technique was found to be the effective blurring of material discontinuities like the moho this is of little consequence in pure shear deformation but when strong simple shear occurs it impairs the resolution the main advantage of the eulerian technique is that very large deformation poses no problem as it does for the lagrangian method the new approach is especially suited to study processes where strong differential flow occurs in part of the modeled system a further possible advantage of the method is its efficiency in terms of computer time

patterns in seismic anisotropy driven by rollback subduction beneath the high lava plains received 24 march 2011 revised 25 may 2011 accepted 31 may 2011 published 8 july 2011 1 we present three dimensional laboratory modeling of the evolution of finite strain and compare these to shear wave splitting observations in the northwest u s under the high lava plains hlp we show that relationships between mantle flow and anisotropy are complicated in subduction zones and factors such as initial orientation of the olivine fast axis style of subduction and time evolving flow are important due to increased horizontal shear systems with a component of rollback subduction have simple trench normal strain alignment within the central region of the backarc mantle wedge while those with more simple longitudinal sinking are often variable and complex in the hlp splitting observations are consistent with rollback driven laboratory results citation druken k a m d long and c kincaid 2011 patterns in seismic anisotropy driven by rollback subduction beneath the high lava plains geophys res lett 38 l13310 doi 10 1029 2011gl047541

okinawa trough backarc basin early tectonic and magmatic evolution the okinawa trough lying between japan and taiwan is a backarc basin formed by extension within the continental lithosphere behind the ryukyu trench arc system stress directions associated with the two last extensional phases in the southwestern okinawa trough have been deduced from a comparison with analog modeling the direction of extension is n 150 o for the pleistocene phase of extension 2 0 1 ma and n 170 o for the late pleistocene holocene phase of extension 0 1 0 ma the present day ryukyu volcanic arc a narrow continuous feature extending from japan to taiwan is located on the eastern side of the okinawa trough 80 100 km above the wadati benioff zone the minimum depth for emplacement of arc magmatism scarce present day backarc volcanism appears in the middle and southern okinawa trough within linear en echelon bathymetric depressions a n045 o oriented seamount volcanic chain cuts across obliquely the southwestern okinawa trough and lies in the direct line of extension of the gagua ridge a n s linear volcanic feature of the philippine sea plate associated with this extension of the gagua ridge a large reentrant located at the base of the ryukyu prism the uplift of part of the nanao forearc basin and the deformation of the sedimentary arc suggest that the voluminous cross backarc volcanism could be tied to the subduction of the gagua ridge located there at a depth of 80 100 km beneath the backarc basin a second area of anomalous volcanism has been identified in the middle okinawa trough in the ene extension of the daito ridge a wnw ese 400 km long volcanic feature of the philippine sea plate we suggest that the gagua and daito ridges initially induced stress at the base of the arc which is still brittle and cracks propagated through the overlying brittle lithosphere allowing magmas with arc affinities to erupt at the seafloor this excessive magmatism reaches the seafloor through conduits which preferentially follow in their shallowest portion the crustal normal faults of the backarc rifts the okinawa trough is consequently still in an early stage of evolving from arc type to backarc activity

subduction dynamics and the origin of andean orogeny and the bolivian orocline the building of the andes results from the subduction of the oceanic nazca plate underneath the south american continent however how and why the andes and their curvature the bolivian orocline formed in the cenozoic era 65 5 million years myr ago to present despite subduction continuing since the mesozoic era 251 0 65 5 myr ago is still unknown three dimensional numerical subduction models demonstrate that variations in slab thickness arising from the nazca plate s age at the trench produce a cordilleran morphology consistent with that observed the age dependent sinking of the slab in the mantle drives traction towards the trench at the base of the upper plate causing it to thicken thus subducting older nazca plate below the central andes can explain the locally thickened crust and higher elevations here we demonstrate that resultant thickening of the south american plate modifies both shear force gradients and migration rates along the trench to produce a concave margin that matches the bolivian orocline additionally the varying forcing along the margin allows stress belts to form in the upper plate interior explaining the widening of the central andes and the different tectonic styles found on their margins the eastern and western cordilleras the rise of the central andes and orocline formation are directly related to the local increase of nazca plate age and an age distribution along the margin similar to that found today the onset of these conditions only occurred in the eocene epoch this may explain the enigmatic delay of the andean orogeny that is the formation of the modern andes

sub slab mantle flow parallel to the caribbean plate boundaries inferences from sks splitting abstract upper mantle deformation near the margins of the caribbean plate is investigated using observations of shear wave splitting in teleseismic and local shear phases the caribbean plate is almost stationary in the hot spot reference frame and is wedged between the north america south america nazca and cocos plates collisional belts and major shear zones encircle the plate data from seismic stations operated by iris geoscope the venezuelan seismological network and the british geological survey have provided nearly 2000 seismic records for analysis analysis of shear wave splitting in teleseismic core phases e g sks at stations reveals fast shear wave polarisations that are conformal to the plate boundary paralleling the major structural features the magnitude of the splitting is in general quite large 1 2 2 1 s in northern venezuela the magnitude of splitting increases towards the caribbean south american collisional front local shear phases from earthquakes up to 200 km deep beneath northeastern venezuela show very small amounts of splitting 0 1 0 3 s analysis of the depth dependence in the magnitude of the splitting suggests that most of the upper mantle wedge is isotropic and the splitting in the local phases is mostly accrued in the crust and uppermost mantle in ne venezuela fast shear wave polarisations in both the local and teleseismic phases closely parallel the fault systems in the region suggesting that the crust and mantle are coupled in this area stations on the island of montserrat and near bucaramanga columbia also show very small magnitudes in local s wave splitting but large amounts of sks splitting in general the small magnitude of the local s phase splitting suggests that the teleseismic phases accrue considerable splitting beneath the top of the slab we interpret the bulk of the splitting in terms of sub slab flow that is forced around the caribbean plate due to the convergence of the surrounding plates

influence of trench width on subduction hinge retreat rates in 3 d models of slab rollback subduction of tectonic plates limited in lateral extent and with a free trailing tail i e free subduction is modeled in a three dimensional 3 d geometry the models use a nonlinear viscoplastic rheology for the subducting plate and exhibit a wide range of behaviors depending on such plate characteristics as strength width and thickness we investigate the time evolution of this progressive rollback subduction measure the accompanying return flow in the upper mantle and quantify the plate kinematics due to the 3 d geometry flow is allowed to accompany slab rollback around the lateral edges of the slab the toroidal component as opposed to 2 d geometry where material is forced to flow underneath the slab tip the poloidal component a simple force balance is provided which relates the speed of backward trench migration to the resistive forces of generating flow and weakening the plate our results indicate most of the gravitational energy of the system i e the negative buoyancy of the subducting slab is converted into a toroidal flow 69 a much smaller amount goes into weakening the plate 18 and the remaining amount goes into driving flow parallel to displacement of the slab 13 for the trench widths w we investigate 1500 km a maximum trench retreat rate occurs for trenches 600 km wide which is attributed to the interaction between a plate of finite width and the induced flow which has a lengthscale in the horizontal direction these numerical results quantitatively agree with comparable 3 d laboratory experiments using analogue models with a purely viscous plate material schellart 2004a 2004b including correlations between increasing retreat rate with increasing plate thickness trench width for maximum retreat rate 500 km and estimated amount of slab buoyancy used to drive rollback induced flow 70 several implications for plate tectonics on earth result from these models such as rollback subduction providing a physical mechanism for ephemeral slab graveyards situated above the more viscous lower mantle and endothermic phase transition prior to a flushing event into the lower mantle mantle avalanche

slab stiffness control of trench motion insights from numerical models 1 subduction zones are not static features but trenches retreat roll back or advance here we investigate the dominant dynamic controls on trench migration by means of two and three dimensional numerical modeling of subduction this investigation has been carried out by systematically varying the geometrical and rheological model parameters our viscoplastic models illustrate that advancing style subduction is promoted by a thick plate a large viscosity ratio between plate and mantle and a small density contrast between plate and mantle or an intermediate width w 1300 km advancing slabs dissipate 45 to 50 of the energy in the system thin plates with relatively low viscosity or relatively high density or wide slabs w 2300 km on the other hand promote subduction in the retreating style i e slab roll back the energy dissipated by a retreating slab is 35 to 40 of the total dissipated energy most of the energy dissipation occurs in the mantle to accommodate the slab motion whereas the lithosphere dissipates the remaining part to bend and unbend with a simple scaling law we illustrate that this complex combination of model parameters influencing trench migration can be reduced to a single one plate stiffness stiffer slabs cause the trench to advance whereas more flexible slabs lead to trench retreat the reason for this is that all slabs will bend into the subduction zone because of their low plastic strength near the surface but stiff slabs have more difficulty unbending at depth when arriving at the 660 km discontinuity those bent slabs tend to cause the trench to advance in a similar way variation of the viscoplasticity parameters in the plate may change the style of subduction a low value of friction coefficient weakens the plate and results in a retreating style while higher values strengthen the plate and promote the advancing subduction style given the fact that also on earth the oldest and therefore probably stiffest plates have the fastest advancing trenches we hypothesize that the ability of slabs to unbend after subduction forms the dominant control on trench migration

global mantle flow and the development of seismic anisotropy differences between the oceanic and continental upper mantle 1 viscous shear in the asthenosphere accommodates relative motion between earth s surface plates and underlying mantle generating lattice preferred orientation lpo in olivine aggregates and a seismically anisotropic fabric because this fabric develops with the evolving mantle flow field observations of seismic anisotropy can constrain asthenospheric flow patterns if the contribution of fossil lithospheric anisotropy is small we use global viscous mantle flow models to characterize the relationship between asthenospheric deformation and lpo and compare the predicted pattern of anisotropy to a global compilation of observed shear wave splitting measurements for asthenosphere 500 km from plate boundaries simple shear rotates the lpo toward the infinite strain axis isa the lpo after infinite deformation faster than the isa changes along flow lines thus we expect the isa to approximate lpo throughout most of the asthenosphere greatly simplifying lpo predictions because strain integration along flow lines is unnecessary approximating lpo with the isa and assuming a type fabric olivine a axis parallel to isa we find that mantle flow driven by both plate motions and mantle density heterogeneity successfully predicts oceanic anisotropy average misfit 13 continental anisotropy is less well fit average misfit 41 but lateral variations in lithospheric thickness improve the fit in some continental areas this suggests that asthenospheric anisotropy contributes to shear wave splitting for both continents and oceans but is overlain by a stronger layer of lithospheric anisotropy for continents the contribution of the oceanic lithosphere is likely smaller because it is thinner younger and less deformed than its continental counterpart

reconciling strong slab pull and weak plate bending the plate motion constraint on the strength of mantle slabs although subducting slabs undergo a bending deformation that resists tectonic plate motions the magnitude of this resistance is not known because of poor constraints on slab strength however because slab bending slows the relatively rapid motions of oceanic plates observed plate motions constrain the importance of bending we estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity predicting plate motions using a global mantle flow model we constrain the viscosity of the bending slab to be at most similar to 300 times more viscous than the upper mantle stronger slabs are intolerably slowed by the bending deformation weaker slabs however cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of similar to 10 15 s 1 the constrained bending viscosity similar to 2 x 10 23 pa s is larger than previous estimates that yielded similar or larger bending resistance here similar to 25 of forces this apparent discrepancy occurs because slabs bend more gently than previously thought with an average radius of curvature of 390 km that permits subduction of strong slabs this gentle bending may ultimately permit plate tectonics on earth c 2008 elsevier b v all rights reserved

the influence of plate motions on three dimensional back arc mantle flow and shear wave splitting both the polarization direction of the fast shear waves and the types of deformation within overriding plates vary between the back arc basins of western pacific subduction zones the goal of this study is to test the possibility that motions of the overriding plates may control the patterns of seismic anisotropy and therefore the observed shear wave splitting we calculated three dimensional models of viscous asthenospheric flow driven by the motions of the subducting slab and overriding plates shear wave splitting was calculated for polymineralic anisotropy within the back arc mantle wedge assuming that the anisotropy was created by flow induced strain predicted splitting may differ substantially depending on whether anisotropy is computed directly using polycrystalline plasticity models or is based on the orientation of finite strain a parameter study shows that both finite strain and textural anisotropy developed within three dimensional plate coupled asthenospheric flow models are very heterogeneous when back arc shearing occurs within the overriding plate therefore predicted shear wave splitting varies strongly as a function of plate motion boundary conditions and with ray path through the back arc asthenosphere flow models incorporating plate motion boundary conditions for the tonga southern kuril and eastern aleutian subduction zones produce splitting parameters consistent with observations from each region trench parallel flow generated by small variations in the dip of the subducting plate may be important in explaining observed fast directions of anisotropy sampled within the innermost corner of the mantle wedge

seismic anisotropy in the upper mantle 2 predictions for current plate boundary flow models 1 the anisotropic seismic structure due to flow induced mineral alignment is investigated for a series of models designed to simulate deformation in the upper mantle within a few hundred kilometers of a plate boundary the orientation distributions of olivine enstatite aggregates evolve along streamlines of each flow model based on each grains plastic response to the local stress strain field the effective elastic tensor for these textured aggregates provides predictions of p wave anisotropy and shear wave splitting throughout the model space p and s travel time delay patterns and fast shear wave polarization angles are found to vary significantly with incidence angle for a given model comparison of predicted fast p direction for our method versus a finite strain based estimate shows that agreement is acceptable for much of the model space but notable differences occur in regions up to several tens of kilometers in size two dimensional models of spreading center flow are presented for slow and fast rates and for several cases in which the ridge migrates over the deeper mantle the effect of flow in the third dimension is addressed in a few calculations for one comparison of flow in the mantle wedge at a subduction zone the introduction of trench parallel flow causes significant changes in the predicted patterns of p wave anisotropy magnitude more than orientation and sks splitting

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

prediction of anisotropy from flow models a comparison of three methods observations of anisotropy in earth are used regularly as constraints for models of deformation using various assumptions about the relationship between deformation and the resulting anisotropic fabric we compare three methods for calculating fabric from velocity fields tracking of finite strain ellipses a kinematic crystallographic code and the evolution of directors we find that the use of the finite strain ellipse provides only limited prediction capabilities as it cannot reproduce experimental observations that involve recrystallization the crystallographic code we tested a variant of the popular code d rex provides a more complete fabric prediction but at a much higher computational cost the directors method provides an intermediate solution while it does not include some of the more complex crystallographic processes that d rex does the results of this method closely resemble those of d rex at a lower computational cost the directors are also capable of tracking anisotropy at much larger strains than d rex we conclude that when computation speed is important for example in self consistent geodynamic models with anisotropic rheology the directors method provides an appropriate approximation

seafloor spreading in the weddell sea and southwest atlantic since the late cretaceous geosat radar altimetry defines a pattern of closely spaced curvilinear gravity anomalies in the northern weddell sea and adjacent parts of the southwest atlantic these anomalies are caused by a fracture zone pattern resulting from seafloor spreading mainly on the southern flank of a ridge system which now exists only to the east of the south sandwich trench flowlines derived from these anomalies are combined with isochrons resulting from an interpretation of all available magnetic anomaly profiles in the region to produce a new tectonic summary chart this chart reveals major changes in spreading direction during i the cretaceous magnetic superchron ii the latest cretaceous paleocene 65 ma and iii the mid eocene 50 ma the latter two changes enclose a period of very slow spreading and define a kink in the flowlines which is discernable throughout most of the weddell sea preliminary modelling indicates that trends from the western weddell sea to near the south atlantic triple junction can be explained by the separation of just two plates there is some disagreement between the observed flowlines and those predicted by rotations published for the south atlantic and southwest indian ridge but abrupt changes in spreading direction and a slowing of rates are observed on all three branches of the south america africa antarctica plate system between the times of anomalies c31 and c21 this suggests that the lack of closure results from errors in the calculated rotations and that the seafloor of the weddell sea and southwest atlantic was created as a direct consequence of south america antarctica plate motion the data do not support the existence of a late cretaceous malvinas plate within the weddell sea although microplates in the northernmost weddell sea do appear to have formed briefly as a result of interactions between the former spreading ridge and a subduction zone located at the south scotia ridge

a laboratory model of subduction zone anisotropy abstract measurements of split seismic shear waves near pacific subduction zones show that fast polarization directions are primarily trench parallel in the seaward side mantle beneath a subducting slab and variable in the mantle wedge beneath back arc basins it is generally accepted that the azimuth of the fast polarization direction defines the average orientation of the crystallographic a axis in olivine here we use laboratory experiments to determine the orientation of the olivine a axis in the upper mantle in the vicinity of a subducting slab a laboratory analog for olivine a axis motion during creep deformation is developed using small cylinders whiskers suspended in a viscous fluid subduction beneath a stationary overriding plate variable dip angles and realistic rollback and down dip plate kinematics are examined plate rollback produces strong plate parallel whisker alignments for all dip angles examined alignment in the wedge occurs in three layers a layer of dip parallel alignment above the downgoing plate a layer of subhorizontal alignment beneath the overriding plate and a core of intermediate orientations the amount of alignment beneath the overriding plate and the wedge core decreases with increasing dip angle our results indicate that olivine a axis orientation in the seaward side mantle is controlled by the amount of slab rollback and that orientations in the mantle wedge depend upon slab dip angle

kinematics of subduction and subduction induced flow in the upper mantle 1 results of fluid dynamical experiments are presented to model the kinematics of lithospheric subduction in the upper mantle the experiments model a dense high viscosity plate subducting lithosphere overlying a less dense low viscosity layer upper mantle the overriding lithosphere is not incorporated several important features of slab behavior were investigated including the temporal variability of hinge line migration the kinematic behavior of the slab and the subduction induced upper mantle flow both fixed and free trailing edge boundary conditions of the subducting plate were investigated results show that hinge line retreat is a natural consequence of subduction of a negatively buoyant slab the migration rate increases until the slab approaches the upper lower mantle discontinuity resulting in a decrease in migration rate followed by a renewed increase and finally approaching a steady state slab retreat results in mantle flow with material initially located underneath the slab flowing around the lateral slab edges toward the mantle wedge experimental results indicate that all rollback induced flow occurs around the lateral slab edges forcing the hinge line to attain a convex shape toward the direction of retreat no signs for poloidal flow underneath the slab tip have been detected only a small component of toroidal type flow was observed underneath slanting slab tips for a fixed trailing edge the slab does not sink vertically downward but sinks at an angle in a regressive manner for a free trailing edge slab sinking is oriented more vertically while the surface part of the subducting plate is pulled into the subduction zone

the southwestern edge of the ryukyu subduction zone a high q mantle wedge abstract the lateral edge of a subduction zone is usually depicted as an opening to the asthenosphere where invigorated dynamics and amplified magmatism take place in this study we present evidence from seismic data suggesting the presence of a cold and dynamically sluggish edge environment at the southwest end of the ryukyu subduction system we measured attenuation or 1 q for p waves from subduction zone events at 100 km depths received by obss in the okinawa trough and land stations in ne taiwan in the okinawa trough 100 200 km from the edge q values are lower than 100 in the vicinity of the edge q values increase from 100 to over 1000 towards taiwan to reconcile arguments from geophysical and geochemical observations we propose that the mantle wedge near the edge has high q values due to low temperatures and probably low water content these may result from coupling of the slab laterally with the thick eurasian lithosphere which inhibits back arc rifting retards subduction and reduces the water supply to the mantle wedge the sw ryukyu subduction system represents a subduction zone edge type distinct from more commonly documented free or warm edges

a model for the termination of the ryukyu subduction zone against taiwan a junction of collision subduction separation and subduction boundaries 1 the nw moving philippine sea plate psp collides with the eurasian plate eup in the vicinity of taiwan and at the same time it subducts toward the north along sw ryukyu the ryukyu subduction zone terminates against eastern taiwan while the ryukyu trench is a linear bathymetric low about 100 km east of taiwan closer to taiwan it cannot be clearly identified bathymetrically owing to the deformation related to the collision making the location of the intersection of the ryukyu with taiwan difficult to decipher we propose a model for this complex of boundaries on the basis of seismicity and 3 d velocity structures in this model the intersection is placed at the latitude of about 23 7 n placing the northern part of the coastal range on eup as psp gets deeper along the subduction zone it collides with eup on the taiwan side only where they are in direct contact thus the eurasian plate on the taiwan side is being pushed and compressed by the nw moving philippine sea plate at increasing depth toward the north offshore of northeastern taiwan the wedge shaped eup on top of the ryukyu subducting plate is connected to the eup on the ryukyu side and coupled to the nw moving psp by friction at the plate interface the two sides of the eup above the western end of the subduction zone are not subjected to the same forces and a difference in motions can be expected the deformation of taiwan as revealed by continuous gps measurements geodetic movement along the east coast of taiwan and the formation of the hoping basin can be understood in terms of the proposed model

scotia sea regional tectonic evolution implications for mantle flow and palaeocirculation abstract the scotia sea and surrounding scotia arc have evolved over the past 40 ma by extension behind an east migrating subduction zone at the boundary between the south american sam and antarctic ant plates the considerable data set now available regional geology and geophysics earthquake seismology satellite altimetry global plate analyses suggest why east migrating subduction began what has been the driving force that has sustained it and what other processes have controlled the mode of back arc extension in the scotia sea a suite of six reconstructions has been developed based on this data set the reconstruction to 40 ma creates a compact cuspate continental connection between south america and the antarctic peninsula at the subducting pacific margin with fragments now dispersed around the scotia arc occupying positions within it compatible with their known geology the driving force has been subduction of south american ocean floor which began as a result of southward migration of the pole of south american antarctic plate rotation and a key modulator of back arc extension has been collision of ridge crest sections of the south american antarctic plate boundary with the east advancing trench cenozoic regional tectonic evolution has two other likely consequences which greatly increase its importance firstly this region saw the tectonic disruption of the final barrier to complete circum antarctic deep water flow that may have had a profound effect on palaeoclimate secondly it is possible that the rapid roll back of the hinge of subduction is related to shallow eastward flow in the sub lithospheric mantle both of these consequences are explored the reconstructions show that rapid roll back of the subduction hinge averaging 50 mm a over the last 40 ma with respect to the south american plate has been a feature of all of scotia sea evolution and provide a history of motion of several oceanic microplates most of which are now welded together within the scotia sea this will guide the location of seismometers and or dredge hauls to test the hypothesis of shallow mantle flow and help interpret the results the reconstructions also allow an assessment of the creation of deep water pathways that would have permitted the development of the present day antarctic circumpolar current acc an early miocene onset within the period 22 17 ma seems likely for the acc depending on the structure and palaeo elevation of davis bank and aurora bank sections of the north scotia ridge however the study shows there was a delay of one or more million years between initial provision of a deep water pathway and the major mid miocene change in global climate involving the general level of antarctic glaciation that may have been related if these changes were related then the delay suggests that other factors possibly rough elevated ocean floor but also non tectonic factors such as atmospheric co2 were important in determining palaeoclimate

induced seismicity at the lg3 reservoir james bay quebec canada the filling of the second reservoir lg3 in the james bay hydroelectric development has resulted in a number of small induced earthquakes up to magnitude 3 7 the activity has occurred in two main areas one of which can be associated with the previously mapped lg3 fault and the other with the intersection of two dominant regional structures the significant induced activity started 7 months after filling commenced when the water depth had reached 64 m and continued for another 11 months after a lull of about 6 months a short burst of events much more energetic than that during filling occurred probably on the lg3 fault after a 2 m lowering of the water level a further episode of seismicity occurred in the spring of 1984 also during the lowering of the water level based on the data available it is suggested that postfilling seismicity permits one to isolate two classes of reservoirs from the general population of reservoirs those with seismicity associated with refilling and those with seismicity associated with lowering water levels lg3 belongs to the latter class where seismicity recurrs during lowering of a few meters of water

small scale convection under the back arc occurring in the low viscosity wedge abstract water released from subducting slabs through a dehydration reaction may lower the viscosity of the mantle significantly thus we may expect a low viscosity wedge lvw above the subducting slabs the lvw coupled with a large scale flow induced by the subducting slabs may allow the existence of roll like small scale convection whose axis is normal to the strike of the plate boundary such a roll structure may explain the origin of along arc variations of mantle temperature proposed recently in northeast japan we study this possibility using both 2d and 3d models with without pressure and temperature dependent viscosity 2d models without pressure and temperature dependence of viscosity show that with a reasonable geometry of the lvw and subduction speed small scale convection is likely to occur when the viscosity of the lvw is less than 10 19 pa s corresponding 3d model studies reveal that the wavelength of rolls depends on the depth of the lvw the inclusion of temperature dependent viscosity requires the existence of further low viscosity in the lvw since temperature dependence suppresses the instability of the cold thermal boundary layer pressure i e depth dependence coupled with temperature dependence of the viscosity promotes short wavelength instabilities the model which shows a relatively moderate viscosity decrease in the lvw most of the lvw viscosity is 10 18 10 19 pa s and a wavelength of roll 80 km has a rather small activation energy and volume 130 kj mol and 4 cm 3 mol of the viscosity this small activation energy and volume may be possible if we regard them as an effective viscosity of non linear rheology

kinematics and flow patterns in deep mantle and upper mantle subduction models influence of the mantle depth and slab to mantle viscosity ratio 1 three dimensional fluid dynamic laboratory simulations are presented that investigate the subduction process in two mantle models an upper mantle model and a deep mantle model and for various subducting plate mantle viscosity ratios sp m 59 1375 the models investigate the mantle flow field geometrical evolution of the slab sinking kinematics and relative contributions of subducting plate motion and trench migration to the total rate of subduction all models show that the subducting plate is always moving trenchward resulting from slab pull furthermore all deep mantle models show trench retreat as do upper mantle models in the initial stage of subduction before slab tip transition zone interaction upper mantle models with a low sp m 66 217 continue to show trench retreat after interaction upper mantle models with a high sp m 378 709 show a period of trench advance after interaction followed by trench retreat upper mantle models with a very high sp m 1375 show continued trench advance after interaction the difference in trench migration behavior and associated slab geometries is attributed to both sp m and the mantle depth to plate thickness ratio tm tsp which both affect the slab bending radius to mantle thickness ratio rb tm four subduction regimes can be defined regime i with rb tm 0 3 trench retreat slab draping and a concave trench regime ii with 0 3 rb tm 0 5 episodic trench migration slab folding and a concave trench regime iii with rb tm 0 5 trench advance slab rollover geometries and minor trench curvature and regime iv with rb tm 0 8 trench retreat slab draping and a rectilinear trench in all models slab parallel downdip motion induces poloidal mantle flow structures in addition trench retreat and rollback motion of the slab induce quasi toroidal return flow around the lateral slab edges toward the mantle wedge rollback induced poloidal flow around the slab tip is not observed in any of the experiments finally comparison between the slab geometries observed in the upper mantle models and slab geometries observed in nature imply that the effective viscosity ratio between slab and ambient upper mantle in nature is less than 103 and of the order 1 7 102 with a best estimate of 1 3 102

a rigid block rotation model for the gps derived velocity field along the ryukyu arc abstract using a gps derived velocity field we try to estimate the kinematic features of the ryukyu arc and discuss its implications to the spreading of the okinawa trough the gps derived strain rate field implies that internal deformation within the ryukyu arc is not significant therefore we adopt a rigid block rotation model for the motion of the ryukyu arc the velocity field in the ryukyu arc is best explained by a model with three blocks rather than one with only one block the estimated block rotations predict the spreading along the okinawa trough the predicted spreading is slow 10 mm per year in the north and fast 50 mm per year in the southwest the direction of the predicted spreading is normal to the okinawa trough n170 e n180 e in the southwest while in the north the direction of the spreading is oblique to the okinawa trough n150 e n180 e