A two-step ionospheric modeling algorithm considering the impact of GLONASS pseudo-range inter-channel biases

The Global Navigation Satellite System presents a plausible and cost-effective way of computing the total electron content (TEC). But TEC estimated value could be seriously affected by the differential code biases (DCB) of frequency-dependent satellites and receivers. Unlike GPS and other satellite systems, GLONASS adopts a frequency-division multiplexing access mode to distinguish different satellites. This strategy leads to different wavelengths and inter-frequency biases (IFBs) for both pseudo-range and carrier phase observations, whose impacts are rarely considered in ionospheric modeling. We obtained observations from four groups of co-stations to analyze the characteristics of the GLONASS receiver P1P2 pseudo-range IFB with a double-difference method. The results showed that the GLONASS P1P2 pseudo-range IFB remained stable for a period of time and could catch up to several meters, which cannot be absorbed by the receiver DCB during ionospheric modeling. Given the characteristics of the GLONASS P1P2 pseudo-range IFB, we proposed a two-step ionosphere modeling method with the priori IFB information. The experimental analysis showed that the new algorithm can effectively eliminate the adverse effects on ionospheric model and hardware delay parameters estimation in different space environments. During high solar activity period, compared to the traditional GPS + GLONASS modeling algorithm, the absolute average deviation of TEC decreased from 2.17 to 2.07 TECu (TEC unit); simultaneously, the average RMS of GPS satellite DCB decreased from 0.225 to 0.219 ns, and the average deviation of GLONASS satellite DCB decreased from 0.253 to 0.113 ns with a great improvement in over 55%.     

The polar form of the spherical harmonic spectrum: implications for filtering <Emphasis Type="BoldSmallCaps">grace</Emphasis> data

Representing the spherical harmonic spectrum of a field on the sphere in terms of its amplitude and phase is termed as its polar form. In this study, we look at how the amplitude and phase are affected by linear low-pass filtering. The impact of filtering on amplitude is well understood, but that on phase has not been studied previously. Here, we demonstrate that a certain class of filters only affect the amplitude of the spherical harmonic spectrum and not the phase, but the others affect both the amplitude and phase. Further, we also demonstrate that the filtered phase helps in ascertaining the efficacy of decorrelation filters used in the grace community.     

Tropospheric refractivity and zenith path delays from least-squares collocation of meteorological and GNSS data

Precise positioning requires an accurate a priori troposphere model to enhance the solution quality. Several empirical models are available, but they may not properly characterize the state of troposphere, especially in severe weather conditions. Another possible solution is to use regional troposphere models based on real-time or near-real time measurements. In this study, we present the total refractivity and zenith total delay (ZTD) models based on a numerical weather prediction (NWP) model, Global Navigation Satellite System (GNSS) data and ground-based meteorological observations. We reconstruct the total refractivity profiles over the western part of Switzerland and the total refractivity profiles as well as ZTDs over Poland using the least-squares collocation software COMEDIE (Collocation of Meteorological Data for Interpretation and Estimation of Tropospheric Pathdelays) developed at ETH Zürich. In these two case studies, profiles of the total refractivity and ZTDs are calculated from different data sets. For Switzerland, the data set with the best agreement with the reference radiosonde (RS) measurements is the combination of ground-based meteorological observations and GNSS ZTDs. Introducing the horizontal gradients does not improve the vertical interpolation, and results in slightly larger biases and standard deviations. For Poland, the data set based on meteorological parameters from the NWP Weather Research and Forecasting (WRF) model and from a combination of the NWP model and GNSS ZTDs shows the best agreement with the reference RS data. In terms of ZTD, the combined NWP-GNSS observations and GNSS-only data set exhibit the best accuracy with an average bias (from all stations) of 3.7 mm and average standard deviations of 17.0 mm w.r.t. the reference GNSS stations.     

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.     

Group delay variations of GPS transmitting and receiving antennas

GPS code pseudorange measurements exhibit group delay variations at the transmitting and the receiving antenna. We calibrated C1 and P2 delay variations with respect to dual-frequency carrier phase observations and obtained nadir-dependent corrections for 32 satellites of the GPS constellation in early 2015 as well as elevation-dependent corrections for 13 receiving antenna models. The combined delay variations reach up to 1.0 m (3.3 ns) in the ionosphere-free linear combination for specific pairs of satellite and receiving antennas. Applying these corrections to the code measurements improves code/carrier single-frequency precise point positioning, ambiguity fixing based on the Melbourne–Wübbena linear combination, and determination of ionospheric total electron content. It also affects fractional cycle biases and differential code biases.     

Error analysis of the 3D similarity coordinate transformation

The 3D similarity coordinate transformation with the Gauss–Helmert error model is investigated. The first-order error analysis of an analytical least-squares solution to this problem is developed in detail. While additive errors are assumed in the translation and scale estimates, a 3 × 1 multiplicative error vector is defined to effectively parameterize the rotation matrix estimation error. The propagation of the errors in the coordinate measurements to the errors in the estimated transformation parameters is derived step-by-step, and the formulae for calculating the variance–covariance matrix of the estimated parameters are presented.     

