Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained 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.

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

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A parametric sensitivity study for single-stage-to-orbit hypersonic vehicles using trajectory optimization

The class of hypersonic vehicle configurations with single stage-to-orbit (SSTO) capability reflect highly integrated airframe and propulsion systems. These designs are also known to exhibit a large degree of interaction between the airframe and engine dynamics. Consequently, even simplified hypersonic models are characterized by tightly coupled nonlinear equations of motion. In addition, hypersonic SSTO vehicles present a major system design challenge; the vehicle's overall mission performance is a function of its subsystem efficiencies including structural, aerodynamic, propulsive, and operational. Further, all subsystem efficiencies are interrelated, hence, independent optimization of the subsystems is not likely to lead to an optimum design. Thus, it is desired to know the effect of various subsystem efficiencies on overall mission performance. For the purposes of this analysis, mission performance will be measured in terms of the payload weight inserted into orbit. In this report, a trajectory optimization problem is formulated for a generic hypersonic lifting body for a specified orbit-injection mission. A solution method is outlined, and results are detailed for the generic vehicle, referred to as the baseline model. After evaluating the performance of the baseline model, a sensitivity study is presented to determine the effect of various subsystem efficiencies on mission performance. This consists of performing a parametric analysis of the basic design parameters, generating a matrix of configurations, and determining the mission performance of each configuration. Also, the performance loss due to constraining the total head load experienced by the vehicle is evaluated. The key results from this analysis include the formulation of the sizing problem for this vehicle class using trajectory optimization, characteristics of the optimal trajectories, and the subsystem design sensitivities.

Lovell, T. Alan↗
NASA NTRS2020

The bioclimatic extent and pattern of the cold edge of the boreal forest: the circumpolar taiga-tundra ecotone

Current configurations of forest structure at the cold edge of the boreal may help understand the future of ecosystem functioning in high northern latitudes. The circumpolar biome boundary at the boreal (taiga) forest and tundra interface is an ecological transition zone (taiga-tundra ecotone; TTE) experiencing changes that affect its forest structure. We accounted for the TTE's horizontal forest structure with an estimate of its extent and pattern as represented by tree canopy cover (TCC). We quantified TCC patterns with an algorithm that describes its spatial gradient, and summarized landscape patterns of structure to represent heterogeneity, capturing abrupt, diffuse, and uniform forest at mesoscales. We used these landscape patterns to constrain the spatial extent of sparse and open canopy forest, and non-forest (forest-adjacent) edge that defines the TTE extent. The resulting map of the TTE extent is based on forest structure spatial patterns resolved at 30 m, highlights structural variability across landscapes, and helps distinguish tundra from boreal domains. We classified 14 594 landscapes as those associated with the TTE within a circumpolar bioclimatic envelope (11.575 million km2), where 44.83% of the area of these landscapes were forest and non-forest edge, yet 36.43% contributed to the TTE extent. We report the overall extent of the TTE (3.032 million km2) across North America and Greenland (53%), and Eurasia (47%), where 0.697 million km2 is non-forest edge, 0.549 million km2 is sparse forest, and 1.787 million km2 is open canopy forest. Diffuse forest landscapes dominate the TTE (79%), and abrupt landscapes (~19%) indicate portions of the TTE where sparse forest and non-forest edge are the prevailing structural patterns. This account of the TTE quantifies the area of the cold edge of the boreal forest where previous global estimates show high discrepancies, and can help target monitoring and prediction of circumpolar dynamics.

Paul M Montesano↗
NASA NTRS2001

Diagnosing Tropical Divergent Circulations From Satellite-Derived Diabatic Forcing and Variational Integral Constraints

