microseismicity and velocity structure in the hawkes bay region new zealand fine structure of the subducting pacific plate summary hawkes bay overlies part of the shallow zone of subduction of the pacific plate beneath the australian plate in the north island of new zealand a temporary dense network of 11 seismograph stations placed in the region recorded more than 350 micro earthquakes in two months providing a data set for velocity and structure inversion the hypocentres define a seismic band 11 12 km thick dipping 6 2 beneath hawke bay and increasing to more than 20 towards the west side of the region this seismic band is thought to delineate the crust of the subducting pacific plate s to p converted phases for deeper events are consistent with a fairly large impedance contrast across the plate interface at the top of the band p and s arrival times were inverted to give a simultaneous determination of hypocentres and velocity structure results indicated a thin layer of low 4 7 ksm 1 p wave velocity in the neighbourhood of the plate interface we interpret this layer as a thickness of subducted sediments beneath the micro earthquake station network this is consistent with seismic reflection work to the south of hawkes bay and with sediment volume calculations the stress distribution as given by focal mechanisms shows that the microseismicity represents primarily intraplate rather than interplate deformation with the main band of microseismicity indicating normal faulting throughout at least the top 10 km of the subducted pacific plate the observation of reverse faulting at depth beneath the main band is interpreted as being due to the flexure of the subducted pacific plate

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

oblique convergence and deformation along the kuril and japan trenches the hypothesis that present day deformation within the southern kuril forearc is driven by oblique subduction of the pacific plate is tested using 397 horizontal slip directions derived from shallow thrust earthquakes from the kuril and japan trenches for the period 1963 1991 a simple two plate model fits the 397 slip vectors significantly worse than a model that permits strike slip motion of the southern kuril forearc relative to the overlying plate weighted mean slip directions along the southern kuril trench are systematically rotated toward the direction orthogonal to the trench which implies that the net convergence is partitioned into less oblique subduction and trench parallel displacement of the southern kuril forearc the angular discrepancy between the observed slip direction and the direction predicted by the nuvel 1 pacific north america euler vector implies that the southern kuril forearc translates 6 11 mm yr 1 to the southwest relative to the overlying north american plate these results are consistent with geologically geodetically and seismologically observed convergence at the leading edge of the forearc sliver in southern hokkaido and with previously inferred extension at the trailing edge of the sliver which is located at the bussol strait at 46 n the percentage of the oblique component of the rigid plate convergence vector that is transferred into forearc strike slip motion is only 30 which may be evidence that the hokkaido honshu bend in the trench hinders southwest forearc motion or that the greater strength of oceanic forearcs impedes motion along forearc strike slip faults

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

changes to the environment of intramontane basins in the light of malacological research of calcareous tufa podhale basin carpathians southern poland the profiles of calcareous tufa occurring in the podhale basin carpathians southern poland were subject to a detailed malacological analysis studies were carried out in 15 sites five of them represent the period of the late glacial and holocene and the remainder is tufa deposited during historical times on the basis of malacological analysis 10 faunistic assemblages were defined each developed under specific environmental and climatic conditions the observed sequences of molluscan assemblages in vertical profiles were the basis for paleoenvironmental reconstruction determining the main phases of climatic changes and the human impact on the natural environment the faunal sequence observed indicates a cold climate and the predominance of open habitats in the late glacial with progressing warming and the increasing significance of shaded habitats in the early and middle holocene as well as anthropogenic pressure in modern times the results of malacological analysis were complemented with determination of the age of deposits by radiocarbon dating and the results of palynological analyses conducted in the sites of peatbogs occurring in podhale the reconstruction indicates similarities and differences in the paleogeographic development of mountain and upland areas of central europe

centroid moment tensor solutions for january march 1992 abstract centroid moment tensor solutions are presented for 160 earthquakes that occurred during the first quarter of 1992 the solutions are obtained using corrections for aspherical earth structure represented by a whole mantle shear velocity model sh8 u4l8 of dziewonski and woodward

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

