scholarly journals Common mode signals and vertical velocities in the great Alpine area from GNSS data

2021 ◽  
Author(s):  
Francesco Pintori ◽  
Enrico Serpelloni ◽  
Adriano Gualandi

Abstract. We study time series of vertical ground displacements from continuous GNSS stations to investigate the spatial and temporal contribution of different geophysical processes to the time-varying displacements that are superimposed on vertical linear trends across the European Alps. We apply a multivariate statistics-based blind source separation algorithm to both GNSS displacement time series and to ground displacements associated with atmospheric and hydrological loading processes, as obtained from global reanalysis models. This allows us to associate each retrieved geodetic vertical deformation signal with a corresponding forcing process. Atmospheric loading is the most important one, reaching amplitudes larger than 2 cm. Besides atmospheric loading, seasonal displacements with amplitudes of about 1 cm are associated with temperature-related processes and with hydrological loading. We find that both temperature and hydrological loading cause peculiar spatial features of GNSS ground displacements. For example, temperature-related seasonal displacements show different behaviour at sites in the plains and in the mountains. Atmospheric and hydrological loading, besides the first-order spatially uniform feature, are associated also with NS and EW displacement gradients. We filter out signals associated with non-tectonic deformation from the raw time series to study their impact on both the estimated noise and linear rates in the vertical direction. While the impact on rates appears rather limited, given also the long-time span of the time-series considered in this work, the uncertainties estimated from filtered time-series assuming a power law + white noise model are significantly reduced, with an important increase in white noise contributions to the total noise budget. Finally, we present the filtered velocity field and show how vertical ground velocities are positively correlated with topographic features of the Alps.

2021 ◽  
Vol 13 (3) ◽  
pp. 431
Author(s):  
Yuefan He ◽  
Guigen Nie ◽  
Shuguang Wu ◽  
Haiyang Li

The displacement of Global Navigation Satellite System (GNSS) station contains the information of surface elastic deformation caused by the variation of land water reserves. This paper selects the long-term coordinate series data of 671 International GNSS Service (IGS) reference stations distributed globally under the framework of World Geodetic System 1984 (WGS84) from 2000 to 2021. Different noise model combinations are used for noise analysis, and the optimal noise model for each station before and after hydrologic loading correction is calculated. The results show that the noise models of global IGS reference stations are diverse, and each component has different optimal noise model characteristics, mainly white noise + flicker noise (WN+FN), generalized Gauss–Markov noise (GGM) and white noise + power law noise (WN+PL). Through specific analysis between the optimal noise model and the time series velocity of the station, it is found that the maximum influence value of the vertical velocity can reach 1.8 mm when hydrological loading is considered. Different complex noise models also have a certain influence on the linear velocity and velocity uncertainty of the station. Among them, the influence of white noise + random walking noise is relatively obvious, and its maximum influence value in the elevation direction can reach over 2 mm/year. When studying the impact of hydrological loading correction on the periodicity of the coordinate series, it is concluded whether the hydrological loading is calculated or not, and the GNSS long-term coordinate series has obvious annual and semi-annual amplitude changes, which are most obvious in the vertical direction, up to 16.48 mm.


Author(s):  
F. C. Çomut ◽  
A. Ustun ◽  
M. Lazecky ◽  
M. M. Aref

The SAR Interferometry (InSAR) application has shown great potential in monitoring of land terrain changes and in detection of land deformations such as subsidence. Longer time analysis can lead to understand longer trends and changes. Using different bands of SAR satellite (C- from ERS 1-2 and Envisat, L- from ALOS) over the study area, we achieve knowledge of movements in long-term and evaluation of its dynamic changes within observed period of time. Results from InSAR processing fit with the position changes in vertical direction based on GPS network established over the basin as an effective geodetic network. Time series (StaMPS PS+SB) of several points over Çumra County in eastern part of Konya City show a general trend of the deformation that is expected to be approximately between -13 to -17 mm/year. Northern part of Karaman is affected by faster subsidence, borders of the subsidence trough were identified from Envisat. <br><br> Presenting InSAR results together with GIS information about locations and time of occurrence of sudden subsidence, urban/industrial growth in time and climate changes helps in better understanding of the situation. This way, the impact of natural and man-made changes will be shown for urban planning thanks to InSAR and GIS comparisons with hydrogeological modeling. In this study we present results of differential and multitemporal InSAR series using different bands and GIS conjunction associated with seasonal and temporal groundwater level changes in Konya Closed Basin.


