A Reappraisal of the H-κ Stacking Technique: Implications for Global Crustal Structure

TY - JOUR

T1 - A Reappraisal of the H-κ Stacking Technique

T2 - Implications for Global Crustal Structure

AU - Ogden, C. S.

AU - Bastow, I. D.

AU - Gilligan, A.

AU - Rondenay, S.

N1 - We thank two anonymous reviewers and editor Michael Ritzwoller for insightful comments which have improved this manuscript. We also thank H. Meek for hard work during the early stages of this project and S. Pilidou, I. Dimitriadis, P. Iosif and their colleagues at the Geological Survey Department of Cyprus for their help establishing the TROODOS network (Bastow et al., 2017). V. Lane and D. Daly (both of SEIS-UK), A. Boyce, M. Liddell and R. Kounoudis were all excellent field assistants in Cyprus. SAC (Helffrich et al., 2013) and GMT (Wessel and Smith, 1991) software were used to process and image seismic data, which were sourced from IRIS DMC and ORFEUS. C.S. Ogden is funded by the Natural Environment Research Council (NERC) Doctoral Training Partnership: Science and Solutions for a Changing Planet, Grant Number NE/L002515/1. S. Rondenay’s contribution to this work was supported by Career Integration Grant 321871 - GLImER from the FP7 Marie Curie Actions of the European Commission, and by the Research Council of Norway FRINATEK programme through SwaMMIS project 231354.

PY - 2019/12

Y1 - 2019/12

N2 - H-κ stacking is used routinely to infer crustal thickness and bulk-crustal VP /VS ratio from teleseismic receiver functions. The method assumes that the largest amplitude P-to-S conversions beneath the seismograph station are generated at the Moho. This is reasonable where the crust is simple and the Moho marks a relatively abrupt transition from crust to mantle, but not if the crust-mantle transition is gradational and/or complex intra-crustal structure exists. We demonstrate via synthetic seismogram analysis that H-κ results can be strongly dependent on the choice of stacking parameters (the relative weights assigned to the Moho P-to-S conversion and its subsequent reverberations, the choice of linear or phase-weighted stacking, input crustal P-wave velocity) and associated data parameters (receiver function frequency content and the sample of receiver functions analyzed). To address this parameter sensitivity issue, we develop an H-κ approach in which cluster analysis selects a final solution from 1000 individual H-κ results, each calculated using randomly-selected receiver functions, and H-κ input parameters. Ten quality control criteria that variously assess the final numerical result, the receiver function dataset, and the extent to which the results are tightly clustered, are used to assess the reliability of H-κ stacking at a station. Analysis of synthetic datasets indicates H-κ works reliably when the Moho is sharp and intra-crustal structure is lacking but is less successful when the Moho is gradational. Limiting the frequency content of receiver functions can improve the H-κ solutions in such settings, provided intra-crustal structure is simple. In cratonic Canada, India and Australia, H-κ solutions generally cluster tightly, indicative of simple crust and a sharp Moho. In contrast, on the Ethiopian plateau, where Paleogene flood-basalts overlie marine sediments, H-κ results are unstable and erroneous. For stations that lie on thinner flood-basalt outcrops, and/or in regions where Blue Nile river incision has eroded through to the sediments below, limiting the receiver function frequency content to longer periods improves the H-κ solution and reveals a 6–10 km gradational Moho, readily interpreted as a lower-crustal intrusion layer at the base of a mafic (VP /VS=1.77–1.87) crust. Moving off the flood-basalt province, H-κ results are reliable and the crust is thinner and more felsic (VP /VS=1.70–1.77), indicating the lower crustal intrusion layer is confined to the region covered by flood-basaltic volcanism. Analysis of data from other tectonically-complex settings (e.g., Japan, Cyprus) shows H-κ stacking results should be treated cautiously. Only in regions of relatively simple crust can H-κ stacking analysis be considered truly reliable.

AB - H-κ stacking is used routinely to infer crustal thickness and bulk-crustal VP /VS ratio from teleseismic receiver functions. The method assumes that the largest amplitude P-to-S conversions beneath the seismograph station are generated at the Moho. This is reasonable where the crust is simple and the Moho marks a relatively abrupt transition from crust to mantle, but not if the crust-mantle transition is gradational and/or complex intra-crustal structure exists. We demonstrate via synthetic seismogram analysis that H-κ results can be strongly dependent on the choice of stacking parameters (the relative weights assigned to the Moho P-to-S conversion and its subsequent reverberations, the choice of linear or phase-weighted stacking, input crustal P-wave velocity) and associated data parameters (receiver function frequency content and the sample of receiver functions analyzed). To address this parameter sensitivity issue, we develop an H-κ approach in which cluster analysis selects a final solution from 1000 individual H-κ results, each calculated using randomly-selected receiver functions, and H-κ input parameters. Ten quality control criteria that variously assess the final numerical result, the receiver function dataset, and the extent to which the results are tightly clustered, are used to assess the reliability of H-κ stacking at a station. Analysis of synthetic datasets indicates H-κ works reliably when the Moho is sharp and intra-crustal structure is lacking but is less successful when the Moho is gradational. Limiting the frequency content of receiver functions can improve the H-κ solutions in such settings, provided intra-crustal structure is simple. In cratonic Canada, India and Australia, H-κ solutions generally cluster tightly, indicative of simple crust and a sharp Moho. In contrast, on the Ethiopian plateau, where Paleogene flood-basalts overlie marine sediments, H-κ results are unstable and erroneous. For stations that lie on thinner flood-basalt outcrops, and/or in regions where Blue Nile river incision has eroded through to the sediments below, limiting the receiver function frequency content to longer periods improves the H-κ solution and reveals a 6–10 km gradational Moho, readily interpreted as a lower-crustal intrusion layer at the base of a mafic (VP /VS=1.77–1.87) crust. Moving off the flood-basalt province, H-κ results are reliable and the crust is thinner and more felsic (VP /VS=1.70–1.77), indicating the lower crustal intrusion layer is confined to the region covered by flood-basaltic volcanism. Analysis of data from other tectonically-complex settings (e.g., Japan, Cyprus) shows H-κ stacking results should be treated cautiously. Only in regions of relatively simple crust can H-κ stacking analysis be considered truly reliable.

KW - Crustal imaging

KW - Crustal structure

KW - Cratons

KW - Large igneous provinces

KW - Statistical methods

KW - Body waves

KW - Ethiopia

UR - https://academic.oup.com/journals/pages/access_purchase/rights_and_permissions/self_archiving_policy_p

UR - http://www.mendeley.com/research/reappraisal-h%CE%BA-stacking-technique-implications-global-crustal-structure

UR - http://www.scopus.com/inward/record.url?scp=85072769027&partnerID=8YFLogxK

U2 - 10.1093/gji/ggz364

DO - 10.1093/gji/ggz364

M3 - Article

VL - 219

SP - 1491

EP - 1513

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 3

ER -