geoid and topography over subduction zones the effect of phase transformations 1 the association between local maxima in the geoid and subduction zones is examined while it is well known that subduction zones are associated with broad local maxima in the geoid at spherical harmonic degrees 4 9 there is an impressive correlation between back arc geoid maxima at the 5000 km length scale i e spherical harmonic degrees 9 and subduction zones with the exception of the scotia and caribbean arcs geoid maxima are also observed in the forearc regions of the aleutian and central american subduction zones numerical modeling of compressible convection with multiple phase transformations indicates that an increase in viscosity at or near the 660 km phase transformation is necessary to explain the pattern of geoid maxima associated with subduction zones the addition of an endothermic phase transformation at 660 km depth does not provide sufficient resistance to the slab to eliminate the need for a viscosity increase at or near 660 km short wavelength geoid and topography profiles might provide an additional constraint on the deformation at subduction zones however it may be necessary to incorporate geologic history into subduction calculations to further improve the results presented here

gravity above subduction zones and forces controlling plate motions 1 short and intermediate wavelength gravity and geoid anomalies are used to provide constraints on the mechanical structure of subduction zones and on the forces involved this study is based on 2d cartesian dynamically self consistent models with newtonian or power law rheologies we show that both strong decoupling of the two convergent plates shear stresses of the order of 107 pa and weakened bending lithosphere are necessary to reproduce the observed geoid and gravity data good fits are found for relatively low failure stresses 30 50 mpa for all models providing reasonable predictions of gravity only a small fraction of the downgoing slab weight is transmitted to the surface plates about 10 of the energy is dissipated in the contact zone between the two plates 10 to 20 in the bending region and more than 70 in the sublithospheric mantle the basal tractions on the order of 1 4 1012 n m induce a net motion of the plates with the subducting lithosphere moving faster than predicted by the no net motion principle a marked positive geoid anomaly is predicted above subduction zones at intermediate wavelength 2000 4000 km in the case of pure whole mantle convection such large geoid highs are not observed introducing partial layering i e mass anomalies hampering mantle flow through the transition zone is necessary to reconcile model predictions and observations for these wavelengths

constraints on the lateral strength of slabs from three dimensional dynamic flow models if the viscosity of subducted lithosphere is described purely by temperature and pressure it should remain considerably more viscous than the surrounding material as it descends through the upper mantle many dynamic models of the long wavelength geoid associated with slabs require an increase in viscosity from the upper to the lower mantle but assume the perturbation to the lateral viscosity structure arising from the slabs themselves can be ignored previous studies indicate that strong and localized viscosity variations should have a dramatic influence on the geoid we present 3d finite element models of the regional geoid of the western pacific subduction zones slab buoyancies and viscosities are defined using the distribution of seismicity the geoid is very sensitive to the lateral strength of the slab very viscous slabs penetrating a low viscosity mantle generate significant 10 50 m long wavelength geoid lows opposite to the geoid high which is observed over slabs to obtain a geoid high comparable to that observed the lower mantle viscosity must be 60 200 times greater than the upper mantle viscosity and the slab must be in contact with the lower mantle these strict requirements suggest that slabs have been weakened and cannot act as stress guides from the deep mantle to the surface

slab detachment during continental collision influence of crustal rheology and interaction with lithospheric delamination abstract collision between continents can lead to the subduction of continental material if the crust remains coupled to the downgoing slab a large buoyancy force is generated this force slows down convergence and promotes slab detachment if the crust resists to subduction it may decouple from the downgoing slab and be subjected to buoyant extrusion we employ two dimensional thermo mechanical modelling to study the importance of crustal rheology on the evolution of subduction collision systems we propose simple quantifications of the mechanical decoupling between lithospheric levels and the potential for buoyant extrusion of the crust the modelling results indicate that a variable crustal rheological structure results in slab detachment delamination or the combination of both mechanisms a strong crust provides coupling at the moho low and remains coherent during subduction low it promotes deep subduction of the crust 180 km and slab detachment exhumation occurs in coherent manners via eduction and thrusting slab detachment triggers the development of topography 4 5 km close to the suture a contrasting style of collision occurs using a weak crustal rheology mechanical decoupling at the moho high promotes the extrusion of the crust high disabling slab detachment ongoing shortening leads to buckling of the crust and development of topography on the lower plate collisions involving rheologically layered crust allow decoupling at mid crustal depths this structure favours both the extrusion of upper crust and the subduction of the lower crust such collisions are successively affected by delamination and slab detachment topography develops together with the buoyant extrusion of crust onto the foreland and is further amplified by slab detachment our results suggest that the occurrence of both delamination apennines and slab detachment himalayas in orogens may indicate differences in the initial crustal structure of subducting continental plates in these regions

the owles iop2b lake effect snowstorm dynamics of the secondary circulation abstractintense lake effect snowfall results from a long lake axis parallel llap precipitation band that often forms when the flow is parallel to the long axis of an elongated body of water such as lake ontario the intensity and persistence of the localized precipitation along the downwind shore and farther inland suggests the presence of a secondary circulation that helps organize such a band and maintain it for some time as the circulation is advected inland unique airborne vertical plane dual doppler radar data are used here to document this secondary circulation in a deep well organized llap band observed during intensive observing period iop 2b of the ontario winter lake effect systems owles field campaign the circulation centered on a convective updraft intensified toward the downwind shore and only gradually weakened inland the question arises as to what sustains such a circulation in the vertical plane across the llap band wrf model simulations indicate that the primary llap band a

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

a numerical study of the interaction between the mantle wedge subducting slab and overriding plate abstract we have formulated a numerical model with strongly temperature dependent viscosity to calculate thermal structure and flow field in subduction zones one important particularity of the model is that the overriding plate is not fixed over its whole thickness in order to allow material exchange between the wedge and the upper lithosphere numerical problems due to very high viscosity contrasts are avoided by coupling a finite difference method and a finite element method for solving the energy conservation equation and the stockes equation respectively in this model a temperature decrease from 1400 to 1300 c increases the viscosity by an order of magnitude we study the temperature structure and the velocity field of the subducting slab and mantle wedge surface heat flow velocity anomalies and geometry of the partial melting zone are also calculated to study the effect that boundary conditions play on the interaction between the mantle wedge overriding plate and subducting plate we examine models with both fixed and free slip conditions applied to the overriding plate when the overriding plate is allowed to move laterally free slip the subducting slab is thick and both the temperature field and the convective motions in the mantle wedge are similar to those observed when using constant viscosity numerical models or analytical corner flow models if the surface of the overriding plate is fixed the subducting slab is thin and the mantle wedge impinges upon the overriding plate forming a high temperature nose between the overriding plate and subducting lithosphere furthermore viscous decoupling occurs implicitly at shallow depth between the slab and the wedge because hot material from the wedge is entrained close to the trench in that case the subducting slab tectonically erodes the lower lithosphere of the overriding plate leading to high temperatures low seismic velocities high attenuation and high heat flow beneath volcanic arc in agreement with geophysical observations

the dynamical control of subduction parameters on surface topography the long wavelength surface deflection of earth s outermost rocky shell is mainly controlled by large scale dynamic processes like isostasy or mantle flow the largest topographic amplitudes are therefore observed at plate boundaries due to the presence of large thermal heterogeneities and strong tectonic forces distinct vertical surface deflections are particularly apparent at convergent plate boundaries mostly due to the convergence and asymmetric sinking of the plates having a mantle convection model with a free surface that is able to reproduce both realistic single sided subduction and long wavelength surface topography self consistently we are now able to better investigate this interaction we separate the topographic signal into distinct features and quantify the individual topographic contribution of several controlling subduction parameters results are diagnosed by splitting the topographic signal into isostatic and residual components and by considering various physical aspects like viscous dissipation during plate bending performing several systematic suites of experiments we are then able to quantify the topographic impact of the buoyancy rheology and geometry of the subduction zone system to each and every topographic feature at a subduction zone and to provide corresponding scaling laws we identify slab dip and slightly less importantly slab buoyancy as the major agents controlling surface topography at subduction zones on earth only the island arc high and the back arc depression extent are mainly controlled by plate strength overall his modeling study sets the basis to better constrain deep seated mantle structures and their physical properties via the observed surface topography on present day earth and back through time