seismicity and geometry of a 110 km long blind thrust fault 2 synthesis of the 1982 1985 california earthquake sequence during 1982 1985 three 5 4 mw 6 5 earthquakes migrated 65 km along the north half of a string of quaternary folds at the east front of the california coast ranges there is no surface fault associated with these earthquakes but the fold chain oriented parallel to the san andreas fault 30 km to the west takes up contraction normal to the san andreas fault identification of active blind thrust faults is hindered by the absence of a fault trace but may be revealed by the presence of growing folds at the surface and by earthquakes at depth because of excellent seismic geodetic and geologic data this earthquake sequence illuminates the geometry of the blind faults and enables us to probe their mechanics we relocate seismicity and compute focal mechanisms in a three dimensional velocity space and model the vertical deformation associated with the 1983 m 6 5 coalinga earthquake the aftershock zones abut at echelon offsets in the fold axes and the mainshocks display reverse slip perpendicular to the axes suggesting that the folds conceal a contiguous segmented thrust fault background seismicity concentrates at bends and breaks in the fold chain sites that may correspond to tears and ramps in the thrust fault at depth seismic reflection profiles reveal thrust and reverse faults dipping toward the san andreas fault at depths of 5 10 km with several kilometers of cumulative slip and high angle reverse faults in the anticlines with several hundred meters cumulative slip coseismic fold uplift accompanied the coalinga and kettleman hills north dome earthquakes suggesting that fold growth is episodic and coupled closely to repeated earthquakes on the underlying thrusts the north half of the fold chain has been the site of several 6 m 6 5 earthquakes since 1885 thus kettleman hills middle dome the next fold segment to the south has an elevated seismic potential the aftershock zones of the three main shocks are diffuse and occupy a region much larger than the site of seismic slip aftershocks occur 5 7 km from the fault regions where on the basis of a boundary element model the shear strain increase caused by the fault slip exceeds 20 ppm equivalent to about 0 7 mpa we argue that the broad aftershock zone is a product of high sustained off fault stress caused by repeated displacement on faults that do not cut the earth s surface the fault tip stresses lead to the formation of secondary faults which can become sites of aftershocks and postseismic creep

properties of aftershock sequences in southern california the temporal behavior of 39 aftershock sequences in southern california 1933 1988 was modeled by the modified omori relation minimum magnitudes for completeness of each sequence catalog were determined and the maximum likelihood estimates of the parameters k p and c with the standard errors on each were determined by the ogata algorithm the b value of each sequence was also calculated many of the active faults in the region both strike slip and thrust were sampled the p values were graded in terms of the size of the standard error relative to the p value itself most of the sequences were modeled well by the omori relation many of the sequences had p values close to the mean of the whole data set 1 11 0 25 but values significantly different from the mean as low as 0 7 and as high as 1 8 exist no correlation of p with either the b value of the sequence or the mainshock magnitude was found the results suggest a direct correlation of p values is with surface heat flow with high values in the salton trough high heat flow and one low value in the san bernardino mountains and on the edge of the ventura basin both low heat flow the large fraction of the sequences with p values near the mean are at locations where the heat flow is near the regional mean 74 mw m2 if the hypothesis that aftershock decay rate is controlled by temperature at depth is valid the effects of other factors such as heterogeneity of the fault zone properties are superimposed on the background rate determined by temperature surface heat flow is taken as an indicator of crustal temperature at hypocentral depths but the effects on heat flow of convective heat transport and variations in near surface thermal conductivity invalidate any simple association of local variations in heat flow with details of the subsurface temperature distribution the interpretation is that higher temperatures in the aftershock source volume caused shortened stress relaxation times in the fault zone materials leading to a faster decay rate higher p value in 1967 mogi proposed such a relation on the basis of japanese data seismic velocity distributions published for southern california generally support the hypothesis with low velocities higher temperatures corresponding to high p values and vice versa

upper mantle p wave velocities in the southeast north island new zealand abstract reversed pn velocities have been determined from regional earthquake sources for the wairarapa region in the southeast north island new zealand values were 8 64 0 13 km s in the easternmost part of the region and 8 90 0 13 km s in the central wairarapa which are consistent with subducted pacific plate lying at shallow depth under the region a value of 8 22 km s between wellington and the wairarapa region is interpreted as a possible difference in the velocity structure above the refractor between the two localities

spatial and temporal correlation between coda q 1 and seismicity and its physical mechanism from the analysis of local earthquake coda waves recorded by the benioff short period vertical seismograph at riverside california we found a surprisingly large and systematic temporal variation in coda q 1 in the southern california during the 55 year period from 1933 to 1987 comparison of the coda q 1 variation with the b value obtained for earthquakes with m 3 occurred in the area within 180 km from riverside corresponding to the area for which q 1 is estimated by the coda method led to an extraordinary result that coda q 1 correlates positively with b value with the correlation coefficient of 0 79 the significance of correlation cannot be rejected at the confidence level of 0 974 since the data used for measuring q 1 namely the coda part of seismograms and the data used for measuring b namely the seismicity catalog are totally independent the observed correlation supports strongly the physical reality of temporal change in both coda q 1 and b value the correlation between coda q 1 and b value has been reported for several other cases but the correlation was positive in some cases and negative in others our current preferred model for explaining the observed correlation is the following aseismic creep is going on from time to time over various sizes of cracks in the ductile part and in some cases in the brittle part of lithosphere at and near the plate boundary creep activities tend to increase the crack density and coda q 1 observed in the region if the creek occurs over distributed cracks with a certain predominant scale length spatial stress concentration may enhance seismicity for earthquakes with a similar size if this characteristic size is in the lower part of the magnitude range for which the b value is estimated the enhanced seismicity will increase b on the other hand if it is in the upper part of the magnitude range it will decrease b thus we can explain the observed strong correlation between coda q 1 and b value for both positive and negative cases the existence of such a characteristic crack size in a seismic region has been supported by various observations on earthquake scaling

