Local Earthquake Tomography by trans-dimensional Monte Carlo Sampling

In this study a new, fully non-linear, approach to Local Earthquake Tomography is presented. Local Earthquakes Tomography (LET) is a non-linear inversion problem that allows the joint determination of earthquakes parameters and velocity structure from arrival times of waves generated by local sources. Since the early developments of seismic tomography several inversion methods have been developed to solve this problem in a linearized way. In the framework of Monte Carlo sampling, we developed a new code based on the Reversible Jump Markov Chain Monte Carlo sampling method (Rj-McMc). It is a trans-dimensional approach in which the number of unknowns, and thus the model parameterization, is treated as one of the unknowns. I show that our new code allows overcoming major limitations of linearized tomography, opening a new perspective in seismic imaging. Synthetic tests demonstrate that our algorithm is able to produce a robust and reliable tomography without the need to make subjective a-priori assumptions about starting models and parameterization. Moreover it provides a more accurate estimate of uncertainties about the model parameters. Therefore, it is very suitable for investigating the velocity structure in regions that lack of accurate a-priori information. Synthetic tests also reveal that the lack of any regularization constraints allows extracting more information from the observed data and that the velocity structure can be detected also in regions where the density of rays is low and standard linearized codes fails. I also present high-resolution Vp and Vp/Vs models in two widespread investigated regions: the Parkfield segment of the San Andreas Fault (California, USA) and the area around the Alto Tiberina fault (Umbria-Marche, Italy). In both the cases, the models obtained with our code show a substantial improvement in the data fit, if compared with the models obtained from the same data set with the linearized inversion codes.

Thermal conductivity at DSDP Leg 60 Holes

A total of 1547 thermal conductivity values were determined by both the NP (needle probe method) and the QTM (quick thermal conductivity meter) on 1319 samples recovered during DSDP Leg 60. The NP method is primarily for the measurement of soft sedimentary samples, and the result is free from the effect of porewater evaporation. Measurement by the QTM method is faster and is applicable to all types of samples-namely, sediments (soft, semilithified, and lithified) and basement rocks. Data from the deep holes at Sites 453, 458, and 459 show that the thermal conductivity increases with depth, the rate of increase ranging from (0.18 mcal/cm s °C)/100 m at Site 459 to (0.72 mcal/cm s °C)/100 m at Site 456. A positive correlation between the sedimentary accumulation rate and the rate of thermal conductivity increase with depth indicates that both compaction and lithification are important factors. Drilled pillow basalts show nearly uniform thermal conductivity. At She 454 the thermal conductivity of one basaltic flow unit was higher near the center of the unit and lower toward the margin, reflecting variable vesicularity. Hydrothermally altered basalts at Site 456 showed higher thermal conductivity than fresh basalt because secondary calcite, quartz, and pyrite are generally more thermally conductive than fresh basalt. The average thermal conductivity in the top 50 meters of sediments correlates inversely with water depth because of dissolution of calcite, a mineral with high thermal conductivity, from the sediments as the water depth exceeds the lysocline and the carbonate compensation depth. Differences between the Mariana Trench data and the Mariana Basin and Trough data may reflect different abundances of terrigenous material in the sediment. There are remarkable correlations between thermal conductivity and other physical properties. The relationship between thermal conductivity and compressional wave velocity can be used to infer the ocean crustal thermal conductivity from the seismic velocity structure.

Horizontal to vertical spectral ratio measurements in Port-au-Prince (Haiti) area damaged by the 2010 Haiti earthquake

In order to evaluate ground shaking characteristics due to surface soil layers in the urban area of Port-au-Prince, short-period ambient noise observation has been performed approximately in a 500x500m grid. The HVSR method was applied to this set of 36 ambient noise measurement points to determine a distribution map of soil predominant periods. This map reveals a general increasing trend in the period values, from the Miocene conglomerates in the northern and southern parts of the town to the central and western zones formed of Pleistocene and Holocene alluvial deposits respectively, where the shallow geological materials that cover the basement increase in thickness. Shorter predominant periods (less than 0.3 s) were found in mountainous and neighbouring zones, where the thickness of sediments is smaller whereas longer periods (greater than 0.5 s) appear in Holocene alluvial fans, where the thickness of sediments is larger. The shallow shear-wave velocity structure have been estimated by means of inversion of Rayleigh wave dispersion data obtained from vertical-component array records of ambient noise. The measurements were carried out at one open space located in Holocene alluvial deposits, using 3 regular pentagonal arrays with 5, 10 and 20m respectively. Reliable dispersion curves were retrieved for frequencies between 4.0 and 14 Hz, with phase velocity values ranging from 420m/s down to 270 m/s. Finally, the average shear-wave velocity of the upper 30 m (VS30) was inverted for characterization of this geological unit.

