Crustal Deformation of the Kendeng Fault Branches Area from GNSS and InSAR Data in Surabaya City, Indonesia

Abstract Slip rates on active faults derived from the interseismic velocity field are critical to understanding seismic hazards in metropolitan cities. This study integrated the data from the Global Navigation Satellite System (GNSS) and Interferometric Synthetic Aperture Radar (InSAR) to evaluate the interseismic velocities in the second-largest city in Indonesia, Surabaya, where branches of the Kendeng fault (the Surabaya and the Waru faults) pass through. Data from 16 campaign-mode GNSS stations collected between 2017 and 2020 from previous research are reprocessed to estimate the velocity field. Horizontal velocities under the ITRF frame range between -23.8 mm/yr and 47.9 mm/yr toward the southeast. Vertical velocities generally range between -1.3 mm/yr and -112.2 mm/yr. Sentinel-1A SAR data of both ascending and descending tracks acquired between November 2014 and July 2020 were used to generate the interferograms with the InSAR Scientific Computing Environment (ISCE) software. Furthermore, cumulative displacement time series were constructed using the Small BAseline Subset (SBAS) technique within the Generic InSAR Analysis Toolbox (GIAnT). This study also carried out the detection of outlier SAR epochs to improve the precision of Line-of-sight (LOS) velocity estimates. The LOS velocities range from -14.8 to 10.8 mm/yr in the ascending track and from -12.7 to 9.5 mm/yr in the descending track. These results will facilitate the detection of coupling behaviors on the Kendeng fault branches, which can improve our understanding of seismic risks in the Surabaya area.

Download Full-text

A strain rate model for Europe from a dense network of GNSS

10.5194/egusphere-egu21-11703 ◽

2021 ◽

Author(s):

Holger Steffen ◽

Rebekka Steffen ◽

Ambrus Kenyeres ◽

Alessandro Caporali ◽

Joaquin Zurutuza ◽

...

Keyword(s):

Velocity Field ◽

Strain Rate ◽

Velocity Model ◽

Satellite System ◽

Least Square ◽

Plate Boundaries ◽

Strain Rates ◽

Rate Model ◽

2D Strain ◽

Global Navigation Satellite

<p>Strain rates are an important factor to find areas that are under stress. Higher strain rates are usually observed along plate boundaries, while lower strain rates are found in intraplate regions. The increased availability of velocity solutions from Global Navigation Satellite System (GNSS) for entire Europe allows a 2D strain rate to be estimated at high resolution. Thus, regions of high and low strain become clearly visible. Here, we will present a new strain rate model, which is based on a recent velocity field solution by the EUREF Permanent Network Densification (EPND2100). This velocity field is obtained by the combination of weekly position SINEX solutions generated by 28 EPND Analysis Centres. More details on EPND can be found in the www.epnd.sgo-penc.hu website. The homogenized and quality checked velocity field is then interpolated via a least-square collocation using a fixed scale length of 135 km. In addition, the effect of known plate boundaries is considered during the interpolation to avoid a smoothing of nearby velocities on different tectonic plates. We also apply a moving variance approach to avoid effects of non-stationarity, which arise due to the variable station densities. The interpolated velocity model is then used to estimate a 2D strain rate covering most of Europe. We will highlight the situation in intraplate areas with very low strain rates but dense GNSS networks.</p>

Download Full-text

The March 2021 Damasi Earthquake Sequence, Central Greece: Reactivation Evidence across the Westward Propagating Tyrnavos Graben

Geosciences ◽

10.3390/geosciences11080328 ◽

2021 ◽

Vol 11 (8) ◽

pp. 328

Author(s):

Ioannis K. Koukouvelas ◽

Konstantinos G. Nikolakopoulos ◽

Aggeliki Kyriou ◽

Riccardo Caputo ◽

Alexandros Belesis ◽

...

Keyword(s):

Active Faults ◽

Satellite System ◽

Epicentral Area ◽

Surface Models ◽

Gnss Receivers ◽

Surface Ruptures ◽

Satellite Radar ◽

Global Navigation Satellite ◽

Geomorphological Analysis ◽

Complex Manner

On 3 March 2021, a strong shallow earthquake affected northern Thessaly, Greece, with an epicenter close to Damasi village causing significant destruction of many stone houses. In this contribution, we provide fieldwork observations, satellite radar interferometry, mapping of the active faults exposed in the epicentral area, liquefactions and coseismic surface ruptures, and preliminary geomorphological analyses of the epicentral area. The geomorphological analysis is based on air photographs, digital surface models analysis, Real-Time Kinematik (RTK) measurements with Global Navigation Satellite System (GNSS) receivers, and data from UAV flight campaigns. Although the seismotectonic setting of the area is complex and there is an apparent mismatch between field and interferometric data, the results of our investigations suggest that at least three fault segments were reactivated by the major shocks of the March seismic sequence. These tectonic structureslikely represent the westward propagation of the Tyrnavos Graben, where newly formed and inherited low-angle faults interplay in a complex manner.

Download Full-text

The characteristics of transforming coordinates from the geocentric system WGS-84 for the Mercator projection under low latitudes conditions

Geodesy and Cartography ◽

10.22389/0016-7126-2018-940-10-2-6 ◽

2018 ◽

Vol 940 (10) ◽

pp. 2-6

Author(s):

J.A. Younes ◽

M.G. Mustafin

Keyword(s):

Coordinate System ◽

Satellite System ◽

Programming Algorithm ◽

Low Latitude ◽

Model Calculations ◽

Mercator Projection ◽

Low Latitudes ◽

Rectangular Coordinates ◽

Global Navigation Satellite ◽

Coordination Methods

The issue of calculating the plane rectangular coordinates using the data obtained by the satellite observations during the creation of the geodetic networks is discussed in the article. The peculiarity of these works is in conversion of the coordinates into the Mercator projection, while the plane coordinate system on the base of Gauss-Kruger projection is used in Russia. When using the technology of global navigation satellite system, this task is relevant for any point (area) of the Earth due to a fundamentally different approach in determining the coordinates. The fact is that satellite determinations are much more precise than the ground coordination methods (triangulation and others). In addition, the conversion to the zonal coordinate system is associated with errors; the value at present can prove to be completely critical. The expediency of using the Mercator projection in the topographic and geodetic works production at low latitudes is shown numerically on the basis of model calculations. To convert the coordinates from the geocentric system with the Mercator projection, a programming algorithm which is widely used in Russia was chosen. For its application under low-latitude conditions, the modification of known formulas to be used in Saudi Arabia is implemented.

Download Full-text

Detecting Apples in the Wild: Potential for Harvest Quantity Estimation

Sustainability ◽

10.3390/su13148054 ◽

2021 ◽

Vol 13 (14) ◽

pp. 8054

Author(s):

Artur Janowski ◽

Rafał Kaźmierczak ◽

Cezary Kowalczyk ◽

Jakub Szulwic

Keyword(s):

Image Analysis ◽

Pilot Study ◽

Production Management ◽

Satellite System ◽

Morphological Operations ◽

Computational Performance ◽

Analysis Methods ◽

Manual Counting ◽

In The Wild ◽

Global Navigation Satellite

Knowing the exact number of fruits and trees helps farmers to make better decisions in their orchard production management. The current practice of crop estimation practice often involves manual counting of fruits (before harvesting), which is an extremely time-consuming and costly process. Additionally, this is not practicable for large orchards. Thanks to the changes that have taken place in recent years in the field of image analysis methods and computational performance, it is possible to create solutions for automatic fruit counting based on registered digital images. The pilot study aims to confirm the state of knowledge in the use of three methods (You Only Look Once—YOLO, Viola–Jones—a method based on the synergy of morphological operations of digital imagesand Hough transformation) of image recognition for apple detecting and counting. The study compared the results of three image analysis methods that can be used for counting apple fruits. They were validated, and their results allowed the recommendation of a method based on the YOLO algorithm for the proposed solution. It was based on the use of mass accessible devices (smartphones equipped with a camera with the required accuracy of image acquisition and accurate Global Navigation Satellite System (GNSS) positioning) for orchard owners to count growing apples. In our pilot study, three methods of counting apples were tested to create an automatic system for estimating apple yields in orchards. The test orchard is located at the University of Warmia and Mazury in Olsztyn. The tests were carried out on four trees located in different parts of the orchard. For the tests used, the dataset contained 1102 apple images and 3800 background images without fruits.

