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 -