Search NASA⌕ Search

SEARCH · Search NASA

Results for “planetary data”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of second Si–O vibrational overtones/combinations on quantifying water in silicate and silica minerals using infrared spectroscopy, and an experimental method for its removal

Infrared spectroscopy (IR) is the most widely used analytical tool to quantify trace water in silicate and silica minerals. A prerequisite for highly accurate IR measurements of trace water is a good understanding of the effect of the second Si–O vibrational overtones/combination bands (2nd Si–O VOCBs) on the water peaks. Silicate and silica minerals can be divided as isolated (Q 0 ), paired (Q 1 ), ring (Q 2 ), chain (Q 1 or Q 2 ), sheet (Q 3 ) and framework (Q 4 ) structures according to the polymerization of their SiO 4 tetrahedral units, and the 2nd Si–O VOCBs of these different structural types attain different vibrational features which are expected to affect the water peaks to different extents. Here, we selected olivine (Q 0 ) and α-quartz (Q 4 ) as two endmember-like structural examples, performed extensive IR measurements on both pristine and heat-treated thin sections prepared for these two minerals, and explored the vibrational features of the 2nd Si–O VOCBs. In this work, we have found that the 2nd Si–O VOCBs are well separated from the water peaks in olivine, but severely overlap with the water peaks in α-quartz, confirming the different roles that the 2nd Si–O VOCBs play in quantifying trace water in silicate and silica minerals with different structural polymerizations. To remove the influence of the 2nd Si–O VOCBs (or any other species rather than water), an experimental protocol has been successfully developed, as approved by some fundamental equations and verified by the data of α-quartz in the literature. This development should lead to significant accuracy improvement in quantifying trace amounts of water in Earth and planetary materials.

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

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Nature of the bonded-to-atomic transition in liquid silica to TPa pressures

First-principles calculations and analysis of the thermodynamic, structural, and electronic properties of liquid SiO 2 characterize the bonded-to-atomic transition at 0.1–1.6 TPa and 10 4 –10 5 K (1–9 eV), the high-energy-density regime relevant to understanding planetary interiors. We find strong ionic bonds that become short-lived due to high kinetics during the transition, with sensitivity of the transition temperature to pressure, and our calculated Hugoniots agree with past experimental data. Furthermore, these results reconcile previous experimental and theoretical findings by clarifying the nature of the bond dissociation process in early Earth and “rocky” (oxide) constituents of large planets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

CORRIGENDUM to Figures S3-S4 of "Quasi-stationary Planetary Scale Waves in Modern Climate Models"

In addition to fixing the color-coding legend, the new figures plot Fourier phase in a more easily recognizable format: polar coordinates rather than Cartesian coordinates. The same phase data is displayed—indicating the position of the first ridge to the east of the prime meridian—but the new plots avoid unphysical jumps, e.g. from +180° to -180° as longitude increases across the International Dateline. A slight disadvantage of the new format is that each hemisphere requires a separate plot. The wavenumber-1 Northern Hemisphere phase plots give, in effect, a view from above the North Pole, with Eurasia at the top and North America at the bottom. Wavenumber-1 Southern Hemisphere plots give an analogous view from the South Pole. Plots for higher wavenumbers, however, wrap around at fractions of 360° (180°, 120°, ...) and cannot be interpreted as a geographical view.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Lithophile Element Budget of Earth's Core

Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO 2 , and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

Chidester, B. A.↗
DOE OSTI2023

A Novel Spatio-Temporal Regime Tracking Method for Impact Simulations

In this proposal, we present a novel method of tracking the rheological regimes activated during impact cratering events that will allow researchers to gain new insights into cratering mechanics. Rheology describes the stress-strain response of rocks to different conditions. Planetary impact cratering events often occur on too large of a scale to be feasibly captured in controlled experiments. Instead, these dynamic events are primarily studied using multi-physics codes equipped with complex material models that enable calculations of the impact event at scale. However, determining which physical processes are required for the problem of interest is challenging. Because the dominant rheological regimes change with space and time during crater formation, it is difficult to link numerical simulations with observable features of craters at the end of the event. The basis of this work is the implementation of numerical flags that track the activation of each rheological regime throughout impact simulations. We demonstrate this with the ‘Rock Model’ implemented in the CTH shock-physics code. We will use this rheology tracking method to ’zoom in’ on a specific region within an event and track the conditions the rock experiences over time. This work will develop community benchmarks to validate and distribute our implemented rheological models. We will focus on improving the melt models used by the planetary impact modeling community by developing an EOS-aware rheological transition from solid to melt. Through analysis of cell and tracer-particle based tracking data, we will study the effects of different rheologies on modeled outcomes, particularly on the volume and distributions of impacts melts. We will also use this method to link observable features with the rheological mechanisms responsible for them. The deliverables (peer-reviewed papers) from the proposed work are (i) tests of the implemented rheologic processes and demonstrations of the tracking flags; (ii) the first calculations of the spatio-temporal evolution of the dominant rheologies during impact cratering events; and (iii) application to delivery of impactor iron during basin-scale impacts.

