Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Uncertainty”

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 937 records · Page 52

THE MEASUREMENT, TREATMENT, AND IMPACT OF SPECTRAL COVARIANCE AND BAYESIAN PRIORS IN INTEGRAL -FIELD SPECTROSCOPY OF EXOPLANETS

The recovery of an exoplanet’s atmospheric parameters from its spectrum requires accurate knowledge of the spectral errors and covariances. Unfortunately, the complex image processing used in high-contrast integral-field spectrograph (IFS) observations generally produces spectral covariances that are poorly understood and often ignored. In this work, we show how to measure the spectral errors and covariances and include them self-consistently in parameter retrievals. By combining model exoplanet spectra with a realistic noise model generated from the Gemini Planet Imager (GPI) early science data, we show that ignoring spectral covariance in high-contrast IFS data can both bias inferred parameters and lead to unreliable confidence regions on those parameters. This problem is made worse by the common practice of scaling the χ {sup 2} per degree of freedom to unity; the input parameters then fall outside the 95% confidence regions in as many as ∼80% of noise realizations. The biases we observe can approach the typical levels of precision achieved in high-contrast spectroscopy. Accounting for realistic priors in fully Bayesian retrievals can also have a significant impact on the inferred parameters. Plausible priors on effective temperature and surface gravity can vary by an order of magnitude across the confidence regions appropriate for objects with weak age constraints; priors for objects with good age constraints are dominated by modeling uncertainties. Our methods are directly applicable to existing high-contrast IFSs including GPI and SPHERE, as well as upcoming instruments like CHARIS and, ultimately, WFIRST-AFTA.

79 ASTRONOMY AND ASTROPHYSICS↗

The IRX–β Relation: Insights from Simulations

We study the relationship between the UV continuum slope and infrared excess (IRX≡L{sub IR}/L{sub FUV}) predicted by performing dust radiative transfer on a suite of hydrodynamical simulations of galaxies. Our suite includes both isolated disk galaxies and mergers intended to be representative of galaxies at both z∼0 and z∼2--3. Our low-redshift systems often populate a region around the locally calibrated Meurer et al. relation but move above the relation during merger-induced starbursts. Our high-redshift systems are blue and IR luminous and therefore lie above the Meurer et al. relation. The value of β strongly depends on the dust type used in the RT simulation: Milky-Way-type dust leads to significantly more negative (bluer) slopes compared with Small-Magellanic-Cloud-type dust. The effect on β due to variations in the dust composition with galaxy properties or redshift is the dominant model uncertainty. The dispersion in β is anticorrelated with specific star formation rate (sSFR) and tends to be higher for the z∼2--3 simulations. In the actively star-forming z∼2--3 simulated galaxies, dust attenuation dominates the dispersion in β, whereas in the z∼0 simulations, the contributions of star formation history (SFH) variations and dust are similar. For low-sSFR systems at both redshifts, SFH variations dominate the dispersion. Finally, the simulated z∼2--3 isolated disks and mergers both occupy a region in the IRX--β plane consistent with observed z∼2--3 dusty star-forming galaxies (DSFGs). Thus, contrary to some claims in the literature, the blue colors of high-z DSFGs do not imply that they are short-lived starbursts.

79 ASTRONOMY AND ASTROPHYSICS↗

Precise Time Delays from Strongly Gravitationally Lensed Type Ia Supernovae with Chromatically Microlensed Images

Time delays between the multiple images of strongly gravitationally lensed Type Ia supernovae (glSNe Ia) have the potential to deliver precise cosmological constraints, but the effects of microlensing on time delay extraction have not been studied in detail. Here we quantify the effect of microlensing on the glSN Ia yield of the Large Synoptic Survey Telescope (LSST) and the effect of microlensing on the precision and accuracy of time delays that can be extracted from LSST glSNe Ia. Microlensing has a negligible effect on the LSST glSN Ia yield, but it can be increased by a factor of ~2 over previous predictions to 930 systems using a novel photometric identification technique based on spectral template fitting. Crucially, the microlensing of glSNe Ia is achromatic until three rest-frame weeks after the explosion, making the early-time color curves microlensing-insensitive time delay indicators. By fitting simulated flux and color observations of microlensed glSNe Ia with their underlying, unlensed spectral templates, we forecast the distribution of absolute time delay error due to microlensing for LSST, which is unbiased at the sub-percent level and peaked at 1% for color curve observations in the achromatic phase, while for light-curve observations it is comparable to state-of-the-art mass modeling uncertainties (4%). About 70% of LSST glSN Ia images should be discovered during the achromatic phase, indicating that microlensing time delay uncertainties can be minimized if prompt multicolor follow-up observations are obtained. Lastly, accounting for microlensing, the 1-2 day time delay on the recently discovered glSN Ia iPTF16geu can be measured to 40% precision, limiting its cosmological utility.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Measurement of the Temperature–density Relation of the IGM from Voigt Profile Fitting