subduction kinematics and dynamic constraints the kinematics of subduction zones shows a variety of settings that can provide clues for dynamic understandings two reference frames are used here to describe the simple 2d kinematics of subduction zones in the first the upper plate is assumed fixed whereas in the second frame upper and lower plates move relative to the mantle relative to a fixed point in the upper plate u the transient subduction hinge h can converge diverge or be stationary similarly the lower plate l can converge diverge or be stationary the subduction rate vs is given by the velocity of the hinge h minus the velocity of the lower plate l vs vh vl the subduction rate 1 increases when h diverges and 2 decreases when h converges combining the different movements at least 14 kinematic settings can be distinguished along the subduction zones variable settings can coexist even along a single subduction zone as shown for the 5 different cases occurring along the apennines subduction zone apart from few exceptions the subduction hinge converges toward the upper plate more frequently along e or ne directed subduction zone whereas it mainly diverges from the upper plate along w directed subduction zones accompanying backarc extension before collision orogen growth occurs mostly at the expenses of the upper plate shortening along e ne directed subduction zones whereas the accretionary prism of w directed subduction zones increases at the expenses of the shallow layers of the lower plate the convergence shortening ratio is 1 along e or ne directed subduction zones whereas it is 1 along accretionary prisms of w directed subduction zones backarc spreading forms in two settings along the w directed subduction zones it is determined by the hinge divergence relative to the upper plate minus the volume of the accretionary prism or in case of scarce or no accretion minus the volume of the asthenospheric intrusion at the subduction hinge since the volume of the accretionary prism is proportional to the depth of the decollement plane the backarc rifting is inversely proportional to the depth of the decollement on the other hand along e or ne directed subduction zones few backarc basins form e g aegean andaman and can be explained by the velocity gradient within the hangingwall lithosphere separated into two plates when referring to the mantle the kinematics of subduction zones can be computed either in the deep or in the shallow hotspot reference frames the subduction hinge is mostly stationary being the slab anchored to the mantle along w directed subduction zones whereas it moves w or sw ward along e or ne directed subduction zones surprisingly along e or ne directed subduction zones the slab moves out of the mantle i e the slab slips relative to the mantle opposite to the subduction direction kinematically this subduction occurs because the upper plate overrides the lower plate pushing it down into the mantle as an example the hellenic slab moves out relative to the mantle i e sw ward opposite to its subduction direction both in the deep and shallow hotspot reference frames in the shallow hotspot reference frame upper and lower plates move westward relative to the mantle along all subduction zones this kinematic observation casts serious doubts on the slab negative buoyancy as the primary driving mechanism of subduction and plate motions w directed subduction zones rather provide about 2 3 times larger volumes of lithosphere re entering into the mantle and the slab is pushed down this opposite behavior is consistent with the down dip extension seismicity along e ne directed subduction zones and the frequent down dip compression along the w directed subduction zones subduction kinematics shows that plate velocity is not dictated by the rate of subduction along the w directed subduction zones the rate of subduction is rather controlled i by the hinge migration due to the slab interaction with the easterly trending horizontal mantle wind along the global tectonic mainstream ii by the far field plate velocities and iii by the value of negative buoyancy of the slab relative to the country mantle alternatively e ne nne directed subduction zones have rates of sinking chiefly determined i by the far field velocity of plates and ii by the value of negative buoyancy of the slab relative to the country mantle along this type of subduction the subduction hinge generally advances e ne ward toward the upper plate decreasing the subduction rate but it moves w sw ward relative to the mantle all this indicates that subduction zones have different origin as a function of their geographic polarity and the subduction process is more a passive feature rather than being the driving mechanism of plate motions a rotational component combined with mantle density and viscosity anisotropies seems more plausible for generating the global tuning in the asymmetry of subduction zones

steady plumes in viscously stratified vigorously convecting three dimensional numerical mantle convection models with mobile plates 1 we present results from three dimensional 3 d cartesian geometry computations featuring vigorous high rayleigh number mantle convection in a system incorporating stiff tectonic plates with dynamically determined time dependent velocities we track plume related hot spot motion while calculating associated synthetic hot spot tracks allowing us to estimate plume migration rates relative to plate motion rates plate like surface motion is achieved explicitly by specifying both the plate geometry and rigidity however plate velocities evolve dynamically in response to the buoyancy distribution within the convecting system including buoyancy within the plate itself we find that convection is characterized by downwelling sheets and varying numbers of 3 d upwelling structures when the lower mantle to upper mantle viscosity ratio is varied from 9 to 90 as the lower mantle viscosity is increased relative to the upper mantle viscosity the upwellings in the calculations evolve into vigorous active plumes characterized by long narrow thermal conduits with broad disk like heads the total number shape and persistence of the plumes are affected by the specified viscosity stratification in the calculations our findings indicate that hot spots associated with mantle plumes drift by smaller distances in comparable amounts of time when the lower mantle viscosity is increased relative to the upper mantle viscosity in cases with a high lower mantle to upper mantle viscosity ratio the depth of the flow aligned with plate motion is diminished and the plate scale return flow is integrated into a deep unorganized sluggish layer as a result plume conduit morphology in the lower mantle of such calculations is less influenced by plate motion than in calculations with a small difference between upper mantle and lower mantle viscosity for lower mantle viscosities 30 times and 90 times greater than the upper mantle viscosity the lower mantle velocity field is dominated locally by active plumes anchored in the lower thermal boundary layer of the convecting system and plume migration rates are typically only about 10 of the overriding plate velocity in such cases we also find that plume related hot spots persist for periods that we estimate to be between 0 55 and 1 1 plume transit times

subducted slabs and the geoid 1 numerical experiments with temperature dependent viscosity one of the most powerful constraints on mantle viscosity comes from the correlation of the long wavelength degree 4 9 geoid with that predicted by a density model for subducted slabs the effect of lateral variations of viscosity should be most pronounced at subduction zones due to the strong effect of temperature on viscosity an idealized slab model with temperature dependent viscosity is considered with various lateral and vertical viscosity structures using a two dimensional finite element formulation the viscosity parameterization affects the amplitude of the long wavelength geoid anomaly but not the sign of the correlation between the geoid and density anomalies depth dependent viscosity models with a high viscosity lithospheric layer do not completely match the temperature dependent laterally varying results suggesting that the rheology of the slab does have an effect on the long wavelength surface topography and geoid however this affect is minor suggesting that the radial models of mantle viscosity inferred from surface observables are correct to first order in contrast the short wavelength features are dramatically affected by the rheology of the region surrounding the downwelling further study of the shorter wavelength geoid over subduction zones may provide better insight into subduction zone processes

overriding plate velocity control on surface topography in 2 d models of subduction zones we study the dynamically induced flexural topography in subduction numerical mechanical models we focus on the topographic changes at the overriding plate op surface induced by variations in op kinematics particularly when the subducting plate sp has a stationary motion after having reached the rigid base of the upper mantle our models consist of two viscoelastic plates with free surfaces and an isoviscous mantle friction is imposed along the planar subduction interface we first characterize the main topographic features at a constant op velocity using spatial definitions based on geometrical estimations of the volcanic arc position the models exhibit the formation of a bulge in the forearc area followed landwards by a depression and a smaller second bulge both bracketing the arc region the steady state distance to the trench of these features increases with op velocity their amplitude is affected by the far field op tectonic regime that depends on kinematics and plates and subduction interface strength we next test the effect of sudden changes in op velocity an op acceleration yields a transient topographic tilt during which the outer forearc quickly subsides whereas the arc region uplifts and that is followed by reverse slower motions an op slowdown induces opposite motions the rates of elevation change during the tilt are approximately proportional to velocity variations and mainly sensitive to the sp strength the rates are higher than 0 1 mm yr for velocity changes higher than 1 cm yr we suggest that topographic accommodations of op velocity changes should be considered when quantifying nonisostatic topography

