the philippine sea new survey results reveal the structure and the history of the marginal basins the hydrographic department of japan initiated a long term geological and geophysical survey of the philippine sea in 1983 we here summarize the bathymetry and magnetic anomalies of this area focusing on the structure and the evolutionary process of the marginal basins the shikoku and parece vela basins developed as one backarc system in the later phase of their formation the significant difference between the two basins however is clearest near the extinct spreading centers the steep oki daito escarpment located in the northern west philippine basin was mapped for the first time during this survey the n s trending spreading fabric north of the escarpment is in striking contrast to the nw se to e w lineament in the south our data reveal a more complicated history of the western philippine sea than has been proposed in previous studies

including anisotropy in 3 d velocity inversion and application to marlborough new zealand summary two types of seismic velocity heterogeneity are often observed in a given region pronounced patterns of anisotropy shown by local earthquake shear wave splitting and seismic velocity structure imaged by local earthquake traveltimes we seek to combine these types of observations by including anisotropy in 3 d velocity inversion the fast polarization directions from local earthquake shear wave splitting are used to define an initial anisotropy model the ray paths are estimated for each of the shear wave splitting observations and the 3 d initial model is determined via an averaging process and a coherence analysis the initial magnitude of anisotropy at a node is dependent on the number of adjacent ray paths and the consistency of the shear wave splitting anisotropy measurements the traveltime inversion is parametrized with an isotropic component and two azimuthal anisotropy parameters for each node the ray paths are defined with approximate 3 d ray tracing and for each point along a ray path the velocity is defined as a function of the ray path azimuth at that point a blended approach is used to damp the perturbations to the magnitude of anisotropy and or perturbations to the fast direction tests with synthetic traveltime data show that the anisotropic velocity inversion method is reliable application to the marlborough region shows variations in anisotropy in the brittle crust ductile lower crust mantle and subducted slab the patterns are consistent with the shear wave splitting observations but show some additional features the maximum anisotropy is approximately 12 per cent oriented northeast between the awatere and wairau faults this high anisotropy may reflect deformation and extensive fracturing in a region of high total shear strain in the region of interaction of the ductile lower crust with the shallow subducting slab moderate anisotropy with east west orientation is imaged this is consistent with ductile deformation roughly parallel to the plate velocity

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

a new mesozoic isochron chart of the northwestern pacific ocean paleomagnetic and tectonic implications a new mesozoic isochron chart of the northwestern pacific ocean has been completed by utilizing all the currently available data from 357 cruises our chart provides the first whole view of the mesozoic magnetic anomaly lineations in the northwestern pacific ocean and yields essential information for understanding the problems of the mesozoic magnetostratigraphic time scale and the evolutional history of the pacific plate the mesozoic magnetostratigraphic time scale must be recalibrated using not only the hawaiian lineation set but also other mesozoic lineation sets the lack of lineation m0 in the lineation sets other than the hawaiian lineation set indicates that a reorganization of the ridges except the pacific farallon ridge occurred between chrons m1 and m0 nearly 120 ma the eruption rate of the shatsky hotspot seems to be comparable to that of the reunion hotspot at their initial stage the growth of the pacific plate began in a small triangular region at about 190 ma

the colima volcanic complex mexico part ii late quaternary cinder cones since the late pleistocene eleven cinder and lava cones have erupted on the floor of the southern colima graben ne and nw of the large active andesitic volcano colima scoria and lava samples from nine of the cones form a completely transitional basic alkalic series including basanites 9 leucite basanites 3 and minettes 15 the commonest variety of mica lamprophyre these samples represent primitive high temperature magmas with 47 6 50 3 sio2 7 4 15 3 mgo 2 5 4 4 k2o and 2 2 9 9 normative nepheline all members of this basic alkalic suite contain mg olivine fo75 94 chromite augite and late plagioclase and titanomagnetite the petrographic transition from basanite to minette is marked by the appearance of sanidine and the volatile bearing phases phlogopite apatite and analcime during late stages of crystallization as these phases increase in abundance presumably reflecting a rise in magmatic volatile content there are corresponding increases in the whole rock concentrations of 16 incompatible elements although these incompatible elements are relatively abundant even in the basanites many are highly concentrated in the minettes ba 4 200 ppm sr 3 100 ppm zr 550 ppm ce 190 ppm hf 18 ppm among the incompatible elements the degrees of enrichment in the minettes relative to the basanites decrease in the order h th ce la nd zr hf u ba sm eu pb p nb sr ti these enrichments may reflect the increasing importance of minor incompatible element rich mantle phases during partial melting the concentrations of alkali metals k na rb and cs do not correlate with these other elemental enrichments the leucite basanties have similar incompatible element contents to many minettes differing from them only in the presence of leucite rather than analcime and ti f rich groundmass phlogopite rather than hydrous phlogopite phenocrysts thus the leucite basanites represent relatively dry equivalents of minettes two of the eleven cinder cones are calc alkaline in nature and do not belong to the basanite minette group the easternmost cone is constructed of high al basalt and the northernmost of basaltic andesite the high al basalt 49 5 sio2 9 3 mgo 221 ppm ni closely approximates a parental magma to the post caldera andesitic suite of v colima 56 5 61 6 sio2 the basaltic andesite is relatively enriched in incompatible elements compared to the high al basalt v colima trend the ne normative basanite minette samples are highly enriched in incompatible elements while the contemporaneous hy qz normative calc alkaline suite encompassing the high al basalt and v colima s andesites is characterized by relatively low abundances of these elements no likely mineral assemblage can relate the alkaline and calc alkaline suites through crystal fractionation they probably represent fundamentally different melting events during the quaternary the main focus of andesitic volcanism in the southern colima graben has migrated southward with time volcan colima marks its present position 5 km south of the pleistocene volcano nevado de colima and another 15 km from the still older volcan cantaro the eruptions of basic alkalic magma probably occurred during the late stages of nevado s life and through the life of v colima they generally migrated from west to east with time towards v cantaro the most recent cone v apaxtepec is the only one east of the andesitic colima cantaro axis the oldest and the two youngest cones produced basanites while minettes dominated at cones of intermediate ages the cinder cone eruptions may have coincided with a phase of lamprophyre dike injection into plutons solidifying beneath the extinct volcanoes north of v colima the southern end of the colima graben can be viewed then as the volcanic analog of many classical post plutonic hypabyssal lamprophyre localities

seismicity and structure in central mexico evidence for a possible slab tear in the south cocos plate the morphology of the transition from flat to normal subduction in eastern central mexico is explored using intraslab earthquakes recorded by temporary and permanent regional seismic arrays observations of a sharp transition in slab dip near the abrupt end of the trans mexican volcanic belt tmvb suggest a possible slab tear located within the subducted south cocos plate the eastern lateral extent of a thin ultra slow velocity layer usl imaged atop the cocos slab in recent studies along the meso america subduction experiment array is examined here using additional data we find an end to this usl which is coincident with the western boundary of a zone of decreased seismicity and the end of the tmvb near the sharp transition in slab dip waveform modeling of the 2 d structure in this region using a finite difference algorithm provides constraints on the velocity and geometry of the slab s seismic structure and confirms the location of the usl analysis of intraslab seismicity patterns reveals clustering sudden increase in depth variable focal mechanism orientations and faulting types and alignment of source mechanisms along the sharp transition in slab dip the seismicity and structural evidence suggests a possible tear in the south cocos slab this potential tear together with the tear along the orozco fracture zone to the northwest indicates a slab rollback mechanism in which separate slab segments move independently allowing for mantle flow between the segments

