Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple cycles”

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 595 records · Page 33

Electrolyte Oxidation Pathways in Lithium-Ion Batteries

The mitigation of decomposition reactions of lithium-ion battery electrolyte solutions is of critical importance in controlling device lifetime and performance. However, due to the complexity of the system, exacerbated by the diverse set of electrolyte compositions, electrode materials, and operating parameters, a clear understanding of the key chemical mechanisms remains elusive. In this work, operando pressure measurements, solution NMR, and electrochemical methods were combined to study electrolyte oxidation and reduction at multiple cell voltages. Two-compartment LiCoO 2 /Li cells were cycled with a lithium-ion conducting glass–ceramic separator so that the species formed at each electrode could be identified separately and further reactions of these species at the opposite electrode prevented. One principal finding is that chemical oxidation (with an onset voltage of ~4.7 V vs Li/Li + for LiCoO 2 ), rather than electrochemical reaction, is the dominant decomposition process at the positive electrode surface in this system. This is ascribed to the well-known release of reactive oxygen at higher states-of-charge, indicating that reactions of the electrolyte at the positive electrode are intrinsically linked to surface reactivity of the active material. Soluble electrolyte decomposition products formed at both electrodes are characterized, and a detailed reaction scheme is constructed to rationalize the formation of the observed species. The insights on electrolyte decomposition through reactions with reactive oxygen species identified through this work have a direct impact on understanding and mitigating degradation in high-voltage/higher-energy-density LiCoO 2 -based cells, and more generally for cells containing nickel-containing cathode materials (e.g., LiNi x Mn y Co z O 2 ; NMCs), as they lose oxygen at lower operating voltages.

25 ENERGY STORAGE↗

Second harmonic generation for estimating state of charge of lithium-ion batteries

This study applied the nonlinear ultrasonic method, second harmonic generation, to precisely estimate the state of charge (SoC) in lithium-ion batteries. The second harmonic of the longitudinal wave is generated on a pouch cell battery at 5 MHz with a through-transmission setup. The relative nonlinear parameter β' is determined by analyzing the amplitudes at the fundamental and second harmonic frequencies. To enhance the nonlinear parameter's measurement accuracy, multiple excitation amplitudes are employed. Two separate charge/discharge tests (four-cycle and eight-cycle) are conducted on the battery at a rate of C/10. The nonlinear parameter is measured periodically during the charge/discharge process, and temperature compensation is applied to the measurement. The correlation curves between the nonlinear parameter and the actual SoC align well for the four-cycle and eight-cycle tests, and a robust linear relationship is observed for both correlation curves. A linear model and a second-order polynomial model are applied to fit the correlation using all data points from both tests. The two models are employed to validate the SoC prediction on a second battery by using a four-cycle test. Furthermore, the results indicate that both models can predict the SoC with an accuracy of approximately 3%, whereas the polynomial model demonstrates smaller errors in the regions near 0% and 100% SoC. Therefore, the nonlinear parameter β', measured through the second harmonic generation, can effectively predict lithium-ion battery SoC with an accuracy of less than 3%.

25 ENERGY STORAGE↗

A Fire Community Observatory: Interdisciplinary, AI-informed Post-Fire Rapid Response for Improved Water Cycle Science at Watershed Scale

Wildfire is an ecological disturbance that disrupts the hydrological cycle. In the past few years, a record number of multiple-and-compounding fires have occurred across urban-wildland gradients in the Western United States. Changes to watershed hydrological partitioning in response to fires (infiltration, runoff, evapotranspiration) presents unprecedented challenges to “Water-in-the-West” through negative impacts to water supply and its quality, and is a direct threat to downstream communities, groundwater, and drinking water supply infrastructure. While much work is being done to advance Artificial Intelligence and Machine Learning (AI/ML) use during fires for emergency response (i.e. predict fire movement, direct evacuations), significant potential exists to use AI/ML to address three scientific grand challenges that are rarely addressed in a convergent science context: 1) how to enhance the potential resiliency of a landscape before fire(s), 2) how to cost-effectively and optimally monitor watershed changes after fires, and 3) how to predict future hydrological and biogeochemical trajectories in fire-impacted watershed given climate change. This whitepaper addresses DOE Focal Area 1: Data acquisition and assimilation enabled by machine learning, AI, and advanced methods.

