  The 3:2 spin-orbit resonance between the rotational and orbital motions of
Mercury results from a functional dependance on a tidal friction adding to a
non-zero eccentricity with a permanent asymmetry in the equatorial plane of the
planet. The upcoming space missions, MESSENGER and BepiColombo with onboard
instrumentation capable of measuring the Mercury's rotational parameters,
stimulate the objective to attempt to an accurate theory of the planet's
rotation. We have used our BJV relativistic model of solar system integration
including the spin-orbit motion of the Moon. This model had been previously
built in accordance with the requirements of the Lunar Laser Ranging
observational accuracy. We extended this model to the spin-orbit couplings of
the terrestrial planets including Mercury; the updated model is called SONYR
(acronym of Spin-Orbit N-BodY Relativistic model). An accurate rotation of
Mercury has been then obtained. Moreover, the conception of the SONYR model is
suitable for analyzing the different families of hermean rotational librations.
We accurately identify the non-linear relations between the rotation of Mercury
and its dynamical figure (\cmr2, $C_{20}$, and $C_{22}$). Notably, for a
variation of 1% on the \cmr2 value, signatures in the $\phi$ hermean libration
in longitude as well as in the $\eta$ obliquity of the planet are respectively
0.45 arcseconds (as) and 2.4 milliarcseconds (mas). These determinations
provide new constraints on the internal structure of Mercury to be discussed
with the expected accuracy forecasted in the BepiColombo mission (respectively
3.2 and 3.7 as according to Milani et al 2001).