the relationship between quaternary volcanism in central mexico and the seismicity and structure of subducted ocean lithosphere late quaternary volcanism in central mexico is related to the subduction of young ocean lithosphere at the middle america trench along arc variations in seismicity volcano structure and composition of volcanic products bear a remarkable correlation with the age and structural framework of the downgoing slab morphological and petrographic characteristics of major volcanoes within the trans mexican volcanic belt tmvb serve to distinguish two calc alkaline subprovinces 1 a western arc averaging 60 km in width associated with aseismic subduction of the rivera plate the main cones of this region are dominated by two pyroxene andesites comprise volumes 70 km 3 and stand less than 3 000 m above sea level 2 a broad central and eastern arc related to subduction of a gently inclined segment of the cocos plate bounded by the rivera transform and the tehuantepec ridge major volcanic edifices possess summit elevations in the range 4 000 to 6 000 m have appropriately larger volumes typically 200 km 3 and are constructed with a high proportion of amphibole bearing lavas the boundary between these subprovinces is marked by a north south oriented structural depression the colima graben and it coincides with a 100 km offset in the volcanic front extensional tectonism in the colima graben accompanied by mixed calc alkaline and alkaline volcanism of potassic affinity is likely related to a hinge type transform fault which marks the cocos rivera plate juncture in the downgoing slab a third segment of ocean floor is presently interacting with continental lithosphere south of the gulf of tehuantepec where quaternary volcanism is weakly developed within a tectonically complex region that marks the diffuse cocos noam caribbean triple junction the northern limit of this triple junction is defined by the seismically active isthmus fault which may be related to alkaline volcanic activity at san andres tuxtla a tectonic reconstruction based on the evolution of oceanic crust reveals that the distribution of intermediate depth earthquakes along the arc is directly dependent upon the age of the subducted slab ocean lithosphere younger than approximately 20 m y is subducted aseismically at convergence rates approaching 9 cm yr the length of the inclined seismic zone indicates that the time constant for thermal relaxation in the slab is approximately 4 m y the tmvb overlies the aseismic extension of this young ocean lithosphere several aspects of this study have a bearing on the segmented nature of converging margins in general 1 the tectonic evolution of the ocean floor may determine the nature of segmentation at the site of subduction 2 the complete record of volcanism in the tmvb over the past million years can be related to the present plate configuration 3 alkaline and calc alkaline volcanism have developed contemporaneously at a converging plate margin 4 lineaments in volcanic arcs may reflect the structural complexity of the crust rather than segment boundaries in the subducted slab

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

experimental constraints on the low pressure evolution of transitional and mildly alkalic basalts the effect of fe ti oxide minerals and the origin of basaltic andesites one atmosphere melting experiments controlled at the fayalite magnetite quartz oxygen buffer on mildly alkalic and transitional basalts from iceland show that these begin to crystallize fe ti oxide minerals magnetite and or ilmenite at 1105 5 c apparently independently of bulk composition and the order of silicate and oxide mineral crystallization most samples crystalline plagioclase and olivine as the first two crystalline phases augite as the third phase and an fe ti oxide mineral as the fourth phase the main effects of fe ti oxide crystallization are a marked decrease in feo and tio2 in the liquid and a notable increase in sio2 and al2o3 and the minor oxides k2o and p2o5 with decreasing temperature the most silicic glasses are compositionally mugearitic and shoshonitic basaltic andesites because the smallest amount of glass that could be analyzed with the microprobe represents 20 55 percent liquid remaining it can be expected that more silicic liquids will occur at lower temperatures on normative pseudoternary projections the general effect of fe ti oxide crystallization for mildly alkalic and transitional basalts is a marked increase in normative quartz this is caused by a strong systematic convergence with the appearance of fe ti oxides of the bulk solid precipitates toward the liquid compositions as projected on the triangle plagioclase diopside olivine for alkalic basalts the bulk solid precipitate shows an increase in normative diopside with falling temperature and fe ti oxide crystallization this causes the liquids to move toward decreasing normative diopside and relatively little variation in nepheline the experimental observations imply that mildly alkalic and transitional magmas without stabilizing a fe ti oxide mineral will not evolve toward early silica saturation

compressional and shear wave anisotropy in the oceanic lithosphere the ngendei seismic refraction experiment summary synthetic seismogram modelling of seismic refraction data from the 1983 ngendei expedition to the south pacific indicate seismic anisotropy both within the top kilometre of oceanic crust and in the uppermost mantle we calculate average p and s wave velocity versus depth functions for two orthogonal azimuths n3ooe approximately aligned with the fast mantle direction and n120 e close to the fast crustal direction probable lateral heterogeneities at the ngendei site prevent perfect matching of data and synthetic waveforms crustal anisotropy is indicated by a 28 pattern of both p and s wave travel times as a function of azimuth and is consistent with an approximate 0 2 km s difference in p wave velocities and 0 1 km s difference in s wave velocities upper mantle anisotropy is characterized by pn velocity differences of 7 95 8 4 km s l and a nearly uniform s velocity of 4 65 km s we use these velocity profiles and traveltime data to calculate bounds on the properties of the elastic constants of the crust and upper mantle the crustal anisotropy can be explained by a model involving aligned cracks parallel to the original spreading ridge resulting in a fast direction perpendicular to the fossil spreading direction the upper mantle anisotropy is consistent with aligned olivine crystal models in which the fast direction is parallel to the fossil spreading direction for both of these models we use bounds from the ngendei data to place constraints on the physical properties of the lithosphere

stratification of anisotropy in the pacific upper mantle 1 on the basis of the use of broadband 25 150 s rayleigh wave group speeds to estimate the 2 component of azimuthal anisotropy we present evidence for a stratification of anisotropy in the uppermost mantle at large scales across the pacific basin we confirm previous surface wave studies that established that the fast axis directions of azimuthal anisotropy for intermediate and long period rayleigh waves approximately align with present day plate motions at shorter periods 25 50 s however fast axes align nearer to the paleospreading or fossil spreading direction which differs from present day plate motion in the old pacific these observations as well as observations of the age dependence of the amplitude of azimuthal anisotropy imply that azimuthal anisotropy in the pacific upper lithosphere 100 km depth is fixed or fossilized on average reflecting the strain conditions extant during the early evolution of the lithosphere rather than the current ambient flow direction in the deeper lithosphere and asthenosphere anisotropic fast axis directions align nearer to present day plate motions apparently having reoriented to conform to the current conditions of mantle flow the mechanism of anisotropy stratification remains unclear but observations are consistent with the anisotropy of the shallow lithosphere being fixed because the shear flows that can produce dislocation creep and a change in anisotropy will occur at increasing depths as the plate ages

p wave anisotropic tomography of the crust and upper mantle under hokkaido japan abstract we determined a 3 d anisotropic velocity structure beneath hokkaido in northern japan using a large number of first p wave arrival times the travel time inversion is parameterized with an isotropic component and two anisotropic parameters for each grid node by assuming p wave azimuthal anisotropy with hexagonal symmetry axis distributed horizontally the geometries of the conrad and moho discontinuities and the upper boundary of the subducting pacific plate are taken into account in the model our results show that low velocity zones exist beneath active arc volcanoes in the crust and in the central portion of mantle wedge above the high velocity pacific slab which are similar to the previous tomographic images p wave anisotropy is revealed in the crust mantle wedge and the subducting pacific slab in the upper crust the anisotropy is possibly caused by microcracks and cracks in the lower crust the anisotropic structure is different from that in the upper crust and it may be affected considerably by the plastic flow deformation in the mantle wedge the fast velocity direction fvd is generally trench normal in most of the study area while it becomes trench parallel in the pacific coast areas close to the trench which agrees with the shear wave splitting results beneath the volcanic front the fvd is also trench parallel which may suggest the existence of complex 3 d mantle flow in the mantle wedge the fvd in the subducting pacific slab is nearly north south which likely keeps the original fossil anisotropy of the pacific plate formed at the mid ocean ridge or is affected by the olivine fabric transition from a type to b type fabric due to the changes in water content stress and temperature

75 plus years of anisotropy in exploration and reservoir seismics a historical review of concepts and methods the idea that the propagation of elastic waves can be anisotropic i e that the velocity may depend on the direction is about 175 years old the first steps are connected with the top scientists of that time people such as cauchy fresnel green and kelvin for most of the 19th century anisotropic wave propagation was studied mainly by mathematical physicists and the only applications were in crystal optics and crystal elasticity during these years important steps in the formal description of the subject were made at the turn of the 20th century rudzki stressed the significance of seismic anisotropy he studied many of its aspects but his ideas were not applied research in seismic anisotropy became stagnant after his death in 1916 beginning about 1950 the significance of seismic anisotropy for exploration seismics was studied mainly in connection with thinly layered media and the resulting transverse isotropy very soon it became clear that the effect of layer induced anisotropy on data ac

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

