The use of shallow seismic reflection technique in near surface exploration of urban sites: An evaluation in the city of SãO Paulo, Brazil

In order to evaluate the use of shallow seismic technique to delineate geological and geotechnical features up to 40 meters depth in noisy urban areas covered with asphalt pavement, five survey lines were conducted in the metropolitan area of São Paulo City. The data were acquired using a 24-bit, 24-channel seismograph, 30 and 100 Hz geophones and a sledgehammer-plate system as seismic source. Seismic reflection data were recorded using a CMP (common mid point) acquisition method. The processing routine consisted of: prestack band-pass filtering (90-250 Hz); automatic gain control (AGC); muting (digital zeroin) of dead/noisy traces, ground roll, air-wave and refracted-wave; CMP sorting; velocity analyses; normal move-out corrections; residual static corrections; f-k filtering; CMP stacking. The near surface is geologically characterized by unconsolidated fill materials and Quaternary sediments with organic material overlying Tertiary sediments with the water table 2 to 5 m below the surface. The basement is composed of granite and gneiss. Reflections were observed from 40 milliseconds to 65 ms two-way traveltime and were related to the silt clay layer and fine sand layer contact of the Tertiary sediments and to the weathered basement. The CMP seismic-reflection technique has been shown to be useful for mapping the sedimentary layers and the bedrock of the São Paulo sedimentary basin for the purposes of shallow investigations related to engineering problems. In spite of the strong cultural noise observed in these urban areas and problems with planting geophones we verified that, with the proper equipment, field parameters and particularly great care in data collection and processing, we can overcome the adverse field conditions and to image reflections from layers as shallow as 20 meters.

Developing a system for down-hole seismic testing together with the CPTU

This paper presents a system for performing down-hole seismic test together with the piezocone test in order to determine the shear wave velocity (Vs) and for calculating the maximum shear modulus (Go); a basic parameter for analyzing the dynamic soil behavior and a reference value of the soil stiffness. The system components are described and tests results for checking the geophone response are also presented, both before and after installation into the probe. The system was used in down-hole tests carried out at three experimental research sites located in the interior of Sao Paulo State, Brazil, where in situ seismic test results are available. The Vs values measured in down-hole tests carried out with this system were consistent with those determined in cross-hole tests and with a commercial seismic piezocone, which enabled to validate the developed system.

The seismic SPT to determine the maximum shear modulus

This paper presents the first results of a system to carry out seismic SPT, which associates the up-hole seismic technique to the currently used SPT test. This hybrid test allows determining the maximum shear modulus (G 0) together with the SPT test. G 0 is the most important dynamic soil parameter and it is calculated based on the shear wave velocity (V S) measurements. The equipment and the method of analysis will be briefly described. The V s values measured using the S-SPT in one experimental research site located inland of Sao Paulo State-Brazil are presented and compared with down-hole and Seismic CPT test data. These preliminary results indicate that the S-SPT was satisfactory and can be used for dynamic site characterization.

CRS seismic data imaging system: a case study for basib reevaluation

Este trabalho apresenta resultados práticos de uma atenção sistemática dada ao processamento e à interpretação sísmica de algumas linhas terrestres do conjunto de dados do gráben do Tacutu (Brasil), sobre os quais foram aplicadas etapas fundamentais do sistema WIT de imageamento do empilhamento CRS (Superfície de Reflexão Comum) vinculado a dados. Como resultado, esperamos estabelecer um fluxograma para a reavaliação sísmica de bacias sedimentares. Fundamentado nos atributos de frente de onda resultantes do empilhamento CRS, um macro-modelo suave de velocidades foi obtido através de inversão tomográfica. Usando este macro-modelo, foi realizado uma migração à profundidade pré- e pós-empilhamento. Além disso, outras técnicas baseadas no empilhamento CRS foram realizadas em paralelo como correção estática residual e migração de abertura-limitada baseada na zona de Fresnel projetada. Uma interpretação geológica sobre as seções empilhadas e migradas foi esboçada. A partir dos detalhes visuais dos painéis é possível interpretar desconformidades, afinamentos, um anticlinal principal falhado com conjuntos de horstes e grábens. Também, uma parte da linha selecionada precisa de processamento mais detalhado para evidenciar melhor qualquer estrutura presente na subsuperfície.

Multiresolution spherical wavelet analysis in global seismic tomography

