趙閃閃, 褚一凡, 姜小玉, 于飛, 靳同霞, 馬劍敏
(河南師范大學生命科學學院,河南新鄉453007)
河南新鄉黃河濕地鳥類國家級自然保護區是典型的河流灘涂濕地,位于我國候鳥遷徙中線的必經之路,鳥類多樣性較高。據最新記載,該地共有鳥類153種,占河南省鳥類總數的42.7%,其中留鳥24.8%,夏候鳥26.1%,冬候鳥28.1%,旅鳥20.9%(內部資料)。然而目前,保護區生境破碎化程度嚴重(徐文茜等,2016),水域面積減少,旅游開發等人為活動日益頻繁,鳥類棲息地受到了威脅。Tilman等(2014)研究表明,在所有影響物種多樣性的因素中,棲息地喪失與破碎化是造成多樣性下降的最主要因素。為了維持保護區鳥類多樣性、了解影響鳥類棲息繁殖的主要因素、進行科學有效的管理,開展鳥類與棲息地之間的關系研究尤為必要。
蒼鷺Ardeacinerea隸屬于鵜形目Pelecaniformes鷺科Ardeidae,不僅在食物鏈中處于較高的營養級,而且是濕地生態系統中重要的指示物種。濕地環境變化勢必對蒼鷺的種群和繁殖產生一定的影響(Sujaketal. ,2016)。根據相關記載(牛紅星等,2001;牛俊英等,2009),蒼鷺常年在新鄉黃河濕地棲息,被視為常見種,是保護區涉禽的代表。為了更好地養育后代,大多數鳥類會在領域內尋找適宜地點,建造不同類型的鳥巢(鄭光美,2012)。巢是鳥類繁育后代的特殊場所,在鳥類繁殖期間具有容納卵或雛鳥、保溫、保護、促進繁殖等功能。巢特征在某種程度上也能反映鳥類的繁殖投入(窩卵數)和雙親的供應能力(Neve & Soler,2002)。另外,巢址選擇是鳥類繁殖的基礎,是鳥類順利完成繁殖的先決條件(劉利等,2015)。適宜的巢址能有效保護卵和雛鳥,降低干擾以及預防天敵捕食,并能從中得到繁殖所需的食物和隱蔽場所,提高鳥類繁殖成功率(Jakubas,2005;陳偉等,2006;龍帥等,2007;王寧等,2007)。巢址選擇取決于當地的環境條件,由雌雄鳥共同選擇或者由雌鳥或雄鳥單獨選擇,不僅要把天敵的捕食壓力降到最低,而且要具備利于胚胎和雛鳥生長發育的熱力學條件(丁長青,鄭光美,1997)。吳慶明等(2008)對扎龍濕地蒼鷺巢址的選擇研究表明,巢距人為活動不頻繁地距離和植被高度是影響蒼鷺巢址選擇的重要因子;劉利等(2015)對內蒙古庫布齊沙漠蒼鷺巢址選擇的研究表明,蒼鷺傾向于距池塘較近、道路較遠、植株密度適中、蓋度較高的地方營巢,食物和人類的干擾條件是巢區選擇的關鍵因素。Lusk等(2017)對百靈雀Chondestegrammacus巢址選擇的研究發現,植被特征是影響其巢址選擇的生態因子。迄今為止,尚未有關于新鄉黃河濕地蒼鷺巢特征和巢址選擇的研究報道。
蒼鷺的繁殖期為每年3—8月,其中,占巢期為3—4月(高欣,劉明玉,2002),巢域固定期為4月末,8月中旬開始集群活動(李曉民,陶宇,1990)。2017年5—8月,采用野外觀察、GPS定位、樣方調查等方法對保護區蒼鷺巢區進行調查研究,以了解蒼鷺巢特征、巢周環境特征以及主要影響因素,為恢復構建濕地鳥類群落、有效管理保護區提供科學依據;為保護區生物多樣性保護提供基礎資料。
新鄉黃河濕地于1988年經河南省人民政府批準建立省級自然保護區,1996年晉升為國家級鳥類自然保護區。保護區(114°13′53″~114°52′30″E,34°53′13″~35°06′21″N)地處河南省新鄉市東部封丘縣和長垣縣境內,北岸是長垣縣魏莊鎮至封丘縣荊隆宮鄉與原陽縣交界處,東岸是河南與山東省界,南岸是開封市蘭考縣至開封市水稻鄉;總面積22 780 hm2,海拔73~80 m。屬大陸性季風氣候,受冷暖空氣交替影響,四季分明,年平均氣溫14 ℃,年平均降水量466.9 mm。據最新的科學考察(內部資料),區內有維管束植物102科392屬903種、兩棲類14種、昆蟲319種、獸類14種、底棲類41種、魚類23種、鳥類153種。水生植物主要有蘆葦Phragmitesaustralis、東方香蒲Typhaorientalis等,陸生植物的優勢科有禾本科Poaceae、菊科Asteraceae、豆科Fabaceae、莎草科Cyperaceae。大天鵝Cygnuscygnus、大鴇Otistarda、灰鶴Grusgrus等珍稀物種會在此越冬。
蒼鷺的繁殖棲息地主要位于陳橋東湖濕地,該地區位于黃河大堤北側,常年有黃河水的側滲和補給。濕地中央的植被有蘆葦、東方香蒲等。在蒼鷺繁殖期,還有黃斑葦Ixobrychussinensis、夜鷺Nycticoraxnycticorax、白琵鷺Platalealeucorodia等在此筑巢繁殖,但數量較少。
2017年5—8月,在保護區的2個固定點用雙筒望遠鏡(10×42)和單筒望遠鏡(Nikon 60ED-A)進行觀察,確定蒼鷺的棲息范圍,鎖定蒼鷺巢穴的大致方位,并租船接近其巢穴,用GPS定位每個巢穴的準確位置,同時在最短的時間內用卷尺對巢的形態特征進行測量,設置樣方對巢周的環境特征進行測量。巢形態特征的測量指標有:巢外徑、巢內徑、巢高、巢深、窩卵數;巢址區環境特征的測量指標有:巢址水深、巢周植被種類、植被高度、植被密度、植被蓋度、巢址距道路距離等。同時,在巢區周圍與巢域大環境一致且沒有蒼鷺巢穴的生境中,選擇可能影響蒼鷺巢址選擇的生態因子進行觀測記錄,主要包括植物種類、植被高度、植被蓋度、植被密度、水深、距人為活動頻繁地的距離、距人為活動不頻繁地的距離、水透明度等,將其記為對照區生態因子。
蒼鷺巢址區與對照區的生態因子數據用Kolmogorov-SmirnovZ檢驗是否符合正態分布,若符合正態分布,使用獨立樣本t檢驗;反之,則使用Mann-WhitneyU檢驗。各變量在描述時采用平均值±標準差表示(王永強,2007)。
為了避免多重共線性的影響,需要對差異有統計學意義的變量進行Spearman相關檢驗(雙尾),若變量間的相關性較高,需要先將生態因子納入主成分分析,然后將保留的主成分分析的因子系數再納入二元Logistic回歸中分析影響蒼鷺巢址選擇的關鍵因子(Dietmaretal. ,2004)。最后選取具有最小Akaike信息標準(Akaike information criteria,AIC)值或AICc值的回歸等式為蒼鷺巢址選擇的最佳模型(Strawdermanetal. ,1998;Boyceetal. ,2002;賈非等,2005;Pan,2015),同時計算最優模型與次優模型的ΔAIC值或ΔAICc值,如果其值<2,則需要對模型進行平均(Symonds & Moussalli,2011;李東來等,2015;吳奕如等,2016)。AIC和AICc的計算公式為:
AIC=-2ln(L)+2K,其中,K為回歸方程中自變量的個數+2,L為似然函數。
AICc=AIC+2K(K+1)/(n-K-1),其中,n為樣本數,當n/K≥40時,采用AIC值作為模型預測準確性的判定標準;當n/K<40時,使用AICc值。
以上所有統計分析在SPSS 17.0和R 3.4.2中完成。
共記錄了蒼鷺20個巢的特征,巢址區域為114°30′28″~114°30′80″E,34°56′30″~34°56′80″N。蒼鷺巢均建在蘆葦叢之上,屬于地棲型巢。近距離觀察發現,筑巢的巢材為蘆葦的莖和葉,長短不等,莖的直徑為0.5~1.0 cm。對巢形態特征各指標排序后,記錄四分位數(Dariusz & Brygida,2011)(表1),巢內徑為(39.00±3.28) cm,巢外徑為(79.10±15.39) cm,巢高為(31.00±19.80) cm,巢深為(5.00±1.43) cm。窩卵數3~5枚,以4枚居多。

