Macroseismology

Macroseismology: the study of felt earthquakes

Nowadays seismology mainly focuses on digital data acquisition gathered by seismometers installed in a seismic network and the detection of earthquakes. Earthquakes, however, happened long before we were able to measure them in the instrumental era. By gathering and mapping witness’ reports on felt earthquakes, the impact of an earthquake can be reliably mapped and the magnitude and depth of historical earthquakes can be estimated based on the distribution of the felt reports. Usually the German mathematician P.N.G. Egen is credited for the development of the first modern intensity scale based on the perception of 28 February 1828 Ms ~4.6 earthquake near Jauche in Belgium. Later pioneering work of Nicholas Ambraseys in the 1970ies showed the enormous value of studying historical earthquakes.

Since 1932, after the ML 4.7 Uden earthquake, when a felt earthquake occurs, the Royal Observatory of Belgium sends out national earthquake inquiries to each Belgian municipality to construct an intensity map of the felt area. Based on the responses, an earthquake intensity value is assigned to each municipality resulting in an Communal Intensity Map. The first success of macroseismology in Belgium was shown by Somville (1939) who gathered the felt reports of the largest instrumentally recorded earthquake on the Belgian territory: i.e. the Zulzeke-Nukerke (near Oudenaarde) 11 June 1938 Ms 5.0 earthquake.

Since the arrival of the internet, the paper inquiry has been complemented by an online internet "Did You Feel It?" (DYFI) macroseismic inquiry. Online DYFI inquiries allow rapid and automatic construction of community internet intensity macroseismic maps with which the authorities and the population can be quickly informed about the damage distribution and felt area of an earthquake, either officially or through social media. This information is complementary to the recordings by the seismic stations to study earthquakes and their impact.

Europe: a problem in gathering macroseismic data?

In the US, the “Did You Feel It?” (DYFI?) system of the USGS asks questions about nationally and internationally felt earthquakes. The advantage of DYFI? is that macroseismic data are collected systematically over the whole North American continent and worldwide by ONE SINGLE INSTITUTE (USGS) through a well-calibrated algorithm (Dengler and Dewey, 1998; Worden et al., 2000).

In Europe, however, the situation is more complex: at least 34 seismological institutes in 24 countries and 2 international institutes (EMSC and USGS) manage and maintain their own web-based macroseismic questionnaire (see table). Germany wins the award of having the most inquiries (8 !) in one country. The benefit of having different European inquiries is that every institute has made the inquiry available in its own national language(s) and hence can deal with specifics in their language. However, disadvantages are that

(1) people can answer multiple national or international inquiries when they felt an earthquake

(2) the data on the perception of internationally-felt seismic events is strongly fragmented across different institutes and countries.

Merging databases carries the risk of duplicate entries in the database but, more importantly, also has the risk of merging intensities that result from different questionnaires, countries and intensity calculation algorithms, which will smooth the mean of the merged intensities. Because of this European fragmentation and intensity variability, performing a proper macroseismic assessment of transfrontier-felt earthquakes remains a complex and not straightforward task.

The table and figure give an overview of European seismological institutes that provide an online macroseismic inquiry to request responses regarding earthquakes that have been felt. Hyperlinks to the online questionnaires are provided in the URL column. QRT? (quasi real time?) asks whether the gathered data are mapped and illustrated online in QRT: y(es), n(o). Source = Van Noten et al. 2017.

A solution?

The current challenge in Europe is to generate denser and more accurate shaking maps by merging all available intensity datasets from national and international seismological institutes. For example, the ROB cooperates with the Cologne University on a trans-frontier online macroseismic data acquisition. Small- to moderate-magnitude (M4 to M5) earthquakes affecting the Belgian and neighbouring territories cause (several thousand) respondents to fill this online questionnaire. To apply this methodology on an European scale, a consensus on a common questionnaire and data gathering is required. Something that will take time and effort.

To provide a simple solution how data can be shared, I have investigated the possibility to store intensity data from different institutes in equally-sized grid cells. First a python software was developed to geocode individual addresses into their geographical coordinates. Afterwards, instead of using classic DYFI representations (average internet intensities per community), I assign average intensities to single 100 km2 grid cells. The resulting macroseismic grid cell distribution of the affected area shows a less subjective and more homogeneous intensity distribution than the classical irregular community. This method was extensively explored in my open access 2017 Solid Earth paper.

Averaging intensities in grid cells is however not ideal. All earthquake effect need to be taken into account for a proper intensity assessment. Determining the intensity of a community, area or grid cell thus requires a more complex approach than just averaging. HOWEVER, because different countries have different questionnaires and are not transparent in their intensity algorithm, averaging intensities is the currently the best we can do. Before questionnaires are homogenised in Europe, this grid cell methodology provides a simple yet rough solution.

