volcanism in the afar rift sustained by decompression melting with minimal plume influence continental breakup is caused by some combination of heating and stretching1 2 the afar rift system in africa is an example of active continental rifting where a mantle plume probably weakened the lithosphere through thermal erosion and magma infiltration however the location and degree of plume influence today are debated2 3 here we use seismic s to p receiver functions to image the mantle structure beneath afar we identify the transition between the lithosphere and underlying asthenosphere at about 75 km depth beneath the flanks of the continental rift however this boundary is absent beneath the rift itself and we instead observe a strong increase in seismic velocities with depth at about 75 km we use geodynamic modelling to show that the velocity increase at this depth is best explained by decompression melting of the mantle in the absence of a strong thermal plume so although the absence of mantle lithosphere beneath the rift implies a plume may have once been active we conclude that the influence of a thermal plume directly beneath afar today is minimal

radial thermo chemical structure beneath western and northern pacific from seismic waveform inversion abstract the earth s deep mantle seismic structure is dominated by two large low shear wave velocity provinces llsvps located beneath africa and the pacific while the existence of these structures has been attested by many studies and data sets their detailed nature purely thermal or thermo chemical is still a matter of debate discriminating between these hypotheses requires constraints independent from seismic velocity structure seismic shear wave attenuation measured by the quality factor q s strongly depends on temperature but not or weakly on composition it may bring key information on temperature resolving in turn the trade off between temperature and composition here we invert seismic waveform data for radial models of shear wave velocity anomalies dln v s and q s at two different locations in the pacific and from a depth of 2000 km down to the core mantle boundary cmb we show that western pacific wp models sampling the western tip of the pacific llsvp and the caroline plume cannot be explained by thermal anomalies alone but require excess in iron of 4 0 from the cmb up to 2600 km by contrast northern pacific models np if unaffected by seismic anisotropy may have a purely thermal origin based on these observations we build radial thermo chemical models at wp and np

sedimentary and structural evolution of the northwest arabian red sea margin abstract an investigation of the northwestern margin of the arabian peninsula midyan region provides new insights into the two stage development of the red sea graben system the early red sea suez stage and the later red sea aqaba stage the early stage is characterized by the passive subsidence of a continental graben since the early oligocene after a period of enhanced subsidence of the basin in the early miocene this stage ended at the langhian serra valian boundary 14 ma the later stage was initiated by a release of stresses by strike slip movement along the newly formed aqaba levant transform creation of the new plate boundary was accompanied by stagnation of extension in the gulf of suez branch of the early rift stage strike slip motion on the aqaba levant structure is younger than 14 ma and is recorded by the post langhian compressional structures in the area of the midyan peninsula which developed by transpression along the curved path of a left lateral strike slip fault

seismic imaging of melt in a displaced hawaiian plume the precise location of the mantle plume upwelling beneath hawaii is debated seismic data reveal a thick layer of melt in the mantle beneath western hawaii implying that the upwelling plume may be deflected around an ancient resistive root beneath the island

deformation of subducted oceanic lithosphere summary the deformation of subducted oceanic lithosphere is revealed by detailed studies of the distribution of seismicity in benioff zones and by the increasingly well resolved images of the upper mantle obtained by seismic p wave tomography by using numerical experiments based on a simplified dynamical model of a subducted lithospheric slab we test the idea that the observed deformation of the subducted slab results from a balance between internal buoyancy forces and the viscous stresses associated with deformation the simplified model assumes that the lithosphere has uniform physical properties it is denser and much more viscous than the upper mantle and its penetration below 670 km depth is resisted by a step like density increase at that level the primary determinant of the style of deformation is the buoyancy number f a ratio of buoyancy stress generated by the mass anomaly in the slab to the stress associated with viscous deformation at a strain rate u l defined in terms of subduction rate u and slab thickness l for a small buoyancy number deformation of the slab is dominated by viscous flexure for large f down dip extension of the slab appears to dominate and a buckling instability of the slab is observed if there is resistance to penetration of the 670 km level the transition value of f at which this change in behaviour is observed is about 0 05 for constant viscosity and about 0 2 for stressdependent viscosity with a stress versus strain rate exponent of n 3 with n 3 a boudinage type instability of the slab is observed for large f penetration of the slab below the 670 km level depends on the density contrast between slab and lower mantle resisted only by the density difference the slab may initially penetrate several hundred kilometres below the 670 km level for plausible density parameters before a buckling instability causes this part of the slab to rotate and ascend comparison of the deforming slab geometry and stress field with seismicity distribution and tomographic images from the tonga subduction zone suggests that the effective buoyancy number f for this slab is approximately equal to the transition value that we determined experimentally we use this constraint to estimate average rheological parameters for the tonga slab which when compared with published creep deformation laws for olivine are consistent with the average temperature of the slab being about 0 4 to 0 45 times the melting temperature based on constitutive laws for olivine at shallow depths 100 km these temperatures correspond to around 590 c for wet olivine or around 700 c for dry olivine

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

global positioning system measurements of present day crustal movements in the arabia africa eurasia plate collision zone we present and interpret global positioning system gps measurements of crustal motions for the period 1988 1994 at 54 sites extending east west from the caucasus mountains of southern russia georgia and armenia to the aegean coast of turkey and north south from the southern edge of the eurasian plate pontus block to the northern edge of the arabian platform viewed from a eurasia fixed reference frame sites on the northern arabian platform move n38 13 w at 20 3 mm yr roughly consistent with the velocity implied by nuvel 1a circuit closure n23 7 w at 24 2 mm yr the motion of arabia appears to be transferred directly to the region of turkey north of the suture however eastern turkey is characterized by distributed deformation while central western turkey is characterized by coherent plate motion involving westward displacement and counterclockwise rotation of the anatolian plate internal deformation within the central part of the anatolian plate is less than 2 mm yr the anatolian plate is decoupled from eurasia along the right lateral strike slip north anatolian fault naf this different response in eastern and western turkey to the collision of arabia may result from the different boundary conditions the hellenic arc forming a free boundary to the west and the asian continent and oceanic lithosphere of the black and caspian seas forming a resistant boundary to the north and east we derive a best fitting euler vector for anatolia eurasia motion of 29 2 0 8 n 32 9 0 4 1 3 0 1 m y the mapped surface trace of the naf corresponds well to a small circle about this pole the new euler vector implies an upper bound for naf slip rate of 30 2 mm yr i e assuming all relative motion is accommodated along the naf using the nuvel 1a euler vector for arabia eurasia and the gps euler vector for anatolia eurasia we determine an arabia anatolia euler vector of 31 2 n 45 2 e 0 9 0 1 m y and an upper bound on the east anatolian fault slip rate of 15 3 mm yr the aegean trough region of western turkey deviates significantly from coherent plate rotation in addition to rotating with anatolia this region shows roughly n s extension at a rate of 14 5 mm yr taken together with satellite laser ranging results along the hellenic arc the contemporary pattern of deformation indicates increasing motions toward the arc suggesting that the westward displacement and counterclockwise rotation of anatolia is driven both by pushing from the arabian plate and by pulling or basal drag associated with the foundering african plate along the hellenic subduction zone

isostatic response of the alula fartak and owen transforms in the eastern gulf of aden and the adjoining arabian sea summary the alula fartak and owen transforms are the active parts of major fracture zones with a distinct topographic expression in the eastern gulf of aden and the arabian sea respectively while the alula fartak transform offsets the sheba ridge by about 180 km and is associated with a broad steep sided valley with a relief of nearly 3 5 km the owen transform offsets the carlsberg ridge by nearly 300 km and is associated with a broad step like valley surrounded by deeper water depths the gravity and topography data along several profiles selected across these two transforms have been analysed using cross spectral analysis in order to investigate their isostatic compensation the observed admittance estimates have been compared with three theoretical isostatic compensation models two local compensation models airy i and ii and one regional compensation plate model comparing the longer wavelength admittance estimates suggests that the regional compensation model gives the best fit for both the alula fartak transform and the owen transform with effective elastic thickness te of 5 km and slightly less than 10 km respectively for the alula fartak transform the airy ii model might also be acceptable though with large scatter in the observed values it suggests a mean value of 9 km for the mantle layer with a 6 km thick crust for the owen transform on the other hand the two local compensation models failed the difference in te estimate between the two transforms could be ascribed to differences in thermal structure arising from their varied tectonic history a comparison with the isostatic response estimates of transform fracture zones along the slow spreading mid atlantic ridge suggests that the regional compensation model is generally applicable for transform fracture zone topography along such mid ocean ridges

dynamics of retreating slabs 2 insights from three dimensional laboratory experiments 1 a laboratory analogue of a three layer linear viscous slab upper mantle lower mantle system is established in a silicone putty honey and crystallized honey tank experiment the same setup as in the numerical investigation part 1 is used we focus on the interaction of the slab with the induced passive mantle flow by widely varying the mantle volume flux boundary conditions in our numerical experiments the lateral volume flux was set to zero in interpreting the results relative to the real earth the base of the box is taken as the bottom of the mantle convection system while the lateral boundaries may be associated with the presence of other nearby slabs dynamic force equilibrium assessed on the basis of an analytical review of forces is described for four different phases 1 the subduction initiation instability 2 the accelerating dynamic free fall phase of the slab 3 the dynamic interaction with the 660 km discontinuity and 4 a final phase of steady state trench retreat phase 3 is an important feature not observed in the numerical experiments this highly dynamic phase of interrupted trench retreat can therefore be attributed to boundary conditions on mantle volume flux on the basis of integration constants of force equilibrium in phases 2 and 4 we identify two different classes of volume flux one in which the lateral boundary can be considered open and the other class where it is closed closed boundary condition cases are obtained if any of the lateral box boundaries are 600 km away from the slab assuming a one to one relation between trench retreat and back arc spreading enigmatic observations of episodic opening of back arc basins can be explained by our experimental observations

constraints on the presence of post perovskite in earth s lowermost mantle from tomographic geodynamic model comparisons abstract lower mantle tomography models consistently feature an increase in the ratio of shear wave velocity v s to compressional wave velocity v p variations and a negative correlation between shear wave and bulk sound velocity v c variations these seismic characteristics also observed in the recent sp12rts model have been interpreted to be indicative of large scale chemical variations other explanations such as the lower mantle post perovskite ppv phase which would not require chemical heterogeneity have been explored less constraining the origin of these seismic features is important as geodynamic simulations predict a fundamentally different style of mantle convection under both scenarios here we investigate to what extent the presence of ppv explains the observed high v s v p ratios and negative v s v c correlation globally we compare the statistical properties of sp12rts with the statistics of synthetic tomography models derived from both thermal and thermochemical models of 3 d global mantle convection we convert the temperature fields of these models into seismic velocity structures using mineral physics lookup tables with and without ppv we account for the limited tomographic resolution of sp12rts using its resolution operator for both v s and v p structures this allows for direct comparisons of the resulting velocity ratios and correlations although the tomographic filtering significantly affects the synthetic tomography images we demonstrate that the effect of ppv remains evident in the ratios and correlations of seismic velocities we find that lateral variations in the presence of ppv have a dominant influence on the v s v p ratio and v s v c correlation which are thus unsuitable measures to constrain the presence of large scale chemical variations in the lowermost mantle to explain the decrease in the v s v p ratio of sp12rts close to the cmb our results favour a ppv bearing cmb region which has implications for the stability field of ppv in the earth s mantle

