Stephen B. DeLong

Oscillations in arid alluvial-channel geometry

Arid alluvial channels on piedmonts and valley floors often exhibit an oscillating pattern of narrow, deeply incised reaches and wide, shallow reaches with a characteristic wavelength. How do these oscillations develop and what controls their wavelengths? To address these questions we developed a two-dimensional numerical model that couples erosion and deposition in a channel bed with cross-sectional widening and narrowing. This model is inherently unstable over a range of spatial scales dependent on the channel width, depth, and slope. In the initial phase of model evolution, wider-than-average channel reaches become zones of distributary flow that aggrade, lose stream power, and further widen in a positive feedback. Simultaneously, narrower-than-average reaches incise, gain stream power, and further narrow. In the second stage of model evolution, this instability is balanced by the diffusive nature of longitudinal profile evolution, and solitary topographic waves propagate in the upstream direction with a characteristic wavelength and amplitude. The model predicts a specific quantitative relationship between the oscillation wavelength and channel width, depth, and slope that is verified by a database of channel geometries in southern Arizona.

Rates and patterns of surface deformation from laser scanning following the South Napa earthquake, California

The A.D. 2014 M6.0 South Napa earthquake, despite its moderate magnitude, caused significant damage to the Napa Valley in northern California (USA). Surface rupture occurred along several mapped and unmapped faults. Field observations following the earthquake indicated that the magnitude of postseismic surface slip was likely to approach or exceed the maximum coseismic surface slip and as such presented ongoing hazard to infrastructure. Using a laser scanner, we monitored postseismic deformation in three dimensions through time along 0.5 km of the main surface rupture. A key component of this study is the demonstration of proper alignment of repeat surveys using point cloud–based methods that minimize error imposed by both local survey errors and global navigation satellite system georeferencing errors. Using solid modeling of natural and cultural features, we quantify dextral postseismic displacement at several hundred points near the main fault trace. We also quantify total dextral displacement of initially straight cultural features. Total dextral displacement from both coseismic displacement and the first 2.5 d of postseismic displacement ranges from 0.22 to 0.29 m. This range increased to 0.33–0.42 m at 59 d post-earthquake. Furthermore, we estimate up to 0.15 m of vertical deformation during the first 2.5 d post-earthquake, which then increased by ∼0.02 m at 59 d post-earthquake. This vertical deformation is not expressed as a distinct step or scarp at the fault trace but rather as a broad up-to-the-west zone of increasing elevation change spanning the fault trace over several tens of meters, challenging common notions about fault scarp development in strike-slip systems. Integrating these analyses provides three-dimensional mapping of surface deformation and identifies spatial variability in slip along the main fault trace that we attribute to distributed slip via subtle block rotation. These results indicate the benefits of laser scanner surveys along active faults and demonstrate that fine-scale variability in fault slip has been missed by traditional earthquake response methods.

Bedrock landscape development modeling: Calibration using field study, geochronology, and digital elevation model analysis

Stream power–based models of bedrock landscape development are effective at producing synthetic topography with realistic fl uvial-network topology and three-dimensional topography, but they are diffi cult to calibrate. This paper examines ways in which fi eld observations, geochronology, and digital elevation model (DEM) data can be used to calibrate a bedrock landscape development model for a specifi c study site. We fi show how uplift rate, bedrock erodibility, and landslide threshold slope are related to steady-state relief, hypsometry, and drainage density for a wide range of synthetic topographies produced by a stream power–based planform landscape development model. Our results indicate that low uplift rates and high erodibility result in low-relief, high drainage density, fl uvially dominated topography, and high uplift rates and low erodibility leads to high-relief, low drainage density, mass wasting–dominated topography. Topography made up of a combination of flchannels and threshold slopes occurs for only a relatively narrow range of model parameters. Using measured values for hypsometric integral, drainage density, and relief, quantitative values of bedrock erodibility can be further constrained, particularly if uplift rates are independently known. We applied these techniques to three sedimentary rock units in the western Transverse Ranges in California that have experienced similar climate, uplift, and incision histories. The 10 Be surface exposure datincision of initially low-relief topography there occurred during the last ~60 k.y. We estimated the relative dependence of drainage area and channel slope on erosion rate in the stream power law from slope-area data, and inferred values for bedrock erodibility ranging from 0.09 to 0.3 m (0.2–0.4) k.y. –1 for the rock types in this study area.

Arroyo channel head evolution in a flash-flood–dominated discontinuous ephemeral stream system

