Search NASA⌕ Search

SEARCH · Search NASA

Results for “systematic errors”

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 775 records · Page 43

Preventing Anomalous Torques in Circumbinary Accretion Simulations

Numerical experiments are the primary method of studying the evolution of circumbinary disks due to the strong nonlinearities involved. Many circumbinary simulations also require the use of numerical mass sinks: source terms that prevent gas from unphysically accumulating around the simulated point masses by removing gas at a given rate. However, special care must be taken when drawing physical conclusions from such simulations to ensure that results are not biased by numerical artifacts. We demonstrate how improved sink methods reduce some of these biases, using simulations of aspect ratio 0.1 accretion disks around binaries with mass ratios between 0.1 and 1. We show that sink terms that do not reduce the angular momentum of gas relative to the accreting object (1) reduce the dependence on the sink rate of quantities such as the torque on the binary, distribution of accretion between binary components, and evolution of the binary semimajor axis; (2) reduce the degree to which the sink rate affects the structure of the accretion disks around each binary component; (3) alter the variability of accretion onto the binary, regularizing it in time. We also investigate other potential sources of systematic error, such as gravitational softening and simplifications to the viscous stress tensor. Because of the strong dependence of binary orbital evolution on both the torque and the distribution of mass between binary components, the sink method influences the orbital evolution of the binary at all mass ratios, with mass ratios below ∼0.3 most affected.

Alexander Dittmann↗

Constraining the Neutron Star Mass–Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes

We describe the X-ray pulse profile models we use and how we use them to analyze Neutron Star Interior Composition Explorer(NICER)observations of rotation-powered millisecond pulsars to obtain information about the mass–radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we analyzed in Miller et al. and Riley et al. and describe a cross-verification of computations of the pulse profiles of a star with R/M 3, in case stars this compact need to be considered in future analyses. We also present our early cross-verification efforts of the parameter estimation procedures used by Miller et al. and Riley et al. by analyzing two distinct synthetic data sets. Both codes yielded credible regions in the mass–radius plane that are statistically consistent with one another, and both gave posterior distributions for model parameter values consistent with the values that were used to generate the data. We also summarize the additional tests of the parameter estimation procedure of Miller et al. that used synthetic pulse profiles and the NICER pulse profile of PSR J0030+0451. We then illustrate how the precision of mass and radius estimates depends on the pulsar’s spin rate and the size of its hot spot by analyzing four different synthetic pulse profiles. Finally, we assess possible sources of systematic error in the estimates made using this technique, some of which may warrant further investigation.

Slavko Bogdanov↗

Application of a new X-ray reflection model to V1223 Sagittarii

In intermediate polars (IPs), the intrinsic thermal emissions from white dwarfs (WDs) have typically been studied. Few reports have analysed X-ray reflections from WDs. We recently developed an elaborate IP-reflection spectral model. Herein, we report the first application of a reflection model with an IP thermal model to the spectra of the brightest typical IP V1223 Sagittarii observed by the Suzaku and Nuclear Spectroscopic Telescope Array satellites. The model reasonably reproduces the spectra within the range of 5–78 keV and estimates the WD mass as 0.92 ± 0.02 Mꙩ. The WD mass estimated by the proposed model is consistent with that measured using an active galactic nucleus reflection model and a partial covering absorption model. However, the choice of incorrect parameter values, such as an unsuitable fitting energy band and an incorrect metal abundance, was found to introduce systematic errors (e.g. ≲0.2 Mꙩ in the WD mass) in the WD mass measurement. Our spin-phase-resolved analysis resulted in discoveries regarding the modulations of the equivalent width of the fluorescent iron K α line and the angle between the post-shock accretion column and the line of sight (viewing angle). The viewing angle anticorrelates approximately with the X-ray flux and has average and semi-amplitude values of 55° and 7°, respectively, which points towards two WD spin axis angles from the line of sight of 55° and 7°, respectively. Both estimated spin axis angles are different from the reported system inclination of 24°.

