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

impact of shelf basin freshwater transport on deep convection in the western labrador sea abstractfreshwater exiting the arctic ocean through the canadian arctic archipelago caa has been shown to affect meridional overturning circulation and thereby the global climate system however because of constraints of spatial resolution in most global ocean models neither the flow of low salinity water through the caa to the labrador sea nor the eddy activity that may transport freshwater from the shelf to areas of open ocean convection can be directly simulated to address these issues this study uses a high resolution ice ocean model of the pan arctic region with a realistic caa and forced with realistic atmospheric data this model resolves conditions in the arctic ocean upstream of the labrador sea and is coupled to a thermodynamic dynamic sea ice model that responds to the atmospheric forcing the major shelf basin exchange of liquid freshwater occurs south of hamilton bank whereas the largest ice flux occurs in the northwest of the basin freshwater flux anomalies entering the labrador sea

absolute transports of mass and temperature for the north atlantic current subpolar front system the flow of subtropical waters carried into the northern north atlantic ocean by the north atlantic current subpolar front system nac spf is an important component of the meridional overturning circulation these waters become colder and denser as they flow through the subpolar region both by mixing with the colder subpolar waters and by atmospheric cooling the relative roles of these two processes remain to be quantified and the mechanisms driving lateral mixing need to be better understood to address those questions a new methodology is developed to estimate the mean absolute transports of mass and heat for the top 1000 dbar in the region of the nac spf for the time period 1993 2000 the transports are obtained by combining historical hydrography with isopycnal rafos float data from the area the mean absolute transport potential field shows an nac spf pipe defined by two bounding transport potential contours this pipe transports 10 0 6 3 5 sv sv 106 m3 s 21 top 1000 dbar from the subtropics into the eastern subpolar north atlantic in contrast to earlier studies the northward flowing nac follows a distinct meandering path with no evidence of permanent branches peeling off the current before reaching the northwest corner as the current enters the northwest corner it loses its tight structure and maybe splits into two or more branches which together constitute the eastward flow along the spf the eastward flow between the northwest corner and the mid atlantic ridge is not as tightly defined because of the meandering and or eddy shedding of the branches constituing the spf as the flow approaches the mid atlantic ridge it converges to cross above the charlie gibbs and faraday fracture zones the mean absolute temperature transport top 1000 dbar by the 10 sv pipe was estimated across 10 transects crossing the nac spf because the mean mass flux is constant in the pipe variations in the mean temperature transports result from lateral exchange and mixing across the pipe s side walls and from air sea fluxes across the surface of the pipe the nac spf current loses 0 18 6 0 05 pw on its transit through the region most of the loss occuring upstream of the northwest corner the heat loss is 10 times the corresponding heat lost to the atmosphere we conclude that cross frontal exchange induced by the steep meanders of the northward flowing nac is the main mechanism by which heat is lost along the current in the region between the tail of the grand banks and the mid atlantic ridge

relations among subduction parameters clues to the dynamics of the subduction process are found in the many measurable parameters of modern subduction zones based on a critical appraisal of the geophysical and geological literature 26 parameters are estimated for each of 39 modern subduction zones to isolate causal relationships among these parameters multivariate analysis is applied to this data set this analysis yields empirical quantitative relations that predict strain regime and strike slip faulting in the overriding plate maximum earthquake magnitude benioff zone length slab dip arc trench gap and maximum trench depth excellent correlation is found between length of the benioff zone and the product of convergence rate and age of the downgoing slab this relationship is consistent with the conductive heating model of molnar et al 1979 if the model is modified in one respect the rate of heating of the slab is not constant it is substantially slower during passage of the slab beneath the accretionary prism and overriding plate the structural style in the overriding plate is determined by its stress state though the stress state of overriding plates cannot be quantified one can classify each individual subduction zone into one of seven semiquantitative strain classes that form a continuum from strongly extensional class 1 back arc spreading to strongly compressional class 7 active folding and thrusting this analysis indicates that strain class is probably determined by a linear combination of convergence rate slab age and shallow slab dip interplate coupling controlled by convergence rate and slab age is an important control on strain regime and the primary control on earthquake magnitude arc parallel strike slip faulting is a common feature of convergent margins forming a forearc sliver between the strike slip fault and trench optimum conditions for the development of forearc slivers are oblique convergence a compressional environment and a continental overriding plate the primary factor controlling presence of strike slip faulting is coupling strongly oblique convergence is not required the rate of strike slip faulting is affected by both convergence obliquity and convergence rate maximum trench depth is a response to flexure of the underthrusting plate the amount of flexural deflection at the trench depends on the vertical component of slab pull force which is very sensitive to slab age and shallow slab dip shallow slab dip conforms to the cross sectional shape of the overriding plate which is controlled by width of the accretionary prism and duration of subduction deep slab dip is affected by the mantle trajectory established at shallow depth but may be modified by mantle flow much of the structural complexity of convergent margins is probably attributable to terrane juxtaposition associated with temporal changes in both forearc strike slip faulting and strain regime empirical equations relating subduction parameters can provide both a focus for future theoretical studies and a conceptual and kinematic link between plate tectonics and the geology of subduction zones

