• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    人工紅松林步行蟲(Coleoptera: Carabidae)群落物種共存格局動態(tài)分析

    2016-11-29 03:22:54高梅香李景科張雪萍朱紀元
    生態(tài)學報 2016年20期
    關鍵詞:紅松林隨機性格局

    高梅香,李景科,石 昊,張雪萍,*,朱紀元

    1 哈爾濱師范大學地理科學學院,哈爾濱 150025 2 黑龍江省普通高等學校地理環(huán)境遙感監(jiān)測重點實驗室,哈爾濱 150025

    ?

    人工紅松林步行蟲(Coleoptera: Carabidae)群落物種共存格局動態(tài)分析

    高梅香1,2,李景科1,2,石 昊1,2,張雪萍1,2,*,朱紀元1,2

    1 哈爾濱師范大學地理科學學院,哈爾濱 150025 2 黑龍江省普通高等學校地理環(huán)境遙感監(jiān)測重點實驗室,哈爾濱 150025

    揭示群落物種共存格局是群落生態(tài)學研究的重點內(nèi)容之一,零模型的應用極大的促進了群落物種共存格局及其調(diào)控機制的進展,然而針對地下生物群落共存格局動態(tài)特征的研究并不多見。在帽兒山森林生態(tài)站的人工紅松林內(nèi),通過5次調(diào)查取樣基于零模型模擬分析小尺度空間(20 m×20 m)步行蟲群落物種共存格局的動態(tài)特征。結果表明:(1)共捕獲步行蟲20種,2278只個體,其中Carabusbillergimaoershanensis為所有調(diào)查季節(jié)數(shù)量最具優(yōu)勢且分布最廣泛的物種,步行蟲群落結構具明顯的時間變異性;(2)2013年6月步行蟲群落為集群性共存格局,而2014年8、10月為競爭性共存格局,基于目前的零模型指標和法則難以準確揭示其他月份是集群性還是競爭性的共存格局,但所有季節(jié)的群落均表現(xiàn)為明顯的非隨機性共存格局,這些共存格局的發(fā)現(xiàn)并不完全支持Diamond的群落構建機制理論;(3)所有調(diào)查季節(jié)均發(fā)現(xiàn)很少的顯著物種對,基于更嚴格的檢驗表明群落中集群性物種對多于隔離性物種對,那些表現(xiàn)為顯著的非隨機性共存關系的物種對往往是群落內(nèi)數(shù)量較大且分布廣泛的優(yōu)勢和常見物種。表明非隨機性共存格局可能是帽兒山人工紅松林小尺度空間步行蟲群落的常見格局,這種非隨機性格局具一定的短期動態(tài)穩(wěn)定性,但不同季節(jié)這種非隨機性共存格局類型表現(xiàn)不同,群落內(nèi)這些較少的顯著物種對可能對群落物種共存格局具有一定的貢獻。

    共存格局;顯著物種對;零模型;步行蟲;人工紅松林;帽兒山

    群落生態(tài)學研究的一個重要任務就是揭示在一定時空尺度內(nèi),群落是否表現(xiàn)為一種非隨機性共存格局并進一步揭示該群落受到何種機制的調(diào)控,如構建機制[1]、鑲嵌性[2]、中心-衛(wèi)星結構[3]、分隔[4]或者擴散限制。Diamond[1]首先提出鳥類的群落組成受到基于種間競爭的構建機制調(diào)控,他進一步指出種間競爭會阻止那些具有生態(tài)學相似性的物種在群落內(nèi)共存,之后物種共存分析被越來越多的用于檢驗群落構建機制是否存在[5]。然而物種共存格局分析一直是群落生態(tài)學研究的難點[6-7],直到零模型(null model)的應用才使構建機制的檢驗得以廣泛開展。零模型是一種格局生成模型(a pattern-generating model)的統(tǒng)計檢驗,該檢驗基于已知或者想象分布的生態(tài)數(shù)據(jù)的隨機性或取樣的隨機性[8],故意排除一種被檢驗的機制而產(chǎn)生隨機性的群落格局,通過隨機格局和調(diào)查群落格局之間的比較來評估生物間相互作用的調(diào)控作用。零模型模擬探討群落物種共存格局時,一個重要的問題就是識別顯著的物種對,即識別那些在群落內(nèi)表現(xiàn)為顯著的正的、負的或者隨機性的物種對,這有助于在群落水平和物種水平共同揭示群落物種共存格局?;诹隳P偷娜郝湮锓N共存格局分析及顯著物種對識別被應用到地下生態(tài)系統(tǒng),包括蚯蚓[9-10]、跳蟲[11-12]、螨類[13]等,但是對地下生物群落物種共存格局時間動態(tài)特征的研究非常少見。

    步行蟲(Coleoptera: Carabidae)是地表活動動物中廣泛分布的物種,也是北半球節(jié)肢動物中研究較深入的物種[14],對環(huán)境條件變化和干擾比較敏感,常常被用來作為生物指示種[15]或生物多樣性指標[16-17],而且步行蟲是相對的長壽命物種,能夠基于簡便可行的陷阱法進行采集[18]?;诖?步行蟲成為評估微生境選擇和小尺度動態(tài)的理想有機體。很多研究探討了步行蟲物種多樣性和豐富度時空變異[19],卻少有研究揭示步行蟲群落物種共存格局的動態(tài)變化。步行蟲物種共存格局及其動態(tài)變化的研究,是揭示群落物種共存機制的基礎,為以步行蟲為指示種的環(huán)境保護和生態(tài)管理提供基礎數(shù)據(jù)。本研究在帽兒山森林生態(tài)系統(tǒng)研究站的人工紅松林內(nèi),通過2013年6月至2014年10月的5次野外調(diào)查,研究小尺度(101—103m)[20]范圍內(nèi)步行蟲科物種共存格局的動態(tài)特征,嘗試性的回答以下問題:(1)步行蟲群落是否為非隨機性的共存格局,(2)這種共存格局是否具有時間穩(wěn)定性,(3)群落內(nèi)是否存在較多顯著的共存物種對,(4)這些共存物種對是否對群落水平物種共存格局有重要貢獻。

    1 研究地區(qū)與研究方法

    1.1 研究地區(qū)

    本實驗設置在帽兒山森林生態(tài)系統(tǒng)研究站內(nèi)(127°30′—127°34′ E,45°20′—45°25′ N)。研究區(qū)位于黑龍江省東北部,平均海拔300m,平均坡度10°—15°,地帶性土壤為暗棕色森林土。為四季分明的大陸性季風氣候,夏季溫暖濕潤而冬季寒冷干燥,年均溫約3.1℃[21],年均降水量為629 mm。研究區(qū)植被屬于長白植物區(qū)系,實驗樣地設置在林齡約為60a的人工紅松林內(nèi),平均坡度12°,西北坡向,平均胸徑12 cm,林分密度3145個/hm2,喬木層以紅松(Pinuskoraiensis)為優(yōu)勢種,同時包括白樺(Betulaplatyplylla)、水曲柳(Fraxinusmandshurica)、楓樺(Betulacostata)、山楊(Populusdavidiana)、胡桃楸(Juglansmandshurica),下木層包括春榆(Ulmuspropiqua)和山梨(Pyrusussuriensis)[22]。

    1.2 樣地設置與調(diào)查方法

    在人工紅松林內(nèi)隨機設置一個20m×20m的樣地,以2m為間隔均分成100個2m×2m的小單元格,在每個小單元格左下角網(wǎng)格交叉點采集樣品(共計100個采樣點)。陷阱法捕獲步行蟲樣品,將誘捕杯(高9 cm,內(nèi)徑6.5 cm)放置在采樣點,杯口與地面齊平,杯口上部離地面約10 cm處支起1個一次性餐盤,以防止凋落物或者穿透雨的進入。2013年6、8、9月內(nèi)置4%福爾馬林,2014年8、10月內(nèi)置飽和NaCl溶液,將陷阱至于野外7天7夜后取回。手撿法分揀步行蟲,將樣品放置于95%醫(yī)用酒精內(nèi),顯微鏡下鑒定種類并計數(shù)。

    1.3數(shù)據(jù)處理分析

    1.3.1 Raunkiaer頻度指數(shù)、方差分析和基于采樣點的稀釋度曲線

    依據(jù)Raunkiaer[23]公式計算物種頻度,劃分標準為:1%—20%為A級,21%—40%為B級,41%—60%為C級,61%—80%為D級,81%—100%為E級。單因素方差分析(one-way ANOVA)比較相同誘捕液條件下,步行蟲物種數(shù)量和個體數(shù)量差異的顯著性?;诓蓸狱c的稀釋度曲線(rarefaction analysis)[24-25]檢驗研究尺度內(nèi)物種-強度調(diào)查關系(species-trapping-effort relationship)。數(shù)量百分比劃分類群,依據(jù)各物種數(shù)量占該群落總捕獲數(shù)量百分比分別劃歸為:>10%為優(yōu)勢物種,1%—10%為常見物種,<1%為稀有物種。Raunkiaer頻度和數(shù)量百分比在Microsoft Excel 2007中計算,單因素方差分析和多重檢驗使用R軟件中的“aov”和“duncan.test”函數(shù)實現(xiàn)[26],基于采樣點的稀釋度曲線使用R軟件vegan軟件包中的“specaccum”函數(shù)實現(xiàn)。

    1.3.2 基于零模型的物種共存格局分析

    采用零模型模擬分析步行蟲群落是否表現(xiàn)出非隨機的共存格局[27-28]。首先將步行蟲豐富度矩陣轉換為“有-無“矩陣(presence-absence matrix)[29],該矩陣每行代表一個步行蟲科物種,每列代表一個不同的采樣點,“0”值代表該物種在特定采樣點不出現(xiàn),“1”值代表該物種在特定采樣點出現(xiàn)[30]。選擇C-score和V-ratio指標量化步行蟲群落的共存格局,其中每個指標都是一個單獨的數(shù)值來量化整個“有-無”矩陣的格局。

    (1)C-score指標 Stone和Roberts[31]引入C-score指標作為量化物種共存程度的一個指標,該指標量化每個物種對能夠被發(fā)現(xiàn)的棋盤單元(checkerboard unit, CU)的數(shù)量。任何物種對的棋盤單元的數(shù)量可以由公式獲得:CU=(Ri-S)(Rj-S),式中Ri和Rj分別為物種i和j的總行數(shù),S是被兩物種共同占據(jù)的采樣點的數(shù)量。C-score是所有獨特物種對(unique pairs of species)計算的CU的均值。C-score的值越大,共存的物種對就越少[32-33],如果C-score的值顯著的比預期機會(expected by chance, EBC)大,那么該集群被認為是競爭性/隔離性群落結構[27],說明種間競爭可能對群落物種共存具有重要的調(diào)控作用;反之,該集群被認為是集群性群落結構,表明種間競爭對群落物種共存的調(diào)控作用可能不明顯。

    (2)V-ratio指標 Schluter[34]提出將V-ratio作為量化群落結構的一個指標,該值是列總和方差與行總和方差的比值,測量每個采樣點物種的數(shù)量的變異性。在一個競爭性結構的群落中,實測的V-ratio應該明顯的比預期機會(EBC)??;而在一個集群性的群落中,實測的V-ratio應該明顯的比預期機會(EBC)大。

    基于零模型的共存格局分析對物種共存頻率的變異敏感,Gotelli建議零模型中的行數(shù)應被限定[27],3種零模型法則被推薦用來模擬物種共存格局,這些法則固定行總和且較少產(chǎn)生Ⅰ類錯誤。本文選擇FF(Fixed row-Fixed column,固定行-固定列), FE (Fixed row-Equiprobable column, 固定行-等概率列)和FP (Fixed row-Probability column, 固定行-或然率列)法則來計算C-score,使用FE和FP算法來計算V-ratio[27-28,35]。

    標準效應量(Standardized Effect Size, SES)用來比較結果的顯著性。SES評估實測指標距離平均模擬指標的標準差[5],對于C-score來說,SES大于2.0表明步行蟲群落共存格局為非隨機性的物種間隔離,說明為競爭性的群落結構;SES小于-2.0表明步行蟲群落共存格局為非隨機性物種間集群,說明群落表現(xiàn)為顯著的非競爭性結構。對于V-ratio指標,SES大于2.0說明共存格局表現(xiàn)為非隨機性物種間集群,而SES小于-2.0說明共存格局表現(xiàn)為非隨機性物種間隔離。當群落內(nèi)產(chǎn)生上百個物種對的時候P值應該被校正[35],本文采用控制發(fā)現(xiàn)錯誤率(the false discovery rate, FDR)對顯著性水平進行校正[36-37]。C-sore和V-ratio均在Ecosim 7.0軟件的“Co-occurrence”模塊計算[38],隨機模擬運行50000次;P值的校正使用R軟件的“p.adjust”函數(shù)實現(xiàn)。

    1.3.3 基于零模型的顯著物種對

    基于零模型的顯著物種對分析會產(chǎn)生大量的物種對,對于具有n個物種的步行蟲矩陣來說,共有[(n)(n-1)/2]個物種對,而這其中的許多物種對可能彼此不是生物學或者統(tǒng)計學相互獨立的[35],在分析物種對間相互作用關系時會在很大程度上增加I類錯誤的產(chǎn)生[35]。本文通過計算每個物種對的C-score值來確定特定物種對之間的關系,該計算過程中C-score被重采樣為0—1之間,C-score=0代表完全重疊,C-score=1代表完全隔離[35]。然后采用Gotelli和Ulrich[35]提出的4種準則來檢驗這些物種對之間為何種共存關系,即集群性、隔離性或隨機性共存物種對。其中CL (the confidence limit criterion)是最簡單和最寬泛的一種檢驗準則,Gotelli和Ulrich[35]同時推薦更嚴格的檢驗準則,即BY (after sequential Bonferroni correction)、BM (empirical Bayes mean based criterion) 和BCL (empirical Bayes confidence limits based criterion)。關于4種檢驗準則的理論基礎和計算過程,請參閱文獻Gotelli和Ulrich[35]。

    最后,計算那些被證明存在顯著相互作用關系的物種對C-score值的標準效應量(SES),SES大于2.0表明物種對為顯著的負相互作用關系(隔離性關系),說明兩個物種間存在顯著的競爭;SES小于-2.0表明物種對為顯著的正相互作用關系,說明物種對不存在顯著的種間相互競爭。顯著物種對的識別在Pairwise軟件中計算[39],隨機模擬運行1 000次。

    1.3.4 PiankaOik生態(tài)位重疊指數(shù)

    和營養(yǎng)資源類似,空間也可以被作為一種資源被物種競爭[40]。Pianka生態(tài)位重疊指數(shù)[41]被用來測度兩個步行蟲科物種利用空間資源而互相重疊的情況:

    (1)

    式中,Oik為物種i和物種k的生態(tài)位重疊度值,Pij和Pkj為物種i和物種k在空間資源梯度級j的數(shù)量特征,r為空間樣點數(shù)。PiankaOik生態(tài)位重疊指數(shù)的值介于0和1之間,該指數(shù)在Ecosim 7.0軟件的“niche overlap”模塊中計算。

    2 結果

    2.1 群落結構特征

    共捕獲步行蟲科20種,2278只個體。2013年6、7、9月均捕獲較多的物種數(shù)量和個體數(shù)量,而2014年10月捕獲最少的物種數(shù)量和個體數(shù)量。2013年各個月份間捕獲類群數(shù)量沒有顯著差異,而個體數(shù)量有顯著差異(P<0.01);2014年2個月份間捕獲類群數(shù)量和個體數(shù)量都有顯著差異(P<0.01)(表1)。群落物種稀釋度曲線表明,在每個調(diào)查季節(jié)內(nèi)步行蟲物種數(shù)量隨著調(diào)查樣點的增加而增加,說明需要在研究尺度內(nèi)適當增加取樣量。

    Carabusbillergimaoershanensis在所有調(diào)查季節(jié)均為優(yōu)勢物種,2013年6月捕獲量為總捕獲量的84.56%,2013年7月和2014年10月所占比例均為60%以上,2013年9月和2014年8月也達到了30%以上。Carabusbillergimaoershanensis為出現(xiàn)頻度最高的物種,2013年6月Raunkiaer頻度為D級,2014年10月為B級,其他月份均為C級。所有群落70%左右的步行蟲科物種的頻度指數(shù)均屬于A級,說明大部分物種在樣地內(nèi)出現(xiàn)頻度不高。2013年6、7、9月和2014年8、10月群落內(nèi)優(yōu)勢、常見和稀有物種的數(shù)量分別為1、4、10;2、3、10;3、4、6;3、5、6;1、5、0。僅2014年10月群落內(nèi)沒有發(fā)現(xiàn)稀有類群,其他季節(jié)稀有類群較多(表1)。

    2.2 群落物種共存格局

    2013年6月,基于FE法則的C-score值明顯的小于期望值(模擬值),而V-ratio值明顯的大于期望值,表明步行蟲群落為集群性共存格局。2013年7月,基于FE法則的C-score和V-ratio表明步行蟲群落為集群性共存格局;而基于FF、FP法則的C-score和基于FP法則的V-ratio表明群落為隔離性格局。2013年9月,基于FE法則的C-score和V-ratio表明群落為集群性格局,而基于FP法則的C-score和V-ratio表明步行蟲群落為隔離性共存格局。2014年8月基于所有法則的C-score和V-ratio均表明步行蟲群落為隔離性共存格局。2014年10月基于FP法則的C-score和V-ratio均表明步行蟲群落為隔離性共存格局(表2)。

    2.3 顯著物種對識別

    基于CL準則,2013年6月,基于FF法則發(fā)現(xiàn)1個具顯著正相互作用關系的物種對,說明二者表現(xiàn)為集群性共存關系;2013年7月,基于FE和FF法則均發(fā)現(xiàn)1個具顯著正相互作用關系的物種對;2013年9月,基于FE法則發(fā)現(xiàn)3個具顯著正相互作用關系的物種對,而基于FF法則分別發(fā)現(xiàn)1個正的和負的相互作用關系物種對;2014年8月,基于FE和FF法則均發(fā)現(xiàn)3個正的和2個負的相互作用關系物種對;2014年10月未發(fā)現(xiàn)顯著物種對?;诟鼑栏竦腂M, BCL或者BY法則,所有月份均發(fā)現(xiàn)較多的正相互作用關系物種對(表3)。這些發(fā)現(xiàn)的顯著物種對之間的PiankaOik生態(tài)位重疊度指標,往往高于群落內(nèi)其他物種對之間的重疊度指標(表4)。

    3 討論

    本調(diào)查共捕獲步行蟲科20種,2278只個體,人工紅松林內(nèi)步行蟲群落存在明顯時間變異性,2013年調(diào)查季節(jié)步行蟲個體數(shù)量彼此間均存在顯著的差異,2014年2個月份類群數(shù)量和個體數(shù)量均有顯著差異。2014年10月為較寒冷的調(diào)查季節(jié),食物資源減少、地表溫濕度等適宜性明顯降低,導致群落的物種數(shù)量和個體數(shù)量明顯下降,稀有物種在樣地內(nèi)消失,群落內(nèi)僅發(fā)現(xiàn)優(yōu)勢和常見物種而沒有稀有物種,這與其他生長季節(jié)明顯不同,其他季節(jié)群落內(nèi)稀有物種的物種數(shù)量多于優(yōu)勢和常見物種。群落內(nèi)物種的Raunkiar頻度指數(shù)多為A級,說明大部分物種不僅數(shù)量稀少,而且僅僅分布在樣地內(nèi)有限的空間內(nèi)。Carabusbillergimaoershanensis是所有季節(jié)最具優(yōu)勢的物種,數(shù)量最豐富且在樣地內(nèi)分布最廣泛。

    表2 基于零模型的步行蟲群落物種共存格局

    表3 基于C-score值的通過4種檢驗準則的步行蟲群落顯著物種對

    物種代碼詳見表1;FE: 固定行-等概率列法則Fixed-Equiprobable algorithm;FF: 固定行-固定列法則Fixed-Fixed algorithm

    表4 2013年6月份步行蟲群落物種的Pianka Oik生態(tài)位重疊度指標

    物種代碼詳見表1

    步行蟲群落表現(xiàn)為非隨機性的共存格局,這種非隨機性的共存格局表現(xiàn)為一定的時間穩(wěn)定性?;贔P法則的V-ratio指標和基于所有法則的C-score和V-ratio指標分別表明,2014年8月和10月步行蟲群落為競爭性-隔離性共存格局,暗示種間競爭在群落共存中可能起到重要的調(diào)控作用,這個結論支持了Diamond群落構建機制理論[1]?;贔E法則的C-score指標表明, 2013年6月步行蟲群落為集群性共存格局,即種間競爭可能不是群落構建的主要調(diào)控機制,此結論并沒有為Diamond群落構建機制理論[1]提供有效證據(jù)?;诓煌▌t的C-score和V-ratio指標則表明,2013年7月和9月步行蟲群落表現(xiàn)為集群性或競爭性共存格局,這些群落存在的大量隨機性物種對會產(chǎn)生稀釋效應(dilution effects)[35, 42],可能導致我們采用的零模型指標和方法難以確定這兩個月份群落非隨機共存格局類型。另外,零模型是基于物種特性或生態(tài)學性狀而隨機構建共存格局的[43- 45],雖然被廣泛的應用在土壤動物群落生態(tài)學研究中[7, 9, 11, 13],但是任何一種方法的應用都有其局限性[28],有研究認為零模型使用的交換法則和Z值(Z-score),容易產(chǎn)生法則和隨后的膨脹I類錯誤水平(inflated type I error levles)帶來的偏差[8, 46]?;诖?在零模型模擬過程中增加隨機模擬運行次數(shù)(5 0000次)以避免這些可能誤差源的產(chǎn)生。雖然依據(jù)目前的零模型指標和方法,難以明確的揭示2013年7、9月步行蟲群落物種共存格局是集群性還是隔離性,但可以在一定程度上說明這兩個月份群落共存格局均為非隨機性格局。綜上,調(diào)查季節(jié)內(nèi)步行蟲群落表現(xiàn)為較穩(wěn)定的非隨機性格局,推測研究區(qū)步行蟲群落非隨機性的共存格局可能是普遍存在的,但非隨機性的表現(xiàn)形式(集群性和隔離性)可能具有季節(jié)性差別。

    本調(diào)查僅發(fā)現(xiàn)較少的顯著物種對,但這些步行蟲物種之間的非隨機性相互作用對整個群落的非隨機性共存格局可能有一定貢獻。在表現(xiàn)為集群性結構的群落中(2013年6月),發(fā)現(xiàn)集群性物種對而沒有發(fā)現(xiàn)隔離性物種對;在表現(xiàn)為隔離性結構的群落中(2014年8月),僅發(fā)現(xiàn)2個隔離性物種對卻發(fā)現(xiàn)3個集群性物種對。本研究的發(fā)現(xiàn)并不完全支持先前的假設,即如果群落內(nèi)大部分物種對都表現(xiàn)為顯著的集群性關系,那么該群落表現(xiàn)為一個集群性的群落結構[42, 47]。但這并不能完全否定這些顯著物種對在群落構建中可能會起到重要的作用。2013年7月份的步行蟲群落,僅發(fā)現(xiàn)1個(“Carabusbillergimaoershanensis-Aulonocarabuscanaliculatushailiniensis”)具正相互作用關系的物種對,僅占該群落物種對數(shù)量的0.95%;這兩個物種均為該群落的優(yōu)勢物種,其個體數(shù)分別占該群落總捕獲量的62.18%和26.78%,兩個物種個體數(shù)合計占該群落總捕獲量的88.96%;且這兩個物種在樣地內(nèi)廣泛出現(xiàn),其Raunkier頻度指數(shù)均為C級;并且兩個物種的PinkaOik生態(tài)位重疊度指數(shù)為0.87,遠遠大于群落內(nèi)其他物種對之間的PinkaOik重疊度指標,表明這兩個物種具有較大的空間生態(tài)位重疊。其他季節(jié)發(fā)現(xiàn)的顯著物種對也多為數(shù)量和分布上占優(yōu)勢的優(yōu)勢和常見物種,在不同季節(jié)表現(xiàn)出不同的競爭性或者集群性共存關系。這些優(yōu)勢和常見物種可能具有相似的環(huán)境喜好,當環(huán)境資源不充足時在局地環(huán)境對食物、空間等資源產(chǎn)生激烈的競爭,會導致較強的物種間競爭排斥現(xiàn)象,當環(huán)境資源比較充足的時候,這種空間生態(tài)位重疊并不一定會導致競爭,可能形成正的集群性的物種間關系[47]。可能恰恰是這些數(shù)量較大、分布廣泛的優(yōu)勢和常見物種間的顯著相互作用,對整個群落的共存格局產(chǎn)生重要的貢獻[35]。但是基于目前的實驗數(shù)據(jù),仍然難以揭示出這些顯著物種對于群落構建起到怎樣的作用,在未來的研究中,需要引入食物網(wǎng)技術和模擬控制實驗,深入分析那些數(shù)量較大且分布廣泛的物種在群落構建中的調(diào)控作用。

    種間競爭通常不被認為是地表甲蟲群落構建的重要調(diào)控機制[48-49],但也有研究表明種間競爭可能導致步行蟲群落出現(xiàn)競爭性共存格局[50]。相關研究證明步行蟲群落的共存格局因物種特性而異,如Ulrich W和Zalewshi M[51]發(fā)現(xiàn)區(qū)域尺度核心步行蟲物種的共存可能受到生物間相互作用和生態(tài)位分化過程的調(diào)控,但是衛(wèi)星步行蟲物種的共存最可能受到隨機擴散的調(diào)控而非種間競爭的調(diào)控,但少有研究揭示小尺度空間群落物種共存格局的時間動態(tài)特征。本研究發(fā)現(xiàn)種間競爭對2014年8月和10月步行蟲的群落構建可能具有調(diào)控作用,但是對2013年6月的群落構建沒有明顯的調(diào)控作用,雖然難以確定2013年7月和9月的非隨機性共存格局的形式,但分析結果表明步行蟲群落的非隨機性格局具有一定的時間穩(wěn)定性。事實上除了生物間相互作用(種間競爭),土壤環(huán)境條件異質(zhì)性和擴散限制[9, 52]等過程也可能導致土壤動物群落的非隨機性共存格局。本實驗空間尺度較小(20m),且人工紅松林的環(huán)境條件相對穩(wěn)定,可能會導致步行蟲群落非隨機性共存格局的相對穩(wěn)定性。本研究沒有考慮生境異質(zhì)性和擴散限制的調(diào)控作用,目前難以明確的揭示步行蟲群落這種非隨機性空間格局穩(wěn)定性的調(diào)控機制,相關問題有待于我們進一步深入研究。

    4 結論

    本研究在黑龍江省帽兒山森林生態(tài)站的人工紅松林內(nèi),研究了小尺度空間步行蟲群落物種共存格局的動態(tài)特征,對本文開端提出的問題做出回答:步行蟲群落具有非隨機性的物種共存格局;這種非隨機性共存格局具有一定的時間穩(wěn)定性,可能是研究區(qū)小尺度空間步行蟲群落的常見格局;群落內(nèi)發(fā)現(xiàn)較少的顯著物種對;這些具有顯著相互作用關系的物種對往往是群落內(nèi)數(shù)量占優(yōu)勢且在樣地分布廣泛的物種,可能對整個群落的非隨機性格局具有一定的貢獻。今后研究中,需要深入揭示這些優(yōu)勢和常見物種在群落共存格局中的調(diào)控作用,更需要探討群落這種非隨機性格局的年度變化及調(diào)控機制。本實驗僅能揭示人工紅松林小尺度空間(20 m×20 m)步行蟲群落的共存格局及其時間動態(tài)性,更大和更小空間有待深入研究。

    致謝:胡媛媛、閆龍、鐘雨晴在野外調(diào)查和樣品制備過程中給予幫助,感謝黑龍江省帽兒山森林生態(tài)系統(tǒng)國家野外科學觀測研究站提供樣地等基礎設施服務和支持,特此致謝。

    [1] Diamond J M. Assembly of Species Communities. Cambridge: Harvard University Press, 1975: 342- 444.

    [2] Patterson B D, Atmar W. Nested subsets and the structure of insular mammalian faunas and archipelagos. Biological Journal of the Linnean Society, 1986, 28(1/2): 65- 82.

    [3] Hanski I. Communities of bumblebees: testing the core-satellite species hypothesis. AnnalesZoologiciFennici, 1982, 19(2): 65- 73.

    [4] Leibold M A, Mikkelson G M. Coherence, species turnover, and boundary clumping: elements of meta-community structure. Oikos, 2002, 97(2): 237- 250.

    [5] Gotelli N J, McCabe D J. Species co-occurrence: a meta-analysis of J. M. diamond′s assembly rules model. Ecology, 2002, 83(8): 2091- 2096.

    [6] Yin X Q, Song B, Dong W H, Xin W D, Wang Y Q. A review on the eco-geography of soil fauna in China. Journal of Geographical Sciences, 2010, 20(3): 333- 346.

    [7] Tiho S, Josens G. Co-occurrence of earthworms in urban surroundings: a null model analysis of community structure. European Journal of Soil Biology, 2007, 43(2): 84- 90.

    [8] Gotelli N J. Research frontiers in null model analysis. Global Ecology and Biogeography, 2001, 10(4): 337- 343.

    [9] Jiménez J J, Deca?ns T, Rossi J P. Soil environmental heterogeneity allows spatial co-occurrence of competitor earthworm species in a gallery forest of the Colombian ‘Llanos’. Oikos, 2012, 121(6): 915- 926.

    [10] Ehouman N M, Tiho S, Dagnogo M. Co-occurrence of earthworms in Lamto savanna: a null model analysis of community structure. European Journal of Soil Biology, 2012, 53: 40- 47.

    [11] Caruso T, Trokhymets V, Bargagli R, Convey P. Biotic interactions as a structuring force in soil communities: evidence from the micro-arthropods of an Antarctic moss model system. Oecologia, 2013, 172(2): 495- 503.

    [12] 高梅香, 何萍, 孫新, 張雪萍, 吳東輝. 環(huán)境篩選、擴散限制和生物間相互作用在溫帶落葉闊葉林土壤跳蟲群落構建中的作用. 科學通報, 2014, 59(24): 2426- 2438.

    [13] Lin L, Gao M X, Liu D, Zhang X P, Wu H T,Wu D H. Co-occurrence patterns of above-ground and below-ground mite communities in farmland of Sanjiang Plain, Northeast China. Chinese Geographical Science, 2014, 24(3): 339- 347.

    [14] L?vei G L, Sunderland K D. Ecology and behavior of ground beetles (Coleoptera: Carabidae). Annual Review of Entomology, 1996, 41: 231- 256.

    [15] Rusch A, Bommarco R, Chiverton P, ?berg S, Wallin H,Wiktelius S, Ekbom B. Response of ground beetle (Coleoptera, Carabidae) communities to changes in agricultural policies in Sweden over two decades. Agriculture, Ecosystems & Environment, 2013, 176: 63- 69.

    [16] Niemel? J, Kotze J, Ashworth A, Brandmayr P, Desender K, New T, Penev L, Samways M, Spence J. The search for common anthropogenic impacts on biodiversity: a global network. Journal of Insect Conservation, 2000, 4(1): 3- 9.

    [17] Koivula M J. Useful model organisms, indicators, or both? Ground beetles (Coleoptera, Carabidae) reflecting environmental conditions. Zookeys, 2011, 100: 287- 317.

    [18] Lindroth C H. Coleoptera: Carabidae. Handbooks for the Identification of British Insects. London: Royal Entomological Society, 1974.

    [19] Günther J, Assmann T. Fluctuations of carabid populations inhabiting an ancient woodland (Coleoptera, Carabidae). Pedobiologia, 2004, 48(2): 159- 164.

    [20] Hortal J, Roura-Pascual N, Sanders N J, Rahbek C. Understanding (insect) species distributions across spatial scales. Ecography, 2010, 33(1): 51- 53.

    [21] 焦振, 王傳寬, 王興昌. 溫帶落葉闊葉林冠層CO2濃度的時空變異. 植物生態(tài)學報, 2011, 35(5): 512- 522.

    [22] 楊金艷, 王傳寬. 東北東部森林生態(tài)系統(tǒng)土壤碳貯量和碳通量. 生態(tài)學報, 2005, 25(11): 2875- 2882.

    [23] Raunkiaer C. The Life Forms of Plants and Statistical Plant Geography.London: Oxford University Press, 1934.

    [24] Si X F, Pimm S L, Russell G J, Ding P. Turnover of breeding bird communities on islands in an inundated lake. Journal of Biogeography, 2014, 41(12): 2283- 2292.

    [25] Gotelli N J, Colwell R K. Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecology Letters, 2001, 4(4): 379- 391.

    [26] 林元震, 陳曉陽. R與ASReml-R統(tǒng)計分析教程. 北京: 中國林業(yè)出版社, 2014.

    [27] Gotelli N J. Null model analysis of species co-occurrence patterns. Ecology, 2000, 81(9): 2606- 2621.

    [28] Gotelli N J, Ulrich W. Statistical challenges in null model analysis. Oikos, 2012, 121(2): 171- 180.

    [29] McCoy E D, Heck Jr K L. Some Observations on the Use of taxonomic similarity in large-scale biogeography. Journal of Biogeography, 1987, 14(1): 79- 87.

    [30] Zhang J, Hao Z Q, Song B, Li B H, Wang X G, Ye J. Fine-scale species co-occurrence patterns in an old-growth temperate forest. Forest Ecology and Management, 2009, 257(10): 2115- 2120.

    [31] Stone L, Roberts A. The checkerboard score and species distributions.Oecologia, 1990, 85(1): 74- 79.

    [32] Haukisalmi V, Henttonen H. Analysing interspecific associations in parasites: alternative methods and effects of sampling heterogeneity. Oecologia, 1998, 116(4): 565- 574.

    [33] Gotelli N J, Arnett A E. Biogeographic effects of red fire ant invasion. Ecology Letters, 2000, 3(4): 257- 261.

    [34] Schluter D.A variance test for detecting species associations, with some example applications. Ecology, 1984, 65(3): 998- 1005.

    [35] Gotelli N J, Ulrich W. The empirical Bayes approach as a tool to identify non-random species associations. Oecologia, 2010, 162(2): 463- 477.

    [36] Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. The Annals of Statistics, 2001, 29(4): 1165- 1188.

    [37] Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society, 1995, 57(1): 289- 300.

    [38] Gotelli N J, Entsminger G L. Ecosim: Null Models Software for Ecology, Version 7. Acquired Intelligence Inc. and Kesey-Bear, Jericho, VT, USA, 2009.

    [39] Ulrich W. Pairs-a FORTRAN program for studying pair-wise species associations in ecological matrices, Version 1.0. 2008.

    [40] Chesson P. Mechanisms of maintenance of species diversity. Annual Review of Ecology and Systematics, 2000, 31: 343- 366.

    [41] Pianka E R. The structure of lizard communities. Annual Review of Ecology and Systematics, 1973, 4: 53- 74.

    [42] Krasnov B R, Shenbrot G I, Khokhlova I S. Aggregative structure is the rule in communities of fleas: null model analysis. Ecography, 2011, 34(5): 751- 761.

    [43] Williams C B. The generic relations of species in small ecological communities. Journal of Animal Ecology, 1947, 16(1): 11- 18.

    [44] Gotelli N J, Graves G R. Null Models in Ecology. Washington: Smithsonian Institution Press,1996.

    [45] Connor E F, Simberloff D. The assembly of species communities: chance or competition?. Ecology, 1979, 60(6): 1132- 1140.

    [46] Sanderson J G, Moulton M P, Selfridge R G. Null matrices and the analysis of species co-occurrences. Oecologia, 1998, 116(1/2): 275- 283.

    [47] Sfenthourakis S, Tzanatos E, Giokas S. Species co-occurrence: the case of congeneric species and a causal approach to patterns of species association. Global Ecology and Biogeography, 2006, 15(1): 39- 49.

    [48] Niemel? J. Interspecific competition in ground-beetle assemblages (Carabidae): what have we learned?. Oikos, 1993, 66(2): 325- 335.

    [49] Loreau M. Species abundance patterns and the structure of ground-beetle communities. AnnalesZoologiciFennici, 1991, 28(2): 49- 56.

    [50] Kaltsas D, EleniP, Tricha A, Chatzaki M, Mylonas M. Co-occurrence patterns and shaping processes of assembly in beetle and spider communities in eastern Mediterranean maquis formations. Journal of Biological Research, 2012, 18: 187- 197.

    [51] Ulrich W, Zalewski M. Abundance and co-occurrence patterns of core and satellite species of ground beetles on small lake islands. Oikos, 2006, 114(2): 338- 348.

    [52] Ulrich W. Species co-occurrences and neutral models: reassessing J.M. diamond′s assembly rules. Oikos, 2004, 107(3): 603- 609.

    Co-occurrence pattern dynamics of a Carabidae community in aPinuskoraiensisplanted forest

    GAO Meixiang1,2, LI Jingke1,2, SHI Hao1,2, ZHANG Xueping1,2,*, ZHU Jiyuan1,2

    1CollegeofGeographicalSciences,HarbinNormalUniversity,Harbin150025,China2KeyLaboratoryofRemoteSensingMonitoringofGeographicEnvironment,CollegeofHeilongjiangProvince,HarbinNormalUniversity,Harbin150025,China

    One of the most important tasks in community ecology is identifying the co-occurrence patterns of the community and identifying the underlying processes that control such patterns. Null model analysis,which has been facilitated by the processes of biodiversity maintenance, has been used to identify the co-occurrence patterns in communities. However, few studies have focused on co-occurrence pattern dynamics of belowground soil communities. In this study, we identify the co-occurrence pattern dynamics of a Carabidae community in aPinuskoraiensisplanted forest on a small scale, based on five investigations from June 2013 to October 2014. A plot (20 m × 20 m) was established and further divided into 100 squares of 2 m × 2 m at the Maoer Mountain Ecosystem Research Station of the Northeast Forestry University in Heilongjiang Province, Northeast China. One hundred pitfall traps were set to collect the Carabidae community for each investigation. To determine non-random patterns of species co-existence and significant species pairs in Carabidae communities, the null model was used to generate expected patterns in the absence of species interactions. In total, 20 species and 2278 individuals were captured.Carabusbillergimaoershanensiswas the most abundant and most widely distributed species among all communities. Carabidae community composition exhibited significant temporal variability among individuals during the experimental period. With the model algorithm of fixed row and column sums, results of theC-score andV-ratio showed that co-occurrence of species in June 2013 was aggregated, a finding that was inconsistent with Diamond′s assembly rules. Moreover, with the model algorithm of fixed row sums and probability column sums and considering both algorithms, results of theC-score andV-ratio indicated that co-occurrence patterns of Carabidae communities in August and October 2014 were segregated, which was consistent with Diamond′s assembly rules. Although the results of theC-score andV-ratio suggest that the co-occurrence patterns of the Carabidae communities in July and September 2013 were aggregated or segregated based on analysis of the null model, it can be inferred that non-random co-occurrence patterns were detected in all Carabidae communities. Furthermore,based on our observations of the co-occurrence patterns in the Carabidae communities, the Diamond′s assembly rules were partly rejected. Few significant species pairs were detected among all Carabidae communities. A greater number of aggregated species pairs were observed in comparison to segregated species pairs among all communities. The species comprising those significant species pairs were dominant or common species, and were comparatively more widely distributed throughout the plot in each season. In addition, the Pianka (Oik) indices between those dominant and common species were relatively greater than that observed between other species pairs. The results suggest that those significant species pairs might also make important contributions to community composition. The present study demonstrated that a non-random co-occurrence pattern might be common in a small-scale Carabidae community in aPinuskoraiensisplanted forest. This non-random co-occurrence pattern also showed short-term stability during the study period;however, the properties of the sepatterns vary among the different seasons. The significant species pairs possibly contribute directly to observed community co-occurrence patterns, and more focus should be placed on the contributions of those dominant and common Carabidae species in the future.

    co-occurrence pattern; significant species pairs; null model; Carabidae community;Pinuskoraiensis; Maoer Mountain

    國家自然科學基金項目(41471037, 41101049, 41371072, 41430857);黑龍江省普通本科高等學校青年創(chuàng)新人才培養(yǎng)計劃(UNPYSCT- 2015054);哈爾濱師范大學青年學術骨干資助計劃項目(KGB201204) ;哈爾濱師范大學碩士研究生創(chuàng)新科研項目(HSDSSCX2015- 10)

    2015- 04- 05;

    日期:2016- 01- 22

    10.5846/stxb201504050680

    *通訊作者Corresponding author.E-mail: hellozxp@163.com

    高梅香,李景科,石昊,張雪萍,朱紀元.人工紅松林步行蟲(Coleoptera: Carabidae)群落物種共存格局動態(tài)分析.生態(tài)學報,2016,36(20):6591- 6601.

    Gao M X, Li J K, Shi H, Zhang X P, Zhu J Y.Co-occurrence pattern dynamics of a Carabidae community in aPinuskoraiensisplanted forest.Acta Ecologica Sinica,2016,36(20):6591- 6601.

    猜你喜歡
    紅松林隨機性格局
    沒有紅松的紅松林
    紅松林采伐更新的確定方式
    格局
    聯(lián)手共建 努力打造大調(diào)解工作格局
    東折棱河楓樺紅松林與柞樹紅松林結構特征對比分析
    防護林科技(2017年7期)2017-09-01 01:08:43
    淺析電網(wǎng)規(guī)劃中的模糊可靠性評估方法
    考慮負荷與分布式電源隨機性的配電網(wǎng)無功優(yōu)化
    適用于隨機性電源即插即用的模塊化儲能電池柜設計
    小人物的大格局
    小說月刊(2015年12期)2015-04-23 08:51:10
    高見的詩
    滿族文學(2014年5期)2014-09-22 23:50:36
    久久久久久人人人人人| 色精品久久人妻99蜜桃| 狠狠精品人妻久久久久久综合| 久久久久精品久久久久真实原创| 99re6热这里在线精品视频| 亚洲熟女毛片儿| 亚洲情色 制服丝袜| 欧美人与性动交α欧美精品济南到| 日本黄色日本黄色录像| 老熟女久久久| 欧美乱码精品一区二区三区| 大码成人一级视频| 久久免费观看电影| 欧美黑人欧美精品刺激| 亚洲成人一二三区av| 久久 成人 亚洲| 十分钟在线观看高清视频www| 天堂中文最新版在线下载| 久久精品aⅴ一区二区三区四区| 免费黄色在线免费观看| 欧美久久黑人一区二区| 美女主播在线视频| 精品国产一区二区三区久久久樱花| 一级爰片在线观看| 少妇人妻精品综合一区二区| 综合色丁香网| 国产欧美日韩综合在线一区二区| 99国产综合亚洲精品| 精品视频人人做人人爽| 亚洲精品在线美女| 嫩草影院入口| 久久久国产一区二区| 汤姆久久久久久久影院中文字幕| 老汉色av国产亚洲站长工具| 日韩 欧美 亚洲 中文字幕| 精品视频人人做人人爽| 五月开心婷婷网| 成人免费观看视频高清| 午夜福利视频在线观看免费| 欧美日韩精品网址| 久久久久久久久免费视频了| 性少妇av在线| 母亲3免费完整高清在线观看| 亚洲男人天堂网一区| 伊人久久国产一区二区| 97人妻天天添夜夜摸| 日本av手机在线免费观看| 国产成人91sexporn| 伊人亚洲综合成人网| 欧美日本中文国产一区发布| 国产亚洲av片在线观看秒播厂| 国产老妇伦熟女老妇高清| 黄片播放在线免费| 最新在线观看一区二区三区 | 色播在线永久视频| 亚洲国产中文字幕在线视频| 成人黄色视频免费在线看| 免费观看人在逋| 亚洲国产精品999| 乱人伦中国视频| 少妇人妻精品综合一区二区| 男男h啪啪无遮挡| 欧美中文综合在线视频| 最黄视频免费看| 国产福利在线免费观看视频| 91精品三级在线观看| 一个人免费看片子| 亚洲精品在线美女| 欧美国产精品va在线观看不卡| 午夜免费鲁丝| 777米奇影视久久| 99香蕉大伊视频| 男女下面插进去视频免费观看| 亚洲在久久综合| 精品人妻一区二区三区麻豆| 亚洲五月色婷婷综合| 最近中文字幕高清免费大全6| 中文精品一卡2卡3卡4更新| 最新的欧美精品一区二区| 叶爱在线成人免费视频播放| 日韩精品免费视频一区二区三区| 欧美日韩亚洲高清精品| 日韩免费高清中文字幕av| 精品一区二区三区av网在线观看 | 在线看a的网站| www.av在线官网国产| 午夜免费观看性视频| 美女视频免费永久观看网站| 中文乱码字字幕精品一区二区三区| 看免费av毛片| www.精华液| 亚洲伊人久久精品综合| 亚洲人成电影观看| 狠狠婷婷综合久久久久久88av| 人人妻人人爽人人添夜夜欢视频| 亚洲成人免费av在线播放| 飞空精品影院首页| 欧美日韩av久久| 日韩欧美精品免费久久| 久久久久久久久免费视频了| 国产男女内射视频| 少妇人妻久久综合中文| 国产 一区精品| 老鸭窝网址在线观看| 国产精品av久久久久免费| 国产在线视频一区二区| 成人午夜精彩视频在线观看| 久久人人97超碰香蕉20202| 国产精品国产三级国产专区5o| 不卡av一区二区三区| 久久韩国三级中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产精品亚洲av一区麻豆 | 亚洲三区欧美一区| 亚洲一区中文字幕在线| av国产久精品久网站免费入址| 91老司机精品| svipshipincom国产片| 亚洲精品乱久久久久久| 亚洲欧洲国产日韩| 国产精品偷伦视频观看了| 亚洲伊人久久精品综合| 日韩av在线免费看完整版不卡| 99精国产麻豆久久婷婷| 亚洲久久久国产精品| 亚洲精品久久成人aⅴ小说| 精品亚洲成a人片在线观看| 美女视频免费永久观看网站| 亚洲五月色婷婷综合| 成年人午夜在线观看视频| av国产久精品久网站免费入址| 久久精品久久精品一区二区三区| 久久影院123| av网站在线播放免费| 最近最新中文字幕免费大全7| bbb黄色大片| 波多野结衣一区麻豆| 国产极品天堂在线| 欧美精品人与动牲交sv欧美| 国产伦人伦偷精品视频| 丝袜美腿诱惑在线| 国产精品.久久久| 亚洲国产日韩一区二区| 人人妻人人爽人人添夜夜欢视频| 国产午夜精品一二区理论片| 久久狼人影院| 国产野战对白在线观看| 久久久久久人妻| 如何舔出高潮| 热re99久久精品国产66热6| 少妇 在线观看| 9热在线视频观看99| 国产老妇伦熟女老妇高清| 我要看黄色一级片免费的| 人妻人人澡人人爽人人| 99九九在线精品视频| 好男人视频免费观看在线| 国产成人a∨麻豆精品| 亚洲,一卡二卡三卡| 日本av手机在线免费观看| 欧美日韩视频高清一区二区三区二| 精品一区在线观看国产| 欧美中文综合在线视频| 久久久国产一区二区| 一级片'在线观看视频| 成人三级做爰电影| 色综合欧美亚洲国产小说| 99热网站在线观看| av片东京热男人的天堂| 性少妇av在线| 国产精品亚洲av一区麻豆 | 免费黄网站久久成人精品| 国产欧美日韩一区二区三区在线| 国产精品久久久久久久久免| 999精品在线视频| 飞空精品影院首页| 亚洲国产成人一精品久久久| 女人被躁到高潮嗷嗷叫费观| 91精品伊人久久大香线蕉| 午夜福利视频精品| 精品少妇内射三级| 欧美日韩国产mv在线观看视频| 精品国产国语对白av| av国产久精品久网站免费入址| 久久精品国产综合久久久| 在线观看免费高清a一片| 少妇猛男粗大的猛烈进出视频| 久久韩国三级中文字幕| 精品卡一卡二卡四卡免费| 天天躁日日躁夜夜躁夜夜| 中文字幕精品免费在线观看视频| 99久国产av精品国产电影| 成年人午夜在线观看视频| 777久久人妻少妇嫩草av网站| 免费黄色在线免费观看| 久久久精品区二区三区| 老司机靠b影院| 成年av动漫网址| 日韩制服丝袜自拍偷拍| 一个人免费看片子| 久久婷婷青草| 亚洲精华国产精华液的使用体验| 欧美激情 高清一区二区三区| 亚洲,欧美,日韩| 国产成人精品福利久久| 亚洲av电影在线进入| 欧美日韩亚洲综合一区二区三区_| 精品国产露脸久久av麻豆| 久久精品国产亚洲av高清一级| 免费观看av网站的网址| 亚洲av欧美aⅴ国产| 欧美亚洲 丝袜 人妻 在线| 国产精品一国产av| 国产麻豆69| 久久99精品国语久久久| 91aial.com中文字幕在线观看| 久久精品久久久久久噜噜老黄| 亚洲国产av影院在线观看| 爱豆传媒免费全集在线观看| 大片电影免费在线观看免费| 国产精品麻豆人妻色哟哟久久| 欧美日韩综合久久久久久| 中文天堂在线官网| 性少妇av在线| 99久久精品国产亚洲精品| 日韩中文字幕视频在线看片| 欧美少妇被猛烈插入视频| 黑人欧美特级aaaaaa片| 我的亚洲天堂| 午夜91福利影院| 日韩 亚洲 欧美在线| 久久久久国产精品人妻一区二区| 妹子高潮喷水视频| 一级黄片播放器| 777久久人妻少妇嫩草av网站| 好男人视频免费观看在线| 亚洲婷婷狠狠爱综合网| 国产片特级美女逼逼视频| 国产亚洲一区二区精品| 亚洲专区中文字幕在线 | 午夜老司机福利片| 日本猛色少妇xxxxx猛交久久| 国产成人av激情在线播放| 美国免费a级毛片| 亚洲欧美激情在线| 伊人久久国产一区二区| 久久毛片免费看一区二区三区| 性高湖久久久久久久久免费观看| 最新的欧美精品一区二区| 色视频在线一区二区三区| 国产亚洲精品第一综合不卡| a级毛片在线看网站| 欧美人与性动交α欧美精品济南到| 久久国产亚洲av麻豆专区| 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 免费av中文字幕在线| 亚洲精品在线美女| 99精国产麻豆久久婷婷| 夫妻午夜视频| 国产爽快片一区二区三区| 免费观看人在逋| 一本色道久久久久久精品综合| 色婷婷久久久亚洲欧美| 亚洲成人免费av在线播放| 最近最新中文字幕大全免费视频 | 黑人猛操日本美女一级片| 免费黄色在线免费观看| 久久久久精品人妻al黑| 欧美精品一区二区免费开放| 亚洲av电影在线观看一区二区三区| 国产精品国产三级国产专区5o| 人人妻人人爽人人添夜夜欢视频| 国产 一区精品| 久久国产亚洲av麻豆专区| 亚洲精品国产av蜜桃| 在线精品无人区一区二区三| 这个男人来自地球电影免费观看 | 午夜福利一区二区在线看| 哪个播放器可以免费观看大片| 久久精品国产综合久久久| 美女午夜性视频免费| 中文天堂在线官网| 亚洲精品国产色婷婷电影| 丰满饥渴人妻一区二区三| 宅男免费午夜| av在线观看视频网站免费| 国产成人免费观看mmmm| 亚洲精品国产av成人精品| 亚洲成人一二三区av| 成人亚洲欧美一区二区av| 成年人免费黄色播放视频| 性少妇av在线| 卡戴珊不雅视频在线播放| 涩涩av久久男人的天堂| 中文字幕另类日韩欧美亚洲嫩草| e午夜精品久久久久久久| 国产老妇伦熟女老妇高清| 久久久欧美国产精品| 成人黄色视频免费在线看| 男男h啪啪无遮挡| 男人舔女人的私密视频| 中文字幕人妻丝袜制服| 亚洲成人手机| 久久久久久久久久久免费av| 久久国产亚洲av麻豆专区| 午夜福利,免费看| 欧美国产精品一级二级三级| 秋霞伦理黄片| 2018国产大陆天天弄谢| 青春草视频在线免费观看| 成人国产av品久久久| 一区二区三区乱码不卡18| 精品亚洲乱码少妇综合久久| 中文字幕另类日韩欧美亚洲嫩草| 成人三级做爰电影| 一区二区三区激情视频| 岛国毛片在线播放| 国产男女内射视频| 自线自在国产av| 免费高清在线观看视频在线观看| 丁香六月欧美| av在线app专区| 久久天躁狠狠躁夜夜2o2o | 午夜日韩欧美国产| 欧美在线黄色| 日韩人妻精品一区2区三区| 久久久久久久久免费视频了| 亚洲国产看品久久| 成人手机av| av卡一久久| 国产精品香港三级国产av潘金莲 | 黄色一级大片看看| 男女边吃奶边做爰视频| 999精品在线视频| 香蕉丝袜av| 久久久久人妻精品一区果冻| 欧美日韩亚洲高清精品| 亚洲av电影在线进入| 又粗又硬又长又爽又黄的视频| 亚洲欧美日韩另类电影网站| a级片在线免费高清观看视频| 久久性视频一级片| 精品国产国语对白av| 亚洲第一av免费看| 亚洲精品aⅴ在线观看| 欧美中文综合在线视频| 欧美成人午夜精品| 欧美日韩国产mv在线观看视频| avwww免费| 国产一区二区三区av在线| 一级毛片我不卡| 欧美黑人精品巨大| 久久精品久久久久久噜噜老黄| 热99国产精品久久久久久7| 老汉色av国产亚洲站长工具| 一边亲一边摸免费视频| 亚洲精品中文字幕在线视频| 免费女性裸体啪啪无遮挡网站| 久久久久久久精品精品| 丝袜喷水一区| 亚洲精品自拍成人| av电影中文网址| 亚洲精品美女久久av网站| 久久久久精品国产欧美久久久 | 国产欧美日韩综合在线一区二区| 热99国产精品久久久久久7| 久久毛片免费看一区二区三区| 亚洲伊人久久精品综合| 欧美亚洲日本最大视频资源| 一级,二级,三级黄色视频| 免费日韩欧美在线观看| 国产一区二区激情短视频 | 国产男女超爽视频在线观看| 丰满乱子伦码专区| 狂野欧美激情性bbbbbb| 久热这里只有精品99| 中文字幕人妻熟女乱码| 男人爽女人下面视频在线观看| 男女床上黄色一级片免费看| 午夜老司机福利片| 建设人人有责人人尽责人人享有的| 你懂的网址亚洲精品在线观看| 色视频在线一区二区三区| 秋霞伦理黄片| 桃花免费在线播放| 老汉色∧v一级毛片| 伊人亚洲综合成人网| 一区二区三区乱码不卡18| 亚洲,欧美精品.| xxx大片免费视频| 久久久久久久国产电影| 欧美日韩成人在线一区二区| xxxhd国产人妻xxx| 亚洲精品在线美女| 亚洲激情五月婷婷啪啪| 99久国产av精品国产电影| 欧美日韩综合久久久久久| 久久综合国产亚洲精品| 久久毛片免费看一区二区三区| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 别揉我奶头~嗯~啊~动态视频 | 日本黄色日本黄色录像| 久久久久视频综合| 精品国产超薄肉色丝袜足j| av有码第一页| 大片电影免费在线观看免费| 欧美成人精品欧美一级黄| 国产福利在线免费观看视频| 人人妻人人添人人爽欧美一区卜| 亚洲国产看品久久| 午夜老司机福利片| 啦啦啦啦在线视频资源| 欧美在线一区亚洲| 一区二区三区乱码不卡18| videosex国产| 国产精品二区激情视频| 观看av在线不卡| 国产精品久久久久成人av| 看十八女毛片水多多多| 久久精品久久精品一区二区三区| 久热这里只有精品99| 成人国产av品久久久| 夫妻午夜视频| 女的被弄到高潮叫床怎么办| 中文字幕另类日韩欧美亚洲嫩草| 哪个播放器可以免费观看大片| 热re99久久精品国产66热6| av在线观看视频网站免费| 老司机在亚洲福利影院| 国产日韩欧美亚洲二区| 日本欧美国产在线视频| 国产有黄有色有爽视频| av国产精品久久久久影院| 免费看不卡的av| 成人国语在线视频| 免费黄网站久久成人精品| 街头女战士在线观看网站| 国产精品熟女久久久久浪| 午夜福利,免费看| 男女边吃奶边做爰视频| 免费黄频网站在线观看国产| 国产精品秋霞免费鲁丝片| 中文字幕制服av| 观看av在线不卡| 一二三四在线观看免费中文在| 亚洲国产精品一区二区三区在线| 国语对白做爰xxxⅹ性视频网站| 伦理电影大哥的女人| 高清不卡的av网站| 宅男免费午夜| 多毛熟女@视频| 国产日韩欧美在线精品| 日本wwww免费看| 中文字幕制服av| 国产国语露脸激情在线看| 国产精品蜜桃在线观看| 欧美亚洲 丝袜 人妻 在线| 久久韩国三级中文字幕| 久久亚洲国产成人精品v| 母亲3免费完整高清在线观看| 国产高清国产精品国产三级| 亚洲第一av免费看| 精品一区二区免费观看| 韩国精品一区二区三区| 丝袜美腿诱惑在线| 中文字幕高清在线视频| 久久精品国产亚洲av涩爱| 久久ye,这里只有精品| 天天影视国产精品| 国产一区亚洲一区在线观看| 亚洲av成人不卡在线观看播放网 | 在线精品无人区一区二区三| 欧美另类一区| 亚洲成人国产一区在线观看 | 免费日韩欧美在线观看| 亚洲一码二码三码区别大吗| 久久久久精品性色| 1024视频免费在线观看| 免费在线观看视频国产中文字幕亚洲 | 欧美日韩一区二区视频在线观看视频在线| 99久国产av精品国产电影| 男的添女的下面高潮视频| 九色亚洲精品在线播放| 国产极品天堂在线| 亚洲精品在线美女| 欧美日韩视频精品一区| 性色av一级| 99香蕉大伊视频| 高清在线视频一区二区三区| 精品久久蜜臀av无| 久久精品aⅴ一区二区三区四区| 久久ye,这里只有精品| 国产精品二区激情视频| xxxhd国产人妻xxx| 五月开心婷婷网| 国语对白做爰xxxⅹ性视频网站| 性高湖久久久久久久久免费观看| 伦理电影免费视频| 久久精品aⅴ一区二区三区四区| 午夜久久久在线观看| 午夜免费男女啪啪视频观看| 大片电影免费在线观看免费| 亚洲国产看品久久| 国产精品亚洲av一区麻豆 | 在线观看免费日韩欧美大片| 丝袜脚勾引网站| 男女之事视频高清在线观看 | 可以免费在线观看a视频的电影网站 | 国产精品国产av在线观看| 在线观看三级黄色| 涩涩av久久男人的天堂| 亚洲欧美成人综合另类久久久| 欧美 日韩 精品 国产| 看非洲黑人一级黄片| 中文字幕亚洲精品专区| 别揉我奶头~嗯~啊~动态视频 | 久久精品久久久久久噜噜老黄| 亚洲精品在线美女| 又大又爽又粗| 777久久人妻少妇嫩草av网站| 欧美激情高清一区二区三区 | 国产欧美亚洲国产| 两个人看的免费小视频| 久热这里只有精品99| 亚洲精品国产一区二区精华液| 天天躁狠狠躁夜夜躁狠狠躁| 操美女的视频在线观看| 久久国产精品男人的天堂亚洲| 国产男人的电影天堂91| 亚洲欧美色中文字幕在线| 国产精品亚洲av一区麻豆 | 国产精品香港三级国产av潘金莲 | 亚洲免费av在线视频| 精品亚洲成国产av| 免费人妻精品一区二区三区视频| 亚洲精品成人av观看孕妇| 街头女战士在线观看网站| 久久久久久久国产电影| 日韩制服骚丝袜av| 少妇精品久久久久久久| 桃花免费在线播放| 国产一区有黄有色的免费视频| 色综合欧美亚洲国产小说| 日韩 欧美 亚洲 中文字幕| 成人三级做爰电影| 久久亚洲精品不卡| 亚洲天堂国产精品一区在线| 欧美黑人欧美精品刺激| 啦啦啦韩国在线观看视频| 国产三级黄色录像| 亚洲av成人一区二区三| 黄频高清免费视频| 97人妻天天添夜夜摸| 亚洲五月天丁香| 日韩中文字幕欧美一区二区| 99在线视频只有这里精品首页| 日韩欧美国产一区二区入口| 十分钟在线观看高清视频www| 巨乳人妻的诱惑在线观看| 国产av一区在线观看免费| 亚洲久久久国产精品| 波多野结衣巨乳人妻| 欧美激情高清一区二区三区| 曰老女人黄片| 亚洲成av人片免费观看| 欧美+亚洲+日韩+国产| 十分钟在线观看高清视频www| 精品国产一区二区三区四区第35| 国产精品永久免费网站| 久久精品国产清高在天天线| 国产精品电影一区二区三区| 极品教师在线免费播放| 亚洲av五月六月丁香网| 亚洲成av人片免费观看| 久久久久亚洲av毛片大全| 少妇的丰满在线观看| 在线永久观看黄色视频| 精品第一国产精品| 国产精品爽爽va在线观看网站 | 精品电影一区二区在线| 伦理电影免费视频| 免费女性裸体啪啪无遮挡网站| 老汉色∧v一级毛片| 法律面前人人平等表现在哪些方面| 精品久久久久久久毛片微露脸| 午夜两性在线视频| 老司机午夜福利在线观看视频| 久久久水蜜桃国产精品网| 97人妻天天添夜夜摸| 极品人妻少妇av视频| 免费在线观看完整版高清| 亚洲va日本ⅴa欧美va伊人久久| 精品熟女少妇八av免费久了| 午夜福利一区二区在线看| 咕卡用的链子| 如日韩欧美国产精品一区二区三区| 国产又色又爽无遮挡免费看| 成人精品一区二区免费| 亚洲人成电影观看| 美女免费视频网站| 麻豆一二三区av精品| 国产成人影院久久av| 搡老岳熟女国产| 真人做人爱边吃奶动态| 十八禁人妻一区二区| 国产成人系列免费观看| 日日摸夜夜添夜夜添小说| www.www免费av| 搡老岳熟女国产| 成人精品一区二区免费| 一区二区三区精品91| 巨乳人妻的诱惑在线观看| 一区二区三区国产精品乱码|