dynamic topography gravity and the role of lateral viscosity variations from inversion of global mantle flow lateral viscosity variations lvvs in the mantle influence geodynamic processes and their surface expressions with the observed long wavelength geoid free air anomaly gravity gradient in three directions and discrete high accuracy residual topography we invert for depth and temperature dependent and tectonically regionalized mantle viscosity with a mantle flow model the inversions suggest that long wavelength gravitational and topographic signals are mainly controlled by the radial viscosity profile the pre cambrian lithosphere viscosity is slightly one order of magnitude higher than that of oceanic and phanerozoic lithosphere plate margins are substantially weaker than plate interiors and viscosity has only a weak apparent dependence on temperature suggesting either a balancing between factors or a smoothing of actual higher amplitude but short wavelength lvvs the predicted large scale lithospheric stress regime compression or extension is consistent with the world stress map thrust or normal faulting both recent compiled high accuracy residual topography and the predicted dynamic topography yield 1 km amplitude long wavelength dynamic topography inconsistent with recent studies suggesting amplitudes of 100 to 500 m such studies use a constant positive admittance transfer function between topography and gravity in contrast to the evidence which shows that the earth has a spatially and wavelength dependent admittance with large negative admittances between 4000 and 10 4 km wavelengths

controls of the structure of subducted slabs numerical simulations of subducting slabs are formulated in which the shape and dip of the slab are determined by the dynamics of the flow rather than imposed a priori the dip of slabs is a function of the time since the initiation of subduction slabs fold develop a kink in dip and thicken on entry into a high viscosity lower mantle comparison of the simulations with seismic observations suggest that the lower mantle is at least 10 30 times more viscous than the upper mantle

stress and flow beneath island arcs summary the flow pattern stress distribution topography and gravity anomalies were computed from numerical models having density and viscosity distributions resemblant to the aleutian arc the results were compatible with the hypothesis that the excess density of the slab drives its descent and that hydrodynamic forces are responsible for topographic and gravity highs over the outer rise seaward of the trench and the frontal arc and lows over the trench in models with simple distributions of rheological parameters the force from the slab was transmitted directly upward producing a negative gravity anomaly over the arc material with low resistance to flow was needed along the fault plane above the slab or within the crust of the frontal arc and within the wedge of asthenosphere above the slab to reduce that force and to allow the horizontal lithosphere to move with the slab models with the resistance to flow thus reduced had outer rises deep trenches horizontal tension seaward of the trench horizontal compression under the trench and downdip tension in the slab free air gravity anomalies which are the sum of between deflections of the free surface due to hydrodynamic forces and direct attractions from the masses driving the flow were not fit excellently by any of the models in part because the coarse grid used precluded accurate representation of the fault zone above the slab and the frontal arc an alternate to the hypothesis that about 5 kb of stress on the fault plane is needed to produce an outer rise is offered by these models shear stress between the slab and the island arc was always below 700 bars in the more successful models if the density distribution was scaled to match the topography of the trench this is much less than the 2000 bars stresses needed if frictional heating causes island arc volcanism

large scale deformation associated with ridge subduction summary continuum models are used to investigate the large scale deformation associated with the subduction of aseismic ridges formulated in the horizontal plane using thin viscous sheet theory these models measure the horizontal transmission of stress through the arc lithosphere accompanying ridge subduction modelling was used to compare the tonga arc and louisville ridge collision with the new hebrides arc and d entrecasteaux ridge collision which have disparate arc ridge intersection speeds but otherwise similar characteristics models of both systems indicate that diffuse deformation low values of the effective stress strain exponent n are required to explain the observed deformation deformation is somewhat insensitive to the vertically integrated strength of the arc inversely proportional to the argand number ar but indicates that the arc lithosphere is not extremely weak ar 100 low values of both ar and n suggest that the thermal structure is typical of cold or normal arcs and that deformation is dominated by flow in the lower crust and mantle in addition low values of n approaching newtonian flow may indicate that specific deformation mechanisms dictate deformation of the arc lithosphere possible mechanisms include low stress grain size dependent creep pyroxenite controlled rheology and mechanisms associated with water weakening changes in the boundary conditions greatly affect deformation within island arcs high rates of arc ridge intersection speed tonga louisville system yield arc parallel tension and crustal thickening in the wake of ridge subduction in contrast low rates of arc ridge intersection speed new hebrides d entrecasteaux system yield compressional deformation directly arcward of the collision zone and transverse strike slip faulting adjacent to the region of compressional deformation localized regions of extensional deformation along the frontal part of the arc adjacent to the collision zone may contribute to the formation of re entrants

viscous flow model of a subduction zone with a faulted lithosphere long and short wavelength topography gravity and geoid the essential features of observed topography free air gravity and the geoid over subduction zones have been reproduced with a dynamically self consistent viscous flow model of a subduction zone with a faulted lithosphere a fault between the overriding and subducting plates is crucial to understanding lithospheric deformation and gravity in a subduction zone the topography for models with a fault is characterized by a deep and narrow trench and a broad and shallow back arc depression quite different from viscous models without faults which only have a broad and large depression over subducting slabs for a model in which a slab is 600 km long and 80 myr old the resulting trench is approximately 3 km deep and 100 km wide and the back arc basin is about 1 5 km deep and 600 km wide the model is characterized by a long wavelength geoid high of about 30 m over the slab a local geoid low of 10 m over the trench and a 150 mgal gravity low over the trench the fact that trench depth and width are reproduced in models with a small resisting shear stress 50 bars suggests that convergent plate boundaries are probably weak

far field tilting of laurentia during the ordovician and constraints on the evolution of a slab under an ancient continent during a brief period of 10 to 15 million years in the middle ordovician the michigan basin departed from its bull s eye subsidence pattern and tilted toward the east opening to the appalachian basin this tilting is observed in maps of tectonic subsidence estimated for the black river and trenton formations and extends over 300 km across the michigan basin and into eastern wisconsin contours of constant tectonic subsidence rate are approximately parallel to the inferred position of the laurentian iapetus convergent margin the distance between the inferred position of the subduction zone to the limit of tilting is approximately 1000 km three alternative models for the tilting are tested two relying on the rigidity of the continental lithosphere and a third on the viscous flow generated by a subducted slab in the first elastic model we assume the edge of the elastic plate is simply loaded from above by a fold and thrust sheet volcanic pile etc this model however cannot simultaneously satisfy the space of tectonic subsidence and subsidence rate even for lithospheres which have a rigidity of 10 28 nm in the second elastic model the laurentian continental margin descends into a trench of an eastward dipping slab that is a slab descending under the iapetus ocean this process cannot generate any significant far field displacements even for extremely rigid plates and must be rejected for the third model we use finite element solutions of a negatively buoyant slab in a viscous medium with a faulted lithosphere such slabs can easily generate not only realistic trenches on the under thrusting plate but also significantly tilt the lithosphere as much as 1000 km from the plate margin the magnitude and distribution of far field displacements depend on the age length and dip angle of the slab in contrast to the elastic models penetration of a west dipping slab beneath the continent can reproduce both the extent magnitude and rate of tectonic subsidence observed in the trenton and black river formations the observed data are best fit by an old slab 140 ma which initially descended at a moderate dip 20 30 for 10 15 m y which then steepened as the slab penetrated deeper into the mantle at the position of the previous ordovician plate margin there are narrow block faulted basins which underwent rapid subsidence 10 m y before the michigan basin tilted toward the east we propose that these earlier subsidence events were caused by the initial descent of the slab under the preexisting cambro ordovician passive margin the time lag of 10 m y may be due to the time it takes the slab to penetrate the upper mantle this result is important for understanding the time evolution of mantle convection and mechanisms for the initiation of subduction

