A GIS-based approach for gully erosion susceptibility modelling: a test in Sicily, Italy

The aim of this study is to analyze the susceptibility conditions to gully erosion phenomena in the Magazzolo River basin and to test a method that allows for driving the factors selection. The study area is one of the largest (225 km²) watershed of southern Sicily and it is mostly characterized by gentle slopes carved into clayey and evaporitic sediments, except for the northern sector where carbonatic rocks give rise to steep slopes. In order to obtain a quantitative evaluation of gully erosion susceptibility, statistical relationships between the spatial distributions of gullies affecting the area and a set of twelve environmental variables were analyzed. Stereoscopic analysis of aerial photographs dated 2000, and field surveys carried out in 2006, allowed us to map about a thousand landforms produced by linear water erosion processes, classifiable as ephemeral and permanent gullies. The linear density of the gullies, computed on each of the factors classes, was assumed as the function expressing the susceptibility level of the latter. A 40-m digital elevation model (DEM) prepared from 1:10,000-scale topographic maps was used to compute the values of nine topographic attributes (primary: slope, aspect, plan curvature, profile curvature, general curvature, tangential curvature; secondary: stream power index; topographic wetness index; LS-USLE factor); from available thematic maps and field checks three other physical attributes (lithology, soil texture, land use) were derived. For each of these variables, a 40-m grid layer was generated, reclassifying the topographic variables according to their standard deviation values. In order to evaluate the controlling role of the selected predictive variables, one-variable susceptibility models, based on the spatial relationships between each single factor and gullies, were produced and submitted to a validation procedure. The latter was carried out by evaluating the predictive performance of models trained on one half of the landform archive and tested on the other. Large differences of accuracy were verified by computing geometric indexes of the validation curves (prediction and success rate curves; ROC curves) drawn for each one-variable model; in particular, soil texture, general curvature and aspect demonstrated a weak or a null influence on the spatial distribution of gullies within the studied area, while, on the contrary, tangential curvature, stream power index and plan curvature showed high predictive skills. Hence, predictive models were produced on a multi-variable basis, by variously combining the one-variable models. The validation of the multi-variables models, which generally indicated quite satisfactory results, were used as a sensitivity analysis tool to evaluate differences in the prediction results produced by changing the set of combined physical attributes. The sensitivity analysis pointed out that by increasing the number of combined environmental variables, an improvement of the susceptibility assessment is produced; this is true with the exception of adding to the multi-variables models a variable, as slope aspect, not correlated to the target variable. The addition of this attribute produces effects on the validation curves that are not distinguishable from noise and, as a consequence, the slope aspect was excluded from the final multi-variables model used to draw the gully erosion susceptibility map of the Magazzolo River basin. In conclusion, the research showed that the validation of one-variable models can be used as a tool for selecting factors to be combined to prepare the best performing multi-variables gully erosion susceptibility model.     

On the evaluation of seismic response of structures by nonlinear static methods

In the most recent seismic codes, the assessment of the seismic response of structures may be carried out by comparing the displacement capacity, provided by nonlinear static analysis, with the displacement demand. In many cases the code approach is based on the N2 method proposed by Fajfar, which evaluates the displacement demand by defining, as an intermediate step, a single degree‐of‐freedom (SDOF) system equivalent to the examined structure. Other codes suggest simpler approaches, which do not require equivalent SDOF systems, but they give slightly different estimation of the seismic displacement demand. The paper points out the differences between the methods and suggests an operative approach that provides the same accuracy as the N2 method without requiring the evaluation of an equivalent SDOF system. A wide parametric investigation allows an accurate comparison of the different methods and demonstrates the effectiveness of the proposed operative approach. Copyright © 2009 John Wiley & Sons, Ltd.     

Effect of Overstrength on the Seismic Behaviour of Multi-Storey Regularly Asymmetric Buildings

In past years, seismic response of asymmetric structures has been frequently analysed by means of single-storey models, because of their simplicity and low computational cost. However, it is widely believed that use of more realistic multi-storey models is needed in order to investigate effects of some system characteristics (such as overstrength, higher modes of vibration, etc.) that make behaviour of multi-storey schemes different from that of single-storey systems. This paper examines effects of the overstrength in element cross-sections on the seismic behaviour of multi-storey asymmetric buildings. It is shown that in actual buildings this characteristic, which is sometimes very variable both in plan and along the height of the building, may lead to distributions of ductility demands different from those expected according to the results from single-storey models. Consequently, torsional provisions, which aim at reducing ductility demands of single-storey asymmetric systems to those of the corresponding torsionally balanced systems, should be re-checked in light of the behaviour of realistic multi-storey buildings.     

