Geology ReportsSearch

SEARCH · Geology Reports

Results for “The Seismological Record”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14Linked to original sources

Orientation-independent, nongeometric-mean measures of seismic intensity from two horizontal components of motion

New measures of spectral intensity based on the horizontal components of ground shaking are introduced. These new measures are independent of the in situ orientation of the recordings and encompass the full range of spectral amplitudes over all possible rotation angles. Unlike previously introduced measures that are also orientation independent, no geometric means are used in the computation of the new measures. The new measures based on fiftieth percentile values of the response spectra show small but systematic increases (to a factor of about 1.07 at a 10 sec period) compared to the comparable geometric-mean measure.

Bulletin of the Seismological Society of America

Application of an iterative least-squares waveform inversion of strong-motion and teleseismic records to the 1978 Tabas, Iran, earthquake

An iterative least-squares technique is used to simultaneously invert the strong-motion records and teleseismic P waveforms for the 1978 Tabas, Iran, earthquake to deduce the rupture history. The effects of using different data sets and different parametrizations of the problem (linear versus nonlinear) are considered. A consensus of all the inversion runs indicates a complex, multiple source for the Tabas earthquake, with four main source regions over a fault length of 90 km and an average rupture velocity of 2.5 km/sec.

Bulletin of the Seismological Society of America

The MW 7.0 Haiti Earthquake of January 12, 2010: USGS/EERI Advance Reconnaissance Team Report

Executive Summary A field reconnaissance in Haiti by a five-member team with expertise in seismology and earthquake engineering has revealed a number of factors that led to catastrophic losses of life and property during the January 12, 2010, Mw 7.0 earthquake. The field study was conducted from January 26 to February 3, 2010, and included investigations in Port-au-Prince and the heavily damaged communities to the west, including Leogane, Grand Goave, Petite Goave, and Oliver. Seismology Despite recent seismic quiescence, Haiti has suffered similar devastating earthquakes in the historical past (1701, 1751, 1770 and 1860). Despite this knowledge of historical seismicity, Haiti had no seismograph stations during the main earthquake, so it is impossible to estimate accurately the intensity of ground motions. Nonetheless, the wide range of buildings damaged by the January 12, 2010, earthquake suggests that the ground motions contained seismic energy over a wide range of frequencies. Another earthquake of similar magnitude could strike at any time on the eastern end of the Enriquillo Fault, directly to the south of Port-au-Prince. Reconstruction must take this hazard into account. The four portable seismographs installed by the team recorded a series of small aftershocks. As expected, the ground motions recorded at a hard-rock site contained a greater proportion of high frequencies than the motions recorded at a soil site. Two of the stations continue to monitor seismic activity. A thorough field investigation of the mapped Enriquillo Fault south of the city of Leogane failed to find any evidence of surface faulting. This led the team to conclude that the earthquake was unlikely to have produced any surface rupture in the study area. Geotechnical Aspects Soil liquefaction, landslides and rockslides in cut slopes, and road embankment failures contributed to extensive damage in Port-au-Prince and elsewhere. A lack of detailed knowledge of the physical conditions of the soils (for example, lithology, stiffness, density, and thickness) made it difficult for us to quantitatively assess the role of ground-motion amplification in the widespread damage. Buildings The Haitian Ministry of Statistics and Informatics reported that one-story buildings represent 73 percent of the building inventory. Most ordinary, one-story houses have roofs made of sheet metal (82 percent), whereas most multistory houses and apartments have roofs made of concrete (71 percent). Walls made of concrete/block/stone predominate both in ordinary houses and apartments. It appears that the widespread damage to residences and commercial and government buildings was attributable to a great extent to the lack of earthquake-resistant design. In many cases, the structural types, member dimensions, and detailing practices were inadequate to resist strong ground motions. These vulnerabilities may have been exacerbated by poor construction practices. Reinforced concrete frames with concrete block masonry infill appeared to perform particularly poorly. Structures with light (timber or sheet metal) roofs performed better compared to structures with concrete roofs and slabs. The seismic performance of some buildings was adequate, and some of the damaged buildings appeared to have had low deformation demands. These observations suggest that structures designed and constructed with adequate stiffness and reinforcing details would have resisted the earthquake without being damaged severely. A damage survey of 107 buildings in downtown Port-au-Prince indicated that 28 percent had collapsed and another 33 percent were damaged enough to require repairs. A similar survey of 52 buildings in Leogane found that 62 percent had collapsed and another 31 percent required repairs. Bridges There was no evidence of bridge collapses attributable to the earthquake. Most bridges in Port-au-Prince are simple box culverts consisting of box girders 2.0 to 2.

