mantle phase changes and deep earthquake faulting in subducting lithosphere inclined zones of earthquakes are the primary expression of lithosphere subduction a distinct deep population of subduction zone earthquakes occurs at depths of 350 to 690 kilometers at those depths ordinary brittle fracture and frictional sliding the faulting processes of shallow earthquakes are not expected a fresh understanding of these deep earthquakes comes from developments in several areas of experimental and theoretical geophysics including the discovery and characterization of transformational faulting a shear instability connected with localized phase transformations under nonhydrostatic stress these developments support the hypothesis that deep earthquakes represent transformational faulting in a wedge of olivine rich peridotite that is likely to persist metastably in coldest plate interiors to depths as great as 690 km predictions based on this deep structure of mantle phase changes are consistent with the global depth distribution of deep earthquakes the maximum depths of earthquakes in individual subductions zones and key source characteristics of deep events

a new self organizing mechanism for deep focus earthquakes the mechanism of deep focus earthquakes has been a puzzle since their discovery almost 70 years ago1 2 because brittle fracture and frictional sliding at depths in excess of 100 200 km would require unrealistic rock strengths3 4 rock strength does increase with pressure but a few hundred mpa equivalent to 10 20 km depth suffices to inhibit most fracture and elevated temperature activates ductile mechanisms that operate at stresses less than the fracture stength a range of mechanisms has been proposed for deep earthquakes including plastic instabilities5 7 shear induced melting8 11 and instabilities accompanying recrystallization12 13 or polymorphic phase transformation14 23 each of these proposed mechanisms has exhibited certain inherent weaknesses see kirby22 for review and discussion here we report experimental observations of high pressure faulting of metastable mg2geo4 olivine as it undergoes incipient transformation to a spinel 24 25 structure we present a model in which this faulting arises from the generation propagation and linking up of spinel filled anticracks26 when applied to the olivine spinel transformation in the earth s mantle the anticrack model27 satisfactorily accounts for the similarities and differences between shallow and deep earthquakes and removes the problems associated with frictional sliding

a view of the lower crustal component of hydrothermal systems at the mid atlantic ridge an extensive suite of hydrothermally altered gabbros was recovered by alvin and dredging from the mid atlantic ridge in the vicinity of the eastern ridge transform intersection rti of the kane fracture zone mark where tectonic extension has provided a window into the lower crustal component of hydrothermal cells four alteration types are distinguished on the basis of metamorphic assemblage mineral composition and deformation textures a conceptual model is presented that places temperature fluid rock ratio relationships in conjunction with styles of deformation into a spatial and temporal framework exsolution of late stage magmatic fluids in the vicinity of the magma hydrothermal interface at temperatures 700 c marks the onset of fluid rock interaction in the mark gabbros kelley et al this issue with subsequent alteration involving seawater derived fluids being strongly influenced by deformation mechanisms in zones of ductile shear hydration was initiated at temperatures between 550 and 700 c and moderate fluid rock ratios elsewhere in the lower crust the onset of seawater penetration occurred at temperatures of 400 550 c and was facilitated by brittle fracturing early vein networks indicate very low initial fluid rock ratios and are interpreted as being related to the solidification of plutons later fractures indicate higher water rock ratios and may be related to the downward propagation of shallow fault systems throughout the crust fluid rock interaction ceased at temperatures 180 300 c subsequent cataclastic deformation associated with the unroofing and emplacement of crustal blocks in the rti massif produced moderate fluid rock ratios within localized zones hydrothermal alteration in the crustal section exposed at mark was initiated at lower temperatures and water rock interaction proceeded to lower temperatures than many other plutonic suites these differences must reflect the temporal and spatial evolution of magmatic and tectonic extension that is particular to each section of sampled crust

stress distribution due to olivine spinel phase transition in descending plate and deep focus earthquakes the stress field due to olivine spinel phase transition in and around a descending plate is studied based on the calculations for an equilibrium phase transition with a constant clapeyron s slope and for a nonequilibrium transition taking a kinetic effect into consideration the material concerned is assumed to be a viscoelastic body whose physical coefficients are functions of temperature and pressures the calculated stress distribution for the equilibrium transition is characterized by a compressional stress near the upper and or lower surfaces of the plate and by a tensional stress at the central part in the depth range from 200 to 550 km the magnitude of principal stress being greater than 0 5 gpa however for the nonequilibrium transition a tensional stress predominates near the upper and lower surfaces and the central part is compressional in the depth range from 300 to 600 km unlike the case of equilibrium transition the maximum principal stress in the latter case is greater than 2 gpa which is much more than that for the equilibrium transition the principal axes in the high stressed region are oriented almost parallel to the descending direction in both cases comparison of the calculated results to the actually observed seismicity depth relation and to focal mechanisms of deep focus earthquakes in descending plates beneath island arcs shows that the nonequilibrium phase transition agrees better with the observations

