experimental calibration of the aluminum in hornblende geobarometer with application to long valley caldera california volcanic rocks a new geobarometer based on the al content of igneous hornblendes in equilibrium with melt fluid biotite quartz sanidine plagioclase sphene and magnetite or ilmenite has been calibrated experimentally the calibration was performed by equilibrating the required phase assemblage over the pressure range 2 8 kbar at 740 780c and then analyzing euhedral hornblendes in equilibrium with glass melt experiments were performed on natural samples of both volcanic and plutonic rocks earlier empirical calibrations of this geobarometer relied on analyzing natural hornblendes from plutons with the required phase assemblage and inferring pressure from nearby metamorphic country rocks the experimental calibration differs from the empirical calibrations especially above 5 kbar and shows that the al content of hornblendes in equilibrium with the required phase assemblage is greater for a given total pressure than previously thought the geobarometers uncertainty is dramatically reduced the derived equation is p plus minus 0 5 kbar 3 minus 3 46 plus minus 0 24 4 23 plus minus 0 13 al sup t the geobarometer is applied to post bishop tuff volcanic rocks from long valley caldera california and reveals that most rhyodacites in this complex erupted from depths of about 6 km these eruptions occurred over 500 000 yr suggesting that the rhyodacitic magma reservoir beneath longmore valley had reached a steady p depth t state less

anelastic response of the earth to a dip slip earthquake the deformation of the earth s surface immediately following an earthquake is well modeled by a dislocation in an elastic half space some time after the event however the elastic stresses begin to relax due to flow in the more fluid regions of the upper mantle the asthenosphere this relaxation causes further slow displacements observable at the earth s surface this paper examines the pattern and timing of these postseismic motions for a vertical dip slip fault and a 30 dip thrust fault a variety of rheological models are studied these models encompass variations in asthenosphere thickness mesosphere viscosity fault depth and most important power law flow both newtonian n 1 and non newtonian n 3 maxwell rheologies are considered for the asthenosphere and mesosphere the results for a vertical dip slip fault show that a peripheral upwarp develops within a distance of a few lithosphere thicknesses from the rupture zone on the downthrown side of the fault a comparable downwarp develops on the upthrown side these peripheral warps grow to a maximum height of about 5 of the total fault slip some time after the earthquake event and then collapse at much greater times the vertical motions for a 30 dip thrust fault show that after the earthquake a general uplift occurs in the region near the rupture zone horizontal displacements propagate asymmetrically away from the rupture zone the largest displacements occurring in the overthrust plate the principal result of the comparison of n 1 and n 3 rheological laws is that the basic pattern of uplifts and horizontal motions is very similar the major difference is the time dependence of the development of the pattern in spite of the complex geometrical structure of the subduction zone the displacements ui 1 t for newtonian rheology are nearly the same as ui 3 t for the non newtonian rheology if t t3 thus although the qualitative pattern of the displacements is not diagnostic of power law flow the great difference in time dependence allows a clear discrimination of the two flow laws

neodymium isotope systematics of the mantle beneath the baltic shield evidence for depleted mantle evolution since the archaean abstract more than 300 published and unpublished nd isotopic analyses of mantle derived rocks from the baltic shield have been compiled the rocks range in age from archaean to phanerozoic within any given age interval the mantle derived rocks range in nd t from depleted mantle values at or above the growth curves of the global depleted mantle reservoirs of depaolo 1981 and depaolo et al 1991 to negative values initial neodymium isotopic compositions below the de paolo curve are best explained by interaction between depleted mantle derived magmas and local crustal contaminants the data now available lend no support to the existence of isolated less depleted or undepleted mantle domains beneath the baltic shield as was suggested by mearns et al 1986 and valbracht 1991a b

did postglacial catastrophic flooding trigger massive changes in the black sea gas hydrate reservoir after the last glacial maximum the semi land locked black sea basin was flooded by warm water from the mediterranean sea this major sea level rise and change of physical water properties had a large impact on the gas hydrate reservoir in the sediments below modelling of the regional response of the gas hydrate stability zone ghsz to the black sea flooding 7100 years ago shows that a strong effect of near bottom temperature increase pushes the gas hydrate reservoir to a large shrinking of 15 62 that may release up to 1 1 4 6 gt of methane this catastrophic scenario is however delayed because of the transient nature of the heat wave propagation the large scale reduction of the ghsz is only to take place within the next thousand years at present widespread hydrate dissociation is only expected to occur where there is a minimum water depth for hydrate stability

numerical comparison of different convergent plate contacts subduction channel and subduction fault summary at convergent plate boundaries the properties of the actual plate contact are important for the overall dynamics convergent plate boundaries both mechanically decouple and link tectonic plates and accommodate large amounts of strain we investigate two fundamental physical states of the subduction contact one based on a fault and the other based on a subduction channel using a finite element method we determine the specific signatures of both states of the subduction contact we pay particular attention to the overriding plate in a tectonic setting of converging plates where the subducting plate is freely moving the subduction channel reduces compression relative to the fault model in a land locked basin setting where the relative motion between the far field of the plates is zero the subduction channel model produces tensile stress regime in the overriding plate even though the amount of slab roll back is small the fault model shows a stronger development of slab roll back and a compressive stress regime in the upper plate based on a consistent comparison of fault and channel numerical models we find that the nature of the plate contact is one of the controlling factors in developing or not of backarc extension we conclude that the type of plate contact plays a decisive role in controlling the backarc state of stress to obtain backarc extension roll back is required as an underling geodynamic process but it is not always a sufficient condition

foredeeps versus subduction zones the geologic characteristics of foredeeps and accretionary wedges suggest that these features are distinguishable on the basis of the direction of the associated subduction east north east dipping subduction related accretionary wedges show high relief and broad outcrops of metamorphic rocks they are associated with shallow foredeeps with low subsidence rates in contrast west dipping subduction related accretionary wedges show low relief and involve mainly sedimentary cover the related foredeeps are deep and have high subsidence rates this differentiation is useful both for oceanic and continental subductions e g eastern vs western pacific subductions or east dipping alpine vs west dipping apenninec subductions in a cross section of the alps the ratio of the area of the orogen to the area of the foredeeps is at least 2 1 whereas this ratio is 0 22 1 for the apennines these ratios explain why foredeeps related to east or northeast dipping subduction are quickly filled and are bypassed by clastic rocks whereas foredeeps related to west dipping subduction maintain a deep water environment longer these differences support the presence of an eastward asthenospheric counterflow relative to the westward drift of the lithosphere detected in the hot spot reference frame even in the mediterranean where no hot spots are present in this interpretation the apennines foredeep was caused by the eastward push of the mantle acting on the subducted slab whereas the foredeeps in the alps were caused by the load of the thrust sheets and the downward component of movement of the upper adriatic plate these forces contrast with the upward component of the eastward mantle flow

numerical modeling of frontal and basal accretion at collisional margins 1 we investigate the deformation of orogenic wedges that form in the early stages of continent continent collisions using a two dimensional numerical model limited to the upper lithosphere our models show that deformation at the plate margins is influenced by rheology surface processes and the balance between inward mass flux and outward subduction flux as controlled by the subduction load which represents the effects of slab pull and resistive forces and flexural downbending we find three characteristic deformation modes 1 near pure subduction with little or no accretion 2 frontal accretion with development of an accretionary wedge built up by offscraping of the sediment layer at shallow depth and 3 independent frontal and basal accretion where a retrothrust allows stacking of basement nappes at crustal to mantle depths near pure subduction is enabled for ordinary rheology materials characterized by brittle and viscous material behavior approximating a christmas tree type depth profile and almost zero friction along the subduction shear zone frontal accretion occurs when slightly increased friction along the subduction shear zone allows offscraping of the sediment layer from the subducting plate independent frontal and basal accretion develops in strong rheology models with an almost fully brittle material behavior major surface erosion or a reduction of the subduction load promote the development of large basement nappes frontal accretion is favored by major sedimentation during convergence a large backstop and in the case of a lateral transition from a strong rheology to an ordinary rheology subducting plate our numerical models develop first order characteristics as observed in natural orogenic wedges for example upper crustal nappe stacks frontal and basal accretion or extension in the core of an orogen frontal and basal accretion are interdependent and tend to stabilize the subduction system

mineralogy and geochemistry of devonian ultramafic minor intrusions of the southern kola peninsula russia implications for the petrogenesis of kimberlites and melilitites minor magmatic intrusions of kimberlite melilitite and cpx melilitite occur in the southern part of the kola peninsula russia on the terskii coast and near the town of kandalaksha they yield k ar ages of 382 14 ma and 365 16 ma similar to the magmatic rocks from the kola alkaline province the terskii coast kimberlites have mineralogical and geochemical affinities with group 1 kimberlites whereas the kandalaksha monticellite kimberlite more closely resembles calcite kimberlites the lower al2o3 content in the kola kimberlites indicates a strongly depleted harzburgitic source while higher al2o3 in the melilitites suggests a lherzolitic source the terskii coast kimberlites are anomalously potassic and significantly enriched in p and ba compared to other group 1 kimberlites in contrast the melilitites are sodic and are anomalously depleted in p compared to worldwide melilitites trace element patterns of the kola kimberlites and melilitites indicate the presence of k and p rich phases in the mantle source to account for the k troughs shown by both magma types a k rich phase such as phlogopite is thought to be residual in their sources however the anomalous k enrichment in the terskii coast kimberlites may indicate that an additional metasomatic k rich phase e g k richterite and or a complex k ba phosphate existed in the kimberlite source the p depletion in the melilitites may suggest that a phosphate phase such as apatite remained residual in the melilititic source however anomalous p enrichment in the kimberlites cannot be explained by complete melting of the same phase because the kimberlites are a smaller degree melt thus it is most likely that another metasomatic phosphate mineral existed in the source of the kimberlites the kola kimberlites and melilitites are all strongly lree enriched but the kimberlites have a steeper ree pattern and are significantly more depleted in hree indicating a higher proportion of garnet in their source higher nb y ratios and lower sio2 values in the kimberlites indicate that they were a smaller degree partial melt than the melilitites the presence of diamonds in the terskii coast kimberlites indicates a relatively deep origin while the melilitites originated from shallower depth the non diamondiferous kandalaksha monticellite kimberlite has lower abundances of all incompatible trace elements suggesting a higher degree of partial melting and or a less enriched and shallower source than the terskii coast kimberlites the 87sr 86sri 143nd 144ndi and pb isotope compositions confirm that the terskii coast kimberlites have close affinities with group 1 kimberlites and were derived from an asthenospheric mantle source while the kandalaksha monticellite kimberlite and terskii coast melilitites were derived from lithospheric mantle impact of a devonian asthenospheric mantle plume on the base of the archaean proterozoic lithosphere of the kola peninsula caused widespread emplacement of kimberlites melilitites ultramafic lamprophyres and other more fractionated alkaline magmas the nature of the mantle affected by metasomatism associated with the plume and in particular the depth of melting and the stability of the metasomatic phases gave rise to the observed differences between kimberlites and the related melilitites and other magmas

