coordinated hydrological regimes in the indo pacific region during the past two millennia instrumental data suggest that major shifts in tropical pacific atmospheric dynamics and hydrology have occurred within the past century potentially in response to anthropogenic warming to better understand these trends we use the hydrogen isotopic ratios of terrestrial higher plant leaf waxes delta d wax in marine sediments from southwest sulawesi indonesia to compile a detailed reconstruction of central indo pacific warm pool ipwp hydrologic variability spanning most of the last two millennia our paleodata are highly correlated with a monsoon reconstruction from southeast asia indicating that intervals of strong east asian summer monsoon easm activity are associated with a weaker indonesian monsoon im furthermore the centennial scale oscillations in our data follow known changes in northern hemisphere climate e g the little ice age and medieval warm period implying a dynamic link between northern hemisphere temperatures and ipwp hydrology the inverse relationship between the easm and im suggests that migrations of the intertropical convergence zone and associated changes in monsoon strength caused synoptic hydrologic shifts in the ipwp throughout most of the past two millennia

mercury s global contraction much greater than earlier estimates observations of compressional structures on mercury have fallen short of accommodating the global contraction that is required owing to cooling of the planet s interior mapping of folds and faults across mercury s surface using messenger spacecraft images reveals deformation consistent with a planet that has contracted radially as much as seven kilometres over its history

initiation of subduction zones along transform and accreting plate boundaries triple junction evolution and forearc spreading centres implications for ophiolitic geology and obduction

mantle convection with a brittle lithosphere thoughts on the global tectonic styles of the earth and venus summary plates are an integral part of the convection system in the fluid mantle but plate boundaries are the product of brittle faulting and plate motions are strongly influenced by the existence of such faults the conditions for plate tectonics are studied by considering brittle behaviour using byerlee s law to limit the maximum stress in the lithosphere in a mantle convection model with temperature dependent viscosity when the yield stress is high convection is confined below a thick stagnant lithosphere at low yield stress brittle deformation mobilizes the lithosphere which becomes a part of the overall circulation surface deformation occurs in localized regions close to upwellings and downwellings in the system at intermediate levels of the yield stress there is a cycling between these two states thick lithosphere episodically mobilizes and collapses into the interior before reforming the mobile lid regime resembles convection of a fluid with temperature dependent viscosity and the boundary layer scalings are found to be analogous this regime has a well defined nusselt number rayleigh number relationship which is in good agreement with scaling theory the surface velocity is nearly independent of the yield stress indicating that the plate motion is resisted by viscous stresses in the mantle analysis suggests that mobilization of the earth s lithosphere can occur if the friction coeycient in the lithosphere is less than 0 03 0 13 lower than laboratory values but consistent with seismic field studies on venus the friction coeycient may be high as a result of the dry conditions and brittle mobilization of the lithosphere would then be episodic and catastrophic

geological consequences of super sized earths 1 the discovery of terrestrial scale extrasolar planets and their calculated abundance in the galaxy has prompted speculation on their surface conditions and thermal structure both are dependent on the tectonic regime of a planet which is itself a function of the balance between driving forces and the resistive strength of the lithosphere here we use mantle convection simulations to show that simply increasing planetary radius acts to decrease the ratio of driving to resisting stresses and thus super sized earths are likely to be in an episodic or stagnant lid regime this effect is robust when associated increases in gravity are included as the more dominant effect is increased fault strength rather than greater buoyancy forces the thermo tectonic evolution of large terrestrial planets is more complex than often assumed and this has implications for the surface and conditions habitability of such worlds

on the initiation of subduction estimates of shear resistance associated with lithospheric thrusting and convergence represent lower bounds on the force necessary to promote trench formation three environments proposed as preferential sites of incipient subduction are investigated passive continental margins transform faults fracture zones and extinct ridges none of these are predicted to convert into subduction zones simply by the accumulation of local gravitational stresses subduction cannot initiate through the foundering of dense oceanic lithosphere immediately adjacent to passive continental margins the attempted subduction of buoyant material at a mature trench can result in large compressional forces in both subducting and overriding plates this is the only tectonic force sufficient to trigger the nucleation of a new subduction zone the ubiquitous distribution of transform faults and fracture zones combined with the common proximity of these features to mature subduction complexes suggests that they may represent the most likely sites of trench formation if they are even marginally weaker than normal oceanic lithosphere

role of a low viscosity zone in stabilizing plate tectonics implications for comparative terrestrial planetology 1 earth s near surface layer its lithosphere is broken into quasi rigid plates that form the upper thermal boundary layer for mantle convection since the discovery of plate tectonics it has been widely conjectured but only recently demonstrated that this peculiar style of convection may be facilitated by an upper mantle low viscosity zone lvz over which the plates glide easily the lvz or asthenosphere concept dates from 19th century investigations of isostatic support of mountain belts and is supported by modern evidence for a seismic low velocity zone and by studies of postglacial rebound and dynamic compensation of the earth s gravity field here we show in both two dimensional 2 d cartesian and 3 d spherical earth models that combining a pronounced lvz and a plastic yield stress to allow localized weakening of the cold thermal boundary layer results in a distinctly plate tectonic style of convection with 30 toroidal surface motion for the 3 d case recycling of water into the upper mantle at subduction zones is a plausible cause of earth s lvz whereas venus is dry and lacks both an lvz and plate tectonics

