comparison of dynamic flow models for the central aleutian and tonga kermadec subduction zones using three dimensional 3 d instantaneous dynamic models of the central aleutian subduction zone we show that a low viscosity low density deltarho 10 kg m3 region in the wedge is needed to match observations of topography the geoid and the stress state in the slab and overriding plate previous models of the tonga kermadec subduction zone also require a low viscosity low density deltarho 20 kg m3 region in the wedge a low viscosity region in the wedge decouples the overriding plate from the downward flow in the wedge caused by the sinking slab reducing the magnitude of negative dynamic topography on the overriding plate and changing the state of stress in the slab from down dip tension to down dip compression if the low viscosity region is large as in the tonga kermadec models then the decoupling may also cause the overriding plate to go into extension these are the first dynamic models to simultaneously match the state of stress in the slab and overriding plate and suggest that the combination of stress state within the slab and overriding plate may provide an additional constraint on the viscosity structure within subduction zones

mesozoic aseismic ridges on the farallon plate and southward migration of shallow subduction during the laramide orogeny we propose that an aseismic ridge of late jurassic and early cretaceous age on the farallon plate was subducted beneath the unites states cordillera during the laramide orogeny the relative buoyancy of the aseismic ridge caused shallow subduction which resulted in a magmatic lull 70 40 ma in the existing near trench magmatic arc and laramide uplift and faulting of crystalline basement 1000 1500 km inland from the trench the timing and latitudinal limits of the subduction of the aseismic ridge are similar to the timing and latitudinal limits of both the laramide orogeny and the magmatic lull in the western cordillera because the north northeast trending aseismic ridge was moving northeast relative to the north american plate the point of collision of the aseismic ridge with north america and therefore the locus of shallow subduction migrated southward with time causing the well documented southward migration of the magmatic lull

mesozoic plate motion history below the northeast pacific ocean from seismic images of the subducted farallon slab the high resolution seismic imaging of subducted oceanic slabs1 2 has become a powerful tool for reconstructing palaeogeography3 the images can now be interpreted quantitatively by comparison with models of the general circulation of the earth s mantle4 here we use a three dimensional spherical computer model of mantle convection5 6 to show that seismic images of the subducted farallon plate provide strong evidence for a mesozoic period of low angle subduction under north america such a period of low angle subduction has been invoked independently to explain rocky mountain uplift far inland from the plate boundary during the laramide orogeny7 the computer simulations also allow us to locate the largely unknown kula farallon spreading plate boundary the location of which is important for inferring the trajectories of suspect terrain across the pacific basin8

localization of the gravity field and the signature of glacial rebound the negative free air gravity anomaly centred on hudson bay canada shows a remarkable correlation with the location of the laurentide ice sheet suggesting that this gravity anomaly is the result of incomplete post glacial rebound this region however is also underlain by higher than average mantle seismic velocities suggesting that the gravity low might result instead from dynamic topography associated with convective downwellings here we analyse the global gravity field as a simultaneous function of geographic location and spectral content we find that the hudson bay gravity low is unique with anomalously high amplitude in the spectral band where the power from the laurentide ice load is greatest and the relaxation times predicted for viable models of viscous relaxation are longest we estimate that about half of the hudson bay gravity anomaly is the result of incomplete post glacial rebound and derive a mantle viscosity model that explains both this gravity signature and the characteristic uplift rates for the central laurentide and fennoscandian regions this model has a jump in viscosity at 670 km depth comparable to that in dynamic models of the geoid highs over subducted slabs but lacks a low viscosity asthenosphere consistent with a higher viscosity in the upper mantle beneath shields than in oceanic regions

long term sea level fluctuations driven by ocean basin dynamics earth9s long term sea level history is characterized by widespread continental flooding in the cretaceous period 145 to 65 million years ago followed by gradual regression of inland seas however published estimates of the late cretaceous sea level high differ by half an order of magnitude from 40 to 250 meters above the present level the low estimate is based on the stratigraphy of the new jersey margin by assimilating marine geophysical data into reconstructions of ancient ocean basins we model a late cretaceous sea level that is 170 85 to 270 meters higher than it is today we use a mantle convection model to suggest that new jersey subsided by 105 to 180 meters in the past 70 million years because of north america9s westward passage over the subducted farallon plate this mechanism reconciles new jersey margin based sea level estimates with ocean basin reconstructions

simultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection through the assimilation of present day mantle seismic structure adjoint methods can be used to constrain the structure of the mantle at earlier times i e mantle initial conditions however the application to geophysical problems is restricted through both the high computational expense from repeated iteration between forward and adjoint models and the need to know mantle properties such as viscosity and the absolute magnitude of temperature or density a priori we propose that an optimal first guess to the initial condition can be obtained through a simple backward integration sbi of the governing equations thus lessening the computational expense given a model with known mantle properties we show that a solution based on an sbi generated first guess has smaller residuals than arbitrary guesses mantle viscosity and the effective rayleigh number are crucial for mantle convection models neither of which is exactly known we place additional constraints on these basic mantle properties when the convection induced dynamic topography on earth s surface is considered within an adjoint inverse method besides assimilating present day seismic structure as a constraint we use dynamic topography and its rate of change in an inverse method that allows simultaneous inversion of the absolute upper and lower mantle viscosities scaling between seismic velocity and thermal anomalies and initial condition the theory is derived from the governing equations of mantle convection and validated by synthetic experiments for both one layer viscosity and two layer viscosity regionally bounded spherical shells for the one layer model at any instant of time the magnitude of dynamic topography is controlled by the temperature scaling while the rate of change of topography is controlled by the absolute value of viscosity for the two layer case the rate of change of topography constrains upper mantle viscosity while the magnitude of dynamic topography determines the temperature scaling lower mantle viscosity when upper mantle lower mantle density anomaly dominates the flow field this two stage scheme minimizes the tradeoff between temperature and lower mantle viscosity for both cases we show that the theory can constrain mantle properties with errors arising through the adjoint recovery of the initial condition for the two layer model this error is manifest as a tradeoff between the temperature scaling and lower mantle viscosity

the ups and downs of north america evaluating the role of mantle dynamic topography since the mesozoic the driving force for transient vertical motions of earth s surface remains an outstanding question a main difficulty lies in the uncertain role of the underlying mantle especially during the geological past here i review previous studies on both observational constraints and physical mechanisms of north american topographic evolution since the mesozoic i first summarize the north american vertical motion history using proxies from structural geology geochronology sedimentary stratigraphy and geomorphology based on which i then discuss the published physical models overall there is a progressive consensus on the contribution of mantle dynamic topography due to buoyancy structures associated with the past subduction at the continental scale a largely west to east migrating deformation pattern suggests an eastward translation of mantle dynamic effects consistent with models involving an eastward subduction and sinking of former farallon slabs since the cretaceous among the existing models the inverse model based on an adjoint algorithm and time dependent data constraints provides the most extensive explanations for the temporal changes of north american topography since the mesozoic at regional scales debates still exist on the predicted surface subsidence and uplift within both the western and eastern united states where discrepancies are likely due to differences in model setup e g mantle dynamic properties and boundary conditions and the amount of time dependent observational constraints toward the development of the next generation predictive geodynamic models new research directions may include 1 development of enhanced data assimilation capabilities 2 exploration of multiscale and multiphysics processes and 3 cross disciplinary code coupling

influence of climate change and uplift on colorado plateau paleotemperatures from carbonate clumped isotope thermometry the elevation history of earth s surface is key to understanding the geodynamic processes responsible for the rise of plateaus we investigate the timing of colorado plateau uplift by estimating depositional temperatures of tertiary lake sediments that blanket the plateau interior and adjacent lowlands using carbonate clumped isotope paleothermometry a measure of the temperature dependent enrichment of 13 c 18 o bonds in carbonates comparison of modern and ancient samples deposited near sea level provides an opportunity to quantify the influence of climate and therefore assess the contribution of changes in elevation to the variations of surface temperature on the plateau analysis of modern lake calcite from 350 to 3300 m elevation in the southwestern united states reveals a lake water carbonate temperature lct lapse rate of 4 2 0 6 c km analysis of miocene deposits from 88 to 1900 m elevation in the colorado river drainage suggests that the ancient lct lapse rate was 4 1 0 7 c km and temperatures were 7 7 2 0 c warmer at any one elevation than predicted by the modern trend the inferred cooling is plausible in light of pliocene temperature estimates off the coast of california and the consistency of lapse rates through time supports the interpretation that there has been little or no elevation change for any of the samples since 6 ma together with previous paleorelief estimates from apatite u th he data from the grand canyon our results suggest most or all of the plateau s lithospheric buoyancy was acquired 80 60 ma and do not support explanations that ascribe most plateau uplift to oligocene or younger disposal of either the farallon or north american mantle lithosphere

