SEARCH · Engineering Papers

Results for “planetary dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Remote sensing of planetary boundary layer from ground and space: Structure, thermodynamics, and boundary layer clouds

The planetary boundary layer (PBL) is a critical interface between Earth's surface and atmosphere, influencing atmospheric convection, weather patterns, and air quality. Recognized by the 2017 National Academies of Sciences, Engineering, and Medicine Earth Science Decadal Survey as an Incubation Targeted Observable, high-quality and effective observation of the PBL has become a priority. This paper reviews surface-based and satellite remote sensing techniques for characterizing PBL features, including PBL height (PBLH), boundary layer thermodynamics, turbulence, and boundary layer clouds (BLCs). These elements are treated as interconnected aspects of the PBL system, while recognizing that different instruments retrieve different physical manifestations of the PBL. The review summarizes recent advances and limitations in Micro-Pulse Lidar (MPL), Doppler lidar, Raman lidar, Differential Absorption Lidar (DIAL), ceilometers, wind profilers, GNSS Radio Occultation, radar, and hyperspectral sounders. These observations help trace dynamic processes within the PBL and link PBL structure to broader weather and climate processes. The review also highlights persistent observational gaps over oceans, remote land, and polar regions, where continuous surface-based profiling is sparse. By comparing instrument capabilities, practical limitations, and interpretation issues, this review emphasizes the need for integrated remote sensing approaches and careful definition of the retrieved PBL quantities.

Boundary layer clouds↗

DOE OSTI2023

Hydrodynamics and Survivability during Post-main-sequence Planetary Engulfment

The engulfment of substellar bodies (SBs), such as brown dwarfs and planets, by giant stars is a possible explanation for rapidly rotating giants, lithium-rich giants, and the presence of SBs in close orbits around subdwarfs and white dwarfs. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an engulfed SB. We model the SB as a rigid body with a reflective surface because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram-pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory inside the star. We find that engulfment can increase the luminosity of a 1 M ⊙ star by up to a few orders of magnitude. The time for the star to return to its original luminosity is up to a few thousand years when the star has evolved to ≈10 R ⊙ and up to a few decades at the tip of the red giant branch (RGB). No SBs can eject the envelope of a 1 M ⊙ star before it evolves to ≈10 R ⊙ if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as ≈10 MJup can eject the envelope at the tip of the RGB. The numerical framework we introduce here can be used to study planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Density of Fe-Ni-C Liquids at High Pressures and Implications for Liquid Cores of Earth and the Moon

The presence of light elements in the metallic cores of the Earth, the Moon, and other rocky planetary bodies has been widely proposed. Carbon is among the top candidates in light of its high cosmic abundance, siderophile nature, and ubiquity in iron meteorites. It is, however, still controversial whether carbon-rich core compositional models can account for the seismic velocity observations within the Earth and lunar cores. In this paper we report the density and elasticity of Fe 90 Ni 10 -3 wt.% C and Fe 90 Ni 10 -5 wt.% C liquid alloys using synchrotron-based X-ray absorption experiments and first-principles molecular dynamics simulations. Our results show that alloying of 3 wt.% and 5 wt.% C lowers the density of Fe 90 Ni 10 liquid by ~2.9-3.1% at 2 GPa, and ~3.4-3.6% at 9 GPa. More intriguingly, our experiments and simulations both demonstrate that the bulk moduli of the Fe-Ni-C liquids are similar to or slightly higher than those of Fe-Ni liquids. Thus, the calculated compressional velocity ($v_p$) of Fe-Ni-C liquids are higher than that of pure Fe-Ni alloy, promoting carbon as a possible candidate to explain the elevated v p in the Earth's outer core. However, the values and slopes of both density and v p of the studied two Fe-Ni-C liquids do not match the outer core seismic models, suggesting that carbon may not be the sole principal light element in Earth's outer core. The high $v_p$ of Fe-Ni-C liquids does not match the presumptive $v_p$ of the lunar outer core well, indicating that carbon is less likely to be its dominant light element.

58 GEOSCIENCES↗