Download Full-text

GPS-BDS-Galileo double-differenced stochastic model refinement based on least-squares variance component estimation

Journal of Navigation ◽

10.1017/s0373463321000564 ◽

2021 ◽

pp. 1-16

Author(s):

Hong Hu ◽

Xuefeng Xie ◽

Jingxiang Gao ◽

Shuanggen Jin ◽

Peng Jiang

Keyword(s):

Stochastic Model ◽

Least Squares ◽

Stochastic Models ◽

Variance Component ◽

Satellite System ◽

Variance Component Estimation ◽

Short Baseline ◽

Component Estimation ◽

Zero Baseline ◽

Global Navigation Satellite

Abstract Stochastic models are essential for precise navigation and positioning of the global navigation satellite system (GNSS). A stochastic model can influence the resolution of ambiguity, which is a key step in GNSS positioning. Most of the existing multi-GNSS stochastic models are based on the GPS empirical model, while differences in the precision of observations among different systems are not considered. In this paper, three refined stochastic models, namely the variance components between systems (RSM1), the variances of different types of observations (RSM2) and the variances of observations for each satellite (RSM3) are proposed based on the least-squares variance component estimation (LS-VCE). Zero-baseline and short-baseline GNSS experimental data were used to verify the proposed three refined stochastic models. The results show that, compared with the traditional elevation-dependent model (EDM), though the proposed models do not significantly improve the ambiguity resolution success rate, the positioning precision of the three proposed models has been improved. RSM3, which is more realistic for the data itself, performs the best, and the precision at elevation mask angles 20°, 30°, 40°, 50° can be improved by 4⋅6%, 7⋅6%, 13⋅2%, 73⋅0% for L1-B1-E1 and 1⋅1%, 4⋅8%, 16⋅3%, 64⋅5% for L2-B2-E5a, respectively.

Download Full-text

Investigation of Pre-Earthquake Ionospheric and Atmospheric Disturbances for Three Large Earthquakes in Mexico

Geosciences ◽

10.3390/geosciences11010016 ◽

2020 ◽

Vol 11 (1) ◽

pp. 16

Author(s):

Christina Oikonomou ◽

Haris Haralambous ◽

Sergey Pulinets ◽

Aakriti Khadka ◽

Shukra R. Paudel ◽

...

Keyword(s):

Chemical Potential ◽

Atmospheric Model ◽

Satellite System ◽

Seismic Events ◽

Electron Content ◽

Total Electron ◽

Atmospheric Disturbances ◽

Global Ionospheric Maps ◽

Global Navigation Satellite ◽

Large Earthquakes

The purpose of the present study is to investigate simultaneously pre-earthquake ionospheric and atmospheric disturbances by the application of different methodologies, with the ultimate aim to detect their possible link with the impending seismic event. Three large earthquakes in Mexico are selected (8.2 Mw, 7.1 Mw and 6.6 Mw during 8 and 19 September 2017 and 21 January 2016 respectively), while ionospheric variations during the entire year 2017 prior to 37 earthquakes are also examined. In particular, Total Electron Content (TEC) retrieved from Global Navigation Satellite System (GNSS) networks and Atmospheric Chemical Potential (ACP) variations extracted from an atmospheric model are analyzed by performing statistical and spectral analysis on TEC measurements with the aid of Global Ionospheric Maps (GIMs), Ionospheric Precursor Mask (IPM) methodology and time series and regional maps of ACP. It is found that both large and short scale ionospheric anomalies occurring from few hours to a few days prior to the seismic events may be linked to the forthcoming events and most of them are nearly concurrent with atmospheric anomalies happening during the same day. This analysis also highlights that even in low-latitude areas it is possible to discern pre-earthquake ionospheric disturbances possibly linked with the imminent seismic events.

Download Full-text

Estimation of fractional cycle bias for GPS/BDS-2/Galileo based on international GNSS monitoring and assessment system observations using the uncombined PPP model

