spatial and temporal distribution of slip for the 1999 chi chi taiwan earthquake we investigated the rupture process of the 1999 chi chi taiwan earth quake using high quality near source strong motion records broadband teleseismic displacement waveforms and well distributed global positioning system gps data the near source strong motion displacement waveforms recorded significant static offsets of up to 8 m the teleseismic displacement records show a significant pulse with duration of about 18 to 20 sec taking into account the surface displacements observed along the chelungpu fault we considered two fault geometries a single planar fault and a two segment fault with a northeast striking section near the north ern end using the finite fault model with variable slip vectors we derived two models of the temporal and spatial slip distribution of the earthquake the gps data provided good surface displacement constraints for the slip distribution determina tion the spatial slip distribution is generally consistent with field observations the results for the simple fault model show a large asperity located in the region about 25 to 55 km north of the hypocenter with maximum slip of about 15 m when we use the two segment model the asperity further extends to the region where the fault bends toward the northeast with a maximum slip of up to 20 m a large amount of right lateral slip beneath station tcu068 is necessary to explain its observed large west movement it implies a local converging slip at the corner where the fault bends to the northeast the slip amplitude near the hypocenter is about 3 to 6 m the seismic moments determined from the various data sets are within the range of 2 to 4 10 27 dyne cm most of the slip concentrated at shallow depths less than 10 km the total rupture duration is about 28 sec and the rupture velocity is 75 to 80 of the shear wave velocity the slip vector shows a clockwise rotation during the fault rupture the static stress drop of the large asperity region is comparable with the dynamic stress drop as observed directly from the slip velocity at the station near the large slip region

crustal deformation in taiwan results from finite source inversions of six mw 5 8 chi chi aftershocks 1 we inverted strong motion data for the finite source parameters of six large aftershocks of the 1999 chi chi taiwan earthquake to investigate seismogenic structure in taiwan for each event we derived a preferred model by testing different focal mechanisms hypocenters rupture velocities and dislocation risetimes as well as different combinations of stations in more than 1000 inversions we documented how the fits between the waveforms and the corresponding synthetics deteriorated as the hypocenter and focal mechanism deviate from those of the preferred model if the deviation in focal mechanisms and hypocenters is less than 20 and 5 km respectively we generally recovered 80 of the preferred model s synthetic waveform fit unlike the dislocation risetime the rupture velocity used in the inversion had a strong influence on the waveform fits in this study we also used the slip models to study fault geometry two of the aftershocks ruptured on the southern extension of the main shock slip area one strike slip aftershock nucleated within the basement but ruptured mainly within the overlying sedimentary strata suggesting that seismogenic deformation in the basement can influence shallow structures p axes of the derived models have azimuths consistent with current plate motion finally gps displacement derived from the six slip models can explain 80 of the postseismic deformation observed in the aftershock regions indicating that studies of postseismic deformation must take into account the cumulative effects of large shallow aftershocks

horizontal coseismic deformation of the 1999 chi chi earthquake measured from spot satellite images implications for the seismic cycle along the western foothills of central taiwan 1 the 1999 chi chi earthquake mw 7 6 broke a major thrust fault along the western foothills of the central range of taiwan we have measured the horizontal coseismic displacement field by correlating optical satellite images acquired before and after the earthquake these data reveal the fault trace and a clockwise rotation of surface displacements toward the north with much larger displacements and strain in the hanging wall at the surface coseismic slip increases from 5 6 m near the epicenter to 10 11 m to the north in the epicentral area we observe a left lateral strike slip zone trending n125 e and farther north a fault zone trending n s with a right lateral component the data were modeled using elastic dislocations the fault geometry consists of a shallow

structural geology of a classic thrust belt earthquake the 1999 chi chi earthquake taiwan mw 7 6 abstract we document the structural context of the 1999 chi chi earthquake m w 7 6 in western taiwan which is one of the best instrumented thrust belt earthquakes the main surface break and large slip 3 10 m is on two segments of the shallow otherwise aseismic bedding parallel chelungpu sanyi thrust system which shows nearly classic ramp flat geometry with shallow detachments 1 6 km in the pliocene chinshui shale and mio pliocene kueichulin tungkeng formations however rupture is complex involving at least six faults including a previously unknown deeper thrust 8 10 km on which the rupture began we compare the coseismic displacements with a new 3d map of the chelungpu sanyi system the displacements are spatially and temporally heterogeneous and well correlated with discrete geometric segments of the 3d shape of the fault system geodetic displacement vectors are statistically parallel to the nearest adjacent fault segment and are parallel to large scale oblique fault corrugations the displacement magnitudes are heterogeneous at several scales which requires in the long term other non chi chi events or significant aseismic deformation the chelungpu thrust has a total displacement of 14 km but the area of largest chi chi slip 10 m is on a newly propagated north chelungpu chinshui detachment 0 3 km total slip which shows abnormally smooth rupture dynamics