s p conversion from the transition zone beneath tonga and the nature of the 670 km discontinuity summary using more than 300 recordings of deep tonga earthquakes from the warramunga seismic array central australia we have searched for short period s p conversions from the 670 km discontinuity arriving in the p wave codas the array seismograms were stacked and plotted in sections to optimize observation of s670p yielding the following results for the northern part of tonga 17 3 18 4 s s670p phases are observed with apparent conversion depths between about 660 and 680 km for the central part 20 2 22 3 s s670p is not clearly observed for the southern part 23 1 25 7 s s670p is observed for conversion depths between about 660 and 700 km intermittent occurrence of s p conversion may result from different characteristic source radiation patterns among the three sections overall uncertainties in conversion depths resulting from unmodelled lateral velocity variations and hypocentral location errors are perhaps a few tens of kilometres the s p rays must travel either within or very close to the subducted slab because of the source receiver slab dip geometry so we conclude that large 50 km deformations of the 670 km discontinuity do not occur beneath tonga our observations appear to be inconsistent with a chemical discontinuity at 670 km but are easily explained if the 670 km discontinuity is a strongly pressure dependent and perhaps moderately endothermic phase change in both the slab and surrounding mantle

intermediate and deep seismicity and lateral structure of subducted lithosphere in the circum pacific region in this paper we present a region by region review of the wadati benioff zone structure of most of the world s seismically active subduction zones focusing primarily on the intermediate and deep seismicity lateral changes in wadati benioff zone structure are common in every major subduction zone in this study we use these changes to define possible boundaries between portions or segments of lithosphere with differing subduction geometries although earthquake data seldom have the resolution to show conclusively whether these boundaries separate independent blocks of lithosphere the available data indicate that the active process at most of these segment boundaries is ductile deformation of the subducting plate rather than tearing we found the strongest evidence for the existence of tears where tears are geometrically necessary such as where a transform boundary terminates a trench as at the new hebrides tonga and south sandwich trenches weak evidence suggesting other tears does exist in some regions such as taiwan the japan izu bonin corner and the philippines the causes of these changes in structure are in most cases unclear only about 34 of the possible segment boundaries coincide with subducting bathymetric features some boundaries occur where there is apparent lateral strain caused by anomalous trench geometry we have designed a simple modeling procedure which incorporates published plate motion and the observed geometry of trenches and wadati benioff zones to estimate the lateral strain in subduction zones throughout the circum pacific region although no observed subduction zone has a perfectly strain free geometry there is a broad range of geometries for which the lateral strain is small indeed the observed geometry of most subduction zones involve relatively little lateral strain comparison with centroid moment tensor focal mechanisms indicates that in zones where the modeling predicted little lateral strain the mechanisms of intermediate and deep earthquakes show no effects of lateral stress and downdip stresses are clearly dominant in regions such as the mariana arc where the model predicts very large lateral extension lateral tension is very evident in the focal mechanisms in regions such as the hokkaido corner where the modeling predicts large compressional strains the plate appears to buckle and bending stresses parallel to the trench are evident in general this study finds that subducted lithosphere is remarkably cohesive and rigid and only rarely deforms by breaking or stretching

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

the state of stress in subducting slabs as revealed by earthquakes analysed by moment tensor inversion abstract we have used a data set of earthquakes whose mechanisms have been studied here or elsewhere by moment tensor inversion to confirm some earlier ideas concerning the state of stress in subducting slabs the moment tensor method offers the advantages of speed and objectivity over earlier methods usually based exclusively on first motion data deep earthquakes below 300 400 km tend to have compression axes aligned with the dip of the subducting slab the behavior of intermediate earthquakes from 70 to 300 km is less clear there is some evidence for down dip tension but the alignment is not as obvious as it is in the case of the compression axes of deep events tonga shows a tendency toward down dip compression at intermediate depths but other deep extending slabs as well as shallow extending slabs slabs with maximum depths shallower than 300 km seem to be closer to down dip tension

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

