SEARCH · USGS Science

Results for “GSA Bulletin”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13Linked to original sources

USGS2012

Faulting within the Mount St. Helens conduit and implications for volcanic earthquakes

The 2004–2008 eruption of Mount St. Helens produced seven dacite spines mantled by cataclastic fault rocks, comprising an outer fault core and an inner damage zone. These fault rocks provide remarkable insights into the mechanical processes that accompany extrusion of degassed magma, insights that are useful in forecasting dome-forming eruptions. The outermost part of the fault core consists of finely comminuted fault gouge that is host to 1- to 3-mm-thick layers of extremely fine-grained slickenside-bearing ultracataclasite. Interior to the fault core, there is an ∼2-m-thick damage zone composed of cataclastic breccia and sheared dacite, and interior to the damage zone, there is massive to flow-banded dacite lava of the spine interior. Structures and microtextures indicate entirely brittle deformation, including rock breakage, tensional dilation, shearing, grain flow, and microfaulting, as well as gas and fluid migration through intergranular pores and fractures in the damage zone. Slickenside lineations and consistent orientations of Riedel shears indicate upward shear of the extruding spines against adjacent conduit wall rocks. Paleomagnetic directions, demagnetization paths, oxide mineralogy, and petrology indicate that cataclasis took place within dacite in a solidified steeply dipping volcanic conduit at temperatures above 500 °C. Low water content of matrix glass is consistent with brittle behavior at these relatively high temperatures, and the presence of tridymite indicates solidification depths of <1 km. Cataclasis was coincident with the eruption’s seismogenic zone at <1.5 km. More than a million small and low-frequency “drumbeat” earthquakes with coda magnitudes (M d ) <2.0 and frequencies <5 Hz occurred during the 2004–2008 eruption. Our field data provide a means with which to estimate slip-patch dimensions for shear planes and to compare these with estimates of slip patches based on seismic moments and shear moduli for dacite rock and granular fault gouge. Based on these comparisons, we find that aseismic creep is achieved by micron-scale displacements on Riedel shears and by granular flow, whereas the drumbeat earthquakes require millimeter to centimeter displacements on relatively large (e.g., ∼1000 m 2 ) slip patches, possibly along observed extensive principal shear zones within the fault core but probably not along the smaller Riedel shears. Although our field and structural data are compatible with stick-slip models, they do not rule out seismic and infrasound models that call on resonance of steam-filled fractures to generate the drumbeat earthquakes. We suggest that stick-slip and gas release processes may be coupled, and that regardless of the source mechanism, the distinctive drumbeat earthquakes are proving to be an effective precursor for dome-forming eruptions. Our data document a continuous cycle of deformation along the conduit margins beginning with episodes of fracture in the damage zone and followed by transfer of motion to the fault core. We illustrate the cycle of deformation using a hypothetical cross section of the Mount St. Helens conduit, extending from the surface to the depth of magmatic solidification.

Washington↗

USGS2013

Origin and lateral migration of linear dunes in the Qaidam Basin of NW China revealed by dune sediments, internal structures, and optically stimulated luminescence ages, with implications for linear dunes on Titan: discussion

Zhou et al. (2012) proposed that longitudinal dunes in the Qaidam Basin, China, formed like yardangs: by erosion into sediment that was not deposited by those dunes. Because erosion occurs on the upwind flanks of most migrating dunes ( Rubin and Hunter, 1982 , 1985 ), the key to demonstrating a yardang-like origin is to show that the dunes did not deposit the strata that they contain. Zhou et al. made this argument by proposing that: (1) The dunes have not deposited cross-strata in the past 810 yr. (2) Cross-bedding within the dunes was not deposited by the dunes on the present-day land surface, but rather by older dunes that had a different morphology. (3) The present dunes are a later generation, “most likely of erosional origin similar to yardangs with orientations controlled by strikes of joints,” (p. 1147). (4) Rates of deflation in the dune field have been extremely high for the past 810–2440 yr. This commentary reviews these conclusions, reviews contradictory observations, and considers alternative interpretations.

GSA Bulletin↗

USGS2014