recent volcanism in the puyehue cordon caulle region southern andes chile 40 5 s petrogenesis of evolved lavas

anisotropic structures of the upper mantle beneath the northern philippine sea region from rayleigh and love wave tomography abstract seismic anisotropy can provide fundamental information on past and present day deformation processes in the upper mantle using rayleigh and love waves recorded by land and seafloor broadband seismometers we analyzed the isotropic and anisotropic shear wave velocity structures in the northern philippine sea region we found that the fast directions of azimuthal anisotropy are parallel to the directions of ancient seafloor spreading in the lithosphere of the shikoku and west philippine basins and pacific ocean whereas they are parallel to the direction of the present day absolute plate motion apm in the asthenosphere of the shikoku basin and oblique to the direction of the apm in the pacific ocean by 430 and in the northern part of the west philippine basin by 455 in the subduction zones around the philippine sea plate the fast direction of azimuthal anisotropy is trench parallel in the ryukyu arc and oriented nw se in the izu ogasawara island arc the philippine sea plate which is a single plate shows very large lateral variations in azimuthal and radial anisotropies compared with the pacific plate

mantle and slab contributions in arc magmas destructive plate margins are major sites of terrestrial magmatism that have long had a key role in models for the generation of continental crust and the development of chemical heterogeneities in the upper mantle in these models it is important to evaluate the nature and size of contributions from preexisting crust in the generation of magmas above subduction zones and to constrain the size of the fluxes of different elements from the subducted crust and the overlying mantle wedge in arc magmas kay 1 980 hawkesworth ellam 1 989 hawkesworth et al 1 99 1 a mcculloch gamble 1 99 1 the more mobile components of the subducted oceanic crust are the products of element fractionation processes in the upper continental crust the ocean basins and hydrothermal systems on the ocean floor they therefore have distinctive isotope and even trace element ratios that in principle should be readily distinguished from those in peridotite in the mantle wedge prior to subduction in practice however element ratios are almost certainly fractionated during shallow level dehydration of sub ducted material and at greater depths where material from the slab appears to contribute in the generation of arc magmas moreover the isotope signal of slab derived material may be heavily diluted by exchange with

tectonic interpretation of the trans mexican volcanic belt discussion discussion sur l article de shurbet d h et cebull s e paru dans tectonophysics nld 1984 101 159 165

contemporaneous calc alkaline and alkaline volcanism at sanganguey volcano nayarit mexico sanganguey volcano located in the northwestern segment of the mexican volcanic belt is a calc alkaline stratovolcano that has erupted a chemically continuous suite of andesites and dacites during the past 200 000 yr alkalic volcanism however appears to have been predominant in the area since the pliocene alkalic basalts and their derivatives erupted along extensional fractures trending northwest through sanganguey and continued to erupt through the flanks of the stratocone following cessation of calc alkaline activity the calc alkaline andesite dacite volcanism was focused beneath sanganguey and appears to have been controlled by crystal fractionation in a shallow magma chamber no basalts that erupted in the area have chemical compositions that could be parental to the sanganguey calc alkaline suite the most basic calc alkaline rocks basaltic andesites have incompatible trace element concentrations nearly identical to the most basic alkaline rocks implying that the two suites were derived from distinct parental magmas alternatively it is suggested that stagnation of mantle derived alkalic basalts provided heat input to mobilize lower crustal rocks resulting in the generation of calc alkaline magmas

stress distribution and subduction of aseismic ridges in the middle america subduction zone the regional distribution of stresses associated with the subduction of the cocos plate is inferred from a synthesis of 190 earthquake focal mechanisms body and surface wave analyses of large earthquakes and seismicity distributions broad patterns of consistent behavior are found across the region from the rivera plate boundary in the northwest to the guatemala el salvador border in the southeast and are used as a framework to evaluate evidence for variations in local stresses due to the subduction of two aseismic ridges the tehuantepec ridge and the orozco fracture zone information which bears on the seismic potential at locations of aseismic ridge subduction is particularly important in that no large ms 7 5 earthquakes have occurred historically we identify three major zones with consistent patterns in focal mechanisms and hypocentral distributions of seismicity the first closest to the trench and reflecting the mechanical interaction of the converging plates is a zone of shallow thrust earthquakes extending 100 150 km inland from the trench the second is a zone of normal faulting beginning at about 200 km inland from the trench h 60 km which extends continuously along the entire length of the descending plate throughout the region the third distinct zone exhibits a relatively low level of activity and separates the zones of thrust and normal faulting at about 150 200 km inland from the trench this zone extends from the rivera plate boundary in the northwest to the guatamala region in the southeast at this point the quiet region pinches out and the thrust and normal faulting zones abut and overlap superimposed on this overall pattern we find locally only minor changes in areas of aseismic ridge subduction aside from the prominent seismic slip gaps furthermore on october 25 1981 the playa azul earthquake ms 7 3 occurred in the midregion of the orozco fracture zone body and surface wave analyses of this event show a simple source rupture and shallow thrust fault mechanism as are found elsewhere in the region the seismic moment is mo 1 3 1020 n m the calculated stress drop is 4 5 mpa not extraordinarily high as might be expected in these pervasive seismic gaps an event in the tehuantepec ridge region on january 24 1983 ms 6 7 was a large normal faulting event but further interpretation is ambiguous due to the proximity of other normal faulting we conclude that while aseismic slip may be occurring in the areas of ridge subduction the possibility of large thrust earthquakes cannot be ruled out due to the overall similarities with adjacent regions

75 plus years of anisotropy in exploration and reservoir seismics a historical review of concepts and methods geophysics 70 9nd 23nd on page 18nd right column the phrase in the valhall oil field in norway was incorrect we apologize for any inconvenience caused by this error

seismic consequences of warm versus cool subduction metamorphism examples from southwest and northeast japan warm and cool subduction zones exhibit differences in seismicity seismic structure and arc magmatism which reflect differences in metamorphic reactions occurring in subducting oceanic crust in southwest japan arc volcanism is sparse and intraslab earthquakes extend to 65 kilometers depth in northeast japan arc volcanism is more common and intraslab earthquakes reach 200 kilometers depth thermal petrologic models predict that oceanic crust subducting beneath southwest japan is 300 to 500 c warmer than beneath northeast japan resulting in shallower eclogite transformation and slab dehydration reactions and possible slab melting

seismotectonics beneath the tokyo metropolitan area japan effect of slab slab contact and overlap on seismicity 1 we first determine the configuration of the upper surface of the pacific pac slab beneath kanto japan from the distribution of interplate earthquakes relocated by an appropriate 1 d velocity model then traveltime tomography is carried out to estimate three dimensional seismic velocity structures around kanto using 735 520 p wave and 444 049 s wave arrival times from 6508 local earthquakes the obtained results suggest that the philippine sea phs slab is subducting to depths of 130 140 km without a gap even to the northwest of the izu collision zone we subsequently define the lateral extent of the contact zone between the bottom of the phs slab and the upper surface of the pac slab phs pac interface and reveal that the slab contact zone underlies a wider area beneath kanto in harmony with the kanto plain the downdip limit of interplate thrust type earthquakes on the pac slab is deepened by 30 km locally under the slab contact zone this deepening is probably caused by a lower temperature environment in the pac slab resulting from the overlap with the phs slab subducting above and consequent thermal shielding by the phs slab from the hot mantle wedge we detect an extremely low velocity anomaly in the easternmost portion of the phs slab which is probably attributable to serpentinization of mantle peridotite interplate earthquakes are almost absent along the phs pac interface overlain by the serpentinized mantle in the phs slab suggesting that ductile deformation takes place along the interface because of low viscosity of the serpentine

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

