Geology topics

Roger Bilham

Publications and source records attributed to Roger Bilham.

11 recordsLinked to original sources

USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗

USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗

USGS2024

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗

USGS2024

On the provenance of field reports of the 1886 Charleston, South Carolina, earthquake: A seismo-historical whodunnit

Much of what is known about the effects of the 1886 Charleston, South Carolina, earthquake throughout the epicentral region can be attributed to meticulous field investigations by an individual with training in geology and engineering, Earle Sloan ( Clendenin, 1926 ). In a recent study, Bilham and Hough (2024) undertook a detailed analysis of the effects of the earthquake on railroads in the Charleston region, drawing heavily from Sloan’s reports. This exercise identified several inconsistencies in Sloan’s field reports, including understandable measurement imprecision, inferred data entry mistakes, and transcription errors. The study also begged the question, where was Sloan at the time of the mainshock and over the following week? And to what extent did he draw from secondhand information in compiling his reports? On this question Sloan’s reports were sometimes enigmatic, lending themselves to misinterpretation in contemporaneous as well as modern interpretations. Beyond the details that were germane for, and briefly summarized by, the studies of Bilham and Hough (2023 , 2024) , in this report we don our historical seismologist caps to chronicle Sloan’s activities following the earthquake. We summarize our inferences here for the benefit of future scholars who might attempt to retrace either Sloan’s footsteps or our own. This study also serves to highlight Sloan’s singular contributions to earthquake science, which were never published separately.

South Carolina↗

USGS2024

The 1886 Charleston, South Carolina, earthquake: Intensities and ground motions

The 1 September 1886 Charleston, South Carolina, earthquake was one of the largest preinstrumental earthquakes in eastern North America for which extensive contemporaneous observations were documented. The distribution of shaking was mapped shortly after the earthquake, and reconsidered by several authors in the late twentieth century, but has not been reconsidered with a modern appreciation for issues associated with macroseismic data interpretation. Detailed contemporary accounts have also never been used to map the distribution of numerical shaking intensities in the near field. In this study we reconsider macroseismic data from far‐field accounts as well as detailed accounts of damage in the near field, estimating modified Mercalli intensity values at 1297 locations including over 200 definite “not felt” reports that delineate the overall felt extent. We compare the results to the suite of ground‐motion models for eastern North America selected by the National Seismic Hazard Model, using a recently proposed mainshock rupture model and an average site condition for the locations at which intensities are estimated. The comparison supports the moment magnitude estimate, 7.3, from a recently proposed rupture model ( Bilham and Hough, 2023 ). A ShakeMap constrained by model predictions and estimated intensities further illustrates this consistency, which we show is insensitive to rupture model details. Given the uncertainty of calibration relations for magnitudes close to 7, the overall intensity distribution provides a good characterization of shaking but cannot improve the independent moment magnitude estimate. We also identify a previously unrecognized early large aftershock that occurred 9–10 min after the mainshock, for which we estimate magnitude ∼5.6.

South Carolina↗

USGS2023

The 1886 Charleston, South Carolina, Earthquake: Relic railroad offset reveals rupture

In the absence of documented surface rupture during the 1 September 1886 Charleston earthquake, there has been considerable speculation about the location and mechanism of the causative fault. We use an inferred coseismic offset of the South Carolina Railroad and additional numerical constraints to develop an elastic deformation model—a west‐dipping fault following strands of two previously identified faults. The constraints are consistent with a blind rupture with 6.5 ± 0.3 m of dextral slip and 2 ± 0.5 m of reverse slip below 450 m depth. We propose that repeated slip on this fault has raised the Penholoway Marine Terrace >6 m since ∼770 ka. The inferred coseismic slip on the fault in an M w 7.3 earthquake is consistent with the distribution of damage in 1886.

South Carolina↗

USGS2020

Active steady-state creep on a nontectonic normal fault in southeast Utah: Implications for strain release in a rapidly deforming salt system

Characterizing short-term temporal variations of fault creep provides insight into the evolution, mechanics, and strength of fault systems. Using spirit leveling and an extensome- ter, we measured surface displacement of a fault southwest of the Needles District, Canyon- lands National Park, Utah, where extension is driven by differential unloading of a subsur- face salt layer due to incision of the Colorado River. Results show continuous creep at max- imum rates of 0.7±0.2 mm/yr without large temporal variations typical of episodic creep events. Occasional, minor transient events in fault slip velocity coincided with water infil- tration; however, we found no significant relationship between precipitation and transient events. Detailed mapping of regionally widespread, fault-parallel sinkholes provide evidence for dilation of faults at shallow depth, a process that lowers fault strength. We propose con- tinuous slip is related to low fault strength and differential unloading, as opposed to other salt systems where dissolution has been linked to episodic slip.

Southeastern Utah↗

USGS2018

Revisiting earthquakes in the Los Angeles, California, basin during the early instrumental period: Evidence for an association with oil production