Seismic retrofitting of braced frame buildings by RC rocking walls and viscous dampers

A design procedure for seismic retrofitting of concentrically and eccentrically braced frame buildings is proposed and validated in this paper. Rocking walls are added to the existing system to ensure an almost uniform distribution of the interstorey displacement in elevation. To achieve direct and efficient control over the seismic performance, the design procedure is founded on the displacement‐based approach and makes use of overdamped elastic response spectra. The top displacement capacity of the building is evaluated based on a rigid lateral deformed configuration of the structure and on the ductility capacity of the dissipative members of the braced frames. The equivalent viscous damping ratio of the braced structure with rocking walls is calculated based on semi‐empirical relationships specifically calibrated in this paper for concentrically and eccentrically braced frames. If the equivalent viscous damping ratio of the structure is lower than the required equivalent viscous damping ratio, viscous dampers are added and arranged between the rocking walls and adjacent reaction columns. The design internal forces of the rocking walls are evaluated considering the contributions of more than one mode of vibration. The proposed design procedure is applied to a large set of archetype braced frame buildings and its effectiveness verified by nonlinear dynamic analysis.     

Corrective eccentricities for assessment by the nonlinear static method of 3D structures subjected to bidirectional ground motions

Nonlinear static methods are reliable in the evaluation of the seismic response of planar structural schemes, but they are not very effective in the assessment of three-dimensional building structures. The authors of this paper have recently proposed a nonlinear static approach for asymmetric structures, which is an improvement on that stipulated by seismic codes. This method is based on the observation that the distribution of the maximum dynamic displacements of the deck can be enveloped by two pushover analyses performed by applying the lateral force with two eccentricities with respect to the center of mass of the deck. These eccentricities, named “corrective eccentricities”, are defined so that the two corresponding pushover analyses provide displacements that are equal to those evaluated by nonlinear dynamic analysis at the two sides of the deck. In this paper, the corrective eccentricities are determined for a wide set of single-story systems. The equations for their analytical evaluation are determined and their reliability is demonstrated. Finally, the analysis of a multi-story structure is conducted to show how the method can be applied to real buildings. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison of nonlinear static methods for the assessment of asymmetric buildings

Traditional nonlinear static methods, e.g. the original version of the N2 method implemented in Eurocode 8, are not always effective in the assessment of asymmetric structures because of the errors committed in the evaluation of the torsional response. To overcome this shortcoming, two methods have recently been suggested by Kreslin and Fajfar (Bull Earthquake Eng 10(2):695–715, 2012) and Bosco et al. (Earthq Eng Struct Dyn 41:1751–1773, 2012). In particular, the method proposed by Kreslin and Fajfar adjusts the results of the nonlinear static analysis by means of those of a standard modal response spectrum analysis. In the method proposed by Bosco et al., the researchers suggested the use of two nonlinear static analyses characterized by lateral forces applied to different points of the deck. In this paper, the two improved nonlinear static methods and the original N2 method are applied to predict the maximum dynamic response of single- and multi-storey systems subjected to artificial and recorded accelerograms. The results highlight that the improved nonlinear static methods provide estimates which are more accurate than those of the original N2 method. Further, the comparison of the results identifies the range of the structural properties within which the original N2 method is still reliable and the range within which one of the two improved methods should be preferred.     

Modelling the effects of soil type and root distribution on shallow groundwater resources

