Seismic Tomography and Earthquake Mechanism Beneath the Central Fennoscandian Shield

The aim of this thesis was to study the seismic tomography structure of the earth s crust together with earthquake distribution and mechanism beneath the central Fennoscandian Shield, mainly in southern and central Finland. The earthquake foci and some fault plane solutions are correlated with 3-D images of the velocity tomography. The results are discussed in relation to the stress field of the Shield and with other geophysical, e.g. geomagnetic, gravimetric, tectonic, and anisotropy studies of the Shield. The earthquake data of the Fennoscandian Shield has been extracted from the Nordic earthquake parameter data base which was founded at the time of inception of the earthquake catalogue for northern Europe. Eight earlier earthquake source mechanisms are included in a pilot study on creating a novel technique for calculating an earthquake fault plane solution. Altogether, eleven source mechanisms of shallow, weak earthquakes are related in the 3-D tomography model to trace stresses of the crust in southern and central Finland. The earthquakes in the eastern part of the Fennoscandian Shield represent low-active, intraplate seismicity. Earthquake mechanisms with NW-SE oriented horizontal compression confirm that the dominant stress field originates from the ridge-push force in the North Atlantic Ocean. Earthquakes accumulate in coastal areas, in intersections of tectonic lineaments, in main fault zones or are bordered by fault lines. The majority of Fennoscandian earthquakes concentrate on the south-western Shield in southern Norway and Sweden. Onwards, epicentres spread via the ridge of the Shield along the west-coast of the Gulf of Bothnia northwards along the Tornio River - Finnmark fault system to the Barents Sea, and branch out north-eastwards via the Kuusamo region to the White Sea Kola Peninsula faults. The local seismic tomographic method was applied to find the terrane distribution within the central parts of the Shield the Svecofennian Orogen. From 300 local explosions a total of 19765 crustal Pg- and Sg-wave arrival times were inverted to create independent 3-D Vp and Vs tomographic models, from which the Vp/Vs ratio was calculated. The 3-D structure of the crust is presented as a P-wave and for the first time as an S-wave velocity model, and also as a Vp/Vs-ratio model of the SVEKALAPKO area that covers 700x800 km2 in southern and central Finland. Also, some P-wave Moho-reflection data was interpolated to image the relief of the crust-mantle boundary (i.e. Moho). In the tomography model, the seismic velocities vary smoothly. The lateral variations are larger for Vp (dVp =0.7 km/s) than for Vs (dVs =0.4 km/s). The Vp/Vs ratio varies spatially more distinctly than P- and S-wave velocities, usually from 1.70 to 1.74 in the upper crust and from 1.72 to 1.78 in the lower crust. Schist belts and their continuations at depth are associated with lower velocities and lower Vp/Vs ratios than in the granitoid areas. The tomography modelling suggests that the Svecofennian Orogen was accreted from crustal blocks ranging in size from 100x100 km2 to 200x200 km2 in cross-sectional area. The intervening sedimentary belts have ca. 0.2 km/s lower P- and S-wave velocities and ca. 0.04 lower Vp/Vs ratios. Thus, the tomographic model supports the concept that the thick Svecofennian crust was accreted from several crustal terranes, some hidden, and that the crust was later modified by intra- and underplating. In conclusion, as a novel approach the earthquake focal mechanism and focal depth distribution is discussed in relation to the 3-D tomography model. The schist belts and the transformation zones between the high- and low-velocity anomaly blocks are characterized by deeper earthquakes than the granitoid areas where shallow events dominate. Although only a few focal mechanisms were solved for southern Finland, there is a trend towards strike-slip and oblique strike-slip movements inside schist areas. The normal dip-slip type earthquakes are typical in the seismically active Kuusamo district in the NE edge of the SVEKALAPKO area, where the Archean crust is ca. 15-20 km thinner than the Proterozoic Svecofennian crust. Two near vertical dip-slip mechanism earthquakes occurred in the NE-SW junction between the Central Finland Granitoid Complex and the Vyborg rapakivi batholith, where high Vp/Vs-ratio deep-set intrusion splits the southern Finland schist belt into two parts in the tomography model.

