Earthquake Report: M 6.0 northeast Pacific Ocean

I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.

Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.

https://earthquake.usgs.gov/earthquakes/eventpage/us7000ilwt/executive

This is possibly one of the most mysterious earthquakes of the year. I forgot to write this up at the time so need to fill in more details after I am done working up my annual summary.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures.

  • In the upper left corner I include a large scale view of the magnetic anomaly data. These anomalies are formed at mid-ocean ridges, so are parallel to the ridges. When transform faults offset these anomalies, the anomalies get offset.
  • In the lower right corner is a map showing the USGS modeled intensity that uses the Modified Mercalli Intensity (MMI) data.
  • In the upper left center and the lower right center I include maps from two papers that show the magnetic anomaly data for the Pacific Ocean.
  • Here is the map with 3 month’s seismicity plotted.

    References:

    Basic & General References

  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
  • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
  • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
  • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
  • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
  • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
  • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
  • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
  • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
  • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
  • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
  • Specific References

Return to the Earthquake Reports page.

Earthquake Report for M 6.9 Earthquake in Taiwan

I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.

There was a magnitude M 6.9 earthquake in Taiwan on 18 September 2022.

https://earthquake.usgs.gov/earthquakes/eventpage/us7000i90q/executive

Taiwan is an interesting place, from a tectonic perspective. There is an intersection of several plate boundary fault systems here. Along the western boundary of Taiwan the Eurasia plate subducts (dives beneath) the Philippine Sea plate forming the Manila trench. This megathrust subduction zone fault system terminates somewhere in central-northern Taiwan.

Intersecting central Taiwan from the east is another subduction zone where the Philippine Sea plate subducts beneath the Eurasia plate, forming the Ryukyu trench.

There was an earthquake in Taiwan in 1999 that has been commemorated by creating a park and museum that preserves some of the evidence of the earthquake. This Chi-Chi earthquake cause lots of damage and, sadly, lots of suffering. In addition, because of the dominance of the computer chip manufacturing industry in Taiwan at the time, the price of computer chips was greatly inflated. The global economy suffered following this earthquake.

This 18 September 2022 M 6.9 earthquake occurred on a crustal fault that strikes (trends) parallel to the coast. Because of the mapped faults, I interpret this to have been a left-lateral strike slip earthquake.

There was a foreshock, a mag M 6.5 earthquake, nearby, the day before.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures.

  • In the upper left corner is a map that shows the plates, their boundaries, and a century of seismicity.
  • In the upper right are two maps that show models of how there may have been landslides or liquefaction because of the earthquake shaking and impacts. Read more about landslides and liquefaction here. I include both the USGS epicenter and the Central Weather Bureau Seismological Center epicenter (which is probably more accurate). However, these ground failure models are based on the USGS epicenter/location.
  • To the left of those two maps is a low angle oblique view of the tectonic plates and how they are oriented relative to each other.
  • Below that figure, in the center, is a map from Chen at al. (2020) that shows the earthquake fault mapping along eastern Taiwan. I place a yellow star in the location of the M 6.9 epicenter (the location of the earthquake on the ground surface).
  • In the lower right corner is a map that shows the ground shaking from the earthquake, with color representing intensity using the Modified Mercalli Intensity (MMI) scale. The closer to the earthquake, the stronger the ground shaking. The colors on the map represent the USGS model of ground shaking. The colored circles represent reports from people who posted information on the USGS Did You Feel It? part of the website for this earthquake. There are things that affect the strength of ground shaking other than distance, which is why the reported intensities are different from the modeled intensities.
  • To the left of the intensity map is a map that shows seismicity from the Central Weather Bureau Seismological Center. The locations of earthquakes from this center are better than those from the USGS since this organization runs a local seismic network (the USGS runs a global network). The local network uses more seismometers than the global network (so can detect more events, in this region).
  • To the left of this seismicity map is a plot that shows how the shaking intensity models and reports relate to each other. The horizontal axis is distance from the earthquake and the vertical axis is shaking intensity (using the MMI scale, just like in the map to the right: these are the same datasets).
  • In the upper left-center is a figure that shows the USGS earthquake slip model. This shows how much the fault slipped in different areas (based on their modeling, not observation). The model shows that there were places that may have slipped over 1.5 meters (5 feet).
  • Here is the map with 3 month’s seismicity plotted.

    Supportive Figures

  • I could not help myself. I am so excited to have this website back up and running, like a fully operational space station, that I include below some additional figures that help us understand the tectonic setting.
  • Here is the low angle oblique view of the plate configuration in Taiwan.

  • Here is the map from Chen at al. (2020) that shows the fault mapping in this area of eastern Taiwan.

  • Geologic map of the Coastal Range on shaded relief (after Wang and Chen, 1993). The Longitudinal Valley Fault (LVF) can be subdivided into the Linding and Juisui locked Fault and the Chihshang and Lichi creeping Fault. Vertical cross-sections of VS perturbation tomography along the AeA0 and BeB0 profiles denote the Central Range, the Coastal Range, and the LVF. EU: Eurasian Plate; PH: Philippine Sea Plate.

  • Here is an oblique view of the plate configuration in this region. This is from Chang (2001).

  • Here is a great interpretation showing how the Island of Taiwan is being uplifted and exhumed. This is from Lin (2002).

  • Needless to say, this is an excellent map showing the complicated faulting of this region. This is from Theunissen et al. (2012).

  • Here is another tectonic interpretation map from here.
  • Here is a great general overview of the tectonics of the region from Shyu et al. (2005). I include their figure caption below the image as a blockquote.

  • A neotectonic snapshot of Taiwan and adjacent regions. (a) Taiwan is currently experiencing a double suturing. In the south the Luzon volcanic arc is colliding with the Hengchun forearc ridge, which is, in turn, colliding with the Eurasian continental margin. In the north both sutures are unstitching. Their disengagement is forming both the Okinawa Trough and the forearc basins of the Ryukyu arc. Thus, in the course of passing through the island, the roles of the volcanic arc and forearc ridge flip along with the flipping of the polarity of subduction. The three gray strips represent the three lithospheric pieces of Taiwan’s tandem suturing and disarticulation: the Eurasian continental margin, the continental sliver, and the Luzon arc. Black arrows indicate the suturing and disarticulation. This concept is discussed in detail by Shyu et al. [2005]. Current velocity vector of the Philippine Sea plate relative to the Eurasian plate is adapted from Yu et al. [1997, 1999]. Current velocity vector of the Ryukyu arc is adapted from Lallemand and Liu [1998]. Black dashed lines are the northern and western limits of the Wadati-Benioff zone of the two subducting systems, taken from the seismicity database of the Central Weather Bureau, Taiwan. DF, deformation front; LCS, Lishan-Chaochou suture; LVS, Longitudinal Valley suture; WF, Western Foothills; CeR, Central Range; CoR, Coastal Range; HP, Hengchun Peninsula. (b) Major tectonic elements around Taiwan. Active structures identified in this study are shown in red. Major inactive faults that form the boundaries of tectonic elements are shown in black: 1, Chiuchih fault; 2, Lishan fault; 3, Laonung fault; 4, Chukou fault. Selected GPS vectors relative to the stable Eurasian continental shelf are adapted from Yu et al. [1997]. A,Western Foothills; B, Hsueshan Range; C, Central Range and Hengchun Peninsula; D, Coastal Range; E, westernmost Ryukyu arc; F, Yaeyama forearc ridge; G, northernmost Luzon arc; H, western Taiwan coastal plains; I, Lanyang Plain; J, Pingtung Plain; K, Longitudinal Valley; L, submarine Hengchun Ridge; M, Ryukyu forearc basins.

  • This figure from Shyu et al. (2005) shows their interpretation of the different tectonic domains in Taiwan. This is a complicated region that includes collision zones in different orientations as the Okinawa Trough, Ryukyu Trench, and Manila Trench (all subduction zones) each intersect beneath and adjacent to Taiwan. I include their figure caption below the image as a blockquote.

  • Map of major active faults and folds of Taiwan (in red) showing that the two sutures are producing separate western and eastern neotectonic belts. Each collision belt matures and then decays progressively from south to north. This occurs in discrete steps, manifested as seven distinct neotectonic domains along the western belt and four along the eastern. A distinctive assemblage of active structures defines each domain. For example, two principal structures dominate the Taichung Domain. Rupture in 1999 of one of these, the Chelungpu fault, caused the disastrous Chi-Chi earthquake. The Lishan fault (dashed black line) is the suture between forearc ridge and continental margin. Thick light green and pink lines are boundaries of domains.

  • This map from Shyu et al. (2005) shows the earthquake slip regions for proposed earthquake scenarios. I include their figure caption below the image as a blockquote.

  • Proposed major sources for future large earthquakes in and around Taiwan. Thick red lines are proposed future ruptures, and the white patches are rupture planes projected to the surface. Here we have selected only a few representative scenarios from Table 1. Earthquake magnitude of each scenario is predicted value from our calculation.

  • This map from here shows the basement geology of Taiwan. Note the accretionary belts, including the forearc basin. This is a compilation from Teng et al. (2001) and Hsiao et al. (1998) as presented in Ustaszewski et al. (2012).

Return to the Earthquake Reports page.

Earthquake Report: M 7.6 Papua New Guinea

I was travelling to southern California to attend the annual meeting for the Southern California Earthquake Center. This was the first in person meeting since 2019, my first SCEC meeting. I had landed and was waiting for the luggage to arrive when I saw the CSEM earthquake app notification that there was a large earthquake in Papua New Guinea (PNG).

I  put together a quick tweet before anything was posted on the USGS earthquake page other than the location and depth.

When I got to my hotel room later, more information was up.

However, due to some problems with Dreamhost (my website hosting company), I am migrating to a different company. For a little while, parts of this website (like the links and the images) will be non-functional. I will not be using Dreamhost again. I don’t have any more to say about this since they have not returned any of my emails in over a week, basically abandoning my website in the middle of the night without any warning.

https://earthquake.usgs.gov/earthquakes/eventpage/us6000iitd/executive

The depth increased to about 90 km and has a normal-oblique sense of motion. This means that the earthquake was the result of a combination of extension (stretching) and strike-slip.

This area is a complicated region from a tectonic perspective. There are old faults and old plate boundaries that may no longer be active and there are known active faults that juxtapose these older structures.

For example, there is a convergent (moving together) plate boundary fault system on the north side of Papua New Guinea. This ‘subduction zone’ formed a deep sea trench called the New Guinea trench, where the Caroline plate subducted south beneath Papua New Guinea. This plate boundary fault system is thought to be inactive on the west side of the island and active, but with a slow convergence rate, on the eastern side of the island.

Then, to the southeast of PNG, there is a deep sea trench formed by the subduction of the Australia plate diving to the north beneath the island. This fault is inactive offshore and turns into the Papua fold and thrust belt (PFTB) onshore. The PFTB onshore is inactive in the western part of the island and has a slow convergence rate on the eastern side of the island. On 25 February 2018 there was a M 7.5 earthquake associated with the PFTB. Here is the earthquake report for that earthquake.

