SEDIMENTARY RECORDS OF VOLCANIC ACTIVITIES IN THE SOUTH CHINA SEA OVER THE PAST 480 ka
-
摘要: 通过对南海西部越南岸外MD05-2901孔的研究,建立了48万年来高分辨率的氧同位素年代地层。针对磁化率结果显示的13层厚度超过4 cm的磁化率峰值层,将其年代和菲律宾巴拉望岛岸外的MD972142孔的火山灰层进行对比,发现其中10层完全对应,证明磁化率峰值层是由火山喷发作用造成的,说明48万年来菲律宾火山带喷发的粉尘可以沉积到南海西部。磁化率峰值层样品的粒度特征明显不同,呈双峰态分布,粗颗粒组分较多,空隙率高,密度和颜色反射率较小。由此证明,磁化率和粒度分析可用于寻找肉眼不可见的火山喷发作用的沉积学记录。研究发现的受火山喷发作用影响的特殊层位,可作为地层对比的参比标志,并为追溯48万年来大气和洋流的变迁提供了参考资料。Abstract: Magnetic susceptibility (MS) and grain size distribution of core MD05-2901, which is located in the upwelling area off eastern Vietnam in the South China Sea, were measured. Oxygen isotope stratigraphy indicates that the core spans a time period of 480 ka. Magnetic susceptibility curve displays glacial-interglacial variations in general, but with prominent peaks. Ten out of thirteen MS peak layers, which are thicker than 4 cm, correlate well with the volcanic ash layers recognized in core MD972142 offshore northwestern Palawan. The result indicates that the volcanic material from the Philippine seismic zone could have reached the western part of the South China Sea during the last 480 ka. The content of size fraction>16μm is relatively high in the layers of high MS values. In addition, grain size distribution of these layers shows distinct double-mode, indicating input from volcanic activities. MS and grain size features of the sediments can be used as an indicator of volcanic activities, which otherwise can't be examined easily by naked-eye. The MS peaks are also a useful tool for chronology and stratigraphic correlations.
-
Keywords:
- volcanic ash /
- magnetic susceptibility /
- grain size /
- late Quaternary /
- South China Sea
-
重磁资料是揭示地质构造信息的重要依据，在利用地球物理资料进行地质解释过程中，重磁数据的处理是一个重要的环节。重力异常是由地下不同规模、不同形态和不同埋深的不均匀地质体的重力作用叠加而成[1]，通常将实测获取的重力异常数据看成是区域异常场和局部异常场的组合[2]。区域场与局部场是一个相对的概念，二者的划分可由所研究的对象不同而变化[3]。区域异常的分布范围大，幅度大，变化平缓，一般由分布广的中深地质因素引起；局部异常的分布有限，幅度小，梯度大，一般由埋深浅、体积小的地质因素引起。依据勘探目标的不同，需要将实测异常数据中区域性异常分离出来，圈定局部异常，从而更好地揭示与矿产、构造等相关地质现象的关系，也能对分离后的异常数据进行更好的定量解释。因此，将重力数据中的区域场和局部场分离对于重力勘探显得尤为重要。
异常分离是一个很早就出现的需求，在20世纪50、60年代，用于分离异常的数字处理方法已经提出，经过几十年来人们不断的探索，异常分离的方法越来越成熟，并且出现了很多新方法[4-13]。这些方法由于数学原理不同，应用前提也不尽相同，具有针对性和选择性[1]，在进行重力异常分离时应根据需求选取合适的方法。
插值切割法在程方道、文百红等[10-17]的研究基础上不断发展成熟，该方法在分离重磁异常数据中具有抗干扰、计算速度快、收敛性好、分辨率和精度较高的特点，是一种非常实用的异常分离方法，并在实践中取得了很多重要的应用。基于上述优点，本文选取插值切割法来对南海的重力场进行分离。
利用插值切割法分离重力场，确定切割半径是一个关键。文百红[13]指出，只有选择了最佳的切割半径，才能取得最佳的分离效果。张虹[17]指出，切割半径及切割精度一定要合理选取，数据间距要小，保证插值切割法发挥作用。由于地质现象的复杂性，所形成的重力异常在幅度变化、波长长短、叠加程度、干扰水平等方面也是错综复杂，因此，在确定最佳切割半径时，必须同时考虑异常和干扰的空间分布特征。一般是根据实际情况，选取一系列的切割半径，通过对比不同切割半径下的异常场分离效果来确定最佳的切割半径。如何从定量上分析确定一个最佳的切割半径，前人很少提及。本文以南沙群岛中北部的一块岛礁区域为例，通过计算比较，认为重力局部异常、区域异常与原始数据及地形数据的相关性系数可以作为确定最佳切割半径的依据，并将确定好的最佳半径应用于整个南海重力场的切割分离。
1. 插值切割法原理
插值切割法的基本原理是构造一个切割算子[13]，重复作用于测点的场值，迭代切割使其收敛达到稳定值，得到区域场，从观测场中减去区域场即得到局部场。
针对矩形网格化的离散数据，汪炳柱等[13]给出了一个多次插值切割法的计算实例，具体算法如下：
设一个测区中某点的重力异常值为G(i, j)，区域异常为R(i, j)，局部异常为L(i, j)，则有
$$ G\left( {i, j} \right) = R\left( {i, j} \right) + L\left( {i, j} \right) $$ (1) 其中(i, j)表示测点的行列坐标(i=0, 1, 2, …, M；j=0, 1, 2, …, N；M+1是x方向的点数；N+1是y方向的点数。)
计算点(i, j)的区域异常是关键的一步，该点的R(i, j)是其重力异常值与其周围4点的重力异常值的加权平均，用(2)式计算。
$$ R\left( {i, j} \right) = \left[{1-a\left( {i, j} \right)/2} \right]B\left( {i, j} \right) + \left[{a\left( {i, j} \right)/2} \right]G\left( {i, j} \right) $$ (2) 公式(2)中a(i, j)为一个加权系数。B(i, j)为4点的重力异常平均值，用公式(3)计算。
$$ \begin{array}{l} B\left( {i, j} \right) = 0.25\left[{G\left( {i + {n_x}, j} \right) + G\left( {i-{n_x}, j} \right) + } \right.\\ \left. {G\left( {i, j + {n_y}} \right) + G\left( {i, j-{n_y}} \right)} \right] \end{array} $$ (3) 这里nx、ny分别为x方向和y方向的切割半径，通常取相同值，为点距的整数倍。
加权系数a(i, j)由公式(4)计算。
$$ a\left( {i, j} \right) = b\left( {i, j} \right) + c\left( {i, j} \right) $$ (4) $$ 其中, b\left( {i, j} \right) = 1/\left\{ {\left[{{F_x}{{\left( {i, j} \right)}^2} + 1} \right]} \right\} $$ (5) 这里Fx(i, j)是x方向的非线性度，由公式(6)—(8)计算。
$$ {F_x}\left( {i, j} \right) = \Delta {B_{x1}}\left( {i, j} \right)/\Delta {B_{x2}}\left( {i, j} \right) $$ (6) $$ \Delta {B_{x1}}\left( {i, j} \right) = G\left( {i, j} \right)- 0.5\left[{G\left( {i + {n_x}, j} \right) + G\left( {i-{n_x}, j} \right)} \right] $$ (7) $$ \Delta {B_{x2}}\left( {i, j} \right) = G\left( {i + {n_x}, j} \right)-G\left( {i-{n_x}, j} \right) $$ (8) 当Bx2(i, j)=0时，若Bx1(i, j)=0，则b(i, j)=1；若Bx1(i, j)≠0，则b(i, j)=0。
c(i, j)和b(i, j)的意义相同，由公式(9)—(12)计算。
$$ c\left( {i, j} \right) = 1/\left\{ {{{\left[{{F_y}\left( {i, j} \right)} \right]}^2} + 1} \right\} $$ (9) 这里Fy(i, j)是y方向的非线性度。
$$ {F_y}\left( {i, j} \right) = \Delta {B_{y1}}\left( {i, j} \right)/\Delta {B_{y2}}\left( {i, j} \right) $$ (10) $$ \Delta {B_{y1}}\left( {i, j} \right) = G\left( {i, j} \right)- 0.5\left[{G\left( {i, j + {n_y}} \right) + G\left( {i, j-{n_y}} \right)} \right] $$ (11) $$ \Delta {B_{y2}}\left( {i, j} \right) = G\left( {i, j + {n_y}} \right)-G\left( {i, j-{n_y}} \right) $$ (12) 当By2(i, j)＝0时，若By1(i, j)＝0，则c(i, j)=1；若By1(i, j)≠0，则c(i, j)=0。
在计算区域异常时，是重复利用公式(2)进行多次插值切割的。设Rk(i, j)，Rk-1(i, j)分别是经过k次和k－1次插值切割得到的点(i, j)的区域异常，当k趋向于∞时，max | Rk(i, j)－Rk-1(i, j) |=0。令V=max | Rk(i, j)-Rk-1(i, j) |，取e为一个很小的正数，若V>e，则令G(i, j)=Rk(i, j)，进行下一次迭代计算；若V≤e，则可将Rk(i, j)作为点(i, j)的区域异常。
将式(1)变换为：L(i, j)=G(i, j)－R(i, j)，即可求出局部异常L(i, j)。
切割半径是一个重要参数，它可以控制划分出不同层次的局部场和区域场。切割半径越大，切割出的局部场反映的地质体的深度和规模越大。
2. 切割半径的确定
本文采用的南海重力数据来源于1:200万“南海地质地球物理图系——空间重力异常图”[21]，该重力数据的主体数据来源于我国的海洋重力实际测量，以广州海洋地质调查局历年调查成果为主，还包括国家海洋局及国外在南海的部分调查成果，空白及资料不足的地方用卫星测量数据补充[18-23]。
陈洁等[18]指出南海的大部分海岛海山等正向高地均与重力高重合。利用南海双子岛礁至九章群礁一带的地形数据和不同切割半径下分离的重力局部异常和区域异常数据，计算它们的相关性，可以从定量上对重力局部异常、区域异常与原始数据及海底地形的对应关系做出探讨，从而确定一个最佳的切割半径来应用于整个南海的重力切割。
图 1为南沙群岛中北部的海底地形，分布着双子群礁、中业群礁、道明群礁、郑和群礁和九章群礁等岛礁以及礼乐西海槽等，水深变化超过4000m。图 2为对应的空间重力异常图，可见重力异常的高低值与岛礁、沟槽有着明显的对应。
![]()
图 1 南沙群岛中北部地形图(引自南海地质地球物理图系)Figure 1. Topographic map of the Nansha Islands (partial)![]()
图 2 南沙群岛中北部空间重力异常图Figure 2. Free-air gravity anomaly map of the Nansha Islands (partial)以一系列不同的切割半径对该区的重力进行局部异常和区域异常的分离，计算这些异常数据与未切割的原始数据及地形数据的相关性系数，如图 3、图 4所示。其中，2~10km之间以2km为间隔，10~100km之间以10km为间隔。未切割的原始重力异常数据与地形的相关性系数为0.692，重力异常与地形有极显著的相关性。
![]()
图 3 重力局部异常及区域异常与原始数据的相关性折线图Figure 3. Correlation of local and regional gravity anomalies with original data![]()
图 4 重力局部异常及区域异常与地形的相关性折线图Figure 4. Correlation of local and regional gravity anomalies with terrain考察局部异常与原始重力数据的相关性(图 3)，当切割半径增大时，相关性系数不断增大。切割半径在2~20km时增长很快；从30km到100km，相关性系数基本保持不变，维持在0.98左右的高位，可以认为切割半径大于30km之后，对分离出局部异常的作用不大。重力与地形密切相关，随切割半径增大，局部异常与地形的相关性系数也不断增大(图 4)。切割半径在2~30km时增长很快，到40km以后，基本保持平稳。
区域异常的相关性系数总体变化与局部异常是相反的，图 3中与原始数据的相关性系数随切割半径的增加而下降，在70km时略有反转，在10km时与局部异常的折线相交；图 4中与地形数据的相关性系数也是随切割半径的增加而下降，在2~30km为0.7~0.76，随后快速下降到50km处的0.493，之后保持稳定。
局部异常反映浅层地质因素的影响，区域异常反映中深部地质因素的影响。从分离角度来说，局部异常与原始数据及地形数据的相关性越低越好，但是切割半径小，分离得到的局部异常规模太小，而相应的区域异常不能反映深部的东西，如图 3中切割半径为2、4km时，区域异常与原始数据的相关性很高，为1和0.993。因此需要在这两者之间权衡。
图 5列出了切割半径为6~30km的局部异常及区域异常分布图，不同的切割半径下均能较好地分离出局部异常及区域异常。切割半径为20、30km的区域异常特征不够明显，切割半径为10km的局部异常比切割半径为6、8km的细节更丰富。综合考虑，可以确定10km为最佳的切割半径，用于对整个南海的重力切割分离。
![]()
图 5 不同切割半径下南沙群岛重力局部异常和区域异常图Figure 5. Local and regional anomalies of the Nansha Islands in different cutting radius3. 南海重力异常分离
前人的研究中，将南海的重力格局划分为北部陆架-陆坡异常区、中-西沙岛礁异常区、海盆异常区、南沙岛礁异常区、西部陆架-陆坡异常区和东部岛弧-海沟异常区[24-26]等，并通过重磁分析探讨南海构造及演化问题[27-30]。
结合之前的分析，以10km作为切割半径对重力数据进行异常分离。图 6为南海空间重力异常图，图 7为以10km切割半径插值切割分离得到的区域异常图，图 8为对应的局部异常图。从插值切割结果中可以更清晰地看出南海重力异常的分布特征。
![]()
图 7 南海空间重力区域异常图(切割半径10km)Figure 7. Regional gravity anomaly map of the South China Sea (cutting radius of 10km)![]()
图 8 南海空间重力局部异常图(切割半径10km)Figure 8. Local gravity anomaly map of the South China Sea (cutting radius of 10km)图 7为区域异常场，去掉了浅层重力异常的影响，某些岛礁异常、岛弧-海沟异常等均不同程度被抹平，只留下规模巨大、分布宽缓的背景场，从而使得该切割半径下的区域异常反映分布更广更深的地质体的特征。区域异常场与图 6的空间重力异常分布基本一致。区域异常数值集中于-40~40mGal，总体显示为北低南高；南海北部的低异常带对应着盆地构造的凹陷区，如莺歌海盆地；台湾岛与吕宋岛之间的巴士海峡及吕宋岛西缘的重力梯度带反映了俯冲带的特征；中沙群岛、南沙群岛及南部的南沙海槽均有明显的展示。
图 8的局部异常中，将分布广、变化缓的-20~20mGal的局部异常标示为白色，便于突出显示。图中，巴士海峡、中部西沙群岛、中沙群岛及黄岩岛等海域的局部异常特征明显，南沙海域10°至12°之间正负相间的异常特征得到了突出显示。
4. 结语
插值切割法是一种非常实用的异常分离方法，在分离重力异常的应用中可以取得较好的效果。该方法中，确定最佳的切割半径是一个重要的步骤。本文通过计算重力局部异常、区域异常与原始数据及地形数据的相关性系数，发现重力异常与地形有极显著的相关性，从相关性结果可以初步确定理想的切割半径；并对比不同切割半径下的分离效果，可以很好地确定最佳切割半径。以10km最佳半径应用于整个南海重力场的切割分离，取得很好的分离效果。
重力异常的插值切割处理只是单纯的分离，要获取更深的认识，还应结合已知的地质资料和其他的地球物理资料来进行解释。
-
[1] 冯文科,薛万俊,杨达源. 南海北部晚第四纪地质环境[M]. 广州:广东科技出版社,1988:156-162.[Feng Wen-ke, XUE Wan-jun, YANG Da-yuan. The Geological Environment of Late Quaternary in the Northern South China Sea[M]. Guangzhou:Guangdong Science and Technology Publishing House, 1988:156 -162.]
[2] Kennett J P. Marine tephrochronology[C]//The Sea, Vol. 7. Wiley, New York, 1981:1373-1436.
[3] Oppenheimer Clive. Limited global change due to the largest known Quaternary eruption, Toba≈74kyr BP?[J].Quaternary Science Reviews, 2002, 21:1593-1609.
[4] Wiesner M G, Wang Y, Zheng L. Fallout of volcanic ash to the deep South China Sea induced by the 1991 eruption of Mount Pinatubo (Philippines)[J]. Geology, 1995, 23:885-888.
[5] Wiesner M G, Wetzel A, Catane S G, et al. Grain size, area thickness distribution and controls on sedimentation of the 1991 Mount Pinatubo tephra layer in the South China Sea[J]. Bulletin of Volcanology, 2004, 66:226-242.
[6] 杨育标,范时清. 南海深海晚第四纪火山沉积物及其起源探讨[J]. 热带海洋,1990,9(1):52-59. [YANG Yu-biao, FAN Shi-qing. Research on volcanic sediments and origion of volcanic substance in South China Sea during late Quaternary[J]. Tropical Oceanography, 1990, 9(1):52-59.]
[7] 杨群慧,林振宏,张富元,等. 南海东部重矿物分布特征及其影响因素[J]. 青岛海洋大学学报自然科学版,2002,32(6):956-964. [YANG Qun-hui, LIN Zhen-hong, ZHANG Fu-yuan, et al. The distribution characteristics of heavy minerals in the east of South China Sea and their controlling factors[J]. Journal of Ocean University of Qingdao, 2002, 32(6):956-964.]
[8] Bühring C, Sarnthein M, Leg 184 Shipboard Scientific Party. Toba ash layers in the South China Sea:evidence of contrasting wind directions during eruption ca. 74ka[J]. Geology, 2000, 28:275-278.
[9] Song S R, Chen C H, Lee M Y, et al. Newly discovered eastern dispersal of the youngest Toba Tuff[J]. Marine Geology, 2000, 167:303-312.
[10] Liang X, Wei G, Shao L, et al. Records of Toba eruptions in the South China Sea:chemical characteristics of the glass shards from ODP 1143A[J]. Science in China (Series D), 2001, 44:871-878.
[11] Liu Zhi-fei, Colin C, Trentesaux A. Major element geochemistry of glass shards and minerals of the Youngest Toba Tephra in the southwestern South China Sea[J]. Journal of Asian Earth Sciences, 2006, 27:99-107.
[12] 陈文斌,周福根. 南海北部晚第四纪火山玻璃的分布及其地层学意义[C]//南海海洋沉积化学过程与地球化学研究. 北京:海洋出版社, 1993:39-45.[CHEN Wen-bin, ZHOU Fu-gen. A study of volcanic glass in northern South China Sea during the last 100 ka[C]//Contributions of Sedimentation Process and Geochemistry of the South China Sea. Beijing:China Ocean Press, 1993 :39-45.]
[13] 王慧中,周福根,翦知湣. 中沙台缘碳酸盐软泥中的火山碎屑及其古环境意义[C]//南海晚第四纪古海洋学研究.青岛:青岛海洋大学出版社,1992:42-55.[WANG Hui-zhong, ZHOU Fu-gen, JIAN Zhi-min. Volcanic clasts in the periplatform carbonate ooze near Zhongsha Islands and their bearing on paleo-environment[C]//Contributions to Late Quaternary Paleoceanography of the South China Sea. Qingdao:Qingdao Ocean University Press, 1992:42 -55.]
[14] Wei Kuo-Yen, Lee Teh-Quei, et al. Late Pleistocene volcanic ash layers in Core MD972142, offshore from northwestern Palawan, South China Sea:A preliminary report[J]. TAO, 1998, 9(1):143-152.
[15] Wei Kuo-Yen, Chiu Tzu-Chien, Chen Yue-Gau. Toward establishing a maritime proxy record of the East Asian summer monsoons for the Late Quaternary[J]. Marine Geology, 2003, 201:67-79.
[16] Ku Y P, Chen C H, Song S R, et al. The volcanic periodicity in central and southern Luzon, Philippines:as inferred from tephra study in the ocean core. Available online at https://gra103.aca.ntu.edu.tw/gdoc/D90224003a.pdf
[17] 时振波,林晓彤,杨慧群,等. 南海东部晚更新世以来的火山沉积特征[J]. 中国海洋大学学报,2004,34(6):1063-1068. [SHI Zhen-bo, LIN Xiao-tong, YANG Hui-qun, et al. Volcanic sediments in the eastern South China Sea during the late Pleistocene epoch[J]. Periodical of Ocean University of China, 2004, 34(6):1063-1068.]
[18] 王汝建. 南沙海区晚更新世以来的火山灰及其地质意义[J]. 海洋地质与第四纪地质,2000, 20(1):51-56. [WANG Ru-jian. Pleistocene volcanic shards and their geological significance in the Nansha area of the South China Sea[J]. Marine Geology and Quaternary Geology, 2000, 20(1):51-56.]
[19] 陈忠,夏斌,颜文,等. 南海火山玻璃的分布特征、化学成分及源区探讨[J]. 海洋学报,2005,27(5):73-81. [CHEN Zhong, XIA Bin, YAN Wen, et al. Distribution, chemical characteristics and source area of volcanic class in the South China Sea[J]. Acta Oceanologica Sinica, 2005, 27(5):73-81.]
[20] Laj C, Wang P, Balut Y. MD147-Marco Polo IMAGESⅫ Cruise Report[R]. France:Institut Paul-Emile Victor (IPEV), 2005.
[21] 谢昕,郑洪波,陈国成,等.古环境研究中深海沉积物粒度测试的预处理方法[J].沉积学报,2007(待刊).[XIE Xin, ZHENG Hong-bo, CHEN Guo-cheng, et al. Pretreatment method of grain size measurement of marine sediments in paleoenvironment research[J]. Acta Sedimentologica Sinica, 2007(in press).] [22] Imbrie J, Hays J D, Martinson D G, et al. The orbital theory of pleistocene climate:support from a revised chronology of the marine δ18O record[C]//In:Berger A, Imbrie J, Hays J D (eds). Milankovitch and Climate. Holland:D. Reidel, 1984:269-307.
[23] 汪品先. 十五万年来的南海[M]. 上海:同济大学出版社,1995.[WANG Pin-xian (eds.). The South China Sea During the Last 150000 Years[M]. Shanghai:Tongji University Press, 1995.] [24] Wang Pin-xian, Wang Lue-jiang, Bian Yun-hua, et al. Late Quqternary paleoceanography of the South China Sea:surface circulation and carbonate cycles[J]. Marine Geology, 1995, 127:145-165.
[25] Barton C E, Lackie M, Peerdeman F M. Environmental control of magnetic properties of upper-slope sediments near the great barrier reef:results from leg 133, site 820[J]. Proceedings of the Ocean Drilling Program, Scientific Results. 1993, 133:543-562.
-
期刊类型引用(5)
1. 孟庆奎，张文志，高维，舒晴，李瑞，徐光晶，张凯淞. 重力位场小波多尺度分解性质的分析与应用. 物探与化探. 2022(04): 946-954 . 
百度学术
2. 沈立祥，王智明，倪冉. 水下磁性目标物空间位置定位研究. 海洋测绘. 2022(04): 5-8 . 百度学术
3. 孟庆奎，舒晴，高维，周坚鑫，张文志，李勇. 一种用于重力场分离的约束插值切割方法. 石油地球物理勘探. 2022(05): 1204-1217+1008-1009 . 百度学术
4. 王亚章，杜晓鹏，王正虎. 基于插值切割法的水下磁性物体深度反演. 安徽地质. 2022(S2): 113-115 . 百度学术
5. 张菲菲，孟祥君，韩波，田振兴. 辽东湾地区重、磁异常特征及其区域构造分析. 海洋地质与第四纪地质. 2019(03): 104-112 . 本站查看
其他类型引用(2)
下载:
计量
- 文章访问数: 1826
- HTML全文浏览量: 229
- PDF下载量: 22
- 被引次数: 7
