Geology topics

Thomas E. Parsons

Publications and source records attributed to Thomas E. Parsons.

At least 37 records · Page 2Linked to original sources

USGS2020

Submarine landslide kinematics derived from high-resolution imaging in Port Valdez, Alaska

Submarine landslides caused by strong ground shaking during the M9.2 1964 Great Alaska earthquake generated a tsunami that destroyed much of the old town of Valdez, Alaska, and was responsible for 32 deaths at that location. We explore structural details of the 1964 landslide deposit, as well as landslide deposits from earlier events, in order to characterize kinematics of the landslide process. We present a new high‐resolution seismic reflection data set that images the 1964 landslide deposit and six pre‐1964 deposits with great detail. These deposits are represented by thick packages (~7–23 m) of debris within >500 m of fjord sedimentation above basement. Internal slide structures are associated with distinctive landslide failure mechanisms, including detailed erosional and depositional features and structures resolved within both landslide blocks and distal debris flow layers. Based on comparisons of deposit volume from subbottom structure and differenced bathymetry, we refine prior interpretations of the source of failed material. New data show evidence for basal erosion and reworking of fjord‐floor sedimentation. Additionally, material comprising the 1964 landslide appears to have been translated and deformed by lateral thrusting, rather than having been sourced entirely from upslope evacuation zones. Taking into account these complexities in depositional patterns, we show variations in slide size through Holocene time and relate the history of landslides to the paleoseismic record. Collectively, these new observations demonstrate that Port Valdez has a repeated history of large submarine landslides, which are likely associated with large megathrust earthquakes.

Alaska↗

USGS2020

The predictive skills of elastic Coulomb rate-and-state aftershock forecasts during the 2019 Ridgecrest, California, earthquake sequence

Operational earthquake forecasting protocols commonly use statistical models for their recognized ease of implementation and robustness in describing the short-term spatiotemporal patterns of triggered seismicity. However, recent advances on physics-based aftershock forecasting reveal comparable performance to the standard statistical counterparts with significantly improved predictive skills when fault and stress field heterogeneities are considered. Here, we perform a pseudo-prospective forecasting experiment during the first month of the 2019 Ridgecrest (California) earthquake sequence. We develop seven Coulomb rate-and-state models that couple static stress change estimates with continuum mechanics expressed by the rate-and-state friction laws. Our model parametrization supports a gradually increasing complexity; we start from a preliminary model implementation with simplified slip distributions and spatially homogeneous receiver faults to reach an enhanced one featuring optimized fault constitutive parameters, finite-fault slip models, secondary triggering effects, and spatially heterogenous planes informed by pre-existing ruptures. The data-rich environment of Southern California allows us to test whether incorporating data collected in near real-time during an unfolding earthquake sequence boosts our predictive power. We assess the absolute and relative performance of the forecasts by means of statistical tests used within the Collaboratory for the Study of Earthquake Predictability (CSEP) and compare their skills against a standard benchmark ETAS model for the short (24 hours after the two Ridgecrest mainshocks) and intermediate-term (one month). Stress-based forecasts expect heightened rates along the whole near-fault region and increased expected seismicity rates in Central Garlock Fault. Our comparative model evaluation supports that faulting heterogeneities coupled with secondary triggering effects are the most critical success components behind physics-based forecasts, but also underlines the importance of model updates incorporating near real-time available aftershock data reaching better performance than ETAS models. We explore the physical basis behind our results by investigating the localized shut down of pre-existing normal faults in the Ridgecrest near-source area.

California↗

USGS2020

The role of seismic and slow slip events in triggering the 2018 M7.1 Anchorage earthquake in the Southcentral Alaska subduction zone

