Search NASA⌕ Search

SEARCH · Search NASA

Results for “global uncertainty analysis”

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 145 records · Page 8

Challenges Facing Design and Analysis Tools

The design and analysis of future aerospace systems will strongly rely on advanced engineering analysis tools used in combination with risk mitigation procedures. The implications of such a trend place increased demands on these tools to assess off-nominal conditions, residual strength, damage propagation, and extreme loading conditions in order to understand and quantify these effects as they affect mission success. Advances in computer hardware such as CPU processing speed, memory, secondary storage, and visualization provide significant resources for the engineer to exploit in engineering design. The challenges facing design and analysis tools fall into three primary areas. The first area involves mechanics needs such as constitutive modeling, contact and penetration simulation, crack growth prediction, damage initiation and progression prediction, transient dynamics and deployment simulations, and solution algorithms. The second area involves computational needs such as fast, robust solvers, adaptivity for model and solution strategies, control processes for concurrent, distributed computing for uncertainty assessments, and immersive technology. Traditional finite element codes still require fast direct solvers which when coupled to current CPU power enables new insight as a result of high-fidelity modeling. The third area involves decision making by the analyst. This area involves the integration and interrogation of vast amounts of information - some global in character while local details are critical and often drive the design. The proposed presentation will describe and illustrate these areas using composite structures, energy-absorbing structures, and inflatable space structures. While certain engineering approximations within the finite element model may be adequate for global response prediction, they generally are inadequate in a design setting or when local response prediction is critical. Pitfalls to be avoided and trends for emerging analysis tools will be described.

Knight, Norman F., Jr.↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Site Distribution and Aliasing Effects in the Inversion for Load Coefficients and Geocenter Motion from GPS Data

Precise GPS measurements of elastic relative site displacements due to surface mass loading offer important constraints on global surface mass transport. We investigate effects of site distribution and aliasing by higher-degree (n greater than or equal 2) loading terms on inversion of GPS data for n = 1 load coefficients and geocenter motion. Covariance and simulation analyses are conducted to assess the sensitivity of the inversion to aliasing and mismodeling errors and possible uncertainties in the n = 1 load coefficient determination. We found that the use of center-of-figure approximation in the inverse formulation could cause 10- 15% errors in the inverted load coefficients. n = 1 load estimates may be contaminated significantly by unknown higher-degree terms, depending on the load scenario and the GPS site distribution. The uncertainty in n = 1 zonal load estimate is at the level of 80 - 95% for two load scenarios.

global positioning system↗

Determination of Interannual to Decadal Changes in Ice Sheet Mass Balance from Satellite Altimetry

A major uncertainty in predicting sea level rise is the sensitivity of ice sheet mass balance to climate change, as well as the uncertainty in present mass balance. Since the annual water exchange is about 8 mm of global sea level equivalent, the +/- 25% uncertainty in current mass balance corresponds to +/- 2 mm/yr in sea level change. Furthermore, estimates of the sensitivity of the mass balance to temperature change range from perhaps as much as - 10% to + 10% per K. Although the overall ice mass balance and seasonal and inter-annual variations can be derived from time-series of ice surface elevations from satellite altimetry, satellite radar altimeters have been limited in spatial coverage and elevation accuracy. Nevertheless, new data analysis shows mixed patterns of ice elevation increases and decreases that are significant in terms of regional-scale mass balances. In addition, observed seasonal and interannual variations in elevation demonstrate the potential for relating the variability in mass balance to changes in precipitation, temperature, and melting. From 2001, NASA's ICESat laser altimeter mission will provide significantly better elevation accuracy and spatial coverage to 86 deg latitude and to the margins of the ice sheets. During 3 to 5 years of ICESat-1 operation, an estimate of the overall ice sheet mass balance and sea level contribution will be obtained. The importance of continued ice monitoring after the first ICESat is illustrated by the variability in the area of Greenland surface melt observed over 17-years and its correlation with temperature. In addition, measurement of ice sheet changes, along with measurements of sea level change by a series of ocean altimeters, should enable direct detection of ice level and global sea level correlations.

Zwally, H. Jay↗

Coherent Evaluation of Aerosol Data Products from Multiple Satellite Sensors

