scale dependence of seismic energy to moment ratio for strike slip earthquakes in japan we analyzed four pairs of a large m w 6 and a small m w 3 5 to 4 shallow strike slip earthquakes to investigate the scale dependence of the seismic energy to moment ratio an important macroscopic parameter which reflects the basic physical process of seismic slip these earthquakes occurred in the south western part of japan and high quality close in records epicentral distance 50 km are available for both the small and large earthquakes the paired events have almost the same focal mechanism and hypocenter location we used the spectral ratio of the paired events in order to remove the effects of attenuation along the wave propagation path and the station site response we then estimated the seismic energy from the source spectra estimated from the spectral ratio the energy to moment ratio increases with the earthquake size this scale dependence is very similar to that found earlier for earthquakes in southern california

geochemical and sr nd isotopic constraints on the genesis of the cenozoic linzizong volcanic successions southern tibet abstract the linzizong volcanic successions that crop out in the lhasa terrane southern tibet have been conventionally regarded as the products of northward subduction of the neotethyan oceanic slab beneath south asia this study reports geochemical data of 100 volcanic rocks from the lhasa terrane to better constrain the temporal spatial distribution and petrogenesis of the linzizong volcanism the linzizong volcanic successions consist dominantly of calc alkaline rocks that erupted from ca 69 to 43 ma and show typical arc lava geochemical features marked with lile enrichment and hfse depletion in the spidergram their sr and nd isotope ratios e nd t 3 3 to 2 4 i sr 0 7048 0 7072 are generally similar to those of the associated gangdese i type granitoids the linzizong volcanism is characterized by a flare up period ca 50 ma that shows significant geochemical variations manifested by the coexistence of five types of volcanic rocks 1 the main suite of calc alkaline rocks sio 2 45 80 wt la 12 45 ppm e nd t 3 8 to 4 9 i sr 0 7037 0 7105 that could be interpreted by partial melting of the mantle wedge followed by assimilation and fractional crystallization afc processes with 10 crustal contamination and or magma mixing 2 the low k low ree suite sio 2 48 61 wt k 2 o 0 5 1 1 wt la 7 10 ppm that has the highest nd isotope ratios e nd t 5 9 to 3 5 suggesting a juvenile mantle origin possibly related to decompression melting of the asthenosphere 3 the shoshonitic suite sio 2 53 71 wt k 2 o 3 8 6 7 wt e nd t 2 8 to 6 1 from small degree melting of the metasomatized lithospheric mantle 4 the high ree suite sio 2 61 75 wt la 36 87 ppm e nd t 2 6 to 3 2 originating from remelting of a newly underplated basaltic lower crust and 5 the evolved suite of rhyolitic flows and ignimbrites sio 2 65 wt that have the lowest e nd t values of 14 to 18 representing remelting products of the basement or continental crust of the lhasa terrane such geochemical heterogeneities are attributed to breakoff of the subducted neotethyan slab under southern tibet that occurred in the early stage of the india asia collision

modelling the generation of gravity waves by a maritime continent thunderstorm observed wind and temperature profiles are used to initialize a model calculation of a maritime continent thunderstorm and the numerical solution is used to explore the effect of tropospheric shear on the gravity wave generation the resultant convective system is qualitatively similar to that observed the modelled gravity waves propagate away from the cloud with wave fronts that are approximately circular implying that the convective clouds do not generate waves which propagate in a preferred direction the gravity wave generation is related to the oscillation of the convective updraughts about their level of neutral buoyancy while the general features of the gravity waves are similar in most respects to a previous study which used an idealized wind profile the frequency of the gravity waves is doppler shifted by the tropospheric wind shear the result is a much broader power spectrum in comparison to the idealized cases

global gravity wave variances from aura mls characteristics and interpretation abstract the gravity wave gw resolving capabilities of 118 ghz saturated thermal radiances acquired throughout the stratosphere by the microwave limb sounder mls on the aura satellite are investigated and initial results presented because the saturated optically thick radiances resolve gw perturbations from a given altitude at different horizontal locations variances are evaluated at 12 pressure altitudes between 21 and 51 km using the 40 saturated radiances found at the bottom of each limb scan forward modeling simulations show that these variances are controlled mostly by gws with vertical wavelengths z 5 km and horizontal along track wavelengths of y 100 200 km the tilted cigar shaped three dimensional weighting functions yield highly selective responses to gws of high intrinsic frequency that propagate toward the instrument the latter property is used to infer the net meridional component of gw propagation by differencing the variances acquired from ascending a and descending d o

vertical velocity characteristics of deep convection over darwin australia abstract continuous vertical velocity measurements using a 50 mhz wind profiler located at darwin in northern australia during periods of active convection have been analyzed this dataset is dominated by continental type convection numerous examples of shallow deep and decaying convection were seen and it is shown that only the deep systems have substantial tilts to the draft structure the most intense updrafts occur above the freezing level but shallow convection also produces large amplitude vertical motions the strength of these updrafts in this dataset is very similar to other tropical oceanic data that observation is consistent with the idea that the magnitude of the updrafts is much less in the tropics than for intense midlatitude convection because the convective available potential energy is distributed over a much deeper layer in the tropics although more intense updrafts may be present at other tropical locations such as the tiwi islands north of darwin the size of the cores however

acceleration spectra for subduction zone earthquakes we estimate the source spectra of shallow earthquakes from digital recordings of teleseismic p wave groups that is p pp sp by making frequency dependent corrections for the attenuation and for the interference of the free surface the correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths we apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from ms 6 2 to 8 1 obtaining corrected p wave acceleration spectra on the frequency band from 0 01 to 2 0 hz seismic moment estimates from surface waves and normal modes are used to extend these p wave spectra to the frequency band from 0 001 to 0 01 hz the acceleration spectra of moderate subduction zone earthquakes that is earthquakes whose seismic moments are less than 1027 dyn cm exhibit square or brune type spectra while the acceleration spectra of large subduction zone earthquakes that is earthquakes whose seismic moments are greater than 1027 dyn cm exhibit intermediate slopes where u 5 4 for frequencies from 0 005 to 0 05 hz for this set of earthquakes spectral shape appears to be a discontinuous function of seismic moment using reasonable assumptions for the phase characteristics we transform the spectral shape observed for large earthquakes into the time domain to fit ekstrom s 1987 moment rate functions for the ms 8 1 michoacan earthquake of september 19 1985 and the ms 7 6 michoacan aftershock of september 21 1985

a global morphology of gravity wave activity in the stratosphere revealed by the gps occultation data gps met using temperature profiles obtained by the gps met gps meteorology experiment from april 1995 to february 1997 we have extracted mesoscale temperature perturbations with vertical wavelengths ranging from 2 to 10 km and background brunt vaisala frequency squared n2 for each occultation event we can evaluate a potential energy ep which is assumed to be caused by atmospheric gravity waves the monthly mean values of ep at 15 20 km around japan showed an annual variation with an enhancement in winter which is consistent with the climatological behavior of the kinetic energy of gravity waves observed with the mu middle and upper atmosphere radar 34 9 n 136 0 e from 1985 to 1989 we have then derived the global distribution of ep at 20 30 km during northern hemisphere winter from november to february our analysis shows that the largest ep values are generally centered around the equator between 25 n and 25 s with considerable longitude variations longitudinal variations of ep at 20 30 km in a latitude range of 30 60 n are also analyzed resulting in larger ep values over the continents than over the pacific ocean using gps met data without antispoofing latitudinal variations of ep are determined in 15 45 km although large ep values are concentrated near the equator at 20 30 km ep tends to become larger at midlatitudes at 30 40 km and higher altitude regions at midlatitudes ep is found to be larger in winter months in both hemispheres height variations of ep indicate a decrease at 25 30 km and a monotonic increase above 30 km

tibetan tectonic evolution inferred from spatial and temporal variations in post collisional magmatism cenozoic magmatism on the tibetan plateau shows systematic variations in space and time that must be considered in models concerning tibetan tectonic evolution after the india asia collision which started in the early tertiary and terminated the gangdese arc magmatism in the lhasa terrane made of the southern tibetan plateau widespread potassium rich lavas and subordinate sodium rich basalts were generated from f50 to 30 ma in the qiangtang terrane of northern tibet subsequent post collisional magmatism migrated southwards producing ultrapotassic and adakitic lavas coevally between f26 and 10 ma in the lhasa terrane then potassic volcanism was renewed to the north and has become extensive and semicontinuous since f13 ma in the western qiangtang and songpan ganze terranes such spatial temporal variations enable us to elaborate a geodynamic evolution model which depicts when and how the indian continental lithospheric mantle started thrusting under asia by involving rollback and breakoff of the subducted neo tethyan slab followed by removal of the thickened lhasa lithospheric root we propose that only after the lithospheric removal which occurred at f26 ma could the indian mantle lithosphere have commenced its northward underthrusting and henceforth served as a pivotal control to the himalayan tibetan orogenesis consequently the tibetan plateau is suggested to have risen diachronously from south to north whereas the southern part of the plateau may have been created and maintained since the late oligocene the northern plateau would have not attained its present day elevation and size until the mid miocene when the lower part of the western qiangtang and songpan ganze lithospheres began to founder owing to the push of the underthrust indian mantle lithosphere d 2004 elsevier b v all rights reserved

ubiquitous weakening of faults due to thermal pressurization faults weaken during earthquakes analysis of the amount of energy released during earthquakes globally suggests that heat induced pressurization of pore fluids can weaken faults during earthquakes of all sizes laboratory simulations of earthquakes show that at high slip rates faults can weaken significantly aiding rupture1 2 3 various mechanisms such as thermal pressurization and flash heating have been proposed to cause this weakening during laboratory experiments1 4 5 6 yet the processes that aid fault slip in nature remain unknown measurements of seismic radiation during an earthquake can be used to estimate the frictional work associated with fault weakening known as an event s fracture energy7 8 9 here we compile new and existing8 9 measurements of fracture energy for earthquakes globally that vary in size from borehole microseismicity to great earthquakes we observe a distinct transition in how fracture energy scales with event size which implies that faults weaken differently during small and large earthquakes and earthquakes are not self similar we use an elastodynamic numerical model of earthquake rupture to explore possible mechanisms we find that thermal pressurization of pore fluid by the rapid shear heating of fault gouge can account for the observed scaling of fracture energy in small and large earthquakes over seven orders of fault slip magnitude we conclude that thermal pressurization is a widespread and prominent process for fault weakening

first report of lower permian basalts in south tibet tholeiitic magmatism during break up and incipient opening of neotethys abstract the neo tethys ocean began to form at early permian times when continental flood basalts were emplaced in various areas of the newly formed indian passive margin exposed today in the so called tibetan sedimentary zone of the himalaya lower permian mafic volcanic rocks which have long been known from various himalayan localities from kashmir to arunachal pradesh are here for the first time reported to occur also in south tibet bhote kosi basalts of the gyirong county the basalts unconformably overlie lowermost permian diamictites with locally intervening black shales and debris flow deposits and are followed in turn by chert bearing quartzarenites and silty to phosphatic marls yielding brachiopods of roadian wordian age the age of the lavas can thus be bracketed as late early permian post sakmarian and pre roadian the geochemistry of these subalkalic tholeiites akin to morbs testifies to their similarity not only with the adjacent nar tsum spilites of central nepal but also with the panjal traps and abor volcanics of the western and eastern himalayas respectively the geochemical signature of lower permian volcanic rocks is in fact uniform all along the himalayan range and markedly different from that of basaltic rhyolitic alkalic products sporadically emplaced during the previous rifting stage rift volcanism in the tethys himalaya began in the early carboniferous and came to an end in sakmarian times in the early permian initial submergence of the rift shoulders and sediment starvation were followed by tholeiitic magmatism which is therefore interpreted as following break up and incipient sea floor spreading in the neotethys ocean roughly contemporaneous emplacement of continental flood basalts of similar geochemical signature along a 2000 km long rift axis would in fact suggest extensive mantle melting at the transition from continental rifting to break up and opening of the neotethys between northern gondwana and the peri gondwanian blocks