58 GEOSCIENCES↗
DOE OSTI2020

OGLE-2017-BLG-0406: Spitzer Microlens Parallax Reveals Saturn-mass Planet Orbiting M-dwarf Host in the Inner Galactic Disk

We report the discovery and analysis of the planetary microlensing event OGLE-2017-BLG-0406, which was observed both from the ground and by the Spitzer satellite in a solar orbit. At high magnification, the anomaly in the light curve was densely observed by ground-based-survey and follow-up groups, and it was found to be explained by a planetary lens with a planet/host mass ratio of q=7.0×10{sup −4} from the light-curve modeling. The ground-only and Spitzer-“only” data each provide very strong one-dimensional (1D) constraints on the 2D microlens parallax vector π{sub E}. When combined, these yield a precise measurement of π{sub E} and of the masses of the host M{sub host}=0.56±0.07 M{sub ⊙} and planet M {sub planet} = 0.41 ± 0.05 M {sub Jup}. The system lies at a distance D {sub L} = 5.2 ± 0.5 kpc from the Sun toward the Galactic bulge, and the host is more likely to be a disk population star according to the kinematics of the lens. The projected separation of the planet from the host is a{sub ⊥}=3.5±0.3 au (i.e., just over twice the snow line). The Galactic-disk kinematics are established in part from a precise measurement of the source proper motion based on OGLE-IV data. By contrast, the Gaia proper-motion measurement of the source suffers from a catastrophic 10σ error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

NEMESIS: Exoplanet Transit Survey of Nearby M-dwarfs in TESS FFIs. I

In this work, we present an analysis of 33,054 M-dwarf stars, located within 100 parsecs, via the Transiting Exoplanet Survey Satellite (TESS) full-frame images (FFIs) of observed sectors 1–5. We present a new pipeline called NEMESIS, developed to extract detrended photometry, and to perform transit searches of single-sector data in TESS FFIs. As many M-dwarfs are faint, and are not observed with a two-minute cadence by TESS, FFI transit surveys can provide an empirical validation of how many planets are missed, using the 30-minute cadence data. In this work, we detect 183 threshold crossing events, and present 29 candidate planets for sectors 1–5, 24 of which are new detections. Our sample contains orbital periods ranging from 1.25 to 6.84 days, and planetary radii from 1.26 to 5.31 R {sub ⊕}. With the addition of our new planet candidate detections, along with detections previously observed in sectors 1–5, we calculate an integrated occurrence rate of 2.49 ± 1.58 planets per star, for the period range ∈ [1, 9] days, and planet radius range ∈ [0.5,11] R {sub ⊕}. We project an estimated yield of 122 ± 11 transit detections of nearby M-dwarfs. Of our new candidates, 23 have signal-to-noise ratios >7, transmission spectroscopy metrics >38, and emission spectroscopy metrics >10. We present all of our data products for our planet candidates via the Filtergraph data visualization service, located at https://filtergraph.com/NEMESIS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Causes of and Solutions to Wind Speed Bias in NREL's 2020 Offshore Wind Resource Assessment for the California Pacific Outer Continental Shelf