double seismic zone for deep earthquakes in the izu bonin subduction zone a double seismic zone for deep earthquakes was found in the izu bonin region an analysis of sp converted phases confirms that the deep seismic zone consists of two layers separated by 20 kilometers numerical modeling of the thermal structure implies that the hypocenters are located along isotherms of 500 to 550 c which is consistent with the hypothesis that deep earthquakes result from the phase transition of metastable olivine to a high pressure phase in the subducting slab

phase transformations earthquakes and the descending lithosphere abstract the energy dissipated from metastable phase transformations is investigated in this paper in addition to the kinetic energy envisaged by randall at least two other major forms of energy latent heat and metastability energy are involved in a metastable phase transformation the seismic energy density for metastable phase transformation is of the order of 10 24 10 25 erg km 3 for most mantle silicates as an example an earthquake of magnitude m 7 5 would require only 0 01 0 03 or 0 04 km 3 of basalt olivine or pyroxene respectively undergoing metastable phase transformation earthquakes caused by both faulting and phase transformation mechanisms have been observed in high pressure laboratory experiments the energy density associated with a faulting mechanism appears to be four to five orders of magnitude less than that associated with metastable phase transformations of mantle silicates confirming theoretical expectations according to all of the thermal models for the descending lithosphere available in the literature metastable phase transformations almost certainly occur in the descending lithosphere enormous energies would be released suddenly in association with these metastable phase transformations thus earthquakes within the descending lithosphere at depths greater than 100 km may be induced by metastable phase transformations however metastable phase transformations in the descending lithosphere have little or an opposing effect on the tectonic engine driving the motion of plates into the mantle

downdip tensional earthquakes beneath the tonga arc a double seismic zone characteristics of seismicity at intermediate depth beneath the tonga arc are studied special attention is paid to newly identified downdip tensional events because the dominance of downdip compression has been known in this area we reexamine the relatively shallow 70 km downdip tensional events reported in the harvard centroid moment tensor solutions by modeling bodywave waveforms of world wide standard seismograph network seismograms and conclude that they are in fact intraplate downdip tensional events furthermore these downdip tensional events are consistently located below and seaward of downdip compressional events and appear to constitute a double seismic zone the result of a relative relocation technique applied to these earthquakes also supports this observation which was originally based on routinely determined epicenters the subducting slab beneath the tonga arc has been considered as a prototype of a compressionally loaded slab the existence of a double seismic zone in tonga therefore would suggest that double seismic zones located at intermediate depth in subducting slabs are more common than previously thought

significance of non double couple components of deep and intermediate depth earthquakes implications from moment tensor inversions of long period seismic waves abstract analysis of long period seismic waves suggests that non double couple components in source moment tensors of intermediate depth and deep earthquakes are significant and appear to respond to the state of predominant strain release within slabs the strain regime is partially or fully induced by the sources themselves or by slab structures near the sources we observe consistency in the non double couple components from three different seismic wave inversions among 21 earthquakes studied 17 events have the same sign in the three non double couple components and the signs show a correlation with the state of strain release within the slabs the effects of unmodeled propagation structure and instability of the inversion procedure are unlikely to be responsible for the consistent non double couple components because the suites of seismic waves traverse different paths and the various inversion schemes have the different resolvabilities of the moment tensors

waveform effects of a metastable olivine tongue in subducting slabs we constructed velocity models of subducting slabs with a kinetically depressed olivine and spinel transition and examined the effect that such structures would have on teleseismic p waveforms using a full wave finite difference method these two dimensional calculations yielded waveforms at a range of distances in the downdip direction the slab models included a wedge shaped low velocity metastable olivine tongue moto to a depth of 670 km as well as a plausible thermal anomaly one model further included a 10 km thick fast layer on the surface of the slab the principal effect of moto is to produce grazing reflections at wide angles off the phase boundary generating a secondary arrival 0 to 4 seconds after the initial arrival depending on the take off angle the amplitude and timing of this feature vary with the lateral location of the seismic source within the slab cross section careful analysis of waveforms from earthquakes with depths near 400 km simple sources and adequate station coverage in appropriate geometries will be required to resolve whether moto is present

localized polymorphic phase transformations in high pressure faults and applications to the physical mechanism of deep earthquakes earthquake mechanisms based on frictional instabilities are widely accepted for relatively shallow earthquakes such mechanisms in unmodified form are highly unlikely for deep earthquakes because of the expected high frictional shear stresses at depths of 300 700 km several alternative mechanisms have been proposed here i discuss a new type of shear instability observed in experiments on ordinary polycrystalline ice and polycrystalline tremolite by my colleagues and me these shear instabilities lead to rapid stress drops and sound emissions and have the following characteristics 1 a single sharp shear zone forms in the maximum shear stress orientation 2 confining pressure does not suppress the instability and under some conditions raising confining pressure pc lowers the differential stress at which it occurs 3 stress drops can occur repeatedly after reloading at about the same differential stress 4 the mean stress vs temperature conditions at which the instabilities occur in ice and tremolite place them in the fields of thermodynamic stability of ice ii and diopside talc respectively although i have not as yet confirmed the existence of these high pressure phases in the shear zones transformation in the shear zone is implicated because such shear instabilities do not occur in similar materials in which polymorphic phase changes do not occur reconstructive polymorphic phase changes localized along shear zones may be a general phenomenon in crystalline materials stressed non hydrostatically near transition boundaries moreover by analogy the reconstructive phase transformations that are thought to occur in deeply subducted lithosphere may lead to localized shear instabilities manifested by deep earthquakes and produce the observed double couple first motion distribution finally if this is the earthquake source mechanism in deeply subducted lithosphere then earthquakes should shut off when all of the polymorphic phase changes have occurred this may explain the maximum depth of earthquakes at about 680 km

