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

extrusion tectonics of the okhotsk plate northeast asia based on focal mechanism geologic and lineament data faulting in the chersky range northeast asia changes from right lateral transpression between north america and eurasia to left lateral strike slip between north america and the okhotsk plate the location of the okhotsk eurasia boundary is poorly constrained available data suggest this boundary is right lateral the okhotsk plate is being compressed and extruded from between north america and eurasia resulting in high levels of microseismic activity along the southern flanks of the chersky range

vertical averages of rheology of the continental lithosphere relation to thin sheet parameters abstract we calculate vertically averaged rheologies for continental lithosphere subject to biaxial states of stress corresponding to extensional compressional and strike slip regimes we assume that brittle deformation occurs by byerlee s law and that ductile deformation occurs by steady state creep mechanisms under these conditions the vertically averaged rheology may be approximated by a single power law over the range of strain rates from 10 13 to 10 17 s 1 this is true not only when either brittle or ductile behavior dominates the strength of the lithospere but also when the contributions of brittle and ductile behavior are comparable the value of the power law exponent n and the strength of the lithosphere depend on the depth and stress difference at the brittle ductile transition the temperature at the moho and the geothermal gradient in the range of vertically integrated stress between 1012 and 1013 n m 1n reaches a minimum of between 5 and 10 under other circumstancesn increases rapidly as brittle behavior dominates over ductile behavior the strength of the lithosphere depends most strongly on the thermal state of the upper mantle but also is influenced by the stress supported in brittle portions of the lithosphere

thermal structure beneath tohoku northeast japan combined with geophysical observations and petrological considerations we have studied the thermal models of the tohoku northeast japan subduction zone the temperature of the deeper part of the mantle wedge is constrained by the petrological requirements that demand a high temperature region 1400 c beneath the volcanic zone to satisfy this condition namely to explain the arc volcanism in the tohoku subduction zone a high temperature upper mantle beneath the japan sea possibly 1400 c at the depth shallower than 100 km is required heat flow distribution which constrains the temperature structure at shallower depth shows that low heat flow values ranging from 20 to 50 mw m2 are distributed between the trench and aseismic front and consistently high heat flow values 80 120 mw m2 prevail in the area landward of the volcanic front including the japan sea recent heat flow measurements indicate that the transition from low to high heat flow occurs between the aseismic front and volcanic front this heat flow distribution restricts the form of the flow in the mantle wedge and shear stress acting at the plate boundary a rigid non moving triangular region bounded by the subducting slab and the vertical plane beneath the aseismic front is necessary to explain the low heat flow zone between the trench and aseismic front shear stress of about 500 to 1000 bar is supposed to act at the interface between the subducting slab and rigid triangular mantle we also suggest that to explain the peak of heat flow near the volcanic front 120 mw m2 magma migration rate of about 0 05 cm3 cm 2 year 1 is necessary a comparison of the thermal structure of the tohoku subduction zone with that of the central south america which was preliminary calculated by the similar scheme shows that the temperature in the mantle wedge of the central south america is probably lower than that of the tohoku subduction zone this difference may result in that of the nature of volcanic rocks

the evolution of forearc structures along an oblique convergent margin central aleutian arc multichannel seismic reflection data were used to determine the evolutionary history of the forearc region of the central aleutian ridge since at least late miocene time this sector of the ridge has been obliquely underthrust 30 west of orthogonal convergence by the northwestward converging pacific plate at a rate of 80 90 km m y our data indicate that prior to late eocene time the forearc region was composed of rocks of the arc massif thinly mantled by slope deposits the forearc region probably lacked both major depositional basins and a tectonically attached accretionary prism of offscraped oceanic deposits beginning in latest miocene or earliest pliocene time a zone of outer arc structural highs and a forearc basin began to form formation of these companion intraarc structures may be linked to the late neogene growth of an accretionary wedge that formed as the result of the deposition of a thick turbidite wedge in the aleutian trench initial structures of the zone of outer arc highs formed as the thickening wedge underran compressively deformed and uplifted the seaward edge of the arc massif above a landward dipping backstop thrust forearc basin strata ponded arcward of the elevating zone of outer arc highs however most younger structures of the zone of outer arc highs cannot be ascribed simply to the orthogonal effects of an underrunning wedge oblique convergence created a major right lateral shear zone the hawley ridge shear zone that longitudinally disrupted the zone of outer arc highs truncating the seaward flank of the forearc basin and shearing the southern limb of hawley ridge an exceptionally large antiformal outer arc high structure slivers of forearc basement rocks and overlying strata have been transported along the shear zone that is flanked by differentially elevated structures attributed to localized transpressive and transtensional processes uplift of hawley ridge may be related to the thickening of the arc massif by westward directed basement duplexes in addition the forearc is disrupted by structures transverse to the margin that occur where unusually high stress accumulations have resulted in the rupture of repeated great earthquakes it is likely that many ancient active margins evolved in tectonic and depositional settings similar to those of the central aleutian ridge great structural complexity including the close juxtaposition of coeval structures recording compression extension differential vertical movements and strike slip displacement should be expected even within areas of generally kindred tectonostratigraphic terranes

fault trends on the seaward slope of the aleutian trench implications for a laterally changing stress field tied to a westward increase in oblique convergence 1 normal faults along the seaward trench slope sts commonly strike parallel to the trench in response to bending of the oceanic plate into the subduction zone this is not the circumstance for the aleutian trench where the direction of convergence gradually changes westward from normal to transform motion gloria side scan sonar images document that the aleutian sts is dominated by faults striking oblique to the trench west of 179 e and east of 172 w these images also show a pattern of east west trending seafloor faults that are aligned parallel to the spreading fabric defined by magnetic anomalies the stress strain field along the sts is divided into two domains west and east respectively of 179 e over the western domain sts faults and nodal planes of earthquakes are oriented oblique 9 46 to the trench axis and 69 90 to the magnetic fabric west of 179 e sts fault strikes change by 36 from the e w trend of sts where the trench parallel slip gets larger than its orthogonal component of convergence this rotation indicates that horizontal stresses along the western domain of the sts are deflected by the increasing obliquity in convergence an analytical model supports the idea that strikes of sts faults result from a superposition of stresses associated with the dextral shear couple of the oblique convergence and stresses caused by plate bending for the eastern domain most nodal planes of earthquakes strike parallel to the outer rise indicating bending as the prevailing mechanism causing normal faulting east of 172 w sts faults strike parallel to the magnetic fabric but oblique 10 26 to the axis of the trench on the basis of a coulomb failure criterion the trench oblique strikes probably result from reactivation of crustal faults generated by spreading

upper limits of power law creep of rocks at depths below a few kilometers most rocks probably flow in the steady state dominantly by diffusion assisted dislocation creep mechanisms whereby strain rate is proportional to a low power of deviatoric stress at high stress levels for metals ceramics and rocks the power law breaks down plb to some other strain rate stress functional relationship and this paper explores four methods by which the breakdown stress b may be estimated for metals alloys and ceramics b 10 3 for non silicate rocks b 5 10 3 and for silicates b 10 2 for silicates the low stress exponent power law breaks down at high stresses to an exponential dependence of strain rate upon stress or to a dorn law the governing flow mechanisms generally being lattice resistance controlled dislocation glide and mechanical twinning plb breakdown stresses for silicates are nearly strain rate insensitive and extrapolated values of b are so high as to suggest that low temperature high stress plasticity may be suppressed entirely in geological deformations plb for olivine rich rocks is discussed in some detail and new three dimensional deformation mechanism surfaces for wet and dry olivine polycrystals are presented

