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

    甘肅敦煌西湖荒漠-濕地生態(tài)系統(tǒng)土壤水分空間異質(zhì)性及其影響因子研究

    2016-03-23 02:29:32孫飛達(dá)陳文業(yè)袁海峰竇英杰
    西北植物學(xué)報(bào) 2016年1期
    關(guān)鍵詞:荒漠土壤水分

    孫飛達(dá),陳文業(yè),袁海峰,竇英杰,

    邴丹琿2,3,4,馮 穎2,3,4,吳 婷4,5

    (1 四川農(nóng)業(yè)大學(xué) 草業(yè)科學(xué)系,成都 611130;2 甘肅省林業(yè)科學(xué)研究院,蘭州 730020;3 甘肅林研科技工程公司,蘭州 730020;4 甘肅敦煌西湖國家級自然保護(hù)區(qū)管理局,甘肅敦煌 736200;5 甘肅敦煌西湖濕地生態(tài)系統(tǒng)國家定位觀測研究站,甘肅敦煌 736200)

    ?

    甘肅敦煌西湖荒漠-濕地生態(tài)系統(tǒng)土壤水分空間異質(zhì)性及其影響因子研究

    孫飛達(dá)1,陳文業(yè)2,3,4*,袁海峰4,5,竇英杰4,5,

    邴丹琿2,3,4,馮穎2,3,4,吳婷4,5

    (1 四川農(nóng)業(yè)大學(xué) 草業(yè)科學(xué)系,成都 611130;2 甘肅省林業(yè)科學(xué)研究院,蘭州 730020;3 甘肅林研科技工程公司,蘭州 730020;4 甘肅敦煌西湖國家級自然保護(hù)區(qū)管理局,甘肅敦煌 736200;5 甘肅敦煌西湖濕地生態(tài)系統(tǒng)國家定位觀測研究站,甘肅敦煌 736200)

    摘要:土壤水分是內(nèi)陸荒漠區(qū)濕地生態(tài)系統(tǒng)中重要的限制因子,為了揭示該區(qū)域土壤水分空間分布特征,采用傳統(tǒng)統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)相結(jié)合的方法,對甘肅敦煌西湖國家級自然保護(hù)區(qū)0~200 cm內(nèi)各層土壤水分的空間變異性及海拔、土壤質(zhì)地和植被對其的影響進(jìn)行了研究,旨在為極干旱區(qū)濕地生態(tài)系統(tǒng)植被修復(fù)和保育提供科學(xué)依據(jù)。結(jié)果表明:(1)本研究所得各變量的變異系數(shù)、塊金方差、基臺(tái)值、變程和結(jié)構(gòu)比分別為36.51%~88.65%、0.007~0.098、0.112~0.549、116~453和76.6%~97.6%,各變量均為中等變異,存在高度異質(zhì)性,具有較強(qiáng)空間自相關(guān)。(2)深層(60~200 cm)土壤水分含量較淺層(0~60 cm)變異大,且不同層次土壤水分含量的空間異質(zhì)性差別也較大,空間變異主要發(fā)生在較小尺度上(分維數(shù)D在1.902~1.989之間)。(3)海拔是影響該區(qū)域深層土壤水分空間變異的主導(dǎo)因子。(4)土壤質(zhì)地與淺層(0~60 cm)土壤水分含量的相關(guān)性大于與深層(60~200 cm)土壤含水量的相關(guān)性,它們與海拔相關(guān)性表現(xiàn)相反;草本植被蓋度與淺層土壤水分含量呈較高的正關(guān)聯(lián)關(guān)系,灌木根量與深層土壤水分含量呈較高的負(fù)關(guān)聯(lián)關(guān)系。

    關(guān)鍵詞:敦煌西湖;荒漠-濕地生態(tài)系統(tǒng);土壤水分;空間異質(zhì)性;環(huán)境因子影響

    空間異質(zhì)性是指系統(tǒng)或系統(tǒng)屬性在空間上的復(fù)雜性和變異程度,包括系統(tǒng)屬性的空間組成,空間構(gòu)型和空間相關(guān)[1-3]。土壤水分空間異質(zhì)性是土壤重要屬性之一[3-4],水在干旱半干旱地區(qū)極具敏感性,是生態(tài)系統(tǒng)中最活躍的因素,是反映土壤特性的重要指標(biāo)[5],對該區(qū)域植被恢復(fù)與建設(shè)以及經(jīng)濟(jì)持續(xù)發(fā)展具有極大的限制性[6]。在不同尺度上研究土壤水分空間異質(zhì)性,對了解植被與土壤水分的關(guān)系具有重要的參考價(jià)值[3,7]。土壤水分的空間異質(zhì)性研究一直是國際上水文學(xué)和土壤學(xué)研究的熱點(diǎn)問題,隨著地統(tǒng)計(jì)學(xué)方法的出現(xiàn),采用該方法對土壤水分空間變異性的研究逐漸增多[8],國內(nèi)外很多學(xué)者對不同類型生態(tài)系統(tǒng)的土壤水分空間異質(zhì)性做了大量的研究工作[9],如高寒區(qū)[10]、荒漠綠洲區(qū)[11]、農(nóng)牧交錯(cuò)帶[12]、黃土區(qū)[1]、荒漠區(qū)[13-14]、喀斯特地貌區(qū)[15]等,但整體來看,研究對象的類型并不多,尺度范圍也較窄,對不同深度層土壤水分變異性的研究較少[9]。尤其是對荒漠區(qū)濕地生態(tài)系統(tǒng)土壤水分空間變異性研究則鮮見報(bào)道。

    甘肅敦煌西湖國家級自然保護(hù)區(qū)(以下簡稱敦煌西湖)地處甘肅河西走廊最西端,西接庫姆塔格沙漠和羅布泊,南接阿克塞哈薩克族自治縣,北連新疆維吾爾自治區(qū)。該區(qū)域濕地屬于沼澤濕地,是內(nèi)陸干旱地區(qū)的典型濕地類型,具有極干旱區(qū)濕地生態(tài)系統(tǒng)和荒漠生態(tài)系統(tǒng)的典型性和代表性,區(qū)位優(yōu)勢明顯、特殊。土壤水分是該區(qū)域植被生長重要的限制因子,調(diào)控著植被的組成、分布、結(jié)構(gòu)及其系統(tǒng)穩(wěn)定性。本研究采用了傳統(tǒng)統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)相結(jié)合的方法,定量研究了敦煌西湖0~200 cm土層不同層次水分的空間異質(zhì)性,以期進(jìn)一步深入研究研究區(qū)植被與土壤水分的相互關(guān)系,旨在了解研究區(qū)土壤水分空間分布特征,為極干旱區(qū)濕地生態(tài)系統(tǒng)植被修復(fù)和保育提供科學(xué)依據(jù)。

    1研究區(qū)自然條件

    敦煌西湖面積為6.6×105hm2,其中濕地面積9.80×104hm2,蘆葦沼澤3.43×104hm2,四周均被沙漠和戈壁所隔絕。地理坐標(biāo)為:92°45′~93°50′E,39°45′~40°36′N,區(qū)內(nèi)海拔820~2 359 m,地勢南高北低,中間為沖積平原。地處北半球暖溫帶干旱氣候區(qū),屬典型的大陸性氣候,年平均氣溫為9.90 ℃,最低氣溫-30 ℃,最高氣溫40 ℃;年均降水量39.90 mm,蒸發(fā)量2 486 mm;年均風(fēng)速2.20 m/s,大風(fēng)日數(shù)15.40 d;年日照時(shí)數(shù)為3 115.0~3 246.7 h,日照率70%~73%;年總輻射量為641.84 kJ/cm2,干燥度大于16。植物區(qū)系歸屬于泛北極植物區(qū)中的亞洲荒漠植物亞區(qū),具有鮮明的溫帶荒漠性質(zhì),敦煌西湖共記錄有種子植物23科61屬80種。土壤主要為沼澤土,另有部分草甸土分布于河漫灘等地。主要保護(hù)對象為濕地生態(tài)系統(tǒng)、荒漠生態(tài)系統(tǒng)及其野生動(dòng)植物[16]。

    2研究方法

    2.1樣地布設(shè)

    從鹽池灣開始,向南設(shè)置長1.6 km、寬300 m的調(diào)查樣帶,植被包括沼澤、鹽沼、草甸、闊葉林荒漠5種植被型組。在樣帶內(nèi),沿樣帶方向設(shè)置3條平行樣線,按南北、東西間隔100 m取樣,共有450個(gè)樣點(diǎn),并以每個(gè)樣點(diǎn)為中心設(shè)置50 m×50 m樣地。按照“五點(diǎn)法”在每個(gè)樣地內(nèi)設(shè)置5根直徑50 mm、深度2.2 m的土壤水分測定管。

    2.2土壤水分含量測定

    土壤含水量采用CNC503B(DR)新型智能中子水分儀測定。按0~10、10~20、20~40、40~60、60~80、80~100、100~120、120~140、140~160、160~180、180~200 cm深度分層測定;每個(gè)樣地每層測定5個(gè)樣點(diǎn),取其平均值作為樣地該次該層的土壤含水量。最后將算得的土壤水分含量和土壤貯水量值用于空間變異性分析。土壤貯水量的計(jì)算公式為:S=0.1hmd,式中:S為土壤儲(chǔ)水量(mm)、h為土層厚度(20 cm)、m為質(zhì)量含水量(%)、d為該土層土壤容重(環(huán)刀法,環(huán)刀規(guī)格:直徑50.46 mm×高50 mm,容積100 cm3)[17]。依據(jù)該公式分別計(jì)算出各層的土壤貯水量,然后累加計(jì)算出2 m土層的水分貯量。

    2.3植被調(diào)查

    每個(gè)樣地采用樣方與樣線相結(jié)合的方法進(jìn)行植物的種類、密度、蓋度、頻度和高度等調(diào)查[16],在交叉樣線處均勻設(shè)置6個(gè)20 m×20 m樣方進(jìn)行灌(喬)木調(diào)查,并設(shè)置9個(gè)5 m×5 m樣方進(jìn)行草本植物調(diào)查。

    2.4地下水埋深調(diào)查

    每個(gè)樣地用土鉆打至地下水水流出,待水位穩(wěn)定后測量地下水埋深。

    2.5根系生物量的調(diào)查

    采用全挖法。①灌木調(diào)查:在樣方中選擇大小相近、生長旺盛的植株5株,離植株20 cm處垂直向下挖一剖面,把根系分布層分為0~60 cm和60~200 cm 2層;②草本調(diào)查:在樣方內(nèi)選擇地勢平坦、植被相對集中的區(qū)域,由上向下按30 cm×50 cm×20 cm挖取土柱,深度以無植株根系為限,直至200 cm深處。剪取各層所有根系,分層裝袋,在室內(nèi)沖洗、晾曬、烘干后稱重。

    2.6土壤粒徑分析

    每個(gè)樣地每層采集5個(gè)樣點(diǎn)土樣充分混合(表層土壤直接取樣,深層采用挖剖面和土鉆法相結(jié)合的方式取樣),用四分法取大約1 kg樣品,室內(nèi)風(fēng)干,對土壤粒徑用MS-S激光粒度分析儀(英國馬爾文Malvem儀器有限公司)進(jìn)行分析[18]。

    為了便于分析比較,主要考慮不受降水影響,每年的9、10月份是當(dāng)?shù)亟邓肯鄬ψ钌俚臅r(shí)期,加之該時(shí)期研究區(qū)植被生長旺盛,所以在2013年的10月5日至10月25日期間對研究區(qū)植被特征與土壤特性調(diào)查及土樣采集同步進(jìn)行。在調(diào)查的同時(shí),用GPS定位,并記錄各樣地的經(jīng)緯度、海拔高度、地貌及土壤類型等生境因子。

    2.7分析方法

    采用GS+進(jìn)行地統(tǒng)計(jì)學(xué)變異函數(shù)數(shù)據(jù)分析。統(tǒng)計(jì)分析分兩步驟:①對土壤水分等數(shù)據(jù)用單樣本方法進(jìn)行正態(tài)分布檢驗(yàn)(SPSS 13.0軟件);②變異函數(shù)的計(jì)算,用于估計(jì)半方差的公式為[19-21]:

    (1)

    式中,r(h)為半方差函數(shù);N(h)是距離等于h時(shí)的點(diǎn)對數(shù),Zi是樣點(diǎn)Z在位置i的實(shí)測值,Z(i+h)是與i距離為h處樣點(diǎn)的值。

    本研究中球狀模型(公式2)和指數(shù)模型(公式3)土壤水分等數(shù)據(jù)符合要求。球狀模型的變程等于相關(guān)距離a;而指數(shù)模型并不表現(xiàn)出有限變程,但實(shí)踐中變程值近似用3a表示。

    (2)

    (3)

    式中,C0為塊金值(nugget),C0+C為基臺(tái)值,a為相關(guān)距離。

    分維數(shù)D的計(jì)算由變異函數(shù)和步長h之間的關(guān)系確定,即:

    2r(h)=h4-2D

    (4)

    上式取雙對數(shù)后再對雙對數(shù)曲線進(jìn)行線性回歸,得到回歸直線的斜率k,分維數(shù)可用斜率k估算:

    D=1/2(4-k)

    (5)

    分維數(shù)D的大小可用于度量生態(tài)變量的復(fù)雜程度,D值越小,變量的空間依賴性越強(qiáng),其空間格局相對簡單。通過分形維數(shù)可分析生態(tài)因子在不同尺度上的差異[19-20]。

    3結(jié)果與分析

    3.1土壤水分的空間變異性特征

    3.1.1土壤水分狀況表1顯示,垂直方向上,各樣點(diǎn)土壤水分含量均隨土層加深而增大,同一樣點(diǎn)相鄰?fù)翆油寥浪趾孔兓淮?土壤水分含量均值為14.873%,介于6.681%~23.004%之間,2 m土層水分貯量平均值為163.598 mm,介于65.588~248.188 mm之間。一般認(rèn)為,CV<0.1為弱變異性,0.1≤CV≤1為中等變異性,CV>1為高度變異性[22],土壤水分含量變異系數(shù)(CV)介于39.28%~88.65%,隨土層加深顯逐漸增大趨勢,說明研究區(qū)各樣點(diǎn)土壤剖面各層土壤水分含量的空間變異性在增大。2 m土層水分貯量變異系數(shù)(CV)為36.51%(表1),變異系數(shù)顯示,研究區(qū)土壤水分含量和2 m土層水分貯量變異相對較高,均為中等變異。深層(60~200 cm)土壤水分含量的變異較淺層(0~60 cm)土壤水分含量的變異大,說明深層土壤具有較高的土壤水分含量和較大的變異程度,而淺層土壤水分狀況較差且變異程度較小。

    3.1.2土壤水分空間結(jié)構(gòu)特征根據(jù)研究區(qū)野外調(diào)查數(shù)據(jù)計(jì)算實(shí)際變異函數(shù),分別用不同類型的模型進(jìn)行擬合,得到最佳擬合模型的參數(shù)值見表2,結(jié)果顯示,0~60 cm各土層土壤水分含量、2 m土層貯水量和植被蓋度的理論模型符合指數(shù)模型,而60~200 cm各土層土壤水分含量的理論模型符合球狀模型。

    垂直方向上,11層土層土壤水分含量及2 m土層貯水量的塊金方差(C0)介于0.007~0.098之間,表明各層土壤水分含量及2 m土層貯水量存在高度的異質(zhì)性,這種異質(zhì)性在各要素塊金值中也可以看出,各要素的基臺(tái)值(C0+C)介于0.112~0.549之間。0~60 cm各土層土壤水分含量的塊金方差和基臺(tái)值較60~200 cm土層的低,且各值均差別不大,而60~200 cm各土層土壤水分含量的塊金方差和基臺(tái)值差別均較大,表明60~200 cm各土層土壤水分含量隨機(jī)部分的空間異質(zhì)性和總的空間異質(zhì)性程度較高,且不同層次土壤含水量的空間異質(zhì)性差別也較大。本研究所得各變量的結(jié)構(gòu)比介于0.766~0.976之間(表2),表明研究區(qū)各土層土壤水分含量和2 m土層貯水量具有較強(qiáng)的空間自相關(guān)。整體上看,研究區(qū)2 m土層貯水量的結(jié)構(gòu)比最小(0.766),說明對于2 m土層貯水量由隨機(jī)因素引起的空間異質(zhì)性占總空間異質(zhì)性的比例比其他變量大[23-25]。

    變程是表示植被特征與土壤水分空間異質(zhì)性的尺度[26-27]。表2顯示,研究區(qū)各土層土壤水分含量和2 m土層貯水量異質(zhì)性尺度不同且值均較小,變程介于116~453,2 m土層貯水量與各層土壤水分含量變程(a)比較,相對較大。

    D表示變異函數(shù)曲線的曲率大小[27]。表2顯示,0~60 cm各土層土壤水分含量、2 m土層貯水量的分形維數(shù)(D)較小(1.883~1.934),說明空間變異主要發(fā)生在較大尺度上,而60~200 cm各土層土壤水分含量的分形維數(shù)(D)較大(1.902~1.989),說明空間變異主要發(fā)生在較小尺度上,即沿樣帶相鄰點(diǎn)間的土壤水分含量差異很大。

    3.2研究區(qū)土壤水分空間異質(zhì)性影響因素分析

    3.2.1海拔用GPS測量研究區(qū)所有樣點(diǎn)的海拔高度,其中最高點(diǎn)海拔為1 050.38 m,最低點(diǎn)海拔為880.43 m。圖1顯示,隨著海拔遞增,各土層土壤水分含量和2 m土層貯水量均呈下降趨勢,2 m土層貯水量變化趨勢較明顯。

    海拔是影響土壤水分含量分布的環(huán)境因子之一,許多研究已表明土壤水分含量與海拔之間具有負(fù)的相關(guān)關(guān)系[18,28-29]。圖2顯示,淺層土壤水分含量(10~60cm)和深層土壤水分含量(60~200 cm)均與海拔有負(fù)的相關(guān)關(guān)系(P<0.01),隨著海拔遞增,深層土壤水分含量與海拔的相關(guān)性呈遞增趨勢且更高,而淺層土壤水分含量則變化相反,且海拔與深層(60~200 cm)土壤水分含量的相關(guān)性大于淺層(0~60 cm)土壤濕度的相關(guān)性。

    表1 研究區(qū)剖面各層土壤水分含量的統(tǒng)計(jì)特征

    表2 土壤水分含量變異函數(shù)理論模型及參數(shù)

    圖1 土壤水分含量隨海拔梯度變化

    圖2 土壤水分含量與海拔間的相關(guān)性特征

    3.2.2土壤質(zhì)地土壤質(zhì)地影響著土壤水分運(yùn)移,進(jìn)而影響著土壤水分異質(zhì)性程度。表3顯示,粘粒和粉粒的變異系數(shù)均大于沙粒,且沿著樣線,隨著海拔增大,沙粒的含量表現(xiàn)出明顯遞增趨勢,而粘粒和粉粒的含量呈現(xiàn)出相反的空間分布特征。土壤質(zhì)地的空間異質(zhì)性能夠產(chǎn)生水文導(dǎo)度的空間變異,從而導(dǎo)致土壤水分分布的空間異質(zhì)性[18]。

    由表3可知,粉粒和粘粒均與各層土壤水分含量之間呈正相關(guān)關(guān)系,且粘粒的相關(guān)性大于粉粒的相關(guān)性;而沙粒與各層土壤水分含量之間呈負(fù)相關(guān)關(guān)系。因?yàn)?土壤粒徑粗細(xì)程度對土壤質(zhì)地產(chǎn)生著重大影響,粒徑均值越小,土壤可塑性和吸濕性越明顯,其透水性較差,持水性越好。土壤質(zhì)地與淺層(0~60 cm)土壤水分含量的相關(guān)性大于深層(60~200 cm)土壤的濕度的相關(guān)性,這結(jié)論與潘顏霞等[18]的研究結(jié)果一致。因?yàn)闇\層土壤中粉粒級以下的顆粒在風(fēng)蝕過程中被吹失,呈現(xiàn)出隨時(shí)間推移變粗的特征,調(diào)查結(jié)果顯示,研究區(qū)淺層(0~60 cm)土壤的沙粒含量為65%,而深層(60~200 cm)土壤的粉粒和粘粒含量接近80%。

    3.2.3植被植被通過遮蓋土壤表面影響蒸發(fā)率,通過根活性影響土壤導(dǎo)度并促進(jìn)土壤表面有機(jī)物質(zhì)的增加而影響土壤水分變化,這些因素對于土壤水分含量變化率的影響隨著植被種類、密度和蓋度變化而變化[18,30]。圖3顯示,研究區(qū)主要植被群落土壤水分含量均相對較低。其中,多枝檉柳沙包群落土壤水分含量最少,均值為6.681%;蘆葦沼澤群落土壤水分含量最大,均值為23.004%。各主要植被群落淺層(0~60 cm)土壤同一層土壤水分含量差異相對較少(多枝檉柳群落和蘆葦沼澤群落;圖3,A),深層(60~200 cm)土壤同一層土壤水分含量差異較大(圖3,B)。不同植被類型及其格局分布影響著土壤水分,植被根系垂直和水平分布格局通過蒸騰耗水影響著其根際層的水分,從而影響到各層土壤水分含量異質(zhì)性[9]。草本植物根系較淺,根系生物量集中分布在40 cm以上土層中,占全部根系生物量的85%以上,對淺層水分影響較大;灌木根系較深,大部分根系分布于60 cm以下土層中,占全部根系生物量的70%以上,改變了土壤質(zhì)地和有機(jī)質(zhì)含量,影響了水文導(dǎo)度,對深層水分影響較大。表4顯示,植被蓋度與土壤水分含量顯正相關(guān),且草本植被蓋度與淺層土壤水分含量顯示較高的關(guān)聯(lián)性,是由于較大的植被蓋度可以阻擋太陽照射,減少土壤的蒸發(fā),改善了表層的土壤有機(jī)質(zhì),往往導(dǎo)致較高的土壤水分;植被根量與土壤水分含量顯負(fù)相關(guān),且灌木根量與深層土壤水分含量顯示較高的關(guān)聯(lián)性,是由于植被根系集中分布層的根系對該層土壤水分含量具有顯著的影響,尤其在植被生長旺盛季節(jié),植被蒸騰作用增強(qiáng),根系吸水強(qiáng)烈,通常導(dǎo)致根系主要分布層的土壤水分含量大幅度降低。

    表3 土壤粒徑與土壤水分含量之間相關(guān)性及其空間變異性

    a.研究區(qū)主要植被淺層土壤水分含量;b.研究區(qū)主要植被深層土壤水分含量及其2 m土層貯水量;

    植被Vegetation土壤層Soillayer/cm0~6060~200草本Herbs蓋度Cover0.67180.3627根重Rootweight-0.5438-0.2567灌木Shrubs蓋度Cover0.42670.3833根重Rootweight-0.5120-0.6633

    表5 不同土層土壤水分含量、2 m土層

    3.2.4地下水埋深調(diào)查發(fā)現(xiàn)地下水埋深較淺的區(qū)域,土壤水分含量較高,否則則反之。表5顯示,土壤含水量與地下水埋深之間存在顯著的負(fù)相關(guān)關(guān)系,自上而下其關(guān)系逐漸增強(qiáng);2 m土層貯水量與地下水埋深的相關(guān)性最強(qiáng),40~60 cm土壤含水量與地下水埋深的相關(guān)性次之,說明研究區(qū)土壤水分含量更依賴于地下水,地下水埋深是影響土壤含水量的重要因素。

    4 討論

    敦煌西湖國家級自然保護(hù)區(qū),四周均被沙漠、戈壁所隔絕,研究區(qū)兼有濕地生態(tài)系統(tǒng)和荒漠生態(tài)系統(tǒng),導(dǎo)致本研究所得各變量均為中等變異,存在高度的異質(zhì)性,具有較強(qiáng)的空間自相關(guān),0~200 cm土層各層土壤水分含量及2 m土層貯水量的變異系數(shù)、塊金方差、基臺(tái)值、變程和結(jié)構(gòu)比分別介于36.51%~88.65%、0.007~0.098、0.112~0.549、116~453和76.6%~97.6%之間。深層(60~200 cm)土壤水分含量的變異較淺層(0~60 cm)土壤水分含量的變異大,且不同層次土壤水分含量的空間異質(zhì)性差別也較大,空間變異主要發(fā)生在較小尺度上。說明深層土壤較高的土壤濕度和較大的變異程度,而淺層土壤水分狀況較差且變異程度較小,這是因?yàn)楸韺虞^底層土壤蒸發(fā)強(qiáng)烈,土壤水分含量較低,且與剖面土壤質(zhì)地不均勻、植被種類、蓋度和多度有關(guān),加之研究區(qū)主要是以旱生、超旱生植被為主,通過其發(fā)達(dá)的根系對深層土壤水分的利用較大,進(jìn)而導(dǎo)致深層土壤水分變異系數(shù)大于淺層。

    參考文獻(xiàn):

    [1]畢華興,李笑吟,劉鑫,等.晉西黃土區(qū)土壤水分空間異質(zhì)性的地統(tǒng)計(jì)學(xué)分析[J].北京林業(yè)大學(xué)學(xué)報(bào),2006,28(5):59-66.

    BI H X,LI X Y,LIU X,etal.Spatial heterogeneity of soil moisture using geological statistics method in the loess region of west Shanxi Province[J].JournalofBeijingForestryUniversity,2006,28(5):59-66.

    [2]MLLLER R E,VEF H J M,FOWLER N L.Spatial heterogeneity in eight central Teaxs gtasslands[J].JournalofEcology,1995,83:919-928.

    [3]郭德亮,樊軍,米美霞.黑河中游綠洲區(qū)不同土地利用類型表層土壤水分空間變異的尺度效應(yīng)[J].應(yīng)用生態(tài)學(xué)報(bào),2013,24(5):1 199-1 208.

    GUO D L,FAN J,MI M X.Scale-dependency of spatial variability of surface soil moisture under different land use types in Heihe Oasis,China [J].ChineseJournalofAppliedEcology,2013,24(5):1 199-1 208.

    [4]ROSSI R E,MULLA D J,JOUREL A G,etal.Geostatistical tools for mdeling and interpreting ecological spatial dependence[J].EcologicalMonographs,1991,62:277-314.

    [5]李朝生,楊曉暉,張克斌,等.沙漠-綠洲系統(tǒng)中降雨、土壤水分與地下水位的響應(yīng)特征[J].北京林業(yè)大學(xué)學(xué)報(bào),2007,29(4):129-135.

    LI C S,YANG X H,ZHANG K B,etal.Response characteristics of precipitation,soil moisture and groundwater level in desert oasis system[J].JournalofBeijingForestryUniversity,2007,29(4):129-135.

    [6]溫從.寧夏鹽池沙地主要植物群落土壤水分動(dòng)態(tài)研究[D].北京:北京林業(yè)大學(xué),2007:2-12.

    [7]LEGENDER P,FORTIN M J.Spatial pattern and ecological analysis[J].Vegetation,1989,80:107-138.

    [8]白永飛,許志信,李德新,等.內(nèi)蒙古高原針茅草原群落土壤水分和碳、氮分布的小尺度空間異質(zhì)性[J].生態(tài)學(xué)報(bào),2002,22(8):1 215-1 223.

    BAI Y F,XU Z X,LI D X,etal.On the small scale spatial heterogeneity of soil moisture,carbon and nitrogen in sipa communities of the Inner Mongolia Plateau[J].ActaEcologicaSinica,2002,22(8):1 215-1 223.

    [9]馬風(fēng)云,李新榮,張景光,等.沙坡頭人工固沙植被土壤水分空間異質(zhì)性[J].應(yīng)用生態(tài)學(xué)報(bào),2006,17(5):789-795.

    MA F Y,LI X R,ZHANG J G,etal.Spatial heterogeneity of soil moisture in Shapotou sand-fixing artificial vegetation area[J].ChineseJournalofAppliedEcology,2006,17(5):789-795.

    [10]李元壽,王根緒,丁永建,等.青藏高原高寒草甸區(qū)土壤水分的空間異質(zhì)性[J].水科學(xué)進(jìn)展,2008,19(1):61-67.

    LI Y S,WANG G X,DING Y J,etal.Spatial heterogeneity of soil moisture in alpine meadow area of the Qinghai Xizang Plateau[J].AdvancesinWaterScience,2008,19(1):61-67.

    [11]何志斌,趙文智.荒漠綠洲區(qū)人工梭梭林土壤水分空間異質(zhì)性的定量研究[J].冰川凍土,2004,26(2):207-211.

    HE Z B,ZHAO W Z.The spatial heterogeneity of soil moisture in artificialHaloxylonammodendron[J].JournalofGlaciologyandGeocryology,2004,26(2):207-211.

    一項(xiàng)建筑工程項(xiàng)目的完成需要涉及到多個(gè)行業(yè)和部門,還需要各個(gè)部門相互之間的溝通交流、共同協(xié)商來調(diào)整因?yàn)樘鞖庖蛩?、季?jié)因素、環(huán)境因素、現(xiàn)場條件而帶來的各個(gè)不良影響。在安全標(biāo)準(zhǔn)化管理的過程中,很容易會(huì)發(fā)生因?yàn)槿藶椴豢煽刂频囊蛩囟o工程施工帶來混亂的現(xiàn)象,增添了不可控制的風(fēng)險(xiǎn)。對此,施工人員可以提前根據(jù)安全標(biāo)準(zhǔn)管理措施做出對應(yīng)的應(yīng)急對策,方便施工人員能夠有效針對突發(fā)的各種情況并及時(shí)處理,降低建筑施工過程中的風(fēng)險(xiǎn)。

    [12]王紅梅,王仲良,王堃,等.華北農(nóng)牧交錯(cuò)帶農(nóng)田-草地景觀鑲嵌體土壤水分空間異質(zhì)性[J].生態(tài)學(xué)報(bào),2013,33(19):6 287-6 294.

    WANG H M,WANG Z L,WANG K,etal.Spatial heterogeneity of soil moisture across a cropland-grassland mosaic:a case study for agro-pastural transition in north of China[J].ActaEcologicaSinica,2013,33(19):6 287-6 294.

    [13]馬全林,魚泳,陳芳,等.干旱區(qū)沙漠化逆轉(zhuǎn)過程土壤水分的空間異質(zhì)性特征[J].干旱區(qū)地理,2010,33(5):716-724.

    MA Q L,YU Y,CHEN F,etal.Spatial heterogeneity of soil water content in the reversion process of desertification in arid area[J].AridLandGeography,2010,33(5):716-724.

    [14]李從娟,雷加強(qiáng),徐新文,等.古爾班通古特沙漠土壤水分與化學(xué)性質(zhì)的空間分布[J].生態(tài)學(xué)報(bào),2014,34(15):4 380-4 389.

    LI C J,LEI J Q,XU X W,etal.Spatial pattern for soil water and chemical properties in Gurbantunggut Desert[J].ActaEcologicaSinica,2014,34(15):4 380-4 389.

    [15]宋同清,彭晚霞,曾馥平,等.喀斯特木論自然保護(hù)區(qū)旱季土壤水分的空間異質(zhì)性[J].應(yīng)用生態(tài)學(xué)報(bào),2009,20(1):98-104.

    SONG T Q,PENG W X,ZENG F P,etal.Spatial heterogeneity of surface soil moisture content in dry season in Mulun National Natural Reserve in Karst area[J].ChineseJournalofAppliedEcology,2009,20(1):98-104.

    [16]陳文業(yè),張繼強(qiáng),趙明,等.甘肅省敦煌西湖荒漠-濕地植物群落物種多樣性特征研究[J].中國沙漠,2012,32(6):1 639-1 646.

    CHEN W Y,ZHANG J Q,ZHAO M,etal.Species diversity characteristics of plant community in Xihu Desert Wetland of Dunhuang,Gansu Province[J].JournalofDesertResearch,2012,32(6):1 639-1 646.

    [17]盧琦,李新榮,肖洪浪,等.荒漠生態(tài)系統(tǒng)觀測方法[M].北京:中國環(huán)境科學(xué)出版社,2004:12-20.

    [18]潘顏霞,王新平,蘇延桂,等.荒漠人工固沙植被區(qū)土壤水分的時(shí)空變異性[J].生態(tài)學(xué)報(bào),2009,29(2):993-1 000.

    PAN Y X,WANG X P,SU Y G,etal.Temporal and spatial variability of surface soil moisture in a re-vegetation desert area in Shapotou[J].ActaEcologicaSinica,2006,29(2):993-1 000.

    [19]邵方麗.冀北山地典型森林植被與土壤成分的空間異質(zhì)性關(guān)系研究[D].北京:北京林業(yè)大學(xué),2012:13-20.

    [20]顏亮,周廣勝,張峰,等.內(nèi)蒙古荒漠草原植被蓋度的空間異質(zhì)性動(dòng)態(tài)分析[J].生態(tài)學(xué)報(bào),2012,32(13):4 017-4 024.

    YAN L,ZHOU G S,ZHANG F,etal.Spatial heterogeneity of vegetation coverage and its temporal dynamics in desert steppe,Inner Mongolia[J].ActaEcologicaSinica,2012,32(13):4 017-4 024.

    [21]王蕙,趙文智,常學(xué)向.黑河中游荒漠綠洲過渡帶土壤水分與植被空間變異[J].生態(tài)學(xué)報(bào),2007,27(5):1 731-1 739.

    WANG H,ZHAO W Z,CHANG X X.Spatial variability of soil moisture and vegetation in desert-oasis ecotone in the middle reaches of Heihe River Basin[J].ActaEcologicaSinica,2007,27(5):1 731-1 739.

    [22]連綱,郭旭東,傅伯杰,等.黃土高原小流域土壤容重及水分空間變異特征[J].生態(tài)學(xué)報(bào),2006,26(3):647-654.

    LIAN G,GUO X D,FU B J,etal.Spatial variability of bulk density and soil water in a small catchment of the Loess Plateau[J].ActaEcologicaSinica,2006,26(3):647-654.

    [23]王政權(quán).地統(tǒng)計(jì)學(xué)及其在生態(tài)學(xué)中的應(yīng)[M].北京:科學(xué)出版社,1999:18-32.

    [24]蘇松錦.格氏栲天然林土壤空間異質(zhì)性及其生長適宜性評價(jià)[D].福州:福建農(nóng)林大學(xué),2012:20-31.

    [25]CAMBARDELLA C A,MOOMAN A T,NOVAK J M,etal.Field-scale variability of soil properties in central Iowa soil[J].SoilSci.Soc.Am.J.,1994,58:1 501-1 511.

    [26]王政權(quán),王慶成,李哈濱.紅松老齡林主要樹種的空間異質(zhì)性特征與比較的定量研究[J].植物生態(tài)學(xué)報(bào),2000,24(6):718-723.

    WANG Z Q,WANG Q C,LI H B.The spatial heterogeneity of soil physical properties in forests[J].ActaEcologicaSinica,2000,20(6):945-950.

    [27]陳鵬,初雨,顧峰雪,等.綠洲一荒漠過渡帶景觀的植被與土壤特征要素的空間異質(zhì)性分析[J].應(yīng)用生態(tài)學(xué)報(bào),2003,14(6):904-908.

    CHEN P,CHU Y,GU F G,etal.Spatial heterogeneity of vegetation and soil characteristics in oasis-desert ecotone[J].ChineseJournalofAppliedEcology,2003,14(6):904-908.

    [28]NYBERG L.Spatial variability of water content in the covered catchment at Gardsjon,Sweden[J].HydrologicalProcesses,1996,10:89-103.

    [29]KRUMBAEH A W.Effects of microrelief on distribution of soil moisture and bulk density[J].JournalofGeophysicalResearch,1959,64:1 587-1 590.

    [30]FAMIGLIETTI J S,RUDNICKI J W,RODELL M.Variability in surface moisture content along a hillslope transeet:Rattlesnake Hill Texa[J].JournalofHydrology,1998,210:259-281.

    (編輯:潘新社)

    Spatial Heterogeneity of Soil Moisture and Related Factors in Desert-Wetland Ecosystem Enclosed in Dunhuang Xihu,Gansu,China

    SUN Feida1,CHEN Wenye2,3,4*,YUAN Haifeng4,5,DOU Yingjie4,5,BING Danhui2,3,4,FENG Ying2,3,4,WU Ting4,5

    (1 Department of Grassland Science,Sichuan Agricultural University,Chengdu 611130,China;2 Gansu Forestry Science and Technology Research Academy,Lanzhou 730020,China;3 Gansu Forestry Science and Technology Engineering Company,Lanzhou 730020,China;4 Administrative Bureau of Dunhuang Xihu National Nature Reserve of Gansu,Dunhuang,Gansu 736200,China;5 National Positioning Observation Station of Xihu Wetland Ecosystem in Dunhuang of Gansu,Dunhuang,Gansu 736200,China)

    Abstract:Soil moisture is the dominant limiting factor of desert-wetland ecosystems in the inland.Based on the traditional and geographical statistics methods,we measured soil moisture ranged from 0 to 200 cm depth to find the rules of their spatial distribution characteristics,and surveyed some associated environmental factors such as altitude,soil texture and related vegetation indices in desert-wetland ecosystem at Xihu district,Dunhuang National Nature Reserve of Gansu Province,which may be helpful to vegetation restoration,protection and sustainable management in wetland ecosystem at extreme arid areas.The main results were as follows:(1)the coefficients of variation,Nugget variance,sill,range and spatial dependence were ranged from 36.51% to 88.65%,0.007 to 0.098,0.112 to 0.549,116 to 453 and 76.6% to 97.6%,respectively,and all their variables were at a medium variation level and a high heterogeneity and strong spatial autocorrelation was appeared.(2)The variability of soil moisture in the deep soil profile at 60 to 200 cm was higher than those in the shallow layer at 0 to 60 cm,which kept much more differences at different layers,especially in some small scales.(3)The altitude was a dominant factor that influenced the variation of deep soil moisture trend.(4)The relationships between soil texture and soil moisture in deep soil at 60 to 200 cm were stronger than those of shallow layers at 0 to 60 cm,otherwise,which did not comply with the factor of altitude;there were a great positive correlation between the herbage coverage and the soil moisture in shallow soil layers,but a negative correlation between the root biomass of shrub and the soil moisture in deep soil layers.

    Key words:Dunhuang Xihu;ecosystem desert and wetland;soilmoisture;spatial heterogeneity;influence factors

    中圖分類號:Q948.11

    文獻(xiàn)標(biāo)志碼:A

    作者簡介:孫飛達(dá)(1978-),男,博士,副教授,主要從事草地資源監(jiān)測與管理研究。E-mail:sunfd08@163.com*通信作者:陳文業(yè),學(xué)士,副研究員,主要從事生態(tài)恢復(fù)研究。E-mail:gschwy@163.com

    基金項(xiàng)目:甘肅省科技重大專項(xiàng)計(jì)劃(1302FKDA035);甘肅省科技支撐計(jì)劃-社會(huì)發(fā)展類項(xiàng)目(1011FKCA136)

    收稿日期:2015-10-27;修改稿收到日期:2016-01-10

    文章編號:1000-4025(2016)01-0165-09

    doi:10.7606/j.issn.1000-4025.2016.01.0165

    猜你喜歡
    荒漠土壤水分
    向荒漠宣戰(zhàn)
    磷素添加對土壤水分一維垂直入滲特性的影響
    中國石化(2021年8期)2021-11-05 07:00:14
    荒漠生態(tài)系統(tǒng)的演變
    綠色中國(2019年14期)2019-11-26 07:11:34
    三代人薪火相傳 誓把荒漠變綠洲
    吉林玉米生長期土壤水分規(guī)律分析
    荒漠水塘
    荒漠之路上的三首情歌
    文史春秋(2016年3期)2016-12-01 05:42:20
    西藏高原土壤水分遙感監(jiān)測方法研究
    不同覆蓋措施對棗園土壤水分和溫度的影響
    丁香六月欧美| 亚洲aⅴ乱码一区二区在线播放 | 两性午夜刺激爽爽歪歪视频在线观看 | 一级毛片女人18水好多| 黄频高清免费视频| 又紧又爽又黄一区二区| 日韩有码中文字幕| 欧美性猛交黑人性爽| 嫩草影视91久久| 久9热在线精品视频| 男人舔奶头视频| 亚洲国产日韩欧美精品在线观看 | 国产精品久久电影中文字幕| 久久亚洲真实| 久久人人精品亚洲av| 亚洲av成人一区二区三| 国产区一区二久久| 在线看三级毛片| 久久午夜亚洲精品久久| 熟女少妇亚洲综合色aaa.| 黄色片一级片一级黄色片| 欧美日韩福利视频一区二区| 午夜激情福利司机影院| 禁无遮挡网站| 日韩高清综合在线| 国产黄a三级三级三级人| 午夜免费激情av| 亚洲激情在线av| 妹子高潮喷水视频| 制服人妻中文乱码| 啪啪无遮挡十八禁网站| 琪琪午夜伦伦电影理论片6080| 色综合亚洲欧美另类图片| 亚洲专区中文字幕在线| 国产成+人综合+亚洲专区| aaaaa片日本免费| 欧美日韩一级在线毛片| 99国产精品一区二区蜜桃av| 成在线人永久免费视频| 亚洲av成人av| 亚洲国产欧洲综合997久久,| 久久久久亚洲av毛片大全| 久热爱精品视频在线9| 久久久久免费精品人妻一区二区| 国内精品久久久久精免费| 欧美成人免费av一区二区三区| 国产精品亚洲av一区麻豆| 久久久国产成人免费| 给我免费播放毛片高清在线观看| 黄频高清免费视频| 国产精品爽爽va在线观看网站| 1024香蕉在线观看| 久久久国产欧美日韩av| 欧美+亚洲+日韩+国产| 妹子高潮喷水视频| 亚洲专区国产一区二区| 窝窝影院91人妻| 中文字幕最新亚洲高清| 国产精品久久电影中文字幕| 99精品久久久久人妻精品| 搡老妇女老女人老熟妇| 国产一区二区激情短视频| 香蕉av资源在线| 久久精品aⅴ一区二区三区四区| 日韩欧美精品v在线| 最好的美女福利视频网| 不卡一级毛片| 999精品在线视频| 国产一区二区三区视频了| 一边摸一边做爽爽视频免费| 少妇熟女aⅴ在线视频| 国产亚洲精品一区二区www| 天堂影院成人在线观看| 可以在线观看的亚洲视频| 亚洲一码二码三码区别大吗| 亚洲精品粉嫩美女一区| 亚洲男人天堂网一区| 色在线成人网| 男人舔女人的私密视频| 全区人妻精品视频| 婷婷亚洲欧美| 激情在线观看视频在线高清| 国产在线精品亚洲第一网站| 亚洲男人的天堂狠狠| 免费在线观看成人毛片| 亚洲人成电影免费在线| 99国产精品一区二区蜜桃av| 他把我摸到了高潮在线观看| 欧美日韩瑟瑟在线播放| 久久草成人影院| 日日干狠狠操夜夜爽| 亚洲人成77777在线视频| 黄片小视频在线播放| 久久久久亚洲av毛片大全| 91国产中文字幕| 两个人视频免费观看高清| 男女床上黄色一级片免费看| 色在线成人网| 一边摸一边做爽爽视频免费| 亚洲专区国产一区二区| 女同久久另类99精品国产91| 国产精品一区二区三区四区免费观看 | 欧美黄色片欧美黄色片| avwww免费| 五月伊人婷婷丁香| 亚洲av成人不卡在线观看播放网| 国产精品亚洲一级av第二区| 欧美不卡视频在线免费观看 | 久久精品夜夜夜夜夜久久蜜豆 | 午夜a级毛片| 精品一区二区三区视频在线观看免费| 日韩精品中文字幕看吧| 制服人妻中文乱码| 最新在线观看一区二区三区| 国产亚洲av嫩草精品影院| 69av精品久久久久久| 久久精品91无色码中文字幕| 一夜夜www| 少妇裸体淫交视频免费看高清 | 亚洲aⅴ乱码一区二区在线播放 | 91老司机精品| 午夜成年电影在线免费观看| 老司机午夜福利在线观看视频| 国产欧美日韩一区二区精品| 亚洲精华国产精华精| 大型av网站在线播放| 两个人视频免费观看高清| 日韩三级视频一区二区三区| 亚洲美女视频黄频| 亚洲成人精品中文字幕电影| 亚洲成人精品中文字幕电影| 神马国产精品三级电影在线观看 | 男女视频在线观看网站免费 | 啦啦啦观看免费观看视频高清| 成年免费大片在线观看| 麻豆成人午夜福利视频| 亚洲精品久久国产高清桃花| 啪啪无遮挡十八禁网站| 性色av乱码一区二区三区2| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av美国av| 天天躁狠狠躁夜夜躁狠狠躁| 少妇人妻一区二区三区视频| 亚洲国产高清在线一区二区三| 欧美绝顶高潮抽搐喷水| 一级毛片精品| 中文字幕最新亚洲高清| 午夜成年电影在线免费观看| 欧美色欧美亚洲另类二区| 男人的好看免费观看在线视频 | 色尼玛亚洲综合影院| 麻豆国产av国片精品| 久久九九热精品免费| 国产精品久久久久久人妻精品电影| 亚洲七黄色美女视频| 亚洲av日韩精品久久久久久密| 精品乱码久久久久久99久播| 精华霜和精华液先用哪个| 99久久久亚洲精品蜜臀av| 久久久久久免费高清国产稀缺| 在线观看66精品国产| 天天一区二区日本电影三级| 亚洲午夜精品一区,二区,三区| 不卡av一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲欧美日韩东京热| 婷婷精品国产亚洲av| 天天一区二区日本电影三级| 国产成+人综合+亚洲专区| 黑人操中国人逼视频| 亚洲中文日韩欧美视频| 麻豆一二三区av精品| 特级一级黄色大片| 小说图片视频综合网站| 又黄又爽又免费观看的视频| 99国产精品一区二区三区| 久久久久久久久中文| 亚洲av片天天在线观看| 岛国在线观看网站| 国产乱人伦免费视频| 亚洲专区字幕在线| 又黄又粗又硬又大视频| 天天躁夜夜躁狠狠躁躁| 啦啦啦观看免费观看视频高清| 亚洲av电影不卡..在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲精品久久国产高清桃花| 香蕉丝袜av| 久久久久久国产a免费观看| 久久精品国产清高在天天线| 无遮挡黄片免费观看| 国产不卡一卡二| 亚洲五月天丁香| 亚洲人与动物交配视频| 夜夜夜夜夜久久久久| 精品一区二区三区四区五区乱码| 亚洲av熟女| 岛国视频午夜一区免费看| 久久中文看片网| 日韩欧美三级三区| 国产精品免费一区二区三区在线| 日韩 欧美 亚洲 中文字幕| 亚洲真实伦在线观看| 久久精品91蜜桃| 久久热在线av| 色综合婷婷激情| 精品久久久久久久人妻蜜臀av| 成人18禁在线播放| 免费在线观看亚洲国产| 制服诱惑二区| 一个人免费在线观看电影 | e午夜精品久久久久久久| 老鸭窝网址在线观看| 丁香欧美五月| 免费看a级黄色片| 老熟妇仑乱视频hdxx| 亚洲专区中文字幕在线| 熟女电影av网| 亚洲国产日韩欧美精品在线观看 | 啪啪无遮挡十八禁网站| 亚洲精品粉嫩美女一区| 无遮挡黄片免费观看| 欧美+亚洲+日韩+国产| 久热爱精品视频在线9| 亚洲av成人一区二区三| 久久热在线av| 国产日本99.免费观看| 久久精品成人免费网站| 我要搜黄色片| 九色成人免费人妻av| 91成年电影在线观看| 精品人妻1区二区| 日韩免费av在线播放| 亚洲精品在线观看二区| 精品人妻1区二区| 亚洲中文字幕一区二区三区有码在线看 | 日韩 欧美 亚洲 中文字幕| 男人舔奶头视频| 欧美日本亚洲视频在线播放| 亚洲国产看品久久| 亚洲av成人av| 欧美午夜高清在线| 日韩成人在线观看一区二区三区| 欧美黄色片欧美黄色片| 欧美一区二区国产精品久久精品 | 女生性感内裤真人,穿戴方法视频| 天天躁夜夜躁狠狠躁躁| 欧美不卡视频在线免费观看 | 欧美乱妇无乱码| 国产成人一区二区三区免费视频网站| 亚洲成av人片免费观看| 一本一本综合久久| 狠狠狠狠99中文字幕| 亚洲av日韩精品久久久久久密| 两个人的视频大全免费| ponron亚洲| 草草在线视频免费看| 亚洲av日韩精品久久久久久密| а√天堂www在线а√下载| 国产亚洲精品久久久久5区| 精品久久久久久久毛片微露脸| 日本 欧美在线| 久久香蕉激情| 精品日产1卡2卡| 欧美日韩乱码在线| 亚洲国产欧美一区二区综合| 亚洲第一欧美日韩一区二区三区| 人妻久久中文字幕网| 淫妇啪啪啪对白视频| xxx96com| 久久伊人香网站| 中文字幕久久专区| 国产亚洲精品av在线| 免费在线观看视频国产中文字幕亚洲| 午夜激情av网站| 午夜福利高清视频| 免费在线观看影片大全网站| 精品人妻1区二区| 午夜福利在线在线| 伊人久久大香线蕉亚洲五| 精品国产亚洲在线| 国产黄片美女视频| 亚洲欧美日韩东京热| 国产熟女xx| 国产亚洲精品av在线| 国产1区2区3区精品| 91老司机精品| 淫妇啪啪啪对白视频| 18禁裸乳无遮挡免费网站照片| 人成视频在线观看免费观看| 51午夜福利影视在线观看| 哪里可以看免费的av片| 久久精品影院6| 色综合亚洲欧美另类图片| 久久久国产欧美日韩av| 国产真人三级小视频在线观看| 男女做爰动态图高潮gif福利片| 国产三级黄色录像| 久久精品91无色码中文字幕| 美女 人体艺术 gogo| 中文字幕精品亚洲无线码一区| 久久精品国产99精品国产亚洲性色| 麻豆一二三区av精品| 亚洲国产高清在线一区二区三| 搡老岳熟女国产| 免费电影在线观看免费观看| 一个人免费在线观看的高清视频| 中国美女看黄片| 老司机午夜福利在线观看视频| 国产探花在线观看一区二区| 国产黄片美女视频| 亚洲熟女毛片儿| 亚洲精品在线观看二区| 色综合欧美亚洲国产小说| 又爽又黄无遮挡网站| 欧美久久黑人一区二区| 一级毛片高清免费大全| 欧美最黄视频在线播放免费| 无限看片的www在线观看| 成人特级黄色片久久久久久久| 亚洲国产精品999在线| 国产精品久久久久久精品电影| 美女高潮喷水抽搐中文字幕| 亚洲成人中文字幕在线播放| 精品一区二区三区视频在线观看免费| 丁香六月欧美| 亚洲在线自拍视频| 欧美成人一区二区免费高清观看 | 久99久视频精品免费| 亚洲第一欧美日韩一区二区三区| 怎么达到女性高潮| 白带黄色成豆腐渣| 国产精品98久久久久久宅男小说| 丁香六月欧美| 亚洲第一电影网av| 男女那种视频在线观看| 一区二区三区国产精品乱码| 亚洲熟女毛片儿| 日韩 欧美 亚洲 中文字幕| 国产真人三级小视频在线观看| 日韩免费av在线播放| 成人特级黄色片久久久久久久| 极品教师在线免费播放| 五月玫瑰六月丁香| 国产成人精品久久二区二区91| a级毛片a级免费在线| www.999成人在线观看| 久久国产乱子伦精品免费另类| 久久久久久九九精品二区国产 | 97人妻精品一区二区三区麻豆| 国产三级黄色录像| bbb黄色大片| 高清在线国产一区| 精品不卡国产一区二区三区| 久久天堂一区二区三区四区| 一边摸一边抽搐一进一小说| 熟女少妇亚洲综合色aaa.| 夜夜躁狠狠躁天天躁| 曰老女人黄片| 精品人妻1区二区| 欧美午夜高清在线| 亚洲欧美一区二区三区黑人| 亚洲五月婷婷丁香| 免费搜索国产男女视频| e午夜精品久久久久久久| 日日爽夜夜爽网站| 变态另类成人亚洲欧美熟女| 亚洲美女黄片视频| 国产黄色小视频在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲国产日韩欧美精品在线观看 | 夜夜躁狠狠躁天天躁| 不卡av一区二区三区| 日韩精品免费视频一区二区三区| 久久久久久久精品吃奶| 亚洲第一电影网av| 一二三四在线观看免费中文在| 啦啦啦免费观看视频1| 国产黄色小视频在线观看| 久久久久国产一级毛片高清牌| 国产精品亚洲美女久久久| 亚洲一区中文字幕在线| 日韩中文字幕欧美一区二区| 岛国视频午夜一区免费看| 大型av网站在线播放| 日韩高清综合在线| 国产高清videossex| 手机成人av网站| 18禁观看日本| 变态另类丝袜制服| 欧美人与性动交α欧美精品济南到| 在线看三级毛片| 91九色精品人成在线观看| 午夜视频精品福利| a级毛片a级免费在线| 少妇粗大呻吟视频| 久久久久久久久久黄片| 一夜夜www| 99热这里只有精品一区 | 男女之事视频高清在线观看| 91麻豆av在线| 欧美 亚洲 国产 日韩一| 久久人人精品亚洲av| 国产亚洲精品第一综合不卡| 国产成人系列免费观看| 午夜福利欧美成人| 欧美日本视频| 男女午夜视频在线观看| 母亲3免费完整高清在线观看| 身体一侧抽搐| 日韩欧美 国产精品| 99精品久久久久人妻精品| 日韩欧美在线乱码| 久久香蕉国产精品| 制服人妻中文乱码| 国内揄拍国产精品人妻在线| 欧美人与性动交α欧美精品济南到| 又黄又粗又硬又大视频| 免费高清视频大片| 麻豆成人午夜福利视频| 免费一级毛片在线播放高清视频| 精品久久久久久,| 在线观看午夜福利视频| 白带黄色成豆腐渣| 亚洲av美国av| 欧美大码av| 人妻夜夜爽99麻豆av| 久久久久久国产a免费观看| 欧美极品一区二区三区四区| 国产亚洲精品久久久久久毛片| 成在线人永久免费视频| 欧美国产日韩亚洲一区| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av成人不卡在线观看播放网| 国产真人三级小视频在线观看| 亚洲av日韩精品久久久久久密| 欧美日韩福利视频一区二区| 淫妇啪啪啪对白视频| 亚洲欧美精品综合一区二区三区| 天天躁夜夜躁狠狠躁躁| www.999成人在线观看| 啦啦啦韩国在线观看视频| svipshipincom国产片| 亚洲va日本ⅴa欧美va伊人久久| 久久人人精品亚洲av| 亚洲中文av在线| 欧美色视频一区免费| 无人区码免费观看不卡| 国产av麻豆久久久久久久| 悠悠久久av| 免费高清视频大片| 欧美日韩福利视频一区二区| 亚洲人成网站高清观看| 老司机午夜福利在线观看视频| 99国产精品一区二区三区| 午夜福利高清视频| 国产精品一区二区三区四区免费观看 | 亚洲av第一区精品v没综合| 丝袜人妻中文字幕| 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 国产成人aa在线观看| 看黄色毛片网站| 国产野战对白在线观看| 变态另类丝袜制服| 精品午夜福利视频在线观看一区| 免费在线观看视频国产中文字幕亚洲| 婷婷六月久久综合丁香| 国产精品免费一区二区三区在线| 男女视频在线观看网站免费 | 欧美一级a爱片免费观看看 | 不卡一级毛片| 午夜福利欧美成人| 国产成人精品久久二区二区91| 1024手机看黄色片| 天天一区二区日本电影三级| 极品教师在线免费播放| 狂野欧美白嫩少妇大欣赏| 亚洲成av人片在线播放无| 又粗又爽又猛毛片免费看| 国内毛片毛片毛片毛片毛片| 免费看日本二区| 免费无遮挡裸体视频| 黄片小视频在线播放| 国产99久久九九免费精品| 亚洲成人精品中文字幕电影| 欧美黑人巨大hd| x7x7x7水蜜桃| 亚洲一区高清亚洲精品| 午夜福利成人在线免费观看| 国产黄片美女视频| 久久中文看片网| av超薄肉色丝袜交足视频| 欧美日韩亚洲国产一区二区在线观看| 国产激情欧美一区二区| 亚洲 国产 在线| 久久中文字幕人妻熟女| 欧美成人免费av一区二区三区| 日韩精品中文字幕看吧| 久久亚洲精品不卡| 一本一本综合久久| 亚洲成av人片在线播放无| 国产av在哪里看| 日本a在线网址| 中文在线观看免费www的网站 | 女同久久另类99精品国产91| 五月玫瑰六月丁香| 免费在线观看亚洲国产| 国产免费男女视频| 1024香蕉在线观看| 三级毛片av免费| 不卡av一区二区三区| 丰满的人妻完整版| 别揉我奶头~嗯~啊~动态视频| 国产1区2区3区精品| ponron亚洲| 国产亚洲精品久久久久5区| 亚洲avbb在线观看| 麻豆久久精品国产亚洲av| av有码第一页| 国产熟女午夜一区二区三区| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 日韩av在线大香蕉| 99热这里只有是精品50| 人人妻,人人澡人人爽秒播| 日韩高清综合在线| 淫妇啪啪啪对白视频| 好男人电影高清在线观看| 久久中文字幕一级| 久久久精品大字幕| 欧美黑人精品巨大| 男人的好看免费观看在线视频 | 国产av在哪里看| 久久精品91无色码中文字幕| 午夜精品一区二区三区免费看| 午夜福利高清视频| а√天堂www在线а√下载| 亚洲欧洲精品一区二区精品久久久| 小说图片视频综合网站| 极品教师在线免费播放| √禁漫天堂资源中文www| 老熟妇仑乱视频hdxx| 搡老妇女老女人老熟妇| 最近最新中文字幕大全电影3| 国产爱豆传媒在线观看 | 黄频高清免费视频| 亚洲精品av麻豆狂野| 岛国视频午夜一区免费看| 最新美女视频免费是黄的| 美女黄网站色视频| 日韩大尺度精品在线看网址| АⅤ资源中文在线天堂| 99久久99久久久精品蜜桃| 淫秽高清视频在线观看| 久久国产精品人妻蜜桃| 黄片小视频在线播放| 国产成人av激情在线播放| 十八禁网站免费在线| 中文字幕久久专区| 麻豆国产97在线/欧美 | 久久久久免费精品人妻一区二区| 久久久国产精品麻豆| 99国产精品99久久久久| 哪里可以看免费的av片| 亚洲在线自拍视频| 欧美午夜高清在线| 亚洲天堂国产精品一区在线| 日韩 欧美 亚洲 中文字幕| 少妇的丰满在线观看| 国产精品久久久久久久电影 | 久久中文字幕人妻熟女| 最近最新中文字幕大全电影3| 国产免费av片在线观看野外av| 国产日本99.免费观看| 在线a可以看的网站| 亚洲18禁久久av| 国产成人欧美在线观看| 妹子高潮喷水视频| 欧美日韩乱码在线| 亚洲天堂国产精品一区在线| 午夜影院日韩av| 嫁个100分男人电影在线观看| 99国产精品一区二区三区| av有码第一页| 窝窝影院91人妻| 男女做爰动态图高潮gif福利片| 精品无人区乱码1区二区| 国内精品一区二区在线观看| 亚洲精品美女久久av网站| 日本撒尿小便嘘嘘汇集6| 国产精品影院久久| 午夜福利视频1000在线观看| 香蕉久久夜色| 亚洲精品国产精品久久久不卡| 亚洲国产精品久久男人天堂| 日韩精品免费视频一区二区三区| 欧美人与性动交α欧美精品济南到| 国产乱人伦免费视频| 听说在线观看完整版免费高清| 日本一本二区三区精品| 国内揄拍国产精品人妻在线| 夜夜夜夜夜久久久久| 亚洲全国av大片| 成人av在线播放网站| 在线a可以看的网站| 久久久久久久精品吃奶| 国产激情欧美一区二区| 亚洲 国产 在线| bbb黄色大片| 亚洲午夜精品一区,二区,三区| 日韩欧美一区二区三区在线观看| 国产成人av教育| 成人国产综合亚洲| 亚洲av日韩精品久久久久久密|