Large Slip Region Research Articles

Updated concepts of seismic gaps and asperities to assess great earthquake hazard along South America

So far in this century, six very large-magnitude earthquakes (MW ≥ 7.8) have ruptured separate portions of the subduction zone plate boundary of western South America along Ecuador, Peru, and Chile. Each source region had last experienced a very large earthquake from 74 to 261 y earlier. This history led to their designation in advance as seismic gaps with potential to host future large earthquakes. Deployments of geodetic and seismic monitoring instruments in several of the seismic gaps enhanced resolution of the subsequent faulting processes, revealing preevent patterns of geodetic slip deficit accumulation and heterogeneous coseismic slip on the megathrust fault. Localized regions of large slip, or asperities, appear to have influenced variability in how each source region ruptured relative to prior events, as repeated ruptures have had similar, but not identical slip distributions. We consider updated perspectives of seismic gaps, asperities, and geodetic locking to assess current very large earthquake hazard along the South American subduction zone, noting regions of particular concern in northern Ecuador and Colombia (1958/1906 rupture zone), southeastern Peru (southeasternmost 1868 rupture zone), north Chile (1877 rupture zone), and north-central Chile (1922 rupture zone) that have large geodetic slip deficit measurements and long intervals (from 64 to 154 y) since prior large events have struck those regions. Expanded geophysical measurements onshore and offshore in these seismic gaps may provide critical information about the strain cycle and fault stress buildup late in the seismic cycle in advance of the future great earthquakes that will eventually strike each region.

High‐Resolution Finite Fault Slip Inversion of the 2019 Ridgecrest Earthquake Using 3D Finite Element Modeling

AbstractThe 2019 Ridgecrest earthquake sequence manifested as one of the most complex fault surface ruptures observed in California in modern times. The M6.4 foreshock and M7.1 mainshock occurred on an intricate network of orthogonal and sub‐parallel faults resulting in observable surface displacement and surface rupture captured by geodetic data. Here we present the application of a high‐resolution 3D finite element model (FEM) approach to invert for the detailed fault slip of the entire sequence using complex rheology and fused coseismic Global Navigation Satellite System (GNSS) data with Sentinel‐1 differential interferometric synthetic aperture radar and pixel offset data. The heterogeneous FEM and the fused geodetic data set of pixel offsets, interferograms, and GNSS data results in our optimal inversion solution. This preferred solution is a complex, high‐resolution non‐planar slip model of both the M6.4 and M7.1 events that features three main regions of large slip (6.9+ m), with depths ranging from 2 to 10 km. The regions of slip are bounded by the mainshock hypocenter and the mainshock aftershocks and appear to be related to spatially varying rheological properties. We successfully reproduce a localized region of observed subsidence in the northern portion of the primary fault through the inclusion of a curved fault strand with a significant dip‐slip component. The curved fault strand is the site of our maximum slip of 7.4 m at a depth of 4.2 km. The results demonstrate a robust fit from a more complete, detailed model for the entire seismogenic zone with reasonable computational cost, providing new insights into the governing rheologic and structural processes.

Fault‐Valve Behavior Estimated From Intensive Foreshocks and Aftershocks of the 2017 M 5.3 Kagoshima Bay Earthquake Sequence, Kyushu, Southern Japan

