a body wave analysis of the 1966 gisborne new zealand earthquake abstract the ml 6 2 gisborne earthquake of march 4 1966 occurred along the east coast of the north island new zealand modeling of p and s body waves shows that the focal mechanism of this event is consistent with northwestward thrusting of the pacific plate beneath the north island 249 25 and 131 the focal depth is constrained to 18 km significantly less than the values of 25 30 km computed from local network data estimates of the scalar moment source duration and stress drop for the event are 4 1024 dyne cm 2 3 s and 20 120 bar respectively cross correlation errors of synthetic to observed waveforms were computed for all possible p and t axis locations and slip vector orientations and contoured on a projection of the focal sphere the error contour at which the synthetic waveforms distinctly diverged from the observed waveforms was established by eye the procedure shows that the analysis of long period body waves at least in this case provides much better constraint on focal mechanism orientation than does first motion data alone

the hawke s bay new zealand earthquake of 1982 september 02 and deformation in the interior of the subducted pacific plate abstract the most damaging earthquake to occur in new zealand in 1982 was the hawke s bay earthquake of september 02 m b 5 6 an aftershock study using portable microearthquake recorders has established that this event occurred some 20 km south of hastings at a depth of about 47 km this depth places the event in the mantle of the subducted pacific plate this earthquake is a member of a recently recognized class of intraplate shocks occurring along the strike of the subducted plate near the east coast of the north island such events are characterized by normal faulting mechanisms when taken together with activity in the crust of the subducted plate this subcrustal activity resembles a double planed seismic zone restricted in the down dip direction of the subducted plate several lines of evidence suggest that deformation in the mantle of the subducted plate is controlled by the state of coupling of the plate interface at shallow depth 1 1 the subcrustal activity occurs directly beneath the source region of major interplate earthquakes along the east coast of the north island and there appears to be a pairing of interplate and subcrustal events 2 2 stress directions determined for the subcrustal events are concordant with geodetic strain directions and 3 3 the amount of subcrustal activity appears to be related to the state of coupling of the plate interface

mechanics of fold and thrust belts and accretionary wedges cohesive coulomb theory a critically tapered fold and thrust belt or submarine accretionary wedge is one that is on the verge of coulomb failure everywhere including its base where frictional sliding along a decollement is assumed to be occurring cohesion within a wedge can add significantly to the overall strength near the toe the effect of this is to decrease the near toe taper leading to a critical topographic profile that is concave upward if the decollement is planar we obtain an approximate self consistent solution for the state of stress within a thin skinned cohesive critical coulomb wedge and determine the relationship between the wedge taper and its strength and basal friction the theory is then applied to the presently deforming fold and thrust belt of western taiwan fitting of theoretical critical wedge shapes to topographic profiles and measurements of the step up angles of thrust faults from the basal decollement are used to constrain the taiwan wedge strength parameters an attractive assertion fully consistent with all the observations is that the mechanics of fold and thrust belts and accretionary wedges is governed by normal frictional and fracture strengths of rocks measured in the laboratory in particular if byerlee s law b 0 85 is adopted as the coefficient of sliding friction on the base we find a coefficient of internal friction 0 9 1 0 in the wedge and a wedge cohesion so 5 20 mpa other solutions having strengths and ambient stresses up to 4 times lower than this can also however satisfy the data

coseismic uplift of holocene marine terraces in the pakarae river area eastern north island new zealand abstract holocene marine terraces along 15 km of the northeastern coast of north island record episodic tectonic uplift a maximum of seven terraces are arranged in staircase fashion and lie about 20 km above the subduction interface between the pacific and australian plates the highest t1 corresponds with the maximum of the holocene marine transgression about 6700 14 c yr b p younger terraces are marine abrasion platforms overlain by thin beach deposits radiocarbon ages of marine shells from beach deposits indicate that uplift above marine conditions occurred ca 6700 t1 5500 t2 3900 t3 2500 4 1600 t5 1000 t6 and slightly younger than 600 t7 yr b p uplift probably occurred coseismically the maximum late holocene uplift rate in the study area is 8 mm yr altitudinal distribution of terraces suggests deformation exists as a ca 20 km elongate dome broken at the southern end by the pakarae fault which trends across the dome rupture on this fault has accompanied the growth of the dome but is not responsible for it bathymetric profiling suggests that an active fault parallel to and about 5 km offshore is probably responsible for the episodic coastal uplift

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

