temporal changes of shallow seismic velocity around the karadere duzce branch of the north anatolian fault and strong ground motion we analyze temporal variations of seismic velocity along the karadere duzce branch of the north anatolian fault using seismograms generated by repeating earthquake clusters in the aftershock zones of the 1999 mw7 4 izmit and mw7 1 duzce earthquakes the analysis employs 36 sets of highly repeating earthquakes each containing 4 18 events the events in each cluster are relocated by detailed multi step analysis and are likely to rupture approximately the same fault patch at different times the decay rates of the repeating events in individual clusters are compatible with the omori s law for the decay rate of regional aftershocks a sliding window waveform cross correlation technique is used to measure travel time differences and evolving decorrelation in waveforms generated by each set of the repeating events we find clear step like delays in the direct s and early s coda waves sharp seismic velocity reduction immediately after the duzce main shock followed by gradual logarithmic type recoveries a gradual increase of seismic velocities is also observed before the duzce main shock probably reflecting post seismic recovery from the earlier izmit main shock the temporal behavior is similar at each station for clusters at various source locations indicating that the temporal changes of material properties occur in the top most portion of the crust the effects are most prominent at stations situated in the immediate vicinity of the recently ruptured fault zones and generally decrease with normal distance from the fault a strong correlation between the co seismic delays and intensities of the strong ground motion generated by the duzce main shock implies that the radiated seismic waves produced the velocity reductions in the shallow material

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

mantle flow in subduction systems the subslab flow field and implications for mantle dynamics 1 the character of the mantle flow field in subduction zones remains poorly understood despite its importance for our understanding of subduction dynamics in particular little attention has been paid to mantle flow beneath subducting slabs in order to identify processes that make first order contributions to the global pattern of subslab mantle flow we have compiled shear wave splitting measurements from subduction zones worldwide from previously published studies and estimated average splitting parameters for the subslab region globally the subslab region is overwhelmingly dominated by trench parallel fast splitting directions we tested for relationships between splitting delay time a measure of the strength of anisotropy and parameters that are indicators of tectonic processes such as trench migration velocity convergence velocity and obliquity age of subducting lithosphere and slab dip curvature seismicity and thickness we used several different plate motion models to describe plate and trench motions and evaluated the differences among the models we find that only one parameter namely the trench migration velocity in a pacific hot spot reference frame appears to correlate well with the strength of the subslab splitting signal this finding supports a model in which the mantle beneath subducting slabs is dominated by three dimensional flow induced by trench migration we explore several implications of our model for various aspects of mantle dynamics including the choice of a suitable reference frame s for mantle flow the existence of a thin decoupling zone between slabs and the subslab mantle and consequences for mass transfer between the upper and lower mantle

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

variability in flow and temperatures within mantle subduction zones 1 a series of laboratory experiments is used to model three dimensional aspects of flow in subduction zones and the consequent temperature variations in the slab and overlying mantle wedge the effects of longitudinal rollback and slab steepening components of motions are considered along with different thicknesses of the over riding lithosphere the results show that the style of plate sinking influences the evolution of subduction zones both in terms of the speed orientation and temperature of flow in the overlying wedge and in terms of the temperatures at the surface of the descending slab in the simplest case of longitudinal sinking without rollback motion velocities in the mantle wedge are 30 40 of the downdip plate speed return flow paths in the shallow mantle wedge are nearly horizontal for slow slab speeds and steepen with higher slab speeds and there is no mass flux around the edges of a slab segment of finite width rollback subduction leads to flow both around and beneath the sinking slab with larger velocities in the wedge up to 150 of the slab speed and flow focused toward the center of the plate segment rollback subduction in which either the trench migrates or the plate steepens with time induces shallow and steep wedge return flow trajectories respectively the thermal evolution of the plate is strongly influenced by sinking style and rate highest temperatures are along the edges of the slab for longitudinal sinking but along the centerline of the slab segment for rollback motion slab surface temperatures when scaled to the mantle are higher than those given by previous two dimensional numerical models and are consistent with recent observational and experimental data on melt compositions