can the okhotsk plate be discriminated from the north american plate the plate geometry in northeast asia has been a long standing question with a major issue being whether the sea of okhotsk and northern japanese islands are better regarded as part of the north american plate or as a separate okhotsk plate this question has been difficult to resolve because earthquake slip vectors along the kuril and japan trenches are consistent with either pacific north america or pacific okhotsk plate motion to circumvent this difficulty we also use slip vectors of earthquakes along sakhalin island and the eastern margin of the japan sea and compare them to the predicted eurasia okhotsk and eurasia north america motions for a model with a separate okhotsk plate we invert 10 eurasia okhotsk and 255 pacific okhotsk slip vectors with pacific north america and eurasia north america nuvel 1 data alternatively for a model without an okhotsk plate those eurasia okhotsk and pacific okhotsk data are regarded as eurasia north america and pacific north america data respectively the model with an okhotsk plate fits the data better than one in which this region is treated as part of the north american plate because the improved fit exceeds that expected purely from the additional plate the data indicate that the okhotsk plate can be resolved from the north american plate the motions on the okhotsk plate s boundaries predicted by the best fitting euler vectors are generally consistent with the recent tectonics the eurasia okhotsk pole is located at northernmost sakhalin island and predicts right lateral strike slip motion on the nne striking fault plane of the may 27 1995 neftegorsk earthquake consistent with the centroid moment tensor focal mechanism and the surface faulting along the northern boundary of the okhotsk plate the north america okhotsk euler vector predicts left lateral strike slip consistent with the observed focal mechanisms on the nw boundary of the okhotsk plate the eurasia okhotsk euler vector predicts e w extension discordant with the limited focal mechanisms and geological data this misfit may imply that another plate is necessary west of the magadan region in southeast siberia but this possibility is hard to confirm without further data such as might be obtained from space based geodesy

geology and k ar dating of the tuxtla volcanic field veracruz mexico the tuxtla volcanic field tvf is located on the coast of the gulf of mexico in the southern part of the state of veracruz mexico volcanism began about 7 my ago in the late miocene and continued to recent times with historical eruptions in ad 1664 and 1793 the oldest rocks occur as highly eroded remnants of lava flows in the area surrounding the historically active cone of san martin tuxtla between about 3 and 1 my ago four large composite volcanoes were built in the eastern part of the area rocks from these structures are hydrothermally altered and covered with lateritic soils and their northern slopes show extensive erosional dissection that has widened preexisting craters to form erosional calderas the eastern volcanoes are composed of alkali basalts hawaiites mugearites and benmoreites with less common calc alkaline basaltic andesites and andesites in the western part of the area san martin tuxtla volcano and its over 250 satellite cinder cones and maars produced about 120 km3 of lava over the last 0 8 my a ridge of flank cinder cones blocked drainage to the north to form laguna catemaco lavas erupted from san martin and its flank vents are restricted to compositions between basanite and alkali basalt the alignment of major volcanoes and flank vents along a n55 w trend suggests an extensional stress field in the crust with a minimum compressional stress orientation of n35 e in total about 800 km3 of lava has been erupted in the tvf in the last 7 my this gives a magma output rate of about 0 1 km3 1000 year a value smaller than most composite cones but similar to cinder cone fields that occur in central mexico individual eruptions over the last 5000 years had volumes on the order of 0 1km3 with average recurrence intervals of 600 years the alkaline compositions of the tvf lavas contrast markedly with the calc alkaline compositions erupted in the subduction related mexican volcanic belt to the west leading previous workers to suggest that the tvf is not related to subduction trace element signatures of tvf lavas indicate however that they are probably related to subduction we suggest that the alkaline character of the tvf lavas is the result of low degrees of melting of a mantle source coupled with a stress regime that allows these small volume melts to reach the surface in the tvf

p wave anisotropy tomography of central japan insight into subduction dynamics abstract a p wave anisotropic tomography under central japan was determined using a large number of arrival time data from 22 244 earthquakes recorded by 770 seismic stations on the japan islands the results show that the anisotropy exists widely beneath the study area except for the central portion of the mantle wedge the fast velocity directions fvds of the mantle wedge are generally trench parallel under the fore arc area and trench normal under the back arc area which may be attributed to the lattice preferred orientation lpo of olivine changing from b type under the fore arc area to a type under the back arc area with the variation of water content and the occurrence of mantle flow the subducting philippine sea plate phs beneath tokai is revealed as a high velocity anomaly dipping northward with nee to ne sw fvds which are consistent with the spreading direction of the phs during its formation and are inferred as the fossil anisotropy of the phs the anisotropy of the upper portion of the pacific slab under central japan is also revealed in this study the fvds of the upper portion of the pacific slab are sub parallel to the japan trench and izu bonin trench at 60 120 km depths and become almost perpendicular to the trenches below 100 120 km depth after comparisons with the mesozoic magnetic anomaly lineations and the principle stresses in the pacific slab we propose that the fossil fvds of the pacific slab beneath central japan were already rebuilt by the principle stresses in the slab

tearing of stagnant slab subducted slabs of oceanic lithosphere below the western pacific tend to be stagnant in the transition zone with poorly known mechanical properties typical examples are the izu bonin and japan slabs that meet each other to form a cusplike junction beneath southwest japan here we show that these two slabs are torn apart at their junction when they bend to flatten over the 660 kilometer discontinuity as is expected from a simple geometric argument we present three lines of evidence for this ongoing slab tear

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

slab segmentation revealed by anisotropic p wave tomography 1 seismic anisotropy is a useful indicator for identifying the physical and chemical condition of the earth s interior such as stress and flow fields and in situ constituent minerals using traveltime tomography we examined three dimensional anisotropic p wave velocity structure of the kii peninsula southwest japan where source regions of megathrust earthquakes along the nankai trough are presumed to exist the tomography revealed that the philippine sea slab beneath the peninsula is segmented i e a high velocity slab with e w anisotropy obtained during seafloor spreading is broken by an anomalous low velocity region with n s anisotropy the anomaly lies along the segmentation boundary of two previous rupture zones the 1944 tonankai and 1946 nankai earthquakes and can be explained by the occurrence of a fracture zone with n s oriented fractures and fault planes revealed by wide angle seismic data it is likely that this anomalous region will form a clear segmentation boundary during future earthquakes

shear wave splitting and waveform complexity for lowermost mantle structures with low velocity lamellae and transverse isotropy 1 shear waves that traverse the lowermost mantle exhibit polarization anomalies and waveform complexities that indicate the presence of complex velocity structure above the core mantle boundary synthetic seismograms for horizontally and vertically polarized shear waves sh and sv respectively are computed using the reflectivity method for structures with low velocity sheets lamellae and for comb like models approximating long wavelength vertical transverse isotropy vti motivated by evidence for partial melt in the deep mantle lamella parameter ranges include 1 vp from 5 to 10 vs 3 vp 2 100 to 300 km thickness of vertical stacks of lamella 3 lamella spacing and thickness varying from 0 5 to 20 km and 4 lamellae concentrated near the top bottom or throughout the d region at the base of the mantle such lamellae represent in effect horizontally emplaced dikes within d excessively complex waveforms are produced when more than 20 of d volume is comprised of low velocity lamellae many lamellae models can match observed sdiff splitting 1 10 s delays of svdiff but typically underpredict scs splitting 1 4 s delays of scsv vti model parameters are selected to address d observations and include 1 0 5 to 3 anisotropy 2 discontinuous d shear velocity increases up to 3 3 d thicknesses from 100 to 300 km and 4 vti concentrated at the top bottom or throughout d vti models readily match observed splits of scs and sdiff we discuss lamellae and vti model attributes in relationship to waveform complexities splitting magnitude triplications from a high velocity d discontinuity and apparently reversed polarity svdiff onsets the possible presence of melt filled lamellae indicates that local chemical or thermal perturbations can produce regions that exceed the solidus within d such melt could occur in the bulk of d because the melt is either close to neutral buoyancy advective velocities exceed percolative velocities or both