Structural controls on geothermal circulation in Surprise Valley, California: A re-evaluation of the Lake City fault zone

Faults and fractures play an important role in the circulation of geothermal fluids in the crust, and the nature of that role varies according to structural setting and state of stress. As a result, detailed geologic and geophysical mapping that relates thermal springs to known structural features is essential to modeling geothermal systems. Published maps of Surprise Valley in northeastern California suggest that the “Lake City fault” or “Lake City fault zone” is a significant structural feature, cutting obliquely across the basin and connecting thermal springs across the valley. Newly acquired geophysical data (audio-magnetotelluric, gravity, and magnetic), combined with existing geochemical and geological data, suggest otherwise. We examine potential field profiles and resistivity models that cross the mapped Lake City fault zone. While there are numerous geophysical anomalies that suggest subsurface structures, they mostly do not coincide with the mapped traces of the Lake City fault zone, nor do they show a consistent signature in gravity, magnetics, or resistivities that would suggest a through-going fault that would promote connectivity through lateral fluid flow. Instead of a single, continuous fault, we propose the presence of a deformation zone associated with the growth of the range-front Surprise Valley fault. The implication for geothermal circulation is that this is a zone of enhanced porosity but lacks length-wise connectivity that could conduct fluids across the valley. Thermal fluid circulation is most likely controlled primarily by interactions between N-S–trending normal faults.

California↗

USGS2014

Geologic and physiographic controls on bed-material yield, transport, and channel morphology for alluvial and bedrock rivers, western Oregon

The rivers of western Oregon have diverse forms and characteristics, with channel substrates ranging from continuous alluvial gravel to bare bedrock. Analysis of several measurable morphologic attributes of 24 valley reaches on 17 rivers provides a basis for comparing nonalluvial and alluvial channels. Key differences are that alluvial reaches have greater bar area, greater migration rates, and show systematic correlation among variables relating grain size to bed-material transport capacity. We relate these differences between channel types to bed-material transport rates as derived from a coupled regional analysis of empirical sediment yield measurements and physical experiments of clast attrition during transport. This sediment supply analysis shows that overall bed-material transport rates for western Oregon are chiefly controlled by (1) lithology and basin slope, which are the key factors for bed-material supply into the stream network, and (2) lithologic control of bed-material attrition from in-transport abrasion and disintegration. This bed-material comminution strongly affects bed-material transport in the study area, reducing transport rates by 50%–90% along the length of the larger rivers in the study area. A comparison of the bed-material transport estimates with the morphologic analyses shows that alluvial gravel-bed channels have systematic and bounding relations between bed-material transport rate and attributes such as bar area and local transport capacity. By contrast, few such relations are evident for nonalluvial rivers with bedrock or mixed-bed substrates, which are apparently more influenced by local controls on channel geometry and sediment supply. At the scale of western Oregon, the physiographic and lithologic controls on the balance between bed-material supply and transport capacity exert far-reaching influence on the distribution of alluvial and nonalluvial channels and their consequently distinctive morphologies and behaviors—differences germane for understanding river response to tectonics and environmental perturbations, as well as for implementing effective restoration and monitoring strategies.

California, Nevada, Oregon, Washington↗

USGS2014

Early to Middle Ordovician back-arc basin in the southern Appalachian Blue Ridge: characteristics, extent, and tectonic significance

