Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high rate”

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

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗
DOE OSTI2025

Self-consistent hardness measurements spanning eleven decades of strain rate on a single material surface

A comprehensive understanding of material strength across strain rates typically requires the combination of results from different methods, which often vary in loading conditions and/or sampled volumes, leading to discrepancies in material behavior. This study presents a microindentation approach to measure hardness covering eleven orders of magnitude in strain rate, from quasi-static to phonon drag-dominated rates, on a single material surface under uniform testing conditions. By engineering the geometry of impactors used in laser induced particle impact testing, we extend the breadth of accessible strain rates, including multiple distinct rates exceeding 10⁵ s⁻¹. This self-consistent approach provides clearer insights into high-rate deformation mechanisms. Our results demonstrate a gradual increase in hardness with strain rate from quasi-static up to ultra-high rates, where a sharp upturn in hardness is observed.

Engineering↗
DOE OSTI2023

Prelithiation of Lithium Peroxide for Silicon Anode: Achieving a High Activation Rate

The use of lithium peroxide (Li2O2) as a cost-effective low-weight prelithiation cathode additive was successfully demonstrated. Through a series of studies on the chemical stability of Li2O2 and the activation process of Li2O2 on the cathode, we revealed that Li2O2 is more compatible with conventional electrolyte and cathode laminate slurry than lithium oxide. Due to the significantlysmaller size of commercial Li2O2, it can beused directly as a cathode additive. Moreover, the activation of Li2O2 on the cathode leads to the impedance growthof the cathode possibly resulting from the release of dioxygen andevacuation of Li2O2 inside the cathode. Withthe introduction of a new Li2O2 spread-coatingtechnique on the cathode, the capacity loss was suppressed. Si||NMCfull cells using Li2O2 spread-coated cathodedemonstrated a highly promising activation rate of Li2O2 and significantly enhanced specific capacity and cyclingstability compared to the uncoated full cells.

cathode additive↗
DOE OSTI2026

Intrinsically Conductive {pi}‑d Conjugated Layers with Co–N4 Active Sites for Efficient Nitrate Electrocatalysis and Zinc-Nitrate Batteries

Electrochemical synthesis of ammonia from nitrate has been extensively investigated as a potential alternative to the energy-intensive Haber-Bosch process. This approach not only operates under ambient conditions but also simultaneously removes nitrate contaminants while producing ammonia as a value-added product. However, the ongoing quest lies in designing an efficient electrocatalyst that achieves a high ammonia yield rate, high selectivity, and long-term stability. Herein, we report the outstanding performance of a Co–N4 coordinated π-d layered Co3(HITP)2 (HITP = 2,3,6,7,10,11-hexaiminotriphenylene) in nitrate electrocatalysis. The unique combination of abundant Co–N4 active sites and superior electrical conductivity enables significant electrocatalytic activity, delivering a maximum ammonia yield rate of 56.8 mg cm–2 h–1 at −0.8 V vs RHE and a Faradaic efficiency of ∼91% at −0.4 V vs RHE. Mechanistic analysis reveals that alkaline conditions accelerate water dissociation to generate adsorbed hydrogen intermediates (H*), which are utilized by Co–N4 sites to drive the stepwise hydrogenation of nitrate to ammonia while suppressing competing hydrogen evolution reaction (HER) pathways. Furthermore, integration of this catalyst into a zinc-nitrate battery resulted in a maximum power density of 5.3 mW cm–2 and an open-circuit potential of ∼1.45 V. These results highlight the potential of π-d conjugated Co–N4 materials as an efficient catalyst for both environmental remediation and energy conversion.

Namvar, shahrirar↗
DOE OSTI2026

Commissioning of the large-scale lead tungstate scintillating calorimeter

Here, we report on the installation and initial commissioning of a large-scale lead tungstate (PbWO4) scintillating crystal calorimeter developed for high-rate photon detection and precise energy measurement. The calorimeter comprises 1596 high-granularity, high-resolution scintillating crystals optimized for electromagnetic-shower detection over a wide energy range. Scintillation light from each crystal is read out by Hamamatsu R4125 photomultiplier tubes equipped with a custom voltage divider and front-end amplifier to ensure stable gain at high rates. All calorimeter modules were fabricated and characterized using a light-emitting diode–based optical test system prior to installation to verify uniformity and photodetector performance. After installation, the electromagnetic calorimeter was fully integrated into the experiment data acquisition and energy-based trigger systems. The optical response of the modules was equalized using the light-monitoring system, cosmic-ray muons, and photons from Compton-scattering events. Commissioning results demonstrate a reliably calibrated optical response and stable detector performance during the first run. These results validate the calorimeter design and commissioning methodology for large-scale scintillator-based photonic instrumentation.

Analog to digital converters↗
DOE OSTI2024