topographic asymmetry of the south atlantic from global models of mantle flow and lithospheric stretching the relief of the south atlantic is characterized by elevated passive continental margins along southern africa and eastern brazil and by the bathymetric asymmetry of the southern oceanic basin where the western flank is much deeper than the eastern flank we investigate the origin of these topographic features in the present and over time since the jurassic with a model of global mantle flow and lithospheric deformation the model progressively assimilates plate kinematics plate boundaries and lithospheric age derived from global tectonic reconstructions with deforming plates and predicts the evolution of mantle temperature continental crustal thickness long wavelength dynamic topography and isostatic topography mantle viscosity and the kinematics of the opening of the south atlantic are adjustable parameters in thirteen model cases model predictions are compared to observables both for the present day and in the past present day predictions are compared to topography mantle tomography and an estimate of residual topography predictions for the past are compared to tectonic subsidence from backstripped borehole data along the south american passive margin and to dynamic uplift as constrained by thermochronology in southern africa comparison between model predictions and observations suggests that the first order features of the topography of the south atlantic are due to long wavelength dynamic topography rather than to asthenospheric processes the uplift of southern africa is best reproduced with a lower mantle that is at least 40 times more viscous than the upper mantle

flexural subsidence and basement tectonics of the cretaceous western interior basin united states the flexural subsidence history recorded in cenomanian to early campanian 97 to 80 ma strata in the cretaceous u s western interior basin was studied with two dimensional flexural backstripping techniques results indicate that the flexural subsidence resulting from thrust loading was superimposed on epeirogenic subsidence in the foreland basin the flexural component exhibits significant spatial and temporal variations along both the strike and dip relative to the sevier thrust belt the greatest cumulative subsidence occurred in southwestern wyoming and northern utah concurrent subsidence in northwestern montana and southern utah was insignificant temporal trends in subsidence also show a distinct regional pattern from the cenomanian to late turonian 97 to 90 ma subsidence rates were high in utah and much lower in wyoming and montana in contrast during the coniacian and santonian 90 to 85 ma subsidence accelerated rapidly in wyoming increased slightly in montana and decreased in utah we suggest that these spatially and temporally varying subsidence patterns reflect the interplay of several geodynamic factors including 1 temporal and spatial variation in emplacement of the thrust loads 2 segmentation of the basement into adjacent blocks with different rheological properties 3 reactivation of basement fault trends and 4 regional dynamic topographic effects

mantle circulation models with variational data assimilation inferring past mantle flow and structure from plate motion histories and seismic tomography summary mantle convection models require an initial condition some time in the past because this initial condition is unknown for earth we cannot infer the geological evolution of mantle flow from forward mantle convection calculations even for the most recent mesozoic and cenozoic geological history of our planet here we introduce a fluid dynamic inverse problem to constrain unknown mantle flow back in time from seismic tomographic observations of the mantle and reconstructions of past plate motions using variational data assimilation we derive the generalized inverse of mantle convection and explore the initial condition problem in high resolution 3 d spherical mantle circulation models for a time period of 100 myr roughly comparable to half a mantle overturn we present a synthetic modelling experiment to demonstrate that mid cretaceous mantle structure can be inferred accurately from fluid dynamic inverse modelling assuming present day mantle structure is well known even if an initial first guess assumption about the mid cretaceous mantle involved only a simple 1 d radial temperature profile we also demonstrate that convecting present day mantle structure back in time by reversing the time stepping of the energy equation is insufficient to model the mantle structure of the past the difficulty arises because such backward convection calculations ignore thermal diffusion effects and therefore cannot account for the generation of thermal buoyancy in boundary layers as we go back in time inverse mantle convection modelling should make it possible to infer a number of flow parameters from observational constraints of the mantle