(Table 1) Comparisons of dynamic and static Young's moduli of basalt at DSDP Holes 69-504B, 70-504B and 83-504B

Compressional- and shear-wave velocity logs (Vp and Vs, respectively) that were run to a sub-basement depth of 1013 m (1287.5 m sub-bottom) in Hole 504B suggest the presence of Layer 2A and document the presence of layers 2B and 2C on the Costa Rica Rift. Layer 2A extends from the mudline to 225 m sub-basement and is characterized by compressional-wave velocities of 4.0 km/s or less. Layer 2B extends from 225 to 900 m and may be divided into two intervals: an upper level from 225 to 600 m in which Vp decreases slowly from 5.0 to 4.8 km/s and a lower level from 600 to about 900 m in which Vp increases slowly to 6.0 km/s. In Layer 2C, which was logged for about 100 m to a depth of 1 km, Vp and Vs appear to be constant at 6.0 and 3.2 km/s, respectively. This velocity structure is consistent with, but more detailed than the structure determined by the oblique seismic experiment in the same hole. Since laboratory measurements of the compressional- and shear-wave velocity of samples from Hole 504B at Pconfining = Pdifferential average 6.0 and 3.2 km/s respectively, and show only slight increases with depth, we conclude that the velocity structure of Layer 2 is controlled almost entirely by variations in porosity and that the crack porosity of Layer 2C approaches zero. A comparison between the compressional-wave velocities determined by logging and the formation porosities calculated from the results of the large-scale resistivity experiment using Archie's Law suggest that the velocity- porosity relation derived by Hyndman et al. (1984) for laboratory samples serves as an upper bound for Vp, and the noninteractive relation derived by Toksöz et al. (1976) for cracks with an aspect ratio a = 1/32 serves as a lower bound.

Estrutura crustal e mantélica da província Borborema através de funções do receptor e dispersão de ondas superficiais

The Borborema Province, located in northeastern Brazil, has a basement of Precambrian age and a tectonic framework structured at the Neoproterozoic (740-560 Ma). After separation between South America and Africa during the Mesozoic, a rift system was formed, giving rise to a number of marginal and inland basins in the Province. After continental breakup, episodes of volcanism and uplift characterized the evolution of the Province. Plateau uplift was initially related to magmatic underplating of mafic material at the base of the crust, perhaps related to the generation of young continental plugs (45-7 Ma) along the Macau-Queimadas Alignment (MQA), due to a small-scale convection at the continental edge. The goal of this study is to investigate the causes of intra-plate uplift and its relationship to MQA volcanism, by using broadband seismology and integrating our results with independent geophysical and geological studies in the Borborema Province. The investigation of the deep structure of the Province with broadband seismic data includes receiver functions and surface-wave dispersion tomography. Both the receiver functions and surface-wave dispersion tomography are methods that use teleseismic events and allow to develop estimates of crustal parameters such as crustal thickness, Vp/Vs ratio, and S-velocity structure. The seismograms used for the receiver function work were obtained from 52 stations in Northeast Brazil: 16 broadband stations from the RSISNE network (Rede Sismográfica do Nordeste do Brasil), and 21 short-period and 6 broadband stations from the INCT-ET network (Instituto Nacional de Ciência e Tecnologia – Estudos Tectônicos). These results add signifi- cantly to previous datasets collected at individual stations in the Province, which include station RCBR (GSN - Global Seismic Network), stations CAUB and AGBL (Brazilian Lithosphere Seismic Project IAG/USP), and 6 other broadband stations that were part of the Projeto Milênio - Estudos geofísicos e tectônicos na Província Borborema/CNPq. For the surface-wave vii tomography, seismograms recorde at 22 broadband stations were utilized: 16 broadband stations from the RSISNE network and 6 broadband stations from the Milênio project. The new constraints developed in this work include: (i) estimates of crustal thickness and bulk Vp/Vs ratio for each station using receiver functions; (ii) new measurements of surfassewave group velocity, which were integrated to existing measurementes from a continental-scale tomography for South America, and (iii) S-wave velocity models (1D) at various locations in the Borborema Province, developed through the simultaneous inversion of receiver functions and surface-wave dispersion velocities. The results display S-wave velocity structure down to the base of the crust that are consistent with the presence of a 5-7.5 km thick mafic layer. The mafic layer was observed only in the southern portion of the Plateau and absent in its northern portion. Another important observation is that our models divide the plateau into a region of thin crust (northern Plateau) and a region of thick crust (southern Plateau), confirming results from independent refraction surveys and receiver function analyses. Existing models of plateau uplift, nonetheless, cannot explain all the new observations. It is proposed that during the Brazilian orogeny a layer of preexisting mafic material was delaminated, as a whole or in part, from the original Brasiliano crust. Partial delamination would have happened in the southern portion of the plateau, where independent studies found evidence of a more resistant rheology. During Mesozoic rifting, thinning of the crust around the southern Plateau would have formed the marginal basins and the Sertaneja depression, which would have included the northern part of the Plateau. In the Cenozoic, uplift of the northern Plateau would have occurred, resulting in a northern Plateau without mafic material at the base of the crust and a southern Plateau with partially delaminated mafic layer.