site amplifications for generic rock sites seismic shear wave velocity as a function of depth for generic rock sites has been estimated from borehole data and studies of crustal velocities and these velocities have been used to compute frequency dependent amplifications for zero attenuation for use in simulations of strong ground motion we define a genetic rock site as one whose velocity at shallow depths equals the average of those from the rock sites sampled by the borehole data most of the boreholes are in populated areas for that reason the rock sites sampled are of particular engineering significance we consider two genetic rock sites rock corresponding to the bulk of the borehole data and very hard rock such as is found in glaciated regions in large areas of eastern north america or in portions of western north america the amplifications on rock sites can be in excess of 3 5 at high frequencies in contrast to the amplifications of less than 1 2 on very hard rock sites the consideration of unattenuated amplification alone is computationally convenient but what matters for ground motion estimation is the combined effect of amplification and attenuation for reasonable values of the attenuation parameter k 0 the combined effect of attenuation and amplification for rock sites peaks between about 2 and 5 hz with a maximum level of less than 1 8 the combined effect is about a factor of 1 5 at 1 hz and is less than unity for frequencies in the range of 10 to 20 hz depending on t using these amplifications we find provisional values of about ao 70 bars and k 0 0 035 sec for rock sites in western north america by fitting our empirically determined response spectra for an m 6 5 event to simulated values the borehole data yield shear velocities 033o of 618 and 306 rrdsec for rock and soil sites respectively when averaged over the upper 30 m from this we recommend that v30 equals 620 and 310 m sec for applications requiring the average velocity for rock and soil sites in western north america by combining the amplifications for rock sites and the site factors obtained from our analysis of strong motion data we derive amplifications for sites with 17 3o 520 m sec nehrp class c corresponding to a mix of rock and soil sites and v3o 310 and 255 m sec average soil and nehrp class d sites respectively for the average soil site the combined effect of amplification and attenuation exceeds a factor of 2 0 for frequencies between 0 4 and about 4 hz with a peak site effect of 2 4 the peak of the nehrp class d site effect is 2 7

strike slip earthquakes in the oceanic lithosphere observations of exceptionally high apparent stress summary the radiated energies e s and seismic moments m0 for 942 globally distributed earthquakes that occurred between 1987 to 1998 are examined to find the earthquakes with the highest apparent stresses a e s m0 where is the modulus of rigidity the globally averaged a for shallow earthquakes in all tectonic environments and seismic regions is 0 3 mpa however the subset of 49 earthquakes with the highest apparent stresses a greater than about 5 0 mpa is dominated almost exclusively by strike slip earthquakes that occur in oceanic environments these earthquakes are all located in the depth range 7 29 km in the upper mantle of the young oceanic lithosphere many of these events occur near plate boundary triple junctions where there appear to be high rates of intraplate deformation indeed the small rapidly deforming gorda plate accounts for 10 of the 49 high a events the depth distribution of a which shows peak values somewhat greater than 25 mpa in the depth range 20 25 km suggests that upper bounds on this parameter are a result of the strength of the oceanic lithosphere a recently proposed envelope for apparent stress derived by taking 6 per cent of the strength inferred from laboratory experiments for young less than 30 ma deforming oceanic lithosphere agrees well with the upper bound envelope of apparent stresses over the depth range 5 30 km the corresponding depth dependent shear strength for young oceanic lithosphere attains a peak value of about 575 mpa at a depth of 21 km and then diminishes rapidly as the depth increases in addition to their high apparent stresses which suggest that the strength of the young oceanic lithosphere is highest in the depth range 10 30 km our set of high a earthquakes show other features that constrain the nature of the forces that cause interplate motion first our set of events is divided roughly equally between intraplate and transform faulting with similar depth distributions of a for the two types secondly many of the intraplate events have focal mechanisms with the t axes that are normal to the nearest ridge crest or subduction zone and p axes that are normal to the proximate transform fault these observations suggest that forces associated with the reorganization of plate boundaries play an important role in causing high a earthquakes inside oceanic plates extant transform boundaries may be misaligned with current plate motion to accommodate current plate motion the pre existing plate boundaries would have to be subjected to large horizontal transform push forces a notable example of this is the triple junction near which the second large aftershock of the 1992 april cape mendocino california sequence occurred alternatively subduction zone resistance may be enhanced by the collision of a buoyant lithosphere a process that also markedly increases the horizontal stress a notable example of this is the aleutian trench near which large events occurred in the gulf of alaska in late 1987 and the 1998 march balleny sea m 8 2 earthquake within the antarctic plate

subsidence history of the tethys himalaya abstract this article reconstructs the evolution of the passive northern margin of peninsular india facing the neotethys from the earliest rifting stages in the carboniferous to final collision with the active southern margin of asia in the early eocene classical techniques of basin analysis were applied to an extensive and coherent stratigraphic data base facies analysis biostratigraphic dating and palaeoenvironmental interpretation of an over 10 km thick sedimentary succession from eight main composite sections allowed us to draw accurate sedimentation and backstripped subsidence curves and to place quantitative constraints on the palaeotectonic scenario inferred from the sedimentary record in each of the investigated composite sections three major stratigraphic gaps are consistently recognised in the rift sequence the rift unconformity mississippian the carboniferous permian hiatus largely kasimovian to asselian and the break up unconformity mid sakmarian two or even all three gaps may merge into a single hiatus spanning up to 70 ma overall gaps are associated with rift shoulder uplift which after calculation of backstripped subsidence sums up to 600 m at least in three of the investigated sections thermal subsidence and low sedimentation rates are documented by the mid permian to middle triassic drift sequence deposited while the neotethys was expanding between northern gondwana and the detached peri gondwanan blocks a marked increase in accumulation rates is recorded in the carnian norian but calculations of uniform lithospheric stretching suggest that extension was minor and associated with tectonic processes affecting distant parts of gondwana the tethys himalayan margin was uplifted around the jurassic cretaceous boundary at the onset of a major volcanic episode finally uplift exceeding in magnitude all previous events is documented by the uppermost units of the stratigraphic column deposited during the earliest stages of collision between india and asia in the latest paleocene to early eocene

effect of directivity on estimates of radiated seismic energy to accurately estimate the seismic energy radiated in an earthquake it is important to use appropriate corrections for path and source effects using slip models obtained by inversion of seismic data we examine the effect of directivity on estimates of radiated energy and develop a method to correct for this source effect from our calculations we suggest that the directivity correction for the earthquakes we studied is less than a factor of three we also conclude that the directivity correction for dip slip earthquakes with rupture propagating along strike is less than a factor of two at teleseismic distances thus teleseismic energy estimates of large subduction zone earthquakes will not be significantly affected by directivity

the andean type gangdese mountains paleoelevation record from the paleocene eocene linzhou basin abstract paleoelevation reconstruction using oxygen isotopes is making a significant contribution to understanding the cenozoic uplift of the himalayas and the tibetan plateau this paper presents new oxygen and carbon isotopic compositions from well dated tertiary paleosols lacustrine calcareous carbonates and marls from the nianbo 60 54 ma and upper pana formations 51 48 ma of the linzizong group in the linzhou penbo basin the sediments of the nianbo formation which are 180 m thick were deposited in alluvial fans braided rivers fan deltas and on nearshore to offshore lacustrine settings whereas those of the upper pana formation are 100 m thick and are comprised predominantly of proximal alluvial fan and braided river deposits correlations between the lithofacies and stable isotopic compositions suggest that the basin was mainly a hydrologically open environment it is confirmed that the 18 o c and 13 c c values from nianbo and pana formations have not yet been reset by late stage diagenesis based on petrographic examination oxygen isotope of the fossil ostracodes and tectonic deformation of strata the paleoelevations are reconstructed using the corrected most negative paleosurface water 18 o psw values these imply that the linzhou area had attained an elevation of 4500 400 m during the period of the indo asian collision i e achieved a near present elevation and may form an andean type mountain range stretching the gangdese arc before collision the gangdese mountains probably maintained high elevations since at least the paleocene and could play a crucial role in the climate change in the interior of the tibetan plateau during the early cenozoic the paleogeomorphic scenario of the eocene tibet is proposed to exist at two high mountains in excess of 4500 m that sandwiched a low elevation basin

late cretaceous to middle tertiary basin evolution in the central tibetan plateau changing environments in response to tectonic partitioning aridification and regional elevation gain located in the bangong suture zone between the lhasa and qiangtang terranes of central tibet the nima basin records cretaceous late miocene sediment accumulation during a period of drastically changing paleogeography and paleoelevation the bangong suture zone originally formed during late jurassic early cretaceous time as the lhasa terrane collided with the qiang tang terrane during early to middle cretaceous time the region lay at northern near equatorial paleolatitudes near sea level by aptian time the nima basin was above sea level and was strongly influenced by nearby volcanic activity and crustal shortening in the reactivated bangong suture zone in the southern nima basin an 50 m y late cretaceous eocene depositional hiatus correlates with major crustal shortening and ensuing voluminous ignimbrite eruptions in the lhasa terrane in the northern nima basin deposition continued during latest cretaceous time recording arid paleoclimate in evaporitic lacustrine and eolian dunefield deposits by oligocene time the nima basin comprised two independent depocenters that accumulated coarse grained alluvial fluvial lacustrine evaporitic and fan delta deposits in close association with reactivated thrusts in the bangong suture zone carbon and oxygen isotope data from oligocene paleosol carbonate reported elsewhere indicate that regional paleoelevation during the late oligocene was 4 6 km as it is today overall the depositional record of the nima basin combined with ongoing structural and geochronologic studies demonstrates that the bangong suture zone was reactivated during middle cretaceous and middle tertiary time that the intervening 50 m y interval was a time of regional upper crustal shortening in the lhasa terrane followed by regional ignimbrite eruptions and that arid paleoclimate and high paleoelevation were established by the late cretaceous and late oligocene respectively within the context of other data sets from the lhasa terrane the record of deposition in the nima basin is consistent with low angle subduction of neotethyan oceanic lithosphere and reactivation of the bangong suture zone during the early cretaceous followed by shortening within the gangdese retroarc and northern lhasa terrane thrust belts during the late cretaceous lithospheric delamination dripping and regional magmatic flare up during latest cretaceous through early tertiary time and underthrusting of indian lower crust and lithosphere as far north as the bangong suture zone during late oligocene time

transition from oceanic to continental lithosphere subduction in southern tibet evidence from the late cretaceous early oligocene 91 30 ma intrusive rocks in the chanang zedong area southern gangdese little is known about the detailed processes associated with the transition from oceanic to continental lithosphere subduction in the gangdese belt of southern tibet gbst here we report zircon u pb age major and trace element and sr nd hf isotopic data for late cretaceous early oligocene similar to 91 30 ma intermediate acid intrusive rocks in the chanang zedong area immediately north of the yarlung tsangpo suture zone these rocks represent five magmatic episodes at similar to 91 similar to 77 similar to 62 similar to 48 and similar to 30 ma respectively the 91 48 ma rocks have slightly lower initial sr 87 sr 86 0 7037 to 0 7047 and higher epsilon nd t 1 8 to 4 3 and epsilon hf t 3 5 to 14 7 values in comparison with those 0 7057 to 0 7062 3 3 to 2 5 and 2 2 to 6 6 of the similar to 30 ma intrusive rocks the similar to 91 similar to 62 and similar to 30 ma rocks are geochemically similar to slab derived adakites the similar to 91 ma somka adakitic granodiorites were likely derived by partial melting of the subducting neo tethyan oceanic crust with minor oceanic sediments and the similar to 91 ma somka dioritic rocks with a geochemical affinity of adakitic magnesian andesites likely resulted from interactions between adakitic magmas and overlying mantle wedge peridotite the similar to 77 ma luomu diorites were probably generated by partial melting of juvenile basaltic lower crust the similar to 62 ma naika and zedong adaldtic diorites and granodiorites were likely generated mainly by partial melting of thickened juvenile mafic lower crust but the source region of the zedong adakitic rocks also contained enriched components corresponding to indian continental crust the similar to 48 ma lamda granites were possibly generated by melting of a juvenile basaltic crust the younger similar to 30 ma chongmuda adakitic quartz monzonites and minor granodiorites were most probably derived by partial melting of early oligocene north ward subducted indian lower crust beneath the southern lhasa block taking into account the regional tectonic and magmatic data we suggest that the gangdese belt of southern tibet gbst underwent a tectonodynamic transition from oceanic subduction to continental subduction between 100 and 30 ma it evolved through four stages 100 65 ma roll back of subducted neo tethyan oceanic lithosphere 65 60 ma initial collision between indian and asian continents 60 40 ma breakoff of subducted neo tethyan oceanic lithosphere and similar to 30 ma northward subduction of the indian continent c 2014 elsevier b v all rights reserved

