downward migration of seismic activity prior to some great shallow earthquakes in japanese subduction zones a possible intermediateterm precursor before the 1944 tonankai earthquake along the nankai trough seismic activity increased in the shallow depths and then the activity gradually migrated downwards when it reached its limit a depth of approximatelty 70 km the main shock occurred several deep earthquakes including one ofm5 3 occurred several months prior to the tonankai earthquake a similar downward migration pattern also can be recognized regarding the 1952 tokachi oki earthquake in this case the deepest earthquakes reached about 400 km this may be one of the intermediate term precursory phenomena of great thrusttype earthquakes in subduction zones recent observations in the tokai district along the suruga trough where a large earthquake is expected to occur in the future suggest a similar downward migration pattern in the land area

pseudotachylytes from corsica fossil earthquakes from a subduction complex subduction zones are the most seismically active tectonic environment on earth here we report fault rock evidence of palaeoearthquakes from a subduction complex pseudotachylyte veins formed by shear heating on faults at seismic strain rates contain dendrites spherulites and acicular minerals characteristic of blueschist facies conditions the veins have been found in both mantle peridotite and gabbro our discovery from corsica opens a new avenue of research towards understanding mechanisms related to subduction zone earthquakes and the conditions at which these are released

seismic evidence for overpressured subducted oceanic crust and megathrust fault sealing water and hydrous minerals play a key role in geodynamic processes at subduction zones however there are few seismological constraints on the presence of water within the subducting plate audet et al model converted teleseismic waves to constrain the seismological properties of subducted oceanic crust from the cascadia continental margin to its intersection with the forearc mantle their observations suggest that water is pervasively present in fluid form at high pore pressures indicating that the megathrust is a low permeability boundary these results may hold important implications for our understanding of seismogenesis subduction zone structure and the mechanism of episodic tremor and slip this study models converted teleseismic waves to constrain the seismological properties of subducted oceanic crust from the cascadia continental margin to its intersection with the forearc mantle the observations suggest that water is pervasively present in fluid form at high pore pressures indicating that the megathrust is a low permeability boundary these results may hold important implications for our understanding of seismogenesis subduction zone structure and the mechanism of episodic tremor and slip water and hydrous minerals play a key part in geodynamic processes at subduction zones1 2 3 by weakening the plate boundary aiding slip and permitting subduction and indeed plate tectonics to occur4 the seismological signature of water within the forearc mantle wedge is evident in anomalies with low seismic shear velocity marking serpentinization5 6 7 however seismological observations bearing on the presence of water within the subducting plate itself are less well documented here we use converted teleseismic waves to obtain observations of anomalously high poisson s ratios within the subducted oceanic crust from the cascadia continental margin to its intersection with forearc mantle on the basis of pressure temperature and compositional considerations the elevated poisson s ratios indicate that water is pervasively present in fluid form at pore pressures near lithostatic values combined with observations of a strong negative velocity contrast at the top of the oceanic crust our results imply that the megathrust is a low permeability boundary the transition from a low to high permeability plate interface downdip into the mantle wedge is explained by hydrofracturing of the seal by volume changes across the interface caused by the onset of crustal eclogitization and mantle serpentinization these results may have important implications for our understanding of seismogenesis subduction zone structure and the mechanism of episodic tremor and slip

a subducting seamount beneath the nankai accretionary prism off shikoku southwestern japan a multi channel seismic mcs reflection survey was conducted to study the structure of the nankai convergent margin off shikoku island in july 1997 based on reflection characteristics we could identify three major seismic reflection units i e units a b and c the mcs data as well as swath bathymetric data reveal a buried circular seamount subducting beneath the nankai accretionary prism the subducting seamount is responsible for deformation of the accretionary wedge resulting in a compressed uplifted sediment knoll and many steep escarpments with ne sw strike a possible thrust fault was identified on the seaward flank of the seamount indicating compressive deformation a tectonic model for subduction of the seamount is presented based on the interpretation of the mcs data

shear zone thickness and the seismicity of chilean and marianas type subduction zones chilean type convergent margins have many large m 7 6 earthquakes whereas marianas type ones do not this dichotomy is enigmatic if the plate interface is viewed as a thin frictional decollement whereas it becomes understandable if it is viewed as a relatively thick sediment filled shear zone which thins or thickens arcward depending on subduction speed and sediment supply chilean type margins have thick trench fills and their shear zones generally thin arcward from inlets as much as several thousand metres high the most pronounced thinning being located near backstops tall up to several kilometres seamounts are subducted essentially intact to relatively great depths and confining pressures before jamming into the roof of the channel and becoming seismogenic asperities their near basal ruptures can generate large thrust type earthquakes mainly concentrated in seismic fronts near backstops marianas type margins in contrast have thin trench fills and their shear zones generally thicken arcward from inlets that can be as little as 300 m high seamounts are truncated near the inlet at low confining pressures and generate only small earthquakes after passing the inlet they do not touch the roof and therefore cannot generate large earthquakes a similar mechanism may explain seismic gaps at sediment poor regions of subduction zones

seismic images of continental lithosphere seismic reflection profiles in usa and west germany over thirty years ago showed strong reflections from the lower crust but technical advances through the 1970s were needed for the method to image the continental lithosphere the british institutions reflection profiling syndicate birps found advantage in working offshore and has obtained over 12 000 km of profiles around britain during the past decade giving the densest regional coverage anywhere reflections are observed from strata in sedimentary basins from low angle faults that can be followed from surface to mid crust from a laterally variable reflective zone in the lower crust from the moho and from fault like features in the upper mantle there is evidence around britain that faults bounding sedimentary basins are older structures which have been reactivated it is argued that lower crustal reflectivity around britain is mainly due to anastomosing ductile shears which accommodate extension in the lower crust commensurate with extension in the upper crust by faulting and block rotation there is also evidence of reflections preserved from older orogenic features so that reflections interpreted as extensional ductile shears in the lower crust may also derive from older features that have been reactivated most reflections in the lower crust arise from rocks that have been highly strained evidence of lateral variability in the reflective character of the moho suggests that the moho has been actively involved in the structural evolution of the lithosphere upper mantle reflections are interpreted as shear zones which transfer crustal deformation through the lithosphere

upper plate tectonic stress state may influence interseismic coupling on subduction megathrusts we document a correlation between along strike variation in interseismic coupling depth on subduction interfaces and the upper plate tectonic stress state in new zealand vanuatu and southwest japan deep interseismic coupling occurs where the upper plate stress regime is contractional to transpressional whereas a shallowing of interseismic coupling occurs where there is an along strike shift to back arc or intra arc extension to explain this relationship we draw on theoretical studies suggesting that the fl uid pressure state within the upper plate and on the subduction interface has a strong control on the depth of the transition from frictional to viscous behavior lower fl uid pressures e g close to hydrostatic are expected where the over riding plate is undergoing tectonic extension whereas higher fl uid pressures e g close to lithostatic are expected where the over riding plate experiences long term e g 10 5 yr tectonic shortening low fl uid pressures within the upper plate may lead to a shallow frictional to viscous transition compared to an upper plate that is highly overpressured we hypothesize that the state of tectonic stress and structural permeability in the upper plate are yet other variables to consider when evaluating which physical mechanisms control interseismic coupling of subduction megathrusts

heterogeneity spectrum and scale anisotropy in the upper crust revealed by the german continental deep drilling ktb holes crustal heterogeneities and their statistical characteristics bear important information about the dynamic processes and evolution of the crust the velocity well log data from the german continental deep drilling project ktb offer a rare opportunity to measure directly the properties of crustal heterogeneities in this study we first estimated the power spectrum of crustal heterogeneities from the p velocity well logs of the two holes for the first time a power law spectrum of crustal heterogeneities was directly observed in the spatial wavelength range from a few meters to a few kilometers the slope of the vertical 1d power spectrum of the random velocity perturbations is about 1 1 corresponding to a flicker noise random medium the scale anisotropy ratio aspect ratio and the horizontal spectral slope of crustal heterogeneities were also estimated by cross correlation analysis between these two holes to be about 1 8 and 2 0 respectively

structure and seismic response of extended continental crust stochastic analysis of the strona ceneri and ivrea zones italy a two dimensional stochastic seismic model of the crustal cross section exposed in the strona ceneri and ivrea zones northern italy was derived from geologic maps and petrophysical data with techniques used in microfabric analysis the model has a layered self affine structure which agrees with the intrinsic scale independence of geologic phenomena synthetic seismograms from this model replicate many pertinent features observed in deep seismic reflection data and allow insight into crustal structure 1 reflectivity depends on the lithological populations present and their distribution within a crustal unit 2 transparent zones within the reflective lower crust can arise from metamorphic reactions 3 the transition from upper to lower crust may be a change in small scale seismic fabric rather than a first order seismic discontinuity

