Guilhem Barruol | Université de la Réunion (original) (raw)

Papers by Guilhem Barruol

The Precambrian rocks of Madagascar were formed and/or modified during continental collision know... more The Precambrian rocks of Madagascar were formed and/or modified during continental collision known as the Pan-African orogeny. Aborted Permo-Triassic Karoo rifting and the subsequent separation from Africa and India resulted in the formation of sedimentary basins in the west and volcanic activity predominantly along the margins. Many geological studies have documented the imprint of these processes, but little was known about the deeper structure. We therefore deployed seismic stations along an SE-NW trending profile spanning nearly all geological domains of southern Madagascar. Here we focus on the crustal structure, which we determined based on joint analysis of receiver functions and surface waves derived from ambient noise measurements. For the sedimentary basin we document a thinning of the underlying crystalline basement by up to ∼60% to 13 km. The crustal velocity structure demonstrates that the thinning was accomplished by removal or exhumation of the lower crust. Both the Proterozoic and Archean crust have a 10 km thick upper crust and 10–12 km thick midcrust. However, in contrast to the typical structure of Proterozoic and Archean aged crust, the Archean lower crust is thicker and faster than the Proterozoic one, indicating possible magmatic intrusions; an underplated layer of 2–8 km thickness is present only below the Archean crust. The Proterozoic mafic lower crust might have been lost during continental collision by delamination or subduction or thinned as a result of extensional collapse. Finally, the Cretaceous volcanics along the east coast are characterized by thin crust (30 km) and very large V P ∕V S ratios.

Keywords: surface wave ambient noise tomography Madagascar intraplate volcanism The crust and upp... more Keywords: surface wave ambient noise tomography Madagascar intraplate volcanism The crust and upper mantle of the Madagascar continental fragment remained largely unexplored until a series of recent broadband seismic experiments. An island-wide deployment of broadband seismic instruments has allowed the first study of phase velocity variations, derived from surface waves, across the entire island. Late Cenozoic alkaline intraplate volcanism has occurred in three separate regions of Madagascar (north, central and southwest), with the north and central volcanism active until <1 Ma, but the sources of which remains uncertain. Combined analysis of three complementary surface wave methods (ambient noise, Rayleigh wave cross-correlations, and two-plane-wave) illuminate the upper mantle down to depths of 150 km. The phase-velocity measurements from the three methods for periods of 8–182 s are combined at each node and interpolated to generate the first 3-D shear-velocity model for sub-Madagascar velocity structure. Shallow (upper 10 km) low-shear-velocity regions correlate well with sedimentary basins along the west coast. Upper mantle low-shear-velocity zones that extend to at least 150 km deep underlie the north and central regions of recent alkali magmatism. These anomalies appear distinct at depths <100 km, suggesting that any connection between the zones lies at depths greater than the resolution of surface-wave tomography. An additional low-shear velocity anomaly is also identified at depths 50–150 km beneath the southwest region of intraplate volcanism. We interpret these three low-velocity regions as upwelling asthenosphere beneath the island, producing high-elevation topography and relatively low-volume magmatism.