a preliminary morphology of precipitation systems in tropical northern australia a preliminary morphology of convective systems observed in the vicinity of darwin australia is presented several types of tropical convection during monsoonal and break season flow are identified using specific examples and compared with a range of systems observed worldwide in the general context of the bulk richardson number a significant spectrum of convective activity ranging from low shear and low convective available potential energy cape to high shear and moderate to high cape typical of mid latitude severe storms is identified significant differences between monsoonal flow and break season systems are apparent the formation of convective lines perpendicular to the low level shear and the apparent balance between advancing shallow cold pools and the shear are ubiquitous features convective scale downdraughts and mesoscale descent appear to be responsible for redistribution of 700 mb momentum to the planetary boundary layer pel the motion of the break season squall lines appears to result from an equilibrium between the cold pool and the undisturbed environmental pel flow deviations from this balance are hypothesized to result in the observed propagating and non propagating or slow moving modes examples of these modes are given showing the classical continuous development of new cells in a balanced state along the leading edge of long lived squalls and in the propagating case the discontinuous development of new cells on the cold pool ahead of the squall

the origin and pre cenozoic evolution of the tibetan plateau different hypotheses have been proposed for the origin and pre cenozoic evolution of the tibetan plateau as a result of several collision events between a series of gondwana derived terranes e g qiangtang lhasa and india and asian continent since the early paleozoic this paper reviews and reevaluates these hypotheses in light of new data from tibet including 1 the distribution of major tectonic boundaries and suture zones 2 basement rocks and their sedimentary covers 3 magmatic suites and 4 detrital zircon constraints from paleozoic metasedimentary rocks the western qiangtang amdo and tethyan himalaya terranes have the indian gondwana origin whereas the lhasa terrane shows an australian gondwana affinity the cambrian magmatic record in the lhasa terrane resulted from the subduction of the proto tethyan ocean lithosphere beneath the australian gondwana the newly identified late devonian granitoids in the southern margin of the lhasa terrane may represent an extensional magmatic event associated with its rifting which ultimately resulted in the opening of the songdo tethyan ocean the lhasa northern australia collision at similar to 263 ma was likely responsible for the initiation of a southward dipping subduction of the bangong nujiang tethyan oceanic lithosphere the yarlung zangbo tethyan ocean opened as a back arc basin in the late triassic leading to the separation of the lhasa terrane from northern australia the subsequent northward subduction of the yarlung zangbo tethyan ocean lithosphere beneath the lhasa terrane may have been triggered by the qiangtang lhasa collision in the earliest cretaceous the mafic dike swarms ca 284 ma in the western qiangtang originated from the panjal plume activity that resulted in continental rifting and its separation from the northern indian continent the subsequent collision of the western qiangtang with the eastern qiangtang in the middle triassic was followed by slab breakoff that led to the exhumation of the qiangtang metamorphic rocks this collision may have caused the northward subduction initiation of the bangong nujiang ocean lithosphere beneath the western qiangtang collision related coeval igneous rocks occurring on both sides of the suture zone and the within plate basalt affinity of associated mafic lithologies suggest slab breakoff induced magmatism in a continent continent collision zone this zone may be the site of net continental crust growth as exemplified by the tibetan plateau c 2012 international association for gondwana research published by elsevier b v all rights reserved

mesozoic and cenozoic volcanic rocks from central and southern tibet 39ar 40ar dating petrological characteristics and geodynamical significance field relations radiometric ages and mineralogical and chemical characteristics of the orogenic volcanic rocks from the lhasa block tibet are studied volcanic activity started before the deposition of the albian aptian sedimentary units takena formation and xigaze group volcanism is volumetrically limited in the northern part of the lhasa block basalts andesites and dacites are present while ignimbrites are scarce radiometric ages of these rocks range between 110 and 80 ma conversely in the southern part of the lhasa block volcanism is largely developed gangdese belt there paroxysmal activity occurred later than in the north between 60 and 50 ma lingzizong formation this volcanism is mainly represented by andesites and ignimbrites a younger episode at about 15 10 ma took place in this southern region the mineralogical and geochemical characteristics indicate that as a whole the volcanic rocks from the lhasa block belong to a calc alkaline continental margin series enriched in k20 and other incompatible elements subordinate potassic alkaline lavas are present in the northern margin of the gangdese belt some basalts exhibit a tholeiitic affinity daquin area southern region the geodynamic significance of these cenozoic volcanic rocks lies in the relative motion of india and the lhasa block after accretion of the latter to the northern qantang block a geodynamic model based on changes of the type of subduction of the indian oceanic crust beneath the lhasa block for the period 120 40 ma is proposed to account for the abundance age and magmatic affinity of the volcanic products

teleseismic estimates of the energy radiated by shallow earthquakes the energy flux contained in the p wave groups p pp sp or the s wave groups s ps ss radiated by a shallow earthquake is modeled assuming that the energy flux in the direct and depth phases adds incoherently by defining generalized radiation patterns which incorporate this neutral interference the wave groups are analyzed as though they were comprised of a single phase measurements of the energy flux in the wave groups are corrected in the frequency domain for both the body wave attenuation and the frequency band of the recording the corrected measurements are then used to estimate the seismic energy radiated by the earthquake this analysis is applied to digital recordings of the teleseismic wave groups radiated by the may 2 1983 coalinga california earthquake and the october 28 1983 borah peak idaho earthquake for the coalinga earthquake an estimate of es 1 6 0 4 1021 dyn cm was determined from six p wave groups while the sh wave group recorded at station col returned an estimate of es 1 2 1021 dyn cm for the borah peak earthquake an estimate ofes 3 2 0 5 1021 dyn cm was determined from seven p wave groups the distribution of the isoseismals for the borah peak earthquake indicate that the energy radiated by this event was focussed to the northwest in the direction of rupture propagation correcting the teleseismic estimate for this focussing gives es 3 6 0 5 1021 dyn cm combining these estimates of the radiated energy with broad band estimates of the seismic moment yields estimates of the apparent stress of 17 and 7 bars for the coalinga and the borah peak earthquakes respectively

generation of sodium rich magmas from newly underplated basaltic crust sodium rich rocks of trondhjemite tonalite dacite ttd or granodiorite ttg suites form much of precambrian continental crust1 they are thought to have formed by partial melting of subducted oceanic crust2 3 a process that would have been much more widespread early in earth history than at present owing to the higher thermal gradients prevailing at that time4 phanerozoic ttd suites do exist however and seem also to relate to subduction zones5 defant and drummond6 proposed that these suites form where young 25 myr hot oceanic lithosphere is subducted and melts thus locally simulating the conditions that led to widespread crustal growth in the archaean here we describe plutonic and volcanic rocks from the cordillera blanca complex in peru which have characteristics of the high al ttd suite but which were produced above a subduction zone containing a 60 myr old slab we present evidence that the complex formed by partial melting of newly underplated basaltic crust and argue that this mechanism should be considered more generally as an additional way of generating sodium rich arc magmas

variability of the australian summer monsoon at darwin 1957 1992 abstract the variability of the australian summer monsoon is reexamined using data covering 35 monsoon seasons a new easily applied objective definition of active and break phases of the monsoon based solely on the zonal wind at darwin is proposed attempts to define wet westerly onsets are shown to be misleading since no clear relationship is found between westerly winds and rainfall on the timescales associated with the transition between active and break phases the resulting dates of monsoon onset at darwin differ from those reported in some recent studies resulting in significantly different relationships with the southern oscillation in particular the date of monsoon onset is shown to be significantly related to and hence predicted from prior values of the southern oscillation index also in contrast to a number of recent studies that have highlighted the so called 40 50 day oscillation in the australian summer monsoon no dominant timescales are found in the length of the active period

zircon u pb age and sr nd hf o isotope geochemistry of the paleocene eocene igneous rocks in western gangdese evidence for the timing of neo tethyan slab breakoff northward neo tethyan oceanic lithosphere subduction beneath southern tibet in the mesozoic early cenozoic produced continental arc magmas in the 1600 km long gangdese belt the most voluminous magmatism occurred in the paleocene eocene and is characterized by extensive i type calc alkaline to high k calc alkaline linzizong volcanic rocks and coeval plutons these rocks have been extensively studied in the eastern gangdese belt east of 89 e but few data exist from the western gangdese belt new data for eleven samples of these rocks combined with existing data from the literature show that they are similar to the eastern gangdese belt rocks with relative depletions in nb ta p and ti and enrichments in rb ba th u k pb zr and hf on a primitive mantle normalized trace element diagrams typical of continental arc related igneous rocks however compared to the east western gangdese igneous rocks range to higher k2o contents up to 6 1 wt higher 87sr 86sr i ratios up to 0 7151 and lower endi values down to 8 1 suggesting that an evolved crustal source was involved in arc magmatism the gangdese arc magmatism lasted to 80 ma and has a gap or quiescent period afterwards starting at 69 ma the arc magmatism initiated and shifted southward from 30 5 n to 29 5 n in southern tibet with an abrupt change of india asia convergence rate the magmas through the whole gangdese belt at 69 53 ma are characterized by intermediate endi values 0 6 to 4 0 positive ehfi values 3 8 to 7 1 intermediate 18o values 5 0 to 6 5 and low th y and la yb ratios 20 these magmas were likely derived from the mantle with crustal contamination by mash process in response to neo tethyan slab rollback as proposed by previous studies the sr nd hf o isotopic compositions of linzizong volcanic rocks and coeval intrusions are relatively homogeneous during 69 53 ma but show a wide range at 53 50 ma zircon temperature values from these rocks are distinctively high during the latter period and whole rock th y and la yb ratios increased significantly after 50 ma we suggest that the high temperature and heterogeneous geochemistry of these magmas at 53 50 ma reflect neo tethyan slab breakoff during india asia collision which triggered asthenospheric mantle upwelling mantle derived magmas resulted in extensive heating induced crustal melting and generated magmas with variable isotopic compositions tibetan crust was significantly thickened after 50 ma

differential energy radiation from two earthquakes in japan with identical mw the kyushu 1996 and tottori 2000 earthquakes we examine two closely located earthquakes in japan that had identical moment magnitudes m w but significantly different energy magnitudes m e we use teleseismic data from the global seismograph network and strong motion data from the national research institute for earth science and disaster prevention s k net to analyze the 19 october 1996 kyushu earthquake m w 6 7 m e 6 6 and the 6 october 2000 tottori earthquake m w 6 7 m e 7 4 to obtain regional estimates of radiated energy e s we apply a spectral technique to regional 200 km waveforms that are dominated by s and lg waves for the thrust fault kyushu earthquake we estimate an average regional attenuation q f 230 f 0 65 for the strike slip tottori earthquake the average regional attenuation is q f 180 f 0 6 these attenuation functions are similar to those derived from studies of both california and japan earthquakes the regional estimate of e s for the kyushu earthquake 3 8 1014 j is significantly smaller than that for the tottori earthquake e s 1 3 1015 j these estimates correspond well with the teleseismic estimates of 3 9 1014 j and 1 8 1015 j respectively the apparent stress a e s m with equal to rigidity for the kyushu earthquake is 4 times smaller than the apparent stress for the tottori earthquake in terms of the fault maturity model the significantly greater release of energy by the strike slip tottori earthquake can be related to strong deformation in an immature intraplate setting the relatively lower energy release of the thrust fault kyushu earthquake can be related to rupture on mature faults at a subduction environment the consistence between teleseismic and regional estimates of e s is particularly significant as teleseismic data for computing e s are routinely available for all large earthquakes whereas often there are no near field data

gravity waves in the tropical lower stratosphere an observational study of seasonal and interannual variability radiosonde observations made at cocos islands 12 s 97 e in the indian ocean between september 1992 and june 1998 are used to study seasonal and interannual variations in gravity wave activity in the lower stratosphere 18 25 km the islands are located in a region of generally strong convection that occurs at all times of the year with the period of strongest convective activity between december and july wet season the prevailing zonal winds during the observational period and height range are westward with a quasi biennial oscillation qbo superimposed time series of wave energy show that largest wave amplitudes occur during the wet season when convection is strongest but a qbo like variation is also apparent maximum energy densities of about 25 j kg 1 occur early 1993 1995 and 1997 at the times when the westward shears are largest wave energy is found to be propagating upward and in the horizontal there is considerable azimuthal anisotropy with predominate eastward propagation against the prevailing wind upward fluxes of zonal momentum flux u w are estimated by combining the temperature and wind information fluxes show a similar temporal behavior to the energy the motion and temperature fields are dominated by waves with vertical wavelengths 2 km using a stokes parameter analysis it is found that the intrinsic frequencies are on average 2 3 times the inertial frequency corresponding to intrinsic periods of 20 25 hours horizontal wavelengths between 200 and 2000 km are inferred with a mean value of about 1000 km the mean intrinsic phase speeds are about 10 ms 1 but ground based phase speeds are centered on 0 ms 1