Open-File Report

Automated detection of clipping in broadband earthquake records

Because the amount of available ground‐motion data has increased over the last decades, the need for automated processing algorithms has also increased. One difficulty with automated processing is to screen clipped records. Clipping occurs when the ground‐motion amplitude exceeds the dynamic range of the linear response of the instrument. Clipped records in which the amplitude exceeds the dynamic range are relatively easy to identify visually yet challenging for automated algorithms. In this article, we seek to identify a reliable and fully automated clipping detection algorithm tailored to near‐real‐time earthquake response needs. We consider multiple alternative algorithms, including (1) an algorithm based on the percentage difference in adjacent data points, (2) the standard deviation of the data within a moving window, (3) the shape of the histogram of the recorded amplitudes, (4) the second derivative of the data, and (5) the amplitude of the data. To quantitatively compare these algorithms, we construct development and holdout datasets from earthquakes across a range of geographic regions, tectonic environments, and instrument types. We manually classify each record for the presence of clipping and use the classified records. We then develop an artificial neural network model that combines all the individual algorithms. Testing on the holdout dataset, the standard deviation and histogram approaches are the most accurate individual algorithms, with an overall accuracy of about 93%. The combined artificial neural network method yields an overall accuracy of 95%, and the choice of classification threshold can balance precision and recall.

Seismological Research Letters

Documentation for a web site to serve ULF-EM (Ultra-Low Frequency Electromagnetic) data to the public

The Stanford Ultra-Low Frequency Electromagnetic (ULF-EM) Monitoring Project is recording naturally varying electromagnetic signals adjacent to active earthquake faults, in an attempt to establish whether there is any variation in these signals associated with earthquakes. Our project is collaborative between Stanford University, the U.S. Geological Survey (USGS), and UC Berkeley. Lead scientists are Simon Klemperer (Stanford University), Jonathan Glen (USGS) and Darcy Karakelian McPhee (USGS). Our initial sites are in the San Francisco Bay Area, monitoring different strands of the San Andreas fault system, at Stanford University's Jasper Ridge Biological Preserve (JRSC), Marin Headlands of the Golden Gate National Recreation Area (MHDL), and the UC Berkeley's Russell Reservation Field Station adjacent to Briones Regional Park (BRIB). In addition, we maintain in conjunction with the Berkeley Seismological Laboratory (BSL) two remote reference stations at the Bear Valley Ranch in Parkfield, Calif., (PKD) and the San Andreas Geophysical Observatory at Hollister, Calif., (SAO). Metadata about our site can be found at http://ulfem-data.stanford.edu/info.html. Site descriptions can be found at the BSL at http://seismo.berkeley.edu/, and seismic data can be obtained from the Northern California Earthquake Data Center at http://www.ncedc.org/. The site http://ulfem-data.stanford.edu/ allows access to data from the Stanford-USGS sites JRSC, MHDL and BRIB, as well as UC Berkeley sites PKD and SAO.

Open-File Report

Regional spectral analysis of three moderate earthquakes in Northeastern North America

We analyze Fourier spectra obtained from the horizontal components of broadband and accelerogram data from the 1997 Cap-Rouge, the 2002 Ausable Forks, and the 2005 Rivière-du-Loup earthquakes, recorded by Canadian and American stations sited on rock at hypocentral distances from 23 to 602 km. We check the recorded spectra closely for anomalies that might result from site resonance or source effects. We use Beresnev and Atkinson’s (1997) near-surface velocity structures and Boore and Joyner’s (1997) quarter-wave method to estimate site response at hard- and soft-rock sites. We revise the Street et al. (1975) model for geometrical spreading, adopting a crossover distance of r o =50 km instead of 100 km. We obtain an average attenuation of Q =410±25f 0.50±0.03 for S+L g +surface waves with ray paths in the Appalachian and southeastern Grenville Provinces. We correct the recorded spectra for attenuation and site response to estimate source spectral shape and radiated energy for these three earthquakes and the 1988 M 5.8 Saguenay earthquake. The Brune stress drops range from 130 to 419 bars, and the apparent stresses range from 39 to 63 bars. The corrected source spectral shapes of these earthquakes are somewhat variable for frequencies from 0.2 to 2 Hz, falling slightly below the fitted Brune spectra.

