Simulation and analysis of the revised Huygens probe entry and descent trajectory and radio link modelling

The Cassini spacecraft will arrive at Saturn in 2004 carrying the Huygens probe. The beginning of the Cassini tour at Saturn has been redesigned to achieve a different relative orbiter/probe geometry in order to compensate for the probe relay receiver design flaw that was discovered during tests in February 2000. This paper presents a numerical simulation of the Huygens atmospheric entry and descent trajectory and the Cassini flyby trajectory during the probe mission. A variety of parameters that are crucial for the probe system and its scientific payload have been calculated and analyzed together with an assessment of their uncertainties. Furthermore the orbiter/probe relay link was simulated in order to assess any potential data loss on the basis of an analytical model of the actual Huygens receiver onboard the Cassini spacecraft. The redesigned Cassini/Huygens mission satisfies all science and engineering requirements and assures the best possible radio link for the entire nominal mission duration.     

Huygens probe entry dynamic model and accelerometer data analysis

During the first phase of Huygens arrival into Titan's atmosphere the probe is subjected to gravitational and aerodynamic forces in aerodynamic hypersonic regime. Atmospheric drag exerts a strong deceleration on the capsule measured by Huygens atmospheric structure instrument (HASI) servo accelerometer. A 6 DOF (Degree of Freedom) model of the Huygens probe entry dynamics has been developed and used for data analysis. The accelerometer data are analysed and the model allows the retrieval of dynamics information of Huygens probe from 1545 km altitude down to end of the entry phase. Probe's initial conditions (velocity and position) were refined to match the measured deceleration profile resulting in a different altitude at interface epoch with respect to those of the Cassini Navigation Team. Velocity and position of probe at interface epoch are compatible with those used by Descent Trajectory Working Group (DTWG).Measurements acquired before atmosphere detection are used to estimate probe's angular rate, bound attitude and characterise the angle of attack profile which results to be lower than 4^° during the whole entry. Probe's spin calculated (6.98 RPM) is slightly different with respect to DTWG of 7.28 RPM but considering a 2% error in the Inertia matrix these results are inside the 1-σ error band.     

Huygens probe entry trajectory and attitude estimated simultaneously with Titan atmospheric structure by Kalman filtering

Space probes entering planetary atmospheres are used for in situ study of their physical structures. During the entry phase aerodynamic forces exerted on the probe depend on atmospheric density. As a consequence accelerations measured by on-board sensors can be used to derive probe trajectory as well as atmospheric density, pressure and temperature profiles. In this work acceleration data acquired by the Huygens Atmospheric Structure Instrument (HASI) have been used to reconstruct the probe trajectory and the Titan's atmospheric structure from down to of altitude. An accurate six degree of freedom model of Huygens during the entry phase has been developed and a new reconstruction technique based on Kalman filtering is presented. This technique estimates simultaneously the probe trajectory, the attitude profile consistent with measured data and the atmospheric density, pressure and temperature.     

Huygens’ entry and descent through Titan's atmosphere—Methodology and results of the trajectory reconstruction

The European Space Agency's Huygens probe separated from the NASA Cassini spacecraft on 25 December 2004, after having been attached for a 7-year interplanetary journey and three orbits around Saturn. The probe reached the predefined NASA/ESA interface point on 14 January 2005 at 09:05:52.523 (UTC) and performed a successful entry and descent sequence. The probe softly impacted on Titan's surface on the same day at 11:38:10.77 (UTC) with a speed of about 4.54 m/s. The probe entry and descent trajectory was reconstructed from the estimated initial state vector provided by the Cassini Navigation team, the probe housekeeping data, and measurements from the scientific payload. This paper presents the methodology and discuss the results of the reconstruction effort. Furthermore the probe roll rate was reconstructed prior to the main entry phase deceleration pulse and throughout the entire descent phase under the main and drogue parachute.     

Reconstruction of the trajectory of the Huygens probe using the Huygens Atmospheric Structure Instrument (HASI)

