SEARCH · USGS Science

Results for “Earth Observer”

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 541 records · Page 30Linked to original sources

USGS1959

Composition of monazites from pegmatites in eastern Minas Gerais, Brazil

Two zoned pegmatites in south-eastern Minas Gerais were sampled in detail for their content of monazite and xenotime and the monazite was analysed for certain of the rare-earth elements and thorium. The ratio of xenotime to monazite increases in both pegmatites from the wall toward the quartz core. The content of the less basic rare-earth elements and of thorium in monazite rises in the same direction. These variation trends suggest that during the crystallization of these pegmatites there was a fractionation of the elements leading to a more or less steady enrichment of the less basic rare-earth elements and of thorium in the residual fluids. One mode of explaining these observed effects postulates that the rare-earth elements and thorium were present in pegmatitic fluids as co-ordination complexes rather than as simple cations.

Minas Gerais↗

USGS2022

Improved resolution across the Global Seismographic Network: A new era in low-frequency seismology

The Global Seismographic Network (GSN)—a global network of ≈150 very broadband stations—is used by researchers to study the free oscillations of the Earth (≈0.3–10 mHz) following large earthquakes. Normal‐mode observations can provide information about the radial density and anisotropic velocity structure of the Earth (including near the core–mantle boundary), but only when signal‐to‐noise ratios at very low frequencies are sufficiently high. Most normal‐mode observations in the past three decades have been made using Streckeisen STS‐1 vault seismometers. However, these sensors are no longer being manufactured or serviced. Candidate replacement sensors, the Streckeisen STS‐6 and the Nanometrics T‐360GSN, have been recently installed in boreholes, postholes, and vaults at several GSN stations and GSN testbeds. In this study, we examine normal‐mode spectra following three M w "> M w 8 earthquakes in 2021 and from one M w "> M w 8.2 earthquake in 2014 to evaluate the change in GSN low‐frequency performance on the vertical component. From this analysis, we conclude that the number of GSN stations capable of resolving normal modes following M w "> M w 8 earthquakes has nearly doubled since 2014. The improved observational capabilities will help better understand the radial velocity and density estimates of the Earth.

The Seismic Record↗

USGS2008

Towards monitoring land-cover and land-use changes at a global scale: the global land survey 2005

