Geology topics

Ryan D. Gold

Publications and source records attributed to Ryan D. Gold.

65 records · Page 4Linked to original sources

USGS2017

Holocene surface-faulting earthquakes at the Spring Lake and North Creek Sites on the Wasatch Fault Zone: Evidence for complex rupture of the Nephi Segment

The Nephi segment of the Wasatch fault zone (WFZ) comprises two fault strands, the northern and southern strands, which have evidence of recurrent late Holocene surface-faulting earthquakes. We excavated paleoseismic trenches across these strands to refine and expand their Holocene earthquake chronologies; improve estimates of earthquake recurrence, displacement, and fault slip rate; and assess whether the strands rupture separately or synchronously in large earthquakes. Paleoseismic data from the Spring Lake site expand the Holocene record of earthquakes on the northern strand: at least five to seven earthquakes ruptured the Spring Lake site at 0.9 ± 0.2 ka (2σ), 2.9 ± 0.7 ka, 4.0 ± 0.5 ka, 4.8 ± 0.8 ka, 5.7 ± 0.8 ka, 6.6 ± 0.7 ka, and 13.1 ± 4.0 ka, yielding a Holocene mean recurrence of ~1.2–1.5 kyr and vertical slip rate of ~0.5–0.8 mm/yr. Paleoseismic data from the North Creek site help refine the Holocene earthquake chronology for the southern strand: at least five earthquakes ruptured the North Creek site at 0.2 ± 0.1 ka (2σ), 1.2 ± 0.1 ka, 2.6 ± 0.9 ka, 4.0 ± 0.1 ka, and 4.7 ± 0.7 ka, yielding a mean recurrence of 1.1–1.3 kyr and vertical slip rate of ~1.9–2.0 mm/yr. We compare these Spring Lake and North Creek data with previous paleoseismic data for the Nephi segment and report late Holocene mean recurrence intervals of ~1.0–1.2 kyr for the northern strand and ~1.1–1.3 kyr for the southern strand. The northern and southern strands have similar late Holocene earthquake histories, which allow for models of both independent and synchronous rupture. However, considering the earthquake timing probabilities and per-event vertical displacements, we have the greatest confidence in the simultaneous rupture of the strands, including rupture of one strand with spillover rupture to the other. Ultimately, our results improve the surface-faulting earthquake history of the Nephi segment and enhance our understanding of how structural barriers influence normal-fault rupture.

Utah↗

USGS2015

A random-walk algorithm for modeling lithospheric density and the role of body forces in the evolution of the Midcontinent Rift

This paper develops a Monte Carlo algorithm for extracting three-dimensional lithospheric density models from geophysical data. Empirical scaling relationships between velocity and density create a 3D starting density model, which is then iteratively refined until it reproduces observed gravity and topography. This approach permits deviations from uniform crustal velocity-density scaling, which provide insight into crustal lithology and prevent spurious mapping of crustal anomalies into the mantle. We test this algorithm on the Proterozoic Midcontinent Rift (MCR), north-central U.S. The MCR provides a challenge because it hosts a gravity high overlying low shear-wave velocity crust in a generally flat region. Our initial density estimates are derived from a seismic velocity/crustal thickness model based on joint inversion of surface-wave dispersion and receiver functions. By adjusting these estimates to reproduce gravity and topography, we generate a lithospheric-scale model that reveals dense middle crust and eclogitized lowermost crust within the rift. Mantle lithospheric density beneath the MCR is not anomalous, consistent with geochemical evidence that lithospheric mantle was not the primary source of rift-related magmas and suggesting that extension occurred in response to far-field stress rather than a hot mantle plume. Similarly, the subsequent inversion of normal faults resulted from changing far-field stress that exploited not only warm, recently faulted crust but also a gravitational potential energy low in the MCR. The success of this density modeling algorithm in the face of such apparently contradictory geophysical properties suggests that it may be applicable to a variety of tectonic and geodynamic problems.