The M 7.1 2018 Anchorage earthquake occurred in the bending part of the subducting North Pacific plate near the geometrical barrier formed by the underthrusting Yakutat terrane. We calculate the triggering potential related with stress redistribution from deformation sources including the M 9.2 1964 earthquake coseismic slip, postseismic deformation, slip from regional M > 5 earthquakes, and the cumulative slip of previously detected slow slip events over the past 55 years. We investigate the deeper shallow depth (20–60 km) seismicity response to these perturbations using an epidemic type aftershock sequence model to describe earthquake‐to‐earthquake interactions. The statistical forecast captures the triggered seismicity during the 1983 M 6+ aftershocks in Columbia Bay but performs poorly during the slow slip event period between 1992.0 and 2004.8 that presents a statistically significant rate change ( β , Z > 2; M < 4.0). We find that stress effects from the 1964 postseismic source and the 12‐year‐long slow slip event (~ M 7.8) contribute to the 2018 Anchorage earthquake occurrence and that slow slip events modulate the deeper shallow depth seismicity patterns in the region.

Alaska↗

USGS2020

A new technique to calculate earthquake stress transfer and to forecast aftershocks

Coseismic stress changes have been the primary physical principle used to explain aftershocks and triggered earthquakes. However, this method does not adequately forecast earthquake rates and diverse rupture populations when subjected to formal testing. We show that earthquake forecasts can be impaired by assumptions made in physics-based models, such as the existence of hypothetical optimal faults, and regional scale invariability of the stress field. We compare calculations made under these assumptions along with different realizations of a new conceptual triggering model that features a complete assay of all possible ruptures. In this concept, there always exists a set of theoretical planes that has positive failure stress conditions under a combination of background and coseismic static stress change. In the Earth, all of these theoretical planes may not exist, and if they do, they may not be ready to fail. Thus the actual aftershock plane may not correspond to the plane with the maximum stress change value. This is consistent with observations that mainshocks commonly activate faults with exotic orientations and rakes. Our testing ground is the M=7.2, 2010 El Mayor-Cucapah earthquake sequence that activated multiple diverse fault populations across the USA-Mexico border in California and Baja California. We carry out a retrospective test involving 748 M≥3.0 triggered earthquakes that occurred during a 3-yr period after the mainshock. We find that a probabilistic expression of possible aftershock planes constrained by pre-mainshock rupture patterns is strongly favoured (89% of aftershocks consistent with static stress triggering) versus an optimal fault implementation (35% consistent). Results show that coseismic stress change magnitudes do not necessarily control earthquake triggering, instead we find that the summed background stress and coseismic stress change promotes diverse ruptures. Our model can thus explain earthquake triggering in regions where optimal plane mapping shows coseismic stress reduction.

Bulletin of the Seismological Society of America↗

USGS2019

A combinatorial approach to determine earthquake magnitude distributions on a variable slip-rate fault

Combinatorial methods are used to determine the spatial distribution of earthquake magnitudes on a fault whose slip rate varies along strike. Input to the problem is a finite sample of earthquake magnitudes that span 5 kyr drawn from a truncated Pareto distribution. The primary constraints to the problem are maximum and minimum values around the target slip-rate function indicating where feasible solutions can occur. Two methods are used to determine the spatial distribution of earthquakes: integer programming and the greedy-sequential algorithm. For the integer-programming method, the binary decision vector includes all possible locations along the fault where each earthquake can occur. Once a set of solutions that satisfy the constraints is found, the cumulative slip misfit on the fault is globally minimized relative to the target slip-rate function. The greedy algorithm sequentially places earthquakes to locally optimize slip accumulation. As a case study, we calculate how earthquakes are distributed along the megathrust of the Nankai subduction zone, in which the slip rate varies significantly along strike. For both methods, the spatial distribution of magnitudes depends on slip rate, except for the largest magnitude earthquakes that span multiple sections of the fault. The greedy-sequential algorithm, previously applied to this fault (Parsons et al., 2012), tends to produce smoother spatial distributions and fewer lower magnitude earthquakes in the low slip-rate section of the fault compared to the integer-programming method. Differences in results from the two methods relate to how much emphasis is placed on minimizing the misfit to the target slip rate (integer programming) compared to finding a solution within the slip-rate constraints (greedy sequential). Specifics of the spatial distribution of magnitudes also depend on the shape of the target slip-rate function: i.e. stepped at the section boundaries versus a smooth function. This study isolates the effects of slip-rate variation along a single fault in determining the spatial distribution of earthquake magnitudes, helping to better interpret results from more complex, interconnected fault systems.

