基于地理探測器的黃海沉積物中值粒徑分布影響因子分析

劉長達，唐秋華*，喬淑卿，周興華

(1.自然資源部 第一海洋研究所,山東 青島 266061；2.山東科技大學(xué) 海洋科學(xué)與工程學(xué)院,山東 青島 266590)

海底沉積物粒度是海洋研究的重要參數(shù),其特征以及空間分布反映了沉積物輸運、沉降的過程,影響海洋物理探測的精度,是進行海洋科學(xué)調(diào)查和研究的重要內(nèi)容。20世紀90年代以來,中外科學(xué)家對中國近海進行了大量科學(xué)調(diào)查和研究,在該區(qū)域積累了詳實的資料。張曉東等[1]采用沉積物粒度端元分析方法,利用南黃海表層沉積物粒度數(shù)據(jù)分離出了4個端元,結(jié)合研究區(qū)沉積動力的格局,給出了評價沉積物粒度對物源和沉積動力環(huán)境響應(yīng)的指標。郭曉陽等[2]采用專業(yè)分析的方法對季風(fēng)邊緣區(qū)內(nèi)陸湖泊的沉積物粒度組分分布特征進行了研究,并按組分分析了影響其特征的形成因素,解釋了水動力條件和海浪二次懸移再沉積等作用在沉積過程中的重要影響。喬淑卿等[3]使用Gao-Collins改進的徑粒趨勢分析方法研究了渤海地區(qū)沉積物粒度特征以及輸運趨勢,解釋了潮流和海洋環(huán)流對沉積物粒度分布格局的控制作用。王華新和線薇微[4]通過分析4個季節(jié)中長江口表層沉積物粒度與沉積物有機碳的分布情況,解釋了沉積物對有機碳的富集作用,證明了沉積物有機碳含量可以指示沉積物粒度。張晉等[5]使用粒度參數(shù)計算和Q 型聚類分析的方法探討了南海南部沉積物粒度分布特征及影響因素,得出了該區(qū)域沉積物粒度主要與物源、水動力條件、水深和地形及洋流等因素有關(guān)的結(jié)論。Alexander等[6]采用同位素方法分析認為,黃河攜帶的泥沙對黃海沉積物粒度的空間分布起著重要影響。目前,針對沉積物粒度空間分布及其影響因素的研究中還存在一些尚未解決的重要問題,如:①研究多局限于水動力環(huán)境因素,未考慮其他因素的影響；②一些研究僅簡單定性地描述了影響沉積物粒度分布的原因,未定量描述因素的影響程度；③現(xiàn)行對于沉積物粒度空間分布的研究多局限于單一影響因素,未綜合考慮不同因素之間的交互作用。

地理探測器可以獲取不同要素間空間分層異質(zhì)性(空間分異性)的顯著程度參數(shù),也可以定量分析要素與某一因子空間分布的相關(guān)性,根據(jù)相關(guān)性結(jié)果分析因子的影響力[7-10]。該方法現(xiàn)已被廣泛應(yīng)用于城鎮(zhèn)擴張、土地利用以及脫貧攻堅等領(lǐng)域,由地理探測器分析結(jié)果可以定量評估并比較各因子的重要性,為專業(yè)決策提供科學(xué)支撐[11-13]。

沉積是一個復(fù)雜過程,其影響因子較多。沉積區(qū)域的水動力環(huán)境、水濁度、化學(xué)因素等都有可能對沉積作用產(chǎn)生影響。沉積物中值粒徑作為表示粒度大小的典型值,可以深刻反映環(huán)境等對沉積作用的影響?；诂F(xiàn)行的沉積物粒徑分布及其影響因素研究中存在的問題,本文應(yīng)用地理探測器的地理學(xué)新方法,探究與黃海沉積物中值粒徑分布有關(guān)的影響因子,旨在找尋與沉積物中值粒徑空間分布有關(guān)的因子,并定量分析不同單因子、多因子對于沉積物空間分布的影響程度,為今后探究沉積物粒度的分布及其影響因素的研究提供新思路。

