hypocentre determination offshore of eastern taiwan using the maximum intersection method summary the maximum intersection maxi method which derives from the master station method msm determines within a 3 d velocity model the absolute hypocentral location based on observed arrival times first the spatial node that better satisfies the arrival time differences computed at all station pairs plus or minus an error tolerance value in seconds is defined as the preliminary hypocentral solution pred second because pred depends neither on the estimate of origin time nor on the residual root mean square rms residual outliers are objectively detected and cleaned out from the original data set without any iterative process or weighting third a statistical minimization residual rms is conducted in a small domain around the pred node which results in a unique final solution the maxi method is applied to the determination of earthquake hypocentres with the proper station correction terms in the southernmost extremity of the ryukyu subduction zone where several dense seismic clusters occur near the seismogenic plate interface the location of earthquakes recorded at both the taiwanese and japanese networks is obtained for about a thousand events between 1992 and 1997 the process uses a detailed 3 d velocity model based on multiple geophysical data sources obtained in the junction area between subduction and collision east of taiwan the earthquake clustering and the significant drop in residual statistics 1 20 0 80 and 0 35 s for taiwanese catalogue msm and maxim solutions respectively indicate the accuracy of the method which can be used to routinely determine absolute hypocentre location based on observed arrival times

is taiwan the result of arc continent or arc arc collision abstract conventionally it has been accepted that the formation of taiwan results from the collision of the luzon arc with the eurasian continental margin we suggest that taiwan results from the collision of the luzon arc with the former ryukyu subduction zone before the collision the latter extended a few hundred kilometres southwest of its present day termination in the early to middle miocene the subduction became inactive in its southwest portion and the manila trench migrated northeastward giving rise to the formation of the luzon arc the collision started in the late miocene by closing the oceanic domain located between the luzon arc and the former ryukyu arc because the luzon arc moved northwestward with respect to the former ne sw trending ryukyu arc the oblique collision first resulted in the indentation on the south ryukyu arc west of 123 5 e and the narrowing of the corresponding portion of the okinawa trough prior to the uplift of taiwan today the hsuehshan and backbone ranges represent the uplifted portion of the backarc basin and the longitudinal valley corresponds to the suture zone between the former ryukyu and the luzon arcs in this model we propose to link although they are not of the same age two truncated backarc basins the active okinawa trough to the northeast and the former active tainan backarc basin to the southwest the former backarc basin extended from southwest of taiwan to southern japan and could have been initiated near the paleo location of taiwan as early as the eocene oligocene the result of the collision is the truncation of the ancient backarc basin associated with the shortening deformation and uplift of the ancient outer arc and related backarc region which provides the fundamental mechanism of the taiwan mountain building this model is in agreement with a counterclockwise rotation of the portion of taiwan located west of the coastal range belonging to the luzon arc and a clockwise rotation of the south ryukyu arc earthquake focal mechanisms geodetic observations and the distribution of the petrologic and structural features are also consistent with the proposed arc arc collision model

crustal deformation at the southernmost part of the ryukyu subduction east taiwan as revealed by new marine seismic experiments abstract the southernmost part of the ryukyu subduction where the philippine sea plate is subducting under the eurasian plate is known to be a very seismically active region of transition from a north dipping subduction along the ryukyu subduction to an se nw collision along the taiwanese orogenic wedge in this paper we will focus on the ryukyu forearc area close to taiwan where the deformation is paroxysmal in order to decipher the nature of the seismic deformation in this region a three month passive experiment combining 22 ocean bottom seismometers and 51 onland stations has been led starting from an a priori heterogeneous model we have obtained 801 well located earthquake hypocenters a precise p wave tomography model and 14 focal mechanisms the seismicity along the ryukyu forearc is mainly located not only in the vicinity of the interplate seismogenic zone isz but also within both the subducting psp and the overriding plate seismicity within the upper plate is essentially localized east of nanao basin where nw se extension occurs and northwest of the hoping basin where strike slip dominates as revealed by both the p wave velocity structure and the newly derived seismicity we argue that a sub vertical step offsetting the subducting psp around 10 km may support the presence of a trench parallel tear the psp also undergoes extension in its upper part that is probably caused by buckling and slab pull the p wave velocity structure reveals three other major features 1 a continuity between the central range and the ryukyu arc with a shallower moho 30 km depth between 122 3 n and 122 5 n along the ryukyu arc 2 high p wave velocities along the eastern side of the central range and 3 two bodies with similar high crustal velocities 6 5 7 0 km s at 12 18 km depths embedded within the ryukyu arc basement just north of hoping basin and north of the nanao basin

hypocentral trend surface analysis probing the geometry of benioff zones a hypocentral trend surface is a continuous function of latitude and longtitude fitted by least squares to a set of hypocenters so that it predicts depth to the middle of a benioff zone in this paper we take a relatively simple approach to hypocentral trend surface analysis the hypocentral trend surface is constructed from a spherical surface harmonic expansion whose coefficients are selected so as to minimize the standard vertical deviation between hypocenters and trend surface the vertical deviation of a hypocenter from the hypocentral trend surface is called its residual consideration of a hypocentral trend surface cannot be divorced from consideration of the associated residuals and their spatial distribution any analysis of benioff zone geometry trades off presumed thickness of the zone with its presumed shape this trade off can be investigated by examining suites of candidate hypocentral trend surfaces together with their associated residual distributions generated by varying the number of degrees of freedom available to the trend surface hypocentral trend surfaces are generated for three high quality sets of hypocenters obtained by local seismic networks in honshu japan cook inlet alaska and south peru hypocenters beneath honshu generate a bimodal distribution of residuals about their trend surface the vertical separation of the upper and lower sheets of this double benioff zone averages about 33 km no double benioff zone configuration is discovered below cook inlet or south peru regional trend surfaces are established for intermediate depth teleseismic data from south and middle america in south america 0 40 s teleseismic data suggest that along strike transitions between flat and moderately steeply dipping sections of the benioff zone are achieved by flexure of a coherent slab rather than by fragmentation of the slab into tear bound flaps or separate tongues of lithosphere hypocentral trend surface analysis of teleseismic data from south peru suggests a geometry similar to that indicated by hypocentral trend surface analysis of the local data of hasegawa and sacks an eastward shallowing benioff zone beneath southwest mexico is closely associated with the geometry of the trans mexican volcanic arc the benioff zone beneath central america is more steeply dipping the flat slab steep slab transition near 96 w may relate to the difference in age of oceanic lithosphere on either side of the tehuantepec ridge while the strong flexure of the subducted nazca plate beneath ecuador is closely associated with the carnegie ridge the nazca ridge extropolates slightly north of the flexure in subducted lithosphere beneath south peru the flexure appears instead more closely associated with the pronounced inflection in the map view shape of the leading edge of the overriding south american plate

transition from interplate slip to double seismic zone along the kuril kamchatka arc during subduction the downgoing slab interacts with a strong overriding plate at shallow depths and with the surrounding asthenosphere at intermediate depths thus the transition from earthquakes representing interplate slip to those reflecting intraplate deformation within subducting lithosphere is important for understanding how the mechanics of subduction varies with depth to characterize this transition we studied 44 large to moderate sized earthquakes that occurred between 1964 and 1990 in two regions along the central portion of the kuril kamchatka arc precise earthquake source parameters are determined by analyzing body waveforms recorded at teleseismic distances at depths of 30 50 km our results delineate a transitional thrust zone characterized by earthquakes showing thrust faulting on moderately dipping nodal planes the northwesterly dipping nodal plane dips approximately 40 some 10 15 steeper than the dip of low angle earthquake faulting along the interplate thrust zone at depths shallower than 30 km for earthquakes along the transitional thrust zone secondary features in observed waveforms suggest that either a thick layer of sediments is present in the forearc region or the sources are embedded in materials with crustal seismic wave speeds in the former case the transitional thrust zone is interpreted as the deep seated portion of the plate interface whose dip increases with depth the alternative implies that the transitional thrust zone occurrs below the plate interface within subducted crust of the pacific plate and that the occurrence of earthquakes along the transitional zone may indicate a locked plate interface between depths of 30 and 50 km the inference of a locked plate interface is consistent with the occurrence of compressional earthquakes beneath the outer rise region however a locked interface does not seem to affect the state of strain in the double seismic zone the upper compressional layer of the double seismic zone abuts the transitional thrust zone at a depth of approximately 50 km at this depth the two layers of seismicity in the double seismic zone are separated by as much as 30 40 km too large a distance to be accommodated by models of seismogenesis invoking either the gabbro eclogite transition or dehydration of subducted oceanic crust

kinematics of convergence deformation and stress distribution in the taiwan collision area 2 d finite element numerical modelling abstract using a 2 d plane stress finite element model with elastic and elasto plastic rheologies appropriate for deformation within the brittle upper crust we analyse the relationship between kinematics of convergence deformation and stress distribution in the present taiwan collision occurring between the ryukyu and luzon subduction zones the distribution of stress trends calculated in our models is compared with a synthetic map of actual stress trajectories based on the most recent data available in the collision zone these data combine present day sources from borehole breakouts and earthquake focal mechanism with the reconstruction of quaternary palaeostress from fault slip data analyses resulting in a complete map of compressional stress trajectories which is used to constrain our models we show that the distribution of stress trajectories in the active taiwan collision is principally controlled by 1 the geometric configuration of the boundary between eurasia and the philippine sea plate 2 the shape and rheological properties of major structural units 3 the direction of convergence of the philippine sea plate relative to eurasia and 4 the influence of the opening of the okinawa trough the study of a two dimensional elastic and elasto plastic finite element modelling of the subduction collision in and around taiwan allows us to estimate the influences of these different parameters in the stress pattern taking into account both the simplifying assumptions of the numerical modelling and the angular uncertainties of field determinations the fit between the calculated stress pattern of the finite element model and that determined based on the geometrical synthesis of field analyses is rather good in general indicating that our model is valid to first approximation misfits remain minor and can be explained by data uncertainties and simplifying modelling assumptions for instance the shape of the corner of the collision zone is critical but is not accurately known also limited decoupling in the longitudinal valley collision zone was not considered in our models although it certainly plays a role some interesting features of our model are 1 the greater influence of the shape of the collision zone in comparison with that of the direction of convergence 2 the requirement for a trench retreat related to suction force in the ryukyu arc and 3 the crucial role of the interaction between okinawa trough opening and collision at the sharp northwestern corner of the philippine sea plate including its influence on the geological evolution of northeastern taiwan

evidence for a role of the downgoing slab in earthquake slip partitioning at oblique subduction zones a new model incorporating shear deformation within a subducting slab is proposed to explain slip partitioning for oblique plate motion at subduction zones on the basis of investigation of 450 interplate earthquakes at 24 subduction zone segments we find that the degree of slip partitioning is largely correlated with the calculated slab pull force such correlation suggests that other than the upper plate deformation the slab pull force plays an important role in controlling oblique subduction our model proposes that the force balance between the slab pull force the interplate coupling resistance and the viscous mantle drag the latter two are passive forces to the former one produces a lateral shear within the slab which causes the slab to deform and change its motion direction gradually toward trench normal as it subducts the amount of direction change which would be observed as the slip partitioning during earthquakes therefore is closely related to the major plate driving force at the subduction zone which is the slab pull force in our model

tectonic analysis and the evolution of a curved collision belt the hsuehshan range northern taiwan abstract detailed studies of fracture patterns and paleostress reconstructions based on analyses of fault slip data sets allow identification of four main patterns of compressional paleostress orientations in the cenozoic rocks of northern taiwan these patterns correspond to five successive tectonic events dominated by 1 n s 2 wsw ene 3 wnw ese 4 nw se and 5 wnw ese directions of compression the first two events are characterized by relatively minor structures in the foreland fold and thrust belt the first event probably reflects major left lateral shear along nne ssw trends which occurred in the adjacent eastern metamorphic belt during the middle late miocene these first two events do not affect the plio pleistocene in the localities studied in contrast the last three events created the major compressional structures of the foreland belt during the plio pleistocene collision in taiwan preliminary studies suggest that the magnetic fabric of the rocks has been significantly influenced by compressional deformation along wnw ese and nw se trends especially in the slate belt thus confirming the importance of these last three events on the basis of paleomagnetic work in progress the northernmost segment of the taiwan belt has been affected by a clockwise rotation of about 20 taking into account the distribution of successive compressional paleostress patterns we conclude that of the 55 angle of the sharply curved belt of northern taiwan about 20 probably correspond to actual clockwise rotation during the plio pleistocene whereas the remaining 35 reflect the initial shape of the curved plio quaternary belt the paleostress history of northern taiwan is examined in terms of geodynamic evolution and kinematic reconstructions the plio quaternary development of the collision between the luzon arc and eurasia followed a late miocene predominantly transform motion from south to north of the philippine sea plate relative to eurasia the succession of plio pleistocene paleostress patterns in northern taiwan is found to be consistent with the longitudinal shift from northeast to southwest and the transverse migration from the inner to the outer zones of compressional phenomena related to the evolution of oblique collision some major changes in compressional trends are also accounted for by variations in the kinematics of plate convergence the most recent charge comprises a counterclockwise change in the direction of philippine sea plate eurasia convergence