Vegetated, shallow groundwater environments typically have high environmental and economic value. A sound understanding of the complex interactions and feedbacks between surface vegetation and groundwater resources is crucial to managing and maintaining healthy ecosystems while responding to human needs. A vegetated shallow groundwater environment was modelled using the software HYDRUS 2D to investigate the effects of several combinations of soil type and root distributions on shallow groundwater resources. Three rainfall regimes coupled to both natural and anthropogenically affected groundwater conditions were used to investigate the effect that combinations of four soil types and five root distributions can have on (a) groundwater level drops, (b) groundwater depletion, (c) groundwater recharge and (d) water stress conditions. Vegetation with roots distributed across the whole unsaturated zone and vegetation with dimorphic root systems (i.e. roots having larger concentrations both near the surface and the capillary fringe) behaved differently from vegetation growing roots mainly near the saturated zone. Specifically, vegetation with roots in the unsaturated zone caused water‐table drops and groundwater depletions that were half the amount due to deep‐rooted vegetation. Vegetation with a large portion of roots near the soil surface benefited from rainfall and was less vulnerable to water‐table lowering; as such, the fraction of the total area of roots affected by water stress conditions could be 40% smaller than in the case with deep‐rooted vegetation. However, roots uniformly distributed in the unsaturated zone could halve groundwater recharge rates observed in bare soils. Our analysis provided insights that can enable the formulation of site‐ and purpose‐specific management plans to respond to both human and ecosystem water requirements. Copyright © 2015 John Wiley & Sons, Ltd.     

Weakening of gneiss surfaces colonized by endolithic lichens in the temperate climate area of northwest Italy

A role of lithobionts in geomorphological processes is increasingly argued, but the spatio‐temporal scale of their impact is largely unexplored in many ecosystems. This study first characterizes in the temperate zone (northwest Italy) the relationships between lithobiontic communities including endolithic lichens and the hardness of their siliceous rock substrate (Villarfocchiardo Gneiss). The communities are characterized, on humid and xeric quarry surfaces exposed for decades and natural outcrops exposed for centuries, in terms of lichen and microbial constituents, using a combined morphological and molecular approach, and with regard to their development on and within the gneiss. A lichen species belonging to Acarosporaceae (Polysporina‐Sarcogyne‐Acarospora group, needing taxonomic revision) chasmoendolithically colonizes both the humid and xeric quarry surfaces, on which epilithic cyanobacterial biofilms and epilithic pioneer lichens respectively occur. Light and electron microscopic observations show the development of the endolithic thalli within rock microcracks and the hyphal penetration along crystal boundaries down to depths of 1 to 3 mm, more pronounced within the humid surfaces. Such colonization patterns are likely related to biogeophysical deterioration, while no chemical alteration characterizes minerals contacted by the endolithic lichen. By contrast, on natural outcrops, where the endolithic colonization is negligible, a reddish rind below epilithic lichens indicates chemical weathering processes. Schmidt Hammer measurements highlight that the endolithic lichens deeply affect the hardness of the gneiss (down to ?60% with respect to fresh controls and surfaces only colonized by cyanobacteria), exerting a significantly higher weakening effect with respect to the associated epilithic lithobionts. The phenomenon is more remarkable on humid than on xeric quarry surfaces and natural outcrops, where epilithic lichens are likely involved in long‐term hardening processes supporting surface stabilization. Endolithic lichens are thus active biogeomorphological agents at the upper millimetric layer of siliceous rocks in temperate areas, exerting their weakening action during the early decade‐scaled stages of surface exposure. Copyright © 2015 John Wiley & Sons, Ltd.     

Source location of long period seismicity at Volcàn de Colima,México

This paper presents an analysis of seismicity associated with the volcanic activity of Volcàn de Colima (México) and recorded in the period November 2005–April 2006 during a field survey by the Istituto Nazionale di Geofisica e Vulcanologia (INGV)–Osservatorio Vesuviano, the Observatorio Vulcanologico de Colima of Colima University and the Instituto Andaluz de Geofisica, University of Granada. Three different types of volcanic earthquakes have been identified on the basis of their spectral properties: Type A (0.3–1 Hz), Type B (1–5 Hz) and Type C (3–4 Hz). Results of polarization analysis applied to Type A events show a predominance of radial motion, indicating that the wavefield comprises compressional waves (P) and shear waves polarized in the vertical plane (SV), while the signal always begins with a negative polarity. Type A, B and C earthquakes have been located using both a flat layered model and a 3D model including topography. Hypocentre distributions indicate that the source of Type A signals is very shallow and confined to a small volume lying about 1 km below the crater. In contrast, the source of Type B and C events is significantly deeper, with most hypocentres located in a volume of about 1 km³ centred at 2.5–3 km depth. A cluster analysis based on the cross-correlation among the waveforms of different events recorded at the same station was applied to Type A earthquakes. Only two clusters, which include only a small percentage of events were found, indicating that earthquake families were uncommon during the period of our survey.