Archaeoseismology: Methodological issues and procedure

Archaeoseismic research contributes important data on past earthquakes. A limitation of the usefulness of archaeoseismology is due to the lack of continuous discussion about the methodology. The methodological issues are particularly important because archaeoseismological investigations of past earthquakes make use of a large variety of methods. Typical in situ investigations include: (1) reconstruction of the local archaeological stratigraphy aimed at defining the correct position and chronology of a destruction layer, presumably related to an earthquake; (2) analysis of the deformations potentially due to seismic shaking or secondary earthquake effects, detectable on walls; (3) analysis of the depositional characteristics of the collapsed material; (4) investigations of the local geology and geomorphology to define possible natural cause(s) of the destruction; (5) investigations of the local factors affecting the ground motion amplifications; and (6) estimation of the dynamic excitation, which affected the site under investigation. Subsequently, a ‘territorial’ approach testing evidence of synchronous destruction in a certain region may delineate the extent of the area struck by the earthquake. The most reliable results of an archaeoseismological investigation are obtained by application of modern geoarchaeological practice (archaeological stratigraphy plus geological–geomorphological data), with the addition of a geophysical-engineering quantitative approach and (if available) historical information. This gives a basic dataset necessary to perform quantitative analyses which, in turn, corroborate the archaeoseismic hypothesis. Since archaeoseismological investigations can reveal the possible natural causes of destruction at a site, they contribute to the wider field of environmental archaeology, that seeks to define the history of the relationship between humans and the environment. Finally, through the improvement of the knowledge on the past seismicity, these studies can contribute to the regional estimation of seismic hazard.     

Benthic foraminifera in sandy (siliciclastic) coastal sediments of the Arabian Gulf (Saudi Arabia): a technical report

The Uromia–Dokhtar Magmatic Arc (UDMA) is a northwest–southeast trending magmatic belt which is formed due to oblique subduction of Neotethys underneath Central Iran and dominantly comprises magmatic rocks. The Jebal-e-Barez Plutonic Complex (JBPC) is located southeast of the UDMA and composed of quartz diorite, granodiorite, granite, and alkali granite. Magmatic enclaves, ranging in composition from felsic to mafic, are abundant in the studied rocks. Based on the whole rock and mineral chemistry study, the granitoids are typically medium-high K calc-alkaline and metaluminous to peraluminous that show characteristics of I-type granitoids. The high field strength (HFS) and large ionic radius lithophile (LIL) element geochemistry suggests fractional crystallization as a major process in the evolution of the JBPC. The tectonomagmatic setting of the granitoids is compatible with the arc-related granitic suite, a pre-plate collision granitic suite, and a syncollision granitic suite. Field observations and petrographic and geochemical studies suggest that the rocks in this area are I-type granitoids and continental collision granitoids (CCG), continental arc granitoids (CAG), and island arc granitoid (IAG) subsections. The geothermobarometry based on the electron probe microanalysis of amphibole, feldspars, and biotite from selected rocks of JBPC implies that the complex formed at high-level depths (i.e., 9–12 km; upper continental crust) and at temperatures ranging from 650 to 750 °C under oxidation conditions. It seems that JBPC is located within a shear zone period, and structural setting of JBPC is extensional shear fractures which are product of transpression tectonic regime. All available data suggested that these granitoids may be derived from a magmatic arc that was formed by northeastern ward subduction of the Neotethyan oceanic crust beneath the Central Iran in Paleogene and subsequent collision between the Arabian and Iranian plates in Miocene.

     

Syn-thrust deformation across a transverse zone: the Grem–Vedra fault system (central Southern Alps, N Italy)

The lateral continuity of the E?CW trending thrust sheets developed within the Lower to Middle Triassic cover of the central Southern Alps (Orobic belt) is disturbed by the occurrence of several N?CS trending transverse zones, such as the poorly known Grem?CVedra Transverse Zone (GVTZ). The GVTZ developed during the emplacement of the up to six S-verging thrust sheets consisting of Lower to Middle Triassic units, occurring immediately south of the Orobic Anticlines. The transverse zone, active during thrust emplacement related to the early Alpine compressions which pre-date the Adamello intrusion, includes three major vertical shear zones, the Grem, Pezzel and Zuccone faults. The major structure of the transverse zone is the dextral Grem fault, forming a deep lateral ramp between thrust sheets 3 and 5. A similar evolution also occurred along the Zuccone and Pezzel faults, which show a left-lateral displacement of syn-thrust folds. The Grem fault was later reactivated as an oblique tear fault during the emplacement of the Orobic Anticlines, due to back-thrusting along out-of-sequence thrust surfaces (Clusone fault). Transpressional deformations along the fault zone are recorded by the rotation of major syn-thrust folds, which also suggest a horizontal offset close to 0.5?km. Records of the first stage of evolution of the Grem fault are better preserved along its northern segment, and structural relationships suggest that it propagated southward and downward in the growing thrust stack. The study of the meso and megascopic structures developed along the GVTZ constrains the evolution of the transverse zone, illustrating the complex deformational phenomena occurring in a transpressional regime. The GVTZ probably reflects the existence of pre-existing tectonic lineaments with a similar orientation. Evidence of pre-existing structures are not preserved in the exposed units, nevertheless the N?CS extensional fault systems that characterize the Norian to Jurassic rifting history of the Lombardian basin are valid candidates.     