58 GEOSCIENCES↗

Integrating Carbon Capture, Utilization, & Sequestration into Chemical Pulp Mills

The U.S. pulp and paper industry presents a unique and largely untapped opportunity for large- scale carbon dioxide removal (CDR). Unlike most industrial sectors, pulp mills rely heavily on biomass, meaning that much of their carbon emissions originate from atmospheric CO₂ that was recently captured by plants. If this biogenic CO₂ can be captured and permanently stored, pulp mills can be transformed from carbon emitters into net carbon removal facilities. This project was motivated by that opportunity and aimed to develop and evaluate integrated, low-cost strategies for capturing, utilizing, and sequestering CO₂ within existing chemical pulping operations. The scope of this work focused on four complementary innovations designed to integrate seamlessly into kraft pulp mill infrastructure: (1) in situ CO₂ capture within the recovery cycle, (2) oxy-fuel retrofitting of the rotary lime kiln to produce a high-purity CO₂ stream, (3) ex situ CO₂ capture and mineralization using pulp mill residues (dregs, grits, and lime mud), and (4) beneficial reuse of these residues as mineral carbonate fertilizers. The project combined process modeling, laboratory experimentation, life cycle assessment (LCA), and field trials to evaluate the technical feasibility, economic viability, and environmental impact of these approaches. The results demonstrate that pulp mills can serve as effective platforms for carbon removal when equipped with integrated carbon capture systems. Process modeling showed that combining sodium spiking with oxy-fuel calcination significantly enhances CO₂ capture efficiency while reducing costs by up to 31% compared to conventional configurations. Experimental work further revealed that calcination behavior in high-CO₂ environments differs substantially from traditional systems, leading to the development of a new kinetic model that predicts reaction rates under these conditions. This model provides essential design guidance for next-generation decarbonized lime kilns. In parallel, the project demonstrated that alkaline mineral residues generated during pulping operations can be repurposed as a sustainable alternative to agricultural lime. Across a wide range of soils in the southeastern United States, these materials performed equivalently to commercial lime in adjusting soil pH while offering lower greenhouse gas emissions and reduced cost. Field and greenhouse studies confirmed that crop and tree growth responses were comparable, supporting their viability as a drop-in replacement. This co-product pathway provides a practical utilization strategy that offsets costs and improves overall system economics. A major contribution of this project is the first comprehensive life cycle assessment of carbon removal in pulp and paper systems across multiple system boundaries. Results show that retrofitted mills can achieve carbon removal efficiencies ranging from 12% to 92%, depending on how the system is defined. This finding highlights a critical issue in carbon accounting: reported performance is highly sensitive to methodological choices. By explicitly quantifying these differences, this work provides valuable guidance for policymakers, carbon registries, and project developers working to standardize carbon removal metrics. From a commercialization perspective, the technologies investigated in this project are well- aligned with existing industrial infrastructure, minimizing the need for entirely new facilities. 3 DE-EE0009413 Industry engagement throughout the project—including collaboration with pulp and paper companies, equipment manufacturers, and carbon removal developers—has accelerated the transition from research to deployment. Notably, a commercial developer is actively pursuing carbon capture projects at pulp mills in the southeastern United States and has cited this research as a contributing foundation. The emergence of voluntary carbon markets and long-term offtake agreements further strengthens the business case for implementation. The broader public benefits of this work are significant. By enabling large-scale carbon removal using existing industrial systems, this approach offers a near-term pathway to reduce atmospheric CO₂ concentrations while supporting domestic manufacturing and rural economies. The reuse of industrial residues as fertilizers reduces reliance on mined materials, lowers costs for farmers, and decreases environmental impacts associated with conventional lime production. In addition, the project has supported workforce development by training graduate students and researchers in carbon capture technologies, helping to build capacity in a critical area of national interest. In conclusion, this project demonstrates that integrated carbon capture, utilization, and sequestration in pulp mills is both technically feasible and economically promising. By combining process innovation, experimental validation, and systems-level analysis, the work advances the understanding of how biomass-based industries can contribute to climate mitigation. The findings provide a strong foundation for commercial deployment and offer a scalable solution for transforming a major U.S. industry into a source of durable carbon removal.

09 BIOMASS FUELS↗