spatial and temporal evolution of stress and slip rate during the 2000 tokai slow earthquake 1 we investigate an ongoing silent thrust event in the tokai seismic gap along the suruga nankai trough central japan prior to the event continuous gps data from april 1996 to the end of 1999 show that this region displaced 2 cm yr to the northwest relative to the landward plate the gps time series show an abrupt change in rate in mid june 2000 that continues as of mid 2005 we model this transient deformation which we refer to as the tokai slow thrust slip event as caused by slip on the interface between the philippine sea and amurian plates the spatial and temporal distribution of slip rate is estimated with kalman filter based inversion methods our inversions reveal two slow subevents the first initiated in late june 2000 slightly before the miyake jima eruption the locus of slip then propagated southeast in the second half of 2000 with maximum slip rates of about 15 cm yr through 2001 a second locus of slip initiated to the northeast in early 2001 the depth of the slip zone is about 25 km which may correspond to the transition zone from a seismogenic to a freely sliding zone the cumulative moment magnitude of the slow slip event up to november 2002 is mw 6 8 we calculate shear stress changes on the plate interface from the slip histories stress change as a function of slip rate shows trajectories similar to that inferred for high speed ruptures however the maximum velocity is 8 orders of magnitude less than in normal earthquakes

conjugate thrust faulting associated with the 1999 chi chi taiwan earthquake sequence 1 the geometrical structure of the responsible faults of the 20 september 1999 chi chi taiwan earthquake ml 7 3 mw 7 6 and its aftershocks can be clearly depicted by well located hypocenters and focal mechanisms of large aftershocks the mainshock and two large aftershocks with ml 6 8 were characterized by thrust faulting along a n s striking fault plane dipping to the east the underground structure of the chelungpu fault which is probably merging with the decollement beneath the western foothills can be clearly associated with the seismicity pattern and the focal mechanisms of the three largest events a group of deeper aftershocks including two moderate events ml 6 3 and 6 0 respectively were located to the southeast of the mainshock along a fault plane dipping steeply to the west down to a depth of about 30 km our results suggest that the spatial pattern of the aftershocks in the southern part of the source area can be interpreted by a conjugate fault system this conjugate fault system is comprised of the gently east dipping chelungpu fault and a steeply west dipping deeper fault zone

high frequency sediment level oscillations in the swash zone sediment level oscillations with heights of about 6 cm and shore normal lengths of order 10 m have been measured in the swash zone of a high energy coarse sand beach crests of oscillations were shore parallel and continuous alongshore the oscillations were of such low steepness height to length ratio approximately 0 006 that they were difficult to detect visually the period of oscillation ranged between 6 and 15 min and decreased landward across the swash zone the sediment level oscillations were progressive landward with an average migration rate in the middle to upper swash zone of 0 8 m min 1 migration was caused mostly by erosion on the seaward flank of the crest of an oscillation during a period of net seaward sediment transport thus the observed migration was a form migration landward rather than a migration involving net landward sediment transport the observed sediment level oscillations were different than sand waves or other swash zone bedforms previously described

imaging of aseismic fault slip transients recorded by dense geodetic networks summary the large scale continuous gps networks that have emerged in the last decade have documented a wide range of aseismic deformation transients that resulted from physical processes such as aseismic fault slip and magma intrusion in particular a new class of m 7 slow earthquakes with durations ranging from days to months have been observed with gps arrays located above the downdip portion of subduction zone thrust interfaces interpretation of the displacement time series resulting from these events is not straightforward owing to the contaminating effects of multiple contributing signals such as fault slip local benchmark motion seasonal effects and reference frame errors we have developed a time dependent inversion algorithm based on the extended kalman filter which can separate the various signals and allow the space time evolution of these slow slip transients to be studied in detail we applied the inversion algorithm to the 1999 cascadia slow earthquake this event had two primary episodes of moment release separated by a two week period in which relatively little moment release occurred the cascadia event and other slow earthquakes share numerous similarities with both ordinary earthquakes and afterslip transients suggesting that they may represent fault slip under a conditionally stable regime

postseismic relaxation driven by brittle creep a possible mechanism to reconcile geodetic measurements and the decay rate of aftershocks application to the chi chi earthquake taiwan we evaluate the effect of coseismic stress changes on the fault slip at midcrustal depth assuming a velocity strengthening brittle creep rheology we show that this model can help reconcile the time evolution of afterslip as measured from geodesy with aftershocks decay we propose an analytical expression for slip of the brittle creeping fault zone bcfz that applies to any dynamic or static stress perturbation including shear stress and normal stress changes the model predicts an initial logarithmic increase of slip with time postseismic slip rate decays over a characteristic time t r a that does not depend on the amplitude of the stress perturbation and it asymptotically joins the long term creep imposed by interseismic stress buildup given that the seismicity rate might be considered proportional to the sliding velocity of the bcfz the model predicts a decay rate of aftershocks that follows omori s law with a mathematical formalism identical to that of dieterich 1994 although based on a different mechanical rationale our model also differs from dieterich s model in that it requires that aftershock sequences and deep afterslip as constrained from geodetic measurements should follow the same temporal evolution we test this for the 1999 chi chi earthquake m w 7 6 and find that both sets of data are consistent with a model of afterslip due to the response of the bcfz the inferred relaxation time t r 8 5 years requires a value for a log v being the coefficient of friction in the range between 1 3 10 3 and 10 2

