USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Holocene”

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

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

At least 1,261 records · Page 70Linked to original sources

Sediment dispersal in the northwestern Adriatic Sea

Sediment dispersal in the Adriatic Sea was evaluated using coupled three-dimensional circulation and sediment transport models, representing conditions from autumn 2002 through spring 2003. The calculations accounted for fluvial sources, resuspension by waves and currents, and suspended transport. Sediment fluxes peaked during southwestward Bora wind conditions that produced energetic waves and strengthened the Western Adriatic Coastal Current. Transport along the western Adriatic continental shelf was nearly always to the south, except during brief periods when northward Sirocco winds reduced the coastal current. Much of the modeled fluvial sediment deposition was near river mouths, such as the Po subaqueous delta. Nearly all Po sediment remained in the northern Adriatic. Material from rivers that drain the Apennine Mountains traveled farther before deposition than Po sediment, because it was modeled with a lower settling velocity. Fluvial sediment delivered to areas with high average bed shear stress was more highly dispersed than material delivered to more quiescent areas. Modeled depositional patterns were similar to observed patterns that have developed over longer timescales. Specifically, modeled Po sediment accumulation was thickest near the river mouth with a very thin deposit extending to the northeast, consistent with patterns of modern sediment texture in the northern Adriatic. Sediment resuspended from the bed and delivered by Apennine Rivers was preferentially deposited on the northern side of the Gargano Peninsula, in the location of thick Holocene accumulation. Deposition here was highest during Bora winds when convergences in current velocities and off-shelf flux enhanced delivery of material to the midshelf. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

Hydrogeologic controls on nitrate transport in a small agricultural catchment, Iowa

