deep mantle forces and the uplift of the colorado plateau 1 we introduce a quantitative model of global mantle convection that reconstructs the detailed motion of a warm mantle upwelling over the last 30 ma towards the interior of the southwestern usa from observed present day mantle heterogeneity the onset and evolution of uplift in the central basin and range province and colorado plateau during this time is determined by tracking the topographic swell due to this mantle upwelling we show that 1 the extension and basaltic volcanism post 25 ma in the central basin and range coincides with the arrival and eastward progression of this upwelling and 2 dynamic uplift of the southern colorado plateau totaling about 1 km transpired in the last 20 ma since 10 ma the center of uplift continued northeastward from the southwestern rim of the plateau consistent with a young grand canyon model and eastward sweep of magmatism in the western colorado plateau

constraints on the venusian crustal thickness variations in the isostatic stagnant lid approximation 1 the isostatic stagnant lid isl approximation is applied to the long wavelength venusian topography and geoid to place constraints on the thickness of the crust and lithosphere previous studies showed that a purely thermal isostasy model of the lithosphere the stagnant lid accounts for a significant portion of the topography and geoid anomalies here we assume that the remaining discrepancy between the observed and model geoid is due to crustal thickness variations and that the crust the lithosphere and the mantle are in a state of double isostatic equilibrium the lithosphere is isostatically floating on top of the mantle and the crust is isostatically floating on top of the lithosphere the addition of another variable the crustal thickness allows an exact fit to the observed topography and geoid this also makes the solution non unique an exact fit to topography and geoid can be achieved in a broad range of average lithospheric and crustal thicknesses additional constraints such as the gabbro eclogite phase transition and the existence of relatively young volcanism limit this range to some extent the inferred variations in the crustal thickness are less sensitive to model assumptions and generally agree with the results obtained in previous studies

seismic evidence for significant crustal thickening beneath jabalpur earthquake 21 may 1997 source region in narmada son lineament central india 1 we investigate crustal properties beneath the deep crustal 35 km jabalpur earthquake of 21 may 1997 narmada son lineament nsl central india in search of a possible cause of stress accumulation in the region teleseismic receiver functions computed at nine digital seismographs along a 250 km long profile suggests moho downwarp to 52 km across the width of the lineament in contrast with an average 40 km depth elsewhere in addition the crust beneath the nsl has higher vp vs of 1 84 compared to 1 73 in the surrounding suggestive of a high density mafic mass at depth that compensates the crustal root also supported by small topographic variation 200 m across the lineament presence of such an anomalous mass in the deep crust may lead to gravitationally induced stresses in the lower crust that contribute to the failure of rock along the pre existing narmada son fault leading to the earthquake in deep crust

shaping mobile belts by small scale convection claudio faccenna and thorsten becker compute global mantle flow based on a recent high resolution seismic tomography model to investigate the role of buoyancy driven and plate motion induced mantle circulation for the mediterranean they show that mantle flow provides an explanation for much of the observed dynamic topography and microplate motion in the region and conclude that vigorous small scale convection in the uppermost mantle may also hold the key for understanding other complex mobile belts such as those in the north american cordillera and the himalayan tibetan collision zone here global mantle flow is computed on the basis of recent high resolution seismic tomography to investigate the role of buoyancy driven and plate motion induced mantle circulation in the mediterranean region the findings show that mantle flow explains much of the observed dynamic topography and microplate motion in the region small scale convection in the uppermost mantle may also hold the key to understanding complex mobile belts elsewhere mobile belts are long lived deformation zones composed of an ensemble of crustal fragments distributed over hundreds of kilometres inside continental convergent margins1 2 the mediterranean represents a remarkable example of this tectonic setting3 the region hosts a diffuse boundary between the nubia and eurasia plates comprised of a mosaic of microplates that move and deform independently from the overall plate convergence4 surface expressions of mediterranean tectonics include deep subsiding backarc basins intraplate plateaux and uplifting orogenic belts although the kinematics of the area are now fairly well defined the dynamical origins of many of these active features are controversial and usually attributed to crustal and lithospheric interactions however the effects of mantle convection well established for continental interiors5 6 7 should be particularly relevant in a mobile belt and modelling may constrain important parameters such as slab coherence and lithospheric strength here we compute global mantle flow on the basis of recent high resolution seismic tomography to investigate the role of buoyancy driven and plate motion induced mantle circulation for the mediterranean we show that mantle flow provides an explanation for much of the observed dynamic topography and microplate motion in the region more generally vigorous small scale convection in the uppermost mantle may also underpin other complex mobile belts such as the north american cordillera or the himalayan tibetan collision zone

active deformation of asia from kinematics to dynamics using estimated strain rates across the actively deforming region of asia to infer variations in viscous stress it is shown that gradients of stress point uphill toward the center of tibet where gravitational potential energy per unit surface area reaches a maximum thus the dynamics of the deformation seem to obey the equation of creeping flow which expresses a balance between gradients in stress and the gravitational body force this balance in the region of the tibetan plateau yields an estimate of 10 22 pascal second for the average viscosity of the tibetan lithosphere which is only about 10 to 100 times greater than the viscosity of the convecting upper mantle

predicting and testing continental vertical motion histories since the paleozoic abstract dynamic topography at the earth s surface caused by mantle convection can affect a range of geophysical and geological observations including bathymetry sea level change continental flooding sedimentation and erosion these observations provide important constraints on and test of mantle dynamic models based on global mantle convection models coupled with the surface plate motion history we compute dynamic topography and its history for the last 400 ma associated with pangea assembly and breakup with particular focus on cratonic regions we propose that burial unroofing histories of cratons inferred from thermochronology data can be used as a new diagnostic to test dynamic topography and mantle dynamic models our models show that there are currently two broad dynamic topography highs in the pacific and africa for the present day earth that are associated with the broad warm structures i e superplumes in the deep mantle consistent with previous proposals of dynamical support for the pacific and african superswells our models reveal that pangea assembly and breakup by affecting subduction and mantle upwelling processes have significant effects on continental vertical motions our models predict that the slave craton in north america subsides before pangea assembly at 330 ma but uplifts significantly from 330 ma to 240 ma in response to pre pangea subduction and post assembly mantle warming the kaapvaal craton of africa is predicted to undergo uplift from 180 ma to 90 ma after pangea breakup but its dynamic topography remains stable for the last 90 ma the predicted histories of elevation change for the slave and kaapvaal cratons compare well with the burial unroofing histories inferred from thermochronology studies thus supporting our dynamic models including the development of the african superplume mantle structure the vertical motion histories for other cratons can provide further tests of and constraints on our mantle dynamic models

origins of topography in the western u s mapping crustal and upper mantle density variations using a uniform seismic velocity model to investigate the physical basis for support of topography in the western u s we construct a subcontinent scale 3 d density model using 1000 estimated crustal thicknesses and s velocity profiles to 150 km depth at each of 947 seismic stations crustal temperature and composition are considered but we assume that mantle velocity variations are thermal in origin from these densities we calculate crustal and mantle topographic contributions typical 2 uncertainty of topography is 500 m and elevations in 84 of the region are reproduced within error remaining deviations from observed elevations are attributed to melt variations in crustal quartz content and dynamic topography compositional variations in the mantle while plausible are not necessary to reproduce topography support for western u s topography is heterogeneous with each province having a unique combination of mechanisms topography due to mantle buoyancy is nearly constant within 250 m across the cordillera relief there 2 km results from variations in crustal chemistry and thickness cold mantle provides 1 5 km of ballast to the thick crust of the great plains and wyoming craton crustal temperature variations and dynamic pressures have smaller magnitude and or more localized impacts positive gravitational potential energy gpe anomalies 2 1012n m calculated from our model promote extension in the northern basin and range and near the sierra nevada negative gpe anomalies 3 1012n m along the western north american margin and yakima fold and thrust belt add compressive stresses stresses derived from lithospheric density variations may strongly modulate tectonic stresses in the western u s continental interior

a tectonic redefinition of the southern rocky mountains abstract the southern rocky mountains of the western united states physiographically defined and described by n m fenneman nearly 60 yrs ago are tectonically redefined and extended they are shown to constitute the crestal range of a mammoth continental arch like feature here named the alvarado ridge its axis trends south from casper wyoming at least as far south as el paso texas beyond which the ridge begins to lose morphological identity maximum elevations along the crest of the ridge exceed 4 2 km the summit ranges are bordered on either side by gently sloping rises that extend outward for at least 1300 km falling to elevations of less than 400 m modest rift structures along the ridge axis continue beyond the southern terminus of the mountains before playing out a major sediment filled axial graben exists over the southern two thirds of the ridge but equivalent parts of it farther north were stripped of their neogene fill by erosion in the headwaters areas of the colorado and north platte rivers the maximum elevation of earlier laramide mountains in this area has been estimated to have been no more than 2 km in colorado half that of the present range related topography was nearly obliterated by erosion prior to late eocene time lateral stream planation produced a southeast sloping major late eocene erosion surface across the region that had what was probably an isostatically adjusted average maximum elevation of less than 900 m present day elevations and relief on the ridge crest are the result of steep crestal normal faulting pronounced block uplift and regional arching with extensional strain limited to an axial corridor less than 200 km wide this episode of mountain building began in middle miocene time 17 to 12 ma culminating in latest miocene and early pliocene time between 7 and 4 ma ago debris from the newly elevated range the southern rocky mountains sensu lato was shed along the full length of the neogene alvarado ridge down parts of both rises as well as into its medial graben the uplifted and tilted east rise the great plains ramp has an extensive aggradational cover that spread eastward to distances of 700 km and more in middle miocene to earliest pliocene time its fluviatile lacustrine and eolian sediments are undergoing dissection today but a sufficient portion remains to define the unmodified profile of the east rise as well as the original length of the ridge it is coincidentally similar in configuration and scale to the west rise of the mid atlantic ridge the topography of the crestal mountain ranges on both structures is also remarkably similar the southern rocky mountains are as high as the alps but are not the product of lithospheric compression and crustal thickening rather they reflect profound epeirogenic uplift driven by lithospheric thinning and accompanied by extensional strain and differential vertical jostling along the crest of the ridge

upper mantle shear structure beneath the colorado rocky mountains a tomographic inversion of teleseismic shear wave delays has been performed using data collected in the 1992 program for array seismic studies of the continental lithosphere rocky mountain front experiment the shear wave residuals used as data were corrected for known crustal sediment and topographic differences and thus were sensitive to mantle shear velocity variations beneath the colorado array large differences in shear wave travel time are observed across colorado and into western kansas with a maximum difference of near 5 sec between the central rockies in colorado and a station in western kansas the shear wave residuals were inverted for shear wave velocity variations in the upper 350 km of the mantle beneath colorado and kansas large differences in mantle shear velocity are found between the rockies and the kansas great plains the average shear velocity in the upper 200 km is about 9 greater beneath kansas than beneath the rockies at greater depths lateral variations are far smaller and not as well correlated with surface tectonics in a qualitative sense there appear to be three mantle structures present beneath the study area beneath the rockies the upper mantle is extremely slow both to the west and east of the rockies we find intermediate mantle shear velocity there is then an abrupt increase in mantle shear velocity near the colorado kansas border we interpret the abrupt increase in mantle shear velocity near the colorado kansas border to mark the boundary of cratonic lithosphere the slow velocity zone beneath the rockies appears to be upwelled asthenosphere as is seen beneath the rio grande rift to the south to match the seismic contrast between the mantle beneath kansas and the colorado rockies requires a temperature contrast of at least 350 c as well as partial melt beneath the rockies the hot mantle beneath the rockies implies that the mantle density is low and the high topography of the rockies is at least partially due to buoyant mantle