Takayuki Hayashi↗

Impact of the Screen Configurations on the Estimated Terra MODIS SD Degradation

The Moderate Resolution Imaging Spectroradiometer (MODIS) instruments on-board the Aqua and Terra spacecraft have collected valuable Earth data for the last 19 and 21 years, respectively. MODIS is equipped with various on-orbit calibrators, including a solar diffuser (SD) and solar diffuser stability monitor (SDSM), that are used to monitor changes in the instrument’s gain over time. Nominally, SD calibrations alternate between two configurations: screen open and screen closed. Terra MODIS, however, experienced an anomaly in 2003, which has left the SD door and screen in a permanent open-closed configuration. This resulted in accelerated degradation of the SD on Terra MODIS due to the direct solar radiation exposure every orbit. It also led to an unexpected divergence between the calibration coefficients, m1, or the inverse of the gain, generated using SD data and those generated using lunar observations at the time of the 2003 anomaly. This paper examines the effect of the screen configuration on the Terra SD degradation by analyzing Terra SDSM data and comparing to Aqua. We generate “pseudo-open” Terra SDSM data by modeling the ratio of Aqua screen open and screen closed SDSM data and calculate the Terra SD degradation using this new screen open data. We then examine the effects on our m1 values calculated using this new pseudo-open SD degradation for Terra. Implementing this new degradation results in a smaller discrepancy between the lunar and SD m1 values after the door anomaly, indicating that there may be a systematic error in the SD degradation calculated with screen closed data.

MODIS↗

Towards Energy Scale Calibration and Drift Correction of TES Detectors for Athena X-IFU

The Athena X-Ray Integral Field Unit (X-IFU) comprises a 2376-pixel array of transition edge sensors (TES) read out with time-division multiplexing (TDM). X-IFU will provide spatially resolved, high-resolution spectroscopy (2.5 eV full-width-half-maximum up to 7 keV) over the energy range 0.2 to 12 keV, with an absolute energy scale accuracy of 0.4 eV. The energy scale function maps the optimally filtered pulse height, in arbitrary engineering units, to real calibrated energy. Uncertainties in the calibration can result from imperfect fitting of the energy scale between the known calibration points. Furthermore, temporal changes in the TES operating environment, such as heat-sink temperature, magnetic field and bias voltage, can cause significant variations in the detector gain function over time. If not properly corrected, this can result in degradation of the energy resolution, and systematic errors in the absolute energy scale. The non-linear nature of TES detectors, coupled with the possibility of multiple simultaneously occurring sources of drift, can make effective corrections over the full bandpass of the instrument extremely challenging. Athena X-IFU will employ an on-board calibration source that provides known reference x-ray lines. This provides real-time monitoring of the gain stability of the detector system and information that can be used to correct for gain drifts. For X-IFU the baseline approach is to measure a series of calibration curves under different environmental conditions, which bound the expected drifts the instrument is predicted to see over the course of the mission. Using the information from the in-flight calibration source, these energy scale functions can be interpolated to generate a new corrected energy scale as a function of time. In this paper we discuss progress towards demonstrating that the X-IFU energy scale requirements can be met. We present measurements on ~ 200 pixels in a prototype X-IFU array read out with 8-column x 32-row TDM. We use a rotating target source containing 12 fluorescent targets to generate x-ray lines covering the energy range 4 keV (Sc-Kα) to 12 keV (Br-Kα). We present measurements of the non-linear energy scale function and show how variations in heat-sink temperature, TES bias voltage and magnetic field affect the shape of TES energy scale differently and introduce different residual gain errors over the bandpass. We explore different drift correction algorithms that use either a single or multiple referential lines to track and correct the gain from these various sources of drift. In addition to the pulse-height, the DC ‘baseline’ level of the TES can contain information about its bias conditions. Thus, we test a multi-parameter gain correction algorithm that attempts to incorporate both the pulse height and the additional baseline information into the algorithm.

