馬永強，馬 超，歐陽芳，杭 佳，石 云

(1寧夏大學 a農(nóng)學院，b資源環(huán)境學院，寧夏 銀川 750021；2 寧夏回族自治區(qū)自然資源廳 土地勘測規(guī)劃院，寧夏 銀川 750000； 3中國科學院動物所 農(nóng)業(yè)蟲害鼠害綜合治理研究國家重點實驗室，北京 100101)

黃土丘陵溝壑區(qū)是中國黃土高原的重要組成部分，該區(qū)域水土流失嚴重，生態(tài)脆弱。近年來實施的生態(tài)恢復舉措使該區(qū)生態(tài)系統(tǒng)功能得到改善，生物多樣性明顯提高。人為設計和改變區(qū)域環(huán)境因子會影響生物群落多樣性的空間分布特征和演變規(guī)律，為此探究群落多樣性與環(huán)境因子變化的關系，維持生態(tài)系統(tǒng)穩(wěn)定性已成為重要的研究方向[1]。

昆蟲作為世界上種類和數(shù)量最多的生物類群，其群落多樣性受局部環(huán)境條件和周圍景觀特征的影響[2-6]。特別是生物環(huán)境中較為多樣化的鞘翅目昆蟲類群，其食性復雜[7]，擴散力較低，易受到地面環(huán)境影響[8-9]，因此適合作為研究物種空間分布規(guī)律的對象。地形、氣候、生物等環(huán)境因子是引起生態(tài)環(huán)境變化的重要因素，這些因子綜合影響物種生境所需光、熱、水等要素的空間分布，同時影響動植物群落分布及種群格局的形成[10-11]。長期以來，國內外學者在研究昆蟲群落與地形關系方面，主要建立在感性認知和定性分析的基礎上[12-13]，定量分析研究較為少見。采用定量分析方法研究區(qū)域地表甲蟲物種多樣性和空間分布規(guī)律，探尋其變化與不同環(huán)境因子的關系，進而可對研究區(qū)生態(tài)效益進行評價。

數(shù)理統(tǒng)計和地理信息系統(tǒng)(GIS)空間分析方法是分析物種分布及其影響因素的重要研究手段。將統(tǒng)計技術和GIS技術引入研究，使量化后的昆蟲數(shù)據(jù)與GIS分析數(shù)據(jù)結合，可定量分析昆蟲群落分布與環(huán)境因子之間的關系。廣義可加模型(GAM)不需要預定義參數(shù)模型[14-16]，數(shù)據(jù)擬合更具靈活性，能應用非參數(shù)的方法檢測數(shù)據(jù)結構，找出其中規(guī)律，得到更好的預測結果[17]。目前，國內學者已將GAM模型應用于海洋漁業(yè)、湖泊微生物[18]、潛在植被預測以及草原昆蟲格局空間關系[19]等方面。但在黃土丘陵溝壑區(qū)僅停留在定性分析的基礎上[20]。

黃土丘陵溝壑區(qū)地形復雜且具有強烈的褶皺隆起，易影響區(qū)域內動植物分布。小流域是黃土丘陵溝壑區(qū)最小的自然集水單元。以小流域作為水土保持治理最小單元，有利于更加精確高效治理生態(tài)脆弱區(qū)域，實現(xiàn)生態(tài)環(huán)境治理和打造清潔型小流域雙目標。研究小流域物種分布與環(huán)境因子之間的關系，有助于進一步了解區(qū)域內生態(tài)效益變化情況。本研究以適應性較強、活動范圍大的鞘翅目地表甲蟲為研究對象，以地處黃土丘陵溝壑區(qū)高建堡典型小流域為研究區(qū)，基于統(tǒng)計方法和地理信息系統(tǒng)(GIS)技術，使量化昆蟲數(shù)據(jù)與GIS分析數(shù)據(jù)相結合，構建GAM，定量分析昆蟲群落與環(huán)境因子的關系，探究其最佳生境，促進生態(tài)脆弱區(qū)域生物多樣性保護及生態(tài)環(huán)境恢復。

1 材料與方法

1.1 研究區(qū)概況

高建堡小流域位于寧夏回族自治區(qū)東南部、六盤山東麓彭陽縣王洼鎮(zhèn)(36°3′44″-36°6′53″N,106°38′15″-106°40′43″E)，總面積1 076.89 hm2，海拔1 297～2 483 m。流域地處黃土高原，具有黃土丘陵溝壑區(qū)第二副區(qū)典型特征，梁、峁、溝壑并存，溝道發(fā)育，農(nóng)業(yè)綜合區(qū)劃屬于農(nóng)牧區(qū)，多年平均氣溫6.8 ℃，年平均降水量420 mm，年平均蒸發(fā)量為1 550 mm，蒸發(fā)強烈，屬半干旱大陸性季風氣候。土壤以黃綿土為主，主要分布在梁峁狀丘陵坡地上，土壤肥力中下，生產(chǎn)力水平低，是寧夏水土流失最嚴重的縣域之一，也是全國重點水土保持治理區(qū)。高建堡小流域面積較小，地形破碎，對地表甲蟲物種多樣性和空間分布有重要的影響。

