Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calibration bias”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper describes the attitude ground system (AGS) design to be used for support of the Magnetospheric MultiScale (MMS) mission. The AGS exists as one component of the mission operations control center. It has responsibility for validating the onboard attitude and accelerometer bias estimates, calibrating the attitude sensors and the spacecraft inertia tensor, and generating a definitive attitude history for use by the science teams. NASA's Goddard Space Flight Center (GSFC) in Greenbelt, Maryland is responsible for developing the MMS spacecraft, for the overall management of the MMS mission, and for mission operations. MMS is scheduled for launch in 2014 for a planned two-year mission. The MMS mission consists of four identical spacecraft flying in a tetrahedral formation in an eccentric Earth orbit. The relatively tight formation, ranging from 10 to 400 km, will provide coordinated observations giving insight into small-scale magnetic field reconnection processes. By varying the size of the tetrahedron and the orbital semi-major axis and eccentricity, and making use of the changing solar phase, this geometry allows for the study of both bow shock and magnetotail plasma physics, including acceleration, reconnection, and turbulence. The mission divides into two phases for science; these phases will have orbit dimensions of 1.2 x 12 Earth radii in the first phase and 1.2x25 Earth radii in the second in order to study the dayside magnetopause and the nightside magnetotail, respectively. The orbital periods are roughly one day and three days for the two mission phases. Each of the four MMS spacecraft will be spin stabilized at 3 revolutions per minute (rpm), with the spin axis oriented near the ecliptic north pole but tipped approximately 2.5 deg towards the Sun line. The main body of each spacecraft will be an eight-sided platform with diameter of 3.4 m and height of 1.2 m. Several booms are attached to this central core: two axial booms of 14.9 m length, two radial magnetometer booms of 5 m length, and four radial wire booms of 60 m length. Attitude and orbit control will use a set of axial and radial thrusters. A four-head star tracker and a slit-type digital Sun sensor (DSS) provide input for attitude determination. In addition, an accelerometer will be used for closed-loop orbit maneuver control. The primary AGS product will be a daily definitive attitude history. Due to power limitations, the star tracker and accelerometer data will not be available at all times. However, tracker data from at least 10 percent of each orbit and continuous DSS data will be provided. An extended Kalman filter (EKF) will be used to estimate the three-axis attitude (i.e., spin axis orientation and spin phase) and rotation rate for all times when the tracker data is valid. For other times, the attitude is generated by assuming a constant angular momentum vector in the inertial frame. The DSS sun pulse will provide a timing signal to maintain an accurate spin phase. There will be times when the Sun is occulted and DSS data is not available. If this occurs at the start or end of a definitive attitude product, then the spin phase will be extrapolated using the mean rate determined by the EKF.

Sedlak, Joseph E.↗
NASA NTRS2011

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper presents an overview of the attitude ground system (AGS) currently under development for the Magnetospheric Multiscale (MMS) mission. The primary responsibilities for the MMS AGS are definitive attitude determination, validation of the onboard attitude filter, and computation of certain parameters needed to improve maneuver performance. For these purposes, the ground support utilities include attitude and rate estimation for validation of the onboard estimates, sensor calibration, inertia tensor calibration, accelerometer bias estimation, center of mass estimation, and production of a definitive attitude history for use by the science teams. Much of the AGS functionality already exists in utilities used at NASA's Goddard Space Flight Center with support heritage from many other missions, but new utilities are being created specifically for the MMS mission, such as for the inertia tensor, accelerometer bias, and center of mass estimation. Algorithms and test results for all the major AGS subsystems are presented here.

Sedlak, Joseph E.↗
NASA NTRS2012

Cassini Thruster Calibration Algorithm Using Reaction Wheel Biasing Data