AbstractDetermining fluid migration and pore pressure change within the Earth is key to understand earthquake occurrences. We investigated the spatiotemporal characteristics of the intense foreshocks and aftershocks of the 2017 ML 5.3 earthquake in Kagoshima Bay, Kyushu, southern Japan, to examine the physical processes governing this earthquake sequence. Our relocated hypocenters show the foreshocks moved on a sharply defined plane with a steep dip. The mainshock rupture initiated at the edge of the foreshock seismic gap. The size of the foreshock seismic gap is comparable to that of the mainshock estimated from the source corner frequency, suggesting this seismic gap corresponds to the large slip region of the mainshock. The aftershocks migrated upward along several steeply dipped planes with a seismicity pattern that deviated from the typical mainshock–aftershock type. This deviation of seismicity pattern, together with the hypocenter migrations, suggests aseismic processes, such as pore pressure migration and aseismic slip, affected this earthquake sequence. We established the following hypothesis. First, fluids originating from the subducting slab migrated upward and intruded into the fault plane, reducing the fault strength and causing the foreshock sequence and potentially aseismic slip. Second, the mainshock rupture occurred due to the decreased fault strength and the increased shear stress in an area with relatively high strength. Third, pore pressure increase associated with post‐failure fluid discharge caused the upward aftershock migration. These observations are consistent with the fault‐valve model and show the importance of fluid movement at depth not only in earthquake swarms but also in foreshock–mainshock–aftershock sequences.

Rupture Process of the Mainshock of the 2019 Ridgecrest Earthquake Sequence from Waveform Inversion with Empirical Green’s Functions

ABSTRACTThe 2019 Mw 7.1 mainshock of the Ridgecrest earthquake sequence, which was the first event exceeding Mw 7.0 in California since the 1999 Hector Mine earthquake, caused near-fault ground motions exceeding 0.5g and 70 cm/s. In this study, the rupture process and the generation mechanism of strong ground motions of the mainshock were investigated through waveform inversions of strong-motion data in the frequency range of 0.2–2.0 Hz using empirical Green’s functions (EGFs). The results suggest that the mainshock involved two large slip regions: the primary one with a maximum slip of approximately 4.4 m was centered ∼3 km northwest of the hypocenter, which was slightly shallower than the hypocenter, and the secondary one was centered ∼25 km southeast of the hypocenter. Outside these regions, the slip was rather small and restricted to deeper parts of the fault. A relatively small rupture velocity of 2.1 km/s was identified. The robustness of the slip model was examined by conducting additional inversion analyses with different combinations of EGF events and near-fault stations. In addition, using the preferred slip model, we synthesized strong motions at stations that were not used in the inversion analyses. The synthetic waveforms captured the timing of the main phases of observed waveforms, indicating the validity of the major spatiotemporal characteristics of the slip model. Our large slip regions are also generally visible in the models proposed by other researchers based on different datasets and focusing on lower frequency ranges (generally lower than 0.5 Hz). In particular, two large slip regions in our model are very consistent with two of the four subevents identified by Ross et al. (2019), which may indicate that part of the large slip regions that generated low-frequency ground motions also generated high-frequency ground motions up to 2.0 Hz during the Ridgecrest mainshock.

The 25 March 2020 MW 7.5 Paramushir, northern Kuril Islands earthquake and major (MW ≥ 7.0) near-trench intraplate compressional faulting

Large compressional-faulting earthquakes located relatively deep in oceanic lithosphere entering subduction zones are primarily caused by plate bending stress, but their timing, depth extent and size can be influenced by temporally-varying shear stress on the plate boundary. The 25 March 2020 MW 7.5 event in the Pacific plate seaward of Paramushir Island (northern Kuril Islands), is among the largest recorded events of this type. Its rupture extends along a large-slip region in the southwestern portion of the 1952 Kamchatka MW 9.0 rupture zone. This region has somewhat lower interplate coupling than the megathrust fault along Kamchatka to the northeast, but there could be 68 yrs of strain accumulation. The 2020 event is considered in the context of the 24 recorded major (MW≥7.0) near-trench intraplate compressional-faulting events. An updated compilation of temporally varying near-trench intraslab faulting relative to major interplate ruptures indicates that the stress cycles on the plate boundary influence both extensional and compressional near-trench faulting caused by plate bending. Particularly noteworthy are such events seaward of areas presumed to be in an advanced stage of their seismic cycle, including relatively shallow compressional events along the 1944 MW 8.1 Tonankai, Japan rupture zone, along with activity along the 1952 Kamchatka and 1922 Chile rupture zones.

Toward comparable relative locations between the mainshock slip and aftershocks via empirical approaches