Stephen J Smith↗

Developing a Spectral Correlation based Method for mitigating Angular Mismatch Effects in Intercalibration Using a Benchmark Hyperspectral Sensor

The CLARREO Pathfinder (CPF) mission will implement a state-of-the-art intercalibration method for transferring CPF’s in-orbit Système Internationale (SI)-traceable reference to the shortwave channel of the Clouds and the Earth’s Radiance Energy System (CERES) and the reflective solar bands of the Visible Infrared Imaging Radiometer Suite (VIIRS) aboard the NOAA-20 satellite platform with a targeted intercalibration methodology uncertainty of 0.3% (k=1). In order to achieve such a high intercalibration accuracy, the CPF intercalibration measurements to be scheduled need to closely match those from CERES and VIIRS in time, space, angles, and wavelength. Despite’s CPF’s two-axis pointing capability to match its boresight line of sight to that of a target sensor, there will be finite residual differences in angular samplings of the two instruments. The potential angular mismatch between CPF and the target sensors can introduce systematic errors in the intercalibration results, and therefore needs to be corrected using a rigorously designed algorithm. The CPF intercalibration team has been developing a correction method for mitigating the impact of these angular anisotropic effects in the CPF-VIIRS and CPF-CERES intercalibration. The method explores the spectral correlation relationship between the reflected solar radiances from the same surface target that would be measured by CPF at two adjacent angles. Our studies have shown that the spectral information based on CPF measurements can be used to accurately predict the spectral radiance or reflectance difference due to a given mismatch in the viewing and solar geometry. The hyper-spectral information can also be used to provide scene stratification without using any auxiliary data, which is critical to reduce the angular correction uncertainty. The angular correction relationship can be well established using simulated CPF-like top-of-atmosphere spectral radiances observed at all sorts of viewing geometry and solar angles and for different scenes. Intensive radiative transfer simulations have been completed using a state-of-art radiative transfer model developed by the CPF team members. The implementation of the algorithm on high-fidelity event simulation data and the characterization for the angular adjustment uncertainty will be presented.

Wan Wu↗

SKYSURF: Constraints on Zodiacal Light and Extragalactic Background Light through Panchromatic HST All-sky Surface-brightness Measurements: II. First Limits on Diffuse Light at 1.25, 1.4, and 1.6 μm

We present the first results from the HST Archival Legacy project "SKYSURF." As described in Windhorst et al., SKYSURF utilizes the large HST archive to study the diffuse UV, optical, and near-IR backgrounds and foregrounds in detail. Here, we utilize SKYSURF's first sky-surface-brightness measurements to constrain the level of near-IR diffuse Extragalactic Background Light (EBL) in three near-IR filters (F125W, F140W, and F160W). This is done by comparing our preliminary sky measurements of >30,000 images to zodiacal light models, carefully selecting the darkest images to avoid contamination from stray light. Our sky-surface-brightness measurements have been verified to an accuracy of better than 1%, which when combined with systematic errors associated with HST, results in sky-brightness uncertainties of ∼ 2%–4% ≃ 0.005 MJy sr−1 in each image. When compared to the Kelsall et al. zodiacal model, an isotropic diffuse background of ∼30 nW m−2 sr−1 remains, whereas using the Wright zodiacal model results in no discernible diffuse background. Based primarily on uncertainties in the foreground model subtraction, we present limits on the amount of diffuse EBL of 29, 40, and 29 nW m−2 sr−1, for F125W, F140W, and F160W, respectively. While this light is generally isotropic, our modeling at this point does not distinguish between a cosmological origin or a solar system origin (such as a dim, diffuse, spherical cloud of cometary dust).

Timothy Carleton↗

Aerothermodynamic Measurements in Hypersonic Non-Equilibrium Flows

