Mariusz Figurski

Results of the application of tropospheric corrections from different troposphere models for precise GPS rapid static positioning

In many surveying applications, determination of accurate heights is of significant interest. The delay caused by the neutral atmosphere is one of the main factors limiting the accuracy of GPS positioning and affecting mainly the height coordinate component rather than horizontal ones. Estimation of the zenith total delay is a commonly used technique for accounting for the tropospheric delay in static positioning. However, in the rapid static positioning mode the estimation of the zenith total delay may fail, since for its reliable estimation longer observing sessions are required. In this paper, several troposphere modeling techniques were applied and tested with three processing scenarios: a single baseline solution with various height differences and a multi-baseline solution. In specific, we introduced external zenith total delays obtained from Modified Hopfield troposphere model with standard atmosphere parameters, UNB3m model, COAMPS numerical weather prediction model and zenith total delays interpolated from a reference network solution. The best results were obtained when tropospheric delays derived from the reference network were applied.

On the Handling of Outliers in the GNSS Time Series by Means of the Noise and Probability Analysis

The data pre-analysis plays a significant role in the noise determination. The most important issue is to find an optimum criterion for outliers removal, since their existence can affect any further analysis. The noises in the GNSS time series are characterized by spectral index and amplitudes that can be determined with a few different methods. In this research, the Maximum Likelihood Estimation (MLE) was used. The noise amplitudes as well as spectral indices were obtained for the topocentric coordinates with daily changes from few selected EPN (EUREF Permanent Network) stations. The data were obtained within the EPN re-processing made by the Military University of Technology Local Analysis Centre (MUT LAC). The outliers were removed from the most noisy 12 EPN stations with the criteria of 3 and 5 times the standard deviations (3σ, 5σ) as well as Median Absolute Deviation (MAD) to investigate how they affect noise parameters. The results show that the removal of outliers is necessary before any further analysis, otherwise one may obtain quite odd and unrealistic values. The probability analysis with skewness and kurtosis was also performed beyond the noise analysis. The values of skewness and kurtosis show that assuming a wrong criterion of outliers removal leads to the wrong results in case of probability distribution. On the basis of the results, we propose to use the MAD method for the outliers removal in the GNSS time series.

Error analysis for European IGS stations

Each of the GPS time series that describes the changes of topocentric components consists of a deterministic and a stochastic part, whose character influences the errors of the deterministic parameters. As to the uncertainties of reliable velocities of permanent satellite station systems, surveys that estimate and take into account any dependencies that may affect subsequent operational efficiency are very important. For this analysis, we used 42 stations from the IGS (International GNSS Service) network from Europe, processed at the Military University of Technology EUREF Permanent Network Local Analysis Centre (MUT LAC). The deterministic part of the GPS time series was removed using the least squares method. The seasonal periods in topocentric components were determined assuming the existence of the residual Chandler oscillation (1.67 cpy), as well as the annual tropical (1 cpy) and draconitic (1.04 cpy) oscillations with their harmonics up to 4th. We assumed the character of the residue as a combination of white and powerlaw noise. The obtained results show, that in the case of the European sub-network of IGS stations we are dealing with the coloured noise between white and flicker noise with the amplitudes between 3 to 6 mm/year-k/4 for horizontal components and between 6 to 15 mm/year-κ/4 for the vertical ones, where κ is a spectral index. Finally, we showed that the amplitudes and spectral indices of noise are reduced after performing a spatio-temporal filtering. All the elicited results referred to the uncertainties of velocities by estimating them before and after filtration and the simulation of their values for different lengths of the time series.

Investigation of the 16-year and 18-year ZTD Time Series Derived from GPS Data Processing

The GPS system can play an important role in activities related to the monitoring of climate. Long time series, coherent strategy, and very high quality of tropospheric parameter Zenith Tropospheric Delay (ZTD) estimated on the basis of GPS data analysis allows to investigate its usefulness for climate research as a direct GPS product. This paper presents results of analysis of 16-year time series derived from EUREF Permanent Network (EPN) reprocessing performed by the Military University of Technology. For 58 stations Lomb-Scargle periodograms were performed in order to obtain information about the oscillations in ZTD time series. Seasonal components and linear trend were estimated using Least Square Estimation (LSE) and Mann—Kendall trend test was used to confirm the presence of a linear trend designated by LSE method. In order to verify the impact of the length of time series on trend value, comparison between 16 and 18 years were performed.

GPS-derived height changes in diurnal and sub-diurnal timescales

