Earthquake Report: Lombok, Indonesia

Earlier today there was a shallow M 6.4 earthquake with an epicenter on the island of Lombok, Indonesia. With a hypocentral depth of about 7.5 km, this size of an earthquake can be quite damaging. The USGS PAGER estimate of impact suggests that there is about a 10% chance that there are more than 10 fatalities. Hopefully there are none. There have been several aftershocks, two M > 5.

This earthquake is probably along a thrust fault associated with the Flores thrust fault, a north vergent (dipping into the earth in a southerly direction) back thrust fault to the Sunda subduction zone fault. The Flores thrust possibly extends from east of Timor on the east to the northern shore of Java (McCaffrey and Nabelek, 1987). Others suggest that the Flores thrust ends at a cross fault just east of Lombok (Hengresh and Whitney, 2016). However, the seismic profiles from Silver et al. (1986) are convincing that there are east-west compressional structures extending between the northern shore of Java to where the Flores thrust is mapped.

Detailed mapping of the seafloor to the east of Lombok, north of the island of Sumbawa, reveals that there are imbricate (overlapping) thrust faults (Silver et al., 1986). I think that it is reasonable to presume that there are similar structures on the northern flank of Lombok.

Lombok is also a volcano complex as part of the Sunda magmatic arc. There may be fault systems associated with the volcanic activity. I include tectonic faults that are included in the global scale fault data set from the Coordinating Committee for Geoscience Programme in East and Southeast Asia. The most active volcano on Lombok is the Rinjani volcano. Here is a great place to learn about this volcano (the Volcano Discovery website).

If the M 6.4 earthquake was on the Flores fault, it would need to dip at about 10°. The Flores thrust fault proposed by Hengesh and Whitney (2016) has a much steeper dip. So this sequence is probably in the upper plate somewhere.

There was a M 6.0 earthquake to the east of the M 6.4, but it was much deeper (almost 600 km), so is unlikely to be genetically related to the M 6.4 sequence.

Magnetic Anomalies

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

Historic Seismicity

  • Below I discuss analogues to today’s M 6.4 earthquake.
  • To the west, between Lombok and Bali, there was a series of earthquakes all in 1979. They happened several months apart, but had a similar magnitude and orientation. The hypocentral depths were in the 25-40 km depth range, so some of these may have been on the Flores thrust system. These alone suggest that the Flores thrust extends at least this far west.
  • To the east, along the eastern part of Sumbawa, there was a series of earthquakes in the first decade of the 21st century, from 2002-2009. These also all share a similar magnitude range and orientation. These earthquakes all happened within a narrow range of depths (18-20 km; though the 2002 earthquake has a default depth on 10 km).
  • Based on earthquakes in the regions to the east and to the west, it is possible that this M 6.4 is the first of a series of mid M 6 earthquakes (either within a year like in Bali or over several years like Sumbawa).

Below is my interpretive poster for this earthquake

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I include earthquake epicenters from 1918-2018 with magnitudes M ≥ 6.0.
I plot the USGS fault plane solutions (moment tensors in blue and focal mechanisms in orange), possibly in addition to some relevant historic earthquakes.

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

  • In the upper right corner is a low angle oblique view of the Sunda subduction zone beneath Java, Bali, Lombok, and Sumbawa (from Earth Observatory Singapore). I place a blue star in the general location of today’s earthquake’s epicenter (as for all figures here). The India-Australia plate is subducting northwards beneath the Sunda plate (part of the Eurasia plate).
  • In the upper left corner is a plate tectonic map showing the major fault systems, volcanic arc islands, and oceanic plateaus and basins of the region (Darman, 2012). The map shows the Flores thrust extending as far west as Lombok. Compare the complicated tectonics in the eastern portion of this region compared to the western portion of this region.
  • To the right of the Darman (2012) map is a cross section of seismicity presented by Hengresh and Whitney (2016). These authors argue for a north vergent Flores thrust in this region, though most of their work was on the subduction/collision zone.
  • In the lower right corner is another, earlier, tectonic map from Silver et al. (1986). These authors use seismic reflection and multibeam bathymetry data to map the Flores thrust as far as Java, west of Bali. The location for the map in the lower left corner of this interpretive poster is outlined here as a dashed line rectangle.
  • In the lower left corner is a map from Silver et al. (1986) that shows the detailed mapping of imbricate north (and some south) vergent thrust faults.
  • Here is the same map but with seismicity from the past month.


  • Here is the same map but with historic seismicity.


