Influence of subdaily polar motion model on nutation offsets estimated by very long baseline interferometry

This paper studies the connection between the subdaily model for polar motion used in the processing of very long baseline interferometry (VLBI) observations and the estimated nutation offsets. By convention accepted by the International Earth Rotation Service, the subdaily model for polar motion recommended for routine processing of geodetic observations does not contain any daily retrograde terms due to their one-to-one correlation with the nutation. Nevertheless, for a 24-h VLBI solution a part of the signal contained in the polar motion given by the used subdaily model is numerically mistaken for a retrograde daily sidereal signal. This fictitious retrograde daily signal contributes to the estimated nutation, leading to systematic differences between the nutation offsets from VLBI solutions computed with different subdaily polar motion models. We demonstrate this effect using solutions for all suitable 24-h VLBI sessions over a time span of 11 years (2000–2011). By changing the amplitudes of one tidal term in the underlying subdaily model for polar motion and comparing the estimated parameters to the solutions computed with the unchanged subdaily model, the paper shows and explains theoretically the effects produced by the individual subdaily terms on the VLBI nutation estimates.     

On the consistency of the current conventional EOP series and the celestial and terrestrial reference frames

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.     

Spherical integral transforms of second-order gravitational tensor components onto third-order gravitational tensor components

New spherical integral formulas between components of the second- and third-order gravitational tensors are formulated in this article. First, we review the nomenclature and basic properties of the second- and third-order gravitational tensors. Initial points of mathematical derivations, i.e., the second- and third-order differential operators defined in the spherical local North-oriented reference frame and the analytical solutions of the gradiometric boundary-value problem, are also summarized. Secondly, we apply the third-order differential operators to the analytical solutions of the gradiometric boundary-value problem which gives 30 new integral formulas transforming (1) vertical-vertical, (2) vertical-horizontal and (3) horizontal-horizontal second-order gravitational tensor components onto their third-order counterparts. Using spherical polar coordinates related sub-integral kernels can efficiently be decomposed into azimuthal and isotropic parts. Both spectral and closed forms of the isotropic kernels are provided and their limits are investigated. Thirdly, numerical experiments are performed to test the consistency of the new integral transforms and to investigate properties of the sub-integral kernels. The new mathematical apparatus is valid for any harmonic potential field and may be exploited, e.g., when gravitational/magnetic second- and third-order tensor components become available in the future. The new integral formulas also extend the well-known Meissl diagram and enrich the theoretical apparatus of geodesy.     

A spatial error model with continuous random effects and an application to growth convergence

We propose a spatial error model with continuous random effects based on Matérn covariance functions and apply this model for the analysis of income convergence processes (\(\beta \)-convergence). The use of a model with continuous random effects permits a clearer visualization and interpretation of the spatial dependency patterns, avoids the problems of defining neighborhoods in spatial econometrics models, and allows projecting the spatial effects for every possible location in the continuous space, circumventing the existing aggregations in discrete lattice representations. We apply this model approach to analyze the economic growth of Brazilian municipalities between 1991 and 2010 using unconditional and conditional formulations and a spatiotemporal model of convergence. The results indicate that the estimated spatial random effects are consistent with the existence of income convergence clubs for Brazilian municipalities in this period.     

Landscape transform and spatial metrics for mapping spatiotemporal land cover dynamics using Earth Observation data-sets

Analysis of Earth observation (EO) data, often combined with geographical information systems (GIS), allows monitoring of land cover dynamics over different ecosystems, including protected or conservation sites. The aim of this study is to use contemporary technologies such as EO and GIS in synergy with fragmentation analysis, to quantify the changes in the landscape of the Rajaji National Park (RNP) during the period of 19 years (1990–2009). Several statistics such as principal component analysis (PCA) and spatial metrics are used to understand the results. PCA analysis has produced two principal components (PC) and explained 84.1% of the total variance, first component (PC1) accounted for the 57.8% of the total variance while the second component (PC2) has accounted for the 26.3% of the total variance calculated from the core area metrics, distance metrics and shape metrics. Our results suggested that notable changes happened in the RNP landscape, evidencing the requirement of taking appropriate measures to conserve this natural ecosystem.     