kinetics of the olivine spinel transformation in subducting lithosphere experimental constraints and implications for deep slab processes abstract the persistence of metastable olivine to depths greater than 400 km in subducting slabs has implications for the generation of deep focus earthquakes the magnitude of buoyancy forces driving plate motion and the state of stress in the slab the depth to which metastable olivine can survive in a subduction zone and the depth interval over which transformation to or mg fe 2sio4 occurs have been evaluated from experimental kinetic data for the mg2geo4 and ni2sio4 transformations and the mg2sio4 transformation the data were extrapolated using a kinetic model for grain boundary nucleation and interface controlled growth under non isothermal and non isobaric conditions the results predict that metastable mg1 8fe0 2sio4 olivine survives to depths greater than 500 km in the cold interior of rapidly subducting slabs of old lithosphere the onset of transformation to mg fe 2sio4 spinel depends only on growth kinetics and coincides with the 550 50 c isotherm including the effects of latent heat production causes the transformation to occur by a runaway process over a very narrow depth interval at the onset of transformation high nucleation rates and low growth rates are consistent with the formation of very fine grained reaction products which are required for the transformational faulting mechanism of deep focus earthquakes when olivine crosses the equilibrium boundaries at higher temperatures e g higher than 700 c at 400 km depth transformation to or occurs much closer to equilibrium at a depth controlled by nucleation kinetics in this case the effect of latent heat production on the transformation kinetics is small and microstructural evolution is unlikely to result in transformational faulting below the depth of olivine breakdown cold slabs are likely to have a complex rheological structure owing to temperature dependent microstructural evolution during the phase transformation

seismicity shear failure and modes of deformation in deep subduction zones abstract planes of shear failure in deep seismic zones are identified as being parallel to focal mechanism nodal planes of deep focus earthquakes we relocate deep seismicity reported in the international seismological centre isc catalogue for earthquakes deeper than 400 km in the honshu bonin marianas java banda and south america subduction zones using the joint hypocentral determination method recent events with harvard centroid moment tensor cmt solutions and nearby seismicity are viewed from the direction of the cmt null axis comparing the two orthogonal nodal planes with the seismicity pattern planar features of seismicity parallel to cmt nodal planes are found in each deep seismic zone and are identified as planes of shear failure the sense of displacement on these planes is one of resistance to deeper penetration honshu and bonin show imbrication and horizontal displacement of the deepest part of the seismogenic slab while marianas and brazil show conjugate shear failure banda and argentina show patterns of vertical and horizontal shear planes for the western pacific slabs we find significant lateral displacements towards the south of the deepest part of the subduction flow tonga has been analyzed previously showing the same pattern as identified here in all benioff zones planes of shear failure are observed which act to shorten the slabs and displace matter horizontally normal to the slab moreover differential faulting acts to rotate the slab of subducted lithosphere to a more horizontal posture indicating that slabs experience resistance to penetration of the 670 km discontinuity

seismological evidence for metastable olivine inside a subducting slab the nature and location of the olivine spinel phase transition inside subducting slabs differ greatly from the situation in the surrounding mantle due to the different temperature distribution inside the slabs two models have been proposed for this phase transition one in which the location of the phase boundary between olivine and modified phase spinel is determined by equilibrium thermodynamics1 and the other including a metastable olivine phase which persists to a depth of 550 km ref 2 the location of the olivine spinel transition in the slab may be relevant to the generation of deep earthquakes3 5 and to the buoyancy forces driving subduction6 here we use travel time residuals from deep earthquakes recorded by the dense seismograph network in japan to investigate the configuration of the olivine spinel phase boundary inside the subducting pacific plate theoretical travel time residuals for the equilibrium model do not fit the observed residuals whereas those for the metastable model do implying the presence of metastable olivine inside the subducting slab