High enthalpy arc-jets are unique facilities particularly suited for simulating complex flows in the aerospace field, such as the aerothermodynamics of a re-entry vehicle. Arc-jets are often used to evaluate important design factors such as heat shield materials and vehicle design. Characterization of these facilities is important, as studies often match specific in-flight environments during experiments. Due to the complex environment produced by an arc-jet, with effects such as thermodynamic and chemical non-equilibrium occurring in the flow, characterization experiments are significantly more difficult than in traditional blow-down wind tunnels. The current work aims to characterize an arc-jet facility in terms of flow temperature, velocity and unsteadiness, in a spatially-resolved analysis. To achieve this goal, a non-intrusive imaging technique called “planar laser-induced fluorescence” was performed in the NASA Langley Hypersonic Materials Environ-mental Test System arc-jet facility. The experimental data was analyzed to produce the quantitative measurements in multiple regions of the flow around a blunt body sample. A quasi-1D low fidelity solver was created to numerically simulate the flow to investigate non-equilibrium effects occur-ring outside the imaging region in the arc-jet nozzle. Unsteadiness in the test section was minimized by analyzing a subset of data assessing the gas injection location. Radial velocity, rotational temperature and translational temperature measurements are provided that can be used to validate future high-fidelity computational studies. The temperature measurements revealed rotational non-equilibrium potentially occurring in the flow. This work provides a methodology for characterizing arc-jet facilities, identifying and correcting sources of systematic error, and estimating the uncertainty associated with the quantitative measurements.

laser induced fluorescence↗

Using Gaia Data 2 to Constrain Local Dark Matter Density and Thin Dark Disk

We use stellar kinematics from the latest Gaia data release (DR2) to measure the local dark matter (DM) density ρDM in a heliocentric cylinder of radius R = 150 pc and half-height z = 200 pc. We also explore the prospect of using our analysis to estimate the DM density in local substructure by setting constraints on the surface density and scale height of a thin dark disk aligned with the baryonic disk and formed due to dissipative dark matter self-interactions. Performing the statistical analysis within a Bayesian framework for three types of tracers, we obtain ρDM = 0.016 ± 0.010 Mꙩ/p c 3 for A stars; early G stars give a similar result, while F stars yield a significantly higher value. For a thin dark disk, A stars set the strongest constraint: excluding surface densities (5–12) Mꙩ/pc 2 for scale heights below 100 pc with 95% confidence. The upper bound of this constraint implies . 1% of the Milky Way DM mass is present in a dissipative dark sector. Comparing our results with those derived using Tycho-Gaia Astrometric Solution (TGAS) data, we find that the uncertainty in our measurements of the local DM content is dominated by systematic errors that arise from assumptions of our dynamical analysis in the low z region. Furthermore, there will only be a marginal reduction in these uncertainties with more data in the Gaia era. We comment on the robustness of our method and discuss potential improvements for future work.

dark matter theory↗

Airborne HSRL Assessment of CALIOP Aerosol Extinction Retrievals Constrained by Column AOD