Influence of rock depth on seismic site classification for shallow bedrock regions

Seismic site classifications are used to represent site effects for estimating hazard parameters (response spectral ordinates) at the soil surface. Seismic site classifications have generally been carried out using average shear wave velocity and/or standard penetration test n-values of top 30-m soil layers, according to the recommendations of the National Earthquake Hazards Reduction Program (NEHRP) or the International Building Code (IBC). The site classification system in the NEHRP and the IBC is based on the studies carried out in the United States where soil layers extend up to several hundred meters before reaching any distinct soil-bedrock interface and may not be directly applicable to other regions, especially in regions having shallow geological deposits. This paper investigates the influence of rock depth on site classes based on the recommendations of the NEHRP and the IBC. For this study, soil sites having a wide range of average shear wave velocities (or standard penetration test n-values) have been collected from different parts of Australia, China, and India. Shear wave velocities of rock layers underneath soil layers have also been collected at depths from a few meters to 180 m. It is shown that a site classification system based on the top 30-m soil layers often represents stiffer site classes for soil sites having shallow rock depths (rock depths less than 25 m from the soil surface). A new site classification system based on average soil thickness up to engineering bedrock has been proposed herein, which is considered more representative for soil sites in shallow bedrock regions. It has been observed that response spectral ordinates, amplification factors, and site periods estimated using one-dimensional shear wave analysis considering the depth of engineering bedrock are different from those obtained considering top 30-m soil layers.

Seismic Bearing Capacity of Shallow Embedded Foundations on a Sloping Ground Surface

By using the lower-bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization, bearing-capacity factors, N-c and N-gamma q, with an inclusion of pseudostatic horizontal seismic body forces, have been determined for a shallow embedded horizontal strip footing placed on sloping ground surface. The variation of N-c and N-gamma q with changes in slope angle (beta) for different values of seismic acceleration coefficient (k(h)) has been obtained. The analysis reveals that irrespective of ground inclination and the embedment depth of the footing, the factors N-c and N-gamma q decrease quite considerably with an increase in k(h). As compared with N-c, the factor N-gamma q is affected more extensively with changes in k(h) and beta. Unlike most of the results reported in literature for the seismic case, the present computational results take into account the shear resistance of soil mass above the footing level. An increase in the depth of the embedment leads to an increase in the magnitudes of both N-c and N-gamma q. (C) 2014 American Society of Civil Engineers.

Seismic analyses of shallow tunnels by dynamic centrifuge tests and finite elements

The increments of internal forces induced in a tunnel lining during earthquakes can be assessed with several procedures at different levels of complexity. However, the substantial lack of well-documented case histories still represents a difficulty in order to validate any of the methods proposed in literature. To bridge this gap, centrifuge model tests were carried out on a circular aluminium tunnel located at two different depths in dense and loose dry sand. Each model has been instrumented for measuring soil motion and internal loads in the lining and tested under several dynamic input signals. The tests performed represented an experimental benchmark to calibrate dynamic analyses with different approaches to account for soil-tunnel kinematic interaction. © 2009 IOS Press.

Isolating shallow foundations from seismic loading

Most modern design codes do not allow for movement between a shallow foundation and the underlying soil during seismic loading. Consequently, the full magnitude of seismic energy is transmitted from the soil to the foundation during an earthquake. This energy either has to be dissipated before reaching the superstructure via engineering solutions such as base isolation systems, or the structure itself must withstand the full impact of the earthquake resulting in high material usage and expensive design. However, the inherent hysteric behaviour of soil can be used to isolate a foundation from the underlying soil. As part of a study into the soil-structure-interaction of shallow foundations, methods to optimise foundation isolation were investigated. In this paper the results from centrifuge tests investigating two of these methods are compared to results when no special foundation layout was implemented and the impact of the proposed isolation methods is discussed. © 2010 Taylor & Francis Group, London.

Centrifuge modeling of seismic liquefaction effects on adjacent shallow foundations

