dynamic modelling of stress accumulation in central italy we develop a viscoelastic dynamic model to calculate the stress distribution in central italy induced by different geodynamic processes likely active in the region africa eurasia convergence adriatic underthrusting subduction and asthenospheric upwelling underneath tuscany model results are compared with the hypocentral distribution of earthquakes and with the orientation of p and t axes available for the region the first order characteristics of the seismotectonic observations in central italy are subcrustal seismicity the pair extension compression in tuscany and in the outer apennines and the flexural behaviour evident in long wavelength orientation of t axes these features are reproduced by the combination of two processes underthrusting of the adriatic plate under the apennines and asthenospheric upwelling underneath tuscany the comparison with available data does not give insights about slab pull as an ongoing process

attenuation of high frequency seismic waves beneath the central andean plateau observed patterns of high frequency seismic wave attenuation suggest that near 22 s the upper mantle structure beneath the central andean plateau changes along strike to the south contrasting regions of high and low upper mantle seismic wave attenuation beneath the plateau are identified based on striking azimuthal variations in the character and frequency content of p and s waves propagating beneath the plateau to a portable seismic network deployed in jujuy province argentina 24 s 65 w ray paths from intermediate depth earthquakes located north and northwest of the network transmit seismic waves with a higher frequency content than ray paths from earthquakes at similar depths and distances but located west and south of the network the estimated apparent q values fall into two categories qp 500 and qs 350 for the high q paths and qp 350 and qs 200 for the low q paths in addition sn phases from regional crustal earthquakes in the subandean foreland fold thrust belt to the north propagate efficiently to the jujuy network while sn is not observed from foreland earthquakes located at similar distances to the south of the network these observations combined with data from previously reported wave propagation studies suggest that south of about 22 s the upper mantle beneath the plateau and its adjacent foreland thrust belt is more highly attenuating than the upper mantle farther north forward modeling of the q measurements made at jujuy indicates that the observations can be explained either by a variable thickness high q upper plate beneath the plateau or by a thin variable width very low q zone in the asthenospheric wedge above the subducted slab this lateral variation in upper mantle structure coincides with two physiographically distinct segments of the central andean plateau and its adjacent foreland thrust belt to the east the bolivian altiplano and subandean ranges in the north and the argentine puna and santa barbara system in the south we interpret the north south change in upper mantle attenuation and the corresponding changes in physiography topography and tectonic style at the surface to be due to a mantle lid that is thicker beneath the altiplano and the subandean belt than beneath the puna and the santa barbara ranges

segmentation and configuration of subducted lithosphere in italy an important control on thrust belt and foredeep basin evolution the apennine foredeep basin system located on continental crust consists of four roughly parallel basin segments each containing a flexural outer rise and an inner trough that deepens continuously toward the thrust belt flexural modeling of subsidence data across the basin indicates that this foredeep basin results mainly from deep subsurface loads acting on the subducted slab from within the subduction zone rather than from the topographic load of the apennine mountains the geometry of the four basin segments observed within the apennine foredeep is thus inferred to reflect subduction of segmented lithosphere where each lithospheric segment corresponds to a segment of the foredeep basin beneath the basin and in the foreland the lithosphere is continuous segmentation occurs only at depth beneath the thrust belt this pattern olf segmentation of the subducted lithosphere at depth is reflected in the large scale geometry of the outer part of the apennine thrust belt and appears to control it similarly the pattern of segmentation also appears to control the position of the eastern limit of back arc type extension in western italy these relations strongly suggest a genetic relation between subduction zone processes foredeep basin geometry and thrust belt evolution in which the configuration of the subducted lithosphere largely controls the large scale surface deformation observed in the apennine system

active intraplate deformation within adria examples from the adriatic region abstract recent seismicity recorded in the adriatic region suggests that active deformation occurs within the adria plate earthquakes are common in the adjoining chains and in particular along the northern apennines in the southern alps and along the dinaric albanian hellenic thrust fronts in addition medium to low grade seismicity also affects the foreland and is clustered along two main transects through the centre of the adriatic region the two zones of active deformation correspond to the mid adriatic ridge and the gargano region and approximately trend nw se and e w respectively several interpretations have been proposed to explain the deformation observed within the adriatic foreland including a variety of contrasting kinematic and geodynamic models we attempted to unravel the complicated structural evolution of the central adriatic region by interpreting industrial and ministerial shallow and deep seismic profiles and the recent seismicity recorded in the adriatic the subsurface study allowed us to better understand the structural style timing and rates of deformation affecting the mid adriatic ridge and the gargano region and to unravel the relationships of active deformation with respect to the adjacent chains the results of our study indicate that the structural setting and type of seismicity of the central adriatic can be explained in terms of active intraplate deformation with the reactivation of inherited discontinuities

p and s wave attenuation in the region of marche italy we have found empirical functions to describe the amplitude decay with distance of p and s waves in the marche region of italy these functions result from the analysis of spectral amplitudes of 161 three component digital records from 26 earthquakes recorded by the marchesan seismograph network rsm jointly run by the italian seismological survey of rome ssn and the geophysical observatory of macerata ogsm the events have magnitudes in the range 1 8 m 3 6 and were recorded at hypocentral distances between 9 and 125 km in order to study the possible variations of the quality factor q with the direction of wave propagation q anisotropy we made independent estimates of the attenuation for p sv and sh waves using the longitudinal vertically polarized and transversely polarized components respectively for the whole distance range and in the frequency band analyzed 1 f 25 hz the frequency dependence of q can be approximated as qp 108 4 f0 70 for the p waves qsv 105 6 f0 65 for sv waves and qsh 86 1 f0 70 for sh waves the small difference between the sh and sv spectral decay is within the statistical uncertainties of the estimates suggesting that q anisotropy is negligible however when only deeper wave paths are used in the analysis we observed that at low frequencies f 4 hz the values of qsv are significantly above the estimates of qsh this observation provides additional evidence of the seismic anisotropy reported by margheriti et al geophys res lett 23 2721 2724 1996 for the region the low q anisotropy is possibly related with the weakly deformed adriatic foreland the high attenuation observed for both p and s waves can be explained by the double low velocity layer reported by ponziani et al tectonophysics 252 391 418 1995 since the seismic energy can be trapped between these layers in addition the low values of q obtained for distance ranges between 5 and 50 km can be due to the presence of hot material in the lower crust beneath the northern apennines as suggested by mele et al j geophys res 102 11863 11875 1997