deep seismicity and modes of deformation in tonga subduction zone the geometry of deep 400 km earthquakes in the subducting western edge of the pacific plate when related to the local distribution of seismic activity reveals a complex pattern of planes of shear failure enabling the modes of deformation by which the descending plate shortens and thickens to be identified the deepest part of the slab is displaced laterally suggesting that it fails to penetrate the discontinuity at 670 km

the north atlantic current and surrounding waters at the crossroads the north atlantic current is a well defined western boundary current that flows north along the east side of the grand banks from 40 to 51 n where it turns sharply to the east and begins its journey across the ocean the current is unique in transporting warm tropical waters to much higher latitudes than any other western boundary current and thus plays a crucial role in ameliorating the climate of the european subcontinent the north atlantic current originates in the gulf stream when the latter curves north around the southeast newfoundland rise a major submarine ridge that stretches se from the grand banks a well defined front delineates the path of the current as long as it flows north as a western boundary current after the current turns east in the north it broadens into a widening band of eastward drift without a sharp or permanent front in the sense of the eastward flowing gulf stream after it separates from cape hatteras the north atlantic current transports more than 40 sv 1 sv 106 m s 1 in the south and about 20 sv by the time it flows east across the mid atlantic ridge the currents along the northward flowing front are quite swift with typical maximum average speeds in the upper 300 m near 1 m s 1 2 knots the current meanders almost as wildly as a snaking river but unlike steep meanders in the gulf stream these meanders appear to be stable and with one exception have not been observed to break off to form pools of warm and or cold waters as frequently occurs in the gulf stream the meanders appear to be induced by major topographic features along the path of the current namely the southeast newfoundland rise the newfoundland seamounts and flemish cap strong recirculations develop on the concave side of the meanders one of these the mann eddy at the first meander crest of the north atlantic current should be regarded as a permanent feature of the north atlantic circulation other meanders also contain recirculations that can persist for months under certain conditions these can merge together to form an extended sw flow recirculation just east of the north atlantic current

the effect of the descending lithosphere beneath the tonga island arc on p wave travel time residuals at the warramunga seismic array abstract p wave travel time residuals at the warramunga seismic array wra in the northern territory australia have been studied from 49 earthquakes with epicenters south of 19 s in the fiji tonga region focal depths are between 42 and 679 km as determined from pp p using the jeffreys bullen and the herrin travel time tables the epicentral parameters have been redetermined by considering only normal seismic stations in the location procedure these are those stations where p wave travel times are probably not affected by lateral heterogeneities caused by the lithosphere descending beneath the tonga trench epicenters of deep earthquakes below 300 km have been relocated by using stations at 25 only epicenters from shallower depth earthquakes have been recalculated without using stations between 35

an alternative hypothesis for the origin of the mediterranean salt lens observed off the bahamas in the fall of 1976 abstract a hypothesis is presented that the original salt lens or meddy observed off the bahamas in the fall of 1976 may have been formed not near the mediterranean outflow but instead in the vicinity of the northwest corner 51 n 43 w of the north atlantic current an eddy was observed near the northwest corner by an isopycnal rafos float deployed during the 1993 95 north atlantic current experiment and had nearly identical temperature salinity properties as those of the bahamas lens hydrographic evidence of thick homogeneous layers with similar properties near the northwest corner suggest a possible formation mechanism by which surface eddies containing warm and saline waters are cooled and subducted a plausible scenario is made whereby a northwest corner eddy might be advected southward in the newfoundland basin by the flow around the high pressure ridge east of the north atlantic current and then enter the recirculation gyre immediately south of the gulf stream such an eddy could be advect

resolution experiments for nw pacific subduction zone tomography we investigate the resolving power of isc neic p travel time data in tomographic inversions for the geometry of the subduction zones in the nw pacific from thermal models for the kurile japan izu bonin mariana and ryukyu slabs we generate three dimensional synthetic velocity anomalies for subducting slabs which are then projected onto a cell model for the uppermost 1400 km of the mantle using these synthetic models we compute synthetic delay times for ray paths corresponding to the source and receiver locations used for the actual data add gaussian noise invert the synthetic data and compare the resulting velocity structure to the initial synthetic models this comparison is illustrated for sections through the kuriles and the mariana arcs a variety of resolution artifacts are observed which in many cases resemble features visible in the tomographic results obtained from inverting the actual isc neic data

