segmentation of the nazca and south american plates along the ecuador subduction zone from wide angle seismic profiles we describe the deep structure of the south colombian northern ecuador convergent margin using travel time inversion of wide angle seismic data recently collected offshore the margin appears segmented into three contrasting zones in the north zone affected by four great subduction earthquakes during the 20th century normal oceanic crust subducts beneath the oceanic cretaceous substratum of the margin underlined by seismic velocities as high as 6 0 6 5 km s in the central zone the subducting oceanic crust is over thickened beneath the carnegie ridge a steeper slope and a well developed high velocity cretaceous oceanic basement characterizes the margin wedge this area coincides with a gap in significant subduction earthquake activity in the south zone the subducting oceanic crust is normal the fore arc is characterized by large sedimentary basins suggesting significant subsidence velocities in the margin wedge are significantly lower and denote a different nature or a higher degree of fracturing even if the distance between the three profiles exceeds 150 km the structural segmentation obtained along the ecuadorian margin correlates well with the distribution of seismic activity and the neotectonic zonation

provenance of late carboniferous bauxite deposits in the north china craton new constraints on marginal arc construction and accretion processes abstract the north china craton ncc is bounded by two paleozoic accretionary arc terranes the north qinling terrane to the south and the bainaimiao terrane to the north the timing of arc accretion to the ncc and the architecture of the bainaimiao arc remain unclear during the building and accretion of the arcs along its margins the ncc experienced a long sedimentary hiatus since the ordovician which ended with the deposition of bauxite bearing sediments in the late carboniferous in this paper we report the u pb and hf isotopes of detrital zircons from the late carboniferous bauxite layer and use these data to constrain the tectonic evolution of the margin of the ncc the detrital zircons yield a minimum u pb age of ca 310 ma and a prominent age peak at ca 450 ma zircon crystals with ages of ca 330 ma and ca 1900 ma are more common in the bauxite samples from the northern part of the ncc than in those from the central part the ehf t values of the ca 450 ma detrital zircon crystals of the bauxite samples from the ncc are similar to those of the contemporaneous detrital zircon crystals from the north qinling arc terrane to the south but different from those of the contemporaneous detrital zircon crystals from the bainaimiao arc terrane to the north the ca 450 ma detrital zircon crystals in the ca 310 ma bauxite deposits are therefore interpreted to have been derived from the north qinling arc terrane the source of the ca 330 ma detrital zircon crystals of the bauxite deposits is interpreted to be the northern margin of the ncc where intermediate felsic plutons formed at ca 330 ma are common the results from this study support the interpretation that the paleozoic continental arc terranes and their concomitant back arc basins were developed along the margins of the ncc before ca 450 ma and the arc complexes were subsequently accreted to the craton in the late carboniferous this was then followed by the formation of a walled continental basin within the ncc

large intermediate depth earthquakes and the subduction process this study provides an overview of intermediate depth earthquake phenomena placing emphasis on the larger tectonically significant events and exploring the relation of intermediate depth earthquakes to shallower seismicity especially we examine whether intermediate depth events reflect the state of interplate coupling at subduction zones and whether this activity exhibits temporal changes associated with the occurrence of large underthrusting earthquakes historic record of large intraplate earthquakes m b 7 0 in this century shows that the new hebrides and tonga subduction zones have the largest number of large intraplate events regions associated with bends in the subducted lithosphere also have many large events e g altiplano and new ireland we compiled a catalog of focal mechanisms for events that occurred between 1960 and 1984 with m 6 and depth between 40 and 200 km the final catalog includes 335 events with 47 new focal mechanisms and is probably complete for earthquakes with m b 6 5 for events with m 6 5 nearly 48 of the events had no aftershocks and only 15 of the events had more than five aftershocks within one week of the mainshock events with more than ten aftershocks are located in regions associated with bends in the subducted slab focal mechanism solutions for intermediate depth earthquakes with m 6 8 can be grouped into four categories 1 normal fault events 44 and 2 reverse fault events 33 both with a strike nearly parallel to the trench axis 3 normal or reverse fault events with a strike significantly oblique to the trench axis 10 and 4 tear faulting events 13 the focal mechanisms of type 1 events occur mainly along strongly or moderately coupled subduction zones where a down dip extensional stress prevails in a gently dipping plate in contrast along decoupled subduction zones great normal fault earthquakes occur at shallow depths e g the 1977 sumbawa earthquake in indonesia type 2 events with strike subparallel to the subduction zone and most of them with a near vertical tension axis occur mainly in regions that have partially coupled or uncoupled subduction zones and the observed continuous seismicity is deeper than 300 km the increased dip of the downgoing slab associated with weakly coupled subduction zones and the weight of the slab may be responsible for the near vertical tensional stress at intermediate depth and consequently the change in focal mechanism from type 1 to type 2 events events of type 3 occur where the trench axis bends sharply causing horizontal parallel to the trench strike extensional or compressional intraplate stress type 4 are hinge faulting events for strongly coupled zones we observed temporal changes of intermediate depth earthquake activity associated with the occurrence of a large underthrusting event after the occurrence of a large underthrusting event the stress axis orientation of intermediate depth earthquakes changes from down dip tensional to down dip compressional e g 1960 chile 1974 peru 1982 tonga and 1952 kamchatka earthquakes or the number of large intermediate events decreases for a few years e g 1964 alaska and 1985 valparaiso earthquakes we conclude that even though the stress changes induced by slab pull and slab distortion control the general pattern of intermediate depth seismicity spatial and temporal variations of the intraplate stress associated with interplate coupling are important in controlling the global occurrence of large intermediate depth events

teleseismic time functions for large shallow subduction zone earthquakes broadband vertical p wave records are analyzed from 63 of the largest shallow subduction zone earthquakes which have occurred in the circum pacific in the last 45 yr most of the records studied come from a common instrument the pasadena california benioff 1 90 seismometer propagation and instrument effects are deconvolved from the p wave records using a damped least squares inversion to obtain the teleseismic source time function the inversion has the additional constraint that the time function be positive everywhere the period band over which the time functions are considered reliable is from 2 5 to 50 sec fourier displacement amplitude spectra computed for each of the 1 90 p wave trains indicate spectral slopes measured between 2 and 50 sec of 1 0 to 2 25 with an average value of 1 5 these values assume an average attenuation of t 1 0 the seismic moments derived from the p wave time functions compare well with other published values for earthquakes having moments smaller than 2 5 10 28 dyne cm m w 8 2 because the 1 90 seismometer has little response at very long periods this technique underestimates the moments of the very largest events the time functions are characterized using five parameters 1 spectral slope between 2 and 50 sec 2 roughness of the time function 3 multiplicity of sources 4 pulse widths of individual sources and 5 overall signal duration the 63 earthquakes studied come from 15 subduction zones with a wide range in the ages of subducted lithosphere convergence rates and maximum size of earthquakes comparing the time function parameters with age rate and m w of the subduction zone does not yield obvious global trends however most of the subduction zones do behave characteristically and can be grouped accordingly