The relative locations between mainshocks and their aftershocks have long been studied to characterize the mainshock–aftershock relationships, yet these comparisons may be subjected to biases inherited from various sources, such as the analysis method, data, and model parameters. Here, we perform both a relocation analysis of interplate events to obtain accurate relative centroid locations and a slip inversion analysis of the mainshock slip relative to the relocated events, with some of the relocated events used as empirical Green’s functions, to retrieve the spatiotemporal slip features of the mainshock relative to all of the relocated events. We perform these analyses on the large (6.0 le M_{W} le 7.3) interplate earthquakes that occurred near four giant (M_{W} ge 8.5) megathrust earthquakes: the 2007 M_{W} ;8.5 Bengkulu (Indonesia), 2005 M_{W} ;8.6 Nias (Indonesia), 2010 M_{W} ;8.8 Maule (Chile), and 2011 M_{W} ;9.1 Tohoku-Oki (Japan) earthquakes. Most of the spatiotemporal slip features of the mainshocks are consistently recovered using different empirical Green’s functions. We qualitatively and quantitatively demonstrate that the large interplate aftershocks within 5 years of the four analyzed mainshocks are largely located on the periphery or outside of the large-slip regions of these four giant megathrust earthquakes.

Why do aftershocks occur? Relationship between mainshock rupture and aftershock sequence based on highly resolved hypocenter and focal mechanism distributions

In order to clarify the origin of aftershocks, we precisely analyze the hypocenters and focal mechanisms of the aftershocks following the 2000 Western Tottori Earthquake, which occurred in the western part of Japan, using data from dense seismic observations. We investigate whether aftershocks occur on the mainshock fault plane on which coseismic slip occurred or they represent the rupture of fractures surrounding the mainshock fault plane. Based on the hypocenter distribution of the aftershocks, the subsurface fault structure of the mainshock is estimated using principal component analysis. As a result, we can obtain the detail fault structure composed of 8 best-fit planes. We demonstrate that the aftershocks around the mainshock fault are distributed within zones of 1.0–1.5 km in thicknesses, and their focal mechanisms are significantly diverse. This result suggests that most of the aftershocks represent the rupture of fractures surrounding the mainshock fault rather than the rerupture of the mainshock fault. The aftershocks have a much wider zone compared with the exhumed fault zone in field observations, suggesting that many aftershocks occur outside the fault damage zone. We find that most aftershocks except in and around the large-slip region are well explained by coseismic stress changes. These results suggest that the thickness of the aftershock distribution may be controlled by the stress changes caused by the heterogeneous slip distribution during the mainshock. The aftershock is also distributed within a much wider zone than the hypocenter distribution observed in swarm activity in the geothermal region, which is thought to be caused by the migration of hydrothermal fluid. This result implies a difference in generation processes: Stress changes due to the mainshock contribute primarily to the occurrence of aftershocks, whereas earthquake swarms in the geothermal region are caused by fluid migration within the localized zone. Graphical abstract.

Source rupture processes of the 2016 Kumamoto, Japan, earthquakes estimated from strong-motion waveforms

The detailed source rupture process of the M 7.3 event (April 16, 2016, 01:25, JST) of the 2016 Kumamoto, Japan, earthquakes was derived from strong-motion waveforms using multiple-time-window linear waveform inversion. Based on the observations of surface ruptures, the spatial distribution of aftershocks, and the geodetic data, a realistic curved fault model was developed for source-process analysis of this event. The seismic moment and maximum slip were estimated as 5.5 × 1019 Nm (M w 7.1) and 3.8 m, respectively. The source model of the M 7.3 event had two significant ruptures. One rupture propagated toward the northeastern shallow region at 4 s after rupture initiation and continued with large slips to approximately 16 s. This rupture caused a large slip region 10–30 km northeast of the hypocenter that reached the caldera of Mt. Aso. Another rupture propagated toward the surface from the hypocenter at 2–6 s and then propagated toward the northeast along the near surface at 6–10 s. A comparison with the result of using a single fault plane model demonstrated that the use of the curved fault model led to improved waveform fit at the stations south of the fault. The source process of the M 6.5 event (April 14, 2016, 21:26, JST) was also estimated. In the source model obtained for the M 6.5 event, the seismic moment was 1.7 × 1018 Nm (M w 6.1), and the rupture with large slips propagated from the hypocenter to the surface along the north-northeast direction at 1–6 s. The results in this study are consistent with observations of the surface ruptures.

