遷徙蒼鷺繁殖期及種群發(fā)展趨勢預測方法探討

劉雙

摘要：遷徙蒼鷺通常是指春季遷來東北等地區(qū)繁殖，秋季遷離繁殖地到南方越冬的蒼鷺，在朝陽為夏候鳥。筆者于2009～2019年，結合“遼寧最大蒼鷺種群生活史研究”項目，對遷徙型蒼鷺繁殖期及種群發(fā)展趨勢預測方法進行了探討。

關鍵詞：蒼鷺；繁殖期；種群；預測

研究地位于朝陽縣七道嶺鎮(zhèn)蘇家營子村，外界稱為蒼鷺村，距朝陽市區(qū)46公里，蒼鷺種群在該地棲息繁殖已達100多年，其間經(jīng)過兩次本地短距離遷移。蘇家營子村四面低山環(huán)繞，植被豐富，植物群落多樣性和物種多樣性較高，居民房前屋后及河灘楊柳科樹木高大。良圖溝河從村北流入，沿村西山根彎肘型向東南流淌，河灘兩側多為楊柳構成的防護林。這條河流是遼西地區(qū)少見的長流水河，從未干枯過，在干旱年份，也只是流量減小、河道變窄而已。該河水質優(yōu)良，常年保持在Ⅰ～Ⅱ級標準，河中魚蝦蛙等資源豐富。

2.1 蒼鷺種群生活史各階段發(fā)生時間預測方法

蒼鷺種群生活史各階段發(fā)生時間預測采用期距法，主要利用當?shù)胤e累多年或多次有關蟲害發(fā)生規(guī)律的資料分析總結當?shù)啬澈οx任何兩種現(xiàn)象（如蟲態(tài)）之間的時間間隔，這種有規(guī)律的帶必然性的時間間隔叫作“期距”。整理多年或多次的平均期距和標準差，對某害蟲下一個現(xiàn)象的日期進行預測。預測公式： f=hi+（xi±sx），式中：f為預測時期出現(xiàn)日期，hi為預測時期前某個時期的實際日期，xi為理論期距值，sx為標準差。

2.2 蒼鷺種群各階段的數(shù)量發(fā)展趨勢預測方法

采用2009～2019年調查數(shù)據(jù)，應用“Microsoft Office Excel2003”做回歸分析。打開“Microsoft Excel”，點擊“圖標向導”，選擇成對的自變量x和因變量y數(shù)據(jù)列，點擊“x、y散點圖”。點擊“添加趨勢線”，在“趨勢預測/回歸分析類型”的線性、對數(shù)、多項式、乘冪、指數(shù)等類型中，采用相關關系系數(shù)r2值大于0.95的方程式作為蒼鷺種群發(fā)展趨勢預測回歸模型。

3.1 蒼鷺種群生活史各階段發(fā)生時間預測

根據(jù)2009～2019年調查的蒼鷺從遷來始見時間到南遷末期時間，各生活時段點時間（見表1），計算出以下生活時段間的期距值和標準差（見表2）。

蘇家營子蒼鷺繁殖群2013年遷來始見日期為3月4日，即“hi”為3月4日，“xi±sx”為35±4=31～39天，在3月4日基礎上加31天和39天，得2013年蒼鷺產(chǎn)卵始見期為4月5～12日，實際調查為4月8日，在預測區(qū)間值內(nèi)。

3.2 蒼鷺種群各階段的數(shù)量發(fā)展趨勢預測

根據(jù)2009～2019年調查數(shù)據(jù)，應用“Excel2013”建立數(shù)學回歸預測模型。

蒼鷺遷來數(shù)量預測回歸模型：y=0.001x2+0.9541x+188.51 r2=0.9874。式中：x為蒼鷺遷來前的鳥巢調查數(shù)據(jù)（個）；y為蒼鷺遷來數(shù)量預測值（只）；r2為相關系數(shù)。

蒼鷺幼鳥離巢數(shù)量預測回歸模型：y=0.0025x2+2.3044x+449.39 r2=0.9898。式中：x為蒼鷺遷來前的鳥巢調查數(shù)據(jù)（個）；y為蒼鷺幼鳥離巢數(shù)量預測值（只）；r2為相關系數(shù)。

幼鳥保存數(shù)量預測回歸模型：y=0.8315x+1.1701 r2=0.9956。式中：x為蒼鷺遷來數(shù)量（只）；y為蒼鷺幼鳥離巢數(shù)量預測值（只）；r2為相關系數(shù)。

蒼鷺南遷數(shù)量預測回歸模型：y=1.88568x+1.1624 r2=0.9912。式中：x為蒼鷺遷來數(shù)量（只）；y為蒼鷺南遷數(shù)量預測值（只）；r2為相關系數(shù)。

預測模型檢驗舉例：蘇家營子蒼鷺繁殖群遷來前的鳥巢調查數(shù)據(jù)為2 5 8個，預測蒼鷺遷來數(shù)量，應用模型y=0.001x2+0.9541x+188.51計算，得y=66.564+246.157+188.51=501 .231（只），實際調查值為498只，相差3只，在允許范圍之內(nèi)。

應用期距預測法對蒼鷺種群某階段發(fā)生時間進行預測和建立數(shù)學回歸模型對蒼鷺種群各階段的數(shù)量發(fā)展趨勢預測，均未見相關報道。經(jīng)檢驗，應用期距法預測時間與調查時間相符；應用回歸模型對蒼鷺種群數(shù)量預測，預測值在允許范圍之內(nèi)。