Earthquake Report: M 7.1 Indonesia

Yesterday on my way home from a Phil & Graham Lesh show, I got a tsunami notification alert from the National Tsunami Warning Center. There was a magnitude M 6.9 earthquake offshore of Indonesia and there was no tsunami threat for the west coast of the US.

I pulled over to investigate and searched the USGS earthquakes page to locate this earthquake. At that time, there were two events spaced closely in time and space. I suspected that there was probably only one earthquake.

Shortly, that turned out to be true (within minutes). There was a M 7.1 earthquake north of the islands of Bali and Lombok, part of Indonesia.

There was a series of earthquakes in this area a few years ago, which came to my mind. However, this M 7.1 was much deeper.

This M 7.1 earthquake was quite deep (over 500 km!). Those earlier events were shallower and appear to have been related to the Flores thrust fault. Read more about these shallower earthquakes here.

This part of the world is geologically dominated by the convergent plate margin between the Australia and Eurasia plates. This convergent plate margin is part of the Alpide belt, a convergent plate margin that spans almost half the globe (from the northern tip of Australia to the western tip of Portugal). The Alpide belt is responsible for building the tallest mountains in the world (the Asian Himalayas and the European Alps).

Here, in Indonesia, the Australia plate dives beneath the Australia plate forming a subduction zone and a deep sea trench (the Java trench). Earthquakes along this megathrust subduction zone fault have generated strong ground shaking, generated tsunami, and triggered landslides in the past.

In this part of the world, the Eurasia plate is subdivided into a sub-plate called the Sunda plate (so one might see maps with either name labeling this plate).

As the Australia plate subducts it starts out dipping shallowly beneath Java, Bali, Lombok, and the other islands.

The oceanic crust has water within it that helps generate melt in the magma that exists above the Australia plate and beneath the Sunda plate. When this magma melts, its density decreases and the magma rises until it erupts forming the volcanoes that comprise these islands.

As the Australia plate subducts further, the angle that it dips down into the Earth gets steeper. During this process the plate bends and gets exposed to higher pressures (and temperatures).

These physical processes change the stresses within and surrounding the plate. These changes in stress can cause earthquakes like the M 7.1.

Even though this earthquake was large in magnitude, it was so deep so that the shaking intensity was smaller.

The shaking intensity is often reported using the Modified Mercalli Intensity (MMI) scale. People (anyone with internet access) can report their observations on the USGS “Did You Feel It?” web page for this earthquake.

These observations are then used to estimate the shaking intensity felt by those people. Reports from this earthquake show that people on these nearby islands felt intensities around MMI 4 to MMI 5 or so.

Below is my interpretive poster for this earthquake

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

    I include some inset figures. Some of the same figures are located in different places on the larger scale map below. In some of these inset figures I place a blue star to locate the M 7.1 earthquake on these figures.

  • In the upper right corner is a map showing historic seismicity and tectonic plate boundaries.
  • Below that map is a low angle oblique view of a cut away of the Earth along the subduction zone in Java, Indonesia from EOS.
  • In the upper left corner are maps that show the seismic hazard and seismic risk for Indonesia. I spend more time explaining this below.
  • To the right of these hazard and risk maps is a map that shows earthquake intensity using the Modified Mercalli Intensity (MMI) Scale.
  • Above the map is a plot that shows the same intensity (both modeled and reported) data as displayed on the map. Note how the intensity gets smaller with distance from the earthquake.
  • In the lower right corner is a cross section showing earthquake hypocenters (3-D locations) from Darman et a. (2012).
  • To the left of the Darman (2012) plot is a cross section of seismicity presented by Hengresh and Whitney (2016).
  • In the lower center I plot USGS seismicity from the past century. I describe this further below.
  • Here is the map with a month’s seismicity plotted.

  • Here is the plot with a century’s seismicity plotted.
  • These data are limited to within the region shown on the map and i highlight the M 7.1 in orange.
  • These are USGS hypocenters and epicenters for earthquakes between 1923 and 2023 for magnitudes M≥5.
  • One may observe the seismicity within the Australia plate as the plate subducts downwards. The top of the crust is above these seismicity trends.
  • If one looks closely, they will notice a horizontal line of earthquakes at about 30km. These are all with a default depth of 33 km. There is also a row of seismicity with a default depth of 10km. These are depths assigned to events prior to a location assigned with greater certainty. It is highly likely that the depths for these events is different than 10 or 33 km.