We decompose the Lyman-α (Lyα) forest of an extensive sample of 75 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-b cutoff in the b- distribution set primarily by the thermal state of the intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range 2.0 ≤ z ≤ 3.4, which allows us to measure the evolution of the IGM temperature-density (T = T 0 (ρ/ρ 0 ) γ-1 ) relation parameters T 0 and γ. We calibrate our measurements against mock Lyα forest data generated using 26 hydrodynamic simulations with different thermal histories from the THERMAL suite, also encompassing different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of T 0 (γ) shows a suggestive peak (dip) at z = 2.9 (z = 3). Our measured evolution of T 0 and γ is generally in good agreement with previous determinations in the literature. Both the peak in the evolution of T 0 at z = 2.8, as well as the high temperatures T 0 ≃ 15,000-20,000 K that we observe at 2.4 < z < 3.4, strongly suggest that a significant episode of heating occurred after the end of H I reionization, which was most likely the cosmic reionization of He II.

79 ASTRONOMY AND ASTROPHYSICS↗

The Challenges Ahead for Multimessenger Analyses of Gravitational Waves and Kilonova: A Case Study on GW190425

In recent years, there have been significant advances in multimessenger astronomy due to the discovery of the first, and so far only confirmed, gravitational wave event with a simultaneous electromagnetic (EM) counterpart, as well as improvements in numerical simulations, gravitational wave (GW) detectors, and transient astronomy. This has led to the exciting possibility of performing joint analyses of the GW and EM data, providing additional constraints on fundamental properties of the binary progenitor and merger remnant. Here, in this work, we present a new Bayesian framework that allows inference of these properties, while taking into account the systematic modeling uncertainties that arise when mapping from GW binary progenitor properties to photometric light curves. We extend the relative binning method presented in Zackay et al. to include extrinsic GW parameters for fast analysis of the GW signal. The focus of our EM framework is on light curves arising from r-process nucleosynthesis in the ejected material during and after merger, the so-called kilonova, and particularly on black hole–neutron star systems. As a case study, we examine the recent detection of GW190425, where the primary object is consistent with being either a black hole or a neutron star. We show quantitatively how improved mapping between binary progenitor and outflow properties, and/or an increase in EM data quantity and quality are required in order to break degeneracies in the fundamental source parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment

We present a detection of 21 cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment. Radio observations acquired over 102 nights are used to construct maps that are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasars (QSOs) from the eBOSS clustering catalogs. We find decisive evidence for a detection when stacking on all three tracers of LSS, with the logarithm of the Bayes factor equal to 18.9 (LRG), 10.8 (ELG), and 56.3 (QSO). An alternative frequentist interpretation, based on the likelihood ratio test, yields a detection significance of 7.1σ (LRG), 5.7σ (ELG), and 11.1σ (QSO). These are the first 21 cm intensity mapping measurements made with an interferometer. We constrain the effective clustering amplitude of neutral hydrogen (H I ), defined as ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}\equiv {10}^{3}\,{{\rm{\Omega }}}_{{\rm{H}}\,{\rm\small{I}}}\left({b}_{{\rm{H}}\,{\rm\small{I}}}+\langle \,f{\mu }^{2}\rangle \right)$, where Ω H I is the cosmic abundance of H I , b H I is the linear bias of H I , and $\langle$fμ 2 $\rangle$ = 0.552 encodes the effect of redshift-space distortions at linear order. We find ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={1.51}_{-0.97}^{+3.60}$ for LRGs (z = 0.84), ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={6.76}_{-3.79}^{+9.04}$ for ELGs (z = 0.96), and ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={1.68}_{-0.67}^{+1.10}$ for QSOs (z = 1.20), with constraints limited by modeling uncertainties at nonlinear scales. We are also sensitive to bias in the spectroscopic redshifts of each tracer, and we find a nonzero bias Δ v = - 66 ± 20 km s -1 for the QSOs. We split the QSO catalog into three redshift bins and have a decisive detection in each, with the upper bin at z = 1.30 producing the highest-redshift 21 cm intensity mapping measurement thus far.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of Polarized Dust Emission Using Data from the First Flight of SPIDER