Toward Urban Water Security: Broadening the Use of Machine Learning Methods for Mitigating Urban Water Hazards

Due to the complex interactions of human activity and the hydrological cycle, achieving urban water security requires comprehensive planning processes that address urban water hazards using a holistic approach. However, the effective implementation of such an approach requires the collection and curation of large amounts of disparate data, and reliable methods for modeling processes that may be co-evolutionary yet traditionally represented in non-integrable ways. In recent decades, many hydrological studies have utilized advanced machine learning and information technologies to approximate and predict physical processes, yet none have synthesized these methods into a comprehensive urban water security plan. In this paper, we review ways in which advanced machine learning techniques have been applied to specific aspects of the hydrological cycle and discuss their potential applications for addressing challenges in mitigating multiple water hazards over urban areas. We also describe a vision that integrates these machine learning applications into a comprehensive watershed-to-community planning workflow for smart-cities management of urban water resources.

54 ENVIRONMENTAL SCIENCES↗

Generating target system specifications from a domain model using CLIPS

The quest for reuse in software engineering is still being pursued and researchers are actively investigating the domain modeling approach to software construction. There are several domain modeling efforts reported in the literature and they all agree that the components that are generated from domain modeling are more conducive to reuse. Once a domain model is created, several target systems can be generated by tailoring the domain model or by evolving the domain model and then tailoring it according to the specified requirements. This paper presents the Evolutionary Domain Life Cycle (EDLC) paradigm in which a domain model is created using multiple views, namely, aggregation hierarchy, generalization/specialization hierarchies, object communication diagrams and state transition diagrams. The architecture of the Knowledge Based Requirements Elicitation Tool (KBRET) which is used to generate target system specifications is also presented. The preliminary version of KBRET is implemented in the C Language Integrated Production System (CLIPS).

Sugumaran, Vijayan↗

Covariability in the Monthly Mean Convective and Radiative Diurnal Cycles in the Amazon

The diurnal cycle of convective clouds greatly influences the radiative energy balance in convectively active regions of Earth, through both direct presence, and the production of anvil and stratiform clouds. Previous studies show that the frequency and properties of convective clouds can vary on monthly timescales as a result of variability in the monthly mean atmospheric state. Furthermore, the radiative budget in convectively active regions also varies by up to 7 Wm-2 in convectively active regions. These facts suggest that convective clouds connect atmospheric state variability and radiation variability beyond clear sky effects alone. Previous research has identified monthly covariability between the diurnal cycle of CERES-observed top-of-atmosphere radiative fluxes and multiple atmospheric state variables from reanalysis over the Amazon region. ASVs that enhance (reduce) deep convection, such as CAPE (LTS), tend to shift the daily OLR and cloud albedo maxima earlier (later) in the day by 2-3 hr. We first test the analysis method using multiple reanalysis products for both the dry and wet seasons to further investigate the robustness of the preliminary results. We then use CloudSat data as an independent cloud observing system to further evaluate the relationships of cloud properties to variability in radiation and atmospheric states. While CERES can decompose OLR variability into clear sky and cloud effects, it cannot determine what variability in cloud properties lead to variability in the radiative cloud effects. Cloud frequency, cloud top height, and cloud microphysics all contribute to the cloud radiative effect, all of which are observable by CloudSat. In addition, CloudSat can also observe the presence and variability of deep convective cores responsible for the production of anvil clouds. We use these capabilities to determine the covariability of convective cloud properties and the radiative diurnal cycle.

Dodson, Jason B.↗

Wintertime losses of ozone in high northern latitudes

Total column ozone data over the past 22 to 30 years from ground-based Dobson and filter ozonometer stations between 30 N and 80 N are analyzed for residual trends remaining after allowance for the known geophysical variations corresponding to: (1) the change of seasons; (2) the quasi-biennial oscillation; and (3) the 11-year solar sunspot cycle. Examination of the data from several ground stations between 45 and 55 N indicated a seasonal difference in the long-term ozone series, with statistically significant losses in several winter months. Accordingly, the data from individual stations were analyzed with multiple regression analysis, seeking trends on a monthly basis after allowance for the known geophysical cycles. Previous statistical analyses were conducted as tests of 1-D model calculations which do not show any differences with the seasons, and any trend toward change in ozone concentrations is expressed in a yearly trend without seasonal variation. Such a model is inappropriate for calculations with a data set which exhibits a pronounced tendency toward seasonal differences in the trends. Comparisons with model calculations then require 2-D models into which seasonal and latitudinal differences can easily be programmed.