1.2 采樣時間與樣地設置

地表甲蟲的取樣采用國際通用的陷阱法[21-22]。用一次性塑料水杯(高9 cm，上口直徑7.5 cm，下口直徑4 cm)作為陷阱容器，引誘劑為醋、糖、醫(yī)用酒精和水混合液(質量比為2∶1∶1∶20)，每個陷阱器放引誘劑40～60 mL。

在高建堡流域選擇6個植物群落分布典型且坡向、坡度和海拔特征分明的地段作為地表甲蟲調查區(qū)域，具體信息見表1。每個調查區(qū)域中設置3個20 m×20 m的樣方，樣方間隔20 m以上，在每個樣方的四角各放置1個陷阱器，每個調查區(qū)域設12個陷阱器，共72個陷阱器，分別收集每個陷阱器的捕獲結果，統(tǒng)計并編號。利用GPS(eXplorist 500 LE,MAGELLAN)確定每個陷阱的精準位置。2018年7月10日-8月25日，每隔5 d對72個陷阱器收集1次，并更換誘液。用體積分數(shù)75%酒精保存帶回實驗室鑒定到種，并統(tǒng)計數(shù)量。

表1 黃土丘陵區(qū)寧夏高建堡流域典型植物群落類型Table 1 Typical plant community types in Gaojianbao watershed of Ningxia in the Loess Hilly Region

1.3 研究方法

廣義可加模型(GAM)是廣義線性模型的半?yún)?shù)擴展[23]，無需定義參數(shù)模型。該模型具有以下優(yōu)點：(1)不受函數(shù)驅動，按數(shù)據(jù)結構對數(shù)據(jù)擬合，更具靈活性；(2)無需定義回歸參數(shù)，可自動選擇；(3)相對于廣義線性模型(GLM)等其他模型，GAM模擬精度更高；(4)更容易發(fā)現(xiàn)數(shù)據(jù)中的線性關系，并可以處理一個響應變量與多個解釋變量之間的非線性關系。假設函數(shù)是相加的，且函數(shù)的組成成分是光滑函數(shù)(smooth functions)，其基本模型如下[24-25]:

g(EY)=a0+f1(X1)+…+fn(Xn)+ε。

(1)

式中：g(…)是連接函數(shù)；Y是收集地表甲蟲種群數(shù)量，EY是Y的期望值；a0是截矩；在fn(Xn)項中，X1,…,Xn分別為各環(huán)境因子，fn(…)是解釋變量Xn的單變量函數(shù)；ε為隨機變量。通常使用光滑樣條函數(shù)(smoothing spline function)擬合，因子的重要性可以通過平滑函數(shù)的自由度來簡單判斷：自由度高表明因子重要性高。本研究分析在R軟件的mgcv包支持下和Arcgis10.5中完成。

1.4 數(shù)據(jù)處理

1.4.1 地形因子 利用ArcGIS數(shù)字化1∶10 000地形圖，獲得研究區(qū)DEM數(shù)據(jù)(分辨率為5 m)，分析提取研究區(qū)海拔、坡度和坡向，用于地表甲蟲物種空間分布的相關性分析。其中坡向和坡度采用3次有限差分方法[26]計算得到，坡向的組別以國際劃分法確定。根據(jù)國家水土保持技術通則將坡度劃分為6個等級[27]?？紤]到研究區(qū)相對高差不大，因此按照50 m間隔將流域劃分為5個高程區(qū)。根據(jù)樣地地形特征，將海拔、坡度、坡向作為模型的預測變量。

1.4.2 土壤水分數(shù)據(jù) 由于土壤含水量受降雨事件影響較大，本研究選擇5～7 d連續(xù)無降雨時段取樣，以每個樣點中央陷阱器為中心，在直徑2 m的范圍內用土鉆取得0～30 cm土壤，帶回實驗室采用烘干法測得。

2 結果與分析

2.1 黃土丘陵區(qū)地表甲蟲種群數(shù)量及空間分布

