Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust estimation”

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 901 records · Page 50

Memory-Efficient Onboard Rock Segmentation

Rockster-MER is an autonomous perception capability that was uploaded to the Mars Exploration Rover Opportunity in December 2009. This software provides the vision front end for a larger software system known as AEGIS (Autonomous Exploration for Gathering Increased Science), which was recently named 2011 NASA Software of the Year. As the first step in AEGIS, Rockster-MER analyzes an image captured by the rover, and detects and automatically identifies the boundary contours of rocks and regions of outcrop present in the scene. This initial segmentation step reduces the data volume from millions of pixels into hundreds (or fewer) of rock contours. Subsequent stages of AEGIS then prioritize the best rocks according to scientist- defined preferences and take high-resolution, follow-up observations. Rockster-MER has performed robustly from the outset on the Mars surface under challenging conditions. Rockster-MER is a specially adapted, embedded version of the original Rockster algorithm ("Rock Segmentation Through Edge Regrouping," (NPO- 44417) Software Tech Briefs, September 2008, p. 25). Although the new version performs the same basic task as the original code, the software has been (1) significantly upgraded to overcome the severe onboard re source limitations (CPU, memory, power, time) and (2) "bulletproofed" through code reviews and extensive testing and profiling to avoid the occurrence of faults. Because of the limited computational power of the RAD6000 flight processor on Opportunity (roughly two orders of magnitude slower than a modern workstation), the algorithm was heavily tuned to improve its speed. Several functional elements of the original algorithm were removed as a result of an extensive cost/benefit analysis conducted on a large set of archived rover images. The algorithm was also required to operate below a stringent 4MB high-water memory ceiling; hence, numerous tricks and strategies were introduced to reduce the memory footprint. Local filtering operations were re-coded to operate on horizontal data stripes across the image. Data types were reduced to smaller sizes where possible. Binary- valued intermediate results were squeezed into a more compact, one-bit-per-pixel representation through bit packing and bit manipulation macros. An estimated 16-fold reduction in memory footprint relative to the original Rockster algorithm was achieved. The resulting memory footprint is less than four times the base image size. Also, memory allocation calls were modified to draw from a static pool and consolidated to reduce memory management overhead and fragmentation. Rockster-MER has now been run onboard Opportunity numerous times as part of AEGIS with exceptional performance. Sample results are available on the AEGIS website at http://aegis.jpl.nasa.gov.

Burl, Michael C.↗

In-Lab Rapid Analytical Detection of Lunar Volatiles By Universal Gas Analyzer With Comparison to GC-MS System

Introduction: The curation of permanently shadowed regions (PSRs) [1] on the lunar surface is centered around studies based upon the observed volatiles from the LCROSS mission [2]. The rapid detection of important volatile gases and vapors present in planetary bodies and Astromaterials by a standalone analytical device is an area of intense research interest in our group and Planetary Exploration & Astromaterials Research Laboratory (PEARL) facility and this work is relevant to the future preparation of viable lunar simulants for testing curation efforts down the road. The groundbreaking results obtained from the LCROSS Mission [2] open the requirements for the direct detection of volatiles present in regolith materials collected from the lunar surface. The mass spectrometry of volatile chemicals is a general technique that utilizes a set of instruments that creates charged ions from a gaseous chemical species and measures the intensities vs. mass-to-charge ratio (m/z) [3]. In this context, we discuss in-lab experimental results and procedures for rapid qualitative analysis of main LCROSS volatiles (water, H2S, NH3, CO2, and CH3OH) by a Universal Gas Analyzer (UGA) instrument. Additionally, the instrument performance was evaluated by measuring the isotopic abundance ratio of atmospheric Ar-40 to Ar-36 present in room air since, argon is a relevant gas in planetary studies as it can provide an insight and reference point to isotope studies [4]. Additional, cross comparisons were attempted and made between the two instruments to develop a robust analytical technique by comparing mass spectral data for H2S headspace samples with a Trace-1310/ ISQ 7000 (ThermoFisher Scientific.) GC-MS system. Background: The benchtop UGA System is equipped with an SRS UGA 300 quadrupole mass spectrometer designed and built by Stanford Research Systems [5]. This system can be configured for several types of gaseous chemical analysis. The inlet line continuously samples gases at low flow rates (several milliliters per minute) through a capillary limiting the intake pressure making the instrument ideal for online analysis of select gases and/or room atmosphere. Moreover, in our current UGA system, a change in composition at the inlet can be detected in about 200 milliseconds and a complete spectrum is acquired (for a range of 1-100 amu) in under 45 sec with masses measured at rates up to 25 msec per point [5]. This system provides a quick upstream analytical data that we can then compare to results obtained by our GC-MS system. Sample Preparation: Small volume (2-4 mL) of analyte sample was taken in a 10 mL glass vial and sealed with a crimped cap and purged with pure Ar or N2 gas to displace air from the top. The headspace sample was scanned by the UGA instrument at analog, histogram, and pressure vs. time modes. The isotopic abundance ratio for 40Ar-to-36Ar was estimated by measuring partial pressure vs. time scans and setting the mass at 40 and 36 respectively. Results and Discussions: In this work, we have investigated the applicability of the UGA system by qualitative analysis of a series of LCROSS volatiles measured individually. Fig. 1 demonstrates a set of vertically offset spectra for the partial pressures measured as a function of mass-to-charge (m/z) ratios. The average acquisition time for each spectrum was less than a minute suggesting that the UGA system is ideal for quick analysis of geochemical volatiles. For the cross-comparison, we analyzed an H2S headspace sample by a Trace-1310/ISQ-7000 system and compared mass spectral data with previously measured UGA histogram scan data (Fig. 2). In both cases, major peak positions are the same, however, the intensities of fragment ions ([1H132S]+ and [32S]+) are higher for UGA suggesting that the fragment ionization process is stronger in UGA compared to that of GC-MS. To investigate how the integrated area under each chromatogram varies with the headspace sample volume, a set of five H2S headspace samples with increasing volumes was analyzed by the GC-MS system (Fig. 3, inset). A small volume (e.g., 200 to 1000 µL) of H2S/H2O vapor was withdrawn from a 20 mL stock sample vial containing ~5 mL of 0.4% H2S in water by a gas-tight syringe and added to another 20 mL vial filled with argon and analyzed by the GC-MS system. Finally, the UGA detector sensitivity was evaluated by calculating the atmospheric 40Ar-to-36Ar isotopic abundance ratio in room air by running a partial pressure vs. time scan with setting the atomic mass at 40, and 36. Fig. 4(a) shows a ~25 min duration “P vs. time” scan for 40Ar (plot for 36Ar is not shown). The partial pressure values (~100 points) were corrected by subtracting the corresponding background pressure value for 37Ar and utilized to calculate 40Ar-to-36Ar isotopic abundance ratios as shown by Fig. 4b. The average isotopic abundance ratio is ~306 with a 2*STDEV ~13. This abundance ratio is significantly close to the previously reported value of 298.56 [6] and the ratio obtained by our GC-MS system (303 for a UHP grade Ar sample). Conclusions: Our study strongly evidenced that the benchtop UGA system is a valuable analytical tool for the detection of major LCROSS volatiles. The rapid scanning capability, the inexpensiveness of the whole system, and impressive detection sensitivity prove its worthiness as an essential device for advanced geochemical applications. Moreover, cross comparisons with the GC-MS provide important bridges into advanced curatorial efforts into the future. References: [1] Bickel, V.T., et al. (2021) Nat Commun 12, 5607. [2] Colaprete, A., et al. (2010) Science, 330, 463-468. [3] Glavin, D. P. et al. (2012) 2012 IEEE Aerospace Conference, 1-11. [4] Willett, C. D., et al. (2022) Geochimica et Cosmochimica Acta 329, 119-134. [5] Operation Manual and Programming Reference. (2018) Universal gas Analyzers, Stanford Research Systems. [6] Lee, J. Y., et al. (2006) Geochimica et Cosmochimica Acta 70, 4507–4512. Notes: (4 figures are attached with text as shown by the attached file)