Iowa, Kansas, Michigan, Minnesota, Nebraska, Wisco↗

USGS2015

Rapid characterization of the 2015 M w 7.8 Gorkha, Nepal, earthquake sequence and its seismotectonic context

Earthquake response and related information products are important for placing recent seismic events into context and particularly for understanding the impact earthquakes can have on the regional community and its infrastructure. These tools are even more useful if they are available quickly, ahead of detailed information from the areas affected by such earthquakes. Here we provide an overview of the response activities and related information products generated and provided by the U.S. Geological Survey National Earthquake Information Center in association with the 2015 M 7.8 Gorkha, Nepal, earthquake. This group monitors global earthquakes 24  hrs/day and 7  days/week to provide rapid information on the location and size of recent events and to characterize the source properties, tectonic setting, and potential fatalities and economic losses associated with significant earthquakes. We present the timeline over which these products became available, discuss what they tell us about the seismotectonics of the Gorkha earthquake and its aftershocks, and examine how their information is used today, and might be used in the future, to help mitigate the impact of such natural disasters.

Seismological Research Letters↗

USGS2015

High‐resolution trench photomosaics from image‐based modeling: Workflow and error analysis

Photomosaics are commonly used to construct maps of paleoseismic trench exposures, but the conventional process of manually using image‐editing software is time consuming and produces undesirable artifacts and distortions. Herein, we document and evaluate the application of image‐based modeling (IBM) for creating photomosaics and 3D models of paleoseismic trench exposures, illustrated with a case‐study trench across the Wasatch fault in Alpine, Utah. Our results include a structure‐from‐motion workflow for the semiautomated creation of seamless, high‐resolution photomosaics designed for rapid implementation in a field setting. Compared with conventional manual methods, the IBM photomosaic method provides a more accurate, continuous, and detailed record of paleoseismic trench exposures in approximately half the processing time and 15%–20% of the user input time. Our error analysis quantifies the effect of the number and spatial distribution of control points on model accuracy. For this case study, an ∼87 m 2 exposure of a benched trench photographed at viewing distances of 1.5–7 m yields a model with <2 cm root mean square error (rmse) with as few as six control points. Rmse decreases as more control points are implemented, but the gains in accuracy are minimal beyond 12 control points. Spreading control points throughout the target area helps to minimize error. We propose that 3D digital models and corresponding photomosaics should be standard practice in paleoseismic exposure archiving. The error analysis serves as a guide for future investigations that seek balance between speed and accuracy during photomosaic and 3D model construction.

Bulletin of the Seismological Society of America↗

USGS2015

Evidence for slip partitioning and bimodal slip behavior on a single fault: Surface slip characteristics of the 2013 Mw7.7 Balochistan, Pakistan earthquake

Deformation is commonly accommodated by strain partitioning on multiple, independent strike-slip and dip-slip faults in continental settings of oblique plate convergence. As a corollary, individual faults tend to exhibit one sense of slip – normal, reverse, or strike-slip – until whole-scale changes in boundary conditions reactivate preexisting faults in a new deformation regime. In this study, we show that a single continental fault may instead partition oblique strain by alternatively slipping in a strike-slip or a dip-slip sense during independent fault slip events. We use 0.5 m resolution optical imagery and sub-pixel correlation analysis of the 200+ km 200 + km "> 2013 M w 7.7 Balochistan, Pakistan earthquake to document co-seismic surface slip characteristics and Quaternary tectonic geomorphology along the causative Hoshab fault. We find that the 2013 earthquake, which involved a ∼6:1 strike-slip to dip-slip ratio, ruptured a structurally segmented fault. Quaternary geomorphic indicators of gross fault-zone morphology reveal both reverse-slip and strike-slip deformation in the rupture area of the 2013 earthquake that varies systematically along fault strike despite nearly pure strike-slip motion in 2013. Observations of along-strike variations in range front relief and geomorphic offsets suggest that the Hoshab fault accommodates a substantial reverse component of fault slip in the Quaternary, especially along the southern section of the 2013 rupture. We surmise that Quaternary bimodal slip along the Hoshab fault is promoted by a combination of the arcuate geometry of the Hoshab fault, the frictional weakness of the Makran accretionary prism, and time variable loading conditions from adjacent earthquakes and plate interactions.