Between these two subduction zones, that dip in opposite directions, are several large strike-slip fault systems, which also have varying levels of activity.

“Yesterday’s” M 7.6 occurred in one of the plates that is/was being subducted. It was probably in the Australia plate that dives beneath PNG (which is responsible for the PFTB).

I updated this page a bit on 13 March 2023 and make some corrections. I note these changes below.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1920-2020 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures. Some of the same figures are located in different places on the larger scale map below.

  • In the lower right corner is a map that shows the earthquake intensity using the modified Mercalli intensity scale. Earthquake intensity is a measure of how strongly the Earth shakes during an earthquake, so gets smaller the further away one is from the earthquake epicenter. The map colors represent a model of what the intensity may be. The USGS has a system called “Did You Feel It?” (DYFI) where people enter their observations from the earthquake and the USGS calculates what the intensity was for that person. The dots with yellow labels show what people actually felt in those different locations.
  • In the lower left center, there is a plate tectonic map showing the tectonic plates and their boundaries (Cloos et al., 2005). I place a yellow star in the general location of this M 7.6 earthquake.
  • On the right, above the intensity map, is a low angle oblique view of the tectonic plate configuration from Cloos et al. (2005). This shows how the Australia plate may be oriented and I place a yellow star to show how the earthquake may be located within the plate (the “slab”).
  • In the upper left corner are two maps showing the probability of earthquake triggered landslides and possibility of earthquake induced liquefaction. I will describe these phenomena below.
  • Here is the map with 3 month’s seismicity plotted.

Other Report Pages

Some Relevant Discussion and Figures

  • Here is the Holm et al. (2016) figure, showing the major plate boundary faults (symbols tell us which plate boundaries are no longer active: dark symbols = active, hollow symbols = inactive).

  • Topography, bathymetry and regional tectonic setting of New Guinea and Solomon Islands. Arrows indicate rate and direction of plate motion of the Australian and Pacific plates (MORVEL, DeMets et al., 2010); Mamberamo thrust belt, Indonesia (MTB); North Fiji Basin (NFB)

  • This is the Cloos et al. (2005) map from the poster.

  • Tectonic map of New Guinea, adapted from Hamilton (1979), Cooper and Taylor (1987), Dow et al. (1988), and Sapiie et al. (1999). AFTB—Aure fold and thrust belt, FTB—fold-and-thrust belt, IOB—Irian Ophiolite Belt, TFB—thrust-and-fold belt, POB—Papuan Ophiolite Belt, BTFZ—Bewani-Torricelli fault zone, MDZ—Mamberamo deformation zone, YFZ—Yapen fault zone, SFZ—Sorong fault zone, WO—Weyland overthrust. Continental basement exposures are concentrated along the southern fl ank of the Central Range: BD—Baupo Dome, MA—Mapenduma anticline, DM—Digul monocline, IDI—Idenberg Inlier, MUA—Mueller anticline, KA—Kubor anticline, LFTB—Legguru fold-and-thrust belt, RMFZ—Ramu-Markham fault zone, TAFZ—Tarera-Aiduna fault zone. The Tasman line separates continental crust that is Paleozoic and younger to the east from Precambrian to the west.

  • This is the Cloos et al. (2005) cross section, showing a different interpretation of the delaminated slab.

  • Lithospheric-scale cross section at 2 Ma. Plate motion is now focused along the Yapen fault zone in the center of the recently extinct arc. This probably occurred because this zone of weakness had a trend that could accommodate the imposed movements as the corner of the Caroline microplate ruptured, forming the Bismarck plate, and the corner of the Australian plate ruptured, forming the Solomon microplate. The collisional delamination-generated magmatic event ends in the highlands as the lower crustal magma chamber solidifies. Upwelled asthenosphere cools and transforms into lithospheric mantle. This drives a slow regional subsidence of the highlands that will continue for tens of millions of years or until other plate-tectonic movements are initiated. Deep erosion is still concentrated on the fl anks of the mountain belt. RMB—Ruffaer Metamorphic Belt, AUS—Australian plate, PAC—Pacific plate.

  • Here is the tectonic map figure from Sappie and Cloos (2004). Their work was focused on western PNG, so their interpretations are more detailed there (and perhaps less relevant for us for these eastern PNG earthquakes).

  • Seismotectonic interpretation of New Guinea. Tectonic features: PTFB—Papuan thrust-and-fold belt; RMFZ—Ramu-Markham fault zone; BTFZ—Bewani-Torricelli fault zone; MTFB—Mamberamo thrust-and-fold belt; SFZ—Sorong fault zone; YFZ—Yapen fault zone; RFZ—Ransiki fault zone; TAFZ—Tarera-Aiduna fault zone; WT—Waipona Trough. After Sapiie et al. (1999).

  • This is the two panel figure from Holm and Richards (2013) that shows how the New Britain trench megathrust splays into three thrust faults as this fault system heads onto PNG. They plot active thrust faults as black triangles (with the triangles on the hanging wall side of the fault) and inactive thrust faults as open triangles. So, either the NG trench subduction zone extends further east than is presented in earlier work or the Bundi Fault Zone is the fault associated with this deep seismicity.

  • Topography, bathymetry and major tectonic elements of the study area. (a) Major tectonic boundaries of Papua New Guinea and the western Solomon Islands; CP, Caroline plate; MB, Manus Basin; NBP, North Bismarck plate; NBT, New Britain trench; NGT, New Guinea trench; NST, North Solomon trench; PFTB, Papuan Fold and Thrust Belt; PT, Pocklington trough; RMF, Ramu-Markham Fault; SBP, South Bismarck plate; SCT, San Cristobal trench; SS, Solomon Sea plate; TT, Trobriand trough; WB,Woodlark Basin; WMT,West Melanesian trench. Study area is indicated by rectangle labelled Figure 1b; the other inset rectangle highlights location for subsequent figures. Present day GPS motions of plates are indicated relative to the Australian plate (from Tregoning et al. 1998, 1999; Tregoning 2002; Wallace et al. 2004). (b) Detailed topography, bathymetry and structural elements significant to the South Bismarck region (terms not in common use are referenced); AFB, Aure Fold Belt (Davies 2012); AT, Adelbert Terrane (e.g. Wallace et al. 2004); BFZ, Bundi Fault Zone (Abbott 1995); BSSL, Bismarck Sea Seismic Lineation; CG, Cape Gloucester; FT, Finisterre Terrane; GF, Gogol Fault (Abbott 1995); GP, Gazelle Peninsula; HP, Huon Peninsula; MB, Manus Basin; NB, New Britain; NI, New Ireland; OSF, Owen Stanley Fault; RMF, Ramu-Markham Fault; SS, Solomon Sea; WMR, Willaumez-Manus Rise (Johnson et al. 1979); WT, Wonga Thrust (Abbott et al. 1994); minor strike-slip faults are shown adjacent to Huon Peninsula (Abers & McCaffrey 1994) and in east New Britain, the Gazelle Peninsula (e.g. Madsen & Lindley 1994). Circles indicate centres of Quaternary volcanism of the Bismarck arc. Filled triangles indicate active thrusting or subduction, empty triangles indicate extinct or negligible thrusting or subduction.

  • Here is the slab interpretation for the New Britain region from Holm and Richards, 2013. I include the figure caption below as a blockquote.

  • 3-D model of the Solomon slab comprising the subducted Solomon Sea plate, and associated crust of the Woodlark Basin and Australian plate subducted at the New Britain and San Cristobal trenches. Depth is in kilometres; the top surface of the slab is contoured at 20 km intervals from the Earth’s surface (black) to termination of slabrelated seismicity at approximately 550 km depth (light brown). Red line indicates the locations of the Ramu-Markham Fault (RMF)–New Britain trench (NBT)–San Cristobal trench (SCT); other major structures are removed for clarity; NB, New Britain; NI, New Ireland; SI, Solomon Islands; SS, Solomon Sea; TLTF, Tabar–Lihir–Tanga–Feni arc. See text for details.

  • Here are the forward models for the slab in the New Britain region from Holm and Richards, 2013. I include the figure caption below as a blockquote.

  • Forward tectonic reconstruction of progressive arc collision and accretion of New Britain to the Papua New Guinea margin. (a) Schematic forward reconstruction of New Britain relative to Papua New Guinea assuming continued northward motion of the Australian plate and clockwise rotation of the South Bismarck plate. (b) Cross-sections illustrate a conceptual interpretation of collision between New Britain and Papua New Guinea.

UPDATE: 2023.03.13

  • Today (13 March 2023) there was an M 6.3 earthquake that made me look back at this M 7.6 event.
  • This area is interesting, in how there are relict plate boundaries surrounded by active ones. The maps shown above make that clear.
  • The series of figures that I want to add are from Holm et al. (2015). They have some excellent cross sections that show their interpretation of the evolution of the plate boundary between the Australia and Pacific plates (including the Bismarck plates as well).
  • Holm et al. (2015) use radiometric ages and geochemical analysis of the rocks in the region to show that there was once a subduction zone, followed by a collision zone, slab foundering, all leading to the existing configuration.
  • This first figure shows (b) the plate boundary faults and topography, (a) their study area and the magmatic arc geology and geochronology, and (c) a large scale map showing their study sites.

  • Panel A: Distribution of mapped magmatic rock units (from Australian Bureau of Mineral Resources (1972)) and major tectonic boundaries of central Papua New Guinea. Published K/Ar (white boxes) and U–Pb ages (grey boxes) for Late Cenozoic magmatic rocks are from Page and McDougall (1972), Grant and Nielsen (1975), Page (1976), Whalen et al. (1982), Rogerson and Williamson (1985), Richards and McDougall (1990), and van Dongen et al. (2010). Panel B: Major tectonic elements of Papua New Guinea; Adelbert Terrane (AT); Aure trough (AuT); Bundi fault zone (BFZ); Bismarck Sea seismic lineation (BSSL); Fly Platform (FP); Finisterre Terrane (FT); Lagaip fault zone (LFZ); Manus Basin (MB); New Britain (NB); New Britain trench (NBT); New Guinea trench (NGT); North Sepik arc (NSA); Owen Stanley fault zone (OSFZ); Papuan Fold and Thrust Belt (PFTB); Papuan Peninsula (PP); Pocklington trough (PT); Ramu–Markham fault zone (RMFZ); South Bismarck plate (SBP); Solomon Sea (SS); Trobriand trough (TT);Woodlark Basin (WB). Panel C: Regional geology of the eastern Papuan Highlands (modified from Australian Bureau of Mineral Resources (1989)).

  • This second figure shows their interpretation of the evolution for this plate boundary in cross sectional view.
  • Time spans 15 million years. The most important figure in this entire report is this one, particularly the lowest panel.
  • The lowest panel shows how there are three plates atop each other and this 7.6 is somewhere within this stack of plates. (Same for the 13 March 2023 earthquake)

  • Simplified geodynamic evolution of the Maramuni arc. Time steps represent: a) north-dipping subduction of the Pocklington Sea slab at the Pocklington trough and associated magmatism of the Maramuni arc intruding the New Guinea Mobile Belt (NGMB); b) onset of continental collision and growth of the New Guinea Orogen from ca. 12 Ma as the Australian continent enters the Pocklington trough. Collision between the Australian continent and the New Guinea Mobile Belt leads to cessation of subduction at the Pocklington trough and associated foundering and steepening of the Pocklington Sea slab in the mantle; c) continued orogenesis results from convergence between the Australian continent and New Guinea Mobile Belt, associated with under thrusting of the leading Australian continental margin at the Pocklington trough. Northward motion of the Australian plate relative to the foundering Pocklington Sea slab leads to buckling and overturning of the slab. A reduction in convergence at the Pocklington trough is accommodated by initiation of north-dipping subduction of the Solomon Sea plate beneath the Finisterre Terrane; d) lithospheric delamination of the Pocklington Sea slab some 6 m.y. after continental collision results in renewed orogenesis in the New Guinea Orogen, and HREE-depleted magmatism at Kokofimpa. Subduction of the Solomon Sea plate continues at the New Britain trench but also begins to underthrust the New Guinea Mobile Belt at the Trobriand trough (Holmet al., unpublished data); e) Maramuni arc magmatism migrates southward into the Papuan Fold and Thrust Belt. Subduction of the Solomon Sea plate at the New Britain and Trobriand subduction systems results in closure of the Solomon Sea and collision and overthrusting of the Finisterre Terrane above the New Guinea Mobile Belt. Magmatism of the West Bismarck arc begins northward of the Finisterre Terrane; and f) convergence between the Finisterre Terrane and New Guinea Mobile Belt continues at the Ramu–Markham fault. Underthrusting of the New Guinea Mobile Belt beneath the Finisterre Terrane results in crustal-derived magmatism in theWest Bismarck arc. The foundered Solomon Sea slab is laterally continuous with the present day Solomon Sea plate to the east (see Fig. 1). Present day cross-section is modified from Holm and Richards (2013).

  • Here is the tweet for the 13 March 2023 M 6.3 earthquake.

