Toward decimeter-level real-time orbit determination: a demonstration using the SAC-C and CHAMP Spacecraft
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
The angular momentum of the atmosphere and oceans change as both the distribution of mass within the atmosphere and oceans changes and as the direction and speed of the winds and currents change.
Explore the source record for details and available documents.
The acronym CHAMP (signifying camera, hand lens, and microscope ) denotes any of several proposed optoelectronic instruments that would be capable of color imaging at working distances that could be varied continuously through a range from infinity down to several millimeters. As in any optical instrument, the magnification, depth of field, and spatial resolution would vary with the working distance. For example, in one CHAMP version, at a working distance of 2.5 m, the instrument would function as an electronic camera with a magnification of 1/100, whereas at a working distance of 7 mm, the instrument would function as a microscope/electronic camera with a magnification of 4.4. Moreover, as described below, when operating at or near the shortest-working-distance/highest-magnification combination, a CHAMP could be made to perform one or more spectral imaging functions. CHAMPs were originally intended to be used in robotic geological exploration of the Moon and Mars. The CHAMP concept also has potential for diverse terrestrial applications that could include remotely controlled or robotic geological exploration, prospecting, field microbiology, environmental surveying, and assembly- line inspection. A CHAMP (see figure) would include two lens cells: (1) a distal cell corresponding to the objective lens assembly of a conventional telescope or microscope and (2) a proximal cell that would contain the focusing camera lens assembly and the camera electronic image-detector chip, which would be of the active-pixel-sensor (APS) type. The distal lens cell would face outward from a housing, while the proximal lens cell would lie in a clean environment inside the housing. The proximal lens cell would contain a beam splitter that would enable simultaneous use of the imaging optics (that is, proximal and distal lens assemblies) for imaging and illumination of the field of view. The APS chip would be mounted on a focal plane on a side face of the beam splitter, while light for illuminating the field of view would enter the imaging optics via the end face of the beam splitter. The proximal lens cell would be mounted on a sled that could be translated along the optical axis for focus adjustment. The position of the CHAMP would initially be chosen at the desired working distance of the distal lens from (corresponding to an approximate desired magnification of) an object to be examined. During subsequent operation, the working distance would ordinarily remain fixed at the chosen value and the position of the proximal lens cell within the instrument would be adjusted for focus as needed.
With the successful launch of the Orsted, SAC-C and CHAMP satellites we are able to make both magnetic and gravity anomaly maps of the Earth's crust; magnetic from all three missions and gravity with CHAMP. We have used these data to study the KMA area of Russia. This is an important region for several reasons: (1) we have already made satellite magnetic anomaly maps of this region and they can be integrated with the gravity data from CHAMP for a comprehensive interpretation; (2) KMA contains the largest know reserves of iron-ore in the world; and (3) there are significant ground truth data available for this region from aeromagnetic, balloon surveys and geophysical mapping, including extensive rock magnetic/paleo-magnetic and geologic studies. Utilizing the gravity observations, collocated with the magnetic data enabled us to make a joint interpretation. While there is a high amplitude magnetic anomaly recorded over the KMA the gravity anomaly at satellite altitude revealed by CHAMP is only around 3-6 mGal but is not centered on the magnetic high. This would indicate that despite the fact that in the region of the KMA the rocks have a higher percentage of iron than in the surrounding formations the entire area is Archean-Proterozoic in age and therefore very dense.
The main problem facing current CAD tools for VLSIs is the large amount of memory required when dealing with large systems, primarily due to the circuit representation used by most current tools. This paper discusses an approach for hierarchical switch-level simulation of digital circuits. The approach exploits the hierarchy to reduce the memory requirements of the simulation, allowing the simulation of circuits that are too large to simulate at one flat level. The approach has been implemented in a hierarchical switch-level simulator, CHAMP, which runs on a SUN workstation. The program performs mixed mode simulation: parts of the circuit can be simulated faster at a behavioral level by supplying a high level software description. CHAMP allows assignable delays, and bidirectional signal flow inside circuit blocks that are represented as transistor networks as well as across the boundaries of higher level blocks. CHAMP is also unique in that it simulates directly from the hierarchical circuit description without flattening to a single level.
Data from the IMP 8 satellite are used here to rule out charged massive particles (CHAMPs) with masses between 2.4 and 56,000 TeV as the source of the dark matter in the Galactic halo. This limit is achieved under the assumption that CHAMPs are virialized.
Regional to continental scale magnetic anomaly maps are becoming increasingly available from airborne, shipborne, and terrestrial surveys. Satellite data are commonly considered to fill the coverage gaps in regional compilations of these near-surface surveys. For the near-surface Antarctic magnetic anomaly map being produced by the Antarctic Digital Magnetic Anomaly Project (ADMAP), we show that near-surface magnetic anomaly estimation is greatly enhanced by the joint inversion of the near-surface data with the satellite observations relative to the conventional technique such as minimum curvature. Orsted observations are especially advantageous relative to the Magsat data that have order-of-magnitude greater measurement errors, albeit at much lower orbital altitudes. CHAMP is observing the geomagnetic field with the same measurement accuracy as the Orsted mission, but at the lower orbital altitudes covered by Magsat. Hence, additional significant improvement in predicting near-surface magnetic anomalies can result as these CHAMP data are available. Our analysis also suggests that considerable new insights on the magnetic properties of the lithosphere may be revealed by a further order-of-magnitude improvement in the accuracy of the magnetometer measurements at minimum orbital altitude.
With the launch of CHAMP and GRACE, we have entered a new phase in the history of satellite geodesy. For the first time, geopotential models are now available based almost exclusively on satellite-satellite tracking either with GPS in the case of the CHAMP-based geopotential models, or co-orbital intersatellite ultra-precise ranging in the case of GRACE. Different groups have analyzed these data, and produced a series of geopotential models (e.g., EIGENlS, EIGEN2, GGM0lS, GGMOlC) that incorporate the new data. We will compare the performance of these "newer" geopotential models with the standard models now used for computations, (e.g., JGM-3, BGM-96, PGS7727, and GRIMS-C1) for TOPEX, JASON, Geosat-Follow-On (GFO), and Envisat using standard metrics such as SLR RMS of fit, altimeter crossovers, and orbit overlaps. Where covariances are available we can evaluate the predicted geographically correlated orbit error. These predicted results can be compared with the Earth-fixed differences between dynamic and reduced-dynamic orbits to test the predictive accuracy of the covariances, as well as to calibrate the error of the solutions.
Spatial and temporal variability of the magnetic field component induced by ocean circulation is investigated on the basis of a standard thin-shell approximation of electro- and magneto-static equations. Well-known difficulties of numerical solution of the governing equations are resolved by reducing the problem to an equation for the electric field potential,(ef) as opposed to a more conventional approach focused on the vertical jump, (psi), of the magnetic field potential across a combined ocean/ marine-sediment-layer spherical shell. The present formulation permits using more realistic input data on ocean currents and ultimately yields much greater (by at least an order of magnitude) values of the magnetic field at sea surface than predicted in earlier studies. Such large values are comparable to, and in some cases exceed, magnetic field variations caused by lithospheric and ionospheric sources on monthly to interannual timescales. At the 400-km altitude (of CHAMP satellite), the field attains 6 nT. The model predictions show favorable comparisons with some in situ measurements as well as with Challenging Minisatellite Payload (CHAMP) satellite magnetometer data.
Currently two gravity field satellite missions, CHAMP and GRACE, are equipped with high sensitivity electrostatic accelerometers, measuring the non-conservative forces acting on the spacecraft in three orthogonal directions. During the gravity field recovery these measurements help to separate gravitational and non-gravitational contributions in the observed orbit perturbations. For precise orbit determination purposes all these missions have a dual-frequency GPS receiver on board. The reduced dynamic technique combines the dense and accurate GPS observations with physical models of the forces acting on the spacecraft, complemented by empirical accelerations, which are stochastic parameters adjusted in the orbit determination process. When the spacecraft carries an accelerometer, these measured accelerations can be used to replace the models of the non-conservative forces, such as air drag and solar radiation pressure. This approach is implemented in a batch least-squares estimator of the GPS High Precision Orbit Determination Software Tools (GHOST), developed at DLR/GSOC and DEOS. It is extensively tested with data of the CHAMP and GRACE satellites. As accelerometer observations typically can be affected by an unknown scale factor and bias in each measurement direction, they require calibration during processing. Therefore the estimated state vector is augmented with six parameters: a scale and bias factor for the three axes. In order to converge efficiently to a good solution, reasonable a priori values for the bias factor are necessary. These are calculated by combining the mean value of the accelerometer observations with the mean value of the non-conservative force models and empirical accelerations, estimated when using these models. When replacing the non-conservative force models with accelerometer observations and still estimating empirical accelerations, a good orbit precision is achieved. 100 days of GRACE B data processing results in a mean orbit fit of a few centimeters with respect to high-quality JPL reference orbits. This shows a slightly better consistency compared to the case when using force models. A purely dynamic orbit, without estimating empirical accelerations thus only adjusting six state parameters and the bias and scale factors, gives an orbit fit for the GRACE B test case below the decimeter level. The in orbit calibrated accelerometer observations can be used to validate the modelled accelerations and estimated empirical accelerations computed with the GHOST tools. In along track direction they show the best resemblance, with a mean correlation coefficient of 93% for the same period. In radial and normal direction the correlation is smaller. During days of high solar activity the benefit of using accelerometer observations is clearly visible. The observations during these days show fluctuations which the modelled and empirical accelerations can not follow.
The interplanetary shock/electric field event of 5-6 November 2001 is analyzed using ACE interplanetary data. The consequential ionospheric effects are studied using GPS receiver data from the CHAMP and SAC-C satellites and altimeter data from the TOPEX/ Poseidon satellite. Data from ~100 ground-based GPS receivers as well as Brazilian Digisonde and Pacific sector magnetometer data are also used. The dawn-to-dusk interplanetary electric field was initially ~33 mV/m just after the forward shock (IMF BZ = -48 nT) and later reached a peak value of ~54 mV/m 1 hour and 40 min later (BZ = -78 nT). The electric field was ~45 mV/m (BZ = -65 nT) 2 hours after the shock. This electric field generated a magnetic storm of intensity DST = -275 nT. The dayside satellite GPS receiver data plus ground-based GPS data indicate that the entire equatorial and midlatitude (up to +/-50(deg) magnetic latitude (MLAT)) dayside ionosphere was uplifted, significantly increasing the electron content (and densities) at altitudes greater than 430 km (CHAMP orbital altitude). This uplift peaked ~2 1/2 hours after the shock passage. The effect of the uplift on the ionospheric total electron content (TEC) lasted for 4 to 5 hours. Our hypothesis is that the interplanetary electric field ''promptly penetrated'' to the ionosphere, and the dayside plasma was convected (by E x B) to higher altitudes. Plasma upward transport/convergence led to a ~55-60% increase in equatorial ionospheric TEC to values above ~430 km (at 1930 LT). This transport/convergence plus photoionization of atmospheric neutrals at lower altitudes caused a 21% TEC increase in equatorial ionospheric TEC at ~1400 LT (from ground-based measurements). During the intense electric field interval, there was a sharp plasma ''shoulder'' detected at midlatitudes by the GPS receiver and altimeter satellites. This shoulder moves equatorward from -54(deg) to -37(deg) MLAT during the development of the main phase of the magnetic storm. We presume this to be an ionospheric signature of the plasmapause and its motion. The total TEC increase of this shoulder is ~80%. Part of this increase may be due to a "superfountain effect." The dayside ionospheric TEC above ~430 km decreased to values ~45% lower than quiet day values 7 to 9 hours after the beginning of the electric field event. The total equatorial ionospheric TEC decrease was ~16%. This decrease occurred both at midlatitudes and at the equator. We presume that thermospheric winds and neutral composition changes produced by the storm-time Joule heating, disturbance dynamo electric fields, and electric fields at auroral and subauroral latitudes are responsible for these decreases.
A computational approach for simulating the effects of rolling element and journal bearings on the vibration and sound transmission through gearboxes has been demonstrated. The approach, using ARL/Penn State s CHAMP methodology, uses Component Mode Synthesis of housing and shafting modes computed using Finite Element (FE) models to allow for rapid adjustment of bearing impedances in gearbox models. The approach has been demonstrated on NASA GRC s test gearbox with three different bearing configurations: in the first condition, traditional rolling element (ball and roller) bearings were installed, and in the second and third conditions, the traditional bearings were replaced with journal and wave bearings (wave bearings are journal bearings with a multi-lobed wave pattern on the bearing surface). A methodology for computing the stiffnesses and damping in journal and wave bearings has been presented, and demonstrated for the journal and wave bearings used in the NASA GRC test gearbox. The FE model of the gearbox, along with the rolling element bearing coupling impedances, was analyzed to compute dynamic transfer functions between forces applied to the meshing gears and accelerations on the gearbox housing, including several locations near the bearings. A Boundary Element (BE) acoustic model was used to compute the sound radiated by the gearbox. Measurements of the Gear Mesh Frequency (GMF) tones were made by NASA GRC at several operational speeds for the rolling element and journal bearing gearbox configurations. Both the measurements and the CHAMP numerical model indicate that the journal bearings reduce vibration and noise for the second harmonic of the gear meshing tones, but show no clear benefit to using journal bearings to reduce the amplitudes of the fundamental gear meshing tones. Also, the numerical model shows that the gearbox vibrations and radiated sound are similar for journal and wave bearing configurations.
Objective quantification of model performance based on metrics helps us evaluate the current state of space physics modeling capability, address differences among various modeling approaches, and track model improvements over time. The Coupling, Energetics, and Dynamics of Atmospheric Regions (CEDAR) Electrodynamics Thermosphere Ionosphere (ETI) Challenge was initiated in 2009 to assess accuracy of various ionosphere/thermosphere models in reproducing ionosphere and thermosphere parameters. A total of nine events and five physical parameters were selected to compare between model outputs and observations. The nine events included two strong and one moderate geomagnetic storm events from GEM Challenge events and three moderate storms and three quiet periods from the first half of the International Polar Year (IPY) campaign, which lasted for 2 years, from March 2007 to March 2009. The five physical parameters selected were NmF2 and hmF2 from ISRs and LEO satellites such as CHAMP and COSMIC, vertical drifts at Jicamarca, and electron and neutral densities along the track of the CHAMP satellite. For this study, four different metrics and up to 10 models were used. In this paper, we focus on preliminary results of the study using ground-based measurements, which include NmF2 and hmF2 from Incoherent Scatter Radars (ISRs), and vertical drifts at Jicamarca. The results show that the model performance strongly depends on the type of metrics used, and thus no model is ranked top for all used metrics. The analysis further indicates that performance of the model also varies with latitude and geomagnetic activity level.
Cycle 23 was exceptional in that it lasted almost two years longer than its predecessors and in that it ended in an extended minimum period that proved all predictions wrong. Comparisons of the International Reference Ionosphere (IRI) with CHAMP and GRACE in-situ measurements of electron density during the minimum have revealed significant discrepancies at 400-500 km altitude. Our study investigates the causes for these discrepancies with the help of ionosonde and Planar Langmuir Probe (PLP) data from the Communications/Navigation Outage Forecasting System (C/NOFS) satellite. Our C/NOFS comparisons confirm the earlier CHAMP and GRACE results. But the ionosonde measurements of the F-peak plasma frequency (foF2) show generally good agreement throughout the whole solar cycle. At mid-latitude stations yearly averages of the data-model difference are within 10% and at low latitudes stations within 20%. The 60-70% differences found at 400-500 km altitude are not seen at the F peak. We will discuss how these seemingly contradicting results from the ionosonde and in situ data-model comparisons can be explained and which parameters need to be corrected in the IRI model.
We investigate satellite orbital drag effects at lowEarth orbit associated with thermosphere heating during magnetic storms caused by coronal mass ejections. CHAllenge Minisatellite Payload (CHAMP) and Gravity Recovery And Climate Experiment (GRACE) neutral density data are used to compute orbital drag. Stormtoquiet density comparisons are performed with background densities obtained by the JacchiaBowman 2008 (JB2008) empirical model. Our storms are grouped in different categories regarding their intensities as indicated by minimum values of the SYMH index. We then perform superposed epoch analyses with storm main phase onset as zero epoch time. In general, we find that orbital drag effects are larger for CHAMP (lower altitudes) in comparison to GRACE (higher altitudes). Results show that storm time drag effects manifest first at high latitudes, but for extreme storms, particularly observed by GRACE, stronger orbital drag effects occur during early main phase at low/equatorial latitudes, probably due to heating propagation from high latitudes. We find that storm time orbital decay along the satellites' path generally increases with storm intensity, being stronger and faster for the most extreme events. For these events, orbital drag effects decrease faster probably due to elevated cooling effects caused by nitric oxide, which introduce modeled density uncertainties during storm recovery phase. Errors associated with total orbit decay introduced by JB2008 are generally the largest for the strongest storms and increase during storm times, particular during recovery phases. We discuss the implication of these uncertainties for the prediction of collision between space objects at lowEarth orbit during magnetic storms.
We examine the accuracy of density prediction by the first principles model Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) developed by the National Center for Atmospheric Research and compare it to the accuracy of three empirical models: Jacchia 71, the Naval Research Laboratory Mass Spectrometer Incoherent Scatter Extended 2000 (NRLMSIS), Jacchia 1971, and Jacchia-Bowman 2008. Comparisons are made for three large storms: the October 2003 storm, the March 2013 storm, and the March 2015 storm. To evaluate the accuracy of these models we use tracking data for nine space objects in low Earth orbit. Additionally, we evaluate the accuracy of the TIEGCM and NRLMSIS with data from high precision accelerometers on the Challenging Minisatellite Payload (CHAMP) and Gravity field and Circulation Explorer (GOCE) satellites. The goal is to assess the use of a first principles model as a potential tool for forecasting satellite drag during large magnetic storms. For the storms considered, we found the TIEGCM, JB2008, and NRLMSIS models to be substantially more accurate than the Jacchia 71 model. The accuracies of the TIEGCM and JB2008 models were similar, but overall, the TIEGCM was more accurate. We found smaller differences for TIEGCM versus CHAMP than for NRLMIS for the Halloween Storm, and smaller differences than results published for JB2008 and the assimilative model HASDM. The empirical models are at present more practical for operational purposes, but the TIEGCM, developed as a research model, with a greater focus on operational use offers the potential for improved utility during stressing conditions.
Ten years of CHAMP magnetic field measurements are integrated into MFACE, a model of field-aligned currents (FACs) using empirical orthogonal functions (EOFs). EOF1 gives the basic Region-1/Region-2 pattern varying mainly with the interplanetary magnetic field Bz component. EOF2 captures separately the cusp current signature and By-related variability. Compared to existing models, MFACE yields significantly better spatial resolution, reproduces typically observed FAC thickness and intensity, improves on the magnetic local time (MLT) distribution, and gives the seasonal dependence of FAC latitudes and the NBZ current signature. MFACE further reveals systematic dependences on By, including 1) Region-1/Region-2 topology modifications around noon; 2) imbalance between upward and downward maximum current density; 3) MLT location of the Harang discontinuity. Furthermore, our procedure allows quantifying response times of FACs to solar wind driving at the bow shock nose: we obtain 20 minutes and 35-40 minutes lags for the FAC density and latitude, respectively.