表1 河南新鄉黃河濕地鳥類國家級自然保護區蒼鷺巢特征Table 1 The nests characteristics of Ardea cinerea inXinxiang Yellow River Wetland Birds National Nature Reserve, Henan
注: Q為巢特征指標的四分位點, Q1為1/4處的蒼鷺巢特征指標, Q3為3/4處的蒼鷺巢特征指標
Notes: Q is the quartile, Q1 is the first quartile ofA.cinereacharacteristic index, Q3 is the third quartile
共設置和調查了20個巢址區樣方和20個對照區樣方。經Kolmogorov-SmirnovZ檢驗,所有變量均符合正態分布(P>0.05)。獨立樣本t檢驗結果顯示(表2),二者植被密度間的差異有統計學意義(P<0.05);植被高度、植被蓋度、水深、距人為活動不頻繁地的距離、距人為活動頻繁地的距離之間的差異有高度統計學意義(P<0.01);而水透明度之間的差異無統計學意義(P>0.05)。說明蒼鷺喜歡在植被高度較低、植被蓋度較低、植被密度適中、水深相對較深、距人為活動干擾地距離較遠的生境中筑巢。

表2 河南新鄉黃河濕地鳥類國家級自然保護區蒼鷺巢址區與對照區生態因子的對比分析Table 2 Comparative analysis of ecological factors in nest sites and control sites of Ardea cinereain Xinxiang Yellow River Wetland Birds National Nature Reserve, Henan
對巢址區和對照區之間差異具有統計學意義的變量給予保留,并對保留變量進行Spearman相關性檢驗(雙尾),結果發現各因子之間的相關系數較高(表3),因此,將巢址區和對照區生態因子納入主成分分析。

