Geology Reports⌕ Search

Geology topics

Adrian Bender

Publications and source records attributed to Adrian Bender.

At least 19 recordsLinked to original sources
USGS2025

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska↗
USGS2024

Oblique contraction along the fastest ocean-continent transform plate boundary focuses rock uplift west of the Fairweather fault, southeast Alaska

Contraction along the Yakutat–North America plate boundary drives 4.6–9.0 mm/year Holocene rock uplift rates along Earth's fastest slipping (≥49 mm/year) ocean–continent transform fault, the Fairweather Fault. Between Icy Point and Lituya Bay, the near-vertical Fairweather fault focuses rock uplift and rapid right-lateral slip by accommodating both vertical and fault-parallel strain during oblique-slip and separate, predominantly strike-slip ruptures. Unusually high uplift rates, indicated by radiocarbon and luminescence dating, result from a 10-km-wide, asymmetric, positive flower structure along a 20°, ∼30-km-long restraining double bend in the Fairweather fault. The principal reverse fault in the flower structure, the offshore, blind Icy Point–Lituya Bay fault, ruptures no more than every 460–1040 years evidenced by uplifted Holocene shorelines. Maximum 3–5 m coseismic uplifts imply 3.1–10 m dip slip per event and earthquake magnitudes of M w 7.0–7.5. The Yakutat block collides obliquely into North America, and our model entails oblique slip on the Fairweather fault with and without corupture on the reverse fault. Oblique slip is evident by vertically offset (>25 m) fluvial and marine terraces and by the primary Fairweather fault strand that strikes >20° to the west of plate-boundary motion.

Alaska↗
USGS2024

Seismicity and anisotropic imaging reveal an active detachment beneath the northern Alaska Range foothills

North of the Denali Fault, the collision between the Yakutat block with North America is accommodated by a fold-thrust belt giving rise to the northern Alaska Range foothills. At the western end, the Kantishna Hills anticline hosts prominent microseismicity and surface deformation, interpreted as active folding of the Kantishna Hills anticline above a midcrustal detachment. We test for this detachment by using anisotropy-aware receiver functions to image fabric contrasts within the crust in context with seismicity. Seismic stations near the crest of the Kantishna Hills anticline and near its southern flank show a single strong contrast in dipping fabric at depths of 12–13 km near microseismicity clustering depths, consistent with a detachment plane beneath the fold. A minimum b -value at 10–13 km depth is consistent with seismicity on the detachment, compatible with the imaged anisotropic contrast, while off-fault seismicity is shallower and deeper with smaller magnitudes. South-dipping imbricate thrusts in schist characterize the northern Alaska Range foothills structure. This supports our interpretation of the observed anisotropy as reflecting SSW–SSE-dipping foliation above a detachment at ∼10–13 km depth that exploits existing crustal weaknesses along subtle fabric contrasts observed in the seismically quiescent region north of the actively deforming belt.

Alaska↗
USGS2024

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon↗
USGS2024

The 2023 Alaska National Seismic Hazard Model

US Geological Survey (USGS) National Seismic Hazard Models (NSHMs) are used extensively for seismic design regulations in the United States and earthquake scenario development, as well as risk assessment and mitigation for both buildings and infrastructure. This 2023 update of the long-term, time-independent Alaska NSHM includes substantial changes to both the earthquake rupture forecast (ERF) and ground motion models (GMMs). The ERF includes numerous additions to the finite-fault model, considers two deformation models, and introduces updated declustering and smoothing algorithms in the gridded background seismicity model. For the Alaska–Aleutian subduction zone, megathrust earthquakes occur on an updated structural and segmentation model, and the moment magnitude (M) 8+ rupture and rate model include a logic tree branch that considers slip rates derived from geodetic models of interface coupling. The megathrust model considers multiple models of down-dip width, and magnitudes are computed using newly developed scaling relations. For subduction intraslab events and subduction interface events with M < 7, the 2023 update uses a smoothed seismicity model with rupture depths derived from Slab2. The 2023 model updates GMMs in all tectonic settings using the recently published Next Generation Attenuation Subduction (NGA-Sub) GMMs for subduction interface and intraslab events, and the NGA-West2 GMMs for active crustal settings. Collectively, additions and updates to the Alaska NSHM result in hazard increases across most of south-central Alaska relative to the previous model, published in 2007. These changes are primarily due to the adoption of updated rate models for the large-magnitude interface events and the NGA-Sub GMMs that have much higher aleatory variability (sigma), consistent with global observations, and that include models of epistemic uncertainty.