Current operational CALIOP aerosol extinction profile retrievals usually rely on accurately specifying the relationship between aerosol extinction and backscattering (i.e. lidar ratio). Uncertainties in the assigned lidar ratios are typically the largest source of systematic error (~30-50%) in the CALIOP retrievals of aerosol extinction, backscatter, and aerosol optical depth. Alternatively, column aerosol optical depth (AOD) can be used to solve the lidar equation and obtain column-equivalent lidar ratios and retrieve aerosol extinction profiles from the CALIOP attenuated backscatter profiles. These derived lidar ratios correspond to the aerosols in the altitude range from 0-7 km. We derive column-equivalent aerosol lidar ratios and aerosol extinction profiles from CALIOP attenuated backscatter profiles constrained using co-located column AODs provided by several instruments/techniques: Synergized Optical Depth of Aerosols (SODA), Ocean-Derived Column Optical Depths (ODCOD), MODIS (dark target), MODIS (Multi-Angle Implementation of Atmospheric Correction-MAIAC), and PARASOL (Generalized Retrieval of Aerosol and Surface Properties-GRASP). These various retrievals of AOD and the column-average aerosol lidar ratios and aerosol extinction profiles derived from CALIOP using these AOD constraints are evaluated using the extensive record of coincident and co-located airborne HSRL measurements of AOD, aerosol lidar ratio, and aerosol extinction profiles acquired during more than 140 HSRL underflights of CALIPSO since 2006. The HSRL technique allows for the independent measurement of extinction and backscatter without the need for external constraints or assumptions regarding the lidar ratios. Initial comparisons with these coincident HSRL aerosol extinction profiles show that the CALIOP aerosol extinction profiles retrieved using these various AOD constraints are generally in better agreement with the HSRL measurements than aerosol extinction profiles from the operational techniques. The quality of agreement depends on the accuracy and magnitude of the AOD constraint. We present these comparisons of AOD, column aerosol lidar ratio, and aerosol extinction profiles for each of the AOD constraints described above.

lidar↗

Hybrid Flush and Synthetic Air Data Filter for Entry Vehicle Atmospheric State Estimation

A hybrid flush/synthetic air data sensing filter for entry vehicle atmosphere estimation is developed. The approach makes use of a Kalman-Schmidt and Rauch-Tung-Streibel smoother. The filter/smoother blends information from pressure sensors distributed on the heatshield with pseudo-measurements of the vehicle aerodynamic forces and moments computed from mass properties and inertial measurement unit data, and prior estimates of the atmosphere. The filter produces estimates of the atmospheric conditions along the entry trajectory, and systematic error estimates to reconcile differences between the pressure and aerodynamic data sources. The filter is applied to data acquired during the Mars 2020 entry, descent, and landing at Jezero crater on February 18th, 2021. The final paper will also include results from the Mars Science Laboratory entry, descent, and landing.

Chris D. Karlgaard↗

Does Collection Time Bias the Ecology of Cleanroom Air Samples?