steady state secondary convection beneath lithospheric plates with temperature and pressure dependent viscosity we have constructed a thermomechanical model for upper mantle convection such that the thickness and structure of the lithosphere are determined self consistently by the heat transported by convection in this study of the interaction between the lithosphere and upper mantle strongly temperature and pressure dependent rheologies for both newtonian and non newtonian creep mechanisms are employed for a strictly temperature dependent rheology an insignificant amount of heat less than 12 5 mw m2 can be transported convectively for an interior viscosity 0 1021 pas compatible with postglacial rebound on the other hand for similar values of the interior viscosity steady heat fluxes between 20 and 40 mw m2 are produced by introducing pressure dependence into the rheology for the temperature and pressure dependent flow law the horizontally averaged interior temperature displays very little variation with the amount of heat evacuated once all of the rheological parameters are fixed this finding may have important ramifications for parameterized convection we employ both the single mode mean field approximation and the complete two dimensional equations using finite elements to obtain solutions for the various types of rheologies from evaluating the geophysically relevant observables such as topography free air gravity anomalies surface heat flow and stress fields in the lithosphere we find that the lateral variations of these quantities predicted by a non newtonian rheology are much smaller than those derived from a linear rheology these results suggest that surface variations of geophysical observables are more compatible with a non newtonian rheology in the upper mantle

convection scaling and subduction on earth and super earths super earths are the smallest class of discovered extra solar planets owing to their relatively small mass some might resemble earth and perhaps be habitable because of the connection to habitability through thermal evolution we investigate the tectonic regime of massive terrestrial planets two independent studies o neill c lenardic a oct 2007 geological consequences of super sized earths grl 34 l19204 valencia d o connell r j sasselov d d nov 2007a inevitability of plate tectonics on super earths apjl 670 l45 l48 have reached opposing conclusions about the likelihood of plate tectonics here we offer possible reasons for the discordant findings and address three key aspects of sustaining plate tectonics deformation on faults negative buoyancy and energy dissipation during subduction we show that in general the ratio of driving force to plate resistance increases with planetary mass this is a consequence of increasing convective stresses thinning plates and similar plate structure we conclude that even though the strength of dry and wet of faults increases for massive terrestrial planets the convective stress increases even more allowing deformation to take place also despite shorter timescales for plate cooling rocky super earths achieve negative buoyancy at subduction zones finally by investigating the effects of energy dissipation during subduction we find that massive terrestrial planets dissipate less energy during subduction and hence provide a positive feedback to sustain active lid tectonics in conclusion rocky super earths have more favorable conditions than earth for the subduction of plates and hence for sustaining plate tectonics

boundary layer theory and subduction numerical models of thermally activated convective flow in earth s mantle do not resemble active plate tectonics because of their inability to model successfully the process of subduction other than by the inclusion of artificial weak zones here we show using a boundary layer argument how the rigid lid style of convection favored by thermoviscous fluids leads to lithospheric stresses which may realistically exceed the yield stress and thus cause subduction to occur through the visoc plastic failure of lithospheric rock an explicit criterion for the failure of the lid is given which is sensitive to the internal viscosity r below the lid for numbers appropriate to earth s mantle this criterion is approximately r a 10 21 pa s

characteristics of winter cyclone activity in the northern north atlantic insights from observations and regional modeling 1 fields from the ncep ncar reanalysis are used in conjunction with simulations from polar mm5 to examine regional processes contributing to cyclone development in the sub arctic arctic sector of the north atlantic strong deepening of systems 6 hpa 6 hours is especially common south of the tip of greenland and in the vicinity of the mean icelandic low il a secondary maximum is located just south of svalbard composite analyses of strongly deepening systems in these areas point to the influence of pronounced preexisting low level temperature gradients linked to ocean currents the sea ice margin and ocean land contrasts the area of preferred deepening near svalbard corresponds to the northernmost penetration of open water in the north atlantic impacts of greenland s orography are prominent in the synoptic complexity of the region typical synoptic situations examined as case studies include splitting or bifurcation of cyclones at the southern tip of greenland orographic cyclogenesis in the lee of greenland at the location of the mean il and deepening of preexisting systems near the il orographic influences are clearly captured in polar mm5 simulations for which control runs of the three cases are compared to runs for which the orography of greenland is removed

plate motions and the viscosity structure of the mantle insights from numerical modelling abstract the viscosity structure of the earth s mantle is likely to play an important role with respect to the motion of the lithospheric plates a zone of low viscosity in the upper mantle as postulated by seismological studies and also by investigations of postglacial uplift has been proposed to facilitate plate motion another inferred feature of the viscosity profile in the mantle is a significant increase at greater depth we have employed a three dimensional mantle convection model to explore the relation between the appearance of plate motion and the viscosity structure of the mantle beneath the model allows for a complex fluid rheology strong temperature pressure and stress dependence of the viscosity and can so account for the self consistent formation of plates at the surface of the convecting mantle plate like motion appears only in a limited parameter range and we demonstrate that this appearance is closely linked to a specific viscosity structure the viscosity is characterized by a local minimum at shallow depth low viscosity zone and by a local maximum at greater depth self consistent plate motion is only observed if such a viscosity structure exists in order to determine the functional relationship between the viscosity structure and plate motion we also employ an approach in which a certain viscosity profile is prescribed these studies again reveal the importance of the local extrema in the viscosity profile taking into account an increased interior viscosity also obtained by the pressure dependence we find stable plate motion as a general result we observe that the balance of the temperature and stress dependence of the viscosity is necessary to obtain plate motion at the same time we find a viscosity profile showing a minimum at shallow depth and a maximum at greater depth

