I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.
There was a magnitude M 7.6 earthquake in Mexico on 1 September 2022.
https://earthquake.usgs.gov/earthquakes/eventpage/us7000i9bw/executive
I am catching up on some Earthquake Reports that I did not yet post since my website was being migrated to a more secure and reliable server (and more expensive).
The tectonics of coastal southwestern Mexico is dominated by the convergent plate boundary between the Cocos plate (to the southwest) and the North America plate (to the northeast). Here, the Cocos plate subducts below (goes underneath) the North America plate.
The fault between these plates is called a megathrust subduction zone fault and the plate boundary forms the Middle America trench.
This M 7.6 earthquake mechanism (the “moment tensor”) shows that this event was a compressional earthquake (reverse or thrust).
Based on it’s location, the event probably happened along the megathrust fault.
This earthquake even generated a tsunami recorded on tide gages in the region!
Below is my interpretive poster for this earthquake
- I plot the seismicity from the past month, with diameter representing magnitude (see legend). I include earthquake epicenters from 1921-2021 with magnitudes M ≥ 3.0 in one version.
- I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
- A review of the basic base map variations and data that I use for the interpretive posters can be found on the Earthquake Reports page. I have improved these posters over time and some of this background information applies to the older posters.
- Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.
- In the upper left corner is a map that shows the plates, their boundaries, and a century of seismicity.
- In the upper right are two maps that show models of how there may have been landslides or liquefaction because of the earthquake shaking and impacts. Read more about landslides and liquefaction here. I include the USGS epicenter as a red circle. However, these ground failure models are based on the USGS epicenter/location.
- On the center right is a map that shows the historic subduction zone earthquake history for the subduction zone offshore of Mexico (National Seismological Service of Mexico).
- To the left of those two maps is a low angle oblique view of the tectonic plates and how they are oriented relative to each other (Manea et al., 2013).
- In the lower right corner is a map that shows the ground shaking from the earthquake, with color representing intensity using the Modified Mercalli Intensity (MMI) scale. The closer to the earthquake, the stronger the ground shaking. The colors on the map represent the USGS model of ground shaking. The colored circles represent reports from people who posted information on the USGS Did You Feel It? part of the website for this earthquake. There are things that affect the strength of ground shaking other than distance, which is why the reported intensities are different from the modeled intensities.
- On the left, below the tectonic setting map is a plot that shows how the shaking intensity models and reports relate to each other. The horizontal axis is distance from the earthquake and the vertical axis is shaking intensity (using the MMI scale, just like in the map to the right: these are the same datasets).
- In the upper left-center is a figure that shows the USGS earthquake slip model. This shows how much the fault slipped in different areas (based on their modeling, not observation). The model shows that there were places that may have slipped over 1.25 meters (~4 feet).
- In the lower left is a series of plots from the tide gages in the region. The location of these gages are shown on the main map. The tsunami wave height (vertical distance between the peak of the wave and the trough of the wave) ranged from 0.6m to 1.7m.
I include some inset figures.
- I could not help myself. I am so excited to have this website back up and running, like a fully operational space station, that I include below some additional figures that help us understand the tectonic setting.
- Here is the low angle oblique view of this tectonic region (Manea et al., 2013).
- Here are the plots from the tide gages operated in the region.
- These tide gages are organized north to south, top to bottom.
- The dark blue line is the tidal forecast. The medium blue line is the tide gage record. The light blue line is the tide gage record minus the tidal forecast (basically the tsunami plus any other influence (like atmospheric pressure influences, or storm surge, etc.).
- The locations for these gages are labeled on the interpretive poster above.
- The earthquake origin time is labeled in orange.
- Time is presented in UTC.
- Here is tectonic map from Franco et al. (2012).
- These figures are from the USGS publication (Benz et al., 2011) that presents an educational poster about the historic seismicity and seismic hazard along the Middle America Trench.
- First is a map showing earthquake depth as color (green depth > red). Seismicity cross section B-B’ is shown on the map. Today’s M=6.6 quake is nearest this section.
- Here are some figures from Manea et al. (2013). First are the map and low angle oblique view of the Cocos plate.
- Here is a map showing the spreading ridge features, along with the plate boundary faults (Mann, 2007). This is similar to the inset map in the interpretive poster.
- Here is the McCann et al. (1979) summary figure, showing the earthquake history of the region.
- This is a more updated figure from Franco et al. (2005) showing the seismic gap.
- Here is a map from Franco et al. (2015) that shows the rupture patches for historic earthquakes in this region.
Supportive Figures
Development of the Tepic–Zacoalco (TZ), Colima, and Chapala rifts. The TZ rift is formed by the Rivera slab rollback, enhanced by the toroidal flow around the slab edges. The Colima rift is probably related with the oblique convergence between Rivera and NAM plates at ~5 Ma.
Tectonic setting of the Caribbean Plate. Grey rectangle shows study area of Fig. 2. Faults are mostly from Feuillet et al. (2002). PMF, Polochic–Motagua faults; EF, Enriquillo Fault; TD, Trinidad Fault; GB, Guatemala Basin. Topography and bathymetry are from Shuttle Radar Topography Mission (Farr&Kobrick 2000) and Smith & Sandwell (1997), respectively. Plate velocities relative to Caribbean Plate are from Nuvel1 (DeMets et al. 1990) for Cocos Plate, DeMets et al. (2000) for North America Plate and Weber et al. (2001) for South America Plate.
A. Geodynamic and tectonic setting alongMiddle America Subduction Zone. JB: Jalisco Block; Ch. Rift—Chapala rift; Co. rift—Colima rift; EGG—El Gordo Graben; EPR: East Pacific Rise; MCVA: Modern Chiapanecan Volcanic Arc; PMFS: Polochic–Motagua Fault System; CR—Cocos Ridge. Themain Quaternary volcanic centers of the TransMexican Volcanic Belt (TMVB) and the Central American Volcanic Arc (CAVA) are shown as blue and red dots, respectively. B. 3-D view of the Pacific, Rivera and Cocos plates’ bathymetrywith geometry of the subducted slab and contours of the depth to theWadati–Benioff zone (every 20 km). Grey arrows are vectors of the present plate convergence along theMAT. The red layer beneath the subducting plate represents the sub-slab asthenosphere.
Marine magnetic anomalies and fracture zones that constrain tectonic reconstructions such as those shown in Figure 4 (ages of anomalies are keyed to colors as explained in the legend; all anomalies shown are from University of Texas Institute for Geophysics PLATES [2000] database): (1) Boxed area in solid blue line is area of anomaly and fracture zone picks by Leroy et al. (2000) and Rosencrantz (1994); (2) boxed area in dashed purple line shows anomalies and fracture zones of Barckhausen et al. (2001) for the Cocos plate; (3) boxed area in dashed green line shows anomalies and fracture zones from Wilson and Hey (1995); and (4) boxed area in red shows anomalies and fracture zones from Wilson (1996). Onland outcrops in green are either the obducted Cretaceous Caribbean large igneous province, including the Siuna belt, or obducted ophiolites unrelated to the large igneous province (Motagua ophiolites). The magnetic anomalies and fracture zones record the Cenozoic relative motions of all divergent plate pairs infl uencing the Central American subduction zone (Caribbean, Nazca, Cocos, North America, and South America). When incorporated into a plate model, these anomalies and fracture zones provide important constraints on the age and thickness of subducted crust, incidence angle of subduction, and rate of subduction for the Central American region. MCSC—Mid-Cayman Spreading Center.
Rupture zones (ellipses) and epicenters (triangles and circles) of large shallow earthquakes (after KELLEHER et al., 1973) and bathymetry (CHASE et al., 1970) along the Middle America arc. Note that six gaps which have earthquake histories have not ruptured for 40 years or more. In contrast, the gap near the intersection of the Tehuantepec ridge has no known history of large shocks. Contours are in fathoms.
The study area encompasses Guerrero and Oaxaca states of Mexico. Shaded ellipse-like areas annotated with the years are rupture areas of the most recent major thrust earthquakes (M≥6.5) in the Mexican subduction zone. Triangles show locations of permanent GPS stations. Small hexagons indicate campaign GPS sites. Arrows are the Cocos-North America convergence vectors from NUVEL-1A model (DeMets et al., 1994). Double head arrow shows the extent of the Guerrero seismic gap. Solid and dashed curves annotated with negative numbers show the depth in km down to the surface of subducting Cocos plate (modified from Pardo and Su´arez, 1995, using the plate interface configuration model for the Central Oaxaca from this study, the model for Guerrero from Kostoglodov et al. (1996), and the last seismological estimates in Chiapas by Bravo et al. (2004). MAT, Middle America trench.
- There are many different ways in which a landslide can be triggered. The first order relations behind slope failure (landslides) is that the “resisting” forces that are preventing slope failure (e.g. the strength of the bedrock or soil) are overcome by the “driving” forces that are pushing this land downwards (e.g. gravity). The ratio of resisting forces to driving forces is called the Factor of Safety (FOS). We can write this ratio like this:
FOS = Resisting Force / Driving Force
- When FOS > 1, the slope is stable and when FOS < 1, the slope fails and we get a landslide. The illustration below shows these relations. Note how the slope angle α can take part in this ratio (the steeper the slope, the greater impact of the mass of the slope can contribute to driving forces). The real world is more complicated than the simplified illustration below.
- Landslide ground shaking can change the Factor of Safety in several ways that might increase the driving force or decrease the resisting force. Keefer (1984) studied a global data set of earthquake triggered landslides and found that larger earthquakes trigger larger and more numerous landslides across a larger area than do smaller earthquakes. Earthquakes can cause landslides because the seismic waves can cause the driving force to increase (the earthquake motions can “push” the land downwards), leading to a landslide. In addition, ground shaking can change the strength of these earth materials (a form of resisting force) with a process called liquefaction.
- Sediment or soil strength is based upon the ability for sediment particles to push against each other without moving. This is a combination of friction and the forces exerted between these particles. This is loosely what we call the “angle of internal friction.” Liquefaction is a process by which pore pressure increases cause water to push out against the sediment particles so that they are no longer touching.
- An analogy that some may be familiar with relates to a visit to the beach. When one is walking on the wet sand near the shoreline, the sand may hold the weight of our body generally pretty well. However, if we stop and vibrate our feet back and forth, this causes pore pressure to increase and we sink into the sand as the sand liquefies. Or, at least our feet sink into the sand.
- Below is a diagram showing how an increase in pore pressure can push against the sediment particles so that they are not touching any more. This allows the particles to move around and this is why our feet sink in the sand in the analogy above. This is also what changes the strength of earth materials such that a landslide can be triggered.
- Below is a diagram based upon a publication designed to educate the public about landslides and the processes that trigger them (USGS, 2004). Additional background information about landslide types can be found in Highland et al. (2008). There was a variety of landslide types that can be observed surrounding the earthquake region. So, this illustration can help people when they observing the landscape response to the earthquake whether they are using aerial imagery, photos in newspaper or website articles, or videos on social media. Will you be able to locate a landslide scarp or the toe of a landslide? This figure shows a rotational landslide, one where the land rotates along a curvilinear failure surface.
- Here is an excellent educational video from IRIS and a variety of organizations. The video helps us learn about how earthquake intensity gets smaller with distance from an earthquake. The concept of liquefaction is reviewed and we learn how different types of bedrock and underlying earth materials can affect the severity of ground shaking in a given location. The intensity map above is based on a model that relates intensity with distance to the earthquake, but does not incorporate changes in material properties as the video below mentions is an important factor that can increase intensity in places.
- If we look at the map at the top of this report, we might imagine that because the areas close to the fault shake more strongly, there may be more landslides in those areas. This is probably true at first order, but the variation in material properties and water content also control where landslides might occur.
- There are landslide slope stability and liquefaction susceptibility models based on empirical data from past earthquakes. The USGS has recently incorporated these types of analyses into their earthquake event pages. More about these USGS models can be found on this page.
Earthquake Triggered Landslides
Social Media:
#EarthquakeReport for the M 7.6 (likely) subduction zone #Earthquake in #Mexico on 19 Sept 2022
catching up on reports that happened after my website went down
generated 0.6-1.7m wave height #Tsunami
probably triggered landslides/induced liquefactionhttps://t.co/9zpN2ZlAhw pic.twitter.com/tMDb4mSwCw— Jason "Jay" R. Patton (@patton_cascadia) November 9, 2022
- 2022.09.19 M 7.6 Mexico
- 2021.09.08 M 7.0 Mexico
- 2021.07.21 M 6.7 Panama
- 2020.06.23 M 7.4 Mexico
- 2019.02.01 M 6.6 Guatemala & Mexico
- 2018.02.16 M 7.2 Oaxaca, Mexico
- 2018.01.19 M 6.3 Gulf of California
- 2017.09.19 M 7.1 Puebla, Mexico
- 2017.09.19 M 7.1 Puebla, Mexico Update #1
- 2017.09.08 M 8.1 Chiapas, Mexico
- 2017.09.08 M 8.1 Chiapas, Mexico Update #1
- 2017.09.23 M 8.1 Chiapas, Mexico Update #2
- 2017.06.22 M 6.8 Guatemala
- 2017.06.14 M 6.9 Guatemala
- 2017.05.12 M 6.2 El Salvador
- 2017.03.29 M 5.7 Gulf of California
- 2016.11.24 M 7.0 El Salvador
- 2016.04.29 M 6.6 East Pacific Rise / MAT
- 2016.01.21 M 6.6 Mexico
- 2015.09.13 M 6.6 Gulf California
- 2015.09.13 M 6.6 Gulf California Update #1
- 2015.01.07 M 6.6 Panama
- 2015.01.31 M 5.5 Panama Update #1
- 2014.05.13 M 6.5 Panama
- 2014.05.13 M 6.5 Panama Update #1
- 2014.10.14 M 7.3 El Salvador
- 2013.10.20 M 6.4 Gulf California
Mexico | Central America
General Overview
Earthquake Reports
- 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
- Franco, A., C. Lasserre H. Lyon-Caen V. Kostoglodov E. Molina M. Guzman-Speziale D. Monterosso V. Robles C. Figueroa W. Amaya E. Barrier L. Chiquin S. Moran O. Flores J. Romero J. A. Santiago M. Manea V. C. Manea, 2012. Fault kinematics in northern Central America and coupling along the subduction interface of the Cocos Plate, from GPS data in Chiapas (Mexico), Guatemala and El Salvador in Geophysical Journal International., v. 189, no. 3, p. 1223-1236. DOI: https://doi.org/10.1111/j.1365-246X.2012.05390.x
- Franco, S.I., Kostoglodov, V., Larson, K.M., Manea, V.C>, Manea, M., and Santiago, J.A., 2005. Propagation of the 2001–2002 silent earthquake and interplate coupling in the Oaxaca subduction zone, Mexico in Earth Planets Space, v. 57., p. 973-985.
- Garcia-Casco, A., Projenza, J.A., Iturralde-Vinent, M.A., 2011. Subduction Zones of the Caribbean: the sedimentary, magmatic, metamorphic and ore-deposit records UNESCO/iugs igcp Project 546 Subduction Zones of the Caribbean in Geologica Acta, v. 9, no., 3-4, p. 217-224
- Benz, H.M., Dart, R.L., Villaseñor, Antonio, Hayes, G.P., Tarr, A.C., Furlong, K.P., and Rhea, Susan, 2011 a. Seismicity of the Earth 1900–2010 Mexico and vicinity: U.S. Geological Survey Open-File Report 2010–1083-F, scale 1:8,000,000.
- Franco, A., Lasserre, C., Lyon-Caen, H., Kostoglodov, V., Molina, E., Guzman-Speziale, M., Monterosso, D., Robles, V., Figueroa, C., Amaya, W., Barrier, E., Chiquin, L., Moran, S., Flores, O., Romero, J., Santiago, J.A., Manea, M., Manea, V.C., 2012. Fault kinematics in northern Central America and coupling along the subduction interface of the Cocos Plate, from GPS data in Chiapas (Mexico), Guatemala and El Salvador in Geophysical Journal International., v. 189, no. 3, p. 1223-1236 https://doi.org/10.1111/j.1365-246X.2012.05390.x
- Manea, V.C., et al., 2013. A geodynamical perspective on the subduction of Cocos and Rivera plates beneath Mexico and Central America in Tectonophysics, http://dx.doi.org/10.1016/j.tecto.2012.12.039
- Manea, M., and Manea, V.C., 2014. On the origin of El Chichón volcano and subduction of Tehuantepec Ridge: A geodynamical perspective in JGVR, v. 175, p. 459-471.
- Mann, P., 2007. Overview of the tectonic history of northern Central America, in Mann, P., ed., Geologic and tectonic development of the Caribbean plate boundary in northern Central America: Geological Society of America Special Paper 428, p. 1–19, doi: 10.1130/2007.2428(01). For
- McCann, W.R., Nishenko S.P., Sykes, L.R., and Krause, J., 1979. Seismic Gaps and Plate Tectonics” Seismic Potential for Major Boundaries in Pageoph, v. 117
References:
Basic & General References
Specific References
Return to the Earthquake Reports page.
- Sorted by Magnitude
- Sorted by Year
- Sorted by Day of the Year
- Sorted By Region
I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.
Well, it has been a busy couple of weeks. https://earthquake.usgs.gov/earthquakes/eventpage/us6000ah9t/executive The west coast coastline of southern Mexico, Central America, and South America is formed by a convergent plate boundary where oceanic tectonic plates dive eastwards beneath the continents. The fault formed at this plate boundary is called a subduction zone and the dynamics of subduction zones form deep sea trenches. I spend a few paragraphs discussing the different faults that form at different plate boundaries here. Offshore of southern Mexico the Middle America trench shows us the location of the subduction zone megathrust fault. This fault system has a long history of damaging earthquakes, including some events that affect areas hundreds of kilometers from the source earthquake (e.g. the 1985 magnitude M 8 Mexico City earthquake). In the past few years, evidences this megathrust is active continue to present themselves. There is a list of some earthquake reports at the bottom of this page. The earthquake generated seismic waves that travelled around the world, including some that caused strong shaking in Mexico City. Mexico City was built where the Aztec Civilization had once constructed a great city. This city was built next to a lake where the Aztec constructed floating gardens. Eventually, these gardens filled the lake and the lake filled with sediment (I am simplifying what happened over a long time). So, Mexico City is built in a sedimentary basin. Sedimentary basins can amplify shaking from seismic waves. These basins can also focus seismic waves and these waves can resonate within the basin, causing further amplification. This is why there was so much damage in Mexico City from the 1985 subduction zone earthquake. The same thing happened a couple years ago for a recent earthquake there. Well, when subduction zone earthquakes happen, the crust around the fault can flex like the elastic on one’s waist band. As the crust moves, if that crust is beneath the water, this crust motion moves the water causing a tsunami. There are a number of organizations that monitor the Earth for earthquakes that may cause tsunami. These organizations alert officials in regions where these tsunami may inundate so that residents and visitors to the coast can take action (e.g. head to high ground). These programs save lives. This M 7.4 earthquake generated a tsunami that was recorded along the coastline, but not at all tide gage stations. The Salina Cruz station has a great record of this tsunami and is located >80 km from the epicenter. The Acapulco station also recorded a tsunami, but those data were not uploaded to the IOC website (they are working this out now). It appeared that the Acapulco data were being streamed in real time, but I noticed that they were the same data as posted for the Salina Cruz station. Here I plot the water surface elevations observed at the Salina Cruz tide gage. I mark the earthquake event time and the tsunami arrival time, then calculate the tsunami travel time. The Wave Height of the tsunami is the vertical distance measured between the peak and the trough. These data show a Maximum Wave Height of 1.4 meters. The strong ground shaking from an earthquake can also cause landslides and liquefaction. I discuss these further down in this report and include maps in the poster.
Tectonic setting of the Caribbean Plate. Grey rectangle shows study area of Fig. 2. Faults are mostly from Feuillet et al. (2002). PMF, Polochic–Motagua faults; EF, Enriquillo Fault; TD, Trinidad Fault; GB, Guatemala Basin. Topography and bathymetry are from Shuttle Radar Topography Mission (Farr&Kobrick 2000) and Smith & Sandwell (1997), respectively. Plate velocities relative to Caribbean Plate are from Nuvel1 (DeMets et al. 1990) for Cocos Plate, DeMets et al. (2000) for North America Plate and Weber et al. (2001) for South America Plate.
A. Geodynamic and tectonic setting alongMiddle America Subduction Zone. JB: Jalisco Block; Ch. Rift—Chapala rift; Co. rift—Colima rift; EGG—El Gordo Graben; EPR: East Pacific Rise; MCVA: Modern Chiapanecan Volcanic Arc; PMFS: Polochic–Motagua Fault System; CR—Cocos Ridge. Themain Quaternary volcanic centers of the TransMexican Volcanic Belt (TMVB) and the Central American Volcanic Arc (CAVA) are shown as blue and red dots, respectively. B. 3-D view of the Pacific, Rivera and Cocos plates’ bathymetrywith geometry of the subducted slab and contours of the depth to theWadati–Benioff zone (every 20 km). Grey arrows are vectors of the present plate convergence along theMAT. The red layer beneath the subducting plate represents the sub-slab asthenosphere.
Kinematic model (mantle reference frame) of the subducting Cocos slab along the MAT in the vicinity of Cocos–Caribbe–North America triple junction since Early Miocene. The evolution of Caribbean–North America tectonic contact is based on the model of Witt et al. (2012). The blue strips represent markers on the Cocos plate. Note how trench roll forward is associated with steep slab in Central America, whereas trench roll back is associated with flat slab in Mexico.
Marine magnetic anomalies and fracture zones that constrain tectonic reconstructions such as those shown in Figure 4 (ages of anomalies are keyed to colors as explained in the legend; all anomalies shown are from University of Texas Institute for Geophysics PLATES [2000] database): (1) Boxed area in solid blue line is area of anomaly and fracture zone picks by Leroy et al. (2000) and Rosencrantz (1994); (2) boxed area in dashed purple line shows anomalies and fracture zones of Barckhausen et al. (2001) for the Cocos plate; (3) boxed area in dashed green line shows anomalies and fracture zones from Wilson and Hey (1995); and (4) boxed area in red shows anomalies and fracture zones from Wilson (1996). Onland outcrops in green are either the obducted Cretaceous Caribbean large igneous province, including the Siuna belt, or obducted ophiolites unrelated to the large igneous province (Motagua ophiolites). The magnetic anomalies and fracture zones record the Cenozoic relative motions of all divergent plate pairs infl uencing the Central American subduction zone (Caribbean, Nazca, Cocos, North America, and South America). When incorporated into a plate model, these anomalies and fracture zones provide important constraints on the age and thickness of subducted crust, incidence angle of subduction, and rate of subduction for the Central American region. MCSC—Mid-Cayman Spreading Center.
Present setting of Central America showing plates, Cocos crust produced at East Pacifi c Rise (EPR), and Cocos-Nazca spreading center (CNS), triple-junction trace (heavy dotted line), volcanoes (open triangles), Middle America Trench (MAT), and rates of relative plate motion (DeMets et al., 2000; DeMets, 2001). East Pacifi c Rise half spreading rates from Wilson (1996) and Barckhausen et al. (2001). Lines 1, 2, and 3 are locations of topographic and tomographic profi les in Figure 6.
(A) Tomographic slices of the P-wave velocity of the mantle at depths of 100, 300, and 500 km beneath Central America. (B) Upper panels show cross sections of topography and bathymetry. Lower panels: tomographic profi les showing Cocos slab detached below northern Central America, upper Cocos slab continuous with subducted plate at Middle America Trench (MAT), and slab gap between 200 and 500 km. Shading indicates anomalies in seismic wave speed as a ±0.8% deviation from average mantle velocities. Darker shading indicates colder, subducted slab material of Cocos plate. Circles are earthquake hypocenters. Grid sizes on profi les correspond to quantity of ray-path data within that cell of model; smaller boxes indicate regions of increased data density. CT—Cayman trough; SL—sea level (modifi ed from Rogers et al., 2002).
Rupture zones (ellipses) and epicenters (triangles and circles) of large shallow earthquakes (after KELLEHER et al., 1973) and bathymetry (CHASE et al., 1970) along the Middle America arc. Note that six gaps which have earthquake histories have not ruptured for 40 years or more. In contrast, the gap near the intersection of the Tehuantepec ridge has no known history of large shocks. Contours are in fathoms.
The study area encompasses Guerrero and Oaxaca states of Mexico. Shaded ellipse-like areas annotated with the years are rupture areas of the most recent major thrust earthquakes (M≥6.5) in the Mexican subduction zone. Triangles show locations of permanent GPS stations. Small hexagons indicate campaign GPS sites. Arrows are the Cocos-North America convergence vectors from NUVEL-1A model (DeMets et al., 1994). Double head arrow shows the extent of the Guerrero seismic gap. Solid and dashed curves annotated with negative numbers show the depth in km down to the surface of subducting Cocos plate (modified from Pardo and Su´arez, 1995, using the plate interface configuration model for the Central Oaxaca from this study, the model for Guerrero from Kostoglodov et al. (1996), and the last seismological estimates in Chiapas by Bravo et al. (2004). MAT, Middle America trench.
Tectonic framework of the Cocos plate convergent margin. Top- General view. Yellow arrows indicate direction and speed (in cm/yr) of plate convergence, calculated from the Euler poles given by DeMets et al. (2010) for CocoeNoam (first three arrows, from left to right), and CocoeCarb (last four arrows). Length of arrow is proportional to speed. Red arrow shows location of the 96 longitude. Box indicates location of lower panel. Bottom- Location of features and places mentioned in text. Triangles indicate volcanoes of the Central American Volcanic Arc (CAVA) with known Holocene eruption (Siebert and Simkin, 2002).
Seismicity along the convergent margin. Top: Map view. Blue circles are shallow (z < 60 km) hypocenters; orange, intermediate-depth (60 < z < 100 km); yellow, deep (z > 100 km). Next three panels: Earthquakes as a function of longitude and magnitude for shallow (blue dots), intermediate (orange), and deep (yellow) hypocenters. Numbers indicate number of events on each convergent margin, with average magnitude in parenthesis. Gray line in this and subsequent figures mark the 96 deg longitude.
Location of hypocentral cross-sections. Hypocentral depths are keyed as in previous figures.
Hypocentral cross-sections. Depths are color-coded as in previous figures. Dashed lines indicate the 60-km and 100-km depths. Tick marks are at 100-km intervals, as shown on the sections. There is no vertical exaggeration and Earth’s curvature is taken into account. Number of sections refers to location on Fig. 3.
Earthquake fault-plane solutions from CMT data. a. Shallow (z < 60 km), thrust-faulting mechanisms. b. Intermediate-depth (60 < z < 100 km) thrust-faulting events. c. Deep (z > 100 km), thrust-faulting earthquakes. d. to f. Normal-faulting events, in same layout as for thrust-faulting events.
FOS = Resisting Force / Driving Force #EarthquakeReport in Mexico for M 7.4 #Earthquake #Terremoto be safe everyonehttps://t.co/jTfZEjQ0H0 more on tectonic background herehttps://t.co/5jgHU1xxbx pic.twitter.com/H1moPjQtqR — Jason "Jay" R. Patton (@patton_cascadia) June 23, 2020 #EarthquakeReport for the M 7.4 #Terremoto #Earthquake #Tsunami in #Oaxaca #Mexico shaking felt in #MexicoCity poster previously posted is updated and explained in report here:https://t.co/jCM8vu0ccX pic.twitter.com/C7OD8UobyL — Jason "Jay" R. Patton (@patton_cascadia) June 25, 2020 Watch the waves from the M7.4 Mexico earthquake roll across seismic stations in North America! (THREAD) pic.twitter.com/7dpPr1Ueuy — IRIS Earthquake Sci (@IRIS_EPO) June 23, 2020 Tide gauge at Salina Cruz, Mexico#SantaMaríaZapotitlán for M7.4 earthquake. pic.twitter.com/RNVYbzxseM — Nick Graehl (@nickgraehl) June 23, 2020 #Terremoto en México: en estas impactantes imágenes, tomadas en Bosques de las Lomas (#CDMX), podemos observar como ambos lados de una grieta se mueven ligeramente durante el seísmo.pic.twitter.com/u50mCx1XFS — eSPAINews (@eSPAINews) June 23, 2020 Heres my first cut at estimating offsets for the event today in Mexico. I will add more sites as they become available. The mechanism is USGS W-phase. pic.twitter.com/PomjZF0JKL — Brendan Crowell (@bwcphd) June 23, 2020 #Earthquake_alert – Source: @Bersa76#CDMX Así se sintió el sismo en el nivel 56 de la Torre Mitikah en la alcaldía Benito Juárez, se trata de un inmueble de 267.3m de altura. pic.twitter.com/J6cpAZfPIU#Mexico #earthquake #alert #warning #mitigation #tsunami — Desianto F. Wibisono (@TDesiantoFW) June 23, 2020 VIDEO: A building swaying in Mexico City’s La Condesa neighborhood, after today’s 7.4 quake in southern Mexico, according to @TVVnoticias pic.twitter.com/BKJrilJZ8s — Manuel Bojorquez (@BojorquezCBS) June 23, 2020 "Nearly" automatic forward model based on the USGS focal mechanism for the M 7.4 #Mexico #earthquake (June 22, 15:29 UTC). And expected InSAR fringe pattern from C-band ascending orbit. Waiting, as usual, for real InSAR data. — Simone Atzori (@SimoneAtzori73) June 23, 2020 #DÚltimaHora #D7 Cerrada la Carretera Oaxaca-Istmo por Derrumbe en las Cercanías de la Comunidad del Camarón, evite Viajar más información en un momento. pic.twitter.com/B3G5b17YZv — Domingo 7 (@domingo7th) June 23, 2020 A 7.4 magnitude #earthquake has hit the southern coastline of #Mexico. #Sismo #Temblor #Oaxaca #CDMX #AlertaSismica pic.twitter.com/doiICbJHuy — DailyNewsf (@DailyNe25683877) June 23, 2020 #sismo2020 #sismo | En la Roma Norte, el movimiento causó daños en la azotea de una edificio, del cual comenzaron a caer chorros de agua pic.twitter.com/LqA4fVKhBy — RedTvo Television (@redtvo) June 23, 2020 El #sismo magnitud 7.1 con epicentro en Crucecita Huatulco, #Oaxaca. Ocasionó daños en el Centro Histórico de la Capital. pic.twitter.com/iD70elgxP1 — Diario El Fortín (@diarioelfortin) June 23, 2020 En la zona platanera de Teapa, Tabasco, así se percibió el sismo de 7.5 #CDMX #oaxaca pic.twitter.com/sgpd4lKvOs — Desde Peninsula (@DesdePeninsula) June 23, 2020 Playa Riscalillo, mar se aleja unos cien metros. Bahía principal de Huatulco, mar se aleja unos 30-40 metros#ReporteCiudadano#TenemosSismo #Oaxaca #Huatulco pic.twitter.com/rZNa5AWDFr — Periodistas Oaxaca (@PeriodistasOax) June 23, 2020 OXUM roughly 75 km away to the west had pretty sizable vertical motions. pic.twitter.com/YIltT3uBE4 — Brendan Crowell (@bwcphd) June 23, 2020 Así se sintió el sismo de 7.5 está mañana en Xochimilco pic.twitter.com/yK0CCKwGNt — Salvador García Soto (@SGarciaSoto) June 23, 2020 Earthquake in Mexico City- still occurring – not sure where epicenter is nor if there is any real damage pic.twitter.com/KDAou9YhPG — (@Andalalucha) June 23, 2020 First rough model of the tsunami from today's M7.5 earthquake in Mexico using the USGS finite fault as input. Clear edgewaves and shelf resonance in the Tehuantepec gulf. https://t.co/2FAV5m81ZX — Diego Melgar (@geosmx) June 23, 2020 Today's Oaxaca tsunami in the Pacific. So far, only on 43413… pic.twitter.com/hDqG2ip3jI — Amir Salaree (@amirsalaree) June 23, 2020 El 23 de junio del 2020, un terremoto M7.4 cerca de Oaxaca, México, ocurrió como resultado de una falla inversa en o cerca del límite de la placa que está entre las placas de Cocos y América del Norte. (Enlaces y detalles: https://t.co/qvcse5eYQe) pic.twitter.com/53IxlOr0VG — USGS (@USGS) June 23, 2020 A section from today's M7.4 earthquake near Santa María Zapotitlán at 2020-06-23 15:29:05 UTC recorded on the worldwide @raspishake network and processed with a wider frequency range than I normally use. See: https://t.co/44McZsCCpE. Uses @obspy & @matplotlib software. pic.twitter.com/kItaJfsJ4W — Mark Vanstone (@wmvanstone) June 23, 2020 GSN and other global surface & body wave record sections for the M7.4 Mexico earthquake https://t.co/FvU7Jlak3R pic.twitter.com/YdvRb3iInc — IRIS Earthquake Sci (@IRIS_EPO) June 23, 2020 How big was that magnitude 7.4 earthquake? According to my calculations (see below), it's as if about 4 trillion college students jumped about 1 ft. And, that's about 500 times the population of planet Earth! (Peer review of my calculations is welcome…) pic.twitter.com/TlrmXhQSLj — Alan Kafka (@Weston_Quakes) June 24, 2020 #Mexico M7.4 #earthquake — Vincenzo De Novellis (@VDN75) June 24, 2020 #ERCC #DailyMap: 2020-06-24 ⦙ Mexico | 7.4M Earthquake of 23 June ▸https://t.co/mdRq0Me6LL pic.twitter.com/z7vgXVrfRG — Copernicus EMS (@CopernicusEMS) June 24, 2020 I find it mightily impressive that the low-cost geophone sensors inside @raspishake seismometers were able to capture the long-period surface wave signals, as well as higher-freq body waves, from yesterday's M7.4 Mexico earthquake. — Stephen Hicks (@seismo_steve) June 24, 2020
This morning (my time) there was a moderately deep earthquake along the coast of southern Mexico and northern Guatemala. Here is my Temblor article about this M=6.6 earthquake and how it might relate to the 2017 M=8.2 quake. I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include earthquake epicenters from 1919-2019 with magnitudes M ≥ 6.5 in one version. There are also some interesting relations between different historic earthquakes.
A. Geodynamic and tectonic setting alongMiddle America Subduction Zone. JB: Jalisco Block; Ch. Rift—Chapala rift; Co. rift—Colima rift; EGG—El Gordo Graben; EPR: East Pacific Rise; MCVA: Modern Chiapanecan Volcanic Arc; PMFS: Polochic–Motagua Fault System; CR—Cocos Ridge. Themain Quaternary volcanic centers of the TransMexican Volcanic Belt (TMVB) and the Central American Volcanic Arc (CAVA) are shown as blue and red dots, respectively. B. 3-D view of the Pacific, Rivera and Cocos plates’ bathymetrywith geometry of the subducted slab and contours of the depth to theWadati–Benioff zone (every 20 km). Grey arrows are vectors of the present plate convergence along theMAT. The red layer beneath the subducting plate represents the sub-slab asthenosphere.
Kinematic model (mantle reference frame) of the subducting Cocos slab along the MAT in the vicinity of Cocos–Caribbe–North America triple junction since Early Miocene. The evolution of Caribbean–North America tectonic contact is based on the model of Witt et al. (2012). The blue strips represent markers on the Cocos plate. Note how trench roll forward is associated with steep slab in Central America, whereas trench roll back is associated with flat slab in Mexico.
Present setting of Central America showing plates, Cocos crust produced at East Pacifi c Rise (EPR), and Cocos-Nazca spreading center (CNS), triple-junction trace (heavy dotted line), volcanoes (open triangles), Middle America Trench (MAT), and rates of relative plate motion (DeMets et al., 2000; DeMets, 2001). East Pacifi c Rise half spreading rates from Wilson (1996) and Barckhausen et al. (2001). Lines 1, 2, and 3 are locations of topographic and tomographic profi les in Figure 6.
(A) Tomographic slices of the P-wave velocity of the mantle at depths of 100, 300, and 500 km beneath Central America. (B) Upper panels show cross sections of topography and bathymetry. Lower panels: tomographic profi les showing Cocos slab detached below northern Central America, upper Cocos slab continuous with subducted plate at Middle America Trench (MAT), and slab gap between 200 and 500 km. Shading indicates anomalies in seismic wave speed as a ±0.8% deviation from average mantle velocities. Darker shading indicates colder, subducted slab material of Cocos plate. Circles are earthquake hypocenters. Grid sizes on profi les correspond to quantity of ray-path data within that cell of model; smaller boxes indicate regions of increased data density. CT—Cayman trough; SL—sea level (modifi ed from Rogers et al., 2002).
Proposed model of faults kinematics and coupling along the Cocos slab interface, revised from Lyon-Caen et al. (2006). Numbers are velocities relative to CA plate in mmyr−1. Focal mechanisms are for crustal earthquakes (depth ≤30 km) since 1976, from CMT Harvard catalogue.
The two beach balls show the stike-slip fault motions for the M6.4 (left) and M6.0 (right) earthquakes. Helena Buurman's primer on reading those symbols is here. pic.twitter.com/aWrrb8I9tj — AK Earthquake Center (@AKearthquake) August 15, 2018
Strike Slip: 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.)
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.
“When a fault slips during an #earthquake, there are changes in stress in the surrounding crust. These changes can either promote or inhibit the subsequent earthquake, depending on the orientation and type of fault on which the stress is imparted.” https://t.co/Qj4iOTQLwR — Dr Lucy Jones Center (@DLJCSS) February 3, 2019
Here I summarize Earth’s significant seismicity for 2018. I limit this summary to earthquakes with magnitude greater than or equal to M 6.5. I am sure that there is a possibility that your favorite earthquake is not included in this review. Happy New Year. One year of #earthquakes recorded by @INGVterremoti in Italy. About 2500 events with magnitude equal or larger than M2, about seven per day. Data source https://t.co/g1RvR2A989) #Italia #terremoto #Italy #earthquake pic.twitter.com/ft8GAsFjKA — iunio iervolino (@iuniervo) December 31, 2018 Earthquakes of 2018: a quick post summarising global seismic activity last year (i.e., the figures I showed you yesterday). https://t.co/ahdwpf1OFv pic.twitter.com/S438okD8QQ — Chris Rowan (@Allochthonous) January 1, 2019 Global #earthquakes by Magnitude (M5+) by year (2000-18), showing remarkable consistency from geologic forcing. Whereas patterns are understood, they do not permit short-term, local predictions; instead, be informed and be prepared. #geohazards @IRIS_EPO @USGS pic.twitter.com/BmtXhhUvWF — Ben van der Pluijm 🌎 (@vdpluijm) January 2, 2019 The pattern of shallow earthquakes (depth < 33 km) is typical, with much of the country susceptible to regular shallow seismicity, with lower rates in Northland/Auckland and southeast Otago. pic.twitter.com/3jip8Lyje9 — John Ristau 🇨🇦 🇳🇿 (@SinistralSeismo) January 3, 2019
Just a couple hours ago there was an earthquake along the Swan fault, which is the transform plate boundary between the North America and Caribbean plates. The Cayman trough (CT) is a region of oceanic crust, formed at the Mid-Cayman Rise (MCR) oceanic spreading center. To the west of the MCR the CT is bound by the left-lateral strike-slip Swan fault. To the east of the MCR, the CT is bound on the north by the Oriente fault. We had a damaging and (sadly) deadly earthquake in southern Peru in the last 24 hours. This is an earthquake, with magnitude M 7.1, that is associated with the subduction zone forming the Peru-Chile trench (PCT). The Nazca plate (NP) is subducting beneath the South America plate (SAP). There are lots of geologic structures on the Nazca plate that tend to affect how the subduction zone responds during earthquakes (e.g. segmentation). This earthquake appears to be located along a reactivated fracture zone in the GA. There have only been a couple earthquakes in this region in the past century, one an M 6.0 to the east (though this M 6.0 was a thrust earthquake). The Gulf of Alaska shear zone is even further to the east and has a more active historic fault history (a pair of earthquakes in 1987-1988). The magnetic anomalies (formed when the Earth’s magnetic polarity flips) reflect a ~north-south oriented spreading ridge (the anomalies are oriented north-south in the region of today’s earthquake). There is a right-lateral offset of these magnetic anomalies located near the M 7.9 epicenter. Interesting that this right-lateral strike-slip fault (?) is also located at the intersection of the Gulf of Alaska shear zone and the 1988 M 7.8 earthquake (probably just a coincidence?). However, the 1988 M 7.8 earthquake fault plane solution can be interpreted for both fault planes (it is probably on the GA shear zone, but I don’t think that we can really tell). As a reminder, if the M 7.9 earthquake fault is E-W oriented, it would be left-lateral. The offset magnetic anomalies show right-lateral offset across these fracture zones. This was perhaps the main reason why I thought that the main fault was not E-W, but N-S. After a day’s worth of aftershocks, the seismicity may reveal some north-south trends. But, as a drama student in 7th grade (1977), my drama teacher (Ms. Naichbor, rest in peace) asked our class to go stand up on stage. We all stood in a line and she mentioned that this is social behavior, that people tend to stand in lines (and to avoid doing this while on stage). Later, when in college, professors often commented about how people tend to seek linear trends in data (lines). I actually see 3-4 N-S trends and ~2 E-W trends in the seismicity data. There was just now an earthquake in Oaxaca, Mexico between the other large earthquakes from last 2017.09.08 (M 8.1) and 2017.09.08 (M 7.1). There has already been a M 5.8 aftershock.Here is the USGS website for today’s M 7.2 earthquake. This morning (local time in California) there was an earthquake in Papua New Guinea with, unfortunately, a high likelihood of having a good number of casualties. I was working on a project, so could not immediately begin work on this report. We had an M 6.8 earthquake near a transform micro-plate boundary fault system north of New Ireland, Papua New Guinea today. Here is the USGS website for this earthquake. The New Britain region is one of the more active regions in the world. See a list of earthquake reports for this region at the bottom of this page, above the reference list. Well, those earthquakes from earlier, one a foreshock to a later one, were foreshocks to an earthquake today! Here is my report from a couple days ago. The M 6.6 and M 6.3 straddle today’s earthquake and all have similar hypocentral depths. A couple days ago there was a deep focus earthquake in the downgoing Nazca plate deep beneath Bolivia. This earthquake has an hypocentral depth of 562 km (~350 miles). There has been a swarm of earthquakes on the southeastern part of the big island, with USGS volcanologists hypothesizing about magma movement and suggesting that an eruption may be imminent. Here is a great place to find official USGS updates on the volcanism in Hawaii (including maps). This version includes earthquakes M ≥ 3.5 (note the seismicity offshore to the south, this is where the youngest Hawaii volcano is). Below are a series of plots from tide gages installed at several sites in the Hawaii Island Chain. These data are all posted online here and here. Yesterday morning, as I was recovering from working on stage crew for the 34th Reggae on the River (fundraiser for the non profit, the Mateel Community Center), I noticed on social media that there was an M 6.9 earthquake in Lombok, Indonesia. This is sad because of the likelihood for casualties and economic damage in this region. Well, yesterday while I was installing the final window in a reconstruction project, there was an earthquake along the Aleutian Island Arc (a subduction zone) in the region of the Andreanof Islands. Here is the USGS website for the M 6.6 earthquake. This earthquake is close to the depth of the megathrust fault, but maybe not close enough. So, this may be on the subduction zone, but may also be on an upper plate fault (I interpret this due to the compressive earthquake fault mechanism). The earthquake has a hypocentral depth of 20 km and the slab model (see Hayes et al., 2013 below and in the poster) is at 40 km at this location. There is uncertainty in both the slab model and the hypocentral depth. We just had a Great Earthquake in the region of the Fiji Islands, in the central-western Pacific. Great Earthquakes are earthquakes with magnitudes M ≥ 8.0. This ongoing sequence began in late July with a Mw 6.4 earthquake. Followed less than 2 weeks later with a Mw 6.9 earthquake. We just had a M 7.3 earthquake in northern Venezuela. Sadly, this large earthquake has the potential to be quite damaging to people and their belongings (buildings, infrastructure). Well, this earthquake, while having a large magnitude, was quite deep. Because earthquake intensity decreases with distance from the earthquake source, the shaking intensity from this earthquake was so low that nobody submitted a single report to the USGS “Did You Feel It?” website for this earthquake. Following the largest typhoon to strike Japan in a very long time, there was an earthquake on the island of Hokkaido, Japan today. There is lots on social media, including some spectacular views of disastrous and deadly landslides triggered by this earthquake (earthquakes are the number 1 source for triggering of landslides). These landslides may have been precipitated (sorry for the pun) by the saturation of hillslopes from the typhoon. Based upon the USGS PAGER estimate, this earthquake has the potential to cause significant economic damages, but hopefully a small number of casualties. As far as I know, this does not incorporate potential losses from earthquake triggered landslides [yet]. Today, there was a large earthquake associated with the subduction zone that forms the Kermadec Trench. Well, around 3 AM my time (northeastern Pacific, northern CA) there was a sequence of earthquakes including a mainshock with a magnitude M = 7.5. This earthquake happened in a highly populated region of Indonesia. Here is a map that shows the updated USGS model of ground shaking. The USGS prepared an updated earthquake fault slip model that was additionally informed by post-earthquake analysis of ground deformation. The original fault model extended from north of the epicenter to the northernmost extent of Palu City. Soon after the earthquake, Dr. Sotiris Valkaniotis prepared a map that showed large horizontal offsets across the ruptured fault along the entire length of the western margin on Palu Valley. This horizontal offset had an estimated ~8 meters of relative displacement. InSAR analyses confirmed that the coseismic ground deformation extended through Palu Valley and into the mountains to the south of the valley. Synthetic Aperture Radar (SAR) is a remote sensing method that uses Radar to make observations of Earth. These observations include the position of the ground surface, along with other information about the material properties of the Earth’s surface. Landslides during and following the M=7.5 earthquake in central Sulawesi, Indonesia possibly caused the majority of casualties from this catastrophic natural disaster. Volunteers (citizen scientists) have used satellite aerial imagery collected after the earthquake to document the spatial extent and magnitude of damage caused by the earthquake, landslides, and tsunami.
Nowicki Jessee and others (2018) is the preferred model for earthquake-triggered landslide hazard. Our primary landslide model is the empirical model of Nowicki Jessee and others (2018). The model was developed by relating 23 inventories of landslides triggered by past earthquakes with different combinations of predictor variables using logistic regression. The output resolution is ~250 m. The model inputs are described below. More details about the model can be found in the original publication. We modify the published model by excluding areas with slopes <5° and changing the coefficient for the lithology layer "unconsolidated sediments" from -3.22 to -1.36, the coefficient for "mixed sedimentary rocks" to better reflect that this unit is expected to be weak (more negative coefficient indicates stronger rock).To exclude areas of insignificantly small probabilities in the computation of aggregate statistics for this model, we use a probability threshold of 0.002.
Zhu and others (2017) is the preferred model for liquefaction hazard. The model was developed by relating 27 inventories of liquefaction triggered by past earthquakes to globally-available geospatial proxies (summarized below) using logistic regression. We have implemented the global version of the model and have added additional modifications proposed by Baise and Rashidian (2017), including a peak ground acceleration (PGA) threshold of 0.1 g and linear interpolation of the input layers. We also exclude areas with slopes >5°. We linearly interpolate the original input layers of ~1 km resolution to 500 m resolution. The model inputs are described below. More details about the model can be found in the original publication.
In this region of the world, the Solomon Sea plate and the South Bismarck plate converge to form a subduction zone, where the Solomon Sea plate is the oceanic crust diving beneath the S.Bismarck plate. This region of the Pacific-North America plate boundary is at the northern end of the Cascadia subduction zone (CSZ). To the east, the Explorer and Juan de Fuca plates subduct beneath the North America plate to form the megathrust subduction zone fault capable of producing earthquakes in the magnitude M = 9 range. The last CSZ earthquake was in January of 1700, just almost 319 years ago. Before I looked more closely, I thought this sequence might be related to the Kefallonia fault. I prepared some earthquake reports for earthquakes here in the past, in 2015 and in 2016. There was a M = 6.8 earthquake along a transform fault connecting segments of the Mid Atlantic Ridge recently. Today’s earthquake occurred along the convergent plate boundary in southern Alaska. This subduction zone fault is famous for the 1964 March 27 M = 9.2 megathrust earthquake. I describe this earthquake in more detail here. There was a sequence of earthquakes along the subduction zone near New Caledonia and the Loyalty Islands. A large earthquake in the region of the Bering Kresla fracture zone, a strike-slip fault system that coincides with the westernmost portion of the Aleutian trench (which is a subduction zone further to the east). This magnitude M = 7.0 earthquake is related to the subduction zone that forms the Philippine trench (where the Philippine Sea plate subducts beneath the Sunda plate). Here is the USGS website for this earthquake.
The two beach balls show the stike-slip fault motions for the M6.4 (left) and M6.0 (right) earthquakes. Helena Buurman's primer on reading those symbols is here. pic.twitter.com/aWrrb8I9tj — AK Earthquake Center (@AKearthquake) August 15, 2018
Strike Slip: 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.)
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.
Today we had an earthquake with magnitude M 6.3 in the Gulf of California (GOC). The GOC is formed by transtension (extension along a strike-slip fault system) along the North-America-Pacific plate boundary. Transtension happens when the plate motion across a fault is not oriented parallel to the fault. This non-optimal relation (plate motion vs fault direction) can generally happen as a result of (1) a bending fault or (2) due to stepped offsets of the fault. The dextral (right-lateral) strike-slip faults here make “right-steps” and pull apart basins form in these locations. 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 ≥ 6.5.
Modified integrated transtensional shear model for the tectonic evolution of the Gulf of California. North America plate fixed. (A) Prior to 28 Ma, the spreading center between Pacific and Vancouver-Farallon tectonic plates approached the subduction zone between North America and Vancouver-Farallon plates. (B) By ca. 20 Ma, contact between the Pacific and North America plates created early dextral transform relative plate motion. The Basin and Range extensional province (light gray) accommodated moderate extension throughout western North America. (C) The Rivera triple junction migrated the full length of the Baja California Peninsula by ca. 12.5 Ma, lengthening the Pacific–North American transform plate boundary. The proto–Gulf of California period commenced (ca.12.5 Ma) with transtensional strain distributed across two distinct transtensional deformation belts, west and east of the stable Baja California microplate. (D) In late proto–Gulf of California time, shear deformation gradually localized within a narrow belt of focused en- echelon dextral shear zones embedded within the greater Mexican Basin and Range extensional province. These shear zones and intervening extensional regions both experienced high-magnitude strain. (E) By ca. 6 Ma, Pacific–North America plate boundary strain was localized and focused crustal thinning and subsidence in transtensional pull-apart basins that formed the Gulf of California. Faults shown represent primary structures active during Quaternary time. RP—Rivera plate, JDFP—Juan de Fuca plate, RTJ—Rivera triple junction, MTJ—Mendocino triple junction, SAF— San Andreas fault, CP—Colorado Plateau, SMO—Sierra Madre Occidental.
The model for the tectonic evolution of the Gulf of California (GOC). (A) The starting point of the evolution of the GOC began 14–12 Ma, where the Magdalena rise stalled off the west coast of Baja California; there was also a marked changed in the style of volcanism. Plate motion was split between the dying spreading ridge, subduction zone, and the new highly oblique extension in the proto–GOC. During this time, the dipping part of the subducted plate appears to have broken off, opening a slab window beneath the southern Baja peninsula. NAM—North American plate; PAC—Pacific plate. (B) Another major change in the system occurs near 8–7 Ma, where the volcanic style changes once again with many lavas of unusual composition deposited. Any minor component of spreading finally ceases and the Tosco-Abreojos fault forms within the borderlands west of Baja. Oblique extension continues in the GOC. (C) Seafloor spreading begins at the Alarcón Rise between 4 and 3 Ma. Small amounts of movement continue along the Tosco-Abreojos fault (TAF); even today the Baja peninsula is not fully transferred to the Pacific plate.
Map showing location of Alarcón transect and the major basins along the profile.
An overview of multichannel seismic transect data presented in this paper. Seismic data are post-stack time migrated. TWTT—two-way travel time.
East Cerralvo basin. The uninterpreted (top) and interpreted (bottom) seismic profi les are shown. TWTT—two-way traveltime. Basement reflections are highlighted in blue, and sedimentary sequence boundaries are separated by green lines, with faults shown in red. Basement has a reflective discontinuous appearance. Unit 1 (divided into 1a and 1b) is a synrift deposit, with a chaotic character; unit 2 (divided into 2a and 2b) exhibits rotation and divergence and appears to be synrift; unit 3 consists of postrift, layered deeper water marine sediments. Two surface-cutting normal faults at the southeast end appear younger than the main basin, although the exact relationship of faulting is unclear and they appear to be overprinted by current-controlled erosion.
(A) Simplified map of the Gulf of California region and Baja California peninsula showing the present plate boundary and some major tectonic features related to the plate-tectonic history since 12 Ma. The Gulf extensional province in gray is bounded by the Main Gulf Escarpment (bold dashed lines), which runs through the Loreto area and is shown in Figure 3. The Salton trough in southern California is merely the northern part of the Gulf extensional province. (B) Map of part of the southern Gulf of California and Baja California peninsula showing bathymetry (in meters), the transform–spreading-ridge plate boundary, and the location of subsequent figures with maps. The bathymetry is after a map in Ness and Lyle (1991) and the transform–spreading-ridge plate boundary is from Lonsdale (1989). The lines with double arrows are the three proposed rift segments modified here after Axen (1995); MS—Mulege´ segment, LS—Loreto segment, TS—Timbabichi
(A) Tectonic map of the southern Baja California microplate (BCM) and Gulf of California extensional province (GEP). The Magdalena fan is deposited on oceanic crust of the Farallon-derived Magdalena microplate located west of Baja California. Deep Sea Drilling Project Site 471 is shown as black dot on the Magdalena fan. Abbreviations: BCT—Baja California trench, BM—Bahia Magdalena, LC—Los Cabos block, T—Trinidad block, LP—La Paz, PV—Puerto Vallarta, SMSLF—Santa Margarita–San Lazaro fault, TAF—Tosco-Abreojos fault, TS—Todos Santos, V—Vizcaino peninsula. Geology is simplifi ed from Muehlberger (1996). Interpretation of marine magnetic anomalies, with numbers denoting the chron of positively magnetized stripes, is from Severinghaus and Atwater (1989) and Lonsdale (1991).
Map-view time slices showing the widely accepted model for the two-phase kinematic evolution of plate margin shearing around the Baja California microplate. (A) Configuration of active ridge segments (pink) west of Baja California just before they became largely abandoned ca. 12.3 Ma. (B) It is thought that plate motion from 12.3 to 6 Ma was kinematically partitioned into dextral strike slip (325 km) on faults west of Baja California and orthogonal rifting in the Gulf of California (90 km). This is known as the protogulf phase of rifting. (C) From 6 to 0 Ma faults west of Baja California are thought to have died and all plate motion was localized in the Gulf of California, which accommodated ~345 km of integrated transtensional shearing. Despite its wide acceptance, our data preclude this kinematic model. In all frames, the modern coastline is blue. Continental crust that accommodated post–12.3 Ma shearing is dark brown. Unfaulted microplates of continental crust are light tan. Farallon-derived microplates are light green. Middle Miocene trench-filling deposits like the Magdalena fan are colored dark green. Deep Sea Drilling Project Site 471 is the black dot on the southern Magdalena microplate. Yellow line (296 km) in the northern Gulf of California connects correlated terranes of Oskin and Stock (2003). Maps have Universal Transverse Mercator zone 12 projection with mainland Mexico fixed in present position.
Eastern Gulf of California contains abandoned rift basins, while active rifting occurs in the western Gulf (Lonsdale, 1989; Fenby and Gastil, 1991; Persaud et al., 2003; Aragón-Arreola et al., 2005; this study). Eastern Gulf constitutes abandoned rift margin (see inset). PA—Pacifi c plate; GC—Gulf of California; GEP—Gulf Extensional Province; B&R—Basin Gulf of California #Earthquake — Ehsan Kosari (@ehsan_seismo) January 19, 2018 M6.3 Earthquake, Max. Intensity V (light to moderate shaking felt) – 77km NNE of Loreto, Mexico (in the Gulf of California). 2018-01-19 09:17 AM local time. PAGER: Green. (Low likelihood of casualties and damage). @USGS Event Page: https://t.co/11DCL51sX8 pic.twitter.com/PlWUe2vq8m — SCEC (@SCEC) January 19, 2018 Mw=6.3, GULF OF CALIFORNIA (Depth: 6 km), 2018/01/19 16:17:44 UTC – Full details here: https://t.co/wIoiiCn8Sg pic.twitter.com/JS5Pmd33K4 — Earthquakes (@geoscope_ipgp) January 19, 2018 Gulf of California Tectonic Setting—Earthquakes & the Spreading Ridge https://t.co/jyEv9N5ly1 #earthquake #GulfofCalifornia pic.twitter.com/7sbuPtoIv1 — IRIS Earthquake Sci (@IRIS_EPO) January 19, 2018 Felt #earthquake (#sismo) M6.1 strikes 147 km SW of Ciudad Obregón (#Mexico) 9 min ago. https://t.co/w8mm58brIO pic.twitter.com/7oCiT5Jesw — EMSC (@LastQuake) January 19, 2018
Well, we had a really interesting earthquake today. There was a M 6.1 earthquake in the North America plate (NAP) to the north of the sequence offshore of Chiapas, with the M 8.1 mainshock. Here is the USGS website for the M 6.1 earthquake. There was also an M 5.8 earthquake that was a more typical aftershock (USGS website). Why is this earthquake interesting? It is outside the region of aftershocks from the M 8.1 earthquake and it is in the upper plate (the NAP). This is not altogether groundbreaking (pardon the pun) as there are many examples of earthquakes in one plate triggering earthquakes in other plates. For example, the recent sequence just to the south of the M 8.1 sequence (which may have led partly to the M 8.1 earthquake). This earthquake also triggered (sorry for the pun, another one) a debate about the difference between triggered earthquakes and aftershocks. This discussion is largely semantic and does not really matter from a natural hazards perspective. The rocks behave to physics, not how we classify them. So, we don’t need to get caught up in this lexicon (as long as we all have a general understanding of what is happening). In the classic sense, I interpret this M 6.1 (and the few nearby earthquakes) to be triggered, but they are in the region that may have an increased coulomb stress. I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I also include USGS epicenters from 1917-2017 for magnitudes M ≥ 8.0. I include fault plane solutions for the 1985 and 1995 earthquakes (along with the MMI contours for those earthquakes, see below for a discussion of MMI contours). Updated seismicity map, with today's two intraplate events: M5.8 within subducting slab, M6.1 at shallower depth within overriding plate pic.twitter.com/90iUgexfdi — Jascha Polet (@CPPGeophysics) September 23, 2017 Cross-section shows different tectonic environments for recent large quakes (excl. Sep 19): extension within slab versus in overriding plate pic.twitter.com/cZFfYUWLJF — Jascha Polet (@CPPGeophysics) September 23, 2017
M6.1 Mexico EQ this morning an upper plate aftershock of the fairly diverse aftershock sequence to the 09-08 M8.1 EQ. Slab xsecs here. pic.twitter.com/X7kCevaNYd — Gavin Hayes (@gph_seismo) September 23, 2017
Stress change at location of M6.1 Mexico EQ this morning as a result of prior 09-08 M8.1 EQ. Stress change on both planes about the same. pic.twitter.com/4GWNURpkex — Gavin Hayes (@gph_seismo) September 23, 2017
Well, the responses of people who are in the midst of a deadly disaster have been inspiring, bringing tears to my eyes often. Watching people searching and helping find survivors. This deadly earthquake brings pause to all who are paying attention. May we learn from this disaster with the hopes that others will suffer less from these lessons. I plot the USGS seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I also include USGS epicenters from 1917-2017 for magnitudes M ≥ 7.0. I include the USGS fault plane solution for the 1985 earthquake. I also include the USGS moment tensor for the 2017.09.08 M 8.1 earthquake. M 7.1 Puebla, Mexico earthquake at 11:14 AM PDT as recorded on the HSU Baby Benioff. The video is showing the surface waves arriving at campus, preceded by the P-waves at the beginning and S-waves immediately prior to the large amplitude waves. Our thoughts are with people in Mexico.
Topographic setting of Mexico City (MC) and the Valley of Mexico. Color scale corresponds to the basin thickness (i.e., the basin contact with the Oligocene volcanics of the Transmexican Volcanic Belt, TMVB). Stars show the epicenters for the vertical body forces applied at the free surface (green) and the magnitude 3.4 earthquake of December 1, 2014 (red). This figure has been created using the Generic Mapping Tools (GMT) Version 5.3.0, http://gmt.soest.hawaii.edu.
Snapshots of the Green’s function for the vertical body force S6 (see Fig. 1) described by the inset time history with flat spectrum up to 1 Hz. Notice the topographic scattering, the generation and propagation of wave trains at different speeds within the basin, and their multiple diffractions. This figure has been created using the Matlab software Version R2016a, http://www.mathworks.com/.
(a,c) Comparison of average eigenfunctions for the 8 sources with standard deviation bars for both elastic (blue solid) and viscoelastic (red solid) simulations at two representative sites, P1 and P2, and different frequencies. Dashed lines show theoretical eigenfunctions for the vertical component of Rayleigh waves in the model of Figure A1a (Table A1) for the fundamental mode (blue) and the first (red) and second (green) overtones. Normalized peak vertical displacements observed in different boreholes (green dots in Fig. 1) are shown with black circles and error bars (after Shapiro et al., 2001). (b) Fourier spectral amplifications (geometric mean of both horizontal components) at 0.5 Hz with respect to the CUIG site (Fig. 1) averaged for the 8 sources. The black contour corresponds to the 2 s dominant-period. (d) Duration of the strong shaking phase for f < 1 Hz averaged for the 8 sources.
Quickly made (disclaimer!) map of lake-bed of Lake Texcoco with crowd-sourced damaged and destroyed buildings map for Mexico City pic.twitter.com/MKCqeUqqJF — Jascha Polet (@CPPGeophysics) September 20, 2017
Tectonic framework of the Cocos plate convergent margin. Top- General view. Yellow arrows indicate direction and speed (in cm/yr) of plate convergence, calculated from the Euler poles given by DeMets et al. (2010) for CocoeNoam (first three arrows, from left to right), and CocoeCarb (last four arrows). Length of arrow is proportional to speed. Red arrow shows location of the 96 longitude. Box indicates location of lower panel. Bottom- Location of features and places mentioned in text. Triangles indicate volcanoes of the Central American Volcanic Arc (CAVA) with known Holocene eruption (Siebert and Simkin, 2002).
Seismicity along the convergent margin. Top: Map view. Blue circles are shallow (z < 60 km) hypocenters; orange, intermediate-depth (60 < z < 100 km); yellow, deep (z > 100 km). Next three panels: Earthquakes as a function of longitude and magnitude for shallow (blue dots), intermediate (orange), and deep (yellow) hypocenters. Numbers indicate number of events on each convergent margin, with average magnitude in parenthesis. Gray line in this and subsequent figures mark the 96 deg longitude.
Location of hypocentral cross-sections. Hypocentral depths are keyed as in previous figures.
Hypocentral cross-sections. Depths are color-coded as in previous figures. Dashed lines indicate the 60-km and 100-km depths. Tick marks are at 100-km intervals, as shown on the sections. There is no vertical exaggeration and Earth’s curvature is taken into account. Number of sections refers to location on Fig. 3.
Earthquake fault-plane solutions from CMT data. a. Shallow (z < 60 km), thrust-faulting mechanisms. b. Intermediate-depth (60 < z < 100 km) thrust-faulting events. c. Deep (z > 100 km), thrust-faulting earthquakes. d. to f. Normal-faulting events, in same layout as for thrust-faulting events.
The 2012 M = 8.6 mainshock and M = 8.2 aftershock fault ruptures and maps of strain duration tstrain at a threshold value of 0.1 microstrain. a, Inferred fault ruptures of the 11 April 2012 M = 8.6 east Indian Ocean earthquake and an M = 8.2 aftershock that occurred 2 h later. Superimposed are the first 20 d of M > 4.5 aftershocks of 0–100-km depth. These earthquakes probably ruptured a complex set of subparallel and conjugate faults with the indicated sense of motion (arrows). Parts of the rupture areas of the 2004 M = 9.2 and 2005 M = 8.7 Nias earthquakes on the Sunda megathrust are indicated. b, c, Global maps of tstrain (colour scale). Superimposed are the epicentres of M>5.5 events that occurred during the 6 d preceding the mainshock (2 epicentres) and following the mainshock (24 epicentres, 16 of which are remote, that is, .1,500km from the mainshock). Focal mechanisms of six post-mainshock events with near-vertical strike-slip mechanisms (plunge of neutral axis, >60 deg) are indicated with red beachballs. The 9:00:09 11 April 2012 M = 5.5 event (in the western Aleutian Islands) occurred 21 min 33 s after the mainshock between the direct P- and S-wave arrivals from the mainshock; all others are delayed by hours to days. The focal mechanism of the mainshock is plotted at its epicentre.
Well, after about 4 hours sleep, my business partner woke me up to talk about the fire alarms we were installing in a rental (#safetyfirst). Now that I have had some breakfast, I here provide some additional observations that people have made since I prepared my initial report. I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I also include USGS epicenters from 1917-2017 for magnitudes M ≥ 8.0. I include fault plane solutions for the 1985 earthquake (along with the MMI contours for that earthquake, see below for a discussion of MMI contours).
Tectonic setting of the Caribbean Plate. Grey rectangle shows study area of Fig. 2. Faults are mostly from Feuillet et al. (2002). PMF, Polochic–Motagua faults; EF, Enriquillo Fault; TD, Trinidad Fault; GB, Guatemala Basin. Topography and bathymetry are from Shuttle Radar Topography Mission (Farr&Kobrick 2000) and Smith & Sandwell (1997), respectively. Plate velocities relative to Caribbean Plate are from Nuvel1 (DeMets et al. 1990) for Cocos Plate, DeMets et al. (2000) for North America Plate and Weber et al. (2001) for South America Plate.
Rupture zones (ellipses) and epicenters (triangles and circles) of large shallow earthquakes (after KELLEHER et al., 1973) and bathymetry (CHASE et al., 1970) along the Middle America arc. Note that six gaps which have earthquake histories have not ruptured for 40 years or more. In contrast, the gap near the intersection of the Tehuantepec ridge has no known history of large shocks. Contours are in fathoms.
The study area encompasses Guerrero and Oaxaca states of Mexico. Shaded ellipse-like areas annotated with the years are rupture areas of the most recent major thrust earthquakes (M≥6.5) in the Mexican subduction zone. Triangles show locations of permanent GPS stations. Small hexagons indicate campaign GPS sites. Arrows are the Cocos-North America convergence vectors from NUVEL-1A model (DeMets et al., 1994). Double head arrow shows the extent of the Guerrero seismic gap. Solid and dashed curves annotated with negative numbers show the depth in km down to the surface of subducting Cocos plate (modified from Pardo and Su´arez, 1995, using the plate interface configuration model for the Central Oaxaca from this study, the model for Guerrero from Kostoglodov et al. (1996), and the last seismological estimates in Chiapas by Bravo et al. (2004). MAT, Middle America trench.
While I was spending time with my friend Steve Tillinghast (he is getting married on Saturday), there was a Great Earthquake offshore of Chiapas, Mexico. This is one of four M 8 or greater earthquakes ever recorded along the subduction zone forming the Middle American Trench. There has recently been some seismic activity to the east of this current M 8.1 earthquake. These earthquakes happened near the boundary between the North America (NAP) and Caribbean (CP) upper plates. I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I also include USGS epicenters from 1917-2017 for magnitudes M ≥ 8.0. I include fault plane solutions for the 1985 and 1995 earthquakes (along with the MMI contours for those earthquakes, see below for a discussion of MMI contours). TSUNAMI THREAT FORECAST ESTIMATED TIMES OF ARRIVAL TSUNAMI OBSERVATIONS
Marine magnetic anomalies and fracture zones that constrain tectonic reconstructions such as those shown in Figure 4 (ages of anomalies are keyed to colors as explained in the legend; all anomalies shown are from University of Texas Institute for Geophysics PLATES [2000] database): (1) Boxed area in solid blue line is area of anomaly and fracture zone picks by Leroy et al. (2000) and Rosencrantz (1994); (2) boxed area in dashed purple line shows anomalies and fracture zones of Barckhausen et al. (2001) for the Cocos plate; (3) boxed area in dashed green line shows anomalies and fracture zones from Wilson and Hey (1995); and (4) boxed area in red shows anomalies and fracture zones from Wilson (1996). Onland outcrops in green are either the obducted Cretaceous Caribbean large igneous province, including the Siuna belt, or obducted ophiolites unrelated to the large igneous province (Motagua ophiolites). The magnetic anomalies and fracture zones record the Cenozoic relative motions of all divergent plate pairs infl uencing the Central American subduction zone (Caribbean, Nazca, Cocos, North America, and South America). When incorporated into a plate model, these anomalies and fracture zones provide important constraints on the age and thickness of subducted crust, incidence angle of subduction, and rate of subduction for the Central American region. MCSC—Mid-Cayman Spreading Center.
Generalized tectonic map of the study area. Transparent red zones show the location of active volcanic belts in México: CMVB — Central Mexican Volcanic Belt, MCVA — Modern Chiapanecan Volcanic Arc. Transparent gray area: the extinct Sierra Madre Miocenic Arc. Orange stars are the El Chichón and San Martin active volcanoes. EPR — East Pacific Rise. MAT — Middle American Trench. Right black dashed line with a question mark is the hypothetical prolongation of Polochic–Montagua fault system which represents the limit between North America (NAM) and Caribbean plates. White dashed line is the onshore prolongation of Tehuantepec Ridge. Onshore white contours represent the slab isodepths (Bravo et al., 2004; Pardo and Suarez, 1995). Arrows show convergence velocities between the Cocos and North American plates (DeMets et al., 1994). TR1 and TR2 are the cross-sections where we calculate the thermal structure across TR (Fig. 10). Cocos plate ages are from Manea et al. (2005). White line dashed squares show the location of magnetic and gravity maps (Figs. 3, 5, 6). Blue dots represent continental heat-flow measurements (mW/m2) from Ziagos et al. (1985).
The study area encompasses Guerrero and Oaxaca states of Mexico. Shaded ellipse-like areas annotated with the years are rupture areas of the most recent major thrust earthquakes (M≥6.5) in the Mexican subduction zone. Triangles show locations of permanent GPS stations. Small hexagons indicate campaign GPS sites. Arrows are the Cocos-North America convergence vectors from NUVEL-1A model (DeMets et al., 1994). Double head arrow shows the extent of the Guerrero seismic gap. Solid and dashed curves annotated with negative numbers show the depth in km down to the surface of subducting Cocos plate (modified from Pardo and Su´arez, 1995, using the plate interface configuration model for the Central Oaxaca from this study, the model for Guerrero from Kostoglodov et al. (1996), and the last seismological estimates in Chiapas by Bravo et al. (2004). MAT, Middle America trench.
Inverted coupling coefficients along the MFe, MFc, MATch and MATgs, and residual velocities for best-fitting 3B model.
Same as Fig. 7 for best-fitting 4B model, with coupling along VAF fixed to 1.
Proposed model of faults kinematics and coupling along the Cocos slab interface, revised from Lyon-Caen et al. (2006). Numbers are velocities relative to CA plate in mmyr−1. Focal mechanisms are for crustal earthquakes (depth ≤30 km) since 1976, from CMT Harvard catalogue.
Earthquake Report Lite: M 7.0 near Acapulco, Mexico
Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here.
Last afternoon (my time) there was an M 7.0 earthquake near Acapulco, Mexico. This event generated a tsunami, landslides, building damage, casualties (one fatality as I write this), and many emotions.
https://earthquake.usgs.gov/earthquakes/eventpage/us7000f93v/executive
I present my interpretive poster and a few figures. Read more about the tectonics of this region here, in a report for an M 7.4 earthquake in 2020.Below is my interpretive poster for this earthquake
I include some inset figures.
Tide Gage Data – Acapulco
Earthquake Intensity
Mexico | Central America
General Overview
Earthquake Reports
References:
Basic & General References
Specific References
Return to the Earthquake Reports page.
Earthquake Report (and Tsunami) Oaxaca, Mexico
The M 7.4 Oaxaca, Mexico Earthquake occurred along the megathrust fault interface (an “interplate” earthquake) based on our knowledge of the location of the fault, our calculation of the earthquake location, and the earthquake mechanisms prepared by seismologists (i.e. focal mechanisms or moment tensors).
I noticed that there is a down-first wave prior to the tsunami. This was observed at both stations (Acapulco and Salina Cruz). Dr. Costas Synolakis (USC) informed me that this is a well known phenomena called a “Leading Depression N-wave.” I mark the location of the Salina Cruz gage on the interpretive poster below.Below is my interpretive poster for this earthquake
I include some inset figures. Some of the same figures are located in different places on the larger scale map below.
Other Report Pages
Some Relevant Discussion and Figures
Earthquake Triggered Landslides
Mexico | Central America
General Overview
Earthquake Reports
Social Media
Calculated with @antandre71 pic.twitter.com/6hQFMV55EX
CERRADA LA CARRETERA OAXACA-ISTMO POR #Sismo 7.5
Waiting #sentinel1…
Surface projection of the slip distribution by @USGS superimposed on the quickly forward model on #ArcGis platform; white lines indicate #Slab2 depth boundaries. The seismogenic fault is also shown in 3D view map (magenta plain). pic.twitter.com/sm62EHfF8Z
This one below is a recording from SE England pic.twitter.com/6RRzPsHAGK
References:
Basic & General References
Specific References
Return to the Earthquake Reports page.
Earthquake Report: Guatemala and Mexico
https://earthquake.usgs.gov/earthquakes/eventpage/us2000jbub/executive
Offshore of Guatemala and Mexico, the Middle America trench is formed by the subduction of the oceanic Cocos plate beneath the North America and Caribbean plates.
To the east of Guatemala and Mexico, the North America and Caribbean plates are separated by a left lateral (sinistral) strike-slip plate boundary fault (that forms the Cayman Trough beneath the Caribbean Sea).
As this plate boundary comes onshore, this fault forms multiple splays, including the Polochi-Montagua fault. As this system trends westwards across Central America, it joins another strike-slip plate boundary associated with the subduction zone (the Volcanic Arc fault).
South of about 15°N, the relative plate motion between the Caribbean and Cocos plates is oblique (they are not moving towards each other in a direction perpendicular to the subduction zone fault). At plate boundaries where plate convergence is oblique (like also found in Sumatra), the strain is partitioned onto the subduction zone (for fault normal component of the relative plate motion) and a forearc sliver fault (for the fault parallel relative motion).
The Tehuantepec fracture zone (TFZ) is a major structure in the Cocos plate. Coincidentally, the strike-slip fault systems trend towards where the TFZ intersects the trench.
There is left-lateral offset of the seafloor across the TFZ so the crust is about 10 million years older on the north side of the eastern TFZ. This age offset changes the depth of the crust across the TFZ and also may affect the megathrust fault properties on either side of the TFZ.
In addition, the TFZ may have geological properties that also affect the fault properties when this part of the plate subducts (affecting where, when, and how the fault slips).
There are so many things going on, but I will mention one more thing. Something that also appears to be happening in this part of the subduction zone is that there may be gaps in the slab beneath the megathrust. If this is true (Mann, 2007), then there may be changes in slab pull tension along strike as a result of different widths of attached downgoing slab.Below is my interpretive poster for this earthquake
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.
Magnetic Anomalies
Age of Oceanic Lithosphere
I include some inset figures. Some of the same figures are located in different places on the larger scale map below.
In 2017 there was a series of large magnitude earthquakes in the region of today’s M=6.6 and further to the south. These quakes are highlighted in the posters above, notable are the 6 Jun M=6.9 and 22 Jun M=6.8. The first quake was a deep extensional event, followed by a thrust event (possibly triggered by the M=6.9). In addition, there was a M=6.9 extensional earthquake in 2014 that also may have been a player.
I presented an interpretive poster showing the zone of aftershocks associated with the June sequence. Later, in Sept, there was a M=8.2 extensional tsunamigenic earthquake to the north of the June sequence. If we look at the aftershock zone for the M=8.2 quake, it looks like a sausage link adjacent to the sausage link formed by the June aftershocks. mmmm veggie sausages.
However there was no megathrust earthquake in the area of the M=8.2 sequence.
Other Report Pages
Some Relevant Discussion and Figures
Geologic Fundamentals
Compressional:
Extensional:
Mexico | Central America
Earthquake Reports
Social Media
References:
in northern Central America: Geological Society of America Special Paper 428, p. 1–19, https://doi.org/10.1130/2007.2428(01).Return to the Earthquake Reports page.
Earthquake Report: 2018 Summary
However, our historic record is very short, so any thoughts about whether this year (or last, or next) has smaller (or larger) magnitude earthquakes than “normal” are limited by this small data set.
Here is a table of the earthquakes M ≥ 6.5.
Here is a plot showing the cumulative release of seismic energy. This summary is imperfect in several ways, but shows how only the largest earthquakes have a significant impact on the tally of energy release from earthquakes. I only include earthquakes M ≥ 6.5. Note how the M 7.5 Sulawesi earthquake and how little energy was released relative to the two M = 7.9 earthquakes.
Below is my summary poster for this earthquake year
This is a video that shuffles through the earthquake report posters of the year
2018 Earthquake Report Pages
Other Annual Summaries
2018 Earthquake Reports
General Overview of how to interact with these summaries
Background on the Earthquake Report posters
Magnetic Anomalies
2018.01.10 M 7.6 Cayman Trough
Based upon our knowledge of the plate tectonics of this region, I can interpret the fault plane solution for this earthquake. The M 7.6 earthquake was most likely a left-lateral strike-slip earthquake associated with the Swan fault.
2018.01.14 M 7.1 Peru
In the region of this M 7.1 earthquake, two large structures in the NP are the Nazca Ridge and the Nazca fracture zone. The Nazca fracture zone is a (probably inactive) strike-slip fault system. The Nazca Ridge is an over-thickened region of the NP, thickened as the NP moved over a hotspot located near Salas y Gomez in the Pacific Ocean east of Easter Island (Ray et al., 2012).
There are many papers that discuss how the ridge affects the shape of the megathrust fault here. The main take-away is that the NR is bull dozing into South America and the dip of the subduction zone is flat here. There is a figure below that shows the deviation of the subducting slab contours at the NR.
Well, I missed looking further into a key update paper and used figures from an older paper on my interpretive poster yesterday. Thanks to Stéphane Baize for pointing this out! Turns out, after their new analyses, the M 7.1 earthquake was in a region of higher seismogenic coupling, rather than low coupling (as was presented in my first poster).
Also, Dr. Robin Lacassin noticed (as did I) the paucity of aftershocks from yesterday’s M 7.1. This was also the case for the carbon copy 2013 M 7.1 earthquake (there was 1 M 4.6 aftershock in the weeks following the M 7.1 earthquake on 2013.09.25; there were a dozen M 1-2 earthquakes in Nov. and Dec. of 2013, but I am not sure how related they are to the M 7.1 then). I present a poster below with this in mind. I also include below a comparison of the MMI modeled estimates. The 2013 seems to have possibly generated more widespread intensities, even though that was a deeper earthquake.
2018.01.23 M 7.9 Gulf of Alaska
This is strange because the USGS fault plane is oriented east-west, leading us to interpret the fault plane solution (moment tensor or focal mechanism) as a left-lateral strike-slip earthquake. So, maybe this earthquake is a little more complicated than first presumed. The USGS fault model is constrained by seismic waves, so this is probably the correct fault (east-west).
I prepared an Earthquake Report for the 1964 Good Friday Earthquake here.
So, that being said, here is the animation I put together. I used the USGS query tool to get earthquakes from 1/22 until now, M ≥ 1.5. I include a couple inset maps presented in my interpretive posters. The music is copyright free. The animations run through twice.
Here is a screenshot of the 14 MB video embedded below. I encourage you to view it in full screen mode (or download it).
2018.02.16 M 7.2 Oaxaca, Mexico
The SSN has a reported depth of 12 km, further supporting evidence that this earthquake was in the North America plate.
This region of the subduction zone dips at a very shallow angle (flat and almost horizontal).
There was also a sequence of earthquakes offshore of Guatemala in June, which could possibly be related to the M 8.1 earthquake. Here is my earthquake report for the Guatemala earthquake.
The poster also shows the seismicity associated with the M 7.6 earthquake along the Swan fault (southern boundary of the Cayman trough). Here is my earthquake report for the Guatemala earthquake.2018.02.25 M 7.5 Papua New Guinea
This M 7.5 earthquake (USGS website) occurred along the Papua Fold and Thrust Belt (PFTB), a (mostly) south vergent sequence of imbricate thrust faults and associated fold (anticlines). The history of this PFTB appears to be related to the collision of the Australia plate with the Caroline and Pacific plates, the delamination of the downgoing oceanic crust, and then associated magmatic effects (from decompression melting where the overriding slab (crust) was exposed to the mantle following the delamination). More about this can be found in Cloos et al. (2005).
The aftershocks are still coming in! We can use these aftershocks to define where the fault may have slipped during this M 7.5 earthquake. As I mentioned yesterday in the original report, it turns out the fault dimension matches pretty well with empirical relations between fault length and magnitude from Wells and Coppersmith (1994).
The mapped faults in the region, as well as interpreted seismic lines, show an imbricate fold and thrust belt that dominates the geomorphology here (as well as some volcanoes, which are probably related to the slab gap produced by crust delamination; see Cloos et al., 2005 for more on this). I found a fault data set and include this in the aftershock update interpretive poster (from the Coordinating Committee for Geoscience Programmes in East and Southeast Asia, CCOP).
I initially thought that this M 7.5 earthquake was on a fault in the Papuan Fold and Thrust Belt (PFTB). Mark Allen pointed out on twitter that the ~35km hypocentral depth is probably too deep to be on one of these “thin skinned” faults (see Social Media below). Abers and McCaffrey (1988) used focal mechanism data to hypothesize that there are deeper crustal faults that are also capable of generating the earthquakes in this region. So, I now align myself with this hypothesis (that the M 7.5 slipped on a crustal fault, beneath the thin skin deformation associated with the PFTB. (thanks Mark! I had downloaded the Abers paper but had not digested it fully.2018.03.08 M 6.8 New Ireland
The main transform fault (Weitin fault) is ~40 km to the west of the USGS epicenter. There was a very similar earthquake on 1982.08.12 (USGS website).
This earthquake is unrelated to the sequence occurring on the island of New Guinea.
Something that I rediscovered is that there were two M 8 earthquakes in 1971 in this region. This testifies that it is possible to have a Great earthquake (M ≥ 8) close in space and time relative to another Great earthquake. These earthquakes do not have USGS fault plane solutions, but I suspect that these are subduction zone earthquakes (based upon their depth).
This transform system is capable of producing Great earthquakes too, as evidenced by the 2000.11.16 M 8.0 earthquake (USGS website). This is another example of two Great earthquakes (or almost 2 Great earthquakes, as the M 7.8 is not quite a Great earthquake) are related. It appears that the M 8.0 earthquake may have triggered teh M 7.8 earthquake about 3 months later (however at first glance, it seemed to me like the strike-slip earthquake might not increase the static coulomb stress on the subduction zone, but I have not spent more than half a minute thinking about this).Main Interpretive Poster with emag2
Earthquakes M≥ 6.5 with emag2
2018.03.26 M 6.6 New Britain
Today’s M 6.6 earthquake happened close in proximity to a M 6.3 from 2 days ago and a M 5.6 from a couple weeks ago. The M 5.6 may be related (may have triggered these other earthquakes), but this region is so active, it might be difficult to distinguish the effects from different earthquakes. The M 5.6 is much deeper and looks like it was in the downgoing Solomon Sea plate. It is much more likely that the M 6.3 and M 6.6 are related (I interpret that the M 6.3 probably triggered the M 6.6, or that M 6.3 was a foreshock to the M 6.6, given they are close in depth). Both M 6.3 and M 6.6 are at depths close to the depth of the subducting slab (the megathrust fault depth) at this location. So, I interpret these to be subduction zone earthquakes.
2018.03.26 M 6.9 New Britain
2018.04.02 M 6.8 Bolivia
We are still unsure what causes an earthquake at such great a depth. The majority of earthquakes happen at shallower depths, caused largely by the frictional between differently moving plates or crustal blocks (where earth materials like the crust behave with brittle behavior and not elastic behavior). Some of these shallow earthquakes are also due to internal deformation within plates or crustal blocks.
As plates dive into the Earth at subduction zones, they undergo a variety of changes (temperature, pressure, stress). However, because people cannot directly observe what is happening at these depths, we must rely on inferences, laboratory analogs, and other indirect methods to estimate what is going on.
So, we don’t really know what causes earthquakes at the depth of this Bolivia M 6.8 earthquake. Below is a review of possible explanations as provided by Thorne Lay (UC Santa Cruz) in an interview in response to the 2013 M 8.3 Okhotsk Earthquake.
2018.05.04 M 6.9 Hawai’i
Hawaii is an active volcanic island formed by hotspot volcanism. The Hawaii-Emperor Seamount Chain is a series of active and inactive volcanoes formed by this process and are in a line because the Pacific plate has been moving over the hotspot for many millions of years.
Southeast of the main Kilauea vent, the Pu‘u ‘Ö‘ö crater saw an elevation of lava into the crater, leading to overtopping of the crater (on 4/30/2018). Seismicity migrated eastward along the ERZ. This morning, there was a M 5.0 earthquake in the region of the Hilina fault zone (HFZ). I was getting ready to write something up, but I had other work that I needed to complete. Then, this evening, there was a M 6.9 earthquake between the ERZ and the HFZ.
There have been earthquakes this large in this region in the past (e.g. the 1975.1.29 M 7.1 earthquake along the HFZ). This earthquake was also most likely related to magma injection (Ando, 1979). The 1975 M 7.1 earthquake generated a small tsunami (Ando, 1979). These earthquakes are generally compressional in nature (including the earthquakes from today).
Today’s earthquake also generated a tsunami as recorded on tide gages throughout Hawaii. There is probably no chance that a tsunami will travel across the Pacific to have a significant impact elsewhere.Temblor Reports:
2018.05.05 Pele, the Hawai’i Goddess of Fire, Lightning, Wind, and Volcanoes
2018.05.06 Pele, la Diosa Hawaiana del Fuego, los Relámpagos, el Viento y los Volcanes de Hawái
2018.08.05 M 6.9 Lombok, Indonesia
However, it is interesting because the earthquake sequence from last week (with a largest earthquake with a magnitude of M 6.4) were all foreshocks to this M 6.9. Now, technically, these were not really foreshocks. The M 6.4 has an hypocentral (3-D location) depth of ~6 km and the M 6.9 has an hypocentral depth of ~31 km. These earthquakes are not on the same fault, so I would interpret that the M 6.9 was triggered by the sequence from last week due to static coulomb changes in stress on the fault that ruptured. Given the large difference in depths, the uncertainty for these depths is probably not sufficient to state that they may be on the same fault (i.e. these depths are sufficiently different that this difference is larger than the uncertainty of their locations).
I present a more comprehensive analysis of the tectonics of this region in my earthquake report for the M 6.4 earthquake here. I especially address the historic seismicity of the region there. This M 6.9 may have been on the Flores thrust system, while the earthquakes from last week were on the imbricate thrust faults overlying the Flores Thrust. See the map from Silver et al. (1986) below. I include the same maps as in my original report, but after those, I include the figures from Koulani et al. (2016) (the paper is available on researchgate).2018.08.15 M 6.6 Aleutians
The Andreanof Islands is one of the most active parts of the Aleutian Arc. There have been many historic earthquakes here, some of which have been tsunamigenic (in fact, the email that notified me of this earthquake was from the ITIC Tsunami Bulletin Board).
Possibly the most significant earthquake was the 1957 Andreanof Islands M 8.6 Great (M ≥ 8.0) earthquake, though the 1986 M 8.0 Great earthquake is also quite significant. As was the 1996 M 7.9 and 2003 M 7.8 earthquakes. Lest we forget smaller earthquakes, like the 2007 M 7.2. So many earthquakes, so little time.2018.08.18 M 8.2 Fiji
This earthquake is one of the largest earthquakes recorded historically in this region. I include the other Large and Great Earthquakes in the posters below for some comparisons.
Today’s earthquake has a Moment Magnitude of M = 8.2. The depth is over 550 km, so is very very deep. This region has an historic record of having deep earthquakes here. Here is the USGS website for this M 8.2 earthquake. While I was writing this, there was an M 6.8 deep earthquake to the northeast of the M 8.2. The M 6.8 is much shallower (about 420 km deep) and also a compressional earthquake, in contrast to the extensional M 8.2.
This M 8.2 earthquake occurred along the Tonga subduction zone, which is a convergent plate boundary where the Pacific plate on the east subducts to the west, beneath the Australia plate. This subduction zone forms the Tonga trench.2018.08.19 M 6.9 Lombok, Indonesia
Today there was an M 6.3 soon followed by an M 6.9 earthquake (and a couple M 5.X quakes).
These earthquakes have been occurring along a thrust fault system along the northern portion of Lombok, Indonesia, an island in the magamatic arc related to the Sunda subduction zone. The Flores thrust fault is a backthrust to the subduction zone. The tectonics are complicated in this region of the world and there are lots of varying views on the tectonic history. However, there has been several decades of work on the Flores thrust (e.g. Silver et al., 1986). The Flores thrust is an east-west striking (oriented) north vergent (dipping to the south) thrust fault that extends from eastern Java towards the Islands of Flores and Timor. Above the main thrust fault are a series of imbricate (overlapping) thrust faults. These imbricate thrust faults are shallower in depth than the main Flores thrust.
The earthquakes that have been happening appear to be on these shallower thrust faults, but there is a possibility that they are activating the Flores thrust itself. Perhaps further research will illuminate the relations between these shallower faults and the main player, the Flores thrust.
2018.08.21 M 7.3 Venezuela
The northeastern part of Venezuela lies a large strike-slip plate boundary fault, the El Pilar fault. This fault is rather complicated as it strikes through the region. There are thrust faults and normal faults forming ocean basins and mountains along strike.
Many of the earthquakes along this fault system are strike-slip earthquakes (e.g. the 1997.07.09 M 7.0 earthquake which is just to the southwest of today’s temblor. However, today’s earthquake broke my immediate expectations for strike-slip tectonics. There is a south vergent (dipping to the north) thrust fault system that strikes (is oriented) east-west along the Península de Paria, just north of highway 9, east of Carupano, Venezuela. Audenard et al. (2000, 2006) compiled a Quaternary Fault database for Venezuela, which helps us interpret today’s earthquake. I suspect that this earthquake occurred on this thrust fault system. I bet those that work in this area even know the name of this fault. However, looking at the epicenter and the location of the thrust fault, this is probably not on this thrust fault. When I initially wrote this report, the depth was much shallower. Currently, the hypocentral (3-D location) depth is 123 km, so cannot be on that thrust fault.
The best alternative might be the subduction zone associated with the Lesser Antilles.2018.08.24 M 7.1 Peru
While doing my lit review, I found the Okal and Bina (1994) paper where they use various methods to determine focal mechanisms for the some deep earthquakes in northern Peru. More about focal mechanisms below. These authors created focal mechanisms for the 1921 and 1922 deep earthquakes so they could lean more about the 1970 deep earthquake. Their seminal work here forms an important record of deep earthquakes globally. These three earthquakes are all extensional earthquakes, similar to the other deep earthquakes in this region. I label the 1921 and 1922 earthquakes a couplet on the poster.
There was also a pair of earthquakes that happened in November, 2015. These two earthquakes happened about 5 minutes apart. They have many similar characteristics, suggest that they slipped similar faults, if not the same fault. I label these as doublets also.
So, there may be a doublet companion to today’s M 7.1 earthquake. However, there may be not. There are examples of both (single and doublet) and it might not really matter for 99.99% of the people on Earth since the seismic hazard from these deep earthquakes is very low.
Other examples of doublets include the 2006 | 2007 Kuril Doublets (Ammon et al., 2008) and the 2011 Kermadec Doublets (Todd and Lay, 2013).2018.09.05 M 6.6 Hokkaido, Japan
This earthquake is in an interesting location. to the east of Hokkaido, there is a subduction zone trench formed by the subduction of the Pacific plate beneath the Okhotsk plate (on the north) and the Eurasia plate (to the south). This trench is called the Kuril Trench offshore and north of Hokkaido and the Japan Trench offshore of Honshu.
One of the interesting things about this region is that there is a collision zone (a convergent plate boundary where two continental plates are colliding) that exists along the southern part of the island of Hokkaido. The Hidaka collision zone is oriented (strikes) in a northwest orientation as a result of northeast-southwest compression. Some suggest that this collision zone is no longer very active, however, there are an abundance of active crustal faults that are spatially coincident with the collision zone.
Today’s M 6.6 earthquake is a thrust or reverse earthquake that responded to northeast-southwest compression, just like the Hidaka collision zone. However, the hypocentral (3-D) depth was about 33 km. This would place this earthquake deeper than what most of the active crustal faults might reach. The depth is also much shallower than where we think that the subduction zone megathrust fault is located at this location (the fault formed between the Pacific and the Okhotsk or Eurasia plates). Based upon the USGS Slab 1.0 model (Hayes et al., 2012), the slab (roughly the top of the Pacific plate) is between 80 and 100 km. So, the depth is too shallow for this hypothesis (Kuril Trench earthquake) and the orientation seems incorrect. Subduction zone earthquakes along the trench are oriented from northwest-southweast compression, a different orientation than today’s M 6.6.
So today’s M 6.6 earthquake appears to have been on a fault deeper than the crustal faults, possibly along a deep fault associated with the collision zone. Though I am not really certain. This region is complicated (e.g. Kita et al., 2010), but there are some interpretations of the crust at this depth range (Iwasaki et al., 2004) shown in an interpreted cross section below.Temblor Reports:
2018.09.06 Violent shaking triggers massive landslides in Sapporo Japan earthquake
2018.09.09 M 6.9 Kermadec
This earthquake was quite deep, so was not expected to generate a significant tsunami (if one at all).
There are several analogies to today’s earthquake. There was a M 7.4 earthquake in a similar location, but much deeper. These are an interesting comparison because the M 7.4 was compressional and the M 6.9 was extensional. There is some debate about what causes ultra deep earthquakes. The earthquakes that are deeper than about 40-50 km are not along subduction zone faults, but within the downgoing plate. This M 6.9 appears to be in a part of the plate that is bending (based on the Benz et al., 2011 cross section). As plates bend downwards, the upper part of the plate gets extended and the lower part of the plate experiences compression.2018.09.28 M 7.5 Sulawesi
This area of Indonesia is dominated by a left-lateral (sinistral) strike-slip plate boundary fault system. Sulawesi is bisected by the Palu-Kola / Matano fault system. These faults appear to be an extension of the Sorong fault, the sinistral strike-slip fault that cuts across the northern part of New Guinea.
There have been a few earthquakes along the Palu-Kola fault system that help inform us about the sense of motion across this fault, but most have maximum magnitudes mid M 6.
GPS and block modeling data suggest that the fault in this area has a slip rate of about 40 mm/yr (Socquet et al., 2006). However, analysis of offset stream channels provides evidence of a lower slip rate for the Holocene (last 12,000 years), a rate of about 35 mm/yr (Bellier et al., 2001). Given the short time period for GPS observations, the GPS rate may include postseismic motion earlier earthquakes, though these numbers are very close.
Using empirical relations for historic earthquakes compiled by Wells and Coppersmith (1994), Socquet et al. (2016) suggest that the Palu-Koro fault system could produce a magnitude M 7 earthquake once per century. However, studies of prehistoric earthquakes along this fault system suggest that, over the past 2000 years, this fault produces a magnitude M 7-8 earthquake every 700 years (Bellier et al., 2006). So, it appears that this is the characteristic earthquake we might expect along this fault.
Most commonly, we associate tsunamigenic earthquakes with subduction zones and thrust faults because these are the types of earthquakes most likely to deform the seafloor, causing the entire water column to be lifted up. Strike-slip earthquakes can generate tsunami if there is sufficient submarine topography that gets offset during the earthquake. Also, if a strike-slip earthquake triggers a landslide, this could cause a tsunami. We will need to wait until people take a deeper look into this before we can make any conclusions about the tsunami and what may have caused it.
My 2018.10.01 BC Newshour Interview
InSAR Analysis
Interferometric SAR (InSAR) utilizes two separate SAR data sets to determine if the ground surface has changed over time, the time between when these 2 data sets were collected. More about InSAR can be found here and here. Explaining the details about how these data are analyzed is beyond the scope of this report. I rely heavily on the expertise of those who do this type of analysis, for example Dr. Eric Fielding.
M 7.5 Landslide Model vs. Observation Comparison
Until these landslides are analyzed and compared with regions that did not fail in slope failure, we will not be able to reconstruct what happened… why some areas failed and some did not.
There are landslide slope stability and liquefaction susceptibility models based on empirical data from past earthquakes. The USGS has recently incorporated these types of analyses into their earthquake event pages. More about these USGS models can be found on this page.
I prepared some maps that compare the USGS landslide and liquefaction probability maps. Below I present these results along with the MMI contours. I also include the faults mapped by Wilkinson and Hall (2017). Shown are the cities of Donggala and Palu. Also shown are the 2 tide gage locations (Pantoloan Port – PP and Mumuju – M). I also used post-earthquake satellite imagery to outline the largest landslides in Palu Valley, ones that appear to be lateral spreads.
Temblor Reports:
2018.09.28 The Palu-Koro fault ruptures in a M=7.5 quake in Sulawesi, Indonesia, triggering a tsunami and likely more shocks
2018.10.03 Tsunami in Sulawesi, Indonesia, triggered by earthquake, landslide, or both
2018.10.16 Coseismic Landslides in Sulawesi, Indonesia
2018.10.10 M 7.0 New Britain, PNG
The subduction zone forms the New Britain Trench with an axis that trends east-northeast. To the east of New Britain, the subduction zone bends to the southeast to form the San Cristobal and South Solomon trenches. Between these two subduction zones is a series of oceanic spreading ridges sequentially offset by transform (strike slip) faults.
Earthquakes along the megathrust at the New Britain trench are oriented with the maximum compressive stress oriented north-northwest (perpendicular to the trench). Likewise, the subduction zone megathrust earthquakes along the S. Solomon trench compress in a northeasterly direction (perpendicular to that trench).
There is also a great strike slip earthquake that shows that the transform faults are active.
This earthquake was too small and too deep to generate a tsunami.Temblor Reports:
2018.10.10 M 7.5 Earthquake in New Britain, Papua New Guinea
2018.10.22 M 6.8 Explorer plate
The Juan de Fuca plate is created at an oceanic spreading center called the Juan de Fuca Ridge. This spreading ridge is offset by several transform (strike-slip) faults. At the southern terminus of the JDF Ridge is the Blanco fault, a transtensional transform fault connecting the JDF and Gorda ridges.
At the northern terminus of the JDF Ridge is the Sovanco transform fault that strikes to the northwest of the JDF Ridge. There are additional fracture zones parallel and south of the Sovanco fault, called the Heck, Heckle, and Springfield fracture zones.
The first earthquake (M = 6.6) appears to have slipped along the Sovanco fault as a right-lateral strike-slip earthquake. Then the M 6.8 earthquake happened and, given the uncertainty of the location for this event, occurred on a fault sub-parallel to the Sovanco fault. Then the M 6.5 earthquake hit, back on the Sovanco fault.2018.10.25 M 6.8 Greece
Both of those earthquakes were right-lateral strike-slip earthquakes associated with the Kefallonia fault.
However, today’s earthquake sequence was further to the south and east of the strike-slip fault, in a region experiencing compression from the Ionian Trench subduction zone. But there is some overlap of these different plate boundaries, so the M 6.8 mainshock is an oblique earthquake (compressional and strike-slip). Based upon the sequence, I interpret this earthquake to be right-lateral oblique. I could be wrong.
Temblor Reports:
2018.10.26 Greek earthquake in a region of high seismic hazard
2018.11.08 M 6.8 Mid Atlantic Ridge (Jan Mayen fracture zone)
North of Iceland, the MAR is offset by many small and several large transform faults. The largest transform fault north of Iceland is called the Jan Mayen fracture zone, which is the location for the 2018.11.08 M = 6.8 earthquake.
2018.11.30 M 7.0 Alaska
During the 1964 earthquake, the downgoing Pacific plate slipped past the North America plate, including slip on “splay faults” (like the Patton fault, no relation, heheh). There was deformation along the seafloor that caused a transoceanic tsunami.
The Pacific plate has pre-existing zones of weakness related to fracture zones and spreading ridges where the plate formed and are offset. There was an earthquake in January 2016 that may have reactivated one of these fracture zones. This earthquake (M = 7.1) was very deep (~130 km), but still caused widespread damage.
The earthquake appears to have a depth of ~40 km and the USGS model for the megathrust fault (slab 2.0) shows the megathrust to be shallower than this earthquake. There are generally 2 ways that may explain the extensional earthquake: slab tension (the downgoing plate is pulling down on the slab, causing extension) or “bending moment” extension (as the plate bends downward, the top of the plate stretches out.Temblor Reports:
2018.11.30 Exotic M=7.0 earthquake strikes beneath Anchorage, Alaska
2018.12.11 What the Anchorage earthquake means for the Bay Area, Southern California, Seattle, and Salt Lake City
2018.12.05 M 7.5 New Caledonia
This part of the plate boundary is quite active and I have a number of earthquake reports from the past few years (see below, a list of earthquake reports for this region).
But the cool thing from a plate tectonics perspective is that there was a series of different types of earthquakes. At first view, it appears that there was a mainshock with a magnitude of M = 7.5. There was a preceding M 6.0 earthquake which may have been a foreshock.
The M 7.5 earthquake was an extensional earthquake. This may be due to either extension from slab pull or due to extension from bending of the plate. More on this later.
Following the M 7.5, there was an M 6.6 earthquake, however, this was a thrust or reverse (compressional) earthquake. The M 6.6 may have been in the upper plate or along the subduction zone megathrust fault, but we won’t know until the earthquake locations are better determined.
A similar sequence happened in October/November 2017. I prepared two reports for this sequence here and here. Albeit, in 2017, the thrust earthquake was first (2017.10.31 vs. 2017.11.19).
There have been some observations of tsunami. Below is from the Pacific Tsunami Warning Center.
2018.12.20 M 7.4 Bering Kresla
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 Okhotsk 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.
UPDATE #1
2018.12.29 M 7.0 Philippines
The earthquake was quite deep, which makes it less likely to cause damage to people and their belongings (e.g. houses and roads) and also less likely that the earthquake will trigger a trans-oceanic tsunami.
Here are the tidal data:
Geologic Fundamentals
Compressional:
Extensional:
Return to the Earthquake Reports page.
Earthquake Report: Gulf of California
The strike-slip faults are offset by oceanic spreading ridges. These spreading ridges are connected the East Pacific Rise to the south (and the Juan de Fuca Ridge and Gorda Rise to the north, via the San Andreas fault).
The geology of this region is much more complicated than this, but this is a good place to start when trying to understand the tectonics here. This M 6.3 earthquake happened on the southern boundary of the Guamas Basin, one of these pull apart basins.
So far there has been a single aftershock (M 4.5).
This M 6.3 earthquake appears to be pretty typical of this part of the Gulf. There have been several M 6 earthquakes in the last century. There have been a couple M 7 earthquakes, to the north.Below is my interpretive poster for this earthquake
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange) for the M 6.3 earthquakes, in addition to some relevant historic earthquakes.
I labeled the pull-apart basins in cyan and the faults in light orange. CdBF – C. de Ballenas fault; GF – Guaymas fault; CF – Carmen fault; FF – Farallon fault; PF – Pescadero fault; AF – Alcaron fault; AR – Alcaron Ridge (from Aragón-Arreola, M. and Martín-Barajas, A., 2007).
I include some inset figures.
segment.
Update
and Range Province; SMO—Sierra Madre Occidental; CP—Colorado Plateau: ITI—Isla Tiburón; IAG—Isla Ángel de la Guarda.
Social Media
Mw 6.3, Depth 10km
Focal Mechanism by @GFZ_Potsdam pic.twitter.com/Tah7HaHGCt
Mexico | Central America
General Overview
Earthquake Reports
References:
Earthquake Report: Chiapas Earthquake Update #2
Below is my interpretive poster for this earthquake
I include some inset figures in the poster.
Mexico | Central America
General Overview
Earthquake Reports
References:
Earthquake Report: Puebla, Mexico Update #1
I have been discussing this earthquake with other experts, both online (i.e. the twitterverse, where most convo happens these days) and offline. Many of these experts are presenting their interpretations of this earthquake as it may help us learn about plate tectonics. While many of us are interested in learning these technical details, I can only hope that we seek a similar goal, to reduce future suffering. Plate tectonics is a young science and we have an ultra short observation period (given that the recurrence of earthquakes can be centuries to millenia, it may take centuries or more to fully understand these processes).
Here I present a review of the material that I have seen in the past day and how I interpret these data. The main focus of the poster is a comparison of ground shaking for three earthquakes. Also of interest is the ongoing discussion about how the 2019.09.08 M 8.1 Chiapas Earthquake and this M 7.1 Puebla Earthquake relate to each other. My initial interpretation holds, that the temporal relations between these earthquakes is coincidental (but we now have the analysis to support this interpretation!).
Here are some additional reports found elsewhere on the internets (President G. W. Bush used to call the internet, “the internets”. heheh). This is the same list of sites from the original report.
Below is my interpretive poster for this earthquake
I include some inset figures in the poster.
Mexico | Central America
General Overview
Earthquake Reports
References:
Earthquake Report: Chiapas, Mexico Update #1
Below I present some figures about the Tehuantepec Seismic Gap (as before, but with additional figures). The impetus for this is two fold: (1) it is interesting for earthquake geologists as they consider earthquake recurrence patterns, globally and (2) that the M 8.1 earthquake was not a subduction zone earthquake and may have loaded the megathrust.
Here are some additional reports found elsewhere on the internets (President G. W. Bush used to call the internet, “the internets”. heheh).
Below is my interpretive poster for this earthquake
I outline the region of seismicity from June 2017. I include posters and links to the reports from that sequence below.
I also outline the region of the megathrust where the Tehuantepec Seismic Gap is located, generally in the region of the M 7.8 1902 megathrust earthquake. The Tehuantepec Ridge is a player in the regional tectonics as I discussed on my report earlier today.
I include some inset figures in the poster.
Seismicity: Mexico | Central America
General Overview
Earthquake Reports
References:
Earthquake Report: Chiapas, Mexico
This M 8.1 earthquake happened in a region of the subduction zone that is interpreted to have a higher coupling ratio than further to the south (higher proportion of the plate convergence rate is accumulated as elastic strain due to seismogenic coupling of the megathrust fault). Faults that are aseismic (fully slipping) have a coupling ratio of zero. The Polochic-Motagua fault zone marks this NAP-CP boundary. The recent seismicity offshore of Guatemala (June 2017) comprised a series of thrust earthquakes along the upper megathrust, along with some down-dip extensional faulting.
Tonight’s earthquake will be a very damaging and deadly earthquake and, based upon the shake map, possibly more damaging than either the 1985 or 1995 earthquakes. The 1985 earthquake caused severe damage in Mexico City. The PAGER alert shows an estimate of 34% probability for between 1000 and 10,000 fatalities. However, please read below about the PAGER alert and go to the USGS website about PAGER alerts (link below). These are just model based estimates of damage, so we won’t really know the damage until this is evaluated with “boots on the ground.” One might consider PAGER alerts to be the “armchair estimate” of damage. Thanks to Dr. Lori Dengler for reviewing my report (though any mistakes are only to be credited to me).
This M 8.1 earthquake is deeper than the megathrust fault and has an extensional moment tensor. This is not a megathrust earthquake, but is related to slip on a fault in the downgoing Cocos plate. At this depth, it may be due to bending in the downgoing oceanic lithosphere.
There is no danger of a tsunami here along the west coast of the U.S. West Coast, British Colombia, or Alaska. There have been some tsunami observations
Here are some additional reports found elsewhere on the internets (President G. W. Bush used to call the internet, “the internets”. heheh).
Below is my interpretive poster for this earthquake
I include some inset figures in the poster.
———————–
* TSUNAMI WAVES REACHING MORE THAN 3 METERS ABOVE THE TIDE
LEVEL ARE POSSIBLE ALONG SOME COASTS OF
MEXICO.
* TSUNAMI WAVES REACHING 0.3 TO 1 METERS ABOVE THE TIDE LEVEL
ARE POSSIBLE FOR SOME COASTS OF
AMERICAN SAMOA… ANTARCTICA… COOK ISLANDS… ECUADOR…
EL SALVADOR… FIJI… FRENCH POLYNESIA… GUATEMALA…
KIRIBATI… NEW ZEALAND… SAMOA… TOKELAU… TUVALU…
VANUATU… AND WALLIS AND FUTUNA.
* TSUNAMI WAVES ARE FORECAST TO BE LESS THAN 0.3 METERS ABOVE
THE TIDE LEVEL FOR THE COASTS OF
AUSTRALIA… CHILE… CHINA… CHUUK… COLOMBIA… COSTA
RICA… GUAM… HAWAII… HONDURAS… HOWLAND AND BAKER…
INDONESIA… JAPAN… JARVIS ISLAND… JOHNSTON ATOLL…
KERMADEC ISLANDS… KOSRAE… MALAYSIA… MARSHALL
ISLANDS… MIDWAY ISLAND… NAURU… NEW CALEDONIA…
NICARAGUA… NIUE… NORTHERN MARIANAS… NORTHWESTERN
HAWAIIAN ISLANDS… PALAU… PALMYRA ISLAND… PANAMA…
PAPUA NEW GUINEA… PERU… PHILIPPINES… PITCAIRN
ISLANDS… POHNPEI… RUSSIA… SOLOMON ISLANDS…
TAIWAN… TONGA… VIETNAM… WAKE ISLAND… AND YAP.
————————–
* ESTIMATED TIMES OF ARRIVAL -ETA- OF THE INITIAL TSUNAMI WAVE
FOR PLACES WITHIN THREATENED REGIONS ARE GIVEN BELOW. ACTUAL
ARRIVAL TIMES MAY DIFFER AND THE INITIAL WAVE MAY NOT BE THE
LARGEST. A TSUNAMI IS A SERIES OF WAVES AND THE TIME BETWEEN
WAVES CAN BE FIVE MINUTES TO ONE HOUR.
——————–
* THE FOLLOWING ARE TSUNAMI WAVE OBSERVATIONS FROM COASTAL
AND/OR DEEP-OCEAN SEA LEVEL GAUGES AT THE INDICATED
LOCATIONS. THE MAXIMUM TSUNAMI HEIGHT IS MEASURED WITH
RESPECT TO THE NORMAL TIDE LEVEL.
Seismicity: Mexico | Central America
General Overview
Earthquake Reports
References: