seismic imaging of subduction zone metamorphism combined analysis of high resolution seismic images of the alaska and cascadia subduction zones reveals where metamorphic fl uids are released both images show the subducted oceanic crust as a dipping low velocity layer with a clear termination depth however in alaska the crust is thicker 15 20 km compared to 8 km and terminates at greater depth 120 km com pared to 40 km than in cascadia based on metamorphic reaction estimates and geodynamic models we demonstrate that the termination depth corresponds to eclogitization of the crust triggered by dehydration of water bearing minerals and that the location of this reaction is dependent on the thermal structure of the subducted slab

space geodetic observations of nazca south america convergence across the central andes space geodetic data recorded rates and directions of motion across the convergent boundary zone between the oceanic nazca and continental south american plates in peru and bolivia roughly half of the overall convergence about 30 to 40 millimeters per year accumulated on the locked plate interface and can be released in future earthquakes about 10 to 15 millimeters per year of crustal shortening occurred inland at the sub andean foreland fold and thrust belt indicating that the andes are continuing to build little 5 to 10 millimeters per year along trench motion of coastal forearc slivers was observed despite the oblique convergence

spatial distribution of focal mechanisms for interplate and intraplate earthquakes associated with the subducting pacific plate beneath the northeastern japan arc a triple planed deep seismic zone the northeastern japan arc is located in one of the most seismically active subduction zones in the world in this study we relocated hypocenters and determined focal mechanisms of small earthquakes m 5 beneath the arc in order to investigate in detail the stress distribution in and around the descending oceanic plate in the hypocenter relocation we adopted a source region station correction method in which station corrections vary with hypocenter locations we have developed a new focal mechanism determination technique named master solution method which is analogous to the master event method in hypocenter determination we applied the method to p and sh wave amplitude data to obtain 1106 focal mechanism solutions from the new mechanism solutions and relocated hypocenters we found that there occur both low angle thrust fault lt type and downdip compression dc type earthquakes at depths from 40 to 70 km near the aseismic front the dc type events are underlying the lt type events almost all the earthquake clusters are composed of lt type events the western limit of the region where lt type events have occurred is subparallel to the trench axis although it undulates considerably it delineates the westernend of the active region of interplate seismicity furthermore we found that normal fault nf type events also occur at depths from 70 km in the upper plane of the double planed deep seismic zone which is characterized mainly by dc type event these nf type events are distributed only in a thin uppermost portion of the slab close to the plate boundary below these in the lower plane are downdip extension de type events this result indicates that the deep seismic zone in the northeastern japan arc is not double planed but triple planed even beneath the land area which cannot be explained by any simple models

mantle phase changes and deep earthquake faulting in subducting lithosphere inclined zones of earthquakes are the primary expression of lithosphere subduction a distinct deep population of subduction zone earthquakes occurs at depths of 350 to 690 kilometers at those depths ordinary brittle fracture and frictional sliding the faulting processes of shallow earthquakes are not expected a fresh understanding of these deep earthquakes comes from developments in several areas of experimental and theoretical geophysics including the discovery and characterization of transformational faulting a shear instability connected with localized phase transformations under nonhydrostatic stress these developments support the hypothesis that deep earthquakes represent transformational faulting in a wedge of olivine rich peridotite that is likely to persist metastably in coldest plate interiors to depths as great as 690 km predictions based on this deep structure of mantle phase changes are consistent with the global depth distribution of deep earthquakes the maximum depths of earthquakes in individual subductions zones and key source characteristics of deep events

evidence for lithospheric detachment in the central andes from local earthquake tomography abstract data from three temporary seismic networks were merged for tomographic inversion although the deployments did not coincide in time spatial overlap was achieved by re occupying existing sites travel times and t operators of about 1600 earthquakes were inverted for 3d models of p p s and p wave attenuation q p 1 all three attributes provide a consistent image of the entire subduction zone on a lithospheric scale the tomographic images reveal low velocities and high attenuation in the crust and mantle underlying the western cordillera and most of the puna plateau indicative of weak rheology and mostly asthenospheric mantle in contrast forearc and eastern foreland are characterized by high q p values corresponding to cold temperatures in accordance with thermal models in the backarc between 23 s and 24 s a high velocity high q p structure beneath the eastern cordillera and eastern puna is interpreted as detaching continental lithosphere that has been thickened in the orogenic process south of this structure the mantle is characterized by low velocities high p s ratios and low q p values here it is believed that lithosphere originally underlying andean crust has already been removed this is supported by new estimates of crustal thickness and volcanic activity

relations among subduction parameters clues to the dynamics of the subduction process are found in the many measurable parameters of modern subduction zones based on a critical appraisal of the geophysical and geological literature 26 parameters are estimated for each of 39 modern subduction zones to isolate causal relationships among these parameters multivariate analysis is applied to this data set this analysis yields empirical quantitative relations that predict strain regime and strike slip faulting in the overriding plate maximum earthquake magnitude benioff zone length slab dip arc trench gap and maximum trench depth excellent correlation is found between length of the benioff zone and the product of convergence rate and age of the downgoing slab this relationship is consistent with the conductive heating model of molnar et al 1979 if the model is modified in one respect the rate of heating of the slab is not constant it is substantially slower during passage of the slab beneath the accretionary prism and overriding plate the structural style in the overriding plate is determined by its stress state though the stress state of overriding plates cannot be quantified one can classify each individual subduction zone into one of seven semiquantitative strain classes that form a continuum from strongly extensional class 1 back arc spreading to strongly compressional class 7 active folding and thrusting this analysis indicates that strain class is probably determined by a linear combination of convergence rate slab age and shallow slab dip interplate coupling controlled by convergence rate and slab age is an important control on strain regime and the primary control on earthquake magnitude arc parallel strike slip faulting is a common feature of convergent margins forming a forearc sliver between the strike slip fault and trench optimum conditions for the development of forearc slivers are oblique convergence a compressional environment and a continental overriding plate the primary factor controlling presence of strike slip faulting is coupling strongly oblique convergence is not required the rate of strike slip faulting is affected by both convergence obliquity and convergence rate maximum trench depth is a response to flexure of the underthrusting plate the amount of flexural deflection at the trench depends on the vertical component of slab pull force which is very sensitive to slab age and shallow slab dip shallow slab dip conforms to the cross sectional shape of the overriding plate which is controlled by width of the accretionary prism and duration of subduction deep slab dip is affected by the mantle trajectory established at shallow depth but may be modified by mantle flow much of the structural complexity of convergent margins is probably attributable to terrane juxtaposition associated with temporal changes in both forearc strike slip faulting and strain regime empirical equations relating subduction parameters can provide both a focus for future theoretical studies and a conceptual and kinematic link between plate tectonics and the geology of subduction zones

seismic structure of the northeastern japan convergent margin a synthesis many studies recently made on the basis of seismic observations have revealed a detailed structure of the crust and the upper mantle beneath the northeastern japan arc and its relationship to seismic and volcanic activity spatial distributions of the depths to the conrad and the moho discontinuities estimated from shallow earthquake data and seismic explosion data show that both discontinuities are deep in the middle of the land area and shallow toward the coastlines of the japan sea and the pacific ocean the pn velocity has a lateral variation it is as low as 7 5 km s beneath the land area while that beneath the japan sea and the pacific ocean is 8 0 8 2 km s it changes abruptly at the transition zones which are located along the coastlines precise structure and location of the subducted pacific plate beneath the land area is inferred from converted or reflected seismic waves at the top or bottom of the plate the pacific plate is composed of a thin 5 km low velocity upper layer and a thick high velocity lower layer its total thickness being 80 90 km the upper plane seismicity of the double seismic zone is confined to the thin low velocity upper layer which probably corresponds to the subducted former oceanic crust the lower plane seismicity lies at the middle of the high velocity lower layer and the lower half of the plate below it is incapable of generating earthquakes the shallower portion of the upper surface of the plate beneath the pacific ocean along which major seismicity with low angle thrust faultings is actually occurring is also located by seismic observations on land and in the sea the pacific plate subducts at an extremely low angle of 5 for the first 25 km depth and then the dip steepens rather abruptly to 30 normal fault type events at the top of the plate have not been detected in the portion where the downward bending is the largest but have been ditected near the trench axis where it is rather small tomographic inversions for seismic velocity structure clearly delineate the inclined high velocity pacific plate with a thickness of 80 90 km and low velocity zones in the crust and the mantle wedge beneath active volcanoes seismic attenuation tomography also shows similar zones of low q value beneath active volcanoes although its spatial resolution is much lower the low velocity zones with 2 6 velocity lows are continuously distributed from the upper crust just beneath active volcanoes to a depth of 100 150 km in the mantle wedge and are approximately parallel to the dip of the underlying pacific plate these low velocity zones probably reflect the pathway of magma ascent from a depth in the mantle wedge to the earth s surface corresponding to a portion of the subduction induced secondary mantle wedge flow

