SEARCH · Engineering Papers

Results for “large deviations”

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.

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 30

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Modeling and Optimizing Pumped Storage in a Multi-stage Large Scale Electricity Market under Portfolio Evolution

To leverage the fast-ramping capability of resources to provide great value to the grid, electricity system operators such as the Midcontinent Independent System Operator (MISO) continue to evolve their approaches for integrating energy storage resources, including pumpedstorage hydro (PSH), into the electricity markets. However, new challenges arise in modeling and optimizing these energy-limited resources across multiple market clearing processes and planning studies with uncertainties and imperfect information. For instance, current market practices of PSH owners specifying pumping/generating hours can result in sub-optimal generation dispatch. Letting grid operators optimize PSH with the consideration of multiple operating modes and energy limitation constraints can potentially bring economic benefits to both the system and the PSH owners. However, in multi-stage clearing process of electricity markets, utilizing the PSH flexibility to deal with realized uncertainties can cause deviation in the multi-stage scheduling processes. The resulting financial risks from the schedule deviation may not be acceptable to PSH owners. In addition, to effectively utilize this energy limited resource, the state of charge (SOC) constraints of PSH needs to be continuously optimized and the marginal cost of deviation need to reflect the expected cost to purchase or sell energy at future times to compensate for deviations. This project aims to develop a prototype enhanced PSH model and improved price signals in the multi-stage market clearing process with proper consideration of the unique characteristics of PSH, in order to better align underlying PSH capabilities with evolving grid needs, particularly including the needs for more frequent and larger cycling to manage variability and uncertainty from renewables.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Evidence of off-shell Higgs boson production from ZZ leptonic decay channels and constraints on its total width with the ATLAS detector

This Letter reports on a search for off-shell production of the Higgs boson using 139 fb —1 of $pp$ collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS detector at the Large Hadron Collider. The signature is a pair of $\overline{Z}$ bosons, with contributions from both the production and subsequent decay of a virtual Higgs boson and the interference of that process with other processes. The two observable final states are $\overline{Z}$ $\overline{Z}$ → 4$\ell$ and $\overline{Z}$ $\overline{Z}$ → 2$\ell$2$ν$ with $\ell$ = $e$ or $μ$. In the $\overline{Z}$ $\overline{Z}$ → 4$\ell$ final state, a dense Neural Network is used to enhance analysis sensitivity with respect to matrix element-based discrimination. The backgroundonly hypothesis is rejected with an observed (expected) significance of 3.3 (2.2) standard deviations, representing experimental evidence for off-shell Higgs boson production. Assuming that no new particles enter the production of the virtual Higgs boson, its total width can be deduced from the measurement of its off-shell production cross-section. The measured total width of the Higgs boson is ${4.5}_{—2.5}^{+3.3}$ MeV, and the observed (expected) upper limit on the total width is found to be 10.5 (10.9) MeV at 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Predicting the neutrinoless double-β-decay matrix element of 136 Xe using a statistical approach

Calculation of the nuclear matrix elements (NMEs) for double-β decay is of paramount importance for guiding experiments and for analyzing and interpreting the experimental data, especially for the search of the neutrinoless double β decay mode (0v⁢ββ). However, there are currently still large differences between the NME values calculated by different methods, hence a quantification of their uncertainties is very much required. Here, in this paper, we propose a statistical analysis of 0⁢v⁢ββ NME for the 136 Xe isotope, based on the interacting shell model, but using three independent effective Hamiltonians, emphasizing the range of the NMEs' most probable values and its correlations with observables that can be obtained from the existing nuclear data. Consequently, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (1.55–2.65) at 90% confidence level, with a mean value of 1.99 and a standard deviation of 0.37.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Let Each Quantum Bit Choose Its Basis Gates