Alaska↗
USGS2024

Boulders modulate hillslope-channel coupling in the northern Alaska Range

Active orogens balance tectonic rock uplift with erosion, commonly via river incision coupled to landslide denudation of “threshold” hillslopes, but sediment’s role in this feedback is unclear. We report fluvial geometry, and sediment size, prevalence, and mobility across two ≤600-m-tall gneissic northern Alaska Range anticlines that sustain steep landslide-clad hillslopes but differ 10× in late Pleistocene−recent rock uplift rate. Enigmatically, the river steepens and narrows prominently across the fold experiencing slow surface uplift (∼0.5 mm/yr) but remains low-gradient and wide downstream across the anticline undergoing rapid differential rock uplift (∼5 mm/yr). Frequent bedload mobilization across both folds implies fluvial equilibration to sediment transport despite discrepant channel forms and similarly prevalent hillslope-derived boulders. Boulder prevalence correlates significantly with channel slope and width on the slowly uplifting anticline, but weakly on the rapidly uplifting anticline. Strong correlations across the tectonically quiescent anticline may reflect local incision-suppressing boulder aggradation that forces the channel to steepen and narrow, consistent with field observations. Conversely, weak correlations across the rapidly uplifting anticline imply that boulders may modulate expected tectonic channel adjustment by preferentially aggrading to subdue slope, and deflecting frequently mobile bedload to drive lateral erosion that maintains channel width, steepens adjacent hillslopes, and perpetuates hillslope-channel coupling. Hence, hillslope-derived boulders may occupy important roles in regulating feedbacks between river incision and landslide erosion that differ fundamentally at high and low tectonic rates.

Geology↗
USGS2024

Repeated coseismic uplift of coastal lagoons above the Patton Bay Splay Fault System, Montague Island, Alaska, USA

Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.

Alaska↗
USGS2023

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra↗
USGS2023

Rapid active thrust faulting at the northern Alaska Range front

Plate convergence rates strongly influence seismicity and mountain building inboard of convergent margins, but the distribution and kinematics of structures accommodating farfield convergence can be elusive. In interior Alaska, Yakutat microplate convergence drives late Pleistocene–recent right slip on the Denali fault, but westward-decreasing slip rates leave substantial residual Yakutat motion unaccounted for. Here, we show that Northern Foothills thrust slip beneath the northern Alaska Range absorbs a modern 4.4 mm/yr geodetic velocity gradient equivalent to ~78% of the 5.6 mm/yr residual Yakutat convergence along the central Denali fault. Infrared-stimulated luminescence ages of strath terrace deposits (67–4 ka; six sites) quantify Totatlanika River bedrock incision across the 1947 M w 7.1 thrust earthquake epicentral region. Incision rates increase abruptly from <1 mm/yr to 4.8–5.6 mm/yr above the blind thrust tip near the range front. Rapid slip at 6.7 mm/yr on a steep thrust ramp beneath the northern Alaska Range front accommodates the geodetic gradient, drives rock uplift at rates matching measured incision rates, and implies that large earthquakes like the 1947 event may recur with 500–1400 yr frequency. Results illuminate focused seismogenic strain inboard of a complex convergent margin and prompt reevaluation of Alaska’s neotectonic framework.

Alaska↗
USGS2022

Yukon River incision drove organic carbon burial in the Bering Sea during global climate changes at 2.6 and 1 Ma