Every seismic event produces seismic waves which travel throughout the Earth. Seismology is the science of interpreting measurements to derive information about the structure of the Earth. Seismic tomography is the most powerful tool for determination of 3D structure of deep Earth's interiors. Tomographic models obtained at the global and regional scales are an underlying tool for determination of geodynamical state of the Earth, showing evident correlation with other geophysical and geological characteristics. The global tomographic images of the Earth can be written as a linear combinations of basis functions from a specifically chosen set, defining the model parameterization. A number of different parameterizations are commonly seen in literature: seismic velocities in the Earth have been expressed, for example, as combinations of spherical harmonics or by means of the simpler characteristic functions of discrete cells. With this work we are interested to focus our attention on this aspect, evaluating a new type of parameterization, performed by means of wavelet functions. It is known from the classical Fourier theory that a signal can be expressed as the sum of a, possibly infinite, series of sines and cosines. This sum is often referred as a Fourier expansion. The big disadvantage of a Fourier expansion is that it has only frequency resolution and no time resolution. The Wavelet Analysis (or Wavelet Transform) is probably the most recent solution to overcome the shortcomings of Fourier analysis. The fundamental idea behind this innovative analysis is to study signal according to scale. Wavelets, in fact, are mathematical functions that cut up data into different frequency components, and then study each component with resolution matched to its scale, so they are especially useful in the analysis of non stationary process that contains multi-scale features, discontinuities and sharp strike. Wavelets are essentially used in two ways when they are applied in geophysical process or signals studies: 1) as a basis for representation or characterization of process; 2) as an integration kernel for analysis to extract information about the process. These two types of applications of wavelets in geophysical field, are object of study of this work. At the beginning we use the wavelets as basis to represent and resolve the Tomographic Inverse Problem. After a briefly introduction to seismic tomography theory, we assess the power of wavelet analysis in the representation of two different type of synthetic models; then we apply it to real data, obtaining surface wave phase velocity maps and evaluating its abilities by means of comparison with an other type of parametrization (i.e., block parametrization). For the second type of wavelet application we analyze the ability of Continuous Wavelet Transform in the spectral analysis, starting again with some synthetic tests to evaluate its sensibility and capability and then apply the same analysis to real data to obtain Local Correlation Maps between different model at same depth or between different profiles of the same model.

Three dimensional seismic imaging and earthquake locations in a complex, segmented fault region in Southern Apennines (Italy)

The southern Apennines of Italy have been experienced several destructive earthquakes both in historic and recent times. The present day seismicity, characterized by small-to-moderate magnitude earthquakes, was used like a probe to obatin a deeper knowledge of the fault structures where the largest earthquakes occurred in the past. With the aim to infer a three dimensional seismic image both the problem of data quality and the selection of a reliable and robust tomographic inversion strategy have been faced. The data quality has been obtained to develop optimized procedures for the measurements of P- and S-wave arrival times, through the use of polarization filtering and to the application of a refined re-picking technique based on cross-correlation of waveforms. A technique of iterative tomographic inversion, linearized, damped combined with a strategy of multiscale inversion type has been adopted. The retrieved P-wave velocity model indicates the presence of a strong velocity variation along a direction orthogonal to the Apenninic chain. This variation defines two domains which are characterized by a relatively low and high velocity values. From the comparison between the inferred P-wave velocity model with a portion of a structural section available in literature, the high velocity body was correlated with the Apulia carbonatic platforms whereas the low velocity bodies was associated to the basinal deposits. The deduced Vp/Vs ratio shows that the ratio is lower than 1.8 in the shallower part of the model, while for depths ranging between 5 km and 12 km the ratio increases up to 2.1 in correspondence to the area of higher seismicity. This confirms that areas characterized by higher values are more prone to generate earthquakes as a response to the presence of fluids and higher pore-pressures.

Deep geometry of subduction below the Andean belt of Colombia as revealed by seismic tomography

In this study new tomographic models of Colombia were calculated. I used the seismicity recorded by the Colombian seismic network during the period 2006-2009. In this time period, the improvement of the seismic network yields more stable hypocentral results with respect to older data set and allows to compute new 3D Vp and Vp/Vs models. The final dataset consists of 10813 P- and 8614 S-arrival times associated to 1405 earthquakes. Tests with synthetic data and resolution analysis indicate that velocity models are well constrained in central, western and southwestern Colombia to a depth of 160 km; the resolution is poor in the northern Colombia and close to Venezuela due to a lack of seismic stations and seismicity. The tomographic models and the relocated seismicity indicate the existence of E-SE subducting Nazca lithosphere beneath central and southern Colombia. The North-South changes in Wadati-Benioff zone, Vp & Vp/Vs pattern and volcanism, show that the downgoing plate is segmented by slab tears E-W directed, suggesting the presence of three sectors. Earthquakes in the northernmost sector represent most of the Colombian seimicity and concentrated on 100-170 km depth interval, beneath the Eastern Cordillera. Here a massive dehydration is inferred, resulting from a delay in the eclogitization of a thickened oceanic crust in a flat-subduction geometry. In this sector a cluster of intermediate-depth seismicity (Bucaramanga Nest) is present beneath the elbow of the Eastern Cordillera, interpreted as the result of massive and highly localized dehydration phenomenon caused by a hyper-hydrous oceanic crust. The central and southern sectors, although different in Vp pattern show, conversely, a continuous, steep and more homogeneous Wadati-Benioff zone with overlying volcanic areas. Here a "normalthickened" oceanic crust is inferred, allowing for a gradual and continuous metamorphic reactions to take place with depth, enabling the fluid migration towards the mantle wedge.