depositional environment reservoir units evolution and hydrocarbon habitat of shuaiba formation lower cretaceous abu dhabi united arab emirates the shuaiba formation aptian is a thick porous shelf carbonate present in the subsurface over much of the arabian gulf it is the uppermost formation of the thamama group lower cretaceous and is overlain by the nahr umr formation albian and overlies the kharaib formation barremian the formation is divided into nine reservoir units a to i from the base upward reservoir units a b and c represent calcareous algal platform sediments showing some open marine influence they are dominated by chalky lime mudstones to wackestones with abundant codiacian green algae bacinella irregularis and lithocodium aggregatum and stylolite seams reservoir units d and f are deeper open marine shelf facies dominated by dense argillaceous lime mudstones to wackestones whic are burrowed and bioturbated reservoir units g and i represent lagoonal facies dominated by dense argillaceous wackestones to packstones which are sparsely fossiliferous and commonly stylolitic reservoir unit h represents a shallow water coarser grained shelf limestone containing abundant rudists algae orbitolinids and other foraminifers reservoir unit ho represents a shoal facies characterized by wackestone to packstone sediment with abundant orbitolina and scattered rudist debris and other foraminifers in northeast and southeast dhabi the basinal shuaiba equivalent bab member was deposited under slightly deeper quieter marine conditions and is dominated in parts of the section by dense argillaceous lime mudstones and shaly wackestone and calcareous shale i tervals the major reefal shuaiba buildup in abu dhabi trends northwest southeast it is in this trend that the major oil fields shah bu hasa and zurarah shaybah were discovered and from which some wells produce in western abu dhabi e g ruwais 2 hydrocarbons also occur outside of the rudist trend in the shuaiba formation in the asab jarn yaphour and sahil onshore oil fields and the zakum and mandous offshore oil fields local reefal anomalies are also encountered in parts of the deep marine basin of eastern abu dhabi these are formed by an unusual pre reefal growth of orbitolinids rudists and algae an isolated buildup patch reef consisting of coralgal boundstones encountered in jarn yaphour and zabbara fields may represent organic growth over either a rising salt structure or an uplifted zone formed by deep seated faults geochemical studies reveal that the main source for the thamama oils are the upper jurassic dukhan formation onshore abu dhabi or diyab formation offshore abu dhabi however some shuaiba source rocks may have existed in intrashelf basinal areas in the upper thamama zone the nahr umr shale acts as a seal for the shuaiba reservoir in this area

megathrust earthquake swarms indicate frictional changes which delimit large earthquake ruptures abstract using a recently constructed catalog of megathrust earthquake swarms we present results that illustrate how fault slip for m w 8 5 earthquakes with a publicly available finite fault rupture model is delimited along strike by swarm regions we further investigate this relationship for all m w 7 5 subduction zone megathrust earthquakes that have occurred over the last 40 years by estimating the rupture area from early aftershocks this extended comparison shows only 2 out of 42 earthquakes propagate through regions of the megathrust that have documented recurring earthquake swarm activity while about 40 17 out of 42 are terminate on at least one side near regions with swarm activity in japan chile sumatra and alaska earthquake swarms correlate with regions of the plate interface that exhibit low interseismic strain accumulation this suggests earthquake swarms are associated with reduced plate coupling and possibly aseismic moment release however earthquakes in chile and alaska ruptured through regions of weak coupling that do not show swarm activity we propose instead that stress heterogeneity on the megathrust can best explain why large earthquake ruptures terminate in swarm regions but not in regions with weak coupling alone earthquake swarms appear to be a proxy for along strike segmentation of subduction megathrusts at least on the time scale of one seismic cycle

landslide deformation behavior influenced by water level fluctuations of the three gorges reservoir china abstract water level fluctuations that occurs during filling drawdown operations in man made reservoirs can result in accelerations or reactivations or even catastrophic failure of large landslides for disaster mitigation a monitoring system comprising manual and automatic gps systems borehole inclinometers and groundwater piezometers was installed to investigate the slope deformation behavior of a large reactivated landslide known as the shuping landslide at the three gorges reservoir china the long term monitoring data 2003 2015 focused on the landslide deformation which was triggered by water level fluctuations is presented in this study the cumulative displacement curves showed a step like behavior suggesting a lower base activity superimposed by phases of acceleration which appeared during reservoir drawdown periods and rainy seasons based on statistical pearson correlation analysis the increase in landslide activity is dominated by reservoir water drawdown rather than rainfall the durations of the slide acceleration phases are much shorter than the reservoir water drawdown durations based on in site piezometer monitoring data and numerical simulation the potential causes are the hydraulic processes in the slope with much smaller groundwater drawdown velocities and hydraulic gradients in the initial stage of the reservoir water drawdown period

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

intense foreshocks and a slow slip event preceded the 2014 iquique mw 8 1 earthquake the subduction zone in northern chile is a well identified seismic gap that last ruptured in 1877 the moment magnitude mw 8 1 iquique earthquake of 1 april 2014 broke a highly coupled portion of this gap to understand the seismicity preceding this event we studied the location and mechanisms of the foreshocks and computed global positioning system gps time series at stations located on shore seismicity off the coast of iquique started to increase in january 2014 after 16 march several mw 6 events occurred near the low coupled zone these events migrated northward for 50 kilometers until the 1 april earthquake occurred on 16 march on shore continuous gps stations detected a westward motion that we model as a slow slip event situated in the same area where the mainshock occurred

tracking repeated subduction of the hikurangi plateau beneath new zealand abstract the hikurangi plateau has been subducted beneath new zealand twice firstly at ca 100 ma during n s convergence with gondwana and currently during e w convergence between the pacific and australian plates here we use a nationwide 3 d seismic velocity model derived from seismic tomography studies and seismicity relocated with this model to determine the location of the subducted western and southern edges of the plateau in the hikurangi subduction zone in the north island and northern south island we identify the western edge in the 37 140 km depth range based on the distribution of seismic velocities in the top of the subducted pacific plate changes in the thickness of the dipping seismic zone and steepening of the subducted plate downdip of the buoyant plateau in the central south island where the hikurangi plateau is subducting for a second time the plateau dips less steeply its southern edge is well defined by the collision of the delaminated lower crust of the plateau with the northern fiordland subduction zone where it bends the younger subducted australian plate to vertical the hikurangi plateau was originally rifted from the ontong java plateau and the subducted portions of the hikurangi plateau identified here fill a gap in reconstructions of this large igneous province

high resolution image of the geometry and thickness of the subducting nazca lithosphere beneath northern chile 1 results obtained from s and p receiver functions produced a clear image of the top and bottom of the subducting nazca lithosphere beneath northern chile using data from the teleseismic events recorded at 15 permanent integrated plate boundary observatory chile ipoc stations we obtained new constraints on the geometry and thickness of the descending nazca lithosphere we observed the subducted crust of the nazca plate at depths ranging from 50 km beneath the coastal cordillera down to 110 km beneath the western cordillera we found significant along strike variations in the geometry of the nazca plate beneath northern chile on closer inspection it appears that the oceanic nazca plate is divided into two distinct segments as it descends beneath the continental south american plate the transition from the relatively steeper 23 and deeper slab to the north of 21 s to the flatter southern segment 19 is shown reasonably clearly by our data this feature could well be associated with variations in the curvature of the plate margin and the geometry of the chile trench which is mainly curved to the north of 21 s we have also mapped the continental moho of the south american plate at depths ranging between 60 and 70 km to the east of the longitudinal valley beneath the coastal cordillera this boundary becomes invisible probably due to the serpentinization of the forearc mantle wedge that reduces the velocity in the uppermost mantle the base of the subducted nazca plate was clearly identified as a sharp boundary in the results obtained from the p and s receiver functions the thickness of the subducted oceanic nazca plate which has an age of 50 my is estimated to be 50 km although this thickness is consistent with that predicted by thermal gradients the explanation of the sharpness of the lithosphere asthenosphere boundary may require another mechanism such as hydration or melting

