Bayesian Spatial Modelling of Gamma Ray Count Data

Gamma ray logging is a method routinely employed by geophysicists and environmental engineers in site geology evaluations. Modelling of gamma ray data from individual boreholes assists in the local identification of major lithological changes; modelling these data from a network of boreholes assists with lithological mapping and spatial stratigraphic correlation. In this paper we employ Bayesian spatial partition models to analyse gamma ray data spatially. In particular, a spatial partition is defined via a Voronoi tessellation and the mean intensity is assumed constant in each cell of the partition. The number of vertices generating the tessellation as well as the locations of vertices are assumed unknown, and uncertainty about these quantities is described via a hierarchical prior distribution. We describe the advantages of the spatial partition modelling approach in the context of smoothing gamma ray count data and describe an implementation that may be extended to the fitting of a more general model than a constant mean within each cell of the partition. As an illustration of the methodology we consider a data set collected from a network of eight boreholes, which is part of a geophysical study to assist in mapping the lithology of a site. Gamma ray logs are linked with geological information from cores and the spatial analysis of log data assists with predicting the lithology at unsampled locations.     

Damage Survey of Water Supply Systems and Fragility Curve of PVC Water Pipelines in the Chi–Chi Taiwan Earthquake

Substantial damage to water supply systems, including water delivery pipelines, water treatment plants, reservoirs, and water storage tanks, was reported after the 1999 Chi–Chi Taiwan Earthquake. This paper first summarizes the damage survey and then presents the results of seismic fragility analysis for underground pipelines. Construction blueprints of the water delivery pipelines and repair work orders of 11 townships and cities in the disastrous area were digitized into a Geographical Information System (GIS) for analysis and assessment. With the aid of the GIS system, we found that PVC pipes made up 86% of water delivery pipelines while steel, cast iron, ductile iron, PE and others took the rest. Therefore, this paper focuses on the fragility analysis of PVC pipes. Three different methods were applied to derive the fragility relations between the PVC water pipes having nominal diameters (approximately inner diameters) greater than or equal to 65 mm and earthquake intensity parameters such as peak ground acceleration and peak ground velocity. The results were then examined with those of other countries. The discrepancy between our results and the empirical equation used by HAZUS, an earthquake loss estimation software developed by the Federal Emergency Management Agency was not significant.     

Assessing the Vibrational Frequencies of the Cathedral of Cologne (Germany) by Means of Ambient Seismic Noise Analysis

Ambient seismic noise measurements were conducted inside the Cathedral of Cologne (Germany) for assessing its frequencies of vibration and for checking whether these occur in the range where soil amplification is expected. If this is the case, damages may increase in case of an earthquake due to an increased structural response of the building. Analysis of the ratio between the horizontal and vertical components of the spectra recorded at stations located inside the building as well as the ratio between the corresponding components of the spectra recorded simultaneously inside the building and at a reference station placed in the basement of the cathedral indicated several modes of vibration. Facilitated by these results an assessment of the seismic vulnerability was attempted for a 2D ground motion scenario using the finite element method.     

Site Effects Study in the Serravalle Village (NE Italy) by Means of Spectral Ratios Analyses Relatively to Multiple Reference Sites

The area of Serravalle, sited in the northern part of the town of Vittorio Veneto (TV), NE Italy, has been the target of a seismic microzonation campaign. 10 seismic stations have been deployed for a 7 months period to record in continuous mode. Three stations were installed on bedrock outcrops and seven on sedimentary sites with variable cover thickness. Spectral analyses have been performed on the collected data-set using the Generalized Inversion Technique (GIT, e.g. Andrews, 1986). In particular, spectral ratios have been calculated for each station relatively to the average of the three reference, bedrock sites. The spectral ratios provide quantitative estimates of the seismic motion amplifications which occur in each of the monitored sites. Two sites show high values of amplification, 5 times larger than signal amplitude at the reference sites, in correspondence of well discernible peak frequencies of 5 Hz. Results for the other stations show smaller amounts of site amplification spreading over a broad range of frequencies. Sites where the highest amplifications were recorded all lie on the left bank of the Meschio River and in areas farther away from its outlet into the plain correlating with the presence of thick layers of Quaternary deposits.     

New Generation of Probabilistic Seismic Hazard Assessment for the Area Cologne/Aachen Considering the Uncertainties of the Input Data