on the choice of boundary conditions in continuum models of continental deformation recent studies of continental deformation have treated the lithosphere as a viscous media and investigated the time evolution of the deformation caused by tectonic and buoyancy forces we examine the differences between continuum models that keep velocity boundary conditions b c constant with time and models that keep stress b c constant with time and demonstrate these differences by using a simple example of a continental lithosphere that is subjected to horizontal compression our results show that in the case of constant stress b c the indentation velocity decreases with time while in the case of constant velocity b c the indentation velocity remains constant with time observations from the tibetan plateau and from the andes indicate that the rate of indentation decreases with time this suggests that when buoyancy forces are comparable in magnitude to tectonic forces constant stress b c are more appropriate for time evolutionary models of continental lithosphere we propose a simple procedure that translates velocity b c which are easy to estimate into stress b c at the initial stage and keep the stress b c constant along the time progression of the calculations

roughness at the base of the seismogenic zone contributing factors the cutout depth of microseismic activity in continental fault zones appears to correspond to the onset of greenschist metamorphic conditions at about 300 c it can generally be modeled as the transition from frictional to quasi plastic behavior in quartzofeldspathic crust shear resistance increases with depth through the frictional regime to peak at the transition beneath which it falls off exponentially with increasing temperature larger earthquake ruptures ml 5 5 nucleate around this transition depth where the highest concentrations of strain energy may accumulate varying depth and amplitude of the peak shear resistance along strike induce fluctuations in strain energy concentration at the base of the seismogenic zone factors affecting the depth of the transition include crustal composition geometry and mode of faulting fluid pressure levels in the frictional regime and water content in the quasi plastic regime quasi plastic strain rate and geothermal gradient evaluation of their relative importance is complicated because several are interdependent however compositional change may cause abrupt irregularities in seismogenic depth and peak shear resistance while regional variations in heat flow look to be particularly effective in creating long wavelength heterogeneities in strain energy concentration affecting faulting style

the role of a heterogeneous inclusion during continental collision abstract we study a simple model of the deformation due to collision between two continental plates the model is patterned after the geometry of the india asia collision but it may be applicable to more general situations the collision is modelled by the indentation of a viscoplastic plate asia by a rigid punch india a finite element approach capable of dealing with finite deformation heterogeneity and non linear rheologies is used to study the role of an inclusion of increased rigidity on the viscoplastic flow of the plate the heterogeneous zone was chosen to approximate the shape and position of the tarim basin behind tibet the lateral boundaries of the plate were defined so as to approximate those that are believed to apply for the asian plate i e fixed boundaries to the west and north free boundaries to the west and south several geometries and contrasts of rheology were analyzed to verify the general validity of our results the purpose was to determine the main kinematic features of the flow not to model the particular case of tarim it is found that the heterogeneity affects the flow in a drastic way even for small contrasts of rigidity and for small size heterogeneities the main effect of the heterogeneity is to create two important shear zones or internal boundary layers the most important is in front of the inclusion separating the main flow from the inclusion as the flow proceeds the inclusion rotates and the shear zone progressively bends in qualitative agreement with the bending of the large strike slip faults like the altyn tagh that separates the tarim from tiber our model of continuous deformation does not allow for the formation of faults or fractures but the shear zones may be understood as the initial stages for the generation of fractures or velocity discontinuities another shear zone develops to the left of the inhomogeneity near the fixed boundary and in the wake of the punch the role of these shear zones is to isolate the heterogeneity from the main flow to facilitate its rigid rotation inside the flow of the softer plate a final observation is that in the presence of the heterogeneity the lower left hand corner rotates almost like a rigid block producing a sharp bend on the left hand side of the plate as the inclusion rotates pushed by the punch this rigid block rotates clockwise to open up the space necessary for the homogeneity

stress modeling of tectonic blocks at cape kamchatka russia using principal stress proxies from high resolution sar new evidence for the komandorskiy block abstract a new data set of 851 lineaments mapped from european remote sensing satellites 1 and 2 full resolution synthetic aperture radar sar data are interpreted as geological proxies for stresses resulting from plate and block collision near the cape kamchatka region of the kamchatka peninsula using a geographical information systems based analysis numerical manifold method nmm analysis is used to model the stress field within the cape kamchatka region resulting from the collision of lithospheric plates and blocks results of our nmm model using different plate motion and plate configuration in the region are compared with orientation data for the mapped set of lineaments these data suggest that the lineaments observed in sar cannot be fully explained by a simple two plate model in this in the cape kamchatka region as an alternative we propose that the data can be explained by the existence of the previously proposed komandorskiy block recent global positioning satellite measurements in the aleutian islands support our lineament derived model and show that the near islands komandorskiy island block of the extreme western aleutians is moving independently of the north american plate

structure of an active arc continent collision area the aleutian kamchatka junction abstract the plate boundary between the north american and pacific plates along the westernmost portion of the aleutian arc is formed by several active dextral strike slip fracture zones fracturing is expressed by fault escarpments visible in reflection seismic profiles and narrow zones of high seismicity the seismicity changes from the aleutian arc where shallow earthquakes occur to intermediate to deep earthquakes in the kamchatka aleutian junction area and to the south in the vicinity of this junction the distribution of earthquakes follows the traces of major aleutian fracture zones down to a depth of 100 km south of the kronotsky peninsula the distribution of earthquakes changes and is similar to other subduction zones the fracture zones cut the western aleutian arc into individual blocks which move along strike slip faults to the west and collide with kamchatka related transtension and transpression has led to the development of pull apart basins and block tilting the fracture zones splay into numerous strike slip faults to the west where they meet the kamchatka margin the main faults of the aleutian fault system can be traced on the kamchatka mys peninsula where deformed pleistocene sediments high uplift rates recent thrust faulting and changes in the drainage system provide evidence for neotectonic activity the fault pattern onshore and offshore in the junction area supports an indentation tectonic model of the collision of aleutian arc slivers with kamchatka

the komandorsky basin as a product of spreading behind a transform plate boundary abstract a long term study of the komandorsky basin in particular by the cruises of r vvulkanolog and r vdmitry mendeleev in the summer of 1988 made it possible to describe in detail the main structural features of this marginal basin several ne and n s trending spreading axes were detected in the basin in the rear of the transform sector of the aleutian arc spreading axes are shifted relative to each other along nw striking transform faults there is variation in style from north to south along the strike of the spreading axes and a less well defined symmetry compared to that typical of mid oceanic ridges the tectonic pattern and heat flow reflect several stages starting from the eocene in the history of the opening of the basin caused by the changes in the velocity and direction of the relative motion between the pacific and north american plates

constraints on the timing of crustal imbrication in the central trans hudson orogen from single zircon 207pb 206pb ages of granitoid rocks from the pelican thrust zone saskatchewan the pelican thrust is a major ductile high strain zone in the reindeer zone trans hudson orogen northern saskatchewan it is interpreted as the main sole thrust separating stacked juvenile paleoproterozoic allochthons and underlying archean microcontinental crust in this central part of the orogen exposed nonmylonitic rocks in the footwall of the thrust consist of the sahli monzocharnockite and the smaller more highly retrograded macmillan point granite protomylonitic to ultramylonitic gneisses in the thrust zone derive from a variety of prethrust protoliths a footwall internal suite mainly comprises quartzofeldspathic orthogneisses q gneisses and high grade migmatitic paragneisses hanging wall external suite mylonitic gneisses include feldspar porphyroclastic hornblendic grey gneisses probably derived from arc plutons and laminated amphibolites derived from volcanic rocks the overlying allochthon mainly comprises protoliths equivalent to those of the porphyroclastic orthogneisses and lam