shear concentration in a collision zone kinematics of the chihshang fault as revealed by outcrop scale quantification of active faulting longitudinal valley eastern taiwan abstract repeated measurements of active deformation were carried out at three sites along the active chihshang fault a segment of the longitudinal valley fault zone of eastern taiwan the present day plate boundary between the philippine sea plate and eurasia reliable annual records of displacement along an active fault were obtained based on detailed surveys of faulted concrete structures along the active chihshang fault striking n18 e we determined average motion vectors trending n37 w with an average shortening of 2 2 cm yr thus the transverse component of motion related to westward thrusting is 1 8 cm yr whereas the left lateral strike slip component of motion is 1 3 cm yr the fault dips 39 45 to the east so that the vertical displacement is 1 5 3 cm yr and the actual oblique offset of the fault increases at a rate of 2 7 3 7 cm yr this is in good agreement with the results of regional geodetic and tectonic analyses in taiwan and consistent with the n54 w trend of convergence between the northernmost luzon arc and south china revealed by gps studies our study provides an example of extreme shear concentration in an oblique collision zone at chihshang the whole horizontal shortening of the longitudinal valley fault 2 2 cm yr on average occurs across a single narrow fault zone so that the whole reverse slip about 2 7 3 7 cm yr depending on fault dip was entirely recorded by walls 20 200 m long where faults are tightly localized this active faulting accounts for more than one fourth 27 of the total shortening between the luzon arc and south china recorded through gps analyses further surveys should indicate whether the decreasing shortening velocity across the fault is significant revealing increasing earthquake risk due to stress accumulation or not revealing continuing fault creep and weak behaviour of the chihshang fault

large earthquakes and crustal deformation near taiwan source parameters of 15 large mb 6 earthquakes which occurred near taiwan during the last 25 years are determined by synthetic seismogram analysis of teleseismic p and s waveforms source depths range from 13 to 88 km three earthquakes of intermediate depth are located within the southernmost limit of the ryukyu arc and show thrust type mechanisms with t axes subparallel to a northward dipping wadati benioff zone shallow earthquakes 30 km are distributed over a broad region along and off the east coast of taiwan and show both high angle reverse and strike slip mechanisms the p axes of the shallow events located off the east coast of northern taiwan are oriented near horizontal trend south to southwest and may reflect relative motion between the okinawa platelet and taiwan or complex deformation at the southwest limit of the ryukyu trench in contrast p axes of shallow earthquakes both onshore and off the east coast of central and southern taiwan consistently trend northwest and are interpreted to reflect relative convergence between the philippine sea and eurasian plates summation of the seismic moment tensors of eight large shallow earthquakes that occurred along and off the central and southern coast of taiwan shows that relative plate motion between the eurasian and philippine sea plates when averaged over the last 25 years has been accommodated by crustal contraction at a rate of 2 6 5 4 cm yr along an azimuth of about 290 o

a synthesis of the geologic evolution of taiwan abstract the island arc of taiwan is composed of cenozoic geosynclinal sediments more than 10 000 m thick lying on a pre tertiary metamorphic basement pleistocene to miocene andesitic islands surround the main island and are related mostly to arc magmatism the penghu island group in the taiwan strait is covered with pleistocene flood basalt neogene shallow marine clastic sediments are exposed mainly in the western foothills with pleistocene andesitic extrusives at the northern tip and the northeastern offshore islands a thick sequence of paleogene to miocene argillitic to slaty metaclastic rocks underlies the western central range and forms the immediate sedimentary cover on the pre tertiary metamorphic complex to the east which represents an older mesozoic arc trench system the coastal range in eastern taiwan is a neogene andesitic magmatic arc including also a large variety of volcaniclastic and turbiditic sediments cenozoic taiwan is the site of arc continent collision where the luzon arc on the philippine sea plate overrides the chinese continental margin on the eurasian plate east and northeast of taiwan the polarity of subduction changes whereby the oceanic philippine sea plate is subducting beneath the ryukyu arc system on the eurasian plate continent arc collision in taiwan island is anomalous and may occur in a broad belt of deformation rather than along a well defined plate boundary or subduction zone

active tectonics of the pamirs and karakorum the source parameters of moderate to large sized earthquakes and the spatial distribution of earthquake hypocenters are used to investigate the active tectonics and three dimensional configuration of the subducted lithosphere in the pamir karakorum intracontinental convergence zone the continental collision has resulted in intense deformation in the upper crust the deformation is distributed over a broad area and is absorbed by thrust and strike slip faulting along the northern margin of the pamirs where the deformation due to continental collision is most intense the northern pamir thrust and several major fault zones mark the present boundary most of the earthquakes that occurred in these fault zones are shallow crustal events with focal mechanism solutions correlated with the tectonic features of the region along the pamir front thrust faulting dominates while on the western and eastern edges of the pamirs the style of deformation is characterized by oblique thrusting with a component of strike slip motion thrust type events beneath the tadjik depression indicate that both the sedimentary rocks and the basement are involved in shortening right lateral strike slip motion is observed on the eastern edge of the pamirs where they border the tarim basin and in the talas fergana kun lun fault zone a southward dipping seismic zone beneath the pamirs and karakorum indicates that the asian lithosphere has been subducted along the pamir front to a depth of at least 200 km this interpretation is consistent with published geological and geophysical data focal mechanism solutions of some intermediate depth events beneath the pamirs show strike slip faulting with approximately n s horizontally oriented p axes indicating that nearly horizontal compression at the intermediate depth is the predominant mode of deformation a 90 km deep event beneath the karakorum is interpreted as occurring in the subducted asian lithosphere the fact that the p axis of this event is oriented parallel to the descending plate suggests downdip compression within the asian lithosphere the overall tectonics of the region is interpreted as a consequence of the underthrusting asian lithosphere being impinged upon by the shallow northward underthrusting indian lithosphere

microseismicity of back arc rifting in the middle okinawa trough a one month microearthquake survey was conducted at iheya knolls in the middle okinawa trough in 1992 to reveal microearthquake characteristics of back arc rifting we deployed five ocean bottom seismometers and observed over 1700 events during 30 days hypocenters were determined for 48 earthquakes with good resolution from this survey we conclude that the survey area with 100km 100km has high seismicity many earthquakes occur in the western part of iheya knolls and few occur in the central part the earthquakes are distributed in an en echelon shape and throughout the entire depth of the crust there is no evidence for the existence of transform faults suggesting that the middle okinawa trough is in a stage of back arc rifting

the configuration of asia prior to the collision of india cretaceous paleomagnetic constraints paleomagnetic data from central asia show that 1700 610 km of shortening of southern asia since cretaceous time have been absorbed by distributed deformation between southern tibet and the siberia craton this result is based on a compilation of cretaceous poles from the junggar tarim tibet indochina south china north china and mongolia blocks complementing the recent compilation of enkin et al 1992a we propose a paleogeographic reconstruction of asia in the cretaceous in which the position of siberia is derived from the synthetic apparent polar wander path of besse and courtillot 1991 the resulting map which likely represents asia as it remained throughout the cretaceous until the collision with india began features an unbent tibet with an east west trending andean margin at tropical latitudes and a rather continuous belt of continental red bed basins extending from sichuan to tarim through tibet the map allows one to estimate continental shortening and rotations between the blocks which are attributed to the collision despite large uncertainties these have amounts and senses which are in all cases compatible with some recent kinematic models such as that of avouac 1991 appendix is available with entire paper on microfiche order from american geophysical union 2000 florida avenue n w washington dc 20009 document b93 007 2 50 payment must accompany order

a thermomechanical wedge model of taiwan constrained by fission track thermochronometry abstract we present 31 new apatite fission track aft ages for the island of taiwan that when combined with existing aft and zircon fission track zft data provide regional spatial coverage of the island with respect to low temperature thermochronometry the overall pattern of zft and aft ages in taiwan exhibits unreset ages in the southern and western portions of the island and reset ages predominantly in the central range and eastern taiwan this pattern supports interpretations of the orogen kinematics as reflecting a crustal scale wedge with a southward propagating collision zone in this model new material is accreted to the wedge from the west and is transferred to the east with the greatest exhumation occurring along the eastern margin as recorded in the reset ages in the east and unreset ages in the west the southward propagating collision is consistent with reset ages in the north where erosional exhumation has been ongoing for longer and unreset ages in the south where the younger collision implies less time for erosional exhumation despite the variation in the age of the collision along the strike of the island the widths of the aft and zft reset age zones remain nearly constant between 23 00 n to 24 00 n and 23 20 n to 24 00 n respectively suggesting that the orogen is in an exhumational steady state over these regions with respect to the aft and zft thermochronometers we use the fission track data in conjunction with observations of crustal structure crustal fabric and heat flow measurements to constrain a time dependent two dimensional thermomechanical model of orogen evolution by accounting for the heat transfer tectonic and erosion processes needed to predict aft and zft ages we are able to investigate the relationship between the measured ages and the tectonic characteristics of the orogen with our model we conclude that 1 roughly half of the material accretion in taiwan occurs through underplating over an approximately 40 km wide region 2 current average erosion rates are 3 3 mm yr in the eastern central range and 2 3 mm yr over the whole island 3 the collision has been propagating southward at a rate between 20 and 51 km ma over the past 2 3 ma and 4 central taiwan is in a topographic thermal and exhumational steady state

a catalog of taiwan earthquakes 1900 2006 with homogenized mw magnitudes taiwan has a relatively complete catalog of earthquakes since the first seismograph was installed in 1897 however due to changes in seismographic characteristics network coverage and observational practice the definition and procedure for magnitude determination were different during different time periods recognizing a complete catalog of earthquakes with consistent magnitudes is essential for delineating seismicity patterns and assessing seismic hazards for taiwan efforts have been made to convert the original magnitudes of earthquakes in taiwan based on various magnitude scales to a common m l or m s magnitude scale unfortunately the m l or m s magnitude scales chosen for previous studies all are subject to a fundamental limitation of saturation toward large earthquakes besides these studies are nearly two decades old in order to avoid this limitation and to follow the current trend we have chosen in this study to convert original magnitudes of various scales to a common m w scale in this study we used two independent methods for magnitude conversion in the first method we converted the original magnitudes to m w through empirical relations between these magnitudes and m w this magnitude is called the old m w in the second method we used the best fitting a and b values to convert the original magnitudes to m w this magnitude is called the new m w the converted m w magnitudes from both methods have resulted in significant improvements over the original magnitudes and are in good agreement with each other nevertheless by examining the log10 n versus m w plots we found better linearity and tighter overlap among different time periods for the new m w than for the old m w thus we chose the new m w as the unified magnitude for the catalog of taiwan earthquakes finally a list of 899 earthquakes from 1900 to 2006 with m w 5 5 is presented graphic 1 see the earthquake lists available in the electronic supplement to this article the list is considered complete for the whole period from 1900 to 2006 1 embed inline graphic 1 gif

a tear fault boundary between the taiwan orogen and the ryukyu subduction zone more than twenty seismic lines were analysed across the intersection of the taiwan orogen and the ryukyu arc forearc system the main results of this study concern the recognition of two superposed sedimentary basins that were previously unknown we have named the recent one the hoping basin and the older tectonized one the suao basin the suao basin containing sediments more than 3 km thick records subsidence of about 3 km on its southwestern portion after being shortened in the direction of plate convergence on the basis of seismicity gravity present day plate kinematics as well as seismic reflection data we relate the dramatic subsidence of the suao basin to the initiation of a wnw trending tear fault within the philippine sea plate this tear fault decouples the subducting philippine sea plate to the northeast and the apparent overthrusting philippine sea plate in the coastal range to the southwest back arc rifting may thus occur in the westernmost section of the okinawa trough above the sinking philippine sea slab while lithospheric shortening prevails south of hualien as a consequence the triangular area located southwest of the tear fault between the taiwan coastline and the toe of the yaeyama ridge no longer belongs to the ryukyu forearc it is now passively transported with the non subducting portion of the philippine sea plate we propose an evolutionary model for the relative chronology of tectono sedimentary events in this region over the last 5 m y

crustal scale seismic profiles across taiwan and the western philippine sea abstract we have used combined onshore and offshore wide angle seismic data sets to model the velocity structure of the taiwan arc continent collision along three cross island transects although taiwan is well known as a collisional orogen relatively few data have been collected that reveal the deeper structure resulting from this lithospheric scale process our southern transect crosses the hengchun peninsula of southernmost taiwan and demonstrates characteristics of incipient collision here 11 km thick transitional crust of the eurasian plate eup subducts beneath a large rapidly growing accretionary prism this prism also overrides the n luzon forearc to the east as it grows just west of the arc axis there is an abrupt discontinuity in the forearc velocity structure because this break is accompanied by intense seismicity we interpret that the forearc block is being detached from the n luzon arc and philippine sea plate psp at this point our middle transect illustrates the structure of the developing collision steep and overturned velocity contours indicate probable large scale thrust boundaries across the orogen the leading edge of the coherent psp appears to extend to beneath the east coast of taiwan deformation of the psp is largely limited to the remnant n luzon arc with no evidence of crustal thickening to the east in the huatung basin our northern transect illustrates slab continent collision the continuing collision of the psp and eup as the psp subducts the collisional contact is below 20 km depths along this transect ne of hualien this transect shows elements of the transition from arc continent collision to ryukyu arc subduction both of our models across the central range suggest that the paleozoic to mesozoic basement rocks there may have been emplaced as thick coherent thrust sheets this suggests a process of partial continental subduction followed by intra crustal detachment and buoyancy aided exhumation although our models provide previously unknown structural information about the taiwan orogen our data do not define the deepest orogenic structure nor the structure of western taiwan additional seismic active and passive geologic and geodynamic modeling work must be done to fully define the structure the active deformation zones and the key geodynamic process of the taiwan arc continent collision