insights on upper mantle melting rheology and anelastic behavior from seismic shear wave tomography in seismic tomography we map the wave speed structure inside the earth but we ultimately seek to interpret those images in terms of physical parameters this is challenging because many parameters can trade off with each other to produce a given wave speed the problem is compounded by the convention of mapping seismic structures as perturbations relative to a 1 d reference model rather than absolute wave speeds using a full waveform tomography model of europe as a case study we quantify the extent to which thermochemical and dynamic properties can be constrained using only s wave speed expressed in absolute values the wave speed distributions of this tomography model are compared with 4 million thermochemical models whose seismic properties are computed via thermodynamic modeling these models sample the full range of realistic mantle compositions including variable water and melt contents and mineral intrinsic anelasticity is taken into account intrinsic anelasticity causes waves to travel more slowly at higher temperatures leading to seismic attenuation but the sensitivity of the wave speed reduction to temperature is in turn controlled by the wave frequency global studies of surface waves indicate an anticorrelation between s wave speed and attenuation we therefore only retain thermochemical models satisfying this anticorrelation our study indicates that the frequency dependence of anelasticity depends on temperature or rheology with 0 1 being most appropriate in cold or lithospheric mantle and 0 3 in warmer regions i e the asthenosphere additionally the slowest regions require specific compositions and or a velocity weakening mechanism such as partial melting elastically accommodated grain boundary sliding or water

the dynamics of back arc extension an experimental approach to the opening of the tyrrhenian sea the e w opening tyrrhenian sea developed after the cretaceous palaeogene alpine collision nearly perpendicular to the motion of the african plate as a back arc of the adria ionian westward subduction three driving mechanisms have been proposed to explain the dynamic evolution of the tyrrhenian apennine system 1 the northward indentation of the african plate 2 the retreating subduction of the adria ionian lithosphere and 3 the gravitational collapse of the alpine post collisional wedge in order to define the relative contribution of each of these mechanisms in the neogene dynamic of the tyrrhenian apennine system we performed 3 d laboratory experiments in which we simulated a retreating subduction process in a compressional regime oriented perpendicularly to the direction of subduction in this framework we also tested the influence of the gravitational collapse of the overriding plate experiments were constructed using dry sand and silicone putties to simulate brittle upper crust and ductile lower crust upper mantle respectively these layers floated on a high density low viscosity glucose syrup which simulated the asthenosphere the main conclusion of our experiments is that large scale continental extension similar to that observed in the tyrrhenian area could be reproduced perpendicular to the shortening direction induced by the indentation of the african plate in this framework extensional processes are indeed possible if the trench retreat velocity is higher than the rate of shortening induced by the advancing african plate our experimental results indicate that this high trench retreat velocity could be explained by the coexistence of the gravitational collapse of the post alpine wedge with a slab pull process linked to the retreating subduction of the adria ionian plate while the first mechanism is predominant in the northern tyrrhenian area the second one seems to be important in the latest stage of extension and oceanic accretion of the southern tyrrhenian area

why did arabia separate from africa insights from 3 d laboratory experiments abstract we have performed 3 d scaled lithospheric experiments to investigate the role of the gravitational force exerted by a subducting slab on the deformation of the subducting plate itself experiments have been constructed using a dense silicone putty plate to simulate a thin viscous lithosphere floating in the middle of a large box filled with glucose syrup simulating the upper mantle we examine three different plate configurations i subduction of a uniform oceanic plate ii subduction of an oceanic continental plate system and iii subduction of a more complex oceanic continental system simulating the asymmetric africa eurasia system each model has been performed with and without the presence of a circular weak zone inside the subducting plate to test the near surface weakening effect of a plume activity our results show that a subducting plate can deform in its interior only if the force distribution varies laterally along the subduction zone i e by the asymmetrical entrance of continental material along the trench in particular extensional deformation of the plate occurs when a portion of the subduction zone is locked by the collisional process the results of this study can be used to analyze the formation of the arabian plate we found that intraplate stresses similar to those that generated the africa arabia break up can be related to the neogene evolution of the northern convergent margin of the african plate where a lateral change from collision mediterranean and bitlis to active subduction makran has been described second intraplate stress and strain localization are favored by the presence of a weakness zone such as the one generated by the afar plume producing a pattern of extensional deformation belts resembling the red sea gulf of aden rift system

the influence of thermodynamic formulation on simulations of subduction zone geometry and history we introduce varying approximations to the thermodynamic properties of a convecting system similar to earth s mantle a realistic geometry rheology and faults are included to insure that the variations we see are robust in a dynamic self buffering environment we find two levels of response in the system first some approximations lead to very different patterns of flow in the simulations second other sets of approximations lead to similar patterns of flow but the time history of the flow varies significantly with the formulation used last we note that only minor changes to an existing code would be needed to integrate a complex thermodynamic model

the somalia plate and the east african rift system present day kinematics summary the motion of the somalia plate relative to the nubia africa arabia and antarctica platcs is re evaluated using a new inversion method based on a monte carlo technique and a least absolute value misfit criterion a subset of the nuvel 1 data set with additional data along the levant fault and in the red sea is used the results confirm that the motion of arabia with respect to africa is significantly different from the motion relative to somalia it is further shown that the data along the sw indian ridge are compatible with a pole of relative motion between africa and somalia located close to the hypothetical diffuse triple junction between the ridge and thc east african rift the resulting africa somalia motion is then compatible with the geological structures and seismological data along the east african rift system assuming a separate somalia plate thus solves kinematic and geological problems around the afar triple junction and along the east african rift

effects of plate bending and fault strength at subduction zones on plate dynamics for subduction to occur plates must bend and slide past overriding plates along fault zones because the lithosphere is strong significant energy is required for this deformation to occur energy that could otherwise be spent deforming the mantle we have developed a finite element representation of a subduction zone in which we parameterize the bending plate and the fault zone using a viscous rheology by increasing the effective viscosity of either the plate or the fault zone we can increase the rates of energy dissipation within these regions and thus decrease the velocity of a plate driven by a given slab buoyancy we have developed a simple physical theory that predicts this slowing by estimating a convecting cell s total energy balance while taking into account the energy required by inelastic deformation of the bending slab and shearing of the fault zone the energy required to bend the slab is proportional to the slab s viscosity and to the cube of the ratio of its thickness to its radius of curvature the distribution of dissipation among the mantle lithosphere and fault zone causes the speed of a plate to depend on its horizontal length scale using the observation that earth s plate velocities are not correlated to plate size we can constrain the lithosphere viscosity to be between 50 and 200 times the mantle viscosity with higher values required if the fault zone can support shear tractions 50 mpa over 300 km these subduction zone strengths imply that the mantle fault zone and lithosphere dissipate about 30 10 and 60 of a descending slab s potential energy release if the slab is 100 km thick the lithospheric component is highly dependent on slab thickness it is smaller for thin plates but may be large enough to prevent bending in slabs that can grow thicker than 100 km ubduction zone strength should be more stable than mantle viscosity to changes in mantle temperature so the controlling influence of subduction zones could serve to stabilize plate velocities over time as the earth cools because the details of convergent plate boundaries are so important to the dynamics of plate motion numerical models of mantle flow should treat subduction zones in a realistic way

cenozoic extension in coastal dhofar southern oman implications on the oblique rifting of the gulf of aden a detailed analysis of fault populations that affect the paleocene to early miocene sedimentary succession of southern dhofar allows the definition of a time sequence of paleostresses and the reconstruction of the rifting history of the eastern gulf of aden two distinct phases of n20j and n160j trending extensions took place during the deposition of the synrift dhofar group from late priabonian 35 ma to middle burdigalian 18 ma a period of tectonic instability preceded these events marked by normal and strike slip faulting as early as the mid eocene the two main directions of extension also recorded in southern yemen developed respectively oblique and normal to the n75e trend of the gulf of aden an anticlockwise succession of the two episodes from n20 to n160 only established in yemen accounts for the observed synrift fault pattern which is consistent with the type and geometry of faults predicted in analogue models of oblique rifting a latest phase of ns compression marked by strike slip faulting may represent a far echo of the arabia eurasia collision d 2002 elsevier science b v all rights reserved

low seismic velocities below mid ocean ridges attenuation versus melt retention the first comprehensive seismic experiment sampling subridge mantle revealed a pronounced low velocity zone between 40 and 100 km depth below the east pacific rise epr that has been attributed to substantial retained melt fractions of 0 3 2 such high melt fractions are at odds with low melt productivity and high melt mobility inferred from petrology and geochemistry here we evaluate whether seismic attenuation can reconcile subridge seismic structure with low melt fractions we start from a dynamic spreading model which includes melt generation and migration and is converted into seismic structure accounting for temperature pressure composition phase and melt dependent anharmonicity and temperature pressure frequency and hydration dependent anelasticity our models predict a double low velocity zone a shallow approximately triangular region due to dry melting and a low velocity channel between 60 and 150 km depth dominantly controlled by solid state high temperature seismic attenuation in a damp mantle with only a minor contribution of 0 1 melt we test how tomographic inversion influences the imaging of our modeled shear velocity features the epr experiment revealed a double low velocity zone but most tomographic studies would only resolve the deeper velocity minimum experimentally constrained anelasticity formulations produce vsas low as observed and can explain lateral variations in near ridge asthenospheric vs with 100 k temperature variations and or zero to high water content furthermore such qs formulations also reproduce low asthenospheric vs below older oceans and continents from basic lithospheric cooling models although these structures are compatible with global qs images they are more attenuating than permitted by epr data

on the origin of the cocos nazca spreading center to investigate the mechanism underlying the break up of the farallon plate into the cocos plate and nazca plate we analyze the state of stress in the farallon plate at about 30 m y b p just prior to the fragmentation to this purpose we use finite element methods and a reconstruction of regional plate boundaries appropriate for 30 m y b p a key role in the model is played by the dependence of two important plate tectonics forces slab pull and ridge push on the age of oceanic lithosphere the results show a highly tensional stress field with maximum principal stresses of 5 to 6 kbar north south tension in the vicinity of present day panama is proposed to have been the cause for the fragmentation of the farallon plate and the inception of spreading along the new cocos nazca plate boundary because the kind of plate interaction that gave rise to the high level of tensional stresses in the farallon plate is not restricted to the area of this study the proposed mechanism seems to shed light on the problem of fragmentation of oceanic plates in general