imaging the ramp decollement geometry of the chelungpu fault using coseismic gps displacements from the 1999 chi chi taiwan earthquake abstract we use coseismic gps data from the 1999 chi chi taiwan earthquake to estimate the subsurface shape of the chelungpu fault that ruptured during the earthquake studies prior to the earthquake suggest a ramp decollement geometry for the chelungpu fault yet many finite source inversions using gps and seismic data assume slip occurred on the down dip extension of the chelungpu ramp rather than on a sub horizontal decollement we test whether slip occurred on the decollement or the down dip extension of the ramp using well established methods of inverting gps data for geometry and slip on faults represented as elastic dislocations we find that a significant portion of the coseismic slip did indeed occur on a sub horizontal decollement located at 8 km depth the slip on the decollement contributes 21 of the total modeled moment release we estimate the fault geometry assuming several different models for the distribution of elastic properties in the earth homogeneous layered and layered with lateral material contrast across the fault it is shown however that heterogeneity has little influence on our estimated fault geometry we also investigate several competing interpretations of deformation within the e w trending rupture zone at the northern end of the 1999 ground ruptures we demonstrate that the gps data require a 22 to 35 km long lateral ramp at the northern end contradicting other investigations that propose deformation is concentrated within 10 km of the chelungpu fault lastly we propose a simple tectonic model for the development of the lateral ramp

response of seismicity to coulomb stress triggers and shadows of the 1999 mw 7 6 chi chi taiwan earthquake 1 the correlation between static coulomb stress increases and aftershocks has thus far provided the strongest evidence that stress changes promote seismicity a correlation that the chi chi earthquake well exhibits several studies have deepened the argument by resolving stress changes on aftershock focal mechanisms which removes the assumption that the aftershocks are optimally oriented for failure here one compares the percentage of planes on which failure is promoted after the main shock relative to the percentage beforehand for chi chi we find a 28 increase for thrust and an 18 increase for strike slip mechanisms commensurate with increases reported for other large main shocks however perhaps the chief criticism of static stress triggering is the difficulty in observing predicted seismicity rate decreases in the stress shadows or sites of coulomb stress decrease detection of sustained drops in seismicity rate demands a long catalog with a low magnitude of completeness and a high seismicity rate conditions that are met at chi chi we find four lobes with statistically significant seismicity rate declines of 40 90 for 50 months and they coincide with the stress shadows calculated for strike slip faults the dominant faulting mechanism the rate drops are evident in uniform cell calculations 100 month time series and by visual inspection of the m 3 seismicity an additional reason why detection of such declines has proven so rare emerges from this study there is a widespread increase in seismicity rate during the first 3 months after chi chi and perhaps many other main shocks that might be associated with a different mechanism

rapid afterslip following the 1999 chi chi taiwan earthquake 1 postseismic displacements of as much as 14 cm were recorded by gps measurements in the 3 months following the mw 7 6 1999 chi chi taiwan earthquake data from 35 continuous and 90 campaign surveyed stations which show continued east over west thrusting are analyzed to estimate the postseismic slip distribution and fault geometry assuming the shallow fault dips 24 e as determined by numerous studies of the mainshock we invert for the deeper fault structure our results show that the fault dip shallows with depth below the hypocenter merging into a nearly horizontal decollement at a depth of 8 12 km the afterslip distribution shows a maximum slip of 25 cm in the hypocentral region at 7 12 km depth as well as significant slip on the lower decollement afterslip is notably absent in the region of maximum coseismic slip consistent with the afterslip being driven by the mainshock stress change