Near-term quantum computers are primarily limited by errors in quantum operations (or gates) between two quantum bits (or qubits). A physical machine typically provides a set of basis gates that include primitive 2-qubit (2Q) and 1-qubit (1Q) gates that can be implemented in a given technology. 2Q entangling gates, coupled with some 1Q gates, allow for universal quantum computation. In superconducting technologies, the current state of the art is to implement the same 2Q gate between every pair of qubits (typically an XX-or XY-type gate). This strict hardware uniformity requirement for 2Q gates in a large quantum computer has made scaling up a time and resource-intensive endeavor in the lab. We propose a radical idea – allow the 2Q basis gate(s) to differ between every pair of qubits, selecting the best entangling gates that can be calibrated between given pairs of qubits. This work aims to give quantum scientists the ability to run meaningful algorithms with qubit systems that are not perfectly uniform. Scientists will also be able to use a much broader variety of novel 2Q gates for quantum computing. We develop a theoretical framework for identifying good 2Q basis gates on “nonstandard” Cartan trajectories that deviate from “standard” trajectories like XX. We then introduce practical methods for calibration and compilation with nonstandard 2Q gates, and discuss possible ways to improve the compilation. To demonstrate our methods in a case study, we simulated both standard XY-type trajectories and faster, nonstandard trajectories using an entangling gate architecture with far-detuned transmon qubits. We identify efficient 2Q basis gates on these nonstandard trajectories and use them to compile a number of standard benchmark circuits such as QFT and QAOA. Furthermore, our results demonstrate an 8x improvement over the baseline 2Q gates with respect to speed and coherence-limited gate fidelity.

quantum computing↗

DOE OSTI2023

Deep Learning Parameterization of Vertical Wind Velocity Variability via Constrained Adversarial Training

Atmospheric models with typical resolution in the tenths of kilometers cannot resolve the dynamics of air parcel ascent, which varies on scales ranging from tens to hundreds of meters. Small-scale wind fluctuations are thus characterized by a subgrid distribution of vertical wind velocity W with standard deviation σ W . The parameterization of σ W is fundamental to the representation of aerosol–cloud interactions, yet it is poorly constrained. Using a novel deep learning technique, this work develops a new parameterization for σ W merging data from global storm-resolving model simulations, high-frequency retrievals of W , and climate reanalysis products. The parameterization reproduces the observed statistics of σ W and leverages learned physical relations from the model simulations to guide extrapolation beyond the observed domain. Incorporating observational data during the training phase was found to be critical for its performance. The parameterization can be applied online within large-scale atmospheric models, or offline using output from weather forecasting and reanalysis products.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

The Challenge of Achieving Site Closure of a Complex Groundwater Plume - 20531

Over the last three decades, DOE has made substantial progress in nearly every area of nuclear waste cleanup and has closed and transitioned 91 of the 107 sites in the EM complex. Several of the remaining groundwater plumes at Hanford, Los Alamos, Paducah, Oak Ridge, and Savannah River are complex and tend to have unique site or contaminant characteristics which provide an almost impossible challenge to remediate to regulatory levels. As a result, these complex sites will require a large investment to remediate. In order for DOE to achieve regulatory closure, a technically sound case must be developed and presented to regulators and stakeholders that clearly demonstrates that remedial goals have been or will be irreversibly achieved or that there is an alternate acceptable end-state. Five lines of evidence are required: - A conceptual site model that is as complete as possible, with the knowledge that it may continue to evolve with time; - Demonstration that trends in contaminant concentrations, or flux, versus time indicate that remedial goals will be achieved on schedule or within a reasonable time thereafter and that trends that deviate from this path are explainable by the conceptual model; - Demonstration that remobilization of attenuated contaminants is unlikely to occur in the future; - Explicit consideration of remedial contingencies, should deviations in the path toward remedial goals occur; - Development of an effective long-term monitoring plan to detect any conditions that might alter the path toward remedial goals. SRS has made significant progress in the remediation of contamination associated with the F-area and H-Area Seepage Basins, and it is currently developing a strategic closure path that addresses each of these lines of evidence and identifies both basic and applied scientific issues that will need to be addressed in order to obtain regulatory closure. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Acoustic focusing of microplastics in microfabricated and steel tube devices: An experimental study on the effects from particle size and medium density

The contamination of aquatic systems with microplastics is becoming a severe threat to living organisms. Microplastics found in the environment consist of a wide range of particle sizes thus give rise to several challenges to current purification methods. Acoustic focusing has been utilized as a particle concentrating method in biological and non-biological applications. In this study, we investigate the feasibility of using acoustic focusing technology to isolate microplastics in aqueous samples and find that a two-stage mechanism is essential to concentrate microplastics of wide size ranges of the same material. We demonstrate for the first time experimentally that when the particle size is comparable to the resonance wavelength, they deviate from typical pressure nodal focusing and show pressure antinodal focusing. Further, we have expanded microplastics isolation beyond microfluidic systems and demonstrate their extraction capabilities relatively at higher flow rates using microplastics of polystyrene, polyethylene, and polymethyl methacrylate ranging from 6 µm to 300 µm diameter and in devices fabricated with cylindrical steel tubes. We could extract 70% of microplastics smaller than 180 µm and 82% of microplastics greater than 180 µm in diameter from an aqueous sample at a flow rate of 12 mL/min. In conclusion, we demonstrate that small microplastics focus on acoustic pressure nodes regardless of the medium density. However, focusing of large microplastics of the same material depends on the medium density.