mantle dynamics isostasy and the support of high terrain dynamic topography is commonly understood to be deflection of the earth s surface that results from convection of the mantle because different authors use the words dynamic topography differently topography designated as dynamic may amount to only a few hundred meters or may exceed 2000 m for most regions however surface heights computed on the assumption that the lithospheric column is in isostatic equilibrium provide good approximations to observed topography the small free air gravity anomalies and still smaller isostatic anomalies 2000 m of dynamic topography moreover simple examples show that if flow induced stresses maintain a given density distribution then the resultant surface deflections can be smaller than those that would exist in isostatic equilibrium in light of confusion over the term dynamic topography we offer readers simple tools that we hope will enable them to diagnose the component of topography that results from flow induced stresses and to distinguish it from topography that is compensated isostatically

tectonics of pliocene removal of lithosphere of the sierra nevada california pliocene ca 3 5 ma removal of dense eclogitic material under the sierra nevada has been proposed from variations in the petrology and geochemistry of neogene volcanic rocks and their entrained xenoliths from the southern sierra the replacement of eclogite by buoyant warm asthenosphere is consistent with present day seismologic and magnetotelluric observations made in the southern sierra a necessary consequence of replacing eclogite with peridotite is that mean surface elevations and gravitational potential energy both increase an increase in potential energy should increase extensional strain rates in the area if these forces are insuffi cient to signifi cantly alter pacifi c north american plate motion then increased extensional strain rates in the vicinity of the sierra must be accompanied by changes in the rate and style of deformation elsewhere changes in deformation in california and westernmost nevada agree well with these predictions existing geologic evidence indicates that a period of rapid uplift along the sierran crest of more than 1 km occurred between 8 and 3 ma most likely as a consequence of removal of lower lithosphere about this same time extensional deformation was initiated within 50 km of the eastern side of the sierra 5 3 ma and regional shortening began to produce the california coast ranges 5 3 ma we suggest that these events were induced by the 1 2 10 12 n m increase of gravitational potential energy generated by the sierran uplift evidence for pliocene uplift adjoining crustal extension and shortening in directly opposing parts of the coast ranges is found along nearly the entire length of the sierra nevada and implies that lithosphere was removed beneath all of the presentday mountain range the uplifted area lies between two large upper mantle high pwave velocity bodies under the south end of the san joaquin valley and the north end of the sacramento valley these high velocity bodies plausibly represent the present position of material removed from the base of the crust lithospheric removal may also be responsible for shifting of the distribution of transform slip from the san andreas fault system to the eastern california shear zone a prediction that awaits better defi ned slip histories on both faults overall the late cenozoic deformational history of the sierra nevada and vicinity illustrates that locally derived forces can infl uence deformation kinematics within plate boundary zones

northward thinning of tibetan crust revealed by virtual seismic profiles 1 a new approach of constructing deep penetrating seismic profiles reveals significant regional variations in crustal thickness under near constant elevation of tibet over distances of hundreds of kilometers the crust is as thick as 75 km in southern tibet but shoals to just over 60 km under the qiangtang terrane in central tibet where the deviation from airy isostasy is equivalent to a thickness of over 10 km in missing crust northward thinning of crust occurs gradually over a distance of about 200 km where mechanical deformation instead of pervasive magmatism also seems to have disrupted the crust mantle interface

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

steady propagation of delamination events delamination of the lithospheric thermal boundary from overlying continental crust propagates laterally from the line of initiation accelerating as the sinking slab of detached lithosphere grows longer this propagation has been numerically modeled with steady state equations in a moving reference frame by matching an interior finite element solution to flexible boundary conditions which represent the mechanical and thermal response of the surroundings the form of the solution depends on the shear coupling of intruding asthenosphere to the top of the sinking slab across a thin layer of crustal material without coupling the tip of the intrusion cools and stiffens to form a wedge dividing the crust cold mode with coupling the intrusion is forced to convect and remains ductile hot mode the cold mode can propagate at all velocities the hot mode has a lower limiting velocity of 1 2 cm year but offers less resistance at higher speeds resistance to delamination includes a constant term from the buoyant crustal downwarp plus a velocity proportional term representing viscous deformation however the proportionality constant of the latter term is only weakly dependent on crust and lithosphere viscosities matching this resistance to loading lines of 100 to 800 km slabs sinking in a mantle of 1022 p velocities of 0 3 8 0 cm year are obtained changes in viscosity affect this rate but cold mode delamination is unstoppable except at continental margins or by failure in the sinking slab the surface expression of delamination is a leading outer rise followed by a submarine trough with a large negative free air anomaly which finally evolves into a 1 km plateau if crustal viscosity and velocity are both low however there is a montonic crustal uplift with no trough thus the present lack of linear supracontinental oceans does not preclude delamination at up to 4 cm year driven by slabs up to 400 km in length

genesis of collision volcanism in eastern anatolia turkey late cenozoic volcanism in eastern anatolia extends in a broad sw ne trending belt across the arabia eurasia collision zone from the arabian foreland basin in the southwest and to the kars plateau and lesser caucasus in the northeast foreland volcanism is dominated by basaltic shield and fissure eruptions of transitional tholeiitic alkaline composition volcanism on the thickened crust north of the bitlis thrust zone varies from the mildly alkaline volcano nemrut and older mus volcanics in the south through the transitional calc alkaline alkaline volcanoes bingol and suphan and the alkaline volcano tendurek to the calc alkaline volcano ararat and older kars plateau volcanics in the north isotope sr nd and trace element systematics indicate that the lavas from the foreland were derived from the mantle lithosphere of the arabian continent which had been enriched by small volumes of asthenospheric melts over a period of time and that lavas from the alkaline volcanic area around mus and the volcanoes nemrut and tendurek north of the bitlis thrust zone were derived from a lithospheric source of similar composition either from the same underthrust arabian continent or from the bitlis massif microcontinent by contrast the transitional lavas from bingol and suphan and the calc alkaline lavas from ararat and kars were derived from lithosphere carrying a distinct subduction signature inherited from pre collision subduction events positive correlations between 87sr 86sr and sio2 and rb nb and sio2 in the alkaline and transitional lavas suggest that combined assimilation and fractional crystallization was an important process within at least part of the thickened crust of the collision zone trace element covariation diagrams such as y rb indicate the importance of hornblende crystallization at depth and orthopyroxene at shallow levels within the calc alkaline provinces in contrast to the consistently anhydrous crystallization sequences of the alkaline lavas trace element diagrams based on the covariation of compatible and incompatible elements point to moderate low degrees of partial melting with residual clinopyroxene throughout and residual garnet in the foreland province consideration of mineral stabilities mantle solidi and geothermal gradients before and after collision suggest that lithospheric thickening should both increase the thickness of metasomatized lithosphere and depress the metasomatized zone to greater depths probably beneath the amphibole and dolomite breakdown curves perturbation of the thickened lithosphere by delamination of the thermal boundary layer perhaps coupled with local stretching associated with pull apart basins on strike slip fault systems is then sufficient to generate melt the composition of that melt being largely dependent on the enrichment history of the lithosphere in question in eastern anatolia volcanism appears to have started between about 8 and 6 ma ago some 5 ma after the start of rapid uplift of the east anatolian plateau

the effective elastic thickness of the india plate from receiver function imaging gravity anomalies and thermomechanical modelling summary the range and the meaning of the effective elastic thickness eet in continental areas have been subject to controversy over the last two decades here we take advantage of the new data set from the hi climb seismological experiment to re estimate the eet of the india plate along a south north profile extending from the ganges basin to central tibet receiver functions give a high resolution image of the base of the foreland basin at 5 km depth and constrain the crustal thickness which increases northwards from 35 km beneath the indo gangetic plain to 70 km in southern tibet together with available data sets including seismic profiles seismological images from both indepth and himnt experiments deep well measurements and bouguer anomaly profiles we interpret this new image with 2 d thermomechanical modelling solutions using different type of crustal and mantle rheologies we find that 1 the eet of the india plate decreases northwards from 60 80 to 20 30 km as it is flexed down beneath himalaya and tibet due to thermal and flexural weakening 2 the only resistant layer of the india plate beneath southern tibet is the upper mantle which serves as a support for the topographic load and 3 the most abrupt drop in the eet located around 200 km south of the mft is associated with a gradual decoupling between the crust and the mantle we show that our geometrical constraints do not allow to determine if the upper and lower crust are coupled or not our results clearly reveal that a rheology with a weak mantle is unable to explain the geometry of the lithosphere in this region and they are in favour of a rheology in which the mantle is strong

gravitational potential energy and active deformation in the apennines we use velocity measurements from a network of continuous gps sites spanning the apennines of peninsular italy to test the hypothesis that the active deformation of the region is explained by variations in gravitational potential energy of the lithosphere the simple geometry of the mountain chain allows us to treat the deformation as two dimensional neglecting gradients of velocity along the strike of the chain under this assumption the integral of gravitational potential energy per unit area of the lithosphere gpe in the direction perpendicular to the chain is related by a simple expression to the velocity in the same direction we show that the observed velocities match this expression with an rms misfit of 0 5mm yr this agreement suggests that deformation of the apennines reflects a balance within the mountain chain itself between lateral variations in gpe and the stresses required to deform the lithosphere forces arising from processes external to the belt are not required to explain the observations

alkylated phenanthrene distribution in artificially matured kerogens from kimmeridge clay and the brent formation north sea abstract kerogen and extractible organic matter from immature kimmeridge clay and brent coal were heated at different temperatures 270 500 c under conditions to mimic natural maturation aromatic fractions have been examined for their phenanthrene compound distribution phenanthrene indices based on methylphenanthrenes mp and dimethylphenanthrenes dmp exhibit very sensible changes in the applied temperature range however at high temperatures about 450 c and more possible thermal degradation of phenanthrene compounds produces a decrease in the values of the indices experimental data obtained in this study are compared and discussed with data obtained on thermal evolution performed on total rocks vitrinite reflectance values measured on artificially matured coals are also tentatively correlated with phenanthrene ratios theoretical thermodynamic calculations on mp compounds confirm also the relative thermal stability of the individual isomers

rapid erosion of the southern african plateau as it climbs over a mantle superswell we present new sedimentary flux data confirming that a large pulse of erosion affected the southern african plateau in the late cretaceous and is likely to be related to a major uplift episode of the plateau this short phase of erosion i e less than 30 myr in duration has commonly been difficult to reconcile with a mantle origin for the plateau anomalous uplift given its size the rise of the african superplume is likely to have lasted much longer here we demonstrate by using a simple model for fluvial erosion that tilting of the continent as it rides over a wide dynamic topography high cannot only cause rapid uplift of the plateau but also trigger continent wide drainage reorganization leading to substantial denudation in a relatively short amount of time the amplitude and short duration of the sedimentary pulse are best reproduced by assuming a strong erodibility contrast between the karoo sedimentary and volcanic rocks and the underlying basement we also present a new compilation of paleoclimate indicators that shows a transition from arid to very humid conditions approximately at the onset of the documented erosional pulse suggesting that climate may have also played a role in triggering the denudation the diachronism of the sedimentary flux between the eastern and western margins of the plateau and the temporal and geographic coincidence between the uplift and kimberlite eruptions are however better explained by our tilt hypothesis driven by the migration of the continent over a fixed source of mantle upwelling

dynamic topography produced by lower crustal flow against rheological strength heterogeneities bordering the tibetan plateau summary dynamic stresses developed in the deep crust as a consequence of flow of weak lower crust may explain anomalously high topography and extensional structures localized along orogenic plateau margins with lubrication equations commonly used to describe viscous flow in a thin gap geometry we model dynamic stresses associated with the obstruction of lower crustal channel flow due to rheological heterogeneity dynamic stresses depend on the mean velocity u viscosity and channel thickness h uniquely through the term u h 2 these stresses are then applied to the base of an elastic upper crust and the deflection of the elastic layer is computed to yield the predicted dynamic topography we compare model calculations with observed topography of the eastern tibetan plateau margin where we interpret channel flow of the deep crust to be inhibited by the rigid sichuan basin model results suggest that as much 1500 m of dynamic topography across a region of several tens to a hundred kilometres wide may be produced for lower crustal material with a viscosity of 2 10 18 p asfl owing in a 15 km thick channel around a rigid cylindrical block at an average rate of 80 mm yr 1