rheological control of wadati benioff zone seismicity intermediate and deep focus earthquakes in wadati benioff zones are thought to occur in the cold interiors of downgoing slabs which are significantly stronger than the warmer mantle given that earthquakes in oceanic lithosphere appear restricted by an isotherm and hence a given value of lithospheric strength we investigate whether a similar formulation is useful for subducting plates strength in downgoing slabs should be affected by both pressure and temperature an effect previously treated using a depth dependent limiting temperature for seismicity wortel 1982 we find this limiting temperature implies that a possible limiting strength increases strongly with depth unless either the temperatures were too low or the activation volume too large comparison of the analytic model used by wortel with numerical thermal models appears to exclude the first possibility we explore the second possibility by using the numerical thermal model to compute strength contours for flow law constants reported from laboratory experiments and find that the expected pressure strengthening is large enough that the slab should have considerable strength well below the deepest seismicity we conclude that if laboratory results are applicable to these conditions either a strongly depth dependent limiting strength exists or factors in addition to strength control the distribution of subduction zone earthquakes

the rheology of faults triggered by the olivine spinel transformation in mg2geo4 and its implications for the mechanism of deep focus earthquakes the relevance of the anticrack faulting mechanism to the origin of deep focus earthquakes depends critically on the physical processes operating during fault generation and movement we show here that the resistance to sliding on fault zones produced by this mechanism in mg2geo4 depends only weakly on confining pressure but strongly on sliding rate both of these characteristics are inconsistent with friction the microstructures of the fault zones and the crystal plastic rheology of the constituent phases indicate that sliding must occur primarily by grain boundary sliding previously we reported that anticrack faulting is accompanied by elastic radiation of energy acoustic emissions and that it can operate during the transformation in mg fe 2sio4 at the pressures and temperatures at which deep earthquakes occur the present results and the recent demonstration that metastable olivine is present in the subducting slab beneath japan provide additional support for the anticrack theory of deep focus earthquakes

relationship of deep seismicity to the thermal structure of subducted lithosphere recent experimental work on silicate olivine polymorphs1 has confirmed earlier observations of mechanical failure in analogous compounds2 3 as a result of phase transformations when non hydrostatic stresses were applied to a metastable phase this transformational faulting 4 mechanism is considered to be a leading candidate5 among others involving olivine transformation6 7 for the cause of deep earthquakes evidence consistent with this mechanism comes from the observation2 4 that worldwide earthquakes become more numerous with increasing depth after a seismicity minimum at about 350 km depth8 this depth lies within the range anticipated for mineralogical transformations in the subducted lithosphere9 but individual subduction zones differ in their thermal structures so if transformational faulting indeed contributes to deep seismicity this should be reflected in the location of the seismicity minimum for each zone the depth at which the minimum occurs should depend on temperature in the same way as do the polymorphic transformations of olivine and should always lie deeper than the depth at which the equilibrium transformation takes place here we test these predictions for eight subduction zones worldwide we find that with one exception the north japan zone the depth of the seismicity minimum decreases linearly with increasing thermal age of the slab a measure of its temperature profile these results support the proposal that deep earthquakes are a consequence of phase transformations in olivine

evidence for transformational faulting from a deep double seismic zone in tonga double seismic zones planes or earthquakes parallel to the dip of a subducting slab and separated by 20 40 km provide important clues about the earthquake generating mechanisms and strain distribution inside subducting slabs double seismic zones have been found at intermediate depths 70 200 km in many subduction zones1 6 but have not been previously reported in deep slabs here by relocating earthquakes with a hypocentroidal decomposition technique7 and visualizing the earthquake positions and uncertainties in three dimensions we identify a double seismic zone at depths of 350 460 km in the tonga subduction zone source parameters of the earthquakes determined by waveform analysis suggest different stress orientations for the two zones with in plane compression in the lower zone and in plane tension in the upper zone the double zone may be due to transformational faulting as olivine along the edges of a metastable olivine wedge becomes warmer and transforms to spinel8 11

first measurements of the ionospheric plasma escape from mars before the phobos mission the magnetosphere of mars was relatively unexplored basically all data on the charged particle environment of mars had come from the soviet mars missions mars 2 3 4 and 5 here we report the results of an ion composition experiment on board the soviet phobos 2 spacecraft which allowed us to determine the loss of plasma from the ionosphere of mars because mars has a very small intrinsic magnetic field the ion outflow had been expected to be similar to that of venus surprisingly there were many similarities between the ionospheric outflow from mars and the earth the ion loss from mars results from both ion pick up that results from mass loading of the solar wind in the martian boundary layer and ionospheric o beams of energies up to several kev possibly from upward acceleration processes similar to those observed above the earth s auroral oval a preliminary estimate of the ionospheric outflow from mars indicates that the planet is losing oxygen at a rate of 3 x 1025 ions s 1 this corresponds to an evacuation of its present total atmospheric oxygen content contained in co2 and o2 in less than 100 million years