the tyrrhenian apennines system structural setting and seismotectogenesis abstract the available geological geophysical and seismological data of the tyrrhenian apennines system are briefly reviewed and synthesized in a number of sketch maps integrating and comparing these data a large scale simple shear kinematic model is proposed for the tyrrhenian apennines system that is here interpreted as an oblique asymmetric passive continental margin developed during the progressive counterclockwise rotation of the adriatic foreland from the late miocene to the present in the frame of this working hypothesis the thinning and stretching of the tyrrhenian crust were accommodated by progressive eastward motion of extensional allochthonous slices bounded by low angle west dipping oblique shear zones rotating above a common detachment that penetrates into the mantle lithosphere at any given time the most external eastward moving slice is bounded by a transtensional shear zone at the top and by a transpressional one at the bottom the displacement within the hanging wall of the transtensional shear zone is accommodated by high angle normal faults and en echelon grabens the displacement within the hanging wall of the transpressional shear zone is instead accommodated by co axial en echelon sets of folds thrusts and high angle strike slip faults thus the development of the late miocene pliocene pleistocene contractional strain field of the apennines of italy might not be a consequence of complex subduction processes but subordinate to the deformation processes that lead to the progressive opening of the tyrrhenian sea the seismicity of the crustal apennines of italy is analyzed in the frame of the above history of progressive deformation and a seismotectogenetic model is outlined

deep crustal roots beneath the northern apennines inferred from teleseismic receiver functions abstract teleseismic waveforms recorded by a regional array crossing the northern italian peninsula and northern corsica are analyzed using the receiver function technique to determine the first order crustal structure the receiver function approach is used to isolate receiver side ps conversions generated at the crust mantle boundary and any major velocity discontinuity beneath the stations we used the time delay between the direct p wave and the ps wave converted at the moho discontinuity to infer crustal thickness beneath the stations the crust mantle boundary is estimated at 25 1 km of depth in northern corsica 20 2 km beneath the elba island and 20 to 24 2 km beneath tuscany in the eastern portion of the array moho depth increases from 28 2 km beneath the adriatic coast to 49 3 km beneath the apennine chain in a distance of about 100 km a double ps conversion produced beneath the val tiberina graben to the west of the apennines corresponds to interfaces as deep as 20 2 and 52 2 km this observation supports the hypothesis of partial overlapping between the shallow tuscan moho and the deeper adriatic moho the westward deepening of the adriatic moho beneath the northern apennines can be explained by lithosphere delamination that has dragged downward the adriatic lower crust the deep crustal root estimated beneath the northern apennines indicates that this portion of the chain is isostatically overcompensated

compressional velocity structure and anisotropy in the uppermost mantle beneath italy and surrounding regions travel times of about 39 000 pn arrivals recorded from regional earthquakes by the italian telemetered seismic network and by stations of nearby countries are inverted to image lateral variations of seismic velocity and anisotropy at subcrustal depth in italy and surrounding regions this method allows simultaneous imaging of variations of pn velocity and anisotropy as well as crustal thickness variations the po plain the adriatic sea and the ionian sea have normal to high pn velocities in contrast lower velocities 7 9 8 0 kin s are imaged in italy beneath the western alps the northern apennines and eastern sicily and nearby calabria as well as in northern albania and beneath the pannonian basin low pn velocities beneath the northern apennines correlate with present day extension and may have resulted from thermal anomalies in the uppermost mantle possibly due to delamination processes low velocities are consistent with the high attenuation zone inferred in the uppermost mantle beneath the internal apennine units and the tyrrhenian margin of the peninsula by mele et al 1996 1997 on the contrary low velocities beneath the western alps may be an apparent effect due to the abrupt thickening of the crustal roots pn velocity is anisotropic in the study area with a maximum amplitude of 0 2 km s the largest anisotropic velocity anomalies are observed along the major arc structures of italy i e the northern apennines and the calabrian arc indicating that these features are controlled by uppermost mantle processes the anisotropy anomaly along the calabrian arc extends as far as albania but ends abruptly north of this area suggesting that a lithospheric discontinuity is present along the northern albanian border

seismotectonics of the calabrian arc abstract relocation of intermediate and deep earthquakes of tyrrhenian sea area through joint hypocenter determination for the period 1962 1979 has allowed a more detailed definition of the geometry of this peculiar benioff zone earthquakes dip along a quasi vertical plane to 250 km depth there is a 50 dip in the 250 340 km depth range and a low dip angle to 480 km depth the structure sketched from the hypocenters is almost continuous but most energy has been released in the 230 340 km depth interval an evaluation of fault plane solutions of intermediate earthquakes in this area indicates predominance of down dip compressions in the central part of the slab at the border strike slip motion occurs independent of depth some earthquakes that occurred at intermediate depth less than 100 km along the ionian margin of calabria show predominance of reverse faulting with the p axis oriented se nw however shallow earthquakes in the calabria sicily region indicate a more complex motion with predominance of normal faulting a possible interpretation of these features according to the available geological history which involves subduction of continental lithosphere is discussed

lateral variations of high frequency seismic wave propagation at regional distances across the turkish and iranian plateaus in this study we investigate crustal and uppermost mantle physical properties that characterize some of the continental plateaus of the middle east this is done as part of a larger effort to map and compare high frequency wave propagation at regional distances across the earth s continental plateaus thousands of short period wwssn seismograms recorded at stations located in the middle east and produced by earthquakes with epicentral distances less than about 20 were examined visually in an effort to study lateral variations of high frequency 0 5 2 hz seismic wave propagation across this area particularly to the north of the zone of continental collision between the africa arabian and eurasian plates variations of frequencies and amplitudes of sn and lg relative to p are examined and mapped throughout the region and this work is supplemented by a study of velocities of pn sn and lg sn amplitude variations are very striking in this area an important observation of this study is that sn propagates efficiently beneath a major part of the turkish and iranian plateaus sn is strongly attenuated however in the northernmost portion of the plateaus south of the black and caspian seas and in an area between the two seas these regions are characterized in general by active tectonism including volcanism faulting and folding however this active tectonism is not restricted to the areas of high sn attenuation but appears to extend beneath other parts of the iranian and turkish plateaus patterns of lateral variations in the propagation of lg are not as consistent as those for sn lg propagates efficiently across turkey iran and adjacent regions but the lg waves that cross the turkish and iranian plateaus are weak and have relatively long predominant periods of about 2 5 s the lg phase is not observed when the path of propagation crosses the southern caspian and the black seas consistent with the evidence of oceanic type crustal structure beneath these seas lg is also not observed from subcrustal events located in the hindu kush region the velocity of pn beneath most of the plateau areas and particularly in the regions of sn attenuation ranges between 8 0 and 8 2 km s velocities for sn and lg throughout the middle east are about 4 5 and 3 4 km s respectively it is possible to interpret the efficient sn propagation and the relatively normal uppermost mantle p and s velocities beneath a major part of the turkish iranian plateau as indicating a partial underthrusting of the arabian continental plate beneath the iranian and turkish continental blocks alternatively the mapped regions of efficient sn propagation and high sn attenuation may result from a differential cooling of a past thermal anomaly in this case the isostatic compensation of the plateaus could be due to an anomalously hot uppermost mantle that developed behind the mid to upper tertiary subduction zone s in this region and to a subsequent crustal shortening after the collision of the arabian and eurasian plates our results cannot uniquely confirm the validity of the above proposed models

