Water-table and discharge changes associated with the 2016–2017 seismic sequence in central Italy: hydrogeological data and a conceptual model for fractured carbonate aquifers

A seismic sequence in central Italy from August 2016 to January 2017 affected groundwater dynamics in fractured carbonate aquifers. Changes in spring discharge, water-table position, and streamflow were recorded for several months following nine Mw 5.0–6.5 seismic events. Data from 22 measurement sites, located within 100 km of the epicentral zones, were analyzed. The intensity of the induced changes were correlated with seismic magnitude and distance to epicenters. The additional post-seismic discharge from rivers and springs was found to be higher than 9 m³/s, totaling more than 0.1 km³ of groundwater release over 6 months. This huge and unexpected contribution increased streamflow in narrow mountainous valleys to previously unmeasured peak values. Analogously to the L’Aquila 2009 post-earthquake phenomenon, these hydrogeological changes might reflect an increase of bulk hydraulic conductivity at the aquifer scale, which would increase hydraulic heads in the discharge zones and lower them in some recharge areas. The observed changes may also be partly due to other mechanisms, such as shaking and/or squeezing effects related to intense subsidence in the core of the affected area, where effects had maximum extent, or breaching of hydraulic barriers.

     

Large-Eddy Simulation of Atmospheric Boundary-Layer Flow Over Fluvial-Like Landscapes Using a Dynamic Roughness Model

A recently developed dynamic surface roughness model (Anderson and Meneveau, J Fluid Mech 679:288–314, 2011) for large-eddy simulation (LES) of atmospheric boundary-layer flow over multi-scale topographies is applied to boundary-layer flow over several types of fluvial-like landscapes. The landscapes are generated numerically with simulation of a modified Kardar–Parisi–Zhang equation (Passalacqua et al., Water Resour Res 42:WOD611, 2006). These surfaces possess the fractal-like channel network and anisotropic features often found in real terrains. The dynamic model is shown to lead to accurate flow predictions when the surface-height distributions exhibit power-law scaling (scale invariance) in the prevalent mean flow direction. In those cases, the LES provide accurate predictions (invariant to resolution) of mean velocity profiles. Conversely, some resolution dependence is found for applications in which the landscape’s streamwise spectra do not exhibit pure power-law scaling near wavenumbers corresponding to the LES grid resolution.     

Resistive instabilities in coronal conditions

Resistive instabilities in a context referring to the solar corona are rigorously investigated. Various equilibrium configurations are considered, differing, among other things, by their behaviour with respect to fast, ideal instabilities. The computations presented cover in a unified scheme all known regimes of resistive modes and allow one to determine the fastest timescale over which resistivity can play a role. Comparisons with previous work as well as possible extensions are presented.     

Magma contamination by direct wall rock interaction: constraints from xenoliths from the walls of a carbonate-hosted magma chamber (Vesuvius 1944 eruption)

During the 1944 eruption of Vesuvius different types of xenoliths were ejected. They represent fragments of the walls of a low volume (<0.5 km³) shallow (3–4 km depth) magma chamber. The study of these xenoliths enables us to estimate the amount of contamination occurring at the boundary of a high-T alkaline magma chamber hosted in carbonate rocks. The process of contamination of the magma by carbonates can be modelled, using isotopic and chemical data, as a mixing between magma and marbles. Mass exchanges occur at the boundary between the crystallizing magma and marble wall rocks, where endoskarn forms. The contamination of the solidification front of the chamber is very limited. The solidification front and the skarn shell effectively isolate the interior of the magma chamber from new inputs of contaminants from the carbonate wall rocks. Therefore, the main volume of magma, hosted in the magma chamber, did not undergo any significant mass exchange with the wall rocks.     

New constraints to the geodynamic evolution of the southern sector of the Calabria–Peloritani Arc (Italy)Nouvelles données sur l'évolution géodynamique du secteur méridional de l'Arc Calabro-Péloritain (Italie)

New studies have been carried out on the Tertiary of the Stilo Unit, the uppermost of the Calabria–Peloritani Arc southern sector, and the Stilo–Capo d'Orlando Formation, sealing the whole nappe stack. The Tertiary terrains linked to the Mesozoic cover of the Stilo Unit consist of the lowermost Oligocene Palizzi Formation and the Late Rupelian–Aquitanian Pignolo Formation. The possibility that they deposited before the emplacement of this unit as the highest tectonic sheet of the sector is suggested. The base of the Stilo–Capo d'Orlando Formation resulted of Burdigalian age in both type areas. This interpretation, together with the existing and new data, allows proposing an age close to the Aquitanian–Burdigalian boundary for the stacking of the whole Calabria–Peloritani Arc southern sector. To cite this article: G. Bonardi et al., C. R. Geoscience 334 (2002) 423–430.     

