粒度指标指示的坝上御道口地区8 000 a B.P.环境变化

黄劭康，贾红娟，武俊叶，张 昊，殷志强

(1.河北地质大学资源学院，河北 石家庄 050031；2.河北地质大学实验实践教学中心，河北 石家庄 050031；3.中国地质环境监测院，北京 100081)

河北坝上地区承担着“为京津涵水源、阻沙源、提供生态支撑”的重大历史使命，是京津冀地区抵御浑善达克沙地南侵的最后一道防线。坝上地区位于季风区与非季风区交汇带，对气候变化非常敏感[1]，可较好地揭示未来可能增温背景下环境和生态系统对气候变化的响应。近几十年坝上地区植被退化、土地沙化、盐碱化面积不断增大，草原生态环境脆弱，自然湿地萎缩干化为碱滩、沼泽，永久湿地变为季节性湿地趋势明显，水源涵养保障能力弱化，旅游经济受到影响。坝上地区以人工林为代表的人工植被建设正面临一个严峻问题——土壤干旱化，认识本区气候变化对未来增温背景下采取应对措施有重要意义。

前人对该区环境有过一些研究，研究结果存在较多分歧，最主要的分歧是全新世湿润期的结束时间及温湿匹配。御道口牧场剖面研究显示此区6 000 a来经历7次明显变化[2]，但深150 cm的剖面只有2个年龄控制，7个阶段年龄划分有欠缺。孙建中等[3]研究显示8 200～5 700 a B.P.为湿润气候，5 700～4 600 a B.P.为干冷气候，4 600～4 000 a B.P.为温湿气候，4 000～3 500 a B.P.为冷干环境，3 500～3 000 a B.P.气候暖湿, 3 000 a B.P.以来气候干旱。齐惠慧等[4]对安固里淖的研究结果表明：8 000～3 200 a B.P.气候偏湿，3 200～0 a B.P.气候干旱。邱维理等[5]研究结果表明6 230～5 300 a B.P.是安固里淖湖面快速下降期，5 300 a B.P.以来湖面波动性下降，此剖面6000年来只有2个年龄。翟秋敏等[6]对安固里淖的研究结果显示：7 453～4 617 a B.P.以湿润气候为主，但是存在相对干旱的时段(5 137～4 617 a B.P.)；4 617～574 a B.P.以干冷气候为主；574 a B.P.以来，前期气候向湿润方向发展，但是在383 a B.P.以后又转为干旱。另外同一剖面不同代用指标结果相互矛盾，阳小兰等[7]研究结果显示安固里淖5 000～1 070 a B.P.是气候相对湿润阶段,中世纪温暖期(1 070～620 a B.P.)气候干旱化明显,小冰期(620～45 a B.P.)干旱化及寒冷程度进一步加剧。刘林敬等[8]对同一剖面的研究结果则显示：在5 030～3 070 a B.P.时期气候相对湿润；在3 070～1 120 a B.P.期间温暖指数逐渐升高，降水减少；1 120～800 a B.P.期间，较高的温暖指数指示气候相对温暖，降水量增加；800～270 a B.P.期间，温暖指数的显著降低，降水减少。

综观这些研究，发现造成这些分歧的主要原因是剖面年龄数据少、分辨率低，不同研究者利用相同指标或不同指标数据解释存在差异，因此在未来可能增温的背景下，环境变化需进一步研究。本研究采用较高分辨率年代数据重建了御道口地区8 000 a B.P.以来的环境变迁。

1 研究区概况

图1 YDK剖面位置Fig.1 Location of YDK section

御道口剖面(YDK剖面：北纬42°07′14.04″东经117°00′47.94″，见图1)位于河北承德市围场县西北部的坝上地区。此区属大陆性季风气候，据塞罕坝气象站多年气象数据，1月平均气温 -21.76 ℃，7月平均气温16.28 ℃；年降水量在 258～652 mm之间。该地区地貌上属内蒙古高原的东南边缘与冀西北山地的交汇带，在地质构造上为内蒙高原背斜面，属大兴安岭余脉[9]。

2019年5月，在水平方向选择第四系沉积物，垂直方向以 2.5 cm 间距，共采样112个。剖面深度及岩性见图2。