2017 ◽  
Vol 10 (9) ◽  
pp. 3589-3607 ◽  
Author(s):  
Jan Dousa ◽  
Pavel Vaclavovic ◽  
Michal Elias

Abstract. In this paper, we present results of the second reprocessing of all data from 1996 to 2014 from all stations in International Association of Geodesy (IAG) Reference Frame Sub-Commission for Europe (EUREF) Permanent Network (EPN) as performed at the Geodetic Observatory Pecný (GOP). While the original goal of this research was to ultimately contribute to the realization of a new European Terrestrial Reference System (ETRS), we also aim to provide a new set of GNSS (Global Navigation Satellite System) tropospheric parameter time series with possible applications to climate research. To achieve these goals, we improved a strategy to guarantee the continuity of these tropospheric parameters and we prepared several variants of troposphere modelling. We then assessed all solutions in terms of the repeatability of coordinates as an internal evaluation of applied models and strategies and in terms of zenith tropospheric delays (ZTDs) and horizontal gradients with those of the ERA-Interim numerical weather model (NWM) reanalysis. When compared to the GOP Repro1 (first EUREF reprocessing) solution, the results of the GOP Repro2 (second EUREF reprocessing) yielded improvements of approximately 50 and 25 % in the repeatability of the horizontal and vertical components, respectively, and of approximately 9 % in tropospheric parameters. Vertical repeatability was reduced from 4.14 to 3.73 mm when using the VMF1 mapping function, a priori ZHD (zenith hydrostatic delay), and non-tidal atmospheric loading corrections from actual weather data. Raising the elevation cut-off angle from 3 to 7° and then to 10° increased RMS from coordinates' repeatability, which was then confirmed by independently comparing GNSS tropospheric parameters with the NWM reanalysis. The assessment of tropospheric horizontal gradients with respect to the ERA-Interim revealed a strong sensitivity of estimated gradients to the quality of GNSS antenna tracking performance. This impact was demonstrated at the Mallorca station, where gradients systematically grew up to 5 mm during the period between 2003 and 2008, before this behaviour disappeared when the antenna at the station was changed. The impact of processing variants on long-term ZTD trend estimates was assessed at 172 EUREF stations with time series longer than 10 years. The most significant site-specific impact was due to the non-tidal atmospheric loading followed by the impact of changing the elevation cut-off angle from 3 to 10°. The other processing strategy had a very small or negligible impact on estimated trends.


2021 ◽  
Author(s):  
Pauline André ◽  
Marie-Pierre Doin ◽  
Marguerite Mathey ◽  
Swann Zerathe ◽  
Riccardo Vassallo ◽  
...  

