ABSTRACT This paper summarizes the hydrological variability in eastern California (central Sierra Nevada) for the past 3000 yr based on three distinct paleoclimate proxies, δ 18 O, total inorganic carbon (TIC), and magnetic susceptibility (chi). These proxies, which are recorded in lake sediments of Pyramid Lake and Walker Lake, Nevada, and Mono Lake and Owens Lake, California, indicate lake-level changes that are mostly due to variations in Sierra Nevada snowpack and rainfall. We evaluated lake-level changes in the four Great Basin lake systems with regard to sediment-core locations and lake-basin morphologies, to the extent that these two factors influence the paleoclimate proxy records. We documented the strengths and weaknesses of each proxy and argue that a systematic study of all three proxies together significantly enhances our ability to characterize the regional pattern, chronology, and resolution of hydrological variability. We used paleomagnetic secular variation (PSV) to develop paleomagnetic chronostratigraphies for all four lakes. We previously published PSV records for three of the lakes (Mono, Owens, Pyramid) and developed a new PSV record herein for Walker Lake. We show that our PSV chronostratigraphies are almost identical to previously established radiocarbon-based chronologies, but that there are differences of 20–200 yr in individual age records. In addition, we used eight of the PSV inclination features to provide isochrons that permit exacting correlations between lake records. We also evaluated the temporal resolution of our proxies. Most can document decadal-scale variability over the past 1000 yr, multidecadal-scale variability for the past 2000 yr, and centennial-scale variability between 2000 and 3000 yr ago. Comparisons among our proxies show a strong coherence in the pattern of lake-level variability for all four lakes. Pyramid Lake and Walker Lake have the longest and highest-resolution records. The δ 18 O and TIC records yield the same pattern of lake-level variability; however, TIC may allow a somewhat higher-frequency resolution. It is not clear, however, which proxy best estimates the absolute amplitude of lake-level variability. Chi is the only available proxy that records lake-level variability in all four lakes prior to 2000 yr ago, and it shows consistent evidence of a large multicentennial period of drought. TIC, chi, and δ 18 O are integrative proxies in that they display the cumulative record of hydrologic variability in each lake basin. Tree-ring estimations of hydrological variability, by contrast, are incremental proxies that estimate annual variability. We compared our integrated proxies with tree-ring incremental proxies and found a strong correspondence among the two groups of proxies if the tree-ring proxies are smoothed to decadal or multidecadal averages. Together, these results indicate a common pattern of wet/dry variability in California (Sierra Nevada snowpack/rainfall) extending from a few years (notable only in the tree-ring data) to perhaps 1000 yr. Notable hydrologic variability has occurred at all time scales and should continue into the future.

Surface, unsaturated-zone, and saturated-zone hydrologic conditions at Yucca Mountain responded to past climate variations and are at least partly preserved by sediment, fossil, and mineral records. Characterizing past hydrologic conditions in surface and subsurface environments helps to constrain hydrologic responses expected under future climate conditions and improve predictions of repository performance. Furthermore, these records provide a better understanding of hydrologic processes that operate at time scales not readily measured by other means. Pleistocene climates in southern Nevada were predominantly wetter and colder than the current interglacial period. Cyclic episodes of aggradation and incision in Fortymile Wash, which drains the eastern slope of Yucca Mountain, are closely linked to Pleistocene climate cycles. Formation of pedogenic cement is favored under wetter Pleistocene climates, consistent with increased soil moisture and vegetation, higher chemical solubility, and greater evapotranspiration relative to Holocene soil conditions. The distribution and geochemistry of secondary minerals in subsurface fractures and cavities reflect unsaturated-zone hydrologic conditions and the response of the hydrogeologic system to changes in temperature and percolation flux over the last 12.8 m.y. Physical and fluid-inclusion evidence indicates that secondary calcite and opal formed in air-filled cavities from fluids percolating downward through connected fracture pathways in the unsaturated zone. Oxygen, strontium, and carbon isotope data from calcite are consistent with a descending meteoric water source but also indicate that water compositions and temperatures evolved through time. Geochronological data indicate that secondary mineral growth rates are less than 1–5 mm/m.y., and have remained approximately uniform over the last 10 m.y. or longer. These data are interpreted as evidence for hydrological stability despite large differences in surface moisture caused by climate shifts between the Miocene and Pleistocene and between Pleistocene glacial-interglacial cycles. Secondary mineral distribution and δ 18 O profiles indicate that evaporation in the shallower welded tuffs reduces infiltration fluxes. Several near-surface and subsurface processes likely are responsible for diverting or dampening infiltration and percolation, resulting in buffering of percolation fluxes to the deeper unsaturated zone. Cooler and wetter Pleistocene climates resulted in increased recharge in upland areas and higher water tables at Yucca Mountain and throughout the region. Discharge deposits in the Amargosa Desert were active during glacial periods, but only in areas where the modern water table is within 7–30 m of the surface. Published groundwater models simulate water-table rises beneath Yucca Mountain of as much as 150 m during glacial climates. However, most evidence from Fortymile Canyon up gradient from Yucca Mountain limits water-table rises to 30 m or less, which is consistent with evidence from discharge sites in the Amargosa Desert. The isotopic compositions of uranium in tuffs spanning the water table in two Yucca Mountain boreholes indicate that Pleistocene water-table rises likely were restricted to 25–50 m above modern positions and are in approximate agreement with water-table rises estimated from zeolitic-to-vitric transitions in the Yucca Mountain tuffs (less than 60 m in the last 11.6 m.y.).