large scale deformation associated with ridge subduction summary continuum models are used to investigate the large scale deformation associated with the subduction of aseismic ridges formulated in the horizontal plane using thin viscous sheet theory these models measure the horizontal transmission of stress through the arc lithosphere accompanying ridge subduction modelling was used to compare the tonga arc and louisville ridge collision with the new hebrides arc and d entrecasteaux ridge collision which have disparate arc ridge intersection speeds but otherwise similar characteristics models of both systems indicate that diffuse deformation low values of the effective stress strain exponent n are required to explain the observed deformation deformation is somewhat insensitive to the vertically integrated strength of the arc inversely proportional to the argand number ar but indicates that the arc lithosphere is not extremely weak ar 100 low values of both ar and n suggest that the thermal structure is typical of cold or normal arcs and that deformation is dominated by flow in the lower crust and mantle in addition low values of n approaching newtonian flow may indicate that specific deformation mechanisms dictate deformation of the arc lithosphere possible mechanisms include low stress grain size dependent creep pyroxenite controlled rheology and mechanisms associated with water weakening changes in the boundary conditions greatly affect deformation within island arcs high rates of arc ridge intersection speed tonga louisville system yield arc parallel tension and crustal thickening in the wake of ridge subduction in contrast low rates of arc ridge intersection speed new hebrides d entrecasteaux system yield compressional deformation directly arcward of the collision zone and transverse strike slip faulting adjacent to the region of compressional deformation localized regions of extensional deformation along the frontal part of the arc adjacent to the collision zone may contribute to the formation of re entrants

new finite element techniques for modeling deformation histories of continents with stratified temperature dependent rheology previous techniques for modeling anelastic deformation of the lithosphere have included plane strain models restricted to two dimensional problems and quasi three dimensional plane stress or thin plate models that did not accurately include the effects of the shallow frictional layer or of kinematic detachment of crust from mantle this paper presents techniques to remedy these deficiencies of thin plate models an iteration strategy in which the rheology is linearized using artificial prestress and a particular effective viscosity tensor causes the calculation of horizontal velocities to converge monotonically even with a frictional layer at the top of the lithosphere a technique using two planar grids allows deformation and displacement to be different at the crust and mantle levels at far less cost than that of a three dimensional grid a finite element technique is developed for computing the changes in thickness of these layers caused by pure shear simple shear and pressure gradients a technique based on relaxation of perturbation eigenfunctions solves the heat equation in the lithosphere during deformation accuracy of component numerical methods is good for simple test problems but in realistic nonlinear problems utilizing all components only the precision can be determined because of the lack of analytic solutions precision of the combined program is tested with a realistic sample problem and presented as a function of the number of iterations in each velocity solution convergence factor 0 73 to 0 88 size of time step in the predictor corrector time integration convergence as t0 8 and number of degrees of freedom in the finite element grid convergence as n 0 5 to 0 8 for most variables overall cost of a simulation varies with the fractional precision as 3 3 a new consequence of kinematic detachment a moving wave of crustal thickness is found unfortunately the form of the wave depends on the finite element size this means that element size must be chosen to approximate the smoothing by flexural rigidity effects that were neglected because of cost

earthquake source parameters and stress distribution in the adak island region of the central aleutian islands alaska source parameters have been systematically determined for all earthquakes with mw 5 that occurred between 172 w and 179 w longitude in the adak island region of the central aleutian islands during 1977 march 1987 we relocate the events using a plate velocity model developed for the region and use two alternative methods of depth determination the first method uses arrival times of direct and prominent reflected phases primarily pwp the second method uses broadband p wave displacement seismograms in an inversion for source depth the analysis of these events provides a well constrained data set for the study of stress release along an active subduction zone three earthquakes which occurred seaward of the trench axis are located just below the crust mantle interface and show extension nearly perpendicular to the trench axis seven events occurred in the wadati benioff zone sixty four events are located in the main thrust zone and except for five unusual events are characterized by thrust mechanisms with one nodal plane dipping north at a shallow angle in cross section the thrust zone appears as a thin 10 15 km thickness interplate region which extends from 15 to 50 km depth a 10 15 average discrepancy between observed slip vector azimuths and the predicted relative motion direction between the pacific and north american plates exists for these events the observed slip vectors are oriented more normal to the trench than is predicted by plate motions five earthquakes which occurred as aftershocks to the may 7 1986 earthquake mw 8 0 in the crust of the overriding plate have strike slip mechanisms consistent with right lateral motion on arc parallel fault planes observations of slip vectors along the whole aleutian arc show a similar trend to that observed in the adak island region the largest differences 30 between observed and predicted slip azimuths occur around 175 e we propose a model of plate interaction in which a portion of the along arc motion occurs along a weak strike slip shear zone in the upper plate near the volcanic line the slip azimuths in the main thrust zone fit this model well if the amount of transcurrent slip occurring in the upper plate is 60 of the arc parallel relative plate motion a consequence of the model is along arc extension of the overriding plate between the accretionary wedge and the volcanic line especially in the western part of the aleutian arc calculations based on a tectonic model by geist et al 1988 for the formation of arc summit basins through block rotation and translation suggest that along arc extension has been significant since late miocene or early pliocene

on the possibility of newtonian flow in the upper mantle abstract the possibility of a linear creep newtonian viscosity upper mantle is reexamined on the basis of present knowledge on flow mechanisms in olivine including differences between activation parameters for creep and for diffusion and revised estimates of grain boundary width results of the comparison between linear superplastic creep and power law creep are presented as crossover temperature between lattice and grain boundary diffusion and crossover stress between non linear and linear creep as a function of temperature depth grain size grain boundary width grain boundary diffusion activation energy and rate controlling species silicon or oxygen for the most realistic values of the parameters linear creep is well within the range of possibilities there is no major objection from rheology to the idea of a newtonian viscosity upper mantle viscosities for the two kinds of creep are compared for laboratory upper mantle and tentatively lower mantle conditions this results in a prediction of the grain size at which coble creep could be observed in olivine polycrystals in the laboratory 10 m or less the upper mantle viscosity is in the 1020 1021 pa s range minimum estimates for the viscosity of the lower mantle are of the same order the pictures of the rheology of the mantle derived from microphysical models of flow and from geophysical observation can therefore be made compatible

terrestrial heat flow in the northwestern pacific abstract on the basis of 150 new heat flow measurements the thermal structure of the uppermost layer of the earth in the transition zone of the northwestern pacific and relation of heat flow values to the age of tectonic magmatic activity have been established an analysis was made of the main factors distorting the geothermal field in the region under study and a map of measured and background heat flow values was constructed calculations were made to determine the temperatures in the deep crust and upper mantle along the dss profiles in the northwestern pacific and a map was constructed showing the depth of the zone of fractional melting 1200 c isotherm identified with the upper boundary of the asthenospheric layer