geophysical evidence for hydration of the crust and mantle of the nazca plate during bending at the north chile trench water transported in subducting oceanic plates plays a key role in a number of phenomena including intraslab seismicity and arc magmatism however the locus of plate hydration and water distribution in crust and mantle of plates entering subduction zones is debated we present evidence for anomalously low seismic velocities and densities of the crust and upper mantle of the nazca plate at the north chile trench crustal seismic velocities at the trench are lower than velocities of mature fast spreading crust and even lower than velocities of highly extended slow spreading crust in addition the nazca plate at the north chile trench may contain an 20 km thick upper mantle layer with 17 serpentine which implies 2 5 wt water these results document pervasive rock alteration by water percolation linked to bending related extensional faulting

the 2007 m7 7 tocopilla northern chile earthquake sequence implications for along strike and downdip rupture segmentation and megathrust frictional behavior 1 in 2007 a m7 7 earthquake occurred near the town of tocopilla within the northern chile seismic gap main shock slip derived from coseismic surface deformation was confined to the depth range between 30 and 55 km we relocated 1100 events during six months before and one week after the main shock aftershock seismicity is first congruent to the main shock slip and then it spreads offshore west and northwest of mejillones peninsula mp waveform modeling for 38 aftershocks reveals source mechanisms that are in the majority similar to the main shock however a few events appear to occur in the upper plate some with extensional mechanisms juxtaposing the tocopilla aftershocks with those following the neighboring 1995 antofagasta earthquake produces a striking symmetry across an ew axis in the center of mp events seem to skirt around mp probably due to a shallower moho there we suggest that the seismogenic coupling zone in northern chile changes its frictional behavior in the downdip direction from unstable to mostly conditionally stable for both earthquake sequences aftershocks agglomerate in the conditionally stable region whereas maximum inter seismic slip deficit and co seismic slip occurs in the unstable region the boundary between the unstable and conditionally stable zones parallels the coastline we identify a similar segmentation for other earthquakes in chile and peru where the offshore segments break in great m 8 earthquakes and the onshore segments in smaller m 8 earthquakes using critical taper analysis we demonstrate a causal relationship between varying slip behavior on the interface and forearc wedge anatomy that can be attributed to spatial variations in the rate dependency of friction

gradual unlocking of plate boundary controlled initiation of the 2014 iquique earthquake a long foreshock series unlocked the south american plate boundary until eventually initiating the m 8 1 iquique chile earthquake two groups publishing in this issue of nature analyse the seismic context of the iquique earthquake that occurred off the coast of northern chile on 1 april 2014 in a seismic zone that had been quiescent since a significant event in 1877 gavin hayes et al identify areas of remaining or elevated earthquake hazard along the megathrust fault in the region and conclude that the 2014 iquique event was not the earthquake that had been anticipated given that significant sections of the northern chile subduction zone have not ruptured in almost 150 years they suggest that it is likely that future megathrust earthquakes will occur south and potentially north of the 2014 iquique sequence bernd schurr et al show that the april 2014 earthquake broke a central fraction of the northern chile seismic gap the last major segment of the south american plate boundary that had yet to rupture in the past century from july 2013 up to the april earthquake they identify three seismic clusters along this part of the plate boundary each lasting a few weeks with earthquakes of increasing peak magnitudes they conclude that these seismic clusters and their slip transients reflect a gradual weakening of the central part of the seismic gap that was instrumental in initiating the final failure on 1 april 2014 northern chile was struck by a magnitude 8 1 earthquake following a protracted series of foreshocks the integrated plate boundary observatory chile monitored the entire sequence of events providing unprecedented resolution of the build up to the main event and its rupture evolution here we show that the iquique earthquake broke a central fraction of the so called northern chile seismic gap the last major segment of the south american plate boundary that had not ruptured in the past century1 2 since july 2013 three seismic clusters each lasting a few weeks hit this part of the plate boundary with earthquakes of increasing peak magnitudes starting with the second cluster geodetic observations show surface displacements that can be associated with slip on the plate interface these seismic clusters and their slip transients occupied a part of the plate interface that was transitional between a fully locked and a creeping portion leading up to this earthquake the b value of the foreshocks gradually decreased during the years before the earthquake reversing its trend a few days before the iquique earthquake the mainshock finally nucleated at the northern end of the foreshock area which skirted a locked patch and ruptured mainly downdip towards higher locking peak slip was attained immediately downdip of the foreshock region and at the margin of the locked patch we conclude that gradual weakening of the central part of the seismic gap accentuated by the foreshock activity in a zone of intermediate seismic coupling was instrumental in causing final failure distinguishing the iquique earthquake from most great earthquakes finally only one third of the gap was broken and the remaining locked segments now pose a significant increased seismic hazard with the potential to host an earthquake with a magnitude of 8 5

two dimensional velocity models of the nazca plate subduction zone between 19 5 s and 25 s from wide angle seismic measurements during the cinca95 project the interpretation of seismic refraction wide angle reflection data from the 1995 crustal investigations off and on shore nazca central andes cinca95 project has resulted in the derivation of nine e w two dimensional 2 d velocity cross sections for the region between the peru chile trench and the coast between 19 5 s and 25 s with three of the cross sections extending a farther 100 200 km inland these sections define the major lithospheric structures of the upper south american plate and the lower nazca plate down to uppermost mantle depths of 30 60 km beneath this part of the present day forearc region in addition to showing the nazca plate subducting at an increasing angle of 9 25 down to 30 50 km depth near the coast these cross sections show a portion of the moho dipping eastward from 43 50 km near the coast to 55 64 km up to about 240 km inland owing to a gap of almost 50 km in the data coverage of the moho near the coast it is uncertain from this data set alone whether the moho defined east of the coast should be correlated with the lower oceanic nazca plate or the upper continental south american plate however a comparison with other seismological data suggests that the moho identified here east of the coast defines the base of the continental crust of the upper south american plate immediately behind the downgoing nazca plate thus the hypothesis of delouis et al 1996 that the zone of seismic coupling between the plates correlates with the region where continental crust is in contact with the plate boundary is supported the cross sections also show between 20 s and 22 5 s a boundary extending from upper crustal levels down to 25 30 km depth subparallel to and 5 10 km above the top surface of the subducting nazca plate the discussion focuses on either that the boundary is a fossil extensional structure or preferably that the wedge below the boundary represents gabbroic material from layer 3 of the oceanic crust tectonically underplated onto the lower part of the upper plate

aftershock seismicity of the 27 february 2010 mw 8 8 maule earthquake rupture zone on 27 february 2010 the mw 8 8 maule earthquake in central chile ruptured a seismic gap where significant strain had accumulated since 1835 shortly after the mainshock a dense network of temporary seismic stations was installed along the whole rupture zone in order to capture the aftershock activity here we present the aftershock distribution and first motion polarity focal mechanisms based on automatic detection algorithms and picking engines by processing the seismic data between 15 march and 30 september 2010 from stations from iris ipgp gfz and university of liverpool we determined 20 205 hypocentres with magnitudes mw between 1 and 5 5 seismic activity occurs in six groups 1 normal faulting outer rise events 2 a shallow group of plate interface seismicity apparent at 25 35 km depth and 50 120 km distance to the trench with some variations between profiles along strike the aftershocks occur largely within the zone of coseismic slip but extend 50 km further north and with predominantly shallowly dipping thrust mechanisms along dip the events are either within the zone of coseismic slip or downdip from it depending on the coseismic slip model used 3 a third band of seismicity is observed further downdip at 40 50 km depth and further inland at 150 160 km trench perpendicular distance with mostly shallow dipping 28 thrust focal mechanisms indicating rupture of the plate interface significantly downdip of the coseismic rupture and presumably above the intersection of the continental moho with the plate interface 4 a deep group of intermediate depth events between 80 and 120 km depth is present north of 36 s within the maule segment a large portion of events during the inter seismic phase originated from this depth range 5 the magmatic arc exhibits a small amount of crustal seismicity but does not appear to show significantly enhanced activity after the mw 8 8 maule 2010 earthquake 6 pronounced crustal aftershock activity with mainly normal faulting mechanisms is found in the region of pichilemu 34 5 s these crustal events occur in a 30 km wide region with sharp inclined boundaries and oriented oblique to the trench the best located events describe a plane dipping to the southwest consistent with one of the focal planes of the large normal faulting aftershock mw 6 9 on 11 march 2010