1 方法與數(shù)據(jù)源

1.1 研究方法

空間分異性(Spatial Stratified Heterogeneity,SSH)指某一屬性值在不同類型或區(qū)域之間存在差異的現(xiàn)象。該性質(zhì)是由自然界中多種因素共同作用的結(jié)果。探究空間分異性背后地規(guī)律,可以更深層次地理解形成該現(xiàn)象的機理。地理探測器正是識別、檢驗、尋找和歸因空間分異的工具,具體包括因子、風(fēng)險、生態(tài)、交互四個探測器。

①因子(X)探測可以定量分析因子對于屬性(Y)的解釋程度,即因子解釋力。解釋力大小通過統(tǒng)計值q(0≤q≤1)度量[7-10]。計算方法:

式中:h為某一個分層；N和N h分別為研究區(qū)總單元數(shù)和層h內(nèi)的單元數(shù)；σ2和分別為樣本Y值在研究區(qū)的總離散方差和層h內(nèi)的離散方差；STS為研究區(qū)總方差；SWS為研究區(qū)所有層內(nèi)方差之和；L為總分層數(shù)。q值越大說明因子X對屬性Y的解釋能力越強。另外,可通過非中心F分布(non-centralF-distribution)來檢驗q值是否統(tǒng)計顯著。

②風(fēng)險探測可以用來判斷2個因子區(qū)域間的屬性均值是否有顯著差異。根據(jù)因子分層原則,可進一步分析屬性均值在不同層間的變化規(guī)律,進一步分析屬性值與因子的關(guān)系。

③生態(tài)探測多用于比較2個因子對屬性值空間分布的影響是否有顯著差異。

④交互探測是地理探測器的關(guān)鍵部分,該部分通過不同要素之間的疊加分析,可以識別不同風(fēng)險因子之間的交互作用。即評估2個因子共同作用時對屬性值的解釋能力,并與因子探測比較解釋力的變化,探究其增加或減弱情況。

1.2 數(shù)據(jù)源

本文通過科考現(xiàn)場海洋采樣以及其他觀測手段得到數(shù)據(jù),主要包括:

①海底沉積物數(shù)據(jù),為“2007—2009年中國科學(xué)院海洋研究所4次黃海科考”所得數(shù)據(jù)[15],該數(shù)據(jù)有明確的采樣記錄,經(jīng)過了嚴格且細致的數(shù)據(jù)篩選?；谠摂?shù)據(jù)可對黃海沉積物情況有一個細致且全面的了解(圖1)。

圖1 黃海海底地形以及采樣點分布Fig.1 Topography of the study area and sampling locations

②海底底層水?dāng)?shù)據(jù),來自于混合坐標海洋模式(Hybrid Coordinate Ocean Model,HYCOM)[16-17],該數(shù)據(jù)集是近些年來主流的全球海洋環(huán)流模式,其再分析數(shù)據(jù)可以得到表層、混合層以及底層水的流速。根據(jù)此數(shù)據(jù)可以派生出底層水矢量流速以及流向數(shù)據(jù)。

③海底底層水化學(xué)指標數(shù)據(jù),來自于CARS2009,該數(shù)據(jù)集為現(xiàn)代海洋測量期間的平均海洋屬性的網(wǎng)格場,且只提供插值或平均之后的數(shù)據(jù),不提供任何特定年份的數(shù)據(jù)[18]。

④海底地形數(shù)據(jù)以及與地形有關(guān)因子,采用全球海底地形(General Bathymetric Chart of the Oceans,GEBCO)數(shù)據(jù)[19],該數(shù)據(jù)是由國際海道測量組織和政府間海洋委員會聯(lián)合發(fā)布的最全面的世界大洋海底地形數(shù)據(jù),包括半分間隔數(shù)據(jù)、一分間隔數(shù)據(jù)及其相關(guān)輔助數(shù)據(jù)[20-21]。

