USGS Science⌕ Search

Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

Follow this topic with RSS
At least 1,027 records · Page 57Linked to original sources
USGS2021

Temporal seismic velocity variations: Recovery following from the 2019 Mw 7.1 Ridgecrest, California earthquake

We investigated seismic velocity changes ( dv / v ) associated with the 2019 Ridgecrest earthquake sequence with high‐frequency autocorrelations of ambient seismic noise data. Daily autocorrelation functions were computed for the entirety of 2019 and the first quarter of 2020 for broadband stations within the region, including the temporary broadband stations installed during the aftershock deployment. Travel time shifts in the daily autocorrelation functions, relative to the mean autocorrelation waveform, were computed to produce dv / v time series, which are sensitive to the evolving material properties of the shallow crust surrounding the Ridgecrest fault zone (RFZ). A short‐term velocity drop follows the M w 7.1 earthquake at stations in the vicinity of the rupture surface, while those greater than 50 km away showed no such drop. The maximum, absolute changes in seismic velocity are proportional to the logarithm of distance from the fault rupture and to the peak dynamic strain experienced during the earthquake. Near the areas of the highest coseismic slip within the RFZ, seismic velocities recovered over 3 months. However, in the vicinity of the nearby Garlock fault, where triggered slip manifested, and north of the RFZ, seismic velocities recovered within a month. We interpret the seismic velocity changes and their recovery to be largely due to changes in the physical properties of the shallow crust, such as fault zone damage recovery caused by the earthquake rupture process and in response to the large dynamic stresses of passing seismic waves from the mainshock.

California↗
USGS2012

b values and ω −γ seismic source models: Implications for tectonic stress variations along active crustal fault zones and the estimation of high-frequency strong ground motion

In this study the tectonic stress along active crustal fault zones is taken to be of the form , where is the average tectonic stress at depth y and Δσ p ( x, y ) is a seismologically observable, essentially random function of both fault plane coordinates; the stress differences arising in the course of crustal faulting are derived from Δσ p ( x, y ). Empirically known frequency of occurrence statistics, moment-magnitude relationships, and the constancy of earthquake stress drops may be used to infer that the number of earthquakes N of dimension ≥ r is of the form N ∼ 1/ r 2 and that the spectral composition of Δσ p ( x, y ) is of the form , where is the two-dimensional Fourier transform of Δσ p ( x, y ) expressed in radial wave number k . The γ = 2 model of the far-field shear wave displacement spectrum is consistent with the spectral composition , provided that the number of contributions to the spectral representation of the radiated field at frequency ƒ goes as ( k / k 0 ) 2 , consistent with the quasi-static frequency of occurrence relation N ∼ 1/ r 2 ; k 0 is a reference wave number associated with the reciprocal source dimension. Separately, a variety of seismologic observations suggests that the γ = 2 model is the one generally, although certainly not always, applicable to the high-frequency spectral decay of the far-field radiation of earthquakes. In this framework, then, b values near 1, the general validity of the γ = 2 model, and the constancy of earthquake stress drops independent of size are all related to the average spectral composition of . Should one of these change as a result of premonitory effects leading to failure, as has been specifically proposed for b values, it seems likely that one or all of the other characteristics will change as well from their normative values. Irrespective of these associations, the far-field, high-frequency shear radiation for the γ = 2 model in the presence of anelastic attenuation may be interpreted as band-limited, finite duration white noise in acceleration. Its rms value, a rms , is given by the expression a rms = 0.85[2 1/2 (2π) 2 /106] (Δσ/ρ R )(ƒ max /ƒ 0 ) 1/2 , where Δσ is the earthquake stress drop, ρ is density, R is hypocentral distance, ƒ 0 is the spectral corner frequency, and ƒ max is determined by R and specific attenuation 1/ Q . For several reasons, one of which is that it may be estimated in the absence of empirically defined ground motion correlations, a rms holds considerable promise as a measure of high-frequency strong ground motion for engineering purposes.

Conference Paper↗
USGS2012

The growth of geological structures by repeated earthquakes 1. Conceptual framework