geodetically determined strain across the southern end of the tonga kermadec hikurangi subduction zone summary a method of simultaneous reduction is presented for determining strain rates from multiple triangulation surveys where common triangulation stations have been used but the angles of the old survey have not necessarily been reobserved this method is applied to triangulation in the northern south island at the southern end of the tonga kermadec hikurangi subduction zone from a profile of shear strain across the indian pacific plate boundary the displacement of the indian plate relative to the pacific is calculated to be 54 9 mm yr at an azimuth of 84 f lo in remarkable agreement with the motion predicted by global plate tectonic models most of this motion occurs within a 15okm wide zone bounded on the east by the hikurangi trough within this zone the motion is partitioned near the hikurangi trough no slip is occurring at the upper surface of the subducting pacific plate the subduction thrust and motion is predominantly thrusting normal to the trough axis to the west is a region of predominantly dextral strike slip faulting this pattern is consistent with fitch s model of oblique subduction to the south of the profile a change is observed in the azimuth of the faulting along a line which marks the southern extent of the subduction slab indicating the end of the partitioned motion

normal faulting through subducted oceanic crust the 19 july 1985 earthquake of hawke s bay new zealand abstract source parameters of the m l 5 7 19 july 1985 earthquake of northern hawke s bay new zealand are determined by modelling the teleseismic p and s body waves a focal depth of 31 1 km scalar moment 5 6 10 17 nm and source duration of 3 0 3 5 s are estimated the focal mechanism strike 224 dip 78 rake 273 shows normal faulting on a nw dipping plane that is parallel to the strike of the subducted pacific plate thirty five aftershocks were recorded by a temporary seismograph network 21 of which were relocated using the jhd method the hypocentres define a northwest dipping plane 20 30 km which extends from the hypocentre of the mainshock to a depth of 47 km the position of the subducted plate in hawke s bay is well defined from previous micro earthquake studies in the region of the mainshock the plate interface is located at a depth of approximately 28 to 30 km placing the hypocentre of the mainshock in the crust of the pacific plate the extent and location of the aftershock distribution suggest that the mainshock ruptured downward and extended through the crust of the plate this event thus differs from other normal faulting events in the region that usually occur in the mantle of the subducted plate the scalar moment at long periods t 45s dziewonski et al 1986 exceeds that determined here by a factor of 2 suggesting a source rich in long period energy an idea supported by the relatively large aftershock area

the eastern khoy metamorphic complex of nw iran a jurassic ophiolite or continuation of the sanandaj sirjan zone the khoy complex in nw iran has been widely regarded to consist of both jurassic and cretaceous ophiolites but whereas the western late cretaceous ophiolite is unequivocal the nature of the eastern jurassic body has been unclear field observations show that the presumed eastern meta ophiolite has no similarities to an ophiolite sequence here we report geological and zircon rutile titanite u pb data for rocks from the presumed jurassic ophiolite obtained to verify whether these are real ophiolitic units and to understand their relations either to the late cretaceous zagros bitlis ophiolites in western iran southern anatolia or to the sevan akera northern armenia and izmir ankara southern pontides complexes the new u pb ages show that the presumed ophiolite is in fact a collage of ediacaran to cambrian c 606 517 ma and jurassic c 160 ma meta igneous rocks similar to ages obtained for igneous rocks of the sanandaj sirjan zone which has been interpreted as a jurassic continental rift the jurassic igneous rocks contain abundant ediacaran ordovician silurian and carboniferous permian inherited zircons further suggesting involvement of pre existing continental crust our results indicate that jurassic continental rifting provided a lithospheric weakness along which a new subduction zone formed in late cretaceous times

