Use of microtremor measurement for assessing site effects in Northern Belgium – interpretation of the observed intensity during the MS = 5.0 June 11 1938 earthquake

To provide quantitative information on site effects in the northern part of Belgium, forty-seven H/V microtremor measurements were performed with 5 second seismometers over an area of about 15.000 km². Most of the results show a northward regular increase of the fundamental period in agreement with the augmentation of the Mesozoic and Cainozoic soft sedimentsthickness from a few meters 40 km south of Brussels to 900 m at the Netherlands-Belgium border. The measured resonance frequency values wereconsistent with theoretical computations performed at different sites onthe basis of existing information and shallow seismic experiments. At oneparticular site (Uccle) where borehole data were available, microtremor measurements using an array of four seismological stations with differentapertures allowed to obtain the low frequency part of the Rayleigh wave dispersion curve, extending the range covered by the analysis of surfacewaves artificially generated. The Vs profile derived from the surface waveinversion corroborates the 1 Hz natural frequency of the site. Comparison of these results with the macroseismic information concerning the M_S = 5.0 1938 earthquake which occurred 50 km west of Brussels,confirmed the hypothesis that the geological structure of the Brabant massifis likely to control damage distribution during such an earthquake. Comparisonbetween the intensity map of the 1938 earthquake and the resonance period ofsediments obtained by our microtremor study shows a clear relation betweenthe two parameters. During the 1938 earthquake, site effects played a prominent role due to the dimension of the source whose corner frequency wasabout 1 Hz.     

On the Limnology of Deep and Saline Lake Burdur in Turkey

Limnological characteristics of Lake Burdur in Lake District in south‐western Turkey are presented. It is a deep, tectonic (estimated max. depth 100 m), athalassic, highly alkaline, and saline lake. A set of physical and chemical variables was monitored, phyto‐ and zooplankton was sampled from surface layer of the lake during 1997. Physico‐chemical variables indicated that the lake is hyposaline and composed of some hydrochemically different water layers formed by groundwater sources located on the bottom of the lake. The phytoplankton composition of Lake Burdur consisted of Cyanophyta, Chlorophyta, Bacillariophyta, Dinophyta, and Chrysophyta. The abundance and number of species of Chlorophyta and Cyanophyta were higher than the other taxa. The zooplankton composition of the Lake consisted of Rhizopoda, Rotifera, and Crustacea. Number of species of Rotifera was higher than the other taxa. The diversities of the phyto‐ and zooplankton were calculated according to the Shannon‐Weaver diversity index. The diversity of each group was found to be low in the lake.     

Some Morphological Aspects and Hydrological Characterization of the Tagus Floods in the Santarém Region, Portugal

In Portugal, few studies have been made of the historical floods of the RiverTagus (the longest in the Iberian Peninsula). This fact led to the study of theLower Tagus (Santarém region) using written historical documents,cartographic documents and hydrological data, consisting mainly of waterlevel records. With the support of the historical documents and the analysisof all the maps, it was possible to verify that the human intervention has beenfundamental in the morphological changes of the Tagus' channel. It becameaware that the river changed from a braided to a single channel with alternatebars. From the hydrological data, return periods for the floods were determined,using the Pearson Type 3 distribution. The selection of the most important floods,from 1855 to 1998, enabled to build a ``flood hazard scale' for this region.     

Indicator Simulation Accounting for Multiple-Point Statistics

Geostatistical simulation aims at reproducing the variability of the real underlying phenomena. When nonlinear features or large-range connectivity is present, the traditional variogram-based simulation approaches do not provide good reproduction of those features. Connectivity of high and low values is often critical for grades in a mineral deposit. Multiple-point statistics can help to characterize these features. The use of multiple-point statistics in geostatistical simulation was proposed more than 10 years ago, on the basis of the use of training images to extract the statistics. This paper proposes the use of multiple-point statistics extracted from actual data. A method is developed to simulate continuous variables. The indicator kriging probabilities used in sequential indicator simulation are modified by probabilities extracted from multiple-point configurations. The correction is done under the assumption of conditional independence. The practical implementation of the method is illustrated with data from a porphyry copper mine.     

An attenuation relationship based on Turkish strong motion data and iso-acceleration map of Turkey

This paper presents an attenuation relationship of peak ground acceleration (PGA) derived from Turkish strong motion data for rock, soil and soft soil sites and an iso-acceleration map of Turkey based on this relationship. For the purpose, among all the three-component accessible records, 221 records from 122 earthquakes that occurred in Turkey between 1976 and November 2003 were selected. The database was compiled for earthquakes with moment magnitudes (M_w) and PGA values ranging between 4.1 and 7.5, and 20 and 806 gal, and distances to epicenter considered in the database were between 5 and 100 km. From the regression analysis of the data, an attenuation equation of PGA considering rock, soil and soft soil conditions was developed. The PGA values predicted from the equation suggested in this study and those both from a few domestic equations and some imported equations were compared. In addition, an iso-acceleration map of Turkey was constructed using the suggested attenuation equation and considering both known active faults and epicenter locations of the earthquakes that have occurred in Turkey.     

Low-frequency 3-D seismic survey of upper crustal magmatic intrusions in the northeastern Pannonian basin of Hungary