Complementary slip distributions of the August 4, 2003 Mw 7.6 and November 17, 2013 Mw 7.8 South Scotia Ridge earthquakes

The South Scotia Ridge Transform (SSRT) plate boundary between the Scotia and Antarctic plates experienced large strike-slip earthquakes on August 4, 2003 (Mw 7.6) and November 17, 2013 (Mw 7.8). These events have overlapping aftershock zones, which is unusual. A 36°–45° southward dipping fault zone ruptured with left-lateral displacements in each event along the northern margin of the South Orkney micro-continent near 60°S. Slip distributions for the two events are determined using teleseismic body and surface wave recordings along with constraints from GPS ground motion recordings at station BORC on Laurie Island (South Orkney Islands), just south of the SSRT. The aftershock distributions, high-frequency back-projections, and unconstrained body wave finite-fault inversions permit significant overlap of the 2003 and 2013 slip zones; however, the GPS static displacements resolve differences in the large-slip regions of the two ruptures. The 2013 earthquake sequence along the SSRT initiated with Mw 6.1 (November 13) and Mw 6.8 (November 16) foreshocks located ∼50 km west of the mainshock hypocenter, and had aftershocks extending ∼250 km eastward. The rupture spread primarily eastward at ∼2.5 km/s with a total rupture duration of about 120 s, with two distinct patches of large-slip located northwest and northeast of the South Orkney Islands. The rupture swept past BORC, with high-rate GPS (HRGPS) ground motion recordings capturing the time-varying slip history of the faulting. Traditional GPS data require that the largest-slip region of the shorter rupture in 2003 is located in the gap NNE of BORC between the two patches that ruptured in 2013. There appears to be some overlap of lower slip regions. The complementary slip distributions comprise a relatively uniform offset along this portion of the SSRT, which is one of the most seismically active regions of the entire Antarctic plate boundary.

Source rupture process of Lushan M S7.0 earthquake, Sichuan, China and its tectonic implications

The source rupture process of the MS7.0 Lushan earthquake was here evaluated using 40 long-period P waveforms with even azimuth coverage of stations. Results reveal that the rupture process of the Lushan MS7.0 event to be simpler than that of the Wenchuan earthquake and also showed significant differences between the two rupture processes. The whole rupture process lasted 36 s and most of the moment was released within the first 13 s. The total released moment is 1.9×1019N m with MW=6.8. Rupture propagated upwards and bilaterally to both sides from the initial point, resulting in a large slip region of 40 km×30 km, with the maximum slip of 1.8 m, located above the initial point. No surface displacement was estimated around the epicenter, but displacement was observed about 20 km NE and SW directions of the epicenter. Both showed slips of less than 40 cm. The rupture suddenly stopped at 20 km NE of the initial point. This was consistent with the aftershock activity. This phenomenon indicates the existence of significant variation of the medium or tectonic structure, which may prevent the propagation of the rupture and aftershock activity. The earthquake risk of the left segment of Qianshan fault is worthy of attention.

The February 6, 2013 Mw 8.0 Santa Cruz Islands earthquake and tsunami