During glacial (pluvial) climatic periods, Death Valley is hypothesized to have episodically been the terminus for the Amargosa, Owens, and Mojave Rivers. Geological and biological studies have tended to support this hypothesis and a hydrological link that included the Colorado River, allowing dispersal of pupfish throughout southeastern California and western Nevada. Recent mitochondrial deoxyribonucleic acid (mtDNA) studies show a common pupfish (Cyprinodontidae) ancestry in this region with divergence beginning 3–2 Ma. We present tephrochronologic and paleomagnetic data in the context of testing the paleohydrologic connections with respect to the common collection point of the Amargosa, Owens, and Mojave Rivers in Death Valley during successive time periods: (1) the late Pliocene to early Pleistocene (3–2 Ma), (2) early to middle Pleistocene (1.2–0.5 Ma), and (3) middle to late Pleistocene (<0.7–0.03 Ma; paleolakes Manly and Mojave). Using the 3.35 Ma Zabriskie Wash tuff and 3.28 Ma Nomlaki Tuff Member of the Tuscan and Tehama Formations, which are prominent marker beds in the region, we conclude that at 3–2 Ma, a narrow lake occupied the ancient Furnace Creek Basin and that Death Valley was not hydrologically connected with the Amargosa or Mojave Rivers. A paucity of data for Panamint Valley does not allow us to evaluate an Owens River connection to Death Valley ca. 3–2 Ma. Studies by others have shown that Death Valley was not hydrologically linked to the Amargosa, Owens, or Mojave Rivers from 1.2 to 0.5 Ma. We found no evidence that Lake Manly flooded back up the Mojave River to pluvial Lake Mojave between 0.18 and 0.12 Ma, although surface water flowed from the Amargosa and Owens Rivers to Death Valley at this time. There is also no evidence for a connection of the Owens, Amargosa, or Mojave Rivers to the Colorado River in the last 3–2 m.y. Therefore, the hypothesis that pupfish dispersed or were isolated in basins throughout southeastern California and western Nevada by such a connection is not supported. Beyond the biologically predicted time frame, however, sparse and disputed data suggest that a fluvial system connected Panamint (Owens River), Death, and Amargosa Valleys, which could account for the dispersal and isolation before 3 Ma.

Pleistocene deposits in Panamint Valley, California, document the changes in pluvial lake level, source water, and elevation of the regional groundwater table associated with climate change. The oxygen isotope stage (OIS) 2 and 6 lacustrine record is well preserved in surficial deposits, whereas the OIS 3–5 lacustrine-paludal and lacustrine record is mainly derived from an archived core sample. Amino acid racemization ratios in ostracodes and gastropods suggest that the shoreline and groundwater-discharge features that lie between ∼600 and 550 m elevation formed during one highstand, probably during OIS 6. A fossiliferous part of the ∼100-m-deep core DH-1, which was drilled in the Ballarat Basin during the late 1950s, was resampled in this study. Comparison of DH-1 with core DH-3 from Panamint Valley and core OL-92 from Owens Lake suggests the 34–78-m-depth interval of DH-1 may span all or much of OIS 4. The microfauna from this depth interval indicate a saline marsh or shallow lacustrine environment, but not a large lake. The ostracode assemblage requires low ratios of alkalinity to calcium (alk/Ca) water likely indicative of solutes in deep regional groundwater sources rather than the high alk/Ca solutes common to the Owens River system. OIS 2–aged sediment from surficial deposits, a shallow auger hole, and core DH-1 contain faunas, including the ostracode Limnocythere sappaensis , which require the high alk/Ca evolved solutes common to the Owens River. The elevation of the lacustrine sediments further indicates a moderate-sized saline lake around 180–200 m depth. In the northern Lake Hill basin, a saline lake persisted until at least 16 ka, and it was succeeded by fresh, groundwater-supported wetlands, which were fully developed by ca. 12,575 14 C yr B.P. and which persisted until around 10,500 14 C yr B.P., when the basin became a dry playa.

