Earthquake Report: 1992.04.25 M 7.1 Petrolia

The 25 April 1992 M 7.1 earthquake was a wake up call for many, like all large magnitude earthquakes are.

    Here is the USGS website for these three large earthquakes.

  • 1992-04-25 18:06:05 UTC 40.335°N 124.229°W 9.9 km depth M 7.2
  • 1992-04-26 07:41:40 UTC 40.433°N 124.566°W 18.8 km depth M 6.5
  • 1992-04-26 11:18:25 UTC 40.383°N 124.555°W 21.7 km depth M 6.6

Below is my interpretive poster for this earthquake.

I plot the seismicity for a week beginning April 25, 1992, with color representing depth and diameter representing magnitude (see legend)..

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely.
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (McCrory et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault.

    I include some inset figures in the poster.

  • In the upper left corner is a map of the Cascadia subduction zone (CSZ) and regional tectonic plate boundary faults. This is modified from several sources (Chaytor et al., 2004; Nelson et al., 2004)
  • Below the CSZ map is an illustration modified from Plafker (1972). This figure shows how a subduction zone deforms between (interseismic) and during (coseismic) earthquakes.
  • In the upper left corner is a figure from Rollins and Stein (2010). In their paper they discuss how static coulomb stress changes from earthquakes may impart (or remove) stress from adjacent crust/faults. To the right of this map are two panels. The upper panel shows the location and orientation of the fault plane used by Rollins and Stein (2010) to model potential changes in coulomb stress following the 1992 M 7.2 earthquake. The Lower panel shows the results from this modeling.
  • In the lower right corner is the map from Stein et al. (1993). This map shows an estimate of coseismic vertical ground motion induced by the 1992 earthquake sequence.
  • In the upper right corner is a series of USGS shakemaps. These plot intensity using the MMI scale.
  • Below the shakemaps is the “Did You Feel It?” map and attenuation relation plot.


  • 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).

  • This figure shows how a subduction zone deforms between (interseismic) and during (coseismic) earthquakes.

  • 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.
  • Following the earthquake, there was lots of work done by local geologists, along with help from those visiting from out of the area. One of the projects included the measurement and modeling of the ground deformation related to the earthquake. The measurements consistend of a first order survey of benchmarks, along with Global Positioning System measurements at GPS monuments. The results from these analyses were presented in a U.S. Geological Survey Open-File Report 93-383 (Stein et al., 1993). Below is a map that shows a modeled estimate of the surface deformation associated with this earthquake.

  • Here is a figure from Oppenheimer et al. (1993) that shows the shaking intensity from this earthquake sequence. Below is a colorized version.


  • Simplified tectonic map in the vicinity of the Cape Mendocino earthquake sequence. Stars, epicenters of three largest earthquakes; contours, Modified Mercalli intensities (values, Roman numerals) of main shock; open circles, strong motion instrument sites (adjacent numbers give peak horizontal accelerations in g). Abbreviations FT Fortuna; F Ferndale; RD, Rio Dell; S, Scotia; P, Petrolia; H, Honeydew; MF, Mendocino fault; CSZ, seaward edge of Cascadia subduction zone; and SAF, San Andreas fault.

  • This map shows an alternate model of earthquake ground deformation (Oppenheimer et al, 1993).

  • Observed and predicted coseismic displacements for the Cape Mendocino main shock (epicenter located at star).

  • This is a figure that shows the tsunami recorded by tide gages in California, Hawaii, and Oregon (Oppenheimer et al., 1993)

  • Here is a map from Rollins and Stein (2010), 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.

  • This figure shows the fault plane and aftershocks used in their analysis of the 1992 earthquake sequence.

  • Source models for earthquakes 25 April 1992, Mw = 6.9, open circles are from Waldhauser and Schaff ’s [2008] earthquake locations for 25 April 1992 (1806 UTC) to 26 April 1992 (0741 UTC)

  • This figure shows the change in coulomb stress imparted by the M 7.1 earthquake onto different faults: (a) the CSZ and (b) the faults that were triggered to generate the two main aftershocks.

  • (a) Coulomb stress changes imparted by the 1992 Mw = 6.9 Cape Mendocino earthquake (J) to the Cascadia subduction zone. Calculation depth is 8 km. Open circles are Waldhauser and Schaff [2008] earthquake locations for 25 April 1992 to 2 May 1992, 0–15 km depth. Seismicity data were cut off at 15 km depth to prevent interference from aftershocks of K and L. Cross section A‐A′ includes seismicity between 40.24°N and 40.36°N. Cross section B‐B′ includes seismicity between 40.36°N and 40.48°N. (b) Coulomb stress changes imparted by the 1992 Mw = 6.9 earthquake (J) to Mw = 6.5 and Mw = 6.6 shocks the next day (K and L). Stress change is resolved on the average of the orientations of K and L (strike 127°/dip 90°/rake 180°). Calculation depth is 21.5 km. (c) Calculated Coulomb stress changes imparted by M ≥ 5.9 shocks in 1983, 1987, and 1992 (C, E, and J) to the epicenters of K and L. The series of three colored numbers represent stress changes imparted by C, E, and J, respectively.

  • Here is a plot of the seismograms from the NCEDC.

    Here is the USGS website for all the earthquakes in this region from 1917-2017 with M ≥ 6.5.

  • 1922.01.31 13:17 M 7.3
  • 1923.01.22 09:04 M 6.9
  • 1934-07-06 22:48 M 6.7
  • 1941-02-09 09:44 M 6.8
  • 1949-03-24 20:56 M 6.5
  • 1954-11-25 11:16 M 6.8
  • 1954-12-21 19:56 M 6.6
  • 1980-11-08 10:27 M 7.2
  • 1984-09-10 03:14 M 6.7
  • 1984-09-10 03:14 M 6.6
  • 1991-07-13 02:50 M 6.9
  • 1991-08-17 22:17 M 7.0
  • 1992-04-25 18:06 M 7.2
  • 1992-04-26 07:41 M 6.5
  • 1992-04-26 11:18 M 6.6
  • 1994-09-01 15:15 M 7.0
  • 1995-02-19 04:03 M 6.6
  • 2005-06-15 02:50 M 7.2
  • 2005-06-17 06:21 M 6.6
  • 2010-01-10 00:27 M 6.5
  • 2014-03-10 05:18 M 6.8
  • 2016-12-08 14:49 M 6.5
  • 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.



  • There are three types of earthquakes, strike-slip, compressional (reverse or thrust, depending upon the dip of the fault), and extensional (normal). Here is are some animations of these three types of earthquake faults. Many of the earthquakes people are familiar with in the Mendocino triple junction region are either compressional or strike slip. The following three animations are from IRIS.
  • Strike Slip:
  • Compressional:
  • Extensional:
  • Here is a primer that helps people learn how to interpret focal mechanisms and moment tensors. Moment tensors are calculated differently from focal mechanisms, but the interpretation of their graphical solution is similar. This is from the USGS.

  • For more on the graphical representation of moment tensors and focal mechnisms, check this IRIS video out:

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.
  • Goldfinger, C., Nelson, C.H., Morey, A., Johnson, J.E., Gutierrez-Pastor, J., Eriksson, A.T., Karabanov, E., Patton, J., Gràcia, E., Enkin, R., Dallimore, A., Dunhill, G., and Vallier, T., 2012 a. Turbidite Event History: Methods and Implications for Holocene Paleoseismicity of the Cascadia Subduction Zone, USGS Professional Paper # 1661F. U.S. Geological Survey, Reston, VA, 184 pp.
  • McCrory, P.A., 2000, Upper plate contraction north of the migrating Mendocino triple junction, northern California: Implications for partitioning of strain: Tectonics, v. 19, p. 11441160.
  • McCrory, P. A., Blair, J. L., Oppenheimer, D. H., and Walter, S. R., 2006, Depth to the Juan de Fuca slab beneath the Cascadia subduction margin; a 3-D model for sorting earthquakes U. S. Geological Survey
  • Nelson, A.R., Kelsey, H.M., Witter, R.C., 2006. Great earthquakes of variable magnitude at the Cascadia subduction zone. Quaternary Research 65, 354-365.
  • Oppenheimer, D., Beroza, G., Carver, G., Dengler, L., Eaton, J., Gee, L., Gonzalez, F., Jayko, A., Ki., W.H., Lisowski, M., Magee, M., Marshall, G., Murray, M., McPherson, R., Romanowicz, B., Satake, K., Simpson, R., Somerille, P., Stein, R., and Valentine, D., The Cape Mendocino, California, Earthquakes of April, 1992: Subduction at the Triple Junction in Science, v. 261, no. 5120, p. 433-438.
  • Patton, J. R., Goldfinger, C., Morey, A. E., Romsos, C., Black, B., Djadjadihardja, Y., and Udrekh, 2013. Seismoturbidite record as preserved at core sites at the Cascadia and Sumatra–Andaman subduction zones, Nat. Hazards Earth Syst. Sci., 13, 833-867, doi:10.5194/nhess-13-833-2013, 2013.
  • Plafker, G., 1972. Alaskan earthquake of 1964 and Chilean earthquake of 1960: Implications for arc tectonics in Journal of Geophysical Research, v. 77, p. 901-925.
  • 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.
  • Stein, R.S., Marshall, G.A., Murray, M.H., Balazs, E., Carver, G.A., Dunklin, T.A>, McLaughlin, R.J., Cyr, K., and Jayko, A., 1993. Permanent Ground Movement Associate with the 1992 M=7 Cape Mendocino, California, Earthquake: Implications for Damage to Infrastructure and Hazards to navigation, U.S. Geological Survey Open-File Report 93-383.
  • Wang, K., Wells, R., Mazzotti, S., Hyndman, R. D., and Sagiya, T., 2003, A revised dislocation model of interseismic deformation of the Cascadia subduction zone Journal of Geophysical Research, B, Solid Earth and Planets v. 108, no. 1.

Earthquake Report: Chile Update #1

Well, I thought more to compare this ongoing earthquake sequence with the 1985 M 8.0 earthquake. This, in context with the 2010 and 2015 earthquakes. My initial report based upon the M ~4-5.9 swarm is here and my report on the “current” M 6.9 mainshock is here. More information about the background for the tectonics along the plate boundary, please refer to those earlier reports.

