1. 引言
气候的频繁波动通过改变降水与地表径流的强度、沉积环境条件,显著影响沉积过程,导致沉积物类型多样,调控全球地貌格局[1] [2]。坡洪积物是坡面重力作用和暴雨径流共同作用下形成的典型沉积物,气候变化、地形条件及物质来源决定着其发育特征[3] [4],是第四纪沉积体系的重要组成部分,不仅记录着气候变化、地表水文过程等环境信息[5]-[7],还可以预防滑坡、洪水等自然灾害事件,在生态环境的维持及环境变化的科学研究中具有不可替代的作用[8],是揭示区域环境演变的关键载体。
国外对洪积物沉积学的研究可以追溯到19世纪,主要聚焦洪积物的形成、概念[9]、位置分布[10]-[12]、沉积特征及地貌识别等方面[13]。20世纪中后期,国外研究者对洪积物进行更深入研究,在地貌提取、形态研究、形成因素及对环境变化的响应等方面取得一定成果[14]-[17]。国内平原、低山丘陵大量发育洪积物[18],近年在野外识别[19] [20]、地貌提取[21] [22]、环境指示[23] [24]及气候灾害等方面具有很大进展。
相较于洪积物,坡积物、坡洪积物在地貌学上较少,一部分是在学者研究古洪水沉积序列时,作为洪积物的伴生成分被分析,其独立性研究相对有限[25]。一部分则集中在区域调查、地层分布、沉积相及沉积特征描述上[26] [27],坡洪积地貌的环境指标定量分析上研究相对薄弱。
粒度组成可以反映沉积物的搬运能力和沉积动力,是判别沉积环境和分析沉积过程的重要依据,测试简单、意义明确、对气候反映敏感,已被国内外广泛地应用到多种沉积环境的研究中[28]-[30]。辽南地区第四纪沉积物类型多样,前人对其研究多集中在区域地质调查,描述地层分布、沉积类型上,为本文剖面沉积类型的确定、沉积特征描述、沉积环境分析提供基础资料,为确定区域沉积环境变化提供坚实的基础[31]-[33]。大连石河地区位于小黑山山麓地带,保存着多层夹杂着坡洪积层。本文在前人对地层分布描述、沉积类型确定的基础上,进行沉积特征的定性描述及沉积粒度的定量分析,明确区域的沉积环境变化过程,揭示其对古气候、水动力条件、搬运距离等信息的响应,对理解辽南地区气候背景、第四纪坡洪积地貌演化及预防地质灾害具有重要意义。
2. 研究区概况
2.1. 研究区域概况
辽南以低山丘陵和河流冲积平原为主,整体地势西北高、东南低[34],石河街道位于辽宁省南部,大连市北部,行政隶属大连市金州区(39˚20'~39˚35'N, 121˚40'~122˚00'E),处于低山丘陵与沿海冲积平原的过渡带上(图1)。研究区主要岩性为洪积相的黄褐色亚砂土与砂砾石,桥头岩组二岩段的石英岩、板岩互层,冲/坡洪积相的上部亚砂土夹砾石层。研究区周围岩性为甘井子组的白云岩、白云质灰岩,南关岭组的砾屑、砂屑灰岩、石灰岩,长岭子组的板岩、钙质板岩夹大理岩等。
注:该图基于自然资源部标准底图服务网站下载的审图号为GS (2024) 0650号的标准地图制作,底图无修改。
Figure 1. Geographical location of the study area
图1. 研究区地理位置
辽南地势北高南低,河流多发源于北部山区,多为放射性水系。石河剖面(39˚19'23"N, 121˚50'25"E)位于连续起伏的低山丘陵的东侧坡麓地带,丘陵沟谷纵横,受气候的影响,夏季降水集中且多暴雨,洪涝灾害频发。剖面北部约1 km处,有一条发源于小黑山北麓名为“石河”的河流,总体自东南向西北流动,经普兰店湾入海,是金州区北部的重要河流。
2.2. 研究剖面
对石河剖面进行多次考察后,于2024年4月采集样品。剖面高5.3 m,整体呈红色褐,细粒层与砾石层交互存在。砾石大小混杂,直径约4~15 cm,大者可达1 m,磨圆度低,分选性差,基岩以石英岩为主,大理岩次之。本研究采用非等间距取样法,依据剖面沉积层序界面及粒度旋回特征,根据剖面砾石层和细粒层的堆积情况,将其自上而下分为⑧层,采集16个样品,标号为SH1至SH16,用于后续粒度分析(图2),有利于识别沉积事件的边界与频率,确保样品能准确反映不同沉积阶段的粒度特征。
其中,第①层砾石与粗细颗粒混杂,磨圆度低、分选极差;第②层含粘土质粉砂,砂含量低且分选差;第③层砾石粒径减小,粗颗粒减少;第④层砾石大小混杂、叠瓦状排列,含直径可达1 m、低磨圆度的大方石,分选极差;第⑤层砾石粒径小,混粉砂透镜体,分选极差;第⑥层砂质粉砂为主,含次棱角状较粗的砾石,分选较差;第⑦层夹薄层粗砂,分选中等;第⑧层砾石少而粗砂增多,顶部覆盖0.1 m表土层。
Figure 2. Sedimentary characteristics of the Shihe section
图2. 石河剖面沉积特征图
3. 研究方法
本文根据粒度组分的百分含量,粒度频率分布曲线、粒度三角图、粒度参数特征及其地层组合方式等方面分析研究区坡洪积物的沉积粒度特征,解译其指示的环境意义。
样品测试仪器采用美国库尔特公司生产的贝克曼LS13 320型激光衍射粒度测试仪,测试范围为0.04~2000 μm,实验在辽宁省自然地理与空间信息科学重点实验室完成。为了准保结果可靠,每样品进行3次重复实验,取平均值。实验步骤如下:
① 筛选:将样品自然风干、碾碎,过2 mm筛。
② 称重:用精度为0.1 mg的电子秤称取0.25~0.30 g样品,置于100 ml小烧杯中,每个样品称取3次。
③ 去除有机质:将配比为1:1的H2O2溶液(100 ml)加入烧杯中,静置48 h,用加热板加热,直至烧杯中不再产生气泡。
④ 去除碳酸盐:在完成第二步的烧杯中加入10 ml配比为1:3的HCl溶液,放置24 h,加热板加热至沸腾,直至无气泡产生。
⑤ 洗酸:冷却至常温之后,加入蒸馏水,放置48 h,倒掉上层清液,重复进行直至溶液呈中性。
⑥ 加分散剂:在测试前烧杯中加入10 ml浓度为0.05 mol/L的六偏磷酸钠并摇匀。
⑦ 测定:在贝克曼LS13 320型激光衍射粒度测试仪进行测试,各样品均采用超声波震荡30 s,测试3次,取平均值。
4. 结果分析
4.1. 粒度组成分析
沉积物粒度分布可以直接反映其机械组分的含量和分布,是揭示古气候演化、沉积动力特征重要指标。由于水动力搬运条件的限制,石河坡洪积物呈现粗颗粒与细颗粒、砾石混杂堆积的特点,本文对筛选(细颗粒,<2 mm)后的样品进行粒度测试。根据Udden-Wentworth的粒度分级标准[35],将沉积物粒径分为6个等级:黏土(<4 μm),粉砂(4~63 μm),极细砂(63~125 μm),细砂(125~250 μm),中砂(250~500 μm),粗砂(>500 μm)。
除砾石外,剖面整体以粉砂为主,变化范围为33.15%~66.09%之间,均值为55.30%。粘土变化范围为8.13%~35.05%,均值17.30%。极细砂变化范围为3.78%~28.30%,均值13.39%。细砂变化范围为1.82%~21.44%,均值8.00%。中砂变化范围为0.00%~13.84%,均值4.43%。粗砂变化范围0.00%~6.51%,均值1.57%。占比顺序为:为粉砂 > 粘土 > 极细砂 > 细砂 > 中砂 > 粗砂(表1)。
剖面粉砂、粘土含量整体呈现自下而上减小的趋势,极细砂、细砂含量整体呈自下而上增大的趋势,但是均具有波动性,在3.5 m处样品波动最为明显。中砂与粗砂百分比含量变化曲线相似,整体波动频繁,在剖面最顶层,中砂、粗砂含量最高。第①层粉砂含量较多、粗砂少量出现,砾石混杂。第②层黏土含量达到峰值,粗粒几乎缺失。第③层粘土含量有所降低,粉砂含量升高,砾石减少。第④层在3.5 m处极细砂与细砂含量突增,粘土含量最低,砾石巨大。第⑤层粉砂含量较高,粗砂少,存在少量砾石。第⑥层粉砂占比下降,极细砂与细砂占比增多。第⑦层粉砂含量又增加。第⑧层粗砂、中砂占比最高,粉砂含量较低,含砾石。
粒度三角图反映沉积物中砂、粉砂和粘土三者之间的比例关系及分布特征。根据三角图解法,石河剖面粉砂与砂含量较大,分别为55.30%、27.39%,粘土的含量最低,为17.30% (图3),剖面除砾石之外,粒度组成包括砂质粉砂、粘土质粉砂与粉砂质砂,且砂质粉砂含量最高,剖面自下而上整体变化趋势为粘土质粉砂–砂质粉砂–粉砂质砂,第④层沉积环境波动较大。
粒度频率曲线不仅可以直观粒度变化特征,还可以通过沉积物的分选性与峰值,分析沉积环境、水动力条件。第①层样品中的SH3与SH8、SH9的粒度频率曲线相似,但SH8第二个峰值波动比SH3大,指示他们沉积环境相似,但SH8样品比SH3样品沉积环境更为复杂,SH8处发生了更强大的搬运作用,扰乱了原本SH3上部沉积物,并在第④层沉积。其余样品呈现不规则双峰、多峰形态,SH4与SH5、SH6与SH7含量最高的峰值粒径数值一样,虽然峰值粒径占比有所差别,但可以说明SH4与SH5、SH6与SH7的沉积环境与动力、来源一致。SH10三峰分布,SH11、SH12峰值极不明显,指示水动力条件复杂且不稳定。SH13近似无峰,与SH14频率曲线为双峰,峰值粒径接近,SH15、SH16为多峰,说明顶层沉积物的沉积环境复杂。
Figure 3. Grain size triangle diagram of the Shihe slope-pluvial deposit profile
图3. 石河坡洪积物剖面粒度三角图
剖面样品粒度频率分布曲线多呈双峰、多峰的形态(图4),说明此剖面沉积环境复杂,应为坡面重力、片流作用与暴雨径流事件的叠加,水动力充足,波动性强。
Figure 4. Grain size frequency distribution of slope-pluvial deposits in the Shihe River
图4. 石河坡洪积物粒度频率分布图
4.2. 粒度参数特征分析
粒度参数是定量表述、分析沉积物粒度特征的常用方法,通常利用平均粒径(Mz)、分选系数(σ)、偏度(SK)、峰态(KG)这四个参数反应颗粒物质粒径的分布特征。本文采用Folk和Wald公式进行计算[36]。
平均粒径(Mz)指示样品所有粒径的平均值,反映整体颗粒的“粗细”程度,Mz越大,颗粒越粗;反之越细。计算公式为:
(1)
分选系数(σ)表示粒径分布的离散程度,σ越小,说明颗粒大小越均匀;反之,颗粒混杂,分选性差。计算公式为:
(2)
式中,Ф5为累积频率曲线5%所对应的粒径值Φ、Ф16、Ф25、Ф50、Ф84、Ф95为16%、25%、50%、84%和95%处所对应的粒径值Φ。
偏度(SK)指示粒径分布对称性的度量,分为正偏(数值0.1~0.3)指示细颗粒后期充填、负偏(数值−0.3~−0.1),反映粗颗粒优势(如洪水沉积)、极负偏(数值−1~−0.3)、近对称(数值−0.1~0.1)、极正偏(数值0.3~1)。计算公式如下:
(3)
峰态(KG)表示粒径分布的集中趋势,高峰度(窄峰)表示粒径集中,低峰度(宽峰)表示粒径分散。计算公式为:
(4)
研究区剖面平均粒径(Mz)最小值与最大值分别为4.14 Φ和6.98 Φ,分别位于第⑧层3、第③层,均值为5.50 Φ,属于粉砂。剖面Mz数值具有波动性,呈现出大幅度上升–大幅度下降–上升–下降的趋势,说明整体粒径经历了变粗–变细–粗–细的阶段。分选系数(σ)变化范围为3.81~7.02,均值为5.36,数值较大,属“极差”等级,表明粗细颗粒混杂,该结论与野外观察的剖面粗细颗粒混杂、砾石磨圆度低且多呈次棱角状、分选性差的结论高度一致,共同指示沉积物质未经历长时间搬运、磨圆作用,呈现典型近源快速堆积特征。
剖面偏度(SK)的变化范围为−0.03~0.38,均值0.11,共有近对称、正偏、极正偏三个等级,占比分别为56.2%、37.5%、6.3%。说明剖面搬运动力以坡面片流为主,在沟谷汇流、暴雨冲刷的作用下,粗颗堆积,洪水退却后悬浮物沉降,使剖面中含有大量细颗粒物质,既有坡积的影响,也有洪水后期细颗粒填充的影响。峰态值(KG)的变化范围为0.86~1.18,均值为0.96,第①~⑦层大部分属于中等峰度,表示粒径不集中,第⑧层属于低峰度,粒径分散。峰态(KG)分为三个等级:中等、宽、尖宽,占比依次为中等(81.3%) > 宽(12.5%) > 尖宽(6.2%) (表1)。
5. 讨论
石河剖面记录的沉积环境变化是构造背景与气候演变共同作用的结果。辽南地区新近纪以来经历了显著的构造抬升,石河地区处于低山丘陵与沿海冲积平原的过渡带上低山丘陵向沿海冲积平原的过渡带上,这种地形为坡洪积物的发育提供了物源条件和动力基础。且辽南地区于新近纪发育的红色风化壳,是高温高湿气候条件下的风化产物,其上覆黄土则说明该区域经历了从热带或亚热带湿热气候向干寒气候的转变[37] http://wfo.lib.ustc.edu.cn/D/Periodical_yndlhjyj201802009.aspx。这种气候转型过程对区域风化作用和沉积物特征产生了深远影响。
粒度参数进一步揭示了沉积环境对气候变化的响应。剖面平均粒径的垂向波动以及分选系数极差,不仅指示了近源快速堆积的特点,也反映水动力条件的频繁变化。偏度以正偏和近对称为主,说明沉积过程中既有洪水粗颗粒沉积,也有后期细颗粒填充,这种特点与辽东半岛全新世气候波动导致的沉积动力变化相一致。频率曲线多峰形态则进一步证实了沉积动力的复杂性。剖面第④层观察到粗颗粒富集,含大砾石,黏土含量最低,这反映了环境湿润、洪水突发、水动力强劲的沉积环境。这一层次很可能对应古气候极端降水事件,强大的水动力条件导致了粗粒物质的快速搬运和堆积。因此,大连石河地区坡洪积物记录的环境突变事件主要受构造–气候耦合作用控制。
Table 1. Grain size composition of the Shihe profile
表1. 石河剖面粒度组成
层位 |
样品编号 |
深度 /m |
粘土 /% |
粉砂 /% |
极细砂 /% |
细砂 /% |
中砂 /% |
粗砂 /% |
平均粒径 |
分选系数 |
偏度 |
峰态 |
⑧ |
SH16 |
0.4 |
9.34 |
38.59 |
17.74 |
14.14 |
13.68 |
6.51 |
4.14 |
3.85 |
0.20 |
0.94 |
SH15 |
0.9 |
13.91 |
42.79 |
13.96 |
12.02 |
12.75 |
4.57 |
4.63 |
4.39 |
0.16 |
0.86 |
⑦ |
SH14 |
1.1 |
10.31 |
53.65 |
15.63 |
10.91 |
6.50 |
3.00 |
4.87 |
4.80 |
0.06 |
0.95 |
SH13 |
1.3 |
16.52 |
61.78 |
11.64 |
6.02 |
3.06 |
0.98 |
5.71 |
5.68 |
0.04 |
1.00 |
⑥ |
SH12 |
2 |
16.71 |
56.83 |
11.55 |
7.83 |
5.44 |
1.63 |
5.51 |
5.58 |
−0.02 |
0.97 |
SH11 |
2.8 |
11.92 |
48.76 |
17.94 |
13.27 |
7.10 |
1.00 |
4.85 |
4.60 |
0.19 |
0.92 |
⑤ |
SH10 |
3 |
19.62 |
63.87 |
10.21 |
4.19 |
1.74 |
0.37 |
6.01 |
5.93 |
0.06 |
0.94 |
SH9 |
3.2 |
12.51 |
53.33 |
19.39 |
11.16 |
3.55 |
0.06 |
5.08 |
4.76 |
0.23 |
0.93 |
④ |
SH8 |
3.5 |
8.13 |
38.05 |
28.31 |
19.21 |
4.65 |
1.65 |
4.25 |
3.81 |
0.38 |
1.18 |
SH7 |
3.7 |
17.77 |
62.39 |
13.65 |
6.03 |
0.16 |
0.00 |
5.83 |
5.75 |
0.09 |
0.91 |
③ |
SH6 |
4 |
12.28 |
65.13 |
14.06 |
4.68 |
2.05 |
1.80 |
5.49 |
5.37 |
0.10 |
1.04 |
SH5 |
4.3 |
31.54 |
62.86 |
3.78 |
1.82 |
0.00 |
0.00 |
6.98 |
7.02 |
−0.02 |
1.00 |
② |
SH4 |
4.5 |
35.05 |
57.36 |
5.31 |
2.28 |
0.00 |
0.00 |
6.97 |
7.02 |
−0.03 |
0.87 |
SH3 |
4.7 |
18.75 |
66.09 |
9.38 |
3.30 |
2.11 |
0.38 |
5.96 |
5.70 |
0.16 |
1.00 |
① |
SH2 |
4.9 |
25.17 |
58.53 |
8.58 |
4.39 |
2.70 |
0.63 |
6.29 |
6.20 |
0.03 |
0.96 |
SH1 |
5.3 |
17.29 |
54.87 |
13.17 |
6.72 |
5.34 |
2.61 |
5.44 |
5.27 |
0.08 |
1.01 |
6. 结论
本文选取大连石河地区坡洪积物沉积剖面,对沉积物进行粒度实验,分析不同地层沉积物的环境指标特征,探讨指示的环境意义,得出以下结论:
(1) 除砾石外,剖面整体以粉砂为主,为砂质粉砂沉积。剖面第①~③层粒径含量无明显上升下降趋势,但粘土含量略高,指示水动力从强到弱,气候略干。第④层粒径变化异常,粗颗粒富集,黏土含量最低,伴随巨型砾石,指示当时沉积气候湿润,水动力强劲。第⑤~⑦层整体粒径变细,但粗砂含量升高,指示其为稳定低能环境、间歇性强水流环境下的产物。第⑧层粗砂与中砂含量最高,粉砂较低,含砾石,说明其在暴雨期,粗粒快速堆积而成。
(2) 平均粒径介于4.14~6.98 Φ之间,垂向波动显著,粒径的粗细旋回,沉积动力多变。分选系数范围为3.81~7.02,均值5.36,属“极差”等级,表明颗粒混杂、物源复杂,与野外观察的粗粒颗粒混杂、砾石磨圆度低且多呈次棱角状、分选差特征高度一致,指示物质未经历长时间搬运磨圆作用,为近源快速堆积。偏度范围处于−0.03~0.38之间,均值0.11,以近对称和正偏为主,指示既有坡积影响,又有洪水后期细颗粒填充作用。峰态介于0.86~1.18之间,均值0.96,81.3%样品属中等峰态,只有第⑧层为低峰态,粗颗粒富集、分选性低,共同指示极端强降水事件。
(3) 石河坡洪积物粒度波动指示气候变化控制沉积动力,粒度频率曲线呈双峰、多峰形态,表明石河剖面坡洪积物水动力充足且沉积条件复杂,是坡面重力、片流作用与突发性洪水沉积事件的叠加。
基金项目
辽宁省社会科学规划基金项目,项目编号:L20AKG001。
NOTES
*通讯作者。