Chemical evolution of thermal waters and changes in the hydrothermal system of Papandayan volcano (West Java,Indonesia) after the November 2002 eruption

Papandayan is a stratovolcano situated in West Java, Indonesia. Since the last magmatic eruption in 1772, only few hydrothermal explosions have occurred. An explosive eruption occurred in November 2002 and ejected ash and altered rocks. The altered rocks show that an advanced argillic alteration took place in the hydrothermal system by interaction between acid fluids and rocks. Four zones of alteration have been defined and are limited in extension and shape along faults or across permeable structures at different levels beneath the active crater of the volcano.     

Landscape response to tectonic and climatic forcing in the foredeep of the southern Apennines,Italy: insights from Quaternary stratigraphy,quantitative geomorphic analysis,and denudation rate proxies

We present new data about the morphological and stratigraphic evolution and the rates of fluvial denudation of the Tavoliere di Puglia plain, a low‐relief landscape representing the northernmost sector of the Pliocene‐Pleistocene foredeep of the southern Apennines. The study area is located between the easternmost part of the southern Apennine chain and the Gargano promontory and it is characterized by several orders of terraced fluvial deposits, disconformably overlying lower Pleistocene marine clay and organized in a staircase geometry, which recorded the emersion and the long‐term incision history of this sector since mid‐Pleistocene times. We used the spatial and altimetric distribution of several orders of middle to late Pleistocene fluvial terraces in order to perform paleotopographic reconstruction and GIS‐aided eroded volumes estimates. Then, we estimated denudation rates on the basis of the terraces chronostratigraphy, supported by published OSL and AAR dating. Middle to upper Pleistocene denudation rates estimated by means of such an approach are slightly lower than 0.1 mm yr^‐1, in good agreement with short‐term data from direct and indirect evaluation of suspended sediment yield. The analysis of longitudinal river profiles using the stream power erosion model provided additional information on the incision rates of the studied area. Middle to late Quaternary uplift rates (about 0.15 mm yr^‐1), calculated on the basis of the elevation above sea level of marine deposits outcropping in the easternmost sector of the study area, are quite similar to the erosion rates average value, thus suggesting a steady‐state fluvial incision. The approach adopted in this work has demonstrated that erosion rates traditionally obtained by quantitative geomorphic analysis and k_sn estimations can be successfully integrated to quantify rates of tectonic or geomorphological processes of a landscape approaching steady‐state equilibrium. Copyright © 2014 John Wiley & Sons, Ltd.     

The role of the seismic trigger in the Calitri landslide (Italy): historical reconstruction and dynamic analysis

The study covered by this paper was focused on the historical case of the Calitri landslide, which was repeatedly reactivated by earthquakes, as reported since 1694. The town of Calitri (Southern Italy) is located on a ridge whose southern slope, from its top to the Ofanto river valley floor, has been historically affected by major landsliding. The last record of recurrence of the Calitri landslide leads back to the 1980 Irpinia earthquake, which caused significant damage to the town and had pervasive and visible ground effects. Based on a detailed historical reconstruction of landsliding and seismicity at Calitri, the study analysed the current static and dynamic stability of the landsliding slope by means of a finite-difference numerical analysis, taking into account the various factors of landslide initiation and kinematics. The results of the numerical analysis: (i) were consistent with the roto-translational mechanism observed upon the latest reactivation of the landslide; (ii) demonstrated that excess pore pressure redistribution caused a lag between the seismic trigger and the initiation of landsliding; and (iii) showed the impact of seismic input frequency on propagation and depth of slope instability.     

Effect of frictional heating and thermal advection on pre-seismic sliding: a numerical simulation using a rate-, state- and temperature-dependent friction law

Laboratory experiments on simulated faults in rocks clearly show the temperature dependence of dynamic rock friction. Since rocks surrounding faults are permeable, we have developed a numerical method to describe the thermo-mechanical evolution of the pre-seismic sliding phase which takes into account both the rate-, state- and temperature-dependent friction law and the heat advection term in the energy equation. We consider a laminar fluid motion perpendicular to a vertical fault plane and assume that fluids move away from the fault plane. A semi-analytical temperature solution which accounts for the variability of slip velocity and stress on the fault has been found. This solution has been generalized to the case of a time varying fluid velocity and then was used to include the thermal pressurization effect. After discretizing the temperature solution, the evolution of the system is obtained by the solution of a system of first order differential equations which allows us to determine the evolution of slip, slip rate, friction coefficient, effective normal stress, temperature and fluid velocity. The numerical solutions are found using a Runge-Kutta method with an adaptative stepsize control in time. When the thermal pressurization effects can be neglected, the heat advection effect gives rise to a delay, with respect to the purely conductive case, of the earthquake occurrence time. This delay increases with increasing permeability H of the system. When the thermal pressurization effects are taken into account the situation is opposite, i.e. the onset of instability tends to precede that of the purely conductive case. The advance in the time of occurrence of instability increases with increasing coefficient of thermal pressurization. In the small permeability range (H ≤ 10^?18 m²), the seismic moment and nucleation length of the pre-seismic phase are significantly smaller than those predicted by the purely conductive model.     

