Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flux analysis”

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

Environmental controls on the light use efficiency of terrestrial gross primary production

Abstract Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf‐level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy‐level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C 3 species; however, additional assumptions are required to “scale up” from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed‐effects model. The same model design was fitted to site‐based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple “goodness of fit” comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

54 ENVIRONMENTAL SCIENCES↗

Metabolic engineering to improve production of 3-hydroxypropionic acid from corn-stover hydrolysate in Aspergillus species

Fuels and chemicals derived from non-fossil sources are needed to lessen human impacts on the environment while providing a healthy and growing economy. 3-hydroxypropionic acid (3-HP) is an important chemical building block that can be used for many products. Biosynthesis of 3-HP is possible; however, low production is typically observed in those natural systems. Biosynthetic pathways have been designed to produce 3-HP from a variety of feedstocks in different microorganisms. In this study, the 3-HP β-alanine pathway consisting of aspartate decarboxylase, β-alanine-pyruvate aminotransferase, and 3-hydroxypropionate dehydrogenase from selected microorganisms were codon optimized for Aspergillus species and placed under the control of constitutive promoters. The pathway was introduced into Aspergillus pseudoterreus and subsequently into Aspergillus niger, and 3-HP production was assessed in both hosts. A. niger produced higher initial 3-HP yields and fewer co-product contaminants and was selected as a suitable host for further engineering. Proteomic and metabolomic analysis of both Aspergillus species during 3-HP production identified genetic targets for improvement of flux toward 3-HP including pyruvate carboxylase, aspartate aminotransferase, malonate semialdehyde dehydrogenase, succinate semialdehyde dehydrogenase, oxaloacetate hydrolase, and a 3-HP transporter. Overexpression of pyruvate carboxylase improved yield in shake-flasks from 0.09 to 0.12 C-mol 3-HP C-mol -1 glucose in the base strain expressing 12 copies of the β-alanine pathway. Deletion or overexpression of individual target genes in the pyruvate carboxylase overexpression strain improved yield to 0.22 C-mol 3-HP C-mol -1 glucose after deletion of the major malonate semialdehyde dehydrogenase. Further incorporation of additional β-alanine pathway genes and optimization of culture conditions (sugars, temperature, nitrogen, phosphate, trace elements) for 3-HP production from deacetylated and mechanically refined corn stover hydrolysate improved yield to 0.48 C-mol 3-HP C-mol -1 sugars and resulted in a final titer of 36.0 g/L 3-HP. The results of this study establish A. niger as a host for 3-HP production from a lignocellulosic feedstock in acidic conditions and demonstrates that 3-HP titer and yield can be improved by a broad metabolic engineering strategy involving identification and modification of genes participated in the synthesis of 3-HP and its precursors, degradation of intermediates, and transport of 3-HP across the plasma membrane.

09 BIOMASS FUELS↗

Software Verification of VARPOW

The VARPOW program is a post-processing utility program for DIF3D, specifically DIF3DVARIANT, and it was developed to provide interface files for the thermal analysis program DASSH. The basic methodology of VARPOW is to take the neutron and gamma flux (moments) calculated by DIF3D (or GAMSOR) and combine them with the heating (coefficient) cross sections to calculate the spatial power distributions using the DIF3DVARIANT spatial basis. VARPOW can use the output from GAMSOR (both steady state neutron and gamma flux calculations) or standard DIF3D/REBUS calculations (neutron flux only). The correct approach for defining the power distribution is to use GAMSOR as its purpose was to properly compute the gamma heating throughout the modeled domain. The power densities calculated by VARPOW are broken into fuel, cladding and coolant terms for which isotope-wise categorization is needed. VARPOW has built in options the user can select for the isotope categorization or VARPOW can import a file that details the isotope categorization. VARPOW can export the solution in the polynomial basis of DIF3D-VARIANT or the monomial basis of DIF3D-VARIANT. The purpose of this work is to verify the power distribution results calculated by VARPOW from both the GAMSOR and DIF3D input options and verify that the input and output options are consistent with the manual. Hand calculation and independent numerical calculation are used for this verification work.