42 ENGINEERING↗

DOE OSTI2022

Calculation of Nuclear Reactor Cooling Tower Performance With Limited Data Streams

Monitoring of cooling tower performance in a nuclear reactor facility is necessary to ensure safe operation; however, instrumentation for measuring performance characteristics can be difficult to install and may malfunction or break down over long duration experiments. This paper describes employing a thermodynamic approach to quantify cooling tower performance, the Merkel model, which requires only five parameters, namely, inlet water temperature, outlet water temperature, liquid mass flowrate, gas mass flowrate, and wet bulb temperature. Using this model, a general method to determine cooling tower operation for a nuclear reactor was developed in situations when neither the outlet water temperature nor gas mass flowrate are available, the former being a critical piece of information to bound the Merkel integral. Furthermore, when multiple cooling tower cells are used in parallel (as would be in the case of large-scale cooling operations), only the average outlet temperature of the cooling system is used as feedback for fan speed control, increasing the difficulty of obtaining the outlet water temperature for each cell. To address these shortcomings, this paper describes a method to obtain individual cell outlet water temperatures for mechanical forced-air cooling towers via parametric analysis and optimization. In this method, the outlet water temperature for an individual cooling tower cell is acquired as a function of the liquid-to-gas ratio (L/G). Leveraging the tight tolerance on the average outlet water temperature, an error function is generated to describe the deviation of the parameterized L/G to the highly controlled average outlet temperature. The method was able to determine the gas flowrate at rated conditions to be within 3.9% from that obtained from the manufacturer’s specification, while the average error for the four individual cooling cell outlet water temperatures were 1.6 °C, -0.5 °C, -1.0 °C, and 0.3 °C.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks

Abstract Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α‐diversity. To calculate the PSS index, we then compare observed partner phylogenetic α‐diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host–parasite, avian seed‐dispersal, lichenized fungi–cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much‐needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss .

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Search for heavy resonances in final states with four leptons and missing transverse momentum or jets in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for a new heavy boson produced via gluon-fusion in the four-lepton channel with missing transverse momentum or jets is performed. The search uses proton-proton collision data equivalent to an integrated luminosity of 139 fb -1 at a centre-of-mass energy of 13 TeV collected by the ATLAS detector between 2015 and 2018 at the Large Hadron Collider. This study explores the decays of heavy bosons: R → SH and A → ZH, where R is a CP-even boson, A is a CP-odd boson, H is a CP-even boson, and S is considered to decay into invisible particles that are candidates for dark matter. In these processes, S → invisible and H → ZZ. The Z boson associated with the heavy scalar boson H decays into all decay channels of the Z boson. The mass range under consideration is 390–1300 (320–1300) GeV for the R (A) boson and 220–1000 GeV for the H boson. No significant deviation from the Standard Model backgrounds is observed. The results are interpreted as upper limits at a 95% confidence level on the cross-section times the branching ratio of the heavy resonances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Detailed study of quark-hadron duality in spin structure functions of the proton and neutron

Background: The response of hadrons, the bound states of the strong force (QCD), to external probes can be described in two different, complementary frameworks: as direct interactions with their fundamental constituents, quarks and gluons, or alternatively as elastic or inelastic coherent scattering that leaves the hadrons in their ground state or in one of their excited (resonance) states. The former picture emerges most clearly in hard processes with high momentum transfer, where the hadron response can be described by the perturbative expansion of QCD, while at lower energy and momentum transfers, the resonant excitations of the hadrons dominate the cross section. The overlap region between these two pictures, where both yield similar predictions, is referred to as quark-hadron duality and has been extensively studied in reactions involving unpolarized hadrons. Some limited information on this phenomenon also exists for polarized protons, deuterons, and 3 He nuclei, but not yet for neutrons. Purpose: In this paper, we present comprehensive and detailed results on the correspondence between the extrapolated deep inelastic structure function g 1 of both the proton and the neutron with the same quantity measured in the nucleon resonance region. Thanks to the fine binning and high precision of our data, and using a well-controlled perturbative QCD (pQCD) fit for the partonic prediction, we can make quantitative statements about the kinematic range of applicability of both local duality and global duality. Method: We use the most updated QCD global analysis results at high x from the Jefferson Lab Angular Momentum Collaboration to extrapolate the spin structure function g 1 into the nucleon resonance region and then integrate over various intervals in the scaling variable x. We compare the results with the large data set collected in the quark-hadron transition region by the CLAS Collaboration, including, for the first time, deconvoluted neutron data, integrated over the same intervals. We present this comparison as a function of the momentum transfer Q 2 . Results: We find that, depending on the integration interval and the minimum momentum transfer chosen, a clear transition to quark-hadron duality can be observed in both nucleon species. Furthermore, we show, for the first time, the approach to scaling behavior for g 1 measured in the resonance region at sufficiently high momentum transfer. Conclusions: Here, our results can be used to quantify the deviations from the applicability of pQCD for data taken at moderate energies and can help with extraction of quark distribution functions from such data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