A total of seven independent M L ≥ 4.0 earthquakes occurred in the Los Angeles, California, basin, during the early instrumental period between 1932 and 1952, the largest of which was the 1933 Long Beach earthquake. Revising available macroseismic and instrumental data for a total of 6 4.0 ≤ M L ≤ 5.1 events between 1938 and 1944, we conclude that early instrumental locations can be grossly inconsistent with detailed macroseismic data. We use available macroseismic data to revisit event locations. We further present evidence that most if not all of these moderate earthquakes may have been induced by oil production. We quantify the predicted stress change associated with production from eight oil fields in the southwestern Los Angeles basin and show that frictional failure would have been encouraged beneath and at the periphery of high-volume fields, with stress changes upward of 0.1 MPa at 5-km depth. The results suggest that if earthquakes are induced by stress changes associated with production, the magnitudes of events might tend to be limited by the limited spatial extent of lobes of increased Coulomb failure stress. It further appears that the advent of fluid injection recovery methods ( water-flooding ) around 1960 mitigated induced earthquake risk considerably.

California↗

USGS2005

Wagon loads of sand blows in White County, Illinois

Several anecdotal accounts provide compelling evidence that liquefaction occurred at several sites in Illinois during the 1811-1812 New Madrid sequence, as much as 250 km north of the New Madrid seismic zone (NMSZ). At one Wabash Valley location, sand blows are still evident near Big Prairie, Illinois, a location described in a particularly detailed and precise historic account. This account includes descriptions of substantial liquefaction (sand blows) as well as a two-mile-long east-west-trending “crack” along which two feet of south-side-down displacement occurred. An offset can no longer be seen at this location, which has been extensively farmed and plowed for decades. Field reconnaissance verifies many of the details provided in the account, however. We conducted a seismic-reflection experiment at this location and observed a modest offset in the Paleozoic strata at this location. The offset is opposite to that described in the historic account, consistent with the hypothesis that large midcontinent earthquakes occur on faults reactivated in a Holocene stress regime different from the one in which they were formed. Only two explanations can account for these observations: Either large NMSZ events triggered substantial liquefaction at distances greater than hitherto realized, or at least one large “New Madrid” event occurred significantly north of the NMSZ. We explore these possibilities and conclude that, while neither one can be ruled out, several disparate lines of evidence suggest that the 23 January 1812 “New Madrid mainshock” occurred in White County, Illinois, near the location of the m b 5.5 1968 southern Illinois earthquake and recent microearthquake activity.

Illinois↗

USGS2003

A media-based assessment of damage and ground motions from the January 26th, 2001 M 7.6 Bhuj, India earthquake

We compiled available news and internet accounts of damage and other effects from the 26th January, 2001, Bhuj earthquake, and interpreted them to obtain modified Mercalli intensities at over 200 locations throughout the Indian subcontinent. These values are used to map the intensity distribution using a simple mathematical interpolation method. The maps reveal several interesting features. Within the Kachchh region, the most heavily damaged villages are concentrated towards the western edge of the inferred fault, consistent with western directivity. Significant sediment-induced amplification is also suggested at a number of locations around the Gulf of Kachchh to the south of the epicenter. Away from the Kachchh region intensities were clearly amplified significantly in areas that are along rivers, within deltas, or on coastal alluvium such as mud flats and salt pans. In addition we use fault rupture parameters inferred from teleseismic data to predict shaking intensity at distances of 0-1000 km. We then convert the predicted hard rock ground motion parameters to MMI using a relationship (derived from internet-based intensity surveys) that assigns MMI based on the average effects in a region. The predicted MMIs are typically lower by 1-2 units than those estimated from news accounts. This discrepancy is generally consistent with the expected effect of sediment response, but it could also reflect other factors such as a tendency for media accounts to focus on the most dramatic damage, rather than the average effects. Our modeling results also suggest, however, that the Bhuj earthquake generated more high-frequency shaking than is expected for earthquakes of similar magnitude in California, and may therefore have been especially damaging.

Bhuj↗

USGS1996

A slow earthquake sequence on the San Andreas fault

EARTHQUAKES typically release stored strain energy on timescales of the order of seconds, limited by the velocity of sound in rock. Over the past 20 years, observations and laboratory experiments have indicated that capture can also occur more slowly, with durations up to hours. Such events may be important in earthquake nucleation and in accounting for the excess of plate convergence over seismic slip in subduction zones. The detection of events with larger timescales requires near-field deformation measurements. In December 1992, two borehole strainmeters close to the San Andreas fault in California recorded a slow strain event of about a week in duration, and we show here that the strain changes were produced by a slow earthquake sequence (equivalent magnitude 4.8) with complexity similar to that of regular earthquakes. The largest earthquakes associated with these slow events were small (local magnitude 3.7) and contributed negligible strain release. The importance of slow earthquakes in the seismogenic process remains an open question, but these observations extend the observed timescale for slow events by two orders of magnitude.

Nature↗