Search NASA⌕ Search

SEARCH · Search NASA

Results for “disruption prediction”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Improbability of Post-Closure Criticality in Compacted Criticality Control Overpacks after Room Closure at Waste Isolation Pilot Plant

As part of its periodic re-certification of the Waste Isolation Pilot Plant (WIPP), an operating repository in bedded salt for the disposal of transuranic (TRU) waste from atomic energy defense activities, the United States Environmental Protection Agency expects a re-evaluation of features, events, and processes, such as post-closure nuclear criticality. Although salt creep beneficially encapsulates the TRU waste in the closed WIPP repository, the spacing between an array of waste packages is disrupted as the salt creep closes disposal rooms and containers lose structural integrity. For most TRU waste, the possibility of post-closure criticality is exceedingly small either because the salt neutronically isolates TRU waste canisters or because closure of a disposal room from salt creep does not sufficiently compact the low mass of fissile material. The criticality evaluation was updated, however, because of the introduction of criticality control overpack (CCO) containers, which may dispose up to 380 fissile gram equivalent plutonium-239 in each container. The criticality potential is evaluated through high-fidelity geomechanical modeling of a disposal room filled with CCO containers during two representative conditions: (1) large salt block fall, and (2) gradual disposal room closure from salt creep. Geomechanical models of roof fall demonstrate three tiers of CCO containers are not greatly disrupted. Geomechanical models of gradual room closure from salt creep (without brine seepage and subsequent gas generation to permit maximum room closure) were used to predict irregular arrays of closely packed CCOs after 1000 years, when room closure has asymptotically approached maximum compaction. Models of spheres or cylinders with 380 fissile gram equivalent of plutonium (as oxide) at the predicted irregular compacted spacing demonstrate that an array of CCO containers is not critical when surrounded by salt and magnesium oxide, provided the mass of hydrogenous material shipped in CCO containers (usually plastics) is controlled or boron carbide (a neutron poison) is mixed with the fissile contents.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Force-free Wave Interaction in Magnetar Magnetospheres: Computational Modeling in Axisymmetry

Abstract Crustal quakes of highly magnetized neutron stars can disrupt their magnetospheres, triggering energetic phenomena like X-ray and fast radio bursts. Understanding plasma wave dynamics in these extreme environments is vital for predicting energy transport across scales to the radiation length. This study models relativistic plasma wave interaction in magnetar magnetospheres with force-free electrodynamics simulations. For propagation along curved magnetic field lines, we observe the continuous conversion of Alfvén waves to fast magnetosonic (FMS) waves. The conversion efficiency can be up to three times higher when counter-propagating Alfvén waves interact in the equatorial region. Alfvén waves generate FMS waves of twice their frequency during their first crossing of the magnetosphere. After the initial transient burst of FMS waves, Alfvén waves convert to FMS waves periodically, generating variations on timescales of the magnetospheric Alfvén wave crossing time. This decaying FMS wave tail carries a significant portion (half) of the total energy emitted. Plastic damping of “bouncing” Alfvén waves by the magnetar crust has minimal impact on the FMS efficiency. We discuss the implications of the identified wave phenomena for magnetar observations. Outgoing FMS waves can develop electric zones, potential sources of coherent radiation. Long wavelength FMS waves could generate FRBs through reconnection beyond the light cylinder.

Astronomy & Astrophysics↗
NASA NTRS2020

POEMMA’s Target-of-Opportunity Sensitivity to Cosmic Neutrino Transient Sources