magma evolution and open system processes at shiveluch volcano insights from phenocryst zoning phenocryst zoning patterns are used to identify open system magmatic processes in the products of the 2001 eruption of shiveluch volcano kamchatka the lavas and pumices studied are hornblende plagioclase andesites with average pre eruptive temperatures of 840c and fo2 of 1 5 2 1 log units above nickel nickel oxide nno plagioclase zoning includes oscillatory and patchy zonation and sieve textures hornblendes are commonly unzoned but some show simple multiple or patchy zoning apatite microphenocrysts display normal and reverse zoning of sulphur the textural similarity of patchy hornblende and plagioclase together with ba sr systematics in patchy plagioclase indicate that the cores of these crystals were derived from cumulate material plagioclase liquid equilibria suggest that the patchy texture develops by resorption during h2o undersaturated decompression when h2o saturated crystallization recommences at lower pressure reduced ph2o results in lower xan in plagioclase causing more al rich hornblende to crystallize plagioclase cores with diffuse oscillatory zoning and unzoned hornblende crystals probably represent a population of crystals resident in the magma chamber for long periods of time in contrast oscillatory zoning in the rims of plagioclase phenocrysts may reflect eruption dynamics during decompression crystallization increasing fe al in oscillatory zoned rims suggests oxidation as a result of degassing of h2o during decompression a general lack of textural overlap between phenocryst types suggests that different phenocryst populations were spatially or temporally isolated during crystallization we present evidence that the host andesite has mixed with both more felsic and more mafic magmas olivine and orthopyroxene xenocrysts with reaction or overgrowth rims and strong normal zoning indicate mixing with basalt sieve textured plagioclase resulted from mixing of a more felsic magma with the host andesite the mineralogy and mineral compositions of a mafic andesite enclave are identical to those of the host magma which implies efficient thermal quenching and thus small volumes of intruding magma mixing of this magma with the host andesite results in phenocryst zoning because of differences in dissolved volatile contents we suggest that small magma pulses differentiated at depth and ascended intermittently into the growing magma chamber producing incremental variations in whole rock compositions

terrestrial heat flow estimation of deep temperature and seismicity of the kamchatka region abstract the results of heat flow measurements in the kamchatka region are discussed in this paper on the basis of more than 60 continental measurements and a number of heat flow measurements made at sea a schematic heat flow map has been compiled an interpretation of the anomalous fraction of heat flow is made and a likely temperature distribution in the lithosphere of kamchatka was calculated based on 3 d modelling of the thermal field an analysis is made of the geothermal regime of the kamchatka sector of the focal zone temperature is shown to be the controlling factor for the maximum depth of occurrence of seismic events the majority of earthquakes occur in a layer 0 100 km thick where the temperature does not exceed 700 c

oblique plate convergence slip vectors and forearc deformation slip vectors from thrust earthquakes at subduction zones where convergence is oblique to the trench often point between the directions of relative plate convergence and normal to the trench axis suggesting that oblique convergence is taken up by partial decoupling decoupling means that a component of arc parallel motion of the leading edge of the upper plate results in less oblique thrusting at the trench partial decoupling is modeled by partitioning of oblique convergence into slip on thrust and strike slip faults that are parallel to the trench and to each other and starting with a force equilibrium condition a relationship between the obliquity and the earthquake slip vector orientation is derived assuming that either fault slips when shear stress on it reaches a yield stress oblique slip parallel to the plate vector should occur on the thrust fault when obliquity is smaller than a critical angle for obliquity at or greater than this angle the stress on the strike slip fault is large enough to start it slipping and when both faults are active the arc parallel motion of the forearc deflects the slip vector back toward the trench normal if we assume that continued slip on either fault occurs at constant stress but the two faults can be at different stresses the slip vector will maintain a constant angle relative to the trench normal even when obliquity is larger than the critical angle this limiting angle of the slip vector called max measured relative to the trench normal is simply the arcsine of the ratio of the shear forces resisting slip on the strike slip and thrust faults a consequence is that when the obliquity exceeds max the slip vectors on the thrust fault are sensitive only to the thrust fault orientation and contain no information about the convergence direction between the plates slip vectors at the java trench southwest of sumatra show the relationship clearly with max 20 5 while slip vectors at the aleutian trench show the relationship less clearly with max 25 to 45 the greater angle at the aleutian trench suggests that the upper plate is stronger in the aleutian arc relative to the thrust fault than in the sumatran arc consistent with the sumatran arc being continental and having a well developed strike slip fault while the aleutian arc is oceanic and without a clear transcurrent fault slip vectors at the philippine trench which like sumatra has a large strike slip fault inboard of it tend to stay within 25 of the trench normal when obliquity is as large as 50 if obliquity exceeds max and continues to increase along a subduction zone the rate of motion of the forearc relative to the upper plate will vary with obliquity in which case the forearc sliver should extend or contract parallel to the arc from the geometry of modern island arcs arc parallel extension should be the more common and has been hypothesized for both sumatra and the aleutians on the basis of earthquake slip vectors and for these and other arcs from geological observations from estimates of max and the arc parallel gradients in obliquity arc parallel strain rates are estimated to be 1 to 3 10 8 yr for the sumatran forearc 2 to 6 10 8 yr for the aleutian forearc and 0 3 to 3 10 8 yr for the philippine forearc oblique convergence and subsequent arc parallel extension if accompanied by crustal thinning may provide an important yet little appreciated mechanism for bringing high grade metamorphic rocks to the surface of subduction complexes

numerical modeling of intraplate deformation simple mechanical models of continental collision we study a model of continental collision in which one of the continents acts as a rigid die indenting the other plate which flows as an incompressible viscoplastic medium we consider two extreme cases of plane deformation 1 plane strain which corresponds to an infinitely thick lithospheric plate and 2 plane stress corresponding to a very thin plate deformation of the lithosphere a thick plate should be intermediate between those extremes we found that the flow in the plane strain case is quite similar to that obtained by slip line theory the plane stress results are quite different since in this case most of the plate shortening is taken up by the thickening of the lithosphere we also explored the role of boundary conditions on the flow in particular the role of the side walls containing the flow of the lithosphere in the case of a free lateral boundary the main feature is a flow of matter toward this free wall and a s like pattern for the horizontal stress field for a rigid wall on the other hand the plane strain and the plane stress results are quite different in the first case there is a large return flow on the sides of the punch the material being extruded along the only free surface available in the plane stress case the return flow disappears and the material displaced by the penetration of the die tends to thicken the plate the role of a nonlinear constitutive relation is studied for power law creep as the power of the flow limit increases the flow retains its general features but the deformation localizes creating sharper contrasts between high and low strain rate areas in plane stress the effect of nonlinearity is to increase the contrasts in vertical motion available data for asia are discussed in the light of the new results

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

centroid moment tensor solutions for october december 1991 centroid moment tensor solutions are presented for 160 earthquakes that occurred during the fourth quarter of 1991 the solutions are obtained using corrections for aspherical earth structure represented by a whole mantle shear velocity model sh8 u4l8

finite strain calculations of continental deformation 2 comparison with the india asia collision zone numerical experiments on a thin viscous sheet model for deformation of continental lithosphere subjected to an indenting boundary condition yield distributions of crustal thickness of stress and strain rate and of latitudinal displacements that may be compared with observations in the india asia collision zone a simple indenting boundary condition applied to initially laterally homogeneous sheets obeying a power law rheology produces results that are in broad agreement with the observations provided that the power law exponent is three or greater and the sheet can support vertically integrated stress differences of 2 1013 5 1012 n m 1 in regions in front of the indenter under these conditions the calculated deformation shows accommodation of convergence primarily by crustal thickening to produce a plateau in front of the indenter palaeomagnetic data from india and tibet and the observed distribution of topography suggest that much of the post eocene convergence of india with asia has been taken up by deformation within asia that involved crustal thickening the principal difference between calculation and observation is the absence from the calculated strain rate fields of east west extension of the plateau in front of the indenting boundary the calculations show that once such a plateau is formed the buoyancy force associated with the crustal thickness contrast inhibits further thickening and the plateau strains at less than half the rate of its immediate surroundings seismically determined regional strain rates exhibit a similar distribution with the tibetan plateau straining at about one quarter the rate of the tien shan and ningxia gansu regions calculated principal compressive stress orientations and regional strain rates agree with the seismically determined quantities in the mongolia baikal tien shan tibet and ningxia gansu regions of asia to within the uncertainty of the latter the vertically integrated stresses that are calculated for the viscous sheet are comparable with those that can be supported by a theologically stratified continental lithosphere obeying laboratory determined flow laws we suggest that the thin viscous sheet model described in this paper and its companion gives a simple and physically plausible description of the observed deformation in central asia in this description the predominant mechanism of accommodation of continental convergence is diffuse crustal thickening with shear on vertical planes playing a subsidiary role once large crustal thickness contrasts have been established