elucidating the geometry of the active shanchiao fault in the taipei metropolis northern taiwan and the reactivation relationship with preexisting orogen structures the shanchiao fault is an active normal fault with documented paleoearthquakes in the taipei metropolis taiwan while posing direct seismic threat on the multimillion population its crustal scale fault plane configuration has not been constrained this study presents the first attempt to resolve the fault plane dip changes of the shanchiao fault within the upper crust by forward modeling late quaternary deformation tectonic subsidence over the last 23 ka is estimated from vertical displacements of a rapidly formed alluvial fan horizon deformed into a dramatic rollover monocline a 2 d profile across the shanchiao fault is chosen for elastic half space dislocation modeling and the results suggest that the fault is listric in the shallow crust with an abrupt change from subvertical ramp 85 75 to near horizontal flat 10 15 at 3 4 km depth consistent with an origin from the inversion of an orogen related thrust detachment given the presence of rift related fabrics in the underthrust chinese continental margin basement beneath the taiwanese orogenic wedge listric ramp flat ramp models with a second deeper bend to 60 dip are also tested reasonable fits with the geological observations are produced when the lower ramp is located at greater than 8 km depth which correlates with the hypocentral location of a moderate earthquake in 2004 joint reactivation of preexisting thrust and rift faults by the shanchiao fault is therefore plausible with implications for seismic hazard in the taipei area

the double seismic zone in kuril kamchatka the tale of two overlapping single zones the classic double seismic zone in the central portion of the kuril kamchatka arc seems to turn into a regular wadati benioff zone with downdip compression and extension in the northeastern and southwestern ends of the arc respectively to delineate this lateral changeover we studied 33 large to moderate sized mb 5 1 earthquakes that occurred at depths shallower than 200 km between 1963 and 1991 these events took place in two regions where the double seismic zone terminates and single seismic zones emerge precise earthquake source parameters including focal depths and mechanisms are determined from a careful analysis of a mixture of analog and digital p and sh waveforms recorded at teleseismic distances our results indicate that the double seismic zone is confined between approximately 46 and 53 n northeast of 53 n intermediate depth earthquakes 70 200 km show only downdip compressional focal mechanisms furthermore the wadati benioff zone there seems to be a continuation of the upper compressional layer of the double seismic zone thus the changeover from a double to a single seismic zone in this case is simply disappearing of the lower extensional sheet of the double seismic zone similarly southwest of 46 n intermediate depth earthquakes show only downdip extensional mechanisms the single seismic zone is unlikely to consist of two layers of extensional seismicity in fact nowhere along the entire arc did we observe two distinct seismogenic layers with a single state of strain in other words these extensional events did not occur in the same seismogenic structure as the upper compressional layer of the double seismic zone however we cannot ascertain whether the single seismic zone in extension is simply a southwestern extension of the lower extensional layer of the double seismic zone lateral variation from a single zone under extension in the southwest to a double seismic zone and finally to a single zone under compression cannot be accounted for by conventional models of a double seismic zone on the other hand the addition of a component of stress transmitted from depths greater than 250 km can explain the observed lateral variation in the state of strain at intermediate depths because the slab between depths of 250 450 km is under extension in the southwestern end of the arc whereas the rest of the arc is under compression vanishing of the lower extensional layer of the double seismic zone toward the northeastern end of the arc is interpreted as a result of compression transmitted from below overshadowing extension caused by effects such as unbending a similar reasoning applies to the southwestern termination of the double seismic zone

deformation history of the southwestern taiwan foreland thrust belt insights from tectono sedimentary analyses and balanced cross sections abstract analysis of the basin scale and local deformation history and estimates of shortening rates are performed in the southwestern foreland thrust belt of the plio pleistocene taiwan orogen timing of deformational events is constrained using both syntectonic unconformities and other markers of syndepositional tectonic activity three cross sections are investigated and balanced intermediate shortening rates are estimated based on balanced and restored cross sections the total amount of shortening along studied balanced sections results from superimposed thin skinned and thick skinned styles a shallow aseismic thick decollement is located to depth of 6 8 km at the base of thick late miocene shales that behaves as ductile layer the occurrence of a deep seismic decollement to depth of 12 15 km is confirmed by seismicity fault plane solutions of focal mechanisms and a seismic reflection profile providing evidence of reactivated and inverted inherited extensional features beneath the coastal plain three major tectonic stages are recognized within the southwestern foothills which are correlated with local thrusting events the first stage begins at 5 ma with the onset of the taiwan collision and ends at 2 1 6 ma it reflects an initial submarine accretionary stage associated with moderate tectonic activity the second major stage occurs during late pliocene early pleistocene this intermediate collision stage sealed the evolution of the southern foreland thrust wedge from submarine to continental environment in relation with exhumation and denudation of the central range the last stage of deformation begins in the middle pleistocene and continues until present it is associated with high rates of sedimentation and most of the thrusts are activated or reactivated at that time and out of sequence thrusting occurs the collision started synchronously near 5 ma throughout the southwestern taiwan total and intermediate shortening rates are more or less constant through time in studied sections that may suggest that long term deformation in the taiwan foreland thrust belt occurred more or less at steady state since 5 ma the southward decrease of rates of shortening from 3 6 to 4 8 mm year in northern sections to 2 2 mm year in southern section might be related to both influence of indentation by the peikang high to the north and southward transition at the scale of the southwestern foothills from collision to subduction setting the effect of the obliquity of the convergence resulting in regional southward migration of the collision cannot be definitely inferred from this analysis our results place new kinematic and time constraints which must be taken into account in forthcoming geodynamic models of the taiwan collision

intracrustal detachment within zones of continental deformation geologic mapping of active to recently active geologic structures in panamint valley western united states and in the haiyuan region and northern tibetan plateau china suggests detachment in the middle and lower crust on a scale of tens of kilometres to at least several hundred kilometres detachment occurs similarly in predominantly extensional panamint valley and in predominantly compressional china environments it involves structures with displacements of more than 10 km and displacement rates of more than 3 mm yr perhaps more than 10 mm yr the steeply dipping strike slip faults present in all three areas hunter mountain fault haiyuan fault and altyn tagh fault zone terminate in zones of extension or compression and geometric relations indicate that all structures including strike slip faults are thin skinned and restricted to the upper crust in panamint valley and in the haiyuan region deformation within these systems can be reconstructed in three dimensions displacement on the strike slip faults is absorbed by extension or compression occurring at the termination of the faults so strike slip displacement is roughly equal to and in the same direction as shortening or extension we propose that left slip on the altyn tagh fault zone in northern tibet is similarly absorbed by shortening southeast of the fault zone within the qaidam basin and the nan shan

improved seismic tomography offshore northeastern taiwan implications for subduction and collision processes between taiwan and the southernmost ryukyu summary to improve the resolution in imaging the 3 d vp and vp vs structures in the taiwan region especially offshore eastern taiwan we combine the arrival times from eleven ocean bottom seismometers obss and from seismic stations of the japan meteorological agency jma with those from the permanent stations of taiwan central weather bureau seismic network cwbsn and taiwan strong motion instrumentation program tsmip by doing so we have obtained a new 3 d tomographic model for taiwan and its surrounding regions with a better resolution especially in the area offshore northeastern taiwan we also used this new tomography model to relocate the hypocentres of the earthquakes in northeastern taiwan and determined the focal mechanisms of relatively large events our tomography results indicate that in the region northeast of taiwan the subducting oceanic philippine sea plate beneath the eurasian continental lithosphere is characterized by a high vp layer surrounded by lower vp areas this p wave velocity characteristics of the subducting slab provides a better constraint on the geometry of the subduction interface especially in its shallower portion in the hangingwall eurasian block above the subduction interface a vertically elongated high vp vs body appears to originate from the interface at depths between 100 and 140 km we suggest that it represents the partially melted materials that are related to the magmatic activity of the ryukyu volcanic arc

brittle frictional mountain building 2 thermal structure and heat budget this paper describes a simple thermal model of an actively deforming critically tapered fold and thrust belt the model determines the steady state temperature distribution and heat flow as well as the pressure temperature time histories of rocks that outcrop at the surface the main parameters controlling the thermal structure are the accretion and erosion rates the undisturbed geothermal gradient at the toe and the amount of frictional heating both shear heating on the decollement fault and internal strain heating within the deforming brittle wedge are incorporated in a mechanically consistent manner and they dominate the effect of radiogenic heating except in fold and thrust belts with significantly overpressured pore fluids the mean stresses temperatures and surface heat flow all increase with an increase in the basal and internal coefficients of friction and this dependence is used to constrain the level of friction on the decollement fault beneath the steady state fold and thrust belt in taiwan rocks outcropping in the core of the central mountain range of taiwan experience maximum theoretical temperatures in excess of 400 c and maximum mean pressures in excess of 500 mpa if the coefficient of basal friction is b 0 5 qualitatively these conditions are in good agreement with the observed high greenschist facies metamorphism the theoretical surface heat flow which increases from 95 mw m2 at the front of the fold and thrust belt to 240 mw m2 at the rear is in excellent agreement with the results of a recent geothermal survey of taiwan and theoretical cooling histories are in good agreement with fission track and other geochronologic studies taken together these results provide strong evidence that sliding on the basal decollement fault beneath taiwan is governed by a coefficient of friction in the range of typical laboratory measurements b 0 5 0 2 approximately 35 of the total surface heat flux of 3 gw is heat conducted into the base of the wedge from the top of the basal decollement fault and somewhat more than 30 is heat advected into the toe by accretion the remaining heat is generated internally about 25 by internal strain heating and about 10 by radiogenic heating either an increase in the coefficient of basal friction b or a reduction in the pore fluid pressure ratio b leads to an increase in the surface heat flow because of the increased frictional heating within the wedge and on the basal decollement fault the overall balance of energy in a steady state fold and thrust belt is described by the equation e w g q where is the rate at which both mechanical and heat energy are added from external sources w g is the rate at which work is performed against gravitational body forces in a reference frame attached to the overriding plate and q is the rate at which waste heat flows out of both the upper and lower boundaries the total power input into the taiwan fold and thrust belt is approximately 4 2 gw the mechanical work being done on the base and front of the fold and thrust belt accounts for 3 of these 4 2 gw in addition 0 9 gw of heat are being advected from the subducting plate into the toe by accretion the remaining 0 3 gw are being supplied by in situ radiogenic heating the outgoing energy is dominated by the 3 gw of heat conducted out the top in the surface heat flow however another 0 8 gw are conducted down beneath the rear portion of the basal decollement fault to heat the underlying subducting slab only 0 4 of the incoming 4 2 gw do useful mechanical work against gravity within the wedge the efficiency of brittle frictional mountain building in taiwan is therefore 10

estimates of modern arc parallel strain rates in fore arcs deflections of slip vectors of interplate thrust earthquakes from expected directions are used to estimate arc parallel strain rates within the overriding plates at the world9s major convergent plate margins arc parallel extension strain rates between 10 8 yr and 10 7 yr and significant at 2 standard deviations are observed in the fore arcs of the aegean aleutian mariana sumatran southern kuriles new hebrides scotia and southern central american lat 8 to 12 n subduction zones and the himalayas northern chile lat 17 to 31 s and central america lat 11 to 18 n show arc parallel compression available geologic and geodetic estimates of fore arc slip and strain rates agree within a factor of two with slip vector estimates arc parallel strain in fore arcs is rapid enough to produce geologically significant effects such as unroofing of high grade metamorphic rocks and disruption of transported fore arc terranes fore arcs deform even where convergence is perpendicular to curved margins demonstrating that head on subduction can produce a three dimensional strain field

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

the taiwan ryukyu subduction collision complex folding of a viscoelastic slab and the double seismic zone 1 the termination of the ryukyu trench against eurasia and the oblique subduction of the philippine sea plate create a subduction collision complex offshore taiwan which has not previously been elucidated in detail we combine traveltime data from the seismic networks in taiwan and japan to better illuminate how the subducting ryukyu slab deforms in this subduction collision zone more than 5000 events recorded by both networks were relocated with the double difference method using an optimal regional one dimensional velocity model the offshore seismicity indicates that the double seismic zone with a gap of 15 20 km exists in the subducting slab in the depth range of 40 80 km focal mechanisms suggest that the double seismic zone is caused by east west compression resulting from oblique convergence the improved hypocentral locations for the first time reveal folding of the slab into a horizontal curvature larger in magnitude than and opposite in sign to that of the ryukyu trench in the depth range 50 100 km the anomalous curvature together with the focal mechanisms suggests that the slab folds against the eurasian lithosphere and that this deformation cannot be fully elastic we model this deformation mode as the developing instability of a viscoelastic maxwell layer embedded in a viscous medium the characteristic wavelength of the instability i e 250 km is consistent with folding of a slab whose viscosity is 100 times higher than that of the surrounding mantle for an along strike elastic membrane strain as small as 0 01 or more than 3 orders of magnitude higher if 5 elastic strain is allowed

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

post collisional magmatism around northern taiwan and its relation with opening of the okinawa trough being part of an active mountain belt formed by oblique collision of the luzon arc with asia northern taiwan and the offshore islets are marked by a series of latest pliocene quaternary volcanoes whose eruptions have been conventionally ascribed to westward propagation of the ryukyu volcanic arc on the basis of new geochemical data along with supporting geologic and geophysical evidence we propose instead that this young volcanism resulted from post collisional lithospheric extension in the northern taiwan mountain belt and the mantle source regions involved in the melt generation have been significantly modified by the nearby ryukyu subduction related processes magmas thus produced through the northeast to the southwest in the northern taiwan volcanic zone ntvz vary from low k to calc alkaline and then shoshonitic compositions such a spatial geochemical variation characterized by southwestward increase in enrichments of potassium and incompatible trace elements runs subparallel to the southwestern part of the present day ryukyu trench the geochemical variation that can be explained by southwestward decrease in degrees of partial melting of the mantle sources is reconciled with a southwestward weakening extensional regime observed in the ntvz the post collisional extension in northern taiwan furthermore might have played a role in reactivation of the opening in the middle okinawa trough and gave way to its rapid southwestward propagation with associated development of the ryukyu subduction zone west of124oe this nascent subduction in turn resulted in abundant submarine volcanoes which delineate an embryo volcanic front along the southern margin of the trough therefore the southwestern okinawa trough is not a fore arc basin as had been previously alleged but represents an atypical back arc basin which developed broadly synchronously or prior to its arc trench system in this particular collision extension subduction tectonic environment 1999 elsevier science b v all rights reserved