isotopic and trace element constraints on the origin and evolution of alkaline and calc alkaline magmas in the northwestern mexican volcanic belt a wide variety of rock types has been produced by pliocene to recent volcanism in the northwestern portion of the mexican volcanic belt composite cones erupted rocks typical of calc alkaline volcanic belts associated with subduction monogenetic cinder cones surrounding the composite cones erupted mildly alkaline basalts and related rocks one larger center las navajas produced basalts trachytes and peralkaline rhyolites such alkaline rocks are typically associated with crustal extension related to rifting both subduction and rifting appear to be taking place in western mexico the most basic calc alkaline rocks although somewhat fractionated show depletion in high field strength elements hfse and enrichment in alkaline earth elements relative to other incompatible elements a characteristic typical of magmas in subduction related volcanic arcs the basic alkaline rocks show no depletion in hfse and show higher concentrations of incompatible elements than the basic calc alkaline rocks the calc alkaline rocks show a tight cluster of 87sr 86sr and 143nd 144nd ratios both of which correlate weakly with sio2 basic alkaline rocks have lower sr and higher nd isotopic ratios whereas the more differentiated alkaline rocks have isotopic ratios that overlap with and extend beyond the range of values found in calc alkaline rocks trace element data indicate that calc alkaline magmas could not have been derived by crystal fractionation from the alkaline magmas nor by assimilation of crustal materials by alkaline magmas the evolution of both calc alkaline and alkaline rocks requires a process of crystal fractionation accompanied by assimilation afc of crustal rocks to account for changes in isotopic ratios and trace element concentrations in the more differentiated members of each suite assuming a contaminant located in the lower or middle crust afc modeling shows that the amounts of assimilation required for alkaline magmas are somewhat higher than those for calc alkaline magmas calc alkaline and alkaline rocks appear to have resulted from melting of two distinct sources alkaline from an oib type source and calc alkaline from a mixed mantle and slab derived source

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

modeling anisotropy and plate driven flow in the tonga subduction zone back arc the goal of this study is to determine whether shear wave splitting observed in subduction zone back arc regions the tonga subduction zone in particular can be quantitatively modeled with flow in the back arc mantle driven by the motions of the subducting slab and the upper back arc plate we calculated two dimensional mantle flow models using known tonga plate motions as boundary conditions and assuming a range of uniform and variable viscosity structures shear wave splitting was predicted for the anisotropy due to lattice preferred orientation lpo of olivine and orthopyroxene in the flow model finite strain fields the predicted shear wave splitting provides a good match to the fast directions parallel to the azimuth of subducting plate motion and splitting times 0 5 1 5 s observed in tonga both for models where lpo anisotropy develops everywhere above 410 km and for models where lpo anisotropy is confined to regions of relatively high stress if lpo anisotropy does develop over the entire upper 410 km of the mantle the strength of anistropy induced by a given amount of shear strain must be relatively weak 4 for shear strains of 1 5 with a maximum value of 6 for very large strains the splitting observations are comparably fit by a wide range of different viscosity models anisotropy due to melt filled cracks aligned by stresses in the back arc flow models predicts fast directions roughly normal to observed values and thus cannot alone explain the observed splitting

subduction of a wedge shaped philippine sea plate beneath kanto central japan estimated from converted waves and small repeating earthquakes 1 we estimate the configuration of the philippine sea plate phs subducted beneath kanto to understand the seismotectonics in the tokyo metropolitan area in addition to small repeating earthquakes that delineate both the upper and lower boundaries of the phs previously unrecognized converted waves at its upper boundary from lower boundary repeater sources provide a new constraint on the position of the boundary the results from 794 s to p and 212 p to s converted waves together with 200 repeater locations show that 1 the phs has a wedge like shape decreasing in thickness from 50 km beneath tokyo to zero at its northeastern limit and 2 the phs is relatively flat in the eastern part of kanto and bends upward near its northeastern limit from the fact that the phs in the izu bonin forearc also has a similar wedge like shape we infer that its shape was formed before its subduction on the other hand the relatively flat phs at its northeastern limit and its upward bend are probably a consequence of ongoing deformation from an interaction with the underlying thick pacific plate pac the history of the phs as a forearc before subduction and its contact with the pac after subduction suggests that the temperature of this part of the phs is colder than other regions which probably causes anomalously deep seismicity on the phs in kanto colder phs temperatures imply that the source regions of large earthquakes may extend deeper than expected an important fact when assessing the locations of disastrous earthquakes for tokyo

seismic anisotropy beneath stable continental interiors a robust knowledge of seismic anisotropy beneath the continents is essential to our understanding of plate tectonic theory as anisotropy provides a unique constraint on the character of past and present deformation in the lithosphere and sublithospheric mantle this review paper summarizes the range of techniques currently available to image seismic anisotropy with passive source seismic data and addresses current issues surrounding the observation and interpretation of continental seismic anisotropy to this end we present case studies for four regions where seismic anisotropy has been extensively investigated in recent years eastern north america the canadian shield australia and southern africa based on this full suite of results we infer that stable continental regions generally contain seismic anisotropy that is clearly located within both the lithosphere and the sublithospheric mantle usually to depths of at least 200 km and perhaps more an implication of these results is that tectonic plates are at most only partially coupled to the underlying mantle the results from these case studies also demonstrate that while remarkable progress in seismic anisotropy imaging has been achieved in recent years it is clear that much more work will be required to adequately understand the origin of continental seismic anisotropy we suggest that a more robust characterization of anisotropic parameters can only be achieved by integrating complementary seismic datasets and by incorporating constraints from key related datasets from mineral physics magnetotellurics gravity and geodesy

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

experimental constraints on the low pressure evolution of transitional and mildly alkalic basalts multisaturated liquids and coexisting augites basaltic liquids saturated at one atmosphere in augite olivine and plagioclase show a strong linear trend in pseudo quaternary projections from nepheline to slightly quartz normative compositions coexisting augites show a corresponding variation from slightly nepheline to hypersthene normative compositions consequently alkalic and mildly alkalic liquids coexist with nepheline normative augites tholeiitic liquids coexist with hypersthene normative augites and transitional liquids coexist with augites compositionally intermediate between those crystallizing in alkalic and tholeciitic melts the variation in augite composition is due to a decrease in wollastonite and non quadrilateral components going from mildly alkalic to transitional and tholeiitic augites the observed liquid and augite trends lack correlations between normative composition and crstallization temperature this is an artifact of plotting unrelated liquid lines of descent individual liquid lines of descent in contrast display good correlation between normative composition and crystllization temperature which for both liquids and coexisting augites are highly oblique to nearly perpendicular to the general trend because the composition of the first augite to appear on the liquidus varies as a function of the bulk system the location of the thermal divide shows corresponding variations as a consequence unifying cotectics and pseudo quaternary normative projections should not without specific experimental knowledge be used to address liquid lines of descent and common magmatic parentages

seismicity and tectonics of the subducted cocos plate we have examined teleseismic earthquake locations reported by the international seismological centre isc for the middle america region and selected 220 as the most reliable these hypocenters and other data are used to delineate the deep structure of the subducted cocos plate the results indicate that the subducted plate consists of three major segments segment i extends from the panama fracture zone to the nicoya peninsula the structure of this segment is poorly defined segment ii is the largest and best defined segment this segment consists of two parts iia and iib part iia extends from the nicoya peninsula to western guatemala and is very well defined and continuous in structure its strike follows the curvature of the trench and dips at about 60 part iib extends from western guatemala to orizaba mexico the dip of this part of the segment decreases slightly toward the northwest and its strike is more northward than that of the trench segment iii extends from orizaba to the rivera fracture zone and is not well defined due to a lack of earthquake activity beneath about 100 km its orientation differs markedly from segment ii and strikes somewhat more westward than the trench between parts iia and iib of segment ii the subducted plate seems to be continuous bending smoothly to accommodate the change in geometry local network data from costa rica suggest there may be a tear between segments i and ii between segments ii and iii there is a gap in the hypocenters which makes it difficult to define the boundary the change in geometry between these two segments indicates that there may be a tear and two strike slip focal mechanisms in the region support this conclusion we find no convincing evidence supporting the existence of segments smaller than the three described above if there is smaller scale segmentation in the shallow part of the subducting plate the plate must still maintain enough continuity to appear continuous at greater depths there is no evidence for any major tear in the subducted plate associated directly with either the tehuantepec ridge or the orozco fracture zone the shallow subduction at the northwestern end of segment ii may be related to the bouyancy of the tehuantepec ridge the cocos ridge is probably directly responsible for the change in geometry between segments i and ii and may even be slowing or stopping subduction in segment i the structure of the subducted plate in segment ii and the changes in the character of volcanism along the arc can be related to the relative motion of the north american and caribbean plates the present geometry of part iib of segment ii is more consistent with the probable configuration of the trench about 7 ma ago than with the present configuration indicating that the north american plate is overriding the subduction zone appendices 2 3 and 4 are available with entire article on microfiche order from american geophysical union 2000 florida avenue n w washington dc 20009 document b84 009 2 50