seismic gap hypothesis ten years after the seismic gap hypothesis states that earthquake hazard increases with time since the last large earthquake on certain faults or plate boundaries one of the earliest and clearest applications of the seismic gap theory to earthquake forecasting was by mccann et al 1979 who postulated zones of high medium and low seismic potential around the pacific rim in the 10 years since there have been over 40 large m 7 0 earthquakes enough to test statistically the earlier forecast we also analyze another forecast of long term earthquake risk that by kelleher et al 1973 the hypothesis of increased earthquake potential after a long quiet period can be rejected with a large confidence the data suggest that contrary to these forecasts places of recent earthquake activity have larger than usual seismic hazard whereas the segments of the circum pacific belt with no large earthquakes in recent decades have remained relatively quiet the clustering of earthquake times does not contradict the plate tectonic model which constrains only the long term average slip rate not the regularity of earthquakes

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

depth of seismic coupling along subduction zones underthrusting at subduction zones can cause large earthquakes at shallow depths but it is always accommodated by aseismic deformation below a certain depth the maximum depth of the seismically coupled zone or seismogenic zone is a transition from unstable to stable sliding along the plate interface we have determined the depth of this stability transition for the circurn pacific subduction zones of honshu kuriles kamchatka aleutians alaska mexico and chile these subduction zones have experienced great interplate earthquakes and the aftershock regions are well located depth estimates of interplate events that are located at the downdip edge of the aftershock regions are used to determine the maximum depth of seismic coupling for an average p wave velocity of 6 7 km s 1 above the plate interface we find that for most subduction zones the stability transition occurs at 40 5 km depth there are however several exceptions at the hokkaido trench junction where the japan trench and the kurile trench intersect seismic coupling is deep and extends down to 52 55 km deep coupling was also found in the coquimbo region in central chile the mexico subduction zone has shallow coupling the transition occurs at 20 30 km depth previous studies of micro earthquakes in honshu hokkaido the aleutians and alaska show that earthquakes within the upper plate extend no deeper than the downdip edge of the coupled zone that we find given our measurements of seismic coupling depth we then explore the mechanism that may determine coupling depth the concept of critical temperature has been used to explain the depth of seismic coupling in other tectonic environments thus we first test whether a critical temperature can explain our results temperatures at the plate interface are dependent on many variables but two that are poorly determined are shear stress and radiogenic heat generation shear stress has been constrained by inversion of heat flow data assuming a crustal radiogenic heat production rate of 3 1 exp z 8 5 wm 3 and a constant coefficient of friction we find two critical temperatures of about 400 c and 550 c the lower critical temperature may be characteristic of regions with a relatively thick continental crust and the higher temperature of regions with a relatively thin continental crust on the other hand one single critical temperature of about 250 c can explain the coupling depths if shear stresses are constant with depth

precursory seismic quiescence past present and future precursory seismic quiescence has played a major role in most of the succesful earthquake predictions made to date in addition to these successes the number of detailed post mainshock documentations of precursory quiescence is steadily growing these facts suggest that precursory quiescence will play an important role in earthquake prediction programs of the future for this reason it is important to critically evaluate the present state of knowledge concerning this phenomenon the history of observations of precursory seismic quiescence includes work on seismic gaps and seismic preconditions as well as actual studies of temporal quiescence these papers demonstrated the importance of quantitative evaluation of seismicity rates and the benefits of systematic analysis during the early 1980 s the impact of man made effects on seismicity rates was demonstrated for the first time despite progress in catalog understanding the identification and correction of man made seismicity changes remains as the major barrier to earthquake prediction using these data effects of man made changes are apparent in many past studies of seismicity patterns making the results difficult to evaluate recent experience with real time anomalies has demonstrated the necessity of determining the false alarm rates associated with quiescence precursors determination of false alarm rates depends on quantitative definitions of anomalies and statistical evaluations of their significance a number of successful predictions which have been made on the basis of seismic quiescence provide important lessons for present and future work there are many presently unanswered questions regarding seismic quiescence which must be answered before we can determine the reliability of this phenomena as a precursor

