长江与黄河黏土粒级沉积物地球化学特征及其物源指示意义

艾丽娜, 韩宗珠, 吴晓, 毕乃双, 王厚杰

艾丽娜, 韩宗珠, 吴晓, 毕乃双, 王厚杰. 长江与黄河黏土粒级沉积物地球化学特征及其物源指示意义[J]. 海洋地质与第四纪地质, 2020, 40(3): 109-118. DOI: 10.16562/j.cnki.0256-1492.2019050603
引用本文: 艾丽娜, 韩宗珠, 吴晓, 毕乃双, 王厚杰. 长江与黄河黏土粒级沉积物地球化学特征及其物源指示意义[J]. 海洋地质与第四纪地质, 2020, 40(3): 109-118. DOI: 10.16562/j.cnki.0256-1492.2019050603
AI Lina, HAN Zongzhu, WU Xiao, BI Naishuang, WANG Houjie. Geochemical characteristics of clay-sized sediments of the Yangtze and Yellow Rivers and their implications for provenance[J]. Marine Geology & Quaternary Geology, 2020, 40(3): 109-118. DOI: 10.16562/j.cnki.0256-1492.2019050603
Citation: AI Lina, HAN Zongzhu, WU Xiao, BI Naishuang, WANG Houjie. Geochemical characteristics of clay-sized sediments of the Yangtze and Yellow Rivers and their implications for provenance[J]. Marine Geology & Quaternary Geology, 2020, 40(3): 109-118. DOI: 10.16562/j.cnki.0256-1492.2019050603

长江与黄河黏土粒级沉积物地球化学特征及其物源指示意义

基金项目: 国家杰出青年科学基金“河口海岸学:现代黄河入海沉积物从源到汇的关键沉积动力过程”(41525021)
详细信息
    作者简介:

    艾丽娜(1989—),女,博士生,海洋地质专业,E-mail:ailinaOUC@163.com

    通讯作者:

    王厚杰(1972—),男,教授,主要从事近海沉积动力学研究,E-mail:hjwang@ouc.edu.cn

  • 中图分类号: P736.4

