Tomographie crustale des Pyrénées et des régions avoisinantes par corrélation de bruit

In this thesis, we applied the ambient noise correlation method in the very heterogeneous region of the Pyrenees and the surrounding areas (mountain belt and thick sedimentary basins). The dataset used is a combination of two temporary broadband arrays from France and Spain (PYROPE and IBERARRAY) and stations of the French and Catalan permanent broadband arrays. Seismic noise recorded over years by the 158 stations was used to calculate correlations in a period range of 5-55 s. Observed Rayleigh and Love wave group velocities between pairs of stations were used as input to a linearized inversion scheme, where we obtained for each period group velocities maps, with a lateral resolution of approximately 40 km. The comparison between the two type of waves demonstrates radial anisotropy at short periods, while little or no radial anisotropy is present at long periods. We developed a new strategy for the inversion of dispersion curves to shear wave velocity models and applied it to the Rayleigh waves group velocity. This approach is based on the combination of a full exploration of the model space and a linearized inversion. The obtained model, validated by the comparison of our results with the results of other methods, is the first complete 3-D crustal Vs model of the region. We in particular note : (1) Atypical S-waves profiles in the Eastern part of the Massif Central which indicate a thinned crust and low velocities in the uppermost mantle. (2) High-velocity anomalies at 25 km depth beneath the Labourd-Mauléon area and the on-land continuation of the Parentis basin. We suggest that they are the traces of the hyper-extension which might have preceded the collision phase which lead to the formation of the Pyrenees. The strong heterogeneity of our study region is also well-adapted to analyze the influence of the ray deviations on the reconstruction of the model. First results show that the reconstructed model, using great-circle paths, does not explain the ray deviations as observed by beamforming. Observed deviations therefore carry the potential of improving the model in a combined inversion scheme.Durant cette thèse, nous appliquons la méthode de la corrélation de bruit ambiant dans la région particulièrement hétérogène des Pyrénées et des ses alentours (socle rocheux affleurant et bassin sédimentaires épais). Le jeu de données utilisé est une combinaison de deux réseaux temporaires large-bande français et espagnol (PYROPE et IBERARRAY) et de stations des réseaux permanents large-bande français et catalan. Le bruit sismique, enregistré pendant un an par les 158 stations est utilisé pour calculer les corrélations dans la gamme de période 5-55 s. Les vitesses de groupe de l'onde de Rayleigh et de l'onde de Love entre paires de stations sont inversées de manière linéarisée et nous obtenons, pour chaque période, des cartes de vitesses de groupe avec une résolution latérale d'environ 40 km. La comparaison entre les deux types d'ondes montre qu'il existe une anisotropie radiale à courtes périodes, alors que peu ou pas d'anisotropie radiale n'est visible aux périodes plus longues. Nous avons développé une nouvelle stratégie d'inversion des courbes de dispersion en modèle de vitesse d'onde S que nous avons appliqué sur les vitesses de groupe de l'onde de Rayleigh. Cette approche d'inversion est basée sur une exploration complète de l'espace des modèles et une inversion linéarisée. Le modèle obtenu, validé par la comparaison avec des résultats provenant d'autres méthodes, est le premier modèle 3-D crustale en vitesse d'onde S de la région et il permet d'apporter des contraintes sur la géodynamique des Pyrénées et de ses alentours. Deux points importants sont soulevés : (1) Des profils de vitesses atypiques sous l'Est du Massif Central, avec une croûte amincie et des vitesses anormalement faibles dans le manteau supérieur. (2) Deux anomalies de vitesses rapides sous la zone du Labourd-Mauléon et dans le prolongement du bassin de Parentis. Ces anomalies, situées à 25 km de profondeur, sont interprétées comme les traces de l'hyper-extension qui aurait précédée la phase de collision amenant à la formation des Pyrénées. La forte hétérogénéité de la zone permet également de faire l'analyse de l'influence de la non prise en compte des déviations de rais lors de l'inversion. Les premiers résultats montrent que le modèle obtenu en utilisant la théorie des rais droits ne permet pas d'expliquer les déviations calculées par la méthode de formation de voie. Ces déviations observées peuvent donc apporter une amélioration du modèle en considérant l'utilisation d'une inversion combinée

ERT and crosswell EM imaging of CO2: Examples from a shallow injection experiment at the Carbon Management Canada CaMI FRS in Southeast Alberta, Canada

Electromagnetic (EM) geophysical techniques offer the possibility of monitoring subsurface CO2 in saline reservoirs due to the fact that CO2 has high electrical resistivity compared to the surrounding geologic materials. In this paper we first discuss the underlying physics of two different borehole-based EM monitoring techniques; electrical resistivity tomography (ERT) and crosswell EM. This discussion is followed by the description of an experiment at the Carbon Management Canada CaMI FRS test site where time-lapse single well ERT and crosswell EM data have been acquired to image CO2 injection into a shallow aquifer. Resistivity imaging results from inversion of the two data types separately and jointly will be compared and contrasted, and interpretation of the extent of the injected CO2 provided

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects

Search for eccentric black hole coalescences during the third observing run of LIGO and Virgo

Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass M>70 M⊙) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities 0<e≤0.3 at 0.33 Gpc−3 yr−1 at 90\% confidence level

Ultralight vector dark matter search using data from the KAGRA O3GK run

Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for U(1)B−L gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the U(1)B−L gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM

Crustal tomography of the Pyrenees and the surrounding regions by noise correlation

Durant cette thèse, nous appliquons la méthode de la corrélation de bruit ambiant dans la région particulièrement hétérogène des Pyrénées et des ses alentours (socle rocheux affleurant et bassin sédimentaires épais). Le jeu de données utilisé est une combinaison de deux réseaux temporaires large-bande français et espagnol (PYROPE et IBERARRAY) et de stations des réseaux permanents large-bande français et catalan. Le bruit sismique, enregistré pendant un an par les 158 stations est utilisé pour calculer les corrélations dans la gamme de période 5-55 s. Les vitesses de groupe de l'onde de Rayleigh et de l'onde de Love entre paires de stations sont inversées de manière linéarisée et nous obtenons, pour chaque période, des cartes de vitesses de groupe avec une résolution latérale d'environ 40 km. La comparaison entre les deux types d'ondes montre qu'il existe une anisotropie radiale à courtes périodes, alors que peu ou pas d'anisotropie radiale n'est visible aux périodes plus longues. Nous avons développé une nouvelle stratégie d'inversion des courbes de dispersion en modèle de vitesse d'onde S que nous avons appliqué sur les vitesses de groupe de l'onde de Rayleigh. Cette approche d'inversion est basée sur une exploration complète de l'espace des modèles et une inversion linéarisée. Le modèle obtenu, validé par la comparaison avec des résultats provenant d'autres méthodes, est le premier modèle 3-D crustale en vitesse d'onde S de la région et il permet d'apporter des contraintes sur la géodynamique des Pyrénées et de ses alentours. Deux points importants sont soulevés : (1) Des profils de vitesses atypiques sous l'Est du Massif Central, avec une croûte amincie et des vitesses anormalement faibles dans le manteau supérieur. (2) Deux anomalies de vitesses rapides sous la zone du Labourd-Mauléon et dans le prolongement du bassin de Parentis. Ces anomalies, situées à 25 km de profondeur, sont interprétées comme les traces de l'hyper-extension qui aurait précédée la phase de collision amenant à la formation des Pyrénées. La forte hétérogénéité de la zone permet également de faire l'analyse de l'influence de la non prise en compte des déviations de rais lors de l'inversion. Les premiers résultats montrent que le modèle obtenu en utilisant la théorie des rais droits ne permet pas d'expliquer les déviations calculées par la méthode de formation de voie. Ces déviations observées peuvent donc apporter une amélioration du modèle en considérant l'utilisation d'une inversion combinée.In this thesis, we applied the ambient noise correlation method in the very heterogeneous region of the Pyrenees and the surrounding areas (mountain belt and thick sedimentary basins). The dataset used is a combination of two temporary broadband arrays from France and Spain (PYROPE and IBERARRAY) and stations of the French and Catalan permanent broadband arrays. Seismic noise recorded over years by the 158 stations was used to calculate correlations in a period range of 5-55 s. Observed Rayleigh and Love wave group velocities between pairs of stations were used as input to a linearized inversion scheme, where we obtained for each period group velocities maps, with a lateral resolution of approximately 40 km. The comparison between the two type of waves demonstrates radial anisotropy at short periods, while little or no radial anisotropy is present at long periods. We developed a new strategy for the inversion of dispersion curves to shear wave velocity models and applied it to the Rayleigh waves group velocity. This approach is based on the combination of a full exploration of the model space and a linearized inversion. The obtained model, validated by the comparison of our results with the results of other methods, is the first complete 3-D crustal Vs model of the region. We in particular note : (1) Atypical S-waves profiles in the Eastern part of the Massif Central which indicate a thinned crust and low velocities in the uppermost mantle. (2) High-velocity anomalies at 25 km depth beneath the Labourd-Mauléon area and the on-land continuation of the Parentis basin. We suggest that they are the traces of the hyper-extension which might have preceded the collision phase which lead to the formation of the Pyrenees. The strong heterogeneity of our study region is also well-adapted to analyze the influence of the ray deviations on the reconstruction of the model. First results show that the reconstructed model, using great-circle paths, does not explain the ray deviations as observed by beamforming. Observed deviations therefore carry the potential of improving the model in a combined inversion scheme

Long-duration transient gravitational-wave search pipeline

