skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: 3-D muographic inversion in the exploration of cavities and low-density fractured zones

Muography is an imaging tool based on the attenuation of cosmic muons for observing density anomalies associated with large objects, such as underground caves or fractured zones. Tomography based on muography measurements, that is, 3-D reconstruction of density distribution from 2-D muon flux maps, brings along special challenges. The detector field of view covering must be as balanced as possible, considering the muon flux drop at high zenith angles and the detector placement possibilities. The inversion from directional muon fluxes to a 3-D density map is usually underdetermined (more voxels than measurements). Therefore, the solution of the inversion can be unstable due to partial coverage. The instability can be solved by geologically relevant Bayesian constraints. However, the Bayesian principle results in parameter bias and artefacts. In this work, linearized (density-length based) inversion is applied by formulating the constraints associated with inversion to ensure the stability of parameter fitting. After testing the procedure on synthetic examples, an actual high-quality muography measurement data set from seven positions is used as input for the inversion. The resulting tomographic imaging provides details on the complicated internal structures of karstic fracture zone. The existence of low density zones in the imaged space was verified by samples from core drills, which consist of altered dolomite powder within the intact high density dolomite.

more » « less
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Medium: X Size: p. 700-710
["p. 700-710"]
Sponsoring Org:
National Science Foundation
More Like this

    The ability to accurately and reliably obtain images of shallow subsurface anomalies within the Earth is important for hazard monitoring and a fundamental understanding of many geologic structures, such as volcanic edifices. In recent years, machine learning (ML) has gained increasing attention as a novel approach for addressing complex problems in the geosciences. Here we present an ML-based inversion method to integrate cosmic-ray muon and gravity data sets for shallow subsurface density imaging at a volcano. Starting with an ensemble of random density anomalies, we use physics-based forward calculations to find the corresponding set of expected gravity and muon attenuation observations. Given a large enough ensemble of synthetic density patterns and observations, the ML algorithm is trained to recognize the expected spatial relations within the synthetic input–output pairs, learning the inherent physical relationships between them. Once trained, the ML algorithm can then interpolate the best-fitting anomalous pattern given data that were not used in training, such as those obtained from field measurements. We test the validity of our ML algorithm using field data from the Showa-Shinzan lava dome (Mt Usu, Japan) and show that our model produces results consistent with those obtained using a more traditional Bayesian joint inversion. Our results are similar to the previously published inversion, and suggest that the Showa-Shinzan lava dome consists of a relatively high-density (2200–2400 km m–3) cylindrical anomaly, about 300 m in diameter. Adding noise to synthetic training and testing data sets shows that, as expected, the ML algorithm is most robust in areas of high sensitivity, as determined by the forward kernels. Overall, we discover that ML offers a viable alternate method to a Bayesian joint inversion when used with gravity and muon data sets for subsurface density imaging.

    more » « less
  2. Abstract

    The Formosa array, with 137 broadband seismometers and ∼5 km station spacing, was deployed recently in Northern Taiwan. Here by using eight months of continuous ambient noise records, we construct the first high‐resolution three‐dimensional (3‐D) shear wave velocity model of the crust in the area. We first calculate multi‐component cross‐correlations to extract robust Rayleigh wave signals. We then determine phase velocity maps between 3 and 10 s periods using Eikonal tomography and measure Rayleigh wave ellipticity at each station location between 2 and 13 s periods. For each location, we jointly invert the two types of Rayleigh wave measurements with a Bayesian‐based inversion method for a one‐dimensional shear wave velocity model. All piecewise continuous one‐dimensional models are then used to construct the final 3‐D model. Our 3‐D model reveals upper crustal structures that correlate well with surface geological features. Near the surface, the model delineates the low‐velocity Taipei and Ilan Basins from the adjacent fast‐velocity mountainous areas, with basin geometries consistent with the results of previous geophysical exploration and geological studies. At a greater depth, low velocity anomalies are observed associated with the Linkou Tableland, Tatun Volcano Group, and a possible dyke intrusion beneath the Southern Ilan Basin. The model also provides new geometrical constraints on the major active fault systems in the area, which are important to understand the basin formation, orogeny dynamics, and regional seismic hazard. The new 3‐D shear wave velocity model allows a comprehensive investigation of shallow geologic structures in the Northern Taiwan.

    more » « less
  3. null (Ed.)
    SUMMARY We have developed a linear 3-D gravity inversion method capable of modelling complex geological regions such as subduction margins. Our procedure inverts satellite gravity to determine the best-fitting differential densities of spatially discretized subsurface prisms in a least-squares sense. We use a Bayesian approach to incorporate both data error and prior constraints based on seismic reflection and refraction data. Based on these data, Gaussian priors are applied to the appropriate model parameters as absolute equality constraints. To stabilize the inversion and provide relative equality constraints on the parameters, we utilize a combination of first and second order Tikhonov regularization, which enforces smoothness in the horizontal direction between seismically constrained regions, while allowing for sharper contacts in the vertical. We apply this method to the nascent Puysegur Trench, south of New Zealand, where oceanic lithosphere of the Australian Plate has underthrust Puysegur Ridge and Solander Basin on the Pacific Plate since the Miocene. These models provide insight into the density contrasts, Moho depth, and crustal thickness in the region. The final model has a mean standard deviation on the model parameters of about 17 kg m–3, and a mean absolute error on the predicted gravity of about 3.9 mGal, demonstrating the success of this method for even complex density distributions like those present at subduction zones. The posterior density distribution versus seismic velocity is diagnostic of compositional and structural changes and shows a thin sliver of oceanic crust emplaced between the nascent thrust and the strike slip Puysegur Fault. However, the northern end of the Puysegur Ridge, at the Snares Zone, is predominantly buoyant continental crust, despite its subsidence with respect to the rest of the ridge. These features highlight the mechanical changes unfolding during subduction initiation. 
    more » « less

    The fiducial cosmological analyses of imaging surveys like DES typically probe the Universe at redshifts z < 1. We present the selection and characterization of high-redshift galaxy samples using DES Year 3 data, and the analysis of their galaxy clustering measurements. In particular, we use galaxies that are fainter than those used in the previous DES Year 3 analyses and a Bayesian redshift scheme to define three tomographic bins with mean redshifts around z ∼ 0.9, 1.2, and 1.5, which extend the redshift coverage of the fiducial DES Year 3 analysis. These samples contain a total of about 9 million galaxies, and their galaxy density is more than 2 times higher than those in the DES Year 3 fiducial case. We characterize the redshift uncertainties of the samples, including the usage of various spectroscopic and high-quality redshift samples, and we develop a machine-learning method to correct for correlations between galaxy density and survey observing conditions. The analysis of galaxy clustering measurements, with a total signal to noise S/N ∼ 70 after scale cuts, yields robust cosmological constraints on a combination of the fraction of matter in the Universe Ωm and the Hubble parameter h, $\Omega _m h = 0.195^{+0.023}_{-0.018}$, and 2–3  per cent measurements of the amplitude of the galaxy clustering signals, probing galaxy bias and the amplitude of matter fluctuations, bσ8. A companion paper (in preparation) will present the cross-correlations of these high-z samples with cosmic microwave background lensing from Planck and South Pole Telescope, and the cosmological analysis of those measurements in combination with the galaxy clustering presented in this work.

    more » « less
  5. Abstract

    A self‐consistent regional‐scale seismic velocity model with resolution from seismogenic depth to the surface is crucial for seismic hazard assessment. Though Southern California is the most seismically imaged region in the world, techniques with high near‐surface sensitivity have been applied only in disparate local areas and have not been incorporated into a unified model with deeper resolution. In the present work, we obtain isotropic values for Rayleigh wave phase velocity and ellipticity in Southern California by cross‐correlating daily time series from the year 2015 across 315 regional stations in period ranges 6 to 18 s. Leveraging the complementary sensitivity of the two Rayleigh wave data sets, we combine H/V and phase velocity measurements to determine a new 3‐D shear velocity model in a Bayesian joint inversion framework. The new model has greatly improved shallow resolution compared to the Southern California Earthquake Center CVMS4.26 reference model. Well‐known large‐scale features common to previous studies are resolved, including velocity contrasts across the San Andreas, San Jacinto, Garlock, and Elsinore faults, midcrustal high‐velocity structure beneath the Mojave Desert, and shallow Moho beneath the Salton Trough. Other prominent features that have previously only been imaged in focused local studies include the correct sedimentary thickness of the southern Central Valley, fold structure of the Ventura and Oak Ridge Anticlines, and velocity contrast across the Newport‐Inglewood fault. The new shallow structure will greatly impact simulation‐based studies of seismic hazard, especially in the near‐surface low‐velocity zones beneath densely populated areas like the Los Angeles, San Bernardino, and Ventura Basins.

    more » « less