Harris, Neil↗

Techno-Economic Evaluation of Electrified Vehicle Options in Drayage Fleets

The electrification of drayage fleets offers potential economic and operational benefits, but the financial viability of electrified vehicles remains sensitive to battery cost, energy price, and fleet usage patterns. While total cost of ownership (TCO) is a useful benchmark, fleet operators and investors are equally concerned with investment performance metrics such as payback period (PB) and Internal Rate of Return (IRR), which better reflect financial risks and investment return timelines. This study develops a unified techno-economic framework that jointly evaluates TCO, PB, and IRR to determine when electrified trucks become cost-effective alternatives to diesel trucks. Building on a previously developed cost modeling tool and using real-world telematics data from a Class 8 drayage fleet at the Port of Savannah, the analysis incorporates projected battery cost trajectories, electricity and diesel price trends, vehicle efficiency improvements, and multiple battery capacities. Parameter ranges reflect widely cited projections and observed drayage-duty-cycle variability. A surrogate-modeling method approximates economic performance across thousands of battery cost–electricity price combinations, enabling high-resolution identification of conditions that achieve TCO parity, acceptable PB thresholds, and target IRR levels. Additionally, the study estimates the evolving share of the fleet that can feasibly electrify over time under multiple economic metrics. This integrated framework offers a novel, data-driven approach to inform risk-aware decision-making for fleet electrification and supports investment planning under evolving cost and operational conditions.

Sun, Ruixiao [ORNL] (ORCID:0000000341768676)↗

Improving the Recyclability of Polymer Composites With Cellulose Nanofibrils

Cellulose nanofibers (CNFs) have been widely studied for their reinforcing potential in high-performance composites. While there are numerous publications on CNF-reinforced composites in a variety of polymer matrices, few have considered the recyclability of such thermoplastic composites and whether the incorporation of CNFs deteriorates or improves their performance upon reprocessing. In this study, two thermoplastic resins, poly(lactic acid) (PLA), and glycol-modified polyethylene terephthalate (PETg), were prepared with CNF reinforcement and thermomechanically recycled to investigate the effect of CNF inclusion on the composite properties after reprocessing as well as their effect on the composites’ number of useful life cycles. Additionally, changes in mechanical, thermal, rheological, molecular, and microstructural properties of the composites and/or base resins were monitored as a function of cycle numbers. As is typical, the polymers’ molecular weight and mechanical performance deteriorated with continued processing. However, the addition of spray dried CNF was found to better maintain the mechanical performance of both polymers throughout multiple recycling steps as compared to neat samples. For example, the tensile strength of PETg with 20 wt% CNF after 6 processing cycles was found to exceed that of virgin neat PETg, and higher loadings of CNF were found to preserve a higher yield strength during multiple rounds of reprocessing compared to PETg composites with lower CNF loadings. Ultimately this study indicates that the addition of CNF to some thermoplastic materials can increase both their sustainability by offsetting the use of high-embodied energy resins and their circularity by enabling performance retention over more use cycles.

42 ENGINEERING↗

Torsional and biaxial (tension-torsion) fatigue damage mechanisms in Waspaloy at room temperature

Strain controlled torsional and biaxial (tension-torsion) low cycle fatigue behavior of Waspaloy was studied at room temperature as a function of heat treatment. Biaxial tests were conducted under proportional and nonproportional cyclic conditions. The deformation behavior under these different cyclic conditions was evaluated by slip trace analysis. For this, a Schmidt-type factor was defined for multiaxial loading conditions, and it was shown that when the slip deformation is predominant, nonproportional cycles are more damaging than proportional or pure axial or torsional cycles. This was attributed to the fact that under nonproportional cyclic conditions, deformation was through multiple slip, as opposed to single slip for other loading conditions, which gave rise to increased hardening. The total life for a given test condition was found to be independent of heat treatment. This was interpreted as being due to the differences in the cycles to initiation and propagation of cracks.

Jayaraman, N.↗

Results of Application Demonstration of Autonomous Offshore Power System Supporting Multiple Mobile and Static Assets