Fault-dismembered segments of a distinctive, extensive, highly allochthonous, and tectonically significant Ordovician (ca. 480–460 Ma) basin, which contains suites of bimodal metavolcanic rocks, associated base metal deposits, and thick immature deep-water (turbiditic) metasediments, occur in parts of the southern Appalachian Talladega belt, eastern Blue Ridge, and Inner Piedmont of Alabama, Georgia, and North and South Carolina. The basin's predominantly metasedimentary strata display geochemical and isotopic evidence of a mixed provenance, including an adjacent active volcanic arc and a provenance of mica (clay)-rich sedimentary and felsic plutonic rocks consistent with Laurentian (Grenvillian) upper-crustal continental rocks and their passive-margin cover sequences. Geochemical characteristics of the subordinate intercalated bimodal metavolcanic rocks indicate formation in a suprasubduction environment, most likely a back-arc basin, whereas characteristics of metasedimentary units suggest deposition above Neoproterozoic rift and outer-margin lower Paleozoic slope and rise sediments within a marginal basin along Ordovician Laurentia's Iapetus margin. This tectonic setting indicates that southernmost Appalachian Ordovician orogenesis (Taconic orogeny) began as an extensional accretionary orogen along the outer margin of Laurentia, rather than in an exotic (non-Laurentian) arc collisional setting. B-type subduction polarity requires that the associated arc-trench system formed southeast of the palinspastic position of the back-arc basin. This scenario can explain several unique features of the southern Appalachian Taconic orogen, including: the palinspastic geographic ordering of key tectonic elements (i.e., back-arc, arc, etc.), and a lack of (1) an obducted arc sensu stricto on the Laurentian margin, (2) widespread Ordovician regional metamorphism, and (3) Taconic klippen to supply detritus to the Taconic foreland basin.

Alabama, Georgia, North Carolina, South Carolina↗

USGS2014

Calibration of a conodont apatite-based Ordovician 87 Sr/ 86 Sr curve to biostratigraphy and geochronology: Implications for stratigraphic resolution

The Ordovician 87 Sr/ 86 Sr isotope seawater curve is well established and shows a decreasing trend until the mid-Katian. However, uncertainties in calibration of this curve to biostratigraphy and geochronology have made it difficult to determine how the rates of 87 Sr/ 86 Sr decrease may have varied, which has implications for both the stratigraphic resolution possible using Sr isotope stratigraphy and efforts to model the effects of Ordovician geologic events. We measured 87 Sr/ 86 Sr in conodont apatite in North American Ordovician sections that are well studied for conodont biostratigraphy, primarily in Nevada, Oklahoma, the Appalachian region, and Ohio Valley. Our results indicate that conodont apatite may provide an accurate medium for Sr isotope stratigraphy and strengthen previous reports that point toward a significant increase in the rate of fall in seawater 87 Sr/ 86 Sr during the Middle Ordovician Darriwilian Stage. Our 87 Sr/ 86 Sr results suggest that Sr isotope stratigraphy will be most useful as a high-resolution tool for global correlation in the mid-Darriwilian to mid-Sandbian, when the maximum rate of fall in 87 Sr/ 86 Sr is estimated at &sim;5.0–10.0 × 10 –5 per m.y. Variable preservation of conodont elements limits the precision for individual stratigraphic horizons. Replicate conodont analyses from the same sample differ by an average of &sim;4.0 × 10 –5 (the 2σ standard deviation is 6.2 × 10 –5 ), which in the best case scenario allows for subdivision of Ordovician time intervals characterized by the highest rates of fall in 87 Sr/ 86 Sr at a maximum resolution of &sim;0.5–1.0 m.y. Links between the increased rate of fall in 87 Sr/ 86 Sr beginning in the mid-late Darriwilian ( Phragmodus polonicus to Pygodus serra conodont zones) and geologic events continue to be investigated, but the coincidence with a long-term rise in sea level (Sauk-Tippecanoe megasequence boundary) and tectonic events (Taconic orogeny) in North America provides a plausible explanation for the changing magnitude and 87 Sr/ 86 Sr of the riverine Sr flux to the oceans.

GSA Bulletin↗

USGS2014

Stratigraphic and microfossil evidence for a 4500-year history of Cascadia subduction zone earthquakes and tsunamis at Yaquina River estuary, Oregon, USA