The Huygens probe returned scientific measurements from the atmosphere and surface of Titan on 14 January 2005. Knowledge of the trajectory of Huygens is necessary for scientific analysis of those measurements. We use measurements from the Huygens Atmospheric Structure Instrument (HASI) to reconstruct the trajectory of Huygens during its mission. The HASI Accelerometer subsystem measured the axial acceleration of the probe with errors of 3E−6 m s⁻², the most accurate measurements ever made by an atmospheric structure instrument on another planetary body. The atmosphere was detected at an altitude of 1498 km. Measurements of the normal acceleration of the probe, which are important for determining the probe's attitude during hypersonic entry, were significantly less accurate and limited by transverse sensitivity of the piezo sensors. Peak acceleration of 121.2 m s⁻² occurred at 234.9 km altitude. The parachute deployment sequence started at 157.1 km and a speed of 342.1 m s⁻¹. Direct measurements of pressure and temperature began shortly afterwards. The measured accelerations and equations of motion have been used to reconstruct the trajectory prior to parachute deployment. Measured pressures and temperatures, together with the equation of hydrostatic equilibrium and the equation of state, have been used to reconstruct the trajectory after parachute deployment. Uncertainties in the entry state of Huygens at the top of the atmosphere are significant, but can be reduced by requiring that the trajectory and atmospheric properties be continuous at parachute deployment.     

Results from the Huygens probe on Titan

The Cassini–Huygens mission, comprising the NASA Saturn Orbiter and the ESA Huygens Probe, arrived at Saturn in late June 2004. The Huygens probe descended under parachute in Titan’s atmosphere on 14 January 2005, 3 weeks after separation from the Orbiter. We discuss here the breakthroughs that the Huygens probe, in conjunction with the Cassini spacecraft, brought to Titan science. We review the achievements ESA’s Huygens probe put forward and the context in which it operated. The findings include new localized information on several aspects of Titan science: the atmospheric structure and chemical composition; the aerosols distribution and content; the surface morphology and composition at the probe’s landing site; the winds, the electrical properties, and the implications on the origin and evolution of the satellite.     

Huygens/HASI 2002 balloon test campaign: Probe trajectory and atmospheric vertical profiles reconstruction

In the framework of the activities going on in preparation for the mission of the Huygens probe in Titan's atmosphere (January 2005), the Huygens Atmospheric Structure Instrument (HASI) team scheduled and performed several balloon campaigns to test the HASI sensors’ performance in flight conditions in the Earth's atmosphere. In particular, pressure conditions reached during each test are similar to those expected in Titan's lower atmosphere. A 1:1 scaled mock-up of the Huygens probe was launched with a stratospheric balloon in 2001 (Br. Assoc. Adv. Sci. 33 (2001) 1109) and in 2002 (Br. Assoc. Adv. Sci. 34 (2002) 911; Adv. Space. Sci. (2003)) from the G. Broglio base of the Italian Space Agency, located in Trapani Milo (Sicily). In both cases the mock-up was dropped from an altitude higher than 27 and , respectively, and recovered on the ground after a parachuted descent. In this paper, we describe the results obtained in reconstructing (i) the probe descent trajectory and (ii) the profiles of the physical quantities characterizing the Earth's atmosphere, on the basis of a complete analysis of the data obtained during the HASI 2002 balloon flight experiment. Using temperature and pressure measurements, we are able to reach an accuracy of the order of 0.5% on the altitude reconstruction during the descent. We validate both the models used for trajectory reconstruction and to check the sensors’ performance. We describe the problems faced in determining the Huygens probe descent trajectory in Titan's atmosphere focusing our discussion on the critical aspects of the descent reconstruction (such as the uncertainties due to measurement errors, limited knowledge of the atmospheric composition, etc.) and the validity of the adopted assumptions.     

New laboratory measurements of CH4 in Titan's conditions and a reanalysis of the DISR near-surface spectra at the Huygens landing site

Laboratory spectra of methane-nitrogen mixtures have been recorded in the near-infrared range (1.0-1.65 μm) in conditions similar to Titan's near surface, to facilitate the interpretation of the DISR/DLIS (DISR—Descent Imager/Spectral Radiometer) spectra taken during the last phase of the descent of the Huygens Probe, when the surface was illuminated by a surface-science lamp. We used a 0.03 cm⁻¹ spectral resolution, adequate to resolve the lines at high pressure (p_N2∼1.5 bar). By comparing the laboratory spectra with synthetic calculations in the well-studied ν₂+2ν₃ band (7515-7620 cm⁻¹), we determine a methane absorption column density of 178±20 cm atm and a temperature of 118±10 K in our experiment. From this, we derive the methane absorption coefficients over 1.0-1.65 μm with a 0.03 cm⁻¹ sampling, allowing for the extrapolation of the results to any other methane column density under the relevant pressure and temperature conditions. We then revisit the calibration and analysis of the Titan “lamp-on” DLIS spectra. We infer a 5.1±0.8% methane-mixing ratio in the first 25 m of Titan's atmosphere. The CH₄ mixing ratio measured 90 s after landing from a distance of 45 cm is found to be 0.92±0.25 times this value, thus showing no post-landing outgassing of methane in excess of ∼20%. Finally, we determine the surface reflectivity as seen between 25 m and 45 cm and find that the 1500 nm absorption band is deeper in the post-landing spectrum as compared to pre-landing.     