(Table 1) Mass loss of the Greenland Ice Sheet

Using satellite radar interferometry observations of Greenland, we detected widespread glacier acceleration below 66° north between 1996 and 2000, which rapidly expanded to 70° north in 2005. Accelerated ice discharge in the west and particularly in the east doubled the ice sheet mass deficit in the last decade from 90 to 220 cubic kilometers per year. As more glaciers accelerate farther north, the contribution of Greenland to sea-level rise will continue to increase.

Constraints on the permeability structure of alluvial aquifers from the poro-elastic inversion of multi-frequency P-wave velocity logs

We have explored the possibility of obtaining first-order permeability estimates for saturated alluvial sediments based on the poro-elastic interpretation of the P-wave velocity dispersion inferred from sonic logs. Modern sonic logging tools designed for environmental and engineering applications allow one for P-wave velocity measurements at multiple emitter frequencies over a bandwidth covering 5 to 10 octaves. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and typical emitter frequencies ranging from approximately 1 to 30 kHz, the observable velocity dispersion should be sufficiently pronounced to allow one for reliable first-order estimations of the permeability structure. The corresponding predictions have been tested on and verified for a borehole penetrating a typical surficial alluvial aquifer. In addition to multifrequency sonic logs, a comprehensive suite of nuclear and electrical logs, an S-wave log, a litholog, and a limited number laboratory measurements of the permeability from retrieved core material were also available. This complementary information was found to be essential for parameterizing the poro-elastic inversion procedure and for assessing the uncertainty and internal consistency of corresponding permeability estimates. Our results indicate that the thus obtained permeability estimates are largely consistent with those expected based on the corresponding granulometric characteristics, as well as with the available evidence form laboratory measurements. These findings are also consistent with evidence from ocean acoustics, which indicate that, over a frequency range of several orders-of-magnitude, the classical theory of poro-elasticity is generally capable of explaining the observed P-wave velocity dispersion in medium- to fine-grained seabed sediments

Structure and dynamics of the Shapley Supercluster - Velocity catalogue, general morphology and mass

We present results of our wide-field redshift survey of galaxies in a 285 square degree region of the Shapley Supercluster (SSC), based on a set of 10 529 velocity measurements (including 1201 new ones) on 8632 galaxies obtained from various telescopes and from the literature. Our data reveal that the main plane of the SSC (v approximate to 14 500 km s(-1)) extends further than previous estimates, filling the whole extent of our survey region of 12 degrees by 30 degrees on the sky (30 x 75 h(-1) Mpc). There is also a connecting structure associated with the slightly nearer Abell 3571 cluster complex (v approximate to 12 000 km s(-1)). These galaxies seem to link two previously identified sheets of galaxies and establish a connection with a third one at v = 15 000 km s(-1) near RA = 13(h). They also tend to fill the gap of galaxies between the foreground Hydra-Centaurus region and the more distant SSC. In the velocity range of the Shapley Supercluster (9000 km s(-1) < cz < 18 000 km s(-1)), we found redshift-space overdensities with b(j) < 17.5 of similar or equal to 5.4 over the 225 square degree central region and similar or equal to 3.8 in a 192 square degree region excluding rich clusters. Over the large region of our survey, we find that the intercluster galaxies make up 48 per cent of the observed galaxies in the SSC region and, accounting for the different completeness, may contribute nearly twice as much mass as the cluster galaxies. In this paper, we discuss the completeness of the velocity catalogue, the morphology of the supercluster, the global overdensity, and some properties of the individual galaxy clusters in the Supercluster.

Which-path problem for one and two particles with two degrees of freedom and a relation between transverse spatial structure and group velocity of light