three dimensional distribution of seismic anisotropy in the hikurangi subduction zone beneath the central north island new zealand 1 we use earthquake arrival time data in the central north island new zealand in inversion for 3 d vp azimuthal anisotropy it is parameterized with an isotropic component and two azimuthal anisotropy parameters for each node and local earthquake shear wave splitting observations define the initial anisotropy model our results suggest that mantle flow either within the mantle wedge or below the subducted slab is not the primary source of anisotropy in the hikurangi subduction zone the largest region of anisotropy that we image is within the subducted slab where anisotropy is consistently trench parallel and 5 9 magnitude from 32 to 185 km depth bending induced yielding in the slab offshore and as it passes through the trench provides an explanation for this slab anisotropy beneath the taupo volcanic zone tvz there is 1 4 trench normal anisotropy within 30 km of the surface of the slab within a region entrained with motion of the slab in contrast only very weak 0 2 anisotropy is imaged in the region of high temperature partial melt arising from active corner flow in the mantle wedge in the crust geologic structure is important with the largest crustal anisotropy 14 related to schist we image significant anisotropy aligned with the margins of the tvz related to extensive fracturing on the margins as the tvz actively extends within the heavily intruded and underplated tvz lower crust the orientation of anisotropy switches to align with the extension direction this orientation is consistent with other geophysical data suggesting the presence of connected melt in this region

seismic structure and crustal magmatism at the mid atlantic ridge 35 n we have determined the three dimensional p wave velocity structure and azimuthal anisotropy within the shallow crust of the inner valley at the mid atlantic ridge mar near 35 n the best fitting one dimensional model is 0 8 km s lower in velocity at a given depth than for fast spreading ridges while the isotropic three dimensional velocity structure shows a strong lateral heterogeneity up to 1 4 km s peak to peak variation at 1 km depth in both the along and cross axis directions an anomalous region with velocities as much as 0 8 km s lower than average and 3 5 km in lateral extent is located beneath a near axis seamount we interpret this feature as a region of near solidus temperatures with possibly a small degree of partial melt from the variation of travel time residuals with azimuth we determine a best fitting p wave anisotropy model that has 4 anisotropy from 500 m to 1 km depth 2 anisotropy in the depth range 1 1 5 km and no anisotropy at greater depth this finding is consistent with vertical cracks aligned parallel to the axial valley we find no evidence for the shallow high velocity ridge parallel feature commonly observed at the axis of fast spreading ridges and attributed to high velocity dikes beneath a thin extrusive layer the absence of this anomaly at the mar implies that the axis of crustal accretion is not at a stable location with respect to the inner valley walls on timescales greater than about 25 000 years an unstable axis of accretion will result in a thicker transition from extrusive rocks to dikes which may account for the generally lower velocities beneath the mar compared with the east pacific rise and a lateral variation in upper crustal lithology

tears or thinning subduction structures in the pacific plate beneath the japanese islands abstract the nature of a subduction zone at depth is affected by the evolution of its tectonic system and the geometry of the trench line can change over time due to slab roll back or the arrival of a distinctive feature with the incoming oceanic lithosphere the configuration of the plate has to accommodate such changes with buckling thinning or the formation of tears depending on the rate of influx to the trench tomographic imaging is commonly used to recognise the presence of such tears through marked reductions in wavespeed anomalies in localised zones a good example is provided by pacific plate subduction beneath the japanese islands a horizontal tear in the plate below 300 km depth can be recognised at the southern end of the izu bonin arc associated with the change in slab morphology to the much steeper mariana arc beneath southern honshu a break in the fast wavespeeds associated with the pacific plate has been described as a tear based on the evidence of converted phases from the edge of the zone and tensional focal mechanisms for seismic events in the tear zone in the north close to the hokkaido bend in the subduction zone the reduction in the shear wavespeed anomaly is just as dramatic but here the characteristics of high frequency guided waves from deep earthquakes indicate continuity of slab material with thinning of the slab the thinned slab has less wavespeed contrast within the affected cells and so appears in the tomographic images as a weakened anomaly the various modes of slab deformation represent different ways in which the subducted material accommodates the strains imposed by the evolution of the geometry of the subduction scenario not all significant reductions in wavespeed anomalies represent tears and thus it is important that such interpretations be checked against the characteristics of wave propagation through the zone

deep structure of japan subduction zone as derived from local regional and teleseismic events we have determined a detailed three dimensional p wave velocity structure of the japan subduction zone to 500 km depth by inverting local regional and teleseismic data simultaneously we used 45 318 p wave arrivals from 1241 shallow and deep earthquakes which occurred in and around the japan islands the arrival times are recorded by the japan university seismic network which covers the entire japan islands densely and uniformly we also used 4211 travel time residuals from 100 teleseismic events which are read from seismograms recorded by seismic stations in northeastern japan in comparison with the previous results obtained from only local and regional events the present result for the area around the lower plate boundary and the mantle below the plate is determined more reliably because of the addition of 7035 data from 100 teleseismic events and 41 very deep earthquakes in the crust and uppermost mantle low velocity zones are clearly visible beneath active volcanoes in the mantle wedge the low velocity zones generally parallel with the slab and exist continuously to a depth of about 200 km which is consistent with the petrological geochemical and geodynamic studies we consider that the existence of volcanism related low velocity anomalies in the mantle wedge is a general seismological characteristic of subduction zones in light of all the available tomographic results for many subduction zones in the world the pacific slab beneath japan is imaged more clearly than in previous studies as a high velocity zone with a thickness of 80 90 km and a p wave velocity 4 6 higher than the normal mantle lower velocity anomalies are found in the mantle below the slab

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

shear wave splitting beneath the southwestern kurile arc and northeastern japan arc a new insight into mantle return flow 1 shear wave splitting in the southwestern part of the kurile arc and the northeastern ne japan arc is investigated using the waveforms from local earthquakes for both arcs observed shear wave splitting shows clear evidence for a striking rotation of fast direction across the arc suggesting the different feature of anisotropy between the fore arc and back arc sides trench parallel fast directions are observed in the fore arc side which is consistent with the anisotropy expected from the deformation of the b type olivine fabric as well as the anisotropy in the slab fast directions observed in the back arc side show approximately e w and ese wnw in the ne japan arc and n s in the southwestern kurile arc which are characterized by the local maximum dip direction of the subducted pacific plate not by the direction of plate convergence this relationship suggests that the direction of the mantle return flow is governed by the local slab geometry