DOE OSTI2021

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Propagation and Periodicity of Mars's Northern Annular Mode Modulates the Dust Cycle

Abstract We document the propagation of annular modes—zonally symmetric patterns of variability—in Mars's atmosphere using a reanalysis dataset. Mars's Northern Annular Mode (MNAM) sees anomalies of zonal‐mean zonal wind emerge near the subtropics and migrate poleward with a period of 150 days, similarly to Earth's Southern Annular Mode. The mechanism of propagation involves the interaction of the two leading empirical orthogonal functions that define the MNAM. Moreover, the propagation encourages alternating bands of surface wind stress to migrate polewards with a 150‐day period. In addition, a 150‐day periodicity in anomalous column dust optical depth most likely emerges in response to extrema of the MNAM. The combination of the impact of the MNAM's internally forced periodicity on the surface wind stress and the seasonal cycle may contribute to the inter‐annual variability of global dust events, as suggested by a Monte Carlo estimate that correctly approximates the observed incidence of global dust events.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

The structure of liquid carbon elucidated by in situ X-ray diffraction

Carbon has a central role in biology and organic chemistry, and its solid allotropes provide the basis of much of our modern technology. However, the liquid form of carbon remains nearly uncharted, and the structure of liquid carbon and most of its physical properties are essentially unknown. But liquid carbon is relevant for modelling planetary interiors and the atmospheres of white dwarfs, as an intermediate state for the synthesis of advanced carbon materials, inertial confinement fusion implosions, hypervelocity impact events on carbon materials and our general understanding of structured fluids at extreme conditions. Here we present a precise structure measurement of liquid carbon at pressures of around 1 million atmospheres obtained by in situ X-ray diffraction at an X-ray free-electron laser. Our results show a complex fluid with transient bonding and approximately four nearest neighbours on average, in agreement with quantum molecular dynamics simulations. The obtained data substantiate the understanding of the liquid state of one of the most abundant elements in the universe and can test models of the melting line. The demonstrated experimental abilities open the path to performing similar studies of the structure of liquids composed of light elements at extreme conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Biological particles and aerosol-cloud interactions in the Southern Great Plains (Final Report)

The United States Southern Great Plains (SGP) is the genesis for much of the warm season precipitation in the central and eastern United States. While atmospheric thermodynamics and large-scale dynamics play an important role in formation of precipitation, precipitation mechanisms are also sensitive to aerosols. Aerosols can suppress precipitation as cloud condensation nuclei (CCN) or act as ice nucleating particles (INP) in deep convective clouds. This project focused on understanding the role of primary biological aerosol particles (PBAP) in the region and its influence on cloud formation. Specifically, we focused on biological aerosol in the form of pollen, one type of PBAP that is emitted in large yet variable quantities from vegetation in the mid latitudes. Field observations provide evidence of pollen in the planetary boundary layer and pollen components in cloud droplets and fine particulate matter. Further, pollen grains can easily rupture when wet, forming smaller, sub pollen particles with sizes less than one micron. We evaluated the potential for PBAP events at the Department of Energy (DoE) SGP Atmospheric Radiation Measurement (ARM) research facility due to the rich dataset available. We also developed model simulations that accounted for pollen emission and the generation of sub pollen particles, which have been shown to be both cloud condensation nuclei and ice nucleating particles. The proposed work was designed to address the following questions: 1. What are the physical and chemical signatures of biological aerosol such as pollen in SGP ARM measurements? 2. What is the role of pollen-derived particles on deep convection and precipitation in the Central US? Using measurements from the DoE SGP ARM site and recent airborne campaigns, we evaluated the signatures of pollen and pollen-derived aerosols on optical properties, cloud properties and precipitation (Subba et al., 2021). We identified pollen-driven events over the data record and used the Weather Research and Forecasting Model with fully coupled chemistry (WRF-Chem) to conduct chemically realistic simulations of aerosols during and summer mesoscale convective events in the Southern Great Plains (Subba et al., in review). This work improved our understanding of the role of biological aerosol on clouds and precipitation in the Central United States and placed these results in context with anthropogenically-driven processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Biological particles and aerosol-cloud interactions in the Southern Great Plains

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