why does venus lack a magnetic field venus and earth have similar radii and estimated bulk compositions and both have an iron core that is at least partially liquid however despite these similarities venus lacks an appreciable dipolar magnetic field here i examine the hypothesis that this absence is due to venus9s also lacking plate tectonics for the past 0 5 b y the generation of a global magnetic field requires core convection which in turn requires extraction of heat from the core into the overlying mantle plate tectonics cools earth9s mantle on the basis of elastic thickness estimates and convection models it is argued here that the mantle temperature on venus is currently increasing this heating will reduce the heat flux out of the core to zero over 1 b y halting core convection and magnetic field generation if plate tectonics was operating on venus prior to ca 0 5 ga a magnetic field may also have existed on earth the geodynamo may be a consequence of plate tectonics this connection between near surface processes and core magnetism may also be relevant to the generation of magnetic fields on mars mercury and ganymede

constraints on plate tectonics initiation from scaling laws for single cell convection abstract the earth is the only planet known to have plate tectonics while other planets are covered with a stagnant lid on the earth the initiation of subduction which is thought to be the fundamental process for plate tectonics initiation is caused not only by the negative buoyancy of the lithosphere but also by the forces from plate motions however for planets which do not have plate tectonics the very first episode of lithospheric failure has to be caused by forces other than plate motions sublithospheric convection has been proposed as a possible mechanism that provides lithospheric instability through inducing stresses in the lithosphere and lithospheric failure can occur when the yield stress is below a critical value we test the applicability of scaling laws for the critical yield stress obtained in single cell convection simulations to strongly time dependent multi cell systems we show that with an appropriate choice of characteristic aspect ratio for the convective system the scaling laws from single cell simulations can be used to evaluate the conditions on the terrestrial planets in the inner solar system for plate tectonics to exist in agreement with previous studies the estimated values for critical yield stress and coefficient of friction are much lower than the expected values for the earth s lithosphere

orographic influence of east greenland on a polar low over the denmark strait we present a numerical study of a polar low which hit western iceland in january 2007 with heavy snowfall and mean wind speeds exceeding 20 m s 1 in several locations the operational models at the time captured the polar low formation rather well but there was a large spread in their predictions of the subsequent evolution and track of the polar low the objective of this study is to investigate possible orographic forcing from greenland as a trigger for the polar low development in addition to an analysis of surface observations and satellite imagery sensitivity studies using hirlam were carried out with various degradations of greenland s orography as well as with modifications to the sea surface temperature sst surface roughness and the data assimilation scheme despite the presence of an upper level trough and weak static stability in all the simulations the polar low development was found to be very sensitive to the presence of the high mountains of eastern greenland whereas the control run captured well the main features of the polar low simulations with parts of east greenland s orography removed gave a southward displaced polar low which moved rapidly eastward resulting in substantially underestimated near surface winds and snowfall amounts setting the orographic heights over all of greenland to zero led to the complete disappearance of the polar low on the other hand artificially increasing the sst by 4 k in the denmark strait reducing the orographic roughness or replacing the four dimensional variational assimilation scheme 4d var by 3d var had only a small effect on the polar low we suggest that hitherto unreported interactions between the high mountains of east greenland and polar low development over the denmark strait may be more important for polar low formation than katabatic flow from valleys in east greenland that was highlighted in earlier studies copyright 2011 royal meteorological society

mechanisms for the generation of plate tectonics by two phase grain damage and pinning abstract weakening and shear localization in the lithosphere are essential ingredients for understanding how and whether plate tectonics is generated from mantle convection on terrestrial planets we present a new theoretical model for the mechanism of lithospheric weakening and shear localization and hence plate generation through damage grain evolution and zener pinning in two phase polycrystalline lithospheric rocks grain size evolves through the competition of coarsening which drives grain growth with damage which drives grain reduction however in a two phase medium the interface between phases induces zener pinning which impedes grain growth and facilitates damage the size of the pinning surfaces is given by the roughness of the interface and damage to the interface causes smaller pinning surfaces which in turn drive down the grain size forcing the rheology into the grain size dependent diffusion creep regime this process allows damage and rheological weakening to co exist which is normally considered impossible in single phase assemblages moreover pinning greatly inhibits grain growth and shear zone healing which is much faster in single phase materials hence the resulting shear localization is rapid less than 1 myr but the healing time for a dormant weak zone is very slow greater than 100 myrs these effects therefore permit rapidly forming and long lived plate boundaries the model therefore provides a key ingredient and predictive theory for the generation of plate tectonics on earth and other planets