late cretaceous gangdese intrusions of adakitic geochemical characteristics se tibet petrogenesis and tectonic implications article i nfo the gangdese batholith emplaced from the cretaceous to eocene in southern tibet has been widely regarded as the major constituent of an andean type convergent margin resulting from northward subduction of the neo tethyan oceanic lithosphere under asia while the gangdese batholith consists predominantly of calc alkaline rocks we identify from the eastern part of the batholith a suite of epidote bearing granodiorites that shows adakitic geochemical characteristics marked with apparently higher la yb and sr y and lower y and hree than other gangdese rocks and common arc magmas shrimp zircon u pb analyses of two of the samples yielded 206 pb 238 u dates of 80 4 1 1 and 82 7 1 6 2 ma which constrain the emplacement ages of the adakitic rocks trace element modeling suggests that these rocks originated from partial melting of a garnet amphibolite source that on the basis of the sr and nd isotope data isr 0 7044 0 7048 end t 3 2 to 0 9 we interpret to be a newly underplated mafic lower crust rather than the subducted neo tethyan oceanic crust this juvenile crust was produced by cretaceous basaltic underplating above the mantle wedge and then thickened by the tectonic contraction owing to flattening of the neo tethyan subduction a process that also led to the adakitic magmatism our interpretation involving a late cretaceous flat slab subduction and related orogenesis in southern tibet is consistent with petrographic data such as the occurrence of magmatic epidote and muscovite rimmed with resorption texture in the granodiorites which indicate deep seated emplacement followed by rapid tectonic exhumation

geochemical criteria for a refined tectonic discrimination of orogenic andesites abstract minor trace element and sr isotope data are compiled for tertiary recent orogenic andesites of three tectonic settings 1 oceanic island arcs divided into low k and other andesites 2 continental island arcs and thin continental margins and 3 thick continental margins andean from 1 to 3 there is a steady overlapping increase in la ce ree la yb la y th th u zr y hf yb and ni co and a decrease in k la p la sc cr and sc ni optimum discriminationis achieved with a plots of la yb vs th and la yb vs sc ni and b cluster analysis using the critical parameters listed above shoshonitic andesites and high mg andesites boninites are considered separately orogenic andesites are distinguished from non orogenic andesites by high contents of al low contents of ti zr ree y nb ta and ga and the absence of pronounced negative eu anomalies

analysis of the ground accelerations radiated by the 1980 livermore valley earthquakes for directivity and dynamic source characteristics the strong motion accelerograph recordings of the 24 january 1980 main shock and the 27 january 1980 aftershock of the livermore valley earthquake sequence are analyzed for systematic variations with azimuth or station location the variation of the peak accelerations with epicentral azimuth is apparently reversed for the two events the main shock accelerations are larger to the south and the aftershock accelerations are larger to the northwest we eliminate the site effects by forming the ratio of the peak accelerations recorded at the same station after correcting for the epicentral distance this analysis indicates that source directivity caused a total variation of a factor of 10 in the peak accelerations comparison of this variation with the spatial extent of the aftershock sequences suggests that the strong directivity in the radiated accelerations is the result of unilateral ruptures in both events the accelerograms recorded at 10 stations within 35 km of the events were digitized to analyze the azimuthal variation of the rms acceleration the peak velocity and the radiated energy flux the variation of rms acceleration correlates almost exactly with the variation of the peak accelerations this correlation is analyzed using both deterministic and stochastic models for the acceleration waveforms the peak velocities corrected for epicentral distance vary with azimuth by a factor of 5 for both events while the radiated energy flux varies by a factor of 30 for the main shock and 15 for the aftershock the peak velocities are strongly correlated with the radiated energy flux the radiated seismic energies are estimated to be 2 6 0 9 1020 dyne cm for the main shock and 1 5 0 3 1020 dyne cm for the aftershock

zircon u pb and hf isotope constraints on the mesozoic tectonics and crustal evolution of southern tibet the first in situ hf and u pb isotope analyses of zircon separates from mesozoic granites in southern tibet identify a significant previously unknown stage of magmatism igneous zircons n 34 from a granite within the gangdese batholith show a weighted mean 206 pb 238 u age of 188 1 1 4 ma and e hf t the parts in 10 4 deviation of initial hf isotope ratios between the zircon sample and the chondritic reservoir values between 10 4 and 16 8 suggesting predominantly early jurassic intrusive activity with a juvenile mantle contribution of 40 inherited zircons from two cretaceous s type granites in the northern magmatic belt 23 delineate a slightly older 206 pb 238 u age cluster between 188 and 210 ma these zircons have e hf t values from 3 9 to 13 7 yielding crustal hf model ages from ca 1 4 to 2 1 ga suggesting a major episode of crustal growth in proterozoic time and remelting of this crust in early jurassic time combining these with literature data we interpret the jurassic gangdese magmatism as an early product of the neo tethyan subduction that played a long lasting role in the tectonic evolution of southern tibet prior to the india asia collision

the gangdese magmatic constraints on a latest cretaceous lithospheric delamination of the lhasa terrane southern tibet abstract this paper presents a systematic study of late cretaceous to early tertiary magmatism in the eastern gangdese batholith lhasa nyingchi segment the southern lhasa terrane southern tibet based on zircon u pb age rock association and geochemistry the magmatism can be divided into three stages 95 86 ma stage 1 85 73 ma stage 2 and 68 60 ma stage 3 with a short magmatic quiescence or gap of 73 69 ma the stage 1 is composed mainly of intermediate mafic rocks including norite gabbro diorite charnockite and subordinate granodiorite whereas the stage 2 consists largely of felsic rocks including tonalite granodiorite monzogranite and granite in contrast the stage 3 comprises a wide range of rock types from gabbro to granite all these samples show i type geochemical characteristics with high zircon hf isotope ratios e hf t 7 1 13 0 medium k calc alkaline and metaluminous to weakly peraluminous features stage 1 granitic rocks show low sr y 31 43 and la yb n 6 3 8 3 ratios in contrast to stage 2 that show gradually increasing sr y 47 450 and la yb n 6 38 ratios their magma temperatures constrained by titanium in zircon thermometer decrease from 740 719 c stage 1 to 685 628 c stage 2 these are consistent with an interpretation of progressive magma source deepening coupled with crustal thickening stage 3 granitoids in turn show reduced sr y 3 40 and la yb n 3 5 17 6 ratios and elevated magma temperatures 755 676 c suggesting a high temperature low pressure melting condition and thus a shallow magma source which may imply the presence of thinner crust the contrasting temperature pressure conditions between stage 2 and stage 3 magmatism suggest an event of crustal thinning which we interpret through a synthesis of other lines of geologic evidence to have resulted from a lithospheric delamination of the lhasa terrane in the latest cretaceous in conclusion the lhasa terrane underwent continuous crustal thickening in the late cretaceous and lithospheric delamination in the latest cretaceous this implies that the overthickened crust of southern tibet is generated in the cenozoic

teleseismic estimates of radiated seismic energy the e m 0 discriminant for tsunami earthquakes we adapt the formalism of boatwright and choy for the computation of radiated seismic energy from broadband records at teleseismic distances to the real time situation when neither the depth nor the focal geometry of the source is known accurately the analysis of a large data set of more than 500 records from 52 large recent earthquakes shows that this procedure yields values of the estimated energy ee in good agreement with values computed from available source parameters for example as published by the national earthquake information center neic the average logarithmic residual being only 0 26 units we analyze the energy to moment ratio by defining log10 ee m0 for regular earthquakes this parameter agrees well with values expected from theoretical models and from the worldwide neic catalogue there is a one to one correspondence between values of that are deficient by one full unit or more and the so called tsunami earthquakes previously identified in the literature as having exceedingly slow sources and believed due to the presence of sedimentary structures in the fault zone our formalism can be applied to single station measurements and its coupling to automated real time measurements of the seismic moment using the mantle magnitude mm should significantly improve real time tsunami warning

determination of earthquake energy release and ml using terrascope we estimated the energy radiated by earthquakes in southern california using on scale very broadband recordings from terrascope the method we used involves time integration of the squared ground motion velocity and empirical determination of the distance attenuation function and the station corrections the time integral is typically taken over a duration of 2 min after the p wave arrival the attenuation curve for the energy integral we obtained is given by q r cr n exp kr r 2 2 h ref 2 with c 0 49710 n 1 0322 k 0 0035 km 1 and h ref 8 km where is the epicentral distance a similar method was used to determine m l using terrascope data the station corrections for m l are determined such that the m l values determined from terrascope agree with those from the traditional optical wood anderson seismographs for 1 5 6 5 m l saturates the ratio e s m 0 m 0 seismic moment a measure of the average stress drop for six earthquakes the 1989 montebello earthquake m l 4 6 the 1989 pasadena earthquake m l 4 9 the 1990 upland earthquake m l 5 2 the 1991 sierra madre earthquake m l 5 8 the 1992 joshua tree earthquake m l 6 1 and the 1992 landers earthquake m w 7 3 are about 10 times larger than those of the others that include the aftershocks of the 1987 whittier narrows earthquake the sierra madre earthquake the joshua tree earthquake and the two earthquakes on the san jacinto fault the difference in the stress drop between the mainshock and their large aftershocks may be similar to that between earthquakes on a fault with long and short repeat times the aftershocks which occurred on the fault plane where the mainshock slippage occurred had a very short time to heal hence a low stress drop the repeat time of the major earthquakes on the frontal fault systems in the transverse ranges in southern california is believed to be very long a few thousand years hence the events in the transverse ranges may have higher stress drops than those of the events occurring on faults with shorter repeat times such as the san andreas fault and the san jacinto fault the observation that very high stress drop events occur in the transverse ranges and the los angeles basin has important implications for the regional seismic potential the occurrence of these high stress drop events near the bottom of the seismogenic zone strongly suggests that these fault systems are capable of supporting high stress that will eventually be released in major seismic events characterization of earthquakes in terms of the e s m 0 ratio using broadband data will help delineate the spatial distribution of seismogenic stresses in the los angeles basin and the transverse ranges

seismic energy for shallow earthquakes in southwest japan we estimated the radiated energy for 88 earthquakes with m w 4 0 that occurred in southwest japan by analyzing direct s waves recorded by the freesia broadband seismic network f net at local regional distances the attenuation factors for each source receiver group are estimated by the coda method using the recently established high sensitivity and high density seismic network hi net data the site amplification for each f net station is calculated by using the coda normalization method we find good agreement between the source spectra calculated by using these propagation effects and those determined by using an empirical green s function method for 4 m w 5 earthquakes the station variance of the energy estimates is a factor of 2 for earthquakes with m w 5 and it is a factor of 3 for the smaller ones the resultant energy estimations show that 1 the overall ratios between the radiated energy and seismic moment e s m 0 lie in a band of 5 10 7 to 5 10 5 for 88 earthquakes with m w 3 9 6 8 for the tottori 2000 earthquake and its 16 aftershocks with m 0 10 14 10 19 n m we found a very slight increase of e s m 0 with increasing seismic moment for the rest of earthquakes e s m 0 are almost constant for earthquakes occurring in the same source area the constants corresponding to each source area are slightly but systematically different indicating a location dependence of e s m 0 earthquakes that occur in volcanic areas show lower e s m 0 values relative to earthquakes occurring in nonvolcanic areas however the weak to non scaledependence is not certain because of the scantiness of large earthquakes in the data set 2 for earthquakes in the volcanic source areas there appears to be a trend of increasing e s m 0 with focal depth we believe this trend is a result of at least in part the decreasing attenuation with depth in volcanic areas for the nonvolcanic area earthquakes however e s m 0 tends to be a constant in a depth range of 12 50 km since the akinada earthquake 2001 m w 6 8 took place in the western part of subducting philippine sea slab and the two deeper earthquakes occurred in the upper mantle below the same subducting slab in its eastern part the constancy of e s m 0 may indicate that the apparent stress at the brittle patches in the ductile zone can be as high as in the brittle part of the crust 3 estimated e s m 0 are inversely related to the surface heat flow approximately such a relationship suggests that the geothermal structure may be the dominating tectonic factor that controls the level of apparent stress