The Santa Cruz Islands region has high seismicity near a 90° bend in the boundary between the Pacific and Australian plates. Southward, along the Vanuatu island arc, the Australian plate under-thrusts the Pacific plate in a well-defined subduction zone. Westward, a transpressional, predominantly transform boundary extends to the southern Solomon Islands subduction zone. The Santa Cruz Islands region has upper plate strike-slip and normal faulting, plate boundary under-thrusting, outer rise extensional faulting, and intraplate compressional faulting. On February 6, 2013 the largest under-thrusting earthquake (Mw 8.0) that has been instrumentally recorded in the region ruptured the megathrust. The epicenter (10.738°S, 165.138°E) is about 1° north of epicenters of prior large shallow under-thrusting events in Vanuatu in 1934 (M ~ 7.8) and 1966 (MS 7.9), and there is overlap of all three events' aftershock zones, but not their large-slip regions. At least 10 lives were lost, with 6 more missing, due to tsunami ~ 1.5 m high striking the town of Lata and several villages on the main Santa Cruz Island of Nendö (Ndeni) and a nearby small island Nibanga. Inundation of 500 m flooded the Lata airport. The tsunami was well-recorded by DART buoys spanning an unusually wide three-quadrant azimuthal aperture. Iterative modeling of teleseismic broadband P waves and the deep-water tsunami recordings resolves the slip distribution. There are two large-slip patches with a southeastward rupture expansion at about 1.5 km/s, with the second patch appearing to have ruptured with large slip at the trench. The event has relatively low short-period seismic wave energy release, but a typical overall moment-scaled total energy. The shallow rupture depth may be associated with the low level of short-period energy, and a near-total stress drop may account for a lack of underthrusting aftershocks among a highly productive aftershock sequence that included three events with Mw ≥ 7.0.

Imaging short‐period seismic radiation from the 27 February 2010 Chile (MW 8.8) earthquake by back‐projection of P, PP, and PKIKP waves

Teleseismic short‐period (0.5–5 s) P waves from the 27 February 2010 Chile earthquake (Mw 8.8) are back projected to the source region to image locations of coherent short‐period seismic wave radiation. Several receiver array configurations are analyzed using different P wave arrivals, including networks of stations in North America (P), Japan (PKIKP), and Europe (PP), as well as a global configuration of stations with a broad azimuthal distribution and longer‐period P waves (5–20 s). Coherent bursts of short‐period radiation from the source are concentrated below the Chilean coastline, along the downdip portion of the megathrust. The short‐period source region expands bilaterally, with significant irregularity in the radiation. Comparison with finite fault slip models inverted from longer‐period seismic waves indicates that the regions of large slip on the megathrust are located updip of the regions of short‐period radiation, a manifestation of frequency‐dependent seismic radiation, similar to observations for the great 2011 Tohoku earthquake (Mw 9.0). Back projection of synthetic P waves generated from the finite fault models demonstrates that if the short‐period energy had radiated with the same space‐time distribution as the long‐period energy, back‐projection analysis would image it in the correct location, updip. We conclude that back‐projection imaging of short‐period signals provides a distinct view of the seismic source that is missed by studies based only on long‐period seismic waves, geodetic data, and/or tsunami observations.

In situ measurement of the hydraulic diffusivity of the active Chelungpu Fault, Taiwan

Hydraulic diffusivity controls fluid pressure and hence affects effective normal stress during rupture. Models suggest a particularly spectacular example of fluid pressurization during the Mw = 7.6 1999 Chichi earthquake when pressurization may have reduced high‐frequency shaking in the regions of large slip if the fault was sufficiently sealed. We investigate in situ hydraulic diffusivity which is the key parameter in such models through a cross‐hole experiment. We find a diffusivity of D = (7 ± 1) × 10−5 m2/s, which is a low value compatible with pressurization of the Chelungpu fault during the earthquake. In most poroelastic media, the hydraulic storativity S lies between 10−7 and 10−5, so that the transmissivity T along the fault zone is comprised between 10−11 m2/s and 10−9 m2/s. The corresponding permeability (10−18–10−16 m2) is at most one hundred times larger than the value obtained on core samples from the host rock. The fault zone is overpressurized by 0.06 to 6 MPa, which is between 0.2% and 20% of the lithostatic pressure.