quantification of interplate coupling in subduction zones and forearc topography the effect of long term coupling in subduction zones on forearc topography is studied using a 2 dimensional finite element model a curved fault representing the interface between subducting slab and overriding plate has been incorporated into a model with an elastic crustal layer a viscoelastic upper mantle an initial topography and an effective coefficient of friction on the fault we assume that traction arising from friction is one of the stresses acting on topography for various friction coefficients the change of topography after some equilibration time is computed with specific geometry and kinematic boundary conditions for three subduction zones northern chile northern japan and tonga while the intrinsic coefficient of friction for small rock samples is high 0 6 0 8 the observed topography of forearc regions are inconsistent with an effective friction coefficient larger than 0 2

numerical modeling of stress and topography coupling during subduction inferences on global vs regional observables interpretation abstract the mechanical strength of the oceanic lithospheres exerts a fundamental control on the dynamics and tectonics of convergent margins here we use numerical models of coupled subducting upper plates and mantle to investigate the impact of the subducting lithosphere properties on the parameters commonly measured at subduction zones the models predict a correlation among parameters such as subducting and convergence velocities slab dip and curvature and stress coupling with the upper plates such signature emerges in a dynamic system through the slab buoyancy similar trends are found when assessing the impact of laboratory constrained non linear rheologies uncertainties on the mechanical strength and mantle flow although some features are more affected when the role of upper plate thickness is tested we find a positive correlation of upper plate subsidence and fore arc topography with the buoyancy of the down going plate however plasticity in the crust limits the development of topography becoming increasingly relevant with thinner upper plates so that no correlation can be established similar control is found on the width of the long wavelength dynamic subsidence due to flexure which increases non linearly with the upper plate thickness the models constrain relations among parameters within the dynamic system driven that allow inferences on the observable parameters the comparisons supports the idea that the subduction dynamics emerges as a global phenomenon in plate motions bending and dynamic subsidence yet andean type margin tectonics is best addressed in a regional context

the mechanical state of the lithosphere in the altiplano puna segment of the andes abstract information about topography the shape of the geoid seismicity neogene deformation and volcanism in the region of altiplano puna of western south america is used to analyse the state of stress across the convergent plate margin in terms of the effects of topography and simple models of its compensation an average elevation near 4 km is consistent with compensation by a yet unresolved combination of crustal root and hot uppermost mantle producing a geoid high of 22 27 meters average horizontal compressive stress in excess of a reference sea level lithostatic value of 390 bars in a 150 km thick lithosphere and an average shear stress of 170 bars along a 30 dipping interplate boundary the basis for these estimates is evidence for a neutral to extensional stress regime within the high plateau contrasted with a compressional regime on the eastern slopes and along the interplate boundary itself comparison with other plateaus in a convergent plate tectonic setting suggests an evolutionary sequence from compressional to extensional tectonics as elevation of the plateau increases

trench migration and subduction zone geometry seismic tomography results show a large variety in upper mantle structure along convergent plate boundaries we numerically investigate the effect of trench migration on the evolution of a slab with temperature dependent viscosity encountering a viscosity interface we find that subduction zone geometry is sensitive to even small rates 1 cm yr of retrograde motion increase in trench migration rate decreases the dip angle of the slab and its ability to penetrate the lower mantle upon including a background mantle flow it turns out that trench migration relative to the upper mantle flow is more decisive than the absolute plate velocities subduction zone geometry appears to be strongly time dependent we conclude that individual tectonic setting and time dependent slab behaviour can account for many different types of observable subduction zone anomalies our model temperature can adequately account for magnitudes and patterns of seismic anomalies as obtained from seismic tomography

petrology and geochemistry of xenoliths from the northern baltic shield evidence for partial melting and metasomatism in the lower crust beneath an archaean terrane lower crustal xenoliths entrained in a paleozoic ultramafic lamprophyre breccia pipe on elovy island kola peninsula russia represent some of the oldest lower crustal material yet investigated from europe the xenoliths vary from feldspar poor garnetrich rocks which resemble eclogites to feldspar rich garnet granulites quartz rich felsic granulites as well as pyroxenites and amphibole rich rocks are also present the mafic granulites eclogites represent a suite of gabbros and norites that is related by olivine fractionation the igneous protoliths may have formed in a manner analogous to lower crustal rocks from most other european xenolith localities i e by basaltic underplating but magmatic cumulates are not in evidence the kola lower crust was subjected to one or more metasomatic events which introduced up to 45 phlogopite and or amphibole into both eclogites granulites and pyroxenites the resulting rocks have strong enrichments in rb ba and k indicating that the lower crust is not uniformly depleted in lil and heat producing elements siliceous 65 sio2 and mafic 50 sio2 lithologies coexist in migmatitic xenoliths which provide evidence for partial melting processes and restite formation in mafic metaigneous lower crust the relationship if any between partial melting and metasomatism is unclear

seismicity and shape of the subducted nazca plate an updated compilation of earthquake locations and focal mechanism solutions from the international seismological centre and preliminary determination of earthquakes is the basis of a comprehensive study of the geometry of the wadati benioff zone beneath western south america the new data support previous mapping of a sharp flexure rather than a tear in the subducted nazca plate beneath southern peru and provide evidence for a similar flexure in the southward transition from nearly horizontal subduction to a slab with 30 dip at latitude 33 s in contrast the transition from 30 slab dip beneath bolivia to a nearly horizontal dip in the region between 28 s 32 s is more gradual occurring over several hundred kilometers of along strike distance between 20 s and 32 s this southward flattening corresponds to a broadening of a horizontal benchlike part of the subducted plate formed between 100 and 125 km depth the transition in continental tectonic style near 27 s 28 s from a wide volcanically active plateau to a narrow nonvolcanic cordillera appears not to be associated with the main slab flattening which begins to the north of these latitudes but with a more abrupt change in curvature of the subducted slab from convex upward to concave upward immediately below the plate boundary interface the concept of gaussian curvature is applied to slab bending to explain how subduction geometry is affected by the shape of the south american plate we hypothesize that the polarity of vertical curvature in the subducting slab is governed by the orientation of lateral curvature of the plate margin focal mechanism solutions for intermediate and deep earthquakes are grouped by geographic region and inverted for the orientation and relative magnitudes of the principal stresses results of the inversion indicate that downdip extension dominates in the slab above 350 km while downdip compression dominates at greater depths

a magmatic hydrothermal transition in arkaroola northern flinders ranges south australia from diopside titanite pegmatites to hematite quartz growth a set of palaeozoic diopside titanite veins are present in mesoproterozoic metagranites and metasediments that constitute the basement mt painter inlier of the adelaide fold belt south aus tralia these massive veins up to 1 m of pegmatitic nature contain large crystals of diopside lree y enriched titanite up to 40 cm in length and minor amounts of quartz they can be used to trace the sys tem s development from a high temperature magmatic stage through to a massive hydrothermal event the pegmatitic origin of these veins is evident from a complex fluid melt inclusion assemblage consisting of a highly saline inhomogeneous fluid and relicts of melt immiscibility of melt and heterogeneous highly saline fluids exceeding 61 eq mass nacl is preserved in primary inclusions in diopside and secondary inclusions in titanite indicating relatively shallow conditions of formation 510 20 c and 130 10 mpa graphic intergrowth of diopside and albite occurs at the contact with granitic pegmatites the system evolved into hydrothermal conditions which can be deduced from a later population of only fluid inclusions homogeneous and less saline 40 eq mass nacl trapped around 350 20 c and 80 10 mpa during quartz crystallization the conditions moved across the halite liquidus resulting in a heterogeneous mixture of brine and halite crystals which were trapped at 200 20 c and 50 10 mpa brecciation and a palaeo geother mal system overprinted the pegmatitic veins with an epithermal hematite quartz assemblage and lesser amounts of bladed calcite and fluorite in an intermit tently boiling hydrothermal system of fairly pure h2o at 100 140 c and 1 5 mpa remobilization of lrees and y from titanite and or the granitic host rock is evidenced by precipitation of apatite allanite and wakefieldite in an intermediate stage occasional incorporation of radioactive elements or minerals presumably u rich in the fluorite is responsible for radiolysis of h2 ot o h 2

overriding plate deformation and variability of fore arc deformation during subduction insight from geodynamic models and application to the calabria subduction zone in nature subducting slabs and overriding plate segments bordering subduction zones are generally embedded within larger plates such large plates can impose far field boundary conditions that influence the style of subduction and overriding plate deformation here we present dynamic laboratory models of progressive subduction in three dimensional space in which the far field boundary conditions at the trailing edges of the subducting plate sp and overriding plate op are varied four configurations are presented free both plates free sp fixed op fixed and sp op fixed we investigate their impact on the kinematics and dynamics of subduction particularly focusing on overriding plate deformation the results indicate that the variation in far field boundary conditions has an influence on the slab geometry subduction partitioning and trench migration partitioning our models also indicate that in natural narrow subduction zones assuming a homogeneous overriding plate the formation of back arc basins e g tyrrhenian sea aegean sea and scotia sea is generally expected to occur at a comparable location 250 700 km from the trench irrespective of the boundary condition in addition our models indicate that the style of fore arc deformation shortening or extension is influenced by the mobility of the overriding plate through controlling the force normal to the subduction zone interface trench suction our geodynamic model that uses the sp op fixed setup is comparable to the calabria subduction zone with respect to subduction kinematics slab geometry trench curvature and accretionary configuration furthermore the model can explain back arc and fore arc extension at the calabria subduction zone since the latest middle miocene as a consequence of subduction of the narrow calabrian slab and the immobility of the subducting african plate and overriding eurasian plate this setting induced strong trench suction driving fore arc extension and forced subduction to be accommodated almost entirely by slab rollback not trenchward subducting plate motion while trench retreat was accommodated almost entirely by back arc and fore arc extension not trenchward overriding plate motion comparable to our sp op fixed model