How Do Uncertainties in Galaxy Formation Physics Impact Field-level Galaxy Bias?

Our ability to extract cosmological information from galaxy surveys is limited by uncertainties in the galaxy–dark matter halo relationship for a given galaxy population, which are governed by the intricacies of galaxy formation. To quantify these uncertainties, we examine quenched and star-forming galaxies using two distinct approaches to modeling galaxy formation: U$\scriptsize{NIVERSE}$M$\scriptsize{ACHINE}$, an empirical semianalytic model, and the I$\scriptsize{llustris}$TNG hydrodynamical simulation. We apply a second-order hybrid N-body perturbative bias expansion to each galaxy sample, enabling direct comparison of modeling approaches and revealing how uncertainties in the galaxy–halo connection affect bias parameters and non-Poisson noise across number densities and redshifts. Notably, we find that quenched and star-forming galaxies occupy distinct parts of the bias parameter space, and that the scatter induced from these different galaxy formation models is small when conditioned on similar selections of galaxies. We also detect the signature of assembly bias in our samples; this leads to small but significant deviations from analytic bias predictions, while assembly bias-removed samples match these predictions well. This work indicates that galaxy samples from a spectrum of reasonable, physically motivated models for galaxy formation give a relatively small range of field-level galaxy bias parameters. We estimate a set of priors from these models that should be useful in extracting cosmological constraints from luminous red galaxy– and emission line galaxy–like samples. Looking forward, careful estimates of the range of impacts of galaxy formation, for a given sample and cosmological analysis, will be an essential ingredient for extracting the most precise cosmological information from current and future large galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Hyperfine-to-rotational energy transfer in ultracold atom–molecule collisions of Rb and KRb

Energy transfer between different mechanical degrees of freedom in atom–molecule collisions has been studied and largely understood. However, systems involving spins remain less explored. Here, in this study, we directly observed energy transfer from atomic hyperfine to molecular rotation in the 87 Rb ($| {F}_{a},{M}_{{F}_{a}}\rangle =| 2,2\rangle$) + 40 K 87 Rb (X 1 Σ + , rotational state N = 0) ⟶ Rb ($| 1,1\rangle$) + KRb ( N = 0, 1, 2) collision with state-to-state precision. We also performed quantum scattering calculations that rigorously included the coupling between spin and rotational degrees of freedom at short range under the assumption of rigid-rotor KRb monomers moving along a single potential energy surface. The calculated product rotational state distribution deviates from the observations even after extensive tuning of the atom–molecule potential energy surface. In addition, our ab initio calculations indicate that spin–rotation coupling is enhanced close to a conical intersection that is energetically accessible at short range. This, together with the deviation, suggests that vibrational degrees of freedom and conical intersections play an important part in the coupling. Our observations confirm that spin is coupled to mechanical rotation at short range and establish a benchmark for future theoretical studies.

chemical physics↗

DOE OSTI2022

Boundary characterization using 3D mapping of geometrically necessary dislocations in AM Ta microstructure

Additive manufacturing (AM) of high strength metallic materials produces microstructures and defects that deviate substantially from those created through conventional manufacturing methods. It has recently been demonstrated that with application of extreme energy densities, a high-temperature refractory metal, tantalum (Ta), can be printed to a fully dense state with exceptionally low porosity. Here we employ a combination of three-dimensional electron backscattered diffraction (EBSD) TriBeam technology and crystallographic geometrically necessary dislocation (GND) theory to characterize the microstructure and defect boundaries of this unusual AM product. Two-dimensional (2D) EBSD and these three-dimensional (3D) measurements indicate that the microstructure of AM Ta is highly oriented <111> along the build direction, yet at the same time contains large crystallographic orientation gradients that span mm’s across the build. Crystallographic GND density analysis of the 3D microstructure reveals that highly misoriented subboundaries exist within this strongly textured microstructure that have large dislocation densities of 1×10 16 m –2 , just as large as those comprising the high-angle grain boundaries (HAGBs) in the same material. The 3D crystallographic GND density mapping reveals these subboundaries are part of a complex, finely spaced network that extends throughout the entire microstructure. Here, the orientation of these boundaries can be related to the scan strategy used during printing. TEM measurements corroborate an extremely high dislocation density at the microscale and indicate a cell-like dislocation network structure existing in the AM Ta at the sub-μm scale.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Self-interacting Dark Matter, Core Collapse, and the Galaxy–Galaxy Strong-lensing Discrepancy