USGS Earthquake Pages

    These are from this current sequence

  • 2018-07-28 17:07:23 UTC M 6.0
  • 2018-07-28 22:47:37 UTC M 6.4
  • 2018-07-28 23:06:49 UTC M 5.4
  • 2018-07-29 01:50:32 UTC M 5.3

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.

  • Here is a map of the same general area from Silver et al. (1986), used here to locate the following large scale map.

  • Location of SeaMARC II survey (Plate 1 and Figures 2) and geographic features discussed in text. Triangles on upper plates of thrust zones.

  • This is the large scale map showing the detailed thrust fault mapping (Silver et al., 1986).

  • Bathymetry, faults, and mud diapirs of the central Flores thrust zone, based on interpretation of SeaMARC II data and seismic reflection profiles. Shown also are locations (circled numbers) of all seismic profiles. Mud diapirs are solid black. Triangles on upper plates of thrust faults.

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

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

  • Here are some focal mechanisms from earthquakes in the region from Hangesh and Whitney (2016). Symbol color represents depth.

  • (a) Focal mechanism solutions for the study region. The focal mechanisms are classified based on depth intervals to illustrate the style of faulting within the different structural domains. Note (b) sinistral reverse motion along Timor trough, (c) subduction related pattern along Java trench, and dextral solutions along the western Australia extended margin (Figure 4a) north of 20°S. Centroid moment tensor (CMT) solutions [Dziewonski et al., 1981] are from the CMT project [Ekström et al., 2012; http://www.globalcmt.org/CMTcite.html] for events of Mw>5.0 for the period 1976 onward.

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

  • This map from Hangesh and Whitney (2016) shows the GPS velocities in this region. Note the termination of the Flores thrust and the north-northeast striking (oriented) cross fault between Lombok and Sumbawa.

  • GPS velocities of Sunda and Banda arc region. Large black and grey arrow shows motion of Australia relative to Eurasia [DeMets et al., 1994]. Thin black arrows show GPS velocities of Sunda and Banda arc regions relative to Australia [Nugroho et al., 2009]. Seismicity from ISC-GEM catalog [Storchak et al., 2013]. Note reduction of station velocities from west to east indicating progressive coupling of the Banda arc to the Australian plate compared to the area along the Sunda arc.

Geologic Fundamentals

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

  • There are three types of earthquakes, strike-slip, compressional (reverse or thrust, depending upon the dip of the fault), and extensional (normal). Here is are some animations of these three types of earthquake faults. The following three animations are from IRIS.
  • Strike Slip:

    Compressional:

    Extensional:

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

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

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

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

    References:

  • Audley-Charles, M.G., 1986. Rates of Neogene and Quaternary tectonic movements in the Southern Banda Arc based on micropalaeontology in: Journal of fhe Geological Society, London, Vol. 143, 1986, pp. 161-175.
  • Audley-Charles, M.G., 2011. Tectonic post-collision processes in Timor, Hall, R., Cottam, M. A. &Wilson, M. E. J. (eds) The SE Asian Gateway: History and Tectonics of the Australia–Asia Collision. Geological Society, London, Special Publications, 355, 241–266.
  • Baldwin, S.L., Fitzgerald, P.G., and Webb, L.E., 2012. Tectonics of the New Guinea Region in Annu. Rev. Earth Planet. Sci., v. 41, p. 485-520.
  • Benz, H.M., Herman, Matthew, Tarr, A.C., Hayes, G.P., Furlong, K.P., Villaseñor, Antonio, Dart, R.L., and Rhea, Susan, 2011. Seismicity of the Earth 1900–2010 New Guinea and vicinity: U.S. Geological Survey Open-File Report 2010–1083-H, scale 1:8,000,000.
  • 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, https://dx.doi.org/10.3133/ofr20101083N.
  • 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

Earthquake Report: Banda Sea

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

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

Below is my interpretive poster for this earthquake.

I plot the seismicity from the past month, with color representing depth and diameter representing magnitude (see legend). I also include seismicity from 1917-2017 for earthquakes with magnitudes M ≥ 7.0. I show the fault plane solutions for some of these earthquakes.

  • moment tensors for 2005, 2012, and 2017 (USGS)
  • focal mechanisms for 1987 (USGS) and 1938 (Okal and Reymond, 2003)
  • The fault plane solutions for the 1963, 1987, 2005, and 2012 earthquakes are all very similar to the 2017 M 5.5. However, these earthquakes are form two depth “populations.” The 1963 and 1987 earthquakes are at ~65 km depth, while the 2005, 2012, and 2017 are between 150-200 km. There are some earthquakes that are much shallower depth eastward along this possibly strike-slip fault. Near where this fault comes on land at the base of the Bird’s Head in New Guinea, there are some earthquakes from the past couple of decades that also have strike-slip fault–plane solutions. The earthquake depths along this Sorong fault (Hall, 2011) appear to show that the Sorong fault is active beneath the Sunda plate. The 1938 earthquake may be the result of some form of strain partitioned faulting(?). Alternatively, these earthquakes may be unrelated to the Sorong fault. There may be some internal structure in the Australia plate that is interacting with the subduction zones or other faults (some preexisting structure that is optimally oriented to reactivate with the .
  • I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely
  • I also include the shaking intensity contours on the map. These use the Modified Mercalli Intensity Scale (MMI; see the legend on the map). This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.
  • I include the slab contours plotted (Hayes et al., 2012), which are contours that represent the depth to the subduction zone fault. These are mostly based upon seismicity. The depths of the earthquakes have considerable error and do not all occur along the subduction zone faults, so these slab contours are simply the best estimate for the location of the fault. The hypocentral depth plots this close to the location of the fault as mapped by Hayes et al. (2012).

    I include some inset figures in the poster.

  • In the upper left corner is a general tectonic map for this part of the world. I placed a green star in the location of this M 5.5 earthquake (Zahirovic et al., 2014).
  • In the upper right corner is a low-angle oblique view of the plate boundaries in the northern part of this region (Hall, 2011). The upper part of the diagram shows the opposing vergent subduction zones along that strike north-south along the Molucca Strait (Halmahera, Philippines). The lower panel shows the downgoing Australia plate along the Timor Trench and Seram Trench. Note the location of the Bird’s Head, the northwestern part of New Guinea. I have also labeled this region in the main map for comparison. The strike-slip fault at the northern boundary of New Guinea is the Sorong fault and this is labeled in this Hall (2011) figure.
  • Below the Hall (2011) figure is a figure from Baldwin et al. (2012) that shows the regional seismicity and faulting as they are related to different geologic types in New Guinea. I placed a green star in the location of this M 5.5 earthquake.
  • In the lower right corner is a part of the USGS Poster that reviews the seismicity of this region for the past century or so (Benz et a., 2011). The map shows seismicity with depth, along with some cross section locations. I place a green star at the location of this M 5.5 earthquake. I present three of the cross sections from this poster, A-A’, B-B’, and C-C.’ Of particular interest is the section B-B’ because this is placed near the M 5.5 earthquake. I have placed a green star that represents the hypocentral location on cross section B-B.’ The hypocentral depth suggests this M 5.5 earthquake is in the downgoing Australia plate slab.
  • In the lower left corner is a diagram showing the subducting Australia plate at the Java Trench (Yves Descatoire).
  • To the right of the Java Trench figure presents a detailed view of the faulting along the eastern Java and western Timor trenches (Hangesh and Whitney, 2016). They present evidence for oblique motion along the Timor trough. And present evidence for a backthrust on the northern side of Timor and the Indonesia islands east of Java.


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

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

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

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

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

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

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

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

References:

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

Earthquake Report: Australia!

Yesterday there was an earthquake in northeastern Australia. Here is the USGS website for this M 5.7 earthquake.

This earthquake occurred in a region that also experienced a similar magnitude earthquake in 2011. Here is the USGS website for that M 5.0 earthquake. More can be found about the larger earthquakes that have occurred in Australia here.

Here is my interpretive map that shows the epicenter, along with the shaking intensity contours. These contours use the Modified Mercalli Intensity (MMI) scale. The MMI is a qualitative measure of shaking intensity. More on the MMI scale can be found here and here. This is based upon a computer model estimate of ground motions, different from the “Did You Feel It?” estimate of ground motions that is actually based on real observations.

I placed a moment tensor / focal mechanism legend on the poster. There is more material from the USGS web sites about moment tensors and focal mechanisms (the beach ball symbols). Both moment tensors and focal mechanisms are solutions to seismologic data that reveal two possible interpretations for fault orientation and sense of motion. One must use other information, like the regional tectonics, to interpret which of the two possibilities is more likely. There are some faults in the region that are northeast striking, so this slightly favors an interpretation of a right-lateral (dextral) strike-slip earthquake.

    I include some inset figures and maps.

  • In the upper right corner is a figure from Matthews et al. (2011) showing the felt region from the 2011.04.16 earthquake near Bowen, AU.
  • To the left of that is a figure that shows the geology, faults, and seismicity in Australia. This map is from Australian Government Geoscience Australia (AGSO).
  • In the lower right corner is a map of the seismic hazard in Australia from 1991 (so it is dated). There is an updated version of this map here, with an online interface here. However, this region of AU still has a low seismic hazard.
  • In the lower left corner is the Rapid Assessment of an Earthquake’s Impact (PAGER) report. More on the PAGER program can be found here. An explanation of a PAGER report can be found here. PAGER reports are modeled estimates of damage. On the top is a histogram showing estimated casualties and on the right is an estimate of possible economic losses.
  • In the upper left corner is a plot of seismicity from the Queensland University Advanced Centre for Earthquake Studies (QUAKES).


Here is more on the 2011.04.16 M 5.0 Earthquake. First the figure and then below is a description from AGSO in blockquote.


On Saturday, 16 April 2011, a magnitude 5.3 earthquake occurred at 3:31pm local time, located 50km west of Bowen in central Queensland, near Mount Abbot. The event was widely felt along the Queensland coast, including Cairns about 400km from the epicentre, and further west in Hughenden about 370km away. Local residents experienced significant ground-shaking, with a maximum intensity of MMI V experienced in Ravenswood and Bowen with reports of slight damage in Guthalungra and Bowen. A focal mechanism produced for the main shock indicates that movement was predominantly strike slip. Four temporary seismic stations were installed around the epicentre, recording over 300 small aftershocks in the six weeks following the main shock. Five aftershocks ranging in magnitude from ML 3.2 to 4.1 were recorded on 16, 17 and 19 April. The Bowen earthquake was the largest recorded in the region since a magnitude 5.7 earthquake occurred north of Ravenswood (about 80km west of Mount Abbot) in December 1913