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

    祁連山東段不同放牧強(qiáng)度下高原鼢鼠棲息地選擇分析

    2016-02-23 02:48:04楚彬花立民周延山劉麗李帥
    草業(yè)學(xué)報(bào) 2016年1期

    楚彬,花立民,周延山,劉麗,李帥

    (甘肅農(nóng)業(yè)大學(xué)草業(yè)學(xué)院,甘肅 蘭州 730070)

    ?

    祁連山東段不同放牧強(qiáng)度下高原鼢鼠棲息地選擇分析

    楚彬,花立民*,周延山,劉麗,李帥

    (甘肅農(nóng)業(yè)大學(xué)草業(yè)學(xué)院,甘肅 蘭州 730070)

    摘要:本研究以祁連山東段高寒草甸棲息的高原鼢鼠為研究對(duì)象,探討不同放牧強(qiáng)度下高原鼢鼠棲息地特征及選擇因素,為合理控制草原鼠害及保護(hù)草地生物多樣性提供科學(xué)依據(jù)。本研究選擇4個(gè)不同放牧強(qiáng)度小區(qū),調(diào)查每個(gè)小區(qū)的高原鼢鼠相對(duì)種群密度、土壤緊實(shí)度和容重、植物學(xué)組成、植物物種豐富度、地上生物量、可食牧草產(chǎn)量及各植物功能群產(chǎn)量,并分析高原鼢鼠相對(duì)種群密度與這些要素之間的關(guān)系。結(jié)果表明,放牧強(qiáng)度與高原鼢鼠相對(duì)種群密度有正相關(guān)關(guān)系但不顯著(P>0.05),重度放牧區(qū)下高原鼢鼠相對(duì)種群密度最高,禁牧區(qū)高原鼢鼠相對(duì)種群密度最低;高原鼢鼠相對(duì)種群密度與土壤緊實(shí)度呈顯著負(fù)相關(guān)(P<0.05),與土壤容重和地上生物量呈負(fù)相關(guān)但不顯著(P>0.05);高原鼢鼠相對(duì)種群密度與雜類草產(chǎn)量有正相關(guān)關(guān)系但不顯著(P>0.05),與可食牧草產(chǎn)量、禾草產(chǎn)量和物種豐富度有負(fù)相關(guān)關(guān)系但不顯著(P>0.05);多元線性回歸方程表明影響高原鼢鼠棲息地選擇的主要因素是土壤緊實(shí)度(Y=205.36-0.145X1,X1為土壤緊實(shí)度)。

    關(guān)鍵詞:高原鼢鼠;放牧強(qiáng)度;棲息地特征;棲息地選擇

    青藏高原高寒草地是世界上海拔最高、面積最大的草地生態(tài)系統(tǒng),在畜牧生產(chǎn)、涵養(yǎng)水源和維持生物多樣性等方面發(fā)揮著重要作用[1]。近年來(lái),由于全球氣候變化及人為因素的影響,高寒草地退化嚴(yán)重,生產(chǎn)力下降以及生態(tài)環(huán)境惡化已經(jīng)嚴(yán)重威脅著區(qū)域乃至全國(guó)社會(huì)和經(jīng)濟(jì)的可持續(xù)發(fā)展[2-3]。在影響草地生態(tài)系統(tǒng)的生物因子中,家畜的過(guò)度放牧以及嚙齒動(dòng)物種群數(shù)量的增加都具有重要的干擾作用[4]。放牧不僅作用于植物群落和土壤[5-7],而且也影響到草原嚙齒動(dòng)物棲息地選擇[8]。嚙齒動(dòng)物棲息地適合度增加可以導(dǎo)致種群數(shù)量增加和棲息地面積擴(kuò)大,加劇草地退化[9]。但是,放牧強(qiáng)度與地下嚙齒動(dòng)物棲息地適合度的關(guān)系,目前尚不清晰。因此,研究放牧干擾對(duì)草原嚙齒動(dòng)物棲息地選擇影響,對(duì)明晰其在草地生態(tài)系統(tǒng)中的地位和作用,以及科學(xué)防治草原鼠害具有重要意義。

    高原鼢鼠(Eospalaxbaileyi)是倉(cāng)鼠科中4種營(yíng)地下生活的嚙齒動(dòng)物之一,廣泛分布于青藏高原,是青藏高原生物多樣性的重要組成部分,在草地生態(tài)系統(tǒng)食物網(wǎng)及其相對(duì)應(yīng)的能量流通和物質(zhì)循環(huán)中有著獨(dú)特的地位[10-12]。但高原鼢鼠采食牧草以及挖掘過(guò)程中向地表推出土丘,對(duì)草地生產(chǎn)和生態(tài)也會(huì)產(chǎn)生負(fù)面影響。同其他動(dòng)物一樣,高原鼢鼠棲息地的選擇受諸多因素影響,以往的研究主要集中在對(duì)高原鼢鼠棲息地植被特征的分析,較少考慮到放牧強(qiáng)度的影響[13-14]。而本研究重點(diǎn)考慮不同放牧強(qiáng)度下,通過(guò)調(diào)查不同放牧強(qiáng)度區(qū)內(nèi)高原鼢鼠相對(duì)種群密度以及土壤和植被因素,重點(diǎn)分析放牧干擾下高原鼢鼠棲息地的選擇。

    1材料與方法

    1.1研究樣地概況

    樣地設(shè)于甘肅省武威市天祝藏族自治縣抓喜秀龍鄉(xiāng)馬營(yíng)灘地區(qū),該區(qū)位于東祁連山的天祝金強(qiáng)河河谷,地理坐標(biāo)N 37°12′,E 102°46′,海拔2937 m,氣候寒冷潮濕,太陽(yáng)輻射強(qiáng)。年均溫-0.1℃,1月平均溫度-18.3℃,7月平均溫度12.7℃,>0℃年積溫1380℃;年降水量416 mm,多為地形雨,集中于7,8,9三個(gè)月。植被類型為高寒草甸,以莎草科、禾本科植物為主,混有其他雜類草,草群蓋度達(dá)到90%左右,土壤為高寒草甸土,主要植物有珠芽蓼(Polygonumviviparum)、鵝絨委陵菜(Potentillaanserina)、早熟禾(Poaannua)、矮嵩草(Kobresiahumilis)、線葉嵩草(Kobresiacapillifoli)、垂穗披堿草(Elymusdahuricus)、球花蒿(Artemisiasmithii)、黃芪(Astragalusmembranaceus)等[15]。

    1.2樣地設(shè)置及實(shí)驗(yàn)時(shí)間

    樣地選取冬季牧場(chǎng),依據(jù)不同放牧強(qiáng)度選擇4塊樣地,每塊樣地面積3 hm2。放牧?xí)r間為當(dāng)年11月至次年5月,放牧期均在3年以上,而禁牧區(qū)為已連續(xù)3年禁牧的草地。放牧強(qiáng)度以家畜頭數(shù)、草地面積和利用時(shí)間而定,依照當(dāng)?shù)仄骄拍翉?qiáng)度,分別劃分為重度放牧40羊單位/hm2(high grazing, HG)、中度放牧27羊單位/hm2(moderate grazing, MG)、輕度放牧10羊單位/hm2(light grazing, LG)、將禁牧區(qū)作為對(duì)照(control check, CK)。樣地皆做相互獨(dú)立的圍欄保護(hù),該圍欄的設(shè)置是防止家畜在另一樣地進(jìn)行采食。所有植被和土壤調(diào)查時(shí)間均在2014年8月。

    1.3高原鼢鼠相對(duì)種群密度調(diào)查

    依據(jù)何俊齡等[16]單位面積新鼠丘數(shù)(當(dāng)年8月以前所形成)表示各樣地高原鼢鼠的相對(duì)種群密度。即相對(duì)種群密度=當(dāng)年8月前新生鼠丘數(shù)量/樣地面積。新土丘數(shù)量調(diào)查采用RTK(real-time kinematic)技術(shù)。RTK是目前使用測(cè)繪、建筑等方面的快速高精度定位技術(shù)[17]。首先設(shè)置RTK基準(zhǔn)站,并將手持移動(dòng)終端與基準(zhǔn)站建立通訊聯(lián)系。然后將手持移動(dòng)終端放置在每個(gè)新土丘上,并記錄其準(zhǔn)確的經(jīng)緯度位置,最后將手持移動(dòng)終端記錄的新鼠丘位置輸入Arcgis軟件,統(tǒng)計(jì)不同樣地的新鼠丘數(shù)作為高原鼢鼠相對(duì)種群密度。根據(jù)高原鼢鼠相對(duì)種群密度,依次劃分為高密度區(qū)(Ⅰ)、次高密度區(qū)(Ⅱ)、中密度區(qū)(Ⅲ)以及低密度區(qū)(Ⅳ),與之相對(duì)應(yīng)的放牧強(qiáng)度區(qū)為HG、LG、MG和CK(圖1)。

    1.4植被性狀調(diào)查

    1.4.1地上生物量及植物學(xué)組成每個(gè)放牧區(qū)設(shè)置3個(gè)重復(fù),每個(gè)重復(fù)面積為1 hm2。利用DWR(dry-weight rank)法測(cè)量地上產(chǎn)草量和植物學(xué)組成,統(tǒng)計(jì)植物學(xué)組成按照其可食牧草所占比例,計(jì)算出可食牧草的產(chǎn)量[18-27]。

    1.4.2植物物種豐富度利用DWR草地快速監(jiān)測(cè)法,記錄每一樣地出現(xiàn)的植物種數(shù)。

    1.4.3植物功能群生物量按照DWR法所記錄的植物學(xué)組成以及地上生物量,將其劃分為禾本科、莎草科、雜類草3個(gè)經(jīng)濟(jì)功能群,按照每個(gè)功能群所占比例計(jì)算各功能群生物量。

    1.5土壤性狀調(diào)查

    1.5.1土壤緊實(shí)度在每個(gè)樣地均勻選擇27個(gè)點(diǎn),每個(gè)點(diǎn)用SC-900土壤緊實(shí)度儀每隔2.5 cm測(cè)量土壤緊實(shí)度,深度為0~20 cm,并在此點(diǎn)附近重復(fù)3次,然后求平均值。

    1.5.2土壤容重每個(gè)樣地內(nèi)均勻選取27個(gè)點(diǎn),每個(gè)點(diǎn)附近挖取1個(gè)剖面,深度20 cm,自表土向下每隔10 cm用100 cm3環(huán)刀取一個(gè)土壤樣品并裝入鋁盒,采用鋁盒烘干法測(cè)定其容重[28]。

    1.6數(shù)據(jù)分析

    用Excel 2007繪制圖表;采用SPSS 17.0軟件對(duì)不同樣地間土壤、植被因子進(jìn)行One-way ANOVA分析,分析樣地間因子差異顯著性;對(duì)高原鼢鼠相對(duì)種群密度和土壤、植被因子進(jìn)行簡(jiǎn)單相關(guān)性分析及多元回歸分析,制作相關(guān)系數(shù)矩陣及最優(yōu)方程。

    2結(jié)果與分析

    2.1高原鼢鼠相對(duì)種群密度與放牧強(qiáng)度的關(guān)系

    圖1 高原鼢鼠相對(duì)種群密度與放牧強(qiáng)度的關(guān)系Fig.1 The relationship between grazing intensity andrelative population density of plateau zokor

    從圖1可以看出,放牧強(qiáng)度(羊單位/hm2)依次是HG(40)>MG(27)>LG(10)>CK(0),種群密度(土丘/hm2)的變化為Ⅰ區(qū)(55.67)>Ⅱ區(qū)(30.00)>Ⅲ區(qū)(16.67)>Ⅳ區(qū)(6.67)。HG區(qū)高原鼢鼠相對(duì)種群密度最高,LG區(qū)和MG區(qū)次之,CK區(qū)高原鼢鼠相對(duì)種群密度最低。高原鼢鼠種群相對(duì)密度與放牧強(qiáng)度有正相關(guān)關(guān)系但不顯著(P>0.05,R=0.646)。

    2.2高原鼢鼠相對(duì)種群密度與地上生物量和可食牧草產(chǎn)量的關(guān)系

    由圖2可以看出,高原鼢鼠低密度區(qū)(Ⅳ)地上生物量與可食牧草產(chǎn)量顯著高于其他3個(gè)區(qū)(P<0.05)。高原鼢鼠高密度區(qū)(Ⅰ)、次高密度區(qū)(Ⅱ)、中密度區(qū)(Ⅲ)三者產(chǎn)草量無(wú)顯著性差異(P>0.05)。

    2.3高原鼢鼠相對(duì)種群密度與植物物種豐富度的關(guān)系

    從圖3可以看出,植物物種豐富度隨著高原鼢鼠相對(duì)種群密度的升高,出現(xiàn)先升后降的現(xiàn)象,其峰值出現(xiàn)在高原鼢鼠次高密度區(qū)(Ⅱ),最低值出現(xiàn)在高原鼢鼠低密度區(qū)(Ⅳ)。

    2.4高原鼢鼠相對(duì)種群密度與植物功能群產(chǎn)量的關(guān)系

    由圖4可以看出,禾草功能群和莎草功能群產(chǎn)量隨著高原鼢鼠相對(duì)種群密度的升高表現(xiàn)出降低的趨勢(shì),低密度區(qū)(Ⅳ)的禾草和莎草產(chǎn)量顯著高于其他3個(gè)區(qū)(P<0.05)。而對(duì)于雜類草功能群產(chǎn)量來(lái)說(shuō),隨著種群密度的上升呈現(xiàn)增大趨勢(shì),最大值出現(xiàn)在高原鼢鼠高密度區(qū)(Ⅰ),最小值出現(xiàn)在高原鼢鼠低密度區(qū)(Ⅳ),低密度區(qū)的雜類草產(chǎn)量顯著低于其他3個(gè)區(qū)(P<0.05)。

    圖2 高原鼢鼠相對(duì)種群密度與地上生物量和可食牧草產(chǎn)量的關(guān)系 Fig.2 The relationship between aboveground biomass,edible herbageyield and relative population density of plateau zokor

    圖3 高原鼢鼠相對(duì)種群密度與物種豐富度的關(guān)系Fig.3 The relationship between plants richnessand relative population density of plateau zokor

    不同字母表示處理間差異顯著(P<0.05),下同。Values with different letters indicate significant difference among treatments atP<0.05, the same below.

    2.5高原鼢鼠相對(duì)種群密度與土壤緊實(shí)度的關(guān)系

    圖4 高原鼢鼠相對(duì)種群密度與植物功能群產(chǎn)量的關(guān)系 Fig.4 The relationship between the field of plant functionalgroups and relative population density of plateau zokor

    從圖5可以看出,土壤深度在0~20 cm,高原鼢鼠中密度區(qū)(Ⅲ)的土壤緊實(shí)度最大,高密度區(qū)(Ⅰ)最小,高原鼢鼠高密度區(qū)(Ⅰ)的土壤緊實(shí)度顯著低于次高密度區(qū)、中密度區(qū)和低密度區(qū)(P<0.05),次高密度區(qū)(Ⅱ)、中密度區(qū)(Ⅲ)和低密度區(qū)(Ⅳ)三者無(wú)顯著差異(P>0.05)。說(shuō)明高原鼢鼠在緊實(shí)度較低的土壤分布相對(duì)較多。

    2.6高原鼢鼠相對(duì)種群密度與土壤容重的關(guān)系

    從圖6可以看出,土壤深度在0~20 cm,高原鼢鼠低密度區(qū)(Ⅳ)的土壤容重最大,高密度區(qū)(Ⅰ)最小,且高原鼢鼠低密度區(qū)(Ⅳ)的土壤容重顯著大于其他3個(gè)區(qū)(P<0.05),而高原鼢鼠高密度區(qū)(Ⅰ)、次高密度區(qū)(Ⅱ)、中密度區(qū)(Ⅲ)三者土壤容重?zé)o顯著性差異(P>0.05)。

    2.7土壤、植被性狀與高原鼢鼠相對(duì)種群密度相關(guān)性及回歸分析

    土壤、植被性狀與高原鼢鼠相對(duì)種群密度相關(guān)系數(shù)表明(表1),高原鼢鼠相對(duì)種群密度與土壤緊實(shí)度有顯著負(fù)相關(guān)關(guān)系(P<0.05);雜類草產(chǎn)量與土壤容重有顯著負(fù)相關(guān)關(guān)系(P<0.05),地上生物量與土壤容重有顯著正相關(guān)關(guān)系(P<0.05),禾草產(chǎn)量與雜類草產(chǎn)量有顯著負(fù)相關(guān)關(guān)系(P<0.05)。

    圖5 高原鼢鼠相對(duì)種群密度與土壤緊實(shí)度的關(guān)系Fig.5 The relationship between soil compaction andrelative population density of plateau zokor

    圖6 高原鼢鼠相對(duì)種群密度與土壤容重的關(guān)系Fig.6 The relationship between soil bulk density andrelative population density of plateau zokor

    項(xiàng)目ItemX0X1X2X3X4X5X6X10.577X2-0.986*-0.667X30.995**0.611-0.976*X40.983*0.631-0.959*0.996**X5-0.5380.1840.389-0.560-0.579X60.995**0.637-0.984*-0.999**0.994**-0.523X7-0.778-0.956*0.831-0.810-0.8290.098-0.827

    X0:土壤容重;X1:土壤緊實(shí)度;X2:雜類草產(chǎn)量;X3:禾草產(chǎn)量;X4:地上生物量;X5:植物物種豐富度;X6:可食牧草產(chǎn)量;X7:高原鼢鼠相對(duì)種群密度?!?”和“**”表示兩種因子的相關(guān)性達(dá)到顯著水平(P<0.05)和極顯著水平(P<0.01)。

    X0:Soil bulk density;X1:Soil compaction;X2:Forb biomass;X3:Grass biomass;X4:Aboveground biomass;X5:Plants richness;X6:Edible forage biomass;X7:Population density of plateau zokor.“*” and “**”indicate significant correlation at the level of 0.05 and 0.01 respectively.

    為了更好說(shuō)明高原鼢鼠相對(duì)種群密度與植被、土壤因子的關(guān)系,采用逐步回歸分析法得到最優(yōu)回歸方程Y=205.36-0.145X1,X1為土壤緊實(shí)度。該方程表明土壤緊實(shí)度對(duì)高原鼢鼠相對(duì)種群密度影響最大。

    3討論與結(jié)論

    國(guó)內(nèi)對(duì)高原鼢鼠棲息地研究開展較早且取得多項(xiàng)成果[29]。但是大部分研究沒有考慮放牧干擾的影響,而是重點(diǎn)研究高原鼢鼠種群密度與土壤及植被特性的關(guān)系[14,30]。本試驗(yàn)通過(guò)劃分4個(gè)不同放牧強(qiáng)度,研究不同放牧強(qiáng)度下高原鼢鼠相對(duì)種群密度與土壤物理性質(zhì)、植被性狀的關(guān)系,發(fā)現(xiàn)高原鼢鼠喜歡棲息地土壤疏松、雜類草產(chǎn)量高的草地。

    棲息于不同放牧區(qū)的高原鼢鼠種群密度不盡相同。隨著放牧強(qiáng)度的增加,優(yōu)良牧草減少而雜類草增多,而雜類草又是高原鼢鼠喜愛的食物資源。本試驗(yàn)重牧區(qū)的雜類草產(chǎn)量最高(157.05 g/m2),高原鼢鼠相對(duì)種群密度也高于其他3個(gè)區(qū)(55.67丘/hm2)。而禁牧區(qū)的雜類草產(chǎn)量最低(95.34 g/m2),相應(yīng)的鼠群密度也最低(6.67丘/hm2),因此高原鼢鼠偏好于雜類草產(chǎn)量較高的棲息地。此研究結(jié)果與劉偉等[30]研究不同放牧強(qiáng)度對(duì)植物及嚙齒動(dòng)物作用所得出的結(jié)果一致。在物種豐富度方面,張軍等[14]在植被性狀與高原鼢鼠棲息地適合度的關(guān)系研究中認(rèn)為適宜的鼠密度有利于增加物種多樣性。而本試驗(yàn)關(guān)于物種豐富度與種群密度的關(guān)系中可以看出其符合“中度干擾理論”[31],即一定的種群密度及放牧強(qiáng)度有利于增加物種的豐富度,高原鼢鼠種群密度較低的禁牧區(qū)物種豐富度反而最低,但是黃彬[32]認(rèn)為植物的物種數(shù)隨著鼢鼠種群密度的遞增而出現(xiàn)下降趨勢(shì),二者呈現(xiàn)顯著負(fù)相關(guān)關(guān)系,本試驗(yàn)的結(jié)果與此結(jié)果有所出入,有可能是因?yàn)楸驹囼?yàn)增加了放牧因素,較低的放牧強(qiáng)度可能有利于增加物種數(shù)。在地上生物量方面,重牧區(qū)的地上生物量顯著低于禁牧區(qū)(P<0.05),高原鼢鼠相對(duì)種群密度與地上生物量有負(fù)相關(guān)關(guān)系,但不顯著(R=-0.829,P>0.05),因?yàn)橹啬羺^(qū)家畜的采食程度遠(yuǎn)遠(yuǎn)大于禁牧區(qū),而且高原鼢鼠采食地下根系從而影響地上植物的生長(zhǎng),因此出現(xiàn)這種現(xiàn)象。

    高原鼢鼠營(yíng)地下生活,其獨(dú)特的生活習(xí)性導(dǎo)致每年都將大量的土壤掘出地面形成土丘,土壤物理性質(zhì)決定著高原鼢鼠的種群密度及分布。不同的放牧強(qiáng)度,家畜的踐踏程度不同,導(dǎo)致土壤物理性質(zhì)的改變。本研究結(jié)果表明,重牧區(qū)的土壤緊實(shí)度、容重顯著低于其他3個(gè)區(qū),高原鼢鼠相對(duì)種群密度與土壤緊實(shí)度呈現(xiàn)顯著負(fù)相關(guān)關(guān)系(R=-0.956,P<0.05),與土壤容重有負(fù)相關(guān)關(guān)系但不顯著(R=-0.778,P>0.05)。說(shuō)明在重牧區(qū)雖然家畜的踐踏程度要高于其他3個(gè)區(qū),使表層土壤更加緊實(shí),但是重牧區(qū)的高原鼢鼠相對(duì)種群密度也較大,因此地下洞道數(shù)量多,結(jié)構(gòu)復(fù)雜,而所研究區(qū)域的高原鼢鼠平均洞道深度在0~20 cm內(nèi),所以會(huì)出現(xiàn)重牧區(qū)的土壤緊實(shí)度顯著低于其他3個(gè)區(qū)(P<0.05)。同時(shí),本試驗(yàn)通過(guò)多元線性回歸分析表明,土壤緊實(shí)度是影響高原鼢鼠棲息地選擇的主要因素。Wang等[33]在研究高原鼢鼠的挖掘行為與土壤硬度的關(guān)系中認(rèn)為高原鼢鼠明顯選擇土壤硬度較小的地方棲息,黃彬[32]研究高原鼢鼠種群密度與土壤特性的關(guān)系中也表明隨著鼠群密度的增加,土壤硬度、土壤容重下降,本試驗(yàn)的研究結(jié)果與此結(jié)果相似。土壤緊實(shí)度是一個(gè)反映土壤容重和水分以及質(zhì)地的綜合指標(biāo)[34],本研究表明適宜高原鼢鼠棲息的生境選擇主要傾向于土壤疏松的草地,其他因子如雜類草影響次之。

    棲息地選擇是動(dòng)物生態(tài)學(xué)的一個(gè)傳統(tǒng)研究領(lǐng)域,也是動(dòng)物生態(tài)學(xué)眾多分支的基礎(chǔ),同時(shí)棲息地選擇研究可為動(dòng)物保護(hù)措施的制定、利用生態(tài)手段防治草原鼠害提供科學(xué)依據(jù)[35]。但是目前關(guān)于棲息地選擇理論的研究還相當(dāng)匱乏,基礎(chǔ)還很薄弱[36]。目前有關(guān)棲息地質(zhì)量評(píng)價(jià)的指標(biāo)存在很多爭(zhēng)議,但國(guó)內(nèi)大多數(shù)學(xué)者認(rèn)為動(dòng)物種群密度的高低依然是評(píng)價(jià)棲息地質(zhì)量高低的指標(biāo)[37]。因此,本研究依然將高原鼢鼠種群密度作為重要指標(biāo)來(lái)確定其棲息地質(zhì)量的好壞,從本研究的結(jié)果可以看出,在放牧干擾下,重度放牧區(qū)鼠群密度最高,棲息地特征表現(xiàn)為土壤緊實(shí)度、容重降低,雜類草產(chǎn)量較大,說(shuō)明高原鼢鼠偏好此類型的棲息地。但是,由于自然條件下各種棲息地特征的非獨(dú)立性,難以從中判斷動(dòng)物個(gè)體棲息地偏好實(shí)際所依據(jù)的信號(hào),而且動(dòng)物偏好某一類型的棲息地既可能由遺傳決定,也可能受后期生長(zhǎng)的影響,還可能受所處生活史階段等生理狀況的影響。因此,對(duì)于高原鼢鼠棲息地選擇研究還需要大量的工作要做。同時(shí),目前關(guān)于棲息地選擇研究大多集中于對(duì)“棲息地特征”的研究[38],多數(shù)研究均建立在對(duì)自然條件下的調(diào)查數(shù)據(jù)的相關(guān)性分析,而相關(guān)性分析的缺陷是相關(guān)現(xiàn)象未必存在必然的因果關(guān)系,因此為了更準(zhǔn)確分析動(dòng)物個(gè)體的棲息地偏好,就需要可控試驗(yàn)的介入,通過(guò)對(duì)條件控制,嚴(yán)格設(shè)計(jì)棲息地選擇試驗(yàn),這也是在今后關(guān)于棲息地研究應(yīng)該參考的方向。

    References:

    [1]Zhou X R, Guo Z G, Guo X H. The role of plateau pika and plateau zokor in alpine meadow. Pratacultural Science, 2010, 27(5): 38-44.

    [2]Li W F, Tang Z M. Environment and safety of Qinghai-Tibet railway. Journal of Safety and Environment, 2003, 3(3): 65-67.

    [3]Fan N C, Shi Y Z. A revision of the zokors of subgenusEospalax. Acta Theriologica Sinica, 1982, 2(1): 183-200.

    [4]Melladoa M, Olverab A, Querob A. Dietary overlap between prairie dog (Cynomysmexicanus) and beef cattle in a desert rangeland of northern Mexico. Journal of Arid Environments, 2005, 62(3): 449-458.

    [5]Rosiere R E. An evaluation of grazing intensity influence on california annual range. Range Manage, 1987, 40(2): 160-165.

    [6]Zhang Y W, Han J G, Li Z Q. A study of the effects of different grazing intensities on soil physical properties. Acta Agrestia Sinica, 2002, 1(4): 74-78.

    [7]Du Y G, Cao G M, Wang Q L,etal. Effect of grazing on surface character and soil physical property in alpine meadow. Journal of Moutain Science, 2007, 25(3): 338-343.

    [8]Harald S, Atle M, Gunnar A. Sleep grazing and rodent populations: Evidence of negative interaction from a landscape scale experiment.Oecologa, 2005, 143(3): 357-364.

    [9]Cui Q H, Jiang Z G, Lian X M,etal. Factors influencing habitat selection of root voles (Microtusoeconomus). Acta Theriologica Sinica, 2005, 25(1): 45-51.

    [10]Fan N C, Wang Q Y, Zhou W Y,etal. The relationship between the plateau zokor population and vegetation damage. Proceedings of Alpine Meadow Ecosystem Congress[C]. Beijing: Science Press, 1988: 109-115.

    [11]Zhang Y M, Liu J K. The biological characteristics of subterranean rodents and their roles in ecosystem. Acta Theriologica Sinica, 2002, 22(5): 144-153.

    [12]Wu R N, Su B D, Bai M L. The influencing of subterranean rodents in agricultural ecosystem. Animal Husbandry and Feed Science, 2004, (6): 105-106.

    [13]Han T H, Zhang W G, Ding L S. Features of vegetation inMyospalaxbaileyihabitat. Acta Prataculturae Sinica, 1999, 8(2): 43-49.

    [14]Zhang J, Ge Q Z, Zhang W G,etal. Correlations between vegetation properties and fitness of plateau zokor habitat. Pratacultural Science, 2011, 28(5): 836-840.

    [15]Li F X, Zhang D G. Indicators and recovery approaches of degenerated grassland in China. Grassland and Turf, 2005, (1): 24-28.

    [16]He J L, Zhang J S, Yang Y B,etal. The study patterns and characteristics of plateau zokor mound. Acta Prataculturae Sinica, 2006, 15(1): 107-112.

    [17]Wu J, Wang J C. Principle and advantages of RTK. China Investigation Design, 2009, (9): 52-53.

    [18]Mannetje L T. The dry-weight-rank method for the botanical analysis of pasture. Grass and Forage Science, 1963, 18(4): 268-275.

    [19]Jones R M, Hargreaves J N G. Improvements to the dry-weight-rank method for measuring botanical composition. Grass and Forage Science, 1979, 34(3): 181-189.

    [20]Gillen R L, Smith E L. Evaluation of the dry-weight-rank method for determining species composition in tall grass prairie. Journal of Range Management, 1986, 5(7): 283-285.

    [21]Barnes D L, Odendaal J J, Beukes B H. Use of the dry-weight-rank method of botanical analysis in the eastern Transvaal Highveld. Proceedings of the Annual Congresses of the Grassland Society of Southern Africa, 1982, 17(1): 79-83.

    [22]Friedel M H, Chewings V H, Bastin G N. The use of comparative yield and dry-weight-rank techniques for monitoring arid rangeland. Journal of Range Management, 1988, 41(5): 430-435.

    [23]Dowhower S L, Teague W R, Ansley R J,etal. Dry-weight-rank method assessment in heterogenous communities. Journal of Range Management, 2001, 16(6): 71-76.

    [24]Dekker B, Kirkman K P, Du Plessis S I. Use of the dry-weight-rank method of botanical analysis in semi-arid savanna communities. African Journal of Range and Forage Science, 2001, 18(1): 63-66.

    [25]Hughes L E. Grazing systems:more thoughts and observations. Rangelands, 2011, 33(1): 35-40.

    [26]Rajput B S, Bhadauria U P S, Rajput R L. Study on specific leaf area and its relation with leaf area, leaf dry weight and seed yield in Mung bean. Journal of Life Sciences, 2014, 11(1): 247-250.

    [27]Haydock K P, Shaw N H. The comparative yield method for estimating dry matter yield of pasture. Animal Production Science, 1975, 15(76): 663-670.

    [28]Lv Y Z, Li B G. Edaphology[M]. Beijing: Chinese Agricultural Press, 2006.

    [29]Zhou J W, Hua L M, Zuo S T,etal. Research progress in habitat selection of plateau zokor. Pratacul Tural Science, 2013, 30(4): 647-653.

    [30]Liu W, Zhou L, Wang X. Responses of plant and rodents to different grazing intensity. Acta Ecologica Sinica, 1999, 19(3): 376-382.

    [31]Wen L Y, Li Z F. The effects of disturbance on maintaining mechanism of species diversity. Journal of Northwest Normal University (Natural Science), 2006, 42(4): 87-91.

    [32]Huang B. The Relationship between Plateau Zokor and Vegetation, Soil Factors[D]. Lanzhou: Lanzhou University, 2014.

    [33]Wang Q Y, Zhou W Y, Wei W H,etal. The burrowing behaviour myospalax bailey and its relation to soil hardness. Acta Theriologica Sinica, 2000, 20(4): 277-283.

    [34]Liu W G, Li L X, Xie H R,etal. Effect of soil bulk density on root morphology and biomass of vetiver grass seedlings. Acta Prataculturae Sinica, 2015, 24(4): 214-220.

    [35]Dai Q, Gu H J, Wang Y Z. Theories and models for habitat selection. Zoological Research, 2007, 28(6): 681-688.

    [36]Zhao Q S, Lou Y Q, Sun Y H. A review of statistical methods for evaluating animal habitat selection. Chinese Journal of Zoology, 2013, 48(5): 732-741.

    [37]Jin L R, Sun K P, He H S,etal. Research advances in habitat suitability index model. Chinese Journal of Ecology, 2008, 27(5): 841-846.

    [38]Bellis L M, Navarro J L, Vignolo P E,etal. Habitat preferences of lesser rheas in Argentine Patagonia. Biodiversity and Conservation, 2006, 15(9): 3065-3075.

    參考文獻(xiàn):

    [1]周雪榮, 郭正剛, 郭興華. 高原鼠兔和高原鼢鼠在高寒草甸中的作用. 草業(yè)科學(xué), 2010, 27(5): 38-44.

    [2]李午峰, 唐禎敏. 青藏鐵路的環(huán)境與安全. 安全與環(huán)境學(xué)報(bào), 2003, 3(3): 65-67.

    [3]樊乃昌, 施銀柱. 中國(guó)鼢鼠(Eospalax)亞屬分類研究. 獸類學(xué)報(bào), 1982, 2(1): 183-200.

    [6]張?zhí)N薇, 韓建國(guó), 李志強(qiáng). 放牧強(qiáng)度對(duì)土壤物理性質(zhì)的影響. 草地學(xué)報(bào), 2002, 1(4): 74-78.

    [7]杜巖功, 曹廣民, 王啟蘭, 等. 放牧對(duì)高寒草甸地表特征和土壤物理性狀的影響. 山地學(xué)報(bào), 2007, 25(3): 338-343.

    [9]崔慶虎, 蔣志剛, 連新明, 等. 田鼠棲息地選擇的影響因素. 獸類學(xué)報(bào), 2005, 25(1): 45-51.

    [10]樊乃昌, 王權(quán)業(yè), 周文揚(yáng), 等. 高原鼢鼠種群數(shù)量與植被破壞程度的關(guān)系-高寒草甸生態(tài)系統(tǒng)國(guó)際學(xué)術(shù)討論會(huì)文集[C]. 北京: 科學(xué)出版社, 1988: 109-115.

    [11]張堰銘, 劉季科. 地下鼠生物學(xué)特征及其在生態(tài)系統(tǒng)中的作用. 獸類學(xué)報(bào), 2002, 22(5): 144-153.

    [12]烏日娜, 蘇布道, 白夢(mèng)蘭. 地下鼠對(duì)農(nóng)業(yè)生態(tài)系統(tǒng)的影響. 畜牧與飼料科學(xué), 2004, (6): 105-106.

    [13]韓天虎, 張衛(wèi)國(guó), 丁連生. 高原鼢鼠棲息地植被特征. 草業(yè)學(xué)報(bào), 1999, 8(2): 43-49.

    [14]張軍, 葛慶征, 張衛(wèi)國(guó), 等. 植被性狀與高原鼢鼠棲息地適合度的關(guān)系. 草業(yè)科學(xué), 2011, 28(5): 836-840.

    [15]李鳳霞, 張德罡. 草地退化指標(biāo)及恢復(fù)措施. 草原與草坪, 2005, (1): 24-28.

    [16]何俊齡, 張金沙, 楊瑩博, 等. 高原鼢鼠土丘空間格局及主要特征研究. 草業(yè)學(xué)報(bào), 2006, 15(1): 107-112.

    [17]吳嬌, 王軍成. RTK工作原理及優(yōu)缺點(diǎn). 中國(guó)勘察設(shè)計(jì), 2009, (9): 52-53.

    [28]呂貽忠, 李保國(guó). 土壤學(xué)[M]. 北京: 中國(guó)農(nóng)業(yè)出版社出版, 2006.

    [29]周建偉, 花立民, 左松濤, 等. 高原鼢鼠棲息地選擇研究進(jìn)展. 草業(yè)科學(xué), 2013, 30(4): 647-653.

    [30]劉偉, 周立, 王溪. 不同放牧強(qiáng)度對(duì)植物及嚙齒動(dòng)物作用的研究. 生態(tài)學(xué)報(bào), 1999, 19(3): 376-382.

    [31]文隴英, 李仲芳. 干擾對(duì)物種多樣性維持機(jī)制的影響. 西北師范大學(xué)學(xué)報(bào), 2006, 42(4): 87-91.

    [32]黃彬. 高原鼢鼠種群密度與植被土壤因素的關(guān)系[D]. 蘭州: 蘭州大學(xué), 2014.

    [34]劉晚茍, 李良賢, 謝海蓉, 等. 土壤容重對(duì)野生香根草幼苗根系形態(tài)及其生物量的影響. 草業(yè)學(xué)報(bào), 2015, 24(4): 214-220.

    [35]戴強(qiáng), 顧海軍, 王躍招. 棲息地選擇的理論與模型. 動(dòng)物學(xué)研究, 2007, 28(6): 681-688.

    [36]趙青山, 樓瑛強(qiáng), 孫悅?cè)A. 動(dòng)物棲息地選擇評(píng)估的常用統(tǒng)計(jì)方法. 動(dòng)物學(xué)雜志, 2013, 48(5): 732-741.

    [37]金龍如, 孫克萍, 賀紅士, 等. 生境適宜都指數(shù)模型研究進(jìn)展. 生態(tài)學(xué)雜志, 2008, 27(5): 841-846.

    *A study of habitat selection by plateau zokor in paddocks with different grazing intensities in the eastern Qilian Mountain region

    CHU Bin, HUA Li-Min*, ZHOU Yan-Shan, LIU Li, LI Shuai

    CollegeofRangelandScienceofGansuAgricultureUniversity,Lanzhou730070,China

    Abstract:The objective of this study was to analyze the habitat selection and habitat characteristics of plateau zokor (Eospalaxbaileyi) in paddocks with different grazing intensities in the eastern Qilian Mountain region, in order to provide scientific information for controlling plateau zokor and conserving biodiversity in alpine meadows. Four paddocks with different grazing intensities were selected to investigate the zokor population density and habitat characteristics, including soil compaction and bulk density (g/cm3), vegetation botanical composition, plant species richness, aboveground plant biomass, and palatable forage biomass. A multiple regression analysis was used to determine the relationship between the plateau zokor population density and habitat characters. The zokor population density was highest in paddocks with the highest grazing intensity, and lowest in ungrazed paddocks (P>0.05). The zokor population density had a significant negative correlation with soil compaction (P<0.05) and a non-significant negative correlation with soil bulk density and aboveground plant biomass (P>0.05). The multiple regression for the relation between plateau zokor density and soil compaction, forb biomass, grass biomass, plants richness indicated that soil compaction was the dominant factor that influenced the plateau zokor habitat selection.

    Key words:plateau zokor; grazing intensity; habitat factors; habitat selection

    *通信作者Corresponding author. E-mail:hualm@gsau.edu.cn

    作者簡(jiǎn)介:楚彬(1991-),男,甘肅蘭州人,在讀碩士。E-mail:ryan_713@163.com

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(No.31460635),農(nóng)業(yè)部公益性行業(yè)科研項(xiàng)目(No.201203041)和甘肅省科技廳國(guó)際合作項(xiàng)目(No.1304WCGA174)資助。

    *收稿日期:2015-05-12;改回日期:2015-07-17

    DOI:10.11686/cyxb2015234

    http://cyxb.lzu.edu.cn

    楚彬,花立民,周延山,劉麗,李帥. 祁連山東段不同放牧強(qiáng)度下高原鼢鼠棲息地選擇分析. 草業(yè)學(xué)報(bào), 2016, 25(1): 179-186.

    CHU Bin, HUA Li-Min, ZHOU Yan-Shan, LIU Li, LI Shuai. A study of habitat selection by plateau zokor in paddocks with different grazing intensities in the eastern Qilian Mountain region. Acta Prataculturae Sinica, 2016, 25(1): 179-186.

    www.av在线官网国产| 欧美3d第一页| 国产国拍精品亚洲av在线观看| 天天操日日干夜夜撸| 少妇被粗大的猛进出69影院 | 我的老师免费观看完整版| 少妇被粗大的猛进出69影院 | 在线观看国产h片| 亚洲欧洲国产日韩| 国产精品一二三区在线看| 亚洲av欧美aⅴ国产| 国产精品不卡视频一区二区| 韩国av在线不卡| 国精品久久久久久国模美| 亚洲av男天堂| 久久免费观看电影| 久久久国产精品麻豆| 国产高清三级在线| 亚洲欧美中文字幕日韩二区| 久久ye,这里只有精品| 人人妻人人澡人人看| 中文天堂在线官网| 免费看光身美女| 乱人伦中国视频| 国产精品久久久久久精品电影小说| 久久久久人妻精品一区果冻| 久久99精品国语久久久| 99精国产麻豆久久婷婷| 91精品三级在线观看| 在线观看人妻少妇| 如何舔出高潮| 美女福利国产在线| 欧美性感艳星| 激情五月婷婷亚洲| 99热国产这里只有精品6| 亚洲精品aⅴ在线观看| 亚洲人成77777在线视频| 国产探花极品一区二区| 免费黄频网站在线观看国产| 在线观看www视频免费| 香蕉精品网在线| 美女视频免费永久观看网站| 日韩欧美精品免费久久| 免费av中文字幕在线| 69精品国产乱码久久久| 2022亚洲国产成人精品| 满18在线观看网站| 久久精品人人爽人人爽视色| 91精品一卡2卡3卡4卡| 中文字幕人妻熟人妻熟丝袜美| 99热全是精品| 春色校园在线视频观看| 久久ye,这里只有精品| 久久久久久伊人网av| 亚洲天堂av无毛| 永久免费av网站大全| 精品酒店卫生间| 欧美老熟妇乱子伦牲交| 国产又色又爽无遮挡免| 91精品国产九色| 在线观看www视频免费| 天堂俺去俺来也www色官网| 啦啦啦啦在线视频资源| 久久99热6这里只有精品| 999精品在线视频| av在线app专区| tube8黄色片| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久久久免| 国产乱人偷精品视频| 欧美日韩亚洲高清精品| 自线自在国产av| 国产高清三级在线| 国产又色又爽无遮挡免| 亚洲,欧美,日韩| 高清视频免费观看一区二区| 日韩电影二区| 尾随美女入室| kizo精华| 夜夜爽夜夜爽视频| 精品亚洲成a人片在线观看| 中文乱码字字幕精品一区二区三区| 制服人妻中文乱码| 国产精品秋霞免费鲁丝片| 成年美女黄网站色视频大全免费 | 三级国产精品欧美在线观看| 欧美最新免费一区二区三区| 国产欧美亚洲国产| 熟女av电影| 国产亚洲精品第一综合不卡 | a 毛片基地| 麻豆乱淫一区二区| 人妻少妇偷人精品九色| 99热6这里只有精品| 另类精品久久| 国产日韩欧美在线精品| 国产爽快片一区二区三区| 韩国av在线不卡| 国产精品一区www在线观看| 丰满少妇做爰视频| 国产精品久久久久久精品电影小说| 亚洲国产精品一区三区| 男女高潮啪啪啪动态图| 又粗又硬又长又爽又黄的视频| 精品一区二区免费观看| 在线观看免费日韩欧美大片 | 18禁在线无遮挡免费观看视频| 精品卡一卡二卡四卡免费| 亚洲,一卡二卡三卡| 日韩大片免费观看网站| 国产精品秋霞免费鲁丝片| 国产欧美日韩综合在线一区二区| 免费观看a级毛片全部| 欧美亚洲 丝袜 人妻 在线| 亚洲av.av天堂| 亚洲一级一片aⅴ在线观看| 亚洲av.av天堂| 十分钟在线观看高清视频www| 女人精品久久久久毛片| 亚洲精品乱久久久久久| 色吧在线观看| 精品午夜福利在线看| 一级片'在线观看视频| 五月玫瑰六月丁香| xxx大片免费视频| 欧美日韩在线观看h| .国产精品久久| 国产男女超爽视频在线观看| 看免费成人av毛片| 人人妻人人爽人人添夜夜欢视频| 一级爰片在线观看| 天天影视国产精品| 涩涩av久久男人的天堂| 在线观看一区二区三区激情| 午夜免费观看性视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 伊人久久国产一区二区| 亚洲欧美日韩另类电影网站| 午夜激情av网站| 午夜激情av网站| 九九爱精品视频在线观看| 国产爽快片一区二区三区| 极品少妇高潮喷水抽搐| 三上悠亚av全集在线观看| 纯流量卡能插随身wifi吗| 性色av一级| 天堂俺去俺来也www色官网| 亚洲人成网站在线播| 毛片一级片免费看久久久久| 国产精品女同一区二区软件| 欧美日韩精品成人综合77777| 国产欧美日韩综合在线一区二区| 精品国产一区二区久久| 日日撸夜夜添| 七月丁香在线播放| 18禁在线播放成人免费| 2018国产大陆天天弄谢| 国产精品秋霞免费鲁丝片| 亚洲综合色惰| 国产成人91sexporn| 一本—道久久a久久精品蜜桃钙片| 久久久久久伊人网av| av又黄又爽大尺度在线免费看| 免费久久久久久久精品成人欧美视频 | 只有这里有精品99| 91aial.com中文字幕在线观看| 亚洲欧美成人综合另类久久久| 最近最新中文字幕免费大全7| 一级a做视频免费观看| 国产精品三级大全| 啦啦啦中文免费视频观看日本| 2021少妇久久久久久久久久久| 97精品久久久久久久久久精品| 在线免费观看不下载黄p国产| 亚洲精品456在线播放app| 免费黄网站久久成人精品| 久久精品国产a三级三级三级| 91成人精品电影| 日韩欧美一区视频在线观看| 午夜福利,免费看| 夜夜看夜夜爽夜夜摸| 三级国产精品欧美在线观看| 国产成人精品无人区| 天天躁夜夜躁狠狠久久av| 大话2 男鬼变身卡| 在线观看免费日韩欧美大片 | 久久99热这里只频精品6学生| 国产亚洲最大av| 国精品久久久久久国模美| av国产久精品久网站免费入址| 99国产综合亚洲精品| 97超碰精品成人国产| 久久久久久久亚洲中文字幕| 精品人妻熟女毛片av久久网站| 大又大粗又爽又黄少妇毛片口| 亚洲av不卡在线观看| 成人综合一区亚洲| 欧美老熟妇乱子伦牲交| 十八禁网站网址无遮挡| av.在线天堂| 国产极品天堂在线| 欧美日韩精品成人综合77777| 我的女老师完整版在线观看| a级毛片免费高清观看在线播放| 日本av手机在线免费观看| 国产免费一区二区三区四区乱码| a级毛片黄视频| 最近的中文字幕免费完整| 中文字幕免费在线视频6| 精品酒店卫生间| av视频免费观看在线观看| 亚洲国产精品999| 国产成人精品在线电影| 亚洲精品视频女| 亚洲国产av影院在线观看| 日本猛色少妇xxxxx猛交久久| 久久青草综合色| 国产亚洲精品久久久com| 国产精品国产三级国产专区5o| 国产深夜福利视频在线观看| 欧美97在线视频| 伦理电影大哥的女人| 在线看a的网站| 亚洲精品日韩av片在线观看| 成人国语在线视频| 人人妻人人添人人爽欧美一区卜| 老司机亚洲免费影院| 青春草国产在线视频| 十八禁高潮呻吟视频| 80岁老熟妇乱子伦牲交| 成年人免费黄色播放视频| 最近中文字幕2019免费版| 久久久亚洲精品成人影院| 最新中文字幕久久久久| 国产男女内射视频| 免费观看无遮挡的男女| 国产精品嫩草影院av在线观看| 亚洲精品成人av观看孕妇| 日本91视频免费播放| 国产不卡av网站在线观看| 一级爰片在线观看| 欧美97在线视频| 国产在线免费精品| 国产免费一级a男人的天堂| 国产国语露脸激情在线看| 成人影院久久| 性色avwww在线观看| 国产欧美日韩一区二区三区在线 | 日韩中字成人| 丰满迷人的少妇在线观看| 国产片特级美女逼逼视频| 极品少妇高潮喷水抽搐| 黄色配什么色好看| 五月伊人婷婷丁香| 国产精品三级大全| 在线观看免费高清a一片| 午夜激情av网站| 欧美亚洲日本最大视频资源| 一个人免费看片子| 国产一区二区三区综合在线观看 | 人人妻人人爽人人添夜夜欢视频| 女人久久www免费人成看片| 成人毛片60女人毛片免费| 少妇人妻久久综合中文| 精品午夜福利在线看| 校园人妻丝袜中文字幕| 国模一区二区三区四区视频| 亚洲人成网站在线观看播放| 国产精品一二三区在线看| 久久精品国产鲁丝片午夜精品| 另类亚洲欧美激情| 日韩熟女老妇一区二区性免费视频| 只有这里有精品99| 亚洲国产av新网站| a级毛片免费高清观看在线播放| 一区二区三区免费毛片| 免费看av在线观看网站| 欧美97在线视频| .国产精品久久| 亚洲av日韩在线播放| 成人黄色视频免费在线看| 亚洲怡红院男人天堂| 18禁裸乳无遮挡动漫免费视频| √禁漫天堂资源中文www| 亚洲av免费高清在线观看| 国产精品一二三区在线看| 午夜影院在线不卡| 一个人看视频在线观看www免费| 久久热精品热| 高清av免费在线| 美女内射精品一级片tv| 久久精品国产亚洲网站| 成人综合一区亚洲| 国产探花极品一区二区| 99热这里只有是精品在线观看| 欧美日韩亚洲高清精品| 国产免费又黄又爽又色| 岛国毛片在线播放| 亚洲av国产av综合av卡| 九色成人免费人妻av| 精品人妻熟女av久视频| 国产一区有黄有色的免费视频| 看十八女毛片水多多多| av在线观看视频网站免费| 国产成人av激情在线播放 | 在线亚洲精品国产二区图片欧美 | 免费观看的影片在线观看| 秋霞伦理黄片| 精品久久久久久久久亚洲| 王馨瑶露胸无遮挡在线观看| 亚洲图色成人| 亚洲精品日韩在线中文字幕| 亚洲欧美日韩卡通动漫| 美女中出高潮动态图| 国产精品无大码| 午夜福利影视在线免费观看| 又黄又爽又刺激的免费视频.| 熟女电影av网| 精品一区二区免费观看| 久久青草综合色| 性色av一级| 热99久久久久精品小说推荐| 国产成人精品久久久久久| 国产一区二区在线观看av| av国产精品久久久久影院| 欧美3d第一页| 另类亚洲欧美激情| 欧美精品一区二区大全| 啦啦啦在线观看免费高清www| 热99国产精品久久久久久7| 十八禁高潮呻吟视频| 大香蕉久久网| 亚洲色图 男人天堂 中文字幕 | 国产白丝娇喘喷水9色精品| 如日韩欧美国产精品一区二区三区 | freevideosex欧美| √禁漫天堂资源中文www| 在现免费观看毛片| 夜夜爽夜夜爽视频| 成人毛片a级毛片在线播放| 国产片内射在线| 国产无遮挡羞羞视频在线观看| 岛国毛片在线播放| 亚洲怡红院男人天堂| 欧美另类一区| 嫩草影院入口| 久久精品夜色国产| 大又大粗又爽又黄少妇毛片口| 免费观看在线日韩| 成人国产av品久久久| 丰满饥渴人妻一区二区三| 亚洲欧美色中文字幕在线| 欧美精品高潮呻吟av久久| 亚洲经典国产精华液单| 国产高清不卡午夜福利| 国产 一区精品| 天堂8中文在线网| 欧美丝袜亚洲另类| 尾随美女入室| 在线观看免费日韩欧美大片 | 亚洲av成人精品一区久久| 国产熟女欧美一区二区| 99久久中文字幕三级久久日本| videossex国产| 精品亚洲乱码少妇综合久久| .国产精品久久| 又大又黄又爽视频免费| 中文字幕av电影在线播放| 国产精品99久久99久久久不卡 | 自线自在国产av| 波野结衣二区三区在线| 精品一区二区免费观看| 丰满迷人的少妇在线观看| 51国产日韩欧美| 久久久欧美国产精品| 亚洲av中文av极速乱| 亚洲国产精品一区三区| 免费看av在线观看网站| 久久久久久久久久久免费av| 国产一区亚洲一区在线观看| 999精品在线视频| 免费高清在线观看日韩| 一边亲一边摸免费视频| 午夜激情福利司机影院| 高清视频免费观看一区二区| 3wmmmm亚洲av在线观看| 毛片一级片免费看久久久久| 日本91视频免费播放| 亚洲精品aⅴ在线观看| 欧美另类一区| 丰满饥渴人妻一区二区三| 一个人看视频在线观看www免费| 久久精品国产自在天天线| 能在线免费看毛片的网站| 亚洲第一av免费看| 三级国产精品欧美在线观看| 有码 亚洲区| 热re99久久精品国产66热6| 五月开心婷婷网| 日日撸夜夜添| 成年人午夜在线观看视频| 91精品一卡2卡3卡4卡| 国产精品一国产av| 久久国产亚洲av麻豆专区| 亚洲图色成人| 成年美女黄网站色视频大全免费 | 久久久久久人妻| 99精国产麻豆久久婷婷| 男女边吃奶边做爰视频| 亚洲精品成人av观看孕妇| 久久久久国产网址| 精品人妻在线不人妻| 欧美xxⅹ黑人| 日韩强制内射视频| 久久久久久久久久久丰满| 欧美xxⅹ黑人| 人人澡人人妻人| videos熟女内射| 18在线观看网站| 国产免费一区二区三区四区乱码| 国产淫语在线视频| 国产精品无大码| 老司机亚洲免费影院| 久久av网站| 日本vs欧美在线观看视频| 51国产日韩欧美| 波野结衣二区三区在线| 国产精品欧美亚洲77777| 自线自在国产av| 亚洲精品亚洲一区二区| 卡戴珊不雅视频在线播放| 精品人妻在线不人妻| 亚洲国产av影院在线观看| 午夜激情福利司机影院| 精品国产露脸久久av麻豆| 欧美性感艳星| 午夜影院在线不卡| 欧美日韩精品成人综合77777| 国产片特级美女逼逼视频| 亚洲国产日韩一区二区| 欧美精品国产亚洲| 国产欧美日韩一区二区三区在线 | 高清不卡的av网站| 丝袜喷水一区| 国产精品无大码| 性色avwww在线观看| 91精品一卡2卡3卡4卡| 女人精品久久久久毛片| 三级国产精品片| 永久网站在线| 亚洲精品久久午夜乱码| tube8黄色片| 欧美日韩精品成人综合77777| 一本一本综合久久| 久久久国产精品麻豆| 精品一区在线观看国产| 久久精品国产鲁丝片午夜精品| 亚洲精品久久久久久婷婷小说| 欧美 亚洲 国产 日韩一| 国产欧美亚洲国产| 欧美少妇被猛烈插入视频| 久久久久久久亚洲中文字幕| 久久精品久久久久久噜噜老黄| 欧美一级a爱片免费观看看| av播播在线观看一区| 欧美日韩国产mv在线观看视频| 日韩成人伦理影院| 国产男人的电影天堂91| 欧美成人精品欧美一级黄| 欧美精品一区二区大全| 热re99久久国产66热| 熟女av电影| 久久人人爽人人爽人人片va| 国产乱人偷精品视频| 亚洲美女黄色视频免费看| 精品人妻熟女av久视频| 国产av精品麻豆| 欧美日韩视频高清一区二区三区二| 伊人久久国产一区二区| 日韩中字成人| 丝袜在线中文字幕| 最新的欧美精品一区二区| 午夜激情av网站| 日本黄大片高清| 国语对白做爰xxxⅹ性视频网站| 免费大片黄手机在线观看| 青青草视频在线视频观看| 亚洲欧美成人精品一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久人妻| 亚洲精品,欧美精品| 黄色视频在线播放观看不卡| 插逼视频在线观看| 久久 成人 亚洲| 人人妻人人爽人人添夜夜欢视频| 日韩强制内射视频| 人妻 亚洲 视频| 午夜激情久久久久久久| 日日啪夜夜爽| 观看av在线不卡| 亚洲精品成人av观看孕妇| 亚洲精品视频女| 久久久久精品性色| 精品熟女少妇av免费看| 亚洲欧美清纯卡通| 性高湖久久久久久久久免费观看| 9色porny在线观看| 国产成人午夜福利电影在线观看| 日韩av免费高清视频| 久久久a久久爽久久v久久| 极品人妻少妇av视频| 18禁观看日本| 久久久久精品久久久久真实原创| 久久99热6这里只有精品| 人成视频在线观看免费观看| 国产在视频线精品| 少妇高潮的动态图| 国产综合精华液| 国产不卡av网站在线观看| 亚洲综合色网址| 中文字幕人妻熟人妻熟丝袜美| 国产极品粉嫩免费观看在线 | 精品一区二区三卡| 日日撸夜夜添| 国产乱人偷精品视频| 黄片无遮挡物在线观看| 曰老女人黄片| 欧美日韩一区二区视频在线观看视频在线| 国产成人免费无遮挡视频| 国产毛片在线视频| 国产视频首页在线观看| 日韩一区二区视频免费看| 在线观看三级黄色| 国产高清三级在线| 亚洲综合色惰| 久久精品国产鲁丝片午夜精品| 最新中文字幕久久久久| 香蕉精品网在线| 亚洲熟女精品中文字幕| 日本黄色片子视频| 亚洲精品久久久久久婷婷小说| 亚洲人与动物交配视频| 男的添女的下面高潮视频| 国产老妇伦熟女老妇高清| 亚洲精品aⅴ在线观看| 色婷婷久久久亚洲欧美| 在线天堂最新版资源| 在线观看免费高清a一片| 一个人看视频在线观看www免费| 国产日韩一区二区三区精品不卡 | 国产精品 国内视频| 免费黄色在线免费观看| 男女啪啪激烈高潮av片| 亚洲欧美一区二区三区国产| 国产黄色免费在线视频| 人妻少妇偷人精品九色| 制服丝袜香蕉在线| 老司机亚洲免费影院| 美女大奶头黄色视频| 99热6这里只有精品| 日产精品乱码卡一卡2卡三| 久久久国产精品麻豆| 一二三四中文在线观看免费高清| av有码第一页| 少妇猛男粗大的猛烈进出视频| 看十八女毛片水多多多| 精品卡一卡二卡四卡免费| 九九在线视频观看精品| 大香蕉97超碰在线| 男的添女的下面高潮视频| 狂野欧美白嫩少妇大欣赏| 亚洲精品久久久久久婷婷小说| 国产视频内射| 日本av免费视频播放| 亚洲国产欧美日韩在线播放| 看非洲黑人一级黄片| 高清在线视频一区二区三区| 在线观看三级黄色| 国产精品国产av在线观看| 久久午夜综合久久蜜桃| 久久久久久久久久久免费av| 能在线免费看毛片的网站| 插逼视频在线观看| 国产精品99久久99久久久不卡 | 亚洲熟女精品中文字幕| 国产欧美日韩一区二区三区在线 | 亚洲少妇的诱惑av| 一级毛片 在线播放| 久久久久久久亚洲中文字幕| √禁漫天堂资源中文www| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 建设人人有责人人尽责人人享有的| 国产熟女午夜一区二区三区 | 亚洲欧洲国产日韩| 丁香六月天网| 又粗又硬又长又爽又黄的视频| 久久久欧美国产精品| 卡戴珊不雅视频在线播放| 人妻夜夜爽99麻豆av| 丰满饥渴人妻一区二区三| 91aial.com中文字幕在线观看| 欧美 日韩 精品 国产| 日本vs欧美在线观看视频| 一级爰片在线观看| 国产精品三级大全| 人妻一区二区av| 黄色毛片三级朝国网站| 国产高清不卡午夜福利| 日韩亚洲欧美综合| 日韩av在线免费看完整版不卡| 晚上一个人看的免费电影| 69精品国产乱码久久久| 国产亚洲av片在线观看秒播厂| 免费av中文字幕在线| 91精品三级在线观看|