strain transfer between disconnected propagating rifts in afar we showed before that both the aden and red sea plate boundaries are currently rifting and propagating along two distinct paths into afar through the opening of a series of disconnected propagating rifts here we use new geochronological tectonic and paleomagnetic data that we acquired mostly in the southeastern part of afar to examine the geometry kinematics and time space evolution of faulting related to strain transfer processes it appears that transfer of strain is accommodated by a bookshelf faulting mechanism wherever rifts or plate boundaries happen to overlap without being connected this mechanism implies the rotation about a vertical axis of small rigid blocks along rift parallel faults that are shown to slip with a left lateral component which is as important as their normal component of slip rates of 2 3 mm yr by contrast where rifts do not overlap either a classic transform fault maskali or an oblique transfer zone mak arrasou kinematically connects them the length of the aden red sea overlap has increased in the last 0 9 myr as the aden plate boundary propagated northward into afar as a consequence the first order blocks that we identify within the overlap did not all rotate during the same time span nor by the same amounts similarly the major faults that bound them did not necessarily initiate and grow as their neighboring faults did despite these variations in strain distribution and kinematics the overlap kept accommodating a constant amount of strain 7 to 15 of the extension amount imposed by plate driving forces which remained distributed on a limited number seven or eight of major faults each one having slipped at constant rates 3 and 2 mm yr for vertical and lateral rates respectively the fault propagation rates and the block rotation rates that we either measure or deduce are so fast 30 130 mm yr and 15 38 myr respectively that they imply that strain transfer processes are transient as has been shown to be the case for the processes of tearing rift propagation and strain jumps in afar

radial profile of mantle viscosity results from the joint inversion of convection and postglacial rebound observables we present new inferences of the radial profile of mantle viscosity that simultaneously fit long wavelength free air gravity harmonics associated with mantle convection and a large set of decay times estimated from the postglacial uplift of sites within previously glaciated regions hudson bay arctic canada and fennoscandia the relative sea level variation at these latter sites is constrained by age height pairs obtained by geological survey rather than the subjective trends which are commonly used in glacial isostatic adjustment gia studies our viscosity inferences are generated using two approaches first we adopt a relative viscosity profile which is known to provide a good fit to the free air gravity harmonics and determine an absolute scaling which yields a best fit to the gia decay time constraints second we perform an iterative nonlinear joint inversion of the two data sets in both cases our inferred profiles are characterized by a significant increase of viscosity 2 orders of magnitude with depth to values of 1022 pa s in the bottom half of the lower mantle the new viscosity profiles are shown to satisfy constraints based on the postglacial uplift of both fennoscandia the classic haskell 1935 number and hudson bay which have commonly been invoked to argue for an isoviscous mantle furthermore the models are used to predict a set of long wavelength signatures of the gia process these include predictions of gia induced variations in 1 the length of day over the late holocene period 2 the earth s precession constant and obliquity over the last 2 6 myr and 3 the present day zonal harmonics of the geopotential jl l 7 the predictions 1 and 3 bound the late holocene and ongoing mass flux between the large polar ice sheets greenland and antarctic and the global oceans to small values 0 4 mm yr equivalent eustatic sea level rise

anisotropic tomography around la reunion island from rayleigh waves in the western indian ocean the reunion hot spot is one of the most active volcanoes on earth temporal interactions between ridges and plumes have shaped the structure of the zone this study investigates the mantle structure using data from the reunion hotspot and upper mantle reunions unterer mantel rhum rum project which significantly increased the seismic coverage of the western part of the indian ocean for more than 1 year 57 ocean bottom seismometer stations and 23 temporary land stations were deployed over this area for each earthquake station path the rayleigh wave fundamental mode phase velocities were measured for the periods from 30 s to 300 s and group velocities for the period from 16 s to 250 s a three dimensional model of the shear velocity in the upper mantle was built in two steps first the path mean phase and group velocities were inverted to obtain regionalized velocity maps for each separate period then all of the phase and group velocity maps were combined and inverted at each grid point to obtain the local s wave velocity as a function of depth using a transdimensional inversion scheme the three dimensional anisotropic s wave velocity model has resolution down to 300 km in depth the tomographic model and surface tectonics are correlated down to 100 km in depth starting at 50 km in depth a slow velocity anomaly beneath rodrigues ridge and the east west orientation of the azimuthal anisotropy show a connection between reunion upwelling and the central indian ridge the slow velocity signature beneath la reunion is connected at greater depths 150 300 km with a large slow velocity zone beneath the entire mascarene basin this develops along a northeast direction following the general motion direction of the african plate these observations indicate nonisotropic spreading of hot plume material and dominant horizontal flow in the upper mantle beneath this area

focal mechanisms and depths of earthquakes in central pakistan a tectonic interpretation focal mechanisms and depths for seven earthquakes in central pakistan were determined from an analysis of rayleigh waves of 20 to 50 s periods in east central pakistan the nodal planes for some solutions strike obliquely to the grain of surface structures this observation supports the contention of other workers that a thin surficial unit decoupled from the basement along a surface of decollement characterizes this region in west central pakistan relative movement between the indian and eurasian plates is at least partially accommodated by seismic slip along the chaman fault other faults which are situated to the east of and lie subparallel to the chaman fault may also take up some of the relative plate motion observed activity within the zone of convergent type structure in central pakistan may be a result of the greater component of convergence across the indian eurasian plate boundary north of quetta pakistan appendix is available with entire article on microfiche order from american geophysical union 2000 florida avenue n w washington d c 20009 document b83 010 2 50 payment must accompany order

mantle conveyor beneath the tethyan collisional belt abstract collisional belts are generated by the arrival of continental lithosphere into a subduction zone the tethyan suture from the bitlis to the himalayas is a prime example where the arabian and indian plates collided with eurasia during the cenozoic while the kinematics of this process are well established its dynamics are more uncertain india and arabia intriguingly keep advancing in spite of large collisional resisting forces and in the absence of a substantial upper mantle slab driving force at present day we perform global mantle circulation computations to test the role of deep mantle flow as a driving force for the kinematics of the tethyan collisional belt evaluating different boundary conditions and mantle density distributions as inferred from seismic tomography or slab models our results show that mantle drag exerted on the base of the lithosphere by a large scale convective conveyor belt with an active upwelling component is likely the main cause for the ongoing indentation of the indian and arabian plates into eurasia

paleomagnetic results along the bitlis zagros suture zone in se anatolia turkey implications for the activation of the dead sea fault zone abstract the tectonic evolution of se anatolia is associated with the collision between the taurides and the arabian platform in the early cenozoic after the final closure of the southern branch of the neotethys ocean during late cretaceous to early eocene times the ongoing deformation was characterized by northward movement of arabia followed by westward extrusion of anatolia having been displaced along the northern and eastern anatolian transform faults the neotectonic deformation history between the anatolian plate and the arabian platform since the closure of the southern neotethyan oceanic basin from late cretaceous to present has not been studied in detail paleomagnetically therefore we carried out a paleomagnetic study of mesozoic and cenozoic rocks at 34 localities along the bitlis zagros suture zone in se anatolia located in the boundary region between the taurides and the arabian platform and delimited by the dead sea fault paleomagnetic results from the area show clockwise rotations of 35 6 14 1 and 26 7 8 0 with respect to the african plate from upper jurassic lower cretaceous and upper cretaceous limestones respectively the paleomagnetic results indicate also an earlier deformation opening due to the dead sea fault zone dsfz farther north on the tauride block counterclockwise rotation of 34 5 12 5 is recorded by upper cretaceous limestones in the se of the study area on the northern margin of arabia however smaller counterclockwise rotations of 12 3 8 3 and 2 5 10 0 in late cretaceous and 12 7 13 6 in middle eocene rocks show relative tectonic stability middle miocene rocks close to the suture zone show counterclockwise rotation by about 39 4 10 9 while farther away from the suture zone counterclockwise rotations of 21 9 9 9 and 7 5 5 8 are observed paleolatitudinal motion from south to north from the late jurassic to the present is in concordance with the apparent polar wander path of eurasia and africa

differential rotation between lithosphere and mantle a consequence of lateral mantle viscosity variations the description of plate motions in the so called hotspot reference frame introduces a global rotation of the lithosphere with respect to the mantle this rotation called toroidal field of degree 1 is roughly westward it reaches an amplitude of about 2 cm yr and has been consistently found in the different generation of plate tectonic models various authors have tried to relate this observation to the deceleration of the earth s rotation to polar wander or to tidal drag however these different physical mechanisms cannot explain the requested amplitude in this paper we compare the values of this rotation vector using different relative plate motion models expressed in the hotspot reference frame in a model earth with lateral viscosity variations a differential rotation is predicted the observed net lithospheric rotation is consistent with the dynamics of a model earth where the asthenospheric viscosity below the oceans is at least one order of magnitude lower than underneath the continents this relative westward drift of the lithosphere may account for the significant structural differences between east or west dipping subduction zones

a thin viscous sheet model for continental deformation summary for the purposes of describing its large scale and long term deformation the continental lithosphere is regarded as a continuum obeying a newtonian or a power law rheology the flow of a thin sheet of power law material overlying an inviscid substrate is studied under the assumption that vertical gradients of the horizontal velocity are negligible a numerical model is used to investigate the deformation of such a sheet under conditions approximating those of continent continent collision the material flows in response to forces applied to its boundaries for example the indenting of one continent by another and to forces in its interior arising from gradients in crustal thickness the horizontal divergence of the flow produces changes in the crustal thickness and hence a time dependent form to the flow itself for a given set of boundary conditions the flow depends on the stress exponent in the power law rheology n and on the argand number which is a measure of the ratio between the stress arising from crustal thickness contrasts and the stress required to deform the material at the ambient strain rates when the effective viscosity of the medium is very high ar o crustal thickness variations do not influence the flow if the material is newtonian n l the deformation associated with an influx of material approximating an indenter is of much greater lateral dimension than the width of the indenter whereas when material has a power law rheology n 3 5 are used the deformation is confined to a region of lateral extent comparable to that of the indenter as ar increases the forces arising from crustal thickness contrasts exert more influence on the flow and the maximum crustal thickness that can be sustained by a given influx of material is related by a simple expression to the effective viscosity of the medium at the ambient strain rates in the limit of a very weak medium ar 10 the lithosphere is unable to sustain appreciable crustal elevation contrasts the results of these numerical experiments show that systems in which the effective viscosities are such that the maximum deviatoric stresses are between 1 kbar

hotspot tracks and the early rifting of the atlantic abstract many hotspot tracks appear to become the locus of later rifting as though the heat of the hotspot weakens the lithosphere and tens of millions of years later the continents are split along these weakened lines examples are the west coast of greenland east coast of labrador madeira hotspot the south coast of mexico north coast of honduras guyana hotspot and the south coast of west africa north coast of brazil st helena hotspot a modern day analog of a possible future rift is the snake river plain where the north american continent is being pre weakened by the yellowstone hotspot track this conclusion is based on reconstructions of the motions of the continents over hotspots for the past 200 million years the relative motions of the plates are determined from magnetic anomaly isochrons in the oceans and the motion of one plate is chosen ad hoc to best fit the motions of the plates over the hotspots however once the motion of this one plate is chosen the motions of all the other plates are prescribed by the relative motion constraints in addition to the correlation between the predicted tracks and sites of later continental breakup exposed continental shields correlate with the tracks their exposure may be the result of hotspot induced uplift which has led to erosion of their former platform sediment cover

