Geometry of the neotectonic stress field in southern Italy: Geological and seismological evidence

The neotectonic regime in southern Italy has been evaluated by making a comparison between all the available structural and seismological data. The area investigated can be subdivided into four distinct zones which are characterized by different stress regimes. In the Southern Apennines the tensile axis of the stress field is oriented approximately NE-SW while the maximum principal stress (σ₁) is subvertical. In Northern Calabria, the tensile axis is ESE-WNW and the σ₁ axis is almost vertical. In the Catanzaro trough both the tensile axis and the σ₁ axis are subhorizontal and act E-W and N-S, respectively. Finally, the Strait of Messina zone is characterized by a tensile axis oriented E-W and by σ₁ being subvertical.     

A review of the geodynamic constraints on the development and evolution of geothermal systems in the Central Andean Volcanic Zone (18–28°Lat.S)

ABSTRACT

Geothermal exploration in the Central Andean Volcanic Zone (CAVZ) focuses on heat capacities of known geothermal systems, yet the role of faults, veins, fractures and folds on the evolution and migration of fluids is far from complete. Here, we present a compilation of He and Sr isotope data and newly generated structural maps to examine if particular tectonic configurations are associated with fluids migrating from different crustal levels. Accordingly, we defined three tectono-geothermal environments (T1–T3) depicting specific structural arrangements and spatial relation with geothermal and volcanic manifestations. T1 is dominated by left-lateral strike-slip NW-striking faults, and geothermal and volcanic manifestations occur along the traces of these structures. T2 is dominated by N-striking thrust faults and parallel fault-propagated folds, cut and displaced by NW-striking faults. Here, geothermal manifestations occur at fault intersections and at fold hinges. T3 is defined by left-lateral/normal NW-striking faults, with geothermal and volcanic manifestations lying along fault traces. Each tectono-geothermal environment yields distinctive isotope ratios and geothermal reservoir temperatures. T1 shows high helium and low strontium ratios, and temperatures between 220° and 310°C. T3 shows low helium and high strontium ratios and temperatures between 260° and 320°C. T2 isotope ratios fall between T1 and T3, with lowest (130°-250°C) reservoir temperatures. We argue that these particular isotope signatures are due to a structural control on reservoir location and orientation. The orientation of the fracture mesh genetically associated with each tectono-geothermal environment is a first-order control on the migration pathway of fluids. Our model shows that T1 allows fluxing of deeper fluids, T2 promotes storage and favors longer residence times and T3 enhances subvertical fluid migration. Results here help to explain the evolution of geothermal systems in a wider context, including fault systems and Sr and He isotope variations, thus providing a framework for geothermal exploration in the CAVZ.     

Tectonic rotations during the Chile Ridge collision and obduction of the Taitao ophiolite (southern Chile)

Abstract The < 6 Ma young Taitao ophiolite, exposed at the westernmost promontory of the Taitao Peninsula, is located approximately 40 km southeast of the Chile triple junction and consists of a complete sequence of oceanic lithosphere. Systematic sampling for paleomagnetic study was performed to understand the complex obduction processes of the ophiolite onto the forearc of the South American Plate. Two representative demagnetization paths of remanent magnetization vectors were observed. One is characterized by stable univectorial demagnetization paths and was observed in volcaniclastic rocks and dyke complexes. Orientations of their remanent magnetization vectors indicate various degrees of counterclockwise rotations. The other is characterized by multivectorial demagnetization paths and was observed in the plutonic units (gabbros and ultramafic rocks). From these, two distinct stable remanent magnetization vectors were isolated; one has high coercivity and the other has low coercivity along the demagnetization paths with little influence of viscous magnetizations. This suggests that the complex deformation history involved at least two rotational events. The clockwise rotation, inferred from high coercivity remanent magnetization vectors, was attributed to a ridge collision event and the counterclockwise rotation, inferred from the low coercivity remanent magnetization vectors, was attributed to an accommodation phase into the South American forearc during obduction and final emplacement of the ophiolite. Folds developed during this period. Paleomagnetic restorations of the internal structures of the plutonic units and dyke complexes suggest that they probably originated in a mid‐oceanic ridge environment near a transform fault. The counterclockwise rotation of the plutonic and dyke complex units during the obduction generated tectonic gaps between these and the basement. The volcaniclastic rocks must have been deposited at nearly their present location, filling the tectonic gaps, as less effect of tectonic rotation was identified on these rocks.     

General and specific growth marks in pegmatite beryls

Growth marks of pegmatite beryls can be divided into general and specific in relation to their types and distributions. Dislocation absence and growth bands have been recognized as general growth marks of the environment while minor element concentrations and optical anomalies as specific to the local growth medium.     

On the mathematics of storage routing

Some mathematical aspects of storage routing are investigated analytically and through numerical simulation. The minimum size of the time interval required to eliminate computed values at the beginning of the outflow hydrograph falling below the initial steady flow rate is given for two routing schemes. The stability behavior of a refined Muskingum scheme is explored; it is found that values of the weighting coefficient up to 1.0 may be used with appropriately specified time intervals. The concept of “relative” storage is analyzed and it is shown that it is unnecessary when the correct initial condition of the storage is considered. The classical best straight-line fit procedure between storage and weighted discharge is investigated. In an example case, it is demonstrated that this method of system identification is not reliable. The possibility of using discharge-dependent parameters is discussed.     

Most similar neighbor imputation of forest attributes using metrics derived from combined airborne LIDAR and multispectral sensors