the frictional hydrologic metamorphic and thermal habitat of shallow slow earthquakes slow earthquakes have recently been discovered on the shallow parts of some subduction zones a review of the conditions under which these quakes form reveals that they could be common in most subduction zones globally

synthetic seismograms for a complex crustal model the algorithm of the original reflectivity method has been vectorized and implemented on a cdc cyber 205 computer calculation times are shortened by a factor of 20 to 30 compared with a general purpose computer with a capacity of several million floating point operations per second mflop the rapid calculation of synthetic seismograms for complex models high frequency sources and all offset ranges is a provision for modeling not only particular phases but the whole observed wavefield as an example we model refraction data of the black forest southwest germany and are able to derive rather tight constraints on the physical properties of the lower crust

crustal deformation dynamics and stress evolution during seamount subduction high resolution 3 d numerical modeling funding was received from the people programme marie curie actions of the european union s seventh framework programme fp7 2007 2013 under rea grant agreement 604713 zip zooming in between plates v s and c r r acknowledge financial support received through the hades project ctm2011 30400 c01 and ctm2011 30400 c02 funded by spanish mineco

deep seismic reflection profiling and the growth of the continental crust abstract deep seismic reflection profiling provides evidence for the style and magnitude of phanerozoic crustal growth by both lateral and vertical accretion profiles across convergent zones have been interpreted as showing tectonic underplating of the continental crust by fragmented oceanic lithosphere profiling across extensional zones provides evidence for magmatic underplating of the crust in addition to studies over active areas deep reflection profiling over large areas of western europe north america and australia shows that the lower crust is strongly reflective and seismically distinct from the upper crust and upper mantle the layering is detected beneath many different terranes of varying crustal type and basement age from archaean to cenozoic and is often apparently younger than the tectonic or metamorphic age of the exposed upper crust thus continental growth curves estimated from upper crustal measurements alone may underestimate phanerozoic crustal additions and mean crustal ages similarly estimated may be too high

wide angle seismic reflections from two dimensional random target zones we have analyzed the shot gather seismic response of a buried zone of isotropic random velocity fluctuations scattering analysis predicts that events at offsets greater than the depth of the target zone will appear more continuous and coherent than events at zero and small offsets due to an effective dip filter inherent to the source target receiver geometry the scattering theory is confirmed with finite difference synthetic seismograms calculated for an idealized crustal model the closely spaced synthetic seismograms show bright events having the appearance of refraction arrivals which suggest a horizontally layered structure in the target zone these effects appear to be independent of the magnitude of the velocity fluctuations provided that the magnitude of the fluctuations is small 3 0 of the mean velocity we have made a second finite difference calculation of a crustal scale model to reinterpret previously published wide angle data from the black forest germany sandmeier and wenzel 1986 the response of the lower crust seen in the conventionally recorded field data can be explained by a lower crustal model of isotropic small amplitude random velocity fluctuations we present this in contrast to a laterally homogeneous model of a finely layered lower crust of alternating high and low velocities proposed by sandmeier and wenzel 1986 the synthetic studies and the scattering analysis demonstrate that the geometry of the wide angle experiment produces signals which appear to be from layered target zones even where no layering exists they also underscore the importance of conducting seismic surveys which make use of closely spaced receivers and which integrate vertical incidence and wide angle recording

small scale heterogeneity and large scale velocity structure of the continental crust since conrad 1925 first identified a midcrustal refraction event in central europe seismologists have developed layered models of the continental crust with velocity steps at various levels modern analysis of crustal refraction data including automated inversion relies on ray theoretical techniques or wave theory for plane layered media which generally parameterize the earth as a series of grossly horizontal layers the bias toward layered models is in part based on physical and geologic intuition but it is to a large measure a result of available interpretation methods for example layer based reflectivity and ray tracing methods in contrast reflection seismology has produced a picture of a highly heterogeneous crust with wavelength scale variations in velocity existing at different levels of the crust in different tectonic regimes finite difference modeling shows that zones of small amplitude random velocity fluctuations not only produce short discontinuous reflection segments as observed on deep seismic reflection data but also produce strong wide angle reflections some of the events identified as crustal wide angle reflections in field data may in fact be the result of scattering from wavelength scale heterogeneities and need not result from reflection from a first or second order velocity discontinuity significant changes in the large scale velocity structure associated with zones of random velocity fluctuations produce relatively small changes in the dynamic properties of the backscattered wave field this observation is in disagreement with the common credo that the seismic refraction method is primarily sensitive to changes in the large scale velocity structure as a consequence the ray theoretical interpretation of such wide angle reflections as first order discontinuities leads to crustal velocity structures that are nonunique and may be erroneous in detail while correctly predicting average crustal velocity and thickness we are not suggesting that crustal velocity does not change with depth on the large scale rather we are proposing that crustal seismic events observed at intermediate and large offsets may be the result of scattering from wavelength scale velocity fabric rather than specular reflection or refraction from first order discontinuities

zagros blueschists episodic underplating and long lived cooling of a subduction zone pressure temperature time p t tp t t trajectories of high pressure rocks provide important constraints to understand the tectonic evolution of convergent margins new field observations and p t tp t t constraints for the evolution of the only known blueschist facies exposure along the zagros suture zone in southern iran are reported here these blueschists now exposed in tectonic windows under the sanandaj sirjan zone upper plate constitute accreted fragments of the tethyan domain during n directed cretaceous subduction three units were identified in the field from top to bottom the ashin unit mafic and felsic gneisses the seghin complex mafic tuffs and ultramafics and the siah kuh massif coherent volcanic edifice microstructural observations p tp t estimates and rb sr deformation ages indicate that the ashin unit possibly underwent burial down to 30 35 km and 550 c along a relatively warm p tp t gradient c 17 km and was ultimately deformed between 85 and 100 ma the seghin complex exhibits remarkably well preserved hp lt assemblages comprising lawsonite glaucophane aragonite omphacite and garnet p t tp t t reconstruction indicates that this slice was subducted down to c 50 km at temperatures of c 500 c along a very cold subduction gradient c 7 km deformation in the seghin complex stopped at around 65 ma close to peak metamorphic conditions field relationships and estimates of the p tp t trajectory followed by the siah kuh volcanic edifice indicate that this massif was lately subducted down to 15 km depth along the same very cold gradient this slice stack represents a well preserved field example i highlighting the existence of transient underplating processes juxtaposing pluri kilometric tectonic slices along the subduction channel and ii imaging the discontinuous down stepping of the active main subduction thrust with ongoing accretion the zagros blueschists also record an apparent cooling of the zagros subduction zone between 90 and 65 ma most likely as a thermal response of a geodynamic perturbation toward the end of the cretaceous

seismic properties of rocks from an exposure of extended continental crust new laboratory measurements from the ivrea zone abstract compressional wave velocities for new samples from the ivrea verbano zone ivz and strona ceneri zone scz have been measured to 600 mpa in order to better understand the petrologic controls on velocity and to provide an expanded base for seismological models of the crust as represented by these terrains these data coupled with previous measurements show that preferred orientation of mica controls both the velocity hysteresis with pressure and seismic anisotropy of scz and ivz schists and gneisses ivz paragneisses form the basis of a linear velocity density relationship for metagraywackes of various metamorphic grade this relationship falls close to that expected for mafic rocks and overlaps velocity density data for ivz mafic gneisses one dimensional near vertical reflection models based on these data and recent structural reconstructions of the terrain show that lower crustal reflectivity can be related to high reflection coefficients associated with interlayered sequences of upper amphibolite granulite facies paragneiss mafic gneiss and ultramafic rocks transparent upper crustal levels in the model result from the overall similarity of velocities for schists gneisses and granites from the structurally and lithologically complex scz a mean lower crustal velocity of 7 1 km s is calculated for the part of the ivz composed of interlayered mafic gneiss and paragneiss and for another area dominated by a mafic layered igneous complex the velocity density analysis coupled with this observation emphasize the ambiguity inherent in interpretation of crustal composition from p wave velocity data in the absence of other constraints