source characteristics of the 12 november 1996 m w 7 7 peru subduction zone earthquake the 12 november 1996 m w 7 7 peru subduction zone earthquake occurred off the coast of southern peru near the intersection of the south american trench and the highest topographical point of the subducting nazca ridge we model the broadband teleseismic p waveforms from stations in the global seismic network to constrain the source characteristics of this subduction zone earthquake we have analyzed the vertical component p waves for this earthquake to constrain the depth source complexity seismic moment and rupture characteristics the seismic moment determined from the nondiffracted p waves is 3 5 1020n m corresponding to a moment magnitude m w of 7 6 7 7 the source time function for the 1996 peru event has three pulses of seismic moment release with a total duration of approximately 45 50 seconds the largest moment release occurs at approximately 35 40 seconds and is located 90km southeast of the rupture initiation approximately 70 of the seismic moment was released in the third pulse

supercycle at the ecuadorian subduction zone revealed after the 2016 pedernales earthquake large earthquakes are usually assumed to release all of the strain accumulated since the previous event implying a reduced seismic hazard after them however long records of seismic history at several subduction zones suggest supercycle behaviour where centuries long accumulated strain is released through clustered large earthquakes resulting in an extended period of enhanced seismic hazard here we combine historical seismology results present day geodesy data and dense local observations of the recent mw 7 8 2016 pedernales earthquake to reconstruct the strain budget at the ecuador subduction zone since the great 1906 earthquake we show that the pedernales earthquake involved the successive rupture of two patches on the plate interface that were locked prior to the earthquake and most probably overlaps the area already ruptured in 1942 by a similar earthquake however we find that coseismic slip in 2016 exceeds the deficit accumulated since 1942 the seismic moment of every large earthquake during the twentieth century further exceeds the moment deficit accumulated since 1906 these results together with the seismic quiescence before 1906 highlighted by historical records and marine palaeoseismology argue for an earthquake supercycle at the ecuador colombia margin this behaviour which has led to an enhanced seismic hazard for 110 years is possibly still going on and may apply to other subduction zones that recently experienced a great earthquake large earthquakes are often assumed to reset the seismic hazard of a region analysis of recent and historical seismicity in ecuador suggests that this region may experience clusters of large earthquakes and extended periods of high seismic hazard

influence of asperities along subduction interfaces on the stressing and seismicity of adjacent areas abstract dmowska r and lovison l c 1992 influence of asperities along subduction interfaces on the stressing and seismicity of adjacent areas in t mikumo k aki m ohnaka l j ruff and p k p spudich editors earthquake source physics and earthquake precursors tectonophysics 211 23 43 we have investigated the influence of large scale fault inhomogeneities in large subduction earthquakes on the style of deformation and seismic behavior of the incoming oceanic plate and slab at intermediate depths during the earthquake cycle the zones of the large subduction events of rat islands 1965 alaska 1964 and valparaiso 1985 have been searched for earthquakes with m b 5 0 if available and for time periods as long as possible it has been found that in general the seismicity in the incoming oceanic plate clusters in front of asperities areas of highest seismic moment release and strongest locking and is positioned relative to them in the direction of plate motion it is usually lacking in areas adjacent to non asperities that is to zones that slip during the main event but with appreciably smaller seismic moment release and possibly slip seismically aseismically during the whole cycle similar behavior occurs in the downgoing slab at intermediate depths where seismicity during the cycle clusters but less strongly than in the oceanic crust next to asperities and down dip from them we infer that the locking of asperities causes higher stresses associated with the earthquake cycle itself to occur in areas adjacent to asperities both up dip and down dip from them along the direction of plate motion and that such stressing is much less pronounced in the areas adjacent to non asperities this opens the possibility of identifying the areas of highest seismic moment release in future subduction earthquakes and carries implications for where the highest deformation and possibly precursory phenomena and or nucleation of a future event might occur

the rupture process and asperity distribution of three great earthquakes from long period diffracted p waves moment tensor source mechanisms for the largest twenty nine shallow earthquakes from 1981 m s 6 5 have been determined using long period rayleigh wave spectra and p wave first motions recorded by the international deployment of accelerograms ida global digital seismic network gdsn and world wide standardised seismic network wwssn these are compared with the centroid moment tensor mechanisms presented by dziewonski and woodhouse which were obtained in a contrasting manner in most cases the source is well represented by a single double couple mechanisms with good agreement of scalar moment and orientation between the two studies apparent mismatches for other cases may be interpreted in terms of source complexity we emphasise the objectivity of comparing moment tensor elements apparent discrepancies in the moment and orientation of double couple fault plane solutions are usually due to differences in the mxz and myz moment tensor elements which are poorly constrained by long period data

the 16 april 2016 mw 7 8 ms 7 5 ecuador earthquake a quasi repeat of the 1942 ms 7 5 earthquake and partial re rupture of the 1906 ms 8 6 colombia ecuador earthquake the 2016 ecuador m w 7 8 earthquake ruptured the subduction zone boundary between the nazca and south american plates joint modeling of seismic and tsunami observations indicates an 120 km long rupture area beneath the coastline north of the 1998 m w 7 2 rupture the slip distribution reveals two discrete asperities near the hypocenter and around the equator their locations and the patchy pattern are consistent with the prior interseismic geodetic strain which showed highly locked patches also beneath the coastline aftershocks cluster along two streaks one aligned nearly parallel to the plate convergence direction up dip of the main slip patches and the other on a trench perpendicular lineation south of the 1958 rupture zone comparisons of seismic waveforms and magnitudes show that the 2016 event and 1942 earthquakes have similar surface wave magnitude m s 7 5 overlapping rupture areas and similar main pulses of moment rate the same area ruptured as the southernmost portion of the larger earthquake of 1906 m w 8 6 m s 8 6 the seismic behavior reflects persistent heterogeneous frictional properties of the colombia ecuador megathrust