⑤水濁度數(shù)據(jù),來自于美國國家航空航天局(National Aeronautics and Space Administration,NASA)水色遙感實驗室的MODIS-Aqua業(yè)務(wù)化三級遙感產(chǎn)品[22]。本研究采用2007—2009年衛(wèi)星490 nm 波長的海水下行輻射通量的漫衰減系數(shù)代表水濁度。由于490 nm 漫衰減系數(shù)與水體濁度對應(yīng)良好,業(yè)界多依此反演水體濁度。胡靜雯等[23]通過現(xiàn)場科考采樣數(shù)據(jù)的比對,印證了MODIS-Aqua遙感數(shù)據(jù)的在東海區(qū)域反演水濁度的可靠性。

1.3 研究區(qū)域概況

黃海(119°～126°E,32°～39°N)北部鄰接渤海,東西部被中國大陸和朝鮮半島包圍,南部與東海相連,是一個半封閉海。其位于東亞大陸架之上,水深較淺,平均水深約44 m,最大水深約140 m。一般認為山東半島成山角與朝鮮長山串角之間的連接線為南北黃海分界線?？傮w來看,黃海海底地勢由北、西、東三面向黃海中央及東南方向傾斜,但坡度不大,總體比較平坦。在黃海中部有一個明顯南北走向的水下洼地,水深約60～80 m[24-25]。

2 地理探測器分析

2.1 因子選擇及預(yù)處理

鑒于地理探測器對于研究問題所需要的假設(shè)條件較少、對因子間共線性免疫,且可以通過查表或相關(guān)軟件進行顯著性檢驗等,在使用地理探測器進行因子分析時可以大膽選取較多因子進行分析,如不滿足顯著性則舍棄該因子即可[26]。

結(jié)合現(xiàn)有的樣本數(shù)據(jù)以及專業(yè)知識,充分選取可能影響沉積物中值粒徑分布的5類因子進行探測研究(表1)。第一,水動力因子。地質(zhì)學(xué)家觀點認為,水動力是影響沉積作用的主要因素。本文選取底層水東向流速(X1)、北向流速(X2)以及派生出的矢量流速(X3)、流向(X4)作為4個水動力因子。第二,化學(xué)因子。普斯托瓦洛夫沉積化學(xué)理論表明:在沉積物搬運和沉積過程中,不同化學(xué)元素的溶解物質(zhì)對應(yīng)不同的沉積順序,會產(chǎn)生化學(xué)沉積分異作用[27]。另外,沉積物顆粒會對水中游離有機物質(zhì)產(chǎn)生吸附作用[28]。本研究主要選取底層水硅酸鹽(X5)、硝酸鹽(X6)、磷酸鹽(X7)、溶解氧(X8)以及底層水鹽度(X9)和溫度(X10)六個化學(xué)因子。第三,生物因素。生物通過新陳代謝和搬運作用會產(chǎn)生某種元素的富集,從而影響沉積作用。本文選取指示生物指標的葉綠素濃度(X11)和沉積物有機碳含量(X12)兩個因子進行分析。第四,地形等因子?？紤]到地形起伏和變化是影響水動力環(huán)境的重要因素,水深(X13)、坡向(X14)、坡度(X15)等的變化必然在某種程度上影響沉積作用。第五,反映水體濁度的因子。本文主要選取490 nm 漫衰減系數(shù)(X16)和沉積物含水量(X17)因子進行分析。圖2為部分因子變量結(jié)果圖。

圖2 部分因子變量圖Fig.2 Variable plots for partial factors

表1 沉積物粒度分異因子探測結(jié)果Table 1 Results of differentiation factors for sediment grain size

地理探測器中為了保證自變量是類型量,需要對連續(xù)變量進行離散化分層[7-10]。為了得到較優(yōu)的離散分層結(jié)果,本文采用專家經(jīng)驗、自然斷點相結(jié)合的方法進行分層。其中,與方向有關(guān)的因子按八方位知識進行分層；其余變量均采用自然斷點法分為5層。

2.2 探測結(jié)果及分析

2.2.1 因子探測

因子探測器揭示了各環(huán)境因子對沉積物中值粒徑影響程度的差異[29]。從地理探測器模型總體來看,通過顯著性檢驗的因子并不多,依次為:底層水矢量流速(以下簡稱底層水流速X3)、流向(X4)、沉積物有機碳含量(X12)以及坡向(X14)、490 nm 漫衰減系數(shù)(X16)、沉積物含水量(X17)等。顯著因子的解釋能力也有較大差別。

水動力因子中底層水流速和流向的因子解釋力(q統(tǒng)計值)為0.076和0.244,在一定程度上印證了水動力因素對沉積過程的重要影響。坡向、490 nm 漫衰減系數(shù)因子解釋力為0.124和0.199,這說明了沉積物中值粒徑不僅與水動力因素有關(guān),還與海底地形、水體渾濁等因素有關(guān)系。另外,指示生物因素的沉積物有機碳含量(X12)通過了顯著檢驗,其解釋力為0.137。沉積物含水量對沉積物中值粒徑有著較強的解釋能力,其因子解釋力為0.408。

造成這種結(jié)果的原因:首先,水動力因素是影響沉積物沉積作用的重要因素,不同粒沉積物在海底的搬運、沉積方式不同；其次,水動力因素與地形有著一定的關(guān)系,復(fù)雜的地形處通常對應(yīng)復(fù)雜的水動力因素；第三,中值粒徑主要反映沉積物水動力環(huán)境以及物源等因素,對于生物搬運沉積和化學(xué)沉積分異作用的指示性不敏感；第四,沉積物顆粒對有機碳含量有一定的富集作用,且該作用與中值粒徑大小有關(guān)；最后,490 nm 漫衰減系數(shù)可對應(yīng)水體渾濁的程度,可以反映出細沉積物懸浮搬運的過程。

2.2.2 風(fēng)險探測

風(fēng)險探測器揭示了各環(huán)境因子內(nèi)部不同類別分區(qū)間的顯著性差異,用于判斷2個子分層間的屬性均值是否具有顯著性差別[30]。通過因子探測顯著性檢驗的因子輸入Geodetector模型的風(fēng)險探測結(jié)果如圖3所示。結(jié)合因子選取以及分層原則,即對于底層水流速(X3)、沉積物有機碳含量(X12)、490 nm 漫衰減系數(shù)(X16)、沉積物含水量(X17)等分為5層的因子,其值均隨著子分層數(shù)的增加而增加。對于流向(X4)和坡向(X14)等因子,具體分層情況見表2。

圖3 沉積物粒度分異風(fēng)險探測結(jié)果Fig.3 Risk detection results of differentiation of sediment grain size

表2 因子分層結(jié)果Table 2 Results of factor stratification

具體而言:沉積物中值粒徑均值隨著底層水流速(X3)的增加呈現(xiàn)先變大后減小的趨勢,在流速為3～4 cm/s時出現(xiàn)粒徑均值的最大值(40.887μm)。沉積物中值粒徑均值與沉積物有機碳含量(X12)、沉積物含水量(X17)均呈現(xiàn)負相關(guān)的關(guān)系,即粒徑均值越大,樣本中有機碳和含水量越低,具體而言,含水量最多時(含水量＞0.55%),粒徑均值最小(7.578μm)；有機碳含量最多時(有機碳含量＞0.95%),粒徑均值最小(13.677 μm),這與王華新和線薇微[4]的研究結(jié)果一致。沉積物中值粒均值隨著490 nm 漫衰減系數(shù)(X17)的增加呈現(xiàn)先變大后減小的趨勢,粒徑均值最大值(44.264μm)出現(xiàn)在490 nm 漫衰減系數(shù)的第4個子分層。最后,沉積物中值粒徑均值在不同坡向、不同流向時有不同粒徑均值,在坡向(X14)為西時粒徑均值最大(50.818 μm),在底層水流向(X4)為北時粒徑均值最大(72.174μm)。另外,以底層水流向分層以及每個子分層對應(yīng)的底層水平均流速進行線性回歸(圖4),兩者回歸性較好,在黃海區(qū)域底層水較高流速往往具有特定方向,流速與流向并不是相互獨立的。

