孫雪瑩,李祎斌,吳慶明, *,徐 卓,2,高曉冬,李玉春,柳 旭,鄒紅菲,李曉民

1 東北林業(yè)大學(xué)野生動物與自然保護地學(xué)院,哈爾濱 150040

2 黑龍江扎龍國家級自然保護區(qū)管理局,齊齊哈爾 161000

3 曲阜師范大學(xué)生命科學(xué)學(xué)院,曲阜 273165

4 煙臺市森林資源監(jiān)測保護服務(wù)中心,煙臺 264001

5 黑龍江省寶清東升自然保護區(qū)管理局,寶清 155600

棲息地是指野生動物的生存環(huán)境空間,這一空間可為野生動物提供完成生活史過程所需要的所有生物與非生物需求,所有這些需求都關(guān)系著野生動物的生存繁衍與種群延續(xù),不僅有宏生境層面數(shù)量的配置供給,更有微生境層面質(zhì)量的精細需要[1-4]。這種情況下,適宜野生動物生存的棲息地則顯得尤為重要,這是野生動物生存的必要條件,對野生動物的生存起決定性作用,鳥類也不例外,對候鳥最為關(guān)鍵[4-6]。因此,基于野生動物現(xiàn)存的棲息地環(huán)境特征,分析其潛在的棲息地并對適宜性分布進行分析,有助于對特定區(qū)域野生動物及其棲息地資源開展針對性的科學(xué)保護[7-8]。

松嫩平原是東北平原的重要組成部分,不僅是我國及世界重要的黑土帶分布區(qū),更是世界九條鳥類遷徙路線中最為重要的東亞-澳大利西亞鳥類遷徙路線的重要組成[9-10]。蒼鷺(Ardeacinerea)是該區(qū)域最為常見的大型濕地鳥類,但繁殖種群數(shù)量遠低于歷史規(guī)模,歷史上的大規(guī)模繁殖群現(xiàn)在已很難再現(xiàn)[11-12]。為了避免蒼鷺這一常見鳥種出現(xiàn)青頭潛鴨(Aythyabaeri)種群驟減且科學(xué)研究資料缺乏的現(xiàn)象,非常有必要對蒼鷺的保護生物學(xué)進行密切關(guān)注和深入思考[13-14]。

而關(guān)于野生動物潛在棲息地分析及棲息地適宜性評價方面的研究,已被眾多學(xué)者認(rèn)為是開展瀕危物種保護工作極為重要的生物學(xué)信息,不僅能指導(dǎo)開展具體的棲息地資源管理與保護,而且能有針對性地挖掘潛在棲息地并進行棲息地資源質(zhì)量的管理預(yù)警[15-17]。目前,已有眾多學(xué)者對丹頂鶴(Grusjapoenesis)、白頭鶴(Grusmonacha)、白鶴(Grusleucogeranus)、白枕鶴(Grusvipio)、灰鶴(Grusgrus)、黑頸鶴(Grusnigricollis)、朱鹮(Nipponianippon)、小鴇(Tetraxtetrax)、鷸類(Scolopacidae)、啄木鳥(Picidae)等鳥類棲息地適宜性開展了一系列的評價研究[1-2,18-27]。而關(guān)于蒼鷺棲息地適宜性評價方面,關(guān)注較少,僅有學(xué)者通過經(jīng)典的生境適合度指數(shù)模型(Habitat suitability index, HSI)對扎龍保護區(qū)蒼鷺營巢最適生境進行了分析[12]。為了了解松嫩平原蒼鷺秋遷期停歇地的適宜性分布信息,本文采用最新的GPS/GSM衛(wèi)星跟蹤技術(shù),并結(jié)合遙感影像和地理信息系統(tǒng)及Maxent模型,對該區(qū)域蒼鷺秋季潛在停歇地進行了分析,并對其適宜性進行了評價,以期為松嫩平原蒼鷺棲息地未來的管理和保護提供科學(xué)依據(jù)。

1 研究地概況

