SEARCH · Search USGS

Results for “Bulletin Seismological Society America”

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

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29Linked to original sources

USGS2016

Noise reduction in long‐period seismograms by way of array summing

Long‐period (>100 s period) seismic data can often be dominated by instrumental noise as well as local site noise. When multiple collocated sensors are installed at a single site, it is possible to improve the overall station noise levels by applying stacking methods to their traces. We look at the noise reduction in long‐period seismic data by applying the time–frequency phase‐weighted stacking method of Schimmel and Gallart (2007) as well as the phase‐weighted stacking (PWS) method of Schimmel and Paulssen (1997) to four collocated broadband sensors installed in the quiet Albuquerque Seismological Laboratory underground vault. We show that such stacking methods can improve vertical noise levels by as much as 10 dB over the mean background noise levels at 400 s period, suggesting that greater improvements could be achieved with an array involving multiple sensors. We also apply this method to reduce local incoherent noise on horizontal seismic records of the 2 March 2016 M w 7.8 Sumatra earthquake, where the incoherent noise levels at very long periods are similar in amplitude to the earthquake signal. To maximize the coherency, we apply the PWS method to horizontal data where relative azimuths between collocated sensors are estimated and compared with a simpler linear stack with no azimuthal rotation. Such methods could help reduce noise levels at various seismic stations where multiple high‐quality sensors have been deployed. Such small arrays may also provide a solution to improving long‐period noise levels at Global Seismographic Network stations.

Bulletin of the Seismological Society of America↗

USGS2016

Continuity of the West Napa–Franklin fault zone inferred from guided waves generated by earthquakes following the 24 August 2014 Mw 6.0 South Napa Earthquake

We measure peak ground velocities from fault‐zone guided waves (FZGWs), generated by on‐fault earthquakes associated with the 24 August 2014 M w 6.0 South Napa earthquake. The data were recorded on three arrays deployed across north and south of the 2014 surface rupture. The observed FZGWs indicate that the West Napa fault zone (WNFZ) and the Franklin fault (FF) are continuous in the subsurface for at least 75 km. Previously published potential‐field data indicate that the WNFZ extends northward to the Maacama fault (MF), and previous geologic mapping indicates that the FF extends southward to the Calaveras fault (CF); this suggests a total length of at least 110 km for the WNFZ–FF. Because the WNFZ–FF appears contiguous with the MF and CF, these faults apparently form a continuous Calaveras–Franklin–WNFZ–Maacama (CFWM) fault that is second only in length (∼300 km) to the San Andreas fault in the San Francisco Bay area. The long distances over which we observe FZGWs, coupled with their high amplitudes (2–10 times the S waves) suggest that strong shaking from large earthquakes on any part of the CFWM fault may cause far‐field amplified fault‐zone shaking. We interpret guided waves and seismicity cross sections to indicate multiple upper crustal splays of the WNFZ–FF, including a northward extension of the Southhampton fault, which may cause strong shaking in the Napa Valley and the Vallejo area. Based on travel times from each earthquake to each recording array, we estimate average P ‐, S ‐, and guided‐wave velocities within the WNFZ–FF (4.8–5.7, 2.2–3.2, and 1.1–2.8 km/s, respectively), with FZGW velocities ranging from 58% to 93% of the average S ‐wave velocities.

California↗

USGS2016

Potentially induced earthquakes during the early twentieth century in the Los Angeles Basin

Recent studies have presented evidence that early to mid‐twentieth‐century earthquakes in Oklahoma and Texas were likely induced by fossil fuel production and/or injection of wastewater ( Hough and Page, 2015 ; Frohlich et al. , 2016 ). Considering seismicity from 1935 onward, Hauksson et al. (2015) concluded that there is no evidence for significant induced activity in the greater Los Angeles region between 1935 and the present. To explore a possible association between earthquakes prior to 1935 and oil and gas production, we first revisit the historical catalog and then review contemporary oil industry activities. Although early industry activities did not induce large numbers of earthquakes, we present evidence for an association between the initial oil boom in the greater Los Angeles area and earthquakes between 1915 and 1932, including the damaging 22 June 1920 Inglewood and 8 July 1929 Whittier earthquakes. We further consider whether the 1933 M w 6.4 Long Beach earthquake might have been induced, and show some evidence that points to a causative relationship between the earthquake and activities in the Huntington Beach oil field. The hypothesis that the Long Beach earthquake was either induced or triggered by an foreshock cannot be ruled out. Our results suggest that significant earthquakes in southern California during the early twentieth century might have been associated with industry practices that are no longer employed (i.e., production without water reinjection), and do not necessarily imply a high likelihood of induced earthquakes at the present time.