River erosion affects the carbon cycle and thus climate by exporting terrigenous carbon to seafloor sediment and by nourishing CO 2 -consuming marine life. The Yukon River–Bering Sea system preserves rare source-to-sink records of these processes across profound changes in global climate during the past 5 million years (Ma). Here, we expand the terrestrial erosion record by dating terraces along the Charley River, Alaska, and explore linkages among previously published Yukon River tributary incision chronologies and Bering Sea sedimentation. Cosmogenic 26 Al / 10 Be "> isochron burial ages of Charley River terraces match previously documented central Yukon River tributary incision from 2.6 to 1.6 Ma during Pliocene–Pleistocene glacial expansion, and at 1.1 Ma during the 1.2–0.7 Ma Middle Pleistocene climate transition. Bering Sea sediments preserve 2–4-fold rate increases of Yukon River-derived continental detritus, terrestrial and marine organic carbon, and silicate microfossil deposition at 2.6–2.1 and 1.1–0.8 Ma. These tightly coupled records demonstrate elevated terrigenous nutrient and carbon export and concomitant Bering Sea productivity in response to climate-forced Yukon River incision. Carbon burial related to accelerated terrestrial erosion may contribute to CO 2 drawdown across the Pliocene–Pleistocene and Middle Pleistocene climate transitions observed in many proxy records worldwide.

Alaska, British Columbia, Northwest Territories, Y↗
USGS2021

Bedrock gorge incision via anthropogenic meander cutoff

Bedrock river-gorge incision represents a fundamental landscape-shaping process, but a dearth of observational data at >10 yr timescales impedes understanding of gorge formation. I quantify 10 2 yr rates and processes of gorge incision using historical records, field observations, and topographic and image analysis of a human-caused bedrock meander cutoff along the North Fork Fortymile River in Alaska (USA). Miners cut off the meander in 1900 CE, abruptly lowering local base level by 6 m and forcing narrowing and steepening of the channel across a knickpoint that rapidly incised upstream. Tectonic quiescence, consistent rock erosivity, and low millennial erosion rates provide ideal boundary conditions for this 10 2 yr gorge-formation experiment. Initial fast knickpoint propagation (23 m/yr; 1900–1903 CE) slowed (4 m/yr; 1903–1981 CE) to diffusion (1981–2019 CE) as knickpoint slope decreased, yielding an ~350-m-long, 6-m-deep gorge within the pre–1900 CE channel. Today, diffusion dominates incision of a 500-m-long knickzone upstream of the gorge, where sediment transport likely limits ongoing adjustments to the anthropogenic cutoff. Results elucidate channel width, slope, discharge, and sediment dynamics consistent with a gradual transition from detachment- to transport-limited incision in fluvial adjustment to local base-level lowering.

Alaska↗
USGS2021

Geomorphic expression and slip rate of the Fairweather fault, southeast Alaska, and evidence for predecessors of the 1958 rupture

Active traces of the southern Fairweather fault were revealed by light detection and ranging (lidar) and show evidence for transpressional deformation between North America and the Yakutat block in southeast Alaska. We map the Holocene geomorphic expression of tectonic deformation along the southern 30 km of the Fairweather fault, which ruptured in the 1958 moment magnitude 7.8 earthquake. Digital maps of surficial geology, geomorphology, and active faults illustrate both strike-slip and dip-slip deformation styles within a 10°–30° double restraining bend where the southern Fairweather fault steps offshore to the Queen Charlotte fault. We measure offset landforms along the fault and calibrate legacy 14 C data to reassess the rate of Holocene strike-slip motion (≥49 mm/yr), which corroborates published estimates that place most of the plate boundary motion on the Fairweather fault. Our slip-rate estimates allow a component of oblique-reverse motion to be accommodated by contractional structures west of the Fairweather fault consistent with geodetic block models. Stratigraphic and structural relations in hand-dug excavations across two active fault strands provide an incomplete paleoseismic record including evidence for up to six surface ruptures in the past 5600 years, and at least two to four events in the past 810 years. The incomplete record suggests an earthquake recurrence interval of ≥270 years—much longer than intervals <100 years implied by published slip rates and expected earthquake displacements. Our paleoseismic observations and map of active traces of the southern Fairweather fault illustrate the complexity of transpressional deformation and seismic potential along one of Earth's fastest strike-slip plate boundaries.

