Analysis of the nonlinear storage–discharge relation for hillslopes through 2D numerical modelling

Storage–discharge curves are widely used in several hydrological applications concerning flow and solute transport in small catchments. This article analyzes the relation Q(S) (where Q is the discharge and S is the saturated storage in the hillslope), as a function of some simple structural parameters. The relation Q(S) is evaluated through two‐dimensional numerical simulations and makes use of dimensionless quantities. The method lies in between simple analytical approaches, like those based on the Boussinesq formulation, and more complex distributed models. After the numerical solution of the dimensionless Richards equation, simple analytical relations for Q(S) are determined in dimensionless form, as a function of a few relevant physical parameters. It was found that the storage–discharge curve can be well approximated by a power law function Q/(LK_s) = a(S/(L²(? ? θ_r)))^b, where L is the length of the hillslope, K_s the saturated conductivity, ? ? θ_r the effective porosity, and a, b two coefficients which mainly depend on the slope. The results confirm the validity of the widely used power law assumption for Q(S). Similar relations can be obtained by performing a standard recession curve analysis. Although simplified, the results obtained in the present work may serve as a preliminary tool for assessing the storage–discharge relation in hillslopes. Copyright © 2012 John Wiley & Sons, Ltd.     

Rapid estimation of ground-shaking maps for seismic emergency management in the Campania Region of southern Italy

Strong ground-shaking mapping soon after a moderate-to-large earthquake is crucial to recognize the areas that have suffered the largest damage and losses. These maps have a fundamental role for emergency services, loss estimation and planning of emergency actions by the Civil Protection Authorities. This is particularly important for areas with high seismic risk levels, such as the Campania-Lucania Region in southern Italy. Taking advantage of the Irpinia Seismic Network (ISNet), a recently installed dense and wide dynamic seismic network, we have developed a procedure for rapid estimation of ground-shaking maps after moderate-to-large earthquakes (GRSmap). This uses an optimal data gridding scheme designed to account for bi-dimensional features of strong ground-motion fields, such as directivity, radiation patterns and focal mechanisms, to which most damage can be correlated. The basis of the mapping technique is a triangulation procedure to locally correct predicted data at the triangle barycentres where their vertices correspond to seismic stations. The method has been tested off-line using a simulated M 6.6 earthquake located at the centre of ISNet and applied to data of the 23 November 1980 Irpina M 6.9 earthquake recorded by a sparse network. This has highlighted its ability to predict peak ground-motion parameters of large magnitude earthquakes with respect to the attenuation relationships.     

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.     

Environmental history since 11,000 C yr B.P. of the northeastern Pampas, Argentina, from alluvial sequences of the Luján River

Sedimentological, malacological, and pollen analyses from ¹⁴C-dated alluvial sections from the Luján River provide a detailed record of environmental changes during the Holocene in the northeastern Pampas of Argentina. From 11,200 to 9000 ¹⁴C yr B.P., both sedimentary and biological components suggest that the depositional environment was eutrophic, alkaline, and freshwater to brackish shallow water bodies without significant water circulation. During this time, bioclastic sedimentation was dominant and the shallow water bodies reached maximum development as the climate became more humid, suggesting an increase in precipitation. Short-term fluctuations in climate during the last stage of this interval may have been sufficient to initiate changes in the water bodies, as reduction of the volume alternated with periods of flooding. The beginning of the evolution of shallow swamps in the wide floodplain or huge wetlands was contemporaneous with a sea level lower than the present one. From 9000 and 7000 ¹⁴C yr B.P., mesotrophic, alkaline, brackish, probably anoxic swamps existed. Between 7000 and 3000 ¹⁴C yr B.P., anoxic calcareous swamps were formed, with subaerial exposure and development of the Puesto Berrondo Soil (3500-2900 ¹⁴C yr B.P.). A trend to a reduction of water bodies is recorded from 9000 to ca. 3000 ¹⁴C yr B.P., with a significant reduction after ca. 7000 ¹⁴C yr B.P. A shift to subhumid-dry climate after 7000 ¹⁴C yr B.P. appears to be the main cause. During this time, an additional external forcing toward higher groundwater levels was caused by Holocene marine transgression causing changes in the water bodies levels. The climate became drier during the late Holocene (ca. 3000 yr B.P.), when clastic sedimentation increased, under subhumid-dry conditions. Flood events increased in frequency during this time. From ca. A.D. 1790 to present, the pollen record reflects widespread disturbance of the vegetation during the European settlement.     

High frequency strong ground motion modelling in the Catania area associated with the Ibleo-Maltese fault system

A hybrid statistical-deterministic approach has been applied to estimate strong ground motion parameters (PGA, spectral ordinates) in South-Eastern Sicily for a M = 7 earthquake. A number of 100 different rupture processes have been simulated along a composite fault system representing two segments of the Ibleo-Maltese fault scarp. Map at regional scale of mean PGA in the 0.5–20 Hz frequency band shows highest values (0.4–0.5 g) nearby and North of Catania, due to a dominant directivity effect. The COV parameter, which expresses the variability of PGA values as a function of source complexity, is higher in the region nearby and South of the town of Augusta, where, depending on the rupture history, rather large PGA values can be observed (>0.4 g). PGA attenuation curves suggest that an azimuthal variation could be related to the source extent and directivity. The response and pseudo acceleration spectra are computed for different sites in the town of Catania including an approximate 1D site response. Ground motion amplification effects at high frequency (5–20 Hz) are produced by thin shallow layer of soft clay, loose pyroclastites and fill. We observe small amplification effects, in the frequency ranges 2–3 Hz and 5–10 Hz, in sites where recent alluvia reach a thickness of some tens of meters. Otherwise, sites located on outcrops of massive lavas show moderate attenuation.     

Statistical analysis of earthquake activity at Etna volcano (March 1981 eruption)

Seismic activity that preceded, accompanied, and followed the 17–23 March 1981 Etnean eruption has been statistically analyzed.On the grounds of both time evolution of seismicity and catalogue completeness, three time intervals have been defined (12 February–2 March, 12–17 March, 19–31 March) and for each of these periods both the b coefficient of theGutenberg-Richter's (1956) relationship and the E parameter of the cluster size (Shlien andToksoz, 1970) have been calculated.No significant variations were observed between the first and third periods, while lower values of bothb andE coefficients were found in the second one. These findings might indicate that changes in the seismicity features occur just before the eruption start.Small but fast variations in the stress field acting on the volcano might originate this type of seismic activity, while the importance of the tectonic control on volcanic phenomena seems to be confirmed.