LEFT: Ramsgate 2015 ML 4.2 earthquake intensity map. Merged grid cell intensity map generated with responses from ROB-BNS, BCSF, BGS and intensity data from EMSC and USGS. The 2015 Ramsgate example clearly shows the influence of the Brabant Massif on good earthquake wave propagation along the structure of this massif.

RIGHT: Goch 2011 ML 4.3 earthquake intensity map. Merged grid cell intensity map generated with geocoded responses from ROB-BNS, KNMI and intensity data from EMSC and USGS. Note the absence of macroseismic entries in NE Belgium. Grid cell size is 100 km2. Rings represent 50 km, 100 km and 150 km epicentral distance. The 2011 Goch example demonstrates that, apart from effective wave propagation, bedrock depth is also important to determine if an earthquake can be felt. E.g. in the NE part of Belgium, the Goch earthquake was not felt because of high-frequency filtering by the sedimentary cover.

Source: Van Noten et al. 2017.

Acknowledgements

Maps created in QGIS.

Want to read more on felt events in Belgium? Please visit this webpage http://seismologie.be/en/research/seismology/macroseismology

Conference Presentations related to Macroseismology

  1. Van Noten, K., de Rubeis, V., Tosi, P., Sbarra, P. 2018. Challenges of the new ESC Working Group Macroseismology – Integrating corrective parameters to merge multiple-sourced online macroseismic data. 6th Colloquium on historical earthquakes & paleoseismology studies, Han-sur-lesse, Belgium. 24-25 October 2018
  2. Tosi, P., De Rubeis, V., Van Noten, K. & Sbarra, P. 2018. Integration of corrective parameters to merge multiple-sourced online macroseismic data. 36th General Assembly of the European Seismological Commission. Malta. 2-7 September 2018.
  3. Van Noten, K. & Lecocq, T. 2017. Macroseismology in Europe. Technical Workshop on internet macroseismology, Ljubljana, Slovenia, 14-15 November 2017. Poster Presentation
  4. Van Noten, K., Lecocq, T., Hinzen, K.-G., Sira, C. & Camelbeeck, T. 2017. Merging transfrontier internet macroseismic data of two M4 earthquakes in NW Europe using a grid cell approach. Technical Workshop on internet macroseismology, Ljubljana, Slovenia, 14-15 November 2017. Presentation
  5. Van Noten, K., Lecocq, T., Hinzen, K.-G., Sira, C. & Camelbeeck, T. 2016. Transfrontier Macroseismic Data Exchange in Europe: Intensity Assessment of M4 Earthquakes by a Grid Cell Approach. American Geophysical Union, AGU2016, 12-16 December 2016, San Francisco, USA.
  6. Consentino, B.C., Van Noten, K., Sintubin, M. 2016. Macroseismic intensity assessment of the 2008-2010 Walloon Brabant Seismic Swarm (Belgium) by a griell procedure. 48th Congresso Brazileiro de Geologia, 9-13 October 2016, Porto Alegre, Brazil.
  7. Van Noten, K., Lecocq, T., Hinzen, K.-G., Sira, C. & Camelbeeck, T. 2016. Transfrontier Macroseismology in Europe. ESC2016, 4-9 September 2016, Trieste, Italy.
  8. Van Noten, K., Lecocq, T., Hinzen, K.-G., Sira, C. & Camelbeeck, T. 2016. Transfrontier macroseismic data exchange in NW Europe: examples of non-circular intensity distributions. Geophysical Research Abstracts 18, EGU2016-4609, session SM3.3, 17-22 April 2016, Vienna, Austria.
  9. Van Noten, K., Lecocq, T. & Camelbeeck, T. 2016. Macroseismic analysis of "Did You Feel It?" responses to ML > 3 earthquakes felt in and around Belgium. 5th International Geologica Belgica Conference, January 2016, Mons, Belgium.
  10. Van Noten, K., Lecocq, T., Hinzen, K. & Camelbeeck, T. 2015. Macroseismic analysis of "Did You Feel It?" responses to ML > 3 earthquakes felt in Belgium and its neighbouring countries. Membach 20-40 Workshop: Earthquake activity and hazard in NW Europe. La Gileppe, Belgium. 15 October 2015.
  11. Lecocq, T., Hinzen, K., Van Noten, K. & Camelbeeck, T. 2015. Nearly 15 years of real-time transfrontier macroseismic data acquisition and mapping. Membach 20-40 Workshop: Earthquake activity and hazard in NW Europe. La Gileppe, Belgium. 15 October 2015.
  12. Van Noten, K., Lecocq, T., Vleminckx, B., Camelbeeck, T., 2013. Macroseismic investigation of the 2008-2010 low magnitude seismic swarm in the Brabant Massif, Belgium. The link between macroseismic intensity and geomorphology. Geophysical Research Abstracts 15, EGU2013-4923, session SM2.8, 7-12 april 2013, Vienna.