Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗
DOE OSTI2021

Intelligent Monitoring Systems and Advanced Well Integrity and Mitigation

Long-term seismic monitoring of carbon capture and storage projects is needed to verify that the injected gas is safely stored in the subsurface until permanence can be assured. Conventional surface seismic monitoring techniques are usually expensive, require highly invasive surface operations, and need significant time investments on the part of personnel for both the field effort and processing the acquired data. For these reasons, permanent reservoir monitoring technologies are preferred, as they can offer a cost-effective solution for long-term monitoring. As part of the monitoring program of the Archer Daniels Midland’s large-scale injection of CO 2 in Decatur, Illinois, USA, a continuous seismic monitoring array was installed using a combination of surface orbital vibrator (SOV) sources and fiber-optic cables for distributed acoustic sensing (DAS) acquisition with the objective to build a continuous monitoring array. The aim of the presented project was to build a monitoring array and platform that integrates real-time seismic data with conventional data streams and provides continuous data analysis using dynamic computational models to deliver a comprehensive real-time assessment of subsurface conditions. It is in this context that the Intelligent Monitoring Systems and Advanced Well Integrity and Mitigation project was proposed with the objective to develop an integrated architecture that utilizes a permanent seismic monitoring network, combines the real-time geophysical and process data with reservoir flow and geomechanical models to create a comprehensive monitoring, visualization, and control system that delivers critical information for process surveillance and optimization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Inline small-angle X-ray scattering-coupled chromatography under extreme hydrostatic pressure

As continuing discoveries highlight the surprising abundance and resilience of deep ocean and subsurface microbial life, the effects of extreme hydrostatic pressure on biological structure and function have attracted renewed interest. Biological small-angle X-ray scattering (BioSAXS) is a widely used method of obtaining structural information from biomolecules in solution under a wide range of solution conditions. Due to its ability to reduce radiation damage, remove aggregates, and separate monodisperse components from complex mixtures, size-exclusion chromatography-coupled SAXS (SEC-SAXS) is now the dominant form of BioSAXS at many synchrotron beamlines. While BioSAXS can currently be performed with some difficulty under pressure with non-flowing samples, it has not been clear how, or even if, continuously flowing SEC-SAXS, with its fragile media-packed columns, might work in an extreme high-pressure environment. Here we show, for the first time, that reproducible chromatographic separations coupled directly to high-pressure BioSAXS can be achieved at pressures up to at least 100 MPa and that pressure-induced changes in folding and oligomeric state and other properties can be observed. Further, the apparatus described here functions at a range of temperatures (0°C–50°C), expanding opportunities for understanding biomolecular rules of life in deep ocean and subsurface environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evaluation of several relevant fractionation processes as possible explanation for radioxenon isotopic activity ratios in samples taken near underground nuclear explosions in shafts and tunnels

Gas samples taken from two historic underground nuclear tests done in 1989 at the Nevada National Security Site (NNSS), formerly the Nevada Test Site (NTS), were examined to determine how xenon isotopes fractionate because of early-time cavity processes, transport through the rock, or dispersal through tunnels. Xenon isotopes are currently being used to distinguish civilian sources of xenon in the atmosphere from sources associated with underground nuclear explosions (UNEs). The two nuclear tests included (1) BARNWELL, a test conducted in a vertical shaft approximately 600 m below ground surface at Pahute Mesa, and (2) DISKO ELM, a horizontal line-of-sight test done in P-tunnel approximately 261 m below the surface of Aqueduct Mesa. Numerical flow and transport models developed for the two sites had mixed success when attempting to match the observed xenon isotope ratios. At the BARNWELL site, the simulated xenon isotope ratios were consistent with measurements from the chimney and ground surface and appeared to have been affected primarily by fractionation during subsurface transport. At the DISKO ELM site, samples taken from two elevations in the chimney failed to show the degree of fractionation predicted by the models during transport and did not show evidence for significant fractionation due to early-time condensation of refractory xenon-precursor radionuclides into the melt glass. Gas samples taken from the adjacent tunnels in the days following the test showed mixed evidence for early-time separation of xenon isotopes from their iodine precursors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The Role of Atmospheric Rivers on Groundwater: Lessons Learned From an Extreme Wet Year