A specially designed 700-km² grid survey, deploying 1000 regularly distributed low-frequency seismic recording systems, successfully investigated one of the most complex geologic environments of the Pannonian basin. The wide-angle signals penetrated through over 1000 m of multi-phase igneous lithology and recognized, for the first time, the underlying enigmatic Permian to Early Triassic basement rocks. Tomographic inversion of the first arrival grid data resulted in determination of an accurate three-dimensional (3-D) velocity field, to a depth of 4 km. The anomalous changes of the spatial velocity data outline the regional extent of the Late Miocene magmatic intrusions, which are covered by over 2000 m of Mid-Miocene to Pleistocene clastics. Complex relationship was found between the surface potential data and the intrusive bodies. This multi-faceted geophysical data analysis established a functional technique for mapping a subsurface with intricate acoustic and structural complexity.     

Organic matter in a low productivity anoxic intraplatform basin in the Triassic Tethys

This paper presents isotopic, bulk geochemical and biomarker data measured on organic matter accumulated in a narrow extensional basin developed at the oceanward margin of the huge Triassic carbonate platform in the Alps–Appenines domain. The integration of isotope signatures, organic petrographical and biomarker evidence together with the composition of kerogen pyrolysates suggests immature organic matter predominantly of algal origin with a minor, but not negligible, higher plant derived and moderate bacterial contribution for the entire sequence. The mineral sources are dominated by platform-derived subtidal Dachstein limestone with a minor palaeosol input and a moderate contribution of autochtonous quartz. Nevertheless, parallel variations observed in the mineral content, as well as in the amount and the quality of the organic matter reflect variations in the palaeoenvironment. The increased humidity, existing in the period of the accumulation of the upper section of the sequence, led to the restriction of dolomitization. A slightly greater higher plant derived contribution, in this section, is evidenced by the composition of bitumen and the results on GC and GC/MS analyses on the non-aromatic hydrocarbon fraction of bitumen. Moreover, the climate-induced weathering enchanced the primary productivity and resulted in a pronounced increase in the TOC content. The average estimated value of the planktonic productivity is about three times higher for the calcite-rich sequence than the dolomite-rich one, being 44.2 and 15.3 tC_org/m²/Ma, respectively. The low to moderate planktonic productivity shows that anoxic conditions, observed for the entire succession, are a consequence of the stagnant water stratification rather than high planktonic productivity. Depth trends in the data measured on kerogens (HI, OI, δ¹³Corg values, composition of pyrolysate) together with the δ¹⁸O excursions and Δδ¹³C values appear to be controlled by sea-level fluctuations. Consistent with the high abundance of alkyltiophenes in kerogen pyrolysates, the high S_org/C ratios (ranging between 0.05 and 0.10) suggest the importance of natural sulfurization in the formation of the sulfur-rich type-II-S kerogen occurring in all of the samples.     

Relationship between solar radiation and the development and morphology of small cirque glaciers (Maladeta Mountain massif, Central Pyrenees, Spain)

Abstract Spatial and temporal variations in radiative fluxes influence glacier mass‐balance in mountain areas. The primary goal of this study was to assess differences in solar radiation on three glacial cirques located in the Maladeta Mountain massif (Central Spanish Pyrenees), and analyse their implications on glacier development and morphology. A quantitative approach is adopted to obtain the values of solar radiation (direct, diffuse and global radiation), combining several field data parameters (measured at 55 control‐points) with the solar radiation modelling package Ecosim. The data obtained confirm that the morphologies of the glacial bodies developed in the three cirques have a good correlation with the spatial variation on solar radiation inputs, favouring also the conservation (Aneto and Coronas cirques) or total vanishing (Llosás cirque) of the glacial remnants analysed here. The study shows how strongly in this Alpine‐Mediterranean context solar radiation — firstly as a function of latitude and time of year, and locally as a function of topographic slope, aspect and shadowing — controls the mass‐balance and the spatial distribution of melting in small glaciers, having an effect on the development of their morphologies.     

Image processing of 2D resistivity data for imaging faults

A methodology to locate automatically limits or boundaries between different geological bodies in 2D electrical tomography is proposed, using a crest line extraction process in gradient images. This method is applied on several synthetic models and on field data set acquired on three experimental sites during the European project PALEOSIS where trenches were dug. The results presented in this work are valid for electrical tomographies data collected with a Wenner-alpha array and computed with an l₁ norm (blocky inversion) as optimization method. For the synthetic cases, three geometric contexts are modelled: a vertical and a dipping fault juxtaposing two different geological formations and a step-like structure. A superficial layer can cover each geological structure. In these three situations, the method locates the synthetic faults and layer boundaries, and determines fault displacement but with several limitations. The estimated fault positions correlate exactly with the synthetic ones if a conductive (or no superficial) layer overlies the studied structure. When a resistive layer with a thickness of 6 m covers the model, faults are positioned with a maximum error of 1 m. Moreover, when a resistive and/or a thick top layer is present, the resolution significantly decreases for the fault displacement estimation (error up to 150%). The tests with the synthetic models for surveys using the Wenner-alpha array indicate that the proposed methodology is best suited to vertical and horizontal contacts. Application of the methodology to real data sets shows that a lateral resistivity contrast of 1:5–1:10 leads to exact faults location. A fault contact with a resistivity contrast of 1:0.75 and overlaid by a resistive layer with a thickness of 1 m gives an error location ranging from 1 to 3 m. Moreover, no result is obtained for a contact with very low contrasts (1:0.85) overlaid by a resistive soil. The method shows poor results when vertical gradients are greater than horizontal ones. This kind of image processing technique should be systematically used for improving the objectiveness of tomography interpretation when looking for limits between geological objects.