lake effect mode and precipitation enhancement over the tug hill plateau during owles iop2b abstractimproved understanding of the influence of orography on lake effect storms is crucial for weather forecasting in many lake effect regions the tug hill plateau of northern new york hereafter tug hill rising 500 m above eastern lake ontario experiences some of the most intense snowstorms in the world herein the authors investigate the enhancement of lake effect snowfall over tug hill during iop2b of the ontario winter lake effect systems owles field campaign during the 24 h study period total liquid precipitation equivalent along the axis of maximum precipitation increased from 33 5 mm at a lowland 145 m msl site to 62 5 mm at an upland 385 m msl site the latter yielding 101 5 cm of snow however the ratio of upland to lowland precipitation or orographic ratio varied with the mode of lake effect precipitation strongly organized long lake axis parallel bands some of which formed in association with the approach or passage of upper level short wave troughs produced the highest pre

subducted slabs and lateral viscosity variations effects on the long wavelength geoid summary the characteristic broad local maxima exhibited by the long wavelength geoid over subduction zones are investigated with a numerical model of mantle flow in a spherical axisymmetric geometry a synthetic model of buoyancy driven subduction is used to test the effects on the geoid caused by the depth of penetration of the lithosphere into the mantle by the viscosity stratification and by lateral viscosity variations lvv in the lithosphere upper and lower mantle the presence of anomalous slab density in the lower mantle guarantees geoid amplitudes comparable with the observations favouring the picture of slabs that penetrate the transition zone and sink into the deep mantle the viscosity of the lower mantle controls the long wavelength geoid to the first order ensuring a clear positive signal when it is at least 30 times greater than the upper mantle viscosity the presence of lvv in the lithosphere in the form of weak plate margins helps to increase the contribution of the surface topography causing a pronounced reduction of the geoid localized lvv associated with the cold slab play a secondary role if they are in the upper mantle on the other hand highly viscous slabs in the lower mantle exert a large influence on the geoid they cause its amplitude to increase dramatically way beyond the values typically observed over subduction zones long wavelength flow becomes less vigorous as the slab viscosity increases deformation in the upper mantle becomes more localized and power is transferred to short wavelengths causing the long wavelength surface topography to diminish and the total geoid to increase slabs may be then weakened in the lower mantle or retain their high viscosity while other mechanisms act to lower the geoid it is shown that a phase change from perovskite to post perovskite above the core mantle boundary can cause the geoid to reduce significantly thereby helping to reconcile models and observations

surface erosion events controlled the evolution of plate tectonics on earth plate tectonics is among the most important geological processes on earth but its emergence and evolution remain unclear here we extrapolate models of present day plate tectonics to the past and propose that since about three billion years ago the rise of continents and the accumulation of sediments at continental edges and in trenches has provided lubrication for the stabilization of subduction and has been crucial in the development of plate tectonics on earth we conclude that the two largest surface erosion and subduction lubrication events occurred after the palaeoproterozoic huronian global glaciations 2 45 to 2 2 billion years ago leading to the formation of the columbia supercontinent and after the neoproterozoic snowball earth glaciations 0 75 to 0 63 billion years ago the snowball earth event followed the boring billion a period of reduced plate tectonic activity about 1 75 to 0 75 billion years ago that was probably caused by a shortfall of sediments in trenches and it kick started the modern episode of active plate tectonics the rise of continents and the accumulation of sediments in trenches since about three billion years ago has had a crucial role in the emergence and evolution of plate tectonics on earth

a comparison of numerical surface topography calculations in geodynamic modelling an evaluation of the sticky air method summary calculating surface topography in geodynamic models is a common numerical problem besides other approaches the so called sticky air approach has gained interest as a free surface proxy at the top boundary the often used free slip condition is thereby vertically extended by introducing a low density low viscosityfluid layer this allows the air crust interface to behave in a similar manner to a true free surface we present here a theoretical analysis that provides the physical conditions under which the sticky air approach is a valid approximation of a true free surface two cases are evaluated that characterize the evolution of topography on different timescales 1 isostatic relaxation of a cosine perturbation and 2 topography changes above a rising plume we quantitatively compare topographies calculated by six different numerical codes using finite difference and finite element techniques using three different topography calculation methods i direct calculation of topography from normal stress ii body fitting methods allowing for meshing the topography and iii lagrangian tracking of the topography on an eulerian grid it is found that the sticky air approach works well as long as the term st ch hst l 3 is sufficiently small where st and hst are the viscosity and thickness of the sticky air layer and ch and l are the characteristic viscosity and length scale of the model respectively spurious lateral fluctuations of topography as observed in some marker based sticky air approaches may effectively be damped by an anisotropic distribution of markers with a higher number of markers per element in the vertical than in the horizontal direction

the midwest snow storm of 8 11 december 1977 abstract conventional surface and upper air aerological data are combined with radar and aircraft measurements to give a description of a major winter storm that deposited over 69 cm of new snow at muskegon michigan between 8 and 11 december 1977 it is shown that most of this snow occurred during four distinct episodes three of which were related to air water temperature contrasts as arctic air flowed over the relatively warm surface of lake michigan these four episodes have been identified as 1 pre cold frontal 2 post cold frontal 3 secondary trough and 4 mesolow phases during the pre cold front phase heavy snow was associated with an advancing synoptic scale trough little direct effect of the lakes was detected following passage of the cold front strong northwesterly winds across lake michigan resulted in strong air mass transformation with light steady snow along the downwind shoreline the evidence suggests the presence of horizontal roll convection as the dominant organizational mode

snowbands over the english channel and irish sea during cold air outbreaks persistent northerly to easterly cold air outbreaks affected the uk during the winters of 2009 10 and 2010 11 with the resulting convection frequently organizing into snowbands over the english channel and irish sea sounding data and composite radar reflectivity images from the met office nimrod precipitation radar network reveal that these bands formed along the major axis of each body of water or sea when the boundary layer flow was roughly parallel to each of those axes along channel for both seas a band was present the majority of times that the 850 hpa flow was along channel of these times of along channel flow the 850 hpa wind speed and surface to 850 hpa temperature difference were significantly greater when bands were present than when they were not for the english channel only the land sea temperature difference was also significantly greater when bands were present than when they were not in a real data weather research and forecasting model wrf control simulation of a typical band over the english channel a trough develops over the water and offshore air streams from either side converge along it in the absence of surface fluxes the trough convergence and organized precipitation fail to develop altogether orography and roughness length variations are less important in band development affecting only the location and morphology

climatology of lake effect snowstorms of the great salt lake abstract characteristics of lake effect snowstorms associated with the great salt lake are described using wsr 88d radar imagery 16 well defined and 18 marginal lake effect events were identified from september 1994 through may 1998 excluding june august with the former used for more detailed analysis precipitation during the well defined events was frequently characterized by the irregular development of radar echoes over and downstream of the great salt lake the most commonly observed precipitation structures were solitary wind parallel bands that developed along or near the major axis of the gsl and broad area precipitation shields with embedded convective elements that formed near the southern shoreline regional scale composite analyses and rawinsonde derived statistics showed that the lake effect events occurred in post frontal westerly to northerly 700 hpa flow following the passage of an upper level trough and associated low level cold front the lake effect environment was characterized by

a simple semi dynamic model of the subduction zone effects of a moving plate boundary on the small scale convection under the island arc summary a simple semi dynamic numerical model of the subduction zone which incorporates the plate velocity and a moving plate boundary as a priori constraints has been developed this model may be useful for studying specific subduction zones where a history of kinematic movements of plates is well constrained as an application of this model the effects of a moving plate boundary on the nature of the small scale convection in the mantle wedge were studied under the restriction of given relative velocity between two converging plates the absolute velocity of each plate was changed and i found that the characteristics of small scale convection are not significantly affected by the movement of the plate boundary at least within the parameter range used in this study i also found a small systematic change which may be related to the difference of the heat supply to the mantle wedge caused by the moving plate boundary thus the model with a fixed continental plate and moving oceanic plate with a given relative velocity may be sufficient for studying the present subject this model can be extended into a simple 3d subduction zone model