application of continuum models to deformation of the aleutian island arc continuum models were constructed to describe large scale deformation of the aleutian island arc over the past 5 m y these models consider the island arc as a continuum in the horizontal plane with the velocity boundary condition at the pacific edge stated as a fraction of pacific plate convergence transferred to the arc first a simple model of uniformly distributed strain is formulated to illustrate the mechanics of continuous deformation lineaments along the arc massif rotated about a vertical axis are matched by small element rotation calculated from the model however this model does not predict across arc variations in deformation and produces an unrealistic amount of crustal thickening after 5 m y a physically more meaningful model of deformation is the thin viscous sheet model based on averages of stress and rheology throughout the lithosphere the amount of motion transferred from the pacific plate to the arc is constrained by the rotated lineaments while the effective stress strain exponent n and the ability the lithosphere has to sustain crustal thickness contrasts the argand number are independent variables primarily bathymetry earthquake focal mechanisms and styles of faulting are used to evaluate the models the preferred model is one where the amount of motion transferred from the pacific plate is greater in an arc parallel direction than in an arc normal direction producing stresses consistent with strike slip faulting at the far western end of the arc and tensional stresses consistent with transverse normal faulting elsewhere in the arc massif this model agrees with observations of slip vectors by ekstr6m and engdahl 1989 who conclude that a portion of the arc parallel component of relative plate motion is taken up in the overriding plate this model implies that compressive stress transferred to the arc is small in comparison to along arc shear stress and that stresses conducive to strikeslip faulting are prevalent throughout the arc

large scale deformation related to the collision of the aleutian arc with kamchatka the far western aleutian island arc is actively colliding with kamchatka westward motion of the aleutian arc is brought about by the tangential relative motion of the pacific plate transferred to major right lateral shear zones north and south of the arc early geologic mapping of cape kamchatka a promontory of kamchatka along strike with the aleutian arc revealed many similarities to the geology of the aleutian islands later studies support the notion that cape kamchatka is the farthest west aleutian island and that il has been accreted to kamchatka by the process of arc continent collision deformation associated with the collision onshore kamchatka includes gravimetrically determined crustal thickening and formation of a narrow thrust belt of intensely deformed rucks directly west of cape kamchatka the trend of the thrust faults is concave toward the collision zone indicating a radial distribution of maximum horizontal compressive stress offshore major crustal faults trend either oblique to the kamchatka margin or parallel to major aleutian shear zones these offshore faults are complex accommodating both strike slip and thrust displacements as documented by focal mechtmisms and seismic reflection data earthquake activity is much higher in the offshore region within a zone bounded to the north by the northernmost aleutian shear zone and to the west by an apparent aseismic front analysis of focal mechanisms in the region indicate that the present day arc continent contact zone is located directly east of cape kamchatka in modeling the dynamics of the coiiision zone using thin viscous sheet theory the rhealogical parameters are only partially constrained to values of n the effective power law exponent 2 3 and ar the argand number s 30 these values are consistent with a fnrcarc thermal profile of kamchatka previously determined from heat flow modeling the thin viscous sheet modeling also indicates that onshore thrust faulting is a consequence not only of compressive stresses resulting from the west directed collision but also of sediment induced coupling of the subducting pacific plate

palaeomagnetism of palaeogene volcanic series of the kamchatskymys peninsula east kamchatka the motion of an active island arc abstract the paleocene tarkhov and eocene baklan formation are the lower and upper members respectively of a continuousisland arc sequence at the kamchatsky mys peninsula on the pacific side of the kamchatka peninsula we studied paleocene tuffs at five localities and eocene tuffs and basalts at nine localities from beds of various attitude thermal demagnetization revealed consistent directions of the characteristic component chrm of natural remanent magnetization in both formations normal and reversed chrm directions in the paleocene rocks are roughly antiparallel thus constituting a positive reversal test at the same time the fold test mcfadden and jones 1981 is not straightforward as locality means display a banana shaped distribution most probably because of local rotations inclination only data however show considerable decrease in dispersion after tilt correction thus constituting a positive fold test and implying a prefolding and most prprobably primary age of the chrm in these rocks the formation mean inclination of 57 5 3 8 calculated with the aid of inclination only statistics corresponds to a palaeolatitude of 38 1 4 1 n chrm directions isolated from the eocene baklan formation are well grouped and pass both the fold and reversal tests the formation mean inclination of 65 0 4 9 corresponds to a palaeolatitude of 47 0 6 4 n geological data point to uninterrupted accumulation of the studied island arc complex and hence to continuous activity of the subduction zone and related island arc the difference between the paleocene and eocene palaeolatitudes of 8 9 6 3 is statistically significant thus implying an absolute motion of the active subduction zone and related island arc during the paleocene to eocene no other active subduction zone existed between the studied area and eurasia we think that the studied island arc belonged to the pacific plate and was moving northwestward while the oceanic part of the eurasian or north american plate was subducted southward underneath this arc this island arc became inactive in the late eocene and was transported with the pacific plate which started subducting under the eurasia margin since that time

rheology of the lithosphere selected topics we review recent results concerning the rheology of the lithosphere with special attention to the following topics 1 the flexure of the oceanic lithosphere 2 deformation of the continental lithosphere resulting from vertical surface loads and forces applied at plate margins 3 the rheological stratification of the continents 4 strain localization and shear zone development and 5 strain induced crystallographic preferred orientations and anisotropies in body wave velocities we conclude with a section citing the 1983 1986 rock mechanics literature by category

the relationship between strain rates crustal thickening palaeomagnetism finite strain and fault movements within a deforming zone the deformation of continental lithosphere generally does not occur on a single fault but is distributed over a zone up to a few hundred kilometers wide if the deformation rate is constant throughout the zone and if the velocity is continuous everywhere and does not change with time the deformation gradient tensor can be obtained analytically this tensor directly determines the amount of crustal thickening and is related in a complicated manner to the finite strain tensor the declination of remnant magnetization measured within the zone can also be rotated by the deformation because palaeomagnetic measurements are generally made on rigid blocks which have not been deformed there is no simple relationship between such rotations and the finite strain expressions for all these quantities are obtained in terms of the velocity gradients which are taken to be constant the history of western north america is used to illustrate the relationship between the present deformation by strike slip and thrust motions the rotations of the remnant magnetization vectors and the shape of accreted terranes at the surface distributed deformation is taken up by motion on faults which separate the zone into blocks the displacement rate on the faults their strike and the finite rotations of the blocks are all controlled by the velocity gradients analytic expressions relating these quantities are obtained a deforming zone containing normal faults which crosses central greece takes up both strike slip and extension and illustrates the relationship between the structures within and the motion across the zone a similar but more complicated zone in northeastern iran involves both strike slip and shortening