Over the last decade, several technology trends have dominated the innovation cycle within the terrestrial economy. Electrification, cloud computing, decarbonization, digitization, and autonomy are representative examples of macro trends that have had a rapid and fundamental impact; an impact felt not only in user-facing technologies, such as electric vehicles and software-as-a-service, but also in the infrastructure required to deliver those technologies. However, the scale of the impact and rollout has not been felt equally between the terrestrial and ocean economies. Here, the innovation cycle in the ocean economy is significantly slower. What accounts for the difference and are there solutions to this issue?

16 TIDAL AND WAVE POWER↗

High-speed ultraviolet photometry of HD 60435

We present the first high-speed ultraviolet photometry of an oscillating Ap star, HD 60435. After removing known orbital effects related to the Hubble Space Telescope, we confirm the presence of a strong pulsation period at a frequency of 123.70 cycles per day. In addition, we find significant amplitude modulation of this frequency that we suggest could be the result of beating of multiple periodicities. In this context, we suggest evidence for the presence of four additional frequencies at nu = 120.56, 126.55, 149.49, and 221.03 cycles per day. Three of these frequencies correspond well to frequencies detected in optical observations of HD 60435. The fourth, at 149.49 cycles per day, if real, is a potentially new pulsation mode that has not been detected in ground-based observations of this star. The amplitude of the 123 cycles per day pulsation is significantly larger in the ultraviolet than it is in the blue.

Taylor, M.↗

Low-Latitude Ethane Rain on Titan

Cassini ISS observed multiple widespread changes in surface brightness in Titan's equatorial regions over the past three years. These brightness variations are attributed to rainfall from cloud systems that appear to form seasonally. Determining the composition of this rainfall is an important step in understanding the "methanological" cycle on Titan. I use data from Cassini VIMS to complete a spectroscopic investigation of multiple rain-wetted areas. I compute "before-and-after" spectral ratios of any areas that show either deposition or evaporation of rain. By comparing these spectral ratios to a model of liquid ethane, I find that the rain is most likely composed of liquid ethane. The spectrum of liquid ethane contains multiple absorption features that fall within the 2-micron and 5-micron spectral windows in Titan's atmosphere. I show that these features are visible in the spectra taken of Titan's surface and that they are characteristically different than those in the spectrum of liquid methane. Furthermore, just as ISS saw the surface brightness reverting to its original state after a period of time, I show that VIMS observations of later flybys show the surface composition in different stages of returning to its initial form.

Cassini↗

Double-multiple streamtube model for Darrieus in turbines

An analytical model is proposed for calculating the rotor performance and aerodynamic blade forces for Darrieus wind turbines with curved blades. The method of analysis uses a multiple-streamtube model, divided into two parts: one modeling the upstream half-cycle of the rotor and the other, the downstream half-cycle. The upwind and downwind components of the induced velocities at each level of the rotor were obtained using the principle of two actuator disks in tandem. Variation of the induced velocities in the two parts of the rotor produces larger forces in the upstream zone and smaller forces in the downstream zone. Comparisons of the overall rotor performance with previous methods and field test data show the important improvement obtained with the present model. The calculations were made using the computer code CARDAA developed at IREQ. The double-multiple streamtube model presented has two major advantages: it requires a much shorter computer time than the three-dimensional vortex model and is more accurate than multiple-streamtube model in predicting the aerodynamic blade loads.

Paraschivoiu, I.↗

Ways to Increase Launch Velocities of 2-Stage Gas Guns