the preservation of seismic anisotropy in the earth s mantle during diffusion creep summary seismic anisotropy in the earth particularly in the mantle is commonly interpreted as the result of solid state deformation by dislocation creep that induces a lattice preferred orientation lpo diffusion creep operates where stress levels are lower and or grain sizes smaller it is often assumed that diffusion creep induces grain rotations that eventually destroy any existing lpo a new numerical test of this assumption shows that it is not necessarily the case diffusion creep will create some relative grain rotations but rotation rates decrease through time hence when microstructural change due to diffusion creep dominates that due to grain growth defined here as type p behaviour the converse being type o behaviour the model indicates that lpo will be weakened but preserved for a variety of strain paths including both pure and simple shear one measure of anisotropy is the proportional difference in shear wave velocities for different polarization vectors avs a model olivine microstructure with equant grains and initial maximum avs of 10 0 percent has this value reduced to 6 7 per cent when rotational steady state is attained other models with different initial maximum avs values exhibit final maximum avs values more than half the initial values if the grains are initially elongate by a factor of 2 maximum avs is reduced just slightly to 8 5 per cent thus when grain growth plays a subordinate role to the deformation diffusion creep weakens seismic anisotropy by a factor of less than 2 using maximum avs as a measure and olivine as an example consequently the link between seismic anisotropy and deformation mechanism in the mantle requires reappraisal regions with lpo may comprise material which once deformed by dislocation creep but is now deforming by diffusion creep in a rotational steady state

perturbation analysis of convective instability of oceanic lithosphere and initiation of subduction zones a linear perturbation analysis is made on convective instability of the oceanic lithosphere to understand the occurrence of subduction on the earth the oceanic lithosphere is modeled by the upper part of top thermal boundary layer in a convecting newtonian temperature dependent viscosity fluid thickening of oceanic lithosphere with its age due to cooling at the surface is modeled by thickening of the thermal boundary layer with time its thickening rate l is a measure of the rate with which temperature in the thermal boundary layer decreases with time and hence is a measure of the rate with which viscosity in the thermal boundary layer increases with time the linear perturbation analysis shows the possible excitation of plural eigenmodes to give rise to an instability in the thermal boundary layer depending on a threshold value rc of the local rayleigh number r of the thermal boundary layer further rc is shown to be a linearly increasing function of l which is treated as a free parameter from estimates of r and l appropriate for the major oceanic lithosphere denoted as robs and l obs respectively i find that rc l 0 robs rc l l obs the flow pattern of mantle convection depends strongly on which eigenmodes are excited by the convective instability in thermal boundary layer 1 when r is lower thanrc only the fundamental and the first higher modes at most are excited the resulting mantle convection manifests itself as the secondary convection 2 when r is higher than rc the modes higher than the first higher mode are also excited the resulting mantle convection involves lithospheric instability leading to large local deformation in the oceanic lithosphere and finally to subduction of the oceanic lithosphere this result implies that 1 subduction does not begin newly in the major oceanic lithosphere at present since lithospheric instability is suppressed by thickening of the major oceanic lithosphere neither by its high viscosity nor by its buoyancy 2 subduction will newly begin in the major oceanic lithosphere due to lithospheric instability in future when the oceanic lithosphere becomes so old that its thickening rate l decreases to sufficiently small value the latter result suggests that lithospheric instability is an important agent that maintains multiplate tectonics on the earth

the effect of rheological parameters on plate behaviour in a self consistent model of mantle convection abstract in a three dimensional numerical model of mantle convection we have explored the self consistent formation and evolution of plates this was done by investigating the influence of temperature stress and pressure dependence of the viscosity on plate like behaviour further the role of a depth dependent thermal expansivity and of internal heating was examined first order self consistent plate tectonics can be obtained for temperature and stress dependent viscosity we observe virtually undeformed moving plates and the phenomenon of trench migration however plate motion only appears during short intervals which are interrupted by long phases characterised by an immobile lid an improvement is achieved by adding pressure dependence of the viscosity and of the thermal expansivity to the model this leads to extended plates with continuous plate motion and the appearance of a low viscosity zone internal heating accelerates subduction and gives an increase in plate velocity the plates are otherwise left intact we find that the surface velocity is largely determined by subduction rather than by upwellings plate like behaviour is restricted to a narrow parameter window the window is bordered by an episodic or a stagnant lid type of flow on the one side and by a domain in which convection is akin to constant viscosity flow on the other side with these different convective regimes the behaviour of various terrestrial planets can be described self consistently occasionally even a change from the plate regime to the stagnant lid regime occurs such a type of transition has possibly changed the style of tectonics on some terrestrial planets such as mars

subduction zone infancy examples from the eocene izu bonin mariana and jurassic california arcs a new model for the earliest stages in the evolution of subduction zones is developed from recent geologic studies of the izu bonin mariana ibm arc system and then applied to late jurassic ophiolites of cailfornia the model accounts for several key observations about the earliest stages in the evolution of the ibm system 1 subduction nucleated along an active transform boundary which separated younger less dense lithosphere in the west from older more dense lithosphere to the east 2 initial arc magmatic activity occupied a much broader zone than existed later 3 initial magmatism extended up to the modern trench over a region now characterized by subnormal heat flow 4 early are magmatism was characterized by depleted tholeiitic and ultra depleted boninitic magmas indicating that melting was more extensive and involved more depleted mantle than is found anywhere else on earth 5 early arc magmatism was strongly extensional with crust forming in a manner similar to slow spreading ridges and 6 crust production rates were 120 to 180 km 3 km ma several times greater than for mature arc systems these observations require that the earliest stages of subduction involve rapid retreat of the trench we infer that this resulted from continuous subsidence of denser lithosphere along the transform fault this resulted in strong extension and thinning of younger more buoyant lithosphere to the west this extension was accompanied by the flow of water from the sinking oceanic lithosphere to the base of the extending lithosphere and the underlying asthenosphere addition of water and asthenospheric upwelling led to catastrophic melting which continued until lithosphere subsidence was replaced by lithospheric subduction application of the subduction zone infancy model to the late jurassic ophiolites of california provides a framework in which to understand the rapid formation of oceanic crust with strong arc affinities between the younger sierran magmatic arc and the franciscan subduction complex provides a mechanism for the formation and subsidence of the great valley forearc basin and explains the limited duration of high t high p metamorphism experienced by franciscan melanges

