Advanced multilateration theory, software development, and data processing: The MICRODOT system

The process of geometric parameter estimation to accuracies of one centimeter, i.e., multilateration, is defined and applications are listed. A brief functional explanation of the theory is presented. Next, various multilateration systems are described in order of increasing system complexity. Expected systems accuracy is discussed from a general point of view and a summary of the errors is listed. An outline of the design of a software processing system for multilateration, called MICRODOT, is presented next. The links of this software, which can be used for multilateration data simulations or operational data reduction, are examined on an individual basis. Functional flow diagrams are presented to aid in understanding the software capability. MICRODOT capability is described with respect to vehicle configurations, interstation coordinate reduction, geophysical parameter estimation, and orbit determination. Numerical results obtained from MICRODOT via data simulations are displayed both for hypothetical and real world vehicle/station configurations such as used in the GEOS-3 Project. These simulations show the inherent power of the multilateration procedure

Three-D multilateration: A precision geodetic measurement system

A technique of satellite geodesy for determining the relative three dimensional coordinates of ground stations within one centimeter over baselines of 20 to 10,000 kilometers is discussed. The system is referred to as 3-D Multilateration and has applications in earthquake hazard assessment, precision surveying, plate tectonics, and orbital mechanics. The accuracy is obtained by using pulsed lasers to obtain simultaneous slant ranges between several ground stations and a moving retroreflector with known trajectory for aiming the lasers

Locating the LCROSS Impact Craters

The Lunar CRater Observations and Sensing Satellite (LCROSS) mission impacted a spent Centaur rocket stage into a permanently shadowed region near the lunar south pole. The Sheperding Spacecraft (SSC) separated \sim9 hours before impact and performed a small braking maneuver in order to observe the Centaur impact plume, looking for evidence of water and other volatiles, before impacting itself. This paper describes the registration of imagery of the LCROSS impact region from the mid- and near-infrared cameras onboard the SSC, as well as from the Goldstone radar. We compare the Centaur impact features, positively identified in the first two, and with a consistent feature in the third, which are interpreted as a 20 m diameter crater surrounded by a 160 m diameter ejecta region. The images are registered to Lunar Reconnaisance Orbiter (LRO) topographical data which allows determination of the impact location. This location is compared with the impact location derived from ground-based tracking and propagation of the spacecraft's trajectory and with locations derived from two hybrid imagery/trajectory methods. The four methods give a weighted average Centaur impact location of -84.6796\circ, -48.7093\circ, with a 1{\sigma} un- certainty of 115 m along latitude, and 44 m along longitude, just 146 m from the target impact site. Meanwhile, the trajectory-derived SSC impact location is -84.719\circ, -49.61\circ, with a 1{\sigma} uncertainty of 3 m along the Earth vector and 75 m orthogonal to that, 766 m from the target location and 2.803 km south-west of the Centaur impact. We also detail the Centaur impact angle and SSC instrument pointing errors. Six high-level LCROSS mission requirements are shown to be met by wide margins. We hope that these results facilitate further analyses of the LCROSS experiment data and follow-up observations of the impact region.Comment: Accepted for publication in Space Science Review. 24 pages, 9 figure

A possible new test of general relativity with Juno