Alaska↗
USGS2021

Extreme Quaternary plate boundary exhumation and strike slip localized along the southern Fairweather fault, Alaska, USA

The Fairweather fault (southeastern Alaska, USA) is Earth’s fastest-slipping intracontinental strike-slip fault, but its long-term role in localizing Yakutat–(Pacific–)North America plate motion is poorly constrained. This plate boundary fault transitions northward from pure strike slip to transpression where it comes onshore and undergoes a <25°, 30-km-long restraining double bend. To the east, apatite (U-Th)/He (AHe) ages indicate that North America exhumation rates increase stepwise from ~0.7 to 1.7 km/m.y. across the bend. In contrast, to the west, AHe age-depth data indicate that extremely rapid 5–10 km/m.y. Yakutat exhumation rates are localized within the bend. Further northwest, Yakutat AHe and zircon (U-Th)/He (ZHe) ages gradually increase from 0.3 to 2.6 Ma over 150 km and depict an interval of extremely rapid >6–8 km/m.y. exhumation rates that increases in age away from the bend. We interpret this migration of rapid, transient exhumation to reflect prolonged advection of the Cenozoic–Cretaceous sedimentary cover of the eastern Yakutat microplate through a stationary restraining bend along the edge of the North America plate. Yakutat cooling ages imply a long-term strike-slip rate (54 ± 6 km/m.y.) that mimics the millennial (53 ± 5 m/k.y.) and decadal (46 mm/yr) rates. Fairweather fault slip can account for all Pacific–North America relative plate motion throughout Quaternary time and indicates stability of highly localized plate boundary strike slip on a single fault where extreme rock uplift rates are persistently localized within a restraining bend.

Alaska↗
USGS2021

Dating by cosmogenic nuclides

Since the 1990s, cosmogenic nuclides have revolutionized the study of Earth surface processes, particularly the understanding of rates and dates. These nuclides, including 3 He, 10 Be, 14 C, 21 Ne, 26 Al, and 36 Cl, enable dating of landforms and the measurement of erosion rates both at the scale of drainage basins and at specific locations on Earth's surface. Cosmogenic nuclides are produced at low rates (several to hundreds of atoms per gram per year) by the interaction of cosmic rays with elements both in the atmosphere and in surficial materials, including in rock and soil. Measuring nuclide concentrations requires elemental separation from source geologic material followed by counting of atoms using sensitive accelerator mass spectrometers. Because nuclide production rates have been quantified, the measured concentration of these nuclides can be interpreted as a near-surface residence time. Here, we review the systematics of commonly used cosmogenic nuclides, describe how they are extracted and measured, and then present case studies focusing on the most commonly measured cosmogenic nuclide, 10 Be. We present common applications such as dating surface features, including moraines and outcrops shaped by glaciation, the use of cosmogenic nuclides for inferring tectonic and erosion processes in drainage basins, and the use of these nuclides to trace sediment sources in drainage basins. When multiple nuclides are measured in one sample, they can be used to model burial and exposure histories in stratigraphic sections. We conclude by exploring what the future might bring in terms of measurements and applications.

Book chapter↗
USGS2020

Late Cenozoic climate change paces landscape adjustments to Yukon River capture