Balochistan↗

USGS2015

On- and off-fault deformation associated with the September 2013 Mw7.7 Balochistan earthquake: Implications for geologic slip rate measurements

The 24 September 2013 M w 7.7 Balochistan, Pakistan earthquake ruptured a ~ 200 km-long stretch of the Hoshab fault in southern Pakistan and produced the second-largest lateral surface displacement observed for a continental strike-slip earthquake. We remotely measured surface deformation associated with this event using high-resolution (0.5 m) pre- and post-event satellite optical imagery. We document left lateral, near-field, on-fault offsets (10 m from fault) using 309 laterally offset piercing points, such as streams, terrace risers, and roads. Peak near-field displacement is 13.6 + 2.5 / − 3.4 m. We characterize off-fault deformation by measuring medium- (< 350 m from fault) and far-field (> 350 m from fault) displacement using manual (259 measurements) and automated image cross-correlation methods, respectively. Off-fault peak lateral displacement values are ~ 15 m and exceed on-fault displacement magnitudes for ~ 85% of the rupture length. Our observations suggest that for this rupture, coseismic surface displacement typically increases with distance away from the surface trace of the fault; however, nearly 100% of total surface displacement occurs within a few hundred meters of the primary fault trace. Furthermore, off-fault displacement accounts for, on average, 28% of the total displacement but exhibits a highly heterogeneous along-strike pattern. The best agreement between near-field and far-field displacements generally corresponds to the narrowest fault zone widths. Our analysis demonstrates significant and heterogeneous mismatches between on- and off-fault coseismic deformation, and we conclude that this phenomenon should be considered in hazard models based on geologically determined on-fault slip rates.

Tectonophysics↗

USGS2014

Latest Quaternary paleoseismology and evidence of distributed dextral shear along the Mohawk Valley fault zone, northern Walker Lane, California

The dextral-slip Mohawk Valley fault zone (MVFZ) strikes northwestward along the eastern margin of the Sierra Nevada in the northern Walker Lane. Geodetic block modeling indicates that the MVFZ may accommodate ~3 mm/yr of regional dextral strain, implying that it is the highest slip-rate strike-slip fault in the region; however, only limited geologic data are available to constrain the system’s slip rate and earthquake history. We mapped the MVFZ using airborne lidar data and field observations and identified a site near Sulphur Creek for paleoseismic investigation. At this site, oblique dextral-normal faulting on the steep valley margin has created a closed depression that floods annually during spring snowmelt to form an ephemeral pond. We excavated three fault-perpendicular trenches at the site and exposed pond sediment that interfingers with multiple colluvial packages eroded from the scarp that bounds the eastern side of the pond. We documented evidence for four surface-rupturing earthquakes on this strand of the MVFZ. OxCal modeling of radiocarbon and luminescence ages indicates that these earthquakes occurred at 14.0 ka, 12.8 ka, 5.7 ka, and 1.9 ka. The mean ~4 kyr recurrence interval is inconsistent with slip rates of ~3 mm/yr; these rates imply surface ruptures of more than 10 m per event, which is geologically implausible for the subdued geomorphic expression and 60 km length of the MVFZ. We propose that unidentified structures not yet incorporated into geodetic models may accommodate significant dextral shear across the northern Walker Lane, highlighting the role of distributed deformation in this region.

California↗

USGS2014

Ball-and-socket tectonic rotation during the 2013 M w 7.7 Balochistan earthquake