Curation↗

A Chlorophyll-a Algorithm for Landsat-8 Based on Mixture Density Networks

Retrieval of aquatic biogeochemical variables, such as the near-surface concentration of chlorophyll-a (Chla) in inland and coastal waters via remote observations, has long been regarded as a challenging task. This manuscript applies Mixture Density Networks (MDN) that use the visible spectral bands available by the Operational Land Imager (OLI) aboard Landsat-8 to estimate Chla. We utilize a database of co-located in situ radiometric and Chla measurements (N = 4,354), referred to as Type A data, to train and test an MDN model (MDN(A)). This algorithm’s performance, having been proven for other satellite missions, is further evaluated against other widely used machine learning models (e.g., support vector machines), as well as other domain-specific solutions (OC3), and shown to offer significant advancements in the field. Our performance assessment using a held-out test data set suggests that a 49% (median) accuracy with near-zero bias can be achieved via the MDN(A) model, offering improvements of 20 to 100% in retrievals with respect to other models. The sensitivity of the MDN(A) model and benchmarking methods to uncertainties from atmospheric correction (AC) methods, is further quantified through a semi-global matchup dataset (N = 3,337), referred to as Type B data. To tackle the increased uncertainties, alternative MDN models (MDN(B)) are developed through various features of the Type B data (e.g., Rayleigh-corrected reflectance spectra ρ(s)). Using held-out data, along with spatial and temporal analyses, we demonstrate that these alternative models show promise in enhancing the retrieval accuracy adversely influenced by the AC process. Results lend support for the adoption of MDN(B) models for regional and potentially global processing of OLI imagery, until a more robust AC method is developed. Index Terms—Chlorophyll-a, coastal water, inland water, Landsat-8, machine learning, ocean color, aquatic remote sensing.

Brandon Smith↗

