The International DORIS Service (IDS): Toward maturity

DORIS is one of the four space-geodetic techniques participating in the Global Geodetic Observing System (GGOS), particularly to maintain and disseminate the Terrestrial Reference Frame as determined by International Earth rotation and Reference frame Service (IERS). A few years ago, under the umbrella of the International Association of Geodesy, a DORIS International Service (IDS) was created in order to foster international cooperation and to provide new scientific products. This paper addresses the organizational aspects of the IDS and presents some recent DORIS scientific results. It is for the first time that, in preparation of the ITRF2008, seven Analysis Centers (AC’s) contributed to derive long-term time series of DORIS stations positions. These solutions were then combined into a homogeneous time series IDS-2 for which a precision of less than 10 mm was obtained. Orbit comparisons between the various AC’s showed an excellent agreement in the radial component, both for the SPOT satellites (e.g. 0.5–2.1 cm RMS for SPOT-2) and Envisat (0.9–2.1 cm RMS), using different software packages, models, corrections and analysis strategies. There is now a wide international participation within IDS that should lead to future improvements in DORIS analysis strategies and DORIS-derived geodetic products.     

Computer modelling of a penetrator thermal sensor

The Philae lander is part of the Rosetta mission to investigate comet 67P/Churyumov-Gerasimenko. It will use a harpoon like device to anchor itself onto the surface. The anchor will perhaps reach depths of 1–2 m. In the anchor is a temperature sensor that will measure the boundary temperature as part of the MUPUS experiment. As the anchor attains thermal equilibrium with the comet ice it may be possible to extract the thermal properties of the surrounding ice, such as the thermal diffusivity, by using the temperature sensor data. The anchor is not an optimal shape for a thermal probe and application of analytical solutions to the heat equation is inappropriate. We prepare a numerical model to fit temperature sensor data and extract the thermal diffusivity. Penetrator probes mechanically compact the material immediately surrounding them as they enter the target. If the thermal properties, composition and dimensions of the penetrator are known, then the thermal properties of this pristine material may be recovered although this will be a challenging measurement. We report on investigations, using a numerical thermal model, to simulate a variety of scenarios that the anchor may encounter and how they will affect the measurement.     

Wolf–Rayet optically thick winds with Alfvén waves

The Wolf–Rayet (WR) stars are hot luminous objects which are suffering an extreme mass loss via a continuous stellar wind. The high values of mass loss rates and high terminal velocities of the WR stellar winds constitute a challenge to the theories of radiation driven winds. Several authors incorporated magnetic forces to the line driven mechanism in order to explain these characteristics of the wind. Observations indicate that the WR stellar winds may reach, at the photosphere, velocities of the order of the terminal values, which means that an important part of the wind acceleration occurs at the optically thick region. The aim of this study is to analyze a model in which the wind in a WR star begins to be accelerated in the optically thick part of the wind. We used as initial conditions stellar parameters taken from the literature and solved the energy, mass and momentum equations. We demonstrate that the acceleration only by radiative forces is prevented by the general behavior of the opacities. Combining radiative forces plus a flux of Alfvén waves, we found in the simulations a fast drop in the wind density profile which strongly reduces the extension of the optically thick region and the wind becomes optically thin too close its base. The understanding how the WR wind initiate is still an open issue.     

Survival and death of the haloarchaeon Natronorubrum strain HG-1 in a simulated martian environment

Halophilic archaea are of interest to astrobiology due to their survival capabilities in desiccated and high salt environments. The detection of remnants of salty pools on Mars stimulated investigations into the response of haloarchaea to martian conditions. Natronorubrum sp. strain HG-1 is an extremely halophilic archaeon with unusual metabolic pathways, growing on acetate and stimulated by tetrathionate. We exposed Natronorubrum strain HG-1 to ultraviolet (UV) radiation, similar to levels currently prevalent on Mars. In addition, the effects of low temperature (4, −20, and −80 °C), desiccation, and exposure to a Mars soil analogue from the Atacama desert on the viability of Natronorubrum strain HG-1 cultures were investigated. The results show that Natronorubrum strain HG-1 cannot survive for more than several hours when exposed to UV radiation equivalent to that at the martian equator. Even when protected from UV radiation, viability is impaired by a combination of desiccation and low temperature. Desiccating Natronorubrum strain HG-1 cells when mixed with a Mars soil analogue impaired growth of the culture to below the detection limit. Overall, we conclude that Natronorubrum strain HG-1 cannot survive the environment currently present on Mars. Since other halophilic microorganisms were reported to survive simulated martian conditions, our results imply that survival capabilities are not necessarily shared between phylogenetically related species.     

IDS contribution to ITRF2008