The September 2013 M w 7.7 Balochistan earthquake ruptured a ∼200-km-long segment of the curved Hoshab fault in southern Pakistan with 10±0.2 m of peak sinistral and ∼1.7±0.8 m of dip slip. This rupture is unusual because the fault dips 60±15° towards the focus of a small circle centered in northwest Pakistan, and, despite a 30° increase in obliquity along strike, the ratios of strike and dip slip remain relatively uniform. Surface displacements and geodetic and teleseismic source inversions quantify a bilateral rupture that propagated rapidly at shallow depths from a transtensional jog near the northern end of the rupture. Static friction prior to rupture was unusually weak (μ<0.05), and friction may have approached zero during dynamic rupture. Here we show that the inward-dipping Hoshab fault defines the northern rim of a structural unit in southeast Makran that rotates – akin to a 2-D ball-and-socket joint – counter-clockwise in response to India's penetration into the Eurasian plate. This rotation accounts for complexity in the Chaman fault system and, in principle, reduces seismic potential near Karachi; nonetheless, these findings highlight deficiencies in strong ground motion equations and tectonic models that invoke Anderson–Byerlee faulting predictions.

Balochistan↗

USGS2014

Holocene and latest Pleistocene paleoseismology of the Salt Lake City segment of the Wasatch Fault Zone, Utah, at the Penrose Drive Trench Site

The Salt Lake City segment (SLCS) of the Wasatch fault zone (WFZ) and the West Valley fault zone (WVFZ) compromise Holocene-active normal faults that bound a large intrabasin graben in northern Salt Lake Valley and have evidence of recurrent, large-magnitude (M ~6-7) surface-faulting earthquakes. However, at the time of this investigation, questions remained regarding the timing, displacement, and recurrence of latest Pleistocene and Holocene earthquakes on the northern SLCS and WVFZ , and whether the WVFZ is seismically independent of, or moves coseismically with, the SLCS. To improve paleoseismic data for the SLCS, we conducted a fault-trench investigation at the Penrose Drive site on the northern SLCS. Two trenches, excavated across an 11-m-high scarp near the northern end of the East Bench fault, exposed colluvial-wedge evidence for fize of six (preferred) surface-faulting earthquakes postdating to Provo-phase shoreline of Lake Bonneville (~14-18 ka). Radiocarbon and luminescence ages support earthquake times at 4.0 ± 0.5 ka (2σ) (PD1), 5.9 ± 0.7 ka (PD2), 7.5 ± 0.8 ka (PD3a), 9.7 ± 1.1 ka (PD3b), 10.9 ± 0.2 ka (PD4), and 12.1 ± 1.6 ka (PD5). At least one additional earthquake occurred at 16.5 ± 1.9 ka (PD6) based on an erosional unconformity that separates deformed Lake Bonneville sily and flat-lying Provo-phase shoreline gravel. Earthquakes PD5-PD1 yield latest Pleistocene (post-Provo) and Holocene mean recurrence intervals of ~1.6 kyr and ~1.7-1.9 kyr, respectively. Using 1.0-1.4 m of per-event vertical displacement for PD5-PD3b corroborate previously identified SLCS earthquakes at 4-10 ka. PD4 and PD5 occurred within an ~8-kyr *17-9 ka) time interval on the SLCS previously interpreted as a period of seismic quiescence, and PD6 possibly corresponds with a previously identified earthquake at ~17 ka (although both events have large timing uncertainties). The Penrose data, when combined with previous paleoseismic results, improve the latest Pleistocene-Holocene earthquake chronology of the SLCS, and demonstrate that the SLCS has been a consistently active source of large-magnitude earthquakes since the latest Pleistocene. At least nine surface-faulting earthquakes (S1-S9) have occurred since the highstand of Lake Bonneville (~18 ka). Where the SLCS earthquake record is most complete (since ~14 ka), per-site estimates of mean recurrence are similar for the latest Pleistocene (post-Provo) (~1.6 kyr), Holocene (~1.6-1.9 kyr), and late Holocene (~1.2-1.4 kyr). These SLCS paleoearthquake data indicate an essentially stable rate of earthquake recurrence since the latest Pleistocene and are important for understanding the earthquake potential of the SLCS, clarifying the seismogenic relation between the SLCS and WVFZ, and forecasting the probabilities of future large-magnitude earthquake in the Wasatch Front region.