In the context of predicting forest attributes using a combination of airborne LIDAR and multispectral (MS) sensors, we suggest the inclusion of normalized difference vegetation index (NDVI) metrics along with the more traditional LIDAR height metrics. Here the data fusion method consists of back-projecting LIDAR returns onto original MS images, avoiding co-registration errors. The prediction method is based on non-parametric imputation (the most similar neighbor). Predictor selection and accuracy assessment include hypothesis tests and over-fitting prevention methods. Results show improvements when using combinations of LIDAR and MS compared to using either of them alone. The MS sensor has little explanatory capacity for forest variables dependent on tree height, already well determined from LIDAR alone. However, there is potential for variables dependent on tree diameters and their density. The combination of LIDAR and MS sensors can be very beneficial for predicting variables describing forests structural heterogeneity, which are best described from synergies between LIDAR heights and NDVI dispersion. Results demonstrate the potential of NDVI metrics to increase prediction accuracy of forest attributes. Their inclusion in the predictor dataset may, however, in a few cases be detrimental to accuracy, and therefore we recommend to carefully assess the possible advantages of data fusion on a case-by-case basis.     

Cross-twinning in a natural spinel from Sri Lanka

A modified cross-twinning growth mechanism is put forward to explain the anomalous morphology of a spinel multiple-twin from Sri Lanka, flattened crosswise the twin planes. Cross-twinning in spinel was found also in other specimens from Pegu (Myanmar), and the results were published in a previous paper. This particular type of twinning is derived from the combination of cyclic twinning with lamellar twinning, so that these samples may be thought of as partial fivelings (cubic cyclic {111} twins with five components sharing a common <110> pseudo-fivefold axis). In the present paper, the sample from Sri Lanka has been suitably cut with the aim of focusing the study on the cross-twinning region. The transformation matrices that link the orientation states of each couple of twin components have been determined by means of White Beam Synchrotron Radiation Topography. They showed that the specimen is made up of four twin components (A, B, C and D), with three twin planes: and They also showed that the cross-twinned individuals (B and D) actually are not twinned to each other, and that a simple crystallographic relationship holds between them. X-ray diffraction topography by conventional source allowed to image the crossing-region and to determine that the cross-twinned individuals are in contact through a semi-coherent boundary, with twinning dislocations contributing to relieve the coherency strains. Electron probe microanalyses with wave dispersive spectroscopy showed that the chemical composition is almost homogeneous, at least within the spatial resolution limit of this technique. The similar growth features observed in the spinel sample from Sri Lanka and in those from Myanmar are interpreted as growth marks, indicators of a similar origin: in both cases they are found in impure dolomitic marbles. In particular, the specimen from Sri Lanka results from the interaction of thermal and metasomatic effects due to contact metamorphism. An unusual stepped morphology of the face close to the twin boundary, possibly due to corrosion and re-growth processes acted preferentially at a re-entrant corner by metasomatic fluids, is interpreted as indicator of a metasomatic event that succeeded to the crystal growth, the latter occurred by thermal effect.

Digital Urban Morphometrics: Automatic Extraction and Assessment of Morphological Properties of Buildings

The aim of this article is to present a method to calculate the morphological properties of the built environment using LiDAR (light detection and ranging) data, geographic information systems (GIS) data and three‐dimensional (3D) models of cities as a source of information. A hybrid approach that takes into account different types of inputs and consequently evaluates the accuracy of each type of used data is presented. This work is intended to give a first response to the lack of a comprehensive and accurate procedure that uses LiDAR data in order to automatically derive precise morphological properties, such as volumes and surfaces (façades and roofs) of buildings. The method was tested on two case‐study areas in the Geneva region with different characteristics, one in the old town along the Rhone River and the other on the CERN campus. A statistical analysis that compares the results of the computation with the 3D model of the built environment was used to validate the results, complemented by significance statistical tests. Outcomes showed that the proposed method to derive morphological properties can reach high levels of accuracy, thus enhancing the potential uses of LiDAR data for numerous applications, typically for the assessment of the urban environmental quality (UEQ) at the city and district scale, such as the estimation of the potential deployment of renewable energies in the built environment and the determination and monitoring of several urban indicators.     

Seismic damage accumulation in multiple mainshock–aftershock sequences

Earthquakes are generally clustered, both in time and space. Conventionally, each cluster is made of foreshocks, the mainshock, and aftershocks. Seismic damage can possibly accumulate because of the effects of multiple earthquakes in one cluster and/or because the structure is unrepaired between different clusters. Typically, the performance-based earthquake engineering (PBEE) framework neglects seismic damage accumulation. This is because (i) probabilistic seismic hazard analysis (PSHA) only refers to mainshocks and (ii) classical fragility curves represent the failure probability in one event, of given intensity, only. However, for life cycle assessment, it can be necessary to account for the build-up of seismic losses because of damage in multiple events. It has been already demonstrated that a Markovian model (i.e., a Markov chain), accounting for damage accumulation in multiple mainshocks, can be calibrated by maintaining PSHA from the classical PBEE framework and replacing structural fragility with a set of state-dependent fragility curves. In fact, the Markov chain also works when damage accumulates in multiple aftershocks from a single mainshock of known magnitude and location, if aftershock PSHA replaces classical PSHA. Herein, this model is extended further, developing a Markovian model that accounts, at the same time, for damage accumulation: (i) within any mainshock–aftershock seismic sequence and (ii) among multiple sequences. The model is illustrated through applications to a series of six-story reinforced concrete moment-resisting frame buildings designed for three sites with different seismic hazard levels in Italy. The time-variant reliability assessment results are compared with the classical PBEE approach and the accumulation model that only considers mainshocks, so as to address the relevance of aftershocks for life cycle assessment.