rupture of a seismic gap in southeastern alaska the 1972 sitka earthquake ms 7 6 abstract the sitka alaska earthquake of july 30 1972 m s 7 6 is one of the largest strike slip events since 1963 this earthquake is thought to have ruptured the pacific north america plate boundary segment between the large strike slip events in 1949 and 1958 and the sitka earthquake was heralded as one of the first successful forecasts of an earthquake based on the seismic gap hypothesis another motivation to study this event is the availability of a strong ground motion recording 50 km from the epicenter which is near the middle of the 180 km long plate boundary segment the teleseismically located aftershocks occur only in the southern half of this segment we have analyzed rayleigh waves and p waves recorded by wwssn long period vertical instruments a total of 25 very long period rayleigh phases r2 r5 define a symmetric four lobed radiation pattern from which we can infer a nearly pure right lateral strike slip fault striking 340 n20 w and a bilateral rupture comparison of synthetic and observed rayleigh waves yields a seismic moment of 3 10 27 dyn cm m w 7 6 the p waves are characterized by a dominant complex pulse in the first 40 s and a succession of smaller pulses after that the source time functions are deconvolved from 32 p waves with good azimuthal distribution to obtain the moment depth extent and rupture history the mean moment of the source time functions is 5 0 2 4 10 27 dyn cm and the depth extent is 0 10 km tomographic inversion of the source time functions produces an image of the earthquake rupture process this procedure maps the moment release of the complex pulse into two high displacement regions i e asperities 1 the epicentral area to 40 km north and 2 between 60 and 90 km south of the epicenter the main asperity just north of the epicenter coincides with a quiescent zone of microaftershocks while the overall average seismic displacement is 6 m the displacement at the main asperity is 7 m the same value as the tectonic displacement accumulated since the previous earthquake in 1848 the fault averaged stress drop is 100 bars while the asperity stress drop is 130 bars tomographic imaging of the source time functions thus resolves several features of the rupture process notably the bilateral propagation and rupture initiation at the largest asperity this is in contrast to the large 1976 guatemala strike slip earthquake where the largest displacements occurred away from rupture initiation

the tsunami signature on a submerged promontory the case study of the atacames promontory ecuador shelf promontories exhibit very specific bathymetric features with regards to tsunamis because of their submerged cape morphology a potential tsunami generated seawards of the promontory will exhibit a specific mode of propagation and coastal impact to identify this peculiar tsunami signature the atacames promontory ecuador was chosen as a case study another example is the shelf of the nile delta egypt the area is tectonically very active hosts earthquakes among the most powerful recorded as well as areas of slope instabilities that have triggered significant submarine landslides in the past several cubic kilometres of volume both types of events are likely to be tsunamigenic to examine the tsunami behaviour at the coastal area of the promontory and at its vicinity we have considered two examples of tsunamigenic landslides of which scars were identified near the base of the continental slope we also took into consideration two earthquake scenarios that are likely to represent most classes of earthquakes possibly occurring in this area depending on their locations and subsequent tsunami directivity that is a sensitivity test investigation we took two distinct earthquake scenarios which are based on the 1942 and 1958 events that stroke the area then we computed their derived tsunamis and analysed their coastal impact we found that significant tsunamis can be generated by either landslides or earthquakes however the maxima of wave amplitude occur offshore but still above the underwater promontory the concave type shape of the bathymetric field often yields a refraction focusing area that is located on the shelf promontory and not at the coast area of the promontory the wave propagates first through the focusing area before striking the considered coast this area may be considered as a sheltered zone besides in the vicinity of the promontory not exactly concerned by the study the city of esmeraldas is relatively sheltered due to the presence of the underwater canyon at its termination and due to diverging waves

study of the source process of the 1992 colombia ms 7 3 earthquake with the empirical green s function method the source time function stf of the october 18 1992 colombia ms 7 3 earthquake has been determined using long period teleseismic surface and body waveform data and the empirical green s function egf method love and sh waveforms of a ms 6 7 foreshock were treated as the egfs and deconvolved from the mainshock waveforms to obtain the relative stf this indicated that the ms 7 3 mainshock was a double event followed 100 seconds later by a large aftershock not previously reported the source durations for the mainshock and the aftershock are about 45 and 25 seconds respectively based on the source radiation directivity the rupture directions for the mainshock and aftershock are about s50 w the total length of the fault for the mainshock and aftershock is about 140 km which agrees with the fault dimension estimated using the aftershock distribution the seismic moment releases for the three asperities on the rupture surface are 1 0 1019 6 9 1019 and 2 3 1019 nm with respect to time

variable rupture mode of the subduction zone along the ecuador colombia coast three large earthquakes occurred within the rupture zone of the 1906 colombia ecuador earthquake m w 8 8 in 1942 m s 7 9 1958 ms 7 8 and 1979 m s 7 7 we compared the size and mechanism of these earthquakes by using long period surface waves tsunami data and macroseismic data the 1979 event is a thrust event with a seismic moment of 2 9 10 28 dyne cm and represents subduction of the nazca plate beneath south america the rupture length and direction are 230 km and n40 e respectively examination of old seismograms indicates that the 1906 event is also a thrust event which ruptured in the northeast direction the seismic moment estimated from the tsunami data and the size of the rupture zone is 2 10 29 dyne cm the 1942 and 1958 events are much smaller about 1 5 to 1 10 of the 1979 event in the seismic moment than the 1979 event we conclude that the sum of the seismic moments of the 1942 1958 and 1979 events is only formula of that of the 1906 event despite the fact that the sequence of the 1942 1958 and 1979 events ruptured approximately the same segment as the 1906 event this difference could be explained by an asperity model in which the fault zone is held by a discrete distribution of asperities with weak zones in between the weak zone normally behaves aseismically but slips abruptly only when it is driven by failure of the asperities a small earthquake represents failure of one asperity and the rupture zone is pinned at both ends by adjacent asperities so that the effective width and the amount of slip are relatively small a great earthquake represents failure of more than one asperity and consequently involves much larger width and slip