stochastic characterization and seismic response of upper and middle crustal rocks based on the lewisian gneiss complex scotland summary we have developed statistical models of upper and middle crustal seismic velocity heterogeneity based on geologic maps and petrophysical data from the lewisian gneiss complex scotland the level of heterogeneity we have measured is relevant to seismic exploration of the crystalline crust using conventional reflection and refraction techniques we digitized two 1 10560 geologic maps of laxfordian archean age lewisian rocks on a 26 8m grid both maps are believed to be representative of the upper and middle crust in north western scotland and both are believed to provide cross sectional views of parts of the crust the digital maps were characterized by the statistics of their lithologic populations and by their 2 d spatial autocorrelation functions different lithologies were assigned seismic velocities appropriate for the mid crust using petrophysical data three lithologies are dominant silicic gneisses vp 6 2 km s 1 mafic amphibolites vp 6 75 kms 1 and intermediate composition schists vp 6 5kms 1 both maps have self affine spatial fabrics the first map covers the core of a syncline its autocorrelation function defines a medium with a fractal dimension of 2 78 a horizontal characteristic length of about 244m and a vertical correlation of about 133m aspect ratio is 2 1 it has an essentially trimodal velocity lithology population consisting of 37 per cent silicicgneiss 43 per cent mafic amphibolites and 20 per cent schists this map is representative of 30 40 per cent of laxfordian rocks the second map is a plan view which can be rotated 90 about an axis perpendicular to strike to give a cross section this map is characterized by a fractal dimension of 2 55 with a horizontal correlation length of about 111m and a vertical correlation of about 38m aspect ratio 3 1 it has a nearly bimodal population consisting of 77 per cent silicic gneisses and 22 per cent mafic amphibotites it is representative of 60 70 per cent of laxfordian rocks lastly we examine the seismic response of an upper to middle crust based on our statistical models using acoustic and elastic 2 d finite difference synthetic seismograms short offset shot records demonstrate that a lewisian upper crust produces scattered waves which significantly disrupt signals reflecting from deeper levels measurements of transmission scattering q and coda decay rates confirm that seismic scattering in lewisian type crust is strong the migrated cmp response of a lewisian crustal model shows the characteristic salt and pepper pattern often observed in the upper crust and described incorrectly as transparent we suggest that translucent is a more appropriate descriptor

gravity interpretation of a unified 2 d acoustic image of the central alpine collision zone summary four individual deep seismic reflection profiles across the eastern and southern swiss alps have been combined along the trend of the bouguer gravity anomalies onto the european geotraverse the profiles were migrated ray theoretically with the migration velocity model derived from a smoothed reinterpretation of the seismic wide angle profiles running parallel to the strike of the alpine arc the contours of this velocity model were ray theoretically depth migrated together with the combined digitized line drawings of the seismic reflection profiles the resulting acoustic image depicts the subduction of parts of the lower crust of the european plate beneath the adriatic promontory of the african plate at a low angle of 15 based on the interpretation of the seismic data orogenic crustal thickening is attributed to the stacking of crystalline nappes onto the upper crust of the european plate and wedging of the european and adriatic middle and lower crusts the latter being ill constrained by the seismic evidence alone the south vergent thrusting of the southern alps can be accounted for by the observed downbending of the adriatic moho and the lower crust in conjunction with the inferred wedging at mid crustal levels using the geometric constraints provided by seismic data gravity modelling of the alpine lithosphere asthenosphere system relative to stable central european platform clearly favours a gently inclined subduction zone reaching down to at least 200 km depth such a gently dipping subduction zone is at odds with prominent models on alpine geodynamics which favour a near vertical orientation of the subducting lithosphere verschluckung in agreement with the wedging hypothesis indirectly inferred from the seismic data the short wavelength part of the alpine gravity anomaly requires a middle crust of anomalously high density in the axial zone of the orogen both seismic and gravimetric evidence therefore suggest that late orogenic lithospheric shortening and crustal thickening was governed by the mechanical decoupling of the upper middle and lower parts of the crust the amount of subducted lower crustal material cannot be constrained by gravity modelling since its gravity effect cannot be separated from the effect of the subducting slab

seismic structure of gneissic granitic upper crust geological and petrophysical evidence from the strona ceneri zone northern italy and implications for crustal seismic exploration summary we present a stochastic seismic model of gneissic granitic upper crust by estimating the amplitude probability distribution and the autocorrelation function of the inferred seismic velocity field of the strona ceneri zone northern italy unlike previous investigations of the adjacent ivrea zone and the lewisian gneiss complex scotland the lithologies mapped in the strona ceneri zone cannot be directly correlated with laboratory measured seismic velocities the amplitude probability distribution of the extensive petrophysical database and the relatively uniform mineralogical composition of the lithologies indicate that the velocity distribution in the strona ceneri zone can be approximated as being continuous and gaussian based on this assumption and using the spatial autocorrelation function measured from a digitized geological cross section we have modelled the fine scale structure of the strona ceneri zone as a self affine medium with a fractal dimension of 2 50 0 15 and horizontal and vertical characteristic scales of 2 6 0 4 km and 1 4 0 3 km respectively on near offset shot gather seismic data strona ceneri type upper crustal heterogeneity produces a reflected field with the appearance of incoherent scattered noise characteristic of non reflective transparent upper crust found in extended phanerozoic provinces this type of heterogeneity has a considerable effect on seismic imaging of underlying reflective lower crust individual reflection segments are disrupted and shifted in time by up to several tenths of a second and the boundaries of the reflective zone become distorted and blurred at far offsets seismic scattering from heterogeneity of this type reproduces the commonly observed rapid decay and scattered distribution of pg amplitudes comparable amplitude decay patterns result from layered deterministic upper crustal velocity structures with velocity inversions

a stochastic view of lower crustal fabric based on evidence from the ivrea zone despite its complicated history the ivrea zone is considered to be a representative surface exposure of extended continental crust we have digitized two standard 1 25 000 geological maps from this area and evaluated their structural statistics because of the subvertical orientation of the ivrea zone these maps can be considered as small scale cross sections through the lower continental crust the autocorrelation functions of the digitized lithologies measured from these maps show a clear self similar or fractal rather than a gaussian or deterministic trend we found that an anisotropic von karman correlation function with an aspect ratio around 4 and a hurst number of 0 3 corresponding to a fractal dimension of 2 7 matches the observed data our results represent an explicit confirmation of previous indirect evidence for the fractal nature of lithospheric heterogeneities and provide the means to construct realistic crustal scale seismic models of ivrea type lower continental crust

migration of synthetic seismograms for crustal structures with random heterogeneities summary fast frequency wavenumber migration is used to migrate synthetic seismograms calculated by a hybrid method for crustal models with complex 2 d scattering structures in the lower crust scatterers with predominantly horizontal or vertical orientation and scatterers without preferred orientation isotropic scatterers are investigated horizontal scatterers or lamellae have frequently been suggested for the lower crust on the basis of modern reflection seismological experiments in none of the cases studied here are the scattering structures imaged correctly the reason is mainly multiple scattering which can produce coherent arrivals in the seismogram sections and hence coherent signatures in the migrated sections which have no relation to structure imaging is generally acceptable for horizontal scatterers but for isotropic or vertical scatterers the migrated signatures are also horizontal and thus do not represent reality it is concluded that the highly detailed line drawings which are popular in crustal reflection seismology are less reliable than believed as far as the internal structure of scattering zones and the scatterer orientations are concerned horizontal or subhorizontal structures in the lower crust may be less common than assumed the paper also briefly addresses a few methodical aspects of fast frequency wavenumber migration for depth dependent background structure in particular the concept of sounding beams their width and the construction of a migrated section from several sounding beams

fractals chaos self organized criticality and tectonics some geological structures have simple geometrical forms and can be analysed using deterministic equations examples include alluvial fans and many sedimentary basins but most geological structures are complex and appear to defy mathematical analyses yet in the complexity there is an order complex geological structures generally obey fractal statistics examples include topography distributions of earthquakes and faults and mineral deposits an unresolved question is whether the fractal order is simply the result of scale invariance or the result of governing equations that yield deterministic chaos in order to try to answer this question a variety of slider block models have been considered the stick slip behaviour of slider block models is a simple analogy to earthquakes a pair of slider blocks has been shown to behave chaotically models that use many slider blocks exhibit self organized criticality and generate fractal statistics similar to the statistics of regional seismicity

accretionary history and crustal evolution of the variscan belt in western europe abstract the variscan belt of western europe is part of a large intra paleozoic belt extending on both sides of the atlantic from the ouachitas in the us and the mauritanides in west africa to the bohemian massif in czechoslovakia and poland the belt was constructed between 500 and 250 ma as a result of convergence and collision between two main continents laurentia baltica and africa after the closure of various oceanic basins iapetus rheic and galicia massif central now recorded as discrete in part cryptic sutures which form the roots of large crystalline nappes in which are found dismembered remnants of oceanic crust and mantle the accretionary history of the variscan belt in western europe is related to the progressive closure of probably two oceanic basins the rheic in the north and another the galicia massif central in the south by respectively southward and northward intraoceanic subduction followed by obduction and intracontinental lithospheric subduction the result is a broad 700 800 km fan like divergent belt which must be considered as a classical obduction collision orogen with a long intracontinental history of about 100 150 ma its arcuate shape is probably due to the northwestward impingement of a promontory of the african continent into laurentia baltica