density of morb eclogite in the upper mantle in order to study the role of eclogite subducted oceanic crust in geodynamics we carried out high pressure experiments on natural mid oceanic ridge basalt morb using a multi anvil apparatus we determined phase relations of morb as a function of pressure and temperature in the temperature range of 1200 1700 c and pressures from 5 to 14 gpa the present results differ from a previous high pressure study of morb earth planet sci lett 77 1986 245 conducted at 1200 c as follows 1 coesite coe or stishovite st appears throughout the studied pressure range and 2 the clinopyroxene garnet cpx gar ratio is smaller particularly in the low pressure range 12 gpa our results show that the calculated bulk density of subducted oceanic crust increases discontinuously at around 9 gpa due to the coesite stishovite phase transition the density is considerably higher than estimates based on the previous experiments at pressures less than 12 gpa thus the role of eclogite in driving subduction should be more important than previously thought on the other hand the density difference between peridotite and eclogite in ascending mantle plumes becomes smaller as temperature increases a result of the increase in clinopyroxene garnet sio2 ratio in eclogite with increasing temperature as eclogite is still denser than peridotite even under subsolidus conditions recycled oceanic crust entrained in ascending plumes must be sustained by the upward flow of ambient peridotite in mantle plumes

mantle provinces under north america from multifrequency p wave tomography this is a survey of mantle provinces large scale seismic anomalies under north america from the surface down to 1500 1800 km depth the underlying p velocity model was obtained by multifrequency tomography a waveform based method that systematically measures and models the frequency dependence of teleseismic body waves a novel kind of three dimensional rendering technique is used to make the considerable structural complexities under north america accessible in the transition zone and below the north american mantle is dominated by seismically fast provinces which represent distinct subduction episodes of the farallon plate i attempt to date and interpret the various slab fragments by reconciling their present positions with paleotrench locations from plate tectonic reconstructions and with major geologic surface episodes differences in vertical sinking velocity have led to large vertical offsets across adjacent coeval slabs some of the mantle provinces have not been discussed much previously including 1 a seismically slow blanket overlying the oldest farallon subduction along the eastern continental margin 2 a transition zone slab coeval with the laramide orogeny ca 80 60 myr which i discuss in analogy to the stagnant slab subduction style commonly found in the western pacific today 3 the lower mantle root of present day cascadia subduction which may have started out as intraoceanic subduction 4 a lower mantle slab under arizona and new mexico the last material to subduct before strike slip motion developed along the san andreas boundary and 5 two narrow plate tears thousands of kilometers long one of which is the subducted conjugate of the mendocino fracture zone

characteristics of instabilities in the mesopause region over maui hawaii 1 characteristics of convective and dynamical instabilities in the mesopause region between 85 and 100 km over maui hawaii 20 7 n 156 3 w are investigated using 19 nights 133 hours of high resolution wind and temperature data obtained by the university of illinois na wind temperature lidar during the maui mesosphere and lower thermosphere maui malt campaigns the mean probabilities of convective and dynamical instabilities are observed to be 3 and 10 respectively but there is considerable night to night variation at any given time the probability that an unstable condition is found at some altitudes in the 85 100 km range is 90 the maui malt data exhibit a distinct trend for n2 to increase with wind shear and vice versa this correlation has important implications in the understanding of the development of instabilities the night of 11 april 2002 is studied in detail in order to investigate the spatial and temporal structures of n2 wind shear and convective and dynamical instabilities a close linkage between instability and the mesosphere inversion layers mils is identified most of the convectively and dynamically unstable regions are located above the mils with a tendency for dynamical instability to develop below convective instability it is found that the vertical variations of n2 are often correlated with those of wind shear but with a phase shift such that the maxima and minima of n2 are located 0 5 1 km below those of wind shear because of this shift dynamical instability tends to develop in the region above the maximum wind shear where relatively small n2 is observed to be associated with large wind shear we also found that the wind shear is dominated by the contribution of the meridional wind especially when the wind shear is strong possible mechanisms for the observed features are discussed