图2 YDK剖面岩性柱和年代-深度曲线Fig.2 Lithology, age-depth curve of YDK section

2 研究方法

2.1 粒度分析测试

根据样品粗细不同称取0.1～1.5 g放入500 mL干净烧杯中，加入10 mL 30%的H2O2摇匀。之后静置1～2 h使其充分反应，去除有机质。另加入10 mL 10%的HCl去除样品中碳酸盐，静置12 h后，抽去上层清液备用。加入5～10 mL1 mol/L NaPO3超声振荡10 min,用英国Malvern公司Mastersizer3000激光粒度仪分析。为保证实验精确性，将每个样品测量时的遮光度都控制在2%～20%之间，并多次进行测量。粒度组分用Grain Analysis 2013-01软件进行提取[10]。

2.2 年龄分析

剖面年龄标尺由5个AMS14C年龄控制，用全岩沉积物有机质进行AMS14C年代测定。测试工作由美国Beta实验室完成，结果见表1，并对获得的年代通过BetaCal3.21:HPD method:INTCAL13进行日历年校正。

表1 YDK剖面年龄测试结果

剖面年龄框架通过采用回归插值法建立(方程为y=177.68e0.014 1x，R2=0.987 3)。校正后YDK剖面年龄跨度范围184～8 000 cal.a B.P。

3 结果与分析

3.1 粒度分布曲线

御道口剖面粒度分布曲线主要有a、b两种类型(图3)。每种分布图由不同的粒度组分组成，C1粒径<2 μm、C2粒径2～15 μm、C3粒径15～70 μm、C4粒径70～700 μm。a类型有C1、C2、C3、C4四种不同的组分；b类型有C2、C3、C4三种不同的组分。

图3 YDK剖面主要粒度曲线类型及组分Fig.3 The main types and five compoments of grain-size distribution of the YDK section

3.2 各粒级标准偏差

在0.01～3 500.00 μm量程的基础上，利用对数划分方法，分出100个粒级，并计算出每个粒级含量在剖面上的标准偏差，绘制出粒级-标准偏差曲线图(图4)。就整体而言，御道口剖面标准偏差曲线图中有三个峰值，分别在0.2 μm、8 μm、309 μm左右，最突出的峰值是309 μm，表明此组分可较好地指示剖面环境变化。

图4 YDK剖面粒级-标准偏差曲线图Fig.4 Grain size-standard deviation curve of YDK section

3.3 粒度对环境的指示意义及对环境变化敏感的粒径范围

粒度是沉积物机械组成的一个重要参数，粒度因其明确的环境指示意义，成为恢复古气候、古环境的重要代用指标之一[11]。不同粒级沉积物和沉降受不同的搬运营力控制，因此实测的沉积物粒度数据可按照粒径区间进行划分，反过来作为沉积动力环境研究的重要依据。

一般情况下，较强的动力条件能搬运较大颗粒，反之只能搬运细颗粒。风成粉尘的大气动力学理论计算和实际观测都表明，几个微米以下的细粒组分可上升到几千米以内的任一高度作长距离的悬移，高空长距离悬浮是它的主要搬运方式。10～70 μm的粉砂组分则主要以短距离悬移方式搬运，砂和粉砂级粗粒组分(>70 μm)每次起动只能上升到近地表的几厘米到几米的高度，并在水平方向上跃移同样量级的距离，形成风成沙[12]。

剖面沉积物粒径分布共有4个组分，分别为C1、C2、C3和C4，其中值粒径范围，C1粒径<2 μm，C2粒径2～15 μm、C3粒径15～70 μm、C4粒径70～700 μm。组分C1、C2代表了在流体介质中长期悬浮搬运组分和中粒悬浮组分，与湍流强度有关；C3是粉尘中重力沉降主导的粗粒悬浮组分；C4属粗粒跳跃组分[13]。C1、C2、C3含量高地面植被好，能捕捉较多的悬浮颗粒，因此反映气候湿润；C4含量高说明悬浮颗粒少、跳跃颗粒多，反映地面植被差，指示气候干旱。