the subduction plate interface rock record and mechanical coupling from long to short timescales abstract short and long term processes at or close to the subduction plate interface e g mineral transformations fluid release seismicity and more generally deformation might be more closely related than previously thought increasing evidence from the fossil rock record suggests that some episodes of their long geological evolution match or are close to timescales of the seismic cycle this contribution uses rocks recovered episodically from subduction zones together with insights from thermomechanical modelling to provide a new dynamic vision of the nature structure and properties of the plate interface and to bridge the gap between the mechanical behavior of active subduction zones e g coupling inferred from geophysical monitoring and fossil ones e g coupling required to detach and recover subducted slab fragments based on critical observations and an exhaustive compilation of worldwide subducted oceanic units for which the presence near the plate interface rock types pressure temperature t p gradients thickness and timing of detachment can be assessed the present study demonstrates how long term mechanical coupling exerts a key control on detachment from the slab and potential rock recovery critical assessment of rock t p characteristics indicates that these fragments can indeed be used as natural probes and provide reliable information on subduction interface dynamics down to 2 8 gpa rock clusters are identified at depths of 30 55 60 and 80 km with some differences between rock types data also reveal a first order evolution with subduction cooling in the first 5 myr which is interpreted as reflecting a systematic trend from strong to weak mechanical coupling after which subduction is lubricated and mostly inhibits rock recovery this contribution places bounds on the plate interface constitution regular thickness

subducted seamount imaged in the rupture zone of the 1946 nankaido earthquake the nankai trough is a vigorous subduction zone where large earthquakes have been recorded since the seventh century with a recurrence time of 100 to 200 years the 1946 nankaido earthquake was unusual with a rupture zone estimated from long period geodetic data that was more than twice as large as that derived from shorter period seismic data in the center of this earthquake rupture zone we used densely deployed ocean bottom seismographs to detect a subducted seamount 13 kilometers thick by 50 kilometers wide at a depth of 10 kilometers we propose that this seamount might work as a barrier inhibiting brittle seismogenic rupture

tectonic slicing of subducting oceanic crust along plate interfaces numerical modeling multikilometer sized slivers of high pressure low temperature metamorphic oceanic crust and mantle are observed in many mountain belts these blueschist and eclogite units were detached from the descending plate during subduction large scale thermo mechanical numerical models based on finite difference marker in cell staggered grid technique are implemented to investigate slicing processes that lead to the detachment of oceanic slivers and their exhumation before the onset of the continental collision phase in particular we investigate the role of the serpentinized subcrustal slab mantle in the mechanisms of shallow and deep crustal slicing results show that spatially homogeneous serpentinization of the sub moho slab mantle leads to complete accretion of oceanic crust within the accretionary wedge spatially discontinuous serpentinization of the slab mantle in form of unconnected patches can lead to shallow slicing of the oceanic crust below the accretionary wedge and to its deep slicing at mantle depths depending on the patch length slab angle convergence velocity and continental geothermal gradient p t paths obtained in this study are compared to natural examples of shallow slicing of the crescent terrane below vancouver island and deeply sliced crust of the lago superiore and saas zermatt units in the western alps

interseismic coupling megathrust earthquakes and seismic swarms along the chilean subduction zone 38 18 s the recent expansion of dense gps networks over plate boundaries allows for remarkably precise mapping of interseismic coupling along active faults the interseismic coupling coefficient is related to the ratio between slipping velocity on the fault during the interseismic period and the long term plates velocity but the interpretation of coupling in terms of mechanical behavior of the fault is still unclear here we investigate the link between coupling and seismicity over the chilean subduction zone that ruptured three times in the last 5 years with major earthquakes maule mw 8 8 in 2010 iquique mw 8 1 in 2014 and illapel mw 8 4 in 2015 we combine recent gps data acquired over the margin 38 18 s with older data to get the first nearly continuous picture of the interseismic coupling variations on the subduction interface here we show that at least six low coupling zones lcz areas where coupling is low relatively to the neighboring highly coupled segments can be identified we also find that for the three most recent mw 8 events co seismic asperities correlate well with highly coupled segments while lczs behaved as barriers and stopped the ruptures the relation between coupling and background seismicity in the interseismic period before the events is less clear however we note that swarm sequences are prone to occur in intermediate coupling areas at the transition between lcz and neighboring segments and that the background seismicity tends to concentrate on the downdip part of the seismogenic locked zone thus highly coupled segments usually exhibit low background seismicity in this overall context the metropolitan segment that partly ruptured during the 2015 illapel earthquake appears as an outlier since both coupling and background seismicity were high before the rupture raising the issue of the remaining seismic hazard in this very densely populated area

zagros orogeny a subduction dominated process this paper presents a synthetic view of the geodynamic evolution of the zagros orogen within the frame of the arabia eurasia collision the zagros orogen and the iranian plateau preserve a record of the long standing convergence history between eurasia and arabia across the neo tethys from subduction obduction processes to present day collision from 150 to 0 ma we herein combine the results obtained on several geodynamic issues namely the location of the oceanic suture zone the age of oceanic closure and collision the magmatic and geochemical evolution of the eurasian upper plate during convergence as testified by the successive sanandaj sirjan kermanshah and urumieh dokhtar magmatic arcs the p t t history of the few zagros blueschists the convergence characteristics across the neo tethys kinematic velocities tomographic constraints subduction zones and obduction processes together with a survey of recent results gathered by others we provide lithospheric scale reconstructions of the zagros orogen from 150 to 0 ma across two sw ne transects the evolution of the zagros orogen is also compared to those of the nearby turkish and himalayan orogens in our geotectonic scenario for the zagros convergence we outline three main periods regimes 1 the mid to late cretaceous 115 85 ma corresponds to a distinctive period of perturbation of subduction processes and interplate mechanical coupling marked by blueschist exhumation and upper plate fragmentation 2 the paleocene eocene 60 40 ma witnesses slab break off major shifts in arc magmatism and distributed extension within the upper plate and 3 from the oligocene onwards 30 0 ma collision develops with a progressive sw migration of deformation and topographic build up sanandaj sirjan zone 20 15 ma high zagros 12 8 ma simply folded belt 5 0 ma and with partial slab tear at depths 10 ma to present our reconstructions underline the key role played by subduction throughout the whole convergence history we finally stress that such a long lasting subduction system with changing boundary conditions also makes the zagros orogen an ideal natural laboratory for subduction processes

introduction seismic wave scattering in three dimensionally heterogeneous earth the classical earth model of spherical symmetry or layered structure is undergoing a revolution the earth has been revealed to be laterally heterogeneous everywhere from the crust mantle to the core with scale from the grain size of rocks to the lowest orders of global spherical harmonics figure 1 shows the strength and scale length of heterogeneities in the crust and mantle of the earth in terms of v the perturbation index of seismic p or s wave speed defined as the r m s fractional variations of wave speed over the measured region the scale length of heterogeneities revealed by seismic waves not including the laboratory measurements of rock samples spans 8 orders of magnitude these heterogeneities with different scales have different effects to seismic waves the velocity and density heterogeneities can cause the change in waveform phase or travel time and amplitude fluctuation as well as apparent attenuation of the direct arrivals they can also generate coda waves such as the p coda s coda and lg coda caused by the lithospheric heterogeneities and precursory waves such as the scattered pkp waves as the precursors to pkikp caused by the heterogeneities near the core mantle boundary the near source or near receiver structures can modify the seismic waveforms by resonance and other effect rough topography or rough interface can cause the coupling between body wave and surface wave aligned cracks in the crust can produce the effective anisotrophy a great complexity arises when heterogeneities have interaction with anisotrophy and nonlinearity

the effect of seamount subduction on seismic coupling it has long been speculated that the subduction of seamounts may have profound mechanical consequences particularly on the generation of large subduction zone earthquakes we consider this problem as a special case of the seismic coupling model of scholz and campos the model predicts that the subduction of a large seamount will increase the normal stress across the subduction interface and hence will enhance seismic coupling in the case of coupled arcs this increased coupling will result in a large increase in the recurrence intervals of earthquakes such events will therefore appear rare in the instrumental record as is observed for decoupled arcs the subduction of a large seamount should cause them to become locally coupled three of the latter localities have been found with satellite altimetry data two on the tonga kermadec trench and one on the izu bonin trench the diagnostic feature is a shallowing or obstruction of the trench axis adjacent to a chain or group of large seamounts on the subducting plate in each case the small to moderate plate interface seismicity typical of decoupled seismic zones is almost completely absent in a 100 200 km arc segment indicating that the subduction zone is locally locked as predicted by the model the occurrence of large interplate earthquakes at both tonga sites confirms that they are seismically coupled