correlated errors in geodetic time series implications for time dependent deformation analysis of frequent trilateration observations from the two color electronic distance measuring networks in california demonstrate that the noise power spectra are dominated by white noise at higher frequencies and power law behavior at lower frequencies in contrast earth scientists typically have assumed that only white noise is present in a geodetic time series since a combination of infrequent measurements and low precision usually preclude identifying the time correlated signature in such data after removing a linear trend from the two color data it becomes evident that there are primarily two recognizable types of time correlated noise present in the residuals the first type is a seasonal variation in displacement which is probably a result of measuring to shallow surface monuments installed in clayey soil which responds to seasonally occurring rainfall this noise is significant only for a small fraction of the sites analyzed the second type of correlated noise becomes evident only after spectral analysis of line length changes and shows a functional relation at long periods between power and frequency of 1 where is frequency and 2 with 2 this type of correlated noise is termed random walk noise and its source is mainly thought to be small random motions of geodetic monuments with respect to the earth s crust though other sources are possible because the line length changes in the two color networks are measured at irregular intervals power spectral techniques cannot reliably estimate the level of 1 noise rather we also use here a maximum likelihood estimation technique which assumes that there are only two sources of noise in the residual time series white noise and random walk noise and estimates the amount of each from this analysis we find that the random walk noise level averages about 1 3 mm yr and that our estimates of the white noise component confirm theoretical limitations of the measurement technique in addition the seasonal noise can be as large as 3 mm in amplitude but typically is less than 0 5 mm because of the presence of random walk noise in these time series modeling and interpretation of the geodetic data must account for this source of error by way of example we show that estimating the time varying strain tensor a form of spatial averaging from geodetic data having both random walk and white noise error components results in seemingly significant variations in the rate of strain accumulation spatial averaging does reduce the size of both noise components but not their relative influence on the resulting strain accumulation model

viscoelastic afterslip concurrence a possible mechanism in the early post seismic deformation of the mw 7 6 1999 chi chi taiwan earthquake summary observed coseismic data as well as 97 days of post seismic gps data for the chi chi earthquake are used as constraints in the modelling of crustal evolution using the 3 d finite element method first the coseismic gps data are used to justify the use of the elastic earth model and the source rupture model subsequently the most likely rheological model is determined by analysing several modelled time dependent displacements for various viscosity structures the range of viscosities of the lower crust in central taiwan is determined in advance from laboratory measurements and the long term strain rate the estimated viscosity of 5 0 10 17 p as seems to be very low and a relaxation time of 116 days seems very short but the latter approximates the gps measurement of 86 days since earlier studies have indicated that both the viscoelastic response model and the afterslip model may affect post seismic deformation we compare theoretical surface displacements for each of the two models that we evaluate the results reveal that there is little doubt that while neither of these models alone is able to predict the gps measurements well in a 97 day period the combination of the two models improves the predictions considerably we conclude that the afterslip mainly dominated chi chi post seismic deformation in the rupture area while the viscoelastic model did so elsewhere

coseismic slip resolution and post seismic relaxation time of the 1999 chi chi taiwan earthquake as constrained by geological observations geodetic measurements and seismicity summary we analysed the surface displacements estimated by global positioning system gps mea surements and systeme probatoire d observation de la terre spot image correlation and the seismicity change induced by the mw 7 6 1999 september 20 chi chi earthquake we first modelled the seismogenic zone for which we introduced a realistic geologically constrained ramp and flat geometry and estimated by joint inversion of the gps and spot data the co seismic slip distribution stressing the resolution of the model parameters the fair adjustment indicates that the ramp and flat geometry is consistent with the geodetic measurement the cal culated coseismic slip increases both towards the surface and northwards with maximum slip of approximately 15 m in good agreement with the observed surface rupture and previously ob tained slip models based on seismic waveforms this distribution may be represented to a good approximation by a linear increase of slip with a constant azimuth in the plate convergence direction special emphasis is given to the analysis of the model resolution which sharply decreases with depth in spite of the exceptionally well distributed and high quality inverted data this resolution decrease expected for any slip distribution based on surface observations only led us then to examine the seismicity as independent data to test our model following the static stress triggering hypothesis the model is accurately defined in terms of coulomb stress change in the central part of the rupture where the fault geometry is well constrained by seismic profile in contrast at the southern and northern ends it is less appropriate however the static stress triggering approach is not precise enough to determine whether fault geometry or slip distribution is less suitable we thus conclude that the analysis of seismicity in terms of coulomb stress change might not allow us to examine small scale features of seismogenic zones finally we discuss post seismic relaxation from the temporal variation of seismicity rate in contrast to the relaxation time of 0 23 0 90 yr inferred from post seismic permanent gps measurements we found that the decay of aftershock sequence follows omori s law with an 8 14 yr relaxation time we suggest that this discrepancy results from the coexistence of several different processes not all associated with the production of aftershocks