anisotropic pn tomography in the western united states pn travel times are affected not only by lateral variations in crust and mantle velocity but also by significant amounts of laterally varying anisotropy to investigate uppermost mantle anisotropy a tomography algorithm was reformulated to include lateral variations in both velocity and horizontal anisotropy and it was applied to pn travel time data from the western united states results show that anisotropy is as important in explaining the travel time residuals as are the velocity variations a detailed resolution study examined the trade off between the velocity variations and the anisotropy variations and showed that both could be resolved for regions with good ray path coverage pn anisotropy beneath the western united states has maximum amplitudes of 0 3 km s 7 5 when resolved on a length scales of around 3 the orientations of pn anisotropy often correlate well with those inferred from shear wave splitting results this correlation suggests that orientation of mantle anisotropy does not change significantly with depth for many regions the orientation of the pn anisotropy can be correlated with some of the tectonic processes which occur within the western united states for example the fast direction of pn anisotropy parallels the northeast direction of subduction of the juan de fuca plate beneath the northwest pacific coast suggesting that there the anisotropy is related to subduction driven deformation the fast direction of pn anisotropy also parallels the strike of the san andreas fault system in central california indicating that shear strains along the plate boundary may be responsible for the anisotropy there within the great basin the pn anisotropy varies substantially and both partial melting and small scale convection within the uppermost mantle could be responsible for these anisotropy variations

subduction of the philippine sea plate beneath southwestern japan slab geometry and its relationship to arc magmatism 1 we carry out high resolution three dimensional seismic tomography of the crust and upper mantle beneath southwestern sw japan using arrival time data obtained from the nationwide seismograph network the tomographic images provide new insights into the configuration of the philippine sea slab and arc magmatism the results confirm the existence of an aseismic portion of the philippine sea slab at greater depths beyond the seismic portion the philippine sea slab is subducting aseismically down to at least a depth of 200 km in chubu 60 80 km in kinki and 60 km in chugoku and is subducting seismically to depths of 150 200 km in kyushu in the chubu district it is subducting seismically at a shallow dip and then bends downward beneath the volcanic area whereas in the kinki district it is subducting subhorizontally as far as the japan sea the differences in the slab geometry between adjacent regions correlate with the differences in volcanic activity suggesting the importance of the geometry of the philippine sea slab on arc magmatism our tomographic images further imply subduction of the philippine sea slab down to a depth of around 70 km northwest of the izu peninsula the philippine sea slab in that region might have split into western and eastern parts separated by a slab tear and collided with the pacific slab at depths of 150 200 km

constraints on the origin of alkaline and calc alkaline magmas from the tuxtla volcanic field veracruz mexico lavas erupted in the tuxtla volcanic field tvf over the last 7 ma include primitive basanites and alkali basalts mildly alkaline hy normative mugearites and benmoreites and calc alkaline basalts and basaltic andesites the primitive lavas are silica undersaturated with high concentrations of both incompatible and compatible trace elements variable la yb with constant yb at 6 to 8 times chondritic and low sr and o and variable pb and nd isotopic ratios the primitive magmas originated by increasing degrees of melting with pressure decreasing from greater than 30 kbar to 20 kbar in the garnet stability field another group of alkali basalts and hawaiites has lower ni and cr concentrations and higher fe mg ratios and was derived from the primitive group by crystal fractionation at pressures of several kbar incompatible trace elements in these silica undersaturated lavas show depletion in high field strength elements hfse relative to large ion lithophile elements similar to subduction related basalts ba nb ratios are nearly constant and thus the hfse depletion cannot be the result of a residual hfse bearing phase in the source but could be the result of generation from a source contaminated by fluids or melts from the subducted lithosphere the silica saturated mugearites and benmoreites and the calc alkaline basalts and basaltic andesites were erupted only between 3 3 and 1 0 ma these have incompatible element concentrations generally lower than in the silica undersaturated lavas and thus could not have been derived by crystal fractionation from the silica undersaturated alkaline magmas magmas parental to the silica saturated magmas originated by higher degrees of melting at lower pressures than the primitive magmas melting may have been promoted by an influx of fluid from the subducted lithosphere trace element and sr nd pb and o isotopic data suggest that three components are involved in the generation of tvf magmas the mantle a fluid from the subducted lithosphere and continental crust tvf alkaline lavas are similar to those erupted in the back arc region of the mvb and japan and show characteristics similar to alkaline magmas erupted in the southern andean volcanic arc these low degree melts reach the surface along with calc alkaline lavas in the tvf due to an extensional stress field that allows their passage to the surface

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

geometry and relative motion of the philippine sea plate and pacific plate beneath the kanto tokai district japan the recent acquisition of high density and improved seismic data by the kanto tokai k t observational network of the national research institute of earth science and disaster prevention requires the revision of conventional plate configuration models for the k t district we propose a new interpretation of the configuration and relative motion of the philippine sea phs plate pacific pac plate and eurasian eur plate on the basis of the distribution of hypocenters velocity structure and focal mechanisms in the k t district the new model for the configuration of the phs and pac slabs clearly delineates the phs slab subducting northwestward beneath the k t district from the sagami and suruga troughs and the pac slab subducting westward beneath the phs slab from the japan trench the thickness of the descending phs slab is estimated to be 30 5 km the proposed configuration of the phs slab is significantly different from those obtained by other researchers including our own previous study the new model is fully consistent with global plate motion data the new model leads to an interpretation of the tectonic processes taking place beneath the k t district for the sagami trough wing of the phs slab the northeastern portion of the slab is dipping northwestward in direct contact with the upper surface of the pac slab the southwestern portion carrying the izu peninsula has been bent upward because of its buoyancy producing the observed contorted shape for the suruga trough wing the dip direction of the phs slab gradually changes from northwestward at the eastern portion to northward at the western portion the new model also indicates the presence of offset within the phs slab between the sagami and suruga trough wings this area is characterized by a low level of seismic activity the absence of a high velocity structure and the presence of quaternary volcanoes it seems reasonable to speculate that the distribution of quaternary volcanic edifice construction is controlled by the location of the offset within the phs slab in view of this new model the tectonic processes of large earthquakes which occurred recently in the studied area are also reexamined

deep structure of the japan subduction zone abstract we determined 3 d p wave velocity structure down to 700 km depth under the japan islands using a large number of arrival time data from local and teleseismic events simultaneously we collected 207 000 arrival times from 7743 shallow and deep earthquakes occurred in and around japan and 34 148 data from 333 teleseismic events which were recorded by over 1000 seismic stations on the japan islands our tomographic model revealed some new features the philippine sea slab is found to subduct down to 500 km depth under southwest japan though the seismicity within the slab ends at 150 200 km depth significant low velocity anomalies are found to exist in the deep portion of the mantle wedge above the pacific slab which may be caused by the deep dehydration process of the slab slow anomalies are detected in the mantle beneath the pacific slab which may be caused by a mantle plume or upwelling of hot mantle materials associated with the deep subduction of the pacific slab and its collapsing down to the lower mantle

primitive calc alkaline and alkaline rock types from the western mexican volcanic belt since the early pliocene small volumes of alkaline magmas have erupted in the western mexican volcanic belt mvb in close association with the volumetrically dominant calc alkaline magmas both suites include relatively rare primitive types characterized by high mg cr and ni contents primitive hypersthene normative basalts parental to the calc alkaline suite are found along the volcanic front but are absent at greater distances from the trench the volcanic front calc alkaline suites are typically hornblende bearing indicating relatively high water contents associated primitive alkaline magmas at volcanic front locations also contain hydrous minerals these nepheline normative suites include basanites phlogopite bearing minettes and other hornblende bearing lamprophyres these are probably the youngest and freshest lamprophyres yet discovered on earth and ideal samples for addressing the origin of this exotic class of rocks on extended morb normalized elemental plots all lamprophyres show patterns similar to the calc alkaline basalts but have overall enrichments of 2 to 25x in ti k p ba sr light rare earth elements and other incompatible elements the lamprophyres show the same strong relative enrichments of ba k and sr and depletions in ti and nb that characterize the calc alkaline rocks of western mexico and all subduction zones these similarities in relative element abundances support a common source region for all primitive magmas of the western mvb the primitive calc alkaline and alkaline rocks of western mexico show narrow overlapping ranges in sr nd and pb isotopic ratios also consistent with derivation of all magmas from a common source the presence of the relatively cold subducted cocos and rivera plates below the volcanic front demands that all primitive magmas originated in the depth range 30 to 75 km these magmas are enriched in ba k sr and other elements probably transported by hydrous fluids rising from the subducted slab several lines of evidence support the presence of phlogopite in the source region of western mvb magmas probably in the form of phlogopite rich veins cutting asthenospheric peridotite the first liquids formed upon partial melting of this veined source region will concentrate the vein component and generate the lamprophyres these lamprophyric melts rarely erupt in arc settings more commonly forming dikes in the arc crust eruption of lamprophyric magmas in the western mvb is favored by through crustal extensional fracture systems related to active rifting of the jalisco block from the n american plate larger volumes of calc alkaline basalt are generated by melting of the same source which dilutes the vein component with larger proportions of peridotitic wall rock