rupture process of an extended earthquake sequence teleseismic analysis of the chilean earthquake of march 3 1985 broadband displacement and velocity records of p waves recorded at teleseismic distances are analyzed to determine the static and dynamic source parameters of the chilean earthquake of march 3 1985 ms 7 8 and seven large mb 5 6 aftershocks besides the usual parameters of depth moment and focal mechanism the analysis provides estimates of radiated energy associated stresses and source dimension for seven of the eight modeled earthquakes to insure that the parameters are describing the same physical aspect of the rupture process the parameters for each shock are computed from the time window in which most of the seismic energy is radiated the hypocenters of an additional 149 small foreshocks and aftershocks of 4 5 mb 5 6 are computed with the method of joint hypocenter determination the spatial and temporal variation of source parameters and earthquake locations are used to infer a description of the rupture process before during and after the main shock the main shock and modeled aftershocks occurred as thrust faulting on the interface between the nazca and south american plates a cluster of foreshocks in the 10 days preceding the main shock also involved thrust faulting on the plate interface the main shock itself was a complex rupture consisting of three events the first two events denoted as ms1 and ms2 released minor amounts of energy and they occurred on the periphery of the small zone defined by the foreshocks the major release of energy occurred with the third event denoted by ms which nucleated downdip of ms1 and ms2 within the time window of major energy release rupture extended approximately 90 km to the south from the point at which ms nucleated the size of the early aftershock zone which far exceeds our inferred dimensions for the major shock and the strong frequency dependence of scalar moment imply that substantial slow slip occurred on the plate interface that was not associated with major energy release the dip of the seismically active interface increases landward from about 15 at 20 km depth to about 35 at 40 km depth where the seismogenic interface reaches its maximum observed depth the vast majority of small and moderate aftershocks occurred in the shallower part of the interface as did the foreshocks ms1 and ms2 it appears however that the asperities that control the rupture of the largest earthquakes are the asperities on the deeper interface ms and the largest 1985 aftershock that of april 9 ms 7 2 occurred on the deeper interface the chilean earthquake of july 9 1971 ms 7 5 which occurred just north of the 1985 earthquake also involved main shock rupture on the deeper interface and was followed by aftershocks on the shallower interface asperities on the deeper interface appear on average to be stronger larger and more uniform in size than asperities on the shallow interface a greater average strength for deeper asperities is suggested by the observation that stress drops and apparent stresses of the main shock and the after shock of april 9 ms 7 2 are higher than those of the large shallow aftershocks and is consistent with arguments that the maximum friction on a fault zone should be near the base of the seismogenic zone one deep interface aftershock with low stress drop and apparent stress may have been situated below the point of maximum friction on the interface the inferred larger and more uniform sizes of the deeper asperities are suggested by the higher proportion of large to small shocks on the deeper interface the downdip change in the asperity size distribution may reflect increased friction on the plate interace or the coalescence of small asperities into large asperities with increased subduction the radiated energies of the main shock and large aftershocks that are computed directly from the broad bandwidth data are smaller than those predicted by standard empirical formulas by factors ranging from 4 to 20

apparent stress fault maturity and seismic hazard for normal fault earthquakes at subduction zones summary the behavior of apparent stress for normal fault earthquakes at subduction zones is derived by examining the apparent stress a es m 0 w here es is radiated energy and m 0 is seismic moment of all globally distributed shallow depth h 1m pa are also generally intraslab but occur where the lithosphere has just begun subduction beneath the overriding plate they usually occur in cold slabs near trenches where the direction of plate motion across the trench is oblique to the trench axis or where there are local contortions or geometrical complexities of the plate boundary lower a 1m pa i sassociated with events occurring at the outer rise or complex between the or and the trench axis as well as with intracrustal events occurring just landward of the trench the average apparent stress of intraslab normal fault earthquakes is considerably higher than the average apparent stress of interplate thrust fault earthquakes in turn the average a of strike slip earthquakes in intraoceanic environments is considerably higher than that of intraslab normal fault earthquakes the variation of average a with focal mechanism and tectonic regime suggests that the level of a is related to fault maturity lower stress drops are needed to rupture mature faults such as those found at plate interfaces that have been smoothed by large cumulative displacements from hundreds to thousands of kilometres in contrast immature faults such as those on which intraslab normal fault earthquakes generally occur are found in cold and intact lithosphere in which total fault displacement has been much less from hundreds of metres to a few kilometres also faults on which high a oceanic strike slip earthquakes occur are predominantly intraplate or at evolving ends of transforms at subduction zones earthquakes occurring on immature faults are likely to be more hazardous as they tend to generate higher amounts of radiated energy per unit of moment than earthquakes occurring on mature faults we have identified earthquake pairs in which an interplate thrust and an intraslab normal earthquake occurred remarkably close in space and time the intraslab normal member of each pair radiated anomalously high amounts of energy compared to its thrust fault counterpart these intraslab earthquakes probably ruptured intact slab mantle and are dramatic examples in which me an energy magnitude is shown to be a far better estimate of the potential for earthquake damage than mw this discovery may help explain why loss of life as a result of intraslab earthquakes was greater in the 20th century in latin america than the fatalities associated with interplate thrust events that represented much higher total moment release

global evaluation of large earthquake energy from 1997 through mid 2010 1 we develop an updated radiated seismic energy e catalog of global earthquakes with seismic moment m0 1019 nm mw 6 7 from 1997 through mid 2010 recording 342 events using 17849 seismograms station specific corrections and event duration dependent total p wave group calculations allow for improved e determinations for large and long duration earthquakes we find the global mean energy to moment ratio log10 e m0 4 59 0 36 robust deviations are found for thrust t 4 74 strike slip ss 4 44 and normal n 4 51 faulting events for two regions with recent energy deficient tsunami earthquakes tse we examine all large thrust earthquakes with m0 5 1017 nm mw 5 7 to identify a regional characterization in their relative radiated energy to moment ratio while thrust events along java indonesia the site of two tse events in 1994 and 2006 are comparable t jv 4 91 to the global thrust average events along the middle america trench mat the site of the 1992 nicaragua tse are consistently deficient t mat 5 15 along the mat a trend of increased deficiency in radiated energy to the southeast occurs in the direction of more rapid plate convergence global thrust mechanisms become increasingly deficient in radiated seismic energy at shallow depths with tse events representing the shallowest and most energetically deficient end member of the events in the study results suggest thrust events are highly variable but tend to increase in apparent stress from about 15 kpa near the surface to 2 mpa near 70 km depth

seismic energy release in mexican subduction zone earthquakes we estimated radiated seismic energy es from 18 shallow thrust mexican subduction zone earthquakes 4 1022 m 1 1 1028 dyne cm 11 h 37 km using digital accelerograms from the guerrero accelerograph array and neglecting stations with large site effects es is computed by integrating squared velocity spectra after applying geometrical spreading and q corrections we discarded epicentral recordings for the largest michoacan event we find that log es m 4 152 0 275 which gives a median value of es m and apparent stress a of 7 1 10 5 and 24 bars respectively the median es m value is in accordance with gutenberg and richter s g r 1956 formula for es in which es m 5 10 5 is implicit worldwide es m data where es is computed from local records mostly fall between 5 10 5 and 5 10 4 for events with m 1022 dyne cm on the other hand es m values generally lie between 5 10 6 and 5 10 5 if es is estimated from teleseismic records especially anomalous are the es m data from kikuchi and fukao 1988 for large and great earthquakes which fall near 5 10 6 thus while es from local data suggests that the g r relation seldom overestimates seismic energy release the teleseismic data point to the contrary the cause of this discrepancy may lie in the difficulty of resolving incoherent radiation from the fault and inappropriate choice of t in the analysis of teleseismic data

tectonic activities in dongshangen polymetallic ore district eastern kunlun mountains qinghai tibet plateau evidences from fission track thermochronology abstract the dongshangen polymetallic ore district is located in east section of the eastern kunlun mountains the south eastern margin of the qaidam basin and is one of the most important polymetallic mining regions however there is a relatively lower research degree on the district especially lack of thermochronology data using zircon and apatite fission track thermochronology this article reveals the tectonic activities experienced in dongshangen polymetallic ore district it is shown that the 6 age groups of geological thermal events i e 197 ma 182 170 ma 160 110 ma 110 80 ma 80 23 ma and 23 0 ma combining with the geological characteristics we consider that the 6 age groups belong to 6 periods of tectonic activities in the dongshangen ore district the first age group of 197 ma mainly resulted from the south kunlun fault strike slip and the qiangtang block collided with the kunlun block the second age group of 182 170 ma was related to the songpan ganzi bayan har basin closure the third age group of 160 110 ma was bound up with the left lateral strike slip took place in the north kunlun fault and the squeezing and collaging activities of both the qiangtang block and the lhasa block toward the eurasian plate the event during the 110 80 ma responded to northward subduction of the neo tethys ocean and strong strike slip in the altyn tagh fault zone the late cretaceous to miocene 80 23 ma was affected by the northward napping of the central kunlun terrane and the northward subduction of neo tethys ocean continued from 23 ma to present the collision of the indian plate toward the eurasian plate occurred leading to the qinghai tibet plateau uplifting as a whole synthesizing the geological characteristics and the tectonothermal events we discriminate 4 metallogenic periods of 197 ma 182 170 ma 155 112 ma and 109 88 ma in the district better reflecting the characteristics of multi stage metallogenic events

zircon u pb geochronology and hf isotopic constraints on petrogenesis of the gangdese batholith southern tibet abstract during the mesozoic cenozoic northward neotethyan subduction and subsequent india asia collision gave rise to the extensive transhimalayan magmatism that stretches from burma and western yunnan through southern tibet to the ladakh and kohistan complexes to understand the age distribution and petrogenesis of the gangdese batholith the largest intrusive exposure along the transhimalayan magmatic belt fifty granitic samples were selected for in situ zircon u pb and hf isotopic analyses the u pb data suggest four discrete stages of magmatic activity i e 205 152 109 80 65 41 and 33 13 ma respectively with the 65 41 ma stage being the most prominent the hf isotopic data indicate that the gangdese batholith is overwhelmed by positive e hf t values which are comparable to those of the kohistan ladakh batholiths in the west but differ markedly from those of the chayu burma batholiths in the east most of the gangdese granites show similar and young hf model ages 1000 300 ma indicating their derivation from juvenile crust however those formed in the 65 41 ma stage exhibit more heterogeneous hf isotopic ratios with negative e hf t values being observed in some granites younger than 50 ma suggesting the involvement of old indian continental crust in their petrogenesis this age may thus mark the onset of the india asia collision the newly established zircon u pb age and hf isotope database of the gangdese batholith can be used as a powerful tracer or fingerprint when studying the source to sink relation of the sediments eroded from the southern tibetan plateau

gravity waves generated by convection in the darwin area during the darwin area wave experiment 1 we describe a model study of gravity wave generation by convection near darwin during the darwin area wave experiment dawex we focus on an 7 hour period on 17 november 2001 for our study which included a major hector event over the tiwi islands followed by a continental convective outbreak that included a squall line to the southeast of darwin our model is a dry version of a three dimensional 3 d mesoscale cloud resolving model with horizontally uniform background wind and stability fields the model is forced with a spatially and temporally varying heating field representative of the convective latent heating in the area we derive this heating field from 3 d volumetric reflectivity from the precipitation radar located just north of darwin at gunn point during dawex the conversion from radar reflectivity to a 3 d latent heating field requires numerous assumptions that do not allow a quantitative heating estimate but which do provide a very realistic measure of the spatial morphology and temporal variations of the latent heating gravity waves generated by convective heating are known to be very sensitive to these characteristics of the heating so our resulting wave field is likely to be the most realistic possible description of the waves emitted from convection in the darwin area uncertainty in the wave amplitudes is large both because the input forcing heating is uncertain and because the use of heating as the sole proxy for wave forcing leads to additional wave amplitude errors we therefore compare the results of the model to other dawex measurements and previous modeling studies both to validate the waves in our model and to calibrate the wave amplitudes our study in turn aids in the interpretation of the other gravity wave measurements during dawex

