This is an executable version of a paper published in Geophysical Research Letters. This notebooks contains the text of the version submitted for review, and does not incorporate the minor revisions made to the text as per the review, which do not affect the computation.
The paper may be found here.
Low-angle normal faults are well-described in the rock record and may serve an important role in crustal extension. However, a significant earthquake on a continental low-angle normal fault has not been observed, and such slip is often interpreted to be in conflict with standard rock mechanical theory. The lack of observed earthquakes with focal mechanisms clearly indicating low-angle normal slip may be an indication that they are not seismically active, or it may be due to the fact that these earthquakes are infrequent compared to the length of focal mechanism catalogs. To address this, we create a compilation of all potentially active continental low-angle normal faults and calculate the likelihood of observing a significant earthquake on them over time windows from 1 to 100 years. We find 20 candidate faults in extensional zones worldwide. We find that the probability of observing a significant low-angle normal fault earthquake is dependent on several factors including the frequency-magnitude distribution, but for either a characteristic or Gutenberg-Richter distribution we calculate a probability of about 0.5 that an earthquake greater than $M6.5$ (and therefore likely to have a known fault scarp and dip angle) will be observed on any low-angle normal fault in a time window of 35 years, which is the length of the Global CMT catalog. We then use Bayes' rule to illustrate how the absence of observed significant low-angle normal fault seismicity over the catalog period moderately decreases the likelihood that the structures generate large earthquakes, but does not reduce the likelihood to zero.
Low-angle normal faults (LANFs), with dips less than 30$^\circ$, are well described in the geologic record. They are thought to play an important role in accommodating large-magnitude continental extension [Howard and John, 1987] and crustal thinning [Lister et al., 1986], and their recognition has been a major development in continental tectonics [Wernicke, 2009]. However, despite widespread field observations of inactive LANFs and their central role in extensional tectonic theory, they remain enigmatic and contentious structures, and it is not clear if they are seismically active at low dip angles in the upper crust. This is for two reasons: because brittle faulting on LANFs is in apparent conflict with standard Andersonian rock mechanical theory as typically applied to the upper crust [Axen, 2004], and because observations of active faulting on LANFs are sparse and at times ambiguous [Wernicke, 1995]. A considerable amount of research has been performed to address the former concern, reconciling LANF slip with rock mechanics [e.g., Axen and Bartley, 1997; Collettini, 2011]. The latter issue is highlighted by studies that have searched the focal mechanism catalogs and found no normal faulting earthquakes with focal mechanisms and surface ruptures clearly indicating slip on planes $\le30^\circ$ [Jackson, 1987; Collettini and Sibson, 2001], which is taken as conclusive evidence that LANFs are inactive or aseismic. However, the lack of observed seismic slip on continental LANFs may be simply be because they are rare structures with long recurrence intervals, so earthquakes on them are very infrequent. Without knowing the likelihood of observing a LANF rupture in a time window of a few decades, it is not clear if an empty search result is strong evidence against LANF seismicity. If this likelihood is known, though, Bayesian probability theory provides a framework for quantifying how the negative search results impact the probability that LANFs are seismogenic.
In this work, we estimate the maximum likelihood of a significant LANF event occurring in time windows from 1 to 100 years, and then we interpret the lack of observed LANF seismicity in a quantified, probabilistic context using Bayesian methods. We estimate the maximum observation likelihood by treating all potentially active LANFs described in the literature as seismically active at their surface dip angles throughout the upper crust. Under these assumptions, we create synthetic earthquake catalogs with both Gutenberg-Richter and `characteristic' frequency--magnitude distributions, using each fault's geometry and slip rate. We then calculate the probability of observing earthquakes on at least one LANF over different observation periods. Then, we use Bayes' rule to incorporate the negative catalog search results and the observance likelihood to show how the negative results reduce the probability that LANFs are seismically active, but do not bring the final probability to zero.
Areas of the crust undergoing active extension are generally assumed to have a subvertical maximum compressive stress. Mohr-Coulomb theory, as applied to the crust, predicts that a fault with a typical coefficient of friction for rocks (0.6--0.8) should lock up if it is oriented at an angle greater than 60$^\circ$ to the maximum compressive stress (i.e., fault dips less than 30$^\circ$), and new, optimally oriented faults should form [Sibson, 1985]. Therefore, for normal faults with dips less than 30$^\circ$, either much lower fault friction or elevated pore fluid pressure is required for fault slip.
Evidence for seismic slip on LANFs is sparse. This is partly due to the ambiguity of the rupture plane in earthquake focal mechanisms, and thus a focal mechanism with a low angle nodal plane will also have a high angle nodal plane. Without ancillary information indicating which nodal plane corresponds to the slip surface, searches of earthquake catalogs cannot yield unique results as to whether they contain LANF events. Several collections of normal fault earthquakes with known surface breaks [Jackson, 1987; Collettini and Sibson, 2001], thereby resolving dip ambiguity, contain no low-angle events, although we note the total number of events in these collections are small ($\le$ 25 events). Some candidate LANF events exist, but they are undersea [e.g., Abers, 2001] or difficult to verify [e.g., Doser, 1987].
We have compiled all potentially active LANFs with known subareal fault traces from a thorough review of the literature; there are twenty total. We have then mapped the approximate fault traces into a GIS file (available at https://github.com/cossatot/LANF_gis), with metadata such as slip rate and source. We then have estimated the probability of observing an earthquake above a given magnitude for each fault individually over some time window, and then calculated the probability of observing a significant earthquake on any of the faults over that same time window.
Over the past decade or so, many field studies have found evidence for LANF activity in orogens throughout the world. These studies typically find arrays of Quaternary normal fault scarps on the fault traces and/or in the hanging walls of mapped or inferred low-angle detachment faults [e.g., Axen et al., 1999] . Some studies also have bedrock thermochronology data from the exhumed detachment footwalls that are suggestive of ongoing rapid exhumation [e.g., Sundell et al. 2013], although this data does not preclude a recent cessation of faulting. In some cases, additional evidence for LANF activity comes from geophysical data such as GPS geodesy [e.g., Hreinsdóttir and Bennett, 2009] and seismic waves [e.g., Doser, 1987].
from IPython.core.display import Image Image(filename='active_lanfs_map_insets.png')