Sub-seismic scale fracture pattern and in situ permeability data in the chalk atop of the Krempe salt ridge at Lägerdorf, NW Germany: Inferences on synfolding stress field evolution and its impact on fracture connectivity

Chalk is exposed in the Heidestrasse quarry at Lägerdorf, at the top of the NE-SW trending Krempe salt ridge. Structural data indicate the presence of two joint sets, striking almost parallel and perpendicular to the salt ridge, respectively, and of a set of conjugate extensional faults and fault zones striking NW-SE, i.e. almost perpendicular to the salt ridge. Within the overall NW-SE trend of joints and faults, strike variations occur from the massive chalk exposed in the lower half of the quarry, to the overlying layered chalk. A large variability characterizes the normalized spacing of both joint sets, which does not show any clear trend with layer dip. In situ measurements indicate that the cross-sectional permeability of tight joints increases 1-2 orders of magnitude with respect to the undeformed chalk. We propose that joint and fault azimuthal variability resulted from changes through time of the stress ellipsoid at the top of the salt ridge, while joint spacing variability is associated with the weak mechanical influence of bedding in chalk. Azimuthal variability improves fracture connectivity and, hence, permeability and fluid flow.     

The effect of composition and cation ordering on the compressibility of columbites up to 7 GPa

The unit-cell parameters of two columbite samples along the (Fe,Mn)Nb₂O₆ solid solution were measured by means of high-pressure single-crystal X-ray diffraction up to pressures of 7 GPa. The compressional behaviour of these minerals was studied as a function of composition and degree of order. The P–V data of all the samples were fitted with a third-order Birch–Murnaghan equation of state. For the two samples with different compositions but identical degree of order the substitution of Mn for Fe causes a decrease of the bulk modulus K _T0, from 153(1) to 146(1) GPa, without any effect on the pressure first derivative K′. For the two samples with the same composition, cation ordering causes an increase of the bulk modulus from 149(1) to 153(1) GPa and of the pressure first derivative from 4.1(2) to 4.8(3). The compressional behaviour is anisotropic with a linear axial compressibility scheme β _b > β _c ≥ β _a for all samples, regardless of composition and degree of order. Such anisotropy increases sligthly with increasing Mn content.     

Estimation procedures for exchangeable Marshall copulas with hydrological application

Complex phenomena in environmental sciences can be conveniently represented by several inter-dependent random variables. In order to describe such situations, copula-based models have been studied during the last year. In this paper, we consider a novel family of bivariate copulas, called exchangeable Marshall copulas. Such copulas describe both positive and (upper) tail association between random variables. Specifically, inference procedures for the family of exchangeable Marshall copulas are introduced, based on the estimation of their (univariate) generator. Moreover, the performance of the proposed methodologies is shown in a simulation study. Finally, an illustration describes how the proposed procedures can be useful in a hydrological application.     

The Calabrian Stage redefined

The name Calabrian was introduced in the geological literature by the French stratigrapher Maurice Gignoux in 1910, and later described in his important monograph （633 pages） ＂Les formations marines pliocknes et quaternaires de l＇ltalie du sud et de la Sicile ＂published in 1913. Detailed data were provided on several sections （Santa Maria di Catanzaro, Caraffa, Monasterace, Palermo） and on their fossil content. The Calabrian Stage has commonly been used for over fifty years as the oldest subdivision of the Qua- ternary, notably in the time scales of Berggren ＆ van Cou- vering （1974） and Haq ＆ Eysinga （1987）. However, after the GSSP for the Pliocene/Pleistocene boundary （P/P） was approved by INQUA in 1982 and ratified by lUGS in 1984 at the Vrica section of Calabria, there was a decline in the usage of the stage name, and an increasing tendency by many Quaternary workers to question the boundary stratotype. This was because there was increasing evidence that it did not correspond to the beginning of the ＂ice age＂. In doing so, they were not complying with the recommendations presented at the 18th International Geological Congress （IGC） in London, 1948 （Oakley, 1950）.     

Hierarchical simulation of aquifer heterogeneity: implications of different simulation settings on solute-transport modeling

The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.