The amount of space debris is rapidly increasing and the debris is distributed over a wide variety of orbits. Satellites, manned space vehicles and space stations will have to pay increasing attention to the dangers of impacts with space debris. Various armoring techniques (i.e., double or triple layer armor) will have to tested extensively to determine the most effective armor per unit weight. Intersecting near-earth orbits can lead to impact velocities up to 15 km/sec. Conventional two-stage light gas guns can launch intact, controlled-shape projectiles with a density of 1.2 gm/cc and length- to-diameter ratios of 0.5-1.0 at velocities up to 8-9 km/sec. Higher velocities (10-11 km/sec) can be obtained' for very light projectiles. The higher launch velocities tend to be very severe on the high pressure coupling and barrel of the gun and lead to short component lifetimes. Clearly, the ability to raise the launch velocity of a gun (for reasonably massive projectile shapes) from 8-9 km/sec to 11-13 km/sec (or higher), without reduction of component lifetimes, would have significant benefits. This would allow much better simulation of the higher velocity debris impacts as well as better simulation of high speed re-entry into planetary atmospheres. Several techniques for increasing the launcher muzzle velocity above 8-9 km/sec have been studied using CFD simulations and appear to offer the potential for significant gains. The first technique is to use multiple compressions, instead of a single compression, in the pump tube of the light gas gun. In a sense, this is a kind of pre-heating of the gas in the pump tube; other types of pre-heating have yielded disappointing results in the past. The dynamics of the multiple compression pump tube is very different, however, from the earlier techniques, where the pump tube was typically heated ohmically before the gun 2 cycle was started. In this paper, we present CFD calculations that show that significant increases in muzzle velocity can be obtained with multiple compressions in the pump tube. With a conventional two-stage gun, an important limitation to obtaining higher velocities is friction and heat transfer to the barrel, which typically has a length- to- diameter ratio of 200-400. These viscous losses greatly reduce the effectiveness of the regions of the barrel far removed from the second stage breech. We have studied computationally the effect of adding an additional breech (or breeches) along the barrel to reduce these viscous losses. Velocity increases from 6.5 to 7.2 km/sec have been obtained using the main breech and one additional breech. In these results, both breeches were operated with hydrogen, heated electrothermally. We have also studied a gun geometry where the main breech is operated in the conventional manner, using piston compression. The additional breech is operated either with electrothermal heating or heating by using a high explosive charge in a novel geometry. The latter option provides very effective compression, heating and acceleration of the hydrogen working gas and is fully reusable. Calculations are presented which show that very substantial increases in muzzle velocity can be obtained this way, without overstressing the projectile or the 'gun. The third technique studied is to add a section of ram accelerator tube after the barrel to further accelerate the projectile. The ram accelerator used here is not the conventional premixed gas ram accelerator, but a new technique using high explosive as the energy source and pure hydrogen as the working gas in a geometry which can be made fully reusable. Preliminary results with this new rain accelerator geometry were presented and showed that stable ram accelerator drive can be established. Herein, detailed calculations axe presented which show that substantial velocity increases can be obtained using this ram accelerator technique in tandem with a conventional light gas gun.

Bogdanoff, David W.↗

Life Cycle Analysis of Thermoelectric Power Generation in the United States

In this article, the basics of performing a life cycle analysis of thermoelectric power generation are discussed, with three examples of life cycle greenhouse gas (GHG) balances for thermoelectric power generation forms: coal, natural gas, and nuclear. The final section compares multiple electricity generation methods in the United States. Results are presented on the basis of 1 megawatt-hour (MWh) of electricity delivered to the end user. Environmental life cycle results for greenhouse gas (GHG) emissions are presented as carbon dioxide equivalents, based on 100-year global warming potentials (GWPs) established by the 6th Assessment Report from the Intergovernmental Panel on Climate Change in 2021, commonly referred to as AR6 GWP values (IPCC, 2021). Additional details on other environmental life cycle impacts from non-GHG emissions to air, emissions to water, solid waste generation, and land use are available on the Department of Energy, National Energy Technology Laboratory’s Life Cycle Analysis website: www.netl.doe.gov/LCA.

Cutshaw, Ashley↗

Plio-Pleistocene time evolution of the 100-kyr cycle in marine paleoclimate records

A new time-series analysis technique (called envelope inversion), which is based on multiple taper spectral analysis, is used to address the question of an abrupt versus a gradual onset of the 100-kyr ice-age periodicity in the middle Pleistocene. Three long (greater than 2.6 m.y.) time series from Deep Sea Drilling Project site 607 (midlatitude Atlantic) and Ocean Drilling Program site 677 (equatorial Pacific) were analyzed using delta(O-18) records. Results do not yield compelling evidence for an abrupt change in the 100-kyr delta(O-18) signal. Rather, the results suggest that the 100-kyr delta(O-18) cycle is phase-locked with the 124-kyr eccentricity cycle some 300-400 kyr prior to its late Pleistocene growth in amplitude and phase lock with the 95-kyr eccentricity cycle.

Park, Jeffrey↗