slow earthquakes and great earthquakes along the nankai trough abstract we have reexamined reports indicating that slow deformation occurred before the great japan earthquakes of 1944 tonankai and 1946 nankaido and find that the observations are well founded although no quantitative models have previously been proposed to explain all of the relevant data we show that they are satisfied by a simple model for both earthquakes the model based on known properties of subduction zones has slow slip on the subduction interface in an area deeper than the seismic rupture if this model is correct and a similar physical situation holds for an anticipated tokai earthquake existing instruments will be able to reveal the pre slip in real time while differences among the deformation time series at different sites will provide strong constraints on the slow rupture propagation these differences could result in delaying the recognition of a coherent event

seismic attenuation in the vicinity of the geothermal anomaly at urach obtained from near vertical reflection profiles1 the attenuation in the vicinity of the geothermal anomaly at urach was determined by means of two near vertical reflection profiles the attenuation in the sediments and in the upper crust 3 4 km depth was estimated by interpretation of the first refracted arrivals for calculating the attenuation the amplitude decay with respect to distance was used corrections for the spread factor i e the geometric amplitude divergence was deduced from the traveltime curves below the anomaly higher attenuation values q 1 0 008 were observed compared with those in the undisturbed crust q 1 0 002 this effect is probably due to the cracks and fissures in the upper part of the crystalline basement the attenuation in the middle and lower crust was determined using near vertical reflections from this depth interval the use of the spectral ratio method leads to higher values of the effective attenuation q 1eff below the heat flow anomaly compared to those of the normal crust this zone of high q 1eff coincides with the low velocity body below the heat flow anomaly both effects the higher attenuation and the lower velocities could be caused by high temperatures cracks and fissures in the crust

earthquake prediction a critical review summary earthquake prediction research has been conducted for over 100 years with no obvious successes claims of breakthroughs have failed to withstand scrutiny extensive searches have failed to find reliable precursors theoretical work suggests that faulting is a non linear process which is highly sensitive to unmeasurably fine details of the state of the earth in a large volume not just in the immediate vicinity of the hypocentre any small earthquake thus has some probability of cascading into a large event reliable issuing of alarms of imminent large earthquakes appears to be effectively impossible key word earthquake prediction

some attributes of wavefields scattered from ivrea type lower crust abstract the ivrea zone in northern italy is a well studied sliver of extended lower continental crust which was brought to the surface in the course of alpine lithospheric shortening based on two 1 25 000 geological maps from the central ivrea zone and petrophysical parameters of the corresponding lithologies we have constructed geologically constrained stochastic seismic models of ivrea type lower crust a stochastic rather than a deterministic approach is applied due to the limited availability of detailed geological maps and the high degree of small scale structural and petrophysical complexity of the ivrea zone the primary characteristic of the resulting model is its layered self affine or fractal structure and its bimodal velocity distribution while we are not suggesting that ivrea type structure is universal for extended lower continental crust synthetic seismograms may help to constrain some pertinent quantitative aspects observed in deep seismic reflection data 1 1 ivrea type lower crust can explain the igh observed reflection coefficients and bright reflectivity in the lower crust 2 2 observed lower crustal q factors can be largely explained by scattering losses and intrinsic q factors of reflective lower crust are likely to be larger than 1000 which is incompatible with the presence of free volatiles 3 3 in agreement with recent broad band studies of the lower crust ivrea type lower crust does not create any tuning effects and 4 4 for ivrea type heterogeneity the lateral correlation of the reflected signal is largely independent of the lateral characteristic scale of the scattering structures which suggests that multiple scattering may be important in the reflective lower crust

lower crustal reflectivity from waveform inversion summary deep seismic reflection profiles show that the lower part of the continental crust often contains many strong subhorizontal reflections based on forward modelling these reflections are interpreted as a sequence of alternating high and low velocity layers a waveform inversion provides an optimum model by minimizing the misfit between observed data and synthetic seismograms however it is very difficult to invert for the whole data set partly because of the large computation time needed in modelling and partly because the high amplitude reflected energy in the upper crust dominates the misfit function to overcome these difficulties we propose a three step strategy for inverting the deeper part of seismic reflection data over a narrow aperture in which the velocity structure can be assumed to be 1 d first both source wavelet and recorded wavefield are propagated in the plane wave domain downwards to a particular depth for a given 1 d medium because the source receiver distances are short for deep reflections it is assumed that these reflections are at normal incidence it is also assumed that the near surface velocity is known or has been estimated using other inversion schemes the downward propagation algorithm takes account of nearsurface effects such as multiple reflections and attenuation second starting from a known background velocity model a non linear iterative waveform inversion is applied to the propagated data set the method is based on minimizing the difference sample by sample between observed and calculated wavefields in a least squares sense the forward modelling which consists of calculating the data from the model uses a reflectivity method to compute the synthetic seismogram in the plane wave domain the model perturbation is calculated with the help of the conjugate gradients method the non linear inversion gives the short wavelengths of impedance contrast which is a function of density and velocity in the lower crust third an interpretational step is tried in order to recover the medium wavelengths of velocity the algorithm is tested on synthetic data and then applied to a part of the wam data set south west of england a good fit to the real data with the optimum model predicts that the velocity contrast of the layers in the lower crust varies from 0 2 km s to 0 5 km s and the thickness from at least 0 075 km to 1 km

aseismic zone and earthquake segmentation associated with a deep subducted seamount in sumatra the subduction of topographic features has been linked both to earthquake generation and the segmentation of faults seismic imaging reveals a seamount subducted to 40 km depth below sumatra that is associated with an aseismic zone suggesting that at this location the seamount reduces coupling of the slab and overriding plate

subduction megathrust creep governed by pressure solution and frictional viscous flow subduction megathrust slip speeds range from slow creep at plate convergence rates centimetres per year to seismic slip rates metres per second in the largest earthquakes on earth the deformation mechanisms controlling whether fast slip or slow creep occurs however remain unclear here we present evidence that pressure solution creep fluid assisted stress driven mass transfer is an important deformation mechanism in megathrust faults we quantify megathrust strength using a laboratory constrained microphysical model for fault friction involving viscous pressure solution and frictional sliding we find that at plate boundary deformation rates aseismic frictional viscous flow is the preferred deformation mechanism at temperatures above 100 c the model thus predicts aseismic creep at temperatures much cooler than the onset of crystal plasticity unless a boundary condition changes within this model framework earthquakes may nucleate when a local increase in strain rate triggers velocity weakening slip and we speculate that slip area and event magnitude increase with increasing spacing of strong topographically derived irregularities in the subduction interface

variations of fluid pressure within the subducting oceanic crust and slow earthquakes 1 we show fine scale variations of seismic velocities and converted teleseismic waves that reveal the presence of zones of high pressure fluids released by progressive metamorphic dehydration reactions in the subducting philippine sea plate in tokai district japan these zones have a strong correlation with the distribution of slow earthquakes including long term slow slip ltss and low frequency earthquakes lfes overpressured fluids in the ltss region appear to be trapped within the oceanic crust by an impermeable cap rock in the fore arc and impede intraslab earthquakes therein in contrast fluid pressures are reduced in the lfe zone which is deeper than the centroid of the ltss because there fluids are able to infiltrate into the narrow corner of the mantle wedge leading to mantle serpentinization the combination of fluids released from the subducting oceanic crust with heterogeneous fluid transport properties in the hanging wall generates variations of fluid pressures along the downgoing plate boundary which in turn control the occurrence of slow earthquakes

global distribution of seamounts from ship track bathymetry data the distribution of submarine volcanoes or seamounts reflects melting within the earth and how the magma generated ascends through the overlying lithosphere globally 60 latitude we use bathymetry data acquired along 39 5 106 km of ship tracks to find 201 055 probable seamounts an order of magnitude more than previous counts across a wider height range 0 1 1 km implying that 24 000 60 remain to be discovered

lithoprobe new perspectives on crustal evolution lithoprobe is canada s national collaborative multidisciplinary earth science research program directed toward an enhanced understanding of how the north american continent evolved research in i

the snow mountain volcanic complex an on land seamount in the franciscan terrain california the snow mountain volcanic complex in the franciscan terrain of northern california is a 1 1 2 to 2 km thick sequence of altered submarine pillow lavas pillow breccias diabase and interbedded sediments forming a klippe that structurally overlies nearby rocks of the franciscan assemblage the coast range ophiolite and great valley sequence all of the snow mountain rocks bear a weak overprint of blue schist facies metamorphism the volcanic rocks were deposited on top of 200 300 m of interbedded turbiditic greywacke shale and radiolarian chert which are interpreted as an interfingering abyssal fan and pelagic sequence the volcanic pile is therefore supracrustal and is not an ophiolitic fragment formed at a mid ocean spreading center rather it is probably the remnant of a former seamount the lavas are mostly titanium rich transitional basalts minor associated rocks include rhyolite and peralkaline rhyolite relict clinopyroxenes in the basaltic rocks are all titanium aluminum rich augites thes