2d coda and direct wave attenuation tomography in northern italy abstract a 1d coda method was proposed by mayeda et al 2003 in order to obtain stable seismic source moment rate spectra using narrowband coda envelope measurements that study took advantage of the averaging nature of coda waves to derive stable amplitude measurements taking into account all propagation site and s to coda transfer function effects recently this methodology was applied to microearthquake data sets from three subregions of northern italy i e western alps northern apennines and eastern alps because the study regions were small ranging between local to near regional distances the simple 1d path assumptions used in the coda method worked very well the lateral complexity of this region would suggest however that a 2d path correction might provide even better results if the data sets were combined especially when paths traverse larger distances and complicated regions the structural heterogeneity of northern italy makes the region ideal to test the extent to which coda variance can be reduced further by using a 2d q tomography technique the approach we use has been developed by phillips et al 2005 and is an extension of previous amplitude ratio techniques to remove source effects from the inversion the method requires some assumptions such as isotropic source radiation which is generally true for coda waves our results are compared against direct s wave inversions for 1 q and results from both share very similar attenuation features that coincide with known geologic structures we compare our results with those derived from direct waves as well as some recent results from northern california obtained by mayeda et al 2005 that tested the same tomographic methodology applied in this study to invert for 1 q we find that 2d coda path corrections for this region significantly improve upon the 1d corrections in contrast to california where only a marginal improvement was observed we attribute this difference to stronger lateral variations in q for northern italy relative to california

geochemical and petrological evidence of the subduction of delaminated adriatic continental lithosphere in the genesis of the neogene quaternary magmatism of central italy abstract serri g innocenti f and manetti p 1993 geochemical and petrological evidence of the subduction of delaminated adriatic continental lithosphere in the genesis of the neogene quaternary magmatism of central italy in m j r wortel u hansen and r sabadini editors relationships between mantle processes and geological processes at or near the earth s surface tectonophysics 223 117 147 the neogene quaternary magmatism of the northern apenninic arc took place in four phases separated in space and time which become progressively younger from west to east phase i 14 ma phase ii 7 3 6 0 ma phase iii 5 1 2 2 ma phase iv 1 3 0 1 ma this magmatism is the result of the activation of three physically separate sources 1 the adriatic continental crust extracted from the mantle in the late proterozoic 2 a strongly refractory recently k enriched harzburgitic mantle located in the mechanical boundary layer mbl of the lithosphere and 3 a recently metasomatized cpx rich mantle compositionally variable from iherzolite to wehrlite clinopyroxenite interpreted as an ephemerally k enriched asthenosphere the adriatic continental crust is the dominant source of the acid plutonic and volcanic rocks of the tuscan region the acid magmatism is mostly found inside an ellipsoidal area about 150 300 km centred on giglio island here defined as the tuscan crustal dome within this area mantle derived magmas unaffected by important crustal contamination processes and mixing with crustal anatectic melts have so far not been found pure crustal magmas are rare but are represented for example by some of the san vincenzo and roccastrada rhyolites virtually all the tuscan acid centres show evidence of mixing with potassic mantle derived magmas major and trace elements as well as 87 sr 86 sr and 143 nd 144 nd data on primitive rocks mg 65 reveal two groups of mantle derived magmas these define two distinct mantle enrichment trends both essentially due to the additions of k rich components which metasomatized separate compositionally diverse upper mantle sectors in both cases the most remarkable mineralogical effect of these enrichment processes is the production of variable amount of phlogopite through reaction between fluids and or melts with the mantle the rocks of group i ol hy and q normative lamproites ultrapotassic high mg latites ultrapotassic shoshonites and shoshonites saturated trend are considered to be derived by partial melting at low pressure 87 sr 86 sr 0 717 ce sr 0 3 and k 2 o na 2 o 6 7 and low ratios of 143 nd 144 nd 0 5121 0 5120 and ba la 87 sr 86 sr 0 712 and k 2 o na 2 o 8 9 values and low 143 nd 144 nd ba la 0 10 ratios these constraints do not allow to exclude a complete identity between the k rich components which metasomatized the mantle sources of the saturated and undersaturated trend magmas the geochemical and isotopic features of the components that metasomatized the mantle sources of the northern apenninic arc magmatism can be explained by a geodynamic process which causes a large amount of crustal materials to be incorporated within the upper mantle we propose that the delamination and subduction of the adriatic continental lithosphere related to the still ongoing northern apennine continental collision provide a viable mechanism to explain the genesis and eastward discontinuous migration of the magmatism in central italy the subduction of delaminated lithospheric mantle with lower crustal slivers would have exposed uppermost mantle adriatic mbl and crustal units previously imbricated in the apennine chain to the heating advected by the upwelling of a recently and ephemerally k enriched asthenospheric mantle wedge and by the underplating of magmas derived from it we consider that the diapiric uprising of a hot crustally contaminated asthenosphere occurs in the wake left above the sinking of the adriatic delaminated subducting continental lithosphere the delamination subduction process of the adriatic lithosphere has probably started in the early middle miocene but earlier than 15 14 ma ago as indicated by the age and petrologic characteristics of the first magmatic episode sisco lamproite of the northern apennine orogenesis

the origin of the apenninic arcs abstract the apennine area is composed of a variety of arcs some of them forming the mountain chain and some of them being submerged below the po plain and the adriatic sea the structure of the apennine mountain belt is that of a succession of northeast east and southeastward directed arcs corresponding to single thrust sheets each arc is limited by n s oriented dextral transcurrent faults data from volcanology gravimetry and reflection seismic surveys show that an important section of the upper mantle is involved in these features the geodynamic reconstruction during neogene times i e the period of apenninic up building suggests that a unique arc existed from oligocene to middle miocene times as a result of the rotation of the sardo corso continental block this arc was successively disrupted in several segments which advanced and shortened differentially as a function of the resistance to the orogenic expansion opposed by the rigid adriatic foreland this new tectonic phase which is responsible for most of the apenninic features has been triggered by the differential evolution of anomalous mantle domes in the tyrrhenian area from post middle miocene times onwards

