Elmira Wan

Importance of groundwater in propagating downward integration of the 6–5 Ma Colorado River system: Geochemistry of springs, travertines, and lacustrine carbonates of the Grand Canyon region over the past 12 Ma

We applied multiple geochemical tracers ( 87 Sr/ 86 Sr, [Sr], δ 13 C, and δ 18 O) to waters and carbonates of the lower Colorado River system to evaluate its paleohydrology over the past 12 Ma. Modern springs in Grand Canyon reflect mixing of deeply derived (endogenic) fluids with meteoric (epigenic) recharge. Travertine ( 87 Sr/ 86 Sr and δ 13 C and δ 18 O values that overlap with associated water values, providing justification for use of carbonates as a proxy for the waters from which they were deposited. The Hualapai Limestone (12–6 Ma) and Bouse Formation (5.6–4.8 Ma) record paleohydrology immediately prior to and during integration of the Colorado River. The Hualapai Limestone was deposited from 12 Ma (new ash age) to 6 Ma; carbonates thicken eastward to ∼210 m toward the Grand Wash fault, suggesting that deposition was synchronous with fault slip. A fanning-dip geometry is suggested by correlation of ashes between subbasins using tephrochronology. New detrital-zircon ages are consistent with the “Muddy Creek constraint, ” which posits that Grand Wash Trough was internally drained prior to 6 Ma, with limited or no Colorado Plateau detritus, and that Grand Wash basin was sedimentologically distinct from Gregg and Temple basins until after 6 Ma. New isotopic data from Hualapai Limestone of Grand Wash basin show values and ranges of 87 Sr/ 86 Sr, δ 13 C, and δ 18 O that are similar to Grand Canyon springs and travertines, suggesting a long-lived springfed lake/marsh system sourced from western Colorado Plateau groundwater. Progressive up-section decrease in 87 Sr/ 86 Sr and δ 13 C and increase in δ 18 O in the uppermost 50 m of the Hualapai Limestone indicate an increase in meteoric water relative to endogenic inputs, which we interpret to record progressively increased input of high-elevation Colorado Plateau groundwater from ca. 8 to 6 Ma. Grand Wash, Hualapai, Gregg, and Temple basins, although potentially connected by groundwater, were hydrochemically distinct basins before ca. 6 Ma. The 87 Sr/ 86 Sr, δ 13 C, and δ 18 O chemostratigraphic trends are compatible with a model for downward integration of Hualapai basins by groundwater sapping and lake spillover. The Bouse Limestone (5.6–4.8 Ma) was also deposited in several hydrochemically distinct basins separated by bedrock divides. Northern Bouse basins (Cottonwood, Mojave, Havasu) have carbonate chemistry that is nonmarine. The 87 Sr/ 86 Sr data suggest that water in these basins was derived from mixing of high- 87 Sr/ 86 Sr Lake Hualapai waters with lower- 87 Sr/ 86 Sr, first-arriving “Colorado River” waters. Covariation trends of δ 13 C and δ 18 O suggest that newly integrated Grand Wash, Gregg, and Temple basin waters were integrated downward to the Cottonwood and Mojave basins at ca. 5–6 Ma. Southern, potentially younger Bouse basins are distinct hydrochemically from each other, which suggests incomplete mixing during continued downward integration of internally drained basins. Bouse carbonates display a southward trend toward less radiogenic 87 Sr/ 86 Sr values, higher [Sr], and heavier δ 18 O that we attribute to an increased proportion of Colorado River water through time plus increased evaporation from north to south. The δ 13 C and δ 18 O trends suggest alternating closed and open systems in progressively lower (southern) basins. We interpret existing data to permit the interpretation that the southernmost Blythe basin may have had intermittent mixing with marine water based on δ 13 C and δ 18 O covariation trends, sedimentology, and paleontology. [Sr] versus 87 Sr/ 86 Sr modeling suggests that southern Blythe basin 87 Sr/ 86 Sr values of ∼0.710–0.711 could be produced by 25%–75% seawater mixed with river water (depending on [Sr] assumptions) in a delta– marine estuary system. We suggest several refinements to the “lake fill-and-spill” downward integration model for the Colorado River: (1) Lake Hualapai was fed by western Colorado Plateau groundwater from 12 to 8 Ma; (2) high-elevation Colorado Plateau groundwater was progressively introduced to Lake Hualapai from ca. 8 to 6 Ma; (3) Colorado River water arrived at ca. 5–6 Ma; and (4) the combined inputs led to downward integration by a combination of groundwater sapping and sequential lake spillover that first delivered Colorado Plateau water and detritus to the Salton Trough at ca. 5.3 Ma. We propose that the groundwater sapping mechanism strongly influenced lake evolution of the Hualapai and Bouse Limestones and that groundwater flow from the Colorado Plateau to Grand Wash Trough led to Colorado River integration.

Review and analysis of the age and origin of the Pliocene Bouse Formation, lower Colorado River Valley, southwestern USA

The lower Pliocene Bouse Formation in the lower Colorado River Valley (southwestern USA) consists of basal marl and dense tufa overlain by siltstone and fine sandstone. It is locally overlain by and interbedded with sands derived from the Colorado River. We briefly review 87 Sr/ 86 Sr analyses of Bouse carbonates and shells and carbonate and gypsum of similar age east of Las Vegas that indicate that all of these strata are isotopically similar to modern Colorado River water. We also review and add new data that are consistent with a step in Bouse Formation maximum elevations from 330 m south of Topock Gorge to 555 m to the north. New geochemical data from glass shards in a volcanic ash bed within the Bouse Formation, and from an ash bed within similar deposits in Bristol Basin west of the Colorado River Valley, indicate correlation of the two ash beds and coeval submergence of both areas. The tuff bed is identified as the 4.83 Ma Lawlor Tuff derived from the San Francisco Bay region. We conclude, as have some others, that the Bouse Formation was deposited in lakes produced by first-arriving Colorado River water that entered closed basins inherited from Basin and Range extension, and estimate that first arrival of river water occurred ca. 4.9 Ma. If this interpretation is correct, addition of Bristol Basin to the Blythe Basin inundation area means that river discharge was sufficient to fill and spill a lake with an area of ∼10,000 km 2 . For spillover to occur, evaporation rates must have been significantly less in early Pliocene time than modern rates of ∼2–4 m/yr, and/or Colorado River discharge was significantly greater than the current ∼15 km 3 /yr. In this lacustrine interpretation, evaporation rates were sufficient to concentrate salts to levels that were hospitable to some marine organisms presumably introduced by birds.

Age, composition, and areal distribution of the Pliocene Lawlor Tuff, and three younger Pliocene tuffs, California and Nevada

The Lawlor Tuff is a widespread dacitic tephra layer produced by Plinian eruptions and ash flows derived from the Sonoma Volcanics, a volcanic area north of San Francisco Bay in the central Coast Ranges of California, USA. The younger, chemically similar Huichica tuff, the tuff of Napa, and the tuff of Monticello Road sequentially overlie the Lawlor Tuff, and were erupted from the same volcanic field. We obtain new laser-fusion and incremental-heating 40 Ar/ 39 Ar isochron and plateau ages of 4.834 ± 0.011, 4.76 ± 0.03, ≤4.70 ± 0.03, and 4.50 ± 0.02 Ma (1 sigma), respectively, for these layers. The ages are concordant with their stratigraphic positions and are significantly older than those determined previously by the K-Ar method on the same tuffs in previous studies. Based on offsets of the ash-flow phase of the Lawlor Tuff by strands of the eastern San Andreas fault system within the northeastern San Francisco Bay area, total offset east of the Rodgers Creek–Healdsburg fault is estimated to be in the range of 36 to 56 km, with corresponding displacement rates between 8.4 and 11.6 mm/yr over the past ∼4.83 Ma. We identify these tuffs by their chemical, petrographic, and magnetic characteristics over a large area in California and western Nevada, and at a number of new localities. They are thus unique chronostratigraphic markers that allow correlation of marine and terrestrial sedimentary and volcanic strata of early Pliocene age for their region of fallout. The tuff of Monticello Road is identified only near its eruptive source.

Incision history of the Black Canyon of Gunnison, Colorado, over the past ∼1 Ma inferred from dating of fluvial gravel deposits

Spatio-temporal variability in fluvial incision rates in bedrock channels provides data regarding uplift and denudation histories of landscapes. The longitudinal profile of the Gunnison River (Colorado), tributary to the Colorado River, contains a prominent knickzone with 800 m of relief across it within the Black Canyon of the Gunnison. Average bedrock incision rates over the last 0.64 Ma surrounding the knickpoint vary from 150 m/Ma (downstream) to 400–550 m/Ma (within) to 90–95 m/Ma (upstream), suggesting it is a transient feature. Lava Creek B ash constrains strath terraces along a paleoprofile of the river. An isochron cosmogenic burial date in the paleo–Bostwick River of 870 ± 220 ka is consistent with the presence of 0.64 Ma Lava Creek B ash in locally derived, stratigraphically younger sediment. With 350 m of incision since deposition, we determine an incision rate of 400–550 m/Ma, reflecting incision through resistant basement rock at 2–3 times regional incision rates. Such contrast is attributed to a wave of transient incision, potentially initiated by downstream base-level fall during abandonment of Unaweep Canyon at ca. 1 Ma. Rate extrapolation indicates that the ∼700 m depth of Black Canyon has been eroded since 1.3–1.75 Ma. The Black Canyon knickpoint overlies a strong gradient between low-velocity mantle under the Colorado Rockies and higher-velocity mantle of the Colorado Plateau. We interpret recent reorganization and transient incision of both the Gunnison River and upper Colorado River systems to be a response to mantle-driven epeirogenic uplift of the southern Rockies in the last 10 Ma.

The story of a Yakima fold and how it informs late Neogene and Quaternary backarc deformation in the Cascadia subduction zone, Manastash anticline, Washington, USA

The Yakima folds of central Washington, USA are prominent anticlines that are the primary tectonic features of the backarc of the northern Cascadia subduction zone. What accounts for their topographic expression and how much strain do they accommodate and over what time period? We investigate Manastash anticline, a north-vergent fault propagation fold typical of structures in the fold province. From retrodeformation of line- and area-balanced cross sections, the crust has horizontally shortened by 11% (0.8-0.9 km). The fold, and by inference all other folds in the fold province, formed no earlier than 15.6 Ma as they developed on a landscape that was reset to negligible relief following voluminous outpouring of Grande Ronde Basalt. Deformation is accommodated on two fault sets including west-northwest-striking frontal thrust faults and shorter north- to northeast-striking faults. The frontal thrust fault system is active with late Quaternary scarps at the base of the range front. The fault-cored Manastash anticline terminates to the east at the Naneum anticline and fault; activity on the north-trending Naneum structures predates emplacement of the Grande Ronde Basalt. The west-trending Yakima folds and west-striking thrust faults, the shorter north to northeast striking faults, and the Naneum fault together constitute the tectonic structures that accommodate deformation in the low-strain-rate environment in the backarc of the Cascadia Subduction Zone.

Late Neogene–Quaternary tephrochronology, stratigraphy, and paleoclimate of Death Valley, California, USA

Sedimentary deposits in midlatitude continental basins often preserve a paleoclimate record complementary to marine-based records. However, deriving that paleoclimate record depends on having well-exposed deposits and establishing a sufficiently robust geochronology. After decades of research, we have been able to correlate 77 tephra beds exposed in multiple stratigraphic sections in the Death Valley area, California, United States. These correlations identify 25 different tephra beds that erupted from at least five different volcanic centers from older than 3.58 Ma to ca. 32 ka. We have informally named and determined the ages for seven previously unrecognized beds: ca. 3.54 Ma tuff of Curry canyon, ca. 3.45 Ma tuff of Furnace Creek, ca. 3.1 Ma tuff of Kit Fox Hills, ca. 3.1 Ma tuff of Mesquite Flat, ca. 3.15 Ma tuff of Texas Spring, 3.117 ± 0.011 Ma tuff of Echo Canyon, and the ca. 1.3 Ma Amargosa ash bed. Several of these tephra beds are found as far northeast as central Utah and could be important marker beds in western North America. Our tephrochronologic data, combined with magnetic polarity data and ^(40)Ar/^(39)Ar age determinations, redefine Neogene sedimentary deposits exposed across 175 km^2 of the Death Valley area. The alluvial/lacustrine Furnace Creek Formation is a time-transgressive sedimentary sequence ranging from ca. 6.0 to 2.5 Ma in age. The ca. 2.5−1.7 Ma Funeral Formation is typically exposed as a proximal alluvial-fan facies overlying the Furnace Creek Formation. We have correlated deposits in the Kit Fox Hills, Salt Creek, Nova Basin, and southern Death Valley with the informally named ca. 1.3−0.5 Ma Mormon Point formation. In addition, our correlation of the late Pleistocene Wilson Creek ash bed 15 in the Lake Rogers deposits represents the first unambiguous sequences deposited during the Last Glacial Maximum (marine isotope stage [MIS] 2) in Death Valley. Based on this new stratigraphic framework, we show that the Pliocene and Pleistocene climate in Death Valley is consistent with the well-established marine tropical/subtropical record. Pluvial lakes in Death Valley and Searles Valley began to form ca. 3.5−3.4 Ma in the late Pliocene during MIS MG5. Initiation of lakes in these two hydrologically separated valleys at the same time at the beginning of a cooling trend in the marine climate record suggests a link to a cooler, wetter (glacial) regional climate in North America. The Death Valley lake persisted until ca. 3.30 Ma, at the peak of the M2 glaciation, after which there is no evidence of Pliocene lacustrine deposition, even at the peak of the Northern Hemisphere Glaciation (ca. 2.75 Ma). If pluvial lakes in the Pliocene are an indirect record of glacial climate conditions, as they are for the Pleistocene, then a glacial climate was present in western North America for ∼200,000 yr during the Pliocene, encompassing MIS MG5−M2. Pleistocene pluvial lakes in Death Valley that formed ca. 1.98−1.78 Ma, 1.3−1.0 Ma, and ca. 0.6 Ma (MIS 16) are consistent with other regional climate records that indicate a regional glacial climate; however, Death Valley was relatively dry at ca. 0.77 Ma (MIS 19), when large lakes existed in other basins. The limited extent of the MIS 2 marsh/shallow lake in the Lake Rogers basin of northern Death Valley reflects the well-known regional glacial climate at that time; however, Death Valley received relatively lower inflow and rainfall in comparison.

Geochronology of the Upper Alturas Formation, northern California: Implications for the Hemphillian-Blancan North American Land Mammal Age boundary

Fossil vertebrates from the Alturas Formation in northern California have previously been considered important for defining the age of the boundary between the Hemphillian and Blancan North American Land Mammal Ages. Diatomaceous mudstone of the upper Alturas Formation contain fossil mammals including the arvicoline rodent Mimomys (Ogmodontomys) sawrockensis that is diagnostic of Blancan faunas. New paleomagnetic and geochemical data from the upper Alturas Formation constrain the age of the first stratigraphic occurrence of M. (O.) sawrockensis at Crowder Flat Road to between 4.5 and 4.6 Ma. This age is approximately 0.2-0.4 Ma younger than previously reported such that the oldest record of Mimomys in North America, south of 55oN, is from Panaca, Nevada, and is constrained geochronologically to be approximately 4.9 Ma. The Hemphillian – Blancan North American Land Mammal Age boundary probably occurs within magnetic polarity Chron C3n.3r at approximately 4.9 Ma.