An Overview of Experiments and Modeling of Polysiloxane-Coated Thermal Protection Systems for Missions to Mars, Titan, and Beyond.

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the historic Stardust mission, where it successfully returned samples from a comet’s tail and has since been instrumental in delivering payloads to the surface of Mars [1-3]. Most recently, PICA enabled the safe return of samples collected from asteroid Bennu as part of the OSIRIS-REx mission. This rich legacy underscores PICA’s critical role in allowing NASA’s most ambitious exploration missions. However, the friable nature of its phenolic phase presents challenges during handling and pre-launch activities. To mitigate this issue, PICA is coated with a polysiloxane resin system, which serves to suppress particulate dispersion and thereby safeguard spacecraft components. A comprehensive understanding of the polysiloxane resin’s behavior is imperative, as it profoundly shapes the material response of PICA during atmospheric entry by influencing its thermal and oxidative stability. This influence extends to thermocouple plugs embedded within thermal protection systems. These plugs have demonstrated their significance in missions such as Mars Science Laboratory (MSL) and Mars 2020, where the MEDLI and MEDLI2 instrumentation suites delivered in-valuable insights into the performance of thermal protection systems during entry into the Martian atmosphere [4]. Looking ahead, missions such as Dragonfly, set to descend into Titan’s dense atmosphere, aim to leverage advanced sensor technologies to further refine our understanding of thermal protection response [5]. Moreover, thermocouple plugs play an essential role in validating cutting-edge material response models, such as those pioneered under NASA’s Entry Systems Modeling Project (ESM), designed, in-part, to predict the operational integrity of thermal protection systems under the extreme stresses of atmospheric entry. To achieve these modeling goals, ground-based experiments are crucial to provide the foundational data necessary for developing and refining these predictive tools. To this end, an extensive test campaign was conducted at the Hypersonic Materials Environmental Test System (HyMETS) to investigate the high-temperature behavior of the polysiloxane resin in an air environment [6]. These experiments revealed critical phenomena, including the formation of a silicon oxycarbide layer that enhances oxidation resistance, moderates surface temperatures, and alters in-depth thermal response. Building on these findings, subsequent tests were designed to simulate atmospheric entry conditions in reactive gases, such as CO2 and N2, to mimic the environments of Mars and Titan, respectively, as well as non-reactive gases representing the atmospheres of the Ice Giants (Neptune and Uranus). A heating rate dependent decomposition mechanism has been identified for the polysiloxane resin under oxidizing conditions (Fig. 1). In the initial stage, the resin and the underlying thermal protection system undergo pyrolysis, rapidly generating a thin amorphous silicon oxycarbide interwoven with carbonaceous char and residual fibers from PICA. During the second stage, the nascent oxide layer establishes a robust, oxidation-resistant thermal barrier coating, which significantly impedes heat transfer to the underlying carbonaceous char, resulting in a stagnation of the surface temperature. A key factor contributing to this thermal resistance is the low recombination efficiency of atomic oxygen (γ), which further diminishes the heat load on the material’s interior layers [7]. Moreover, as the surface temperature stagnates, the silicon oxycarbide phase separates into distinct regions of silica and free graphite. Ultimately, when the heat flux reaches a critical threshold, a third stage is triggered, leading to the breakdown of the coating through carbothermal reduction, exposing the underlying char layer. This exposure leads to a dramatic surface temperature spike, driven by highly exothermic reactions between atomic oxygen and the char layer, further accelerating material degradation. A detailed mass and heat transfer model of PICA coated with polysiloxane resin was implemented in the Porous material Analysis Toolbox based on OpenFOAM, PATO [8]. The initial stage was considered negligible in this model because the resin decomposition occurs rapidly within a thin surface layer. Instead, the coating was directly considered as an oxygen-resistant thermal barrier coating. For the second stage, the thin amorphous silicon oxycarbide was treated as a pure silica surface to simplify the thermochemical behavior. The model ac-counts for surface equilibrium processes using representative elements of the coating-environment system. For the third stage, specific boundary conditions were developed to estimate the onset and progression of the coating removal. Two-dimensional material response simulations were conducted to compare uncoated and coated PICA using boundary conditions calibrated with HyMETS data. Fig. 2 illustrates that the simulations closely align with experimental data, successfully reproducing measured temperature profiles. This work will include the latest advancements in the coating model, including the calibration of recombination of atomic oxygen at the surface during the second phase. These simulated results will be further validated against additional CO2 data points from HyMETS, reinforcing the models’ predictive capabilities. These mechanisms and their effects on thermal protection systems, including thermochemical behavior and thermocouple probe performance in extreme environments, provide crucial insights for optimizing spacecraft designs that safeguard scientific payload and ensure mission success in future planetary exploration endeavors.

Active Oxidation↗

Plans for NASA Contributions to Position, Navigation, and Timing at the Moon

Dating back to the Apollo missions, navigation in cislunar space has traditionally been performed with radiometric tracking and orbit estimation in Earth-based systems. The selected landing zones for Apollo enabled direct to Earth line of sight. Up to the present day, lunar orbiting science missions continue to rely on Earth sensors and processing. The upcoming Artemis and science missions seek different objectives, including landing and traversing in lunar locations occulted from view by Earth. Sustaining a burgeoning and diverse lunar mission set requires in-situ navigation capabilities that do not rely exclusively on Earth-based tracking and provide accurate and timely position, navigation, and timing (PNT) knowledge defined in a lunar reference system. These missions prefer a seamless transition from Earth-centric to lunar-centric PNT systems to reduce complexity and size, weight, and power (SWaP) on their missions. NASA’s heritage in successful high-altitude Global Navigation Satellite System (GNSS) reception is informing an expanded Space Service Volume to the moon and the development of lower SWaP GNSS receivers for the lunar missions. To overcome the geometric restrictions from GNSS –from both limited Dilution of Precision at the moon and lunar occultation – NASA is pursuing PNT services from the Lunar Communications Relay and Navigation System (LCRNS) as part of the interoperable LunaNet architecture to provide lunar-centric navigation services. For situations that require additional sensors, e.g. hazard detection, NASA’s optical navigation algorithms and close-range topographic sensors combine with observables from LunaNet in autonomous navigation software that brings robustness to the PNT solution. This presentation will provide an overview of NASA’s LCRNS PNT services, high-altitude GNSS and close-range sensors, processing algorithms and performance analysis, as well as efforts underway to define the underpinning Lunar Reference System.

Cheryl J. Gramling↗

Simulated Impacts of Climate Change on Water Use and Yield of Irrigated Sugarcane in South Africa

