Earthquake Report: Italy

Italy continues to shake following the Armatrice Earthquake series in August 2016. Here is my report for that series of earthquakes. Today’s earthquakes occurred along the northern end of the earthquakes that happened a few months ago.

    Today’s series started with a M 5.5, which was a foreshock to a M 6.1 earthquake. Here are the USGS websites for these two earthquakes. I plot the USGS moment tensors for each of these earthquakes on the interpretive poster below.

  • 2016.10.26 M 5.5 Italy
  • 2016.10.26 M 6.1 Italy

This earthquake is north of the region that had an M 6.3 earthquake in 2009 that led to an interesting (putting it nicely) interaction between scientists, public employees/politicians, and the legal system. Basically, several seismologists were sentenced to prison. More on this is found online, for example, here and here.
Below 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. There is a legend for MMI intensities in the upper part of the interpretive poster below.
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. The tectonics of this region has many normal (extensional) faults, which explain the extensional moment tensor. However, I do not know enough of this region to interpret is this is an east or west dipping fault that ruptured (depends upon which side of which basin experience this earthquake; see below).

    I include some inset figures and maps.

  • In the lower right corner I include 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. This PAGER report suggests that there will be quite a bit of damage from this earthquake (and casualties). This earthquake has a high probability of damage to people and their belongings.
  • Above that I show the Seismic Hazard for Europe as prepared by the SHARE Consortium (Giardini et al., 2013).
  • In the upper left corner I include a basic tectonic map of this region (Woudloper, 2009). Maps with local (larger) scale have much more detailed views of the faulting.
  • In the lower left corner is a map showing the faulting mapped in the region surrounding and including Italy (Billi et al., 2006). There is a convergent plate boundary along the eastern part of Italy (part of the Alpide belt, a convergent boundary that extends from the Straits of Gibraltar to Australia). This fault system dips westward and is onshore in the south, but extends offshore into the Adriatic Sea in central-northern Italy. In the central part of Italy is a series of north-northwest striking extensional faults. It is these extensional (normal) faults that are responsible for the damaging seismicity in this region of central Italy. This includes the 1915, 1997, 2009, and 2016 earthquakes.
  • To the right of the Billi et al. (2006) fault map is a plot showing the seismicity from the last year. Today’s earthquakes are plotted as orange circles and the epicenters from August are plotted as gray circles.


  • Bonio et a. (2004) present a remarkable assessment of the seismic hazard in this region based on a 3-D model for seismogenic sources. I present some of their figures below. I include their original figure captions as blockquotes.
  • This map shows a detailed view of the normal faults in the region. Today’s earthquake is in the region shown in box 3, east of the Umbra Valley.

  • Digital elevation model of central Italy with active normal faults of the Umbria-Marche-Abruzzo Apennines and parameters of active stress tensors obtained by inversion of focal mechanisms of background microseismicity (1), aftershock sequences (2, 3, 4, 5) or striated active faults in seismic areas (6); stress data from Brozzetti and Lavecchia (1994), Boncio and Lavecchia (2000a) and Pace et al. (2002a); the stress axes are given as trend (first three numbers) and plunge (last two numbers).

  • This map shows an even more detailed and large scale view of the faults and seismicity in this region. Today’s earthquakes align to the north of Norcia, approximately along the cross section labeled “sec B.” The two cross sections are in the lower right part of the figure, with section B the lowermost cross section. Today’s earthquake may be on the AF2, the C-NFs (Colfiorito-Norcia fault systems), or MVf (Mt. Vettore fault). The AF2 fault is a proposed low angle detachment fault. These types of faults are controversial in that there are arguments about whether they are seismogenic or not. This year’s Pacific Cell Friends of the Pleistocene field trip in Panamint Valley presented research results that attempted to address this question. In Panamint Valley there are faults that have similar configurations as these faults in Italy.

  • Geological cross sections from seismic reflection profiles across the Gubbio, Gualdo T. and Colfiorito seismic areas (from Boncio et al., 1998; Boncio and Lavecchia, 2000b); epi- and hypocentral distribution of back-ground microseismicity recorded in the Umbria-Marche Apennines and rheological profiles (strength envelopes in critical stress difference, σ1–σ3) built for two different thermal contexts (50 and 40 mW/m2 surface heat flow, see Figure 3 for location); the depth of the brittle-plastic transition on rheological profiles is indicated by arrows; the used rheological parameters are indicated: crustal layering is from DSS data; A (empirical material constant), n (stress exponent) and E (activation energy) are creep parameters; ´ε = longitudinal strain rate (calculated by balancing of a regional geologic section; Figure 5 in Boncio et al., 2000); see text for further details.

  • This map shows a smaller scaled view (than the above figure) with focal mechanisms and cross sections (with structural interpretations). Hypocenters are also plotted on these cross sections. Today’s earthquakes are just south of cross section b. (earthquakes happened here in 1997)

  • Epicentres of the major seismic sequences of the last twenty years (Gubbio, 1984; Colfiorito, 1997; Norcia, 1979; Sangro, 1984) plus three small seismic sequences in the L’Aquila area (1992, 1994 and 1996); seismotectonic sections and rheological profiles built according to the local thermal context. The dashed line (sections ‘a’ and ‘b’) represents the AF low-angle extensional detachment; arrows in seismotectonic sections indicate the maximum depth-extent of the activated seismogenic faults as suggested by the best defined aftershock volume; rheological parameters as in Figure 7; in the southern Abruzzo area, creep strengths for geologic and geodetic longitudinal strain rates are compared (geologic strain rate calculated from data of Galadini and Galli, 2000; geodetic strain rate from D’Agostino et al., 2001); seismological data from Amato et al. (1998); Boncio (1998); Boncio et al. (2004); Cattaneo et al. (2000); De Luca et al. (2000); Deschamps et al. (1984); Ekstrom et al. (1998); Haessler et al. (1988); Harvard CMT database at www.seismology.harvard.edu.

  • Here is a poster that shows the seismic hazard for Europe (Giardini et al., 2013).

  • This is the PAGER alert for the M 6.1 earthquake.

  • Here is a more detailed seismic hazard map for Italy (Stucchi et al., 2011). This shows the median PGA (Peak Ground Acceleration, units of g where g = 9.8 m/s2) that has a 10% probability of exceedance (PE) in the next 50 years. This model assumes a Vs30 greater than 800 m/s. Vs is the average seismic velocity in the upper 30 meters. Vs30 is a proxy used for global to regional estimates of seismic hazard.I include their original figure captions as blockquote.

  • The seismic hazard map showing the PGA distribution with 10% probability of exceedance in 50 years, computed on hard ground (VS30 > 800 m=s).

3 thoughts on “Earthquake Report: Italy

Leave a Reply

Your email address will not be published. Required fields are marked *