International audienceAs the sensitivity and observing time of gravitational-wave detectors increase, a more diverse range of signals is expected to be observed from a variety of sources. Especially, long-lived gravitational-wave transients have received interest in the last decade. Because most of long-duration signals are poorly modeled, detection must rely on generic search algorithms, which make few or no assumption on the nature of the signal. However, the computational cost of those searches remains a limiting factor, which leads to sub-optimal sensitivity. Several detection algorithms have been developed to cope with this issue. In this paper, we present a new data analysis pipeline to search for un-modeled long-lived transient gravitational-wave signals with duration between 10 and 1000 s, based on an excess cross-power statistic in a network of detectors. The pipeline implements several new features that are intended to reduce computational cost and increase detection sensitivity for a wide range of signal morphologies. The method is generalized to a network of an arbitrary number of detectors and aims to provide a stable interface for further improvements. Comparisons with a previous implementation of a similar method on simulated and real gravitational-wave data show an overall increase in detection efficiency depending on the signal morphology, and a computing time reduced by at least a factor 10

Ambient noise tomography of the Pyrenees and the surrounding regions: Inversion for a 3-D Vs model in the presence of a very heterogeneous crust

The lithospheric architecture of the Pyrenees is still uncertain and highly debated. Here, we provide new constraints from a high-resolution 3-D S-wave velocity model of the Pyrenees and the adjacent foreland basins. This model is obtained from ambient noise tomography on records of temporary and permanent seismic arrays installed in southwestern France and northern Spain. We first computed group velocity maps for Rayleigh waves in the 5 to 55 s period range using noise correlation stacks at 1500-8500 station pairs. As the crust is very heterogeneous, poor results were obtained using a single starting model in a linearized inversion of group velocity dispersion curves for the shear wave structure. We therefore built a starting model for each grid node by full exploration of the model space. The resulting 3-D shear wave velocity model is compared to data from previous geophysical studies as a validation test. Despite the poor sensitivity of surface waves to seismic discontinuities, the geometry of the top of the basement and the Moho depth are retrieved well, except along the Cantabrian coast. Major reflectors of the ECORS deep seismic sounding profiles in the central and western Pyrenees coincide with sharp velocity gradients in our velocity model.We retrieve the difference between the thicker Iberian crust and the thinner European crust, the presence of low-velocity material of the Iberian crust underthrust beneath the European crust in the central Pyrenees, and the structural dissymmetry between the South Pyrenean Zone and the North Pyrenean Zone at the shallow crustal level. In the Labourd-Mauléon-Arzacq region (western Pyrenees), there is a high S-wave velocity anomaly at 20-30 km in depth, which might explain the positive Bouguer anomaly of the Labourd Massif. This high-velocity lower crust, which is also detected beneath the Parentis area, might be an imprint of the Albian-Aptian rifting phase. The southeastern part of the Massif Central has an unusual velocity structure, with a very shallowMoho (21-25 km) above an uppermost mantle with anomalously low shear wave velocity.We thank all participants in the fieldwork, and the municipalities and landlords that hosted a PYROPE temporary station. We also acknowledge SISMOB, the French seismic mobile pool (a com- ponent of the RESIF Research Facility), for providing us with the seismological instrumentation for the temporary deployments. We used data from the FR and RD (RESIF), G (Geoscope) and CA (‘Institut Cartogr ` afic i Geol ` ogic de Catalunya’) permanent net- works. RESIF ( http://portal.resif.fr/ ) is a national Research Infras- tructure, recognized as such by the French Ministry of Higher Education and Research. RESIF is managed by the RESIF Con- sortium, composed of 18 Research Institutions and Universities in France. RESIF additionally supported by a public grant overseen by the French National Research Agency (ANR) as part of the «Investissements d’Avenir» program (reference: ANR-11-EQPX- 0040)andtheFrenchMinistryofEcology,SustainableDevelopment and Energy. The PYROPE experiment was supported by the French Research Agency “ANR blanc” program (project PYROPE, ANR-09-BLAN-0229). This is a contribution of the Team Consolider-Ingenio 2010 TOPO-IBERIA (CSD2006–00041). CPeer reviewe

Searches for long-duration gravitational wave transients in LIGO and Virgo data

International audienceLong-lived gravitational wave (GW) transients have received interest in the last decade, as the sensitivity of LIGO and Virgo increases. Such signals, lasting between 10 and 1000s, can come from a variety of sources, including accretion disk instabilities around black holes, binary neutron stars post-merger, core-collapse supernovae, non-axisymmetric deformations in isolated neutron stars, and magnetar giant flares. Given the large parameter space and the lack of precisely modeled waveforms, searches must rely on robust detection algorithms, which make few or no assumptions on the nature of the signal. Here we present a new data analysis pipeline to search for long-lived transient GW signals, based on an excess cross-power statistic computed over a network of detectors. It uses a hierarchical strategy that allows to estimate the background quickly and implements several features aimed to increase detection sensitivity by 30% for a wide range of signal morphology compared to an older implementation. We also report upper limits on the GW energy emitted from a search conducted with the pipeline for GW emission around a sample of nearby magnetar giant flares, and discuss detection potential of such sources with second and third-generation detectors