the influence of post perovskite strength on the earth s mantle thermal and chemical evolution abstract we have investigated the influence of post perovskite ppv viscosity on mantle convective dynamics and stirring efficiency using numerical modeling and simple analytical theory our results show that the strength of ppv has a dramatic influence on convective dynamics the presence of a weak ppv enhances heat transfer across the bottom thermal boundary layer resulting in higher temperatures lower mantle viscosities and considerably larger convective velocities this leads to a significant increase in stirring efficiencies with decreasing the ppv strength by at least one order of magnitude in addition using a simple parameterized convection evolution that includes the influence of ppv coupled to a mixing model we show that during the long term history of the earth s mantle the presence of ppv yields systematically hotter thermal evolution and more efficient convective stirring such a strong effect of ppv strength on mantle stirring efficiency suggests that the influence of ppv phase must be considered when interpreting both geochemical and geophysical observations

aspherical structure of the core mantle boundary from pkp travel times maps of aspherical structure in the vicinity of the core mantle boundary have been constructed from a tomographic analysis of p df and p ab travel times in the distance range 150 180 these maps reveal long wavelength features not previously detected by mantle tomography the peak to peak amplitude of these features is on the order of 1 s in vertical travel time too large to be explained by conventional thermal boundary layer models or dynamically supported topography on the core mantle boundary we hypothesize the existence of one or more chemical boundary layers whose large scale accumulations on the surface of the core may in some ways be analogous to continents at the earth s surface

a simple model of plate generation from mantle flow summary a simple model of non newtoniancreeping ow is used to evaluate classes of rheologieswhichallow viscous mantle ow to becomeplate like the model describes shallow layerlithosphericmotion driven by sources and sinks the sources represent spreading ridges while the sinks rep resent subduction zones the sources and sinks thus also prescribe the poloidal component ofthe surface ow eld the toroidal strike slip component of the ow eld is found via thesolution of the stokes equation with non newtonian rheology as a rst basic investigation ofthe model the horizontal divergence from the two dimensional rectangular velocity eld ofolson bercovici 1991 is used for the source sink eld the degree to which the induced uid ow reproducesthe rectangularplate is used to measure the success of differentrheologiesin generating plate like ows results indicate that power law rheologies even in the limit ofvery high power law index can only produce modest plate like ow for example the ratioof toroidal to poloidal kinetic energy for a source sink eld derived from a square plate is atbest 0 65 whereas a perfect square plate has a ratio of 1 0 moreover the power law rheologyappears to reach an asymptotic limit in its ability to produce plate like behavior this impliesthat plate tectonics is unlikely to arise from a power law rheology even in the limit of very high a class of rheologies that yield signi cantly more promisi ng results arise from the carreaupseudo plastic rheology with the power law index taken to be 0 one rheology in thisclass is the continuum model for stick slip earthquake behavior of whitehead gans 1974 which is essentially the carreau equation with 1 this class of rheologies referred to asthe stick slip rheologies induces a toroidal to poloidal kinetic energy ratio for a square plates ssource sink function which can be as high as0 9 the viscosity or strength distribution forthis class of rheologies also appears more plate like showing fairly uniform high viscosity re gions pseudo plates and sharply de ned low viscosity zones pseudo margins in contrast even the most nonlinear power law rheology produces spatially varying high viscosity regionsand relatively smooth low viscosity margins the greater success of the stick slip rheologiesin producing plates is attributed to a self lubricating mechanism in which the transfer of mo mentum from regions of high shear to low shear is inhibited in contrast even in the limit ofin nite power law index a power law rheology can retard but never prohibit momentum trans fer this feature is essential to the sharpening of velocity pro les into plate like pro les whichis illustrated with a simple boundary layertheory key words plate tectonics mantle convection non newtonian ow toroidal poloidal cou pling

glacial isostatic adjustment and mantle viscosity the determination of mantle viscosity conditions from observations of relative sea level changes associated with glacial unloading of the laurentide and fennoscandia ice sheets has been reexamined the conclusion is reached that the adjustment has been an upper mantle process associated with the low seismic velocity low viscosity asthenospheric layer rather than a whole mantle process of essentially uniform mantle viscosity the critical data for resolving this dichotomy are those from the zone of peripheral subsidence adjacent to the region of central uplift of the laurentide and fennoscandia ice sheets

inferring the viscosity and the 3 d density structure of the mantle from geoid topography and plate velocities summary in a dynamic earth mantle mass heterogeneities induce gravity anomalies surface velocities and surface topography these lateral density heterogeneities can be estimated on the basis of seismic tomographic models recent papers have described a realistic circulation model that takes into account the observed plate geometry and is able to predict the rotation vectors of the present plates the relationship between the surface observables and the heterogeneities is sensitive to the viscosity stratification of the mantle here we use this model combined with a generalized least squares method in order to infer the viscosity profile of the earth from the surface observations and to get some new insight into the 3 d density structure of the mantle the computed radial viscosity profile presents a continuous increase of more than two orders of magnitude the asthenosphere has a viscosity close to 2 x lo2 pa s no sharp discontinuity is requested at the upper lower mantle interface the largest viscosity 7 x 10 pas is reached in the middle of the lower mantle at greater depth approaching the core mantle boundary the viscosity decreases by one order of magnitude the model suggests that the well known degree 2 and order 2 anomaly in the transition zone of the upper mantle is merely the signature of the slabs it also slightly increases the degree 2 and order 0 in the lower mantle and decreases it in the upper mantle in other words the inversion requests a hotter lower mantle beneath the equator and a colder upper mantle at the same latitudes

influence of lithospheric and mantle stratification on global post seismic deformation a new class of multilayered viscoelastic earth models based on prem with an incompressible linear viscoelastic maxwell rheology is applied to the modeling of co seismic and post seismic deformation following large earthquakes the novelty of our approach stands on the usage for the first time in post seismic studies of a fully analytical scheme based on normal mode theory that allows to deal with some of the complexities of the real earth such as lithospheric and mantle layering sphericity and self gravitation for a dip slip source used as a test case elastic layering of the lithosphere is necessary for the modeling of realistic patterns of post seismic deformation although only a few layers namely three are sufficient to reproduce correct deformation patterns in consequence of the fact that below the crust the lithosphere inferred from prem is rather homogeneous the viscosity structure of the mantle has a major influence on post seismic deformation in the far field at distances from the epicenter of the order of 100 km at least for lithospheric sources precise modeling of post seismic deformation does not require for the mantle the same fineness of layering that is necessary for the lithosphere our findings are key for future analyses of the co seismic and post seismic effects of large earthquakes in view of a correct interpretation of vlbi gps and satellite differential radar altimetry data

thermal expansion measurements at very high pressure systematics and a case for a chemically homogeneous mantle new data on thermal expansivity of iron to over 400 kbar in laser and electrically heated diamond anvil cells and on mgo and forsterite from spectroscopic data to over 200 kbar and previously measured adiabats boehler 1982 imply that at high pressure may be described by a constant value of 1n 1nv t 5 5 0 5 this behavior is supported by previous measurements at lower pressures on a large variety of materials densities of magnesio wustite and silicate perovskite in a pyrolitic ratio calculated at lower mantle conditions using a finite strain equation and at high pressures using the above systematics match those in prem through the entire lower mantle making a strong case for a uniform mantle