Electron conductivity and density profiles derived from the mutual impedance probe measurements performed during the descent of Huygens through the atmosphere of Titan

During the descent of the Huygens probe through the atmosphere of Titan, on January 14th, 2005, the permittivity, waves and altimetry (PWA) subsystem, a component of the Huygens atmospheric structure instrument (HASI), detected an ionized layer at altitudes around 63 km with two different instruments, the relaxation probe (RP) and the mutual impedance probe (MIP). A very detailed analysis of both data sets is required, in order to correct for environmental effects and compare the two independent estimates of the electrical conductivity. The present work is dedicated to the MIP data analysis. New laboratory tests have been performed to validate or improve the available calibration results. Temperature effects have been included and numerical models of the MIP sensors and electric circuitry have been developed to take into account the proximity of the Huygens probe body. The effect of the vertical motion of the vessel in the ionized atmosphere is estimated in both analytical and numerical ways. The peculiar performance of the instrument in the altitude range 100–140 km is scrutinized. The existence of a prominent ionized layer, and of enhancements in the conductivity and electron density profiles at 63 km, are discussed in the light of previous theoretical predictions.     

A technique to determine the mean molecular mass of a planetary atmosphere using pressure and temperature measurements made by an entry probe: Demonstration using Huygens data

It is possible to determine the mean molecular mass of a planetary atmosphere using pressure and temperature measurements made by an entry probe descending at terminal velocity. The descent trajectory of an entry probe can be determined from pressure, temperature, and mean molecular mass data. This technique offers redundancy for large entry probes in the event of a mass spectrometer failure and increases the potential scientific yield of small entry probes that do not carry mass spectrometers. This technique is demonstrated on Huygens atmospheric structure instrument (HASI) data from Titan. Accurate knowledge of entry probe and parachute drag coefficients is required for this technique to be useful.     

The photochemical formation of a titan haze analog. Structural analysis by x-ray photoelectron and infrared spectroscopy

The photochemical flow reactor (D.W. Clarke et al., 2000, Icarus 147, 282-291) has been modified to minimize the incorporation of oxygen and other impurities in the photoproducts. A mixture of gases that approximate their mixing ratios on Titan (N₂, CH₄, H₂, C₂H₂, C₂H₄, and HC₃N) (0.98, 0.018, 0.002, 3.5 × 10⁻⁴, 3 × 10⁻⁴, 1.7 × 10⁻⁵, respectively) was irradiated in the flow photochemical reactor using a 185-nm source to give a Titan haze analog as a solid product. X-ray photoelectron spectroscopy (XPS) gave a composition of 93.3% C, 5.3% N, and 1.4% O. Of the 93.3% carbon, high-resolution XPS revealed that 81.2% was present as CH, CC, and CC groups, 12.1% may be CO, CN, CN, CN, and/or CN groups, 5.3% as a CN group. The peak for N was symmetrical and was assigned to the CN while that for oxygen was assigned to the CO and/or the CO group. Some of these assignments were confirmed by FTIR spectroscopy. The polymeric product had a C:N ratio of 17.6, which is significantly greater than that for Titan haze analogs prepared in discharge reactions. When the polymer was exposed to air for seven days the oxygen content increased by 6% along with an increase in the infrared absorption at 1710 cm⁻¹ assigned to the CO group of a ketone. The oxidation is attributed to the reaction of oxygen with free radicals trapped in the polymer matrix. It is proposed that the photochemical initiation of Titan haze formation from compounds formed from starting materials formed high in Titan’s atmosphere is a more plausible model than haze formed in reactions initiated by solely by discharges. These data will be helpful in the interpretation of the data returned from the Huygens probe of the Cassini mission.     

Coupling photochemistry with haze formation in Titan's atmosphere, Part II: Results and validation with Cassini/Huygens data