Peak runoff contributing area as hydrological signature of the probability distribution of floods

For the analysis of hydrological extremes and particularly in flood prediction, deeper investigation is needed on the relative effects of different hydrological processes acting at the basin scale in different hydroclimatic areas of the world. In this framework, the theoretical derivation of flood distribution shows a great potential for development and knowledge advancement. In addition, another promising path of investigation is represented by the use of distributed hydrological models via simulation modelling (including Monte Carlo, discrete event and continuous simulation). In this paper results of a theoretically derived flood frequency distribution are analyzed and compared with the results of a simulation scheme that uses a distributed hydrological model (DREAM) in cascade with a rainfall generator (IRP). The numerical simulation allows the reproduction of a large number of extreme events and provides insight into the main control for flood generation mechanisms with particular emphasis to the peak runoff contributing areas, highlighting the relevance of soil texture and morphology in different climatic environments. The proposed methodology is applied here to the Agri and the Bradano basin, in Southern Italy.     

A performance-based adaptive methodology for the seismic evaluation of multi-span simply supported deck bridges

A performance-based adaptive methodology for the seismic assessment of highway bridges is proposed. The proposed methodology is based on an Inverse (I), Adaptive (A) application of the Capacity Spectrum Method (CSM), with the capacity curve of the bridge derived through a Displacement-based Adaptive Pushover (DAP) analysis. For this reason, the acronym IACSM is used to identify the proposed methodology. A number of Performance Levels (PLs), for which the seismic vulnerability and seismic risk of the bridge shall be evaluated, are identified. Each PL is associated to a number of Damage States (DSs) of the critical members of the bridge (piers, abutments, joints and bearing devices). The IACSM provides the earthquake intensity level (PGA) corresponding to the attainment of the selected DSs, using over-damped elastic response spectra as demand curves. The seismic vulnerability of the bridge is described by means of fragility curves, derived based on the PGA values associated to each DS. The seismic risk of the bridge is evaluated as convolution integral of the product between the fragility curves and the seismic hazard curve of the bridge site. In this paper, the key aspects and basic assumptions of the proposed methodology are presented first. The IACSM is then applied to nine existing simply supported deck bridges, characterized by different types of piers and bearing devices. Finally, the IACSM predictions are compared with the results of nonlinear response time-history analysis, carried out using a set of seven ground motions scaled to the expected PGA values.     

Mt. Melbourne volcano,antarctica: Evidence of seismicity related to volcanic activity

Characterization of the microseismic activity (M _L <2.0) has been performed at Mt. Melbourne since 1990. We recorded a group of low frequency events with common morphological characteristics, i.e., an emerging onset, an unclear second phase and a sharply dropping coda. Spectral analysis of events recorded at more than one station indicates that the seismogram characteristics and spectral content are largely due to source effects. A polarization filter applied to a set of three component data revealed a first phase made up ofP waves followed (after about 0.9–1.4 sec.) by a second phase probably composed ofSH-type waves. Particle motion analysis detected a seismic ray angle direction mainly between N70°E and N110°E and apparent angle of incidence between 35° and 48° for the first phase. The studied seismicity was localized in an area on the eastern slope of Mt. Melbourne Volcano which presents a surface temperature anomaly (Mazzarini andSalvini, 1994). We formulate two hypotheses for the type of earthquakes recorded: 1) long-period events involving active presence of magmatic fluids in the source processes; 2) or the result of fracturing processes (shear?) in a medium characterized by transition between brittle and plastic behaviors. In the latter hypothesis the superficial thermal anomaly may be a symptom of this behavior at depth and is confirmed by the lown values observed for the exponential fit in the codaQ analysis.     

Modeling of Rayleigh waves dispersion in the Sannio region (Southern Italy) from seismic active refraction data

Short period surface waves, recorded during a seismic refractionsurvey in the Sannio region (Southern Italy), have been modeled to infera shallow velocity model for the area. Based on the decrease of resolutionwith depth, due to the bias on group velocity estimates arising frominterference of the Rayleigh waves with higher modes, we carried out aprocedure of fitting, with synthetic seismograms, of selected filtered traceswith a gaussian filter, having a width at half height equal to 1 Hz and acentral frequency lying in the range [1,4] Hz. We estimated the likelihoodbetween synthetic and observed seismograms by measuring their semblance.In this way we were able to infer a more refined local velocity modelcharacterized by a high Vp and Vs vertical gradient in the sedimentarycover. Two ad hoc resolution studies, based on group velocity andamplitude data respectively, indicate that the local velocity model is a goodvelocity model also for the entire studied area. The increase in the numberof available data when using amplitude information allows us to make amore selective choice in the model parameter space (Vp and Vs of eachlayer) and to solve for the Vp/Vs ratio. The inferred Vp velocity in thehalf-space is equal to 2.8 km/s. This value is in excellent agreement withthat inferred by other authors (3 km/s) by modeling P-wave travel timevs. distance. The best-fit model furnish low Vp/Vs for the sedimentarycover so indicating a high degree of the sediment's compaction in thestudied area. The inferred shallow high-velocity gradient indicates thatthe shallow sedimentary layer in the area could trap and focus the energytraveling into it.