Laser-driven ion and electron acceleration from near-critical density gas targets: Towards high-repetition rate operation in the 1 PW, sub-100 fs laser interaction regime

Ion acceleration from gaseous targets driven by relativistic-intensity lasers was demonstrated as early as the late 1990s, yet most of the experiments conducted to date have involved picosecond-duration, Nd:glass lasers operating at low repetition rate. Here, we present measurements on the interaction of ultraintense ( ≈ 10 20 W cm − 2 , 1 PW ) , ultrashort ( ≈ 70 fs ) Ti:Sa laser pulses with near-critical ( ≈ 10 20 cm − 3 ) helium gas jets, a debris-free targetry with the potential for future compatibility with high ( ≈ 1 Hz ) repetition rate operation. We provide evidence of α particles being forward accelerated up to ≈ 2.7 − MeV energy with a total flux of ≈ 10 11 sr − 1 as integrated over > 0.1 − MeV energies and detected within a 0.5 − mrad solid angle. We also report on on-axis emission of relativistic electrons with an exponentially decaying spectrum characterized by a ≈ 10 − MeV slope, i.e., five times larger than the standard ponderomotive scaling. The total charge of these electrons with energy above 2 MeV is estimated to be of ≈ 1 nC , corresponding to ≈ 0.1 % of the laser drive energy. In addition, we observe the formation of a plasma channel, extending longitudinally across the gas density maximum and expanding radially with time. These results are well captured by large-scale particle-in-cell simulations, which reveal that the detected fast ions most likely originate from reflection off the rapidly expanding channel walls. The latter process is predicted to yield ion energies in the MeV range, which compare well with the measurements. Finally, direct laser acceleration is shown to be the dominant mechanism behind the observed electron energization. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neutron Detectors for Scattering Applications

The book covers, for the first time, the most common neutron detectors used in neutron scattering facilities and all of those in use at Oak Ridge National Lab. The first two chapters describe the facilities, instruments and the critical detector parameters needed by various instruments. The third chapter focuses on the key components of the 3He-based linear position sensitive detectors as well as on their electronics, which require particular attention to signal-processing and noise reduction. The fourth chapter is dedicated to the 3He alternatives where scintillators play a critical role. Finally, chapter 5 covers emerging neutron detection technologies including semiconductors, vacuum-based devices and their associated readouts, which will be required in the future for high rate and high-resolution neutron detectors.

Diawara, Yacouba↗
DOE OSTI2024

Modeling the Enceladus dust plume based on in situ measurements performed with the Cassini Cosmic Dust Analyzer

We analyzed data recorded by the Cosmic Dust Analyzer on board the Cassini spacecraft during Enceladus dust plume traversals. Our focus was on profiles of relative abundances of grains of different compositional types derived from mass spectra recorded with the Dust Analyzer subsystem during the Cassini flybys E5 and E17. The E5 profile, corresponding to a steep and fast traversal of the plume, has already been analyzed. In this paper, we included a second profile from the E17 flyby involving a nearly horizontal traversal of the south polar terrain at a significantly lower velocity. Additionally, we incorporated dust detection rates from the High Rate Detector subsystem during flybys E7 and E21. We derived grain size ranges in the different observational data sets and used these data to constrain parameters for a new dust plume model. This model was constructed using a mathematical description of dust ejection implemented in the software package DUDI. Further constraints included published velocities of gas ejection, positions of gas and dust jets, and the mass production rate of the plume. Our model employs two different types of sources: diffuse sources of dust ejected with a lower velocity and jets with a faster and more colimated emission. From our model, we derived dust mass production rates for different compositional grain types, amounting to at least 28 kg s –1 . Previously, salt-rich dust was believed to dominate the plume mass based on E5 data alone. The E17 profile shows a dominance of organic-enriched grains over the south polar terrain, a region not well constrained by E5 data. By including both E5 and E17 profiles, we find the salt-rich dust contribution to be at most 1% by mass. This revision also results from an improved understanding of grain masses of various compositional types that implies smaller sizes for salt-rich grains. Our new model can predict grain numbers and masses for future mission detectors during plume traversals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Highly Reversible Plating/Stripping of Porous Zinc Anodes for Multivalent Zinc Batteries

Zinc continues to garner immense interest due to its versatility as an anode material in several configurations utilizing either alkaline or mild-pH electrolytes. Current research on using mild-pH electrolytes has improved the rechargeability aspect of Zn-based batteries since Zn2+ is solely utilized for plating/stripping of Zn. Several studies have incorporated Zn metal foils, yet, dramatic improvements can be achieved by expressing Zn as a porous structure. Herein, we use a quasi-pulsed electrodeposition process to prepare a conformal Zn coating onto 3D porous copper foam. By tuning the electrodeposition parameters, we achieved an optimal Zn coating that undergoes reversible plating/stripping when tested in symmetric Zn cells, which supported a low overpotential of ~60 mV for up to 100 cycles. We further investigated changes in the surface morphology by studying the Zn surface of both foil and 3D structure using scanning electron microscopy and X-ray micro-computed tomography. Both techniques showed that the Zn foil undergoes dramatic alterations at the surface, which results in inhomogeneous deposition of Zn, whereas the 3D form exhibited minimal changes. Lastly, we paired both Zn foil and 3D Zn with vanadium oxide and demonstrated that the porous structure supports high rate capability and high specific capacity.

25 ENERGY STORAGE↗
DOE OSTI2021

Hydrogen isotope fractionation using graphene and related 2-D materials. Final report

This report provides a final technical report and listing of products from the research on a research project from Clemson university, focused on use of two-dimensional materials within membranes to accomplish hydrogen ion separations. The project used CVD graphene and related 2D materials to accomplish hydrogen / deuterium / tritium (H/D/T) separation during water electrolysis in proton-exchange-membrane (PEM)-style cells. 2D materials were incorporated into PEM cells at varying places with the aim of accomplishing isotopic sieving through 2D materials by selective hydron transport at high rates and with high selectivity. The project work scope included studies on H/D separation with small-scale cells at Clemson, collaborative work between Clemson and SRNL on computational modeling and tritium separations, and issues associated with scaleup.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

CMS High Granularity Calorimeter ECON-D ASIC overview and radiation testing results

The Compact Muon Solenoid (CMS) Experiment’s High Granularity Calorimeter (HGCAL) upgrade replaces the CMS electromagnetic and hadronic endcap calorimeters in preparation for the high-rate and high-radiation environment of the High Luminosity LHC. To effectively use the over 6 million channels of this highly-segmented “imaging” calorimeter, CMS is pioneering very front-end data compression with the Endcap Concentrator (ECON) ASICs – the ECON-T for the trigger path and the ECON-D for the data path. These 65 nm CMOS ASICs are radiation tolerant (200 Mrad) and low-power (< 2.5 mW/channel). In June 2023, we received the first full-functionality prototype of the data path concentrator ASIC, the ECON-D-P1. This talk will present an overview of the ECON-D-P1, summarize functionality and system testing, and present results from both Total Ionizing Dose (TID) and Single Event Effect (SEE) testing campaigns completed in summer 2023, validating the ECON-D radiation tolerant performance.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

28nm front end ASIC and 12” LGADs for 3D integration

The 3DIntSenS Collaboration—a joint effort between SLAC, Fermilab, and LLNL—is developing enabling technologies for next-generation radiation imaging detectors that combine ultra-fine spatial resolution (≈10 μm) with precision timing (<20 ps), while maintaining low power <1 W/cm2 and high data throughput. The approach leverages 3D integration between advanced CMOS readout ASICs and finely pixelated LGAD sensors to achieve the performance and scalability required for large-area, high-rate applications. High-granularity, precision-timing detectors are essential for scientific advances in HEP, NP, BES, and FES, but widespread adoption is limited by the cost and complexity of 3D integration. To close this gap, the collaboration is developing LGAD sensors compatible with 12-inch commercial CMOS processes, enabling cost-effective integration with high-performance ASICs under development. We present the design and results from a 28 nm CMOS ASIC prototype, including a low-jitter front end, and in-pixel TDC demonstrating sub-10 ps timing resolution. We also report on the co-design and characterization of reticle-scale LGAD sensors with 50 μm and 100 μm pixels and introduce the next 10k-pixel ASIC designed for full 3D integration. These advances represent a critical step toward scalable, high-resolution radiation imaging systems for future scientific instrumentation.

England, Troy [Fermilab] (ORCID:0000000154405255)↗
DOE OSTI2026

Two 28-nm front-end ASICs for ultra-fine spatial resolution and precision timing to be 3D integrated with 12 LGADs

The 3DIntSenS Collaboration—a joint effort between SLAC, Fermilab, and LLNL—is developing enabling technologies for next-generation radiation imaging detectors that combine ultra-fine spatial resolution (about 10 µm) with precision timing (<20 ps), while maintaining low power <1 W/cm2 and high data throughput. The approach leverages 3D integration between advanced CMOS readout ASICs and finely pixelated LGAD sensors to achieve the performance and scalability required for large-area, high-rate applications. High-granularity, precision-timing detectors are essential for scientific advances in HEP, NP, BES, and FES, but widespread adoption is limited by the cost and complexity of 3D integration. To close this gap, the collaboration is developing LGAD sensors compatible with 12-inch commercial CMOS processes, enabling cost-effective integration with high-performance ASICs under development. We present two 28 nm CMOS ASIC prototypes, including a low-jitter front end, and in-pixel TDC demonstrating sub-10 ps timing resolution. These advances represent a critical step toward scalable, high-resolution radiation imaging systems for future scientific instrumentation.