historical 1942 ecuador and 1942 peru subduction earthquakes and earthquake cycles along colombia ecuador and peru subduction segments two large shallow earthquakes occurred in 1942 along the south american subduction zone inclose proximity to subducting oceanic ridges the 14 may event occurred near the subducting carnegie ridge off the coast of ecuador and the 24 august event occurred off the coast of southwestern peru near the southern flank of the subducting nazca ridge source parameters for these for these two historic events have been determined using long periodp waveforms p wave first motions intensities and local tsunami data we have analyzed thep waves for these two earthquakes to constrain the focal mechanism depth source complexity and seismic moment modeling of thep waveform for both events yields a range of acceptable focal mechanisms and depths all of which are consistent with underthrusting of the nazca plate beneath the south american plate the source time function for the 1942 ecuador event has one simple pulse of moment release with a duration of 22 suconds suggesting that most of the moment release occurred near the epicenter the seismic moment determined from thep waves is 6 8 1020n m corresponding ot a moment magnitude of 7 8 7 9 the reported location of the maximum intensities ix for this event is south of the main shock epicenter the relocated aftershcks are in an area that is approximately 200 km by 90 km elongated parallel to the trench with the majority of aftershocks north of the epicenter in contrast the 1942 peru event has a much longer duration and higher degree of complexity than the ecuador earthquake suggesting a heterogeneous rupture seismic moment is released in three distinct pulses over approximately 74 seconds the largest moment release occurs 32 seconds after rupture initiation the seismic moment as determined from thep waves for the 1942 peru event is 10 25 1020n m corresponding to a moment magnitude of 7 9 8 2 aftershock locations reported by the iss occur over a broad area surrounding the main shock the reported locations of the maximum intensities ix are concentrated south of the epicenter suggesting that at least part of the rupture was to the south we have also examined great historic earthquakes along the colombia ecuador and peru segments of the south american subduction zone we find that the size and rupture length of the underthrusting earthquakes vary between successive earthquake cycles this suggests that the segmentation of the plate boundary as defined by earthquakes this century is not constant

inelastic models of lithospheric stress 11 implications for outer rise seismicity and dynamics summary outer rise seismicity and dynamics are examined using inelastic models of lithospheric deformation which allow a more realistic characterization of stress distributions and failure behaviour we conclude that thrust and normal faulting outer rise earthquakes represent substantially different states of stress within the oceanic lithosphere specifically the normal faulting events occur in response to downward plate bending which establishes the standard bending dominated state of outer rise stress and the thrust faulting events occur in response to an elevated level of in plane compression which develops only in response to exceptional circumstances this interpretation accounts for the observation that normal faulting outer rise earthquakes occur more frequently and are more widely distributed than their thrust faulting counterparts an observation for which the simple bending model offers no explanation in addition attributing both thrust and normal faulting outer rise earthquakes to plate bending implies that both classes of events should occur within relatively close lateral proximity to one another because both are allegedly a manifestation of the same bendingdominated stress distribution whereas in reality this is not observed we propose that the tendency for thrust faulting outer rise earthquakes to exhibit greater source depths than their normal faulting counterparts an observation that is frequently cited in support of the bending interpretation of the former is merely a consequence of the fact that bending induced tension is confined to the upper lithosphere our model predicts that outer rise in plane force variations may promote thrust faulting outerrise activity prior to an underthrusting interplate subduction earthquake and normalfaulting outer rise activity following such an earthquake but that both forms of outerrise activity are unlikely to be associated with the same subduction earthquake a corollary implication of our model is that subduction earthquakes are likely to be either preceded by or followed by an absence of large outer rise earthquakes levels of in plane compression necessary to generate thrust faulting outer rise earthquakes are attributed to stress concentrations within the subducting plate that are induced by relatively localized resistance to regionally distributed plate driving forces resistance of this nature may result from either the attempted subduction of relatively buoyant i e isostatically compensated bathymetric features or the existence of strong interplate

great earthquakes seismicity gaps and potential for earthquake disaster along the himalaya plate boundary abstract analysis of the space time patterns of seismicity in the himalaya plate boundary has established the existence of three seismic gaps 1 1 the kashmir gap lying west of the 1905 kangra earthquake 2 2 the central gap situated between the 1905 kangra and the 1934 bihar earthquakes 3 3 the assam gap between the 1897 and 1950 assam earthquakes this study has shown that the above great earthquakes were preceded as well as followed by long periods 19 years of decreased levels of seismic activity in the epicentral regions remarkable decrease in the seismicity following the year 1970 has been observed in the western half of the central gap as well as in the assam gap local seismic investigation in the assam gap confirms this feature and the seismicity suggests the existence there of an asperity the local seismic investigations in garhwal himalaya have shown that the small earthquakes are confined to the upper 6 8 km of the crust and may have strike slip motions these earthquakes occur in a region where teleseismically recorded events were few

fault size and depth extent of the ecuador earthquake m w 7 8 of 16 april 2016 from teleseismic and tsunami data the april 2016 ecuador mw 7 8 earthquake was the first megathrust tsunamigenic earthquake along the ecuador colombia subduction zone since 1979 mw 8 2 with 200 deaths from tsunami while there was no tsunami damage from the 2016 earthquake small tsunamis were recorded at deep ocean assessment and reporting of tsunami and tide gauges here we designed various fault models with and without shallow slip area and compared the computed teleseismic and tsunami waveforms with the observations while teleseismic inversions were indifferent about inclusion or exclusion of the shallow slip tsunami waveforms strongly favored the slip model without shallow slip our final slip model has a depth range of 15 44 km and its western shallowest limit is located at the distance of 60 km from the trench maximum and average slips were 2 5 and 0 7 m respectively the large slip area was 80 km along strike 60 km along dip in the depth range of 15 35 km

directivity analysis of the deep bolivian earthquake of june 9 1994 the june 9 1994 deep bolivian earthquake m w 8 3 d 636 km occurred on a nearly horizontal rupture plane with an approximate area of 2000 km 2 the rupture plane is constrained by directivity analysis of the p wave displacement pulses recorded at local regional and teleseismic distances the fact that there are local recordings of the displacement provides a strong constraint on the vertical extent of faulting which is less than 10 km the moment release had at least four distinct episodes or subevents the best overall explanation for the subevents is rupture which initiates on the western margin of the fault plane propagates to the east and then continues bilaterally along a 350 azimuth with most of the moment release occurring to the north the apparent rupture velocity between subevents is very low approximately 1 2 km s

a new gps velocity field for the pacific plate part 1 constraints on plate motion intraplate deformation and the viscosity of pacific basin asthenosphere