magnitudes of large shallow earthquakes from 1904 to 1980 abstract in order to obtain a uniform magnitude catalogue surface wave magnitudes ms and broad band body wave magnitudes mb have been determined for large shallow earthquakes from 1904 to 1980 in making the catalogue homogeneous the author consistently adheres to the original definitions of ms and mb given by gutenberg 1945 and gutenberg and richter 1956 the determinations of ms and mb are all based on the amplitude and period data listed in gutenberg and richter s unpublished notes bulletins from stations worldwide and other basic information mb is measured on broad band instruments in periods of 8 s consistency of the magnitude determinations from these different sources is carefully checked in detail more than 900 shallow shocks of magnitude 7 and over are catalogued the meaning of the magnitude scales in various catalogues is examined in terms of ms and mb most of the magnitudes listed by gutenberg and richter 1954 in their seismicity of the earth are basically ms for large shocks shallower than 40 km but are basically mb for large shocks at depths of 40 60 km the surface wave magnitudes given by earthquake data reports are higher than ms by 0 2 unit whether the combined horizontal amplitude or the vertical amplitude is used mb and the currently used 1 s body wave magnitude are measured at different periods and should not be directly compared

cause and effect among thrust and normal faulting anatectic melting and exhumation in the himalaya the transition from compression to extension in a regime of overall convergence in a narrow mountain belt requires a change with time in at least one of the following horizontal normal stress applied to the belt by its surroundings the distribution of shear traction on the base of the belt the constitutive relation between stress and strain rate of material in the belt itself the surface height of the belt the following processes could have caused one or other of these changes in the himalaya in miocene time and thus could have triggered the onset of extension by normal faulting on planes parallel to thrust faults that were active at the same time the shear stress on the thrust fault at the base of the himalaya could have decreased for instance because of melting on the fault plane alternatively the creation of a new deeper thrust fault and hence a new locus of convergence could have initiated extensional faulting if the distribution of shear stress on the new thrust fault were unable to support the himalaya above it the himalaya could have been weakened by an increase in pore fluid pressure for instance owing to metamorphic dehydration reactions within the wedge the himalaya could have been elevated for instance following the detachment of a tongue of mantle lithosphere from beneath it or by the underthrusting of a lens of material beneath its base the apparent simultaneity of thrust and normal faulting on nearly parallel planes of anatexis and of exhumation in the himalaya in miocene time can be explained if normal faulting was responsible for most of the exhumation and through decompression melting for anatexis existing data suggest that the most probable cause for miocene normal faulting was either uplift of the himalaya or the cessation of slip on the main central thrust at the base of the orogenic wedge and the creation of a new deeper thrust fault and a new base to the deforming himalaya 1993 the geological society

quaternary extension in southern tibet field observations and tectonic implications we summarize evidence for quaternary and active faulting collected in the field during three sino french expeditions to southeastern tibet 1980 1982 detailed mapping of quaternary and active faults as well as microtectonic measurements indicate that normal faulting has been the dominant tectonic regime north of the himalayas in the last 2 0 5 m y the maximum horizontal principal stress in south tibet appears to be only the intermediate principal stress 2 1 being vertical south of the chord joining the eastern and western syntaxes of the himalayan arc extensional strains are principally localized within seven regularly spaced rift zones three of which have been studied in some detail the extension direction is determined to be n96 e 7 mainly from statistical averaging of strikes of newly formed normal faults throw rates on normal faults are evaluated for different time spans 2 0 5 m y 60 40 kyr and 10 2 kyr b p using structural and topographic reliefs as well as synglacial and postglacial vertical offsets the rate of quaternary extension is about 1 m y 1 along an 1100 km long ese traverse across south tibet this corresponds to a spreading rate of 1 0 6 cm yr this rate and the divergent horizontal projections of slip vectors of earthquakes along the himalayan front constrain the rate at which rigid india underthrusts southern tibet to be 2 1 cm yr although most of the normal faults appear to be independent of and nearly orthogonal to the e w mesozoic tertiary tectonic fabric the yadong gulu rift appears to be guided for over 130 km by the older oblique ne sw nyainqentanglha range and fault zones along it this reactivated zone is the most prominent left lateral strike slip fault system in se tibet excepting this zone and the vicinity of the se extremity of the karakorum fault quaternary strike slip faulting is rare in south tibet i e south of the chord between the syntaxes of the himalayan arc north of the chord the tectonic style is different there minor conjugate strike slip faulting is widespread and appears to control quaternary normal faulting which is more diffuse and subdued than in south tibet along the chord the presence of a major zone of active right lateral en echelon strike slip faults karakorum jiali fault zone probably reflects the greater facility of eastward extrusion in north central tibet in response to the northward push of india the eastern himalayan syntaxis may be an obstacle to such extrusion movements south of the chord

double seismic zones and stresses of intermediate depth earthquakes data from japanese local seismograph networks suggest that the stresses in double seismic zones are in plate compression for the upper zone and in plate tension for the lower zone the stresses do not necessarily appear to be down dip it may therefore be possible to identify other double seismic zones on the basis of data which indicate that events with differing orientations of in plate stresses occur in a given segment of slab a global survey of published focal mechanisms for intermediate depth earthquakes suggests that the stress in the slab is controlled at least in part by the age of the slab and the rate of convergence old and slow slabs are under in plate tensile stresses and the amount of in plate compression in the slab increases with increasing convergence rate or decreasing slab age young and fast slabs are an exception to this trend all such slabs are down dip tensile since these slabs all subduct under continents they may be bent by continental loading double seismic zones are not a feature common to all subduction zones and are only observed in slabs which are not dominated by tensile or compressive stresses unbending of the lithosphere and upper mantle phase changes are unlikely to be the causes of the major features of double zones although they may contribute to producing some of their characteristics sagging or thermal effects possibly aided by asthenospheric relative motion may produce the local deviatoric stresses that cause double zones

numerical modeling of intraplate deformation simple mechanical models of continental collision we study a model of continental collision in which one of the continents acts as a rigid die indenting the other plate which flows as an incompressible viscoplastic medium we consider two extreme cases of plane deformation 1 plane strain which corresponds to an infinitely thick lithospheric plate and 2 plane stress corresponding to a very thin plate deformation of the lithosphere a thick plate should be intermediate between those extremes we found that the flow in the plane strain case is quite similar to that obtained by slip line theory the plane stress results are quite different since in this case most of the plate shortening is taken up by the thickening of the lithosphere we also explored the role of boundary conditions on the flow in particular the role of the side walls containing the flow of the lithosphere in the case of a free lateral boundary the main feature is a flow of matter toward this free wall and a s like pattern for the horizontal stress field for a rigid wall on the other hand the plane strain and the plane stress results are quite different in the first case there is a large return flow on the sides of the punch the material being extruded along the only free surface available in the plane stress case the return flow disappears and the material displaced by the penetration of the die tends to thicken the plate the role of a nonlinear constitutive relation is studied for power law creep as the power of the flow limit increases the flow retains its general features but the deformation localizes creating sharper contrasts between high and low strain rate areas in plane stress the effect of nonlinearity is to increase the contrasts in vertical motion available data for asia are discussed in the light of the new results

three dimensional patterns of seismic deformation in the taiwan region with special implication from the 1999 chi chi earthquake sequence abstract we compile a comprehensive list of hypocentral locations and source parameters for earthquakes occurred in the taiwan region the collected database is used to construct 3d patterns of seismic deformation that are associated with various tectonic processes such as subduction and collision between the philippine sea plate and eurasia specifically we divide the entire taiwan region into a 3d grid of 10km interval and the amount of seismic moment released at each point is the sum of all individual events in the vicinity in case when the physical dimension of an earthquake is larger than 10km the corresponding seismic moment is distributed across multiple grid points this quantitative approach gives a better depiction on many first order features in the region the most intriguing one is that a significant deficit in the total amount of released seismic moment can be clearly identified around the chi chi source region before the big earthquake occurred in 1999 the deficit trough is then filled by the seismic moment of the chi chi earthquake sequence this deficit then fill pattern in a region s seismic moment distribution can be used as an indicator to pinpoint the locations of large earthquakes in the foreseeable future following this argument we suggest that disastrous earthquakes with magnitudes comparable to that of the chi chi earthquake are due for the two regions adjacent to the chi chi source area i e the miaoli hsinchu domain to the north and the chiayi domain to the south

layered anisotropy within the crust and lithospheric mantle beneath the sea of japan abstract continental rifting during the oligocene to mid miocene caused the opening of the sea of japan and the separation between the japanese islands and the eurasian plate the tectonic evolution in the sea of japan is important for understanding the evolution of back arc regions in active convergent margins here we use data from the seismic stations surrounding the sea of japan to map the rayleigh wave azimuthal anisotropy in the crust and lithospheric mantle beneath the sea of japan we explore the variations of rayleigh wave phase velocity beneath the sea of japan in a broad period range 30 80 s rayleigh wave dispersion curves are measured by the two station technique for a total of 231 interstation paths using vertical component broad band waveforms at 22 seismic stations around the sea of japan from 1411 global earthquakes the resulting maps of rayleigh wave phase velocity and azimuthal anisotropy allow the examination of azimuthal anisotropy at specific periods they exhibit several regions with different isotropic and anisotropic patterns the japan basin displays fast velocities at shorter periods 30 and 40 s with nne ssw anisotropy whereas at 60 s and longer the velocities become slow even if the anisotropy remains ne sw the east china sea shows fast velocities at all periods 30 80 s with constant nw se anisotropy trench normal anisotropy beneath the japanese islands is found at short periods 30 40 s and become trench parallel at periods of 60 s and longer overall our model resolves two layers of anisotropy the shallowest and deepest layers being potentially related to frozen deformation due to recent geodynamic events and asthenospheric flow respectively

earthquake focal mechanisms deformation state and seismotectonics of the pamir tien shan region central asia this paper presents a review of the seismotectonics of the pamir tien shan collision zone in the garm region tajikistan based on geological structure seismicity and focal mechanism solutions the region is dominated by horizontal compression manifested by imbricate low angle thrust faults that separate the upper crust into a series of tectonic sheets these thrust systems verge northward from the tajik depression toward the southern tien shan and southward toward the northern pamir the pattern of seismicity across the region suggests that similar low angle thrust faults exist within the crystalline basement as well in order to reconstruct the present day stress strain state of the region we used data from over 15 000 earthquake focal mechanism solutions for small earthquakes m 1 0 gathered over the 27 year period 1963 1989 the method of reconstruction of the seismotectonic deformation std field involves dividing the data set into small spatial windows and summing the individual focal mechanism solutions to form an average earthquake mechanism tensor for each cell the std state of that cell can then be represented by the orientation of the principal axes of compression and tension the intensity or relative uniformity of the average mechanism and the lode nadai coefficient which defines the relation between the magnitudes of the principal strain components we observe throughout the study area a general predominance of subhorizontal compression manifested in a mixture of thrust and strike slip deformation normal fault deformation is observed in only two small portions of the study area the orientation of principal compression varies from nearly n s in the northern pamir and the southern tien shan to nw se in the peter the first range in general the std field has a clearly developed mosaic structure defined by spatial groupings of nearly uniform orientations of compression and tension axes the boundaries of these groups cannot in most cases be directly correlated with known geological contacts the std structure also shows some variation with depth with more or less uniform std orientation in the upper crust 0 14 km depth changing significantly at greater depths examination of various subsets of the focal mechanism catalog shows the spatial structure of the std field to be largely stable with respect to both time and magnitude the reliability of the reconstructed std field was verified using three approaches 1 comparison of individual focal mechanism determinations for a large subset of the data using both soviet and u s algorithms 2 comparison of std reconstruction results using the two independent focal mechanism catalogs and 3 examination of the effect of spatial sampling on the results the dominance of generally north oriented subhorizontal compression is interpreted to be primarily the result of convergence between the pamir and tien shan ranges in turn caused by the ongoing collision of india and eurasia

active hydrothermal mounds in the okinawa trough backarc basin japan abstract dive studies by the shinkai 2000 were made in september 1984 and july 1986 in the central axial rift of the middle okinawa trough a possible active backarc basin in a nascent stage hydrothermally active mounds were discovered on a small knoll 27 34 4 n 127 08 6 e directly south of the area where a high heat flow higher than 1 w m2 had been observed previously the largest mound 5 6 m high and 15 25 m wide has several peaks and ridges along which yellow materials are exposed shimmering water was noted at several places on the yellow materials the rest of the mound is thinly covered by brown materials the temperature of the ejecting water was higher than that of the ambient seawater by 2 3 c and a 40 cm long thermometer inserted into the mound recorded a temperature of 20 50 c analysis of the ejected water showed a methane content of about 200 nl kg it is possible that further dive studies will reveal larger scale hydrothermal activity in the area