Shaking Intensity

Luckily I updated this page because I noticed that the interpretive figure below was incorrect (it was for a different earthquake).

  • Here is a figure that shows a more detailed comparison between the modeled intensity and the reported intensity. Both data use the same color scale, the Modified Mercalli Intensity Scale (MMI). More about this can be found here. The colors and contours on the map are results from the USGS modeled intensity. The DYFI data are plotted as colored dots (color = MMI, diameter = number of reports).
  • In the upper panel is the USGS Did You Feel It reports map, showing reports as colored dots using the MMI color scale. Underlain on this map are colored areas showing the USGS modeled estimate for shaking intensity (MMI scale).
  • In the lower panel is a plot showing MMI intensity (vertical axis) relative to distance from the earthquake (horizontal axis). The models are represented by the green and orange lines. The DYFI data are plotted as light blue dots. The mean and median (different types of “average”) are plotted as orange and purple dots. Note how well the reports fit the green line (the model that represents how MMI works based on quakes in California).
  • Below the lower plot is the USGS MMI Intensity scale, which lists the level of damage for each level of intensity, along with approximate measures of how strongly the ground shakes at these intensities, showing levels in acceleration (Peak Ground Acceleration, PGA) and velocity (Peak Ground Velocity, PGV).

Potential for Ground Failure

Luckily I updated this page because I noticed that the interpretive figure below was incorrect (it was for a different earthquake).

  • Below are a series of maps that show the potential for landslides and liquefaction. These are all USGS data products.
    There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:

    FOS = Resisting Force / Driving Force

  • When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.

  • Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
  • Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
  • An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
  • Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.

  • Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.

  • Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides.
  • I use the same color scheme that the USGS uses on their website. Note how the areas that are more likely to have experienced earthquake induced liquefaction are in the valleys. Learn more about how the USGS prepares these model results here.

    References:

    Basic & General References

  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
  • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
  • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
  • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
  • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
  • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
  • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
  • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
  • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
  • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
  • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
  • Specific References

  • Abers, G. and McCaffrey, R., 1988. Active Deformation in the New Guinea Fold-and-Thrust Belt: Seismological Evidence for Strike-Slip Faulting and Basement-Involved Thrusting in JGR, v. 93, no. B11, p. 13,332-13,354
  • Baldwin, S.L., Monteleone, B.D., Webb, L.E., Fitzgerald, P.G., Grove, M., and Hill, E.J., 2004. Pliocene eclogite exhumation at plate tectonic rates in eastern Papua New Guinea in Nature, v. 431, p/ 263-267, doi:10.1038/nature02846.
  • Baldwin, S.L., Fitzgerald, P.G., and Webb, L.E., 2012. Tectonics of the New Guinea Region, Annu. Rev. Earth Planet. Sci., v. 40, pp. 495-520.
  • Cloos, M., Sapiie, B., Quarles van Ufford, A., Weiland, R.J., Warren, P.Q., and McMahon, T.P., 2005. Collisional delamination in New Guinea: The geotectonics of subducting slab breakoff: Geological Society of America Special Paper 400, 51 p., doi: 10.1130/2005.2400.
  • Dow, D.B., 1977. A Geological Synthesis of Papua New Guinea, Bureau of Mineral Resources, Geology, and Geophysics, Bulltein 201, Australian Government Publishing Sevice, Canberra, 1977, 58 pp.
  • ..

  • Hamilton, W.B., 1979. Tectonics of the Indonesian Region, USGS Professional Paper 1078.
  • Holm, R. and Richards, S.W., 2013. A re-evaluation of arc-continent collision and along-arc variation in the Bismarck Sea region, Papua New Guinea in Australian Journal of Earth Sciences, v. 60, p. 605-619.
  • Holm, R.J., Richards, S.W., Rosenbaum, G., and Spandler, C., 2015. Disparate Tectonic Settings for Mineralisation in an Active Arc, Eastern Papua New Guinea and the Solomon Islands in proceedings from PACRIM 2015 Congress, Hong Kong ,18-21 March, 2015, pp. 7.
  • Holm, R.J., Rosenbaum, G., Richards, S.W., 2016. Post 8 Ma reconstruction of Papua New Guinea and Solomon Islands: Microplate tectonics in a convergent plate boundary setting in Eartth Science Reviews, v. 156, p. 66-81.
  • Johnson, R.W., 1976, Late Cainozoic volcanism and plate tectonics at the southern margin of the Bismarck Sea, Papua New Guinea, in Johnson, R.W., ed., 1976, Volcanism in Australia: Amsterdam, Elsevier, p. 101-116
  • Koulali, A., tregoning, P., McClusky, S., Stanaway, R., Wallace, L., and Lister, G., 2015. New Insights into the present-day kinematics of the central and western Papua New Guinea from GPS in GJI, v. 202, p. 993-1004, doi: 10.1093/gji/ggv200
  • Sapiie, B., and Cloos, M., 2004. Strike-slip faulting in the core of the Central Range of west New Guinea: Ertsberg Mining District, Indonesia in GSA Bulletin, v. 116; no. 3/4; p. 277–293
  • Tregoning, P., McQueen, H., Lambeck, K., Jackson, R. Little, T., Saunders, S., and Rosa, R., 2000. Present-day crustal motion in Papua New Guinea, Earth Planets and Space, v. 52, pp. 727-730.

Return to the Earthquake Reports page.

Earthquake Report: M 6.9 Mid Atlantic Ridge

There have not been that many large earthquakes this year. This is good for one main reason, there is a lower potential for human suffering.

Therefore, there are fewer Earthquake Reports for this year.

This morning (my time) there was a magnitude M 6.9 earthquake along the Romanche transform fault, a right-lateral strike-slip fault system that offsets the Mid Atlantic Ridge in the equatorial Atlantic Ocean. The fault is part of the Romanche fracture zone.

https://earthquake.usgs.gov/earthquakes/eventpage/us7000i53f/executive

The transform faults in this part of the Mid Atlantic Ridge plate boundary have a pattern of earthquakes that seem to max out in the lower 7 magnitudes. This may be (at least partly) due to the maximum length of these faults (?).

The Romanche fault is about 900 kilometers long. The Chain fault is about 250 km long. The St. Paul fault is about 350 km long.

    Earthquake magnitude is controlled by three things:

  1. the size of the earthquake slip area, for most events, this is basically the length of the fault (since the width of the fault is controlled by the thickness of the lithosphere, or the crust)
  2. the amount that the fault slipped
  3. physical properties of the lithosphere or crust on either side of the fault (how “elastic” the Earth is)

Using empirical (data) based relations between earthquake subsurface rupture length and earthquake magnitude (Wells and Coppersmith, 1994), I calculate the maximum earthquake magnitude we may get on these three faults listed above.

Here are the data that Wells and Coppersmith use to establish these relations.


(a) Regression of subsurface rupture length on magnitude (M). Regression line shown for all-slip-type relationship. Short dashed line indicates 95% confidence interval. (b) Regression lines for strike-slip relationships. See Table 2 for regression coefficients. Length of regression lines shows the range of data for each relationship.

Here are the magnitude estimates for each of these fault systems.


Looking at the interpretive poster, we can see that there have not been any temblors that approach the sizes listed in this table. The largest historic earthquake was M 7.1 (there were several).

So, we may ask ourselves one of the most common questions people ask regarding earthquakes. Was this M 6.9 a foreshock to a larger earthquake?

Obviously, we cannot yet know this. Nobody can predict the future (at least not yet).

However, based on the incredibly short historic record of earthquakes, we may answer this question: “no, probably not.” This answer is tempered by the very short seismic record. If magnitude 8 earthquakes occur, on average, every 1000 years, then our ~100 year record might be too short to “notice” one of these M 8 events.
If we continue to look at the historic record, we will see that there appear to be three instances where one of these M 6.5-7 earthquakes had a later earthquake of a similar magnitude.
When an earthquake fault slips, the crust surrounding the fault squishes and expands, deforming elastically (like in one’s underwear). These changes in shape of the crust cause earthquake fault stresses to change. These changes in stress can either increase or decrease the chance of another earthquake.
I wrote more about this type of earthquake triggering for Temblor here. Head over there to learn more about “static coulomb stress triggering.”
In the poster, I label these earthquakes as “Linked Earthquakes.” Perhaps the later of each earthquake pair (or triple) was triggered by the change in static coulomb stress.

    Here are the three sets of “Linked Earthquakes:”

  • In 1992, along the Chain fault, the 16 Feb M 6.6 appears to have triggered the 18 Feb M 6.6. More speculatively, about 6 months later, it seems that there was a triggered M 6.9 on 18 Aug. Static Coulomb triggering typically has a limit of about 2-3 times the rupture length (and this depends of the pre-existing stress on the receiver fault, the fault that may be having triggered slip). A M 6.6 may have a rupture length of 50 km, so could possibly affect faults as far as 100-150 km away. The M 6.9 is about 70 km from the easternmost M 6.6, so it seems possible that the M 6.9 was triggered by the M 6.6.
  • In 2003, along the Romanche fault, there were two M 6.6 earthquakes separated by about 6 weeks. These quakes are about 100 km apart, possibly close enough to be triggered.
  • In 2020, along the St. Paul transform fault, there was a pair of quakes about 3 weeks and 340 km apart. The first quake was M 6.5, so this pair of events seems to far apart to be related.

