Recent Advances in Prediction and Processing of Strong Ground Motions

We present an overview of our recent results on utilizing small earthquakes in the earthquake engineering practice. Site-specific ground motion time-histories of large earthquakes can be successfully simulated using recordings of small earthquakes which are often referred to as 'empirical Green's functions' in seismology. Another important practical problem is whether and how these observations can be used in seismic risk studies which are based on empirical attenuation relations for ground motion parameters. We study a possibility of extrapolating attenuation relations for small earthquakes, to larger magnitudes using the data from the Garner Valley downhole array in Southern California. Finally we introduce efficient ground motion processing techniques in frequency- and time-domains and apply them to site response estimation.     

Seismic Ground Motion Expected for the Eastern District of Naples

The seismic ground motion of a test area in the eastern district of Naples is computed with a hybrid technique based on the mode summation and the finite difference methods. This technique allows us the realistic modelling of source and propagation effects, including local soil conditions. In the modelling, we consider the 1980 Irpinia earthquake, a good example of strong shaking for the area of Naples, which is located about 90 km from the epicenter.The detailed geological setting is reconstructed from a large number of drillings. The sub-soil is mainly formed by alluvial (ash, stratified sand and peat) and pyroclastic materials overlying a pyroclastic rock (yellow neapolitan tuff), representing the neapolitan bedrock. The detailed information available on mechanical properties of the sub-soil and its geometry warrants the application of the sophisticated hybrid technique.As expected, the sedimentary cover causes an increase of the signal's amplitudes and duration. If thin peat layers are present, the amplification effects are reduced, and the peak ground accelerations are similar to those observed for the bedrock model. This can be explained by the backscattering of wave energy at such layers, that tend to seismically decouple the upper from the lower part of the structure.For SH-waves, the influence of the variations of the S-wave velocities on the spectral amplification is studied, by considering locally measured velocities and values determined from near-by down-hole measurements. The comparison between the computed spectral amplifications confirms the key role of an accurate determination of the seismic velocities of the different layers.The comparison performed between a realistic 2-D seismic response and a standard 1-D response, based on the vertical propagation of waves in a plane layered structure, shows considerable difference, from which it is evident that serious caution must be taken in the modelling of expected ground motion at a specific site.     

Thermodynamic Modelling of Ore-Forming Mechanism of the Changkeng Gold-Silver Deposits in Guangdong Province

The Changkeng gold-silver deposits consist of a sediment-hosted, disseminated gold deposit and a replacement-type silver deposit. The mineralizations of gold and silver are zoned and closely related to the silicification of carbonate and clastic rocks, so that siliceous ores dominate in the deposit. The mineralizing temperature ranges mainly from 300 to 170℃, and K+, Na+, Ca2+, Mg2+, and Cl- are the major ions in the ore-forming fluid. Calculations of distribution of metal complexes show that gold is mainly transported by hydrosulphide complexes, but chloride complexes of silver, iron, lead, and zinc, which are transformed into hydroxyl and hydrosulphide complexes under neutral to weak-alkaline circumstances in the late stage, predominate in the ore-forming solutions. Water-rock interaction is confirmed to be the effective mechanism for the formation of silver ores by computer modelling of reaction of hydrothermal solution with carbonate rocks. The solubility analyses demonstrate that the precipitation     

In-plane and anti-plane strong shaking of soil systems and structures

The concept of in-plane and anti-plane shaking is introduced with a rigid block on a plane surface with Coulomb friction. Using a hypoplastic constitutive relation to model the mechanical behaviour of the soil, numerical solutions for a rigid block on a thin dry or saturated soil layer are obtained. The coupled nature of dynamic problems involving granular materials is shown, i.e. the motion of the block changes the soil state—skeleton stresses and density—which in turn affects the block motion. Motions of the block as well as soil response can be more realistically calculated by the new model. The same constitutive equation is applied to the numerical simulation of the propagation of plane waves in homogeneous and layered level soil deposits induced by a wave coming from below. Experiments with a novel laminar shake box as well as real seismic records from well-documented sites during strong earthquakes are used to verify the adequacy of the hypoplasticity-based numerical model for the prediction of soil response during strong earthquakes. The response of a homogeneous earth dam subjected to in-plane and anti-plane shaking is investigated numerically. In-plane and anti-plane shaking is shown to cause nearly the same spreading of a sand dam under drained conditions, whereas under undrained conditions anti-plane shaking causes stronger spreading of the dam. The dynamic behaviour of a breakwater founded on rockfill and soft clay during the 1995 Kobe earthquake is back-calculated to show the good performance of the proposed numerical model also with a structure. Section 9 deals with buildings on mattresses of densified cohesionless soils or fine-grained soils with granular columns, slopes with ‘hidden’ dams and structures on piles traversing clayey slopes to show the suitability of hypoplasticity-based models for the earthquake-resistant design and safety assessment of geotechnical systems.     

Tectonic evolution of the Lachlan Fold Belt, southeastern Australia: constraints from coupled numerical models of crustal deformation and surface erosion driven by subduction of the underlying mantle