表3 保留變量的Spearman相關分析結果(雙尾)Table 3 The remained variables by Spearman correlation analysis (two-tailed)
注Note:*P<0.05,**P<0.01
主成分分析獲得3個主成分(表4),根據各主成分的特征向量分別命名為植被狀態、干擾因子和水深,累計貢獻率達到81.49%,說明其基本包含所有生態因子的大部分信息。將3個主成分的因子系數納入二元Logistic回歸,共構建8個模型。計算AICc值結果(表5)發現,最優模型和次優模型的ΔAICc值<2,經模型平均獲得三者的相對重要值分別為1.00、0.99和0.57(表6),系數估計分別為-0.381±0.047、0.168±0.048和0.045±0.046。由此可見,影響蒼鷺巢址選擇的生態因子的重要性順序依次為:植被狀態、干擾因子和水深。

表4 蒼鷺巢址生態因子的主成分分析Table 4 Principal component analysis of Ardea cinerea’snest-site ecological factors
注: F1. 植被狀態, F2. 干擾因子, F3. 水深; 下同
Notes: F1. vegetation status, F2. interference factor, F3. water depth; the same below

表5 巢址選擇的模型選擇Table 5 Models for nest-site selection
注: —. 沒有變量入選, 只有常數項
Notes: —. No varible was selected, only a constant term