Thrust force estimates for the reaction control thrusters on-board Cassini spacecraft are presented in this paper. Cassini consists of two thruster branches (A and B) each with eight thrusters. The four Z-thrusters control the X and Y-axes, while the four Y-thrusters control the Z-axis. It is important to track the thrust force estimates in order to detect any thruster degradation and for supporting various activities in spacecraft operations (Titan flyby, spacecraft maneuvers). The Euler equation, which describes the rotational motion of the spacecraft during a reaction wheel bias event, is used to develop the algorithm. The thrust estimates are obtained from the pseudo inverse solution using flight telemetry during the bias. Results show that the A-branch Z3A and Z4A thrusters exhibited degraded thrust in November 2008. Due to the degraded thrust performance of Z3A and Z4A, A-branch usage was discontinued and prime branch was swapped to B-branch in March 2009. The thrust estimates from the B-branch do not show any degradation to date. The algorithm is used to trend the B-branch thrust force estimates as the mission continues.

attitude control↗
NASA NTRS2011

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper describes the attitude ground system (AGS) design to be used for support of the Magnetospheric MultiScale (MMS) mission. The AGS exists as one component of the mission operations control center. It has responsibility for validating the onboard attitude and accelerometer bias estimates, calibrating the attitude sensors and the spacecraft inertia tensor, and generating a definitive attitude history for use by the science teams. NASA's Goddard Space Flight Center (GSFC) in Greenbelt, Maryland is responsible for developing the MMS spacecraft, for the overall management of the MMS mission, and for mission operations. MMS is scheduled for launch in 2014 for a planned two-year mission. The MMS mission consists of four identical spacecraft flying in a tetrahedral formation in an eccentric Earth orbit. The relatively tight formation, ranging from 10 to 400 km, will provide coordinated observations giving insight into small-scale magnetic field reconnection processes. By varying the size of the tetrahedron and the orbital semi-major axis and eccentricity, and making use of the changing solar phase, this geometry allows for the study of both bow shock and magnetotail plasma physics, including acceleration, reconnection, and turbulence. The mission divides into two phases for science; these phases will have orbit dimensions of l.2xl2 Earth radii in the first phase and l.2x25 Earth radii in the second in order to study the dayside magnetopause and the nightside magnetotail, respectively. The orbital periods are roughly one day and three days for the two mission phases. Each of the four MMS spacecraft will be spin stabilized at 3 revolutions per minute (rpm), with the spin axis oriented near the ecliptic north pole but tipped approximately 2.5 deg towards the Sun line. The main body of each spacecraft will be an eight-sided platform with diameter of 3.4 m and height of 1.2 m. Several booms are attached to this central core: two axial booms of 14.9 m length, two radial magnetometer booms of 5 m length, and four radial -wire booms of 60 m length. Attitude and orbit control will use a set of axial and radial thrusters. A four-head star tracker and a slit-type digital Sun sensor (DSS) provide input for attitude determination. In addition, an accelerometer will be· used for closed-loop orbit maneuver control. The primary AGS product will be a daily definitive attitude history. Due to power limitations; the star tracker and accelerometer data will not be available at all times. However, tracker data from at least 10 percent of each orbit and continuous DSS data will be provided. An extended Kalman filter (EKF) will be used to estimate the three-axis attitude (i.e., spin axis orientation and spin phase) and rotation rate for all times when the tracker data is valid. For other times, the attitude is generated by assuming a constant angular momentum vector in the inertial frame. The DSS sun pulse will provide a timing signal to maintain an accurate spin phase. There will be times when the Sun is occulted and DSS data is not available. If this occurs at the start or end of a definitive attitude product, then the spin phase will be extrapolated using the mean rate determined by the EKF.

Sedlak, Joseph E.↗
NASA NTRS1999

Uncertainty Analysis of Instrument Calibration and Application