Seismic sequences analysis for estimation of earthquake source parameters: corner frequency, stress drop, and seismic moment observations

The present study has been carried out with the following objectives: i) To investigate the attributes of source parameters of local and regional earthquakes; ii) To estimate, as accurately as possible, M0, fc, Δσ and their standard errors to infer their relationship with source size; iii) To quantify high-frequency earthquake ground motion and to study the source scaling. This work is based on observational data of micro, small and moderate -earthquakes for three selected seismic sequences, namely Parkfield (CA, USA), Maule (Chile) and Ferrara (Italy). For the Parkfield seismic sequence (CA), a data set of 757 (42 clusters) repeating micro-earthquakes (0 ≤ MW ≤ 2), collected using borehole High Resolution Seismic Network (HRSN), have been analyzed and interpreted. We used the coda methodology to compute spectral ratios to obtain accurate values of fc , Δσ, and M0 for three target clusters (San Francisco, Los Angeles, and Hawaii) of our data. We also performed a general regression on peak ground velocities to obtain reliable seismic spectra of all earthquakes. For the Maule seismic sequence, a data set of 172 aftershocks of the 2010 MW 8.8 earthquake (3.7 ≤ MW ≤ 6.2), recorded by more than 100 temporary broadband stations, have been analyzed and interpreted to quantify high-frequency earthquake ground motion in this subduction zone. We completely calibrated the excitation and attenuation of the ground motion in Central Chile. For the Ferrara sequence, we calculated moment tensor solutions for 20 events from MW 5.63 (the largest main event occurred on May 20 2012), down to MW 3.2 by a 1-D velocity model for the crust beneath the Pianura Padana, using all the geophysical and geological information available for the area. The PADANIA model allowed a numerical study on the characteristics of the ground motion in the thick sediments of the flood plain.

Seismic tomographic full-waveform inversion for the Vrancea sinking lithosphere structure using the adjoint method.

The Vrancea region, at the south-eastern bend of the Carpathian Mountains in Romania, represents one of the most puzzling seismically active zones of Europe. Beside some shallow seismicity spread across the whole Romanian territory, Vrancea is the place of an intense seismicity with the presence of a cluster of intermediate-depth foci placed in a narrow nearly vertical volume. Although large-scale mantle seismic tomographic studies have revealed the presence of a narrow, almost vertical, high-velocity body in the upper mantle, the nature and the geodynamic of this deep intra-continental seismicity is still questioned. High-resolution seismic tomography could help to reveal more details in the subcrustal structure of Vrancea. Recent developments in computational seismology as well as the availability of parallel computing now allow to potentially retrieve more information out of seismic waveforms and to reach such high-resolution models. This study was aimed to evaluate the application of a full waveform inversion tomography at regional scale for the Vrancea lithosphere using data from the 1999 six months temporary local network CALIXTO. Starting from a detailed 3D Vp, Vs and density model, built on classical travel-time tomography together with gravity data, I evaluated the improvements obtained with the full waveform inversion approach. The latter proved to be highly problem dependent and highly computational expensive. The model retrieved after the first two iterations does not show large variations with respect to the initial model but remains in agreement with previous tomographic models. It presents a well-defined downgoing slab shape high velocity anomaly, composed of a N-S horizontal anomaly in the depths between 40 and 70km linked to a nearly vertical NE-SW anomaly from 70 to 180km.

Seismic properties of sheeted dikes from ODP Hole 111-504B

Seismic velocities have been measured at confining pressures to 600 MPa for eight samples of sheeted dike rock obtained from Hole 504B during Leg ODP 111. The compressional- and shear-wave velocities are, in general, higher than the velocities measured in overlying dike rocks obtained from the hole during DSDP Leg 83. The velocity gradients observed in Layer 2C result from decreasing porosity with depth and increasing metamorphic grade. The laboratory-measured velocities of the Leg 111 dike rocks are similar to those of dike rocks reported for the Bay of Islands, Samail, and Troodos ophiolites.