space geodetic estimation of the nazca south america euler vector abstract gps data from four sites in the nazca plate easter island galapagos robinson crusoe and san felix islands and from five sites in the stable core of the south american plate enabled us to estimate the euler vector of the nazca plate with respect to south america the observed velocities of easter island 6 6 cm yr at 102 3 galapagos 5 1 cm yr at 90 0 robinson crusoe 6 6 cm yr at 80 1 and san felix 6 0 cm yr at 82 1 are significantly slower than the global plate model nuvel 1a velocites for those four sites our estimated euler pole is located at 48 8 n 91 7 w with a rate of 0 59 m y the nuvel 1a and earlier space geodetic studies give rotation rates that are 20 faster

existence of interplane earthquakes and neutral stress boundary between the upper and lower planes of the double seismic zone beneath tohoku and hokkaido northeastern japan abstract using data from the recently constructed nationwide dense seismic network we determined the hypocenters and focal mechanisms of many intermediate depth intraslab earthquakes within the pacific slab beneath northeastern japan the results show that in addition to the upper and lower planes of the double seismic zone a considerable number of intraslab earthquakes also occur between the two planes this interplane earthquake activity is not homogeneously distributed in space being high beneath eastern hokkaido and the fore arc regions of southeast and central tohoku the focal mechanisms of the interplane earthquakes tend to be of the down dip compressional type dc in tohoku and the hokkaido corner but down dip tensional type de in eastern hokkaido upper plane earthquakes are characterized by dc type stress while lower plane earthquakes are de type in both tohoku and eastern hokkaido the existence of interplane earthquakes enables estimation of the position of the neutral plane between the upper plane dc stress and the lower plane de stress we did so by applying stress tensor inversions to focal mechanism data obtained in the present study and japan meteorological agency jma data the results show that the neutral plane is located about 22 km beneath the upper plate interface beneath tohoku but only about 10 km beneath the upper surface beneath eastern hokkaido this difference in the location of the neutral plane may be due to the difference in buoyancy force exerted by the less dense metastable olivine wedge that is a result of the oblique plate subduction beneath hokkaido comparison of large intraslab earthquakes beneath the two regions shows that their aftershock areas are limited by the neutral plane suggesting that large earthquake ruptures are confined to either the dc or de stress field and do not go beyond the neutral plane

continuing megathrust earthquake potential in chile after the 2014 iquique earthquake the 2014 iquique event was not the earthquake that had been expected to fill the regional seismic gap given that significant sections of the northern chile subduction zone have not ruptured in almost 150 years it is likely that future megathrust earthquakes will occur south and potentially north of the 2014 iquique sequence two groups publishing in this issue of nature analyse the seismic context of the iquique earthquake that occurred off the coast of northern chile on 1 april 2014 in a seismic zone that had been quiescent since a significant event in 1877 gavin hayes et al identify areas of remaining or elevated earthquake hazard along the megathrust fault in the region and conclude that the 2014 iquique event was not the earthquake that had been anticipated given that significant sections of the northern chile subduction zone have not ruptured in almost 150 years they suggest that it is likely that future megathrust earthquakes will occur south and potentially north of the 2014 iquique sequence bernd schurr et al show that the april 2014 earthquake broke a central fraction of the northern chile seismic gap the last major segment of the south american plate boundary that had yet to rupture in the past century from july 2013 up to the april earthquake they identify three seismic clusters along this part of the plate boundary each lasting a few weeks with earthquakes of increasing peak magnitudes they conclude that these seismic clusters and their slip transients reflect a gradual weakening of the central part of the seismic gap that was instrumental in initiating the final failure the seismic gap theory1 identifies regions of elevated hazard based on a lack of recent seismicity in comparison with other portions of a fault it has successfully explained past earthquakes see for example ref 2 and is useful for qualitatively describing where large earthquakes might occur a large earthquake had been expected in the subduction zone adjacent to northern chile3 4 5 6 which had not ruptured in a megathrust earthquake since a m 8 8 event in 1877 on 1 april 2014 a m 8 2 earthquake occurred within this seismic gap here we present an assessment of the seismotectonics of the march april 2014 iquique sequence including analyses of earthquake relocations moment tensors finite fault models moment deficit calculations and cumulative coulomb stress transfer this ensemble of information allows us to place the sequence within the context of regional seismicity and to identify areas of remaining and or elevated hazard our results constrain the size and spatial extent of rupture and indicate that this was not the earthquake that had been anticipated significant sections of the northern chile subduction zone have not ruptured in almost 150 years so it is likely that future megathrust earthquakes will occur to the south and potentially to the north of the 2014 iquique sequence

accurate hypocentre determination in the seismogenic zone of the subducting nazca plate in northern chile using a combined on offshore network the coupled plate interface of subduction zones commonly called the seismogenic zone has been recognized as the origin of fatal earthquakes a subset of the after shock series of the great antofagasta thrust type event 1995 july 30 mw 8 0 has been used to study the extent of the seismogenic zone in northern chile to achieve reliable and precise hypocentre locations we applied the concept of the minimum 1 d model which incorporates iterative simultaneous inversion of velocity and hypocentre parameters the minimum 1 d model is complemented by station corrections which are influenced by near surface velocity heterogeneity and by the individual station elevations by relocating mine blasts which were not included in the inversion we obtain absolute location errors of 1 km in epicentre and 2 km in focal depth a study of the resolution parameters ale and dspr documents the importance of offshore stations on location accuracy for offshore events based on precisely determined hypo centres we calculate a depth of 46 km for the lower limit of the seismogenic zone which is in good agreement with previous studies for this area for the upper limit we found a depth of 20 km our results of an aseismic zone between the upper limit of the seismogenic zone and the surface correlates with a detachment zone proposed by other studies the results are also in agreement with thermal studies for the antofagasta forearc region

the seismic coupling of subduction zones revisited 1 the nature of seismic coupling for many of the world s subduction zones has been reevaluated geodetic estimates of seismic coupling obtained from gps measurements of upper plate deformation during the interseismic period are summarized we compared those with new estimates of seismic coupling obtained from seismological data the results show that with a few notable exceptions the two methods agree to within about 10 the seismological estimates have been greatly improved over those made 20 30 years ago because of an abundance of paleoseismological data that greatly extend the temporal record of great subduction earthquakes and by the occurrence in the intervening years of an unusual number of great and giant earthquakes that have filled in some of the most critical holes in the seismic record the data also again with a few notable exceptions support the frictional instability theory of seismic coupling and in particular the test of that theory made by scholz and campos 1995 overall the results support their prediction that high coupling occurs for subduction zones subjected to high normal forces with a switch to low coupling occurring fairly abruptly as the normal force decreases below a critical value there is also considerable variation of coupling within individual subduction zones earthquake asperities correlate with areas of high coupling and hence have a semblance of permanence but the rupture zones and asperity distributions of great earthquakes may differ greatly between seismic cycles because of differences in the phase of seismic flux accumulation