Abstract In the coastal regions of the western United States, atmospheric rivers (ARs) are associated with the largest precipitation generating storms and contribute up to half of annual precipitation, but the impact of ARs on the integrated hydrologic cycle, specifically on groundwater storage and hydrodynamics, is largely unknown. To better explore the hydrologic behavior of AR versus non‐AR event precipitation, we present a novel combination of two water tracking methods (one in the atmosphere and one in the subsurface) to explicitly track the full lifecycle of water parcels generated by ARs. Simulations of northern California's Cosumnes River watershed during the record wet 2017 water year are performed via the coupling of a high‐resolution regional climate model and a land surface‐groundwater model accounting for lateral groundwater flow. Despite ARs contributing more precipitation than non‐AR storms, we find less AR water is preferentially stored in aquifers by year end. Fractionally, ARs result in 300% less snow derived groundwater‐recharged compared to non‐AR precipitation. Rain‐on‐snow (RoS) plays an important role in AR‐driven discharge, where over 50% of total discharge from ARs snow is from RoS events. Finally, despite record‐breaking annual precipitation, simulated groundwater depletion occurs by year end due to estimates of groundwater pumping activities. The results from these simulations serve as a partial analogue of future hydrologic conditions where ARs are expected to intensify and provide a greater fraction of annual precipitation due to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2023

Combining synchrotron X-ray diffraction, mechanistic modeling and machine learning for in situ subsurface temperature quantification during laser melting

Laser melting, such as that encountered during additive manufacturing, produces extreme gradients of temperature in both space and time, which in turn influence microstructural development in the material. Qualification and model validation of the process itself and the resulting material necessitate the ability to characterize these temperature fields. However, well established means to directly probe the material temperature below the surface of an alloy while it is being processed are limited. To address this gap in characterization capabilities, a novel means is presented to extract subsurface temperature-distribution metrics, with uncertainty, from in situ synchrotron X-ray diffraction measurements to provide quantitative temperature evolution data during laser melting. Temperature-distribution metrics are determined using Gaussian process regression supervised machine-learning surrogate models trained with a combination of mechanistic modeling (heat transfer and fluid flow) and X-ray diffraction simulation. The trained surrogate model uncertainties are found to range from 5 to 15% depending on the metric and current temperature. The surrogate models are then applied to experimental data to extract temperature metrics from an Inconel 625 nickel superalloy wall specimen during laser melting. The maximum temperatures of the solid phase in the diffraction volume through melting and cooling are found to reach the solidus temperature as expected, with the mean and minimum temperatures found to be several hundred degrees less. The extracted temperature metrics near melting are determined to be more accurate because of the lower relative levels of mechanical elastic strains. However, uncertainties for temperature metrics during cooling are increased due to the effects of thermomechanical stress.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Using AI to build a hydrobiogeochemical soil model

Soil water content is a function of inputs from precipitation and outputs via evaporation, transpiration, lateral flow, and vertical percolation, and is sensitive to biogeochemical processes. As such, soils serve as an ideal integrator of atmospheric, hydrological, and biogeochemical processes affecting the water cycle. In addition, soil water retention capacity, infiltration rates, and hydraulic conductivity can buffer or exacerbate the effects of extreme precipitation events (e.g., flooding, runoff, subsurface transport, erosion, greenhouse gas emissions) and mitigate the impact of droughts and heat waves on land systems (e.g., fire, crop failure). However, integrating water cycle measurements spanning different land atmosphere compartments across scales is a fundamental barrier for numerical model predictability. A significant challenge is that each domain (soil, hydrology, biology, and atmosphere) typically collects different sets of data at different temporal and spatial frequencies/scales, and even different dimensionalities (2D vs 3D). To implement soil as an integrator of the water cycle in land models, we suggest that novel machine learning (ML) tools can be developed to effectively simulate complex landscapes across various domains and scales, extended to regions with sparse or no data. The ultimate goals are to improve predictive understanding of land-atmosphere interactions and to extend the predictability of current Earth System Models (ESMs) through better integration of hydrological and biogeochemical data. We envision a framework in which: (1) ML-aided data reconstructions enable the merger of data sources into a unified geospatial product; (2) automated detection techniques are used to improve the knowledge of complex soil processes and interactions; and (3) this knowledge is leveraged and incorporated into models through AI-based emulators to distinctly connect the land and atmospheric compartments of the water cycle in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