Seismic properties of volcanic rocks from DSDP Hole 62-465A

Hess Rise, a major structural feature in the northern Pacific, is one of several oceanic plateaus which apparently possess anomalous seismic structures (Sutton et al., 1971). Hence, Laboratory measurements of compressional- and shear-wave velocities in rocks from oceanic plateau regions are of considerable interest. Several questions come to mind: (1) Are compressional - wave velocities of volcanic rocks from oceanic plateaus similar to basalts of equivalent age from normal oceanic crust? (2) Do velocity-density relations for plateau rocks fit the well-established trends for Layer 2 basalts? (3) How do Poisson's ratios, determined from compressional- and shear-wave velocities, of oceanic plateau rocks compare with those of normal sea-floor basalts? To answer these questions, we have selected for velocity measurements five volcanic samples from different depths from Hole 465A, on southern Hess Rise. It is particularly significant that all the rocks are vesicular and have been highly altered, which significantly affects their elastic properties.

Seismic investigations of the Greenland ice sheet at EGIG points

The compressional and the shear wave velocities in the Greenland ice sheet are derived from seismic records of the EGIG 1959. Further the variation of velocities in the firn and the dependance of Poisson's ratio from depth are determined. At Station Centrale, two P-waves are recorded from underground layers. Their velocities show that the ice basement consists of crystalline rocks. The P-wave velocities derived from reflections agree well with those obtained by refraction shooting. From this agreement results that the ice is ± homogenous and ± isotropic for Pwaves. The elastic constants for isotropic ice are calculated. Finally the temperature dependence of the velocities is discussed.

Multibeam and seismic survey, physical oceanography, and sea-bottom videos from the Bay of Biscay

We report the northernmost and deepest known occurrence of deep-water pycnodontine oysters, based on two surveys along the French Atlantic continental margin to the La Chapelle continental slope (2006) and the Guilvinec Canyon (2008). The combined use of multibeam bathymetry, seismic profiling, CTD casts and a remotely operated vehicle (ROV) made it possible to describe the physical habitat and to assess the oceanographic control for the recently described species Neopycnodonte zibrowii. These oysters have been observed in vivo in depths from 540 to 846 m, colonizing overhanging banks or escarpments protruding from steep canyon flanks. Especially in the Bay of Biscay, such physical habitats may only be observed within canyons, where they are created by both long-term turbiditic and contouritic processes. Frequent observations of sand ripples on the seabed indicate the presence of a steady, but enhanced bottom current of about 40 cm/s. The occurrence of oysters also coincides with the interface between the Eastern North Atlantic Water and the Mediterranean Outflow Water. A combination of this water mass mixing, internal tide generation and a strong primary surface productivity may generate an enhanced nutrient flux, which is funnelled through the canyon. When the ideal environmental conditions are met, up to 100 individuals per m² may be observed. These deep-water oysters require a vertical habitat, which is often incompatible with the requirements of other sessile organisms, and are only sparsely distributed along the continental margins. The discovery of these giant oyster banks illustrates the rich biodiversity of deep-sea canyons and their underestimation as true ecosystem hotspots.

(Table 1) Bulk and grain density, porosity, and p- and s-wave velocity characteristics of ODP Hole 163-990A sediments

Seismic velocities in rocks are influenced by the properties of the solid, the pore fluid, and the pore space. Cracks dramatically affect seismic velocities in rocks; their influence on the effective elastic moduli of rocks depends on their shape and concentration. Thin cracks (or fractures) substantially lower the moduli of a rock relative to the effect of spherical voids (or vesicles), and lower moduli are reflected by lower P- and S-wave velocities. The objective of this research is to determine the types and concentrations of cracks and their influence on the seismic properties of subaerially erupted basalts drilled from Hole 990A on the Southeast Greenland margin during Ocean Drilling Program Leg 163. Ellipsoidal cracks are used to model the voids in the rocks. The elastic moduli of the solid (grains) are also free parameters in the inverse modeling procedure. The apparent grain moduli reflect a weighted average of the moduli of the constituent minerals (e.g., plagioclase, augite, and clay minerals). The results indicate that (1) there is a strong relationship between P-wave velocity and porosity, suggesting a similarity of pore shape distributions, (2) the distribution of crack types within the massive, central region of aa flows from Hole 990A is independent of total porosity, (3) thin cracks are the first to be effectively sealed by alteration products, and (4) grain densities (an alteration index) and apparent grain moduli of the basalt samples are directly related.