the relationship between seismic velocity and density in the continental crust a useful constraint summary plots of seismic velocity and density of rock samples show that a range of densities is possible for rocks of each seismic velocity and vice versa although a single linear relationship is often assumed in crustal gravity calculations because of the scatter whenever rocks of known seismic velocity are converted to density using this relationship a reduction is made to the resolving power of the resulting gravity calculation if these rocks reach thicknesses of more than a few kilometres then the uncertainties become significant when compared with the size of commonly observed gravity anomalies examples are considered from the north sea mississippi and carolina trough it is concluded that the use of a seismic velocity measurement as the only indication of rock density does not provide a useful constraint when attempting to reproduce observed gravity variations an appropriate model for isostatic compensation is probably the most important factor for successful predictions of crustal structure on the basis of gravity data

mantle dynamics uplift of the tibetan plateau and the indian monsoon convective removal of lower lithosphere beneath the tibetan plateau can account for a rapid increase in the mean elevation of the tibetan plateau of 1000 m or more in a few million years such uplift seems to be required by abrupt tectonic and environmental changes in asia and the indian ocean in late cenozoic time the composition of basaltic volcanism in northern tibet which apparently began at about 13 ma implies melting of lithosphere not asthenosphere the most plausible mechanism for rapid heat transfer to the midlithosphere is by convective removal of deeper lithosphere and its replacement by hotter asthenosphere the initiation of normal faulting in tibet at about 8 3 ma suggests that the plateau underwent an appreciable increase in elevation at that time an increase due solely to the isostatic response to crustal thickening caused by india s penetration into eurasia should have been slow and could not have triggered normal faulting another process such as removal of relatively cold dense lower lithosphere must have caused a supplemental uplift of the surface folding and faulting of the indo australian plate south of india the most prominent oceanic intraplate deformation on earth began between about 7 5 and 8 ma and indicates an increased north south compressional stress within the indo australian plate a tibetan uplift of only 1000 m if the result of removal of lower lithosphere should have increased the compressional stress that the plateau applies to india and that resists india s northward movement from an amount too small to fold oceanic lithosphere to one sufficient to do so the climate of the equatorial indian ocean and southern asia changed at about 6 9 ma monsoonal winds apparently strengthened northern pakistan became more arid but weathering of rock in the eastern himalaya apparently increased because of its high altitude and lateral extent the tibetan plateau provides a heat source at midlatitudes that should oppose classical symmetric hadley circulation between the equator and temperate latitudes and that should help to drive an essentially opposite circulation characteristic of summer monsoons for the simple case of axisymmetric heating no dependence on longitude of an atmosphere without dissipation theoretical analyses by hou lindzen and plumb show that an axisymmetric heat source displaced from the equator can drive a much stronger meridianal monsoonlike circulation than such a source centered on the equator but only if heating exceeds a threshold whose level increases with the latitude of the heat source because heating of the atmosphere over tibet should increase monotonically with elevation of the plateau a modest uplift 1000 2500 m of tibet already of substantial extent and height might have been sufficient to exceed a threshold necessary for a strong monsoon the virtual simultaneity of these phenomena suggests that uplift was rapid approximately 1000 m to 2500 m in a few million years moreover nearly simultaneously with the late miocene strengthening of the monsoon the calcite compensation depth in the oceans dropped plants using the relatively efficient c4 pathway for photosynthesis evolved rapidly and atmospheric co2 seems to have decreased suggesting causal relationships and positive feedbacks among these phenomena both a supplemental uplift of the himalaya the southern edge of tibet and a strengthened monsoon may have accelerated erosion and weathering of silicate rock in the himalaya that in turn enhanced extraction of co2 from the atmosphere thus these correlations offer some support for links between plateau uplift a downdrawing of co2 from the atmosphere and global climate change as proposed by raymo ruddiman and froehlich mantle dynamics beneath mountain belts not only may profoundly affect tectonic processes near and far from the belts but might also play an important role in altering regional and global climates

the many surface expressions of mantle dynamics plate tectonic theory suggests that present day topography can be explained by the repeated interactions between the tectonic plates moving along earth s surface however mounting evidence indicates that a significant proportion of earth s topography results from the viscous stresses created by flow within the underlying mantle rather than by the moving plates this dynamic topography is transient varying as mantle flow changes and is characterized by small amplitudes and long wavelengths it is therefore often hidden by or confused with the more obvious topographic anomalies resulting from horizontal tectonic movements however dynamic topography can influence surface processes and thus enter the geological record it has for example played a role in the establishment of amazon drainage patterns in turn surface processes such as the erosion of topographical anomalies could affect mantle flow this emerging view of dynamic topography suggests that the concept of plate tectonics as the driver of surface deformation needs to be extended to include the vertical coupling between the mantle and the surface unravelling this coupling back in time with the help of models and the geological record can potentially provide unprecedented insights into past mantle dynamics

rayleigh taylor instability lithospheric dynamics surface topography at convergent mountain belts and gravity anomalies 1 surface topography and associated gravity anomalies above a layer resembling continental lithosphere whose mantle part is gravitationally unstable depend strongly on the ratio of viscosities of the lower density crustal part to that of the mantle part for linear stability analysis growth rates of rayleigh taylor instabilities depend largely on the wave number or wavelength of the perturbation to the base of the lithosphere and weakly on this viscosity ratio on plausible density differences among crust mantle lithosphere and asthenosphere and on ratios of crustal to total lithospheric thicknesses for all likely densities viscosities and thicknesses the moho is drawn down pushed up where the base of the lithosphere subsides rises for large viscosities of crust compared to mantle lithosphere ratios 30 a sinking and thickening mantle lithosphere also pulls the surface down for smaller viscosity ratios crustal thickening overwhelms the descent of the moho and the surface rises subsides above regions where mantle lithosphere thickens and descends thins and rises ignoring vertical variations of viscosity within the crust and mantle lithosphere we find that the maximum surface height occurs for approximately equal viscosities of crust and mantle lithosphere for large crust mantle lithosphere viscosity ratios gravity anomalies follow those of surface topography with negative positive free air anomalies over regions of descent ascent in this case topography anomalies are smaller than those that would occur if the lithosphere were in isostatic equilibrium hence flow induced stresses dynamic pressure and deviatoric stress create smaller topography than that expected for an isostatic state for small crust mantle viscosity ratios 10 however calculated surface topography at long wavelengths is greater than it would be if the lithospheric column were in isostatic equilibrium and at short wavelengths local isostasy predicts surface deflections of the wrong sign for the range of wavelengths appropriate for convergent mountain belts 150 600 km calculated gravity anomalies are negative over regions of lithospheric thickening especially when allowance for flexural rigidity of a surface layer is included correspondingly calculated values of admittance the ratio of fourier transforms of surface topography and free air gravity anomalies are also negative for wave numbers relevant to mountain belts for essentially all mountain belts however measured free air anomalies and admittance are positive whether gravitational instability of the lithosphere affects the structure of convergent belts or not its contribution to the topography of mountain belts seems to be small compared to that predicted for isostatic balance of crustal thickness variations

constraints on the structure of the himalaya from an analysis of gravity anomalies and a flexural model of the lithosphere the intracontinental subduction of india beneath the himalaya presents several similarities to that occurring at island arcs we study one of those similarities by analyzing gravity anomalies across the himalaya assuming that the topography is supported by the indian elastic plate flexed under the weight of both the overthrust mountains and the sediments in the ganga basin we first examine in detail the effects of each of the following parameters on the configuration of the elastic plate and on the gravity anomalies the flexural rigidity the position of the northern end of the elastic plate the amount of underthrusting of such a plate beneath the range and the density contrasts between the crust and mantle and between the sediments and the crust a plate with a constant flexural rigidity of about 0 7 1025 n m between 0 2 and 2 0 1025 n m allows a good fit to the data from the lesser himalaya and the ganga basin such a plate however cannot underthrust the entire himalaya instead the gravity anomalies show that the moho steepens from only about 3 beneath the lesser himalaya to about 15 beneath the greater himalaya this implies a smaller flexural rigidity beneath the greater himalaya 0 1 to 1 0 1023 n m than beneath the ganga basin and the lesser himalaya even with a thin weak plate beneath the greater himalaya the weight of the mountains depresses the plate too much unless an additional force or moment is applied to the plate the application of a bending moment unit length to the end of the plate of about 0 6 1018 n m is adequate to elevate the indian plate and to bring the calculated gravity anomalies in agreement with those observed both the smaller flexural rigidity and the bending moment can be understood if we assume that part or all of the indian crust has been detached from the lower lithosphere that underthrusts the greater himalaya we study the tectonic implications of these results by means of a series of idealized balanced cross sections from the collision to the present that reproduce several important features of the geology of the himalaya and predict an amount of eroded material comparable to that in the ganga basin and the bay of bengal these cross sections include high grade metamorphic rocks near the main central thrust and a steeper dip of it there than in the lesser himalaya they predict rapid uplift only in the greater himalaya and at the foot of the lesser himalaya

the isostatic stagnant lid approximation and global variations in the venusian lithospheric thickness we use the long wavelength venusian topography and geoid to constrain the global lithospheric thickness by assuming that the stagnant lid overlying the convective mantle is in isostatic equilibrium the isostatic stagnant lid isl approximation two dimensional and three dimensional numerical simulations of convection show that the isl approximation is reasonably accurate in the tested parameter range of temperature dependant viscosity convection with bottom heating and internal heating the errors in predicting the thickness of the stagnant lid i e the lithosphere are somewhat larger for internally heated convection we attribute this to a flatter lid and as a result to a larger relative contribution to topography and geoid anomalies from the deeper parts of the layer application of the isl approximation to venus suggests that the average lithospheric thickness might be as large as 600 km which is on the high end of previously proposed estimates

gravity anomalies and flexure of the lithosphere at mountain ranges the bouguer gravity anomaly over the himalayan alpine and appalachian mountains is characterized by a generally asymmetric gravity low which spans the mountains and associated foreland basins the minimum of the gravity low is generally systematically displaced from the region of greatest topographic relief and shows no obvious relationship to surface geology in addition the alps and appalachians are associated with a generally symmetric gravity high that is unrelated to the topographic relief together the gravity low and high form a characteristic positive negative anomaly couple the steep gravity gradient between the positive and negative couple correlates with the insubric line and its equivalents in the alps and the brevard zone in the southern appalachians this gravity anomaly couple is interpreted as evidence for flexure of the continental lithosphere by subsurface buried and surface topographic loading the magnitude of the subsurface load is estimated from the integrated excess mass represented by the positive anomaly and the magnitude of the surface load from the topography by combining these loads the flexure of the lithosphere and the associated bouguer gravity anomaly are calculated for different assumed values of the elastic thickness te of continental lithosphere the best fitting value of te for the himalayas is 80 100 km for the alps 25 50 km and for the appalachians 80 130 km the model calculations suggest that the main contribution to the gravity couple in the alps and appalachians arises not from surface loads such as thrust sheets and nappes but from subsurface loads however in the case of the himalayas the topography appears to be a sufficient load to explain the observed asymmetric gravity low when subsurface loads do exist therefore we would not expect a close correlation between mountain topographic relief and depth to the m discontinuity as required by classical models of isostasy rather the m discontinuity is expected to relate to the center of mass of all the loads acting on the lithosphere the surface and subsurface loads inferred from the gravity anomaly play a major role in the development of mountains and foreland basins the emplacement of large subsurface loads represents the primary event in the alps and appalachians while in the himalayas the emplacement of the surface load is the primary event the emplacement of the subsurface load probably marks the initiation of foreland basin development the regional characteristics of the basin are controlled by the primary load while local variations are controlled by the secondary thrust sheet loads subsequent movement of either load will produce migration of the basin depocenter the origin of the subsurface load is not clear but its association with the insubric line in the alps and the brevard zone in the appalachians suggests that it may be related to the obduction of crustal blocks flakes onto the underlying plate during continental collision or possibly to the preloaded crustal structure of the underlying plate