influence of lithosphere asthenosphere interaction on the overriding lithosphere in a subduction zone numerical modeling 1 several factors influence the topography of the overriding plate in a subduction zone mantle wedge flow caused by motion of the descending slab lithosphere elasticity and interaction between the solid lithosphere and the viscous asthenosphere this paper presents the results of new subduction modeling that incorporates both elastic elastoplastic deformation in the overriding plate and viscous deformation in the asthenosphere the effects of fluid solid coupling and interaction were considered using an arbitrary lagrangian eulerian technique that solves for the continuous deformation of fluid flow in contact with the rigid plate modeling results show that thicker lithosphere which is harder to deform contributes to the development of broadly shallower back arc basins as lithospheric strength decreases its ability to withstand outside force fluid load resulting from flow in the mantle wedge also declines leading to the development of deep and narrow basins large scale deformation of the lithosphere modifies both the mantle wedge geometry and the wedge flow pattern which in turn increases the fluid load and thereby enhances topographic variation generally the effect of interaction is enhanced with lithospheric deformation the elastoplastic model also leads to a weaker lithosphere and a correspondingly larger deflection in the lithosphere similarly to previous studies back arc basins are deeper when mantle viscosity is higher for a 50 km thick elastic lithosphere the lithosphere asthenosphere interaction is negligible when the asthenospheric viscosity is less than 4 1019 pa s increasing the rate of subduction or decreasing the dip angle will also lead to the development of deeper basins

holocene climatic change and the development of the lake effect snowbelt in michigan usa lake effect snow is an important constraint on ecological and socio economic systems near the north american great lakes little is known about the holocene history of lake effect snowbelts and it is difficult to decipher how lake effect snowfall abundance affected ecosystem development we conducted oxygen isotope analysis of calcite in lake sediment cores from northern lower michigan to infer holocene climatic variation and assess snowbelt development the two lakes experience the same synoptic scale climatic systems but only one of them huffman lake receives a significant amount of lake effect snow a 177 cm difference in annual snowfall causes groundwater inflow at huffman lake to be 18o depleted by 2 3 relative to o brien lake to assess when the lake effect snowbelt became established we compared calcite 18o profiles of the last 11 500 years from these two sites the chronologies are based on accelerator mass spectrometry 14c ages of 11 and 17 terrestrial plant samples from huffman and o brien lakes respectively the values of 18o are low at both sites from 11 500 to 9500 cal yr bp when the laurentide ice sheet lis exerted a dominant control over the regional climate and provided periodic pulses of meltwater to the great lakes basin carbonate 18o increases by 2 6 at o brien lake and by 1 4 at huffman lake between 9500 and 7000 cal yr bp suggesting a regional decline in the proportion of runoff derived from winter precipitation the great lakes snowbelt probably developed between 9500 and 5500 cal yr bp as inferred from the progressive 18o depletion at huffman lake relative to o brien lake with the largest increase of lake effect snow around 7000 cal yr bp lake effect snow became possible at this time because of increasing contact between the great lakes and frigid arctic air these changes resulted from enhanced westerly flow over the great lakes as the lis collapsed and from rapidly rising great lakes levels during the nipissing transgression the 18o difference between huffman and o brien lakes declines after 5500 cal yr bp probably because of a northward shift of the polar vortex that brought increasing winter precipitation to the entire region however 18o remains depleted at huffman lake relative to o brien lake because of the continued production of lake effect snow

subduction zones observations and geodynamic models abstract this review of subduction and geodynamic models is organized around three central questions 1 why is subduction asymmetric 2 are subducted slabs strong or weak 3 how do subducted slabs interact with phase transformations changes in mantle rheology and possibly chemical boundaries in the mantle based on laboratory measurements of the temperature dependence of olivine one would conclude that the core of a subducting slab is at least 10 000 times more viscous than ambient mantle however there are a number of complementary but independent observations that suggest that slabs are much weaker than this slabs undergo significant deformation in the upper mantle and may thicken to twice their original width by the time they reach the base of the transition zone the lack of a clear correlation between the observed dip angle of deep slabs and plate velocity rate of trench migration and slab age in modern subduction zones is consistent with hypothesis that subduction is a time dependent phenomenon both tank and numerical convection experiments with plates conclude that subduction is not a steady phenomenon but that slabs bend thicken stretch and change dip through time this is at odds with the assumptions used in steady state slab thermal models where slab deformation is not considered

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

vertical motions of the puerto rico trench and puerto rico and their cause 1 the puerto rico trench exhibits great water depth an extremely low gravity anomaly and a tilted carbonate platform between reconstructed elevations of 1300 m and 4000 m i argue that these features are manifestations of large vertical movements of a segment of the puerto rico trench its forearc and the island of puerto rico that took place 3 3 m y ago over a time period as short as 14 40 kyr i explain these vertical movements by a sudden increase in the slab s descent angle that caused the trench to subside and the island to rise the increased dip could have been caused by shearing or even by a complete tear of the descending north american slab although the exact nature of this deformation is unknown the rapid 14 40 kyr and uniform tilt along a 250 km long section of the trench is compatible with scales of mantle flow and plate bending the proposed shear zone or tear is inferred from seismic morphological and gravity observations to start at the trench at 64 5 w and trend southwestwardly toward eastern puerto rico the tensile stresses necessary to deform or tear the slab could have been generated by increased curvature of the trench following a counterclockwise rotation of the upper plate and by the subduction of a large seamount

porphyry copper deposit formation by sub volcanic sulphur dioxide flux and chemisorption the processes that create economic grade accumulations of metals above magma chambers are unclear high temperature laboratory experiments show that rapid reactions between magmatic gases and earth s crust can trigger efficient metal deposition

a low viscosity wedge in subduction zones geochemical petrologic and seismological observations indicate that there may be high concentrations of water in the region above a subducting slab the mantle wedge which could decrease the viscosity of the mantle locally by several orders of magnitude using numerical models we demonstrate that a low viscosity wedge has a dramatic influence on the force balance in a subduction zone and leads to an observable signal in the topography gravity and geoid a regional dynamic model of the tonga kermadec subduction zone shows that the viscosity of the wedge is at least a factor of 10 smaller than surrounding mantle lithosphere and asthenosphere consistent with estimates from seismic dissipation and deformation experiments

crustal structure of the kermadec arc from mango seismic refraction profiles three active source seismic refraction profiles are integrated with morphological and potential field data to place the first regional constraints on the structure of the kermadec subduction zone these observations are used to test contrasting tectonic models for an along strike transition in margin structure previously known as the 32 s boundary we use residual bathymetry to constrain the geometry of this boundary and propose the name central kermadec discontinuity ckd north of the ckd the buried tonga ridge occupies the forearc with vp 6 5 7 3 km s 1 and residual free air gravity anomalies constrain its latitudinal extent north of 30 5 s width 110 20 km and strike 005 south of 25 s south of the ckd the forearc is structurally homogeneous down dip with vp 5 7 7 3 km s 1 in the havre trough backarc crustal thickness south of the ckd is 8 9 km which is up to 4 km thinner than the northern havre trough and at least 1 km thinner than the southern havre trough we suggest that the eocene arc did not extend along the current length of the tonga kermadec trench the eocene arc was originally connected to the three kings ridge and the ckd was likely formed during separation and easterly translation of an eocene arc substrate during the early oligocene we suggest that the first order crustal thickness variations along the kermadec arc were inherited from before the neogene and reflect mesozoic crustal structure the cenozoic evolution of the tonga kermadec hikurangi margin and along strike variations in the duration of arc volcanism

the stress heat flow paradox and thermal results from cajon pass conventional friction models predict a substantial thermal anomaly associated with active traces of strike slip faults but no such anomaly is observed from over 100 heat flow determinations along 1 000 km of the san andreas fault the cajon pass well is being drilled to bring deep heat flow and stress data to bear on this paradox preliminary heat flows to 1 7 km in the neighboring arkoma well indicated a positive anomaly of about 30 as would be expected if frictional resistance of the fault were 50 mpa or so but the pliocene and pleistocene history of the transverse ranges particularly high rates of uplift and erosion could also explain the anomaly preliminary stress results from cajon pass and a new interpretation of regional data by mark d zoback and colleagues suggests that the maximum compressive stress near the fault is almost normal to the trace and hence the resolved shear stress is low and the fault weak although this is consistent with long standing evidence from low heat flow vexing questions remain regarding the physics of strike slip faulting the heat flow data show large variability with depth probably from three dimensional structure and an overall decrease from over 90 mw m in the upper kilometer to less than 80 mw m in the lower 300 m with no evidence of advective heat transfer whether this decreasing trend is significant and whether it continues downward to the background value 70 mw m less the radiogenic heat from overlying rocks as expected for a weak fault must be verified with data from greater depth

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

