Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust statistics”

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 163 records · Page 9

Four Years of Type Ia Supernovae Observed by TESS: Early-time Light-curve Shapes and Constraints on Companion Interaction Models

Abstract We present 307 type Ia supernova (SN) light curves from the first 4 yr of the Transiting Exoplanet Survey Satellite mission. We use this sample to characterize the shapes of the early-time light curves, measure the rise times from first light to peak, and search for companion star interactions. Using simulations, we show that light curves must have noise <10% of the peak flux to avoid biases in the early-time light-curve shape, restricting our quantitative analysis to 74 light curves. We find that the mean power-law index t β 1 of the early-time light curves isβ 1 = 1.93 ± 0.57, and the mean rise time to peak is 15.7 ± 3.5 days. The underlying population distribution forβ 1 may instead consist of a Gaussian component with mean 2.29, width 0.34, and a long tail extending to values less than 1.0. We find that the data can rarely distinguish between models with and without companion interaction models. Nevertheless, we find three high-quality light curves that tentatively prefer the addition of a companion interaction model, but the statistical evidence for the companion interactions is not robust. We also find two SNe that disfavor the addition of a companion interaction model to a curved power-law model. Taking the 74 SNe together, we calculate 3σupper limits on the presence of companion signatures to control for orientation effects that can hide companions in individual light curves. Our results rule out common progenitor systems with companions having Roche lobe radii >31R ⊙ (separations >5.7 × 10 12 cm, 99.9% confidence level) and disfavor companions having Roche lobe radii >10R ⊙ (separations >1.9 × 10 12 cm, 95% confidence level). Lastly, we discuss the implications of our results for the intrinsic fraction of single degenerate progenitor systems.

Astronomy & Astrophysics↗

Synthesis of ARM User Facility Surface Rainfall Datasets to Construct a Best Estimate Value Added Product (PrecipBE)