double seismic zone and dehydration embrittlement of the subducting slab 1 dehydration embrittlement of metamorphosed oceanic crust and mantle in the subducting slab may be responsible for the occurrence of intermediate depth earthquakes we explore the possibility that this hypothesis can explain the morphology of the double seismic zones observed in northeast japan southwest japan northeast taiwan northern chile cape mendocino and eastern aleutians we calculate transient temperature structures of slabs based on geologically estimated subduction histories of these regions we then determine dehydration loci of metamorphosed oceanic crust and serpentinized mantle using experimentally derived phase diagrams the depth range of the dehydration loci of metamorphosed oceanic crust and serpentine is dependent on slab age the dehydration loci of serpentine produce a double layered structure because the upper dehydration loci of serpentine are mostly located in the wedge mantle above the slab we regard the upper plane seismicity representing dehydration embrittlement in the oceanic crust and we fix the slab geometry so that the upper plane seismicity is just below the upper surface of the slab we find that the lower plane seismicity is located at the lower dehydration loci of serpentine which indicates that the morphology of the double seismic zones is consistent with the dehydration embrittlement

double seismic zone of the nazca plate in northern chile high resolution velocity structure petrological implications and thermomechanical modeling this paper presents an interdisciplinary study of the northern chile double seismic zone first a high resolution velocity structure of the subducting nazca plate has been obtained by the tomodd double difference tomography method the double seismic zone dsz is observed between 80 and 140 km depth and the two seismic planes is 20 km apart then the chemical and petrologic characteristics of the oceanic lithosphere associated with this dsz are deduced by using current thermal petrological seismological models and are compared to pressure temperature conditions provided by a numerical thermomechanical model our results agree with the common hypothesis that seismicity in both upper and lower planes is related to fluid releases associated with metamorphic dehydration reactions in the seismic upper plane located within the upper crust these reactions would affect material of basaltic morb composition and document different metamorphic reactions occurring within high p 2 4 gpa and low t 130 km lawsonite amphibole eclogite conditions the lower plane lying in the oceanic mantle can be associated with serpentinite dehydration reactions the vp and vs characteristics of the region in between both planes are consistent with a partially 25 30 vol antigorite 0 10 vol brucite and 4 10 vol chlorite hydrated harzburgitic material discrepancies persist that we attribute to complexities inherent to heterogeneous structural compositions while various geophysical indicators evidence particularly cold conditions in both the descending nazca plate and the continental fore arc thermomechanical models indicate that both seismic planes delimit the inner slab compressional zone around the 400 c 50 c isotherm lower plane earthquakes are predicted to occur in the slab s flexural neutral plane where fluids released from surrounding metamorphic reactions could accumulate and trigger seismicity fluids migrating upward from the tensile zone below could be blocked in their ascension by the compressive zone above this plane thus producing a sheeted layer of free fluids or a serpentinized layer therefore earthquakes may present either downdip compression and downdip tensile characteristics numerical tests indicate that the slab s thermal structure is not the only factor that controls the occurrence of inner slab compression 1 a weak ductile subduction channel and 2 a cold mantle fore arc both favor inner slab compression by facilitating transmission of compressional stresses from the continental lithosphere into the slab 3 decreasing the radius of curvature of the slab broadens the depth of inner slab compression whereas 4 decreasing upper plate convergence diminishes its intensity all these factors indicate that if dszs indeed contour inner slab compression they cannot be linked only to slab unbending but also to the transmission of high compressional stresses from the upper plate into the slab

revisiting the north chile seismic gap segmentation using gps derived interseismic coupling no major earthquake occurred in north chile since the 1877 m w 8 6 subduction earthquake that produced a huge tsunami however geodetic measurements conducted over the last decade in this area show that the upper plate is actually deforming which reveals some degree of locking on the subduction interface this accumulation of elastic deformation is likely to be released in a future earthquake because of the long elapsed time since 1877 and the rapid accumulation of deformation thought to be 6 7 cm yr many consider this area is a mature seismic gap where a major earthquake is due and seismic hazard is high we present a new global positioning system gps velocity field acquired between 2008 and 2012 that describes in some detail the interseismic deformation between 18 s and 24 s we invert for coupling distribution on the nazca south america subduction interface using elastic modelling our measurements require that at these latitudes 10 to 12 mm yr 1 i e 15 per cent of the whole convergence rate are accommodated by the clockwise rotation of an andean block bounded to the east by the subandean fold and thrust belt this reduces the accumulation rate on the subduction interface to 56 mm yr 1 in this area coupling variations on the subduction interface both along strike and along dip are described we find that the north chile seismic gap is segmented in at least two highly locked segments bounded by narrow areas of weak coupling this coupling segmentation is consistent with our knowledge of the historical ruptures and of the instrumental seismicity of the region intersegment zones iquique mejillones correlate with high background seismic rate and local tectonic complexities on the upper or downgoing plates the rupture of either the paranal or the loa segment alone could easily produce a mw 8 0 8 3 rupture and we propose that the loa segment from 22 5 s to 20 8 s may be the one that ruptured in 1877

three dimensional models of p wave velocity and p to s velocity ratio in the southern central andes by simultaneous inversion of local earthquake data the pisco 94 proyecto de investigation sismologica de la cordillera occidental 1994 seismological network of 31 digital broad band and short period three component seismometers was deployed in northern chile between the coastal cordillera and the western cordillera more than 5300 local seismic events were observed in a 100 day period a subset of high quality p and s arrival time data was used to invert simultaneously for hypocenters and velocity structure additional data from two other networks in the region could be included the velocity models show a number of prominent anomalies outlining an extremely thickened crust about 70 km beneath the forearc region an anomalous crustal structure beneath the recent magmatic arc western cordillera characterized by very low velocities and a high velocity slab a region of an increased vp vs ratio has been found directly above the wadati benioff zone which might be caused by hydration processes a zone of lower than average velocities and a high vp vs ratio might correspond to the asthenospheric wedge the upper edge of the wadati benioff zone is sharply defined by intermediate depth hypocenters while evidence for a double seismic zone can hardly be seen crustal events between the precordillera and the western cordillera have been observed for the first time and are mainly located in the vicinity of the salar de atacama down to depths of about 40 km

reflection image spectroscopy across the andean subduction zone abstract this paper presents new insights into the south american subduction zone from reprocessed seismic images we applied a 3d kirchhoff prestack depth migration scheme to data sets containing different narrow frequency ranges in order to extract additional details from seismic reflection images this approach accounts for the effects of scattering on the seismic image especially for structures below a heterogeneous overburden the reflection image in such environments will differ significantly when focusing on different frequency ranges due to the frequency dependence of scattering that is likely to be present the narrow frequency range images uncover reflectors in one frequency range that are masked in another range furthermore the images enable for instance the characterization of the medium in terms of scatterer concentration and thus improve the structural interpretation the analysis of these images might help to distinguish between small scale structures in the high frequency band and large scale structures in the low frequency band we call this imaging approach reflection image spectroscopy ris we applied the ris approach to the ancorp 96 data set an onshore deep seismic reflection profile across the south american central andes the narrow frequency range images revealed additional details that we interpret as features directly linked to fluid migration processes in the subduction zone furthermore structural details of the oceanic crust and the overlying mantle and crust are revealed from the narrow frequency range images we interpret that the top of the so called nazca reflector at 70 km depth marks the upper limit of the hydrated mantle wedge whereas the bottom of the reflector represents the top of the subducted oceanic crust the compilation of our results with local earthquake data confirms this interpretation similar features observed in another deep seismic profile precorp 95 support this interpretation too furthermore the ris images show a highly reflective heterogeneous zone between the nazca reflector at 70 km depth and a prominent mid crustal bright spot quebrada blanca bright spot at about 30 km depth we associate this zone with a complex network of ascending fluids or partial melts initiated by ascending fluids released from the subducted oceanic plate this observation links the quebrada bright spot directly to the dehydrating oceanic plate

