Major-ion chemistry, strontium isotope ratios ( 87 Sr/ 86 Sr), stable isotope ratios (δ 18 O, δ 2 H), and tritium were analyzed for water samples from the southern Bear Lake Valley, Utah and Idaho, to characterize the types and distribution of groundwater sources and their relation to Bear Lake’s pre-diversion chemistry. Four ground-water types were identified: (1) Ca-Mg-HCO 3 water with 87 Sr/ 86 Sr values of ~0.71050 and modern tritium concentrations was found in the mountainous carbonate terrain of the Bear River Range. Magnesium (Mg) and bicarbonate (HCO 3 ) concentrations at Swan Creek Spring are discharge dependent and result from differential carbonate bedrock dissolution within the Bear River Range. (2) Cl-rich groundwater with elevated barium and strontium concentrations and 87 Sr/ 86 Sr values between 0.71021 and 0.71322 was found in the southwestern part of the valley. This groundwater discharges at several small, fault-controlled springs along the margin of the lake and contains solutes derived from the Wasatch Formation. (3) SO 4 -rich groundwater with 87 Sr/ 86 Sr values of ~0.70865, and lacking detectable tritium, discharges from two springs in the northeast quadrant of the study area and along the East Bear Lake fault. (4) Ca-Mg-HCO 3 -SO 4 -Cl water with 87 Sr/ 86 Sr values of ~0.71060 and sub-modern tritium concentrations discharges from several small springs emanating from the Wasatch Formation on the Bear Lake Plateau. The δ 18 O and δ 2 H values from springs and streams discharging in the Bear River Range fall along the Global Meteoric Water Line (GMWL), but are more negative at the southern end of the valley and at lower elevations. The δ 18 O and δ 2 H values from springs discharging on the Bear Lake Plateau plot on an evaporation line slightly below the GMWL. Stable isotope data suggest that precipitation falling in Bear Lake Valley is affected by orographic effects as storms pass over the Bear River Range, and by evaporation prior to recharging the Bear Lake Plateau aquifers. Approximately 99% of the solutes constituting Bear Lake’s pre-diversion chemistry were derived from stream discharge and shallow groundwater sources located within the Bear River Range. Lake-marginal springs exposed during the recent low lake levels and springs and streams draining the Bear Lake Plateau did not contribute significantly to the pre-diversion chemistry of Bear Lake.

Bear Lake, Utah and Idaho, is one of only a few lakes worldwide with endemic ostracode species. In most lakes, ostracode species distributions vary systematically with depth, but in Bear Lake, there is a distinct boundary in the abundances of cosmopolitan and endemic valves in surface sediments at ~7 m water depth. This boundary seems to coincide with the depth distribution of endemic fish, indicating a biological rather than environmental control on ostracode species distributions. The cosmopolitan versus endemic ostracode species distribution persisted through time in Bear Lake and in a neighboring wetland. The endemic ostracode fauna in Bear Lake implies a complex ecosystem that evolved in a hydrologically stable, but not invariant, environmental setting that was long lived. Long-lived (geologic time scale) hydrologic stability implies the lake persisted for hundreds of thousands of years despite climate variability that likely involved times when effective moisture and lake levels were lower than today. The hydrologic budget of the lake is dominated by snowpack meltwater, as it likely was during past climates. The fractured and karstic bedrock in the Bear Lake catchment sustains local stream flow through the dry summer and sustains stream and ground-water flow to the lake during dry years, buffering the lake hydrology from climate variability and providing a stable environment for the evolution of endemic species.

A continuous, 120-m-long core (BL00-1) from Bear Lake, Utah and Idaho, contains evidence of hydrologic and environmental change over the last two glacial-interglacial cycles. The core was taken at 41.95°N, 111.31°W, near the depocenter of the 60-m-deep, spring-fed, alkaline lake, where carbonate-bearing sediment has accumulated continuously. Chronological control is poor but indicates an average sedimentation rate of 0.54 mm yr ‒1 . Analyses have been completed at multi-centennial to millennial scales, including (in order of decreasing temporal resolution) sediment magnetic properties, oxygen and carbon isotopes on bulk-sediment carbonate, organic- and inorganic- carbon contents, palynology; mineralogy (X-ray diffraction), strontium isotopes on bulk carbonate, ostracode taxonomy, oxygen and carbon isotopes on ostracodes, and diatom assemblages. Massive silty clay and marl constitute most of the core, with variable carbonate content (average = 31 ± 19%) and oxygen-isotopic values (δ 18 O ranging from ‒18‰ to ‒5‰ in bulk carbonate). These variations, as well as fluctuations of biological indicators, reflect changes in the water and sediment discharged from the glaciated headwaters of the dominant tributary, Bear River, and the processes that influenced sediment delivery to the core site, including lake-level changes. Although its influence has varied, Bear River has remained a tributary to Bear Lake during most of the last quarter-million years. The lake disconnected from the river and, except for a few brief excursions, retracted into a topographically closed basin during global interglaciations (during parts of marine isotope stages 7, 5, and 1). These intervals contain up to 80% endogenic aragonite with high δ 18 O values (average = ‒5.8 ± 1.7‰), indicative of strongly evaporitic conditions. Interglacial intervals also are dominated by small, benthic/tychoplanktic fragilarioid species indicative of reduced habitat availability associated with low lake levels, and they contain increased high-desert shrub and Juniperus pollen and decreased forest and forest-woodland pollen. The 87 Sr/ 86 Sr values (>0.7100) also increase, and the ratio of quartz to dolomite decreases, as expected in the absence of Bear River inflow. The changing paleoenvironments inferred from BL00-1 generally are consistent with other regional and global records of glacial-interglacial fluctuations; the diversity of paleoenvironmental conditions inferred from BL00-1 also reflects the influence of catchment-scale processes.