Utah↗

USGS2012

Interactive terrain visualization enables virtual field work during rapid scientific response to the 2010 Haiti earthquake

The moment magnitude (M w ) 7.0 12 January 2010 Haiti earthquake is the first major earthquake for which a large-footprint LiDAR (light detection and ranging) survey was acquired within several weeks of the event. Here, we describe the use of virtual reality data visualization to analyze massive amounts (67 GB on disk) of multiresolution terrain data during the rapid scientific response to a major natural disaster. In particular, we describe a method for conducting virtual field work using both desktop computers and a 4-sided, 22 m 3 CAVE immersive virtual reality environment, along with KeckCAVES (Keck Center for Active Visualization in the Earth Sciences) software tools LiDAR Viewer, to analyze LiDAR point-cloud data, and Crusta, for 2.5 dimensional surficial geologic mapping on a bare-earth digital elevation model. This system enabled virtual field work that yielded remote observations of the topographic expression of active faulting within an ∼75-km-long section of the eastern Enriquillo–Plantain Garden fault spanning the 2010 epicenter. Virtual field observations indicated that the geomorphic evidence of active faulting and ancient surface rupture varies along strike. Landform offsets of 6–50 m along the Enriquillo–Plantain Garden fault east of the 2010 epicenter and closest to Port-au-Prince attest to repeated recent surface-rupturing earthquakes there. In the west, the fault trace is well defined by displaced landforms, but it is not as clear as in the east. The 2010 epicenter is within a transition zone between these sections that extends from Grand Goâve in the west to Fayette in the east. Within this transition, between L'Acul (lat 72°40′W) and the Rouillone River (lat 72°35′W), the Enriquillo–Plantain Garden fault is undefined along an embayed low-relief range front, with little evidence of recent surface rupture. Based on the geometry of the eastern and western faults that show evidence of recent surface rupture, we propose that the 2010 event occurred within a stepover that appears to have served as a long-lived boundary between rupture segments, explaining the lack of 2010 surface rupture. This study demonstrates how virtual reality–based data visualization has the potential to transform rapid scientific response by enabling virtual field studies and real-time interactive analysis of massive terrain data sets.

Geosphere↗

USGS2012

A terrestrial lidar-based workflow for determining three-dimensional slip vectors and associated uncertainties

Three-dimensional (3D) slip vectors recorded by displaced landforms are difficult to constrain across complex fault zones, and the uncertainties associated with such measurements become increasingly challenging to assess as landforms degrade over time. We approach this problem from a remote sensing perspective by using terrestrial laser scanning (TLS) and 3D structural analysis. We have developed an integrated TLS data collection and point-based analysis workflow that incorporates accurate assessments of aleatoric and epistemic uncertainties using experimental surveys, Monte Carlo simulations, and iterative site reconstructions. Our scanning workflow and equipment requirements are optimized for single-operator surveying, and our data analysis process is largely completed using new point-based computing tools in an immersive 3D virtual reality environment. In a case study, we measured slip vector orientations at two sites along the rupture trace of the 1954 Dixie Valley earthquake (central Nevada, United States), yielding measurements that are the first direct constraints on the 3D slip vector for this event. These observations are consistent with a previous approximation of net extension direction for this event. We find that errors introduced by variables in our survey method result in <2.5 cm of variability in components of displacement, and are eclipsed by the 10–60 cm epistemic errors introduced by reconstructing the field sites to their pre-erosion geometries. Although the higher resolution TLS data sets enabled visualization and data interactivity critical for reconstructing the 3D slip vector and for assessing uncertainties, dense topographic constraints alone were not sufficient to significantly narrow the wide (<26°) range of allowable slip vector orientations that resulted from accounting for epistemic uncertainties.

Geosphere↗