borobudur a basin under volcanic influence 361 000 years bp to present abstract borobudur basin is located in central java indonesia 30 km to the north of yogyakarta city the basin is famous for its unesco world heritage temple and for one of the world s most active volcanoes merapi located to the east of borobudur basin merapi is one of the three andesitic volcanoes that surround the basin merapi merbabu and sumbing volcanoes therefore volcanic activity has strongly influenced the evolution of borobudur basin the object of this contribution is to present the evolution of borobudur basin over the last 161 000 years in the light of volcanic influence the methodology and tools developed for this research span over different areas of expertise from geochemistry geology and geomorphology to remote sensing gis and archeology results highlight the following points 1 two major volcanic events deposited volcaniclastic materials up to tens of meters thick 119 000 years bp and 31 000 years bp in the southern part of the borobudur basin the second volcanic event could correspond to the collapse of the older ancient merapi camus et al 2000 or proto merapi stage newhall et al 2000 2 there is no trace in the borobudur basin of a large debris avalanche 3 there are several generations of paleolakes in the borobudur basin the latest one has shrunk until historical times corroborating the theory of newhall et al 2000 and murwanto et al 2004 that borobudur temple was standing by a water body most of these paleolakes were impounded following volcanic events paleolakes most probably originated from the blockage of the hydrographic network by volcanic material 4 borobudur temple was never buried under volcanic material during historic times

finite strain calculations of continental deformation 1 method and general results for convergent zones a finite element method is used to calculate stresses and strain rates in a thin viscous sheet representing the continental lithosphere that is bordered in one part by an indenting boundary we calculate solutions for velocity fields crustal thickness distributions strain rates rotation time integrated vorticity and finite strain ellipsoids to show the dependence of these quantities on the controlling parameters the stress strain exponent n and the strength of the sheet relative to the gravitational forces which may be expressed as a dimensionless number the argand number ar for n 3 and ar between 1 and 10 a plateau of thickened crust is formed in front of the indenter the strain histories of individual elements within the viscous sheet may be quite complex and the finite strain at the end of deformation may bear no obvious relation to the principal stress orientations during the deformation the strain rate fields are interpreted in terms of the style of faulting that would be expected in the brittle upper crust if it were coupled to the ductile but stronger upper mantle in general the length scale of deformation increases with ar and time since the start of the collision and decreases with n for ar 3 the increase in length scale is accompanied by a change in the style of deformation from predominantly crustal thickening with thrust faulting in the brittle layer to predominantly transcurrent deformation strike slip faulting in the brittle layer

a continuum model of continental deformation above subduction zones application to the andes and the aegean large scale continental tectonics of back arc extensional and andean type compressional environments are investigated by using the thin viscous sheet model to calculate the deformation within continental lithosphere that is subjected to horizontal forces on its plate boundaries and to basal drag from the asthenospheric flow beneath the shear tractions acting on the base of a deformable lithosphere are derived from a corner flow model that assumes a rigid subducting plate and a deformable overlying plate because the calculated shear tractions and the deformation within the overlying plate are interdependent the corner flow and the thin viscous sheet models are solved simultaneously we use a perturbation method to obtain analytical solutions for the velocity and strain rate fields within the overlying continental lithosphere the solutions depend on the angle of subduction the dimensionless thickness of the lithosphere and the ratio of asthenospheric to lithospheric viscosities which governs the viscous coupling between the asthenosphere and the lithosphere calculations are compared with observations from the andes and the aegean our results explain some of the features of the deformation in these regions that have heretofore not been explained by other models our model predicts that in a compressional environment a broad region of uplifted topography will tend to develop above a more steeply dipping slab 30o rather than above a shallower slab 10o 15o this is in accord with observations in the various segments of the central andes for an extensional environment the model predicts that a zone of compression can develop near the trench and that extensional strain rate can increase with distance from the trench as is observed in the aegean we also estimate the effective viscosities of 102o pa s for the asthenosphere 2x 1021 pa s for the aegean lithosphere and 1022 pa s for the andean lithosphere

provenance of the palaeocene eocene clastic rocks of the komandorsky islands the palaeocene eocene sequence on the komandorsky islands is represented by sediments and minor volcanics with a total thickness of more than 3000 m sandstone siltstone mudstone siliceous rocks and conglomerate which have accumulated in the upper and middle parts of a deep water fan prevail all the debris was delivered from the north east from the side of the contemporary komandorsky basin in the bering sea products of volcanic denudation and fresh tephra dominate its composition however at some levels the sandstone contains fragments of jasper quartzite slate schist gabbro and granite as well as non volcanic quartz along with a volcanic component the heavy minerals fraction of the sandstone contains considerable quantities of apatite zircon garnet and chromite as well as pyroxenes amphiboles and magnetite the chemical composition of the garnet is specific to metamorphic rocks a significant admixture of chromite is indicative of denudation from an ophiolite complex the folded and partly metamorphosed pre palaeogenic rocks were a source of these fragments their most likely source terrane was the shirshov bowers chain cretaceous palaeoarc the komandorsky block originally was the forearc part of the aleutian arc where the latter was superimposed on the shirshov bauers chain a subsequent north west displacement of forearc structures in response to the transform fault at the boundary between the pacific and north american plates could have resulted in the separation of the komandorsky block from its early palaeogene debris source area

the origin of summit basins of the aleutian ridge implications for block rotation of an arc massif it is proposed that many summit basins along the aleutian arc form from the clockwise rotation of blocks of the arc massif summit basins are arc parallel grabens or half grabens formed within the arc massif and are commonly located near or along the axis of late cenozoic volcanism geomorphically the aleutian arc appears to consist of contiguous rhombic blocks of varying size tens to hundreds of kilometers in length the boundaries between adjacent blocks are delineated by fault controlled canyons that cut the southern slope of the arc transverse to its regional trend evidence that these blocks have rotated clockwise is provided by the triangular shaped summit basins bordering the blocks to the north oblique physiographic trends offsets in the summit platform and broad deflections in the southern slope of the arc we present a model for block rotation that involves translation of blocks parallel to an arc it is suggested that block rotation which appears to have accelerated in late cenozoic time is linked to 1 a shift in the euler pole for the pacific plate 2 the consequential start up of late cenozoic volcanism 3 improved interplate coupling instigated by sediment flooding of the aleutian trench and 4 westward subduction of northeast striking segments of the inactive kula pacific ridge

orogeny in arc continent collision the banda arc and western new guinea eastern indonesia contains a modern continent island arc collision that demonstrates how a complex juxtaposition of continental and oceanic elements can occur in an orogenic belt shallow earthquakes show that strike slip faulting transports large crustal blocks into out of and along the collision zone while crustal shortening and thickening occur by steep angle thrust faulting near the strike slip faulting whereas strike slip faulting is paramount in positioning australian continental rocks so as to enclose the oceanic and island arc rocks of eastern indonesia its role may be overshadowed by the contemporaneous thrusting and may confuse geologic interpretations of the resulting orogenic belt

initial rupture and the main fault of earthquakes a comparison of the body wave first arrivals and cmt data for the kamchatka commander region abstract centroid moment tensor cmt data by dziewonski and coworkers and body wave first arrival earthquake parameters obtained from local networks were simultaneously investigated for 55 earthquakes in the kamchatka commander region which occurred from 1977 to 1986 the magnitude range considered was 4 5 7 5 the focal depths were 0 190 km the following parameters were studied origin time hypocenter coordinates focal mechanism and seismic moment the first arrival data were referred to the initial rupture in the earthquake focus and cmt estimates to the main earthquake fault the principal results obtained are as follows shifts of the centroid location with respect to initial hypocenters and the change or stability of the type of faulting from the initial rupture to the main fault depend on the tectonic position of the hypocentral area the character of the source time function depends on the sign of the main fault propagation i e downwards upwards horizontally