vertical gps ground motion rates in the euro mediterranean region new evidence of velocity gradients at different spatial scales along the nubia eurasia plate boundary we use 2 5 to 14 years long position time series from 800 continuous global positioning system gps stations to study vertical deformation rates in the euro mediterranean region we estimate and remove common mode errors in position time series using a principal component analysis obtaining a significant gain in the signal to noise ratio of the displacements data following the results of a maximum likelihood estimation analysis which gives a mean spectral index 0 7 we adopt a power law white noise stochastic model in estimating the final vertical rates and find 95 of the velocities within 2 mm yr with uncertainties from filtered time series 40 smaller than from the unfiltered ones we highlight the presence of statistically significant velocity gradients where the stations density is higher we find undulations of the vertical velocity field at different spatial scales both in tectonically active regions like eastern alps apennines and eastern mediterranean and in regions characterized by a low or negligible tectonic activity like central iberia and western alps a correlation between smooth vertical velocities and topographic features is apparent in many sectors of the study area glacial isostatic adjustment and weathering processes do not completely explain the measured rates and a combination of active tectonics and deep seated geodynamic processes must be invoked excluding areas where localized processes are likely or where subduction processes may be active mantle dynamics is the most likely process but regional mantle modeling is required for a better understanding

plio quaternary vertical motion of the northern apennines insights from dynamic modeling we test the effects of different geodynamic mechanisms on the pliocene to present day dynamics and in particular on the vertical motions of the northern apennines system by means of two dimensional finite element models we show that the pliocene features of the northern apennines exhumation of deep rocks in western italy and fast subsidence rates in the adriatic foredeep can be adequately reproduced by a model slab characterized by the passive sinking of a negatively buoyant adriatic slab and by the upwelling of buoyant asthenosphere beneath western italy model slab is however not able to account for the quaternary and present day pattern of vertical motion of eastern italy and of the adriatic basin a deep configuration of the system characterized by a detached slab model detach may explain the vertical motions of these areas constrained by uplift of the eastern portions of the apenninic chain by an eastward shift of the foredeep depocenter and by lower subsidence rates in the adriatic basin the tectonic and stratigraphic data showing major differences between the tertiary and the quaternary evolution of the northern apennines are discussed

uppermost mantle velocities and anisotropy beneath europe pn data collected within southern europe and the mediterranean are used to tomographically image variations in both seismic velocity and seismic anisotropy seismic anisotropy is an essential part of the inversion and without it several low velocity features within the uppermost mantle could not be properly imaged the technically active mantle of southern europe has much lower seismic velocities 7 6 8 1 km s than the more stable mantle of the sub african plate of the adriatic sea 8 3 km s however the most dramatic features within europe s uppermost mantle relate to the apennine dinaride and hellenide arcs these arcs all have extremely low 5 amounts of arc parallel anisotropy the tyrrhenian and aegean back arc regions also show low velocities 7 7 7 9 km s but less anisotropy the same may be true for the pannonian basin but the tomography has poorer resolution there a model explaining these observations focuses on subducted water metasomatizing the mantle wedge the addition of water causes melting creates arc volcanism lowers the seismic velocity and enhances the formation of anisotropy due to preferential olivine orientation for collisional arcs of the northern mediterranean arc parallel anisotropy has formed in response to compression across them and extension along them within back arc regions water is no longer a major factor instead convection associated with subduction and back arc extension controls the anisotropy

shear wave attenuation tomography of the lithosphere asthenosphere system beneath the mediterranean region abstract 3 d tomographic images of the shear wave attenuation structure of the lithosphere asthenosphere system beneath the mediterranean region are derived from the inversion of anelastic attenuation coefficients obtained for the fundamental mode of nearly 3100 rayleigh wavetrains recorded at very broad stations placed in the mediterranean region the two station method is used to compute 124 interstation path averaged attenuation coefficients within the 15 120 s period range an adaptation to attenuation coefficients of the yanovskaya ditmar formulation specifically devoted to path averaged surface wave dispersion data and body waves permits to obtain the local attenuation coefficients and their uncertainties for all the period range at every node of a grid covering the studied region maps of stretching azimuths and mean radii of the resolution ellipses are also obtained for every period all of these parameters are strongly related to the spatial path coverage data quality and data coherency a local attenuation dispersion curve is defined then in each grid node and the stochastic inversion of these curves yields 1 d models of inverse shear wave intrinsic quality factor q 1 the shear wave velocity models obtained for every grid node in an elastic tomography developed in a previous work for the same domain are taken as reference models for the inversion process the kriging contouring algorithm is used to represent the resulting 3 d q 1 structure by means of horizontal and vertical cross sections down to 160 km depth the main patterns of this anelastic structure are mostly in agreement with the structural and tectonic features obtained in the reference shear velocity tomography the local attenuation maps manifest notable lateral changes for periods of 20 and 40 s whereas more moderate spatial variation is observed for longer periods the 3 d anelastic models show that at shallow levels the attenuation in the western mediterranean is usually slightly higher than in the central eastern part the area from corsica sardinia islands to italy presents high q 1 values of 40 0 70 0 10 3 down to 75 km depth whereas q 1 does not exceed 20 0 10 3 in the rest of the mediterranean basin except for the south eastern part where slightly higher values are obtained 20 0 30 0 10 3 those high q 1 values suggest the presence of asthenospheric material at shallow levels which would be related to distensive tectonic processes responsible for the recent opening of extensional basins as the tyrrhenian sea and the ligurian provencal basin in the eastern mediterranean the vertical cross sections at different latitudes and longitudes reveal at deep levels some structural patterns consistent with signs of subduction which are linked to the convergence between the african and eurasian plates