Land cover is a critical component of the Earth system, infl uencing land-atmosphere interactions, greenhouse gas fl uxes, ecosystem health, and availability of food, fi ber, and energy for human populations. The recent Integrated Global Observations of Land (IGOL) report calls for the generation of maps documenting global land cover at resolutions between 10m and 30m at least every fi ve years (Townshend et al., in press). Moreover, despite 35 years of Landsat observations, there has not been a unifi ed global analysis of land-cover trends nor has there been a global assessment of land-cover change at Landsat-like resolution. Since the 1990s, the National Aeronautics and Space Administration (NASA) and the U.S. Geological Survey (USGS) have supported development of data sets based on global Landsat observations (Tucker et al., 2004). These land survey data sets, usually referred to as GeoCover ™, provide global, orthorectifi ed, typically cloud-free Landsat imagery centered on the years 1975, 1990, and 2000, with a preference for leaf-on conditions. Collectively, these data sets provided a consistent set of observations to assess land-cover changes at a decadal scale. These data are freely available via the Internet from the USGS Center for Earth Resources Observation and Science (EROS) (see http://earthexplorer.usgs.gov or http://glovis.usgs.gov). This has resulted in unprecedented downloads of data, which are widely used in scientifi c studies of land-cover change (e.g., Boone et al., 2007; Harris et al., 2005; Hilbert, 2006; Huang et al. 2007; Jantz et al., 2005, Kim et al., 2007; Leimgruber, 2005; Masek et al., 2006). NASA and USGS are continuing to support land-cover change research through the development of GLS2005 - an additional global Landsat assessment circa 20051 . Going beyond the earlier initiatives, this data set will establish a baseline for monitoring changes on a 5-year interval and will pave the way toward continuous global land-cover monitoring at Landsat-like resolution in the next decade.

Photogrammetric Engineering and Remote Sensing↗

USGS2016

Helium as a tracer for fluids released from Juan de Fuca lithosphere beneath the Cascadia forearc

The ratio between helium isotopes (3He/4He) provides an excellent geochemical tracer for investigating the sources of fluids sampled at the Earth's surface. 3He/4He values observed in 25 mineral springs and wells above the Cascadia forearc document a significant component of mantle-derived helium above Juan de Fuca lithosphere, as well as variability in 3He enrichment across the forearc. Sample sites arcward of the forearc mantle corner (FMC) generally yield significantly higher ratios (1.2-4.0 RA) than those seaward of the corner (0.03-0.7 RA). The highest ratios in the Cascadia forearc coincide with slab depths (40-45 km) where metamorphic dehydration of young oceanic lithosphere is expected to release significant fluid and where tectonic tremor occurs, whereas little fluid is expected to be released from the slab depths (25-30 km) beneath sites seaward of the corner.Tremor (considered a marker for high fluid pressure) and high RA values in the forearc are spatially correlated. The Cascadia tremor band is centered on its FMC, and we tentatively postulate that hydrated forearc mantle beneath Cascadia deflects a significant portion of slab-derived fluids updip along the subduction interface, to vent in the vicinity of its corner. Furthermore, high RA values within the tremor band just arcward of the FMC, suggest that the innermost mantle wedge is relatively permeable.Conceptual models require: (1) a deep fluid source as a medium to transport primordial 3He; (2) conduits through the lithosphere which serve to speed fluid ascent to the surface before significant dilution from radiogenic 4He can occur; and (3) near lithostatic fluid pressure to keep conduits open. Our spatial correlation between high RA values and tectonic tremor provides independent evidence that tremor is associated with deep fluids, and it further suggests that high pore pressures associated with tremor may serve to keep fractures open for 3He migration through ductile upper mantle and lower crust.

Geochemistry International↗

USGS2008

Data-resolution matrix and model-resolution matrix for Rayleigh-wave inversion using a damped least-squares method

Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (>2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. We employed a data-resolution matrix to select data that would be well predicted and we find that there are advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher-mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher-mode data are normally more accurately predicted than fundamental-mode data because of restrictions on the data kernel for the inversion system. We used synthetic and real-world examples to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher-mode data in inversion can provide better results. We also calculated model-resolution matrices in these examples to show the potential of increasing model resolution with selected surface-wave data. ?? Birkhaueser 2008.

Pure and Applied Geophysics↗

USGS2018

Powerful turbidity currents driven by dense basal layers

Seafloor sediment flows (turbidity currents) are among the volumetrically most important yet least documented sediment transport processes on Earth. A scarcity of direct observations means that basic characteristics, such as whether flows are entirely dilute or driven by a dense basal layer, remain equivocal. Here we present the most detailed direct observations yet from oceanic turbidity currents. These powerful events in Monterey Canyon have frontal speeds of up to 7.2 m s −1 , and carry heavy (800 kg) objects at speeds of ≥4 m s −1 . We infer they consist of fast and dense near-bed layers, caused by remobilization of the seafloor, overlain by dilute clouds that outrun the dense layer. Seabed remobilization probably results from disturbance and liquefaction of loose-packed canyon-floor sand. Surprisingly, not all flows correlate with major perturbations such as storms, floods or earthquakes. We therefore provide a new view of sediment transport through submarine canyons into the deep-sea.

California↗

USGS2010

Marine-target craters on Mars? An assessment study

Observations of impact craters on Earth show that a water column at the target strongly influences lithology and morphology of the resultant crater. The degree of influence varies with the target water depth and impactor diameter. Morphological features detectable in satellite imagery include a concentric shape with an inner crater inset within a shallower outer crater, which is cut by gullies excavated by the resurge of water. In this study, we show that if oceans, large seas, and lakes existed on Mars for periods of time, marine-target craters must have formed. We make an assessment of the minimum and maximum amounts of such craters based on published data on water depths, extent, and duration of putative oceans within “contacts 1 and 2,” cratering rate during the different oceanic phases, and computer modeling of minimum impactor diameters required to form long-lasting craters in the seafloor of the oceans. We also discuss the influence of erosion and sedimentation on the preservation and exposure of the craters. For an ocean within the smaller “contact 2” with a duration of 100,000 yr and the low present crater formation rate, only ˜1–2 detectable marine-target craters would have formed. In a maximum estimate with a duration of 0.8 Gyr, as many as 1400 craters may have formed. An ocean within the larger “contact 1-Meridiani,” with a duration of 100,000 yr, would not have received any seafloor craters despite the higher crater formation rate estimated before 3.5 Gyr. On the other hand, with a maximum duration of 0.8 Gyr, about 160 seafloor craters may have formed. However, terrestrial examples show that most marine-target craters may be covered by thick sediments. Ground penetrating radar surveys planned for the ESA Mars Express and NASA 2005 missions may reveal buried craters, though it is uncertain if the resolution will allow the detection of diagnostic features of marine-target craters. The implications regarding the discovery of marine-target craters on Mars is not without significance, as such discoveries would help address the ongoing debate of whether large water bodies occupied the northern plains of Mars and would help constrain future paleoclimatic reconstructions.

Meteoritics and Planetary Science↗

USGS2016

Leaf development and demography explain photosynthetic seasonality in Amazon evergreen forests

In evergreen tropical forests, the extent, magnitude, and controls on photosynthetic seasonality are poorly resolved and inadequately represented in Earth system models. Combining camera observations with ecosystem carbon dioxide fluxes at forests across rainfall gradients in Amazônia, we show that aggregate canopy phenology, not seasonality of climate drivers, is the primary cause of photosynthetic seasonality in these forests. Specifically, synchronization of new leaf growth with dry season litterfall shifts canopy composition toward younger, more light-use efficient leaves, explaining large seasonal increases (~27%) in ecosystem photosynthesis. Coordinated leaf development and demography thus reconcile seemingly disparate observations at different scales and indicate that accounting for leaf-level phenology is critical for accurately simulating ecosystem-scale responses to climate change.

Science↗

USGS2009

Deep drilling in the Chesapeake Bay impact structure - An overview

The late Eocene Chesapeake Bay impact structure lies buried at moderate depths below Chesapeake Bay and surrounding landmasses in southeastern Virginia, USA. Numerous characteristics made this impact structure an inviting target for scientific drilling, including the location of the impact on the Eocene continental shelf, its threelayer target structure, its large size (??85 km diameter), its status as the source of the North American tektite strewn field, its temporal association with other late Eocene terrestrial impacts, its documented effects on the regional groundwater system, and its previously unstudied effects on the deep microbial biosphere. The Chesapeake Bay Impact Structure Deep Drilling Project was designed to drill a deep, continuously cored test hole into the central part of the structure. A project workshop, funding proposals, and the acceptance of those proposals occurred during 2003-2005. Initial drilling funds were provided by the International Continental Scientific Drilling Program (ICDP) and the U.S. Geological Survey (USGS). Supplementary funds were provided by the National Aeronautics and Space Administration (NASA) Science Mission Directorate, ICDP, and USGS. Field operations were conducted at Eyreville Farm, Northampton County, Virginia, by Drilling, Observation, and Sampling of the Earth's Continental Crust (DOSECC) and the project staff during September-December 2005, resulting in two continuously cored, deep holes. The USGS and Rutgers University cored a shallow hole to 140 m in April-May 2006 to complete the recovered section from land surface to 1766 m depth. The recovered section consists of 1322 m of crater materials and 444 m of overlying postimpact Eocene to Pleistocene sediments. The crater section consists of, from base to top: basement-derived blocks of crystalline rocks (215 m); a section of suevite, impact melt rock, lithic impact breccia, and cataclasites (154 m); a thin interval of quartz sand and lithic blocks (26 m); a granite megablock (275 m); and sediment blocks and boulders, polymict, sediment-clast-dominated sedimentary breccias, and a thin upper section of stratified sediments (652 m). The cored postimpact sediments provide insight into the effects of a large continental-margin impact on subsequent coastal-plain sedimentation. This volume contains the first results of multidisciplinary studies of the Eyreville cores and related topics. The volume is divided into these sections: geologic column; borehole geophysical studies; regional geophysical studies; crystalline rocks, impactites, and impact models; sedimentary breccias; postimpact sediments; hydrologic and geothermal studies; and microbiologic studies. ?? 2009 The Geological Society of America.

Special Paper of the Geological Society of America↗

USGS2008

A trade-off between model resolution and variance with selected Rayleigh-wave data

Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (??? 2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. First, we employed a data-resolution matrix to select data that would be well predicted and to explain advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher mode data are normally more accurately predicted than fundamental mode data because of restrictions on the data kernel for the inversion system. Second, we obtained an optimal damping vector in a vicinity of an inverted model by the singular value decomposition of a trade-off function of model resolution and variance. In the end of the paper, we used a real-world example to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher mode data in inversion can provide better results. We also calculated model-resolution matrices of these examples to show the potential of increasing model resolution with selected surface-wave data. With the optimal damping vector, we can improve and assess an inverted model obtained by a damped least-square method.

Conference Paper↗

USGS2010

Seismic evidence for widespread western-US deep-crustal deformation caused by extension

Laboratory experiments have established that many of the materials comprising the Earth are strongly anisotropic in terms of seismic-wave speeds. Observations of azimuthal and radial anisotropy in the upper mantle are attributed to the lattice-preferred orientation of olivine caused by the shear strains associated with deformation, and provide some of the most direct evidence for deformation and flow within the Earths interior. Although observations of crustal radial anisotropy would improve our understanding of crustal deformation and flow patterns resulting from tectonic processes, large-scale observations have been limited to regions of particularly thick crust. Here we show that observations from ambient noise tomography in the western United States reveal strong deep (middle to lower)-crustal radial anisotropy that is confined mainly to the geological provinces that have undergone significant extension during the Cenozoic Era (since 65 Myr ago). The coincidence of crustal radial anisotropy with the extensional provinces of the western United States suggests that the radial anisotropy results from the lattice-preferred orientation of anisotropic crustal minerals caused by extensional deformation. These observations also provide support for the hypothesis that the deep crust within these regions has undergone widespread and relatively uniform strain in response to crustal thinning and extension. ?? 2010 Macmillan Publishers Limited. All rights reserved.

Nature↗

USGS2016

Potential evapotranspiration and continental drying

By various measures (drought area and intensity , climatic aridity index , and climatic water deficits ), some observational analyses have suggested that much of the Earth ’ s land has been drying during recent decades, but such drying seems inconsistent with observations of dryland greening and decreasing pan evaporation . ‘Offline ’ analyses of climate-model outputs from anthropogenic climate change (ACC) experiments portend continuation of putative drying through the twenty-first century , despite an expected increase in global land precipitation . A ubiquitous increase in estimates of potential evapotranspiration (PET), driven by atmospheric warming , underlies the drying trends , but may be a methodological artefact . Here we show that the PET estimator commonly used (the Penman–Monteith PET for either an open-water surface or a reference crop ) severely overpredicts the changes in non-water-stressed evapotranspiration computed in the climate models themselves in ACC experiments. This overprediction is partially due to neglect of stomatal conductance reductions commonly induced by increasing atmospheric CO 2 concentrations in climate models . Our findings imply that historical and future tendencies towards continental drying, as characterized by offline-computed runoff, as well as other PET-dependent metrics, may be considerably weaker and less extensive than previously thought.

Nature Climate Change↗

USGS2012

State of the Earth’s cryosphere at the beginning of the 21st century: Glaciers, global snow cover, floating ice, and permafrost and periglacial environments

This chapter is the tenth in a series of 11 book-length chapters, collectively referred to as “this volume,” in the series U.S. Geological Survey Professional Paper 1386, Satellite Image Atlas of Glaciers of the World. In the other 10 chapters, each of which concerns a specific glacierized region of Earth, the authors used remotely sensed images, primarily from the Landsat 1, 2, and 3 series of spacecraft, in order to analyze that glacierized region and to monitor changes in its glaciers. Landsat images, acquired primarily during the period 1972 through 1981, were used by an international team of glaciologists and other scientists to study the various glacierized regions and (or) to discuss related glaciological topics. In each glacierized region, the present distribution of glaciers within its geographic area is compared, wherever possible, with historical information about their past areal extent. The atlas provides an accurate regional inventory of the areal extent of glacier ice on our planet during the 1970s as part of an expanding international scientific effort to measure global environmental change on the Earth’s surface. However, this chapter differs from the other 10 in its discussion of observed changes in all four elements of the Earth’s cryosphere (glaciers, snow cover, floating ice, and permafrost) in the context of documented changes in all components of the Earth System. Human impact on the planet at the beginning of the 21st century is pervasive. The focus of Chapter A is on changes in the cryosphere and the importance of long-term monitoring by a variety of sensors carried on Earth-orbiting satellites or by a ground-based network of observatories in the case of permafrost. The chapter consists of five parts. The first part provides an introduction to the Earth System, including the interrelationships of the geosphere (cryosphere, hydrosphere, lithosphere, and atmosphere), the biosphere, climate processes, biogeochemical cycles, and the critically important hydrologic cycle, in which glacier ice is the second largest reservoir of water after the oceans. The second part assesses the state of glaciers in all of the glacierized regions of the planet, primarily as drawn in the other 10 chapters. It includes sections on ice cores and the climate record they contain, volumetric changes in glaciers, harnessing spaceborne sensors to measure changes in glaciers, and related topics. The third part summarizes trends in global snow cover. The fourth part summarizes long-term changes in area and thickness of floating ice, including polar sea ice and freshwater (lake and river) ice. The fifth part assesses the loss of permafrost and changes in periglacial environments at high latitudes and high altitudes.

Professional Paper↗

USGS2012

A dislocation model of strain accumulation and release at a subduction zone

Strain accumulation and release at a subduction zone are attributed to stick slip on the main thrust zone and steady aseismic slip on the remainder of the plate interface. This process can be described as a superposition of steady state subduction and a repetitive cycle of slip on the main thrust zone, consisting of steady normal slip at the plate convergence rate plus occasional thrust events that recover the accumulated normal slip. Because steady state subduction does not contribute to the deformation at the free surface, deformation observed there is completely equivalent to that produced by the slip cycle alone. The response to that slip is simply the response of a particular earth model to embedded dislocations. For a purely elastic earth model, the deformation cycle consists of a coseismic offset followed by a linear‐in‐time recovery to the initial value during the interval between earthquakes. For an elastic‐viscoelastic earth model (elastic lithosphere over a viscoelastic asthenosphere), the postearthquake recovery is not linear in time. Records of local uplift as a function of time indicate that the long‐term postseismic recovery is approximately linear, suggesting that elastic earth models are adequate to describe the deformation cycle. However, the deformation predicted for a simple elastic half‐space earth model does not reproduce the deformation observed along the subduction zones in Japan at all well if stick slip is restricted to the main thrust zone. As recognized earlier by Shimazaki, Seno, and Kato, the uplift profiles could be explained if stick slip were postulated to extend along the plate interface beyond the main thrust zone to a depth of perhaps 100 km, but independent evidence suggests that stick slip at such depths is unlikely.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

Damage and restoration of geodetic infrastructure caused by the 1994 Northridge, California, earthquake

We seek to restore the integrity of the geodetic network in the San Fernando, Simi, Santa Clarita Valleys and in the northern Los Angeles Basin by remeasurement of the network and identification of BMs which experienced non-tectonic displacements associated with the Northridge earthquake. We then use the observed displacement of BMs in the network to portray or predict the permanent vertical and horizontal deformation associated with the 1994 Northridge earthquake throughout the area, including sites where we lack geodetic measurements. To accomplish this, we find the fault geometry and earthquake slip that are most compatible with the geodetic and independent seismic observations of the earthquake. We then use that fault model to predict the deformation everywhere at the earth's surface, both at locations where geodetic observations exist and also where they are absent. We compare displacements predicted for a large number of numerical models of the earthquake faulting to the coseismic displacements, treating the earthquake fault as a cut or discontinuity embedded in a stiff elastic solid. This comparison is made after non-tectonic deformation has been removed from the measured elevation changes. The fault slip produces strain in the medium and deforms the ground surface. The model compatible with seismic observations that best fits the geodetic data within their uncertainties is selected. The acceptable model fault bisects the mainshock focus, and the earthquake size , magnitude, is compatible with the earthquake size measured seismically. Our fault model was used to identify geodetic monuments on engineered structures that were anomalously displaced by the earthquake.

California↗

USGS2026

The use of fluorite geochemistry and machine learning to identify critical mineral systems

Fluorite (CaF 2 ) is a potential pathfinder to critical mineral and rare earth element (REE) deposits but its application has been limited to a narrow range of mineralization types. I show that fluorite is a robust recorder of mineralization fertility by applying statistical and machine-learning methods to a new global fluorite geochemical database. Distinct median rare earth and trace element patterns are observed among deposit types and genetic environments. Fluorite associated with carbonatites and REE deposits are relatively enriched in Sr and have minimal Eu anomalies. These characteristics define new bivariate discrimination diagrams that correctly identify 78% of carbonatite-related fluorite and 88% of fluorite from REE deposits. Random forest classifiers were developed for a wide range of mineralization types and genetic settings. Trained solely on rare earth element patterns, these models achieve accuracies of 77–79%. Higher classification accuracies (up to 88–96%) are obtained when including elements such as Sr, highlighting the significance of trace elements for optimal fluorite classification. The recognition of diagnostic fluorite compositional fingerprints, particularly in REE-fertile systems, underscores its potential as a pathfinder and indicator for critical mineral exploration in F-bearing environments.

Mineralium Deposita↗

USGS2023

Shallow fault slip of the 2020 M5.1 Sparta, North Carolina, earthquake

The 2020 M 5.1 Sparta, North Carolina, earthquake is the largest in the eastern United States since the 2011 M 5.8 Mineral, Virginia, earthquake and produced a ∼2.5‐km‐long surface rupture, unusual for an event of this magnitude. A geological field study conducted soon after the event indicates oblique slip along a east‐southeast‐trending fault with a consistently observed thrust component. My analysis of regional seismic waveforms, Interferometric Synthetic Aperture Radar, and Global Positioning System survey data yields a compact shallow rupture extending from Earth’s surface down‐dip to the southwest over a ∼3 km fault length. The inferred kinematic rupture is primarily toward the up‐dip and eastward along‐strike directions and has predominantly thrust motion in the west, transitioning to roughly equal thrust and left‐lateral strike‐slip motion in the east. No normal faulting component, as proposed in an earlier geophysical study, is necessary to explain the data. The prevalence of only dip‐slip motions observed at Earth’s surface may demand slip partitioning between dip slip and lateral motions at depth.

North Carolina, Virginia, West Virginia↗

USGS2019

Most Earth-surface calcites precipitate out of isotopic equilibrium

Oxygen-isotope thermometry played a critical role in the rise of modern geochemistry and remains extensively used in (bio-)geoscience. Its theoretical foundations rest on the assumption that 18 O/ 16 O partitioning among water and carbonate minerals primarily reflects thermodynamic equilibrium. However, after decades of research, there is no consensus on the true equilibrium 18 O/ 16 O fractionation between calcite and water ( 18 α cc/w ). Here, we constrain the equilibrium relations linking temperature, 18 α cc/w , and clumped isotopes ( Δ 47 ) based on the composition of extremely slow-growing calcites from Devils Hole and Laghetto Basso (Corchia Cave). Equilibrium 18 α cc/w values are systematically ~1.5‰ greater than those in biogenic and synthetic calcite traditionally considered to approach oxygen-isotope equilibrium. We further demonstrate that subtle disequilibria also affect Δ 47 in biogenic calcite. These observations provide evidence that most Earth-surface calcites fail to achieve isotopic equilibrium, highlighting the need to improve our quantitative understanding of non-equilibrium isotope fractionation effects instead of relying on phenomenological calibrations.

Nature Communications↗