eroding dynamic topography 1 geological observations of mantle flow driven dynamic topography are numerous especially in the stratigraphy of sedimentary basins on the contrary when it leads to subaerial exposure of rocks dynamic topography must be substantially eroded to leave a noticeable trace in the geological record here we demonstrate that despite its low amplitude and long wavelength and thus very low slopes dynamic topography is efficiently eroded by fluvial erosion providing that drainage is strongly perturbed by the mantle flow driven surface uplift using simple scaling arguments as well as a very efficient surface processes model we show that dynamic topography erodes in direct proportion to its wavelength we demonstrate that the recent deep erosion experienced in the colorado plateau and in central patagonia is likely to be related to the passage of a wave of dynamic topography generated by mantle upwelling

an inversion for radial viscosity structure using seismic tomography using a uniform velocity density scaling we have inverted the seismically inferred 3 dimensional structure of the whole manfie for the radial viscosity structure which best fits the geoid we are able to obtain 60 72 variance reductions for three different s wave tomographic models the resulting viscosity structures are remarkably similar showing a high viscosity lid a low viscosity zone in the transition region and a high viscosity lower manfie a resolution analysis indcates that the viscosity structure in the upper manfie is well resolved by the data however the resolution in the lower manfie is poorer our models axe in general agreement with previous studies except that our inversions prefer a low viscosity layer at 400 670 km as opposed to 100 400 kin

complex subduction and small scale convection revealed by body wave tomography of the western united states upper mantle abstract new high resolution p and s wave tomography of the united states upper mantle from the pacific coast to the great plains reveals strong multi scale heterogeneity closely correlated with tectonic and magmatic activity we invert teleseismic travel time residuals from the earthscope transportable array and more than 1700 additional temporary and permanent stations for 3 d velocity perturbations to a depth of 1000 km the inversion uses recent advances in western u s crust models to better isolate the mantle component of travel time residuals and frequency dependent 3 d sensitivity kernels to map travel time residuals measured in multiple frequency bands into velocity structure in addition to separate v p and v s models we jointly invert the two datasets for v p v s perturbations by imposing a smoothness constraint on the ln v s ln v p field the joint inversion helps us identify regions where partial melt is probable the amplitude of v p v s and v p v s variations is greatest in the upper 200 km of the mantle and the form of velocity anomalies suggests a provincially heterogeneous lithosphere and the occurrence of widespread small scale convection partially molten mantle is inferred beneath yellowstone and the eastern snake river plain srp the salton trough and the clear lake volcanic field the inferred depth extent of partial melt is consistent with a generally hydrated upper mantle and elevated temperatures beneath the eastern srp and yellowstone despite continuous subduction since the cretaceous the distribution of sub lithospheric high velocity anomalies is dissected similar to other recent studies based on our new tomography models western u s geologic history and plate tectonic reconstructions we infer patchy and incomplete removal of the flat subducting laramide slab and slab tearing associated with eocene accretion in the northwestern u s

laramide orogenic influence on late mesozoic cenozoic subsidence history western deep gulf of mexico basin the deep part of the gulf of mexico basin which is underlain by oceanic crust has generally been regarded as a tectonically stable part of the basin since its formation by sea floor spreading and thermal subsidence during the late jurassic early cretaceous subsequent subsidence is believed to have been caused mainly by thermal cooling and sediment loading a backstripping study however has revealed a tectonic event in the western deep basin with tectonic loading effect today one possible explanation is laramide fold thrust loading along the western flank of the basin during 66 40 ma