seismic imaging of a convergent continental margin and plateau in the central andes andean continental research project 1996 ancorp 96 1 a 400 km long seismic reflection profile andean continental research project 1996 ancorp 96 and integrated geophysical experiments wide angle seismology passive seismology gravity and magnetotelluric depth sounding across the central andes 21 s observed subduction of the nazca plate under the south american continent an east dipping reflector nazca reflector is linked to the down going oceanic crust and shows increasing downdip intensity before gradual breakdown below 80 km we interpret parts of the nazca reflector as a fluid trap located at the front of recent hydration and shearing of the mantle the fluids being supplied by dehydration of the oceanic plate patches of bright quebrada blanca bright spot to more diffuse reflectivity underlie the plateau domain at 15 30 km depth this reflectivity is associated with a low velocity zone p to s wave conversions the upper limits of high conductivity and high vp vs ratios and to the occurrence of neogene volcanic rocks at surface we interpret this feature as evidence of widespread partial melting of the plateau crust causing decoupling of the upper and lower crust during neogene shortening and plateau growth the imaging properties of the continental moho beneath the andes indicate a broad transitional character of the crust mantle boundary owing to active processes like hydration of mantle rocks in the cooler parts of the plate margin system magmatic underplating and intraplating under and into the lowermost crust mechanical instability at moho etc hence all first order features appear to be related to fluid assisted processes in a subduction setting

unbending combined with dehydration embrittlement as a cause for double and triple seismic zones 1 using a model of elastic perfectly plastic plate i demonstrate that plate unbending combined with dehydration embrittlement can account for the morphology and stresses of double or triple seismic zones dsz s or tsz s observed in many subducting slabs densifying metamorphic reactions may induce tensile stresses as well as high pore fluid pressures in the slab crust and thus generate the upper seismicity band of the tsz limited dehydration may moderately elevate the fluid pressure in the slab mantle but unbending usually dominates mantle stresses and results in another two seismicity bands dsz s are underdeveloped tsz s

tectonic and magmatic responses to the subduction of high bathymetric relief subduction of high bathymetric relief such as aseismic ridges and magmatic plateaus is considered to be responsible for dramatic changes in the dynamics and kinematics of the subduction zone for example the buoyancy of high bathymetric relief is thought to flatten the dip of the subducting slab modifying the structural and magmatic evolution of the overriding plate and terminating arc volcanism in addition the effect of ridge subduction in retreating plate boundaries can inhibit subduction rollback a process that could locally pin the subduction hinge and lead to the development of cusps and slab tearing here we discuss the tectonic response to subduction of high bathymetric relief using examples from the circum pacific subduction systems we demonstrate that flattening of the subduction dip angle is only significant in the eastern pacific where the average slab dip angle is relatively shallow in the western pacific in contrast the average subduction dip angle is steeper and there is no significant flattening of the dip angle in areas of ridge subduction subduction of high bathymetric relief in the circum pacific is commonly associated with reduced arc volcanism and in many cases the area of ridge subduction coincides with a volcanic gap in the overriding plate ridge subduction is associated with pronounced changes in the style of deformation involving uplift reactivation of basement thrusts development of orogen perpendicular tear faults and block rotations leading to oroclinal bending the discussed characteristic patterns associated with ridge subduction provide important guidelines for reconstructing past plate tectonic processes and could help constraining the geodynamics of ancient subduction systems

exposed plate interface in the european alps reveals fabric styles and gradients related to an ancient seismogenic coupling zone 1 we present observations from a continuous exposure of an ancient plate interface in the depth range of its former seismogenic zone in the central alps of europe related to late cretaceous early tertiary subduction and accretion of the south penninic lower plate underneath the adriatic upper plate the material forming the exposed plate interface zone has experienced flow and fracturing over an extended period of time followed by syncollisional exhumation thus reflecting a multistage evolution fabric formation and metamorphism however chiefly record the deformation conditions of the precollisional setting along the plate interface we identify an unstable slip domain from pseudotachylytes occurring in the temperature range between 200 and 300 c this zone coincides with a domain of intense veining in the subduction melange with mineral growth into open cavities indicating fast possibly seismic rupture evidence for transient near lithostatic fluid pressure as well as brittle fractures competing with mylonitic shear zones continues into the region below the occurrence of pseudotachylytes possibly reflecting a zone of conditionally stable slip the zone above the unstable slip area is devoid of veins but displays ample evidence of fluid assisted processes similar to the deeper zone solution precipitation creep and dehydration reactions in the melange matrix hydration and sealing of the base of the upper plate seismic rupture here is possibly expressed by ubiquitous localized deformation zones we hypothesize that trenchward sealing of parts of the plate interface as well as reaction enhanced destruction of upper plate permeability is an important component localizing the unstable slip zone this relation may result from the competition of the pervasive presumably interseismic pressure solution creep destroying permeability and building elevated fluid pressure until the strength threshold is reached with seismic failure

splay fault branching along the nankai subduction zone seismic reflection profiles reveal steeply landward dipping splay faults in the rupture area of the magnitude m 8 1 tonankai earthquake in the nankai subduction zone these splay faults branch upward from the plate boundary interface that is the subduction zone at a depth of approximately 10 kilometers approximately 50 to 55 kilometers landward of the trough axis breaking through the upper crustal plate slip on the active splay fault may be an important mechanism that accommodates the elastic strain caused by relative plate motion

intermediate depth earthquake faulting by dehydration embrittlement with negative volume change earthquakes are observed to occur in subduction zones to depths of approximately 680 km even though unassisted brittle failure is inhibited at depths greater than about 50 km owing to the high pressures and temperatures it is thought that such earthquakes particularly those at intermediate depths of 50 300 km may instead be triggered by embrittlement accompanying dehydration of hydrous minerals principally serpentine a problem with failure by serpentine dehydration is that the volume change accompanying dehydration becomes negative at pressures of 2 4 gpa 60 120 km depth above which brittle fracture mechanics predicts that the instability should be quenched here we show that dehydration of antigorite serpentinite under stress results in faults delineated by ultrafine grained solid reaction products formed during dehydration this phenomenon was observed under all conditions tested pressures of 1 6 gpa temperatures of 650 820 degrees c independent of the sign of the volume change of reaction although this result contradicts expectations from fracture mechanics it can be explained by separation of fluid from solid residue before and during faulting a hypothesis supported by our observations these observations confirm that dehydration embrittlement is a viable mechanism for nucleating earthquakes independent of depth as long as there are hydrous minerals breaking down under a differential stress

locking of the chile subduction zone controlled by fluid pressure before the 2010 earthquake large subduction zone earthquakes are thought to occur where the down going and overriding tectonic plates are strongly locked analysis of geodetic and seismic data collected in the decade before the 2010 chile earthquake shows that variations in pore fluid pressure correlate with the degree of plate interface locking and may therefore control earthquake rupture

long term persistence of subduction earthquake segment boundaries evidence from mejillones peninsula northern chile 1 the mejillones peninsula in northern chile has been recognized as the surface expression of a segment boundary for large subduction zone earthquakes the sharp contact between the rupture planes of two instrumentally recorded earthquakes the mw 8 0 antofagasta 1995 and the mw 7 7 tocopilla 2007 events is located beneath the central part of mejillones peninsula we present new chronostratigraphic and structural data that allow reconstructing the evolution of the peninsula at the surface and correlation of the latter with seismic cycle deformation on the plate interface uplift commenced after 3 4 myr as recorded in the western highland the central graben area on the peninsula started uplifting above sea level as an anticlinal hinge zone prior to 400 kyr ago most probably 790 kyr ago the resulting e w trending hinge exactly overlies the limit between the rupture planes of the antofagasta and tocopilla earthquakes by correlating the uplift data with the slip distribution of the above earthquakes we demonstrate that deformation and uplift is focused during the postseismic and interseismic periods of the megathrust seismic cycle with coseismic deformation opposed to the long term motion additionally the slip deficit beneath the peninsula accumulating between events is probably largely recovered by creep hence we suggest that mejillones peninsula owes its existence to the lateral variation of the propensity for unstable slip at the interface since the latter is a material property the long term spatial stability of the peninsula as a barrier to rupture propagation since at least the middle pleistocene is a necessary consequence