Thermal behaviour of davyne-group minerals

The thermal behaviour of microsommite (MC), davyne from Vesuvius (DV) and from Zabargad (DZ) was determined by X-ray single crystal data obtained employing a microfurnace connected to a four-circle diffractometer. Upon heating, the a parameter increased linearly, with similar thermal expansion rates for the three samples: the mean linear expansion coefficients, _a , were 10.2(3)·10^-6, 13.4(7)·10^-6, 15.1(8)·10^-1 K^-1 for MC, DV and DZ respectively.At about 473 K both MC and DZ showed a discontinuity in the expansion of the c parameter. The mean linear expansion coefficient, _c , changed abruptly from 16(4)·10^-6 K^-1 for both minerals below the discontinuity to 2(1)·10^-6 and 3(1)·10^-6 K^-1 for MC and DZ, respectively, above the discontinuity. In DV, however, the _c coefficient was constant between 293 und 827 K and equal to 1(2)·10^-6 K^-1.     

Complexity in a cellular model of river avulsion

We propose a new model of river avulsion that emphasizes simplicity, self-organization, and unprogrammed behavior rather than detailed simulation. The model runs on a fixed cellular grid and tracks two elevations in each cell, a high elevation representing the channel (levee) top and a low one representing the channel bottom. The channel aggrades in place until a superelevation threshold for avulsion is met. After an avulsion is triggered a new flow path is selected by steepest descent based on the low values of elevation. Flow path depends sensitively on floodplain topography, particularly the presence of former abandoned channels. Several behavioral characteristics emerge consistently from this simple model: (1) a tendency of the active flow to switch among a small number of channel paths, which we term the active channel set, over extended periods, leading to clustering and formation of multistory sand bodies in the resulting deposits; (2) a tendency for avulsed channels to return to their previous paths, so that new channel length tends to be generated in relatively short segments; and (3) avulsion-related sediment storage and release, leading to pulsed sediment output even for constant input. Each of these behaviors is consistent with observations from depositional river systems. A single-valued threshold produces a wide variety of avulsion sizes and styles. Larger “nodal” avulsions are rarer because pre-existing floodplain topography acts to steer flow back to the active channel. Channel stacking pattern is very sensitive to floodplain deposition. This work highlights the need to develop models of floodplain evolution at large time and space scales to complement the improving models of river channel evolution.     

Ground-Motion Scaling in the Western Alps

In order to empirically obtain the scaling relationships for the high-frequency ground motion in the Western Alps (NW Italy), regressions are carried out on more than 7500 seismograms from 957 regional earthquakes. The waveforms were selected from the database of 6 three-component stations of the RSNI (Regional Seismic network of Northwestern Italy). The events, M _W ranging between 1.2 and 4.8, were recorded within a hypocentral distance of 200 km during the time period: 1996–2001. The peak ground velocities are measured in selected narrow-frequency bands, between 0.5 and 14 Hz. Results are presented in terms of a regional attenuation function for the vertical ground motion, a set of vertical excitation terms at the reference station STV2 (hard-rock), and a set of site terms (vertical and horizontal), all relative to the vertical component of station STV2.The regional propagation of the ground motion is modeled after quantifying the expected duration of the seismic motion as a function of frequency and hypocentral distance. A simple functional form is used to take into account both the geometrical and the anelastic attenuation: a multi-variable grid search yielded a quality factor Q(f) = 310f ^0.20, together with a quadri-linear geometrical spreading at low frequency. A simpler, bi-linear geometrical spreading seems to be more appropriate at higher frequencies (f > 1.0 Hz). Excitation terms are matched by using a Brune spectral model with variable, magnitude-dependent stress drop: at M _w 4.8, we used Δσ = 50 MPa. A regional distance-independent attenuation parameter is obtained (κ₀ = 0.012 s) by modelling the average spectral decay at high frequency of small earthquakes.In order to predict the absolute levels of ground shaking in the region, the excitation/attenuation model is used through the Random Vibration Theory (RVT) with a stochastic point-source model. The expected peak-ground accelerations (PGA) are compared with the ones derived by Ambraseys et al. (1996) for the Mediterranean region and by Sabetta and Pugliese (1996) for the Italian territory.