Large-scale divergent circulations are part of the atmospheric dynamic response to diabatic heating from condensation, radiative processes, and surface energy fluxes. Vertical motion and the associated divergent wind is thus intimately tied to the hydrologic cycle and the global beat balance. Despite its importance, the divergent circulation is too small in comparison to the rotational flow to measure directly with any accuracy. Vertical motions are recovered diagnostically from reanalyses and, as such, are subject to shortcomings in model physics, numerics, and data availability. While reanalysis estimates of tropical divergent circulations are much improved over those from the Global Weather Experiment, there are still substantial differences between products from operational centers. This is because these circulations are still forced largely by model physics and only secondarily by observations. In order to produce a refined estimate of tropical divergence and its interannual variability we have used a number of remotely-sensed data sets along with variational constraints to improve upon reanalysis estimates. Among these are: precipitation from SSM/I and GPCP, TOVS Path-A vertical cloud distributions; ISCCP radiative cooling rates; TOA radiative fluxes from ERBS, surface radiative fluxes from the SRB project, and surface latent and sensible flux estimates from SSM/I. The TOVS Path-A data constrain the divergent outflow in precipitating regions to have the same vertical structure as observed cloudiness. Using integral constraints for moisture, heat, and mass balance, we retrieve consistent divergent wind flows. We examine the ability of this type analysis to capture regional details of ENSO related perturbations to the divergent wind and associated tropical energy balance. Precipitation from satellite is found to be the major constraint in supplying this horizontal structure. We also consider the ability of this analysis to quantify integrated (land vs. ocean) fluctuations in the global monsoonal flow and energy balance. This depends, of course, on the error characteristics of the data constraints. Some estimate of these properties and their relative importance is provided.

Robertson, Franklin R.↗
NASA NTRS2022

The breakup of a long-period comet is not a likely match to the Chicxulub impactor

Since the discovery of Ir in the clay layer at the K-Pg boundary [1], scientists have sought to constrain the origin of the extraterrestrial impactor that triggered the end-Cretaceous mass extinction of the non-avian dinosaurs and other species. While the first proposal was for an asteroid, for a while some theories invoked a cometary impactor to explain perceived periodicities in mass extinctions. Such models have long been disfavored by the mass of Ir in the layer, inferred to be2.0−2.8×1011g. The size of Chicxulub crater leads to an estimated asteroid impactor diameter,D≈10 km. Comets typically impact at higher speeds, reducing the impactor mass for the same impact energy. Although it is increasingly recognized that a continuum exists between comets and asteroids, ’comets’ are considered to be more ice-rich (estimates for67P/ Churyumov-Gerasimenko are about 20%;), implying lower Ir contents per impactor mass. A carbonaceous chondrite-like asteroid of the appropriate size would likely deliver≈2.3×1011gof Ir, in the center of the estimated mass range of the global Ir layer; but a comet would only deliver∼0.1×1011g, because it would be less massive. Although these conclusions are long standing, Siraj & Loeb have recently argued anew in favor of a comet over an asteroid, based on dynamical and geochemical evidence. Here we demonstrate that their arguments are based on misinterpretations of the literature, and that an asteroid is in fact still highly favored over a comet.

Steven J Desch↗
NASA NTRS2014

The Interior and Orbital Evolution of Charon as Preserved in Its Geologic Record

Pluto and its largest satellite, Charon, currently orbit in a mutually synchronous state; both bodies continuously show the same face to one another. This orbital configuration is a natural end-state for bodies that have undergone tidal dissipation. In order to achieve this state, both bodies would have experienced tidal heating and stress, with the extent of tidal activity controlled by the orbital evolution of Pluto and Charon and by the interior structure and rheology of each body. As the secondary, Charon would have experienced a larger tidal response than Pluto, which may have manifested as observable tectonism. Unfortunately, there are few constraints on the interiors of Pluto and Charon. In addition, the pathway by which Charon came to occupy its present orbital state is uncertain. If Charon's orbit experienced a high-eccentricity phase, as suggested by some orbital evolution models, tidal effects would have likely been more significant. Therefore, we determine the conditions under which Charon could have experienced tidally-driven geologic activity and the extent to which upcoming New Horizons spacecraft observations could be used to constrain Charon's internal structure and orbital evolution. Using plausible interior structure models that include an ocean layer, we find that tidally-driven tensile fractures would likely have formed on Charon if its eccentricity were on the order of 0.01, especially if Charon were orbiting closer to Pluto than at present. Such fractures could display a variety of azimuths near the equator and near the poles, with the range of azimuths in a given region dependent on longitude; east-west-trending fractures should dominate at mid-latitudes. The fracture patterns we predict indicate that Charon's surface geology could provide constraints on the thickness and viscosity of Charon's ice shell at the time of fracture formation.

Tectonics↗
NASA NTRS2016

Detection and Measurement of Micrometeoroids with LISA Pathfinder