magellan and comparative planetology magellan data have been widely used to provide explanations for the tectonics and volcanism on venus however it is also important to compare our understanding of these processes on venus with the earth moon and mars one essential question concerning the behavior of these bodies is how does the heat get out internal heat generation and planetary cooling are the energy sources responsible for planetary tectonics and volcanism on the earth the plates of plate tectonics are the thermal boundary layers of mantle convection cells the subduction of these cold plates at ocean trenches contributes about 75 of the mantle heat flux it is now quite clear that active plate tectonics is not occurring on venus the surface of venus is a single rigid plate as is the case for the moon mars and mercury prior to magellan the general consensus was that venus was in a near uniformatarian state with conductive heat losses through the lithosphere balancing internal heat generation and secular cooling of the planet however crater studies from magellan images provide direct evidence for a catastrophic or episodic evolution of venus one explanation is that venus died tectonically 500 m y ago plate tectonics occurred during the early evolution of the planet but ceased as the radioactive heat production decreased an alternative explanation is that venus experiences episodic global subduction events two major questions concerning the earth remain largely unanswered particularly when compared with venus 1 why does the earth have continents 2 why does the earth have plate tectonics a major difference between the earth and venus is the presence of liquid water oceans during much of the evolution of the earth it has been hypothesized that water is an essential ingredient for the silicic volcanism that produces continental rocks on the earth and the absence of well defined continental masses on venus supports this hypothesis it has long been recognized that plate tectonics cannot evolve in time without internal deformation of plates plate evolution in time results in geometrical incompatibilities that require overlaps or gaps on the earth these incompatibilities lead to broad zones of continental deformation i e the western united states and china the absence of plate tectonics on venus suggests that rheologically soft continents are necessary for plate tectonics a reasonable conclusion is that the earth is the remarkable planet in terms of tectonics and volcanism not venus

a community benchmark for 2 d cartesian compressible convection in the earth s mantle benchmark comparisons are an essential tool to verify the accuracy and validity of computational approaches to mantle convection six 2 d cartesian compressible convection codes are compared for steady state constant and temperature dependent viscosity cases as well as time dependent constant viscosity cases in general we find good agreement between all codes when comparing average flow characteristics such as nusselt number and rms velocity at rayleigh numbers near 106 and dissipation numbers between 0 and 2 the results differ by approximately 1 per cent differences in discretization and use of finite volumes versus finite elements dominate the differences there is a small systematic difference between the use of the anelastic liquid approximation ala compared to that of the truncated ala in determining the onset of time dependence there was less agreement between the codes with a spread in the rayleigh number where the first bifurcation occurs ranging from 7 79 10 5 to 1 05 10 6

viscosity of oceanic asthenosphere inferred from remote triggering of earthquakes a sequence of large interplate earthquakes from 1952 to 1965 along the aleutian arc and kurile kamchatka trench released accumulated stresses along nearly the entire northern portion of the pacific plate boundary the postseismic stress evolution across the northern pacific and arctic basins calculated from a viscoelastic coupling model with an asthenospheric viscosity of 5 x 10 17 pascal seconds is consistent with triggering of oceanic intraplate earthquakes temporal patterns in seismicity at remote plate boundaries and space based geodetic measurements of anomalous velocity over an area 7000 by 7000 kilometers square during the 30 year period after the sequence

dynamics of subduction initiation at preexisting fault zones we show that a new subduction zone can initiate at a preexisting dipping fault zone with reasonable plate forces consistent with what is observed in the western pacific the dynamics of subduction initiation within a viscoelastic medium has been systematically explored with the finite element method we investigate the compression of oceanic lithosphere with both force and velocity boundary conditions and track the thermal structure in which heat is transported by both advection and diffusion the viscosity of the medium is non newtonian and temperature dependent we also examine the influence of a perfectly plastic yield stress a new method to model the preexisting weak zone as a graded fault with remeshing is introduced the first 400 km of plate convergence is examined during which time there are profound changes in plate boundary dynamics initially topography across the plate boundary is characterized by over 1000 m of uplift close to the trench subsequently a bathymetric depression develops on the overriding plate with an initial subsidence rate of 125 m myr for a slab subducting at 2 cm yr tectonic subsidence of the overriding plate is primarily dependent upon the depth of slab penetration for a plate force of 4 10 12 n m subduction initiates even with a fault shear stress of 3 mpa the resistance most consistent with seismological observations however subduction will not initiate when this fault shear stress is 5 mpa our results are consistent with plate reconstructions which predict the initiation of subduction across preexisting weak zones within oceanic basins