The new one-dimensional radiative-convective/photochemical/microphysical model described in Part I is applied to the study of Titan's atmospheric processes that lead to haze formation. Our model generates the haze structure from the gaseous species photochemistry. Model results are presented for the species vertical concentration profiles, haze formation and its radiative properties, vertical temperature/density profiles and geometric albedo. These are validated against Cassini/Huygens observations and other ground-based and space-borne measurements. The model reproduces well most of the latest measurements from the Cassini/Huygens instruments for the chemical composition of Titan's atmosphere and the vertical profiles of the observed species. For the haze production we have included pathways that are based on pure hydrocarbons, pure nitriles and hydrocarbon/nitrile copolymers. From these, the nitrile and copolymer pathways provide the stronger contribution, in agreement with the results from the ACP instrument, which support the incorporation of nitrogen in the pyrolized haze structures. Our haze model reveals a new second major peak in the vertical profile of haze production rate between 500 and 900 km. This peak is produced by the copolymer family used and has important ramifications for the vertical atmospheric temperature profile and geometric albedo. In particular, the existence of this second peak determines the vertical profile of haze extinction. Our model results have been compared with the DISR retrieved haze extinction profiles and are found to be in very good agreement. We have also incorporated in our model heterogeneous chemistry on the haze particles that converts atomic hydrogen to molecular hydrogen. The resultant H₂ profile is closer to the INMS measurements, while the vertical profile of the diacetylene formed is found to be closer to that of the CIRS profile when this heterogenous chemistry is included.     

Biomization and quantitative climate reconstruction techniques in northwestern Mexico—With an application to four Holocene pollen sequences

New paleovegetation and paleoclimatic reconstructions from the Sierra Madre Occidental (SMO) in northwestern Mexico are presented. This work involves climate and biome reconstruction using Plant Functional Types (PFT) assigned to pollen taxa. We used fossil pollen data from four Holocene peat bogs located at different altitudes (1500‑2000 m) at the border region of Sonora and Chihuahua at around 28° N latitude (Ortega-Rosas, C.I. 2003. Palinología de la Ciénega de Camilo: datos para la historia de la vegetación y el clima del Holoceno medio y superior en el NW de la Sierra Madre Occidental, Sonora, Mexico. Master Thesis, Universidad Nacional Autónoma de México, México D.F.; Ortega-Rosas, C.I., Peñalba, M.C., Guiot, J. Holocene altitudinal shifts in vegetation belts and environmental changes in the Sierra Madre Occidental, Northwestern Mexico. Submitted for publication of Palaeobotany and Palynology). The closest modern pollen data come from pollen analysis across an altitudinal transect from the Sonoran Desert towards the highlands of the temperate SMO at the same latitude (Ortega-Rosas, C.I. 2003. Palinología de la Ciénega de Camilo: datos para la historia de la vegetación y el clima del Holoceno medio y superior en el NW de la Sierra Madre Occidental, Sonora, Mexico. Master Thesis, Universidad Nacional Autónoma de México, México D.F.). An additional modern pollen dataset of 400 sites across NW Mexico and the SW United States was compiled from different sources (Davis, O.K., 1995. Climate and vegetation pattern in surface samples from arid western U.S.A.: application to Holocene climatic reconstruction. Palynology 19, 95–119, North American Pollen Database, Latin-American Pollen Database, personal data, and different scientific papers). For the biomization method (Prentice, I.C., Guiot, J., Huntley, B., Jolly, D., Cheddadi, R., 1996. Reconstructing biomes from paleoecological data: a general method and its application to European pollen data at 0 and 6 ka. Climate Dynamics 12, 185–194), we modified the pollen-PFT and PFT-biomes assignation of Thompson and Anderson (Thompson, R.S., Anderson, K.H., 2000. Biomes of western North America at 18,000; 6000 and 0 14C yr BP reconstructed from pollen and packrat midden data. Journal of Biogeography 27, 555–584) for a better representation of the modern vegetation of NW Mexico. The biome reconstruction method was validated with the modern pollen sites and applied to the fossil sites. Our results show that, during the early Holocene, a cool conifer forest extended at least down to 1700 m, while today this biome is present above 2000 m in the Chihuahua state. The Younger Dryas event was recorded in one site with cold and dry conditions. The reconstructed annual temperature for this period was 3°–6 °C colder than today, and annual precipitation was 250 mm lower than at present (900 mm/yr). The middle Holocene after 9200 cal yr BP was marked by a warming trend, reaching temperatures 2 °C warmer than today at 7000 cal yr BP, and by the installation of a warm mixed forest, the present day biome, at 1700 m elevation, while at higher elevations (1900 m) the cool conifer forest was still present. Summer precipitation was 200 mm/yr above the early Holocene values, suggesting that monsoon-like conditions strengthened since 9200 cal yr BP at this region. During the last 4000 yr, the same warm mixed forest was reconstructed below 1700 m and a conifer forest above 1700 m. A great variability of vegetation and climate patterns was recorded for the last 3000 yr particularly at high elevation sites, where warming and cooling trends would be coeval of the Medieval warm period and Little Ice Age, likely related to ENSO variability.     

Extraterrestrial samples from low Earth orbits: techniques for their collection and analysis