Aerosol retrieval from satellite has practically become routine, especially during the last decade. However, there is often disagreement between similar aerosol parameters retrieved from different sensors, thereby leaving users confused as to which sensors to trust for answering important science questions about the distribution, properties, and impacts of aerosols. As long as there is no consensus, and the inconsistencies are not well characterized and understood, there will be no way of developing reliable model inputs and climate data records from satellite aerosol measurements. Fortunately, the Aerosol Robotic Network (AERONET) is providing well-calibrated globally representative ground-based aerosol measurements corresponding to the satellite-retrieved products. Through a recently developed web-based Multi-sensor Aerosol Products Sampling System (MAPSS), we are utilizing the advantages offered by collocated AERONET and satellite products to characterize and evaluate aerosol retrieval from multiple sensors. Indeed, MAPSS and its companion statistical tool AeroStat are facilitating detailed comparative uncertainty analysis of satellite aerosol measurements from Terra-MODIS, Aqua-MODIS, Terra-MISR, Aura-OMI, Parasol-POLDER, and Calipso-CALIOP. In this presentation, we will describe the strategy of the MAPSS system, its potential advantages for the aerosol community, and the preliminary results of an integrated comparative uncertainly analysis of aerosol products from multiple satellite sensors.

Ichoku, Charles↗

Uncertainties in Predicting Rice Yield by Current Crop Models Under a Wide Range of Climatic Conditions

Predicting rice (Oryza sativa) productivity under future climates is important for global food security. Ecophysiological crop models in combination with climate model outputs are commonly used in yield prediction, but uncertainties associated with crop models remain largely unquantified. We evaluated 13 rice models against multi-year experimental yield data at four sites with diverse climatic conditions in Asia and examined whether different modeling approaches on major physiological processes attribute to the uncertainties of prediction to field measured yields and to the uncertainties of sensitivity to changes in temperature and CO2 concentration [CO2]. We also examined whether a use of an ensemble of crop models can reduce the uncertainties. Individual models did not consistently reproduce both experimental and regional yields well, and uncertainty was larger at the warmest and coolest sites. The variation in yield projections was larger among crop models than variation resulting from 16 global climate model-based scenarios. However, the mean of predictions of all crop models reproduced experimental data, with an uncertainty of less than 10 percent of measured yields. Using an ensemble of eight models calibrated only for phenology or five models calibrated in detail resulted in the uncertainty equivalent to that of the measured yield in well-controlled agronomic field experiments. Sensitivity analysis indicates the necessity to improve the accuracy in predicting both biomass and harvest index in response to increasing [CO2] and temperature.

crop-model ensembles↗

Synergy of Observations and Dynamo Models to Understand and Predict Solar Activity Cycles

The long-standing problem of understanding the evolution of the global magnetic fields that drive solar activity through different temporal scales is becoming more tractable because, in addition to 400 years of sunspot records, we now have almost 4 solar cycles of magnetic field observations. These observations allow us to discern physical connections between dynamo model variables and observations using data assimilation analysis. In particular, the Ensemble Kalman Filter approach takes into account uncertainties in both observations and modeling and allows us to make reliable forecasts of solar cycle activity cycles by using a relatively simple non-linear dynamical model of the solar dynamo. To expand this approach for more complex 2D and 3D dynamo modeling, it is necessary to decompose the observed synoptic magnetograms into poloidal and toroidal field components. In this presentation I will present initial results on magnetogram decomposition and assimilation of magnetogram data into dynamo modeling.

Synergy↗

Synergy of Observations and Dynamo Models to Understand and Predict Solar Activity Cycles

The long-standing problem of understanding the evolution of the global magnetic fields that drive solar activity through different temporal scales is becoming more tractable because, in addition to 400 years of sunspot records, we now have almost 4 solar cycles of magnetic field observations. These observations allow us to discern physical connections between dynamo model variables and observations using data assimilation analysis. In particular, the Ensemble Kalman Filter approach takes into account uncertainties in both observations and modeling and allows us to make reliable forecasts of solar cycle activity by using a relatively simple non-linear dynamical model of the solar dynamo. To expand this approach for more complex 2D and 3D dynamo modeling, it is necessary to decompose the observed synoptic magnetograms into poloidal and toroidal field components. In this presentation I will present initial results on magnetogram decomposition and assimilation of magnetogram data into dynamo modeling.

Kitiashvili, Irina N.↗

Challenges to Aboveground Biomass Prediction from Waveform Lidar