Reliable predictions of climate change impacts on water use, irrigation requirements and yields of irrigated sugarcane in South Africa (a water-scarce country) are necessary to plan adaptation strategies. Although previous work has been done in this regard, methodologies and results vary considerably. The objectives were (1) to estimate likely impacts of climate change on sugarcane yields, water use and irrigation demand at three irrigated sugarcane production sites in South Africa (Malelane, Pongola and La Mercy) for current (1980-2010) and future (2070-2100) climate scenarios, using an approach based on the Agricultural Model Inter-comparison and Improvement Project (AgMIP) protocols; and (2) to assess the suitability of this methodology for investigating climate change impacts on sugarcane production. Future climate datasets were generated using the Delta downscaling method and three Global Circulation Models (GCMs) assuming atmospheric CO2 concentration [CO2] of 734 ppm(A2 emissions scenario). Yield and water use were simulated using the DSSAT-Canegro v4.5 model. Irrigated cane yields are expected to increase at all three sites (between 11 and 14%), primarily due to increased interception of radiation as a result of accelerated canopy development. Evapotranspiration and irrigation requirements increased by 11% due to increased canopy cover and evaporative demand. Sucrose yields are expected to decline because of increased consumption of photo-assimilate for structural growth and maintenance respiration. Crop responses in canopy development and yield formation differed markedly between the crop cycles investigated. Possible agronomic implications of these results include reduced weed control costs due to shortened periods of partial canopy, a need for improved efficiency of irrigation to counter increased demands, and adjustments to ripening and harvest practices to counter decreased cane quality and optimize productivity. Although the Delta climate data downscaling method is considered robust, accurate and easily-understood, it does not change the future number of rain-days per month. The impacts of this and other climate data simplifications ought to be explored in future work. Shortcomings of the DSSAT-Canegro model include the simulated responses of phenological development, photosynthesis and respiration processes to high temperatures, and the disconnect between simulated biomass accumulation and expansive growth. Proposed methodology refinements should improve the reliability of predicted climate change impacts on sugarcane yield.

Climate change↗

Science Benefits of Onboard Spacecraft Navigation

Primitive bodies (asteroids and comets), which have remained relatively unaltered since their formation, are important targets for scientific missions that seek to understand the evolution of the solar system. Often the first step is to fly by these bodies with robotic spacecraft. The key to maximizing data returns from these flybys is to determine the spacecraft trajectory relative to the target body-in short, navigate the spacecraft- with sufficient accuracy so that the target is guaranteed to be in the instruments' field of view. The most powerful navigation data in these scenarios are images taken by the spacecraft of the target against a known star field (onboard astrometry). Traditionally, the relative trajectory of the spacecraft must be estimated hours to days in advance using images collected by the spacecraft. This is because of (1)!the long round-trip light times between the spacecraft and the Earth and (2)!the time needed to downlink and process navigation data on the ground, make decisions based on the result, and build and uplink instrument pointing sequences from the results. The light time and processing time compromise navigation accuracy considerably, because there is not enough time to use more accurate data collected closer to the target-such data are more accurate because the angular capability of the onboard astrometry is essentially constant as the distance to the target decreases, resulting in better "plane-of- sky" knowledge of the target. Excellent examples of these timing limitations are high-speed comet encounters. Comets are difficult to observe up close; their orbits often limit scientists to brief, rapid flybys, and their coma further restricts viewers from seeing the nucleus in any detail, unless they can view the nucleus at close range. Comet nuclei details are typically discernable for much shorter durations than the roundtrip light time to Earth, so robotic spacecraft must be able to perform onboard navigation. This onboard navigation can be accomplished through a self- contained system that by eliminating light time restrictions dramatically improves the relative trajectory knowledge and control and subsequently increases the amount of quality data collected. Flybys are one-time events, so the system's underlying algorithms and software must be extremely robust. The autonomous software must also be able to cope with the unknown size, shape, and orientation of the previously unseen comet nucleus. Furthermore, algorithms must be reliable in the presence of imperfections and/or damage to onboard cameras accrued after many years of deep-space operations. The AutoNav operational flight software packages, developed by scientists at the Jet Propulsion Laboratory (JPL) under contract with NASA, meet all these requirements. They have been directly responsible for the successful encounters on all of NASA's close-up comet-imaging missions (see Figure !1). AutoNav is the only system to date that has autonomously tracked comet nuclei during encounters and performed autonomous interplanetary navigation. AutoNav has enabled five cometary flyby missions (Table!1) residing on four NASA spacecraft provided by three different spacecraft builders. Using this software, missions were able to process a combined total of nearly 1000 images previously unseen by humans. By eliminating the need to navigate spacecraft from Earth, the accuracy gained by AutoNav during flybys compared to ground-based navigation is about 1!order of magnitude in targeting and 2!orders of magnitude in time of flight. These benefits ensure that pointing errors do not compromise data gathered during flybys. In addition, these benefits can be applied to flybys of other solar system objects, flybys at much slower relative velocities, mosaic imaging campaigns, and other proximity activities (e.g., orbiting, hovering, and descent/ascent).

Autonomy↗

Spectral Range Within Global aCDOM(440) Algorithms for Oceanic, Coastal, and Inland Waters with Application to Airborne Measurements