England, Troy [Fermilab] (ORCID:0000000154405255)↗
DOE OSTI2026

Two 28-nm front-end ASICs for ultra-fine spatial resolution and precision timing to be 3D integrated with 12 LGADs

The 3DIntSenS Collaboration—a joint effort between SLAC, Fermilab, and LLNL—is developing enabling technologies for next-generation radiation imaging detectors that combine ultra-fine spatial resolution (about 10 µm) with precision timing (<20 ps), while maintaining low power <1 W/cm2 and high data throughput. The approach leverages 3D integration between advanced CMOS readout ASICs and finely pixelated LGAD sensors to achieve the performance and scalability required for large-area, high-rate applications. High-granularity, precision-timing detectors are essential for scientific advances in HEP, NP, BES, and FES, but widespread adoption is limited by the cost and complexity of 3D integration. To close this gap, the collaboration is developing LGAD sensors compatible with 12-inch commercial CMOS processes, enabling cost-effective integration with high-performance ASICs under development. We present two 28 nm CMOS ASIC prototypes, including a low-jitter front end, and in-pixel TDC demonstrating sub-10 ps timing resolution. These advances represent a critical step toward scalable, high-resolution radiation imaging systems for future scientific instrumentation.

England, Troy [Fermilab] (ORCID:0000000154405255)↗
DOE OSTI2023

Rapid Fabrication of High Frame Rate Multichannel FTIR Spectrometers

Spectrally resolved signals in the short- to mid-wave infrared (SWIR/MWIR) bands at high-temporal resolution are critical for many national security remote sensing missions. Currently available off the shelf technology can achieve either high temporal resolution or high spectral resolution, but rugged instruments that can achieve both simultaneously remain mostly in the realm of one-off R&D projects. This report documents efforts to demonstrate a new technique for designing and building high resolution, high framerate multichannel FTIR (MC-FTIR) spectrometers that operate in the SWIR/MWIR bands. The core optical element in a MC-FTIR spectrometer is an array of statically-tuned lamellar grating interferometers (LGI). In the original MC-FTIR work these arrays were fabricated using a synchrotron x-ray lithography method. We proposed to instead fabricate these LGI arrays using multiphoton lithography (MPL), a 3D printing technique that can fabricate meso-scale structures with sub-micron precision. Although we were able to fabricate LGI arrays of sufficient size using MPL, the realized optical surfaces had unsuitably high optical form errors, precluding their use in a fieldable instrument. Further advancement in MPL technology may eventually enable fabrication of interferometer-grade LGI arrays.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dates and rates of loess deposition and soil formation on the Snake River Plain, Idaho, USA

Loess-paleosol sequences in eolian deposits on the Snake River Plain, Idaho, western United States, preserve records of late Pleistocene to Holocene glacial-interglacial cycles. Here, we examine climate-driven changes in the timing and rate of loess accretion, soil formation, and pedogenic carbonate accumulation in the eastern and central Snake River Plain. High rates of loess deposition often correspond with dry, cold, and windy conditions and/or high sediment supply, while soil formation indicates landscape stability. Optically stimulated luminescence dating of four loess-soil sequences provides a record of the timing of loess deposition, and radiocarbon dating ( 14 C) of soil inorganic carbon provides ages for soil formation and pedogenic carbonate precipitation. On the central Snake River Plain, marine isotope stage (MIS) 3 is marked by loess deposition between ca. 50 ka and 33 ka, and pedogenic carbonate accumulation between ca. 42 ka calibrated years before present (cal. B.P.) and 31 ka cal. B.P. During MIS 2, loess deposition between ca. 21 ka and 15 ka is contemporaneous with pedogenic carbonate precipitation. MIS 2 is characterized by diffuse pedogenic carbonate formation and rapid loess accumulation. Both the eastern and central Snake River Plain study sites experienced high rates of loess deposition during MIS 2; however, the eastern Snake River Plain had almost three times faster rates of loess accumulation (0.56−0.63 m/k.y. from ca. 25 ka to 17 ka) as compared to the central Snake River Plain (0.18−0.26 m/k.y. from ca. 21 ka to 15 ka). During MIS 1, loess deposition on the central Snake River Plain ca. 11 ka was coincident with pedogenic carbonate formation ca. 10 ka cal. B.P. MIS 3 and MIS 1 are characterized by A horizons and well-developed Bk horizons, indicating soil stability and lower rates of loess accumulation during these intervals. Finally, this study highlights the role of late Quaternary−scale climatic variability on rates and amounts of soil inorganic carbon storage.

Geosciences↗