upper mantle seismic wave velocity effects of realistic partial melt geometries we investigate seismic wave velocity reduction resulting from the presence of partial melt in the upper mantle the amount of shear and bulk modulus reduction produced by the presence of a connected network of realistically shaped and naturally organized melt inclusions is found using finite element calculations the geometries of the inclusions are taken directly from laboratory experiments of mantle melting with finite element meshes constructed to conform to these shapes the shear and bulk moduli of the composite material are found for both the unrelaxed isolated inclusions and relaxed pressure equalized inclusions cases by assigning appropriate material properties to the fluid modulus reduction from deformation simulations of a solid containing realistically shaped and ellipse shaped melt inclusions quantify the effect of melt pocket cuspateness and melt pocket organization on seismic velocity reduction the three dimensional response is estimated from two dimensional distributions of the melt phase by determining the mode ii and mode iii components of elastic modulus reduction separately and summing their effects in general cuspate and naturally organized melt inclusions cause greater velocity reduction it is shown that v and vs reduction per percent partial melt are at least 3 6 and 7 9 respectively even higher values for velocity reduction are possible above 1 melt fraction if melt exists only in tubules below 1 melt fraction the lower more conservative values of velocity reduction are 070 greater for vp and 84 greater for vs than the analytically determined values for ellipsoidal inclusions somewhat greater effects are possible if nonrandom organization of melt occurs on scales greater than our model

the development of slabs in the upper mantle insights from numerical and laboratory experiments we have performed numerical and laboratory experiments to model subduction of oceanic lithosphere in the upper mantle from its beginnings as a gravitational instability to the fully developed slab a two dimensional finite element code is applied to model newtonian creep in the numerical experiments scaled analog media are used in the laboratory a sand mixture models the brittle crust silicone putty simulates creep in the lower crust and mantle lithosphere and glucose syrup is the asthenosphere analog both model approaches show similar results and reproduce first order observations of the subduction process in nature based on density and viscosity heterogeneities in a stokes flow model subduction nucleates slowly and a pronounced slab forms only when the viscosity contrast between oceanic plate and mantle is below a threshold we find that the subduction velocity and angle are time dependent and increase roughly exponentially over tens of millions of years before the slab reaches the 670 km discontinuity the style of subduction is controlled by the prescribed velocity of convergence the density contrast between the plates and the viscosity contrast between the oceanic plate and the mantle these factors can be combined in the buoyancy number f which expresses the ratio between driving slab pull and resisting viscous dissipation in the oceanic plate variations in f control the stress in the plates the speed and the dip of subduction and the rate of trench retreat reproducing the contrasting styles of subduction observed in nature the subduction rate is strongly influenced by the work of bending the lithosphere as it subducts

propagation of rifting along the arabia somalia plate boundary the gulfs of aden and tadjoura the localization and propagation of rifting between arabia and somalia are investigated by assessing the deformation geometry and kinematics at different scales between the eastern gulf of aden and the gulf of tadjoura using bathymetric magnetic seismological and structural evidence large scale southwestward propagation of the aden ridge markedly oblique to the arabia somalia relative motion vector began about 30 myr ago between the error and sharbithat ridges it was an episodic process with stages of rapid propagation mostly at rates 10 cm yr interrupted by million year pauses on transverse discontinuities coinciding with rheological boundaries between different crustal provinces of the arabia somalia plate the longest pause was at the shukra el sheik discontinuity 45 e where the ridge tip stalled for 13 myr between 17 and 4 ma west of that discontinuity rifting and spreading took place at an azimuth n25 10 e and rate 1 2 0 3 cm yr different from those of the global arabia somalia motion vector n39 1 73 cm yr implying an additional component of movement n65 10 e 0 7 0 2 cm yr due to rotation of the danakil microplate at shukra el sheik the typical oceanic ridge gives way to a narrow wsw trending axial trough resembling a large fissure across a shallow shelf this trough is composed of about eight rift segments which result from normal faulting and fissuring along n110 n130 e trends all the segments step to the left southwestward mostly through oblique transfer zones with en echelon normal faults only two segments show clear significant overlap there is one clear transform the maskali fault between the obock and tadjoura segments the latter segment which encroaches onland is composed of two parallel subrifts iboli ambabbo that propagated northwestward and formed in succession the most recent southwestern subrift ambabbo represents the current tip of the aden ridge we propose a mechanical model in which the large scale propagation of the ridge followed a wsw trending zone of maximum tensile stress while the small scale propagation of its nw trending segments was dictated by the orientation of that stress oblique propagation was a consequence of passive lithospheric necking of the arabia somalia plate along its narrow section in map view between socotra and the kink of the red sea ethiopian rift above the afar plume individual ridge segments oriented roughly perpendicular to plate motion like lithospheric cracks were forced to jump southward because of confinement within the necking zone self sustaining plate scale necking may explain why the aden ridge did not connect with the red sea through bab el mandeb but continued straight into afar

displacements strains and rotations within the afar depression djibouti from restoration in map view within the afar depression djibouti the stratiform basaltic lava flows emplaced between 4 ma and 1 ma cover the two thirds of the depression because of the desert climate they are almost non eroded the topography can then be used as a reference surface to measure the vertical offset of younger faults using a detailed map the southeastern part of the depression showing the fault pattern and the vertical offset on each fault obtained from stereoparts of spot images we have restored the topography in map view we have computed the fields of finite displacement rotation and strain the field of horizontal displacements shows vectors trending to the ne on average this direction is compatible with the displacement to the ne of the arabian plate with reference to the african plate in detail however displacement vector trend from nne ssw to ne sw major axes of strain ellipse trend ne sw suggesting that major faults n120 are dominantely extensional strain intensities and magnitudes of block rotations are increasing from southwest to northeast of the restored area magnitudes of block rotations are very small compared with those obtained from paleomagnetic measurements therefore we computed a second restoration that covers a smaller area where the true vertical offset of the major faults taking into account the sedimentary fill of basins were avalaible this second model gives rotation magnitudes in better agreement with paleomagnetic data but still underestimates we discuss these results with the available kinematic models of the southeastern part of the depression

comparison of dike intrusions in an incipient seafloor spreading segment in afar ethiopia seismicity perspectives 1 oceanic crust is accreted through the emplacement of dikes at spreading ridges but the role of dike intrusion in plate boundary deformation during continental rupture remains poorly understood between 2005 and 2009 the 70 km long dabbahu manda hararo rift segment in ethiopia has experienced 14 large volume dike intrusions 9 of which were recorded on temporary seismic arrays a detailed comparison of the seismic characteristics of the seismically monitored dikes is presented with implications for dike intrusion processes and magmatic plumbing systems all of the migrating swarms of earthquakes started from a 80 of energy is released during the propagation phase with minimal seismic energy release after the dike propagation ceases we interpret that faulting and graben formation above the dikes occurs hours after the passage of the dike tip coincident with the onset of low frequency earthquakes dike lengths show no systematic reduction in length with time suggesting that topographic loading and stress barriers influence dike length as well as changes in tectonic stress the propagation velocities of all the dikes follow a decaying exponential northward propagating dikes had faster average velocities than those that propagated southward suggesting preconditioning by the 2005 megadike or ongoing heating from a subcrustal magma source north of the midsegment

coupling of south american and african plate motion and plate deformation although the african plate s northeastward absolute motion slowed abruptly 30 million years ago the south atlantic s spreading velocity has remained roughly constant over the past 80 million years thus requiring a simultaneous westward acceleration of the south american plate this plate velocity correlation occurs because the two plates are coupled to general mantle circulation the deceleration of the african plate due to its collision with the eurasian plate diverts mantle flow westward increasing the net basal driving torque and westward velocity of the south american plate one result of south america s higher plate velocity is the increased cordilleran activity along its western edge beginning at about 30 million years ago

episodic back arc extension during restricted mantle convection in the central mediterranean abstract in the central mediterranean during the last 80 myr africa has been slowly converging toward eurasia with an average relative velocity of less than 1 cm yr 1 for 50 myr this slow convergence led to the build up of the alpine orogenic belt and to the initiation of consumption of the former tethys ocean in the last 30 myr however the whole central western mediterranean has been swept by rapid episodes of trench migration up to 6 cm yr 1 and back arc opening which consumed the whole former tethys ocean we combine plate tectonic history geological timing and laboratory modelling to reconstruct the subduction history and the evolution of the central mediterranean over the last 80 myr we find that the dynamic evolution of the subducting oceanic lithosphere can reconcile the rapid episodic back arc migration with the slow african convergence but only if the subduction process is restricted to the upper mantle

three dimensional full waveform inversion of short period teleseismic wavefields based upon the sem dsm hybrid method we present a method for high resolution imaging of lithospheric structures based on full waveform inversion of teleseismic waveforms we model the propagation of seismic waves using our recently developed direct solution method spectral element method hybrid technique which allows us to simulate the propagation of short period teleseismic waves through a regional 3 d model we implement an iterative quasi newton method based upon the l bfgs algorithm where the gradient of the misfit function is computed using the adjoint state method compared to gradient or conjugate gradient methods the l bfgs algorithm has a much faster convergence rate we illustrate the potential of this method on a synthetic test case that consists of a crustal model with a crustal discontinuity at 25 km depth and a sharp moho jump this model contains short and long wavelength heterogeneities along the lateral and vertical directions the iterative inversion starts from a smooth 1 d model derived from the iasp91 reference earth model we invert both radial and vertical component waveforms starting from long period signals filtered at 10 s and gradually decreasing the cut off period down to 1 25 s this multiscale algorithm quickly converges towards a model that is very close to the true model in contrast to inversions involving short period waveforms only which always get trapped into a local minimum of the cost function

