S. Jason Zhao, Researcher in geodesy and geophysics
E-mail address: subduction99@netscape.com or subduction99@yahoo.com
(If you could not get the pictures, it is because someone wants you to
read the interesting advertisements. This is why it is free).
Stress change before large
earthquakes
Horizontal deformation around
Japanese islands
Tectonic stress change in
Japan
Numerical modelling of pore-fluid
flow
Free software
(fault modelling, 3D finite
elements)
(p4)
1. Stress change before the 1976 M=7.8 Tangshan earthquake (China)
Levelling results during 1971 and 1973 show significant height changes in the Tangshan region. The vertical deformation several years before the 1976 Tangshan earthquake could be caused mainly by the movement of the active faults in the region. Inversion analysis of the levelling data indicates that the observed vertical deformation can be well interpreted by the movement of the Jiyunhe fault which truncates the Tangshan earthquake to the southwest. Further stress analysis reveals that shear stress concentration occurred around the Jiyunhe fault at least 3~4 years before the 1976 Tangshan earthquake.
Implications: (1). Aseismic fault movement occurred before the 1976 Tangshan earthquake; (2). The activity of the faults was mainly aseismic and could not be detected by present seismometers; (3). The active fault segments were found beneath the seismogenic layer.
P5
2. Stress change before the 1995 M=7.2 Kobe earthquake (Japan)
Significant horizontal deformation were detected by GPS survey before the 1995 Kobe earthquake in a larger area surrounding the Kobe region. The heterogeneous pattern of the observed horizontal deformation could not be interpreted by a uniform strain model. An inversion analysis of the observed horizontal deformation indicates that an active blind-fault might be present in the region. The blind-fault (model) can explain the observed horizontal deformation reasonably well, and its activity could lead to significant stress localization in the area close to the Kobe earthquake fault. The activity of the blind-fault could be related to the reactivation of the pre-existing rupture planes of historical earthquakes and the aftershocks following the 1946 M=8.2 Nankaido earthquake, or alternatively, the blind-fault could be a portion of the upper face of the subducted plate-slab in the region.
P6
3. Stress change before the 1996 M=7.0 Lijiang earthquake (China)
Inversion analysis of observed GPS baseline and height changes between 1985 and 1992 reveals two active fault segments at the Red River fault zone, Yunnan Province of Southwest China. The modelled active fault segments lead to a clear interpretation of the observed GPS baseline and height changes, and the geometric characteristics of the fault segments are in good agreement with those from geological investigations. The significant activity of the fault segments during this period might be transient, and may reflect the macroscopic, ductile deformation in the mid- to lower crust, which could be responsible for the stress localization in this region. The 3 February 1996 ML = 7.0 Lijiang earthquake and its aftershocks occurred in a north-south band about 45-km long and 16-km wide, truncated by the Jianchuan-Lijiang fault to the south. The aseismic slip of the Red River and Jianchuan-Lijiang faults during 1988-1991 resulted in significant Coulomb stress concentration (at a depth of 10km) around the 1996 earthquake fault (up to 6.0 bar) and southwestern portion of the Jianchuan-Lijiang fault segment (up to 9.5 bar), and a rise in the Coulomb failure stress by 2.6 bar at the hypocenter of the 1996 Lijiang earthquake. These results demonstrate that the 1996 Lijiang earthquake was possibly triggered or induced by aseismic movement of the Jianchuan-Lijiang fault.
.
The Red Rive fault zone and the 1996 Lijiang earthquake
Map showing the geological faults (dashed curves), epicenters of the Lijiang main shock (star, ML =7.0) and aftershocks (dots for ML ³ 4.0; the plus sign for the ML = 6.0 aftershock), the rupture segment (thin solid-line marked by EF), active segments (bold lines for lower margins) of the Red River (RR) and Jianchuan-Lijiang (JL) faults estimated from the joint inversion of GPS and levelling data. The lower-left inset is a sketch map of the Chinese mainland which shows the location (star) of the 1996 Lijiang earthquake.
See (Dx2ya.gif)
(p13)
Distribution of GPS baselines and levelling stations. (also shown are the Red River (RR) and Jianchuan-Lijiang (JL) faults estimated from the joint inversion of GPS and levelling data (bold lines). The stars denote the epicenters of the 1996 earthquakes (ML³ 6.0). Also shown are the levelling points (dots), GPS stations (circles), and baselines (solid lines).
See dx21b.gif
Observed (left) and predicted (right) height changes at the Red River fault zone (1985 - 1992). The contour interval is 15 mm. Also shown are the epicenters (stars) of the 1996 earthquakes (ML ³ 6.0). The squares (left) denote levelling points, and the straight lines (right) represent the two active segments (marked by RR and JL, respectively) estimated from the joint inversion.
See hsp2ar.gif
Locations of the active segments at depth for (a) the Red River fault and (b) the Jianchuan-Lijiang fault. Also shown in (b) are hypocenters (circles) of the ML³ 4.0 aftershocks in the 1996 Lijiang earthquake sequence which occurred within 10 km of the JL fault trace. The star denotes the hypocenter of the ML = 6.0 aftershock which occurred on the JL fault. The standard error in hypocentral locations is estimated to be less than ± 5km.
See dxdep1.gif
Coulomb stress (CCFF) generated by aseismic movement of the Red River (RR) and Jianchuan-Lijiang (JL) faults during 1988-1991 (at a layer of 10km deep). The 1996 earthquake fault (straight line marked by EF) is represented as a north-south trending dislocation (normal dip-slip) extending from the surface to a depth of 10 km. Also shown are the epicenters of the main shock (star) and aftershocks (triangles, ML ³ 4.0), and active fault segments estimated from the joint inversion (straight lines marked by RR and JL, respectively). The plus sign denotes the forecasted epicenter in 1993 (Zhao et al., 1993). The contour interval is 1.5 bar.
See dxcl10r.gif
Annual horizontal deformation around Japanese islands
Some remarks:
The horizontal deformation along the Japanese
islands has been well determined by GPS. The data set can be directly used
to evaluate the plate subduction models which were commonly used for the
Japan subduction zone. In the past, some kinematic models for simulating plate
subduction along the Japanese islands were mainly constrained by vertical
deformation observed along the Japan subduction zone. Recent studies show
that neither could these models interpret the observed horizontal deformation,
nor predict correctly the stress change assiciated with the plate subduction
process. Future kinematic or mechanic models for simulating plate subduction
along the Japanese islands should be able to account for the pattern of the
horizontal deformation observed by GPS.
View horizontal deformation. Please select an area:
Back to topor Go to
stress map
Annual horizontal deformation -whole Japan(t1)
Observational data were from the Japanese permanent GPS network (www.gsi-mc.go.jp). Note that the station Tsusima was used as a fixed point.
See (jpg97b.gif)
Annual horizontal deformation - Kyushu region (t2)
Note: The short solid-lines denote the displacements computed with a plate subduction model. See the Section of Stress Change for more details.
Annual horizontal deformation - Chugoku and Kinki region (t3)
Note: The short dashed-lines denote the displacements computed with a plate subduction model. See the Section of Stress Change for more details.
Annual horizontal deformation - Tokai and Kanto region(t4)
See (tk7a.gif)
Note: The short solid-lines denote the displacements computed with a plate subduction model. See the Section of Stress Change for more details.
Annual horizontal deformation - northern Honshu (Tohoku)(t5)
See (hok12.gif)
Note: The short dotted-lines denote the displacements computed with a plate subduction model. See the Section of Stress Change for more details.
Annual horizontal deformation - Hokkaido(t6)
See (bht5.gif)
Note: The short dotted-lines denote the displacements computed with a plate subduction model. See the Section of Stress Change for more details.
(Annual stress change around Japanese islands)
By S. Zhao, 1998.
Some remarks:
Efforts have been paid to establish a 3D-plate subduction model to explain the horizontal deformation observed along the Japanese islands. The model was constructed based on the geometry of the subducted slabs and the rates of the plate motion in the region. The parameters of the model were estimated by inversion of the observed horizontal deformation. The model has provided a reasonable fit (within ± 5mm on an average) to the horizontal deformationcomponents observed at 504 GPS stations and is in good agreement with the principle of the present plate tectonic theory, which justify the model as well as the resultant stress estimates. However, the geometry of the subducted slabs at depth, which were inferred mainly from the depth distribution of earthquakes, could be subject to some large uncertainties, hence the stress changes evaluated by the model as well as the stress maps displayed below may only be taken as semi-quantitative results.
In general, the stress map may be used as a preliminary guide to the assessment of the long-term seismic hazards along the Japan subduction zone. In addition, the stress map would be further improved and updated with input of more GPS data as well as refined geometry of the subducted slabs at depth.
Note that the activity of inland-faults was not considered in the present study, and hence one should be cautious when trying to use the results to explain the local stress changes or seismicity.
Acknowledgements: Part of the research work was done in Kyoto University during the author held a visiting fellowship from the Japan Society for the Promotion of Science (JSPS). The author would like to thank Prof. Shuzo Takemoto and other colleagues in the Department of Geophysics, Kyoto University for their kind assistance.
View stress map. Please select a region:
Sorry, this subject is under construction.
Zhao S. et al., Deformation and stress localization in southwest Japan, Earth Planet. Sci. Lett., 2003, in press.
Zhao S. & R. D. Mueller, Three-dimensional
finite element modelling of tectonic stress field in continental Australia,
Joint Publication of Geological Societies of America and Australia,
2003, in press.
Zhao S. et al., 3D finite element modelling of
the Northwest Australian stress field, Proceedings of West Australian
Basins Symposium (WABS) III, October 20-23, Perth, 2002.
Zhao S. & R. D. Mueller, Effect of crustal
heterogeneities on deformation and stress change associated with faulting,
Proceedings of the International Workshop on Physics of Active Fault,
E. Fukuyama and R. Ikeda (Eds.), 26-27 Feb., 2002, Tsukuba, Japan, pp. 333-340,
2002.
Zhao S. & Takemoto S., Deformation and stress change associated with plate interaction at subduction zones, Geophys. J. Int, 142, 300-318, 2000.
Zhao S. & S. Takemoto, Aseismic fault movement before the 1976 Tangshan earthquake detected by levelling survey: Implication for preseismic stress localization? Geophys. J. Int., 136, 68-82, 1999.
Zhao S., D. Chao, X. Lai, J. Xu, and G. Seeber, The 1996 ML=7.0 Lijiang earthquake, Yunnan, China: An anticipated event, J. Geodynamics, 27, 529-546, 1999.
Zhao S. & S. Takemoto, Aseismic fault movement before
the 1995 Kobe earthquake detected by GPS survey:
Implication for preseismic stress localization? Geophys.
J. Int., 135, 595-606, 1998.
Zhao S. & S. Takemoto, Fault geometry and stress
release pattern of the 1995 Hyogo-ken Nanbu earthquake,
Japan estimated by inversion of GPS data, J. Geod.
Soc. Japan, 43, 33-43, 1997.
Zhao S., Joint inversion of observed gravity and GPS
baseline changes for the detection of the active fault segment at the Red
River fault zone. Geophys. J. Int., 122, 70-88, 1995.
Observed and modelled vertical deformation in the Tangshan region between 1971-1972 (Rate in unit of mm/year; sequares denote the levelling stations; JF=Jiyunhe fault; TF=Tangshan fault; stars=earthquake epcenters)
Observed vertical deformation of 1971-72 (dl712r.gif)
Modelled vertical deformation of 1971-72 (el712r.gif) (the short-line, marked J1, is the found active segment on the Jiyunhe fault; note comparing the observed and modelled rates)
------------------------------------------------------------------------
Observed and modelled vertical deformation in the Tangshan region between 1972-1973 (Rate in unit of mm/year; sequares denote the levelling stations; JF=Jiyunhe fault; TF=Tangshan fault; stars=earthquake epcenters)
Observed vertical deformation of 1972-73 (dl723r.gif)
Modelled vertical deformation of 1972-73 (el723r.gif)(The bold-line, marked J2, is the found active segment on the Jiyunhe fault; note comparing the observed and modelled rates).
p10
Active fault (Jiyunhe fault) at depth
See (tangv1.gif)
p2
Shear stress concentration before the 1976 Tangshan earthquake (depth=20km) (unit in bar/year; TF=Tangshan fault; JF=Jiyunhe fault; J2 = the active segment of JF for the period of 1972-73)
See (ssh20r.gif)
P7.
Horizontal deformation before the 1995 Kobe earthquake (The dotted-line denotes the 1995 Kobe earthquake fault; the bold-line is the assumed active blind-fault; the arrows with error ellipses are the observed horizintal deformation vectors (the dashed-lines for the modelled vectors).
See (kco8a.gif)
p8
The active blind-fault detected (The dashed rectangle denotes the location of the active blind- fault at depth; the circles are the locations of the micro-earthquakes (M>2; Depth<30km; 1994/01-12); the star denotes the hypocenter of the 1995 Kobe earthquake).
See kco9.gif
p9
Shear stress change before the 1995 Kobe earthquake (depth=20km) (Unit in bar/yr; the bold-line denotes the active blind-fault; the shear stress = sigma xy, x = south-north, and y = west-east)
See kcos2.gif
P12
The aftershocks (circles) of the 1946 M=8.2 Nankaido earthquake and their possible relationship with the active blind-fault (straight-line)
See (nankdo1.gif)
Annual stress change - the whole Japan (s1)
(Depth=15km)
See (jstr7r.gif)
Captions: Yellow curves = plate boundaries; Blue-curves = coastal lines; EU= Eurasia Plate; PH = Philippine Sea Plate; PA= Pacific Plate; and NA = North America Plate.
Implication: Seismic risks in the areas with the stress change rate > 0.5 bar/yr could be relatively higher than in other areas.
Annual stress change - Kyushu region (s2)
(Depth=30 km)
See (jzh11a.gif)
Captions: Stars = historical earthquakes (M>=7.0); Yellow-line = Nankai trough; Triangles = Recent earthquakes (M>6.0, 1996/12-1997/12).
Note the areas with the stress change rate > 5.0 (or 4.0) bar/yr.
Want to see the subducted slab? (jzh9596.gif).
See horizontal defomation produced by the subducted slabs.
Annual stress change - Chugoku and Kinki region (s3)
(Depth=20km)
See (nank29a.gif)
Caption: Stars = historical earthquakes (M>=7.0).
Note the areas with the stress change rate > 0.6 bar/yr.
Want to see the subducted slab? Kick (zshi3a.gif)
See horizontal defomation produced by the subducted slabs.
Annual stress change - Tokai and Kanto region (s4)
(Depth=15km)
See tkp1.gif
Captions: Stars = historical earthquakes (M>=7.0); Yellow-curves = plate boundaries.
Note the area with the stress change rate > 0.9 bar/yr.
Want to see the subducted slab? Kick (tk4.gif)
Note: The geometry of the subducted slabs in the Tokai and Kanto region is quite complex, and could not be well displayed here. The geometry of the subducted slabs in the region estimated by Ishida (1992, MODEL-91) was used in the computation.
See horizontal defomation produced by the subducted slabs.
Annual stress change - Northern Honshu (Tohoku) (s5)
(Depth=15 km)
See (Th15r.gif)
Captions: Stars = historical earthquakes (M>=7.0); Yellow-curve = plate boundary.
Note the area with the stress change rate > 0.3 bar/yr.
Want see the subducted slab? Kick (hok7a.gif).
See horizontal defomation produced by the subducted slabs.
Annual stress change - Hokkaido
(Depth=15 km)
See (bht15r.gif)
Captions: Stars = historical earthquakes (M>=7.0); Yellow-curve =plate boundary; PA = Pacific Plate; NA = North America Plate.
Note the area with the stress change rate > 0.4 bar/yr.
Want to see the subducted slab? Kick (hok6a.gif).
See horizontal defomation produced by the subducted slabs.
Ando M. (1975): Source mechanisms and tectonic significance of historical earthquakes along the Nankai trough, Japan, Tectonophysics, 27, 119-140.
DeMets, C.R., R.G. Gordon, D. Argus, and S. Stein (1990): Current plate motions, Geophys. J. Int., 101, 425-478.
Forsyth D. and S. Uyeda (1975): On the relative importance of the driving forces of plate motion, Geophys. J. R. astr. Soc., 43, 163-200.
Geographic Survey Institute (ed.) (1997): Report of the Coordinating Committee for earthquake prediction, 58, 90-99 (in Japanese).
Geographic Survey Institute (ed.) (1998): Report of the Coordinating Committee for earthquake prediction, 59, 123-125 (in Japanese).
Hyndman R.D., K. Wang, and M. Yamano (1995): Thermal constraints on the seismogenic portion of the southwestern Japan subduction thrust, J. Geophys. Res., 100, 15373-15392.
Ishida M. (1992): Geometry and relative motion of the Philippine sea plate and Pacific plate beneath the Kanto-Tokai district, Japan, J. Geophys. Res., 97, 489-513.
Miyazaki S., T. Saito, M. Sasaki, Y. Hatanaka and Y. Iimura, Expansion of GSI’s nationwide GPS array, Bull. Geograph. Surv. Inst., 43., 23-34, 1997.
Seno T., S. Stein, and A.E. Gripp (1993): A model for the motion of the Philippine sea plate consistent with NUVEL-1 and geological data, J. Geophys. Res., 98, 17941-17948.
Seno T., T. Sakurai, and S. Stein, Can the Okhotsk plate be discriminated from the North American plate? J. Geophys. Res., 101, 11304-11315, 1996.
Usami T., Study of historical earthquakes in Japan, Bull. Earthq. Res.
Inst., Tokyo Uni., 54, 399-439, 1979.
Free software for geodynamic or mechanical
modelling
(faulting, fault motion, elastic, visco-elastic, viscous
modelling, and finite element method)
0. Introduction to FEM and programs (pdf file)
1. Fortran program for modelling faulting/dislocation
(based on elastic dislocation theory)
a. Forward modelling (readme, program's zip files)
b. Inverse modelling (readme, program's zip files)
Applications: crustal deformation and stress
change associated with fault motion...
c. Forward and inverse modelling of gravity
changes due to dislocation (faulting)
(readme, program's zip files)
2. Finite element program for 3D deformation and stress modelling
(elastic, viscoelastic, viscous/Newtonian flow)
a. visco-elastic model (readme, program's zip files)
b. viscous flow model (readme, program's zip files)
Applications: Crustal and lithosphere deformation,
mantle flow, and earthquake/faulting...
3. others
Last Revised: 2002-09-23