Earthquake Report: M 7.0 Vanuatu

Early this morning (my time) I got a notification from the Pacific Tsunami Warning Center that there was no tsunami threat from an M 7.2 earthquake in the Vanuatu Islands.

Later, as I woke up I checked the USGS website to see that there was an M7.0 earthquake offshore of the Vanuatu Islands.

Based on the depth of the hypocenter (the 3-D location of the earthquake) it appears that this M 7.0 ruptured a thrust fault within the Australia plate. Given the uncertainty of the location of the megathrust fault, it is possible that this actually was on the megathrust subduction zone fault (so is what we call an “interface” event). I don’t think that the USGS finite fault model is correct (it seems unlikely that this earthquake ruptured a fault within the Australia plate and slipped up into the upper plate). But I could be wrong (which is quite common).

I don’t always have the time to write a proper Earthquake Report. However, I prepare interpretive posters for these events.

Because of this, I present Earthquake Report Lite. (but it is more than just water, like the adult beverage that claims otherwise). I will try to describe the figures included in the poster, but sometimes I will simply post the poster here. (UPDATE: I could not resist spending a little time looking at updated papers from this region, so have included some figures below.)

Below is my interpretive poster for this earthquake

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

    I include some inset figures.

  • Here is the map with 3 month’s seismicity plotted.

