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 an...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.展开更多
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,...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.展开更多
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. ...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).展开更多
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. ...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.展开更多
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 level...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.展开更多
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...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.展开更多
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 o...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.展开更多
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 l...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.展开更多
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 pola...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.展开更多
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 orb...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.展开更多
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 estimat...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.展开更多
基金the University of French Polynesiafunding by several successive“Decision Aide a la Recherche”(DAR)grants to the Geodesy Observatory of Tahiti from the French Space Agency(CNES)+2 种基金fundings from the local government of French Polynesia(Observatoire Polynesien du Rechauffement Climatique)funding by“National Natural Science Foundation of China”(Grand No.41931075)funding by“the Fundamental Research Funds for the Central Universities"(Grand No.2042022kf1198)。
文摘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.
基金funding from the“Talent Introduction Scientific Research Start-Up Fund”of Shandong University of Science and Technology(Grant number 0104060510217)the“Open Fund of State Key Laboratory of Geodesy and Earth’s Dynamics”(Grant number SKLGED2021-3-5)。
文摘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.
基金the innovation carrier project by Zhejiang provincial science and Technology Department (2017F10008)the French Space Agency (CNES) for their funding, through a DAR grant to the Geodesy Observatory of Tahiti
文摘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).
基金supported by the National Natural Science Foundation of China (Nos. U1531136, U1831132 and U1531104)Innovation Group of Natural Fund of Hubei Province(2018CFA087)+1 种基金Open Funding of Macao University of Science and Technology (FDCT 119/2017/A3)Open Funding of Guizhou Provincial Key Laboratory of Radio Astronomy and Data Processing (KF201813)
文摘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.
基金supported by the pre-research Project on Civil Aerospace Technologies(No.D020103)supported by the National Natural Science Foundation of China(Grant No.42030110)。
文摘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.
基金the National Natural Science Foundation of China(U1831132 and 41804025)grant of Hubei Province Natural Science(2018CFA087)+2 种基金supported by the grant from Key Laboratory of Lunar and Deep Space Exploration,CASLIESMARS Special Research Fundingfunded through a DAR grant in planetology from the French Space Agency(CNES)。
文摘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.
基金supported by the National Natural Science Foundation of China (Grant Nos.42030110,41874010 and U1831132)Hubei Province Foundation innovation group project (2018CFA087)+1 种基金funded by a DAR grant in planetology from the French Space Agency (CNES)funded by Deutsches Zentrum für Luft-und Raumfahrt,Bonn-Oberkassel (Grant Nos.50QM1704 and 50QM1901)。
文摘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.
基金supported by the National Natural Science Foundation of China (42241116 and U1831132)supported by the Macao Science and Technology Development Fund (No. SKL-LPS(MUST)-20212023)+2 种基金the 2022 Project of Xinjiang Uygur Autonomous Region of China for Heaven Lake Talent Programsupported by the Chinese Academy of Sciences Foundation of the young scholars of western (2020-XBQNXZ-019)supported by a DAR grant in planetology from the French Space Agency (CNES).
文摘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.
基金supported by the National Natural Science Foundation of China(NSFC)under Nos.12203002 and 42241116Key Laboratory of Geospace Environment and Geodesy,Ministry of Education,Wuhan University under No.21-01-01funded by a DAR grant in planetology from the French Space Agency(CNES)。
文摘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.
基金the National Natural Science Foundation of China(U1831132 and 41874010)Innovation Group of Natural Fund of Hubei Province(2018CFA087)supported by a DAR fund in planetology from the French Space Agency(CNES)。
文摘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.
基金Open Access funding enabled and organized by CAUL and its Member InstitutionsJO’L acknowledges support from the Australian Research Council Centre of Excellence for Gravitational Wave Discovery(OzGrav)(Grant No.CE170100004)JPB was funded by a DAR grant in planetology from the French Space Agency(CNES).
文摘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.