Geophysical Journal International↗

USGS2018

Characteristic earthquake magnitude frequency distributions on faults calculated from consensus data in California

An estimate of the expected earthquake rate at all possible magnitudes is needed for seismic hazard forecasts. Regional earthquake magnitude frequency distributions obey a negative exponential law (Gutenberg‐Richter), but it's unclear if individual faults do. We add three new methods to calculate long‐term California earthquake rupture rates to the existing Uniform California Earthquake Rupture Forecast (UCERF3) efforts to assess method and parameter dependence on magnitude frequency results for individual faults. All solutions show strongly characteristic magnitude‐frequency distributions on the San Andreas and other faults, with higher rates of large earthquakes than would be expected from a Gutenberg‐Richter distribution. This is a necessary outcome that results from fitting high fault slip rates under the overall statewide earthquake rate budget. We find that input data choices can affect the nucleation magnitude‐frequency distribution shape for the San Andreas fault; solutions are closer to a Gutenberg‐Richter distribution if the maximum magnitude allowed for earthquakes that occur away from mapped faults (background events) is raised above the consensus threshold of M =7.6, if the moment rate for background events is reduced, or if the overall maximum magnitude is reduced from M =8.5. We also find that participation magnitude‐frequency distribution shapes can be strongly affected by slip‐rate discontinuities along faults that may be artifacts related to segment boundaries.

California↗

USGS2018

Testing earthquake links in Mexico from 1978 up to the 2017 M=8.1 Chiapas and M=7.1 Puebla shocks

The M = 8.1 Chiapas and the M = 7.1 Puebla earthquakes occurred in the bending part of the subducting Cocos plate 11 days and ~600 km apart, a range that puts them well outside the typical aftershock zone. We find this to be a relatively common occurrence in Mexico, with 14% of M > 7.0 earthquakes since 1900 striking more than 300 km apart and within a 2 week interval, not different from a randomized catalog. We calculate the triggering potential caused by crustal stress redistribution from large subduction earthquakes over the last 40 years. There is no evidence that static stress transfer or dynamic triggering from the 8 September Chiapas earthquake promoted the 19 September earthquake. Both recent earthquakes were promoted by past thrust events instead, including delayed afterslip from the 2012 M = 7.5 Oaxaca earthquake. A repeated pattern of shallow thrust events promoting deep intraslab earthquakes is observed over the past 40 years.

Geophysical Research Letters↗

USGS2018

Determining on-fault earthquake magnitude distributions from integer programming

Earthquake magnitude distributions among faults within a fault system are determined from regional seismicity and fault slip rates using binary integer programming. A synthetic earthquake catalog (i.e., list of randomly sampled magnitudes) that spans millennia is first formed, assuming that regional seismicity follows a Gutenberg-Richter relation. Each earthquake in the synthetic catalog can occur on any fault and at any location. The objective is to minimize misfits in the target slip rate for each fault, where slip for each earthquake is scaled from its magnitude. The decision vector consists of binary variables indicating which locations are optimal among all possibilities. Uncertainty estimates in fault slip rates provide explicit upper and lower bounding constraints to the problem. An implicit constraint is that an earthquake can only be located on a fault if it is long enough to contain that earthquake. A general mixed-integer programming solver, consisting of a number of different algorithms, is used to determine the optimal decision vector. A case study is presented for the State of California, where a 4 kyr synthetic earthquake catalog is created and faults with slip ≥3 mm/yr are considered, resulting in >10 6 variables. The optimal magnitude distributions for each of the faults in the system span a rich diversity of shapes, ranging from characteristic to power-law distributions.

Computers & Geosciences↗

USGS2017

Probabilistic tsunami hazard analysis: Multiple sources and global applications