The Sallys Bend swamp and marsh area on the central Oregon coast onshore of the Cascadia subduction zone contains a sequence of buried coastal wetland soils that extends back ∼4500 yr B.P. The upper 10 of the 12 soils are represented in multiple cores. Each soil is abruptly overlain by a sandy deposit and then, in most cases, by greater than 10 cm of mud. For eight of the 10 buried soils, times of soil burial are constrained through radiocarbon ages on fine, delicate detritus from the top of the buried soil; for two of the buried soils, diatom and foraminifera data constrain paleoenvironment at the time of soil burial. We infer that each buried soil represents a Cascadia subduction zone earthquake because the soils are laterally extensive and abruptly overlain by sandy deposits and mud. Preservation of coseismically buried soils occurred from 4500 yr ago until ∼500–600 yr ago, after which preservation was compromised by cessation of gradual relative sea-level rise, which in turn precluded drowning of marsh soils during instances of coseismic subsidence. Based on grain-size and microfossil data, sandy deposits overlying buried soils accumulated immediately after a subduction zone earthquake, during tsunami incursion into Sallys Bend. The possibility that the sandy deposits were sourced directly from landslides triggered upstream in the Yaquina River basin by seismic shaking was discounted based on sedimentologic, microfossil, and depositional site characteristics of the sandy deposits, which were inconsistent with a fluvial origin. Biostratigraphic analyses of sediment above two buried soils—in the case of two earthquakes, one occurring shortly after 1541–1708 cal. yr B.P. and the other occurring shortly after 3227–3444 cal. yr B.P.—provide estimates that coseismic subsidence was a minimum of 0.4 m. The average recurrence interval of subduction zone earthquakes is 420–580 yr, based on an ∼3750–4050-yr-long record and seven to nine interearthquake intervals. The comparison of the Yaquina Bay earthquake record to similar records at other Cascadia coastal sites helps to define potential patterns of rupture for different earthquakes, although inherent uncertainty in dating precludes definitive statements about rupture length during earthquakes. We infer that in the first half of the last millennia, the northern Oregon part of the subduction zone had a different rupture history than the southern Oregon part of the subduction zone, and we also infer that at ca. 1.6 ka, two earthquakes closely spaced in time together ruptured a length of the megathrust that extends at least from southwestern Washington to southern Oregon.

Oregon↗

USGS2016

Active faulting on the Wallula fault zone within the Olympic-Wallowa lineament, Washington State, USA

The Wallula fault zone is an integral feature of the Olympic-Wallowa lineament, an ∼500-km-long topographic lineament oblique to the Cascadia plate boundary, extending from Vancouver Island, British Columbia, to Walla Walla, Washington. The structure and past earthquake activity of the Wallula fault zone are important because of nearby infrastructure, and also because the fault zone defines part of the Olympic-Wallowa lineament in south-central Washington and suggests that the Olympic-Wallowa lineament may have a structural origin. We used aeromagnetic and ground magnetic data to locate the trace of the Wallula fault zone in the subsurface and map a quarry exposure of the Wallula fault zone near Finley, Washington, to investigate past earthquakes along the fault. We mapped three main packages of rocks and unconsolidated sediments in an ∼10-m-high quarry exposure. Our mapping suggests at least three late Pleistocene earthquakes with surface rupture, and an episode of liquefaction in the Holocene along the Wallula fault zone. Faint striae on the master fault surface are subhorizontal and suggest reverse dextral oblique motion for these earthquakes, consistent with dextral offset on the Wallula fault zone inferred from offset aeromagnetic anomalies associated with ca. 8.5 Ma basalt dikes. Magnetic surveys show that the Wallula fault actually lies 350 m to the southwest of the trace shown on published maps, passes directly through deformed late Pleistocene or younger deposits exposed at Finley quarry, and extends uninterrupted over 120 km.

Washington↗

USGS2016

Linking silicate weathering to riverine geochemistry—A case study from a mountainous tropical setting in west-central Panama

Chemical analyses from 71 watersheds across an ∼450 km transect in west-central Panama provide insight into controls on weathering and rates of chemical denudation and CO 2 consumption across an igneous arc terrain in the tropics. Stream and river compositions across this region of Panama are generally dilute, having a total dissolved solute value = 118 ± 91 mg/L, with bicarbonate and silica being the predominant dissolved species. Solute, stable isotope, and radiogenic isotope compositions are consistent with dissolution of igneous rocks present in Panama by meteoric precipitation, with geochemical signatures of rivers largely acquired in their upstream regions. Comparison of a headwater basin with its entire watershed observed considerably more runoff production from the high-elevation upstream portion of the catchment than in its much more spatially extensive downstream region. Rock alteration profiles document that weathering proceeds primarily by dissolution of feldspar and pyroxene, with base cations effectively leached in the following sequence: Na > Ca > Mg > K. Control on water chemistry by bedrock lithology is indicated through a linking of elevated ([Na + K]/[Ca + Mg]) ratios in waters to a high proportion of catchment area silicic bedrock and low ratios to mafic bedrock. Sr-isotope ratios are dominated by basement-derived Sr, with only very minor, if any, contribution from other sources. Cation weathering of Ca sil + Mg sil + Na + K spans about an order in magnitude, from 3 to 32 tons/km 2 /yr. Strong positive correlations of chemical denudation and CO 2 consumption are observed with precipitation, mean watershed elevation, extent of land surface forest cover, and physical erosion rate.