seismic imaging of upper mantle structure with new evidence for a 520 km discontinuity stacked images of long period seismic data reveal many phases some not previously observed which are caused by reflections and conversions at discontinuities in the earth s upper mantle these images support the existence of discontinuities at 410 520 and 660km depth but no evidence is found for a significant interface at 220 km

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

observations in the northwest corner of the north atlantic current abstract the north atlantic current a northern branch of the gulf stream flows in a northerly direction to the east of the grand banks of newfoundland and turns sharply to the right north of flemish cap in the northwest corner the currents and water masses in and near this turn are described using data from batfish satellite tracked drifters ctds and current meters drifter trajectories and water mass properties indicate that the northward current toward the northwest corner is sometimes well defined narrow converging and accelerating at other times it is poorly defined weaker less coherent and diverging when it is well defined the current turns sharply eastward between 51 and 52 n but rather than continuing east it forms 100 km anticyclonic eddies which persist for at least 1 2 months records from 16 current meters placed on four moorings across the current show the flow direction to be variable in time but nearly constant throughout the water column speed is also constant with depth

earthquake generating stresses in a descending slab abstract a numerical study of stress within descending slabs beneath island arcs provides a general interpretation for the occurrence of deep and intermediate focus earthquakes the computation of stress due to several factors which have been proposed for the interpretation of the double planed seismicity of intermediate depth events shows that thermal expansion and olivine spinel phase change play an important role in earthquake generation the stress field due to the combination of these two causes gives the following results in the depth range 50 300 km the compressional and tensional stress predominates in the upper and central part of the slab respectively maximum shear stress amounting to 1 6 gpa in compressional and 0 4 gpa in tensional fields occurs at a depth of 200 km the stress is generally very small at around 300 350 km below 400 km compressional stress exists at the central part of the slab and tensional stress appears near the upper and lower surfaces of the slab maximum shear stress is 2 8 gpa in compression and 1 6 gpa in tension at a depth of 500 km the principal axes in these high stress fields align almost parallel to the dip direction of the slab these stress features explain well such observed seismological facts beneath island arcs as the seismicity depth relationship planer configuration of seismic zone existence or non existence of double planed seismicity and focal mechanisms of intermediate and deep focus earthquakes

spherical shell tectonics effects of sphericity and inextensibility on the geometry of the descending lithosphere the shape of a deep seismic zone is thought to represent that of the descending slab of lithosphere the lithosphere before subduction is a spherical shell and the shape of the descending slab is the result of the deformation of the spherical lithosphere at the subduction zone upon bending a spherical shell often deforms in a very different way from a simple plate we examine whether the actual shape of the descending slab can be explained by a simple bending of an inextensible spherical shell which shows little surface deformation under moderate stress this examination is made region by region for most of the subduction zones in the world by means of an analogue method the lithosphere is simulated by an inextensible spherical shell made of polyvinyl chloride resin the wadati benioff zone is shaped by plaster by referring to the reliable hypocentral data selected from the international seismological centre isc bulletins the spherical shell is forced to fit the miniature of the wadati benioff zone fitting is first attempted only by bending if a good fit is not attainable and if a discontinuity or gap in seismic activity is observed in the relevant region the spherical shell is torn along this discontinuity or gap and the goodness of fit is reexamined the results of the analysis are summarized as follows 1 the shape of the wadati benioff zone can be simulated largely by a simple bending of a spherical shell without surface deformation 2 in almost all of the regions of poor fit with bending a good fit can be achieved by tearing the spherical shell along the trace of low seismicity the sphericity of the lithosphere and the inextensibility upon deformation are the two essential factors in controlling the slab shape this means that the lateral constraint is most important for understanding the geometry of the downgoing slab of lithosphere and the stress state within it further several problems related to the deformational characteristics of the spherical lithosphere are also reviewed and discussed in connection with subduction tectonics