&lt;p&gt;Based on geomorphological criteria, large-scale slow gravitational deformation affecting entire mountain flank, often being referred as Deep-Seated Gravitational Slope Deformation (DSGSD), have been shown to affect most of the reliefs worldwide. For instance in the European Alps, these deformation patterns were identified in several areas such as the Aosta Valley (Martinotti et al., 2011) or the Mercantour massif (Jomard, 2006). DSGSD inventories based on visual interpretation of scarps and field mapping were then compiled (e.g. Crosta et al., 2013) revealing the widespread occurrence of DSGSD. However, many aspects of these large-scale gravitational processes remain unclear and in particular their present-day activity and temporal evolution remain largely unknown.&lt;/p&gt;&lt;p&gt;The present study aims at characterizing the spatial extent of DSGSD, and their velocity, at the scale of Western Alps through InSAR time series analysis using NSBAS processing chain (Doin et al., 2001). We used the whole SAR Sentinel-1 archive, between 2014 and 2018, with an acquisition every 6 days, on an ascending track. The processing was adapted to fit the specific conditions of the Alps (seasonal snow cover, strong local relief, vegetation and strong atmospheric heterogeneities). In particular we implemented a correction using the ERA 5 weather model and we used snow masks in winter allowing to select long temporal baseline interferograms with as little snow as possible. As we specifically aim to study deformation patterns at the scale of valley flanks, an average high-pass filter on moving subwindows has been applied to the interferograms prior to the implementation of time-serie inversions. This step strongly reduced the impact of residual atmospheric delays.&lt;/p&gt;&lt;p&gt;The resulting velocity map in the line of sight (LOS) of the satellite reveals ubiquitous gravitational deformation patterns over the whole Western Alps, with localized patches of moving slopes showing sharp discontinuities with stable surrounding areas. We used radar geometry and InSAR measurement quality factors as indicators to identify the most trusted areas and to extract an inventory of potential DSGSD with their spatial extent. Doing so, we identified more than two thousands slowly deforming areas characterized by LOS velocities from 4 to 20 mm/year. We then compared the geometries of our &amp;#8220;InSAR-detected-deforming-slopes&amp;#8221; with previously published DSGSD inventories. Good agreements were found for example in the Aosta valley where most of the deforming areas from our velocity map are falling into the DSGSD outlines of Crosta et al. (2013). Currently, we continue to investigate the potential of this large-scale velocity map for DSGSD understanding and we plan to use artificial intelligence to search for possible generic properties between the detected sites.&lt;/p&gt;


Author(s):  
Richard McCleary ◽  
David McDowall ◽  
Bradley J. Bartos

Chapter 3 introduces the Box-Jenkins AutoRegressive Integrated Moving Average (ARIMA) noise modeling strategy. The strategy begins with a test of the Normality assumption using a Kolomogov-Smirnov (KS) statistic. Non-Normal time series are transformed with a Box-Cox procedure is applied. A tentative ARIMA noise model is then identified from a sample AutoCorrelation function (ACF). If the sample ACF identifies a nonstationary model, the time series is differenced. Integer orders p and q of the underlying autoregressive and moving average structures are then identified from the ACF and partial autocorrelation function (PACF). Parameters of the tentative ARIMA noise model are estimated with maximum likelihood methods. If the estimates lie within the stationary-invertible bounds and are statistically significant, the residuals of the tentative model are diagnosed to determine whether the model’s residuals are not different than white noise. If the tentative model’s residuals satisfy this assumption, the statistically adequate model is accepted. Otherwise, the identification-estimation-diagnosis ARIMA noise model-building strategy continues iteratively until it yields a statistically adequate model. The Box-Jenkins ARIMA noise modeling strategy is illustrated with detailed analyses of twelve time series. The example analyses include non-Normal time series, stationary white noise, autoregressive and moving average time series, nonstationary time series, and seasonal time series. The time series models built in Chapter 3 are re-introduced in later chapters. Chapter 3 concludes with a discussion and demonstration of auxiliary modeling procedures that are not part of the Box-Jenkins strategy. These auxiliary procedures include the use of information criteria to compare models, unit root tests of stationarity, and co-integration.


2020 ◽  
Vol 12 (6) ◽  
pp. 915 ◽  
Author(s):  
Benjamin Brede ◽  
Jochem Verrelst ◽  
Jean-Philippe Gastellu-Etchegorry ◽  
Jan G. P. W. Clevers ◽  
Leo Goudzwaard ◽  
...  