glacial isostatic adjustment i the forward problem summary the isostatic adjustment of a radially stratified visco elastic spheroid is treated using space time green functions for the associated surface mass load boundary value problem these impulse response functions are convolved with a heaviside function to give the time dependent deformation of the planet which would be produced by a unit point mass brought up from infinity at t 0 and allowed to remain on the surface the resulting heaviside green functions can be employed to simulate all of the important signatures of glacial isostatic adjustment given a space time dependent surface mass load consisting of ice sheet ablation histories and a model of the simultaneous filling of the ocean basins these source terms are simply convolved with the heaviside green function appropriate to a specific response signature a realistic model of the spatial distribution of the main late pleistocene ice loads and of their temporal disintegration is constructed the response of two visco elastic earth models is computed and compared to a global set of relaxation data relative sea level curves on the basis of this initial comparison of theory and observation the possibility that the lower mantle has a viscosity which is significantly in excess of the viscosity of the upper mantle is excluded in addition clear evidence of the presence of the lithosphere has been found in relaxation data from sites which were near the edge of the laurentide ice sheet such data should therefore prove useful as a basis for analysis of lateral variations in lithosphere thickness further extensions of the calculation are suggested

effects of lateral viscosity variations on long wavelength geoid anomalies and topography models for the earth s geoid and dynamically supported topography usually account only for radial vertical variations in mantle rheology and flow structure however expected lateral variations in viscosity due to temperature or stress dependent rheology can degrade our ability to infer mantle structure from geoid models based on seismically determined density contrasts in the mantel theoretically the main difficulty is that the harmonic components of the flow stress deformation field are no longer independent and exact analytical methods cannot account for this mode coupling for large thermal density contrasts of small lateral dimension e g subducted slabs it is difficult to form any general theory for these effects however we believe that the largest geoid anomalies which occur at the lowest harmonic degrees l 2 3 are due to a global pattern of relatively mild temperature variations for convection with mainly broadscale thermal and rheological changes some relatively straightforward ideas can be developed combining results from simple perturbation theory and numerical models for convective flow we find the following 1 if we are correct in inferring viscosity variations of less than one order of magnitude globally the longest wavelength geoid anomalies l 2 3 are not seriously contaminated by lateral viscosity variations 2 contamination of the higher degrees l 4 is expected and in particular self coupling between l 2 heterogeneity and viscosity variations may be observable at degree 4 the most severe contamination occurs at the doubled harmonic l 4 of the dominant heterogeneity harmonic l 2 3 perturbation techniques are adequate to model global scale lateral viscosity variations of an order of magnitude or less 4 the steady state surface deformation problem at the heart of our geoid models may be more robust with respect to broadscale lateral viscosity variations than are models for transient deformations e g postglacial rebound 5 for loading by broadscale density heterogeneity the effects of stress dependent rheology are probably secondary to the effects of radial viscosity stratification

pleistocene deglaciation and the earth s rotation implications for mantle viscosity summary recent results from the analysis of postglacial rebound data suggest that the viscosity of the earth s mantle increases through the transition region models which fit both relative sea level and free air gravity data have viscosities which increase from a value near 1022 poise in the upper mantle beneath the lithosphere to a value of about 1023 poise in the lower mantle in this paper we analyse the effect of deglaciation upon the earth s rotation and thereby show that the observed secular trend polar wander evident in the ils ipms pole path and measurements of the non tidal acceleration of the length of day are both consistent with the viscosity profile deduced from postglacial rebound the two analyses are therefore mutually reinforcing

the impact of variability in the rheological activation parameters on lower mantle viscosity stratification and its dynamics the recently discovered spin state crossover in iron in major mantle minerals at high pressures should exert dramatic influences on transport properties such as the activation creep parameters in the deep mantle wentzcovitch et al 2009 have computed the elastic parameter of ferropericlase that is affected by this crossover the bulk modulus at high realistic conditions using first principles density functional theory plus hubbard u as a consequence of the spin crossover there is a significant softening of the bulk modulus and an attendant reduction of the activation energy in the creep law using a parameterized model capturing the basic physics we have studied the dynamical consequences within the framework of a 2 d cartesian convection model our models reveal a series of asthenospheres and are characterized first by a low viscosity channel below the 670 km boundary second asthenosphere caused by an increase of temperature and decrease in grain size reduction superplascity with the post spinel transition which results in partially layered convection further the variability of the rheological parameters leads to the prevalent formation of viscosity minimum at a depth of about 1600 km third asthenosphere caused by spin crossovers followed by a viscosity hill in the bottom half of the lower mantle the latter due to the increase of the activation energy in some low spin irons at the bottom of the mantle our numerical simulations reveal a tendency to the formation of small scale convection below the 670 km boundary and bigger stabilized plumes and or plume clustering in the deep lower mantle other material properties at the bottom of the mantle decrease of thermal expansivity radiative thermal conductivity also exert significant influence on the multi scale lower mantle plume dynamics

observational hints for a plume fed suboceanic asthenosphere and its role in mantle convection an asthenosphere layer which is entirely fed from below by plumes and which loses equal mass by accretion to the overlying oceanic lithosphere and at subduction zones may play a critical role in shaping the form of mantle convection in this study we discuss geochemical seismic and geoid depth evidence for lateral flow within this type of asthenosphere in particular we suggest that there are large scale layered horizontal flow structures that connect upward plume input beneath hotspots to near ridge regions of increased asthenosphere accretion into the growing oceanic lithosphere lateral asthenosphere flow is also shaped by oceanic subduction zones with a partial return flow from trenches and by deep continental roots that are migrating barriers to asthenosphere flow this alternative paradigm offers relatively simple explanations for several puzzles about mantle convection for example the low mantle heat flow beneath continents it also offers an explanation for why mid ocean ridges appear to be passive features that migrate with little geochemical or morphological change with respect to the lower mantle and seem to be uncoupled from large scale mantle flow while in contrast trenches appear to be strongly coupled to mantle thick regions of fast colder seismic velocity anomalies we also discuss several implications of this paradigm that should be testable in future studies such as the prediction of cogenetic off axis seamount volcanism that is created between an off axis hotspot and its neighboring ridge axis

viscosity undulations in the lower mantle the dynamical role of iron spin transition a b s t r a c t a proper determination of the lower mantle viscosity profile is fundamental to understanding earth geodynamics based on results coming from different sources several models have been proposed to constrain the variations of viscosity as a function of pressure stress and temperature while some models have proposed a relatively modest viscosity variation across the lower mantle others have proposed variations of several orders of magnitude here we have determined the viscosity of ferropericlase a major mantle mineral and explored the role of the iron high to low spin transition viscosity was described within the elastic strain energy model in which the activation parameters are obtained from the bulk and shear wave velocities those velocities were computed combining first principles total energy calculations and the quasi harmonic approximation as ar esult of a strong elasticity softening across the spin transition there is a large reduction in the activation free energies of the materials creep properties leading to viscosity undulations these results suggest that the variations of the viscosity across the lower mantle resulting from geoid inversion and postglacial rebound studies may be caused by the iron spin transition in mantle minerals implications of the undulated lower mantle viscosity profile exist for both down and up wellings in the mantle we find that a viscosity profile characterized by an activation free energy of g z0 300 400 kj mol based on diffusion creep and dilation factor 0 5b etter fits the observed high velocity layer at mid mantle depths which can be explained by the stagnation and mixing of mantle material our model also accounts for the growth of mantle plume heads up to the size necessary to explain the large igneous provinces that characterize the start of most plume tracks