the tectonics of mercury the view after messenger s first flyby article i nfo during its first flyby of mercury messenger imaged many tectonic landforms most of which are contractional in nature and include lobate scarps high relief ridges and wrinkle ridges wrinkle ridges are found on the smooth plains that partially fill the interior and surround the exterior of the caloris basin and also on smooth plains that fill the interiors of smaller impact basins and larger craters messenger revealed a radial graben complex pantheon fossae nearlyco centered with the caloris basin pantheon fossae and a network of mostly basin concentric graben intheouterportionsofthecalorisbasininteriorformapatternofextensionnotseenelsewhereonmercury the first clear example of extensional faults outside of the caloris basin has been documented on smooth plains inside the peak ring of the relatively young raditladi basin a map of the distribution of tectonic landforms imaged by messenger and mariner 10 shows that lobate scarps are the most widespread type of feature estimates of the horizontal shortening associated with lobate scarps that crosscut and overthrust small impact craters imaged by messenger range from 1 to 3 km previously unrecognized lobate scarps detected in areas imaged by mariner 10 indicate thatearlierestimatesof contractional strainarelow anewestimate of theaverage arealcontractionalstrain 0 06 accommodated byscarps is at leastonethirdgreater thancomparableprevious estimates andcorresponds to a decrease in mercury s radius of at least 0 8 km since the end of heavy impact bombardment of the inner solar system these figures are lower bounds because of the likelihood that not all lobate scarps have been identified even in areas imaged to date crosscutting and embayment relations revealed by messenger suggest that lobate scarps formed before the end of smooth plains emplacement and continued to be active after the emplacement of the youngest smooth plains deposits relatively recent activity on lobate scarps may be the expression of slow but continuous radial contraction that accompanied cooling of mercury s interior and the growth of the planet s solid inner core published by elsevier b v

subduction origin on early earth a hypothesis i propose a hypothesis for the origin of subduction on early earth that is directly coupled with mantle upwelling and the formation of mafi c crust the hypothesis invokes spatial and temporal overlap of both endogenic and exogenic processes broad quasi cylindrical mantle upwelling and large bolide impact respectively leading to subduction and spreading two signature processes of modern plate tectonics the spatial and temporal intersection of these processes could have occurred variably across early earth and thus subduction could have begun at different locations and times globally the hypothesis postulates that the ability of a terrestrial planet to evolve plate tectonics results from a balance between the strength of its lithosphere and the size bolide it can attract and survive both factors to a fi rst order are a function of planet size

a tensile flexural model for the initiation of subduction we argue that subduction may be initiated at passive continental margins without shortening the lithosphere overcoming the lithosphere s high compressive strength requires special circumstances and these make it difficult to explain the nearly complete recycling of old sea floor instead we present a model that predicts tensile decoupling of the continental and oceanic lithosphere passive rifting and foundering of the seafloor beneath material welling up in the rift this occurs because the lithosphere in a new ocean basin establishes mechanical continuity with the continent at a depth comparable to mid ocean ridges later subsidence at the margin is therefore inhibited by flexure which implies shear stresses that promote fault slip and tensile stresses necessary to balance the component of the plate s weight directed down the margin slope we show that this tension can more than offset ridge push in our model an important additional tension arises from basal shear tractions resisting the plate s motion away from the mid ocean ridge although these tractions cannot be evaluated with confidence slip on a high angle fault decouples the oceanic and continental lithosphere when shear stresses arising from flexure and the applied tension exceed the lithosphere s shear strength under these loads a passive rift then forms allowing a mantle column to rise to the height of mid ocean ridges over 3 km above the old seafloor and flow onto the surface in a gravity current this load flexes the plate downwards which enhances the flow and lets the old oceanic lithosphere founder this model is consistent with the presence of oceanic material in continental forearcs and the youth of ophiolites when they are obducted as both might be explained by rifting a margin and underthrusting the juvenile crust formed there boninites in ophiolite complexes and tectonites at their bases show depleted and hydrated geochemistries consistent with the melting of rifted mantle lithosphere that receives volatiles from foundering oceanic crust

heat flow maintained in ocean basins of great age investigations in the north equatorial west pacific summary multipenetration heat flow measurements have been made at four sites in deep basins of the west central pacific ocean the west mariana basin central mariana basin nauru basin and central pacific basin the final heat flows are respectively 46 6 f 0 5 49 4 f 0 2 44 2 f 0 9 and 49 5 f 1 1 mw m each site was surveyed by single channel seismic reflection profiling and provided a gravity core the instrument measured thermal conductivity in situ over the entire depth intervals used for determination of the gradients and the reduction scheme iterated conductivity and heat capacity changes into the fitting procedure both of entry frictional decays and of conductivity heat pulse decays the absolute accuracy of the instrument should approach 2 per cent and the first site would make a good intercalibration standard for heat flow measurement the heat flow variation between the sites is real and there is also a significant variation in the isostatically reduced depths of the sites there is no age progression of either depth or heat flow and when five other good multidata points are included the relationship between depth and heat flow conforms to that expected from simple cooling models only in an average sense for the whole group the most plausible explanation for the variations is that heat flow and thermal elevation are dependent on different levels of deep lithosphere reheating at different times between 70 and 120myr ago it is suggested that additional topographic variation is caused by the different accumulations of sediment and lava flows at each site and to errors in the isostatically reduced depths due to incomplete knowledge of the stratigraphy down to the crust mantle interface these explanations of the topographic variation could be tested by seismic refraction measurements