source mechanisms of colombia aftershocks using digital surface wave data love and rayleigh wave signals recorded at stations of the global digital seismograph network were used to retrieve focal mechanism information for moderate magnitude mb 5 0 aftershocks that followed the large m s 7 7 colombia earthquake of 12 december 1979 source mechanisms compatible with the observed love rayleigh spectral amplitude ratios were identified following a comparison against theoretical values calculated for a suite of source models p wave first motions observed at worldwide stations were then used to remove ambiguities in the sense of slip and strike orientation for most of the events examined the focal depths of the events were fixed at values determined from depth phase arrival times in the surface wave comparison a reference earthquake with known focal mechanism and depth was used to calibrate the procedure and to minimize the path and size effects aftershocks whose epicenters lie within the asperity region from which most of the moment release originated during the large colombia earthquake denote thrust faulting similar to that of the main shock consistent with occurrence along the interplate zone the few near trench events in the vicinity of the asperity may have resulted from secondary deformation away from the major plane of rupture rather than from a propagation of interplate rupture seaward to the trench axis these near trench aftershocks which are seaward of the region of maximum energy release do not provide strong evidence for seismogenic failure of the plate interface at least one of the events is a normal faulting earthquake that occurred either above the interplate boundary or in the subducting nazca plate source mechanisms computed for aftershocks located northeast of the asperity but within the length of the rupture reported for the main shock are also consistent with interplate slip the aftershock activity observed north of this distance is associated with a separate and perhaps offset segment of the nazca south america plate boundary

great earthquakes and subduction along the peru trench subduction along the peru trench between 9 and 15 s involves both large interplate underthrusting earthquakes and intraplate normal fault earthquakes the four largest earthquakes along the peru trench are from north to south the 1970 m 7 9 intraplate normal fault earthquake and the interplate underthrusting earthquakes in 1966 m 8 0 1940 m 8 and 1974 m 8 0 we have studied the rupture process of these earthquakes and can locate spatial concentrations of moment release through directivity analysis of source time functions deconvolved from long period p wave seismograms the 1966 earthquake has a source duration of 45 s with most of the moment release concentrated near the epicenter two intraplate normal fault events occurred in 1963 m 6 7 and 7 0 at the down dip edge of the 1966 dominant asperity the 1940 earthquake is an underthrusting event with a simple source time function of 30 s duration that represents the rupture of a single asperity near the epicenter the 1974 earthquake has a source duration of 45 50s and two pulses of moment release this earthquake has a bilateral rupture with the first pulse of moment release located northwest of the epicenter and the second pulse of moment release located southeast of the epicenter both pulses of moment release occur on the northern half of the aftershock area the 1970 earthquake is one of the largest intraplate normal fault earthquakes to occur in a subduction zone and has a moment release comparable with many large underthrusting events the aftershocks for the 1970 earthquake form two distinct clusters the smaller cluster near the epicenter has focal mechanisms characterized by down dip tension but the second aftershock cluster located 80 km southeast of the epicenter has focal mechanisms characterized by down dip compression the p waves for the main shock can be modeled as two sources with different focal mechanisms and depths similar to the two clusters of aftershocks the first event has a down dip tensional focal mechanism and is followed 40 s later by a distinct second event located 80 km southeast of the epicenter with a down dip compressional focal mechanism and a somewhat shallower depth than the first event the observable directivity indicates that the second source is located at the second cluster of aftershocks that have down dip compressional focal mechanisms the occurrence of both down dip tensional and compressional focal mechanisms may be explained by extreme unbending stresses associated with the anomalous slab geometry the unusually large size of the 1970 earthquake may also be related to the subduction of the mendana fracture zone the historic earthquake record along the peru trench indicates that the previous event in 1746 was much larger than any of the three underthrusting earthquakes this century the 1746 earthquake may have ruptured the entire segment in a multiple asperity earthquake thus the mode of rupture along the peru coast has changed between successive earthquake cycles

seismicity associated with the great colombia ecuador earthquakes of 1942 1958 and 1979 implications for barrier models of earthquake rupture abstract we have relocated earthquakes along the colombia ecuador coast with the method of joint hypocenter determination in order to examine the nature of fault heterogeneities that controlled northward propagation of plate boundary rupture from the source region of the earthquake of 14 may 1942 to the source region of the 19 january 1958 earthquake and eventually to the source region of the earthquake of 12 december 1979 the entire region of the 1942 1958 and 1979 shocks had previously ruptured in the single great earthquake of 31 january 1906 the relocated hypocenters lie on the same plane to within the approximately 20 km uncertainty of the focal depths also the main shocks nucleated at nearly the same distance from the colombia trench these observations suggest that the heterogeneities between the 1942 and 1958 ruptures and between the 1958 and 1979 ruptures do not correspond to a major distortion of the downgoing slab but rather to either minor distortions of the slab or to regions of high friction or low available strain energy on a continuous fault surface in particular the heterogeneity between the 1958 and 1979 rupture zones seems to have been a high strength barrier with dimensions much smaller than the dimensions of either of the rupture zones the 1942 and 1958 earthquakes had source dimensions no larger than the 1979 main shock but they had stronger aftershock sequences than the 1979 earthquake this observation suggests that the arrest of earthquake rupture in 1979 left the plate boundary segment that had ruptured in 1906 in a state of lower stress than it had been following the 1942 and 1958 earthquakes long term seismicity in the decades preceding the 1979 earthquake occurred mostly outside or on the boundaries of the rupture area defined by the distribution of 1979 aftershocks the intense aftershock activity that followed the 1958 main shock within tens of kilometers of the eventual 1979 hypocenter may correspond to a long term precursory seismic swarm for the 1979 earthquake

geodetic teleseismic and strong motion constraints on slip from recent southern peru subduction zone earthquakes we use seismic and geodetic data both jointly and separately to constrain coseismic slip from the 12 november 1996 m w 7 7 and 23 june 2001 m w 8 5 southern peru subduction zone earthquakes as well as two large aftershocks following the 2001 earthquake on 26 june and 7 july 2001 we use all available data in our inversions gps interferometric synthetic aperture radar insar from the ers 1 ers 2 jers and radarsat 1 satellites and seismic data from teleseismic and strong motion stations our two dimensional slip models derived from only teleseismic body waves from south american subduction zone earthquakes with m w 7 5 do not reliably predict available geodetic data in particular we find significant differences in the distribution of slip for the 2001 earthquake from models that use only seismic teleseismic and two strong motion stations or geodetic insar and gps data the differences might be related to postseismic deformation or more likely the different sensitivities of the teleseismic and geodetic data to coseismic rupture properties the earthquakes studied here follow the pattern of earthquake directivity along the coast of western south america north of 5 s earthquakes rupture to the north south of about 12 s directivity is southerly and in between earthquakes are bilateral the predicted deformation at the arequipa gps station from the seismic only slip model for the 7 july 2001 aftershock is not consistent with significant preseismic motion

