Search NASA⌕ Search

SEARCH · Search NASA

Results for “Earth Impact”

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 OSTI2024

The Utility of Infrasound in Global Monitoring of Extraterrestrial Impacts: A Case Study of the 2008 July 23 Tajikistan Bolide

Abstract Among the various observational techniques used for the detection of large bolides on a global scale is a low-frequency sound known as infrasound. Infrasound, which is also one of the four sensing modalities used by the International Monitoring System, offers continuous global monitoring and can be leveraged for planetary defense. Infrasonic records can provide an additional dimension for event characterization and a distinct perspective that might not be available through any other observational method. This paper describes the infrasonic detection and characterization of the bolide that disintegrated over Tajikistan on 2008 July 23. This event was detected by two infrasound stations at distances of 1530 and 2130 km. Propagation paths to one of the stations were not predicted by the model, despite being clearly detected. The presence of the signal is attributed to the acoustic energy being trapped in a weak but leaky stratospheric AtmoSOFAR channel. The infrasound signal analysis indicates that the shock originated at the point of the main breakup at an altitude of 35 km. The primary mode of the shock production of the signal detected at the two stations was a spherical blast resulting from the main gross fragmentation episode. The energy estimate, based on the signal period, is 0.17–0.51 kt of TNT equivalent, suggesting a mass of 6.6–23.5 tons. The corresponding object radius, assuming the chondritic origin, was 0.78–1.18 m.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Machine learning reveals strong grid-scale dependence in the satellite N d –LWP relationship

The relationship between cloud droplet number concentration ( N d ) and liquid water path (LWP) is highly uncertain yet crucial for determining the impact of aerosol-cloud interactions (ACI) on Earth's radiation budget. The N d -LWP relationship is examined using a machine learning (ML) random forest model applied to five years of satellite data at grid resolutions ranging from 10° to 0.05° in 12 distinct regions. In the subtropics, the shape of the N d -LWP relationship switches from an inverted-V at 1° grid-resolution to an “M” shape at 0.1° resolution with decreased $\frac{\textrm{dln⁡LWP}}{\textrm{dln⁡}N_d}$ sensitivity. Tropical and midlatitude regions generally show a more positive sensitivity. Cloud sampling and filtering also influence this slope, wherein the exclusion of thin clouds, as commonly performed to reduce retrieval uncertainty, leads to strongly negative sensitivity across all regions. Precipitation is primarily responsible for driving the strength of the sensitivity, with strong positive slopes in raining clouds and negative and/or neutral responses found in non-raining clouds. A new method to compute radiative forcing from the ML model shows a robust Twomey radiative forcing across all regions and grid resolutions. However, LWP and cloud fraction adjustments to the radiative forcing, which are ∼50 % or smaller than the Twomey effect, decrease to negligible values with higher spatial resolution data. As Earth system models move toward higher spatial resolutions in the future, evaluating the LWP and CF adjustment contributions to the radiative forcing budget at these finer resolutions will be essential for evaluation and model development.

Aerosol-Cloud Interactions↗
DOE OSTI2024

Accurate assessment of land–atmosphere coupling in climate models requires high-frequency data output

Land–atmosphere (L–A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centered climate anomalies. Local L–A coupling (LoCo) metrics capture relevant L–A processes, highlighting the impact of soil and vegetation states on surface flux partitioning and the impact of surface fluxes on boundary layer (BL) growth and development and the entrainment of air above the BL. A primary goal of the Climate Process Team in the Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) project is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional Earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L–A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land–atmosphere interactions span timescales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability in behavioral regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here, we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyze, and better understand L–A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behavior. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant timescales, (ii) enable demonstration of inter-model spread in L–A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ultrafast structural response of shock‐compressed plagioclase