a local magnitude scale for crustal earthquakes in italy abstract in this study a new local magnitude m l scale is derived for the italian region from the analysis of seismic signals recorded by a dense broadband network between 2003 and 2009 the m l computation is performed by measuring peak amplitudes of synthetic wood anderson seismograms based on the original definition by richter 1935 these amplitudes are inverted for earthquake magnitudes site corrections and the distance dependent term that corrects amplitudes for geometrical spreading and scattering anelastic attenuation the latter is expressed as a function of r hypocentral distance and log r by means of a two parameter linear form and is constrained at the distance of 100 km in accordance with richter s definition using the newly derived distance dependent term and site corrections magnitudes are consistently estimated for the stations recording each earthquake within an overall root mean square residual of 0 18 no anomalous trend with distance is appreciable in the magnitude residuals for the italian region as a whole or for smaller areas along the italian peninsula the only exception is represented by northeastern italy where increasingly negative residuals are observed for distances decreasing below 100 km the spatial pattern of site corrections reveals contrasting attenuation properties of the crust between the tyrrhenian side of the italian peninsula high attenuation west of the apennines and the adriatic side low attenuation east of the apennines as previously found in other studies finally the magnitudes computed with the new m l scale are compared with global estimates of body wave magnitude m b and with local and regional estimates of moment magnitude m w they display a satisfactory agreement in most of the magnitude range considered roughly from 3 to 5 5 online material tables of station parameters maximum amplitudes and magnitude

interactions between mantle upwelling drainage evolution and active normal faulting an example from the central apennines italy summary in this paper we show that the processes that have shaped the quaternary surface development of the apennines in central italy are all consequences of a single subcrustal process the upwelling of the mantle the relationship between gravity and topography shows that mantle convection is responsible for a long wavelength 150 200 km topographic bulge over the central apennines and stratigraphic evidence suggests this bulge developed in the quaternary active normal faulting is localized at the crest of this bulge and produces internally draining fault bounded basins these basins have been progressively captured by the aggressive headward erosion of major streams that cut down to the sea on the flanks of the regional bulge the only surviving closed basins are those on the apennine watershed most distant from the marine base level where continued normal faulting is still able to provide local subsidence that defeats their capture by the regional drainage network understanding the competition between regional capture and local fault related subsidence of intermontane basins is crucial for recognizing potentially hazardous active faults in the landscape and also for interpreting the sediment supply to adjacent offshore regions central italy provides a good modern analogue for processes that are probably common in the geological record particularly on rifted margins and intracontinental rifts but may not have been fully appreciated

the influence of iron sulfide oxidation on the sulfur isotope analysis of miocene phosphorites from florida usa abstract comparison of 34 s values obtained using an oxic and anoxic gravimetric method on miocene phosphorite from florida indicates that oxidation of fe sulfides can significantly lower s isotope values the lowest s isotope values are obtained from phosphorite grains sampled from cores recovered from northeastern florida nef 34 s values of nef phosphorite grains range from 11 9 to 19 2 cdt when prepared under oxic conditions and 17 4 to 22 8 cdt when prepared under anoxic conditions phosphorites from nef have a mean terrigenous insoluble residue content of 29 wt and fe sulfide contents of 0 5 to 3 wt the chemical composition of the carbonate fluorapatite from nef phosphorite indicates minimal alteration by postdepositional weathering most of the fe sulfide oxidation is believed to occur rapidly during drying of bulk powdered samples after washing to remove soluble sulfate and possibly during filtration and acid digestion phosphorites sampled from shallow open pit mines of the central florida phosphate district cfpd have low terrigenous insoluble residue contents mean of 7 wt and pyrite contents less than 0 5 wt 34 s values of cfpd phosphorite grains range from 22 4 to 25 5 cdt and show no difference for samples prepared under oxic and anoxic conditions 34 s values of cfpd phosphorites are slightly enriched relative to miocene seawater and indicate formation immediately below the suboxic anoxic sediment boundary regional differences in pyrite content and s isotope composition of the phosphorite are related to proximity to terrigenous sediment sources and postdepositional weathering lower than miocene seawater 34 s values of nef phosphorites obtained using the anoxic gravimetric method may reflect microbial oxidation of sulfide at the suboxic anoxic sediment interface or extremely rapid oxidation of pyrite during sample filtration in the lab

orogens and slabs vs their direction of subduction abstract subduction zones appear primarily controlled by the polarity of their direction i e w directed or e to nne directed probably due to the westward drift of the lithosphere relative to the asthenosphere the decollement planes behave differently in the two end members in the w directed subduction zone the decollement of the plate to the east is warped and subducted whereas in the e to nne directed it is ramping upward at the surface there are w directed subduction zones that work also in absence of active convergence like the carpathians or the apennines w directed subduction zones have shorter life 30 40 ma than e or ne directed subduction zones even longer than 100 ma the different decollements in the two end members of subduction should control different ptt paths and therefore generate variable metamorphic assemblages in the associated accretionary wedges and orogens these asymmetries also determine different topographic and structural evolutions that are marked by low topography and a fast eastward migrating structural wave along w directed subduction zones whereas the topography and the structure are rapidly growing upward and expanding laterally along the opposite subduction zones the magmatic pair calc alkaline and alkaline tholeiitic volcanic products of the island arc and the back arc basin characterise the w directed subduction zones magmatic rocks associated with e or ne directed subduction zones have higher abundances of incompatible elements and mainly consist of calc alkaline shoshonitic suites with large volumes of batholithic intrusions and porphyry copper ore deposits the subduction zones surrounding the adriatic plate in the central mediterranean confirm the differences among subduction zones as primarily controlled by the geographic polarity of the main direction of the slab the western margin of the adriatic plate contemporaneously overridden and underthrust europe toward the west to generate respectively the alps and the apennines while the eastern margin subducted under the dinarides hellenides these belts confirm the characters of the end members of subduction zones as a function of their geographic polarity similarly to the pacific subduction zones

radiogenic helium degassing and rock fracturing a case study of the southern apennines active tectonic region gas from mud volcanoes dry mofettes springs and wells were sampled in a region of active tectonics and high seismicity in the southern apennines italy where there is a long history of disastrous earthquakes with the latest ms 6 9 occurring in 1980 the fluids consist of a mixture of mantle derived and crust derived volatiles with a low atmosphere derived contribution as identified by the he isotope signature and he ne ratio measurements one year of monthly monitoring of the he concentrations and he isotopes revealed no seasonal modifications or variations induced by low seismicity there are extraordinary high outputs of 4he produced in the crust in the area up to 2 5 1028 atoms yr 1 these outputs cannot be solely due to the whole rock production rate and a long lasting diffusion degassing through the crust of the produced 4he this study explored the relation between the volume of fractured rock and the related release of he the results support that crustal degassing can be controlled by tectonic events resulting in earthquakes the high seismicity in this sector of the apennines provides the conditions necessary for a massive release of he that has accumulated in the rock over a long time period we identified that the assessed high crustal 4he output can be attributed to an intense fracturing of a calculable volume of rock which gives new constraints on the volume of rock involved in high magnitude earthquakes in the region