We present two independent, automated methods for estimating the absolute horizontal misori-entat... more We present two independent, automated methods for estimating the absolute horizontal misori-entation of seismic sensors. We apply both methods to 44 free-fall ocean-bottom seismometers (OBSs) of the RHUM-RUM experiment (http://www.rhum-rum.net/). The techniques measure the 3-D directions of particle motion of (1) P-waves and (2) Rayleigh waves of earthquake recordings. For P-waves, we used a principal component analysis to determine the directions of particle motions (polarizations) in multiple frequency passbands. We correct for polarization deviations due to seismic anisotropy and dipping discontinuities using a simple fit equation, which yields significantly more accurate OBS orientations. For Rayleigh waves, we evaluated the degree of elliptical polarization in the vertical plane in the time and frequency domain. The results obtained for the RHUM-RUM OBS stations differed, on average, by 3.1 • and 3.7 • between the methods, using circular mean and median statistics, which is within the methods' estimate uncertainties. Using P-waves, we obtained orientation estimates for 31 ocean-bottom seismometers with an average uncertainty (95 per cent confidence interval) of 11 • per station. For 7 of these OBS, data coverage was sufficient to correct polarization measurements for underlying seismic anisotropy and dipping discontinuities, improving their average orientation uncertainty from 11 • to 6 • per station. Using Rayleigh waves, we obtained misorientation estimates for 40 OBS, with an average uncertainty of 16 • per station. The good agreement of results obtained using the two methods indicates that they should also be useful for detecting misorientations of terrestrial seismic stations.

S U M M A R Y Ocean wave activity excites seismic waves that propagate through the solid earth, k... more S U M M A R Y Ocean wave activity excites seismic waves that propagate through the solid earth, known as microseisms, which, once recorded on oceanic islands, can be used to analyse the swell. Here, we analyse the microseismic noise recorded in different period ranges by the permanent seismic station RER on La Réunion Island and by a temporary network of 10 broad-band seismic stations deployed on the island to analyse extreme swell events. We perform a comparative analysis of cyclonic and austral swell events by analysing not only the primary (PM, ∼10– 20 s period) and secondary (SM, ∼3–10 s) microseisms but also the long-period secondary microseisms (LPSMs, ∼ 7–10 s), which may result from the interaction between incident ocean waves and the reflected waves off the coast. We compare the microseismic observations with buoy data when available and with hindcasts from numerical ocean wave models. We show that each cyclone is characterized by its own individual signature in the SM, which depends not only on its distance and intensity but also on its dynamics and trajectory. Thus, the SM contains relevant information for cyclone detection and monitoring. Analysing the PM and the LPSM, and comparing it to direct buoy observations and/or wave numerical models allows characterizing the local impact of the swell with the island in terms of amplitude, period, and sometimes, direction of propagation, making possible to use a seismic station as an ocean wave gauge. The microseisms, which link the atmosphere, the ocean and the solid Earth, can thus provide valuable observations on extreme swell events, in addition to oceanic and meteorological data.

RHUM-RUM is a German-French seismological experiment based on the sea floor surrounding the islan... more RHUM-RUM is a German-French seismological experiment based on the sea floor surrounding the island of La Réunion, western Indian Ocean (Barruol and Sigloch, 2013). Its primary objective is to clarify the presence or absence of a mantle plume beneath the Reunion volcanic hotspot. RHUM-RUM's central component is a 13-month deployment (October 2012 to November 2013) of 57 broad-band ocean bottom seismometers (OBS) and hydrophones over an area of 2000 ⇥ 2000 km 2 surrounding the hotspot. The array contained 48 wideband OBS from the German DE-PAS pool and 9 broadband OBS from the French INSU pool. It is the largest deployment of DEPAS and INSU OBS so far, and the first joint experiment. This article reviews network performance and data quality: of the 57 stations, 46 and 53 yielded good seismome-ter and hydrophone recordings, respectively. The 19 751 total deployment days yielded 18 735 days of hydrophone recordings and 15 941 days of seismometer recordings, which are 94 and 80 % of the theoretically possible yields. The INSU seismic sensors stand away from their OBS frames, whereas the DEPAS sensors are integrated into their frames. At long periods (> 10 s), the DEPAS seismome-ters are affected by significantly stronger noise than the INSU seismometers. On the horizontal components, this can be explained by tilting of the frame and buoy assemblage, e.g. through the action of ocean-bottom currents, but in addition the DEPAS intruments are affected by significant self-noise at long periods, including on the vertical channels. By comparison, the INSU instruments are much quieter at periods > 30 s and hence better suited for long-period signals studies. The trade-off of the instrument design is that the integrated DEPAS setup is easier to deploy and recover, especially when large numbers of stations are involved. Additionally , the wideband sensor has only half the power consumption of the broadband INSU seismometers. For the first time, this article publishes response information of the DE-PAS instruments, which is necessary for any project where true ground displacement is of interest. The data will become publicly available at the end of 2017.

Ocean waves activity is a major source of microvibrations that travel through the solid Earth, kn... more Ocean waves activity is a major source of microvibrations that travel through the solid Earth, known as microseismic noise and recorded worldwide by broadband seismometers. Analysis of microseismic noise in continuous seismic records can be used to investigate noise sources in the oceans such as storms, and their variations in space and time, making possible the regional and global-scale monitoring of the wave climate. In order to complete the knowledge of the Atlantic and Pacific oceans microseismic noise sources, we analyse 1 yr of continuous data recorded by permanent seismic stations located in the Indian Ocean basin.We primarily focus on secondary microseisms (SM) that are dominated by Rayleigh waves between 6 and 11 s of period. Continuous polarization analyses in this frequency band at 15 individual seismic stations allow us to quantify the number of polarized signal corresponding to Rayleigh waves, and to retrieve their backazimuths (BAZ) in the time–frequency domain. We observe clear seasonal variations in the number of polarized signals and in their frequencies, but not in their BAZ that consistently point towards the Southern part of the basin throughout the year. This property is very peculiar to the Indian Ocean that is closed on its Northern side, and therefore not affected by large ocean storms during Northern Hemisphere winters. We show that the noise amplitude seasonal variations and the backazimuth directions are consistent with the source areas computed from ocean wave models.

Geophysical Research Letters, 2014

Eos, Transactions American Geophysical Union, 2013

The Piton de la Fournaise basaltic volcano, on La Reunion Island in the western Indian Ocean, is ... more The Piton de la Fournaise basaltic volcano, on La Reunion Island in the western Indian Ocean, is one of the most active volcanoes in the world. This volcano is classically considered as the surface expression of an upwelling mantle plume and its activity is continuously monitored, providing detailed information on its superficial dynamics and on the edifice structure. Deeper crustal and upper mantle structure under La Reunion Island is surprisingly poorly constrained,
motivating this study.We used receiver function techniques to determine a shear wave velocity profile through the crust and uppermost mantle beneath La Reunion, but also at other seismic stations located on the hotspot track, to investigate the plume and lithosphere interaction and
its evolution through time. Receiver functions (RFs) were computed at permanent broadband seismic stations from the GEOSCOPE network (on La Reunion and Rodrigues), at
IRIS stations MRIV and DGAR installed on Mauritius and Diego Garcia islands, and at the GEOFON stations KAAM and HMDM on the Maldives. We performed non-linear inversions
of RFs through modelling of P-to-S conversions at various crustal and upper mantle interfaces. Joint inversion of RF and surface wave dispersion data suggests a much deeper Mohorovic discontinuity (Moho) beneath Mauritius (∼21 km) compared to La Reunion (∼12 km). A magmatic underplated body may be present under La Reunion as a thin layer (≤3 km thick), as suggested by a previous seismic refraction study, and as a much thicker layer beneath other
stations located on the hotspot track, suggesting that underplating is an important process resulting from the plume–lithosphere interaction.We find evidence for a strikingly lowvelocity layer starting at about 33 km depth beneath La Reunion that we interpret as a zone of partial melt beneath the active volcano. We finally observe low velocities below 70 km beneath La Reunion and below 50 km beneath Mauritius that could represent the base of the oceanic lithosphere.

Geophysical Research Letters, 2013

S U M M A R Y A temporary network of 10 broad-band seismic stations has been installed in French ... more S U M M A R Y A temporary network of 10 broad-band seismic stations has been installed in French Polynesia for the Polynesian Lithosphere and Upper Mantle Experiment (PLUME). All the seismic stations were installed either on volcanic islands or on atolls of the various archipelagos of French Polynesia in a manner which complements the geographic coverage provided by the regional permanent stations. The primary aim of PLUME is to image the upper mantle structures related to plate motion and hotspot activity. However, because of its proximity to all sites, the ocean is responsible for a high level of noise in the seismic data and we show that these data can also be used to analyse ocean wave activity. The power spectral density (PSD) analyses of the seismic data recorded in French Polynesia show clear peaks in the 0.05– 0.10 Hz band (periods between 10 and 20 s), which corresponds to swell frequencies. Clear peaks in this frequency band are also observed in infrasonic data recorded on Tahiti. Ground motion analysis shows that the swell-related seismic noise (SRSN) is linearly polarized in the horizontal plane and its amplitude decreases rapidly with the distance from the shore. The microseismic and the infrasonic 'noise' amplitudes show very similar variations from station to station and both are strongly correlated with the swell amplitudes predicted by the National Oceanic and Atmospheric Administration (NOAA), wind-forced, 'WaveWatch' models. The swell direction can be estimated from SRSN polarization analysis but this has to be done with care since, for some cases, the ground motions are strongly controlled by the islands' anisometric shapes and by swell refraction processes. We find cases, however, such as Tahiti or roughly circular Tuamotu atolls, where the azimuth of the swell is in good agreement with the seismic estimates. We, therefore, demonstrate that the SRSN and the infrasonic signal observed in French Polynesia can be used in such cases as a proxy for swell amplitude and azimuth. From the continuous analysis of the data recorded in 2003 at the permanent seismic station PPTL in Tahiti, transfer functions have been obtained. This could provide a way to quantify the swell activity during the last two decades and, therefore, assist in the investigation of climate changes.

Seismic records from La Réunion Island very broadband Geoscope station are investigated to constr... more Seismic records from La Réunion Island very broadband Geoscope station are investigated to constrain the link between the 2007 eruptive sequence and the related caldera collapse of the Piton de la Fournaise volcano. Tilt estimated from seismic records reveals that the three 2007 eruptions belong to a single inflation-deflation cycle. Tilt trend indicates that the small-volume summit eruption of 18 February occurred during a phase of continuous inflation that started in January 2007. Inflation decelerated 24 days before a second short-lived, small-volume eruption on 30 March, almost simultaneous with a sudden, large-scale deflation of the volcano. Deflation rate, which had stabilized at relatively low level, increased anew on 1 April while no magma was erupted, followed on 2 April by a major distal eruption and on 5 April by a summit caldera collapse. Long-term tilt variation suggests that the 2007 eruptive succession was triggered by a deep magma input.

Geophysical Monograph Series, 2000

Eos, Transactions American Geophysical Union, 2002

Geophysical Journal International

The Terre Adélie and George V Land (East Antarctica) represent key areas for understand-ing tecto... more The Terre Adélie and George V Land (East Antarctica) represent key areas for understand-ing tectonic relationships between terranes forming the Neoarchean-Palaeoproterozoic Terre Adélie Craton (TAC) and the neighbouring lithospheric blocks, together with the nature of its boundary. This region that represents the eastern border of the TAC is limited on its eastern side by the Mertz shear zone (MSZ) separating more recent Palaeozoic units from the craton. The MSZ, that recorded dextral strike-slip movement at 1.7 and 1.5 Ga, is likely correlated with the Kalinjala or Coorong shear zone in South Australia, east of the Gawler Craton and may therefore represent a frozen lithospheric-scale structure. In order to investigate the litho-spheric structure of the TAC and the MSZ, we deployed from 2009 October to 2011 October four temporary seismic stations, which sampled the various lithospheric units of the TAC and of the neighbouring Palaeozoic block, together with the MSZ. We used receiver...

Eos, Transactions American Geophysical Union, 2005