a transient subduction zone slip episode in southwest japan observed by the nationwide gps array 1 we investigate the spatiotemporal variation of fault slip on the subduction interface in southwest japan following the two 1996 hyuganada earthquakes m 6 7 using continuous gps observations a silent thrust slip event has been reported by hirose et al 1999 following these earthquakes approximately 100 km to the north beneath the bungo channel we employ the network inversion filter to rigorously invert for the fault slip distribution in space and time the gps time series is modeled as a sum of 1 secular deformation 2 coseismic station displacements 3 transient fault slip 4 random benchmark motion 5 reference frame shifts and 6 observational errors nonnegativity constraints effectively prohibit normal faulting on the subduction interface spatial and temporal smoothing parameters are estimated by the extended kalman filter the estimated fault slip rate images postseismic slip following the two hyuganada earthquakes slip associated with the silent thrust event began roughly 1 month after the second earthquake but initiated 100 km to the northeast of the hyuganada earthquakes the slow event propagated southwest and downdip beneath the bungo channel lasting approximately 1 year with maximum slip rates of 0 6 m yr the equivalent seismic moment was mw 7 2 we can rule out the possibility that slip propagated from the source region of the hyuganada earthquakes to the location of the silent earthquake however the timing of the events suggests a causal connection perhaps due to static stress transfer this kind of transient event or silent earthquake may reflect different fault zone characteristics and may be responsible for accommodating strain that would otherwise be released in large earthquakes

space time distribution of afterslip following the 2003 tokachi oki earthquake implications for variations in fault zone frictional properties 1 gps time series following the september 25 2003 ut tokachi oki earthquake m w 8 0 reveal significant deformation we use the network inversion filter to invert for the spatial and temporal evolution of afterslip in the 30 days following the earthquake afterslip is concentrated adjacent to the coseismic rupture zone between the 2003 rupture and the inferred source regions of the 1968 tokachi oki and the 1973 nemuro oki earthquakes the inversion shows a rapid decay of slip rate in the 2003 rupture zone with slower deceleration in the surrounding regions the stress velocity paths for the afterslip regions approximately follow d d in v 0 6 mpa suggestive of steady state velocity strengthening friction the complementary spatial pattern of coseismic rupture and afterslip may indicate along strike spatial variations of frictional properties although other interpretations can not be ruled out

ten second love wave propagation and strong ground motions in taiwan applying the surface wave gaussian beam method to a three dimensional 3 d structure of taiwan we have studied the generation and propagation of short period 10 s surface waves in a region of pronounced crustal heterogeneity especially in a region defined to be at an epicentral distance larger than a few focal depths away where high frequency near field motions are largely attenuated and surface waves are adequately developed by perturbing the source model as well as the crustal model in an iterative inversion process we have achieved an excellent fit to the recent strong motion observations recorded during the tung ao earthquake june 5 1994 ml 6 0 in taiwan this allows us to 1 obtain a refined 3 d crustal model of taiwan from an initial model that was derived from recent tomographic results and more important 2 explain the distribution patterns of the strong shaking reported in terms of intensity maps for historical large earthquakes and predict the long period strong motion distribution patterns for future large earthquakes in taiwan with an improved 3 d structure we have gained better insight into short period surface wave propagations as they are modified by crustal lateral heterogeneity focusing and defocusing are clearly demonstrated of particular interest is the central mountain range a massive tertiary metamorphic body forming the backbone of taiwan that seems to serve as a divergent lens for the propagating surface waves whereas the large sedimentary basin in southwestern taiwan seems to focus the wavefronts instead the computed distortions of the propagating surface wave field due to lateral heterogeneity are correctly reflected by the observations

temporal and spatial variations of post seismic deformation following the 1999 chi chi taiwan earthquake summary we use gps displacements collected in the 15 months after the 1999 chi chi taiwan earthquake m w 7 6 to evaluate whether post seismic deformation is better explained by afterslip or viscoelastic relaxation of the lower crust and upper mantle we find that all viscoelastic models tested fail to fit the general features in the post seismic gps displacements in contrast to the satisfactory fit obtained with afterslip models we conclude that afterslip is the dominant mechanism in the 15 month period and invert for the space time distribution of afterslip using the extended network inversion filter our results show high slip rates surrounding the region of greatest coseismic slip the slip rate distribution remains roughly stationary over the 15 month period in contrast to the limited coseismic slip on the

slip distribution and tectonic implication of the 1999 chi chi taiwan earthquake we report on the fault complexity of the large m w 7 6 chi chi earthquake obtained by inverting densely and well distributed static measurements consisting of 119 gps and 23 doubly integrated strong motion records we show that the slip of the chi chi earthquake was concentrated on the surface of a wedge shaped block the inferred geometric complexity explains the difference between the strike of the fault plane determined by long period seismic data and surface break observations when combined with other geophysical and geological observations the result provides a unique snapshot of tectonic deformation taking place in the form of very large 10m displacements of a massive wedge shaped crustal block which may relate to the changeover from over thrusting to subducting motion between the philippine sea and the eurasian plates