本次采樣共捕獲甲蟲標本共5 347個，分屬于11科60余種，地表甲蟲種群數(shù)量見表2。由表2可知，步甲科(Carabidae)昆蟲在高建堡流域內屬于廣布種和優(yōu)勢種，共收集到3 844個樣本，占地表甲蟲總數(shù)的71.90%，因此其適合作為分析地表甲蟲與環(huán)境因子關系的研究對象；其次為金龜科(Scarabaeidae)和虎甲科(Cicindelidae)，所占比例分別為13.62%和8.64%。

由表2還可知，在高建堡小流域的6種典型生境中，灌草混交林地、喬灌混交林地和生態(tài)薪炭林地是地表甲蟲生存、棲息的較佳生境，地表甲蟲個體數(shù)量分布較多；其他生境地表甲蟲種群數(shù)量排序依次為水平農(nóng)田地>生態(tài)經(jīng)濟林地>天然封育草地。由此可知,不同生境地表甲蟲種群分布存在很大的差異性。步甲科(Carabidae)昆蟲在生態(tài)薪炭林地分布最多，其次為喬灌混交林地，灌草混交林地、生態(tài)經(jīng)濟林地和水平農(nóng)田地中的數(shù)量相近，天然封育草地最少；金龜科(Scarabaeidae)甲蟲主要分布于灌草混交林地和喬灌混交林地；虎甲科(Cicindelidae)昆蟲在喬灌混交林地、天然封育草地和灌草混交林地中分布較多；其他科類地表甲蟲個體數(shù)量較少，在6種生境中的分布無明顯規(guī)律。

表2 黃土丘陵區(qū)寧夏高建堡流域地表甲蟲種群數(shù)量Table 2 Surface beetles quantity in the Loess Hilly Region in Ningxia

2.2 黃土丘陵區(qū)環(huán)境因子的t檢驗

環(huán)境因子能否作為解釋變量，首先需要對環(huán)境因子進行t檢驗，確定它們之間是否存在共線性相關。經(jīng)t檢驗后，如果兩個因子之間的Pearson相關系數(shù)R<0.5，說明兩個因子間不存在共線性相關，因此這兩個因子均可作為解釋變量。黃土丘陵區(qū)寧夏高建堡流域環(huán)境因子的t檢驗結果如表3所示。由表3可知，本研究所選環(huán)境因子間的Pearson相關系數(shù)R均小于0.5，說明所選環(huán)境因子均可以作為解釋變量，可全部納入模型。

表3 黃土丘陵區(qū)寧夏高建堡流域環(huán)境因子t檢驗結果Table 3 t test results of environmental factors in the Loess Hilly Region in Ningxia

2.3 GAM計算結果

通過對比不同解釋變量組成函數(shù)的GAM結果，選取模型方差解釋量最大、廣義交叉驗證值最小、F檢驗概率值(P)精度最高的模型作為最優(yōu)模型。一般情況下，當F檢驗概率值P<0.05，表明解釋變量對響應變量影響顯著；如果調整擬合系數(shù)(R2)大于0.5，說明模型具有良好的穩(wěn)定性，易于解釋響應變量與解釋變量之間關系。由表4可知，本研究所涉及環(huán)境因子均對步甲科甲蟲分布有顯著影響。

選取黃土丘陵區(qū)對地表甲蟲生境影響最直接的地形因子(海拔、坡度、坡向)和土壤水分因子作為解釋變量，利用GAM分析和獲取其相關性，結果顯示，地表步甲科昆蟲種群的個體數(shù)量與海拔、坡度、坡向、土壤水分均顯著相關(表4)。

表4 基于不同環(huán)境因子的黃土丘陵區(qū)寧夏高建堡小流域地表步甲科甲蟲生境的GAM運行結果Table 4 GAM results for surface beetle habitat in the Gaojianbao watershed of Ningxia in the Loess Hilly Region based on different environmental factors

黃土丘陵區(qū)寧夏高建堡流域步甲科甲蟲種群個體數(shù)量與環(huán)境因子的關系如圖1所示。由圖1-A可知，步甲科昆蟲種群的個體數(shù)量與海拔總體呈非線性正相關(F=9.449，P<0.001)。在海拔分別為1 650～1 680，1 700～1 720，1 750～1 780，1 820～1 850 m區(qū)域，步甲科甲蟲群落個體數(shù)量與海拔之間主要呈正相關；在海拔1 680～1 700和1 780～1 820 m區(qū)域，二者之間主要呈負相關。

由圖1-B可知，步甲科昆蟲種群的個體數(shù)量與坡度總體呈非線性負相關(F=3.144，P<0.001)。在坡度為5°～8°，12°～15°和19°～23°區(qū)域內，二者之間主要呈正相關；在坡度為8°～12°、15°～19°和23°～26°區(qū)域，二者之間主要呈負相關。