剖面所有样品同一粒级百分含量的标准偏差可指示该粒级含量的波动强度，标准偏差越大，表明该粒级颗粒含量在整个剖面中的变化幅度越大，可作为指示沉积特征变化的敏感组分[14-15]。整体而言，剖面标准偏差曲线图中有三个峰值，分别在0.2，8，309 μm左右，最突出的峰值是309 μm，表明此组分可较好地指示剖面环境变化。剖面粗颗粒含量高，指示风力较大气候干旱，反之细颗粒含量高，指示风力较小气候湿润。

4 御道口地区8 000 a B.P.环境变化与区域对比

根据中值粒径(Md)、粒级-标准偏差、粒度分布、各粒级百分含量变化可将剖面划分为三个阶段(图5)：

图5 御道口YDK剖面沉积物粒度指标特征随年代的变化Fig.5 Changes of the characteristics of sediment grain index with age in YDK section

第1阶段(8 000～6 100 aB.P.)气候干旱：此时段粒度分布曲线主要为a型，Md值较大，波动较小，<0.87 μm百分含量较少，0.87～45.6 μm和>45.6 μm百分含量较多，各粒级百分含量变化不大。此时段C1、 C3含量少，C2、C4含量高。此时段动力较大，粉尘距源区较近，风沙活动较强，携带较多粗颗粒物质造成粒径较粗，气候干燥，区域植被较差。

此时段区域研究结果存在较大分歧，孙建中等通过研究坝上地区赋存的黄土与古土壤，古沙丘与古土壤和湖泊沉积物中的环境变化信息得出8 200～5 700 a B.P.气候冷湿润[3]。安固里淖研究结果显示8 000～6 000 a B.P.以森林草原植被为主，气候湿润[4]。呼伦湖重建降水量较高(图6e)[16]。此阶段有些记录结果相反，达里湖δ13C呈增加趋势,指示湖区大气湿度逐渐降低(图6g)[17]。三江平原风沙含量8 000～6 200 a B.P.较高，指示相对干旱[18](图6h)。翁牛特磁化率值低，呈上升趋势，指示降水量少[19](图6i)。

图6 YDK剖面中值粒径和其他记录对比Fig.6 Comparison of the Md of this study with other records

第2阶段(6 100～1 700 a B.P.)为湿期：此阶段粒度分布线类型为a型，与第Ⅰ阶段相比，Md减小，为剖面最低阶段，且减小的幅度较大，从90～100 μm减小到20 μm左右；<0.87 μm百分含量增加，从2%增长到7%左右;0.87～45.6 μm百分含量增加，为剖面最大阶段；>45.6 μm百分含量减小，含量减少了约10%，为剖面最小阶段。此时段C1、C2、C3含量高，C4含量低。此时段动力较小，粉尘距源区较远，粗颗粒物质减少，估计此时段气候湿润，区域植被较好，风沙活动减弱。此时段以3 700 a B.P为界可划分为两个亚段，2a比2b小于0.87 μm、C1、C3含量高，Md较小，应该更湿润。

此时段刘永慧等[2]研究显示御道口地区5 692～5 356 a B.P.气候凉略湿，5 356～4 684 a B.P.气候特征为凉干，4 684～947 a B.P.气候变得相对略湿或湿。阳小兰等[7]研究结果显示安固里淖5 000～1 070 a B.P.是气候相对湿润阶段(图6d)。孙建中等[3]研究结果显示此时段气候有多次波动，5 700～4 600 a B.P.为干冷气候，4 600～4 000 a B.P.为温湿气候，4 000～3 500 a B.P.为冷干环境，3 500～3 000 a B.P.气候暖湿，3 000 a B.P以来干旱。辽西地区 5 000 a B.P.以来环境演变显示此时段乔木含量较高相对湿润，期间存在气候波动，5 000～4 400 a B.P.环境为冷干，4 400～2 750 a B.P环境以温暖湿润为主，2 750～1 200 a B.P.环境为冷干[20](图6f)。内蒙古呼伦湖定量重建降水量6 100～1 700 a B.P.大部分时间较高，4 400～3 350 a B.P.，湖区耐旱的藜科植物大量生长，荒漠化十分严重，降水和温度都明显降低(图6e)[16]。内蒙古达里湖δ13C指示湖区大气湿度较高，自高位波动下降[17](图6g)。三江平原风沙含量 6 200～1 700 a B.P.较低，指示此时段相对湿润[20](图6h)。翁牛特6 000～3 900 a B.P.磁化率值高，指示植被覆盖好，气候相对湿润，3 900 a B.P.之后磁化率降低[19](图6i)。内蒙古辉腾锡勒5 390～3 300 a B.P.气候偏湿，3 300～2 380 a B.P气候偏干，2 380～1 160 a B.P.气候凉湿[21]。内蒙古巴汗淖湖6 100～3 400 a B.P.气候温暖湿润，3 400～2 700 a B.P.气候向干旱化方向发展,2 700 a B.P.以来气候进一步变干[22]。