The expansion in multipoles of the gravitational potential of a rotating body affects the orbital motion of a test particle orbiting it with long-term perturbations both at a classical and at a relativistic level. In this preliminary sensitivity analysis, we show that, for the first time, the J2 c^-2 effects could be measured by the ongoing Juno mission in the gravitational field of Jupiter during its yearlong science phase (10 November 2016-5 October 2017) thanks to its high eccentricity (e=0.947) and to the huge oblateness of Jupiter (J2=1.47 10^-2). The semi-major axis a and the perijove \omega\ of Juno are expected to be shifted by \Delta a =700-900 m and \Delta\omega = 50-60 milliarcseconds, respectively, over 1-2 yr. A numerical analysis shows also that the expected J2c^-2 range-rate signal for Juno should be as large as 280 microns per second during a typical 6 h pass at its closest approach. Independent analyses previously performed by other researchers about the measurability of the Lense-Thirring effect showed that the radio science apparatus of Juno should reach an accuracy in Doppler range-rate measurements of 1-5 microns per second over such passes. The range-rate signature of the classical even zonal perturbations is different from the 1PN one. Thus, further investigations, based on covariance analyses of simulated Doppler data and dedicated parameters estimation, are worth of further consideration. It turns out that the J2 c^-2 effects cannot be responsible of the flyby anomaly in the gravitational field of the Earth. A dedicated spacecraft in a 6678 km X 57103 km polar orbit would experience a geocentric J2 c^-2 range-rate shift of 0.4 mm s^-1.Comment: LaTex2e, 16 pages, no tables, 2 figures, 39 references. Version published in Classical and Quantum Gravity (CQG

LIBRA: An inexpensive geodetic network densification system

A description is given of the Libra (Locations Interposed by Ranging Aircraft) system, by which geodesy and earth strain measurements can be performed rapidly and inexpensively to several hundred auxiliary points with respect to a few fundamental control points established by any other technique, such as radio interferometry or satellite ranging. This low-cost means of extending the accuracy of space age geodesy to local surveys provides speed and spatial resolution useful, for example, for earthquake hazards estimation. Libra may be combined with an existing system, Aries (Astronomical Radio Interferometric Earth Surveying) to provide a balanced system adequate to meet the geophysical needs, and applicable to conventional surveying. The basic hardware design was outlined and specifications were defined. Then need for network densification was described. The following activities required to implement the proposed Libra system are also described: hardware development, data reduction, tropospheric calibrations, schedule of development and estimated costs

First In-orbit Experience of TerraSAR-X Flight Dynamics Operations

TerraSAR-X is an advanced synthetic aperture radar satellite system for scientific and commercial applications that is realized in a public-private partnership between the German Aerospace Center (DLR) and the Astrium GmbH. TerraSAR-X was launched at June 15, 2007 on top of a Russian DNEPR-1 rocket into a 514 km sun-synchronous dusk-dawn orbit with an 11-day repeat cycle and will be operated for a period of at least 5 years during which it will provide high resolution SAR-data in the X-band. Due to the objectives of the interferometric campaigns the satellite has to comply to tight orbit control requirements, which are formulated in the form of a 250 m toroidal tube around a pre-flight determined reference trajectory. The acquisition of the reference orbit was one of the main and key activities during the Launch and Early Orbit Phase (LEOP) and had to compensate for both injection errors and spacecraft safe mode attitude control thruster activities. The paper summarizes the activities of GSOC flight dynamics team during both LEOP and early Commissioning Phase, where the main tasks have been 1) the first-acquisition support via angle-tracking and GPS-based orbit determination, 2) maneuver planning for target orbit acquisition and maintenance, and 3) precise orbit and attitude determination for SAR processing support. Furthermore, a presentation on the achieved results and encountered problems will be addressed

Improved solution accuracy for TDRSS-based TOPEX/Poseidon orbit determination

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using a batch-least-squares estimator available in the Goddard Trajectory Determination System (GTDS) and an extended Kalman filter estimation system to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. GTDS is the operational orbit determination system used by the FDD in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. The extended Kalman filter was implemented in an orbit determination analysis prototype system, closely related to the Real-Time Orbit Determination System/Enhanced (RTOD/E) system. In addition, the Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generated an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the geodynamics (GEODYN) orbit determination system with laser ranging and Doppler Orbitography and Radiopositioning integrated by satellite (DORIS) tracking measurements. The TOPEX/Poseidon trajectories were estimated for November 7 through November 11, 1992, the timeframe under study. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch-least-squares solutions were assessed based on the solution residuals, while the sequential solutions were assessed based on primarily the estimated covariances. The batch-least-squares and sequential orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 2 meters for the batch-least-squares and less than 13 meters for the sequential estimation solutions. After the sequential estimation solutions were processed with a smoother algorithm, position differences with POD orbit solutions of less than 7 meters were obtained. The differences among the POD, GTDS, and filter/smoother solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail

The GLAS Algorithm Theoretical Basis Document for Precision Orbit Determination (POD)

The Geoscience Laser Altimeter System (GLAS) was the sole instrument for NASA's Ice, Cloud and land Elevation Satellite (ICESat) laser altimetry mission. The primary purpose of the ICESat mission was to make ice sheet elevation measurements of the polar regions. Additional goals were to measure the global distribution of clouds and aerosols and to map sea ice, land topography and vegetation. ICESat was the benchmark Earth Observing System (EOS) mission to be used to determine the mass balance of the ice sheets, as well as for providing cloud property information, especially for stratospheric clouds common over polar areas. The GLAS instrument operated from 2003 to 2009 and provided multi-year elevation data needed to determine changes in sea ice freeboard, land topography and vegetation around the globe, in addition to elevation changes of the Greenland and Antarctic ice sheets. This document describes the Precision Orbit Determination (POD) algorithm for the ICESat mission. The problem of determining an accurate ephemeris for an orbiting satellite involves estimating the position and velocity of the satellite from a sequence of observations. The ICESatGLAS elevation measurements must be very accurately geolocated, combining precise orbit information with precision pointing information. The ICESat mission POD requirement states that the position of the instrument should be determined with an accuracy of 5 and 20 cm (1-s) in radial and horizontal components, respectively, to meet the science requirements for determining elevation change

Historical forest biomass dynamics modelled with Landsat spectral trajectories

Acknowledgements National Forest Inventory data are available online, provided by Ministerio de Agricultura, Alimentación y Medio Ambiente (España). Landsat images are available online, provided by the USGS.Peer reviewedPostprin