California↗

USGS2016

Ground motion in the presence of complex Topography II: Earthquake sources and 3D simulations

Eight seismic stations were placed in a linear array with a topographic relief of 222 m over Mission Peak in the east San Francisco Bay region for a period of one year to study topographic effects. Seventy‐two well‐recorded local earthquakes are used to calculate spectral amplitude ratios relative to a reference site. A well‐defined fundamental resonance peak is observed with individual station amplitudes following the theoretically predicted progression of larger amplitudes in the upslope direction. Favored directions of vibration are also seen that are related to the trapping of shear waves within the primary ridge dimensions. Spectral peaks above the fundamental one are also related to topographic effects but follow a more complex pattern. Theoretical predictions using a 3D velocity model and accurate topography reproduce many of the general frequency and time‐domain features of the data. Shifts in spectral frequencies and amplitude differences, however, are related to deficiencies of the model and point out the importance of contributing factors, including the shear‐wave velocity under the topographic feature, near‐surface velocity gradients, and source parameters.

California↗

USGS2017

Uncertainty, variability, and earthquake physics in ground‐motion prediction equations

Residuals between ground‐motion data and ground‐motion prediction equations (GMPEs) can be decomposed into terms representing earthquake source, path, and site effects. These terms can be cast in terms of repeatable (epistemic) residuals and the random (aleatory) components. Identifying the repeatable residuals leads to a GMPE with reduced uncertainty for a specific source, site, or path location, which in turn can yield a lower hazard level at small probabilities of exceedance. We illustrate a schematic framework for this residual partitioning with a dataset from the ANZA network, which straddles the central San Jacinto fault in southern California. The dataset consists of more than 3200 1.15≤ M ≤3 earthquakes and their peak ground accelerations (PGAs), recorded at close distances ( R ≤20 km). We construct a small‐magnitude GMPE for these PGA data, incorporating V S 30 site conditions and geometrical spreading. Identification and removal of the repeatable source, path, and site terms yield an overall reduction in the standard deviation from 0.97 (in ln units) to 0.44, for a nonergodic assumption, that is, for a single‐source location, single site, and single path. We give examples of relationships between independent seismological observables and the repeatable terms. We find a correlation between location‐based source terms and stress drops in the San Jacinto fault zone region; an explanation of the site term as a function of kappa, the near‐site attenuation parameter; and a suggestion that the path component can be related directly to elastic structure. These correlations allow the repeatable source location, site, and path terms to be determined a priori using independent geophysical relationships. Those terms could be incorporated into location‐specific GMPEs for more accurate and precise ground‐motion prediction.

Bulletin of the Seismological Society of America↗

USGS2016

Toppling analysis of the Echo Cliffs precariously balanced rock

Toppling analysis of a precariously balanced rock (PBR) can provide insight into the nature of ground motion that has not occurred at that location in the past and, by extension, can constrain peak ground motions for use in engineering design. Earlier approaches have targeted 2D models of the rock or modeled the rock–pedestal contact using spring‐damper assemblies that require recalibration for each rock. Here, a method to model PBRs in 3D is presented through a case study of the Echo Cliffs PBR. The 3D model is created from a point cloud of the rock, the pedestal, and their interface, obtained using terrestrial laser scanning. The dynamic response of the model under earthquake excitation is simulated using a rigid‐body dynamics algorithm. The veracity of this approach is demonstrated through comparisons against data from shake‐table experiments. Fragility maps for toppling probability of the Echo Cliffs PBR as a function of various ground‐motion parameters, rock–pedestal interface friction coefficient, and excitation direction are presented. These fragility maps indicate that the toppling probability of this rock is low (less than 0.2) for peak ground acceleration (PGA) and peak ground velocity (PGV) lower than 3 m/s 2 and 0.75 m/s, respectively, suggesting that the ground‐motion intensities at this location from earthquakes on nearby faults have most probably not exceeded the above‐mentioned PGA and PGV during the age of the PBR. Additionally, the fragility maps generated from this methodology can also be directly coupled with existing probabilistic frameworks to obtain direct constraints on unexceeded ground motion at a PBR’s location.