Quantum mechanics, optics and indeed any wave theory exhibits the phenomenon of interference. In this thesis we present two problems investigating interference due to indistinguishable alternatives and a mostly unrelated investigation into the free space propagation speed of light pulses in particular spatial modes. In chapter 1 we introduce the basic properties of the electromagnetic field needed for the subsequent chapters. In chapter 2 we review the properties of interference using the beam splitter and the Mach-Zehnder interferometer. In particular we review what happens when one of the paths of the interferometer is marked in some way so that the particle having traversed it contains information as to which path it went down (to be followed up in chapter 3) and we review Hong-Ou-Mandel interference at a beam splitter (to be followed up in chapter 5). In chapter 3 we present the first of the interference problems. This consists of a nested Mach-Zehnder interferometer in which each of the free space propagation segments are weakly marked by mirrors vibrating at different frequencies [1]. The original experiment drew the conclusions that the photons followed disconnected paths. We partition the description of the light in the interferometer according to the number of paths it contains which-way information about and reinterpret the results reported in [1] in terms of the interference of paths spatially connected from source to detector. In chapter 4 we briefly review optical angular momentum, entanglement and spontaneous parametric down conversion. These concepts feed into chapter 5 in which we present the second of the interference problems namely Hong-Ou-Mandel interference with particles possessing two degrees of freedom. We analyse the problem in terms of exchange symmetry for both boson and fermion pairs and show that the particle statistics at a beam splitter can be controlled for suitably chosen states. We propose an experimental test of these ideas using orbital angular momentum entangled photons. In chapter 6 we look at the effect that the transverse spatial structure of the mode that a pulse of light is excited in has on its group velocity. We show that the resulting group velocity is slower than the speed of light in vacuum for plane waves and that this reduction in the group velocity is related to the spread in the wave vectors required to create the transverse spatial structure. We present experimental results of the measurement of this slowing down using Hong-Ou-Mandel interference.

The structure and dynamics of Abell 1942

Aims. We present a dynamical analysis of the galaxy cluster Abell 1942 based on a set of 128 velocities obtained at the European Southern Observatory. Methods. Data on individual galaxies are presented and the accuracy of the determined velocities as some properties of the cluster are discussed. We have also made use of publicly available Chandra X-ray data. Results. We obtained an improved mean redshift value z = 0.22513 +/- 0.0008 and velocity dispersion sigma = 908(139)(+147) km s(-1). Our analysis indicates that inside a radius of similar to 1.5 h(70)(-1) Mpc (similar to 7 arcmin) the cluster is well relaxed, without any remarkable features and the X-ray emission traces the galaxy distribution fairly well. Two possible optical substructures are seen at similar to 5 arcmin from the centre in the northwest and the southwest directions, but are not confirmed by the velocity field. These clumps are, however, kinematically bound to the main structure of Abell 1942. X-ray spectroscopic analysis of Chandra data resulted in a temperature kT = 5.5+/-0.5 keV and metal abundance Z = 0.33 +/- 0.15 Z(circle dot). The velocity dispersion corresponding to this temperature using the T(X-sigma) scaling relation is in good agreement with the measured galaxy velocities. Our photometric redshift analysis suggests that the weak lensing signal observed to the south of the cluster and previously attributed to a ""dark clump"" is produced by background sources, possibly distributed as a filamentary structure.

Solution structure of the human signaling protein RACK1

Background: The adaptor protein RACK1 (receptor of activated kinase 1) was originally identified as an anchoring protein for protein kinase C. RACK1 is a 36 kDa protein, and is composed of seven WD repeats which mediate its protein-protein interactions. RACK1 is ubiquitously expressed and has been implicated in diverse cellular processes involving: protein translation regulation, neuropathological processes, cellular stress, and tissue development. Results: In this study we performed a biophysical analysis of human RACK1 with the aim of obtaining low resolution structural information. Small angle X-ray scattering (SAXS) experiments demonstrated that human RACK1 is globular and monomeric in solution and its low resolution structure is strikingly similar to that of an homology model previously calculated by us and to the crystallographic structure of RACK1 isoform A from Arabidopsis thaliana. Both sedimentation velocity and sedimentation equilibrium analytical ultracentrifugation techniques showed that RACK1 is predominantly a monomer of around 37 kDa in solution, but also presents small amounts of oligomeric species. Moreover, hydrodynamic data suggested that RACK1 has a slightly asymmetric shape. The interaction of RACK1 and Ki1/57 was tested by sedimentation equilibrium. The results suggested that the association between RACK1 and Ki-1/57(122-413) follows a stoichiometry of 1:1. The binding constant (KB) observed for RACK1-Ki-1/57(122-413) interaction was of around (1.5 +/- 0.2) x 10(6) M(-1) and resulted in a dissociation constant (KD) of (0.7 +/- 0.1) x 10(-6) M. Moreover, the fluorescence data also suggests that the interaction may occur in a cooperative fashion. Conclusion: Our SAXS and analytical ultracentrifugation experiments indicated that RACK1 is predominantly a monomer in solution. RACK1 and Ki-1/57(122-413) interact strongly under the tested conditions.