This paper describes the research concerning precise short-time GPS solutions conducted in the Centre of Applied Geomatics, Military University of Technology, Warsaw, Poland. The data from ASG-EUPOS (Polish Active Geodetic Network) was processed using Bernese 5.0 software and EPN (EUREF Permanent Network) standards and models. In this study, the adapted 3-hour observation window is shifted every hour for obtaining hourly geocentric coordinates in ITRF2005 reference frame. The adjusted network consisted of over 130 stations from Poland and the neighbouring countries, the period covered observations collected from June 2008 through June 2010. These two years of observations allowed for examining short-period oscillations which were found to be closely related to the tidal (dynamic) frequencies. The analysis of the residua from the IERS2003 tidal model was performed using the least squares method with the Eterna software. It confirmed the existence of significant (several millimetres) oscillations in the frequencies corresponding to S1, K1 and K2. The paper describes the idea of data processing and analysis, presents the results of vertical (the Up component) oscillations in main tidal frequency bands, and also includes discussion of possible explanation of the existence of short-period oscillations in the GPS precise solutions and the possibility of propagation of short-period oscillations into long-period spurious changes in the daily (standard) GPS solutions.

A Bridge Deflection Monitoring with GPS

A Bridge Deflection Monitoring with GPS This paper introduces results of investigation carried on by The Applied Geomatics Section in Military University of Technology. Research includes possibilities of monitoring dynamic behavior of a bridge using high rate GPS data. Whole event was executed with collaboration of The Road and Bridge Management and The Warsaw Geodesy Company. Interdisciplinary approach with this project allows authors to get reliable information about investigating constructions and their respond for true traffic loading detected by GPS receivers. Way of compute data and used software (TRACK) are also shown in this paper.

Efficient Usage of Dense GNSS Networks in Central Europe for the Visualization and Investigation of Ionospheric TEC Variations

The technique of the orthogonal projection of ionosphere electronic content variations for mapping total electron content (TEC) allows us to visualize ionospheric irregularities. For the reconstruction of global ionospheric characteristics, numerous global navigation satellite system (GNSS) receivers located in different regions of the Earth are used as sensors. We used dense GNSS networks in central Europe to detect and investigate a special type of plasma inhomogeneities, called travelling ionospheric disturbances (TID). Such use of GNSS sensors allows us to reconstruct the main TID parameters, such as spatial dimensions, velocities, and directions of their movement. The paper gives examples of the restoration of dynamic characteristics of ionospheric irregularities for quiet and disturbed geophysical conditions. Special attention is paid to the dynamics of ionospheric disturbances stimulated by the magnetic storms of two St. Patrick’s Days (17 March 2013 and 2015). Additional opportunities for the remote sensing of the ionosphere with the use of dense regional networks of GNSS receiving sensors have been noted too.

Irregular variations in GPS time series by probability and noise analysis

Abstract The character of the topocentric components in ETRF2000(R08) from the Polish ASG-EUPOS system was analysed using skewness and kurtosis derived from the data probability density function. The data from 115 permanent GPS stations with a time span of more than 5 years were used. The main goal of this research was to show that any unmodelled systematics can disrupt the results. The obtained median values of skewness and kurtosis clearly indicate the discrepancies between the assumed normality of the GPS time series distribution and the reality, mainly due to stochastic and/or deterministic parts that are still present in the data. The quadratic relationship between skewness and kurtosis was developed with the empirically determined constants. The noise analysis with the Maximum Likelihood Estimation was also performed with the assumptions of white and white plus power-law noise. The estimated spectral indices for power law noises are close to −1 (flicker noise). The uncertainties of the intraplate velocities with white and white plus power law noise assumptions were calculated. It was shown that these uncertainties can be underestimated up to 5 mm/year. Additionally we made 1000 simulation which were aimed at showing how the values of skewness and kurtosis are changed when some mismodelled part of the data remains in the time series. We assumed different values of trend (1 and 6 mm/year), seasonal amplitudes (1 and 6 mm) and offsets (3 and 10 mm) with flicker noise of 1 mm amplitude to see how the metrics we use are biased.

Investigation on Tidal Components in GPS Coordinates

Investigation on Tidal Components in GPS Coordinates This paper presents analyses on the GPS coordinates from sub-diurnal solutions of EPN data provided by Warsaw Military University of Technology. The aim of this research is to investigate the way the tidal models used in Bernese software (solid Earth and ocean tides as well) fit to the individual conditions of EPN stations. The 1-hour solution technique of GPS data processing was utilized to obtain coordinates of above 70 EPN stations. Additionally several Polish permanent sites with clearly seen oscillations were examined. This processing technique allowed us to recognize diurnal and sub-diurnal residual oscillations which could be next utilized for validation of the tidal models.

ASG-EUPOS Monitoring with Reference to EPN

ASG-EUPOS Monitoring with Reference to EPN This paper concerns methods and results of processing of data collected by ASG-EUPOS (Active Geodetic Network - European Position Determination System) sites since the beginning of the system. The elaboration was made by Centre of Applied Geomatics from Military University of Technology and it was based on the experience gained during EPN (European Permanent Network) processing. Strategy was similar to strategy used by EPN Local Analysis Centers for European network elaborations and EPN test reprocessing (Kenyeres et al., 2008). As the outcome daily and weekly solutions (North, East, Up component for each station) were obtained. Basing on weekly solutions local velocities and vertical movements were determined. Preliminary tests of combining national and regional networks were also made and results were compared to official EPN solutions (differences below 1 cm). Daily solutions were analyzed in the context of stability and reliability basing on time series standard deviation.