Surya Pachhai portrait
  • Research Assistant Professor, Geology & Geophysics

Research Statement

The analysis of seismological observations (e.g., seismic waves, Earth normal modes – the vibrational pattern due to large earthquakes) is a key component to inferring the structure and dynamic processes of the Earth. My research interests include the quantitative study of Earth structure and dynamic processes by inversion of seismological observations and the development of computational methods to elucidate Earth structure. Currently, my specific focus is on developing and applying probabilistic (Bayesian) inversion to achieve better understanding of the structure, composition, and dynamics of the solid Earth interior. While increased data quality, new inference techniques, and increased computational power have improved our understanding of the Earth’s interior, the interpretation of inversion results remains challenging because of inherent nonuniqueness (many models explain the data equally well) and nonlinearities. I assert that rigorous model parameter estimates and uncertainties are required for meaningful physical interpretation. In Bayesian inversion, parameters are considered as random variables and their values are represented by probability densities that express degrees of belief. Parameters are assigned prior information which is independent of the data and updated by data information to produce the inversion result (posterior probability density) from which rigorous parameter values and uncertainties can be inferred.

Core-reflected waves (i.e., PcP and ScP) have been extensively utilized to constrain small-scale structures that exhibit a strong decrease in shear (S)- and compressional (P)- wave velocity above the CMB, known as Ultra Low Velocity Zones (ULVZs). However, the details of extent and fine structure are poorly understood because of the nonuniqueness and nonlinearity of the inverse problem. Therefore, I developed a novel Bayesian approach to constrain ULVZ parameter values and uncertainties. I applied this method to ScP waves recorded by the Hi-Net array (Japan) that sample the CMB beneath the east of Philippines. This study found a complex ULVZ with dense material sitting at the bottom of two layers. An extension of this study generalizes the inversion to consider an unknown model parameterization (i.e. the number of layers and data noise are unknown — a trans-dimensional model). The new study produces more general, but consistent, results than our previous work and, in particular, more meaningful uncertainties which improves comparisons between data sets. Recently, I also applied the trans-dimensional (trans-D) approach to additional waveforms recorded by the Australian transportable array WOMBAT to study an extensive CMB area beneath the east of Australia.

Another project improves the measurement of splitting functions (a visualization of the radial average of the normal-mode sensed structure) for inner core sensitive modes, which are relatively well isolated in frequency, by inversion of normal mode data. Splitting functions depend linearly on isotropic and anisotropic velocities in the Earth. Inner core structure has been extensively studied using normal modes and body waves but is still controversial because of strong non-linearity. Therefore, my project applies nonlinear optimization (neighborhood algorithm — NA) to compute the splitting functions of the 13S2 (spheroidal mode with radial order 13 and angular order 2) mode. For the observed data, the splitting functions of this mode are dominated by zonal structure with positive frequency anomalies near the poles and negative frequency anomalies along the equator. One of the main benefits of the NA approach is that the splitting functions of this mode can be constrained with significantly fewer events in comparison to those using traditional approaches. Limitations of the NA approach are under sampling for higher orders (billions of samples are required) and lack of rigorous uncertainty estimation. To address the limitations, I have adapted the Bayesian ULVZ inversion to this problem and successfully applied for 19 inner-core sensitive modes. I am currently working on the application of this approach to the mantle sensitive modes.

I am also working on the crustal and mantle structure using both radial and transverse components of receiver functions and surface wave dispersions, mainly focused on Himalaya, Afar, and Iran.

Not only the deep Earth structure, I am also involved in the estimation of subsurface structure using spatial autocorrelation (SPAC) of ambient noise. SPAC represents the azimuthal average of normalized cross-correlation of a station at the center of the circle to the stations at the circumference of a circle. The 2-D shape of station geometry is preferred to reduce the directionality of noise. The SPAC method helps to constrain high-resolution subsurface structures (top 10s of meters) depending on the size of the array. In addition to the estimation of the S-wave velocity structure, I am also working on the modeling of electrical resistivity and gravity data to estimate the resistivity and density, respectively. I mostly develop and apply Bayesian inversion of these geophysical measurements to not only constrain the structure but also its uncertainties. Estimation of rigorous uncertainty in the Bayesian framework allows us to make meaningful interpretations of the results.