So, given the historic record, it sure seems likely that there may be another M6-7 earthquakes in the region of the fault sometime in the next couple of months. And, given our lack of knowledge about the long term behavior of these faults, it is also possible that there could be a larger M 8 event.
Since we cannot yet know the real answer to this question, we are reminded of the advice that educators and emergency response people provide: If one lives in Earthquake Country, get earthquake prepared. Just a little effort to get better prepared makes a major difference in the outcome.
Head over to Earthquake Alliance where there are some excellent brochures about how to be better prepared and more resilient to earthquake and tsunami hazards. Living on Shaky Ground is one of my favorites!

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1920-2020 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures. Some of the same figures are located in different places on the larger scale map below.

  • In the lower right corner I include a map that shows the age of the oceanic crust in the Atlantic Ocean. Oceanic crust (or lithosphere) is created at mid ocean ridges, where there is extension that allows upward movement of magma, leading to the formation of oceanic crust. The Mid Atlantic Ridge system is one of these types of plate boundaries.
  • Above the crust age map is an illustration showing the how the crust moving away from the ocean ridges leaves behind oceanic crust. The Earth’s magnetic polarity changes at times and the oceanic crust records these changes in magnetic polarity. These changes are the main reason why we know that the crust is formed along these ridge systems. Read more here.
  • In the upper left corner is a small scale map that shows the historic seismicity, the plate boundary fault systems, and the magnetic anomalies. Places with crust formed when the magnetic field is like today, is colored red (a.k.a. normal polarity) and crust formed when the poles were reversed relative to today is blue (i.e., reversed polarity).
  • In the upper right corner is a map that shows the earthquake intensity from this earthquake (using the modified Mercalli Intensity Scale). Intensity is a measure of how strongly the shaking is felt, not a measure of the earthquake size. So, the intensity gets smaller with distance (see how the highest intensity is nearest the earthquake epicenter).
  • In the lower left center there is a map from Heezen et al. (1964). Heezen was an oceanographer that contributed greatly to our knowledge of the oceans. In this study, one of the things that they were studying is the flow of deep water (deep water flows largely because of changes in density of the seawater, controlled by salinity and temperature). Because of this, they were mapping the shape of the seafloor to see where this deep water could flow. Ths location of this map is outlined by a dashed rectangle in the main map.
  • Here is the map with a month’s seismicity plotted.

Some Relevant Discussion and Figures

  • Here is the Müller et al. (2008) figure from the interpretive poster above.

  • Here a the Bonatti et al. (2001) figure showing the bathymetry of this area. I include the figure caption as a blockquote below.

  • A: Multibeam topography of Romanche region, showing north-south profiles where sampling was carried out. Black dots and red numbers indicate estimated age (in million years) of lithosphere south of Romanche Transform, assuming spreading half-rate of 17 mm/yr within present-day ridge and transform geometry. White dots indicate epicenters of teleseismically recorded 1970–1995 events (magnitude . 4). FZ is fracture zone. B: Topography and petrology at eastern intersection of Romanche Fracture Zone with Mid-Atlantic Ridge. Data were obtained during expeditions S-16, S-19, and G-96 (Bonatti et al., 1994, 1996). C: Location of A along Mid-Atlantic Ridge.

  • Dr. Stephen Hicks and their colleagues conducted a fascinating study of the 2016 M 7.1 earthquake. They hypothesize that the Romanche fault slipped in different parts of the fault at different times (during the earthquake).
  • This map shows the historic seismicity of the region.

  • Seismotectonic context. The map location is given by the red rectangle on the inset globe. Focal mechanisms are shown for events with Mw > 6 (ref. 30). Mw > 7.0 events are labelled. Stations of the PI-LAB ocean bottom seismometer network are indicated by triangles. Our relocated hypocentre and low-frequency RMT of the 2016 earthquake are shown by the red star and red beach ball, respectively. The orange beach ball is a colocated Mw 5.8 used for the Mach cone analysis. The black rectangle shows the location of the map in Fig. 2. ISC Bulletin, Bulletin of the International Seismological Centre.

  • Here is where Hicks et al. (2020) hypothesize that the slip slipped.

  • Interpretation of rupture dynamics for the 2016 Romanche earthquake. Top: perspective view of bathymetry along the Romanche FZ. Bottom: interpretive cross-section along the ruptured fault plane. Colours show a thermal profile based on half-space cooling. The green line denotes the predicted transition between velocity-strengthening and velocity-weakening frictional regimes (as expressed by the a – b friction rate parameter) from Gabbro data35. The numbers show the key stages of rupture evolution: (1) rupture initiation (star) in the oceanic mantle, (2) initiation phase has sufficient fracture energy to propagate upwards to the locked section of fault, (3) weak subshear rupture front travels east in the lower crust and/or upper mantle, (4) rupture reaches the locked, thinner crustal segment close to the weaker RTI (SE1), (5) sufficient fracture energy for a westward supershear rupture in the crust along the strongly coupled fault segment (SE2) and (6) rupture possibly terminated by a serpentinized and hydrothermally altered fault segment.

    References:

    Basic & General References

  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
  • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
  • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
  • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
  • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
  • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
  • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
  • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
  • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
  • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
  • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
  • Specific References

  • Abercrombie, R.E. and Ekstrom, G., 2001. Earthquake slip on oceanic transform faults in Nature, v. 410, p. 74-77
  • Bonatti, E., Brunello, D., Fabretti, P., Ligi, M., Porcaro, R.A., and Sealer, M., 2001. Steady-state creation of crust-free lithosphere at cold spots in mid-ocean ridges in Geology, v. 29, no. 11, p. 979-982.
  • Hicks, S.P., Okuwaki, R., Steinberg, A., Rychert, C.A., Harmon, N. Abercrombie, R.E., Bogiatzis, P., Cataphors, D., Zahradnik, J., Kendall, J-M., Yagi, Y., Shimizu, K., and Sudhaus, H., 2020. Back-propagating supershear rupture in the 2016 Mw 7.1 Romanche transform fault earthquake in Nature Geoscience, v. 13, p. 647-653, https://doi.org/10.1038/s41561-020-0619-9
  • Heezen, B.C., Bunce, E.T., Hersey, J.B., and Tharp, M., 1964. Chain and Romanche fracture zones in Deep-Sea research, v. 11, p. 11-33
  • Müller, R.D., Sdrolias, M., Gaina, C., and Roest, W.R., 2008. Age, spreading rates and spreading symmetry of the world’s ocean crust in Geochem. Geophys. Geosyst., 9, Q04006, doi:10.1029/2007GC001743
  • Torsvik, T.H., Tousse, S., Labaila, C., and Smethurst, M.A., 2009. A new scheme for the opening of the South Atlantic Ocean and the dissection of an Aptian salt basin in Geophysical Journal International, v. 177, p. 1315-1333.

Return to the Earthquake Reports page.


Earthquake Report M 7.0 Philippines

I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.

https://earthquake.usgs.gov/earthquakes/eventpage/us6000i5rd/executive

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures.

  • Here is the map with a month’s seismicity plotted.

20220727_philippines_interpretation.pdf 16 MB pdf

    References:

    Basic & General References

  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
  • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
  • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
  • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
  • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
  • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
  • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
  • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
  • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
  • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
  • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
  • Specific References

Social Media: Here is my thread for this event.

Return to the Earthquake Reports page.

Earthquake Report M 6.9 Macquarie Island

I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.

Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.

https://earthquake.usgs.gov/earthquakes/eventpage/us6000hm9j/executive

I will fill in the details later, after I complete the annual summary. I have not even made a poster for this earthquake yet. Woops.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures.

  • Here is the map with 3 month’s seismicity plotted.

    References:

    Basic & General References

  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
  • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
  • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
  • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
  • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
  • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
  • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
  • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
  • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
  • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
  • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
  • Specific References

Return to the Earthquake Reports page.

Earthquake Report M 6.7 in Taiwan

I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.

Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.

https://earthquake.usgs.gov/earthquakes/eventpage/us6000h6nr/executive

I will fill in these details a bit more after I complete the annual summary.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures.

  • Here is the map with 3 month’s seismicity plotted.

Earthquake Report M 7.3 offshore of Japan

I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.

Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.

https://earthquake.usgs.gov/earthquakes/eventpage/us6000h519/executive

I forgot to write this up when the earthquake happened. I will fill in some details after I have completed my annual summary.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures.

  • Here is the map with a month’s seismicity plotted.

    Ground Failure Poster


    References:

    Basic & General References

  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
  • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
  • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
  • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
  • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
  • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
  • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
  • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
  • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
  • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
  • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
  • Specific References

Return to the Earthquake Reports page.

Earthquake Report: M 6.2 along the Great Sumatra fault

There was a magnitude M 6.2 Gempa or Earthquake on 25 February 2022.

https://earthquake.usgs.gov/earthquakes/eventpage/us6000gzyg/executive

The plate boundary fault system that dominates the tectonics of the Island of Sumatera, Indonesia, is the complicated.

The oceanic India-Australia plate converges with the Eurasia plate to form the Sunda trench. This convergent plate boundary forms a subduction zone where the oceanic plate subducts beneath the continental plate.
Below is a low-angle oblique view cut into the Earth showing this plate configuration from the Earth Observatory Singapore.


However, the direction of plate convergence is not perpendicular to the plate boundary fault (the megathrust subduction zone). Why does this matter?
Because the convergence is at an angle oblique to the plate boundary, we can imagine that this convergence can be subdivided into two components of motion:

  1. the fault perpendicular motion
  2. and the fault parallel motion

The amount of plate convergence that is perpendicular to the plate boundary is accommodated by earthquake fault slip on the megathrust.

The amount of plate convergence that is parallel to the plate boundary is accommodated by earthquake fault slip on a different series of faults that we call sliver faults. The Great Sumatra fault is one of these [forearc] sliver faults.

Here is a figure from Lange et al. (2008) that shows how oblique plate convergence forms both a subduction zone and a forearc sliver fault system.


The M 6.2 earthquake is a strike-slip earthquake along the Great Sumatra fault, one of these forearc sliver faults.

Based on our knowledge of this fault system and the earthquake mechanism, we can easily interpret this to be a right-lateral strike-slip fault.

There are numerous historical analogies from the past century. Most of the events in the past few decades have been in the M 6-7 range, though there have been events of larger magnitude in the past centuries.