long term elasticity in the continental lithosphere modelling the aden ridge propagation and the anatolian extrusion process summary the evolution of the gulf of aden and the anatolian fault systems are modelled using the principles of elastic fracture mechanics usually applied to smaller scale cracks or faults the lithosphere is treated as a plate and simple boundary conditions are applied that correspond to the known plate boundary geometry and slip vectors the models provide a simple explanation for many observed geological features for the gulf of aden the model predicts why the ridge propagated from east to west from the owen fracture zone towards the afar and the overall form of its path the smaller en echelon offsets can be explained by upward propagation from the initially created mantle dyke while the larger ones may be attributed to the propagating rupture interacting with pre existing structures for anatolia the modelling suggests that the east anatolian fault was created before the north anatolian fault could form once both faults were formed however activity could switch between them the time scales over which this should take place are not known but evidence for switching can be found in the historical seismicity for aden and anatolia pre existing structures or inhomogeneous stress fields left from earlier orogenic events have modified the processes of propagation and without an understanding of the existence of such features the propagation processes cannot be fully understood furthermore a propagating fault can extend into an active region where it would not have initiated the north anatolian fault encountered slow but active extension when it entered the aegean about 5 ma and the stress field associated with the extending fault has progressively modified aegean extension in the central aegean activity has been reduced while to the north west on features such as the gulfs of evvia and corinth activity has been increased the field observation that major structures propagate and the success of simple elastic models suggest that the continental crust behaves in an elastic brittle or elastic plastic fashion even though laboratory tests may be interpreted to suggest viscous behaviour there are major problems in scaling from the behaviour of small homogeneous samples to the large heterogeneous mantle and large scale observations should be treated more seriously than extrapolations of the behaviour of laboratory experiments over many orders of magnitude in space and time the retention of long term elasticity and localised failure suggests a similar gross rheology for the oceanic and continental lithospheres even though it is incorrect to attribute differences in behaviour to the former being rigid i e elastic and the latter viscous oceanic and continental lithosphere behave in different ways unlike oceanic crust continental crust is buoyant and cannot be simply created or destroyed the process of thickening or thinning works against gravity preventing large displacements on extensional or contractional features in the upper mantle the equivalents of ridge or subduction systems are suppressed before they can accommodate large displacements and activity must shift elsewhere on the other hand strike slip boundaries and extrusion processes are favoured

reappraisal of the arabia india somalia triple junction kinematics abstract we propose alternative kinematics for the arabia india somalia triple junction based on a re interpretation of seismological and magnetic data the new triple junction of the ridge ridge ridge type is located at the bend of the sheba ridge in the eastern gulf of aden at 14 5 n and 56 4 e the owen fracture zone arabia india boundary is connected to the sheba ridge by an ultra slow divergent boundary trending n80 e 10 marked by diffuse seismicity the location of the arabia india rotation pole is constrained at 14 1 n and 71 2 e by fitting the active part of the owen fracture zone with a small circle the finite kinematics of the triple junction is inferred from the present day kinematics since the inception of the accretion 15 18 ma ago the sheba ridge has probably receded 300 km at the expense of the carlsberg ridge which propagated northwestward in the gulf of aden while an ultra slow divergent plate boundary developed between the arabian and indian plates the overall geometry of the new triple junction is very similar to that of the azores triple junction

role of ocean continent contrast and continental keels on plate motion net rotation of lithosphere and the geoid three dimensional spherical shell models of mantle convection have been formulated to investigate the effects of ocean continent contrast i e different viscosity for continental and oceanic lithospheres and continental keels on plate motion net rotation of lithosphere i e degree 1 toroidal plate motion and the geoid the models include relatively realistic plate rheology i e strong plate interiors and weak plate margins and continental keel structure based on seismic models the mantle flow in the models is driven by slab buoyancy the models demonstrate that the ocean continent contrast and continental keels have minor effects on the long wavelength geoid for degrees 2 5 however continental keels overriding subducted slabs for example in north america can produce large negative gravity anomalies at a regional scale this may have important implications to interpreting the gravity anomalies in north america our models show that when plates bounded by weak plate margins have the same thickness neither weak plate margins nor the ocean continent contrast efficiently excite net rotation of lithosphere although they excite toroidal motion at higher harmonic degrees however plate thickness variations like continental keels excite the net rotation in spite of their ability to excite the net rotation continental keels as thick as 300 km cannot provide the necessary coupling to the deep mantle to produce the observed net rotation of lithosphere and slow continental motion if there is a weak asthenosphere underlying them if the thermomechanic structure of north american continental upper mantle derived from seismic and heat flow studies is representative our models suggest that the temperature and pressure dependent mantle rheology may not produce sufficiently high viscosity below continental keels that is needed to explain the observed plate motion our study hints the necessity to include more realistic treatment of plate boundaries in future studies to assess possible effects of plate plate coupling on plate motion

receiver function imaging of lithospheric structure and the onset of melting beneath the galapagos archipelago abstract the galapagos archipelago represents an opportunity to investigate the properties of young oceanic lithosphere the effects of a hotspot anomaly on lithospheric thickness and melting dynamics in a hotspot ridge interaction here we use data recorded by the signet array and permanent station payg on the islands santa cruz and isabela respectively we used p to s ps and s to p sp receiver functions to constrain crust and mantle structure a simultaneous deconvolution method was used to constrain 1 d structure and also for the modeling of robust features a migrated extended multitaper method was used to investigate 3 d structural variations ps images a velocity increase with depth at 11 7 km probably the base of the pre plume crust or old moho sp imaging and modeling images a second deeper velocity increase at 37 7 km depth a velocity decrease with depth is imaged on average at 75 12 km likely associated with the lithosphere asthenosphere boundary this discontinuity is imaged deeper 82 km in the southwest and shallower 66 km in the northeast near the spreading ridge although the trend is consistent with lithospheric thickening with age the thickness is much larger than predicted by conductive cooling models of 0 10 my oceanic lithosphere we infer a compositional contribution to velocity variations finally a velocity increase with depth is imaged at 125 to 145 15 km depth that is likely associated with the onset of melting the discontinuity is imaged deeper in 3 sectors of the galapagos platform ridge region all coincident with the slowest surface wave shear velocity anomalies in the upper 100 km one is located in the southwest in a hypothesized plume location the other two are to the northwest and northeast possibly illuminating multiple plume diversions related to complex plume ridge interactions

closing the gap between regional and global travel time tomography recent global travel time tomography studies by zhou 1996 and van der hilst et al 1997 have been performed with cell parameterizations of the order of those frequently used in regional tomography studies i e with cell sizes of 1 2 these new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations 5 cells the inferred structures are however of larger scale than is usually obtained in regional models and it is not clear where and if individual cells are actually resolved this study aims at resolving lateral heterogeneity on scales as small as 0 6 in the upper mantle and 1 2 3 in the lower mantle this allows for the adequate mapping of expected small scale structures induced by for example lithosphere subduction deep mantle upwellings and mid ocean ridges there are three major contributions that allow for this advancement first we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the earth s mantle by seismic waves spakman and bijwaard 1998 second we exploit the global data set of engdahl et al 1998 which is a reprocessed version of the global data set of the international seismological centre their reprocessing included hypocenter redetermination and phase reidentification finally we combine all data used p pp and pwp phases into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal to noise ratio is improved in the approximate solution of the huge inverse problem we obtain a variance reduction of 57 1 synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells 0 6 or 1 2 in the shallow parts of subduction zones decreasing to approximately 2 3 resolution in well sampled regions in the lower mantle vertical resolution can be worse up to several hundreds of kilometers in subduction zones with rays predominantly pointing along dip important features of the solution are as follows 100 200 km thick high velocity slabs beneath all major subduction zones sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle large high velocity anomalies in the lower mantle that have been attributed to subduction of the tethys ocean and the farallon plate and low velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under iceland east africa the canary islands yellowstone and the society islands our findings corroborate that the 660 km boundary may resist but not prevent present day large scale mass transfer from upper to lower mantle or vice versa this observation confirms the results of previous global mantle studies that employed coarser parameterizations

building the zagros collisional orogen timing strain distribution and the dynamics of arabia eurasia plate convergence he zagros mountains are the result of the arabia eurasia collision initiated at 35 ma as the rifted arabian lithosphere was underthrusted beneath the iranian plate due to its negative buoyancy the onset of crustal thickening started at 25 ma as recorded by the hinterland exhumation and foreland clastic deposition deformation throughout the arabia eurasia collision zone and the uplift of the iranian plateau occurred after 15 12 ma as a result of shortening thickening of the thin iranian crust we emphasize that only 42 of the post 35 ma convergence is partitioned by shortening within central iran tomographic constraints show ongoing slab steepening or breakoff in the nw zagros whereas underthrusting of the arabian plate is observed beneath central zagros the current subduction dynamics can be explained by the original lateral difference in the buoyancy of the distal margin that promoted slab sinking in nw zagros and underthrusting in central zagros critical wedge approach applied to the zagros favors the hypothesis of strong brittle crust detached above a viscous lower crust in contrast the weak sedimentary cover deforms by buckling of a thick multilayered cover thrust faulting associated with folding occurs in the competent layers and is responsible for most of the earthquakes there is evidence that the role of the slab pull force in driving the arabian plate motion was reduced after 12 ma large scale mantle flow induced by mantle upwelling at the afar plume appears to be the main driver of the arabia plate motion we stress that the main kinematic change in the zagros region occurred at 15 12 ma as the zagros uplifted before the arabian slab detached the zagros appears key to investigate coupling between continental rheology plate driving forces and mountain building in which the role of rift inheritance appears to be central

effects of mantle heat source distribution on supercontinent stability a two dimensional cartesian geometry numerical model is used to study the influence of internal heating on the ability of mantle convection to rift and disperse an assembled supercontinent our wide aspect ratio models incorporate migrating plate boundaries and finite thickness rigid mobile oceanic and continental plates 12 000 km across we examine models characterized by both whole mantle and upper mantle convection parameters specifying different degrees of internal heating in the mantle we compare the patterns of convection that develop below stationary continental plates that form as a result of the aggregation of smaller mobile continents following continental collisions we find that a reorganization of mantle flow may result in the development of deviatoric stresses at the base of the continental lithosphere which are capable of initiating continental breakup our results suggest that continental basal stresses are largely controlled by the wavelength of the subcontinental convection pattern and that the inclusion of internal heating in the mantle generally increases this wavelength for models with 80 internal heating in which active upwellings are intrinsically absent basal stresses arise primarily from subduction at the continental margins and result in a geologically reasonable duration for our model supercontinent assemblages 200 myr for models with 40 or less internal heating continental breakup is controlled by mantle upwellings and our supercontinents remain intact much longer 600 700 myr

subducted slabs and the geoid constraints on mantle rheology and flow geoid anomalies are primarily the result of the density contrasts driving mantle convection and plate motion the total geoid anomaly resulting from a given density contrast in a convecting earth is affected by the mass anomalies associated with the flow induced deformation of the upper surface and internal compositional boundaries as well as by the density contrast itself these boundary deformations and hence the total gravity field depend on the radial distribution of effective viscosity if the internal density contrasts can be estimated as is the case for subducted slabs useful constraints can be placed on the depth and on the variation of viscosity with depth of the convecting system the degree 4 9 components of the observed long wavelength geoid are highly correlated with those predicted by a density model for seismically active subducted slabs the positive sign of the correlation requires that the effective viscosity increases with depth by a factor of 30 or more the amplitude of the correlation cannot be explained by the density contrasts associated with just the seismically active parts of subducted slabs however the amplitude can be explained if the density contrasts associated with subduction extend into the lower mantle or if old lithosphere is piled up at the base of the upper mantle beneath subduction zones to a thickness in excess of 350 km over horizontal distances of thousands of kilometers mantlewide convection in a mantle that has a viscosity increasing with depth provides a simple explanation of the long wavelength geoid anomalies over subduction zones