The Solar System contains a population of dust and small particles originating from asteroids, comets, and other bodies. These particles have been studied using a number of techniques ranging from in-situ satellite detectors to analysis of lunar microcraters to ground-based observations of zodiacal light. In this paper, we describe an approach for using the LISA Pathfinder (LPF) mission as an instrument to detect and characterize the dynamics of dust particles in the vicinity of Earth-Sun L1. Launched on Dec. 3rd, 2015, LPF is a dedicated technology demonstrator mission that will validate several key technologies for a future space-based gravitational-wave observatory. The primary science instrument aboard LPF is a precision accelerometer which we show will be capable of sensing discrete momentum impulses as small as 4 x 10(exp -8) N s.We then estimate the rate of such impulses resulting from impacts of micrometeoroids based on standard models of the micrometeoroid environment in the inner solar system. We find that LPF may detect dozens to hundreds of individual events corresponding to impacts of particles with masses greater than 10(exp -9) g during LPF's roughly six-month science operations phase in a 5 x 10(exp 5) km by 8 x 10(exp 5) km Lissajous orbit around L1. In addition, we estimate the ability of LPF to characterize individual impacts by measuring quantities such as total momentum transferred, direction of impact, and location of impact on the spacecraft. Information on flux and direction provided by LPF may provide insight as to the nature and origin of the individual impact and help constrain models of the interplanetary dust complex in general. Additionally, this direct in situ measurement of micrometeoroid impacts will be valuable to designers of future spacecraft targeting the environment around L1.

meteoroids – instrumentation: miscellaneous↗
NASA NTRSDate not supplied

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain-Corrected Satellite Sea Surface Salinity Using the GMOA S2s Forecast System

The ENSO phenomenon has a significant global socio-economic impact and has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help improve the initialization of the thermocline, while satellite SST aids in constraining surface heat-fluxes, leading to improved coupled system sub-seasonal to seasonal forecasts. However, few studies have focused on improving the near-surface density and mixing through satellite sea surface salinity (SSS) assimilation. The few ocean models that assimilate satellite SSS, bias correct to normalize towards the near-surface Argo data for expediency. This assumption is likely inadequate in rainy regions, where buoyant water forms a fresh surface lens. In previous work, we showed that adjusting SSS to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 improves the near-surface density and mixed layer depth, leading to deeper thermocline and improved NINO3.4 SST forecasts. We now utilize the Soil Moisture and Ocean Salinity rain-corrected (SMOS_RC) SSS, available in SMOS-CATDS products, to represent Sb more accurately at the first model layer (e.g., 5 m). Rather than a diffusivity model as RIM, SMOS_RC uses a statistical correction dependent on Integrated Multi-satellitE Retrievals for GPM (IMERG) rain rates, established on observed SMOS SSS decreases related to Sb in the presence of rain (Supply et al., 2020). For all experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is, and a separate reanalysis assimilates SMOS_RC. We assess the impact on near-surface and subsurface dynamics by validating against observations and explore how SSS assimilation (SMOS vs SMOS_RC) impacts ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (S2S-v3, Molod et al. 2020). We show that improved estimates of density and near-surface mixing led to more accurate coupled air/sea interaction and better ENSO forecasts. The increased SSS, resulting from the removal of the instantaneous rain effect, modifies the ocean state by enhancing mixing and deepening the thermocline.

Veronica Ruiz Xomchuk↗
NASA NTRS2008

Neutron Stars and Black Holes Seen with the Rossi X-Ray Timing Explorer (RXTE)