slab1 0 a three dimensional model of global subduction zone geometries 1 we describe and present a new model of global subduction zone geometries called slab1 0 an extension of previous efforts to constrain the two dimensional non planar geometry of subduction zones around the focus of large earthquakes slab1 0 describes the detailed non planar three dimensional geometry of approximately 85 of subduction zones worldwide while the model focuses on the detailed form of each slab from their trenches through the seismogenic zone where it combines data sets from active source and passive seismology it also continues to the limits of their seismic extent in the upper mid mantle providing a uniform approach to the definition of the entire seismically active slab geometry examples are shown for two well constrained global locations models for many other regions are available and can be freely downloaded in several formats from our new slab1 0 website http on doi gov d9arbs we describe improvements in our two dimensional geometry constraint inversion including the use of average active source seismic data profiles in the shallow trench regions where data are otherwise lacking derived from the interpolation between other active source seismic data along strike in the same subduction zone we include several analyses of the uncertainty and robustness of our three dimensional interpolation methods in addition we use the filtered subduction related earthquake data sets compiled to build slab1 0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far concluding that the width of these seismogenic zones is on average 30 larger than previous studies have suggested

the kinematic and geodynamic significance of the atacama fault zone northern chile abstract the atacama fault zone afz is the dominant feature in the structure of the north chilean coastal cordillera left lateral displacement took place along its system of longitudinal faults during the jurassic and early cretaceous this development was contemporaneous with arc magmatism and was later reactivated resulting in a steep normal fault strike slip movements along the afz consist of two sets of ductile shear zones of different ages one jurassic formed under amphibolite facies conditions the other early cretaceous with greenschist facies mylonites structural asymmetries point to a sinistral sense of shear in both sets the afz can be interpreted as a magmatic arc structure which accommodated the oblique subduction of an oceanic plate trench linked strike slip fault the sinistral sense of shear is consistent with reconstructions of late jurassic to early cretaceous plate configurations in the se pacific

are the lower planes of double seismic zones caused by serpentine dehydration in subducting oceanic mantle in a number of subduction zones earthquakes at 50 200 km depth define two dipping planes separated by 20 40 km that appear to merge downdip upper plane earthquakes are inferred to occur within the subducting oceanic crust whereas lower plane earthquakes occur in the subducting oceanic mantle beneath northeast japan the results of a new finite element heat transfer model suggest that the lower seismic plane cuts across isotherms at a shallow angle lower plane earthquakes occur at 550 800 c at 100 km depth and at 350 600 c at 160 km depth these conditions coincide with the dehydration reaction antigorite serpentine forsterite enstatite h 2 o which suggests that lower plane earthquakes may be triggered by dehydration embrittlement which in turn suggests that the subducting oceanic mantle is partially hydrated serpentinization may occur in the trench outer rise region where faulting may promote infiltration of seawater several tens of kilometers into the oceanic lithosphere if this hypothesis is correct current subduction zone h 2 o budgets may significantly underestimate the amount of bound h 2 o entering the subduction factory

seismicity and average velocities beneath the argentine puna plateau a network of 60 seismographs was deployed across the andes at 23 5 s the array was centered in the backarc atop the puna high plateau in nw argentina p and s arrival times of 426 intermediate depth earthquakes were inverted for 1 d velocity structure and hypocentral coordinates average velocities and v p v s in the crust are low average mantle velocities are high but difficult to interpret because of the presence of a fast velocity slab at depth although the hypocenters sharply define a 35 dipping benioff zone seismicity in the slab is not continuous the spatial clustering of earthquakes is thought to reflect inherited heterogeneties of the subducted oceanic lithosphere additionally 57 crustal earthquakes were located seismicity concentrates in the fold and thrust belt of the foreland and eastern cordillera and along and south of the el toro olacapato calama lineament tocl focal mechanisms of two earthquakes at this structure exhibit left lateral strike slip mechanisms similar to the suggested kinematics of the tocl we believe that the puna north of the tocl behaves like a rigid block with little internal deformation whereas the area south of the tocl is weaker and currently deforming

the mw 8 1 antofagasta north chile earthquake of july 30 1995 first results from teleseismic and geodetic data a strong mw 8 1 subduction earthquake occurred on july 30 1995 in antofagasta northern chile this is one of the largest events during this century in the region it ruptured the southernmost portion of a seismic gap between 18 s and 25 s in 1992 we had used gps to survey a network with about 50 benchmarks covering a region nearly 500 km long n s and 200 km wide e w part of these marks were re surveyed with gps after the 1995 earthquake comparison with 1992 positions indicate relative horizontal displacement towards the trench reaching 0 7 m the inland subsided several decimeters the mejillones peninsula was uplifted by more than 15 cm teleseismic body wave modelling of vbb records gives a subduction focal mechanism and source time function with three distinct episodes of moment release and southward directivity modelling the displacement field using a dislocation with uniform slip in elastic half space suggests a rupture zone with 19 24 eastward dip extending to a depth no greater than 50 km with n s length of 180 km and an average slip of about 5 m the component of right lateral slip inferred both from the teleseismic and geodetic data does not require slip partitioning at the plate boundary that the well constrained northern end of the 1995 rupture zone lies under the southern part of the mejillones peninsula increases the probability for a next rupture in the gap north of it

revisiting viscoelastic effects on interseismic deformation and locking degree a case study of the peru north chile subduction zone viscoelastic effects potentially play an important role during all phases of the earthquake cycle in subduction zones however most current models neglect such effects in the interseismic deformation pattern here we use finite element method fem models to investigate the control of viscoelasticity on interseismic deformation and to highlight the pitfalls of interpreting the data with purely elastic models for both the forward and inverse problems our results confirm that elastic models are prone to overestimating the interseismic locking depth a crucial parameter for estimating the maximum possible earthquake magnitude the application of the viscoelastic model improves the fit to the interseismic deformation especially in the inland area additionally we construct 3 d fem models constrained by geophysical and gps data and apply our methodology to the peru north chile subduction zone our results indicate that viscoelastic effects contribute significantly to the observed gps data the signals interpreted as back arc shortening in the elastic model can be alternatively explained by viscoelastic deformation which in turn dramatically refines the interseismic locking pattern in both dip and strike directions our viscoelastic locking map exhibits excellent correlation with the slip distributions of previous earthquakes especially the recent 2014 mw 8 1 iquique earthquake the incorrect elastic assumptions affect the analysis of interseismic deformation with respect to slip deficit calculations our results thus suggest that it is necessary to thoroughly reevaluate existing locking models that are based on purely elastic models some of which attribute viscoelastic deformation to different sources such as microplate sliver motions

three dimensional density model of the nazca plate and the andean continental margin 1 we forward modeled the bouguer anomaly in a region encompassing the pacific ocean 85 w and the andean margin 60 w between northern peru 5 s and patagonia 45 s the three dimensional density model that reproduces the gravity field is a continental scale representation of density structure to 410 km depth that characterizes the mantle and crust of the oceanic nazca plate subducted slab and continental margin with a minimum number of bodies we predefined the density of each body after studying the dependency of density on composition of crustal and mantle materials and pressure temperature conditions appropriate for the andean setting a database of independent geophysical information constrains the geometry of the top of the subducted slab locally the moho of the oceanic and continental crusts and indirectly the lithosphere asthenosphere boundary underneath the continental plate other boundaries notably the intracrustal density discontinuity separating upper from lower crust below the continent were not constrained and their geometry is the result of fitting the observed and calculated bouguer anomaly during forward modeling this contribution presents the model to the andean geoscientific community and contains some tools like a sensitivity analysis that helps potential users of the model to interpret its results we describe and discuss some of these results in order to illustrate the application of the model to the study of a wide range of phenomena e g modification of oceanic plate structure by hot spots shape of the subducted slab thermal structure of the continental lithosphere compensation mechanism and formation of orogenic relieve causes of andean segmentation

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