The chronologic history of pluvial Owens Lake along the eastern Sierra Nevada in Owens Valley, California, has previously been reported for the interval of time from ca. 25 calibrated ka to the present. However, the age, distribution, and paleoclimatic context of higher-elevation shoreline features have not been formally documented. We describe the location and characteristics of wave-formed erosional and depositional features, as well as fluvial strath terraces that grade into an older shoreline of pluvial Owens Lake. These pluvial-lacustrine features are described between the Olancha area to the south and Poverty Hills area to the north, and they appear to be vertically deformed ∼20 ± 4 m across the active oblique-dextral Owens Valley fault zone. They occur at elevations from 1176 to 1182 m along the lower flanks of the Inyo Mountains and Coso Range east of the fault zone to as high as ∼1204 m west of the fault zone. This relict shoreline, referred to as the 1180 m shoreline, lies ∼20–40 m higher than the previously documented Last Glacial Maximum shoreline at ∼1160 m, which occupied the valley during marine isotope stage 2 (MIS 2). Crosscutting relations of wave-formed platforms, notches, and sandy beach deposits, as well as strath terraces on lava flows of the Big Pine volcanic field, bracket the age of the 1180 m shoreline to the time interval between ca. 340 ∼ 60 ka and ca. 130 ∼ 50 ka. This interval includes marine oxygen isotope stages 8–6 (MIS 8–6), corresponding to 260–240 ka and 185–130 ka, respectively. An additional age estimate for this shoreline is provided by a cosmogenic 36 Cl model age of ca. 160 ∼ 32 ka on reefal tufa at ∼1170 m elevation from the southeastern margin of the valley. This 36 Cl model age corroborates the constraining ages based on dated lava flows and refines the lake age to the MIS 6 interval. Documentation of this larger pluvial Owens Lake offers insight to the hydrologic balance along the east side of the southern Sierra Nevada and will assist with future regional paleoclimatic models within the western Basin and Range.

Owens Lake has existed for most of the past 800,000 yr, but the sequence of interconnected lakes and streams of which it was often part, the Owens River cascade, last flourished during late Pleistocene time. A fluctuating, increasingly saline, terminal lake survived into the late Holocene until upstream water diversions to the Los Angeles Aqueduct began in 1913. Shoreline fragments and beach stratigraphy indicate that the lake reached its highest late Pleistocene level around 23.5 ka, during the Last Glacial Maximum, when it was fed by meltwaters from Sierra Nevada glaciers and spilled southward to Searles Lake and beyond. The lake then fell to relatively low levels after 16.5 ka before experiencing terminal Pleistocene oscillations related to hydroclimatic forcing, which involved changing regional precipitation regimes rather than major inputs from Sierra Nevada glaciers. Two major transgressions occurred. The first culminated around 14.3 ka and was probably related to a cooler, wetter regional climate. The second culminated around 12.8 ka and was linked to the earlier wetter phase of the Younger Dryas cold event. However, the high late Pleistocene shoreline is deformed, and the highest beach ranges in elevation from 1140 m to 1167 m above sea level. If the terminal Pleistocene lake overflowed, as suggested here, then its outlet has also been raised since 12.8 ka. This deformation appears to have involved uplift of the Coso Range magmatic complex relative to subsidence and faulting within the Owens Lake graben between the Sierra Nevada and Inyo Mountains frontal faults. Such deformation confounds simple hydroclimatic explanations of lake behavior and must be incorporated into models that seek to interpret the changing form and geochemistry of Owens Lake and the frequency of its spillage southward to Searles Lake.