automatic measurement of shear wave splitting and applications to time varying anisotropy at mount ruapehu volcano new zealand 1 we present an automatic shear wave splitting measurement tool for local earthquakes with the sole manual step of choosing an s arrival time we apply the technique to three data sets recorded on mount ruapehu volcano in new zealand that have previously been determined to have fast polarizations that vary in time and with earthquake depth the technique uses an eigenvalue minimization technique applied over multiple measurement windows the dominant period of each waveform sets minimum and maximum window lengths cluster analysis determines the best solution among all the windows and quality grading criteria assess the results automatically when the same filters are used for events determined to be high quality from manual studies the automatic technique returns virtually identical results scatter increases when the automatic technique chooses the best filter but the average automatic results remain consistent with the manual results when the automatic technique is used on sets that include data previously judged as poor quality some stations yield distributions of fast polarizations that include peaks that were not present in previously published results the difference may stem from two factors automatic grading lets through some measurements that independent analysts consider poor quality but also unconscious bias in the manual selection process may downgrade measurements that do not fit expectations nonetheless the new objective analysis confirms changes in the average fast polarizations between 1994 and 2002 and between shallow and deep events therefore this new technique is valuable for objective assessment of anisotropy and its variation in time

shear wave splitting beneath southern korea and its tectonic implication abstract seismic anisotropy beneath the southern korean peninsula skp was investigated based on the scs waves of deep focused earthquakes recorded at 35 broadband stations the most striking feature of the observed splitting patterns is the preference of the fast directions nw se that are nearly parallel to the direction of the absolute plate motion in the region however the splitting patterns over the region show significant variation in splitting parameters indicating a complex anisotropic structure variations of the splitting directions and dissimilarity in the source domains of basaltic volcanisms suggest that the asthenospheric mantle flow since at least the late cenozoic cannot explain the seismic anisotropy beneath the region comparison to shear wave splitting measurements from eastern china revealed that the nw se fast direction of splitting measurements in the skp is close to that in the north china block ncb while the ne sw fast direction might be related to that in the south china block scb the shallow mantle lithosphere beneath the skp retains the fossil anisotropy amalgamated prior to the late paleozoic before the collision between the ncb and scb and the anisotropic structure was not completely realigned by the major orogenic events during the late paleozoic to mesozoic eras

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

acoustic measurements of surface gravity wave spectra in monterey bay using mode travel time fluctuations earlier work by miller et al j acoust soc am 86 326 345 1989 suggested that the surface gravity wave frequency direction spectrum could be measured using acoustic tomography signals an experiment was performed in monterey bay california from 12 16 december 1988 to validate this earlier work during that period an acoustic tomography source with a center frequency of 224 hz and bandwidth of 16 hz was deployed continuously repeating its coded transmissions every 1 9375 s the transmissions were received on a number of modified sonobuoys placed on the continental shelf at a range of approximately 55 km although the acoustic propagation was complicated by the presence of the monterey submarine canyon a hybrid ray and mode interpretation explains the observations the main effect of the surface waves on the acoustic signals is determined to occur in the shallow continental shelf portion of the path on the shelf it is shown that the resolved acoustic normal modes exhibit the fluctuations in travel time described in our earlier work the acoustic measurements of the surface wave spectrum compared well with measurements made by a pitch roll wave buoy maintained in the area by the national data buoy center ndbc the acoustic and buoy spectra agreed well on the location of the peak swell frequency the root mean square rms wave amplitudes from the two measures agreed to within 7