Using data from the first flight of Spider and from the Planck High Frequency Instrument, we probe the properties of polarized emission from interstellar dust in the Spider observing region. Component-separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses of diffuse Galactic dust emission spanning the full Spider region demonstrate (i) a spectral energy distribution that is broadly consistent with a modified-blackbody (MBB) model with a spectral index of β d = 1.45 ± 0.05 (1.47 ± 0.06) for E (B)-mode polarization, slightly lower than that reported by Planck for the full sky; (ii) an angular power spectrum broadly consistent with a power law; and (iii) no significant detection of line-of-sight polarization decorrelation. Tests of several modeling uncertainties find only a modest impact (~10% in σ r ) on Spider's sensitivity to the cosmological tensor-to-scalar ratio. The size of the Spider region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature T d = 19.6 K, an analysis of two independent subregions of that field results in inferred values of β d = 1.52 ± 0.06 and β d = 1.09 ± 0.09, which are inconsistent at the 3.9σ level. Furthermore, a joint analysis of Spider and Planck 217 and 353 GHz data within one subregion is inconsistent with a simple MBB at more than 3σ, assuming a common morphology of polarized dust emission over the full range of frequencies. This evidence of variation may inform the component-separation approaches of future cosmic microwave background polarization experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

The Impact of White Dwarf Luminosity Profiles on Oscillation Frequencies

KIC 08626021 is a pulsating DB white dwarf (WD) of considerable recent interest, and the first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based WD models choose a luminosity profile where the luminosity is proportional to the enclosed mass, ${L}_{r}\,\propto \,{M}_{r}$, independent of the effective temperature T eff . Evolutionary models of young WDs with T eff ≳ 25,000 K suggest that neutrino emission gives rise to luminosity profiles with L r ${/}\!\!\!\!\!\!{\propto }$ M r . We explore this contrast by comparing the oscillation frequencies between two nearly identical WD models: one with an enforced ${L}_{r}\propto {M}_{r}$ luminosity profile, and the other with a luminosity profile determined by the star's previous evolution history. We find that the low-order g-mode frequencies differ by up to ≃70 μHz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an ${L}_{r}\propto {M}_{r}$ profile may have uncorrected, effectively random errors at the level of tens of μHz. A mean frequency difference of 30 μHz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of ≃12% in WD mass, ≃9% in the radius, and ≃3% in the central oxygen mass fraction.

79 ASTRONOMY AND ASTROPHYSICS↗

Assimilating partial observation to enhance feedback control of stochastic dynamical systems

Here, in this paper, we present a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation, which is often not feasible in real-world fluid dynamics control problems. Our approach underscores the significance of utilizing observational data to inform control policy design. Specifically, we introduce a kernel learning backward stochastic differential equation (SDE) filter to enhance data assimilation efficiency and propose a sample-wise stochastic optimization method within the stochastic maximum principle framework. We demonstrate the efficacy and accuracy of our method in the control of advection-diffusion-reaction flow problem and the Dubins airplane maneuvering problem with model uncertainty.

data driven↗

Important role of stratospheric injection height for the distribution and radiative forcing of smoke aerosol from the 2019–2020 Australian wildfires

More than 1 Tg smoke aerosol was emitted into the atmosphere by the exceptional 2019–2020 southeastern Australian wildfires. Triggered by the extreme fire heat, several deep pyroconvective events carried the smoke directly into the stratosphere. Once there, smoke aerosol remained airborne considerably longer than in lower atmospheric layers. The thick plumes traveled eastward, thereby being distributed across the high and mid-latitudes in the Southern Hemisphere, enhancing the atmospheric opacity. Due to the increased atmospheric lifetime of the smoke plume, its radiative effect increased compared to smoke that remains in lower altitudes. Global models describing aerosol-climate impacts lack adequate descriptions of the emission height of aerosols from intense wildfires. Here, we demonstrate, by a combination of aerosol-climate modeling and lidar observations, the importance of the representation of those high-altitude fire smoke layers for estimating the atmospheric energy budget. Through observation-based input into the simulations, the Australian wildfire emissions by pyroconvection are explicitly prescribed to the lower stratosphere in different scenarios. Based on our simulations, the 2019–2020 Australian fires caused a significant top-of-atmosphere (TOA) hemispheric instantaneous direct radiative forcing signal that reached a magnitude comparable to the radiative forcing induced by anthropogenic absorbing aerosol. Up to +0.50 W m -2 instantaneous direct radiative forcing was modeled at TOA, averaged for the Southern Hemisphere (+0.25 W m -2 globally) from January to March 2020 under all-sky conditions. At the surface, on the other hand, an instantaneous solar radiative forcing of up to -0.81 W m -2 was found for clear-sky conditions, with the respective estimates depending on the model configuration and subject to the model uncertainties in the smoke optical properties. Since extreme wildfires are expected to occur more frequently in the rapidly changing climate, our findings suggest that high-altitude wildfire plumes must be adequately considered in climate projections in order to obtain reasonable estimates of atmospheric energy budget changes.

