Ultramafic intrusion triggers hydrothermal explosions at Colle Fabbri (Spoleto, Umbria), Italy

The eroded Colle Fabbri volcano comprises intrusive and extrusive rocks which cover an area of about 10,000 m² in the Umbria region, Central Italy. The outcrop is located at the SW boundary of the Umbria Valley Graben (Umbria Region) on an N140 normal fault inside the Intramountain Ultra-Alkaline Province of carbonatite and melilitite rocks of central–southern Italy. A field survey of the outcrop allowed a reconstruction of igneous activity events of this unusual small-scale volcano. It is younger than 0.7 m.y. and rests on Lower Pleistocene conglomerates with intercalated clay beds. A palaeosoil marks the base of the volcanic sequence. It follows several metres of extrusive breccias composed of fragments of thermometamorphosed clay and travertine. Key features of these breccias are mud shells on blocks, plastic mud lumps, slumps, and mechanical injection from one layer into another, as well as plastic and vesiculated, micro-brecciated matrix. The breccias are cemented by a variety of silicate, sulphate, and sulphate–carbonate minerals deposited by intense hydrothermal circulation. The breccias are related to phreatic explosions triggered by the eruption of a superficial hydrothermal system. Up to 1 m thick sill and a dykelet swarm intrude the breccias. In the southern part of the outcrop there is a plug of melilitolitic composition which intrudes the breccias and deforms them upward. A variety of contact breccias is also scattered around the sill and the plug. In some places, hyaloclastite formed when melt invaded water pockets contained in the encasing rocks under hydrostatic pressure conditions. A thermal aureole, which moves up to 10 m away from the contact, is characterised by cordierite–trydimite association, thus indicating high-temperature (>>1000 °C) contact phenomena.     

The Effect of Velocity Inversions on H/V

We analyzed the phenomenology of microtremor H/V curves under inversions in the shear-wave velocity (V_s) profile in the subsoil. Under no V_s inversion the spectral signature of the H/V peaks is found to be ‘eye-shaped’ with the horizontal components higher than the vertical. Conversely, under negative velocity gradients, numerous of differences emerge. I) A H/V ratio below 1 is observed for a wide range of frequencies, due to the decrease of the horizontal components below the vertical one. II) In the presence of persistent H/V < 1, small bumps in the H/V ratio given by local minima in the vertical spectral component may represent the relics of the peaks indicating resonances and stratigraphic discontinuities. As a consequence, in the presence of velocity inversions the H/V > 2 SESAME (2004) criterion fails but a stratigraphic interpretation may still be possible. III) The H/V curves should always be interpreted together with the single component spectra. IV) Microtremor H/V measurements for stratigraphic/microzonation purposes on stiff artificial soils, (asphalt, concrete, cement, pavements) should always be avoided since the latter often produce velocity inversions. This may have consequences in the intermediate to high frequency domain ( > 1 Hz) also in the application of reference site methods, like H_site/H_bedrock, to microtremor. Theoretical modeling confirms these experimental findings.     

Hydraulic Tomography: Continuity and Discontinuity of High‐K and Low‐K Zones

Hydraulic tomography is an emerging field and modeling method that provides a continuous hydraulic conductivity (K) distribution for an investigated region. Characterization approaches that rely on interpolation between one‐dimensional (1D) profiles have limited ability to accurately identify high‐K channels, juxtapositions of lenses with high K contrast, and breaches in layers or channels between such profiles. However, locating these features is especially important for groundwater flow and transport modeling, and for design and operation of in situ remediation in complex hydrogeologic environments. We use transient hydraulic tomography to estimate 3D K in a volume of 15‐m diameter by 20‐m saturated thickness in a highly heterogeneous unconfined alluvial (clay to sand‐and‐gravel) aquifer with a K range of approximately seven orders of magnitude at an active industrial site in Assemini, Sardinia, Italy. A modified Levenberg‐Marquardt algorithm was used for geostatistical inversion to deal with the nonlinear nature of the highly heterogeneous system. The imaging results are validated with pumping tests not used in the tomographic inversion. These tests were conducted from three of five clusters of continuous multichannel tubing (CMTs) installed for observation in the tomographic testing. Locations of high‐K continuity and discontinuity, juxtaposition of very high‐K and very low‐K lenses, and low‐K “plugs” are evident in regions of the investigated volume where they likely would not have been identified with interpolation from 1D profiles at the positions of the pumping well and five CMT clusters. Quality assessment methods identified a suspect high‐K feature between the tested volume and a lateral boundary of the model.     

Bedload hysteresis in a glacier‐fed mountain river

Sediment transport during flood events often reveals hysteretic patterns because flow discharge can peak before (counterclockwise hysteresis) or after (clockwise hysteresis) the peak of bedload. Hysteresis in sediment transport has been used in the literature to infer the degree of sediment availability. Counterclockwise and clockwise hysteresis have been in fact interpreted as limited and unlimited sediment supply conditions, respectively. Hysteresis has been mainly explored for the case of suspended sediment transport, but it was rarely reported for bedload transport in mountain streams. This work focuses on the temporal variability of bedload transport in an alpine catchment (Saldur basin, 18.6 km², Italian Alps) where bedload transport was monitored by means of an acoustic pipe sensor which detects the acoustic vibrations induced by particles hitting a 0.5m‐long steel pipe. Runoff dynamics are dominated by snowmelt in late spring/early summer, mostly by glacier melt in late summer/early autumn, and by a combination of the snow and glacier melt in mid‐summer. The results indicate that hysteretic patterns during daily discharge fluctuations are predominantly clockwise during the snowmelt period, likely due to the ready availability of unpacked sediments within the channel or through bank erosion in the lower part of the basin. On the contrary, counterclockwise hysteresis tend to be more frequent during late glacier melting period, possibly due to the time lag needed for sediment provided by the glacial and peri‐glacial area to be transported to the monitoring section. However, intense rainfall events occurring during the glacier melt period generated predominantly clockwise hysteresis, thus indicating the activation of different sediment sources. These results indicate that runoff generation processes play a crucial role on sediment supply and temporal availability in mountain streams. Copyright © 2014 John Wiley & Sons, Ltd.     