Chagras and Pacora watersheds↗

USGS2017

Rare earth mineral potential in the southeastern U.S. Coastal Plain from integrated geophysical, geochemical, and geological approaches

We combined geophysical, geochemical, mineralogical, and geological data to evaluate the regional presence of rare earth element (REE)−bearing minerals in heavy mineral sand deposits of the southeastern U.S. Coastal Plain. We also analyzed regional differences in these data to determine probable sedimentary provenance. Analyses of heavy mineral separates covering the region show strong correlations between thorium, monazite, and xenotime, suggesting that radiometric equivalent thorium (eTh) can be used as a geophysical proxy for those REE-bearing minerals. Airborne radiometric data collected during the National Uranium Resource Evaluation (NURE) program cover the southeastern United States with line spacing varying from ∼2 to 10 km. These data show eTh highs over Cretaceous and Tertiary Coastal Plain sediments from the Cape Fear arch in North Carolina to eastern Alabama; these highs decrease with distance from the Piedmont. Quaternary sediments along the modern coasts show weaker eTh anomalies, except near coast-parallel ridges from South Carolina to northern Florida. Prominent eTh anomalies are also observed over large riverbeds and their floodplains, even north of the Cape Fear arch where surrounding areas are relatively low. These variations were verified using ground geophysical measurements and sample analyses, indicating that radiometric methods are a useful exploration tool at varying scales. Further analyses of heavy mineral separates showed regional differences, not only in concentrations of monazite, but also of rutile and staurolite, and in magnetic susceptibility. The combined properties suggest the presence of subregions where heavy mineral sediments are primarily sourced from high-grade metamorphic, low-grade metamorphic, or igneous terrains, or where they represent a mixing of these sources. Comparisons between interpreted sources of heavy mineral sands near the Fall Line and igneous and metamorphic Piedmont and Blue Ridge units showed a strong correspondence with rocks closest to the Fall Line and poor correspondence with rocks farther inland. This strongly suggests that the primary source of those heavy minerals, especially monazite, is the rocks that formed the rocky coast that was present during opening of the Atlantic Ocean, which in turn indicates the importance of coastal processes in forming heavy mineral sand concentrations. Furthermore, narrow radiometric eTh and K anomalies are associated with major rivers, indicating limited spatial influence of fluvial processes. Later coastal plain sediment deposition appears to have involved reworking of sediments, providing an “inheritance” of the rocky coast composition that persists for some distance from the Fall Line. However, this inheritance is reduced with distance, and sediments within ∼100 km of the coast in Georgia and Florida exhibit properties indicative of mixing from multiple sources.

GSA Bulletin↗

USGS2017

The Valmy thrust sheet: A regional structure formed during the protracted assembly of the Roberts Mountains allochthon, Nevada, USA