slip partitioning along major convergent plate boundaries along plate boundaries characterized by oblique convergence earthquake slip vectors are commonly rotated toward the normal of the trench with respect to predicted plate motion vectors consequently relative plate motion along such convergent margins must be partitioned between displacements along the thrust plate interface and deformation within the forearc and back arc regions the deformation behind the trench may take the form of strike slip motion back arc extension or some combination of both we observe from our analysis of the harvard moment tensor catalog that convergent arcs characterized by back arc spreading specifically the marianas and new hebrides are characterized by a large degree of slip partitioning however the observed rates directions and location of back arc spreading are not sufficient to account for degree of partitioning observed along the respective arcs implying that the oblique component of subduction is also accommodated in part by shearing of the overriding plate in the case of the sumatran arc where partitioning is accommodated by strike slip faulting in the overriding plate the degree of partitioning is similar to that observed along the marianas but the result is viewed with caution because it is based on a predicted plate motion vector that is based on locally derived earthquake slip vectors in the case of the alaskan aleutian arc where back arc spreading is also absent the degree of partitioning is less and rotation of slip vectors toward the trench normal appears to increase linearly as a function of the obliquity of convergence if partitioning in the alaskan aleutian arc is accommodated by strike slip faulting within the upper plate the positive relationship between obliquity of convergence and the rotation of earthquake slip vectors to the trench normal may reflect that either 1 the ratio of the depth extent of strike slip faults behind the trench z s to the subduction thrust z t increases westward along the arc 2 the dip of the subduction thrust increases westward along the arc or 3 the strength of the subduction thrust decreases westward along the arc

slab pull and the seismotectonics of subducting lithosphere this synthesis links many seismic and tectonic processes at subduction zones including great subduction earthquakes to the sinking of subducted plate earthquake data and tectonic modeling for subduction zones indicate that the slab pull force is much larger than the ridge push force interactions between the forces that drive and resist plate motions cause spatially and temporally localized stresses that lead to characteristic earthquake activity providing details on how subduction occurs compression is localized across a locked interface thrust zone because both the ridge push and the slab pull forces are resisted there the slab pull force increases with increasing plate age thus because the slab pull force tends to bend subducted plate downward and decrease the force acting normal to the interface thrust zone the characteristic maximum earthquake at a given interface thrust zone is inversely related to the age of the subducted plate the 1960 chile earthquake mw 9 5 the largest earthquake to occur in historic times began its rupture at an interface bounding oceanic plate 30 m y old however this rupture initiation was associated with the locally oldest subducting lithosphere weakest coupling the rupture propagated southward along an interface bounding progressively younger oceanic lithosphere terminating near the subducting chile rise prior to a great subduction earthquake the sinking subducted slab will cause increased tension at depths of 50 200 km with greatest tension near the shallow zone resisting plate subduction plate sinking not only leads to compressional stresses at a locked interface thrust zone but may load compressional stresses at plate depths of 260 350 km provided that the shallow sinking occurs faster than the relaxation time of the deeper mantle this explains k mogi s observations of m 7 thrust earthquakes at depths of 260 350 km immediately downdip and within 3 years prior to five great shallow earthquakes of northern japan the slab pull model explains the lower layer of double seismic zones as due to tension from the deeper sinking plate and the upper layer as due to localized in plate compression as plate motion is resisted by the bounding mantle just downdip of the interface thrust zone there occurs an aseismic 20 50 dip increase of subducted plate this slab bend reflects the summed slab pull force of deeper plate and probably is at the crustal basalt to eclogite phase change resistance to subduction provided by a continually developing slab bend may be an important factor in the size of slab pull force delivered to an interface thrust zone

global variability in subduction thrust zone forearc systems deviations of slip vector azimuths of interplate thrust earthquakes from expected plate convergence directions at oblique subduction zones provide kinematic information about the deformation of forearcs and indirect evidence on the dynamics of the plate boundary a global survey of slip vectors at major trenches of the world reveals a large variability in the kinematic response of forearcs to shear produced by oblique convergence the variability in forearc deformation inferred from slip vector deflections is suggested to be caused by variations in forearc rheology rather than in the stresses acting on subduction zone thrust faults estimated apparent macroscopic rheologies range from elastic to perfectly plastic or viscous forearc rheologies inferred from slip vectors do not correlate with age of the subducting lithosphere but continental forearcs or old arcs appear to deform less than oceanic or young arcs the inferred absence of forearc deformation at continental arcs from this study is counter to inferences drawn from compiled geologic information on forearc faults correlations of the apparent forearc rheology with backarc spreading convergence rate slab dip arc curvature and downdip length of the thrust contact are poor however great subduction zone earthquakes occur where forearcs are apparently more elastic i e less deformed by oblique convergence which suggests that the mechanical properties of forearcs rather than stress magnitude on thrust faults control both the kinematic behavior of forearcs and where great subduction zone earthquakes occur