phanerozoic marine inundation of continents driven by dynamic topography above subducting slabs a spherical model of mantle flow constrained by the locations of trenches can be used to predict the dynamic topography of the earth s surface and hence the marine inundation of continents for past periods of high sea level the predicted geographical pattern of flooding correlates well with the geological record the high spatial correlation may result from increased plate velocities at these times leading to increased rates of subduction subsidence and inundation at convergent margins

mantle convection and the recent evolution of the colorado plateau and the rio grande rift valley the colorado plateau contains late cretaceous marine strata that are at a mean elevation of 2 km the timing and amount of uplift since the cretaceous has generated considerable debate with the exception of a few studies topography supported by vertical stresses generated by viscous fl ow in the mantle has not been explicitly considered to contribute to the elevation of this region herein we compute the viscous fl ow beneath north america that is driven by density anomalies inferred from joint seismic geodynamic modeling we fithat the colorado plateau overlies a strong mantle upwelling that is coupled to the sinking farallon slab currently beneath the eastern united states consequently the colorado plateau is currently a focused dynamic topography high within the western u s cordillera moreover this strong upwelling impacts the base of the lithosphere at an oblique angle east of the plateau directly below the rio grande rift we attribute this fl ow as being responsible for some of the recent magmatic activity along the jemez lineament as well as contributing to the recent rifting process in the rio grande rift valley

gravity anomalies flexure of the indian plate and the structure support and evolution of the himalaya and ganga basin bouguer gravity anomalies along four profiles across the western himalaya and ganga basin show large deviations from local isostatic equilibrium a deficit of mass characterizes the ganga basin and an excess underlies the lesser himalaya both can be understood if the indian plate is flexed down by the distributed load of part of the mountains the cross sectional shape of the ganga basin seems to be controlled by the deflection of the indian plate which we compute assuming the indian plate to overlie an inviscid fluid from the shapes of both bouguer anomaly profiles and the basement topography we place bounds on the flexural rigidity of such a plate if the ganga basin is a steady state feature then the age of the basal sediments in a given locality should be proportional to the distance of that locality from the southern edge of the basin if the rate of convergence of india and the himalaya were constant that rate should equal the distance divided by the corresponding age we find a rate of 10 to 15 mm a for the last 15 to 20 ma which is consistent with a large part of the 50 mm a rate of convergence between india and eurasia being absorbed by the eastward extrusion of parts of tibet profiles of bouguer gravity anomalies show only a small peak or plateau over the southern edge of the lesser himalaya implying that the boundary between the light sediments of the ganga basin and the heavier crustal rocks of the lesser himalaya is not sharp and that there exists some light material beneath the range we infer that some sediment deposited in the ganga basin has been underthrust beneath the lesser himalaya but the quantity is small most of this sediment probably is scraped off the indian plate to make the foothills of the range we find that the load of the high himalaya is too large to be supported solely by elastic stress in the indian plate if the flexural rigidity of the plate is constant and if no other external forces act on the plate the observed gradient in bouguer gravity anomalies increases from an average of about 1 mgal km over the ganga basin and lesser himalaya to about 2 mgal km over the high himalaya this increase in the gravity gradient implies that the moho dips more steeply 10 15 beneath the high himalaya than beneath the lesser himalaya 2 3 we interpret this steepening of the moho to be due to a weakening of the plate which allows it to bend more sharply beneath the high himalaya than farther south with the inclusion of a weak segment of indian plate beneath the high himalaya calculated anomalies show a somewhat increased gradient beneath the high himalaya but when the weight of the entire himalaya is used as a load on the plate calculated anomalies are more negative than observed therefore an external force system is needed to support much of the weight of the high himalaya as well as to bend the plate sufficiently beneath the high himalaya the magnitudes of the bending moment and the force per unit length that must be applied to the end of the plate are compatible with their sources being gravity acting on a part of india s mantle lithosphere stripped of its crust and underthrust beneath southern tibet

high resolution surface wave tomography beneath the aegean anatolia region constraints on upper mantle structure summary this study provides new constraints on the upper mantle structure from western greece to central anatolia using seismic data of permanent broad band networks recently installed in greece and turkey and from a two year temporary array simbaad experiment we used 200 seismic events recorded at 146 broad band stations with a typical interstation distance of 60 100 km across the study area the high resolution 3 d shear wave velocity model of the mantle is obtained by inversion of fundamental mode rayleigh wave phase velocity maps for periods between 20 and 195 s the tomography is based on ray tracing in heterogeneous media taking into account external propagation effects the horizontal resolution is approximately 100 km however small heterogeneities may suffer from some horizontal smearing and damping the vertical resolution is approximately 100 km the vertical smoothing is necessary to avoid unresolved spurious shear wave velocity oscillations in the upper mantle the errors on shear wave velocities in our 3 d model 0 02 0 1 km s 1 are significantly smaller than the amplitude of vs variations 0 3 0 5 km s 1 in spite of the vertical and horizontal smoothing our model shows details in the upper mantle structure never reached at regional scale in the area the overall structure is characterized by a low velocity zone 80 200 km depth reflecting a slow and warm asthenosphere underlying a thin lithosphere the southwesternmost termination of the low velocity anomaly corresponds to the northward dipping hellenic slab the detailed shear velocity structure of the upper mantle beneath anatolia appears to be far more geometrically complex than revealed in previous tomographic studies of the area at depths larger than or equal to 160 km velocities are overall high beneath anatolia partly due to the presence of dipping high velocity anomalies which we tentatively interpret as remnant slabs the southernmost high velocity anomaly beneath anatolia is separated from the eastern edge of the hellenic slab by a major low velocity anomaly which we interpret as the trace of asthenospheric mantle material rising inside a vertical slab tear beneath southwestern anatolia

anomalous uplift of the apennines and subsidence of the adriatic the result of active mantle flow 1 we consider the potential for interpreted mantle flow to generate anomalous topography across the north central mediterranean a seismic tomography section was translated into the temperature domain and used as input into a series of 2d thermo mechanical numerical models of crust mantle deformation the generated dynamic topography signal was then compared with a derived residual topography that showed a distinct residual low across the adriatic sea and easternmost apennines and a residual high to the west of the apennine range and into the tyrrhenian sea a good spatial correlation was found between the modeled dynamic effects of inferred underlying mantle structure and the observed residual topography the results of this study suggest that underlying mantle structure may play a role in supporting the dynamic anomalous surface topography along the investigated profile and within the north central mediterranean

u th he thermochronometry constraints on unroofing of the eastern kaapvaal craton and significance for uplift of the southern african plateau the timing and causes of the 1 0 km elevation gain of the southern african plateau since paleozoic time are widely debated we report the first apatite and titanite u th he thermochronometry data for southern africa to resolve the unroofing history across a classic portion of the major escarpment that encircles the plateau the study area encompasses 1500 m of relief within archean basement of the barberton greenstone belt region of the eastern kaapvaal craton titanite dates are neoproterozoic apatite dates are cretaceous with most results clustering at ca 100 ma thermal history simulations confirm mesozoic heating followed by accelerated cooling in mid to late cretaceous time the lower temperature sensitivity of the apatite u th he method relative to previous thermochronometry in southern africa allows tighter constraints on the cenozoic thermal history than past work the data limit cenozoic temperatures east of the escarpment to 35 c and appear best explained by temperatures within a few degrees of the modern surface temperature these results restrict cenozoic unroofing to less than 850 m and permit negligible erosion since the cretaceous if substantial uplift of the southern african plateau occurred in the cenozoic as advocated by some workers then it was not responsible for the majority of post paleozoic unroofing across the eastern escarpment significant mesozoic unroofing is coincident with large igneous province activity kimberlite magmatism and continental rifting within and along the margins of southern africa compatible with a phase of plateau elevation gain due to mantle buoyancy sources associated with these events

the nature of the crust in southern india implications for precambrian crustal evolution 1 we present crustal thickness and poisson s ratio determinations from receiver function analyzes at 32 sites on the archaean and proterozoic terrains of south india the crustal thickness in the late archaean 2 5 ga eastern dharwar craton varies from 34 39 km similar crustal thickness is observed beneath the deccan volcanic province and the cuddapah basin the most unexpected result is the anomalous present day crustal thickness of 42 51 km beneath the mid archaean 3 4 3 0 ga segment of the western dharwar craton since the amphibolite grade metamorphic mineral assemblages 5 7 kbar paleopressures in this part of western dharwar craton equilibrated at the depths of 15 20 km our observations suggest the existence of an exceptionally thick 57 70 km crust 3 0 ga ago beneath the exhumed granulite terrain in southernmost india the crustal thickness varies between 42 60 km the poisson s ratio ranges between 0 24 0 28 beneath the precambrian terrains indicating the presence of intermediate rock type in the lower crust these observations of thickened crust suggest significant crustal shortening in south india during the archaean

mantle structure and dynamic topography in the mediterranean basin 1 we study the contribution of mantle flow to surface deformation within the mediterranean basin flow is modeled numerically based on lateral changes in mantle temperature estimated from tomography models we find that modeling results are significantly affected by the properties of the selected tomography models shear velocity models based on surface wave observations achieve the highest resolution of upper mantle structure and as a result are most successful in predicting microplate motion and dynamic topography

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

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

quest for dynamic topography observations from southeast asia transient dynamic topography is maintained by time varying density anomalies deep within earth9s mantle there is considerable interest in predicting and measuring dynamic topography because its spatial and temporal variation will help to track changes in the pattern of mantle convection despite the availability of increasingly elaborate predictive models there is little agreement about the amplitude wavelength and rate of change of dynamic topography here we analyze a hybrid data set of subsidence observations from continental and oceanic basins in southeast asia where published dynamic models predict 1 2 km of regional subsidence residual subsidence was isolated by removing the well defined effects of lithospheric extension and or cooling at long wavelengths 10 3 km the amplitude of residual subsidence i e maximum permissible dynamic subsidence has an upper bound of only 300 m the spatial distribution of permissible dynamic topography has no simple spatial or temporal relationship to the pattern of subduction over the past 65 m y anomalous subsidence commenced 5 10 m y ago accelerating to a present day rate of 100 m m y our results show that density anomalies within the lower mantle have little apparent influence at earth9s surface instead the distribution of residual subsidence reported here can be accounted for by minor variations in lithospheric thickness or by transient density variations beneath the lithospheric plate the absence of measurable dynamic topography is surprising and indicates that mantle mass anomalies are supported elsewhere presumably at internal boundaries within earth our hybrid data set can thus be used as an independent constraint on the viscosity structure of the mantle

dynamics of the india eurasia collision zone we present simple new dynamic calculations of a vertically averaged deviatoric stress field over a depth average of 100 km for asia from geodetic geologic topographic and seismic data a first estimate of the minimum absolute magnitudes and directions of vertically averaged deviatoric stress is obtained by solving force balance equations for deviatoric stresses associated with gravitational potential energy differences within the lithosphere plus a first order contribution of deviatoric stresses associated with stress boundary conditions this initial estimate of the vertically averaged deviatoric stress field is obtained independent of assumptions about the rheology of the lithosphere absolute magnitudes of vertically averaged deviatoric stresses vary between 5 and 40 mpa assuming bulk viscous behavior for the lithosphere the magnitudes of deviatoric stresses together with the magnitudes of strain rates inferred from quaternary fault slip rate and gps data yield vertically averaged effective viscosities for tibet of 0 5 5 1022 pa s compared with 1 2 5 1023 pa s in more rigid areas elsewhere in the region a forward modeling method that solves force balance equations using velocity boundary conditions allows us to refine our estimates of the vertically averaged effective viscosity distribution and deviatoric stress field the total vertically averaged deviatoric stress and effective viscosity field are consistent with a weak lower crust in tibet they are consistent with some eastward motion of tibet and south china lithosphere relative to eurasia and they confirm that gravitational potential energy differences have a profound effect on the spatially varying style and magnitude of strain rate around the tibetan plateau our results for the vertically averaged deviatoric stress argue for a large portion of the strength of the lithosphere to reside within the seismogenic upper crust to get deviatoric stress magnitudes there to be as high as 100 300 mpa in accord with laboratory and theoretical friction experiments indicating that stress drops in earthquakes are small fractions of the total deviatoric stress