Applying probabilistic methods to infrequent but devastating natural events is intrinsically challenging. For tsunami analyses, a suite of geophysical assessments should be in principle evaluated because of the different causes generating tsunamis (earthquakes, landslides, volcanic activity, meteorological events, and asteroid impacts) with varying mean recurrence rates. Probabilistic Tsunami Hazard Analyses (PTHAs) are conducted in different areas of the world at global, regional, and local scales with the aim of understanding tsunami hazard to inform tsunami risk reduction activities. PTHAs enhance knowledge of the potential tsunamigenic threat by estimating the probability of exceeding specific levels of tsunami intensity metrics (e.g., run-up or maximum inundation heights) within a certain period of time (exposure time) at given locations (target sites); these estimates can be summarized in hazard maps or hazard curves. This discussion presents a broad overview of PTHA, including (i) sources and mechanisms of tsunami generation, emphasizing the variety and complexity of the tsunami sources and their generation mechanisms, (ii) developments in modeling the propagation and impact of tsunami waves, and (iii) statistical procedures for tsunami hazard estimates that include the associated epistemic and aleatoric uncertainties. Key elements in understanding the potential tsunami hazard are discussed, in light of the rapid development of PTHA methods during the last decade and the globally distributed applications, including the importance of considering multiple sources, their relative intensities, probabilities of occurrence, and uncertainties in an integrated and consistent probabilistic framework.

Reviews of Geophysics↗

USGS2017

From coseismic offsets to fault-block mountains

In the Basin and Range extensional province of the western United States, coseismic offsets, under the influence of gravity, display predominantly subsidence of the basin side (fault hanging wall), with comparatively little or no uplift of the mountainside (fault footwall). A few decades later, geodetic measurements [GPS and interferometric synthetic aperture radar (InSAR)] show broad (∼100 km) aseismic uplift symmetrically spanning the fault zone. Finally, after millions of years and hundreds of fault offsets, the mountain blocks display large uplift and tilting over a breadth of only about 10 km. These sparse but robust observations pose a problem in that the coesismic uplifts of the footwall are small and inadequate to raise the mountain blocks. To address this paradox we develop finite-element models subjected to extensional and gravitational forces to study time-varying deformation associated with normal faulting. Stretching the model under gravity demonstrates that asymmetric slip via collapse of the hanging wall is a natural consequence of coseismic deformation. Focused flow in the upper mantle imposed by deformation of the lower crust localizes uplift, which is predicted to take place within one to two decades after each large earthquake. Thus, the best-preserved topographic signature of earthquakes is expected to occur early in the postseismic period.

Proceedings of the National Academy of Sciences of↗

USGS2017

A synoptic view of the Third Uniform California Earthquake Rupture Forecast (UCERF3)

Probabilistic forecasting of earthquake‐producing fault ruptures informs all major decisions aimed at reducing seismic risk and improving earthquake resilience. Earthquake forecasting models rely on two scales of hazard evolution: long‐term (decades to centuries) probabilities of fault rupture, constrained by stress renewal statistics, and short‐term (hours to years) probabilities of distributed seismicity, constrained by earthquake‐clustering statistics. Comprehensive datasets on both hazard scales have been integrated into the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3). UCERF3 is the first model to provide self‐consistent rupture probabilities over forecasting intervals from less than an hour to more than a century, and it is the first capable of evaluating the short‐term hazards that result from multievent sequences of complex faulting. This article gives an overview of UCERF3, illustrates the short‐term probabilities with aftershock scenarios, and draws some valuable scientific conclusions from the modeling results. In particular, seismic, geologic, and geodetic data, when combined in the UCERF3 framework, reject two types of fault‐based models: long‐term forecasts constrained to have local Gutenberg–Richter scaling, and short‐term forecasts that lack stress relaxation by elastic rebound.

Seismological Research Letters↗

USGS2017

A physics-based earthquake simulator and its application to seismic hazard assessment in Calabria (Southern Italy) region

