摘要:現(xiàn)代黃河三角洲位于海陸交互地帶,為我國(guó)重要的生態(tài)保護(hù)區(qū),同時(shí)也是石油開采地。受河流、潮流、人類活動(dòng)影響,現(xiàn)代黃河三角洲沉積物粒度組成復(fù)雜,物質(zhì)來(lái)源多樣。本文共采集了黃河三角洲表層沉積物樣品219個(gè),首先采用粒度分析方法分析了現(xiàn)代黃河三角洲沉積物的粒度特征以及沉積物來(lái)源,然后采用參數(shù)端元模型分析方法對(duì)沉積物的物質(zhì)來(lái)源進(jìn)行了劃分,最后采用聚類分析方法分析了表層沉積物類型。結(jié)果表明:現(xiàn)代黃河三角洲沉積物粒度組成以粉砂(體積分?jǐn)?shù)71.10%)和砂(27.62%)為主,大部分沉積物粒度分選中等,粒度分布曲線呈近對(duì)稱分布。沉積物劃分為EM1、EM2、EM3三個(gè)沉積端元,分別代表三種不同物質(zhì)來(lái)源,其中EM1代表舊河流沉積作用,EM2代表海洋潮流沉積作用,EM3代表現(xiàn)黃河河流沉積作用。將表層沉積物分成三種沉積物類型,結(jié)合聚類分析結(jié)果將黃河三角洲分為三類沉積環(huán)境,其中:第一類主要分布在黃河三角洲等西部河流的中上游,包括馬新河、沾利河、草橋溝、挑河,主要為舊河流沉積,水動(dòng)力較弱;第二類主要分布在刁口河、神仙溝、現(xiàn)黃河(清水溝)等東部河流,現(xiàn)黃河河流沉積作用較強(qiáng);第三類主要分布在黃河三角洲北部沿岸及潮灘,受海洋潮流作用較強(qiáng),海岸受到侵蝕。此外,人類活動(dòng)對(duì)黃河三角洲岸線變化也產(chǎn)生了一定影響。
關(guān)鍵詞:黃河三角洲;表層沉積物;粒度;端元;物源分析
doi:10.13278/j.cnki.jjuese.2023020
中圖分類號(hào):P736.21
文獻(xiàn)標(biāo)志碼:A
Supported by the National Natural Science Foundation of China (42271162,41971100),the Natural Science Foundation of Shandong Province (ZR2020MD063) and the Joint Fund Key Support Project of NSFCShandong Province (U1706214)
Grain Size Characteristic and Source of Sediments
in Modern Yellow River Delta
Yu Jihuai1,Jie Dongmei 3,4,5,Li Ping6,7
1. School of Geographical Sciences, Northeast Normal University, Changchun 130024, China
2. Institute for Peat and Mire Research, Northeast Normal University, Changchun 130024, China
3. State Environmental Protection Key Laboratory of Wetland Ecology and Vegetation Restoration, Changchun 130024, China
4. Key Laboratory of Geographical Processes and Ecological Security of Changbai Mountains (Northeast Normal University), "Ministry of Education, Changchun 130024, China
5. Key Laboratory of Vegetation Ecology, Institute of Grassland (Northeast Normal University), Ministry of Education, "Changchun 130024, China
6. First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, Shandong,China
7. Key Laboratory of Coastal Zone Science and Integrated Management, Ministry of Natural Resources, Qingdao 266061, "Shandong,China
Abstract: There are some important ecological reserves and oil exploitations in modern Yellow River delta, which located at the junction of sea and land. Influenced by rivers, currents and human activities, the modern Yellow River delta has a complex composition of sediment grain size and sources. A total of 219 sediment samples from the Yellow River delta were collected in this paper. Firstly, the grain size distribution and material source of the modern Yellow River delta sediment were analyzed by particle size analysis method. Then the parametric endmember model analysis method was used to divide the material sources of sediments. Finally, the cluster analysis method was used to analyze the surface sediment types. The results show that the grain size composition of delta sediments is mainly silt (71.10%) and sand (27.62%), most of the sediments are of medium grain size sorted and the curve of the grain size distribution is nearly symmetrical. The delta sediments can be divided into three types of end-mumber such as EM1, EM2 and EM3, which represent three different material sources. Among them, EM1 is from fluvial sediments, EM2 is from tidal current sediments and EM3 is from the current Yellow River fluvial sediments. According to the results of cluster analysis, the Yellow River delta is divided into three types of sedimentary environments. The first type is mainly distributed in the middle and upper reaches of western part of the Yellow River delta rivers, including Maxin River, Zhanli River, Caoqiaogou River and Tiaohe River, which is mainly deposited by old river with weak hydrodynamic force. The second type is mainly distributed in the eastern part of rivers such as Diaokou River, Shenxiangou and the current Yellow River (Qingshuigou), and the current Yellow River has strong sedimentary dynamics. The third type is mainly distributed along the northern coast and tidal flat of the Yellow River delta, which is strongly affected by the marine tidal and the coastline is eroded. In addition, human activities also had a certain impact on the changes of the coastline in the Yellow River delta.
Key words: Yellow River delta; surface sediments; grain size; end mumber; source analysis
0 引言
黃河作為我國(guó)第二大河,具有水少沙多、水沙異源、快速淤積的特點(diǎn)。黃河平均每年向渤海輸送8.7億t泥沙,其中70%沉積在河口附近并形成新的陸地[1]。1855年,黃河從銅瓦廂決口奪大清河匯入渤海,形成現(xiàn)代黃河三角洲[2]。近年來(lái),受氣候變化以及人類活動(dòng)的影響,黃河三角洲下游部分河流發(fā)生淤積,有的地段岸線受海洋侵蝕作用影響,發(fā)生蝕退[3]。為保護(hù)海岸生態(tài)環(huán)境,保障沿海居民生產(chǎn)生活安全,學(xué)者們對(duì)現(xiàn)代黃河三角洲潮間帶及水下三角洲地區(qū)展開了較多研究,研究?jī)?nèi)容包括沉積物粒度特征及物質(zhì)來(lái)源[48]、沉積特征與沉積環(huán)境[913]、黃河泥沙輸送對(duì)地形地貌影響[1416]、人類活動(dòng)影響[1718]、岸線變化[1921]、濕地健康評(píng)價(jià)[22]等。相比于黃河三角洲潮間帶及水下三角洲,黃河三角洲陸上部分沉積環(huán)境的研究工作較少,而且學(xué)界對(duì)于黃河三角洲地區(qū)沉積動(dòng)力與物質(zhì)來(lái)源的關(guān)系認(rèn)識(shí)也不清楚,因此,有必要在黃河三角洲開展表層沉積物水動(dòng)力條件及物質(zhì)來(lái)源分析。
粒度作為沉積學(xué)研究中的重要指標(biāo),可以反映沉積物質(zhì)的搬運(yùn)方式、物質(zhì)來(lái)源、沉積動(dòng)力等多重信息,對(duì)研究沉積物輸運(yùn)過(guò)程以及沉積環(huán)境演變具有重要意義[23]。由于河流沉積物粒度組成及分布受水動(dòng)力環(huán)境和不同物質(zhì)來(lái)源的綜合影響,因此利用粒度數(shù)據(jù)分析沉積物特征及物質(zhì)來(lái)源時(shí),常常需要對(duì)粒度組分進(jìn)行分離,其中粒度端元分析方法是區(qū)分粒度組分的有效手段。粒度端元分析通過(guò)運(yùn)用數(shù)學(xué)算法將多種沉積動(dòng)力作用下形成的沉積物粒度數(shù)據(jù)有效分離成多個(gè)獨(dú)立的具有某特定特征的粒度組分,從而判斷沉積物來(lái)源和識(shí)別沉積動(dòng)力條件[24]。端元分析方法自提出以來(lái),模型經(jīng)過(guò)了不斷地改進(jìn),在沉積學(xué)分析上應(yīng)用越來(lái)越廣泛[2526]。比如:張曉東等[27]對(duì)長(zhǎng)江口臨近海域表層沉積物進(jìn)行了粒度端元分析,并對(duì)該地區(qū)沉積動(dòng)力環(huán)境進(jìn)行了描述;林鎮(zhèn)坤等[28]對(duì)南流江河口水下三角洲進(jìn)行端元分析并判斷了當(dāng)?shù)氐某练e動(dòng)力環(huán)境;薛成鳳等[29]運(yùn)用端元分析推測(cè)了中小河流對(duì)長(zhǎng)江水下三角洲遠(yuǎn)端泥沉積的貢獻(xiàn)。另外,在進(jìn)行沉積動(dòng)力分區(qū)和探討沉積環(huán)境方面[3031],也常常運(yùn)用聚類分析的方法,結(jié)合粒度特征參數(shù)(粒度參數(shù)、粒度組分)對(duì)沿海地區(qū)沉積物進(jìn)行沉積環(huán)境劃分[3233]。因此,粒度指標(biāo)可用于判斷黃河三角洲表層沉積物的物質(zhì)來(lái)源和動(dòng)力環(huán)境。
本文通過(guò)對(duì)黃河三角洲陸上沉積物的粒度特征及物質(zhì)來(lái)源進(jìn)行研究,探討不同沉積動(dòng)力對(duì)黃河三角洲岸線及環(huán)境的影響,以期為黃河三角洲地區(qū)岸線防護(hù)及環(huán)境治理提供參考。
1 區(qū)域概況與研究方法
1.1 區(qū)域概況
黃河三角洲位于渤海盆地東部的濟(jì)陽(yáng)斷陷和埕寧隆起之間,基底構(gòu)造以強(qiáng)烈的褶皺和斷層為特征。黃河三角洲地勢(shì)總體西南高,東北低(圖1a)。1855年以來(lái),黃河尾閭發(fā)生多次改流,并在不同時(shí)期形成多條分流河流。近幾十年來(lái)受人工改流控制,黃河三角洲存在多條河流并行現(xiàn)象。除現(xiàn)黃河外,河流多為南北走向,河流自西向東分別為馬新河、沾利河、草橋溝、挑河、刁口河、神仙溝、現(xiàn)黃河(圖1b),其中馬新河、沾利河、草橋溝為人工開挖的河流。黃河三角洲地區(qū)屬于暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候,四季差異分明,年均溫12 ℃,年降水量776 mm。海岸潮汐類型多樣,大部分為不正規(guī)半日潮[2]。植被組成以草本植物為主,木本植物所占比例較少。黃河三角洲植被的主要建群種和優(yōu)勢(shì)種為濕生植物和鹽生植物。植被形成時(shí)間較短,群落穩(wěn)定性較差[34]。
1.2 研究方法
研究區(qū)位于黃河下游、黃河口附近。采樣點(diǎn)分布于現(xiàn)代黃河三角洲的陸上部分及潮間帶區(qū)域(圖1c)。沿現(xiàn)黃河河流以及神仙溝、刁口河、挑河、草橋溝、沾利河、馬新河等黃河尾閭擺蕩形成較早的舊河流地區(qū)進(jìn)行等距采樣,采樣相鄰點(diǎn)間隔約為5 km,選取表層深度0~2 cm的沉積物樣品共219個(gè)。舊河流基本信息見表1。
粒度測(cè)定方法:首先將樣品烘干,稱取約1 g樣品;然后量取10%的稀鹽酸溶液以除去樣品中的鈣質(zhì)膠結(jié),待反應(yīng)完全后加入蒸餾水并離心至中性;接著量取約5 mL雙氧水并進(jìn)行水浴加熱,以去除有機(jī)質(zhì),待反應(yīng)完全后加入蒸餾水進(jìn)行水洗并離心直至中性;再向離心管中加入0.05 mol/L的(NaPO3)6約5 mL,用超聲波清洗儀清洗樣品約10 min,使樣品充分分散;最后使用粒度分析儀Microtrac S3500進(jìn)行測(cè)試,測(cè)量范圍為0.025~2 000 μm,重復(fù)測(cè)量3次使相對(duì)誤差在±3%以內(nèi)。
數(shù)據(jù)分析方法:首先,在Matlab中使用Paterson and Heslop提供的AnalySize程序包,采用Gen.Weibull參數(shù)模型對(duì)黃河三角洲表層沉積物粒度各粒級(jí)的體積分?jǐn)?shù)進(jìn)行分析,并對(duì)粒度數(shù)據(jù)進(jìn)行EMA(端元分析)矩陣計(jì)算。然后,用Folk和Ward公式計(jì)算各個(gè)樣品中的平均粒度、峰度、分選系數(shù)、偏度[36]。最后,運(yùn)用IBM SPSS Statistics 中的K均值聚類方法,根據(jù)粒度參數(shù)Mz對(duì)樣品進(jìn)行聚類分析。
2 結(jié)果與分析
2.1 黃河三角洲表層沉積物粒度組成及分布
按國(guó)際通用標(biāo)準(zhǔn)進(jìn)行分類,Φ≤4為砂,4<Φ≤8為粉砂,Φ>8為黏土(其中Φ=-log2D,D為顆粒直徑,mm)[37],得到沉積物組成按照體積分?jǐn)?shù)由高到低分別為粉砂(71.10%)、砂(27.62%)、黏土(1.28%)。從不同河流位置上表層沉積物的粒度組成中粉砂、砂和黏土的平均體積分?jǐn)?shù)來(lái)看,河流間粒度組成差異較明顯(表2)。從表2可知:整體上看,黃河三角洲表層沉積物粉砂體積分?jǐn)?shù)最高,范圍在59.93%~80.29%之間;刁口河、神仙溝、現(xiàn)黃河等東部河流地區(qū)砂體積分?jǐn)?shù)較高,分布于26.56%~38.84%之間;馬新河、沾利河、草橋溝、挑河等舊河流地區(qū)黏土體積分?jǐn)?shù)較高,在 1.82%~2.99%之間。
黃河三角洲地區(qū)表層沉積物砂、粉砂、黏土粒度分布特征(圖2)顯示,現(xiàn)代黃河三角洲研究區(qū)范圍內(nèi)的樣品中,砂和粉砂的體積分?jǐn)?shù)在空間上的分布變化范圍較大,黏土的空間變化較小。不同河流表層沉積物粒度特征具體表現(xiàn)為:黃河三角洲表層沉積物粒度峰度平均值為1.09,近對(duì)稱;分選系數(shù)平均值為1.11,分選中等;偏度系數(shù)為0.10,近對(duì)稱。
2.2 粒度分布曲線及概率累積曲線
沉積物的粒度主要受搬運(yùn)介質(zhì)、搬運(yùn)方式、沉積環(huán)境等多種因素影響[23]。粒度概率累積曲線可以反映沉積物的搬運(yùn)方式。根據(jù)累積曲線的截點(diǎn)以及分段可以判斷物質(zhì)的搬運(yùn)是以懸移、躍移、底移中的哪種搬運(yùn)方式為主,以及不同搬運(yùn)方式所占的比例如何。原則上根據(jù)現(xiàn)代黃河三角洲研究區(qū)的大多數(shù)點(diǎn)的粒度分布曲線和概率累積曲線的形態(tài),將粒度分布曲線及概率累積曲線劃分為三組類型:類型一為馬新河、沾利河、草橋溝、挑河等舊河流;類型二為刁口河、神仙溝、現(xiàn)黃河等東部河流;類型三為沿海地區(qū)海岸及潮灘。將這三組不同類型的沉積物粒度組合取平均值,并做出相應(yīng)地區(qū)的概率累積曲線以及粒度分布曲線圖(圖3)。其中,粒度分布曲線可以反映粒度的分布,進(jìn)而可以推測(cè)沉積物的物質(zhì)來(lái)源。
馬新河、沾利河、草橋溝、挑河等距離海岸有一定距離的舊河流平均粒度為5.26Φ(圖3a),粒級(jí)為細(xì)粉砂,分選系數(shù)為1.41,分選較差,偏度為0.09,粒度分布曲線呈近對(duì)稱分布,峰態(tài)為0.94,峰態(tài)中等。概率累積曲線為二段式(圖3a),可知物質(zhì)搬運(yùn)方式主要為躍移和懸移。
刁口河、神仙溝、現(xiàn)黃河等東部河流平均粒度為3.75Φ(圖3b),粒級(jí)為極細(xì)砂,分選系數(shù)為0.85,分選中等,偏度為0.08,粒度分布曲線呈近對(duì)稱分布,峰態(tài)為1.24,峰態(tài)窄。概率累積曲線為二段式(圖3b),物質(zhì)搬運(yùn)方式主要為躍移和懸移?,F(xiàn)黃河在10Φ左右存在明顯截點(diǎn),說(shuō)明在沉積物的細(xì)端部分物質(zhì)的沉積方式發(fā)生了明顯變化。
在黃河三角洲北部沿海地區(qū),平均粒度為4.62Φ(圖3c),粒級(jí)為粗粉砂,分選系數(shù)為1.03,分選較差,偏度為0.16,粒度分布曲線為正偏,峰態(tài)為1.19,峰態(tài)窄。概率累積曲線為四段式(圖3c),在6Φ和11Φ位置處存在拐點(diǎn);其中躍移組分的斜率較高,反映了沉積物在能量相對(duì)較低的河口上快速堆積的特征,該處沉積物受海洋潮流的影響較為強(qiáng)烈。
2.3 端元分析
對(duì)黃河三角洲表層沉積物粒級(jí)數(shù)據(jù)進(jìn)行EMA(end member analysis)矩陣計(jì)算并得到端元數(shù)與線性相關(guān)性、角度偏差的關(guān)系(圖4),其中R2為各端元之間的判定系數(shù),角度分布曲線之間的偏離程度[24]。評(píng)估端元分解效果時(shí),需要盡量保證 R2在0.9以上,角度偏差盡量小,同時(shí)選擇盡可能少的端元數(shù)目。隨著端元數(shù)量增加,線性相關(guān)增強(qiáng),角度偏差下降。最終選擇3個(gè)端元(EM1、EM2和EM3)對(duì)該組粒度數(shù)據(jù)進(jìn)行反演并分析。
沉積物粒度端元的性質(zhì)通常受沉積物物質(zhì)來(lái)源和沉積環(huán)境(括沉積物的運(yùn)移特征以及水動(dòng)力環(huán)境條件[27])的影響。3個(gè)不同端元所代表的粒度參數(shù)(表3)反映了黃河三角洲整體沉積環(huán)境受3種不同沉積動(dòng)力來(lái)源的影響,結(jié)合每個(gè)樣品各端元的空間分布特征(圖5)可以分析沉積物粒度端元所指示的物質(zhì)來(lái)源及水動(dòng)力環(huán)境。
EM1代表沉積物的平均粒度和眾數(shù)粒度分別為6.13Φ和6.01Φ(表3),屬于細(xì)粉砂粒級(jí),其中黏土體積分?jǐn)?shù)為4.07%,粉砂體積分?jǐn)?shù)為93.89%。分選系數(shù)為1.08,分選程度較中等;偏度為0.02,近對(duì)稱;峰度為0.97,峰態(tài)為中等。從空間位置分布上看,EM1物質(zhì)來(lái)源主要分布在現(xiàn)代黃河三角洲西北部沿岸潮灘以及馬新河、沾利河、草橋溝、挑河等舊河流地區(qū)。因此推測(cè)EM1為舊河流沉積動(dòng)力,代表舊河流在較弱的水動(dòng)力條件下逐漸沉積形成的河流沉積物。
EM2代表沉積物的平均粒度和眾數(shù)粒度分別為4.74Φ和4.51Φ(表3),屬于粗粉砂粒級(jí);分選系數(shù)為0.59,分選較好;偏度為0.11,偏度為正偏,顆粒在粗粒部分集中;峰度為1.00,峰態(tài)為中等。EM2主要分布在黃河三角洲北部沿海地區(qū)。北部廢棄河口岸線由于缺少河流泥沙補(bǔ)給,受波浪、潮流、風(fēng)暴潮等多種海洋因素作用,在刁口河口無(wú)潮點(diǎn)附近的海岸段岸線強(qiáng)烈蝕退[19]。因此推斷EM2受海洋影響強(qiáng)烈,為海洋潮流及波浪動(dòng)力作用下形成的沉積物。
EM3代表沉積物的平均粒度和眾數(shù)粒度分別為3.68Φ和3.76Φ(表3),屬于粗粉砂粒級(jí);分選系數(shù)為0.56,分選較好;偏度為-0.13,為負(fù)偏,顆粒在細(xì)粒部分集中;峰度為0.96,峰態(tài)為中等。與EM1和EM2相比,EM3粒度更粗,且為負(fù)偏,在空間位置分布上看,EM3主要分布在黃河三角洲研究區(qū)的東北側(cè),以及現(xiàn)黃河。在神仙溝及刁口河流靠近陸地地區(qū)EM3含量較高,這些河流曾經(jīng)為黃河流經(jīng)河流,由于在這些地區(qū)黃河改流時(shí)間較短,且河流之間相互連通,因此仍受到黃河沉積作用影響。因此推斷EM3主要來(lái)源為現(xiàn)黃河沉積。
從端元組成上看,EM2和EM3的占比較高,而EM1占比較低(表3)??傮w上是粉砂體積分?jǐn)?shù)高,黏土體積分?jǐn)?shù)低,沉積物粒度粗;說(shuō)明黃河三角洲研究區(qū)沉積環(huán)境水動(dòng)力強(qiáng),沉積環(huán)境受河流和潮流等多種因素影響。近年來(lái)隨著入海泥沙減少,除了現(xiàn)黃河(清水溝)以外,其余地方受到海洋侵蝕作用加強(qiáng),因而EM2占比較高(圖5)。
3 討論
3.1 沉積動(dòng)力分區(qū)及沉積環(huán)境探討
由于現(xiàn)代黃河三角洲粒度分布是與現(xiàn)代河口動(dòng)力環(huán)境相適應(yīng)的,因此,分析沉積物的粒度分布特征可以對(duì)河口動(dòng)力環(huán)境進(jìn)行推斷。
結(jié)合現(xiàn)代黃河三角洲的粒度特征,運(yùn)用K均值聚類方法對(duì)黃河三角洲表層沉積物的平均粒度進(jìn)行分類,共分為三種沉積類型,根據(jù)這三種類型反映的不同沉積特征將沉積物的空間分布劃分為三種沉積動(dòng)力環(huán)境區(qū)(圖6)。對(duì)不同沉積動(dòng)力環(huán)境區(qū)聚類占比及粒度特征進(jìn)行計(jì)算,結(jié)果見表4。
沉積Ⅰ區(qū)主要分布在黃河三角洲西部河流,包括馬新河、沾利河、草橋溝、挑河等河流的上游及中游部分。從聚類組成上來(lái)看,沉積Ⅰ區(qū)中聚類1占比38.64%,聚類2占比13.64%,聚類3占比47.73%;聚類3占比最高,其平均粒度為5.46Φ(表4)。由于河流改流時(shí)間較早,河流廢棄時(shí)間較長(zhǎng),受人類活動(dòng)用水增加以及人工堤壩的修建等影響,沉積Ⅰ區(qū)的沉積物粒度較細(xì),該地區(qū)物質(zhì)搬運(yùn)方式主要以躍移和懸移為主。從岸線變化上看,結(jié)合前人[3]分析認(rèn)為該地區(qū)由于連接港口的公路修建,攔截了潮灘,岸線較為穩(wěn)定。沉積Ⅰ區(qū)的沉積環(huán)境以舊河流沉積為主。
沉積Ⅱ區(qū)主要分布在刁口河、神仙溝、現(xiàn)黃河(清水溝)等東部河流。從聚類組成上來(lái)看,沉積Ⅱ區(qū)聚類1占比44.83%,聚類2占比48.28%,聚類3占比0.07%,主要以聚類2和聚類1為主。聚類2和聚類1平均粒度分別為4.03Φ和4.79Φ(表4)。該區(qū)域沉積環(huán)境主要為黃河河流沉積作用?,F(xiàn)行河口區(qū)主要受黃河水動(dòng)力影響。一般來(lái)說(shuō),水動(dòng)力環(huán)境越強(qiáng),粒度越粗。由于黃河徑流量較大,流速較快,因而沉積物粒度較粗。其中,徑流區(qū)域影響下砂的含量較高。這可能是由于河流流速較快,細(xì)顆粒物質(zhì)隨水流被搬運(yùn)至遠(yuǎn)處,粗顆粒物質(zhì)搬運(yùn)距離較近,多在原地沉積。袁萍等[6]認(rèn)為現(xiàn)行河口三角洲葉瓣周圍的表層沉積物粒度較遠(yuǎn)離河口位置的沉積物粒度粗,且現(xiàn)黃河口沉積物粒度有粗化的趨勢(shì)。從岸線變化上看,結(jié)合前人[11]分析認(rèn)為黃河三角洲岸線在孤東以北和清水溝流路南側(cè)的區(qū)域岸線作用以侵蝕為主,而在孤東海域和清水溝流路河口處則以淤積作用為主。
沉積Ⅲ區(qū)主要分布在黃河三角洲北部沿岸及潮灘地區(qū)。從聚類組成上來(lái)說(shuō),沉積Ⅲ區(qū)聚類1占比67.82%,聚類2占比12.64%,聚類3占比19.54%。主要以聚類1為主,其平均粒度為4.79Φ(表4)。該區(qū)域沉積環(huán)境主要為海洋沉積作用,岸線蝕退。廢棄刁口河口和神仙溝屬于波流控制型[16]。自1976年黃河改流清水溝流路以來(lái),三角洲海岸線總的變化趨勢(shì)是現(xiàn)行河口海岸線明顯向海淤進(jìn),而北部廢棄河口岸線則明顯侵蝕后退[19]。由于潮灘地區(qū)地勢(shì)平坦,河流流速迅速減弱,水動(dòng)力環(huán)境變?nèi)?,因而其粒度較河流地區(qū)細(xì)。黃河改流后,黃河三角洲北部潮間帶部分岸灘波浪作用加強(qiáng),侵蝕效應(yīng)亦增強(qiáng)[4]。在黃河三角洲北部沿岸沉積物主要受海洋潮流作用影響。其中在近岸,漲落潮的往復(fù)運(yùn)動(dòng)對(duì)細(xì)顆粒泥沙的輸移和分布有重要影響。當(dāng)來(lái)自海洋的潮流與受風(fēng)驅(qū)動(dòng)下的水流、河流和海洋環(huán)流相互作用時(shí),沉積物質(zhì)的凈運(yùn)輸方向會(huì)發(fā)生相應(yīng)改變。隨著水深沿海岸變化,海底摩擦力增大,潮流速度向岸邊減小。細(xì)粒沉積物的搬運(yùn)受侵蝕和沉積滯后效應(yīng)的控制,呈現(xiàn)出向岸邊搬運(yùn)和沉積的趨勢(shì)。受潮流頂托作用,沉積物細(xì)顆粒物質(zhì)在海岸沉積,因而在黃河三角洲近岸粒度較細(xì)。近年來(lái),黃河徑流減少導(dǎo)致海水潮流動(dòng)力增強(qiáng)[15]。因而,在黃河三角洲北部海岸線有侵蝕的趨勢(shì)。
3.2 物質(zhì)來(lái)源探討
黃河三角洲沉積物物質(zhì)來(lái)源主要為現(xiàn)黃河河流沉積、舊河流沉積和海洋潮流沉積三個(gè)來(lái)源。黃河三角洲地區(qū)沉積物的粒度特征與物質(zhì)來(lái)源有密切關(guān)系。在黃河三角洲地區(qū),受現(xiàn)黃河河流沉積和受海洋潮流沉積影響的區(qū)域,沉積物粒度較粗,而受舊河流沉積影響的地區(qū)沉積物粒度較細(xì)。從黃河三角洲沉積物端元組成上看,現(xiàn)代黃河三角洲不同位置的端元構(gòu)成存在差異:在馬新河、沾利河、草橋溝、挑河等西部河流,端元組成以EM1為主,沉積物主要來(lái)源為舊河流沉積作用下的泥沙沉積;在刁口河、神仙溝、現(xiàn)黃河等東部河流,端元組成以EM3為主,沉積物主要物質(zhì)來(lái)源為現(xiàn)黃河河流沖刷下來(lái)的黃土高原的黃土[8];在黃河三角洲沿海地區(qū),端元組成以EM2為主,沉積物主要物質(zhì)來(lái)源可能為海洋和潮流挾沙堆積而成。這說(shuō)明該區(qū)域主要受黃河河流沉積作用影響,海洋潮流僅影響黃河入??诒辈康貐^(qū)。
本文所討論的陸上黃河三角洲沉積動(dòng)力形式與前人關(guān)于黃河三角洲潮間帶及水下三角洲沉積動(dòng)力形式較為一致,如陳小英等[13]認(rèn)為黃河三角洲物質(zhì)來(lái)源受到徑流、潮流以及海洋風(fēng)浪影響,喬淑卿等[7]認(rèn)為黃河三角洲沉積物物質(zhì)來(lái)源主要為河流沉積,同時(shí)受到海洋潮流動(dòng)力影響;而且本研究發(fā)現(xiàn)黃河三角洲不同沉積區(qū)沉積動(dòng)力形式具有明顯差別,其中本研究沉積Ⅲ區(qū)與陳小英等[13]劃分的廢棄三角洲沉積區(qū)距離較近,兩項(xiàng)研究均揭示了該地區(qū)主要受海洋沉積作用的影響。此外,本研究中黃河與陳小英等[13]
研究的現(xiàn)行河口區(qū)位置一致,主要受河流作用以及沙嘴前方強(qiáng)潮流帶的作用??傊?,黃河三角洲沉積動(dòng)力同時(shí)受河流和海洋潮流影響,在黃河三角洲不同沉積區(qū)域沉積物主導(dǎo)的動(dòng)力因素有所不同。
黃河三角洲沉積物的粒度變化主要受到黃河影響,為河控型三角洲。1855年以來(lái),黃河頻繁改流導(dǎo)致黃河三角洲不同河流上沉積物特征均不相同。黃學(xué)勇等[12]認(rèn)為黃河徑流是現(xiàn)代黃河口南岸沉積物的主要物質(zhì)來(lái)源,黃河改流對(duì)黃河三角洲沉積環(huán)境產(chǎn)生重要影響。對(duì)比黃河三角洲不同位置的黃河河流沉積物粒度特征并結(jié)合黃河改流歷史發(fā)現(xiàn),較早時(shí)期廢棄的河流粒度較細(xì),而現(xiàn)黃河以及廢棄時(shí)間較晚的河流由于仍主要受現(xiàn)河流流水搬運(yùn)作用影響,河流流速較快,水動(dòng)力環(huán)境較強(qiáng),因而粒度更粗。在沾利河、草橋溝和挑河入??谔幱猩倭看憩F(xiàn)黃河的沉積物。這說(shuō)明黃河改流后,仍有少量黃河沉積物通過(guò)西部舊河流運(yùn)輸至入???。
近年來(lái),人類工程活動(dòng)、修建堤防、人工改流和開挖溝槽(尤其是河口疏浚)對(duì)黃河三角洲和其河口產(chǎn)生了巨大的影響。隨著河口向海延伸,人們?cè)邳S河三角洲陸地上種植樹木,用來(lái)防止海浪沖刷侵蝕,在油田附近建造防浪墻來(lái)防止海岸侵蝕[17]。在黃河三角洲北部,受修建堤壩、鞏固岸灘、石油開采和人工養(yǎng)殖業(yè)發(fā)展等人類活動(dòng)的干擾,人工岸線向海推進(jìn)。許多島嶼并入陸地,海灣被攔截,原來(lái)曲折的天然韻律型海岸被平直的人工海岸取代[21]。因此,對(duì)黃河三角洲地區(qū)進(jìn)行開發(fā)時(shí)應(yīng)注意生態(tài)防護(hù)與經(jīng)濟(jì)建設(shè)協(xié)調(diào)發(fā)展,以減少黃河三角洲北部岸線侵蝕對(duì)生態(tài)環(huán)境的不利影響。
4 結(jié)論與建議
1)現(xiàn)代黃河三角洲沉積物組成以粉砂(體積分?jǐn)?shù)71.10%)、砂(27.62%)為主,黏土(1.28%)較少。在馬新河、沾利河、草橋溝、挑河等大部分舊河流地區(qū),黏土體積分?jǐn)?shù)相對(duì)其他河道較高;在刁口河、神仙溝、現(xiàn)黃河以及部分沿海地區(qū),砂體積分?jǐn)?shù)較高。在馬新河、沾利河、草橋溝、挑河等距離海岸有一定距離的舊河流,粒級(jí)為細(xì)粉砂,分選較差,粒度分布曲線呈近對(duì)稱分布,峰態(tài)中等,概率累積曲線為兩段式,物質(zhì)搬運(yùn)方式主要為躍移和懸移;在刁口河、神仙溝、現(xiàn)黃河等東部河流,粒級(jí)為極細(xì)砂,分選中等,粒度分布曲線呈近對(duì)稱分布,峰態(tài)窄,概率累積曲線為二段式,物質(zhì)搬運(yùn)方式主要為躍移和懸移;在黃河三角洲北部沿海地帶,粒級(jí)為粗粉砂,分選較差,粒度分布曲線為正偏,峰態(tài)窄,概率累積曲線為四段式,其中躍移組分的斜率較高。
2)黃河三角洲表層沉積物由三個(gè)端元構(gòu)成,其中EM1為舊河流沉積動(dòng)力,代表舊河流在較弱的水動(dòng)力條件下逐漸沉積形成的河流沉積物,EM2為海洋潮流及波浪動(dòng)力作用下形成的沉積物,EM3主要來(lái)源為現(xiàn)黃河沉積物。
3)運(yùn)用聚類分析方法將黃河三角洲分為三類沉積環(huán)境:沉積Ⅰ區(qū)主要分布在黃河三角洲等西部河流的中上游,包括馬新河、沾利河、草橋溝、挑河,主要為舊河流沉積,水動(dòng)力較弱;沉積Ⅱ區(qū)主要分布在刁口河、神仙溝、現(xiàn)黃河(清水溝)等東部河流,現(xiàn)黃河河流搬運(yùn)作用較強(qiáng);沉積Ⅲ區(qū)主要分布在黃河三角洲北部沿岸及潮灘,受海洋潮流作用較強(qiáng),海岸受到侵蝕。
人類活動(dòng)對(duì)黃河三角洲和其河口產(chǎn)生了一定的影響,在此背景下,黃河三角洲岸線發(fā)生動(dòng)態(tài)改變。在黃河三角洲神仙溝以東部分岸線保存良好,但在刁口河以西岸線仍受到海洋侵蝕,應(yīng)加強(qiáng)岸線防護(hù)。
參考文獻(xiàn)(References):
[1] Zhang J, Huang W W, Shi M C. Huanghe (Yellow River) and Its Estuary: Sediment Origin, Transport and Deposition[J]. Journal of Hydrology, 1990, 120(1/2/3/4): 203223.
[2] 成國(guó)棟. 黃河三角洲現(xiàn)代沉積作用及模式[M]. 北京:地質(zhì)出版社, 1991.
Cheng Guodong. Modern Sedimentation and Model of Yellow River Delta [M]. Beijing: Geological Publishing House, 1991.
[3] 尹延鴻, 周永青, 丁東. 現(xiàn)代黃河三角洲海岸演化研究[J]. 海洋通報(bào), 2004, 23(2): 3240.
Yin Yanhong, Zhou Yongqing, Ding Dong. Evolultion of Modern Yellow River Delt Coast[J]. Marine Science Bulletin, 2004, 23(2): 3240.
[4] 薛允傳, 尹延鴻, 高抒. 黃河三角洲北部潮間帶沉積物的粒度特征[J]. 海洋科學(xué), 2001, 25(5): 5054.
Xue Yunchuan, Yin Yanhong, Gao Shu. Grain Size Characteristics of Tidal Flat Sediments, Northern Yellow River Delta[J]. Marine Sciences, 2001, 25(5): 5054.
[5] 李九發(fā), 李為華, 應(yīng)銘, 等. 黃河三角洲飛雁灘沉積物顆粒度分布和粒度參數(shù)特征及水動(dòng)力解釋[J]. 海洋通報(bào), 2006,25(3): 3844.
Li Jiufa, Li Weihua, Ying Ming, et al. Characteristic and Hydrodynamic Explanation of Distribution and Parameters of Sediment Granularity in Feiyan Shoal of Yellow River Delta[J]. Marine Science Bulletin, 2006, 25(3): 3844.
[6] 袁萍, 畢乃雙, 吳曉, 等. 現(xiàn)代黃河三角洲表層沉積物的空間分布特征[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2016, 36(2): 4957.
Yuan Ping, Bi Naishuang, Wu Xiao, et al. Surface Sediments at the Subaqueous Yellow River Delta: Classification and Distribution[J]. Marine Geology amp; Quaternary Geology, 2016, 36(2): 4957.
[7] 喬淑卿, 方習(xí)生, 石學(xué)法, 等. 黃河口及鄰近渤海海域表層沉積物中CaO和蒙皂石分布及其對(duì)黃河入海物質(zhì)運(yùn)移的指示[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2010, 30(1): 1723.
Qiao Shuqing, Fang Xisheng, Shi Xuefa, et al. Distribution of CaO and Smectite in Surface Sediments off the Yellow River Mouth and in the Nearby Bohai Sea and the Implications for Dispersion of the River Sediments to the Sea[J]. Marine Geology amp; Quaternary Geology, 2010, 30(1): 1723.
[8] 王昆山, 石學(xué)法, 蔡善武, 等. 黃河口及萊州灣表層沉積物中重礦物分布與來(lái)源[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2010, 30(6): 18.
Wang Kunshan, Shi Xuefa, Cai Shanwu, et al. Distribution and Provenance of the Surface Sediments of the Yellow River Mouth and Laizhou Bay Deduced from Heavy Minerals[J]. Marine Geology amp; Quaternary Geology, 2010, 30(6): 18.
[9] 李栓科. 近代黃河三角洲的沉積特征[J]. 地理研究, 1989, 8(4): 4555.
Li Shuanke. Sedimentary Characteristics of the Modern Yellow River Delta[J]. Geographical Research, 1989,8(4): 4555.
[10] 任韌希子. 黃河三角洲沉積特征與環(huán)境演變研究[D]. 上海: 華東師范大學(xué), 2012.
Ren Renxizi. Sedimentary Characteristics and Environment Evolution of the Yellow River Delta[D]. Shanghai: East China Normal University, 2012.
[11] 趙廣明, 葉青, 薛春汀, 等. 現(xiàn)代黃河三角洲陸上表層沉積物類型與沉積環(huán)境分區(qū)及岸線演變[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2013, 33(5): 4752.
Zhao Guangming, Ye Qing, Xue Chunting, et al. Surface Sediments, Sedimentary Subenvironments and Shoreline Evolution of Modern Yellow River Delta[J]. Marine Geology amp; Quaternary Geology, 2013, 33(5): 4752.
[12] 黃學(xué)勇, 張戈, 高茂生, 等. 現(xiàn)代黃河口南岸海洋沉積環(huán)境分析[J]. 海洋地質(zhì)前沿, 2019, 35(1): 1221.
Huang Xueyong, Zhang Ge, Gao Maosheng, et al. Analysis of Marine Sedimentary Environment on the South Bank of the Modern Yellow River Estuary[J]. Marine Geology Frontier, 2019, 35(1): 1221.
[13] 陳小英, 陳沈良, 劉勇勝. 黃河三角洲濱海區(qū)沉積物的分異特征與規(guī)律[J]. 沉積學(xué)報(bào), 2006, 24(5): 714721.
Chen Xiaoying, Chen Shenliang, Liu Yongsheng. Sediment Differentiation Along Nearshore Zone of the Yellow River Delta[J]. Acta Sedimentologica Sinica, 2006, 24(5): 714721.
[14] 王一鳴. 黃河水沙變化及其對(duì)黃河三角洲變化的影響分析[D]. 北京:中國(guó)科學(xué)院大學(xué), 2018.
Wang Yiming. Analysis on the Water and Sediment Variation in Yellow River and Its Influence on the Change of Yellow River Delta[D]. Beijing:University of Chinese Academy of Sciences, 2018.
[15] 樊輝. 黃河口泥沙輸移及三角洲的近期演變[D].青島:中國(guó)科學(xué)院研究生院, 2005.
Fan Hui. Sediment Transport in the Huanghe River Mouth and the Recent Deltaic Evolution [D]. Qingdao: Graduate School of the Chinese Academy of Sciences, 2005.
[16] 劉鋒. 黃河口及其鄰近海域泥沙輸運(yùn)及其動(dòng)力地貌過(guò)程[D]. 上海:華東師范大學(xué), 2012.
Liu Feng. Sediment Transport and Dynamic Geomorphology Process in the Yellow River Estuary and Its Adjacent Sea[D]. Shanghai: East China Normal University, 2012.
[17] 王兆印, 程?hào)|升, 劉成. 人類活動(dòng)對(duì)典型三角洲演變的影響:Ⅱ:黃河和海河三角洲[J]. 泥沙研究,2006,25(1): 7680.
Wang Zhaoyin, Cheng Dongsheng, Liu Cheng. Delta Processes and Management Strategies in China:Ⅱ :The Yellow and Haihe River Deltas[J]. Journal of Sediment Research, 2006,25(1): 7680.
[18] 繆紅兵, 喬璐璐, 仲毅, 等. 人類活動(dòng)和自然演變共同驅(qū)動(dòng)下黃河三角洲海域潮波及物質(zhì)輸運(yùn)變化[J]. 海洋學(xué)報(bào), 2022, 44(9): 7386.
Miao Hongbing, Qiao Lulu, Zhong Yi, et al. Evolution of Tidal System and Material Transport off the Huanghe River Delta Induced by Human Activities and Natural Evolution[J]. Haiyang Xuebao, 2022, 44(9): 7386.
[19] 楊偉. 現(xiàn)代黃河三角洲海岸線變遷及灘涂演化[J].海洋地質(zhì)前沿, 2012, 28(7): 1723.
Yang Wei. Shifting of Coastline and Evolution of Tidal Flat in Modern Yellow River Delta[J]. Marine Geology Frontiers, 2012, 28(7): 1723.
[20] 黃波, 郭慧敏, 孫吉?jiǎng)偅?等. 黃河三角洲刁口河海岸侵蝕時(shí)空演變與防護(hù)對(duì)策[M]. 鄭州: 黃河水利出版社, 2020.
Huang Bo,Guo Huimin, Sun Jigang, et al. Spatial-Temporal Evolution of Coastal Erosion and Protection Countermeasures of the Diaokou River in the Yellow River Delta[M]. Zhengzhou: Yellow River Water Conservancy Press, 2020.
[21] 黃海軍, 樊輝. 1976年黃河改流以來(lái)三角洲近岸區(qū)變化遙感監(jiān)測(cè)[J]. 海洋與湖沼, 2004, 35(4): 306314.
Huang Haijun, Fan Hui. Monitoring Changes of Nearshore Zones in the Huanghe (Yellow River) Delta Since 1976[J]. Oceanologia et Limnologia Sinica, 2004, 35(4): 306314.
[22] 安樂生, 劉貫群, 葉思源, 等. 黃河三角洲濱海濕地健康條件評(píng)價(jià)[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2011, 41(4): 11571165.
An Lesheng, Liu Guanqun, Ye Siyuan, et al. Health Assessment of Coastal Wetlands in the Yellow River Delta[J]. Journal of Jilin University(Earth Science Edition), 2011, 41(4): 11571165.
[23] 成都地質(zhì)學(xué)院陜北隊(duì). 沉積巖(物)粒度分析及其應(yīng)用[M]. 北京: 地質(zhì)出版社,1978.
Shanbei Team,
Chengdu College of Geology. Sedimentary Rock (Material) Particle-Size Analysis and Its Application[M]. Beijing: Geological Publishing House, 1978.
[24] Weltje G J. End-Member Modeling of Compositional Data: Numerical-Statistical Algorithms for Solving the Explicit Mixing Problem[J]. Mathematical Geology, 1997, 29(4): 503549.
[25] 白敏, 魯瑞潔, 丁之勇, 等.青海湖湖東沙地粒度端元分析及其指示意義[J]. 第四紀(jì)研究, 2020, 40(5): 12031215.
Bai Min, Lu Ruijie, Ding Zhiyong, et al. End-Member Analysis of GrainSize in the East of Qinghai Lake and Its Environmental Implications[J]. Quaternary Sciences, 2020, 40(5): 12031215.
[26] Paterson G A, Heslop D. New Methods for Unmixing Sediment Grain Size Data[J]. Technical Reports : Methods, 2015, 16(12): 44944506.
[27] 張曉東, 翟世奎, 許淑梅. 端元分析模型在長(zhǎng)江口鄰近海域沉積物粒度數(shù)據(jù)反演方面的應(yīng)用[J]. 海洋學(xué)報(bào), 2006, 28(4): 159166.
Zhang Xiaodong, Zhai Shikui, Xu Shumei. The Application of Grain-Size End-Member Modeling to the Shelf Near the Estuary of Changjiang River in China[J]. Haiyang Xuebao, 2006, 28(4): 159166.
[28] 林鎮(zhèn)坤, 王愛軍, 葉翔. 南流江河口水下三角洲表層沉積物端元分析及其沉積動(dòng)力環(huán)境意義[J]. 沉積學(xué)報(bào), 2019, 37(1): 124134.
Lin Zhenkun, Wang Aijun, Ye Xiang. End-Member Analysis for Surficial Sediment of Nanliujiang River Subaqueous Delta and Associated Sediment Dynamic Environmental Significance[J]. Acta Geologica Sinica, 2019, 37(1): 124134.
[29] 薛成鳳, 賈建軍, 高抒, 等. 中小河流對(duì)長(zhǎng)江水下三角洲遠(yuǎn)端泥沉積的貢獻(xiàn):以椒江和甌江為例[J]. 海洋學(xué)報(bào), 2018, 40(5): 7589.
Xue Chengfeng, Jia Jianjun, Gao Shu, et al. The Contribution of Middle and Small Rivers to the Distal Mud of Subaqueous Changjiang Delta: Results from Jiaojiang River and Oujiang River[J]. Haiyang Xuebao, 2018, 40(5): 7589.
[30] 杜新強(qiáng),王鈺升,冶雪艷,等. 大清河流域平原區(qū)地下水人工補(bǔ)給潛力與補(bǔ)給方式分析[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2022, 52(2):535549.
Du Xinqiang, Wang Yusheng, Ye Xueyan,et al. Artificial Aquifer Recharge Potential and Methods in Plain Area of Daqing River Basin[J]. Journal of Jilin University (Earth Science Edition), 2022, 52(2):535549.
[31] 洪梅,趙明明,魏濤,等.基于SWAT模型的翔陽(yáng)河流域斷流河道復(fù)流研究[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2024,54(1);253263.
Hong Mei,Zhao Mingming,Wei Tao,et al. Study on Water Restoration in the Dried-up River of the Fuyang River Basin Based on the SWAT Model[J]. Journal of Jilin University (Earth Science Edition),2024,54( 1):253263.
[32] 張盼, 吳建政, 胡日軍, 等. 萊州灣西南部表層沉積物粒度分布特征及其現(xiàn)代沉積環(huán)境分區(qū)[J]. 海洋地質(zhì)前沿, 2014, 30(9): 1117.
Zhang Pan, Wu Jianzheng, Hu Rijun, et al. Distribution of Surface Sediments and Modern Sedimentary Environment in the Southwestern Laizhou Bay[J]. Marine Geology Frontiers, 2014, 30(9): 1117.
[33] 楊陽(yáng), 高抒, 周亮, 等. 海南新村港潟湖表層沉積物粒度特征及其沉積環(huán)境[J]. 海洋學(xué)報(bào), 2016, 38(1): 94105.
Yang Yang, Gao Shu, Zhou Liang, et al. Grain Size Distribution of Surface Sediments and Sedimentary Environment in the Lagoon of Xincun, Hainan Island[J]. Haiyang Xuebao, 2016, 38(1): 94105.
[34] 邢尚軍, 郗金標(biāo), 張建鋒, 等. 黃河三角洲植被基本特征及其主要類型[J]. 東北林業(yè)大學(xué)學(xué)報(bào), 2003, 31(6): 8586.
Xing Shangjun, Xi Jinbiao, Zhang Jianfeng, et al. Basic Characteristics and Main Types of Vegetation in the Yellow River Delta[J]. Journal of Northeast Forestry University, 2003, 31(6): 8586.
[35] 東營(yíng)市人民政府. 自然環(huán)境 [EB/OL]. (20230514).http://www.dongying.gov.cn/col/col307764/index.html.
Dongying Municipal People’s Government. Natural Environment [EB/OL]. (20230514)." http://www.dongying.gov.cn/col/col307764/index.html.
[36] Folk R L, Ward W C. Brazos River Bar (Texas): A Study in the Significance of Grain Size Parameters[J]. Journal of Sedimentary Research, 1957, 27(1): 326.
[37] Wentworth C K. A Scale of Grade and Class Terms for Clastic Sediments[J]. Journal of Geology, 1922, 30(5): 377392.