California↗

USGS2017

A spatiotemporal clustering model for the Third Uniform California Earthquake Rupture Forecast (UCERF3‐ETAS): Toward an operational earthquake forecast

We, the ongoing Working Group on California Earthquake Probabilities, present a spatiotemporal clustering model for the Third Uniform California Earthquake Rupture Forecast (UCERF3), with the goal being to represent aftershocks, induced seismicity, and otherwise triggered events as a potential basis for operational earthquake forecasting (OEF). Specifically, we add an epidemic‐type aftershock sequence (ETAS) component to the previously published time‐independent and long‐term time‐dependent forecasts. This combined model, referred to as UCERF3‐ETAS, collectively represents a relaxation of segmentation assumptions, the inclusion of multifault ruptures, an elastic‐rebound model for fault‐based ruptures, and a state‐of‐the‐art spatiotemporal clustering component. It also represents an attempt to merge fault‐based forecasts with statistical seismology models, such that information on fault proximity, activity rate, and time since last event are considered in OEF. We describe several unanticipated challenges that were encountered, including a need for elastic rebound and characteristic magnitude–frequency distributions (MFDs) on faults, both of which are required to get realistic triggering behavior. UCERF3‐ETAS produces synthetic catalogs of M ≥2.5 events, conditioned on any prior M ≥2.5 events that are input to the model. We evaluate results with respect to both long‐term (1000 year) simulations as well as for 10‐year time periods following a variety of hypothetical scenario mainshocks. Although the results are very plausible, they are not always consistent with the simple notion that triggering probabilities should be greater if a mainshock is located near a fault. Important factors include whether the MFD near faults includes a significant characteristic earthquake component, as well as whether large triggered events can nucleate from within the rupture zone of the mainshock. Because UCERF3‐ETAS has many sources of uncertainty, as will any subsequent version or competing model, potential usefulness needs to be considered in the context of actual applications.

California↗

USGS2017

Broadband seismic noise attenuation versus depth at the Albuquerque Seismological Laboratory

Seismic noise induced by atmospheric processes such as wind and pressure changes can be a major contributor to the background noise observed in many seismograph stations, especially those installed at or near the surface. Cultural noise such as vehicle traffic or nearby buildings with air handling equipment also contributes to seismic background noise. Such noise sources fundamentally limit our ability to resolve earthquake‐generated signals. Many previous seismic noise versus depth studies focused separately on either high‐frequency ( > 1    Hz "> > 1 Hz ) or low‐frequency ( < 0.05    Hz "> < 0.05 Hz ) bands. In this study, we use modern high‐quality broadband (BB) and very broadband (VBB) seismometers installed at depths ranging from 1.5 to 188 m at the Albuquerque Seismological Laboratory to evaluate noise attenuation as a function of depth over a broad range of frequencies (0.002–50 Hz). Many modern seismometer deployments use BB or VBB seismometers installed at various depths, depending on the application. These depths range from one‐half meter or less in aftershock study deployments, to one or two meters in the Incorporated Research Institutions for Seismology Transportable Array (TA), to a few meters (shallow surface vaults) up to 100 m or more (boreholes) in the permanent observatories of the Global Seismographic Network (GSN). It is important for managers and planners of these and similar arrays and networks of seismograph stations to understand the attenuation of surface‐generated noise versus depth so that they can achieve desired performance goals within their budgets as well as their frequency band of focus. The results of this study will assist in decisions regarding BB and VBB seismometer installation depths. In general, we find that greater installation depths are better and seismometer emplacement in hard rock is better than in soil. Attenuation for any given depth varies with frequency. More specifically, we find that the dependence of depth will be application dependent based on the frequency band and sensitive axes of interest. For quick deployments (like aftershock studies), 1 m may be deep enough to produce good data, especially when the focus is on vertical data where temperature stability fundamentally limits the low‐frequency noise levels and little low‐frequency data will be used. For temporary (medium‐term) deployments (e.g., TA) where low cost can be very important, 2–3 m should be sufficient, but such shallow installations will limit the ability to resolve low‐frequency signals, especially on horizontal components. Of course, one should try for maximum burial depth within the budget when there is interest in using the data for low‐frequency applications. For long‐term deployments like the permanent observatories of the GSN and similar networks, 100–200 m depth in hard rock is desirable to achieve lowest noise, although 30–60 m may be acceptable.

