Search NASA⌕ Search

SEARCH · Search NASA

Results for “POTENTIAL GRADIENT”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Non-linear Plasma Wake Growth of Electron Holes

An object’s wake in a plasma with small Debye length that drifts across the magnetic field is subject to electrostatic electron instabilities. Such situations include, for example, the moon in the solar wind wake and probes in magnetized laboratory plasmas. The instability drive mechanism can equivalently be considered drift down the potential-energy gradient or drift up the density gradient. The gradients arise because the plasma wake has a region of depressed density and electrostatic potential into which ions are attracted along the field. The non-linear consequences of the instability are analysed in this paper. At

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Core-excited states of SF 6 probed with soft-x-ray femtosecond transient absorption of vibrational wave packets

A vibrational wavepacket in SF 6 is created by impulsive stimulated Raman scattering with a few-cycle infrared pulse and mapped simultaneously onto five sulfur core-excited states using table-top soft x-ray transient absorption spectroscopy between 170 to 200 eV. The femtosecond vibrations induce real-time energy shifts of the x-ray absorption, whose amplitude depend strongly on the nature of the core-excited state. The pump laser intensity is used to control the number of vibrational states in the superposition, thereby accessing core-excited levels for various extensions of the S-F stretching motion. This enables the determination of the relative core-level potential energy gradients for the symmetric stretching mode, in good agreement with TDDFT calculations. This experiment demonstrates a new means of characterizing core-excited potential energy curves.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Use of diffusive gradients in thin films (DGT) to measure potentially bioavailable metals in southeastern USA blackwater streams

Here, we used diffusive gradients in thin films (DGT) to measure potentially bioavailable metals in coastal plain streams in the southeastern USA that exhibited strong to moderate blackwater characteristics. Metals were partitioned into particulate metals, DGT-inert metals (i.e., colloidal and refractory organic complexes not accumulated by DGT), and DGT-labile metals (i.e., free metal ions, small inorganic complexes, and labile organic complexes). We also examined the influence of different DGT deployment times using data collected from the field and a follow-up laboratory study. The DGT-measured fraction of dissolved metals in the streams was 15% for Cd, 21% for Zn, 33% for Cu, 37% for Pb, and 98% for Mn. Metals bound to particulates predominated only for Pb. Most of the Cd, Pb, Zn, and Cu were associated with colloids, refractory organic complexes, or particles. Relatively small amounts were in free ion or labile complexes likely to be bioavailable through respiratory surfaces. Modeled concentrations of free and inorganically bound Cu and Pb were lower than the DGT fraction indicating that DGT accumulated some organically bound Cu and Pb that might not have been bioavailable. DGT-exposure times in excess of 5 days may have contributed to the accumulation of partly labile organic–metal complexes and were associated with substantial biofouling that caused metal uptake by DGT to depart from linearity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sheath constraints on turbulent magnetised plasmas

Abstract A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale λ D . In magnetised plasmas where the magnetic field is inclined at an oblique angle α with the target, some of the sheath potential variation occurs also on the ion sound gyroradius length scale ρ S cos ⁡ α , caused by finite ion gyro-orbit distortion and losses. We consider a collisionless and steady-state magnetised plasma sheath whose thickness l ms ∼ max ( λ D , ρ S cos ⁡ α ) is smaller than the characteristic length scale L of spatial fluctuations in the bulk plasma, such that the limit l ms / L → 0 is appropriate. Spatial structures are assumed to be magnetic field-aligned. In the case of small magnetic field angle α ∼ δ ≡ ρ S / L ≪ 1 , electric fields tangential to the target transport ions towards the target via E × B drifts at a rate comparable to the one from parallel streaming. A generalised form of the kinetic Bohm–Chodura criterion at the sheath entrance is derived by requiring that the sheath electric field have a monotonic spatial decay far from the target. The criterion depends on tangential gradients of potential and ion distribution function, with additional nontrivial conditions.

kinetic Bohm–Chodura condition↗
DOE OSTI2022

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Magnetic field-induced non-trivial electronic topology in Fe3− x GeTe2