Microbial monitoring of astromaterials collections has taken on increased importance with the return of biologically sensitive samples from the asteroids Ryugu and Bennu and the initiation of the Mars Sample Return Program. Terrestrial bacteria and fungi can alter the mineralogy and organic composition of our collections causing irreversible contamination of pristine samples and increasing the risk of false positives for life detection measurements. NASA has conducted routine microbial monitoring of its existing collections since 20181. Initial monitoring focused on surface samples collected with foam swabs. Although, airborne microbiology is often decoupled from surface microbiology in the built environment2 culture-based air sampling techniques like impactors were not compliant with existing contamination control requirements. Bringing organic rich media, gelatin or liquids into curation cleanrooms presents an unacceptable risk to pristine samples. In 2022 NASA purchased a materials complaint air sampler and began collecting air samples from the cleanrooms in addition to surface samples3. The new instrument uses an electret filter to collect samples that are suitable for cultivating organisms or for direct DNA sequencing. Preliminary DNA sequencing results appeared to indicate that longer sampling times biased the microbial community in favor of hearty, spore-forming bacteria3. We present the results of a study comparing overnight sampling (17 hours) to short (1 hour) sampling of unoccupied curation cleanrooms. The results will help us optimize our monitoring protocols and develop a more detailed inventory of the ecology of astromaterials curation cleanrooms. Methods: We analyzed 72 paired air samples from six different cleanrooms including the meteorite processing lab (ISO 7 equivalent, 16 samples), the lunar lab (ISO 6 equivalent, 10 samples), the stardust lab (ISO 5 equivalent 14 samples), the OSIRIS-REx lab (ISO 5 equivalent, 12 samples), the Hayabusa2 lab (ISO 5 equivalent, 14 samples), and the Genesis lab (ISO 4 equivalent, 6 samples). All the samples were collected with an InnovaPrep Bobcat air sampler operating at a sampling rate of 200 L/min. The sampler operates for 5 minutes out of every 20 minute period. Half of the samples were collected by filtering 3,000L (15 min. of active sampling) of air across an electret filter for one hour. The rest of the samples were collected by filtering approximately 51,000 L air across the filter overnight (~17 hours, 255 min. of active sampling). Cells were eluted from the filter using 6-7 ml of pressurized 0.15% tween 20 in PBS (phosphate buffered saline). This liquid was used to cultivate bacteria according to previously published methods1,4,5 and for DNA extraction and next generation sequencing. DNA was extracted with a Qiagen MagAttract PowerMicrobiome kit6. To identify bacteria and archaea, the 16S rRNA gene was amplified using Earth Microbiome primers for the V4 region 7. The amplified DNA was sequenced on an Illumina MiSeq using a V3 reagent kit. The resulting sequences were processed using DADA2 and QIIME2 as implemented on the EDGE bioinformatics platform8–10. Results: Only two of the 72 samples had no amplifiable DNA. Amplified DNA concentrations ranged from 2.67 – 0.272 ng/µl. The median concentration of amplified DNA for the 1 hour samples was 0.770 ± 0.368 ng/µl. The median concentration of amplified DNA for the overnight samples was 0.877 ± 0.434 ng/µl. On average the overnight samples had slightly more sequences (58,960 vs. 59,456) and ASV’s (amplicon sequence variants) (60 vs 64.5) than the one hour samples, but these differences are not statistically significant. The most abundant ASV in every sample mapped to the genus Cupravidus. ASV’s mapping to the genuses Bacillus, Schlegelella, Thermus, and Staphylococcus were also common. Discussion and Future Work: Alpha diversity statistics like Shannon Entropy and Faith Phylogenetic Diversity are used to describe the diversity of organisms in a single sample. If a longer sampling time was biasing the data, we would expect to see a change in these diversity statistics vs. sample time. However, we did not observe this in our data. The median Shannon entropy was slightly higher for the overnight samples (3.773 vs 3.611) as was the Faith Phylogenetic Diversity (4.042 vs 3.596), but both values were within a standard deviation of each other for the two sampling times (Fig. 1). It is unlikely, that the longer sampling time is introducing bias into our data. We do observe a significant decrease in diversity when comparing the air samples by lab. The Genesis lab (ISO 4 equivalent) has a lower median number of ASV’s (45.5) than the other labs (62). Median values for Shannon Entropy (3.717 vs. 3.430) and Faith Phylogenetic Diversity (3.796 vs. 3.548) are also lower for Genesis, but those values are with one standard deviation of each other for the different sampling times. This is consistent with previous culture-based results suggesting that the environment in cleanrooms tends to select for a core group of organisms capable of surviving under dry, low nutrient, conditions. The presence of the ASV’s mapping to Cupravidus and Thermus in our sequencing blanks and controls suggests that several of the most common organisms in our samples represent contaminants from the reagents used to perform the DNA extractions and sequencing. Further work is needed to identify these contaminants, remove them from our data and recalculate the diversity statistics. This is a systematic error. Therefore, we do not expect removing the sequencing contaminants to change our conclusions. Longer air sample collection times appear to result in slightly higher diversity and do not bias the results towards “hardy” bacteria like spore-formers. Based on these preliminary results we conclude that sampling at least 3,000 liters of air is sufficient to capture the microbial diversity of cleanrooms, and that air samples can also be collected overnight without negatively impacting diversity. These results allow us to be flexible when designing microbial monitoring plans so that they do not interfere with routine lab activity. References: 1. Regberg, A. B. et al. 49th Lunar and Planetary Science Conference (2018). 2. The United States Pharmacopeial Convention. USP General Chapter <1116> (2013). 3. Regberg, A. B., et al. 54th Lunar and Planetary Science Conference (2023). 4. Regberg, A. B. et al. 53rd Lunar and Planetary Science Conference ( 2022). 5. Davis, R. E.,et al. 50th Lunar and Planetary Science Conference (2019). 6. Qiagen. MagAttract® PowerMicrobiome® DNA/RNA EP Kit Handbook. (2018). 7. Walters, W. et al. mSystems 1, (2015). 8. Callahan, B. J. et al. Nat. Methods 13, 581–583 (2016). 9. Hall, M. & Beiko, R. G. Microbiome Analysis: Methods and Protocols113–129 (Springer, 2018). 10. Philipson, C. et al. Bio-Protoc. 7, e2622 (2017).