exploring the recovery lakes region and interior dronning maud land east antarctica with airborne gravity magnetic and radar measurements long range airborne geophysical measurements were carried out in the icegrav campaigns covering hitherto unexplored parts of interior east antarctica and part of the antarctic peninsula the airborne surveys provided a regional coverage of gravity magnetic and ice penetrating radar measurements for major dronning maud land ice stream systems from the grounding lines up to the recovery lakes drainage basin and filled in major data voids in antarctic data compilations such as antgp for gravity data admap for magnetic data and bedmap2 for ice thickness data and the sub ice topography we present the first maps of gravity magnetic and ice thickness data and bedrock topography for the region and show examples of bedrock topography and basal reflectivity patterns the 2013 recovery lakes campaign was carried out with a british antarctic survey twin otter aircraft operating from the halley and belgrano ii stations as well as a remote field camp located at the recovery subglacial lake b site gravity measurements were the primary driver for the survey with two airborne gravimeters lacoste and romberg and chekan am providing measurements at an accuracy level of around 2 mgal r m s supplementing goce gravity field and steady state ocean circulation explorer satellite data and confirming an excellent sub milligal agreement between satellite and airborne data at longer wavelengths

source retrieval for mantle earthquakes by iterative deconvolution of long period p waves abstract long period p waves of eighteen deep and intermediate depth earthquakes in a moment range 1 10 26 3 10 28 dyn cm are analysed by the method of kikuchi and fukao this issue to iteratively derive point dislocations on a plane the fault plane is defined as the plane of point dislocations that gives the best fit between the synthetic and observed waveforms for dip slip type earthquakes the fault plane is found to dip more steeply than the auxiliary plane in general point dislocations tend to cluster into one or a few subevents we characterize the largest subevent by its maximum moment rate m 0 m and its half duration from which the associated stress drop can be roughly estimated for earthquakes shallower than 400 km half duration is approximately constant 5 s so that stress drop increases from 1000 bar as m 0 m increases the value of for the great banda sea earthquake of 1963 h 100 km is exceptionally large 10 s for which 700 bar for earthquakes deeper than 400 km increases from 3 s to 9 s as m 0 m increases so that is approximately constant 700 bar stress drop so obtained approximates to an overall stress drop with a possible exception of the several largest earthquakes the distribution of point dislocations and the moment rate function indicate that rupture most often expands in a fan shaped fault area and that a large part of the rupture front ceases propagation in a coherent manner a sharp contrast in stress drop behavior across a depth of 400 km may be related to the olivine spinel transition within the downgoing slab of lithosphere

subducting plate boundary beneath the northeastern japan arc estimated from sp converted waves abstract in the tohoku district the northeastern part of honshu japan a conspicuous phase is frequently observed between p and s phases in seismograms for events in the lower seismic plane of the double planed deep seismic zone this phase is dominant in the vertical component and has an apparent velocity almost equal to or slightly smaller than the p phase we have never found this phase in seismograms for events in the upper seismic plane from these characteristics the phase is interpreted as sp converted from s to p at the upper boundary of the descending oceanic plate ubp using the travel times of the sp waves we estimated the location of the ubp beneath the tohoku district to be just above the upper seismic plane of the double planed deep seismic zone this position is consistent with those determined from the analyses of the travel times of the scsp and ps waves which are respectively the seismic wave converted from scs to p and that from p to s at the ubp

deep focus earthquakes and recycling of water into the earth s mantle for more than 50 years observations of earthquakes to depths of 100 to 650 kilometers inside earth have been enigmatic at these depths rocks are expected to deform by ductile flow rather than brittle fracturing or frictional sliding on fault surfaces laboratory experiments and detailed calculations of the pressures and temperatures in seismically active subduction zones indicate that this deep focus seismicity could originate from dehydration and high pressure structural instabilities occurring in the hydrated part of the lithosphere that sinks into the upper mantle thus seismologists may be mapping the recirculation of water from the oceans back into the deep interior of our planet

the relationship between wadati benioff zone geometry and p t and b axes of intermediate and deep focus earthquakes for 865 earthquakes with focal depths exceeding 90 km this paper presents a concise global summary of the distributions of p t and b axes of their focal mechanisms with respect to the downdip along strike and normal axes of wadati benioff zones the focal mechanisms are best double couple solutions for the centroid moment tensors reported by the harvard group with the p t and b axes corresponding to the principal axes of the individual moment tensors we find as did previous investigators that events below 300 km have predominantly downdip p axes with 50 of all p axes of deep focus events lying within 28 of downdip for events between 90 and 300 km depth there are regional differences but most regions possess events with nearly downdip p axes or events with nearly downdip t axes in almost all regions having downdip p or t axes bingham statistics show that the greatest concentration of p or t axes lies somewhat below the plane of the wadati benioff zone for events of all focal depths in nearly all regions b axes cluster along the strike direction with 50 occurring within 38 of the strike from these regional results we define a typical focal mechanism as one having the b axis within 30 of the along strike and having either the t or p axis within 30 of the downdip such a focal mechanism might occur in response to a stress tensor having its greatest or least principal stresses channeled approximately along the downdip direction of the subducting slab with the intermediate principal stress lying horizontally within the slab controlled possibly by slab bending however by this criterion only 29 of the focal mechanisms for deep and intermediate events are typical thus the majority of the individual events cannot be explained simply as occurring in response to a downdip and slab normal greatest and least principal stress system with a horizontal along strike intermediate principal stress a few regions such as java and the marianas are clearly anomalous in that they possess significantly fewer typical focal mechanisms than do most other regions