Publications

  • Sisay Alemayehu, Abdelkrim Aoudia, Atalay Ayele, Surya Pachhai, Hari Ram Thapa, CJ Ebinger, Radia Kherchouche, Mariangela Guidarelli & Seongryong Kim (2023). Structure of the crust-uppermost mantle beneath the Ethiopian volcanic province using ambient seismic noise and teleseismic P wave coda autocorrelation. Tectonophysics. Vol. 869, 230092. Published, 12/20/2023.
  • Surya Pachhai, Michael S. Thorne & Sebastian Rost (2023). Improved characterization of ultralow-velocity zones through advances in Bayesian inversion of ScP waveforms. Journal of Geophysical Research: Solid Earth. Published, 06/2023.
    https://doi.org/10.1029/2023JB026415
  • Hari Ram Thapa, Surya Pachhai, Abdelkrim Aoudia, Daniel Manu-Marfo, Keith Priestley & Supriyo Mitra (2023). The Main Himalayan Thrust beneath Nepal and Southern Tibet illuminated by seismic ambient noise and teleseismic P wave coda autocorrelation. Journal of Geophysical Research: Solid Earth. Published, 06/2023.
    https://doi.org/10.1029/2022JB026195
  • Bladimir Moreno Toiran, Abdelkrim Aoudia, Daniel Manu-Marfo, Radia Kherchouche & Surya Pachhai (2023). Crust-Uppermost Mantle Structure beneath the Caribbean Region from Seismic Ambient Noise Tomography. Bulletin of Seismological Society of America. Published, 04/2023.
    https://doi.org/10.1785/0120220062
  • Surya Pachhai, Michael S. Thorne & Tarje Nissen-Meyer (2022). Quantification of small-scale heterogeneity at the core-mantle boundary using sample entropy of SKS and SPdKS synthetic waveforms. Minerals. Published, 06/2022.
    https://www.mdpi.com/2075-163X/12/7/813
  • Najmieh Mohammadi, Habib Rahimi, Ali Gholami, Surya Pachhai & Abdelkrim Aoudia (2022). Shear-wave velocity structure of upper mantle along the Zagros collision zone. Tectonophysics. Published, 06/2022.
    https://www.sciencedirect.com/science/article/abs/...
  • Surya Pachhai, Mingming Li, Michael S. Thorne, Jan Dettmer & Hrvoje Tkalcic (2022). Internal structure of ultralow-velocity zones consistent with origin from a basal magma ocean. Nature Geoscience. Published, 01/2022.
    https://www.nature.com/articles/s41561-021-00871-5
  • Michael Thorne, Kuangdai Leng, Surya Pachhai, Sebastian Rost, June Wicks & Tarje Nissen-Meyer (2021). The most parsimonious ultralow‐velocity zone distribution from highly anomalous SPdKS waveforms. Geochemistry, Geophysics, Geosystems. Published, 01/14/2021.
    https://doi.org/10.1029/2020GC009467
  • Amir Sadeghi-Bagherabadi, Farhad Sobouti, Surya Pachhai & Abdelkrim Aoudia (2020). Estimation of Geometrical Spreading, Quality Factor and Kappa in the Zagros Region. Soil Dynamics and Earthquake Engineering. Published, 03/2020.
    https://doi.org/10.1016/j.soildyn.2020.106110
  • Surya Pachhai, Guy Masters & Gabi Laske (2020). Probabilistic estimation of structure coefficients and their uncertainties, for inner-core sensitive modes, using matrix autoregression. Geophysical Journal International. Published, 02/2020.
    https://doi.org/10.1093/gji/ggaa077
  • Michael Thorne, Surya Pachhai, Kuangdai Leng, June Wicks & Tarje Nissen-Meyer (2020). New candidate ultralow-velocity zone locations from highly anomalous . Minerals. Published, 02/2020.
    https://doi.org/10.3390/min10030211
  • Daniel Manu-Marfo, Abdelkrim Aoudia, Surya Pachhai & Radia Kherchouche (2019). 3D shear wave velocity model of the crust and uppermost mantle beneath the Tyrrhenian basin and margins. Scientific Reports. Published, 03/2019.
    https://doi.org/10.1038/s41598-019-40510-z
  • Surya Pachhai, Hrvoje Tkalcic & Guy Masters (2016). Estimation of splitting functions from normal mode spectra using the neighbourhood algorithm. Geophysical Journal International. Published, 01/2016.
    https://doi.org/10.1093/gji/ggv414
  • Surya Pachhai (2015). Bayesian inference for deep Earth structure using body waves and free oscillations of the Earth. The Australian National University. Published, 11/2015.
    http://hdl.handle.net/1885/149922
  • Surya Pachhai, Jan Dettmer & Hrvoje Tkalcic (2015). Ultra-low velocity zones beneath the Philippine and Tasman Seas revealed by trans-dimensional Bayesian waveform inversion. Geophysical Journal International. Published, 11/2015.
    https://doi.org/10.1093/gji/ggv368
  • Surya Pachhai, Hrvoje Tkalcic & Jan Dettmer (2014). Bayesian inference for ultralow velocity zones in the Earth’s lowermost mantle: complex ULVZ beneath the east of Philippines. Journal of Geophysical Research: Solid Earth. Published, 10/2014.
    https://doi.org/10.1002/2014JB011067