We study whether arroyo channel head retreat in dryland discontinuous ephemeral streams is driven by surface runoff, seepage erosion, mass wasting, or some combination of these hydrogeomorphic processes. We monitored precipitation, overland flow, soil moisture, and headcut migration over several seasonal cycles at two adjacent rangeland channel heads in southern Arizona. Erosion occurred by headward retreat of vertical to overhanging faces, driven dominantly by surface runoff. No evidence exists for erosion caused by shallow-groundwater–related processes, even though similar theater-headed morphologies are sometimes attributed to seepage erosion by emerging groundwater. At our field site, vertical variation in soil shear strength influenced the persistence of the characteristic theater-head form. The dominant processes of erosion included removal of grains and soil aggregates during even very shallow (1–3 cm) overland flow events by runoff on vertical to overhanging channel headwalls, plunge-pool erosion during higher-discharge runoff events, immediate postrunoff wet mass wasting, and minor intra-event dry mass wasting on soil tension fractures developing subparallel to the headwall. Multiple stepwise linear regression indicates that the migration rate is most strongly correlated with flow duration and total precipitation and is poorly correlated with peak flow depth or time-integrated flow depth. The studied channel heads migrated upslope with a self-similar morphologic form under a wide range of hydrological conditions, and the most powerful flash floods were not always responsible for the largest changes in landscape form in this environment.

Tearing the terroir: Details and implications of surface rupture and deformation from the 24 August 2014 M6.0 South Napa earthquake, California

The Mw 6.0 South Napa earthquake of 24 August 2014 caused slip on several active fault strands within the West Napa Fault Zone (WNFZ). Field mapping identified 12.5 km of surface rupture. These field observations, near-field geodesy and space geodesy together provide evidence for more than ~30 km of surface deformation with a relatively complex distribution across a number of subparallel lineaments. Along a ~7 km section north of the epicenter, the surface rupture is confined to a single trace that cuts alluvial deposits, reoccupying a low-slope scarp. The rupture continued northward onto at least four other traces through subparallel ridges and valleys. Postseismic slip exceeded coseismic slip along much of the southern part of the main rupture trace with total slip one year post-event approaching 0.5 meters at locations where only a few centimeters were measured the day of the earthquake. Analysis of airborne interferometric synthetic aperture radar data provides slip distributions along fault traces, indicates connectivity and extent of secondary traces, and confirms that postseismic slip only occurred on the main trace of the fault, perhaps indicating secondary structures ruptured as coseismic triggered slip. Previous mapping identified the WNFZ as a zone of distributed faulting, and this was generally borne out by the complex 2014 rupture pattern. Implications for hazard analysis in similar settings include the need to consider the possibility of complex surface rupture in areas of complex topography, especially where multiple potentially Quaternary-active fault strands can be mapped. This article is protected by copyright. All rights reserved.

Late holocene alluvial history of the Cuyama River, California, USA

Dryland river deposits are archives of past changes in fluvial-system form and process. Chronostratigraphic reconstruction of the late Holocene alluvial history of the Cuyama River in west-central California reveals past spatial and temporal variation in dryland channel form and process. The modern Cuyama River consists of a wide braided reach in the upper drainage basin, a narrower arroyo reach in the middle of the drainage basin, and a bedrock canyon reach in the lower drainage basin that drains to the Santa Maria coastal plain. Along the arroyo reach, late Holocene stratigraphy is well exposed and is the focus of this study. Between ca. 1700 and 350 yr B.P., two widespread deposits of tabular-bedded clay, silt, and fine sand were deposited, separated by a buried soil formed between 950 and 700 yr B.P. Channel incision occurred between 550 and 350 yr B.P. Next, deposition of massive to bedded sands occurred in a pattern alternating between poorly confined deposition on top of older deposits onto the broad valley floor, and channelized deposition along the valley axis inset into older deposits. The superposed deposits now underlie the main valley terrace, and the valley-axis deposits are preserved as inset fill terraces. Historical arroyo cutting then formed a 65-km-long arroyo ca. 150 yr B.P. Based on correlation to regional paleoclimate records, channel aggradation occurred during periods of relative aridity and low peak discharge events, while wet periods, possibly floods after drought, led to fluvial incision. These cycles are superimposed on a transition from a wide, silt- and fine-sand-dominated fluvial system to a modern, narrow, sand– and gravel-dominated arroyo channel. The relationship between times of fluvial process change and climate change in Cuyama Valley bears considerable similarity to other well-studied dryland rivers in the southwestern United States; however, the complex sedimentology and geometry of preserved fluvial deposits suggest that a wider range of fluvial modes occurred along the Cuyama River than has been reported in simpler cut-and-fill–dominated channels elsewhere. These spatial and sedimentological complexities underscore the need to link fluvial deposits and their bounding unconformities along a channel in order to fully understand the spatial and temporal evolution of ancient fluvial systems.