Shallow foundations built on saturated deposits of granular soils in seismically active areas are, regardless of their static bearing capacity, critical structures during seismic events. A single centrifuge experiment involving shallow foundations situated atop a liquefiable soil deposit has been performed to identify the mechanisms involved in the interaction between liquefaction-induced effects on neighboring shallow foundations. Centrifuge test results indicate that liquefaction causes significant settlements of footings, which are affected by the presence of neighboring foundations and can be extremely damaging to the superstructure. The understanding of these interaction effects is very important, mainly in densely populated urban areas. The development of high excess pore-pressures, localized drainage in response to the high transient hydraulic gradients, and earthquake-induced vertical motions to the footings are also important effects that are discussed to assist in enhancing current understanding and ability to predict liquefaction effects on shallow foundations. © 2014 Taylor & Francis Group.

TUNING, AVO, AND FLAT-SPOT EFFECTS IN A SEISMIC ANALYSIS OF NORTH SEA BLOCK F3

Reflection seismic data from the F3 block in the Dutch North Sea exhibits many large-amplitude reflections at shallow horizons, typically categorized as “brightspots ” (Schroot and Schuttenhelm, 2003), mainly because of their bright appearance. In most cases, these bright reflections show a significant “flatness” contrasting with local structural trends. While flatspots are often easily identified in thick reservoirs, we have often occasionally observed apparent flatspot tuning effects at fluid contacts near reservoir edges and in thin reservoir beds, while only poorly understanding them. We conclude that many of the shallow large-amplitude reflections in block F3 are dominated by flatspots, and we investigate the thin-bed tuning effects that such flatspots cause as they interact with the reflection from the reservoir’s upper boundary. There are two possible effects to be considered: (1) the “wedge-model” tuning effects of the flatspot and overlying brightspots, dimspots, or polarity-reversals; and (2) the stacking effects that result from possible inclusion of post-critical flatspot reflections in these shallow sands. We modeled the effects of these two phenomena for the particular stratigraphic sequence in block F3. Our results suggest that stacking of post-critical flatspot reflections can cause similar large-amplitude but flat reflections, in some cases even causing an interface expected to produce a ‘dimspot’ to appear as a ‘brightspot’. Analysis of NMO stretch and muting shows the likely exclusion of critical offset data in stacked output. If post-critical reflections are included in stacking, unusual results will be observed. In the North Sea case, we conclude the tuning effect was the primary reason causing for the brightness and flatness of these reflections. However, it is still important to note that care should be taken while applying muting on reflections with wide range of incidence angles and the inclusion of critical offset data may cause some spurious features in the stacked section.

Ocean bottom seismometer sgy-files of refraction seismic profiles from Professor Logachev cruise TTR16/3

Over the last decade pockmarks have proven to be important seabed features that provide information about fluid flow on continental margins. Their formation and dynamics are still poorly constrained due to the lack of proper three dimensional imaging of their internal structure. Numerous fluid escape features provide evidence for an active fluid-flow system on the Norwegian margin, specifically in the Nyegga region. In June-July 2006 a high-resolution seismic experiment using Ocean Bottom Seismometers (OBS) was carried out to investigate the detailed 3D structure of a pockmark named G11 in the region. An array of 14 OBS was deployed across the pockmark with 1 m location accuracy. Shots fired from surface towed mini GI guns were also recorded on a near surface hydrophone streamer. Several reflectors of high amplitude and reverse polarity are observed on the profiles indicating the presence of gas. Gas hydrates were recovered with gravity cores from less than a meter below the seafloor during the cruise. Indications of gas at shallow depths in the hydrate stability field show that methane is able to escape through the water-saturated sediments in the chimney without being entirely converted into gas hydrate. An initial 2D raytraced forward model of some of the P wave data along a line running NE-SW across the G11 pockmark shows, a gradual increase in velocity between the seafloor and a gas charged zone lying at ~300 m depth below the seabed. The traveltime fit is improved if the pockmark is underlain by velocities higher than in the surrounding layer corresponding to a pipe which ascends from the gas zone, to where it terminates in the pockmark as seen in the reflection profiles. This could be due to the presence of hydrates or carbonates within the sediments.