small scale sublithospheric continental mantle deformation constraints from sks splitting observations summary it is generally believed that azimuthal seismic anisotropy is an important property of the subcontinental mantle and that this is a relatively small scale phenomenon compared with anisotropy of the oceanic mantle during recent years shear wave splitting has come to be viewed as the most direct diagnostic of azimuthal anisotropy the magnitude of the splitting delay times can often rule out a crustal source of anisotropy however it is not yet clear whether the source of anisotropy lies within the lithosphere or the sublithosphere it is essential to know the scale of lateral variations and the distribution of anisotropy with depth in order to understand the origin of this phenomenon here we describe a way to constrain the depth of the anisotropy by studying the lateral variation of the splitting parameters at neighbouring seismic stations we use sks and skks phases recorded at the nars nl array a relatively dense network of broad band stations separated by about 50 km fresnel zones at different depths are calculated for these phases the depth of the anisotropy is constrained by the criterion that fresnel zones corresponding to different splitting observations should not overlap variability in the splitting measurements for events with different directions of approach recorded at one station and for single events recorded at various stations provides evidence for a non homogeneously distributed source of anisotropy located below the array at different depths we show that this variability in the measurements cannot be due to effects of crustal scattering our results indicate that the anisotropy is partly located in the sublithospheric mantle at about 400 km depth this suggests that the deformation in this region is subject to significant lateral variations

distribution of seismic anisotropy in the subduction zone beneath the wellington region new zealand summary shear wave splitting measurements from s arrivals of local earthquakes recorded at the incorporated research institutions for seismology iris broadband sensor snzo are used to determine a basic anisotropic structure for the subduction zone in the wellington region with the use of high frequency filters fast anisotropic polarization w and splitting time dt measurements typical of crustal anisotropy are evident but the larger splitting expected from the mantle is often not resolved the small splitting seen agrees well with the results of previous studies concerning shallow crustal anisotropy with the use of lower frequency filters measurements more consistent with mantle anisotropy are made anisotropy of 4 4 0 9 per cent with a fast polarization of 29 38 is calculated for the subducting slab from 20 to 70 km depth using this result in addition to the results of previous studies a model is proposed the model requires a frequency dependent anisotropy of less than 1 4 per cent when measured with a period of 2 s to be present in the sub slab mantle separate from this population a band of events in northern cook strait with an 86 10 fast polarization is seen this is at about 40 from the strike of the hikurangi margin and suggests a source of shear strain 40 removed from that found in the majority of the region the cause of this is probably a deformation in the subducting slab in this region as it moves towards a greater incline to the south

seismicity structure and tectonics of the wellington region new zealand summary the wellington region lies on the border between the pacific and australian plates the former being subducted the plate interface lies at shallow depth under the region 10 60 km so that detailed and long term monitoring of the seismicity of the zone where large thrust earthquakes are expected to occur is possible hypocentres for events recorded by a 12 station seismic network for the period 1978 1982 7231 events have been refined by the use of a threedimensional velocity model derived from a simultaneous inversion of arrival time data from 129 local events for both hypocentres and velocity the hypocentres show that the subducting pacific plate is defined by a nw dipping zone of relatively intense activity the plate interface as given by the upper envelope of this activity is offset vertically by 7 km on average along a nw se line through cook strait parallel to the direction of absolute motion of the pacific plate rather than to the direction of convergence the strike of the interface at deeper levels also changes across this line this distortion may partly reflect effects of the transition to continent continent type convergence further south within the subducting pacific plate the seismic activity defines two parallel zones separated by about 20 km the lower zone being much less active and perhaps merging with the upper zone at a depth of 80 90 km composite focal mechanisms of events in the upper zone imply predominantly down dip tension with a compression axis perpendicular to the plate interface the mechanism in the lower zone also implies down dip tension but the compression axis has rotated to be horizontal and along strike indicating a component of lateral shearing within the plate these two types of mechanism correlate well with those found previously for events at depths of 150 km or more thus it seems that the large negative buoyancy of the subducted lithosphere dominates the internal stress field throughout the length of the seismically active region the lithosphere acting as a stress guide there is no evidence for a significant number of events with mechanisms indicating thrusting along the plate interface the seismic velocity in the region corresponding to the upper seismic band is close to that typical for oceanic crust which here is probably considerably thicker than normal there is no evidence for a substantial increase in velocity down

frequency dependent anisotropy in wellington new zealand shear wave splitting measurements are made using teleseismic s scs and sks waveforms recorded at the gdsn broadband station snzo situated in south karori new zealand the average s and sks delay times are around 2 to 3 s among the highest in the world the average scs delay time is significantly smaller around 1 s this discrepancy appears to be due to differences in dominant frequency the most likely cause of frequency dependent anisotropy is oriented heterogeneities with a scale length much smaller than the wavelength the fast polarizations range between 21 and 79 with an average ne sw direction which is sub parallel to the trend of the local geologic structure and the strike of the hikurangi subduction zone azimuthal variations in delay time which cannot be explained by differences in period may be due to a dipping axis of symmetry or laterally varying anisotropy or a more complicated symmetry system