ARM SGP PBLH and MLH datasets from Raman lidar and Doppler lidar

The planetary boundary layer (PBL) plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including the daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary for studying PBL characterization and processes. Here, new approaches are developed to provide reliable PBLH and MLH estimates to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) User Facility. Compared to widely used lidar aerosol measurements for PBLH, WVMR is a better tracer for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to reconstruct 2D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm-season RL and ML measurements at the SGP site for five years to study warm-season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in a warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable data set of PBLH and MLH values will be valuable for studying PBL processes, model evolution, and PBL parameterization improvements. The MLH dataset includes the MLH in values of km above ground level. The PBLH dataset includes the PBLH in values of km above ground level, along with a flag ("situation_PBLH") to determine the state of the PBL (1 = Cloudy Condition, 2 = Stable Layer, 3 = Multi-layer WVMR structure, 4 = Well-Mixed PBL, 5 = A de-coupled layer, 6 = Other).

mixing layer height↗

DOE OSTI2023

The microphysics of the warm-rain and ice crystal processes of precipitation in simulated continental convective storms

Precipitation in clouds can form by either warm-rain or ice crystal processes, referred to as warm and cold formation pathways, respectively. Here, we investigate the warm and cold pathway contributions to surface precipitation in simulated continental convective storms. We analyze three contrasting convective storms that are cold-based, slightly warm-based and very warm-based. We apply tracer-tagging techniques in our aerosol-cloud model to determine simulated microphysical pathways that lead to precipitation. We find cold components of graupel and rain mass were higher than warm components in cold- and slightly warm-based clouds. By contrast, in very warm-based clouds nearly 80% of surface precipitation was formed via warm-rain processes. Lowering of cloud base altitude to levels about 10–20 K warmer switched surface precipitation to being mostly warm, due to enhanced moisture content in the planetary boundary layer and larger cloud droplets aloft intensifying raindrop freezing. Our simulations indicate that warm and cold processes co-exist in any storm and the balance between them is determined by cloud base temperature and solute aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Capture and migration of Jupiter and Saturn in mean motion resonance in a gaseous protoplanetary disc

Here, we study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets’ orbits, Δa SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inwards or outwards). We also provide an assessment of the planet’s orbital dynamics at different epochs of Saturn’s growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Δa SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the time-scales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low-probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

A new framework to map fine resolution cropping intensity across the globe: Algorithm, validation, and implication

We report accurate estimation of cropping intensity (CI), an indicator of food production, is well aligned with the ongoing efforts to achieve sustainable development goals (SDGs) under diminishing natural resources. The advancement in satellite remote sensing provides unprecedented opportunities for capturing CI information in a spatially continuous manner. However, challenges remain due to the lack of generalizable algorithms for accurately and efficiently mapping global CI with a fine spatial resolution. In this study, we developed a 30-m planetary-scale CI mapping framework with the reconstructed time series of Normalized Difference Vegetation Index (NDVI) from multiple satellite images. Using a binary crop phenophase profile indicating growing and non-growing periods, we estimated pixel-by-pixel CI by enumerating the total number of valid cropping cycles during the study years. Based on the Google Earth Engine cloud computing platform, we implemented the framework to estimate CI during 2016–2018 in eight geographic regions across continents that are representative of global cropping system diversity. Comparison with PhenoCam network data in four cropland sites suggests that the proposed framework is capable of capturing the seasonal dynamics of cropping practices. Spatially, overall accuracies based on validation samples range from 80.0% to 98.9% across different regions worldwide. Regarding the CI classes, single cropping systems are associated with more robust and less biased estimations than multiple cropping systems. Finally, our CI estimates reveal high agreement with two widely used land surface phenology products, including Vegetation Index and Phenology V004 (VIP4) and Moderate Resolution Imaging Spectroradiometer Land Cover Dynamics (MCD12Q2), meanwhile providing much more spatial details. Due to its robustness, the developed CI framework can be potentially generalized to produce global fine resolution CI products for food security and other applications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

