Impact of global seismicity on sea level change assessment

We analyze the effect of seismic activity on sealevel variations, by computing the time-dependent vertical crustal movement and geoid change due to coseismic deformations and postseismic relaxation effects. Seismic activity can affect both the absolute sealevel, by changing the Earth gravity field and hence the geoid height, and the relative sealevel, i.e. the radial distance between seafloor and geoid level. By using comprehensive seismic catalogues we assess the net effect of seismicity on tidal relative sealevel measurements as well as on the global oceanic surfaces, and we obtain an estimate of absolute sealevel variations of seismic origin. Our results confirm that, on a global scale, most of the signal is associated with few giant thrust events, and that RSL estimates obtained using tide-gauge data can be sensibly affected by the seismic driven sealevel signal. The recent measures of sealevel obtained by satellite altimetry show a wide regional variation of sealevel trends over the oceanic surfaces, with the largest deviations from the mean trend occurring in tectonically active regions. While our estimates of average absolute sealevel variations turn out to be orders of magnitude smaller than the satellite measured variations, we can still argue that mass redistribution associated with aseismic tectonic processes may contribute to the observed regional variability of sealevel variations.Comment: 34 pages, submitted to Journal of Geophysical Researc

Evolution, nucleosynthesis and yields of AGB stars at different metallicities (III): intermediate mass models, revised low mass models and the ph-FRUITY interface

We present a new set of models for intermediate mass AGB stars (4.0, 5.0 and, 6.0 Msun) at different metallicities (-2.15<=Fe/H]<=+0.15). This integrates the existing set of models for low mass AGB stars (1.3<=M/M<=3.0) already included in the FRUITY database. We describe the physical and chemical evolution of the computed models from the Main Sequence up to the end of the AGB phase. Due to less efficient third dredge up episodes, models with large core masses show modest surface enhancements. The latter is due to the fact that the interpulse phases are short and, then, Thermal Pulses are weak. Moreover, the high temperature at the base of the convective envelope prevents it to deeply penetrate the radiative underlying layers. Depending on the initial stellar mass, the heavy elements nucleosynthesis is dominated by different neutron sources. In particular, the s-process distributions of the more massive models are dominated by the \nean~reaction, which is efficiently activated during Thermal Pulses. At low metallicities, our models undergo hot bottom burning and hot third dredge up. We compare our theoretical final core masses to available white dwarf observations. Moreover, we quantify the weight that intermediate mass models have on the carbon stars luminosity function. Finally, we present the upgrade of the FRUITY web interface, now also including the physical quantities of the TP-AGB phase of all the models included in the database (ph-FRUITY).Comment: Accepted for publication on ApJ

Core-mantle boundary deformations and J2 variations resulting from the 2004 Sumatra earthquake

The deformation at the core-mantle boundary produced by the 2004 Sumatra earthquake is investigated by means of a semi-analytic theoretical model of global coseismic and postseismic deformation, predicting a millimetric coseismic perturbation over a large portion of the core-mantle boundary. Spectral features of such deformations are analysed and discussed. The time-dependent postseismic evolution of the elliptical part of the gravity field (J2) is also computed for different asthenosphere viscosity models. Our results show that, for asthenospheric viscosities smaller than 10^18 Pa s, the postseismic J2 variation in the next years is expected to leave a detectable signal in geodetic observations.Comment: 14 pages, 8 figures, 1 table. It will appear in Geophysical Journal Internationa

A Sea Level Equation for seismic perturbations

Large earthquakes are a potentially important source of relative sea level variations, since they can drive global deformation and simultaneously perturb the gravity field of the Earth. For the first time, we formalize a gravitationally self-consistent, integral sea level equation suitable for earthquakes, in which we account both for direct effects by the seismic dislocation and for the feedback from water loading associated with sea level changes. Our approach builds upon the well-established theory first proposed in the realm of glacio-isostatic adjustment modelling. The seismic sea level equation is numerically implemented to model sea level signals following the 2004 Sumatra–Andaman earthquake, showing that surface loading from ocean water redistribution (so far ignored in post-seismic deformation modelling) may account for a significant fraction of the total computed post-seismic sea level variatio

New insights on the Messina 1908 seismic source from post-seismic sea level change