the crustal deformation of the ilan plain acted as a westernmost extension of the okinawa trough abstract we analyze the current deformation in ilan plain of northeastern taiwan considered as the western extension of okinawa trough based on 34 geodetic sites observed between 2002 and 2006 the station velocities are between 1 9 mm yr and 43 0 mm yr in a direction between 38 and 143 relative to stable continent margin of eurasian plate this largest velocity is moving southeastward 143 at a rate of 43 mm yr was enhanced by a seismic doublet that occurred on march 6 2005 the station velocities at the northern margin of ilan plain are insignificant with a magnitude of 1 9 mm yr to 7 5 mm yr on the contrary the station velocities increase approximately southward along the southern margin of ilan plain the stations are moving roughly e w with a magnitude of 5 8 mm yr to 21 8 mm yr further southward the station velocities increase to a value of 43 0 mm yr in a direction of 143 the most prominent feature of the strain distribution patterns in ilan plain are the remarkable extensional strain rates observed in central and southern flank of the ilan plain the largest extensional rate is found 2 66 strain yr in 147 with a shortening rate in the perpendicular direction 1 28 strain yr in 57 indicating transtensional deformation mode in the northeastern flank of ilan plain insignificant or minor deformations are observed overall the extensional rates in this area clearly increase from north to south the extension directions in study area mostly trend e w to nw se directions using the velocity field of 34 gps stations and its strain rate filed we suggest the lateral extrusion process towards mechanically weak domain adjacent ryukyu subduction zone the lateral extrusion is facilitated by the opening of the okinawa trough we conclude that the northeast taiwan is subjected to both the compressional force exerted by the indentation of philippine sea plate and the tensile force induced by trench retreat related to the suction force at the ryukyu trench due to the interaction of two forces the stress regime rapidly changes from pure compression in the taiwan collision belt to transtension near ilan plain and its offshore area

delineating the rupture planes of an earthquake doublet using source scanning algorithm application to the 2005 march 3 ilan doublet northeast taiwan summary correct identification of the fault plane s associated with an earthquake doublet is a very challenging problem because the pair of events often occurs in close space and time with almost the same magnitude most long period waveforms of an earthquake doublet are severely tangled and thus unsuitable for conventional waveform inversion methods in this study we try to resolve this issue by utilizing the recently developed source scanning algorithm ssa the ssa systematically searches the model space for seismic sources whose times and locations are most compatible with the observed arrivals of large amplitudes on seismograms the identification of a seismic source is based on the brightness function which is defined as the summation of the normalized waveform amplitudes at the predicted arrival times at all stations by illuminating the spatiotemporal distribution of asperities during an earthquake s source process we are able to constrain the orientation of the rupture propagation that in turn leads to the identification of the fault plane a series of synthetic experiments are performed to test ssa s resolution under various scenarios including different directions of rupture propagation imperfect station coverage and short origin time difference between the two events of a doublet because only short period records are needed in the analysis the proposed method is best suited for an earthquake doublet with a short time gap between the two events using the 2005 ilan doublet the origin time difference is only 70 s that occurred in northeast taiwan as an example we show that the trace of the brightest spots moves towards the west and infer the e w striking plane to be the actual fault plane

crustal structure of the southernmost ryukyu subduction zone obs mcs and gravity modelling summary in 1995 a combined ocean bottom seismometric obs and multichannel seismic mcs survey with strong air gun shots was carried out in the southernmost ryukyu subduction zone a crustal velocity structure constructed from the layer stripping monte carlo inversion of three obs mcs profiles and the associated density models inverted from gravity data in the sw end of the ryukyu arc trench system are presented parallel to the arc in the southernmost ryukyu subduction system the obs mcs profiles show sedimentary layers of the hoping nanao and east nanao forearc basins from west to east warping of the arc basement and buckling of the subducted slab beneath the hoping basement rise the arc parallel variation of the crustal structure may result from increasing lateral compression westward due to oblique subduction of the philippine sea plate and collision with the luzon arc near the northwestern edge of the forearc region northward subduction and arc parallel compression of the slab also have generated thrust faulting along the subduction interface and strike slip faulting within the subducted slab respectively on 2002 march 31 an earthquake with a moment magnitude of 6 84 was induced by buckling of the subducted slab and strongly affected cities within an epicentral distance of 100 km the velocity interface models the density models and the focal mechanisms presented in this paper therefore suggest that earthquakes induced by slab buckling or arc parallel compression have been stronger but less frequent than those generated by northward subduction in the ryukyu seismogenic zone off taiwan

a model for the termination of the ryukyu subduction zone against taiwan a junction of collision subduction separation and subduction boundaries 1 the nw moving philippine sea plate psp collides with the eurasian plate eup in the vicinity of taiwan and at the same time it subducts toward the north along sw ryukyu the ryukyu subduction zone terminates against eastern taiwan while the ryukyu trench is a linear bathymetric low about 100 km east of taiwan closer to taiwan it cannot be clearly identified bathymetrically owing to the deformation related to the collision making the location of the intersection of the ryukyu with taiwan difficult to decipher we propose a model for this complex of boundaries on the basis of seismicity and 3 d velocity structures in this model the intersection is placed at the latitude of about 23 7 n placing the northern part of the coastal range on eup as psp gets deeper along the subduction zone it collides with eup on the taiwan side only where they are in direct contact thus the eurasian plate on the taiwan side is being pushed and compressed by the nw moving philippine sea plate at increasing depth toward the north offshore of northeastern taiwan the wedge shaped eup on top of the ryukyu subducting plate is connected to the eup on the ryukyu side and coupled to the nw moving psp by friction at the plate interface the two sides of the eup above the western end of the subduction zone are not subjected to the same forces and a difference in motions can be expected the deformation of taiwan as revealed by continuous gps measurements geodetic movement along the east coast of taiwan and the formation of the hoping basin can be understood in terms of the proposed model

heat flow in the southwestern okinawa trough abstract lying to the southeast of the east china sea continental shelf the okinawa trough is believed to have been formed by crustal spreading behind the ryukyu island arc it has the highest mean heat flow among the pacific marginal basins the mean of 18 reliable measurements is3 34 2 75hfu the next highest mean heat flow for a marginal basin is2 53 1 60hfu for the fiji plateau the high mean value for the okinawa trough may be partially attributable to a sampling bias because most measurements were made in or near the well sedimented central rift or graben the reported values range from 0 36 to 10 49 hfu with low values probably resulting from hydrothermal circulation in an area of large relief the high heat flow suggests that the age of the okinawa trough may be as young as 2 ma such an estimate is in gross agreement with the inferred plio pleistocene age for the drifting in the okinawa trough the opening of the okinawa trough may continue onto taiwan but the relation between the opening of the trough and the occurrence of a nearby onshore geothermal resource area remains unsolved

seismotectonics of taiwan abstract high quality seismicity data and focal mechanism solutions obtained during 1973 1983 by the permanent taiwan telemetered seismographic network and several temporary local seismographic networks are used for a detailed study of the seismotectonics of the taiwan area seismicity distribution in southern taiwan clearly reveals an east dipping benioff zone which has a thickness of about 30 km and begins to deepen along 121 e at a dip angle of 55 60 the leading edge of this benioff zone reaches a depth of about 180 km between 21 n and 22 n but tapers off to a shallower depth of about 100 km from 22 n to 23 n the presence of this seismic zone implies that subduction of the south china sea plate under the philippine sea plate extends from luzon northward to about 23 n the position of the northern boundary of the south china sea plate as tentatively defined according to the seismicity distribution passes through southern taiwan from the offshore area in the taiwan strait west of kaohsiung in an east northeast direction to the taitung area where a triple junction probably lies seismicity is found to disappear abruptly below a certain depth in many parts of taiwan this phenomenon may be attributed to the frictional to quasiplastic transition in the crust or upper mantle comparison of shallow seismicity with surface faults and fractures shows that all areas of active shallow seismicity are marked by densely developed faults and fractures however the converse is not necessarily true this may be partly due to the relatively short duration of seismicity data and partly due to excessive weakening of some of the severely faulted and fractured areas finally focal mechanism solutions for west central taiwan and the kuangfu fuli area in eastern taiwan predominantly show a maximum horizontal compression in the se nw direction which can be related to collision between the eurasian and philippine sea plates however focal mechanism solutions for both the hualien area in eastern taiwan and the tainan area in southwestern taiwan show remarkable irregularities which may result from local tectonic complexities

tomographic imaging of lithospheric structures under taiwan abstract tomographic images of the crustal and mantle velocity structures under taiwan are obtained by simultaneous inversion of local earthquake p wave arrival times for hypocenters and p wave velocity structures in northern taiwan a high velocity zone coinciding with the wadati benioff zone can readily be identified as the subducted philippine sea plate the imaged zone dips toward the north at an angle of 40 from a depth range of 20 55 to 100 130 km an upper mantle low velocity wedge ranging from 40 to 80 km in depth exists above the subducted slab above this wedge is the ilan plain of northern taiwan which lies at the west end of the okinawa trough a well recognized back arc basin the crustal velocities under the plain are also relatively low the well defined high velocity zone and the low velocity wedge provide some constraints on the thermal structures of the subduction system under northern taiwan in tomographic images across the central section of taiwan thickening of the crust and up arching of the lower crustal materials under the central range are commonly observed the crust under the western foothills region is clearly thinner and the near surface low velocity layers are well developed the structures under the central range show that although the taiwan orogeny is quite young a root deeper in the north and shallowing to the south has formed the results of our tomography show that a significant portion of the lithosphere is involved in the taiwan orogeny

strength of plate coupling in the southern ryukyu subduction zone abstract understanding the strength of a plate coupling is critical for assessing potential seismic and tsunamic hazards in subduction zones the interaction between an overriding plate and the associated subducting plate can be used to evaluate the strength of plate coupling by examining the mantle lithospheric buoyancy here we calculate the mantle lithosphere buoyancy across the northern portion of the southern ryukyu subduction zone based on gravity modeling with the constraints from a newly derived p wave seismic velocity model the result indicates that the strength of the plate coupling in the study area is relatively strong which is consistent with previous observations in the southernmost ryukyu subduction zone because few large earthquakes mw 7 have occurred in the southern ryukyu subduction zone a large amount of energy is locked and accumulated by plate coupling that could be released in the near future

present day interseismic surface deformation along the longitudinal valley eastern taiwan from a ps insar analysis of the ers satellite archives 1 in taiwan about one third of the lithospheric plate convergence between eurasia and the philippine sea plate is accommodated on the eastern coast across the narrow longitudinal valley lv the longitudinal valley fault lvf is the main seismically active fault zone in this region however the spatial distribution of ground deformation due to interseismic loading on locked fault segments remains widely unknown to address this problem we use a permanent scatterer ps interferometric synthetic aperture radar insar approach for processing the ers satellite archives the ps mean velocity maps allow accurately mapping and quantifying the ground deformation all along the lv assuming that close to the surface the fault plane along the central lvf creeping segment dips 55 to the east we derive a slip vector at the surface of 25 mm yr with a mean rake of 70 in agreement with our ps observations shallow seismicity indicates a clear change in the fault mechanical behavior of the creeping segment around fuli town finally we propose a combination of the ps velocity fields for estimating a 3 d map of ground deformation in the southern part of the lv the vertical component reveals an uplift of the coastal range with respect to the lv of 10 mm yr moreover between taitung and luyeh the vertical deformation is essentially accommodated along the luyeh strand finally north of rueisuei the coastal range uniformly subsides at 15 mm yr relative to the western side of the lv however a local uplift pattern is detected around latitude 23 42

kinematic model of active deformation in central asia the velocity field of present day deformation in central asia is modelled using a set of four rotating blocks siberia tarim tibet india on a spherical earth a best fit is inverted on the basis of estimated shortening rates across the main thrust zones himalayas tien shan and measured slip rates along the principal strike slip faults altyn tagh and karakorum separating those blocks the fit to the data implies that nearly all the present convergence between india and asia can be accounted for by slip partitioning on these four zones with as much as 50 absorbed by northeastwards extrusion of tibet this suggests that localised deformation governs the present mechanical behaviour of the central asian lithosphere

tectonics of western tibet between the tarim and the indus a tarim to indus traverse provides insight into the tectonics of western tibet the kunlun was the site of a mid paleozoic collision at least three phases of post paleozoic accretion have thickened the blanket of sediments that covers western tibet sizeable parts of western qiangtang have remained stable however since the mid mesozoic since the neogene deformation and volcanism have been localized near the edges of the plateau strike slip motion along the karakorum and altyn tagh faults has been coeval with overthrusting in the himalayas and kunlun such slip partitioning and the volcanism appear to result simply from northward subduction of india and southward subduction of the tarim as tibet is extruded eastwards by india s penetration into asia

seismogenic patterns in the taiwan region insights from source parameter inversion of bats data abstract we systematically invert source parameters of 96 earthquakes that occurred in the taiwan region between july 1995 and june 1998 using waveforms recorded by the newly established broadband array in taiwan for seismology bats the results are utilized to delineate the seismogenic patterns associated with the regional tectonic processes between the eurasia and philippine sea plates in general high seismicity is observed at five locations including the nan ao basin 24 n 122 4 e near the hualien area 24 2 n 121 7 e east of the longitudinal valley and within the philippine sea plate the forearc region to the west of lanshu and the okinawa trough events with focal depths less than 25 km scatter across a wide region from the taiwan strait to the ryukyu and luzon arcs whereas events between 25 and 50 km concentrate offshore east of taiwan most earthquakes that occurred beneath the nan ao basin are consistent with the northward subduction of the philippine sea plate along the plate interface events offshore east of the longitudinal valley and within the philippine sea plate are dominated by a compressional regime along the nw nww directions whereas the lanshu lutao forearc is characterized by e w compression normal faults showing n s extension are observed in the okinawa trough where we observed a successive rotation of t axes from n s to nw se as the epicenters approach taiwan indicating a possible interaction between the extensional and compressional strain regimes there most events between 50 and 65 km depth within the subducted philippine sea plate beneath ryukyu are in downdip extension with lateral compression approximately parallel to the strike of the slab the depth where the state of strain switches from downdip extension to downdip compression is constrained at 110 125 km for the southernmost ryukyu slab based on our results we propose that the subducted philippine sea plate beneath ne taiwan plays a significant role in the overall orogenic process of taiwan the lithospheric collision in taiwan should consist of two major components one is the well recognized arc continent collision that dominates the central and southern taiwan and another the slab continent collision that reigns in northern taiwan a cluster of earthquakes near hualien suggests that an incipient west or northwestward subduction of the philippine sea plate is currently taking place there alternatively they can be interpreted to suggest that the steeply east dipping longitudinal valley fault may have transformed near hualien into a diffuse thrust zone which eventually connects itself to the plate interplate beneath the southernmost ryukyu arc earthquakes beneath taiwan and the taiwan strait are compatible with the overall collision convergence and lateral extrusion nonetheless the lateral extrusion is probably only a secondary feature in the tectonic processes of taiwan