a record of appalachian denudation in postrift mesozoic and cenozoic sedimentary deposits of the u s middle atlantic continental margin abstract the complex interplay between source terrain uplift basin subsidence paleoclimatic shifts and sea level change left an extensive sedimentary record in the contiguous offshore basins of the u s middle atlantic margin salisbury embayment baltimore canyon trough and hatteras basin isopach maps of 23 postrift lower jurassic to quaternary a allostratigraphic units coupled with a revised stratigraphic framework reveal that tectonism by regulating sediment supply accumulation rate dominated the interplay of forcing mechanisms tectonic pulses are evidenced by abruptly accelerated sediment accumulation marked latitudinal shifts in the location of depocenters and regional changes in lithofacies relatively rapid tectonic subsidence during the early and middle jurassic history of the basins may have enhanced sediment accumulation rates beginning in the late jurassic however subsidence rates decreased significantly though occasional short pulses of subsidence may have effected relative sea level rises sea level change heavily influenced the distribution and redistribution of sediments one they reached the basins and paleoclimate regulated the relative abundance of carbonates and evaporites in the basins we conclude that source terrains of the central appalachian highlands were tectonically uplifted intensely weathered and rapidly eroded three times since the late triassic 1 early to middle jurassic aalenian to callovian 2 mid early cretaceous barremian and 3 late cenozoic middle miocene intervals of tectonic quiescence following these three tectonic pulses provided conditions suitable for the formation of regional erosion surfaces geomorphic features commonly reported to characterize the central appalachian highlands this series of three irregularly spaced tectonic quiescent cycles does not however match the traditional four cycle concept of post triassic appalachian peneplanation

segmentation of the farallon slab abstract recent tomography images reveal a complex 3d mantle structure beneath western united states with feature morphology varying rapidly with depth by assimilating plate motion history paleo age of sea floor and paleo geography of plate boundaries in a 3 d numerical model we simulate the farallon juan de fuca subduction during the past 40 ma we find that the highly segmented upper mantle structure of western u s is a direct result of the farallon subduction we show that the tilted horseshoe shaped fast seismic anomaly beneath nevada and utah at 300 600 km depth range is in fact a segment of curled slab subducted since 15 ma and the shallower linear slab beneath the cascades is younger than 5 ma the distinct morphology between these two parts of the subduction system indicates the strong influence of the fast trench rollback since 20 ma the northward migrating jf pa na triple junction and the toroidal flow around slab edges the observed mantle structures are used to constrain the rheology of the upper mantle through matching the shape depth and location of modeled subducted slab segments the inferred viscosity for the asthenosphere is 5 10 19 pa s and those for the transition zone and lower mantle are 1 5 10 21 pa s and 2 10 22 pa s respectively the slab is found to be about 2 orders of magnitude stronger than the ambient mantle above 410 km depth but of similar order of magnitude viscosity in the transition zone

mass balanced paleogeographic reconstructions paleogeography paleotopography and paleobathymetry of a closed erosion depositional system can be reconstructed by restoring sedimentary masses to elevated surfaces in a drainage basin based on the inverse of present erosion equations and adjusting for isostasy sea level changes sediment compaction and thermal subsidence the erosion deposition history of the northwestern gulf of mexico margin and its western central north american source area during the cenozoic is used to explore the sensitivity of mass balance reconstructions and changes in assumptions concerning erosion rate parameters and sea level

intra oceanic subduction shaped the assembly of cordilleran north america the western quarter of north america consists of accreted terranes crustal blocks added over the past 200 million years but the reason for this is unclear the widely accepted explanation posits that the oceanic farallon plate acted as a conveyor belt sweeping terranes into the continental margin while subducting under it here we show that this hypothesis which fails to explain many terrane complexities is also inconsistent with new tomographic images of lower mantle slabs and with their locations relative to plate reconstructions we offer a reinterpretation of north american palaeogeography and test it quantitatively collision events are clearly recorded by slab geometry and can be time calibrated and reconciled with plate reconstructions and surface geology the seas west of cretaceous north america must have resembled today s western pacific strung with island arcs all proto pacific plates initially subducted into almost stationary intra oceanic trenches and accumulated below as massive vertical slab walls above the slabs long lived volcanic archipelagos and subduction complexes grew crustal accretion occurred when north america overrode the archipelagos causing major episodes of cordilleran mountain building a new explanation for the origin of the accreted terranes that form the mountainous cordillera of western north america is proposed and tested stationary intra oceanic subduction deposited massive slab walls in the mantle and grew volcanic archipelagos at the surface which were overridden by and accreted to north america during cretaceous times the generally accepted account of how the geology of western north america developed involves the oceanic farallon plate acting as a conveyor belt sweeping crustal terranes onto the continental margin as it moves beneath the north american plate this model fails to explain many aspects of this accreted terrane and here karin sigloch and mitchell mihalynuk present an alternative model based on new tomographic images of lower mantle slabs and modelling of deep mantle structure in the region they propose that during the cretaceous period there were a series of island arcs west of north america beneath which proto pacific plates subducted into almost stationary trenches long lived volcanic archipelagos and subduction complexes grew above the slabs and then accreted onto north america as it overrode the archipelagos causing episodes of cordilleran mountain building