effects of plate bending and fault strength at subduction zones on plate dynamics for subduction to occur plates must bend and slide past overriding plates along fault zones because the lithosphere is strong significant energy is required for this deformation to occur energy that could otherwise be spent deforming the mantle we have developed a finite element representation of a subduction zone in which we parameterize the bending plate and the fault zone using a viscous rheology by increasing the effective viscosity of either the plate or the fault zone we can increase the rates of energy dissipation within these regions and thus decrease the velocity of a plate driven by a given slab buoyancy we have developed a simple physical theory that predicts this slowing by estimating a convecting cell s total energy balance while taking into account the energy required by inelastic deformation of the bending slab and shearing of the fault zone the energy required to bend the slab is proportional to the slab s viscosity and to the cube of the ratio of its thickness to its radius of curvature the distribution of dissipation among the mantle lithosphere and fault zone causes the speed of a plate to depend on its horizontal length scale using the observation that earth s plate velocities are not correlated to plate size we can constrain the lithosphere viscosity to be between 50 and 200 times the mantle viscosity with higher values required if the fault zone can support shear tractions 50 mpa over 300 km these subduction zone strengths imply that the mantle fault zone and lithosphere dissipate about 30 10 and 60 of a descending slab s potential energy release if the slab is 100 km thick the lithospheric component is highly dependent on slab thickness it is smaller for thin plates but may be large enough to prevent bending in slabs that can grow thicker than 100 km ubduction zone strength should be more stable than mantle viscosity to changes in mantle temperature so the controlling influence of subduction zones could serve to stabilize plate velocities over time as the earth cools because the details of convergent plate boundaries are so important to the dynamics of plate motion numerical models of mantle flow should treat subduction zones in a realistic way

current plate velocities relative to the hotspots incorporating the nuvel 1 global plate motion model nuvel 1 is a new global model of current relative plate velocities demets et al 1990 which differ significantly from those of prior models here we incorporate nuvel 1 into hs2 nuvel1 a new global model of plate velocities relative to the hotspots hs2 nuvel1 was determined from the hotspot data and errors used by minster and jordan 1978 to determine am1 2 which is their model of plate velocities relative to the hotspots am1 2 is consistent with minster and jordan s relative plate velocity model rm2 here we compare hs2 nuvel1 with am1 2 and examine how their differences relate to differences between nuvel 1 and rm2 hs2 nuvel1 plate velocities relative to the hotspots are mainly similar to those of am1 2 minor differences between the two models include the following 1 in hs2 nuvel1 the speed of the partly continental apparently non subducting indian plate is greater than that of the purely oceanic subducting nazca plate 2 in places the direction of motion of the african antarctic arabian australian caribbean cocos eurasian north american and south american plates differs between models by more than 10 3 in places the speed of the australian caribbean cocos indian and nazca plates differs between models by more than 8 mm yr although 27 of the 30 rm2 euler vectors differ with 95 confidence from those of nuvel 1 only the am1 2 arabia hotspot and india hotspot euler vectors differ with 95 confidence from those of hs2 nuvel1 thus substituting nuvel 1 for rm2 in the inversion for plate velocities relative to the hotspots changes few euler vectors significantly presumably because the uncertainty in the velocity of a plate relative to the hotspots is much greater than the uncertainty in its velocity relative to other plates

nature of the plate contact and subduction zones diversity abstract in recent studies we showed that the nature of the plate contact in subduction zones is an important physical feature in both oceanic lithospheric subduction and continental collision we investigated two fundamental states of the plate contact one based on a fault and the other based on a subduction channel using geodynamic modeling we determined the speci c signatures of both states of the subduction contact we established that the nature of the plate contact influences the dynamic response of the overriding and subducting plate and is one of the controlling factors whether back arc extension develops or not in the present study we combine results of our previous numerical experiments with a re analysis of published observations overall our synthesis connects seismic moment release with back arc deformation and tectonic processes at the margin it leads us to identify four classes of subduction zones the rst two classes result directly from our numerical experiments in class 1 subduction zones are characterized by a plate contact that is largely fault like with an accretionary margin in class 2 the plate contacts are largely channel type and have an erosive margin class 3 where the plate contact is entirely channel like consists of accretionary margins with a high sediment supply subduction zones of class 4 mostly characterized by an erosive convergent margin northern chile peru honshu and kuril are more complicated they can be explained by incorporating regional observations

metastability of gas hydrate marine occurrences of gas hydrate are normally confined to the top few hundred meters of sediments along deep continental margins the zone of stability for gas hydrate is limited in depth by increases in temperature below the seafloor we use thermodynamic calculations to show that gas hydrate can exist in a metastable state below the usual base of the stability zone we estimate that gas hydrate can be overheated by several degrees and that it may persist in this metastable state in the seafloor for as long as 106 years sudden decomposition of metastable hydrate should produce substantial pore pressure in the sediments contributing to slope failure in locations where gas hydrate is found such a mechanism might help to explain why slumping appears to be more frequent than average during the interval around the last glacial maximum

dynamic models of subduction geophysical and geological evidence in the tyrrhenian sea summary predictions based on a 2 d finite element model for subduction underneath the calabrian arc in southern italy are compared with a variety of geophysical and geological data such as the present day stress pattern within the slab uplift from the elevation of marine terraces in calabria and subsidence in the tyrrhenian marsili basin from odp leg 107 we model the behaviour of the slab driven by slab pull in agreement with the present tectonic style in this part of the mediterranean as suggested by several investigators the model accounts for the crustal lithospheric and mantle structures in a vertical cross section perpendicular to the calabrian subduction zone the shape of the slab is constrained on the basis of new tomographic images in the southern tyrrhenian sea which were obtained from the regional seismic stations of the istituto nazionale di geofisica while the rheological properties of the mantle are taken from global dynamic models density contrasts between the subducted slab and the surrounding mantle based on petrological models drive the flow in our viscoelastic model stress values displacements and vertical velocities at the surface are sampled at different times after loading until dynamic equilibrium is reached our estimates are appropriate for a time window of 100 kyr the validity of our comparison with the geological record is based on the assumption that the tectonic configuration in the past was not substantially different from that of the present day two families of models with unlocked and locked subduction faults are considered the unlocked models allow for roll back of the trench of about 20 mm yr in agreement with some geological estimates the same family of models predicts uplift of the calabrian arc of about 1 mm yr l and subsidence in the marsili basin of 1 2 mm yr in agreement with geological surveys the deviatoric stress obtained from the unlocked model is consistent with the continuous distribution of deep seismicity in the southern tyrrhenian sea with minor concentration within the lithospheric wedge locked models fail to reproduce these geophysical and geological observations predictions derived from a detached slab model are not consistent with the continuous hypocentral distribution of deep seismicity in the southern tyrrhenian sea deformation at the surface and the stress patterns at depth for a detached slab differ substantially from those of a continuous plate dynamic topography and horizontal motions are reduced when compared with the continuous plate with deviatoric stresses concentrated within the relict slab our results indicate that subduction is a major tectonic process in the southern tyrrhenian sea

two dimensional simulations of surface deformation caused by slab detachment detachment of the deeper part of subducted lithosphere causes changes in a subduction zone system which may be observed on the earth s surface constraints on the expected magnitudes of these surface effects can aid in the interpretation of geological observations near convergent plate margins where detachment is expected in this study we quantify surface deformation caused by detachment of subducted lithosphere we determine the range of displacement magnitudes which can be associated with slab detachment using numerical models the lithospheric plates in our models have an effective elastic thickness which provides an upper bound for rapid processes like slab detachment to the surface deformation of lithosphere with a more realistic rheology the surface topography which develops during subduction is compared with the topography shortly after detachment is imposed subduction with a non migrating trench system followed by detachment leads to a maximum surface uplift of 2 6 km while this may be higher for the case of roll back preceding detachment in the latter situation the back arc basin may experience a phase of compression after detachment within the context of our elastic model the surface uplift resulting from slab detachment is sensitive to the depth of detachment a change in friction on the subduction fault during detachment and viscous stresses generated by sinking of the detached part of the slab overall surface uplift of these magnitudes is not diagnostic of slab detachment since variations during ongoing subduction may result in similar vertical surface displacements d 2002 elsevier science b v all rights reserved

a 40ar 39ar and u pb isotopic study of the ilimaussaq complex south greenland implications for the 40k decay constant and for the duration of magmatic activity in a peralkaline complex abstract magmatism in the gardar province south greenland is related to two main rifting events at 1280 ma and 1180 to 1140 ma little is known about the duration of the magmatic activity in a specific complex the ilimaussaq intrusion belongs to the second period of rifting and comprises an extraordinary diversity of granitic and syenitic rock types which intruded and fractionated in three successive magmatic events as the intrusion contains some of the most evolved incompatible element rich rocks on earth it was chosen for a detailed geochronological study to quantify the duration of melt production intrusion fractionation and cooling amphiboles which are abundant in all rock types including pegmatites and late magmatic veins were dated using the 40 ar 39 ar technique since the solidus temperature of the most evolved melts is below the closure temperature of amphibole and no later heating event occurred the 40 ar 39 ar ages reflect the magmatic crystallisation and can be used to determine the duration of igneous differentiation the 40 ar 39 ar plateau ages range between 1142 6 2 2 ma and 1152 3 3 7 ma using the steiger and jager steiger r h and jager e 1977 subcommision on geochronology convention on the use of decay constants in geo and cosmochronology earth planet sci lett 36 359 362 40 k decay constant these ages are younger than an u pb age of 1160 5 ma on baddeleyite from the first magma batch our results indicate that the 40 k decay constant of steiger and jager may be too high a lambda similar to that proposed by kwon et al kwon j min k bickel p and renne p r 2002 statistical methods for jointly estimating decay constant of 40 k and age of a dating standard math geol 34 457 474 is required to make the 40 ar 39 ar match the u pb age

