Correlation between damage distribution and ambient noise H/V spectral ratio: the SESAME project results

In the framework of the SESAME project one of the tasks was the compilation of all available ambient noise measurements within urban environments affected by historical or/and recent strong earthquakes in Europe. The aim of such a task was to give an answer to the question; “How does horizontal-to-vertical ambient noise spectral ratio compare with damage in modern cities?”. For this purpose five European urban areas, namely, Angra do Heroismo (Portugal), Fabriano and Palermo (Italy), Thessaloniki and Kalamata (Greece) were selected for which spatial damage information was available either in terms of modified Mercalli intensity or in EMS98 damage grades. The geological setting of the examined sites as well as the causative earthquakes are satisfactorily known. Ambient noise recordings compiled for all examined sites have been homogeneously processed by a technique developed and agreed upon SESAME project. Using a standard multivariate statistical analysis, namely, factor analysis and canonical correlation, the horizontal-to-vertical ambient noise spectral ratio (HVNSR) is correlated with damage pattern observed within examined urban areas. Results show that, in some cases (Thessaloniki, Palermo), the HVNSR seems to be able to differentiate between areas previously shown to be associated with higher damage. In other cases however (Angra do Heroismo, Fabriano, Kalamata), the correlation is not statistically significant indicating thus the complex character of the parameters involved, implying that currently there is no a straightforward way that a value of HVNSR can correctly predict the extent to which a given region will be associated with increased damage.     

Influence of instruments on the H/V spectral ratios of ambient vibrations

For an optimal analysis of the H/V curve, it appears necessary to check the instrument signal to noise ratio in the studied frequency band, to ensure that the signal from the ground noise is well above the internal noise. We assess the reliability and accuracy of various digitizers, sensors and/or digitizer-sensor couples. Although this study is of general interest for any kind of seismological study, we emphasize the influence of equipment on H/V analysis results. To display the impact of the instrumental part on the H/V behavior, some series of tests have been carried out following a step-by-step procedure: first, the digitizers have been tested in the lab (sensitivity, internal noise...), then the three components sensors, still in the lab, and finally the usual user digitizers-sensors couple in lab and outdoors. In general, the digitizer characteristics, verified during this test, correspond well to the manufacturer specifications, however, depending on the digitizer, the quality of the digitized waveform can be very good to very poor, with variation from a channel to another channel (gain, time difference etc.). It appears very clearly that digitizers need a warming up time before the recording to avoid problems in the low-frequency range. Regarding the sensors, we recommend strongly to avoid the use of “classical” accelerometers (i.e., usual force balance technology). The majority of tested seismometers (broadband and short period, even 4.5 Hz) can be used without problems from 0.4 to 25 Hz. In all cases, the instrumentation should be checked first to verify that it works well for the defined study aim, but also to define its limit of use (frequency, sensitivity...).     

Polyphase thrusting and dyke emplacement in the central Southern Alps (Northern Italy)

The Triassic succession of the central Southern Alps (Italy) is stacked into several units bounded by south-verging low-angle thrust faults, which are related to two successive steps of crustal shortening. The thrust surfaces are cut by high-angle extensional and strike-slip faults, which controlled the emplacement of hypabissal magmatic intrusions that post-date thrusts motions. Intrusion ages based on SHRIMP U–Pb zircon dating span between 42 ± 1 and 39 ± 1 Ma, suggesting close time relationships with the earliest Adamello intrusion stages and, more in general, with the widespread calc-alkaline magmatism described in the Southern Alps. Fission-track ages of magmatic apatites are indistinguishable from U–Pb crystallization ages of zircons, suggesting that the intrusion occurred in country rocks already exhumed above the partial annealing zone of apatite (depth < 2–4 km). These data indicate that the central Southern Alps were already structured and largely exhumed in the Middle Eocene. Although we describe minor faults affecting magmatic bodies and local reactivations of older structures, no major internal deformations have occurred in the area after the Bartonian. Neogene deformations were instead concentrated farther south, along the frontal part of the belt.     

A roadmap for a dedicated Earth Science Grid platform

Due to its intensive data processing and highly distributed organization, the multidisciplinary Earth Science applications community is uniquely positioned for the uptake and exploitation of Grid technologies. Currently Enabling Grids for E-sciencE, and other large Grid infrastructures are already deployed and capable of operational services. So far however, the adoption and exploitation of Grid technology throughout the Earth Science community has been slower than expected. The Dissemination and Exploitation of GRids in Earth sciencE project, proposed by the European Commission to assist and accelerate this process in a number of different ways, had between its main goals the creation of a roadmap towards Earth Science Grid platform. This paper presents the resulting roadmap.     

Catchment properties,function, and conceptual model representation: is there a correspondence?