We investigate the capability of the Probe Of Extreme Multi-Messenger Astrophysics (POEMMA)in performing Target-of-Opportunity (ToO) neutrino observations. POEMMA is a proposed space-based probe-class mission for ultra high-energy cosmic ray and very-high-energy neutrino detection using two spacecraft, each equipped with a large Schmidt telescope to detect optical and near-ultraviolet signals generated by extensive air showers (EASs). POEMMA will be sensitive to Cherenkov radiation from upward-moving EASs initiated by tau neutrinos interacting in the Earth. POEMMA will be able to quickly re-point (90◦in 500 s) each of the two spacecraft to the direction of an astrophysical source, which in combination with its orbital speed will provide it with unparalleled capability to follow up transient alerts. We calculate POEMMA’s transient sensitivity for two observational configurations for the satellites (ToO-stereo and ToO-dual for smaller and larger satellite separations, respectively) and investigate the impact of variations arising due to POEMMA’s orbital characteristics on its sensitivity to tau neutrinos in various regions of the sky. We explore separate scenarios for long (∼105−6s) and short (∼103s) duration events, accounting for intrusion from the Sun and the Moon in the long-duration scenario. We compare the sensitivity and sky coverage of POEMMA for ToO observations with those for existing experiments (e.g., IceCube, ANTARES, and the Pierre Auger Observatory) and other proposed future experiments (e.g., GRAND200k).For long bursts, we find that POEMMA will provide a factor of&7 improvement in average neutrino sensitivity above 300 PeV with respect to existing experiments, reaching the level of model predictions for neutrino fluences at these energies and above from several types of long-duration astrophysical transients (e.g., binary neutron star mergers and tidal disruption events). For short bursts, POEMMA will improve the sensitivity over existing experiments by at least an order of magnitude forEν∼>100 PeV in the “best-case” scenario. POEMMA’s orbital characteristics and rapid re-pointing capability will provide it access to the full celestial sky, including regions that will not be accessible to ground-based neutrino experiments. Finally, we discuss the prospects for POEMMA to detect neutrinos from candidate astrophysical neutrino sources in the nearby universe. Our results demonstrate that with its improved neutrino sensitivity at ultra-high energies and unique full-sky coverage, POEMMA will be an essential, complementary component in a rapidly expanding multi-messenger network.

Tonia M Venters↗
NASA NTRS2011

Consequence and Resilience Modeling for Chemical Supply Chains

The U.S. chemical sector produces more than 70,000 chemicals that are essential material inputs to critical infrastructure systems, such as the energy, public health, and food and agriculture sectors. Disruptions to the chemical sector can potentially cascade to other dependent sectors, resulting in serious national consequences. To address this concern, the U.S. Department of Homeland Security (DHS) tasked Sandia National Laboratories to develop a predictive consequence modeling and simulation capability for global chemical supply chains. This paper describes that capability , which includes a dynamic supply chain simulation platform called N_ABLE(tm). The paper also presents results from a case study that simulates the consequences of a Gulf Coast hurricane on selected segments of the U.S. chemical sector. The case study identified consequences that include impacted chemical facilities, cascading impacts to other parts of the chemical sector. and estimates of the lengths of chemical shortages and recovery . Overall. these simulation results can DHS prepare for and respond to actual disruptions.

Stamber, Kevin L.↗
DOE OSTI2025

Multi-device analysis of energy loss duration and pellet penetration with implications for shattered pellet injection in ITER

A robust disruption mitigation system (DMS) requires accurate characterization of key disruption timescales, one of the most notable being the thermal quench (TQ). Recent modeling of shattered pellet injection (SPI) into ITER plasmas, using JOREK and INDEX, suggests long TQ durations (6–10 ms) and slow cold front propagation due to the large plasma size. If validated, these predictions would have an impact on the desired pellet parameters and mitigation strategies for the ITER DMS. To resolve these questions, a database of SPI experiments from several small-to-large sized devices (J-TEXT, KSTAR, AUG, DIII-D, and JET) has been compiled under the auspices of the International Tokamak Physics Activity MHD, disruptions, and control topical group. Analysis of the energy loss duration (proxy for the TQ duration) with machine size is presented for both mixed neon/deuterium (Ne/D) SPI and pure deuterium (D) SPI. Several metrics for the energy loss onset (e.g. soft x-ray signal drop, I p dip, and radiation flash) were considered as the conventional metric, electron cyclotron emission, is often cut-off during SPI. Several scalings with different onset metrics showed an increase in energy loss duration with machine size. The energy loss duration was additionally shown to be a function of the ratio between the number of SPI neon atoms injected and the stored energy. Analysis of the pellet shard position relative to the cold front found that in larger devices, pellets are typically found inboard of the q = 2 surface at the energy loss onset. Lastly, the delay between the pellet shards hitting the q = 2 surface and the energy loss onset was additionally found to increase with machine size. This suggests that the pellet shards in large devices will penetrate faster and further than the cooling front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

STOCHASTIC OPTIMAL POWER FLOW FOR REAL-TIME MANAGEMENT OF DISTRIBUTED RENEWABLE GENERATION AND DEMAND RESPONSE (Final Report)

To meet the grand challenge of a sustainable energy future, there has been a surge of interest in renewable energy. Today, the uncertainty associated with renewable resources is handled by using operating reserves. The high penetration of renewable resources, however, introduces difficult-to-control dynamics and challenges for power system operation. Decision support tools are necessary at the bulk system operational level to recognize and efficiently utilize renewable resources and distributed demand response products in concert with traditional grid resources. It is envisaged that responsive load can potentially have very significant cost advantages over either spinning or non-spinning ramping reserve. Critical decisions are made during hour(s)-ahead and real-time power system operation regarding the commitment and dispatch of generators to ensure power delivery is both reliable and economic. These decisions are typically made by a security constrained optimal flow, which determines future generator commitments, dispatches, and ensures adequate reserves are available in the event of a contingency (unexpected outage) or if future system conditions deviate from forecasts. However, security has been always based on a pre-specified subset of contingency constraints whose enforcement does not guarantee security under all possible future possibilities while also giving little or no weight to the likelihood of each contingent event or the severity of its consequences. Existing tools, which are based exclusively on deterministic optimization models, do not yield optimal operational decisions to address these new challenges, in terms of both reliability and cost-effectiveness. This project has focused on developing a stochastic optimal power flow (SOPF) framework, which integrates renewable resource uncertainty, load uncertainty, distributed storage (DS), demand response (DR) products, in a holistic manner to address the uncertainty associated with ever-increasing renewable resources, along with the inclusion of distributed demand response products in future power systems. A proof-of-concept problem was created using the Pennsylvania-Jersey-Maryland (PJM) power system network. Synthetic wind generation was added to the system to simulate 50% wind penetration. A 1-hour test of SOPF operation indicated more than 6% operational cost savings. The project continued by adding the Midwestern Independent System Operator (MISO) as a partner, with focus shifting from SOPF to Stochastic Look-Ahead Unit Commitment (SLAC). Unlike PJM, MISO is faced with significant renewable energy resources within its footprint and is challenged with substantial uncertainty in its operations. The SLAC distinguishes itself from existing tools that operators use. At best, today’s tools solve two to three cases independently, where one or two system parameters, such as forecasted load level (e.g., a low, base, and high forecast), are varied and the resulting scenarios are analyzed independently. The stochastic-based optimization of SLAC leverages statistical information from an ensemble of potential operational scenarios and their respective likelihood. The SLAC output can be translated into valuable information to the operator such as suggested commitments, optimal scheduling and dispatch of resources, reserve requirements at both locational and zonal resolutions, ramping availability and requirements, availability of demand response including operational guidance concerning the near-term and real-time coordination between distributed energy resources, and utilization of distributed storage resources. The developed SOPF/SLAC tool, a stand-alone tool compatible with existing EMSs, will provide system operators with unprecedented visibility, flexibility and predictability to these resources and operational guidance concerning the real-time coordination between DERs and DR/DS products. The game changing and practical impact of this disruptive technology will be dramatic and will usher in a new era in the electric power industry, wherein green energy concepts are fully embraced, and electric power costs are lowered throughout the nation.

42 ENGINEERING↗
NASA NTRS2024

How Many Vesta-Like Bodies Existed in the Asteroid Belt?