The use of a newly developed earthquake simulator has allowed the production of catalogs lasting 100 kyr and containing more than 100,000 events of magnitudes ≥4.5. The model of the fault system upon which we applied the simulator code was obtained from the DISS 3.2.0 database, selecting all the faults that are recognized on the Calabria region, for a total of 22 fault segments. The application of our simulation algorithm provides typical features in time, space and magnitude behavior of the seismicity, which can be compared with those of the real observations. The results of the physics-based simulator algorithm were compared with those obtained by an alternative method using a slip-rate balanced technique. Finally, as an example of a possible use of synthetic catalogs, an attenuation law has been applied to all the events reported in the synthetic catalog for the production of maps showing the exceedance probability of given values of PGA on the territory under investigation.

Calabria region↗

USGS2017

Nucleation speed limit on remote fluid induced earthquakes

Earthquakes triggered by other remote seismic events are explained as a response to long-traveling seismic waves that temporarily stress the crust. However, delays of hours or days after seismic waves pass through are reported by several studies, which are difficult to reconcile with the transient stresses imparted by seismic waves. We show that these delays are proportional to magnitude and that nucleation times are best fit to a fluid diffusion process if the governing rupture process involves unlocking a magnitude-dependent critical nucleation zone. It is well established that distant earthquakes can strongly affect the pressure and distribution of crustal pore fluids. Earth’s crust contains hydraulically isolated, pressurized compartments in which fluids are contained within low-permeability walls. We know that strong shaking induced by seismic waves from large earthquakes can change the permeability of rocks. Thus, the boundary of a pressurized compartment may see its permeability rise. Previously confined, overpressurized pore fluids may then diffuse away, infiltrate faults, decrease their strength, and induce earthquakes. Magnitude-dependent delays and critical nucleation zone conclusions can also be applied to human-induced earthquakes.

Science Advances↗

USGS2016

Prospective earthquake forecasts at the Himalayan Front after the 25 April 2015 M 7.8 Gorkha Mainshock

When a major earthquake strikes, the resulting devastation can be compounded or even exceeded by the subsequent cascade of triggered seismicity. As the Nepalese recover from the 25 April 2015 shock, knowledge of what comes next is essential. We calculate the redistribution of crustal stresses and implied earthquake probabilities for different periods, from daily to 30 years into the future. An initial forecast was completed before an M 7.3 earthquake struck on 12 May 2015 that enables a preliminary assessment; postforecast seismicity has so far occurred within a zone of fivefold probability gain. Evaluation of the forecast performance, using two months of seismic data, reveals that stress‐based approaches present improved skill in higher‐magnitude triggered seismicity. Our results suggest that considering the total stress field, rather than only the coseismic one, improves the spatial performance of the model based on the estimation of a wide range of potential triggered faults following a mainshock.

Seismological Research Letters↗

USGS2016

Reconstruction of far-field tsunami amplitude distributions from earthquake sources

The probability distribution of far-field tsunami amplitudes is explained in relation to the distribution of seismic moment at subduction zones. Tsunami amplitude distributions at tide gauge stations follow a similar functional form, well described by a tapered Pareto distribution that is parameterized by a power-law exponent and a corner amplitude. Distribution parameters are first established for eight tide gauge stations in the Pacific, using maximum likelihood estimation. A procedure is then developed to reconstruct the tsunami amplitude distribution that consists of four steps: (1) define the distribution of seismic moment at subduction zones; (2) establish a source-station scaling relation from regression analysis; (3) transform the seismic moment distribution to a tsunami amplitude distribution for each subduction zone; and (4) mix the transformed distribution for all subduction zones to an aggregate tsunami amplitude distribution specific to the tide gauge station. The tsunami amplitude distribution is adequately reconstructed for four tide gauge stations using globally constant seismic moment distribution parameters established in previous studies. In comparisons to empirical tsunami amplitude distributions from maximum likelihood estimation, the reconstructed distributions consistently exhibit higher corner amplitude values, implying that in most cases, the empirical catalogs are too short to include the largest amplitudes. Because the reconstructed distribution is based on a catalog of earthquakes that is much larger than the tsunami catalog, it is less susceptible to the effects of record-breaking events and more indicative of the actual distribution of tsunami amplitudes.