松嫩平原是中國東北三大平原之一,地理位置119°45′—129°36′E,42°50′—49°18′N,西連大興安嶺東麓,北部和東部與小興安嶺和長白山外緣山麓臺地相鄰,南以松遼分水嶺為界,跨黑龍江和吉林兩省,行政區(qū)域涵蓋哈爾濱、大慶、白城、齊齊哈爾、綏化、扶余、長春、四平9個地級行政區(qū)內(nèi)的55個縣級行政區(qū)[28-29]。松嫩平原地處溫帶半濕潤、半干旱的森林草甸與草甸草原地帶,大陸性季風(fēng)氣候明顯,四季分明,年均溫2—6℃,年均降水量可達400—600mm,降水多集中在6—9月[28-29]。松嫩平原源于松花江和嫩江兩大河流,其間分布的湖泡等濕地分布廣泛,這些成為丹頂鶴、白枕鶴、白頭鶴、白鶴、蒼鷺等許多大型水鳥的繁殖地和停歇地[9,11-13],也因此建立了吉林莫莫格國家級自然保護區(qū)(下文簡稱莫莫格保護區(qū))、黑龍江扎龍國家級自然保護區(qū)(下文簡稱扎龍保護區(qū))、吉林向海國家級自然保護區(qū)(下文簡稱向海保護區(qū))、吉林查干湖國家級自然保護區(qū)(下文簡稱查干湖保護區(qū))、吉林大布蘇國家級自然保護區(qū)(下文簡稱大布蘇保護區(qū))、黑龍江烏裕爾河國家級自然保護區(qū)(下文簡稱烏裕爾河保護區(qū))等重要保護地,并成為東亞-澳大利西亞鳥類遷徙路線的重要組成[13]。

2 研究方法

2.1 數(shù)據(jù)收集與篩選

2018年6月下旬,在黑龍江省齊齊哈市泰來縣泰湖國家濕地公園為五只蒼鷺雛鳥佩戴了衛(wèi)星跟蹤器(黑林護許準(zhǔn)[2018]32號)。衛(wèi)星跟蹤器生產(chǎn)于湖南環(huán)球信士科技有限公司,型號為鳥類背負(fù)式追蹤器HQBG3527,重量為40g,體積為65mm×35mm×50mm。跟蹤器佩戴個體的選擇原則上選取蒼鷺雛鳥各個器官基本發(fā)育完善的體型較大的個體,該階段雛鳥已能短距離飛行,能獨立覓食但捕食能力不強,每日早晚仍需親鳥喂食。通過網(wǎng)捕對捕捉到的符合體重要求的蒼鷺雛鳥個體佩戴衛(wèi)星跟蹤器。

關(guān)于本文分析所使用的數(shù)據(jù),選取的是蒼鷺在松嫩平原區(qū)域內(nèi)離開繁殖地后的活動位點數(shù)據(jù)。利用GPS/GSM衛(wèi)星跟蹤技術(shù)結(jié)合實地調(diào)查的方法獲取了2018年7月17日—10月24日間蒼鷺的空間位點數(shù)據(jù)。數(shù)據(jù)每小時采集一組,每組數(shù)據(jù)包括時間、溫度、定位精度(最高5m)、經(jīng)緯度、飛行速度、飛行高度等信息。為了真正反映蒼鷺對空間的棲息利用情況并提升數(shù)據(jù)精度,對采集的5976組數(shù)據(jù)進行了篩選,通過剔除重復(fù)位點、選取定位精度為5m且連續(xù)停留時間大于1h的棲息位點作為蒼鷺選擇的棲息位點,共獲得棲息位點231個(圖1)。這些棲息位點被認(rèn)為是被蒼鷺選擇有繼續(xù)利用概率的棲息點。其中,重復(fù)位點的考慮原則以距離超過10km為準(zhǔn)。

圖1 蒼鷺在松嫩平原的棲息位點

2.2 環(huán)境變量數(shù)據(jù)