Accurate accounting of aboveground biomass density (AGBD) is crucial for carbon cycle, biodiversity, and climate change science. The Global Ecosystem Dynamics Investigation (GEDI), which maps global AGBD from waveform lidar, is the first of a new generation of Earth observation missions designed to improve carbon accounting. This paper explores the possibility that lidar waveforms may not be unique to AGBD—that forest stands with different AGBD may produce highly similar waveforms—and we hypothesize that non-uniqueness may contribute to the large uncertainties in AGBD predictions. Our analysis integrates simulated GEDI waveforms from 428 in situ stem maps with output from an individual-based forest gap model, which we use to generate a database of potential forest stands and simulate GEDI waveforms from those stands. We use this database to predict the AGBD of the 428 in situ stem maps via two different methods: a linear regression from waveform metrics, and a waveform-matching approach that accounts for waveform-AGBD non-uniqueness. We find that some in situ waveforms are more unique to AGBD than others, which notably impacts AGBD prediction uncertainty (7–411 Mg ha−1, average of 167 Mg ha−1). We also find that forest structure complexity may influence the non-uniqueness effect; stands with low structural complexity are more unique to AGBD than more mature stands with multiple cohorts and canopy layers. These findings suggest that the non-uniqueness phenomena may be introduced by the measuring characteristics of waveform lidar in combination with how forest structure manifests at small scales, and we discuss how this complexity may complicate uncertainty estimation in AGBD prediction. This analysis suggests a limit to the accuracy and precision of AGBD predictions from lidar waveforms seen in empirical studies, and underscores the need for further exploration of the relationships between lidar remote sensing measurements, forest structure, and AGBD.

GEDI↗

Study of geopotential error models used in orbit determination error analysis

The uncertainty in the geopotential model is currently one of the major error sources in the orbit determination of low-altitude Earth-orbiting spacecraft. The results of an investigation of different geopotential error models and modeling approaches currently used for operational orbit error analysis support at the Goddard Space Flight Center (GSFC) are presented, with emphasis placed on sequential orbit error analysis using a Kalman filtering algorithm. Several geopotential models, known as the Goddard Earth Models (GEMs), were developed and used at GSFC for orbit determination. The errors in the geopotential models arise from the truncation errors that result from the omission of higher order terms (omission errors) and the errors in the spherical harmonic coefficients themselves (commission errors). At GSFC, two error modeling approaches were operationally used to analyze the effects of geopotential uncertainties on the accuracy of spacecraft orbit determination - the lumped error modeling and uncorrelated error modeling. The lumped error modeling approach computes the orbit determination errors on the basis of either the calibrated standard deviations of a geopotential model's coefficients or the weighted difference between two independently derived geopotential models. The uncorrelated error modeling approach treats the errors in the individual spherical harmonic components as uncorrelated error sources and computes the aggregate effect using a combination of individual coefficient effects. This study assesses the reasonableness of the two error modeling approaches in terms of global error distribution characteristics and orbit error analysis results. Specifically, this study presents the global distribution of geopotential acceleration errors for several gravity error models and assesses the orbit determination errors resulting from these error models for three types of spacecraft - the Gamma Ray Observatory, the Ocean Topography Experiment, and the Cosmic Background Explorer.

Yee, C.↗

Satellite Remote Sensing of Aerosol Forcing

The role of aerosol forcing remains one of the largest uncertainties in estimating man's impact on the global climate system. One school of thought suggests that remote sensing by satellite sensors will provide the data necessary to narrow these uncertainties. While satellite measurements of direct aerosol forcing appear to be straightforward, satellite measurements of aerosol indirect forcing will be more complicated. Pioneering studies identified indirect aerosol forcing using AVHRR data in the biomass burning regions of Brazil. We have expanded this analysis with AVHRR to include an additional year of data and assimilated water vapor fields. The results show similar latitudinal dependence as reported by Kaufman and Fraser, but by using water vapor observations we conclude that latitude is not a proxy for water vapor and the strength of the indirect effect is not correlated to water vapor amounts. In addition to the AVHRR study we have identified indirect aerosol forcing in Brazil at much smaller spatial scales using the MODIS Airborne Simulator. The strength of the indirect effect appears to be related to cloud type and cloud dynamics. There is a suggestion that some of the cloud dynamics may be influenced by smoke destabilization of the atmospheric column. Finally, this study attempts to quantify remote sensing limitations due to the accuracy limits of the retrieval algorithms. We use a combination of numerical aerosol transport models, ground-based AERONET data and ISCCP cloud climatology to determine how much of the forcing occurs in regions too clean to determine from satellite retrievals.