由圖1-C可知，步甲科昆蟲種群的個體數(shù)量與坡向總體呈現(xiàn)“w”形(F=2.459，P<0.001)。在坡向0°～50°和155°～225°呈負相關；在坡向51°～154°和226°～360°呈正相關。

由圖1-D可知，步甲科昆蟲種群的個體數(shù)量與土壤水分總體呈非線性正相關(F=31.62，P<0.001)。土壤含水量在10%～25%區(qū)域，二者之間呈負相關；在土壤含水量為25%～45%區(qū)域，二者之間呈正相關。

縱坐標為各自變量(環(huán)境因子)非線性解釋函數(shù)中的殘差值。曲線的峰值表示自變量與因變量間的相關性，陰影代表樣本數(shù)量The ordinate is the residual value in the nonlinear interpretation function of the respective variable (environmental factor).The peaks and troughs of the curve represent the positive and negative correlations between independent and dependent variables,respectively,and the shade represents the number of samples圖1 黃土丘陵區(qū)寧夏高建堡流域步甲科甲蟲種群個體數(shù)量與環(huán)境因子的關系Fig.1 Nonlinear relationship between number of carabid beetles and environmental factors in the Loess Hilly Region in Ningxia

3 討 論

3.1 不同生境地表甲蟲個體數(shù)量差異

棲息地的環(huán)境因子對地表甲蟲的生存、繁衍、個體組成和分布規(guī)律有很大的影響，鑒于地表甲蟲物種對環(huán)境因子有趨向性，故其數(shù)量能夠作為生物學指標在環(huán)境監(jiān)測中應用[28]。本研究結果顯示，高建堡小流域地表甲蟲群落的總體特點為：類群分布不均勻，多樣性低，優(yōu)勢度高，且不同生境差異顯著，其中步甲科甲蟲的種類和數(shù)量在流域內占據(jù)絕對優(yōu)勢，屬于廣布種和優(yōu)勢種；其次為金龜科，這與前人在華北地區(qū)[29-30]、西北荒漠地區(qū)[31]和南方地區(qū)[32-33]的研究結果不一致，可能是不同區(qū)域地理位置和地表甲蟲棲息環(huán)境存在差異所致。

3.2 環(huán)境因子對地表甲蟲種群的影響

本研究所選取的坡度、坡向、海拔和土壤水分均對步甲科甲蟲數(shù)量有顯著影響，其中地表步甲科昆蟲種群數(shù)量與海拔、土壤水分總體呈現(xiàn)非線性正相關，說明土壤水分對步甲科甲蟲空間分布造成的影響最大，直接影響其生境。已有研究表明，不同生境內土壤含水量與甲蟲群落特征值具有較高的相關性，其是影響地表甲蟲物種豐富度的關鍵因子[31,34-35]，這與本研究結論一致。也有學者對不同研究區(qū)土壤含水量與地表甲蟲的關系進行研究，結果卻不盡相同[30,36-37]，這說明昆蟲物種多樣性隨環(huán)境的變化在不同地區(qū)、不同尺度、不同生態(tài)系統(tǒng)類型等的結果存在差異, 其具體機制有待進一步研究。

本研究中，步甲科甲蟲的個體數(shù)量與坡度總體呈現(xiàn)非線性負相關，且相關性較弱。當坡度大于25°時地表甲蟲與坡度之間基本保持穩(wěn)定負相關，說明在坡度小于25°區(qū)域適宜地表甲蟲生存和繁衍。本研究中由于不同生境冠層蓋度和草本蓋度不同，直接改變地表甲蟲棲息環(huán)境的光熱濕等條件，且地表步甲科昆蟲多在陰濕地表生存，因此其種群的個體數(shù)量隨坡向變化總體呈現(xiàn)“w”形。地表步甲科昆蟲種群的個體數(shù)量與海拔總體呈現(xiàn)非線性正相關，但海拔是決定地表甲蟲棲息環(huán)境的直接因素，它影響區(qū)域內土壤含水量、溫度等因子，進而影響地表甲蟲空間分布和個體數(shù)量特征。

4 結 論

在研究區(qū)6種典型生境中，灌草混交林地是地表甲蟲生存、棲息的最佳生境，步甲科甲蟲在流域內屬于廣布種和優(yōu)勢種。研究所選取的4個環(huán)境因子對地表步甲科昆蟲種群數(shù)量都具有顯著的相關性，且不同環(huán)境因子對地表甲蟲數(shù)量和空間分布影響存在明顯差異。其中，土壤水分對地表甲蟲空間分布影響最大，但海拔是決定地表甲蟲棲息環(huán)境的直接因素。