A synthesis of field, biostratigraphic, detrital zircon geochronologic, and remote sensing data across north-central Nevada, United States, defines a thick, regionally extensive sheet of Middle–Upper Ordovician Valmy Formation quartzite that structurally overlies deformed early Paleozoic units of the Roberts Mountains allochthon. Late Paleozoic regional unconformities that record tectonic disruptions have been recognized in the foreland of central and eastern Nevada and locally within the Roberts Mountains allochthon; these identify multiple, regional tectonic events between the Devonian–Mississippian initiation of the Antler orogeny and the Permian–Triassic Sonoma orogeny. However, few studies have documented the regional kinematic history of the Robert Mountains allochthon sensu stricto. In the Independence Mountains of northern Nevada, emplacement of the Roberts Mountains allochthon is restricted to the Mississippian. In the Tuscarora Mountains, the range west and southwest of the Independence Mountains, several deformation events have been identified, and emplacement of the thrust sheet containing the Valmy Formation is restricted to the Late Pennsylvanian–Early Permian. These structural and temporal relations, reflected in the Antler foreland basin adjacent to the Roberts Mountains allochthon and overlap sequences, suggest that the Roberts Mountains allochthon is a composite stratigraphic terrane assembled along the Cordilleran margin during two or more late Paleozoic contractional events. Valmy Formation deposits likely represent the development of coalescing submarine fans below or within bypass channels in a deep slope or rise environment. Petrographic characteristics, biostratigraphy, and detrital zircon U-Pb age populations of the Valmy Formation link it to coeval slope and rise turbidites of the Vinini Formation and shelfal Eureka Quartzite; Valmy Formation detrital zircon age populations are dissimilar to the rift-to-drift facies of the Neoproterozoic–Cambrian Prospect Mountain Quartzite. Throughout north-central Nevada, the Valmy Formation is in fault contact with units of the Roberts Mountains allochthon, including the Devonian–Mississippian Slaven Chert, Silurian–Devonian Elder Sandstone, and Cambrian(?)–Ordovician Vinini Formation, which were deformed prior to, or during, emplacement of the thrust sheet containing Valmy Formation quartzite. Our mapping and data synthesis, guided by regional quartz maps based on remote sensing (Advanced Spaceborne Thermal Emission and Reflection Radiometer [ASTER]) data, delineate similar structural relationships discontinuously for >200 km along strike of the Roberts Mountains allochthon. Exploration for concealed gold deposits within reach of drilling requires knowledge of the relative thicknesses of the Roberts Mountains allochthon and the Valmy Formation. Overall thicknesses of deformed Roberts Mountains allochthon units between the Valmy Formation and underlying carbonate rocks, which host large, world-class Carlin-type gold deposits, vary by hundreds of meters, but are generally less than 700 m in three of the areas studied here. Recognition of windows through and klippen of the Roberts Mountains allochthon is essential for identification of areas where deposits may be at or near the surface. Correspondingly, most ongoing exploration for Carlin-type gold deposits subjacent to the Roberts Mountains allochthon targets concealed deposits. The model proposed in this study is applicable to determining depth to rocks prospective for undiscovered deposits.

Nevada↗

USGS2017

High-resolution seismic profiling reveals faulting associated with the 1934 Ms 6.6 Hansel Valley earthquake (Utah, USA)

The 1934 Ms 6.6 Hansel Valley, Utah, earthquake produced an 8-km-long by 3-km-wide zone of north-south−trending surface deformation in an extensional basin within the easternmost Basin and Range Province. Less than 0.5 m of purely vertical displacement was measured at the surface, although seismologic data suggest mostly strike-slip faulting at depth. Characterization of the origin and kinematics of faulting in the Hansel Valley earthquake is important to understand how complex fault ruptures accommodate regions of continental extension and transtension. Here, we address three questions: (1) How does the 1934 surface rupture compare with faults in the subsurface? (2) Are the 1934 fault scarps tectonic or secondary features? (3) Did the 1934 earthquake have components of both strike-slip and dip-slip motion? To address these questions, we acquired a 6.6-km-long, high-resolution seismic profile across Hansel Valley, including the 1934 ruptures. We observed numerous east- and west-dipping normal faults that dip 40°−70° and offset late Quaternary strata from within a few tens of meters of the surface down to a depth of ∼1 km. Spatial correspondence between the 1934 surface ruptures and subsurface faults suggests that ruptures associated with the earthquake are of tectonic origin. Our data clearly show complex basin faulting that is most consistent with transtensional tectonics. Although the kinematics of the 1934 earthquake remain underconstrained, we interpret the disagreement between surface (normal) and subsurface (strike-slip) kinematics as due to slip partitioning during fault propagation and to the effect of preexisting structural complexities. We infer that the 1934 earthquake occurred along an ∼3-km wide, off-fault damage zone characterized by distributed deformation along small-displacement faults that may be alternatively activated during different earthquake episodes.

