Search NASA⌕ Search

SEARCH · Search NASA

Results for “Columbia Investigation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Evaluating ecosystem water use efficiency and recovery dynamics during flash droughts: insights from observations and model simulations

Flash droughts (FD), rapidly emerging in a warming future, disrupt ecosystems, agriculture, and water security. Ecosystem water use efficiency (WUE), the ratio of gross primary production (GPP) to actual evapotranspiration (AET), balances carbon assimilation and water loss. FD rapidly disrupts this balance, making WUE critical for assessing plant stress and recovery. Here, this study investigates the dynamics of landscape-scale WUE, and the components of GPP and AET under FD utilizing both observed data from the Missouri Ozark AmeriFlux site (US-MOz) and version 2 of the U.S. Department of Energy’s Earth, Energy, Exascale System Model (E3SM) Land Model (ELMv2). Observations and simulations reveal GPP as dominant for WUE during earlier FD events (2005, 2007, 2012), shifting to AET in recent events (2014, 2018). This agreement indicates that the ELM can capture the shifting dynamics of GPP and AET in regulating WUE under FD conditions. However, the ELM systematically underestimates both GPP and AET and does so in a manner that does not preserve their ratio. As a result, WUE is also underestimated, suggesting that GPP is more strongly underestimated than AET. Furthermore, the ELM also underestimates the speed of GPP recovery, producing an artificially prolonged GPP recovery time following FD events. Observed environmental drivers such as vapor pressure deficit (VPD), soil moisture (SM), and predawn leaf water potential (PLWP) effectively predict WUE, but ELM primarily highlights SM, underestimating VPD’s role. This study demonstrates that relying solely on soil moisture fails to capture the rapid hydraulic recovery observed in PLWP, underscoring the necessity of integrating plant hydraulics into land surface models to improve flash drought predictability.

Evapotranspiration↗
DOE OSTI2025

Evaluation of common trap types for capturing wild pigs

Abstract Wild pigs ( Sus scrofa ) are one of the most ecologically harmful invasive vertebrates globally. Due to their negative impacts, wild pig populations are managed via trapping, yet studies investigating various trap designs have been limited in duration, location, and number of trappers. Using data from professional trappers across the southeastern U.S., we evaluated the performance of 3 wild pig trap types (corral, drop, passive net), and the effect of forage season (barren, green‐up, masting) on trapping success. Specifically, we compared trap effectiveness (estimated proportion of targeted wild pigs caught per capture event), trap night effort (number of days to first catch), and the amount of bait and time used per individual take over the lifetime of each trap location across the 3 trap types and forage seasons. Across >850 capture events by 31 trappers, all 3 trap types evaluated were effective at capturing entire social groups of wild pigs. Few differences were observed among trap types, although drop traps in masting season captured a greater estimated proportion of wild pigs per capture event compared to other trap type‐season combinations. Our data also revealed that passive net traps required slightly more bait/wild pig captured than drop traps but there was no difference in the time/wild pig captured among traps. Our findings demonstrate that when operated by wildlife professionals, multiple designs of wild pig traps can be effective at removing targeted groups of wild pigs. Further, our results offer insights regarding the performance of wild pig trap types across seasons that can be used to inform wild pig trapping programs.