Abstract Meteor impacts can induce unique pressure‐dependent structural changes in minerals due to the propagation of shock waves. Plagioclase—ubiquitous throughout the Earth’s crust, extraterrestrial bodies, and meteorites—is commonly used for reconstructing the impact history and conditions of the parent bodies. However, there have been unresolved inconsistencies in the interpretation of shock transformations across previous studies: The pressure at which amorphization begins and the process by which it occurs is the subject of ongoing debate. Here, we utilize time‐resolved in situ X‐ray diffraction (XRD) to probe the phase transformation pathway of plagioclase during shock compression at a sub‐nanosecond timescale. Direct amorphization begins at pressures much lower than what was previously assumed, just above the Hugoniot elastic limit of 5 GPa, with full amorphization to a high‐density amorphous phase, observed at 32(10) GPa and 20 ns. Upon release, the material partially recrystallizes back into the original structure, demonstrating a memory effect.

58 GEOSCIENCES↗
DOE OSTI2026

Machine Learning Calibration of Groundwater Table Depth in ELM: Impact on Land Surface Hydrology and Land‐Atmosphere Fluxes

Accurate representation of groundwater table depth (GWTD) is crucial for simulating hydrological cycling in Earth system models (ESM). Nevertheless, there is a notable gap in the literature regarding the validation of GWTD simulations in ESMs and their subsequent impact on downstream hydrological components. This study explores the calibration of parameterization of global GWTD using machine learning within the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Despite achieving significant gains in simulating GWTD through calibration, offline ELM simulations unexpectedly show that these improvements do not translate to substantial enhancements in model performance for other key hydrological variables, including soil moisture (SM), runoff, groundwater contribution to runoff or base flow index (BFI), and evapotranspiration and its partitioning. The performance in SM and runoff was even degraded in some regions, while BFI was mostly overestimated. Although there is significant improvement in GWTD within the critical range of 1–5 m, where groundwater traditionally influences land surface energy fluxes, these improvements occurred mostly in humid areas where the impact of GWTD on surface processes is minimal. Although the impacts of model calibration are generally small in offline ELM simulations, coupled land-atmosphere simulations exhibit much stronger responses to GWTD calibration, highlighting the role of land-atmosphere feedbacks in Earth system modeling. These findings underscore the need for integrated calibration strategies that simultaneously optimize multiple hydrological variables. However, if a single-variable approach is necessary, it is crucial to establish clear priorities for calibration, identifying the most critical variables that have the greatest impact on overall model performance.