The optically active component of dissolved organic material in aquatic ecosystems, or colored dissolved organic matter (CDOM), is represented by the coefficient of absorption due to the dissolved aquatic constituents at 440 nm, aCDOM(440). Remote sensing of aCDOM(440) enables characterization of ecosystem processes and aids in retrieval of chlorophyll a, a proxy for phytoplankton biomass. Spectrally adjacent band-ratio domains, e.g., blue to green, have previously been applied for remote sensing of aCDOM(440) in coastal and oceanic waters with similar results compared to more complex semi-analytical algorithms. Estimation of aCDOM(440) from ratios of the most spectrally separated ocean color wavebands (end members), e.g., ultraviolet (UV) to near-infrared (NIR), termed end-member analysis (EMA), has previously been shown to increase the accuracy of global aCDOM(440) retrievals from in-water observations of diffuse attenuation and to enable a unified algorithmic perspective without requiring regional adjustment of internal bio-optical parameters. EMA of above-water observations is evaluated herein, with a focus on coastal and inland waters in which increasing optical complexity and likelihood of bottom reflectance challenge the oceanic algorithms developed for deep and optically simple (case-1) waters. Analysis herein of three independent, in situ, bio-optical datasets indicates significant correlation between aCDOM(440) and end-member band ratios (next-generation 320 and 780 nm or legacy 412 and 670 nm ratios) with a coefficient of determination, R2, of 0.87 (log-scale) or higher based on a dataset spanning the dynamic range of global, conservative water bodies. For applicable wavelengths, EMA algorithms are shown to agree with case-1 relationships and to produce consistent log-scale uncertainties across more than three orders of magnitude in aCDOM(440) values (0.001–2.305 m􀀀1). EMA using UV and NIR wavelengths (320 and 780 nm) is applied to low-altitude airborne observations and satisfies 25% uncertainty based on unbiased percent differences (UPDs) within each of three dissimilar match-up sites ranging in aCDOM(440) from 0.02–0.57 m􀀀1. Results demonstrate that EMA is a useful and robust approach for the remote sensing of aCDOM(440) in coastal and inland waters, which are generally shallower, contain more optically complex environments, and span a greater range in aCDOM(440) than oceanic waters.

CDOM↗

Aerosol-Induced Changes of Convective Cloud Anvils Produce Strong Climate Warming

The effect of aerosol on clouds poses one of the largest uncertainties in estimating the anthropogenic contribution to climate change. Small human-induced perturbations to cloud characteristics via aerosol pathways can create a change in the top-of-atmosphere radiative forcing of hundreds of Wm(exp-2) . Here we focus on links between aerosol and deep convective clouds of the Atlantic and Pacific Intertropical Convergence Zones, noting that the aerosol environment in each region is entirely different. The tops of these vertically developed clouds consisting of mostly ice can reach high levels of the atmosphere, overshooting the lower stratosphere and reaching altitudes greater than 16 km. We show a link between aerosol, clouds and the free atmosphere wind profile that can change the magnitude and sign of the overall climate radiative forcing. We find that increased aerosol loading is associated with taller cloud towers and anvils. The taller clouds reach levels of enhanced wind speeds that act to spread and thin the anvi1 clouds, increasing areal coverage and decreasing cloud optical depth. The radiative effect of this transition is to create a positive radiative forcing (warming) at top-of-atmosphere. Furthermore we introduce the cloud optical depth (r), cloud height (Z) forcing space and show that underestimation of radiative forcing is likely to occur in cases of non homogenous clouds. Specifically, the mean radiative forcing of towers and anvils in the same scene can be several times greater than simply calculating the forcing from the mean cloud optical depth in the scene. Limitations of the method are discussed, alternative sources of aerosol loading are tested and meteorological variance is restricted, but the trend of taller clouds; increased and thinner anvils associated with increased aerosol loading remains robust through all the different tests and perturbations.

Koren, I.↗

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z = 4 to 8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Hubble Frontier Field deep parallel observations near the Abell 2744 and MACS J0416.1- 2403 clusters. The combination of these surveys provides an effective volume of 0.6-1.2 ×10(exp 6) Mpc(exp 3) over this epoch, allowing us to perform a robust search for bright (M(sub UV) less than −21) and faint (M(sub UV) = −18) galaxies. We select galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5 less than z less than 8.5, with more than 1000 galaxies at z of approximately 6 - 8. We measure both a stepwise luminosity function for galaxies in our redshift samples, as well as a Schechter function, using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our UV luminosity functions agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z of greater than or equal to 6. Our bestfit value of the characteristic magnitude M* is consistent with −21 at z of greater than or equal to 5, different than that inferred based on previous trends at lower redshift. At z = 8, a single power-law provides an equally good fit to the UV luminosity function, while at z = 6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare our luminosity functions to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift, though a decreasing impact of feedback may also be possible. We measure the evolution of the cosmic star-formation rate (SFR) density by integrating our observed luminosity functions to M(sub UV) = −17, correcting for dust attenuation, and find that the SFR density declines proportionally to (1 + z)((exp −4.3)(+/-)(0.5)) at z greater than 4, consistent with observations at z greater than or equal to 9. Our observed luminosity functions are consistent with a reionization history that starts at redshift of approximately greater than 10, completes at z greater than 6, and reaches a midpoint (x(sub HII) = 0.5) at 6.7 less than z less than 9.4. Finally, using a constant cumulative number density selection and an empirically derived rising star-formation history, our observations predict that the abundance of bright z = 9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z similar to 10 galaxies.

Galaxy↗

Constraints on the Martian cratering rate imposed by the SNC meteorites and Vallis Marineris layered deposits

