Jacques Marteau

Improvement of density models of geological structures by fusion of gravity data and cosmic muon radiographies

Both muon tomography and gravimetry are geophysical methods that provide information on the density structure of the Earth’s subsurface. Muon tomography measures the natural flux of cosmic muons and its attenuation produced by the screening effect of the rock mass to image. Gravimetry generally consists in measurements of the vertical component of the local gravity field. Both methods are linearly linked to density, but their spatial sensitivity is very different. Muon tomography essentially works like medical X-ray scan and integrates density information along elongated narrow conical volumes while gravimetry measurements are linked to density by a 3-dimensional integral encompassing the whole studied domain. We develop the mathematical expressions of these integration formulas – called acquisition kernels – to express resolving kernels that act as spatial filters relating the true unknown density structure to the density distribution actually recoverable from the available data. The resolving kernels provide a tool to quantitatively describe the resolution of the density models and to evaluate the resolution improvement expected by adding new data in the inversion. The resolving kernels derived in the joined muon/gravimetry case indicate that gravity data are almost useless to constrain the density structure in regions sampled by more than two muon tomography acquisitions. Interestingly the resolution in deeper regions not sampled by muon tomography is significantly improved by joining the two techniques. Examples taken from field experiments performed on La Soufrière of Guadeloupe volcano are discussed.

Muon dynamic radiography of density changes induced by hydrothermal activity at the La Soufrière of Guadeloupe volcano

Imaging geological structures through cosmic muon radiography is a newly developed technique which shows a great potential in volcanology. Here we demonstrate that muon radiography permits to detect and characterize mass movements in shallow hydrothermal systems of low-energy active volcanoes like the La Soufrière lava dome. We present an experiment conducted on this volcano during the Summer 2014 and bring evidence that very important density changes occurred in three domains of the lava dome. Depending on their position and on the medium porosity the volumes of these domains vary from 1u2009×u2009106 m3 to 7u2009×u2009106 m3. However, the total mass budget remains approximately constant : two domains show a mass loss (Δm∈ [−0.8;−0.4]u2009×u2009109u2009kg) and the third one a mass gain (Δm∈ [1.5; 2.5]u2009×u2009109u2009kg). We attribute the negative mass changes to the formation of steam in shallow hydrothermal reservoir previously partly filled with liquid water. This coincides with the emergence of new fumaroles on top of the volcano. The positive mass change is synchronized with the negative mass changes indicating that liquid water probably flowed from the two reservoirs invaded by steam toward the third reservoir.

Monitoring temporal opacity fluctuations of large structures with muon tomography : a calibration experiment using a water tower tank

Usage of secondary cosmic muons to image the geological structures density distribution significantly developed during the past ten years. Recent applications demonstrate the method interest to monitor magma ascent and volcanic gas movements inside volcanoes. Muon radiography could be used to monitor density variations in aquifers and the critical zone in the near surface. However, the time resolution achievable by muon radiography monitoring remains poorly studied. It is biased by fluctuation sources exterior to the target, and statistically affected by the limited number of particles detected during the experiment. The present study documents these two issues within a simple and well constrained experimental context: a water tower. We use the data to discuss the influence of atmospheric variability that perturbs the signal, and propose correction formulas to extract the muon flux variations related to the water level changes. Statistical developments establish the feasibility domain of muon radiography monitoring as a function of target thickness (i.e. opacity). Objects with a thickness comprised between ≈50u2009±u200930u2009m water equivalent correspond to the best time resolution. Thinner objects have a degraded time resolution that strongly depends on the zenith angle, whereas thicker objects (like volcanoes) time resolution does not.

Muon tomography applied to activevolcanoes

Muon tomography is a generic imaging method using the differential absorption of cosmic muons by matter. The measured contrast in the muons flux reflects the matter density contrast as it does in conventional medical imaging. The applications to volcanology present may advantadges induced by the features of the target itself: limited access to dangerous zones, impossible use of standard boreholes information, harsh environmental conditions etc. The Diaphane project is one of the largest and leading collaboration in the field and the present article summarizes recent results collected on the Lesser Antilles, with a special emphasis on the Soufri`ere of Guadeloupe.

Feasibility study of archaeological structures scanning by muon tomography

One of the main concerns in archaeology is to find of a method to study precisely archaeological structures in the least invasive way possible to avoid damage. The requirement of preserving the structures integrity prevents, in the case of pyramids or tumuli, the study of any internal structure (halls or tombs) which are not reachable by existing corridors. n nOne non-invasive method is the muon tomography. By placing a detector which allows to register the muon direction after the structure, it is possible to have an idea of its composition based on the attenuation of the muon flux, which depends on the material length and density that muons have crossed. This technique, alone or together with other exploration techniques as seismic tomography or electrical resistivity tomography, can provide useful information about the internal structure of the archaeological form that can not be obtained by conventional archaeological methods. n nIn this work, the time measurement necessary to obtain a significant result about the composition of an archaeological structure is estimated. To do that, a Monte Carlo simulation framework based on the MUSIC software, properly tuned for this study, has been developed. The particular case of the Kastas Amfipoli Macedonian tumulus has been considered to perform the simulations.