high rates of arc consumption by subduction processes some consequences new estimates of long term consumption of japanese and peruvian continental fore arc margins together with recent drilling data ocean drilling program on some other erosional island arcs such as izu bonin mariana or tonga allow the calculation of a mean rate of 7 3 km m y landward migration of the arc trench system on both east and west sides of the pacific basin during cenozoic time such relative arcward retreat counterbalances seaward migration of the slab hinge line to accommodate opening of the atlantic ocean retreat can thus be interpreted as a delay in seaward or retrograde migration due to upper mantle resistance to contraction of the pacific basin mass transfers in subduction zones must be reconsidered in the light of these new estimates especially the balance between net growth and loss of continental crust in recent geologic time a 50 m y period is long enough to erode the initial width of the volcanic arc and fore arc massif of any subduction zone assuming a mean trench volcanic arc landward migration of 5 km m y and a mean trench volcanic arc distance of 250 km this circumstance thus raises the question of whether the izu bonin mariana and tonga arcs are older than their reported eocene age

smnd evidence for late archaean crust formation in the lapland kola mobile belt kola peninsula russia and norway abstract a sm nd model age traverse across the kola peninsula has been undertaken to gain insight into the major episodes of crustal growth and to test the possible presence of early archaean saamian crust in the region granitoid gneisses and subordinate metasediments have been sampled from the main archaean geological units i e from northeast to southwest the murmansk domain model ages of 2 68 2 96 ga central kola domain including samples from the superdeep borehole 2 74 2 92 ga belomorian domain including the eastern part of the tersk area 2 70 2 93 ga the narrow range of depleted mantle model ages from 2 68 to 2 96 ga for archaean granitoid orthogneisses close to the oldest crystallization ages from the region indicate a major episode of late archaean crustal growth and reveal no indication of ancient saamian material in keeping with other geochronological evidence similar results 2 86 to 2 93 ga for archaean metasediments support this conclusion proterozoic model ages for supracrustals from the western part of the tersk area 2 16 to 2 23 ga indicate that the crust in this area is much younger than previously recognized initial nd values for late archaean samples with reasonable age constraints fall within a range of 2 3 to 2 5 two well dated samples from the belomorian domain have the highest initial nd values nd 2 81 ga 2 0 and nd 2 78 ga 2 3 and suggest that a depleted mantle source was present in the kola region in late archaean times dated samples from the central kola domain have lower initial nd between 2 5 and 1 2 suggesting a possible regional difference in mantle character in contrast to the depleted signature of the archaean orthogneisses between 2 7 ga and 2 8 ga previous studies of earliest palaeoproterozoic rift related majic intrusives in the kola and karelia regions report negative nd values at 2 4 2 5 ga possibly these data reflect a change to more enriched mantle values

experimental data for a proposed empirical igneous geobarometer aluminum in hornblende at 10 kbar pressure knowledge of the depth of solidification of talc alkalic plutons is critical to understanding and unraveling the complex spatial temporal and chemical evolution of orogenic belts an empirical geobarometer has recently been proposed based on the total al content alt of near solidus calcic hornblendes in tonalite and granodiorite intrusive rocks in the pressure range 1 to 8 kbar pressures were estimated from mineral assemblages in associated metamorphic country rocks we report experimental data on newly crystallized near solidus amphiboles in partly melted vapor absent garnet tonalite at 10 kbar that provide some experimental support for the proposed geobarometer consistent with the extrapolated empirical result and confirm its limited applicability to specific multimineral assemblages

a three dimensional model of the early precambrian crust under the southeastern fennoscandian shield karelia craton and belomorian tectonic province abstract the early precambrian structures of the archaean karelia craton and belomorian tectonic province have been investigated by reflection seismology integration of the existing geological database for these two entities and adjacent tectonic units with the seismic data for the crust and upper mantle along the 1 eu geotraverse and cross traverse 4b totaling c 1300 km allows development of a three dimensional model of the deep structure of the southeastern fennoscandian shield interpretation of the seismic data provides evidence that the early precambrian crust is characterized by inclined layering the crustal architecture is a combination of thrust nappe and thrust underthrust structural assemblages consisting of neoarchaean and palaeoproterozoic rocks that were subsequently further deformed during uplift and emplacement of granite gneiss domes the new model substantially modifies previous conceptions of this part of the shield as being constituted of blocks with sub vertical margins and individual internal layering it is apparent that the surface traces of gently dipping boundaries are defined by the erosion level and consequently these boundaries are not those of particular tectonic blocks in the traditional sense

narrow subducting slabs and the origin of backarc basins abstract both tectonic and seismic data suggest that horizontally wide in the direction parallel to the axis of the trench dipping subducting slabs typically move stably at fixed angles to the horizontal and that their associated trenches migrate very slowly or not at all in contrast where backarc basins babs exist the associated trenches typically migrate at rates of several cm year this behavior may be related to the finite width of the slabs under babs analysis shows that hydrodynamic suction created by the downdip motion of a wide slab is sufficient to balance the pull of gravity on the heavy slab however for a narrow slab sideways asthenospheric flow into the slab overriding plate corner effectively eliminates the suction which results in the rapid sinking of the slab hydrodynamic modeling is used to quantify this hypothesis modeling is used also to show that wide slabs are stable in their inclined position not because of the gravitational wing like sliding but due to the lifting force produced by the negative pressure suction in the slab overriding plate corner finally the model results confirm that for relatively narrow slabs this lifting force is small and cannot counterbalance gravity therefore when bending moments in a narrow slab exceed the plastic yield limit the slab becomes unstable this instability can cause the rapid migration of the slab s hinge and the associated arc at rates of 10 20 cm year and thus form the bab

effects of slab mineralogy on subduction rates although velocities of subducting slabs should be controlled primarily by their negative thermal buoyan cies their mineralogy can also have signicant eects we explore this by using thermo kinetic modeling to predict mineralogy and compare the resultant buoyant driving force to the opposing viscous drag phase transitions of mg fe 2sio4 in subducting slabs depend on thermal struc ture in two ways first equilibrium phase boundaries should be deflected causing local buoyancy anomalies whose sign depends upon that of the clapeyron slope as slabs rst enter the transition zone negative anomalies should acceler ate them but positive anomalies that arise when they fully penetrate the transition zone should slow them such ef fects may induce geologically abrupt changes in plate mo tions second olivine that persists metastably in slabs will form regions of positive buoyancy which should reduce slab velocities the coldest and fastest slabs should be slowed more greatly thus narrowing the range of feasible subduc tion rates decreased descent rates however allow slabs to warm and metastable wedges to thermally erode such neg ative feedback mechanisms may serve to regulate subduction rates

the tectonic expression slab pull at continental convergent boundaries examination of five thrust belt systems developed at continental subduction boundaries suggests that they comprise two distinct groups that display pronounced and systematic differences in structural style topographic elevation denudation metamorphism postcollisional convergence and foredeep basin geometry and facies the distinctive geological features developed within each thrust belt group appear to be causally linked to the relative rates of subduction and convergence via the magnitude of horizontal compressional stress transmitted across the subduction boundary at subduction boundaries where the rate of overall plate convergence is less than the rate of subduction termed here retreating subduction boundaries the transmission of horizontal compressive stress across the plate boundary is small and regional deformation of the overriding plate is by horizontal extension the tectonic expression of these retreating subduction boundaries includes topographically low mountains little erosion or denudation low grade to no metamorphism little to no involvement of crystalline basement in shortening little to no postcollisional convergence anomalously deep foredeep basins and a protracted history of flysch deposition within the adjacent foredeep basin analysis of deflection and gravity data across three retreating subduction boundaries apennine carpathian and hellenic systems shows that subduction is driven by gravitational forces acting on dense subducted slabs at depths between about 40 and 80 km carpathians 50 and 150 km apennines and 50 and 250 km hellenides the total mass anomalies represented by the slabs are approximately 3 1012 6 1012 and 12 1012 n m respectively the slabs are partially supported by flexural stresses transmitted through the subducted lithosphere to the foreland and partially supported by dynamic viscous stresses in the asthenosphere at subduction boundaries where the rate of overall plate convergence is greater than the rate of subduction termed here advancing subduction boundaries the transmission of horizontal compressive stress across the plate boundary is large and regional deformation of the overriding plate is by horizontal shortening the tectonic expression of these advancing subduction boundaries includes topographically high mountains antithetic thrust belts large amounts of erosion and denudation exposure of high grade metamorphic rocks at the surface extensive deformation of crystalline basement to midcrustal depths protracted postcollisional convergence tens of millions of years and a protracted history of molasse deposition within the adjacent foredeep basins analysis of gravity and deflection data across two advancing subduction boundaries developed within the continental lithosphere western to eastern and southern alps and himalayas shows that the thrust sheets have been translated for great distances over the foreland lithosphere relative to the point at which the subduction forces are applied thus obscuring any flexural and gravity signals from the subducted slab however it appears that far field stresses presumably related to global plate motions drive most of the convergent motion across these subduction boundaries the concept that orogenic belts formed above retreating subduction boundaries have recognizable tectonic signatures that differ from those of orogenic belts formed above advancing subduction boundaries suggests that it may be possible to interpret the plate boundary settings in which ancient orogenic belts evolved appendix b is available with entire article on microfiche order from the american geophysical union 2000 floridaavenue n w washington d c 20009 document t92 004 2 50 payment must accompany order

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

stresses and plate boundary forces associated with subduction plate margins primary tectonic stress in the lithosphere is predominantly caused by lateral density variations within the earth and associated topographical loading when such a stress system caused by major sublithospheric density anomalies is intersected by a weak zone which cuts across the lithosphere plate boundary forces develop and modify the plate interior stresses in this paper unite element analysis is used to model the stresses and plate boundary forces associated with subduction plate boundaries with dipping and vertical slabs extending to about 270 and 400 km depths having the subduction fault both locked and unlocked in such regions several types of horizontal deviatoric stress may occur including 1 local compression in the trench arc region caused by the dense sinking slab and the associated surface downflexure 2 plate interior tension which occurs when this compression is intersected by a weak subduction fault 3 local tension associated with thickening of the crust at the arc and elsewhere 4 local tension in the back arc region produced by the underlying low density upper mantle and 5 downbending stresses in the subducting slab thermal stresses and transmitted ridge push which are not included in the modelling here a gradient from compression in the forearc to tension in the back arc can be modelled in terms of these stresses when the fault is partially locked it is however the intersection of the local compression 1 by an unlocked or partially locked subduction fault that modifies the plate interior stresses and gives rise to the slab pull and trench suction plate boundary forces the state of stress in the interior of the overriding plate is also crucially influenced by back arc spreading where this occurs plate boundary forces have been evaluated for each of the models it is shown that slab pull and trench suction may be of comparable magnitude