shallow morphology of the subducted pacific plate along the hikurangi margin new zealand abstract along the east coast of the north island new zealand the pacific plate is subducted beneath the australian plate the seismicity associated with this subduction has been determined using the results from a number of temporary and permanent microearthquake networks a synthesis of the microearthquake seismicity the depths of thrust events the depth of s to p conversions and deep reflection seismic studies leads to a coherent picture of the subducting slab both in terms of its mechanical structure and in terms of the general morphology the slab is initially subducted at a low angle of a few degrees is bent through a circular arc of few hundred kilometre radius and then becomes planar again below 120 km depth the radius of curvature varies along strike from about 280 km just north of wellington to 240 km north of hawke bay the slab shape between depths of 15 45 km is thus approximately conical it also approximates within the errors of hypocentre location a surface of constant gaussian curvature as would be geometrically expected from the simple bending of a thin spherical lithospheric shell observed seismicity is primarily concentrated within what is interpreted as the subducted crust of the pacific plate with a lower level of activity within the subducted mantle normal fault mechanisms at the top of the plate probably reflect tensional stress associated with the plate bending whereas thrust events observed at the lower bound of the observed seismicity probably reflect compression below the neutral plane

focal mechanisms of large earthquakes in the north island of new zealand slip partitioning at an oblique active margin we have used body wave modelling to determine the source parameters of 22 moderate to large earthquakes that have occurred along the hikurangi subduction margin and elsewhere in the north island of new zealand since 1964 we have also included events from the harvard cmt catalogue since 1977 as that catalogue contains smaller events than can be modelled with body waves we have found that shallow earthquakes occurring in the back arc taupo volcanic zone and its extension to the north show predominantly normal faulting with nodal planes parallel to the regional fabric and t axes indicating extension in a direction in agreement with geodetic measurements a normal faulting solution for the 1974 opunake earthquake is compatible with the mapping of active faults in the offshore taranaki region and the initiation depth of 17 km is close to the expected brittle ductile transition for that area earthquakes within the subducted plate at the southwestern end of the margin are dominated by down dip tension that is rotated towards the deepest part of the subducting plate further north the events in the upper part of the wadati benioff zone 200 km show pure down dip tension while this association is less obvious for the deeper events suggesting that the aseismic part of the subducting slab inferred from plate reconstructions has reached the 670 km discontinuity normal faulting earthquakes associated with plate bending below the region of interplate contact occur along the length of the margin while shallow normal faulting events also occur in the trench indicating that some slab pull propagates to that region we infer that large interplate earthquakes are not imminent because in other similar areas they are preceded by outer rise compressional events that have not been recorded in the region in historical times slip vectors of interplate thrust earthquakes indicate that the oblique plate motion across the margin is fully partitioned that is the plate interface accommodates very little transcurrent motion which is instead accommodated by faults in the overlying australian plate extensional geodetic strain has been measured across the margin in the northeast which has been regarded as unusual in view of the plate convergence we suggest that the topography at this part of the margin is too steep to be supported by the current horizontal forces due to friction on the plate interface this may have occurred due to a reduction in friction because of the subduction of seamounts and the subsequent entrapment of sediments this also being consistent with the idea of underplating at this part of the margin or due to the subduction of less buoyant crust other margins such as japan and mid america where there is known to be subduction of sediments also show extension in the forearc compressional strain observed across the margin in the southwest is more usual for a convergent margin but it does not show the partitioning we see in focal mechanisms these data might only agree after the long term compressional strain accumulation has been averaged over a complete seismic cycle large historical and palaeo earthquakes have occurred along the length of the margin associated with deformation in both the axial and the coastal ranges in addition along the length of the margin the rate of recent interplate thrust activity is low compared to margins that have a high proportion of aseismic slip these factors indicate that the entire length of this margin is capable of producing large earthquakes