I used the USGS epicenters for earthquakes with magnitudes M ≥ 2.5. For each earthquake (1985, 2010, and 2015) I chose a month of seismicity beginning 3 days before the mainshock. Then I digitized the general outline of the earthquakes. This is a rough approximation for the slip patch for each of these earthquakes. I separated the interface earthquakes from the triggered outer rise earthquakes into separate polygons for the 2010 and 2015 earthquakes (they both appear to have triggered earthquakes in the downgoing Nazca plate to the west of the subduction zone fault, where it flexes in response to subduction here.

This current sequence is about the same magnitude and along-strike size as the 1971 earthquake (M 7.0). This sequence also lies within the 1985 earthquake aftershock region (and also within the northernmost area of the 2010 aftershock region). The M 6.9 could still be a foreshock of a larger earthquake. The 1985 earthquake was preceded by 3 earthquakes in the M 4-5.5 range. But, looking into the past, there are instances when this part of the fault only ruptures a small patch (1971, 1873, 1851). Given that this part of the fault slipped recently (2010), it seems more probable that there won’t be a larger earthquake (M > 8.0). This is difficult to know because we don’t really know the state of stress on the fault (how ready it is to rupture in an earthquake). I still cannot stop thinking about the Juan Fernandez Ridge and how this plays a part in this story.

I include the moment tensors from each of the Great Earthquakes, as well as the 2017 M 6.9 earthquake.

  • In the interpretive poster below
    • I outline the 1985 aftershock region in black dashed lines
    • I outline the 2010 aftershock region in blue dashed lines
    • I outline the 2015 aftershock region in white dashed lines
    • I outline the 2017 aftershock region in red dashed lines

Below is my interpretive poster for this earthquake. Click on the map to enlarge.

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely.
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault.

    I include one inset figure in the poster.

  • In the upper right corner I present Figure 2 from Beck et al. (1998 ) on the map, the space-time plot of historic and prehistoric earthquakes associated with the Chile subduction zone. I add a green line showing my interpretation for the strike length of the 2015 M 8.3 earthquake. Originally it appeared to match the 1943 and 1880 earthquakes, though it appears to extend further along strike. The 1922 and 1880 strike lengths are not well constrained, so this 2015 earthquake may indeed be slipping the same patch of this part of the subduction zone. Indeed, Juan Fernandez Ridge may be a structural boundary that may cause segmentation in this part of the subduction zone. If it does, it does not do so every time, as evidenced by the strike-length of the 1730 AD and 1647 AD earthquakes. I placed a green triangle at the approximate location of this 2017 swarm. This M 6.9 appears to be correlative in space with the 1985 earthquake (albeit a much smaller magnitude, closer to the 1971 in size).


Some background about the heterogeneous megathrust in this region

  • Here is the first of two figures from Moreno et al., 2010. Note that the M 6.9 is close in space to the 1985 earthquake. Also note the along strike heterogeneous seismogenic coupling. I include the figure caption below in blockquote.

  • Tectonic setting of the study area, data, observations and results. a, Shaded relief map of the Andean subduction zone in South- Central Chile. Earthquake segmentation along the margin is indicated by ellipses that enclose the approximate rupture areas of historic earthquakes (updated from refs 4–6). The inset shows the location of panel a (rectangle) relative to the South American continent. b, Compilation of GPS-observed surface velocities (1996–2008) with respect to stable South America before the 2010 Maule earthquake (for references see online-only Methods). Ellipses attached to the arrows represent 95% confidence limits. c, GPS 1 FEM modelled interface locking (fraction of plate convergence) distribution along the Andean subduction zone megathrust in the decade before the 2010 Maule earthquake. The epicentre (white star, USGS NEIC) and focal mechanism (beach ball, GCMT, http://www.globalcmt.org) of the 2010 Maule earthquake are shown in panels a and c.

  • Here is the second of the two figures from Moreno et al. (2010).

  • Relationship between pre, co- and postseismic deformation patterns. a, Coseismic slip distribution during the 2010 (blue contours; USGS slip model26) and 1960 (green contours; from ref. 30) earthquakes overlain onto pre-seismic locking pattern (red shading $0.75), as well as early (during the first 48 h post-shock) M$5 aftershock locations (the grey circle sizes scale with magnitude; GEOFON data29). b, Histograms of early (first 48 h; total number of events, 80) and late (first 3 months; total number of events, 168) aftershock density along a north–south profile (GEOFON data29, M$5). c, Residual slip deficits since 1835 as observed after the 2010 earthquake along a north–south profile (left column, based on the USGS slip model26). The middle and right columns show the effects on slip deficit of overlapping twentieth-century earthquakes (the black lines are polynomial fits to the data). Coloured data points and dates indicate earthquakes by year of occurrence.

References:

Earthquake Report: Chile!

Well, we had another earthquake in the region of a recent (yesterday and the day before) swarm offshore of Valparaiso, Chile (almost due west of Santiago, one of the largest cities in Chile). My previous report on the M 4-5 earthquakes can be found here. The earlier swarm was a series of shallower earthquakes (though some were of intermediate depth and some were deeper). The M 6.9 earthquake, in contrast, is deeper and likely on the megathrust. The slab contours are at 20 km and the hypocentral depth is 25 km (pretty good match considering the uncertainty with the location of the megathrust). Another difference is that the M 6.9 has a greater potential (likelihood, or chance) to damage people or their belongings.

Here are the USGS websites for these earthquakes

  • 2017.04.22 22:46 M 4.9
  • 2017.04.23 01:49 M 4.5
  • 2017.04.23 02:36 M 5.9 (mainshock)
  • 2017.04.23 02:43 M 4.8
  • 2017.04.23 02:52 M 4.8
  • 2017.04.23 03:00 M 4.8
  • 2017.04.23 03:02 M 4.9
  • 2017.04.23 19:40 M 5.6
  • 2017.04.24 21:38 M 6.9 (triggered mainshock)

I took a look at the seismicity from the past century. Here are Google Earth kml files from the USGS website for earthquakes from 1917-2017 with magnitudes M ≥ 5.0, M ≥ 6.0, and M ≥ 7.0.

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include the USGS epicenters for earthquakes from 1917-2017 with magnitudes M ≥ 6.0. I outline the regions of the subduction zone that have participated in earthquake slip during the 21st century (in white dashed polygons). I include USGS moment tensors from the largest earthquakes. I plot the focal mechanism for the 1960 earthquake from Moreno et al. (2011). Note the gap in seismicity in the region of the 1960 M 9.5 earthquake, except for the 2016 M 7.6 earthquake. Also, note how the 1960 and 2010 earthquake slip patches overlap.

Much of the subduction zone has ruptured, except for some spots between the 2001 and 2015 earthquakes. In 2015, I speculated that the region north of the 2015 earthquakes constituted a seismic gap. This region may get filled by a Great subduction zone earthquake or may continue to slip in moderate sized earthquakes (or be aseismic). There was an earthquake in 1877 that spanned 19-23 degrees (overlapping with the 2014 earthquake). This is shown on the Schurr et al. (2014) figure below).

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely.
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault.

    I include some inset figures in the poster.

  • In the lower left corner, I include a map and a cross section of the subduction zone just to the south of this Sept/Nov 2015 swarm (Melnick et al., 2006). I placed a green triangle at the approximate location of this 2017 swarm.
  • In the upper right corner I present Figure 2 from Beck et al. (1998 ) on the map, the space-time plot of historic and prehistoric earthquakes associated with the Chile subduction zone. I add a green line showing my interpretation for the strike length of the 2015 M 8.3 earthquake. Originally it appeared to match the 1943 and 1880 earthquakes, though it appears to extend further along strike. The 1922 and 1880 strike lengths are not well constrained, so this 2015 earthquake may indeed be slipping the same patch of this part of the subduction zone. Indeed, Juan Fernandez Ridge may be a structural boundary that may cause segmentation in this part of the subduction zone. If it does, it does not do so every time, as evidenced by the strike-length of the 1730 AD and 1647 AD earthquakes. I placed a green triangle at the approximate location of this 2017 swarm. This M 6.9 appears to be correlative in space with the 1985 earthquake (albeit a much smaller magnitude, closer to the 1971 in size).
  • In the lower right corner I include two figures from Moreno et al. (2010). The upper one shows the spatial extent of historic subduction zone earthquakes in this region, the GPS velocities, and the fraction of plate convergence attributed to fault seismogenic coupling. The lower panel shows the amount of slip that is attributed to the 1960 and 2010 earthquakes (on the left) and various measures of seismicity and slip deficit (on the right). I place a green star in the general location of the M 6.9 and a green horizontal bar that matches the latitude of this M 6.9 earthquake.
  • In the upper left corner, I include a local map showing the MMI contours for the M 6.9 earthquake. I include the USGS moment tensors from most of the earthquakes in this swarm, including the M 6.9 earthquake.


  • As mentioned above, this earthquake has the potential to cause more harm than the earlier earthquakes due to its larger magnitude. Below is the USGS report that includes estimates of damage to people (possible fatalities) and their belongings from the Rapid Assessment of an Earthquake’s Impact (PAGER) report. More on the PAGER program can be found here. An explanation of a PAGER report can be found here. PAGER reports are modeled estimates of damage. On the top is a histogram showing estimated casualties and on the right is an estimate of possible economic losses. This PAGER report suggests that there will be quite a bit of damage from this earthquake (and casualties). This earthquake has a high probability of damage to people and their belongings.

  • UPDATE: Below are some observations of the tsunami. This comes from the Pacific Tsunami Warning Center.

  • Here is the figure from Lin et al. (2013) that shows the tectonic context of the 2010 Maule earthquake. I include the figure captions as blockquote.

  • (a) Regional tectonic map showing slab isodepth contours (blue lines) [Cahill and Isacks, 1992], M>=4 earthquakes from the National Earthquake Information Center catalog between 1976 and 2011 (yellow circles for depths less than 50 km, and blue circles for depths greater than 50 km), active volcanoes (red triangles), and the approximate extent of large megathrust earthquakes during the past hundred years (red ellipses) modified from Campos et al. [2002]. The large white vector represents the direction of Nazca Plate with respect to stable South America [Kendrick et al., 2003]. (b) Simplified seismo-tectonic map of the study area. Major Quaternary faults are modified after Melnick et al. [2009] (black lines). The Neogene Deformation Front is modified from Folguera et al. [2004]. The west-vergent thrust fault that bounds the west of the Andes between 32 and 38S is modified from Melnick et al. [2009]. (c) Schematic cross-section along line A–A0 (Figure 1b), modified from Folguera and Ramos [2009]. The upper bound of the coseismic slip coincides with the boundary between the frontal accretionary prism and the paleo-accretionary prism [Contreras-Reyes et al., 2010], whereas the contact between the coseismic and postseismic patch is from this study. The thick solid red line and dashed red line on top of the slab represent the approximate coseismic and postseismic plus interseismic slip section of the subduction interface. The thin red and grey lines within the overriding plate are active and inactive structures in the retroarc, modified from Folguera and Ramos [2009]. The red dashed line underneath the Andean Block represents the regional décollement. Background seismicity is from the TIPTEQ catalog, recorded between November 2004 and October 2005 [Rietbrock et al., 2005; Haberland et al., 2009].

  • Here is a cross section of the subduction zone just to the south of this Sept/Nov 2015 swarm (Melnick et al., 2006). Below I include the text from the Melnick et al. (2006) figure caption as block text.

  • (A) Seismotectonic segments, rupture zones of historical subduction earthquakes, and main tectonic features of the south-central Andean convergent margin. Earthquakes were compiled from Lomnitz (1970, 2004), Kelleher (1972), Comte et al. (1986), Cifuentes (1989), Beck et al. (1998), and Campos et al. (2002). Nazca plate and trench are from Bangs and Cande (1997) and Tebbens and Cande (1997). Maximum extension of glaciers is from Rabassa and Clapperton (1990). F.Z.—fracture zone. (B) Regional morphotectonic units, Quaternary faults, and location of the study area. Trench and slope have been interpreted from multibeam bathymetry and seismic-reflection profiles (Reichert et al., 2002). (C) Profile of the offshore Chile margin at ~37°S, indicated by thick stippled line on the map and based on seismic-reflection profiles SO161-24 and ENAP-017. Integrated Seismological experiment in the Southern Andes (ISSA) local network seismicity (Bohm et al., 2002) is shown by dots; focal mechanism is from Bruhn (2003). Updip limit of seismogenic coupling zone from heat-fl ow measurements (Grevemeyer et al., 2003). Basal accretion of trench sediments from sandbox models (Lohrmann, 2002; Glodny et al., 2005). Convergence parameters from Somoza (1998 ).

  • Here is the first of two figures from Moreno et al., 2010. Note that the M 6.9 is close in space to the 1985 earthquake. I include the figure caption below in blockquote.

  • Tectonic setting of the study area, data, observations and results. a, Shaded relief map of the Andean subduction zone in South- Central Chile. Earthquake segmentation along the margin is indicated by ellipses that enclose the approximate rupture areas of historic earthquakes (updated from refs 4–6). The inset shows the location of panel a (rectangle) relative to the South American continent. b, Compilation of GPS-observed surface velocities (1996–2008) with respect to stable South America before the 2010 Maule earthquake (for references see online-only Methods). Ellipses attached to the arrows represent 95% confidence limits. c, GPS 1 FEM modelled interface locking (fraction of plate convergence) distribution along the Andean subduction zone megathrust in the decade before the 2010 Maule earthquake. The epicentre (white star, USGS NEIC) and focal mechanism (beach ball, GCMT, http://www.globalcmt.org) of the 2010 Maule earthquake are shown in panels a and c.

  • Here is the second of the two figures from Moreno et al. (2010).

  • Relationship between pre, co- and postseismic deformation patterns. a, Coseismic slip distribution during the 2010 (blue contours; USGS slip model26) and 1960 (green contours; from ref. 30) earthquakes overlain onto pre-seismic locking pattern (red shading $0.75), as well as early (during the first 48 h post-shock) M$5 aftershock locations (the grey circle sizes scale with magnitude; GEOFON data29). b, Histograms of early (first 48 h; total number of events, 80) and late (first 3 months; total number of events, 168) aftershock density along a north–south profile (GEOFON data29, M$5). c, Residual slip deficits since 1835 as observed after the 2010 earthquake along a north–south profile (left column, based on the USGS slip model26). The middle and right columns show the effects on slip deficit of overlapping twentieth-century earthquakes (the black lines are polynomial fits to the data). Coloured data points and dates indicate earthquakes by year of occurrence.

  • Here is the Beck et al. (1998) space time diagram.

Here is an animation of seismicity from the 21st century

Useful Resources

References:

Earthquake Report: Chile!

There have been a number of earthquakes along the subduction zone offshore of Chile. These have happened near the boundary of two Great Earthquakes from 2010 and 2015. This region may be a segment boundary along the subduction zone, albeit possibly a non persistent one. The Juan Ferndandez ridge may control this segmentation.

The earthquakes from today and yesterday form a range of about 1 1/2 magnitudes (M 4.2- M 5.9). This may be considered a swarm (when there are a series of earthquakes along a fault with similar magnitudes), though there is an M 5.9 that could be considered the mainshock. But, I would not get hung up on terminology as that is not very important. However, there is a great page with a discussion about swarms, including some good examples.

Here are the USGS websites for these earthquakes

I took a look at the seismicity from the past century. Here are Google Earth kml files from the USGS website for earthquakes from 1917-2017 with magnitudes M ≥ 5.0, M ≥ 6.0, and M ≥ 7.0.

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include the USGS epicenters for earthquakes from 1917-2017 with magnitudes M ≥ 5.0. I outline the regions of the subduction zone that have participated in earthquake slip during the 21st century (in white dashed polygons). I include USGS moment tensors from the largest earthquakes. I plot the focal mechanism for the 1960 earthquake from Moreno et al. (2011). Note the gap in seismicity in the region of the 1960 M 9.5 earthquake, except for the 2016 M 7.6 earthquake. Also, note how the 1960 and 2010 earthquake slip patches overlap.

Much of the subduction zone has ruptured, except for some spots between the 2001 and 2015 earthquakes. In 2015, I speculated that the region north of the 2015 earthquakes constituted a seismic gap. This region may get filled by a Great subduction zone earthquake or may continue to slip in moderate sized earthquakes (or be aseismic). There was an earthquake in 1877 that spanned 19-23 degrees (overlapping with the 2014 earthquake). This is shown on the Schurr et al. (2014) figure below).

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely.
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault.

    I include some inset figures in the poster.

  • In the lower left corner, I include a figure from Lin et al. (2013) that shows the tectonic context of the 2010 Maule earthquake. On the map are plotted extents of historic earthquakes along this convergent plate margin. On the right is a large scale map showing the active magmatic arc volcanoes associated with this subduction zone. Finally, there is a cross section showing where the coseismic slip and postseismic slip occurred as part of the 2010 earthquake sequence. I placed a green triangle at the approximate location of this 2017 swarm.
  • In the lower right corner, I include a time-space diagram from Moernaut et al. (2010). There is also a map showing the fracture zones. I placed a green triangle at the approximate location of this 2017 swarm.
  • Above the Moernaut et al. (2010) figure, I present Figure 2 from Beck et al. (1998 ) on the map, the space-time plot of historic and prehistoric earthquakes associated with the Chile subduction zone. This space-time plot overlaps slightly with the Moernaut figure. I add a green line showing my interpretation for the strike length of the 2015 M 8.3 earthquake. Originally it appeared to match the 1943 and 1880 earthquakes, though it appears to extend further along strike. The 1922 and 1880 strike lengths are not well constrained, so this 2015 earthquake may indeed be slipping the same patch of this part of the subduction zone. Indeed, Juan Fernandez Ridge may be a structural boundary that may cause segmentation in this part of the subduction zone. If it does, it does not do so every time, as evidenced by the strike-length of the 1730 AD and 1647 AD earthquakes. I placed a green triangle at the approximate location of this 2017 swarm.
  • In the upper right corner is a space-time figure showing earthquakes for the past few centuries. This diagram does not overlap with the Beck figure. This figure shows the outline of some subduction zone earthquakes and shows how the 2014 earthquake is composed of two earthquakes (an M 8.1 and an M 7.6) that ruptured different but adjacent patches of the subduction zone.
  • In the upper left corner, I include a local map showing the MMI contours for the M 5.9 earthquake. I include the USGS moment tensors from most of the earthquakes in this swarm.


  • Here is the figure from Lin et al. (2013) that shows the tectonic context of the 2010 Maule earthquake. I include the figure captions as blockquote.

  • (a) Regional tectonic map showing slab isodepth contours (blue lines) [Cahill and Isacks, 1992], M>=4 earthquakes from the National Earthquake Information Center catalog between 1976 and 2011 (yellow circles for depths less than 50 km, and blue circles for depths greater than 50 km), active volcanoes (red triangles), and the approximate extent of large megathrust earthquakes during the past hundred years (red ellipses) modified from Campos et al. [2002]. The large white vector represents the direction of Nazca Plate with respect to stable South America [Kendrick et al., 2003]. (b) Simplified seismo-tectonic map of the study area. Major Quaternary faults are modified after Melnick et al. [2009] (black lines). The Neogene Deformation Front is modified from Folguera et al. [2004]. The west-vergent thrust fault that bounds the west of the Andes between 32 and 38S is modified from Melnick et al. [2009]. (c) Schematic cross-section along line A–A0 (Figure 1b), modified from Folguera and Ramos [2009]. The upper bound of the coseismic slip coincides with the boundary between the frontal accretionary prism and the paleo-accretionary prism [Contreras-Reyes et al., 2010], whereas the contact between the coseismic and postseismic patch is from this study. The thick solid red line and dashed red line on top of the slab represent the approximate coseismic and postseismic plus interseismic slip section of the subduction interface. The thin red and grey lines within the overriding plate are active and inactive structures in the retroarc, modified from Folguera and Ramos [2009]. The red dashed line underneath the Andean Block represents the regional décollement. Background seismicity is from the TIPTEQ catalog, recorded between November 2004 and October 2005 [Rietbrock et al., 2005; Haberland et al., 2009].

  • Here is a cross section of the subduction zone just to the south of this Sept/Nov 2015 swarm (Melnick et al., 2006). Below I include the text from the Melnick et al. (2006) figure caption as block text.

  • (A) Seismotectonic segments, rupture zones of historical subduction earthquakes, and main tectonic features of the south-central Andean convergent margin. Earthquakes were compiled from Lomnitz (1970, 2004), Kelleher (1972), Comte et al. (1986), Cifuentes (1989), Beck et al. (1998), and Campos et al. (2002). Nazca plate and trench are from Bangs and Cande (1997) and Tebbens and Cande (1997). Maximum extension of glaciers is from Rabassa and Clapperton (1990). F.Z.—fracture zone. (B) Regional morphotectonic units, Quaternary faults, and location of the study area. Trench and slope have been interpreted from multibeam bathymetry and seismic-reflection profiles (Reichert et al., 2002). (C) Profile of the offshore Chile margin at ~37°S, indicated by thick stippled line on the map and based on seismic-reflection profiles SO161-24 and ENAP-017. Integrated Seismological experiment in the Southern Andes (ISSA) local network seismicity (Bohm et al., 2002) is shown by dots; focal mechanism is from Bruhn (2003). Updip limit of seismogenic coupling zone from heat-fl ow measurements (Grevemeyer et al., 2003). Basal accretion of trench sediments from sandbox models (Lohrmann, 2002; Glodny et al., 2005). Convergence parameters from Somoza (1998 ).

  • In March 2015, there was some seismicity in this September/November 2015 earthquake slip region. I put together an earthquake report about those earthquake of magnitudes M = 5.0-5.3. I speculate that the 1922 earthquake region is a seismic gap. Note that this September/November 2015 earthquake region is along the southern portion of the seismic gap that I labeled on the map below.
  • Here is a map that shows the recent swarm of ~M = 5 earthquakes. There are moment tensors for the earthquakes listed below, some recent historic subduction zone earthquakes. I placed the general along-strike distance for older historic earthquakes in green (and labeled their years). The largest earthquake ever recorded, the Mw = 9.5 Chile earthquake, had a slip patch that extends from the south of the map to just south of the 2010 earthquake swarm. The 2010 and 2014 earthquake swarm epicenters are plotted as colored circles, while most other historic earthquake epicenters are plotted as gray circles. Note how this March 2015 swarm is at the northern end of the 1922/11/11 M 8.3 earthquake. At the bottom of this page, I put a USGS graphic about what these moment tensor plots (beach balls) tell us about the earthquakes.

  • Here is the first space-time figure from Schurr et al., 2014. I include their caption as blockquote below.

  • Map of Northern Chile and Southern Peru showing historical earthquakes and instrumentally recorded megathrust ruptures. IPOC instruments used in the present study (BB, broadband; SM, strong motion) are shown as blue symbols. Left: historical1,2 and instrumental earthquake record. Centre: rupture length was calculated using the regression suggested in ref. 28, with grey lines for earthquakes M .7 and red lines for Mw .8. The slip distribution of the 2014 Iquique event and its largest aftershock derived in this study are colour coded, with contour intervals of 0.5 m. The green and black vectors are the observed and modelled horizontal surface displacements of the mainshock. The slip areas of the most recent other large ruptures4,5,7 are also plotted. Right: moment deficit per kilometre along strike left along the plate boundary after the Iquique event for moment accumulated since 1877, assuming current locking (Fig. 3a). The total accumulated moment since 1877 from 17u S to 25u S (red solid line) is 8.97; the remaining moment after subtracting all earthquake events with Mw .7 (grey dotted line) is 8.91 for the entire northern Chile–southern Peru seismic gap

  • Here is the Beck et al. (1998) space time diagram.

  • Finally, here is the southernmost space-time diagram from Moernaut et al. (2010). These data are largely derived from Melnick et al. (2009).

  • Setting and historical earthquakes in South-Central Chile. Data derived from Barrientos (2007); Campos et al. (2002); Melnick et al.(2009).

Here is an animation of seismicity from the 21st century

References:

Earthquake Report: Phillipines

I put these together earlier this week for me classes and finally have a moment to write about these earthquakes. The Philippines region has been quite active lately, as it frequently is. I show below a series of earthquakes from the past ~30 days. These earthquakes occurred in 4 different regions and 3 different tectonic settings. These are probably unrelated to each other, but it is difficult to really know without further analyses. After I made these posters, there was an earthquake with a magnitude of M 5.8 on the island of Mindanao (I include the USGS link below), possibly associated with the Davao River fault (the closest fault mapped in this region).

Here are the USGS websites for these earthquakes.

I took a look at the seismicity from the past century. Here are Google Earth kml files from the USGS website for earthquakes from 1917-2017 with magnitudes M ≥ 7.0 and M ≥ 7.5.

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include the USGS epicenters for earthquakes from 1917-2017 with magnitudes M ≥ 7.5.

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely. The
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault.

    I include some inset figures in the poster.

  • In the upper left corner is a figure from Hall, 2011. This shows the plate tectonic configuration in the equatorial Pacific. Note how the upper panel shows a west dipping slab on the east side of the Philippines. Note the contrast in the center panel (Halmahera), where the eastern fault is dipping to the east (westward vergent) and the western fault is dipping to the west (eastward vergent). This region near Halmahera forms the Molucca Strait, one of the most tectonically active areas in this region.
  • In the upper right corner is a map showing the regional faults from Noda (2013). The Philippine and
  • In the lower left corner I include a map showing the seismicity and tectonic plate boundary faults for this region (Smoczyk et al., 2013). Earthquakes are plotted with color representing depth and diameter representing magnitude (see legend).


  • I choose 4 earthquakes for which I plot the MMI intensity.
  • The easternmost two are earthquakes related to the subduction of the Philippine Sea plate.
  • The northwesternmost earthquake series (largest M = 5.9) is interesting and mysterious. There is a left-lateral forearc sliver fault (Philippine fault) that parallels the Philippine trench. Internal deformation in the upper plate is accommodated by a complicated series of other strike-slip faults. In the region near Manila, there is a north-south striking right-lateral Mirikina Valley fault system. The MV fault trends south from Manila into the region of the M 5.9 earthquake. This fault seems to terminate in a northeast-southwest trending zone of “extension and young volcanism” (Nelson et al., 2000). Further to the south is an east-west striking left-lateral Lubang fault. This may be an extension of the Sibuyan Sea fault (Noda, 2013), a splay of the Philippine fault. If we look at the moment tensor for the M 5.5 and 5.9 earthquakes, the nodal planes suggest either NE striking left-lateral or NW striking right-lateral motion. This does not fit the orientation nor sense of motion for any of the mapped faults in the region. It is possible that these earthquakes are related to the extensional rifting instead. The moment tensor for the M 5.1 earthquake here shows an extensional beach ball. The orientation is not completely correct based upon the Nelson et al. (2000) figure, but there are no faults shown on their map.
  • The southwestern earthquakes from March (both M 5.6) are also interesting. The epicenters show locations on north Sulawesi, part of Indonesia. There is a NW striking left-lateral strike-slip faults mapped to the west and east of these earthquakes. However, the moment tensor shows that this would be a right-lateral fault in that orientation. The NW striking reverse fault is also equally challenging to interpret. Needless to say, this is a tectonically complicated region.

Below is my interpretive poster for the M 5.9 earthquake.


  • This is the low-angle oblique view of the region (Hall, 2011).

  • 3D cartoon of plate boundaries in the Molucca Sea region modified from Hall et al. (1995). Although seismicity identifies a number of plates there are no continuous boundaries, and the Cotobato, North Sulawesi and Philippine Trenches are all intraplate features. The apparent distinction between different crust types, such as Australian continental crust and oceanic crust of the Philippine and Molucca Sea, is partly a boundary inactive since the Early Miocene (east Sulawesi) and partly a younger but now probably inactive boundary of the Sorong Fault. The upper crust of this entire region is deforming in a much more continuous way than suggested by this cartoon.

  • Here is the fault map for the region on the Island of Luzon (where Manila is located) from Nelson et al. (2000). The panel on the right (B) shows the Marikina Valley fault system (MV) and the Lubang fault. The MV is a right-lateral strike slip fault and the Lubang fault is a left-lateral strike-slip fault.

  • Tectonic setting of the Marikina Valley fault system (MV) in central Luzon, the Philippines. Diagram A shows subduction zone trenches by barbed lines, other faults with high rates of Quaternary activity by heavy black lines. White dots show locations of recent earthquakes on the Philippine fault in Luzon (M 7.8; 1990) and the Aglubang River fault in Mindoro (M 7.1; 1994). Diagram B shows how the Marikina Valley pull-apart basin (MV) may have been formed through extension caused by clockwise rotation (dashed circle) and shearing of central Luzon, which is caught between two active left-lateral strike-slip faults—the Philippine fault (Nakata et al., 1977; Barrier et al., 1991; Ringenbach et al., 1993; Aurelio et al., 1993) and the Lubang fault. A zone of extension and young volcanism south of the fault system has also influenced the structural development of the valley (Fo¨rster et al., 1990; Defant et al., 1988).

  • Here is the map from Noda (2013) that shows various strike slip faults associated with subduction zones.

  • Modern examples of trench-linked strike-slip faults. (A) The Median Tectonic Line (MTL) active fault system in southwestern Japan, related to oblique subduction of the Philippine Sea Plate (PS) along the Nankai Trough (NT). (B) The Great Sumatra Fault system (GSF) along the Java–Sumatra Trench (JST). (C) Strike-slip faults in Alaska. Fault names: DF, Denali; BRF, Boarder Ranges; CSEF, Chugach St. Elias; FF, Fairweather; TF, Transition. (D) The Philippine Fault system (PF). Abbreviations: SSF, Sibuyan Sea Fault; MT, Manila Trench; PT, Philippine Trench; ELT, East Luzon Trough. Plate names: AM, Amur; OK, Okhotsk; PS, Philippine Sea; AU, Australian; SU, Sundaland; NA, North American; PA, Pacific; YMC, Yukutat microcontinent. Black and purple lines are subduction zones and trench-linked strike-slip faults, respectively. All maps were drawn using SRTM and GEBCO with plate boundary data [30]. Blue arrows indicate the direction and velocity of relative plate motion (mm yr-1) based on [31].

The USGS Maps and Cross-Sections

  • Here is the map from Smocyk et al., 2013, followed by the legend. The entire poster is here (92 MB pdf).


  • Below are two cross sections that show the subduction zone seismicity, followed by the legend. The location of these cross sections are labeled on the map above.



  • This map shows the seismic hazard for this region. The color represents the likelihood of any region experiencing ground shaking of a particular magnitude. The scale is “Peak Ground Acceleration.” Units are m/s^2. Purple represents gravitational acceleration of 1 g, gravity at Earth’s surface. Note how most of the earthquakes were in the region of higher likely ground shaking, except for the Sulawesi earthquakes.

  • In January of this year, there was an M 7.3 earthquake in the Celebes Sea south of the Philippines. Below is my interpretive map for that earthquake. I also present the same poster with 1917-2017 seismicity for earthquakes M ≥ 6.5. Here is my earthquake report for this M 7.3 earthquake. I include more background information for the Molucca Strait region on this page.


References:

  • Bock et al., 2003. Crustal motion in Indonesia from Global Positioning System measurements in JGR, v./ 108, no. B8, 2367, doi:10.1029/2001JB000324
  • Hall, R., 2011. Australia–SE Asia collision: plate tectonics and crustal flow in Hall, R., Cottam, M. A. &Wilson, M. E. J. (eds) The SE Asian Gateway: History and Tectonics of the Australia–Asia Collision. Geological Society, London, Special Publications, 355, 75–109.
  • 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
  • McCaffrey, R., Silver, E.A., and Raitt, R.W., 1980. Crustal Structure of the Molucca Sea Collision Zone, Indonesia in The Tectonic and Geologic Evolution of Southeast Asian Seas and Islands-Geophysical Monograph 23, p. 161-177.
  • Nelson, A.R., Personius, S.F., Rimando, R.E>, Punongbayan, R.S., Tungol, N, Mirabueno, H., and Rasdas, A., 2000. Multiple Large Earthquakes in the Past 1500 Years on a Fault in Metropolitan Manila, the Philippines in BSSA vol. 90, p. 73-85.
  • Noda, A., 2013. Strike-Slip Basin – Its Configuration and Sedimentary Facies in Mechanism of Sedimentary Basin Formation – Multidisciplinary Approach on Active Plate Margins http://www.intechopen.com/books/mechanism-of-sedimentarybasin-formation-multidisciplinary-approach-on-active-plate-margins http://dx.doi.org/10.5772/56593
  • Smoczyk, G.M., Hayes, G.P., Hamburger, M.W., Benz, H.M., Villaseñor, Antonio, and Furlong, K.P., 2013. Seismicity of the Earth 1900–2012 Philippine Sea plate and vicinity: U.S. Geological Survey Open-File Report 2010–1083-M, 1 sheet, scale 1:10,000,000.
  • Waltham et al., 2008. Basin formation by volcanic arc loading in GSA Special Papers 2008, v. 436, p. 11-26.
  • Zahirovic et al., 2014. The Cretaceous and Cenozoic tectonic evolution of Southeast Asia in Solid Earth, v. 5, p. 227-273, doi:10.5194/se-5-227-2014.

Earthquake Report: Botswana!

This is a very interesting M 6.5 earthquake, which was preceded by a probably unrelated M 5.2 earthquake.

    Here are the USGS web pages for these earthquakes


  • 2017.04.03 M 5.2
  • 2017.04.03 M 6.5

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past century, with color representing depth and diameter representing magnitude (see legend). I include the USGS epicenters for earthquakes from 1917-2017 with magnitudes M ≥ 2.5. The M 5.2 earthquake happened in a region that is seismically active and this preceded the M 6.5 earthquake. They are at a large distance and are unlikely related to each other.

I also include the generalized location of the East Africa Rift (EAR) in this region as yellow bands with white dashed lines. These are the Eastern Branch and Southwestern Branch of the EAR.

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely. This earthquake is a normal fault event that strikes to the northwest and dips to the northeast or southwest. Due to the paucity of seismic data and mapped faults here, it is difficult to tell which is the principal fault plane. Someone online suggested that this may not be an earthquake, but an alien invasion force. Just joking. However, it is possible that this might have been an underground explosion. But, at 25 km, this seems highly unlikely.
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.

    I include some inset figures in the poster.

  • In the lower right corner I include a geologic map for southern Africa and Botswana (Leseane et al., 2015). Note that the Southwestern Branch of the EAR extends into northwestern Botswana (the white fault lines in the red colored geologic unit). I place a red star in the general location of the M 6.5 earthquake. This earthquake happened in the Kaapvaal Craton, a region of low historical seismicity.
  • In the upper right corner I include a map that shows the seismic hazard for the EAR region of Africa (Hayes et al., 2014). Today’s M 6.5 earthquake happened in the lower left corner of the map. I place a red star in the general location of the M 6.5 earthquake. Note that this earthquake happened in a region of low seismic hazard.
  • In the upper left corner I include two maps and their associated plots. The map on the left is the shaking intensity map as modeled by the USGS, which uses the MMI scale. The plot below shows the results from the Ground Motion Prediction Equation (GMPE) used to generate the map. This plot shows how ground shaking (MMI) attenuates (diminishes as the seismic waves are absorbed by the earth) with distance from the earthquake. The map on the right is the “Did You Feel It?” map generated by the online responses from people who observed the ground shaking. This also uses the MMI scale. The plot below that shows how the reported observations match the GMPE relations for this earthquake. The GMPE relations are plotted as orange and green lines, which represent GMPE models developed for different geologic settings (e.g. hard rock like granite vs. soft rock like accreted terrane).
  • In the lower left corner I include the Rapid Assessment of an Earthquake’s Impact (PAGER) report. More on the PAGER program can be found here. An explanation of a PAGER report can be found here. PAGER reports are modeled estimates of damage. On the top left is a histogram showing estimated casualties and on the top right is an estimate of possible economic losses. There is a list of cities in the lower right corner which shows their populations and the MMI that they were likely exposed to.


  • This is the geologic map from the poster (Leseane et al., 2015). I include their caption below in blockquote.

  • Precambrian tectonic map of (a) southern Africa and (b) Botswana outlining the spatial extent of Archean cratons and Proterozoic orogenic belts. White lines represent the fault system of the Okavango Rift Zone. Modified after Singletary et al. [2003] and Begg et al. [2009].

  • Here is the USGS “Seismicity of the Earth” poster for this region (Hayes et al., 2014).

  • This is the latest geologic maps of Africa (Thieblemont, D., 2016). Click on the map for a 67 MB pdf version.

References:

  • Hayes, G.P., Jones, E.S., Stadler, T.J., Barnhart, W.D., McNamara, D.E., Benz, H.M., Furlong, K.P., and Villaseñor, Antonio, 2014, Seismicity of the Earth 1900–2013 East African Rift: U.S. Geological Survey Open-File Report 2010–1083-P, 1 sheet, scale 1:8,500,000 http://dx.doi.org/10.3133/of20101083P
  • Leseane, K., Atekwana, E.A., Mickus, K.L., Abdelsalam, M.G., Shemanq, E.M., and Atekwana, E.A., 2015. Thermal perturbations beneath the incipient Okavango Rift Zone, northwest Botswana in JGR: Solid Earth, v. 120, doi:10.1002/2014JB011029.
  • Thieblemont, D. (ed.), 2016. Geological Map of Africa et 1:10M scale, CGMW-BRGM 2016

Earthquake Report: Gulf California!

There was an earthquake yesterday in the Gulf of California nearby a series of earthquakes that happened in 2015 and earlier in 2013. The 2017 and 2013 earthquakes are happening along a fault that forms the Carmen Basin and the 2015 earthquakes are rupturing a fault that appears to be in the middle of the Farallon Basin. Here is my Earthquake Report for the 2013 earthquake (an early report, so it is rather basic). Here is my Earthquake Report for the 2015 earthquake sequence. This is an update to the initial 2015 report.

    Here are the USGS web pages for these earthquakes


    2013

  • 2013.10.19 M 6.4

  • 2015

  • 2015.09.13 M 6.6
  • 2015.09.13 M 5.3
  • 2015.09.13 M 4.9
  • 2015.09.13 M 5.2

  • 2017

  • 2017.03.29 M 5.7

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include the USGS epicenters for earthquakes from 1917-2017 with magnitudes M ≥ 6.5.

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely. I suspect that the fault that ruptured is eastward vergent (dipping to the west), so the west dipping nodal plane is probably the primary fault plane. However, this region of Kamchatka has numerous upper plate thrust and reverse faults (so the primary fault plane could be the other one, dipping to the east).
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.

    I include some inset figures in the poster.

  • In the upper left corner I include a map that shows the tectonic setting of this region, with the geological units colored relative to their age and type (marine or continental).This is from a paper that discusses the interaction between spreading ridges and subduction trenches (Fletcher et al., 2007).
  • In the upper right corner I include a larger saled version of this map with the same seismicity plotted, but without the MMI contours. One can see the shape of the seafloor and how this is formed due to plate tectonics. Note the spreading ridge in the lower right corner and the parallel ridges formed as the plates extend from the ridge. The magnitude scale is slightly different than the main map.
  • In the lower left corner I include a map from the 2015 earthquake series. I include this because I labeled the Basins formed by the enechelon steps in this plate boundary. In my 2015 report, I provide more maps that include the names of some of these fracture zones (the transform plate boundary strike-slip faults).


  • Here is a great diagram showing the major faults in the region. I include their figure caption below.

  • (A) Simplified map of the Gulf of California region and Baja California peninsula showing the present plate boundary and some major tectonic features related to the plate-tectonic history since 12 Ma. The Gulf extensional province in gray is bounded by the Main Gulf Escarpment (bold dashed lines), which runs through the Loreto area and is shown in Figure 3. The Salton trough in southern California is merely the northern part of the Gulf extensional province. (B) Map of part of the southern Gulf of California and Baja California peninsula showing bathymetry (in meters), the transform–spreading-ridge plate boundary, and the location of subsequent figures with maps. The bathymetry is after a map in Ness and Lyle (1991) and the transform–spreading-ridge plate boundary is from Lonsdale (1989). The lines with double arrows are the three proposed rift segments modified here after Axen (1995); MS—Mulege´ segment, LS—Loreto segment, TS—Timbabichi
    segment.

  • This map shows the magnetic anomalies and the geologic map for the land and the youngest oceanic crust.

  • (A) Tectonic map of the southern Baja California microplate (BCM) and Gulf of California extensional province (GEP). The Magdalena fan is deposited on oceanic crust of the Farallon-derived Magdalena microplate located west of Baja California. Deep Sea Drilling Project Site 471 is shown as black dot on the Magdalena fan. Abbreviations: BCT—Baja California trench, BM—Bahia Magdalena, LC—Los Cabos block, T—Trinidad block, LP—La Paz, PV—Puerto Vallarta, SMSLF—Santa Margarita–San Lazaro fault, TAF—Tosco-Abreojos fault, TS—Todos Santos, V—Vizcaino peninsula. Geology is simplifi ed from Muehlberger (1996). Interpretation of marine magnetic anomalies, with numbers denoting the chron of positively magnetized stripes, is from Severinghaus and Atwater (1989) and Lonsdale (1991).

  • This map shows a more broad view of the magnetic anomalies through time.

  • Map-view time slices showing the widely accepted model for the two-phase kinematic evolution of plate margin shearing around the Baja California microplate. (A) Configuration of active ridge segments (pink) west of Baja California just before they became largely abandoned ca. 12.3 Ma. (B) It is thought that plate motion from 12.3 to 6 Ma was kinematically partitioned into dextral strike slip (325 km) on faults west of Baja California and orthogonal rifting in the Gulf of California (90 km). This is known as the protogulf phase of rifting. (C) From 6 to 0 Ma faults west of Baja California are thought to have died and all plate motion was localized in the Gulf of California, which accommodated ~345 km of integrated transtensional shearing. Despite its wide acceptance, our data preclude this kinematic model. In all frames, the modern coastline is blue. Continental crust that accommodated post–12.3 Ma shearing is dark brown. Unfaulted microplates of continental crust are light tan. Farallon-derived microplates are light green. Middle Miocene trench-filling deposits like the Magdalena fan are colored dark green. Deep Sea Drilling Project Site 471 is the black dot on the southern Magdalena microplate. Yellow line (296 km) in the northern Gulf of California connects correlated terranes of Oskin and Stock (2003). Maps have Universal Transverse Mercator zone 12 projection with mainland Mexico fixed in present position.

  • Here is the Earthquake Report Poster from the 2015 sequence.

This is a nice simple figure, from the University of Sydney here, showing the terminology of strike slip faulting. It may help with the following figures.

Here is a fault block diagram showing how strike-slip step overs can create localized compression (positive flower) or extension (negative flower). More on strike-slip tectonics (and the source of this image) here.


Here is another great figure showing how sedimentary basins can be developed as a result of step overs in strike slip fault systems (source: Becky Dorsey, University of Oregon, Dept. of Geological Sciences).


I also put together an animation of seismicity from 1065 – 2015. First, here is a map that shows the spatial extent of this animation.


Here is the animation link (2 MB mp4 file) if you cannot view the embedded video below. Note how the animation begins in 1965, but has the recent seismicity plotted for reference.

    There have been two large magnitude earthquakes in this region over the past 50 years.

  • 2007.09.01 M 6.1
  • 2010.10.21 M 6.7

This is an animation from Tanya Atwater. Click on this link to take you to yt (if the embedded video below does not work).

Here is an animation from IRIS. This link takes you to yt (if you cannot view the embedded version below). Here is a link to download the 21 MB mp4 vile file.

This is a link to a tectonic summary map from the USGS. Click on the map below to download the 20 MB pdf file.


References:

  • Fletcher, J.M., Grove, M., Kimbrough, D., Lovera, O., and Gehrels, G.E., 2007. Ridge-trench interactions and the Neogene tectonic evolution of the Magdalena shelf and southern Gulf of California: Insights from detrital zircon U-Pb ages from the Magdalena fan and adjacent areas in GSA Bulletin, v. 119, no. 11/12, p. 1313-1336.
  • Umhoefer, P.J., Mayer, L., and Dorse, R.J., 2002. Evolution of the margin of the Gulf of California near Loreto, Baja California Peninsula, Mexico in GSA Bulletin, v. 114, no. 7, p. 849-868.

Earthquake Report: Kamchatka!

This earthquake happened last night as I was preparing course materials for this morning. Initially it was a magnitude 6.9, but later modified to be M 6.9.

This earthquake happened in an interesting region of the world where there is a junction between two plate boundaries, the Kamchatka subduction zone with the Aleutian subduction zone / Bering-Kresla Shear Zone. The Kamchatka Trench (KT) is formed by the subduction (a convergent plate boundary) beneath the Okhtosk plate (part of North America). The Aleutian Trench (AT) and Bering-Kresla Shear Zone (BKSZ) are formed by the oblique subduction of the Pacific plate beneath the Pacific plate. There is a deflection in the Kamchatka subduction zone north of the BKSZ, where the subduction trench is offset to the west. Some papers suggest the subduction zone to the north is a fossil (inactive) plate boundary fault system. There are also several strike-slip faults subparallel to the BKSZ to the north of the BKSZ. These are shown in two of the inset maps below.

    Here is the USGS website for this earthquake.

  • 2017.03.29 M 6.6

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include the USGS epicenters for earthquakes from 1917-2017 with magnitudes M ≥ 7.0 (the search window is limited to the region west of the Amlia fracture zone).

  • I also include moment tensors for earthquakes associated with the Kamchatka-Kuril subduction zone. There are some interesting earthquakes plotted here:
    • The pair of earthquakes 2008.07.05 M 7.7 and 2013.05.24 M 8.3 are very deep earthquakes (the M 8.3 is one of the largest and deepest earthquake ever recorded by modern seismometers) and may be due to bending of the downgoing slab. Here is my report for the M 8.3 (it was an early report of mine).
    • The pair of earthquakes 2006.11.15 M 833 and 2007.01.13 M 8.1 are directly related to each other. Lay et al. (2009) discussed this earthquake sequence and how the 2006 subduction zone earthquake led to the 2007 outer-rise earthquake in the Pacific plate. Dr. Erica Emry studied this earthquake pair for her Ph.D. research. I also wrote a little about these earthquakes in my earthquake report here.
    • The largest earthquake plotted for this region is the 1952 M 9.0 earthquake (the large epicenter between teh 1993 and 1997 earthquakes. This earthquake is the 5th largest earthquake recorded by modern seismometers. However, there is no USGS moment tensor (I couldn’t find a focal mechanism either). This earthquake generated a tsunami that traveled across the Pacific.
  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely. I suspect that the fault that ruptured is eastward vergent (dipping to the west), so the west dipping nodal plane is probably the primary fault plane. However, this region of Kamchatka has numerous upper plate thrust and reverse faults (so the primary fault plane could be the other one, dipping to the east).
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault. The hypocentral depth of the M 5.5 plots this close to the location of the fault as mapped by Hayes et al. (2012).

    I include some inset figures in the poster.

  • In the upper left corner I include a map that shows the tectonic setting of this region, with the major plate boundary faults and volcanic arc designated by triangles (Bindeman et al., 2002). I placed an orange circle in the general location of the M 6.6 earthquake (sized relative to the magnitude range in the main map). Note the reverse fault mapped to the northeast of the epicenter. This fault dips to the southeast, supporting the east dipping solution for the M 6.6 moment tensor. I post this figure and their figure caption below.
  • To the right of this figure, I include a figure from Portnyagin and Manea (2008 ) that shows a low angle oblique view of the downgoing Pacific plate slab. I post this figure and their figure caption below.
  • In the lower left corner I include a map from the USGS Open File report (Rhea et al., 2010) that explains the historic seismicity for this region. I also plot the epicenter (orange dot).
  • In the upper right corner I include a map that shows more details about the faulting in the region. I place the epicenter for the M 6.6 as an orange circle. The location of the cross section I-I’ (plotted in the lower right corner) is designated by a dashed purple line.


  • Here is the tectonic map from Bindeman et al., 2002. The original figure caption is below in blockquote.

  • Tectonic setting of the Sredinny and Ganal Massifs in Kamchatka. Kamchatka/Aleutian junction is modified after Gaedicke et al. (2000). Onland geology is after Bogdanov and Khain (2000). 1, Active volcanoes (a) and Holocene monogenic vents (b). 2, Trench (a) and pull-apart basin in the Aleutian transform zone (b). 3, Thrust (a) and normal (b) faults. 4, Strike-slip faults. 5–6, Sredinny Massif. 5, Amphibolite-grade felsic paragneisses of the Kolpakovskaya series. 6, Allochthonous metasedimentary and metavolcanic rocks of the Malkinskaya series. 7, The Kvakhona arc. 8, Amphibolites and gabbro (solid circle) of the Ganal Massif. Lower inset shows the global position of Kamchatka. Upper inset shows main Cretaceous-Eocene tectonic units (Bogdanov and Khain 2000): Western Kamchatka (WK) composite unit including the Sredinny Massif, the Kvakhona arc, and the thick pile of Upper Cretaceous marine clastic rocks; Eastern Kamchatka (EK) arc, and Eastern Peninsulas terranes (EPT). Eastern Kamchatka is also known as the Olyutorka-Kamchatka arc (Nokleberg et al. 1998) or the Achaivayam-Valaginskaya arc (Konstantinovskaya 2000), while Eastern Peninsulas terranes are also called Kronotskaya arc (Levashova et al. 2000).

  • This map shows the configuration of the subducting slab. The original figure caption is below in blockquote.

  • Kamchatka subduction zone. A: Major geologic structures at the Kamchatka–Aleutian Arc junction. Thin dashed lines show isodepths to subducting Pacific plate (Gorbatov et al., 1997). Inset illustrates major volcanic zones in Kamchatka: EVB—Eastern Volcanic Belt; CKD—Central
    Kamchatka Depression (rift-like tectonic structure, which accommodates the northern end of EVB); SR—Sredinny Range. Distribution of Quaternary volcanic rocks in EVB and SR is shown in orange and green, respectively. Small dots are active vol canoes. Large circles denote CKD volcanoes: T—Tolbachik; K l — K l y u c h e v s k o y ; Z—Zarechny; Kh—Kharchinsky; Sh—Shiveluch; Shs—Shisheisky Complex; N—Nachikinsky. Location of profiles shown in Figures 2 and 3 is indicated. B: Three dimensional visualization of the Kamchatka subduction zone from the north. Surface relief is shown as semi-transparent layer. Labeled dashed lines and color (blue to red) gradation of subducting plate denote depths to the plate from the earth surface (in km). Bold arrow shows direction of Pacific Plate movement.

  • Here is the more detailed tectonic map from Konstantinovskaia et al. (2001).


  • This is the cross section associated with the above map.


  • Here is the Rhea et al. (2010) poster.

References:

  • Bindeman, I.N., Vinogradov, V.I., Valley, J.W., Wooden, J.L., and Natal’in, B.A., 2002. Archean Protolith and Accretion of Crust in Kamchatka: SHRIMP Dating of Zircons from Sredinny and Ganal Massifs in The Journal of Geology, v. 110, p. 271-289.
  • Hayes, G. P., D. J. Wald, and R. L. Johnson (2012), Slab1.0: A three-dimensional model of global subduction zone geometries, J. Geophys. Res., 117, B01302, doi:10.1029/2011JB008524.
  • Konstantinovskaia, E.A., 2001. Arc-continent collision and subduction reversal in the Cenozoic evolution of the Northwest Pacific: and example from Kamchatka (NE Russia) in Tectonophysics, v. 333, p. 75-94.
  • Koulakov, I.Y., Dobretsov, N.L., Bushenkova, N.A., and Yakovlev, A.V., 2011. Slab shape in subduction zones beneath the Kurile–Kamchatka and Aleutian arcs based on regional tomography results in Russian Geology and Geophysics, v. 52, p. 650-667.
  • Krutikov, L., et al., 2008. Active Tectonics and Seismic Potential of Alaska, Geophysical Monograph Series 179, doi:10.1029/179GM07
  • Lay, T., H. Kanamori, C. J. Ammon, A. R. Hutko, K. Furlong, and L. Rivera, 2009. The 2006 – 2007 Kuril Islands great earthquake sequence in J. Geophys. Res., 114, B11308, doi:10.1029/2008JB006280.
  • Portnyagin, M. and Manea, V.C., 2008. Mantle temperature control on composition of arc magmas along the Central Kamchatka Depression in Geology, v. 36, no. 7, p. 519-522.
  • Rhea, Susan, Tarr, A.C., Hayes, Gavin, Villaseñor, Antonio, Furlong, K.P., and Benz, H.M., 2010, Seismicity of the earth 1900–2007, Kuril-Kamchatka arc and vicinity: U.S. Geological Survey Open-File Report 2010–1083-C, scale 1:5,000,000.

Good Friday Earthquake 27 March 1964

In March of 1964, plate tectonics was still a hotly debated topic at scientific meetings worldwide. Some people still do not accept this theory (some Russian geologists favor alternative hypotheses; Shevchenko et al., 2006). At the time, there was some debate about whether the M 9.2 earthquake (the 2nd largest earthquake recorded with modern seismometers) was from a strike-slip or from a revers/thrust earthquake. Plafker and his colleagues found the evidence to put that debate to rest (see USGS video below).

I have prepared a new map showing the 1964 earthquake in context to the plate boundary using the same methods I have been using for my other earthquake reports. I also found a focal mechanism for this M 9.2 earthquake and included this on the map (Stauder and Bollinger, 1966).

    Here is the USGS website for this earthquake.

  • 1964.03.27 M 9.2

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include a focal mechanism for the M 9.2 earthquake determined by Stauder and Bollinger (1966). I include the USGS epicenters for earthquakes with magnitudes M ≥ 7.0.

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely.
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault. The hypocentral depth of the M 5.5 plots this close to the location of the fault as mapped by Hayes et al. (2012).

    I include some inset figures in the poster.

  • In the upper left corner I include two maps from the USGS, both using the MMI scale of shaking intensity mentioned above. The map on the left is the USGS Shakemap. This is a map that shows an estimate of how strongly the ground would shake during this earthquake. This is based upon a numerical model using Ground Motion Prediction Equations (GMPE), which are empirical relations between fault types, earthquake magnitude, distance from the fault, and shaking intensity. The map on the right is based upon peoples’ direct observations. Below each map are plots that show how these models demonstrate that the MMI attenuates (diminishes) with distance. The lines are the empirical relations. The dots are the data points.
  • To the right of those maps and figures is a map produced by Dr. Peter Haeussler from the USGS Alaska Science Center (pheuslr at usgs.gov) that shows the historic earthquakes along the Aleutian-Alaska subduction zone.
  • In the lower right corner I include an inset map from the USGS Seismicity History poster for this region (Benz et al., 2011). There is one seismicity cross section with its locations plotted on the map. The USGS plot these hypocenters along this cross section E-E’ (in green).
  • In the upper right corner, I include a figure that shows the measurements of uplift and subsidence observed by Plafker and his colleagues following the earthquake (Plafker, 1969). This is shown in map view and as a cross section.


Below is an educational video from the USGS that presents material about subduction zones and the 1964 earthquake and tsunami in particular.
Youtube Source IRIS

mp4 file for downloading.

    Credits:

  • Animation & graphics by Jenda Johnson, geologist
  • Directed by Robert F. Butler, University of Portland
  • U.S. Geological Survey consultants: Robert C. Witter, Alaska Science Center Peter J. Haeussler, Alaska Science Center
  • Narrated by Roger Groom, Mount Tabor Middle School

This is a map from Haeussler et al. (2014). The region in red shows the area that subsided and the area in blue shows the region that uplifted during the earthquake. These regions were originally measured in the field by George Plafker and published in several documents, including this USGS Professional Paper (Plafker, 1969).


Here is a cross section showing the differences of vertical deformation between the coseismic (during the earthquake) and interseismic (between earthquakes).


This figure, from Atwater et al. (2005) shows the earthquake deformation cycle and includes the aspect that the uplift deformation of the seafloor can cause a tsunami.


Here is a figure recently published in the 5th International Conference of IGCP 588 by the Division of Geological and Geophysical Surveys, Dept. of Natural Resources, State of Alaska (State of Alaska, 2015). This is derived from a figure published originally by Plafker (1969). There is a cross section included that shows how the slip was distributed along upper plate faults (e.g. the Patton Bay and Middleton Island faults).


Here is a graphic showing the sediment-stratigraphic evidence of earthquakes in Cascadia, but the analogy works for Alaska also. Atwater et al., 2005. There are 3 panels on the left, showing times of (1) prior to earthquake, (2) several years following the earthquake, and (3) centuries after the earthquake. Before the earthquake, the ground is sufficiently above sea level that trees can grow without fear of being inundated with salt water. During the earthquake, the ground subsides (lowers) so that the area is now inundated during high tides. The salt water kills the trees and other plants. Tidal sediment (like mud) starts to be deposited above the pre-earthquake ground surface. This sediment has organisms within it that reflect the tidal environment. Eventually, the sediment builds up and the crust deforms interseismically until the ground surface is again above sea level. Now plants that can survive in this environment start growing again. There are stumps and tree snags that were rooted in the pre-earthquake soil that can be used to estimate the age of the earthquake using radiocarbon age determinations. The tree snags form “ghost forests.


This is a photo that I took along the Seward HWY 1, that runs east of Anchorage along the Turnagain Arm. I attended the 2014 Seismological Society of America Meeting that was located in Anchorage to commemorate the anniversary of the Good Friday Earthquake. This is a ghost forest of trees that perished as a result of coseismic subsidence during the earthquake. Copyright Jason R. Patton (2014). This region subsided coseismically during the 1964 earthquake. Here are some photos from the paleoseismology field trip. (Please contact me for a higher resolution version of this image: quakejay at gmail.com)


Here is the USGS shakemap for this earthquake. The USGS used a fault model, delineated as black rectangles, to model ground shaking at the surface. The color scale refers to the Modified Mercalli Intensity scale, shown at the bottom.


http://earthjay.com/earthquakes/19640327_alaska/intensity.jpg

There is a great USGS Open File Report that summarizes the tectonics of Alaska and the Aleutian Islands (Benz et al., 2011). I include a section of their poster here. Below is the map legend.




Most recently, there was an earthquake along the Alaska Peninsula, a M 7.1 on 2016.01.24. Here is my earthquake report for this earthquake. Here is a map for the earthquakes of magnitude greater than or equal to M 7.0 between 1900 and today. This is the USGS query that I used to make this map. One may locate the USGS web pages for all the earthquakes on this map by following that link.


Here is an interesting map from Atwater et al., 2001. This figure shows how the estuarine setting in Portage, Alaska (along Turnagain Arm, southeast of Anchorage) had recovered its ground surface elevation in a short time following the earthquake. Within a decade, the region that had coseismically subsided was supporting a meadow with shrubs. By 1980, a spruce tree was growing here. This recovery was largely due to sedimentation, but an unreconciled amount of postseismic tectonic uplift contributed also. I include their figure caption as a blockquote.


(A and B) Tectonic setting of the 1964 Alaska earthquake. Subsidence from Plafker (1969). (C) Postearthquake deposits and their geologic setting in the early 1970s. (D–F) Area around Portage outlined in C, showing the landscape two years before the earthquake (D), two years after the earthquake (E), and nine years after the earthquake (F). In F, location of benchmark P 73 is from http://www.ngs.noaa.gov/cgi-bin/ds2.prl and the Seward (D-6) SE 7.5-minute quadrangle, provisional edition of 1984.

    Here is the tsunami forecast animation from the National Tsunami Warning Center. Below the animation, I include their caption as a blockquote. This includes information about the earthquake and the formation of the warning center.

  • Here is a link to the file for the embedded video below (22 MB 720 mp4)
  • Here is a link to the higher resolution file for the embedded video below (44 MB 1080 mp4)
  • At 5:36 pm on Friday, March 27, 1964 (28 March, 03:36Z UTC) the largest earthquake ever measured in North America, and the second-largest recorded anywhere, struck 40 miles west of Valdez, Alaska in Prince William Sound with a moment magnitude we now know to be 9.2. Almost an hour and a half later the Honolulu Magnetic and Seismic Observatory (later renamed the Pacific Tsunami Warning Center, or PTWC) was able to issue its first “tidal wave advisory” that noted that a tsunami was possible and that it could arrive in the Hawaiian Islands five hours later. Upon learning of a tsunami observation in Kodiak Island, Alaska, an hour and a half later the Honolulu Observatory issued a formal “tidal wave/seismic sea-wave warning” cautioning that damage was possible in Hawaii and throughout the Pacific Ocean but that it was not possible to predict the intensity of the tsunami. The earthquake did in fact generate a tsunami that killed 124 people (106 in Alaska, 13 in California, and 5 in Oregon) and caused about $2.3 billion (2016 dollars) in property loss all along the Pacific coast of North America from Alaska to southern California and in Hawaii. The greatest wave heights were in Alaska at over 67 m or 220 ft. and waves almost 10 m or 32 ft high struck British Columbia, Canada. In the “lower 48” waves as high as 4.5 m or 15 ft. struck Washington, as high as 3.7 m or 12 ft. struck Oregon, and as high as 4.8 m or over 15 ft. struck California. Waves of similar size struck Hawaii at nearly 5 m or over 16 ft. high. Waves over 1 m or 3 ft. high also struck Mexico, Chile, and even New Zealand.
  • As part of its response to this event the United States government created a second tsunami warning facility in 1967 at the Palmer Observatory, Alaska–now called the National Tsunami Warning Center (NTWC, http://ntwc.arh.noaa.gov/ )–to help mitigate future tsunami threats to Alaska, Canada, and the U.S. Mainland.
  • Today, more than 50 years since the Great Alaska Earthquake, PTWC and NTWC issue tsunami warnings in minutes, not hours, after a major earthquake occurs, and will also forecast how large any resulting tsunami will be as it is still crossing the ocean. PTWC can also create an animation of a historical tsunami with the same tool that it uses to determine tsunami hazards in real time for any tsunami today: the Real-Time Forecasting of Tsunamis (RIFT) forecast model. The RIFT model takes earthquake information as input and calculates how the waves move through the world’s oceans, predicting their speed, wavelength, and amplitude. This animation shows these values through the simulated motion of the waves and as they travel through the world’s oceans one can also see the distance between successive wave crests (wavelength) as well as their height (half-amplitude) indicated by their color. More importantly, the model also shows what happens when these tsunami waves strike land, the very information that PTWC needs to issue tsunami hazard guidance for impacted coastlines. From the beginning the animation shows all coastlines covered by colored points. These are initially a blue color like the undisturbed ocean to indicate normal sea level, but as the tsunami waves reach them they will change color to represent the height of the waves coming ashore, and often these values are higher than they were in the deeper waters offshore. The color scheme is based on PTWC’s warning criteria, with blue-to-green representing no hazard (less than 30 cm or ~1 ft.), yellow-to-orange indicating low hazard with a stay-off-the-beach recommendation (30 to 100 cm or ~1 to 3 ft.), light red-to-bright red indicating significant hazard requiring evacuation (1 to 3 m or ~3 to 10 ft.), and dark red indicating a severe hazard possibly requiring a second-tier evacuation (greater than 3 m or ~10 ft.).
  • Toward the end of this simulated 24 hours of activity the wave animation will transition to the “energy map” of a mathematical surface representing the maximum rise in sea-level on the open ocean caused by the tsunami, a pattern that indicates that the kinetic energy of the tsunami was not distributed evenly across the oceans but instead forms a highly directional “beam” such that the tsunami was far more severe in the middle of the “beam” of energy than on its sides. This pattern also generally correlates to the coastal impacts; note how those coastlines directly in the “beam” are hit by larger waves than those to either side of it.

References:

Earthquake Report: Bougainville and Solomons

In the past couple of weeks, there were some earthquakes in the equatorial western Pacific.

    Here are the USGS websites for the two largest earthquakes in the two regions.

  • 2017.03.04 M 6.1
  • 2017.03.19 M 6.0

The M 6.1 was a well behaved subduction zone earthquake associated with subduction of the Solomon Sea plate beneath the Solomon Islands, an island arc formed between opposing subduction zones. The M 6.0 earthquake and related earthquakes are more interesting. Baldwin et al. (2012) and Holm et al. (2016) both consider the North Solomon Trench (NST) to be inactive subduction zone, while the South Solomon Trench (SST) is an active subduction zone fault. There were 9 earthquakes larger than magnitude 2.5. The depths ranged between 4.2 and 53 km. There does not appear to be any direction that favors a deepening trend (i.e. getting deeper in the direction of the downgoing subduction zone fault). It is possible that these depths are not well constrained (the M 6.0 has a 4.2 km depth). These NNT earthquakes are either related to the NST, related to crustal faults, or reactivation of NST fault structures. Regardless, these are some interesting and mysterious earthquakes.

Below is my interpretive poster for this earthquake.


I plot the seismicity from the past year, with color representing depth and diameter representing magnitude (see legend).

  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely. The moment tensor shows northeast-southwest compression, perpendicular to the convergence at these plate boundaries.
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault.
  • I include some inset figures.

  • In the upper right corner is a map from Holm et al. (2016) that shows various key elements of the plate configuration in this region. Colors representing the age of the crust shows how the Woodlark Basin has an active spreading ridge system with ridges striking East-West. Large Igneous Provinces are shown here in dark gray. A relevant cross section is C-D, with the location highlighted in magenta on the map. The cross section shows how the North Solomon Trench is configured.
  • In the lower left corner is a generalized tectonic map of the region from Holm et al., 2016. This map shows the major plate boundary faults. Active subduction zones have shaded triangle fault symbols, while inactive subduction zones have un-shaded triangle fault line symbols.
  • In the upper left corner a figure from Oregon State University, which are based upon Hamilton (1979). “Tectonic microplates of the Melanesian region. Arrows show net plate motion relative to the Australian Plate.” This is from Johnson, 1976.
  • In the lower right corner is a figure from Baldwin et al. (2012). This figure shows a series of cross sections along this convergent plate boundary from the Solomon Islands in the east to Papua New Guinea in the west. Cross sections ‘A’ and ‘B’ are the most representative for the earthquakes today. I placed the general location of these cross sections on the main map as dashed orange lines. I present the map and this figure again below, with their original captions.


  • The M 6.1 is probably related to a series of earthquakes, including two M 7.9 megathrust events, further north along the San Crostobal Trench. Below are two earthquake interpretive posters for these earthquakes. This first poster is from my earthquake report on 2017.01.22.

  • In earlier earthquake reports, I discussed seismicity from 2000-2015 here. The seismicity on the west of this region appears aligned with north-south shortening along the New Britain trench, while seismicity on the east of this region appears aligned with more east-west shortening. Here is a map that I put together where I show these two tectonic domains with the seismicity from this time period (today’s earthquakes are not plotted on this map, but one may see where they might plot).


Background Figures

  • This figure shows an interpretation of the regional tectonics (Holm et al., 2016). I include the figure caption below as a blockquote.

  • Tectonic setting of Papua New Guinea and Solomon Islands. a) Regional plate boundaries and tectonic elements. Light grey shading illustrates bathymetry b 2000 m below sea level indicative of continental or arc crust, and oceanic plateaus; 1000 m depth contour is also shown. Adelbert Terrane (AT); Bismarck Sea fault (BSF); Bundi fault zone (BFZ); Feni Deep (FD); Finisterre Terrane (FT); Gazelle Peninsula (GP); Kia-Kaipito-Korigole fault zone (KKKF); Lagaip fault zone (LFZ); Mamberamo thrust belt (MTB); Manus Island (MI); New Britain (NB); New Ireland (NI); North Sepik arc (NSA); Ramu-Markham fault (RMF); Weitin Fault (WF);West Bismarck fault (WBF); Willaumez-Manus Rise (WMR).

  • This figure shows details of the regional tectonics (Holm et al., 2016). I include the figure caption below as a blockquote.

  • a) Present day tectonic features of the Papua New Guinea and Solomon Islands region as shown in plate reconstructions. Sea floor magnetic anomalies are shown for the Caroline plate (Gaina and Müller, 2007), Solomon Sea plate (Gaina and Müller, 2007) and Coral Sea (Weissel and Watts, 1979). Outline of the reconstructed Solomon Sea slab (SSP) and Vanuatu slab (VS)models are as indicated. b) Cross-sections related to the present day tectonic setting. Section locations are as indicated. Bismarck Sea fault (BSF); Feni Deep (FD); Louisiade Plateau
    (LP); Manus Basin (MB); New Britain trench (NBT); North Bismarck microplate (NBP); North Solomon trench (NST); Ontong Java Plateau (OJP); Ramu-Markham fault (RMF); San Cristobal trench (SCT); Solomon Sea plate (SSP); South Bismarck microplate (SBP); Trobriand trough (TT); projected Vanuatu slab (VS); West Bismarck fault (WBF); West Torres Plateau (WTP); Woodlark Basin (WB).

  • This map shows plate velocities and euler poles for different blocks. Note the counterclockwise motion of the plate that underlies the Solomon Sea (Baldwin et al., 2012). I include the figure caption below as a blockquote.

  • Tectonic maps of the New Guinea region. (a) Seismicity, volcanoes, and plate motion vectors. Plate motion vectors relative to the Australian plate are surface velocity models based on GPS data, fault slip rates, and earthquake focal mechanisms (UNAVCO, http://jules.unavco.org/Voyager/Earth). Earthquake data are sourced from the International Seismological Center EHB Bulletin (http://www.isc.ac.uk); data represent events from January 1994 through January 2009 with constrained focal depths. Background image is generated from http://www.geomapapp.org. Abbreviations: AB, Arafura Basin; AT, Aure Trough; AyT, Ayu Trough; BA, Banda arc; BSSL, Bismarck Sea seismic lineation; BH, Bird’s Head; BT, Banda Trench; BTFZ, Bewani-Torricelli fault zone; DD, Dayman Dome; DEI, D’Entrecasteaux Islands; FP, Fly Platform; GOP, Gulf of Papua; HP, Huon peninsula; LA, Louisiade Archipelago; LFZ, Lowlands fault zone; MaT, Manus Trench; ML, Mt. Lamington; MT, Mt. Trafalgar; MuT, Mussau Trough; MV, Mt. Victory; MTB, Mamberamo thrust belt; MVF, Managalase Plateau volcanic field; NBT, New Britain Trench; NBA, New Britain arc; NF, Nubara fault; NGT, New Guinea Trench; OJP, Ontong Java Plateau; OSF, Owen Stanley fault zone; PFTB, Papuan fold-and-thrust belt; PP, Papuan peninsula; PRi, Pocklington Rise; PT, Pocklington Trough; RMF, Ramu-Markham fault; SST, South Solomons Trench; SA, Solomon arc; SFZ, Sorong fault zone; ST, Seram Trench; TFZ, Tarera-Aiduna fault zone; TJ, AUS-WDKPAC triple junction; TL, Tasman line; TT, Trobriand Trough;WD, Weber Deep;WB, Woodlark Basin;WFTB, Western (Irian) fold-and-thrust belt; WR,Woodlark Rift; WRi, Woodlark Rise; WTB, Weyland thrust; YFZ, Yapen fault zone.White box indicates the location shown in Figure 3. (b) Map of plates, microplates, and tectonic blocks and elements of the New Guinea region. Tectonic elements modified after Hill & Hall (2003). Abbreviations: ADB, Adelbert block; AOB, April ultramafics; AUS, Australian plate; BHB, Bird’s Head block; CM, Cyclops Mountains; CWB, Cendrawasih block; CAR, Caroline microplate; EMD, Ertsberg Mining District; FA, Finisterre arc; IOB, Irian ophiolite belt; KBB, Kubor & Bena blocks (including Bena Bena terrane); LFTB, Lengguru fold-and-thrust belt; MA, Mapenduma anticline; MB, Mamberamo Basin block; MO, Marum ophiolite belt; MHS, Manus hotspot; NBS, North Bismarck plate; NGH, New Guinea highlands block; NNG, Northern New Guinea block; OKT, Ok Tedi mining district; PAC, Pacific plate; PIC, Porgera intrusive complex; PSP, Philippine Sea plate; PUB, Papuan Ultramafic Belt ophiolite; SB, Sepik Basin block; SDB, Sunda block; SBS, South Bismarck plate; SIB, Solomon Islands block; WP, Wandamen peninsula; WDK, Woodlark microplate; YQ, Yeleme quarries.

  • This figure incorporates cross sections and map views of various parts of the regional tectonics (Baldwin et al., 2012). The New Britain region is in the map near the A and B sections. I include the figure caption below as a blockquote.

  • Oblique block diagram of New Guinea from the northeast with schematic cross sections showing the present-day plate tectonic setting. Digital elevation model was generated from http://www.geomapapp.org. Oceanic crust in tectonic cross sections is shown by thick black-and-white hatched lines, with arrows indicating active subduction; thick gray-and-white hatched lines indicate uncertain former subduction. Continental crust, transitional continental crust, and arc-related crust are shown without pattern. Representative geologic cross sections across parts of slices C and D are marked with transparent red ovals and within slices B and E are shown by dotted lines. (i ) Cross section of the Papuan peninsula and D’Entrecasteaux Islands modified from Little et al. (2011), showing the obducted ophiolite belt due to collision of the Australian (AUS) plate with an arc in the Paleogene, with later Pliocene extension and exhumation to form the D’Entrecasteaux Islands. (ii ) Cross section of the Papuan peninsula after Davies & Jaques (1984) shows the Papuan ophiolite thrust over metamorphic rocks of AUS margin affinity. (iii ) Across the Papuan mainland, the cross section after Crowhurst et al. (1996) shows the obducted Marum ophiolite and complex folding and thrusting due to collision of the Melanesian arc (the Adelbert, Finisterre, and Huon blocks) in the Late Miocene to recent. (iv) Across the Bird’s Head, the cross section after Bailly et al. (2009) illustrates deformation in the Lengguru fold-and-thrust belt as a result of Late Miocene–Early Pliocene northeast-southwest shortening, followed by Late Pliocene–Quaternary extension. Abbreviations as in Figure 2, in addition to NI, New Ireland; SI, Solomon Islands; SS, Solomon Sea; (U)HP, (ultra)high-pressure.


Background Videos

  • I put together an animation that shows the seismicity for this region from 1900-2016 for earthquakes with magnitude of M ≥ 6.5. Here is the USGS query I used to search for these data used in this animation. I show the location of the Benz et al. (2011) cross section E-E’ as a yellow line on the main map.
  • In this animation, I include some figures from the interpretive poster above. I also include a tectonic map based upon Hamilton (1979). Music is from copyright free music online and is entitled “Sub Strut.” Above the video I present a map showing all the earthquakes presented in the video.
  • Here is a link to the embedded video below (5 MB mp4)


References:

Course Material and Educational Resources