Effects of subsurface deposits on nitrate loss in stream riparian zones are recognized, but little attention has been focused on similar processes occurring in upland agricultural settings. In this paper, we evaluated hydrogeologic controls on nitrate transport processes occurring in a small 7.6 ha Iowa catchment. Subsurface deposits in the catchment consisted of upland areas of loess overlying weathered pre-Illinoian till, drained by two ephemeral drainageways that consisted of Holocene-age silty and organic rich alluvium. Water tables in upland areas fluctuated more than 4 m per year compared to less than 0.3 m in the drainageway. Water quality patterns showed a distinct spatial pattern, with groundwater in the drainageways having lower nitrate concentrations (<0.5 mg L-1 compared to upland areas (>10 mg L-1) as wells as lower pH, dissolved oxygen and redox, and higher ammonium and dissolved organic carbon levels. Several lines of evidence suggested that conditions are conducive for denitrification of groundwater flowing from uplands through the drainageways. Field-measured nitrate decay rates in the drainageways (???0.02 day-1) were consistent with other laboratory studies and regional patterns. Results from MODFLOW and MT3DMS simulations indicated that soils in the ephemeral drainageways could process all upland groundwater nitrate flowing through them. However, model-simulated tile drainage increased both water flux and nitrate loss from the upland catchment. Study results suggest that ephemeral drainageways can provide a natural nitrate treatment system in our upland glaciated catchments, offering management opportunities to reduce nitrate delivery to streams. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗

Quickly erupted volcanic sections of the Steens Basalt, Columbia River Basalt Group: Secular variation, tectonic rotation, and the Steens Mountain reversal

The Steens Basalt, now considered part of the Columbia River Basalt Group (CRBG), contains the earliest eruptions of this magmatic episode. Lava flows of the Steens Basalt cover about 50,000 km 2 of the Oregon Plateau in sections up to 1000 m thick. The large number of continuously exposed, quickly erupted lava flows (some sections contain over 200 flows) allows for small loops in the magnetic field direction paths to be detected. For volcanic rocks, this detail and fidelity are rarely found outside of the Holocene and yield estimates of eruption durations at our four sections of ∼2.5 ka for 260 m at Pueblo Mountains, 0.5 to 1.5 ka for 190 m at Summit Springs, 1–3 ka for 170 m at North Mickey, and ∼3 ka for 160 m at Guano Rim. That only one reversal of the geomagnetic field occurred during the eruption of the Steens Basalt (the Steens reversal at approximately 16.6 Ma) is supported by comparing 40 Ar/ 39 Ar ages and magnetic polarities to the geomagnetic polarity timescale. At Summit Springs two 40 Ar/ 39 Ar ages from normal polarity flows (16.72 ± ± 0.29 Ma (16.61) and 16.92 ± ± 0.52 Ma (16.82); ± ± equals 2 σ error) place their eruptions after the Steens reversal, while at Pueblo Mountains an 40 Ar/ 39 Ar age of 16.72 ± ± 0.21 Ma (16.61) from a reverse polarity flow places its eruption before the Steens reversal. Paleomagnetic field directions yielded 50 nontransitional directional-group poles which, combined with 26 from Steens Mountain, provide a paleomagnetic pole for the Oregon Plateau of 85.7°N, 318.4°E, K = 15.1, A 95 = 4.3. Comparison of this new pole with a reference pole derived from CRBG flows from eastern Washington and a synthetic reference pole for North America derived from global data implies relative clockwise rotation of the Oregon Plateau of 7.4 ± 5.0° or 14.5 ± 5.4°, respectively, probably due to northward decreasing extension of the basin and range.

Geochemistry, Geophysics, Geosystems↗

Was pre–twentieth century sea level stable?

Sea level rise (SLR) ranks high on the list of climate change issues because the expected acceleration from the current rate (about 3.1 millimeters per year) poses threats to coastal regions. Tide gauge, salt marsh, and archaeological records, and modeling of glacioisostatic adjustment (GIA) have led to the widely accepted idea that late Holocene (the past ∼2000 years) sea level was stable prior to acceleration beginning around 1850–1900 C.E. For instance, according to the Intergovernmental Panel on Climate Change Fourth Assessment Report, before the last century, sea level had “stabilized” over the past 2000 years, rising at a mean rate of 0–0.2 millimeter per year [ Bindoff et al. , 2007]. Others maintain that sea level was “nearly stable” over the past few thousand years [ Nicholls and Cazenave , 2010], pre–twentieth century rates were “close to zero” [ Church et al. , 2008], or “stable from at least BC 100 until AD 950” and “stable, or slightly falling” from 1350 until the nineteenth century [ Kemp et al. , 2011].

Eos, Transactions, American Geophysical Union↗

Peat accumulation in drained thermokarst lake basins in continuous, ice-rich permafrost, northern Seward Peninsula, Alaska

Thermokarst lakes and peat-accumulating drained lake basins cover a substantial portion of Arctic lowland landscapes, yet the role of thermokarst lake drainage and ensuing peat formation in landscape-scale carbon (C) budgets remains understudied. Here we use measurements of terrestrial peat thickness, bulk density, organic matter content, and basal radiocarbon age from permafrost cores, soil pits, and exposures in vegetated, drained lake basins to characterize regional lake drainage chronology, C accumulation rates, and the role of thermokarst-lake cycling in carbon dynamics throughout the Holocene on the northern Seward Peninsula, Alaska. Most detectable lake drainage events occurred within the last 4,000 years with the highest drainage frequency during the medieval climate anomaly. Peat accumulation rates were highest in young (50–500 years) drained lake basins (35.2 g C m −2 yr −1 ) and decreased exponentially with time since drainage to 9 g C m −2 yr −1 in the oldest basins. Spatial analyses of terrestrial peat depth, basal peat radiocarbon ages, basin geomorphology, and satellite-derived land surface properties (Normalized Difference Vegetation Index (NDVI); Minimum Noise Fraction (MNF)) from Landsat satellite data revealed significant relationships between peat thickness and mean basin NDVI or MNF. By upscaling observed relationships, we infer that drained thermokarst lake basins, covering 391 km 2 (76%) of the 515 km 2 study region, store 6.4–6.6 Tg organic C in drained lake basin terrestrial peat. Peat accumulation in drained lake basins likely serves to offset greenhouse gas release from thermokarst-impacted landscapes and should be incorporated in landscape-scale C budgets.

Alaska↗

Chronostratigraphic framework for the IODP Expedition 318 cores from the Wilkes Land Margin: Constraints for paleoceanographic reconstruction

The Integrated Ocean Drilling Program Expedition 318 to the Wilkes Land margin of Antarctica recovered a sedimentary succession ranging in age from lower Eocene to the Holocene. Excellent stratigraphic control is key to understanding the timing of paleoceanographic events through critical climate intervals. Drill sites recovered the lower and middle Eocene, nearly the entire Oligocene, the Miocene from about 17 Ma, the entire Pliocene and much of the Pleistocene. The paleomagnetic properties are generally suitable for magnetostratigraphic interpretation, with well-behaved demagnetization diagrams, uniform distribution of declinations, and a clear separation into two inclination modes. Although the sequences were discontinuously recovered with many gaps due to coring, and there are hiatuses from sedimentary and tectonic processes, the magnetostratigraphic patterns are in general readily interpretable. Our interpretations are integrated with the diatom, radiolarian, calcareous nannofossils and dinoflagellate cyst (dinocyst) biostratigraphy. The magnetostratigraphy significantly improves the resolution of the chronostratigraphy, particularly in intervals with poor biostratigraphic control. However, Southern Ocean records with reliable magnetostratigraphies are notably scarce, and the data reported here provide an opportunity for improved calibration of the biostratigraphic records. In particular, we provide a rare magnetostratigraphic calibration for dinocyst biostratigraphy in the Paleogene and a substantially improved diatom calibration for the Pliocene. This paper presents the stratigraphic framework for future paleoceanographic proxy records which are being developed for the Wilkes Land margin cores. It further provides tight constraints on the duration of regional hiatuses inferred from seismic surveys of the region.

Paleoceanography↗

Origin of primitive tholeiitic and calc-alkaline basalts at Newberry Volcano, Oregon

Major and trace element and radiogenic isotopic characteristics of primitive mafic Pleistocene and Holocene lavas from Newberry Volcano, Oregon, define two groups. The first consists of dry tholeiitic high-alumina olivine basalts that are slightly enriched in highly incompatible elements. The second group consists of calc-alkaline basalts that contained 2–4 wt % H 2 O prior to eruption and shows strong enrichment in the light rare earth elements, Ba, and Sr, and deficits in Nb, Ta, Hf, and Zr. The tholeiitic basalts reflect 6–11% anhydrous adiabatic decompression melting of spinel peridotite. The calc-alkaline basalts derived from compositionally distinct sources with strong LIL enrichment and relative depletion in HFSE, but with Sr, Nd, Hf, and Pb isotopic composition only slightly distinct from the sources of the tholeiitic magmas. Radiogenic Os correlates with LREE enrichment in the calc-alkaline magmas, which indicates that their source materials include a contribution from a mafic component that was melted in the garnet stability field. The calc-alkaline magmas were derived by melting of peridotite metasomatized by a fluid/melt that originated by melting of a mixture of the sediment plus MORB basalt/mantle in the underlying subducting oceanic plate. While the trace element characteristics of the calc-alkaline magmas were determined by the subduction component, their isotopic characteristics were modified during transit through the mantle by interaction with the highly magmatically processed mantle wedge beneath Newberry Volcano that, without the slab component, serves as the source of the tholeiitic magmas.

Oregon↗

Constraining the early eruptive history of the Mono Craters rhyolites, California, based on 238U–230Th isochron dating of their explosive and effusive products

The Mono Craters are an overlapping chain of at least 28 domes and coulees located south of Mono Lake, east central California, and represent the most recent eruptions of high‐silica rhyolite magma in the Mono Lake‐Long Valley volcanic region. Regionally widespread tephra fall deposits from the Mono Craters serve as important chronostratigraphic markers for correlations of late Quaternary terrestrial archives in the western United States. A well‐resolved eruption chronology for the Mono Craters is thus not only critical for volcanic hazard considerations but also for paleoclimatic and paleomagnetic studies. Here we constrain the timing of early Mono Craters volcanism using ion microprobe 238 U‐ 230 Th dating of unpolished crystal faces of autocrystic zircon and allanite microphenocrysts, which samples the final stage of crystallization prior to eruption. The stratigraphically oldest domes (biotite‐bearing rhyolites) yield 238 U‐ 230 Th isochron dates between ca. 42 and ca. 19 ka and are unambiguously linked to dated tephras in the Wilson Creek formation of ancestral Mono Lake via coupled 238 U‐ 230 Th geochronology and titanomagnetite geochemistry. The second oldest set of domes (orthopyroxene‐bearing rhyolites) have overlapping 238 U‐ 230 Th isochron dates that are within uncertainty of published K/Ar and 40 Ar/ 39 Ar dates for the fayalite‐bearing rhyolite domes, suggesting a period of intense effusive activity for the Mono Craters near the Pleistocene‐Holocene boundary. Our new and previously published 238 U‐ 230 Th dates for tephras in the Wilson Creek formation provide robust geochronologic constraints for the controversial geomagnetic excursion recorded in the Wilson Creek formation.

California↗

The mighty Susquehanna—extreme floods in Eastern North America during the past two millennia

The hazards posed by infrequent major floods to communities along the Susquehanna River and the ecological health of Chesapeake Bay remain largely unconstrained due to the short length of streamgage records. Here we develop a history of high‐flow events on the Susquehanna River during the late Holocene from flood deposits contained in MD99‐2209, a sediment core recovered in 26 m of water from Chesapeake Bay near Annapolis, Maryland, United States. We identify coarse‐grained deposits left by Hurricane Agnes (1972) and the Great Flood of 1936, as well as during three intervals that predate instrumental flood records (~1800–1500, 1300–1100, and 400–0 CE). Comparison to sedimentary proxy data (pollen and ostracode Mg/Ca ratios) from the same core site indicates that prehistoric flooding on the Susquehanna often accompanied cooler‐than‐usual winter/spring temperatures near Chesapeake Bay—typical of negative phases of the North Atlantic Oscillation and conditions thought to foster hurricane landfalls along the East Coast.

Massachusetts↗

Tsunamigenic splay faults imply a long-term asperity in southern Prince William Sound, Alaska

Coseismic slip partitioning and uplift over multiple earthquake cycles is critical to understanding upper‐plate fault development. Bathymetric and seismic reflection data from the 1964 M w 9.2 Great Alaska earthquake rupture area reveal sea floor scarps along the tsunamigenic Patton Bay/Cape Cleare/Middleton Island fault system. The faults splay from a megathrust where duplexing and underplating produced rapid exhumation. Trenchward of the duplex region, the faults produce a complex deformation pattern from oblique, south‐directed shortening at the Yakutat‐Pacific plate boundary. Spatial and temporal fault patterns suggest that Holocene megathrust earthquakes had similar relative motions and thus similar tsunami sources as in 1964. Tsunamis during future earthquakes will likely produce similar run‐up patterns and travel times. Splay fault surface expressions thus relate to plate boundary conditions, indicating millennial‐scale persistence of this asperity. We suggest structure of the subducted slab directly influences splay fault and tsunami generation landward of the frontal subduction zone prism.

Alaska↗

Thermal, deformation, and degassing remote sensing time-series (A.D. 2000-2017) at the 47 most active volcanoes in Latin America: Implications for volcanic systems

Volcanoes are hazardous to local and global populations, but only a fraction are continuously monitored by ground-based sensors. For example, in Latin America, more than 60% of Holocene volcanoes are unmonitored, meaning long-term multi-parameter datasets of volcanic activity are rare and sparse. We use satellite observations of degassing, thermal anomalies, and surface deformation spanning 17 years at 47 of the most active volcanoes in Latin America, and compare these datasets to ground-based observations archived by the Global Volcanism Program (GVP). This first comparison of multi-satellite time-series on a regional scale provides information regarding volcanic behavior during, non-, pre-, syn- and post-eruptive periods. For example, at Copahue volcano, deviations from background activity in all three types of satellite measurements were manifested months to years in advance of renewed eruptive activity in 2012. By quantifying the amount of degassing, thermal output, and deformation measured at each of these volcanoes, we test the classification of these volcanoes as open or closed volcanic systems. We find that ~28% of the volcanoes do not fall into either classification and the rest show elements of both, demonstrating a dynamic range of behavior that can change over time. Finally, we recommend how volcano monitoring could be improved through better coordination of available satellite-based capabilities and new instruments.

Journal of Geophysical Research↗

Submarine landslide kinematics derived from high-resolution imaging in Port Valdez, Alaska

Submarine landslides caused by strong ground shaking during the M9.2 1964 Great Alaska earthquake generated a tsunami that destroyed much of the old town of Valdez, Alaska, and was responsible for 32 deaths at that location. We explore structural details of the 1964 landslide deposit, as well as landslide deposits from earlier events, in order to characterize kinematics of the landslide process. We present a new high‐resolution seismic reflection data set that images the 1964 landslide deposit and six pre‐1964 deposits with great detail. These deposits are represented by thick packages (~7–23 m) of debris within >500 m of fjord sedimentation above basement. Internal slide structures are associated with distinctive landslide failure mechanisms, including detailed erosional and depositional features and structures resolved within both landslide blocks and distal debris flow layers. Based on comparisons of deposit volume from subbottom structure and differenced bathymetry, we refine prior interpretations of the source of failed material. New data show evidence for basal erosion and reworking of fjord‐floor sedimentation. Additionally, material comprising the 1964 landslide appears to have been translated and deformed by lateral thrusting, rather than having been sourced entirely from upslope evacuation zones. Taking into account these complexities in depositional patterns, we show variations in slide size through Holocene time and relate the history of landslides to the paleoseismic record. Collectively, these new observations demonstrate that Port Valdez has a repeated history of large submarine landslides, which are likely associated with large megathrust earthquakes.

Alaska↗

Interglacial paleoclimate in the Arctic

Marine Isotope Stage 11 from ~424 to 374 ka experienced peak interglacial warmth and highest global sea level ~410–400 ka. MIS 11 has received extensive study on the causes of its long duration and warmer than Holocene climate, which is anomalous in the last half million years. However, a major geographic gap in MIS 11 proxy records exists in the Arctic Ocean where fragmentary evidence exists for a seasonally sea ice-free summers and high sea-surface temperatures (SST; ~8–10 °C near the Mendeleev Ridge). We investigated MIS 11 in the western and central Arctic Ocean using 12 piston cores and several shorter cores using proxies for surface productivity (microfossil density), bottom water temperature (magnesium/calcium ratios), the proportion of Arctic Ocean Deep Water versus Arctic Intermediate Water (key ostracode species), sea ice (epipelagic sea ice dwelling ostracode abundance), and SST (planktic foraminifers). We produced a new benthic foraminiferal δ 18 O curve, which signifies changes in global ice volume, Arctic Ocean bottom temperature, and perhaps local oceanographic changes. Results indicate that peak warmth occurred in the Amerasian Basin during the middle of MIS 11 roughly from 410 to 400 ka. SST were as high as 8–10 °C for peak interglacial warmth, and sea ice was absent in summers. Evidence also exists for abrupt suborbital events punctuating the MIS 12-MIS 11-MIS 10 interval. These fluctuations in productivity, bottom water temperature, and deep and intermediate water masses (Arctic Ocean Deep Water and Arctic Intermediate Water) may represent Heinrich-like events possibly involving extensive ice shelves extending off Laurentide and Fennoscandian Ice Sheets bordering the Arctic.

Paleoceanography and Paleoclimatology↗

Radiocarbon dating of silica sinter and postglacial hydrothermal activity in the El Tatio geyser field

The El Tatio geothermal field in the Chilean Altiplano contains hydrothermal silica sinter deposits overlaying glacial and volcanic units, providing an opportunity to constrain the timing of deglaciation and volcanic activity in an area with sparse absolute chronologies. We obtained 51 new radiocarbon ages and δ 13 C values on the organic material trapped in these sinter deposits. Based on the δ 13 C values, we exclude 29 samples for possible contamination with bacterial mats that incorporate old carbon. We infer that hydrothermal activity initiated ~27 ka ago and has been nearly continuous ever since. The ages of the oldest sinter deposits coincide with ages of moraines that stabilized after the most recent deglaciation. Whereas late Pleistocene sinters are broadly distributed in the field, Holocene deposits are found around active hydrothermal features. Although recent volcanism is absent in the vicinity of El Tatio, persistent hydrothermal discharge implies a long‐lived magmatic heat source.

El Tatio Geothermal Field↗

Geophysical constraints on the crustal architecture of the transtensional Warm Springs Valley fault zone, northern Walker Lane, western Nevada, USA

The Walker Lane is a zone of distributed transtension where normal faults are overprinted by strike-slip motion. We use two newly-acquired high-resolution seismic reflection profiles and a reprocessed Consortium for Continental Reflection Profiling (COCORP) deep crustal reflection profile to assess the subsurface geometry of the Holocene-active, transtensional Warm Springs Valley fault zone (WSVFZ) near Reno, Nevada, USA. Our multi-scale observations extend to 12 km depth and suggest that the WSVFZ is more complex in the subsurface than implied by late Pleistocene surface fault traces. Two ~4-km-long high-resolution profiles image to a depth of ~2 km and reveal moderately dipping reflections and truncations, some of which project to mapped scarps formed in late Pleistocene surfaces. The shallow lines are co-located with COCORP profile NV 08 along ~40° N latitude. Re-analysis of the COCORP data reveals previously unidentified coherent reflections to a depth of ~12 km and a previously mapped ~30 west-dipping fault at 8-12 km. From these seismic profiles, the WSVFZ is not a simple, sub-vertical fault zone extending through the entire seismogenic crust. Instead, the reflections are consistent with a zone of steep- and moderately-dipping faults that simplify and steepen with depth before intersecting a mid-crustal, low angle (~25-30°) fault. The complex fault geometry of the WSVFZ implies that crustal shear is accommodated by a mix of dipping and subvertical faults in the transtensional northern Walker Lane. If so, transtensional fault zones may present challenges to paleoseismic and geodetic studies and require careful treatment when included in seismic hazard analyses.

Nevada↗

Predicted vulnerability of carbon in permafrost peatlands With future climate change and permafrost thaw in western Canada

Climate warming in high-latitude regions is thawing carbon-rich permafrost soils, which can release carbon to the atmosphere and enhance climate warming. Using a coupled model of long-term peatland dynamics (Holocene Peat Model, HPM-Arctic), we quantify the potential loss of carbon with future climate warming for six sites with differing climates and permafrost histories in Northwestern Canada. We compared the net carbon balance at 2100 CE resulting from new productivity and the decomposition of active layer and newly thawed permafrost peats under RCP8.5 as a high-end constraint. Modeled net carbon losses ranged from −3.0 kg C m −2 (net loss) to +0.1 kg C m −2 (net gain) between 2015 and 2100. Losses of newly thawed permafrost peat comprised 0.2%–25% (median: 1.6%) of “old” C loss, which were related to the residence time of peat in the active layer before being incorporated into the permafrost, peat temperature, and presence of permafrost. The largest C loss was from the permafrost-free site, not from permafrost sites. C losses were greatest from depths of 0.2–1.0 m. New C added to the profile through net primary productivity between 2015 and 2100 offset ∼40% to >100% of old C losses across the sites. Differences between modeled active layer deepening and flooding following permafrost thaw resulted in very small differences in net C loss by 2100, illustrating the important role of present-day conditions and permafrost aggradation history in controlling net C loss.

western Canada↗

A comparison of the CMIP6 midHolocene and lig127k simulations in CESM2

Results are presented and compared for the Community Earth System Model version 2 (CESM2) simulations of the middle Holocene (MH, 6 ka) and Last Interglacial (LIG, 127 ka). These simulations are designated as Tier 1 experiments ( midHolocene and lig127k ) for the Coupled Model Intercomparison Project phase 6 (CMIP6) and the Paleoclimate Modeling Intercomparison Project phase 4 (PMIP4). They use the low-top, standard 1° version of CESM2 contributing to CMIP6 DECK, historical, and future projection simulations, and to other modeling intercomparison projects. The midHolocene and lig127k provide the opportunity to examine the responses in CESM2 to the orbitally induced changes in the seasonal and latitudinal distribution of insolation. The insolation anomalies result in summer warming over the Northern Hemisphere continents, reduced Arctic summer minimum sea ice, and increased areal extent of the North African monsoon. The Arctic remains warm throughout the year. These changes are greater in the lig127k than midHolocene simulation. Other notable changes are reduction of the Niño3.4 variability and Drake Passage transport and a small increase in the Atlantic Meridional Overturning Circulation from the piControl to midHolocene to lig127k simulation. Comparisons to paleo-data and to simulations from previous model versions are discussed. Possible reasons for mismatches with the paleo-observations are proposed, including missing processes in CESM2, simplifications in the CMIP6 protocols for these experiments, and dating and calibration uncertainties in the data reconstructions.

Paleoceanography and Paleoclimatology↗

Devils Hole calcite was precipitated at ±1°C stable aquifer temperatures during the last half million years

Subaqueous carbonates from the Devils Hole caves (southwestern USA) provide a continuous Holocene to Pleistocene North American paleoclimate record. The accuracy of this record relies on two assumptions: That carbonates precipitated close to isotope equilibrium and that groundwater temperature did not change significantly in the last 570 thousand years. Here, we investigate these assumptions using dual clumped isotope thermometry. This method relies on simultaneous analyses of carbonate ∆ 47 and ∆ 48 values and provides information on the existence and extent of kinetic isotope fractionation. Our results confirm the hypothesis that calcite precipitation occurred close to oxygen and clumped isotope equilibrium during the last half million years in Devils Hole. In addition, we provide evidence that aquifer temperatures varied by less than ±1°C during this interval. Thus, the Devils Hole calcite δ 18 O time series exclusively represents changes in groundwater δ 18 O values.

Nevada↗