Experimental aerodynamic researchers require estimated precision and bias uncertainties of measured physical quantities, typically at 95 percent confidence levels. Uncertainties of final computed aerodynamic parameters are obtained by propagation of individual measurement uncertainties through the defining functional expressions. In this paper, rigorous mathematical techniques are extended to determine precision and bias uncertainties of any instrument-sensor system. Through this analysis, instrument uncertainties determined through calibration are now expressed as functions of the corresponding measurement for linear and nonlinear univariate and multivariate processes. Treatment of correlated measurement precision error is developed. During laboratory calibration, calibration standard uncertainties are assumed to be an order of magnitude less than those of the instrument being calibrated. Often calibration standards do not satisfy this assumption. This paper applies rigorous statistical methods for inclusion of calibration standard uncertainty and covariance due to the order of their application. The effects of mathematical modeling error on calibration bias uncertainty are quantified. The effects of experimental design on uncertainty are analyzed. The importance of replication is emphasized, techniques for estimation of both bias and precision uncertainties using replication are developed. Statistical tests for stationarity of calibration parameters over time are obtained.

Tripp, John S.↗
NASA NTRS2012

The Fixed-bias Langmuir Probe on the Communication-navigation Outage Forecast System Satellite: Calibration and Validation

A fixed-bias spherical Langmuir probe is included as part of the Vector Electric Field Instrument (VEFI) suite on the Communication Navigation Outage Forecast System (CNOFS) satellite.CNOFS gathers data in the equatorial ionosphere between 400 and 860 km, where the primary constituent ions are H+ and O+. The ion current collected by the probe surface per unit plasmadensity is found to be a strong function of ion composition. The calibration of the collected current to an absolute density is discussed, and the performance of the spherical probe is compared to other in situ instruments on board the CNOFS satellite. The application of the calibration is discussed with respect to future xed-bias probes; in particular, it is demonstrated that some density fluctuations will be suppressed in the collected current if the plasma composition rapidly changes along with density. This is illustrated in the observation of plasma density enhancements on CNOFS.

VEFI↗
NASA NTRS2012

The Fixed-Bias Langmuir Probe on the Communication-Navigation Outage Forecast System Satellite: Calibration and Validation

A fixed-bias spherical Langmuir probe is included as part of the Vector Electric Field Instrument (VEFI) suite on the Communication Navigation Outage Forecast System (CNOFS) satellite.CNOFS gathers data in the equatorial ionosphere between 400 and 860 km, where the primary constituent ions are H+ and O+. The ion current collected by the probe surface per unit plasma density is found to be a strong function of ion composition. The calibration of the collected current to an absolute density is discussed, and the performance of the spherical probe is compared to other in situ instruments on board the CNOFS satellite. The application of the calibration is discussed with respect to future fixed-bias probes; in particular, it is demonstrated that some density fluctuations will be suppressed in the collected current if the plasma composition rapidly changes along with density. This is illustrated in the observation of plasma density enhancements on CNOFS.

Vector Electric Field Instrument↗
DOE OSTI2023

An empirical approach to model selection: weak lensing and intrinsic alignments

ABSTRACT In cosmology, we routinely choose between models to describe our data, and can incur biases due to insufficient models or lose constraining power with overly complex models. In this paper, we propose an empirical approach to model selection that explicitly balances parameter bias against model complexity. Our method uses synthetic data to calibrate the relation between bias and the χ2 difference between models. This allows us to interpret χ2 values obtained from real data (even if catalogues are blinded) and choose a model accordingly. We apply our method to the problem of intrinsic alignments – one of the most significant weak lensing systematics, and a major contributor to the error budget in modern lensing surveys. Specifically, we consider the example of the Dark Energy Survey Year 3 (DES Y3), and compare the commonly used non-linear alignment (NLA) and tidal alignment and tidal torque (TATT) models. The models are calibrated against bias in the Ωm–S8 plane. Once noise is accounted for, we find that it is possible to set a threshold Δχ2 that guarantees an analysis using NLA is unbiased at some specified level Nσ and confidence level. By contrast, we find that theoretically defined thresholds (based on, e.g. p-values for χ2) tend to be overly optimistic, and do not reliably rule out cosmological biases up to ∼1–2σ. Considering the real DES Y3 cosmic shear results, based on the reported difference in χ2 from NLA and TATT analyses, we find a roughly $30{{\ \rm per\ cent}}$ chance that were NLA to be the fiducial model, the results would be biased (in the Ωm–S8 plane) by more than 0.3σ. More broadly, the method we propose here is simple and general, and requires a relatively low level of resources. We foresee applications to future analyses as a model selection tool in many contexts.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2011