Global ionosphere map constructed by using total electron content from ground-based GNSS receiver and FORMOSAT-3/COSMIC GPS occultation experiment

Effects of rapidly changing ionospheric weather are critical in high accuracy positioning, navigation, and communication applications. A system used to construct the global total electron content (TEC) distribution for monitoring the ionospheric weather in near-real time is needed in the modern society. Here we build the TEC map named Taiwan Ionosphere Group for Education and Research (TIGER) Global Ionospheric Map (GIM) from observations of ground-based GNSS receivers and space-based FORMOSAT-3/COSMIC (F3/C) GPS radio occultation observations using the spherical harmonic expansion and Kalman filter update formula. The TIGER GIM (TGIM) will be published in near-real time of 4-h delay with a spatial resolution of 2.5° in latitude and 5° in longitude and a high temporal resolution of every 5 min. The F3/C TEC results in an improvement on the GIM of about 15.5%, especially over the ocean areas. The TGIM highly correlates with the GIMs published by other international organizations. Therefore, the routinely published TGIM in near-real time is not only for communication, positioning, and navigation applications but also for monitoring and scientific study of ionospheric weathers, such as magnetic storms and seismo-ionospheric anomalies.     

A modified geometry- and ionospheric-free combination for static three-carrier ambiguity resolution

The reliability of the classical geometry- and ionospheric-free (GIF) three-carrier ambiguity resolution (TCAR) degrades when applied to long baselines of hundreds of kilometers. To overcome this deficiency, we propose two new models, which are used sequentially to resolve wide-lane (WL) and narrow-lane (NL) ambiguities and form a stepwise ambiguity resolution (AR) strategy. In the first model, after a successful extra-wide-lane AR, the pseudorange and phase observations are combined to estimate WL ambiguities, in which the residual ionospheric delays and geometry effects are eliminated. In the second model, using the resolved ambiguities from the first step, the two WL ambiguities are combined to remove ionospheric and geometry effects. The unknown coefficients in the two models are determined in such that they minimize the formal errors in the ambiguity estimates to optimize the ambiguity estimation. Using experimental BeiDou triple-frequency observations, we evaluate our method and identify three advantages. First, the two models use double-differenced phase observations that are not differences across frequency. Second, the two models are entirely free from ionospheric delay and geometry effects. Third, the unknown estimates in the two models satisfy the minimum noise condition, which makes the formal errors in the float NL ambiguity estimates much lower than those obtained with common GIF TCAR methods, thereby directly and significantly increasing the success rate of AR compared to the cascaded integer resolution method and two other GIF combinations.     

Multi‐scale approach to mining significant spatial co‐location patterns

Spatial co‐location pattern mining aims to discover a collection of Boolean spatial features, which are frequently located in close geographic proximity to each other. Existing methods for identifying spatial co‐location patterns usually require users to specify two thresholds, i.e. the prevalence threshold for measuring the prevalence of candidate co‐location patterns and distance threshold to search the spatial co‐location patterns. However, these two thresholds are difficult to determine in practice, and improper thresholds may lead to the misidentification of useful patterns and the incorrect reporting of meaningless patterns. The multi‐scale approach proposed in this study overcomes this limitation. Initially, the prevalence of candidate co‐location patterns is measured statistically by using a significance test, and a non‐parametric model is developed to construct the null distribution of features with the consideration of spatial auto‐correlation. Next, the spatial co‐location patterns are explored at multi‐scales instead of single scale (or distance threshold) discovery. The validity of the co‐location patterns is evaluated based on the concept of lifetime. Experiments on both synthetic and ecological datasets show that spatial co‐location patterns are discovered correctly and completely by using the proposed method; on the other hand, the subjectivity in discovery of spatial co‐location patterns is reduced significantly.     

A Novel Approach of Despeckling SAR Images Using Nonlocal Means Filtering

Synthetic aperture radar (SAR) is a day and night, all weather satellite imaging technology. Inherent property of SAR image is speckle noise which produces granular patterns in the image. Speckle noise occurs due to the interference of backscattered echo from earth’s rough surface. There are various speckle reduction techniques in spatial domain and transform domain. Non local means filtering (NLMF) is the technique used for denoising which uses Gaussian weights. In NLMF algorithm, the filtering is performed by taking the weighted mean of all the pixels in a selected search area. The weight given to the pixel is based on the similarity measure calculated as the weighted Euclidean distance over the two windows. Non local means filtering smoothes out homogeneous areas but edges are not preserved. So a discontinuity adaptive weight is used in order to preserve heterogeneous areas like edges. This technique is called as discontinuity adaptive non local means filtering and is well-adapted and robust in the case of Additive White Gaussian Noise (AWGN) model. But speckle is a multiplicative random noise and hence Euclidean distance is not a good choice. This paper presents evaluation results of using different distance measures for improving the accuracy of the Non local means filtering technique. The results are verified using real and synthetic images and from the results it can be concluded that the usage of Manhattan distance improves the accuracy of NLMF technique. Non local approach is used as a preprocessing or post processing technique for many denoising algorithms. So improving NLMF technique would help improving many of the existing denoising techniques.