Earth gravity field recovery using GPS, GLONASS, and SLR satellites

The time variable Earth’s gravity field provides the information about mass transport within the system Earth, i.e., the relationship of mass transport between atmosphere, oceans, and land hydrology. We recover the low-degree parameters of the time variable gravity field using microwave observations from GPS and GLONASS satellites and from SLR data to five geodetic satellites, namely LAGEOS-1/2, Starlette, Stella, and AJISAI. GPS satellites are particularly sensitive to specific coefficients of the Earth's gravity field, because of the deep 2:1 orbital resonance with Earth rotation (two revolutions of the GPS satellites per sidereal day). The resonant coefficients cause, among other, a “secular” drift (actually periodic variations of very long periods) of the semi-major axes of up to 5.3 m/day of GPS satellites. We processed 10 years of GPS and GLONASS data using the standard orbit models from the Center of Orbit Determination in Europe (CODE) with a simultaneous estimation of the Earth gravity field coefficients and other parameters, e.g., satellite orbit parameters, station coordinates, Earth rotation parameters, troposphere delays, etc. The weekly GNSS gravity solutions up to degree and order 4/4 are compared to the weekly SLR gravity field solutions. The SLR-derived geopotential coefficients are compared to monthly GRACE and CHAMP results.

Earth’s Rotation: A Challenging Problem in Mathematics and Physics

A suitable knowledge of the orientation and motion of the Earth in space is a common need in various fields. That knowledge has been ever necessary to carry out astronomical observations, but with the advent of the space age, it became essential for making observations of satellites and predicting and determining their orbits, and for observing the Earth from space as well. Given the relevant role it plays in Space Geodesy, Earth rotation is considered as one of the three pillars of Geodesy, the other two being geometry and gravity. Besides, research on Earth rotation has fostered advances in many fields, such as Mathematics, Astronomy and Geophysics, for centuries. One remarkable feature of the problem is in the extreme requirements of accuracy that must be fulfilled in the near future, about a millimetre on the tangent plane to the planet surface, roughly speaking. That challenges all of the theories that have been devised and used to-date; the paper makes a short review of some of the most relevant methods, which can be envisaged as milestones in Earth rotation research, emphasizing the Hamiltonian approach developed by the authors. Some contemporary problems are presented, as well as the main lines of future research prospected by the International Astronomical Union/International Association of Geodesy Joint Working Group on Theory of Earth Rotation, created in 2013.

Rotation de l'ellipsoide hétérogène et figure exacte de la terre.

Thèse - Faculté des sciences de Paris.

TeV scale implications of non commutative space time in laboratory frame with polarized beams

We analyze e(+)e(-) -> gamma gamma, e(-)gamma -> e(-)gamma and gamma gamma -> e(+)e(-) processes within the Seiberg-Witten expanded noncommutative scenario using polarized beams. With unpolarized beams the leading order effects of non commutativity starts from second order in non commutative(NC) parameter i.e. O(Theta(2)), while with polarized beams these corrections appear at first order (O(Theta')) in cross section. The corrections in Compton case can probe the magnetic component(Theta(B)) while in Pair production and Pair annihilation probe the electric component((Theta) over right arrow (E)) of NC parameter. We include the effects of earth rotation in our analysis. This study is done by investigating the effects of non commutativity on different time averaged cross section observables. The results which also depends on the position of the collider, can provide clear and distinct signatures of the model testable at the International Linear Collider(ILC).

Transformação de Helmert generalizada no posicionamento de alta precisão: fundamentação teórica e exemplificações

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A construção de conceitos relacionados com os movimentos Terra-Lua-Sol por alunos da E.J.A. à luz da teoria histórico-cultural

A construção de conceitos científicos no âmbito escolar ainda precisa ser melhor compreendida. No caso de conceitos relacionados à astronomia, as pesquisas ainda são escassas no Brasil, principalmente quando se trata da Educação de Jovens e Adultos (EJA). O presente estudo objetivou compreender, à luz da Teoria Histórico-Cultural, as formas através das quais estudantes da EJA constroem conceitos relacionados com os movimentos do sistema Terra-Lua-Sol, em suas interações com o professor e os colegas em sala de aula. A pesquisa foi realizada em uma turma do ensino médio da EJA da Escola de Aplicação da Universidade Federal do Pará. A turma era formada por 19 estudantes, com idades variando entre 16 e 37 anos. A coleta de informações foi feita durante um semestre letivo, inicialmente através de questionários. Eles continham perguntas abertas sobre a temática, a fim de investigar as concepções prévias dos estudantes. Posteriormente, as aulas em que o assunto foi ensinado foram gravadas em áudio e vídeo. Nestas aulas os alunos elaboraram individualmente e coletivamente explicações para a sucessão dia-noite na terra. Os grupos foram formados espontaneamente pelos alunos e, em seguida, foram recombinados pelo professor. Os registros foram transcritos e analisados microgeneticamente. As respostas dos estudantes ao questionamento inicial que tratava sobre a sucessão do dia e da noite na Terra foram categorizadas em quatro níveis A, B, C e D desde o mais afastado até o mais próximo do conceito escolar cientificamente aceito. Os resultados obtidos mostraram que 13 estudantes melhoram o perfil conceitual, pois migraram dos níveis A, B ou C para o nível D da categoria de respostas, ou seja, estes estudantes entendiam que a sucessão dianoite era decorrente do movimento de rotação da terra. Os outros seis estudantes, que já se encontravam no nível D, permaneceram nesse nível, porém melhoraram suas explicações em relação as suas respostas iniciais. Foram selecionados três estudantes para a análise microgenética dos percursos da construção de seus conhecimentos. Eles tinham suas respostas escritas iniciais classificadas nos níveis A ou B e durante suas interações com o professor e com os colegas incorporaram elementos do discurso científico, conseguindo elaborar explicações teóricas para o fenômeno observado. Os resultados ilustram diferentes mecanismos de ajuste da ajuda educacional oferecida pelo professor e pelos colegas, que salientam a importância de uma abordagem dialógica e do trabalho com diferentes formatos de grupos em sala de aula.

Impact of loading displacements on SLR-derived parameters and on the consistency between GNSS and SLR results

Displacements of the Earth’s surface caused by tidal and non-tidal loading forces are relevant in high-precision space geodesy. Some of the corrections are recommended by the international scientific community to be applied at the observation level, e.g., ocean tidal loading (OTL) and atmospheric tidal loading (ATL). Non-tidal displacement corrections are in general recommended not to be applied in the products of the International Earth Rotation and Reference Systems Service, in particular atmospheric non-tidal loading (ANTL), oceanic and hydrological non-tidal corrections. We assess and compare the impact of OTL, ATL and ANTL on SLR-derived parameters by reprocessing 12 years of SLR data considering and ignoring individual corrections. We show that loading displacements have an influence not only on station long-term stability, but also on geocenter coordinates, Earth Rotation Parameters, and satellite orbits. Applying the loading corrections reduces the amplitudes of annual signals in the time series of geocenter and station coordinates. The general improvement of the SLR station 3D coordinate repeatability when applying OTL, ATL and ANTL corrections are 19.5 %, 0.2 % and 3.3 % respectively, w.r.t. the solutions without loading corrections. ANTL corrections play a crucial role in the combination of optical (SLR) and microwave (GNSS, VLBI, DORIS) space geodetic observation techniques, because of the so-called Blue-Sky effect: SLR measurements can be carried out only under cloudless sky conditions—typically during high air pressure conditions, when the Earth’s crust is deformed, whereas microwave observations are weather-independent. Thus, applying the loading corrections at the observation level improves SLR-derived products as well as the consistency with microwave-based results. We assess the Blue-Sky effect on SLR stations and the consistency improvement between GNSS and SLR solutions when ANTL corrections are included. The omission of ANTL corrections may lead to inconsistencies between SLR and GNSS solutions of up to 2.5 mm for inland stations. As a result, the estimated GNSS–SLR coordinate differences correspond better to the local ties at the co-located stations when applying ANTL corrections.

SLR-derived terrestrial reference frame using the observations to LAGEOS-1/2, Starlette, Stella, and AJISAI

Currently, the contributions of Starlette, Stella, and AJISAI are not taken into account when defining the International Terrestrial Reference Frame (ITRF), despite the large amount of data collected in a long time-span. Consequently, the SLR-derived parameters and the SLR part of the ITRF are almost exclusively defined by LAGEOS-1 and LAGEOS-2. We investigate the potential of combining the observations to several SLR satellites with different orbital characteristics. Ten years of SLR data are homogeneously processed using the development version 5.3 of the Bernese GNSS Software. Special emphasis is put on orbit parameterization and the impact of LEO data on the estimation of the geocenter coordinates, Earth rotation parameters, Earth gravity field coefficients, and the station coordinates in one common adjustment procedure. We find that the parameters derived from the multi-satellite solutions are of better quality than those obtained in single satellite solutions or solutions based on the two LAGEOS satellites. A spectral analysis of the SLR network scale w.r.t. SLRF2008 shows that artifacts related to orbit perturbations in the LAGEOS-1/2 solutions, i.e., periods related to the draconitic years of the LAGEOS satellites, are greatly reduced in the combined solutions.

Homogeneous reprocessing of GPS, GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from −35 and −38 mm to −12 and −13 mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.

CODE’s new ultra-rapid orbit and ERP products for the IGS

The International GNSS Service (IGS) issues four sets of so-called ultra-rapid products per day, which are based on the contributions of the IGS Analysis Centers. The traditional (“old”) ultra-rapid orbit and earth rotation parameters (ERP) solution of the Center for Orbit Determination in Europe (CODE) was based on the output of three consecutive 3-day long-arc rapid solutions. Information from the IERS Bulletin A was required to generate the predicted part of the old CODE ultra-rapid product. The current (“new”) product, activated in November 2013, is based on the output of exactly one multi-day solution. A priori information from the IERS Bulletin A is no longer required for generating and predicting the orbits and ERPs. This article discusses the transition from the old to the new CODE ultra-rapid orbit and ERP products and the associated improvement in reliability and performance. All solutions used in this article were generated with the development version of the Bernese GNSS Software. The package was slightly extended to meet the needs of the new CODE ultra-rapid generation.

CODE’s new solar radiation pressure model for GNSS orbit determination

The Empirical CODE Orbit Model (ECOM) of the Center for Orbit Determination in Europe (CODE), which was developed in the early 1990s, is widely used in the International GNSS Service (IGS) community. For a rather long time, spurious spectral lines are known to exist in geophysical parameters, in particular in the Earth Rotation Parameters (ERPs) and in the estimated geocenter coordinates, which could recently be attributed to the ECOM. These effects grew creepingly with the increasing influence of the GLONASS system in recent years in the CODE analysis, which is based on a rigorous combination of GPS and GLONASS since May 2003. In a first step we show that the problems associated with the ECOM are to the largest extent caused by the GLONASS, which was reaching full deployment by the end of 2011. GPS-only, GLONASS-only, and combined GPS/GLONASS solutions using the observations in the years 2009–2011 of a global network of 92 combined GPS/GLONASS receivers were analyzed for this purpose. In a second step we review direct solar radiation pressure (SRP) models for GNSS satellites. We demonstrate that only even-order short-period harmonic perturbations acting along the direction Sun-satellite occur for GPS and GLONASS satellites, and only odd-order perturbations acting along the direction perpendicular to both, the vector Sun-satellite and the spacecraft’s solar panel axis. Based on this insight we assess in the third step the performance of four candidate orbit models for the future ECOM. The geocenter coordinates, the ERP differences w. r. t. the IERS 08 C04 series of ERPs, the misclosures for the midnight epochs of the daily orbital arcs, and scale parameters of Helmert transformations for station coordinates serve as quality criteria. The old and updated ECOM are validated in addition with satellite laser ranging (SLR) observations and by comparing the orbits to those of the IGS and other analysis centers. Based on all tests, we present a new extended ECOM which substantially reduces the spurious signals in the geocenter coordinate z (by about a factor of 2–6), reduces the orbit misclosures at the day boundaries by about 10 %, slightly improves the consistency of the estimated ERPs with those of the IERS 08 C04 Earth rotation series, and substantially reduces the systematics in the SLR validation of the GNSS orbits.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Testing a new Free Core Nutation empirical model

The Free Core Nutation (FCN) is a free mode of the Earth's rotation caused by the different material characteristics of the Earth's core and mantle. This causes the rotational axes of those layers to slightly diverge from each other, resulting in a wobble of the Earth's rotation axis comparable to nutations. In this paper we focus on estimating empirical FCN models using the observed nutations derived from the VLBI sessions between 1993 and 2013. Assuming a fixed value for the oscillation period, the time-variable amplitudes and phases are estimated by means of multiple sliding window analyses. The effects of using different a priori Earth Rotation Parameters (ERP) in the derivation of models are also addressed. The optimal choice of the fundamental parameters of the model, namely the window width and step-size of its shift, is searched by performing a thorough experimental analysis using real data. The former analyses lead to the derivation of a model with a temporal resolution higher than the one used in the models currently available, with a sliding window reduced to 400 days and a day-by-day shift. It is shown that this new model increases the accuracy of the modeling of the observed Earth's rotation. Besides, empirical models determined from USNO Finals as a priori ERP present a slightly lower Weighted Root Mean Square (WRMS) of residuals than IERS 08 C04 along the whole period of VLBI observations, according to our computations. The model is also validated through comparisons with other recognized models. The level of agreement among them is satisfactory. Let us remark that our estimates give rise to the lowest residuals and seem to reproduce the FCN signal in more detail.

Satellite laser ranging and the ETALON Geodetic satellites

The technique of Satellite Laser Ranging is today a mature, important tool with applications in many area of geodynamics, geodesy and satellite dynamics. A global network of some 40 stations regularly obtains range observations with sub-cm precision to more than twelve orbiting spacecraft. At such levels of precision it is important to minimise potential sources of range bias in the observations, and part of the thesis is a study of subtle effects caused by the extended nature of the arrays of retro-reflectors on the satellites. We develop models that give a precise correction of the range measurements to the centres of mass of the geodetic satellites Lageos and Etalon, appropriate to a variety of different ranging systems, and use the Etalon values, which were not determined during pre-launch tests, in an extended orbital analysis. We have fitted continuous 2.5 year orbits to range observations of the Etalons from the global network of stations, and analysed the results by mapping the range residuals from these orbits into equivalent corrections to orbital elements over short time intervals. From these residuals we have detected and studied large un-modelled along-track accelerations associated with periods during which the satellites are undergoing eclipse by the Earth's shadow. We also find that the eccentricity residuals are significantly different for the two satellites, with Etalon-2 undergoing a year-long eccentricity anomaly similar in character to that experienced at intervals by Lageos-1. The nodal residuals show that the satellites define a very stable reference frame for Earth rotation determination, with very little drift-off during the 2.5 year period. We show that an analysis of more than about eight years of tracking data would be required to derive a significant value for 2. The reference frame defined by the station coordinates derived from the analyses shows very good agreement with that of ITRF93.

Shape and Rotation Modeling and Thermophysical Analysis of Near-Earth Asteroid (1917) Cuyo

We are conducting an ESO Large Program that includes optical photometry, thermal-IR observations, and optical-NIR spectroscopy of selected NEAs. Among the principal goals of the program are shape and spin-state modeling, and searching for YORP-induced changes in rotation periods. One of our targets is asteroid (1917) Cuyo, a near-Earth asteroid from the Amor group. We carried out an extensive observing campaign on Cuyo between April 2010 and April 2013, operating primarily at the ESO 3.6m NTT for optical photometry, and the 8.2m VLT at Paranal for thermal-IR imaging. Further optical observations were acquired at the ESO 2.2m telescope, the Palomar 200" Hale telescope (California), JPL’s Table Mountain Observatory (California) and the Faulkes Telescope South (Australia). We obtained optical imaging data for rotational lightcurves throughout this period, as the asteroid passed through a wide range of observational geometries, conducive to producing a good shape model and spin state solution. The preliminary shape and spin state model indicates a nearly spherical shape and a rotation pole at ecliptic longitude λ = 53° ± 20° and latitude β = -37° ± 10° (1-sigma error bars are approximate). The sidereal rotation period was measured to be 2.6899522 ± (3 × 10^-7) hours. Linkage with earlier lightcurve data shows possible evidence of a small change in rotation rate during the period 1989-2013. We applied the NEATM thermal model (Harris A., Icarus 131, 291, 1998) to our VLT thermal-IR measurements (8-19.6 μm), obtained in September and December 2011. The derived effective diameter ranges from 3.4 to 4.2 km, and the geometric albedo is 0.16 (+0.07, -0.04). Using the shape model and thermal fluxes we will perform a detailed thermophysical analysis using the new Advanced Thermophysical Model (Rozitis, B. & Green, S.F., MNRAS 415, 2042, 2011; Rozitis, B. & Green, S.F., MNRAS 423, 367, 2012). This work was performed in part at the Jet Propulsion Laboratory under a contract with NASA.