The anomalous Hall, Nernst, and thermal Hall coefficients of the itinerant ferromagnet Fe3−xGeTe2 display anomalies upon cooling that are consistent with a topological transition that could induce deviations with respect to the Wiedemann–Franz (WF) law. This law has not yet been validated for the anomalous transport variables, with recent experimental studies yielding material-dependent results. Nevertheless, the anomalous Hall and thermal Hall coefficients of Fe3−xGeTe2 are found, within our experimental accuracy, to satisfy the WF law for magnetic fields μ0H applied along its c axis. Remarkably, large anomalous transport is also observed for μ0H||a axis with the field aligned along the gradient of the chemical potential generated by thermal gradients or electrical currents, a configuration that should not lead to their observation. These anomalous planar quantities are found to not scale with the component of the planar magnetization (M||), showing instead a sharp decrease beyond μ0H||= 4 T or the field required to align the magnetic moments along μ0H||. We argue that chiral spin structures associated with Bloch domain walls lead to a field-dependent spin chirality that produces a novel type of topological transport in the absence of interaction between the magnetic field and electrical or thermal currents. Locally chiral spin structures are captured by our Monte Carlo simulations incorporating small Dzyaloshinskii–Moriya and biquadratic exchange interactions. These observations reveal not only a new way to detect and expose topological excitations, but also a new configuration for heat conversion that expands the current technological horizon for thermoelectric energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Scaling Behavior of a Turbulent Kinetic Energy Closure Scheme for the Stably Stratified Atmosphere: A Steady-State Analysis

Abstract We present a turbulent kinetic energy (TKE) closure scheme for the stably stratified atmosphere in which the mixing lengths for momentum and heat are not parameterized in the same manner. The key difference is that, while the mixing length for heat tends toward the stability independent mixing length for momentum in neutrally stratified conditions, it tends toward one based on the Brunt–Väisälä time scale and square root of the TKE in the limit of large stability. This enables a unique steady-state solution for TKE to be obtained, which we demonstrate would otherwise be impossible if the mixing lengths were the same. Despite the model’s relative simplicity, it is shown to perform reasonably well with observational data from the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) using commonly employed model constants. Analyzing the scaling behavior of the nondimensional velocity and potential temperature gradients, or of the stability (correction) functions, reveals that for large stability the present model scales in the same manner as the first-order operational scheme of Viterbo et al. Alternatively, it appears as a blend of two cases of the TKE closure scheme of Baas et al. Critically, because a unique steady-state TKE can be obtained, the present model avoids the nonphysical behavior identified in one of the cases of Baas et al.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Nonlinear Response of Atmospheric Blocking to Early Winter Barents–Kara Seas Warming: An Idealized Model Study

Wintertime Ural blocking (UB) has been shown to play an important role in cold extremes over Eurasia, and thus it is useful to investigate the impact of warming over the Barents–Kara Seas (BKS) on the behavior of Ural blocking. In our study the response of UB to stepwise tropospheric warming over the BKS is examined using a dry dynamic core model. Nonlinear responses are found in the frequency and local persistence of UB. The frequency and local persistence of the UB increase with the strength of BKS warming in a less strong range and decrease with the further increase of BKS warming, which is linked to the UB propagation influenced by upstream background atmospheric circulation. For a weak BKS warming, the UB becomes more persistent due to its less westward movement associated with intensified upstream zonal wind and meridional potential vorticity gradient (PVy) in the North Atlantic mid-high latitudes, which corresponds to a negative height response over the North Atlantic high latitudes. When BKS warming is strong, a positive height response appears in the early winter stratosphere, and its subsequent downward propagation leads to a negative NAO response or increased Greenland blocking events, which reduces zonal wind and PVy in the high latitudes from North Atlantic to Europe, thus enhancing the westward propagation of UB and reducing its local persistence. The transition to the negative NAO phase and the retrogression of UB are not found when numerically suppressing the downward influence of weakened stratospheric polar vortex, suggesting a crucial role of the stratospheric pathway in nonlinear responses of UB to the early winter BKS warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Subtropical Jet in Reanalysis Data from STJ_PV

Subtropical jet position from a new method for locating the subtropical jet, called the tropopause gradient method. It is based on the peak gradient in potential temperature along the dynamic tropopause. This data has the identified subtropical jet latitude, level, and intensity across four different reanalysis products (CFSR-2, ERA-Interim, JRA-55, and MERRA-2), at both daily and monthly output frequency.

58 GEOSCIENCES↗
DOE OSTI2022

Microstructural Characterization of the Second High Fluence Baffle-Former Bolt Retrieved from a Westinghouse Two-loop Downflow Type PWR