sks splitting in southern italy anisotropy variations in a fragmented subduction zone abstract in this paper we present a collection of good quality shear wave splitting measurements in southern italy in addition to a large amount of previous splitting measurements we present new data from 15 teleseisms recorded from 2003 to 2006 at the 40 stations of the cat scan temporary network these new measurements provide additional constraints on the anisotropic behaviour of the study region and better define the fast directions in the southern part of the apulian platform for our analysis we have selected well recorded sks phases and we have used the method of silver and chan to obtain the splitting parameters the azimuth of the fast polarized shear wave and delay time t shear wave splitting results reveal the presence of a strong seismic anisotropy in the subduction system below the region three different geological and geodynamic regions are characterized by different anisotropic parameters the calabrian arc domain has fast directions oriented nne ssw and the southern apennines domain has fast directions oriented nnw sse this rotation of fast axes following the arcuate shape of the slab is marked by a lack of resolved measurements which occurs at the transition zone between those two domains the third domain is identified in the apulian platform here fast directions are oriented almost n s in the northern part and nne ssw to ene wsw in the southern one the large number of splitting parameters evaluated for events coming from different back azimuth allows us to hypothesize the presence of a depth dependent anisotropic structure which should be more complicated than a simple 2 layer model below the southern apennines and the calabrian arc domains and to constrain at 50 km depth the upper limit of the anisotropic layer at least at the edge of southern apennines and apulian platform we interpret the variability in fast directions as related to the fragmented subduction system in the mantle of this region the trench parallel observed in calabrian arc and in southern apennines has its main source in the asthenospheric flow below the slab likely due to the pressure induced by the retrograde motion of the slab itself the pattern of in the apulian platform does not appear to be the direct result of the rollback motion of the slab whose influence is limited to about 100 km from the slab the anisotropy in the apulian platform may be related to an asthenospheric flow deflected by the complicated structure of the adriatic microplate or may also be explained as frozen in lithospheric anisotropy

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

the magnitude symmetry and origin of upper mantle anisotropy based on fabric analyses of ultramafic tectonites summary seismic anisotropy within the upper mantle originates from the preferred orientation of highly anisotropic single crystals the symmetry and magnitude of anisotropy depend upon 1 the volume percentages of the minerals constituting the upper mantle 2 the degree and symmetry of preferred orientation of each mineral and 3 the alignment of the minerals crystallographic axes relative to one another the nature of upper mantle anisotropy can be examined by studying mineral orientations within ultramafic rocks which were once part of the mantle petrofabric data for olivine and pyroxene have been used to obtain velocity anisotropy patterns over large regions of ultramafic rocks from the samail ophiolite oman the troodos ophiolite cyprus the bay of islands ophiolife newfoundland the twin sisters ultramafic washington usa the dun mountain ophiolite new zealand the red hills ophiolite new zealand and the red mountain ophiolite new zealand the compressional wave anisotropy calculated for these massifs ranges from 3 to 8 per cent in excellent agreement with observed seismic anisotropy in the upper continental and oceanic mantle the symmetry varies from orthorhombic to axial with the axial symmetry axis corresponding to the olivine a axes maxima and subparallel to spreading directions in oceanic upper mantle pyroxene a b and c axes maxima generally parallel olivine b c and a axes respectively and anisotropy decreases with increasing pyroxene content shear wave splitting is predicted for all propagation directions within the upper mantle symmetry is also orthorhombic or axial with the minimum difference in velocity between the two shear waves parallel to the maximum compressional wave velocity