Asteroid 4 Vesta is typically thought to be the parent body of the HED (howardite, eucrite, and diogenite) meteorites due to spectral similarities. The discovery of asteroids far from Vesta with HED-like spectra like (1459) Magnya and HED-like meteorites (e.g., NWA 011) with anomalous oxygen isotopic values compared to typical HEDs is evidence that other Vesta-like bodies formed. We broadly define a Vesta-like body as a differentiated object with a crust composed primarily of low-Ca pyroxene and plagioclase feldspar. We estimate the number of Vesta-like bodies that did form by looking at the astronomical evidence; the oxygen isotopic, chemical, and petrologic evidence; and the iron meteorite evidence. Assuming that fragments of Vesta were scattered from Vesta by giant planet migration, we conservatively estimate that at least two Vesta-like bodies (Vesta and the Magnya parent bodies) existed. From the oxygen isotopic, chemical, and petrologic evidence, we also conservatively estimate that seven Vesta-like bodies formed. Analyses of iron meteorites indicate that there may be as many as 23 Vesta-like bodies (Vesta, 10 magmatic iron groups, South Byron trio, Emsland/Mbosi duo, 10 ungrouped irons). This estimate from iron meteorites is most certainly an overestimation due to the existence of a number of non-HED crustal/mantle fragments that potentially originated from bodies with magmatic iron cores. Using our three estimates as a guide, we predict that there were ~10 Vesta-like bodies (including Vesta) that formed in the early solar system. Only Vesta remains intact with the others being disrupted early in solar system history.

T. H. Burbine↗
DOE OSTI2021

Topological Machine Learning Methods for Power System Responses to Contingencies

While deep learning tools, coupled with the emerging machinery of topological data analysis, are proven to deliver various performance gains in a broad range of applications, from image classification to biosurveillance to blockchain fraud detection, their utility in areas of high societal importance such as power system modeling and, particularly, resilience quantification in the energy sector yet remains untapped. To provide fast acting synthetic regulation and contingency reserve services to the grid while having minimal disruptions on customer quality of service, we propose a new topology-based system that depends on a neural network architecture for impact metric classification and prediction in power systems. This novel topology-based system allows one to evaluate the impact of three power system contingency types, in conjunction with transmission lines, transformers, and transmission lines combined with transformers. We show that the proposed new neural network architecture equipped with local topological measures facilitates more accurate classification of unserved load as well as the amount of unserved load. In addition, we are able to learn more about the complex relationships between electrical properties and local topological measurements on their simulated response to contingencies for the NREL-SIIP power system.

contingency analysis↗
DOE OSTI2023

Model predictive current profile control in tokamaks by exploiting spatially moving electron cyclotron current drives