Following two independent lines of evidence -- estimates of the age and formation time of a portion of the Martian geologic column exposed in the layered deposits and the crystallization and ejection ages of the SNC meteorites -- it appears that the Martian cratering rate must be double the lunar rate or even higher. This means models such as NHII or NHIII (Neukum and Hiller models II and III), which estimate the Martian cratering rate as being several times lunar are probably far closer to reality on Mars than lunar rates. The effect of such a shift is profound: Mars is transformed from a rather Moon-like place into a planet with vigorous dynamics, multiple large impacts, erosion, floods, and volcanism throughout its history. A strong shift upward in cratering rates on Mars apparently solves some glaring problems; however, it creates others. The period of time during which Earth-like atmospheric conditions existed, the liquid water era on Mars, persists in NHIII up to only 0.5 b.y. ago. Scenarios of extended Earth-like conditions on Mars have been discounted in the past because they would have removed many of the craters from the early bombardment era found in the south. It does appear that some process of crater removal was quite vigorous in the north during Mars' past. Evidence exists that the northern plains may have been the home of long-lived seas or perhaps even a paleo-ocean, so models exist for highly localized destruction of craters in the north. However, the question of how the ancient crater population could be preserved in the south under a long liquid-water era found in any high-cratering-rate models is a serious question that must be addressed. It does appear to be a higher-order problem because it involves low-energy dynamics acting in localized areas, i.e., erosion of craters in the south of Mars, whereas the two problems with the low-cratering-rate models involve high-energy events acting over large areas: the formation of the Vallis Marineris, the SNC ejecting impacts, and the global atmospheric pressure and temperature conditions that allow liquid water to exist as a robust entity anywhere on the Martian surface. In any case, it appears Mars is a more complex and dynamic planet than previously supposed. It has canyons dating from the middle to late period of its history that contain apparent lake sediments bedded deeper than most sediments on Earth. Recent multiple, violent impacts on Mars have apparently provided us with multiple random samples of its surface that all crystallized less than 1.5 b.y. ago. These things cannot be accommodated in our present cratering chronologies of Mars, based on 1x lunar cratering rates, without great difficulties. These difficulties suggest that a new chronology, probably based on NHII or even NHIII, should be adopted; this new chronology will provide us with a new view of Mars as a dynamic planet of rich history.

Brandenburg, J. E.↗

Distributed Prognostics and Health Management with a Wireless Network Architecture

A heterogeneous set of system components monitored by a varied suite of sensors and a particle-filtering (PF) framework, with the power and the flexibility to adapt to the different diagnostic and prognostic needs, has been developed. Both the diagnostic and prognostic tasks are formulated as a particle-filtering problem in order to explicitly represent and manage uncertainties in state estimation and remaining life estimation. Current state-of-the-art prognostic health management (PHM) systems are mostly centralized in nature, where all the processing is reliant on a single processor. This can lead to a loss in functionality in case of a crash of the central processor or monitor. Furthermore, with increases in the volume of sensor data as well as the complexity of algorithms, traditional centralized systems become for a number of reasons somewhat ungainly for successful deployment, and efficient distributed architectures can be more beneficial. The distributed health management architecture is comprised of a network of smart sensor devices. These devices monitor the health of various subsystems or modules. They perform diagnostics operations and trigger prognostics operations based on user-defined thresholds and rules. The sensor devices, called computing elements (CEs), consist of a sensor, or set of sensors, and a communication device (i.e., a wireless transceiver beside an embedded processing element). The CE runs in either a diagnostic or prognostic operating mode. The diagnostic mode is the default mode where a CE monitors a given subsystem or component through a low-weight diagnostic algorithm. If a CE detects a critical condition during monitoring, it raises a flag. Depending on availability of resources, a networked local cluster of CEs is formed that then carries out prognostics and fault mitigation by efficient distribution of the tasks. It should be noted that the CEs are expected not to suspend their previous tasks in the prognostic mode. When the prognostics task is over, and after appropriate actions have been taken, all CEs return to their original default configuration. Wireless technology-based implementation would ensure more flexibility in terms of sensor placement. It would also allow more sensors to be deployed because the overhead related to weights of wired systems is not present. Distributed architectures are furthermore generally robust with regard to recovery from node failures.

Goebel, Kai↗

Updated MISR Over-Water Research Aerosol Retrieval Algorithm - Part 2: A Multi-Angle Aerosol Retrieval Algorithm for Shallow, Turbid, Oligotrophic, and Eutrophic Waters

Coastal waters serve as transport pathways to the ocean for all agricultural and other runoff from terrestrial sources, and many are the sites for upwelling of nutrient rich, deep water; they are also some of the most biologically productive on Earth. Estimating the impact coastal waters have on the global carbon budget requires relating satellite-based remote-sensing retrievals of biological productivity (e.g., chlorophyll alpha concentration) to in situ measurements taken in near-surface waters. The Multi-angle Imaging SpectroRadiometer (MISR) can uniquely constrain the "atmospheric correction" needed to derive ocean color from remote-sensing imagers. Here, we retrieve aerosol amount and type from MISR over all types of water. The primary limitation is an upper bound on aerosol optical depth (AOD), as the algorithm must be able to distinguish the surface. This updated MISR research aerosol retrieval algorithm (RA) also assumes that light reflection by the underlying ocean surface is Lambertian. The RA computes the ocean surface reflectance (R (sub rs)) analytically for a given AOD, aerosol optical model, and wind speed. We provide retrieval examples over shallow, turbid, and eutrophic waters and introduce a productivity and turbidity index (PTI), calculated from retrieved spectral R (sub rs), that distinguished water types (similar to the the normalized difference vegetation index, NDVI, over land). We also validate the new algorithm by comparing spectral AOD and Angstrom exponent (ANG) results with 2419 collocated AErosol RObotic NETwork (AERONET) observations. For AERONET 558 nanometer-interpolated AOD less than 1.0, the root-mean square error (RMSE) is 0.04 and linear correlation coefficient is 0.95. For the 502 cloud-free MISR and AERONET collocations with an AERONET AOD greater than 0.20, the ANG RMSE is 0.25 and r is 0.89. Although MISR RA AOD retrieval quality does not appear to be substantially impacted by the presence of turbid water, the MISR-RA-retrieved Angstrom exponent seems to suffer from increased uncertainty under such conditions. MISR supplements current ocean color sources in regions where sunglint precludes retrievals from single-view-angle instruments. MISR atmospheric correction should also be more robust than that derived from single-view instruments such as the Moderate Resolution Imaging Spectroradiometer (MODIS). This is especially true in regions of shallow, turbid, and eutrophic waters, locations where biological productivity can be high, and single-view-angle retrieval algorithms struggle to separate atmospheric from oceanic features.