For the first time, the International DORIS Service (IDS) has produced a technique level combination based on the contributions of seven analysis centers (ACs), including the European Space Operations Center (ESOC), Geodetic Observatory Pecny (GOP), Geoscience Australia (GAU), the NASA Goddard Space Flight Center (GSFC), the Institut Géographique National (IGN), the Institute of Astronomy, Russian Academy of Sciences (INASAN, named as INA), and CNES/CLS (named as LCA). The ACs used five different software packages to process the DORIS data from 1992 to 2008, including NAPEOS (ESA), Bernese (GOP), GEODYN (GAU, GSC), GIPSY/OASIS (INA), and GINS (LCA). The data from seven DORIS satellites, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5, Envisat and Jason-1 were processed and all the analysis centers produced weekly SINEX files in either variance–covariance or normal equation format. The processing by the analysis centers used the latest GRACE-derived gravity models, forward modelling of atmospheric gravity, updates to the radiation pressure modelling to improve the DORIS geocenter solutions, denser parameterization of empirically determined drag coefficients to improve station and EOP solutions, especially near the solar maximum in 2001–2002, updated troposphere mapping functions, and an ITRF2005-derived station set for orbit determination, DPOD2005. The CATREF software was used to process the weekly AC solutions, and produce three iterations of an IDS global weekly combination. Between the development of the initial solution IDS-1, and the final solution, IDS-3, the ACs improved their analysis strategies and submitted updated solutions to eliminate troposphere-derived biases in the solution scale, to reduce drag-related degradations in station positioning, and to refine the estimation strategy to improve the combination geocenter solution. An analysis of the frequency content of the individual AC geocenter and scale solutions was used as the basis to define the scale and geocenter of the IDS-3 combination. The final IDS-3 combination has an internal position consistency (WRMS) that is 15 to 20 mm before 2002 and 8 to 10 mm after 2002, when 4 or 5 satellites contribute to the weekly solutions. The final IDS-3 combination includes solutions for 130 DORIS stations on 67 different sites of which 35 have occupations over 16 years (1993.0–2009.0). The EOPs from the IDS-3 combination were compared with the IERS 05 C04 time series and the RMS agreement was 0.24 mas and 0.35 mas for the X and Y components of polar motion. The comparison to ITRF2005 in station position shows an agreement of 6 to 8 mm RMS in horizontal and 10.3 mm in height. The RMS comparison to ITRF2005 in station velocity is at 1.8 mm/year on the East component, to 1.2 mm/year in North component and 1.6 mm/year in height.     

A Far Infrared Photometer (FIRP) for the Infrared Telescope in Space (IRTS)

We describe the design and calibration of the Far-Infrared Photometer (FIRP), one of four focal plane instruments on the Infrared Telescope in Space (IRTS). The FIRP will provide absolute photometry in four bands centered at 150, 250, 400, and 700 μm with spectral resolution λ/Δλ ≈ 3 and spatial resolution ΔΘ = 0.5 degrees. High sensitivity is achieved by using bolometric detectors operated at 300 mK in an AC bridge circuit. The closed-cycle ³He refrigerator can be recycled in orbit. A 2 K shutter provides a zero reference for each field of view. More than 10% of the sky will be surveyed during the ≈3 week mission lifetime with a sensitivity of <10⁻¹³ W·cm⁻²·sr⁻¹ per 0.5 degree pixel.     

Studies of the geocenter motion using 16-years DORIS data

An accuracy of geocenter motion estimation is strongly dependent on the geodetic network size and stations distribution over the Earth’s surface. From this point of view DORIS system has an advantage, as its ground network of beacons consists of more than 50 sites, equally distributed over the Earth’s surface. Aiming to study variations of the geocenter movements, the results of DORIS data analysis for the time span 1993.0–2009.0 (inawd06.snx series), performed at the Analysis Centre of the Institute of astronomy of the Russian Academy of Sciences, have been used. DORIS data processing was made with GIPSY/OASIS II software, developed by Jet Propulsion Laboratory and modified for DORIS data processing by Institute Géographique National. Standard deviations of stations coordinates are estimated at the level 0.5–4.0 cm (internal consistency), depending on the number of satellites used in the solution. RMS of estimated components of the DORIS satellites orbits, compared with the solutions of other IDS analysis centres, do not exceed 1–2 cm. Weekly solutions for coordinates have been transformed from free network solutions (inawd06.snx series) to a well defined terrestrial reference frame ITRF2005 with the use of seven parameters of Helmert transformation, which were examined with a view to study variations of the geocenter movements (ina05wd01.geoc time series). In order to estimate linear trend, amplitudes, periods and phases of geocenter variation a method of linear regression was applied. The evaluated amplitudes of annual variations are of the order of 5–7 mm for X and Y components and 27–29 mm for Z component. Semi-annual amplitudes are also noticeable in all components (1–34 mm for X, Y and Z components). Secular trends in the DORIS geocenter coordinates are: −1.2, −0.1 and −0.3 mm/year for X, Y and Z directions respectively.     