the interplay between global tectonic processes and the seismic cycle in the umbria marche seismogenic region summary for the central apennines peninsular italy a series of tectonic mechanisms are reproduced by means of finite element models in order to study the effects of active tectonics on the seismic cycle in the umbria marche seismogenic zone continental extension and rift push effects induced by small scale convection are modelled within 2 d viscoelastic models of the crust lithosphere system in vertical cross sections perpendicular to the strike of the major tectonic structures under study namely the apennines and the colfiorito fault zone where the 1997 seismic sequence took place with the aim of constraining the active tectonic mechanisms at the regional scale and the behaviour of the fault in the seismogenic zone at the local scale modelled baseline rate of change are compared with newly acquired gps data retrieved from the two permanent gps receivers of camerino came and elba elba deliberately installed along the modelled transect these receivers are located at both edges of the continental extension in the front of the apennines close to the adriatic plate in the east and in the rear of the chain in the tyrrhenian domain the deformation pattern inferred from seismicity and from the geodetic data is consistent with small scale convection in the tyrrhenian domain which reproduces extension in the rear of the apennines and compression in the front of the chain a convective mechanism associated with backarc opening and doming of the asthenosphere provides an extensional rate along a baseline connecting two sites in the front of the chain camerino and in its rear elba comparable to the observed baseline rate of change the viscosity of the lower crust plays a fundamental role in determining the style of stress in the crust lithosphere system once constrained by means of the extensional baseline rate inferred from gps the modelled slip across the colfiorito fault and the modelled earthquake recurrence time are consistent with the 1997 normal fault event and with palaeoseismicity respectively

geothermal structure of the tyrrhenian sea over 200 heat flow measurements have been carried out in the tyrrhenian sea mostly since 1977 as part of the heat flow map for italy and surrounding seas italian program as well as part of the french italian tyrrhenia program all data have been corrected for the influence of the plio quaternary sedimentation measured heat flow values were evaluated according to the sedimentary environment following criteria outlined by sclater et al 1980 in this way areas affected mainly by shallow thermal effects in the zones with maximum density of data have been displayed filtering out these shallow contributions the major regional geothermal trends were obtained the tyrrhenian basin can be divided in two distinct regions with different geothermal characteristics the northern part of the basin appears to be connected with the continental structures involved in the interaction between the corsica and tuscan margins the maximum heat flow area located offshore civitavecchia is related to the tensile structures of the tuscan geothermal province the most important thermal features of the central and southern part are located along the central belt trending nw se the lowest heat flow values are located in the peri tyrrhenian areas and are related to the continental margins of the basin while the maxima are connected with the stretched and thinned area localized in the bathyal plain taking into account the heat flow data as well as other geological and geophysical information an evolutionary model for the tyrrhenian sea is presented suggesting a generalized stretching of an initially continental crust accompanied by intensive rifting phases which may have favoured the onset of oceanic type crust

the deep seismicity of the tyrrhenian sea the deep seismicity of the tyrrhenian sea is analysed using data from a new instrumental catalogue of the seismicity of the italian area we use algorithms for the determination of absolute and relative hypocentral locations and for the evaluation of the geometry and coherence of the state of stress within the subducting slab the structure of the benioff zone although simpler than previously indicated reveals anomalous traits both in the seismicity distribution and in the stress geometry confirming that standard subduction models cannot be applied in the tyrrhenian region the velocity anomaly and the location of few isolated events indicate that the subducted slabs extend to the north along the apenninic chain approximately to the latitude of the irpinia region in central italy 42 n

linking the alps and apennines subduction systems new constraints revealed by high resolution teleseismic tomography abstract we report a new model of the upper mantle structure beneath italy obtained by means of p wave teleseismic tomography besides the recent and remarkable development of the italian seismic network a high model resolution has been achieved improving the inversion method upon the ach method used in previous investigations and picking high quality arrival times with the multi channel cross correlation technique the finer details of our vp model yield new insights into the heterogeneous structure of the adria continental lithosphere involved in the collision between the africa and europe plates a wide low vp anomaly located in the northern adria mantle facing the alpine high vp slab supports the idea that the adria lithosphere has been hydrated and thinned during the alpine subduction we argue that this mantle softening may have played a key role in favoring the subsequent delamination of the adria lithosphere in the northern apennines we hypothesize that delamination of continental lithosphere previously thinned in a back arc setting may be considered a key process to favor subduction polarity reversal and recycling of continental material into the mantle circulation conversely in the central southern apennines the velocity structure is consistent with the existence of a deeper oceanic slab that flattens at the base of the upper mantle in agreement with the widely accepted geodynamic evolution of the central mediterranean by slab retreat and back arc spreading the oceanic slab is discontinuously detached from the surface plate suggesting a different structure of the adria lithosphere which resists subduction instead of favoring delamination

contribution of numeric dynamic modelling to the understanding of the seismotectonic regime of the northern apennines abstract in this paper we investigate the forces possibly active in the area of the northern apennines by means of two dimensional finite element modelling assuming a viscoelastic rheology the forces included in the models are related to continental convergence between africa and eurasia to negative buoyancy of the subducted adriatic lithosphere and to positive buoyancy of anomalously hot mantle material underneath tuscany model predicted stress distribution is compared with seismotectonic data available for the area of the northern apennines and with earthquakes distribution our results indicate that continental convergence cannot be the only mechanism acting in the study area since it causes insignificant subcrustal stresses in disagreement with intermediate seismicity observed under the apennines asthenospheric upwelling under tuscany back arc region is shown to play a crucial role in the present day dynamics of the tyrrhenian apennines positive buoyancy causes an upwards and eastwards flow that generates extensional bending stresses in tuscany and at the top of the slab under the apennines horizontal pushing of the slab is proposed as an alternative mechanism to slab pull to generate subcrustal stresses activation of slab pull does not result in a significant change in orientation of the principal stress axes at shallow depths with respect to the model in which only asthenospheric diapirism is active therefore the existence of a continuous gravitationally sinking slab cannot be ruled out on the basis of comparison with seismotectonic observations both models successfully reproduce extension in tuscany and compression in the outer margin of the apenninic belt but fail to reproduce extension along the inner portion of the chain our modelling highlights the importance of better constraining the lateral extent of the asthenospheric diapir in order to find out whether the extension and tectonic uplift in the apennines is caused by asthenospheric upwelling as in the area of tuscany