anisotropic growth in the olivine spinel transformation of mg2geo4 under nonhydrostatic stress abstract the olivine spinel phase transformation in mg 2 geo 4 does not occur by a martensitic mechanism the evidence from samples transformed in a griggs type solid medium deformation apparatus are 1 1 lack of microstructural features in the olivine phase which can be specifically associated with a martensitic mechanism 2 2 the orientation relationship between the two phases that is predicted by the martensitic mechanism does not occur nor is there any apparent consistency of relative orientations 3 3 application of a differential stress to the transforming sample resulted in an anisotropic growth rate for the spinel phase indicating that growth was externally controlled rather than crystallographically controlled anisotropic growth of the spinel phase results in elongation of the residual olivine phase grains in the plane normal to the direction of maximum principal compressive stress a velocity ratio of 1 7 0 7 5 4 has been determined for the growth rate of the spinel from measurements on residual olivine grains the interphase grain boundary in samples transformed under stress has cusp shaped fingers of spinel with a blunt end separated by thin spikes of olivine samples transformed isostatically do not exhibit this feature providing further confirmation of anisotropic growth of the spinel the preferred growth of the spinel is consistent with a theory of phase transformation under nonhydrostatic stress the predicted spinel finger shape based on this theory is generally consistent with observed shapes except for the blunt end the discrepancy may be due to surface energy which has not been considered here or to local deviations of the applied macroscopic stress

subducting slabs stagnant in the mantle transition zone the p wave velocity structure beneath the western pacific is found from the international seismological center first arrival data special attention was paid to the deep structure beneath the wadati benioff zone we discretized the whole mantle into blocks with finer blocks in the region of interest to obtain the velocities of all the blocks this way of discretization minimizes a problem with tomographic studies of regional scale difficulty in making corrections for the effects outside the region of interest our solution is iterative with the alternate step of the relocation of earthquakes using the whole mantle model of inoue et al 1990 as a starting model a first order smoothness constraint was imposed to suppress the possible fluctuation of the solution around the initial model the essential result depends little on whether the reference spherical model is smooth or discontinuous near 400 and 670 km depths we examined the resolution by calculating the resolving kernels for selected blocks and by reconstructing the checkerboard test patterns of velocity perturbation and the test structures of subducting lithosphere the resolution is depth dependent but in general good enough to see the slab configuration beneath the southern kurile japan izu bonin arcs and the java arc it is relatively poor beneath the northern kurile and mariana arcs the seismic image of subducting slab beneath the southern kurile to bonin arcs bends to subhorizontal near the leading edge of the wadati benioff zone and extends continentward over a distance of more than 1000 km the subhorizontal portion of the slab connects a high velocity blob to the bottom that reaches a depth of at least 800 km across the 670 km discontinuity under the japan arc although the image of the java slab directly penetrates the 670 km discontinuity it then bends to a shallow dip with a considerable spread reaching a depth of about 1200 km these results suggest that descending slabs of lithosphere in the western pacific tend to be stagnant in the transition zone under a subtle control of the 670 km discontinuity although stagnant slab materials eventually descend into the lower mantle they no longer maintain their original configuration below the 670 km discontinuity

anticrack associated faulting at very high pressure in natural olivine shallow earthquakes are produced by brittle shear fracture of rock and or fictional sliding on pre existing fault surfaces1 at very high pressures however brittle fracture and frictional sliding are impossible because frictional resistance to movement on any potential fault surface is so great that ductile stress relief processes accomodate strain at lower shear stresses than those necessary to activate faulting nevertheless more than 20 of earthquakes with magnitudes greater than five occur at depths greater than 300km where the pressure exceeds l0 gpa ref 2 possible explanations for this paradox have been offered by several recent experimental studies of shearing instabilities associated with phase transformations at high 3 8 one of these mechanisms is associated with anticrack development during the olivine spinel phase transformation in mg2geo4 at pressures of 1 2 gpa refs 4 6 it is particularly attractive because it operates between mineral structures known to occur in the mantle we show here that this mechanism also can operate in natural silicate olivine mg fe 2sio4 during onset of the transformation at the much higher pressures at which deep earthquakes occur these results lend strong support to the hypothesis that the anticrack mechanism is responsible for such earthquakes4

faulting associated with the olivine to spinel transformation in mg2geo4 and its implications for deep focus earthquakes samples of metastable mg2geo4 olivine fail by faulting when loaded at temperatures where the incoherent nucleation and growth of spinel the stable phase just becomes possible at lower temperatures where the transformation is inhibited metastable olivine is ductile and strong while at higher temperatures where the transformation happens very quickly samples are ductile and weak microfaults that have just been created and have slipped very little contain only very fine grained spinel those that have slipped more contain a mixture of fine grained spinel and larger broken olivine fragments a new type of transformation microstructure that develops only under nonhydrostatic stress occurs in all specimens that failed the spinel nucleates on grain boundaries and on subgrain boundaries within olivine and grows as fine grained lens shaped bodies with a strong preferred orientation normal to maximum compression the correlation between the onset of transformation the new microstructure and the mechanical instability suggests a cause and effect relationship the grain size of the spinel within the microfaults suggests that these zones have superplastic properties in this case the instability should not be suppressed by increasing pressure as normal frictional processes are suppressed analogous conditions exist in down going slabs of lithosphere in the mantle and we postulate that the same mechanism that causes faulting in our samples operates in the mantle causing deep focus earthquakes