the viscosity of earth s lower mantle inferred from sinking speed of subducted lithosphere abstract the viscosity of the mantle is indispensable for predicting earth s mechanical behavior at scales ranging from deep mantle material flow to local stress accumulation in earthquakes zones but mantle viscosity is not well determined for the lower mantle particularly only few constraints result from elaborate high pressure experiments karato 2008 and a variety of viscosity depth profiles result from joint inversion of the geoid and postglacial rebound data forte and mitrovica 1996 kaufmann and lambeck 2000 mitrovica and forte 2004 here we use inferred lower mantle sinking speed of lithosphere subduction remnants as a unique internal constraint on modeling the viscosity profile this entails a series of elaborate dynamic subduction calculations spanning a range of viscosity profiles from which we select profiles that predict the inferred sinking speed of 12 3 mm yr van der meer et al 2010 our modeling shows that sinking speed is very sensitive to lower mantle viscosity good predictions of sinking speed are obtained for nearly constant lower mantle viscosity of about 3 4 1022 pa s viscosity profiles incorporating a viscosity maximum in the deep lower mantle as proposed in numerous studies only lead to a good prediction in combination with a weak postperovskite layer at the bottom of the lower mantle and only for a depth average viscosity of 5 1022 pa s

time dependent convection models of mantle thermal structure constrained by seismic tomography and geodynamics implications for mantle plume dynamics and cmb heat flux summary one of the outstanding problems in modern geodynamics is the development of thermal convection models that are consistent with the present day flow dynamics in the earth s mantle in accord with seismic tomographic images of 3 d earth structure and that are also capable of providing a time dependent evolution of the mantle thermal structure that is as realistic earth like as possible a successful realization of this objective would provide a realistic model of 3 d mantle convection that has optimal consistency with a wide suite of seismic geodynamic and mineral physical constraints on mantle structure and thermodynamic properties to address this challenge we have constructed a time dependent compressible convection model in 3 d spherical geometry that is consistent with tomography based instantaneous flow dynamics using an updated and revised pseudo spectral numerical method the novel feature of our numerical solutions is that the equations of conservation of mass and momentum are solved only once in terms of spectral green s functions we initially focus on the theory and numerical methods employed to solve the equation of thermal energy conservation using the green s function solutions for the equation of motion with special attention placed on the numerical accuracy and stability of the convection solutions a particular concern is the verification of the global energy balance in the dissipative compressible mantle formulation we adopt such validation is essential because we then present geodynamically constrained convection solutions over billion year timescales starting from present day seismically constrained thermal images of the mantle the use of geodynamically constrained spectral green s functions facilitates the modelling of the dynamic impact on the mantle evolution of 1 depth dependent thermal conductivity profiles 2 extreme variations of viscosity over depth and 3 different surface boundary conditions in this case mobile surface plates and a rigid surface the thermal interpretation of seismic tomography models does not provide a radial profile of the horizontally averaged temperature i e the geotherm in the mantle one important goal of this study is to obtain a steady state geotherm with boundary layers which satisfies energy balance of the system and provides the starting point for more realistic numerical simulations of the earth s evolution we obtain surface heat flux in the range of earth like values 37 tw for a rigid surface and 44 tw for a surface with tectonic plates coupled to the mantle flow also our convection simulations deliver cmb heat flux that is on the high end of previously estimated values namely 13 tw and 20 tw for rigid and plate like surface boundary conditions respectively we finally employ these two end member surface boundary conditions to explore the very long time scale evolution of convection over billion year time windows these billion year scale simulations will allow us to determine the extent to which a memory of the starting tomography based thermal structure is preserved and hence to explore the longevity of the structures in the present day mantle the two surface boundary conditions along with the geodynamically inferred radial viscosity profiles yield steady state convective flows that are dominated by long wavelengths throughout the lower mantle the rigid surface condition yields a spectrum of mantle heterogeneity dominated by spherical harmonic degree 3 and 4 and the plate like surface condition yields a pattern dominated by degree 1 our exploration of the time dependence of the spatial heterogeneity shows that for both types of surface boundary condition deep mantle hot upwellings resolved in the present day tomography model are durable and stable features these deeply rooted mantle plumes show remarkable longevity over very long geological time spans mainly owing to the geodynamically inferred high viscosity in the lower mantle

new perspectives on mantle dynamics from high resolution seismic tomographic model p1200 recently a high resolution tomographic model the p1200 based on p wave travel times was developed which allowed for detailed imaging of the top 1200 km of the mantle this model was used in diverse ways to study mantle viscosity structure and geodynamical processes in the spatial domain there are lateral variations in the transition zone suggesting interaction between the lower mantle plumes and the region from 600 km to 1000 km some examples shown here include the continental region underneath manchuria ukraine and south africa where horizontal structures lie above or below the 660 km discontinuity the blockage of upwelling is observed under central africa and the interaction between the upwelling and the transition zone under the slow icelandic region appears to be complex an expansion of the aspherical seismic velocities has been taken out to spherical harmonics of degree 60 for degrees exceeding around 10 the spectra at various depths decay with a power law like dependence on the degree with the logarithmic slopes in the asymptotic portion of the spectra containing values between 2 and 2 6 these spectral results may suggest the time dependent nature of mantle convection details of the viscosity structure in the top 1200 km of the mantle have been inferred both from global and regional geoid data and from the high resolution tomographic model we have considered only the intermediate degrees l 12 25 in the nonlinear inversion with a genetic algorithm approach several families of acceptable viscosity profiles are found for both oceanic and global data the families of solutions for the two data sets have different characteristics most of the solutions asociated with the global geoid data show the presence of asthenosphere below the lithosphere in other families a low viscosity zone between 400 and 600 km depth is found to lie atop a viscosity jump other families evidence a viscosity decrease across the 660 km discontinuity solutions from oceanic geoid show basically two low viscosity zones one lying right below the lithosphere the other right under 660 km depth all of these results bespeak clearly the plausible existence of strong vertical viscosity stratification in the top 1000 km of the mantle the presence of the second asthenosphere may have important dynamical ramifications on issues pertaining to layered mantle convection numerical modelling of mantle convection with two phase transitions and a realistic temperature and pressure dependent viscosity demonstrates that a low viscosity region under the endothermic phase transition can indeed be generated self consistently in time dependent situations involving a partially layered configuration in an axisymmetric spherical shell model

the driving mechanism of plate tectonics abstract plate velocities reflect the global mantle circulation which carries the heat from the deep part of the earth to the surface three other measurable parameters put strong constraints on the physical processes involved in mantle dynamics the surface topography the geoid and the lateral density heterogeneities mapped throughout the mantle by seismic tomography using these internal loads and the navier stokes equations a first generation of dynamic earth models have been developed in the last five years they have produced satisfactory predictions for the divergence of the surface velocities and for the geoid these predictions are quite sensitive to the assumed radial viscosity variations the best match with observations has therefore helped to discriminate between various possible viscosity structures a new set of models taking the existence of plates into account has been particularly successful in predicting the plate velocities including their vorticity the driving mechanism of plate tectonics is thus put on a firm physical basis and the computed solutions are quantitatively reasonable