(Project 18-15502) Reducing Uncertainty in Radionuclide Transport Prediction Using Multiple Environmental Tracers (Final Report)

In order to successfully site and design a nuclear waste disposal facility, DOE scientists are required to show safe containment of the radioactive waste for up to 1 million years. A significant hurdle to accurately and convincingly demonstrating disposal safety is predicting the fate of radioactive elements once they enter the groundwater system surrounding the repository. These predictions are often made with computer models, which contain accurate physics of groundwater movement and chemical reactions that occur during groundwater flow. As computer power increases, the physics and chemistry of these computer simulators can become more and more realistic and the physically based error decreases. However, the ability of these computer models to provide accurate predictions in a specific place, over long time periods, requires the scientists and engineers to know the subsurface properties of the Earth that control groundwater movement. In particular, groundwater scientists and engineers need to know the groundwater fluid velocity, which can change over time and strongly vary with location within the groundwater system. Because the Earth’s subsurface cannot be directly seen, and can only be sampled at drilling locations, the properties of the groundwater system are never known completely, and computer models of groundwater transport will always have some amount of uncertainty. This project’s principal goal was to use chemicals and isotope “tracers”, which have been introduced to the groundwater system by natural processes over long time periods, to help inform computer models of the groundwater velocity and subsurface properties. The goal was to calculate how much better predictions of groundwater transport were when these “tracers” were used to inform the computer models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗
DOE OSTI2025

Addressing Critical Challenges of Coal Mines for Underground Pumped Storage Hydropower Using Numerical Calculations

Using existing coal mines to build underground pumped storage hydropower (UPSH) power plants could help reduce the costs associated with starting a new PSH operation. This innovative approach can leverage the existing infrastructure of coal mine facilities available in many regions of the United States. However, technical challenges associated with water quality, mineral chemistry and integrity of subsurface topologies have thus far prevented any of these facilities from going online. Rye Development is investigating the feasibility of UPSH at existing coal mine facilities. In this project, Rye Development received the technical assistance from Oak Ridge National Laboratory to investigate the impacts of minerals, metal particles, and chemical substances in water flow on hydraulic equipment and assess the overall structural integrity of hydraulic equipment in former coal mines.

13 HYDRO ENERGY↗
DOE OSTI2021

Investigating Fracture Network Deformation Using Noble Gas Release

We investigate deformation mechanics of fracture networks in unsaturated fractured rocks from subsurface conventional detonation using dynamic noble gas measurements and changes in air permeability. We dynamically measured the noble gas isotopic composition and helium exhalation of downhole gas before and after a large subsurface conventional detonation. These noble gas measurements were combined with measurements of the subsurface permeability field from 64 discrete sampling intervals before and after the detonation and subsurface mapping of fractures in borehole walls before well completion. We saw no observable increase in radiogenic noble gas release from either an isotopic composition or a helium exhalation point of view. Large increases in permeability were observed in 13 of 64 discrete sampling intervals. Of the sampling intervals which saw large increases in flow, only two locations did not have preexisting fractures mapped at the site. Given the lack of noble gas release and a clear increase in permeability, we infer that most of the strain accommodation of the fractured media occurred along previously existing fractures, rather than the creation of new fractures, even for a high strain rate event. These results have significant implications for how we conceptualize the deformation of rocks with fracture networks above the percolation threshold, with application to a variety of geologic and geological engineering problems.