creep mechanism in mg2geo4 effects of a phase transition we deformed polycrystalline mg2geo4 olivine and spinel in a griggs type solid medium apparatus flow of mg2geo4 olivine can be represented by 6 5 l07 3 5 exp 105 kcal mol rt where kbar is the differential stress and s 1 is the natural strain rate for mg2geo4 spinel the flow law is 2 6 104 2exp 73 kcal mol rt the low stress exponent and activation enthalpy coupled with fine grain size 3 m suggest that mg2geo4 spinel deformed by a superplastic mechanism flow parameters for the olivine phase suggest a dislocation creep mechanism comparison of theoretical superplastic flow laws for mg2geo4 olivine with the spinel phase data suggests that strain rates in mg2geo4 spinel are only about a factor of 3 lower than for mg2geo4 olivine of the same grain size a similar estimate holds for dislocation creep of the two phases if it is controlled by diffusion transformation from mg2geo4 olivine to spinel reduced grain size to approximately 3 m thus we might expect a similar reduction in grain size in the earth s transition zone which could result in superplastic deformation of the transformed phase and cause a weak decoupling zone at the transition boundary in the mantle superplasticity brought about by transformation induced reduction in grain size may also provide a mechanism for deep focus earthquakes and an explanation for the correlation observed between their distribution with depth and the depths of phase transitions within the mantle

the double seismic zone in kuril kamchatka the tale of two overlapping single zones the classic double seismic zone in the central portion of the kuril kamchatka arc seems to turn into a regular wadati benioff zone with downdip compression and extension in the northeastern and southwestern ends of the arc respectively to delineate this lateral changeover we studied 33 large to moderate sized mb 5 1 earthquakes that occurred at depths shallower than 200 km between 1963 and 1991 these events took place in two regions where the double seismic zone terminates and single seismic zones emerge precise earthquake source parameters including focal depths and mechanisms are determined from a careful analysis of a mixture of analog and digital p and sh waveforms recorded at teleseismic distances our results indicate that the double seismic zone is confined between approximately 46 and 53 n northeast of 53 n intermediate depth earthquakes 70 200 km show only downdip compressional focal mechanisms furthermore the wadati benioff zone there seems to be a continuation of the upper compressional layer of the double seismic zone thus the changeover from a double to a single seismic zone in this case is simply disappearing of the lower extensional sheet of the double seismic zone similarly southwest of 46 n intermediate depth earthquakes show only downdip extensional mechanisms the single seismic zone is unlikely to consist of two layers of extensional seismicity in fact nowhere along the entire arc did we observe two distinct seismogenic layers with a single state of strain in other words these extensional events did not occur in the same seismogenic structure as the upper compressional layer of the double seismic zone however we cannot ascertain whether the single seismic zone in extension is simply a southwestern extension of the lower extensional layer of the double seismic zone lateral variation from a single zone under extension in the southwest to a double seismic zone and finally to a single zone under compression cannot be accounted for by conventional models of a double seismic zone on the other hand the addition of a component of stress transmitted from depths greater than 250 km can explain the observed lateral variation in the state of strain at intermediate depths because the slab between depths of 250 450 km is under extension in the southwestern end of the arc whereas the rest of the arc is under compression vanishing of the lower extensional layer of the double seismic zone toward the northeastern end of the arc is interpreted as a result of compression transmitted from below overshadowing extension caused by effects such as unbending a similar reasoning applies to the southwestern termination of the double seismic zone

deep structure of the northeastern japan arc and its relationship to seismic and volcanic activity evidence for possible deep seated magmatic activity beneath northeastern japan can be obtained from seismic observations tomographic inversions of p wave velocity data show low velocity zones distributed in the crust and upper mantle beneath active volcanoes in or around these zones anomalously low frequency micro earthquakes perhaps caused by magmatic activity are found at depths of 25 40 km distinct s wave reflectors corresponding to the upper surfaces of magma bodies are found at depths of 10 18 km large crustal earthquakes also occur around the low velocity zones these observations reveal the distribution of magma reservoirs at depths and elucidate its relation with shallow seismicity beneath volcanic arcs