mantle dynamics constrained by degree 6 surface topography seismic tomography and geoid inference on the origin of the south pacific superswell abstract the relationship between global geophysical observables i e seismic tomography geoid and earth surface topography is investigated this study focuses on the high correlation at degree 6 which in addition to degree 2 characterizes these global fields surface topography is found to be highly correlated at degree 6 to seismic velocity anomalies between the 150 km and 450 km depths in the upper mantle several highs of the degree 6 topography coincide with particular provinces of the earth surface among others the south pacific superswell and the afar region these provinces overlie low velocity hot upper mantle by developing a global circulation model in a compressible radially viscosity stratified mantle we show that the degree 6 global fields can be explained as the result of the mantle convective flow driven by density anomalies the viscosity structure that best fits the degree 6 observables also explains the low degree 2 and 3 components of the dynamical topography and geoid the best fit model has a viscosity increase of factor 30 at the upper mantle lower mantle discontinuity and a low viscosity layer extending from the 100 to the 400 km depth in the upper mantle with a viscosity drop of factor 30 with respect to the transition zone these values are in good agreement with most previously published results based on the low degree geoid alone from this study we infer that the south pacific superswell which lies over a broad zone of hot upper mantle extending down to the 450 km depth is generated by the degree 6 flow of the mantle circulation

mapping the upper mantle three dimensional modeling of earth structure by inversion of seismic waveforms a method is presented for the inversion of waveform data for the three dimensional distribution of seismic wave velocities the method is applied to data from the global digital networks international deployment of accelerometers global digital seismograph network the selected data set consists of some 2000 seismograms corresponding to 53 events and 870 paths the moment tensors of the earthquakes are determined through an iterative procedure which minimizes the corrupting influence of lateral heterogeneity a global model is constructed for shear wave velocity expanded up to degree and order 8 in spherical harmonics and described by a cubic polynomial in depth for the upper 670 km of the earth s mantle although no a priori information is incorporated the model predictions reproduce much of what is known about the dispersion of mantle waves for example high phase velocities for shields low velocities at ridges and a strong degree 2 pattern for rayleigh waves since the method makes use of complete waveforms overtone data are also included it is shown that the model is reproducible in that substantially the same model can be constructed from each half of the total data set considered independently the model shows that shields and ridges are major features in the depth interval 25 250 km the ridges of the southern pacific and the larger shields persist to 350 km but the southeast indian rise is underlain by a high velocity anomaly at this depth as is much of the mid atlantic ridge at 450 650 km the major features are a broad region of high velocities incorporating south america much of the south atlantic and parts of west africa a broad region of low velocities in the central and eastern pacific high velocities in the western pacific and a low velocity anomaly beneath the red sea and the gulf of aden in the absence of a crustal correction degrees 2 and 3 show a high positive correlation with the geoid paradoxically this is largely destroyed when the distribution in crustal thickness is taken into account spherical harmonic degrees 4 7 show a significant negative correlation

inferences of viscosity from the oceanic geoid indication of a low viscosity zone below the 660 km discontinuity abstract we have attempted to infer details of the viscosity structure in the top 1000 km of the mantle from the geoid and tomographic structure beneath the oceans in order to eliminate the gravity signal from problematic masses located below the subduction zones and the continents we have considered only the intermediate degrees of the oceanic geoid l 12 25 a genetic algorithm has been used to determine the family of viscosity models which give the best correlation with the observed geoid our inversion clearly identifies the asthenosphere just below the lithosphere and also confirms the viscosity increase in the lower mantle predicted by previous inferences but suggests that the main viscosity jump occurs at a depth of about 1000 km and not at the usually stated 660 km boundary somewhere in the depth range of 400 1000 km a low viscosity zone may exist where the viscosity decreases to a value comparable with the asthenosphere existence of such a low viscosity zone is supported by recent analysis of deep mantle anisotropy which favours a flow pattern with a strong horizontal component in the top part of the lower mantle unfortunately the resolution of the inversion as well as the quality of recent seismic tomographic models are not sufficient to localize the depth and to come up with a higher accuracy for the viscosity of this low viscosity channel

glacial isostatic adjustment and the free air gravity anomaly as a constraint on deep mantle viscosity summary a gravitationally self consistent theory for relative sea level variations forced by pleistocene deglaciation events is employed to explore the extent to which rsl and free air gravity observations together constrain the viscosity of the mantle beneath the seismic discontinuity at 670 km depth a trade off is revealed between errors in the assumed deglaciation history and errors in the inferred value of the viscosity of the lower mantle taking full account of such uncertainty plausible values of the viscosity beneath the transition region are bounded above by 1023 poise cgs units the preferred value is at least a factor of 2 lower than this and is strongly constrained by the observed free air gravity anomaly over hudson bay the calculations described in detail here show for the first time that the relatively large gravity anomalies observed over sites of wurm wisconsin deglaciation do not require any pronounced increase of mantle viscosity with depth in order to explain them

various influences on plumes and dynamics in time dependent compressible mantle convection in 3 d spherical shell the effects on the dynamics of compressible mantle convection from increasing rayleigh number viscosity stratification and internal heating are investigated in a time dependent 3 d spherical shell model resolution of up to spherical harmonics of degree 192 and 150 higher order finite difference grid points unevenly distributed in the radial direction have been employed surface rayleigh number varies from 5 104 to o 107 the viscosity profiles used range from constant to a depth dependent viscosity profile the adams williamson adiabatic assumption is adopted for describing the equation of state both depth dependences of thermal expansivity and thermal conductivity are included we have focused on the viscosity stratification in the lower mantle the mean temperature profiles are very sensitive to the input parameters such as the amount of internal heating and the dissipation number in the isoviscous model the number of hot plumes decreases for low rayleigh number 5 104 and increases dramatically with higher rayleigh number 106 in the model with a high viscosity zone in the lower mantle the number of plumes is small and does not change too much when the rayleigh number varies from 3 106 to 8 106 in the purely base heated viscously stratified models well organized triple junctions of the downwelling sheets are developed near the surface and interconnected cold walls of downwellings are produced around upwelling plumes the presence of internal heating increases the number of upwelling plumes and breaks up the interconnection of downwelling sheets strong viscous heating in the compressible solutions is found near the top and bottom at the higher values of rayleigh number the laterally averaged value of viscous dissipation is comparable with the chondritic strength of internal heating whether the maximum viscous heating is found in upwelling plumes or downwelling flows depends on the presence of internal heating and the mechanical boundary conditions the effects of compressibility on local structures and plume morphologies can be significant however these compressible effects are moderate on the global patterns