Frictional heat from faulting of the 1999 Chi‐Chi, Taiwan earthquake

The frictional heat generated during an earthquake is thought to be a major portion of the total energy release. However, there have been no direct measurements of the heat generated by a large earthquake. We present an estimate of the frictional heat produced by the 1999 Chi‐Chi, Taiwan earthquake (Mw 7.6), from temperature anomaly that is measured in a borehole penetrating a large slip region (∼10 m) of the fault. A local increase in the temperature profile across the fault is interpreted to be the residual heat generated during the earthquake. Analyses of this temperature anomaly lead to a low estimate of the dynamic shear stress, 0.5 to 0.9 MPa, and the heat produced in the large slip region of the fault, 0.68 to 1.32 × 1017 J. The total frictional energy is estimated as 2.4 to 6.1 × 1017 J, indicating that seismic efficiency is 1 to 3%.

The Estimation of Minimum-Misfit Stochastic Models from Empirical Ground-Motion Prediction Equations

In areas of moderate to low seismic activity there is commonly a lack of recorded strong ground motion. As a consequence, the prediction of ground motion expected for hypothetical future earthquakes is often performed by employing em- pirical models from other regions. In this context, Campbell's hybrid empirical ap- proach (Campbell, 2003, 2004) provides a methodological framework to adapt ground-motion prediction equations to arbitrary target regions by using response spectral host-to-target-region-conversion filters. For this purpose, the empirical ground-motion prediction equation has to be quantified in terms of a stochastic model. The problem we address here is how to do this in a systematic way and how to assess the corresponding uncertainties. For the determination of the model param- eters we use a genetic algorithm search. The stochastic model spectra were calculated by using a speed-optimized version of SMSIM (Boore, 2000). For most of the em- pirical ground-motion models, we obtain sets of stochastic models that match the empirical models within the full magnitude and distance ranges of their generating data sets fairly well. The overall quality of fit and the resulting model parameter sets strongly depend on the particular choice of the distance metric used for the stochastic model. We suggest the use of the hypocentral distance metric for the stochastic simulation of strong ground motion because it provides the lowest-misfit stochastic models for most empirical equations. This is in agreement with the results of two recent studies of hypocenter locations in finite-source models which indicate that hypocenters are often located close to regions of large slip (Mai et al., 2005; Mani- ghetti et al., 2005). Because essentially all empirical ground-motion prediction equa- tions contain data from different geographical regions, the model parameters corre- sponding to the lowest-misfit stochastic models cannot necessarily be expected to represent single, physically realizable host regions but to model the generating data sets in an average way. In addition, the differences between the lowest-misfit sto- chastic models and the empirical ground-motion prediction equation are strongly distance, magnitude, and frequency dependent, which, according to the laws of un- certainty propagation, will increase the variance of the corresponding hybrid empir- ical model predictions (Scherbaum et al., 2005). As a consequence, the selection of empirical ground-motion models for host-to-target-region conversions requires con- siderable judgment of the ground-motion analyst.

Hypocenter Locations in Finite-Source Rupture Models

We use a database of more than 80 finite-source rupture models for more than 50 earthquakes ( M w 4.1–8.1) with different faulting styles occurring in both tectonic and subduction environments to analyze the location of the hypocenter within the fault and to consider the correlation between hypocenter location and regions of large slip. Rupture in strike-slip and crustal dip-slip earthquakes tends to nucleate in the deeper sections of the fault; subduction earthquakes do not show this tendency. Ratios of the hypocentral slip to either the average or the maximum slip show that rupture can nucleate at locations with any level of relative displacement. Rupture nucleates in regions of very large slip ( D ≥ 2/3 D max ) in only 16% of the events, in regions of large slip (1/3 D max D D max ) in 35% of the events, and in regions of low slip ( D ≤ 1/3 D max ) in 48% of the events. These percentages significantly exceed the percentages of fault area with very large (∼7%) and large (∼28%) slip. Ruptures that nucleate in regions of low slip, however, tend to nucleate close to regions of large slip and encounter a zone of very large slip within half the total rupture length. Applying several statistical tests we conclude that hypocenters are not randomly located on a fault but are located either within or close to regions of large slip.

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.

