The EUREF Permanent Network Densification is a collaborative effort of 26 European GNSS analysis centers providing series of daily or weekly station position estimates of dense national and regional GNSS networks, in order to combine them into one homogenized set of station positions and velocities. During the combination, the station meta-data, including station names, DOMES numbers, and position offset definitions were carefully homogenized, position outliers were efficiently eliminated, and the results were cross-checked for any remaining inconsistencies. The results cover the period from March 1999 to January 2017 (GPS week 1000-1933) and include 31 networks with positions and velocities for 3192 stations, well covering Europe. The positions and velocities are expressed in ITRF2014 and ETRF2014 reference frames based on the Minimum Constraint approach using a selected set of ITRF2014 reference stations. The position alignment with the ITRF2014 is at the level of 1.5, 1.2, and 3.2 mm RMS for the East, North, Up components, respectively, while the velocity RMS values are 0.17, 0.14, and 0.38 mm/year for the East, North, and Up components, respectively. The high quality of the combined solution is also reflected by the 1.1, 1.1, and 3.5 mm weighted RMS values for the East, North, and Up components, respectively.

The seismic sequence of November 2019 in Albania culminating with the Mw = 6.4 event of 26 November 2019 was examined from the geodetic (InSAR and GNSS), structural, and historical viewpoints, with some ideas on possible areas of greater hazard. We present accurate estimates of the coseismic displacements using permanent GNSS stations active before and after the sequence, as well as SAR interferograms with Sentinel-1 in ascending and descending mode. When compared with the displacements predicted by a dislocation model on an elastic half space using the moment tensor information of a reverse fault mechanism, the InSAR and GNSS data fit at the mm level provided the hypocentral depth is set to 8 ± 2 km. Next, we examined the elastic stress generated by the Mw = 7.2 Montenegro earthquake of 1979, with the Albania 2019 event as receiver fault, to conclude that the Coulomb stress transfer, at least for the elastic component, was too small to have influenced the 2019 Albania event. A somewhat different picture emerges from the combined elastic deformation resulting after the two (1979 and 2019) events: we investigated the fault geometries where the Coulomb stress is maximized and concluded that the geometry with highest induced Coulomb stress, of the order of ca. 2–3 bar (0.2–0.3 MPa), is that of a vertical, dextral strike slip fault, striking SW to NE. This optimal receiver fault is located between the faults activated in 1979 and 2019, and very closely resembles the Lezhe fault, which marks the transition between the Dinarides and the Albanides.

The Central European GNSS Research Network (CEGRN) collects GNSS data since 1994 from contributors which today include 42 Institutions in 33 Countries. CEGRN returns a dataset of coordinates and velocities computed according to international standards and the most recent processing procedures and recommendations. We provide a dataset of 1229 positions and velocities resulting from 3 or more repetitions of coordinate measurements of each site over 4 or more years. The velocity data result from a combination of eight multiyear, partially overlapping networks, using 234 stations of class A of the European Permanent Network (EPN) for alignment to the 'European Fixed' ETRF2000 Reference Frame. The rms (root mean square) of the 8 individual contributions to the combined solution, after a 7 - parameter Helmert transformation, is less than 5 mm in the observation period 1996-2017. This combined CEGRN network maintains the origin coincident with that of the ETRF2000 reference frame to within 1.8 mm rms for the entire period of analysis. The mean positions and velocities of common EPN Class A and CEGRN stations differ by 0.0 ± 1.1, 0.5 ± 1.0 and 0.1 ± 2.7 mm for the coordinates and 0.06 ± 0.13, -0.07 ± 0.12, 0.38 ± 0.28 mm/yr for the velocities respectively for the North, East and Up components at epoch 2010.0.

To investigate the transfer of elastic energy between a regional stress field and a set of localized faults, we project the stress rate tensor inferred from the Italian GNSS (Global Navigation Satellite Systems) velocity field onto faults selected from the Database of Individual Seismogenic Sources (DISS 3.2.0). For given Lamé constants and friction coefficient, we compute the loading rate on each fault in terms of the Coulomb failure function (CFF) rate. By varying the strike, dip and rake angles around the nominal DISS values, we also estimate the geometry of planes that are optimally oriented for maximal CFF rate. Out of 86 Individual Seismogenic Sources (ISSs), all well covered by GNSS data, 78–81 (depending on the assumed friction coefficient) load energy at a rate of 0–4 kPa yr⁻¹. The faults displaying larger CFF rates (4–6 ± 1 kPa yr⁻¹) are located in the central Apennines and are all characterized by a significant strike-slip component. We also find that the loading rate of 75% of the examined sources is less than 1 kPa yr⁻¹ lower than that of optimally oriented faults. We also analysed 2016 August 24 and October 30 central Apennines earthquakes (M_w 6.0–6.5, respectively). The strike of their causative faults based on seismological and tectonic data and the geodetically inferred strike differ by <30°. Some sources exhibit a strike oblique to the direction of maximum strain rate, suggesting that in some instances the present-day stress acts on inherited faults. The choice of the friction coefficient only marginally affects this result.

Here we test the capability of the Broadcast Ephemeris Message, in both its GPS-like (Keplerian ellipse with secular and periodic perturbations) and Glonass-like (numerical integration of a 9D state vector) formats, to reproduce a corresponding precise ephemeris. We start from a daily Rinex 3.04 navigation file for multiple GNSS and the corresponding SP3 precise orbits computed by CNES (Centre National d’Etudes Spatiales) for GPS, Glonass, Galileo and CODE (Center for Orbit Determination in Europe) for Beidou, and compute broadcast ECEF coordinates and clocks. The pre-fit discrepancies are converted by least squares to corrections to the broadcast ephemeris parameters in two-hour consecutive arcs (for GPS, Galileo and Beidou) and to a set of seven Helmert parameters for the entire day, to align in origin, orientation and scale to the common GNSS IGS14 Reference Frame. The test cases suggest that the Broadcast Ephemeris Message, complemented with Reference Frame information, can reproduce the precise ephemeris and clocks with centimetric accuracy for intervals at least equal to the respective validity times, typically 2 h. The broadcast ephemeris of Glonass consists of three initial positions and velocities at epoch, three constant Lunisolar accelerations for the satellite position, and of three polynomial coefficients for the satellite clock. The 9D vector of state is numerically integrated to generate position and velocity data within the validity time (0.5 h) of the message. To test the capability of this model to reproduce the corresponding values of a precise ephemeris, the 9D vector of state and clock polynomials are adjusted until the rms (root mean squared spread) of the post-fit residuals relative to a precise orbit (CNES’s in our case) is minimum. We show in one example (one satellite for one day) that the Glonass type of message can reproduce a precise ephemeris and clock with a rms spread of 0.025 m over one-hour arcs. Volume computations on one month of data with all available satellites confirm the test results. For GPS, Glonass, Galileo and Beidou, the best fitting clock values predicted by our second order polynomials, based on a 15 min sampling, are shown to fit the corresponding high rate clocks (30 s sampling) of MGEX with zero bias and a rms spread of 0.062 ns (GPS G01), 0.023 ns (Galileo E01), 0.43 ns (Glonass R01), 0.086 ns (Beidou C07) and 0.086 ns (Beidou C12). Modifications to the GPS-like message structure and Glonass algorithm are proposed to increase the validity time by including the effect of the 3rd zonal harmonic of the Earth’s gravity field. The potential of the RTCM messages for broadcasting the improved navigation message is reviewed.