Earthquake Report: Mendocino triple junction

Well, I was on the road for 1.5 days (work party for the Community Village at the Oregon Country Fair). As I was driving home, there was a magnitude M 5.6 earthquake in coastal northern California.
https://earthquake.usgs.gov/earthquakes/eventpage/nc73201181/executive
I didn’t realize this until I was almost home (finally hit the sack around 4 am).
This earthquake follows a sequence of quakes further to the northwest, however their timing is merely a coincidence. Let me repeat this. The M 5.6 earthquake is not related to the sequence of earthquakes along the Blanco fracture zone.
Contrary to what people have posted on social media, there was but a single earthquake. This earthquake happened beneath the area of Petrolia, nearby the 1991 Honeydew Earthquake. More about the Honeydew Earthquake can be found here.
This region also had a good sized shaker in 1992, the Cape Mendocino Earthquake, which led to the development of the National Tsunami Hazard Mitigation Program. More about the Cape Mendocino Earthquake can be found on the 25th anniversary page here and in my earthquake report here.
The regional tectonics in coastal northern California are dominated by the Pacific-North America plate boundary. North of Cape Mendocino, this plate boundary is convergent and forms the Cascadia subduction zone (CSZ). To the south of Cape Mendocino, the plate boundary is the right-lateral (dextral) San Andreas fault (SAF). Where these 2 fault systems meet, there is another plate boundary system, the right-lateral strike-slip Mendocino fault (don’t write Mendocino fracture zone on your maps!). Where these 3 systems meet is called the Mendocino triple junction (MTJ).
The MTJ is a complicated region as these plate boundaries overlap in ways that we still do not fully understand. Geologic mapping in the mid- to late-20th century provides some basic understanding of the long term history. However, recent discoveries have proven that this early work needs to be revisited as there are many unanswered questions (and some of this early work has been demonstrated to be incorrect). Long live science!
Last night’s M 5.6 temblor happened where one strand of the MF trends onshore (another strand bends towards the south). But, it also is where the SAF trends onshore. At this point, I am associating this earthquake with the MF (so, a right-lateral strike-slip earthquake). The mechanism suggest that this is not a SAF related earthquake. However, it is oriented in a way that it could be in the Gorda plate (making it a left-lateral strike-slip earthquake). However, this quake is at the southern edge of the Gorda plate (sedge), so it is unlikely this is a Gorda plate event.

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 earthquake epicenters from 1918-2018 with magnitudes M ≥ 5.0 in one version.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.

  • 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 2.0 contours plotted (Hayes, 2018), 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.

    Magnetic Anomalies

  • In the map below, I include a transparent overlay of the magnetic anomaly data from EMAG2 (Meyer et al., 2017). As oceanic crust is formed, it inherits the magnetic field at the time. At different points through time, the magnetic polarity (north vs. south) flips, the North Pole becomes the South Pole. These changes in polarity can be seen when measuring the magnetic field above oceanic plates. This is one of the fundamental evidences for plate spreading at oceanic spreading ridges (like the Gorda rise).
  • Regions with magnetic fields aligned like today’s magnetic polarity are colored red in the EMAG2 data, while reversed polarity regions are colored blue. Regions of intermediate magnetic field are colored light purple.
  • We can see the roughly ~north-south trends of these red and blue stripes in the Pacific plate. These lines are parallel to the ocean spreading ridges from where they were formed. The stripes disappear at the subduction zone because the oceanic crust with these anomalies is diving deep beneath the North America plate, so the magnetic anomalies from the overlying Sunda plate mask the evidence for the Juan de Fuca and Gorda plates.

    Global Strain

  • In a map below, I include a transparent overlay of the Global Strain Rate Map (Kreemer et al., 2014).
  • The mission of the Global Strain Rate Map (GSRM) project is to determine a globally self-consistent strain rate and velocity field model, consistent with geodetic and geologic field observations. The overall mission also includes:
    1. contributions of global, regional, and local models by individual researchers
    2. archive existing data sets of geologic, geodetic, and seismic information that can contribute toward a greater understanding of strain phenomena
    3. archive existing methods for modeling strain rates and strain transients
  • The completed global strain rate map will provide a large amount of information that is vital for our understanding of continental dynamics and for the quantification of seismic hazards.
  • The version used in the poster(s) below is an update to the original 2004 map (Kreemer et al., 2000, 2003; Holt et al., 2005).

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

  • n 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 lower right corner is a map that shows a comparison between the USGS Did You Feel It? reports and the USGS Modified Mercalli Intensity shakemap model. This comparison shows that the model is a decent fit for the reports from real people. If you felt the earthquake, please submit a report to the USGS here.
  • In the upper right corner I include a larger scale view of seismicity for this area. I highlight the important historic events (e.g. the 1991 Honeydew Earthquake and the 1992 Cape Mendocino Earthquake sequence.
  • Here is the map with a month’s seismicity plotted.

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

  • Here is the map with a century’s seismicity plotted along with the Global Strain Map with a 30% transparency.

  • Here is the educational interpretive poster from the 1992 Cape Mendocino Earthquake (report here).

  • The USGS has been increasing the list of products that are produced in association with their earthquake pages. One of these products is an earthquake forecast (not a prediction as nobody can predict earthquakes yet) that lists the chance of an earthquake with a given magnitude over a certain period of time. The forecast for the M 5.6 earthquake is found here. These forecasts are updated periodically, so the information will change with time. Below is a table where I present the forecast as it was when I checked the page this morning (would be nice if the USGS would produce an easy to read table).
  • From the USGS:

    Be ready for more earthquakes

    • More earthquakes than usual (called aftershocks) will continue to occur near the mainshock.
    • When there are more earthquakes, the chance of a large earthquake is greater which means that the chance of damage is greater.
    • The USGS advises everyone to be aware of the possibility of aftershocks, especially when in or around vulnerable structures such as unreinforced masonry buildings.
    • This earthquake could be part of a sequence. An earthquake sequence may have larger and potentially damaging earthquakes in the future, so remember to: Drop, Cover, and Hold on.

    What we think will happen next

    • According to our forecast, over the next 1 Week there is a < 1 % chance of one or more aftershocks that are larger than magnitude 5.6. It is likely that there will be smaller earthquakes over the next 1 Week, with 0 to 11 magnitude 3 or higher aftershocks. Magnitude 3 and above are large enough to be felt near the epicenter. The number of aftershocks will drop off over time, but a large aftershock can increase the numbers again, temporarily.

    About our earthquake forecasts

    • No one can predict the exact time or place of any earthquake, including aftershocks. Our earthquake forecasts give us an understanding of the chances of having more earthquakes within a given time period in the affected area. We calculate this earthquake forecast using a statistical analysis based on past earthquakes.
    • Our forecast changes as time passes due to decline in the frequency of aftershocks, larger aftershocks that may trigger further earthquakes, and changes in forecast modeling based on the data collected for this earthquake sequence.


  • Gosh, almost forgot to include this photo of the seismic waves recorded on the Humboldt State University Department of Geology Baby Benioff seismometer. Photo Credit: Amanda Admire.

USGS Landslide and Liquefaction Ground Failure data products

  • Below I present a series of maps that are intended to address the excellent ‘new’ products included in the USGS earthquake pages: landslide probability and liquefaction susceptibility (a.k.a. the Ground Failure data products).
  • First I present the landslide probability model. This is a GIS data product that relates a variety of factors to the probability (the chance of) landslides as triggered by this earthquake. There are a number of assumptions that are made in order to be able to produce this model across such a large region, though this is still of great value (like other aspects from teh USGS, e.g. the PAGER alert). Learn more about all of these Ground Failure products here.
  • There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). I spend more time discussing landslides and liquefaction in this recent earthquake report.
  • This model, like all landslide computer models, uses similar inputs. I review these here:
    1. Some information about ground shaking. Often, people use Peak Ground Acceleration, though in the past decade+, it has been recognized that the parameter “Arias Intensity” is a better measure of the energy imparted by the earthquake across the land and seascape. Instead of simply accounting for the peak accelerations, AI integrates the entire energy (duration) during the earthquake. That being said, PGA is a more common parameter that is available for people to use. For example, when I was modeling slope stability for the 2004 Sumatra-Andaman subduction zone earthquake, the only model that was calibrated to observational data were in units of PGA. The first order control to shaking intensity (energy observed at any particular location) is distance to the earthquake fault that slipped.
    2. Some information about the strength of the materials (e.g. angle of internal friction (the strength) and cohesion (the resistance).
    3. Information about the slope. Steeper slopes, with all other things being equal, are more likely to fail than are shallower slopes. Think about skiing. Beginners (like me) often choose shallower slopes to ski because they will go down the slope slower, while experts choose steeper slopes.
  • Areas that are red are more likely to experience landslides than areas that are colored blue. I include a coarse resolution topographic/bathymetric dataset to help us identify where the mountains are relative to the coastal plain and continental shelf (submarine).

  • Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
  • Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
  • An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
  • The liquefaction susceptibility map for the M 5.6 earthquake did not suggest that there would be possibly much liquefaction from this earthquake (probably due to the small magnitude). I discuss liquefaction more in my earthquake report on the 28 September 20018 Sulawesi, Indonesia earthquake, landslide, and tsunami here.
  • Here is a map that shows shaking intensity using the MMI scale (mentioned and plotted in the main earthquake poster maps). I present this here in the same format as the ground failure model maps so we can compare these other maps with the ground shaking model (which is a first order control on slope failure).

Other Report Pages

Some Relevant Discussion and Figures

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

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


Geologic Fundamentals

  • For more on the graphical representation of moment tensors and focal mechanisms, check this IRIS video out:
  • Here is a fantastic infographic from Frisch et al. (2011). This figure shows some examples of earthquakes in different plate tectonic settings, and what their fault plane solutions are. There is a cross section showing these focal mechanisms for a thrust or reverse earthquake. The upper right corner includes my favorite figure of all time. This shows the first motion (up or down) for each of the four quadrants. This figure also shows how the amplitude of the seismic waves are greatest (generally) in the middle of the quadrant and decrease to zero at the nodal planes (the boundary of each quadrant).

  • Here is another way to look at these beach balls.
  • 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. The following three animations are from IRIS.
  • Strike Slip:

    Compressional:

    Extensional:

  • This is an image from the USGS that shows how, when an oceanic plate moves over a hotspot, the volcanoes formed over the hotspot form a series of volcanoes that increase in age in the direction of plate motion. The presumption is that the hotspot is stable and stays in one location. Torsvik et al. (2017) use various methods to evaluate why this is a false presumption for the Hawaii Hotspot.

  • A cutaway view along the Hawaiian island chain showing the inferred mantle plume that has fed the Hawaiian hot spot on the overriding Pacific Plate. The geologic ages of the oldest volcano on each island (Ma = millions of years ago) are progressively older to the northwest, consistent with the hot spot model for the origin of the Hawaiian Ridge-Emperor Seamount Chain. (Modified from image of Joel E. Robinson, USGS, in “This Dynamic Planet” map of Simkin and others, 2006.)

  • Here is a map from Torsvik et al. (2017) that shows the age of volcanic rocks at different locations along the Hawaii-Emperor Seamount Chain.

  • Hawaiian-Emperor Chain. White dots are the locations of radiometrically dated seamounts, atolls and islands, based on compilations of Doubrovine et al. and O’Connor et al. Features encircled with larger white circles are discussed in the text and Fig. 2. Marine gravity anomaly map is from Sandwell and Smith.

  • Here is a great tweet that discusses the different parts of a seismogram and how the internal structures of the Earth help control seismic waves as they propagate in the Earth.

    Social Media

    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.
  • Dengler, L.A., and McPherson, R.C., 1993. The 17 August 1991 Honeydew Earthquake, North Coast California: A Case for Revising the Modified Mercalli Scale in Sparsely Populated Areas in BSSA, v. 83, no. 4, pp. 1081-1094
  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release, https://doi.org/10.5066/F7PV6JNV.
  • Holt, W. E., C. Kreemer, A. J. Haines, L. Estey, C. Meertens, G. Blewitt, and D. Lavallee (2005), Project helps constrain continental dynamics and seismic hazards, Eos Trans. AGU, 86(41), 383–387, , https://doi.org/10.1029/2005EO410002. /li>
  • Kreemer, C., J. Haines, W. Holt, G. Blewitt, and D. Lavallee (2000), On the determination of a global strain rate model, Geophys. J. Int., 52(10), 765–770.
  • Kreemer, C., W. E. Holt, and A. J. Haines (2003), An integrated global model of present-day plate motions and plate boundary deformation, Geophys. J. Int., 154(1), 8–34, , https://doi.org/10.1046/j.1365-246X.2003.01917.x.
  • Kreemer, C., G. Blewitt, E.C. Klein, 2014. A geodetic plate motion and Global Strain Rate Model in Geochemistry, Geophysics, Geosystems, v. 15, p. 3849-3889, https://doi.org/10.1002/2014GC005407.
  • 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
  • Meyer, B., Saltus, R., Chulliat, a., 2017. EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. https://doi.org/10.7289/V5H70CVX
  • Müller, R.D., Sdrolias, M., Gaina, C. and Roest, W.R., 2008, Age spreading rates and spreading asymmetry of the world’s ocean crust in Geochemistry, Geophysics, Geosystems, 9, Q04006, https://doi.org/10.1029/2007GC001743
  • 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.

Return to the Earthquake Reports page.


Earthquake Report: Mendocino fault! (northern California)

I was driving around Eureka today, running to the appliance center to get an appliance (heheh). I got a message from a long time held friend (who lives in Salinas, CA). They asked me if I was OK, given that there was an earthquake up here. I thought I had not felt it because I was driving around. However, after looking at the USGS website, I learned the earthquake happened earlier, while I was back working on my house. The main reason I did not feel it is because it was too far away.

Once I got home, after work, I noticed that lots of people were discussing how they were confused about the earthquake notifications from the USGS. Apparently, there were two M 5.X earthquakes in the USGS earthquake online system for a while. Then there was one. This is a common occurrence and I prepared an explanation for some people Here is what I wrote for these people on social media:

this happens regularly. earthquake notifications are automatic as epicenter locations are identified from incoming seismic waves in the seismic network. sometimes the named arrivals (eg. p wave, s wave, and the many other arrivals) are miss-correlated between stations. this miss-correlation then leads to earthquakes in the database that are not real.

seismologists are monitoring the process and review these data for quality, looking for mistakes, and refining magnitude estimates, moment tensor and focal mechanism solutions, location estimates, casualty estimages (PAGER alerts), and all the derivative data products (intensity, PGA, PGV, etc. maps and data).

sometimes these earthquakes are from data in the same location as the real earthquake (like today) and sometimes they are “picked” from seismic data from remote earthquakes.
some of these earthquakes are listed here:

USGS Why do some earthquakes disappear from the map/list?

Today’s M 5.7 earthquake was along the western part of the Mendocino fault (MF), a right-lateral (dextral) transform plate boundary. This plate boundary connects the Gorda ridge and Juan de Fuca rise spreading centers with their counterparts in the Gulf of California, with the San Andreas strike-slip fault system. Transform plate boundaries are defined that they are strike-slip and that they connect spreading ridges. In this sense of the definition, the Mendocino fault and the San Andreas fault are part of the same system. Here is the USGS website for this earthquake.

See the figures from Rollins and Stein (2010) below. More on earthquakes in this region can be found in Earthquake Reports listed at the bottom of this page above the appendices.

The San Andreas fault is a right-lateral strike-slip transform plate boundary between the Pacific and North America plates. The plate boundary is composed of faults that are parallel to sub-parallel to the SAF and extend from the west coast of CA to the Wasatch fault (WF) system in central Utah (the WF runs through Salt Lake City and is expressed by the mountain range on the east side of the basin that Salt Lake City is built within).

The three main faults in the region north of San Francisco are the SAF, the MF, and the Bartlett Springs fault (BSF). I also place a graphical depiction of the USGS moment tensor for this earthquake. The SAF, MF, and BSF are all right lateral strike-slip fault systems. There are no active faults mapped in the region of Sunday’s epicenter, but I interpret this earthquake to have right-lateral slip. Without more seismicity or mapped faults to suggest otherwise, this is a reasonable interpretation.

The Cascadia subduction zone is a convergent plate boundary where the Juan de Fuca and Gorda plates (JDFP and GP, respectively) subduct norteastwardly beneath the North America plate at rates ranging from 29- to 45-mm/yr. The Juan de Fuca and Gorda plates are formed at the Juan de Fuca Ridge and Gorda Rise spreading centers respectively. More about the CSZ can be found here.

There was a good sized (M 6.5) MF earthquake late last year on 2016.12.08. I present my poster for that earthquake below. Here is my report for that earthquake. Here is the updated report.

Below I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I use the USGS Quaternary fault and fold database for the faults.

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

This is a preliminary report and I hope to prepare some updates as I collect more information.

    I have placed several inset figures.

  • In the upper right 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). I placed a blue star in the general location of today’s M 5.7 earthquake.
  • 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. Today’s earthquake did not occur along the CSZ, so did not produce crustal deformation like this. However, it is useful to know this when studying the CSZ.
  • In the lower 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. This map shows the major earthquakes that have occurred in this region, prior to their publication in 2010. I place a blue star in the general location of today’s earthquake.
  • Above the Rollins and Stein (2010) map are two illustrations showing the difference between a right-lateral and a left-lateral strike slip fault. This is from California Institute of Technology (Caltech).
  • To the right of the Rollins and Stein (2010) map, is a generalized illustration showing an interpretation of the results from these authors. They suggest that, for a variety of earthquake sources in this region, which types of faults have inhibited or promoted earthquake likelihood. The relevant part is C, which tests whether there is an increased or decreased likelihood (chance) of an earthquake on the left-lateral strike-slip faults in the Gorda plate. Based upon today’s M 5.7, there is a slight increase in the chance of a Gorda plate earthquake to the northwest of today’s M 5.7 earthquake. This is the distant side of the M 5.7 earthquake, so any potential GP earthquake would be further away.
  • In the upper right corner is a figure that many people in Humboldt and Del Norte counties might be interested in (the two most northwesterly counties in CA). These two panels both show the same general result (as relevant to this discussion), the increased or decreased chance of an earthquake on two types of faults (north of the dashed line, the chance on GP left-lateral faults; south of dashed line, the chance on the MF. The region of this figure is outlined in dashed white transparent box on the main poster. We can see that the CSZ is just to the east of this figure. People always want to know if there is an increased chance of a megathrust earthquake on the CSZ. This M 5.7 will not have a direct impact upon the CSZ. Over time, earthquakes like this actually bring the CSZ closer to an earthquake (they do not relieve stress, but increase it). But the deformation of the Gorda and Pacific plates is localized near the earthquake. So, it does not change the stress on the megathrust. But, hundreds of earthquakes like this, over time, do increase the stress on the megathrust.
  • The figure here helps us evaluate this concept for this M 5.7 earthquake. The 1994 earthquake, represented in this figure, caused an increase in stress along faults generally in the region of this figure (extending outwards more a little to the south, less more to the west, and very little more to the north and east. The take away is that the 1994 did not change the stress on faults very much in the region of the megathrust. Because today’s M 5.7 earthquake is even further to the west, there is not a possibility that this M 5.7 had any affect on the megathrust.


  • Here is the 2016.12.08 earthquake report poster from this report.

    • For more on the graphical representation of moment tensors and focal mechnisms, check this IRIS video out:
    • Here is a fantastic infographic from Frisch et al. (2011). This figure shows some examples of earthquakes in different plate tectonic settings, and what their fault plane solutions are. There is a cross section showing these focal mechanisms for a thrust or reverse earthquake. The upper right corner includes my favorite figure of all time. This shows the first motion (up or down) for each of the four quadrants. This figure also shows how the amplitude of the seismic waves are greatest (generally) in the middle of the quadrant and decrease to zero at the nodal planes (the boundary of each quadrant).

    • Here is a map from Rollins and Stein, showing their interpretations of different historic earthquakes in the region. This was published in response to the January 2010 Gorda plate earthquake. The faults are from Chaytor et al. (2004). The 1980, 1992, 1994, 2005, and 2010 earthquakes are plotted and labeled. I did not mention the 2010 earthquake, but it most likely was just like 1980 and 2005, a left-lateral strike-slip earthquake on a northeast striking fault.

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

    • Here is the Rollins and Stein (2010) figure that is in the report above. I include their figure caption as blockquote below.

    • Coulomb stress changes imparted by our models of (a) a bilateral rupture and (b) a unilateral eastward rupture for the 1994 Mw = 7.0 Mendocino Fault Zone earthquake to the epicenters of the 1995 Mw = 6.6 southern Gorda zone earthquake (N) and the 2000 Mw = 5.9 Mendocino Fault Zone earthquake (O). Calculation depth is 5 km.

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

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

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

    • The Gorda and Juan de Fuca plates subduct beneath the North America plate to form the Cascadia subduction zone fault system. In 1992 there was a swarm of earthquakes with the magnitude Mw 7.2 Mainshock on 4/25. Initially this earthquake was interpreted to have been on the Cascadia subduction zone (CSZ). The moment tensor shows a compressional mechanism. However the two largest aftershocks on 4/26/1992 (Mw 6.5 and Mw 6.7), had strike-slip moment tensors. These two aftershocks align on what may be the eastern extension of the Mendocino fault.
    • There have been several series of intra-plate earthquakes in the Gorda plate. Two main shocks that I plot of this type of earthquake are the 1980 (Mw 7.2) and 2005 (Mw 7.2) earthquakes. I place orange lines approximately where the faults are that ruptured in 1980 and 2005. These are also plotted in the Rollins and Stein (2010) figure above. The Gorda plate is being deformed due to compression between the Pacific plate to the south and the Juan de Fuca plate to the north. Due to this north-south compression, the plate is deforming internally so that normal faults that formed at the spreading center (the Gorda Rise) are reactivated as left-lateral strike-slip faults. In 2014, there was another swarm of left-lateral earthquakes in the Gorda plate. I posted some material about the Gorda plate setting on this page.
    • 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:

    • This figure shows what a transform plate boundary fault is. Looking down from outer space, the crust on either side of the fault moves side-by-side. When one is standing on the ground, on one side of the fault, looking across the fault as it moves… If the crust on the other side of the fault moves to the right, the fault is a “right lateral” strike slip fault. The Mendocino and San Andreas faults are right-lateral (dextral) strike-slip faults. I believe this is from Pearson Higher Ed.

    Update

    • Abbreviations
    • BSF – Bartlett Springs fault
    • CA – California
    • CSZ – Cascadia subduction zone
    • GP – Gorda plate
    • JDFP – Juan de Fuca plate
    • MF – Mendocino fault
    • MMI – Modified Mercalli Intensity Scale
    • SAF – San Andreas fault
    • USGS – U.S. Geological Survey
    • WF – Wasatch fault

      References

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

    Earthquake Report: Philippines

    I put these together earlier this week for my 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.

    • A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
    • Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.
    • I 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 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 Sibuyan Sea faults are “forearc sliver” faults that accommodate oblique convergence between the Philippine Sea and Sunda plates. The M 5.9 epicenter is designated by a red star. The Manila trench (MT; northernmost east-dipping subduction zone) and the Philippine trench (PT; easternmost west-dipping subduction zone) are major players here. The high velocity plate motion vectors (in mm/yr) show the relative plate motions and the possible source of hi frequency seismicity.
    • 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).


    • In the above poster, I choose 4 earthquakes for which to 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 the interpretive poster for the 08 April 2017 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:

      Basic & General References

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

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

    Return to the Earthquake Reports page.

    Earthquake Report: Botswana!

    This is a very interesting M 6.5 earthquake, which was preceded by a probably unrelated M 5.2 earthquake. Last September, there was an M 5.7 earthquake in Tanzania along the western shores of Lake Victoria. Here is my report for that 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: 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:

    Earthquake Report: New Britain

    We just had an earthquake with a USGS magnitude of M 6.5 along the subduction zone formed by the convergence of the Solomon Sea plate on the south and the South Bismarck plate on the north.
    Here is the USGS website for this M 6.5 earthquake.

    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 also include seismicity from 1917-2017 for earthquakes with magnitudes M ≥ 7.0. I include the USGS moment tensor from this earthquake.

    • 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 this plate boundary.
    • 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 plots this close to the location of the fault as mapped by Hayes et al. (2012). This earthquake appears to be in the upper plate based upon the hypocentral depth of <30 km and the slab contour of 60 km at this location, though this is unclear.

      I include some inset figures in the poster.

    • At the top of the poster are two figures 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.” To the right of the map is a cross section showing how the Solomon Sea plate is subducting beneath New Britain. This is from Johnson, 1976.
    • In the upper right corner are two figures from Holm and Richards (2013). Their paper discusses the back-arc spreading in the Bismarck Sea. They use hypocenter data to construct this 3-D model of the slab. On the right is a forecast of how the slab will be consumed along these subduction zones in the future.
    • In the lower right corner is a cross sectional view of the New Britain subduction zone for the past 10 million years (Holm and Richards, 2013). The approximate location of for this cross section is designated by an orange dashed line on the main map, with N and S labels.
    • In the upper left 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 section ‘C’ is the most representative for the earthquake today. I place the general location of the C-C’ section on the main map as an orange dashed line. I present the map and this figure again below, with their original captions.
    • In the lower left corner is the map from Baldwin et al. (2012) that goes along with the cross sections above.


    • This is a map showing the seismicity of this region since 2000 A.D.

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

    • Here is the generalized tectonic map of the region from Holm et al., 2015. I include the figure caption below as a blockquote.

    • Tectonic setting and mineral deposits of eastern Papua New Guinea and Solomon Islands. The modern arc setting related to formation of the mineral deposits comprises, from west to east, the West Bismarck arc, the New Britain arc, the Tabar-Lihir-Tanga-Feni Chain and the Solomon arc, associated with north-dipping subduction/underthrusting at the Ramu-Markham fault zone, New Britain trench and San Cristobal trench respectively. Arrows denote plate motion direction of the Australian and Pacific plates. Filled triangles denote active subduction. Outlined triangles denote slow or extinct subduction. NBP: North Bismarck plate; SBP: South Bismarck plate; AT: Adelbert Terrane; FT: Finisterre Terrane; RMF: Ramu-Markham fault zone; NBT: New Britain trench.

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

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

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

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

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

    References:

    Earthquake Report: Banda Sea

    Earlier this week there was a moderate earthquake along a strike-slip fault that appears to adjoin the Banda/Timor/Java Arc with New Guinea. This strike-slip fault appears to cross oblique to the subduction zone that forms the Timor Trench to the south and the Seram Trench to the north. Various researchers portray the faulting in this region differently. Given earthquake moment tensor and focal mechanism from the USGS, this earthquake supports the interpretation that this fault system is left-lateral (synistral) strike-slip. A focal mechanism from an earthquake in 1938 (magnitude M 8.5) provides evidence that is a little more confusing. But, this region is a complicated region.
    UPDATE: 2017.03.05 (23:00 local time): Interesting, I was reading my tweet after Lila Lisle noticed a mistake. I fixed that, but later realized that the possible fault in the downgoing plate is not the Sorog fault, but a fault that might intersect with the Sorong fault. I did not delete the tweet since this mistake is a topic of conversation and not really a key part of the story.
    Here is the USGS earthquake website for this M 5.5 earthquake.

    Here are the USGS websites for the major earthquakes in this region from the past century.

    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 also include seismicity from 1917-2017 for earthquakes with magnitudes M ≥ 7.0. I show the fault plane solutions for some of these earthquakes.

    • moment tensors for 2005, 2012, and 2017 (USGS)
    • focal mechanisms for 1987 (USGS) and 1938 (Okal and Reymond, 2003)
    • The fault plane solutions for the 1963, 1987, 2005, and 2012 earthquakes are all very similar to the 2017 M 5.5. However, these earthquakes are form two depth “populations.” The 1963 and 1987 earthquakes are at ~65 km depth, while the 2005, 2012, and 2017 are between 150-200 km. There are some earthquakes that are much shallower depth eastward along this possibly strike-slip fault. Near where this fault comes on land at the base of the Bird’s Head in New Guinea, there are some earthquakes from the past couple of decades that also have strike-slip fault–plane solutions. The earthquake depths along this Sorong fault (Hall, 2011) appear to show that the Sorong fault is active beneath the Sunda plate. The 1938 earthquake may be the result of some form of strain partitioned faulting(?). Alternatively, these earthquakes may be unrelated to the Sorong fault. There may be some internal structure in the Australia plate that is interacting with the subduction zones or other faults (some preexisting structure that is optimally oriented to reactivate with the .
    • 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 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 is a general tectonic map for this part of the world. I placed a green star in the location of this M 5.5 earthquake (Zahirovic et al., 2014).
    • In the upper right corner is a low-angle oblique view of the plate boundaries in the northern part of this region (Hall, 2011). The upper part of the diagram shows the opposing vergent subduction zones along that strike north-south along the Molucca Strait (Halmahera, Philippines). The lower panel shows the downgoing Australia plate along the Timor Trench and Seram Trench. Note the location of the Bird’s Head, the northwestern part of New Guinea. I have also labeled this region in the main map for comparison. The strike-slip fault at the northern boundary of New Guinea is the Sorong fault and this is labeled in this Hall (2011) figure.
    • Below the Hall (2011) figure is a figure from Baldwin et al. (2012) that shows the regional seismicity and faulting as they are related to different geologic types in New Guinea. I placed a green star in the location of this M 5.5 earthquake.
    • In the lower right corner is a part of the USGS Poster that reviews the seismicity of this region for the past century or so (Benz et a., 2011). The map shows seismicity with depth, along with some cross section locations. I place a green star at the location of this M 5.5 earthquake. I present three of the cross sections from this poster, A-A’, B-B’, and C-C.’ Of particular interest is the section B-B’ because this is placed near the M 5.5 earthquake. I have placed a green star that represents the hypocentral location on cross section B-B.’ The hypocentral depth suggests this M 5.5 earthquake is in the downgoing Australia plate slab.
    • In the lower left corner is a diagram showing the subducting Australia plate at the Java Trench (Yves Descatoire).
    • To the right of the Java Trench figure presents a detailed view of the faulting along the eastern Java and western Timor trenches (Hangesh and Whitney, 2016). They present evidence for oblique motion along the Timor trough. And present evidence for a backthrust on the northern side of Timor and the Indonesia islands east of Java.


    • Here is the map from Baldwin et al. (2012)

    • 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 phiolite; SB, Sepik Basin block; SDB, Sunda block; SBS, South Bismarck plate; SIB, Solomon Islands block; WP, Wandamen p ninsula; WDK, Woodlark microplate; YQ, Yeleme quarries.

    • Here is the tectonic map from Hangesh and Whitney (2016)

    • Illustration of major tectonic elements in triple junction geometry: tectonic features labeled per Figure 1; seismicity from ISC-GEM catalog [Storchak et al., 2013]; faults in Savu basin from Rigg and Hall [2011] and Harris et al. [2009]. Purple line is edge of Australian continental basement and fore arc [Rigg and Hall, 2011]. Abbreviations: AR = Ashmore Reef; SR = Scott Reef; RS = Rowley Shoals; TCZ = Timor Collision Zone; ST = Savu thrust; SB = Savu Basin; TT = Timor thrust; WT =Wetar thrust; WASZ = Western Australia Shear Zone. Open arrows indicate relative direction of motion; solid arrows direction of vergence.

    • Here is the Audley (2011) cross section showing how the backthrust relates to the subduction zone beneath Timor. I include their figure caption in blockquote below.

    • Cartoon cross section of Timor today, (cf. Richardson & Blundell 1996, their BIRPS figs 3b, 4b & 7; and their fig. 6 gravity model 2 after Woodside et al. 1989; and Snyder et al. 1996 their fig. 6a). Dimensions of the filled 40 km deep present-day Timor Tectonic Collision Zone are based on BIRPS seismic, earthquake seismicity and gravity data all re-interpreted here from Richardson & Blundell (1996) and from Snyder et al. (1996). NB. The Bobonaro Melange, its broken formation and other facies are not indicated, but they are included with the Gondwana mega-sequence. Note defunct Banda Trench, now the Timor TCZ, filled with Australian continental crust and Asian nappes that occupy all space between Wetar Suture and the 2–3 km deep deformation front north of the axis of the Timor Trough. Note the much younger decollement D5 used exactly the same part of the Jurassic lithology of the Gondwana mega-sequence in the older D1 decollement that produced what appears to be much stronger deformation.

    • Here is a figure showing the regional geodetic motions (Bock et al., 2003). I include their figure caption below as a blockquote.

    • Topographic and tectonic map of the Indonesian archipelago and surrounding region. Labeled, shaded arrows show motion (NUVEL-1A model) of the first-named tectonic plate relative to the second. Solid arrows are velocity vectors derived from GPS surveys from 1991 through 2001, in ITRF2000. For clarity, only a few of the vectors for Sumatra are included. The detailed velocity field for Sumatra is shown in Figure 5. Velocity vector ellipses indicate 2-D 95% confidence levels based on the formal (white noise only) uncertainty estimates. NGT, New Guinea Trench; NST, North Sulawesi Trench; SF, Sumatran Fault; TAF, Tarera-Aiduna Fault. Bathymetry [Smith and Sandwell, 1997] in this and all subsequent figures contoured at 2 km intervals.

    References:

    • Audley-Charles, M.G., 1986. Rates of Neogene and Quaternary tectonic movements in the Southern Banda Arc based on micropalaeontology in: Journal of fhe Geological Society, London, Vol. 143, 1986, pp. 161-175.
    • Audley-Charles, M.G., 2011. Tectonic post-collision processes in Timor, 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, 241–266.
    • Baldwin, S.L., Fitzgerald, P.G., and Webb, L.E., 2012. Tectonics of the New Guinea Region in Annu. Rev. Earth Planet. Sci., v. 41, p. 485-520.
    • Benz, H.M., Herman, Matthew, Tarr, A.C., Hayes, G.P., Furlong, K.P., Villaseñor, Antonio, Dart, R.L., and Rhea, Susan, 2011. Seismicity of the Earth 1900–2010 New Guinea and vicinity: U.S. Geological Survey Open-File Report 2010–1083-H, scale 1:8,000,000.
    • Hall, R., 2011. Australia-SE Asia collision: plate tectonics and crustal flow in Geological Society, London, Special Publications 2011; v. 355; p. 75-109 doi: 10.1144/SP355.5
    • Hangesh, J. and Whitney, B., 2014. Quaternary Reactivation of Australia’s Western Passive Margin: Inception of a New Plate Boundary? in: 5th International INQUA Meeting on Paleoseismology, Active Tectonics and Archeoseismology (PATA), 21-27 September 2014, Busan, Korea, 4 pp.
    • Hayes, G.P., Wald, D.J., and Johnson, R.L., 2012. Slab1.0: A three-dimensional model of global subduction zone geometries in, J. Geophys. Res., 117, B01302, doi:10.1029/2011JB008524
    • Hayes, G.P., Smoczyk, G.M., Benz, H.M., Villaseñor, Antonio, and Furlong, K.P., 2015. Seismicity of the Earth 1900–2013, Seismotectonics of South America (Nazca Plate Region): U.S. Geological Survey Open-File Report 2015–1031–E, 1 sheet, scale 1:14,000,000, http://dx.doi.org/10.3133/ofr20151031E.
    • Okal, E. A., & Reymond, D., 2003. The mechanism of great Banda Sea earthquake of 1 February 1938: applying the method of preliminary determination of focal mechanism to a historical event in EPSL, v. 216, p. 1-15.
    • Zahirovic, S., Seton, M., and Müller, R.D., 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: Alaska

    We had an earthquake a few days ago along the Cook Strait west of Anchorage, Alaska. This earthquake happened nearby a couple earthquakes from the past 2 years that have similar senses of motion along faults that seem to be oriented the same. Here is my report for the 2016 M 7.1 earthquake.

      The USGS websites for the three large earthquakes with moment tensors plotted on the poster are here

    • 2015.07.29 M 6.3
    • 2016.01.24 M 7.1
    • 2017.03.02 M 5.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 also include seismicity from 2015-2017 for earthquakes with magnitudes M ≥ 4.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. The moment tensors for all three of these earthquakes are very similar.
    • 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 plots this close to the location of the fault as mapped by Hayes et al. (2012). All three earthquakes listed above plot clearly within the downgoing Pacific plate slab.

      I include some inset figures in the poster.

    • In the upper right corner 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 upper left corner I include an inset map from the USGS Seismicity History poster for this region (Benz et al., 2010). There is one seismicity cross section with its locations plotted on the map. The USGS plot these hypocenters along this cross section and I include that below (with the legend). I placed orange circles on the map and cross section showing the general location of the M 7.1 and M 5.5 earthquakes. The M 6.3 earthquake would plot almost in the same location as the M 7.1.
    • To the right of the Benz et al. (2010) figures is a map showing seismicity plotted as dots colored vs. depth. This map is from the Alaska Earthquake Center as presented by IRIS.
    • In the lower right corner are the MMI intensity maps for the two earthquakes listed above: 2015 M 6.3, 2016 M 7.1, and 2017 M 5.5. These figures were created by the USGS, but were made at different sizes, so they don’t match perfectly (don’t ask me why they keep changing the sizes of their figures and maps, I don’t know the answer). Below each map are plotted the reports from the Did You Feel It? USGS website for each earthquake. These reports are plotted as green dots with intensity on the vertical axes and distance on the horizontal axes. There are comparisons with Ground Motion Prediction Equation (attenuation relations) results (the orange model uses empirical data from central and eastern US earthquakes; the green model uses empirical data from earthquakes in California). Neither model seems appropriate given the DYFI results, though the California model works slightly better for the M 5.5 earthquake.



    Here is the same poster with only seismicity from the past 30 days plotted.

    • Dr. Peter Haeussler produced a cross section for this region, as prepared for the 2016 M 7.1 earthquake. Below is his description of this figure.

    • I made up a quick diagram (thanks Alaska Earthquake Center tools) showing the tectonic setting of the earthquake. This was a “Benioff zone” event, which means that the earthquake is related to bending of the subducting Pacific Plate as it slides into the mantle.

    • Here is a map for the earthquakes of magnitude greater than or equal to M 7.0 between 1900 and 2016. 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 a cross section showing the differences of vertical deformation between the coseismic (during the earthquake) and interseismic (between earthquakes).

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

    References: