Ocean bottom seismograph orientation and crustal structure of the Woodlark Rift
-
摘要: 伍德拉克裂谷位于巴布亚新几内亚东南部,是发育在澳大利亚板块和西南太平洋板块碰撞带中的年轻大陆裂谷,为研究汇聚构造背景下裂谷起始演化的地壳结构提供了理想场所。伍德拉克裂谷海域地区海水层的存在使得获取高质量地震数据成为难题,而数据的获取主要依赖海底地震仪(ocean bottom seismograph,缩写为OBS)。OBS的布放一般是自由下落式,其地震计的北向水平分量方位与地理北向通常不一致,这使得利用三分量波形数据的反演结果产生较大误差甚至失效,例如接收函数方法。为确定伍德拉克裂谷地区OBS水平分量的方位偏转角度,本文同时引入纵波和瑞雷面波偏振分析方法进行方位校正,并利用校正后的三分量波形数据开展接收函数研究,进而约束该裂谷海域地区的地壳结构。结果分析表明,OBS方位校正后其获得的可用接收函数波形数量显著增多,并且利用纵波偏振分析校正后的数据处理获得了更加合理的地壳结构。基于在该裂谷地区获得的地壳构造结果,基里比斯盆地和裂谷扩张轴所在的古迪纳夫盆地呈现对比鲜明的地壳结构特征:古迪纳夫盆地的地壳厚度朝着裂谷扩张轴处减薄,其平均值为(33.3±2.42) km;基里比斯盆地的地壳厚度更薄,平均值为(24.1±5.44) km。此外,研究区域内所有OBS处均观测到了较高的地壳纵横波速比值,这可能是巴布亚超镁铁质岩体富集和古俯冲残片脱水熔融共同作用的结果。Abstract: The Woodlark Rift in southeastern Papua New Guinea is a young continental rift and develops within the collision zone between the Australian and SW Pacific Plates, which offers an ideal location to explore the crustal structure beneath the incipient rift under a convergent setting. However, the sea water layer makes it difficult to collect high-quality seismic data, and the common step is to deploy the ocean bottom seismographs (OBSs) in a free-fall way. Therefore, mis-orientation of horizontal components of the OBS usually leads to failure of applying the inversion techniques such as the receiver function to the three-component waveforms. In this study, we employed both P-wave and Rayleigh-wave polarization analyses to determine all available OBS orientations, and then used the recorded teleseismic waveforms to conduct a receiver function study on the crustal structure beneath the Woodlark rift. The number of the receiver function traces has greatly increased after the orientation corrections and the crustal structures can be better constrained based on the results from P-wave polarization analysis. Contrasting crustal structures were revealed beneath the Kiribishi Basin and the Goodenough Basin where the rift axis is located. The crust beneath the Goodenough Basin is deciphered to thin towards the rift axis with an average of (33.3±2.42) km, while a much thinner crust is observed beneath the Kiribisi Basin with a mean value of (24.1±5.44) km. High vP/vS ratios were determined at all stations, which may be attributed to the Papuan ultramafic body and dehydration melting of subducted slab segments.
-
Key words:
- Woodlark rift /
- crustal structure /
- OBS orientation /
- receiver function
-
图 1 伍德拉克裂谷地区地形图
紫色和黑色的圆圈分别表示含有效和无效数据的海底地震仪;红色三角为新生代火山;黑色虚线代表裂谷扩张轴,红色实线为欧文—斯坦利断裂带,绿色实线为地震测线(Fitz,Mann,2013)。右上插图代表了纵波(蓝色圆圈)和瑞雷波偏振分析(红色圆圈)所用的地震事件分布图,其中绿色三角形表示研究区域的中心位置。底部右侧插图中的黄色方块显示了研究区域的具体位置,其中红线为板块边界(Bird,2003)。AU:澳大利亚板块,SB:南俾斯麦板块,SS:所罗门海板块,WL:伍德拉克板块
Figure 1. Topographic map of the Woodlark rift showing the locations of the ocean bottom seismographs (OBSs)
The purple and black stars represent OBSs with and without valid data,respectively. Red triangles denote Cenozoic volcanos. The black dashed lines indicate the rift axis. The red line represents the Owen-Stanley fault. The green lines are seismic profiles from previous study (Fitz,Mann,2013). The top-right inset presents the events used for the polarization analysis of Rayleigh wave (red circles)and P-wave (blue circles). The green star marks the center of the study area. The bottom right inset displays the location of the study area highlighted by the yellow rectangle. The red lines denote the plate boundaries (Bird,2003). AU:Australia;SB:South Bismarck;SS:Solomon Sea;WL:Woodlark
图 2 (a) OBS水平分量的方位偏转示意图,χ和θ分别为OBS方位偏转角度和相对于北向水平分量的地震后方位角;(b) 地震台站B的瑞雷波偏振分析结果,图中蓝色圆圈为对应单事件最优方位偏转角度,红色虚线为该地震台站最终方位偏转角度(115.1°);(c) 地震台站B的纵波偏振分析结果,蓝色虚线表示最佳方位偏转角度,为108°
Figure 2. (a) A schematic map of coordinate system exhibiting the relationship between the geographical North direction,the North component and the back-azimuth (BAZ) of an event.
$ \chi $ and θ indicate the sensor orientation and the BAZ of the event relative to the North component,respectively. (b) Example of the Rayleigh-wave polarization analysis from the station B. The blue circles represent the optimal sensor orientations determined from each event based on the maximum correlation coefficient between the vertical and the Hilbert transformed radial components. The red dash line depicts the resulting station orientation (115.1°). (c) The P-wave polarization analysis for determining orientation of station B. The blue dashed line depicts the optimal sensor orientation,which is 108°1 台站B瑞雷波(a)与P波(b)径向分量波形随不同台站方位角的变化
瑞雷波和P波的波形数据分别来自于2010年3月4日和2010年4月30日发生的地震事件,图中红色波形对应最佳台站方位的情况
1. Waveform variations of Rayleigh wave (a) and P wave (b) in the radial component with sensor orientations
Data for Rayleigh wave and P wave are from events occuring on March 4,2010 and April 30,2010,respectively. The red waveforms in (a) and (b) represent the condition with the optimal sensor orientation
图 3 地震台站B方位偏转校正前(a)和后(b)的h-κ叠加结果
上图中黑色和红色线分别为每条接收函数和时间域所有接收函数简单叠加后的波形;下图左侧表示归一化后的h-κ叠加能量图,图中红点(为叠加能量最强点)表示最终确定的地壳厚度h和纵横波速比值κ (vP/vS),NRFs为接收函数数量。右下图表示h-κ叠加能量图中不同纵横波速比值所对应的最大能量点连线。(b)中蓝色线为校正前接收函数的简单叠加波形
Figure 3. h-κ stacking results from representative station B before (a) and after (b) misorientation correction
The upper panel shows individual (black trace) receiver functions (RFs) and a simple time-domain stack (red trace) of all RFs. The lower left panel illustrates the h-κ plot in which the maximum stacking amplitude (red dot) determines the optimal pair of crustal thickness h and vP/vS ratio κ. NRFs,number of receiver functions. The lower right panel displays the maxim stacking amplitude for each candidate vP/vS ratio in the h-κ plot. The blue trace in (b) represents the simple stacked RF trace in time domain before misorientation correction
图 4 瑞雷波和纵波偏振分析得到的OBS方位结果对比图(a)以及本文和前人的地壳厚度结果(Abers et al,2016)对比图(b)。其中*表示前人研究的结果。
Figure 4. Comparisons of OBS orientations from the Rayleigh-wave and P-wave polarization analyses (a) and crustal thickness from this and previous (indicated by *) studies (Abers et al,2016), respectively
6 地震台站D纵波偏振分析方位偏转校正前(a)后(b)的h-κ叠加结果对比图
上图中黑色和红色线分别为每条接收函数和时间域所有接收函数简单叠加后的波形。下图左侧表示归一化后的h-κ叠加能量图,途中红点(为叠加能量最强点)表示最终确定的地壳厚度(Depth,缩写为Dep)和纵横波速比值(vP/vS),NRFs为接收函数数量。右下图表示h-κ叠加能量图中不同纵横波速比值随对应的最大能量点连线
6. h-κ stacking results from station D before (a) and after (b) misorientation correction from the P-wave polarization analysis
The upper panel shows individual (black trace) receiver functions (RFs) and a simple time-domain stack (red trace) of all RFs. The lower left panel illustrates the h-κ plot in which the maximum stacking amplitude (red dot) determines the optimal pair of crustal thickness and vP/vS ratio. The lower right panel displays the maxim stacking amplitude for each candidate vP/vS ratio in the h-κ plot
图 5 方位校正后地壳的平面(a,b)与剖面(c,d)结果图
(a) 地壳厚度平面分布图;(b) 纵横波速比值平面分布图;(c) 地壳厚度(蓝色三角号)沿图(a)中AA′测线的剖面结果图,其中上图为地形起伏,下图红色波形为时深转换后的叠加接收函数,KB:基里比斯盆地,GB:古迪纳夫盆地,DI:当特尔卡斯托群岛);(d)地壳纵横波速比值(红色圆圈)沿(b)中AA′测线的剖面结果图
Figure 5. Planar (a,b) and vertical (c,d) display of the resulting crustal measurements after misorientation corrections
(a,b) Crustal thickness and vP/vS ratio in planar view,respectively;(c) Vertical display of crustal thickness (blue triangles) along the AA′ profile in (a). The upper panel shows the topography. Red traces in the lower panel are the stacked RF traces after time-depth conversion. KB (Kiribisi Basin),GB (Goodenough Basin),DI (D’ Entrecasteaux Island);(d) Vertical display of crustal vP/vS ratio (red circles) along the AA′ profile in (b)
4 瑞雷波偏振分析得到的地震台站D,E和F的方位偏转角度
蓝色圆圈为通过每个地震的垂直分量和希尔伯特变换后的径向分量之间的最大相关系数确定的对应最优方位偏转角度,红色虚线为该地震台站的最终方位偏转角度
4. The Rayleigh-wave polarization analysis for determining orientations of stations D,E and F,respectively
The blue dots represent the optimal sensor orientations determined from each event based on the maximum correlation coefficient between the vertical and the Hilbert transformed radial components. The red dash line depicts the resulting station orientation
表 1 纵波和瑞雷波偏振分析得到的每台OBS方位偏转角χ
Table 1. The resulting OBS orientations χ from analysis of the Rayleigh-wave and P-wave polarization analyses for each station
台站 瑞雷波偏振分析 纵波偏振分析 χ*/° χ/° STD/° χ/° STD/° B 115.1 0.55 108 3.03 116.8 D 3.0 1.48 356 2.68 0.3 E 176.0 1.65 168 3.17 183.1 F 293.1 1.78 310 5.26 301.0 G 44.2 1.92 29 7.68 47.2 H 323.0 1.88 319 1.12 324.5 J 131.8 2.53 144 12.97 143.9 *STD为标准差;χ*为Eilon等(2014)基于瑞雷波偏振分析所得。 表 2 本文和前人获得的地壳厚度h和纵横波速比值κ
Table 2. Crustal thickness h and vP/vS ratios κ from this and previous studies
台站 西经/° 北纬/° h /km h*/km κ B 9.749 150.350 30.4±0.21 27.8±0.53 1.94±0.028 D 9.943 150.707 36.7±0.16 44.9±0.73 1.89±0.005 E 10.080 150.621 34.2±2.69 1.82±0.059 F 9.950 150.200 31.7±0.96 36.2±0.61 1.90±0.072 G 9.333 149.667 25.8±0.73 31.5±0.84 2.04±0.039 H 9.000 149.668 16.8±0.71 36.8±0.74 1.94±0.085 J 8.870 149.934 29.8±0.60 39.2±0.38 1.99±0.043 * h*为Abers等(2016)的结果。 1 不同后方位角(Baz)组分别求得的方位偏转角度(
$ \chi $ )1. Resulting sensor orientations (χ) for different backazimuth (Baz) ranges
台站 Baz:0°— 90° Baz:90°—180° Baz:180°—270° Baz:270°—360° χ/° N χ/° N χ/° N χ/° N B 118 2 101 4 110 17 D 8 2 7 4 354 19 E 168 1 169 2 168 18 F 305 1 300 1 300 1 311 11 G 19 1 19 6 39 1 30 15 H 309 1 309 7 321 15 J 134 3 148 11 注:N表示地震事件的个数。 2 瑞雷波偏振分析方位校正前后的地壳厚度h和纵横波速比κ的值
2. Observations of crustal thickness h and Vp/Vs ratios κ before and after misorientation corrections from the Rayleigh-wave polarization (R-pol) analysis
台站 西经/° 北纬/° h /km h*/km κ κ* B 9.749 150.350 21.5±0.62 30.3±0.20 1.85±0.026 1.95±0.010 D 9.943 150.707 36.7±0.59 36.7±0.34 1.89±0.049 1.89±0.022 E 10.080 150.621 32.1±1.56 24.5±6.91 1.90±0.025 2.14±0.073 F 9.950 150.200 27.3±0.69 29.7±1.90 2.11±0.47 2.01±0.161 G 9.333 149.667 27.5±0.15 24.4±1.09 1.91±0.013 2.14±0.031 H 9.000 149.668 22.8±0.67 26.4±2.27 1.74±0.081 1.65±0.028 J 8.870 149.934 28.2±1.76 29.8±0.49 1.71±0.073 2.02±0.036 *Results after misorientation corrections from the R-pol analysis
*瑞雷波偏振分析校正后得到的地壳结果。 -
胡昊,阮爱国,游庆瑜,李家彪,郝天珧,龙江平. 2016. 海底地震仪远震记录接收函数反演v:以南海西南次海盆为例[J]. 地球物理学报,59(4):1426–1434. doi: 10.6038/cjg20160423 Hu H,Ruan A G,You Q Y,Li J B,Hao T Y,Long J P. 2016. Using OBS teleseismic receiver functions to invert lithospheric structure:A case study of the southwestern subbasin in the South China Sea[J]. Chinese Journal of Geophysics,59(4):1426–1434 (in Chinese). 黄海波,丘学林,胥颐,曾钢平. 2011. 利用远震接收函数方法研究南海西沙群岛下方地壳结构[J]. 地球物理学报,54(11):2788–2798. doi: 10.3969/j.issn.0001-5733.2011.11.009 Huang H B,Qiu X L,Xu Y,Zeng G P. 2011. Crustal structure beneath the Xisha Islands of the South China Sea simulated by the teleseismic receiver function method[J]. Chinese Journal of Geophysics,54(11):2788–2798 (in Chinese). 丘学林,曾钢平,胥颐,郝天珧,李志雄,Priestley K,McKenzie D. 2006. 南海西沙石岛地震台下的地壳结构研究[J]. 地球物理学报,49(6):1720–1729. doi: 10.3321/j.issn:0001-5733.2006.06.019 Qiu X L,Zeng G P,Xu Y,Hao T Y,Li Z X,Priestley K,McKenzie D. 2006. The crustal structure beneath the Shidao station on Xisha Islands of South China Sea[J]. Chinese Journal of Geophysics,49(6):1720–1729 (in Chinese). Abers G A,Ferris A,Craig M,Davies H,Lerner-Lam A L,Mutter J C,Taylor B. 2002. Mantle compensation of active metamorphic core complexes at Woodlark rift in Papua New Guinea[J]. Nature,418(6900):862–865. doi: 10.1038/nature00990 Abers G A,Eilon Z,Gaherty J B,Jin G,Kim Y,Obrebski M,Dieck C. 2016. Southeast Papuan crustal tectonics:Imaging extension and buoyancy of an active rift[J]. J Geophys Res:Solid Earth,121(2):951–971. doi: 10.1002/2015JB012621 Ammon C J. 1991. The isolation of receiver effects from teleseismic P waveforms[J]. Bull Seismol Soc Am,81(6):2504–2510. doi: 10.1785/BSSA0810062504 Baker G E,Stevens J L. 2004. Backazimuth estimation reliability using surface wave polarization[J]. Geophys Res Lett,31(9):L09611. Baldwin S L,Lister G S,Hill E J,Foster D A,McDougall I. 1993. Thermochronologic constraints on the tectonic evolution of active metamorphic core complexes,D'Entrecasteaux Islands,Papua New Guinea[J]. Tectonics,12(3):611–628. doi: 10.1029/93TC00235 Baldwin S L,Monteleone B D,Webb L E,Fitzgerald P G,Grove M,June Hill E. 2004. Pliocene eclogite exhumation at plate tectonic rates in eastern Papua New Guinea[J]. Nature,431(7006):263–267. doi: 10.1038/nature02846 Biemiller J,Boulton C,Wallace L,Ellis S,Little T,Mizera M,Niemeijer A,Lavier L. 2020. Mechanical implications of creep and partial coupling on the world's fastest slipping low-angle normal fault in Southeastern Papua New Guinea[J]. J Geophys Res:Solid Earth,125(10):e2020JB020117. Bird P. 2003. An updated digital model of plate boundaries[J]. Geochem, Geophys, Geosyst, 4(3): 1027. Davies H L,Smith I E. 1971. Geology of Eastern Papua[J]. GSA Bull,82(12):3299–3312. doi: 10.1130/0016-7606(1971)82[3299:GOEP]2.0.CO;2 Davies H L,Jaques A L. 1984. Emplacement of ophiolite in Papua New Guinea[J]. Geol Soc Lond Spec Publ,13(1):341–349. doi: 10.1144/GSL.SP.1984.013.01.27 Efron B,Tibshirani R. 1986. Bootstrap methods for standard errors,confidence intervals,and other measures of statistical accuracy[J]. Statist Sci,1(1):54–75. Eilon Z, Abers G A, Jin G, Gaherty J B. 2014. Anisotropy beneath a highly extended continental rift[J]. Geochem, Geophys, Geosyst, 15(3): 545-564. Eilon Z, Abers G A, Gaherty J B, Jin G. 2015. Imaging continental breakup using teleseismic body waves: The Woodlark Rift, Papua New Guinea[J]. Geochem, Geophys, Geosyst, 16(8): 2529-2548. Ellis,S. M.,Little T. A.,Wallace L. M.,Hacker B. R.,Buiter S. J. H 2011. Feedback between rifting and diapirism can exhume ultrahigh-pressure rocks[J]. Earth Planet. Sci. Lett.,311:427–438. doi: 10.1016/j.jpgl.2011.09.031 Ferris A,Abers G A,Zelt B,Taylor B,Roecker S. 2006. Crustal structure across the transition from rifting to spreading:The Woodlark rift system of Papua New Guinea[J]. Geophys J Int,166(2):622–634. doi: 10.1111/j.1365-246X.2006.02970.x Finlayson D M,Drummond B J,Collins C D M,Connelly J B. 1977. Crustal structures in the region of the Papuan Ultramafic Belt[J]. Phys Earth Planet Inter,14(1):13–29. doi: 10.1016/0031-9201(77)90043-7 Finlayson D M,Muirhead K J,Webb J P,Gibson G,Furumoto A S,Cooke R J S,Russell A J. 1976. Seismic investigation of the Papuan Ultramafic Belt[J]. Geophys J R Astron Soc,44(1):45–59. doi: 10.1111/j.1365-246X.1976.tb00274.x Fitz G, Mann P. 2013. Evaluating upper versus lower crustal extension through structural reconstructions and subsidence analysis of basins adjacent to the D'Entrecasteaux Islands, eastern Papua New Guinea[J]. Geochem, Geophys, Geosyst, 14(6): 1800-1818. Gordon S M,Little T A,Hacker B R,Bowring S A,Korchinski M,Baldwin S L,Kylander-Clark A R C. 2012. Multi-stage exhumation of young UHP-HP rocks:Timescales of melt crystallization in the D'Entrecasteaux Islands,southeastern Papua New Guinea[J]. Earth Planet Sci Lett,351-352:237–246. doi: 10.1016/j.jpgl.2012.07.014 Hill E J,Baldwin S L,Lister G S. 1992. Unroofing of active metamorphic core complexes in the D'Entrecasteaux Islands,Papua New Guinea[J]. Geology,20(10):907–910. doi: 10.1130/0091-7613(1992)020<0907:UOAMCC>2.3.CO;2 Holbrook W S, Mooney W D, Christensen N I. 1992. The seismic velocity structure of the deep continental crust[G]//Continental Lower Crust. New York: Elsevier: 21–43. Hung T D,Yang T,Le B M,Yu Y Q,Xue M,Liu B H,Liu C G,Wang J,Pan M H,Huong P T,Liu F,Morgan J P. 2021. Crustal structure across the extinct mid-ocean ridge in South China Sea from OBS receiver functions:Insights into the spreading rate and magma supply prior to the ridge cessation[J]. Geophys Res Lett,48(3):e2020GL089755. Jin G,Gaherty J B,Abers G A,Kim Y,Eilon Z,Buck W R. 2015. Crust and upper mantle structure associated with extension in the Woodlark Rift,Papua New Guinea from Rayleigh-wave tomography[J]. Geochem,Geophys,Geosyst,16(11):3808–3824. Kennett B L N,Engdahl E R. 1991. Traveltimes for global earthquake location and phase identification[J]. Geophys J Int,105(2):429–465. doi: 10.1111/j.1365-246X.1991.tb06724.x Laske G. 1995. Global observation of off-great-circle propagation of Long-Period surface waves[J]. Geophys J Int,123(1):245–259. doi: 10.1111/j.1365-246X.1995.tb06673.x Lim H,Kim Y,Song T R A,Shen X Z. 2018. Measurement of seismometer orientation using the tangential P-wave receiver function based on harmonic decomposition[J]. Geophys J Int,212(3):1747–1765. doi: 10.1093/gji/ggx515 Little T A,Hacker B R,Gordon S M,Baldwin S L,Fitzgerald P G,Ellis S,Korchinski M. 2011. Diapiric exhumation of Earth's youngest (UHP) eclogites in the gneiss domes of the D'Entrecasteaux Islands,Papua New Guinea[J]. Tectonophysics,510(1/2):39–68. Liu K H, Gao S S. 2010. Spatial variations of crustal characteristics beneath the Hoggar swell, Algeria, revealed by systematic analyses of receiver functions from a single seismic station[J]. Geochem, Geophys, Geosyst, 11(8): Q08011. Lus W Y,McDougall I,Davies H L. 2004. Age of the metamorphic sole of the Papuan Ultramafic Belt ophiolite,Papua New Guinea[J]. Tectonophysics,392(1/2/3/4):85–101. Martinez F,Goodliffe A M,Taylor B. 2001. Metamorphic core complex formation by density inversion and lower-crust extrusion[J]. Nature,411(6840):930–934. doi: 10.1038/35082042 Monteleone B D,Baldwin S L,Webb L E,Fitzgerald P G,Grove M,Schmitt A K. 2007. Late Miocene-Pliocene eclogite facies metamorphism,D'Entrecasteaux Islands,SE Papua New Guinea[J]. J Metamorph Geol,25(2):245–265. doi: 10.1111/j.1525-1314.2006.00685.x Niu F L,Li J. 2011. Component azimuths of the CEArray stations estimated from P-wave particle motion[J]. Earthq Sci,24(1):3–13. doi: 10.1007/s11589-011-0764-8 Selby N D. 2001. Association of Rayleigh waves using backazimuth measurements:Application to test ban verification[J]. Bull Seismol Soc Am,91(3):580–593. doi: 10.1785/0120000068 Stachnik J C,Sheehan A F,Zietlow D W,Yang Z H,Collins J,Ferris A. 2012. Determination of New Zealand ocean bottom seismometer orientation via Rayleigh-wave polarization[J]. Seismol Res Lett,83(4):704–713. doi: 10.1785/0220110128 Tregoning P,Lambeck K,Stolz A,Morgan P,Mcclusky S C,van der Beek P,Mcqueen H,Jackson R J,Little R P,Laing A,Murphy B. 1998. Estimation of current plate motions in Papua New Guinea from Global Positioning System observations[J]. J Geophys Res:Solid Earth,103(B6):12181–12203. doi: 10.1029/97JB03676 van Ufford A Q,Cloos M. 2005. Cenozoic tectonics of New Guinea[J]. AAPG Bull,89(1):119–140. doi: 10.1306/08300403073 Wallace L M,Stevens C,Silver E,McCaffrey R,Loratung W,Hasiata S,Stanaway R,Curley R,Rosa R,Taugaloidi J. 2004. GPS and seismological constraints on active tectonics and arc-continent collision in Papua New Guinea:Implications for mechanics of microplate rotations in a plate boundary zone[J]. J Geophys Res:Solid Earth,109(B5):B05404. Watanabe T. 1993. Effects of water and melt on seismic velocities and their application to characterization of seismic reflectors[J]. Geophys Res Lett,20(24):2933–2936. doi: 10.1029/93GL03170 Webb L E,Baldwin S L,Little T A,Fitzgerald P G. 2008. Can microplate rotation drive subduction inversion?[J]. Geology,36(10):823–826. doi: 10.1130/G25134A.1 Yu Y Q,Tilmann F,Zhao D P,Gao S S,Liu K H. 2022. Continental break-up under a convergent setting:Insights from P wave radial anisotropy tomography of the Woodlark rift in Papua New Guinea[J]. Geophys Res Lett,49(5):e2022GL098086. Zelt B C,Taylor B,Goodliffe A M. 2001. 3-D crustal velocity structure at the rift tip in the western Woodlark Basin[J]. Geophys Res Lett,28(15):3015–3018. doi: 10.1029/2001GL012864 Zeng S J,Zheng Y,Niu F L,Ai S X. 2021. Measurements of seismometer orientation of the first phase CHINArray and their implications on vector-recording-based seismic studies[J]. Bull Seismol Soc Am,111(1):36–49. doi: 10.1785/0120200129 Zha Y,Webb S C,Menke W. 2013. Determining the orientations of ocean bottom seismometers using ambient noise correlation[J]. Geophys Res Lett,40(14):3585–3590. doi: 10.1002/grl.50698 Zheng Y F. 2009. Fluid regime in continental subduction zones:Petrological insights from ultrahigh-pressure metamorphic rocks[J]. J Geol Soc,166(4):763–782. doi: 10.1144/0016-76492008-016R Zhu G H,Yang H F,Lin J,You Q Y. 2020. Determining the orientation of ocean-bottom seismometers on the seafloor and correcting for polarity flipping via polarization analysis and waveform modeling[J]. Seismol Res Lett,91(2A):814–825. doi: 10.1785/0220190239 Zhu L P,Kanamori H. 2000. Moho depth variation in southern California from teleseismic receiver functions[J]. J Geophys Res:Solid Earth,105(B2):2969–2980. doi: 10.1029/1999JB900322 -