effects of multiple phase transitions in a three dimensional spherical model of convection in earth s mantle numerical models of mantle convection that incorporate the major mantle phase changes of the transition zone reveal an inherently three dimensional flow pattern with cylindrical features and linear features that behave differently in their ability to penetrate the 670 km discontinuity the dynamics are dominated by accumulation of cold linear downwellings into rounded pools above the endothermic phase change at 670 km depth resulting in frequent avalanches of upper mantle material into the lower mantle the effect of the exothermic phase transition at 400 km depth is to reduce the overall degree of layering by pushing material through the 670 km phase change resulting in smaller and more frequent avalanches and a wider range of morphologies large quantities of avalanched cold material accumulate above the coremantle boundary cmb resulting in a region of strongly depressed mean temperature at the base of the mantle the 670 km phase change has a strong effect on the temperature field with three distinct regions being visible 1 the upper mantle containing linear downwellings and pools of cold material in the transition zone and characterized by a high amplitude long wavelength spectrum 2 the midmantle containing quasi cylindrical avalanche conduits and characterized by a low amplitude broad spectrum and 3 the deep mantle containing large pools of cold avalanched material and characterized by a high amplitude ultra red i e long wavelength spectrum the effect on the velocity field is very different flow penetration across the 670 km phase change is strongly wavelength dependent with easy penetration at long wavelengths but strong inhibition at short wavelengths thus when comparing numerical models with long wavelength seismic tomography diagnostics based on the density field such as the radial correlation function are much more sensitive to the effects of phase transitions than those based on the velocity field the amplitude of the geoid is not significantly affected by the partial layering because the contribution from the strong heterogeneity in the transition zone is almost perfectly balanced by the contribution from deflection of the 670 km discontinuity avalanches are associated with geoid lows however a more complex viscosity structure is required to correctly match the sign of the geoid over slabs in earth

melting temperature of the earth s mantle and core earth s thermal structure abstract although the temperature at the top of the lower mantle is well constrained by phase equilibrium data for the transformation of transition zone minerals to the denser perovskite polymorphs the temperature distribution in the lower mantle is poorly known models depend strongly on the assumptions of the amount of internal heating and the viscosity profile new melting data on iron to pressures of the outer core 2 mbar and the observed strong decrease of eutectic melting depression in the fe feo fes system with increasing pressure however tightly constrain the temperature at the inner core boundary to slightly less than 5000 k this estimate can be reconciled with all recent static melting measurements on iron and lays within the uncertainty of shock temperature measurements the resulting temperature at the core mantle boundary of about 4000 k then requires a large temperature gradient at the bottom of the lower mantle of about 1500 k recent findings of the very high melting temperatures o

a simple global model of plate dynamics and mantle convection cooling and thickening of lithospheric plates with age and subduction result in large scale horizontal density contrasts tending to drive plate motions and mantle flow we quantify the driving forces associated with these density contrasts to determine if they can drive the observed plate motions first two dimensional models are computed to evaluate the effects of assumed rheologies and boundary conditions we are unable to obtain platelike behavior in viscous models with traction free boundary conditions the piecewise uniform velocities distinctive of plate motion can be imposed as boundary conditions and the dynamic consistency of the models evaluated by determining if the net force on each vanishes if the lithosphere has a newtonian viscous rheology the net force on any plate is a strong function of the effective grid spacing used leading to ambiguities in interpretation incorporating a rigid plastic lithosphere which fails at a critical yield stress into the otherwise viscous model removes these ambiguities the model is extended to the actual three dimensional spherical plate geometry the observed velocities of rigid plastic plates are matched to the solution of the viscous stokes equation at the lithosphere asthenosphere boundary body forces from the seismically observed slabs from the thickening of the lithosphere obtained from the actual lithospheric ages and from the differences in structure between continents and oceans are included interior density contrasts such as those resulting from upwellings from a hot bottom boundary layer are assumed to occur on a scale small compared to plate dimensions and are not included the driving forces from the density contrasts within the plates are calculated and compared to resisting forces resulting from viscous drag computed from the three dimensional global return flow and resistance to deformation at converging boundaries the rms residual torque is 30 of the driving torque the density contrasts within the plates themselves can reasonably account for plate motions body forces from convection in the interior may provide only a small net force on the plates at converging boundaries the lithosphere has a yield stress of 100 bars drag at the base of the plates is 5 bars and resists plate motion the net driving forces from subducting slabs and collisional resistance are localized and approximately balance driving forces from lithospheric thickening are distributed over the areas of the plates as is viscous drag the approximate balance of these two forces predicts plate velocities uncorrelated with plate area as observed the model represents a specific case of boundary layer convection the dynamical results are consistent with either upper mantle or mantle wide convection

topography of the ocean floor thermal evolution of the lithosphere and interaction of deep mantle heterogeneities with the lithosphere topography associated with large wave length geoid is derived by two methods first the bathymetry corrected for sediment loading is decomposed using a least squares method as a function of age plus terms proportional to the geoid the topography at old ages is flatter predicted by the law at degrees 2 and 3 very little topography is associated with the geoid the transfer function is not even certain in sign but is never smaller than 0 6 more topography seems associated with higher degree harmonics 6 and 7 second the function of age alone best fitting the topography is determined the correlation coefficients of the residual topography with the geoid are then computed a slight correlation appears for degree 6 but not for degrees 2 and 3 these results suggest a mechanical decoupling between the lower mantle mass heterogeneities and the lithosphere

internal loading of a compressible earth effect of a solid lithosphere abstract this paper shows the importance of including a lithosphere in the earth model when inferring a mantle viscosity profile or density velocity conversion factor from mantle circulation the free slip and no slip conditions are investigated at the mantle lithosphere boundary results obtained for earth models that include a lithosphere are compared with those obtained for models where no lithosphere is considered in particular it is shown that models without a lithosphere are inconsistent with the no slip condition at the mantle lithosphere interface because they introduce a tangential stress discontinuity at that level with regard to the estimations of the mantle viscosity profile from the comparison between observed and computed geoids and plate velocities the principal results obtained previously using the no slip condition without a lithosphere i e the best fit is obtained for a one cell convection with a large viscosity jump at a depth of 670 km are still valid nevertheless to infer a precise viscosity profile or density velocity conversion factor when using the no slip condition a lithosphere must be accounted for because the results are different for models with and without a lithosphere

geodynamical implications from the correlation of surface geology and seismic tomographic structure abstract we have investigated the relationship between the lateral seismic structure in the mantle and the geological features associated with major mantle upwellings the oceanic ridge system and hotspots using the technique of correlation analysis we have identified that these features extend quite deeply into the mantle to a depth of around 1300 km this finding is surprising for the ridges because it indicates that the style of upwelling is rather active and is somewhat two dimensional in character the confidence level of the correlation between ridges and slow seismic regions is maintained at a level of greater than 90 to a depth of 1300 km this is followed by a sharp decrease to 50 at 1700 km a similar correlation pattern is found for the hotspots in both cases the correlation maxima lie close to a depth of 1000 km where remnants of old slabs have recently been revealed by correlation analysis these findings suggest that dynamical processes at a depth of around 1000 km can be linked with surface tectonic activity this correlation disappears at a depth of around 1500 1700 km certain features of the surface tectonic pattern can again be recognized around a depth of 2000 km and in the deepest mantle where major hotspot plumes and the remnant old lithosphere are the dominant features there may be some dynamical relationship between the ridge system and the two megaplumes we have proposed a physical mechanism based on viscous heating within megaplumes to explain the two peaks in the correlation between the hotspots and the seismic heterogeneities