54 ENVIRONMENTAL SCIENCES↗

Indirect climate impacts of the Hunga eruption

Injection of sulfur and water vapour by the Hunga volcanic eruption significantly altered chemical composition and radiative budget of the stratosphere. Yet, whether the eruption could also affect surface climate, especially via indirect pathways, remains poorly understood. Here we investigate these effects using large ensembles of simulations with the CESM2(WACCM6) Earth system model, incorporating interactive chemistry and aerosols in both coupled ocean and atmosphere-only configurations. We find some statistically significant extratropical regional climate responses to the eruption driven by circulation changes; these are partially linked to the modulation of El Niño–Southern Oscillation, and its associated teleconnections, and to perturbations of the stratospheric polar vortex in both hemispheres. The stratospheric anomalies affect surface climate through modulating the North Atlantic Oscillation in the Northern Hemisphere (up to three boreal winters following the eruption) and the Southern Annular Mode in the Southern Hemisphere in late 2023. The latter is partly related to a concurrent reduction in Antarctic ozone, as increased stratospheric aerosols and water vapor reach the polar vortex. Our study suggests that the eruption could have had a non-negligible influence on regional surface climate, and discusses the mechanisms via which such an influence could occur. However, the results also highlight that this forcing is relatively weak compared to interannual variability, and is subject to model uncertainties in the representation of key processes. More research is thus needed before definitive statements on the role of the eruption in contributing to surface climate and weather events in the following years are made.

Geosciences↗

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗

A new test statistic for climate models that includes field and spatial dependencies using Gaussian Markov random fields

Abstract. A new test statistic for climate model evaluation has been developed that potentially mitigates some of the limitations that exist for observing and representing field and space dependencies of climate phenomena. Traditionally such dependencies have been ignored when climate models have been evaluated against observational data, which makes it difficult to assess whether any given model is simulating observed climate for the right reasons. The new statistic uses Gaussian Markov random fields for estimating field and space dependencies within a first-order grid point neighborhood structure. We illustrate the ability of Gaussian Markov random fields to represent empirical estimates of field and space covariances using "witch hat" graphs. We further use the new statistic to evaluate the tropical response of a climate model (CAM3.1) to changes in two parameters important to its representation of cloud and precipitation physics. Overall, the inclusion of dependency information did not alter significantly the recognition of those regions of parameter space that best approximated observations. However, there were some qualitative differences in the shape of the response surface that suggest how such a measure could affect estimates of model uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Hydrologic calibration of paired watersheds using a MOSUM approach

Paired watershed studies have historically been used to quantify hydrologic effects of land use and management practices by concurrently monitoring two neighboring watersheds (a control and a treatment) during the calibration (pre-treatment) and post-treatment periods. This study characterizes seasonal water table and flow response to rainfall during the calibration period and tests a change detection technique of moving sums of recursive residuals (MOSUM) to select calibration periods for each control-treatment watershed pair when the regression coefficients for daily water table elevation (WTE) were most stable to reduce regression model uncertainty. The control and treatment watersheds included 1–3 year intensively managed loblolly pine ( Pinus taeda L. ) with natural understory, same age loblolly pine intercropped with switchgrass ( Panicum virgatum ), 14–15 year thinned loblolly pine with natural understory (control), and switchgrass only. Although monitoring during the calibration period spanned 2009 to 2012, silvicultural operational practices that occurred during this period such as harvesting of existing stand and site preparation for pine and switchgrass establishment may have acted as external factors, potentially shifting hydrologic calibration relationships between control and treatment watersheds. Results indicated that MOSUM was able to detect significant changes in regression parameters for WTE due to silvicultural operations. This approach also minimized uncertainty of calibration relationships which could otherwise mask marginal treatment effects. All calibration relationships developed using this MOSUM method were quantifiable, strong, and consistent with Nash–Sutcliffe Efficiency (NSE) greater than 0.97 for WTE and NSE greater than 0.92 for daily flow, indicating its applicability for choosing calibration periods of paired watershed studies.