Satellite Remote Sensing↗

Atmospheric corrections for TIMS estimated emittance

The estimated temperature of the average of 500 lines of Thermal Infrared Multispectral Scanner (TIMS) data of the Pacific Ocean, from flight line 94, collected on 30 Sep. 1988, at 1931 GMT is shown. With no atmospheric corrections, estimated temperature decreases away from nadir (the center of the scan line). A LOWTRAN modeled correction, using local radiosonde data and instrument scan angle information, results in reversed limb darkening effects for most bands, and does not adequately correct all bands to the same temperature. The atmosphere tends to re-radiate energy at the wavelengths at which it most absorbs, and thus the overall difference between corrected and uncorrected temperatures is approximately 40 C, despite the average LOWTRAN calculated transmittance of only 60 percent between 8.1 and 11.6 microns. An alternative approach to atmospheric correction is a black body normalization. This is done by calculating a normalization factor for each pixel position and wavelength, which when applied results in a single calculated temperature, as would be expected for a gray body with near uniform emittance. The black body adjustment is based on the atmospheric conditions over the sea. The ground elevation profile along the remaining 3520 scan lines (approximately 10 km) of flight line 94, up the slopes of Kilauea, determined from aircraft pressure and laser altimeter data is shown. This flight line includes a large amount of vegetation that is clearly discernible on the radiance image, being much cooler than the surrounding rocks. For each of the 3520 scan lines, pixels were classified as vegetation or 'other'. A moving average of 51 lines was applied to the composite vegetation emittance for each scan line, to reduce noise. Assuming vegetation to be like water, and to act as gray body with an emittance of 0.986 across the spectrum, it is shown that that the LOWTRAN induced artifacts are severe, and other than for the 0.9.9 micron channel, not significantly different from applying no corrections at all. As expected, with increasing elevation atmospheric effects are slightly reduced, because moisture tends to be concentrated in the lowermost part of the atmosphere. The black body adjustment is highly robust, and even at elevations nearly 600 meters above the sea, remains an alternative procedure for use in calculating emittance.

Warner, T. A.↗

Spacecraft Modeling, Attitude Determination, and Control: Quaternion-Based Approach

Spacecraft attitude determination and control is an important part of a spacecraft to achieve its designed mission. As of today, many spacecrafts have been successfully launched, and most of them have performed well as they were designed. Many research papers have been published to address the attitude determination and control design problems. Several text books are available for students to learn the technology and for engineers to use as references. The most popular spacecraft models for attitude determination algorithms and control design methods are the Euler angle models and the quaternion models . The Euler angle models have been proved very efficient as the linearized models are controllable, and all standard linear control system design methods are directly applicable. The drawbacks related to the Euler angle methods are (a) the designs based on linearized models may not globally stabilize the original nonlinear spacecraft, i.e., the design may not work when the attitude of the spacecraft is far away from the point where the linearization is performed; (b) the models depend on the rotational sequences, this can be error prone if several teams work on the same project and they use different rotational sequences; (c) for any rotational sequence, there is a singular point where the model is not applicable; and (d) since most attitude determination methods use quaternion to represent the spacecraft attitude, there is a need to transform quaternion into Euler angles. On the other hand, for quaternion models, people have found controllers that can globally stabilize nonlinear spacecraft systems; the models do not depend on rotational sequences and they have no singular point; and the quaternion is provided by attitude determination system and ready to use. The main problem with the quaternion model based control system design is that the linearized quaternion model is not controllable. Therefore, most published design methods heavily rely on Lyapunov functions for the nonlinear spacecraft system. But there is no systematic way to obtain a desired Lyapunov functions. Moreover, the Lyapunov function based designs focus on the closed-loop system stability but pay little attention to the closed-loop system performance. In a series of papers, the author proposed some reduced quaternion models which lead to some controllable linearized spacecraft models. Therefore, all standard linear system theory can be directly applied to analyze and design the spacecraft control systems. We showed that, in some cases, the designed control system is not only optimal for the linearized system, but also globally stabilize the original nonlinear system . Clearly, the reduced quaternion models do not depend on rotational sequences. Due to the special structure of the linearized spacecraft model, some most important design methods, such as LQR design and robust pole assignment design are very simple, enjoy the analytical solutions for some problems, have direct connection to the performance measures, such as settling time, rising time , and percentage of overshoot . All these features are attractive for high quality control system designs. The idea mentioned above is then extended to more spacecraft control problems using specific actuators such as magnetic torque bars and control momentum gyroscopes. These types of actuators may not provide exactly desired torques. Most existing methods use different conversions to get approximate solutions, meaning that these actuators may generate a torque close to but not equal to the desired one. Using the reduced quaternion models that incorporate the actuators into the system model, the control inputs are not torques but the operational parameters. The main benefit of this idea is that the control actions are not approximate but accurate. As all actuators have their operational limit, design with input constraints are also considered in this book by using recently developed interior-point optimization techniques. This book grows up from my research on the spacecraft attitude determination and control design methods in more than a decade which is focused on using reduced quaternion models because of their merits stated above. The book provides all necessary background materials on orbital dynamics, rotations and quaternion, frequently used reference frames, transformations between reference frames, space environment and disturbance torques, ephemeris astronomical vector calculations and measurement instruments, spacecraft control actuators and their models, so that the readers will get a global picture and can apply all these information into the spacecraft system modeling, attitude determination, and spacecraft control system designs, which is the main purpose of this book. This book is different from existing books in that we focus on quaternion based spacecraft control system designs and we consider only attitude control system design related problems, from spacecraft modeling, to attitude determination and estimation, to control system design method selection, to control algorithm development, and to the simulation of the control system designs. Moreover, this book addresses different attitude control tasks in the spacecraft life cycle, including spacecraft maneuver, orbit raising, attitude control, and rendezvous. Finally, this book emphasizes the state space design methods rather than the classical frequency design methods.