free surface formulation of mantle convection ii implication for subduction zone observables summary viscous and viscoelastic models for a subduction zone with a faulted lithosphere and internal buoyancy can self consistently and simultaneously predict long wavelength geoid highs over slabs short wavelength gravity lows over trenches trench forebulge morphology and explain the high apparent strength of oceanic lithosphere in trench environments the models use two different free surface formulations of buoyancydriven flows see for example part i lagrangian viscoelastic and pseudo free surface viscous formulations the lower mantle must be stronger than the upper in order to obtain geoid highs at long wavelengths trenches are a simple consequence of the negative buoyancy of slabs and a large thrust fault decoupling the overriding from underthrusting plates the lower oceanic lithosphere must have a viscosity of less than pa s in order to be consistent with the flexural wavelength of forebulges forebulges are dynamically maintained by viscous flow in the lower lithosphere and mantle and give rise to apparently stiffer oceanic lithosphere at trenches with purely viscous models using a pseudo free surface formulation we find that viscous relaxation of oceanic lithosphere in the presence of rapid trench rollback leads to wider and shallower back arc basins when compared to cases without viscous relaxation moreover in agreement with earlier studies the stresses necessary to generate forebulges are small 100 bars compared to the unrealistically high stresses needed in classic thin elastic plate models

the role of cloud top generating cells and boundary layer circulations in the finescale radar structure of a winter cyclone over the great lakes abstractdata from airborne w band radar are used in conjunction with thermodynamic fields from the weather research and forecasting model and air parcel back trajectories from the hysplit model to investigate the finescale reflectivity vertical motion and airmass structure of the comma head of a winter cyclone in the vicinity of the great lakes cloud top generating cells formed along an upper level frontal boundary vertically separating dry air which 48 h earlier was located in the upper troposphere over south central canada from moist air which was located in the lower troposphere over the southeast united states the stronger updrafts within the generating cells had vertical velocities ranging from 1 to 3 m s 1 the generating cells were important to precipitation production within the comma head precipitation trails formed within the generating cells could sometimes be followed to the boundary layer before merging boundary layer air beneath the cyclone s comma head exhibited convective circulati

plate tectonics on the early earth limitations imposed by strength and buoyancy of subducted lithosphere the tectonic style and viability of modern plate tectonics in the early earth is still debated field observations and theoretical arguments both in favor and against the uniformitarian view of plate tectonics back until the archean continue to accumulate here we present the first numerical modeling results that address for a hotter earth the viability of subduction one of the main requirements for plate tectonics a hotter mantle has mainly two effects 1 viscosity is lower and 2 more melt is produced which in a plate tectonic setting will lead to a thicker oceanic crust and harzburgite layer although compositional buoyancy resulting from these thick crust and harzburgite might be a serious limitation for subduction initiation our modeling results show that eclogitization significantly relaxes this limitation for a developed ongoing subduction process furthermore the lower viscosity leads to more frequent slab breakoff and sometimes to crustal separation from the mantle lithosphere unlike earlier propositions not compositional buoyancy considerations but this lithospheric weakness could be the principle limitation to the viability of plate tectonics in a hotter earth these results suggest a new explanation for the absence of ultrahigh pressure metamorphism uhpm and blueschists in most of the precambrian early slabs were not too buoyant but too weak to provide a mechanism for uhpm and exhumation

large trench parallel gravity variations predict seismogenic behavior in subduction zones we demonstrate that great earthquakes occur predominantly in regions with a strongly negative trench parallel gravity anomaly tpga whereas regions with strongly positive tpga are relatively aseismic these observations suggest that over time scales up to at least 1 million years spatial variations of seismogenic behavior within a given subduction zone are stationary and linked to the geological structure of the fore arc the correlations we observe are consistent with a model in which spatial variations in frictional properties on the plate interface control trench parellel variations in fore arc topography gravity and seismogenic behavior

modeling the dynamics of subducting slabs cold dense subducting lithosphere provides the primary force driving tectonic plates at earth s surface the force available to drive the plates depends on a balance between the buoyancy forces driving subduction and the mechanical and buoyancy forces resisting subduction because both the buoyancy and rheology of the slab and mantle depend on temperature composition grain size water content and melt fraction unraveling which of these processes exert a first order control on slab dynamics and under what circumstances other processes become first order effects can be challenging laboratory and numerical models of slab dynamics provide a powerful method for testing the combined effects of buoyancy and strength changes that accompany the slab evolution in the upper mantle transition zone and lower mantle recent studies have focused on understanding how rheologic variations newtonian versus non newtonian viscosity or water content geometry 2d versus 3d and plate motions trench roll back or advance influence the evolution of slabs in the upper mantle and how they sink into the lower mantle these models suggest that spatial and temporal variations in slab strength and the history of subduction determine whether slabs sink directly into the lower mantle or are trapped in the transition zone

regional compensation of subducted lithosphere effects on geoid gravity and topography from a preliminary model abstract it is proposed that the anomalous mass of subducted lithospheric slabs is partially compensated by two mechanisms a by stresses transmitted up the slab to the surface lithosphere causing the trench and b by stresses transmitted through the mantle wedge overlying the slab which cause a broad shallow depression results calculated from a preliminary two dimensional model show that it can give a first order account of observed topography gravity and geoid anomalies in particular the model can reconcile the large mass anomaly predicted for the slab by thermal conduction theory with the smaller than expected gravity and geoid anomalies preliminary results indicate substantial compensation by mechanism a implying rather large stresses at the top of the slab these could be important for theories of earthquake mechanisms and plate driving mechanisms roughly 50 of the possible compensation by mechanism b is indicated the resulting depression the order of 500 m deep and over 500 km wide could explain anomalously deep marginal basins and some marine transgressions of continents

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

climatological characteristics and orographic enhancement of lake effect precipitation east of lake ontario and over the tug hill plateau abstractlake effect snowstorms east of lake ontario are frequently intense and contribute to substantial seasonal accumulations especially over the tug hill plateau hereafter tug hill which rises at a gentle 1 25 slope to 500 m above lake level using a variety of datasets including radar imagery from the ktyx fort drum new york wsr 88d this paper examines the characteristics of lake effect precipitation east of lake ontario over 13 cool seasons 16 september 2001 15 may 2014 during this period days with at least 2 h of lake effect account for 61 76 of the mean cool season snowfall and 24 37 of the mean cool season liquid precipitation mean monthly lake effect frequency and snowfall peak in december and january the highest lake effect frequency and snowfall occur over the western and upper tug hill with an arm of relatively high lake effect frequency and snowfall extending to the southeast shore of lake ontario to the east lee lake effect frequency and snowfall decrease abruptly ov

subduction back arc spreading and global mantle flow the shapes and orientations of benioff zones beneath island arcs interpreted as marking the location of subducted lithosphere provide the best presently available constraints on the global convective flow pattern associated with plate motions this global flow influences the dynamics of subduction subduction zone phenomena therefore provide powerful tests for models of mantle flow we compute global flow models which while simple include those features which are best constrained namely the observed plate velocities applied as boundary conditions and the density contrasts given by thermal models of the lithosphere and subducted slabs two viscosity structures are used for one flow is confined to the upper mantle while for the other flow extends throughout the mantle instantaneous flow velocity vectors match observed benioff zone dips and shapes for the model which allows mantle wide flow but not for the upper mantle model which has a highly contorted flow pattern the effect of trench migration on particle trajectories is calculated it is not important if subduction velocities are greater than migration rates two dimensional finite element models show that including a coherent high viscosity slab does not change these conclusions a coherent high viscosity slab extending deep into the upper mantle would significantly slow subduction if flow were confined to the upper mantle the maximum earthquake magnitude m w for island arcs correlates well with the age of the subducted slab and pressure gradient between the trench and back arc region for the whole mantle but not the upper mantle flow model the correlations with orientations of benioff zones and seismic coupling strongly suggest that the global return flow associated with plate motions extends below 700 km for both models regions of back arc spreading have asthenospheric shear pulling the back arc toward the trench regions without back arc spreading have the opposite sense of shear suggesting global flow strongly influences back arc spreading

