-
水系演化是地球动力系统的重要组成部分，在地球表层系统中扮演着非常重要的角色[1]。水系重建对于理解全球变化和区域响应至关重要[2⁃3]。发源于青藏高原东北缘的东亚大水系，长江和黄河的水系演化得到广泛的关注，尽管仍存在争议，但已取得许多进展[4⁃9]。
作为中国七大水系之一，东北地区最重要的一条河流，松花江水系重建对理解该区域的构造—地貌—季风演化具有重要意义。松花江水系演化研究极为薄弱，较为集中在20世纪70—90年代[10⁃14]，但研究方法主要集中在地貌学方面，同位素地球化学物源示踪研究较为薄弱。在水系演化的关键科学问题上，第四纪松花江水系是否发生倒转一直存在争议，且缺乏直接的指示证据。
沉积物作为河流沉积作用的直接产物，是研究水系演化的首选对象。国际上已展开大量的Sr⁃Nd同位素地球化学研究[15⁃18]，由于其在表生环境中的稳定性，已广泛应用于陆地和海洋沉积物物源追踪[19⁃23]。尤其是Nd同位素在地球表生物质循环过程中具有较高的稳定性而更可靠的运用于沉积物的物源研究中。哈尔滨荒山位于松花江右岸的二级阶地[24]，其沉积物物源演化信息可以为松花江形成和演化研究尤其是对于松花江水系是否发生倒转提供直接证据。本文拟通过荒山第四纪岩芯沉积物的磁化率、古地磁和Sr⁃Nd同位素组成，分析岩芯沉积物的Sr⁃Nd同位素随深度的变化特征，以期揭示荒山岩芯沉积物源的变化，重建第四纪松花江水系演化过程。这项研究不仅有助于加深理解松花江水系演化对构造—地貌—气候变化的响应过程，为东亚水系演化贡献新成果，而且也为探索东亚水系演化的一致性和区域差异性提供研究实证。
-
哈尔滨地区位于松嫩平原东部，松花江中游地区，属于干旱—半干旱温带大陆性气候，处于松辽拗陷的中部东南端、东南隆起区边缘地带，由松花江及其支流的河漫滩、河流阶地以及东部山地洪积—冲积台地构成[25]，多为沉积地层，较少发育侵入岩。其一级阶地[25]分布在市区火车站、康庄桥以及郊区的顾乡屯、王家店、四方台等地，平均海拔高度为140~150 m；二级阶地[25]分布于市区香坊、南岗以及郊区的顾乡屯砖瓦厂、东长岭子和西长岭子等地，海拔达155~175 m。
荒山（又称天恒山）位于哈尔滨市道外区团结镇东郊（图1），北望松花江，西临阿什河，属于松花江二级河流堆积阶地，最高海拔约180 m，主要由黄土和河湖相沉积物组成。荒山剖面第四纪沉积物出露厚度约53.41 m，约46.3 m被埋藏，沉积厚度较深，沉积序列较完整、连续，基底为中生代白垩纪泥岩、砂岩风化壳[26]。在荒山第四纪地层中发现有古冰缘现象以及新构造运动形成的特殊现代地貌[25]，自下而上依次划分为白土山组、荒山组、哈尔滨组，是东北地区第四纪研究的典型剖面。
图 1 松嫩平原松花江水系及钻孔、取样点分布图（红色箭头表示现在的水流方向）
Figure 1. DEM (Digital Elevation Model) map of the Songnen Plain, showing the distribution of the Songhua River drainage system and the borehole sampling point location in this study (the red arrows indicate the current flow direction)
-
荒山钻孔岩芯（45°47′32.73″ N；126°47′43.75″ E）位于天恒山公园，采用双管单动内衬塑料套管取芯技术钻取，该岩芯已打穿整个第四纪直至白垩系基岩，钻井深度101.11 m，最终获得岩芯长度93.3 m，取芯率在92%以上。对钻取的岩芯使用链锯锯开岩芯外的塑料外壳，并去除最外层被钻井泥浆污染的部分。对处理完的岩芯先记录岩性并完成对古地磁样品的定向取样。岩芯柱在实验室按照2.5 cm间距进行分割，共得到3 673个样品。
磁化率样品的测试间距为10 cm，共得到912个样品。样品在室温下风干后用手轻轻地粉碎，放入2 cm×2 cm×2 cm立方的弱磁盒。磁化率试验样品送至中国科学院地球环境研究所黄土与第四纪地质国家重点实验室，利用Bartington MS2仪器在470 Hz的低频下测得低频磁化率。
对整个岩芯以40~100 cm（岩芯0~52 m）或10~50 cm（岩芯52~101 m）为取样间隔进行古地磁定向取样两套（共406块）。剩磁测量采用安装在零磁空间（<300 nT）的2G760三轴超导磁力仪上完成。对荒山钻孔的191件样品进行系统的交变退磁，采用不包括天然剩磁（NRM）的16 步退磁，最大退磁场为80 mT。
以40~60 cm间距共获取Sr⁃Nd同位素组成样品138件。我们没有采用分粒级处理，主要是考虑分粒级样品可能会产生有偏差的结果，并可能丢失某些源区信息[27]。样品在室内自然风干后经玛瑙研钵研磨成粉末，然后过200目标准分样筛，获取<73 μm组分以供测试。尽管沉积物的Sr同位素组成易受化学风化的影响，并明显影响其的物源解释[27]，但我们并没有对样品进行酸处理，主要基于以下考虑[27]：酸处理会溶解来自母岩的碎屑方解石，从而丢失一些物源信息，影响物源解释；样品的长时间酸浸泡也会对黏土矿物造成破坏。
Sr和Nd的分离采用标准离子交换方法。测定的样品Sr和Nd同位素比值分别用86Sr/88Sr=0.119 4 和146Nd/144Nd=0.721 9进行了标准化以进行质量分馏效应的校正。仪器的准确度分别采用国际标样NBS 987和JNdi⁃1进行检测，测定值分别为87Sr/86Sr= 0.710 275±10(2σ)和143Nd/144Nd=0.512 119±8(2σ)。化学分析空白：对Sr为<1 ng，对Nd为<50 pg。Nd同位素组成用εNd(0)值表示，采用143Nd/144Nd=0.512 638的球粒陨石均一源（CHUR）现代值。Sr⁃Nd同位素的测定在中国地质大学(武汉)地质过程与矿产资源国家重点实验室完成。
为了更好分析岩芯沉积物源，在松花江哈尔滨段上、下游的主流及支流进行了样品采样。上游河段在嫩江河流阶地采样6件，雅鲁河采样1件，拉林河采样11件，松花江哈尔滨段采样3件，松花江吉林段采样8件，库勒河采样4件；下游河段在岔林河、蚂蜒河、牡丹江和通河松花江T3河流阶地各采样1件。
-
以30.4 m为界，岩芯岩性整体可分为两个部分，以上为现代土壤和黄土—古土壤序列，以下为河湖相沉积物。基于岩芯沉积物的颜色、粒度、结构和构造特征，荒山岩芯的岩性描述如下（图2）：
0~0.98 m：暗褐色的现代土壤，具备现代植物根茎，受现代人类活动影响较大。
0.98~30.4 m：黄土—古土壤层序列，包括5个黄土层（L1~L5）和5个古土壤层（S0~S4）。古土壤层为深灰褐色—灰黑色，结构致密，发育白色菌丝体。在S0古土壤层中发现炭屑、铁锈色斑和铁锈色的颗粒。黄土层为浅黄褐色，结构疏松，无黏性，含有白色菌丝体。其中26.5~30.4 m为受流水改造的L5次生黄土。
30.4~57.08 m：以浅黄褐色—灰褐色泥和粉砂质泥为主的沉积，结构致密，有黏性，微细层理发育充分，具有铁锈颗粒和黑色铁锰颗粒。顶部为粉砂质泥，中间夹杂有褐色泥沉积，在36.07~36.89 m深处发育泥球、泥粒。52 m处发现碳化明显的动物骨骼。是曲流河河漫滩相堆积。
57.08~62.3 m：以棕黄色细—极细砂为主的沉积，顶部发育泥质砂夹层，结构疏松，杂砂结构，铁染明显。是曲流河河床边滩亚相沉积。
62.3~66.13 m：以青灰色的泥—粉砂质泥为主的沉积，具有明显铁染，顶部为砂质泥过渡层，底部为泥质砂，而65.08 m以下粒度呈现逐渐变粗的变化趋势。为河漫湖泊沉积[31]。
66.13~73.39 m：以青灰色、灰绿色中—粗砂为主的沉积，夹杂有铁染严重的泥质砂，含有少量的砾石（a=8 mm，次棱角状，成分为脉石英）；70.38 m以下为粒度逐渐变粗的黄白色粗砂。为曲流河主流线附近向河漫滩过渡的边滩堆积[31⁃32]。
73.39~95.18 m：以灰黄、灰绿色、青灰色含砾粗砂为主的沉积，夹杂有两层灰绿色、青灰色粉砂质泥—泥质砂。砾石粒径多集中在10~20 mm，最大可达40 mm。是曲流河主流线附近的心滩堆积[31⁃32]。
95.18~101.11 m：以青灰、灰绿色细砂—极细砂为主的沉积，局部夹杂粉砂质泥、细砂—极细砂、灰黑色泥发育。为水动力较弱的浅湖相沉积。
综上所述，荒山岩芯沉积相反映了其由湖泊相（95.18~101.11 m）至河流相（30.4~95.18 m）再到风尘黄土堆积（0.98~30.4 m）的变化，而30.4~95.18 m河流相沉积进一步可以划分为两个完整的沉积旋回：第一沉积旋回（62.3~95.18 m）为河床滞留堆积（73.39~95.18 m）向边滩（66.13~73.39 m）再向河漫湖泊（62.3~66.13 m）转变的沉积环境；第二个旋回（30.4~62.3 m）依次为边滩（57.08~62.25 m）和河漫滩（30.4~57.08 m）的沉积环境。
-
根据荒山岩芯磁化率随深度的变化特征和岩芯岩性，我们将岩芯地层进一步划分为4段：
（1） 0~30.4 m（Ⅰ段）：该段地层为现代土壤和黄土—古土壤序列，其中黄土—古土壤序列磁化率呈现周期性髙低值变化；而现代土壤层磁化率呈现为高值，最高达94.70×10-8 m3kg-1。该段地层磁化率平均值为32.22×10-8 m3kg-1。
（2） 30.4~57.1 m（Ⅱ段）：该段地层磁化率变化特征表现为自下至上小幅度波动并逐渐升高；至50 m以上表现为周期性的髙低值变化。该段地层磁化率平均值为14.64×10-8 m3kg-1。
（3） 57.1~62.3 m（Ⅲ段）：该段地层磁化率急剧增大，波动剧烈，最高甚至达到99.67×10-8 m3kg-1。平均值为35.69×10-8 m3kg-1。
（4） 62.3~101 m（Ⅳ段）：该段地层磁化率值较低，多次出现0值；波动幅度较小。平均值为3.60×10-8 m3kg-1。
由于钻孔的岩芯已经发生转动，磁偏角并无实际指导意义，故使用磁倾角建立磁极性柱。结果表明（图2），荒山岩芯钻孔具有较为规律的极性顺序，和标准地磁极性柱[28⁃30]有良好的对比性。钻孔的Brunhes正极性带和Matuyama负极性带界限比较清晰，位于钻孔51.4 m处，时间界限为0.78 Ma。根据磁极性序列该钻孔自上至下依次划分为8个极性带，包括4个正极性和4个负极性。51.4 m以上为Brunhes正极性带；65.71~69.82 m为Jaramillo正极性带；77.46~77.88 m为Cobb Mountain正极性带；95.24~95.44 m为Gilsa正极性带。通过对岩芯极性柱的线性内插推得，荒山岩芯62.3 m处的地质年龄为0.94 Ma B.P.。
-
荒山岩芯的87Sr/86Sr在0.709 043~0.714 077变化，平均值为0.711 557；εNd(0)在-8.7~-4.03变化，平均值为-6.09。根据岩芯Sr⁃Nd同位素组成随深度变化曲线和岩芯的岩性（图3），可以将岩芯Sr⁃Nd同位素组成进一步划分为四段：
图 3 荒山岩芯Sr⁃Nd同位素组成随深度变化曲线
Figure 3. Down⁃core variations of Sr⁃Nd isotope composition for the HS core;
（1） 0~30.4 m（Ⅰ段）：黄土—古土壤序列的87Sr/86Sr值为0.710 465~0.712 213，平均值为0.711 399；εNd（0）为-6.9~-3.7，平均值为-5.9。该地层15.1 m以上地层Sr⁃Nd同位素组成均表现为在岩芯均值上下小幅度的变化特征；15.1 m以下地层的Sr同位素表现为小于岩芯均值的周期性变化特征；Nd同位素表现为大于岩芯均值的周期性变化特征。
（2） 30.4~57.1 m（Ⅱ段）：该段地层87Sr/86Sr值介于0.711 963~0.713 528，平均值为0.712 738。εNd(0)值介于-7.2~-5.9，平均值为-6.6。Sr⁃Nd同位素组成变化幅度较小，其中Sr同位素表现为大于岩芯均值；Nd同位素表现为小于岩芯均值。
（3） 57.1~62.3 m（Ⅲ段）：该段地层Ⅲ段地层87Sr/86Sr值介于0.709 043~0.712 852，平均值为0.710 048；εNd(0)值介于-8.7~-4.4，平均值为-6.0。其Sr⁃Nd同位素组成变化特征表现为在岩芯均值上下剧烈波动，且具有Ⅱ段和Ⅳ段的Sr⁃Nd同位素组成特征。
（4） 62.3~101 m（Ⅳ段）：该段地层87Sr/86Sr值组成介于0.709 049~0.714 077，平均值为0.711 056；εNd(0)值变化范围为-8.5~-4.0，平均值为-6.0。其Sr⁃Nd同位素组成变化特征表现为在岩芯均值上下大幅度变化，且Sr同位素整体低于岩芯均值；Nd同位素整体高于岩芯均值。
-
松花江水系河流冲积物Sr⁃Nd同位素组成测试结果（表1）显示，松花江松原方向支流的样品87Sr/86Sr值变化范围为0.707 394~0.711 336，平均值为0.709 784；Nd同位素组成εNd(0)值变化范围为-11.0~-5.1，平均值为-8.1。
表 1 荒山岩芯和松花江水系主要河流沉积物Sr⁃Nd同位素组成
Table 1.
Comparison of Sr⁃Nd isotope composition between the HS core and sediments from the Songhua River drainage 流域 87Sr/86Sr 143Nd/144Nd εNd(0) 荒山 岩芯 荒山Ⅱ段 0.711 963~0.713 528 0.512 27~0.512 334 -7.2~-5.9 荒山Ⅲ段 0.709 043~0.712 852 0.512 194~0.512 410 -8.7~-4.4 荒山Ⅳ段 0.709 049~0.714 077 0.512 204~0.512 431 -8.5~-4.0 上游 嫩江肇源 0.708 014~0.709 607 0.512 239~0.512 359 -7.8~-5.4 第二松花江 0.710 479~0.711 243 0.512 113~0.512 222 -10.2~-8.1 雅鲁河 0.708 322 0.512 300 -6.6 拉林河 0.710 158~0.711 336 0.512 076~0.512 278 -11.0~-7.0 松花江哈尔滨段 0.709 521~0.709 665 0.512 084~0.512 204 -10.8~-8.5 库勒河 0.707 394~0.710 277 0.512 274~0.512 379 -7.1~-5.1 中游 牡丹江 0.712 026 0.512 14 -9.7 岔林河 0.709 731 0.512 448 -3.7 蚂蜒河 0.709 018 0.512 569 -1.3 通河T3阶地 0.715 166 0.512 408 -4.5 松花江依兰方向支流样品87Sr/86Sr值变化范围为0.709 018~0.712 026，平均值为0.711 485；Nd同位素组成εNd(0)值变化范围为-9.7~-1.3，平均值为-4.8。
-
沉积物磁化率值的大小主要决定于铁磁性矿物的浓度、颗粒的大小和磁性矿物的成分，反映了沉积物物源、沉积动力大小及次生条件等沉积环境的变化[33]。荒山岩芯62.3 m以下地层磁化率与62.3 m以上地层存在较大差异，Ⅳ段（62.3~101 m）地层与Ⅲ地层（57.1~62.3 m）的磁化率均值差异达到10倍；和Ⅱ段（30.4~57.1 m）地层磁化率平均值差异也达到5倍之多。沉积动力以及次生条件的变化难以解释如此急剧的磁化率波动。沉积物的磁学特征与物源区出露的岩石类型有着密切的联系，荒山岩芯过渡段磁化率的剧烈波动应是由物源的变化引起的，而非次生条件的变化所导致。
在物源示踪方面，Sr⁃Nd同位素组成已成功地应用于长江、恒河等流域物源的约束[2,34⁃35]，这不仅因为源区不同年代地质体具有不同Sr⁃Nd同位素组成，更是由于Sr⁃Nd同位素组成在风化和搬运等过程中所受影响不大，在地球表生循环过程中比较稳定[19⁃23]。尤其是Nd同位素，在地球表生循环过程中基本只受Sm元素在壳幔分异过程的影响，可以更可靠地运用于沉积物的物源指示中[19,35⁃37]。
由于岩芯30.4 m以上为现代土壤和风成黄土—古土壤系列，其Sr⁃Nd同位素组成变化不在我们水系演化的讨论之中。岩芯Ⅱ段（30.4~57.1 m）段和Ⅳ段（62.3~101 m）沉积物具有明显不同的Sr⁃Nd同位素组成。岩芯62.3 m以下地层Sr⁃Nd同位素组成在岩芯均值上下剧烈变化，Sr同位素整体表现为低于岩芯均值，Nd同位素整体则表现为高于岩芯均值。57.08 m以上地层Sr⁃Nd同位素组成变幅较小，Sr同位素高于岩芯均值，呈现为递减的变化趋势，Nd同位素低于岩芯均值，呈现递增的变化趋势。岩芯62.3 m上、下地层Sr⁃Nd同位素组成的这种差异反映了荒山地区0.94 Ma前后沉积物源并不一致，其物源区存在较大变化。
-
第四纪松花江水系演化模式一直存在争议。杨秉赓、初本君等人认为早更新世松花江水流模式为西流模式[11⁃13,38]，反转时间笼统的指向早更新世末或中更新世初[11,13]，但是反转的原因存在争议。杨秉赓等人认为是向源侵蚀引起的河流袭夺和水系倒转[11,39]，初本君等[12]认为是由三江平原的沉降改使松嫩古湖东外泄至三江盆地[12]。然而，裘善文等认为第四纪以来松花江一直是东流模式[10,40⁃41]，由松嫩湖盆的东北出口东流至三江平原。更早的古近纪和新近纪时期，上松花江和牡丹江流向达连河和祥顺地区的一古大湖[42]，即松花江从未西流向松嫩平原。尽管各家观点不同，但是均认为松花江水系在早更新世晚期和中更新世之前在依兰和佳木斯区域存在佳依（佳木斯—依兰）分水岭[11⁃13,43⁃44]。
从荒山岩芯的岩性来看，73.39~95.18 m地层岩性为含砾粗砂，66.13~73.39 m为中—粗砂，62.25~66.13 m为亚黏土和黏土，表现出自下而上沉积物粒度由粗变细的沉积过程，这反映了受构造和地貌控制的水动力条件逐渐减弱的变化过程。已有研究表明，晚第三纪末松辽平原受新构造运动的影响，山区抬升较平原上升显著[11]。佳依分水岭也随之抬升，导致河流向源侵蚀加剧，产生大量新鲜碎屑物质。砾石多为造山作用的产物[45⁃46]，砾石的出现说明该时期荒山岩芯沉积物来自构造抬升的山地。哈尔滨以西的松原一带地处松嫩平原腹地，地势平坦，且发源于大兴安岭的嫩江也无法将砾石远距离搬运至哈尔滨地区进行堆积，松花江哈尔滨段以及松原至哈尔滨段的现代河床及河流阶地沉积物皆为中—细砂，并无砾石堆积。所以荒山岩芯底部的砾石不可能由松原方向搬运过来，只能来自哈尔滨以东的小兴安岭和张广才岭。同时，由于佳依分水岭的构造抬升使得松花江依兰—哈尔滨河段河流纵比降加大，水流搬运能力强，河流将大量砾石搬运至哈尔滨地区，所以岩芯底部出现砾石沉积。随着河流沉积物在下游的堆积，下游河床被抬高，而上游河床由于河流的下切侵蚀则不断被蚀低，导致依兰—哈尔滨河段河床的相对高差逐渐减小，河流纵比降减小，水流搬运能力下降，河流所搬运的沉积物粒度逐渐变细，造成岩芯下部沉积物的粒度由下至上逐渐变细。因此荒山岩芯下部沉积物的粒度变化规律在一定程度上证明了松花江的西流模式。
荒山岩芯沉积物Sr⁃Nd地球化学组成随深度的变化显示（图3），岩芯Ⅳ段（62.3 m以下）Sr⁃Nd同位素组成呈现高频高幅的变化特征。造成这种现象的原因是，河流的持续下蚀，使河谷两侧及底部的基岩不断暴露并受到流水侵蚀，导致大量新鲜碎屑物质加入河流载荷，这些新生碎屑物质具有低的再循环特征和较高的Nd同位素组成。“新”与“老”碎屑物质的混合导致Sr⁃Nd同位素组成在岩芯中的剧烈波动变化。河流的下切侵蚀往往发育在纵比降较大的河流上游的构造隆升地区。现今松花江上游地区流经松嫩平原，地势较为平坦，并不具备下切侵蚀的条件，而依兰—佳木斯一带在早更新世为构造抬升地区，被小兴安岭和张广才岭所围限，地势较高，地层广泛出露，可以为河流下蚀作用提供条件并提供大量碎屑物质进入河流沉积物。
荒山岩芯Ⅳ段的Nd同位素组成εNd(0)的最大值为-4.0（表1），达到岩芯剖面的最大值。松花江水系松原方向冲积物的Nd同位素组成εNd(0)偏负，最大值仅为-5.1，远小于岩芯Ⅳ段的Nd同位素组成的最大值，而依兰方向河流冲积物的Nd同位素组成εNd(0)明显偏正，除了牡丹江外，其他河流冲积物在-1.3~-4.5之间。由此可见，松原方向水系冲积物不可能为岩芯Ⅳ段提供如此高εNd(0)值的碎屑物质，而依兰方向水系冲积物可以为岩芯Ⅳ段提供高εNd(0)值的碎屑物质。因此我们认为，岩芯62.3 m以下地层沉积物来源于依兰方向，而62.3 m以上地层沉积物来源于松原方向。
松嫩平原属于周边向中心倾斜的半封闭不对称盆地，四周被大兴安岭、小兴安岭和长白山包围。早更新世佳依分水岭被切穿前，水流顺应地貌梯度从盆地边缘向心状地流向盆地中心。所以从地貌角度看，在佳依分水岭被切穿前，依兰方向应该是松花江的上游。另外，在通河县野外考察中发现的松花江T3堆积阶地具有明显的二元结构，上部为厚度约2.6 m棕红色河漫滩相含砾粉砂质泥，下部为厚度约5.6 m的灰白色河床相砂砾石层（图4a）。砾石层粒径多集中在20~50 mm，最大粒径达260 mm（图4b），磨圆度较好，有较好的定向排列。这些砾石是佳依分水岭被切穿之前河流上游的产物。首先，在松嫩平原只有山区由于受新构造抬升影响而存在T3阶地，且多数为堆积阶地而少数为基座阶地；而哈尔滨地区地势平坦，受松嫩平原沉降作用的影响，只能形成T1（全新世）和T2（形成于晚更新世，大约在40 ka B.P.，未出版数据）河流堆积阶地。其次，砾石的大规模堆积往往反映了造山运动的山体隆升，河流持续下蚀将砾石向下游输送堆积形成[45⁃46]。若这套砾石层来源于哈尔滨方向，那么位于通河上游的哈尔滨地区应该存在粒径更大的砾石堆积，然而松原至哈尔滨一带的地表不见砾石堆积，且荒山岩芯底部仅存在细颗粒的砾石。最重要的是，通河T3阶地砾石的排列方向显示（图4c），该砾石层的砾石ab面倾向东南，倾向基本在120°~150°之间（图4c），表明了自东向西的古水流流向。所以该套砾石层的存在也直接证明了当时松花江是从依兰方向向西流至哈尔滨方向。
图 4 通河松花江T3阶地的砾石堆积
Figure 4. Gravel deposit accumulation on the Songhua River T3 terraces in Tonghe county.
尽管裘善文等[10,40⁃41]认为松花江水流只能是向东注入通河—依兰之间的古大湖，但他们承认佳依分水岭的存在。然而，裘善文等的水系模式存在一个问题，即佳依分水岭被切穿之前，东流至通河—依兰古大湖的松花江水流的出口问题。由于佳依分水岭的存在，松花江东流注入通河—依兰古大湖后必然要有出口，但是，在佳依分水岭被切穿之前这个出口是不存在的。因此，在佳依分水岭被切穿前松花江根本不可能向东流。
综上所述，我们提出如下的松花江水系演化模式（图5）。在早更新世早—中期，松花江中上游与下游尚未贯通。以佳依分水岭为界，松花江下游从佳木斯向东流经三江平原，松花江中上游（即，肇源—依兰河段）自依兰向西流入松嫩湖盆。同时，嫩江、松花江吉林段、东辽河和西辽河分别注入松嫩古湖，古湖南侧有一个缺口，通过辽河，水流流经下辽河裂谷外泄至渤海。早更新世晚期，佳依分水岭构造抬升，分水岭两侧的松嫩平原和三江平原明显沉降，河流纵比降不断增大，分水岭两侧河流不断向源侵蚀；同时受0.9 Ma中更新世气候转型的影响[47⁃49]，冬季风加强，气温降低，全球冰量增加，侵蚀基准面下降，河流的向源侵蚀加剧。最终导致佳依分水岭在0.94 Ma B.P.被切穿，松花江中上游水系被下游水系所袭夺，流向发生反转，松花江中上游河段自西向东流经佳依峡谷进入三江平原，松花江中上游与下游河段得以贯通。与此同时，受松花江反向以及松辽分水岭抬升的影响，嫩江和松花江吉林段河道被迫改道最终并入松花江。而由于下辽河的持续向源侵蚀以及松辽分水岭的抬升，下辽河袭夺东、西辽河，松嫩古湖不再有河流水流注入，逐渐消失。
图 5 松花江水系演化示意图
Figure 5. Reconstruction of the Songhua River system evolution
-
通过对荒山岩芯的沉积相、磁化率、古地磁和Sr⁃Nd同位素的测定对第四纪以来松花江水系的演化，得到以下结论：
（1） 荒山岩芯沉积相反映了其由湖泊相（95.18~101.11 m）至河流相（30.4~95.18 m）再到风尘黄土堆积（0.98~30.4 m）的变化，而河流相沉积进一步可以划分为由河床滞留堆积向边滩再向河漫湖转变的第一个沉积旋回和由边滩向河漫滩沉积环境转变的第二个沉积旋回。
（2） 荒山岩芯磁化率以62.3 m为界限，上下两段地层磁化率存在明显差异。62.3 m以下地层磁化率值较低，多处出现0值，波动较小；62.3 m以上地层磁化率在急剧升高，并呈现周期性的髙低变化。
（3） 荒山岩芯钻孔的B/M界限位于钻孔51.4 m处，时间界限为0.78 Ma。自51.4 m以下包含有Jaramillo、Cobb Mountain和Gilsa正极性带，通过对岩芯极性柱的线性内插推得，荒山岩芯62.3 m处的地质年龄为0.94 Ma B.P.。
（4） 荒山岩芯62.3 m以下地层的Sr⁃Nd同位素组成剧烈变化，其中Sr同位素组成整体表现为低于岩芯均值，Nd同位素整体表现为高于岩芯均值；岩芯30.4~57.1 m段地层的Sr⁃Nd同位素组成变幅较小，其中Sr同位素高于岩芯均值，Nd同位素低于岩芯均值。表明了岩芯上下两段沉积物源发生了明显变化。
（5） 荒山岩芯62.3 m以下地层的Sr⁃Nd同位素组成与依兰方向水系沉积物相近，而62.3 m以上地层的Sr⁃Nd同位素组成与松原方向水系沉积物相近。表明62.3 m（0.94 Ma B.P.）前后沉积物物源和河流流向发生了显著变化。
（6） 早更新世晚期之前，以佳依分水岭为界，松花江下游向东流经三江平原，而松花江中上游河段向西流入松嫩湖盆。早更新世晚期，佳依分水岭构造抬升，松嫩平原与三江平原不断下沉，同时受中更新世气候转型的影响，全球气温降低，侵蚀基准面下降，河流向源侵蚀加剧。在构造—气候耦合作用的影响下，最终导致佳依分水岭在0.94 Ma B.P.左右被切穿，松花江中上游水系被下游水系所袭夺，河流发生反转，松花江水系得以贯通，逐渐形成现代松花江水系格局。
The Inversion of the Songhua River System in the Early Pleistocene: Implications from Sr⁃Nd isotopic composition in the Harbin Huangshan cores
-
摘要:
水系重建对于理解全球变化和区域响应至关重要。松花江水系演化研究极为薄弱，尤其对第四纪早期松花江水系是否发生倒转一直存在争议，且无明确证据。河流沉积物是水系演化最直观的证据。为此，选择哈尔滨荒山（HS）钻孔岩芯沉积物作为研究对象，对其进行了磁化率、古地磁和Sr⁃Nd同位素组成分析。结果表明，岩芯62.3 m（0.94 Ma B.P.）深度上、下地层的磁化率和Sr⁃Nd同位素特征存在较大差异。62.3 m以下地层磁化率极低，多次出现0值，且变幅较小。Sr⁃Nd同位素组成与依兰方向现代水系的Nd同位素特征相近；然而62.3 m以上地层磁化率急剧增大，并呈现周期性的高低变化。Sr⁃Nd同位素组成小幅度变化，表现出与松原方向现代水系相近的Nd同位素特征。表明了早更新世晚期松花江中上游（肇源—依兰河段）水系流向曾发生反转。早更新世早—中期，松花江中上游与下游尚未贯通。以佳依（佳木斯—依兰）分水岭为界，松花江下游向东流经三江平原，而松花江中上游向西流入松嫩湖盆；早更新世晚期，由于构造—气候变化的共同作用，佳依分水岭不断抬升，而三江平原和松嫩平原持续下降，导致分水岭两侧河流向源侵蚀，在0.94 Ma B.P.，佳依分水岭被切穿，松花江中上游水系被下游水系所袭夺，流向开始倒转，中上游与下游贯通，现代松花江水系逐渐建立。 Abstract:Drainage evolution reconstruction is crucial for understanding global change and the regional response. As the most important river in northeast China， the study of the evolution of the Songhua River system is extremely weak. Whether the water system of the Songhua River reversed in the early Pleistocene has been especially controversial and without clear evidence. As a direct product of river geological processes， sediment is the most direct evidence of the water system evolution. Located in the east Songnen Plain and the middle reaches of the Songhua River， Harbin has an arid⁃semiarid temperate continental climate. Located in the eastern suburb of Tuanjie town， Daowai district， Harbin city， the Huangshan （HS） belongs to the second⁃order terrace on the right bank of the Songhua River geomorphology The evolution information of its sediment provenance can provide direct evidence for the study of the formation and evolution of the Songhua River， especially for the inversion of the Songhua River system. Therefore， magnetic susceptibility， paleomagnetic， and Sr⁃Nd isotope composition analysis of sediments from boreholes core in the HS are performed. The results showed that the magnetic susceptibility and Sr⁃Nd isotope characteristics of the upper and lower strata at a core depth of 62.3 m were significantly different. The magnetic susceptibility of the strata below 62.3 m is extremely low， and 0 value appears many times with only a small variation in range. The Sr⁃Nd isotope composition differ sharply above and below the mean value of the core. However， the magnetic susceptibility of the strata above 62.3 m increased sharply， and showed periodic high and low changes. And the isotopic composition of Sr⁃Nd varied slightly. Through the comparison with the Nd isotope composition of the upper and lower strata of 62.3 m in the rock core of Huangshan Mountain， we observe that the Nd isotope characteristics of the strata below 62.3 m are similar to those of the modern river system in the direction of Yilan⁃Tonghe. The Nd isotope characteristics of the Nenjiang River， second Songhua River， Lalin River， and Kule River systems are similar to those of the 62.3 m above. It indicates that the river direction in the middle and upper reaches of the Songhua River （Zhaoyuan⁃Yilan reach） was reversed in the late early Pleistocene.To summarize， the following evolution model of the Songhua River system is proposed： In the early and middle early Pleistocene， the upper and lower reaches of Songhua River were not connected. Bounded by the Jiayi （Jiamusi⁃Yilan） watershed， the lower reaches of the Songhua River flow eastward through the Sanjiang Plain， while the upper reaches of the Songhua River flow westward into the Songnen Lake basin. In the late early Pleistocene， due to the combined effect of tectonic⁃climate change， the Sanjiang and Songnen Plains continued to decline， while the Jiayi watershed continued to rise， leading to the erosion of rivers on both sides of the watershed to the source. At the same time， affected by 0.9 Ma in the Pleistocene climate transformation， headward erosion at 0.94 Ma B.P.， the watershed is cut to wear. The Songhua River shelter⁃forest system was captured by the downstream drainage and began to reverse the Songhua River shelter⁃forest and downstream， the Songhua River shelter⁃forest river from west to east through better in canyon in Sanjiang plain， and the Songhua River shelter⁃forest and downstream river， thus linking modern Songhua River water system gradually established. -
Key words:
- Songhua River /
- drainage evolution /
- Huangshan core /
- Sr⁃Nd isotope /
- source tracing /
- Jiayi watershed
-
图 1 松嫩平原松花江水系及钻孔、取样点分布图（红色箭头表示现在的水流方向）
Figure 1. DEM (Digital Elevation Model) map of the Songnen Plain, showing the distribution of the Songhua River drainage system and the borehole sampling point location in this study (the red arrows indicate the current flow direction)
图 3 荒山岩芯Sr⁃Nd同位素组成随深度变化曲线
实线a，b分别为荒山岩芯εNd（0）和87Sr/86Sr的平均值
Figure 3. Down⁃core variations of Sr⁃Nd isotope composition for the HS core;
lines a and b are the mean values of Nd(0) and 87Sr/86Sr of the HS core
图 4 通河松花江T3阶地的砾石堆积
（a）通河T3堆积阶地；（b）砾石a轴长度频率分布直方图；（c）古流向玫瑰花图
Figure 4. Gravel deposit accumulation on the Songhua River T3 terraces in Tonghe county.
(a) T3 constructional terrace in Tonghe; (b) gravel a axis length frequency distribution histogram; (c) paleocurrent rosette map
图 5 松花江水系演化示意图
（a）贯通前水系格局 （b）贯通后水系格局
Figure 5. Reconstruction of the Songhua River system evolution
表 1 荒山岩芯和松花江水系主要河流沉积物Sr⁃Nd同位素组成
Table 1.
Comparison of Sr⁃Nd isotope composition between the HS core and sediments from the Songhua River drainage 流域 87Sr/86Sr 143Nd/144Nd εNd(0) 荒山 岩芯 荒山Ⅱ段 0.711 963~0.713 528 0.512 27~0.512 334 -7.2~-5.9 荒山Ⅲ段 0.709 043~0.712 852 0.512 194~0.512 410 -8.7~-4.4 荒山Ⅳ段 0.709 049~0.714 077 0.512 204~0.512 431 -8.5~-4.0 上游 嫩江肇源 0.708 014~0.709 607 0.512 239~0.512 359 -7.8~-5.4 第二松花江 0.710 479~0.711 243 0.512 113~0.512 222 -10.2~-8.1 雅鲁河 0.708 322 0.512 300 -6.6 拉林河 0.710 158~0.711 336 0.512 076~0.512 278 -11.0~-7.0 松花江哈尔滨段 0.709 521~0.709 665 0.512 084~0.512 204 -10.8~-8.5 库勒河 0.707 394~0.710 277 0.512 274~0.512 379 -7.1~-5.1 中游 牡丹江 0.712 026 0.512 14 -9.7 岔林河 0.709 731 0.512 448 -3.7 蚂蜒河 0.709 018 0.512 569 -1.3 通河T3阶地 0.715 166 0.512 408 -4.5 
下载: 导出CSV
-
[1] Jia J T， Zheng H B， Huang X T， et al. Detrital zircon U-Pb ages of Late Cenozoic sediments from the Yangtze delta： Implication for the evolution of the Yangtze River［J］. Chinese Science Bulletin， 2010， 55（15）： 1520-1528. [2] Shao L， Li C A， Yuan S Y， et al. Neodymium isotopic variations of the Late Cenozoic sediments in the Jianghan Basin： Implications for sediment source and evolution of the Yangtze River［J］. Journal of Asian Earth Sciences， 2012， 45： 57-64. [3] Hu Z B， Pan B T， Bridgland D， et al. The linking of the upper-middle and lower reaches of the Yellow River as a result of fluvial entrenchment［J］. Quaternary Science Reviews， 2017， 166： 324-338. [4] Clark M K， Schoenbohm L M， Royden L H， et al. Surface uplift， tectonics， and erosion of eastern Tibet from large-scale drainage patterns［J］. Tectonics， 2004， 23（1）： TC1006. [5] Richardson N J， Densmore A L， Seward D， et al. Did incision of the three gorges begin in the Eocene？［J］. Geology， 2010， 38（6）： 551-554. [6] McPhillips D， Hoke G D， Jing L Z， et al. Dating the incision of the Yangtze River gorge at the first bend using three-nuclide burial ages［J］. Geophysical Research Letters， 2016， 43（1）： 101-110. [7] Kong P， Jia J， Zheng Y. Time constraints for the Yellow River traversing the Sanmen Gorge［J］. Geochemistry， Geophysics， Geosystems， 2014， 15（2）： 395-407. [8] Li B F， Sun D H， Xu W H， et al. Paleomagnetic chronology and paleoenvironmental records from drill cores from the Hetao Basin and their implications for the formation of the Hobq Desert and the Yellow River［J］. Quaternary Science Reviews， 2017， 156： 69-89. [9] Lin A M， Yang Z Y， Sun Z M， et al. How and when did the Yellow River develop its square bend？［J］. Geology， 2001， 29（10）： 951-954. [10] 裘善文，孙广友，李卫东，等. 三江平原松花江古水文网遗迹的发现［J］. 地理学报，1979，34（3）：265-274. Qiu Shanwen， Sun Guangyou， Li Weidong， et al. The discovery of the fossil hydrographic net of Songhua Jiang in the three-river plain［J］. Acta Geographica Sinica， 1979， 34（3）： 265-274. [11] 杨秉赓，孙肇春，吕金福. 松辽水系的变迁［J］. 地理研究，1983，2（1）：48-56. Yang Binggeng， Sun Zhaochun， Jinfu Lü. The changes of Songhua-Liao river system［J］. Geographical Research， 1983， 2（1）： 48-56. [12] 初本君，高振操，杨世生，等. 黑龙江省第四纪地质与环境［M］. 北京：海洋出版社，1988：1-146. Chu Benjun， Gao Zhencao， Yang Shisheng， et al. Quaternary geology and environment of Heilongjiang province［M］. Beijing： Ocean Press， 1988： 1-146. [13] 孙广友. 松辽平原中部第四纪地壳运动与平原发育：兼论松辽分水岭的形成［M］//裘善文. 中国东北平原第四纪自然环境形成与演化. 哈尔滨：哈尔滨地图出版社，1990：44-50. Sun Guangyou. The Quaternary crustal movement and the plain development in the Central Songliao Plain， the formation of Songliao divide［M］//Qiu Shanwen. The formation and evolution of quaternary natural environment in North East China plain. Harbin： Harbin Map Press， 1990： 44-50.］ [14] 张庆云. 松辽分水岭的形成及其演化［M］//裘善文. 中国东北平原第四纪自然环境形成与演化. 哈尔滨：哈尔滨地图出版社，1990：30-33. ［ Zhang Qingyun. The formation of Songliao divide［M］//Qiu Shanwen. The formation and evolution of Quaternary natural environment in North East China plain. Harbin： Harbin Map Press， 1990： 30-33. [15] Krishnaswami S， Trivedi J R， Sarin M M， et al. Strontium isotopes and rubidium in the Ganga-Brahmaputra river system： Weathering in the Himalaya， fluxes to the Bay of Bengal and contributions to the evolution of oceanic 87Sr/86Sr［J］. Earth and Planetary Science Letters， 1992， 109（1/2）： 243-253. [16] Pande K， Sarin M M， Trivedi J R， et al. The Indus River system （India-Pakistan）： Major-ion chemistry， uranium and strontium isotopes［J］. Chemical Geology， 1994， 116（3/4）： 245-259. [17] Derry L A， France-Lanord C. Neogene Himalayan weathering history and river87Sr/86Sr： Impact on the marine Sr record［J］. Earth and Planetary Science Letters， 1996， 142（1/2）： 59-74. [18] Quade J， Roe L， DeCelles P G， et al. The Late Neogene 87Sr/86Sr record of lowland Himalayan rivers［J］. Science， 1997， 276（5320）： 1828-1831. [19] Goldstein S L， O’Nions R K， Hamilton P J. A Sm-Nd isotopic study of atmospheric dusts and particulates from major river systems［J］. Earth and Planetary Science Letters， 1984， 70（2）： 221-236. [20] Goldstein S J， Jacobsen S B. Nd and Sr isotopic systematics of river water suspended material： Implications for crustal evolution［J］. Earth and Planetary Science Letters， 1988， 87（3）： 249-265. [21] Grousset F E， Biscaye P E， Zindler A， et al. Neodymium isotopes as tracers in marine sediments and aerosols： North Atlantic［J］. Earth and Planetary Science Letters， 1988， 87（4）： 367-378. [22] Revel M， Cremer M， Grousset F E， et al. Grain-size and Sr-Nd isotopes as tracer of paleo-bottom current strength， Northeast Atlantic Ocean［J］. Marine Geology， 1996， 131（3/4）： 233-249. [23] Grousset F E， Biscaye P E. Tracing dust sources and transport patterns using Sr， Nd and Pb isotopes［J］. Chemical Geology， 2005， 222（3/4）： 149-167. [24] 叶启晓，吴金城，魏正一. 哈尔滨荒山近期的上升运动［J］. 地理科学，1984，4（4）：383-388. Ye Qixiao， Wu Jincheng， Wei Zhengyi. The tectonic uplift occured in Huang Shan in the last thousants of year［J］. Scientia Geographica Sinica， 1984， 4（4）： 383-388. [25] 朱景湖. 哈尔滨地区地貌和第四纪沉积及其在农业利用上的意义［J］. 哈尔滨师范学院学报（自然科学版），1963：188-205. Zhu Jinghu. Geomorphology and quaternary sediments in Harbin area and their significance in agricultural utilization［J］. Natural Science Journal of Harbin Normal University， 1963： 188-205. [26] 张序强. 哈尔滨荒山地区100万年来地理环境分析［J］. 哈尔滨师范大学自然科学学报，1995（3）：96-97. Zhang Xuqiang. The differentite on palaeoenvironment of Huang Shan area since 100 million years［J］. Natural Science Journal of Harbin Normal University， 1995（3）： 96-97. [27] Újvári G， Varga A， Ramos F C， et al. Evaluating the use of clay mineralogy， Sr–Nd isotopes and zircon U–Pb ages in tracking dust provenance： An example from loess of the Carpathian Basin［J］. Chemical Geology， 2012， 304-305： 83-96. [28] Cande S C， Kent D V. Revised calibration of the geomagnetic polarity timescale for the Late Cretaceous and Cenozoic［J］. Journal of Geophysical Research， 1995， 100（B4）： 6093-6095. [29] Gradstein F M， Ogg J G， Smith A G. A geologic time scale 2004［M］. Cambridge： Cambridge University Press， 2004： 1-589. [30] Laj C， Channell J E T. Geomagnetic excursions［M］//Kono M. Treatise on geophysics， geomagnetism. Amsterdam： Elsevier， 2007： 373-416. [31] Reading H G. Sedimentary environments： Processes， facies and stratigraphy［M］. 3rd ed. Oxford： Blackwell Science， 1996： 1-231. [32] Pan B T， Chen D B， Hu X F， et al. Drainage evolution of the Heihe River in western Hexi Corridor， China， derived from sedimentary and magnetostratigraphic results［J］. Quaternary Science Reviews， 2016， 150： 250-263. [33] 杨晓强，李华梅. 泥河湾盆地典型剖面沉积物磁化率特征及其意义［J］. 海洋地质与第四纪地质，1999，19（1）：75-84. Yang Xiaoqiang， Li Huamei. The susceptibility characteristics of the typical sections of the Nihewan Basin and its implication［J］. Marine Geology & Quaternary Geology， 1999， 19（1）： 75-84. [34] Awasthi N， Ray E， Paul D. Sr and Nd isotope compositions of alluvial sediments from the Ganga Basin and their use as potential proxies for source identification and apportionment［J］. Chemical Geology， 2018， 476： 327-339. [35] Li X H， Wei G J， Shao L， et al. Geochemical and Nd isotopic variations in sediments of the South China Sea： A response to Cenozoic tectonism in SE Asia［J］. Earth and Planetary Science Letters， 2003， 211（3/4）： 207-220. [36] Clift P D， Lee J I， Hildebrand P， et al. Nd and Pb isotope variability in the Indus River system： Implications for sediment provenance and crustal heterogeneity in the western Himalaya［J］. Earth and Planetary Science Letters， 2002， 200（1/2）： 91-106. [37] Weldeab S， Emeis K C， Hemleben C， et al. Provenance of lithogenic surface sediments and pathways of riverine suspended matter in the eastern Mediterranean Sea： Evidence from 143Nd/144Nd and 87Sr/86Sr ratios［J］. Chemical Geology， 2002， 186（1/2）： 139-149. [38] 刘祥，孙文丽. 西辽河水系变迁［J］. 内蒙古水利，2001（4）：70-72. Liu Xiang， Sun Wenli. Changes of Xiliao River system［J］. Inner Mongolia Water Resourcas， 2001（4）： 70-72. [39] 赵济. 中国自然地理［M］. 北京：高等教育出版社，1995：1-343. Zhao Ji. Atlas of Chinese physical geography［M］. Beijing： Higher Education Press， 1995： 1-343. [40] 裘善文，王锡魁，张淑芹，等. 松辽平原古大湖演变及其平原的形成［J］. 第四纪研究，2012，32（5）：1011-1021. Qiu Shanwen， Wang Xikui， Zhang Shuqin， et al. The evolution of the large paleolake in Songliao Plain and its formation［J］. Quaternary Sciences， 2012， 32（5）： 1011-1021. [41] 裘善文，夏玉梅，李凤华，等. 松辽平原第四纪中期古地理研究［J］. 科学通报，1984（3）：172-174. Qiu Shanwen， Xia Yumei， Li Fenghua， et al. Palaeogeography of the Songliao Plain in the Middle Quaternary［J］. Chinese Science Bulletin， 1984（3）： 172-174. [42] 裘善文，王锡魁，Makhinov A N，等. 中国东北平原及毗邻地区古水文网变迁研究综述［J］. 地理学报，2014，69（11）：1604-1614. Qiu Shanwen， Wang Xikui， Makhinov A N， et al. Summary of the paleodrainage pattern changes in the Northeast China Plain and its neighboring areas［J］. Acta Geographica Sinica， 2014， 69（11）： 1604-1614. [43] 斯米尔诺夫A M. 东北平原地形与水文网的发育［J］. 中国第四纪研究，1958，1（2）：159-162. Smirnov A M. The development of the landforms and drainage patterns in the North East China Plain［J］. Quaternaria Sinica， 1958， 1（2）： 159-162. [44] 孙肇春，钟金岳，胡俭彬，等. 东北河流阶地与河流劫夺［J］. 吉林师大学报，1964（2）：135-152. Sun Zhaochun， Zhong Jinyue， Hu Jianbin， et al. River terraces and river robberies in Northeast China［J］. Journal of Jilin Normal University， 1964（2）： 135-152. [45] 傅开道，方小敏，高军平，等. 青藏高原北部砾石粒径变化对气候和构造演化的响应［J］. 中国科学（D辑）：地球科学，2006，36（8）：733-742. Fu Kaidao， Fang Xiaomin， Gao Junping， et al. Response of gravel size change to climate and tectonic evolution in the northern Qinghai-Tibet Plateau China［J］. Science China （Seri. D）： Earth Sciences， 2006， 36（8）： 733-742. [46] 蔡观强，郭锋，刘显太，等. 碎屑沉积物地球化学：物源属性、构造环境和影响因素［J］. 地球与环境，2006，34（4）：75-83. Cai Guanqiang， Guo Feng， Liu Xiantai， et al. Clastic sediment geochemistry： Implications for provenance and tectonic setting and its influential factors［J］. Earth and Environment， 2006， 34（4）： 75-83. [47] Strani F， DeMiguel D， Alba D M， et al. The effects of the “0.9 Ma event” on the Mediterranean ecosystems during the Early-Middle Pleistocene transition as revealed by dental wear patterns of fossil ungulates［J］. Quaternary Science Reviews， 2019， 210： 80-89. [48] 王婷，孙有斌，刘星星. 中更新世气候转型：特征、机制和展望［J］. 科学通报，2017，62（33）：3861-3872. Wang Ting， Sun Youbin， Liu Xingxing. Mid-Pleistocene climate transition： Characteristic， mechanism and perspective［J］. Chinese Science Bulletin， 2017， 62（33）： 3861-3872. [49] 王海粟，党皓文，翦知湣. 中更新世转型时期南海北部上层水体结构演化特征：ODP1146站浮游有孔虫稳定同位素记录［J］. 第四纪研究，2019，39（2）：316-327. Wang Haisu， Dang Haowen， Jian Zhimin. Variations in the upper water structure of northern South China Sea during the Mid-Pleistocene climate transition period： Planktonic foraminifera oxygen isotope records of ODP site 1146［J］. Quaternary Sciences， 2019， 39（2）： 316-327. -
点击查看大图
图(5) / 表 (1)
计量
- 文章访问数: 1082
- HTML全文浏览量: 203
- PDF下载量: 94
- 被引次数: 0
