Ryan D. Gold

Low Quaternary slip rate reconciles geodetic and geologic rates along the Altyn Tagh fault, northwestern Tibet

For more than two decades the slip rate along the active, left-slip Altyn Tagh fault of northwestern Tibet has been disputed, with millennial rates reported to be as much as three times faster than those determined geodetically. This problem is significant because the total offset, plate-boundary length, and age of the Altyn Tagh fault make it the most important single structure accommodating India-Asia convergence north of the Himalayas. Here we show that the central Altyn Tagh fault slipped at only 14–9 mm/a over the past 4–6 ka by tightly bracketing the age of a displaced fluvial terrace riser at Yuemake (88.51°E, 38.19°N). This result contradicts previous latest Quaternary rates and is consistent with those derived from geodetic, paleoseismic, and geologic measurements, and thus resolves the long-standing dispute over the latest Quaternary slip rate along the longest active strike-slip fault in Tibet.

Faulted terrace risers place new constraints on the late Quaternary slip rate for the central Altyn Tagh fault, northwest Tibet

The active, left-lateral Altyn Tagh fault defines the northwestern margin of the Tibetan Plateau in western China. To clarify late Quaternary temporal and spatial variations in slip rate along the central portion of this fault system (85°–90°E), we have more than doubled the number of dated offset markers along the central Altyn Tagh fault. In particular, we determined offset-age relations for seven left-laterally faulted terrace risers at three sites (Kelutelage, Yukuang, and Keke Qiapu) spanning a 140-km-long fault reach by integrating surficial geologic mapping, topographic surveys (total station and tripod–light detection and ranging [T-LiDAR]), and geochronology (radiocarbon dating of organic samples, 230 Th/U dating of pedogenic carbonate coatings on buried clasts, and terrestrial cosmogenic radionuclide exposure age dating applied to quartz-rich gravels). At Kelutelage, which is the westernmost site (37.72°N, 86.67°E), two faulted terrace risers are offset 58 ± 3 m and 48 ± 4 m, and formed at 6.2–6.1 ka and 5.9–3.7 ka, respectively. At the Yukuang site (38.00°N, 87.87°E), four faulted terrace risers are offset 92 ± 12 m, 68 ± 6 m, 55 ± 13 m, and 59 ± 9 m and formed at 24.2–9.5 ka, 6.4–5.0 ka, 5.1–3.9 ka, and 24.2–6.4 ka, respectively. At the easternmost site, Keke Qiapu (38.08°N, 88.12°E), a faulted terrace riser is offset 33 ± 6 m and has an age of 17.1–2.2 ka. The displacement-age relationships derived from these markers can be satisfied by an approximately uniform slip rate of 8– 12 mm/yr. However, additional analysis is required to test how much temporal variability in slip rate is permitted by this data set.

Enabling scientific workflows in virtual reality

To advance research and improve the scientific return on data collection and interpretation efforts in the geosciences, we have developed methods of interactive visualization, with a special focus on immersive virtual reality (VR) environments. Earth sciences employ a strongly visual approach to the measurement and analysis of geologic data due to the spatial and temporal scales over which such data ranges. As observations and simulations increase in size and complexity, the Earth sciences are challenged to manage and interpret increasing amounts of data. Reaping the full intellectual benefits of immersive VR requires us to tailor exploratory approaches to scientific problems. These applications build on the visualization methods strengths, using both 3D perception and interaction with data and models, to take advantage of the skills and training of the geological scientists exploring their data in the VR environment. This interactive approach has enabled us to develop a suite of tools that are adaptable to a range of problems in the geosciences and beyond.

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

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

The moment magnitude (Mw) 7.0 12 January 2010 Haiti earthquake is the fi rst 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 scientifi c response to a major natural disaster. In particular, we describe a method for conducting virtual fi eld 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 pointcloud data, and Crusta, for 2.5 dimensional surfi cial geologic mapping on a bare-earth digital elevation model. This system enabled virtual fi eld 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 fi eld 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 defi ned by displaced landforms, but it is not as clear as in the east. The 2010 epi center 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 un defi ned 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 scientifi c response by enabling virtual fi eld studies and real-time interactive analysis of massive terrain data sets.

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. This article is protected by copyright. All rights reserved.

Interactive mapping on 3‐D terrain models

We present an interactive, real-time mapping system for use with digital elevation models and remotely sensed multispectral imagery that aids geoscientists in the creation and interpretation of geologic/neotectonic maps at length scales of 10 m to 1000 km. Our system provides a terrain visualization of the surface of the Earth or other terrestrial planets by displaying a virtual terrain model generated from a digital elevation model overlain by a color texture generated from orthophotos or satellite imagery. We use a quadtree-based, multiresolution display method to render in real time high-resolution virtual terrain models that span large spatial regions. The system allows users to measure the orientations of geologic surfaces and record their observations by drawing lines directly on the virtual terrain model. In addition, interpretive surfaces can be generated from these drawings and displayed to facilitate understanding of the three-dimensional geometry of geologic surfaces. The main strength of our system is the combination of real-time rendering and interactive mapping performed directly on the virtual terrain model with the ability to navigate the scene while changing viewpoints arbitrarily during mapping. User studies and comparisons with commercially available mapping software show that our system improves mapping accuracy and efficiency and also yields observations that cannot be made with existing systems.

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 systems 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.

Real-time Terrain Mapping

We present an interactive, real-time mapping system for digital elevation maps (DEMs), which allows Earth scientists to map and therefore understand the deformation of the continental crust at length scales of 10m to 1000km. Our system visualizes the surface of the Earth as a 3D surface generated from a DEM, with a color texture generated from a registered multispectral image and vector-based mapping elements draped over it. We use a quadtree-based multiresolution method to be able to render high-resolution terrain mapping data sets of large spatial regions in real time. The main strength of our system is the combination of interactive rendering and interactive mapping directly onto the 3D surface, with the ability to navigate the terrain and to change viewpoints arbitrarily during mapping. User studies and comparisons with commercially available mapping software show that our system improves mapping accuracy and eciency, and also enables qualitatively dierent observations that are not possible to make with existing systems.

An updated stress map of the continental United States reveals heterogeneous intraplate stress

Knowledge of the state of stress in Earth’s crust is key to understanding the forces and processes responsible for earthquakes. Historically, low rates of natural seismicity in the central and eastern United States have complicated efforts to understand intraplate stress, but recent improvements in seismic networks and the spread of human-induced seismicity have greatly improved data coverage. Here, we compile a nationwide stress map based on formal inversions of focal mechanisms that challenges the idea that deformation in continental interiors is driven primarily by broad, uniform stress fields derived from distant plate boundaries. Despite plate-boundary compression, extension dominates roughly half of the continent, and second-order forces related to lithospheric structure appear to control extension directions. We also show that the states of stress in several active eastern United States seismic zones differ significantly from those of surrounding areas and that these anomalies cannot be explained by transient processes, suggesting that earthquakes are focused by persistent, locally derived sources of stress. Such spatially variable intraplate stress appears to justify the current, spatially variable estimates of seismic hazard. Future work to quantify sources of stress, stressing-rate magnitudes and their relationship with strain and earthquake rates could allow prospective mapping of intraplate hazard.Crustal stress in the interior of the United States is spatially variable and largely controlled by local forces, rather than those transmitted from tectonic plate boundaries, according to a map of the continental stress field.