orographic influences on a great salt lake effect snowstorm abstractalthough several mountain ranges surround the great salt lake gsl of northern utah the extent to which orography modifies gsl effect precipitation remains largely unknown here the authors use observational and numerical modeling approaches to examine the influence of orography on the gsl effect snowstorm of 27 october 2010 which generated 6 10 mm of precipitation snow water equivalent in the salt lake valley and up to 30 cm of snow in the wasatch mountains the authors find that the primary orographic influences on the event are 1 foehnlike flow over the upstream orography that warms and dries the incipient low level air mass and reduces precipitation coverage and intensity 2 orographically forced convergence that extends downstream from the upstream orography is enhanced by blocking windward of the promontory mountains and affects the structure and evolution of the lake effect precipitation band and 3 blocking by the wasatch and oquirrh mountains which funnels the flow into the sal

controls on trench topography from dynamic models of subducted slabs a finite element method with constrained elements and lagrange multipliers is used to study tectonic faults in a viscous medium a fault representing the interface between overriding and subducting plates has been incorporated into a viscous flow model of a subduction zone in which both thermal buoyancy and the buoyancy associated with the phase change from olivine to spinel are included the fault causes stress to concentrate in its vicinity giving rise to a weak plate margin and a mobile plate if a power law rheology is used surface dynamic topography with either a newtonian or a power law rheology and with typical subduction zone parameters is characterized by a narrow and deep trench and a broadly depressed back arc basin this suggests that oceanic trenches and back arc basins over subduction zones are dynamically compensated by viscous flow our models show that trench depth increases with fault dip angle slab dip angle slab length and age of oceanic lithosphere just prior to subduction the influence of fault dip angle and age of lithosphere on trench depth is greater than the influence of slab dip angle and slab length these relationships of trench depth versus subduction zone parameters explain well the statistics of observed trench depths for those subduction zones with oceanic lithosphere on both sides of the trench observed trench depths have been corrected for fault and slab dip angles based on the relationships from the dynamic models after correction to a common set of parameters trench depth correlates linearly with age of lithosphere prior to subduction with a slope which is close to what models having high viscosities within the transition zone and lower mantle predict comparison between the trench depths corrected for fault and slab dip angles and model trench depths suggests that the resisting tangential stress on faults in subduction zones may range from 15 mpa to 30 mpa depending on model details

model simulations examining the relationship of lake effect morphology to lake shape wind direction and wind speed abstract idealized model simulations with an isolated elliptical lake and prescribed winter lake effect environmental conditions were used to examine the influences of lake shape wind speed and wind direction on the mesoscale morphology this study presents the first systematic examination of variations in lake shape and the interplay between these three parameters the array of 21 model simulations produced cases containing each of the three classic lake effect morphologies i e vortices shoreline bands and widespread coverage and in some instances the mesoscale circulations were composed of coexisting morphologies located over the lake near the downwind shoreline or inland from the downwind shore as with lake effect circulations simulated over circular lakes the ratio of wind speed u to maximum fetch distance l was found to be a valuable parameter for determining the morphology of a lake effect circulation when variations of lake shape wind speed and wind direction were introduced f

comparison of observations with idealized model results for a method to resolve winter lake effect mesoscale morphology abstract forecasters in the great lakes region have for several decades recognized a general relationship of wind speed and overlake fetch to lake effect snowstorm morphology a recent study using idealized mesoscale model simulations of lake effect conditions over circular and elliptical lakes showed the ratio of wind speed to maximum fetch distance u l may be used to effectively predict lake effect snowstorm morphology the current investigation provides an assessment of the u l criteria using observational datasets previously published great lakes lake effect snowstorm observational studies were used to identify events of known mesoscale morphology hindcasts of nearly 640 lake effect events were performed using historical observations with u l as the predictor results show that the quantity u l contains important information on the different mesoscale lake effect morphologies however it provides only a limited benefit when being used to predict mesoscale morphology in real lake effect situations

multiscale dynamics of the tonga kermadec subduction zone our understanding of mantle convection and the motion of plates depends intimately on our understanding of the viscosity structure of the mantle while geoid and gravity observations have provided fundamental constraints on the radial viscosity structure of the mantle the influence of short wavelength variations in viscosity is still poorly understood we present 2 d and 3 d finite element models of mantle flow including strong lateral viscosity variations and local sources of buoyancy owing to both thermal and compositional effects we first use generic 2 d models of a subduction zone to investigate how different observations depend on various aspects of the viscosity structure in particular the slab and lower mantle viscosity and the presence of a low viscosity region in the mantle wedge above the slab we find that 1 the strain rate provides a strong constraint on the absolute viscosity of the slab 10 23 pa s 2 stress orientation within the slab is sensitive to the relative viscosity of the slab lower mantle and the wedge and 3 the stress state and topography of the overriding plate depend on the wedge viscosity and local sources of buoyancy in particular the state of stress in the overriding plate changes from compression to extension with the addition of a low viscosity wedge we then use observations of strain rate stress orientation dynamic topography and the geoid for the tonga kermadec subduction zone as simultaneous constraints on the viscosity and buoyancy in a 3 d regional dynamic model together these observations are used to develop a self consistent model of the viscosity and buoyancy by taking advantage of the sensitivity of each observation to different aspects of the dynamics over a broad range of length scales the presence of a low viscosity wedge makes it possible to match observations of shallow dynamic topography and horizontal extension within the backarc and down dip compression in the shallow portion of the slab these results suggest that a low viscosity wedge plays an important role in controlling the presence of backarc spreading however for a model with a low viscosity and low density region that provides a good fit to the observed topography we find that a reduction of the slab density by a factor of 1 3 relative to the reference density model is required to match the observed geoid these results suggest that compensation of the slab by dynamic topography may be a much smaller effect at short to intermediate wavelengths than predicted by long wavelength modelling of the geoid

the influence of temperature and depth dependent viscosity on geoid and topography profiles from models of mantle convection abstract we consider the geoid and topography profiles from convection models with depth and temperature dependent viscosity models in a 2 d cartesian geometry we use a very low viscosity region in the top center of the geometry to generate plate like viscosity models we focus on subduction regions which have narrow 100 200 km trenches the narrow trench is only a minor feature in the regional geoid in depth dependent viscosity models we change the geometry of the weak zone around the ocean trench to study the impact of the weak zone on trench topography and the geoid we also vary the vertical viscosity structure of the models the weak zone geometry has very little effect on the surface velocity of the ocean plate or the large scale pattern of convection the width of the weak zone within a reasonable range has little effect on topography and geoid profiles there is a narrow trough in the geoid and topography profiles this trough becomes wider and deeper with a high viscosity upper mantle this may suggest that the high viscosity upper mantle is not suitable for our models in temperature dependent viscosity models we vary the viscosity values for different layers temperature dependent viscosity is also important to producing realistic subduction models however the activation energy must be weaker than laboratory estimates of olivine under mantle pressures and temperatures from both depth and temperature dependent viscosity models the viscosity in the lower mantle should be at least ten times greater than in the upper mantle consistent with previous studies of the geoid in regions of subduction the topography and geoid profiles can match the observations when the viscosity increases with depth