Late Cenozoic cooling and changes in glacial–interglacial cycle tempo are thought to increase global rates of erosion starting ~3 million years ago (Ma). Bedrock rivers set rates and patterns of erosion in most landscapes, but constraints on river response to late Cenozoic climate change remain elusive. Here, we determine cosmogenic isotope and luminescence ages of well-preserved bedrock terraces along the Fortymile River (Yukon River basin) to reconstruct an ~5 Myr history of fluvial adjustment to late Cenozoic climate and Yukon River headwater capture at 2.6 Ma. Post-capture Yukon River downcutting lowered the Fortymile River outlet, forcing subsequent bedrock incision throughout the Fortymile basin in two pulses, from 2.4 to 1.8 Ma and at ~1 Ma. These pulses of incision disrupted longer intervals of slow river channel sedimentation under near-consistent climate forcing from 4.8 to 2.4 Ma and from 1.8 to ~1 Ma. The Fortymile River delivers sediment to the Bering Sea, where provenance and accumulation rate changes since 4.3 Ma match observed variations in incision. Our results link alluviation and incision to late Cenozoic climate steadiness and change, respectively, and support the hypothesis that climate-forced changes in precipitation and runoff fundamentally control the pace of river incision and landscape erosion.

Alaska↗
USGS2019

Ground failure from the Anchorage, Alaska, earthquake of 30 November 2018

Investigation of ground failure triggered by the 2018 M w "> M w Mw 7.1 Anchorage earthquake showed that landslides, liquefaction, and ground cracking all occurred and caused significant damage. Shallow rock falls and rock slides were the most abundant types of landslides, but they occurred in smaller numbers than global models that are based on earthquake magnitude predict; this might result from the 2018 earthquake being an intraslab event. Liquefaction was common in alluvial and intertidal areas; ground deformation probably related to liquefaction damaged numerous houses and port facilities in Anchorage. Ground cracking was pervasive near the edges of slopes in hilly areas and caused perhaps the most significant property damage of all types of ground failure. A complex of slump–earth flows was triggered along coastal bluffs in southern Anchorage where slides also occurred in 1964; the 2018 slides involved both mobilization of new landside material and reactivation of parts of the 1964 landslide deposits. Large translational slides that formed during the 1964 Alaska earthquake showed evidence of deformation along pre‐existing failure surfaces but did not reactivate with new net downslope displacement. Modeling suggests that ground motion in 2018 was of insufficient duration and too high frequency to trigger reactivation of the deep landslides.

Alaska↗
USGS2019

The 30 November 2018 Mw7.1 Anchorage Earthquake

The M w "> M w 7.1 47 km deep earthquake that occurred on 30 November 2018 had deep societal impacts across southcentral Alaska and exhibited phenomena of broad scientific interest. We document observations that point to future directions of research and hazard mitigation. The rupture mechanism, aftershocks, and deformation of the mainshock are consistent with extension inside the Pacific plate near the down‐dip limit of flat‐slab subduction. Peak ground motions &gt; 25 % g "> g >25%g were observed across more than 8000 &#x2009;&#x2009; km 2 "> 8000 km 2 ⁠ , though the most violent near‐fault shaking was avoided because the hypocenter was nearly 50 km below the surface. The ground motions show substantial variation, highlighting the influence of regional geology and near‐surface soil conditions. Aftershock activity was vigorous with roughly 300 felt events in the first six months, including two dozen aftershocks exceeding M 4.5. Broad subsidence of up to 5 cm across the region is consistent with the rupture mechanism. The passage of seismic waves and possibly the coseismic subsidence mobilized ground waters, resulting in temporary increases in stream flow. Although there were many failures of natural slopes and soils, the shaking was insufficient to reactivate many of the failures observed during the 1964 M 9.2 earthquake. This is explained by the much shorter duration of shaking as well as the lower amplitude long‐period motions in 2018. The majority of observed soil failures were in anthropogenically placed fill soils. Structural damage is attributed to both the failure of these emplaced soils as well as to the ground motion, which shows some spatial correlation to damage. However, the paucity of instrumental ground‐motion recordings outside of downtown Anchorage makes these comparisons challenging. The earthquake demonstrated the challenge of issuing tsunami warnings in complex coastal geographies and highlights the need for a targeted tsunami hazard evaluation of the region. The event also demonstrates the challenge of estimating the probabilistic hazard posed by intraslab earthquakes.

Alaska↗