seismic rupture associated with subduction of the cocos ridge the april 3 1983 costa rica earthquake mw 7 4 ms 7 3 occurred where the aseismic cocos ridge has been partially subducting beneath the caribbean plate the epicenter is at the eastern edge of the 70 110 km initial aftershock area which later shows considerable expansion p wave directivity analysis indicates that the rupture front propagated westward 63 km from the epicenter to where the major moment release was truncated the source time history of the 1941 event ms 7 5 which occurred in the same epicentral area appears to be similar to the 1983 event but with more complexity and a slower initiation of moment release the low stress drop of the 1983 event 2 3 bar and low ratio of seismic slip to tectonic slip 0 15 for the 1983 event suggest that the coupling between the descending slab and the overlying plate is weak yet a recurrence time for large events mag 7 0 in this zone of about 40 years is implied by the 1904 1941 and 1983 events suggesting that the cocos ridge causes increased partial coupling compared with that in the adjacent subduction zone the slow source rupture process may reflect the nature of the contact area between the costa rica arc and the subducting plate in the vicinity of the bouyant cocos ridge the normal faulting event of july 1 1979 mb 5 3 is consistent with present day uplift of the burica peninsula which previous studies suggest is related to the subduction of the cocos ridge the burica and osa peninsulas may behave as detached blocks which influence the distribution of partial coupling reflected in the aftershock expansion pattern

strain budget of the ecuador colombia subduction zone a stochastic view the 2016 pedernales earthquake m w 7 8 ruptured a portion of the colombia ecuador subduction interface where several large historical earthquakes have been documented since the great 1906 earthquake m 8 6 considering all significant ruptures that occurred in the region it has been suggested that the cumulative moment generated co seismically along this part of the subduction over the last century exceeds the moment deficit accumulated inter seismically since 1906 such an excess challenges simple models with earthquakes resetting the elastic strain accumulated inter seismically in locked asperities these inferences are however associated with large uncertainties that are generally unknown the impact of spatial smoothing constraints on co seismic and inter seismic models also prevents any robust assessment of the strain budget we propose a bayesian kinematic slip model of the 2016 pedernales earthquake using the most comprehensive dataset to date including insar and gps offsets tsunami waveforms and kinematic records from high rate gps and strong motions in addition we use inter seismic geodetic velocities to produce a probabilistic inter seismic coupling model of the subduction interface our stochastic co seismic and inter seismic solutions include the ensemble of all plausible models consistent with our prior information and that fit the observations within uncertainties the analysis of these model ensembles indicates that an excess of co seismic moment during the 1906 2016 period is likely in central ecuador only if we assume that 1942 and 2016 earthquakes are colocated if this assumption is relaxed we show that this conclusion no longer holds given uncertainties in co and inter seismic processes the comparison of 1942 and 2016 teleseismic records reveals large uncertainties in the location of the 1942 event hampering our ability to draw strong conclusions on the unbalanced moment budget in the region our results also show a heterogeneous coupling of the subduction interface that coincides with two slip asperities in our co seismic model for the 2016 pedernales earthquake and with the location of historical ruptures in 1958 1979 and 1998 the spatial variability in coupling and complexity in earthquake history suggest strong heterogeneities in frictional properties of the subduction megathrust

are rupture zone limits of great subduction earthquakes controlled by upper plate structures evidence from multichannel seismic reflection data acquired across the northern ecuador southwest colombia margin 1 subduction of the nazca plate beneath the ecuador colombia margin has produced four megathrust earthquakes during the last century the 500 km long rupture zone of the 1906 mw 8 8 event was partially reactivated by three thrust events in 1942 mw 7 8 1958 mw 7 7 and 1979 mw 8 2 whose rupture zones abut one another multichannel seismic reflection and bathymetric data acquired during the sisteur cruise show evidence that the margin wedge is segmented by transverse crustal faults that potentially correlate with the limits of the earthquake coseismic slip zones the paleogene neogene jama quininde and esmeraldas crustal faults define a 200 km long margin crustal block that coincides with the 1942 earthquake rupture zone subduction of the buoyant carnegie ridge is inferred to partially lock the plate interface along central ecuador however coseismic slip during the 1942 and 1906 earthquakes may have terminated against the subducted northern flank of the ridge we report on a newly identified manglares crustal fault that cuts transversally through the margin wedge and correlates with the limit between the 1958 and 1979 rupture zones during the earthquake cycle the fault is associated with high stress concentration on the plate interface an outer basement high which bounds the margin seaward of the 1958 rupture zone may act as a deformable buttress to seaward propagation of coseismic slip along a megathrust splay fault coseismic uplift of the basement high is interpreted as the cause for the 1958 tsunami we propose a model of weak transverse faults which reduce coupling between adjacent margin segments together with a splay fault and an asperity along the plate interface as controlling the seismogenic rupture of the 1958 earthquake

thermal regime from bottom simulating reflectors along the north ecuador south colombia margin relation to margin segmentation and great subduction earthquakes 1 the north ecuador south colombia nesc margin has three transverse morphotectonic segments the manglares tumaco and patia segments each with different tectonic and structural patterns following the 1906 subduction earthquake the nesc margin has been the site of three megathrust events in 1942 1958 and 1979 for which the rupture zones abut one another we first investigated variations in heat flow derived from bottom simulating reflectors bsr observed along multichannel seismic lines sampling the three morphotectonic segments strong along strike variations of the bsr derived heat flow just landward of the deformation front suggest that each morphotectonic segment has a specific thermal regime finite element thermal models show that these variations are mainly produced by changes in the age and the dip of the oceanic plate and local hydrothermal cooling we then examined the relationship between seismogenesis and thermal structure along the plate boundary the updip limits of the 1942 and 1979 seismogenic zones estimated from the aftershock area and rupture zone reach the trench where the temperature at the top of the subducting plate appears to be 50 60 c between these events the seismogenic zone of the 1958 earthquake is restricted to a region landward of a prominent outer basement high where the temperature appears to be 100 120 c we propose that on the nesc margin the updip limit of the seismogenic zone is primarily controlled by low temperature processes except for the 1958 event for which the seismogenic updip limit appears to be related to a structural feature in the upper plate

