Application of ray-traced tropospheric slant delays to geodetic VLBI analysis

The correction of tropospheric influences via so-called path delays is critical for the analysis of observations from space geodetic techniques like the very long baseline interferometry (VLBI). In standard VLBI analysis, the a priori slant path delays are determined using the concept of zenith delays, mapping functions and gradients. The a priori use of ray-traced delays, i.e., tropospheric slant path delays determined with the technique of ray-tracing through the meteorological data of numerical weather models (NWM), serves as an alternative way of correcting the influences of the troposphere on the VLBI observations within the analysis. In the presented research, the application of ray-traced delays to the VLBI analysis of sessions in a time span of 16.5 years is investigated. Ray-traced delays have been determined with program RADIATE (see Hofmeister in Ph.D. thesis, Department of Geodesy and Geophysics, Faculty of Mathematics and Geoinformation, Technische Universität Wien. http://resolver.obvsg.at/urn:nbn:at:at-ubtuw:1-3444, 2016) utilizing meteorological data provided by NWM of the European Centre for Medium-Range Weather Forecasts (ECMWF). In comparison with a standard VLBI analysis, which includes the tropospheric gradient estimation, the application of the ray-traced delays to an analysis, which uses the same parameterization except for the a priori slant path delay handling and the used wet mapping factors for the zenith wet delay (ZWD) estimation, improves the baseline length repeatability (BLR) at 55.9% of the baselines at sub-mm level. If no tropospheric gradients are estimated within the compared analyses, 90.6% of all baselines benefit from the application of the ray-traced delays, which leads to an average improvement of the BLR of 1 mm. The effects of the ray-traced delays on the terrestrial reference frame are also investigated. A separate assessment of the RADIATE ray-traced delays is carried out by comparison to the ray-traced delays from the National Aeronautics and Space Administration Goddard Space Flight Center (NASA GSFC) (Eriksson and MacMillan in http://lacerta.gsfc.nasa.gov/tropodelays, 2016) with respect to the analysis performances in terms of BLR results. If tropospheric gradient estimation is included in the analysis, 51.3% of the baselines benefit from the RADIATE ray-traced delays at sub-mm difference level. If no tropospheric gradients are estimated within the analysis, the RADIATE ray-traced delays deliver a better BLR at 63% of the baselines compared to the NASA GSFC ray-traced delays.     

Co-location of space geodetic techniques carried out at the Geodetic Observatory Wettzell using a closure in time and a multi-technique reference target

The quality of the links between the different space geodetic techniques (VLBI, SLR, GNSS and DORIS) is still one of the major limiting factors for the realization of a unique global terrestrial reference frame that is accurate enough to allow the monitoring of the Earth system, i.e., of processes like sea level change, postglacial rebound and silent earthquakes. According to the specifications of the global geodetic observing system of the International Association of Geodesy, such a reference frame should be accurate to 1 mm over decades, with rates of change stable at the level of 0.1 mm/year. The deficiencies arise from inaccurate or incomplete local ties at many fundamental sites as well as from systematic instrumental biases in the individual space geodetic techniques. Frequently repeated surveys, the continuous monitoring of antenna heights and the geometrical mount stability (Lösler et al. in J Geod 90:467–486, 2016.  https://doi.org/10.1007/s00190-016-0887-8) have not provided evidence for insufficient antenna stability. Therefore, we have investigated variations in the respective system delays caused by electronic circuits, which is not adequately captured by the calibration process, either because of subtle differences in the circuitry between geodetic measurement and calibration, high temporal variability or because of lacking resolving bandwidth. The measured system delay variations in the electric chain of both VLBI- and SLR systems reach the order of 100 ps, which is equivalent to 3 cm of path length. Most of this variability is usually removed by the calibrations but by far not all. This paper focuses on the development of new technologies and procedures for co-located geodetic instrumentation in order to identify and remove systematic measurement biases within and between the individual measurement techniques. A closed-loop optical time and frequency distribution system and a common inter-technique reference target provide the possibility to remove variable system delays. The main motivation for the newly established central reference target, locked to the station clock, is the combination of all space geodetic instruments at a single reference point at the observatory. On top of that it provides the unique capability to perform a closure measurement based on the observation of time.     

Determination of the optimized single-layer ionospheric height for electron content measurements over China

The ionosphere effective height (IEH) is a very important parameter in total electron content (TEC) measurements under the widely used single-layer model assumption. To overcome the requirement of a large amount of simultaneous vertical and slant ionospheric observations or dense “coinciding” pierce points data, a new approach comparing the converted vertical TEC (VTEC) value using mapping function based on a given IEH with the “ground truth” VTEC value provided by the combined International GNSS Service Global Ionospheric Maps is proposed for the determination of the optimal IEH. The optimal IEH in the Chinese region is determined using three different methods based on GNSS data. Based on the ionosonde data from three different locations in China, the altitude variation of the peak electron density (hmF2) is found to have clear diurnal, seasonal and latitudinal dependences, and the diurnal variation of hmF2 varies from approximately 210 to 520 km in Hainan. The determination of the optimal IEH employing the inverse method suggested by Birch et al. (Radio Sci 37, 2002. doi: 10.1029/2000rs002601) did not yield a consistent altitude in the Chinese region. Tests of the method minimizing the mapping function errors suggested by Nava et al. (Adv Space Res 39:1292–1297, 2007) indicate that the optimal IEH ranges from 400 to 600 km, and the height of 450 km is the most frequent IEH at both high and low solar activities. It is also confirmed that the IEH of 450–550 km is preferred for the Chinese region instead of the commonly adopted 350–450 km using the determination method of the optimal IEH proposed in this paper.     

Improved efficiency of maximum likelihood analysis of time series with temporally correlated errors

Most time series of geophysical phenomena have temporally correlated errors. From these measurements, various parameters are estimated. For instance, from geodetic measurements of positions, the rates and changes in rates are often estimated and are used to model tectonic processes. Along with the estimates of the size of the parameters, the error in these parameters needs to be assessed. If temporal correlations are not taken into account, or each observation is assumed to be independent, it is likely that any estimate of the error of these parameters will be too low and the estimated value of the parameter will be biased. Inclusion of better estimates of uncertainties is limited by several factors, including selection of the correct model for the background noise and the computational requirements to estimate the parameters of the selected noise model for cases where there are numerous observations. Here, I address the second problem of computational efficiency using maximum likelihood estimates (MLE). Most geophysical time series have background noise processes that can be represented as a combination of white and power-law noise, \(1/f^{\alpha }\) with frequency, f. With missing data, standard spectral techniques involving FFTs are not appropriate. Instead, time domain techniques involving construction and inversion of large data covariance matrices are employed. Bos et al. (J Geod, 2013. doi: 10.1007/s00190-012-0605-0) demonstrate one technique that substantially increases the efficiency of the MLE methods, yet is only an approximate solution for power-law indices >1.0 since they require the data covariance matrix to be Toeplitz. That restriction can be removed by simply forming a data filter that adds noise processes rather than combining them in quadrature. Consequently, the inversion of the data covariance matrix is simplified yet provides robust results for a wider range of power-law indices.     

On the equivalence of spherical splines with least-squares collocation and Stokes’s formula for regional geoid computation

This work is an investigation of three methods for regional geoid computation: Stokes’s formula, least-squares collocation (LSC), and spherical radial base functions (RBFs) using the spline kernel (SK). It is a first attempt to compare the three methods theoretically and numerically in a unified framework. While Stokes integration and LSC may be regarded as classic methods for regional geoid computation, RBFs may still be regarded as a modern approach. All methods are theoretically equal when applied globally, and we therefore expect them to give comparable results in regional applications. However, it has been shown by de Min (Bull Géod 69:223–232, 1995. doi: 10.1007/BF00806734) that the equivalence of Stokes’s formula and LSC does not hold in regional applications without modifying the cross-covariance function. In order to make all methods comparable in regional applications, the corresponding modification has been introduced also in the SK. Ultimately, we present numerical examples comparing Stokes’s formula, LSC, and SKs in a closed-loop environment using synthetic noise-free data, to verify their equivalence. All agree on the millimeter level.     

Exploring the spatially varying innovation capacity of the US counties in the framework of Griliches’ knowledge production function: a mixed GWR approach

Griliches’ knowledge production function has been increasingly adopted at the regional level where location-specific conditions drive the spatial differences in knowledge creation dynamics. However, the large majority of such studies rely on a traditional regression approach that assumes spatially homogenous marginal effects of knowledge input factors. This paper extends the authors’ previous work (Kang and Dall’erba in Int Reg Sci Rev, 2015. doi: 10.1177/0160017615572888) to investigate the spatial heterogeneity in the marginal effects by using nonparametric local modeling approaches such as geographically weighted regression (GWR) and mixed GWR with two distinct samples of the US Metropolitan Statistical Area (MSA) and non-MSA counties. The results indicate a high degree of spatial heterogeneity in the marginal effects of the knowledge input variables, more specifically for the local and distant spillovers of private knowledge measured across MSA counties. On the other hand, local academic knowledge spillovers are found to display spatially homogenous elasticities in both MSA and non-MSA counties. Our results highlight the strengths and weaknesses of each county’s innovation capacity and suggest policy implications for regional innovation strategies.     

VMF3/GPT3: refined discrete and empirical troposphere mapping functions

Incorrect modeling of troposphere delays is one of the major error sources for space geodetic techniques such as Global Navigation Satellite Systems (GNSS) or Very Long Baseline Interferometry (VLBI). Over the years, many approaches have been devised which aim at mapping the delay of radio waves from zenith direction down to the observed elevation angle, so-called mapping functions. This paper contains a new approach intended to refine the currently most important discrete mapping function, the Vienna Mapping Functions 1 (VMF1), which is successively referred to as Vienna Mapping Functions 3 (VMF3). It is designed in such a way as to eliminate shortcomings in the empirical coefficients b and c and in the tuning for the specific elevation angle of \(3^{\circ }\). Ray-traced delays of the ray-tracer RADIATE serve as the basis for the calculation of new mapping function coefficients. Comparisons of modeled slant delays demonstrate the ability of VMF3 to approximate the underlying ray-traced delays more accurately than VMF1 does, in particular at low elevation angles. In other words, when requiring highest precision, VMF3 is to be preferable to VMF1. Aside from revising the discrete form of mapping functions, we also present a new empirical model named Global Pressure and Temperature 3 (GPT3) on a \(5^{\circ }\times 5^{\circ }\) as well as a \(1^{\circ }\times 1^{\circ }\) global grid, which is generally based on the same data. Its main components are hydrostatic and wet empirical mapping function coefficients derived from special averaging techniques of the respective (discrete) VMF3 data. In addition, GPT3 also contains a set of meteorological quantities which are adopted as they stand from their predecessor, Global Pressure and Temperature 2 wet. Thus, GPT3 represents a very comprehensive troposphere model which can be used for a series of geodetic as well as meteorological and climatological purposes and is fully consistent with VMF3.     

Detecting communities with the multi-scale Louvain method: robustness test on the metropolitan area of Brussels

Detecting communities in large networks has become a common practice in socio-spatial analyses and has led to the development of numerous dedicated mathematical algorithms. Nowadays, however, researchers face a deluge of data and algorithms, and great care must be taken regarding methodological questions such as the values of the parameters and the geographical characteristics of the data. We aim here at testing the sensitivity of multi-scale modularity optimized by the Louvain method to the value of the resolution parameter (introduced by Reichardt and Bornholdt (Phys Rev Lett 93(21):218701, 2004.  https://doi.org/10.1103/PhysRevLett.93.218701) and controlling the size of the communities) and to a number of spatial issues such as the inclusion of internal loops and the delineation of the study area. We compare the community structures with those found by another well-known community detection algorithm (Infomap), and we further interpret the final results in terms of urban geography. Sensitivity analyses are conducted for commuting movements in and around Brussels. Results reveal slight effects of spatial issues (inclusion of the internal loops, definition of the study area) on the partition into job basins, while the resolution parameter plays a major role in the final results and their interpretation in terms of urban geography. Community detection methods seem to reveal a surprisingly strong spatial effect of commuting patterns: Similar partitions are obtained with different methods. This paper highlights the advantages and sensitivities of the multi-scale Louvain method and more particularly of defining communities of places. Despite these sensitivities, the method proves to be a valuable tool for geographers and planners.     

Short note: the experimental geopotential model XGM2016

As a precursor study for the upcoming combined Earth Gravitational Model 2020 (EGM2020), the Experimental Gravity Field Model XGM2016, parameterized as a spherical harmonic series up to degree and order 719, is computed. XGM2016 shares the same combination methodology as its predecessor model GOCO05c (Fecher et al. in Surv Geophys 38(3): 571–590, 2017. doi: 10.1007/s10712-016-9406-y). The main difference between these models is that XGM2016 is supported by an improved terrestrial data set of \(15^\prime \times 15^\prime \) gravity anomaly area-means provided by the United States National Geospatial-Intelligence Agency (NGA), resulting in significant upgrades compared to existing combined gravity field models, especially in continental areas such as South America, Africa, parts of Asia, and Antarctica. A combination strategy of relative regional weighting provides for improved performance in near-coastal ocean regions, including regions where the altimetric data are mostly unchanged from previous models. Comparing cumulative height anomalies, from both EGM2008 and XGM2016 at degree/order 719, yields differences of 26 cm in Africa and 40 cm in South America. These differences result from including additional information of satellite data, as well as from the improved ground data in these regions. XGM2016 also yields a smoother Mean Dynamic Topography with significantly reduced artifacts, which indicates an improved modeling of the ocean areas.     

Automatic Building Extraction from Image and LiDAR Data with Active Contour Segmentation

Automatic building extraction is an important topic for many applications such as urban planning, disaster management, 3D building modeling and updating GIS databases. Its approaches mainly depend on two data sources: light detection and ranging (LiDAR) point cloud and aerial imagery both of which have advantages and disadvantages of their own. In this study, in order to benefit from the advantages of each data sources, LiDAR and image data combined together. And then, the building boundaries were extracted with the automated active contour algorithm implemented in MATLAB. Active contour algorithm uses initial contour positions to segment an object in the image. Initial contour positions were detected without user interaction by a series of image enhancements, band ratio and morphological operations. Four test areas with varying building and background levels of detail were selected from ISPRS’s benchmark Vaihingen and Istanbul datasets. Vegetation and shadows were removed from all the datasets by band ratio to improve segmentation quality. Subsequently, LiDAR point cloud data was converted to raster format and added to the aerial imagery as an extra band. Resulting merged image and initial contour positions were given to the active contour algorithm to extract building boundaries. In order to compare the contribution of LiDAR to the proposed method, the boundaries of the buildings were extracted from the input image before and after adding LiDAR data to the image as a layer. Finally extracted building boundaries were smoothed by the Awrangjeb (Int J Remote Sen 37(3): 551–579.  https://doi.org/10.1080/01431161.2015.1131868, 2016) boundary regularization algorithm. Correctness (Corr), completeness (Comp) and accuracy (Q) metrics were used to assess accuracy of segmented building boundaries by comparing extracted building boundaries with manually digitized building boundaries. Proposed approach shows the promising results with over 93% correctness, 92% completeness and 89% quality.     

Global empirical model for mapping zenith wet delays onto precipitable water

We can map zenith wet delays onto precipitable water with a conversion factor, but in order to calculate the exact conversion factor, we must precisely calculate its key variable $T_\mathrm{m}$ . Yao et al. (J Geod 86:1125–1135, 2012. doi:10.1007/s00190-012-0568-1) established the first generation of global $T_\mathrm{m}$ model (GTm-I) with ground-based radiosonde data, but due to the lack of radiosonde data at sea, the model appears to be abnormal in some areas. Given that sea surface temperature varies less than that on land, and the GPT model and the Bevis $T_\mathrm{m}$ – $T_\mathrm{s}$ relationship are accurate enough to describe the surface temperature and $T_\mathrm{m}$ , this paper capitalizes on the GPT model and the Bevis $T_\mathrm{m}$ – $T_\mathrm{s}$ relationship to provide simulated $T_\mathrm{m}$ at sea, as a compensation for the lack of data. Combined with the $T_\mathrm{m}$ from radiosonde data, we recalculated the GTm model coefficients. The results show that this method not only improves the accuracy of the GTm model significantly at sea but also improves that on land, making the GTm model more stable and practically applicable.     

A complete VLBI delay model for deforming radio telescopes: the Effelsberg case

Deformations of radio telescopes used in geodetic and astrometric very long baseline interferometry (VLBI) observations belong to the class of systematic error sources which require correction in data analysis. In this paper we present a model for all path length variations in the geometrical optics of radio telescopes which are due to gravitational deformation. The Effelsberg 100 m radio telescope of the Max Planck Institute for Radio Astronomy, Bonn, Germany, has been surveyed by various terrestrial methods. Thus, all necessary information that is needed to model the path length variations is available. Additionally, a ray tracing program has been developed which uses as input the parameters of the measured deformations to produce an independent check of the theoretical model. In this program as well as in the theoretical model, the illumination function plays an important role because it serves as the weighting function for the individual path lengths depending on the distance from the optical axis. For the Effelsberg telescope, the biggest contribution to the total path length variations is the bending of the main beam located along the elevation axis which partly carries the weight of the paraboloid at its vertex. The difference in total path length is almost \(-\) 100 mm when comparing observations at 90 \(^\circ \) and at 0 \(^\circ \) elevation angle. The impact of the path length corrections is validated in a global VLBI analysis. The application of the correction model leads to a change in the vertical position of \(+120\)  mm. This is more than the maximum path length, but the effect can be explained by the shape of the correction function.     

The AUSTRAL VLBI observing program

The AUSTRAL observing program was started in 2011, performing geodetic and astrometric very long baseline interferometry (VLBI) sessions using the new Australian AuScope VLBI antennas at Hobart, Katherine, and Yarragadee, with contribution from the Warkworth (New Zealand) 12 m and Hartebeesthoek (South Africa) 15 m antennas to make a southern hemisphere array of telescopes with similar design and capability. Designed in the style of the next-generation VLBI system, these small and fast antennas allow for a new way of observing, comprising higher data rates and more observations than the standard observing sessions coordinated by the International VLBI Service for Geodesy and Astrometry (IVS). In this contribution, the continuous development of the AUSTRAL sessions is described, leading to an improvement of the results in terms of baseline length repeatabilities by a factor of two since the start of this program. The focus is on the scheduling strategy and increased number of observations, aspects of automated operation, and data logistics, as well as results of the 151 AUSTRAL sessions performed so far. The high number of the AUSTRAL sessions makes them an important contributor to VLBI end-products, such as the terrestrial and celestial reference frames and Earth orientation parameters. We compare AUSTRAL results with other IVS sessions and discuss their suitability for the determination of baselines, station coordinates, source coordinates, and Earth orientation parameters.     

Adding source positions to the IVS combination—First results

The consistent estimation of terrestrial reference frames (TRF), celestial reference frames (CRF) and Earth orientation parameters (EOP) is still an open subject and offers a large field of investigations. Until now, source positions resulting from Very Long Baseline Interferometry (VLBI) observations are not routinely combined on the level of normal equations in the same way as it is a common process for station coordinates and EOPs. The combination of source positions based on VLBI observations is now integrated in the IVS combination process. We present the studies carried out to evaluate the benefit of the combination compared to individual solutions. On the level of source time series, improved statistics regarding weighted root mean square have been found for the combination in comparison with the individual contributions. In total, 67 stations and 907 sources (including 291 ICRF2 defining sources) are included in the consistently generated CRF and TRF covering 30 years of VLBI contributions. The rotation angles \(A_1\), \(A_2\) and \(A_3\) relative to ICRF2 are ?12.7, 51.7 and 1.8 \({\upmu }\) as, the drifts \(D_\alpha \) and \(D_\delta \) are ?67.2 and 19.1 \(\upmu \) as/rad and the bias \(B_\delta \) is 26.1 \(\upmu \) as. The comparison of the TRF solution with the IVS routinely combined quarterly TRF solution shows no significant impact on the TRF, when the CRF is estimated consistently with the TRF. The root mean square value of the post-fit station coordinate residuals is 0.9 cm.     

Graph theory for analyzing pair-wise data: application to geophysical model parameters estimated from interferometric synthetic aperture radar data at Okmok volcano,Alaska

Graph theory is useful for analyzing time-dependent model parameters estimated from interferometric synthetic aperture radar (InSAR) data in the temporal domain. Plotting acquisition dates (epochs) as vertices and pair-wise interferometric combinations as edges defines an incidence graph. The edge-vertex incidence matrix and the normalized edge Laplacian matrix are factors in the covariance matrix for the pair-wise data. Using empirical measures of residual scatter in the pair-wise observations, we estimate the relative variance at each epoch by inverting the covariance of the pair-wise data. We evaluate the rank deficiency of the corresponding least-squares problem via the edge-vertex incidence matrix. We implement our method in a MATLAB software package called GraphTreeTA available on GitHub (https://github.com/feigl/gipht). We apply temporal adjustment to the data set described in Lu et al. (Geophys Res Solid Earth 110, 2005) at Okmok volcano, Alaska, which erupted most recently in 1997 and 2008. The data set contains 44 differential volumetric changes and uncertainties estimated from interferograms between 1997 and 2004. Estimates show that approximately half of the magma volume lost during the 1997 eruption was recovered by the summer of 2003. Between June 2002 and September 2003, the estimated rate of volumetric increase is \((6.2 \, \pm \, 0.6) \times 10^6~\mathrm{m}^3/\mathrm{year} \). Our preferred model provides a reasonable fit that is compatible with viscoelastic relaxation in the five years following the 1997 eruption. Although we demonstrate the approach using volumetric rates of change, our formulation in terms of incidence graphs applies to any quantity derived from pair-wise differences, such as range change, range gradient, or atmospheric delay.     

Continental hydrology loading observed by VLBI measurements

Variations in continental water storage lead to loading deformation of the crust with typical peak-to-peak variations at very long baseline interferometry (VLBI) sites of 3–15 mm in the vertical component and 1–2 mm in the horizontal component. The hydrology signal at VLBI sites has annual and semi-annual components and clear interannual variations. We have calculated the hydrology loading series using mass loading distributions derived from the global land data assimilation system (GLDAS) hydrology model and alternatively from a global grid of equal-area gravity recovery and climate experiment (GRACE) mascons. In the analysis of the two weekly VLBI 24-h R1 and R4 network sessions from 2003 to 2010 the baseline length repeatabilities are reduced in 79 % (80 %) of baselines when GLDAS (GRACE) loading corrections are applied. Site vertical coordinate repeatabilities are reduced in about 80 % of the sites when either GLDAS or GRACE loading is used. In the horizontal components, reduction occurs in 70–80 % of the sites. Estimates of the annual site vertical amplitudes were reduced for 16 out of 18 sites if either loading series was applied. We estimated loading admittance factors for each site and found that the average admittances were 1.01 \(\pm \) 0.05 for GRACE and 1.39 \(\pm \) 0.07 for GLDAS. The standard deviations of the GRACE admittances and GLDAS admittances were 0.31 and 0.68, respectively. For sites that have been observed in a set of sufficiently temporally dense daily sessions, the average correlation between VLBI vertical monthly averaged series and GLDAS or GRACE loading series was 0.47 and 0.43, respectively.     

Geographic information system–based morphometric characterization of sub-watersheds of Meenachil river basin,Kottayam district,Kerala, India

The study areas Tikovil and Payppara sub-watersheds of Meenachil river cover 158.9 and 111.9 km², respectively. These watersheds are parts of Western Ghats, which is an ecologically sensitive region. The drainage network of the sub-watersheds was delineated using SOI topographical maps on 1:50,000 scale using the Arc GIS software. The stream orders were calculated using the method proposed by Strahler's (1964 Strahler, A. N. 1964. “Quantitative geomorphology of drainage basins and channel networks”. In Hand book of applied hydrology. Vol. 4, Edited by: Chow, V. T. Vol. 4, 39–44.  [Google Scholar]). The drainage network shows that the terrain exhibits dendritic to sub-dendritic drainage pattern. Stream order ranges from the fifth to the sixth order. Drainage density varies between 1.69 and 2.62 km/km². The drainage texture of the drainage basins are 2.3 km^–1 and 6.98 km^–1 and categorized as coarse to very fine texture. Stream frequency is low in the case of Payappara sub-watershed (1.78 km^–2). Payappara sub-watershed has the highest constant of channel maintenance value of 0.59 indicating much fewer structural disturbances and fewer runoff conditions. The form factor value varies in between 0.42 and 0.55 suggesting elongated shape formed for Payappara sub-watershed and a rather more circular shape for Tikovil sub-watershed. The mean bifurcation ratio (3.5) indicates that both the sub-watersheds are within the natural stream system. Hence from the study it can be concluded that GIS techniques prove to be a competent tool in morphometric analysis.     

Earth orientation parameters estimated from VLBI during the CONT11 campaign

In this paper we investigate the accuracy of the earth orientation parameters (EOP) estimated from the continuous VLBI campaign CONT11. We first estimated EOP with daily resolution and compared these to EOP estimated from GNSS data. We find that the WRMS differences are about 31  $\upmu $ as for polar motion and 7  $\upmu $ s for length of day. This is about the precision we could expect, based on Monte Carlo simulations and the results of the previous CONT campaigns. We also estimated EOP with hourly resolution to study the sub-diurnal variations. The results confirm the results of previous studies, showing that the current IERS model for high-frequency EOP variations does not explain all the sub-diurnal variations seen in the estimated time series. We then compared our results to various empirical high-frequency EOP models. However, we did not find that any of these gave any unambiguous improvement. Several simulations testing the impact of various aspects of, e.g. the observing network were also made. For example, we made simulations assuming that all CONT11 stations were equipped with fast VLBI2010 antennas. We found that the WRMS error decreased by about a factor five compared to the current VLBI system. Furthermore, the simulations showed that it is very important to have a homogenous global distribution of the stations for achieving the highest precision for the EOP.     

Evaluation of gravitational curvatures of a tesseroid in spherical integral kernels

Proper understanding of how the Earth’s mass distributions and redistributions influence the Earth’s gravity field-related functionals is crucial for numerous applications in geodesy, geophysics and related geosciences. Calculations of the gravitational curvatures (GC) have been proposed in geodesy in recent years. In view of future satellite missions, the sixth-order developments of the gradients are becoming requisite. In this paper, a set of 3D integral GC formulas of a tesseroid mass body have been provided by spherical integral kernels in the spatial domain. Based on the Taylor series expansion approach, the numerical expressions of the 3D GC formulas are provided up to sixth order. Moreover, numerical experiments demonstrate the correctness of the 3D Taylor series approach for the GC formulas with order as high as sixth order. Analogous to other gravitational effects (e.g., gravitational potential, gravity vector, gravity gradient tensor), numerically it is found that there exist the very-near-area problem and polar singularity problem in the GC east–east–radial, north–north–radial and radial–radial–radial components in spatial domain, and compared to the other gravitational effects, the relative approximation errors of the GC components are larger due to not only the influence of the geocentric distance but also the influence of the latitude. This study shows that the magnitude of each term for the nonzero GC functionals by a grid resolution 15\(^{{\prime } }\,\times \) 15\(^{{\prime }}\) at GOCE satellite height can reach of about 10\(^{-16}\) m\(^{-1}\) s\(^{2}\) for zero order, 10\(^{-24 }\) or 10\(^{-23}\) m\(^{-1}\) s\(^{2}\) for second order, 10\(^{-29}\) m\(^{-1}\) s\(^{2}\) for fourth order and 10\(^{-35}\) or 10\(^{-34}\) m\(^{-1}\) s\(^{2}\) for sixth order, respectively.