sensitivity of the lateral correlation function in deep seismic reflection data in the absence of multiple scattering the lateral correlation function of the reflected wavefield can be related to the lateral correlation function of the medium we have analyzed the lateral correlation function of a series of synthetic seismograms from stochastic models of reflective lower crust with the structural and petrophysical statistics of the ivrea zone our results show that in the presence of ivrea type lower crustal heterogeneity the lateral correlation function of the backscattered wavefield is insensitive to the lateral correlation function of the scattering medium which suggests that multiple scattering is important in the reflective lower crust this interpretation is consistent with lower crustal reflections commonly being shorter than a fresnel zone the absence of clear diffractions and the poor performance of standard migration algorithms on deep seismic reflection data

fossil thermal structure of the southern sanandaj sirjan zone sw iran implications for regional scale tectonics abstract the sanandaj sirjan zone ssz is one of the cimmerian blocks accreted to eurasia during the late triassic and the only tectonic unit of the zagros orogen affected by barrovian metamorphism both the age and absolute temperature of this metamorphic event are still poorly defined in the present study we use the raman spectrometry of carbonaceous material rscm to estimate maximum temperatures attained in the southern ssz along two sw ne cross sections while maximum obtained temperature for precambrian to jurassic rocks ranges from 200 to 600 c all precambrian and early paleozoic samples record temperatures 500 c results thus reveal the existence of a major temperature contrast between precambrian to devonian samples showing tmax 500 c and late devonian to cretaceous ones tmax

crustal structure of the southern alps beneath the intersection with the european geotraverse abstract newly digitized and amplitude controlled record sections from the 1977 southern alps refraction campaign permitted a reinterpretation of the crustal structure in the area between western lombardy and the giudicaria fault the resulting model exhibits considerable lateral heterogeneity in the west below 7 5 km of sediments of the lombardy basin the crust reaches a depth of only 31 km whereas it thickens towards the more mountainous area in the east reaching a depth of 46 km below the adamello massif although the signal character of the corresponding reflections is somewhat erratic the data are satisfied best by models with a low velocity zone in the upper crust an additional small velocity discontinuity from 6 2 to 6 4 km s was found in the middle crust at around 20 km earlier interpretations based on travel times alone included a layer with a velocity of about 7 km s at this depth this high velocity layer was then interpreted as lower crustal material of the adriatic african plate which had been overthrust onto the european plate during the alpine orogeny thus explaining the uplift of the southern alps however this model of crustal doubling is questionable because such a mid crustal high velocity layer is not in agreement with the amplitude data the relatively thin crystalline part of the crust under the lombardy basin is interpreted in accordance with geological evidence as a relic of a late hercynian rifting event

recent seismic activity in the tokai japan region where a large earthquake is expected in the near future abstract a large shallow earthquake is expected to occur in the tokai region of central honshu japan in the near future and various intensive observations are being carried out in the region with the goal of predicting its occurrence so far there have been no reports of any changes which could be regarded as precursors of the expected large earthquake in this area but the following noteworthy changes have been observed in the seismic activity over the past decade 1 1 there has been a spate of large earthquakes of m 6 5 7 in a wide area surrounding the focal region of the expected tokai earthquake 2 2 there has been a significant decline in seismic activity in the area along the suruga trough and the nankai trough which includes the focal region of the anticipated earthquake 3 3 the hypocenters of earthquakes that have occurred in the tokai district have migrated downward from shallow depths since these changes resemble the changes that have appeared before other large shallow earthquakes along the trench in the past it is conceivable that they are long term precursors of a large earthquake future changes in seismic activity in this area demand close attention

link between plate fabric hydration and subduction zone seismicity in alaska subduction carries water into the earth where it can influence seismicity analysis of the structure of the alaskan subduction zone suggests fluid delivery is influenced by faults in the oceanic plate that formed at the mid ocean ridge

seismic refraction investigation of the black forest abstract within the framework of site surveys of the deep continental drilling program ktb deep seismic sounding experiments were carried out in the black forest southwest germany in august 1984 these site surveys comprised seismic refraction as well as seismic reflection measurements this paper will concentrate on the interpretation of the refraction data which were observed along a 240 km long line approximately following the morphological axis of the black forest in a n s direction and extending in the south into the swiss molasse basin from the available seismic refraction data a detailed velocity depth model of the crust and uppermost mantle was derived a clear phase refracted in the basement with a velocity of 5 9 km s is observed for all shots this phase is followed by a reflection from the upper boundary of a low velocity zone this zone is very accentuated under the northern and central black forest and is characterized by a velocity decrease from 6 2 to 5 4 km s between depths of about 8 and 14 km towards the south the velocity inversion gradually loses its intensity under the center of the southern black forest velocities of about 6 km s prevail throughout the upper crust the lower crust which starts at the lower boundary of the low velocity zone shows a very complicated response in the record sections due to long reverberations between travel times of 1 5 3 5 s reduced by 6 0 km s at an epicentral distance of approximately 30 70 km a classical phase correlation is not possible based on a combined interpretation of near angle and wide angle reflection data the lower crust is characterized by a stack of thin layers with alternating high and low velocities the average velocity of the laminated zone is 6 7 km s the reflection from the moho at a depth of 25 26 km is the most prominent phase at epicentral distances greater than 60 km

do subducting seamounts generate or stop large earthquakes seamount subduction is a common process in subduction zone tectonics contradicting a widely held expectation that subducting seamounts generate large earthquakes seamounts subduct largely aseismically producing numerous small earthquakes on rare occasions when they do produce relatively large events the ruptures tend to be complex suggesting multiple rupture patches or faults we explain that the seismogenic behavior of these seamounts is controlled by the development and evolution of an adjacent fracture network during subduction and cannot be described using the frictional behavior of a single fault the complex structure and heterogeneous stresses of this network provide a favorable condition for aseismic creep and small earthquakes but an unfavorable condition for the generation and propagation of large ruptures

the crust as a heterogeneous optical medium or crocodiles in the mist abstract based on petrophysical data geologic maps and a well log we present statistical descriptions of likely upper middle and lower crustal rocks to characterize the fine scale heterogeneity observed in crustal exposures and inferred from deep crustal seismic data the statistical models developed for granitic and metamorphic upper crust and for an extended metamorphic lower crust are used to construct whole crustal models of seismic velocity heterogenity we present finite difference synthetic cmp data from several models which compare favorably with field data the statistical models also permit classification of the seismic reflection experiment and the crustal heterogeneity according to scattering regime the optical or scattering properties of importance for classification are the velocity fluctuation intensity the horizontal and vertical correlation lengths of the medium the correlation function of the medium and the velocity population function for the crustal properties we measured the bandwidth of a typical deep crustal experiment overlaps from the weak to the strong scattering regime with implications for crustal seismic data processing and imaging notably deep crustal signals are likely to have experienced multiple scattering making common seismic imaging techniques of questionable value moreover the details of the unmigrated cmp stacked section bears little resemblance to the underlying medium

the autocorrelation function an image processing tool for fabric analysis abstract the use of the autocorrelation function acf is proposed as a method for fabric and or strain analysis its chief purpose is to extend fabric analysis to fine grained materials such as clays micrites etc or to materials where the determination of individual grain shapes is tedious or impossible the main advantage of the method is that no segmentation of the image prior to analysis is necessary in other words any bias associated with identifying grain boundaries and the like is avoided in as much as the recognition of grains or grain boundaries is a matter of scale the acf can be used without loss of generality whenever the grain size is small in comparison with the size of the analysed area the method is computer based it requires a scanner for the digitization of the images and sufficient storage space and computing power for digital image processing the method uses half tone images as input these can be obtained from thin sections acetate peels polished surfaces photographic negatives or positives etc the acf is a two dimensional function that is related to the average shape that occurs in the fabric it is scaled with respect to the original image and allows simultaneous analysis of small and large features the acf enables a a quantitative description of the fabric in general terms b an efficient grain size determination and c if strain is a valid interpretation of the fabric the derivation of the finite strains in the plane of observation in the first part of the paper the acf method is introduced by way of a few numerically generated fabrics and compared to other methods of fabric and strain analysis in the second part an application to three samples of experimentally deformed solnhofen limestone is demonstrated