increased continental margin slumping frequency during sea level lowstands above gas hydrate bearing sediments we present sup 14 c data on sediment samples from cores of the upper 7 m of the sediment column overlying a major continental rise gas hydrate field on the southern carolina rise and inner blake ridge offshore the southeastern united states the data show that glacial age sediments are underrepresented in the cores the observation is consistent with a previously predicted association between sea level lowstands and increased frequency of sea floor slumping on continental margins containing gas hydrates 26 refs 3 figs

controls of subduction geometry location of magmatic arcs and tectonics of arc and back arc regions most variation in geometry and angle of inclination of subducted oceanic lithosphere is caused by four interdependent factors combinations of 1 rapid absolute upper plate motion toward the trench and active overriding of the subducted plate 2 rapid relative plate convergence and 3 subduction of intraplate island seamount chains aseismic ridges and oceanic plateaus anomalously low density oceanic lithosphere cause low angle subduction under conditions of low angle subduction the upper surface of the subducted plate is in contact with the base of the overlying plate the wedge of low density asthenosphere is replaced by subducted lithosphere and the width of the arc trench gap either is significantly increased or a magmatic arc is not developed within the overlying plate the fourth factor is age of the subducting lithosphere subduction of young lithosphere produces two opposing tendencies 1 low angle subduction and increased arc trench distance owing to its low density and 2 decreased arc trench distance owing to its higher temperature two factors of secondary importance contribute to variation in subduction zone geometry and arc trench distance accretion of sediment in trenches depresses the upper portion of the subducting oceanic plate and causes the trench axis to migrate seaward prolonged subduction thickens the upper plate depresses the isotherms in the subducted plate and may create a broader arc both factors increase the arc trench gap the four primary factors also control development of other tectonic elements such as regional subsidence for example the amazon basin and a portion of the cretaceous interior seaway of western united states intra arc extension for example the basin and range province foreland fold and thrust belts and laramide style tectonics

the development of slabs in the upper mantle insights from numerical and laboratory experiments we have performed numerical and laboratory experiments to model subduction of oceanic lithosphere in the upper mantle from its beginnings as a gravitational instability to the fully developed slab a two dimensional finite element code is applied to model newtonian creep in the numerical experiments scaled analog media are used in the laboratory a sand mixture models the brittle crust silicone putty simulates creep in the lower crust and mantle lithosphere and glucose syrup is the asthenosphere analog both model approaches show similar results and reproduce first order observations of the subduction process in nature based on density and viscosity heterogeneities in a stokes flow model subduction nucleates slowly and a pronounced slab forms only when the viscosity contrast between oceanic plate and mantle is below a threshold we find that the subduction velocity and angle are time dependent and increase roughly exponentially over tens of millions of years before the slab reaches the 670 km discontinuity the style of subduction is controlled by the prescribed velocity of convergence the density contrast between the plates and the viscosity contrast between the oceanic plate and the mantle these factors can be combined in the buoyancy number f which expresses the ratio between driving slab pull and resisting viscous dissipation in the oceanic plate variations in f control the stress in the plates the speed and the dip of subduction and the rate of trench retreat reproducing the contrasting styles of subduction observed in nature the subduction rate is strongly influenced by the work of bending the lithosphere as it subducts

holocene emerged notches and tectonic uplift along the jalisco coast southwest mexico this paper presents the preliminary results from a study of holocene emerged shorelines marine notches and their tectonic implications along the jalisco coast the pacific coast of jalisco sw mexico is an active tectonic margin this coast has been the site of two of the largest earthquakes to occur in mexico this century the 1932 mw 8 2 jalisco earthquake and the 1995 mw 8 0 colima earthquake measurement and preliminary radiocarbon dating of emergent paleoshorelines along the jalisco coast provide the first constraints upon the timing for tectonic uplift along this coastline uplifted holocene marine notches and wave cut platforms occur at elevations ranging from ca 1 to 4 5 m amsl in situ intertidal organisms dated with radiocarbon the first ever reported for the jalisco area provide preliminary results that record tectonic uplift during at least the past 1300 years bp at an average rate of about 3 mm year we propose a model in which coseismic subsidence produced by offshore earthquakes is rapidly recovered during the postseismic and interseismic periods the long term period is characterized by slow tectonic uplift of the jalisco coast we found no evidence of coastal interseismic and long term subsidence along the jalisco coast

buckling of a pervasively faulted lithosphere we investigated the buckling response of a faulted elastic plate under horizontal compression using the finite element technique to better understand the effect of faults on the elastic behavior of a plate we studied the effect of changes in fault spacing depth and dip on the effective young s modulus buckling stress and wavelength our model consists of a thick elastic plate whose entire upper surface is cut by evenly spaced faults we impose either an initial sinusoidal deformation with a fixed wavelength or a random deformation to the grid a fault is represented as a free surface with no resolved shear stress and is allowed to slip in a specified direction using the method of slippery nodes with the assumption of free slip on the faults our model results represent an end member case in which the buckling wavelength and buckling stress are minimized by the presence of the faults

role of naturally occurring gas hydrates in sediment transport naturally occurring gas hydrates have the potential to store enormous volumes of both gas and water in semi solid form in ocean bottom sediments and then to release that gas and water when the hydrate s equilibrium conditions are disturbed therefore hydrates provide a potential mechanism for transporting large volumes of sediments under the combined low bottom water temperatures and moderate hydrostatic pressures that exist over most of the continental slopes and all of the continental rises and abyssal plains hydrocarbon gases at or near saturation in the interstitial waters of the near bottom sediments will form hydrates the gas can either be autochthonous microbially produced gas or allochthonous catagenic gas from deeper sediments equilibrium conditions that stabilize hydrated sediments may be disturbed for example by continued sedimentation or by lowering of sea level in either case some of the solid gas water matrix decomposes released gas and water volume exceeds the volume occupied by the hydrate so the internal pressure rises drastically if large volumes of hydrate are decomposed part of the once rigid sediment is converted to a gas and water rich relatively low density mud when the internal pressure due to the presence of the compressed gas or to buoyancy is sufficiently high the overlying sediment may be lifted and or breached and the less dense gas cut mud may break through such hydrate related phenomena can cause mud diapirs mud volcanos mud slides or turbidite flows depending on s diment configuration and bottom topography

mechanisms for the origin of mid ocean ridge axial topography implications for the thermal and mechanical structure of accreting plate boundaries to better understand the implications of ridge axis morphology for the thermal and mechanical structure of accreting plate boundaries we examine mechanisms for the creation of stress supported axial topography and its dependence on parameters such as spreading rate asthenosphere viscosity and plate thickness near the ridge axis three basic mechanisms are considered 1 asthenospheric upwelling in a narrow steady state conduit beneath a ridge axis 2 mantle flow stresses due to plate spreading extending ridge axis lithosphere or melt segregation and associated mantle compaction and 3 topography supported by moments due to horizontal extensional stresses in a thickening lithosphere if mid ocean ridge topography is due to asthenosphere rising in a narrow conduit at the ridge axis then the conduit structure cannot be thermally controlled as the associated thermal structure has relatively broad flat isotherms instead of a conduitlike shape mantle flow associated with plate spreading alone induces no stress supported topography thus the intuitive notion that diverging flow at the ridge axis due solely to plate spreading produces the axial valley is incorrect vertical mantle flow at the base of a stretching accreting lithosphere and mantle flow associated with melt segregation can however lead to stress supported axial topography however for these mantle flow related processes to produce significant amounts of axial relief the mantle viscosity beneath the ridge axis must be on the order of 1020 pa s 1021 p poise significantly higher than generally accepted values of 1018 pa s 1019 p horizontal extensional stresses in a strong brittle lithosphere that thickens with distance from the ridge axis can produce axial topography of the form and amplitude observed on slow spreading ridges a continuum idealization of deformation due to faulting is used to examine this mechanism this model predicts the axial valley width is controlled by the plate thickness near the ridge axis an 8 km thick plate near the ridge axis would produce the 30 km wide axial valley observed at slow spreading ridges an increased lithosphere thickness of only several kilometers over the half width of the axial valley can produce kilometer scale relief to examine whether a plate 8 10 km thick at the ridge axis that thickens by several kilometers within 15 km of the axis is plausible at a slow spreading ridge thermal models that include both the heat of magmatic crustal accretion on the ridge axis and upwelling mantle flow beneath the axis are formulated hydrothermal cooling that increases the effective conductivity by about a factor of 10 is required for an 8 km thick plate to exist at a slow spreading ridge axis however even with this enhanced hydrothermal cooling the plate thickness at a fast spreading ridge axis and the resulting horizontal force in the plate are too small to result in appreciable topography while mid ocean ridge axial topography is strongly influenced by spreading rate other factors also shape axial structure for example the anomalously shallow slow spreading reykjanes ridge has an axial structure more like that of a fast spreading ridge and the anomalously deep intermediate spreading australian antarctic discordance has a median valley axial structure typical of a slower spreading ridge in addition within a given spreading segment there is often a systematic deepening of the ridge axis toward transform offsets these non spreading rate dependent axial topographic variations may be produced by variations in magmatic heat input associated with crustal emplacement which shapes the thermal structure and lithosphere thickness near the ridge axis

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

study of plate deformation and stress in subduction processes using two dimensional numerical models two dimensional finite element modeling is used to model subduction of an oceanic lithospheric plate beneath continental lithosphere the subduction process is initiated along a preexisting inclined fault and continues until reaching 400 km of total convergence the lithosphere is assumed to be underlain by an in viscid asthenosphere different rheological laws have been considered for the lithosphere including elasticity and elastoplasticity the modeling shows that both the stress system in the plates and the surface topography are strongly dependent on two main parameters the density contrast between lithosphere and asthenosphere l a and the coefficient of friction along the subduction plane varying these two parameters allows explanation of the main characteristics of real subduction zones and results in two major regimes manifested by extension or compression in the arc back arc system extension and back arc rifting corresponds to a positive density contrast and a low coefficient of friction while negative values and or high friction leads to a compressional regime the coexistence of trench arc compression and back arc tension is only possible for a coefficient of friction lower than 0 1 the results of the numerical experiments agree with those of experimental modeling conducted under similar physical assumptions