Other Report Pages

    Some Relevant Discussion and Figures

    • Below is a map showing historic seismicity (Jones et al., 2014). Cross sections B-B’ and C-C’ are shown. The seismicity for the cross sections below are sourced from within each respective rectangle.

    • Here are the seismcity cross sections.

    • Here is the map from McCaffrey and Nabelek (1987). They used seismic reflection profiles, gravity modeling along these profiles, seismicity, and earthquake source mechanism analyses to support their interpretations of the structures in this region.

    • Tectonic and geographic map of the eastern Sunda arc and vicinity. Active volcanoes are represented by triangles, and bathymetric contours are in kilometers. Thrust faults are shown with teeth on the upper plate. The dashed box encloses the study area.

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

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

    • This are the seismicity cross sections from Hangesh and Whitney (2016). These are shown to compare the subduction zone offshore of Java and the collision zone in the Timor region.

    • Comparison of hypocentral profiles across the (a) Java subduction zone and (b) Timor collision zone (paleo-Banda trench). Catalog compiled from multiple reporting agencies listed in Table 1. Events of Mw>4.0 are shown for period 1815 to 2015.

    • Below are the maps and cross sections from Darman et al., 2012.
    • Here is the map in the interpretive poster above.

    • Tectonic map of the Lesser Sunda Islands, showing the main tectonic units, main faults, bathymetry and location of seismic sections discussed in this paper.

    • Here is the seismicity cross section in the interpretive poster above.

    • This plot shows the earthquake localizations on a South-North cross section for the lat -14°/-4° long 114°/124° quadrant corresponding to the Lesser Sunda Islands region. The localizations are extracted from the USGS database and corresponds to magnitude greater than 4.5 in the 1973-2004 time period (shallow earthquakes with undetermined depth have been omitted.

    • Here is their interpretations of seismic data used to interpret the tectonics of the subduction zone and Flores thrust.

    • Six 15 km deep seismic sections acquired by BGR from west to east traversing oceanic crust, deep sea trench, accretionary prism, outer arc high and fore-arc basin, derived from Kirchoff prestack depth migration (PreSDM) with a frequency range of 4-60 Hz. Profile BGR06-313 shows exemplarily a velocity-depth model according to refraction/wide-angle
      seismic tomography on coincident profile P31 (modified after Lüschen et al, 2011).

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

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

    • In addition to the orientation of relative plate motion (that controls seismogenic zone and strain partitioning), the Indo Australia plate varies in crustal age (Lasitha et al., 2006). I include their figure caption below as a blockquote.

    • Tectonic sketch map of the Sumatra–Java trench-arc region in eastern Indian Ocean Benioff Zone configuration. Hatched line with numbers indicates depth to the top of the Benioff Zone (after Newcomb and McCann13). Magnetic anomaly identifications have been considered from Liu et al.14 and Krishna et al.15. Magnitude and direction of the plate motion is obtained from Sieh and Natawidjaja11. O indicates the location of the recent major earthquakes of 26 December 2004, i.e. the devastating tsunamigenic earthquake (Mw = 9.3) and the 28 March 2005 earthquake (Mw = 8.6).

    • Below are the 4 figures from Koulani et al., 2016. First is the plate tectonic map. I include their figure captions in block quote.

    • Seismotectonic setting of the Sunda-Banda arc-continent collision, East Indonesia. Major faults (thick black lines) [Hamilton, 1979]. Topography and bathymetry are from Shuttle Radar Topography Mission ( Focal mechanisms are from the Global Centroid Moment Tensor. Blue mechanisms correspond to earthquakes with Mw>7 (brown transparent ellipses are the corresponding rupture areas for Flores 1992 and Alor 2004 earthquakes), while the green focal mechanism shows the highest magnitude recorded in Sumbawa. Red dots indicate the locations of major historical earthquakes [Musson, 2012].

    • This figure shows their estimates for plate motion relative velocities as derived from GPS data, constrained by the fault geometry in their block modeling.

    • GPS velocities determined in this study with respect to Sunda Block. Uncertainty ellipses represent 95% confidence level. The inset figure corresponds to the area of the dashed rectangle in the map. Light blue arrows show the velocities for East and West Makassar Blocks.

    • This figure shows their estimates of slip rate deficit along all the plate boundary faults in this region.

    • Relative slip vectors across block boundaries, derived from our best fit model. Arrows show motion of the hanging wall (moving block) relative to the footwall (fixed block) with 95% confidence ellipses. The tails of arrows is located within the “moving” block. Black thick lines show well-defined boundaries we use as active faults in our model and dashed lines show less well-defined boundaries (green : free-slipping boundaries and black: fixed locked faults) . Principal axes of the horizontal strain tensor estimated for the SUMB, EMAK, and EJAV are shown in pink. The thick pink arrow shows the relative motion of Australia with respect to Sunda (AUST/SUND). Abbreviations are Sumba Block (SUMB), West Makassar Block (WMAK), East Makassar Block (EMAK), East Java Block (EJAV), and Timor Block (TIMO). The background seismicity is from the International Seismological Centre catalog with magnitudes ≥5.5 and depths <40 km.

    • Here is their figure that shows the slip deficit along the plate boundary faults.

    • Fault slip rate components: (a) fault normal (extension positive) and (b) fault parallel (right-lateral positive).

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

      • The GEM Seismic Risk Map:

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

    Tsunami Hazard

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

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


    Basic & General References

  • Frisch, W., Meschede, M., Blakey, R., 2011. Plate Tectonics, Springer-Verlag, London, 213 pp.
  • Hayes, G., 2018, Slab2 – A Comprehensive Subduction Zone Geometry Model: U.S. Geological Survey data release,
  • 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, , /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.
  • 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, ,
  • 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,
  • 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.
  • 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,
  • 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).
  • Storchak, D. A., D. Di Giacomo, I. Bondár, E. R. Engdahl, J. Harris, W. H. K. Lee, A. Villaseñor, and P. Bormann (2013), Public release of the ISC-GEM global instrumental earthquake catalogue (1900–2009), Seismol. Res. Lett., 84(5), 810–815, doi:10.1785/0220130034.
  • 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,
  • Specific References

  • Audley-Charles, M.G., 1986. Rates of Neogene and Quaternary tectonic movements in the Southern Banda Arc based on micropalaeontology in: Journal of the Geological Society, London, Vol. 143, 1986, pp. 161-175.
  • Audley-Charles, M.G., 2011. Tectonic post-collision processes in Timor, Hall, R., Cottam, M. A. &Wilson, M. E. J. (eds) The SE Asian Gateway: History and Tectonics of the Australia–Asia Collision. Geological Society, London, Special Publications, 355, 241–266.
  • Baldwin, S.L., Fitzgerald, P.G., and Webb, L.E., 2012. Tectonics of the New Guinea Region in Annu. Rev. Earth Planet. Sci., v. 41, p. 485-520.
  • Benz, H.M., Herman, Matthew, Tarr, A.C., Hayes, G.P., Furlong, K.P., Villaseñor, Antonio, Dart, R.L., and Rhea, Susan, 2011. Seismicity of the Earth 1900–2010 New Guinea and vicinity: U.S. Geological Survey Open-File Report 2010–1083-H, scale 1:8,000,000.
  • Darman, H., 2012. Seismic Expression of Tectonic Features in the Lesser Sunda Islands, Indonesia in Berita Sedimentologi, Indonesian Journal of Sedimentary Geology, no. 25, po. 16-25.
  • Hall, R., 2011. Australia-SE Asia collision: plate tectonics and crustal flow in Geological Society, London, Special Publications 2011; v. 355; p. 75-109 doi: 10.1144/SP355.5
  • Hangesh, J. and Whitney, B., 2014. Quaternary Reactivation of Australia’s Western Passive Margin: Inception of a New Plate Boundary? in: 5th International INQUA Meeting on Paleoseismology, Active Tectonics and Archeoseismology (PATA), 21-27 September 2014, Busan, Korea, 4 pp.
  • Hayes, G.P., Wald, D.J., and Johnson, R.L., 2012. Slab1.0: A three-dimensional model of global subduction zone geometries in, J. Geophys. Res., 117, B01302, doi:10.1029/2011JB008524
  • Jones, E.S., Hayes, G.P., Bernardino, Melissa, Dannemann, F.K., Furlong, K.P., Benz, H.M., and Villaseñor, Antonio, 2014. Seismicity of the Earth 1900–2012 Java and vicinity: U.S. Geological Survey Open-File Report 2010–1083-N, 1 sheet, scale 1:5,000,000,
  • Koulali, A., S. Susilo, S. McClusky, I. Meilano, P. Cummins, P. Tregoning, G. Lister, J. Efendi, and M. A. Syafi’i, 2016. Crustal strain partitioning and the associated earthquake hazard in the eastern Sunda-Banda Arc in Geophys. Res. Lett., 43, 1943–1949, doi:10.1002/2016GL067941
  • Krabbenhoeft, A., Weinrebe, R.W., Kopp, H., Flueh, E.R., Ladage, S., Papenberg, C., Planert, L., and Djajadihardja, Y., 2010. Bathymetry of the Indonesian Sunda margin-relating morphological features of the upper plate slopes to the location and extent of the seismogenic zone in NHESS, v. 10, p. 1899-1911, doi:10.5194/nhess-10-1899-2010
  • McCaffrey, R., and Nabelek, J.L., 1984. The geometry of back arc thrusting along the Eastern Sunda Arc, Indonesia: Constraints from earthquake and gravity data in JGR, Atm., vol., 925, no. B1, p. 441-4620, DOI: 10.1029/JB089iB07p06171
  • Okal, E. A., & Reymond, D., 2003. The mechanism of great Banda Sea earthquake of 1 February 1938: applying the method of preliminary determination of focal mechanism to a historical event in EPSL, v. 216, p. 1-15.
  • Silver, E.A., Breen, N.A., and Prastyo, H., 1986. Multibeam Study of the Flores Backarc Thrust Belt, Indonesia, in JGR., vol. 91, no. B3, p. 3489-3500
  • Zahirovic, S., Seton, M., and Müller, R.D., 2014. The Cretaceous and Cenozoic tectonic evolution of Southeast Asia in Solid Earth, v. 5, p. 227-273, doi:10.5194/se-5-227-2014

Return to the Earthquake Reports page.