Below is my interpretive poster for this earthquake

  • I plot the seismicity from the past 3 months, with diameter representing magnitude (see legend). I include earthquake epicenters from 1922-2022 with magnitudes M ≥ 6.0.
  • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
  • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
  • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

    I include some inset figures. Some of the same figures are located in different places on the larger scale map below.

  • In the upper right corner is a map that shows the major plate boundaries with tectonic strain shown in red and blue. Areas that are red have a higher rate of tectonic deformation due to the motion of these plates and the orientation of the plate boundaries.
  • In the lower center is a low angle oblique view of the Sumatra subduction zone that forms the Sunda Trench. I placed a red circle in the location of the M 6.2 earthquake.
  • In the upper left center there is a map that shows the earthquake shaking intensity. Read more about this further down in this report.
  • In the upper right center is a plot showing earthquake shaking intensity (vertical axis) relative to distance from the earthquake (horizontal axis). This shows a comparison between the USGS shaking models as colored lines (also shown on the map to the left) relative to real reports from real people (the Did You Feel It? (dyfi) points).
  • Below the earthquake intensity map is a map that shows the mapped active faults, along with their slip rates from Natawidjaja (2018).
  • On the right margin are two maps that show models of earthquake triggered landslides and earthquake induced liquefaction. I describe these phenomena later in this report.
  • In the upper left corner are two maps from the Global Earthquake Model program: Seismic Hazard and Seismic Risk. Read more about this later in this report.
  • Here is the map with 3 month’s seismicity plotted.

Some Relevant Discussion and Figures

  • Here is a part of the poster “Seismicity of the Earth 1900-2012, Sumatra and Vicinity” (Hayes et al., 2013). Note the location of Padang, which is southwest of the M 6.2 earthquake.
  • The map shows the location of the seismicity cross section (the next figure). The cross section includes earthquakes from locations within the rectangle. These events are plotted along the line C-C.’
  • See that the Sumatra fault crosses this cross section just to the east of the center.

  • Here is the seismicity cross section C-C’. Many of the earthquakes plotted here follow the subducting slab of the India-Australia plate.
  • However, there are some shallow earthquakes in the upper plate that represent slip along Sumatra fault zone faults (like yesterday’s M 6.2).

  • Further to the north there was a pair of subduction zone earthquakes in 2004 and 2005. The 2004 Sumatra-Andaman subduction zone earthquake was devastating and the earthquake and tsunami led to almost a quarter of million deaths.
  • There are several earthquake report pages on these two events. Here is one of them:/li>
  • This is my poster from that report.

  • Here is a cross section showing where the earthquake hypocenter is compared to where we think the mantle exists. We have not been here, so nobody actually knows… These interpretations are based on industry deep seismic data (Singh et al., 2008 ).

  • Here is the historic rupture map showing the locations of historic subduction zone earthquakes. I include a figure caption below that I wrote as blockquote.

  • India-Australia plate subducts northeastwardly beneath the Sunda plate (part of Eurasia) at modern rates (GPS velocities are based on regional modeling of Bock et al, 2003 as plotted in Subarya et al., 2006). Historic earthquake ruptures (Bilham, 2005; Malik et al., 2011) are plotted in orange. 2004 earthquake and 2005 earthquake 5 meter slip contours are plotted in orange and green respectively (Chlieh et al., 2007, 2008). Bengal and Nicobar fans cover structures of the India-Australia plate in the northern part of the map. RR0705 cores are plotted as light blue. SRTM bathymetry and topography is in shaded relief and colored vs. depth/elevation (Smith and Sandwell, 1997).

  • One of my favorite interpretive posters for this part of the world is from the 2012 outer rise sequence offshore of northern Sumatra. This poster includes additional details about the structure of the India-Australia plate.
  • The India-Australia plate is important as it appears to control some of the tectonics of the megathrust, as well as for the upper plate faults like the Sumatra fault.
  • Read more about the 2012 sequence here.

  • This is a video showing a visualization of the seismic waves transmitted from the 2004 SASZ earthquake from IRIS and others.
  • This movie illustrates simulation of seismic wave propagation generated by Dec. 26 Sumatra earthquake. Colors indicate amplitude of vertical displacement at the surface of the Earth. Red is upward and blue is downward. Total duration of this simulation is 20 minutes. Source model we used is that of Chen Ji of Caltech. Simulation was performed by using the Earth Simulator of JAMSTEC.

  • Here is a great figure from Philobosian et al. (2014) that shows the slip patches from the subduction zone earthquakes in this region.

  • Map of Southeast Asia showing recent and selected historical ruptures of the Sunda megathrust. Black lines with sense of motion are major plate-bounding faults, and gray lines are seafloor fracture zones. Motions of Australian and Indian plates relative to Sunda plate are from the MORVEL-1 global model [DeMets et al., 2010]. The fore-arc sliver between the Sunda megathrust and the strike-slip Sumatran Fault becomes the Burma microplate farther north, but this long, thin strip of crust does not necessarily all behave as a rigid block. Sim = Simeulue, Ni = Nias, Bt = Batu Islands, and Eng = Enggano. Brown rectangle centered at 2°S, 99°E delineates the area of Figure 3, highlighting the Mentawai Islands. Figure adapted from Meltzner et al. [2012] with rupture areas and magnitudes from Briggs et al. [2006], Konca et al. [2008], Meltzner et al. [2010], Hill et al. [2012], and references therein.

  • Now let’s take a closer look at the Sumatra fault. Here is a map that is from Natawidjaja (2017). Dr. Natawidjaja worked with Dr. Kerry Sieh on the Sumatra fault for his Ph.D. research. This map is an updated version showing the different fault segments of the Sumatra fault system. I inlcude their figure caption in blackquote.

  • New revised (simplified) active fault map of the Sumatran Fault Zone (SFZ) according to the PuSGeN Team for Updating Indonesia Seismic Hazard Map (2016) with new slip rates from geological and geodetical (GPS) recent studies.

  • Here is a figure from their dissertation showing these fault segments in greater detail (Natawidjaja , 2002). The M 6.2 earthquake is near the southern boundary of the Barumun segment.

  • Map of 20 geometrically defined segments of the Sumatran fault system and their spatial relationships to active volcanoes, major graben, and lakes.

  • This map shows what the slip rate would be on an hypothetical forearc sliver fault in the location of the Sumatra fault, given plate convergence rates and coral uplift rates (Natawidjaja, 2017).

  • Tectonic modelling based on continuous GPS – SuGAr 9 Sumatran GPS Array) and coral uplift rates,

  • This figure shows a comparison of fault slip rates for different parts of the Sumatra fault (Natawidjaja, 2017).

  • Comparison of GPS velocity profiles across the Sumatran fore arc inferred from (left) kinematic block models (right) with previously published velocity profiles. Modeling all fore-arc site velocities with a single strike-slip fault results in anomalously high inferred slip-rates (>22mm/yr) and missing the Sumatran Fault trace by up to 40km. Incorporating the effect of oblique locking of the Sunda megathrust results in lower inferred slip – rates for the Sumatran Fault (~15mm/yr) that are more consistent with updated geological slip rates.

  • And finally here is an interpretive figure showing how Natawidjaja(2002) interpret the formation of the Sumatra fault system.

  • A plausible (but nonunique) history of deformation along the obliquely convergent Sumatran plate margin, based upon our work and consistent with GPS results and the timing of deformation in the forearc region. (a) By about 4 Ma, the outer-arc ridge has formed. The former deformation front and the Mentawai homocline provide a set of reference features for assessing later deformations. From 4 to 2 Ma, partitioning of oblique plate convergence occurs only north of the equator. Dextral-slip faults on the northeast flank of the forearc sliver plate parallel the trench in northern Sumatra but swing south and disarticulate the forearc basin and outer-arc ridge north of the equator. (b) Slip partitioning begins south of the equator about 2 Ma, with the creation of the Mentawai and Sumatran faults. Transtension continues in the forearc north of the equator. (c) In perhaps just the past 100 yr, the Mentawai fault has become inactive, and the rate of slip on the Sumatran fault north of 2°N has more than doubled. This difference in slip rate may be accommodated by a new zone of transtension between the Sumatran fault and the deformation front in the forearc and outer-arc regions.

  • What about further back in time? Hurukawa et al., 2014 prepared a summary of the earthquakes along the Sumatra fault. Below is a map showing the location of these ruptures.

  • Relocated MJHD epicenters. (a) Northern Sumatra. (b) Central Sumatra. (c) Southern Sumatra. Solid lines with names indicate segments of the Sumatran fault (Sieh and Natawidjaja, 2000). Symbols are as in Figure 2. The thick solid line (see Fig. 4c) indicates the Ranau–Suwoh area, which was severely damaged by the 1933 Liwa earthquake (Berlage, 1934;Widiwijayanti et al., 1996). The slip rates of the Sumatran fault in northern, central, and southern Sumatra are taken from Ito et al. (2012) and Genrich et al. (2000) for Global Positioning System (GPS) and Bellier and Sebrier (1995) for Satellite Pour l’Observation de la Terre (SPOT).

  • This shows these records on a space-time diagram.

  • Earthquake history along the Sumatran fault since 1892. Fault planes estimated in this study are shown by thick lines. SG: Seismic gap.

    Earthquake Stress Triggering

    • When an earthquake fault slips, the crust surrounding the fault squishes and expands, deforming elastically (like in one’s underwear). These changes in shape of the crust cause earthquake fault stresses to change. These changes in stress can either increase or decrease the chance of another earthquake.
    • I wrote more about this type of earthquake triggering for Temblor here. Head over there to learn more about “static coulomb stress triggering.”
    • There are two kinds of earthquake triggering.
      1. Dynamic Triggering – When seismic waves travel through the Earth, they change the stresses in the crust. IF the faults are “locked and loaded” (i.e. they are just about ready to slip in an earthquake), there may be an earthquake on the “receiver” fault. Generally, once the seismic waves are done travelling, this effect is over. Though, some suggest that this affect on the stress changes may last longer (but not much longer).
      2. Static Triggering – When an earthquake fault slips, it deforms (changes the shape) of the crust surrounding that earthquake. These changes can cause increases and decreases in the stress on faults (either increasing or decreasing the chance for an earthquake). Just like for dynamic triggering, the fault needs to be about ready to slip. The effect on fault slip changes in “static coulomb stress” generally extend a distance of about 2-3 times the fault length of the “source” fault.
    • Raffie et al. (2021) calculated static coulomb stress changes on the central Sumatra fault segments as imposed by several megathrust subduction zone earthquakes.
    • Below is a series of maps that show the results from their analyses.

    • Coulomb stress models resolved on receiver faults of central part of GSF from coseismic slip model of each large interplate earthquakes. The color represents the maximum stress changes at 10 km depth with a scale saturated at 1 bar.

    • Here are the results when they consider all sources in a cumulative manner.

    • Cumulative ΔCFF of each earthquake listed in Table 1 (a) and cumulative ΔCFF of 1797, 1833, and 1861 earthquakes (b). The cyan ellipses are the damage area of large intraplate earthquakes marked as green star. The ΔCFF is calculated at 10 km depth with a scale saturated at 1 bar.

    Shaking Intensity

    • Here is a figure that shows a more detailed comparison between the modeled intensity and the reported intensity. Both data use the same color scale, the Modified Mercalli Intensity Scale (MMI). More about this can be found here. The colors and contours on the map are results from the USGS modeled intensity. The DYFI data are plotted as colored dots (color = MMI, diameter = number of reports).
    • In the upper panel is the USGS Did You Feel It reports map, showing reports as colored dots using the MMI color scale. Underlain on this map are colored areas showing the USGS modeled estimate for shaking intensity (MMI scale).
    • In the lower panel is a plot showing MMI intensity (vertical axis) relative to distance from the earthquake (horizontal axis). The models are represented by the green and orange lines. The DYFI data are plotted as light blue dots. The mean and median (different types of “average”) are plotted as orange and purple dots. Note how well the reports fit the green line (the model that represents how MMI works based on quakes in California).
    • Below the lower plot is the USGS MMI Intensity scale, which lists the level of damage for each level of intensity, along with approximate measures of how strongly the ground shakes at these intensities, showing levels in acceleration (Peak Ground Acceleration, PGA) and velocity (Peak Ground Velocity, PGV).

    Potential for Ground Failure

    • Below are a series of maps that show the potential for landslides and liquefaction. These are all USGS data products.
      There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:

      FOS = Resisting Force / Driving Force

    • When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.

    • Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
    • Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
    • An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
    • Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.

    • Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.

    • Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides.
    • I use the same color scheme that the USGS uses on their website. Note how the areas that are more likely to have experienced earthquake induced liquefaction are in the valleys. Learn more about how the USGS prepares these model results here.

    Seismic Hazard and Seismic Risk

    • Here is a map that shows the seismic hazard in southeast Asia, including Sumatra (Hayes et al., 2013). The plate convergence vectors, showing the direction of plate convergence and the rate of plate convergence in mm per year. Note how the plate convergence vectors are not perpendicular to the plate boundary.

    • These are the two maps shown in the map above, the GEM Seismic Hazard and the GEM Seismic Risk maps from Pagani et al. (2018) and Silva et al. (2018).
      • The GEM Seismic Hazard Map:


      • The Global Earthquake Model (GEM) Global Seismic Hazard Map (version 2018.1) depicts the geographic distribution of the Peak Ground Acceleration (PGA) with a 10% probability of being exceeded in 50 years, computed for reference rock conditions (shear wave velocity, VS30, of 760-800 m/s). The map was created by collating maps computed using national and regional probabilistic seismic hazard models developed by various institutions and projects, and by GEM Foundation scientists. The OpenQuake engine, an open-source seismic hazard and risk calculation software developed principally by the GEM Foundation, was used to calculate the hazard values. A smoothing methodology was applied to homogenise hazard values along the model borders. The map is based on a database of hazard models described using the OpenQuake engine data format (NRML). Due to possible model limitations, regions portrayed with low hazard may still experience potentially damaging earthquakes.
      • Here is a view of the GEM seismic hazard map for Indonesia.

      • The GEM Seismic Risk Map:


      • The Global Seismic Risk Map (v2018.1) presents the geographic distribution of average annual loss (USD) normalised by the average construction costs of the respective country (USD/m2) due to ground shaking in the residential, commercial and industrial building stock, considering contents, structural and non-structural components. The normalised metric allows a direct comparison of the risk between countries with widely different construction costs. It does not consider the effects of tsunamis, liquefaction, landslides, and fires following earthquakes. The loss estimates are from direct physical damage to buildings due to shaking, and thus damage to infrastructure or indirect losses due to business interruption are not included. The average annual losses are presented on a hexagonal grid, with a spacing of 0.30 x 0.34 decimal degrees (approximately 1,000 km2 at the equator). The average annual losses were computed using the event-based calculator of the OpenQuake engine, an open-source software for seismic hazard and risk analysis developed by the GEM Foundation. The seismic hazard, exposure and vulnerability models employed in these calculations were provided by national institutions, or developed within the scope of regional programs or bilateral collaborations.
    • Here is a view of the GEM seismic risk map for Indonesia.

    Tsunami Hazard

    • Here are two maps that show the results of probabilistic tsunami modeling for the nation of Indonesia (Horspool et al., 2014). These results are similar to results from seismic hazards analysis and maps. The color represents the chance that a given area will experience a certain size tsunami (or larger).
    • The first map shows the annual chance of a tsunami with a height of at least 0.5 m (1.5 feet). The second map shows the chance that there will be a tsunami at least 3 meters (10 feet) high at the coast.

    • Annual probability of experiencing a tsunami with a height at the coast of (a) 0.5m (a tsunami warning) and (b) 3m (a major tsunami warning).

      References:

      Basic & General References

    • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
    • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
    • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
    • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
    • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
    • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
    • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
    • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
    • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
    • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
    • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
    • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
    • Specific References

    • Abercrombie, R.E., Antolik, M., Ekstrom, G., 2003. The June 2000 Mw 7.9 earthquakes south of Sumatra: Deformation in the India–Australia Plate. Journal of Geophysical Research 108, 16.
    • Bassin, C., Laske, G. and Masters, G., The Current Limits of Resolution for Surface Wave Tomography in North America, EOS Trans AGU, 81, F897, 2000.
    • Bock, Y., Prawirodirdjo, L., Genrich, J.F., Stevens, C.W., McCaffrey, R., Subarya, C., Puntodewo, S.S.O., Calais, E., 2003. Crustal motion in Indonesia from Global Positioning System measurements: Journal of Geophysical Research, v. 108, no. B8, 2367, doi: 10.1029/2001JB000324.
    • Bothara, J., Beetham, R.D., Brunston, D., Stannard, M., Brown, R., Hyland, C., Lewis, W., Miller, S., Sanders, R., Sulistio, Y., 2010. General observations of effects of the 30th September 2009 Padang earthquake, Indonesia. Bulletin of the New Zealand Society for Earthquake Engineering 43, 143-173.
    • Chlieh, M., Avouac, J.-P., Hjorleifsdottir, V., Song, T.-R.A., Ji, C., Sieh, K., Sladen, A., Hebert, H., Prawirodirdjo, L., Bock, Y., Galetzka, J., 2007. Coseismic Slip and Afterslip of the Great (Mw 9.15) Sumatra-Andaman Earthquake of 2004. Bulletin of the Seismological Society of America 97, S152-S173.
    • Chlieh, M., Avouac, J.P., Sieh, K., Natawidjaja, D.H., Galetzka, J., 2008. Heterogeneous coupling of the Sumatran megathrust constrained by geodetic and paleogeodetic measurements: Journal of Geophysical Research, v. 113, B05305, doi: 10.1029/2007JB004981.
    • DEPLUS, C. et al., 1998 – Direct evidence of active derormation in the eastern Indian oceanic plate, Geology.
    • DYMENT, J., CANDE, S.C. & SINGH, S., 2007 – Oceanic lithosphere subducting beneath the Sunda Trench: the Wharton Basin revisited. European Geosciences Union General Assembly, Vienna, 15-20/05.
    • Hayes, G. P., Wald, D. J., and Johnson, R. L., 2012. Slab1.0: A three-dimensional model of global subduction zone geometries in J. Geophys. Res., 117, B01302, doi:10.1029/2011JB008524.
    • Hayes, G.P., Bernardino, Melissa, Dannemann, Fransiska, Smoczyk, Gregory, Briggs, Richard, Benz, H.M., Furlong, K.P., and Villaseñor, Antonio, 2013. Seismicity of the Earth 1900–2012 Sumatra and vicinity: U.S. Geological Survey Open-File Report 2010–1083-L, scale 1:6,000,000, https://pubs.usgs.gov/of/2010/1083/l/.
    • JACOB, J., DYMENT, J., YATHEESH, V. & BHATTACHARYA, G.C., 2009 – Marine magnetic anomalies in the NE Indian Ocean: the Wharton and Central Indian basins revisited. European Geosciences Union General Assembly, Vienna, 19-24/04.
    • Ji, C., D.J. Wald, and D.V. Helmberger, Source description of the 1999 Hector Mine, California earthquake; Part I: Wavelet domain inversion theory and resolution analysis, Bull. Seism. Soc. Am., Vol 92, No. 4. pp. 1192-1207, 2002.
    • Hurukawa, N., Wulandari, B.R., and Kasahara, M., 2014. Earthquake History of the Sumatran Fault, Indonesia, since 1892, Derived from Relocation of Large Earthquakes in BSSA,v. 104, no. 4, p. 1750-1762, doi: 10.1785/0120130201
    • Ishii, M., Shearer, P.M., Houston, H., Vidale, J.E., 2005. Extent, duration and speed of the 2004 Sumatra-Andaman earthquake imaged by the Hi-Net array. Nature 435, 933.
    • Kanamori, H., Rivera, L., Lee, W.H.K., 2010. Historical seismograms for unravelling a mysterious earthquake: The 1907 Sumatra Earthquake. Geophysical Journal International 183, 358-374.
    • Konca, A.O., Avouac, J., Sladen, A., Meltzner, A.J., Sieh, K., Fang, P., Li, Z., Galetzka, J., Genrich, J., Chlieh, M., Natawidjaja, D.H., Bock, Y., Fielding, E.J., Ji, C., Helmberger, D., 2008. Partial Rupture of a Locked Patch of the Sumatra Megathrust During the 2007 Earthquake Sequence. Nature 456, 631-635.
    • Lange, D., Cembrano, J., Rietbrock, A., Haberland, C., Dahm, T., and Bataille, K., 2008. First seismic record for intra-arc strike-slip tectonics along the Liquiñe-Ofqui fault zone at the obliquely convergent plate margin of the southern Andes in Tectonophysics, v. 455, p. 14-24
    • Maus, S., et al., 2009. EMAG2: A 2–arc min resolution Earth Magnetic Anomaly Grid compiled from satellite, airborne, and marine magnetic measurements, Geochem. Geophys. Geosyst., 10, Q08005, doi:10.1029/2009GC002471.
    • Malik, J.N., Shishikura, M., Echigo, T., Ikeda, Y., Satake, K., Kayanne, H., Sawai, Y., Murty, C.V.R., Dikshit, D., 2011. Geologic evidence for two pre-2004 earthquakes during recent centuries near Port Blair, South Andaman Island, India: Geology, v. 39, p. 559-562.
    • Meltzner, A.J., Sieh, K., Chiang, H., Shen, C., Suwargadi, B.W., Natawidjaja, D.H., Philobosian, B., Briggs, R.W., Galetzka, J., 2010. Coral evidence for earthquake recurrence and an A.D. 1390–1455 cluster at the south end of the 2004 Aceh–Andaman rupture. Journal of Geophysical Research 115, 1-46.
    • Meng, L., Ampuero, J.-P., Stock, J., Duputel, Z., Luo, Y., and Tsai, V.C., 2012. Earthquake in a Maze: Compressional Rupture Branching During the 2012 Mw 8.6 Sumatra Earthquake in Science, v. 337, p. 724-726.
    • Natawidjaja, D.H., Sieh, K., Chlieh, M., Galetzka, J., Suwargadi, B., Cheng, H., Edwards, R.L., Avouac, J., Ward, S.N., 2006. Source parameters of the great Sumatran megathrust earthquakes of 1797 and 1833 inferred from coral microatolls. Journal of Geophysical Research 111, 37.
    • Natawidjaja, D.H., 2018. Updating active fault maps and slip rates along the Sumatran Fault Zone, Indonesia in IOP Conf. Ser.: Earth Environ. Sci. 118 012001
    • Natawidjaja, D.H., 2002. Neotectonics of the Sumatran Fault And Paleogeodesy of the Sumatran Subduction Zone [thesis] Californiat Institute of Technilogy, Pasadena, California, 289 pp.
    • Newcomb, K.R., McCann, W.R., 1987. Seismic History and Seismotectonics of the Sunda Arc. Journal of Geophysical Research 92, 421-439.
    • Philibosian, B., Sieh, K., Natawidjaja, D.H., Chiang, H., Shen, C., Suwargadi, B., Hill, E.M., Edwards, R.L., 2012. An ancient shallow slip event on the Mentawai segment of the Sunda megathrust, Sumatra. Journal of Geophysical Research 117, 12.
    • Prawirodirdjo, P., McCaffrey,R., Chadwell, D., Bock, Y, and Subarya, C., 2010. Geodetic observations of an earthquake cycle at the Sumatra subduction zone: Role of interseismic strain segmentation, JOURNAL OF GEOPHYSICAL RESEARCH, v. 115, B03414, doi:10.1029/2008JB006139
    • Rafie, M.T., Cummins, P.R., Sahara, D.P., Widiyantoro, S., Triyoso, W., and Nugraha, A.D., 2021. Preliminary study of stress changes evolution on central part of Sumatran Fault influenced by large interplate earthquakes and tectonic stress rates, IOP Conf. Ser.: Earth Environ. Sci. 873 012011
    • Rivera, L., Sieh, K., Helmberger, D., Natawidjaja, D.H., 2002. A Comparative Study of the Sumatran Subduction-Zone Earthquakes of 1935 and 1984. BSSA 92, 1721-1736.
    • Shearer, P., and Burgmann, R., 2010. Lessons Learned from the 2004 Sumatra-Andaman Megathrust Rupture, Annu. Rev. Earth Planet. Sci. v. 38, pp. 103–31
    • SATISH C. S, CARTON H, CHAUHAN A.S., et al., 2011 – Extremely thin crust in the Indian Ocean possibly resulting from Plume-Ridge Interaction, Geophysical Journal International.
    • Sieh, K., Natawidjaja, D.H., Meltzner, A.J., Shen, C., Cheng, H., Li, K., Suwargadi, B.W., Galetzka, J., Philobosian, B., Edwards, R.L., 2008. Earthquake Supercycles Inferred from Sea-Level Changes Recorded in the Corals of West Sumatra. Science 322, 1674-1678.
    • Singh, S.C., Carton, H.L., Tapponnier, P, Hananto, N.D., Chauhan, A.P.S., Hartoyo, D., Bayly, M., Moeljopranoto, S., Bunting, T., Christie, P., Lubis, H., and Martin, J., 2008. Seismic evidence for broken oceanic crust in the 2004 Sumatra earthquake epicentral region, Nature Geoscience, v. 1, pp. 5.
    • Smith, W.H.F., Sandwell, D.T., 1997. Global seafloor topography from satellite altimetry and ship depth soundings: Science, v. 277, p. 1,957-1,962.
    • Sorensen, M.B., Atakan, K., Pulido, N., 2007. Simulated Strong Ground Motions for the Great M 9.3 Sumatra–Andaman Earthquake of 26 December 2004. BSSA 97, S139-S151.
    • Subarya, C., Chlieh, M., Prawirodirdjo, L., Avouac, J., Bock, Y., Sieh, K., Meltzner, A.J., Natawidjaja, D.H., McCaffrey, R., 2006. Plate-boundary deformation associated with the great Sumatra–Andaman earthquake: Nature, v. 440, p. 46-51.
    • Tolstoy, M., Bohnenstiehl, D.R., 2006. Hydroacoustic contributions to understanding the December 26th 2004 great Sumatra–Andaman Earthquake. Survey of Geophysics 27, 633-646.
    • Zhu, Lupei, and Donald V. Helmberger. “Advancement in source estimation techniques using broadband regional seismograms.” Bulletin of the Seismological Society of America 86.5 (1996): 1634-1641.

    Return to the Earthquake Reports page.


    Earthquake Report: M 5.7 & 6.2 Mendocino triple junction

    I was returning from New Orleans where I was attending the American Geophysical Union Fall Meeting. There was a short layover in Denver and I had a short time to find some food, which is challenging with my dietary restrictions. I cannot recall precisely, but I got some notification from my CGS crew about a magnitude M 6.2 earthquake offshore of the Mendocino triple junction. One of these notifications was from Cindy as we both collaborate to prepare quick reports for earthquakes in California. These reports are sent upstream to management in our organization and others. I was unavailable to contribute this time.

    Needless to say, I was sad to have missed experiencing this good sized shaker for myself. This is the first earthquake of this size that I have missed (in Humboldt) since I moved here in 1991.
    I got home about 3 am the next morning and did not have energy to prepare an earthjay report. Though I started working on it the next day. However, I soon learned that this was a complicated earthquake and I decided to await additional analyses by the Berkeley Seismmo Lab and the USGS.

    Last week or so, their analyses were produced publicly and the earthquake catalog was updated. What we discovered is that there were two closely spaced (in time but not space) earthquakes, an M 5.7 and and M 6.2.

    https://earthquake.usgs.gov/earthquakes/eventpage/nc71127029/executive
    https://earthquake.usgs.gov/earthquakes/eventpage/nc73666231/executive

    It was complicated for the seismologists to work out because the seismic waves of the two events overlapped in time. i.e., the waves from the first quake were still passing through the Earth when the waves from the second quake started.

    Basically, there was initially an M 5.7 strike-slip earthquake along the Mendocino transform fault zone about 20 km (12.5 miles) offshore. About 10 or 11 seconds later, there was an M 6.2 strike-slip earthquake within the Gorda plate, below the megathrust fault.

    Here is a plot from the USGS. Each horizontal squiggly line is the seismograph record from an individual seismometer. They are plotted with the seismometer closest to the earthquake on the bottom row and the furthest seismometer on the uppermost row.

    The P wave (primary wave) is the first of four major types of seismic waves. Next comes the S (secondary) wave, then the Love waves, and finally the Raleigh waves.

    The P wave arrives at closer seismometers before it arrives at more distant seismometers. Because of this, we generally call this type of plot a travel time plot.


    In the above plot we can see how the M 6.1 P waves are arriving while the M 5.7 S waves are still being transmitted.

    The M 5.7 is clearly a right-lateral strike-slip event given the aftershock pattern and the known location and type of the Mendocino fault system (a right-lateral strike-slip fault.
    The interpretation for the type of earthquake for the M 6.2 is a little more complicated.

    Earthquake mechanisms (the “beach balls”) show two possible ways that the earthquake could have slipped. We use aftershock patterns and existing mapped faults to help us interpret which of these [nodal] fault planes is the more likely one.

    If we look at the earthquake poster below, we see that the M 6.2 earthquake is an almost pure strike-slip earthquake. The two possible fault planes are one that is oriented in the northwest direction (would be right-lateral) and one that is in the northeast direction (would be left-lateral).

      There are two reasons why I interpret the M 6.2 to be right-lateral (of course, I could be wrong).

    1. There are a series of aftershocks that appear to align along a northwest trajectory. See the 6 mechanisms for the earthquakes just north of the M 6.2. The epicenters for these M 3.8, 3.0, 3.2, 3.0, 3.3, and 3.8 earthquakes align with a northwest oriented trend.
    2. The 1992 Cape Mendocino Earthquake main events began with an M 7.2 thrust type earthquake mainshock followed by two triggered events with magnitudes of M 6.6 and M 6.5. These triggered events are strike-slip and within the Gorda plate (“intraplate” events). If we take a look at the earthquake catalog after it has been modeled using a “double differencing” technique (read more here), we would notice several northwest trends in seismicity. These trends appeared after the 1992 sequence.

    So, while most of our experience with the Gorda plate is with northeast oriented (striking) left-lateral strike-slip faults (e.g., 1980, 2010, 2014, etc.) it is possible that there are other faults, sub-parallel to the post-1992 seismicity trends, where the M 6.2 and other aftershocks were hosted. I mention these northwest trending faults in a recent Earthquake Report here.

    Something that is interesting is that the onshore events from this 20 Dec 2021 sequence are just to the north of the aftershocks from the 1992 sequence. They are at similar depths as those ’92 quakes and have similar earthquake mechanisms. As Spock would say, Fascinating.

    Dr. Anthony Lomax, famous for his work locating the hypocenter for the 1906 San Francisco Earthquake, has been developing excellent tools for seismologists ever since. He recently applied one of his new tools to locate earthquakes to the Mendocino triple junction region. I present some of his figures below.

    Below is my interpretive poster for this earthquake

    • I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1920-2020 with magnitudes M ≥ 3.0 in one version.
    • I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
    • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
    • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.

      I include some inset figures. Some of the same figures are located in different places on the larger scale map below.

    • In the upper left corner is a small scale view of the Mendocino triple junction where the southern Cascadia subduction zone meets the right-lateral Mendocino and San Andreas strike-slip faults.
    • In the lower center is a plot of earthquake intensity (vertical axis) relative to distance from the earthquake (horizontal axis). The blue and orange dots represent USGS Did You Feel It? reports, observations from real people. The green and orange lines show the plots from the USGS [empirical] models of shaking intensity.
    • To the left of the plot is a map that shows these same data. The colored areas are the average intensity reported by people (dyfi numbers). The colored lines represent the USGS modeled intensities. Both the plot and the map show that the shaking intensity gets smaller with distance from the earthquake.
    • In the left upper center is a map that shows liquefaction susceptibility. This is a model from the USGS that uses ground shaking data to estimate where there may be liquefaction. I drove around after the earthquake and could not locate any evidence for liquefaction but my search was far from comprehensive.
    • Here is the map with a month’s seismicity plotted.

    • This one helps us compare the mainshock and two main triggered earthquakes.

    • Here is a poster that shows a comparison between the 1991 Honeydew and 1992 Cape Mendocino earthquakes..

    Other Report Pages

    Shaking Intensity and Potential for Ground Failure

    • Below are a series of maps that show the shaking intensity and potential for landslides and liquefaction. These are all USGS data products.
    • There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:

      FOS = Resisting Force / Driving Force

      When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.


      Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.

      Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.

      An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.

      Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.


      Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.


    • Below is the liquefaction susceptibility and landslide probability map (Jessee et al., 2017; Zhu et al., 2017). Please head over to that report for more information about the USGS Ground Failure products (landslides and liquefaction). Basically, earthquakes shake the ground and this ground shaking can cause landslides. We can see that there is a low probability for landslides. However, we have already seen photographic evidence for landslides and the lower limit for earthquake triggered landslides is magnitude M 5.5 (from Keefer 1984)
    • I use the same color scheme that the USGS uses on their website. Note how the areas that are more likely to have experienced earthquake induced liquefaction are in the valleys. Learn more about how the USGS prepares these model results here.

    Some Relevant Discussion and Figures

    • Here is a map of the Cascadia subduction zone, modified from Nelson et al. (2006). The Juan de Fuca and Gorda plates subduct norteastwardly beneath the North America plate at rates ranging from 29- to 45-mm/yr. Sites where evidence of past earthquakes (paleoseismology) are denoted by white dots. Where there is also evidence for past CSZ tsunami, there are black dots. These paleoseismology sites are labeled (e.g. Humboldt Bay). Some submarine paleoseismology core sites are also shown as grey dots. The two main spreading ridges are not labeled, but the northern one is the Juan de Fuca ridge (where oceanic crust is formed for the Juan de Fuca plate) and the southern one is the Gorda rise (where the oceanic crust is formed for the Gorda plate).

    • Here is a version of the CSZ cross section alone (Plafker, 1972). This shows two parts of the earthquake cycle: the interseismic part (between earthquakes) and the coseismic part (during earthquakes). Regions that experience uplift during the interseismic period tend to experience subsidence during the coseismic period.

    • This figure shows how a subduction zone deforms between (interseismic) and during (coseismic) earthquakes. We also can see how a subduction zone generates a tsunami. Atwater et al., 2005.

    • Here is an animation produced by the folks at Cal Tech following the 2004 Sumatra-Andaman subduction zone earthquake. I have several posts about that earthquake here and here. One may learn more about this animation, as well as download this animation here.

    The Gorda and Juan de Fuca plates subduct beneath the North America plate to form the Cascadia subduction zone fault system. In 1992 there was a swarm of earthquakes with the magnitude Mw 7.2 Mainshock on 4/25. Initially this earthquake was interpreted to have been on the Cascadia subduction zone (CSZ). The moment tensor shows a compressional mechanism. However the two largest aftershocks on 4/26/1992 (Mw 6.5 and Mw 6.7), had strike-slip moment tensors. In my mind, these two aftershocks aligned on what may be the eastern extension of the Mendocino fault. However, looking at their locations, my mind was incorrect. These two earthquakes were not aftershocks, but were either left-lateral or right-lateral strike-slip Gorda plate earthquakes triggered by the M 7.1 thrust event.

    These two quakes appear to be aligned with the two northwest trends in seismicity and the 18 March 2020 M 5.2. The orientation of the mechanisms are not as perfectly well aligned, but there are lots of reasons for this (perhaps the faults were formed in a slightly different orientation, but have rotated slightly).

    There have been several series of intra-plate earthquakes in the Gorda plate. Two main shocks that I plot of this type of earthquake are the 1980 (Mw 7.2) and 2005 (Mw 7.2) earthquakes. I place orange lines approximately where the faults are that ruptured in 1980 and 2005. These are also plotted in the Rollins and Stein (2010) figure above. The Gorda plate is being deformed due to compression between the Pacific plate to the south and the Juan de Fuca plate to the north. Due to this north-south compression, the plate is deforming internally so that normal faults that formed at the spreading center (the Gorda Rise) are reactivated as left-lateral strike-slip faults. In 2014, there was another swarm of left-lateral earthquakes in the Gorda plate. I posted some material about the Gorda plate setting on this page.

    • This is the map used in the animation below. Earthquake epicenters are plotted (some with USGS moment tensors) for this region from 1917-2017 with M ≥ 6.5. I labeled the plates and shaded their general location in different colors.
    • I include some inset maps.
      • In the upper right corner is a map of the Cascadia subduction zone (Chaytor et al., 2004; Nelson et al., 2004).
      • In the upper left corner is a map from Rollins and Stein (2010). They plot epicenters and fault lines involved in earthquakes between 1976 and 2010.


      • Here is a map from Rollins and Stein, showing their interpretations of different historic earthquakes in the region. This was published in response to the Januray 2010 Gorda plate earthquake. The faults are from Chaytor et al. (2004).

      • Tectonic configuration of the Gorda deformation zone and locations and source models for 1976–2010 M ≥ 5.9 earthquakes. Letters designate chronological order of earthquakes (Table 1 and Appendix A). Plate motion vectors relative to the Pacific Plate (gray arrows in main diagram) are from Wilson [1989], with Cande and Kent’s [1995] timescale correction.

      • Here is a large scale map of the 1994 earthquake swarm. The mainshock epicenter is a black star and epicenters are denoted as white circles.

      • Here is a plot of focal mechanisms from the Dengler et al. (1995) paper in California Geology.

        • In this map below, I label a number of other significant earthquakes in this Mendocino triple junction region. Another historic right-lateral earthquake on the Mendocino fault system was in 1994. There was a series of earthquakes possibly along the easternmost section of the Mendocino fault system in late January 2015, here is my post about that earthquake series.

        • Here is a map from Chaytor et al. (2004) that shows some details of the faulting in the region. The moment tensor (at the moment i write this) shows a north-south striking fault with a reverse or thrust faulting mechanism. While this region of faulting is dominated by strike slip faults (and most all prior earthquake moment tensors showed strike slip earthquakes), when strike slip faults bend, they can create compression (transpression) and extension (transtension). This transpressive or transtentional deformation may produce thrust/reverse earthquakes or normal fault earthquakes, respectively. The transverse ranges north of Los Angeles are an example of uplift/transpression due to the bend in the San Andreas fault in that region.

        • A: Mapped faults and fault-related ridges within Gorda plate based on basement structure and surface morphology, overlain on bathymetric contours (gray lines—250 m interval). Approximate boundaries of three structural segments are also shown. Black arrows indicated approximate location of possible northwest- trending large-scale folds. B, C: uninterpreted and interpreted enlargements of center of plate showing location of interpreted second-generation strike-slip faults and features that they appear to offset. OSC—overlapping spreading center.

        • These are the models for tectonic deformation within the Gorda plate as presented by Jason Chaytor in 2004.
        • Mw = 5 Trinidad Chaytor

          Models of brittle deformation for Gorda plate overlain on magnetic anomalies modified from Raff and Mason (1961). Models A–F were proposed prior to collection and analysis of full-plate multibeam data. Deformation model of Gulick et al. (2001) is included in model A. Model G represents modification of Stoddard’s (1987) flexural-slip model proposed in this paper.


      Social Media

      References:

      Basic & General References

    • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
    • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
    • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
    • Jessee, M.A.N., Hamburger, M. W., Allstadt, K., Wald, D. J., Robeson, S. M., Tanyas, H., et al. (2018). A global empirical model for near-real-time assessment of seismically induced landslides. Journal of Geophysical Research: Earth Surface, 123, 1835–1859. https://doi.org/10.1029/2017JF004494
    • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
    • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
    • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
    • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
    • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
    • Pagani,M. , J. Garcia-Pelaez, R. Gee, K. Johnson, V. Poggi, R. Styron, G. Weatherill, M. Simionato, D. Viganò, L. Danciu, D. Monelli (2018). Global Earthquake Model (GEM) Seismic Hazard Map (version 2018.1 – December 2018), DOI: 10.13117/GEM-GLOBAL-SEISMIC-HAZARD-MAP-2018.1
    • Silva, V ., D Amo-Oduro, A Calderon, J Dabbeek, V Despotaki, L Martins, A Rao, M Simionato, D Viganò, C Yepes, A Acevedo, N Horspool, H Crowley, K Jaiswal, M Journeay, M Pittore, 2018. Global Earthquake Model (GEM) Seismic Risk Map (version 2018.1). https://doi.org/10.13117/GEM-GLOBAL-SEISMIC-RISK-MAP-2018.1
    • Zhu, J., Baise, L. G., Thompson, E. M., 2017, An Updated Geospatial Liquefaction Model for Global Application, Bulletin of the Seismological Society of America, 107, p 1365-1385, https://doi.org/0.1785/0120160198
    • Specific References

    • Atwater, B.F., Musumi-Rokkaku, S., Satake, K., Tsuju, Y., Eueda, K., and Yamaguchi, D.K., 2005. The Orphan Tsunami of 1700—Japanese Clues to a Parent Earthquake in North America, USGS Professional Paper 1707, USGS, Reston, VA, 144 pp.
    • Chaytor, J.D., Goldfinger, C., Dziak, R.P., and Fox, C.G., 2004. Active deformation of the Gorda plate: Constraining deformation models with new geophysical data: Geology v. 32, p. 353-356.
    • Dengler, L.A., Moley, K.M., McPherson, R.C., Pasyanos, M., Dewey, J.W., and Murray, M., 1995. The September 1, 1994 Mendocino Fault Earthquake, California Geology, Marc/April 1995, p. 43-53.
    • Geist, E.L. and Andrews D.J., 2000. Slip rates on San Francisco Bay area faults from anelastic deformation of the continental lithosphere, Journal of Geophysical Research, v. 105, no. B11, p. 25,543-25,552.
    • Irwin, W.P., 1990. Quaternary deformation, in Wallace, R.E. (ed.), 1990, The San Andreas Fault system, California: U.S. Geological Survey Professional Paper 1515, online at: http://pubs.usgs.gov/pp/1990/1515/
    • McCrory, P.A.,. Blair, J.L., Waldhauser, F., kand Oppenheimer, D.H., 2012. Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity in JGR, v. 117, B09306, doi:10.1029/2012JB009407.
    • McLaughlin, R.J., Sarna-Wojcicki, A.M., Wagner, D.L., Fleck, R.J., Langenheim, V.E., Jachens, R.C., Clahan, K., and Allen, J.R., 2012. Evolution of the Rodgers Creek–Maacama right-lateral fault system and associated basins east of the northward-migrating Mendocino Triple Junction, northern California in Geosphere, v. 8, no. 2., p. 342-373.
    • Nelson, A.R., Asquith, A.C., and Grant, W.C., 2004. Great Earthquakes and Tsunamis of the Past 2000 Years at the Salmon River Estuary, Central Oregon Coast, USA: Bulletin of the Seismological Society of America, Vol. 94, No. 4, pp. 1276–1292
    • Rollins, J.C. and Stein, R.S., 2010. Coulomb stress interactions among M ≥ 5.9 earthquakes in the Gorda deformation zone and on the Mendocino Fault Zone, Cascadia subduction zone, and northern San Andreas Fault: Journal of Geophysical Research, v. 115, B12306, doi:10.1029/2009JB007117, 2010.
    • Stoffer, P.W., 2006, Where’s the San Andreas Fault? A guidebook to tracing the fault on public lands in the San Francisco Bay region: U.S. Geological Survey General Interest Publication 16, 123 p., online at http://pubs.usgs.gov/gip/2006/16/
    • Wallace, Robert E., ed., 1990, The San Andreas fault system, California: U.S. Geological Survey Professional Paper 1515, 283 p. [http://pubs.usgs.gov/pp/1988/1434/].
    • Wells, D.L., and Coopersmith, K.J., 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement in BSSA, v. 84, no. 4, p. 974-1002

    Return to the Earthquake Reports page.