outer rise earthquakes and seismic coupling outer rise earthquakes can be separated into two basic categories based on focal mechanism compressional and tensional events reflecting compressional and tensional outer rise stress regimes we propose a model that connects the stress regime in the outer rise to the behavior of large subduction earthquakes and explains the focal mechanisms and both the temporal and spatial occurrence of outer rise earthquakes we find that compressional outer rise events tend to occur oceanward of seismic gaps in seismically coupled subduction zones or equivalently just prior to large subduction zone earthquakes the compressional stress regime is probably caused by slip in adjacent subduction zone segments there are two modes of occurrence for tensional outer rise earthquakes depending on whether they occur in coupled or uncoupled subduction zones in coupled subduction zones tensional events tend to occur after large subduction events as tensional stress from slab pull is temporarily transmitted to the outer rise in subduction zones that are inherently uncoupled tensional events can occur at any time therefore in uncoupled subduction zones the outer rise is always in tension and only tensional outer rise events occur there while in coupled subduction zones the outer rise is in a compressional state prior to large subduction earthquakes and a tensional state afterwards due to the large displacements

taiwan orogeny thin skinned or lithospheric collision abstract the taiwan orogeny is young and presently very active it provides an excellent environment for studying ongoing orogenic processes especially since the region is monitored intensively with dense seismological and geodetic networks and new studies aiming at deciphering shallow and deep structures in and around taiwan have been recently conducted or are being planned the available data can be used continually to test critically hypotheses of the taiwan orogeny hypotheses dealing with the mechanics of mountain building are basic to the understanding of taiwan orogeny and are particularly amenable to testing the widely cited thin skinned tectonics hypothesis was formulated to explain mainly the geologic and relatively shallow

focal mechanism determination in taiwan by genetic algorithm we determined the focal mechanism solutions for earthquakes with mag nitude ml 4 0 that occurred in the taiwan region between 1991 and 2005 first motion polarities of p waves recorded at over 700 seismic stations in taiwan were used because of the large number of events and stations involved we implemented the genetic algorithm in a nonlinear global search for the focal mechanism solutions the algorithm was tuned and validated through synthetic tests we finally determined the focal mechanisms of 1635 events with good qualities among 4188 earthquakes focal mechanism solutions for a majority of the earthquakes display a dominant pat tern of thrust fault type reflecting the compressive stress field due to the plate colli sion normal fault events occurred at intermediate depths in subduction zones which is likely the result of the bending of the subducting slabs strike slip faults are also found within the eurasia plate around the peikang basement high and in collision zones near ilan where the geometry of the colliding plates is complex our study pro vides a database of focal mechanisms for studying seismogenic structures and plate tectonics this database can also be used by structural seismologists to compute syn thetics for waveform tomography studies

oblique plate convergence slip vectors and forearc deformation slip vectors from thrust earthquakes at subduction zones where convergence is oblique to the trench often point between the directions of relative plate convergence and normal to the trench axis suggesting that oblique convergence is taken up by partial decoupling decoupling means that a component of arc parallel motion of the leading edge of the upper plate results in less oblique thrusting at the trench partial decoupling is modeled by partitioning of oblique convergence into slip on thrust and strike slip faults that are parallel to the trench and to each other and starting with a force equilibrium condition a relationship between the obliquity and the earthquake slip vector orientation is derived assuming that either fault slips when shear stress on it reaches a yield stress oblique slip parallel to the plate vector should occur on the thrust fault when obliquity is smaller than a critical angle for obliquity at or greater than this angle the stress on the strike slip fault is large enough to start it slipping and when both faults are active the arc parallel motion of the forearc deflects the slip vector back toward the trench normal if we assume that continued slip on either fault occurs at constant stress but the two faults can be at different stresses the slip vector will maintain a constant angle relative to the trench normal even when obliquity is larger than the critical angle this limiting angle of the slip vector called max measured relative to the trench normal is simply the arcsine of the ratio of the shear forces resisting slip on the strike slip and thrust faults a consequence is that when the obliquity exceeds max the slip vectors on the thrust fault are sensitive only to the thrust fault orientation and contain no information about the convergence direction between the plates slip vectors at the java trench southwest of sumatra show the relationship clearly with max 20 5 while slip vectors at the aleutian trench show the relationship less clearly with max 25 to 45 the greater angle at the aleutian trench suggests that the upper plate is stronger in the aleutian arc relative to the thrust fault than in the sumatran arc consistent with the sumatran arc being continental and having a well developed strike slip fault while the aleutian arc is oceanic and without a clear transcurrent fault slip vectors at the philippine trench which like sumatra has a large strike slip fault inboard of it tend to stay within 25 of the trench normal when obliquity is as large as 50 if obliquity exceeds max and continues to increase along a subduction zone the rate of motion of the forearc relative to the upper plate will vary with obliquity in which case the forearc sliver should extend or contract parallel to the arc from the geometry of modern island arcs arc parallel extension should be the more common and has been hypothesized for both sumatra and the aleutians on the basis of earthquake slip vectors and for these and other arcs from geological observations from estimates of max and the arc parallel gradients in obliquity arc parallel strain rates are estimated to be 1 to 3 10 8 yr for the sumatran forearc 2 to 6 10 8 yr for the aleutian forearc and 0 3 to 3 10 8 yr for the philippine forearc oblique convergence and subsequent arc parallel extension if accompanied by crustal thinning may provide an important yet little appreciated mechanism for bringing high grade metamorphic rocks to the surface of subduction complexes

evidence for early cretaceous oceanic crust trapped in the philippine sea plate abstract the huatung basin is a small oceanic basin located east of taiwan previous age estimates from magnetic lineation studies indicated an eocene age for the basin and formation from the central basin spreading center of the west philippine basin during the last phase of spreading in middle eocene new ar ar ages obtained on gabbros dredged on oceanic basement highs of the huatung basin are early cretaceous these old ages are consistent with early cretaceous ages determined on radiolarian assemblages from lanyu island luzon arc we have performed magnetic anomalies modeling for an early cretaceous oceanic crust our results are in good agreement with new ar ar ages determinations the best fit is indeed obtained with an opening of the huatung basin during the early cretaceous from 131 to 119 ma with a half spreading rate varying between 25 and 30 mm yr the spreading center appears to be located south of the actual basin the abnormal depth 5500 m instead of 5900 m and thickness 12 km instead of 6 km of the crust beneath the basin indicate that there was probably an excess supply of magma during its formation we propose that the basin is a fragment of the former proto south china sea or possibly the new guinea basin that has been trapped by the philippine sea plate

thirty year land elevation change from subsidence to uplift following the termination of groundwater pumping and its geological implications in the metropolitan taipei basin northern taiwan several levelling routes in the metropolitan taipei basin have been repeatedly conducted during the past decades mainly in order to monitor the anthropogenic ground subsidence due to massive pumping of groundwater we analysed the levelling data released from government and investigated the rate of ground level change from 1975 to 2003 which postdate the massive groundwater exploitation in taipei area based on the contour maps created from the levelling data of 406 benchmarks the overall subsidence rate in the taipei basin gradually decreased since 1975 and around 1989 the basin switched to slight uplift throughout a large part of the basin three mechanisms are proposed to be responsible for the observed land elevation changes including shallow soil compaction deformation within aquifer and tectonic subsidence the trend of the ground level change in 1975 2003 essentially demonstrated the effects of natural recharge to previously depleted aquifers and is explained by the hydro mechanical coupling of aquifer materials i e elastic rebound to the rising piezometric level the rate of shallow soil compaction is estimated about 1 8 mm yr throughout the basin according primarily to the shallow clay thickness asymmetric tectonic subsidence related to the shanchiao fault was estimated to be 1 75 mm yr and 0 9 mm yr in the western part and the central part of the basin respectively by subtracting the components of the soil compaction and tectonic subsidence from the surface land elevation change the rebound of aquifer strata was estimated to be about 6 7 cm and 16 cm in western margin and central taipei respectively the amount of rebound is approximately 10 in magnitude comparing to the amount of previous anthropogenic subsidence due to massive groundwater pumping totally about 2 m 2007 elsevier b v all rights reserved

plate tearing in the northwestern corner of the subducting philippine sea plate the philippine sea plate phs simultaneously subducts northwestward and collides eastward with the eurasian plate eu in northeast taiwan these two tectonic events induce high seismic activity which makes northeastern taiwan one of the most seismically active zones in the world to understand the mechanical processes at work we used existing geophysical data and the aftershocks recorded following a recent large strike slip event occurring within the phs oceanic crust during this event a nw se trending left lateral sub parallel to the phs eu convergence vector was active as a consequence of the collision subduction plate geometry we show that the lithosphere of the northwestern corner of the phs has been torn in a nw se orientation this tectonic feature is associated with an abrupt tectonic stress boundary and could generate large intra plate earthquakes

cenozoic plate reconstruction of the south china sea region abstract reconstructions of the south china sea region at 60 ma 40 ma 30 ma 20 ma 10 ma and 5 ma are presented we have attempted to place the south china sea basin in a regional tectonic framework the tectonic evolution of the major blocks surrounding the south china sea were analyzed as well as the relative motions of the indian and australian plates we have tried to correct the tectonic models available in this region a 3 d graphics terminal was used to derive rotation poles for the different tectonic blocks and our model was then tested to determine its self consistency when the model conflicted with previous interpretations the input data were evaluated for alternative explanations at least two and possibly three stages of extension can be recognized in this region the earliest one active in the late cretaceous to eocene involved nw se extension the second one active from the late eocene to early miocene involved north south extension the third stage of extension which probably trended nw se can be dated as post oligocene the first extensional event produced the ne sw trending proto south china sea and a series of sedimentary basins along the south china margin following the southeastward extrusion of indochina the proto south china sea was mostly consumed at the palawan trough renewed north south extension in the south china continental margin started the present day south china sea spreading in the oligocene the southeastward extrusion of indochina blocked by sundaland resulted in the nw se trending opening of the south china sea basin in the early miocene collision of the north palawan microcontinental block with the west philippines block stopped the opening of the south china sea at the end of early miocene spreading activity switched to the sulu sea basin in the middle miocene but collision between the sulu ridge and the west philippines at mindanao halted the opening of the sulu sea at the end of the middle miocene in the late miocene greater india continued its northward path and seems to have ripped open the andaman sea in the pliocene subduction along the northern manila trench placed the north luzon arc on a collision path with the east asia continental margin at taiwan our reconstructions along with detailed geological and geophysical information may be used as a predictive tool for basin evolution models and block interactions in this region the development of the south china sea basin the gulf of thailand the malay basin and the central thailand basins are the result of collision induced extensional forces the sulu celebes and sumatra basins were formed as a consequence of prolonged subduction the opening of the pearl river mouth west natuna south china sea sulu and possibly celebes basins were terminated by various plate collisions during the course of plate reorganizations major boundary faults have changed their slip senses during different stages of evolution

terrane motion by strike slip faulting of forearc slivers forearc slivers bounded by a trench and an active strike slip fault occur in about 50 of modern subduction zones analysis of earthquake slip vectors indicates that modern sliver terranes typically migrate at rates of 1 2 cm yr tertiary transport of some forearc slivers by 1000 km or more is therefore expected as a consequence of normal subduction active arc parallel strike slip faulting occurs whenever interplate coupling is strong and convergence is somewhat oblique strike slip rate increases with greater convergence obliquity

gravity recovery and climate experiment grace alias error from ocean tides 1 gravity recovery and climate experiment grace is able to observe interannual and longer period mass redistribution driven by climate changes however alias errors due to inaccurate ocean tide models may contaminate these signals to quantify this problem we estimate alias contamination due to expected errors in eight tidal constituents q1 o1 p1 k1 n2 m2 s2 and k2 for the period august 2002 to february 2006 differences between got00 2 and tpxo6 2 tide models are used as estimates of tide model errors the alias estimate is based on an approximate method employing least squares fits of spherical harmonics to potential differences between the two grace satellites results from this method are similar to previous estimates confirming this is a useful and efficient way to study the alias problem new findings are 1 the period of the m2 tide alias is about 140 d much longer than the predicted period 13 5 d because of monthly sampling of the 13 5 d period 2 there are two alias periods for k1 one about 90 d and one longer than 7 years and 3 the slow decay in the grace orbit causes tide component aliases to differ from pure frequencies as the orbit evolves we also find that aliases can contaminate in both space and time and thus may affect estimates of mass changes over land finally we compare the simulated s2 alias with the estimated s2 alias in grace data they are similar but our approximate method appears to overestimate alias error over high latitude oceans and along certain coastlines while underestimating error over low latitude oceans

