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.     

Stochastic models in the DORIS position time series: estimates for IDS contribution to ITRF2014

This paper focuses on the investigation of the deterministic and stochastic parts of the Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) weekly time series aligned to the newest release of ITRF2014. A set of 90 stations was divided into three groups depending on when the data were collected at an individual station. To reliably describe the DORIS time series, we employed a mathematical model that included the long-term nonlinear signal, linear trend, seasonal oscillations and a stochastic part, all being estimated with maximum likelihood estimation. We proved that the values of the parameters delivered for DORIS data are strictly correlated with the time span of the observations. The quality of the most recent data has significantly improved. Not only did the seasonal amplitudes decrease over the years, but also, and most importantly, the noise level and its type changed significantly. Among several tested models, the power-law process may be chosen as the preferred one for most of the DORIS data. Moreover, the preferred noise model has changed through the years from an autoregressive process to pure power-law noise with few stations characterised by a positive spectral index. For the latest observations, the medians of the velocity errors were equal to 0.3, 0.3 and 0.4 mm/year, respectively, for the North, East and Up components. In the best cases, a velocity uncertainty of DORIS sites of 0.1 mm/year is achievable when the appropriate coloured noise model is taken into consideration.     

A grid-based tropospheric product for China using a GNSS network

Tropospheric delay accounts for one source of error in global navigation satellite systems (GNSS). To better characterize the tropospheric delays in the temporal and spatial domain and facilitate the safety-critical use of GNSS across China, a method is proposed to generate a grid-based tropospheric product (GTP) using the GNSS network with an empirical tropospheric model, known as IGGtrop. The prototype system generates the GTPs in post-processing and real-time modes and is based on the undifferenced and uncombined precise point positioning (UU-PPP) technique. GTPs are constructed for a grid form (\(2.0{^{\circ }}\times 2.5{^{\circ }}\) latitude–longitude) over China with a time resolution of 5 min. The real-time GTP messages are encoded in a self-defined RTCM3 format and broadcast to users using NTRIP (networked transport of RTCM via internet protocol), which enables efficient and safe transmission to real-time users. Our approach for GTP generation consists of three sequential steps. In the first step, GNSS-derived zenith tropospheric delays (ZTDs) for a network of GNSS stations are estimated using UU-PPP. In the second step, vertical adjustments for the GNSS-derived ZTDs are applied to address the height differences between the GNSS stations and grid points. The ZTD height corrections are provided by the IGGtrop model. Finally, an inverse distance weighting method is used to interpolate the GNSS-derived ZTDs from the surrounding GNSS stations to the location of the grid point. A total of 210 global positioning system (GPS) stations from the crustal movement observation network of China are used to generate the GTPs in both post-processing and real-time modes. The accuracies of the GTPs are assessed against with ERA-Interim-derived ZTDs and the GPS-derived ZTDs at 12 test GPS stations, respectively. The results show that the post-processing and real-time GTPs can provide the ZTDs with accuracies of 1.4 and 1.8 cm, respectively. We also apply the GTPs in real-time kinematic GPS PPP, and the results show that the convergence time of the PPP solutions is shortened. These results confirm that the GTPs can act as an efficient information source to augment GNSS positioning over China.     

Determination of Phobos’ rotational parameters by an inertial frame bundle block adjustment

A functional model for a bundle block adjustment in the inertial reference frame was developed, implemented and tested. This approach enables the determination of rotation parameters of planetary bodies on the basis of photogrammetric observations. Tests with a self-consistent synthetic data set showed that the implementation converges reliably toward the expected values of the introduced unknown parameters of the adjustment, e.g., spin pole orientation, and that it can cope with typical observational errors in the data. We applied the model to a data set of Phobos using images from the Mars Express and the Viking mission. With Phobos being in a locked rotation, we computed a forced libration amplitude of \(1.14^\circ \pm 0.03^\circ \) together with a control point network of 685 points.     

What governs the presence of residual vegetation in boreal wildfires?

Wildfires are frequent boreal forest disturbances in Canada, and emulating their patterns with forest harvesting has emerged as a common forest management goal. Wildfires contain many patches of residual vegetation of various size, shape, and composition; understanding their characteristics provides insights for improved emulation criteria. We studied the occurrence of residual vegetation within eleven boreal wildfire events in a natural setting; fires ignited by lightning, no suppression efforts, and no prior anthropogenic disturbances. Relative importance of the measurable geo-environmental factors and their marginal effects on residual presence are studied using Random Forests. These factors included distance from natural firebreaks (wetland, bedrock and non-vegetated areas, and water), land cover, and topographic variables (elevation, slope, and ruggedness index). We present results at spatial resolutions ranging from four to 64 m while emphasizing four and 32 m since they mimic IKONOS- and Landsat-type images. Natural firebreak features, especially the proximity to wetlands, are among the most important variables that explain the likelihood residual occurrence. The majority of residual vegetation areas are concentrated within 100 m of wetlands. Topographic variables, typically important in rugged terrain, are less important in explaining the presence of residuals within our study fires.     