on some consequences and possible causes of layered mantle convection laboratory and numerical experiments are used to study convection confined to discrete layers the laboratory experiments which use the water content of a predominantly glycerine solution as a way of establishing chemical density variations show that convection will remain confined to superimposed layers when the chemical density contrast between layers is greater than the density change associated with the greatest temperature difference within the convecting system the numerical experiments show a double thermal boundary layer at the interface between layers a feature that would result in an average temperature increase of about 500 c in the mantle near 700 km depth if the lower mantle is isolated as suggested by recent geochemical models the possible causes of a layered mantle convection system include major element chemistry variations the combined effect of major element chemistry and a phase change or a phase change alone if it has a sufficiently negative clapeyron slope

evolving force balance during incipient subduction nearly half of all active subduction zones initiated during the cenozoic all subduction zones associated with active back arc extension have initiated since the eocene hinting that back arc extension may be intimately associated with an interval several tens of myr following subduction initiation that such a large proportion of subduction zones are young indicates that subduction initiation is a continuous process in which the net resisting force associated with forming a new subduction zone can be overcome during the normal evolution of plates subduction initiation is known to have occurred in a variety of tectonic settings old fracture zones transform faults and extinct spreading centers and through polarity reversal behind active subduction zones although occurring within different tectonic settings four known subduction initiation events izu bonin mariana ibm along a fracture zone tonga kermadec along an extinct subduction boundary new hebrides within a back arc and puysegur fiordland along a spreading center were typified by rapid uplift within the forearc followed by sudden subsidence other constraints corroborate the compressive nature of ibm and tonga kermadec during initiation using an explicit finite element method within a two dimensional domain we explore the evolving force balance during initiation in which elastic flexure viscous flow plastic failure and heat transport are all considered in order to tie theory with observation known tectonic settings of subduction initiation are used as initial and boundary conditions we systematically explore incipient compression of a homogeneous plate a former spreading center and a fracture zone the force balance is typified by a rapid growth in resisting force as the plate begins bending reaching a maximum value dependent on plate thickness but typically ranging from 2 to 3 1012 n m for cases that become self sustaining this is followed by a drop in stress once a shear zone extends through the plate the formation of a throughgoing fault is associated with rapid uplift on the hanging wall and subsidence on the footwall cumulative convergence not the rate of convergence is the dominant control on the force balance viscous tractions influence the force balance only if the viscosity of the asthenosphere is 1020 pa s and then only after plate failure following plate failure buoyancy of the oceanic crust leads to a linear increase with crustal thickness in the work required to initiate subduction the total work done is also influenced by the rate of lithospheric failure a self sustaining subduction zone does not form from a homogeneous plate a ridge placed under compression localizes subduction initiation but the resisting ridge push force is not nearly as large as the force required to bend the subducting plate the large initial bending resistance can be entirely eliminated in ridge models explaining the propensity for new subduction zones to form through polarity reversals a fracture zone fz placed in compression leads to subduction initiation with rapid extension of the overriding plate a fz must be underthrust by the older plate for 100 150 km before a transition from forced to self sustaining states is reached in fz models the change in force during transition is reflected by a shift from forearc uplift to subsidence subduction initiation is followed by trench retreat and back arc extension moderate resisting forces associated with modeled subduction initiation are consistent with the observed youth of pacific subduction zones the models provide an explanation for the compressive state of western pacific margins before and during subduction initiation including ibm and tonga kermadec in the eocene and the association of active back arcs with young subduction zones on the basis of our dynamic models and the relative poles of rotation between pacific and australia during the eocene we predict that the northern segment of the tonga kermadec convergent margin would have initiated earlier with a progressive southern migration of the transition between forced and self sustaining states

small scale convective instability and upper mantle viscosity under california thermal calculations and convection analysis constrained by seismic tomography results suggest that a small scale convective instability developed in the upper 200 kilometers of the mantle under california after the upwelling and cooling of asthenosphere into the slab window associated with the formation of the san andreas transform boundary the upper bound for the upper mantle viscosity in the slab window 5 x 10 19 pascal seconds is similar to independent estimates for the asthenosphere beneath young oceanic and tectonically active continental regions these model calculations suggest that many tectonically active continental regions characterized by low upper mantle seismic velocities may be affected by time dependent small scale convection that can generate localized areas of uplift and subsidence

large scale atmospheric circulation patterns during polar low events over the nordic seas 1 the large scale atmospheric circulation during polar low events over the nordic seas the north atlantic between greenland and novaya zemlya is investigated on the basis of a recently compiled data set that was derived purely from satellite observations a classification system is applied that divides the polar lows into four types western polar lows eastern polar lows greenland lee polar lows and storm track polar lows type specific large scale atmospheric circulation patterns are identified from national centers for environmental prediction reanalysis data these are distinct in sea level pressure upper level geopotential height and the difference between the skin temperature of the ocean and upper level temperature eastern polar lows are found to be associated with a strong blocking situation caused by anomalously high pressure over iceland and a synoptic scale low pressure anomaly over the barents sea a weaker blocking situation with an anomalous ridge over greenland reaching into the irminger sea and a low pressure anomaly over the norwegian sea favors the development of western polar lows typical values for polar low genesis are a geopotential height of 5030 geopotential meters at 500 hpa and a difference of 48 k between the skin temperature of the ocean and the temperature at 500 hpa the locations of upper level low pressure anomalies relative to the locations of the related anomalies at the sea level show that western and eastern polar lows form on average in a much less baroclinic large scale environment than greenland lee and storm track polar lows