terrestrial and gps measurements of deformation across the taupo back arc and hikurangi forearc regions in new zealand geodetic measurements of deformation across the taupo back arc and hikurangi forearc regions of new zealand are derived from global positioning system gps measurements made in 1990 and 1991 and triangulation observations made in the 1920s 1950s and 1970s the gps horizontal coordinate differences have precisions of 5 6 mm and the triangulation observations have precisions of 0 7 1 2 arc sec these different kinds of observation allow simultaneous estimation of strain parameters and single epoch coordinates for a total of 184 stations at 68 of which gps observations were made under the assumption that no shift of the scale or orientation of the gps reference frame is embodied in the data from 13 gps stations within the taupo volcanic zone tvz common to the 1990 and 1991 surveys the rotational dilatational and shear components of the deformation rate tensor within the tvz can be estimated for that 1 year interval the principal extension rate from 1990 to 1991 is 0 21 0 09 10 6 yr 68 confidence at an azimuth of 124 13 corresponding to 8 4 mm yr extension over 40 km neither the orthogonal principal extension rate 0 02 0 03 10 6 yr nor the rotation rate 0 06 0 04 rad yr is significant the dilatation rate of 0 23 0 09 10 6 yr is therefore produced by uniaxial extension the data from 75 stations distributed from behind the back arc region and well into the forearc region which were surveyed in at least two of the triangulation or gps epochs between the 1920s and 1991 allow the spatial variation of the shear components of the deformation rate tensor to be estimated the maximum engineering shear rates lie in the range 0 1 0 2 10 6 yr within the tvz and maintain these values with similar orientations in the forearc region to the east the western and southern margins of active extension are reasonably well determined by changes in orientation and magnitude of the shear component these results provide confirmation of the previously less well determined deformation field and in comparison with subduction models indicate that the surface deformation is reflecting variation from north to south of coupling of the subduction plate interface and that both trench suction and gravitational collapse probably contribute to the extension in the forearc region

body waveforms and source parameters of some moderate sized earthquakes near north island new zealand p wave first motion and synthetic seismogram analysis of p and sh waveforms recorded at teleseismic distances on the wwssn are used to estimate source parameters of seven of the largest earthquakes 6 1 mb 6 3 that occurred in the vicinity of north island new zealand since 1965 the source parameters of three other mb 6 1 events determined outside of this study are included and considered in the final analysis four of the earthquakes occurred at shallow depths 20 km of which three were located within and to the north of north island two of the shallow events show strike slip and normal focal mechanisms with t axes oriented in a manner consistent with their location in an area of known back arc extension one of the shallow events occurred in northern south island and shows a reverse type mechanism indicating horizontal contraction of the crust in an easterly azimuth six events occurred at intermediate depths h 39 to 195 km of which five exhibit thrust mechanisms with t axes consistently oriented near vertical in the light of previously published plate tectonic models the near vertical orientation of t axes of the intermediate depth events may be used to infer that the southern kermadec plate boundary immediately north of north island is not strongly coupled and hence not likely capable of producing great earthquakes a similar inference cannot be made for the section of the hikurangi margin adjacent to north island since the intermediate depth events considered in this study lie to the north of this segment of the plate boundary

holocene paleoseismicity in the fold and thrust belt of the hikurangi subduction zone eastern north island new zealand abstract along the 500 km long eastern coastline of the north island new zealand there are at most seven distinct holocene marine terraces aged 7 ka b p or less the highest of the terraces is 27 m above present day mean sea level the coastal region is subdivided into fourteen distinct subregions based on radiocarbon ages of marine deposits overlying wave cut shore platforms and geographic distribution of similarly aged terraces holocene marine terraces are the result of uplift associated with large earthquakes co seismic deformation this conclusion is based on characteristic stepped terrace morphology clustering of ages of terrace deposits within subrogions and the occurrence of co seismic uplift in historic time differential uplift across structures and distinct age variations at subregion boundaries are also characteristic recurrence intervals of uplift in any one terrace sequence vary from 0 4 ka to 2 0 ka and individual amounts of uplift vary from 1 0 m to 4 0 m in the past 2 5 ka ages from all terraces within subrogions indicate at least 21 paleoseismic events affecting coastal areas of eastern north island in that period these events cluster in time and in separate parts of eastern north island several earthquakes occurred 0 3 0 6 1 0 1 5 2 1 and 2 3 ka b p our data strongly support the concept of segmentation of deformation along the subduction margin the likely cause of coastal uplift is movement on steep reverse faults local structures rather than the subduction thrust that propagate from or near the subduction thrust some 20 25 km below the region earthquakes of moment magnitude 7 3 8 0 are estimated to be associated with these fault movements