The plasma control systems in next-generation tokamaks like ITER will balance competing control objectives to achieve the desired level of performance in advanced scenarios while preventing magnetohydrodynamic instabilities and disruptions. During normal tokamak operation, the points of incidence of the electromagnetic waves generated by the electron cyclotron heating and current drives (EC H&CDs) are usually fixed in space. However, the points of incidence can be modified in real-time by changing the angles of the mirrors that reflect the EC H&CD waves. Altering the points of incidence, in turn, varies the ability of the plasma control system to regulate a plasma property. For instance, changing the EC H&CD wave incidence location may place the power demands necessary to achieve a particular plasma target within saturation limits. Therefore, using the EC H&CD deposition location, which is related to the EC H&CD mirror angle, as a supplementary controllable variable may facilitate access to a given target scenario. However, active scenario-control algorithms have not been designed so far to fully exploit this capability in real time. In this work, a model predictive controller that can handle actuation locations as control inputs is developed. In particular, the controller is designed to regulate both the auxiliary powers and the EC H&CD deposition locations in a pre-defined optimal sense to achieve the control objective of attaining and sustaining a target current profile. Further, the proposed controller is tested for a DIII-D tokamak scenario in nonlinear simulations using the Control Oriented Transport SIMulator (COTSIM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

[Kelvin-Helmholtz instability in protostellar jets]

NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.

Stone, James↗
DOE OSTI2026

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗
NASA NTRS1993

(abstract) Cometary Particles as a Tracer of Jupiter's Stratospheric Circulation

The impact of fragments of comet Shoemaker-Levy 9 on Jupiter's atmosphere in July 1994 may provide an unprecedented opportunity to study Jupiter's stratospheric circulation. Recent calculations by Z. Sekanina predict that much of the comet material will be deposited in Jupiter's stratosphere. If so, and if the material is deposited in a confined region (10 000 km or less, horizontally) we can expect a situation analogous to an El Chichon or Pinatubo event for the terrestrial stratosphere. Initially the volatile material will be vaporized and will rapidly recondense. The large ice crystals and dust particles will rain out and be lost to the troposphere. The cloud of small particles which remain may have settling times of more than a year. These submicron to micron particles would probably be easily seen in methane filter images in the near-IR, and possibly in the ultraviolet. An observational program to monitor the dispersal of this cloud or clouds would reveal much about the nature of the circulation. Some predictions about the meridional evolution of the clouds can be made already, based on the meridional circulation model of West et al. unless the impact itself significantly disrupts the annual average circulation well after the initial transients die away.

comet Shoemaker-Levy 9 impact Jupiter stratosphere↗
NASA NTRS2009

Icy Satellites of Saturn: Impact Cratering and Age Determination

Saturn is the first giant planet to be visited by an orbiting spacecraft that can transmit large amounts of data to Earth. Crater counts on satellites from Phoebe inward to the regular satellites and ring moons are providing unprecedented insights into the origin and time histories of the impacting populations. Many Voyager-era scientists concluded that the satellites had been struck by at least two populations of impactors. In this view, the Population I impactors, which were generally judged to be comets orbiting the Sun, formed most of the larger and older craters, while Population II impactors, interpreted as Saturn-orbiting ejecta from impacts on satellites, produced most of the smaller and younger craters. Voyager data also implied that all of the ring moons, and probably some of the midsized classical moons, had been catastrophically disrupted and reaccreted since they formed. We examine models of the primary impactor populations in the Saturn system. At the present time, ecliptic comets, which likely originate in the Kuiper belt/scattered disk, are predicted to dominate impacts on the regular satellites and ring moons, but the models require extrapolations in size (from the observed Kuiper belt objects to the much smaller bodies that produce the craters) or in distance (from the known active Jupiter family comets to 9.5 AU). Phoebe, Iapetus, and perhaps even moons closer to Saturn have been struck by irregular satellites as well. We describe the Nice model, which provides a plausible mechanism by which the entire Solar System might have experienced an era of heavy bombardment long after the planets formed. We then discuss the three cratering chronologies, including one based upon the Nice model, that have been used to infer surface ages from crater densities on the saturnian satellites. After reviewing scaling relations between the properties of impactors and the craters they produce, we provide model estimates of the present-day rate at which comets impact, and catastrophically disrupt, the saturnian moons. Finally, we present crater counts on the satellites from two different groups. Many of the heavily cratered terrains appear to be nearly saturated, so it is difficult to infer the provenance of the impactors from crater counts alone. More large craters have been found on Iapetus than on any other satellite. Enceladus displays an enormous range of surface ages, ranging from the old mid-latitude plains to the extremely young South Polar Terrain. Cassini images provide some evidence for the reality of Population II. Most of the observed craters may have formed in one or more cataclysms, but more work is needed to determine the roles of heliocentric and planetocentric bodies in creating the craters.

Dones, L.↗
DOE OSTI2022

Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water–water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water–water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. Here, the ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water–water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average “hydrogen–acceptor-oxygen–donor-oxygen” angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Forest Structural Complexity and Biomass Predict First-Year Carbon Cycling Responses to Disturbance

The pre-disturbance vegetation characteristics that predict carbon (C) cycling responses to disturbance are not well known. To address this gap, we initiated the Forest Resilience Threshold Experiment, a manipulative study in which more than 3600 trees were stem girdled to achieve replicated factorial combinations of four levels (control, 45, 65, and 85% gross defoliation) of disturbance severity and two disturbance types (targeting upper or lower canopy strata). Applying a standardized stability framework in which initial C cycling resistance to disturbance was calculated as the first-year natural log response ratio of disturbance and control treatments, we investigated to what extent pre-disturbance levels of species diversity, aboveground woody biomass, leaf area index, and canopy rugosity—a measure of structural complexity—predict the initial responses of subcanopy light-saturated leaf CO 2 assimilation (A sat ), aboveground wood NPP (ANPP w ), and soil respiration (R s ) to phloem-disrupting disturbance. In the year following stem girdling, we found that above-ground C cycling processes, A sat and ANPP w , were highly resistant to increases in disturbance severity, while Rs resistance declined as severity increased. Disturbance type had no effect on first-year resistance. Pre-disturbance aboveground woody biomass, and canopy rugosity were positive predictors of ANPP w resistance and, conversely, negatively related to R s resistance. Subcanopy A sat resistance was not related to pre-disturbance vegetation characteristics. Stability of C uptake processes along with R s declines suggest the net C sink was sustained in the initial months following disturbance. Therefore, we conclude that biomass and complexity are significant, but not universal, predictors of initial C cycling resistance to disturbance. Moreover, our findings highlight the utility of standardized stability measures when comparing functional responses to disturbance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Contributions of anoxic microsites to soil carbon protection across soil textures

Anoxic microsites, zones of oxygen depletion in otherwise oxic soils, may slow soil C turnover. However, the abundance of anoxic microsites and their contribution to soil C protection is yet undefined. In this study, we determine the contribution of anoxic microsites to soil C protection in soils of three distinct textures (clay loam, loam, sandy loam) across a range of soil moistures. We examined the influence of soil oxygen supply by increasing oxygen content in the incubation atmosphere (oxygen enrichment) and through disaggregation. We attributed increases in CO 2 efflux to the aeration of anoxic microsites. The contribution of anoxic microsites to soil C protection increased with decreasing clay content. Clay loam CO 2 efflux was relatively unaffected by aeration. Moderately moist, loam soils had CO 2 effluxes that did not increase with oxygen enrichment but increased by 375% upon disaggregation. Sandy loam soil CO 2 efflux increased by 50–75% with oxygen enrichment and 250% with disaggregation. Geochemical and microbial data reveal that anoxic microsite abundance also increased with decreasing clay content. The proportion of acid extractable Fe present as Fe(II) increased with decreasing clay content, and methanogens were more abundant in loam and sandy loam soils. Our results suggest that oxygen demand, rather than supply, can regulate anoxic microsite formation and that anoxic protection of soil C can be diminished through physical disruption of soil structure. Our findings further illustrate that anoxic microsites should be included in conceptual models of soil C protection to avoid soil C loss and improve predictions of soil C response to disturbance.

58 GEOSCIENCES↗
DOE OSTI2026

Advances in 3D transient plasma dynamics and control through MHD and hybrid fluid-kinetic simulations with JOREK

Transient phenomena and their control are of high relevance in magnetic confinement fusion plasmas to guarantee a stable and safe plasma operation. Interpretative simulations can maximize the insights gained from experiments on present machines and predictive simulations can help in the preparation of design, mitigation techniques and operational scenarios for future devices. In this article, we provide an overview of recent advances and novel scientific results obtained with the 3D non-linear hybrid fluid-kinetic code JOREK, covering physics of plasma transients from the core to the scrape-off layer (SOL) both for tokamak and stellarator devices. Substantial progress was made in the physics understanding, model validation with experiments and experiment interpretation, thus, giving confidence for predictions to devices like DTT, ITER and DEMO. The topics addressed comprise a wide range: the edge physics of new operation scenarios and edge localized mode suppression; major disruptions with a focus on runaway electrons and vertical displacement events as well as disruption mitigation by shattered pellet injection; the physics mechanisms and operational limits of the flux pumping regime for sawtooth control; MHD limits of stellarators and work towards incorporating advanced edge/SOL/exhaust dynamics; continuing improvements of the code for more efficient hybrid simulations on conventional and accelerated high performance computing architectures.

disruptions↗