global patterns of radiated seismic energy and apparent stress radiated energies from shallow earthquakes with magnitudes 5 8 that occurred between 1986 and 1991 are used to examine global patterns of energy release and apparent stress in contrast to traditional methods which have relied upon empirical formulas these energies are computed through direct spectral analysis of broadband seismic waveforms velocity squared spectra of body waves are integrated after they have been corrected for effects arising from depth phases frequency dependent attenuation and focal mechanism the least squares regression fit of energy es to surface wave magnitude ms for a global set of 397 earthquakes yields log es 4 4 1 5ms which implies that the gutenberg richter relationship overestimates the energies of earthquakes the least squares fit between es and seismic moment m0 is given by the relationship es 1 6 10 5 m0 which yields 0 47 mpa as the average global value of apparent stress however the regression lines of both es ms and es m0 yield poor empirical predictors for the actual energy radiated by any given earthquake the scatter of data is more than an order of magnitude about each of the regression lines on the other hand global variations between es and mo while large are not random when subsets of es m0 are plotted as a function of seismic region and faulting type the scatter in data is substantially reduced the es m0 fits for many seismic regions and tectonic environments are very distinctive and a characteristic apparent stress c can be derived the lowest apparent stresses 3 0 mpa are associated with strike slip earthquakes that occur at oceanic ridge ridge transforms and in intraplate environments seaward of island arcs intermediate values of apparent stress 1 5 a 3 0 mpa are associated with strike slip earthquakes at incipient or transitional plate boundaries in general the dominant mode of failure for a tectonic environment is associated with the faulting type that has the lowest apparent stress an energy magnitude me can complement moment magnitude mw in describing the size of an earthquake me being derived from velocity power spectra is a measure of seismic potential for damage mw being derived from the low frequency asymptote of displacement spectra is more physically related to the final static displacement of an earthquake when earthquake size is ranked by moment a list of the largest events is dominated by earthquakes with thrust mechanisms when earthquake size is ranked by energy the list of the largest events is dominated by strike slip earthquakes

the lhasa terrane record of a microcontinent and its histories of drift and growth article i nfo the lhasa terrane in southern tibet has long been accepted as the last geological block accreted to eurasia before its collision with the northward drifting indian continent in the cenozoic but its lithospheric architecture drift and growth histories and the nature of its northern suture with eurasia via the qiangtang terrane remain enigmatic using zircon in situ u pb and lu hf isotopic and bulk rock geochemical data of mesozoic early tertiary magmatic rocks sampled along four north south traverses across the lhasa terrane we show that the lhasa terrane has ancient basement rocks of proterozoic and archean ages up to 2870 ma in its centre with younger and juvenile crust phanerozoic accreted towards its both northern and southern edges this finding proves that the central lhasa subterrane was once a microcontinent this continent has survived from its long journey across the paleo tethyan ocean basins and has grown at the edges through magmatism resulting from oceanic lithosphere subduction towards beneath it during its journey and subsequent collisions with the qiangtang terrane to the north and with the indian continent to the south zircon hf isotope data indicate significant mantle source contributions to the generation of these granitoid rocks e g 50 90 0 70 and 30 100 to the mesozoic magmatism in the southern central and northern lhasa subterranes respectively we suggest that much of the mesozoic magmatism in the lhasa terrane may be associated with the southward bangong nujiang tethyan seafloor subduction beneath the lhasa terrane which likely began in the middle permian or earlier and ceased in the late early cretaceous and that the significant changes of zircon ehf t at 113 and 52 ma record tectonomagmatic activities as a result of slab break off and related mantle melting events following the qiangtang lhasa amalgamation and india lhasa amalgamation respectively these results manifest the efficacy of zircons as a chronometer u pb dating and a geochemical tracer hf isotopes in understanding the origin and histories of lithospheric plates and in revealing the tectonic evolution of old orogenies in the context of plate tectonics

evolution of the northward subduction of the neo tethys implications of geochemistry of cretaceous arc volcanics in qinghai tibetan plateau abstract the subduction history of the neo tethys is critical for understanding the early stage tectonic evolution of the qinghai tibetan plateau arc volcanics which usually record the subduction process thus can help trace the evolution of the neo tethys in this paper we report for the first time zircon la icp ms u pb ages whole rock major and trace element geochemistry as well as sr nd hf isotopic data for the volcanic rocks from the danshiting formation at the southern edge of the lhasa terrane two samples collected from the danshiting formation volcanic rocks yield zircon u pb ages of 97 1 0 8 ma and 90 1 0 7 ma respectively the volcanic rocks of the danshiting formation are composed dominantly of basalts and andesites characterized by depletion of high field strength elements low 87sr 86sr i ratios 0 7042 0 7045 high end t 2 1 to 2 7 and high ehf t values mean 11 6 these rocks show geochemical affinity of arc volcanic rocks the volcanic rocks of the danshiting formation may have formed during the northward subduction of the neo tethys this paper combined with the previous data briefly divide the northward subduction of the neo tethys into four periods the jurassic 195 150 ma when the magmatic activity was controlled by the northward high angle subduction of the neo tethys that produced volcanics of the yeba and bima formations the late jurassic to early cretaceous 150 110 ma when the neo tethys subducted at low angles the late cretaceous 110 80 ma when subduction was of the andean type and produced volcanic rocks of the danshiting formation the paleogene 70 40 ma when there was transition from ocean to continent subduction to continent continent collision that produced volcanic rocks of the linzizong group

spectral analysis of k net and kik net data in japan part ii on attenuation characteristics source spectra and site response of borehole and surface stations abstract in this study we apply a nonparametric spectral inversion scheme to a data set of accelerograms recorded by the k net and kik net networks in japan in order to derive attenuation characteristics source spectra and site response for this purpose we use a total of more than 67 000 s wave records from 2178 earthquakes m jma 2 7 8 obtained at 1555 stations at the earth s surface and more than 29 000 records from 1826 events recorded at 637 borehole stations at depths of 100 to 3000 m attenuation characteristics are investigated in five separate regions showing that crustal q depicts lower values in central compared to southern japan and a significant frequency dependence is observed in every region the source spectra follow the 2 model with higher stress drops for subcrustal earthquakes as compared with crustal ones while strong amplification effects dominate the site contributions for the surface sensors those for the borehole sensors are characterized by smaller variability nevertheless consistent with observations from deconvolution of borehole surface recording pairs downgoing wave effects are visible in the site contributions for many borehole stations finally the site amplification functions obtained at the surface are compared with surface to borehole s b and horizontal to vertical h v spectral ratios showing that the s b ratios generally provide better estimates of the horizontal amplification than the h v ratios due to amplification of the vertical component of ground motion

petrogenesis of the early eocene adakitic rocks in the napuri area southern lhasa partial melting of thickened lower crust during slab break off and implications for crustal thickening in southern tibet abstract cenozoic adakitic rocks in the lhasa block southern tibet have been widely used to trace the lateral extent of crustal thickening however their petrogenesis remains controversial here we report geochronological and geochemical data for the napuri intrusive rocks in the core area of the quxu batholith southern lhasa zircon u pb dating suggests that they were generated at approximately 48 ma the studied samples show significant geochemical variations manifested by the coexistence of three types of igneous rocks groups i and ii rocks exhibit variable and high sio 2 66 4 73 9 wt high al 2 o 3 14 0 17 4 wt k 2 o 3 9 5 3 wt sr 273 718 ppm and sr y 18 3 to 81 3 values and low y 3 6 to 16 ppm heavy rare earth element ree e g yb 0 48 to 1 8 ppm mgo 0 4 1 0 wt cr 2 9 7 4 ppm and ni 1 6 4 5 ppm contents which are similar to those of thickened lower crust derived adakitic rocks the group i rocks show higher sr y 77 5 81 3 ratios and lower total ree 55 5 63 2 ppm contents with clearly positive eu and sr anomalies whereas the group ii rocks have relatively lower sr y 18 3 65 7 ratios and higher total ree 115 375 ppm contents with negligible or slightly negative eu and sr anomalies group iii rocks have the highest sio 2 74 5 76 0 wt y 17 0 23 7 ppm and yb 2 91 3 30 ppm contents and the lowest al 2 o 3 12 5 13 2 wt sr 81 3 141 ppm and sr y 4 8 5 9 values with distinctly negative eu and sr anomalies compared with the jurassic cretaceous granitoids in southern lhasa the relative enrichment in sr nd hf isotopic compositions 87 sr 86 sr i 0 7049 0 7055 e nd t 0 3 0 7 and e hf t zircon 3 6 11 4 for the napuri intrusive rocks indicates that they likely contained indian continental components the group i and group ii rocks most probably originated from thickened mafic lower crust amphibolite eclogites or garnet amphibolites with garnet rutile plagioclase as residual minerals in the source at 1 5 gpa corresponding to depths of 50 km and group iii rocks were probably generated by fractional crystallization of plagioclase from the adakitic magmas taking into account the narrow linear nature of the eocene magmatic belt and reported synchronous asthenosphere derived basaltic rocks in southern lhasa we suggest that upwelling asthenosphere triggered by the break off of subducted neo tethyan slab probably provided the required thermal conditions for lower crustal melting the identification of indian continental components in the napuri intrusive rocks probably indicates that the asia india collision had taken place prior to their emplacement the dramatic changes in the la yb n ratios and e nd t and e hf t values of magmatic rocks in the gangdese area at ca 51 46 ma indicate that the cenozoic crustal thickening associated with the indentation of the indian continent began in the early eocene ca 51 46 ma at the latest

broadband source spectrum seismic energy and stress drop of the 1989 macquarie ridge earthquake i compute the broadband source spectrum at periods from 1 to 50 seconds using teleseismic p body waves of the may 23 1989 macquarie ridge earthquake mw 8 1 recorded by the gdsn geoscope and ida networks the average source spectrum is obtained by windowing tapering and fourier transforming p waves removing from the spectra the effects of attenuation geometrical spreading and radiation pattern and averaging logarithmically over the stations the source spectrum for the strike slip macquarie ridge earthquake is higher than an average source spectrum of seven recent large earthquakes scaled to be comparable to a mw 8 1 earthquake by a factor of 2 to 3 at periods of 1 to 20 seconds these other earthquakes were underthrusting events in subduction zones using haskell s formulation assuming a point source with no directivity i estimate the seismically radiated energy from the source spectrum by integrating the square of the source spectrum in velocity and scaling the result the seismic energy thus estimated for the macquarie ridge earthquake is 3 to 8 1023 ergs an orowan stress drop can be obtained from the seismic energy and moment the orowan stress drop for the macquarie ridge earthquake is about 20 to 50 bars much higher than similarly determined stress drops of other recent large earthquakes there is a correlation between the orowan stress drops and time since the last earthquake of comparable or larger magnitude for seven recent large earthquakes this correlation suggests that a healing process operates that may control the mechanical strength of the fault and is important on time scales of tens to hundreds of years

metagabbros of the gangdese arc root south tibet implications for the growth of continental crust abstract the gangdese batholith of the southern tibet is regarded as a principal component of the pre himalayan collision andean type continental margin however the gangdese arc root remains poorly understood because most of previous studies concentrated on granitoids from the middle and upper crust of the arc in this paper we report new petrochemical and geochronological characteristics of metagabbros occurred as a part of the late cretaceous lilong batholith in the eastern gangdese arc the gabbros with protolith ages of 82 95 ma underwent slightly late 68 77 ma granulite facies metamorphism and nearly coeval anatexis typical of lower arc crusts the metagabbros and diorites of the lilong batholith show geochemical features of juvenile crust in continental magmatic arcs and an original magmatic differentiation trend with gabbros formed as the cumulates whereas diorites as the crystalline products of evolved magmas the lilong batholith forms an arc crustal section with a possible thickness of 20 30 km and may have an overall andesitic composition providing a support for the andesite model of continental crust growth

vp and vs tomography of kyushu japan new insight into arc magmatism and forearc seismotectonics high resolution 3d tomographic images of the crust and upper mantle are determined under the entire arc region of kyushu by inverting a large number of arrival time data of p wave s wave and sp depth phases the hypocenters of the offshore earthquakes are accurately located by applying a combined method that makes use of sp depth phases and double difference of the p wave and s wave arrival times our tomographic results confirmed the major features delineated by previous studies and revealed some new features of structural heterogeneity related to arc magmatism beneath the kyushu arc high velocity anomalies correspond to the subducting philippine sea slab whilst low velocity zones in the mantle wedge are clearly imaged coinciding well with the velocity heterogeneities determined by previous studies strong low velocity anomalies are distributed extensively along the volcanic front and extend to the back arc side in the crust and upper mantle our present seismic images together with the previous results suggest that the fluids supplied by the dehydration processes of the descending philippine sea slab and the hot upwelling materials related to the opening of the okinawa trough might have contributed to the unzen volcanism under the forearc region of the kyushu subduction zone prominent low velocity anomalies are visible at depths of 10 40 km along the upper boundary of the subducting philippine sea slab and large interplate earthquakes are located around the low velocity zones which is similar to the features revealed in the northeast japan subduction zone the spatial distribution of large thrust earthquakes below the offshore region reflects the existence of weak coupled or decoupled portions on the interface of the subducting philippine sea slab