strength of stick slip and creeping subduction megathrusts from heat flow observations subduction faults called megathrusts can generate large and hazardous earthquakes the mode of slip and seismicity of a megathrust is controlled by the structural complexity of the fault zone however the relative strength of a megathrust based on the mode of slip is far from clear the fault strength affects surface heat flow by frictional heating during slip we model heat flow data for a number of subduction zones to determine the fault strength we find that smooth megathrusts that produce great earthquakes tend to be weaker and therefore dissipate less heat than geometrically rough megathrusts that slip mainly by creeping

the crustal structure of western europe a thorough understanding of the tectonic processes that formed the con tinental lithosphere under changing geometric and physical conditions through time requires detailed knowledge of the structure of the crust one of the best locations to study this fundamental problem is the northern and western part of europe because it is made up of a number of tectonic provinces ranging in succession from the oldest precambrian areas of fennoscandia to the currently active regions of the mediterranean sea the geology of europe a region extending from the northern tip of scandinavia north cape southward to north africa encompasses the pre 2500 ma archean craton in the northernmost part of the baltic or fennoscandian shield the proterozoic paleozoic and cenozoic pro vinces of northern and central europe that have been added on to this nucleus and the active transition zone between the eurasian and african plates in the mediterranean area cf figure 1 from a tectonic geologic point of view europe can be divided into four main units that cover the whole time span of the history of the continental crust

recent earthquake prediction research in japan japan has experienced many major earthquake disasters in the past early in this century research began that was aimed at predicting the occurrence of earthquakes and in 1965 an earthquake prediction program was started as a national project in 1978 a program for constant monitoring and assessment was formally inaugurated with the goal of forecasting the major earthquake that is expected to occur in the near future in the tokai district of central honshu island the issue of predicting the anticipated tokai earthquake is discussed in this article as well as the results of research on major recent earthquakes in japan the izu earthquakes 1978 and 1980 and the japan sea earthquake 1983

seismic anisotropy of metapelites from the ivrea verbano zone and serie dei laghi northern italy abstract laboratory measurements of compressional and shear wave velocities up to 600 mpa confining pressure were carried out on amphibolite to granulite facies metapelites from the serie dei laghi and ivrea verbano zone northern italy the metapelite samples show moderate to high compressional wave anisotropy and significant shear wave splitting for propagation directions parallel to foliation serie dei laghi samples are generally transversely isotropic for compressional wave propagation a pattern that based on calculated velocities is the consequence of strong preferred orientation of micas ivrea verbano samples generally have an orthorhombic anisotropy pattern for compressional waves that is caused by the preferred orientation of both mica and sillimanite shear wave anisotropy and shear wave splitting patterns differ significantly for the two fabric types intense static recrystallization can obliterate seismic anisotropy but partial static recrystallization does not appreciably reduce anisotropy we conclude that seismic anisotropy may occur in mid to lower continental crust where it is predominantly composed of metapelitic rocks and if there is a consistent regional fabric this seismic anisotropy may be diminished if there is no consistent regional fabric or if postmetamorphic intrusions have caused static recrystallization

stochastic modeling of seafloor morphology inversion of sea beam data for second order statistics at scale lengths less than 100 km or so statistical descriptions of seafloor morphology can be usefully employed to characterize ridge crest processes off ridge tectonics and vulcanism sedimentation and postdepositional transport we seek to develop methods for the estimation of seafloor statistics that take into account the finite precision resolution and sampling obtained by actual echo sounding systems in this initial paper we restrict our attention to the problem of recovering second order statistics from data sets collected by multibeam devices such as sea beam the seafloor is modeled as a stationary zero mean gaussian random field completely specified by its two point covariance function we introduce an anisotropic two point covariance function that has five free parameters describing the amplitude orientation characteristic wave numbers and hausdorff fractal dimension of seafloor topography we formulate the general forward problem relating this model to the statistics of an ideal multibeam echo sounder in particular the along track autocovariance functions of individual beams and the cross covariance functions between beams of arbitrary separation using these second moments as data functionals we then pose the inverse problem of estimating the seafloor parameters from realistic noisy data sets with finite sampling and beamwidth and we solve this inverse problem by an iterative linearized least squares method the inversion method is applied to sea beam transit data from both the pacific and atlantic oceans sea beam system noise stands out as a sharp spike on the along track autocovariance function and can be modeled as a white noise process whose amplitude generally increases with beam angle the five parameters in our second order model can be estimated from the inversion of data sets comprising 100 200 km of track length in general the cross track wave number is the most poorly determined although uncertainties in the assumed sea beam response may bias the values of the fractal dimension using the assumed beamwidth the measured noise values and the seafloor parameters recovered from the inversion we generate sea beam synthetics whose statistical character can be directly compared with raw sea beam data for most of the track segments we have processed thus far the synthetics are similar to the data in the case of one atlantic profile however the comparison clearly indicates the necessity of incorporating higher order statistics the space domain procedures described in this paper can be extended for this purpose

the updip and downdip limits to great subduction earthquakes thermal and structural models of cascadia south alaska sw japan and chile this study examines thermal and structural controls of the updip and downdip rupture limits of great subduction thrust earthquakes data on past great earthquake seismic limits have been compiled for four continental subduction zones cascadia sw japan south alaska and chile these limits have then been compared to the predictions of several models for what constrains great earthquake rupture temperatures on the subduction thrusts have been estimated by finite element numerical models the landward limits of the observed updip aseismic zones correspond to the position where the thrust temperature reaches about 100 c that is depths of about 2 to 10 km for the subduction zones studied this temperature agrees with the dehydration of stable sliding smectite clays to illite chlorite the temperatures in this region are controlled mainly by the thickness of sediment on the incoming crust and by the crustal age and thus heat flow the downdip limits correspond to the depth on the thrust where either 1 the temperature reaches about 350 c which corresponds to thermally activated stable sliding behavior for crustal rocks with a transition to 450 c or 2 about 40 km depth if 350 c is reached at greater depth depths of about 40 km approximately correspond to the intersection of the thrust with the continental forearc moho and this downdip limit may be a consequence of stable sliding serpentinite or talc and other hydrated forearc mantle rocks the primary temperature controls on the downdip region are the age of the subducting oceanic plate and the thrust dip profile secondary control comes from the thickness of incoming sediment the convergence rate and the radioactive heat generation in the overlying forearc the 100 c updip limit occurs near the trench for young subducting plates with a thick sediment section such as cascadia 6 8 ma and up to 80 km landward for older oceanic crust such as south alaska 50 ma the 350 c downdip thermal limit is applicable for young oceanic plates e g cascadia and nankai whereas the forearc mantle limit applies for older plates e g south alaska and chile except near the chile rise for the margins studied that have experienced great earthquakes there is generally good agreement between the postulated thermal and moho limits and the rupture or seismogenic zone as defined by the distribution of aftershocks and by waveform tsunami and dislocation modeling the downdip limit of the interseismic locked zone from dislocation modeling is also in agreement with these limits

features of the water mass front in the northwestern north pacific analysis of the hydrographic data taken in the northern north pacific has revealed a water mass front defined as the abrupt change of mixing rates between the two typical water masses typical water masses are situated above the main thermocline in the subtropical and subpolar gyres this front termed water mass front in the present study can be regarded as the boundary between the subtropical gyre and the subpolar gyre the index for detecting the frontal position from the temperature field is given as the 5 c isotherm at the 300 m level comparisons with data from previous studies are also made and this front appears to correspond to the subarctic front in almost every instance