variation with spreading rate of oceanic crustal thickness and geochemistry abstract away from the influence of fracture zones hotspots and marginal basins oceanic crust exhibits remarkably uniform crustal thickness rare earth element concentrations and bulk composition at all but the most slowly spreading ridges at full spreading rates below 15 mm a however all these observed parameters show marked and abrupt changes the seismically determined crustal thicknesses decrease sharply from their average of 7 1 km as do melt thicknesses inferred from rare earth element inversions of basalts bulk element compositions of basalts show an increase in the percentage of na2o and decreases in the feo content and cao al2o3 ratio as the full spreading rate drops below 15 mm a we present a model of the melting beneath oceanic spreading centres which reproduces the observed variations with spreading rate of crustal thickness and the main geochemical changes observed in basalts changes in the volume and composition of melt generated at very slow spreading rates are caused by conductive heat loss from the mantle welling up beneath rifts we model the observed behaviour with normal mantle potential temperatures of 1300 20 c greater crustal thicknesses are generated at the north atlantic spreading centre adjacent to the iceland mantle plume suggesting mantle temperatures about 80 c hotter than normal in this area

existence of a seismic belt in the upper plane of the double seismic zone extending in the along arc direction at depths of 70 100 km beneath ne japan 1 we relocated microearthquakes using data obtained via a dense seismic network and systematically detected the characteristic distribution of the upper seismic plane seismicity within the pacific slab beneath ne japan we found a seismic belt parallel to the iso depth contour of the plate interface that is beneath the forearc area at depths of 70 100 km indicating that the distribution of the upper plane earthquakes is non uniform the location of the deeper limits of this belt and seismicity of the upper seismic plane appear to correspond respectively to two facies boundaries where h2o contents change in the slab crust events in the upper seismic plane have mainly down dip compression type focal mechanisms but several events have normal fault type nf type versions whose spatial distribution appears to correspond to these boundaries these nf events might be induced by the tensional stress field that is caused by volume reduction due to dehydration reactions

the geometry of the wadati benioff zone under southern central america and its tectonic significance results from a high resolution local seismographic network abstract we present here a detailed geometry of the wadati benioff zone under costa rica obtained from seismicity recorded by a dense local seismographic network jointly operated by the costa rica volcanological and seismological observatory national university and the charles f richter seismological laboratory university of california santa cruz underneath the nicaragua costa rica border the wadati benioff zone smoothly contorts from steep to shallow dip angles nw to se but does not show evidence of a brittle tear as postulated by others however further to the se ne of puerto caldera costa rica the wadati benioff zone does show a segmentation the quesada sharp contortion at intermediate depths h 70km nw of this sharp contortion the deepest portion of the seismically active slab dips at about 80 and reaches maximum depths ranging from 200 km near the nicaragua costa rica border to 135 km under ciudad quesada to the se the deeper portion of the wadati benioff zone dips at about 60 and the seismicity does not extend below depths ranging from 125 km behind the volcanic arc to 50 km east of quepos in southern costa rica east of 83 55 w we find no evidence of the wadati benioff zone deeper than 50 km the obtained geometry and other known tectonic features related to the subduction of the cocos plate under the caribbean plate along the southern terminus of the middle america trench nicaragua and costa rica correlate well with along trench variations in age of the subducted cocos plate some of these tectonic features are 1 the shallowing of middle america trench bathymetry from nw to se 2 variations in the energy release within the subducted slab 3 differences in coupling between cocos and caribbean plates 4 the termination of the central america volcanic chain in central costa rica 5 distinct stress field variations on the overriding caribbean plate the subduction of the cocos ridge under southern costa rica is partially responsible for some of these features shallowing of the middle america trench in southern costa rica and for the high uplift rates of the outer arc however as the presence of the panama fracture zone limits the subducted extension of the cocos ridge to less than 100 km from the trench we propose that the overall geometry of the wadati benioff zone is controlled by abrupt along trench changes in the age of the subducted cocos plate

geometry of the pamir hindu kush intermediate depth earthquake zone from local seismic data 1 we present new seismicity images based on a two year seismic deployment in the pamir and sw tien shan a total of 9532 earthquakes were detected located and rigorously assessed in a multistage automatic procedure utilizing state of the art picking algorithms waveform cross correlation and multi event relocation the obtained catalog provides new information on crustal seismicity and reveals the geometry and internal structure of the pamir hindu kush intermediate depth seismic zone with improved detail and resolution the relocated seismicity clearly defines at least two distinct planes one beneath the pamir and the other beneath the hindu kush separated by a gap across which strike and dip directions change abruptly the pamir seismic zone forms a thin approximately 10 km width curviplanar arc that strikes east west and dips south at its eastern end and then progressively turns by 90 to reach a north south strike and a due eastward dip at its southwestern termination pamir deep seismicity outlines several streaks at depths between 70 and 240 km with the deepest events occurring at its southwestern end intermediate depth earthquakes are clearly separated from shallow crustal seismicity which is confined to the uppermost 20 25 km the hindu kush seismic zone extends from 40 to 250 km depth and generally strikes east west yet bends northeast toward the pamir at its eastern end it may be divided vertically into upper and lower parts separated by a gap at approximately 150 km depth in the upper part events form a plane that is 15 25 km thick in cross section and dips sub vertically north to northwest seismic activity is more virile in the lower part where several distinct clusters form a complex pattern of sub parallel planes the observed geometry could be reconciled either with a model of two sided subduction of eurasian and previously underthrusted indian continental lithosphere or by a purely eurasian origin of both pamir and hindu kush seismic zones which necessitates a contortion and oversteepening of the latter

deep slab subduction and dehydration and their geodynamic consequences evidence from seismology and mineral physics abstract we present new pieces of evidence from seismology and mineral physics for the existence of low velocity zones in the deep part of the upper mantle wedge and the mantle transition zone that are caused by fluids from the deep subduction and deep dehydration of the pacific and philippine sea slabs under western pacific and east asia the pacific slab is subducting beneath the japan islands and japan sea with intermediate depth and deep earthquakes down to 600 km depth under the east asia margin and the slab becomes stagnant in the mantle transition zone under east china the western edge of the stagnant pacific slab is roughly coincident with the ne sw daxing anling taihangshan gravity lineament located west of beijing approximately 2000 km away from the japan trench the upper mantle above the stagnant slab under east asia forms a big mantle wedge bmw corner flow in the bmw and deep slab dehydration may have caused asthenospheric upwelling lithospheric thinning continental rift systems and intraplate volcanism in northeast asia the philippine sea slab has subducted down to the mantle transition zone depth under western japan and ryukyu back arc though the seismicity within the slab occurs only down to 200 300 km depths combining with the corner flow in the mantle wedge deep dehydration of the subducting pacific slab has affected the morphology of the subducting philippine sea slab and its seismicity under southwest japan slow anomalies are also found in the mantle under the subducting pacific slab which may represent small mantle plumes or hot upwelling associated with the deep slab subduction slab dehydration may also take place after a continental plate subducts into the mantle

seismic consequences of warm versus cool subduction metamorphism examples from southwest and northeast japan warm and cool subduction zones exhibit differences in seismicity seismic structure and arc magmatism which reflect differences in metamorphic reactions occurring in subducting oceanic crust in southwest japan arc volcanism is sparse and intraslab earthquakes extend to 65 kilometers depth in northeast japan arc volcanism is more common and intraslab earthquakes reach 200 kilometers depth thermal petrologic models predict that oceanic crust subducting beneath southwest japan is 300 to 500 c warmer than beneath northeast japan resulting in shallower eclogite transformation and slab dehydration reactions and possible slab melting

young displacements on the atacama fault system northern chile from field observations and cosmogenic 21ne concentrations 1 we present the first numerical age constraint for young deformation of the atacama fault system afs in northern chile the young activity of the afs is expressed by several fault scarps which affects alluvial fan sediments of the eastern side of the coastal cordillera 23 30 23 42 s detailed mapping of alluvial fans reveals a complex relationship between fault motion erosion and alluvial fan development an older group of alluvial fans became inactive prior to the scarp formation younger alluvial fans arising directly from feeder channels and entrenched in the fault scarps posts date the scarp formation the youngest slip on the afs is recorded by headward eroding channels entrenched across the scarp which are in turn displaced vertically 0 3 0 5 m by the fault quartz fragments in four sites on the older inactive fan group were analyzed for cosmogenic 21ne concentrations yielding an average age of 424 151 ka the upper limit for the recent activity of the fault combined with the height of fault scarp we calculate a 0 01 mm yr minimum vertical fault slip rate thus young displacement on the afs is quaternary in age and confined to the late pleistocene