tectonic setting origin and obduction of the oman ophiolite the semail ophiolite in the oman mountains is the world9s largest and best preserved thrust sheet of oceanic crust and upper mantle 10 000 km 2 550 km long 150 km wide it was emplaced onto the arabian continental margin during late cretaceous time the ophiolite originated 96 94 ma at a spreading center above a northeast dipping subduction zone associated with initiation of immature island arc tholeiitic lavas lasail arc at the highest levels of the ophiolite simultaneous underthrusting of triassic and jurassic mid oceanic ridge basalt and alkalic volcanic rocks beneath 12 km of upper mantle depleted harzburgites produced garnet clinopyroxene amphibolites formed at temperatures of 850 c dated as 95 93 ma subduction cannot have been initiated at a mid oceanic ridge otherwise the protolith of the amphibolites in the metamorphic sole would be the same age and composition as the ophiolite volcanic rocks above in the northern part of the oman mountains in the bani hamid area united arab emirates 870 m of granulite facies rocks enstatite spinel diopside quartzites garnet diopside wollastonite calc silicate marbles clinopyroxene bearing amphibolites were formed at temperatures similar to those of the garnet diopside amphibolites of the oman sole 800 850 c but at slightly higher pressures as much as 9 kbar they are interpreted as deeper level metamorphosed continental margin sedimentary rocks exhumed by out of sequence thrusting placing granulites over mantle sequence harzburgites during the later stages of obduction subduction of the arabian continental crust beneath the obducting semail ophiolite to 78 90 km depth has been proven by thermobarometry of the as sifah eclogites to 20 23 kbar in the eastern sector in the united arab emirates the subducted continental crust began to partially melt producing unusual biotite muscovite garnet tourmaline cordierite andalusite bearing granites that intrude the uppermost mantle sequence harzburgites and lowermost crustal sequence cumulate gabbros of the ophiolite we suggest that the entire leading northeast edge of the arabian plate was subducted beneath the ophiolite during the final stages of obduction leading to eclogitization of the crustal rocks higher temperatures and pressures in the united arab emirates sector possibly due to a thicker or double thickness ophiolite section led to blueschist amphibolite and granulite facies conditions in the metamorphic sole and crustal melting in the subophiolite basement produced leucocratic granites that intruded up as dikes through the obducted ophiolite a model for ophiolite obduction is presented which accounts for all the structural and metamorphic conditions reported from the oman mountains

propagation of an oblique spreading centre the western gulf of aden abstract the east west trending gulf of aden rift between arabia and somalia is spreading along a n35 e direction it propagates westward inside the african continent toward the afar triple junction oblique to the spreading direction these combined propagation and oblique rifting generate a special faulting pattern associated with an irregular distribution of volcanism we analyse this pattern using a structural map derived from the analysis of multibeam bathymetric data and backscattering images from the 1995 tadjouraden cruise in the western gulf of aden the rift valley is bounded by faults trending from n110 e to n130 e in the eastern part of the surveyed area where the valley trends n90 e and by large faults in the central western part where the valley is oriented n70 e the direction of extension deduced from the fault histogram is n37 e 10 over the whole area whatever the rift trend therefore this bend of the plate boundary does not result from a change in the stretching direction this bending increases the rifting obliquity westward and therefore modifies the faulting pattern faults display two distinct patterns either en echelon on the walls indicating a right lateral component of shear or purely normal and parallel to the rift trend inside the valley volcanoes are mostly located in the eastern part of the rift where they are aligned into swarms trending roughly perpendicular to the extension direction deduced from fault directions backscattering images show that recent lava flows fill the axial valley of the eastern part whereas in the central and western parts volcanoes are scattered in the middle part of the valley this distribution is interpreted as a westward propagation of the volcanism and as the direction of deformation migration the thickness of the upper brittle layer of the lithosphere is estimated from the width if the inner graben is inside the rift it indicates a deepening of the brittle ductile transition from east 6 km to west 16 km these values are typical of oceanic and continental systems respectively and consistent with the depth of earthquakes in the western part this points out the existence of a transitional area between 44 e and 45 20 e characterised by the lack of recent volcanism the nature of this mechanical transition is discussed in terms of propagation of an spreading centre into a continental or transitional afar type lithosphere

no net rotation model of current plate velocities incorporating plate motion model nuvel 1 we present nnr nuvel1 a model of plate velocities relative to the unique reference frame defined by requiring no net rotation of the lithosphere while constraining relative plate velocities to equal those in global plate motion model nuvel 1 demets et al 1990 differences between nnr nuvel1 and no net rotation plate motion model am0 2 minster and jordan 1978 are as large as 15 mm yr in nnr nuvel1 the pacific plate rotates in a right handed sense relative to the no net rotation reference frame at 0 67 m y about 63 s 107 e this rotation nearly parallels but is 32 slower than the pacific hotspot euler vector of hs2 nuvel1 which is a global model of plate hotspot velocities constrained to consistency with nuvel 1 gripp and gordon 1990 at hawaii the pacific plate moves relative to the no net rotation reference frame at 70 mm yr which is 25 mm yr slower than the pacific plate moves relative to the hotspots differences between nnr nuvel1 and hs2 nuvel1 are described completely by a right handed rotation of 0 33 m y about 49 s 65 e the 99 confidence ellipsoids for plate hotspot motion in model hs2 nuvel1 exclude the corresponding euler vectors from model nnr nuvell thus the no net rotation reference frame differs significantly from the hotspot reference frame if the difference between reference frames is caused by motion of the hotspots relative to a mean mantle reference frame then hotspots beneath the pacific plate move with coherent motion towards the east southeast alternatively the difference between reference frames may show that the uniform drag no net torque reference frame which is kinematically equivalent to the no net rotation reference frame is based on a dynamically incorrect premise possible exceptions to the assumption of uniform drag include a net torque on the lithosphere due to attached subducting slabs and greater resistance to plate motion beneath continents than beneath oceans

how preexisting weaknesses control the style of continental breakup continental breakup is a highly variable process differences occur in the relative timing and extent of volcanism uplift and graben formation as well as in the mode and amount of continental extension before breakup we propose a model which reconciles this variability with the previously recognized tendency for breakup to occur along preexisting weak trends continental lithosphere is viewed as a composite material composed of two strong layers one in the upper mantle and one in the middle crust finite element simulation indicates that extensional failure at weaknesses in the mantle causes concentrated extension in the mantle and diffuse extension in the crust this leads to early melt segregation and volcanism margin uplift during the late stages of the extension process and relatively narrow symmetrical extended margins in contrast failure at weaknesses in the crustal strong zone causes focused extension in the crust and diffuse extension in the mantle this produces initial graben formation cooling in the lower crust and upper mantle and broad asymmetrical extended margins volcanism only occurs late in the process failure at laterally offset weaknesses within both strong layers leads to a deformation pattern dominated simple shear thus differences in the prerift configuration of the continental lithosphere can control the overall style of continental breakup

thickening of young pacific lithosphere from high resolution rayleigh wave tomography a test of the conductive cooling model direct seismic measurements of the thickening of oceanic lithosphere away from the spreading axis are rare due to the remoteness of mid ocean ridges on the east pacific rise at 17 s there were two long term broadband ocean bottom seismometer deployments the melt and glimpse experiments these arrays spanned young pacific plate seafloor ranging in age from 0 to 8 ma combining these two data sets we describe the increase in rayleigh wave phase velocities for 16 100 s period as function of distance from the ridge with two parameterizations an arbitrary function of seafloor age and a simple polynomial function of age on the pacific and nazca plates although resolution analysis shows that 10 to 15 independent pieces of information about the age variation of phase velocity on the pacific plate can be resolved at periods of 50 s the three parameter polynomial model fits the data almost as well indicating a simple pattern of the evolution of structure with age to compare our observations to the predictions for the conductive cooling of the lithosphere we convert a thermal half space cooling model to shear velocity using anharmonic and anelastic contributions the patterns in the velocities are not consistent with simple conductive cooling the conductive cooling model under predicts the changes in phase velocity with age observed in the 20 to 78 s period range we observe significant changes in shear velocity in the asthenosphere at depths greater than conductive cooling should extend the asthenospheric low velocity zone is asymmetric dipping to the west with the lowest velocities beneath the pacific plate west of the spreading center the conductive cooling model also over predicts the shear velocities in the lithosphere by 0 2 km s these anomalies indicate that more complicated mantle flow and significant mantle heterogeneities such as melt are required

the rotation of arabia and the levant fault system abstract we establish the kinematics of the rotation of arabia with respect to africa and show that the rotation rate of africa increased by a factor of about four in late serravallian time about 12 13 ma ago we discuss evidence that oceanic accretion started simultaneously in the gulf of aden and in the red sea at this time it is also at this time that the levant and anatolian fault systems were initiated our initial reconstruction of the red sea rift in late oligocene early miocene time approximately superposes the present shore lines which coincide with the mid point of the zones of rapid crustal thinning the future continental slopes the amount of extension during the early phase of continental rifting is obtained by assuming that the total volume of continental crust is conserved during extension and by checking these estimates through closure around the sinai triple junction in the late miocene the southern evant fault zone was a zone of distributed shear deformation through counterclockwise block rotation extending soul h ward within the northernmost red sea since the early pliocene the shear motion has been localized within the gulf of aqaba changing the tectonic pattern in the northern red sea

the influence of trench migration on slab penetration into the lower mantle abstract a two dimensional numerical convection model in cartesian geometry is used to study the influence of trench migration on the ability of subducted slabs to penetrate an endothermic phase boundary at 660 km depth the transient subduction history of an oceanic plate is modelled by imposing plate and trench motion at the surface the viscosity depends on temperature and depth a variety of styles of slab behaviour is found depending predominantly on the trench velocity when trench retreat is faster than 2 4 cm a the descending slab flattens above the phase boundary at slower rates it penetrates straight into the lower mantle although flattening in the transition zone may occur later leading to a complex slab morphology the slab can buckle independent of whether it penetrates or not especially when there is a localised increase in viscosity at the phase boundary flattened slabs are only temporarily arrested in the transition zone and sink ultimately into the lower mantle the results offer a framework for understanding the variety in slab geometry revealed by seismic tomography