表6 蒼鷺巢址選擇生態因子的重要值和模型平均后的系數估計Table 6 Relative importance of Ardea cinerea nest-siteselection’s factors and weighted average parameter estimates
通過查閱國內6個地區蒼鷺巢特征、窩卵數、緯度等信息,與本調查數據對比發現,7個地區蒼鷺巢特征指標數值不一樣(表7)。整體來看,新鄉黃河濕地蒼鷺的巢特征與陜西省勉縣銅錢壩鄉(史東仇,路寶忠,1991)的數值相近。雖然后者的蒼鷺巢型為樹棲型,但二者的巢外徑、巢內徑、巢高、巢深、窩卵數均相近,可能是兩地區的緯度相近,在蒼鷺繁殖期間,兩地的氣候條件相近,從而形成相似的巢特征。而同為樹棲型、緯度更高的遼寧省朝陽縣蘇家營子(袁宏宇等,2003)的蒼鷺巢,則具有巢較小、較深、窩卵數較多的特點。
蒼鷺巢深在7個地區間不同,新鄉黃河濕地的數值與陜西省勉縣銅錢壩鄉的相近,但明顯小于另外5個地區。對巢深和緯度進行相關分析,發現二者呈顯著正相關關系(圖1)。可能與不同緯度所導致的氣候差異有關。在高緯度地區,氣溫相對較低,大巢深有利于保溫,從而使孵卵親鳥減少對卵加熱的能量支出。

表7 新鄉黃河濕地、扎龍濕地等7個地區蒼鷺巢特征的比較Table 7 Comparison of Ardea cinerea nest characteristics in Xinxiang Yellow River Wetland, Zhalong Wetland and other 5 regions