rheological control of wadati benioff zone seismicity intermediate and deep focus earthquakes in wadati benioff zones are thought to occur in the cold interiors of downgoing slabs which are significantly stronger than the warmer mantle given that earthquakes in oceanic lithosphere appear restricted by an isotherm and hence a given value of lithospheric strength we investigate whether a similar formulation is useful for subducting plates strength in downgoing slabs should be affected by both pressure and temperature an effect previously treated using a depth dependent limiting temperature for seismicity wortel 1982 we find this limiting temperature implies that a possible limiting strength increases strongly with depth unless either the temperatures were too low or the activation volume too large comparison of the analytic model used by wortel with numerical thermal models appears to exclude the first possibility we explore the second possibility by using the numerical thermal model to compute strength contours for flow law constants reported from laboratory experiments and find that the expected pressure strengthening is large enough that the slab should have considerable strength well below the deepest seismicity we conclude that if laboratory results are applicable to these conditions either a strongly depth dependent limiting strength exists or factors in addition to strength control the distribution of subduction zone earthquakes

effects of near source bathymetry on teleseismic p waveforms teleseismic p waveforms from submarine earthquakes show water reverberations that are often poorly modeled by synthetic seismograms this study suggests that considerable misfit results from the non horizontal nature of the seafloor in the epicentral region and examines this effect for a planar dipping water crust interface geometric ray theory synthetics calculated for a double couple source show that p waveforms are extremely sensitive to seafloor orientation because of the large impedance contrast across the boundary seafloor dips as small as 0 5 significantly modify the ray paths affecting the take off angle at the source reflection and transmission coefficients at the interface geometric spreading and travel time to test the applicability of such modeling p waveforms from the march 29 1976 cocos ridge intraplate earthquake were studied the waveforms many of which show prominent water multiples poorly modeled with flat lying bathymetry were inverted for the best fitting planar ocean bottom orientation using a grid search technique synthetics computed for the best fitting seafloor orientation match the data much better than synthetics for a flat lying structure and reproduce the observed azimuthal dependence of water multiple amplitudes the derived ocean bottom dip angle and direction are in agreement with bathymetric charts suggesting p waveforms provide constraints on the seafloor orientation in the epicentral region green s functions incorporating the effects of bathymetry should provide better resolution of the source parameters of submarine earthquakes

loci and maximum size of thrust earthquakes and the mechanics of the shallow region of subduction zones earthquakes do not extend updip along the plate interface at most subduction zones all of the way to the plate boundary at the trench axis or deformation front rather the shallowest part of that interface moves mainly through stable aseismic slip that part of the plate boundary referred to here as the aseismic zone occurs along the base of an accreted wedge of young sediments the probable primary cause for the existence of this aseismic zone is the stable slip properties of the unconsolidated and semiconsolidated sediments in that zone subducted sediment is progressively dewatered and underplated to the base of the overriding plate through this process more consolidated rocks eventually come into contact at depth across both sides of the plate boundary from the point where sufficiently hard rock is found across the plate interface that interface will change its slip behavior to unstable stick slip sliding which is characteristic of consolidated material under most conditions this type of motion is accommodated seismically as episodic slip in large earthquakes the location of this transition to seismic behavior referred to as the seismic front marks the deep end of the aseismic zone and the top of the seismogenic zone i e that part of the plate interface that moves primarily in thrust earthquakes several convergent plate margins are discussed to illustrate the seismic front and the aseismic zone we find that the seismic front defined by smaller earthquakes that occur during the interval between large events is nearly the same as that for large events as inferred from the locations of their aftershocks the location of the seismic front is important for making estimates of the maximum possible size of thrust earthquakes along the plate boundary because it delimits the trenchward limit of the potential rupture area of these events the size of large thrust earthquakes is proportional to rupture area and to at least the cube of the downdip width w of that area thus by defining the location of the seismic front for convergent margins it is possible to make better estimates of w and hence to deduce the maximum size of future interplate earthquakes the average repeat time of such events is also related to w and so better estimates of average repeat time are possible from improved knowledge of w the hypothesis that the location of the seismic front is related to the maximum depth of subduction of unconsolidated sediment has implications for forearc mechanics we use a mechanical analysis laboratory modeling and multichannel seismic information to develop a simple model explaining the growth of forearcs the outer arc high or trench slope break that occurs arcward of most accretionary wedges represents an abrupt change in the critical taper of the accretionary wedge we argue that that change is caused by a large arcward increase in the strength of the material within the overriding plate that stronger material is called the backstop our laboratory modeling experiments indicate that a backstop with a trenchward dipping upper surface results in the development of an overlying structure that includes all of the primary morphologic features observed in modern forearcs in our modeling an accretionary wedge develops trenchward of the backstop an outer arc high develops above the trenchward toe of the backstop and farther arcward a passive forearc basin forms above the stronger material of the backstop this model is consistent with the observation that earthquakes do not extend updip along the plate interface all the way to the trench axis we hypothesize that plate motion is accommodated seismically along the base of the consolidated backstop and mostly or entirely aseismically trenchward of the toe of the backstop along the base of the accretionary wedge at several margins the seismic front is approximately coincident with the outer arc high supporting this interpretation

