Using GPS multipath to measure soil moisture fluctuations: initial results

Measurements of soil moisture are important for studies of climate and weather forecasting, flood prediction, and aquifer recharge studies. Although soil moisture measurement networks exist, most are sparsely distributed and lack standardized instrumentation. Measurements of soil moisture from satellites have extremely large spatial footprints (40–60 km). A methodology is described here that uses existing networks of continuously-operating GPS receivers to measure soil moisture fluctuations. In this technique, incoming signals are reflected off and attenuated by the ground before reception by the GPS receiver. These multipath reflections directly affect signal-to-noise ratio (SNR) data routinely collected by GPS receivers, creating amplitude variations that are a function of ground reflectivity and therefore soil moisture content. After describing this technique, multipath reflection amplitudes at a GPS site in Tashkent, Uzbekistan are compared to estimates of soil moisture from the Noah land surface model. Although the GPS multipath amplitudes and the land surface model are uncalibrated, over the 70-day period studied, they both rise sharply following each rainfall event and slowly decrease over a period of ∼10 days.     

Bandlimited implicit Runge–Kutta integration for Astrodynamics

We describe a new method for numerical integration, dubbed bandlimited collocation implicit Runge–Kutta (BLC-IRK), and compare its efficiency in propagating orbits to existing techniques commonly used in Astrodynamics. The BLC-IRK scheme uses generalized Gaussian quadratures for bandlimited functions. This new method allows us to use significantly fewer force function evaluations than explicit Runge–Kutta schemes. In particular, we use a low-fidelity force model for most of the iterations, thus minimizing the number of high-fidelity force model evaluations. We also investigate the dense output capability of the new scheme, quantifying its accuracy for Earth orbits. We demonstrate that this numerical integration technique is faster than explicit methods of Dormand and Prince 5(4) and 8(7), Runge–Kutta–Fehlberg 7(8), and approaches the efficiency of the 8th-order Gauss–Jackson multistep method. We anticipate a significant acceleration of the scheme in a multiprocessor environment.     

Parallel initial orbit determination using angles-only observation pairs

We propose a type of admissible-region analysis for track initiation in multi-satellite problems when angles are the primary observable. Pairs of optical observations are used to calculate candidate orbits via a Lambert solver by hypothesizing range values. The method is attractive because it allows multiple levels of parallelization of the track-initiation process. Orbital element partitions are introduced to divide the admissible region into smaller search spaces to be processed on individual computer nodes. For a specified rectangular partition in the space of orbital elements, constraints are developed to bound the values of range that will lead to initial orbit hypotheses (data association hypotheses) associated with that partition. These bounds allow us to parallelize the generation of candidate orbits, because each element-space partition can be handled independently of the others. Several constraints are developed and shown to limit the range pair hypotheses effectively to the constrained admissible region based on the orbital element partitions. Examples are provided to highlight the topology of the proposed constraints.     

Satellite clock bias estimation for iGPS

The High Integrity GPS program seeks to provide enhanced navigation performance by combining conventional GPS with a communications and ranging broadcast from the Iridium^® Communications System. Through clock and message aiding, it would enable existing GPS receivers to acquire and track in more challenging environments. As is the case for standard GPS, accurate and precise timing is key to performance. An approach is presented for estimating the bias of each Iridium satellite clock using satellite-to-ground and satellite-to-satellite measurements. The satellite clock bias estimates are based on a Kalman filter that incorporates code-type observations from the measurements at 10 s intervals. Filter parameters are set based on the expected behavior of the clocks, allowing for discontinuous bias and frequency adjustments due to ground commands. Typical results show the current filter to be accurate to within 200 ns while always meeting the initial system specification of half a microsecond.