time dependent distributed afterslip on and deep below the izmit earthquake rupture surface deformation transients measured with the global positioning system during the 87 days between the 17 august 1999 izmit earthquake and the 12 november 1999 duzce earthquake indicate rapidly decaying aseismic fault slip on and well below the coseismic rupture elastic model inversions for time dependent distributed fault slip using a network inversion filter approach show that afterslip was highest between and below the regions of maximum coseismic slip and propa gated downward to or even below the base of the crust maximum afterslip rates decayed from greater than 2 m yr immediately after the i zmit earthquake to about 1 2 m yr just prior to the duzce earthquake maximum afterslip occurred below the eastern karadere rupture segment and near the i zmit hypocenter afterslip in the upper 16 km decayed more rapidly than that below the seismogenic zone these observations are consistent with a phase of rapid aseismic fault slip concentrated near the base of the seismogenic zone continued loading from the rapid deep afterslip along the eastern rupture zone is a plausible mechanism that helped trigger the nearby mw 7 2 12 november duzce earthquake

a two dimensional dislocation model for interseismic deformation of the taiwan mountain belt we use a global positioning system gps derived surface velocity field of taiwan for the time period between 1993 and 1999 to infer interseismic slip rates on subsurface faults we adopt a composite elastic half space dislocation model constrained by the observed horizontal velocities projected into the direction of plate motion 306 the gps data are divided into northern and southern regions and the velocities in each region are projected into single profiles the model fault geometry includes a shallowly dipping decollement based on the balanced geological cross sections in the coastal plain and western foothills and a two segment fault representing the longitudinal valley fault lvf in eastern taiwan the decollement is composed of two fault segments one extending west under the central range cr and one extending east of the lvf with estimated slip rates of about 35 and 80 mm yr respectively the optimal geometry of decollement is subhorizontal 2 v11 at a depth of 8v9 km the inferred surface location of the western end point of dislocation in the northern profile is located 15 km east of the chelungpu fault while in the southern section it is located beneath the chukou fault the elastic dislocation model successfully matches the horizontal velocity data and predicts elastic strain accumulation in the western foothills that will presumably be released in future earthquakes however considered over multiple earthquake cycles our model cannot explain the topography of the cr and thus fails to predict the active mountain building process in taiwan this failure indicates that both horizontal and vertical velocity fields require a more complex rheological model that incorporates inelastic behavior

viscoelastic shear zone model of a strike slip earthquake cycle i examine the behavior of a two dimensional 2 d strike slip fault system embedded in a 1 d elastic layer schizosphere overlying a uniform viscoelastic half space plastosphere and within the boundaries of a finite width shear zone the viscoelastic coupling model of savage and prescott 1978 considers the viscoelastic response of this system in the absence of the shear zone boundaries to an earthquake occurring within the upper elastic layer steady slip beneath a prescribed depth and the superposition of the responses of multiple earthquakes with characteristic slip occurring at regular intervals so formulated the viscoelastic coupling model predicts that sufficiently long after initiation of the system 1 average fault parallel velocity at any point is the average slip rate of that side of the fault and 2 far field velocities equal the same constant rate because of the sensitivity to the mechanical properties of the schizosphere plastosphere system i e elastic layer thickness plastosphere viscosity this model has been used to infer such properties from measurements of interseismic velocity such inferences exploit the predicted behavior at a known time within the earthquake cycle by modifying the viscoelastic coupling model to satisfy the additional constraint that the absolute velocity at prescribed shear zone boundaries is constant i find that even though the time averaged behavior remains the same the spatiotemporal pattern of surface deformation particularly its temporal variation within an earthquake cycle is markedly different from that predicted by the conventional viscoelastic coupling model these differences are magnified as plastosphere viscosity is reduced or as the recurrence interval of periodic earthquakes is lengthened application to the interseismic velocity field along the mojave section of the san andreas fault suggests that the region behaves mechanically like a 600 km wide shear zone accommodating 50 mrn yr fault parallel motion distributed between the san andreas fault system and eastern california shear zone in the viscoelastic coupling model savage and prescott 1978 time dependent deformation around a strike slip fault zone is realized by loading subsequent rupture and relaxation of a fault present in a schizosphere elastic upper crust overlying a ductile plastosphere assuming a two dimensional 2 d geometry periodic rupture of the fault and identical slip interseismic deformation is obtained as the superposition of viscoelastic relaxation of the plastosphere summed over the infinity of past periodic earthquakes this model yields far field interseismic motions equal to 1 2 of the average slip rate v0 imposed on the fault thus even when no external mechanism is invoked to load the fault i e when the locked portion of the fault penetrates the entire schizosphere the model satisfies the expectation of rigid plate like behavior in the far field in addition when periodic coseismic offsets are factored in the model predicts that the time averaged fault parallel velocity of each point is vol2 depending on location