gravity anomalies flexure and the elastic thickness structure of the india eurasia collisional system wehaveusedbouguergravityanomalyandtopographydatatodeterminetheequivalentelasticthicknessofthelithosphere te in the region of the india eurasia collisional system comparison of observed and modelled gravity anomalies along 1dimensional profiles suggest there are significant variations in te along strike of the himalaya foreland estimates decrease from a high of 70 km in the central region to 30 50 km in the east and west we have verified these inferences of spatial variations using a 2 dimensional non spectral interative flexure and gravity anomaly modelling technique the himalaya foreland forms a high te 40bteb100 km rigid block with a well defined edge as shown by the localisation of faulting and deformation along its northern margin other highteblocks occur to the north beneath the qaidam and sichuan basins the tibetan plateau forms a low te 0bteb20 km weak region that extends from the central part of the plateau into south western china tectonic styles in the india eurasia collisional system therefore involve both drigidt and dnon rigidt blocks where high te rigid blocks are present the stylesdominatedbyunderthrustingofthemorerigidblock wherethecollisionalzoneisnotconstrainedbyrigidblocks however the style appears to be dominated by lower crustal flow and a more continuous style of deformation d 2005 elsevier b v all rights reserved

the influence of dynamically supported topography on estimates of te abstract the effective elastic thickness t e of continental lithosphere can be estimated from the relationship between gravity and topography in the spectral domain two methods have been used one of which depends on the coherence between bouguer gravity anomalies and topography whereas the other uses the transfer function commonly known as the admittance between the free air gravity and topography the two methods give estimates of t e which differ by as much as an order of magnitude in those continental regions where variations in elevation are small this problem has led to much controversy an important concern is the extent to which estimates of t e are affected by dynamically maintained gravity and topography arising from mantle convection and postglacial recovery unlike elastically supported anomalies these processes can generate gravity and topographic long wavelength 500 km anomalies if such anomalies are modelled as being supported by elastic forces the resulting values of t e are overestimated often by a large amount

continental lithosphere folding in central asia part ii constraints from gravity and topography abstract periodical sub parallel fold like structures observed in the western gobi have two characteristic dominant wavelengths of approximately 50 and 300 360 km we explain these observations in terms of independent quasi viscous folding of the crustal and upper mantle parts of the lithosphere due to horizontal transpressional stresses in central asia induced by the collision with the indian plate we derive a model for the mechanical response of the continental lithosphere to horizontal stresses for a rheologically layered plate with non newtonian power law rheology overlying a low viscosity asthenosphere the differentiation in the effective viscosity and thickness of the strong upper crust and upper mantle lithosphere along with the presence of a low viscosity lower crust between them leads to their partial decoupling during compressional deformation these features result in the appearance of different wavelengths of the folds in the western gobi region we derive simple semi analytical estimates for the dominant wavelengths and rates of growth of surface undulations constrained by data from experimental rock mechanics topography and gravity

mantle dynamics in the mediterranean the mediterranean offers a unique opportunity to study the driving forces of tectonic deformation within a complex mobile belt lithospheric dynamics are affected by slab rollback and collision of two large slowly moving plates forcing fragments of continental and oceanic lithosphere to interact this paper reviews the rich and growing set of constraints from geological reconstructions geodetic data and crustal and upper mantle heterogeneity imaged by structural seismology we proceed to discuss a conceptual and quantitative framework for the causes of surface deformation exploring existing and newly developed tectonic and numerical geodynamic models we illustrate the role of mantle convection on surface geology a coherent picture emerges which can be outlined by two almost symmetric upper mantle convection cells the downwellings are found in the center of the mediterranean and are associated with the descent of the tyrrhenian and the hellenic slabs during plate convergence these slabs migrated backward with respect to the eurasian upper plate inducing a return flow of the asthenosphere from the backarc regions towards the subduction zones this flow can be found at large distance from the subduction zones and is at present expressed in two upwellings beneath anatolia and eastern iberia this convection system provides an explanation for the general pattern of seismic anisotropy in the mediterranean first order anatolia and adria microplate kinematics and may contribute to the high elevation of scarcely deformed areas such as anatolia and eastern iberia more generally the mediterranean is an illustration of how upper mantle small scale convection leads to intraplate deformation and complex plate boundary reconfiguration at the westernmost terminus of the tethyan collision

joint seismic geodynamic mineral physical modelling of african geodynamics a reconciliation of deep mantle convection with surface geophysical constraints recent progress in seismic tomography provides the first complete 3 d images of the combined thermal and chemical anomalies that characterise the unique deep mantle structure below the african continent we present a tomography based model of mantle convection that provides an excellent match to fundamental surface geodynamic constraints on internal density heterogeneity that includes both compositional and thermal contributions where the latter are constrained by mineral physics the application of this thermochemical convection model to the problem of african mantle dynamics yields a reconciliation of both surface gravity and topography anomalies to deep seated mantle flow under the african plate over a wider range of wavelengths than has been possible to date on the basis of these results we predict flow in the african asthenosphere characterised by a clear pattern of focussed upwellings below the major centres of late cenozoic volcanism including the kenya domes hoggar massif cameroon volcanicline cape verde and canary islands the flow predictions also reveal a deep seated large scale active hot upwelling below the western margin of africa under the cape verde islands that extends down to the core mantle boundary the scale and dynamical intensity of this west african superplume is comparable to the south african superplume that has long been assumed to dominate the large scale flow dynamics in the deep mantle under africa we evaluate the plausibility of the predicted asthenospheric flow patterns through a comparison with seismic azimuthal anisotropy derived from independent analyses of african shear wave splitting data

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

geothermal anomalies in the crust and upper mantle along southern rocky mountain transitions heat flow temperature gradient and gravity data are analyzed along five profiles crossing from the colorado plateau or great plains to the high topographic areas of the southern rocky mountains these profiles cross large negative bouguer gravity anomalies associated with the high topography in the region and are within the area studied by two recent seismic investigations the heatflow and gravity profiles typically cross from maximum to minimum values in the region allowing an fft fast fourier transform algorithm to be applied to the data so that high frequencies can be filtered the resulting smoothed data can be fitted by the first few terms of a sine series to estimate mid depths of crustal and upper mantle step sources although anomalies could also be caused by different source geometries at much shallower or much deeper depths the analysis provides a depth comparison for thermal and gravity anomaly sources modeled with similar geometries and a comparison of thermal source depths with seismic velocity anomalies there is considerable similarity along the profiles between the gravity and heat flow data and anomaly sources appear to be both in the crust and upper mantle there often appears to be spatial coincidence between thermal and gravity anomaly sources and seismic velocity anomalies proposed in the upper mantle implying thermal influence to these gravity and seismic anomalies thermal anomalies proposed to have sources in the crust in some areas are consistent with gravity and shear wave anomalies and support studies in those areas suggesting high crustal radiogenic concentrations in relatively less dense granitic crustal rocks

detecting the limit of slab break off in central turkey new high resolution pn tomography results summary inversion of pn traveltime residuals from a 39 station broad band array provides a highresolution image of the velocity structure in the uppermost mantle beneath central turkey individually picked pn phase arrivals from events recorded by the north anatolian fault passive seismic experiment and the kandilli observatory were combined with additional events associated with the eastern turkey seismic experiment tomography results show no change in pn velocity across the north anatolian fault although longitudinal variations are evident a region of very low pn velocities 8 1 km s 1 w est of the fault the sharp transition along the cafz which follows the palaeotectonic inner tauride suture may represent the location of the edge of the slab window created when the oceanic slab broke off along the bitlis zagros suture around 11 ma as the arabian plate collided with the eurasian plate

accurate measurements of residual topography from the oceanic realm in the oceans our understanding of plate subsidence as a function of age permits residual depth anomalies to be identified and mapped these anomalies may reflect dynamic topography and could be an important means for constraining convective circulation of the sublithospheric mantle here we analyze a global database of seismic reflection and wide angle profiles from heavily sedimented oceanic crust which abuts continental lithosphere at 449 locations we calculated water loaded subsidence compared it with a reference age depth relationship and determined residual depth we then combined these spot measurements of residual depth with observations from mid oceanic ridges and from selected ship track bathymetry to construct a global map of residual depth our results suggest that the amplitude of residual depth varies by up to 1 km with wavelengths of order 103 km we compare our residual depths with free air gravity and seismic tomographic anomalies our results show that residual depths correlate with long wavelength gravity anomalies in contrast correlations between residual depths and vertically averaged shear velocity anomalies within the upper and or the lower mantle are weaker the largest discrepancies occur at short 1000 km wavelengths these combined observations suggest that residual depth anomalies could be generate by density variations within a thin 102 km low viscosity layer beneath the lithosphere our global compilation should play a significant role in helping to refine predictive geodynamical models

dynamic topography above retreating subduction zones dynamic topography provides a measure of stresses within the earth s interior dense slabs induce an upper mantle flow that deflects the surface of the earth downward above them by combining a simple theoretical stokeslet model of subduction gravity modeling and seismic tomography i show that a significant fraction as much as 2000 m of the topographic variations observed above the scotia mariana and hellenic subduction systems appears to be dynamically induced by stresses related to the underlying subduction

measuring dynamic topography an analysis of southeast asia 1 theoretical models of dynamic topography generated by subducting slabs and consistent with the observed geoid predict 1 2 km of subsidence on wavelengths of 102 103 km the demonstrable existence of such subsidence has important implications for understanding sedimentary basin formation relative sea level changes as well as the earth s viscosity and density structure here we present the results of a combined analysis of subsidence data from oceanic marginal basins and their fringing continental shelves we have chosen to analyze a large area of southeast asia where dynamical predictions are believed to be well constrained by the history of subduction our joint analysis constrains the maximum amplitude of dynamic subsidence to be 300 m with a range of 0 500 m significantly less than predicted actively spreading marginal basins have no resolvable dynamic subsidence but after removing the effects of crustal thickness variation a few inactive marginal basins are consistently 200 500 m deeper than normal adjacent continental sedimentary basins record a possible dynamic subsidence of 0 500 m although wavelengths are shorter than expected 102 km there is also disagreement between the predicted and observed distribution of anomalous subsidence as well as inconsistencies in initiation and duration these results suggest that mass excess within the mantle represented by subducted slabs has little discernible effect on surface topography instead subsidence and marine inundations throughout southeast asia are predominantly controlled by lithospheric extension and or by foreland loading the spatial and temporal distribution of anomalous subsidence suggests that it might not be caused by dynamic topography associated with subducting slabs a large discrepancy between predicted and observed dynamic topography could be reconciled by changing the density structure of subducted slabs by supporting mass excess with a marked increase in viscosity at depth or by decoupling slabs from the surface with a very low viscosity layer

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

dynamic subsidence and uplift of the colorado plateau we use inverse models of mantle convection to explore the vertical evolution of the colorado plateau by satisfying multiple constraints seismic tomography stratigraphy in the western united states and great plains and other structural and volcanic data adjacent to the plateau the model provides predictions on the continuous history of colorado plateau vertical motion since 100 ma with the arrival of the flat lying farallon slab dynamic subsidence swept from west to east over the plateau and reached a maximum ca 86 ma two stages of uplift followed the removal of the farallon slab below the plateau one in the latest cretaceous and the other in the eocene with a cumulative uplift of 1 2 km both the descent of the slab and buoyant upwellings raised the plateau to its current elevation during the oligocene a locally thick plateau lithosphere enhances the coupling to the upper mantle so that the plateau has a higher topography with sharp edges the models predict that the plateau tilted downward to the northeast before the oligocene caused by northeast trending subduction of the farallon slab and that this northeast tilting diminished and reversed to the southwest during the miocene in response to buoyant upwellings