尽管有些研究年代相差较远，但较高分别率研究显示此时段大部分时间相对湿润。

第3阶段(1 700 a B.P.至今)为干期：Md为剖面最大阶段，Md及各粒级百分含量波动明显，可以分为两个子阶段：

3a阶段：(1 700～900 a B.P.)：此时段粒度分布曲线主要为b型，1 700 a B.P. Md迅速增大，由 20 μm 增长到356 μm左右，<0.87 μm百分含量为剖面最低，0.87～45.6 μm由10%下降到3%左右，>45.6 μm百分含量在90%～96%之间，波动幅度不大。此时段C1、C2、C3含量减少，C4含量增加。

3b阶段(900 a B.P.至今)：此时段粒度分布曲线主要为a型，Md与上一阶段比明显下降，总体在130 μm～245 μm之间波动。<0.87 μm百分含量仍较低，大部分都在2%以内，0.87～45.6 μm和 >45.6 μm 两个粒级百分含量变化呈负相关，含量比上一阶段增加。此时段C1含量少，C2、C3含量增加，C4高位波动。1 700 a B.P.至今动力较大，粉尘距源区较近，估计此时段气候变干，区域植被较差，风沙活动强。

从区域研究结果看虽然年代略有差异但总体来说均显示此时段干旱。刘永慧等[2]御道口牧场的研究结果显示947 a B.P.至今的气候总体为凉略干。阳小兰等[7]认为安固里淖1 000 a B.P.以来，干旱程度加剧(图6d)。辽西地区1 200 a B.P.～现在环境干旱[20](图6f)。内蒙古辉腾锡勒1160 a B.P.以来气候向温凉偏干转变[21]。呼伦湖最近1 700年里降水先减少后增加，孢粉谱中伴人植物花粉藜科和禾本科同时大量出现，而最近的500年里湖区气候逐渐向暖湿的方向发展[16](图6e)。藜科孢粉指示干旱，其增长说明气候变干，御道口剖面900年以来中值粒径比1 700～900 a B.P.小，说明这一阶段是干旱时段里湿度略增阶段。达里湖δ13C指示湖区大气湿度呈先低后高趋势(图6g)和YDK剖面中值粒径变化趋势一致。三江平原风沙含量先减少后呈增加趋势指示气候逐渐变干。

5 结论

本文利用5个 AMS14C年龄建立御道口剖面年代框架，根据112个样品粒度数据，重建了河北坝上御道口地区8 000 a B.P.以来的古环境古气候，并把这 8 000 年分成3个阶段：

第一阶段：8 000～6 100 a B.P.，此时段中值粒径较大，细颗粒物百分含量低，粗颗粒百分含量高，各粒级百分含量变化不大。粒度指标指示此时段动力较大，粉尘距源区较近，风沙活动较强，携带较多粗颗粒物质造成粒径较粗，估计此时段气候干燥，区域植被较差。

第二阶段：6 100～1 700 a B.P.,此时段中值粒径值为剖面最小阶段，悬浮组分含量高，粗颗粒物质减少，指示动力较小，粉尘距源区较远，估计此时段气候湿润，区域植被较好，风沙活动减弱。区域对比结果看此时段大部分时间湿润。

第三阶段：1 700 a B.P.至今，中值粒径为剖面最大阶段，C4含量高，悬浮颗粒少、跳跃颗粒多，动力较大，粉尘距源区较近，估计此时段气候变干，区域植被较差，风沙活动强。从区域研究结果看虽然年代略有差异但总体来说近年来研究结果显示气候趋向于干旱。