In many places, earthquakes with similar characteristics have been shown to recur. If this is common, then relatively small deformations associated with individual earthquake cycles should accumulate over time to create geological structures. Following this paradigm, we show that existing models developed to describe leveling line changes associated with the seismic cycle can be adapted to explain geological features associated with a fault. In these models an elastic layer containing the fault overlies a viscous half-space with a different density. Fault motion associated with an earthquake results in immediate deformation followed by a long period of readjustment as stresses relax in the viscous layer and isostatic equilibrium is restored. Deformation is also caused as a result of the loading and unloading due to sediment deposition and erosion. In this paper, the parameters that control the growth of dip-slip structures are identified. We find that the flexural rigidity of the crust (or the apparent elastic thickness) provides the main control of the width of a structure. The loading due to erosion and deposition of sediment determines the ratio of uplift to subsidence between the two sides of the fault. The flexure due to sediment load is much more important in this respect than whether the fault is normal or reverse in character. We find that, in general, real structures are associated with apparent elastic thicknesses of 4 km or less and thus with very low flexural rigidities.

Journal of Geophysical Research Solid Earth↗
USGS2012

Equivalent strike-slip earthquake cycles in half-space and lithosphere-asthenosphere earth models

By virtue of the images used in the dislocation solution, the deformation at the free surface produced throughout the earthquake cycle by slippage on a long strike-slip fault in an Earth model consisting of an elastic plate (lithosphere) overlying a viscoelastic half-space (asthenosphere) can be duplicated by prescribed slip on a vertical fault embedded in an elastic half-space. For the case in which each earthquake ruptures the entire lithosphere (thickness H ), the half-space equivalent slip rate is as follows: Depth interval 0- H , slip identical to that in lithosphere-asthenosphere model (i.e., abrupt coseismic slip and no subsequent slip); depth interval (2 n −1) H to (2 n +1)H ( n = 1,2,…), slip rate uniform in space and dependent upon time as F n ( t ) exp (− t /τ a ) where F n is a ( n - 1) degree polynomial in t , τ a is twice the asthenosphere relaxation time (η/μ), and t is measured from the instant after the preceding earthquake. The slip rate averaged over the seismic cycle in each depth interval equals the secular rate of relative plate motion. For reasonable values of τ a , slip rates below 5 H do not vary much from that mean value and can be treated as constant. Thus the surface deformation due to the earthquake cycle in the lithosphere-asthenosphere model can be calculated very simply from the half-space model with time-dependent slip in the two depth intervals H −3 H and 3 H −5 H , and uniform slip at a rate equal to the secular relative plate velocity below depth 5 H . Inversion of 1973–1988 geodetic measurements of deformation across the segment of the San Andreas fault in the Transverse Ranges north of Los Angeles for the half-space equivalent slip distribution suggests no significant slip on the fault above 30 km and a uniform slip rate of 36 mm/yr below 30 km. One equivalent lithosphere-asthenosphere model would have a 30-km thick lithosphere and an asthenosphere relaxation time greater than 33 years, but other models are possible.

Journal of Geophysical Research↗
USGS1990

Decrease in deformation rate as a possible precursor to the next Parkfield earthquake

Crustal deformations along fault zones are monitored because they show the accumulation and possible premonitory release of strain energy before major earthquakes. On the Parkfield segment of the San Andreas fault the next mainshock (local magnitude M L =5.7) is expected to occur between 1983 and 1993, based on the relatively regular recurrence of characteristic earthquakes in 1857, 1881, 1901, 1922, 1934 and 1966 1 . A concerted monitoring effort is therefore underway at Parkfield, with the aim of detecting precursory changes in observable quantities that can be measured reliably 2 . Here we report data from two geodetic lines crossing the San Andreas fault near Parkfield which show that the rate of shortening decreased by ˜20% starting in August 1986. This change coincided with a rate decrease of 45% for earthquakes of M L ⩾ 2.1 in the area 3,4 . Both phenomena may have been caused by temporary locking of the fault surface, or a decrease in strain accumulation rate, and they are interpreted as intermediate-term precursors to the next Parkfield earthquake.

California↗
USGS2015

Spatial-temporal variation of low-frequency earthquake bursts near Parkfield, California