dynamic surface topography a new interpretation based upon mantle flow models derived from seismic tomography the very long wavelength i e in the degree range l 1 8 dynamic topography at the earth s surface estimated by isostatic reduction of the observed topography is shown to be in good agreement with the topography predicted by a model of viscous flow in the mantle this flow model was previously derived by independent consideration of the observed nonhydrostatic geoid and is based upon a recent model of large scale seismic shear velocity heterogeneity in the mantle our estimate and prediction of the long wavelength dynamic surface topography reveals significant depressions 2 3 km of the continental shields relative to their hydrostatic reference positions

controls of mantle plumes and lithospheric folding on modes of intraplate continental tectonics differences and similarities summary mantle plume activity and lithospheric folding by far field stresses exerted from plate boundaries are two important end members as mechanisms for continental intraplate deformation the topographic expression of mantle plume impingement on continental lithosphere and lithospheric folding has some striking similarities observations from a number of areas in europe s intraplate lithosphere demonstrate that these mechanisms commonly interact in space and time we present the results of thermomechanical modelling addressing the role of factors such as the presence of a hot upper mantle the spatial dimensions of the plume and the time constants involved in the temporal succession of plume activity and lithospheric folding by stress accumulation in intraplate continental lithosphere the results demonstrate that both the processes plume lithosphere interactions and folding may interact resulting either in strong amplification attenuation or modification of their surface expression these inferences are compatible with a number of key observations on the nature and the temporal succession of topography evolution in the alpine foreland the pannonian basin the scandinavian continental margin and the iberian peninsula

mafic dykes intrusive into pre cambrian rocks of the sao luis cratonic fragment and gurupi belt parnaiba province north northeastern brazil geochemistry sr nd pb o isotopes 40ar 39ar geochronology and relationships to camp magmatism abstract dykes of diabase and microgabbro intruded into pre cambrian rocks of the sao luis cratonic fragment and gurupi belt which are tectonic and erosive windows of the parnaiba basin in north northeastern brazil ar ar ages were determined and major trace element and nd sr pb o isotopic compositions of these dykes were measured to provide insights into their age and into the nature of their mantle sources and petrogenetic processes the data have also been used to compare the chemical and isotopic signatures of the dykes with those of the central atlantic magmatic province camp four chemical groups of mafic dykes have been identified these comprise two subtypes of high ti rocks i hti 1 tio 2 2 47 wt ii hti 2 tio 2 2 7 wt sio 2 47 wt in addition to iii evolved high ti tio 2 4 wt sio 2 of 46 wt and iv low si tio 2 2 2 wt sio 2 40 ar 39 ar geochronology of plagioclase returned ages of 201 4 ma and 193 10 ma for the hti 2 subtype and of 201 2 ma and 207 9 ma for the evolved high ti group the hti 1 and low si groups presented highly disturbed age spectra and did not allow the definition of their emplacement ages the argon data indicate an age 200 ma for the low si group and are dubious with respect to the age of thehti 1 subtype if coeval with i e 200 ma or older than the hti 2 and evolved high ti types all groups present 18 o values of pyroxene that are compatible with uncontaminated mantle derived magmas the hti 1 subtype average 143 nd 144 nd 200 0 512644 87 sr 86 sr 200 0 7035 206 pb 204 pb of 17 86 shows the less enriched and less fractionated more primitive trace element distribution of all groups the hti 2 subtype shows enriched trace element pattern and depleted nd sr signature 143 nd 144 nd 200 0 512610 87 sr 86 sr 200 0 7037 and average 206 pb 204 pb ratios of 17 23 the evolved high ti chemical group shows average ratios of 143 nd 144 nd 200 0 512558 87 sr 86 sr 200 0 7035 and 206 pb 204 pb of 16 88 and the more enriched trace elements signature among the four groups the chemical and isotopic compositions and trends of the hti 1 hti 2 ehti types are consistent with their derivation from an asthenosphere derived parental magma further modified by differentiation and minimal crust contamination higher in the hti 2 and ehti types and by the derivation of one type from another via fractional crystallization these high and evolved high tio 2 types show ages and some chemical and isotopic features that are consistent with those of the camp magmatism some differences found are ascribed to petrogenetic processes such as magma differentiation a combination of warming of the mantle and edge driven convection beneath the pangea supercontinent after the closure of the neoproterozoic brasiliano pan african orogenies in the ediacaran cambrian boundary might have triggered the magmatic event the low si type shows paired ta nb and zr hf depletions and depleted sr nd average 143 nd 144 nd 200 0 512687 87 sr 86 sr 200 0 703 and enriched pb 206 pb 204 pb of 18 66 isotopic compositions that may be interpreted to result either from interaction of a subcontinental lithospheric mantle with products of an earlier subduction or by contamination of the mantle derived magma during ascent and emplacement in the continental crust it is hypothesized that these dykes were emplaced in the ediacaran cambrian boundary after the neoproterozoic orogeny that built up the gurupi belt and in the early extensional stages that preceded the formation of the parnaiba basin

effects of a low viscosity lower crust on topography and gravity at convergent mountain belts during gravitational instability of mantle lithosphere we examine the gravitational stability of a stratified structure consisting of a low density crust through which viscosity decreases exponentially over a denser mantle lithosphere of constant viscosity which in turn overlies an inviscid slightly less dense layer like the asthenosphere the most important aspect of the viscosity structure is the contrast in viscosity at the moho surface uplift above a mantle downwelling is greatest when the viscosity contrast at the moho is negligible the surface is depressed when the viscosity of the lowermost crust is much greater than mantle viscosity for a wide range of viscosity structures calculated gravity anomalies and in particular admittances ratios of the fourier transform of gravity to that of topography produced by rayleigh taylor instability differ markedly from those observed across convergent mountain belts in particular whereas values of calculated admittance are negative or very small observed admittances are consistently positive thus if downwelling of mantle lithosphere directly beneath mountain ranges plays an important role in the construction of mountain belts it cannot be the main process that induces crustal thickening

laramide crustal thickening event in the rocky mountain foreland and great plains the rocky mountain foreland and great plains of the western united states were formerly part of a continental platform adjusted by erosion and deposition in cambrian through jurassic time to near mean sea level and therefore to a near uniform crustal thickness of approximately 33 km today the region stands at regional elevations up to 2 km isostatically supported by a crust exceeding 50 km in thickness reasonable estimates of tertiary sedimentation and laramide strain do not account for more than 15 of the implied thickening however from approximately 70 40 m y b p this region was underlain by a horizontally subducting slab of farallon plate lithosphere moving northeast this slab could have been the cause of the thickening finite difference thermal models with specified kinematics show only a temporary cooling of the base of the north american lithosphere by this slab the excess weight of the slab would have depressed the region plate bending calculations show a quantitative agreement of predicted depression with upper cretaceous isopachs since depression by this slab lasted until the eocene at least the latest cretaceous regression was probably caused by a laramide crustal thickening event the farallon slab might have caused crustal thickening in two ways its excess weight would have drawn in ductile lower crust from surrounding regions however calculations show that this effect is too slow too local and too reversible to explain most of the crustal thickening therefore it seems likely that ductile lower crust was transported from sw to ne by shear stresses which the farallon plate exerted on the base of the north american lithosphere a preliminary finite element calculation based on this hypothesis shows the correct general pattern of crustal thickening an unexpected but encouraging result is that predicted principal compression directions are orthogonal to many laramide basement uplifts

eucrust 07 a new reference model for the european crust 1 we present a new digital model eucrust 07 for the crust of western and central europe and surroundings 35n 71n 25w 35e available results of seismic reflection refraction and receiver functions studies are assembled in an integrated model at a uniform grid 15 0 15 0 the model consists of three layers sediments and two layers of the crystalline crust besides depth to the boundaries we provide average p wave velocities in the upper and lower parts of the crystalline crust the new model demonstrates large differences in the moho depth compared to previous compilations over 10 km in some specific areas e g the baltic shield furthermore the velocity structure of the crust is much more heterogeneous than in previous maps eucrust 07 offers a starting point for numerical modeling of deeper structures by allowing correction for crustal effects beforehand and to resolve trade off with mantle heterogeneities citation tesauro m m k kaban and s a p l cloetingh 2008 eucrust 07 a new reference model for the european crust geophys res lett 35 l05313 doi 10 1029 2007gl032244

lith5 0 a revised crustal model for canada based on lithoprobe results summary during the last 15 yr numerous seismic refraction experiments have been conducted overthe canadian landmass most as part of the lithoprobe initiative the new results providethe basis for revising existing global models of crustal structure e g crust5 1 mooney et al 1998 which do not incorporate the new lithoprobe data the crustal model discussedhere lith5 0 combines the lithoprobe and crust5 1 data to provide better constraints oncrustal structure in canada the objectives of the study are to 1 compile available new datainto a single consistent database of crustal thickness and velocity including the uppermostmantle pn velocity and 2 to construct a revised 5 5 crustal model of the canadianlandmass in a format that is generally compatible with existing global models the lithoproberesults suggest a 5 5 km average root mean square rms crustal thickness difference betweenlith5 0andcrust5 1 evenwiththeincorporationofthenewdata thegeographiccoveragein north and central canada remains sparse the new model does not explicitly incorporategeological constraints but instead relies on spherical spline interpolation for extrapolatingvalues into regions that lack control to illustrate the implications of this new crustal model we use lith5 0 to correct for the isostatically compensated crustal heterogeneity in order toobtain improved estimates of the dynamic surface topography the latter provides importantconstraints on the density structure of the uppermost mantle

extension during continental convergence with application to the tibetan plateau the tibetan plateau is the product of crustal thickening caused by the collision between india and asia and is the largest active example of extensional tectonics in a zone of continental collision throughout most of the tertiary the tectonics of the plateau were dominated by north south shortening a significant proportion of which took place on east west striking thrust faults for the last 5 m y or so the plateau has been thinning by the mechanism of extension on north south trending normal faults previous investigations of the collision have been able to account for the large scale features of the tertiary deformation but have failed to explain the transition in the late tertiary to quaternary strain rate field of the plateau from north south compression to east west extension the transition could in principle be effected by a reduction in the rate of convergence between india and asia or by a uniform reduction in strength of the whole continental lithosphere of asia it could not however be effected by a reduction in strength of the elevated region alone this produces increased compressional strain rates in the weakened zone an alternative explanation for the transition to extension comes from considering the thermal evolution of thickened continental lithosphere the lower part of the lithosphere consists of a thermal boundary layer which when thickened by horizontal shortening is colder and denser than its surroundings convective instability of the thickened thermal boundary layer and its replacement by hot asthenosphere would rapidly raise the surface elevation and gravitational potential energy of the overlying part of the lithosphere the convective instability would happen in a time brief compared with the collision time scale 50 m y in the case of india and asia but would only occur after there had already been substantial thickening of the lithosphere numerical experiments show that for a range of lithospheric parameters the increase in surface height as much as 2 km and of potential energy 5 to 10 1012 n m 1 resulting from convective instability of the lower lithosphere are sufficient for east west extension to replace north south compression as the dominant feature of the stress field of the tibetan plateau

flexure of the india plate underneath the bhutan himalaya 1 we investigate flexural geometry and rheology of the india plate beneath the eastern himalaya from a new gravity data set acquired in bhutan compared to the well studied nepal himalaya the obtained bouguer anomaly profiles across the range show shorter wavelength flexure of the lithosphere with a narrower and shallower foreland basin this new data set and seismic moho constraints are interpreted together in terms of lithospheric flexure using a 2 d thermomechanical model our results suggest that the strongest layer of the continental lithosphere beneath bhutan is the upper mantle as it is beneath nepal the observed west to east decrease in flexural wavelength is associated with weakening mantle rheology the simulations show that this decrease can be related to ductile mantle behavior either hydrated wet dunite or dry and hot olivine rheology both models display decoupled lithospheric layers leading to an eastward decrease of flexural rigidity from 10 24 to 5 10 22 n m in nepal and bhutan respectively