Bulletin of the Seismological Society of America

The 30 November 2018 Mw7.1 Anchorage Earthquake

The M w "> M w 7.1 47 km deep earthquake that occurred on 30 November 2018 had deep societal impacts across southcentral Alaska and exhibited phenomena of broad scientific interest. We document observations that point to future directions of research and hazard mitigation. The rupture mechanism, aftershocks, and deformation of the mainshock are consistent with extension inside the Pacific plate near the down‐dip limit of flat‐slab subduction. Peak ground motions > 25 % g "> g >25%g were observed across more than 8000    km 2 "> 8000 km 2 ⁠ , though the most violent near‐fault shaking was avoided because the hypocenter was nearly 50 km below the surface. The ground motions show substantial variation, highlighting the influence of regional geology and near‐surface soil conditions. Aftershock activity was vigorous with roughly 300 felt events in the first six months, including two dozen aftershocks exceeding M 4.5. Broad subsidence of up to 5 cm across the region is consistent with the rupture mechanism. The passage of seismic waves and possibly the coseismic subsidence mobilized ground waters, resulting in temporary increases in stream flow. Although there were many failures of natural slopes and soils, the shaking was insufficient to reactivate many of the failures observed during the 1964 M 9.2 earthquake. This is explained by the much shorter duration of shaking as well as the lower amplitude long‐period motions in 2018. The majority of observed soil failures were in anthropogenically placed fill soils. Structural damage is attributed to both the failure of these emplaced soils as well as to the ground motion, which shows some spatial correlation to damage. However, the paucity of instrumental ground‐motion recordings outside of downtown Anchorage makes these comparisons challenging. The earthquake demonstrated the challenge of issuing tsunami warnings in complex coastal geographies and highlights the need for a targeted tsunami hazard evaluation of the region. The event also demonstrates the challenge of estimating the probabilistic hazard posed by intraslab earthquakes.

Alaska

A seismic nodal deployment to understand magmatic structure in the vicinity of the Pahala earthquake swarm

In summer-fall 2022, 80 three-component SmartSolo IGU-BD3C-5 nodal seismometers were deployed surrounding the Pāhala seismic swarm on the Island of Hawaiʻi, with the goal of improving seismicity catalogs, and seismic velocity images of the crust and upper mantle in this region. The Pāhala swarm, located south of Mauna Loa and Kīlauea, has been the site of a multiyear sustained swarm of seismicity at depths of ~ 25 – 40 km, with order of magnitude increases in rate in 2015, and then again in 2019. This seismicity is possibly related to the input of magma from the mantle plume below, which may then be subsequently transported to volcanic edifices. However, these processes remain enigmatic, in part due to a lack of precise earthquake locations and seismic velocity models in this region. Here we provide an overview of the deployment, an assessment of the quality of the collected data, and discuss the viability of the dataset for local earthquake relocation, tomography, and teleseismic receiver functions. Through comparisons with proximal permanent broadband and short period instruments, we find that the nodes produce high quality data, particularly at periods shorter than 5 s, although we find, document, and correct discrepancies with the gain and polarities of the instruments. We successfully record signals from teleseismic earthquakes, even at periods longer than 5 s (the corner of the flat response of the nodes). We also record local earthquakes, including details related to source characteristics. This indicates that the data is likely to prove useful for investigations using both local and teleseismic earthquake signals to better understand the connections between the deep and shallow magmatic systems of Hawaiʻi. While this deployment provides a snapshot in time, its success may provide a useful benchmark for future studies as the volcanic systems of Hawaiʻi continue to evolve in the future.

Hawaii

GEOS seismograms recorded for aftershocks of the earthquakes of December 7, 1988, near Spitak, Armenia SSR, during the time period 3 January 1989 through 2 February 1988 (UTC)