The European Space Agency (ESA)’s Sentinel-2A (S2A) mission is providing time series that allow the characterisation of dynamic vegetation, especially when combined with the National Aeronautics and Space Administration (NASA)/United States Geological Survey (USGS) Landsat 7 (L7) and Landsat 8 (L8) missions. Hybrid retrieval workflows combining non-parametric Machine Learning Regression Algorithms (MLRAs) and vegetation Radiative Transfer Models (RTMs) were proposed as fast and accurate methods to infer biophysical parameters such as Leaf Area Index (LAI) from these data streams. However, the exact design of optimal retrieval workflows is rarely discussed. In this study, the impact of five retrieval workflow features on LAI prediction performance of MultiSpectral Instrument (MSI), Enhanced Thematic Mapper Plus (ETM+) and Operational Land Imager (OLI) observations was analysed over a Dutch beech forest site for a one-year period. The retrieval workflow features were the (1) addition of prior knowledge of leaf chemistry (two alternatives), (2) the choice of RTM (two alternatives), (3) the addition of Gaussian noise to RTM produced training data (four and five alternatives), (4) possibility of using Sun Zenith Angle (SZA) as an additional MLRA training feature (two alternatives), and (5) the choice of MLRA (six alternatives). The features were varied in a full grid resulting in 960 inversion models in order to find the overall impact on performance as well as possible interactions among the features. A combination of a Terrestrial Laser Scanning (TLS) time series with litter-trap derived LAI served as independent validation. The addition of absolute noise had the most significant impact on prediction performance. It improved the median prediction Root Mean Square Error (RMSE) by 1.08 m2 m−2 when 5% noise was added compared to inversions with 0% absolute noise. The choice of the MLRA was second most important in terms of median prediction performance, which differed by 0.52 m2 m−2 between the best and worst model. The best inversion model achieved an RMSE of 0.91 m2 m−2 and explained 84.9% of the variance of the reference time series. The results underline the need to explicitly describe the used noise model in future studies. Similar studies should be conducted in other study areas, both forest and crop systems, in order to test the noise model as an integral part of hybrid retrieval workflows.


2016 ◽  
Author(s):  
Anna Klos ◽  
Addisu Hunegnaw ◽  
Felix Norman Teferle ◽  
Kibrom Ebuy Abraha ◽  
Furqan Ahmed ◽  
...  

Abstract. Zenith Total Delay (ZTD) time series, derived from the re-processing of Global Positioning System (GPS) data, provide valuable information for the evaluation of global atmospheric reanalysis products such as ERA-Interim. Identifying the correct noise characteristics in the ZTD time series is an important step to assess the "true" magnitude of ZTD trend uncertainties. The ZTD residual time series for 1995–2015 are generated from our homogeneously re-processed and homogenized GPS time series from over 700 globally distributed stations classified into five major climate zones. The annual peak of ZTD data ranges between 10 and 150 mm with the smallest values for the polar and Alpine zone. The amplitudes of daily curve fall between 0 and 12 mm with the greatest variations for the dry zone. The autoregressive process of fourth order plus white noise model were found to be optimal for ZTD series. The tropical zone has the largest amplitude of autoregressive noise (9.59 mm) and the greatest amplitudes of white noise (13.00 mm). All climate zones have similar median coefficients of AR(1) (0.80 ± 0.05) with a minimum for polar and Alpine, which has the highest coefficients of AR(2) (0.27 ± 0.01) and AR(3) (0.11 ± 0.01) and clearly different from the other zones considered. We show that 53 of 120 examined trends became insignificant, when the optimum noise model was employed, compared to 11 insignificant trends for pure white noise. The uncertainty of the ZTD trends may be underestimated by a factor of 3 to 12 compared to the white noise only assumption.


2015 ◽  
Vol 5 (1) ◽  
Author(s):  
M. A. Goudarzi ◽  
M. Cocard ◽  
R. Santerre

