Well, last weekend I was working on a house, so did not have the time to write this up until now. (2023 update: the magnitude is now M 6.5)
https://earthquake.usgs.gov/earthquakes/eventpage/us700098qd/executive
The eastern Mediterranean Sea region is dominated by plate tectonics (no surprise, right?). The plate boundary fault system that is responsible for this earthquake near Crete is a convergent plate boundary called a subduction zone.
Convergent means that one plate is moving towards another plate. One of the largest plate boundary systems in the world is a convergent plate boundary that extends from between the north side of Australia and Indonesia, through southern Asia forming the Himalayan Mountains, through the Middle East, into Europe and west past the Mediterranean.
Near Crete the Africa plate is diving (northwards) beneath the Anatolia plate (a sliver of the Eurasia plate). The 2 May magnitude M 6.6 earthquake appears to have been an earthquake on the subduction zone megathrust fault interface (a subduction zone earthquake).
The earthquake was felt across the region with intensity as high as MMI 6 in Crete, to around MMI 4 in Cairo, Egypt.
The earthquake even caused a tsunami that was recorded at teh Lerapetra tide gage in Crete, Greece. The wave was small at about 40 cm peak to trough (measured vertically from the highest part of the wave, the peak, to the lowest part of the wave, the trough).
Here are the tide gage data downloaded from the IOC website here. The tsunami starts at around 13:00 hours.
Below is my interpretive poster for this earthquake
- I plot the seismicity from the past 6 months, with diameter representing magnitude (see legend). I include earthquake epicenters from 1920-2020 with magnitudes M ≥ 6.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.
- Some basic fundamentals of earthquake geology and plate tectonics can be found on the Earthquake Plate Tectonic Fundamentals page.
- In the center left is an inset map from Dilek and Sandovol (2009) that shows the tectonic plates and the plate boundary faults in the region. There is a blue star in the general location of the M 6.6 earthquake.
- In the upper right corner is a smaller scale view of the region with 6 months of seismicity plotted.
- In the lower right corner is a map that shows a model estimate of the shaking intensity from this M 6.6 earthquake.
- Above the intensity map is a map that shows earthquake mechanisms for historic earthquakes in the region.
- In the bottom center are seismic hazard and seismic risk maps for the European area. There is more about hazard and risk later in this report.
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
- Here is the tectonic map from Dilek and Sandvol (2009).
Tectonic map of the Aegean and eastern Mediterranean region showing the main plate boundaries, major suture zones, fault systems and tectonic units. Thick, white arrows depict the direction and magnitude (mm a21) of plate convergence; grey arrows mark the direction of extension (Miocene–Recent). Orange and purple delineate Eurasian and African plate affinities, respectively. Key to lettering: BF, Burdur fault; CACC, Central Anatolian Crystalline Complex; DKF, Datc¸a–Kale fault (part of the SW Anatolian Shear Zone); EAFZ, East Anatolian fault zone; EF, Ecemis fault; EKP, Erzurum–Kars Plateau; IASZ, Izmir–Ankara suture zone; IPS, Intra–Pontide suture zone; ITS, Inner–Tauride suture; KF, Kefalonia fault; KOTJ, Karliova triple junction; MM, Menderes massif; MS, Marmara Sea; MTR, Maras triple junction; NAFZ, North Anatolian fault zone; OF, Ovacik fault; PSF, Pampak–Sevan fault; TF, Tutak fault; TGF, Tuzgo¨lu¨ fault; TIP, Turkish–Iranian plateau (modified from Dilek 2006).
- This is the Wouldloper (2009) tectonic map of the Mediterranean Sea.
- Here is the large scale tectonic setting map (Taymaz et al., 2007) with their figure below.
Summary sketch map of the faulting and bathymetry in the Eastern Mediterranean region, compiled from our observations and those of Le Pichon & Angelier (1981), Taymaz (1990), Taymaz et al. (1990, 1991a, b); S¸arogˇlu et al. (1992), Papazachos et al. (1998), McClusky et al. (2000) and Tan & Taymaz (2006). Large black arrows show relative motions of plates with respect to Eurasia (McClusky et al. 2003). Bathymetry data are derived from GEBCO/97–BODC, provided by GEBCO (1997) and Smith & Sandwell (1997a, b). Shaded relief map derived from the GTOPO-30 Global Topography Data taken after USGS. NAF, North Anatolian Fault; EAF, East Anatolian Fault; DSF, Dead Sea Fault; NEAF, North East Anatolian Fault; EPF, Ezinepazarı Fault; PTF, Paphos Transform Fault; CTF, Cephalonia Transform Fault; PSF, Pampak–Sevan Fault; AS, Apsheron Sill; GF, Garni Fault; OF, Ovacık Fault; MT, Mus¸ Thrust Zone; TuF, Tutak Fault; TF, Tebriz Fault; KBF, Kavakbas¸ı Fault; MRF, Main Recent Fault; KF, Kagˇızman Fault; IF, Igˇdır Fault; BF, Bozova Fault; EF, Elbistan Fault; SaF, Salmas Fault; SuF, Su¨rgu¨ Fault; G, Go¨kova; BMG, Bu¨yu¨k Menderes Graben; Ge, Gediz Graben; Si, Simav Graben; BuF, Burdur Fault; BGF, Beys¸ehir Go¨lu¨ Fault; TF, Tatarlı Fault; SuF, Sultandagˇ Fault; TGF, Tuz Go¨lu¨ Fault; EcF, Ecemis¸ Fau; ErF, Erciyes Fault; DF, Deliler Fault; MF, Malatya Fault; KFZ, Karatas¸–Osmaniye Fault Zone.
- This figure shows GPS velocities in the region (Taymaz et al., 2007).
GPS horizontal velocities and their 95% confidence ellipses in a Eurasia-fixed reference frame for the period 1988–1997 superimposed on a shaded relief map derived from the GTOPO-30 Global Topography Data taken after USGS. Bathymetry data are derived from GEBCO/97–BODC, provided by GEBCO (1997) and Smith & Sandwell (1997a, b). Large arrows designate generalized relative motions of plates with respect to Eurasia (in mm a21) (recompiled after McClusky et al. 2000). NAF, North Anatolian Fault; EAF, East Anatolian Fault; DSF, Dead Sea Fault; NEAF, North East Anatolian Fault; EPF, Ezinepazarı Fault; CTF, Cephalonia Transform Fault; PTF, Paphos Transform Fault; CMT, Caucasus Main Thrust; MRF, Main Recent Fault.
- Finally their summary figure showing the tectonic regimes (Taymaz et al., 2007).
Schematic map of the principal tectonic settings in the Eastern Mediterranean. Hatching shows areas of coherent motion and zones of distributed deformation. Large arrows designate generalized regional motion (in mm a21) and errors (recompiled after McClusky et al. (2000, 2003). NAF, North Anatolian Fault; EAF, East Anatolian Fault; DSF, Dead Sea Fault; NEAF, North East Anatolian Fault; EPF, Ezinepazarı Fault; CTF, Cephalonia Transform Fault; PTF, Paphos Transform Fault.
- This is a tectonic summary figure from Kokkalas et al. (2006).
Simplified map showing the main structural features along the Hellenic arc and trench system, as well as the main active structures in the Aegean area. The mean GPS horizontal velocities in the Aegean plate, with respect to a Eurasia-fixed reference frame, are shown (after Kahle et al., 1998; McClusky et al., 2000). The lengths of vectors are
proportional to the amount of movement. The thick black arrows indicate the mean motion vectors of the plates. The polygonal areas on the map (dashed lines) define the approximate borders of the five different structural regions discussed in the text. The borders between structural regions are not straightforward, and wide transitional zones probably exist between them. The inset shows a schematic map with the geodynamic framework in the eastern Mediterranean area (modified from McClusky et al., 2000). DSF—Dead Sea fault; EAF—East Anatolia fault; HT—Hellenic trench; KFZ— Kefallonia fault zone; MRAC—Mediterranean Ridge accretionary complex; NAF—North Anatolia fault; NAT—North Aegean trough.
- The following three figures are from Dilek and Sandvol, 2006. The locations of the cross sections are shown on the map as orange lines. Cross section G-G’ is located in the region of today’s earthquake.
- Here is the map (Dilek and Sandvol, 2006). I include the figure caption below in blockquote.
- Here are cross sections A-D (Dilek and Sandvol, 2006). I include the figure caption below in blockquote.
- (A) Eastern Alps. The collision of Adria with Europe produced a bidivergent crustal architecture with both NNW- and SSE-directed nappe structures that involved Tertiary molasse deposits, with deep-seated thrust faults that exhumed lower crustal rocks. The Austro-Alpine units north of the Peri-Adriatic lineament represent the allochthonous outliers of the Adriatic upper crust tectonically resting on the underplating European crust. The Penninic ophiolites mark the remnants of the Mesozoic ocean basin (Meliata). The Oligocene granitoids between the Tauern window and the Peri-Adriatic lineament represent the postcollisional intrusions in the eastern Alps. Modified from Castellarin et al. (2006), with additional data from Coward and Dietrich (1989); Lüschen et al. (2006); Ortner et al. (2006).
- (B) Northern Apennines. Following the collision of Adria with the Apenninic platform and Europe in the late Miocene, the westward subduction of the Adriatic lithosphere and the slab roll-back (eastward) produced a broad extensional regime in the west (Apenninic back-arc extension) affecting the Alpine orogenic crust, and also a frontal thrust belt to the east. Lithospheric-scale extension in this broad back-arc environment above the west-dipping Adria lithosphere resulted in the development of a large boudinage structure in the European (Alpine) lithosphere. Modified from Doglioni et al. (1999), with data from Spakman and Wortel (2004); Zeck (1999).
- (C) Western Mediterranean–Southern Apennines–Calabria. The westward subduction of the Ionian seafloor as part of Adria since ca. 23 Ma and the associated slab roll-back have induced eastward-progressing extension and lithospheric necking through time, producing a series of basins. Rifting of Sardinia from continental Europe developed the Gulf of Lion passive margin and the Algero-Provencal basin (ca. 15–10 Ma), then the Vavilov and Marsili sub-basins in the broader Tyrrhenian basin to the east (ca. 5 Ma to present). Eastward-migrating lithospheric-scale extension and
necking and asthenospheric upwelling have produced locally well-developed alkaline volcanism (e.g., Sardinia). Slab tear or detachment in the Calabria segment of Adria, as imaged through seismic tomography (Spakman and Wortel, 2004), is probably responsible for asthenospheric upwelling and alkaline volcanism in southern Calabria and eastern Sicily (e.g., Mount Etna). Modified from Séranne (1999), with additional data from Spakman et al. (1993); Doglioni et al. (1999); Spakman and Wortel (2004); Lentini et al. (this volume). - (D) Southern Apennines–Albanides–Hellenides. Note the break where the Adriatic Sea is located between the western and eastern sections along this traverse. The Adria plate and the remnant Ionian oceanic lithosphere underlie the Apenninic-Maghrebian orogenic belt. The Alpine-Tethyan and Apulian platform units are telescoped along ENE-vergent thrust faults. The Tyrrhenian Sea opened up in the latest Miocene as a back-arc basin behind the Apenninic-Maghrebian mountain belt. The Aeolian volcanoes in the Tyrrhenian Sea represent the volcanic arc system in this subduction-collision zone environment. Modified from Lentini et al. (this volume). The eastern section of this traverse across the Albanides-Hellenides in the northern Balkan Peninsula shows a bidivergent crustal architecture, with the Jurassic Tethyan ophiolites (Mirdita ophiolites in Albania and Western Hellenic ophiolites in Greece) forming the highest tectonic nappe, resting on the Cretaceous and younger flysch deposits of the Adria affinity to the west and the Pelagonia affinity to the east. Following the emplacement of the Mirdita- Hellenic ophiolites onto the Pelagonian ribbon continent in the Early Cretaceous, the Adria plate collided with Pelagonia-Europe obliquely starting around ca. 55 Ma. WSW-directed thrusting, developed as a result of this oblique collision, has been migrating westward into the peri-Adriatic depression. Modified from Dilek et al. (2005).
- (E) Dinarides–Pannonian basin–Carpathians. The Carpathians developed as a result of the diachronous collision of the Alcapa and Tsia lithospheric blocks, respectively, with the southern edge of the East European platform during the early to middle Miocene (Nemcok et al., 1998; Seghedi et al., 2004). The Pannonian basin evolved as a back-arc basin above the eastward retreating European platform slab (Royden, 1988). Lithospheric-scale necking and boudinage development occurred synchronously with this extension and resulted in the isolation of continental fragments (e.g., the Apuseni mountains) within a broadly extensional Pannonian basin separating the Great Hungarian Plain and the Transylvanian subbasin. Steepening and tearing of the west-dipping slab may have caused asthenospheric flow and upwelling, decompressional melting, and alkaline volcanism (with an ocean island basalt–like mantle source) in the Eastern Carpathians. Modified from Royden (1988), with additional data from Linzer (1996); Nemcok et al. (1998); Doglioni et al. (1999); Seghedi et al. (2004).
- (F) Arabia-Eurasia collision zone and the Turkish-Iranian plateau. The collision of Arabia with Eurasia around 13 Ma resulted in (1) development of a thick orogenic crust via intracontinental convergence and shortening and a high plateau and (2) westward escape of a lithospheric block (the Anatolian microplate) away from the collision front. The Arabia plate and the Bitlis-Pütürge ribbon continent were probably amalgamated earlier (ca. the Eocene) via a separate collision event within the Neo-Tethyan realm. BSZ—Bitlis suture zone; EKP—Erzurum-Kars plateau. A slab break-off and the subsequent removal of the lithospheric mantle (lithospheric delamination) beneath the eastern Anatolian accretionary complex caused asthenospheric upwelling and extensive melting, leading to continental volcanism and regional uplift, which has contributed to the high mean elevation of the Turkish-Iranian plateau. The Eastern Turkey Seismic Experiment results have shown that the crustal thickness here is ~ 45–48 km and that the Turkish-Iranian plateau is devoid of mantle lithosphere. The collision-induced convergence has been accommodated by active diffuse north-south shortening and oblique-slip faults dispersing crustal blocks both to the west and the east. The late Miocene through Plio-Quaternary volcanism appears to have become more alkaline toward the south in time. The Pleistocene Karacadag shield volcano in the Arabian foreland represents a local fissure eruption associated with intraplate extension. Data from Pearce et al. (1990); Keskin (2003); Sandvol et al. (2003); S¸engör et al. (2003).
- (G) Africa-Eurasia collision zone and the Aegean extensional province. The African lithosphere is subducting beneath Eurasia at the Hellenic trench. The Mediterranean Ridge represents a lithospheric block between the Africa and Eurasian plate (Hsü, 1995). The Aegean extensional province straddles the Anatolide-Tauride and Sakarya continental blocks, which collided in the Eocene. NAF—North Anatolian fault. South-transported Tethyan ophiolite nappes were derived from the suture zone between these two continental blocks. Postcollisional granitic intrusions (Eocone and Oligo-Miocene, shown in red) occur mainly north of the suture zone and at the southern edge of the Sakarya continent. Postcollisional volcanism during the Eocene–Quaternary appears to have migrated southward and to have changed from calc-alkaline to alkaline in composition through time. Lithospheric-scale necking, reminiscent of the Europe-Apennine-Adria collision system, and associated extension are also important processes beneath the Aegean and have resulted in the exhumation of core complexes, widespread upper crustal attenuation, and alkaline and mid-ocean ridge basalt volcanism. Slab steepening and slab roll-back appear to have been at work resulting in subduction zone magmatism along the Hellenic arc.
Simplified tectonic map of the Mediterranean region showing the plate boundaries, collisional zones, and directions of extension and tectonic transport. Red lines A through G show the approximate profile lines for the geological traverses depicted in Figure 2. MHSZ—mid-Hungarian shear zone; MP—Moesian platform; RM—Rhodope massif; IAESZ— Izmir-Ankara-Erzincan suture zone; IPS—Intra-Pontide suture zone; ITS—inner Tauride suture zone; NAFZ—north Anatolian fault zone; KB—Kirsehir block; EKP—Erzurum-Kars plateau; TIP—Turkish-Iranian plateau.
Simplified tectonic cross-sections across various segments of the broader Alpine orogenic belt.
- Here is another cross section that shows the temporal evolution of the tectonics of this region in the area of cross section G-G’ above (Dilek and Sandvol, 2009).
Late Mesozoic–Cenozoic geodynamic evolution of the western Anatolian orogenic belt as a result of collisional
and extensional processes in the upper plate of north-dipping subduction zone(s) within the Tethyan realm. See text
for discussion.
- Here is the map showing the historic earthquake mechanisms from Jolivet et al. (2013).
Focal mechanisms of earthquakes over the Aegean Anatolian region.
Seismic Hazard and Seismic Risk
- These are the two maps shown in the map above, the GEM Seismic Hazard and the GEM Seismic Risk maps from Pagani et al. (2018) and Silva et al. (2018).
- The GEM Seismic Hazard Map:
- The Global Earthquake Model (GEM) Global Seismic Hazard Map (version 2018.1) depicts the geographic distribution of the Peak Ground Acceleration (PGA) with a 10% probability of being exceeded in 50 years, computed for reference rock conditions (shear wave velocity, VS30, of 760-800 m/s). The map was created by collating maps computed using national and regional probabilistic seismic hazard models developed by various institutions and projects, and by GEM Foundation scientists. The OpenQuake engine, an open-source seismic hazard and risk calculation software developed principally by the GEM Foundation, was used to calculate the hazard values. A smoothing methodology was applied to homogenise hazard values along the model borders. The map is based on a database of hazard models described using the OpenQuake engine data format (NRML). Due to possible model limitations, regions portrayed with low hazard may still experience potentially damaging earthquakes.
- Here is a view of the GEM seismic hazard map for Europe, the western Middle East, and Northern Africa.
- The GEM Seismic Risk Map:
- The Global Seismic Risk Map (v2018.1) presents the geographic distribution of average annual loss (USD) normalised by the average construction costs of the respective country (USD/m2) due to ground shaking in the residential, commercial and industrial building stock, considering contents, structural and non-structural components. The normalised metric allows a direct comparison of the risk between countries with widely different construction costs. It does not consider the effects of tsunamis, liquefaction, landslides, and fires following earthquakes. The loss estimates are from direct physical damage to buildings due to shaking, and thus damage to infrastructure or indirect losses due to business interruption are not included. The average annual losses are presented on a hexagonal grid, with a spacing of 0.30 x 0.34 decimal degrees (approximately 1,000 km2 at the equator). The average annual losses were computed using the event-based calculator of the OpenQuake engine, an open-source software for seismic hazard and risk analysis developed by the GEM Foundation. The seismic hazard, exposure and vulnerability models employed in these calculations were provided by national institutions, or developed within the scope of regional programs or bilateral collaborations.
- Here is a view of the GEM seismic risk map for Europe, the western Middle East, and Northern Africa.
- 2020.05.02 M 6.6 Crete, Greece
- 2020.01.24 M 6.7 Turkey
- 2019.11.26 M 6.4 Albania
- 2018.10.25 M 6.8 Greece
- 2017.07.20 M 6.7 Turkey
- 2017.06.12 M 6.3 Turkey/Greece
- 2016.10.30 M 6.6 Italy
- 2016.10.30 M 6.6 Italy Update #1
- 2016.10.28 M 5.8 Tyrrhenian Sea
- 2016.10.26 M 6.1 Italy
- 2016.10.16 M 5.3 Greece/Albania
- 2016.08.23 M 6.2 Italy
- 2016.01.24 M 6.1 Mediterranean
- 2015.11.17 M 6.5 Greece
- 2015.04.16 M 6.0 Crete
Europe
General Overview
Earthquake Reports
Social Media
- 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
- Basili R., G. Valensise, P. Vannoli, P. Burrato, U. Fracassi, S. Mariano, M.M. Tiberti, E. Boschi (2008), The Database of Individual Seismogenic Sources (DISS), version 3: summarizing 20 years of research on Italy’s earthquake geology, Tectonophysics, doi:10.1016/j.tecto.2007.04.014
- Brun, J.-P., Sokoutis, D., 2012. 45 m.y. of Aegean crust and mantle flow driven by trench retreat. Geol. Soc. Am., v. 38, p. 815–818.
- Caputo, R., Chatzipetros, A., Pavlides, S., and Sboras, S., 2012. The Greek Database of Seismogenic Sources (GreDaSS): state-of-the-art for northern Greece in Annals of Geophysics, v. 55, no. 5, doi: 10.4401/ag-5168
- Dilek, Y. and Sandvol, E., 2006. Collision tectonics of the Mediterranean region: Causes and consequences in Dilek, Y., and Pavlides, S., eds., Postcollisional tectonics and magmatism in the Mediterranean region and Asia: Geological Society of America Special Paper 409, p. 1–13
- DISS Working Group (2015). Database of Individual Seismogenic Sources (DISS), Version 3.2.0: A compilation of potential sources for earthquakes larger than M 5.5 in Italy and surrounding areas. http://diss.rm.ingv.it/diss/, Istituto Nazionale di Geofisica e Vulcanologia; DOI:10.6092/INGV.IT-DISS3.2.0.
- Ersoy, E.Y., Cemen, I., Helvaci, C., and Billor, Z., 2014. Tectono-stratigraphy of the Neogene basins in Western Turkey: Implications for tectonic evolution of the Aegean Extended Region in Tectonophysics v. 635, p. 33-58.
- Ganas, A., and T. Parsons (2009), Three-dimensional model of Hellenic Arc deformation and origin of the Cretan uplift, J. Geophys. Res., 114, B06404, doi:10.1029/2008JB005599
- Ganas, A., Oikonomou, I.A., and Tsimi, C., 2013. NOAFAULTS: A Digital Database for Active Faults in Greece in Bulletin of the Geological Society of Greece, v. XLVII, Proceedings fo the 13th International Cogfress, Chania, Sept, 2013
- Kokkalas, S., Xypolias, P., Koukouvelas, I., and Doutsos, T., 2006, Postcollisional contractional and extensional deformation in the Aegean region, in Dilek, Y., and Pavlides, S., eds., Postcollisional tectonics and magmatism in the Mediterranean region and Asia: Geological Society of America Special Paper 409, p. 97–123, doi: 10.1130/2006.2409(06)
- Papazachos, B.C., Papadimitrious, E.E., Kiratzi, A.A., Papazachos, C.B., and Louvari, E.k., 1998. Fault Plane Solutions in the Aegean Sea and the Surrounding Area and their Tectonic Implication, in Bollettino Di Geofisica Terorica Ed Applicata, v. 39, no. 3, p. 199-218.
- Taymaz, T. , Yilmaz, Y., and Dilek, Y., 2007. The geodynamics of the Aegean and Anatolia: introduction in TAYMAZ, T., YILMAZ, Y. & DILEK, Y. (eds) The Geodynamics of the Aegean and Anatolia. Geological Society, London, Special Publications, 291, 1–16. DOI: 10.1144/SP291.1 0305-8719/07
- Wouldloper, 2009. Tectonic map of southern Europe and the Middle East, showing tectonic structures of the western Alpide mountain belt. Only Alpine (tertiary) structures are shown.
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
This M 6.7 earthquake was the result of slip probably along a left-lateral strike-slip fault associated with the East Anatolia fault zone (EAF). The event was shallow and produced strong ground shaking in the region. https://earthquake.usgs.gov/earthquakes/eventpage/us60007ewc/executive As I write this, there have been about 5 building collapses and 22 deaths. The high number of deaths may be due to the building design used in the region. The EAF accommodates the relative plate motion between the Anatolia and Arabia plates. Because the northern motion of the Arabia plate is oblique to the plate boundary, the tectonic strain (deformation of the Earth) is proportioned on different fault types. We call this strain partitioning. The lateral strain is localized along the EAF in the form of strike-slip faults. The compressive strain formed the Southeast Anatolia fault zone, a series of imbricate thrust faults south and east of the EAF. Further to the west, this north-south compression results in the subduction of the Africa plate northwards beneath the Anatolia and Eurasia plates. This subduction forms the Hellenic trench. On the northern part of Turkey is bordered by a right-lateral strike-slip fault, the North Anatolia fault. Last year (2019) was the 20 year commemoration of the 1999 Izmit M 7.6 earthquake. The M 6.7 earthquake may have caused landslides or liquefaction in places, but the chances of this are modest at best. Geologists have studied the EAF and subdivided the fault into segments based on their mapping efforts. This M 6.7 is within the Pütürge segment of the EAF. If we look at the historic record of the EAF here, we find that the M 6.7 happened in a part of the fault that does not have an historic rupture. There was an earthquake in 1875 that appears to end to the north of the M 6.7 and there is an earthquake in 1893 that appears to terminate just to the south of the M 6.7.
Tectonic setting of continental extrusion in eastern Mediterranean. Anatolia-Aegean block escapes westward from Arabia-Eurasia collision zone, toward Hellenic subduction zone. Current motion relative to Eurasia (GPS [Global Positioning System] and SLR [Satellite Laser Ranging] velocity vectors, in mm/yr, from Reilinger et al., 1997). In Aegean, two deformation regimes are superimposed (Armijo et al., 1996): widespread, slow extension starting earlier (orange stripes, white diverging arrows), and more localized, fast transtension associated with later, westward propagation of North Anatolian fault (NAF). EAF—East Anatolian fault, K—Karliova triple junction, DSF—Dead Sea fault,NAT—North Aegean Trough, CR—Corinth Rift.Box outlines Marmara pull-apart region, where North Anatolian fault enters Aegean.
Tectonic map of the Aegean and eastern Mediterranean region showing the main plate boundaries, major suture zones, fault systems and tectonic units. Thick, white arrows depict the direction and magnitude (mm a21) of plate convergence; grey arrows mark the direction of extension (Miocene–Recent). Orange and purple delineate Eurasian and African plate affinities, respectively. Key to lettering: BF, Burdur fault; CACC, Central Anatolian Crystalline Complex; DKF, Datc¸a–Kale fault (part of the SW Anatolian Shear Zone); EAFZ, East Anatolian fault zone; EF, Ecemis fault; EKP, Erzurum–Kars Plateau; IASZ, Izmir–Ankara suture zone; IPS, Intra–Pontide suture zone; ITS, Inner–Tauride suture; KF, Kefalonia fault; KOTJ, Karliova triple junction; MM, Menderes massif; MS, Marmara Sea; MTR, Maras triple junction; NAFZ, North Anatolian fault zone; OF, Ovacik fault; PSF, Pampak–Sevan fault; TF, Tutak fault; TGF, Tuzgo¨lu¨ fault; TIP, Turkish–Iranian plateau (modified from Dilek 2006).
A: Tectonic map of the Aegean and Anatolian region showing the main active structures
C: GPS velocity field with a fixed Eurasia after Reilinger et al. (2010) D: the domain affected by distributed post-orogenic extension in the Oligocene and the Miocene and the stretching lineations in the exhumed metamorphic complexes.
E: The thick blue lines illustrate the schematized position of the slab at ~150 km according to the tomographic model of Piromallo and Morelli (2003), and show the disruption of the slab at three positions and possible ages of these tears discussed in the text. Velocity anomalies are displayed in percentages with respect to the reference model sp6 (Morelli and Dziewonski, 1993). Coloured symbols represent the volcanic centres between 0 and 3 Ma after Pe-Piper and Piper (2006). F: Seismic anisotropy obtained from SKS waves (blue bars, Paul et al., 2010) and Rayleigh waves (green and orange bars, Endrun et al., 2011). See also Sandvol et al. (2003). Blue lines show the direction of stretching in the asthenosphere, green bars represent the stretching in the lithospheric mantle and orange bars in the lower crust.
G: Focal mechanisms of earthquakes over the Aegean Anatolian region.
Schematic map of the principal tectonic settings in the Eastern Mediterranean. Hatching shows areas of coherent motion and zones of distributed deformation. Large arrows designate generalized regional motion (in mm a21) and errors (recompiled after McClusky et al. (2000, 2003). NAF, North Anatolian Fault; EAF, East Anatolian Fault; DSF, Dead Sea Fault; NEAF, North East Anatolian Fault; EPF, Ezinepazarı Fault; CTF, Cephalonia Transform Fault; PTF, Paphos Transform Fault.
The main fault systems of the AN–AR and TR–AF plate boundaries (modified from Sengor & Yılmaz 1981; Saroglu et al. 1992a, b; Westaway 2003; Emre et al. 2011a, b, c). Arrows indicate relative plate motions (McClusky et al. 2000). Abbreviations: AN, Anatolian microplate; AF, African plate; AR, Arabian plate; EU, Eurasian plate; NAFZ, North Anatolian Fault Zone; EAFZ, East Anatolian Fault Zone; DSFZ, Dead Sea Fault Zone; MF; Malatya Fault, TF, Tuzgo¨lu¨ fault; EF, Ecemis¸ fault; SATZ, Southeast Anatolian Thrust Zone; SS, southern strand of the EAFZ; NS, northern strand of the EAFZ.
Map of the East Anatolian strike-slip fault system showing strands, segments and fault jogs. Abbreviations: FS, fault Segment; RB, releasing bend; RS, releasing stepover; RDB, restraining double bend; RSB, restraining bend; PB, paired bend; (1) Du¨zic¸i–Osmaniye fault segment; (2) Erzin fault segment; (3) Payas fault segment; (4) Yakapınar fault segment; (5) C¸ okak fault segment; (6) Islahiye releasing bend; (7) Demrek restraining stepover; (8) Engizek fault zone; (9) Maras¸ fault zone.
Map of the (a) Palu and (b) Puturge segments of the East Anatolian fault. Abbreviations: LHRB, Lake Hazar releasing bend; PS, Palu segment; ES, Erkenek segment; H, hill; M, mountain; C, creek; (1) left lateral strike-slip fault; (2) normal fault; (3) reverse or thrust fault; (4) East Anatolian Fault; (5) Southeastern Anatolian Thrust Zone; (6) syncline;(7) anticline; (8) undifferentiated Holocene deposits; (9) undifferentiated Quaternary deposits; (10) landslide.
Surface ruptures produced by large earthquakes during the 19th and 20th centuries along the EAF. Data from Arpat (1971), Arpat and S¸arog˘lu (1972), Seymen and Aydın (1972), Ambraseys (1988), Ambraseys and Jackson (1998), Cetin et al. (2003), Herece (2008), Karabacak et al. (2011) and this study. Ruptured fault segments are highlighted.
The Global Seismic Hazard Map. Peak ground acceleration (pga) with a 10% chance of exceedance in 50 years is depicted in m/s2. The site classification is rock everywhere except Canada and the United States, which assume rock/firm soil site classifications. White and green correspond to low seismicity hazard (0%-8%g), yellow and orange correspond to moderate seismic hazard (8%-24%g), pink and dark pink correspond to high seismicity hazard (24%-40%g), and red and brown correspond to very high seismic hazard (greater than 40%g).
#EarthquakeReport for #depremElazig #Deprem #Earthquake in #Turkey Posters here I have tweeted before, now there is a report to explain them. I will update this report over time — Jason "Jay" R. Patton (@patton_cascadia) January 25, 2020 #EarthquakeReport pages can now be translated on https://t.co/fGEEJoACJA by simply clicking the language on the upper right. Uses Google Translate plugin "GTranslate: for WordPress" pic.twitter.com/hkLvy8VlAq — Jason "Jay" R. Patton (@patton_cascadia) January 26, 2020 Mw=6.7, TURKEY (Depth: 23 km), 2020/01/24 17:55:14 UTC – Full details here: https://t.co/rL50XD3kRs pic.twitter.com/uBFWy5q1hN — Earthquakes (@geoscope_ipgp) January 24, 2020 very strong and dangerous #earthquake #deprem along border of #Turkey and #Syria near #Elazig and #Malatya, felt widely including as far as #Istanbul and all over the Middle East Due to crisis situation in this region, impact estimates are difficult@ShakingEarth pic.twitter.com/Z1EPX0byHi — CATnews | Andreas M. Schäfer (@CATnewsDE) January 24, 2020 Map of felt reports received so far following the #earthquake M6.9 in Eastern Turkey 40 min ago pic.twitter.com/RXQx9Vkbw3 — EMSC (@LastQuake) January 24, 2020 Updated aftershock map of Jan 24 Mw6.8 #earthquake near Sivrice, Elazığ (eastern Turkey); >60km rupture along the Pütürge Segment of East Anatolia Fault Zone. Epicenters & focal mechanism from AFAD, active fault traces from MTA. Latest #Landsat8 image from Jan 22. pic.twitter.com/YR45rY6C1z — Sotiris Valkaniotis (@SotisValkan) January 25, 2020 Watch the waves from the M6.7 Turkey #earthquake roll across the USArray Transportable Array seismic network (https://t.co/RIcNz4bgWq ). #TurkeyEarthquake (THREAD) pic.twitter.com/UoxtaFVOQ6 — IRIS Earthquake Sci (@IRIS_EPO) January 24, 2020 Mw6.9 #earthquake in eastern #Turkey. Possible activation of Pütürge Segment, East Anatolia Fault. Preliminary locations of the epicenter, and GFZ focal mech. Fault map from Duman & Emre (2013) pic.twitter.com/iot7Rb58EA — Sotiris Valkaniotis (@SotisValkan) January 24, 2020 The @USGS generates finite fault models for larger EQs, where they try to reconstruct the quake. For the #TurkeyEarthquake, their prelim estimate is ~20 seconds of rupture, or 20 seconds just for the fault(s) to break. This is different from how long the earth shook. https://t.co/ln4bqol3OU — Alka Tripathy-Lang, PhD (@DrAlkaTrip) January 25, 2020 🗺 New map: [#EMSR423] Elazig: Grading Product, version 1, release 1, RTP Map #01 [v1, 1:30000] — Copernicus EMS (@CopernicusEMS) January 25, 2020 29 dead, 1,466 injured as massive quake of magnitude 6.8 rocks Turkey's Elazığ https://t.co/HaUzX5C8hR — Jason "Jay" R. Patton (@patton_cascadia) January 25, 2020 The Elazig earthquake from a 6-day Sentinel-1 pair (@ESA_EO @CopernicusData), processed using ISCE. The colour gradient from red to blue (instead of a sharp discontinuity) suggests that most of the slip occurred at depth and did not make it to the surface (intriguing…) pic.twitter.com/gZcqhs1wsq — Gareth Funning (@gfun) January 28, 2020
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.
Well, I was about to head to town and noticed a magnitude M = 5.0 earthquake in Greece. I thought to myself, I wonder if that is a foreshock. It was. 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 in one version.
Summary sketch map of the faulting and bathymetry in the Eastern Mediterranean region, compiled from our observations and those of Le Pichon & Angelier (1981), Taymaz (1990), Taymaz et al. (1990, 1991a, b); S¸arogˇlu et al. (1992), Papazachos et al. (1998), McClusky et al. (2000) and Tan & Taymaz (2006). Large black arrows show relative motions of plates with respect to Eurasia (McClusky et al. 2003). Bathymetry data are derived from GEBCO/97–BODC, provided by GEBCO (1997) and Smith & Sandwell (1997a, b). Shaded relief map derived from the GTOPO-30 Global Topography Data taken after USGS. NAF, North Anatolian Fault; EAF, East Anatolian Fault; DSF, Dead Sea Fault; NEAF, North East Anatolian Fault; EPF, Ezinepazarı Fault; PTF, Paphos Transform Fault; CTF, Cephalonia Transform Fault; PSF, Pampak–Sevan Fault; AS, Apsheron Sill; GF, Garni Fault; OF, Ovacık Fault; MT, Mus¸ Thrust Zone; TuF, Tutak Fault; TF, Tebriz Fault; KBF, Kavakbas¸ı Fault; MRF, Main Recent Fault; KF, Kagˇızman Fault; IF, Igˇdır Fault; BF, Bozova Fault; EF, Elbistan Fault; SaF, Salmas Fault; SuF, Su¨rgu¨ Fault; G, Go¨kova; BMG, Bu¨yu¨k Menderes Graben; Ge, Gediz Graben; Si, Simav Graben; BuF, Burdur Fault; BGF, Beys¸ehir Go¨lu¨ Fault; TF, Tatarlı Fault; SuF, Sultandagˇ Fault; TGF, Tuz Go¨lu¨ Fault; EcF, Ecemis¸ Fau; ErF, Erciyes Fault; DF, Deliler Fault; MF, Malatya Fault; KFZ, Karatas¸–Osmaniye Fault Zone.
GPS horizontal velocities and their 95% confidence ellipses in a Eurasia-fixed reference frame for the period 1988–1997 superimposed on a shaded relief map derived from the GTOPO-30 Global Topography Data taken after USGS. Bathymetry data are derived from GEBCO/97–BODC, provided by GEBCO (1997) and Smith & Sandwell (1997a, b). Large arrows designate generalized relative motions of plates with respect to Eurasia (in mm a21) (recompiled after McClusky et al. 2000). NAF, North Anatolian Fault; EAF, East Anatolian Fault; DSF, Dead Sea Fault; NEAF, North East Anatolian Fault; EPF, Ezinepazarı Fault; CTF, Cephalonia Transform Fault; PTF, Paphos Transform Fault; CMT, Caucasus Main Thrust; MRF, Main Recent Fault.
Schematic map of the principal tectonic settings in the Eastern Mediterranean. Hatching shows areas of coherent motion and zones of distributed deformation. Large arrows designate generalized regional motion (in mm a21) and errors (recompiled after McClusky et al. (2000, 2003). NAF, North Anatolian Fault; EAF, East Anatolian Fault; DSF, Dead Sea Fault; NEAF, North East Anatolian Fault; EPF, Ezinepazarı Fault; CTF, Cephalonia Transform Fault; PTF, Paphos Transform Fault.
Simplified map showing the main structural features along the Hellenic arc and trench system, as well as the main active structures in the Aegean area. The mean GPS horizontal velocities in the Aegean plate, with respect to a Eurasia-fixed reference frame, are shown (after Kahle et al., 1998; McClusky et al., 2000). The lengths of vectors are
General simplified structural map of Greece showing the main currently active structures in the five structural regions along the Hellenic Arc, as well as some main pre-existing lineaments. Insets illustrate the main structural features of each region and the period of activity of these structures (for further details see discussion). KFZ—Kefallonia Fault zone; MCL—Mid-Cycladic lineament; NAFZ—North Anatolia fault zone; NAT—North Aegean trough; PF—Pelagonian fault.
Schematic structural map of the central Hellenic Peninsula (Region II), with stress nets showing the orientation of principal stress axes. Stress net explanation as for Figure 3. Also included are cross-sections showing the geometry and kinematics of the External Hellenides in the area (A-A′) and the evolution of the synorogenic basin in the Paleros area (B-B′-B′′). AG—Abelon graben; ALG—Almyros graben; AMG—Amvrakikos graben; CG—Corinth graben; KB—Kymi basin; KF—Klenia fault zone; KFZ—Kefalonia fault zone; LF—Lapithas fault; MG—Megara graben; NG—Nedas graben; P—Parnitha area; PG—Pyrgos graben; PLB—Paleros basin; PTG—Patras graben; RG—Rio graben; S-A.G— Sperchios-Atalanti graben; SEG—South Evoikos graben; TB—Thiva basin; TG—Tithorea graben; TRG—Trihonis graben; VF—Vounargos fault.
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.
M6.8 #earthquake was a truly international quake felt up to 800 km in at least 8 countries (Malta, Libya, Italy, Greece, Macedonia, Albania, Bosnia, Bulgaria, Turkey) pic.twitter.com/gt2oVn5tLS — EMSC (@LastQuake) October 26, 2018 Todays strong #earthquake in Ionian Sea, #Greece, in an area of high seismicity near the junction of Hellenic Arc subduction and Cephalonia Transform Fault. NW-SE/WNW-ESE thrusting in accord. with local tectonic structures. Revised epicenter location (ML 6.4) from NOA. pic.twitter.com/0aqZhsTwOb — Sotiris Valkaniotis (@SotisValkan) October 26, 2018 Sea level changes following M6.8 earthquake are still being observed. Their amplitude remain hopefully limited. Do not go checking your boats there may be strong currents in harbors pic.twitter.com/0WGNjZGnd2 — EMSC (@LastQuake) October 26, 2018 Mw=7.0, IONIAN SEA (Depth: 15 km), 2018/10/25 22:54:51 UTC – Full details here: https://t.co/jXsx6jWEC9 pic.twitter.com/pZiAZvPeKx — Earthquakes (@geoscope_ipgp) October 25, 2018 strong #earthquake offshore Western #Greece near #Zakynthos, also felt in Southern #Italy and #Albania @LastQuake @Quake_Tracker @JuskisErdbeben pic.twitter.com/BF8W44NOjq — CATnews (@CATnewsDE) October 25, 2018 potential #tsunami after #earthquake offshore #Greece, still uncertain, yet a run-up of up to 1.5m can be possible locally, but mostly expected peak coastal wave heights of 20-50cm in Western Greece @LastQuake @JuskisErdbeben @Quake_Tracker4 pic.twitter.com/cejAwCoGgD — CATnews (@CATnewsDE) October 25, 2018 A 6.9 magnitude earthquake took place in Greece today. The epicenter of the quake was reported to be in the Iyon Sea, 250 km west to Athens. Users shared footage on social media of the quake, that happened in a depth of 16.6 km. pic.twitter.com/p4VrGUXT9y — EHA News (@eha_news) October 25, 2018 Faulting mechanism for the M6.8 Greece earthquake (blue beach ball; from the USGS) has a complex oblique mechanism. Past quakes in the area related to the Hellenic subduction thrust. The Cephalonia strike-slip fault lies a little further to the north. Tectonically complex region. pic.twitter.com/asklvRpB7V — Stephen Hicks 🇪🇺 (@seismo_steve) October 25, 2018
Return to the Earthquake Reports page.
Earthquake Report: East Anatolia fault zone
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
(black lines), the main sutures zones (thick violet or blue lines), the main thrusts in the Hellenides where they have not been reworked by later extension (thin blue lines), the North Cycladic Detachment (NCDS, in red) and its extension in the Simav Detachment (SD), the main metamorphic units and their contacts; AlW: Almyropotamos window; BD: Bey Daglari; CB: Cycladic Basement; CBBT: Cycladic Basement basal thrust; CBS: Cycladic Blueschists; CHSZ: Central Hellenic Shear Zone; CR: Corinth Rift; CRMC: Central Rhodope Metamorphic Complex; GT: Gavrovo–Tripolitza Nappe; KD: Kazdag dome; KeD: Kerdylion Detachment; KKD: Kesebir–Kardamos dome; KT: Kephalonia Transform Fault; LN: Lycian Nappes; LNBT: Lycian Nappes Basal Thrust; MCC: Metamorphic Core Complex; MG: Menderes Grabens; NAT: North Aegean Trough; NCDS: North Cycladic Detachment System; NSZ: Nestos Shear Zone; OlW: Olympos Window; OsW: Ossa Window; OSZ: Ören Shear Zone; Pel.: Peloponnese; ÖU: Ören Unit; PQN: Phyllite–Quartzite Nappe; SiD: Simav Detachment; SRCC: South Rhodope Core Complex; StD: Strymon Detachment; WCDS: West Cycladic Detachment System; ZD: Zaroukla Detachment. B: Seismicity. Earthquakes are taken from the USGS-NEIC database. Colour of symbols gives the depth (blue for shallow depths) and size gives the magnitude (from 4.5 to 7.6).
Seismic Hazard and Seismic Risk
Europe Earthquake Reports
General Overview
Earthquake Reports
Middle East Earthquake Reports
General Overview
Earthquake Reports
Social Media
https://t.co/EuFXnSeqlc#TurkeyEarthquake #elazığdepremi #elazığdadeprem #Elazıg #Malatya #Ergani pic.twitter.com/OXwkvfjDJP
🔗 https://t.co/uJIGwWyfBc — #earthquake #grading in #Turkey#Copernicus #CEMS #RapidMapping #EUCivPro
References:
Basic & General References
Specific References
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: Greece
Then, the M 6.8 mainshock hit while i was out and about, followed by a M = 5.2 aftershock.
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.
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.
There are records of tsunami observed on tide gage data.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.
The poster below includes earthquakes that represent the different plate boundaries and tectonic regimes.
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
proportional to the amount of movement. The thick black arrows indicate the mean motion vectors of the plates. The polygonal areas on the map (dashed lines) define the approximate borders of the five different structural regions discussed in the text. The borders between structural regions are not straightforward, and wide transitional zones probably exist between them. The inset shows a schematic map with the geodynamic framework in the eastern Mediterranean area (modified from McClusky et al., 2000). DSF—Dead Sea fault; EAF—East Anatolia fault; HT—Hellenic trench; KFZ— Kefallonia fault zone; MRAC—Mediterranean Ridge accretionary complex; NAF—North Anatolia fault; NAT—North Aegean trough.
Geologic Fundamentals
Compressional:
Extensional:
Europe
General Overview
Earthquake Reports
Social Media
References: