Monitoring absolute vertical land motions and absolute sea-level changes from GPS and tide gauges data over French Polynesia

In this study,we estimate the absolute vertical land motions at three tidal stations with collocated Global Navigation Satellite System(GNSS)receivers over French Polynesia during the period 2007-2020,and obtain,as ancillary results,estimates of the absolute changes in sea level at the same locations.To verify our processing approach to determining vertical motion,we first modeled vertical motion at the International GNSS Service(IGS)THTI station located in the capital island of Tahiti and compared our estimate with previous independent determinations,with a good agreement.We obtained the following estimates for the vertical land motions at the tide gauges:Tubuai island,Austral Archipelago-0.92±0.17 mm/yr,Vairao village,Tahiti Iti:-0.49±0.39 mm/yr,Rikitea,Gambier Archipelago-0.43±0.17 mm/yr.The absolute variations of the sea level are:Tubuai island,Austral Archipelago 5.25±0.60 mm/yr,Vairao village,Tahiti Iti:3.62±0.52 mm/yr,Rikitea,Gambier Archipelago 1.52±0.23 mm/yr.We discuss these absolute values in light of the values obtained from altimetric measurements and other means in French Polynesia.

Combination of Tsoft and ET34-ANA-V80 software for the preprocessing and analysis of tide gauge data in French Polynesia

Since 2008 a network of five sea-level monitoring stations was progressively installed in French Polynesia.The stations are autonomous and data,collected at a sampling rate of 1 or 2 min,are not only recorded locally,but also transferred in real time by a radio-link to the NOAA through the GOES satellite.The new ET34-ANA-V80 version of ETERNA,initially developed for Earth Tides analysis,is now able to analyze ocean tides records.Through a two-step validation scheme,we took advantage of the flexibility of this new version,operated in conjunction with the preprocessing facilities of the Tsoft software,to recover co rrected data series able to model sea-level variations after elimination of the ocean tides signal.We performed the tidal analysis of the tide gauge data with the highest possible selectivity(optimal wave grouping)and a maximum of additional terms(shallow water constituents).Our goal was to provide corrected data series and modelled ocean tides signal to compute tide-free sea-level variations as well as tidal prediction models with centimeter precision.We also present in this study the characteristics of the ocean tides in French Polynesia and preliminary results concerning the non-tidal variations of the sea level concerning the tide gauge setting.

Monitoring Zenithal Total Delays over the three different climatic zones from IGS GPS final products:A comparison between the use of the VMF1 and GMF mapping functions

The International GNSS Service(IGS) final products(ephemeris and clocks-correction) have made the GNSS an indispensable low-cost tool for scientific research, for example sub-daily atmospheric water vapor monitoring. In this study, we investigate if there is a systematic difference coming from the choice between the Vienna Mapping Function 1(VMF1) and the Global Mapping Function(GMF) for the modeling of Zenith Total Delay(ZTD) estimates, as well as the Integrated Precipitable Water Vapor(IPWV) estimates that are deduced from them. As ZTD estimates cannot be fully separated from coordinate estimates, we also investigated the coordinate repeatability between subsequent measurements.For this purpose, we monitored twelve GNSS stations on a global scale, for each of the three climatic zones(polar, mid-latitudes and tropical), with four stations on each zone. We used an automated processing based on the Bernese GNSS Software Version 5.2 by applying the Precise Point Positioning(PPP)approach, L3 Ionosphere-free linear combination, 7 cutoff elevation angle and 2 h sampling. We noticed an excellent agreement with the ZTD estimates and coordinate repeatability for all the stations w.r.t to CODE(the Center for Orbit Determination in Europe) and USNO(US Naval Observatory) products, except for the Antarctic station(Davis) which shows systematic biases for the GMF related results. As a final step, we investigated the effect of using two mapping functions(VMF1 and GMF) to estimate the IPWV,w.r.t the IPWV estimates provided by the Integrated Global Radiosonde Archive(IGRA). The GPS-derived IPWV estimates are very close to the radiosonde-derived IPWV estimates, except for one station in the tropics(Tahiti).

A new method to determine the gravity field of small bodies from line-of-sight acceleration data

We present a new method to derive line-of-sight acceleration observables from spacecraft radio tracking data. The observables can be used to estimate the mass and gravity of a natural satellite as a spacecraft flyby. The corresponding observation model adapts to one-way and two/three-way tracking modes. As a test case for method validation and application, we estimated the mass and degree two gravity field for the Martian moon Phobos using simulated tracking data when the spacecraft Mars Express flew by Phobos on 2013 December 29. We have a few real tracking data during flyby and they will be used to confirm raw data simulation. The main purpose of this paper is to demonstrate the method of line-of-sight acceleration reduction from raw tracking data and the feasibility to estimate mass and gravity of a natural satellite using this type of observable. This novel method is potentially applicable to planet and asteroid gravity field studies combined with Doppler tracking data.

Strong Spatial Aggregation of Martian Surface Temperature Shaped by Spatial and Seasonal Variations in Meteorological and Environmental Factors

Spatio-temporal variation in the Martian surface temperature(MST)is an indicator of ground level thermal processes and hence a building block for climate models.However,the distribution of MST exhibits different levels of spatial aggregation or heterogeneity,and varies in space and time.Furthermore,the effect of regional differences in meteorological or environmental factors on the MST is not well understood.Thus,we investigated the degree of spatial autocorrelation of MST across the surface of Mars globally by Moran’s I,and identified the hot spots by GetisOrd G;*.We also estimated the regional differences in the influence of seasonally dominant factors including thermal inertia(TI),albedo,surface pressure,latitude,dust and slope on MST by a geographically weighted regression model.The results indicate(1)that MST is spatially aggregated and hot and cold spots varied over time and space.(2)Hemispheric differences in topography,surface TI and albedo were primarily responsible for the hemispheric asymmetry of hot spots.(3)The dominant factors varied by geographical locations and seasons.For example,the seasonal Hadley circulation dominates at the low-latitudes and CO;circulation at the high-latitudes.(4)Regions with extreme variations in topography and low TI were sensitive to meteorological and environmental factors such as dust and CO_(2)ice.We conclude that the spatial autocorrelation of MST and the spatial and seasonal heterogeneity of influencing factors must be considered when simulating Martian climate models.This work provides a reference for further exploration of Martian climatic processes.

Computation of the atmosphere-less light intensity curve during a total solar eclipse by using Lunar Reconnaissance Orbiter topography data and the DE430 astronomical ephemeris

Observations of the sky irradiation intensity in the visible wavelengths during a solar eclipse permit to model the Sun diameter,a key number to constrain the internal structure of our star.In this paper,we present an algorithm that takes advantage of the precise Moon topography from Lunar Reconnaissance Orbiter to compute,with a high resolution in time,the geometrical part(i.e.top-of-atmosphere,and for a given wavelength)of the sky irradiation at any given location on the Earth during these events.The algorithm is also able to model the Baily’s beads.We give as an application the theoretical computation of the light curve corresponding to the solar eclipse observed at Lakeland(Queensland,North Australia)on 2012 November 13.The application to real data,with the introduction of atmospheric and instrumental passbands,will be considered in a forthcoming paper.

Reducing loss of significance in the computation of Earth-based two-way Doppler observables for small body missions

Two-way Doppler measurement is a typical Earth-based radiometric technique for interplanetary spacecraft navigation and gravity science investigation.The most widely used model for the computation of two-way Doppler observables is Moyer’s differenced-range Doppler(DRD) formula,which is based on a Schwarzschild approximation of the Solar-System space-time.However,the computation of range difference in DRD formula is sensitive to round-off errors due to approximate numbers defined by the norm IEEE754 in all PCs.This paper presented two updated models and their corresponding detailed instructions for the computation of the two-way Doppler observables so as to impair the effects of this type of numerical error.These two models were validated by two case studies related to the Rosetta mission—asteroid Lutetia flyby and comet 67 P/Churyumov-Gerasimenko orbiting case.In these two cases,the numerical noise from the updated models can be reduced by two orders-of-magnitude in the computed two-way Doppler observables.The results showed an accuracy from better than 6 × 10^(-3) mm s^(-1) at 1 s counting time interval to better than 3 × 10^(-5) mm s^(-1) at 60 s counting time interval.

Simulated Gravity Field Estimation for the Main Belt Comet 133P/Elst-Pizarro Based on a Satellite-to-satellite Tracking Mode

Most asteroids and comets are formed in the early stages of the solar system and therefore contain a wealth of information about their birth.The asteroid exploration mission planned in the coming years by China will likely target the celestial body named 133P/Elst-Pizarro(estimated diameter of about 4 km).The orbit of this asteroid stays within the asteroid belt,but nevertheless,it displays a comet-like dust tail.In this study,we used differential tracking data between two simulated probes and the data from an Earth station to estimate 133P gravity field model.This observation mode is similar to how the gravity field was estimated for large celestial objects in the GRAIL and GRACE missions,but here the object is the very small 133P asteroid.We compared the estimated gravity fields obtained for 133P from the satellite-to-satellite combined with the Earth-based two-way range-rate observation mode,with only the Earth-based two-way range rate mode.The results show that the accuracy of the low-degree(4 degree and order)estimate of the gravity field is improved by one order of magnitude by using the satellite-to-satellite combined with the Earth-based two-way range-rate observation mode with respect to the Earthonly tracking.Furthermore,another order of magnitude improvement in the gravity field solution is gained by decreasing the orbit altitude from 12 to 8 km.