paleoseismology of the rueisuei segment of the longitudinal valley fault eastern taiwan paleoseismic trench excavations in the dafu area of eastern taiwan have provided data on the rupture history of the rueisuei segment of the longitudinal valley fault during the late holocene the 1951 taitung earthquake ruptured the dafu site which is characterized by several terrace raises with late holocene sediments uplifted by the westward thrust fault trenches across the northwest facing fault scarp exposed fluvial and alluvial deposit sediments although nearly all of the 23 radiocarbon ages vary somewhat within each layer the overall age determined for each layer is in good accord with the stratigraphic ordering of these layers based on stratigraphic ordering and a statistical comparison of radiocarbon dates using the oxcal program we estimate that two pre 1951 earthquake surface ruptures at the dafu site occurred in the periods a d 1736 1898 and a d 1564 1680 the same oxcal model constrains the past two recurrence intervals to about 165 and 140 yr although with sizable uncertainties 55 215 and 90 260 yr respectively which are 95 ranges through the correlation of three trenches across the longitudinal valley fault we are able to identify evidence for at least three earthquakes with moment magnitudes of about 7 0 7 2 that occurred up to 390 yr prior to and during 1951 furthermore based on the radiocarbon dates the mean recurrence interval is roughly 150 yr uncertainty is indeterminate with a minimum vertical uplift rate of 8 5 12 2 mm yr

seismicity in the outer rise forearc region and configuration of the subducting lithosphere with special reference to the japan trench along nine subduction zones where the largest mw or ms 7 6 normal faulting earthquakes have occurred beneath the outer rise there is an apparent correlation between a large curvature of the subducting slab and regions of low background seismicity in the forearc in six cases a gentle inner trench is also observed taking advantage of the excellent coverage of teleseismic stations and extensive previous studies using local data we investigate the nature of this spatial correlation between seismicity and bathymetry along the honshu forearc in some detail an aseismic region in the forearc lies directly west of the rupture zone of the 1933 sanriku earthquake mw 8 4 a great normal faulting earthquake that most likely has ruptured the entire brittle portion of the lithosphere in the outer rise along this particular cross section of the japan trench large to moderate sized interplate earthquakes span a short interval of depths between approximately 14 and 31 km about half the interval observed along an adjacent profile to the south 12 45 km the aseismic portion of the plate interface between depths of approximately 30 and 45 km has persisted for over 100 years the occurrence of lithospheric normal faulting earthquakes in the outer rise region indicates that bending moment within the subducting slab limited by yield strength of the lithosphere has nearly reached an asymptotic value moment saturation at this state the curvature of the plate can increase with little increase in bending moment as such any stresses transmitted across the plate interface will be relaxed by a permanent increase in the plate s curvature and little elastic strain can accumulate across the interplate thrust zone to produce large earthquakes thus the saturation of bending moment in the subducting lithosphere offers a straightforward interpretation for the observed correlation between seismicity in the outer rise forearc region and configuration of the subducting lithosphere

too low magnetic inclinations in central asia an indication of a long term tertiary non dipole field observed tertiary palaeolatitudes in central asia are more southerly by about 1600 km on average than those predicted from the eurasian reference palaeopoles subsequent northward displacements of the central asian terranes are unlikely to have been this large in this study we analyse to what extent non dipole fields especially octupole fields can explain this phenomenon a global zonal octupole field manifests itself in two ways 1 because the reference apwp is based mostly on results from the uk and north america its palaeopoles will be far sided as seen from the north atlantic but near sided as seen from eastern asia giving predicted palaeolatitudes that are too high 2 an octupole field contribution produces observed palaeolatitudes as calculated with the dipole formula for central asia that are too low both effects therefore increase the palaeolatitude anomalies in asia we find that an octupole dipole field ratio g3 of 0 06 or greater will reduce the palaeolatitude discrepancies significantly and is of the same magnitude as the g3 estimate of a recent analysis of early tertiary european and north american data

spatial variation of seismogenic depths of crustal earthquakes in the taiwan region implications for seismic hazard assessment abstract this paper presents the first whole taiwan scale spatial variation of the seismogenic zone using a high quality crustal seismicity catalog the seismicity onset and cutoff depths i e seismogenic depths are determined by the earthquake depth moment distribution and used to define the upper and lower boundaries of the seismogenic zone respectively together with the published fault geometries and fault area moment magnitude relations the depth difference in the onset and cutoff depths i e seismogenic thickness is used as the fault width to determine the moment magnitudes of potential earthquakes for the major seismogenic faults results show that the largest mw7 9 8 0 potential earthquake may occur along the changhua fault in western taiwan where the seismic risk is relatively high and seismic hazard mitigation should be a matter of urgent concern in addition the first motion focal mechanism catalog is used to examine the relation between the seismogenic depths and earthquake source parameters for crustal earthquakes 50 km the shallowest onset and cutoff depths are observed for normal and strike slip events respectively this observation is different from the prediction of the conventional continental rheology model which states that thrust events have the shallowest cutoff depth thus a more sophisticated rheology model is necessary to explain our observed dependence of the seismogenic depths on faulting types meanwhile for intermediate to large crustal mw 4 depth 50 km earthquakes thrust events tend to occur at the bottom region of the seismogenic zone but normal and strike slip events distribute at a large depth range

sequential thrusting model for mountain building constraints from geology and heat flow of taiwan in this study we analyze the thermal evolution of the taiwan fold and thrust belt based on the available geological information a finite element program is used to simulate the complex kinematics of the observed sequential thrusting and the interplay of uplift erosion and sedimentation processes the discrete sequential thrusting model successfully reproduces the observed decrease in heat flow at the deformation front in northern taiwan this reveals that the heat flow drop may be due to the recent stacking of cold thrust sheets at the deformation front this inference is consistent with the high seismicity and dense quaternary fault distribution in western taiwan the model also successfully reproduces the observed high heat flow in the central ranges suggesting that the observed high heat flow may be a result of tectonic unroofing due to active underplating uplift and erosion these conclusions are insensitive to uncertainties in the frictional coefficient on the thrust plane

collisional tectonics between the eurasian and philippine sea plates from tomography evidences in southeast china abstract the upper mantle structure of southeast china is important for us to understand the deformation and mantle dynamics process associated with the interaction between the eurasian plate and philippine sea phs slab we determined a detailed three dimensional p wave velocity vp structure of the crust and upper mantle down to 400 km depth beneath southeast china by applying teleseismic tomography to 6869 high quality p wave arrival times the data were collected very carefully from the original seismograms of 635 teleseismic events recorded by 65 broadband stations deployed in southeast china our images show that the high vp phs slab subducts toward the north along the ryukyu trench at the latitude of about 24 n and extends down to 350 km depth and even more high vp anomalies are imaged in the upper mantle under central and southern taiwan which represent the subducted eurasian plate break off eurasian plate at a big angle subducting eastward is revealed under central taiwan at depths from the upper mantle to 400 km while continuous eurasian plate under south taiwan is mainly imaged from the moho down to 400 km depth a torn mantle window within the eurasian continent beneath central and northern taiwan created by the northward motion of the philippine sea plate is the upwelling path of the asthenosphere the tomographic images also show the low vp anomalies spread widely under the coastal areas of mainland china and taiwan strait the structure of the crust and upper mantle suggests that the mountain building process in the central part of taiwan is mainly attributed to the subduction collision tectonics at the boundary between the eurasian continental lithosphere and the subducting oceanic lithosphere of the phs slab

extension during continental convergence with application to the tibetan plateau the tibetan plateau is the product of crustal thickening caused by the collision between india and asia and is the largest active example of extensional tectonics in a zone of continental collision throughout most of the tertiary the tectonics of the plateau were dominated by north south shortening a significant proportion of which took place on east west striking thrust faults for the last 5 m y or so the plateau has been thinning by the mechanism of extension on north south trending normal faults previous investigations of the collision have been able to account for the large scale features of the tertiary deformation but have failed to explain the transition in the late tertiary to quaternary strain rate field of the plateau from north south compression to east west extension the transition could in principle be effected by a reduction in the rate of convergence between india and asia or by a uniform reduction in strength of the whole continental lithosphere of asia it could not however be effected by a reduction in strength of the elevated region alone this produces increased compressional strain rates in the weakened zone an alternative explanation for the transition to extension comes from considering the thermal evolution of thickened continental lithosphere the lower part of the lithosphere consists of a thermal boundary layer which when thickened by horizontal shortening is colder and denser than its surroundings convective instability of the thickened thermal boundary layer and its replacement by hot asthenosphere would rapidly raise the surface elevation and gravitational potential energy of the overlying part of the lithosphere the convective instability would happen in a time brief compared with the collision time scale 50 m y in the case of india and asia but would only occur after there had already been substantial thickening of the lithosphere numerical experiments show that for a range of lithospheric parameters the increase in surface height as much as 2 km and of potential energy 5 to 10 1012 n m 1 resulting from convective instability of the lower lithosphere are sufficient for east west extension to replace north south compression as the dominant feature of the stress field of the tibetan plateau

determinations of focal depths of earthquakes associated with the bending of oceanic plates at trenches abstract the earthquakes examined in this paper are all within the oceanic lithosphere and are associated with the bending of plates before subduction accurate determinations of the depth of these earthquakes are needed to study the stress pattern within a bending plate routinely determined depths of shallow sub oceanic earthquakes published in bulletins are unreliable the depths can be accurately determined to within a few kilometers if the original seismograms from these events are studied in some cases the reflected phases pp and pwp can be clearly identified there exists the possibility that the wave reflected at the water air interface pwp may be misidentified as pp leading to erroneous estimates of depth additional methods of analysis such as surface wave radiation patterns or the apparent frequency dependence of reflection at the crust water interface can remove this possible source of confusion one of the most powerful techniques for depth analysis is the modelling of long period waveforms the pattern of stresses within the bending oceanic lithosphere revealed by the depths and focal mechanisms of these intraplate earthquakes is one of horizontal deviatoric tension down to a depth of about 25 km with horizontal compression at greater depths

back arc rifting in the okinawa trough abstract geological and geophysical data reveal that the okinawa trough shows incipient continental rifting and crustal separation started from about 2 ma the early extensional movements in the trough are probably of miocene age in addition to the miocene phase two main periods of extension are recognized a pleistocene phase between 1 9 and 0 5 ma and the present day phase during the stage short central rifts central grabens were formed the opening however may have occurred only in the southern part of the trough basin having an average half spreading rate of 2 cm yr 1 since early pleistocene time producing its present width of several tens of kilometres these activities were well represented by igneous intrusions sedimentary facies and sedimentary structures in and around the okinawa trough the width of the zone affected by back arc extension defined as greater okinawa trough is larger than the present okinawa trough whose width is 200 250 km the present form of the greater okinawa trough started to form at the same time as that of the okinawa trough

back arc extension in the okinawa trough the okinawa trough lying to the east of china is a back arc basin formed by extension within continental lithosphere behind the ryukyu trench arc system middle to late miocene uplift associated with normal faulting of the initially adjacent ryukyu nonvolcanic arc and the taiwan sinzi folded belt corresponds to the first rifting phase the timing of rifting is supported by the presence of marine sediments of corresponding age drilled in the northern okinawa trough the rifting occurred after a major early miocene change in the motion of the philippine plate with respect to eurasia and ceased during the pliocene a second rifting phase started about 2 m y ago at the plio pleistocene boundary and has continued until the present time it has proceeded to a more advanced stage in the middle and southern okinawa trough than it has farther north detailed bathymetric sea beam seismic reflection and magnetics data collected during the pop 1 cruise of the r v jean charcot reveal the principal features of the extensional processes the back arc spreading phase started very recently in the southern and middle okinawa trough as exemplified by several en echelon and in some cases overlapping active central graben oriented n70 e n80 e some of these depressions are intruded by volcanic ridges of fresh back arc basalt with associated large magnetic anomalies transform faults between these en echelon active rifts are not obvious we suggest that the major part of the southern okinawa trough is underlain by a thinned continental crust and that except for the system of en echelon rifts of the southern okinawa trough the back arc basin oceanic domain is limited to a width of a few tens of kilometers or less in the axial portion of the trough the system of axial back arc volcanic ridges that occur in the rifts ends at the latitude of okinawa island whereas active volcanoes in the ryukyu arc occur only north of okinawa island we refer to this transition between active arc and back arc volcanism as the volcanic arc rift migration phenomenon vamp globally back arc volcanism propagated from the southern okinawa trough to the okinawa vamp area rifting continues to occur in the northern okinawa trough but is not yet accompanied by associated volcanism the okinawa vamp area is characterized by a series of parallel basaltic ridges oriented n75 e with associated linear magnetic anomalies characteristic of dyke intrusions we suggest that the formation of the back arc oceanic domain took place along the axial back arc extensional zone trending n75 e and that this zone presently ends at the southern extremity of the active volcanic chain the initial phase of formation of back arc basin oceanic crust is non steady state and is characterized by the lack of a developed fracture zone pattern the termination of the vamp area in the direction of the volcanic zone of the arc is consistent with the suggestion of molnar and atwater that the volcanic arc is a fundamental line of weakness which determines where initial back arc spreading occurs apparently back arc extension initially occurred within the continental lithosphere located westward of the ryukyu arc along its whole length but the subsequent back arc volcanism was initiated in the southernmost portion of the region and then moved northward this migration was accompanied by the shutting down of volcanic activity along the abandoned portions of the arc it is this transfer of volcanism that we call the vamp process thus arc and back arc basin volcanism seem to be associated in such a manner that spreading tends to migrate simultaneously with a cessation in volcanic activity along the arc this interplay of arc and back arc activity is probably linked to changes in the parameters of plate convergence since the plate motion in the phillippine sea is oblique to the trench at least in the southern part of the okinawa trough we suggest that the oblique resisting force applied to the edge of the overriding plate engenders the development of en echelon extensional features within and behind the arc the motion of the ryukyu platelet with respect to eurasia is consequently an extensional strike slip motion trending roughly parallel to the okinawa trough