As one of the pressurized water reactor (PWR) internal components, baffle-former bolts (BFBs) are subjected to significant mechanical stress and neutron irradiation from the reactor core during the plant operation. Over the long operation period, these conditions lead to potential degradation and reduced load-carrying capacity of the bolts. In support of evaluating long-term operational performance of materials used in core internal components, the Oak Ridge National Laboratory (ORNL), through the Department of Energy (DOE), Light Water Reactor Sustainability (LWRS) Program, Materials Research Pathway (MRP) has harvested two high fluence BFBs from a commercial Westinghouse two-loop downflow type PWR. The two bolts of interest, i.e. bolts # 4412 and 4416, were withdrawn from service in 2011 as part of a preventative replacement plan. No identification of cracking or potential damage was found for these bolts during their removal in 2011. However, the bolts required a lower torque for removal from the baffle structure than the original torque specified during installation. Irradiation displacement damage levels in the bolts range from 15 to 41 displacements per atom. The goal of this project is to perform detailed microstructural and mechanical property characterization of BFBs following in-service exposures. The information from these bolts will be integral to the LWRS program initiatives in evaluating end of life microstructure and properties. Furthermore, valuable data will be obtained that can be incorporated into model predictions of long-term irradiation behavior and compared to results obtained in high flux experimental reactor conditions. In this report, we present our latest study in FY22 on microstructural characterizations of the second high fluence baffle-former bolt, i.e., bolt # 4412. Analytical electron microscopy and atom probe tomography characterization were performed. The radiation-induced defects in the material add to the large wealth of knowledge for neutron-induced defects in 304/316 grades of stainless steels, specifically for radiation-induced precipitation after high fluence commercial PWR irradiation. The main findings are summarized as follows: 1) The cavity size was considerably larger in the bolt thread section than in the bolt head, with the bolt thread section having a bimodal distribution of cavities greater than ~6 nm in diameter and less than ~3 nm in diameter. The bolt head only had the small-sized cavities. In addition, there was a denuded zone of large cavities near grain boundaries in the thread section of the bolt. 2) Radiation-induced precipitation in the BFB #4412 was highly complex, with the volume fraction, size, and number density of Ni/Si and Cu-rich precipitates depending strongly on the radiation temperature/dose. In many cases, co-precipitates of adjoined clusters were found with Ni/Si-rich precipitates sandwiched between Cu-rich clusters and Mo/Cr/P-rich clusters. 3) Solute segregation out of solution was highest for most solutes in the thread section of the bolt #4412 with the exception of Cu, which experienced more separation out of solution into Curich clusters in the bolt head section. This highlights the difference in the mechanisms for precipitation of Ni/Si clusters, which have the Ni 3 Si phase composition, and precipitation of Cu-rich clusters. 4) There appear to be multiple simultaneous influences that affect the microstructural variation along the length of the bolt that overcomes the ~2X difference in irradiation dose between the bolt head and the bolt thread. The irradiation temperature, thermal/fast neutron ratio variation, potential strain gradient, and exposure to PWR coolant water that each section of the bolt sees may have more influence on the microstructural evolution than the total irradiation dose. The bolt thread and shank, with higher temperature, higher relative fast neutron flux, higher strain, and exposure to coolant but lower dose, underwent more enhanced cavity formation, precipitation, and solute segregation than the bolt head section.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

The characteristics of atmospheric boundary layer height over the Arctic Ocean during MOSAiC

The important roles that the atmospheric boundary layer (ABL) plays in the central Arctic climate system have been recognized, but the atmospheric boundary layer height (ABLH), defined as the layer of continuous turbulence adjacent to the surface, has rarely been investigated. Using a year-round radiosonde dataset during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, we improve a Richardson-number-based algorithm that takes cloud effects into consideration and subsequently analyze the characteristics and variability of the ABLH over the Arctic Ocean. The results reveal that the annual cycle is clearly characterized by a distinct peak in May and two respective minima in January and July. This annual variation in the ABLH is primarily controlled by the evolution of the ABL thermal structure. Temperature inversions in the winter and summer are intensified by seasonal radiative cooling and warm-air advection with the surface temperature constrained by melting, respectively, leading to the low ABLH at these times. Meteorological and turbulence variables also play a significant role in ABLH variation, including the near-surface potential temperature gradient, friction velocity, and turbulent kinetic energy (TKE) dissipation rate. In addition, the MOSAiC ABLH is more suppressed than the ABLH during the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment in the summer, which indicates that there is large variability in the Arctic ABL structure during the summer melting season.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Boundary Layer Height Variation over the Savannah River Site during a Heat Wave

With the number of heat waves expect to increase in both in their severity and frequency in the coming years, an insight into their impact on planetary boundary layer (PBL) heights is essential for one development of future emergency and safety plans. Using High-Resolution Rapid Refresh (HRRR) atmospheric model output, an analysis is done on the effects a heat wave has on PBL heights and temperature inversions across the Savannah River Site (SRS). This study showed that due to its adjacency to anthropogenic land use, SRS experiences a non-negligible gradient in its PBL height, with higher heights to the northwest. The singular location PBL heigh measurement and forecasts for SRS do not account for this gradient, which has potential implications on pollutant and smoke dispersion onsite. During a heat wave this gradient is present but PBL heights are less impacted than surrounding developed areas, possibly attributable to the woody wetlands that surround the Savannah River's banks. Additionally, the heat wave appears to strongly enhance the nocturnal temperature inversions that the site and surrounding areas experience. This is possibly due to increased stability from higher pressure conditions, which are often associated with heat waves. This influence on nocturnal inversions, however, is limited to their intensity, leaving their timing relatively unaltered. The impact a heat wave has on PBL heights and temperature inversions across SRS is seen to be more dependent on adjacent land use and features rather than the unique land use types onsite.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