Attitude determination and control↗

Diurnal Variability of Vertical Structure from a TRMM Passive Microwave "Virtual Radar" Retrieval

Robust description of the diurnal cycle from TRMM observations is complicated by the limitations of Low Earth Orbit (LEO) sampling; from a 'climatological' perspective, sufficient sampling must exist to control for both spatial and seasonal variability, before tackling an additional diurnal component (e.g., with 8 additional 3-hourly or 24 1-hourly bins). For documentation of vertical structure, the narrow sample swath of the TRMM Precipitation Radar limits the resolution of any of these components. A neural-network based 'virtual radar" retrieval has been trained and internally validated, using multifrequency / multipolarization passive microwave(TM1) brightness temperatures and textures parameters and lightning (LIS) observations, as inputs, and PR volumetric reflectivity as targets (outputs). By training the algorithms (essentially highly multivariate, nonlinear regressions) on a very large sample of high-quality co-located data from the center of the TRMM swath, 3D radar reflectivity and derived parameters (VIL, IWC, Echo Tops, etc.) can be retrieved across the entire TMI swath, good to 8-9% over the dynamic range of parameters. As a step in the retrieval (and as an output of the process), each TMI multifrequency pixel (at 85 GHz resolution) is classified into one of the 25 archetypal radar profile vertical structure "types", previously identified using cluster analysis. The dynamic range of retrieved vertical structure appears to have higher fidelity than the current (Version 6) experimental GPROF hydrometeor vertical structure retrievals. This is attributable to correct representation of the prior probabilities of vertical structure variability in the neural network training data, unlike the GPROF cloud-resolving model training dataset used in the V6 algorithms. The LIS lightning inputs are supplementary inputs, and a separate offline neural network has been trained to impute (predict) LIS lightning from passive-microwave-only data. The virtual radar retrieval is thus, in principle, extensible to Aqua/AMSR-E and NPOESS/CMIS passive microwave instruments. The virtual radar approach yields a threefold increase in effective sampling from the mission, albeit of lower-quality "retrieved" data, reducing the variance of local estimates by one third (or the standard deviation by-0.57). In this talk, the variance reduction is leveraged to more finely resolve global diurnal variability in both space and time (local hour).

Boccippio, Dennis J.↗

Progress in Development of the Axel Rovers

Progress has been made in the development of a family of robotic land vehicles having modular and minimalist design features chosen to impart a combination of robustness, reliability, and versatility. These vehicles at earlier stages of development were described in two previous NASA Tech Briefs articles: "Reconfigurable Exploratory Robotic Vehicles" (NPO-20944), Vol. 25, No. 7 (July 2001), page 56; and "More About Reconfigurable Exploratory Robotic Vehicles" (NPO-30890), Vol. 33, No. 8 (August 2009), page 40. Conceived for use in exploration of the surfaces of Mars and other remote planets, these vehicles could also be adapted to terrestrial applications, including exploration of volcanic craters or other hostile terrain, military reconnaissance, inspection of hazardous sites, and searching for victims of earthquakes, landslides, avalanches, or mining accidents. In addition, simplified versions of these vehicles might be marketable as toys. The most basic module in this family of reconfigurable robots is the Axel rover, which has a cylindrical body with two main wheels and a trailing link. Inside its body are three motors and associated mechanisms for driving the two wheels and for rotating the link 360 around its symmetrical body. The actuated link serves several purposes: It is used as a lever arm to react to the wheels thrust to move Axel in multiple directions. It is used to rotate the Axel housing in order to tilt, to the desired angle, any sensors and instruments mounted on or in the Axel housing. It provides an alternative mobility mode, which is primarily used in its tethered configuration. Turn ing the link into the ground in lieu of driving the wheels causes the Axel housing and wheels to roll as a unit and thereby leads to a tumbling motion along the ground. With a tether mounted around Axel s cylindrical body, the link serves as a winch mechanism to reel and unreel the tether raising and lowering Axel over steep and vertical surfaces (Figure 1). Sensors, computation, and communication modules are also housed inside Axel s body. A pair of stereo vision cameras provides three-dimensional view for autonomous navigation and avoiding obstacles. Inertial sensors determine the tilt of the robot and are used for estimating its motion. In a fully developed version, power would be supplied by rechargeable batteries aboard Axel; at the time of reporting the information for this article, power was supplied from an external source via a cable. In and of itself, the Axel rover is fully capable of traversing and sampling terrains on planetary surfaces. By use of only the two main wheel actuators and the caster link actuator, Axel can be made to follow an arbitrary path, turn in place, and operate upside- down or right-side-up. If operated in a tethered configuration, as shown in Figure 1, it can be made to move down and up a steep crater wall, descend from an overhang to a cave, and ascend from the cave back to the overhang, all by use of the same three actuators. Such tethered operation could be useful in searching for accident victims or missing persons in mines, caves, and rubble piles. Running the tether through the caster link enhances the stability of Axel and provides a restoring force that keeps the link off the ground for the most part during operation on a steep slope. In its extended configuration, two Axel modules can dock to either side of a payload module to form the four wheeled Axel2 rover (Figure 2). Additional payload and Axel modules can dock to either side of the Axel2 to form the Axel3 rover, extending its payload capacity and its mobility capabilities.

Nesnas, Issa A.↗