anisotropic loci in the mantle beneath central peru abstract seismic anisotropy of the upper mantle beneath the central part of peru is examined by analyzing shear waves observed at broad band stations and a temporary array of short period stations shear wave splitting is seen on various shear phases such as direct s waves from local intermediate to deep earthquakes scs waves from regional deep earthquakes and sks waves from teleseismic earthquakes it is inferred that the shear wave anisotropy in the uppermost 100 km of the subcontinental mantle overlying the subducted nazca plate is 0 5 at most while the anisotropy in the subslab asthenosphere depth range of about 150 350 km is 2 or greater the fast shear wave polarization direction in this depth range as observed at two broad band stations is 30 40 different from the absolute motion of the nazca plate this does not fit simple two dimensional 2 d models of olivine alignment caused by slab induced mantle flow and implies either the existence of a complex flow pattern in the asthenosphere underneath the nazca plate or the presence of an unknown mechanism for the anisotropy formation the lower mantle beneath central peru is found to be effectively isotropic for nearly vertically propagated shear waves

azimuthal shear wave velocity anisotropy in the basin and range province using moho ps converted phases an eight station large aperture seismometer array was deployed in the basin and range province of west central nevada for the 1988 1989 program for array seismic studies of the continental lithosphere passcal passive source seismic experiment during the 10 months of data collection a total of 100 teleseismic events were recorded 30 source equalized p wave receiver functions show waveform variations in timing amplitude and polarity in the moho ps phase that is diagnostic of shear wave velocity anisotropy information regarding the crustal anisotropic component can be recovered through the ps phase since it is generated at the moho and is sensitive to shear wave splitting in the crust only we utilize single event stacked and array beam formed receiver functions to minimize the influence of structural variations and scattering on the observed converted phases analysis of shear wave splitting of the ps conversion from the crust mantle boundary indicates a fast azimuth of anisotropy oriented approximately nw se with observed time delays of 0 20 s this strongly suggests coherent crustal anisotropy with a fast direction perpendicular to the strike of fault block mountain ranges this fabric is most likely due to the preferred crystallographic orientation of seismically anisotropic minerals in the middle to lower crust rather than the distribution of basin and range upper crustal fracture systems

shear wave velocity estimation in porous rocks theoretical formulation preliminary verification and applications shear wave velocity logs are useful for various seismic interpretation applications including bright spot analyses amplitude versus offset analyses and multicomponent seismic interpretations measured shear wave velocity logs are however often unavailable we developed a general method to predict shear wave velocity in porous rocks if reliable compressional wave velocity lithology porosity and water saturation data are available the precision and accuracy of shear wave velocity prediction are 9 and 3 respectively the success of our method depends on 1 robust relationships between compressional and shear wave velocities for water saturated pure porous lithologies 2 nearly linear mixing laws for solid rock constituents 3 first order applicability of the biot gassmann theory to real rocks we verified these concepts with laboratory measurements and full waveform sonic logs shear wave velocities estimated by our method can improve formation evaluation our method has been successfully tested with data from several locations

seismic anisotropy inferred from wide angle reflections off lofoten norway indicative of shear aligned minerals in the upper mantle abstract data from six three component ocean bottom seismographs situated on a 175 km long profile on the continental shelf off lofoten northern norway have been modelled using 2 d seismic ray tracing and p and s wave velocity models for the upper mantle are presented a dipping structure about 10 km below the moho causes strong p and s reflections as well as p to s conversions this upper mantle reflection might be related to a possible seaward dipping master fault and or the presence of partially hydrated peridotite the modelling of these arrivals reveals significant seismic anisotropy in the upper mantle the maximum p wave anisotropy is inferred to be in excess of 15 the maximum s wave anisotropy in the order of 5 10 and the axis of the fast velocity is inferred to be inclined at least 45 from the horizontal the strong anisotropy for p waves indicates that its cause is related to alignment of minerals we postulate that the inferred steep inclination of anisotropic minerals might be a result of shearing along faults and fractures during post caledonian extensional episodes it must be emphasized however that it cannot be excluded that more complicated anisotropic models with other geological implications can also explain the observations satisfactorily the study demonstrates the importance of incorporating s waves in analysis of seismic data because the strong p wave anisotropy would not have been discovered without the constraints given by the s waves