A flume experiment on wood storage and remobilization in braided river systems

This work investigates wood dynamics in braided streams through physical modelling in a mobile bed laboratory flume, with the specific objective to characterize wood storage and turnover as a function of wood input rate and of wood element type. Three parallel channels (1.7 m wide, 10 m long) filled with uniform sand were used to reproduce braided networks with constant water discharge and sediment feeding. Wood dowels with and without simplified root wads were regularly added at the upstream end of each flume at different input rates, with a 1:2:3 ratio between the three flumes. Temporal evolution of wood deposition patterns and remobilization rates were monitored by a series of vertical images that permitted the recognition of individual logs. Results show that wood tends to disperse in generally small accumulations (< 5 logs), with higher spatial density on top of sediment bars, and is frequently remobilized due to the intense morphological changes. The amount of wood stored in the channel depends on log input rate through a non‐linear relationship, and input rates exceeding approximately 100 logs/hour determine a sharp change in wood dynamics, with higher storage volume and augmented formation of large jams (> 10 elements) that are less prone to remobilization. Presence of root wads seems to play a minor role in wood deposition, but it reduces the average travel distance of logs. Turnover rates of logs were similar in the three flumes, independently of wood input rate and largely resembling the turnover rate of exposed bars. For the simulated conditions, significant effects of wood on bed morphology were not observed, suggesting that interactions with fine sediments and living vegetation are crucial to form large, stable wood jams able to bring about relevant morphological changes. Copyright © 2014 John Wiley & Sons, Ltd.     

A relocatable ocean model in support of environmental emergencies

During the Costa Concordia emergency case, regional, subregional, and relocatable ocean models have been used together with the oil spill model, MEDSLIK-II, to provide ocean currents forecasts, possible oil spill scenarios, and drifters trajectories simulations. The models results together with the evaluation of their performances are presented in this paper. In particular, we focused this work on the implementation of the Interactive Relocatable Nested Ocean Model (IRENOM), based on the Harvard Ocean Prediction System (HOPS), for the Costa Concordia emergency and on its validation using drifters released in the area of the accident. It is shown that thanks to the capability of improving easily and quickly its configuration, the IRENOM results are of greater accuracy than the results achieved using regional or subregional model products. The model topography, and to the initialization procedures, and the horizontal resolution are the key model settings to be configured. Furthermore, the IRENOM currents and the MEDSLIK-II simulated trajectories showed to be sensitive to the spatial resolution of the meteorological fields used, providing higher prediction skills with higher resolution wind forcing.     

Site effects in the Aterno River Valley (L’Aquila,Italy): comparison between empirical and 2D numerical modelling starting from April 6th 2009 Mw 6.3 earthquake

In this paper, we focused our attention on a cross-section of the Aterno River Valley where a good quality geological and geophysical dataset allowed to reconstruct accurately the geometry and the Vs profiles along all the plane of the section. Its trace is deliberately aligned close to the strong motion stations that recorded the Mw 6.3 (April 6th 2009) L’Aquila earthquake. We analysed strong and weak motion data available at these latter stations as well as at one of the temporary stations installed during the Microzonation activities and located on outcropping bedrock, in proximity of the cross-section. We used the H/V technique to select a reliable reference site and once we found it, we applied the SSR technique to compute amplification functions in correspondence of two strong motion stations. In turn, for both sites we performed a site response numerical modelling with two different 2D codes and we compared simulated versus experimental transfer functions. We found that the cross-section is well constrained based on the very reasonable agreement between results of numerical modelling and earthquake data analysis. We pointed out also a strong amplification of the deposit at the centre of the valley due to the constructive interference of S and surface waves, not predictable by means of 1D numerical modelling. We also compared the H/V as well as the SSR obtained from strong motion data with the ones computed from weak motion finding evidences of non-linearity in soil behaviour.     

Simplified seismic soil classification: the Vfz matrix

Site effect assessment studies aim at predicting the effect of seismic shaking on structures by modeling the subsoil as an oscillator coupled to another oscillator representing the construction. The resulting amplification functions and response spectra depend on so many strong assumptions and parameters that, in the standard engineering practice, simplified seismic classifications appear preferable to complex modeling procedures which can only offer an illusory better accuracy. Since stratigraphic seismic amplification is not properly related to the absolute rigidity of subsoil but to impedance contrasts, the standard simplified approaches based on the ‘average’ rigidity of subsoil in the first few meters (e.g. Vs30) can hardly be effective. Here it is proposed a simplified soil classification approach that takes into account the basic Physics of seismic amplification and its parameters, i.e. the average shear wave velocity of the cover layer, the resonance frequency and the impedance contrast between the cover and the bedrock, which we summarize as VfZ. A possible classification approach is illustrated through a set of examples.