influence of peruvian flat subduction dynamics on the evolution of western amazonia convection in the earth s mantle is mainly driven by cold dense subducting slabs but relatively little is known about how 3d variations in slab morphology and buoyancy affect mantle flow or how the surface above deforms in response i e dynamic topography we investigate this problem by studying the dynamics of an active region of flat slab subduction located in peru in south america here the slab geometry is well known based on the regional seismicity and we have observations from the local geological record to validate our models of particular interest is the widespread subsidence and deposition of the solimoes formation across western amazonia that coincided with the development of the peruvian flat slab during the mid late miocene this formation covers an extensive area from the foredeep to the purus arch located view the mathml source 2000km away from the trench close to the andes the preservation of several kilometers of sedimentary thicknesses can be easily accounted for by flexure based on an estimate of the andean loading we predict 2 8 to 3 6 km of accommodation space that spans 100 km the spatial and temporal history of the solimoes formation however particularly the thick distal foreland accumulations up to 1 2 km deep can only be matched with the addition of a longer wavelength dynamic source of topography following the transition from normal to flat subduction we predict over 1 km of dynamic subsidence view the mathml source 1500km wide that propagates over 1000 km away from the trench tracking the subduction leading edge this is followed by a pulse of dynamic uplift over the flat segment behind it we therefore propose that a combination of uplift flexure and dynamic topography during slab flattening in peru is responsible for the sedimentation history and landscape evolution of western amazonia that eventually led to the configuration of the amazon drainage basin we know today

unraveling topography around subduction zones from laboratory models the relief around subduction zones results from the interplay of dynamic processes that may locally exceed the iso static contributions the viscous dissipation of the energy in and around subduction zones is capable of generating kilometer scale vertical ground movements in order to evaluate dynamic topography in a selfconsistent subduction system we carried out a set of laboratory experiments wherein the lithosphere and mantle are simulated by means of newtonian viscous materials namely silicone putty and glucose syrup models are kept in their most simple form and are made of negative buoyancy plates of variable width and thickness freely plunging into the syrup the surface of the model and the top of the slab are scanned in three dimensions a forebulge systematically emerges from the bending of the viscous plate adjacent to the trench with a large wavelength dynamic pressure offsets the foreside and backside of the slab by 500 m on average the suction that accompanies the vertical descent of the slab depresses the surface on both sides at a distance equal to the half width of the slab the topographic depression amounts to 500 m on average and becomes negligible at a distance that equals the width of the slab in order to explore the impact of slab rollback on the topography the trailing edge of the plates is alternatively fixed to fixed mode and freed from free mode the end wall of the tank both the pressure and suction components of the topography are 30 lower in the free mode indicating that slab rollback fosters the dynamic subsidence of upper plates our models are compatible with first order observations of the topography around the east scotia tonga kermadec and banda subduction zones which exhibit anomalous depths of nearly 1 km as compared to adjacent sea floor of comparable age

convective evolution across lake michigan during a widespread lake effect snow event abstract lake effect snowstorms generally develop within convective boundary layers which are induced when cold air flows over relatively warm lakes in fall and winter mesoscale circulations within the boundary layers largely control which communities near the downwind shores of the lakes receive the most intense snow the lack of quantitative observations over the lakes during lake effect storms limits the ability to fully understand and predict these mesoscale circulations this study provides the first observations of the concurrent spatial and temporal evolution of the thermodynamic and microphysical boundary layer structure and mesoscale convective patterns across lake michigan during an intense lake effect event observations analyzed in this study were taken during the lake induced convection experiment lake ice aircraft and sounding observations indicate that the lake effect snows of 13 january 1998 developed within a convective boundary layer that grew rapidly across lake michigan boundary

topography and geoid induced by a convecting mantle beneath an elastic lithosphere summary in the absence of seismological measurements observations of the topography and gravity fields of solid planets are the primary constraints on their internal structure to compute the synthetic geoid and topography induced by the dynamics of planetary interiors we introduce a 3 d numerical tool describing mantle convection beneath an elastic lithosphere although the energy conservation is treated in the whole spherical domain the deformation aspect is solved using a hybrid technique finite volume method for the viscous flow spectral method for elastic deformation the mechanical coupling is achieved via the imposition of the traction at the surface of the viscous flow as a basal boundary condition for the elastic deformation we present both response functions and full thermal convection cases computed with our new method for planetary bodies of varying dimensions the filtering effect of the lithosphere on the dynamic topography and geoid is specific for each planetary body justifying the importance of such a tool furthermore since our approach specifically focuses on the mechanical coupling at the base of the lithosphere it will permit future more elaborate rheological treatments it also enables to discriminate between the radial and tangential components of the viscous traction the latter is found to have a significant influence on the elastic deformation the effect on geoid is prominent more specifically while a thin elastic lithosphere is usually considered to play little role on the dynamic topography and geoid of venus a 35 per cent reduction is obtained for geoid height in the numerical example we propose on a planet with thicker elastic lithosphere such as mars the consequence of this filtering effect is to rule out the possibility of a dynamical support for the tharsis rise even for the lowest admissible values of elastic thickness in this region

formation of the maturin foreland basin eastern venezuela thrust sheet loading or subduction dynamic topography 1 the maturin basin in eastern venezuela is considered a good example of a peripheral foreland basin earthquake and tomographic data indicate that eastern venezuela is affected by the oblique subduction of the south american plate underneath the caribbean plate new forward flexural isostatic modeling of eastern venezuela has been carried out in order to determine whether the maturin basin was generated purely by thrust sheet loading from the serrania and monagas foreland thrust belts a sequence of forward models from middle miocene to present was generated for 3 profiles across the serrania del interior thrust belt the monagas foreland thrust belt and the maturin foreland basin the predictions of these models are constrained using seismic reflection and well data the flexural isostatic modeling shows that thrust sheet loading associated with the serrania del interior and monagas foreland thrust belts is insufficient to generate the observed subsidence within the maturin basin dynamic fluid flow modeling of subduction related dynamic topography of eastern venezuela has been used to investigate the influence of south american plate subduction on the generation of the accommodation space observed in the maturin basin fluid flow modeling of subduction related dynamic topography suggests that the subduction of the south american lithospheric mantle caused downward deflection of the south american crust affecting the maturin basin and the serrania thrust belt this modeling suggests that the maturin basin subsidence has two components 55 related to thrust sheet loading and 45 driven by continental subduction

topography of the overriding plate during progressive subduction a dynamic model to explain forearc subsidence overriding plate topography provides constraints on subduction zone geodynamics we investigate its evolution using fully dynamic laboratory models of subduction with techniques of stereoscopic photogrammetry and particle image velocimetry model results show that the topography is characterized by an area of forearc dynamic subsidence with a magnitude scaling to 1 44 3 97 km in nature and a local topographic high between the forearc subsided region and the trench these topographic features rapidly develop during the slab free sinking phase and gradually decrease during the steady state slab rollback phase we propose that they result from the variation of the vertical component of the trench suction force along the subduction zone interface which gradually increases with depth and results from the gradual slab steepening during the initial transient slab sinking phase the downward mantle flow in the nose of the mantle wedge plays a minor role in driving forearc subsidence

sensitivity of the short to intermediate wavelength geoid to rheologic structure in subduction zones 1 the rheologic structure of subduction zones determines how the negative buoyancy of sinking lithosphere slabs couples to the surface to cause surface deformation and plate motion however the rheologic structures determined using different methods are inconsistent analysis of the geoid and some dynamical models conclude that slabs must be weak whereas other dynamical models and experimentally determined rheology predict that slabs are strong using idealized three dimensional models of a subduction zone we show that lateral variations in viscosity lvvs within a slab caused by yielding render the geoid insensitive to the viscosity structure outside of the bending region of the slab spectral analysis of the geoid up to harmonic degree 360 features at wavelengths down to 110 km shows that lvvs cause significant mode coupling with changes to the long wavelength geoid 1000 km caused by short wavelength 500 km changes to dynamic topography models with either a uniformly weak slab or a strong slab with yielding provide a good match to the observed geoid spectrum and characteristic features of the geoid over subduction zones therefore observations other than the geoid are needed to constrain the strength profile of subducting slabs