the dynamics of plate tectonics and mantle flow from local to global scales plate tectonics is regulated by driving and resisting forces concentrated at plate boundaries but observationally constrained high resolution models of global mantle flow remain a computational challenge we capitalized on advances in adaptive mesh refinement algorithms on parallel computers to simulate global mantle flow by incorporating plate motions with individual plate margins resolved down to a scale of 1 kilometer back arc extension and slab rollback are emergent consequences of slab descent in the upper mantle cold thermal anomalies within the lower mantle couple into oceanic plates through narrow high viscosity slabs altering the velocity of oceanic plates viscous dissipation within the bending lithosphere at trenches amounts to 5 to 20 of the total dissipation through the entire lithosphere and mantle

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

testing a coupled ice mixed layer model under subarctic conditions abstract a one dimensional oceanic mixed layer model has been coupled with a thermodynamic sea ice model in order to study the seasonal cycle of ice ocean interaction in the subarctic ocean the ice thickness is assumed constant and only variations of ice compactness under the effect of the atmospheric and oceanic heat fluxes are considered the mixed layer model allows prediction of the rate of entrainment mainly as the result of the penetrative buoyancy driven convection both heat and salt conservation equations are included so that the respective efficiency with regard to vertical mixing of brine rejection below the ice and intense surface cooling in open water can be compared the model has been run for ten years on a zonal vertical section across the greenland sea the annual entrainment retreat cycle of the mixed layer is well reproduced with maximum depths occurring in late winter even down to the bottom at some locations near the ice edge surface cooling is found to be a more efficient mecha

constraining cretaceous subduction polarity in eastern pacific from seismic tomography and geodynamic modeling interpretation of recent mantle seismic images below the america ignited a debate on the cretaceous subduction polarity in the eastern pacific ocean the traditional view is that the massive vertical slab wall under eastern north america resulted from an eastward farallon subduction an alternative interpretation attributes this prominent seismic structure to a westward subduction of the north american plate against a stationary intraoceanic trench here i design quantitative subduction models to test these two scenarios using their implied plate kinematics as velocity boundary conditions modeling results suggest that the westward subduction scenario could not produce enough slab volume as seismic images reveal as is due to the overall slow subduction rate 2 5 cm yr the results favor the continuous eastward farallon subduction scenario which with an average convergence rate of 10 cm yr prior to the eocene can properly generate both the volume and the geometry of the imaged lower mantle slab pile the eastward subduction model is also consistent with most cretaceous geological records along the west coast of north america

reconstructing earth history in three dimensions an inverse model elucidates the connection between plate tectonics at earth s surface and the dynamics of the underlying mantle

p and s wave tomography of the mantle beneath the united states mantle seismic structure beneath the united states spanning from the active western plate margin to the passive eastern margin was imaged with teleseismic p and s wave traveltime tomography including usarray data up to may 2014 to mitigate artifacts from crustal structure 5 40 s rayleigh wave phase velocities were used to create a 3 d starting model major features of the final p and s models include two distinct low velocity anomalies at depths of 60 300 km beneath the central and northern appalachians and passive margin the central appalachian low velocity anomaly coincides with eocene basaltic magmatism and the northern anomaly is located along the cretaceous track of the great meteor hot spot at depths of 300 700 km beneath the central and eastern u s large high velocity anomalies are inferred to be remnants of the farallon slab that subducted prior to 40 ma during the laramide orogeny