棲息地環(huán)境變量包括人為環(huán)境變量和自然環(huán)境變量兩大類。其中,人為環(huán)境變量通過距道路距離(Distance to road,Dr)、人類足跡指數(shù)(Human footprint index,Hf)進行體現(xiàn)；而自然環(huán)境變量通過距水體距離(Distance to water body,Dw)、土地覆蓋類型(Land cover type,Lc)、亮度指數(shù)(Brightness index,Bi)、濕度指數(shù)(Wetness index,Wi)和綠度指數(shù)(Greenness index,Gi)等進行體現(xiàn)。

其中,人為環(huán)境變量中的道路矢量數(shù)據(jù)、人類足跡指數(shù)數(shù)據(jù)來源于哥倫比亞大學(xué)社會經(jīng)濟數(shù)據(jù)和應(yīng)用中心數(shù)據(jù)庫(http://sedac.ciesin.columbia.edu/data)的Global Roads Open Access Data Set(Version 1)和Last of the Wild Data Version 3模塊。

水體、土地覆蓋類型數(shù)據(jù)來自國家綜合地球觀測數(shù)據(jù)共享平臺(http://www.chinageoss.org/dsp/home/index.jsp)中空間分辨率為30m×30m的Landsat8和國產(chǎn)高分衛(wèi)星,下載于LandCover_China2017數(shù)據(jù)集,通過ArcGIS 10.6平臺提取土地覆蓋類型柵格數(shù)據(jù)中的水體,并用空間分析工具(Spatial Analyst Tools)對提取的道路和水體的矢量數(shù)據(jù)通過歐式距離分析方法進行柵格化處理。

亮度指數(shù)、濕度指數(shù)和綠度指數(shù)數(shù)據(jù)來自美國地質(zhì)勘探局官方網(wǎng)站(https://earthexplorer.usgs.gov/)上Landsat 8 OLI/TIRS C1 Level-1數(shù)據(jù)集中17景云量低的多波段衛(wèi)星影像數(shù)據(jù)(2018/9/26,LC81200262018269LGN00,LC81200272018269LGN00,LC81200282018269LGN00,LC81200292018269LGN00；2018/9/19,LC81190272018262LGN00；2018/8/18,LC81190262018230LGN00；2018/8/2,LC81190292018214LGN00,LC81190302018214LGN00；2018/7/1,LC81190282018182LGN00；2018/6/24,LC81180262018175LGN00,LC81180302018175LGN00；2018/6/2,LC81170282018152LGN00；2018/6/1,LC81170292018152LGN00,LC81170272018152LGN00；2017/9/25,LC81180272017268LGN00,LC81180282017268LGN00, LC81180292017268LGN00),利用ENVI 5.3對松嫩平原邊界裁剪后的拼接矢量數(shù)據(jù)進行纓帽變換(Kauth-Thomas Transformation),提取亮度指數(shù)、綠度指數(shù)、濕度指數(shù),這些分別反映了地物目標(biāo)總輻射能量水平的亮度信息、影像中植被覆蓋葉面積指數(shù)和生物量等綠色植被、影像的濕度信息[30]。

2.3 模型構(gòu)建與檢驗

用ArcGIS 10.6將所有柵格數(shù)據(jù)變量的行列數(shù)、像元大小進行一致化處理,轉(zhuǎn)換為模型所需的ASC格式,導(dǎo)入Maxent 3.3.3k軟件中[18,24-25,27,31]。在模型設(shè)置中,隨機選取70%的蒼鷺棲息位點用于模型建立,其余30%的棲息位點用于模型驗證。在模型檢驗過程中,利用刀切法(Jackknife test)AUC值對模型進行評價。AUC值為模型自帶的受試者工作特征曲線ROC(Receiver Operator Characteristic)下的面積(Area Under Curve,簡稱AUC)。AUC值越大,表示環(huán)境變量與模型的相關(guān)性越大,模型的預(yù)測結(jié)果就越準(zhǔn)確：當(dāng)AUC值為介于0.5—0.6之間時,意味著所建模型的預(yù)測效果不及格；當(dāng)AUC值為介于0.6—0.7之間時,意味著所建模型的預(yù)測效果較差；當(dāng)AUC值為介于0.7—0.8之間時,意味著所建模型的預(yù)測效果一般；當(dāng)AUC值為介于0.8—0.9之間時,意味著所建模型的預(yù)測效果良好；當(dāng)AUC值為介于0.9—1.0之間時,意味著所建模型的預(yù)測效果優(yōu)秀。

2.4 適宜棲息地分布

模型預(yù)測的棲息地適宜性分級標(biāo)準(zhǔn)值為0—1,值為0意味著棲息地適宜性最低,值為1意味著棲息地適宜性最高。本研究中,將值介于0.0—0.2之間的棲息地定義為不適宜棲息地,將值介于0.2—0.5之間的棲息地定義為低適宜棲息地,將值介于0.5—0.7之間的棲息地定義為較適宜棲息地,將值介于0.7—1.0之間的棲息地定義為適宜棲息地。通過ArcGIS 10.6將模型預(yù)測的結(jié)果進行重分類,用空間分析工具中的區(qū)域分析提取松嫩平原55個縣區(qū)和6個國家級自然保護區(qū)內(nèi)蒼鷺秋季適宜棲息地面積,分析和統(tǒng)計各縣區(qū)的分布情況及受保護情況。

3 結(jié)果與分析

3.1 Maxent模型預(yù)測結(jié)果檢測

ROC曲線評價結(jié)果顯示(圖2)：訓(xùn)練數(shù)據(jù)集(Training data)AUC值為0.946,測試數(shù)據(jù)集(Test data)AUC值為0.933,預(yù)測結(jié)果達到優(yōu)秀水平,表明建立Maxent模型對松嫩平原蒼鷺秋季棲息地適宜性預(yù)測的結(jié)果可信度高。

圖2 蒼鷺秋季棲息地適宜性評價模型的ROC曲線

3.2 蒼鷺秋季棲息地環(huán)境變量分析

從刀切法測定的蒼鷺秋季棲息地變量的貢獻率可以看出(圖3)：距水體距離(Dw)和綠度指數(shù)(Gi)是影響蒼鷺秋季棲息地利用的重要環(huán)境變量,亮度指數(shù)(Bi)、人類足跡指數(shù)(Hf)、土地覆蓋類型(Lc)、濕度指數(shù)(Wi)、距離道路的距離(Dr)等環(huán)境變量對蒼鷺秋季棲息地利用的影響較小。

圖3 蒼鷺秋季棲息地環(huán)境變量刀切法檢驗

3.3 蒼鷺秋季棲息地適宜性分布

松嫩平原蒼鷺秋季棲息地適宜性分布結(jié)果顯示(圖4)：蒼鷺在松嫩平原秋季的適宜棲息地面積為2761.06 km2,僅占研究區(qū)總面積的1.24%；較適宜棲息地7049.57 km2,占比為3.16%；低適宜棲息地36125.15 km2,占比16.19%；不適宜棲息地177129.60 km2,占比79.41%。

圖4 松嫩平原蒼鷺秋季棲息地適宜分布

其中,適宜棲息地在松嫩平原9個地/州的分布各不相同,大慶、白城、齊齊哈爾3個行政區(qū)是主要的分布區(qū)(占適宜棲息地面積的62.78%),其次是哈爾濱和綏化(20.23%),黑河、扶余、長春、四平區(qū)域內(nèi)的適宜棲息地面積較少(16.99%)。適宜棲息地在松嫩平原55個縣/區(qū)的分布也各不相同,杜爾伯特蒙古族自治縣區(qū)域內(nèi)適宜棲息地面積最多(429.90km2,占研究區(qū)適宜棲息地總面積的15.57%),其次是鎮(zhèn)賚縣(334.92km2,12.13%)和肇源縣(185.54km2,6.72%),五常縣、德都縣、大安市區(qū)域內(nèi)蒼鷺適宜棲息地面積均超過100 km2,這些是蒼鷺秋季適宜棲息地面積較多的行政區(qū)域。

3.4 蒼鷺秋季適宜棲息地的保護現(xiàn)狀

由上述分析可以看出,松嫩平原蒼鷺秋季適宜棲息地面積(圖5)的15.79%受到保護,分布在6個國家級自然保護區(qū)內(nèi),分別為莫莫格保護區(qū)(285.62 km2,占秋季適宜棲息地面積的10.34%)、扎龍保護區(qū)(95.80 km2,3.47%)、向海保護區(qū)(18.57 km2,0.67%)、查干湖保護區(qū)(14.89 km2,0.54%)、大布蘇保護區(qū)(11.26 km2,0.41%)、烏裕爾河保護區(qū)(9.95 km2,0.36%)；而五大連池山口湖水庫、綏化東北部、五?？h龍鳳山國家森林公園是處于未保護狀態(tài)的、面積連續(xù)分布較大的區(qū)域。

圖5 蒼鷺適宜棲息地的保護性分布及主要潛在棲息地

4 討論

4.1 蒼鷺秋季棲息地環(huán)境變量

關(guān)于蒼鷺棲息地環(huán)境特征方面的分析,僅有7篇零星的研究報道,分別分散在中國、俄羅斯、日本、波蘭、法國、希臘,在蒼鷺最適營巢生境、蒼鷺生境與仔魚群落關(guān)系、生境特征和不同空間尺度下蒼鷺繁殖力、蒼鷺對村園區(qū)域的季節(jié)性利用、水文網(wǎng)絡(luò)特征與蒼鷺群體大小及空間分布的關(guān)系、針對棲息地蒼鷺形成的覓食策略等方面進行了分析研究[11-12,32-36]。

本文關(guān)于蒼鷺秋季棲息地環(huán)境變量方面的研究,Maxent模型得出影響蒼鷺秋季棲息地的主要環(huán)境變量是距水體距離和綠度指數(shù)。這與蒼鷺生活史期間的食物需求和行為相符合[11-12,35],蒼鷺主要捕食魚類,是典型的濕地鳥類,水域是其偏好的主要生境類型；同時,該研究結(jié)果也與同域分布的其他大型水鳥如白琵鷺(Platalealeucorodia)、丹頂鶴的結(jié)果一致[27,37],蘆葦塘和開闊水面及綠度是其偏好的生境及主要環(huán)境特征。這種既配置有充足食物的水體、也配置有大尺度開闊空間而利于隱蔽警戒逃逸的環(huán)境,不僅能保證蒼鷺棲息所需要的的食物、水、隱蔽等生境需求,也能為其提供警戒、逃逸等安全的空間保障。

4.2 蒼鷺秋季棲息地適宜分布

關(guān)于松嫩平原水禽棲息地適宜性評價方面,已有學(xué)者通過目視解譯、監(jiān)督分類、層次分析等方法結(jié)合遙感影像和GIS系統(tǒng)對松嫩平原水禽棲息地進行了適宜性分析[38-39],也有學(xué)者通過最大熵模型對松嫩平原丹頂鶴秋季生境適宜性進行了評價,分析了適宜性分布的環(huán)境影響因素、適宜性空間分布及景觀格局[40]。本文蒼鷺秋季棲息地適宜性分布的內(nèi)容與丹頂鶴有些類似。

本文分析的蒼鷺秋季適宜棲息地在松嫩平原不同區(qū)域的分布與受保護情況各不相同。其中,西部區(qū)域分布的適宜棲息地較多,集中于嫩江中下游地區(qū),這也是松嫩平原生物多樣性保護優(yōu)先區(qū),受保護性和受關(guān)注度也較高,連片面積較大的區(qū)域主要分布在莫莫格保護區(qū)(占適宜棲息地面積的10.34%)、扎龍保護區(qū)(3.47%)、向海保護區(qū)(0.67%)、查干湖保護區(qū)(0.54%)、大布蘇保護區(qū)(0.41%)、烏裕爾河保護區(qū)(0.36%)等國家級自然保護區(qū)中。

其中,莫莫格保護區(qū)內(nèi)分布的適宜棲息地面積最大、連通性最好。保護好莫莫格保護區(qū)蒼鷺的棲息地對于松嫩平原的蒼鷺種群而言至關(guān)重要,應(yīng)該進一步摸清莫莫格保護區(qū)內(nèi)蒼鷺的種群數(shù)量和時空分布,這不僅能提前預(yù)警蒼鷺種群的動態(tài),更能有助于深入認(rèn)識蒼鷺棲息的濕地資源質(zhì)量情況,利于進而開展及時的資源管護工作。

扎龍保護區(qū)蒼鷺適宜棲息地相對莫莫格保護區(qū)較分散。關(guān)于扎龍保護區(qū)的蒼鷺種群與分布,零星的研究記載顯示出,不僅種群數(shù)量較歷史上有所下降,曾經(jīng)常見的成千上萬只群體已下降為不足千只,甚至幾百只也很難見,而且分布區(qū)斑塊也變少,由原來的成片相連變?yōu)槟壳暗陌邏K狀分布[12,14]。而土地覆蓋類型方面的研究表明,扎龍保護區(qū)內(nèi)水體發(fā)生了明顯的空間分布轉(zhuǎn)移[41],隨之而來的環(huán)境特征等信息也必然會發(fā)生變化,這是否是蒼鷺群體發(fā)生時空變化的根本原因,非常有必要進行深入的大尺度相關(guān)性分析。

向海保護區(qū)蒼鷺適宜棲息地分布相對其他保護區(qū)較少。眾所周知,由于降水量等自然驅(qū)動力和農(nóng)業(yè)開墾、過度放牧等人為驅(qū)動力的影響,向海保護區(qū)的濕地面積呈現(xiàn)出年度間波動較大的下降趨勢,時空均發(fā)生了變化、濕地破碎化增加[42],這些因素極大程度地限制了以濕地資源為生活史需求的各種鳥類的時空分布,蒼鷺也不例外。

松嫩平原東部區(qū)域面積較大的蒼鷺適宜棲息地主要分布在五大連池市的山口湖水庫、綏化地區(qū)東北部的濕地、五常龍鳳山國家森林公園及阿什河流域,此外還有一些零星的適宜棲息地分布在松嫩平原的中東部地區(qū),所有這些區(qū)域均缺少蒼鷺的分布信息報道和管護關(guān)注,如何對這些分布小區(qū)進行科學(xué)管理和積極的必要關(guān)注,對蒼鷺種群的遷徙管理極為必要,應(yīng)加強對這些區(qū)域的關(guān)注和監(jiān)測。

4.3 GPS/GSM技術(shù)

GSM是全球移動通信系統(tǒng)(Global System for Mobile communications)的簡稱,是數(shù)字制式移動通信(2G)網(wǎng)絡(luò)的主導(dǎo)技術(shù),其與GPS結(jié)合,用以傳輸目標(biāo)動物目標(biāo)位置的系列信息,具有信號覆蓋范圍廣、定位精度高、傳輸效率高的特點,這彌補了原來僅通過衛(wèi)星進行信號傳播、導(dǎo)致部分地區(qū)信號弱的缺點。結(jié)合GPS/GSM技術(shù)和Manent模型對野生動物棲息地適宜性、空間擴散、停留規(guī)律等大尺度空間生態(tài)學(xué)進行研究,這種研究融合了衛(wèi)星跟蹤器采集的最為精準(zhǔn)的GPS分布位點、衛(wèi)星拍攝的最為真實的遙感影像、不受樣本點數(shù)量限制的Maxent模型,相對于傳統(tǒng)方法,這種耦合多種先進方法的科學(xué)研究提供的數(shù)據(jù)更為精準(zhǔn)、分析的結(jié)果更為科學(xué),其研究結(jié)果不僅有利于評價瀕危物種棲息地適宜性、發(fā)現(xiàn)潛在棲息地、制定科學(xué)的保護對策,更有助于更為深入地了解瀕危物種的生物學(xué)規(guī)律,這已經(jīng)成為眾多學(xué)者認(rèn)可的最先進的研究技術(shù)[43-45]。本文用該技術(shù)研究了松嫩平原蒼鷺秋季棲息地的適宜性分布,更為深入地掌握了松嫩平原秋季蒼鷺適宜棲息地的空間分布及受保護情況,將能推動我們制定進一步地科學(xué)保護對策,更能促進未受保護的小區(qū)采取相應(yīng)的保護措施和科學(xué)規(guī)劃。