Astrophysical X-rays bring information about location, energy, time, and polarization. X-rays from compact objects were seen in the first explorations to vary in time. Eclipses and pulsations have simple explanations that identified the importance of X-ray binaries and magnetic neutron stars in the first decade of X-ray astronomy. The dynamics of accretion onto stellar and supermassive black holes and onto neutron stars with relatively low magnetic fields shows up as more complex variations, quasi-periodic oscillations, noise with characteristic frequency spectra, broad-band changes in the energy spectra. To study these variations, RXTE instruments needed to have large area and operational flexibility to find transient activity and observe when it was present. Proportional counters and Phoswich scintillators provided it in a modest mission that has made textbook level contributions to understanding of compact objects. The first seen, and the brightest known, X-ray binary, Sco X-1 is one of a class of neutron stars with low mass companions. Before RXTE, none of these had been seen to show pulsations, though they were hypothesized to be the precursors of radio pulsars with millisecond periods and low magnetic fields. RXTE's large area led to identifying coherent millisecond pulsars in a subset which are relatively faint transients. It also led to identifying short episodes of pulsation during thermonuclear bursts, in sources where a steady signal is not seen. The X-ray stage verifies the evolution that produces millisecond radio pulsars.Masses and radii of neutron stars are being determined by various techniques, constraining the equation of state of matter at nuclear densities. Accretion should lead to a range of neutron star masses. An early stage of superstrong magnetic field neutron stars is now known to produce X-ray and gamma-ray bursts in crust quakes and magnetic field reconnection releases of energy. Soft Gamma Repeaters, Anomolous X-ray Pulsars, and high magnetic field rotation-powered pulsars are all now called magnetars, because they have pulse periods indicating they are slowing down as they would with magnetic dipole radiation for a surface field above 5x1013 gauss. The accretion disk has been connected to the launching of radio jets from black holes, and even from neutron stars. Estimates of the angular momenta of black holes are being made from different approaches, modelling a high frequency oscillation that may be related to how close the inner part of the accretion disk is to the black hole, modelling the continua spectra of the X-ray emission, and modeling the emission of red-shifted iron that may be emitted from the accretion disk. These investigations require early discovery of the black hole transient with the All Sky Monitor on RXTE or other monitoring information, frequent extended observations, and coordinated observations with missions that give higher energy resolution, or radio and infrared information.

Swank, Jean↗
NASA NTRS2018

Impact of Assimilating AIRS Cloud-Cleared Radiances on Atmospheric Dynamics and Polar Low Representation at High Latitudes

This study explores the sensitivity of planetary boundary layer height and related atmospheric dynamics to the assimilation of cloud-cleared AIRS (Atmospheric Infrared Sounder) radiances in the Goddard Earth Observing System (GEOS, version 5) data assimilation and forecast system during the boreal fall 2014 season using observing system experiments (OSEs). Examined here are comparisons between the current, operational approach of assimilating AIRS clear-sky radiances against the assimilation of cloud-cleared radiances (CCR). In polar regions, assimilation of AIRS CCRs is particularly beneficial because of the sparsity of conventional observations and the prevalence of extended low-level stratus cloud cover, which limit the ingestion of clear-sky data. Assimilation of hyperspectral infrared information from AIRS over the Arctic region slightly modifies the lower midtropospheric temperature structure, which in turn contributes to adjustments in geopotential height, affecting the baroclinic instability properties over the entire hemisphere and explaining the overall improvement in global forecast skill. In addition, it is shown that the assimilation of CCRs benefits the representation of convectively-driven small-scale cyclones at high latitudes in the same way as previously noted for tropical cyclones. Specifically, assimilation of CCRs create a temperature dipole over the top of meteorologically active and strongly convective systems such as polar, arctic, and antarctic lows, which helps constrain the analyzed representation of their scale and vertical structure.

Mcgrath-Spangler, E. L.↗
NASA NTRS2022

The breakup of a long-period comet is not a likely match to the Chicxulub impactor

Since the discovery of Ir in the clay layer at the K-Pg boundary, scientists have sought to constrain the origin of the extraterrestrial impactor that triggered the end-Cretaceous mass extinction of the non-avian dinosaurs and other species. While the first proposal was for an asteroid, for a while some theories invoked a cometary impactor to explain perceived periodicities in mass extinctions. Such models have long been disfavored by the mass of Ir in the layer, inferred to be 2.0−2.8×10^(11)cu.g. The size of the Chicxulub crater leads to an estimated asteroid impactor diameter, D≈10km. Comets typically impact at higher speeds, reducing the impactor mass for the same impact energy. Although it is increasingly recognized that a continuum exists between comets and asteroids, ’comets’ are considered to be more ice-rich (estimates for 67P/Churyumov-Gerasimenko are about 20%;), implying lower Ir contents per impactor mass. A carbonaceous chondrite-like asteroid of the appropriate size would likely deliver ≈2.3×10^(11)g of Ir, in the center of the estimated mass range of the global Ir layer; but a comet would only deliver ∼0.1×10^(11)g, because it would be less massive. Although these conclusions are long standing, Siraj and Loeb have recently argued anew in favor of a comet over an asteroid, based on dynamical and geochemical evidence. Here we demonstrate that their arguments are based on misinterpretations of the literature, and that an asteroid is in fact still highly favored over a comet.

Steven J. Desch↗
NASA NTRS2026

Physical Parameters of Space Mission Asteroid Targets