detecting rock uplift across southern taiwan mountain belt by integrated gps and leveling data abstract rock uplift on the earth surface is a key observation for studies of tectonics and geodynamic processes geodetic measurements from global positioning system gps and leveling are able to track vertical deformation over a wide range of spatial and temporal scale taiwan mountain belt is subject to a rapid uplift rate of 20 mm yr as revealed by gps and leveling measurements in previous studies the extremely high rates motivate us to analyze a set of newly processed gps and leveling data along a nw se transect across southern taiwan by considering the effects of the earthquake related deformation as well as hydrological and surface processes our estimates of rock uplift rates from gps and leveling are from 12 to 14 mm yr across the southern central range of taiwan vertical velocity changes on the small spatial scale are primary related to fault locking in the seismic cycle the large spatial feature showing the mean uplift rate of 6 mm yr between the eastern central range and the frontal thrust faults can be explained by a large scale pop up structure assuming the depth of 15 20 km in the orogenic wedge and the horizontal shortening of 30 mm yr our results draw attention to carefully study geodetic vertical deformation in the mountain belts and seek relevant causes

fault rupture process of the 20 september 1999 chi chi taiwan earthquake we modeled the source rupture process of the 20 september 1999 chi chi taiwan earthquake mw 7 6 using an integrated dataset of near field strong ground motion and global positioning system observations this large inland thrust earthquake produced a surface break more than 80 km long along the chelungpu fault the event triggered over 700 strong motion instruments on the island of tai wan this large ground motion dataset provides the best opportunity to date to un derstand earthquake source processes in a major earthquake first we develop a 3d curved fault model consistent with the surface rupture data and aftershock distribu tion the subsurface fault plane has an average dip of 30 to the east and measures 96 km along strike and 40 km down dip next we applied a genetic algorithm to reconstruct the source rupture process the inversion results indicate a complex slip distribution with a generally broad triangular shaped slip zone that points down dip most of the displacement occurred over a 15 km depth range in the northern segment of the fault with the maximum slip exceeding 20 m the average slip is 3 8 m and the average rupture velocity is 2 6 km sec the mainshock slip distribution clearly complements the aftershock location distribution slip rise times range from 8 to 6 sec over most of the primary slip area and decrease slightly as rupture propagates south to north this relatively long rise time resulted in a low dynamic stress drop and low peak ground accelerations observed at surface stations the rake angle rotates across the rupture plane from predominantly dip slip in the southern segment to mostly left lateral oblique slip in the north we find that geometrical irregularities over the fault plane played an important role in controlling the temporal rupture behavior at fault bends rupture tends to decelerate the rake angle rotates to a near dip slip orientation and the shear stress increases as rupture tries to overcome the bending barriers the total seismic moment from our inversion is 2 9 10 27 dyne cm in good agreement with the harvard and the u s geological survey moment estimation

preseismic deformation and coseismic displacements associated with the 1999 chi chi taiwan earthquake the destructive 1999 chi chi taiwan earthquake m w 7 5 produced an approximately 100 km long surface rupture mostly along the previously recognized north south trending chelungpu fault preseismic deformation in central taiwan is realized from annually repeated global positioning system gps data acquired during the 1992 1999 period the total wnw ese shortening rate in the vicinity of the epicentral region that is from the west coast to the western boundary of the central range is up to 25 mm yr the crustal deformation before the chi chi earthquake was essentially a uniaxial compressional strain of 0 36 strain yr in the direction of 114 the gps measurements taken 0 2 2 7 yr before and within 3 months after the mainshock were utilized to estimate the coseismic displacements horizontal movements of 1 1 9 1 m in the nw nnw directions are observed on the hanging wall eastern side of the fault there is a northward increasing trend in the magnitude of the displacement vectors and a dramatic change in the direction of about 50 toward the east along the fault strike in contrast much smaller se see movements of 0 1 1 5 m are found on the footwall western side of the fault the gps data show 2 4 10 1 m of total horizontal offsets across the chelungpu fault vertical offsets of 1 2 4 4 m with the eastern side up are also observed along the surface rupture the uplift on the hanging wall decreases rapidly toward the east it becomes subsidence at sun moon lake and in the puli wushe area the stations on the footwall show subsidence of 0 02 0 26 m the width of the uplift zone increases from about 10 km in the south to approximately 30 km in the north manuscript received 13 october 2000

time dependent triggered afterslip following the 1989 loma prieta earthquake we employ a modified version of the network inversion filter to investigate time dependent slip following the 1989 loma prieta earthquake previous analysis of global positioning system gps and leveling data suggests afterslip on the loma prieta rupture as well as aseismic slip on a thrust fault northeast of the san andreas fault which we identify with the foothills thrust belt we analyzed 173 daily gps solution files at 62 stations collected from 1989 8 to 1998 3 a total of 1 134 three dimensional relative baseline determinations the observed position changes are assumed to result from secular deformation random benchmark motions and temporally varying fault slip the data reveal temporal variations in slip rate but poorly resolve spatial variations in fault slip the amount of temporal smoothing is estimated by maximum likelihood conditional on this estimate reverse slip on the foothills thrust decays from 45 12 mm yr immediately after the earthquake to zero by 1992 reverse slip on the loma prieta rupture surface decays from 57 11 mm yr to zero by 1994 right lateral slip on the loma prieta rupture surface decays monotonically from 30 10 mm yr to zero by 1994 these results suggest that 1 triggered afterslip can occur off the main rupture zone on adjacent faults 2 shallow afterslip dominated the postseismic deformation for the 8 years following the earthquake and 3 postseismic slip on the foothills thrust may account for a significant portion of its total slip budget