圖1 7個蒼鷺筑巢地區的緯度與巢深的關系
Fig. 1 Relationship between the latitude and nest depth ofArdeacinereanests in 7 regions
新鄉黃河濕地蒼鷺窩卵數以4枚居多,與陜西省勉縣銅錢壩鄉、寧夏永寧王太村的數據相近,但小于七星河濕地和遼寧省朝陽縣蘇家營子。窩卵數整體上隨著緯度的增加呈增加趨勢,此現象符合鄭光美等(2012)得出的繁殖地從熱帶向南、北兩極移動時,同種鳥類的窩卵數逐漸增多的結論。Ricklef(1990)認為,自然選擇使生活在高緯度地區的鳥類具有較多的窩卵數。這屬于動物的繁殖策略。
對比7個地區蒼鷺的巢材發現:新鄉黃河濕地、七星河濕地、扎龍濕地、黑龍江牡丹江地區(地棲型)4個濕地均為蘆葦;遼寧省朝陽縣蘇家營子為小楊樹Populussimonii、旱柳Salixmatsudana、刺槐Robiniapseudoacacia等;陜西省勉縣銅錢壩鄉為樹枝、藤本等;黑龍江牡丹江地區(樹棲型)為榆樹枝、蒿莖、桑科莖葉等;寧夏永寧王太村地區為樹枝、雜草、樹葉等。綜合分析,蒼鷺在選擇巢材時偏向于就地取材,與周圍的植物種類有很大關系。
巢址選擇是鳥類繁殖的重要環節,涉及食物的可獲得性、天敵的捕食、人類活動的干擾和遮蔽物的狀態等重要因素(丁未等,2011)。與對照區相比,蒼鷺偏好在高度適中(1.51 cm±0.36 cm)、蓋度較低(32.00%±26.43%)、密度適中(111.00株·m-2±77.04株·m-2)的蘆葦叢中筑巢(相較于當地許多連片生長、密度可達300株·m-2、蓋度幾乎100%的茂密蘆葦蕩,其巢區的蘆葦蓋度較低,密度和高度適中)。蒼鷺選擇此環境筑巢的原因是:(1)巢區蘆葦常年生長,無收割現象,硬度足以支撐巢重,蘆葦的高度、密度可以起到遮蔽作用;(2)與蒼鷺食物的可獲得性有關。蘆葦等水生植物密度、蓋度過高不易獲取水中的食物。另據觀察,在巢區正東50 m處,長有4棵高5~6 m的檉柳Tamarixchinehsis,但其上以及周圍其他樹上未見蒼鷺巢穴。在新鄉黃河濕地,蒼鷺選擇蘆葦叢筑巢,而不選擇樹棲的原因需進一步研究。
本研究發現,巢下水深是另一個影響蒼鷺巢址選擇的重要因素。巢下水深平均為89.45 cm,與對照區差異顯著。走訪村民發現,蒼鷺所在的巢區是濕地中心地勢最低的一個水塘,不管濕地整體水位上升或者下降,該地區水深受影響較小。說明蒼鷺喜歡選擇有一定水深的地方筑巢,此現象與吳慶明等(2008)的報道一致。
調查還發現,距人為活動頻繁地的距離和距人為活動不頻繁地的距離對蒼鷺巢址的選擇也有重要影響。實地走訪發現,在蒼鷺筑巢所在區域內有小的土路,當地村民有時會駕駛機動車輛進出,在調查期間定義為人為活動不頻繁地。在蒼鷺巢區的外部邊緣修有水泥路,車輛出入較頻繁,且在巢區的正北方向有一度假村,人員活動頻繁,視其為人為活動頻繁地。這些人為活動對蒼鷺的繁殖和棲息產生了較大的干擾。
模型平均的結果顯示,巢周植被狀態和干擾因子是影響蒼鷺巢址選擇的關鍵因素。總的來說,蒼鷺巢址選擇是對巢周植被狀態、干擾因子、水深等權衡的結果。
鳥類在絕大多數濕地生態系統中處于頂級群落,其群落研究對物種多樣性保護具有重要的意義。蒼鷺為大型涉禽,在保護區鳥類群落中具有獨特地位,其繁殖成效和種群數量對保護區的鳥類棲息繁殖環境具有指示作用。因此,本研究對于評估保護區生態環境狀況、有效管理保護區具有重要價值。近幾年,在鳥類棲息地的周邊建有曹崗萬畝油菜花基地和宋源湖風景區等旅游開發項目,對鳥類棲息環境和濕地生態系統造成了一定影響。為了更好地保護濕地生態系統,維護鳥類的多樣性,建議:(1)適當擴大水域面積,增加灘地,擴大鳥類的可棲息空間和取食區域;(2)適當封鎖濕地中的道路,將人為活動的干擾降到最低;(3)加強濕地保護區的有效管理,定期蓄水,做好防火等工作;(4)適當擴大蘆葦、香蒲等水生植物分布面積,為鳥類棲息繁殖和取食提供遮蔽物和食材;(5)處理好濕地管理與旅游開發之間的關系,優化旅游規劃項目,在確保不會對生態環境造成明顯損害的前提下,適當發展生態旅游。
致謝:本調查的野外工作得到河南新鄉黃河濕地鳥類國家級自然保護區管理處班永田站長和陳橋村民劉立新、劉東洋的大力支持與幫助,誠表謝意。
:
陳偉, 吳映梅, 胡錦矗, 等. 2006. 四川南充市城市白頭鵯的巢址選擇和巢材研究[J]. 四川動物, 25(3): 590-593.
丁未, 劉迺發, 王亮, 等. 2011. 黑頂麻雀的巢址選擇[J]. 四川動物, 30(6): 928-931.
丁長青, 鄭光美. 1997. 黃腹角雉的巢址選擇[J]. 動物學報, (1): 27-33.
費殿金, 車仁富, 楊正明. 1985. 扎龍保護區草鷺、蒼鷺繁殖習性的初步研究[J]. 動物學雜志, (2): 12-16.
高欣, 劉明玉. 2002. 蒼鷺營巢習性的初步觀察[J]. 四川動物, 21(2): 96-97.
高中信, 賈競波, 閆文, 等. 1991. 蒼鷺繁殖生態研究[J]. 東北林業大學學報, (3): 35-40.
賈非, 王楠, 鄭光美. 2005. 冬季白馬雞群體夜棲地特征分析[J]. 生態學雜志, 24(2): 153-158.
李東來, 魏宏偉, 孫興海, 等. 2015. 震旦鴉雀在鑲嵌型蘆葦收割生境中的巢址選擇[J]. 生態學報, 35(15): 5009-5017.
李曉民, 陶宇. 1990. 蒼鷺年周期的研究[J]. 國土與自然資源研究, (4): 68-69.
劉利, 劉曉光, 張樂, 等. 2015. 內蒙古庫布齊沙漠蒼鷺巢址選擇研究[J]. 干旱區資源與環境, 29(5): 144-148.
龍帥, 周材權, 王維奎, 等. 2007. 南充雉雞的巢址選擇和春夏季棲息地選擇[J]. 動物學研究, 28(3): 249-254.
牛紅星, 卜艷珍, 路紀琪, 等. 2001. 豫北黃河故道濕地鳥類自然保護區鳥類區系調查[J]. 生物多樣性, 9(3): 260-264.
牛俊英, 馬朝紅, 馬書釗, 等. 2009. 河南黃河濕地國家級自然保護區鳥類資源調查[J]. 四川動物, 28(3): 462-467.
任青峰. 1994. 寧夏蒼鷺和夜鷺的數量及生態研究[J]. 寧夏大學學報(自然科學版), 15(2): 71-75.
史東仇, 路寶忠. 1991. 蒼鷺的繁殖生物學觀察研究[J]. 陜西師范大學學報(自然科學版), 19(3): 52-55.
王寧, 張雁云, 鄭光美. 2007. 北京小龍門林場黃眉姬鹟的巢與巢址特征[J]. 動物學研究, 28(4): 337-343.
王永強. 2007. 內蒙古東部草原狼(Canislupus)越冬期生境選擇[D]. 曲阜: 曲阜師范大學.
吳慶明, 繆克傳, 鄒紅菲, 等. 2008. 扎龍濕地蒼鷺巢址選擇[J]. 野生動物學報, 29(1): 26-28.
吳奕如, 斯幸峰, 陳傳武, 等. 2016. 千島湖陸橋島嶼繁殖鳥類的擴散能力差異對群落動態的影響[J]. 生物多樣性, 24(10): 1135-1145.
徐文茜, 湯茜, 丁圣彥. 2016. 河南新鄉黃河濕地鳥類國家級自然保護區景觀格局動態分析[J]. 濕地科學, 14(2): 235-241.
袁宏宇, 周正, 任秀奇, 等. 2003. 蒼鷺繁殖生態習性研究[J]. 遼寧林業科技, (6): 10-12.
鄭光美. 2012. 鳥類學[M]. 北京: 北京師范大學出版社.
朱井麗, 吳慶明. 2012. 七星河保護區蒼鷺巢特征及巢周環境[J]. 野生動物學報, 33(3): 143-145.
Boyce MS, Vernier PR, Nielsen SE,etal. 2002. Evaluating resource selection functions[J]. Ecological Modelling, 157(2): 281-300.
Dariusz J, Brygida M. 2011. The response of grey heronsArdeacinereato changes in prey abundance[J]. Bird Study, 58(4): 487-494.
Dietmar Z, Frieder K, Zellner GE. 2004. Variable selection in logistic regression models[J]. Communications in Statistics-Simulation and Computation, 33(3): 787-805.
Jakubas D. 2005. Factors affecting the breeding success of the grey heron (Ardeacinerea) in northern Poland[J]. Journal of Ornithology, 146(1): 27-33.
Lusk JJ, Wells KS,Guthery FS,etal. 2017. Lark sparrow (Chondestesgrammacus) nest-site selection and success in a mixed-grass prairie[J]. Auk, 120(Jan 2003): 120-129.
Neve LD, Soler JJ. 2002. Nest-building activity and laying date influence female reproductive investment in magpies: an experimental study[J]. Animal Behaviour, 63(5): 975-980.
Pan W. 2015. Akaike’s information criterion in generalized estimating equations[J]. Biometrics, 57(1): 120-125.
Ricklefs RE. 1980. Geographical variation in clutch size among passerine birds: ashmole’s hypothesis[J]. Auk, 97(1): 38-49.
Strawderman RL, Burnham KP, Anderson DR. 1998. Model selection and inference: a practical information-theoretic approach[M]. New York: Springer: 76-81.
Sujak A, Kusz A, Rymarz M,etal. 2016. Environmental bioindication studies by Bayesian network with use of grey heron as model species[J]. Environmental Modeling & Assessment, 22(2): 103-113.
Symonds MRE, Moussalli A. 2011. A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using akaike’s information criterion[J]. Behavioral Ecology and Sociobiology, 65(1): 13-21.
Tilman D, Fargione J, Wolff B,etal. 2001. Forecasting agriculturally driven global environmental change[J]. Science, 292: 281-284.