Remer, Lorraine↗

Multi-Factor Impact Analysis of Agricultural Production in Bangladesh with Climate Change

Diverse vulnerabilities of Bangladesh's agricultural sector in 16 sub-regions are assessed using experiments designed to investigate climate impact factors in isolation and in combination. Climate information from a suite of global climate models (GCMs) is used to drive models assessing the agricultural impact of changes in temperature, precipitation, carbon dioxide concentrations, river floods, and sea level rise for the 2040-2069 period in comparison to a historical baseline. Using the multi-factor impacts analysis framework developed in Yu et al. (2010), this study provides new sub-regional vulnerability analyses and quantifies key uncertainties in climate and production. Rice (aman, boro, and aus seasons) and wheat production are simulated in each sub-region using the biophysical Crop Environment REsource Synthesis (CERES) models. These simulations are then combined with the MIKE BASIN hydrologic model for river floods in the Ganges-Brahmaputra-Meghna (GBM) Basins, and the MIKE21Two-Dimensional Estuary Model to determine coastal inundation under conditions of higher mean sea level. The impacts of each factor depend on GCM configurations, emissions pathways, sub-regions, and particular seasons and crops. Temperature increases generally reduce production across all scenarios. Precipitation changes can have either a positive or a negative impact, with a high degree of uncertainty across GCMs. Carbon dioxide impacts on crop production are positive and depend on the emissions pathway. Increasing river flood areas reduce production in affected sub-regions. Precipitation uncertainties from different GCMs and emissions scenarios are reduced when integrated across the large GBM Basins' hydrology. Agriculture in Southern Bangladesh is severely affected by sea level rise even when cyclonic surges are not fully considered, with impacts increasing under the higher emissions scenario.

Ruane, Alex C.↗

Lidar Measurements of Atmospheric CO2 From Regional to Global Scales

Atmospheric CO2 is a critical forcing for the Earth's climate and the knowledge on its distributions and variations influences predictions of the Earth's future climate. Large uncertainties in the predictions persist due to limited observations. This study uses the airborne Intensity-Modulated Continuous-Wave (IMCW) lidar developed at NASA Langley Research Center to measure regional atmospheric CO2 spatio-temporal variations. Further lidar development and demonstration will provide the capability of global atmospheric CO2 estimations from space, which will significantly advances our knowledge on atmospheric CO2 and reduce the uncertainties in the predictions of future climate. In this presentation, atmospheric CO2 column measurements from airborne flight campaigns and lidar system simulations for space missions will be discussed. A measurement precision of approx.0.3 ppmv for a 10-s average over desert and vegetated surfaces has been achieved. Data analysis also shows that airborne lidar CO2 column measurements over these surfaces agree well with in-situ measurements. Even when thin cirrus clouds present, consistent CO2 column measurements between clear and thin cirrus cloudy skies are obtained. Airborne flight campaigns have demonstrated that precise atmospheric column CO2 values can be measured from current IM-CW lidar systems, which will lead to use this airborne technique in monitoring CO2 sinks and sources in regional and continental scales as proposed by the NASA Atmospheric Carbon and Transport â€" America project. Furthermore, analyses of space CO2 measurements shows that applying the current IM-CW lidar technology and approach to space, the CO2 science goals of space missions will be achieved, and uncertainties in CO2 distributions and variations will be reduced.

Lin, Bing↗

Nominal 30-M Cropland Extent Map of Continental Africa by Integrating Pixel-Based and Object-Based Algorithms Using Sentinel-2 and Landsat-8 Data on Google Earth Engine