The identification of a source model for the catastrophic 1908 December 28 Messina earth- quake (Mw = 7.2) has been the subject of many papers in the last decades. Several authors proposed different models on the basis of seismological, macroseismic and geodetic data sets; among these models, remarkable differences exist with regard to almost all parameters. We selected a subset of six models among those most cited in literature and used them to model the post-seismic sea level variation recorded at the tide gauge station of Messina (until 1923), to attempt an independent discrimination among them. For each model, we assumed a simple rheological structure and carried out a direct-search inversion of upper crust thickness and lower crust viscosity to fit the post-seismic sea level signal. This approach enabled us to iden- tify a class of fault geometries which is consistent with the post-seismic signal at the Messina tide gauge and with the known structural and rheological features of the Messina strai

FEMSA: A finite element simulation tool for quasi-static seismic deformation

Modeling postseismic deformation is an increasingly valuable tool in earthquake seismology. In particular, the Finite Element (FE) numerical method allows accurate modeling of complex faulting geometry, inhomogeneous materials and realistic viscous flow, appearing an excellent tool to investigate a lot of specific phenomena related with earthquakes

FEMSA: a finite element simulation tool for quasi-static seismic deformation modeling

We set up a computational tool to numerically model static and quasi-static deformation generated by faulting sources embedded in plane or spherical domains. We use a Finite Element (FE) approach to automatically implement arbitrary faulting sources and calculate displacement and stress fields induced by slip on the fault. The package makes use of the capabilities of CalculiX, a non commercial FE software designed to solve field problems (see <http://www.calculix.de> for details), and is freely distributed by request

The effects of a revised 7^7Be e−^--capture rate on solar neutrino fluxes

The electron-capture rate on $^7$Be is the main production channel for $^7$Li in several astrophysical environments. Theoretical evaluations have to account for not only the nuclear interaction, but also the processes in the plasma where $^7$Be ions and electrons interact. In the past decades several estimates were presented, pointing out that the theoretical uncertainty in the rate is in general of few percents. In the framework of fundamental solar physics, we consider here a recent evaluation for the $^7$Be+e$^-$ rate, not used up to now in the estimate of neutrino fluxes. We analysed the effects of the new assumptions on Standard Solar Models (SSMs) and compared the results obtained by adopting the revised $^7$Be+e$^-$ rate to those obtained by the one reported in a widely used compilation of reaction rates (ADE11). We found that new SSMs yield a maximum difference in the efficiency of the $^7$Be channel of about -4\% with respect to what is obtained with the previously adopted rate. This fact affects the production of neutrinos from $^8$B, increasing the relative flux up to a maximum of 2.7\%. Negligible variations are found for the physical and chemical properties of the computed solar models. The agreement with the SNO measurements of the neutral current component of the $^8$B neutrino flux is improved.Comment: 7 pages, 3 figures, 4 tables. Accepted for the publication on A&

Finite Element Modeling of the 2004 giant Sumatra Earthquake Postseismic Displacement Field

The 26 December 2004 Sumatra-Andaman earthquake is one of the largest earthquakes ever recorded since 1900. The earthquake resulted from complex slip on the fault where the oceanic portion of the Indian Plate slides under the Eurasian Plate, by the Indonesian Island of Sumatra. The particular features of the detected quasi-static displacement field has been previously attributed to the heterogeneous distribution of moment release on the fault plane. In the present work, we use a new computational FEM strategy to model the co- and postseismic displacement field associated with the Sumatra earthquake. For the first time we can study the joint effects of sphericity and 3D mechanical and rheological heterogeneities on the investigated observables. The comparison between our synthetic results and the available deformation data allows us to ascertain if also lateral heterogeneities in the physical properties of the medium could have played a role in assessing the deformation field

Trade--Off Between Seismic Source Detail and Crustal Heterogeneities in Spherical 3D Finite Element Modeling: the 2004 Sumatra Earthquake case-study

Finite Element methods (FEMs) are a powerful numerical simulation tool for modeling seismic events as they allow to solve three-dimensional complex models. We used a 3D Finite Element approach to evaluate the co-seismic displacement eld produced by the devastating 2004 Sumatra Andaman earthquake, which caused permanent deformations recorded by continuously operating GPS networks in a region of unprecedented extent. Previous analysis of the static displacement fi eld focused on the heterogeneous distribution of moment release on the fault plane; our intention here is to investigate how much the presence of crustal heterogeneities trades off seismic source details. To this aim, we adopted a quite simple source model in modeling the event. The key feature of our analysis is the generation of a complex three dimensional spherical domain. More-over, we also made an accurate analysis concerning boundary conditions, which are crucial for FE simulations