Utah↗

USGS2016

Rapid fluvial incision of a late Holocene lava flow: Insights from LiDAR, alluvial stratigraphy, and numerical modeling

Volcanic eruptions fundamentally alter landscapes, paving over channels, decimating biota, and emplacing fresh, unweathered material. The fluvial incision of blocky lava flows is a geomorphic puzzle. First, high surface permeability and lack of sediment should preclude geomorphically effective surface runoff and dissection. Furthermore, past work has demonstrated the importance of extreme floods in driving incision via column toppling and plucking in columnar basalt, but it is unclear how incision occurs in systems where surface blocks are readily mobile. We examine rapid fluvial incision of the Collier lava flow, an andesitic Holocene lava flow in the High Cascades of Oregon. Since lava flow emplacement ∼1600 yr ago, White Branch Creek has incised bedrock gorges up to 8 m deep into the coherent core of the lava flow and deposited >0.2 km 3 of sediment on the lava flow surface. Field observation points to a bimodal discharge regime in the channel, with evidence for both annual snowmelt runoff and outburst floods from Collier glacier, as well as historical evidence of vigorous glacial meltwater. To determine the range of discharge events capable of incision in White Branch Creek, we used a mechanistic model of fluvial abrasion. We show that the observed incision implies that moderate flows are capable of both initiating channel formation and sustaining incision. Our results have implications for the evolution of volcanic systems worldwide, where glaciation and/or mass wasting may accelerate fluvial processes by providing large amounts of sediment to otherwise porous, sediment-starved landscapes.

GSA Bulletin↗

USGS2017

Luminescence dating of paleolake deltas and glacial deposits in Garwood Valley, Antarctica: Implications for climate, Ross ice sheet dynamics, and paleolake duration

The formation of perched deltas and other lacustrine deposits in the McMurdo Dry Valleys of Antarctica is widely considered to be evidence of valley-filling lakes dammed by the grounded Ross Sea ice sheet during the local Last Glacial Maximum, with lake drainage interpreted as a record of grounding line retreat. We used luminescence dating to determine the age of paleolake deltas and glacial tills in Garwood Valley, a coastal dry valley that opens to the Ross Sea. Luminescence ages are stratigraphically consistent with radiocarbon results from algal mats within the same delta deposits but suggest radiocarbon dates from lacustrine carbonates may overestimate deposit ages by thousands of years. Results suggest that late Holocene delta deposition into paleolake Howard in Garwood Valley persisted until ca. 3.5 ka. This is significantly younger than the date when grounded ice is thought to have retreated from the Ross Sea. Our evidence suggests that the local, stranded ice-cored till topography in Garwood Valley, rather than regional ice-sheet dynamics, may have controlled lake levels for some McMurdo Dry Valleys paleolakes. Age control from the supraglacial Ross Sea drift suggests grounding and up-valley advance of the Ross Sea ice sheet into Garwood valley during marine oxygen isotope stage (MIS) 4 (71–78 ka) and the local Last Glacial Maximum (9–10 ka). This work demonstrates the power of combining luminescence dating with existing radiocarbon data sets to improve understanding of the relationships among paleolake formation, glacial position, and stream discharge in response to climate change.

GSA Bulletin↗

USGS2017

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.

California↗

USGS2017

New insight into the origin of manganese oxide ore deposits in the Appalachian Valley and Ridge of northeastern Tennessee and northern Virginia, USA