This report provides the results of a detailed analysis into the causes of high wind speed bias in the 20-year wind resource data set for offshore California the National Renewable Energy Laboratory (NREL) released in 2020, herein called CA20. The data set was developed using the state-of-the-art Weather Research and Forecasting (WRF) model. Notably, no floating lidars were available at the time in offshore California to validate offshore hub-height wind speeds. In late 2020, the Pacific Northwest National Laboratory (PNNL) deployed two floating lidars in the California outer continental shelf (OCS), near the Bureau of Ocean Energy Management (BOEM) call areas of Humboldt and Morro Bay. Using these observations through 2021, NREL found considerable bias in modeled hub-height winds at both locations: up to +2 m/s at Humboldt over a 6-month period, and up to +1 m/s at Morro Bay over a one-year period. Upon the discovery of this bias, the Department of Energy (DOE) and BOEM funded NREL and PNNL to investigate the causes of, impacts of, and solutions to the bias in the CA20 data set. This report summarizes the findings of this research. We first investigated whether different WRF model setups could lead to reduced bias. We found that the choice of planetary boundary layer (PBL) scheme - which controls the vertical turbulent mixing of momentum, heat, and moisture in the lowermost part of the atmosphere - greatly affected hub-height wind speeds in the region. Specifically, switching from the Mellor-Yamada-Nakanishi-Niino (MYNN) scheme used in CA20 (and widely used across a range of operational and research weather models) to the less common Yonsei University (YSU) scheme nearly eliminated the bias at both the Humboldt and Morro Bay lidar locations. The large discrepancy between the MYNN- and YSU-modeled hub-height winds pointed towards the role of atmospheric stability. In general, PBL schemes agree well in conditions of high turbulence and mixing, normally referred to as "unstable" conditions. By contrast, PBL schemes start to diverge in "stable" conditions, where turbulence is low and thermal stratification (i.e., higher temperature air sitting on top of colder air) greatly suppresses vertical mixing. Under such conditions, winds aloft can decouple from surface effects and greatly accelerate, causing high wind speeds at hub-height and frequent low-level jets (LLJs). We determined that these stable conditions are in fact dominant in offshore California. The region is characterized by moderate-to-extreme stable stratification with a LLJ on average around 200 meters above sea-level. To our knowledge, no wind energy area globally has as strongly stable stratification as offshore California. Under these extreme conditions, we determined that the MYNN scheme models higher stability than YSU, resulting in less vertical turbulent mixing than YSU, allowing for the acceleration of hub-height winds, more intense LLJs, and higher-amplitude inertial oscillations. Using surface observations, we found that MYNN overestimates near-surface stability, whereas YSU tends to model stability better. We then considered several short-term case studies to assess additional meteorological drivers of the bias at Humboldt. We found that during synoptic scale northerly flows driven by the North Pacific High and inland thermal low, a coastal warm bias in the MYNN case studies contributes to the modeled wind speed bias by altering the boundary layer thermodynamics via a thermal wind mechanism. Given the strong performance of the YSU-based runs in offshore California, NREL has produced and published an updated version of the CA20 data set with YSU as the PBL scheme. This updated data set is now part of NREL's 2023 National Offshore Wind (NOW-23) data set, which covers all the U.S. offshore waters. The development and final validation of the NOW-23 data set in offshore California is documented in this report.

17 WIND ENERGY↗
DOE OSTI2021

A Machine Learns to Predict the Stability of Highly Diverse Multi-planetary Systems

Machine learning has proven to be an invaluable tool for characterizing the stability of planets in simplified planetary systems. In this work, we investigate the performance of a machine learning classifier on tightlypacked systems containing a rich diversity of planets, from Earths to Jupiters. Using information derived from short numerical simulations about a planet’s early orbital evolution and its relationship with the most massive planets in the system, we train a random forest classifier to predict instability with a > 88 percent accuracy. Our classifier relies on relative planet masses and the standard deviation of eccentricity for much of its predictive power. Most misclassified planets lie along a multi-dimensional boundary between stable and unstable planets, indicating that their early orbital evolution is ambiguous. The major reason for misclassification in this work is timescale: because our classifier uses information from only the first 137 years of simulation data, it is blind to late time interactions that cause or prevent instability. Machine learning methods like those utilized in this work provide powerful tools to complement numerical simulations across a wide range of planetary architectures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Effect of Surface Heating Heterogeneity on Boundary Layer Height and Its Dependence on Background Wind Speed