the conditions for plate tectonics on super earths inferences from convection models with damage abstract numerical simulations of mantle convection with a damage grainsize feedback are used to develop scaling laws for predicting conditions at which super earths would have plate tectonics in particular the numerical simulations are used to determine how large a viscosity ratio between pristine lithosphere and mantle l m can be offset by damage to allow mobile plate like convection regime diagrams of l m versus the damage number d show that the transition from stagnant lid to mobile convection occurs for higher l m as d increases a similar trend occurs for increasing rayleigh number we hypothesize a new criterion for the onset of plate tectonics on terrestrial planets that damage must reduce the viscosity of shear zones in the lithosphere to a critical value equivalent to the underlying mantle viscosity a scaling law based on this hypothesis reproduces the numerical results for the earth damage is efficient in the lithosphere and provides a viable mechanism for the operation of plate tectonics we scale our theory to super earths and map out the transition between plate like and stagnant lid convection with a planetary plate tectonic phase diagram in planet size surface temperature space both size and surface conditions are found to be important with plate tectonics being favored for larger cooler planets this gives a natural explanation for earth venus and mars and implies that plate tectonics on exoplanets should correlate with size incident solar radiation and atmospheric composition

conditions for the onset of plate tectonics on terrestrial planets and moons abstract plate tectonics on earth is driven by the subduction and stirring of dense oceanic lithosphere into the underlying mantle for such a regime to exist on any planet stresses associated with mantle convection must exceed the strength of the lithosphere this condition is sufficiently restrictive that plate tectonics currently operates only on earth and mantle convection in most terrestrial planets and moons is probably in a stagnant lid regime convective stresses on the lithosphere depend on the viscosity and velocity of underlying cold downwellings the lithospheric yield stress is controlled by its friction coefficient and elastic thickness the depth to the brittle ductile transition or bdt both convective stresses and the plate s yield strength depend critically on the size thermal state and cooling history of a planet accordingly here we use numerical simulations and scaling theory to identify conditions in which mantle convection leads to lithospheric failure for a range of conditions relevant to the terrestrial planets whereas earth is expected to be in a plate tectonic regime over its full thermal evolution the moon and mercury are expected to have always remained in a stagnant lid regime venus io and europa currently fall on the transition between the two regimes which is consistent with an episodic style of mantle convection for venus a tectonic component to deformation on io and the resurfacing history and lithospheric evolution of europa our results suggest that venus may have been in a plate tectonic regime in the past while stagnant now it is plausible that mars may have also been in an active lid regime depending on whether there was liquid water on the surface

cenozoic plate driving forces past studies of plate driving forces have concluded that the forces due to subducted slabs in the upper mantle and those due to the thickening of the oceanic lithosphere are the principal driving forces we reexamine the balance of driving forces for the present day and extend our analysis through the cenozoic using an analytical torque balance method which accounts for interactions between plates via viscous coupling to the induced mantle flow we use an evolving mantle density heterogeneity field based on the last 200 myr of subduction to drive plate motions an approach which has proven successful in predicting the present day mantle heterogeneity field we find that for plausible upper mantle viscosities the forces due to subducted slabs in the cenozoic and mesozoic account for in excess of 90 of plate driving forces and those due to lithospheric thickening for less than 10

post glacial rebound and asthenosphere viscosity in iceland during the weichselian glaciation iceland was covered with an ice cap which caused downward flexure of the earth s surface the post glacial rebound in iceland was very rapid being completed in about 1 000 years the length of this time interval constrains the maximum value of asthenosphere viscosity in iceland to be 1 times 10 sup 19 pa s or less further clarification of the ice retreat and uplift history may reveal lower viscosity current changes in the mass balance of icelandic glaciers must lead to measurable elevation changes considering this low viscosity expected current elevation changes around the vatnajoekull ice cap are of the order of 1 cm per year due to mass balance change in this century

initiation of subduction by small scale convection 1 the possibility of initiation of subduction by sublithospheric small scale convection is investigated with the help of fowler s theory numerical constraints on the lithospheric stresses generated by convection and direct numerical simulations of subduction the derived scaling laws show that the critical yield stress for initiation of subduction by small scale convection is of the order of 30 mpa this value is comparable with the critical yield stress required for initiation of subduction by plate forces and with the estimates of stresses in the subduction zone this implies that small scale convection can play an important role in the initiation of subduction on the present day earth this also implies that small scale convection can initiate subduction in the absence of already occurring plate tectonics and thus can be a trigger mechanism for plate tectonics on the early earth and other terrestrial planets the critical yield stress for initiation of subduction by either mechanism is at least one order of magnitude smaller than the estimates of the lithospheric strength suggested by laboratory experiments on rock deformation such a low lithospheric strength is most likely caused by a relatively high abundance of water in the interior of the earth since the critical yield stress depends only weakly on planetary parameters water seems to be the key condition required for plate tectonics to occur on a terrestrial planet