Manganese oxide deposits have long been observed in association with carbonates within the Appalachian Mountains, but their origin has remained enigmatic for well over a century. Ore deposits of Mn oxides from several productive sites located in eastern Tennessee and northern Virginia display morphologies that include botryoidal and branching forms, massive nodules, breccia matrix cements, and fracture fills. The primary ore minerals include hollandite, cryptomelane, and romanèchite. Samples of Mn oxides from multiple localities in these regions were analyzed using electron microscopy, X-ray analysis, Fourier transform infrared spectroscopy, and trace and rare earth element (REE) geochemistry. The samples from eastern Tennessee have biological morphologies, contain residual biopolymers, and exhibit REE signatures that suggest the ore formation was due to supergene enrichment (likely coupled with microbial activity). In contrast, several northern Virginia ores hosted within quartz-sandstone breccias exhibit petrographic relations, mineral morphologies, and REE signatures indicating inorganic precipitation, and a likely hydrothermal origin with supergene overprinting. Nodular accumulations of Mn oxides within weathered alluvial deposits that occur close to breccia-hosted Mn deposits in Virginia show geochemical signatures that are distinct from the breccia matrices and appear to reflect remobilization of earlier-emplaced Mn and concentration within supergene traps. Based on the proximity of all of the productive ore deposits to mapped faults or other zones of deformation, we suggest that the primary source of all of the Mn may have been deep seated, and that Mn oxides with supergene and/or biological characteristics resulted from the local remobilization and concentration of this primary Mn.

Tennessee, Virginia↗

USGS2018

Stress concentrations at structural discontinuities in active fault zones in the western United States: Implications for permeability and fluid flow in geothermal fields

Slip can induce concentration of stresses at discontinuities along fault systems. These structural discontinuities, i.e., fault terminations, fault step-overs, intersections, bends, and other fault interaction areas, are known to host fluid flow in ore deposition systems, oil and gas reservoirs, and geothermal systems. We modeled stress transfer associated with slip on faults with Holocene-to-historic slip histories at the Salt Wells and Bradys geothermal systems in western Nevada, United States. Results show discrete locations of stress perturbation within discontinuities along these fault systems. Well field data, surface geothermal manifestations, and subsurface temperature data, each a proxy for modern fluid circulation in the fields, indicate that geothermal fluid flow is focused in these same areas where stresses are most highly perturbed. These results suggest that submeter- to meter-scale slip on these fault systems generates stress perturbations that are sufficiently large to promote slip on an array of secondary structures spanning the footprint of the modern geothermal activity. Slip on these secondary faults and fractures generates permeability through kinematic deformation and allows for transmission of fluids. Still, mineralization is expected to seal permeability along faults and fractures over time scales that are generally shorter than either earthquake recurrence intervals or the estimated life span of geothermal fields. This suggests that though stress perturbations resulting from fault slip are broadly important for defining the location and spatial extent of enhanced permeability at structural discontinuities, continual generation and maintenance of flow conduits throughout these areas are probably dependent on the deformation mechanism(s) affecting individual structures.

GSA Bulletin↗

USGS2018

Contemporary fluvial geomorphology and suspended sediment budget of the partly confined, mixed bedrock-alluvial South River, Virginia, USA

We developed a conceptual model and suspended sediment budget for a 38 km reach of the fifth-order South River, Virginia, for the past 75 yr. Bedrock, terraces, and alluvial fans confine 64% of the channel’s lateral boundaries, while bedrock exposures impose vertical confinement along 37% of the channel. Bedrock exposures in the bed separate pools and riffles developed in gravelly bed material, create unusual kilometer-long pools, and divide the study area into a gently sloping upstream reach and a steeply sloping downstream reach. Bedrock exposures upstream and downstream of an alluvial monitoring site limit changes in bed elevation (documented by scour chains and repeat surveys) by flows with up to 10 yr return periods. Fifty-seven islands (features rarely mentioned in previous studies), mostly created by avulsive floodplain incision, occur in the study reach. Rates of bank retreat, likely moderated by bedrock exposures, have modal values of only a few centimeters per year, while floodplain growth by lateral accretion is negligible. Overbank deposition dominates the sediment budget, but the areal of the extent of the floodplain is currently being reduced by bank erosion and channel widening. The South River stores 2.5% of its annual suspended sediment load per kilometer of downstream transport, demonstrating that suspended sediment storage along partly confined, mixed bedrock-alluvial rivers can be equivalent to storage along fully alluvial rivers. The future evolution of the South River will likely be controlled by bank stabilization designed to control mercury loading into the channel from erosion of contaminated floodplain sediments.

Virginia↗