large scale crustal heterogeneities and lithospheric strength in cratons abstract the rheology and thermal structure of the continental lithosphere are intimately linked in old cratons the effective elastic thickness of the lithosphere has been estimated by various spectral inverse methods based on the correlation between topography and gravity anomalies estimates vary within a very large range from 40 km to 120 km depending on the method used in this paper we use forward models to account for lateral variations in mechanical properties and their effect on the equivalent elastic thickness eet of the lithosphere from these models which allow brittle elastic ductile rheologies and mechanical discontinuities faults we have calculated the strain stress distributions and displacement fields vertical integration of the stress permits a local determination of the effective elastic thickness the computed displacements were used to calculate related bouguer and free air gravity anomalies and compare them with the observations the analysis is applied to the 2000 ma kapuskasing uplift in the superior province of the canadian shield where the presence of a high density block in the upper crust is due to the upthrusting of midcrustal rocks along a major thrust fault the study shows that the stability of this structure on geologic time scales requires a strong lower crustal rheology a cold geotherm and the fault to be healed this study also shows that because of stress dependence of the non linear rheology crustal heterogeneities may cause significant 40 local reductions of the lithospheric strength away from the kapuskasing structure the average strength of the lithosphere remains high eet 100 km conventional methods for estimating the elastic thickness would not resolve such local strength reductions in cratons but would predict depending on the method used highly overestimated or instead underestimated eet

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

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

buildup of a dynamically supported orogenic plateau numerical modeling of the zagros central iran case study the iranian plateau is a vast inland region with a smooth average elevation of c 1 5 km formed at the rear of the zagros orogen as a result of the arabia eurasia collision i e over the last 30 35 myr this collision zone is of particular interest due to its disputed resemblance to the faster himalayan collision which gave birth to the tibetan plateau around 50 myr ago recent studies have suggested that a recent 10 5 ma slab break off event below central iran caused the formation of the iranian plateau here we test several hypotheses through large scale 3082 3 590 km numerical models of continental subduction mod els that incorporate a free upper surface erosion rheological stratification brittle elastic ductile rheologies and metamorphic phase changes density and physical properties and account for the specific crustal and thermal structure of the arabian and iranian continental lithospheres we test the impact of the transition from oceanic to continental subduction and the topographic consequences of the progressive slowdown of the convergence rate during continental subduction our results demonstrate the role of mantle flow beneath the overriding plate initiated as an indirect consequence of slab break off this flow creates a dynamic topography support during continental subduction and results in delamination of the overriding plate lithospheric mantle followed by isostatic readjustment hence of further uplift and maintenance of a plateau like topography without significant crustal thickening the slowdown of the convergence rate dur ing the development of the continental subduction collision phase largely contributes to this process by controlling the timing and depth of slab break off

the geoid constraint in global geodynamics viscosity structure mantle heterogeneity models and boundary conditions summary the earth s non hydrostatic gravity field or geoid provides a first order constraint on mantle density structure and dynamics geodynamic models for the geoid have proliferated since the advent of seismic mapping of mantle heterogeneity structure tomography because the geoid offers perhaps the best measured independent constraint on mantle density heterogeneity however dynamic geoid models involve a number of questionable physical assumptions and uncertainties whose effects need to be evaluated before geodynamic inferences based upon the geoid can be considered sound troubling issues include the appropriate surface boundary conditions free slip no slip plates and parametrization of radial viscosity variations how many layers can be resolved in addition to lateral viscosity variations possible chemical layering of the mantle phase transitions etc there are also uncertainties in the density heterogeneity models used as input to dynamic geoid models most of which are derived from seismic tomography and require weakly constrained empirical conversion factors to go from seismic velocity variations to density variations here we address several of the most straightforward problems inherent in geoid modelling namely the issues of viscosity structure resolution uncertainties in appropriate boundary conditions and differences among mantle heterogeneity models a robust feature of all models is a lower mantle viscosity at least a factor of 30 greater than that of the upper mantle but there is little resolution with regard to finer details such as lithospheric or uppermost mantle low viscosity zone viscosity ironically free slip boundary conditions result in the best fits to the geoid in all cases but all boundary conditions exhibit predictable trade offs with the uppermost mantle viscosity models with a single viscosity layer representing the lower mantle yield similar dynamic topography estimates of the order of 700 1000 m in amplitude regardless of the finer details of upper mantle viscosity structure boundary conditions or input heterogeneity models comparing mantle heterogeneity models based on two independent seismological determinations harvard and berkeley models and on the history of subduction we find that these models are virtually indistinguishable regarding inferences of mantle viscosity structure and amplitude of dynamic topography and in terms of the effects of different boundary conditions uncertainties concerning which type of boundary condition is appropriate are much more important than which mantle heterogeneity model is chosen given other uncertainties in modelling the geoid particularly the strong effects due to lateral viscosity variations for intermediate 10 000 km wavelengths we conclude that the class of dynamic geoid models explored so far cannot reliably elucidate the details of upper mantle viscosity structure

parallel thrust and normal faulting in peru and constraints on the state of stress abstract we present evidence for recent normal faulting in the cordillera blanca of peru along planes trending approximately n40 w parallel to the chain this type of faulting suggests that there is an important component of extension perpendicular to the chain and contrasts markedly with the underthrusting of the nazca plate on the west side of peru and with the folding and thrusting observed in the sub andes on the eastern flank of the chain these contrasting styles of deformation and associated stress fields can be explained by buoyancy forces due to the high mountains and crustal roots in so far as this explanation is correct the difference in elevation between zones of normal and thrust faulting can be used to place an upper bound of 500 bars on the average stress difference in the crust

lithospheric velocity structure of the anatolian plateau caucasus caspian region 1 the anatolian plateau caucasus caspian region is an area of complex lithospheric structure accompanied by large variations in seismic wave velocities despite the complexity of the region little is known about the detailed lithospheric structure using data from 31 new permanent broadband seismic stations along with results from a previous 29 temporary seismic stations and 3 existing global seismic stations in the region a 3 d velocity model is developed using joint inversion of teleseismic receiver functions and surface waves both group and phase dispersion curves love and rayleigh were derived from regional and teleseismic events additional rayleigh wave group dispersion curves were determined using ambient noise correlation receiver functions were calculated using p arrivals from 789 teleseismic 30 90 earthquakes the stacked receiver functions and surface wave dispersion curves were jointly inverted to yield the absolute shear wave velocity to a depth of 100 km at each station the depths of major discontinuities sediment basement crust mantle and lithosphere asthenosphere were inferred from the velocity depth profiles at the location of each station distinct spatial variations in crustal and upper mantle shear velocities were observed the kura basin showed slow 2 7 2 9 km s upper crustal 0 11 km velocities but elevated 3 8 3 9 km s velocities in the lower crust the anatolian plateau varied from 3 1 3 2 in the upper crust to 3 5 3 7 in the lower crust while velocities in the arabian plate south of the bitlis suture were slightly faster upper crust between 3 3 and 3 4 km s and lower crust between 3 8 and 3 9 km s the depth of the moho which was estimated from the shear velocity profiles was 35 km in the arabian plate and increased northward to 54 km at the southern edge of the greater caucasus moho depths in the kura and at the edge of the caspian showed more spatial variability but ranged between 35 and 45 km upper mantle velocities were slow under the anatolian plateau but increased to the south under the arabian plate and to the east 4 3 4 4 km s under the kura basin and greater caucasus the areas of slow mantle coincided with the locations of holocene volcanoes differences between rayleigh and love dispersions at long wavelengths reveal a pronounced variation in anisotropy between the anatolian plateau and the kura basin

topography and geoid due to lithospheric mass anomalies summary a model of lithospheric thickness and a recent compilation of moho depths are used to compute the earth s isostatic surface topography and associated gravity anomalies the results are strongly influenced by the uncertainties in lithospheric depth and crustal density profiles the preferred models explain most of the observed topography and are highly correlated with observed gravity anomalies for harmonic degrees larger than 10 comparisons of our residual topography with geodynamical calculations of dynamic topography based on mantle circulation are rather poor

images of crustal variations in the intermountain west 1 we develop a map of crustal thickness variations across the great basin colorado plateau rocky mountain and great plains provinces of the western united states using common conversion point stacking of teleseismic receiver functions below the rocky mountains and high plains in colorado we find the thickest crust in the region at 45 50 km thick beneath the basin and range thinner between 30 and 40 km crust is found thin 30 km thick crust is present in the northern portion of nevada and utah despite elevations similar to those farther south crustal thickness across the colorado plateau can be characterized as a broad transitional region between the thin crust of basin and range to the thicker crust of the rocky mountains the impedance contrast across the mohorovicic discontinuity decreases below the colorado plateau as converted arrivals recorded in this region appear weak compared to surrounding areas variations in vp vs across the region indicate higher values along the western boundary of the basin and range in the rocky mountains and in the western great plains we are not able to characterize vp vs in the colorado plateau we find that crustal thickness does not closely correlate with surface topography within each region or across the region as a whole differences in crustal thickness in each tectonic province indicate the need for a mantle component to support the high elevations across the western united states

lithosphere of the dharwar craton by joint inversion of p and s receiver functions summary the archean dharwar craton in south india is known for long time to be different from most other cratons specifically at station hyderabad hyb the ps converted phases from the 410and 660 km mantle discontinuities arrive up to 2 s later than in other cratons of comparable age which implies lower upper mantle velocities to resolve the unique lithosphere asthenosphere system of the dharwar craton we inverted jointly p and s receiver functions and teleseismic p and s traveltime residuals at 10 seismograph stations this method operates in the same depth range as long period surface waves but differs by much higher lateral and radial resolution we observe striking differences in crustal structures between the eastern and western dharwar craton edc and wdc respectively crustal thickness is of around 31 km with predominantly felsic velocities in the edc and of around 55 km with predominantly mafic velocities in the wdc in the mantle we observe significant variations in the p velocity with depth practically without accompanying variations in the s velocity in the mantle s velocity there are azimuthdependent indications of the hales discontinuity at a depth of 100 km the most conspicuous feature of our models is the lack of the high velocity mantle keel with the s velocity of 4 7 km s 1 typical of other archean cratons the s velocity in our models is close to 4 5 km s 1 from the moho to a depth of 250 km there are indications of a similar upper mantle structure in the northeast of the indian craton and of a partial recovery of the normal shield structure in the northwest a division between the high s velocity western tibet and low s velocity eastern tibet may be related to a similar division between the northeastern and northwestern indian craton

lithospheric structure of the aegean obtained from p and s receiver functions combined p and s receiver functions from seismograms of teleseismic events recorded at 65 temporary and permanent stations in the aegean region are used to map the geometry of the subducted african and the overriding aegean plates we image the moho of the subducting african plate at depths ranging from 40 km beneath southern crete and the western peloponnesus to 160 km beneath the volcanic arc and 220 km beneath northern greece however the dip of the moho of the subducting african plate is shallower beneath the peloponnesus than beneath crete and rhodes and flattens out beneath the northern aegean observed p to s conversions at stations located in the forearc indicate a reversed velocity contrast at the moho boundary of the aegean plate whereas this boundary is observed as a normal velocity contrast by the s to p conversions our modeling suggests that the presence of a large amount of serpentinite more than 30 in the forearc mantle wedge which generally occurs in the subduction zones may be the reason for the reverse sign of the p to s conversion coefficient moho depths for the aegean plate show that the southern part of the aegean crustal thickness of 20 22 km has been strongly influenced by extension while the northern aegean sea which at present undergoes the highest crustal deformation shows a relatively thicker crust 25 28 km this may imply a recent initiation of the present kinematics in the aegean western greece crustal thickness of 32 40 km is unaffected by the recent extension but underwent crustal thickening during the hellenides mountains building event the depths of the aegean moho beneath the margin of the peloponnesus and crete 25 28 and 25 33 km respectively show that these areas are also likely to be affected by the aegean extension even though the cyclades crustal thickness of 26 30 km were not significantly involved in this episode the aegean lithosphere asthenosphere boundary lab mapped with s receiver functions is about 150 km deep beneath mainland greece whereas the lab of the subducted african plate dips from 100 km beneath crete and the southern aegean sea to about 225 km under the volcanic arc this implies a thickness of 60 65 km for the subducted african lithosphere suggesting that the aegean lithosphere was not significantly affected by the extensional process associated with the exhumation of metamorphic core complexes in the cyclades