97 MATHEMATICS AND COMPUTING↗

Software Verification of VARPOW

The VARPOW program is a post-processing utility program for DIF3D, specifically DIF3DVARIANT, and it was developed to provide interface files for the thermal analysis program DASSH. The basic methodology of VARPOW is to take the neutron and gamma flux (moments) calculated by DIF3D (or GAMSOR) and combine them with the heating (coefficient) cross sections to calculate the spatial power distributions using the DIF3DVARIANT spatial basis. VARPOW can use the output from GAMSOR (both steady state neutron and gamma flux calculations) or standard DIF3D/REBUS calculations (neutron flux only). The correct approach for defining the power distribution is to use GAMSOR as its purpose was to properly compute the gamma heating throughout the modeled domain. The power desnities calculated by VARPOW are broken into fuel, cladding and coolant terms for which isotope-wise categorization is needed. VARPOW has built in options the user can select for the isotope categorization or VARPOW can import a file that details the isotope categorization. VARPOW can export the solution in the polynomial basis of DIF3D-VARIANT or the monomial basis of DIF3D-VARIANT. The purpose of this work is to verify the power distribution results calculated by VARPOW from both the GAMSOR and DIF3D input options and verify that the input and output options are consistent with the manual. Hand calculation and independent numerical calculation are used for this verification work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Design of Critical Experiments Involving Graphite elements at the IPEN/MB-01 Plate-Type Core

A series of critical experiments were designed involving the use of characterized graphite elements in the new plate-type core of the IPEN/MB-01 reactor. Eleven configurations with different graphite quantities are considered, each producing distinct effects on neutron fluxes and offering opportunities to analyze varying sensitivities to k eff in the system. The analysis shows that the experiments can be performed with acceptably low k eff uncertainties considering an adequate characterization of the graphite elements to be used. The sensitivity results obtained using SCALE TSUNAMI calculations show that the experiments have a notable k eff sensitivity coefficient to graphite. In the configurations tested, the k eff sensitivity to graphite TSL is low. This experiment will be useful for criticality safety validation of applications using uranium fuel enriched to around 19.75 wt % 235 U, light water and graphite. It will also be useful to gain more insight into graphite material properties and their effect on k eff as a result of using the well-characterized graphite elements. The execution of the experiments is planned for the summer of 2025. The 11 designed configurations serve as a basis for the final design, and fewer configurations will be executed. Once executed, these critical experiments will be evaluated and submitted for publication in the International Criticality Safety Benchmark Evaluation Project Handbook, supporting the DOE/NRC Collaboration for Criticality Safety Support for Commercial-Scale HALEU for Fuel Cycles and Transportation project goal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Software Verification of VARPOW

The VARPOW program is a post-processing utility program for DIF3D, specifically DIF3DVARIANT, and it was developed to provide interface files for the thermal analysis program DASSH. The basic methodology of VARPOW is to take the neutron and gamma flux (moments) calculated by DIF3D (or GAMSOR) and combine them with the heating (coefficient) cross sections to calculate the spatial power distributions using the DIF3DVARIANT spatial basis. VARPOW can use the output from GAMSOR (both steady state neutron and gamma flux calculations) or standard DIF3D/REBUS calculations (neutron flux only). The correct approach for defining the power distribution is to use GAMSOR as its purpose was to properly compute the gamma heating throughout the modeled domain. The power desnities calculated by VARPOW are broken into fuel, cladding and coolant terms for which isotope-wise categorization is needed. VARPOW has built in options the user can select for the isotope categorization or VARPOW can import a file that details the isotope categorization. VARPOW can export the solution in the polynomial basis of DIF3D-VARIANT or the monomial basis of DIF3D-VARIANT. Recently the clad neutron damage calculation capability was introduced into VARPOW, and it now allows a pathway to get distributions of clad DPA. The purpose of this work is to verify the power distribution, fast neutron flux and clad DPA results calculated by VARPOW from both the GAMSOR and DIF3D input options and verify that the input and output options are consistent with the manual. Hand calculation and independent numerical calculation are used for this verification work.