seismic structure and magmatism of the young kyushu subduction zone seismic tomography shows that slow velocity anomalies exist in the mantle wedge and extend to the forearc region down to the subducting philippine sea slab under northern kyushu we also conducted numerical simulations with petrologic data to estimate the fluid distribution in the mantle wedge the seismic and simulation results indicate that regimes of melting and magmatism in a subduction zone with a young slab are different from those with old slabs dehydration and melting occur beneath the arc and forearc above a young and warm slab such as in northern kyushu 26 ma while fluids aqueous solution and melt occur mainly beneath the back arc above old slabs

early cretaceous subduction related adakite like rocks of the gangdese belt southern tibet products of slab melting and subsequent melt peridotite interaction abstract the limited geochronology and geochemistry data available for the early cretaceous igneous rocks of the southern gangdese belt southern tibet has resulted in the proposal of conflicting geodynamic models for the generation of the widespread cretaceous igneous rocks in the middle and northern parts of the belt to explore this issue we present shrimp u pb zircon data and geochemical and sr nd pb hf isotopic data for the mamen andesites from the southern margin of the gangdese belt the mamen andesites emplaced at 136 5 ma are sodic na2o k2o 1 2 2 3 and have geochemical characteristics typical of adakites i e high al2o3 high la yb ratios and sr contents low y and hree contents and positive eu anomalies except for high cr ni and mgo contents the andesites have initial 87sr 86sr t ratios of 0 70413 0 70513 positive end t values of 3 7 5 8 and 206pb 204pb t ratios of 18 37 18 51 207pb 204pb t ratios of 15 59 15 65 and 208pb 204pb t ratios of 38 43 38 72 in situ hf isotopic analyses of zircons that had previously been dated by shrimp yielded positive initial ehf t values ranging from 11 0 to 15 5 a model calculation using trace element and sr nd pb isotopic data indicates that several percent of subducted sediment is required to generate the mamen andesites which were derived via the partial melting of subducted neo tethyan slab morb sediment fluid and subsequently hybridized by peridotite in the mantle wedge our data indicate that the neo tethyan oceanic crust was subducted northward beneath the gangdese belt during the early cretaceous at a high angle our results are inconsistent with a tectonic model that advocates the low angle or flat slab subduction of neo tethyan oceanic crust in generating the widespread cretaceous magmatism recorded in the gangdese belt

early late cretaceous ca 93 ma norites and hornblendites in the milin area eastern gangdese lithosphere asthenosphere interaction during slab roll back and an insight into early late cretaceous ca 100 80 ma magmatic flare up in southern lhasa tibet abstract at more than 500 km in length the mainly jurassic early eocene gangdese batholith is one of the most important constituents of the southern lhasa sub block and provides an ideal site for study of tibetan orogenesis recent studies on gangdese intermediate felsic intrusive rocks mainly granites demonstrate that remarkable crustal growth as well as an early late cretaceous ca 100 80 ma magmatic flare up event occurred in southern tibet however the mechanism that drove this magmatic event and its relationship to crustal growth event are not yet clear here we report detailed petrological geochronological geochemical and sr nd hf isotopic data for recently identified norites and hypersthene bearing hornblendites in the milin area southern lhasa sub block these mafic rocks are dated at early late cretaceous ca 93 ma and are characterized by relatively uniform sr nd hf isotopic compositions 87sr 86sr i 0 7042 to 0 7047 end t 2 9 to 3 6 and ehf t zircon 10 9 to 17 0 suggesting that they evolved from similar parental magmas the milin norites and hornblendites are likely to be the products of mineral fractionation and accumulation from a common parental magma during the early and late stages of the magma evolution respectively thermometric calculations indicate that pyroxenes from the milin norites have high crystallization temperatures 1240 1349 c the parental magmatic compositions calculated from pyroxene trace element compositions in the milin norites show slightly flat to enriched light rare earth element lree patterns la yb n 2 9 3 4 la gd n 2 1 2 9 with variable negative nb anomalies nb la n 0 18 0 81 indicating their dual or hybrid geochemical characteristics we suggest that their parental magmas may have been generated by the interaction of upwelling asthenospheric and metasomatized lithospheric mantle taking into account the spatial and temporal distribution of the mesozoic magmatic rocks and regional paleomagnetic data we further suggest that the early late cretaceous magmatic flare up in the southern lhasa sub block was also triggered by the asthenospheric upwelling which resulted from roll back of subducted tethyan oceanic slab

cenozoic volcanism in tibet evidence for a transition from oceanic to continental subduction geochronological k ar or ar ar major and trace element sr nd pb isotopic and mineral chemical data are presented for newly discovered cenozoic volcanic rocks in the western qiangtang and central lhasa terranes of tibet alkali basalts of 65 45ma occur in the western qiangtang terrane and represent primitive mantle melts as indicated by high mgnumbers 100 mg mg fe 54 65 cr 204 839 ppm and ni 94 218 ppm contents and relatively low ratios of sr sr 0 7046 0 7061 pb pb 18 21 18 89 pb pb 15 49 15 61 and pb pb 38 42 38 89 and high ratios of nd nd 0 5124 0 5127 in contrast younger volcanic rocks in the western qiangtang terrane 30ma and the central lhasa terrane 23 13 and 8ma are potassic to ultrapotassic and interpreted to have been derived from an enriched mantle source they are characterized by very high contents of incompatible trace elements negative ta nb and ti anomalies and radiogenic pb isotopic compositions pb pb 18 43 19 10 pb pb 15 64 15 83 pb pb 39 14 39 67 sr sr 0 7088 0 7092 and nd nd 0 5122 ratios of the western qiangtang terrane potassic lavas are similar to those of 45 29ma potassic volcanic rocks in the north central qiangtang terrane whereas sr sr 0 7167 0 7243 and nd nd 0 5119 ratios of central lhasa terrane lavas are similar to those of 25 16ma ultrapotassic volcanic rocks in the western lhasa terrane the 65 45ma alkali basalts in the western qiangtang terrane along with widespread calc alkaline volcanic rocks of this age in the lhasa terrane may be related to roll back of a previously shallow north dipping slab of tethyan oceanic lithosphere beneath tibet subduction as opposed to convective thinning of continental lithosphere is favored to explain potassic volcanism in tibet because of its occurrence in distinct east west trending belts 45 29ma in the qiangtang terrane 25 17ma in the northern lhasa terrane 16 8ma in the southern lhasa terrane and temporal and spatial relationships with major thrust systems

characteristics of gravity waves with short vertical wavelengths observed with radiosonde and gps occultation during dawex darwin area wave experiment 1 we conducted radiosonde soundings during three intensive observation periods iop of the darwin area wave experiment dawex campaign in october to december 2001 and measured humidity temperature and wind velocity up to 30 35 km every 3 hours for 40 times at three sites in each iop we analyzed height profiles of kinetic ek and potential ep energy per unit mass caused by gravity waves with vertical wavelengths less than 3 km the wave energy was clearly enhanced between 15 20 km and 25 30 km and it was considerably depressed at 20 25 km between the two height regions this feature was seen at all the sites during the three iops different types of wave activity seem to contribute to the enhancement of ek and ep below 20 km and at 25 30 km height distribution of the wave energy at 20 30 km seems to correlate with the structure of the mean zonal wind shear we also analyzed a latitude height section of ep in october to december 2001 using gps occultation data with the champ satellite collected around the dawex sites 80 180 e and 30 n to 30 s compared with the latitude distribution of ep with gps the dawex results were slightly smaller and larger at 20 25 km and 25 30 km respectively the longitude variation of ep is also analyzed from the gps data at 10 15 s which generally agreed well with the dawex results at both 20 25 km and 25 30 km

zircon u pb dating and in situ hf isotopic analysis of permian peraluminous granite in the lhasa terrane southern tibet implications for permian collisional orogeny and paleogeography article i nfo the lhasa terrane has long been interpreted as a simple tectonic block rifted from gondwana during the late paleozoic and then drifted northward before finally amalgamating with the qiangtang terrane during the early cretaceous in this paper we document permian peraluminous granites near pikang in the southern margin of the central lhasa terrane close to the recently documented songdo eclogite of permian age zircon shrimp and la icpms u pb dating for a pikang granite sample gives an identical crystallization age of about 263 ma and a wide age range of inherited zircons 283 2141 ma in situ hf isotopic analyses for 20 zircons of 263 ma yielded ehf t values of 4 5 to 1 9 the pikang granites have high a cnk values 1 08 and high normative corundum 1 3 2 0 indicative of peraluminous s type granite they are characterized by moderately negative eu anomalies eu eu 0 48 0 61 and strongly negative ba nb sr p and ti anomalies the granites have high end t values 6 4 to 6 0 and low initial 87 sr 86 sr ratios 0 7082 0 7096 relative to melts derived from mature continental crust these rocks are interpreted to have been generated by mixing between mantle melts and their induced melting of mature crustal materials we interpret the pikang peraluminous granite magmatism the regional angular unconformity between the middle and the upper permian and the eclogite of the same metamorphic age 262 ma from the same geotectonic location to represent different products of a common event in time and space we tentatively term this common event as syncollisional orogeny i e the permian gangdese orogeny we further hypothesize that the orogeny may be genetically associated with the collision between the lhasa terrane and the northern margin of australia following the closure of the paleo tethyan ocean south of the lhasa terrane

mantle contributions to crustal thickening during continental collision evidence from cenozoic igneous rocks in southern tibet the tibetan crust is twice as thick as average continental crust crustal compression and shortening as a result of indian asian collision isoftenconsideredtobetheprimarycauseforthecrustalthickening inthispaper weshowthatmagmaticcontributionisalsoimportant we come to this conclusion by documenting the paleogene linzizong volcanic succession lvs its coeval granitoid batholiths and the miocene adakitic rocks along the gangdese magmatic belt in southern tibet it has been widely accepted that the indian asian collision proceeded from a soft phase at 65 70 ma to a hard phase at 45 40 ma followed by continued post collisional convergence to the present in response to the collision and post collision convergence are a series of tectono magmatic events recorded in the gangdese magmaticbelt theseinclude 1 thesyn collisionallvsvolcanism 65 40ma andtheemplacementofsoutherngangdesebatholiths a peakageof 50ma 2 aperiod 40mato25ma thatismagmaticallyquiescent yettectonicallydominatedbyactivecompressionand crustal shortening and 3 the emplacements of post collisional adakitic rocks 25 12 ma potassic ultrapotassic volcanics 25 10 ma and peraluminous muscovite bearing granites between 24 and 18 ma these three major events contribute in different ways to the crustal thickness phase i formation of the lower juvenile crust from 65 ma to 50 ma with crustal thickening largely concentrated at 50 40maviainputofmantle derivedmagmas phaseii crustalthickeningbytectonicshorteningat 40 25ma andphaseiii retaining crustal thickness but thinning of the lithospheric mantle since 25 ma in response to crustal extension and upwelling and lateral flow of asthenosphericmantle weemphasizethatcollision inducedcrustalthickeningtookplacemainlyintheperiodof 50 40maand 25ma i e the period between the late stage of the lvs volcanism and the beginning of the adakitic rock emplacement most of the lvs rocks and the collision related granitoids in southern gangdese have end 0 attesting to the significance of mantle input most likely through remelting of mantle derived basaltic rocks including the subducted neo tethyan ocean crust the petrologic and geochemical characteristics of the miocene potassic adakitic rocks support the idea that the lower portion of the thickened tibetan crust is mafic and is genetically associated with the earlier lvs magmatism we estimate that the mantle material input contributed about 30 of the total thickness of the present day tibetan crust by assuming a pre collision crustal thickness of 35 km then the tectonic contribution would be about 20 km