54 ENVIRONMENTAL SCIENCES↗

Lidar measurements of yawed-wind-turbine wakes: characterization and validation of analytical models

Abstract. Wake measurements of a scanning Doppler lidar mounted on the nacelle of a full-scale wind turbine during a wake-steering experiment were used for the characterization of the wake flow, the evaluation of the wake-steering set-up, and the validation of analytical wake models. Inflow-scanning Doppler lidars, a meteorological mast, and the supervisory control and data acquisition (SCADA) system of the wind turbine complemented the set-up. Results from the wake-scanning Doppler lidar showed an increase in the wake deflection with the yaw angle and that the wake deflection was not in all cases beneficial for the power output of a downstream turbine due to a bias of the inflow wind direction perceived by the yawed wind turbine and the wake-steering design implemented. Both observations could be reproduced with an analytical model that was initialized with the inflow measurements. Error propagation from the inflow measurements that were used as model input and the power coefficient of a waked wind turbine contributed significantly to the model uncertainty. Lastly, the span-wise cross section of the wake was strongly affected by wind veer, masking the effects of the yawed wind turbine on the wake cross sections.

17 WIND ENERGY↗

MSD CoP Webinar: Applied Science for Decision-Making at the Water-Energy Nexus

Context: This webinar was hosted by the MultiSector Dynamics Community of Practice (MSD CoP; https://multisectordynamics.org). Abstract: The energy transition is a game-changer across multiple sectors, requiring researchers and practitioners to re-evaluate the evolving connections between climate, water, and power systems, as well as the associated co-management of resources and decision-making. In anticipation of new technologies, policies, business models, and other solutions, applied science for decision-making at the climate-water-energy nexus is needed to support transitions with the expected climate resilience, maintained and enhanced energy security, thriving economies, and equity considerations. Through panelists' applied research and operational experience, the webinar will provide lesson-learnt in how to increase the readiness of water-energy research and become actionable in the context of energy transitions. A number of themes will be discussed, including the connection between growing computational resources and complexity of models, and the implications on the intended users, the actual actionable products, the availability of data for validation of the models, uncertainty characterization and decision-making under uncertainty, and the communication of the novelty to decision-makers and the public. Presenters : David McCollum (Oak Ridge National Laboratory; Co-Chair), and Gokul Iyer (Pacific Northwest National Laboratory; Co-Chair), Curt Jawdy (Tennessee Valley Authority), Nathalie Voisin (Pacific Northwest National Laboratory), and Andrew D. Jones (Lawrence Berkeley National Laboratory) Moderator: Pat M. Reed (MSD CoP Facilitation Team) This webinar was held on: April 30, 2024 from 1-2 PM ET

Energy↗

Tests and comparisons of gravity models.

Optical observations of the GEOS satellites were used to obtain orbital solutions with different sets of geopotential coefficients. The solutions were compared before and after modification to high order terms (necessary because of resonance) and were then analyzed by comparing subsequent observations with predicted trajectories. The most important source of error in orbit determination and prediction for the GEOS satellites is the effect of resonance found in most published sets of geopotential coefficients. Modifications to the sets yield greatly improved orbits in most cases. The results of these comparisons suggest that with the best optical tracking systems and gravity models, satellite position error due to gravity model uncertainty can reach 50-100 m during a heavily observed 5-6 day orbital arc. If resonant coefficients are estimated, the uncertainty is reduced considerably.

Marsh, J. G.↗

Mars lander position estimation in the presence of ephemeris biases.

The process of estimating the location of a spacecraft landed on the surface of Mars is investigated through the application of statistical estimation techniques to earth-based radio tracking data. The spacecraft location and the tracking geometry and schedule are consistent with Viking-type mission constraints. With mission control requirements in mind, the investigation is restricted to analysis of a short data arc (approximately 3 days). Statistics of the spacecraft location are obtained through analysis of (direct-link) tracking data for the landed spacecraft and through simultaneous analysis of tracking data for both a landed and an orbiting spacecraft. These estimates include the effects of model uncertainties in the ephemeris of Mars, tracking station locations, the Mars rotational period, the Mars gravity field, and the orientation of Mars axis of rotation. The most significant of these effects is shown to be due to the Mars ephemeris uncertainty. A dual spacecraft tracking technique is presented for substantially reducing these ephemeris effects.

Blackshear, W. T.↗