<p>In Next Generation Gravity Missions (NGGM) the Laser Ranging Interferometer (LRI) is applied to measure inter-satellite range rate with nanometer-level precision. Thereby the precision of numerical orbit integration must be higher or at least same as that of LRI and the currently widely-used double-precision orbit integration technique cannot meet the numerical requirements of LRI measurements. Considering quadruple-precision orbit integration arithmetic is time consuming, we propose a hybrid-precision numerical orbit integration technique, in which the double- and quadruple-precision arithmetic is employed in the increment calculation part and orbit propagation part, respectively. Since the round-off errors are not sensitive to the time-demanding increment calculation but to the least time-consuming orbit propagation, the proposed hybrid-precision numerical orbit integration technique is as efficient as the double-precision orbit integration technique, and as precise as the quadruple-precision orbit integration. By using hybrid-precision orbit integration technique, the range rate precision is easily achieved at 10-12m/s in either nominal or Encke form, and furthermore the sub-nanometer-level range precision is obtainable in the Encke form with reference orbit selected as the best-fit one. Therefore, the hybrid-precision orbit integration technique is suggested to be used in the gravity field solutions for NGGM.</p>