frequency distribution and quantification of deep earthquakes the frequency distribution of deep earthquakes d 350 km is analyzed on a regional scale using magnitude data from 4660 events reported in the international seismological center bulletin 1964 1982 and 391 centroid moment tensor solutions 1977 1987 we use the moment magnitude mw to define b values comparable with studies of shallow seismicity scaling laws relating mw mb and mo are derived mw 1 6mb 3 2 mb log mo 4 3 2 4 which compare very favorably with published laws uncertainties and biases of the mb and mo data sets are considered the slope of the frequency distribution is evaluated over a restricted magnitude range the results are very robust and independent from the magnitude range used from the mb magnitude data we obtain regional b values in the range 0 4 b 1 22 with a worldwide average value b 0 87 these results are confirmed by the analysis of the seismic moment data the steepest frequency distribution is observed in the tonga islands region b 1 22 although 60 of the worldwide deep seismicity occurs here only few events since 1977 had seismic moment exceeding 1019 n m the high fractal dimension d 2 4 implies a volume partially filled with fault planes in agreement with the observations on the seismic morphology of the region at the other extreme south america is characterized by the largest recorded deep events and an almost complete absence of small ones b 0 4

shear instability in a viscoelastic material as the cause of deep focus earthquakes the mechanism of deep focus earthquakes has been examined by numerical and linear analysis of shear instability in subducting slabs we assume subducting slabs deform such that the spatially averaged strain rate becomes time independent furthermore we assume that rheology depends on temperature in the manner of an activation process using this model we investigate quantitatively the conditions under which shear instability takes place i e the deformation in the material concentrates in a thin layer whose temperature increases explosively until melting occurs first we study instability in a material whose rheological properties are spatially homogeneous shear instability takes place in a homogeneous material when exceeds a certain critical level e c which depends strongly on the stress in the slabs but does not depend on the detail of the creep law the value of the activation energy next we study instability in a material that contains shear zones similar to those often observed by field geologists i e a thin layer made of material with lower vicosity than the surrounding region we find that shear zones trigger shear instability even for the case of provided that g the ratio of characteristic time of thermal diffusion to that of shear heating in shear zones exceeds 1 and that f the ratio of stiffness of elastic deformation in the surrounding region to that of ductile deformation in shear zones becomes lower than 1 applying this condition to subducting slabs with shear zones we find that shear zones trigger shear instability even when the average strain rate is lower than by an order of magnitude from these results and from an estimate of in subducting slabs we conclude that shear instability is a mechanism likely to induce deep focus earthquakes

the temporal distribution of seismic radiation during deep earthquake rupture the time history of energy release during earthquakes illuminates the process of failure which remains enigmatic for events deeper than about 100 kilometers stacks of teleseismic records from regional arrays for 122 intermediate depths of 100 to 350 kilometers and deep depths of 350 to 700 kilometers earthquakes show that the temporal pattern of short period seismic radiation has a systematic variation with depth on average for intermediate depth events more radiation is released toward the beginning of the rupture than near the end whereas for deep events radiation is released symmetrically over the duration of the event with an abrupt beginning and end of rupture these findings suggest a variation in the style of rupture related to decreasing fault heterogeneity with depth

plastic instabilities implications for the origin of intermediate and deep focus earthquakes adiabatic or catastrophic plastic shear has been reported in metals polymers and metallic glasses the phenomenon is associated with rapid stress drops and audible pings or clicks as the material deforms in a plastic manner the driving force for the plastic instability is the stored elastic strain energy of the loading system and in many respects the behavior is reminiscent of the shear stress response arising from stick slip events during unstable frictional sliding although the precise mechanism is different it is shown here that adiabatic plastic shear is capable of explaining the detailed distribution of intermediate and deep focus earthquakes within subduction zones the earthquake events being the result of instabilities in material undergoing plastic flow it is argued that for a particular strain rate there exists a critical temperature tc which is depth dependent for temperatures below tc the material is strain rate softening and for a soft enough loading system may undergo catastrophic plastic shear for temperatures above tc the material is strain rate hardening and is always stable during plastic shear the cutoff depth for deep focus earthquakes then corresponds to the transition from strain rate softening to strain rate hardening material and for commonly accepted geothermal gradients within the slab corresponds to approximately 800 km the precise distribution of earthquakes within the slab is a function of the subtle interplay between the geothermal gradient and the tc gradient in particular a decrease in seismic activity is to be expected below about 300 km in the slab with total stress drops decreasing from a maximum of 700 mpa above 300 km to a maximum of 50 mpa below 300 km the differences in foci distribution between subduction zones such as tonga new hebrides and peru result from minor differences in the geothermal gradients within the slabs the model predicts the development of triple seismic zones high in the slab double seismic zones down to approximately 300 km and single seismic zones down to 800 km such a distribution is to be expected of relatively young cool slabs as the slab heats up the seismic activity retreats up the slab the paper only proposes a deformation mechanism for earthquake generation it does not address the stress field within the slab but only the distribution of strength thus the distribution of focal plane mechanisms is not considered only the locations where earthquakes due to plastic instabilities are possible the absence of earthquakes does not necessarily mean that the slab does not exist it only means that the slab is too hot to undergo plastic instability this means that aseismic subduction is a distinct possibility in many regions of high geothermal gradient within the slab i e circa 3 c km 1