DORIS geodesy: A dynamic determination of geocentre location

Geoscience Australia contributed a multi-satellite, multi-year weekly time series to the International DORIS Service combined submission for the construction of International Terrestrial Reference Frame 2008 (ITRF2008). This contributing solution was extended to a study of the capability of DORIS to dynamically estimate the variation in the geocentre location. Two solutions, comprising different constraint configurations of the tracking network, were undertaken. The respective DORIS satellite orbit solutions (SPOT-2, SPOT-4, SPOT-5 and Envisat) were verified and validated by comparison with those produced at the Goddard Space Flight Center (GSFC), DORIS Analysis Centre, for computational consistency and standards. In addition, in the case of Envisat, the trajectories from the GA determined SLR and DORIS orbits were compared. The results for weekly dynamic geocentre estimates from the two constraint configurations were benchmarked against the geometric geocentre estimates from the IDS-2 combined solution. This established that DORIS is capable of determining the dynamic geocentre variation by estimating the degree one spherical harmonic coefficients of the Earth’s gravity potential. It was established that constrained configurations produced similar results for the geocentre location and consequently similar annual amplitudes. For the minimally constrained configuration Greenbelt–Kitab, the mean of the uncertainties of the geocentre location were 2.3, 2.3 and 7.6 mm and RMS of the mean uncertainties were 1.9, 1.2 and 3.5 mm for the X, Y and Z components, respectively. For GA_IDS-2_Datum constrained configuration, the mean of the uncertainties of the geocentre location were 1.7, 1.7 and 6.2 mm and RMS of the mean uncertainties were 0.9, 0.7 and 2.9 mm for the X, Y and Z components, respectively. The mean of the differences of the two DORIS dynamic geocentre solutions with respect to the IDS-2 combination were 1.6, 4.0 and 5.1 mm with an RMS of the mean 21.2, 14.0 and 31.5 mm for the Greenbelt–Kitab configuration and 4.1, 3.9 and 4.3 mm with an RMS 8.1, 9.0 and 28.6 mm for the GA_IDS-2_Datum constraint configuration. The annual amplitudes for each component were estimated to be 5.3, 10.8 and 11.0 mm for the Greenbelt–Kitab configuration and 5.3, 9.3 and 9.4 mm for the GA_IDS-2_Datum constraint configuration. The two DORIS determined dynamic geocentre solutions were compared to the SLR determined dynamic solution (which was determined from the same process of the GA contribution to the ITRF2008 ILRS combination) gave mean differences of 3.3, −4.7 and 2.5 mm with an RMS of 20.7, 17.5 and 28.0 mm for the X, Y and Z components, respectively for the Greenbelt–Kitab configuration and 1.1, −5.4 and 4.4 mm with an RMS of 9.7, 13.3 and 24.9 mm for the GA_IDS-2_Datum configuration. The larger variability is reflected in the respective amplitudes. As a comparison, the annual amplitudes of the SLR determined dynamic geocentre are 0.9, 1.0 and 6.8 mm in the X, Y and Z components. The results from this study indicate that there is potential to achieve precise dynamically determined geocentre from DORIS.     

High resolution radar clutter statistics

The generalized compound probability density function (GC-pdf) is presented for modeling high resolution radar clutter. In particular, the model is used to describe deviation of the speckle component from the Rayleigh to Weibull or other pdfs with longer tails. The GC-pdf is formed using the generalized gamma (GΓ) pdf to describe both the speckle and the modulation component of the radar clutter. The proposed model is analyzed and thermal noise is incorporated into it. The validation of the GC-pdf with real data is carried out employing the statistical moments as well as goodness-of-fit tests. A large variety of experimental data is used for this purpose. The GC-pdf outperforms the K-pdf in modeling high resolution radar clutter and reveals its structural characteristics     

Intelligent control of a planning system for astronaut training

This work intends to design, analyze and solve, from the systems control perspective, a complex, dynamic, and multiconstrained planning system for generating training plans for crew members of the NASA-led International Space Station. Various intelligent planning systems have been developed within the framework of artificial intelligence. These planning systems generally lack a rigorous mathematical formalism to allow a reliable and flexible methodology for their design, modeling, and performance analysis in a dynamical, time-critical, and multiconstrained environment. Formulating the planning problem in the domain of discrete-event systems under a unified framework such that it can be modeled, designed, and analyzed as a control system will provide a self-contained theory for such planning systems. This will also provide a means to certify various planning systems for operations in the dynamical and complex environments in space. The work presented here completes the design, development, and analysis of an intricate, large-scale, and representative mathematical formulation for intelligent control of a real planning system for Space Station crew training. This planning system has been tested and used at NASA-Johnson Space Center