In the context of dust samples collections in space, the COMET experiment (Collecte en Orbite de Matière ExtraTerrestre) was proposed for the first time in 1982. The idea of such an experiment was to collect grains with identified cometary parent body, instead of mixing all extraterrestrial contributions present in low Earth orbit. It was thus proposed to install collectors inside hermetic boxes, to have these boxes mounted outside a space station, orbiting the Earth and to have the capability of choosing the time and duration of the collection. Since 1985, the COMET experiment has been exposed three times to space (COMET-1, in October 1985 during the encounter of the Earth with the Draconid meteor stream; the EUROMIR-95 instrument, exposing collectors, during the crossing by the Earth of the Orionid meteor stream associated to comet P/Halley and, in November 1998, during the crossing by the Earth of the Leonid meteor stream associated to comet Temple-Tuttle, COMET-99). Specific collection techniques, and corresponding analytical procedures have been developed. The collected particles are the only ones accessible in the laboratory with a known cometary origin, before the return to Earth (2006) of the Stardust mission, which will collect cometary grains in the tails of comet Wild 2. Such a challenge justifies the tremendous efforts brought into play, and that are summarized here.     

Phytolith analysis and paleoenvironmental reconstruction from Lake Poukawa Core, Hawkes Bay, New Zealand

As part of a multiproxy investigation, phytoliths were extracted from sediments in a 197-m core in Hawkes Bay, New Zealand. They provide a continuous vegetation–climate record spanning the time period from at least the last interglacial (marine oxygen isotope stage 5) to the present. The phytolith record demonstrates that grass/cyperaceae grew during warmer periods and woody taxa dominated the site during colder periods.During the present interglacial, the Poukawa basin is occupied by a shallow lake surrounded by an extensive fen. During colder–drier periods, the floor of the basin dried out and woody taxa occupied the basin floor. This contrasts with the pollen record, which demonstrates a converse pattern. The apparent discrepancy reflects the purely local provenance of the phytolith assemblage.Significant changes in phytolith assemblages occur at the same depth as major tephras, indicating a sharp decline in trees and shrubs and a surge in grass and cyperaceae. A series of successional changes follow each major tephra fall. Initially, the woody taxa are killed off and replaced by grass and cyperaceae that rapidly colonise the fresh surface. Trees and or shrubs succeed the grass and cyperaceae after a significant lag.     

New confirmation of image-processing techniques for astrometry of Saturn and its satellites

The image-processing techniques used by Peng et al. are further improved to measure precisely the positions of Saturn and its satellites. 495 CCD images taken with the 1-m telescope at the Yunnan Observatory during the years 2002–2004 are processed with these techniques. These measured pixel positions are compared to their theoretical positions computed with the ephemerides of TASS1.7 for the satellites and JPL DE405 for Saturn itself. Analysis of the data for the intersatellite positions among four bright Saturnian satellites (S3–S6) and for Saturn–satellite (i.e. Saturn–Titan) positions shows that these measured positions have the same dispersions, i.e. about 0.05 and 0.06 arcsec in right ascension and declination, respectively. However, for the fainter satellites, Enceladus and Mimas, poorer residuals up to 0.1 and 0.2 arcsec, respectively, in both directions are found mainly due to their small separations from the primary planet and short exposure time in order to obtain useful images of Saturn.     

A catalogue and analysis of local galaxy ages and metallicities

We have assembled a catalogue of relative ages, metallicities and abundance ratios for about 150 local galaxies in field, group and cluster environments. The galaxies span morphological types from cD and ellipticals, to late-type spirals. Ages and metallicities were estimated from high-quality published spectral line indices using Worthey & Ottaviani (1997) single stellar population evolutionary models.
The identification of galaxy age as a fourth parameter in the fundamental plane ( Forbes, Ponman & Brown 1998 ) is confirmed by our larger sample of ages. We investigate trends between age and metallicity, and with other physical parameters of the galaxies, such as ellipticity, luminosity and kinematic anisotropy. We demonstrate the existence of a galaxy age–metallicity relation similar to that seen for local galactic disc stars, whereby young galaxies have high metallicity, while old galaxies span a large range in metallicities.
We also investigate the influence of environment and morphology on the galaxy age and metallicity, especially the predictions made by semi-analytic hierarchical clustering models (HCM). We confirm that non-cluster ellipticals are indeed younger on average than cluster ellipticals as predicted by the HCM models. However we also find a trend for the more luminous galaxies to have a higher [Mg/Fe] ratio than the lower luminosity galaxies, which is opposite to the expectation from HCM models.