subduction and back arc activity at the hikurangi convergent margin new zealand the hikurangi margin is a region of oblique subduction with northwest dipping intermediate depth seismicity extending southwest from the kermadec system to about 42 s the current episode of subduction is at least 16 20 ma old the plate convergence rate varies along the margin from about 60 mm a at the south end of the kermadec trench to about 45 mm a at 42 s the age of the pacific lithosphere adjacent to the hikurangi trench is not known

shallow subduction tectonics and fault interaction the weber new zealand earthquake sequence of 1990 1992 a series of four large earthquakes with epicenters very close to one another occurred in 1990 1992 in the weber region of the southeast coast of the north island new zealand the region is one of oblique plate convergence and subduction the plate interface being at about 20 km depth the first event february 19 1990 ml 6 1 mw 6 2 occurred within the upper part of the subducting pacific plate on a steeply northwest dipping normal fault as defined by relocated aftershocks the mechanism of this event reflects downdip tension within the subducting plate due to slab pull forces as is common in new zealand the second event may 13 1990 ml 6 2 mw 6 4 occurred on an imbricate fault in the overlying australian plate dipping shallowly northwest with components of both thrusting and dextral strike slip motion the slip vector is similar to the direction of plate convergence showing no arc normal rotation as observed for subduction thrust events in other areas of oblique convergence it is unlikely that the spatial and temporal closeness of these events is coincidence a similar pair occurred to the south in 1942 but the coupling mechanism is not clear perhaps being related to aseismic slip downdip on the subduction interface two subsequent events both of ml 5 5 occurred on august 15 1990 and march 2 1992 and extended the aftershock zone of the deeper first main shock to the northeast matching the extent of the shallower second main shock mechanisms of aftershocks of the two larger events are quite variable but as a whole are similar to the relevant mainshock the time history of the aftershocks of the first deeper event shows a period of quiescence beginning 35 days before the second shallow event the mean magnitude of those events that did occur was high this observation cannot be explained by station biases weather or changes in analysis methods and may be due to stress changes induced by precursory aseismic slip which would discourage aftershock occurrence

varied responses to subduction in nankai trough and japan trench forearcs responses of sediment entering trenches range from accretion at the foot of the inner slope to subduction into the mantle with inner trench slopes variously displaying uplift subsidence growth and retreat leg 87 is the latest of the dsdp voyages which have provided many clues to the behaviour of subduction zones and was devoted to the study of subduction processes in two trenches off japan the nankai trough fig 1 on the basis of very high quality reflection profiles1 2 is accreting most or all its trench fill which is being stripped off the underlying shikoku basin haemipelagites along a decollement near their contact fig 2a because these seismic profiles define the large scale geometry so clearly the objective of the leg 87 drilling programme was to investigate the structural fabric and physical properties of sediments as they enter the subduction zone using these data which we now report we hoped to deduce or constrain the stress conditions accompanying deformation