Satellite Navigation ◽

10.1186/s43020-021-00039-x ◽

2021 ◽

Vol 2 (1) ◽

Author(s):

Jin Wang ◽

Qin Zhang ◽

Guanwen Huang

Keyword(s):

Estimation Method ◽

Satellite Orbit ◽

Satellite System ◽

Assessment System ◽

Navigation Satellite ◽

Dual Frequency ◽

Fractional Cycle Bias ◽

The Stability ◽

High Consistency ◽

Global Navigation Satellite

AbstractThe Fractional Cycle Bias (FCB) product is crucial for the Ambiguity Resolution (AR) in Precise Point Positioning (PPP). Different from the traditional method using the ionospheric-free ambiguity which is formed by the Wide Lane (WL) and Narrow Lane (NL) combinations, the uncombined PPP model is flexible and effective to generate the FCB products. This study presents the FCB estimation method based on the multi-Global Navigation Satellite System (GNSS) precise satellite orbit and clock corrections from the international GNSS Monitoring and Assessment System (iGMAS) observations using the uncombined PPP model. The dual-frequency raw ambiguities are combined by the integer coefficients (4,− 3) and (1,− 1) to directly estimate the FCBs. The details of FCB estimation are described with the Global Positioning System (GPS), BeiDou-2 Navigation Satellite System (BDS-2) and Galileo Navigation Satellite System (Galileo). For the estimated FCBs, the Root Mean Squares (RMSs) of the posterior residuals are smaller than 0.1 cycles, which indicates a high consistency for the float ambiguities. The stability of the WL FCBs series is better than 0.02 cycles for the three GNSS systems, while the STandard Deviation (STD) of the NL FCBs for BDS-2 is larger than 0.139 cycles. The combined FCBs have better stability than the raw series. With the multi-GNSS FCB products, the PPP AR for GPS/BDS-2/Galileo is demonstrated using the raw observations. For hourly static positioning results, the performance of the PPP AR with the three-system observations is improved by 42.6%, but only 13.1% for kinematic positioning results. The results indicate that precise and reliable positioning can be achieved with the PPP AR of GPS/BDS-2/Galileo, supported by multi-GNSS satellite orbit, clock, and FCB products based on iGMAS.

Download Full-text

Improving DGNSS Performance through the Use of Network RTK Corrections

Remote Sensing ◽

10.3390/rs13091621 ◽

2021 ◽

Vol 13 (9) ◽

pp. 1621

Author(s):

Duojie Weng ◽

Shengyue Ji ◽

Yangwei Lu ◽

Wu Chen ◽

Zhihua Li

Keyword(s):

Carrier Phase ◽

Local Area ◽

Satellite System ◽

High Accuracy ◽

Single Frequency ◽

Baseline Length ◽

Low Latitude ◽

Phase Measurements ◽

Network Rtk ◽

Global Navigation Satellite

The differential global navigation satellite system (DGNSS) is an enhancement system that is widely used to improve the accuracy of single-frequency receivers. However, distance-dependent errors are not considered in conventional DGNSS, and DGNSS accuracy decreases when baseline length increases. In network real-time kinematic (RTK) positioning, distance-dependent errors are accurately modelled to enable ambiguity resolution on the user side, and standard Radio Technical Commission for Maritime Services (RTCM) formats have also been developed to describe the spatial characteristics of distance-dependent errors. However, the network RTK service was mainly developed for carrier-phase measurements on professional user receivers. The purpose of this study was to modify the local-area DGNSS through the use of network RTK corrections. Distance-dependent errors can be reduced, and accuracy for a longer baseline length can be improved. The results in the low-latitude areas showed that the accuracy of the modified DGNSS could be improved by more than 50% for a 17.9 km baseline during solar active years. The method in this paper extends the use of available network RTK corrections with high accuracy to normal local-area DGNSS applications.

Download Full-text

Galileo L10 Satellites: Orbit, Clock and Signal-in-Space Performance Analysis

Sensors ◽

10.3390/s21051695 ◽

2021 ◽

Vol 21 (5) ◽

pp. 1695

Author(s):

Constantin-Octavian Andrei ◽