Surface precipitation measurements are essential for Earth system model (ESM) evaluation and understanding cloud processes. An ever-growing need for robust, temporally evolving, and easy-to-use statistical datasets provides motivation for a baseline ground-based precipitation properties data product. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility operates an extensive suite of precipitation instruments with various sensitivities and operating mechanisms, which render the decision of which instrument to use based on one or more fixed thresholds challenging and prone to errors and bias. Using a long-term instrument inter-comparison from a unique per-precipitation event perspective, rather than instantaneous sample comparison, we demonstrate that ARM rainfall-measuring instruments are generally consistent with each other at the statistical level. Inter-instrument deviations at the single event level can be large, especially for specific rainfall event properties such as maximum precipitation rates. A machine-learning (ML) analysis using a random forest regressor indicates that in some cases, depending on instrument, local site climatology, and/or specific deployment configuration, certain atmospheric state variables influence the measured quantities in an unpredictable manner. Thus, a-priori weighting of different instruments does not necessarily lead to more accurate and less biased synthesis of instrument data. These results motivate the design of the ARM precipitation best-estimate (PrecipBE) value-added product, which incorporates all valid precipitation data while considering data quality and other instrument limitations. PrecipBE consists of time series and tabular statistics datasets in an easy-to-use and insightful per-precipitation event format. It provides a large set of precipitation event properties supplemented with ancillary data from ARM datasets that correspond to the detected precipitation events. We describe the PrecipBE algorithm and demonstrate its use via the examination of a single-day output as well as a long-term trend analysis of precipitation events at the ARM Southern Great Plains (SGP) site, covering more than 30 years of data. The trend analysis tentatively suggests a long-term temporal tendency for mainly shorter and less intense precipitation events at the SGP site, but a long-term increase in annual rainfall by more than 36 mm (5 %) per decade. This rainfall trend is catalyzed primarily by more extreme event properties of relatively rare, intense precipitation events, with event total and 1 min maximum precipitation rate at a 1 year timeframe increasing up to 5 mm and 9 mm h −1 (several percent) per decade, respectively. While the currently available PrecipBE datasets (at https://adc.arm.gov/discovery/, last access: 8 December 2025) cover rainfall from multiple ARM deployments up to March 2025, PrecipBE is planned to be expanded to include solid-phase precipitation and will soon become an operational product with a several-day lag from real-time. We invite the ARM user community to leverage this new product and welcome user feedback to further enhance the dataset.

Silber, Israel [Pacific Northwest National Laborat↗

Stratospheric Assimilation of Chemical Tracer Observations Using a Kalman Filter: Chi-Square Validated Results and Analysis of Variance and Correlation Dynamics - Pt. 2

A Kalman filter system designed for the assimilation of limb-sounding observations of stratospheric chemical tracers, which has four tunable covariance parameters, was developed in Part I (Menard et al. 1998) The assimilation results of CH4 observations from the Cryogenic Limb Array Etalon Sounder instrument (CLAES) and the Halogen Observation Experiment instrument (HALOE) on board of the Upper Atmosphere Research Satellite are described in this paper. A robust (chi)(sup 2) criterion, which provides a statistical validation of the forecast and observational error covariances, was used to estimate the tunable variance parameters of the system. In particular, an estimate of the model error variance was obtained. The effect of model error on the forecast error variance became critical after only three days of assimilation of CLAES observations, although it took 14 days of forecast to double the initial error variance. We further found that the model error due to numerical discretization as arising in the standard Kalman filter algorithm, is comparable in size to the physical model error due to wind and transport modeling errors together. Separate assimilations of CLAES and HALOE observations were compared to validate the state estimate away from the observed locations. A wave-breaking event that took place several thousands of kilometers away from the HALOE observation locations was well captured by the Kalman filter due to highly anisotropic forecast error correlations. The forecast error correlation in the assimilation of the CLAES observations was found to have a structure similar to that in pure forecast mode except for smaller length scales. Finally, we have conducted an analysis of the variance and correlation dynamics to determine their relative importance in chemical tracer assimilation problems. Results show that the optimality of a tracer assimilation system depends, for the most part, on having flow-dependent error correlation rather than on evolving the error variance.

Menard, Richard↗

Contemporary Trends in the Acid-Base Status of the Two Acid-Sensitive Streams in Western Maryland

Recovery of streamwater acid neutralizing capacity (ANC) resulting from declines in regional acid deposition was examined using contemporary (1990-2005) data from two long-term monitoring stations located on the Appalachian Plateau in western Maryland, U.S. Two computational methods were used to estimate daily, monthly, and annual fluxes and discharge-weighted concentrations of ANC, sulfate, nitrate, and base cations over the period of record, and two statistical methods were used to evaluate long-term trends in fluxes and concentrations. The methods used to estimate concentrations, as well as the, statistical techniques, produced very similar results, underlining the robustness of the identified trends. We found clear evidence that streamwater sulfate concentrations have declined at an average rate of about 3 (microns)eq L(exp -1) yr(exp -1) at the two sites due to a 34% reduction in wet atmospheric sulfur deposition. Trends in nitrate concentrations appear to be related to other watershed factors, especially forest disturbance. The best evidence of recovery is based on a doubling of ANC (from 21 to 42 (microns)eq L(exp -1) at the more acid-sensitive site over the 16-year period. A slowing, or possible reversal, in the sulfate, nitrate, and SBC trends is evident in our data and may portend a decline in the rate of--or end to--further recovery.

Eshleman, Keith N.↗

Optimized O'Neill/Glaser Model for Human Population of Space and its Impact on Survival Probabilities

Two contemporary issues foretell a shift from our historical Earth based industrial economy and habitation to a solar system based society. The first is the limits to Earth's carrying capacity, that is the maximum number of people that the Earth can support before a catastrophic impact to the health of the planet and human species occurs. The simple example of carrying capacity is that of a bacterial colony in a Petri dish with a limited amount of nutrient. The colony experiences exponential population growth until the carrying capacity is reached after which catastrophic depopulation often results. Estimates of the Earth s carrying capacity vary between 14 and 40 billion people. Although at current population growth rates we may have over a century before we reach Earth s carrying limit our influence on climate and resources on the planetary scale is becoming scientifically established. The second issue is the exponential growth of knowledge and technological power. The exponential growth of technology interacts with the exponential growth of population in a manner that is unique to a highly intelligent species. Thus, the predicted consequences (world famines etc.) of the limits to growth have been largely avoided due to technological advances. However, at the mid twentieth century a critical coincidence occurred in these two trends humanity obtained the technological ability to extinguish life on the planetary scale (by nuclear, chemical, biological means) and attained the ability to expand human life beyond Earth. This paper examines an optimized O Neill/Glaser model (O Neill 1975; Curreri 2007; Detweiler and Curreri 2008) for the economic human population of space. Critical to this model is the utilization of extraterrestrial resources, solar power and spaced based labor. A simple statistical analysis is then performed which predicts the robustness of a single planet based technological society versus that of multiple world (independent habitats) society.

Curreri, Peter A.↗

Optimized O'Neill/Glaser Model for Human Population of Space and its Impact on Survival Probabilities

Two contemporary issues foretell a shift from our historical Earth based industrial economy and habitation to a solar system based society. The first is the limits to Earth s carrying capacity, that is the maximum number of people that the Earth can support before a catastrophic impact to the health of the planet and human species occurs. The simple example of carrying capacity is that of a bacterial colony in a Petri dish with a limited amount of nutrient. The colony experiences exponential population growth until the carrying capacity is reached after which catastrophic depopulation often results. Estimates of the Earth s carrying capacity vary between 14 and 40 billion people. Although at current population growth rates we may have over a century before we reach Earth s carrying limit our influence on climate and resources on the planetary scale is becoming scientifically established. The second issue is the exponential growth of knowledge and technological power. The exponential growth of technology interacts with the exponential growth of population in a manner that is unique to a highly intelligent species. Thus, the predicted consequences (world famines etc.) of the limits to growth have been largely avoided due to technological advances. However, at the mid twentieth century a critical coincidence occurred in these two trends humanity obtained the technological ability to extinguish life on the planetary scale (by nuclear, chemical, biological means) and attained the ability to expand human life beyond Earth. This paper examines an optimized O Neill/Glaser model (O Neill 1975; Curreri 2007; Detweiler and Curreri 2008) for the economic human population of space. Critical to this model is the utilization of extraterrestrial resources, solar power and spaced based labor. A simple statistical analysis is then performed which predicts the robustness of a single planet based technological society versus that of multiple world (independent habitats) society.

Curreri, Peter A.↗

Optimized O'Neill/Glaser Model for Human Population of Space and its Impact on Survival Probabilities

Two contemporary issues foretell a shift from our historical Earth based industrial economy and habitation to a solar system based society. The first is the limits to Earth s carrying capacity, that is the maximum number of people that the Earth can support before a catastrophic impact to the health of the planet and human species occurs. The simple example of carrying capacity is that of a bacterial colony in a Petri dish with a limited amount of nutrient. The colony experiences exponential population growth until the carrying capacity is reached after which catastrophic depopulation often results. Estimates of the Earth s carrying capacity vary between 14 and 40 billion people. Although at current population growth rates we may have over a century before we reach Earth s carrying limit our influence on climate and resources on the planetary scale is becoming scientifically established. The second issue is the exponential growth of knowledge and technological power. The exponential growth of technology interacts with the exponential growth of population in a manner that is unique to a highly intelligent species. Thus, the predicted consequences (world famines etc.) of the limits to growth have been largely avoided due to technological advances. However, at the mid twentieth century a critical coincidence occurred in these two trends humanity obtained the technological ability to extinguish life on the planetary scale (by nuclear, chemical, biological means) and attained the ability to expand human life beyond Earth. This paper examines an optimized O'Neill/Glaser model (O Neill 1975; Curreri 2007; Detweiler and Curreri 2008) for the economic human population of space. Critical to this model is the utilization of extraterrestrial resources, solar power and spaced based labor. A simple statistical analysis is then performed which predicts the robustness of a single planet based technological society versus that of multiple world (independent habitats) society.

Curreri, Peter A.↗

Fe3 Traffic Simulator Overview

NASA Ames developed a first of its kind air traffic simulation tool known as Flexible engine for Fast time evaluation of Flight environments (Fe3). The patent-pending technology is a robust simulation tool that provides the capability of statistically analyzing the high-density and low-altitude traffic system. With this capability, stakeholders can study the impacts of critical components (such as wind, surveillance, communication, collision, avoidance, traffic rules, energy consumption, etc.) in the low-altitude high-density traffic system, gain insights and help define requirements, policies, and protocols for a safe and efficient traffic system, and assess operational risks and optimize fight schedules.

Flexible engine for Fast time evaluation of Flight↗

Dark Energy Survey Year 6 Results: Galaxy-galaxy lensing

We present galaxy--galaxy lensing (GGL) measurements from the full six years of data from the Dark Energy Survey (DES Y6), covering $4031\,\mathrm{deg}^2$ and used in the DES Y6 $3\times2$pt cosmological analysis. We use the MagLim++ lens sample, containing $\sim 9$ million galaxies divided into six redshift bins, and the Metadetection source catalog, including $\sim 140$ million galaxies divided into four redshift bins. The mean tangential shear signal achieves a total signal-to-noise ratio (S/N) of $173$, corresponding to a $17\%$ improvement over DES Y3. After applying the scale cuts used in the cosmological analysis, with $R_{\min}=6\,\mathrm{Mpc}/h$ ($4\,\mathrm{Mpc}/h$) for the linear (nonlinear) galaxy-bias model, the S/N is reduced to $75$ (90). A comprehensive suite of validation tests demonstrates that the measurement is robust against observational and astrophysical systematics at the statistical precision required for the DES Y6 analysis. Although not used in the main cosmological analysis, we extract high--signal-to-noise geometric shear-ratio measurements from the galaxy--galaxy lensing signal on small angular scales. These measurements provide an internal consistency check on the photometric redshift distributions and shear calibration used in the $3\times2$pt analysis.

Giannini, G. [Barcelona, IEEC; Chicago U., KICP; B↗

Characterizing the Joint Effect of Diverse Test-Statistic Correlation Structures and Effect Size on False Discovery Rates in a Multiple-Comparison Study of Many Outcome Measures

In their 2009 Annals of Statistics paper, Gavrilov, Benjamini, and Sarkar report the results of a simulation assessing the robustness of their adaptive step-down procedure (GBS) for controlling the false discovery rate (FDR) when normally distributed test statistics are serially correlated. In this study we extend the investigation to the case of multiple comparisons involving correlated non-central t-statistics, in particular when several treatments or time periods are being compared to a control in a repeated-measures design with many dependent outcome measures. In addition, we consider several dependence structures other than serial correlation and illustrate how the FDR depends on the interaction between effect size and the type of correlation structure as indexed by Foerstner s distance metric from an identity. The relationship between the correlation matrix R of the original dependent variables and R~, the correlation matrix of associated t-statistics is also studied. In general R~ depends not only on R, but also on sample size and the signed effect sizes for the multiple comparisons.

Feiveson, Alan H.↗

Data and code for Daily and Multi-Day Extreme Rainfall Analysis Under Future Climates Using Stochastic Storm Transposition and NEX-GDDP-CMIP6 Over CONUS

This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.

54 ENVIRONMENTAL SCIENCES↗

Frequency domain identification for robust large space structure control design

A methodology is demonstrated for frequency domain identification of large space structures which systematically transforms experimental raw data into a form required for synthesizing H(infinity) controllers using modern robust control design software (e.g., Matlab Toolboxes). A unique feature of this approach is that the additive uncertainty is characterized to a specified statistic confidence rather than with hard bounds. In this study, the difference in robust performance is minimal between the two levels of confidence. In general cases, the present methodology provides a tool for performance/confidence level tradeoff studies. For simplicity, the additive uncertainty on a frequency grid is considered and the interpolation error in between grid points is neglected.

Yam, Y.↗

Statistical Decomposition and Machine Learning to Clean In Situ Spaceflight Magnetic Field Measurements

Robust in situ magnetic field measurements are critical to understanding the various mechanisms that couple mass, momentum, and energy throughout our solar system. However, the spacecraft on which magnetometers are often deployed contaminate the magnetic field measurements via onboard subsystems including reaction wheels and magnetorquers. Two magnetometers can be deployed at different distances from the spacecraft to determine an approximation of the interfering field for subsequent removal, but constant data streams from both magnetometers can be impractical due to power and telemetry limitations. Here we propose a method to identify and remove time-varying magnetic interference from sources such as reaction wheels using statistical decomposition and convolutional neural networks, providing high-fidelity magnetic field data even in cases where dual-sensor measurements are not constantly available. For example, a measurement interval from the Parker Solar Probe outboard magnetometer experienced a 95.1% reduction in reaction wheel interference following application of the proposed technique.

magnetic field data↗

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets such as sex or age of the model organism used. In the present study, NASA GeneLab-hosted RNAseq datasets from rodent liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC, to determine statistical differences between datasets before and after correction, Principal Component Analysis, to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the standard approach. Thus, the most robust standard correction will be implemented in the GeneLab Visualization 2.0 platform when datasets are combined.

GeneLab, RNA-seq, Batch Correction↗

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the respective standard approach. Of the methods tested, standard ComBat and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

GeneLab↗

Combining RNA-SEQ Datasets from NASA GENELAB: An Evaluation of Correction Methods

Background: Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. Methods: In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, the median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. Results: The results showed that the reference-based approach introduced several additional (and likely artificial) differentially expressed genes when compared with the respective standard approach. Conclusions: Of the methods tested, standard ComBat_seq and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

Finsam Samson↗

Vacuum Neutral Transport Model in UEDGE for Tokamak Far Scrape‐Off Layer

A model for neutral transport in the far scrape-off layer (SOL) vacuum region (vacuum neutral model) has been developed and implemented in UEDGE. Free-streaming neutral trajectories between the outermost UEDGE boundary and the vessel wall are preprocessed using DEGAS2 to construct a tele-transport matrix that captures non-local neutral relocation through the vacuum. Here, this matrix is then used in UEDGE as a non-local boundary condition for the neutral equations, preserving the robustness and convergence of the implicit solver without introducing statistical noise. Simulations of a DIII-D lower single-null configuration show that the vacuum neutral model relocates neutrals from the divertor to the upstream region, increasing the outer midplane separatrix density required for detachment onset by about 30%. However, the characteristic target temperature at detachment (T e, osp ~ 3 - 4eV) and the radiation front behavior remain unchanged.

DEGAS2↗

Mars Phoenix Entry, Descent, and Landing Simulation Design and Modelling Analysis

The 2007 Mars Phoenix Lander was launched in August of 2007 on a ten month cruise to reach the northern plains of Mars in May 2008. Its mission continues NASA s pursuit to find evidence of water on Mars. Phoenix carries upon it a slew of science instruments to study soil and ice samples from the northern region of the planet, an area previously undiscovered by robotic landers. In order for these science instruments to be useful, it was necessary for Phoenix to perform a safe entry, descent, and landing (EDL) onto the surface of Mars. The EDL design was defined through simulation and analysis of the various phases of the descent. An overview of the simulation and various models developed to characterize the EDL performance is provided. Monte Carlo statistical analysis was performed to assess the performance and robustness of the Phoenix EDL system and are presented in this paper. Using these simulation and modelling tools throughout the design and into the operations phase, the Mars Phoenix EDL was a success on May 25, 2008.

Prince, Jill L.↗