Ground-based characterization of asteroids that are planned or potential targets of space missions provides important data on their physical parameters and properties. Knowledge of the properties of mission targets is important especially during mission preparation and planning, as it serves to select best or suitable targets for specific purpose of a given mission, to prepare mission plans, to constrain possible mission scenarios, and to design mission experiments. Such characterization efforts may be particularly critical for flyby missions that take only limited data during the high-speed flybys of their target asteroids, but they also provide very crucial data for targets of rendezvous missions. Moreover, long-term observations taken from Earth also allows the modelling of, or constraining, processes acting on the asteroids over extended time scales. Over the past years and decades we have obtained rich data on physical parameters of 82 asteroids that are planned or potential targets of space missions. Our primary observing technique is time-resolved (lightcurve) photometry, but we use also data obtained with other techniques, such as spectroscopy, thermal or radar observations. 15 of the 82 characterized asteroids are planned or possible targets of several space missions that are in flight or preparation, such as ESA’s Hera, RAMSES and PRIAMOS, NASA’s OSIRIS-APEX, JAXA’s Hayabusa2#, DESTINY+ and Next Generation Sample Return (NGSR), the Emirates Mission to Asteroids (EMA), and Karman+’s High Frontier, but we have also characterized 67 asteroids that are potential targets of space missions for their low delta-V’s and were announced as being “of interest to NASA” in the Small-Bodies-Observations- NASA mailing list or listed on the NHATS page of “Accessible NEAs”. The sample of asteroid targets span 3 orders of magnitude in size, with absolute magnitudes H from 12.57 to 26.8, corresponding to diameters from about 10 m to about 10 km. The sample contains a variety of taxonomic types and physically or dynamically interesting objects. Among them, we have identified 5 binary asteroids and 20 tumblers (i.e., asteroids in excited, non-principal axis rotation states). Rotation periods of the characterized asteroids range from 1.45 min to 280 h, reflecting diversity of their properties and formation/evolution paths. We will present an overview of the data set and highlight several representative cases.

Petr Pravec↗
NASA NTRS1996

Quality and Control of Water Vapor Winds

Water vapor imagery from the geostationary satellites such as GOES, Meteosat, and GMS provides synoptic views of dynamical events on a continual basis. Because the imagery represents a non-linear combination of mid- and upper-tropospheric thermodynamic parameters (three-dimensional variations in temperature and humidity), video loops of these image products provide enlightening views of regional flow fields, the movement of tropical and extratropical storm systems, the transfer of moisture between hemispheres and from the tropics to the mid- latitudes, and the dominance of high pressure systems over particular regions of the Earth. Despite the obvious larger scale features, the water vapor imagery contains significant image variability down to the single 8 km GOES pixel. These features can be quantitatively identified and tracked from one time to the next using various image processing techniques. Merrill et al. (1991), Hayden and Schmidt (1992), and Laurent (1993) have documented the operational procedures and capabilities of NOAA and ESOC to produce cloud and water vapor winds. These techniques employ standard correlation and template matching approaches to wind tracking and use qualitative and quantitative procedures to eliminate bad wind vectors from the wind data set. Techniques have also been developed to improve the quality of the operational winds though robust editing procedures (Hayden and Veldon 1991). These quality and control approaches have limitations, are often subjective, and constrain wind variability to be consistent with model derived wind fields. This paper describes research focused on the refinement of objective quality and control parameters for water vapor wind vector data sets. New quality and control measures are developed and employed to provide a more robust wind data set for climate analysis, data assimilation studies, as well as operational weather forecasting. The parameters are applicable to cloud-tracked winds as well with minor modifications. The improvement in winds through use of these new quality and control parameters is measured without the use of rawinsonde or modeled wind field data and compared with other approaches.

Jedlovec, Gary J.↗
NASA NTRS2016

Virtual Flight Demonstration of the Stratospheric Dual-Aircraft Platform

