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

劉長(zhǎng)達(dá)，唐秋華*，喬淑卿，周興華

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

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

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

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

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

1.1 研究方法

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

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

式中:h為某一個(gè)分層；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對(duì)屬性Y的解釋能力越強(qiáng)。另外,可通過非中心F分布(non-centralF-distribution)來檢驗(yàn)q值是否統(tǒng)計(jì)顯著。

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

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

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

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

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

①海底沉積物數(shù)據(jù),為“2007—2009年中國(guó)科學(xué)院海洋研究所4次黃?？瓶肌彼脭?shù)據(jù)[15],該數(shù)據(jù)有明確的采樣記錄,經(jīng)過了嚴(yán)格且細(xì)致的數(shù)據(jù)篩選?；谠摂?shù)據(jù)可對(duì)黃海沉積物情況有一個(gè)細(xì)致且全面的了解(圖1)。

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

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

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

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

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

1.3 研究區(qū)域概況

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

2 地理探測(cè)器分析

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

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

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

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

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

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

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

2.2.1 因子探測(cè)

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

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

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

2.2.2 風(fēng)險(xiǎn)探測(cè)

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

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

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

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

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

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

2.2.3 交互探測(cè)

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

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

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

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

3 結(jié)論

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

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

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

③交互作用均為非線性交互增強(qiáng)。進(jìn)行交互作用時(shí)底層水流速影響力最大,底層水流向和底層水流速的解釋力最高。

④進(jìn)行因子探測(cè)時(shí)解釋力最強(qiáng)的2個(gè)因子分別為沉積物含水量(0.408)、底層水流向(0.244)；進(jìn)行交互探測(cè)時(shí)解釋力最強(qiáng)的2個(gè)因子分別為底層水流速與流向(0.668)、底層水流向與沉積物含水量(0.563)。