GRAIN SIZE DISTRIBUTION IN THE TIDAL ZONE AND SHALLOW SEA ON THE NORTH SIDE OF THE HUANGHUA HARBOR
-
摘要: 在黄骅港北侧潮间带和浅海区,分别获得1个柱状样品(岩心长度132 cm,共分得44个样品)和50个底质样品,通过激光粒度分析仪对沉积物粒度进行了测试,获得了研究区底质和柱状样品的沉积物粒度特征,结果显示:研究区底质样品只有3种类型,即黏土质粉砂、粉砂和砂质粉砂;砂质粉砂分布在港池防沙堤附近的潮下带地区,粉砂主要分布在大辛堡村东水深-2~-5 m的浅海区及站位13,其余区域均为黏土质粉砂;黄骅港防沙堤截断了沿岸流,并改变了其方向,使其携带的沉积物在潮间带和防沙堤附近沉积,形成了黏土质粉砂在潮间带和港池防沙堤附近的沉积;与其相邻的砂质粉砂区是废弃三角洲海岸自然演化的证据;表层沉积物的粒度特征表面,黄骅港防沙堤的影响仅限于潮间带和其附近海区;根据沉积物粒度参数特征,HHG2孔可分为4个演化阶段,从下至上分别为:(1)132~100.5 cm的沉积物粗细小幅波动段;(2)100.5~56 cm的沉积物缓慢粗化段;(3)56~16 cm的事件影响型沉积物粗细快速波动段;(4)16~0 cm的沉积物快速粗化段。结合137Cs示踪定年和黄河在大口河附近入海历史,可知16~0 cm的沉积物快速粗化段,是对1917年以来黄河入海位置逐渐东移引起的河流作用减弱以及海平面上升导致的海洋作用增强的响应。Abstract: 1 core (132 cm long and 44 samples) and 50 bottom sediments samples were collected in the offshore area on the north side of the Huanghua Harbor, and grain size analysis was carried out with Mastersizer2000. Results indicate that there are three types of bottom sediments, namely sandy silt, silt and clayey silt, in the area. Sandy silt is distributed in the subtidal zone on the north side of the harbor, silt mainly in the offshore area -2m~-5m in water depth at Daxinpu and sampling site 13, and clayey silt in the other area. The Harbor has cut off the coastal current, and changed it to the direction along the dam of the Harbor, that caused the currents-carried sediments deposited on the tidal flat and the area near the dam. As an evidence of the abandoned deltaic coast, sandy silt near the Harbor suggests that the impact of the construction of the Harbor is limited to the tidal flat and its nearshore area. According to the vertical grain size distribution pattern of the sediments, the Hole HHG2 can be divided into four segments from bottom to top:(1) the segment with slight fluctuation of coarse-fine grains in 132~100.5 cm, (2) the segment slowly coarsening up in 100.5~56 cm, (3) the segment with rapidly coarse-fine fluctuation caused by episodic events in 56~16 cm, and (4) the segment rapidly coarsening up in 16~0 cm. Furthermore, the forth segment of 16~0 cm is well correlated to the period from 1917 up to present defined by 137Cs tracing and dating. During the period, the fluvial process was weakened due to the shift of the Yellow River eastward, and the enhancement of marine process by sea level rise.
-
Keywords:
- grain size /
- bottom sediment /
- core /
- Huanghua Harbor
-
自然情况下流域的水沙通量主要受气候变化的影响[1]；但近100年来，流域人类活动加剧(主要包括水库建设、土地利用变化、矿山开采、水土保持等)导致的流域变化已经逐渐成为影响流域入海输沙量变化的最主要因素[2, 3]。统计表明，在全球受人类活动影响较大的142条代表性河流中，有近50%的河流输沙量呈现下降趋势；其中，中国最主要的10条河流近十年来的年均入海输沙量比1950—2006年的多年平均值低50%[4]，并因此对河口水下三角洲及邻近海域的沉积环境和生态系统产生了巨大影响[5]。研究表明，流域产沙受到流域的地貌、地质特征、气候变化和人类活动的共同影响[6]。但目前针对水沙通量变化的研究，多为定性地综合因素分析，缺少不同控制因子贡献的定量估算[2, 3]。流域水文模型可从正演的角度模拟气候变化和人类活动对于流域水沙的影响，从而确定不同因子的贡献及其对流域水文过程的影响，是解决上述问题的有效手段[7, 8]。
鸭绿江作为中国东部陆架区北部一条重要的河流，其入海水沙通量变化与辽东半岛东部全新世泥质沉积体系发育规模和形成演化有密切关系[9]。鸭绿江流域在有水沙观测数据之前的1940年就已经修建了大型水库(水丰水库，库容约占鸭绿江年均径流量的56%)。因此，仅依据现有的入海水沙通量数据(1956年开始)无法解释清水河流的鸭绿江为何形成高悬沙浓度的河口，且较低的入海输沙量何以形成不断淤涨的宽广的潮间带浅滩，乃至体积巨大的辽东半岛东部泥质区[10]。测年数据显示，鸭绿江西岸潮滩的沉积速率在水丰水库修建后减少了90%[11]。这也间接证明了，鸭绿江流域有水文记录以来的入海输沙量要远小于水库修建前。因此弄清楚大坝修建以前鸭绿江流域的入海水沙通量对于研究辽东半岛东部泥质沉积体系的源汇关系以及气候变化和人类活动影响下的源汇过程具有重要的指导意义。
全世界绝大多数河流的可靠水文泥沙观测记录始于1950年，导致目前有关水沙通量变化及其影响机制的研究仅限于较短时间尺度[12]，而长周期、可靠的水沙数据是研究河口-陆架区沉积环境变化和地貌演化的基础。模型可模拟长时间尺度的水沙通量，还可通过模拟不同情景模式下流域的水沙通量变化，估算不同影响因子的贡献量，进而分析其对流域水文过程的影响。HydroTrend3.0.4是以气候为驱动因素的概念性集总式流域水文模型，通过输入不同的影响水沙通量变化的因子，模拟不同驱动模式下的流域水沙通量[7, 8]，另外，该模型还可模拟高纬度流域冰雪融水对水沙通量的贡献，因此在全世界许多不同类型河流入海水沙通量变化的模拟中已得到了成功的应用[13-15]。本文拟应用HydroTrend3.0.4水文模型模拟鸭绿江流域入海的水沙通量，将实测水沙数据作为验证资料，通过改变模型的输入参数，模拟出气候变化、人类活动以及两者共同驱动下的流域水沙通量变化，从而估算出主要人类活动(植被覆盖变化、大坝建设)对流域入海水沙通量的贡献量。
1. 鸭绿江流域概况
鸭绿江是一条沿断裂发育而成的河流(图 1)，发源于海拔高度2500m的长白山，沿中朝边境流向西南，于辽宁省丹东市汇入北黄海。河流全长790km，流域面积6.45×104km2，有浑江、虚川江、秃噜江、瑷河等多条支流[16]。鸭绿江具有山溪性河流的特征，除上游小片区域为玄武岩外，整个河床由前寒武纪岩石构成，流域内较大的河床比降易产生暴涨、暴落的地表径流，多年(1958—2012)平均径流量为259.2×108m3/a，年入海输沙量为1.48Mt/a[10]。鸭绿江流域属于典型的温带季风气候，多年平均气温为5.2~7.6℃，每年的12月初至次年4月份为江面冰封期，多年平均降雨量约800~1200mm[16]。荒沟站是鸭绿江下游的控口水文站，集水面积为55420km2；梨树沟站是瑷河的控口水文站，集水面积为5629km2。两个水文站的集水面积占整个鸭绿江流域面积的95%，因此可以将两者的水沙通量作为整个鸭绿江流域的入海水沙通量。鸭绿江流程短流域面积小，因而对于气候变化和人类活动的响应更加迅速[17]，更加适合于模拟气候变化和人类活动共同影响下流域的水沙过程。
2. 材料与方法
2.1 数据来源
2.1.1 气象和水文数据
荒沟和梨树沟1958—2012年逐月水沙通量数据来源于中国水利部。鸭绿江流域上、中和下游的长白、集安和丹东站1958—2012年的月均降雨和温度数据来源于中国气象数据共享网(http://data.cma.cn/)(图 2)。
![]()
图 2 长白、集安和丹东1958—2012年月均降雨量、气温数据Figure 2. 1958-2012 monthly precipitation and temperature data at the Chang Bai, Ji An and Dan Dong stations2.1.2 植被覆盖度和人口密度
林地属于不易发生水土流失的区域，而被灌木、草地以及农作物/经济果木覆盖的区域属于易侵蚀区域[18]，因此本文依据鸭绿江流域林地覆盖面积的变化来估算水土流失状况变化。鸭绿江是中国与朝鲜的界河，朝鲜境内的流域面积占整个流域面积的45.6%。而流域人类活动的形式和强度与国家政策有着紧密的联系，所以不能够简单的使用我国一侧的流域林地覆盖率来代表整个流域的林地覆盖状况。本文获取鸭绿江流域(1985、1988、1994、1998、2003、2006、2009、2012)的遥感影像，使用像元二分模型获取流域1985—2012年的林地覆盖率[19](图 3)。依据1985—2000年间东北三省的林地覆盖度变化(中国经济大数据研究平台http://data.cnki.net/Home/Index/)和鸭绿江流域林地覆盖度变化的关系(图 4a)，虽然，2000—2012年间，我国一侧在鸭绿江流域继续实施大规模的水土保持，流域的林地覆盖度不断增加，但朝鲜一侧鸭绿江流域由于大饥荒却进行了大规模的土地开垦，流域的林地覆盖度降低。导致了2000—2012年间东北三省的林地覆盖度变化和鸭绿江流域林地覆盖度变化不具有变化一致性(图 4a，4b)。但是，1985—2000年间，东北三省的林地覆盖度变化和鸭绿江流域林地覆盖度变化，具有空间差异和时间变化一致性(图 4a，4b)。因此，根据1985—2000年间东北三省的林地覆盖度和鸭绿江流域林地覆盖度的关系，对1940—1985年间的鸭绿江流域林地覆盖度进行估算。最终获取鸭绿江流域1940—2012年的林地覆盖率(图 4b)。
![]()
图 4 东北三省林地覆盖率与鸭绿江流域林地覆盖率变化关系(a)，1940—2012年东北三省、鸭绿江流域和瑷河流域林地覆盖率变化(b)、1950—2012年鸭绿江流域和瑷河流域人口变化(c)和流域累积库容指数变化(d)(水库库容指数指水库库容量与流域年均入海径流量的比值)Figure 4. The relationship between vegetation coverage in northeast three provinces and Yalu catchment during 1985-2000(a), vegetation change of northeast three provinces, Yalu catchment and Ai catchment during 1940-2012 (b), the variations in population density during 1950-2012 (c), and the total reservoir storage capacity index (RSCI) of Yalu rivers during 1940-2012 (d).Here, the reservoir storage capacity index is defined as the ratio of the reservoir storage capacity to average annual water discharge of the contributed catchment.鸭绿江流域中国一侧1949—2012年人口数据来自于中华人民共和国统计年鉴。由于研究区朝鲜的人口数据无法直接获取，故根据世界银行获取朝鲜1961—2012年以来的全国人口数据，通过(DIVA-GIS)获取朝鲜全国人口分布图，并获取鸭绿江朝鲜一侧人口占朝鲜全国人口的比例，因此，计算出1961—2012年朝鲜鸭绿江流域人口数据(图 4c)。
2.1.3 流域高程与大坝
流域的高程数据来自中国科学院地理数据云(http://www.gscloud.cn/)的GDEM30m分辨率数据。1940年鸭绿江流域下游修建了水丰水库，其库容量占整个流域年均径流量的56%。1965年以来，鸭绿江流域又陆续修建了云峰、渭源、太平湾、临江、桓仁和回龙山水库，从上游到下游形成了梯级水库(图 4d)。
2.2 模型概述
Hydro Trend3.0.4集总式概念水文模型是以气候作为驱动因素，依据水量平衡的基本原理，综合考虑流域地形地质特征(高程、河网、坡降和岩性)、物理参数(温度、降雨、蒸发蒸腾、冰雪覆盖)和人类活动影响(土壤侵蚀和大坝)模拟出流域出口的水沙通量[7, 8]。该模型最初由Nicholson于1992年开发，最新版本3.0.4由Kettner于2010年改进。HydroTrend3.0.4水文模型在相关文献中已有大量的介绍和应用实例[7, 8]，本研究中主要介绍模型计算流量和输沙率的控制方程和基本原理。
Hydro Trend3.0.4水文模型使用经典的水量平衡原理来模拟流域出口的流量数据，河流出口的流量是由流域面积(A)、降雨量(P)、蒸散量(Ev)、地下水存储蓄(Sr)所控制的。与此同时模型综合考虑降雨(Qr)、雪融水(Qn)、冰融水(Qice)、蒸发(QEv)和地下水通量(Qg)五个基本的径流过程来模拟河流出口的流量。
$$ Q = A\sum\nolimits_{i = 1}^{ne} {\left( {{P_i} - E{\upsilon _i} \pm S{r_i}} \right)} $$ (1) $$ Q = {Q_r} + {Q_n} + {Q_{ice}} - {Q_{E\upsilon }} \pm {Q_g} $$ (2) (1) 式中：ne代表模拟的时间段，i为每日的时间序列。
Hydro Trend3.0.4水文模型使用半经验关系模型来模拟流域出口的悬移质和推移质通量，且模型较原版本而言，综合考虑了冰雪面积和沉积物输运量之间的经验关系，将冰雪消融导致的沉积物输运量计算进入总的沉积物输运量之中。
$$ Q{s_T} = Qs + Q{s_G} $$ (3) (3) 式中：Qs是流域气候特征和流域属性所决定的悬移质输运量, QsG是流域冰雪覆盖所产生的悬移质输运量，QsT为流域出口总的悬移质输运量。
模拟流域长时间尺度沉积物通量可以选择基于流域面积、流量、最大高差和温度的ART和QRT模式，也可以选择综合包含流域岩性和人类活动的BQART模式。基于研究目标本文选用BQART模式模拟流域出口的沉积物通量。第一步，根据半经验关系模型计算长周期输沙率平均值(30years:Qs[kg/s])。
$$ {Q_S} = \overline \omega B{\overline Q ^{0.31}}{A^{0.5}}R{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} T \ge 2^\circ {\rm{C}} $$ (4a) $$ {Q_S} = 2\overline \omega B{\overline Q ^{0.31}}{A^{0.5}}R{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} T＜2^\circ {\rm{C}} $$ (4b) $$ B = IL\left( {1 - {T_E}} \right){E_h} $$ (5) $$ Q{s_G} = \left( {1 - \frac{{V{s_g}}}{{{P_g}}}} \right)\frac{{\sum\nolimits_{a = 1}^n {Q{s_{G\left( a \right)}}} }}{n} $$ (6) 式中：ω=0.02kg/(s.km2.℃)，Q为流域多年平均流量，A为流域面积，R流域内最大高程差，T为流域多年平均气温，L流域岩性因子，Eh人类活动侵蚀因子，TE为水库对悬浮物的捕获系数。(6)中Vsg为每年的冰川存储量(降雨和降雪转化为冰的数量)，Pg为直接降落在冰川区域的降雨量。
第二步，确定长周期时间段输沙率平均Qs后，基于流域出口水沙经验比值关系生成随机模型(Psi)下的日输沙率序列Qs(i):
$$ \frac{{{Q_{s\left( i \right)}}}}{{{Q_{sT}}}} = \varphi \left( i \right){\left( {\frac{{{Q_{\left( i \right)}}}}{{\overline Q }}} \right)^{c\left( a \right)}} $$ (7) (7) 式中：QsT为多年平均输沙率，Q为多年平均流量。φ(i)从随机分布中采样且随机分布的均值为1，C(a)也是从随机分布中采集，该随机分布的均值为流域出口水沙关系的相关指数。且φ(i)、C(a)的均值和方差的计算方法决定了Psi模型模拟的年内和年间的水沙变化。
3. 结果和讨论
3.1 模拟结果验证
实测数据与模拟结果进行对比可以发现(图 5，表 1)，流域6—9月份实测入海径流量占总径流量的48%，而同期模拟值为49%；6—9月实测输沙量占年输沙量的90%，而同期模拟值为86%。统计分析的结果也表明，荒沟站和梨树沟站的月均实测流量和模拟值的偏差分别为1.5%和4%，而两站的月均输沙率的偏差则分别为7.4%和7.7%。逐日或逐月尺度的径流和输沙的模拟将会出现一定幅度的波动，该波动累计到年或十年尺度的量，才可以更好地反映模拟结果的有效性, 而通过对比两个水文站模拟累积水沙通量与实测累积水沙通量的结果也可发现，无论在入海水沙通量的逐年变化幅度方面，还是在一定时间段内的累积入海水沙通量方面都具有较高的一致性，因此模型的模拟结果较好的反应了流域入海水沙通量的实际变化。
![]()
图 5 1958—2012月均流量(及累计)，输沙率(及累计)模拟和实测对比结果Figure 5. Visual comparison of monthly measured value versus simulated value from water discharge(accumulation) and sediment load(accumulation) during 1958-2012 in two hydrographic stations表 1 不同模式下流域年均入海水沙通量变化(气候变化驱动M1、人类活动驱动M2、气候与人类共同影响驱动M3只考虑水土流失、气候与人类共同影响驱动M3综合考虑水土流失+水库拦截)Table 1. Comparison of annual water discharge and sediment load variations forced with different scenarios (climate change-driven M1, human activities-driven M2, combined climate change and human activities-driven M3 soil erosion, combined climate change and human activities-driven M3-soil erosion and dams)流域 年均径流量(108m3/a) 年均输沙量(Mt/a) M1 M2 M3 M1 M2 M3 M3(水土流失) M3(水土流失+水库拦截) 瑷河梨树沟 32.9 33.1 32.9 0.59 0.83 0.80 0.80 鸭绿江荒沟 231.6 231.2 230.8 4.94 0.79 5.93 0.79 全流域 264.5 264.3 263.7 5.53 1.62 6.73 1.59 3.2 不同情景模式下的流域水沙通量变化模拟
3.2.1 不同情景模式下的流域入海径流量变化模拟
本文设定了三种情景模式来模拟鸭绿江流域入海水沙通量的变化，即，气候变化驱动(M1)、人类活动驱动(M2)以及气候变化和人类活动共同驱动(M3)。如图 6所示，M1模式与M3模式下的荒沟站、梨树沟站和整个流域径流量的变化曲线基本重合，且进一步对比实测数据与M1、M3驱动下整个流域的径流量变化可得出相同的结果：1958—2012年M1、M3和实测数据的径流量分别为264.5×108、263.7×108和259.2×108m3/a。由此可见，人类活动对于鸭绿江流域的径流量变化影响很小，其径流量变化是由气候变化所主导。
![]()
图 6 1958—2012年不同模式下的流域年入海径流量和输沙量曲线变化Figure 6. Annual water discharge and sediment load variations forced with three scenarios in the period from 1958 to 2012作为典型的山溪性中小型河流，鸭绿江入海径流量对于流域降水量的变化响应更为敏感和迅速。而流域的降水存在丰水和枯水的周期变化，例如辽河流域的丰、枯水期变化周期约为8~14年[20]。因此，本文拟通过对比不同丰水期和枯水期内降雨量、径流量和输沙量的关系，来进一步研究气候变化对流域入海水沙通量的影响。结果表明(表 2)，1958—2012年荒沟站丰水期与枯水期降雨量分别为913mm和810mm，M1驱动下两时期的径流量分别为270.1×108m3/a和221.8×108m3/a；梨树沟站具有相同的变化趋势，其丰水期与枯水期降雨量分别为1018mm和817mm，径流量分别为37.7×108m3/a和25.7×108m3/a。相较于丰水期，气候变化导致荒沟站与梨树沟站降雨量分别减少了11.3%和19.8%；径流量分别减少了18.8%和31.6%。由此可以看出，相较于长江等大河流域，小流域对于气候的变化(降雨量)响应更加迅速。高抒也得出了相同的结论，认为相对于大河流域及其河口-海岸-陆架区规模巨大的沉积体系，中小型河流的流域面积，以及河口-陆架区沉积体系的规模相对要小，对流域气候变化和人类活动的响应更加迅速[5]。综合上述分析可以看出，人类活动对于流域的径流量影响较小，流域的径流量是由气候变化所主导，且小流域对于气候变化的响应更加迅速。
3.2.2 不同情景模式下的流域入海输沙量变化模拟
气候变化和人类活动可通过改变流域水土流失状况与泥沙输运过程影响流域的入海输沙量[21]。三种模式驱动下的鸭绿江流域年均入海输沙量变化如图 6和表 1所示。在M2驱动下，荒沟站输沙量的变化曲线与M3驱动下输沙量变化曲线接近重合，两种模式驱动下整个流域1958—2012年的输沙量分别为1.62Mt/a和1.59Mt/a。M1驱动下荒沟站输沙量与M3驱动下输沙量相差巨大：1958—2012年，两种模式下整个流域的输沙量分别为5.53Mt/a和1.59Mt/a。通过上述分析可以得出，人类活动导致整个流域的输沙量减少了3.94Mt/a, 比M1驱动模式下入海输沙量减少了71.2%。Zhang等[2]研究也表明，近50年，人类活动对于长江流域产沙的影响越来越大，且长江流域的人类活动(植被的减少、大坝的建设等)对于入海输沙量的影响要远大于径流量产生的影响。通过进一步的计算，若鸭绿江流域未修建大坝，整个流域的入海输沙量为6.73Mt/a，约为现在鸭绿江年均入海输沙量的5倍。
对于山溪性的中小型河流而言，即使气候变化没有明显的改变趋势，不同年际间的降雨量变化(即暴雨和干旱的频率和强度)也是影响其入海输沙量的一个重要因素[22]。本文的研究结果也表明(表 2)，M1驱动下1958—2012年整个流域丰水期与枯水期的输沙量分别达到7.96Mt/a和4.68Mt/a, 枯水期的年均输沙量仅为丰水期的58.8%，而对比M3驱动下整个流域丰水期和枯水期的输沙量分别为2.17Mt/a和1.13Mt/a，枯水期的年均输沙量占丰水期的52.1%。综合上述分析可知，近100年来，虽然人类活动是导致鸭绿江流域入海输沙量变化的最主要原因，但气候变化仍可通过降雨的分布模式对入海输沙量产生影响。
表 2 不同模式下(气候变化驱动、气候与人类共同影响驱动)流域丰、枯水期年均水沙通量变化Table 2. Comparison of annual water discharge and sediment load variations forced with different scenarios (climate change-driven, combined climate change and human activities-driven) in high precipitation year to that in low precipitation year水文站 降雨丰、枯期 时间 降雨(mm) 年径流量(108m3/a) 年输沙量(Mt/a) M1 M3 M1 M3 瑷河 丰水期 1958-1974 1058.7 39.7 39.8 0.77 0.96 枯水期 1975-1982 854.5 25.6 23.3 0.52 0.85 梨树沟 丰水期 1983-1998 978.8 35.6 35.2 0.62 0.88 枯水期 1999-2012 779.3 25.8 26.0 0.37 0.48 鸭绿江 丰水期 1958-1964 917.1 301.2 300.9 10.08 2.36 枯水期 1965-1979 800.0 240.2 242.8 4.81 0.57 荒沟 枯水期 1980-2002 820.5 203.3 198.0 3.65 0.55 丰水期 2003-2012 910.3 239.0 239.3 4.45 0.59 3.3 不同人类活动强度下的子流域水沙通量对比分析
荒沟站与梨树沟站具有相同的气候、地理与地质特征，但两水文站的人类活动强度和人类活动形式存在巨大差别。荒沟水文站以上流域自1940年就已建立大型水库，进入1960年代后又陆续在各支流和干流修建水库，流域存在梯级水库。梨树沟水文站所代表的瑷河流域在2010年前没有大型水库修建，且两个子流域林地覆盖度相差不大(图 4b)。而实测数据也表明虽然瑷河流域面积仅为荒沟集水面积的1/10，但两者的年均输沙量却相同(0.74Mt/a)，因此，两者人类活动影响的差异主要体现在水库对入海沉积物的拦截效应。
径流模数是单位面积上单位时间所产生的径流量，可消除流域面积大小的影响，说明与自然地理相联系的径流特征；输沙模数指河流某断面以上单位面积上所输移的泥沙量，输沙模数的高低可反映流域水土流失的强弱[4]。M1驱动下，1958—2012年荒沟站与梨树沟站的径流模数分别为0.013m3/km2s和0.018m3/km2s，而输沙模数分别为88t/km2和104t/km2，两者没有显著的差异。但在M3共同驱动模式下，虽然荒沟站与梨树沟站径流模数仍相差不大(0.013m3/km2·s和0.018m3/km2·s)，但荒沟站的输沙模数却远低于梨树沟(13t/km2和131t/km2)。与单纯的M1驱动模式相比，在M3驱动模式下的荒沟站输沙模数减少了85%，而梨树沟的输沙模数却增加了26%。流域内人类活动强度和形式的差异，是导致这种差异的最主要原因。
鸭绿江流域近90%左右的输沙集中在洪季，自然状态下流域输沙的主导事件是流域频发的山洪。已有研究结果表明梯级水库平均削减洪峰流量约50%，其中水丰水库削减洪峰介于33%~89.5%之间[23]，削减洪峰的同时梯级水库拦截了大量的入海沉积物。相对于M1，在M3驱动下，荒沟水文站以上鸭绿江干流流域人类活动导致输沙量减少4.15Mt/a，分析M1与M3(水土流失)数据可知，流域水土流失造成输沙量增加了0.99Mt/a，而通过分析M3(水土流失)和M3(水土流失+水库拦截)可得出，大坝拦截造成输沙量减少5.14Mt/a，即，水库拦截导致荒沟站的年均输沙量减少了86.7%(表 1)。而梨树沟站所代表的瑷河流域的输沙量变化主要由林地覆盖率的变化所主导，相对于M1，在M3驱动下水土流失导致其输沙量增加了0.21Mt/a。从整个鸭绿江流域来看，人类活动导致流域的入海输沙量减少了3.94Mt/a，其中，水库拦截效应导致输沙量减少了5.14Mt/a，而水土流失使流域输沙量增加了1.20Mt/a，由于水库修建，导致鸭绿江流域的年均入海输沙量减少了76.4%(表 1)。据此可以估算，假设水库拦截的沉积物全部输运到河口，则鸭绿江流域未修建水库前(1940年水丰水库修建)入海年均输沙量约为6.73Mt/a。
4. 结论
模拟结果验证表明，HydroTrend3.0.4水文模型可以较好的模拟流域年间和年内的水沙特征变化，模拟值与实测值的偏差在7%以内。研究表明，1958—2012年间流域实测入海径流量与输沙量分别为259.2×108m3/a和1.48Mt/a，而模拟气候变化、人类活动以及气候变化和人类活动共同驱动下流域年入海径流量分别为264.5×108、264.3×108和263.7×108m3/a，年入海输沙量分别达到5.53、1.62和1.59Mt/a。
进一步分析发现，不同年际间的降雨量变化(即暴雨和干旱的频率和强度)是导致鸭绿江流域入海径流量变化的最主要原因，且对于山溪性的中小河流，气候变化可通过降雨的分布模式对输沙量变化产生一定的影响。
人类活动是鸭绿江流域入海输沙量变化的主要原因，流域人类活动对年入海输沙量变化的贡献量为3.94Mt/a，其中流域林地覆盖变化导致水土流失量增加1.20Mt/a, 而流域大坝拦截量为5.14Mt/a，鸭绿江流域的梯级水库拦截了流域76.4%左右的入海输沙量。考虑到流域入海输沙量主要受人类活动所主导，因此，1940年鸭绿江流域修建水库前，其年入海输沙量达到6.73Mt/a，约为现在鸭绿江年入海输沙量的5倍。
-
[1] Folk R L,Ward W C.Brazos river bar:a study in the significance of grain size parameters[J].Journal of Sedimentary Petrology,1957,(27):3-27.
[2] Visher G S.Grain size distribution and depositional processes[J].Journal of Sedimentary Petrology,1969,(39):1074-1106.
[3] Gao S,Collins M.Net sediment transport patterns inferred from grain size trends,based upon definition of "transport vectors"[J].Sedimentary Geology,1992,81:47-60.
[4] 高抒,Collins M.沉积物粒径趋势与海洋沉积动力学[J].中国科学基金,1998(4):241-246.[GAO Shu,Collins M.The use of grain size trends in marine sediment dynamics[J].Science Foundation in China,1998 (4):241-246.]
[5] 成都地质学院陕北队.沉积岩(物)粒度分析及其应用[M].北京:地质出版社,1978.[Chengdu College of Geology,Shanbei team.Sedimentary rock (material) particle-size analysis and its application[M].Beijing:Geological Publishing House,1978.] [6] 孙有斌,高抒,鹿化煜.前处理方法对北黄海沉积物粒度的影响[J].海洋与湖沼,2001,32(6):665-671. [SUN Youbin,GAO Shu,LU Huayu.Influence of different pretreatment procedures on the particle-size distribution of surficial sediments in the northern Yellow Sea[J].Oceanologia et Limnologia Sinica,2001,32(6):665-671.]
[7] 王德杰,范代读,李从先.不同预处理对沉积物粒度分析结果的影响[J].同济大学学报,2003,31(3):314-318. [WANG Dejei,FAN Dadu,LI Congxian.Influence of different pretreatments on size analysis and its implication[J].Journal of Tongji University,2003,31(3):314-318.]
[8] 中华人民共和国国家质量监督检验检疫局,中国国家标准化管理委员会.海洋地质调查规范第8部分:海洋地质地球物理调查(GB/T 12763.8-2007 )[S].2007:1-79.[General Administration of Quality Supervision,Inspection and Quarantine of the People's Republic of China,Standardization administration of the People's Republic of China.Specifications for Oceanographic Survey-Part8:Marine Geology and Geoiphysics survey.(GB/T 12763.8-2007)[S].2007:1-79.]
[9] 贾建军,高抒,薛允传.图解法与矩法沉积物粒度参数的对比[J].海洋与湖沼,2002,33(6):577-582. [JIA Jianjun,GAO Shu,XUE Yunchuan.Grain-size parameters derived from graphic and moment methods:A comparative study[J].Oceanologia et Limnologia Sinica,2002,33(6):577-582.]
[10] 李春初.大口河海岸特性和黄骅港的扩建问题[J].港口工程,1993(5):6-14.[LI Chunchu.Coastal characteristic of Dakouhe River and the expansion of Huanghua Harbor[J].Harbor Engineering,1993 (5):6-14.]
[11] McManus J.Grain size determination and interpretation[M]//Techniques in Sedimentlogy.Oxford,Blackwell,1988:63-85.
[12] 徐家声,刘鹤智.黄骅海岸带开发环境概论[M].北京:海洋出版社,1991.[XU Jiasheng,LIU Hezhi.Coastal Zone Development Environment of Huanghua[M].Beijing:China Ocean prass,1991.] [13] 张庆河,王崇贤,杨华,等.黄骅港海域表层泥砂特性及其影响[J].中国港湾建设,2004(4):14-17.[ZHANG Qinghe,WANG Chongxian,YANG Hua,et al.Property and influence of sediments on seabed surface around Huanghua Harbor[J].China Harbour Engineering,2004 (4):14-17.]
[14] 张庆河,林全泓.黄骅港海域沿堤流现象遥感图像分析[J].中国港湾建设,2004(1):29-31.[ZHANG Qinghe,LIN Quanhong.Investigation of currents along breakwater for Huanghua Port with help of satellite romote sensing images[J].China Harbor Engineering,2004 (1):29-31.]
[15] 王宏,李建芬,康慧,等.泥质海岸带现代地质作用(沉积,剥蚀及岸线变迁)与精细测年[R].天津地质矿产研究所地调项目研究报告,2002:1-62.[WANG Hong,LI Jianfen,KANG Hui,et al.Modern Sedimentation of the Muddy Coast (Accumulation Rate,Erosion and Coastline Change) and High Resolution Dating Report[R].Tianjin Institute of Geology and Mineral Resources,China,2002,1 -82.]
[16] 李建芬,王宏,夏威岚,等.渤海湾西岸210Pbexc和137Cs测年与现代沉积速率[J].地质调查与研究,2003,26(2):114-128. [LI Jiafen,WANG Hong,XIA Weilan,et al.210Pbexc and 137Cs dating and modern sedimentation rate on the western coast of Bohai Bay[J].Geological Survey and Research,2003,26(2):114-128.]
[17] 王福.渤海湾海岸带210Pb、137Cs示踪与测年研究:现代沉积及环境意义[D].中国地质科学院博士学位论文,2009.[WANG Fu.210Pbexc,137Cs tracing and dating on Bohai Bay Coast:modern sedimentation and implication[D].PhD thesis,Chinese Academy of Geological Sciences,China,2009. ]
[18] 王福,裴艳东,刘志广,等.137Cs示踪的黄骅港两侧海区的沉积特征[J].地质通报,2008,27(7):1054-1059. [WANG Fu,PEI Yangdong,LIU Zhiguang,et al.Sedimentary characteristics of Huanghua Harbour traced by 137Cs[J].Geological Bulletin of China,2008,27(7):1054-1059.]
[19] 王福,王宏,王云生,等.黄骅港地区岸线变迁及滩面(浅海区形态变化研究[J].海洋地质与第四纪地质,2010,30(5):25-32. [WANG Fu,WANG Hong,WANG Yunsheng,et al.Shoreline changes and morphology changes of tidal flat-shallow sea area of Huanghua Harbor[J].Marine Geology and Quaternary Geology,2010,30(5):25-32.]
[20] 王福,王宏,王云生,等.黄骅港地质环境研究报告[R].河北省国土资源厅/海洋局、中国地质调查局天津地质调查中心,2008.[WANG Fu,WANG Hong,WANG Yunsheng,et al.Study report of Geological environment in Huanghua Harbor[R].Hebei Province Department Land&Resources/SOA,Tianjin Center of Geological Survey,CGS,2008.]
下载:
计量
- 文章访问数: 1832
- HTML全文浏览量: 308
- PDF下载量: 16