Uncertainty Assessment of the SeaWiFS On-Orbit Calibration

Ocean color climate data records require water-leaving radiances with 5% absolute and 1% relative accuracies as input. Because of the amplification of any sensor calibration errors by the atmospheric correction, the 1% relative accuracy requirement translates into a 0.1% long-term radiometric stability requirement for top-of-the atmosphere radiances. The rigorous on-orbit calibration program developed and implemented for SeaWiFS by the NASA Ocean Biology Processing Group (OBPG) Calibration and Validation Team (CVT) has allowed the CVT to maintain the stability of the radiometric calibration of SeaWiFS at 0.13% or better over the mission. The uncertainties in the resulting calibrated top-of-the-atmosphere (TOA) radiances can be addressed in terms of accuracy (biases in the measurements), precision (scatter in the measurements), and stability (repeatability of the measurements). The calibration biases of lunar observations relative to the USGS RObotic Lunar Observatory (ROLO) photometric model of the Moon are 2-3%. The biases from the vicarious calibration against the Marine Optical Buoy (MOBY) are 1-2%. The precision of the calibration derived from the solar calibration signal-tonoise ratios are 0.16%, from the lunar residuals are 0.13%, and from the vicarious gains are 0.10%. The long-term stability of the TOA radiances, derived from the lunar time series, is 0.13%. The stability of the vicariouslycalibrated TOA radiances, incorporating the uncertainties in the MOBY measurements and the atmospheric correction, is 0.30%. These results allow the OBPG to produce climate data records from the SeaWiFS ocean color data.

Eplee, Robert E., Jr.↗
NASA NTRS2011

A Comparison of Methods for a Priori Bias Correction in Soil Moisture Data Assimilation

Data assimilation is being increasingly used to merge remotely sensed land surface variables such as soil moisture, snow and skin temperature with estimates from land models. Its success, however, depends on unbiased model predictions and unbiased observations. Here, a suite of continental-scale, synthetic soil moisture assimilation experiments is used to compare two approaches that address typical biases in soil moisture prior to data assimilation: (i) parameter estimation to calibrate the land model to the climatology of the soil moisture observations, and (ii) scaling of the observations to the model s soil moisture climatology. To enable this research, an optimization infrastructure was added to the NASA Land Information System (LIS) that includes gradient-based optimization methods and global, heuristic search algorithms. The land model calibration eliminates the bias but does not necessarily result in more realistic model parameters. Nevertheless, the experiments confirm that model calibration yields assimilation estimates of surface and root zone soil moisture that are as skillful as those obtained through scaling of the observations to the model s climatology. Analysis of innovation diagnostics underlines the importance of addressing bias in soil moisture assimilation and confirms that both approaches adequately address the issue.

Kumar, Sujay V.↗
DOE OSTI2026

Automatic Calibration and Health Monitoring of Infrastructure Sensors