圖4 底層水平均流速與流向關(guān)系Fig.4 The relationship between bottom horizontal mean velocity and flow direction

分析造成這種結(jié)果的原因:首先,海底碎屑顆粒物質(zhì)的搬運、移動和沉積過程與其對應(yīng)的起動流速有關(guān),沉積物中值粒徑的變化可以反映出相應(yīng)水動力環(huán)境的變化；其次,黃海區(qū)域底層水高速水流對應(yīng)著特定方向,流速與流向共同形成了水動力參數(shù),進而對沉積作用產(chǎn)生影響；第三,水體中顆粒的沉降受限于顆粒半徑以及水體流速,其中粒徑小于100μm 的碎屑細顆粒沉降速度與顆粒的半徑的平方成正比,粗顆粒沉降速度則與半徑的平方根成正比,并且水流速度遠大于沉降速度時會產(chǎn)生懸浮作用,表現(xiàn)為水中細顆粒增加,濁度上升；最后,細顆粒沉積物對有機碳有著富集作用,因此粒徑均值較小部分對應(yīng)著沉積物有機碳含量的高值。

2.2.3 交互探測

地理探測器的優(yōu)勢之一是可以探測2種因子之間的交互作用,通過因子之間的疊加,分析疊加之后的解釋力變化并判斷出因子之間的交互作用類型[31-32]。將通過顯著性檢驗的因子輸入Geodetector模型,其交互結(jié)果如表3所示。

由表3可見,任意2種因子進行交互時,解釋力均有了一定提高,全部大于對應(yīng)的2種因子的解釋力之和,因此任意2種因子均具有非線性交互增強作用。交互作用中對其他因子解釋力提升最顯著的因子為底層水流速。具體看,因子交互之后,解釋力多在0.24～0.68的范圍內(nèi)。底層水流速與底層水流向的交互之后,解釋力提高0.668,這與傳統(tǒng)研究結(jié)果一致[1-4,6]。造成這種結(jié)果的原因在于:首先,底層水流速是水動力中最主要的因素,沉積物在不同流動強度的水中會產(chǎn)生不同的沉積結(jié)果；其次,底層水流速、流向共同描述了水動力因素,水動力深刻影響著沉積物的沉降過程,對沉積物的物質(zhì)組成、粒徑大小和分選效果都有一定影響。

表3 沉積物粒度分異交互探測結(jié)果Table 3 Interactive detection results of differentiation of sediment grain size

總之,沉積物沉積是一個復(fù)雜的過程,水動力因素在其中扮演著重要角色,但沉積作用不僅與水動力因素相關(guān),海底地形、490 nm 漫衰減系數(shù)、沉積物的物質(zhì)組成等不同因子都不同程度地影響著沉積作用。

3 結(jié)論

本文分析了影響黃海海底沉積物中值粒徑空間分布的因子,并探究了各因子與沉積物中值粒徑的對應(yīng)關(guān)系,定量刻畫了各因子的影響力程度,得出結(jié)論如下:

①影響黃海沉積物中值粒徑空間分布的主要因子有:底層水流速、底層水流向、沉積物有機碳含量、坡向、490 nm 漫衰減系數(shù)和沉積物含水量。

②沉積物中值粒徑均值與沉積物含水量、有機碳含量之間呈負相關(guān)；不同坡向、流向?qū)?yīng)著不同的粒徑均值；沉積物在底層水流速、漫衰減系數(shù)的不同子分層中有不同的粒徑均值。另外,黃海區(qū)域的底層水流向與流速之間有著一定的對應(yīng)關(guān)系。

③交互作用均為非線性交互增強。進行交互作用時底層水流速影響力最大,底層水流向和底層水流速的解釋力最高。

④進行因子探測時解釋力最強的2個因子分別為沉積物含水量(0.408)、底層水流向(0.244)；進行交互探測時解釋力最強的2個因子分別為底層水流速與流向(0.668)、底層水流向與沉積物含水量(0.563)。