Sonja Lahtinen ◽

Markku Poutanen ◽

Hannu Koivula ◽

Jan Johansson

Keyword(s):

Performance Indicators ◽

Satellite System ◽

Short Term ◽

Orbital Inclination ◽

Navigation Data ◽

Periodic Signals ◽

Global Navigation Satellite ◽

Minor Correction ◽

Very High

The tenth launch (L10) of the European Global Navigation Satellite System Galileo filled in all orbital slots in the constellation. The launch carried four Galileo satellites and took place in July 2018. The satellites were declared operational in February 2019. In this study, we report on the performance of the Galileo L10 satellites in terms of orbital inclination and repeat period parameters, broadcast satellite clocks and signal in space (SiS) performance indicators. We used all available broadcast navigation data from the IGS consolidated navigation files. These satellites have not been reported in the previous studies. First, the orbital inclination (56.7±0.15°) and repeat period (50680.7±0.22 s) for all four satellites are within the nominal values. The data analysis reveals also 13.5-, 27-, 177- and 354-days periodic signals. Second, the broadcast satellite clocks show different correction magnitude due to different trends in the bias component. One clock switch and several other minor correction jumps have occurred since the satellites were declared operational. Short-term discontinuities are within ±1 ps/s, whereas clock accuracy values are constantly below 0.20 m (root-mean-square—rms). Finally, the SiS performance has been very high in terms of availability and accuracy. Monthly SiS availability has been constantly above the target value of 87% and much higher in 2020 as compared to 2019. Monthly SiS accuracy has been below 0.20 m (95th percentile) and below 0.40 m (99th percentile). The performance figures depend on the content and quality of the consolidated navigation files as well as the precise reference products. Nevertheless, these levels of accuracy are well below the 7 m threshold (95th percentile) specified in the Galileo service definition document.

Download Full-text

Cloud Radar Observations of Diurnal and Seasonal Cloudiness over Reunion Island

Atmosphere ◽

10.3390/atmos12070868 ◽

2021 ◽

Vol 12 (7) ◽

pp. 868

Author(s):

Jonathan Durand ◽

Edouard Lees ◽

Olivier Bousquet ◽

Julien Delanoë ◽

François Bonnardot

Keyword(s):

Sea Breeze ◽

Satellite System ◽

Research Infrastructure ◽

Weather Data ◽

Reunion Island ◽

Wet Season ◽

Trade Winds ◽

Cloud Radar ◽

Réunion Island ◽

Global Navigation Satellite

In November 2016, a 95 GHz cloud radar was permanently deployed in Reunion Island to investigate the vertical distribution of tropical clouds and monitor the temporal variability of cloudiness in the frame of the pan-European research infrastructure Aerosol, Clouds and Trace gases Research InfraStructure (ACTRIS). In the present study, reflectivity observations collected during the two first years of operation (2016–2018) of this vertically pointing cloud radar are relied upon to investigate the diurnal and seasonal cycle of cloudiness in the northern part of this island. During the wet season (December–March), cloudiness is particularly pronounced between 1–3 km above sea level (with a frequency of cloud occurrence of 45% between 12:00–19:00 LST) and 8–12 km (with a frequency of cloud occurrence of 15% between 14:00–19:00 LST). During the dry season (June–September), this bimodal vertical mode is no longer observed and the vertical cloud extension is essentially limited to a height of 3 km due to both the drop-in humidity resulting from the northward migration of the ITCZ and the capping effect of the trade winds inversion. The frequency of cloud occurrence is at its maximum between 13:00–18:00 LST, with a probability of 35% at 15 LST near an altitude of 2 km. The analysis of global navigation satellite system (GNSS)-derived weather data also shows that the diurnal cycle of low- (1–3 km) and mid-to-high level (5–10 km) clouds is strongly correlated with the diurnal evolution of tropospheric humidity, suggesting that additional moisture is advected towards the island by the sea breeze regime. The detailed analysis of cloudiness observations collected during the four seasons sampled in 2017 and 2018 also shows substantial differences between the two years, possibly associated with a strong positive Indian Ocean Southern Dipole (IOSD) event extending throughout the year 2017.

Download Full-text