Zhong, Zhaopeng [Argonne National Laboratory (ANL)↗

Carbon dioxide (CO2) flux, amplicon sequencing and liquid chromatography–mass spectrometry (LC-MS) data for pyrogenic organic matter (PYOM) amended soil incubations in lab, Madison, WI, 2021-22

This dataset comprises soil-pyrogenic organic matter (PyOM) incubation experiments conducted in a laboratory setting. The focus was on measuring carbon dioxide (CO2) flux emitted in the headspace of incubation jars over a one-month period. The aim was to estimate carbon mineralization from two different carbon fractions in PyOM, as well as carbon mineralization from soil organic carbon (SOC). Bacterial and fungal community profiles were tracked at various time points during the incubation to observe changes in the overall community and specific responders to PyOM. LC-MS analysis of soil samples was also performed to monitor changes in the soil's chemical composition.The provided dataset includes processed CO2 flux measurements obtained from the Picarro-multiplexer measurement setup for each sample incubation jar. This data is available in the "CO2Flux.csv" file. Additionally, the dataset includes CO2 flux data partitioned to estimate mean carbon mineralization from the two PyOM carbon fractions, bulk PyOM, and SOC. The cumulative carbon mineralization data is provided in "cml_respired.csv," while the rate of carbon mineralization data is in "rate_respired.csv." The file "cml_respired_lastcyc.csv" contains cumulative carbon mineralization data for the PyOM carbon fractions, bulk PyOM, and SOC throughout the entire incubation period.For the LC-MS analysis, the raw data can be found in the "20220721_NZ_PeakHeights_FinalAnalysis.csv" file, while sample metadata is provided in "20220721_NZpilot_Metadata.csv". Processed LC-MS data is available in the file "20220721_NZ_FINAL_msdat_sub3.csv". Lastly, the dataset includes the relative abundance data for bacterial genera that exhibited a significant positive response to the addition of PyOM produced at 350 degrees Celsius, which can be found in "resp_genera_350.csv".

54 ENVIRONMENTAL SCIENCES↗

NSTX-U theory, modeling and analysis results

Here, the mission of the low aspect ratio spherical tokamak NSTX-U is to advance the physics basis and technical solutions required for optimizing the configuration of next-step steady-state tokamak fusion devices. NSTX-U will ultimately operate at up to 2 MA of plasma current and 1 T toroidal field on axis for 5 s, and has available up to 15 MW of neutral beam injection power at different tangency radii and 6 MW of high harmonic fast wave heating. With these capabilities NSTX-U will develop the physics understanding and control tools to ramp-up and sustain high performance fully non-inductive plasmas with large bootstrap fraction and enhanced confinement enabled via the low aspect ratio, high beta configuration. With its unique capabilities, NSTX-U research also supports ITER and other critical fusion development needs. Super-Alfvénic ions in beam-heated NSTX-U plasmas access energetic particle (EP) parameter space that is relevant for both α-heated conventional and low aspect ratio burning plasmas. NSTX-U can also generate very large target heat fluxes to test conventional and innovative plasma exhaust and plasma facing component solutions. This paper summarizes recent analysis, theory and modelling progress to advance the tokamak physics basis in the areas of macrostability and 3D fields, EP stability and fast ion transport, thermal transport and pedestal structure, boundary and plasma material interaction, RF heating, scenario optimization and real-time control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global mapping of urban thermal anisotropy reveals substantial potential biases for remotely sensed urban climates

Urban thermal anisotropy (UTA) drastically biases satellite-derived urban surface temperatures and fluxes, and consequently inhibits a more comprehensive understanding of global urban climates. However, UTA patterns and associated biases in observed urban climate variables have not been investigated across an adequate number of global cities with diverse contexts; nor is it known whether there are globally measurable factors that are closely related to the UTA intensity (UTAI, quantified as the maximum difference between the nadir and off-nadir urban thermal radiation). Here we investigate the UTAI over more than 5500 cities worldwide using the moderate resolution imaging spectroradiometer (MODIS) multi-angle land surface temperature (LST) observations from 2003 to 2021. The results show that the global mean UTAI can reach 5.1, 2.7, 2.4, and 1.7 K during summer daytime, winter daytime, summer nighttime, and winter nighttime, respectively. Using nadir LST observations as a reference, our analysis reveals that UTA can lead to an underestimation of satellite-based urban surface sensible heat fluxes ( H ) by 45.4% and surface urban heat island intensity ( I s ) by 43.0% when using LST observations obtained from sensor viewing zenith angles (VZAs) of ±60°. However, practitioners can limit the biases of H and I s within ±10% by using LSTs from sensor VZAs within ±30°. We find that UTAI is closely related to urban impervious surface percentage and surface air temperature across global cities. In conclusion, these findings have implications for angular normalization of satellite-retrieved instantaneous LST observations across cities worldwide.

54 ENVIRONMENTAL SCIENCES↗

Steady-State Mechanical Analysis for Target Assembly in the Material Plasma Exposure eXperiment Facility

The Material Plasma Exposure eXperiment (MPEX) steady-state linear plasma facility is currently under design at Oak Ridge National Laboratory to expose target specimens to fusion divertor regimes. The neutron-irradiated target is actively cooled and remote handled in the MPEX facility for conducting plasma-material–interaction (PMI) experiments. Here, the steady-state stresses in the target and target assembly system are investigated using two-dimensional (2-D) and three-dimensional (3-D) models to provide expected stresses/strains under the heat loads to which various system components would be exposed during MPEX operation. The calculated temperatures from the 2-D axisymmetric mechanical model were found to be in excellent agreement with those from the full 3-D thermohydraulic model, providing a strong model validation. Numerical simulation results for the steady-state mechanical model indicate nonuniform distributions for the temperature, stress, and deformation within the critical components. For the initial design, the deformation results indicate possible gap openings between contacting surfaces below the plasma-facing materials. To reduce the possibility of interfacial gap opening, the target assembly was slightly changed and evaluated using the 2-D stress model. Numerical simulation results indicate that the interfacial gap openings can be minimized without drastically changing the entire target assembly. The stress-strain conditions for the target will be further used to assess the appropriate operation during MPEX experiments and gain insight into materials science phenomena during PMI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a Fast Spectrum Test Loop at the Advanced Test Reactor

The United States currently has no infrastructure to provide a fast neutron flux for research purposes. With the increased interest in Sodium Fast Reactor (SFR) technology, this type of infrastructure is necessary for further development of key components such as Versatile Test Reactor. The purpose of this research is to characterize the thermal hydraulics of a test loop capable of providing SFR like operating conditions. This includes showcasing feasibility in mechanical design, reaching a high fast to thermal neutron flux ratio ideally near 15, an inlet temperature of 600 K, a temperature increase of 50 K from inlet to outlet, and test chamber radial temperature profiles similar to SFR fuel rods. Existing infrastructure must be used to provide the neutron flux for the test loop, and it was decided that the Advanced Test Reactor (ATR) was the optimal choice due to its availability and previous work conducted similar to this project. To achieve the goals of this project, several mechanical designs were put forward with the most feasible design being chosen. The neutronics and thermal hydraulics analysis were conducted in tandem, with an iterative process used to maximize efficiency of both systems. Analysis of the neutronics system was done using Attila and MCNP while the thermal hydraulic analysis was done using RELAP5. This document outlines the thermal hydraulic analysis, though it includes the results from the mechanical design and neutronics analysis. The loop design includes a primary system with molten sodium, a set of “booster fuel plates” to increase the fast flux within the test section, a water system to cool the booster fuel plates, a helium system to separate the sodium and water systems, and a thermal neutron shield to decrease the thermal flux. Each of these systems is modeled in the neutronics and thermal hydraulic analyses. The RELAP5 model showcased the feasibility of the design proposed by this report. For all designs tested, the desired 50 K temperature increase from inlet to outlet was reached. The fuel rod radial temperature profile was lower than was expected, likely due to a lower rod thermal power than expected. A parametric study was also conducted to characterize the system with three fuel pins. This study showed the effect of thermal power, fluid pressure, fluid velocity, and inlet temperature on the system. ©Copyright by Anton Higgins December 11, 2019 All Rights Reserved

99 GENERAL AND MISCELLANEOUS↗

The STAT7 Code for Statistical Propagation of Uncertainties In Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

STAT7 was written to automate many of the steady-state thermal hydraulic safety calculations for the MIT research reactor, both for conversion of the reactor from highly enriched uranium fuel to low-enriched uranium fuel and for future fuel re-loads after the conversion. A Monte-Carlo statistical propagation approach is used to treat uncertainties in important parameters in the analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to ONB for each axial node of each stripe of each plate of each element in the core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Effect of Cell Geometry and Trigger Method on the Risks Associated with Thermal Runaway of Lithium-Ion Batteries

Consideration of thermal runaway heat output variability is paramount for the development of safe lithium-ion battery assemblies. This study utilizes data gathered from fractional thermal runaway calorimetry (FTRC) experiments to conduct a comparative analysis of thermal runaway heat output for three cell formats (18650, 21700, and 33600) as a function of trigger method (heaters, internal short-circuiting device, and nail penetration). The analysis is based on comparisons for the calculated total energy yield, fractional energy yield, heat rate, and heat flux. This study reveals that nail penetration tends to result in higher thermal runaway heat output for larger cells (21700 & 33600); these experiments also tended to result in higher fractions of the total energy being released through the cell body. The smaller cells (18650) did not appear to have significant variation in heat output as a function of trigger method. This finding suggests that, for this cell type, worst-case scenario heat output could be achievable in assembly level testing regardless of the utilized trigger method. This study also demonstrates successful translation of FTRC results, as recorded in the Battery Failure Databank, into meaningful analysis that breaks down the influence of specific conditions on thermal runaway heat output.

calorimetry↗

Extension of the high-resolution thermal-hydraulics code ESCOT to hexagonal core geometries for multi-physics calculations

The extension of the capabilities of the pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT to analyze hexagonal fueled cores and its performance are presented. ESCOT is an accurate yet fast core thermal-hydraulics solution aiming at high-fidelity and high-resolution multi-physics core analysis in the framework of massively parallel computing platforms. Its algorithm solution is based on the four-equation drift-flux model for two-phase calculations, these are numerically solved by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models such as turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analysis. ESCOT is parallelized by a double (radial and axial) domain decomposition that enables its highly parallelized execution. The coupling of the code with the neutronics whole core solver for hexagonal geometries nTRACER is described. The newly implemented ESCOT features are validated by comparing single assembly and full core steady state nTRACER-ESCOT solutions with nTRACER standalone internal one-dimensional T/H solver results. The validation problems are based on the VVER 440 and VVER 1000 cores. ESCOT results show differences within an acceptable range with respect to the simple 1D nTRACER built-in solver. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Design and analysis of actively-cooled, edge-transport diagnostic for long-pulsed operation in WEST

Next step fusion devices that will operate in steady-state will require complex plasma-facing components (PFCs) that can survive the harsh environment over long timescales not common in current devices. This will require robust plasma facing surfaces that are integrated with active cooling systems. In a collaboration between CEA and ORNL, a plasma-interacting diagnostic is being designed for the W Environment in Steady-state Tokamak (WEST) in Cadarache, France which requires plasma-facing protection like those needed for steady-state PFCs. This integrated diagnostic studies edge transport and impurity migration within WEST, and is planned to include imbedded temperature Langmuir probe sensors, as well as removeable sample slots for ex-situ surface analysis of plasma-material interactions. The entire assembly is expected to move into the plasma edge for periods up to 1000 s having an energy removeable capability of ∼6 kW and seeing a peak heat flux of more than 7 MW/m2. The assembly compliments WEST high fluence campaigns that plan for multiple 1000 s pulses. Within these specifications, the assembly will require a refractory metal plasma facing surface and integral cooling in order to function within the limited space allotted for such diagnostics. Because of the limited space and linear actuator needs, additive manufacturing of the high heat flux working end of the assembly is being considered which could allow for precision cooling-channels and lighter weight designs. Conceptual design along with simulation and analysis results will be presented for this complex diagnostic with novel PFCs.

Lumsdaine, Arnold↗

Analyzing the upscaling potential and geospatial siting of calcination-free calcium hydroxide production in the United States

This study evaluates the techno-economic feasibility and the embodied carbon dioxide intensity (eCI) of a novel process for producing nominally pure (>95 mass %) calcium hydroxide without the need for the thermal calcination of limestone. The process relies on the aqueous extraction of calcium from alkaline industrial wastes following which portlandite (Ca(OH) 2 : CH, a.k.a. slaked lime or hydrated lime) is precipitated by application of a waste-heat based thermal swing. This approach takes advantage of the temperature dependent solubility of CH at ambient pressure. We evaluated the feasibility of implementing this process in the U.S. based on the geospatial availability of waste heat and slags as a Ca-source. For the base case, the cost of production of “Low-Temperature Portlandite (LTP)” is 2-to-3 times that of traditional portlandite (~$180/tonne). The main driver of cost is the electricity demand for reverse osmosis (RO) which is used to concentrate Ca-ions in solution, and the costs of membrane replacement. Our sensitivity analysis showed that parity with the cost of production of traditional portlandite is readily achievable by selecting membranes with better durability (i.e., better pH resistance) and flux (i.e., higher permeability) without sacrificing selectivity. Significantly, LTP features an eCI that is between 40%- and - 80 % lower than its traditional counterpart when electricity is sourced from natural gas combustion or wind power, respectively. Finally, our geospatial analysis reveals that there are three areas in the U.S. with the potential for implementation of industrial-scale facilities that could produce at least 50 tonnes of pure Ca(OH) 2 per day, while achieving a production cost of ~$270 per tonne of Ca(OH) 2 , owing to the proximity between slag feedstocks and waste heat sources.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Constraints on the IBEX Ribbon’s Origin from Its Evolution over a Solar Cycle

Abstract In 2009, the Interstellar Boundary Explorer (IBEX) discovered a narrow “ribbon” of energetic neutral atom emissions across the sky with properties correlated with the solar wind latitudinal structure and the interstellar magnetic field draped around the heliosphere. It is widely believed that the ribbon is formed from the escape of heliospheric ENAs into the local interstellar medium and their eventual return as secondary ENAs. However, there is no consensus on the rate of pitch angle scattering of these PUIs before they become secondary ENAs. We test two opposing limits of scattering rates (“weak” versus “strong”) by solving a time-dependent model of the ribbon that evolves with the solar cycle, and we compare them to IBEX observations over 2009–2019. First, we find that both models qualitatively reproduce the evolution of IBEX fluxes for most of the data set, with a few exceptions, although the strong (or “spatial retention”) scattering model greatly underestimates the observed fluxes. Regardless, time dependence of fluxes cannot distinguish these models. Second, the ribbon’s geometric properties, i.e., its center and radius, are significantly different between the models. The spatial retention model reproduces the observed ribbon centers as a function of energy and time slightly better than the weak scattering model, and the spatial retention model reproduces the observed ribbon radius over energy and time almost perfectly, whereas the weak scattering model compares poorly. Our analysis favors the spatial retention mechanism as the source of the IBEX ribbon, but it requires modification to increase the flux of ENAs observed at 1 au.

Zirnstein, E. J. (ORCID:0000000172400618)↗

Estimates of Southern Hemispheric Gravity Wave Momentum Fluxes across Observations, Reanalyses, and Kilometer-Scale Numerical Weather Prediction Model

Abstract Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis. Significance Statement Gravity wave–induced momentum forcing forms a key component of the middle atmospheric circulation. However, complete knowledge of gravity waves, their atmospheric effects, and their long-term trends are obscured due to limited global observations, and the inability of current climate models to fully resolve them. This study combines a kilometer-scale forecast model, modern reanalysis, and a coarse-resolution climate model to first compare the resolved and parameterized momentum fluxes by gravity waves generated over the Andes, and then evaluate the fluxes using a state-of-the-art ground-based Rayleigh lidar. Our analysis reveals shortcomings in current model parameterizations of gravity waves in the middle atmosphere and highlights the sensitivity of the estimated flux to the formulation used.

Meteorology & Atmospheric Sciences↗