generation and trapping of gravity waves from convection with comparison to parameterization some parameterizations of gravity wave mean flow forcing in global circulation models gcms add realism by describing wave generation by tropospheric convection because the convection in gcms is itself a parameterized process these convectively generated wave parameterizations necessarily use many simplifying assumptions in this work the authors use a realistic simulation of wave generation by convection described in previous work which was validated by observations from the darwin area wave experiment dawex to test these assumptions and to suggest some possible improvements to the parameterizations in particular the authors find that wave trapping in the troposphere significantly modifies the spectrum of vertically propagating waves entering the stratosphere above convective wave sources and offer a linear method for computing wave transmission and reflection effects on the spectrum suitable for inclusion in the parameterizations the wave fluxes originate from both a time varying heating mechanism and an obstacle effect mechanism acting in the simulation methods for including both mechanisms in the parameterizations are described waves emanating from the obstacle effect remain very sensitive to the depth of penetration of latent heating cells into an overlying shear zone which will continue to make it difficult to accurately parameterize in a gcm where the convective cells are not resolved

darwin area wave experiment dawex field campaign to study gravity wave generation and propagation 1 a field campaign studying aspects of the atmospheric dynamics from the ground to the lower thermosphere was conducted in australia during october to december 2001 this darwin area wave experiment dawex involved participation from nine institutions in the usa japan and australia and focused on studying wave perturbations in the middle atmosphere in conjunction with detailed observations of tropospheric moist convection the experiment involved extended deployment of six airglow imagers medium frequency radars and a boundary layer wind profiler tropospheric convection was observed by meteorological radars including the australian bureau of meteorology research centre c band polarized radar this paper provides an overview of the experimental setup and also discusses the basic meteorological fields during the experiment including some examples of convective activity near darwin wind and temperature fields extending up into the lower thermosphere were constructed on the basis of campaign observations and the ukmo global meteorological analysis these fields can be used for the modeling of gravity wave propagation through the middle atmosphere

vertical velocity in oceanic convection off tropical australia abstract time series of 1 hz vertical velocity data collected during aircraft penetrations of oceanic cumulonimbus clouds over the western pacific warm pool as part of the equatorial mesoscale experiment emex are analyzed for updraft and downdraft events called cores an updraft core is defined as occurring whenever the vertical velocity exceeds 1 m s 1 for at least 500 m a downdraft core is defined analogously over 19 000 km of straight and level flight legs are used in the analysis five hundred eleven updraft cores and 253 downdraft cores are included in the dataset core properties are summarized as distributions of average and maximum vertical velocity diameter and mass flux in four altitude intervals between 0 2 and 5 8 km distributions are approximately lognormal at all levels examination of the variation of the statistics with height suggests a maximum in vertical velocity between 2 and 3 km slightly lower or equal vertical velocity is indicated at 5 km near the freezing level virtual t

can slab melting be caused by flat subduction slab melting has been suggested as a likely source of adakitic arc magmas i e andesitic and dacitic magmas strongly depleted in y and heavy rare earth elements existing numerical and petrologic models however restrict partial melting to very young 5 ma oceanic crust typically at 60 80 km depth paradoxically most of the known pliocene quaternary adakite occurrences are related to subduction of 10 45 ma lithosphere which should not be able to melt under normal subduction zone thermal gradients we propose an unusual mode of subduction known as flat subduction occurring in 10 of the world s convergent margins that can produce the temperature and pressure conditions necessary for fusion of moderately old oceanic crust of the 10 known flat subduction regions worldwide eight are linked to present or recent 6 ma occurrences of adakitic magmas observations from chile ecuador and costa rica suggest a three stage evolution 1 steep subduction produces a narrow calc alkaline arc typically 300 km from the trench above the asthenospheric wedge 2 once flat subduction begins the lower plate travels several hundred kilometers at nearly the same depth thus remaining in a pressure temperature window allowing slab melting over this broad distance and 3 once flat subduction continues for several million years the asthenospheric wedge disappears and a volcanic gap results as in modern day central chile or peru the proposed hypothesis which reconciles thermal models with geochemical observations has broad implications for the study of arc magmatism and for the thermal evolution of convergent margins

late cretaceous crustal growth in the gangdese area southern tibet petrological and sr nd hf o isotopic evidence from zhengga diorite gabbro abstract recent studies of gangdese granitic magmatism demonstrate a mantle contribution to crustal growth in southern tibet during the jurassic early eocene however the specific mechanism for adding such juvenile crust has been disputed owing to a lack of reliable evidence for contemporaneous deep mantle geodynamic processes here we report on the zhengga diorite gabbro suite from the gangdese area they consist of plagioclase labradorite and anorthite amphibole clinopyroxene biotite and minor magnetite epidote zircon and apatite la icp ms zircon u pb dating for two samples gives a ca 94 ma age for the zhengga intrusive rocks i e the late cretaceous apart from one high sio 2 52 2 wt diorite sample with slightly high k 2 o 1 72 wt and initial 87 sr 86 sr 0 7068 and low e nd t 5 6 values the gabbro samples are geochemically characterized by low sio 2 39 8 50 1 wt and k 2 o 0 3 1 1 wt strongly negative nb ta and positive sr anomalies and uniform initial 87 sr 86 sr 0 7043 0 7048 the gabbros can be divided into two groups group i gabbros with relatively low total rare earth element ree rb and sio 2 contents and positive eu anomalies and group ii gabbros with slightly higher total ree rb and sio 2 contents and negative eu anomalies the group i gabbros have e nd t 1 7 to 4 1 and zircon e hf t 6 5 to 11 1 and 18 o 5 89 to 7 24 values which are slightly different to those of the group ii gabbros 0 2 to 2 0 2 9 to 6 5 and 6 24 to 7 05 trace element compositions of amphibole and clinopyroxene grains suggest that the zhengga mafic magmas contained a significant fluid transported component probably released from subducted oceanic lithosphere we suggest that the parental magmas of the zhengga gabbros were generated by the hydrous partial melting of lithospheric mantle metasomatized by sediment melts fluids the group i gabbros were likely generated by the fractional crystallization of olivine or clinopyroxene from such parental magmas with insignificant crustal contamination whereas the group ii gabbros were probably produced by assimilation and fractional crystallization afc processes from mafic magmas that were geochemically similar to the group i gabbros pre collisional underplating of mantle derived mafic magmas likely played an important role in crustal growth and supplied the source materials for some late late cretaceous cenozoic granitoids of the gangdese batholiths this study also demonstrates that the hydrous partial melting of mantle wedge triggered by the dehydration in a subduction setting has a capacity to create significant volumes of juvenile continental crust

crustal thickening prior to 38 ma in southern tibet evidence from lower crust derived adakitic magmatism in the gangdese batholith article thepetrogenesisandgeodynamicimplicationsofthecenozoicadakitesinsoutherntibetremaintopicsofdebate here we report geochronological and geochemical data for host granites and mafi ce nclaves from wolong in the easterngangdesebatholith southerntibet zirconla icp ms datingindicatesthatthe wolonghostgranitesand enclaves were synchronously emplaced at ca 38 ma the host granites are medium to high k calc alkaline metaluminous a cnk 0 93 0 96 with high al2o3 15 47 17 68 low mgo 0 67 1 18 very low abundances of compatible elements e g cr 3 87 8 36 ppm ni 3 04 5 71 ppm and high sr y ratios 127 217 similar to those typical of adakite the mafi ce nclaves sio2 51 08 56 29 have 3 83 5 02 mgo

cretaceous tertiary shortening basin development and volcanism in central tibet the geologic map pattern of the qiangtang terrane in central tibet defines a 600 km long and up to 270 km wide east plunging structural culmination it is characterized by early mesozoic blueschist bearing melange and upper paleozoic strata in the core and mainly triassic jurassic strata along the limbs in the western qiangtang terrane 84 e the culmination is unconformably overlain by weakly deformed mid cretaceous volcanic flows and tuffs along the bangong suture to the south 32 n 84 e mid cretaceous nonmarine red beds and volcanic rocks lie unconformably on jurassic suture zone rocks these relationships demonstrate that west central tibet was above sea level during the mid cretaceous and experienced significant denudation prior to mid cretaceous time growth of the qiangtang culmination is inferred to have initiated during southward emplacement of a thrust sheet of early mesozoic melange and upper paleozoic strata during the early cretaceous lhasa qiangtang collision the north south width of the inferred thrust sheet provides a minimum slip estimate of 150 km at 84 e decreasing eastward to 70 km at 87 e paleogene deformation in the qiangtang terrane is characterized by widely distributed mainly north dipping thrust faults that cut eocene oligocene red beds and volcanic rocks in their footwall along the bangong suture the north dipping shiquanhe gaize amdo thrust system cuts 64 and 43 m y old volcanic tuffs in its footwall and accommodated 40 km of post mid cretaceous shortening the tertiary south dipping gaize siling co backthrust bounds the southern margin of the bangong suture and marks the northernmost limit of mid cretaceous marine strata in central tibet cretaceous deformation and denudation in central tibet is attributed to northward underthrusting of the lhasa terrane beneath the qiangtang terrane along the bangong suture this model implies that 1 cretaceous strata along the bangong suture and in the northern lhasa terrane were deposited in a flexural foreland basin system and derived at least in part from the qiangtang terrane and 2 the central tibetan crust was thickened substantially prior to the indo asian collision although its magnitude is poorly known tertiary shortening in the qiangtang terrane is more prevalent than in the lhasa terrane this difference may be attributed to the presence of underthrust melange in the deeper central tibetan crust which would have made it weaker than the lhasa terrane during the indo asian collision

detailed fault structure of the 2000 western tottori japan earthquake sequence we investigate the faulting process of the aftershock region of the 2000 western tottori earthquake m w 6 6 by combining aftershock hypocenters and moment tensor solutions aftershock locations were precisely determined by the double difference method using p and s phase arrival data of the japan meteorological agency unified catalog by combining the relocated hypocenters and moment tensor solutions of aftershocks by broadband waveform inversion of freesia f net we successfully resolved very detailed fault structures activated by the mainshock the estimated fault model resolves 15 individual fault segments that are consistent with both aftershock distribution and focal mechanism solutions rupture in the mainshock was principally confined to the three fault elements in the southern half of the zone which is also where the earliest aftershocks concentrate with time the northern part of the zone becomes activated which is also reflected in the postseismic deformation field from the stress tensor analysis of aftershock focal mechanisms we found a rather uniform stress field in the aftershock region although fault strikes were scattered the maximum stress direction is n107 e which is consistent with the tectonic stress field in this region in the northern part of the fault where no slip occurred during the mainshock but postseismic slip was observed the maximum stress direction of n130 e was possible as an alternative solution of stress tensor inversion

the 1978 tabas iran earthquake an interpretation of the strong motion records the 1978 tabas iran earthquake ms 7 4 was registered on 11 accelerograph stations with source distances ranging from 3 to 350 km and corresponding recorded peak accelerations ranging from 0 95 down to 0 01 g respectively at least four subevents are recognized at dayhook the station nearest the epicenter at the southern limit of the fault the rupture propagated to the northwest with a rupture velocity of 2 7 km sec the subevents occurred along the fault as the rupture proceeded curves which describe the distance dependence of peak values of acceleration velocity and band passed acceleration are constructed for this earthquake although attenuation of seismic waves is similar in iran and the western united states the prediction curves based on peak strong motion acceleration data from the western united states systematically underestimate the peak values of the tabas earthquake and the similar 1952 kern county earthquake prediction curves based on peak velocity in the western united states overestimate peak velocity for distances less than about 60 km site effects are recognized at some recording stations the various estimates for stress operating during the earthquake range from 25 to 90 bars the radiated energy of the main event estimated from accelerograms by two methods is about 2 3 1023 ergs larger than the energy calculated by gutenberg and richter energy formula with ms 7 4 each of the four subevents released approximately the same amount of energy corresponding to ms 7 teleseismic energy estimates for this earthquake are one to two orders of magnitude smaller indicating large systematic errors exist in one or both of the energy estimation procedures

comparison of global distributions of zonal mean gravity wave variance inferred from different satellite instruments gravity wave temperature fluctuations acquired by the crista instrument are compared to previous estimates of zonal mean gravity wave temperature variance inferred from the lims mls and gps met satellite instruments during northern winter careful attention is paid to the range of vertical wavelengths resolved by each instrument good agreement between crista data and previously published results from lims mls and gps met are found key latitudinal features in these variances are consistent with previous findings from ground based measurements and some simple models we conclude that all four satellite instruments provide reliable global data on zonal mean gravity wave temperature fluctuations throughout the middle atmosphere