The earthquakes of December 7, 1988, near Spitak, Armenia SSR, serve as another grim reminder of the serious hazard that earthquakes pose throughout the world. We extend our heartfelt sympathies to the families of the earthquake victims and intend that our cooperative scientific endeavours will help reduce losses in future earthquakes. Only through a better understanding of earthquake hazards can earthquake losses be reduced for all peoples in seismically active regions of the world. The tragic consequences of these earthquakes remind scientists and public officials alike of their urgent responsibilities to understand and mitigate the effects of earthquakes. On behalf of the U.S. Geological Survey, I would like to express appreciation to our Soviet colleagues for their kind invitation to participate in joint scientific and engineering studies. Without their cooperation and generous assistance, the conduct of these studies would not have been possible. This report provides seismologic and geologic data collected during the time period December 21, 1988, through February 2, 1989. These data are presented in their entirety to expedite analysis of the data set for inferences regarding hazard mitigation actions, applicable not only in Armenia but other regions of the world exposed to high seismic risk

Open-File Report

Empirical ground-motion relations for subduction-zone earthquakes and their application to Cascadia and other regions

Ground-motion relations for earthquakes that occur in subduction zones are an important input to seismic-hazard analyses in many parts of the world. In the Cascadia region (Washington, Oregon, northern California, and British Columbia), for example, there is a significant hazard from megathrust earthquakes along the subduction interface and from large events within the subducting slab. These hazards are in addition to the hazard from shallow earthquakes in the overlying crust. We have compiled a response spectra database from thousands of strong-motion recordings from events of moment magnitude ( M ) 5-8.3 occurring in subduction zones around the world, including both interface and in-slab events. The 2001 M 6.8 Nisqually and 1999 M 5.9 Satsop earthquakes are included in the database, as are many records from subduction zones in Japan (Kyoshin-Net data), Mexico (Guerrero data), and Central America. The size of the database is four times larger than that available for previous empirical regressions to determine ground-motion relations for subduction-zone earthquakes. The large dataset enables improved determination of attenuation parameters and magnitude scaling, for both interface and in-slab events. Soil response parameters are also better determined by the data.

British Columbia, California, Oregon, Washington

Determination of earthquake magnitude for early warning from the time-dependence of P-wave amplitudes

We propose a method that utilizes the time dependence of P ‐wave displacement amplitudes to estimate the final magnitude ( ⁠M⁠ ) for earthquake early warning (EEW) before the arrival of the peak amplitude. A relation between M and P ‐wave displacement amplitude is employed for the method. Its value is set as a function of time from the P arrival, and is determined using a K‐NET dataset recorded in Japan from a scaling relation between M and the time dependence of P ‐wave displacement. A test to check the performance of the proposed equation demonstrates in a statistical sense that this technique enables us to estimate M more rapidly than conventional methods without loss of accuracy. We conclude that the approach proposed in this article effectively gains a longer lead time as well as reduces the blind zone for EEW.

Bulletin of the Seismological Society of America

GEOS seismograms recorded for aftershocks of the earthquakes of December 7, 1988, near Spitak, Armenia SSR, during the time period 26 December 1988 14:00 through 29 December 1988 (UTC)

The earthquakes of December 7, 1988, near Spitak, Armenia SSR, serve as another grim reminder of the serious hazard that earthquakes pose throughout the world. We extend our heartfelt sympathies to the families of the earthquake victims and intend that our cooperative scientific endeavours will help reduce losses in future earthquakes. Only through a better understanding of earthquake hazards can earthquake losses be reduced for all peoples in seismically active regions of the world. The tragic consequences of these earthquakes remind scientists and public officials alike of their urgent responsibilities to understand and mitigate the effects of earthquakes. On behalf of the U.S. Geological Survey, I would like to express appreciation to our Soviet colleagues for their kind invitation to participate in joint scientific and engineering studies. Without their cooperation and generous assistance, the conduct of these studies would not have been possible. This report provides seismologic and geologic data collected during the time period December 21, 1988, through February 2, 1989. These data are presented in their entirety to expedite analysis of the data set for inferences regarding hazard mitigation actions, applicable not only in Armenia but other regions of the world exposed to high seismic risk.

Open-File Report

Ground motions recorded in Rome during the April 2009 L’Aquila seismic sequence: site response and comparison with ground‐motion predictions based on a global dataset