extensive dilatancy anisotropy relative information in vsps and reflection surveys1 yardley g and crampin s 1990 extensive dilatancy anisotropy relative information in vsps and reflection surveys geophysical prospecting 39 337 355 shear waves have complicated interactions with the free surface particularly in the presence of low velocity surface layers topographic irregularities and the expected near surface crack and stress anomalies consequently it has been suggested that shear waves should be recorded subsurface in vertical seismic profiles vsps in order to extract accurate information about the in situ crack and stress geometry contained in shear wave splitting this paper compares the information in synthetic shear waves in reflection gathers and vsps in order to assess the relative merits of the two techniques for investigating shear wave splitting synthetic seismograms demonstrate that in the presence of even very simple surface layers shearwaves recorded in reflection surveys at the surface have polarizations which may not indicate crack and stress geometry at depth in contrast shear waves recorded in vsps are relatively unaffected by surface layers and near surface stress and crack anomalies and the behaviour of shear wave splitting is dominated by the structure of the rock mass in the vicinity of subsurface geophones matrix rotations of multicomponen t multisource shear wave reflection data to extract the information contained in the split shear waves are found to be directly meaningful only in situations where crack orientations do not change with depth

continental subduction and three dimensional crustal structure the northern south island new zealand the three dimensional vp and vp vs structure of a region where subduction transitions to oblique transform faulting has been determined using arrival times from 579 local earthquakes recorded during a temporary deployment and 3146 earthquakes have been relocated between 40 km and 100 km depth the subducted plate is imaged as a relatively low velocity feature in the uppermost mantle reflecting the continental nature of the subducted crust in this region an increase in amplitude of this low velocity feature from northeast to southwest can be related to an increase in the thickness of the crust of the subducted plate in this direction velocity variations within the subducted and overlying plates show some spatial correlation this suggests an interaction between the plates which extends well beyond the plate interface and is consistent with other geophysical and geological evidence that the plate interface beneath marlborough is currently not accommodating much active subduction in the overlying plate the awatere fault is a major structural feature associated with a low velocity zone extending to 23 km depth there is a marked change in structure near this fault with seismic velocities being lower to the southeast a relatively high level of seismicity occurs in this region of lower seismic velocities suggesting a relationship between the two a possible explanation for this is elevated pore pressures caused by fluids derived from dehydration of the continental subducted crust the low velocity region in the overlying plate coincides with the region of most intense active deformation suggesting it is relatively weak

anisotropic structure under a back arc spreading region the taupo volcanic zone new zealand 1 we investigate anisotropy surrounding a continental back arc spreading region the central volcanic region cvr in the north island of new zealand and present a simple method for spatial averaging to examine heterogeneous anisotropy s phases from local earthquakes yield consistent trench parallel fast directions in the southern compressional region but varied directions in the northern extensional region in the forearc east of the cvr is nearly trench parallel suggesting shear in the mantle or trench parallel flow in the western cvr local phases from 50 to 80 km depth give parallel to extension 120 and phases from 150 to 220 km depth yield similar 150 these results suggest asthenospheric flow with olivine a axes oriented in the extension direction down to 100 150 km depth in the eastern cvr local phases show mixed suggesting that it is a transition between the western cvr and the forearc trench parallel for shallow local events within the cvr may be caused by fluid filled cracks trench parallel for deeper local events and for sks phases can be explained by standard olivine fabrics developed by lithospheric shearing by trench parallel flow and by fossil anisotropy in the subducting oceanic lithosphere alternatively in the cvr they could be caused in part by hydrous minerals at the slab mantle upper interface which could cause olivine a axes to orient parallel to the trench due to a maximum shear stress oriented perpendicular to the trench along the slab we suggest that trench perpendicular flow in the extensional region helps to drive fluids away from the slab while in the compressional region trench parallel flow fails to distribute the fluids explaining the location of changes in geophysical properties

internal deformation of the subducted nazca slab inferred from seismic anisotropy subducting oceanic plates are often considered as cold rigid slabs analysis of seismic anisotropy in the subducted nazca plate beneath peru suggests that the plate has deformed internally during subduction

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