A baseline configuration for the dual-aircraft platform (DAP) concept is described and evaluated in a physics-based flight dynamics simulations for two month-long missions as a communications relay in the lower stratosphere above central Florida, within 150-miles of downtown Orlando.The DAP configuration features two large glider-like (130 ft wing span) unmanned aerial vehicles connected via a long adjustable cable (total extendible length of 3000 ft) which effectively sail without propulsion using available wind shear. Use of onboard LiDAR wind profilers to forecast wind distributions are found to be necessary to enable the platform to efficiently adjust flight conditions to remain sailing by finding sufficient wind shear across the platform. The aircraft derive power from solar cells, like a conventional solar aircraft, but also extract wind power using the propeller as a turbine when there is an excess of wind shear available.Month-long atmospheric profiles (at 3-5 min intervals) in the vicinity of 60,000-ft are derived from archived data measured by the 50-Mhz Doppler Radar Wind Profiler at Cape Canaveral and used in the DAP flight simulations. A cursory evaluation of these datasets show that sufficient wind shear for DAP sailing is persistent, suggesting that DAP could potentially sail over 90% of the month-long durations even when limited by modest ascent/descent rates.DAP's novel guidance software uses a non-linear constrained optimization technique to define waypoints such that sailing mode of flight is maintained where possible, and minimal thrust is required where sailing is not practical. A set of constraints are identified which result in waypoints that enable efficient flight (i.e., minimal use of propulsion) over the two month-long flight simulations. Waypoint solutions may need to be tabulated for a wide range of potential atmospheric conditions and stored onboard for quick retrieval on a real DAP.DAP's flight control software uses an unconventional mixture of spacecraft and aircraft control techniques. Flight simulations confirms that this controls approach enables the platform to consistently reach successive waypoints over the month-long flight simulations.The ability of DAP to transition between the sailing mode (i.e., cable tension is high) and standard formation flight (i.e., cable tension is low) is a vital capability (e.g., to enable intermittent turns while stationkeeping). A new method to perform these transitions has been identified and characterized with flight simulation which requires special aircraft modifications.The energy-usage of the DAP configuration during two month-long stationkeeping missions over central Florida (i.e., stationkeeping over Orlando) is evaluated and compared to that of a pure solar aircraft of the same weight and aerodynamic performance. DAP is shown to consistently reduce net propulsion usage while simultaneously increasing solar energy capture.A baseline 700 GHz communications system is described and its performance evaluated for the proposed mission over central Florida. It is found that the variable roll orientation of the aircraft would increase the power required to maintain coverage over the stationkeeping radius of 150 miles (e.g., by as much as 100% when DAP is 150 miles from Orlando), compared to level flight. This effect can be mitigated via additional antenna design complexity or a more restricted stationkeeping radius.

Demonstrations↗
NASA NTRS2022

The Warm Neptune GJ 3470b Has a Polar Orbit