A. B. Regberg↗

Geodetic Reference Instrument Transponder for Small Satellites (Gritss)

The accuracy of the Terrestrial Reference Frame is limited by systematic errors in tying the contributions from the different geodetic techniques. Local survey uncertainties are limited due to the inaccessibility of the instrument reference points. NASA is implementing a technology demonstration mission called Geodetic Reference Instrument Transponder for Small Satellites(GRITSS) that uses a small satellite as a space-based ref-refence point for tying together collocated geodetic stations. The technology being demonstrated uses a novel idea of upconverting the Global Positioning System (GPS) signals received at the satellite and transponding them to a VLBI Global Observing System (VGOS) antenna ground station. This approach does not require the satellite to be in view of more than one VLBI station at a time, allowing the use of Low Earth Orbits and an inexpensive CubeSat. The demonstration will initially involve only the NASA VGOS stations in Maryland, Hawaii, and Texas, but may be expanded later in the mission to include other international VGOS stations. This presentation will provide an overview of the GRITSS mission and the measurement concept.

S M Merkowitz↗

Geodetic Reference Instrument Transponder for Small Satellites (GRITTS)

The accuracy of the Terrestrial Reference Frame is limited by systematic errors in tying the contributions from the different geodetic techniques. Local survey uncertainties are limited due to the inaccessibility of the instrument reference points. NASA is implementing a technology demonstration mission called Geodetic Reference Instrument Transponder for Small Satellites(GRITSS) that uses a small satellite as a space-based ref-refence point for tying together collocated geodetic stations. The technology being demonstrated uses a novel idea of upconverting the Global Positioning System (GPS) signals received at the satellite and transponding them to a VLBI Global Observing System (VGOS) antenna ground station. This approach does not require the satellite to be in view of more than one VLBI station at a time, allowing the use of Low Earth Orbits and an inexpensive CubeSat. The demonstration will initially involve only the NASA VGOS stations in Maryland, Hawaii, and Texas, but may be expanded later in the mission to include other international VGOS stations. This presentation will provide an overview of the GRITSS mission and the measurement concept.

S. M. Merkowitz↗

Geodetic Reference Instrument Transponder for Small Satellites (GRITTS)

The accuracy of the Terrestrial Reference Frame is limited by systematic errors in tying the contributions from the different geodetic techniques. Local survey uncertainties are limited due to the inaccessibility of the instrument reference points. NASA is implementing a technology demonstration mission called Geodetic Reference Instrument Transponder for Small Satellites(GRITSS) that uses a small satellite as a space-based ref-refence point for tying together collocated geodetic stations. The technology being demonstrated uses a novel idea of upconverting the Global Positioning System (GPS) signals received at the satellite and transponding them to a VLBI Global Observing System (VGOS) antenna ground station. This approach does not require the satellite to be in view of more than one VLBI station at a time, allowing the use of Low Earth Orbits and an inexpensive CubeSat. The demonstration will initially involve only the NASA VGOS stations in Maryland, Hawaii, and Texas, but may be expanded later in the mission to include other international VGOS stations. This presentation will provide an overview of the GRITSS mission and the measurement concept.

S M Merkowitz↗

Error analysis of bent ray radio occultation measurments