We have used a coupled thermo-mechanical finite-element (FE) model of crustal deformation driven by mantle/oceanic subduction to demonstrate that the tectonic evolution of the Lachlan Fold Belt (LFB) during the Mid-Palaeozoic (Late Ordovician to Early Carboniferous) can be linked to continuous subduction along a single subduction zone. This contrasts with most models proposed to date which assume that separate subduction zones were active beneath the western, central and eastern sections of the Lachlan Orogen. We demonstrate how the existing data on the structural, volcanic and erosional evolution of the Lachlan Fold Belt can be accounted for by our model. We focus particularly on the timing of fault movement in the various sectors of the orogen. We demonstrate that the presence of the weak basal decollement on which most of the Lachlan Fold Belt is constructed effectively decouples crustal structures from those in the underlying mantle. The patterns of faulting in the upper crust appears therefore to be controlled by lateral strength contrasts inherited from previous orogenic events rather than the location of one or several subduction zones. The model also predicts that the uplift and deep exhumation of the Wagga-Omeo Metamorphic Belt (WOMB) is associated with the advection of this terrane above the subduction point and is the only tectonic event that gives us direct constraints on the location of the subduction zone. We also discuss the implications of our model for the nature of the basement underlying the present-day orogen.     

Bearing capacity of rock over mined cavities in Nottingham

A significant geohazard is created in Nottingham, UK, by hundreds of man-made caves cut in the weak sandstone beneath the city centre. Stability of the caves has been assessed by a single full-scale loading test, by numerical modelling with FLAC and by physical modelling in plaster. For typical caves 4 m wide, bearing capacity of the rock roof rises from 2 MPa where it is 1 m thick to 8 MPa where 3 m thick. Stability decreases over wider caves and where the loading pad edge is over the edge of the cave. Numerical modelling of a very wide cave revealed the failure mechanisms and also showed that an internal support wall increased roof bearing capacity by 50%. Local building regulations that require 3-5 m of rock cover over the sandstone caves appear to be conservative. In stronger rocks, including karstic limestone, a guideline that cover thickness exceeds 70% of the cave width appears to be appropriate.     

The Significance of Multiple Saturation Points in the Context of Polybaric Near-fractional Melting

Experimental petrologists have successfully located basalticliquid compositions parental to mid-ocean ridge basalt thatare, within experimental resolution, multiply saturated withthree-phase harzburgite or four-phase lherzolite assemblageson their liquidus at some elevated pressure. Such an experimentalresult is a necessary consequence of any paradigm in which eruptedbasalts derive from single-batch primary liquids that equilibratewith a mantle residue and undergo no subsequent magma mixingbefore differentiation and eruption. Here we investigate whether,conversely, such evidence of multiple saturation is sufficientto exclude dynamic melting models wherein increments of meltare mixed after segregation from residues, during melt transportor in magma chambers. Using two independent models of crystal–liquidequilibria to simulate polybaric near-fractional peridotitemelting, we find that aggregate liquids from such melting processescan display near-intersections of liquidus surfaces too closeto distinguish experimentally from exact multiple saturationpoints. Given uncertainties in glass compositions, fractionationcorrections, experimental temperature and pressure conditions,and achievement of equilibrium, these results suggest that polybaricmixtures can in fact masquerade as mantle-equilibrated single-batchprimary liquids. Multiple saturation points on the liquidussurfaces of primitive basalts do, however, preserve informationabout the average pressure of extraction of their constituentincrements of liquid. KEY WORDS: mantle melting; basaltic volcanism; experimental igneous petrology; thermodynamic modelling; inverse method     

Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations

Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature‐induced grain coarsening. Despite this large‐scale trend, a considerable grain size scatter occurs on the outcrop‐scale indicating local influence of second‐order effects such as thermal perturbations, fluid flow and second‐phase particles. Second‐phase particles, whose sizes range from nano‐ to the micron‐scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron‐scale second‐phase particles, or to samples containing a large number of nano‐scale particles. The obtained data set suggests that the second phases induce a temperature‐controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D = C₂ e^Q*/RT(d_p/f_p)^m*, where C₂ is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio d_p/f_p, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C₂ and Q*. Besides variations in the average grain size, the presence of second‐phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian‐type distributions of the pure samples. In contrast, fluid‐enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable δ¹⁸O and δ¹³C isotope ratios in fluid‐affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid‐induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1. To summarize, considerable grain size variations of up to one order of magnitude can locally result from second‐order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies.     

Analysis of geomagnetic field along seismic profile P4 of the International Project POLONAISE'97

This paper presents relative secular variations of the total intensity of the geomagnetic field against a background of results of magnetic anomaly interpretation along seismic profile P4. Profile P4 crosses a Variscan folding zone in the Paleozoic Platform (PLZ), the Trans-European Suture Zone (TESZ), and the Polish part of the East European Craton (EEC). Secular geomagnetic field variations measured in 1966–2000 along a line adjacent to seismic profile P4 were analysed. The study of secular variations, reduced to the base recordings at the Belsk Magnetic Observatory, showed that the growth of geomagnetic field at the East European Craton was slower than in the Trans-European Suture Zone and the Paleozoic Platform.A 2D crustal magnetic model was interpreted as a result of magnetic modelling, in which seismic, geological and geothermal data were also used. The modelling showed that there were significant differences in the magnetic model for geotectonic units, which had been earlier determined based on deep seismic survey data. It should be noted that a fundamental change of trend of the relative secular variations was observed at the slope of the Precambrian Platform. After analysing the geomagnetic field observed along profile P4, the hypothesis that the contact between Phanerozoic and Precambrian Europe lies in Poland's territory can be proven.