trench motion slab geometry and viscous stresses in subduction systems summary a semi analytic 3 d model for subduction within a newtonian viscous upper mantle provides a dynamically consistent means of computing viscous stress trench motion and slab geometry in subduction systems although negative slab buoyancy provides the basic driving force for subduction slabs that extend from the surface to the base of the upper mantle are oversupported by viscous stresses in the shallow 100 150 km mantle and undersupported by viscous stresses at greater depth in the upper mantle these deeper parts of the subduction system act as an engine for subduction while shallower parts act as a brake on trench motion trench migration rates and slab geometry reflect a competition between these two effects during steady state subduction trench migration rates vary approximately linearly with slab buoyancy and model rates of trench motion are in good agreement with the range of observed rates for a two layer upper mantle viscosity of 2 10 20 pa s above 300 km and 5 10 20 pa s below steady state slab dip increases as slab density decreases especially for very low density slabs which dip significantly more steeply than high density slabs the horizontal velocity at the top of the lower mantle measured relative to the foreland has a very large effect on trench migration rates rivalling or even exceeding that of slab buoyancy slab width parallel to the trench also has a significant effect on trench migration rates due to the viscous pressure of toroidal flow around the slab the stiffness of the subducting lithosphere does not exert a significant effect on trench migration rates or slab geometries for rigidities compatible with oceanic lithosphere very stiff slabs with elastic plate thicknesses more that 40 km or viscosities in the range of 10 25 10 26 pa s subduct significantly more slowly than weak slabs with trench migration rates in the range of half to a third that of weak slabs large unexpected effects on trench migration rates and slab geometry are exerted by the structure and density of the frontal prism and overriding plate indicating that local geology can exert important constraints on slab dynamics during non steady state subduction rates of trench migration respond rapidly as variably buoyant lithosphere penetrates into the asthenosphere in the absence of other driving forces for convergence trench migration rates can change by a factor of two or more in as little as 2 3 myr for example from 35 to 70 mm yr 1 when an oceanic piece of slab follows a continental one into the subduction system subduction of variable buoyancy lithosphere is accompanied by changes in slab dip with depth and through time

a methodology for inversion of broadband seismic waveforms for elastic and anelastic structure and its application to the mantle transition zone beneath the northwestern pacific abstract we develop a quantitative and objective methodology for inversion of body wave waveform data for elastic and anelastic structure and apply it to the upper mantle and mantle transition zone mtz beneath the northwestern pacific using prem as the starting model we use the singular value decomposition svd as the basis for our inversion and use akaike s information criterion aic to evaluate the statistical significance of the models although the resulting models are preliminary we obtain an average positive s wave velocity anomaly in the mtz beneath the philippine sea plate and the japan sea we also find a negative s wave velocity anomaly in the depth range 300 400 km beneath the subduction zone to the east of northern japan and its southern extension beneath the izu bonin marianas arc our inversions obtained a lower q as compared to prem for all of the subregions

does partial melting explain geophysical anomalies abstract the existence of partial melt is frequently invoked to explain geophysical anomalies such as low seismic wave velocity and high electrical conductivity i review various experimental and theoretical studies to evaluate the plausibility of this explanation in order for a partial melt model to work not only the presence of melt but also the presence of appropriate amount of melt needs to be explained using the mineral physics observations on the influence of melt on physical properties and the physics and chemistry of melt generation and transport i conclude that partial melt model for the asthenosphere with homogeneous melt distribution does not work one needs to invoke inhomogeneous distribution of melt if one wishes to explain observed geophysical anomalies by partial melting however most of models with inhomogeneous melt distribution are either inconsistent with some geophysical observations or the assumed structures are geodynamically unstable and or implausible therefore partial melt models for the geophysical anomalies of the asthenosphere are unlikely to be valid and some solid state mechanisms must be invoked the situation is different in the deep upper mantle where melt could completely wet grain boundaries and continuous production of melt is likely by dehydration melting at around 410 km in the ultralow velocity zone in the d layer where continuous production of melt is unlikely easy separation of melt from solid precludes the partial melt model for low velocities and high electrical conductivity unless the melt density is extremely close to the density of co existing solid minerals or if there is a strong convective current to support the topography of the ulvz region compositional variation such as fe enrichment is an alternative cause for the anomalies in the d layer

spatial distribution of hotspot material added to the lithosphere under la reunion from wide angle seismic data wide angle seismic lines recorded by ocean bottom and land seismometers provide a pseudo three dimensional investigation of the crust and upper mantle structure around the volcanically active hotspot island of la reunion the submarine part of the edifice has fairly low seismic velocities without evidence for intrusives an upper unit with a velocity depth gradient is interpreted as made of material erupted subaerially then transported and compacted downslope between this unit and the top of the oceanic plate imaged by normal incidence seismic reflection a more homogeneous unit indicated by shadow zones on several wide angle sections may correspond to lavas of a different nature extruded underwater in the earlier phase of volcanism coincident wide angle and normal incidence reflections document that the oceanic plate is not generally downwarping toward the island but doming instead toward its southeastern part with limited evidence for some intracrustal intrusion deeper in the lithosphere the presence of a layer of intermediate velocity between the crust and mantle is firmly established it is interpreted as resulting from the advection of hotspot magmatic products possibly partially molten and of a composition for which the crust is a density barrier the extensive wide angle coverage constrains the extent of this body it does not show the elongated shape expected from plate drift above a steady hotspot supply alternative propositions can hence be considered for example that la reunion is caused by a solitary wave of hotspot material or by a young hotspot the size of the underplate 140 km wide and up to 3 km thick corresponds to less than half the volume of the edifice on top of the plate

structure and deformation history of the bitlis suture near lake hazar southeastern turkey detailed mapping of the bitlis suture zone near lake hazar shows that it is composed of three tectonostratigraphic units which crop out as distinct north dipping thrust slices from south to north these are the puturge metamorphic complex the maden melange and the elazig igneous complex the puturge metamorphic complex consists of pre tertiary continental margin sediments metamorphosed to the greenschist facies during the campanian maastrichtian and deformed by four generations of structures these are in sequence 1 isoclinal folds and a transposition foliation 2 open folds crenulation cleavage and north dipping thrust faults 3 kink bands and 4 small displacement faults the maden melange represents middle eocene back arc basin sediments and volcanics metamorphosed to the greenschist facies and deformed by three generations of structures 1 a north dipping cleavage 2 kink bands and 3 small displacement faults the elazig igneous complex comprises an imbricated maastrichtian early eocene island arc and young marginal basin terrain thrust faults between units are north dipping listric and form a thin skinned system these features suggest the following deformational and tectonic history the puturge metamorphic complex was generated by isoclinal folding and metamorphism of sediments composing the arabian continental margin during campanian maastrichtian ophiolite obduction to the south the elazig arc developed on the deformed margin because of subsequent southward subduction of oceanic lithosphere the arc migrated to the north opening a back arc basin by the early eocene the back arc basin was filling with volcaniclastics of the maden melange and the elazig arc collided with a continent to the north during middle to late eocene convergence caused thrust stacking of the elazig maden and puturge complexes second generation deformation in the puturge complex and first generation deformation in the maden melange continued convergence with the arabian plate in the late miocene caused kink banding in the puturge and maden complexes since the late miocene this convergence has been accommodated by shortening and thickening along numerous internal faults

paleogeographic maps of the continents bordering the indian ocean since the early jurassic paleomagnetic data from the african plate from the early jurassic onward are critically re evaluated and updated resulting in a set of 30 paleopoles from which an apparent polar wander path is derived using 20 ma intervals this path is then used in conjunction with kinematic models of the indian and atlantic oceans to place other continents bordering the indian ocean eurasia india islands and plateaus madagascar seychelles mascarene and blocks involved in the collision of india with eurasia south tibet indochina in an absolute paleogeographic frame fixed with respect to africa paleogeographic maps are presented for the early jurassic pre breakup the jurassic cretaceous boundary anomaly ml5 the early cretaceous m2 the late cretaceous a 34 the cretaceous tertiary boundary a 29 the early eocene a 24 the middle eocene a 21 and the oligocene miocene boundary a 6 observed and predicted paleomagnetic directions for landmasses other than africa i e for western europe china korea india and australia which are not used in deriving the maps are found to be in very satisfactory agreement the quality of reconstructions as measured by rms discrepancies between observed and predicted paleolatitudes is of the order of 5 the main features of the paleogeographic maps are discussed with the aim of relating oceanic and on land geodynamic events the constraints provided by tibetan data on the india asia collision and the 20 ma reorganization of geology and plate boundaries in a band extending from the south china sea through tibet to the red sea system are discussed in more detail we emphasize the episodic nature of the breakup of gondwanaland and formation of the indian ocean and subsequent interaction with eurasia in which six major discrete events are identified separated by 30 ma phases of rather uniform plate motion

water rich sublithospheric melt channel in the equatorial atlantic ocean the lithosphere asthenosphere boundary is the most extensive boundary on earth separating the mobile plate above from the convecting mantle below but its nature remains a matter of debate using an ultra deep seismic reflection technique here we show a systematic seismic image of two deep reflectors that we interpret as the upper and lower limits of the lithosphere asthenosphere boundary beneath a 40 70 million year old oceanic lithosphere in the atlantic ocean these two reflections correspond to 1 260 c and 1 355 c isotherms and bound a low velocity channel suggesting that the lithosphere asthenosphere boundary is thermally controlled we observe a clear age dependency of this sublithospheric channel its depth increases with age from 72 km where it is 40 myr old to 88 km where it is 70 myr old whereas its thickness decreases with age from 18 km to 12 km we suggest that partial melting facilitated by water is the main mechanism responsible for the low velocity channel the required water concentration for melting increases with age nevertheless its corresponding total mass remains relatively constant suggesting that most of the volatiles in the oceanic sublithospheric channel originate from a horizontal flux near the ridge axis

implications for melt transport and source heterogeneity in upwelling mantle from the magnitude of sp converted phases generated at the onset of melting recently detected converted phases in upwelling mantle have been attributed to seismic velocity gradients at the onset of melting in this study we investigate conditions required to generate a melting onset phase by combining melt migration models with synthetic receiver functions we find that increasing upwelling velocity increasing mantle viscosity and decreasing water content in the mantle source increases the predicted strength of a melting onset phase differences in these parameters between study regions may explain observed variations in converted phase magnitude for a wet mantle source the amplitude of synthetic receiver functions calculated using a standard relationship for dependence of seismic velocity on melt fraction is lower than observed amplitudes one possibility is that the observed receiver functions indicate a heterogeneous mantle source in which wet and dry components melt independently alternatively the dependence of seismic velocity on melt fraction at very low melt fractions may be much stronger than that used here as suggested in recent studies

current plate motions across the red sea summary a new objective method is used to estimate precisely 64 rates of seafloor spreading since chron 2a 3 2 ma from all the red sea magnetic profiles available from 15 9 to 26 n the fastest spreading rate 16 mm yr 1 occurs near 18 n whereas the slowest rate 10 mm yr 1 occurs at 25 5 n and is consistent with the rate predicted from arabia nubia data to the south the standard deviation of the spreading rates is 0 8 mm yr 1 much smaller than the median standard deviation of 4 mm yr 1 previously assigned to spreading rates in the global plate motion model nuvel 1 demets et al 1990 the fit to the spreading rates as well as the locations of earthquakes in and near the red sea indicate that spreading south of approximately 17 7 n is less than the full rate of spreading between the arabian and nubian plates the red sea spreading centre is instead usefully interpreted as the boundary between the arabian plate and a danakil microplate that includes the subaerial danakil block and a larger oceanic portion of lithosphere despite the absence of reliable azimuths of transform faults in the red sea all components of the angular velocity of arabia relative to nubia are usefully constrained from just the 45 relevant spreading rates the new compact 95 per cent confidence region of the angular velocity excludes prior estimates based on only four and two spreading rates by chase 1978 and jestin et al 1994 respectively 19 spreading rates in the southern red sea are used to estimate the angular velocity between the danakil microplate and arabia an approach based on singular value decomposition shows that without slip vectors only two of the three components of the angular velocity of the danakil microplate relative to arabia or nubia are usefully constrained but that all three components are usefully constrained if one earthquake slip vector is included