There are two types of experimental measurement errors of the Doppler data associated with the radio occultation, random and systematic. Random errors are due to thermal noise in the transmission channel, and the phase lock loop, and quantization error in the digital circuitry. These are called noise type errors. The systematic errors are due to geometric uncertainty and equipment phase instability. Considered is the amount of uncertainty, due to random measurement errors, in the refractivity profiles reconstructed by this type of indirect sensing experiment. A class of refractivity profiles is defined which approximately fit the set of measured data. Bounds are placed on the extent of this class of solution profiles. To accomplish this, the sensivity of the reconstructed refractivity profiles to errors in the measured quantity and the statistics of the errors in the measurement are examined.

Saintgermain, R. L.↗

Effects of vertical distribution of water vapor and temperature on total column water vapor retrieval error

Results are presented of a test of the physically based total column water vapor retrieval algorithm of Wentz (1992) for sensitivity to realistic vertical distributions of temperature and water vapor. The ECMWF monthly averaged temperature and humidity fields are used to simulate the spatial pattern of systematic retrieval error of total column water vapor due to this sensitivity. The estimated systematic error is within 0.1 g/sq cm over about 70 percent of the global ocean area; systematic errors greater than 0.3 g/sq cm are expected to exist only over a few well-defined regions, about 3 percent of the global oceans, assuming that the global mean value is unbiased.

Sun, Jielun↗

Elevation Change of the Southern Greenland Ice Sheet from Satellite Radar Altimeter Data

Long-term changes in the thickness of the polar ice sheets are important indicators of climate change. Understanding the contributions to the global water mass balance from the accumulation or ablation of grounded ice in Greenland and Antarctica is considered crucial for determining the source of the about 2 mm/yr sea-level rise in the last century. Though the Antarctic ice sheet is much larger than its northern counterpart, the Greenland ice sheet is more likely to undergo dramatic changes in response to a warming trend. This can be attributed to the warmer Greenland climate, as well as a potential for amplification of a global warming trend in the polar regions of the Northern Hemisphere. In collaboration with Drs. Curt Davis and Craig Kluever of the University of Missouri, we are using data from satellite radar altimeters to measure changes in the elevation of the Southern Greenland ice sheet from 1978 to the present. Difficulties with systematic altimeter measurement errors, particularly in intersatellite comparisons, beset earlier studies of the Greenland ice sheet thickness. We use altimeter data collected contemporaneously over the global ocean to establish a reference for correcting ice-sheet data. In addition, the waveform data from the ice-sheet radar returns are reprocessed to better determine the range from the satellite to the ice surface. At JPL, we are focusing our efforts principally on the reduction of orbit errors and range biases in the measurement systems on the various altimeter missions. Our approach emphasizes global characterization and reduction of the long-period orbit errors and range biases using altimeter data from NASA's Ocean Pathfinder program. Along-track sea-height residuals are sequentially filtered and backwards smoothed, and the radial orbit errors are modeled as sinusoids with a wavelength equal to one revolution of the satellite. The amplitudes of the sinusoids are treated as exponentially-correlated noise processes with a time-constant of six days. Measurement errors (e.g., altimeter range bias) are simultaneously recovered as constant parameters. The corrections derived from the global ocean analysis are then applied over the Greenland ice sheet. The orbit error and measurement bias corrections for different missions are developed in a single framework to enable robust linkage of ice-sheet measurements from 1978 to the present. In 1998, we completed our re-evaluation of the 1978 Seasat and 1985-1989 Geosat Exact Repeat Mission data. The estimates of ice thickness over Southern Greenland (south of 72N and above 2000 m) from 1978 to 1988 show large regional variations (+/-18 cm/yr), but yield an overall rate of +1.5 +/- 0.5 cm/yr (one standard error). Accounting for systematic errors, the estimate may not be significantly different from the null growth rate. The average elevation change from 1978 to 1988 is too small to assess whether the Greenland ice sheet is undergoing a long-term change.

Haines, Bruce J.↗