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.

Velocity-density systematics for basalts at DSPD Leg 54 Holes

Precise velocity and density measurements at atmospheric and elevated pressures have been obtained on basalts drilled in 13 holes during Leg 54. The vp-sigma data show systematics which are controlled by the degree of crystallinity (or conversely, glassiness), microstructure, and original mineralogy and chemistry. Extensively fractionated basalts with marked iron enrichment produce anomalously low measured velocities at densities above 2.90 g/cm**3. Also, the effective in situ pressure acting on Leg 54 basalts is less than hydrostatic, and perhaps close to zero. At zero effective pressure, the measured velocities average 2.5 km/s higher than East Pacific Rise upper crustal velocities determined by seismic refraction. This implies that the in situ velocities are undoubtedly a result of the highly fragmented nature of East Pacific Rise crust.

Seismic velocities and elastic properties of gabbroic rocks from ODP Hole 118-735B

The nearly continuous recovery of 0.5 km of generally fresh, layer 3 gabbroic rocks at Hole 735B, especially near the bottom of the section, presents scientists an unusual opportunity to study the detailed elastic properties of the lower oceanic crust. Extending compressional-wave and density shipboard measurements at room pressure, Vp and Vs were measured at pressures from 20 to 200 MPa using the pulse transmission method. All of the rocks exhibit significant increases in velocity with increasing pressure up to about 150 MPa, a feature attributed to the closing of microcrack porosity. Measured velocities reflect the mineralogical makeup and microstructures acquired during the tectonic history of Hole 735B. Most of the undeformed and unaltered gabbros are approximately 65:35 plagioclase/clinopyroxene rocks plus olivine or oxide minerals, and the observed densities and velocities are fully consistent with the Voigt-Reuss-Hill (VRH) averages of the component minerals and their proportions. Depending on their olivine content, the predominant olivine gabbros at 200 MPa have average Vp = 7.1 ± 0.2 km/s, Vs = 3.9 ± 0.1 km/s, and grain densities of 2.95 ± 0.5 g/cm3. The less abundant iron-titanium (Fe-Ti) oxide gabbros average Vp = 6.75 ± 0.15 km/s, Vs = 3.70 ± 0.1 km/s, and grain densities of 3.22 ± 0.05 g/cm3, reflecting the higher densities and lower velocities of oxide minerals compared to olivine. About 30% of the core is plastically deformed, and the densities and directionally averaged velocities of these shear-zone tectonites are generally consistent with those of the gabbros, their protoliths. Three sets of observations indicate that the shear-zone metagabbros are elastically anisotropic: (1) directional variations in Vp, both vertical and horizontal and with respect to foliation and lineation; (2) discrepancies among Vp values for the horizontal cores and the VRH averages of the component minerals and their mineral proportions, suggesting preferred crystallographic orientations of anisotropic minerals; and (3) variations of Vs of up to 7%, with polarization directions parallel and perpendicular to foliation. Optical inspection of thin sections of the same samples indicates that plagioclase feldspar, clinopyroxene, and amphibole typically display crystallographic-preferred orientations, and this, plus the elastic anisotropy of these minerals, suggests that preferred orientations are responsible for much of the observed anisotropy, particularly at high pressure. Alteration tends to be localized to brittle faults and brecciated zones, and typical alteration minerals are amphibole and secondary plagioclase, which do not significantly change the velocity-density relationships.

P-wave velocities and applied bed thickness and velocity changes for synthetic seismograms, ODP Leg 188 and Leg 119 sites

Synthetic seismograms provide a crucial link between lithologic variations within a drill hole and reflectors on seismic profiles crossing the site. In essence, they provide a ground-truth for the interpretation of seismic data. Using a combination of core and logging data, we created synthetic seismograms for Ocean Drilling Program Sites 1165 and 1166, drilled during Leg 188, and Site 742, drilled during Leg 119, all in Prydz Bay, Antarctica. Results from Site 1165 suggest that coring penetrated a target reflector initially thought to represent the onset of drift sedimentation, but the lithologic change across the boundary does not show a change from predrift to drift sediments. The origin of a shallow reflector packet in the seismic line across Site 1166 and a line connecting Sites 1166 and 742 was resolved into its constituent sources, as this reflector occurs in a region of large-scale, narrowly spaced impedance changes. Furthermore, Site 1166 was situated in a fluvio-deltaic system with widely variable geology, and bed thickness changes were estimated between the site and both seismic lines.

Physical properties, velocity and attenuation characteristics, and mineralogy of gabbros from ODP Hole 176-735B

Laboratory compressional wave (Vp) and shear wave (Vs) velocities were measured as a function of confining pressure for the gabbros from Hole 735B and compared to results from Leg 118. The upper 500 m of the hole has a Vp mean value of 6895 m/s measured at 200 MPa, and at 500 meters below seafloor (mbsf), Vp measurements show a mean value of 7036 m/s. Vs mean values in the same intervals are 3840 m/s and 3857 m/s, respectively. The mean Vp and Vs values obtained from log data in the upper 600 m are 6520 and 3518 m/s, respectively. These results show a general increase in velocity with depth and the velocity gradients estimate an upper mantle depth of 3.32 km. This value agrees with previous work based on dredged samples and inversion of rare element concentrations in basalts dredged from the conjugate site to the north of the Atlantis Bank. Laboratory measurements show Vp anisotropy ranging between 0.4% and 8.8%, with the majority of the samples having values less than 3.8%. Measurements of velocity anisotropy seem to be associated with zones of high crystal-plastic deformation with predominant preferred mineral orientations of plagioclase, amphiboles, and pyroxenes. These findings are consistent with results on gabbros from the Hess Deep area and suggest that plastic deformation may play an important role in the seismic properties of the lower oceanic crust. In contrast to ophiolite studies, many of the olivine gabbros show a small degree of anisotropy. Log derived Vs anisotropy shows an average of 5.8% for the upper 600 m of Hole 735B and tends to decrease with depth where the overburden pressure and the age of the crustal section suggests closure of cracks and infilling of fractures by alteration minerals. Overall the results indicate that the average shear wave splitting in Hole 735B might be influenced by preferred structural orientations and the average value of shear wave splitting may not be a maximum because structural dips are <90°. The maximum fast-wave orientation values could be influenced by structural features striking slightly oblique to this orientation or by near-field stress concentrations. However, flexural wave dispersion analyses have not been performed to confirm this hypothesis or to indicate to what extent the near-field stresses may be influencing shear wave propagation. Acoustic impedance contrasts calculated from laboratory and logging data were used to generate synthetic seismograms that aid in the interpretation of reflection profiles. Several prominent reflections produced by these calculations suggest that Fe-Ti oxides and shear zones may contribute to the reflective nature of the lower oceanic crust. Laboratory velocity attenuation (Q) measurements from below 500 m have a mean value of 35.1, which is consistent with previous vertical seismic profile (VSP) and laboratory measurements on the upper 500 m.