Doppler Orbit Determination of Deep Space Probes by the Bernese GNSS Software: First Results of the Combined Orbit Determination from DSN and Inter-Satellite Ka-Band Data from the Grail Mission

Navigation of deep space probes is most commonly operated using the spacecraft Doppler tracking technique. Orbital parameters are determined from a series of repeated measurements of the frequency shift of a microwave carrier over a given integration time. Currently, both ESA and NASA operate antennas at several sites around the world to ensure the tracking of deep space probes. Just a small number of software packages are nowadays used to process Doppler observations. The Astronomical Institute of the University of Bern (AIUB) has recently started the development of Doppler data processing capabilities within the Bernese GNSS Software. This software has been extensively used for Precise Orbit Determination of Earth orbiting satellites using GPS data collected by on-board receivers and for subsequent determination of the Earth gravity field. In this paper, we present the currently achieved status of the Doppler data modeling and orbit determination capabilities in the Bernese GNSS Software using GRAIL data. In particular we will focus on the implemented orbit determination procedure used for the combined analysis of Doppler and intersatellite Ka-band data. We show that even at this earlier stage of the development we can achieve an accuracy of few mHz on two-way S-band Doppler observation and of 2 µm/s on KBRR data from the GRAIL primary mission phase.

Kinetic simulation of neutral∕ionized gas and electrically charged dust in the coma of comet 67P∕Churyumov-Gerasimenko

The cometary coma is a unique phenomenon in the solar system being a planetary atmosphere influenced by little or no gravity. As a comet approaches the sun, the water vapor with some fraction of other gases sublimate, generating a cloud of gas, ice and other refractory materials (rocky and organic dust) ejected from the surface of the nucleus. Sublimating gas molecules undergo frequent collisions and photochemical processes in the near‐nucleus region. Owing to its negligible gravity, comets produce a large and highly variable extensive dusty coma with a size much larger than the characteristic size of the cometary nucleus. The Rosetta spacecraft is en route to comet 67P/Churyumov‐Gerasimenko for a rendezvous, landing, and extensive orbital phase beginning in 2014. Both, interpretation of measurements and safety consideration of the spacecraft require modeling of the comet’s dusty gas environment. In this work we present results of a numerical study of multispecies gaseous and electrically charged dust environment of comet Chyuryumov‐Gerasimenko. Both, gas and dust phases of the coma are simulated kinetically. Photolytic reactions are taken into account. Parameters of the ambient plasma as well as the distribution of electric/magnetic fields are obtained from an MHD simulation [1] of the coma connected to the solar wind. Trajectories of ions and electrically charged dust grains are simulated by accounting for the Lorentz force and the nucleus gravity.

Viscous-Flow Approach to In Situ Infiltration and In Vitro Saturated Hydraulic Conductivity Determination

Infiltration is dominantly gravity driven, and a viscous-flow approach was developed. Laminar film flow equilibrates gravity with the viscous force and a constant flow velocity evolves during a period lasting 3/2 times the duration of a constant input rate, qS. Film thickness F and the specific contact area L of the film per unit soil volume are the key parameters. Sprinkler irrigation produced in situ time series of volumetric water contents, θ(z,t), as determined with TDR probes. The wetting front velocity v and the time series of the mobile water content, w(z,t) were deduced from θ(z,t). In vitro steady flow in a core of saturated soil provided volume flux density, q(z,t), and flow velocity, v, as determined from a heat front velocity. The F and L parameters of the in situ and the in vitro experiments were compared. The macropore-flow restriction states that, for a particular permeable medium, the specific contact area L must be independent from qS i.e., dL/dqS = 0. If true, then the relationship of qS ∝ v3/2 could scale a wide range of input rates 0 ≤ qS ≤ saturated hydraulic conductivity, Ksat, into a permeable medium, and kinematic-wave theory would become a versatile tool to deal with non-equilibrium flow. The viscous-flow approach is based on hydromechanical principles similar to Darcy’s law, but currently it is not suited to deduce flow properties from specified individual spatial structures of permeable media.

Gravitational slopes, geomorphology, and material strengths of the nucleus of comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Aims. We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material (tensile, shear, and compressive strengths). Methods. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface (Imhotep, Ash, Seth, Hathor, and Agilkia), using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies (overhangs, collapsed structures, boulders, cliffs, and Philae's footprint) and mechanical considerations. Results. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 degrees) are covered by a fine material and contain a few large (>10 m) and isolated boulders; ii) intermediate-slope terrains (20-45 degrees) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from <1m to 10m for the majority of them; and iii) high-slope terrains (45-90 degrees) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

Two independent and primitive envelopes of the bilobate nucleus of comet 67P

The factors shaping cometary nuclei are still largely unknown, but could be the result of concurrent effects of evolutionary(1,2) and primordial processes(3,4). The peculiar bilobed shape of comet 67P/Churyumov-Gerasimenko may be the result of the fusion of two objects that were once separate or the result of a localized excavation by outgassing at the interface between the two lobes(5). Here we report that the comet's major lobe is enveloped by a nearly continuous set of strata, up to 650 metres thick, which are independent of an analogous stratified envelope on the minor lobe. Gravity vectors computed for the two lobes separately are closer to perpendicular to the strata than those calculated for the entire nucleus and adjacent to the neck separating the two lobes. Therefore comet 67P/Churyumov-Gerasimenko is an accreted body of two distinct objects with 'onion-like' stratification, which formed before they merged. We conclude that gentle, low-velocity collisions occurred between two fully formed kilometre-sized cometesimals in the early stages of the Solar System. The notable structural similarities between the two lobes of comet 67P/Churyumov-Gerasimenko indicate that the early-forming cometesimals experienced similar primordial stratified accretion, even though they formed independently.