The planetary boundary layer height (PBLH) is a fundamental variable of the planetary boundary layer (PBL). The effect of surface heterogeneity is challenging in the PBL parameterization for numerical weather and climate models. Here we use large eddy simulation data to investigate the impact of surface heating heterogeneity on the PBLH spatial variations and mean values over a domain of typical mesoscale model grid size. It is found that this impact depends on background wind condition. Variance decomposition and Fourier spectra both reveal that for the fully developed sheared convective PBL, surface heterogeneity contributes little to the spatial variation of PBLH, except for the case with heterogeneity scale of 14.4 km. This suggests that surface heterogeneity with length scales less than ~10 km, which is typical for mesoscale modeling, does not induce significant subgrid spatial variation of PBLH. The domain-averaged PBLH is found to generally increase as surface heterogeneity scale decreases, until reaching the scale on the order of PBLH. The findings can be used to inform the parameterization development of PBL processes over heterogeneous surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High-precision shock equation of state measurements for metallic fluid carbon between 15 and 20 Mbar

Diamond is an efficient ablator material to convert the energy of high-power giant lasers into ablation pressure with applications for High-Energy-Density (HED) science, planetary science, and Inertial Confinement Fusion (ICF) research at the National Ignition Facility (NIF). Unfortunately, current theoretical equation of state models cannot reproduce all the observed experimental data in the multi-megabar regime particularly relevant for HED and ICF research. New experimental data on the behavior of carbon at extreme pressures and temperatures are, therefore, essential to improve our predictive capability to design and analyze dynamic compression experiments for HED or ICF research and build improved equation of state models in the future. In this paper, we report high-precision laser-driven shock compression measurements on diamond single crystals at the Omega Laser Facility. Using ultrafast Doppler optical Velocimetry Interferometer System for Any Reflector (VISAR) to track the leading shock front and a quartz plate as an in situ reference, we obtain relative pressure-density shock equation-of-state measurements between 15 and 20 Mbar with an impedance-matching procedure. We also report shock-and-release measurements in a spherical geometry at the NIF. The new data provide tight constraints on the compressibility of warm dense carbon along the Hugoniot of full density diamond, allowing us to discriminate between existing theoretical equation-of-state models. We find that both LLNL LEOS 9061 and LANL Sesame 7835 models capture well the shock compressibility in the explored range. LANL Sesame 7835 also reproduces well the observed shock-and-release behavior of diamond near 10–20 Mbar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Toward a Climate OSSE Framework for Satellite Mission Design

The rich history of observing system simulation experiments (OSSEs) does not yet include a well-established framework for using climate models. The need for a climate OSSE is triggered by the need to quantify the value of a particular measurement for reducing the uncertainty in climate predictions, which differ from numerical weather predictions in that they depend on future atmospheric composition rather than the current state of the weather. However, both weather and climate modeling communities share a need for motivating major observing system investments. Here, we outline a new framework for climate OSSEs that leverages the use of machine learning to calibrate climate model physics against existing satellite data. We demonstrate its application using NASA’s GISS-E3 model to objectively quantify the value of potential future improvements in spaceborne measurements of Earth’s planetary boundary layer. A mature climate OSSE framework should be able to quantitatively compare the ability of proposed observing system architectures to answer a climate-related question, thus offering added value throughout the mission design process, which is subject to increasingly rapid advances in instrument and satellite technology. Technical considerations include selection of observational benchmarks and climate projection metrics, approaches to pinpoint the sources of model physics uncertainty that dominate uncertainty in projections, and the use of instrument simulators. Community and policy-making considerations include the potential to interface with an established culture of model intercomparison projects and a growing need to economically assess the value-driven efficiency of social spending on Earth observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

TESS Discovery of a Super-Earth and Three Sub-Neptunes Hosted by the Bright, Sun-like Star HD 108236