The warm Neptune GJ 3470b transits a nearby (d = 29 pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5 m Telescope at Kitt Peak Observatory, we model the classical Rossiter–McLaughlin effect, yielding a sky-projected obliquity of λ=98 +15 −12 ˚and a v sin i = 0.85 +0.27 −0.33 kms -1 . Leveraging information about the rotation period and size of the host star, our analysis yields a true obliquity of ψ=95 +9 −8 ◦ , revealing that GJ 3470b is on a polar orbit. Using radial velocities from HIRES, HARPS, and the Habitable-zone Planet Finder, we show that the data are compatible with a long-term radial velocity (RV) slope of 𝛾̀=-0.0022±0.0011 ms -1 day -1 over a baseline of 12.9 yr. If the RV slope is due to acceleration from another companion in the system, we show that such a companion is capable of explaining the polar and mildly eccentric orbit of GJ 3470b using two different secular excitation models. The existence of an outer companion can be further constrained with additional RV observations, Gaia astrometry, and future high-contrast imaging observations. Lastly, we show that tidal heating from GJ 3470b’s mild eccentricity has most likely inflated the radius of GJ 3470b by a factor of ∼1.5–1.7, which could help account for its evaporating atmosphere.

Exoplanet astronomy↗
NASA NTRS2014

Humic Acid Complexation of Th, Hf and Zr in Ligand Competition Experiments: Metal Loading and Ph Effects

The mobility of metals in soils and subsurface aquifers is strongly affected by sorption and complexation with dissolved organic matter, oxyhydroxides, clay minerals, and inorganic ligands. Humic substances (HS) are organic macromolecules with functional groups that have a strong affinity for binding metals, such as actinides. Thorium, often studied as an analog for tetravalent actinides, has also been shown to strongly associate with dissolved and colloidal HS in natural waters. The effects of HS on the mobilization dynamics of actinides are of particular interest in risk assessment of nuclear waste repositories. Here, we present conditional equilibrium binding constants (Kc, MHA) of thorium, hafnium, and zirconium-humic acid complexes from ligand competition experiments using capillary electrophoresis coupled with ICP-MS (CE- ICP-MS). Equilibrium dialysis ligand exchange (EDLE) experiments using size exclusion via a 1000 Damembrane were also performed to validate the CE-ICP-MS analysis. Experiments were performed at pH 3.5-7 with solutions containing one tetravalent metal (Th, Hf, or Zr), Elliot soil humic acid (EHA) or Pahokee peat humic acid (PHA), and EDTA. CE-ICP-MS and EDLE experiments yielded nearly identical binding constants for the metal- humic acid complexes, indicating that both methods are appropriate for examining metal speciation at conditions lower than neutral pH. We find that tetravalent metals form strong complexes with humic acids, with Kc, MHA several orders of magnitude above REE-humic complexes. Experiments were conducted at a range of dissolved HA concentrations to examine the effect of [HA]/[Th] molar ratio on Kc, MHA. At low metal loading conditions (i.e. elevated [HA]/[Th] ratios) the ThHA binding constant reached values that were not affected by the relative abundance of humic acid and thorium. The importance of [HA]/[Th] molar ratios on constraining the equilibrium of MHA complexation is apparent when our estimated Kc, MHA values attained at very low metal loading conditions are compared to existing literature data. Overall, experimental data suggest that the tetravalent transition metal/-actinide-humic acid complexation is important over a wide range of pH values, including mildly acidic conditions, and thus, these complexes should be included in speciation models.

humic acid↗
NASA NTRS2018

GEOS S2S-2_1 File Specification: GMAO Seasonal and Sub-Seasonal Forecast Output

The NASA GMAO seasonal (9 months) and subseasonal (45 days) forecasts are produced with the Goddard Earth Observing System (GEOS) Atmosphere-Ocean General Circulation Model and Data Assimilation System Version S2S-2_1. The new system replaces version S2S-1.0 described in Borovikov et al (2017), and includes upgrades to many components of the system. The atmospheric model includes an upgrade from a pre-MERRA-2 version running on a latitude-longitude grid at approx. 1 degree resolution to a current version running on a cubed sphere grid at approximately 1/2 degree resolution. The important developments are related to the dynamical core (Putman et al., 2011), the moist physics (''two-moment microphysics'' of Barahona et al., 2014) and the cryosphere (Cullather et al., 2014). As in the previous GMAO S2S system, the land model is that of Koster et al (2000). GMAO S2S-2_1 now includes the Goddard Chemistry Aerosol Radiation and Transport (GOCART, Colarco et al., 2010) single moment interactive aerosol model that includes predictive aerosols including dust, sea salt and several species of carbon and sulfate. The previous version of GMAO S2S specified aerosol amounts from climatology, which were used to inform the atmospheric radiation only. The ocean model includes an upgrade from MOM4 to MOM5 (Griffies 2012), and continues to be run on the tripolar grid at approximately 1/2 degree resolution in the tropics with 40 vertical levels. As in S2S-1.0, the sea ice model is from the Los Alamos Sea Ice model (CICE4, Hunke and Lipscomb 2010). The Ocean Data Assimilation System (ODAS) has been upgraded from the one described in Borovikov et al., 2017 to one that uses a modified version of the Penny, 2014 Local Ensemble Transform Kalman Filter (LETKF), and now assimilates along-track altimetry. The ODAS also does a nudging to MERRA-2 SST and sea ice boundary conditions. The atmospheric data assimilation fields used to constrain the atmosphere in the ODAS have been upgraded from MERRA to a MERRA-2 like system. The system is initialized using a MERRA-2-like atmospheric reanalysis (Gelaro et al. 2017) and the GMAO S2S-2_1 ocean analysis. Additional ensemble members for forecasts are produced with initial states at 5-day intervals, with additional members based on perturbations of the atmospheric and ocean states. Both subseasonal and seasonal forecasts are submitted to the National MultiModel Ensemble (NMME) project, and are part of the US/Canada multimodel seasonal forecasts (http://www.cpc.ncep.noaa.gov/products/NMME/). A large suite of retrospective forecasts (''hindcasts'') have been completed, and contribute to the calculation of the model's baseline climatology and drift, anomalies from which are the basis of the seasonal forecasts.

Seasonal Forecast↗
NASA NTRS2017

OpenSatKit Enables Quick Startup for CubeSat Missions

The software required to develop, integrate, and operate a spacecraft is substantial regardless of whether its a large or small satellite. Even getting started can be a monumental task. To solve this problem, NASAs Core Flight System (cFS), NASA's 42 spacecraft dynamics simulator, and Ball Aerospaces COSMOS ground system have been integrated together into a kit called OpenSatKit that provides a complete and open source software solution for starting a new satellite mission. Users can have a working system with flight software, dynamics simulation, and a ground command and control system up and running within hours.Every satellite mission requires three primary categories of software to function. The first is Flight Software (FSW) which provides the onboard control of the satellites and its payload(s). NASA's cFS provides a great platform for developing this software. Second, while developing a satellite on earth, it is necessary to simulate the satellites orbit, attitude, and actuators, to ensure that the systems that control these aspects will work correctly in the real environment. NASAs 42 simulator provides these functionalities. Finally, the ground has to be able to communicate with the satellite, monitor its performance and health, and display its data. Additionally, test scripts have to be written to verify the system on the ground. Ball Aerospace's COSMOS command and control system provides this functionality. Once the OpenSatKit is up and running, the next step is to customize the platform and get it running on the end target. Starting from a fully working system makes porting the cFS from Linux to a users platform much easier. An example Raspberry Pi target is included in the kit so users can gain experience working with a low cost hardware target. All users can benefit from OpenSatKit but the greatest impact and benefits will be to SmallSat missions with constrained budgets and small software teams. This paper describes OpenSatKits system design, the steps necessary to run the system to target the Raspberry Pi, and future plans. OpenSatKit is a free fully functional spacecraft software system that we hope will greatly benefit the SmallSat community.

Software↗
NASA NTRS2013

Comet C2012 S1 (ISON): Observations of the Dust Grains From SOFIA and of the Atomic Gas From NSO Dunn and Mcmath-Pierce Solar Telescopes

Comet C/2012 S1 (ISON) is unique in that it is a dynamically new comet derived from the Oort cloud reservoir of comets with a sun-grazing orbit. Infrared (IR) and visible wavelength observing campaigns were planned on NASA's Stratospheric Observatory For Infrared Astronomy (SOFIA) and on National Solar Observatory Dunn (DST) and McMath-Pierce Solar Telescopes, respectively. We highlight our SOFIA (+FORCAST) mid- to far-IR images and spectroscopy (approx. 5-35 microns) of the dust in the coma of ISON are to be obtained by the ISON-SOFIA Team during a flight window 2013 Oct 21-23 UT (r_h approx. = 1.18 AU). Dust characteristics, identified through the 10 micron silicate emission feature and its strength, as well as spectral features from cometary crystalline silicates (Forsterite) at 11.05-11.2 microns, and near 16, 19, 23.5, 27.5, and 33 microns are compared with other Oort cloud comets that span the range of small and/or highly porous grains (e.g., C/1995 O1 (Hale-Bopp) and C/2001 Q4 (NEAT) to large and/or compact grains (e.g., C/2007 N4 (Lulin) and C/2006 P1 (McNaught)). Measurement of the crystalline peaks in contrast to the broad 10 and 20 micron amorphous silicate features yields the cometary silicate crystalline mass fraction, which is a benchmark for radial transport in our protoplanetary disk. The central wavelength positions, relative intensities, and feature asymmetries for the crystalline peaks may constrain the shapes of the crystals. Only SOFIA can look for cometary organics in the 5-8 micron region. Spatially resolved measurements of atoms and simple molecules from when comet ISON is near the Sun (r_h< 0.4 AU, near Nov-20-Dec-03 UT) were proposed for by the ISON-DST Team. Comet ISON is the first comet since comet Ikeya-Seki (1965f) suitable for studying the alkalai metals Na and K and the atoms specifically attributed to dust grains including Mg, Si, Fe, as well as Ca. DST's Horizontal Grating Spectrometer (HGS) measures 4 settings: Na I, K, C2 to sample cometary organics (along with Mg I), and [OI] as a proxy for activity from water (along with Si I and Fe I). State-of-the-art instruments that will also be employed include IBIS, which is a Fabry-Perot spectral imaging system that concurrently measures lines of Na, K, Ca II, or Fe, and ROSA (CSUN/QUB), which is a rapid imager that simultaneously monitors Ca II or CN. From McMath-Pierce, the Solar-Stellar Spectrograph also will target ISON (320-900 nm, R approx. 21,000, r_h<0.3 AU). Assuming survival, the intent is to target ISON over r_h<0.4 AU, characteristic of prior Na detections.

Oort cloud comets↗