Smart transportation infrastructure relies on networks of heterogeneous sensors - cameras, radars, and lidars - continuously monitoring traffic conditions. However, executing the initial spatial calibration of multiple sensors and the subsequent health monitoring presents significant operational challenges. Environmental factors, mechanical vibrations, and gradual drift cause spatial misalignment, degrading fusion performance and tracking accuracy. Traditional calibration approaches require manual intervention with specialized targets or survey equipment, resulting in service interruptions and high maintenance costs. This work presents an automated framework for initial calibration and continuous health monitoring without human intervention or service disruption. Our approach addresses two critical problems: (1) detecting when sensors become miscalibrated during operation, and (2) automatically re-establishing spatial alignment using only operational traffic data. The health monitoring component analyzes measurement innovations - differences between sensor observations and predicted object states - to detect systematic biases indicative of calibration drift. By computing bias magnitude, directional consistency, and rejection rates, the system identifies miscalibrations as small as 0.5 meters. Unlike traditional methods requiring known calibration targets, our diagnostic operates continuously on live traffic observations, enabling early detection before fusion quality degrades. The automatic recalibration algorithm leverages overlapping sensor fields-of-view and temporal correlation of vehicle observations. Using graph-based optimization, the system automatically discovers which sensor pairs observe common regions, estimates pairwise spatial transformations using RANSAC-based robust estimation, and jointly optimizes all sensor poses through bundle adjustment. The framework handles practical deployment challenges, including different sensor sampling rates (1-10 Hz), varying installation positions, unknown orientations, and limited overlap regions (>10%). When approximate sensor positions are available from installation surveys (+/-1m accuracy), the algorithm additionally estimates sensor orientations, refining both position and rotation to sub-meter and sub-degree accuracy. We validate the framework on multi-hour traffic datasets from six heterogeneous sensors with sampling rates ranging from 1 Hz to 10 Hz. Results demonstrate successful calibration even with sparse overlap (<20%) and automatic detection of miscalibrations exceeding 0.8 meters. This work enables a "deploy-and-forget" sensor infrastructure that maintains calibration autonomously, reducing maintenance costs while improving tracking accuracy. The techniques generalize beyond transportation to any multi-sensor monitoring application requiring robust spatial alignment, including smart cities, industrial monitoring, and surveillance systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Understanding Biases in Sample Preparation Techniques for Coupled Scanning Electron Microscopy and MAMA PuO 2 Morphological Analysis

In this project, the scanning electron microscopy (SEM) sampling method used during the statistical design study (SDS) was investigated to determine if any sampling biases were present in the analyzed data. Using standard particle size distribution powders from the National Institute of Standards and Technology (NIST 1984 standard reference material) with the origin wet dispersion method, it was determined that a bias to smaller particles was present. This was supported by theoretical calculations using Stokes’ law to determine the settling rate of spherical particles of roughly the same size and mass as those found in the SDS. Based on the theoretical calculations, it was determined that the settling rate for each of the 76 powder sets in the SDS could be unique based on specific particle shape and mass distributions, making a universal correction factor/formula not applicable. Therefore, priority shifted to developing an improved wet dispersion method that significantly reduced the particle settling rate for all particle size and shapes. This was achieved by replacing the original solvent (isopropyl alcohol) with a heavy liquid (lithium heteropolytungstates), which dramatically slowed the settling rate and allowed for the capture of a suitable homogeneous aliquot. SEM imaging and Morphological Analysis for Material Attribution (MAMA) software analysis were conducted on the NIST standard, and the SEM/MAMA data were compared to data captured by a dynamic image analysis particle size analyzer. The resulting data confirmed that the new wet dispersion method does indeed deliver an improved representative aliquot to the SEM stub. For instance, in the NIST certificate, the average particle size is ~17.1 µm ± 2.2 µm with a normal distribution. The initial wet dispersion method resulted in a drastically reduced average particle size of 6.1 µm in addition to a non-representative heavy bi-modal distribution whereas the improved LST wet dispersion method resulting in an average particle size that was much closer to the NIST certificate (12.7 µm) with a similar normal distribution. Although the improved method was still short of the NIST certificate average, atomic force microscopy analysis determined that the resulting ~20-25% reduction in size was due to particles sinking into the carbon sticky tape used for SEM imaging. It is believed that that this bias can be calibrated in a much more predicable manner than the original settling rate bias. In addition, the matching normal distribution curves between the NIST certificate and the heavy liquid method indicate a much-improved representative aliquot has been sampled and imaged. A surrogate CeO 2 powder was used to reflect PuO 2 more accurately and to aid in implementing radiological controls and shielding. The resulting data sets from the SEM/MAMA method and the particle size analyzer give almost identical average particle sizes and particle distribution statistics. Future work will re-analyze several select runs from the SDS to determine if morphological signatures can be found with the improved sampling method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