Some Supporting Information

  • Here is the USGS poster showing the seismicity for this region from 1900-2010 (Benz et al., 2011). Below I include the legend (not the correct scale; click on this link for the entire poster (65 MB pdf)). Note the cross section F-F’ which I plot on the poster above.

  • Here is the cross section F-F’ again, with the legend below.

  • Here are some figures from Bergeot et al., 2009.
  • This first figure shows the tectonic setting and the plate convergence rates (the rates, in mm/year, that the Australia plate is converging relative to the North Fiji Basin.

  • (a) Geodynamic setting of the VSZ, with block motions relative to the North Fiji Basin [from Calmant et al., 2003]. The Vanuatu arc is split into three blocks, with anticlockwise rotation (north), convergence (center), and clockwise rotation (south). Dashed line is the BATB; solid lines are the spreading ridge; bold line is the VSZ. Bathymetry data are from Calmant et al. [2002]. The black rectangle is the central part of the Vanuatu arc. White arrows are velocities (millimeters per year) with respect to the Australian plate (AP); black arrows are block motion with respect to the North Fiji Basin. Dotted line is the cross section of Figure 2b. (b) Schematic of the central part of the VSZ [from Lagabrielle et al., 2003]. The direction of this cross section is west to east, and it intersects the Santo and Maewo Islands (dotted line in Figure 2a). Abbreviations are as follows: IAB, Aoba Intra-arc Basin; BATB, back-arc thrust belt; NFB, North Fiji Basin.

  • This figure shows the horizontal motion rates (in mm/year) for the GPS sites in the region.

  • Horizontal interseismic GPS velocities for the VSZ in an Australia-fixed reference frame. The Australian motion is estimated as a rigid rotation from our GPS results with a least squares inversion. Abbreviations are as follows: WTP, West Torres Plateau; DER, D’Entrecasteaux Ridge. Lines are (1) BATB, (2) spreading ridge, (3) VSZ, (4) discontinuity supposed between TGOA and Epi island, and (5) transition zone.

  • This figure shows a map where they plot, in the next figure, a comparison between their modeled vertical velocities with the observed vertical velocities.

  • Transects and GPS stations used to assess the locked zone parameters in this study. Shaded triangles represent the A-A0 (TNMR, LVMP, LMBU, WLRN, SWBY, VMVS, NSUP, RNSR, and AMBR) transect GPS stations, and solid triangles represent the B-B0 transect GPS stations (LISB, TASM, AVNA, RATA, RATU, SANC, AOBA, PNCT, and MAWO). The bold lines represent the A-A0 and B-B0 transects. The white arrows show the convergence direction. Abbreviations are as follows: DER, D’Entrecasteaux Ridge; WTP, Wet Torres Plateau. The stars indicate the edge of the locked zone as deduced from the GPS velocity interpretation (Figure 12). Lines are (1) BATB and (2) VSZ.

  • This figure shows the comparison between their modeled velocities and the observed velocities.

  • (top) Vertical and (bottom) horizontal (perpendicular to the trench) velocity profiles for the GPS stations of the A-A0 (open circle) and B-B0 (filled circle) transects. Distances are given with respect to the trench. The bold curves represent the best fit of the locked zone and long-term convergence rate model (dip, 20; width, 50 km; slip, 54 mm a1) estimated from observed velocities. Lines 2 and 3 represent the effect of the width variation in the model (45 and 60 km, respectively). See Figure 11 for the transect location.

  • Here are the two figures from Cleveland et al. (2014).
  • Figure 1. I include the figure caption below as a blockquote.

  • (left) Seismicity of the northern Vanuatu subduction zone, displaying all USGS-NEIC earthquake hypocenters since 1973. The Australian plate subducts beneath the Pacific in nearly trench-orthogonal convergence along the Vanuatu subduction zone. The largest events are displayed with dotted outlines of the magnitude-scaled circle. Convergence rates are calculated using the MORVEL model for Australia Plate relative to Pacific Plate [DeMets et al., 2010]. (right) All GCMT moment tensor solutions and centroids for Mw ≥ 5 since 1976, scaled with moment. This region experiences abundant moderate and large earthquakes but lacks any events with Mw >8 since at least 1900.

  • Figure 17. I include the figure caption below as a blockquote.

  • One hundred day aftershock distributions of all earthquakes listed in the ISC catalog for the 1966 sequence and in the USGS-NEIC catalog for the 1980, 1997, 2009, and 2013 sequences in northern Vanuatu. The 1966 main shocks are plotted at locations listed by Tajima et al. [1990]. Events of the 1997 and 2009 sequences were relocated using the double difference method [Waldhauser and Ellsworth, 2000] for P wave first arrivals based on EDR picks. The event symbol areas are scaled relative to the earthquake magnitudes based on a method developed by Utsu and Seki [1954]. Hypocenters of most aftershock events occurred at <50 km depth.

  • Figure 17. I include the figure caption below as a blockquote.

  • (right) Space-time plot of shallow (≤ 70 km) seismicity M ≥ 5.0 in northern Vanuatu recorded in the NEIC catalog as a function of distance south of ~10°N, 165.25°E. (left) The location of the seismicity on a map rotated to orient the trench vertically.

  • Craig et al. (2014) evaluated the historic record of seismicity for subduction zones globally. In particular, the evaluated the relations between upper and lower plate stresses and earthquake types (cogent for the southern New Hebrides trench). Below is a figure from their paper for this part of the world. I include their figure caption below in blockquote.

  • Outer-rise seismicity along the New Hebrides arc. (a) Seismicity and focal mechanisms. Seismicity at the southern end of the arc is dominated by two major outer-rise normal faulting events, and MW 7.6 on 1995 May 16 and an MW 7.1 on 2004 January 3. Earthquakes are included from Chapple & Forsyth (1979); Chinn & Isacks (1983); Liu & McNally (1993). (b) Time versus latitude plot.

  • Here is a summary figure from Craig et al. (2014) that shows different stress configurations possibly existing along subduction zones.

  • Schematic diagram for the factors influencing the depth of the transition from horizontal extension to horizontal compression beneath the outer rise. Slab pull, the interaction of the descending slab with the 660 km discontinuity (or increasing drag from the surround mantle), and variations in the interface stress influence both the bending moment and the in-plane stress. Increases in the angle of slab dip increases the dominance of the bending moment relative to the in-plane stress, and hence moves the depth of transition towards the middle of the mechanical plate from either an shallower or a deeper position. A decrease in slab dip enhances the influence of the in-plane stress, and hence moves the transition further from the middle of the mechanical plate, either deeper for an extensional in-plane stress, or shallower for a compressional in-plane stress. Increased plate age of the incoming plate leads to increases in the magnitude of ridge push and intraplate thermal contraction, increasing the in-plane compressional stress in the plate prior to bending. Dynamic topography of the oceanic plate seawards of the trench can result in either in-plane extension or compression prior to the application of the bending stresses.

  • Here are the figures from Richards et al. (2011) with their figure captions below in blockquote.
  • Here is the map showing the current configuration of the slabs in the region.

  • Map showing distribution of slab segments beneath the Tonga-Vanuatu region. West-dipping Pacifi c slab is shown in gray; northeast-dipping Australian slab is shown in red. Three detached segments of Australian slab lie below the North Fiji Basin (NFB). HFZ—Hunter Fracture Zone. Contour interval is 100 km. Detached segments of Australian plate form sub-horizontal sheets located at ~600 km depth. White dashed line shows outline of the subducted slab fragments when reconstructed from 660 km depth to the surface. When all subducted components are brought to the surface, the geometry closely approximates that of the North Fiji Basin.

  • This is the cross section showing the megathrust fault configuration based on seismic tomography and seismicity.

  • Previous interpretation of combined P-wave tomography and seismicity from van der Hilst (1995). Earthquake hypocenters are shown in blue. The previous interpretation of slab structure is contained within the black dashed lines. Solid red lines mark the surface of the Pacifi c slab (1), the still attached subducting Australian slab (2a), and the detached segment of the Australian plate (2b). UM—upper mantle;
    TZ—transition zone; LM—lower mantle.

  • Here is their time step interpretation of the slabs that resulted in the second figure above.

  • Simplifi ed plate tectonic reconstruction showing the progressive geometric evolution of the Vanuatu and Tonga subduction systems in plan view and in cross section. Initiation of the Vanuatu subduction system begins by 10 Ma. Initial detachment of the basal part of the Australian slab begins at ca. 5–4 Ma and then sinking and collision between the detached segment and the Pacifi c slab occur by 3–4 Ma. Initial opening of the Lau backarc also occurred at this time. Between 3 Ma and the present, both slabs have been sinking progressively to their current position. VT—Vitiaz trench; dER—d’Entrecasteaux Ridge.

  • Here are two great figures from Deng et al. (2022). This article focuses on the influence of the D’Entrecasteaux Ridge on subduction here in Vanuatu. They focus on geochemical data from magmatic rocks in the area.
  • This one shows a variety of processes going on in this area.

  • Geological and geophysical constraints regarding the evolution of the Vanuatu arc. (a) Bathymetric map showing the locations of islands for which samples were included in our geochemical compilation. Slab dip contours below the Vanuatu arc are displayed every 20° (from Hayes et al., 2018). (b) Bathymetric map of the Vanuatu arc and an inset showing depth-to-slab versus distance-from-trench for each of the sample localities included in our compilation (Table S1 in Supporting Information S1). Slab depth contours beneath the Vanuatu arc are displayed every 20 km (from Hayes et al., 2018). The orange lines show the chosen cross sections (i.e., Sections A, B, C) across the different blocks of the Vanuatu arc, which were used to estimate slab dips. Orange dots denote the location of Deep Sea Drilling Project Site 286 and Ocean Drilling Program Legs 134 Sites 828 and 831. (c) Interpreted geodynamic setting of the Vanuatu arc based on modern global positioning system velocity measurements (observed, black arrows; modeled, white arrows; from Bergeot et al., 2009). The Vanuatu arc can be divided into three tectonic blocks that are separated by two strike-slip faults (magenta dashed lines; Calmant et al., 2003; Taylor et al., 1995), which are the counterclockwise rotated Northern Block, the eastward migrated Central Block and the clockwise rotated Southern Block. Orange arrows indicate plate convergence velocities (in mm/year) with respect to the Australian plate (Bergeot et al., 2009). (d) Intermediate-depth seismicity distribution (50–170 km) since 1972 with magnitudes in the range of 4–7, from USGS Earthquake Catalog ( The seismic gap is highlighted by a solid polygon. The wide red arrow depicts the influx of hot sub-slab mantle to the forearc mantle wedge through a slab tear.

  • Schematic of the Vanuatu subduction zone to illustrate the model proposed by this study. The conceptual model highlights the role that the subducting buoyant D’Entrecasteaux ridge plays in the dynamic evolution of the Vanuatu arc. The introduction of D’Entrecasteauz Ridge causes shallow subduction and the development of a slab tear south of the ridge and the segmentation of the Vanuatu arc into the Northern Block, Central Block and Southern Block. Shallow slab subduction beneath the Central Block results in (a) squeezing out of the asthenospheric mantle; (b) scraping off the bottom of the ancient continental lithospheric mantle beneath the forearc, which then migrates ahead of the advancing slab and forms a bulldozed keel underneath the main-arc and (c) transmitting compressional stresses in the over-riding plate, which inhibits the formation of backarc spreading and instead produces a backarc thrust belt. Additionally, the ingression of hot subslab mantle causes partial melting of the cold forearc mantle and produces magmatism anomalously close to the trench (i.e., the Efate, Nguna, and Pele volcanoes that are situated in the forearc).

  • Here is a video that shows a simulation from Deng et al. (2022).
    • Baillard et al. (2015) provide some insight into the geometry of the Australia plate slab. I am adjusting my hypothesis to be that this M 7.0 probably was on the megathrust, based on their work.
    • This first figure shows a map and some cross sections.
    • Note how cross section 4 is really close to the M 7.0. This geometry places the top of the slab below the M 7.0 hypocenter.

    • Geometry of the subduction interface and updip/downdip extents of the seismogenic zone. (a) Map view. The green contour is at 27 km depth and marks the intersection with the fore-arc Moho. The dashed contours present the updip and downdip extents of the seismogenic zone. The numbered lines showthe location of cross sections plotted to the right. NDR: north d’Entrecasteaux ridge; BS: Bougainville seamount. (b) Geometric cross sections of the subduction interface (depth as a function of distance from the subduction front). (c) Dip cross sections (dip angle as a function of distance from the subduction front).

    • This figure shows how they interpret the seismicity in this area.

    • Cross section of seismic activity through the center of our total catalog (only events with residuals <0.2 s are plotted). Three clusters of activity are observed: (1) around the subduction interface (green), (2) within the subducting plate beneath the subduction interface (red), and (3) at intermediate depths (blue). The dotted line is our interpretation of the subduction interface.

    • This figure shows some specific earthquakes they used in their analyses. Earthquake 1 (pink) is really close to the M 7.0. They interpret that event to be in the upper part of the Australia plate. At this point, I suggest that it is equivocal, about whether or not the M 7.0 was an interface event or a slab or crustal event.

    • Clusters and focal mechanisms in the local catalog. Simple focal mechanisms are illustrated in black, composite focal mechanisms in colors corresponding to the cluster events (circles). P axes indicated in red. (a) Map view. The boxes indicate the orientation and dimensions of the cross sections. (b) Cross section beneath Santo Island. (c) Cross section between Santo and Malekula Islands. The cross sections also show the picked subduction interface (thick black curve), the Australian Plate Moho (dotted line, assuming a 8 km thick crust), and the North Fiji Basin Moho (dotted line, assuming a 27 km thick fore-arc crust).


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

  • Baillard, C., W. C. Crawford, V. Ballu, M. Régnier, B. Pelletier, and E. Garaebiti (2015), Seismicity and shallow slab geometry in the central Vanuatu subduction zone, J. Geophys. Res. Solid Earth,120,5606–5623,
  • Benz, H.M., Herman, Matthew, Tarr, A.C., Furlong, K.P., Hayes, G.P., Villaseñor, Antonio, Dart, R.L., and Rhea, Susan, 2011. Seismicity of the Earth 1900–2010 eastern margin of the Australia plate: U.S. Geological Survey Open-File Report 2010–1083-I, scale 1:8,000,000.
  • Bergeot, N., M. N. Bouin, M. Diament, B. Pelletier, M. Re´gnier, S. Calmant, and V. Ballu (2009), Horizontal and vertical interseismic velocity fields in the Vanuatu subduction zone from GPS measurements: Evidence for a central Vanuatu locked zone, J. Geophys. Res., 114, B06405,
  • Cleveland, K.M., Ammon, C.J., and Lay, T., 2014. Large earthquake processes in the northern Vanuatu subduction zone in Journal of Geophysical Research: Solid Earth, v. 119, p. 8866-8883, doi:10.1002/2014JB011289.
  • Deng, C., Jenner, F. E., Wan, B., & Li, J.-L. (2022). The influence of ridge subduction on the geochemistry of Vanuatu arc magmas. Journal of Geophysical Research: Solid Earth, 127, e2021JB022833.
  • Hayes, G. P., D. J. Wald, and R. L. Johnson, 2012. Slab1.0: A three-dimensional model of global subduction zone geometries, J. Geophys. Res., 117, B01302, doi:10.1029/2011JB008524.
  • Richards, S., Holm., R., Barber, G., 2011. When slabs collide: A tectonic assessment of deep earthquakes in the Tonga-Vanuatu region, Geology, v. 39, pp. 787-790.
  • Oceanic-Oceanic Subduction Zone Figure
    Music Reference (in 1900-2016 seismicity video)

  • Bumba Crossing Kevin MacLeod ( | Licensed under Creative Commons: By Attribution 3.0 License |

Return to the Earthquake Reports page.