Bulletin of the Seismological Society of America↗

USGS2016

A note on adding viscoelasticity to earthquake simulators

Here, I describe how time‐dependent quasi‐static stress transfer can be implemented in an earthquake simulator code that is used to generate long synthetic seismicity catalogs. Most existing seismicity simulators use precomputed static stress interaction coefficients to rapidly implement static stress transfer in fault networks with typically tens of thousands of fault patches. The extension to quasi‐static deformation, which accounts for viscoelasticity of Earth’s ductile lower crust and mantle, involves the precomputation of additional interaction coefficients that represent time‐dependent stress transfer among the model fault patches, combined with defining and evolving additional state variables that track this stress transfer. The new approach is illustrated with application to a California‐wide synthetic fault network.

Bulletin of the Seismological Society of America↗

USGS2017

Apparent late Quaternary fault slip rate increase in the southwestern Lower Rhine Graben, central Europe

In regions of low strain, long earthquake recurrence intervals (10 4 –10 6 yrs) and erosive processes limit preservation of Quaternary markers suitable for distinguishing whether faults slip at uniform or secularly varying rates. The Lower Rhine graben in the border region of Germany, The Netherlands, and Belgium provides a unique opportunity to explore Quaternary slip‐rate variations in a region of low strain using the basal (2.29±0.29 Ma) and surface (700±80 ka) contacts of the regionally extensive main terrace (“Hauptterrasse”), deposited by the Rhine and Maas Rivers. These surfaces are vertically offset 3–140 m and 0–68 m, respectively, across individual fault strands within a distributed network of northwest‐trending, slow‐slipping (<0.1 mm/yr) normal faults. In this investigation, we construct Quaternary slip histories for the southern Lower Rhine graben faults using new main terrace surface vertical offset measurements made from light detection and ranging (lidar)‐derived bare‐earth digital terrain models, which we synthesize with existing constraints on the offset basal contact of this fluvial deposit ( n =91 collocated sites with displacement constraints). We find that >80% of the sites record an apparent increase in slip rate for the more recent interval from 700 ka to present, which corresponds to a period of increased uplift of the nearby Rhenish Massif and regional volcanism. However, the apparent increase in slip rate could result, in part, from erosion of the footwall surface below the main terrace, leading to an apparent displacement that is smaller than the total vertical offset since the start of the Quaternary. Prior work focused on characterization of these faults as seismic sources in the Lower Rhine graben has preferentially relied on the average fault‐slip rate constrained using the base of the main terrace. We suggest that average fault‐slip rates calculated using the ∼700 ka main terrace surface are subjected to fewer uncertainties and sample a time interval that is more relevant for seismic‐hazard analysis.

Lower Rhine graben↗

USGS2017

An updated geospatial liquefaction model for global application