A new earthquake catalogue for central, northern and northwestern Europe with unified M_w magnitudes, in part derived from chi-square maximum likelihood regressions, forms the basis for seismic hazard calculations for the Lower Rhine Embayment. Uncertainties in the various input parameters are introduced, a detailed seismic zonation is performed and a recently developed technique for maximum expected magnitude estimation is adopted and quantified. Applying the logic tree algorithm, resulting hazard values with error estimates are obtained as fractile curves (median, 16% and 84% fractiles and mean) plotted for pga (peak ground acceleration; median values for Cologne 0.7 and 1.2 m/s² for probabilities of exceedence of 10% and 2%, respectively, in 50 years), 0.4 s (0.8 and 1.5 m/s²) and 1.0 s (0.3 and 0.5 m/s²) pseudoacclerations, and intensity (I₀ = 6.5 and 7.2). For the ground motion parameters, rock foundation is assumed. For the area near Cologne and Aachen, maps show the median and 84% fractile hazard for 2% probability of exceedence in 50 years based on pga (maximum median value about 1.5 m/s²), and 0.4 s (>2 m/s²) and 1.0 s (about 0.8 m/s²) pseudoaccelerations, all for rock. The pga 84% fractile map also has a maximum value above 2 m/s² and shows similarities with the median map for 0.4 s. In all maps, the maximum values fall within the area 6.2–6.3° E and 50.8–50.9° N, i.e., east of Aachen.     

Crustal structure across the Grand Banks–Newfoundland Basin Continental Margin – II. Results from a seismic reflection profile

New multichannel seismic reflection data were collected over a 565 km transect covering the non-volcanic rifted margin of the central eastern Grand Banks and the Newfoundland Basin in the northwestern Atlantic. Three major crustal zones are interpreted from west to east over the seaward 350 km of the profile: (1) continental crust; (2) transitional basement and (3) oceanic crust. Continental crust thins over a wide zone (∼160 km) by forming a large rift basin (Carson Basin) and seaward fault block, together with a series of smaller fault blocks eastwards beneath the Salar and Newfoundland basins. Analysis of selected previous reflection profiles (Lithoprobe 85-4, 85-2 and Conrad NB-1) indicates that prominent landward-dipping reflections observed under the continental slope are a regional phenomenon. They define the landward edge of a deep serpentinized mantle layer, which underlies both extended continental crust and transitional basement. The 80-km-wide transitional basement is defined landwards by a basement high that may consist of serpentinized peridotite and seawards by a pair of basement highs of unknown crustal origin. Flat and unreflective transitional basement most likely is exhumed, serpentinized mantle, although our results do not exclude the possibility of anomalously thinned oceanic crust. A Moho reflection below interpreted oceanic crust is first observed landwards of magnetic anomaly M4, 230 km from the shelf break. Extrapolation of ages from chron M0 to the edge of interpreted oceanic crust suggests that the onset of seafloor spreading was ∼138 Ma (Valanginian) in the south (southern Newfoundland Basin) to ∼125 Ma (Barremian–Aptian boundary) in the north (Flemish Cap), comparable to those proposed for the conjugate margins.     

Interpretation of subhorizontal crustal reflections by metamorphic and rheologic effects in the eastern part of the Pannonian Basin

The geologic origin of subhorizontal reflections, often observed in crustal seismic sections, was investigated by establishing metamorphic facies and strength of rocks in depth, and correlating these properties to seismic reflection sections from eastern Hungary. Estimation of the depths of metamorphic mineral stability zones utilized the principles developed by Fyfe et al. and known geothermal data of the area. The strength versus depth profile was derived by relating local seismic P -wave interval velocities to Meissner et al. 's activation energy. The results show that the series of subhorizontal reflections, observed in the Pannonian Basin, are a consequence of combined metamorphic and rheologic changes in depths. The synthesis of the integrated data set suggests that the retrograde alteration of the pre-Tertiary basement above the percolation threshold was made possible by the softening effect of shear zones and their water-conducting capacity. The subhorizontal reflections of highest energy, of the consolidated crust below the percolation threshold, originate in the depths of greenschist, amphibolite and granulite metamorphic mineral facies, which were formed in geothermal and pressure conditions similar to those existing today. These results imply the overprint of earlier (Variscan) metamorphic sequences of the crust by more recent retrograde metamorphic processes.