AbstractWe analyzed the noise characteristics of 112 continuously operating GPS stations in eastern North America using the Spectral Analysis and the Maximum Likelihood Estimation (MLE) methods. Results of both methods show that the combination ofwhite plus flicker noise is the best model for describing the stochastic part of the position time series. We explored this further using the MLE in the time domain by testing noise models of (a) powerlaw, (b)white, (c)white plus flicker, (d)white plus randomwalk, and (e) white plus flicker plus random-walk. The results show that amplitudes of all noise models are smallest in the north direction and largest in the vertical direction. While amplitudes of white noise model in (c–e) are almost equal across the study area, they are prevailed by the flicker and Random-walk noise for all directions. Assuming flicker noise model increases uncertainties of the estimated velocities by a factor of 5–38 compared to the white noise model.


2018 ◽  
Vol 24 (4) ◽  
pp. 545-563
Author(s):  
Christian Gonzalo Pilapanta Amagua ◽  
Claudia Pereira Krueger ◽  
Alfonso Rodrigo Tierra Criollo

Abstract It is well known that daily estimates of GPS coordinates are highly temporally correlated and that the knowledge and understanding of this correlation allows to establish more realistic uncertainties of the parameters estimated from the data. Despite this, there are currently no studies related to the analysis and calculation of the noise sources in geodetic time series in Brazil. In this context, this paper focuses on the investigation of the stochastic properties of a total of 486 coordinates time series from 159 GPS stations belonging to the Brazilian Network for Continuous Monitoring of GNSS (RBMC) using the maximum likelihood estimation approach. To reliably describe the GPS time series, we evaluate 4 possible stochastic models as models of each time series: 3 models with integer spectral indices (white noise, flicker plus white noise and random-walk plus white noise model) and 1 with fractional spectral index (fractional power-law plus white noise). By comparing the calculated noise content values for each model, it is possible to demonstrate a stepwise increase of the noise content, being the combination of a fractional power-law process and white noise process, the model with smaller values and the combination of random walk process with white noise process, the model with greater values. The analysis of the spatial distribution of the noise values of the processes allow demonstrate that the GPS sites with the highest accumulated noise values, coincide with sites located in coastal zones and river basins and that their stochastic properties can be aliased by the occurrence of different physical signals typical of this type of zones, as the case of the hydrological loading effect.


2019 ◽  
Vol 11 (4) ◽  
pp. 386 ◽  
Author(s):  
Wenhao Li ◽  
Fei Li ◽  
Shengkai Zhang ◽  
Jintao Lei ◽  
Qingchuan Zhang ◽  
...  

The common mode error (CME) and optimal noise model are the two most important factors affecting the accuracy of time series in regional Global Navigation Satellite System (GNSS) networks. Removing the CME and selecting the optimal noise model can effectively improve the accuracy of GNSS coordinate time series. The CME, a major source of error, is related to the spatiotemporal distribution; hence, its detrimental effects on time series can be effectively reduced through spatial filtering. Independent component analysis (ICA) is used to filter the time series recorded by 79 GPS stations in Antarctica from 2010 to 2018. After removing stations exhibiting strong local effects using their spatial responses, the filtering results of residual time series derived from principal component analysis (PCA) and ICA are compared and analyzed. The Akaike information criterion (AIC) is then used to determine the optimal noise model of the GPS time series before and after ICA/PCA filtering. The results show that ICA is superior to PCA regarding both the filter results and the consistency of the optimal noise model. In terms of the filtering results, ICA can extract multisource error signals. After ICA filtering, the root mean square (RMS) values of the residual time series are reduced by 14.45%, 8.97%, and 13.27% in the east (E), north (N), and vertical (U) components, respectively, and the associated speed uncertainties are reduced by 13.50%, 8.06% and 11.82%, respectively. Furthermore, different GNSS time series in Antarctica have different optimal noise models with different noise characteristics in different components. The main noise models are the white noise plus flicker noise (WN+FN) and white noise plus power law noise (WN+PN) models. Additionally, the spectrum index of most PN is close to that of FN. Finally, there are more stations with consistent optimal noise models after ICA filtering than there are after PCA filtering.


Sign in / Sign up

Export Citation Format

Share Document