an uplift history of the colorado plateau and its surroundings from inverse modeling of longitudinal river profiles 1 it is generally agreed that a region encompassing the colorado plateau has been uplifted by sub crustal processes admittance calculations tomographic studies and receiver function analyses suggest that dynamic support is generated by some combination of convective upwelling and lithospheric thickness changes notwithstanding advances in our understanding of present day setting uplift rate histories are poorly constrained and debated an improved history will aid discrimination between proposed models here we show that a regional uplift rate history can be obtained by inverting longitudinal river profiles we assume that the shape of a river profile is controlled by uplift rate and moderated by erosion in our model uplift rate is allowed to vary smoothly as a function of space and time upstream drainage area is invariant with time simultaneous inversion of river profiles from the colorado rio grande columbia and mississippi catchments shows that three phases of regional uplift occurred the first phase occurred between 80 and 50 myrs when 1 km of uplift was generated at a rate of 0 03 mm yr a second phase occurred between 35 and 15 myrs when 1 5 km of uplift was generated at a faster rate of 0 06 mm yr a final phase of uplift commenced 5 myrs ago these distinct phases of late cretaceous and oligocene uplift are corroborated by stratigraphic considerations by thermochronometric data and by stratigraphic evidence of periodic clastic efflux delivered into the gulf of mexico an episodic uplift history is consistent with staged removal of thick lithospheric mantle beneath a large region which is currently centered on yellowstone

apatite 4he 3he and u th he evidence for an ancient grand canyon the grand canyon is one of the most dramatic features on earth yet when and why it was carved have been controversial topics for more than 150 years here we present apatite 4he 3he thermochronometry data from the grand canyon basement that tightly constrain the near surface cooling history associated with canyon incision 4he 3he spectra for eastern grand canyon apatites of differing he date radiation damage and u th zonation yield a self consistent cooling history that substantially validates the he diffusion kinetic model applied here similar data for the western grand canyon provide evidence that it was excavated to within a few hundred meters of modern depths by 70 million years ago ma in contrast to the conventional model in which the entire canyon was carved since 5 to 6 ma

the role of oceanic plateau subduction in the laramide orogeny the cause of the laramide phase of mountain building remains uncertain conceptual models implicate the subduction of either ocean ridges or conjugates of the buoyant hess or shatsky oceanic plateaux independent verification of these models has remained elusive because the putative ridges or plateaux are no longer at the earth s surface inverse convection models have identified two prominent seismic anomalies on the recovered farallon plate here we combine inverse convection models with reconstructions of plate motions to show that these seismic anomalies coincide palaeogeographically with the restored positions of the shatsky and hess conjugate plateaux as they subducted beneath north america specifically the distribution of laramide crustal shortening events tracked the passage of the shatsky conjugate beneath north america whereas the effects of the hess conjugate subduction were restricted to the northern mexico foreland belt we propose that continued subduction caused the oceanic crust to undergo the basalt eclogite phase transformation during which the shatsky conjugate lost its extra buoyancy and was effectively removed increases in slab density and coupling between the overriding and subducting plates initially dragged the surface downward followed by regional scale surface rebound we conclude that laramide uplift resulted from the removal rather than emplacement of the shatsky conjugate

subsidence in intracontinental basins due to dynamic topography the origin of anomalous tectonic subsidence ats of large intracontinental basins long after their most recent phase of extension and last thermal perturbation is the subject of a long standing debate we show that deep earth processes may contribute to the subsidence of these tectonically stable basins by analysing the tertiary mantle convection driven topography history of a global set of more than 220 intracontinental basins integrated into a plate kinematic framework most basins are affected by increasing negative dynamic topography over the last 70 myr due to the motion of many continents away from large mantle upwellings and towards downwellings previously located along the perimeter of the supercontinent pangaea during continental dispersal increasing negative dynamic topography causes dynamic subsidence of the basins creating additional accommodation space we utilise a parameter from a global crustal analysis of intraplate basins termed anomalous tectonic subsidence to quantify sediment accumulation not related to crustal stretching we propose that dynamic subsidence due to plate motions relative to the underlying mantle as well as variations in the large scale convection patterns can significantly contribute to the creation and destruction of accommodation space in intraplate basins