The mainshock and moderate‐magnitude aftershocks of the 6 April 2009 M 6.3 L’Aquila seismic sequence, about 90 km northeast of Rome, provided the first earthquake ground‐motion recordings in the urban area of Rome. Before those recordings were obtained, the assessments of the seismic hazard in Rome were based on intensity observations and theoretical considerations. The L’Aquila recordings offer an unprecedented opportunity to calibrate the city response to central Apennine earthquakes—earthquakes that have been responsible for the largest damage to Rome in historical times. Using the data recorded in Rome in April 2009, we show that (1) published theoretical predictions of a 1 s resonance in the Tiber valley are confirmed by observations showing a significant amplitude increase in response spectra at that period, (2) the empirical soil‐transfer functions inferred from spectral ratios are satisfactorily fit through 1D models using the available geological, geophysical, and laboratory data, but local variability can be large for individual events, (3) response spectra for the motions recorded in Rome from the L’Aquila earthquakes are significantly amplified in the radial component at periods near 1 s, even at a firm site on volcanic rocks, and (4) short‐period response spectra are smaller than expected when compared to ground‐motion predictions from equations based on a global dataset, whereas the observed response spectra are higher than expected for periods near 1 s.

Rome

Testing small-aperture array analysis on well-located earthquakes, and application to the location of deep tremor

We have here analyzed local and regional earthquakes using array techniques with the double aim of quantifying the errors associated with the estimation of propagation parameters of seismic signals and testing the suitability of a probabilistic location method for the analysis of nonimpulsive signals. We have applied the zero-lag cross-correlation method to earthquakes recorded by three dense arrays in Puget Sound and Vancouver Island to estimate the slowness and back azimuth of direct P waves and S waves. The results are compared with the slowness and back azimuth computed from the source location obtained by the analysis of data recorded by the Pacific Northwest seismic network (PNSN). This comparison has allowed a quantification of the errors associated with the estimation of slowness and back azimuth obtained through the analysis of array data. The statistical analysis gives ??BP = 10?? and ??BS = 8?? as standard deviations for the back azimuth and ??SP = 0.021 sec/km and ??SS = 0.033 sec /km for the slowness results of the P and S phases, respectively. These values are consistent with the theoretical relationship between slowness and back azimuth and their uncertainties. We have tested a probabilistic source location method on the local earthquakes based on the use of the slowness estimated for two or three arrays without taking into account travel-time information. Then we applied the probabilistic method to the deep, nonvolcanic tremor recorded by the arrays during July 2004. The results of the tremor location using the probabilistic method are in good agreement with those obtained by other techniques. The wide depth range, of between 10 and 70 km, and the source migration with time are evident in our results. The method is useful for locating the source of signals characterized by the absence of pickable seismic phases.

Bulletin of the Seismological Society of America

Rupture history of the 2011 M 9 Tohoku Japan earthquake determined from strong‐motion and high‐rate GPS recordings: Subevents radiating energy in different frequency bands

Strong‐motion records from KiK‐net and K‐NET, along with 1 sample/s Global Positioning System (GPS) records from GEONET, were analyzed to determine the location, timing, and slip of subevents of the M 9 2011 Tohoku earthquake. Timing of arrivals on stations along the coast shows that the first subevent was located closer to the coast than subevent (2), which produced the largest slip. A waveform inversion of data from 0 to 0.2 Hz indicates that the first subevent primarily ruptured down‐dip and north of the hypocenter and had an M of 8.5. The areas of this subevent that generated the low (<0.2 Hz) and high (>0.2 Hz) frequency energy are located in the same vicinity. The inversion result for the second subevent ( M 9.0) has large slip on the shallow part of the fault with peak slip of about 65 m above about 25 km depth. This slip generated the tsunami. The preferred inversion has initiation of subevent 2 on the shallow portion of the fault so that rupture proceeded down‐dip and mainly to the south. Subevent 2 started about 35 s after subevent 1, which allows for the possibility of dynamic triggering from subevent 1. The slip model predicts displacements comparable to those found from ocean‐bottom transducers near the epicenter. At frequencies that most affect tall buildings (0.1–0.5 Hz), there is a strong pulse (subevent 3) in the strong‐motion records that arrives after the near‐field ramp from subevent 2. High‐frequency subevent 3 was located down‐dip and south of the high‐slip portion of subevent 2 and was initiated as rupture from subevent 2 proceeded down‐dip. The compact pulse for subevent 3 is modeled with an M 8.0 source in a 75 by 30 km area that ruptured down‐dip and to the south with a high slip velocity, indicating high stress drop.