insights from scs s measurements on deep mantle attenuation abstract we apply a recently developed method based on the instantaneous frequency to analyze broadband seismic data recorded by the transportable usarray we measure in the frequency band 0 018 0 2 hz about 700 high quality differential scs s anelastic delay times t scs s sampling the mantle below central america and below alaska that we compare to elastic delay times t scs s obtained by cross correlating the s and scs signals we confirm that the instantaneous frequency matching method is more robust than the classical spectral ratio method by a series of careful analyses of the effects of signal to noise ratio source mechanism characteristics and possible phase interferences on measurements of differential anelastic delay times we demonstrate that in order to obtain accurate values of t scs s the seismic records must be rigorously selected in spite of the limited number of data that satisfy our quality criteria we recover using an additional stacking procedure a clear dependence of t scs s on the epicentral distance in the two regions the absence of correlation between the obtained anelastic and elastic delay times indicates a complex compositional thermal origin of the attenuation structure or effects of scattering by small scale structure in accordance with possible presence of subducted material the regional 1 d inversions of our measurements indicate a non uniform lower mantle attenuation structure a zone with high attenuation in the mid lower mantle q 250 and a low attenuation layer at its base q 450 a comparison of our results with low frequency normal model q models is consistent with frequency dependent attenuation with q and 0 1 0 2 i e less attenuation at higher frequencies although possible effects of lateral variations in q in the deep mantle add some uncertainty to these values

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

ice age cycles earth s rotation instabilities and sea level changes a new class of multilayered viscoelastic earth models based on prem is applied to the modeling of earth s rotation instabilities and associated sea level changes induced by the occurrence of pleistocene ice age cycles that match the oxygen isotope records over the last 0 8 myr the novelty of our approach stands on the usage for the first time in post glacial rebound induced sea level studies of a fully analytical scheme based on normal mode theory that allows to deal with complexities of the real earth such as lithospheric and mantle layering sphericity and self gravitation the ice models are based on ice 3g our results show that differences in true polar wander tpw and tpw induced sea level changes between earth models with a few layers and models having enough layers so that saturated continuum limits are reached can amount factors 2 to 3 this may change conclusions derived from some earlier studies on polar wander induced changes in climate that were based on earth models with a limited amount of layers the results indicate that models containing about 15 layers with a stratified lower mantle and transition zone have reached saturated continuum responses for both tpw and tpw induced sea level changes stratification of the lithosphere is not important this is in contrast with the sensitivity of post seismic deformation models this difference in sensitivity on the radial profile of the earth model can be explained with simple arguments on the dependence on zonal degree

late pleistocene and holocene sea level change in the australian region and mantle rheology summary spatial and temporal variations in sea level are produced by the melting of the late pleistocene ice and by the earth s response to the redistribution in surface loads by examining different parts of the sea level curves of the past 20 000 yr from geographically widely distributed regions it becomes possible to constrain models of the melting history of the ice sheets and of the earth s rheology observations from sites away from the former arctic ice sheets such as the australian and south pacific region are particularly important for constraining the total meltwater volumes added into the oceans in the past 20 000 yr and the rates at which this occurred these observations indicate that the antarctic ice sheets provided a significant contribution to the sea level rise at a rate that was approximately synchronous with the melting of the laurentide ice sheet except for the interval 9000 6000 yr ago when it may have lagged behind minor melting of the antarctic ice sheet appears to have continued throughout the late holocene differential observations of the late holocene sea level change recorded at sites in the same region are particularly useful for estimating parameters describing the earth s non elastic response to surface loading the effective parameters used here are a lithospheric thickness an upper mantle viscosity and a lower mantle viscosity describing the response below 670 km depth with the observations used here it is not possible to separate the lithospheric thickness h from the upper mantle viscosity and the viscosity results are based on the assumption that 50 h 100 km neither do these observations provide a resolution of the depth dependence of viscosity in the upper mantle and the resulting estimates are effective parameters only differential sea levels along continental margins and along the shores of large gulfs and bays constrain the effective upper mantle viscosity to be about 1 2 1020 pa s 1 while differential values from islands of different sizes are suggestive of a somewhat lower value the lower mantle taken to be below 670 km depth viscosity is about two orders of magnitude greater than this the estimated ice and rheological models explain many of the holocene sea level observations throughout the australian and southern pacific region the tilting of continental margins as exemplified by observations of variable amplitudes of the holocene high stands and the variable times at which sea levels first reached their present level along the north queensland coast and great barrier reef is well represented by these models differential holocene sea levels observed along the narrow spencers gulf of south australia are also well explained by the models and in neither case is it necessary to invoke tectonic motions predicted late holocene sea levels at small to medium sized islands are characterized by small amplitude high stands that reached their maximum values about 4000 2000 yr ago consistent with observations from the society cook and tuamotu islands

will earthquake activity in eastern canada increase in the next few thousand years knowledge of whether earthquake activity will increase in the next few thousand years is important for the planning of nuclear waste repositories assuming that fault instability portends earthquak

early response of the northeast subarctic pacific plankton assemblage to volcanic ash fertilization fe poor water collected at sta p20 in the gulf of alaska in june 2011 was enriched with different concentrations of volcanic ash 0 12 1 2 and 10 mg l 1 from two subduction zone volcanoes kasatochi and chaiten and incubated onboard under in situ conditions for 6 d the experimental setup also included a control no addition and a positive control addition of 0 6 nmol l 1 feso4 following a 4 d lag period there were increases in carbon fixation rates up to a factor of 10 and chlorophyll a concentrations up to a factor of 3 in the positive control and in the ash enriched 1 2 and 10 mg l 1 treatments diatoms dominated at the end of the incubations but cyanobacteria dinoflagellates pelagophytes and haptophytes were also stimulated by the presence of ash deposition of 1 mg ash l 1 which is in the lower range of those estimated to have caused the august 2008 bloom following the eruption of the kasatochi volcano in the aleutian islands would suffice to trigger a bloom in the gulf of alaska

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

initial value and modal approaches for transient viscoelastic responses with complex viscosity profiles summary an initial value theory is proposed for studying viscoelastic responses of compressible earth models with complex viscosity profiles continuous spectra are caused by both compressibility and the viscosity stratification and they can cause numerical difficulties for the modal approach we have studied the different responses for various types of viscosity profiles and found that there are differences in the responses at short wavelengths for viscosity profiles with sharp low viscosity zones we found that many modes are required to match the full initial value solutions in the presence of sharp viscosity stratifications in both the upper and lower mantle the normal mode approach is best suited for simple layered models long wavelengths and timescales greater than several thousand years while the initial value approach is indispensable in treating short timescale problems with sharp low viscosity zones in the upper mantle and viscosity stratification in the lower mantle

vhf radar observations of the dynamics of the summer polar mesopause region observations of the dynamics of the mesopause region over the norwegian island of andoya 69 n 16 e made during 2 days in summer 1987 using the mobile sousy vhf 53 5 mhz doppler radar are reported measurements of the three mean wind components and the 12 and 24 hour horizontal tidal components as well as preliminary estimates of the component of the total upward flux of density normalized horizontal momentum in the sw ne plane and the anisotropy of the horizontal wind are presented for the 83 to 90 km height interval values of the component of the density normalized flux in the sw ne plane for the 10 min to 48 hour period range after the 12 and 24 hour tidal components have been removed vary from 1 m2 s 2 near 84 km to 16 m2 s 2 near 89 km the major contribution is from the 10 min to 12 hour period range the corresponding mean flow acceleration takes a mean value of 136 m s 1 d 1over the 83 90 km height interval the anisotropy is considerable and may suggest the presence of standing waves