magmatic gradients in the cretaceous caleu pluton central chile injections of pulses from a stratified magma reservoir the caleu pluton central chile extending over 338 km2 and with more than 1 400 m of vertical relief intrudes the n s trending lower cretaceous volcano sedimentary and volcanic successions at a depth equivalent to a pressure of 2 kb the host stratified volcanic successions are tilted about 30 40 e whereas the pluton shows paleomagnetic evidence of either tilting of 15 e or clockwise rotation by few degrees a gradient of westward increasing sio2 content is recognized across the pluton giving rise to three n s elongated zones gabbro diorite zone gdz tonalite zone tz and granodiorite zone gz biotite and hornblende compositions also exhibit a westward decreasing gradient in mg mg fe indicating that the more mafic the zone is the more oxidizing is its crystallization condition horizontal inward gradients of progressively less evolved rocks are recognized across gdz and tz whereas no horizontal gradients were found in the gz vertical compositional gradients are recognized in the gdz and tz which consist of an upward increase in sio2 and decrease in mgo feo fe2o3 and compatible trace elements a vertical compositional boundary was recognized along a traverse across the tz separating two magma pulses with similar trends of compositional variations the three zones of the caleu pluton were derived from a common isotopically sr nd depleted source each zone probably evolved independently as their compositional characteristics would have not been acquired in situ the resulting compositional characteristics of the zones would have been developed prior to the intrusion in a subjacent stratified reservoir placed at about seven kilometers below the pluton

rupture process of the february 4 1965 rat islands earthquake the great rat islands underthrusting earthquake mw 8 7 of february 4 1965 represents subduction of the pacific plate beneath the north american plate along a 600 km segment of the western end of the aleutian islands body wave inversion techniques are used to determine the spatial and temporal heterogeneities associated with the rat islands earthquake we have deconvolved world wide standard seismograph network long period teleseismic p wave seismograms to obtain source time functions directivity associated with the three major pulses of moment release in the source time functions indicates a total source duration of 160 s unilateral rupture in the direction 300 fault length of 420 km and average rupture velocity of 2 5 km s the three pulses of moment release are located along the fault and these regions of high moment release are interpreted as asperities the first asperity extends from the epicenter to 100 km to the wnw this is the largest asperity and corresponds to a smooth pulse of moment release in the source time function with a duration of 50 s the second pulse of moment release is very jagged is less coherent between stations and is centered 200 km wnw of the epicenter the third pulse of moment release extends from 360 to 420 km wnw of the epicenter although the aftershock area is 600 km in length we can not resolve any moment release from the p waves beyond 420 km wnw of the epicenter the rat islands event was closely followed by a large tensional outer rise event on march 30 1965 ms 7 5 which is located oceanward of the largest moment release associated with the rat islands mainshock rupture detailed analysis of the p waves for this large outer rise event confine the depth extent to the upper 30 5 km of the crust the spatial and temporal association between the february 4 mainshock and the march 30 tensional outer rise event suggests the tensional event may have been triggered by the large displacement near the mainshock epicenter the overriding plate along the western aleutian subduction zone is laterally segmented into a series of rigid tectonic blocks separated by fault controlled canyons and extensional basins geist et al 1988 we suggest that the central undeformed parts of the blocks have the strongest coupling with the down going plate and hence are the sites of the largest moment release during an underthrusting earthquake the three asperities determined from the p waves correspond to the rat buldir and near tectonic blocks respectively hence the p wave seismic moment release of the rat islands earthquake is controlled by the lateral segmentation of the overriding plate

asperity distribution and earthquake occurrence in the southern kurile islands arc the seismicity in the southern kurile islands arc is characterized by a complex overlapping of aftershock zones of a few great and many large earthquakes ms 7 0 the spatial relationships of these aftershock zones can be explained in terms of a simple pattern of fault plane heterogeneity the four largest earthquakes in this region are from north to south the 1963 kurile islands mw 8 5 the 1958 etorofu mw 8 3 the 1969 kurile islands mw 8 2 and the 1973 nemuro oki mw 7 8 events detailed body wave studies indicate that the moment release associated with each event is concentrated in a small portion of the aftershock area with scale lengths 100 km these regions of high moment release are the largest asperities the areas that surround the large asperities slip in two modes 1 coseismically during the great events and 2 as distinct large ms 7 0 earthquakes both before and after the great events the second mode involves either rupture of relatively weak regions which results in low stress drop events or rupture of smaller asperities the earthquakes produced by failure of small asperities tend to repeatedly rupture the same regions and they have impulsive source time functions indicating rupture of highly stressed regions of the fault surface by analysis of all ms 7 0 events since 1946 in this region we find a segmentation of the southern kurile islands thrust zone into strongly coupled regions with two scale lengths dominant and secondary asperities and weakly coupled regions the surrounding areas generating the low stress drop events the largest earthquakes and the smaller asperity earthquakes occur at the down dip portion of the lithospheric interface while the low stress drop events occur much closer to the trench axis this implies a difference in the degree of mechanical coupling between the subducting and overriding plates with the deeper portions having greater strength than the shallower parts

distribution of discrete seismic asperities and aseismic slip along the ecuadorian megathrust abstract a dense gps network deployed in ecuador reveals a highly heterogeneous pattern of interseismic coupling confined in the first 35 km depth of the contact between the subducting oceanic nazca plate and the north andean sliver interseismic models indicate that the coupling is weak and very shallow 0 15 km in south ecuador and increases northward with maximum found in the rupture areas of large mw 7 0 megathrust earthquakes that occurred during the 20th century since the great 1906 mw 8 8 colombia ecuador earthquake may have involved the simultaneous rupture of three to six asperities only one or two asperities were reactivated during the large seismic sequence of 1942 mw 7 8 1958 mw 7 7 1979 mw 8 2 and 1998 mw 7 1 the axis of the carnegie ridge which is entering the subduction zone south of the equator coincides well with the location of a 50 km wide creeping corridor that may have acted as persistent barrier to large seismic ruptures south of this creeping region a highly locked asperity is found right below la plata island while this asperity may have the potential to generate an mw 7 0 7 5 earthquake and a local tsunami until now it is unknown to have produced any similar events that region is characterized by the presence of slow slip events that may contribute significantly to reduce the long term moment deficit accumulated there and postpone the failure of that asperity at the actual accumulation rate a characteristic recurrence time for events such as those in 1942 1958 and 1979 is 140 30 yr 90 20 yr 153 80 yr respectively for the great 1906 event we find a recurrence time of at least 575 100 yr making the great 1906 earthquake a rare super cycle event