dynamics of retreating slabs 1 insights from two dimensional numerical experiments from the effects of mantle flux part 2 therefore in this paper we apply forces to the slab using simple analytical functions related to buoyancy and viscous forces in order to isolate the role of rheology on slab dynamics we analyze parameters for simplified elastic viscous and nonlinear viscoelastoplastic single layer models of slabs and compare them with a stratified thermomechanical viscoelastoplastic slab embedded in a thermal solution the near surface behavior of slabs is summarized by assessing the amplitude and wavelength of forebulge uplift for each rheology in the complete thermomechanical solutions vastly contrasting styles of slab dynamics and force balance are observed at top and bottom bends however we find that slab subduction can be modeled using simplified rheologies characterized by a narrow range of selected benchmark parameters the best fit linear viscosity ranges between 5 10 22 pa s and 5 10 23 pa s the closeness of the numerical solution to nature can be characterized by a deborah number 0 5 indicating that elasticity is an important ingredient in subduction index terms 8120 tectonophysics dynamics of lithosphere and mantle general 8159 tectonophysics rheology crust and lithosphere 8160 tectonophysics rheology general keywords subduction numerical models lithospheric rheology

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

thermal and visual time series at a seafloor gas hydrate deposit on the gulf of mexico slope ambient temperature is a critical factor determining the stability of gas hydrate deposits on continental margins to study this process directly under varying conditions a monitoring array comprising a time lapse camera and in situ temperature probes was deployed at a hydrocarbon seep known as bush hill where gas hydrates deposits are exposed at the seafloor in a water depth of 570 m for intervals of 91 days and 30 days the digital camera recorded several daily images of a prominent gas hydrate mound consisting of structure ii gas hydrate the temperature probes were constructed with one autonomous thermistor at each end of a 50 cm pvc wand and recorded temperatures with precision of better than 0 1 8c at 30 min intervals over 327 days one probe was implanted with a tight seal into a drill hole about 7 cm deep in the top of the gas hydrate mound the second was inserted about 50 cm deep into the adjacent sediments for each probe the top thermistor recorded the ambient water temperature while the bottom thermistor synoptically recorded the internal temperature of the hydrate or sediment the bottom water temperatures ranged from 6 64 to 9 73 8c with a mean of 7 90 8c and standard deviation of 0 437 photographic results showed no dramatic changes in the size shape or gas venting from the mound by comparing the temperatures recorded at the tips of the probes with the synoptic water temperature we estimate that thermal diffusivity for the gas hydrate deposit was 7 14d 10 8 to 8 33d 10 8 m 2 s 1 and was 1 77d 10 7 to 3 01d 10 7 m 2 s 1 for the sediments at the sites where temperatures were recorded the diffusivity measured in gas hydrate was lower than that measured under laboratory conditions stability of gas hydrate in this setting is not likely to be affected by short term changes in bottom water temperature within the range observed d 2005 elsevier b v all rights reserved

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

quantifying the energy dissipation of overriding plate deformation in three dimensional subduction models in a subduction system the force and the energy required to deform the overriding plate are generally thought to come from the negative buoyancy of the subducted slab and its potential energy respectively such deformation might involve extension and back arc basin formation or shortening and mountain building how much of the slab s potential energy is consumed during overriding plate deformation remains unknown in this work we present dynamic three dimensional laboratory experiments of progressive subduction with an overriding plate to quantify the force fopd that drives overriding plate deformation and the associated energy dissipation rate opd and we compare them with the negative buoyancy fbu of the subducted slab and its total potential energy release rate bu respectively we varied the viscosity ratio between the plates and the sublithospheric upper mantle with sp um 157 560 and the thickness of the overriding plate with top 0 5 2 5 cm scaling to 25 125 km in nature the results show that fopd fbu has average values of 0 5 2 0 with a maximum of 5 3 and opd bu has average values of 0 05 0 30 with a maximum of 0 41 the results indicate that only a small portion of the negative buoyancy of the slab and its potential energy are used to deform the overriding plate our models also suggest that the force required to deform the overriding plate is of comparable magnitude as the ridge push force furthermore we show that in subduction models with an overriding plate bending dissipation at the subduction zone hinge remains low 3 15 during steady state subduction

rates of neogene and quaternary tectonic movements in the southern banda arc based on micropalaeontology it is generally considered that the outer banda arc islands are composed in part of strongly deformed sediments that accumulated at the north australian continental margin the onset and completion of folding nappe emplacement imbrication and uplift that produced the islands of the outer banda arc have been dated by means of planktonic foraminifera the dating of the onset of uplift of timor from a deep submarine position at the end of the neogene nappe emplacement to mountains now reaching nearly 3 km high indicates the rate of post collision uplift to have been initially 3 mm yr then it slowed to about 1 5mm yr where the australian continental margin meets the eastern end of the present java trench there are indications that the australian margin has overridden the trench in the timor region by 240 km micropalaeontological dating of the neogene collision between the australian continental slope proximal rise and the trench provides constraints for checking the overriding hypothesis by calculating rates of convergence derived from plate motions 7 5cm yr against the rates based on the tectonic model 7 4cm yr after nappe emplacement possibly faster than 6 25 cm yr shortening of the continental crust appears to have migrated towards the australian continent so that the shelf itself has become involved in the imbrication with shortening of the cover rocks between the nappes and the present shelf edge amounting to about 40 km during the last 2 ma

thermobarometry and geotectonic significance of high pressure granulites examples from the moldanubian zone of the bohemian massif in lower austria

mechanics of graben formation in crustal rocks a finite element analysis small 1 10 km wide grabens are ubiquitous tectonic structures that have been observed on the surfaces of all the terrestrial planets and satellites in spite of their frequent occurrence there is no fundamental understanding of how they form other than that they are characteristic of extensional terranes this paper examines the mechanics of the initial stages of graben formation and demonstrates that the configuration of a pair of antithetically dipping normal faults a graben is the energetically most favorable fault configuration in elastic brittle rocks subjected to pure extension a corollary to this hypothesis is that the occurrence of synthetic normal faults is indicative of a stress field other than purely extensional or of nonhomogeneous characteristics of the rock the analysis proceeds by computing the stress field in the vicinity of a single initial normal fault using a two dimensional finite element model it is shown that slip on this fault substantially lowers the proximity to failure of the surrounding rocks except for two small regions one near the surface of the downthrown block and the other near the fault tip the region on the downthrown block is located a horizontal distance from the initial fault slightly less than the fault s depth the proximity to failure actually increases here and a secondary fault develops an antithetic secondary fault would link this region with the other failure zone near the fault tip to demonstrate the preference for antithetic secondary faults more rigorously we show that this fault configuration results in a smaller total strain energy than the synthetic configuration in addition a two layer model is used to show that the graben s width is only weakly dependent upon the presence of mechanically distinct layers in the crust both weak over strong and strong over weak configurations are investigated the depth of the initial fault and the thickness of the top layer are varied and the width of the resulting graben is evaluated by computing the distance between the initial fault and the region of highest proximity to failure near the surface this analysis indicates that the depth of the initial fault is the major factor controlling graben width

subduction of the lithosphere and back arc dynamics insights from physical modeling physical modeling of the subduction of an elastico plastic lithosphere underlained by a low viscosity asthenosphere and driven by both horizontal compressional force and gravitational sinking of the subducting plate has shown the subduction regime to be strongly affected by 1 a where a and 1 are the asthenosphere and lithosphere densities at 0 a steady state process develops which requires a horizontal lithosphere compression of h 0 3 0 4 s 108 pa where s is the lithosphere yield limit the central part of the leading wedge of the overriding plate is nonisostatically elevated in this case especially in the region of the frontal arc which corresponds to a free air gravity anomaly maximum of 200 400 mgal a decrease in 0 compression of the overriding plate is considerably less from the onset of subduction and then decreases during subduction at some point compression changes into tension tensional nonhydrostatic stresses grow during the process to a value near h 0 1 0 2 s 0 4 108 pa both subsidence of the overriding wedge below the isostatic equilibrium level and an increase in the benioif dip angle accompany this process the gravity anomaly maximum is displaced into the volcanic arc in this subduction regime where it is normally less than 150 mgal if some weaknesses are present in the overriding plate tensional stresses result in a failure of the lithosphere at these locations and the broken lithosphere segment the arc plate starts to move counter to the subduction plate the arc plate slides along the underthrusting plate s surface without loss of contact and causes back arc opening hydrostatic suction prevents the plates from being separated this effect is caused by the peculiarity of a mechanical system whereby a low viscosity asthenosphere and a plastic lithosphere which effective yield limit is considerably lower than the hydrostatic stresses in the plate if the overriding lithosphere is not weakened enough to fail under the existing tension then the hanging slab breaks down under its own weight when it reaches a specific length both interacting plates undergo isostatic recovery in this case and the whole cycle is repeated again the hydrostatic suction effect also provides back arc tension under other conditions such as when asthenospheric flow forces the subducted plate to move away from the island arc

capillary effects on hydrate stability in marine sediments 1 we study the three phase liquid gas hydrate stability of the methane hydrate system in marine sediments by considering the capillary effects on both hydrate and free gas phases the capillary pressure a measure of the pressure difference across a curved phase interface exerts a key control on the methane solubility in liquid hydrate l h and liquid gas l g systems by calculating the l h and l g solubilities as a function of water depth pressure and pore size interface curvature we show how the solubility requirements for forming both gas hydrate and free gas can be met in a three phase zone the top of the three phase zone shifts upward in sediments as the water depth increases and the mean pore size decreases the thickness of the three phase zone increases as the distribution of pore sizes widens the top of the three phase zone can overlie or underlie the bulk three phase equilibrium depth at blake ridge we predict that the three phase zone is 27 7 m thick and that the top of the three phase zone lies 13 m above the predicted bulk equilibrium depth this reconciles the observation of the bottom simulating reflector bsr at blake ridge that is shallower than the predicted bulk equilibrium depth in contrast at hydrate ridge where water depth is shallower we predict that the three phase zone is 20 4 m thick and that the top of the three phase zone lies 0 7 m below the predicted bulk equilibrium depth our model which predicts an upward shift in the top of free gas occurrence with increasing water depth pressure is compatible with worldwide observations that the bsr is systematically shifted upward relative to the bulk equilibrium depth as water depth pressure is increased