Taylor, Charles R. [Warnell School of Forestry and↗
DOE OSTI2025

Advanced spectroscopic studies of (PPh 4 ) 2 [Co(N 3 ) 4 ], a field-induced single-ion magnet

The high-spin Co II complex (PPh 4 ) 2 [Co(N 3 ) 4 ] (Co-N 3 )has been investigated using advanced spectroscopic techniques [far-IR magneto-spectroscopy (FIRMS), high-frequency and high-field EPR (HFEPR), and inelastic neutron scattering (INS)] to study its zero-field-splitting (ZFS), giving spin-Hamiltonian (SH) parameters. The analysis of multi-frequency HFEPR reveals the easy-axis anisotropy with a D value of −10.39(5) cm −1 and a rhombic ratio (E/D) of 0.21(1). The magnetic properties have also been probed by direct-current (DC) magnetometry, suggesting minor differences in anisotropy from the previously reported polymorph (Co-N 3 ′). Ligand-field theory (LFT) analysis indicates that the structures of Co-N 3 and Co-N 3 ′ are closer to D 2d symmetry than other symmetries considered. Alternate-current (AC) susceptibility reveals slow magnetic relaxation under an applied field, indicating that Co-N 3 is a field-induced single-ion magnet (SIM). Here, while both Co-N 3 and Co-N 3 ′ were studied by DC magnetometry, one unusual aspect of the current work on Co-N 3 is that advanced spectroscopies HFEPR, FIRMS, and INS were used to directly observe transitions between ZFS split states, giving accurate SH parameters.

Hand, Adam T. [Univ. of Tennessee, Knoxville, TN (↗
DOE OSTI2025

Tuning the Properties of Lignin‐Derived Deep Eutectic Solvents for Biomass Processing

Lignin-derived deep eutectic solvents (DESs) have been investigated as sustainable green media for biomass processing. However, the properties and processability of DESs have not been fully understood with the chemical structures of their constituents for biomass fractionation. Here, in this article, the properties of the phenolic DESs are discussed with different numbers of functional groups, such as –OCH 3 and –CHO in their hydrogen bond donor (HBD) structures. The formation of DES is significantly related to the hydrogen bond between its constituents, identified by nuclear magnetic resonance (NMR) analysis and density functional theory calculation (DFT). Lower viscosity and net basicity of DES are achieved with fewer –OCH 3 on HBD structures, resulting in enhanced processability and fractionation efficiency. The thermal stability of the DES is also influenced by the –OCH 3 and –CHO of HBD, as indicated by its onset temperature. The recyclability of the phenolic DES is confirmed by the fractionation performance of the recycled DES. Understanding the structural impacts of DES constituents on the properties and performance is crucial for designing solvents in biorefinery applications.

Ryu, Jiae [State University of New York (SUNY), Sy↗
DOE OSTI2025

Resilient stellarator divertor characteristics in the Helically Symmetric eXperiment

Resilient divertor features connected to open chaotic edge structures in the Helically Symmetric eXperiment are investigated. For the first time, an expanded vessel wall was considered that would give space for implementation of a physical divertor target structure. The analysis was done for four different magnetic configurations with very different chaotic plasma edges. A resilient plasma wall interaction pattern was identified across all configurations. This manifests as qualitatively very similar footprint behavior across the different plasma equilibria. Overall, the resilient field lines of interest with high connection length L C lie within a helical band along the wall for all configurations. This resiliency can be used to identify the best location of a divertor. The details of the magnetic footprint’s resilient helical band is subject to specific field line structures which are linked to the penetration depth of field lines into the plasma and directly influence the heat and particle flux patterns. The differences arising from these details are characterized by introducing a new metric, the minimum radial connection min $(\delta_N)$ of a field line from the last closed flux surface. The relationship, namely the deviation from a scaling law, between min $(\delta_N)$ and L C of the field lines in the plasma edge field line behavior suggests that the field lines are associated with structures such as resonant islands, cantori, and turnstiles. This helps determine the relevant magnetic flux channels based on the radial location of these chaotic edge structures and the divertor target footprint. These details will need to be taken into account for resilient divertor design.

divertor↗
DOE OSTI2025

Constraining Erosion Rates and Landscape Evolution With In Situ 10 Be and 26 Al Cosmogenic Nuclides at Table Mountain, Antarctica

Abstract This study investigates surface weathering and sediment preservation at Table Mountain, a high‐elevation, hyperarid, polar landscape in the Transantarctic Mountains. We report cosmogenic nuclide concentrations ( 10 Be and 26 Al) in quartz from bedrock surfaces, erratic boulder lag, and cobbles embedded within Sirius Group sediments to quantify erosion rates. In situ 10 Be and 26 Al depth profiles from a 2.95 m permafrost core in the Sirius Group further constrain surface erosion rates and elucidate landscape stability. Measured 10 Be and 26 Al concentrations from two sandstone bedrock surfaces adjacent to Sirius Group sediments give erosion rates of 0.18–0.28 m/Myr. An erratic sandstone boulder within the lag above the Sirius Group yields erosion rates of ∼0.42 ± 0.03 m/Myr, whereas two cobbles embedded within the Sirius Group yield higher rates of 0.81–1.12 m/Myr. Depth profiles of in situ 10 Be and 26 Al indicate no vertical mixing of Sirius Group permafrost since deposition. Depth profile models are best explained by erosion rates of 0.53 +0.13 / −0.12 m/Myr, and an exposure age of 0.78 +0.06 / −0.08 Ma. We view the model “age” to represent the ∼0.8‐million‐year time‐scale for surface lowering equivalent to one attenuation length of cosmic ray production to achieve steady‐state conditions. Continual exhumation of embedded clasts from within the Sirius Group results in an accumulation of clasts forming the observed erosional lag deposit covering the landscape. Our erosion rates of the Sirius Group surface based on in situ 10 Be and 26 Al depth profiles are an order‐of‐magnitude larger than those based on meteoric 10 Be infiltration and further clarification is required.

58 GEOSCIENCES↗
DOE OSTI2026

Data for An Orphan Gene BOOSTER Enhances Photosynthetic Efficiency and Plant Productivity

Seeds of Col-0 wild type, sig6 T-DNA mutants (CS877785, ABRC), PRL-1-OE, and sig6 T-DNA mutants transfected with PRL-1 (sig6::PRL-1) were planted in 1/2 MS media. Seedlings growth including chlorophyll development defects were investigated across the genotypes. Four-days-old-post-light exposure seedlings were harvested and performed RNAseq analysis with four biological replicates.

Biomass Analytics↗
DOE OSTI2025

Multiplicity dependent 𝐽/𝜓 and 𝜓⁡(2⁢𝑆) production at forward and backward rapidity in 𝑝 + 𝑝 collisions at $\sqrt{𝑠}$ = 200 GeV

Recent measurements of 𝐽/𝜓 production as a function of event charged-particle multiplicity at the collision energies of both the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) show enhanced 𝐽/𝜓 production yields with increasing multiplicity. One potential explanation for this type of dependence is multiparton interactions (MPI). We present the first study of potential autocorrelations at RHIC energies and forward and backward rapidity of self-normalized 𝐽/𝜓 yields and 𝜓⁡(2⁢𝑆) to 𝐽/𝜓 ratio, as a function of self-normalized multiplicity in 𝑝 + 𝑝 collisions. In addition, detailed pythia studies tuned to RHIC energies were performed to investigate the MPI impacts. We find that the PHENIX data at RHIC are consistent with recent LHC measurements and can only be described by pythia calculations that include MPI effects. The forward and backward 𝜓⁡(2⁢𝑆) to 𝐽/𝜓 ratio is found to be less dependent on the charged-particle multiplicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Ion Manipulation from Liquid Xe to Vacuum: Ba-Tagging for a nEXO Upgrade and Future 0 νββ Experiments

Neutrinoless double beta decay (0𝜈𝛽𝛽) provides a way to probe physics beyond the Standard Model of particle physics. The upcoming nEXO experiment will search for 0𝜈𝛽𝛽 decay in 136 Xe with a projected half-life sensitivity exceeding 10 28 years at the 90% confidence level using a liquid xenon (LXe) Time Projection Chamber (TPC) filled with 5 tonnes of Xe enriched to ∼90% in the 𝛽𝛽-decaying isotope 136 Xe. In parallel, a potential future upgrade to nEXO is being investigated with the aim to further suppress radioactive backgrounds and to confirm 𝛽𝛽-decay events. This technique, known as Ba-tagging, comprises extracting and identifying the 𝛽𝛽-decay daughter 136 Ba ion. One tagging approach being pursued involves extracting a small volume of LXe in the vicinity of a potential 𝛽𝛽-decay using a capillary tube and facilitating a liquid-to-gas phase transition by heating the capillary exit. The Ba ion is then separated from the accompanying Xe gas using a radio-frequency (RF) carpet and RF funnel, conclusively identifying the ion as 136 Ba via laser-fluorescence spectroscopy and mass spectrometry. Simultaneously, an accelerator-driven Ba ion source is being developed to validate and optimize this technique. The motivation for the project, the development of the different aspects, along with the current status and results, are discussed here.

a-tagging↗
DOE OSTI2025

Search for an Anomalous Production of Charged-Current 𝜈 𝑒 Interactions without Visible Pions across Multiple Kinematic Observables in MicroBooNE

This Letter presents an investigation of low-energy electron-neutrino interactions in the Fermilab Booster Neutrino Beam by the MicroBooNE experiment, motivated by the excess of electron-neutrino-like events observed by the MiniBooNE experiment. This is the first measurement to use data from all five years of operation of the MicroBooNE experiment, corresponding to an exposure of 1.11 × 10 21 protons on target, a 70% increase on past results. Two samples of electron neutrino interactions without visible pions are used, one with visible protons and one without any visible protons. The MicroBooNE data show reasonable agreement with the nominal prediction, with 𝑝 values ≥26.7% when the two 𝜈 𝑒 samples are combined, though the prediction exceeds the data in limited regions of phase space. The data are further compared to two empirical models that modify the predicted rate of electron-neutrino interactions in different variables in the simulation to match the unfolded MiniBooNE low energy excess. In the first model, this unfolding is performed as a function of electron neutrino energy, while the second model aims to match the observed shower energy and angle distributions of the MiniBooNE excess. This measurement excludes an electronlike interpretation of the MiniBooNE excess based on these models at >99% CLs in all kinematic variables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Thermal stress in degraded forests in the Brazilian Amazon Arc of Deforestation

Understanding thermal stress in tropical forests has taken on new urgency in light of accelerating climate change and expansion of deforestation and forest degradation. Degraded tropical forests in particular may be approaching critical temperature thresholds even more rapidly than intact forests, with implications for tree survival and ecosystem recovery. We investigate thermal stress in degraded tropical forests within the Brazilian Amazon Arc of Deforestation. Using land surface temperature data from the ECOsystem Spaceborne Thermal Radiometer Experiment on the international Space Station (ECOSTRESS), we compared canopy temperatures of intact, selectively logged, and burned forests in Feliz Natal, Mato Grosso, Brazil. Upper canopy temperatures in previously burned forests were 4.1% higher (mean = 36.5 °C) and 50.9% more variable compared to intact and logged forests, which showed remarkably similar temperature distributions (means of 34.9 °C and 35.1 °C, respectively). Modeled leaf temperature distributions based on canopy temperature measurements from one of the warmest days in a 2 year dry season record indicated that 87% of leaves in the warmest burned forest patches exceeded the temperature threshold where respiration surpasses photosynthesis, compared to approximately 72%–74% in intact and logged forests, respectively, in the same time period. After controlling for environmental factors, burned forests were predicted to be 2.6 °C warmer on average [1.39 °C–3.96 °C, 95% credible interval] than intact forests across a 5–40 m canopy height range. Burned forests showed modest thermal recovery over time, with temperatures decreasing by approximately 1.2 °C over a 30 year recovery period. In contrast, logged forests showed minimal thermal differences from intact forests (−0.353 on average, [−1.125–0.274, 95% credible interval]) and negligible change (0.15 °C) over the same timeframe. While the absolute probabilities of exceeding damaging thermal thresholds remain low across all forest types under current climate conditions, the probability of leaves exceeding temperatures that cause permanent leaf damage was ten times higher in burned forests, with implications for the future of burned forest regeneration in water-limited regions of the Amazon basin. In particular, the combination of higher mean temperatures, greater temperature variability, and more frequent exposure to damaging thermal thresholds implies that burned tropical forests will experience substantially higher mortality rates and slower biomass recovery compared to intact and selectively logged forests, especially in water-limited regions where trees cannot rely on evaporative cooling to moderate canopy temperatures.

Brazilian Amazon Arc of Deforestation↗
DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

First Search for Dark Sector 𝑒 + ⁢𝑒 − Explanations of the MiniBooNE Anomaly at MicroBooNE

We present MicroBooNE’s first search for dark sector 𝑒 + ⁢𝑒 − explanations of the long-standing MiniBooNE anomaly. The MiniBooNE anomaly has garnered significant attention over the past 20 years including previous MicroBooNE investigations into both anomalous electron and photon excesses, but its origin still remains unclear. In this Letter, we provide the first direct test of dark sector models in which dark neutrinos, produced through neutrino-induced scattering, decay into missing energy and visible 𝑒 + ⁢𝑒 − pairs comprising the MiniBooNE anomaly. Many such models have recently gained traction as a viable solution to the anomaly while evading past bounds. Using an exposure of 6.87 × 10 20 protons-on-target in the Booster Neutrino Beam, we implement a selection targeting forward-going, coherently produced 𝑒 + ⁢𝑒 − events. After unblinding, we observe 95 events, which we compare with the constrained background-only prediction of 69.7 ±17.3. This analysis sets the world’s first direct limits on these dark sector models and, at the 95% confidence level, excludes the entirety of the single dark neutrino and majority of the dual dark neutrino, parameter space that is viable as a solution to the MiniBooNE anomaly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Impact of Urbanization on Convection, Lightning, and Precipitation Over the Houston Metropolitan Area: Case Study Simulation From the TRACER Campaign

This study investigates the effects of urbanization, specifically land use change and anthropogenic emissions (AE), on convection, lightning, and surface precipitation for a case of summertime sea‐breeze convection observed over the Houston metropolitan area. The unique capabilities of the NASA‐Unified Weather Research and Forecasting model allows us to conduct a series of sensitivity experiments with complex configurations, in particular including multi‐year land model spin‐up simulations, treatment of aerosols and their precursors, and explicit cloud charging and lightning. The simulation results show that urban land use primarily alters the temporal evolution of convection, lightning, and surface precipitation, leading to late afternoon thunderstorm development. The decrease in latent heat flux from the land surface caused by urbanization weakens convection in the early afternoon, while a condition suitable for convection development is maintained in the late afternoon due to less stabilization of the lower troposphere by the weaker convection development and high sensible heat flux from the surface. On the other hand, anthropogenic aerosols directly enhance convection, lightning, and surface precipitation by increasing convective updrafts due to the aerosol‐induced convective invigoration. The combined effects of urban land use and AE lead to even stronger thunderstorms in the late afternoon, mostly consistent with observations. These results indicate that urbanization increases the probability of late afternoon thunderstorms over the Houston area during the summer season. Advanced weather forecasting models that incorporate these urbanization effects might support sustainable urban planning to better mitigate the impacts of urbanization on local weather and public safety.

Iguchi, T. [Univ. of Maryland, College Park, MD (U↗
DOE OSTI2025

Physical properties, internal structure, and the three‐dimensional petrography of CI chondrites

physical properties and the nature of their breccation, we investigated nine samples of the Ivuna and Orgueil CI chondrites ranging in size from 1 mm to 4 cm in approximate diameter. The combined mass of unique material investigated in this work is 113 g. For our investigations, we use ideal gas pycnometry, 3-D laser scanning, x-ray computed microtomography (μCT), and accompanying digital data extraction techniques. We found that the bulk density of the samples ranged from 1.61 to 2.10 g cm −3 . Larger samples tend to have a lower bulk density. Grain density (ranging from 2.44 to 2.55 g cm −3 ) is significantly less variable than the bulk density in our samples and the quantity of porosity (ranging from 14.6% to 33.8%) is the dominant factor in determining the bulk density of CI chondrite material. Our μCT results show that the visible porosity across all sizes of our CI chondrite samples is in the form of cracks, but these cracks can account for less than two-thirds of the porosity in the CI chondrites. Other porosity is not visible, even at μCT resolutions of 2.7 μm voxel edge −1 and we conclude that it is sub-micron in nature. It is not clear if the cracks seen in our samples are indigenous to the chondrites or are a result of terrestrial processes. We also find that the CI chondrites are excellent examples of the fractal-like nature of brecciation, where clasts can be observed at all scales we imaged. The breccias are composed of sub-equant-shaped and sub-rounded-textured clasts like melt-free impact breccias on other solar system bodies. From our μCT volume and digital data extraction, we determine that the Ivuna CI chondrite breccia is organized: the mostly sub-equant clasts within our ~2 cm chunk of Ivuna have a mean diameter of 1.33 mm and their aligned longest axes define a lineation structure. We speculate that the lineation was imparted after fragmentation of the clasts by slight shear on the parent asteroid which could be the result of seismic-related granular flow or mild non-axial impact-related compaction. These data will help to place returned asteroidal material from asteroids 162173 Ryugu and 101955 Bennu and the CI chondrites into a mutual geological context.

CI chondrite↗
DOE OSTI2025

Coulomb excitation of 80 Sr and the limits of the 𝑁 ≈ 𝑍 ≈ 40 island of deformation

The region of 𝑁 ≈ 𝑍 ≈ 40 has long been associated with strongly deformed nuclear configurations. The presence of this strong deformation was recently confirmed through lifetime measurements in 𝑁 ≈ 𝑍 Sr and Zr nuclei. Theoretically, however, these nuclei present a challenge due to the vast valence space required to incorporate all necessary orbitals. Recent state-of-the-art predictions indicate a near axial prolate deformation for 𝑁 = 𝑍 and 𝑁 = 𝑍+2 nuclei between 𝑁 = 𝑍 = 36 and 𝑁 = 𝑍 = 40. In this work we investigate the shores of this island of deformation through a sub-barrier Coulomb excitation study of the 𝑁 = 𝑍+4 nucleus, 80 Sr . Extracting a spectroscopic quadrupole moment of 𝑄 𝑠 ⁡($2^+_1$)= $0.5^{+0.8}_{−0.9}$𝑒⁢b, we find that 80 Sr is inconsistent with significant axial prolate deformation with a significance of 1.5⁢𝜎. This result, albeit with a large uncertainty, indicates that the predicted region of strong prolate deformation around 𝑁 = 𝑍 = 40 is tightly constrained to the quartet of nuclei: 76,78 Sr and 78,80 Zr .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Comparing multi-model ensemble simulations with observations and decadal projections of upper atmospheric variations following the Hunga eruption

The Hunga Tonga–Hunga Ha'apai Model–Observation Comparison (HTHH–MOC) project aims to comprehensively investigate the evolution of volcanic water vapor and sulfur emissions and their subsequent atmospheric impacts and underlying response mechanisms using state-of-the-art global climate models. This study evaluates multi-model ensemble simulations participating in the HTHH–MOC free-run experiment with climate projections for 10 years (2022–2032). Model results are evaluated against satellite observations to assess their ability to reproduce the observed evolution of stratospheric water vapor, aerosols, temperature, and ozone from 2022 to 2024. The participating models accurately capture the observed distribution patterns and associated upper atmospheric responses, providing confidence for their future projections. Model simulations suggest that the Hunga eruption-induced stratospheric water vapor anomaly lasts 4–7 years, with a water vapor e-folding time of 31–43 months. This prolonged water vapor perturbation leads to significant stratospheric and mesospheric cooling, resulting in significant ozone loss in the upper stratosphere and lower mesosphere for 7–10 years. Comparisons between simulations with both SO 2 and H 2 O emissions and those with H 2 O-only emissions indicate that the pronounced dipole response with upper-stratospheric cooling and lower-stratospheric warming is driven by the combined effects of SO 2 and H 2 O injections. These results highlight the prolonged atmospheric impacts of the Hunga eruption and the potential critical role of stratospheric water vapor in modulating long-term atmospheric chemistry and dynamics.

Zhuo, Zhihong [Univ. of Quebec, Montreal, QC (Cana↗