evidence for precursory changes in the frequency magnitude b value summary a clear case has been established for a significant increase in the slope of the frequency magnitude relation for periods of up to six years before a number of earthquakes in new zealand and california in a region of about 200 km diameter in central new zealand the b value increased from near unity to more than 2 0 at the end of 1968 and remained abnormally high until the beginning of 1975 an earthquake of magnitude 6 0 followed in 1977 january within 50 km of the location which had registered the highest b value the anomaly was detected in contour maps of b values windowed in space and time then resolved more precisely in time using the cusum technique it was statistically significant at a level higher than 99 99 per cent compared with the average rate of moment release in this region there was a deficit of 1 3 x 10 nm during the anomalous period which is equivalent to an earthquake of magnitude 6 08 this phenomenon of a high b value has been found before all shallow new zealand earthquakes of magnitude 6 0 or more since 1965 wherever the seismicity rate and detection threshold were adequate to resolve it further the only statistically significant occurrences of high b value were those which were followed by large earthquakes two of magnitude 6 8 and 6 0 and a swarm of equivalent magnitude 5 3 the phenomenon has also been found before the san fernando and coalinga earthquakes in california both of which were thrust events but it is not apparent before two large californian earthquakes with strike slip mechanisms observations of high b value are consistent with reported instances of quiescence in background seismicity in fact what has been interpreted as quiescence may probably be more accurately described as a localized high b value the duration of the anomaly is not a simple function of magnitude different strain rates in the regions examined may account for this the available data suggest that careful monitoring of the b value may provide a tool for medium term earthquake prediction although perhaps not in all tectonic environments

an improved method for determining the regional stress tensor using earthquake focal mechanism data application to the san fernando earthquake sequence the orientations of fault planes and slip directions indicated by a population of earthquake focal mechanisms can be used to determine best fit regional principal stress directions and r 2 1 3 1 a measure of relative stress magnitudes under the assumption of uniform stress in the source region this analysis allows for the possibility that failure occurs on preexisting zones of weakness of any orientation in the inversion we perform a grid search of stress models to find the one which requires the smallest total rotation of all the fault planes that is needed to match the observed and predicted slip directions the method allows for errors in orientations of both the fault planes and slip directions we have an objective means for identifying the more likely of the two possible fault planes from each focal mechanism relative to a given stress model thus we do not face the problem of ambiguity of nodal planes which plagues other analyses of this kind by using a grid search of stress models rather than a linearization scheme we are able to perform a realistic error analysis and thus establish confidence limits for the preferred regional stresses the method can be used to investigate possible stress inhomogeneities during earthquake sequences by analyzing subsets of the data population the technique has been applied to 76 events from the san fernando earthquake sequence for which we have found best fit stresses plunge and azimuth 1 7 187 2 27 281 3 62 84 and r 0 65 the average misfit between this stress model and all the data is about 8 and all but eight of the aftershocks have misfits of less than 20 these values are considerably less than the uncertainty of the focal mechanism determinations therefore significant stress inhomogeneities are not required by the data our analysis does not support the suggestion of a change in stresses during the aftershock sequence as proposed by others on the basis of an apparent change in focal mechanisms

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

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

phosphorus burial and diagenesis in the central bering sea bowers ridge iodp site u1341 perspectives on the marine p cycle abstract to reconstruct the cycling of reactive phosphorus p in the bering sea a p speciation record covering the last 4 ma was generated from sediments recovered during integrated ocean drilling program iodp expedition 323 at site u1341 bowers ridge a chemical extraction procedure distinguishing between different operationally defined p fractions provides new insight into reactive p input burial and diagenetic transformations reactive p mass accumulation rates mars are 20 110 mol cm2 ka which is comparable to other open ocean locations but orders of magnitude lower than most upwelling settings we find that authigenic carbonate fluorapatite cfa and opal bound p are the dominant p fractions at site u1341 an overall increasing contribution of cfa to total p with sediment depth is consistent with a gradual sink switching from more labile p fractions fish remains fe oxides organic matter to stable authigenic cfa however the positive correlation of cfa with al content implies that a significant portion of the supposedly reactive cfa is non reactive detrital contamination by eolian and or riverine cfa in contrast to cfa opal bound p has rarely been studied in marine sediments we find for the first time that opal bound p directly correlates with excess silica contents this p fraction was apparently available to biosiliceous phytoplankton at the time of sediment deposition and is a long term sink for reactive p in the ocean despite the likelihood for diagenetic re mobilisation of this p at depth indicated by increasing ratios of excess silica to opal bound p average reactive p mars at site u1341 increase by 25 if opal bound p is accounted for but decrease by 25 if 50 of the extracted cfa fraction based on the lowest cfa value at site u1341 is assumed to be detrital combining our results with literature data we present a qualitative perspective of terrestrial cfa and opal bound p deposition in the modern ocean riverine cfa input has mostly been reported from continental shelves and margins draining p rich lithologies while eolian cfa input is found across wide ocean regions underlying the northern hemispheric dust belt opal bound p burial is important in the southern ocean north pacific and likely in upwelling areas shifts in detrital cfa and opal bound p deposition across ocean basins likely occurred over time responding to changing weathering patterns sea level and biogenic opal deposition