the end of the ice age methane gas hydrates may have contributed to the rapid rise in atmospheric ch4 co2 and global temperatures at the end of the last major glaciation about 13 500 years ago given suitable orbital conditions and insolation at high latitude a small triggering event such as the release of one or more arctic gas pools may have initiated massive release of methane from hydrate under ice and in permafrost the consequent greenhouse warming would have provided strong positive feedback amplifying emission this warming driven by ch4 may have induced the release of co2 from the oceans to the biosphere stabilizing the interglacial carbon cycle

global trench migration velocities and slab migration induced upper mantle volume fluxes constraints to find an earth reference frame based on minimizing viscous dissipation sincetheadventofplatetectonicsdifferentglobalreferenceframeshavebeenusedtodescribethemotionofplatesandtrenches thedifferencein plate motion and trench migration between different reference frames can be substantial up to 4 cm yr this study presents an overview of trench migration velocities for all the mature and incipient subduction zones on earth as calculated in eight different global reference frames calculations show that irrespective of the reference frame 1 trench retreat always dominates over trench advance with 62 78 of the 244 trench segments retreating 2 the mean and median trench velocity are always positive retreating and within the range 1 3 1 5 cm yr and 0 9 1 3 cm yr respectively 3 rapid trench retreat is only observed close to lateral slab edges b1500 km and 4 trench retreat is always slow far from slab edges n2000km thesecalculationsarepredictedbygeodynamicmodelswithavaryingslabwidth inwhichplatemotion trenchmotionandmantleflow result from subduction ofdense slabs suggesting that trenchmotion is indeedprimarily drivenby slab buoyancyforces andthat proximity to alateral slab edge exerts a dominant control on the trench migration velocity despite these four general conclusions significant differences in velocities between such reference frames remain it is therefore important to determine which reference frame most likely describes the true absolute velocities togetanunderstandingoftheforcesdrivingplatetectonicsandmantleconvection itishereproposedthat basedonfluiddynamicconsiderationsand predictions from geodynamic modelling the best candidate is the one which optimises the number of trench segments that retreat minimizes the trench perpendicular trench migration velocity vt in the centre of wide n4000 km subduction zones maximizes the number of retreating trench segments located within 2000 km of the closest lateral slab edge minimizes the average of the absolute of the trench perpendicular trench migration velocity vt forallsubductionzonesonearth andminimizestheglobaluppermantletoroidalvolumeflux to thatresultsfromtrenchmigration and associated lateral slab migration i e slab rollback or slab roll forward calculations show that these conditions are best met in one particular indo atlantic hotspot reference frame where 75 of the subduction zones retreat vt in the centre of wide subduction zones ranges between 3 5 and 1 8 cm yr 83 of the trench segments located within 2000 km of the closest lateral slab edge retreat the average of vt is 2 1 cm yr and to 456 km 3 yr lower limit and 539 km 3 yr upper limit inclusion of all the incipient subduction zones on earth results in slightly greater fluxes of 465 km 3 yr lower limit and 569 km 3 yr upper limit it is also found that this reference frame is close to minimizing the total sub lithospheric uppermantlevolumeflux k associatedwithmotionofcontinentalkeelslocatedbelowthemajorcratons itisstressed however that kisanorder of magnitude smaller than to and thus of subordinate importance in conclusion the indo atlantic hotspot reference frame appears preferable for calculating plate velocities and plate boundary velocities

migration of methane gas through the hydrate stability zone in a low flux hydrate province new high resolution seismic data show clear evidence for upward injection of methane gas well into the hydrate stability zone at the stable low methane flux blake ridge crest this movement of gaseous methane through a thermo dynamic regime where it should be trapped as hydrate suggests that dynamic migrations of gas play an important role in the interaction of subseafloor methane with the ocean in the study area none of the seismic amplitude anomalies that provide evidence for gas migration reaches the seafloor instead they terminate at the base of a highly reflective unfaulted capping layer seismic inversions of anomalous regions show 1 increased velocities beneath the hydrate stability zone suggesting less gas and 2 increased velocities within the hydrate stability zone associated with observed low amplitude chimneys and bright spots indicating increased hydrate concentrations these observations and analyses indicate that methane migrates upward as free gas hundreds of meters into the hydrate stability zone before forming hydrate the observations strongly imply that given appropriate permeable pathways free gas can escape into the ocean even in a low flux environment the hydrate stability zone is not an impermeable barrier to free gas migration

molecular and stable isotope compositions of natural gas hydrates a revised global dataset and basic interpretations in the context of geological settings abstract a global dataset of molecular and stable isotope compositions of gases released from 209 different specimens of natural gas hydrate is presented and discussed the 26 hydrate bearing areas from 21 geographic regions are grouped into high gas flux hgf settings low gas flux lgf settings and hydrated gas accumulations hga methane ch4 is the most abundant hydrate bound gas while co2 and c2 hydrocarbon gases are frequently present in small amounts non hydrocarbon gases such as h2s are uncommon the stable isotope composition of hydrate bound gases varies significantly e g 13c of ch4 from 74 7 to 39 6 suggesting that gases of both microbial and thermogenic origin form gas hydrates hydrate bound gases are derived from autochthonous located predominantly within the gas hydrate stability zone ghsz and allochthonous located in deep sediments sources the occurrence and concentration of gas hydrates in sediments are controlled not by the origin of gases but rather by their sources which in turn strongly depend on geological setting allochthonous gases microbial and or thermogenic dominate in hgf and hga settings where they are focused in the shallow ghsz along faults within mud volcanoes in permeable carrier beds and other geological features from underlying petroleum systems relatively high concentrations of gas hydrate average 5 15 of pores and locally up to 100 of volume occurring over small areas are typical of hgf and hga settings in contrast autochthonous and diffuse allochthonous gases mostly microbial occur in stratigraphically and structurally simple lgf settings and result in relatively low concentrations of gas hydrate average 2 of pores in the ghsz spread over large areas the major implication of this finding is that successful prediction of resource and geohazard potential of gas hydrates is possible only if regional petroleum systems extending well below the ghsz are properly evaluated in addition to the diagnosis of origins and sources of hydrate bound gases molecular and isotopic data help to better identify hydrate bearing intervals and provide valuable insights into the dynamics of hydrate bearing sites

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

the influence of trench migration on slab penetration into the lower mantle abstract a two dimensional numerical convection model in cartesian geometry is used to study the influence of trench migration on the ability of subducted slabs to penetrate an endothermic phase boundary at 660 km depth the transient subduction history of an oceanic plate is modelled by imposing plate and trench motion at the surface the viscosity depends on temperature and depth a variety of styles of slab behaviour is found depending predominantly on the trench velocity when trench retreat is faster than 2 4 cm a the descending slab flattens above the phase boundary at slower rates it penetrates straight into the lower mantle although flattening in the transition zone may occur later leading to a complex slab morphology the slab can buckle independent of whether it penetrates or not especially when there is a localised increase in viscosity at the phase boundary flattened slabs are only temporarily arrested in the transition zone and sink ultimately into the lower mantle the results offer a framework for understanding the variety in slab geometry revealed by seismic tomography

the kola alkaline province of the cis and finland precise rbsr ages define 380 360 ma age range for all magmatism abstract twenty one new rb sr mineral and biotite whole rock isochrons are presented for intrusive rock complexes of the palaeozoic kola alkaline province in the cis and finland samples from the ultrabasic alkaline and nepheline syenite complexes khibina lovozero afrikanda turiy mys ozernaya varaka sokli and iivaara all yield ages in the very narrow range from 361 7 to 377 3 ma i e late middle devonian to upper devonian these data represent crystallization rather than intrusion ages but are nevertheless argued as marking very closely the time of intrusion there is no indication of bimodal age distribution the emplacement of the alkaline magmas is related temporally to the subsidence of the kontozero graben of the central kola peninsula sr initial ratios of the agpaitic nepheline syenites of the giant khibina and lovozero intrusions scatter between 0 703 and 0 704 because of high rb sr ratios of the rocks these isotope ratios exclude long term pre enrichment of alkalis in the sources of the magmas prior to segregation and intrusion and thus episodic enrichment during earlier rifting events in the northern baltic shield the magma formation may be related to wrench movements in the surroundings of the major trollfjord komagelv fault which forms part of the late caledonian arctic north atlantic megashear system

trench motion slab geometry and viscous stresses in subduction systems summary a semi analytic 3 d model for subduction within a newtonian viscous upper mantle provides a dynamically consistent means of computing viscous stress trench motion and slab geometry in subduction systems although negative slab buoyancy provides the basic driving force for subduction slabs that extend from the surface to the base of the upper mantle are oversupported by viscous stresses in the shallow 100 150 km mantle and undersupported by viscous stresses at greater depth in the upper mantle these deeper parts of the subduction system act as an engine for subduction while shallower parts act as a brake on trench motion trench migration rates and slab geometry reflect a competition between these two effects during steady state subduction trench migration rates vary approximately linearly with slab buoyancy and model rates of trench motion are in good agreement with the range of observed rates for a two layer upper mantle viscosity of 2 10 20 pa s above 300 km and 5 10 20 pa s below steady state slab dip increases as slab density decreases especially for very low density slabs which dip significantly more steeply than high density slabs the horizontal velocity at the top of the lower mantle measured relative to the foreland has a very large effect on trench migration rates rivalling or even exceeding that of slab buoyancy slab width parallel to the trench also has a significant effect on trench migration rates due to the viscous pressure of toroidal flow around the slab the stiffness of the subducting lithosphere does not exert a significant effect on trench migration rates or slab geometries for rigidities compatible with oceanic lithosphere very stiff slabs with elastic plate thicknesses more that 40 km or viscosities in the range of 10 25 10 26 pa s subduct significantly more slowly than weak slabs with trench migration rates in the range of half to a third that of weak slabs large unexpected effects on trench migration rates and slab geometry are exerted by the structure and density of the frontal prism and overriding plate indicating that local geology can exert important constraints on slab dynamics during non steady state subduction rates of trench migration respond rapidly as variably buoyant lithosphere penetrates into the asthenosphere in the absence of other driving forces for convergence trench migration rates can change by a factor of two or more in as little as 2 3 myr for example from 35 to 70 mm yr 1 when an oceanic piece of slab follows a continental one into the subduction system subduction of variable buoyancy lithosphere is accompanied by changes in slab dip with depth and through time