OARE flight maneuvers and calibration measurements on STS-58

The Orbital Acceleration Research Experiment (OARE), which has flown on STS-40, STS-50, and STS-58, contains a three axis accelerometer with a single, nonpendulous, electrostatically suspended proofmass which can resolve accelerations to the nano-g level. The experiment also contains a full calibration station to permit in situ bias and scale factor calibration. This on-orbit calibration capability eliminates the large uncertainty of ground-based calibrations encountered with accelerometers flown in the past on the orbiter, thus providing absolute acceleration measurement accuracy heretofore unachievable. This is the first time accelerometer scale factor measurements have been performed on orbit. A detailed analysis of the calibration process is given along with results of the calibration factors from the on-orbit OARE flight measurements on STS-58. In addition, the analysis of OARE flight maneuver data used to validate the scale factor measurements in the sensor's most sensitive range is also presented. Estimates on calibration uncertainties are discussed. This provides bounds on the STS-58 absolute acceleration measurements for future applications.

Blanchard, Robert C.↗
NASA NTRS1995

Summary of OARE flight calibration measurements

To date, the Orbital Acceleration Research Experiment (OARE) has flown on the shuttle orbiter for five missions; namely, STS-40, STS-50, STS-58, STS-62, and STS-65. The OARE instrument system contains a three-axis accelerometer which can resolve accelerations to the nano-g (10(exp -9) g) level and a full calibration station to permit in situ bias and scale factor calibration measurements. This calibration capability eliminates the large uncertainty encountered with accelerometers flown in the past on the orbiter which use ground-based calibrations to provide absolute acceleration measurements. A detailed flight data report presentation is given for the OARE calibration measurements from all missions, along with an estimate of the calibration errors. The main aim is to collect, process, and present the calibration data in one archival report. These calibration data are the necessary key ingredient to produce the absolute acceleration levels from the OARE acceleration flight data.

Blanchard, Robert C.↗
NASA NTRS2004

Ar-Ar and I-Xe Ages and the Thermal History of IAB Meteorites

Studies of several samples of the large Caddo County IAB iron meteorite reveal andesitic material, enriched in Si, Na, Al and Ca, which is essentially unique among meteorites. This material is believed to have formed from a chondritic source by partial melting and to have further segregated by grain coarsening. Such an origin implies extended metamorphism of the IAB parent body. New Ar-39-Ar-40 ages for silicate from three different Caddo samples are consistent with a common age of 4.50-4.51 Gyr ago. Less well defined Ar-Ar degassing ages for inclusions from two other IABs, EET8333 and Udei Station, are approx.4.32 Gyr, whereas the age for Campo del Cielo varies considerably over approx.3.23-4.56 Gyr. New I-129-Xe-129 ages for Caddo County and EET8333 are 4561.9+/-0.1 Myr and 4560- 4563 Myr, respectively, relative to an age of 4566 Myr for Shallowater. Considering all reported Ar-Ar ages for IABs and related winonaites, the range is approx.4.32-4.53 Gyr, but several IABs give similar Ar ages of 4.50-4.52 Gyr. We interpret these older ages to represent cooling after the time of last significant metamorphism on the parent body, and the younger ages to represent later Ar-40 diffusion loss. These older Ar-Ar ages are similar to Sm-Nd and Rb-Sr isochron ages reported in the literature for Caddo County. Considering the possibility that IAB parent body formation was followed by impact disruption, reassembly, and metamorphism (e.g., Benedix et al. 2000), the time of the post-assembly metamorphism may have been as late as approx.4.53 Gyr ago. However, precise I-Xe ages reported for some IABs define a range of ages of approx.4560 to approx.4576 Myr. The older I-Xe ages exceed the oldest precise radiometric ages of meteorites, appear unrealistic, and s,uggest a bias in the calibration of all I-Xe ages. But even with such a bias, the I-Xe ages of IABs cannot easily be reconciled with the much younger Ar-Ar and Sm-Nd ages and with cooling rates deduced from Ni concentration profiles in IAB metal (Herpfer et al., 1994). An explanation for the difference in radiometric ages of IABs may reside in combinations of the following: a) I-Xe ages have very high closure temperatures and were not reset during metamorphism; b) a bias exists in the K-40 decay constants; c) the reported Sm-Nd and Rb-Sr ages for Caddo are in error by amounts equal to or exceeding their reported 2-sigma uncertainties; and 4) the IAB parent body may have experienced a mild metamorphism approx.30 Myr after the initial heating that produced differentiation of Caddo silicate and mixing of silicate and metal.