aftershock area expansion and mechanical heterogeneity of fault zone within subduction zones we investigate global variations in the pattern of aftershock area expansion associated with large subduction zone earthquakes the expansion of aftershock areas with time is relatively small in the subduction zones of alaska the aleutians mexico and parts of the kuriles and south america these subduction zones are interpreted to be characterized by moderate to strong interplate coupling in contrast the extent of aftershock expansion is greater for events within the subduction zones of northeast japan and western and south pacific where interplate coupling is inferred to be weak these results are interpreted in terms of an asperity model where a fault zone is represented by a distribution of strong spots asperities and weaker zones the immediate aftershock area is assumed to extend over an area occupied by asperities and the subsequent expansion to occur into the surrounding weaker zones if large asperities abut each other on a fault plane interplate coupling is strong and little expansion of aftershock area occurs smaller and more sparsely distributed asperities cause weak interplate coupling and significant expansion of aftershock area small and densely distributed asperities result in moderate interplate coupling and little aftershock expansion

the rupture process of the great 1979 colombia earthquake evidence for the asperity model we have determined the asperity distribution in the segment of the ecuador colombia subduction zone that was ruptured by the 1979 mw 8 2 and 1958 mw 7 7 earthquakes for the 1979 earthquake we are able to document directivity in the p wave source time functions and consequently have spatially located the moment release by combining results from both p wave and surface wave directivity we have defined three regions 1 from 0 to 56 km northeast of the epicenter with low moment release and a small seismic displacement 2 the asperity from 56 to 116 km northeast of the epicenter with most of the moment release and a displacement of 6 m and 3 a small displacement region from 116 to 180 240 km northeast of the epicenter thus the 1979 rupture zone is composed of an asperity 60 km in length bounded by two weak small displacement regions the displacement at the asperity is equal to the cumulative tectonic displacement between 1979 and 1906 the year of the previous great earthquake in this segment the 1958 earthquake is characterized by the rupture of a small asperity length scale approximately 25 km in the epicentral region the 1958 event continued to rupture past this small epicentral asperity and stopped near the epicenter of the following 1979 event this detailed study of the spatial variation in moment release offers strong support for the asperity model of large earthquake occurrence

do trench sediments affect great earthquake occurrence in subduction zones seismic energy release is dominated by the underthrusting earthquakes in subduction zones and this energy release is further concentrated in a few subduction zones while some subduction zones are characterized by the occurrence of great earthquakes others are relatively aseismic this variation in maximum earthquake size between subduction zones is one of the most important features of global seismicity previous work has shown that the variation in maximum earthquake size is correlated with the variation in two other subduction zone properties age of the subducting lithosphere and convergence rate these two properties do not explain all the variance in maximum earthquake size i propose that a third subduction zone property trench sediments explains part of the remaining variance in maximum earthquake size subduction zones are divided into two groups 1 those with excess trench sediments and 2 those with horst and graben structure at the trench thirteen of the 19 largest subduction zone events including the three largest occur in zones with excess trench sediments about half the zones with excess trench sediments are characterized by great earthquake occurrence most of the other zones with excess trench sediments but without great earthquakes are predicted to have small earthquakes by the age rate correlation two notable exceptions are the oregon washington and middle america zones overall the presence of excess trench sediments appears to enhance great earthquake occurrence one speculative physical mechanism that connects trench sediments and earthquake size is that excess trench sediments are associated with the subduction of a coherent sedimentary layer which at elevated temperature and pressure forms a homogeneous and strong contact zone between the plates

earthquakes of the holocene geologic studies of earthquakes involve mapping of coseismic surface faulting and secondary deformation from historical events trenching and geomorphic analysis to define the timing and size of past earthquakes and investigations of fault zone structure and geometry in both unconsolidated sediments and bedrock this research is now being referred to as paleoseismology seismic geology and earthquake geology since the mid 1970s it has led to some of the most exciting and important contributions to the understanding of earthquake behavior in space and time hanks 1985 allen 1986 the present report is the first of what is hoped will become a regularly contributed summary of the geologic aspects of the study of earthquakes

seamount subduction at the north ecuadorian convergent margin effects on structures inter seismic coupling and seismogenesis at the north ecuadorian convergent margin 1 s 1 5 n the subduction of the rough nazca oceanic plate leads to tectonic erosion of the upper plate and complex seismogenic behavior of the megathrust we used three selected pre stack depth migrated multi channel seismic reflection lines collected during the sisteur cruise to investigate the margin structure and decipher the impact of the subducted atacames seamounts on tectonic erosion interseismic coupling and seismogenesis in the region of the 1942 mw7 8 earthquake this dataset highlights a subducted 30 40 km 30 40 km double peak seamount that belongs to the atacames seamount chain and that is associated with a deep morphologic re entrant containing mass transport deposits the seamount subduction uplifted the margin basement by 1 6 km and pervasively broke the margin by deep and intense reverse faulting ahead of the seamount a process that is likely to weaken considerably the margin in the seamount wake the basement reverse fault system rotated counter clockwise this faulted basement is overlain with slope sediment sliding along listric normal faults that sole out onto the bsr this superposition of deep tectonic contraction within the basement and shallow gravitational extension deformation within the sediment highlights the key role of gas hydrate on outer slope erosion in addition to long term regional basal erosion the margin basement has thinned locally by an extra 0 8 1 km in response to the subduction of the atacames seamount chain and hydrofracturing by overpressured fluids at the margin toe this pervasively and deeply fractured margin segment is associated with a seismically quiet and gps modeled low interseismic coupling corridor that terminates downdip near the 1942 epicenter and locked zone we suggest that the deeply buried double peak atacames seamount triggered the 1942 earthquake ahead of its leading flank this result supports previous studies proposing that subducted seamounts provide unfavorable conditions for locking the updip segment of the plate boundary limiting the updip extent of seismogenic zones but may favor large subduction earthquakes at greater depths