Abstract Gravitational lensing by galaxy clusters has emerged as a powerful tool to probe the standard cold dark matter (CDM) paradigm of structure formation in the Universe. Despite the remarkable explanatory power of CDM on large scales, tensions with observations on small scales have been reported. Recent studies find that the observational cross section for galaxy–galaxy strong lensing (GGSL) in clusters exceeds the CDM prediction by more than an order of magnitude and persists even after rigorous examination of various possible systematics. We investigate the impact of modifying the internal structure of cluster dark matter subhalos on GGSL and report that altering the inner density profile, given by r γ , to steeper slopes with γ > 2.5 can alleviate the GGSL discrepancy. This is steeper than slopes obtained with the inclusion of the contribution of baryons to the inner regions of these subhalos. Deviating from the γ ~ 1.0 cusps that CDM predicts, these steeper slopes could arise in models of self-interacting dark matter undergoing core collapse. Our results motivate additional study of subhalo core collapse in dense cluster environments.

Dutra, Isaque (ORCID:0000000170404930)↗

DOE OSTI2024

Resilience of the slow component in timescale separated synchronized oscillators

Physiological networks are usually made of a large number of biological oscillators evolving on a multitude of different timescales. Phase oscillators are particularly useful in the modelling of the synchronization dynamics of such systems. If the coupling is strong enough compared to the heterogeneity of the internal parameters, synchronized states might emerge where phase oscillators start to behave coherently. Here, we focus on the case where synchronized oscillators are divided into a fast and a slow component so that the two subsets evolve on separated timescales. We assess the resilience of the slow component by, first, reducing the dynamics of the fast one using Mori-Zwanzig formalism. Second, we evaluate the variance of the phase deviations when the oscillators in the two components are subject to noise with possibly distinct correlation times. From the general expression for the variance, we consider specific network structures and show how the noise transmission between the fast and slow components is affected. Interestingly, we find that oscillators that are among the most robust when there is only a single timescale, might become the most vulnerable when the system undergoes a timescale separation. We also find that layered networks seem to be insensitive to such timescale separations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Resilience of the slow component in timescale-separated synchronized oscillators

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Multiscale Temporal Variability of the Global Air‐Sea CO 2 Flux Anomaly

Abstract The global air‐sea CO 2 flux (F) impacts and is impacted by a plethora of climate‐related processes operating at multiple time scales. In bulk mass transfer formulations, F is driven by physico‐ and bio‐chemical factors such as the air‐sea partial pressure difference (∆pCO 2 ), gas transfer velocity, sea surface temperature, and salinity–all varying at multiple time scales. To de‐convolve the impact of these factors on variability in F at different time scales, time‐resolved estimates of F were computed using a global data set assembled between 1988 and 2015. The F anomalies were defined as temporal deviations from the 28‐year time‐averaged value. Spectral analysis revealed four dominant timescales of variability in F–subseasonal, seasonal, interannual, and decadal with relative amplitude differences varying across regions. A second‐order Taylor series expansion was then conducted along these four timescales to separate drivers across differing regions. The analysis showed that on subseasonal timescales, wind speed variability explains some 66% of the global F anomaly and is the dominant driver. On seasonal, interannual, and decadal timescales, the ∆pCO 2 effect controlled by the ∆pCO 2 anomaly, explained much of the F anomaly. On decadal timescales, the F anomaly was almost entirely governed by the ∆pCO 2 effect with large contributions from high latitudes. The main drivers across timescales also dominate the regional F anomaly, particularly in the mid‐high latitude regions. Finally, the driver of the ∆pCO 2 effect was closely connected with the relative strength of atmospheric pCO 2 and the nonthermal component of oceanic pCO 2 anomaly associated with dissolved inorganic carbon and alkalinity.

Environmental Sciences & Ecology↗