A satellite-derived cropland extent map at high spatial resolution (30-m or better) is a must for food and water security analysis. Precise and accurate global cropland extent maps, indicating cropland and non-cropland areas, is a starting point to develop high-level products such as crop watering methods (irrigated or rainfed), cropping intensities (e.g., single, double, or continuous cropping), crop types, cropland fallows, as well as assessment of cropland productivity (productivity per unit of land), and crop water productivity (productivity per unit of water). Uncertainties associated with the cropland extent map have cascading effects on all higher-level cropland products. However, precise and accurate cropland extent maps at high spatial resolution over large areas (e.g., continents or the globe) are challenging to produce due to the small-holder dominant agricultural systems like those found in most of Africa and Asia. Cloud-based Geospatial computing platforms and multi-date, multi-sensor satellite image inventories on Google Earth Engine offer opportunities for mapping croplands with precision and accuracy over large areas that satisfy the requirements of broad range of applications. Such maps are expected to provide highly significant improvements compared to existing products, which tend to be coarser in resolution, and often fail to capture fragmented small-holder farms especially in regions with high dynamic change within and across years. To overcome these limitations, in this research we present an approach for cropland extent mapping at high spatial resolution (30-m or better) using the 10-day, 10 to 20-m, Sentinel-2 data in combination with 16-day, 30-m, Landsat-8 data on Google Earth Engine (GEE). First, nominal 30-m resolution satellite imagery composites were created from 36,924 scenes of Sentinel-2 and Landsat-8 images for the entire African continent in 2015-2016. These composites were generated using a median-mosaic of five bands (blue, green, red, near-infrared, NDVI) during each of the two periods (period 1: January-June 2016 and period 2: July-December 2015) plus a 30-m slope layer derived from the Shuttle Radar Topographic Mission (SRTM) elevation dataset. Second, we selected Cropland/Non-cropland training samples (sample size 9791) from various sources in GEE to create pixel-based classifications. As supervised classification algorithm, Random Forest (RF) was used as the primary classifier because of its efficiency, and when over-fitting issues of RF happened due to the noise of input training data, Support Vector Machine (SVM) was applied to compensate for such defects in specific areas. Third, the Recursive Hierarchical Segmentation (RHSeg) algorithm was employed to generate an object-oriented segmentation layer based on spectral and spatial properties from the same input data. This layer was merged with the pixel-based classification to improve segmentation accuracy. Accuracies of the merged 30-m crop extent product were computed using an error matrix approach in which 1754 independent validation samples were used. In addition, a comparison was performed with other available cropland maps as well as with LULC maps to show spatial similarity. Finally, the cropland area results derived from the map were compared with UN FAO statistics. The independent accuracy assessment showed a weighted overall accuracy of 94, with a producers accuracy of 85.9 (or omission error of 14.1), and users accuracy of 68.5 (commission error of 31.5) for the cropland class. The total net cropland area (TNCA) of Africa was estimated as 313 Mha for the nominal year 2015.

Cropland mapping; cropland areas; 30-m; Landsat-8;↗

U.S. Eastern Continental Shelf Carbon Cycling (USECoS): Modeling, Data Assimilation, and Analysis

Although the oceans play a major role in the uptake of fossil fuel CO2 from the atmosphere, there is much debate about the contribution from continental shelves, since many key shelf fluxes are not yet well quantified: the exchange of carbon across the land-ocean and shelf-slope interfaces, air-sea exchange of CO2, burial, and biological processes including productivity. Our goal is to quantify these carbon fluxes along the eastern U.S. coast using models quantitatively verified by comparison to observations, and to establish a framework for predicting how these fluxes may be modified as a result of climate and land use change. Our research questions build on those addressed with previous NASA funding for the USECoS (U.S. Eastern Continental Shelf Carbon Cycling) project. We have developed a coupled biogeochemical ocean circulation model configured for this study region and have extensively evaluated this model with both in situ and remotely-sensed data. Results indicate that to further reduce uncertainties in the shelf component of the global carbon cycle, future efforts must be directed towards 1) increasing the resolution of the physical model via nesting and 2) making refinements to the biogeochemical model and quantitatively evaluating these via the assimilation of biogeochemical data (in situ and remotely-sensed). These model improvements are essential for better understanding and reducing estimates of uncertainties in current and future carbon transformations and cycling in continental shelf systems. Our approach and science questions are particularly germane to the carbon cycle science goals of the NASA Earth Science Research Program as well as the U.S. Climate Change Research Program and the North American Carbon Program. Our interdisciplinary research team consists of scientists who have expertise in the physics and biogeochemistry of the U.S. eastern continental shelf, remote-sensing data analysis and data assimilative numerical models.

Mannino, Antonio↗

Estimating Dust and Water Ice Content of the Martian Atmosphere From THEMIS Data