Tectonic tremor (TT) and low-frequency earthquakes (LFEs) have been found in the deeper crust of various tectonic environments globally in the last decade. The spatial-temporal behaviour of LFEs provides insight into deep fault zone processes. In this study, we examine recurrence times from a 12-yr catalogue of 88 LFE families with ∼730 000 LFEs in the vicinity of the Parkfield section of the San Andreas Fault (SAF) in central California. We apply an automatic burst detection algorithm to the LFE recurrence times to identify the clustering behaviour of LFEs (LFE bursts) in each family. We find that the burst behaviours in the northern and southern LFE groups differ. Generally, the northern group has longer burst duration but fewer LFEs per burst, while the southern group has shorter burst duration but more LFEs per burst. The southern group LFE bursts are generally more correlated than the northern group, suggesting more coherent deep fault slip and relatively simpler deep fault structure beneath the locked section of SAF. We also found that the 2004 Parkfield earthquake clearly increased the number of LFEs per burst and average burst duration for both the northern and the southern groups, with a relatively larger effect on the northern group. This could be due to the weakness of northern part of the fault, or the northwesterly rupture direction of the Parkfield earthquake.

California↗
USGS1997

Seismic source study of the Racha-Dzhava (Georgia) earthquake from aftershocks and broad-band teleseismic body-wave records: An example of active nappe tectonics

The Racha-Dzhava earthquake ( Ms =7.0) that occurred on 1991 April 29 at 09:12:48.1 GMT in the southern border of the Great Caucasus is the biggest event ever recorded in the region, stronger than the Spitak earthquake ( Ms =6.9) of 1988. A field expedition to the epicentral area was organised and a temporary seismic network of 37 stations was deployed to record the aftershock activity. A very precise image of the aftershock distribution is obtained, showing an elongated cloud oriented N105°, with one branch trending N310° in the western part. The southernmost part extends over 80 km, with the depth ranging from 0 to 15 km, and dips north. The northern branch, which is about 30 km long, shows activity that ranges in depth from 5 to 15 km. The complex thrust dips northwards. A stress-tensor inversion from P -wave first-motion polarities shows a state of triaxial compression, with the major principal axis oriented roughly N-S, the minor principal axis being vertical. Body-waveform inversion of teleseismic seismograms was performed for the main shock, which can be divided into four subevents with a total rupture-time duration of 22 s. The most important part of the seismic moment was released by a gentle northerly dipping thrust. The model is consistent with the compressive tectonics of the region and is in agreement with the aftershock distribution and the stress tensor deduced from the aftershocks. The focal mechanisms of the three largest aftershocks were also inverted from body-wave records. The April 29th ( Ms =6.1) and May 5th ( Ms =5.4) aftershocks have thrust mechanisms on roughly E-W-oriented planes, similar to the main shock. Surprisingly, the June 15th ( Ms =6.2) aftershock shows a thrust fault striking N-S. This mechanism is explained by the structural control of the rupture along the east-dipping geometry of the Dzirula Massif close to the Borzhomi-Kazbeg strike-slip fault. In fact, the orientation and shape of the stress tensor produce a thrust on a N-S oriented plane. Nappe tectonics has been identified as an important feature in the Caucasus, and the source mechanism is consistent with this observation. A hidden fault is present below the nappe, and no large surface breaks were observed due to the main shock. The epicentral region is characterized by sediments that are trapped between two crystalline basements: the Dzirula Massif, which crops out south of Chiatoura, and the Caucasus Main Range north of Oni. Most, if not all, of the rupture is controlled by the thrusting of overlapping, deformed and folded sediments over the Dzirula Massif. This event is another example of blind active faults, with the distinctive feature that the fault plane dips at a gentle angle. The Racha Range is one of the surface expressions of this blind thrust, and its growth is the consequence and evidence of similar earthquakes in the past.

Geophysical Journal International↗
USGS2017

Discriminating between natural vs induced seismicity from long-term deformation history of intraplate faults

To assess whether recent seismicity is induced by human activity or is of natural origin, we analyze fault displacements on high-resolution seismic reflection profiles for two regions in the central United States (CUS): the Fort Worth Basin (FWB) of Texas, and the northern Mississippi embayment (NME). Since 2009 earthquake activity in the CUS has increased dramatically, and numerous publications suggest that this increase is primarily due to induced earthquakes caused by deep-well injection of wastewater, both flowback water from hydrofracturing operations and produced water accompanying hydrocarbon production. Alternatively, some argue that these earthquakes are natural, and that the seismicity increase is a normal variation that occurs over millions of years. Our analysis shows that within the NME, faults deform both Quaternary alluvium and underlying sediments dating from Paleozoic through Tertiary, with displacement increasing with geologic unit age, documenting a long history of natural activity. In the FWB, a region of ongoing wastewater injection, basement faults show deformation of the Proterozoic and Paleozoic units, but little or no deformation of younger strata. Specifically, vertical displacements in the post-Pennsylvanian formations, if any, are below the resolution (~15 m) of the seismic data, far less than expected had these faults accumulated deformation over millions of years. Our results support the assertion that recent FWB earthquakes are of induced origin; this conclusion is entirely independent of analyses correlating seismicity and wastewater injection practices. To our knowledge this is the first study to discriminate natural and induced seismicity using classical structural geology analysis techniques.