Fang, Yilin [Pacific Northwest National Laboratory↗
DOE OSTI2023

Grounding our Understanding of the Impacts of Boreal Forest Expansion on Shallow Cumulus Clouds with a Simple Modeling Framework

Abstract The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns. Significance Statement Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Hunga Tonga–Hunga Ha′apai Volcano Impact Model Observation Comparison (HTHH-MOC) project: experiment protocol and model descriptions

The 2022 Hunga volcanic eruption injected a significant amount of water vapor and a moderate amount of sulfur dioxide into the stratosphere, causing observable responses in the climate system. We have developed a model–observation comparison project to investigate the evolution of volcanic water and aerosols and their impacts on atmospheric dynamics, chemistry, and climate, using several state-of-the-art chemistry climate models. The project goals are (1) to evaluate the current chemistry–climate models to quantify their performance in comparison to observations and (2) to understand atmospheric responses in the Earth system after this exceptional event and investigate the potential impacts in the projected future. To achieve these goals, we designed specific experiments for direct comparisons to observations, for example from balloons and the Microwave Limb Sounder satellite instrument. Experiment 1 consists of two sets of free-running ensemble experiments from 2022 to 2031: one with fixed sea-surface temperatures and sea ice and one with coupled ocean. These experiments will help to understand the long-term evolution of water vapor and aerosols; quantify HTHH effects on stratospheric and mesospheric temperatures, dynamics, and transport; understand the impact of dynamic changes on ozone chemistry; quantify the net radiative forcings; and evaluate any surface climate impact. Experiment 2 is a nudged-run experiment from 2022 to 2023 using observed meteorology. To allow participation of more climate models with varying complexities of aerosol simulation, we include two sets of simulations in Experiment 2: Experiment 2a is designed for models with internally generated aerosol, while Experiment 2b is designed for models using prescribed aerosol surface area density. This experiment will help to analyze H 2 O and aerosol evolution, quantify the net radiative forcings, understand the impacts on mid-latitude and polar O 3 chemistry, and allow close comparisons with observations.

Zhu, Yunqian [Univ. of Colorado, Boulder, CO (Unit↗
DOE OSTI2026

Arctic shipping under global change: A case study of offshore oil exports

We explore impacts of sea ice thinning and evolutions in the energy sector on future use of the Northern Sea Route (NSR) versus the Suez Canal Route (SCR), using a case study of shipping oil extracted from the offshore Russian Arctic to China. We combine an integrated human-Earth system model with a shipping cost model to incorporate impacts on both oil production and shipping costs under internally consistent scenarios. We find that the NSR could become cost-competitive with the SCR as sea ice thickness declines, especially in an RCP8.5 scenario, due to decreasing fuel and icebreaker escort costs. In a global energy evolution scenario consistent with RCP2.6, high emissions costs on the longer SCR may outweigh the costs associated with thicker sea ice on the NSR. Our novel framework provides integrated projections of NSR shipping traffic driven by a specific commodity likely to be shipped through the Arctic.

Arctic↗
DOE OSTI2020

Global land use for 2015–2100 at 0.05° resolution under diverse socioeconomic and climate scenarios

Abstract Global future land use (LU) is an important input for Earth system models for projecting Earth system dynamics and is critical for many modeling studies on future global change. Here we generated a new global gridded LU dataset using the Global Change Analysis Model (GCAM) and a land use spatial downscaling model, named Demeter, under the five Shared Socioeconomic Pathways (SSPs) and four Representative Concentration Pathways (RCPs) scenarios. Compared to existing similar datasets, the presented dataset has a higher spatial resolution (0.05° × 0.05°) and spreads under a more comprehensive set of SSP-RCP scenarios (in total 15 scenarios), and considers uncertainties from the forcing climates. We compared our dataset with the Land Use Harmonization version 2 (LUH2) dataset and found our results are in general spatially consistent with LUH2. The presented dataset will be useful for global Earth system modeling studies, especially for the analysis of the impacts of land use and land cover change and socioeconomics, as well as the characterizing the uncertainties associated with these impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modeling the Air Quality Impacts of Future Energy Scenarios

Many greenhouse gas (GHG) emission reduction measures achieve simultaneous reductions in air pollutants. Human-Earth system models can estimate such emission changes in the energy system but using them in chemistry-transport models (CTMs) to study their air quality impacts involves resource-intensive emissions processing. This is greatly simplified by an emissions scaling approach linking state-level emissions estimated by a human-Earth system model to a CTM. A scenario continuing pre-2022 energy policy in the U.S. to 2050 shows widespread air quality improvements over the 2015 baseline from SO 2 and NO x emission reductions of 50–80% from electricity generation and light-duty vehicles. Scenarios of GHG mitigation and vehicle electrification at the state and national level add further benefits. However, PM 2.5 increases from increased use of wood heating and bioenergy suggest that additional PM 2.5 management may be needed when using biofuels. In conclusion, this approach helps assess multiple future energy scenarios efficiently without sacrificing chemical detail in the air quality simulations.

air quality↗
DOE OSTI2023

Barotropic tides in MPAS-Ocean (E3SM V2): impact of ice shelf cavities

Abstract. Oceanic tides are seldom represented in Earth system models (ESMs) owing to the need for high horizontal resolution to accurately represent the associated barotropic waves close to coasts. This paper presents results of tides implemented in the Model for Prediction Across Scales–Ocean or MPAS-Ocean, which is the ocean component within the U.S. Department of Energy developed Energy Exascale Earth System Model (E3SM). MPAS-Ocean circumvents the limitation of low resolution using unstructured global meshing. We are at this stage simulating the largest semidiurnal (M2, S2, N2) and diurnal (K1, O1) tidal constituents in a single-layer version of MPAS-O. First, we show that the tidal constituents calculated using MPAS-Ocean closely agree with the results of the global tidal prediction model TPXO8 when suitably tuned topographic wave drag and bottom drag coefficients are employed. Thereafter, we present the sensitivity of global tidal evolution due to the presence of Antarctic ice shelf cavities. The effect of ice shelves on the amplitude and phase of tidal constituents are presented. Lower values of complex errors (with respect to TPXO8 results) for the M2 tidal constituents are observed when the ice shelf is added in the simulations, with particularly strong improvement in the Southern Ocean. Our work points towards future research with varying Antarctic ice shelf geometries and sea ice coupling that might lead to better comparison and prediction of tides and thus better prediction of sea-level rise and also the future climate variability.

58 GEOSCIENCES↗
DOE OSTI2025

The DOE E3SM version 2.1: overview and assessment of the impacts of parameterized ocean submesoscales

The U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 2.1 builds on E3SMv2 with several changes, with the most notable being the addition of the Fox-Kemper et al. (2011) mixed-layer eddy parameterization. This parameterization captures the effect of finite-amplitude, mixed-layer eddies as an overturning streamfunction and has the primary function of restratification. Herein, we outline the changes to the mean climate state of E3SM that were introduced by the addition of this parameterization. Overall, the presence of the submesoscale parameterization improves the fidelity of the v2.1 simulation by reducing the ocean surface biases in the North Atlantic present in v2, as illustrated by changes in the climatological sea surface temperature and salinity and the Arctic sea-ice extent. Other impacts include a slight shoaling of the mixed-layer depths in the North Atlantic and a small improvement in the Atlantic Meridional Overturning Circulation (AMOC). We note that the expected shoaling due to the parameterization is regionally dependent in our coupled configuration. In addition, we investigate why the parameterization and its impacts on mixed-layer depth have little impact on the simulated AMOC: despite increased dense-water formation in the Norwegian Sea, only a small fraction of the water formed makes its way south into the North Atlantic basin. Version 2.1 also exhibits small improvements in the atmospheric climatology, with smaller biases in many notable quantities and modes of variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Automatic detection of impact craters on Al foils from the Stardust interstellar dust collector using convolutional neural networks

NASA's Stardust mission utilized a sample collector composed of aerogel and aluminum foil to return cometary and interstellar particles to Earth. Analysis of the aluminum foil begins with locating craters produced by hypervelocity impacts of cometary and interstellar dust. Interstellar dust craters are typically less than one micrometer in size and are sparsely distributed, making them difficult to find. In this paper, we describe a convolutional neural network based on the VGG16 architecture that achieves high specificity and sensitivity in locating impact craters in the Stardust interstellar collector foils. Finally, we evaluate its implications for current and future analyses of Stardust samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Effects of Realistic Impact-processed Asteroid Topology on Atmospheric Airbursts

When small asteroids (∼10–200 m) enter Earth’s atmosphere they can break up before hitting the ground, releasing a considerable amount of energy in what is known as an atmospheric airburst. To date, numerous studies have examined how an asteroid’s strength, size, shape, and entry characteristics affect the breakup process; but few have assessed the effects of heterogenous internal topologies. We present a two-step high-fidelity simulation approach, which we use to model the atmospheric breakup of asteroids with realistic impact-processed internal topologies and their associated ground effects. In the first step, we model the asteroid as it is impacted to generate an “impact-processed” topology. We then map the resulting fractures (i.e., damage field) from that simulation to the initial state of an atmospheric breakup simulation. For both phases, we use the smoothed particle hydrodynamics code Spheral++. As a demonstrative example, we apply our new simulation pipeline to a Tunguska-scale impactor (50 m diameter, 20 km s −1 entry velocity, and 45° angle) in a case study comparing impact-processed and fully intact asteroids in both planar 2D and full 3D. For sufficiently strong asteroids, we find that impact processing generally causes the asteroid to break up at higher altitudes. We also observe that impact-processed asteroids break up into smaller fragments than objects with equivalent properties that have not been impact-processed.

Astronomy and AstroPhysics↗
DOE OSTI2024

Fire–precipitation interactions amplify the quasi-biennial variability in fires over southern Mexico and Central America

Fires have great ecological, social, and economic impact. However, fire prediction and management remain challenges due to a limited understanding of their roles in the Earth system. Fires over southern Mexico and Central America (SMCA) are a good example of this, greatly impacting local air quality and regional climate. Here we report that the spring peak (April–May) of fire activities in this region has a distinct quasi-biennial signal based on multiple satellite datasets measuring different fire characteristics. The variability is initially driven by quasi-biennial variations in precipitation. Composite analysis indicates that strong fire years correspond to suppressed ascending motion and weakened precipitation over the SMCA. The anomalous precipitation over the SMCA is further found to be mostly related to the East Pacific–North Pacific (EP-NP) pattern 2 months prior to the fire season. The positive phase of the EP-NP leads to enhanced precipitation over the eastern US but suppressed precipitation over the SMCA, similar to the spatial pattern of precipitation differences between strong and weak fire years. Meanwhile, the quasi-biennial signals in precipitation and fires appear to be amplified by their interactions through a positive feedback loop at short timescales. Model simulations show that in strong fire years, more aerosol particles are released and transported downstream over the Gulf of Mexico and the eastern US, where suspended light-absorbing aerosols warm the atmosphere and cause the ascending motion of the air aloft. Subsequently, a compensating downward motion is formed over the region of the fire source and ultimately suppresses precipitation and intensifies fires. Statistical analysis shows the different durations of the two-way interaction, where the fire suppression effect of precipitation lasts for more than 20 d, while fire leads to a decrease in precipitation at shorter timescales (3–5 d). This study demonstrates the importance of fire–climate interactions in shaping the fire activities on an interannual scale and highlights how precipitation–fire interactions at short timescales contribute to the interannual variability in both fire and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Global-Scale Convergence Obscures Inconsistencies in Soil Carbon Change Predicted by Earth System Models

Soil carbon (C) responses to environmental change represent a major source of uncertainty in the global C cycle. Feedbacks between soil C stocks and climate drivers could impact atmospheric CO 2 levels, further altering the climate. Here, we assessed the reliability of Earth system model (ESM) predictions of soil C change using the Coupled Model Intercomparison Project phases 5 and 6 (CMIP5 and CMIP6). ESMs predicted global soil C gains under the high emission scenario, with soils taking up 43.9 Pg (95% CI: 9.2–78.5 Pg) C on average during the 21st century. The variation in global soil C change declined significantly from CMIP5 (with average of 48.4 Pg [95% CI: 2.0–94.9 Pg] C) to CMIP6 models (with average of 39.3 Pg [95% CI: 23.9–54.7 Pg] C). For some models, a small C increase in all biomes contributed to this convergence. For other models, offsetting responses between cold and warm biomes contributed to convergence. Although soil C predictions appeared to converge in CMIP6, the dominant processes driving soil C change at global or biome scales differed among models and in many cases between earlier and later versions of the same model. Random Forest models, for soil carbon dynamics, accounted for more than 63% variation of the global soil C change predicted by CMIP5 ESMs, but only 36% for CMIP6 models. Although most CMIP6 models apparently agree on increased soil C storage during the 21st century, this consensus obscures substantial model disagreement on the mechanisms underlying soil C response, calling into question the reliability of model predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