tomographic filtering of geodynamic models implications for model interpretation and large scale mantle structure 1 the resolution operator r is a critical accompaniment to tomographic models of the mantle r facilitates the comparison between conceptual three dimensional velocity models and tomographic models because it can filter these theoretical models to the spatial resolution of the tomographic model we compute r for the tomographic model s20rts ritsema et al 1999 2004 and two companion models that are based on the same data but derived with different norm damping values the three models explain within measurement uncertainty s sks and s skks travel times equally well to demonstrate how artifacts distort tomographic images and complicate model interpretation we apply r to 1 a thermochemical and 2 an isochemical model of convection in the mantle that feature different patterns of shear velocity heterogeneity in the deep mantle if we assume that shear velocity heterogeneity is caused by temperature variations only r suppresses short wavelength structures removes strong velocity gradients and introduces artificial stretching and tilting of velocity anomalies temperature anomalies in the thermochemical model resemble the spatial extent of low seismic velocity anomalies and the shear velocity spectrum in the d region better than the isochemical model however the thermochemical model overpredicts the amplitude of shear velocity variation and places the african and pacific anomalies imperfectly we suspect that inaccurate velocity scaling laws and uncertain initial conditions control these mismatches extensive hypothesis testing is required to identify successful models

analysis of tide and offshore storm induced water table fluctuations for structural characterization of a coastal island aquifer analysis of water table fluctuations can provide important insight into the hydraulic properties and structure of a coastal aquifer system including the connectivity between the aquifer and ocean this study presents an improved approach for characterizing a permeable heterogeneous coastal aquifer system through analysis of the propagation of the tidal signal as well as offshore storm pulse signals through a coastal aquifer offshore storms produce high wave activity but are not necessarily linked to significant onshore precipitation in this study we focused on offshore storm events during which no onshore precipitation occurred extensive groundwater level data collected on a sand barrier island sable island ns canada show nonuniform discontinuous propagation of the tide and offshore storm pulse signals through the aquifer with isolated inland areas showing enhanced response to both oceanic forcing signals propagation analysis suggests that isolated inland water table fluctuations may be caused by localized leakage from a confined aquifer that is connected to the ocean offshore but within the wave setup zone two dimensional groundwater flow simulations were conducted to test the leaky confined unconfined aquifer conceptualization and to identify the effect of key parameters on tidal signal propagation in leaky confined unconfined coastal aquifers this study illustrates that analysis of offshore storm signal propagation in addition to tidal signal propagation provides a valuable and low resource approach for large scale characterization of permeable heterogeneous coastal aquifers such an approach is needed for the effective management of coastal environments where water resources are threatened by human activities and the changing climate

cenozoic dynamics of the african plate with emphasis on the africa eurasia collision we investigate the forces acting on the african plate from about 56 ma to the present day we aim 1 to establish a basis for quantification of the overall dynamics of the mediterranean middle east region and 2 to understand how africa reacts to temporal variation in the forces associated with the africa eurasia collision we use the assumption of dynamical equilibrium which implies that at any given time the net torque of forces on africa s northern margin must have been balanced by the other torques acting on the plate moreover temporal variation in the forces on africa s northern margin must have been balanced by contemporaneous changes in one or more of the other torques by quantifying those other torques we may derive a constraint on the forces at the northern margin we expect the africa eurasia collision to have been associated with changes in the associated forces and examine how these changes are balanced using alternative parameterizations we calculate the net torque of ridge push transform fault resistance and basal drag we also quantify changes in the torque of forces associated with the continent ocean transition the results constitute a framework for detailed dynamic models of the africa eurasia convergence zone temporal changes in the ridge push torque prove to be small the torque of transform resistance does change but these changes do not relate to forces at the northern margin evidence for a relation with collision is found only for the case of basal drag the results point to a mechanism where the plate maintains its dynamical equilibrium mainly through collision induced changes in absolute motion

the role of inherited crustal and lithospheric architecture during the evolution of the red sea insights from three dimensional analogue experiments abstract the red sea represents an incipient ocean basin that is in transition from continental rifting to ocean initiation in this paper we assess how pre existing mantle and crustal heterogeneities within the african and arabian lithospheres influenced rift development and propagation using three dimensional 3d brittle ductile analogue laboratory experiments with a rotational boundary condition here we investigate the effect of three different configurations i pre existing rheological heterogeneities in the crust ii in the lithospheric mantle and iii in both the experimental outcomes show that a linear weak zone in the lithospheric mantle layer exerts a first order control on the overall mechanical behaviour of extending lithosphere rift architecture and style of rifting crustal heterogeneities control second order features such as the trends of the main boundary escarpments along the rift axis our experiments reproduce the geological history of the red sea and provide key insights on its deformation history in space and time we provide for the first time 3d dynamic constraints on how pre existing weak zones in the lithosphere have controlled the temporal evolution of the red sea and determined its distinctive present day morphology the results find that certain morphological elements of the red sea can only be explained by the combined effects of crustal and mantle lithosphere heterogeneities when applied to other natural examples our findings provide significant new insights on the role of crustal and lithospheric heterogeneities during continental rifting

modelling melting rates in upwelling mantle abstract upwelling regions of the mantle can undergo partial melting as a result of decompression many models for the dynamics of these regions have largely ignored the actual melting process or have prescribed a uniform melting rate proportional to the upwelling velocity this paper uses a simple model for an upwelling column to calculate the melting rate from conservation principles the model rock comprises two chemical components and is assumed to be in thermodynamic equilibrium for idealized linear phase constraints the melting rate can be calculated analytically and is found to be proportional to the average upwelling velocity of both the matrix and melt a secondary aim is to discuss reactive instabilities the model predicts that the one dimensional state will be linearly stable whereas previous models have suggested that reactive infiltration instability should occur this is argued to be a result of the background melting rate which has not usually been fully accounted for but which has a stabilizing effect the model here can also be applied to a column in which some melt is already present and in that case it does exhibit a channeling instability it is concluded that accounting for melt production consistently in mid ocean ridge models is important when assessing the likely modes of melt transport

stretching instabilities and lithospheric boudinage a competent layer with a nonlinear rheology can under extension exhibit pinch and swell instabilities such instabilities can explain small scale regular deformations of rock recently they have also been invoked in relation to the distribution of basins and ranges fletcher and hallet 1983 and the undulation of the bouguer map froidevaux 1986 in the western part of north america for the case of a simple stratified structure different authors have expressed mathematical solutions describing these instabilities in this paper we elucidate the general physical meaning of these nonhomogeneous deformations in the simplest case where the competent layer is fully embedded in a softer material a fundamental mode develops with a preferred wavelength equal to 4 times the thickness of a competent layer the novelty is the occurrence of overtones with alternate symmetrical and antisymmetrical deflections of the interfaces between layers these modes have a large amplitude when the rheology tends to perfect plasticity the presence of a free surface like the earth s surface for the upper crust is shown to double the preferred wavelength of the fundamental mode and to forbid one overtone out of two this new effect is caused by gravity which inhibits the growth of surface topography another shift is of opposite sign to that caused by a free surface numerical solutions for multilayer structures can elucidate the nature of couplings between various layers the case of a three layer lithosphere having a soft lower crust and overlying a mantle asthenosphere has been analyzed in some detail it shows that in the range of wavelengths larger than 30 km up to five preferred modes may be superimposed this should encourage further analysis of the western u s topography and gravity pattern which is known to contain at least two wavelengths about 40 km for the individual basins and ranges and about 200 km for the bouguer troughs in nevada and their associated broad topographies

the influence of olivine metastability on deep subduction of oceanic lithosphere abstract the dynamics of subduction of oceanic lithosphere is mainly controlled by thermal buoyancy forces and mineralogy in self consistent models the influence of metastable olivine mo on the dynamics and evolution of subducting oceanic lithosphere is investigated we consider both phase boundaries olivine ol to spinel sp 410 km and sp to perovskite pv magnesiowustite mw 660 km treating the kinetics of the phase transition in a simplified way while models do not allow the formation of mo penetrate straight into the lower mantle slabs with a sufficient amount of mo flatten out at the 660 km boundary before penetrating at a later stage furthermore our models show that the effect of mo to slow down subduction velocities within the transition zone is rather moderate

perturbation to the lithosphere along the hotspot track of la reunion from an offshore onshore seismic transect a 250 km long ne sw lithospheric transect spanning the 40 km wide island of la reunion and its submarine edifice is derived from lines of air gun shots at sea on either side along the assumed hotspot trace seismic records were obtained from an array spanning the whole transect and including sea bottom and land receivers providing a system of reversed and overlapping observations low seismic velocity and hence density is found on average for the whole edifice above the oceanic plate we attribute high velocity anomalies within the edifice to an intrusive core confined under the central northern quarter of the island crossing segment unexpectedly the main seismic interfaces top and bottom of the prevolcanic crust do not show significant flexural downwarping under the island in addition clear multipathing in the recorded wave field requires the presence of a body with a seismic velocity intermediate between the prevolcanic crustal material and the normal mantle this lithospheric structure provides the first example where underplating occur s beneath an active volcanic island suggesting a genetic relationship the underplated body could represent residues of the evolution of primary picritic melts that yield erupted basalts evidence for reflectors deeper in the lithosphere may indicate further related heterogeneity in the plate hotspot model commonly assumed the structural variation along the transect could be interpreted as a variation with time of the amount and physical state of underplated material

the relationship between plate velocity and trench viscosity in newtonian and power law subduction calculations convection with a newtonian temperature dependent rheology leads to little or no surface velocity unless zones of weakness are introduced plate like features are observed in calculations both with newtonian rheology employing imposed weak zones and with power law non newtonian rheology where high stresses at the trench reduce the effective viscosity since deformation at subduction zones involves faulting both of these parameterizations should be treated with some skepticism it is important to understand how the parameterizations affect the model results we study the relationship between trench viscosity and plate velocity using a newtonian rheology by varying the viscosity at the trench the plate velocity is a function of the trench viscosity for fixed rayleigh number and plate slab viscosity slab velocities for non newtonian rheology calculations are significantly different from slab velocities from newtonian rheology calculations at the same effective rayleigh number both models give reasonable strain rates for the slab when compared with estimates of seismic strain rate non newtonian rheology eliminates the need for imposed weak zones and provides a self consistent fluid dynamical mechanism for subduction in numerical convection models