On the Sampling-Based Computation of Nash Equilibria Under Uncertainty via the Nikaido–Isoda Function

We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be L 0 -Lipschitz continuous and convex, given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function defined using the Nikaido–Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected-gradient-response method, reliant on inexact resolution of a player-level best-response subproblem. Here, by deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of generated iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.

Nikaido-Isoda function↗
DOE OSTI2024

Plasma potential shaping using end-electrodes in the Large Plasma Device

We perform experiments in the Large Plasma Device at the University of California, Los Angeles, studying how different end-electrode biasing schemes modify the radial potential profile in the machine. We impose biasing profiles of different polarities and gradient signs on a set of five concentric electrodes placed 12 m downstream from the plasma source. We find that imposing concave-down profiles (negative potential radial gradient) on the electrodes creates radial potential profiles halfway up the plasma column that are comparable to those imposed on the electrodes and a few electron temperature in height, regardless of the biasing polarity. On the other hand, imposing concave-up profiles (positive potential radial gradient) leads to non-monotonic radial potential profiles. This observation can be explained by the current drawn through the electrodes and the parallel plasma resistivity, highlighting their important role in controlling the rotation of plasma. Concave-down plasma potential profiles, obtained by drawing electrons on the axis, are predicted to drive azimuthal drift velocities that can approach significant fractions of the ion sound speed in the central region of the plasma column.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

CODE-TO-CODE BENCHMARK STUDY FOR THERMAL STRESS MODELING AND PRELIMINARY ANALYSIS OF THE HIGH-TEMPERATURE SINGLE HEAT PIPE EXPERIMENT

Microreactors are very small nuclear reactors with typical thermal-energy output of up to 20 MWth. The microreactor concept is gaining more and more attention for safe, robust, and reliable supply of electricity for remote locations and to meet industrial process heat needs. Heat-pipe-cooled microreactor is one of the microreactors being investigated at Idaho National Laboratory. In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios for the heat-pipe-cooled microreactor designs. This paper discusses finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory, using two commercial software packages, Abaqus and STAR-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the locations where the largest temperature gradients appeared between heat pipe and electric heater. The temperature fields have good agreements between Abaqus and STAR-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment

In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Seawater pH reconstruction using boron isotopes in multiple planktonic foraminifera species with different depth habitats and their potential to constrain pH and <i>p</i>CO<sub>2</sub> gradients

Boron isotope systematics of planktonic foraminifera from core-top sediments and culture experiments have been studied to investigate the sensitivity ofδ11B of calcite tests to seawater pH. However, our knowledge of the relationship between δ11B and pH remains incomplete for many taxa. Thus, to expand the potential scope of application of this proxy, we report δ11B data for seven different species of planktonic foraminifera from sediment core tops. We utilize a method for the measurement of small samples of foraminifera and calculate the δ11B-calcite sensitivity to pH for Globigerinoides ruber, Trilobus sacculifer (sacc or without sacc), Orbulina universa, Pulleniatina obliquiloculata, Neogloboquadrina dutertrei, Globorotalia menardii, and Globorotalia tumida, including for unstudied core tops and species. These taxa have diverse ecological preferences and are from sites that span a range of oceanographic regimes, including some that are in regions of air–sea equilibrium and others that are out of equilibrium with the atmosphere. The sensitivity of δ11Bcarbonate to δ11Bborate (e.g., Δδ11Bcarbonate/Δδ11Bborate) in core tops is consistent with previous studies for T. sacculifer and G. ruber and close to unity for N. dutertrei, O. universa, and combined deep-dwelling species. Deep-dwelling species closely follow the core-top calibration fo rO. universa, which is attributed to respiration-driven microenvironments likely caused by light limitation and/or symbiont–host interactions. Our data support the premise that utilizing boron isotope measurements of multiple species within a sediment core can be utilized to constrain vertical profiles of pH andpCO2 at sites spanning different oceanic regimes, thereby constraining changes in vertical pH gradients and yielding insights into the past behavior of the oceanic carbon pumps.

54 ENVIRONMENTAL SCIENCES↗