constraints on extension of continental lithosphere the extension of the continental lithosphere is considered in terms of the deformation of a thin viscous sheet with a temperature dependent rheology which is governed primarily by the strength of the upper mantle the lithosphere is subjected either to a constant strain rate or to a strain rate decreasing with time corresponding to the stretching of a region of fixed volume with a constant velocity boundary condition the progress of the extension is influenced by the decrease in strength of the lithosphere owing to its attenuation and by an increase in strength resulting from cooling the relative importance of these two effects depends on two parameters a dimensionless strain rate and the total strain for a given strain rate the extending lithosphere goes through two phases as the strain increases in the first the effect of thinning the lithosphere is predominant and the average strength of the lithosphere remains close to and perhaps slightly below its original strength in the second the effect of cooling of the lithosphere leads to very rapid increase in strength as a function of strain resulting presumably in an effective limit on the degree of stretching that the lithosphere experiences for the case of strain rate decreasing with time when the initial strain rates are of order 10 15s 1 this maximum extension is between 100 and 150 for a strain rate of 10 l4 s 1 the maximum lies around 400 to 500 at higher strain rates than this very large extension may occur before cooling significantly affects the strength of the lithosphere this model accords well with the constraints available from the continental margins of the north atlantic

an experiment in systematic study of global seismicity centroid moment tensor solutions for 201 moderate and large earthquakes of 1981 data from the global digital seismograph network were used to obtain centroid moment tensor solutions using the method of dziewonski et al 1981 results were obtained for 201 earthquakes ranging in seismic moment from 7 1023 to 3 1027 dyne cm the wide dynamic range of the sro asro stations allows us to investigate using the same algorithm series of events among which the smallest and the largest may differ in moment by a factor as large as 1000 among the events studied is a particularly interesting series in january 1981 off the island of honshu the main shock was preceded by three foreshocks and followed by four aftershocks all falling within the range of our analysis there is an indication that the source mechanism changes with the size of the event statistical analysis of the results for all earthquakes reveals that the half duration of shallow and intermediate events is approximated by 1 6 10 8 m01 3 the shifts in origin time for deep focus earthquakes are indicative of the general complexity of these events which depends only weakly on the seismic moment the extreme departures from a double couple mechanism seem to depend both on the moment and on the focal depth while for shallow earthquakes these departures decrease from a maximum at about 1025 dyne cm with increasing moment the reverse appears to be true with respect to deep focus earthquakes shallow earthquakes north of new guinea and along the solomon islands show systematic and substantial departures from the double couple mechanism

cold moho boundary beneath island arc systems an example from the north kamchatka arc a suite of amphibole bearing ultramafic and mafic xenoliths representing the crust mantle boundary was sampled in the northern segment of the kamchatka arc these nodules exhibit low equilibration temperatures at elevated pressures suggesting relatively cold conditions at moho depths which are uncommon in modern island arcs a model including successive stages of subduction initiation subduction cessation resulting in cooling of the lithospheric wedge followed by intra arc rifting due to upwelling of astenospheric diapirs is proposed to explain the observed textural petrologic and chemical characteristics of sampled rock assemblage

paleomagnetism of paleogene rocks of the central east kamchatka and komanorsky islands tectonic implications abstract paleomagnetic investigation of paleocene eocene to miocene rocks has been carried out at five localities in the eastern kamchatka and komandorsky islands after thermal and af cleanings consistent paleomagnetic directions were isolated for a number of formations data from the kronotsky peninsula and medny island passed both the fold and reversal tests and thus can be considered as reliable other results are of lower quality all the results were derived from sedimentary rocks and thus the corresponding paleolatitudes may be biased due to inclination errors a limited collection of basalts yielded inclination values systematically lower than those in sediments this anomaly was tentatively explained by shape anisotropy the correctness of paleolatitude estimations was thus left unproved several models of the northern pacific tectonic evolution were suggested and analyzed but the authors failed to find a model that was fully compatible with paleomagnetic kinematic and geological data

effects of a temperature dependent rheology on large scale continental extension the influence of a temperature dependent rheology on large scale continental extension is investigated using a thin viscous sheet model a vertically averaged rheology is used that is consistent with laboratory experiments on power law creep of olivine and that depends exponentially on temperature the behavior of the model depends principally on two parameters the peclet number which describes the relative rates of advection and diffusion of heat and a dimensionless activation energy which controls the temperature dependence of the rheology numerical calculations of such a thin sheet subjected to an extensional boundary condition show the following results at short times following the beginning of extension deformation occurs with negligible change in temperature so that only small changes in lithospheric strength occur due to attenuation of the lithosphere and maximum rates of deformation are located close to the extensional boundary however after a certain time interval thermal diffusion results in lowered temperatures in the lithosphere strongly increasing lithospheric strength and slowing the rate of extension this critical time depends principally on the peclet number and is short compared with the thermal time constant of the lithosphere the changes in strength cause the locus of high extensional strain rates to shift with time from regions of high strain to regions of low strain results of the calculations are compared with observations from the aegean where maximum extensional strains are found in the south near crete but the largest present day extensional strain rates occur about 300 km further north the comparisons support the hypothesis that the observed separation of the loci of maximum strain and maximum present day strain rates in the aegean may be the consequence of changes in lithospheric strength arising from the temperature dependent mechanical properties of lithospheric materials

origin transport and emplacement of an exotic island arc terrane exposed in eastern kamchatka russia the regional stratigraphy of eastern kamchatka includes an exotic early late cretaceous ophiolite and late cretaceous island arc volcanic sequence integrating the existing geologic and geophysical data we examine the origin transport emplacement and postemplacement deformation of the island arc terrane which is named the olyutorsky island arc results from several paleomagnetic studies consistently indicate that the island arc terrane originated 1000 km to the south of where it is presently exposed although the formative paleolatitudes of the island arc rocks approximately correspond to the location of the izanagi farallon subduction zone the age of the volcanic rocks postdates the cessation of izanagi farallon convergence thus indicating that an unnamed plate or back arc basin existed in the northwest pacific during late cretaceous time we examine two possible models for northward transport of the island arc terrane to kamchatka 1 infra oceanic transport with the pacific or kula plates and 2 coastwise translation of the island arc terrane after accretion to the eurasian margin far to the south of kamchatka for both models the dominant eocene and miocene deformation ages observed in eastern kamchatka are used as two possible age limits for the cessation of northward transport although the observed paleolatitudes from paleomagnetic data correspond best with the infra oceanic transport model the provenance of the paleogene transport stratigraphy indicates a near shore sediment supply our preferred interpretation is that the island arc terrane 1 accreted onto the eurasian margin concurrent with cessation of island arc volcanism maastrichtian danian and 2 underwent northward coastwise translation along a major strike slip fault zone ending by middle late eocene time 43 50 ma it is unclear whether the ophiolite was exposed during arc continent collision or whether the ophiolite was obducted onto the island arc prior to collision a consequence of either infra oceanic transport or coastwise translation is that an open corridor between the western terminus of the aleutian arc and kamchatka must have existed until middle to late eocene time spreading within the komandorsky basin subduction of sea mounts and collision of the aleutian arc with kamchatka are proposed to have instigated the second miocene phase of deformation which uplifted and reexposed the island arc terrane