2003年5月26日に宮城県沖で発生したスラブ内地震の震源モデルと強震動シミュレーション

A large intraslab earthquake (MJMA 7.1) with the focal depth of about 70km occurred in the Pacific slab below off Miyagi Prefecture, northeastern Japan on May 26, 2003. This earthquake source brought strong high-frequency ground motion to the near-source region. At first, a source model is estimated by the forward modeling through strong motion simulations using the empirical Green's function method. Strong ground motions are simulated in the broadband frequency range (0.3-10Hz). Four borehole station data of KiK-net are used. A simple source model having several rectangular “strong motion generation areas” (Miyake et al., 2003) is assumed. The fittings of acceleration, velocity, and displacement waveforms are fairly well. Next, strong ground motion simulation is conducted at all the stations of K-NET and KiK-net (surface) within about 100km distance from the epicenter to validate the obtained source model. The spatial distributions of peak horizontal velocity and acceleration are well reproduced, which gives an evidence that our obtained source model can describe principal characteristics of the actual source process of this earthquake. However, synthetic peak horizontal accelerations are much larger than observed ones at several stations, which may be caused by nonlinear effect of the superficial layers during the mainshock. The obtained source model contains three strong motion generation areas. The spatial distribution of the strong motion generation areas is almost corresponding to the large slip regions from the kinematic waveform inversion analysis. The total size of the strong motion generation areas against a given seismic moment follows the empirical relationship with source-depth for intraslab earthquakes by Asano et al. (2003). Each strong motion generation area of this source model has 105MPa of stress drop. The varieties and uncertainties in estimation of the stress drop value on strong motion generation areas are also discussed.

Rupture behavior of the 1999 Chi-Chi, Taiwan, earthquake—slips on a curved fault in response to the regional plate convergence

Simple kinematic modeling of particle motion along a curved fault similar to the ruptured Chelungpu fault, Taiwan, indicates a unique spatial slip pattern. Specifically, we find that the large convergent slip on the curve is a result of the minimum deformation scenario of the fault geometry and regional northwestern movement of the Philippine Sea Plate (PSP). The modeled deformation regime portrays an accumulation of deformation in the curved region, which coincides very well with a long-term observed NW-SE-trending seismogenic zone in the central Taiwan. This consistence suggests that the Chelungpu fault is a preexisting curved fault. This is further evidenced by geological and geophysical observations. Because the spatial slip pattern is locally and regionally tectonically controlled, it indicates that the rupture behavior of the Chi-Chi earthquake is repeatable. Better knowledge of the fault geometry and the regional plate motion may help us to predict the possible spatial slip distribution of large earthquakes. This discovery is important for avoiding large buildings and constructions near predicted large slip regions.

Numerical Analysis of Dusty Hypersonic Viscous Gas Flow over a Flat Plate

We present an analysis of dilute gas-particle flows at hypersonic Mach number in which the flow parameters are obtained by solving complete two-dimensional Navier-Stokes equations numerically for a Bat-plate geometry. The conservation equations for both the gas and particle phases have been solved using finite volume upwind schemes, Time marching is carried out using the Euler-implicit scheme for the gas phase and the Euler-explicit scheme for the particle phase. Roe's flux-difference-splitting scheme is used for the gas phase, and the Steger-Warming Bu-vector-splitting scheme is applied for the particle phase. The detailed Bow structures of the gas and particle phases are given in three distinct regions: the large-slip region near the leading edge, the moderate-slip region, and the small-slip region far downstream. It is found that the presence of solid particles in the flow has an appreciable effect on the flowfield parameters. Also the presence of particles enhance the skin friction and heal transfer along the wall and decreases the boundary-layer thickness compared with the pure-gas case.