Arkansas, Illinois, Mississippi, Missouri, Tenness↗
USGS2026

Predictable seismic cycles result from structural rupture barriers on oceanic transform faults

Earthquakes of magnitude ( M ) >5.5 on oceanic transform faults (OTFs) repeatedly rupture the same locked patches, sometimes quasiperiodically. These patches are separated by “barriers” that halt earthquake propagation and slip mostly aseismically. However, the physical processes governing this systematic behavior remain unclear. We analyzed two barriers along the Gofar transform fault that have arrested ~15 M 6 earthquakes over the past three decades. Ocean bottom seismometer data indicate that the barriers hosted intense microseismicity before the mainshocks and comprise multistrand faults and transtensional stepovers with 100- to 400-m lateral offset. These characteristics contradict earthquake rupture termination models invoking velocity-strengthening friction or large geometric steps and instead point to damage-enhanced porosity and dilatancy-strengthening mechanisms. By isolating rupture segments, the barriers regulate the quasiperiodic recurrence of OTF earthquakes.

Science↗
USGS2004

Land-level changes from a late Holocene earthquake in the northern Puget lowland, Washington

An earthquake, probably generated on the southern Whidbey Island fault zone, caused 1-2 m of ground-surface uplift on central Whidbey Island ???2800-3200 yr ago. The cause of the uplift is a fold that grew coseismically above a blind fault that was the earthquake source. Both the fault and the fold at the fault's tip are imaged on multichannel seismic refection profiles in Puget Sound immediately east of the central Whidbey Island site. Uplift is documented through contrasting histories of relative sea level at two coastal marshes on either side of the fault. Late Holocene shallow-crustal earthquakes of Mw = 6.5-7 pose substantial seismic hazard to the northern Puget Lowland. ?? 2004 Geological Society of America.

Geology↗
USGS2018

Combining conflicting Bayesian models to develop paleoseismic records—An example from the Wasatch Fault Zone, Utah