Bulletin of the Seismological Society of America

Results of seismological monitoring in the Cascade Range 1962-1989: earthquakes, eruptions, avalanches and other curiosities

Modern monitoring of seismic activity at Cascade Range volcanoes began at Longmire on Mount Rainier in 1958. Since then, there has been an expansion of the regional seismic networks in Washington, northern Oregon and northern California. Now, the Cascade Range from Lassen Peak to Mount Shasta in the south and Newberry Volcano to Mount Baker in the north is being monitored for earthquakes as small as magnitude 2.0, and many of the stratovolcanoes are monitored for non-earthquake seismic activity. This monitoring has yielded three major observations. First, tectonic earthquakes are concentrated in two segments of the Cascade Range between Mount Rainier and Mount Hood and between Mount Shasta and Lassen Peak, whereas little seismicity occurs between Mount Hood and Mount Shasta. Second, the volcanic activity and associated phenomena at Mount St. Helens have produced intense and widely varied seismicity. And third, at the northern stratovolcanoes, signals generated by surficial events such as debris flows, icequakes, steam emissions, rockfalls and icefalls are seismically recorded. Such records have been used to alert authorities of dangerous events in progress. -Authors

Geoscience Canada

The uses and limitations of the square‐root‐impedance method for computing site amplification

The square‐root‐impedance (SRI) method is a fast way of computing approximate site amplification that does not depend on the details from velocity models. The SRI method underestimates the peak response of models with large impedance contrasts near their base, but the amplifications for those models is often close to or equal to the root mean square of the theoretical full resonant (FR) response of the higher modes. On the other hand, for velocity models made up of gradients, with no significant impedance changes across small ranges of depth, the SRI method systematically underestimates the theoretical FR response over a wide frequency range. For commonly used gradient models for generic rock sites, the SRI method underestimates the FR response by about 20%&ndash;30%. Notwithstanding the persistent underestimation of amplifications from theoretical FR calculations, however, amplifications from the SRI method may often provide more useful estimates of amplifications than the FR method, because the SRI amplifications are not sensitive to details of the models and will not exhibit the many peaks and valleys characteristic of theoretical full resonant amplifications (jaggedness sometimes not seen in amplifications based on averages of site response from multiple recordings at a given site). The lack of sensitivity to details of the velocity models also makes the SRI method useful in comparing the response of various velocity models, in spite of any systematic underestimation of the response. The quarter‐wavelength average velocity, which is fundamental to the SRI method, is useful by itself in site characterization, and as such, is the fundamental parameter used to characterize the site response in a number of recent ground‐motion prediction equations.

Bulletin of the Seismological Society of America

How does the onset of offset influence geologic slip rates?

Geologic slip rates are typically based on the displacement accrued by a geomorphic or stratigraphic feature and the age of the offset feature. Because slip rates are commonly calculated by dividing the displacement of a faulted marker by its age, they contain two open time intervals: the elapsed time between the age of an offset feature and the age of the earthquake that displaced the feature, and the time between the present‐day and the most recent earthquake. Here, we explore the influence of including unconstrained open intervals in geologic slip rate calculations. We test the degree to which these open intervals affect geologic slip rates and their uncertainties, and we find that their influence depends primarily on mean earthquake recurrence intervals (RIs). Slip rates on faults with longer RIs, such as the Wasatch fault, can be greatly influenced by an increase of up to 20% when accounting for open intervals. In contrast, slip rates on faults with shorter RIs, such as the San Andreas fault, are only slightly influenced by the assumption that slip rates calculated over open intervals approximate those calculated over closed intervals. Our analyses indicate that faults with moderate slip rates (∼0.2–5 mm/yr) are sensitive to both open interval effects themselves, as well as methods to quantify and account for these effects. We re‐evaluate how slip rates are calculated and defined in displacement–time space using published deformation records. We explore the utility of assigning a probability distribution to the initiation of offset of the oldest faulted feature and the timing of the most recent earthquake (MRE). We find that calculating geologic slip rates without using probability distributions that capture the timing of the MRE and the onset of offset of the oldest faulted feature, especially on slow‐to‐moderate slip rate faults, can lead to systematic underestimation of average geologic slip rates.

Seismological Research Letters