Gravitational features are a fundamental source of information to learn more about the interior structure and composition of planets, moons, asteroids and comets. Gravitational field modeling typically approximates the target body with a sphere, leading to a representation in spherical harmonics. However, small celestial bodies are often irregular in shape, and hence poorly approximated by a sphere. A much better suited geometrical fit is achieved by a tri-axial ellipsoid. This is also mirrored in the fact that the associated harmonic expansion (ellipsoidal harmonics) shows a significantly better convergence behavior as opposed to spherical harmonics. Unfortunately, complex mathematics and numerical problems (arithmetic overflow) so far severely limited the applicability of ellipsoidal harmonics. In this paper, we present a method that allows expanding ellipsoidal harmonics to a considerably higher degree compared to existing techniques. We apply this novel approach to model the gravitational field of comet 67P, the final target of the Rosetta mission. The comparison of results based on the ellipsoidal parameterization with those based on the spheroidal and spherical approximations reveals that the latter is clearly inferior; the spheroidal solution, on the other hand, is virtually just as accurate as the ellipsoidal one. Finally, in order to generalize our findings, we assess the gravitational field modeling performance for some 400 small bodies in the solar system. From this investigation we generally conclude that the spheroidal representation is an attractive alternative to the complex ellipsoidal parameterization on the one hand, and the inadequate spherical representation on the other hand.

The aim of this study is first to determine the gravity field of the comet 67P Churyumov-Gerasimenko and second to derive the solar component of the precession rate and nutation coefficients of the spin axis of the comet nucleus, i.e. without the direct, usually larger, effect of outgassing. The gravity field, and related moments of inertia, are obtained from two polyhedra, that are provided by the OSIRIS and NAV-CAM experiments on Rosetta, and are based on the assumption of uniform density for the comet nucleus. We also calculate the forced precession rate as well as the nutation coefficients on the basis of Kinoshita's theory of rotation of the rigid Earth and adapted it to be able to indirectly include the effect of outgassing on the rotational parameters. The 2nd degree denormalized Stokes coefficients of comet 67P C-G turn out to be (bracketed numbers refer to second shape model) C20 = -6.74 [-7.93] * 10^-2, C22 = 2.60 [2.71] * 10^-2 consistent with normalized principal moments of inertia A/MR^2 = 0.13 [0.11], B/MR^2 = 0.23 [0.22], with polar moment c = C/MR^2 = 0.25, depending on the choice of the polyhedron model. The obliquity between the rotation axis and the mean orbit normal is 52 degree, and the precession rate only due to solar torques becomes 20-30 "/y. Oscillations in longitude caused by the gravitational pull of the Sun turn out to be of the order of 1', oscillations in obliquity can be estimated to be of the order of 0.5'.

For the period July 2003 to August 2010, the interplanetary coronal mass ejection (ICME) catalogue maintained by Richardson and Cane lists 106 Earth-directed events, which have been measured in-situ by plasma and field instruments onboard the ACE satellite. We present a statistical investigation of the Earth's thermospheric neutral density response by means of accelerometer measurements collected by the GRACE satellites, which are available for 104 ICMEs in the data set, and its relation to various geomagnetic indices and characteristic ICME parameters such as the impact speed, southward magnetic field strength (Bz). The majority of ICMEs causes a distinct density enhancement in the thermosphere, with up to a factor of eight compared to the pre-event level. We find high correlations between ICME Bz and thermospheric density enhancements (~0.9), while the correlation with the ICME impact speed is somewhat smaller (~0.7). The geomagnetic indices revealing the highest correlations are Dst and SYM-H (~0.9), the lowest correlations are obtained for kp and AE (~0.7), which show a nonlinear relation with the thermospheric density enhancements. Separating the response for the shock sheath region and the magnetic structure of the ICME, we find that the Dst and SYM-H reveal a tighter relation to the Bz minimum in the magnetic structure of the ICME, whereas the polar cap indices show higher correlations with the Bz minimum in the shock sheath region. Since the strength of the Bz component - either in the sheath or the magnetic structure of the ICME - is highly correlated (~0.9) with the neutral density enhancement, we discuss the possibility of satellite orbital decay estimates based on magnetic field measurements at L1, i.e. before the ICME hits the Earth's magnetosphere. This will further stimulate progress in space weather understanding and applications regarding satellite operations.