subduction zone earthquakes and stress in slabs the pattern of seismicity as a function of depth in the world and the orientation of stress axes of deep and intermediate earthquakes are explained using viscous fluid models of subducting slabs with a barrier in the mantle at 670 km 670 km is the depth of a seismic discontinuity and also the depth below which earthquakes do not occur the barrier in the models can be a viscosity increase of an order of magnitude or more or a chemical discontinuity where vertical velocity is zero longn versus depth wheren is the number of earthquakes shows 1 a linear decrease to about 250 300 km depth 2 a minimum near that depth and 3 an increase thereafter stress magnitude in a subducting slab versus depth for a wide variety of models shows the same pattern since there is some experimental evidence thatn is proportional toe where is a constant and is the stress magnitude the agreement is encouraging in addition the models predict down dip compression in the slab at depths below 400 km this has been observed in earlier studies of earthquake stress axes and we have confirmed it via a survey of events occurring since 1977 which have been analysed by moment tensor inversion at intermediate depths the models predict an approximate but not precise state of down dip tension when the slab is dipping the observations do not show an unambiguous state of down dip tension at intermediate depths but in the majority of regions the state of stress is decidedly closer to down dip tension than it is to down dip compression chemical discontinuities above 670 km or phase transitions with an elevation of the boundary in the slab predict when incorporated into the models stress peaks which are not mirrored in the profile of seismicity versus depth models with an asthenosphere and mesosphere of appropriate viscosity can not only explain the state of stress observed in double benioff zones but also yield stress magnitude profiles consistent with observed seismicity models where a nonlinear rheology is used are qualitatively consistent with the linear models

current broadening as a mechanism for anticyclogenesis at the northwest corner of the north atlantic current 1 the warm north atlantic current flowing north just east of the grand banks makes an abrupt turn to the east near 51 n as the waters flow through this region known as the northwest corner the current broadens and weakens effectively splitting into two zonal yet variable flows forming what is generally known as the subpolar front the northwest corner also has the curious property that it expands and expels anticyclonic vortices on a quasi regular basis approximately every 2 1 2 months using an energy budget approach we show that this anticyclogenesis can be accounted for in terms of a kinetic energy convergence due to the broadening of the flow as it passes through the region the principal export route for the newly formed eddies appears to be density driven subduction i e sinking driven by severe heat loss due to the exposed position of these warm waters in the southern labrador sea

mantle dynamics and basalt petrogenesis abstract differentiation at mid ocean ridges generates a layered lithosphere consisting of a basaltic crust immediately underlain by harzburgite and further underlain by pyrolite which has experienced depletion only of highly incompatible elements the body forces driving subduction are concentrated mainly in the upper half of the lithosphere which is relatively cool and brittle during subduction the lower layer of relatively ductile slightly depleted pyrolite is stripped off and resorbed into the upper mantle thereby providing a future source region for morb magmas the slab which sinks to 600 km is comprised mainly of differentiated former basalt and harzburgite which undergo a different series of phase transformations to those experienced by mantle pyrolite in consequence the former basaltic crust remains denser than surrounding mantle whereas former harzburgite becomes relatively buoyant below the 650 km seismic discontinuity the resulting non uniformity in stress distribution causes the slab to buckle at this depth and accrete to form a large relatively cool ovoid megalith of mixed former harzburgite and basaltic crust heating of the megalith occurs over 1 2 b y leading to partial melting of the former basaltic crust the resultant liquids contaminate regions of former harzburgite rendering them fertile in the sense of future capacity to produce basaltic magmas after thermal equilibration the newly fertile former harzburgite becomes buoyant leading to the separation of diapirs which rise into the upper mantle such diapirs rising beneath sub oceanic lithosphere experience small degrees of partial melting to produce ocean island basalts mainly of the alkaline suite diapirs of fertile former harzburgite rising beneath continents become incorporated into the sub continental lithosphere this is a cumulative process and is ultimately responsible for the development of the chemical physical and isotopic characteristics of the sub continental lithosphere including its capacity to produce alkaline magmas when subjected to small degrees of partial melting

waveform effects of a high velocity subducted slab creager and jordan 1986 propose that penetration of subducting slabs under the kurile islands and other subduction zones to depths of at least 1000 km is necessary to explain travel time anomalies of deep earthquakes such penetration would also be expected to affect the amplitudes and waveforms of the body waves from earthquakes synthetic seismograms appropriate for a record section in a plane perpendicular to the strike of the slab are presented using a coupled finite difference and kirchhoff method our shear wave version of the compressional wave velocity structure of creager and jordan 1986 produces an amplitude decrease up to a factor of four and waveform broadening up to 20 seconds for sh arrivals with a take off angle pointing straight down the slab slabs that extend only 300 km below the earthquake but are half as thick and twice as anomalously fast as creager and jordan s 1986 velocity model will roughly preserve the travel time variation pattern and show less waveform broadening but produce first arrivals that are emergent slabs that become thicker with depth show less waveform broadening reconciliation of the amplitude waveform distortion and timing of body waves from deep events is necessary to understand the geometry of slabs in the upper mantle and near and below the 650 km discontinuity