We present an updated geospatial approach to estimation of earthquake-induced liquefaction from globally available geospatial proxies. Our previous iteration of the geospatial liquefaction model was based on mapped liquefaction surface effects from four earthquakes in Christchurch, New Zealand, and Kobe, Japan, paired with geospatial explanatory variables including slope-derived VS30, compound topographic index, and magnitude-adjusted peak ground acceleration from ShakeMap. The updated geospatial liquefaction model presented herein improves the performance and the generality of the model. The updates include (1) expanding the liquefaction database to 27 earthquake events across 6 countries, (2) addressing the sampling of nonliquefaction for incomplete liquefaction inventories, (3) testing interaction effects between explanatory variables, and (4) overall improving model performance. While we test 14 geospatial proxies for soil density and soil saturation, the most promising geospatial parameters are slope-derived VS30, modeled water table depth, distance to coast, distance to river, distance to closest water body, and precipitation. We found that peak ground velocity (PGV) performs better than peak ground acceleration (PGA) as the shaking intensity parameter. We present two models which offer improved performance over prior models. We evaluate model performance using the area under the curve under the Receiver Operating Characteristic (ROC) curve (AUC) and the Brier score. The best-performing model in a coastal setting uses distance to coast but is problematic for regions away from the coast. The second best model, using PGV, VS30, water table depth, distance to closest water body, and precipitation, performs better in noncoastal regions and thus is the model we recommend for global implementation.

Christchurch, Kobe↗

USGS2017

Modeling of high‐frequency seismic‐wave scattering and propagation using radiative transfer theory

This is a study of the nonisotropic scattering process based on radiative transfer theory and its application to the observation of the M 4.3 aftershock recording of the 2008 Wells earthquake sequence in Nevada. Given a wide range of recording distances from 29 to 320 km, the data provide a unique opportunity to discriminate scattering models based on their distance‐dependent behaviors. First, we develop a stable numerical procedure to simulate nonisotropic scattering waves based on the 3D nonisotropic scattering theory proposed by Sato (1995). By applying the simulation method to the inversion of M 4.3 Wells aftershock recordings, we find that a nonisotropic scattering model, dominated by forward scattering, provides the best fit to the observed high‐frequency direct S waves and S ‐wave coda velocity envelopes. The scattering process is governed by a Gaussian autocorrelation function, suggesting a Gaussian random heterogeneous structure for the Nevada crust. The model successfully explains the common decay of seismic coda independent of source–station locations as a result of energy leaking from multiple strong forward scattering, instead of backscattering governed by the diffusion solution at large lapse times. The model also explains the pulse‐broadening effect in the high‐frequency direct and early arriving S waves, as other studies have found, and could be very important to applications of high‐frequency wave simulation in which scattering has a strong effect. We also find that regardless of its physical implications, the isotropic scattering model provides the same effective scattering coefficient and intrinsic attenuation estimates as the forward scattering model, suggesting that the isotropic scattering model is still a viable tool for the study of seismic scattering and intrinsic attenuation coefficients in the Earth.

Bulletin of the Seismological Society of America↗

USGS2016

Relations between some horizontal‐component ground‐motion intensity measures used in practice

Various measures using the two horizontal components of recorded ground motions have been used in a number of studies that derive ground‐motion prediction equations and construct maps of shaking intensity. We update relations between a number of these measures, including those in Boore et al. (2006) and Boore (2010) , using the large and carefully constructed global database of ground motions from crustal earthquakes in active tectonic regions developed as part of the Pacific Earthquake Engineering Research Center–Next Generation Attenuation‐West2 project. The ratios from the expanded datasets generally agree to within a few percent of the previously published ratios. We also provide some ratios that were not considered before, some of which will be useful in applications such as constructing ShakeMaps. Finally, we compare two important ratios with those from a large central and eastern North American database and from many records from subduction earthquakes in Japan and Taiwan. In general, the ratios from these regions are within several percent of those from crustal earthquakes in active tectonic regions.

Bulletin of the Seismological Society of America↗

USGS2017

Alternative rupture-scaling relationships for subduction interface and other offshore environments