lithospheric flexure and thermal anomalies plots of effective elastic thickness te of the oceanic lithosphere as a function of age when elastically deformed show a general increase in te with the square root of age although there is much scatter some of the scatter can be eliminated by correcting for the effect of finite yield strength in rocks which makes te for more sharply bent plates underestimate the true depth tm to the rheological boundary at the base of the high strength mechanical lithosphere after converting effective values te to true thicknesses tm it appears that the mechanical lithosphere is always thinner beneath islands and seamounts compared to lithosphere of similar age flexed at subduction zones from several other lines of evidence such as depth anomalies seamount subsidence and heat flow data it has been suggested that the thermal structure in plates subjected to hot spot type volcanism is reset to values appropriate for much younger lithosphere the reduction in mechanical plate thickness beneath hot spot volcanoes is also consistent with this concept of thermal rejuvenation since the strength of rocks decreases exponentially with rising temperature this dependence implies that variations in mechanical plate thickness can be used to map temperature anomalies in the lithosphere since the mechanical plate thickness depends principally upon the temperature at one depth specifically the depth at which the temperature reaches 550 600 c depending on strain rate flexural data are particularly sensitive indicators of the vertical structure of lithospheric reheating compared to quantities dependent on depth integrals of temperature such as swell height and subsidence the thermal rejuvenation hypothesis predicts that the flexural profile should mimic viscoelastic relaxation in the first 10 ma after loading if initially reheating is confined to depths close to and below the 550 c isotherm the flexural data combined with depth anomaly information for hawaii point to extremely rapid reheating localized in the lower lithosphere rejuvenation in the society islands on the other hand is compatible with slower methods of heat injection involving the entire lithospheric column

subcrustal earthquakes in the northern apennines italy evidence for a still active subduction previously unreported subcrustal earthquakes have been located up to 90 km depth beneath the northern apennines italy these earthquakes occur beneath a zone of abundant upper crustal seismicity mostly confined to the upper 20 km of the crust although there are relatively few well located subcrustal earthquakes there appears to be a general southwestward deepening of the hypocenters from the adriatic to the tyrrhenian sea consistent with other geophysical data which suggest that the adriatic lithosphere is presently subducting beneath the northern apennines

velocities and propagation characteristics of pn and sn beneath the himalayan arc and tibetan plateau possible evidence for underthrusting of indian continental lithosphere beneath tibet new seismological observations on velocities and propagation characteristics of pn and sn waves beneath tibet can be interpreted although not uniquely to indicate the shallow angle underthrusting of the greater indian continental lithosphere beneath the tibetan plateau this is inferred from available geological and geophysical data and is now supported by new seismological observations the most significant observation is that high frequency sn waves shear waves that travel in the mantle lithosphere propagate efficiently in the uppermost mantle beneath the tibetan plateau except beneath the north central part of tibet the chang thang terrane furthermore sn propagates efficiently across the stable block of the tarim basin along the tien shan and himalayan mountains and across the indian shield apparent velocities of pn and sn waves that traverse the uppermost mantle beneath the tibetan plateau are 8 42 and 4 73 km s respectively these velocities are very similar to those beneath the himalayan mountains and the indian shield the efficient propagation and the velocities of pn and sn waves beneath tibet are also similar to those commonly observed beneath shield and stable continental regions these results are consistent with the underthrusting model and do not favor the class of alternative models in which the tibetan plateau is formed by shortening and thickening of hot crustal and uppermost mantle material in response to the convergence of india and asia

the role of membrane like stresses in determining the stability and sensitivity of the antarctic ice sheets back pressure and grounding line motion membrane stresses act along thin bodies which are relatively well lubricated on both surfaces they operate in ice sheets because the bottom is either sliding or is much less viscous than the top owing to stress and heat softening of the basal ice ice streams flow over very well lubricated beds and are restrained at their sides the ideal of the perfectly slippery bed is considered in this paper and the propagation of mechanical effects along an ice stream considered by applying spatially varying horizontal body forces propagation distances depend sensitively on the rheological index and can be very large for ice type rheologies a new analytical solution for ice shelf profiles and grounded tractionless stream profiles is presented which show blow up of the profile in a finite distance upstream at locations where the flux is non zero this is a feature of an earlier analytical solution for a floating shelf the length scale of decay of membrane stresses from the grounding line is investigated through scale analysis in ice sheets such effects decay over distances of several tens of kilometres creating a vertical boundary layer between sheet flow and shelf flow where membrane stresses adjust bounded physically reasonable steady surface profiles only exist conditionally in this boundary layer where bounded steady profiles exist adjacent profile equilibria for the whole ice sheet corresponding to different grounded areas occur neutral equilibrium if no solution in the boundary layer can exist the ice sheet profile must change the conditions for existence can be written in terms of whether the basal rate factor sliding or internal deformation is too large to permit a steady solution the critical value depends extremely sensitively on ice velocity and the back stress applied at the grounding line high ice velocity and high stress both favour the existence of solutions and stability changes in these parameters can cause the steady solution existence criterion to be traversed and the ice sheet dynamics to change a finite difference model which represents both neutral equilibrium and the dynamical transition is presented and preliminary investigations into its numerical sensitivity are carried out evidence for the existence of a long wavelength instability is presented through the solution of a numerical eigenproblem which will hamper predictability

shear wave attenuation in the lithosphere beneath italy and surrounding regions tectonic implications more than 700 waveforms produced by 51 shallow earthquakes and recorded at regional distances 250 1400 km by the italian seismic network have been analyzed to characterize the amplitude and frequency contents of the crustal and uppermost mantle shear waves lg and sn respectively the crustal phase lg efficiently propagates through the relatively stable adriatic continental crust while it is not observed along propagation paths crossing major physiographic features like the apennine chain and the basinal domain of the tyrrhenian and ionian seas similar to lg efficient sn propagation is observed in the uppermost mantle beneath the po plain and the adriatic sea efficient sn transmission is also observed across the northern ionian sea and sicily and in the area between sardinia and the northern coasts of africa sn are efficiently transmitted across the sicily channel and rather efficient sn propagate beneath the ligurian sea on the contrary inefficient sn transmission characterizes the uppermost mantle beneath the apennines the western margin of the italian peninsula and the southern tyrrhenian sea shear wave attenuation suggests the presence of asthenospheric material in the uppermost mantle probably related to the present day extension along the apennine chain and in the tyrrhenian basin this interpretation is consistent with the presence of extensive neogene and quaternary volcanic activity in these areas and related high heat flow proposed lithospheric delaminatioh processes beneath the apennines and subduction beneath the tyrrhenian sea can reasonably explain the observed high attenuation zones in the uppermost mantle in contrast a high strength mantle lid is inferred to underlay the po plain the adriatic sea and the northern ionian sea the available waveforms also indicate that a continuous mantle lid is present beneath sicily and the extensional domain of the sicily channel as well as in the marine area south of sardinia