CO2 diffusion into pore spaces limits weathering rate of an experimental basalt landscape

Basalt weathering is a key control over the global carbon cycle, though in situ measurements of carbon cycling are lacking. In an experimental, vegetation-free hillslope containing 330 m 3 of ground basalt scoria, we measured real-time inorganic carbon dynamics within the porous media and seepage flow. The hillslope carbon flux (0.6–5.1 mg C m –2 h –1 ) matched weathering rates of natural basalt landscapes (0.4–8.8 mg C m –2 h –1 ) despite lacking the expected field-based impediments to weathering. After rainfall, a decrease in CO 2 concentration ([CO 2 ]) in pore spaces into solution suggested rapid carbon sequestration but slow reactant supply. Persistent low soil [CO 2 ] implied that diffusion limited CO 2 supply, while when sufficiently dry, reaction product concentrations limited further weathering. Strong influence of diffusion could cause spatial heterogeneity of weathering even in natural settings, implying that modeling studies need to include variable soil [CO 2 ] to improve carbon cycling estimates associated with potential carbon sequestration methods.

Geomorphology, denudation rates, and stream channel profiles reveal patterns of mountain building adjacent to the San Andreas fault in northern California, USA

Relative horizontal motion along strike-slip faults can build mountains when motion is oblique to the trend of the strike-slip boundary. The resulting contraction and uplift pose off-fault seismic hazards, which are often difficult to detect because of the poor vertical resolution of satellite geodesy and difficulty of locating offset datable landforms in active mountain ranges. Sparse geomorphic markers, topographic analyses, and measurement of denudation allow us to map spatiotemporal patterns of uplift along the northern San Andreas fault. Between Jenner and Mendocino, California, emergent marine terraces found southwest of the San Andreas fault record late Pleistocene uplift rates between 0.20 and 0.45 mm yr−1 along much of the coast. However, on the northeast side of the San Andreas fault, a zone of rapid uplift (0.6−1.0 mm yr−1) exists adjacent to the San Andreas fault, but rates decay northeastward as the coast becomes more distant from the San Andreas fault. A newly dated 4.5 Ma shallow-marine deposit located at ∼500 m above sea level (masl) adjacent to the San Andreas fault is warped down to just 150 masl 15 km northeast of the San Andreas fault, and it is exposed at just 60−110 masl to the west of the fault. Landscape denudation rates calculated from abundance of cosmogenic radionuclides in fluvial sediment northeast of, and adjacent to, the San Andreas fault are 0.16−0.29 mm yr−1, but they are only 0.03−0.07 mm yr−1 west of the fault. Basin-average channel steepness and the denudation rates can be used to infer the erosive properties of the underlying bedrock. Calibrated erosion rates can then be estimated across the entire landscape using the spatial distribution of channel steepness with these erosive properties. The lower-elevation areas of this landscape that show high channel steepness (and hence calibrated erosion rate) are distinct from higher-elevation areas with systematically lower channel steepness and denudation rates. These two areas do not appear to be coincident with lithologic contacts. Assuming that changes in rock uplift rates are manifest in channel steepness values as an upstream-propagating kinematic wave that separates high and low channel steepness values, the distance that this transition has migrated vertically provides an estimate of the timing of rock uplift rate increase. This analysis suggests that rock uplift rates along the coast changed from 0.3 to 0.75 mm yr−1 between 450 and 350 ka. This zone of recent, relatively rapid crustal deformation along the plate boundary may be a result of the impingement of relatively strong crust underlying the Gualala block into the thinner, weaker oceanic crust left at the western margin of the North American plate by the westward migration of the subduction zone prior to establishment of the current transform plate boundary. The warped Pliocene marine deposits and the presence of a topographic ridge support the patterns indicated by the channel steepness analyses, and further indicate that the zone of rapid uplift may herald elevated off-fault seismic hazard if this uplift is created by periodic stick-slip motion on contractional structures.

Imaging of earthquake faults using small UAVs as a pathfinder for air and space observations

Large earthquakes cause billions of dollars in damage and extensive loss of life and property. Geodetic and topographic imaging provide measurements of transient and long-term crustal deformation needed to monitor fault zones and understand earthquakes. Earthquake-induced strain and rupture characteristics are expressed in topographic features imprinted on the landscapes of fault zones. Small UAVs provide an efficient and flexible means to collect multi-angle imagery to reconstruct fine scale fault zone topography and provide surrogate data to determine requirements for and to simulate future platforms for air- and space-based multi-angle imaging.