Hawaii↗

USGS2016

M≥7 Earthquake rupture forecast and time-dependent probability for the Sea of Marmara region, Turkey

We forecast time-independent and time-dependent earthquake ruptures in the Marmara region of Turkey for the next 30 years using a new fault-segmentation model. We also augment time-dependent Brownian Passage Time (BPT) probability with static Coulomb stress changes (ΔCFF) from interacting faults. We calculate M w > 6.5 probability from 26 individual fault sources in the Marmara region. We also consider a multisegment rupture model that allows higher-magnitude ruptures over some segments of the Northern branch of the North Anatolian Fault Zone (NNAF) beneath the Marmara Sea. A total of 10 different M w =7.0 to M w =8.0 multisegment ruptures are combined with the other regional faults at rates that balance the overall moment accumulation. We use Gaussian random distributions to treat parameter uncertainties (e.g., aperiodicity, maximum expected magnitude, slip rate, and consequently mean recurrence time) of the statistical distributions associated with each fault source. We then estimate uncertainties of the 30-year probability values for the next characteristic event obtained from three different models (Poisson, BPT, and BPT+ΔCFF) using a Monte Carlo procedure. The Gerede fault segment located at the eastern end of the Marmara region shows the highest 30-yr probability, with a Poisson value of 29%, and a time-dependent interaction probability of 48%. We find an aggregated 30-yr Poisson probability of M >7.3 earthquakes at Istanbul of 35%, which increases to 47% if time dependence and stress transfer are considered. We calculate a 2-fold probability gain (ratio time-dependent to time-independent) on the southern strands of the North Anatolian Fault Zone.

Sea if Marmara↗

USGS2016

Vertical deformation associated with normal fault systems evolved over coseismic, postseismic, and multiseismic periods

Vertical deformation of extensional provinces varies significantly and in seemingly contradictory ways. Sparse but robust geodetic, seismic, and geologic observations in the Basin and Range province of the western United States indicate that immediately after an earthquake, vertical change primarily occurs as subsidence of the normal fault hanging wall. A few decades later, a ±100 km wide zone is symmetrically uplifted. The preserved topography of long-term rifting shows bent and tilted footwall flanks rising high above deep basins. We develop finite element models subjected to extensional and gravitational forces to study time-varying deformation associated with normal faulting. We replicate observations with a model that has a weak upper mantle overlain by a stronger lower crust and a breakable elastic upper crust. A 60° dipping normal fault cuts through the upper crust and extends through the lower crust to simulate an underlying shear zone. Stretching the model under gravity demonstrates that asymmetric slip via collapse of the hanging wall is a natural consequence of coseismic deformation. Focused flow in the upper mantle imposed by deformation of the lower crust localizes uplift under the footwall; the breakable upper crust is a necessary model feature to replicate footwall bending over the observed width ( < 10 km), which is predicted to take place within 1-2 decades after each large earthquake. Thus the best-preserved topographic signature of rifting is expected to occur early in the postseismic period. The relatively stronger lower crust in our models is necessary to replicate broader postseismic uplift that is observed geodetically in subsequent decades.

Journal of Geophysical Research↗

USGS2016

Missing link between the Hayward and Rodgers Creek faults

The next major earthquake to strike the ~7 million residents of the San Francisco Bay Area will most likely result from rupture of the Hayward or Rodgers Creek faults. Until now, the relationship between these two faults beneath San Pablo Bay has been a mystery. Detailed subsurface imaging provides definitive evidence of active faulting along the Hayward fault as it traverses San Pablo Bay and bends ~10° to the right toward the Rodgers Creek fault. Integrated geophysical interpretation and kinematic modeling show that the Hayward and Rodgers Creek faults are directly connected at the surface—a geometric relationship that has significant implications for earthquake dynamics and seismic hazard. A direct link enables simultaneous rupture of the Hayward and Rodgers Creek faults, a scenario that could result in a major earthquake ( M = 7.4) that would cause extensive damage and loss of life with global economic impact.

Science Advances↗