mantle stratification and long term polar wander there is accumulating evidence1 2 that the average location of the earth s magnetic pole has moved relative to the deep mantle over the past 200 myr this phenomenon of true polar wander may be caused by the advance and retreat of ice sheets3 or by mass redistribution in the earth s interior due to changes in the pattern of mantle convection4 5 new analyses1 2 6 of polar wander data show a significant shift of the pole in the late cretaceous but this period is thought to have been too warm for glaciation to have occurred thus most of the true polar wander must be due to mass movements in the mantle here we show by analysing the appropriate equations for polar wander that both viscosity and chemical stratification in the mantle are important in determining the rate of polar wander the dynamical effects of chemical stratification and high viscosity 1023 poise in the lower mantle are to decrease polar wander speed to a level consistent with the averaged velocities inferred from palaeomagnetic data whole mantle convection models without non adiabatic density jumps at 670 km depth and with viscosities less than 1023 p would produce wandering rates in excess of 5 per myr

mapping the lower mantle determination of lateral heterogeneity in p velocity up to degree and order 6 the data from the international seismological centre bulletins for the years 1964 1979 are used to derive a three dimensional model of lateral variations of the p velocity in the lower mantle unlike previous studies in which such perturbations were represented by a blocklike parameterization we seek the solution in the form r k 0k l 0l m 0l k r kalm cos m kblm sin m plm cos some 500 000 travel time residuals for teleseismic distances from 5000 earthquakes are used in an iterative procedure to derive the coefficients a and b the maximum k and l used are 4 and 6 respectively in each iteration in addition to solving for the three dimensional structure we also relocate all the earthquakes perturb the average travel time curve and determine station corrections for over 1000 stations particular attention is given to the problem of weighting the individual observations in order to avoid as much as possible the bias due to the uneven distribution of sources and receivers the resulting model shows a high level of perturbations just below the 670 km discontinuity and just above the core mantle boundary where the maximum perturbations reach 1 1 5 of the average velocity even for this highly smoothed model at a depth of 2000 km the rms perturbations are some 3 to 4 times lower than at either of the two boundaries the model predicts well the large scale pattern of observed travel time residuals for various source regions except for the distinct effects of the subduction zones the most striking large scale three dimensional feature of the model is a ring of high velocities circumscribing the pacific basin within a depth range from 1000 km to the core mantle boundary if the proportionality between perturbations in the velocity and temperature is assumed the derived pattern of anomalies could be consistent with a two boundary layer model of convection in the lower mantle some caution however should be exercised in interpreting the results near the 670 km discontinuity where the resolving power of our data set is limited by the lack of adequate control of the upper mantle structure at the same time a significant global increase in the level of perturbations in the lowermost mantle is established beyond question the results of dziewonski et al 1977 on the correlation between the equivalent geoid and the observed geoid for the angular order numbers 2 and 3 is confirmed but there is no significant correlation for degrees from 4 to 6 these calculations are made for the rigid earth richards and hager this issue show that introduction of finite viscosity can change the sign of a geoid anomaly in addition recent studies of lateral variations in the shear velocity reveal the presence of large scale anomalies in the upper mantle thus while the inferences made here with respect to the origin of the gravest orders of the geopotential field are not conclusive they point the way in which results from seismology can be used to address some of the basic questions of geodynamics

layered mantle convection a model for geoid and topography the long wavelength geoid and topography are dynamic effects of a convecting mantle the long wavelength geoid of the earth is controlled by density variations in the mantle and has been explained by circulation models involving whole mantle flow however the relationship of long wavelength topography to mantle circulation has been a puzzling problem in geodynamics we show that the dynamic topography is mainly due to density variations in the upper mantle even after the effects of lithospheric cooling and crustal thickness variation are taken into account layered mantle convection with a shallow origin for surface dynamic topography is consistent with the spectrum small amplitude and pattern of the topography layered mantle convection with a barrier about 250 km deeper than the 670 km phase boundary provides a self consistent geodynamic model for the amplitude and pattern of both the long wavelength geoid and surface topography

on transient rheology and glacial isostasy theoretical predictions from analytical self gravitating layered earth models with maxwell burgers body and standard linear solid rheologies are employed to explore the extent to which transient creep may have influenced the inference of long term mantle viscosity from analysis of the various signatures derived from the process of glacial isostatic adjustment the compatibility of the various models is checked by their ability to fit simultaneously the relative sea level free air gravity and the rotational data we conclude that transient creep cannot operate throughout the entire mantle on account of the extreme sensitivity of the uplift and gravity data to any significant departures of the steady upper mantle viscosity from 1021 pa s transient rheology would most likely be unimportant in the upper mantle restricting transient rheology to the lower mantle we find that models with an adiabatic density jump at 670 km depth are not satisfactory because of their inability to produce sufficiently large gravity anomalies and polar wander speeds our calculations show that the preferred models with transient creep in the lower mantle are chemically stratified with a nonadiabatic density change of around 10 the long term viscosities in the lower mantle associated with compatible transient rheological models lie between 5 1022 and 1023 pa s these values are consistent with viscosity estimates derived from correlation of long wavelength geoid anomalies and the density perturbations in the deep mantle

mantle heterogeneities geoid and plate motion a monte carlo inversion seismic tomography in both the upper and the lower mantle as well as subducting oceanic slabs defined by seismicity has been translated into density heterogeneities to generate models of mantle circulation these models can predict both the surface velocities and the geoid which can be compared with plate tectonics and gravity data a given model is specified by 6 parameters related to the viscosities of 3 mantle layers and the absolute amplitudes of density variations in the upper and lower mantle as well as in the slabs the values of these parameters are chosen at random within an acceptable range each model is submitted to an appropriate test comparing observations with predictions the results of the most successful models selected by this monte carlo inversion are displayed they yield preferred mantle viscosity structures exhibiting large variations at depth with a physical interface between upper and lower mantle i e with the possibility for the circulation to penetrate the 650 km discontinuity two classes of viscosity profiles stand out the first one implies a regular increase of the viscosity in the sublithospheric mantle with reasonable values for the density parameters the second one is unexpected in the sense that it predicts a very stiff bottom for the upper mantle it also requires vanishingly small amplitudes for the upper mantle density heterogeneities defined by tomography which would thus have to be of lithological rather than thermal origin with a chemical interface at 650 km the outcome is very similar the same two classes of viscosity structures do yield a satisfactory geoid prediction however only the class of models with a stiff layer at midmantle depths predicts acceptable surface velocities altogether the best models out of some 60 000 which have been tested only explain one third of the geoid and two thirds of the surface divergence for spherical harmonic degrees 1 to 6 nevertheless the main features of these two observed patterns are present in the computed maps and 4 out of 6 correlation coefficients lie close to the 90 confidence level this is true for the geoid as well as for the surface divergence of the displacement velocity however as the internal viscosity structure has been assumed to have spherical symmetry the rotational component of the surface velocities cannot be predicted

is the lower mantle rheology newtonian today the issue concerning the dominant creep mechanism in mantle convection has been studied numerically with a rheology incorporating both linear and non linear elements we have employed a rheology suggested by the recently obtained melting temperature of lower mantle constituents for an effective rayleigh number between 105 and 106 this type of strongly temperature and pressure dependent rheology induces a circulation in the earth s mantle which is characterized by a sluggish type of flow in the lower mantle with a prevailing newtonian character while in the upper mantle much more vigorous flows are developed in a predominantly non newtonian rheology