study of a novel range dependent propagation effect with application to the axial injection of signals from the kaneohe source a novel range dependent propagation effect occurs when a source is placed on the seafloor in shallow water with a downward refracting sound speed profile and sound waves propagate down a slope into deep water under these conditions small grazing angle sound waves slide along the bottom downward and outward from the source until they reach the depth of the sound channel axis in deep water where they are detached from the sloping bottom and continue to propagate outward near the sound channel axis this mudslide effect is one of a few robust and predictable acoustic propagation effects that occur in range dependent ocean environments as a consequence of this effect a bottom mounted source in shallow water can inject a significant amount of acoustic energy into the axis of the deep ocean sound channel that can then propagate to very long ranges numerical simulations with a full wave range dependent acoustic model show that the kaneohe experiment had the appropriate source bathymetry and sound speed profiles that allows this effect to operate efficiently this supports the interpretation that some of the near axial acoustic signals received near the coast of california from the bottom mounted source located in shallow water in kaneohe bay oahu hawaii were injected into the sound channel of the deep pacific ocean by this mechanism numerical simulations suggest that the mudslide effect is robust

calibration of shear wave splitting in the subcontinental upper mantle beneath active orogenic belts using ultramafic xenoliths from the canadian cordillera and alaska abstract despite the abundance of measurements of shear wave splitting around the world some fundamental questions about its geological interpretation have not yet been answered in order to constrain 1 i the orientation and magnitude of s wave anisotropy 2 ii the thickness of the anisotropic layer and 3 iii the possible variation of the anisotropy with depth in the subcontinental upper mantle beneath active orogenic belts we have carried out a systematic investigation on shear wave properties of mantle xenoliths from recent alkaline basalts of the canadian cordillera and alaska it is found that the polarization direction of the fast s wave as well as the time delay t between the two arrivals are strongly dependent on the propagation direction with respect to the structural frame at a single station the variations of measured and t values for sks phases from different events can result from this dependence because the average grain size of olivine in upper mantle samples is commonly larger than 0 5 mm dislocation creep prevails over a thickness of at least 250 km and results in olivine lpo and hence anisotropy and shear wave splitting a single thick anisotropic layer in mountain parallel subvertical shear zones may cause high values of tgd while multiple subhorizontal layers will produce multiple splitting making the measurement and interpretation of s wave splitting difficult our results also suggest that the s wave splitting observed in cold old stable cratons such as the canadian shield is likely dominated by fossil anisotropy whereas that observed in hot active orogenic regions such as the northwest american cordillera and the tibetan plateau may be dominated by modern tectonic deformation in the lithosphere and flow in the asthenosphere

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

evidence for shallow and pervasive seismic anisotropy in the wellington region new zealand the shear waves of local earthquakes were recorded during a 5 month deployment of seven three component digital seismographs on the wellington peninsula new zealand the seismographs were spaced an average of less than 5 km apart and over 300 local earthquakes were recorded with phase arrivals within the shear wave window a significant number 37 of the earthquakes recorded showed clear evidence of shear wave splitting identifiable fast and slow shear wave arrivals with similar pulse shapes consistent polarization directions at particular stations were also observed even when poor signal to noise ratio or scattering meant that no slow shear wave arrival could be identified however there were large station to station differences in the polarization directions correcting for the observed shear wave splitting improved the fit between the measured shear wave polarizations and those calculated assuming a double couple focal mechanism the cause of the observed shear wave splitting is therefore most likely to be seismic anisotropy large station to station differences in the polarization alignments ranging from 61 24 to 137 18 suggest that most anisotropy is confined to the top 2 3 km of the crust however there is evidence from one station for a small amount of pervasive anisotropy if such a trend existed on the other stations it could not be identified because of the large scatter in the data points the measured delay times between split shear waves vary from 0 02 to 0 22 s with a mean value of 0 1 0 06 s this translates to a near surface shear wave velocity anisotropy of about 10 with up to 2 pervasive anisotropy possible throughout the crust this data set indicates that extensive dilatancy anisotropy cannot be the sole cause of crustal seismic anisotropy and that foliations in the rock fabric and the fracture zones of active faults may also be important there is no evidence for temporal change in the shear wave splitting parameters during the period of the experiment