microseismicity in the middle okinawa trough we conducted a microearthquake survey with ocean bottom seismometers obs in the middle okinawa trough in 1984 during about one month of operation very active microseismicity was observed the microearthquake activity was characterized by the ceaseless occurrence of numerous small earthquakes in the vicinity of obs stations this activity was recorded by almost every obs station but not by land based instruments swarm activities also occurred at many places average activity levels of earthquakes whose durations are longer than 30 s recorded at individual obss are from 4 to 20 events per day which seem higher than those of ordinary mid oceanic ridge fracture zone systems and other back arc basins apparently the number of smaller size earthquakes is fewer than expected this might indicate the rather high attenuation of seismic waves or low q of the crust around obs stations which makes it difficult to detect small size events at a distant range attenuation of seismic waves was large beyond the epicentral distance of about 120 km this seems to indicate low q in the lower crust also a particularly low q area was found at the northern most part of the research area where high heat flow measurements were obtained this suggests that igneous activity reaches to the surface in that vicinity from the bottom of the crust where magmatic activity has been implied in other studies at least two types of swarm activities are indicated for this area some have no dominant events show high b value properties and seem to be volcanic others have dominant events and appear to be tectonic in nature among the swarm activity two kinds of signals closely resembling volcanic tremors were observed from these features we conclude that the crust of this area is highly fractured and subject to intense tectonic activities as well as some kind of igneous processes

global seismicity of 1981 centroid moment tensor solutions for 542 earthquakes abstract centroid moment tensor solutions are presented for 542 earthquakes that occurred during 1981

earthquakes along the ryukyu kyushu arc strain segmentation lateral compression and the thermomechanical state of the plate interface we systematically determined the focal depths and mechanisms of 49 large to moderate sized earthquakes mb 5 4 that occurred along the ryukyu kyushu arc since 1963 by inverting the waveforms and amplitudes of p and sh wave trains at teleseismic distances the results are sufficiently precise to delineate seismogenic structures near the plate interface in detail revealing features not predicted by plate kinematics in contrast to previous studies of this arc and those along other subduction zones shallow seismicity along the plate interface is systematic showing two distinct layers of activities the second layer of seismicity is delineated by a few earthquakes that occurred at depths between 50 and 65 km some 10 20 km directly beneath the seismogenic portion of the interplate thrust zone these earthquakes indicate lateral compressional strain within the subducted slab as their p axes are subparallel to the local strike of the arc not perpendicular to the arc as one would expect from the direction of plate convergence the occurrence of these events cannot be accounted for by membrane stress due to the geometry of the subducted slab to our knowledge similar earthquakes occurred only beneath the northern indoburman ranges and along the banda arc where subducted slabs as part of the indian plate are being dragged northward with their northern edges bumping into east west trending collision zones nearby by drawing an analogy between the tectonic settings of these three regions we interpret events beneath the plate interface along the ryukyu kyushu arc as a consequence of ongoing collision between the philippine sea plate and eurasia near taiwan the interplate thrust zone is largely aseismic down to a depth of approximately 30 km a large number of earthquakes showing low angle thrust faulting commence at this depth and are accompanied by two events that show antithetic thrust faulting at a slightly shallower depth of 20 25 km these observations suggest that the strength of the plate interface increases significantly below the depth of proximately 25 5 km based upon available heat flow measurements between the trench axis and the volcanic arc we estimated the temperature field and magnitude of shear traction along the interplate thrust zone the deepest interplate earthquakes 40 50 km correspond to temperatures of approximately 730 980 c comparable to the limiting temperature of intraplate mantle earthquakes 800 c and that of intermediate and deep focus earthquakes at wadati benioff zones potential temperature 900 k therefore the cessation of seismicity in general is probably controlled by similar temperature conditions regardless of the tectonic settings of source regions the magnitude of shear traction along the seismogenic portion 40 50 km depth of the interplate thrust zone seems to be of the order of 100 mpa which in turn implies a low average coefficient of friction of only 0 10 0 05 considerably less than those observed for laboratory specimens the inception of intermediate focus earthquakes within the subducted slab is at a depth of about 100 km whereas shallow earthquakes show no apparent variations along the entire arc intermediate depth earthquakes indicate downdip extension along the northern end of the arc near kyushu but abruptly change to downdip compression along the rest of the arc all available evidence indicates that this is a spatial pattern unrelated to earthquake cycles at the moment tectonic interpretation of such a distinct strain segmentation along the arc is enigmatic without the presence of clear discontinuities in the age of the subducted slab the sudden switch in the strain field of the slab is difficult to explain by the continuously varying rate of subduction along the arc in any case this pattern of strain segmentation requires the presence of a major accommodation structure such as a tear fault in the subducted slab near the tokara channel appendix is available with entire article on microfiche order from american geophysical union 2000 florida ave n w washington dc 20009 document b91 004 2 50 payment must accompany order

variations of b values at the western edge of the ryukyu subduction zone north east taiwan using earthquakes relocated in north east taiwan we estimated b value distribution along a cross section located near the ryukyu slab edge and four b value anomalous areas are evidenced 1 a high b value body lying on top of a low vp low vs and high vp vs sausage like body was considered as a region of enhanced partial melt or water supply above which seismicity occurs 2 beneath the ilan plain an anomalous area characterized by b values slightly higher than 1 1 might give evidence to the magma conduits to the kueishantao island 3 above the ryukyu wadati benioff zone at depths ranging from 90 to 110 km a high b value anomaly might correspond to the depth where dehydration occurs in the subducting oceanic plate and 4 a low b value area located within the ryukyu slab at depths ranging from 70 to 90 km might be linked to the compressive mechanisms shown by focal mechanisms and the bending of the subducting plate

seismic coupling and outer rise earthquakes variation of interplate seismic coupling at subduction zones is a major factor controlling the size of the largest underthrusting events this variation also has a profound effect on the regional intraplate stresses in the vicinity of the subduction zone outer rise seismicity is strongly correlated with variations in interplate coupling reflecting the stress state of the interplate coupled zone over 200 outer rise earthquakes with known focal mechanisms are used to investigate the relationship between stresses in the outer rise and interplate seismic coupling these events occur within the downgoing i e oceanic plate near the bathymetric trench axes and generally fall into the categories of tensional normal or compressional thrust with their tensional or compressional stress axes oriented approximately horizontal and perpendicular to the trench in uncoupled subduction zones only tensional outer rise earthquakes occur which indicates that the outer rise is dominated by tensional stresses associated with plate bending and or slab pull forces in strongly coupled subduction zones both tensional and compressional outer rise events are found these events are related both spatially and temporally to the distribution of large underthrusting earthquakes and are thus an integral part of the earthquake cycle in the strongly coupled regions tensional outer rise events follow large underthrusting events as the outer rise is temporarily in tension due to the underthrusting motion compressional outer rise events take place as compressional stress slowly accumulates oceanward of locked sections of the interplate zone in four instances compressional outer rise earthquakes have been followed by large underthrusting events which have occurred 2 4 7 and 19 years after the associated outer rise event the remaining compressional outer rise events are located in regions that are either known seismic gaps or in regions where the seismic potential is unknown the occurrence of compressional outer rise earthquakes suggests that compressional stress is accumulating in the adjacent interplate region and that there is the potential for a future large underthrusting event in the region thirty compressional outer rise events have been located in trench segments of middle and south america the kurile islands the tonga and kermadec islands the new hebrides arc and the solomon islands regions in both the southern kamchatka and northern new hebrides regions the outer rise seismicity indicates that the stress regimes in the outer rise have changed with time from tensional following a previous large underthrusting event to compressional at present thus three stages of the cycle from underthrusting to tensional outer rise regime to compressional outer rise regime are present requiring only the occurrence of the next underthrusting event to complete the cycle the occurrence of compressional outer rise events is useful for assessing the seismic potential of a region on an intermediate time scale

shear wave anisotropy beneath the taiwan orogen abstract shear wave splitting in regional and teleseismic shear waves scs and s has been analyzed for fast polarization directions and delay times at 33 short period and 17 broad band permanent and temporary seismic stations in the taiwan region the orientations of anisotropy for the taiwan stations are generally parallel to the strike of the mountain belt with a clockwise rotation following the trend of regional geologic fabric in north taiwan the delay times on the other hand vary greatly from 0 5 to 2 1 s the largest values are from the southwestern taiwan foothills and an abrupt change in delay times is found in south taiwan the main source of shear wave splitting in taiwan is compatible with an anisotropic region between 25 and 230 km depths the observed splitting parameters are consistent with upper mantle anisotropy resulting from the collisional tectonics that built the taiwan orogen

transition from oblique subduction to collision earthquakes in the southernmost ryukyu arc taiwan region tectonic characteristics of the region between taiwan and the southernmost ryukyu arc are inferred from a detailed analysis of local seismicity and source parameters of 62 recent earthquakes of 5 5 mb 6 6 five major seismogenic structures can be delineated the collision seismic zone csz the interface seismic zone isz the wadati benioff seismic zone wbsz the lateral compression seismic zone lcsz and the okinawa seismic zone osz in the csz located along the east coast of taiwan and offshore earthquake focal mechanisms show horizontal p axes distributed in two directions 287 10 and 333 16 possibly reflecting a strain partition associated with the relative plate convergence between the eurasia plate and the philippine sea plate the corresponding seismic strain tensor indicates a maximum compressive strain rate of 1 2 10 7 yr 1 along 293 and a comparable extension in vertical direction presumably resulted from plate collision in the region the geometry of the isz which is distorted significantly at its westernmost end can be approximated by a north dipping plane that is gradually pushed northward with increasing dip the seismogenic portion of the interface spans a short depth range from 10 km to 35 km a clear pattern of earthquake slip partition is observed the average slip vector residual is as large as 35 seismic strain patterns within the subducted philippine sea slab show predominantly downdip extension between 80 and 120 km and downdip compression at 270 km different from the pattern of strain segmentation observed for the rest of the ryukyu arc where the northern and southern portions are dominated by downdip extension and compression respectively owing to the large convergence obliquity the slab is descending at a rate significantly slower near taiwan than in the southern ryukyu thus we interpret the appearance of downdip extension within the subducted lithosphere as a combined result of oblique subduction and the slab s negative buoyancy a number of thrust or oblique strike slip earthquakes between isz and wbsz show a consistent pattern of lateral compression with p axes oriented roughly parallel to the local strike of the trench arc system they are probably due to the compressive strain originated from the collision and transmitted laterally within the lithosphere shallow normal faulting earthquakes show successive rotation of t axes from approximately n s in the okinawa trough to approximately e w in northeast taiwan possibly as a result of interaction between the extension from the opening of the okinawa trough and the compression from collision one normal event january 18 1991 mb 5 9 occurred in the central range with t axis roughly parallel to the structural trend of taiwan implying that the nature of the orogeny in taiwan has changed from thin skinned deformation to lithospheric collision involving the whole crust and uppermost mantle

trench parallel stretching and folding of forearc basins and lateral migration of the accretionary wedge in the southern ryukyus a case of strain partition caused by oblique convergence detailed seafloor mapping in the area east of taiwan revealed trench parallel stretching and folding of the ryukyu forearc and lateral motion of the accretionary wedge under oblique convergence east of 122 40 e a steep accretionary wedge is elongated in an e w direction a major transcurrent right lateral strike slip fault accommodates the strain partitioning caused by an oblique convergence of 40 a spectacular out of sequence thrust may be related to the subduction of a structural high lying in the axis of the n s trending gagua ridge this asperity is likely responsible for the uplift of the accretionary wedge and forearc basement and may have augmented strain partitioning by increasing the coupling between the two plates west of 122 40 e the low taper accretionary wedge is sheared in a direction subparallel to the convergence vector with respect to the ryukyu arc the bayonet shape of the southern ryukyu arc slope partly results from the recent re opening of the southern okinawa trough at a rate of about 2 to 4 cm yr right lateral shearing of the sedimentary forearc with respect to the nonlinear ryukyu backstop generates trench parallel extension in the forearc sediment sequence at dilational jogs and trench parallel folding at compressive jogs the hoping basin lies above a diffuse trench trench fault ttf or tff unstable triple junction moving toward the south along a n s transform zone which accommodates the southward drift of the ryukyu arc with respect to eurasia

origin and evolution of the west philippine basin a new interpretation abstract combined analyses of the magnetic lineations and seafloor structures of the west philippine basin show that it was formed by spreading from the central basin spreading center in two distinctly different spreading phases before and after subduction was initiated along the palau kyushu trend at 45 ma b p these episodes are recognized by differences in the strike of the magnetic lineations spreading rates density and strike of fracture zones and basement relief from 60 45 ma b p spreading was ne sw relative to present orientations at a half rate of 44 mm a at about 45 ma b p the spreading direction changed to a more n s direction with a reconfiguration of the central basin spreading center into numerous short e w segments offset by closely spaced n s transform faults spreading slowed to a half rate of 18 mm a and ceased at 35 ma b p thus the west philippine basin originated at 45 ma b p by the trapping of normal ocean crust west of the initial subduction along the palau kyushu trend the 45 35 ma b p period represents the dying phase of spreading on the central basin spreading center following isolation of the west philippine basin from the plate driving forces of the pacific