the relationship between topography and gravity on earth and venus abstract the measurement of the admittance between the topography and the gravity of a planet as a function of wavelength provides important constraints on its interior processes a detailed discussion of the influence of noise on such measurements shows that at long wavelengths the admittance is best determined in the frequency domain using gravity as input and topography as output this procedure is used to estimate the admittance in the pacific and indian oceans the observed values are reduced by volcanism above plumes which produces a contribution to the long wavelength compensated topography when this effect is taken into account there is little evidence that the admittance is reduced by the presence of a low viscosity zone in the upper mantle a similar study of venus shows an excellent correlation between gravity g and topography h for wavelengths greater than 700 km with an admittance of about 50 mgal km 1 this value agrees well with that expected from numerical axisymmetric models of plumes and can be used to construct maps of residual topography by calculating h g 50 where h is in kilometers and g in milligals the power spectrum of the gravity field and its skewness with positive anomalies covering a smaller area but being larger than the negative ones are also in agreement with those of the numerical models of the major topographic features only aphrodite and ishtar terrae behave differently with a value of admittance of about one third of other large structures this difference requires their topography to be partly supported by variations in crustal thickness by analogy with similar structures on earth both are likely to be supported by volcanically thickened crust on top of rising plumes though the observed relationships do not rule out other models the behavior of the admittance at long wavelengths is important because it can constrain the thickness of the convecting layer at wavelengths of more than 3000 km the gravity data provides some evidence that the layer depth is about 1000 km over the center of beta and phoebe regiones there are short wavelength anomalies of residual topography produced by crust that is thinner than that elsewhere flanked by positive anomalies of thicker crust the regions of thin crust occur beneath rift zones this behavior can be explained if crust flows down the flanks of elevated regions above plumes and can generate coronae when the plume support disappears

the role of gravitational potential energy in active deformation in the southwestern united states in orogenic belts intrinsic buoyancy and extrinsic boundary forces act on the continental lithosphere determining the relative importance of these forces is fundamental to understanding the dynamics of continental deformation in the western united states estimates of gravitational potential energy show that given available constraints on rheology buoyancy forces are sufficient to produce the wide variety of rates and styles of observed deformation including absence of deformation this underscores the importance of gravitational potential energy for deformation in orogenic belts

an analog experiment for the aegean to describe the contribution of gravitational potential energy the southern aegean seafloor exhibits clear evidence of internal deformation stretching as shown by tectonics seismology and space geodesy we use an analog three layer laboratory experiment of sand silicone putty and honey to investigate the deformation of the southern aegean lithosphere the model is installed in a box and confined by a vertical wall we open a gate in the wall and observe the deformation of the two upper layers due to buoyancy forces the general pattern of the deformation of the southern aegean is found in the analog model we observe the formation of an arc spreading outward with time the extension is radial in the inner part but parallel to the arc in the external part and of comparable importance at both ends of the gate we observe strike slip motion dextral in the western part sinistral in the eastern part rotation clockwise in the western part counterclockwise in the eastern part of up to 40 is seen on both sides of the gate but is also present with a smaller amplitude far in the internal region partially due to distributed shear the spreading is associated with the thinning of the two upper layers and affects a region of dimensions comparable to the length of the free boundary this spreading does not propagate inward with time some pieces of material located near the active boundary remain undeformed during the experiment

a new isostatic residual gravity map of the conterminous united states with a discussion on the significance of isostatic residual anomalies to display more clearly the gravity anomalies caused by geologic bodies in the upper parts of the crust a new colored isostatic residual gravity map of the conterminous united states has been prepared using the gravity data set compiled for the gravity anomaly map of the united states society of exploration geophysicists 1982 the new isostatic residual gravity map is based on an airy heiskanen model of local compensation in which the surface load requiring compensation is defined by 5 min topographic and bathymetric data sets a colored first vertical derivative map of isostatic residual gravity further enhances the short wavelength anomalies produced by bodies at or near the surface and emphasizes the regional fabrics and trends in the gravity field for the purpose of displaying gravity anomalies caused by shallow bodies of geologic significance the nature of the isostatic model and the values of its parameters are of lesser importance than the application of an isostatic correction of some sort most isostatic models result in residual gravity maps that appear nearly identical in their main patterns and features anomalies on isostatic residual gravity maps should not be casually interpreted in terms of undercompensation or overcompensation because large amplitude anomalies can be produced by crustal bodies in complete local isostatic equilibrium many isostatic residual gravity anomalies less than several hundred kilometers wide can be related to known geologic bodies we present here a classification scheme that attempts to categorize such anomalies on the basis of tectonic environment in general highs correlate with intruded or accreted mafic material or with upthrusted crustal sections whereas lows occur over low density sedimentary or volcanic sections felsic intrusive bodies or down warped crustal sections although some longer wavelength anomalies such as the broad gravity high centered over montana could be manifestations of density contrasts deep in the mantle many such anomalies can also be modeled by geologically reasonable density contrasts that are isostatically compensated and confined to depths of less than several hundred kilometers so that their source bodies need not be deep the fact that certain of these broader anomalies have well defined boundaries which correlate with near surface geologic features increases the likelihood that their sources lie entirely within the lithosphere if so then the density contrasts required to explain the gravity data imply fundamental anomalies in chemical composition or thermal state for the crustal and upper mantle columns under these regions we have investigated spectral analysis as a method to quantitatively characterize regional anomaly patterns contoured plots of normalized amplitude spectra were prepared for various areas of the isostatic residual gravity field of the united states these fourier domain representations show characteristic patterns that can be interpreted in terms of the trends and wavelengths of anomalies and may help to more objectively distinguish geologic basements with different origins or tectonic histories

the interaction of plume heads with compositional discontinuities in the earth s mantle the effects of compositional discontinuities of density and viscosity in the earth s mantle on the ascent of mantle plume heads is studied using a boundary integral numerical technique three specific problems are considered 1 a plume head rising away from a deformable interface 2 a plume head passing through an interface and 3 a plume head approaching the surface of the earth for the case of a plume attached to a free surface the calculated time dependent plume shapes are compared with experimental results two principal modes of plume head deformation are observed plume head elongation or the formation of a cavity inside the plume head the inferred structure of mantle plumes namely a large plume head with a long tail is characteristic of plumes attached to their source region and also of buoyant material moving away from an interface and of buoyant material moving through an interface from a high to low viscosity region as a rising plume head approaches the upper mantle most of the lower mantle will quickly drain from the gap between the plume head and the upper mantle if the plume head enters the upper mantle if the plume head moves from a high to low viscosity region the plume head becomes significantly elongated and for the viscosity contrasts thought to exist in the earth could extend from the 670 km discontinuity to the surface plume heads that are extended owing to a viscosity decrease in the upper mantle have a cylindrical geometry the dynamic surface topography induced by plume heads is bell shaped when the top of the plume head is at depths greater than about 0 1 plume head radii as the plume head approaches the surface and spreads the dynamic topography becomes plateau shaped the largest stresses are produced in the early stages of plume spreading when the plume head is still nearly spherical and the surface expression of these stresses is likely to be dominated by radial extension as the plume spreads compressional stresses on the surface are produced beyond the edges of the plume consequently extensional features will be produced above the plume head and may be surrounded by a ring of compressional features

estimating the compensation depth of the hawaiian swell with linear filters broad areas of anomalously shallow seafloor surround regions of active or recently active hot spot volcanism the cause of these midplate swells has been modeled both as elevated temperatures in a convecting layer and as thermal expansion within the conducting portion of the lithosphere according to the former explanation we would expect the low density material compensating the swell to lie beneath the conducting lid whereas the latter model allows compensation at shallower depths we present a technique based on the application of linear filters for determining the compensation depth of midplate swells given observations of topography and geoid or gravity anomalies our method does not require an a priori estimate of swell height or width nor does it assume complete spatial or spectral separation between the swell topography and the signal from the hot spot volcanoes which are compensated by flexure of the oceanic moho rather we exploit the difference in the predicted sign and amplitude of the moho deflection beneath elevated topography produced by loading above versus loading below the elastic plate in order to separate the two effects an application of these linear filters to residual topography and geoid data from the region surrounding the hawaiian hot spot yields a compensation depth of 70 10 km within the lower lithosphere

mantle flow uplift of western anatolia and the aegean interpretations from geophysical analyses and geodynamic modeling 1 the western anatolian and aegean region demonstrates a complex geologic history of horizontal and vertical tectonics active normal faulting and exhumation zones indicate that western anatolia has experienced significant extension since the oligocene early miocene 30 ma our geophysical analyses demonstrate that the region is also uplifted relative to an elevation that would be expected given an isostatic response to the lithosphere structure namely topography residuals indicate a residual uplift of about 1500 m over 200 km sections of western anatolia and the aegean admittance functions between free air gravity and topography indicate that the regional topography is isostatically uncompensated and as it approaches 50 mgal km at the longest wavelengths the uncompensated topography is likely owing to an underlying mantle flow component using forward geodynamic modelling we consider an idealized section of western anatolian lithosphere based on tomographic inversions and examine the magnitude and pattern of surface topography to reconcile with the geophysical observables the models consistently show a plateau type uplift and horizontal extension through western anatolia with an amplitude and wavelength consistent with the residual topography calculations together the geophysical analyses and modelling provide independent quantitative evidence that the thin anatolian aegean lithosphere is being buoyed upwards by underlying mantle flow the mantle flow may be associated with active lithosphere delamination beneath the region a process that would also explain the ongoing crustal extension

three dimensional modeling of plume lithosphere interaction in order to understand better the dynamics of hot spots such as hawaii we present a three dimensional numerical model for the interaction of a thermal plume with a moving lithosphere the model domain is a rectangular box filled with fluid whose viscosity depends upon temperature and pressure the lithosphere is represented by a layer of cold highly viscous fluid moving with an imposed horizontal velocity u in the x direction and a thermal plume is generated by a circular temperature anomaly on the bottom of the box the steady flow is determined numerically using a hybrid spectral finite difference technique the flow is characterized by a stagnation streamline of width y x1 5 that represents the edge of the spreading plume material we illustrate the detailed behavior of the model using the example of the hawaiian plume our best fitting hawaiian model is obtained by adjusting the plume buoyancy flux b until the predicted topography anomaly matches the observed hawaiian swell topography we find b 4100 kg s 1 which implies a plume radius of 90 km for an assumed plume mantle temperature contrast of 300 c the predicted topography is supported primarily by density anomalies beneath the lithosphere and cannot be explained by lithospheric erosion we therefore conclude that the classical lithospheric reheating model is unable to account for hotspot swells the horizontal flux of buoyancy associated with the swell exceeds b by up to 80 suggesting that current estimates of b for mantle plumes are too high empirically derived scaling laws for the width of the stagnation streamline and for the topography anomaly exhibit power law dependencies on b and u that agree well with those predicted by the refracted plume model of olson 1990 the principal weakness of the model is that the predicted geoid topography ratio of 0 010 for the hawaiian swell is about twice the observed value of 0 004 0 006