fault geometry and slip distribution of the 1999 chi chi taiwan earthquake imaged from inversion of gps data gps measurements of coseismic displacements from the 1999 chi chi taiwan earthquake are modeled using elastic dislocation theory we find that a single fault plane cannot fit the data but rather a curved fault surface consisting of multiple segments dipping 20 25 best fits the observations the model fault exhibits reverse and left lateral slip on a 75 km long n s trending segment and reverse and right lateral slip on a 25 km e w trending segment at the northern end of the rupture the 21 dipping e w segment is inconsistent with previous interpretations of high angle tear faulting

thermal modeling of continental subduction and exhumation constrained by heat flow and seismicity in taiwan abstract the thermal evolution of crustal exhumation subsequent to the subduction of the continental margin is modeled through numerical solutions to the two dimensional heat conduction equation the boundary conditions used in the modeling are basically constrained by available geophysical and geological observations in the taiwan area temperature distributions are calculated at one million year intervals over a 5 ma duration although this study presents a tentative estimation based on some reasonable assumptions the preliminary results not only explain the sharp change in seismicity within the crust but also agree with the general heat flow pattern on the surface of taiwan more specifically both the high heat flow and the aseismic belt beneath the eastern central range which are probably attributed to higher temperatures within the crust are modeled by crustal exhumation after the subduction of the continental margin however the occurrence of deeper earthquakes east of the central range is a result of the lower thermal gradient which is due to the subduction of the cold continental crust into the uppermost mantle in summary preliminary thermal modeling supports the possibility of a tectonic model of active continental subduction and crustal exhumation for the taiwan orogen

radiometric dating of ash partings in coal of the eocene puget group washington implications for paleobotanical stages new 40 k 40 ar and fission track ages from volcanic ash partings in coal beds of the eocene puget group of western washington indicate a time span of about 41 2 1 8 to 45 0 2 1 2 m y for the 1 890 m section of sediments exposed in the green river area these age data do not entirely support the previous early eocene through early oligocene paleobotanical age estimates for this section and for the four paleobotanical stages defined within it radiometric dating of floras assigned to the same stages outside the type section appears to be partially inconsistent with radiometric ages from the type section

mapping the northern portion of the chelungpu fault taiwan by shallow reflection seismics 1 the chelungpu fault was activated by the 1999 chi chi earthquake mw 7 6 taiwan this fault exhibited extraordinarily large surface ruptures up to 9 8 m as well as underground fault slippages up to 12 m during the earthquake these large displacements were concentrated along the northern portion of the fault 40 km north of the epicenter to prepare data for the future drilling of deep wells in this area many shallow seismic reflection surveys were conducted to investigate the sites an approximate 3d structure of the fault surface can be deduced by this cost effective approach although the depth penetration may be limited e g 3 km the method still provides reliable information to study large ruptures and to better plan future deep wells

gps measurement of postseismic deformation following the 1999 chi chi taiwan earthquake 1 in the 15 month period september 1999 to december 2000 after the 1999 chi chi taiwan earthquake mw 7 6 about 80 campaign surveyed gps stations in central taiwan were occupied up to seven times furthermore seven continuous gps stations were put into operation in the epicentral region mostly within three weeks after the main shock gps data from these near field stations and some pre existing far field stations were utilized to study the postseismic deformation following the chi chi earthquake the postseismic displacements of the gps stations were approximated by a combination of an exponential transient decay function with a relaxation time of 86 days and a postseismic linear rate change stations on the hanging wall displayed west northwesterly horizontal displacements of up to 252 mm and uplift as large as 229 mm postseismic displacements are much larger to the south near the epicenter of chi chi main shock than in the north where the maximum coseismic displacement occurred stations on the footwall show only small displacements the postseismic gps data were inverted to infer the deeper fault geometry and afterslip distribution based on a four segment fault model with its shallow part constrained by coseismic fault geometry deduced from gps and seismological studies the optimal model requires the lower fault segment to be a horizontal plane at a depth of 10 4 km consistent with a model based on gps data taken during the first three months after the main shock and other geological or geophysical studies the maximum afterslip 459 mm occurs in the hypocentral region and in the northern part of the decollement slip of about 430 mm also occurs in the southern part of the decollement the afterslip moment inferred from 15 months of gps data is 4 7 1019 n m approximately 2 3 times the seismic moment released by aftershocks this implies that a major part of the postseismic deformation is aseismic

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