scs wave splitting of deep earthquakes around japan summary scs wave splitting of five deep earthquakes in subduction zones near japan is investigated using horizontal seismograms recorded al jma stations for each earthquake we clearly observe uniform scs wave splitting in all stations over japan especially for the events located south of honshu in 1982 1984 and 1993 however the directions of fast polarized waves of these events differed by a maximum of about 50 from one another the orientation of fast polarized waves in the 1982 event was nnw sse those in the two later events wnw ese we also recognize this discrepancy in the results of the analysis of the 1971 sea of okhotsk event reported by fukao 1984 the sakhalin islands event in 1990 reveals a linear particle motion without such a change in direction of the second arrivals implying no anisotropy these observations are interpreted as indicating an anisotropic region within the slab near the earthquake sources but not beneath the receivers since the orientations of fast polarized waves recorded at each station are not common to all the earthquakes furthermore we consider that anisotropy exists non uniformly within the slab the event in 1982 which occurred in almost the same area as those in 1984 and 1993 showed a fast direction different from the events in 1984 and 1993 the 1982 event was 179 km deep but the two later events were at 398 km and 360 km respectively the fast direction observed from the 1982 event is parallel to the fossil plate motion whereas those from the events in 1984 and 1993 are parallel to the compression axis within the subducting slab the depth of 400 km is a phase boundary where olivine changes to spinel we consider that the most likely cause of the change in anisotropy direction is the re orientation of crystals associated with the phase change of olivine to spinel due to subduction of the slab

numerical modeling of acoustic tomography in the straits of florida sensitivity to bathymetry an acoustic model based on the split step parabolic equation technique is used to study the effects of irregularities in bathymetry on the transmission loss or equivalently sound energy distribution for continuous wave cw acoustic signals in particular contours of transmission loss in the range depth plane are predicted for the propagation of a 510 hz cw signal across the straits of florida between miami florida and bimini bahamas the predicted intensity contours exhibit recognizable energy bundles that correspond to actual ray bundles the stability of these bundles is studied as random bottom irregularities are introduced the distribution of acoustic energy is found to be insensitive to random variations in bottom bathymetry of up to 5 m rms the same model is used to investigate the feasibility of transport and depth averaged temperature measurements of the florida current by acoustic method energy contour plots are predicted for several source receiver geometries a system of five transceivers positioned to form two triangles sharing a common vertex in the eastern edge of the miami terrace is evaluated as a monitoring strategy

the fracture criticality of crustal rocks summary the shear wave splitting observed along almost all shear wave ray paths in the earth s crust is interpreted as the effects of stress aligned fluid filled cracks microcracks and preferentially oriented pore space once away from the free surface where open joints and fractures may lead to strong anisotropy of 10 per cent or greater intact ostensibly unfractured crustal rock exhibits a limited range of shear wave splitting from about 1 5 to 4 5 per cent differential shear wave velocity anisotropy interpreting this velocity anisotropy as normalized crack densities a factor of less than two in crack radius covers the range from the minimum 1 5 per cent anisotropy observed in intact rock to the 10 per cent observed in heavily cracked almost disaggregated near surface rocks this narrow range of crack dimensions and the pronounced effect on rock cohesion suggests that there is a state of fracture criticality at some level of anisotropy between 4 5 and 10 per cent marking the boundary between essentially intact and heavily fractured rock when the level of fracture criticality is exceeded cracking is so severe that there is a breakdown in shear strength the likelihood of progressive fracturing and the dispersal of pore fluids through enhanced permeability the range of normalized crack dimensions below fracture criticality is so small in intact rock that any modification to the crack geometry by even minor changes of conditions or minor deformation particularly in the presence of high pore fluid pressures may change rock from being essentially intact below fracture criticality to heavily fractured above fracture criticality this recognition of the essential compliance of most crustal rocks and its effect on shear wave splitting has implications for monitoring changes in any conditions affecting the rock mass these include monitoring changes in reservoir evolution during hydrocarbon production and enhanced oil recovery and in monitoring changes before and after earthquakes amongst others