An Independent Degree-eight Mars Gravity Field Model and the Expected Results from the Tianwen-1 Mission

Tianwen-1 is China's first independent interplanetary exploration mission,targeting Mars,and includes orbiting,landing,and rover phases.Similar to previous Mars missions,the Tianwen-1 orbiter was designed for polar orbits during the scientific mission period but has an exceptional eccentricity of approximately 0.59.We provide the first independent eight-degree Martian gravity field model in this paper,which was developed exclusively by a team working in China with our independent software as well,based on about two months of radiometric Doppler and range data from only the Tianwen-1 mission.This model is independent from the models created by the groups at NASA Jet Propulsion Laboratory and Goddard Space Flight Center in the United States,as well as the Centre National d'Etudes Spatiales in France.Furthermore,in order to optimize the engineering and scientific benefits,we proposed a number of potential orbits for the extended Tianwen-1 mission.In order to solve a higher-degree independent Mars gravity field model,the viability of modifying the perigee height was investigated,with the priority considerations of fuel savings and implementation hazards being controlled.

A simulation experiment on the GM estimation for Comet 133P/Elst-Pizarro

In China’s asteroid mission to be launched around 2025,(7968)133 P/Elst-Pizarro(hereafter 133 P)will be the second target,after a visit to asteroid(469219)Kamo’oalewa.This paper describes a simulation of precise orbit determination for the spacecraft around comet 133 P,as well as estimation of its gravitational parameter(GM)value and the solar radiation pressure coefficient Cr for the spacecraft.Different cometocentric distances of 200,150 and 100 km orbits are considered,as well as two tracking modes:exclusive two-way range-rate mode(Earth station to spacecraft)and combinations of two-way range-rate and local spacecraft onboard ranging to the comet.Compared to exclusive two-way range-rate,the introduction of local ranging observables improves the final GM uncertainties by up to one order of magnitude.An ephemeris error in the orbit of 133 P is also considered,and we show that,to obtain a reliable estimate of the GM for 133 P,this error cannot exceed a one km range.

Reconstructing the cruise-phase trajectory of deep-space probes in a general relativistic framework:An application to the Cassini gravitational wave experiment

Einstein’s theory of general relativity is playing an increasingly important role in fields such as interplanetary navigation,astrometry,and metrology.Modern spacecraft and interplanetary probe prediction and estimation platforms employ a perturbed Newtonian framework,supplemented with the Einstein-Infeld-Hoffmann n-body equations of motion.While time in Newtonian mechanics is formally universal,the accuracy of modern radiometric tracking systems necessitate linear corrections via increasingly complex and error-prone post-Newtonian techniques—to account for light deflection due to the solar system bodies.With flagship projects such as the ESA/JAXA BepiColombo mission now operating at unprecedented levels of accuracy,we believe the standard corrected Newtonian paradigm is approaching its limits in terms of complexity.In this paper,we employ a novel prototype software,General Relativistic Accelerometer-based Propagation Environment,to reconstruct the Cassini cruise-phase trajectory during its first gravitational wave experiment in a fully relativistic framework.The results presented herein agree with post-processed trajectory information obtained from NASA’s SPICE kernels at the order of centimetres.

基于改进彗星星历的67P/Churyumov-Gerasimenko的三阶重力场模型

本文解算了Rosetta彗星探测任务的主要目标——彗星67P/Churyumov-Gerasimenko的重力场模型.与以往模型的主要区别是,本文使用了JPL团队最新解算的67P星历模型,该彗星星历在很大程度上得到了改进.67P的改进星历能够对其GM值的解算产生一定的影响,这一影响与67P通过近日点由于彗星喷发造成的质量损失处于同一水平.新的彗星重力场解算结果表明,67P重力异常比过去的模型更加显著.从2.65 km的参考球上看,新的重力场模型最大重力异常值约为0.9±0.03 mGal,较大的异常值表明了该彗星内部质量分布不均一的特性,这可能是由于高孔隙率和分层的彗核结构导致的.