stress change and effective friction coefficient along the sumatra andaman sagaing fault system after the 26 december 2004 m w 9 2 and the 28 march 2005 m w 8 7 earthquakes the 2004 aceh and 2005 nias events are the two greatest earthquakes of the past 40 years with a total rupture of 1700 km long and a coseismic slip reaching up to 25 m these two earthquakes have caused large stress perturbations which significantly altered seismic activity in the sumatra andaman region using both detailed mapping of failure planes and various slip distributions we calculate this stress change along the sumatra andaman sagaing fault system from central sumatra to southern myanmar the static coulomb stress change delta cff and the observed seismic activity are in very good agreement with a coulomb index similar to 20 greater than the one obtained for random events compared to previous studies this high coulomb index confirms two important issues on the use of static stress change criterion unsuited to study near field aftershocks and only relevant for aftershocks analysis on large and mature faults at a time scale of several months the calculated delta cff distribution suggests that the 2004 and 2005 earthquakes inhibit failure on the north andaman rift and on the sagaing fault while failure is encouraged along the transform andaman zone the central andaman rift the west andaman fault the sumatra fault system and the offshore thrust faults west of sumatra island the maximum value of similar to 15 20 bar 1 5 2 mpa for delta cff is reached in the northern part of the sumatra fault system this high value together with the lack of major earthquake in the last 170 years result in a high seismic hazard for this region our results are also consistent with temporal evolution of both earthquakes location and focal mechanism prior to and after the events in particular we explain the occurrence and the mechanism of seismic swarms observed in the central andaman rift and along the west andaman fault finally our calculations reveal that the seismicity in the andaman rift zone can only be explained if mu 0 5 this result leads to two end member models one with a constant and high fault friction and one with spatial variations for which friction may depend on either the nature of the lithosphere oceanic versus continental or the fault type

detachment of part of the downgoing slab and uplift of the new hebrides vanuatu islands several seisniological observations suggest that there is a gap within the downgoing slab of australian lithosphere plunging beneath the new hebrides and ages of elevated coral tenaces on the new hebrides islands suggest that the islands arc rising rapidly we suggest that the creation of the gap in the slab within the last 1 m y occurred by part of the slab detaching from the rest and sinking rapidly into the underlying asthenosphere this lowered the downward force on the slab at shallower depths and therefore on the overriding plate and the continued decrease of this force as the deeper slab sinks freely has allowed the islands above it to rise although this suggestion cannot account for all aspects of the uplift of the islands it provides a simple mechanism for explaining the fairly young uplift and its distribution along the new hebrides arc

mixing behavior in quaternary garnet solid solution and an extended ellis and green garnet clinopyroxene geothermometer applications of the different formulations of garnet clinopyroxene geothermometers to crustal granulites and particularly to metamorphosed manganese formations show a wide scatter of estimates this is primarily related to the uncertainty associated with the formulation of the activity composition relationship in garnet solid solution besides that in the clinopyroxene solid solution the largest uncertainty appears to have resulted from the adopted mg mn interaction parameter in garnet which is yet to be experimentally determined a fresh statistical regression of eleven pairs of garnet and clinopyroxene equilibrated at identical p t conditions from a manganese formation in india yields wmg mngt wmn mg value of 1600 500 cal mols an extension of the ellis and green 1979 geothermometer known to work well in normal granulites incorporating terms associated with the activity coefficients of the exchangeable cations now yields consistent temperature estimation from normal granulites as well as from metamorphosed manganese formations the geothermometer however will not work well for garnet having xfe 0 65 and clinopyroxenes containing appreciable na al and or fe3

interplay between subduction and continental convergence a three dimensional dynamic model for the central mediterranean summary the mediterranean region is attracting considerable attention due to the complexities of its tectonic setting which is considered worldwide a unique natural laboratory for studying the occurrence of extensional tectonics in a general context of continental convergence the tyrrhenian apennine system is controlled by the west dipping subduction of the adria ionian lithosphere and by the near north south convergence between the african and eurasian plates we provide the first 3 d dynamic model of the central mediterranean that quantifies the effects of subduction and convergence on surface deformation in simplified geometry the axis of the model extends from sicily to the alps along the subduction hinge line a convergence rate of 1 cm yr 1parallel to the subduction hinge has been applied to the tyrrhenian block in agreement with global plate motion models density contrasts within the slab cause the gravitational sinking and roll back of the slab in the southern tyrrhenian domain modelling results show a gradual decrease of hinge retreat from south to north with values ranging between 8 and 2 mm yr 1 indicating that the arculate geometry of the hinge line along the italian peninsula is ultimately controlled by the interplay between subduction and convergence the pattern of vertical velocity along directions perpendicular to the hinge with subsidence in the foredeeps and uplift at the eastern border of the tyrrhenian domain is maintained along the whole italian peninsula with higher values in the southern areas

some remarks on the origin of foredeeps abstract geological observations show that there are two main types of foredeep basins on the base of their subsidence rate and their relationship with the associated orogen or accretionary wedge high rates of subsidence up to 1600 m myr and a ratio lower than 1 between the area of the elevated belt and the area of the basin characterize the foredeeps associated with w dipping subduction low rates of subsidence up to 300 m myr and a ratio higher than 1 between the area of the orogen and the area of the basin characterize the foredeeps associated with e or ne dipping subductions this observation enables us to interpret the slow filling of foredeeps of the first type and the much faster filling of foredeeps of the second type moreover w dipping subductions have a steeper and faster monocline bounding the base of the deep foredeep in contrast to e or ne dipping subductions which have a shallower monocline this differentiation supports the notion of the eastward mantle flow with respect to the lithosphere detected in the hot spot reference frame in this view the foredeep depth in w dipping subductions is mainly controlled by the roll back of the subduction hinge pushed by the relative eastward mantle flow while foredeep depth in e or ne dipping subductions is instead mainly generated by the load of the thrust sheets and by the roll back of the subduction hinge due to the advancing upper plate contrasting the upward push of the mantle the shape of the foredeeps is regularly arcuate in case of w dipping subduction while is linear or following the shape of the inherited continental margin in case of e or ne dipping subduction fold development in the two opposite foredeeps is significantly different in w dipping subduction the folds are transported down in subduction while they are forming and consequently they are poorly eroded syntectonic sediments drape a preserved fold in the e ne dipping subduction folds and thrust sheets are instead uplifted and deeply eroded this is also predicted by the envelope of the crests in the two distinct accretionary wedges in the w dipping case the envelope may dip toward the subduction while it is rising toward the hinterland in the other cases this approach may explain the differences in terms of geometries amount and rates of subsidence between the foredeeps around the adriatic plate for instance the foredeep of the southern alps back thrust belt foredeep is deformed in the western part by the later apenninic foredeep development associated with the sinistral oblique ramp of the w dipping adriatic plate subduction the relationship in the frontal parts of a thrust belt between accretion and eustacy has to be analyzed in terms of when where and how fast thrusts propagation took place with respect to the randomly distributed sea level third order fluctuations and by the type of subduction which is connected to the foredeep and thrust belt or accretionary wedge therefore the stratigraphic patterns in foredeep basins are peculiar and casual successions that record the combination of regional tectonic evolution with independent wider scale sea level changes

agpaitic magmatism in the northeastern baltic shield a study of the niva intrusion kola peninsula russia abstract the recently discovered niva agpaitic syenite intrusion is located 50 km southwest of the khibina complex central kola peninsula it forms a lense like body 1 5 2 km in size that outcrops over an area of 300 200 m on the southern bank of the imandra lake rb sr 379 21 ma and sm nd 383 58 ma isochron ages of the intrusion coincide with the age of the khibina and lovozero complexes a characteristic feature of the niva rocks is their extremely high agpaicity and enrichment in sr ba y zr nb rb lree and volatiles the mineral association in the agpaitic syenite is k feldspar aegirine diopside biotite lamprophyllite aenigmatite and shcherbakovite isotopic compositions indicate the absence of contamination by continental crust during formation of the niva alkaline magmas a depleted asthenospheric mantle reservoir mixed with a slightly enriched mantle component is suggested to be the main source of the niva magma a metasomatically altered garnet granulite xenolith from the western contact of the niva intrusion contains a primary paragenesis of diopside garnet and plagioclase with chemical characteristics similar to those in the granulite suite of the adjacent lapland zone and lower crustal granulite xenoliths from devonian pipes in the kandalaksha area additionally the granulite contains vinogradovite neptunite pectolite and natrolite which partially replace both primary minerals of the granulite assemblage and secondary minerals which were formed during a later regional metamorphism pt estimates support the lower crustal origin of the xenolith the u pb zircon age of the niva garnet granulite is 2711 7 ma

a modelling study of vertical surface displacements at convergent plate margins summary during the evolution of a subduction zone system variations are likely to occur in for example surface plate velocities and buoyancy of the subducting lithosphere we quantify vertical surface displacements at convergent plate margins resulting from such imposed variations for this purpose we use a 2 d numerical model in which the lithospheric plates have an effective elastic thickness we first define a model in which the subducting plate is driven by its negative buoyancy and a velocity at its surface side boundary its equilibrium topography after around 2 myr is the reference level for examination of surface displacements resulting from variations in buoyancy velocity of the surface plates friction along the interplate contact and subduction zone roll back we find that a decrease increase in buoyancy of the subducting material leads to a deepening uplift of the plate margins an increase in friction along the subduction fault deepens the overriding plate margin subduction zone roll back due to sinking of the negatively buoyant subducting plate induces subsidence of the overriding plate margin this subsidence is reduced when roll back takes place in a land locked basin setting trench retreat forced by the motion of the overriding plate is characterized by higher topography of the overriding plate margin than the case of retreat due to the sinking of the negatively buoyant slab in the first case in plane stress in the back arc region is compressive while it is tensional for roll back due to the sinking of the slab we conclude that vertical surface displacements during ongoing subduction may reach a magnitude of a few kilometres on the overriding and subducting plate margins