Geochemical characteristics of clay-sized sediments of the Yangtze and Yellow Rivers and their implications for provenance

  • 摘要: 本文对长江与黄河口黏土粒级沉积物主微量元素的地球化学特征进行了研究。结果表明:长江相对富集Al、K、Fe、Ti等常量元素以及Cr、V、Li、Zn、Ni和Rb等微量元素,黄河以高Ca、Sr和Ba为特征;长江沉积物中稀土元素含量高于黄河沉积物,长江与黄河黏土粒级沉积物中稀土元素的分馏程度相同,均具有轻稀土元素富集、重稀土元素亏损的右倾的球粒陨石标准化配分模式,上陆壳标准化配分模式为中稀土元素富集,黄河沉积物的Ce负异常和Eu正异常程度较长江沉积物偏弱。长江沉积物中元素含量变化较大,而黄河沉积物中元素含量较为稳定。黏土粒级沉积物的稀土元素更接近其物源区,<2 μm的黏土粒级沉积物中的REE可以作为判识长江与黄河沉积物的地球化学指标。ΣREE、δCe可以作为长江、黄河入海沉积物的判别指标。长江与黄河黏土粒级沉积物的地球化学特征受源区岩石类型、化学风化、水动力分选等因素的控制。化学风化引起河流沉积物相对其源岩发生地球化学分异,水动力分选使矿物在不同粒级沉积物中富集,从而引起地球化学组成在不同粒级沉积物中的分异。
    Abstract: In this paper, the element compositions and the geochemical characteristics of the clay-sized sediments of the Yangtze and Yellow Rivers estuaries were studied. It is found that the sediments in the Yangtze River is relatively rich in major elements Al, K, Fe, Ti and trace elements Cr, V, Li, Zn, Ni, Rb, whereas the Yellow River sediments characterized by Ca, Sr and Ba. As far as the rare earth elements are considered, the sediments from the Yangtze River contain more rare earth elements (REEs) with same fractionation degree comparing to those from the Yellow River. Chondrite standardized curve of rare earth elements reveals a pattern of LREE enrichment and HREE heavy loss in a slightly right-wing type. In the standardized distribution model by the upper continental crust, medium REEs are enriched. However, the negative Ce anomaly and the positive Eu anomaly of the Yellow River sediments are smaller than those in the Yangtze River sediments. The content of elements in the sediments of the Yangtze River varies greatly. Instead it is relatively concentrated in the Yellow River sediments. RERs, ΣREE, δCe of clay-sized sediments are closer to the source rocks and can be used as a geochemical indicator to distinguish the sediments of the Yangtze River from the sediments of the Yellow River. The geochemical characteristics of clay-sized sediments from the Yangtze and Yellow Rivers are controlled by lithology, chemical weathering and hydrodynamic sorting in the source area. Chemical weathering caused the geochemical differentiation of river sediments starting from the source, and the hydrodynamic environment contributed to the differentiation and separation of minerals of different sizes in the transportation.
  • 西太平洋地区因拥有世界上最独特的双俯冲系统、最完整的“沟-弧-盆”体系、最宽的大陆架、最大的边缘海、最深的海沟、最老的洋壳、最多的海山等而成为地学界最关注的区域[1]。西太平洋中部地区恰是西太平洋板块边缘沟-弧-盆体系构造演化的关键区域,其地质特征与构造演化一直是地学家关注的焦点问题之一。区内发育有海盆、海槽、海沟、不活跃和活跃的洋中脊、裂谷以及活跃的岛弧,其地球动力学背景和构造过程极为复杂。尽管通过一些地球物理资料,对于西太平洋中部的部分区域的地壳构造和年龄进行了研究[2-4],但对整体的岩石圈结构和流变性的了解相对较少。

    岩石圈有效弹性厚度(Te)是一个表征岩石圈长期强度和动力学响应的重要参数,其对应于一个理想的弹性板的厚度,在相同的荷载作用下它可以产生与岩石圈相同的弯曲程度[5]。岩石圈有效弹性厚度主要与幂律蠕变参数有关,例如温度、岩石圈物质组成、地壳厚度、应力状态等[6-9]。因此,有效弹性厚度能提供岩石圈的流变性和应力状态的重要信息[7, 10-11],并为定量研究岩石圈强度变化和标定其屈服强度提供必要数据[12]。前人对于西太平洋中部岩石圈的有效弹性厚度进行了一定程度的研究[13-17],但是大都没有考虑内部荷载的影响,这样会导致Te值被低估[5]。表面荷载作用主要由地形变化和地表密度的大尺度变化(如山地和沉积盆地)引起,内部荷载主要包括侵入的铁镁质物质、增生的下地壳物质、热异常和地壳物质成分的变化,这些会造成岩石圈内部密度的横向变化。本文在挠曲模型中将表面荷载和内部荷载同时考虑在内,得到的有效弹性厚度将更加准确地反映西太平洋中部岩石圈的强度变化。

    目前有效弹性厚度的研究方法很多[18-20],其中通过分析自由空间重力异常和地形之间的导纳来计算岩石圈Te值的滑动窗口导纳法,是当前最为常见的一种方法。它被广泛应用于计算不同地质背景下的海洋岩石圈和大陆岩石圈的Te值,如欧洲[21]、太平洋[14, 16-17]以及俯冲带[22]等。鉴于此,本文拟利用西太平洋中部的水深数据、自由空间重力异常数据和地壳结构数据,并在挠曲模型中考虑表面荷载和内部荷载同时存在的情况,采用滑动窗口导纳技术计算其岩石圈的有效弹性厚度值,以进一步探讨该区的岩石圈结构和地质构造过程。

    西太平洋中部主要由太平洋板块、菲律宾板块和卡罗琳板块组成,它们之间以俯冲带为边界(图1)。太平洋板块是中生代海底的重要组成部分,存在大量的海底火山。太平洋板块向西北运动,俯冲到菲律宾板块底部,其以25°N为界由北向南大致可以分为两部分:呈线型的伊豆–小笠原海沟,呈弓型的马里亚纳海沟。

    图  1  西太平洋中部的水深地形图
    Figure  1.  The bathymetry map of the Central Western Pacific

    菲律宾板块是一个相对较小的板块,被东部的太平洋板块和西部的欧亚板块夹在中间,南部以印度-澳大利亚板块和卡罗琳板块为边界。沿着东部边缘,菲律宾板块与太平洋板块形成了伊豆-小笠原-马里亚纳海沟。在西部边缘,菲律宾沿着南海海槽、琉球海沟和菲律宾海沟俯冲到欧亚板块底部,而在中部的马尼拉海沟却为俯冲带的上覆板块[23]。菲律宾板块由海底扩张形成的西菲律宾海盆、四国海盆和帕里西维拉海盆组成。位于菲律宾板块中部呈南北走向的非震海岭九州-帕劳海脊(KPR)将西菲律宾海盆与四国海盆和帕里西维拉海盆分隔开。由于四国海盆和帕里西维拉海盆的打开(约25 Ma),KPR从IBM火山弧分离出来[24]。在3~4 Ma,第二次海底扩张作用形成了马里亚纳海槽。位于西菲律宾海盆中部的中央盆地断陷被认为是一个古洋中脊,通过海底扩张作用形成了西菲律宾海盆[25-26]

    卡罗琳板块位于太平洋、菲律宾和印度-澳大利亚板块之间。卡罗琳板块曾被认为是太平洋板块的一部分,但是基于对其海底地形、重力、地震数据的分析,发现卡罗琳板块可能是一个独立的板块[27]。卡罗琳板块东部以穆索海沟为界,北部以雅浦海沟为界,南部以新几内亚海沟和马努斯海沟为界。卡罗琳板块被位于中部的欧里皮克海隆分为西卡罗琳次海盆和东卡罗琳次海盆。Gaina和Müller[28]通过对海底磁异常的研究,发现东卡罗琳次海盆的地壳以中部的古洋中脊对称分布,而西卡罗琳次海盆的古洋中脊靠近它的北部边缘,呈现明显的非对称性,认为可能是受到北部马努斯热点的影响使得洋中脊向北跳跃而导致的结果。卡罗琳海岭位于雅浦海沟的东部以及马里亚纳海沟的南部,被认为是由地幔柱的火山作用形成[4],索罗尔海槽将其分为东卡罗琳海岭和西卡罗琳海岭。在早中新世,卡罗琳海岭与雅浦火山弧相撞,导致了雅浦火山弧火山作用停止[29],帕里西维拉海盆逆冲到雅浦火山弧上面[30],马里亚纳火山弧弯曲[31]

    地形和测深数据来源于网格分辨率为30"×30"的GEBCO_2014数字高程模型,其数据主要源于船测数据和卫星测高反演数据(图1)。刘洋[32]在马里亚纳海沟“挑战者深渊”,以多波束条带数据的中央波束测深结果为基准,对GEBCO_2014数据进行精度评估,结果表明其精度在限差容许范围内(平均差值为−10.7 m)。自由空气重力异常数据来源于网格分辨率为2.5ˊ×2.5ˊ的地球重力模型EGM2008[33]图2a)。EGM2008重力模型的球谐展开阶数达到了2159,其采用了卫星跟踪数据、卫星测高数据、地面重力数据和海洋重力数据等。张精明等[34]在中国大陆对EGM2008模型精度进行了分析,结果显示各区域内EGM2008模型计算高程与中国高程系统水准高程对比误差精度均在1 m 以内,表明其具有较高的精度。此外,计算岩石圈Te值的过程需要了解地壳和上地幔的内部结构,为了确定内部不同界面的密度分布和横向变化,我们使用了全球地壳模型CRUST1.0[35],CRUST1.0包含冰层、水层、3层沉积物和3层地壳,第8层为Moho面深度(图2b)。根据Forsyth[5]的假设,我们认为初始内部荷载发生在Moho面。尽管Te和假定的荷载深度存在一种权衡关系,但是Pére-Gussinyé和Watts[21]通过测试Te对荷载深度的敏感度,发现将荷载深度从中地壳更改到Moho面,Te只改变了约5 km且变化模式保持不变,表明荷载深度的差别对于Te计算结果的影响较小。

    图  2  西太平洋中部的自由空气重力异常图(a)和Moho面深度图(b)
    Figure  2.  Free-air gravity anomaly(a)and depth of the Moho(b)of the Central Western Pacific

    为了消除曲线坐标的平面处理而产生的误差,我们使用墨卡托投影将所有数据投影到笛卡尔坐标系中,并将所有数据重采样为网格分辨率为8 km的数据。此外,为了消除边缘效应,我们选取了比研究区域更大的网格数据用于计算Te结果。

    岩石圈有效弹性厚度的计算需要根据模型,最简单的模型是一个位于无黏性流体上的均匀弹性板模型,它在荷载作用下发生挠曲[5]。对于一个给定的Te值,我们可以从观测重力异常和地形中得到初始的表面荷载和内部荷载,从而根据挠曲模型重构重力异常和地形,称为荷载反褶积法[5]

    观测导纳函数Zobsk)为:

    $$ {Z}_{\rm{obs}}\left(k\right)=<\frac{{S}_{\rm{fh}}\left(k\right)}{{S}_{\rm{hh}}\left(k\right)}> $$ (1)

    式中,k为二维波数,Sfhk)和Shhk)分别表示地形和自由空间重力异常的互功率谱密度和地形的自功率谱密度,<>表示一个环形波数带上的平均值。

    岩石圈在荷载作用下会产生挠曲,其弯曲程度受岩石圈挠曲强度D的控制,公式如下:

    $$ D=\frac{E{\rm{T}}_{\rm{e}}^{3}}{12(1-{v}^{2})} $$ (2)

    式中,E为杨氏模量,v为泊松比,Te为岩石圈有效弹性厚度(具体数值见表1)。在假定Ev保持不变的情况下,Te的大小决定了岩石圈的挠曲强度,较大的Te对应于一个较强的岩石圈,表明岩石圈的很大一部分荷载由弹性应力支撑,岩石圈弯曲程度较小;而较小的Te对应于一个较弱的岩石圈,表明弹性应力支撑了较小部分的荷载,岩石圈弯曲程度较大。

    表  1  常量的符号和数值
    Table  1.  Symbols and values of constants
    常量符号数值单位
    杨氏模量E100GPa
    牛顿引力常数G6.6726m3·kg−1· s−2
    泊松比v0.25
    重力加速度g9.79m·s−2
    海水密度ρw1030kg·m−3
    地壳密度ρc2800kg·m−3
    地幔密度ρm3300kg·m−3
    下载: 导出CSV 
    | 显示表格

    弹性板块在荷载作用下的挠曲方程[36]

    $$ D{\nabla }^{4}w+\left({\rho }_{\rm{m}}-{\rho }_{\rm{fill}}\right)gw=q\left(x\right) $$ (3)

    式中,D为岩石圈的挠曲强度,w为弹性挠曲,ρmρfill分别为地幔密度和填充挠曲的物质密度,g为重力加速度,qx)为弹性板块受到的垂向荷载(具体数值见表1)。

    假定岩石圈厚度在挠曲过程中保持不变,对于一个作用在岩石圈表面且密度为$ \rho $c的初始表面荷载(HI),其导致的最终地形和自由空气重力异常分别为HTGT[22]

    $$ {H}_{\rm{T}}={H}_{\rm{I}}\frac{D{\left|k\right|}^{4}+({\rho }_{\rm{m}}-{\rho }_{\rm{c}})g}{D{\left|k\right|}^{4}+({\rho }_{\rm{m}}-{\rho }_{\rm{w}})g} $$ (4)
    $$ \begin{split} {G_{\rm{T}}} = &{H_I}\frac{{2\pi G({\rho _{\rm{c}}} - {\rho _{\rm{w}}})}}{{\left( {{\rho _{\rm{m}}} - {\rho _{\rm{w}}}} \right) + D{{\left| k \right|}^4}/g}}\left[ {\frac{{D{{\left| k \right|}^4}}}{g} + } \right.\\ & {({\rho _{\rm{m}}} - {\rho _{\rm{c}}})(1 - {\rm{exp}}\left( { - kt} \right))]{\rm{exp}}\left( { - kd} \right)} \Biggr] \end{split}$$ (5)

    式中,G为牛顿引力常数,ρcρw分别为地壳密度和海水密度,td分别为平均地壳厚度和平均水深。

    对于一个作用在Moho面且密度差为(ρmρc)的初始内部荷载(WI),其导致的最终地形和自由空气重力异常分别为HBGB[22];

    $$ {H}_{\rm{B}}={W}_{\rm{I}}\frac{-({\rho }_{\rm{m}}-{\rho }_{\rm{c}})g}{D{\left|k\right|}^{4}+({\rho }_{\rm{m}}-{\rho }_{\rm{w}})g} $$ (6)
    $$\begin{split} &{G_{\rm{B}}} = {W_{\rm{I}}}\frac{{2\pi G({\rho _{\rm{m}}} - {\rho _{\rm{c}}})}}{{\left( {{\rho _{\rm{m}}} - {\rho _{\rm{w}}}} \right) + D{{\left| k \right|}^4}/g}}\\ &\left[ {{\rm{exp}}\left( { - kt} \right)\left( {\frac{{D{{\left| k \right|}^4}}}{g} + {\rho _{\rm{c}}} - {\rho _{\rm{w}}}} \right) - ({\rho _{\rm{c}}} - {\rho _{\rm{w}}})} \right]{\rm{exp}}\left( { - kd} \right) \end{split} $$ (7)

    当表面荷载和底部荷载同时作用于海洋岩石圈时,最终的地形(H)和自由空气重力异常(G)是两者各自单独作用于岩石圈所产生的地形和自由空气重力异常的累加,公式如下:

    $$ H={H}_{\rm{T}}+{H}_{\rm{B}} $$ (8)
    $$ G={G}_{\rm{T}}+{G}_{\rm{B}} $$ (9)

    理论导纳函数Zprek)为:

    $$ {Z}_{\rm{pre}}\left(k\right)=<\frac{{S}_{{f}_{t}{h}_{t}}\left(k\right)+{S}_{{f}_{b}{h}_{b}}\left(k\right)}{{S}_{{h}_{t}{h}_{t}}\left(k\right)+{S}_{{h}_{b}{h}_{b}}\left(k\right)}> $$ (10)

    式中,Sfthtk)和Sfbhbk)分别表示表面荷载和底部荷载作用下的地形和自由空间重力异常的互功率谱密度,Shthtk)和Shbhbk)分别表示表面荷载和底部荷载作用下的地形的自功率谱密度。

    对于一个给定的Te值,我们首先利用观测的地形和自由空间重力异常,通过公式(8)和(9)求得表面荷载(HI)和底部荷载(WI);然后利用HIWI并通过公式(4)—(7)求得最终地形和自由空间重力异常的分量(HTGTHBGB);最后通过公式(10)计算地形和自由空间重力异常之间的理论导纳曲线。我们将观测导纳曲线与不同Te值的一系列理论导纳曲线进行对比,将匹配程度最高的理论导纳曲线所对应的Te值作为研究区的Te值。

    为了获得研究区有效弹性厚度的空间变化值,我们将其划分为多个重叠的窗口,计算每个窗口的观测导纳曲线和理论导纳曲线,将获得的岩石圈有效弹性厚度作为窗口中心处的Te值,然后将窗口移动40 km进行新的Te估计。在滑动窗口导纳法中,窗口大小的选择至关重要,它控制着Te结果的分辨率和方差之间的权衡[20, 37]。较小的窗口能更好地反映短波长的特征,提高数据分辨率,但会丢失长波长的信息;较大的窗口能突出长波长特征,提高Te的反演精度,但会丢失短波长信息,降低数据的分辨率[20, 37]。为了获得研究区较为精确的Te估计值和较好的Te空间变化性,我们选取了比较广泛的窗口范围,分别为200 km×200 km、400 km×400 km、600 km×600 km和800 km×800 km。

    图3展示了4个不同窗口计算的西太平洋中部的Te空间变化结果。虽然4个不同窗口的Te结果在绝对值和短波长变化方面存在差异,但是总体的Te变化模式看起来非常相似。计算结果显示,研究区的Te值整体上处于0~50 km,这与前人的计算结果相一致[14-15, 17],Te变化基本上与构造单元相吻合,且与主要的构造边界密切相关。通常情况下,Te的变化相对平缓,尤其在海洋岩石圈,然而,我们的结果显示出部分区域的Te值在很小的空间间隔内发生了很大的变化,从而产生了很大的Te梯度,尤其是200 km×200 km窗口和400 km×400 km窗口的Te结果,这种现象在使用频谱法反演有效弹性厚度的研究中比较常见,属于反演的人为误差[38-40],在滑动窗口导纳法中主要取决于窗口大小的选择。Pérez-Gussinyé等[39]通过对不同窗口大小的反演结果进行测试,认为600 km×600 km的窗口对于0~50 km的岩石圈有效弹性厚度具有较好的反演效果。然而,由于600 km×600 km的窗口相对较大,对于一些较小的构造单元的反演效果较差,如马里亚纳火山弧。结果显示,200 km×200 km窗口的计算结果能够较好地反映研究区的构造单元之间的Te值变化。因此,本文主要利用600 km×600 km窗口(精确的Te估计值)和200 km×200 km窗口(较好的Te空间变化性)计算的Te结果,对研究区的岩石圈强度和地质构造过程进行分析。

    图  3  西太平洋中部岩石圈有效弹性厚度(Te)的空间变化图
    a-d分别代表200 km×200 km窗口、400 km×400 km窗口、600 km×600 km窗口和800 km×800 km窗口的计算结果。
    Figure  3.  Te estimates in the Central Western Pacific
    Panels(a)to(d)show Te results for four different window sizes(200 km×200 km, 400 km×400 km, 600 km×600 km, and 800 km×800 km)

    除海底火山具有相对较小的Te值(15~20 km)外,太平洋板块整体上具有较强的岩石圈强度(25~30 km)。马里亚纳海沟东侧的岩石圈具有西太平洋中部最大的Te值(30~50 km),这与Zhang等[15]计算的结果相对一致。Zhang等[15]对马里亚纳海沟附近的岩石圈进行弹性挠曲研究,结果显示外隆起背离海沟方向的岩石圈有效弹性厚度为45~53 km。马里亚纳海沟东侧的Te值具有较大的空间变化性,由北向南整体上表现为‘高-低-高-低’的Te变化特征(图3a, b)。马里亚纳海沟西侧的马里亚纳火山弧在大窗口(图3b-d)的计算结果中表现为高Te值特征(30~50 km),这可能是由于窗口较大导致的平滑作用,使得马里亚纳火山弧的岩石圈短波长信息被马里亚纳海沟东侧岩石圈的长波长信息所覆盖,没有凸显出来。与之相对,在200 km×200 km的窗口中(图3a),马里亚纳海沟和马里亚纳火山弧均表现为明显的较低的Te值特征(约 12 km)。马里亚纳海槽具有异常弱的岩石圈强度,Te值为3~5 km。帕里西维拉海盆中部具有较强的岩石圈强度,Te值约为30 km,其西部岩石圈强度明显弱于东部岩石圈(东部的Te值比西部的Te值高约10 km)。西菲律宾海盆具有与帕里西维拉海盆相似的Te变化特征。西菲律宾海盆的中央盆地断陷的岩石圈表现为高Te值特征(约18 km),其南北两侧次海盆的岩石圈表现为低Te值特征(<10 km)。靠近菲律宾俯冲带的区域,岩石圈的Te值明显增大(约25 km)。卡罗琳板块的岩石圈整体上表现为相对均一的低Te值特征(<15 km),尤其是欧里皮克海隆、卡罗琳海岭和索罗尔海槽,它们的岩石圈Te值为3 km,这与杨安等[17]的计算结果相一致。卡罗琳海岭东侧的卡罗琳火山链在200 km×200 km的窗口中表现为明显的低Te值特征(约5 km)(图3a)。卡罗琳板块北部的雅浦海沟和东部的穆索海沟具有较高的Te值(25~35 km),表明它们具有相对较强的岩石圈强度。位于卡罗琳板块西部的阿尤海槽和南部的马努斯海沟具有较弱的岩石圈强度(约3 km)。

    大家普遍认为岩石圈强度与岩石圈温度和年龄之间存在一定的相关性[7, 41-43],岩石圈会随着年龄的增加而不断变冷、变强。然而,后期的构造-岩浆活动会导致岩石圈的热再生和/或岩石圈组分的改变,使得岩石圈强度降低[39, 41],如热点或火山作用引起的热扰动以及伸展或断层作用引起的岩石圈减薄或破裂。因此,Te的空间变化可以反映出构造单元的地质构造过程。

    太平洋板块虽然是一个较老的板块,但是我们的结果显示,它并没有表现为相对均一的高岩石圈强度,表明其岩石圈流变强度是不均匀的。前人利用滑动窗口导纳法对于太平洋的岩石圈有效弹性厚度进行了计算[14, 16-17],结果显示Te值主要为8~15 km,明显小于我们的计算结果。在计算Te值时,他们使用的模型只考虑了表面荷载,没有考虑内部荷载的情况,这可能导致低估了Te值的大小[5]。太平洋板块后期受到火山作用的影响形成了众多的海底火山,强烈的岩浆活动导致岩石圈的热再生和/或岩石圈组分的改变,使得岩石圈的强度显著降低。我们的结果显示,马里亚纳海沟和菲律宾海沟的岩石圈强度从外隆起到海沟方向表现为明显的减弱,表明岩石圈由外隆起向海沟发生了弱化,这可能是由于岩石圈发生挠曲形变时,产生了大量断层,断层作用使得岩石圈破裂,降低了岩石圈的强度。海洋岩石圈的弹性厚度受到岩石圈温度的控制,Te值与岩石圈温度呈负相关关系[41-43]:岩石圈的温度越高,其强度会越低。马里亚纳火山弧由许多火山组成,岩石圈处于较热状态,因此它具有较弱的岩石圈强度(约12 km)。马里亚纳海槽是最年轻的海洋盆地,形成于6~8 Ma[44-45],岩石圈应该处于均衡状态,具有异常弱的岩石圈强度(Te<5 km)。

    帕里西维拉海盆属于KPR和西马里亚纳残留弧之间的弧间盆地,在西菲律宾海盆停止扩张后(约30 Ma),裂谷作用使得四国海盆和帕里西维拉海盆开始形成,然后在29~26 Ma和15 Ma之间发生海底扩张作用[46]。我们的结果显示,帕里西维拉海盆中部的扩张中心表现为高Te值特征,表明岩石圈处于非均衡状态。张正一等[3]通过对重力数据、磁力数据和地震数据的分析,发现帕里西维拉海盆西部具有比东部更高的磁异常和更厚的磁性体。因此,我们推测,帕里西维拉海盆西部相较于东部具有较弱的岩石圈强度,这可能与海盆的非对称扩张有关,扩张过程中岩浆可能更加集中于海盆的西侧,导致岩浆作用更加强烈,岩石圈更热,从而使其具有更低的岩石圈Te值。西菲律宾海盆被中央盆地断陷分为南北两个次海盆,除南部次海盆靠近菲律宾海沟的岩石圈具有较大的Te值外,两个次海盆整体上表现为相同的低Te值特征,表明西菲律宾海盆在海底扩张过程中岩石圈处于相对均衡状态。

    我们的结果显示,欧里皮克海隆具有异常低的Te值,表明它的岩石圈处于均衡状态,这与Altis[13]的观点相一致,Altis通过对测深数据和地震数据的分析,认为欧里皮克海隆是由强烈的火山作用所形成,岩石圈处于Airy均衡状态,因此其具有较弱的岩石圈强度。Zhang 等[4]认为,卡罗琳海岭代表了一个中等大小的海底高原,其火山活动来源于地幔柱(卡罗琳热点)。强烈的火山活动导致岩石圈的热再生,使得卡罗琳海岭和索罗尔海槽的岩石圈处于较热状态,岩石圈强度显著降低(3 km)。阿尤海槽被认为是一个慢速扩张洋中脊,且自25 Ma以来一直处于活动状态[27]。一般情况下,洋中脊地形和重力异常之间具有较好的相关性,表明岩石圈处于Airy均衡,岩石圈应该较弱[47]。此外,我们的结果显示,位于卡罗琳板块南部的马努斯海沟具有较弱的岩石圈强度(3 km),杨安等[17]认为,此处可能是由于处于非弹性均衡导致的结果。

    (1)除海底火山具有相对较小的Te值(15~20 km)外,太平洋板块整体上具有较强的岩石圈强度(25~30 km)。马里亚纳海沟和菲律宾海沟的岩石圈强度从外隆起到海沟方向表现为明显的减弱,可能是由于岩石圈发生挠曲形变时,产生了大量断层,降低了岩石圈的强度。

    (2)帕里西维拉海盆西部相较于东部具有较弱的岩石圈强度,这可能与海盆的非对称扩张有关,扩张过程中岩浆可能更加集中于海盆的西侧,导致岩浆作用更加强烈,岩石圈更热。

    (3)卡罗琳板块的岩石圈整体上表现为相对均一的低Te值特征(<15 km)。欧里皮克海隆、卡罗琳海岭和索罗尔海槽的Te值为3 km,强烈的火山活动导致岩石圈的热再生,使得它们的岩石圈处于较热状态,岩石圈强度显著降低。

  • 图  1   取样站位

    a.黄河口,b.长江口。

    Figure  1.   Sampling positions

    a. Yellow River estuary, b. Yangtze River estuary.

    图  2   长江与黄河黏土粒级沉积物球粒陨石标准化曲线(a)与上陆壳标准化曲线(b)

    Figure  2.   Chondrite standardized curves (a) and upper continental crust standardized curves (b) for rare elements of clay-sized sediments in the Yangtze and Yellow Rivers

    图  3   长江与黄河不同粒级沉积物稀土元素配分曲线

    a. 球粒陨石标准化曲线,b. 上陆壳标准化曲线。

    Figure  3.   REE distribution curves of different grain-sized sediments of Yangtze and Yellow Rivers

    a. Chondrite standardized curves, b. Upper Continental Crust standardized curves.

    图  4   REE-δCe判别图

    Figure  4.   Discrimination diagram of REE-δCe

    图  5   长江干流沉积物Sr-Nd同位素、(La/Yb)UCC、(Gd/Yb)UCC变化、上陆壳标准化曲线 (长江沉积物同位素数据来自[30-31],黄河沉积物同位素数据来自[32],长江沉积物REE数据来自[33])

    Figure  5.   Sr-Nd isotope,(La/Yb)UCC,(Gd/Yb)UCC, and upper continental crust standardized curves of sediments from the main stream of Yangtze (Sr-Nd isotope data of the Yangtze and Yellow Rivers sediments from [30-31] and [32], respectively, REE of Yangtze Sediment from [33])

    表  1   长江与黄河黏土粒级沉积物的常量元素含量(%)

    Table  1   Major element contents in the clay-sized sediments of the Yangtze and Yellow Rivers(%)

    Al2O3CaOTFe2O3K2OMgOMnONa2OP2O5TiO2
    黄河平均值19.6110.038.593.383.820.170.410.290.63
    标准偏差0.520.650.240.080.090.010.020.020.02
    变异系数2.646.462.762.362.377.405.106.833.43
    长江平均值23.501.999.813.623.240.190.370.300.90
    标准偏差0.500.600.470.130.110.030.040.020.07
    变异系数2.1130.014.823.613.2817.569.547.297.77
    下载: 导出CSV

    表  2   长江与黄河黏土粒级沉积物的微量元素含量(μg/g)

    Table  2   Trace element contents in the clay-sized sediments of the Yangtze and Yellow Rivers(μg/g)

    SrCrVZnBaLiBeScCo
    黄河平均值226.81106.33140.96149.13578.0076.913.4824.7223.48
    标准偏差8.952.425.176.5722.833.990.140.710.68
    变异系数3.942.273.674.403.955.193.892.862.91
    长江平均值105.72130.12178.84193.58542.67105.774.1826.3427.88
    标准偏差11.098.429.2939.3443.523.200.170.852.19
    变异系数10.496.475.1920.328.023.034.163.247.86
    NiCuGaRbMoCsPbThU
    黄河平均值59.1960.5527.45169.601.3215.9453.8220.512.80
    标准偏差1.634.040.754.850.120.452.810.550.11
    变异系数2.766.682.732.869.102.835.212.684.05
    长江平均值70.3565.5732.74218.730.9520.5260.4722.223.58
    标准偏差4.1413.790.952.610.290.8613.301.580.33
    变异系数5.8921.032.901.1930.314.1822.007.129.16
    下载: 导出CSV

    表  3   长江与黄河黏土粒级沉积物的稀土元素含量及参数

    Table  3   REE contents and parameters in the clay-sized sediments of the Yangtze and Yellow Rivers

    ΣREE/(μg/g)LREE/HREEδEuδCe(La/Yb)UCC(La/Sm)UCC(Gd/Yb)UCC
    黄河平均值226.448.760.660.970.900.981.13
    标准偏差6.690.100.010.010.010.010.02
    变异系数2.951.090.880.950.571.391.62
    长江平均值243.578.750.680.930.910.991.12
    标准偏差25.780.100.010.020.030.020.03
    变异系数10.591.181.081.753.171.743.14
    下载: 导出CSV

    表  4   长江与黄河不同粒级沉积物稀土元素组成及参数

    Table  4   REE contents and parameters in different sized sediments from the Yangtze and Yellow Rivers

    ΣREE/(μg/g)LREE/HREEδEuδCe(La/Yb)UCC(La/Sm)UCC(Gd/Yb)UCC数据来源
    长江(全样)170.958.430.660.841.170.891.44[20]
    长江(<63 μm)167.108.480.670.851.190.891.45[21]
    长江(<2 μm)243.578.750.680.930.990.911.12本文
    黄河(全样)139.118.150.650.861.040.871.32[20]
    黄河(<63 μm)137.768.120.650.861.040.871.31[21]
    黄河(<2 μm)226.448.760.660.970.980.901.13本文
    下载: 导出CSV
  • [1] 胡邦琦, 李军, 李国刚, 等. 长江和黄河入海沉积物的物源识别研究进展[J]. 海洋地质与第四纪地质, 2011, 31(6):147-156. [HU Bangqi, LI Jun, LI Guogang, et al. Distinguishing the Changjiang and Huanghe sediments: a review [J]. Marine Geology & Quaternary Geology, 2011, 31(6): 147-156.
    [2]

    Yang S Y, Jung H S, Lim D I, et al. A review on the provenance discrimination of sediments in the Yellow Sea [J]. Earth-Science Reviews, 2003, 63(1-2): 93-120. doi: 10.1016/S0012-8252(03)00033-3

    [3] 杨守业. 亚洲主要河流的沉积地球化学示踪研究进展[J]. 地球科学进展, 2006, 21(6):648-655. [YANG Shouye. Advances in sedimentary geochemistry and tracing applications of Asian rivers [J]. Advances in Earth Science, 2006, 21(6): 648-655. doi: 10.3321/j.issn:1001-8166.2006.06.013
    [4] 杨守业, 韦刚健, 石学法. 地球化学方法示踪东亚大陆边缘源汇沉积过程与环境演变[J]. 矿物岩石地球化学通报, 2015, 34(5):902-910. [YANG Shouye, WEI Gangjian, SHI Xuefa. Geochemical approaches of tracing source-to-sink sediment processes and environmental changes at the East Asian continental margin [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2015, 34(5): 902-910. doi: 10.3969/j.issn.1007-2802.2015.05.003
    [5] 赵一阳, 鄢明才. 黄河、长江、中国浅海沉积物化学元素丰度比较[J]. 科学通报, 1992, 37(23):1991-1994. [ZHAO Yiyang, YAN Mingcai. Abundance of chemical elements in sediments from the Huanghe river, the Changjiang river and the continental shelf of China [J]. Chinese Science Bulletin, 1992, 37(23): 1991-1994.
    [6] 赵一阳, 鄢明才. 中国浅海沉积物地球化学[M]. 北京: 科学出版社, 1994.

    ZHAO Yiyang, YAN Mingcai. Geochemistry of sediments of the China shelf sea[M]. Beijing: Science Press, 1994.

    [7] 蓝先洪, 王红霞, 李日辉, 等. 南黄海沉积物常量元素组成及物源分析[J]. 地学前缘, 2007, 14(4):197-203. [LAN Xianhong, WANG Hongxia, LI Rihui, et al. Major elements composition and provenance analysis in the sediments of the South Yellow Sea [J]. Earth Science Frontiers, 2007, 14(4): 197-203. doi: 10.3321/j.issn:1005-2321.2007.04.021
    [8] 蓝先洪, 梅西, 李日辉, 等. 晚更新世以来南黄海北部泥质区沉积物元素的分布及来源分析[J]. 现代地质, 2015, 29(4):777-788. [LAN Xianhong, MEI Xi, LI Rihui, et al. Distribution and source analysis of elements from sediments in the Northern South Yellow Sea since the late Pleistocene [J]. Geoscience, 2015, 29(4): 777-788. doi: 10.3969/j.issn.1000-8527.2015.04.007
    [9] 赵一阳, 王金土, 秦朝阳, 等. 中国大陆架海底沉积物中的稀土元素[J]. 沉积学报, 1990, 8(1):37-43. [ZHAN Yiyang, WANG Jintu, QIN Chaoyang, et al. Rare-earth elements in continental shelf sediments of the China seas [J]. Acta Sedimentologica Sinica, 1990, 8(1): 37-43.
    [10] 鄢明才, 迟清华, 顾铁新, 等. 中国东部上地壳化学组成[J]. 中国科学(D辑), 1997, 40(5):530-539. [YAN Mingcai, CHI Qinghua, GU Tiexin, et al. Chemical composition of upper crust in eastern China [J]. Science in China Series D: Earth Sciences, 1997, 40(5): 530-539. doi: 10.1007/BF02877620
    [11] 杨守业, 李从先. 长江与黄河沉积物元素组成及地质背景[J]. 海洋地质与第四纪地质, 1999, 19(2):19-26. [YANG Shouye, LI Congxian. Characteristic element compositions of the Yangtze and the Yellow River sediments and their geological background [J]. Marine Geology & Quaternary Geology, 1999, 19(2): 19-26.
    [12] 杨守业, 李从先. REE示踪沉积物物源研究进展[J]. 地球科学进展, 1999, 14(2):164-167. [YANG Shouye, LI Congxian. Research progress in REE tracer for sediment source [J]. Advance in Earth Sciences, 1999, 14(2): 164-167. doi: 10.3321/j.issn:1001-8166.1999.02.010
    [13]

    McLennan S M. Weathering and global denudation [J]. The Journal of Geology, 1993, 101(2): 295-303. doi: 10.1086/648222

    [14]

    Bouchez J, Gaillardet J, France C, et al. Grain size control of river suspended sediment geochemistry: Clues from Amazon River depth profiles [J]. Geochemistry, Geophysics, Geosystems, 2011, 12(3): Q03008.

    [15]

    Garzanti E, Andó S, France-Lanord C, et al. Mineralogical and chemical variability of fluvial sediments 2. Suspended-load silt (Ganga–Brahmaputra, Bangladesh) [J]. Earth and Planetary Science Letters, 2011, 302(1-2): 107-120. doi: 10.1016/j.jpgl.2010.11.043

    [16]

    Resentini A, Malusà M G, Garzanti E. MinSORTING: An Excel® worksheet for modelling mineral grain-size distribution in sediments, with application to detrital geochronology and provenance studies [J]. Computers & Geosciences, 2013, 59: 90-97.

    [17]

    Garzanti E, Resentini A. Provenance control on chemical indices of weathering (Taiwan river sands) [J]. Sedimentary Geology, 2015, 336: 81-95.

    [18]

    Garzanti E, Andò S, Padoan M, et al. The modern Nile sediment system: Processes and products [J]. Quaternary Science Reviews, 2015, 130: 9-56. doi: 10.1016/j.quascirev.2015.07.011

    [19] 韩宗珠, 艾丽娜, 陈筱林, 等. 南黄海泥质区西北缘B01孔黏土粒级沉积物地球化学特征及其物质来源的识别[J]. 中国海洋大学学报: 自然科学版, 2016, 46(10):82-91. [HAN Zongzhu, AI Li’na, CHEN Xiaolin, et al. Geochemical characteristics of sediments and provenance of B01 core in northwest margin of South Yellow Sea mud area [J]. Periodical of Ocean University of China, 2016, 46(10): 82-91.
    [20] 杨守业, 李从先. 长江与黄河沉积物REE地球化学及示踪作用[J]. 地球化学, 1999, 28(4):374-380. [YANG Shouye, LI Congxian. REE geochemistry and tracing application in the Yangtze River and the Yellow River sediments [J]. Geochimica, 1999, 28(4): 374-380. doi: 10.3321/j.issn:0379-1726.1999.04.008
    [21] 杨守业, 李从先. 元素地球化学特征的多元统计方法研究——长江与黄河沉积物元素地球化学研究[J]. 矿物岩石, 1999, 19(1):63-67. [YANG Shouye, LI Congxian. Multiple statistic study of element geochemical characteristics element geochemical study on the Changjiang and Huanghe sediments [J]. Journal of Mineralogy and Petrology, 1999, 19(1): 63-67.
    [22]

    Jung H S, Lim D, Jeong D H, et al. Discrimination of sediment provenance in the Yellow Sea: Secondary grain-size effect and REE proxy [J]. Journal of Asian Earth Sciences, 2016, 123: 78-84. doi: 10.1016/j.jseaes.2016.03.020

    [23]

    Jung H S, Lim D, Choi J Y, et al. Rare earth element compositions of core sediments from the shelf of the South Sea, Korea: Their controls and origins [J]. Continental Shelf Research, 2012, 48: 75-86. doi: 10.1016/j.csr.2012.08.008

    [24]

    Lim D, Choi J Y, Shin H H, et al. Multielement geochemistry of offshore sediments in the southeastern Yellow Sea and implications for sediment origin and dispersal [J]. Quaternary International, 2013, 298: 196-206. doi: 10.1016/j.quaint.2013.01.004

    [25]

    Lim D, Jung H S, Choi J Y. REE partitioning in riverine sediments around the Yellow Sea and its importance in shelf sediment provenance [J]. Marine Geology, 2014, 357: 12-24. doi: 10.1016/j.margeo.2014.07.002

    [26] 王金土. 黄海表层沉积物稀土元素地球化学[J]. 地球化学, 1990(1):44-53. [WANG Jintu. Ree geochemistry of surficial sediments from the Yellow sea of China [J]. Geochimica, 1990(1): 44-53. doi: 10.3321/j.issn:0379-1726.1990.01.005
    [27]

    Yang S Y, Jung H S, Choi M S, et al. The rare earth element compositions of the Changjiang (Yangtze) and Huanghe (Yellow) river sediments [J]. Earth and Planetary Science Letters, 2002, 201(2): 407-419. doi: 10.1016/S0012-821X(02)00715-X

    [28] 茅昌平, 陈骏, 袁旭音, 等. 长江下游悬浮物Sr-Nd同位素组成的季节性变化与物源示踪[J]. 科学通报, 2011, 56(22):2371-2378. [MAO Changping, CHEN Jun, YUAN Xuyin, et al. Seasonal variations in the Sr-Nd isotopic compositions of suspended particulate matter in the lower Changjiang River: Provenance and erosion constraints [J]. Chinese Science Bulletin, 2011, 56(22): 2371-2378. doi: 10.1007/s11434-011-4589-6
    [29] 何梦颖, 郑洪波, 黄湘通, 等. 长江流域沉积物黏土矿物组合特征及物源指示意义[J]. 沉积学报, 2011, 29(3):544-551. [HE Mengying, ZHENG Hongbo, HUANG Xiangtong, et al. Clay mineral assemblages in the Yangtze drainage and provenance implications [J]. Acta Sedimentologica Sinica, 2011, 29(3): 544-551.
    [30] 杨守业, 蒋少涌, 凌洪飞, 等. 长江河流沉积物Sr-Nd同位素组成与物源示踪[J]. 中国科学 D辑: 地球科学, 2007, 50(10):1556-1565. [YANG Shouye, JIANG Shaoyong, LING Hongfei, et al. Sr-Nd isotopic compositions of the Changjiang sediments: Implications for tracing sediment sources [J]. Science in China Series D: Earth Sciences, 2007, 50(10): 1556-1565. doi: 10.1007/s11430-007-0052-6
    [31]

    Meng X W, Liu Y G, Shi X F, et al. Nd and Sr isotopic compositions of sediments from the Yellow and Yangtze Rivers: Implications for partitioning tectonic terranes and crust weathering of the Central and Southeast China [J]. Frontiers of Earth Science in China, 2008, 2(4): 418-426. doi: 10.1007/s11707-008-0054-5

    [32] 郭玉龙. 地球化学方法示踪中国东部入海河流沉积物物质组成的不均一性[D]. 上海: 同济大学, 2018.

    GUO Yulong. Geochemical approaches for compositional heterogeneity of sediments from rivers in East China[D]. Shanghai: Tongji University, 2018.

    [33] 杨守业, 王中波. 长江主要支流与干流沉积物的REE组成[J]. 矿物岩石地球化学通报, 2011, 30(1):31-39. [YANG Shouye, WANG Zhongbo. Rare Earth Element compositions of the sediments from the major tributaries and the main stream of the Changjiang River [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2011, 30(1): 31-39. doi: 10.3969/j.issn.1007-2802.2011.01.005
    [34] 罗超, 郑洪波, 吴卫华, 等. 长江河水87Sr/86Sr值的季节性变化及其指示意义: 以长江大通站为例[J]. 地球科学进展, 2014, 29(7):835-843. [LUO Chao, ZHEGN Hongbo, WU Weihua, et al. Temporal variation in Sr and 87Sr/86Sr of Yangtze river: an example from datong hydrological station [J]. Advances in Earth Science, 2014, 29(7): 835-843.
    [35] 张玉泉, 谢应雯, 李献华, 等. 青藏高原东部钾玄岩系岩浆岩同位素特征: 岩石成因及其构造意义[J]. 中国科学(D辑), 2000, 30(5):493-498. [ZHAGN Yuquan, XIE Yingwen, LI Xianhua, et al. Isotopic characteristics of magmatic rocks of the potassium basaltic series in the eastern Qinghai-Tibet Plateau: petrogenesis and tectonic significance [J]. Science in China: Geoscience, 2000, 30(5): 493-498. doi: 10.3321/j.issn:1006-9267.2000.05.007
    [36] 张招崇, 王福生. 峨眉山玄武岩Sr、Nd、Pb同位素特征及其物源探讨[J]. 地球科学—中国地质大学学报, 2003, 28(4):431-439. [ZHANG Zhaochong, WANG Fusheng. Sr, Nd and Pb isotopic characteristics of Emeishan basalt province and discussion on their source region [J]. Earth Science——Journal of China University of Geosciences, 2003, 28(4): 431-439.
    [37] 孟宪伟, 杜德文, 陈志华, 等. 长江、黄河流域泛滥平原细粒沉积物87Sr/86Sr空间变异的制约因素及其物源示踪意义[J]. 地球化学, 2000, 29(6):562-570. [MENG Xianwei, DU Dewen, CHEN Zhihua, et al. Factors controlling spatial variation of 87Sr/86Sr in the fine-grained sediments from the overbanks of the Yellow River and Yangtze River and its implication for provenance of marine sediments [J]. Geochimica, 2000, 29(6): 562-570. doi: 10.3321/j.issn:0379-1726.2000.06.008
    [38]

    Nesbitt H W, Markovics G, Price R C. Chemical processes affecting alkalis and alkaline earths during continental weathering [J]. Geochimica et Cosmochimica Acta, 1980, 44(11): 1659-1666. doi: 10.1016/0016-7037(80)90218-5

    [39]

    Bi L, Yang S Y, Zhao Y, et al. Provenance study of the holocene sediments in the Changjiang (Yangtze River) Estuary and inner shelf of the East China Sea [J]. Quaternary International, 2016, 441: 147-161.

    [40]

    Bentley S J Jr, Blum M D, Maloney J, et al. The Mississippi river source-to-sink system: perspectives on tectonic, climatic, and anthropogenic influences, Miocene to anthropocene [J]. Earth-Science Reviews, 2016, 153: 139-174. doi: 10.1016/j.earscirev.2015.11.001

    [41]

    Hoffmann T, Thorndycraft V R, Brown A G, et al. Human impact on fluvial regimes and sediment flux during the Holocene: Review and future research agenda [J]. Global and Planetary Change, 2010, 72(3): 87-98. doi: 10.1016/j.gloplacha.2010.04.008

    [42]

    Walling D E, Fang D. Recent trends in the suspended sediment loads of the world's rivers [J]. Global and Planetary Change, 2003, 39(1-2): 111-126. doi: 10.1016/S0921-8181(03)00020-1

    [43] 赵一阳, 鄢明才, 李安春, 等. 中国近海沿岸泥的地球化学特征及其指示意义[J]. 中国地质, 2002, 29(2):181-185. [ZHAO Yiyang, YAN Mingcai, LI Anchun, et al. Geochemistry of muds along the coast of China and their significance [J]. Geology in China, 2002, 29(2): 181-185. doi: 10.3969/j.issn.1000-3657.2002.02.014
  • 期刊类型引用(1)

    1. 潘峰,蒋裕强,蒋立伟,王玥. 川南大安龙马溪组电成像测井资料中大型滑塌体的发现及地质意义. 古地理学报. 2024(04): 972-984 . 百度学术

    其他类型引用(0)

图(5)  /  表(4)
计量
  • 文章访问数:  2859
  • HTML全文浏览量:  527
  • PDF下载量:  64
  • 被引次数: 1
出版历程
  • 收稿日期:  2019-05-05
  • 修回日期:  2019-05-20
  • 网络出版日期:  2020-04-28
  • 刊出日期:  2020-05-31

目录

/

返回文章
返回