Bayesian statistical analyses of paleoseismic data result in the probabilistic determination of earthquake times using geochronological data evaluated in the context of a stratigraphic model. However, a fundamental problem in paleoseismology is how to use the Bayesian approach to model sparse and/or conflicting geochronological datasets, such as those derived from sites exhibiting episodic sedimentary and pedogenic processes in moderate‐ to high‐energy environments (e.g., a normal‐faulted alluvial fan). Using paleoseismic data for the Corner Canyon site on the Salt Lake City segment of the Wasatch fault zone (Utah), we develop an approach by which multiple Bayesian models are combined to generate an earthquake history at a site. This approach accommodates mutually exclusive interpretations of the geochronological data and thereby limits the influence of sparse data, stratigraphically inconsistent ages, or a single, subjective model interpretation. For the Corner Canyon site, we integrate four OxCal Bayesian models to generate a chronology of six events between ∼ 4.8 "> ∼ 4.8 and ∼ 0.5    ka "> ∼ 0.5 ka ⁠ . Late Holocene (post‐5 ka) mean recurrence and vertical slip‐rate estimates are ∼ 0.9    ky "> ∼ 0.9 ky (0.7–1.0 ky; 95% confidence) and 1.1    mm / yr "> 1.1 mm / yr ( ⁠ 0.8 – 1.7    mm / yr "> 0.8 – 1.7 mm / yr range), respectively. Although our method increases the uncertainty in the timing of individual earthquakes, it more objectively accounts for potential geochronological errors and different interpretations of stratigraphic age control. By relaxing the need to select a single age model, our approach yields more accurate earthquake‐timing results that will better facilitate evaluations of along‐fault event correlation and earthquake rupture length.

Utah↗
USGS2020

A stress-similarity triggering model for aftershocks of the MW6.4 and MW7.1 Ridgecrest earthquakes

The July 2019 M w "> M w 6.4 and 7.1 Ridgecrest earthquakes triggered numerous aftershocks, including clusters of off‐fault aftershocks in an extensional stepover of the Garlock fault, near the town of Olancha, and near Panamint Valley. The locations of the off‐fault aftershocks are consistent with the stress‐similarity model of triggering, which hypothesizes that aftershocks preferentially occur in areas where the mainshock static stress change tensor is similar in orientation to the background stress tensor. The background stress field is determined from the inversion of earthquake focal mechanisms, with the spatial resolution adapted to the local density of earthquakes. The mainshock static stress change is computed using finite‐source models for the M w "> M w 6.4 foreshock and M w "> M w 7.1 mainshock. I quantify the similarity between these two stress fields using the tensor dot product of the normalized deviatoric stress tensors. The off‐fault aftershocks in the Garlock stepover and the Olancha area fall within lobes of positive stress similarity, whereas the aftershocks near Panamint Valley are partially within a lobe. The cluster in the Garlock fault stepover and the smaller of two clusters near Olancha occur in regions of locally anomalous background stress that results in higher stress similarity. I compute the spatial density of M ≥ 2.0 "> M≥2.0 aftershocks and find that the aftershock density increases as a function of stress similarity, with a factor of ∼ 15 "> ∼15 difference between high stress‐similarity and low stress‐similarity areas. This result is robust with respect to the choice of mainshock model and the uncertainty of the background stress field. The aftershock density varies substantially inside the high stress‐similarity lobes, however, indicating that other variable background conditions, such as material properties, temperature, and fluid pressure, may also be playing a role. Specifically, temperature and fluid pressure conditions might help explain the low rate of aftershocks in the Coso geothermal field.

California↗
USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗
USGS1994

Triggering of the Ms = 5.4 Little Skull Mountain, Nevada, earthquake with dynamic strains

We have developed an approach to test the viability of dynamic strains as a triggering mechanism by quantifying the dynamic strain tensor at seismogenic depths. We focus on the dynamic strains at the hypocenter of the M s = 5.4 Little Skull Mountain ( LSM ), Nevada, earthquake. This event is noteworthy because it is the largest earthquake demonstrably triggered at remote distances (∼280 km) by the M s = 7.4 Landers, California, earthquake and because of its ambiguous association with magmatic activity. Our analysis shows that, if dynamic strains initiate remote triggering, the orientation and modes of faulting most favorable for being triggered by a given strain transient change with depth. The geometry of the most probable LSM fault plane was favorably oriented with respect to the geometry of the dynamic strain tensor. We estimate that the magnitude of the peak dynamic strains at the hypocentral depth of the LSM earthquake were ∼4 μstrain (∼.2 MPa) which are ∼50% smaller than those estimated from velocity seismograms recorded at the surface. We suggest that these strains are too small to cause Mohr-Coulomb style failure unless the fault was prestrained to near failure levels, the fault was exceptionally weak, and/or the dynamic strains trigger other processes that lead to failure.

Nevada↗
USGS2002

Interpreting the earthquake source of the Wabash Valley seismic zone (Illinois, Indiana, and Kentucky) from seismic-reflection, gravity, and magnetic-intensity data

Reprocessing of seismic-reflection data reveals new images of upper- to middle-crustal structures beneath the Wabash Valley seismic zone, located north of the New Madrid seismic zone within the seismically active southern Illinois basin. Four intersecting deep seismic profiles (243 km total) indicate an anomalous, 5–10-km-wide zone of dipping reflections and diffractions below the western flank of the Wabash Valley fault system (WVFS). The zone corresponds in places to gently arched regions of Paleozoic strata. The reflector zone can be interpreted as a result of either (or a combination of) magmatic intrusion or structural deformation. The area encompassing the reflection profiles has experienced several moderate magnitude (3.0 ≤ m bLg ≤ 5.5) earthquakes during the past 50 years, defining the central part of the Wabash Valley seismic zone. The hypocenter of the largest 20th-century earthquake in the central USA midcontinent (9 November 1968, m bLg 5.5) corresponds to the most prominent zone of dipping middle-crustal reflections, just west of the WVFS. Both the focal mechanism (moderately dipping reverse fault) and the expected rupture zone size (∼2.9 km fault length) of this earthquake are consistent with the orientation and size of observed reflectors. Dipping reflector patterns in the Precambrian crust are not collinear with fault surfaces updip in the Paleozoic sedimentary section. This indicates that shallow Paleozoic structures are effectively “decoupled” from deeper, possibly seismogenic structure, which suggests that understanding Paleozoic structure is not the key to understanding the earthquake source. The complex dipping crustal reflectivity beneath the WVFS is typical of Paleozoic continental convergent zones observed elsewhere ( e.g. , Appalachian orogen) and thus may suggest a preserved Proterozoic suture, possibly associated with the distal Grenville orogeny or an older event. Although magnetic intensity, Bouguer gravity, and seismic-reflection data present different means of understanding the deep geology of the area, their integration aids in limiting the number of admissible interpretations. The reflection profiles indicate a variable zone of anomalous crustal structure, including the dipping reflector zone, along a trend of northeast-trending gravity and magnetic highs locally defining the Commerce geophysical lineament (CGL), which is a suspected source of seismic hazard in the central USA midcontinent. Three-dimensional inverse modeling of the residual isostatic gravity anomaly values indicates that the upper part of the dipping reflector zone beneath the CGL lies near an important density boundary in the upper Precambrian crust. The results of our study suggest that the seismogenic source just north of the New Madrid seismic zone consists, in part, of a pre-existing fabric of blind thrusts localized along pre-existing igneous intrusions, locally coincident with the CGL. This suggests that the CGL may be seismogenic in places and thus a potential seismic hazard. The variation in the expression of the CGL using reflection and potential-field data sets is probably partly related to the differing geologic features by which it is expressed, but would be also consistent with its reactivation numerous times under varying stress regimes.

Illinois, Indiana, Kentucky↗
USGS2012

QuakeCaster, an earthquake physics demonstration and exploration tool

A fundamental riddle of earthquake occurrence is that tectonic motions at plate interiors are steady, changing only subtly over millions of years, but at plate boundary faults, the plates are stuck for hundreds of years and then suddenly jerk forward in earthquakes. Why does this happen? The answer, as formulated by Harry F. Reid (Reid 1910, 192) is that the Earth’s crust is elastic, behaving like a very stiff slab of rubber sliding over a substrate of “honey”-like asthenosphere, and that faults are restrained by friction. The crust near the faults—zones of weakness that separate the plates—slowly deforms, building up stress until frictional resistance on the fault is overcome and the fault suddenly slips. For the past century, scientists have sought ways to use this knowledge to forecast earthquakes.

Seismological Research Letters↗
USGS2017

Sculpted by water, elevated by earthquakes—The coastal landscape of Glacier Bay National Park, Alaska

Within Glacier Bay National Park in southeastern Alaska, the Fairweather Fault represents the onshore boundary between two of Earth’s constantly moving tectonic plates: the North American Plate and the Yakutat microplate. Satellite measurements indicate that during the past few decades the Yakutat microplate has moved northwest at a rate of nearly 5 centimeters per year relative to the North American Plate. Motion between the tectonic plates results in earthquakes on the Fairweather Fault during time intervals spanning one or more centuries. For example, in 1958, a 260-kilometer section of the Fairweather Fault ruptured during a magnitude 7.8 earthquake, causing permanent horizontal (as much as 6.5 meters) and vertical (as much as 1 meter) displacement of the ground surface across the fault. Thousands to millions of years of tectonic plate motion, including earthquakes like the one in 1958, raised and shifted the ground surface across the Fairweather Fault, while rivers, glaciers, and ocean waves eroded and sculpted the surrounding landscape along the Gulf of Alaska coast in Glacier Bay National Park.

Alaska↗
USGS1976

A preliminary evaluation of selected earthquake-related geologic hazards in the Kenai Lowland, Alaska

Several major faults exist beneath the Kenai Lowland, and others may be inferred. No surface evidence was found to indicate that any of the faults have been active in Holocene time. The sparse shallow seismicity thus far recorded does not correlate with known faults, nor does it define linear trends suggestive of faulting. Tidal flats at the mouths of the Kenai, Kasilof, and Chickaloon Rivers have been uplifted in the Holocene, however, and the uplift may reflect growing anticlines at depth. Some faults associated with folding may be active, but they are probably too limited in size to generate destructive earthquakes. They could, however, constitute a hazard to manmade structures in their immediate vicinity. Vibrational damage and ground failure during the 1964 earthquake were most extensive north of Tustumena Lake in areas of water-saturated unconsolidated deposits with uneven topography. South of Tustumena Lake, damage resulted mainly from landslides formed along unconfined bluffs and riverbanks, and from submarine sliding and subsidence at Homer Spit. A comparable distribution of surficial effects may be anticipated in future large earthquakes on the Kenai Lowland.

Alaska↗