An Open Benchmark of One Million High-Fidelity Cislunar Trajectories

Cislunar space spans from geosynchronous altitudes to beyond the Moon and will underpin future exploration, science, and security operations. We describe and release an open dataset of one million numerically propagated cislunar trajectories generated with the open-source Space Situational Awareness Python package (SSAPy). The model includes high-degree Earth/Moon gravity, solar gravity, and Earth/Sun radiation pressure; other planetary gravities are omitted by design for computational efficiency. Initial conditions uniformly sample commonly used osculating-element ranges, and each trajectory is propagated for up to six years under a single, fixed start epoch. The dataset is intended as a reusable benchmark for method development (e.g., space domain awareness, navigation, and machine-learning pipelines), a reference library for statistical studies of orbit families, and a starting point for community-driven extensions (e.g., alternative epochs). We report empirically observed stability trends (e.g., a band near ~5 GEO and persistence of some co-orbital classes including L4/L5 librators) as dataset descriptors rather than new dynamical results. The chief contribution is the scale, fidelity, organization (CSV/HDF5 with full state time series and metadata), and open availability, which together lower the barrier to comparative and data-driven studies in the cislunar regime.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

BLADE: An Automated Framework for Classifying Light Curves from the Center for Near-Earth Object Studies Fireball Database

Fireballs (bolides) are high-energy luminous phenomena produced when meteoroids and small asteroids enter Earth’s atmosphere at hypersonic speeds, often resulting in fragmentation or complete disintegration accompanied by significant energy release. The resulting bolide light curves capture temporal brightness variations as these objects traverse increasingly dense atmospheric layers, providing essential information on meteoroid entry dynamics, fragmentation behavior, and atmospheric energy deposition processes. The Center for Near-Earth Object Studies’ (CNEOS) continuously expanding fireball database offers a globally comprehensive archive of bolide events, including light curves and associated metadata. Events associated with infrasound detections allow direct correlations between acoustic signatures and light curve features, therefore enabling detailed analyses of fragmentation dynamics and energy deposition. Here, we introduce Bolide Light-curve Analysis and Discrimination Explorer (BLADE), a robust and high-fidelity framework specifically designed to analyze bolide light curves for objects detected from space. BLADE incorporates a processing pipeline integrating Savitzky–Golay filtering, prominence-based peak detection, and gradient analysis, enabling systematic identification and classification of fragmentation events and their associated energy release characteristics. Preliminary results demonstrate that BLADE reliably distinguishes distinct bolide behaviors, providing an objective, scalable methodology for characterization and analysis of large bolide light curve data sets. This foundational work establishes a novel pathway for advanced bolide research, with promising applications in planetary defense and global atmospheric monitoring. Future research should adopt an integrative approach combining CNEOS optical data with complementary infrasound measurements, further clarifying relationships between bolide energy deposition and acoustic signatures, thus refining our understanding of meteoroid and asteroid atmospheric entry processes.

Asteroids↗

DOE OSTI2023

Plasma Wave and Particle Dynamics During Interchange Events in the Jovian Magnetosphere Using Juno Observations

Abstract Interchange instability is known to drive fast radial transport of particles in Jupiter's inner magnetosphere. Magnetic flux tubes associated with the interchange instability often coincide with changes in particle distributions and plasma waves, but further investigations are required to understand their detailed characteristics. We analyze representative interchange events observed by Juno, which exhibit intriguing features of particle distributions and plasma waves, including Z‐mode and whistler‐mode waves. These events occurred at an equatorial radial distance of ∼9 Jovian radii on the nightside, with Z‐mode waves observed at mid‐latitude and whistler‐mode waves near the equator. We calculate the linear growth rate of whistler‐mode and Z‐mode waves based on the observed plasma parameters and electron distributions and find that both waves can be locally generated within the interchanged flux tube. Our findings are important for understanding particle transport and generation of plasma waves in the magnetospheres of Jupiter and other planetary systems.

Geology↗