Bogard, Donald D.↗
NASA NTRS2004

Ar-Ar and I-Xe Ages and the Thermal History of IAB Meteorites

Studies of several samples of the large Caddo County IAB iron meteorite reveal andesitic material, enriched in Si, Nay Al and Ca, which is essentially unique among meteorites. This material is believed to have formed from a chondritic source by partial melting and to have further segregated by grain coarsening. Such an origin implies extended metamorphism of the IAB parent body. New Ar-39-Ar-40 ages for silicate from three different Caddo samples are consistent with a common age of 4.50- 4.51 Gyr ago. Less well defined Ar-Ar degassing ages for inclusions from two other IABs, EET8333 and Udei Station, are approx.4.32 Gyr, whereas the age for Campo del Cielo varies considerably over approx.3.23-4.56 Gyr. New I-129-Xe-129 ages for Caddo County and EET8333 are 4561.9 +/-0.1 Myr and 4560-4563 Myr, respectively, relative to an age of 4566 Myr for Shallowater. Considering all reported Ar-Ar ages for IABs and related winonaites, the range is approx.4.32-4.53 Gyr, but several IABs give similar Ar ages of 4.50-4.52 Gyr. We interpret these older ages to represent cooling after the time of last significant metamorphism on the parent body, and the younger ages to represent later Ar-40 diffusion loss. These older Ar-Ar ages are similar to Sm-Nd and Rb-Sr isochron ages reported in the literature for Caddo County. Considering the possibility that IAB parent body formation was followed by impact disruption, reassembly, and metamorphism (e.g., Benedix et al. 2000), the time of the post-assembly metamorphism may have been as late as approx.4.53 Gyr ago. However, precise I-Xe ages reported for some IABs define a range of ages of approx.4560 to approx.4576 My. The older I-Xe ages exceed the oldest precise radiometric ages of meteorites, appear unrealistic, and suggest a bias in the calibration of all I-Xe ages. But even with such a bias, the I-Xe ages of IABs cannot easily be reconciled with the much younger Ar-Ar and Sm-Nd ages and with cooling rates deduced from Ni concentration profiles in IAB metal (Herpfer et al., 1994). An explanation for the difference in radiometric ages of IABs may reside in combinations of the following: a) I-Xe ages have very high closure temperatures and were not reset during metamorphism; b) a bias exists in the K-40 decay constants; c) the reported Sm-Nd and Rb-Sr ages for Caddo are in error by amounts equal to or exceeding their reported 2-sigma uncertainties; and 4) the IAB parent body may have experienced a mild metamorphism approx.30 Myr after the initial heating that produced differentiation of Caddo silicate and mixing of silicate and metal. 129 129

Bogard, Donald D.↗