Researchers at JPL and Arizona State University conducted a comparative study of three candidate algorithms for estimating components of the Martian atmosphere, using raw (uncalibrated) data collected by the Thermal Emission Imaging System (THEMIS). THEMIS is an instrument onboard the Mars Odyssey spacecraft that acquires image data in five visible and nine infrared (IR) wavelength bands. The algorithms under study used data collected from eight of the nine IR bands to estimate the dust and water ice content of the atmosphere. Such an algorithm could be used in onboard data processing to trigger other algorithms that search for features of scientific interest and to reduce the volume of data transmitted to Earth. The algorithms studied were based on regression models. In the study, the optical depths estimated by these algorithms were compared with optical depths estimated in ground-based processing using fully calibrated data from both THEMIS and the Thermal Emission Spectrometer (TES). TES is an instrument onboard the Mars Global Surveyor spacecraft that also observes the planet at infrared wavelengths, but at a lower spatial resolution than THEMIS does. Of the algorithms studied, the one that performed best was based on a Gaussian Support Vector Machine regression model. The test results indicated that this algorithm, operating on the raw data, had error rates that were within the uncertainty associated with the estimates obtained by the groundbased analysis of the fully calibrated data. This level of fidelity demonstrates that these algorithms are sufficiently accurate for use in an onboard setting.

Bandfield, Joshua↗

P21C-2113: Constraining Hesperian Martian PCO2 from Mineral Analysis at Gale Crater

Carbon dioxide is an essential atmospheric component in martian climate models that attempt to reconcile a faint young sun with widespread evidence of liquid water at the planet's surface in the Noachian and Early Hesperian. Current estimates of ancient martian CO levels, derived from global inventories of carbon, and orbital detections of Noachian and Early Hesperian clay mineralbearing terrains indicate CO levels that are unable to support warm and wet conditions. These estimates are subject to various sources of uncertainty however. Mineral and contextual sedimentary environmental data collected by the Mars Science Laboratory rover Curiosity in Gale Crater provide a more direct means of estimating the atmospheric partial pressure of CO (P ) coinciding with a long-lived lake system in Gale crater at approximately 3.5 Ga. Results from a reaction transport model, which simulates mineralogy observed within the Sheepbed member at Yellowknife Bay by coupling mineral equilibria with carbonate precipitation kinetics and rates of sedimentation, indicate atmospheric levels in the 10's mbar range. At such low P levels, climate models are unable to warm Hesperian Mars anywhere near the freezing point of water and other gases are required to raise atmospheric pressure to prevent lakes from boiling away. Thus, lacustrine features of Gale formed in a cold environment by a mechanism yet to be determined, or the climate models lack an essential component that would serve to elevate surface temperatures, at least temporally and/or locally, on Hesperian Mars. Our results also impose restrictions on the potential role of atmospheric CO in inferred warmer conditions of the Noachian.

Bristow, Thomas↗

Constraining Warm Dark Matter and Population III Stars with the Global 21 cm Signal

Upcoming ground- and space-based experiments may have sufficient accuracy to place significant constraints upon high-redshift star formation, Reionization, and dark matter (DM) using the global 21 cm signal of the intergalactic medium. In the early universe, when the relative abundance of low-mass DM halos was important, measuring the global signal would place constraints on the damping of structure formation caused by DM having a higher relic velocity (warm dark matter, or WDM) than in cold dark matter. Such damping, however, can be mimicked by altering the star formation efficiency (SFE) and can be difficult to detect because of the presence of Population III stars with unknown properties. We study these various cases and their degeneracies with the WDM mass parameter mX using a Fisher matrix analysis. We study the mX = 7 keV case and a star formation model that parameterizes the SFE as a strong function of halo mass and include several variations of this model along with three different input noise levels for the likelihood; we also use a minimum halo virial temperature for collapse near the molecular cooling threshold. We find that when the likelihood includes only Population II stars, mX is constrained to an uncertainty of ∼0.4 keV for all models and noise levels at the 68% confidence level. When the likelihood includes weak Population III stars, mX ∼ 0.3 keV, and if Population III star formation is relatively efficient, mX ∼ 0.1 keV uncertainty, with tight Population III star formation parameter constraints. Our results show that the global 21 cm signal is a promising test-bed for WDM models, even in the presence of strong degeneracies with astrophysical parameters.

Cosmology↗