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

    石羊河流域自然植被對(duì)生態(tài)系統(tǒng)服務(wù)的約束效應(yīng)

    2020-06-11 11:15:58王玉純金貴峻王建邦胡秀芳
    生態(tài)學(xué)報(bào) 2020年9期
    關(guān)鍵詞:石羊河產(chǎn)水量覆蓋度

    高 超,趙 軍,王玉純,金貴峻,王建邦,胡秀芳

    1 西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,蘭州 735300 2 甘肅林業(yè)職業(yè)技術(shù)學(xué)院,天水 741020 3 南通大學(xué)地理科學(xué)學(xué)院,南通 226007

    生態(tài)系統(tǒng)服務(wù)是指生態(tài)系統(tǒng)所形成和所維持的人類賴以生存的自然條件及效用[1],是人類直接或間接從生態(tài)系統(tǒng)中所得到的收益[2]。生態(tài)系統(tǒng)服務(wù)功能的可持續(xù)性是人類社會(huì)可持續(xù)發(fā)展的重要基礎(chǔ)[2- 3],人類為滿足不斷增長(zhǎng)的服務(wù)需求和對(duì)服務(wù)的選擇性使用,常通過(guò)改造生態(tài)系統(tǒng)來(lái)增加供給服務(wù),造成其他類型服務(wù)不斷減少,調(diào)節(jié)服務(wù)尤為顯著。生態(tài)系統(tǒng)服務(wù)功能退化不僅危及當(dāng)代社會(huì)福祉,且將極大影響人類后代從生態(tài)系統(tǒng)中謀求利益,使人類社會(huì)可持續(xù)發(fā)展面臨巨大挑戰(zhàn)[4- 5]。

    植被生態(tài)系統(tǒng)的涵養(yǎng)水源、土壤保持及固碳釋氧等功能對(duì)維持區(qū)域經(jīng)濟(jì)、生態(tài)平衡具有重要地理價(jià)值。植被是生態(tài)系統(tǒng)服務(wù)產(chǎn)生和傳遞的重要基礎(chǔ),氣候變化和人為擾動(dòng)的影響結(jié)果也直接體現(xiàn)在植被上,植被的生長(zhǎng)狀況勢(shì)必對(duì)區(qū)域內(nèi)生態(tài)系統(tǒng)服務(wù)供給產(chǎn)生重要影響[6]。精確認(rèn)識(shí)植被變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響,對(duì)于深化生態(tài)系統(tǒng)服務(wù)研究具有重要意義。

    生態(tài)系統(tǒng)服務(wù)仍是學(xué)術(shù)界熱點(diǎn)之一,梳理國(guó)內(nèi)外研究成果發(fā)現(xiàn):評(píng)估的主流模型仍為InVEST模型,研究尺度主要以流域尺度或特定的自然保護(hù)區(qū)、城市(縣區(qū))[7];研究的生態(tài)功能以水源涵養(yǎng)、土壤保持、碳固定、水質(zhì)凈化及生物多樣性為主;研究方向從最初的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估逐步深化[8- 9],向生態(tài)服務(wù)的時(shí)空格局及驅(qū)動(dòng)因素[10- 11]、權(quán)衡與協(xié)同[12- 13]、生態(tài)恢復(fù)與補(bǔ)償[14]、生態(tài)服務(wù)供需流[15]等多方向延伸。在生態(tài)系統(tǒng)服務(wù)間的關(guān)系中,由于供給服務(wù)和調(diào)節(jié)服務(wù)之間存在矛盾[16],有學(xué)者以情境模擬法來(lái)實(shí)現(xiàn)生態(tài)服務(wù)的優(yōu)化;由于生態(tài)系統(tǒng)的復(fù)雜性和生態(tài)服務(wù)形成機(jī)制的研究不足,已有的情境模擬研究均通過(guò)改變某一個(gè)參數(shù)、其余參數(shù)不變進(jìn)行討論,較難準(zhǔn)確模擬該情境下的生態(tài)系統(tǒng)服務(wù)[17]。因此,介于植被在生態(tài)服務(wù)中的特殊地位,加強(qiáng)植被變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響是十分必要的。

    石羊河流域生態(tài)問(wèn)題受廣大學(xué)者關(guān)注,有關(guān)生態(tài)環(huán)境及生態(tài)系統(tǒng)服務(wù)研究已有成果發(fā)表[17- 21],在植被-生態(tài)服務(wù)、生態(tài)服務(wù)-生態(tài)服務(wù)間關(guān)系的研究中,多以相關(guān)分析進(jìn)行探討,缺乏非線性關(guān)系研究。約束線法在解決復(fù)雜生態(tài)學(xué)問(wèn)題中具有巨大潛力,為研究生態(tài)系統(tǒng)中的非線性關(guān)系提供了可能。因此,本文引入約束線法來(lái)研究石羊河流域植被覆蓋度對(duì)其產(chǎn)水量、土壤保持和固碳服務(wù)的約束關(guān)系,為生態(tài)系統(tǒng)服務(wù)變化的驅(qū)動(dòng)因素定量化研究提供參考。

    1 數(shù)據(jù)來(lái)源與評(píng)估方法

    1.1 研究區(qū)概況

    石羊河流域介于101°41′—104°16′E、36°29′—39°27′N之間,總面積4.16×104km2,平均海拔約2000 m,屬于典型的大陸溫帶干旱氣候,其上中下游的氣候特征具有明顯垂直分異:上游祁連山高寒半干旱濕潤(rùn)區(qū);中游走廊平原干旱區(qū),干旱程度明顯大于上游區(qū)域;下游溫暖干旱區(qū),主要為民勤綠洲和荒漠,屬于重度缺水區(qū)[22]。全流域共有8條子流域,按水文地質(zhì)單元又可分為3個(gè)獨(dú)立子水系[23]。

    1.2 數(shù)據(jù)來(lái)源

    植被覆蓋度(fractional vegetation cover,FVC)以MODIS-NDVI產(chǎn)品數(shù)據(jù)計(jì)算獲得,歸一化植被指數(shù)(normalized difference vegetation index,NDVI)產(chǎn)品空間分辨率為250 m,時(shí)間分辨率為16 d,以最大值合成法得到年數(shù)據(jù)。固碳服務(wù)以MODIS-NPP產(chǎn)品數(shù)據(jù)獲得,NPP產(chǎn)品空間分辨率為1000 m,時(shí)間分辨率為1 a。以上MODIS數(shù)據(jù)均來(lái)自https://modis.gsfc.nasa.gov/。產(chǎn)水服務(wù)與土壤服務(wù)采用InVEST模型計(jì)算,所使用的基礎(chǔ)數(shù)據(jù)見(jiàn)表1所列,對(duì)獲得的原始數(shù)據(jù)進(jìn)行處理和計(jì)算得到模型的輸入?yún)?shù)。

    表1 InVEST模型參數(shù)來(lái)源說(shuō)明

    1.3 評(píng)估方法

    本研究以上述數(shù)據(jù),定量估算2000、2005、2010、2015年石羊河流域自然植被覆蓋度、產(chǎn)水量、土壤保持及固碳服務(wù),提取植被覆蓋度與生態(tài)系統(tǒng)服務(wù)的約束線,在此基礎(chǔ)上總結(jié)植被覆蓋度對(duì)生態(tài)系統(tǒng)服務(wù)的約束效應(yīng),定量認(rèn)識(shí)植被覆蓋變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響,探索在約束效應(yīng)下石羊河流域生態(tài)系統(tǒng)服務(wù)優(yōu)化的理論。

    本研究使用多尺度數(shù)據(jù),其中NDVI數(shù)據(jù)重采樣為1 km;InVEST模型參數(shù)除DEM外均為1 km空間分辨率,經(jīng)驗(yàn)參數(shù)以相近研究區(qū)、模型可借鑒性及文獻(xiàn)的新近性為原則,以引文分析工具分析文獻(xiàn)獲得。Zhang系數(shù)、匯水累計(jì)閾值及坡度閾值需多次模擬,使其與研究區(qū)實(shí)際情況相符。

    1.3.1植被覆蓋度估算

    植被覆蓋度采用像元二分模型,以NDVI數(shù)據(jù)計(jì)算植被覆蓋度,公式為:

    FVC=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

    式中,FVC為植被覆蓋度,NDVIsoil和NDVIveg分別為純裸土和純植被對(duì)應(yīng)的NDVI值。結(jié)合前人研究和研究區(qū)實(shí)際情況,選取生長(zhǎng)季內(nèi)植被NDVI最小值0.076和最大值0.8179代替NDVIsoil和NDVIveg[32]。

    1.3.2生態(tài)系統(tǒng)服務(wù)評(píng)估

    (1)產(chǎn)水量。以水量平衡公式計(jì)算,公式如下:

    式中,AET(x)表示柵格單元x的年實(shí)際蒸散量,P(x)表示柵格單元x的年降水量。

    (2)土壤保持。土壤保持服務(wù)以土壤侵蝕和土壤保持量為指標(biāo),土壤保持量定義為裸地條件下的土壤侵蝕量RKLS與植被覆蓋情況下的土壤侵蝕量USLE之差。公式如下:

    RKLS=R·K·LS

    USLE=R·K·LS·C·P

    SD=RKLS-USLE

    式中,RKLS為潛在土壤侵蝕,USLE為實(shí)際土壤侵蝕,SD為土壤保持量;R為降水可蝕量因子,K為土壤可蝕性因子,LS為坡度坡長(zhǎng)因子,C為植被覆蓋和管理因子,P為土壤保持措施因子。

    (3)固碳服務(wù)。

    6CO2+12H2O=C6H12O6+6O2+6H2O

    根據(jù)上式,植物干物質(zhì)生產(chǎn)與CO2消耗比例為1∶1.62,因此固碳量為:

    Vc=NPP·1.62

    1.3.3約束線法

    圖1 約束線示意圖Fig.1 Constraint line diagram

    約束線法是以分段分位數(shù)回歸刻畫的,約束線上的點(diǎn)代表研究?jī)蓚€(gè)生態(tài)變量時(shí),受到其他變量的影響最小。由于植被的生長(zhǎng)情況是限制其提供生態(tài)系統(tǒng)服務(wù)量的主要因素,因此將表征植被生長(zhǎng)狀況的植被覆蓋度作為坐標(biāo)橫軸,提供的生態(tài)系統(tǒng)服務(wù)作為坐標(biāo)縱軸,構(gòu)成二維坐標(biāo)系(圖1),其中限制變量為植被覆蓋度(X),響應(yīng)變量為生態(tài)系統(tǒng)服務(wù)(Y)。約束線上的點(diǎn)表示變量Y受變量X影響時(shí),受其他因素的影響最小,即變量Y的理論最大值。

    約束線法以O(shè)rigin 2017擬合,所有約束線均通過(guò)95%置信度檢驗(yàn);結(jié)合研究對(duì)象間的關(guān)系及前人研究[33],將約束線歸納為以下幾種類型(圖2)。

    圖2 約束線類型Fig.2 Constraint line type正線型:變量X對(duì)變量Y的約束效應(yīng)成比例下降;正凸型:變量X能夠有效的約束變量Y,且成比例增加約束效應(yīng);指數(shù)型:變量X對(duì)變量Y的約束效應(yīng)在整個(gè)范圍內(nèi)持續(xù)下降;S型:隨著變量X的增加,其對(duì)變量Y的約束效應(yīng)在X的某個(gè)范圍內(nèi)急劇下降;開(kāi)口向下拋物線型:隨著變量X的增加,其對(duì)變量Y的約束效應(yīng)先逐漸減少,后逐漸增大;凸波型:變量X對(duì)變量Y的約束效應(yīng)表現(xiàn)出波動(dòng)特征

    1.3.4生態(tài)系統(tǒng)服務(wù)潛力估算

    約束線上的點(diǎn)表示在其他因素的最小約束下,一個(gè)生態(tài)變量受到另一個(gè)生態(tài)變量約束,因此代表了效率前沿,類似于經(jīng)濟(jì)學(xué)中帕累托效率的含義[34-35]。在生態(tài)系統(tǒng)管理中,實(shí)際生態(tài)系統(tǒng)服務(wù)量與效率前沿的差值就是對(duì)應(yīng)生態(tài)系統(tǒng)服務(wù)的提升潛力,也是生態(tài)系統(tǒng)服務(wù)優(yōu)化的重點(diǎn)。本文以2015年石羊河流域生態(tài)系統(tǒng)服務(wù)為例,分別計(jì)算產(chǎn)水、土壤保持、固碳等三種服務(wù)潛力的空間分布,并對(duì)3種服務(wù)潛力進(jìn)行加權(quán)分析,確定綜合潛力分布態(tài)勢(shì)。生態(tài)系統(tǒng)服務(wù)潛力表示為:

    Pi=Li-Ai

    式中,Pi表示每個(gè)柵格單元的潛力值,Li表示每個(gè)柵格單元的理想值,Ai表示每個(gè)柵格單元的實(shí)際值。

    在干旱區(qū),水資源可用性是維持各項(xiàng)生態(tài)過(guò)程最重要的制約因素。在測(cè)算綜合潛力時(shí),產(chǎn)水量的權(quán)重應(yīng)高于土壤保持和固碳服務(wù)的權(quán)重。綜合潛力表示為:

    CPi=0.4×Wi+0.3×Si+0.3×Ni

    式中,CPi為每個(gè)柵格單元的綜合潛力值,Wi為每個(gè)柵格單元的產(chǎn)水量潛力值,Si為每個(gè)柵格單元的土壤保持潛力值,Ni為每個(gè)柵格單元的固碳服務(wù)潛力值。

    2 結(jié)果與分析

    2.1 植被覆蓋與生態(tài)系統(tǒng)服務(wù)的約束效應(yīng)

    根據(jù)USGS分類體系,本文從整體植被和分林草植被兩個(gè)層次提取約束線,并在此基礎(chǔ)上進(jìn)行討論,整體植被是將所有自然植被作為整體提取約束線,包括所有的林地和草地,分林草植被是將林地和草地分開(kāi),分別提取約束線。

    2.1.1植被覆蓋與產(chǎn)水量的約束效應(yīng)

    從整體來(lái)看,植被覆蓋度與產(chǎn)水量的約束線基本呈開(kāi)口向下拋物線類型,但不同植被類型的約束線略有差別(圖3)。在閾值左側(cè),隨著植被覆蓋度不斷增加,植被覆蓋度對(duì)產(chǎn)水量的約束效應(yīng)逐漸減小,其他因素對(duì)產(chǎn)水量的影響不斷增加。當(dāng)植被覆蓋度不斷增加時(shí),植被對(duì)降水的截留能力在不斷增加,降水量也在不斷增加,表明植被可利用水分在增加,但總體處于缺水狀況,在植被可用水分不能充分滿足植被生長(zhǎng)的情況下,更多的水分被植被吸收用于自身生長(zhǎng),因此這一階段水分成為限制植被提供產(chǎn)水服務(wù)的主要因素。在閾值右側(cè),隨著植被覆蓋度不斷增加,植被覆蓋度對(duì)產(chǎn)水量的約束效應(yīng)逐漸增大,其他因素對(duì)產(chǎn)水量的影響不斷減小。石羊河流域高植被覆蓋區(qū)位于流域上游,為典型的山地區(qū)。高植被覆蓋度下枝葉的凋落量大,分解速率不高,枯枝存量大,造成枯枝落葉層對(duì)于水分的大量截留,枯枝落葉層的截留水量達(dá)到飽和后部分水分將滲入土中,增加了土壤含水量。同時(shí),高植被覆蓋從一個(gè)側(cè)面也反映出該區(qū)有更大的降水量,充足的水分在滿足了植被自身的生長(zhǎng)需求后產(chǎn)生徑流,加之山區(qū)坡度較大處的土壤更易形成土壤侵蝕,影響植被生長(zhǎng),進(jìn)而影響到植被提供產(chǎn)水服務(wù)的能力。

    圖3 植被覆蓋度與產(chǎn)水量約束線Fig.3 Vegetation coverage and water production constraint line

    從時(shí)間尺度來(lái)看,同類植被的約束線類型基本保持不變,但表現(xiàn)出一定的時(shí)間分異特征。分析來(lái)看,造成時(shí)間分異的本質(zhì)原因是降水量、潛在蒸散量等其他影響因素的不同組合導(dǎo)致。其中,土壤深度、植物可利用含水率、植物蒸散系數(shù)等在一定時(shí)間內(nèi)基本保持穩(wěn)定,因此年降水量和潛在蒸散量的變化是造成約束線具有時(shí)間分異特征的主要原因。

    從擬合優(yōu)度R2來(lái)看,整體植被的約束線和草地的約束線擬合效果極佳,在不同時(shí)間尺度下擬合優(yōu)度較為接近。林地的約束線擬合效果總體較低。由于石羊河流域林地?cái)?shù)量較少,且主要分布在上游地區(qū),植被覆蓋度較高,低植被覆蓋度地區(qū)林地?cái)?shù)量較少,這可能是影響擬合結(jié)果的主要原因。

    從約束線的閾值來(lái)看,草地閾值明顯高于林地閾值,表明在理想條件下,同等植被覆蓋度時(shí)草地能提供更多的產(chǎn)水服務(wù)。結(jié)合林地與草地提供產(chǎn)水量的實(shí)際能力與理論最大值(圖3)發(fā)現(xiàn),在現(xiàn)實(shí)情況下,林地的產(chǎn)水能力強(qiáng)于草地;在理想情況下,草地的產(chǎn)水能力強(qiáng)于林地。

    2.1.2植被覆蓋與土壤保持的約束效應(yīng)

    植被覆蓋度與土壤保持的約束線均呈開(kāi)口向下拋物線類型,與產(chǎn)水量的約束線相似(圖4)。在閾值左側(cè),植被覆蓋度總體處于中、低覆蓋水平,主要分布在流域中下游,屬于干旱區(qū)。由于降水較少,蒸散量較大,總體處于缺水環(huán)境,加之植被覆蓋度不高,對(duì)于降水的保留能力不強(qiáng),大部分降水被直接蒸散散失。有研究表明,在干旱區(qū)小于13 mm的降雨屬于無(wú)效降雨,全部降水都將被蒸發(fā)消耗掉[36]。土壤保持主要受降水可蝕性因子、土壤可蝕性因子、坡度坡長(zhǎng)因子、土壤飽和導(dǎo)水率因子、土壤保持因子及植被覆蓋和管理因子影響。在降水較少的干旱區(qū),降水可蝕性因子和土壤飽和導(dǎo)水率可以忽略;中下游地勢(shì)較為平坦,因此坡度坡長(zhǎng)因子也可以不予考慮;植被覆蓋和管理因子和土壤保持因子對(duì)土壤保持具有促進(jìn)作用,因此土壤可蝕性因子成為限制植被提供土壤保持的主要因素。土壤可蝕性因子大小取決于土壤的特性指標(biāo),包括砂粒含量、土壤有機(jī)質(zhì)含量等,結(jié)合石羊河流域中下游自然地理特征,可以認(rèn)為風(fēng)蝕是該區(qū)土壤侵蝕的主要類型。另外,由于石羊河流域上游是產(chǎn)水量的高值區(qū),上游來(lái)水中攜帶的礦物質(zhì)隨徑流來(lái)到中下游,而中下游蒸發(fā)強(qiáng)烈,可能出現(xiàn)土壤鹽堿化,限制植被生長(zhǎng),進(jìn)而限制植被提供土壤保持服務(wù)的能力。在閾值右側(cè),植被覆蓋度總體處于高覆蓋水平,主要分布在流域上游,區(qū)域內(nèi)降水量較大,降水強(qiáng)度較高,地形起伏較大,這些因素不但會(huì)對(duì)土壤造成水蝕,同時(shí)形成的徑流會(huì)帶走土壤中的部分礦物質(zhì)和有機(jī)鹽等,影響土壤質(zhì)量。

    在時(shí)間尺度、擬合優(yōu)度及閾值上,土壤保持約束線的結(jié)果與產(chǎn)水量相似。

    圖4 植被覆蓋度與土壤保持的約束線Fig.4 Vegetation coverage and soil retention constraints

    2.1.3植被覆蓋與固碳服務(wù)的約束效應(yīng)

    植被覆蓋與固碳服務(wù)的約束線呈類直線的正凸型(圖5),隨著植被覆蓋度的不斷增加,植被提供的固碳服務(wù)量也不斷增加。由于正凸型函數(shù)類型不存在閾值,所以植被覆蓋度越高,提供的固碳服務(wù)量越高。

    固碳服務(wù)的約束線同樣存在時(shí)間分異特征。從整體植被的角度看,2000年的固碳量最大值在950左右,為最大值,2005年的固碳量最大值在750以下,其余兩年的對(duì)應(yīng)最大值在800—850之間。

    在擬合優(yōu)度上,整體植被與分林草植被的約束線的擬合效果均較為理想,R2值較為接近。

    對(duì)比不同時(shí)期林地和草地約束線發(fā)現(xiàn),同等植被覆蓋度情況下,林地的固碳量理論最大值高于草地的固碳量最大值;結(jié)合林地與草地的實(shí)際固碳能力,發(fā)現(xiàn)在實(shí)際與理想情況下,同等植被覆蓋度時(shí)林地的固碳能力均大于草地的固碳能力。

    圖5 植被覆蓋與固碳服務(wù)的約束線Fig.5 Constraint line for vegetation cover and carbon sequestration services

    2.1.4模型評(píng)價(jià)

    在約束線上,一個(gè)生態(tài)變量受到另一個(gè)生態(tài)變量的約束時(shí),受到其他影響因素的影響最小。以植被覆蓋對(duì)產(chǎn)水量的約束線為例,產(chǎn)水量受植被、氣象、土壤及地形等多因素影響,約束線上的點(diǎn)表示當(dāng)產(chǎn)水量受植被覆蓋度影響時(shí),受氣象、土壤、地形等其他因素的最小影響。所以約束線的意義在于,在研究某兩種生態(tài)變量的關(guān)系時(shí),能夠最大限度地減少局部因素的影響。因此,為盡可能的減少局部因素的影響,研究對(duì)象應(yīng)盡可能的廣泛分布,以涵蓋不同的氣候、土壤及地形等類型。

    應(yīng)用約束線分析生態(tài)問(wèn)題時(shí),在研究區(qū)分布廣泛的變量有較好的擬合效果,更適合約束線理論分析。由于石羊河流域林地面積較少且集中分布,主要集中在上游地區(qū),應(yīng)用約束線法減小局部因素的影響效果較差,因此擬合優(yōu)度比草地的擬合優(yōu)度差。

    2.2 植被與生態(tài)系統(tǒng)服務(wù)的響應(yīng)關(guān)系

    通過(guò)不同植被類型下植被覆蓋與生態(tài)系統(tǒng)服務(wù)的約束線,發(fā)現(xiàn)在現(xiàn)實(shí)、理想情況下草地和林地的服務(wù)供給能力差異,因此,有必要討論不同植被類型下植被與生態(tài)系統(tǒng)服務(wù)的響應(yīng)關(guān)系,以期為植被恢復(fù)類型選擇提供科學(xué)依據(jù)。從現(xiàn)實(shí)情況看,林地提供產(chǎn)水、土壤保持及固碳服務(wù)的能力強(qiáng)于草地,因此增加林地是提高服務(wù)供給的更佳選擇。從理想情況來(lái)看,草地的產(chǎn)水能力和土壤保持能力均大于草地,與實(shí)際情況相反。這是由于理想條件下降水量、潛在蒸散量及土壤可蝕性等限制服務(wù)供給的主要因素被最大限度削弱導(dǎo)致;結(jié)合石羊河流域具體情況及林地的分布情況,林地對(duì)生長(zhǎng)環(huán)境的要求相對(duì)更高。另外,荒漠生態(tài)系統(tǒng)服務(wù)功能監(jiān)測(cè)與評(píng)估技術(shù)研究項(xiàng)目組指出,植被冠層對(duì)降水量截留率在2.7%—33.1%間,最大可達(dá)50 mm,這部分降水多以蒸發(fā)形式直接散失[37];尹立河等[38]指出,植被覆蓋度增加會(huì)減少地下水補(bǔ)給,導(dǎo)致河流流量降低;李小英等[39]指出,相同地區(qū)人工植被耗水量大于天然植被耗水量。因此,選擇自然恢復(fù)或耗水量較小的當(dāng)?shù)刂脖豢蓪?shí)現(xiàn)干旱區(qū)水與生態(tài)的和諧發(fā)展。

    2.3 生態(tài)系統(tǒng)服務(wù)潛力估算

    2.3.1產(chǎn)水量服務(wù)潛力

    圖6 石羊河流域產(chǎn)水服務(wù)潛力分布圖 Fig.6 Distribution map of water production service potential in Shiyang River Basin

    根據(jù)1.3.4中提出的潛力估算模型計(jì)算得到石羊河流域產(chǎn)水量服務(wù)潛力分布圖(圖6)。由圖6可知,研究區(qū)產(chǎn)水量服務(wù)潛力總體呈現(xiàn)西南高東北低、具有流域分層的分布特征。流域上游是產(chǎn)水服務(wù)潛力值的高值區(qū),植被茂密,水熱條件較好,也是產(chǎn)水服務(wù)的高值區(qū)。流域中游西部山區(qū)的產(chǎn)水量服務(wù)潛力較高,流域中游中部和東部是產(chǎn)水量服務(wù)潛力的中低值區(qū)。流域下游是產(chǎn)水量服務(wù)潛力的低值區(qū),提升空間不大,僅有零星區(qū)域有一定提升潛力。

    綜合石羊河流域2015年實(shí)際產(chǎn)水量和產(chǎn)水量潛力分析發(fā)現(xiàn),上游產(chǎn)水量服務(wù)的高產(chǎn)值和高潛力并存。石羊河流域上游生態(tài)條件較好,是產(chǎn)水量服務(wù)的高值區(qū);流域中游是人類活動(dòng)的主要集中地,水資源使用和消耗量全流域最大,短缺也最為嚴(yán)重,中游地區(qū)人類活動(dòng)對(duì)水資源的超量利用也是造成下游水資源短缺的重要原因。正是因?yàn)橹邢掠螄?yán)重的生態(tài)問(wèn)題和水資源短缺,社會(huì)各界在管理和關(guān)注生態(tài)問(wèn)題時(shí)更側(cè)重于中下游地區(qū)。

    面對(duì)石羊河流域的水資源問(wèn)題,可持續(xù)發(fā)展和合理節(jié)約利用水資源是我們關(guān)注的重點(diǎn),但不是解決水資源問(wèn)題的唯一途徑。通過(guò)對(duì)石羊河流域產(chǎn)水服務(wù)潛力高值區(qū)的科學(xué)管理,也是解決石羊河流域水資源問(wèn)題的一種行之有效的策略。

    2.3.2土壤保持服務(wù)潛力

    由圖7可知,研究區(qū)土壤保持服務(wù)潛力總體呈現(xiàn)西南高東北低、流域分層明顯的分布特征。石羊河流域上游地形起伏度較大,降水較多,容易造成土壤流失,但上游茂密的植被在一定程度上減少了地形和降水造成的土壤流失,是土壤保持的高值區(qū),也是土壤保持服務(wù)潛力的高值區(qū)。流域中游主要為土壤保持服務(wù)潛力中低值分布區(qū),有一定的提升潛力。流域下游是土壤保持服務(wù)潛力的低值區(qū),提升空間不大??傮w來(lái)說(shuō),流域西南部及西部山區(qū)地帶土壤保持服務(wù)有較大提升空間,中游中部區(qū)域土壤保持服務(wù)也有一定提升空間。

    綜合石羊河流域2015年實(shí)際土壤保持量和土壤保持潛力分析發(fā)現(xiàn),石羊河流域上游土壤保持服務(wù)的高產(chǎn)值和高潛力并存;流域中下游土壤保持服務(wù)的低產(chǎn)值和低潛力并存。石羊河流域上游生態(tài)條件較好,但降水量、降水強(qiáng)度、地形等因素同樣限制土壤保持服務(wù)的供給,合理開(kāi)發(fā)流域上游的土壤保持服務(wù)潛力,是改善區(qū)域生態(tài)問(wèn)題的有效途徑。石羊河流域中下游生態(tài)問(wèn)題突出,提供土壤保持服務(wù)的能力較弱,提升潛力較小,改善中下游的生態(tài)問(wèn)題是當(dāng)下的難題。

    2.3.3固碳服務(wù)潛力

    由圖8可知,研究區(qū)固碳服務(wù)潛力總體呈現(xiàn)西南高東北低、流域分層分布特征。石羊河流域固碳服務(wù)高值區(qū)和服務(wù)潛力高值區(qū)分布在上游;流域固碳服務(wù)的中值區(qū)和服務(wù)潛力的中值區(qū)主要分布在中游;流域固碳服務(wù)的低值區(qū)和服務(wù)潛力的低值區(qū)主要分布在下游。

    石羊河流域的植被覆蓋度及植被分布范圍總體呈從上游向下游遞減趨勢(shì),結(jié)合上述分析推測(cè),增加植被覆蓋度及植被分布范圍可能提高固碳量。有研究指出,石羊河流域下游的生態(tài)輸水能提高植被覆蓋度,有利于改善區(qū)域生態(tài)環(huán)境[7]。因此,通過(guò)人工增加植被、生態(tài)輸水等方式,改善區(qū)域生態(tài)環(huán)境,是增加固碳服務(wù)的有效方式。

    圖7 石羊河流域土壤保持服務(wù)潛力分布圖 Fig.7 Distribution map of soil conservation service potential in Shiyang River Basin

    圖8 石羊河流域固碳服務(wù)潛力分布圖 Fig.8 Distribution map of carbon sequestration service potential in Shiyang River Basin

    2.3.4綜合潛力

    圖9 石羊河流域綜合潛力分布圖 Fig.9 Shiyang River Basin comprehensive potential distribution map

    根據(jù)產(chǎn)水、土壤保持、固碳等3種服務(wù)的潛力值柵格圖,利用1.3.4中的綜合潛力估算模型計(jì)算得到石羊河流域生態(tài)系統(tǒng)服務(wù)綜合潛力分布圖(圖9)。分析圖9發(fā)現(xiàn),生態(tài)系統(tǒng)服務(wù)綜合潛力高值區(qū)主要分布在上游西南部,表明石羊河流域上游具有較高的生態(tài)系統(tǒng)服務(wù)能力提升空間;中下游也有不同程度的提升空間。

    結(jié)合產(chǎn)水量服務(wù)、土壤保持服務(wù)及固碳服務(wù)空間分布圖發(fā)現(xiàn),3種服務(wù)潛力的主要提升空間均集中在流域上游,與綜合潛力空間分布相似。因此,有必要從子流域尺度進(jìn)行分析。

    由表2可知,在上游的8條子流域中,總潛力前4位的子流域按總潛力值大小依次為:東大河>西大河>西營(yíng)河>古浪河。結(jié)合產(chǎn)水量、土壤保持及固碳服務(wù)的子流域占比情況發(fā)現(xiàn),產(chǎn)水量較高的子流域?yàn)椋汗爬撕?東大河>西營(yíng)河>雜木河>西大河;土壤保持較高的子流域?yàn)椋何鳡I(yíng)河>東大河>雜木河>黃羊河>古浪河>西大河;固碳量較高的子流域?yàn)椋簴|大河>西大河>西營(yíng)河>古浪河。分析發(fā)現(xiàn),東大河、西大河、西營(yíng)河及古浪河提供的3種服務(wù)較多,且總潛力較大,說(shuō)明這4條子流域不僅是流域內(nèi)各項(xiàng)生態(tài)系統(tǒng)服務(wù)的主要供給者,也是流域內(nèi)生態(tài)系統(tǒng)服務(wù)潛力提升的主要貢獻(xiàn)者。

    表2 2015年各子流域自然植被綜合潛力占流域總值的比例/%

    分析來(lái)看,植被的分布狀況是造成子流域綜合潛力差異的原因之一。另外,綜合潛力由產(chǎn)水量潛力、土壤保持潛力及固碳潛力計(jì)算得出,因此,有必要從三項(xiàng)服務(wù)潛力分別分析。從產(chǎn)水量潛力來(lái)看,水熱條件組合是影響產(chǎn)水量潛力的主要原因,對(duì)應(yīng)到輸入?yún)?shù)為年降水量及潛在蒸散量。在較大尺度,通過(guò)人為影響的方式改善年降水量及潛在蒸散量,從而達(dá)到生態(tài)恢復(fù)是不現(xiàn)實(shí)的。從土壤保持潛力來(lái)看,土壤可蝕性是影響土壤保持潛力的主要原因,但不論是降水可蝕性,還是土壤可蝕性,這些因素都較難通過(guò)人工來(lái)大范圍改變。從固碳潛力來(lái)看,提高植被覆蓋度有利用增加固碳服務(wù),生態(tài)輸水能提高植被覆蓋度,因此,通過(guò)生態(tài)輸水,植樹(shù)植草等方式,可增加區(qū)域的固碳服務(wù)。結(jié)合植被與生態(tài)系統(tǒng)服務(wù)的響應(yīng)關(guān)系分析發(fā)現(xiàn),根據(jù)各子流域的實(shí)際情況,在自然恢復(fù)的基礎(chǔ)上,增加耗水量較少的當(dāng)?shù)刂脖荒茉谝欢ǔ潭壬细纳粕鷳B(tài)環(huán)境,提高服務(wù)供給。

    3 結(jié)論

    本文基于分段分位數(shù)回歸法提取了植被覆蓋與生態(tài)系統(tǒng)服務(wù)之間的約束線,探討了植被覆蓋對(duì)生態(tài)系統(tǒng)服務(wù)的約束效應(yīng),并以2015年數(shù)據(jù)為例,分析了石羊河流域生態(tài)系統(tǒng)服務(wù)潛力,得出以下結(jié)論:

    (1)植被覆蓋度與產(chǎn)水量、土壤保持服務(wù)的約束線均呈開(kāi)口向下拋物線類型,隨著植被覆蓋度的增加,兩種服務(wù)均呈先增加后減小的變化趨勢(shì),表明植被覆蓋度對(duì)兩種服務(wù)的約束效應(yīng)均呈先減小后增加的變化特征;植被覆蓋度與固碳服務(wù)間的約束線呈類直線的正凸型,表明隨著植被覆蓋度增加,其對(duì)固碳服務(wù)的約束效應(yīng)逐漸減小。

    (2)植被覆蓋對(duì)產(chǎn)水量、土壤保持及固碳服務(wù)的約束效應(yīng)均存在時(shí)空異質(zhì)性。在理想情況下,草地的產(chǎn)水能力與土壤保持能力均強(qiáng)于林地,與實(shí)際情況相反;草地的固碳能力弱于林地,與實(shí)際情況相同,但差距不大。

    (3)約束線法代表效率前沿,可用于分析生態(tài)系統(tǒng)服務(wù)能力的可提升潛力。石羊河流域產(chǎn)水、土壤保持、固碳3種服務(wù)的提升潛力均呈西南高東北低、高服務(wù)值和高潛力并存的分布特征,東大河、西大河、西營(yíng)河及古浪河流域不僅是流域內(nèi)各項(xiàng)生態(tài)系統(tǒng)服務(wù)的主要供給者,也是流域內(nèi)生態(tài)系統(tǒng)服務(wù)潛力提升的主要貢獻(xiàn)者。

    (4)約束線法在生態(tài)系統(tǒng)服務(wù)領(lǐng)域有較好的應(yīng)用前景,其存在時(shí)空異質(zhì)性的根本原因在于其他因素的時(shí)空異質(zhì)性。受理論和方法限制,具有在研究區(qū)廣泛分布特征的要素可達(dá)到較好的擬合效果,更適合約束線法的推廣應(yīng)用。

    (5)石羊河流域的生態(tài)系統(tǒng)服務(wù)優(yōu)化通過(guò)生態(tài)恢復(fù)進(jìn)行,應(yīng)結(jié)合子流域的實(shí)際情況,在自然恢復(fù)的基礎(chǔ)上,通過(guò)增加耗水量較少的當(dāng)?shù)刂脖桓纳粕鷳B(tài)環(huán)境,提高服務(wù)供給,實(shí)現(xiàn)干旱區(qū)生態(tài)與服務(wù)的和諧發(fā)展。

    猜你喜歡
    石羊河產(chǎn)水量覆蓋度
    人民黃河(2023年7期)2023-08-27 15:41:53
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測(cè)
    基于INVEST 模型的資水流域產(chǎn)水量及其對(duì)環(huán)境響應(yīng)的評(píng)估分析
    贛江流域產(chǎn)水功能對(duì)土地利用變化的響應(yīng)
    雅魯藏布江下游產(chǎn)水量時(shí)空演變及對(duì)氣候和土地利用變化的響應(yīng)
    基于不同旱情指數(shù)的石羊河流域春旱監(jiān)測(cè)研究
    基于NDVI的晉州市植被覆蓋信息提取
    石羊河流域永昌縣地下水及水資源供需平衡分析
    低覆蓋度CO分子在Ni(110)面的吸附研究
    煤層氣井長(zhǎng)沖程、大泵徑排采設(shè)備的研究及應(yīng)用
    日韩,欧美,国产一区二区三区| 久久久久精品国产欧美久久久 | 国产淫语在线视频| 久久人妻熟女aⅴ| 日韩有码中文字幕| 一区二区三区四区激情视频| 飞空精品影院首页| 18禁国产床啪视频网站| 日韩大码丰满熟妇| 久久影院123| 久久精品亚洲熟妇少妇任你| tube8黄色片| 90打野战视频偷拍视频| 色婷婷av一区二区三区视频| 丝袜美足系列| av在线app专区| 亚洲精品美女久久久久99蜜臀| 午夜福利在线观看吧| e午夜精品久久久久久久| 啦啦啦视频在线资源免费观看| 丰满饥渴人妻一区二区三| 日韩精品免费视频一区二区三区| 亚洲黑人精品在线| 亚洲欧洲日产国产| 一级黄色大片毛片| 欧美老熟妇乱子伦牲交| 成人手机av| 欧美少妇被猛烈插入视频| 美国免费a级毛片| 欧美+亚洲+日韩+国产| 欧美日本中文国产一区发布| 一进一出抽搐动态| 亚洲av成人一区二区三| 欧美亚洲 丝袜 人妻 在线| 日本wwww免费看| 国产高清国产精品国产三级| 国产亚洲av高清不卡| av视频免费观看在线观看| 精品福利永久在线观看| 亚洲av男天堂| av电影中文网址| 99久久人妻综合| 国产三级黄色录像| 国产欧美日韩综合在线一区二区| 日日摸夜夜添夜夜添小说| 国产精品久久久久久精品古装| 欧美激情 高清一区二区三区| 最近最新免费中文字幕在线| 两个人看的免费小视频| 亚洲专区中文字幕在线| 性高湖久久久久久久久免费观看| 日本撒尿小便嘘嘘汇集6| 大片免费播放器 马上看| 亚洲 国产 在线| av有码第一页| 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 国产一区有黄有色的免费视频| 狠狠婷婷综合久久久久久88av| 老司机午夜十八禁免费视频| 欧美 亚洲 国产 日韩一| 国产亚洲精品第一综合不卡| 国产男女超爽视频在线观看| 国产麻豆69| 真人做人爱边吃奶动态| 国产在线观看jvid| 亚洲一区二区三区欧美精品| 91老司机精品| 韩国精品一区二区三区| 肉色欧美久久久久久久蜜桃| 午夜成年电影在线免费观看| 亚洲熟女毛片儿| 欧美老熟妇乱子伦牲交| 18禁黄网站禁片午夜丰满| 日日爽夜夜爽网站| av在线app专区| 中文字幕人妻丝袜制服| 夫妻午夜视频| 欧美日韩黄片免| 中文字幕人妻丝袜一区二区| 午夜老司机福利片| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成av片中文字幕在线观看| 嫁个100分男人电影在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲精品在线美女| 午夜福利免费观看在线| 成人国产av品久久久| 免费在线观看影片大全网站| 国产一区二区三区在线臀色熟女 | 久热这里只有精品99| 久久久久久久久久久久大奶| 亚洲欧美日韩高清在线视频 | 日韩,欧美,国产一区二区三区| 操出白浆在线播放| 亚洲成人手机| 久久精品熟女亚洲av麻豆精品| 天天操日日干夜夜撸| 日本猛色少妇xxxxx猛交久久| 一级黄色大片毛片| 女性生殖器流出的白浆| 精品一区二区三区四区五区乱码| 久久久久网色| 欧美黑人欧美精品刺激| 视频区图区小说| netflix在线观看网站| 女人高潮潮喷娇喘18禁视频| 满18在线观看网站| 高清在线国产一区| 老司机靠b影院| 亚洲中文av在线| 亚洲国产av影院在线观看| 亚洲精品国产av成人精品| 国产成人系列免费观看| 最新的欧美精品一区二区| 久久热在线av| 精品人妻一区二区三区麻豆| 交换朋友夫妻互换小说| 国产人伦9x9x在线观看| 狠狠精品人妻久久久久久综合| 成人免费观看视频高清| 99久久精品国产亚洲精品| 不卡av一区二区三区| avwww免费| 亚洲精品自拍成人| 日本av手机在线免费观看| 老熟妇仑乱视频hdxx| 国产精品一区二区在线观看99| 亚洲成国产人片在线观看| 女性被躁到高潮视频| 一区二区日韩欧美中文字幕| 国产精品偷伦视频观看了| 看免费av毛片| 国产日韩欧美在线精品| 青草久久国产| 人妻人人澡人人爽人人| 国产不卡av网站在线观看| 国产色视频综合| 午夜福利视频精品| 欧美黑人精品巨大| 99精品久久久久人妻精品| 午夜日韩欧美国产| 久久精品熟女亚洲av麻豆精品| 丝袜人妻中文字幕| 午夜日韩欧美国产| 日韩欧美一区视频在线观看| 精品国产乱码久久久久久小说| 国产免费一区二区三区四区乱码| 蜜桃国产av成人99| 婷婷色av中文字幕| 亚洲人成电影免费在线| 狂野欧美激情性xxxx| 99久久综合免费| 国产成人欧美| 午夜影院在线不卡| 90打野战视频偷拍视频| 日韩制服骚丝袜av| 欧美大码av| 久久久精品区二区三区| 极品少妇高潮喷水抽搐| 1024香蕉在线观看| 亚洲欧美精品自产自拍| 两人在一起打扑克的视频| 在线 av 中文字幕| 国精品久久久久久国模美| 精品亚洲成a人片在线观看| 99国产精品一区二区蜜桃av | 国产黄色免费在线视频| 亚洲欧美一区二区三区久久| 国产又爽黄色视频| 欧美日韩成人在线一区二区| 国产精品久久久久成人av| 80岁老熟妇乱子伦牲交| 操美女的视频在线观看| 啦啦啦啦在线视频资源| 国产免费福利视频在线观看| 看免费av毛片| 在线永久观看黄色视频| 色婷婷av一区二区三区视频| 亚洲精品中文字幕在线视频| 久久99一区二区三区| 国产成人av激情在线播放| 12—13女人毛片做爰片一| 国产一区二区三区在线臀色熟女 | 捣出白浆h1v1| 熟女少妇亚洲综合色aaa.| 久久久精品区二区三区| 高清黄色对白视频在线免费看| 天堂8中文在线网| av在线app专区| 久久精品成人免费网站| 亚洲欧美清纯卡通| 黄片小视频在线播放| 亚洲欧美成人综合另类久久久| 欧美精品啪啪一区二区三区 | 美女扒开内裤让男人捅视频| 99久久人妻综合| av超薄肉色丝袜交足视频| 一级毛片精品| 午夜福利乱码中文字幕| 日韩一卡2卡3卡4卡2021年| 777久久人妻少妇嫩草av网站| 十八禁网站网址无遮挡| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一二三区在线看| 国产一区二区三区在线臀色熟女 | 午夜福利免费观看在线| 精品第一国产精品| 久久国产精品男人的天堂亚洲| 91国产中文字幕| 欧美乱码精品一区二区三区| 99精品欧美一区二区三区四区| 日韩一卡2卡3卡4卡2021年| 久久人人爽人人片av| 亚洲人成77777在线视频| 午夜成年电影在线免费观看| 一个人免费看片子| bbb黄色大片| 水蜜桃什么品种好| 久久中文字幕一级| 黄色毛片三级朝国网站| 国产一卡二卡三卡精品| 欧美另类亚洲清纯唯美| 91字幕亚洲| 男女边摸边吃奶| 男男h啪啪无遮挡| 日韩欧美一区二区三区在线观看 | 国产精品一区二区在线观看99| 精品一区在线观看国产| tube8黄色片| 波多野结衣一区麻豆| 亚洲一区中文字幕在线| 成人免费观看视频高清| 50天的宝宝边吃奶边哭怎么回事| 久热这里只有精品99| 精品国产一区二区久久| 久久久久久免费高清国产稀缺| 肉色欧美久久久久久久蜜桃| 久久久久久久精品精品| 另类亚洲欧美激情| 欧美一级毛片孕妇| 国产野战对白在线观看| 精品国产一区二区久久| 建设人人有责人人尽责人人享有的| 俄罗斯特黄特色一大片| 国产免费视频播放在线视频| 久久亚洲国产成人精品v| 亚洲免费av在线视频| 亚洲中文字幕日韩| 一级片免费观看大全| 午夜91福利影院| 久久中文看片网| 国产精品亚洲av一区麻豆| 欧美日韩黄片免| 欧美精品一区二区大全| 亚洲欧美清纯卡通| 在线天堂中文资源库| 久久久精品94久久精品| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| a级毛片黄视频| 免费高清在线观看日韩| 亚洲色图综合在线观看| 在线精品无人区一区二区三| 欧美国产精品一级二级三级| 亚洲第一av免费看| 99热网站在线观看| 日韩制服丝袜自拍偷拍| 免费在线观看黄色视频的| 国产视频一区二区在线看| 免费观看av网站的网址| 9191精品国产免费久久| 国产免费现黄频在线看| 脱女人内裤的视频| 91精品三级在线观看| svipshipincom国产片| 国产日韩欧美亚洲二区| 69精品国产乱码久久久| 激情视频va一区二区三区| 最近最新中文字幕大全免费视频| 又黄又粗又硬又大视频| 免费高清在线观看日韩| 国产成人a∨麻豆精品| 亚洲国产精品999| 老司机福利观看| 日本五十路高清| 国产av又大| 国产精品免费视频内射| 天堂中文最新版在线下载| 午夜福利视频精品| 国产一区二区三区综合在线观看| 久久久久久久久免费视频了| 黄片小视频在线播放| 国产男女内射视频| 成年动漫av网址| 亚洲成国产人片在线观看| 日韩制服丝袜自拍偷拍| 热99re8久久精品国产| 999久久久精品免费观看国产| 一级,二级,三级黄色视频| 热99国产精品久久久久久7| 免费高清在线观看日韩| 男男h啪啪无遮挡| 亚洲精品美女久久久久99蜜臀| 国产高清videossex| 国产一区有黄有色的免费视频| 国产精品秋霞免费鲁丝片| 在线观看www视频免费| 最近最新免费中文字幕在线| 欧美精品亚洲一区二区| 亚洲一码二码三码区别大吗| 人妻人人澡人人爽人人| 国内毛片毛片毛片毛片毛片| 亚洲av成人一区二区三| 国产男女超爽视频在线观看| 99九九在线精品视频| 亚洲中文av在线| 精品少妇久久久久久888优播| 在线永久观看黄色视频| 多毛熟女@视频| 免费少妇av软件| 久久久久久亚洲精品国产蜜桃av| 日本91视频免费播放| 日韩一卡2卡3卡4卡2021年| 999久久久国产精品视频| 亚洲色图 男人天堂 中文字幕| 欧美成人午夜精品| 日本精品一区二区三区蜜桃| 久久午夜综合久久蜜桃| 久久av网站| 国产高清国产精品国产三级| 国产91精品成人一区二区三区 | 侵犯人妻中文字幕一二三四区| 男女免费视频国产| 99国产精品一区二区三区| 日韩视频一区二区在线观看| 精品久久久精品久久久| 欧美激情高清一区二区三区| 欧美日韩国产mv在线观看视频| 久久久国产精品麻豆| 一级毛片女人18水好多| 亚洲精品美女久久久久99蜜臀| 国产一区二区 视频在线| 国产av一区二区精品久久| 美女扒开内裤让男人捅视频| 两个人免费观看高清视频| 自拍欧美九色日韩亚洲蝌蚪91| av在线老鸭窝| 亚洲国产欧美在线一区| 精品一品国产午夜福利视频| √禁漫天堂资源中文www| 久久精品熟女亚洲av麻豆精品| 国产男人的电影天堂91| 国产精品香港三级国产av潘金莲| 2018国产大陆天天弄谢| 丁香六月天网| 777米奇影视久久| 日本91视频免费播放| 亚洲免费av在线视频| 男女午夜视频在线观看| 99热网站在线观看| 国产在线观看jvid| 国产精品免费大片| 久久久欧美国产精品| 久久ye,这里只有精品| 亚洲精品久久成人aⅴ小说| 777米奇影视久久| 国产又爽黄色视频| 欧美中文综合在线视频| 两个人免费观看高清视频| 国产精品 欧美亚洲| 国产精品香港三级国产av潘金莲| 欧美另类一区| 亚洲国产精品999| 亚洲情色 制服丝袜| 高清视频免费观看一区二区| 啦啦啦在线免费观看视频4| 在线天堂中文资源库| 窝窝影院91人妻| 国产一区二区激情短视频 | 性色av乱码一区二区三区2| 天天躁日日躁夜夜躁夜夜| 午夜福利在线免费观看网站| 亚洲欧美一区二区三区久久| 亚洲欧美精品自产自拍| 亚洲一区中文字幕在线| 中国美女看黄片| 麻豆乱淫一区二区| 汤姆久久久久久久影院中文字幕| 不卡一级毛片| 午夜老司机福利片| 一个人免费看片子| 日本猛色少妇xxxxx猛交久久| 日韩欧美一区视频在线观看| 亚洲av电影在线观看一区二区三区| 久久天堂一区二区三区四区| 国产成人av激情在线播放| 精品人妻熟女毛片av久久网站| 淫妇啪啪啪对白视频 | 丰满人妻熟妇乱又伦精品不卡| av欧美777| 黑人猛操日本美女一级片| 国产又色又爽无遮挡免| 我的亚洲天堂| 色94色欧美一区二区| 在线观看舔阴道视频| 丝袜美腿诱惑在线| 色老头精品视频在线观看| 久久久久久久久久久久大奶| 国产亚洲一区二区精品| 亚洲一卡2卡3卡4卡5卡精品中文| 99国产极品粉嫩在线观看| 成人国语在线视频| 欧美97在线视频| √禁漫天堂资源中文www| 蜜桃在线观看..| 淫妇啪啪啪对白视频 | 国产一区二区 视频在线| 亚洲国产精品999| 免费在线观看日本一区| videosex国产| 亚洲欧美清纯卡通| 日韩人妻精品一区2区三区| 午夜两性在线视频| 国产极品粉嫩免费观看在线| 人人妻,人人澡人人爽秒播| e午夜精品久久久久久久| 亚洲第一av免费看| 成人影院久久| 久久久久精品人妻al黑| 欧美乱码精品一区二区三区| 久久99一区二区三区| 久久久久国产精品人妻一区二区| 桃红色精品国产亚洲av| 男人舔女人的私密视频| 在线av久久热| 久久亚洲精品不卡| 制服人妻中文乱码| 日韩,欧美,国产一区二区三区| 亚洲国产欧美网| 免费观看人在逋| a 毛片基地| 亚洲精品久久久久久婷婷小说| 国产一级毛片在线| 久久久久久亚洲精品国产蜜桃av| 制服人妻中文乱码| 老汉色av国产亚洲站长工具| 亚洲av片天天在线观看| 69av精品久久久久久 | 美女高潮喷水抽搐中文字幕| 啦啦啦视频在线资源免费观看| 美女脱内裤让男人舔精品视频| 成年人午夜在线观看视频| 中亚洲国语对白在线视频| 亚洲三区欧美一区| 女警被强在线播放| 国产xxxxx性猛交| 成人影院久久| 久久久国产精品麻豆| 国产精品久久久久成人av| 久久免费观看电影| 国产高清视频在线播放一区 | 亚洲视频免费观看视频| 国产免费视频播放在线视频| 最近最新中文字幕大全免费视频| 国产亚洲一区二区精品| 日韩一卡2卡3卡4卡2021年| 精品一区二区三卡| 淫妇啪啪啪对白视频 | av在线播放精品| 国内毛片毛片毛片毛片毛片| 亚洲五月色婷婷综合| 后天国语完整版免费观看| 午夜免费成人在线视频| 黑人猛操日本美女一级片| 老汉色av国产亚洲站长工具| 一本大道久久a久久精品| 一级,二级,三级黄色视频| 12—13女人毛片做爰片一| 夜夜夜夜夜久久久久| 日韩精品免费视频一区二区三区| 黄色a级毛片大全视频| 午夜福利在线免费观看网站| 国产精品99久久99久久久不卡| 国产欧美日韩一区二区精品| 午夜福利视频在线观看免费| 正在播放国产对白刺激| 色94色欧美一区二区| 色老头精品视频在线观看| 欧美成狂野欧美在线观看| 97在线人人人人妻| 精品一区在线观看国产| 色94色欧美一区二区| 亚洲国产精品成人久久小说| 国产熟女午夜一区二区三区| 国产成人系列免费观看| 久久综合国产亚洲精品| 久久中文字幕一级| 亚洲成人免费电影在线观看| 日本vs欧美在线观看视频| 欧美黑人精品巨大| 窝窝影院91人妻| 亚洲七黄色美女视频| 亚洲成人免费电影在线观看| 亚洲av男天堂| 亚洲成人国产一区在线观看| 欧美日韩国产mv在线观看视频| 男女边摸边吃奶| 男人舔女人的私密视频| 欧美中文综合在线视频| 五月开心婷婷网| 在线十欧美十亚洲十日本专区| 大香蕉久久成人网| 极品少妇高潮喷水抽搐| 青春草亚洲视频在线观看| 日韩熟女老妇一区二区性免费视频| 美女国产高潮福利片在线看| 99精品久久久久人妻精品| 欧美在线黄色| av在线播放精品| 在线看a的网站| 黄色 视频免费看| 可以免费在线观看a视频的电影网站| 国产欧美日韩一区二区三 | 国产欧美日韩精品亚洲av| 亚洲精品国产av成人精品| 桃红色精品国产亚洲av| 91成人精品电影| 男女国产视频网站| videosex国产| 一级黄色大片毛片| 日韩一卡2卡3卡4卡2021年| 他把我摸到了高潮在线观看 | 亚洲va日本ⅴa欧美va伊人久久 | 欧美日韩精品网址| 国产国语露脸激情在线看| 久久精品熟女亚洲av麻豆精品| 国产av国产精品国产| 69av精品久久久久久 | 色94色欧美一区二区| 国产精品免费视频内射| 极品人妻少妇av视频| 久久毛片免费看一区二区三区| 亚洲中文日韩欧美视频| 纵有疾风起免费观看全集完整版| 午夜精品国产一区二区电影| 多毛熟女@视频| 日韩三级视频一区二区三区| 下体分泌物呈黄色| 一区在线观看完整版| 美女中出高潮动态图| 中文字幕av电影在线播放| 永久免费av网站大全| 亚洲专区中文字幕在线| 男女高潮啪啪啪动态图| 国产精品一区二区在线观看99| 又紧又爽又黄一区二区| 亚洲国产av影院在线观看| 国产91精品成人一区二区三区 | tube8黄色片| 午夜免费鲁丝| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 久久午夜综合久久蜜桃| 性色av乱码一区二区三区2| 悠悠久久av| 国产成人啪精品午夜网站| 女人久久www免费人成看片| 亚洲专区国产一区二区| 女人精品久久久久毛片| www.999成人在线观看| 99久久99久久久精品蜜桃| 视频在线观看一区二区三区| 精品亚洲成国产av| 国产精品一区二区免费欧美 | 韩国高清视频一区二区三区| 国产精品亚洲av一区麻豆| 老汉色av国产亚洲站长工具| 97人妻天天添夜夜摸| 亚洲伊人久久精品综合| 欧美黑人欧美精品刺激| 一边摸一边抽搐一进一出视频| 午夜老司机福利片| 欧美激情 高清一区二区三区| 欧美性长视频在线观看| 后天国语完整版免费观看| 午夜福利乱码中文字幕| 亚洲欧美一区二区三区久久| 制服诱惑二区| avwww免费| netflix在线观看网站| 免费人妻精品一区二区三区视频| 波多野结衣av一区二区av| 日日夜夜操网爽| 中国美女看黄片| 国产精品麻豆人妻色哟哟久久| 久久99热这里只频精品6学生| av超薄肉色丝袜交足视频| 涩涩av久久男人的天堂| 免费在线观看日本一区| 男女边摸边吃奶| 亚洲国产毛片av蜜桃av| 国产1区2区3区精品| 视频在线观看一区二区三区| 欧美精品高潮呻吟av久久| 中文字幕av电影在线播放| 汤姆久久久久久久影院中文字幕| 巨乳人妻的诱惑在线观看| 黄色毛片三级朝国网站| 十八禁人妻一区二区| 国产xxxxx性猛交| 99久久99久久久精品蜜桃| 成人黄色视频免费在线看| 日韩欧美一区视频在线观看| 日韩免费高清中文字幕av|