We report the discovery and validation of four extrasolar planets hosted by the nearby, bright, Sun-like (G3V) star HD 108236 using data from the Transiting Exoplanet Survey Satellite (TESS). We present transit photometry, reconnaissance, and precise Doppler spectroscopy, as well as high-resolution imaging, to validate the planetary nature of the objects transiting HD 108236, also known as the TESS Object of Interest (TOI) 1233. The innermost planet is a possibly rocky super-Earth with a period of 3.79523{sub −0.00044}{sup +0.00047} days and has a radius of 1.586 ± 0.098 R {sub ⊕}. The outer planets are sub-Neptunes, with potential gaseous envelopes, having radii of 2.068{sub −0.091}{sup +0.10} R {sub ⊕}, 2.72 ± 0.11 R {sub ⊕}, and 3.12{sub −0.12}{sup +0.13} R {sub ⊕} and periods of 6.20370{sub −0.00052}{sup +0.00064} days, 14.17555{sub −0.0011}{sup +0.00099} days, and 19.5917{sub −0.0020}{sup +0.0022} days, respectively. With V and K {sub s} magnitudes of 9.2 and 7.6, respectively, the bright host star makes the transiting planets favorable targets for mass measurements and, potentially, for atmospheric characterization via transmission spectroscopy. HD 108236 is the brightest Sun-like star in the visual (V) band known to host four or more transiting exoplanets. The discovered planets span a broad range of planetary radii and equilibrium temperatures and share a common history of insolation from a Sun-like star (R {sub ⋆} = 0.888 ± 0.017 R {sub ⊙}, T {sub eff} = 5730 ± 50 K), making HD 108236 an exciting, opportune cosmic laboratory for testing models of planet formation and evolution.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Seismic constraints from a Mars impact experiment using InSight and Perseverance

NASA’s InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission has operated a sophisticated suite of seismology and geophysics instruments on the surface of Mars since its arrival in 2018. On 18 February 2021, we attempted to detect the seismic and acoustic waves produced by the entry, descent and landing of the Perseverance rover using the sensors onboard the InSight lander. Similar observations have been made on Earth using data from both crewed and uncrewed spacecraft, and on the Moon during the Apollo era, but never before on Mars or another planet. This was the only seismic event to occur on Mars since InSight began operations that had an a priori known and independently constrained timing and location. It therefore had the potential to be used as a calibration for other marsquakes recorded by InSight. Here we report that no signal from Perseverance’s entry, descent and landing is identifiable in the InSight data. Nonetheless, measurements made during the landing window enable us to place constraints on the distance–amplitude relationships used to predict the amplitude of seismic waves produced by planetary impacts and place in situ constraints on Martian impact seismic efficiency (the fraction of the impactor kinetic energy converted into seismic energy).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multiple Explanations for the Single Transit of KIC 5951458 Based on Radial Velocity Measurements Extracted with a Novel Matched-template Technique

Planetary systems that show single-transit events are a critical pathway to increasing the yield of long-period exoplanets from transit surveys. From the primary Kepler mission, KIC 5951458 b (Kepler-456b) was thought to be a single-transit giant planet with an orbital period of 1310 days. However, radial velocity (RV) observations of KIC 5951458 from the HIRES instrument on the Keck telescope suggest that the system is far more complicated. To extract precise RVs for this V≈13 star, we develop a novel matched-template technique that takes advantage of a broad library of template spectra acquired with HIRES. We validate this technique and measure its noise floor to be 4–8 m s{sup −1} (in addition to internal RV error) for most stars that would be targeted for precision RVs. For KIC 5951458, we detect a long-term RV trend that suggests the existence of a stellar companion with an orbital period greater than a few thousand days. We also detect an additional signal in the RVs that is possibly caused by a planetary or brown dwarf companion with mass in the range of 0.6–82 M {sub Jup} and orbital period below a few thousand days. Curiously, from just the data on hand, it is not possible to determine which object caused the single “transit” event. We demonstrate how a modest set of RVs allows us to update the properties of this unusual system and predict the optimal timing for future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Search for L4 Earth Trojan Asteroids Using a Novel Track-before-detect Multiepoch Pipeline

We report Earth Trojan asteroids are an important but elusive population that co-orbit with Earth at the L4 and L5 Lagrange points. There is only one known, but a large population is theoretically stable and could provide insight into our solar system’s past and present, as well as planetary defense. In this paper, we present the results of an Earth Trojan survey that uses a novel shift-and-stack detection method on 2 nights of data from the Dark Energy Camera. We find no new Earth Trojan asteroids. We calculate an upper limit on the population that is consistent with previous searches despite much less sky coverage. Additionally, we elaborate on previous upper-limit calculations using current asteroid population statistics and an extensive asteroid simulation to provide the most up-to-date population constraints. We find an L4 Earth Trojan population of N ET < 1 for H = 13.93, N ET < 10 for H = 16, and N ET < 938 for H = 22.

79 ASTRONOMY AND ASTROPHYSICS↗