Using low-density discrete Airborne Laser Scanning data to assess the potential carbon dioxide emission in case of a fire event in a Mediterranean pine forest

The aim of study is to map the carbon dioxide (CO₂) emission of the aboveground tree biomass (AGB) in case of a fire event. The suitability of low point density, discrete, multiple-return, Airborne Laser Scanning (ALS) data and the influence of several characteristics of these data and the study area on the results obtained have been evaluated. A sample of 45 circular plots representative of Pinus halepensis Miller stands were used to fit and validate the model of AGB. The ALS point clouds were processed to obtain the independent variables and a multivariate linear regression analysis between field data and ALS-derived variables allowed estimation of AGB. Then, the influence of several characteristics on the residuals of the model was analyzed. Finally, conversion factors were applied to obtain the CO₂ values. The AGB model presented a R² value of 0.84 with a relative root-mean-square error of 27.35%. This model included ALS variables related to vegetation height variability and to canopy density. Terrain slope, aspect, canopy cover, scan angle and the number of laser returns did not influence AGB estimations at plot level.     

SVLBI用于大地测量的卫星轨道及其跟踪网设计

SVLBI （space very long baseline interferometry） has some important potential applications in geodesy and geodynamics, for which one of the most difficult tasks is to precisely determine the orbit of an SVLBI satellite. This work studies several technologies that will possibly be able to determine the orbit of a space VLBI satellite. Then, according to the types and charac- teristics of the satellite and the requirements for geodetic study and the geometry of the GNSS （GPS, GALILEO） satellite to track the space VLBI satellite, the six Keplerian elements of the SVLBI satellite （TEST-SVLBI） are determined. A program is designed to analyze the coverage area of space of different altitudes by the stations of the network, with which the tracking network of TEST-SVLBI is designed. The efficiency of tracking TEST-SVLBI by the network is studied, and the results are presented.     

Spin rate estimation of sounding rockets using GPS wind-up

Carrier phase wind-up is a well-known effect that arises from the relative rotation between a transmitting and receiving antenna. In GPS measurements at L1 frequency, this effect translates into an error of 19.029 cm per full relative rotation of antennas. Since this effect is independent of the satellite elevation for pure rotation about the antenna boresight axis, it is usually absorbed by the clock estimation in navigation algorithms. Therefore, the impact of wind-up is usually neglected for applications that do not require accuracies to the cm level like RTK. However, in receiving platforms with high rotation rate, the accumulated wind-up value can be important and actually be larger than receiver noise or even ionospheric variations. Therefore, in such scenarios, the wind-up contribution can be isolated and used as a source of information to compute the spin rate of such platforms using an appropriate combination of GPS observables. This work shows some results of a coarse, yet simple, approach to monitor the rotation angle and spin-rate of spin stabilized sounding rockets flown by DLR.     

The ionosphere: effects,GPS modeling and the benefits for space geodetic techniques

The main goal of this paper is to provide a summary of our current knowledge of the ionosphere as it relates to space geodetic techniques, especially the most informative technology, global navigation satellite systems (GNSS), specifically the fully deployed and operational global positioning system (GPS). As such, the main relevant modeling points are discussed, and the corresponding results of ionospheric monitoring are related, which were mostly computed using GPS data and based on the direct experience of the authors. We address various phenomena such as horizontal and vertical ionospheric morphology in quiet conditions, traveling ionospheric disturbances, solar flares, ionospheric storms and scintillation. Finally, we also tackle the question of how improved knowledge of ionospheric conditions, especially in terms of an accurate understanding of the distribution of free electrons, can improve space geodetic techniques at different levels, such as higher-order ionospheric effects, precise GNSS navigation, single-antenna GNSS orientation and real-time GNSS meteorology.