Alternative fault-rupture-scaling relationships are developed for Mw 7.1– 9.5 subduction interface earthquakes using a new database of consistently derived finitefault rupture models from teleseismic inversion. Scaling relationships are derived for rupture area, rupture length, rupture width, maximum slip, and average slip. These relationships apply width saturation for large-magnitude interface earthquakes (approximately Mw >8:6) for which the physical characteristics of subduction zones limit the depth extent of seismogenic rupture, and consequently, the down-dip limit of strong ground motion generation. On average, the down-dip rupture width for interface earthquakes saturates near 200 km (196 km on average). Accordingly, the reinterpretation of rupture-area scaling for subduction interface earthquakes through the use of a bilinear scaling model suggests that rupture asperity area is less well correlated with magnitude for earthquakes Mw >8:6. Consequently, the size of great-magnitude earthquakes appears to be more strongly controlled by the average slip across asperities. The sensitivity of the interface scaling relationships is evaluated against geographic region (or subduction zone) and average dip along the rupture interface to assess the need for correction factors. Although regional perturbations in fault-rupture scaling could be identified, statistical significance analyses suggest there is little rationale for implementing regional correction factors based on the limited number of interface rupture models available for each region. Fault-rupture-scaling relationships are also developed for intraslab (within the subducting slab), extensional outer-rise and offshore strike-slip environments. For these environments, the rupture width and area scaling properties yield smaller dimensions than interface ruptures for the corresponding magnitude. However, average and maximum slip metrics yield larger values than interface events. These observations reflect both the narrower fault widths and higher stress drops in these faulting environments. Although expressing significantly different rupture-scaling properties from earthquakes in subduction environments, the characteristics of offshore strike-slip earthquake ruptures compare similarly to commonly used rupture-scaling relationships for onshore strike-slip earthquakes.

Bulletin of the Seismological Society of America↗

USGS2017

Investigation of late Pleistocene and Holocene activity in the San Gregorio fault zone on the continental slope north of Monterey Canyon, offshore central California

We provide an extensive high‐resolution geophysical, sediment core, and radiocarbon dataset to address late Pleistocene and Holocene fault activity of the San Gregorio fault zone (SGFZ), offshore central California. The SGFZ occurs primarily offshore in the San Andreas fault system and has been accommodating dextral strike‐slip motion between the Pacific and North American plates since the mid‐Miocene. Our study focuses on the SGFZ where it has been mapped through the continental slope north of Monterey Canyon. From 2009 to 2015, the Monterey Bay Aquarium Research Institute collected high‐resolution multibeam bathymetry and chirp sub‐bottom profiles using an autonomous underwater vehicle (AUV). Targeted samples were collected using a remotely operated vehicle (ROV) to provide radiocarbon age constraints. We integrate the high‐resolution geophysical data with radiocarbon dates to reveal Pleistocene seismic horizons vertically offset less than 5 m on nearly vertical faults. These faults are buried by continuous reflections deposited after ∼17.5 ka and likely following erosion during the last sea‐level lowstand ∼21 ka, bracketing the age of faulting to ∼32–21 ka. Clearly faulted horizons are only detected in a small area where mass wasting exhumed older strata to within ∼25 m of the seafloor. The lack of clearly faulted Holocene deposits and possible highly distributed faulting in the study area are consistent with previous interpretations that late Pleistocene and Holocene activity along the SGFZ may decrease to the south. This study illustrates the complexity of the SGFZ, offshore central California, and demonstrates the utility of very high‐resolution data from combined AUV (geophysical)–ROV (seabed sampling) surveys in offshore studies of fault activity.

Bulletin of the Seismological Society of America↗

USGS2017

Finite‐fault Bayesian inversion of teleseismic body waves

Inverting geophysical data has provided fundamental information about the behavior of earthquake rupture. However, inferring kinematic source model parameters for finite‐fault ruptures is an intrinsically underdetermined problem (the problem of nonuniqueness), because we are restricted to finite noisy observations. Although many studies use least‐squares techniques to make the finite‐fault problem tractable, these methods generally lack the ability to apply non‐Gaussian error analysis and the imposition of nonlinear constraints. However, the Bayesian approach can be employed to find a Gaussian or non‐Gaussian distribution of all probable model parameters, while utilizing nonlinear constraints. We present case studies to quantify the resolving power and associated uncertainties using only teleseismic body waves in a Bayesian framework to infer the slip history for a synthetic case and two earthquakes: the 2011 M w 7.1 Van, east Turkey, earthquake and the 2010 M w 7.2 El Mayor–Cucapah, Baja California, earthquake. In implementing the Bayesian method, we further present two distinct solutions to investigate the uncertainties by performing the inversion with and without velocity structure perturbations. We find that the posterior ensemble becomes broader when including velocity structure variability and introduces a spatial smearing of slip. Using the Bayesian framework solely on teleseismic body waves, we find rake is poorly constrained by the observations and rise time is poorly resolved when slip amplitude is low.