models of the core mantle boundary and the travel times of internally reflected core phases topographic maps of the core mantle boundary cmb obtained by inverting different subsets of international seismological centre arrival time data involving either waves refracted through the boundary pkp or waves reflected from the topside pcp do not agree we report on the implications for models of cmb structure when data from bottomside reflections notably pkkp are included in the inversion travel time residuals with respect to a reference earth model prem of all core phases are positively biased but the mean residual is different for different phases and must be accounted for prior to inversion the new data illustrate that models of cmb topography can explain only a small fraction of the data variance the new cmb model involves both topography and a laterally varying boundary layer the cmb itself is smooth compared to previous results

investigation of upper mantle structure under young regions of the southeast pacific using long period rayleigh waves abstract long period rayleigh wave group velocity dispersion curves are presented for paths across the east pacific ocean the records are from the ipg site at pamatai tahiti and the ida station at nana for the first time direct path observations of group velocities up to 300 s have been obtained this study shows that for young oceanic regions group velocities are low even at long periods the observations are interpreted in terms of an s wave velocity model by a generalized inversion scheme in the models for young ages the low velocity zone under the asthenosphere lid is well developed with a strong velocity gradient at the bottom of this zone followed at larger depths by a plateau representing a lower velocity zone and a marked gradient at 400 km

plate tectonics and aspherical earth structure the importance of poloidal toroidal coupling we derive poloidal flow internal loading green functions for incompressible fluid models of the mantle consisting either of a single constant viscosity spherical shell or of two adjacent spherical shells having different viscosities from the green function we obtain kernels connecting the surface divergence the geoid and the surface topography fields to the lateral density heterogeneity inferred on the basis of the application of seismic tomographic imaging techniques and with these kernels we argue that both the surface divergence and geoid fields consisting of harmonic degrees 2 5 may be reasonably fit with only a factor of 8 viscosity increase at a depth of 1200 km we point out however that the coupling from poloidal to toroidal flow which is required to understand surface velocity spectra may allow these observations to be understood in terms of a viscosity increase at depth which is smaller than required by nonhydrostatic geoid data in the context of pure poloidal models using the observed kinetic energy in the surface plate motions as a constraint allows us to infer a value for the steady state upper mantle viscosity of 2 0 0 5 1021 pa s we comment on the implication of the difference between this number and the somewhat lower value 1 1021 pa s which has been derived on the basis of analyses of signatures of the glacial isostatic adjustment process

a new method for the computation of global viscoelastic post seismic deformation in a realistic earth model i vertical displacement and gravity variation summary we introduce a new method by which to compute global post seismic deformation psd in a spherically symmetric self gravitating viscoelastic earth model previous methods are based on simplified earth models that neglect compressibility and or the continuous variation of the radial structure of earth this is because the previous mode summation technique cannot avoid intrinsic numerical difficulties caused by the innumerable poles that appear in a realistic earth model that considers such effects in contrast the proposed method enables both of these effects to be taken into account simultaneously we carry out numerical inverse laplace integration which allows evaluation of the contribution from all of the innumerable modes of the realistic earth model using this method a complete set of green s functions is obtained including functions of the time variation of the displacement gravity change and the geoid height change at the surface for strike slip dip slip horizontal and vertical tensile point dislocations as an earth model we employ the preliminary reference earth model prem and a convex viscosity profile further we investigate the effects of fine layering of the viscoelastic structure and compressibility on a time series of psd using the green s function for a dip slip fault the result indicates that the effect of increasing number of layers is saturated at several tens of layers even when compressibility is taken into account and that the effect of compressibility is detectable with modern observational techniques for a shallower large earthquake mw 8 as an application the psd due to the sumatra andaman islands earthquake mw 9 3 is estimated we show that the rate of post seismic vertical displacement and gravity change is possibly detected in the far field where the epicentral distance exceeds 400 km

lower mantle thermal structure deduced from seismic tomography mineral physics and numerical modelling abstract the long wavelength thermal anomalies in the lower mantle have been mapped out using several seismic tomographic models in conjunction with thermodynamic parameters derived from high pressure mineral physics experiments these parameters are the depth variations of thermal expansivity and of the proportionality factor between changes in density and seismic velocity the giant plume like structures in the lower mantle under the pacific ocean and africa have outer fringes with thermal anomalies around 300 400 k but very high temperatures are found in the center of the plumes near the base of the core mantle boundary these extreme values can exceed 1500 k and may reflect large hot thermal anomalies in the lower mantle which are supported by recent measurements of high melting temperatures of perovskite and iron extremely cold anomalies around 1500 k are found for anomalies in the deep mantle around the pacific rim and under south america numerical simulations show that large negative thermal anomalies of this magnitude can be produced in the lower mantle following a catastrophic flushing event cold anomalies in the mid lower mantle have modest magnitudes of around 500 k a correlation pattern exists between the present day locations of cold masses in the lower mantle and the sites of past subduction since the cretaceous results from correlation analysis show that the slab mass flux in the lower mantle did not conform to a steady state nature but exhibited time dependent behavior

rotational instabilities of the earth induced by mantle density anomalies from the seismic tomographic models the authors estimate the offset between the geographic pole and the main inertial axis this offset which is large for a static earth is drastically reduced when the dynamical behavior of the mantle is taken into account for a mantle stratification with a viscosity increase by a factor of about 50 this offset is comparable to what can be estimated on the basis of remanent post glacial depression the authors also discuss the rate of polar drift induced by both internal mass heterogeneities and pleistocenic deglaciation contributions from internal sources and ice sheets disintegration are consistent with the observed secular trend for high viscosity contrasts the results are quite insensitive to the phase transition of chemical nature of the upper lower mantle interface this study shows that there is no contradiction between the tomographic models of the earth and their implications for the rotational stability of our planet the findings also agree with recent analyses of paleomagnetic data

geoid anomalies in a dynamic earth in order to obtain a dynamically consistent relationship between the geoid and the earth s response to internal buoyancy forces we have calculated potential and surface deformation love numbers for internal loading these quantities depend on the depth and harmonic degree of loading they can be integrated as green functions to obtain the dynamic response due to an arbitrary distribution of internal density contrasts spherically symmetric self gravitating flow models are constructed for a variety of radial newtonian viscosity variations and flow configurations including both whole mantle and layered convection we demonstrate that boundary deformation due to internal loading reaches its equilibrium value on the same time scale as postglacial rebound much less than the time scale for significant change in the convective flow pattern by calculating relaxation times for a series of spherically symmetric viscous earth models for uniform mantle viscosity the geoid signature due to boundary deformations is larger than that due to internal loads resulting in net negative geoid anomalies for positive density contrasts geoid anomalies from intermediate wavelength density contrasts are amplified by up to an order of magnitude geoid anomalies are primarily the result of density contrasts in the interior of convecting layers density contrasts near layer boundaries are almost completely compensated layered mantle convection results in smaller geoid anomalies than mantle wide flow for a given density contrast viscosity stratification leads to more complicated spectral signatures because of the sensitivity of the dynamic response functions to model parameters forward models for the geoid can be used to combine several sources of geophysical data e g subducted slab locations seismic velocity anomalies surface topography to constrain better the structure and viscosity of the mantle