thermal evolution of a rift basin the tyrrhenian sea existing geological and geophysical data support the view that the opening of the tyrrhenian sea was by rifting of a formerly continuous continental lithosphere which began in upper miocene on the upper sardinian margin and migrated southeastward to the calabria margin heat flow across the tyrrhenian sea shows a pronounced asymmetry from about 50 70 mw m 2 over the sardinian margin to more than 150 mw m 2 over its southeastern margin off southern italy and calabria in this study we use the newly acquired data from ocean drilling program leg 107 heat flow seismic velocity and bathymetric data to constrain a group of models that supposes the tyrrhenian sea opened by detachment of a formerly continuous lithosphere along a southeastern facing low angle normal fault specific models tested include 1 the wernicke model in which thinning and extension of the continental lithosphere are supposed to be accomplished by simple shear along a large scale gently dipping detachment fault zone which cut through the entire lithosphere and 2 the delamination model in which the low angle detachment fault is supposed to cut only the upper and middle crust but to merge subhorizontally into the lower crust below which concurrent pure shear may occur a versatile two dimensional finite element procedure is used to evaluate the tectonic and thermal evolution of these models and the results are compared with observations to test the validity of the models delamination models with large pure shear components in the lower lithosphere and the wernicke model predict heat flow within 1 of the measurements the wernicke model however predicts a basin configuration much too deep in comparison with bathymetry and a subbasinal lithospheric thickness substantially different from rayleigh wave dispersion results the delamination model with a large component of pure shear in the subbasinal lower lithosphere on the other hand predicts satisfactorily both the bathymetry and the lithospheric thickness beneath the basin

lateral variations of pn propagation in italy evidence for a high attenuation zone beneath the apennines pn phases recorded by 40 stations of the italian seismic network are analyzed using the spectral ratio technique to estimate the q structure of the uppermost mantle beneath the italian peninsula and nearby adriatic sea a total of 344 digital waveforms are analyzed from 22 earthquakes that occurred within distances of 300 to 1600 km the calculated apparent q values are grouped into two categories q 800 characterizes the adriatic side of the italian peninsula indicating that the adriatic lithosphere is very efficient in propagating pn phases through the uppermost mantle q 600 characterizes the uppermost mantle beneath the apennines and western italy indicating less efficient wave propagation the presence of asthenospheric mantle material at shallow depths beneath the apennines can explain the observed q

a transportable mb lg scale for central europe and implications for low magnitude ms mb discrimination summary we have extended the mb lg method of nuttli using root mean square rms amplitudes corrected for noise and a 1 dependence for geometrical spreading lg waves recorded on the german regional seismic network grsn for earthquakes in south central europe were used to develop an mb lg formula requiring a new calibration constant crms to keep rms mb lg on the same baseline as nuttli s traditional formula based on 3rd peak amplitudes grsn stations had to be calibrated for site terms and for lg attenuation lateral variations in lgq appear to be significant across the study area and a regional q model consisting of constant q partitions north south and in the central alps was developed using measurements based on interstation and two event single station methods when plotted against surface wave estimates of mw rms mb lg measurements in central europe are found to be consistent with mw mb lg relationships for north america and southern asia thus supporting the transportability of our mb lg formula frequency wavenumber processing of grafenberg array data enabled us to extract rayleigh waves for small events and regional ms were measured using the marshall and basham formula our ms mb lg relationship extends to ms 2 5 and agrees well with observations in other regions including the western united states the discrimination potential of ms mb lg observations was examined under realistic monitoring conditions where path corrections were inferred from earthquake data and applied uniformly to natural sources and explosions under these conditions mb pn p are greater than mb lg for large nts explosions however ms mb scaling slopes are steeper for p waves than they are for lg and ms mb observations for nts explosions converge near mb 4 thus allowing for measurement errors and additional uncertainty in mb pn due to regional bias there is little difference in the discrimination potential for pn and lg waves at small magnitudes as such a regional ms mb discriminant based on lg might be preferred owing to the better detectability of lg waves for small earthquakes these results need to be confirmed for explosions at other test sites compared to teleseismic experience regional ms mb observations extend the discrimination capability to lower magnitudes by at least one ms unit

high frequency seismic wave propagation beneath the indian shield himalayan arc tibetan plateau and surrounding regions high uppermost mantle velocities and efficient sn propagation beneath tibet summary this paper reports the regional propagation characteristics of high frequhcy p p s and l waves 0 5 to 2hz and p pg and s velocities for several regions in central asia tibet and india an important result is that s waves propagate efficiently in the uppermost mantle beneath most of the tibetan plateau also s transmission in tibet is not affected by crossing many physiographic and tectonic boundaries such as the indus tsangpo suture these observations are similar to observations of s propagation in stable continental and shield regions efficient s propagation is also observed across the stable tarim basin along the tien shan and himalayan mountains and across the indian shield s is strongly attenuated across the northern tibetan plateau chang thang terrane such strong attenuation suggests the existence of a low q zone in the uppermost mantle beneath that region this interpretation is supported by recent geological observations of widespread quaternary to recent basaltic volcanism in the chang thang terrane l propagates efficiently across the tarim basin the indian shield and along the strike of the himalayan mountains l is also observed although less pronounced for paths crossing large portions of the complicated structures of tien shan karakoram and pamir l is not observed for paths crossing the tibetan plateau the southern boundary of the region of inefficient l propagation in tibet lies in the vicinity of the indus tsangpo suture zone the lack of l can be attributed to attenuation within the tibetan crust however other possibilities including scattering by structural discontinuities at the margins of tibet are equally plausible pg and s are the most prominent features on most short period wwssn seismograms at epicentral distances less than about 4 apparent pg velocity determined for the area of the himalayan mountains is 6 0 0 2 km s l apparent velocities of p and s waves that have at least 50 per cent of their paths in the uppermost mantle beneath the tibetan plateau are 8 42 it 0 10 and 4 73 f 0 06 km s respectively these velocities are quite similar to those determined for the himalayan mountains p 8 45 f 0 08 km s i s 4 74 0 05 km s and the indian shield p 8 40 0 08 km s s 4 69 k 0 04 km s p