correction to isotopic and trace element constraints on the origin and evolution of alkaline and calc alkaline magmas in the northwestern mexican volcanic belt by surendra p verma and stephen a nelson

seismic anisotropy of the pacific slab and mantle wedge beneath the japanese islands 1 the mantle beneath the japanese islands is complex because of subduction of the pacific and philippine sea plates and the deformation associated with it detection of seismic anisotropy should be useful for understanding the processes occurring in such mantle here we resolve seismic anisotropy of the mantle wedge from that of the underlying subducted pacific slab using a large number of measurements of shear wave splitting for a family of core reflected shear phases scs sscs scs2 and sscs2 the anisotropy of the mantle wedge changes sharply across the volcanic front on the pacific side of the volcanic front a vertically propagating shear wave is polarized with the fast direction approximately parallel to the trench whereas on the back arc side it is polarized with the fast direction approximately parallel to the plate convergence direction the pacific slab is uniformly anisotropic with the nnw fast direction consistent with the paleospreading direction of the northwestern pacific seafloor the preferential alignment of anisotropic crystals frozen in the plate at the time of its formation appears to be preserved down to depths beyond 400 km within the slab

p wave anisotropic tomography beneath northeast japan we determined a 3 d p wave anisotropic tomography beneath northeast japan by using first p arrival times the travel time inversion is parameterized with an isotropic component and two anisotropic parameters for each grid node by assuming p wave azimuthal anisotropy with hexagonal symmetry axis distributed horizontally the geometries of the conrad and moho discontinuities and the upper boundary of the subducting pacific plate are taken into account in the model our extensive synthetic tests show that the anisotropic velocity inversion method is reliable our results show that low velocity zones exist below active volcanoes in the crust in the fore arc mantle wedge and in the central portion of the mantle wedge above the high velocity pacific slab which are similar to the previous tomographic images anisotropic velocity variations are revealed in the crust mantle wedge and the subducting pacific slab the anisotropic patterns are complex and different in the upper and lower crust in the mantle wedge the fast direction of p velocity is generally trench parallel in the fore arc area while it becomes trench normal in the back arc area which agree well with the shear wave splitting results the trench parallel fast direction is revealed beneath the volcanic front which may suggest the existence of complex 3 d mantle flows in the mantle wedge the fast velocity direction in the subducting pacific slab is mostly n s which is likely due to the original fossil anisotropy of the pacific plate formed at the mid ocean ridge our result also suggests that the anisotropy is weak in the upper crust and the average amplitude of the velocity anisotropy in the mantle is no larger than 4

novel anisotropic teleseismic body wave tomography code anitomo to illuminate heterogeneous anisotropic upper mantle part i theory and inversion tuning with realistic synthetic data

mapping p wave anisotropy of the honshu arc from japan trench to the back arc abstract we determined a 3 d p wave anisotropic tomography beneath the entire honshu arc using about 448 000 high quality p wave arrival times from 18 335 local earthquakes that occurred beneath the northeast ne japan land area and the fore arc area under the pacific ocean our results show that low velocity low v zones exist beneath the active arc volcanoes in the crust and in the central portion of the mantle wedge above the subducting pacific slab the low v anomalies are related to the arc magmatism low v zones are also revealed in the fore arc area which are probably caused by large volumes of water releasing upwards from dehydration of the subducting oceanic crust and sediments the anisotropic amplitude in the upper crust is weaker than that in other portions under ne japan in the mantle wedge the fast velocity direction fvd is generally trench normal in back arc area which may reflect that the olivine a axis aligns with the transport direction induced by the slab driven corner flow the fvd becomes trench parallel in the central portion of the fore arc mantle wedge which is possibly induced by the olivine of b type fabric in the slab driven corner flow the fvd shows trench parallel in the low v zones in the fore arc mantle wedge close to the upper boundary of the pacific slab which may reflect the b type olivine fabric dominating in those areas or it may be induced by the dextral shearing of the overlying crust the trench parallel fvd is also revealed beneath the volcanic front indicating complex 3 d mantle flows in the mantle wedge the fvd in the subducting pacific slab is mostly trench parallel which may reflect one or more of the following possibilities 1 the original fossil anisotropy of the pacific plate formed at the mid ocean ridge 2 cracks within the slab and 3 the olivine fabric transition due to the changes in water content stress and temperature

seismic wave conversion near the upper boundary of the pacific plate beneath the kanto district japan summary a detailed structure of the pacific plate beneath the kanto district japan was studied by the use of conspicuous sp converted waves that are generated near the upper boundary of the pacific plate these sp waves are observed for earthquakes occurring in the lower plane of the double seismic zone of the subducted pacific plate the location of the conversion interface is estimated by inverting the observed tsp p times 293 seismic traces recorded at 42 stations for 49 earthquakes are used for analysis the conversion interface is located above the upper seismic plane of the double seismic zone and is equivalent to the upper boundary of the pacific plate within the estimated uncertainty in general the conversion interface coincides well with the zone of low angle thrust fault type earthquakes in the entire area studied comparison with the result of 3 d tomographic inversion reveals that the conversion interface is probably located in the low velocity layer on the top of the pacific plate which is interpreted as the oceanic crust of the pacific plate

anomalous low velocity zone and linear alignment of seismicity along it in the subducted pacific slab beneath kanto japan reactivation of subducted fracture zone 1 a detailed investigation of the hypocenter distribution beneath kanto japan reveals a nw se trending linear alignment of seismicity within the subducted pacific slab we estimate the 3d seismic velocity structure in the pacific slab to understand the factors controlling the genesis of such intraslab earthquakes a narrow low velocity zone is imaged within the subducted slab over a length of 150 km which partly penetrates into the mantle portion of the slab the low velocity zone correlates in space with the nw se trending earthquake cluster a reactivation of hydrated fracture zone formed prior to subduction is probably related to the low velocity anomaly dehydration reactions of the hydrated oceanic mantle as well as the oceanic crust might lower the seismic velocity along the fossil fracture zone accompanied by intraslab earthquakes these observations support the hypothesis of dehydration embrittlement as the most viable mechanism for generating intraslab earthquakes

exhumed upper mantle within the ocean continent transition on the northern west iberia margin evidence from prestack depth migration and total tectonic subsidence analyses 1 the west iberia continental margin contains an enigmatic north south ridge of serpentinized peridotite located within the ocean continent transition we used prestack depth migration and total tectonic subsidence analysis to constrain the basement type adjacent to the ridge in order to learn more about its emplacement prestack depth migration was performed on segments of the nine east west lines from the iberia seismic experiment ise 97 and total tectonic subsidence analyses were performed on the depth sections the seismic data and subsidence analyses are consistent with exhumed upper mantle occurring immediately landward of the ridge suggesting that the ridge is located within a broad zone of exhumed mantle that has been serpentinized the subsidence analyses also provide minimum bounds on the thickness of the serpentinized layer which extends from the zone of exhumed mantle landward under thinned continental crust exhumed upper mantle has previously been identified in the wider ocean continent transition to the south in the southern iberia abyssal plain and we conclude that similar but narrower exposures occur to the north in the galicia bank area while upper mantle appears to have been tectonically exhumed along all the ise 97 lines it does not always form a distinct basement ridge there is no well defined ridge on three seismic profiles in the middle of the ise 97 survey where it is well developed the peridotite ridge parallels a deeply penetrating west dipping normal fault and the three middle lines appear to represent a transfer zone where this normal fault is not well developed