subduction zone recycling processes and the rock record of crustal suture zonesthis article is one of a series of papers published in this special issue on the theme lithoprobe parameters processes and the evolution of a continent offshore observations at modern ocean margin subduction zones omszs reveal that bodies of accreted material are commonly volumetrically small or missing that crustal thinning and subsidence 3 5 km has occurred and that most trench axes lie close 5 30 km to the seaward tapering edge of coastal basement rock onshore mapping commonly documents missing or only narrow terranes of former forearc rock and the inboard migration of the arc magmatic front these observations are evidence that subduction is accompanied by the removal of sediment and crustal material from the submerged forearc by the kindred tectonic processes respectively of sediment subduction and subduction erosion subduction erosion truncates the margin migrates the trench inboard at 2 5 km ma onshore observations at ancient crust suturing subduction zones csszs imply that collisional suturing is accompanied by sediment subduction and truncation of both upper and lower plates during a protracted period of suturing 20 50 millio

deformation of the central andes 15 27 s derived from a flow model of subduction zones a simple viscous flow model of a subduction zone is used to calculate the deformation within continental lithosphere above a subducting slab our formulation accounts for two forces that dominate the deformation in the overriding lithosphere tectonic forces and buoyancy forces the tectonic forces arise from the subduction of oceanic lithosphere and produce minor extension in the lithosphere near the trench and compression farther inland buoyancy forces can induce horizontal pressure gradients due to variations in crustal thickness which cause lithospheric extension where the crust is thick numerical solutions obtained by using a finite element technique are compared with observations from the central andes 15 27 s the model predicts the observed deformation pattern of extension in the forearc compression in the western monocline corresponding to magmatic activity extension in the altiplano compression in the eastern monocline and subandes and no deformation in the brazilian shield by comparing the calculated solutions with the large scale tectonic observations we are able to evaluate the forces that govern the deformation in the central andes the approximately constant subduction velocity in the past 26 m y suggests that the rate of crustal shortening in the andes has decreased with time due to the thickening of the crust

formation of the rocky mountains western united states a continuum computer model one hypothesis for the information of the rocky mountain structures in late cretaceous through eocene time is that plate of oceanic lithosphere was underthrust horizontally along the base of the north american lithosphere the horizontal components of the motion of this plate are known from paleomagnetism and the edge of the region of flat slab can estimated from reconstructed patterns of volcanism new techniques of finite element modeling allow prediction of the thermal and mechanical effects of horizontal subduction on the north american plate a model that has a realistic temperature dependent rheology and a simple plane layered initial condition is used to compute the consequences of horizontal underthrusting in the time interval 75 million to 30 million years before present successful prediction of this model include i the location amount and direction of horizontal shortening that has been inferred from laramide structures ii massive transport of lower crust from southwest to northeast iii the location and timing of the subsequent extension in metamorphic core complexes and the rio grande rift and iv the total area eventually involved in basin and range style extension in a broad sense this model has predicted the belt of laramide structures the transport of crust from the coastal region to the continental interior the subsequent extension in metamorphic core complexes and the rio grande rift and the geographic region of late tertiary basin and range extension its principal defects are that i many events are predicted about 5 million to 10 million years too late and ii the wave of crustal thickening does not travel far enough to the east reasonable modifications to the oceanic plate kinematics and rheologies that were assumed may correct these defects the correspondence of model predictions to actual geology is already sufficiently close to show that the hypothesis that horizontal subduction caused the laramide orogeny is probably correct the rocky mountain thrust and reverse faults formed in an environment of east west to northeast southwest compressive stress that was caused by the viscous coupling between the oceanic plate and the base of the north american crust nonuniform crustal thickening by simple shear transport also caused relative uplifts therefore this model is consistent with both of the range forming mechanisms that have been inferred 1 a new proposal that arises from this simulation is that horizontal subduction also caused the subsequent extensional basin and range taphrogeny by stripping away the mantle lithosphere so that the crust was exposed to hot asthenosphere after the oceanic slab dropped away

a test of present day plate geometries for northeast asia and japan alternative geometries for the present day configuration of plate boundaries in northeast asia and japan are tested using nuvel 1 and 256 horizontal earthquake slip vectors from the japan and northern kuril trenches statistical analysis of the slip vectors is used to determine whether the north american eurasian or okhotsk plate overlies the trench along the northern kuril trench sup vectors are well fit by the nuvel 1 pacific north america euler pole but are poorly fit by the pacific eurasia euler pole results for the japan trench are less conclusive but suggest that much of honshu and hokkaido are also part of the north american plate the simplest geometry consistent with the trench slip vectors is a geometry in which the north american plate extends south to 41 n and possibly includes northern honshu and southern hokkaido although these results imply that the diffuse seismicity that connects the lena river delta to sakhalin island and the eastern sea of japan records motion between eurasia and north america onshore geologic and seismic data define an additional belt of seismicity in siberia that cannot be explained with this geometry assuming that these two seismic belts constitute evidence for an okhotsk block two published kinematic models for motion of the okhotsk block are tested the first model which predicts motion of up to 15 mm yr 1 relative to north america is rejected because kuril and japan trench slip vectors are fit more poorly than for the simpler geometry described above the second model gives a good fit to the trench slip vectors but only if okhotsk north america motion is slower than 5 mm yr 1

trenching studies of active faults in kamchatka eastern russia palaeoseismic tectonic and hazard implications the central part of the kamchatka peninsula is characterized by a well defined depression associated with active volcanism aligned ne sw on the east the depression is bounded by a prominent system of active faults known as the east kamchatka fault zone ekfz in order to improve understanding of the behaviour and kinematic role of this fault zone a fieldwork programme including study of trenches was conducted in the north central part of this system aerial photograph analysis ground truthed indicates a westward fault dip with predominantly normal slip while lateral offsets of river terraces and stream channels demonstrate a combined dextral component over 20 excavated pits and natural exposures were examined to confirm a detailed tephra succession extending from the early holocene to recent historic eruptions this chronological framework then provided age control on five past faulting events recognised in three trenches these events took place at about 10 5 6 0 4 5 and in a two event succession within a short time span at 3 3 3 2 ka bp event clustering may be characteristic and fault length displacement values suggest earthquakes of m6 5 thus representing a significant new element in regional seismic hazard evaluations additional to events generated at the subduction interface the relatively long gap in faulting since the two most recent events may also be significant for hazard scenarios and there is a possible link between the faulting and volcanic activity in the depression overall the ekfz together with the nachiki transverse zone farther south is thought to define a regional scale block that is extending eastwards independently from the rest of kamchatka

continuum calculations of continental deformation in transcurrent environments a thin viscous sheet approximation is used to investigate continental deformation near a strike slip boundary the vertically averaged velocity field is calculated for a medium characterized by a power law rheology with stress exponent n driving stresses include those applied along boundaries of the sheet and those arising from buoyancy forces related to lateral differences in crustal thickness calculations are performed with either velocity or stress prescribed on the transcurrent boundary exact and approximate analytic solutions for a region with a sinusoidal strike slip boundary condition are compared with solutions for more geologically relevant boundary conditions obtained using a finite element technique the across strike length scale of the deformation is approximately 1 4 n 1 2 times the dominant wavelength of the imposed strike slip boundary condition for both the analytic and the numerical solutions this result is consistent with length scales observed in continental regions of large scale transcurrent faulting an approximate linear relationship between displacement and rotation is found that depends only on the deformation length scale and the rheology calculated displacements finite rotations and distribution of crustal thicknesses are consistent with those observed in the region of the pacific north america plate boundary in california