Geographical impacts on social networks from perspectives of space and place: an empirical study using mobile phone data

Space and place are two fundamental concepts in geography. Geographical factors have long been known as drivers of many aspects of people’s social networks. But whether and how space and place affect social networks differently are still unclear. The widespread use of location-aware devices provides a novel source for distinguishing the mechanisms of their impacts on social networks. Using mobile phone data, this paper explores the effects of space and place on social networks. From the perspective of space, we confirm the distance decay effect in social networks, based on a comparison between synthetic social ties generated by a null model and actual social ties derived from real-world data. From the perspective of place, we introduce several measures to evaluate interactions between individuals and inspect the trio relationship including distance, spatio-temporal co-occurrence, and social ties. We found that people’s interaction is a more important factor than spatial proximity, indicating that the spatial factor has a stronger impact on social networks in place compared to that in space. Furthermore, we verify the hypothesis that interactions play an important role in strengthening friendships.     

Spatial representation in the social interaction potential metric: an analysis of scale and parameter sensitivity

The social interaction potential (SIP) metric measures urban structural constraints on social interaction opportunities of a metropolitan region based on the time geographic concept of joint accessibility. Previous implementations of the metric used an interaction surface based on census tracts and the locations of their centroids. This has been shown to be a shortcoming, as the metric strongly depends on the scale of the zoning system in the region, making it difficult to compare the SIP metric between metropolitan regions. This research explores the role of spatial representation in the SIP metric and identifies a suitable grid-based representation that allows for comparison between regions while retaining cost-effectiveness with respect to computational burden. We also report on findings from an extensive sensitivity analysis investigating the SIP metric’s input parameters such as a travel flow congestion factor and the length of the allowable time budget for social activities. The results provide new insights on the role of the modifiable areal unit problem in the computation of time geographic measures of accessibility.     

Weak GPS signal tracking using FFT discriminator in open loop receiver

An open loop tracking architecture, which tracks GPS signals under weak and challenging conditions, is analyzed. The in-phase and quadrature-phase integration pair is regarded as a single tone complex signal. An FFT-based method is used as a frequency discriminator to estimate the Doppler frequency residual of the single tone signal. Another FFT-based method applies complex squaring to eliminate the effect of the navigation data bits polarities. The performance of the FFT-based discriminators is assessed in three criteria. Those criteria are the signal strength and dynamic range that can be tracked and the accuracy of the estimated Doppler frequency. In addition, the performance of the discriminators is analyzed to provide the theoretical and simulated peak detection probability. The results indicate that the FFT discriminator can track signals about 5 dB weaker than the signals that can be tracked by the complex squared FFT discriminator. In a quasi-static environment, the Doppler frequency residual can be assumed to be around zero, which can enable the FFT-based discriminators to track signals with approximately 2 dB less power. Moreover, the performance of the FFT-based discriminators is compared with the performance of two other frequency discriminators, namely the fast–slow and power-based. The comparison results indicate that these two frequency discriminators give higher frequency estimation accuracy, but they have a narrower dynamic range.     

Precise orbit determination of BeiDou constellation: method comparison

Chinese BeiDou navigation satellite system is in official service as a regional constellation with five geostationary earth orbit (GEO) satellites, five inclined geosynchronous satellite orbit (IGSO) satellites and four medium earth orbit (MEO) satellites. There are mainly two methods for precise orbit determination of the BeiDou constellation found in the current literatures. One is the independent single-system method, where only BeiDou observations are used without help from other GNSS systems. The other is the two-step GPS-assisted method where in the first step, GPS data are used to resolve some common parameters, such as station coordinates, receiver clocks and zenith tropospheric delay parameters, which are then introduced as known quantities in BeiDou processing in the second step. We conduct a thorough performance comparison between the two methods. Observations from the BeiDou experimental tracking stations and the IGS Multi-GNSS Experiment network from January 1 to March 31, 2013, are processed with the Positioning and Navigation Data Analyst (PANDA) software. The results show that for BeiDou IGSO and MEO satellites, the two-step GPS-assisted method outperforms the independent single-system method in both internal orbit overlap precision and external satellite laser ranging validation. For BeiDou GEO satellites, the two methods show close performances. Zenith tropospheric delays estimated from the first method are very close to those estimated from GPS precise point positioning in the second method, with differences of several millimeters. Satellite clock estimates from the two methods show similar performances when assessing the stability of the BeiDou on board clocks.