the water masses of the central north atlantic in 1983 84 abstract hydrographic surveys carried out in 1983 84 along both sides of the mid atlantic ridge between 24 and 53 n provide a detailed description of the well known north atlantic water masses with particular emphasis on their meridional distribution and zonal dissymetry in the upper layers the dense horizontal sampling resolves the several narrow gulf stream extensions into the ocean interior giving the image in the central water density range of a mosaic of mode waters separated by fronts at intermediate depths a vertical shear in the distributions of the mediterranean water and labrador sea water stands out both water masses having their lower part displaced southwards relative to their upper parts bottom waters containing nearly 20 percent of pure antarctic bottom water are observed at 50 n in the eastern basin in contrast with the western basin where proportions greater than 10 percent were found only south of 36 n along our section this water mass analysis also gives indications that stron

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

modeling of residual spheres for subduction zone earthquakes 1 apparent slab penetration signatures in the nw pacific caused by deep diffuse mantle anomalies we have computed focal residual spheres for 145 subduction zone earthquakes along the northwest edge of the pacific using regional and global mantle velocity models from tomographic inversions the mantle models explain much of the observed residual sphere data and to a certain extent suggest the location of mantle velocity heterogeneities which are responsible for various residual sphere patterns for most deep events considered the fast slablike residual sphere anomalies are caused by diffuse heterogeneities mainly of deep lower mantle and receiver mantle origin rather than by an extension of the slab the region immediately below the deepest earthquakes depths of 650 1500 km has an effect usually smaller than or comparable to the effect of other regions of the mantle without a proper account of the teleseismic effect attributing the long wavelength anomalies of the residual sphere to near source slab effects alone or even primarily is not valid the fast bands in many observed residual spheres agree with seismicity trends once the deep mantle and receiver mantle effects are removed these may give the approximate orientation but not the depth extent of near source fast velocities for most deep earthquakes under japan the predominant fast band is subhorizontal rather than near vertical this type feature would be overlooked in conventional residual sphere studies using only steeply diving rays and cosine weighting of the data

horizontal shear flow in the mantle beneath the tonga arc in the tonga subduction zone a systematic shear deformation in which deep material moves relatively to the south is superimposed on the mode of down dip shortening which is characteristic of the region we interpret this observation together with a tectonic reconstruction in terms of the interaction between the subducting lithospheric slab and a horizontal shear flow in the mantle

search for deep slabs in the northwest pacific mantle a residual sphere is formed by projecting seismic ray travel time anomalies relative to a reference earth model onto an imaginary sphere around an earthquake any dominant slab like fast band can be determined with spherical harmonic expansion the technique is useful in detecting trends associated with high velocity slabs beneath deep earthquakes after deep mantle and near receiver effects are removed two types of corrections are used the first uses a tomographic global mantle model the second uses teleseismic station averages of residuals from many events over a large area centered on the events of interest under the mariana izu bonin and japan trenches the dominant fast bands are generally consistent with seismicity trends the results are unstable and differ from the seismicity trend for kurile events the predominant fast band for most deep earthquakes under japan is subhorizontal rather than near vertical we find little support for the deep slab penetration hypothesis

on the dynamics of magma mixing by reintrusion implications for pluton assembly processes many plutons have formed by repeated intrusion with complex internal contacts characterized by sharp and diause kinematic and compositional domains we performed numerical experiments of mixing following magma chamber recharge that explicitly considers the fluid dynamics of multiphase mixtures in the limit of zero diausivity of intensive scalar quantities and low stokes number mingling by fluid instability is equivalent to deformation and persistent fluid structures are kinematic attractors three distinct regimes are exemplified and can be described by their multiphase reynolds or grashof number for a reynolds number greater than about 100 an internal intrusive contact will collapse by internal wave breaking and chaotic magma mingling and mixing yield a nearly chamber wide stratification this flow may scour mushy regions at the walls and widely distribute previously crystallized material for reynolds numbers from 10 to 100 the internal wave does not break and the stratification occupies less of the chamber leaving islands of unmixed material for reynolds numbers around 1 internal slumping and folding can occur and is the most likely to be preserved by a mineral fabric internal sub vertical contacts require a high absolute viscosity in the resident magma representing a time break between episodes of reintrusion the shape of perched mafic domains captures the progress of sinking from higher levels into regions of increasing strength and crystallinity 7 2000 elsevier science ltd all rights reserved