This study investigates the possible correspondence between catchment structure, as represented by perceptual hydrological models developed from fieldwork investigations, and mathematical model structures, selected on the basis of reproducing observed catchment hydrographs. Three Luxembourgish headwater catchments are considered, where previous fieldwork suggested distinct flow‐generating mechanisms and hydrological dynamics. A set of lumped conceptual model structures are hypothesized and implemented using the SUPERFLEX framework. Following parameter calibration, the model performance is examined in terms of predictive accuracy, quantification of uncertainty, and the ability to reproduce the flow–duration curve signature. Our key research question is whether differences in the performance of the conceptual model structures can be interpreted based on the dominant catchment processes suggested from fieldwork investigations. For example, we propose that the permeable bedrock and the presence of multiple aquifers in the Huewelerbach catchment may explain the superior performance of model structures with storage elements connected in parallel. Conversely, model structures with serial connections perform better in the Weierbach and Wollefsbach catchments, which are characterized by impermeable bedrock and dominated by lateral flow. The presence of threshold dynamics in the Weierbach and Wollefsbach catchments may favour nonlinear models, while the smoother dynamics of the larger Huewelerbach catchment were suitably reproduced by linear models. It is also shown how hydrologically distinct processes can be effectively described by the same mathematical model components. Major research questions are reviewed, including the correspondence between hydrological processes at different levels of scale and how best to synthesize the experimentalist's and modeller's perspectives. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

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.     

Interplay of S and As in Mekong Delta sediments during redox oscillations

The cumulative effects of periodic redox cycling on the mobility of As,Fe,and S from alluvial sediment to groundwater were investigated in bioreactor experiments.Two particular sediments from the alluvial floodplain of the Mekong Delta River were investigated:Matrix A(14 m deep)had a higher pyrite concentration than matrix B(7 m deep)sediments.Gypsum was present in matrix B but absent in matrix A.In the reactors,the sediment suspensions were supplemented with As(Ⅲ)and SO_4~(2-),and were subjected to three full-redox cycles entailing phases of nitrogen/CO_2,compressed air sparging,and cellobiose addition.Major differences in As concentration and speciation were observed upon redox cycling.Evidences support the fact that initial sediment composition is the main factor controlling arsenic release and its speciation during the redox cycles.Indeed,a high pyrite content associated with a low SO_4~(2-)content resulted in an increase in dissolved As concentrations,mainly in the form of As(Ⅲ),after anoxic half-cycles;whereas a decrease in As concentrations mainly in the form of As(Ⅴ),was instead observed after oxic half-cycles.In addition,oxic conditions were found to be responsible for pyrite and arsenian pyrite oxidation,increasing the As pool available for mobilization.The same processes seem to occur in sediment with the presence of gypsum,but,in this case,dissolved As were sequestered by biotic or abiotic redox reactions occurring in the Fe—S system,and by specific physico-chemical condition(e.g.pH).The contrasting results obtained for two sediments sampled from the same core show that many complexes and entangled factors are at work,and further refinement is needed to explain the spatial and temporal variability of As release to groundwater of the Mekong River Delta(Vietnam).     

HTP21/c–C2/c phase transition and kinetics of Fe2+–Mg order–disorder of an Fe-poor pigeonite: implications for the cooling history of ureilites

A natural Ca-poor pigeonite (Wo₆En₇₆Fs₁₈) from the ureilite meteorite sample PCA82506-3, free of exsolved augite, was studied by in situ high-temperature single-crystal X-ray diffraction. The sample, monoclinic P2₁/c, was annealed up to 1,093°C to induce a phase transition from P2₁/c to C2/c symmetry. The variation with increasing temperature of the lattice parameters and of the intensity of the b-type reflections (h + k = 2n + 1, present only in the P2₁/c phase) showed a displacive phase transition P2₁/c to C2/c at a transition temperature T _Tr = 944°C, first order in character. The Fe–Mg exchange kinetics was studied by ex situ single-crystal X-ray diffraction in a range of temperatures between the closure temperature of the Fe–Mg exchange reaction and the transition temperature. Isothermal disordering annealing experiments, using the IW buffer, were performed on three crystals at 790, 840 and 865°C. Linear regression of ln k _D versus 1/T yielded the following equation: ln k_\textD = - 3717( ±416)/T(K) + 1.290( ±0.378);    (R² = 0.988) \ln \,k_{\text{D}} = - 3717( \pm 416)/T(K) + 1.290( \pm 0.378);\quad (R^{2} = 0.988) . The closure temperature (T _c) calculated using this equation was ∼740(±30)°C. Analysis of the kinetic data carried out taking into account the e.s.d.'s of the atomic fractions used to define the Fe–Mg degree of order, performed according to Mueller’s model, allowed us to retrieve the disordering rate constants C ₀ K _dis⁺ for all three temperatures yielding the following Arrhenius relation: ln( C₀ K_\textdis ⁺ ) = ln K₀ - Q/(RT) = 20.99( ±3.74) - 26406( ±4165)/T(K);    (R² = 0.988) \ln \left( {C_{0} K_{\text{dis}}^{ + } } \right) = \ln \,K_{0} - Q/(RT) = 20.99( \pm 3.74) - 26406( \pm 4165)/T(K);\quad (R^{2} = 0.988) . An activation energy of 52.5(±4) kcal/mol for the Fe–Mg exchange process was obtained. The above relation was used to calculate the following Arrhenius relation modified as a function of X _Fe (in the range of X _Fe = 0.20–0.50): ln( C₀ K_\textdis ⁺ ) = (21.185 - 1.47X_\textFe ) - \frac(27267 - 4170X_\textFe )T(K) \ln \left( {C_{0} K_{\text{dis}}^{ + } } \right) = (21.185 - 1.47X_{\text{Fe}} ) - {\frac{{(27267 - 4170X_{\text{Fe}} )}}{T(K)}} . The cooling time constant, η = 6 × 10⁻¹ K⁻¹ year⁻¹ calculated on the PCA82506-3 sample, provided a cooling rate of the order of 1°C/min consistent with the extremely fast late cooling history of the ureilite parent body after impact excavation.