twelve years of geophysical research in the charlevoix seismic zone abstract various and extensive geophysical monitoring began in 1974 in the charlevoix seismic zone csz the historical record and a major earthquake in 1925 characterize this zone as the likely site of a future damaging earthquake i e m 5 5 or larger and the aim of the monitoring was to develop the capability of predicting such earthquakes the parameters studied included microseismicity seismic travel times electrical impedance vertical movement horizontal movement tilt gravity change and strain through water well level changes although no clear precursor was detected in the months before the largest seismic event that occurred during the observation period m 5 0 the multiparameter experiment provided new insights into the structure and the mechanics of this active region observations of several hundred well located microearthquakes have permitted the delineation of ancient rift faults that are now being activated by high horizontal compression the geographical extent of the zone is sharply defined by the microseismicity which occurs in an area of 35 km 80 km straddling the st lawrence river observation of p waves from accurately timed explosions and s waves from earthquakes support the hypothesis of extensive dilatancy anisotropy eda involving saturated parallel cracks in the crust the data have permitted an estimate of both crack density and the direction for the principal compressive stress in the region n 30 e changes in electrical parameters with periods of months to years were detected at five of the six magnetotelluric stations installed in the region changes in the polarization angle of the earth s electrical telluric field correlated well for seven years with calculated groundwater accumulation levels at one site electric impedance changes showed no clear correlation from station to station but changes at one station appear to correlate with a seismic p wave travel time change anomaly and with the disappearance of a negative gravity anomaly at the time of an unusual number of earthquakes in the range m 2 to 3 from mid 1976 to mid 1977 it is suggested that these correlated changes could be related to aseismic stress redistribution and fluid flow along one of the ne striking rift faults about two years before the m 5 0 earthquake in 1979 continuous monitoring of tilt surface and borehole and well water levels was carried out at a single site charlevoix geophysical observatory although there were unexplained near surface tilts of the order of the earth tide amplitude within two or three days before some earthquakes they were either not spatially coherent or were related to small unexplained temperature changes in the tiltmeter vault no intermediate term one week to one year precursors were seen against a noise level of 1 3 rad rms for near surface tilt observations and a noise level of 0 2 0 3 rad rms for borehole tilt measurements at a depth of 50 100 m the threshold of detectable tidal tilt admittance changes for the borehole tiltmeters was lowered to 2 but no significant changes were seen four coseismic changes in the levels of water in boreholes were observed all were drops in water level which were much larger than expected on the basis of elastic dislocation theory it is suggested that these could be explained either by the collapse of a dilatant zone around the respective hypocentres or by triggered aseisnuc movement on a nearby active fault borehole water level tides observed from 1980 to 1986 exhibit significant variations in amplitude and phase which are thought to be the result of time varying elastic properties in a nearby structure some constraints can now be placed on the long term deformation of the csz as a whole from tiltmeter levelling and horizontal surveys data from a levelling array observed from 1976 to 1986 puts a limit of 0 54 rad yr on the secular tilt rate at the observatory a further limit of 0 1 rad yr is provided by borehole tilt measurements over the period 1983 1986 a limit of 2 cm is placed on vertical movements northwest of the st lawrence river by levelling surveys in 1965 1978 1980 and 1982 no significant 90 confidence deviatoric strain has been detected across the csz from horizontal control surveys carried out in 1912 1919 1965 1978 and 1983 the maximum allowable long term strain rate across the csz from 1919 to 1983 is 0 16 strain yr

weak interplate coupling by seamounts and repeating m 7 earthquakes subducting seamounts are thought to increase the normal stress between subducting and overriding plates however recent seismic surveys and laboratory experiments suggest that interplate coupling is weak a seismic survey in the japan trench shows that a large seamount is being subducted near a region of repeating earthquakes of magnitude m 7 both observed seismicity and the pattern of rupture propagation during the 1982 m 7 0 event imply that interplate coupling was weak over the seamount a large rupture area with small slip occurred in front of the seamount its northern bound could be determined by a trace of multiple subducted seamounts whereas a subducted seamount itself may not define the rupture area its width may be influenced by that of the seamount

thrust type subduction zone earthquakes and seamount asperities a physical model for seismic rupture a thrust type subduction zone earthquake of m sub w 7 6 ruptures an area of approximately 6 000 km sup 2 has a seismic slip of approximately 1 m and is nucleated by the rupture of an asperity approximately 25km across a model for thrust type subduction zone seismicity is proposed in which basaltic seamounts jammed against the base of the overriding plate act as strong asperities that rupture by stick slip faulting a m sub w 7 6 event would correspond to the near basal rupture of a approximately 2 km tall seamount the base of the seamount is surrounded by a low shear strength layer composed of subducting sediment that also deforms between seismic events by distributed strain viscous flow planar faults form in this layer as the seismic rupture propagates out of the seamount at speeds of kilometers per second the faults in the shear zone are disrupted after the event by aseismic slow viscous flow of the subducting sediment layer consequently the extent of fault rupture varies for different earthquakes nucleated at the same seamount asperity because new fault surfaces form in the surrounding subducting sediment layer during each fast seismic rupture

seismic reflection character of the hikurangi subduction interface new zealand in the region of repeated gisborne slow slip events summary we use seismic reflection data to map the geometry and character of the subduction interface in the gisborne area of the hikurangi subduction margin new zealand which experiences repeated shallow 15 km slow slip events the reflection character and geometry in this area is highly variable which we interpret to be related to the subduction of seamounts and underthrust sediments three zones of high amplitude interface reflectivity hrz 1 2 and 3 are interpreted to be the result of fluid rich sediments that have been entrained with subducting seamounts the interface above the hrz zones is shallower than the surrounding areas by 2 4 km due to the warping of the interface to accommodate seamount subduction these zones of high amplitude reflectivity and shallower interface geometry correlate broadly with locations of recorded slow slip events from 2002 to 2008 we hypothesize that effective stresses on the interface may be lower along the northeast margin in areas of high amplitude reflectivity due to 1 the enhanced underthrusting of fluid rich sediment 2 reduced overburden stresses where the interface has been warped to shallower depths to accommodate seamount subduction and 3 potential fluid flow concentration effects leading to overpressure along these shallower interface corrugations from our observations we propose localized reductions in effective stress caused by interface structural relief may be a potential factor in promoting shallow slow slip events

seismic attenuation in the continental crust sw of england summary significant variability in crustal attenuation exists off england s sw coast as shown by birps deep seismic reflection lines the presence of the thick cornubian batholith with no sedimentary cover acting as a low attenuating region compared to normal upper crust with about 300 m of sedimentary cover can account for the observed differences in the strength of signal returning from the lower crust application of the spectral amplitude ratio method to determine the effective quality factor q of the lower crust beneath the granite yields a value of 400 200 a modification of the spectral amplitude ratio method to measure relative q s at different locations indicates that the average q of the upper and middle crust containing the granite is twice that elsewhere although the hypothesis that the lower crustal composition differs significantly beneath the granite cannot be ruled out geological and seismic observations suggest that differential attenuation is the main factor causing variation in the observed reflectivity thus thick granites act as low attenuating windows to the lower crust

earthquake rupture dynamics frozen in exhumed ancient faults just by looking at the traces of solidified friction induced melt produced by ancient earthquakes it is possible to learn about fractures that lasted a few seconds during seismic events that took place millions of years ago in the italian alps large sections of an ancient fault have been uprooted by tectonic movements brought to the surface almost intact and polished by the abrasion of a glacier secondary fractures created during the earthquakes tend to lie to one side of the fault in a particular direction suggesting that rupture propagated in the same direction during repeated earthquakes most of our knowledge about co seismic rupture propagation is derived from inversion and interpretation of strong ground motion seismograms1 2 3 laboratory experiments on rock4 5 and rock analogue material6 or inferred from theoretical and numerical elastodynamic models7 8 9 however additional information on dynamic rupture processes can be provided by direct observation of faults exhumed at the earth s surface10 pseudotachylytes solidified friction induced melts11 12 are the most certain fault rock indicator of seismicity on ancient faults13 here we show how the asymmetry in distribution and the orientation of pseudotachylyte filled secondary fractures around an exhumed fault can be used to reconstruct the earthquake rupture directivity rupture velocity and fracture energy by comparison with the theoretical dynamic stress field computed around propagating fractures in particular the studied natural network of pseudotachylytes is consistent with a dominant propagation direction during repeated seismic events and subsonic rupture propagation close to the rayleigh wave velocity

temporal variation of crustal deformation during the days preceding a thrust type great earthquake the 1944 tonankai earthquake of magnitude 8 1 japan the temporal variation in precursory ground tilt prior to the 1944 tonankai japan earthquake which is a great thrust type earthquake along the nankai trough is discussed using the analysis of data from repeated surveys along short distance leveling routes sato 1970 pointed out that an anomalous tilt occurred one day before the earthquake at kakegawa near the northern end of the focal region of the earthquake from the analysis of additional leveling data sato s result is re examined and the temporal change in the ground tilt is deduced for the period of about ten days beginning six days before the earthquake a remarkable precursory tilt started two or three days before the earthquake the direction of the precursory tilt was up towards the south uplift on the southern nankai trough side but the coseismic tilt was up towards the southeast perpendicular to the strike of the main thrust fault of the tonankai earthquake the postseismic tilt was probably opposite of the coseismic tilt the preseismic tilt is attributed to precursory slip on part of the main fault if similar precursory deformation occurs before a future earthquake expected to occur in the adjacent tokai region the deformation may help predict the time of the tokai earthquake

lower crustal reflectivity patterns in wide angle seismic recordings many features of wide angle reflection data have been attributed to finely layered structure in the lower crust reflections with similar character however can be produced by small amplitude 2 d random velocity variations finite difference seismograms for a lower crustal model with isotropic random velocity variations match field data recorded at the black forest germany at source receiver offsets beyond about twice the target depth waves scattered from a random zone produce the appearance of horizontal layering in wide angle data lower crustal reflectivity patterns can thus be explained by models that have lateral and vertical random fluctutations rather than fine scale layering