# Geophysics (physics.geo-ph)

• We investigate the accuracy and robustness of one of the most common methods used in glaciology for the discretization of the $\mathfrak{p}$-Stokes equations: equal order finite elements with Galerkin Least-Squares (GLS) stabilization. Furthermore we compare the results to other stabilized methods. We find that the vertical velocity component is more sensitive to the choice of GLS stabilization parameter than horizontal velocity. Additionally, the accuracy of the vertical velocity component is especially important since errors in this component can cause ice surface instabilities and propagate into future ice volume predictions. If the element cell size is set to the minimum edge length and the stabilization parameter is allowed to vary non-linearly with viscosity, the GLS stabilization parameter found in literature is a good choice on simple domains. However, near ice margins the standard parameter choice may result in significant oscillations in the vertical component of the surface velocity. For these cases, other stabilization techniques, such as the interior penalty method, result in better accuracy and are less sensitive to the choice of the stabilization parameter. During this work we also discovered that the manufactured solutions often used to evaluate errors in glaciology are not reliable due to high artificial surface forces at singularities. We perform our numerical experiments in both FEniCS and Elmer/Ice.
• Many measurements of geophysical interest are contaminated by tidal effects, which can have amplitudes greater than the ones in the events intended to be investigated. One of these measurements is the precise monitoring of the Earth's gravity field, regarded as a powerful tool for investigation of seismically induced deformations, hydrological effects on the ground and ice cap variations. Superconducting gravimeters have sensitivity below the 1nm/s$^2$ level, able to investigate those phenomena, and data from instruments in 33 locations globally distributed are available from the International Geodynamics and Earth Tides Service. However, fluctuations due to the Earth tides, with amplitudes typically in the order of 2000nm/s$^2$, dominate the signals, and must be removed before proceeding to the specific applications. In this work, we compare, for 5 locations(8 instruments) representative of the network, three different conceptual approaches that could be used for tide removal: frequency filtering, physical modelling and data-based modelling. Each approach applied to gravity time-series has shown a different limitation to be considered depending on the intended application. Results reveal that vestiges of tides remain in the residues for the modelling procedures, whereas the signal was distorted in the filtering procedures. The classical method with the best performance was the least-squares data-based modelling. However, the validation tests were not fulfilled, issue commonly related to overfitting. Although the tides could not be completely eliminated by the methods, they were sufficiently reduced to allow observation of geophysical events of interest above the 10nm/s$^2$ level, exemplified by a hydrology-related event of 60nm/s$^2$. The implementations adopted for each conceptual approach are general, so that the principles could be applied to other kinds of signals affected by tides.
• Isostatic equilibrium is commonly understood to be the state of equilibrium--neglecting mantle dynamics and the slow relaxation of the crust--achieved when there are no lateral gradients in hydrostatic or lithostatic pressure, and thus no lateral flow, at depth within the lower viscosity mantle that underlies the outer crust of a planetary body. In a constant-gravity Cartesian framework, this definition is equivalent to the requirement that columns of equal width contain equal masses. Here we show, however, that this equivalence breaks down when the spherical geometry of the problem is taken into account. Imposing the 'equal masses' requirement in a spherical geometry, as is commonly done in the literature, leads to significant lateral pressure gradients along internal equipotential surfaces, and thus corresponds to a state of disequilibrium. Compared with the 'equal pressures' model we present here, the 'equal masses' model always leads to an overestimate of the compensation depth. The magnitude of the discrepancy depends on the density structure of the body and the wavelength of the relevant topography, and is most pronounced when the compensation depth is a substantial fraction of the body's radius. Compared with the 'equal pressures' model, we show that analyses incorporating the 'equal masses' model may overestimate crustal thicknesses by as much as ~27% in the case of the lunar highlands, by ~10% in the case of the Martian highlands, and by nearly a factor of two in the case of Saturn's small icy moon Enceladus.
• Earthquakes cannot be predicted with precision, but algorithms exist for intermediate-term middle range prediction of main shocks above a pre-assigned threshold, based on seismicity patterns. Few years ago, a first attempt was made in the framework of project SISMA, funded by Italian Space Agency, to jointly use seismological tools, like CN algorithm and scenario earthquakes, and geodetic methods and techniques, like GPS and SAR monitoring, in order to effectively constrain priority areas where to concentrate prevention and seismic risk mitigation. We present a further development of integration of seismological and geodetic information, clearly showing the contribution of geodesy to the understanding and prediction of earthquakes. As a relevant application, the seismic crisis that started in Central Italy in August 2016 is considered in a retrospective analysis. Differently from the much more common approach, here GPS data are not used to estimate the standard 2D velocity and strain field in the area, but to reconstruct the velocity and strain pattern along transects, which are properly oriented according to the a priori information about the known tectonic setting. Overall, the analysis of the available geodetic data indicates that it is possible to highlight the velocity variation and the related strain accumulation in the area of Amatrice event, within the area alarmed by CN since November 1st, 2012. The considered counter examples, across CN alarmed and not-alarmed areas, do not show any comparable spatial acceleration localized trend. Therefore, we show that the combined analysis of the results of CN prediction algorithms, with those from the processing of adequately dense and permanent GNSS network data, may allow the routine highlight in advance of the strain accumulation. Thus it is possible to significantly reduce the size of the CN alarmed areas.