Gardner, W. Payton↗
DOE OSTI2020

Understanding the distinct impacts of MCS and non-MCS rainfall on the surface water balance in the central US using a numerical water-tagging technique

Warm-season rainfall associated with mesoscale convective systems (MCSs) in the central US is characterized by higher intensity and nocturnal timing compared to rainfall from non-MCS systems, suggesting their potentially different footprints on the land surface. To differentiate the impacts of MCS and non-MCS rainfall on the surface water balance, a water tracer tool embedded in the Noah land surface model with multi-parameterization options (WT-Noah-MP) is used to numerically “tag” water from MCS and non-MCS rainfall separately during April to August (1997-2018) and track their transit in the terrestrial system. From the water-tagging results, over 50% of warm-season rainfall leaves the surface-subsurface system through evapotranspiration by the end of August, but non-MCS rainfall contributes a larger fraction. However, MCS rainfall plays a more important role in generating surface runoff. These differences are mostly attributed to the rainfall intensity differences. The higher intensity MCS rainfall tends to produce more surface runoff through infiltration excess flow and drives a deeper penetration of the rainwater into the soil. Over 70% of the top 10 percentile runoff is contributed by MCS rainfall, demonstrating its important contribution to local flooding. In contrast, lighter intensity non-MCS rainfall resides mostly in the top layer and contributes more to evapotranspiration through soil evaporation. Diurnal timing of rainfall has negligible effects on the flux partitioning for both MCS and non-MCS rainfall. Differences in soil moisture profiles for MCS and non-MCS rainfall and the resultant evapotranspiration suggest differences in their roles in soil moisture-precipitation feedbacks and ecohydrology.

Hu, Huancui↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2022

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗
DOE OSTI2023

Model-parallel Fourier neural operators as learned surrogates for large-scale parametric PDEs

Fourier neural operators (FNOs) are a recently introduced neural network architecture for learning solution operators of partial differential equations (PDEs), which have been shown to perform significantly better than comparable deep learning approaches. Once trained, FNOs can achieve speed-ups of multiple orders of magnitude over conventional numerical PDE solvers. However, due to the high dimensionality of their input data and network weights, FNOs have so far only been applied to two-dimensional or small three-dimensional problems. To remove this limited problem-size barrier, we propose a model-parallel version of FNOs based on domain-decomposition of both the input data and network weights. Here, we demonstrate that our model-parallel FNO is able to predict time-varying PDE solutions of over 2.6 billion variables on Perlmutter using up to 512 A100 GPUs and show an example of training a distributed FNO on the Azure cloud for simulating multiphase CO 2 dynamics in the Earth’s subsurface.

58 GEOSCIENCES↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The role of Northeast Pacific meltwater events in deglacial climate change

Columbia River megafloods occurred repeatedly during the last deglaciation, but the impacts of this fresh water on Pacific hydrography are largely unknown. To reconstruct changes in ocean circulation during this period, we used a numerical model to simulate the flow trajectory of Columbia River megafloods and compiled records of sea surface temperature, paleo-salinity, and deep-water radiocarbon from marine sediment cores in the Northeast Pacific. The North Pacific sea surface cooled and freshened during the early deglacial (19.0-16.5 ka) and Younger Dryas (12.9-11.7 ka) intervals, coincident with the appearance of subsurface water masses depleted in radiocarbon relative to the sea surface. We infer that Pacific meltwater fluxes contributed to net Northern Hemisphere cooling prior to North Atlantic Heinrich Events, and again during the Younger Dryas stadial. Abrupt warming in the Northeast Pacific similarly contributed to hemispheric warming during the Bølling and Holocene transitions. These findings underscore the importance of changes in North Pacific freshwater fluxes and circulation in deglacial climate events.

54 ENVIRONMENTAL SCIENCES↗