Bulletin of the Seismological Society of America↗

USGS2017

Recurrent Holocene movement on the Susitna Glacier Thrust Fault: The structure that initiated the Mw 7.9 Denali Fault earthquake, central Alaska

We conducted a trench investigation and analyzed pre‐ and postearthquake topography to determine the timing and size of prehistoric surface ruptures on the Susitna Glacier fault (SGF), the thrust fault that initiated the 2002 M w 7.9 Denali fault earthquake sequence in central Alaska. In two of our three hand‐excavated trenches, we found clear evidence for a single pre‐2002 earthquake (penultimate earthquake [PE]) and determined an age of 2210±420 cal. B.P. (2 σ ) for this event. We used structure‐from‐motion software to create a pre‐2002‐earthquake digital surface model (DSM) from 1:62,800‐scale aerial photography taken in 1980 and compared this DSM with postearthquake 5‐m/pixel Interferometric Synthetic Aperature Radar topography taken in 2010. Topographic profiles measured from the pre‐earthquake DSM show features that we interpret as fault and fold scarps. These landforms were about the same size as those formed in 2002, so we infer that the PE was similar in size to the initial ( M w 7.2) subevent of the 2002 sequence. A recurrence interval of 2270 yrs and dip slip of ∼4.8 m yield a single‐interval slip rate of ∼1.8 mm/yr. The lack of evidence for pre‐PE deformation indicates probable episodic (clustering) behavior on the SGF that may be related to strain migration among other similarly oriented thrust faults that together accommodate shortening south of the Denali fault. We suspect that slip‐partitioned thrust‐triggered earthquakes may be a common occurrence on the Denali fault system, but documenting the frequency of such events will be very difficult, given the lack of long‐term paleoseismic records, the number of potential thrust‐earthquake sources, and the pervasive glacial erosion in the region.

Alaska↗

USGS2018

On the depth extent of co-seismic rupture

We investigate the implications of deformation experiments for the coseismic down‐dip extent of rupture in quasi‐dynamic, whole‐cycle earthquake models of a fault for which the depth of the transition between seismic and aseisimic fault slip depends on strain rate. The calculations use a dislocation fault model from Tse and Rice (1986) with a vertical strike‐slip orientation, mode III rupture, and variable along‐strike length. Our reference calculation is the original rheological representation of Tse and Rice with a strain‐rate‐independent transition. The primary calculations use two different representations of a strain‐rate‐dependent transition: (1) between rate‐weakening friction and dislocation creep and (2) between rate‐weakening and rate‐strengthening frictions. For both these cases, when fault strength is high (friction between 0.5 and 0.6) and the transition is sharp, coseismic slip extends a small distance (1–2 km) below the fixed temperature (depth) that is commonly used to define the rheological transition at the plate‐motion rate. Thus, coseismic slip occurs below the depth assumed in seismic hazard models using microseismicity or a chosen fixed‐temperature contour. Though significant coseismic slip occurs below the plate‐rate transition depth, the added moment is < 10 % "> < 10 % of the total. The deep extension is a region that is rheologically distinct; for example, deep coseismic slip can produce a stress increase rather than a stress drop. If friction is smaller, the deepening effect and its contribution to moment are larger. For all representations of the transition, average and surface slip increase with the along‐strike rupture length in a manner consistent with the limited data from natural observations. However, this property is not controlled by the assumed fault rheology; instead, it arises because the stiffness of the fault decreases weakly with fault length, an intrinsic and unrealistic property of the particular crustal scale fault model used.

Bulletin of the Seismological Society of America↗