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

    黃河中游綠水系數(shù)變化及其生態(tài)環(huán)境意義

    2015-02-05 06:03:40許炯心
    生態(tài)學報 2015年22期
    關鍵詞:藍水綠水梯田

    許炯心

    中國科學院地理科學與資源研究所, 中國科學院陸地水循環(huán)及地表過程重點實驗室, 北京 100101

    黃河中游綠水系數(shù)變化及其生態(tài)環(huán)境意義

    許炯心*

    中國科學院地理科學與資源研究所, 中國科學院陸地水循環(huán)及地表過程重點實驗室, 北京 100101

    綠水和藍水的概念和理論對于半干旱、半濕潤區(qū)水資源的評價和管理有重要意義。引入流域尺度上綠水系數(shù)指標,定義為某一流域內由降水到綠水的轉換系數(shù)。研究發(fā)現(xiàn),1950—2011年間,黃河中游河口鎮(zhèn)至龍門區(qū)間(河龍區(qū)間)的綠水系數(shù)呈現(xiàn)增大的趨勢。除了氣溫、降水的影響外,水土保持是一個重要因素。水土保持措施減少了降雨到徑流(藍水)的轉化率,增大了降雨到綠水的轉化率。綠水系數(shù)的增大意味著坡面徑流減弱和河流徑流的減弱,從而減少了坡面侵蝕和流域產沙。同時,綠水系數(shù)的增大意味著植被蒸騰作用的增強,說明植被對地表的保護作用增強,這也會導致坡面侵蝕的減弱。河龍區(qū)間產沙量與流域綠水系數(shù)之間呈顯著的負相關關系,產沙量變化的53.7%可以用流域綠水系數(shù)的變化來解釋。水土保持實施后,梯田和壩地使得生產性綠水系數(shù)增大,從而增加了糧食產量。梯田、壩地面積與綠水系數(shù)和糧食產量之間都存在著顯著的正相關關系,糧食產量與綠水系數(shù)之間也存在顯著的正相關關系。綠水系數(shù)具有流域生態(tài)環(huán)境質量變化的指示意義,在年降水可比的情況下,流域綠水系數(shù)的減小意味著集水區(qū)生態(tài)環(huán)境環(huán)境質量降低,流域綠水系數(shù)的增大意味著集水區(qū)生態(tài)環(huán)境質量提高。依照綠水系數(shù)的變化,可以將近60余年來河龍區(qū)間生態(tài)環(huán)境的變化過程劃分為3個階段。

    綠水;綠水系數(shù);氣候變化;人類活動;黃河中游

    為了更好地進行半干旱、半濕潤區(qū)農業(yè)水資源的評價,F(xiàn)alkenmark等于1995年提出綠水和藍水的概念。藍水是指儲存在河流、湖泊以及含水層中的水,而綠水是指直接來源于降水并用于蒸散的水[1]。此后,綠水的研究在世界上受到重視。綠水可以被定義為蒸散流,是流向大氣圈的水汽流,包括灌溉農田、濕地、水面和自然植被等不同地表的蒸散發(fā)產生的水汽流[2-3],也可以被定義為具體的資源,即綠水是源于降水、存儲于土壤并被植被蒸散發(fā)消耗的水資源,某一地區(qū)在一段時間內能夠獲得的總的綠水資源量等于該時期內蒸散發(fā)累計量[4-5]。上述理論的提出,豐富了水資源的內涵,為更加科學地進行水資源管理提供了理論依據(jù)?;谶@一理論,水資源的管理應分為藍水管理和綠水管理。研究表明,從水循環(huán)的角度分析,全球尺度上總降水的65%通過森林、草地、濕地和雨養(yǎng)農田的蒸散返回到大氣中,成為綠水,僅有35%的降水儲存于河流、湖泊以及含水層中,成為藍水[5]。綠水可分為兩個分量:一部分為植物蒸騰量,與陸地生態(tài)系統(tǒng)中的生物量生產密切相關,可稱為生產性綠水;另一部分為陸面蒸發(fā)量,可稱為非生產性藍水[6]。國際上對于綠水的研究,已取得很多成果[5- 10]。近年來,劉昌明、程國棟等在我國倡導展開綠水的研究[11-12]。從某種意義上說,傳統(tǒng)的水資源管理屬于藍水管理,目前亟待加強對于綠水的管理。這對于提高我國水資源管理水平具有重要意義。綠水、藍水的概念引入我國之后,已產生了一些研究成果[13- 15],但尚待展開深入的研究。

    已有研究大部分集中于農田尺度。為了進行流域綠水管理,必須在流域尺度上深入研究綠水形成的機理,研究降水-綠水的轉化率的變化如何受到自然與人為因素的影響。流域系統(tǒng)既是一個水文地貌系統(tǒng),又是一個生態(tài)系統(tǒng)。目前國內外研究者對于綠水在水文學和水資源管理上的意義進行了比較充分的研究,但對于其在流域尺度上生態(tài)意義的揭示還較少。近60余年黃河中游河口鎮(zhèn)至龍門區(qū)間流域綠水系數(shù)有明顯的減小趨勢,并從氣候變化和人類活動的影響揭示了減小的原因[16]。本文是已有工作[16]的延伸,將在前文研究的基礎上,進一步探討綠水系數(shù)增大的生態(tài)環(huán)境意義及其在流域水資源管理和水土流失治理效益評價中的應用意義。其目的是深化對于綠水系數(shù)科學內涵的認識,為更科學、更確切地評價水土保持效益和更好地開展流域水資源管理和水土保持規(guī)劃提供新的知識。

    1 研究區(qū)概況

    黃河中游河口鎮(zhèn)至龍門區(qū)間(以下簡稱河龍區(qū)間)(圖1)是黃土高原的主體。屬中溫帶和暖溫帶半干旱區(qū),自然植被類型區(qū)為溫帶干草原和暖溫帶森林草原。年降水量為350—450mm,多以暴雨形式集中下落,降雨侵蝕力較強。地表物質主要為黃土,黃土層深厚,滲透性強,雨季入滲的雨水只有一小部分補給到地下水中,其余部分在旱季中就地蒸發(fā),因而徑流系數(shù)較低,徑流深度較小。黃河在河口鎮(zhèn)流入該區(qū),在龍門流出該區(qū),河龍區(qū)間的流域面積為112000km2。按照1950—1999年間的統(tǒng)計,河口鎮(zhèn)水文站和龍門水文站的多年平均實測徑流量分別229.1×108m3和283.3×108m3,河龍區(qū)間產生的徑流量為54.2×108m3。該區(qū)產流量占花園口以上流域產流量的13.9%。河龍區(qū)間是黃河中游多沙粗沙區(qū)的主體,位于黃土高原丘陵溝壑區(qū),黃土層厚度大,溝谷密度大、切割深,侵蝕強烈,是黃河泥沙特別是粗泥沙的主要來源區(qū)。由于黃土層深厚,滲透性強,故徑流深度較小,對黃河徑流量的貢獻相對較小。該區(qū)的面積只占花園口站以上流域面積的18.4%,年產流量占花園口以上流域產流量的14.9%,年產沙量卻占花園口以上年產沙量的55.8%[17]。本區(qū)是黃河流域水土流失重點治理區(qū),從20世紀60年代以來梯田、林草和淤地壩措施的大規(guī)模實施以及1998年以后退耕還林(草)的開展,使流域下墊面特性發(fā)生了很大的變化。與此同時,氣溫升高,降水量減少,蒸發(fā)特性改變,使得流域產流特性發(fā)生了很大的變化,進而導致了綠水、藍水關系的變化。由于特殊的自然地理條件,該地區(qū)的產流過程對于環(huán)境因素的變化十分敏感,是研究水文過程對氣候變化和人類活動響應的理想地區(qū)。

    圖1 黃河流域示意圖Fig.1 A map of the Yellow River basin

    2 研究方法

    如前所述,綠水包括生產性綠水(指被蒸騰量)和非生產性綠水(陸面蒸發(fā)量)。其中,生產性綠水對于水資源管理有更為重要的意義。然而,目前在廣大流域面積上對于植被蒸騰量的觀測資料很少,不足以進行深入研究。因此,本文研究包括蒸騰量和蒸發(fā)量在內的廣義綠水。按水量平衡方程:

    降水量=蒸騰蒸發(fā)量+徑流量+蓄水變量

    (1)

    在較長的時間尺度上,流域中的蓄水變量可以忽略。降水量轉化為徑流量和蒸散發(fā)量,前者為藍水,后者為綠水。由此可寫出:

    降水量=蒸騰蒸發(fā)量+徑流量=藍水量+綠水量

    (2)

    引入流域尺度上綠水系數(shù)(Cgw)的指標,定義為從降水到綠水的轉化率,即綠水量與降水量之比。因此,

    綠水系數(shù)=綠水量/降水量=(降水量-藍水量)/降水量=(降水量-徑流量)/降水量

    (3)

    必須指出,在用上式研究綠水系數(shù)時,應該采用天然徑流量而不是實測徑流量,因為人類所引用的水量屬于藍水,引水導致實測徑流量及藍水減小,對于這一部分必須進行還原計算。因此:

    綠水系數(shù)=綠水量/降水量 =(降水量-天然徑流量)/降水量

    (4)

    本文在年尺度上研究綠水系數(shù)。應該指出,在年尺度上水庫對徑流的調蓄是一種重要的蓄水變量,不宜忽略。在進行年天然徑流還原計算時,已扣除了水庫蓄水變量。流域土壤的調蓄量對于水量平衡也有一定影響,但目前尚難以得到長系列、大尺度的歷年土壤蓄水變量的資料。作為近似,本文暫時忽略了這一因素。

    本文采用統(tǒng)計方法與時間序列分析方法,研究綠水系數(shù)的變化及其與影響因素的關系,揭示綠水系數(shù)變化所包含的生態(tài)環(huán)境意義。同時還采用Mann-Kendall方法研究了綠水系數(shù)的變化趨勢。Mann-KendallU值隨時間的變化反映變量變化趨勢的改變,可以用來探測由增到減(或由減到增)的轉折點與突變點[18]。

    3 資料來源

    河口鎮(zhèn)至龍門區(qū)間的年徑流量由龍門站的年徑流量減去河口鎮(zhèn)站的年徑流量而得到,時間系列為1950年至20011年。所依據(jù)的徑流資料來自河口鎮(zhèn)水文站和龍門水文站。水文站實測的徑流量(Qwm),并不代表流域實際產生的徑流量即天然徑流量。黃河水利委員會有關部門統(tǒng)計出了不同年份某一站點以上流域中的灌溉引水量、工業(yè)及城鎮(zhèn)用水量、水庫調蓄量,并對人類引水后回歸到河道中的水量進行了估算,求出了歷年全流域的人類凈引水量Qw,div,計算出了天然徑流量(Qwn):Qwn=Qwm-Qw,div+ΔSr。由此,可得到:

    Qwn,H-L=Qwm,H-L+Qw,div,H-L=Qwm,L-Qwm,H+Qw,div,H-L

    (5)

    式中,Qwn,H-L為河龍區(qū)間天然徑流量,Qwm,H-L為河龍區(qū)間實測徑流量,Qwm,L為龍門站實測徑流量,Qwm,H為河口鎮(zhèn)站實測徑流量,ΔSr為水庫調蓄量,Qw,div,H-L為河龍區(qū)間人類凈引水量。以上均按年系列計算。河龍區(qū)間的綠水系數(shù)Cgw由下式計算:

    Cgw=Qwn,H-L/(AH-L×PH-L)

    (6)

    式中,Qwn,H-L為河龍區(qū)間天然年徑流量(108m3/a),AH-L為河龍區(qū)間的流域面積(km2),PH-L為河龍區(qū)間的面平均年降水量(mm)。河龍區(qū)間年降水量、天然徑流量、凈引水量的資料來自歷年《黃河水資源公報》;各項水土保持面積的資料來自黃河水利科學研究院;氣溫資料來自中國氣象局。由于河龍區(qū)間的干流河道位于陜西、山西間的峽谷中,河段沖淤量很少,可以忽略,河龍區(qū)間的年產沙量按龍門站的年輸沙量減去河口鎮(zhèn)(頭道拐)站的年輸沙量來計算,兩站年輸沙量的資料來自歷年《黃河泥沙公報》。實測徑流量、輸沙量和氣象資料在觀測時嚴格遵照國家技術標準,并進行了精度檢驗;水土保持面積在統(tǒng)計數(shù)據(jù)的基礎上進行了抽樣核實,具有一定的準確性。文中還利用了陜西、山西兩省歷年糧食產量的數(shù)據(jù),該數(shù)據(jù)來自中國統(tǒng)計年鑒[19];利用了山西、陜西兩省2000—2006年間退耕及封禁治理面積的資料,數(shù)據(jù)來自中國統(tǒng)計年鑒[19]和黃河年鑒[20]。

    4 結果與分析

    4.1 綠水系數(shù)和影響因素的時間變化趨勢

    流域綠水系數(shù)與氣候因素(如降水量和氣溫)和人類活動因素(如各種水土保持措施)有密切的關系。圖2點繪了河龍區(qū)間綠水系數(shù)隨時間的變化。圖中顯示,綠水系數(shù)呈現(xiàn)增大的趨勢,與時間(年份)的決定系數(shù)R2=0.5276,顯著性概率P<0.001。圖2還點繪了河龍區(qū)間流域因素如年降水量和年平均氣溫以及梯田面積、造林面積、種草面積和壩地面積等隨時間的變化。梯田、造林、種草和壩地面積均表現(xiàn)出顯著的增大趨勢,R2分別為0.9801、0.9114、0.7598和0.8515,顯著性概率P<0.001。年降水量略有減小的趨勢,R2為0.061 (P=0.064,P<0.10);年平均氣溫呈現(xiàn)顯著的增大趨勢,R2為0.4522(P<0.001)。在先前的研究中對于河龍區(qū)間綠水系數(shù)增大的原因進行了較深入的研究[16],表明綠水系數(shù)的增大與夏季風強度減弱、降水量減小、年均氣溫升高和水土保持面積增大有密切的相關關系。建立了綠水系數(shù)(Cgw)與水土保持措施總面積(Asw)、5a滑動平均年降水量(Pm,5m)和年平均氣溫(Tm,5m)之間的多元線性回歸方程,在此基礎上估算出Asw、Pm,5m、Tm,5m對Cgw的貢獻率分別為51%, 37% 和 12%。本文進一步揭示綠水系數(shù)增大的生態(tài)環(huán)境意義。

    圖2 河龍區(qū)間綠水系數(shù)(1950—2011年)和影響因素(1950—2006 年)及隨時間的變化Fig.2 Temporal variations in green water coefficient from 1950 to 2011 and in the influencing factors (1950—2006) including annual precipitation, temperature and the areas of 4 types of soil and water conservation measures

    4.2 綠水系數(shù)與流域侵蝕產沙的關系

    從實質上說,水土保持措施的作用是對于地表徑流進行調節(jié)。這種調節(jié)包括兩個方面。一是對于降雨-徑流轉換率即地表徑流系數(shù)的調節(jié)。梯田、造林種草等水土保持措施會增加降雨入滲,增加土壤水,減少地表徑流。由于黃土高原黃土層厚度很大,地下水埋藏很深,入滲的雨水量經(jīng)由土壤水和地下水最終轉化為河流基流的比例不大,最后大部分都消耗于蒸發(fā)與蒸騰作用,重新回到大氣,即轉化為綠水。因此,水土保持會顯著地增大綠水系數(shù)。二是對于徑流的影響。坡地改為水平梯田或緩坡梯田之后,地表坡度大大減小;林草植被則增大了坡面地表的糙率系數(shù),這兩方面的作用都會減慢坡面流的流速、延長其歷時,從而增加徑流運動過程中的下滲量。淤地壩以及前期淤成的壩地大量攔蓄來自溝道流域的洪水徑流,這些徑流大部分都滲入土壤之中,轉化為土壤水,最后消耗于蒸散發(fā)而成為綠水。因此,水土保持措施對于徑流的調節(jié),最后會表現(xiàn)為增大綠水系數(shù)[16]。圖3中點繪的綠水系數(shù)與各項水土保持措施面積的關系,證明了這一點。圖中的R2表明,梯田面積、造林面積、種草面積和淤地壩造地面積等變量可以解釋綠水系數(shù)變化的39.5%—49.1%。

    圖3 河龍區(qū)間水土保持措施面積與綠水系數(shù)之間的關系Fig.3 Relationship between green water coefficient and the area of soil and water conservation measures

    從本質上說,水土保持措施減少了降雨到徑流(藍水)的轉化率,增大了降雨到綠水的轉化率。綠水系數(shù)的減小意味著坡面徑流減弱,河流徑流也減弱,前者可以減少坡面侵蝕,后者則可以減少河道侵蝕。同時,綠水系數(shù)的增大意味著植被蒸騰作用的增強,說明植被對地表的保護作用增強,這也會導致坡面侵蝕的減弱。因此,河流的產沙量與綠水系數(shù)之間存在密切的關系。圖4顯示,河龍區(qū)間產沙量與綠水系數(shù)之間呈顯著的負相關關系(R2=0.5372,P<0.001),意味著河龍區(qū)間產沙量變化的53.7%可以用綠水系數(shù)的變化來解釋。圖4中對于二者的時間變化進行了比較。圖中顯示,前者具有減小趨勢,后者則有增大趨勢。由于河龍區(qū)間降水存在3—5a的波動,年產沙量與綠水系數(shù)也有3—5a的波動。為了更好的體現(xiàn)出變化的趨勢,計算出年產沙量與綠水系數(shù)的5a滑動平均值,其變化曲線也繪在圖中??梢钥吹?,兩條曲線具有很好的鏡像關系。產沙量的5a滑動平均值與綠水系數(shù)的5a滑動平均值具有很強的負相關,R2=0.8266,比圖4的R2=0.5372要高得多。這一現(xiàn)象表明,植被通過對增大綠水系數(shù)來影響侵蝕產沙,具有一定的時間尺度效應。在較長的時間(例如5a)尺度上,如果綠水系數(shù)有增大趨勢,則意味著植被狀況持續(xù)改善,植被的抗蝕力也會持續(xù)增加,從而使侵蝕產沙量穩(wěn)定地減少。因此,在較長的時間尺度上(例如5a)綠水系數(shù)的變化對于產沙量變化的解釋能力要比年尺度關系的解釋能力強得多。

    圖4 河龍區(qū)間年產沙量與綠水系數(shù)之間的關系Fig.4 Relationships between the annual green water coefficient and sediment yield in the Hekouzhen-Longmen drainage area

    4.3 綠水系數(shù)增大在農業(yè)生產上的意義

    黃河流域綠水系數(shù)的增大意味著生產性綠水系數(shù)的增大。因此,在綠水系數(shù)增大的同時,植被覆蓋度和糧食產量均增大。由于梯田與坡耕地相比,降雨入滲量大大增加;壩地攔蓄的洪水,大部分最后都滲入地下,因此梯田和壩地土壤的含水率要大大高于坡耕地。據(jù)黃河水利委員會綏德水土保持試驗站實測資料,坡地、梯田、壩地土壤含水率分別為9.47%、10.72%和17.61%,梯田和壩地分別是是坡地的1.13倍和1.86倍[21]。在其他條件可比時,這是梯田和壩地的糧食產量要大大高于坡耕地的原因。據(jù)黃土高原典型調查資料,壩地糧食產量為3700—4500 kg/hm2,高的可達7500 kg以上,是坡耕地的4—6倍、梯田的2—3倍,梯田每公頃產量則為坡耕地的2倍以上。據(jù)綏德水土保持試驗站對王茂溝1960—1996年37年間的觀測,壩地、水平梯田、坡耕地年均糧食產量分別為每公頃4750、1606和566 kg,即壩地、水平梯田分別為坡耕地的8.39倍和2.83倍[21]。農作物耗水量體現(xiàn)為蒸騰作用,是流域綠水量的重要組成部分。在其他條件可比時,單位面積糧食產量越高,農作物的耗水量越大,梯田和壩地糧食的增產意味著綠水系數(shù)的增大。河龍區(qū)間流域絕大部分位于山西和陜西兩省,以這兩省的糧食產量的相對變化來近似反映河龍區(qū)間糧食產量的相對變化。以實施水土保持以前的20世紀50年代平均糧食產量作為基準,計算出歷年的產量與基準產量之比值,稱為糧食產量指數(shù)。以歷年數(shù)據(jù)為基礎,圖5中分別點繪了糧食產量指數(shù)與梯田面積與壩地面積的關系,表現(xiàn)出很強的正相關。圖5還顯示,糧食產量指數(shù)與河龍區(qū)間綠水系數(shù)也呈較顯著的正相關,說明水土保持實施后,生產性綠水占降水的比率增大,這是糧食產量顯著提高的原因之一。

    圖5 河龍區(qū)間糧食產量指數(shù)與梯田面積、壩地面積和綠水系數(shù)的關系(糧食產量數(shù)據(jù)來自文獻[20])Fig.5 Index of grain yield plotted against the area of terrace land, the area of the land created by check-dams and Cwg[20]

    4.4 綠水系數(shù)對生態(tài)環(huán)境的指示意義

    在蒸散發(fā)量中,蒸發(fā)量為非生產性綠水流,散發(fā)量(植物蒸騰量)為生產性綠水流[6]。其中,農作物散發(fā)量是生產糧食所消耗的綠水,非農作物植被的散發(fā)量則是維持生態(tài)系統(tǒng)的所消耗的綠水量,這兩部分綠水量都對生態(tài)環(huán)境有利。為了對這兩部分進行區(qū)分,可以稱前一部分為農業(yè)生產性綠水流,后一部分為生態(tài)性綠水流。因土壤蒸發(fā)而消耗的綠水量可以增加空氣的濕度,就廣義而言對于生態(tài)環(huán)境的維持也是有利的。因此,綠水流對于生態(tài)環(huán)境是有利的。從這一意義出發(fā),可以將綠水系數(shù)作為衡量生態(tài)環(huán)境變化的指標。在年降水可比的情況下,如果綠水系數(shù)減小,說明生態(tài)系統(tǒng)中以徑流的方式流失的水量即藍水的比率增大,用于其自身維持的水量(即綠水)的比率減小,在這一過程中與徑流流失相伴隨的物質流失(土壤顆粒、土壤有機質和營養(yǎng)元素)也增加??梢哉J為,這會導致生態(tài)環(huán)境質量的下降。反之,綠水系數(shù)增大,則說明生態(tài)系統(tǒng)中以徑流的方式流失的水量即藍水的比率減小(由此導致物質流失減少)、綠水的比率增大(這意味著生態(tài)系統(tǒng)的生產力提高),可以認為生態(tài)環(huán)境趨于好轉。從這一思路出發(fā),可以評價黃河中游河口鎮(zhèn)至龍門區(qū)間的綠水系數(shù)的變化及其生態(tài)環(huán)境意義。

    采用Mann-Kendall方法研究了綠水系數(shù)的變化趨勢,探測了變化過程中的轉折點。圖6點繪了河龍區(qū)間綠水系數(shù)的Mann-KendallU值隨時間的變化。從圖中可以看到,U值變化曲線有兩個明顯的轉折點,發(fā)生于1970年和1999年;1990年還有一個次一級的轉折點??梢詫?950—2011年間河龍區(qū)間綠水系數(shù)的變化分為4個階段:(1)階段Ⅰ(1950—1969年):綠水系數(shù)呈減小趨勢。這一階段中,水土保持措施尚未生效,人類對環(huán)境的破壞大于治理,流域水分流失相對嚴重,因而綠水系數(shù)減小。與水分流失相伴隨的土壤流失與土壤養(yǎng)分流失均較嚴重,這意味著生態(tài)環(huán)境的質量下降;(2)階段Ⅱ(1970—1998年):綠水系數(shù)在波動中呈增大趨勢。水土保持措施生效,流域水分流失大大減少,與水分流失相伴隨的土壤流失與土壤養(yǎng)分流失也大大減少,這意味著生態(tài)環(huán)境質量的提高。這一階段又可細分為兩個時段。1970—1989年(Ⅱa),水土保持的生效使得綠水系數(shù)增大。然而,1990—1998年(Ⅱb),淤地壩的攔沙效應明顯衰減。河龍區(qū)間的淤地壩,絕大部分是20世紀70年代修建的,80年代以后淤地壩修建量大為減少,而淤地壩的攔沙壽命為10—20a,70年代修建的淤地壩與攔沙庫到這一時期已大部失效[22],因此流域水分流失增加,綠水系數(shù)減?。?3)階段Ⅲ(1999—2011年):綠水系數(shù)迅速增大。這一階段中,除了原有的梯田、林草和淤地壩措施得到加強外,大面積退耕還林還草和以自然封禁為主的生態(tài)恢復也在這一地區(qū)廣泛開展,據(jù)中國統(tǒng)計年鑒資料[19],山西、陜西兩省2000—2006年間累積退耕耕地面積達300.4萬hm2;據(jù)黃河年鑒資料[20],河龍區(qū)間1996—2005年間累積封禁治理133.2萬hm2,這使得植被進一步改善,生態(tài)環(huán)境質量進一步提高。

    圖6 河龍區(qū)間綠水系數(shù)的Mann-Kendall U值、綠水系數(shù)距平及累積距平隨時間的變化Fig.6 Temporal variations in Mann-Kendall U of Cgw, the departure of Cgw to the average and the cumulative departure of Cgw to the average

    為了對于圖6中的變化階段作進一步論證,在圖6中分別點繪了綠水系數(shù)的距平百分數(shù)及其累積值的變化??梢钥吹剑?950—1970年間,綠水系數(shù)距平百分數(shù)以負值為主,21a中正值只出現(xiàn)了1次;1971—1999年間,綠水系數(shù)有所增大,28a中距平百分數(shù)正值出現(xiàn)了16次,負值出現(xiàn)了12次;2000—2011年間,綠水系數(shù)距平百分數(shù)進一步增大,11a中正值出現(xiàn)了10次,負值只出現(xiàn)了1次。綠水系數(shù)累積距平百分比曲線則顯示出完全相似的階段性,在1970年和1999年分別出現(xiàn)兩個轉折點,將整個變化過程分為3個階段。1950—1970年,綠水系數(shù)累積距平百分數(shù)減??;1971—1999年在總體上緩慢增大,但可以以1990年為界劃分為兩個次一級的階段,1970—1989年增大,1990—1999年減小;2000—2011年,綠水系數(shù)累積距平百分數(shù)呈現(xiàn)加速增大的趨勢。

    5 關于綠水系數(shù)應用意義的討論

    本文所揭示的黃河中游綠水系數(shù)變化的生態(tài)環(huán)境意義,有助于更全面地評價水土保持的水資源效應,從而為流域水資源的管理提供新的認識。這體現(xiàn)在4個方面:第一,綠水系數(shù)的增大意味著綠水(蒸發(fā)蒸騰量)在降水中所占比例的增大。蒸騰量比率的增大表明植物消耗水量的增加,即林草植被和農作物耗水的增加。前者屬于生態(tài)耗水量,對于生態(tài)環(huán)境的改善有利;后者是糧食等農作物生產的耗水量,有益于社會經(jīng)濟發(fā)展。第二,綠水系數(shù)和天然徑流系數(shù)之和是一個常數(shù),二者互為消長。水土保持措施的生效是這些措施對于徑流進行調節(jié)的結果,梯田、林草、淤地壩的減沙效益來自減水效應。與這種減水效應相伴隨的是土壤水資源的增加和生態(tài)性、生產性綠水資源的增加。因此,水土保持在減少了下游的河川徑流資源的同時,增加了流域中的土壤水和生態(tài)性、生產性綠水資源。第三,水土保持的水資源效應的含義應包括藍水和綠水在內。水土保持對藍水資源的影響體現(xiàn)為對徑流的影響,是一種異地(off-site)或下游效應,一般表現(xiàn)為減少年徑流。水土保持對綠水資源的影響則體現(xiàn)為對土壤水資源的影響,是一種當?shù)?in-site)效應。綠水系數(shù)的增大意味著當?shù)厣鷳B(tài)耗水量和農作物耗水量的增大,這兩部分耗水量是綠水資源的消耗量,是由土壤水資源轉化而來的。下游河川徑流的減少不一定意味著水資源的凈減少;徑流資源的減少是以土壤水資源(或生態(tài)性、生產性綠水資源)的增加為補償?shù)模m不利于下游河道,但有利于實施水土保持的小尺度、中尺度流域。水土保持的水資源效應包括兩部分,即減少藍水和增加生產性、生態(tài)性綠水,不能只看到前一方面而忽略了后一方面。第四,水土保持的水資源效應該是藍水減少量與生產性、生態(tài)性綠水增加量的代數(shù)和。對于水保措施的水資源效應的評價必須從包含土壤水和生態(tài)生產性綠水資源在內的廣義水資源的概念出發(fā),才能得到全面、正確的認識。在廣義水資源的管理規(guī)劃中如何科學地確定河川徑流資源和生態(tài)生產性綠水資源之間的合理比例,在理論和實踐上都有重要意義。

    6 結論

    (1)本文以大量實測資料為基礎,研究了黃河中游河龍區(qū)間綠水系數(shù)的變化及其影響因素,發(fā)現(xiàn)該地區(qū)綠水系數(shù)呈現(xiàn)增大的趨勢。

    (2)系統(tǒng)地研究了水土保持措施與綠水系數(shù)的關系,發(fā)現(xiàn)水土保持措施的大規(guī)模實施是導致綠水系數(shù)增大的重要因素,梯田面積、造林面積、種草面積和淤地壩造地面積等變量可以解釋綠水系數(shù)變化的39.5%—49.1%。

    (3)首次發(fā)現(xiàn)河龍區(qū)間產沙量與綠水系數(shù)之間呈顯著的負相關關系,綠水系數(shù)的變化可以解釋產沙量變化的53.7%;梯田、壩地面積與綠水系數(shù)和糧食產量之間都存在著顯著的正相關關系,糧食產量與綠水系數(shù)之間也存在顯著的正相關關系;水土保持的實施通過增大梯田和壩地生產性綠水系數(shù),增加了糧食產量。

    (4)綠水系數(shù)具有生態(tài)環(huán)境指示意義,可以作為評價流域生態(tài)環(huán)境變化的指標之一。在年降水可比的情況下綠水系數(shù)的減小意味著生態(tài)環(huán)境環(huán)境質量降低,綠水系數(shù)的增大意味著生態(tài)環(huán)境質量提高。綠水系數(shù)對于廣義的流域水資源評價和水土保持措施水資源效應的確切評價,都有重要的應用意義。

    [1] Falkenmark M. Coping with water scarcity under rapid population growth. Conference of SADC Minister, Pretoria, 1995: 23- 24.

    [2] Falkenmark M. Land-water linkages: A synopsis. Land and Water Bulletin No.1. Rome: FAO, 1995: 15- 17.

    [3] Falkenmark M, Lannerstad M. Consumptive water use to feed humanity-curing a blind spot. Hydrology and Earth System Sciences Discussions, 2004, (1): 7- 40.

    [4] Savenije H H G. Water scarcity indicators: The deception of the numbers. Physics and Chemistry of the Earth (B), 2000, 25(3): 199- 204.

    [5] Ringersma J, Batjes N, Dent D. Green Water: Definitions and Data for Assessment (ISRIC Report). Wageningen, 2003.

    [6] Falkenmark M, Rockstrom J. The new blue and green water paradigm: Breaking new ground for water resources planning and management. Water Resource Planning and Management, American Society of Civil Engineers, 2006, 132(3): 129- 132.

    [7] Postel S L, Daily G C, Ehrlich P R. Human appropriation of renewable fresh water. Science, 1996, 271(5250): 785- 788.

    [8] Rockstrom J, Gordon L. Assessment of green water flows to sustain major biomes of the world: Implications for future ecohydrological landscape management. Physics and Chemistry of the Earth (B), 2001, 26(11- 12): 843- 851.

    [9] Gerten D, Hoff H, Bondeau A, Lucht W, Smith P, Zaehle S. Contemporary "green" water flows: Simulations with a dynamic global vegetation and water balance model. Physics and Chemistry of the Earth, 2005, 30(6-7): 334- 338.

    [10] Jewitt G. Integrating blue and green water flows for water resources management and planning. Physics and Chemistry of the Earth, 2006, 31(15-16): 753- 762.

    [11] 程國棟, 趙文智. 綠水及其研究進展. 地球科學進展, 2006, 21(3): 221- 227.

    [12] 劉昌明, 李云成. “綠水”與節(jié)水: 中國水資源內涵問題討論. 科學對社會的影響, 2006, (1): 16- 20.

    [13] 吳洪濤, 武春友, 郝芳華, 金英學. “綠水”的多角度評估及其管理研究. 中國人口·資源與環(huán)境, 2008, 18(6): 61- 67.

    [14] 王玉娟, 楊勝天, 劉昌明, 戴東, 鄭東海, 曾紅娟. 植被生態(tài)用水結構及綠水資源消耗效用—以黃河三門峽地區(qū)為例. 地理研究, 2009, 28(1): 74- 84.

    [15] 吳洪濤, 武春友, 郝芳華, 金英學. 綠水的多角度評估及其在碧流河上游地區(qū)的應用. 資源科學, 2009, 31(3): 420- 428.

    [16] Xu J X. Effects of climate and land-use change on green-water variations in the middle Yellow River, China. Hydrological Sciences Journal, 2013, 58(1): 1- 12.

    [17] 趙文林. 黃河泥沙. 鄭州: 黃河水利出版社, 1996: 798- 799.

    [18] Demaree G R, Nicolis C. Onset of Sahelian drought viewed as a fluctuation-induced transition. Quarterly Journal of the Royal Meteorological Society, 1990, 116(491): 221-238.

    [19] 中國國家統(tǒng)計局. 中國統(tǒng)計年鑒(2007). 北京: 中國統(tǒng)計出版社, 2008.

    [20] 黃河水利委員會. 黃河年鑒(2006). 鄭州: 黃河水利出版社, 2006.

    [21] 常茂德, 鄭新民, 柏躍勤. 黃土丘陵溝壑區(qū)小流域壩系相對穩(wěn)定及水土資源開發(fā)利用研究. 鄭州: 黃河水利出版社, 2007: 219- 230, 292- 293.

    [22] 許炯心. 黃河中游多沙粗沙區(qū)水土保持減沙的近期趨勢及其成因. 泥沙研究, 2004, (2): 5- 10.

    Increasing trend of green water coefficient in the middle Yellow River basin and the eco-environmental implications

    XU Jiongxin*

    InstituteofGeographicSciencesandNaturalResourcesResearch,KeyLaboratoryforWaterCycleandRelatedLandSurfaceProcesses,ChineseAcademyofSciences,Beijing100101,China

    The concept and theory of green and blue waters are important for water resources assessment and management, especially in semi-arid and sub-humid zones. Most previous studies have been at farm-field scales; to better manage green water at drainage-basin scales, it is important to study how the transformation of precipitation to green water is influenced by climate and human activity. An index of the green water coefficient (Cgw) at the drainage-basin scale is adopted in the present study, defined as the ratio of the annual amount of green water to the annual precipitation within a drainage basin. Based on hydrological and meteorological data from the Yellow River basin and using statistical methods, we analyzed the temporal variation ofCgwin relation to climate change and human activity. During 1950—2011, theCgwover the Hekouzhen-Longmen drainage basin showed an increasing trend. Apart from changes in air temperature and precipitation, the implementation of large-scale soil and water conservation measures is an important factor. Soil and water conservation measures reduce the transformation rate from precipitation to runoff (blue water) and increase the transformation rate from precipitation to green water. The increase in the green water coefficient implies a weakening of runoff on hill slopes and flows in the river, thus reducing soil erosion and sediment yield. On the other hand, the increased green water also indicates higher transpiration by the improved vegetation, and the latter may increase the vegetation′s protection of land surfaces against erosion, also reducing sediment yield. The sediment yield (Qs) over the Hekouzhen-Longmen drainage area has negatively correlated withCgw(R2=0.53), indicating that 53% of the decrease inQscan be explained by the increase inCgw. After the implementation of soil and water conservation measures, the increased terracing land and the land created by check-dams increased the productive green water, increasing grain yield. We found that grain yield was positively correlated withCgwand with the areas of terrace land and the land created by check-dams. The increased proportion of productive green water was an important factor for the increased grain yields. The green water coefficient may be used as an indicator of the environmental quality of a drainage basin. At the same precipitation, an increasedCgwmeans that the environmental quality has improved, and vice versa. According to the variation inCgw, the changes in the environmental quality of the Hekouzhen-Longmen drainage area has been be divided into three stages. In Stage 1 (1950—1969),Cgwshowed a decreasing trend, meaning that the environmental quality was lowered due to serious soil and water losses. In Stage 2,Cgwshowed an increasing trend, meaning that the environmental quality improved due to the large-scale practice of soil and water conservation measures. In Stage 3,Cgwincreased more rapidly, indicating that the environmental quality was further improved due to the large-scale ecological restoration focused on “returning farmland to forests (and grasses)”. Therefore, the application of the green water coefficient is useful for the assessment of generalized water resources at the drainage-basin scale and for a better understanding of the water-resource effect induced by soil and water conservation measures.

    green water; green-water coefficient; climate change; human activity; the middle Yellow River

    國家自然科學基金項目(41371037, 41071016)

    2014- 04- 04;

    日期:2015- 04- 20

    10.5846/stxb201404040646

    *通訊作者Corresponding author.E-mail: xujx@igsnrr.ac.cn

    許炯心.黃河中游綠水系數(shù)變化及其生態(tài)環(huán)境意義.生態(tài)學報,2015,35(22):7298- 7307.

    Xu J X.Increasing trend of green water coefficient in the middle Yellow River basin and the eco-environmental implications.Acta Ecologica Sinica,2015,35(22):7298- 7307.

    猜你喜歡
    藍水綠水梯田
    趕著青山和綠水
    《山區(qū)修梯田》
    藝術品鑒(2019年12期)2020-01-18 08:46:52
    與時俱進樹節(jié)水新觀念
    秀美梯田
    太行山區(qū)藍水綠水沿垂直帶演變規(guī)律及其歸因分析
    梯田之戀
    民族音樂(2017年4期)2017-09-22 03:02:50
    萬綠水 東江請原創(chuàng)歌曲|全國征集
    嶺南音樂(2016年4期)2016-11-21 09:25:24
    為了天藍地綠水清
    綠水繞家園(國畫)
    文藝論壇(2016年3期)2016-03-28 06:54:11
    從江加榜梯田
    www.www免费av| 三级国产精品欧美在线观看| 精品一区二区三区人妻视频| 色吧在线观看| 97超视频在线观看视频| 18禁裸乳无遮挡免费网站照片| 女警被强在线播放| 99久久精品国产亚洲精品| 亚洲avbb在线观看| 夜夜躁狠狠躁天天躁| 欧美乱码精品一区二区三区| 夜夜看夜夜爽夜夜摸| 最好的美女福利视频网| 深夜精品福利| 最近视频中文字幕2019在线8| 男女之事视频高清在线观看| 欧美日韩福利视频一区二区| 亚洲久久久久久中文字幕| 日韩大尺度精品在线看网址| 成人国产综合亚洲| 国产高清激情床上av| 19禁男女啪啪无遮挡网站| 女人十人毛片免费观看3o分钟| 真人一进一出gif抽搐免费| 狂野欧美激情性xxxx| 亚洲在线自拍视频| 免费av不卡在线播放| 国产av在哪里看| 久久精品国产清高在天天线| 最近最新中文字幕大全免费视频| 欧美色欧美亚洲另类二区| 99久久成人亚洲精品观看| 色综合婷婷激情| 黄色日韩在线| 国产男靠女视频免费网站| av欧美777| 在线免费观看的www视频| 我要搜黄色片| 99热6这里只有精品| 久久伊人香网站| 免费av观看视频| 夜夜看夜夜爽夜夜摸| 精品久久久久久,| 一级a爱片免费观看的视频| 亚洲熟妇中文字幕五十中出| 丝袜美腿在线中文| 热99re8久久精品国产| 亚洲国产精品久久男人天堂| 小蜜桃在线观看免费完整版高清| 精品不卡国产一区二区三区| www日本黄色视频网| 亚洲av免费在线观看| 亚洲av电影在线进入| 国产成人影院久久av| 女人被狂操c到高潮| 亚洲在线观看片| 国内精品一区二区在线观看| 国产精品综合久久久久久久免费| 99精品欧美一区二区三区四区| 久久久精品欧美日韩精品| 色精品久久人妻99蜜桃| 精品日产1卡2卡| 欧美一级a爱片免费观看看| 日韩国内少妇激情av| 悠悠久久av| 亚洲成人中文字幕在线播放| 哪里可以看免费的av片| 欧美成人一区二区免费高清观看| 美女高潮的动态| 亚洲在线观看片| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产精品综合久久久久久久免费| 岛国视频午夜一区免费看| а√天堂www在线а√下载| 久久久色成人| 真实男女啪啪啪动态图| 亚洲电影在线观看av| 亚洲av美国av| 非洲黑人性xxxx精品又粗又长| 看黄色毛片网站| 亚洲精品在线观看二区| 国产成人福利小说| 久久久久国内视频| 特大巨黑吊av在线直播| 校园春色视频在线观看| bbb黄色大片| 亚洲欧美一区二区三区黑人| 国产私拍福利视频在线观看| 神马国产精品三级电影在线观看| 国产免费一级a男人的天堂| 91久久精品电影网| 国产黄片美女视频| 成年版毛片免费区| 国产色爽女视频免费观看| aaaaa片日本免费| 日韩亚洲欧美综合| 欧美成人一区二区免费高清观看| 国产精品久久视频播放| 69av精品久久久久久| 日日干狠狠操夜夜爽| 天美传媒精品一区二区| 午夜亚洲福利在线播放| av女优亚洲男人天堂| 母亲3免费完整高清在线观看| 日本a在线网址| 日韩欧美免费精品| 日本 欧美在线| 亚洲成a人片在线一区二区| 怎么达到女性高潮| 精品不卡国产一区二区三区| 一进一出好大好爽视频| 99国产综合亚洲精品| 好男人在线观看高清免费视频| 国内久久婷婷六月综合欲色啪| 国产熟女xx| 91久久精品电影网| 琪琪午夜伦伦电影理论片6080| 母亲3免费完整高清在线观看| e午夜精品久久久久久久| 午夜精品一区二区三区免费看| 内射极品少妇av片p| 日韩欧美 国产精品| 午夜免费观看网址| 午夜免费激情av| 欧洲精品卡2卡3卡4卡5卡区| 国内精品一区二区在线观看| 国产不卡一卡二| 日韩欧美免费精品| 老司机午夜福利在线观看视频| 啦啦啦韩国在线观看视频| 变态另类成人亚洲欧美熟女| 3wmmmm亚洲av在线观看| 日韩欧美三级三区| 日本一二三区视频观看| 国产精品一区二区免费欧美| 久久香蕉精品热| 国产一区在线观看成人免费| 国产高清有码在线观看视频| 亚洲精品456在线播放app | 又粗又爽又猛毛片免费看| 久久久久久久久大av| 精品乱码久久久久久99久播| 综合色av麻豆| a在线观看视频网站| 免费人成在线观看视频色| 国产极品精品免费视频能看的| 午夜a级毛片| 国内毛片毛片毛片毛片毛片| 欧美黑人巨大hd| 性欧美人与动物交配| 国产亚洲精品综合一区在线观看| 男人舔奶头视频| 18禁国产床啪视频网站| 丝袜美腿在线中文| 黄色丝袜av网址大全| 免费无遮挡裸体视频| 亚洲成人中文字幕在线播放| 美女黄网站色视频| 国产亚洲精品久久久com| 精品日产1卡2卡| 两性午夜刺激爽爽歪歪视频在线观看| 日韩成人在线观看一区二区三区| 国产黄色小视频在线观看| 国产97色在线日韩免费| 内射极品少妇av片p| 国产真实伦视频高清在线观看 | 成年女人看的毛片在线观看| 天堂√8在线中文| 免费看十八禁软件| 色尼玛亚洲综合影院| 又黄又爽又免费观看的视频| 国产欧美日韩精品亚洲av| 久久午夜亚洲精品久久| 成人一区二区视频在线观看| 波多野结衣高清无吗| 久久久久久大精品| 亚洲精品456在线播放app | 99热6这里只有精品| 一级毛片高清免费大全| 人人妻,人人澡人人爽秒播| 一个人免费在线观看的高清视频| 国产成人欧美在线观看| h日本视频在线播放| 在线观看一区二区三区| 午夜激情福利司机影院| 欧美日韩综合久久久久久 | 欧美区成人在线视频| 成人特级av手机在线观看| 午夜福利成人在线免费观看| 97碰自拍视频| 久久久久九九精品影院| 搡老熟女国产l中国老女人| 香蕉av资源在线| 精品久久久久久久末码| a级一级毛片免费在线观看| 久久久久久人人人人人| 精品久久久久久久毛片微露脸| 女警被强在线播放| 91九色精品人成在线观看| 婷婷精品国产亚洲av在线| 亚洲黑人精品在线| 一级毛片女人18水好多| 嫩草影视91久久| 精品国产三级普通话版| 麻豆成人午夜福利视频| 免费一级毛片在线播放高清视频| 国产乱人视频| 天天添夜夜摸| 精品日产1卡2卡| 亚洲av二区三区四区| 日韩欧美一区二区三区在线观看| 嫁个100分男人电影在线观看| 熟女少妇亚洲综合色aaa.| 99久久精品一区二区三区| 一个人看的www免费观看视频| 欧美在线黄色| 欧美在线一区亚洲| 在线播放无遮挡| 国产三级在线视频| 岛国在线免费视频观看| 国产成人av激情在线播放| 人妻夜夜爽99麻豆av| 老汉色av国产亚洲站长工具| 熟女人妻精品中文字幕| 欧美日本视频| 美女免费视频网站| 非洲黑人性xxxx精品又粗又长| 日本成人三级电影网站| 国产成人福利小说| 日本撒尿小便嘘嘘汇集6| 精品不卡国产一区二区三区| 一二三四社区在线视频社区8| 黄色片一级片一级黄色片| 在线播放国产精品三级| 亚洲 欧美 日韩 在线 免费| 久久久久久大精品| 99久久综合精品五月天人人| 麻豆久久精品国产亚洲av| 国产一区二区亚洲精品在线观看| 国产真实乱freesex| 国产探花极品一区二区| 最新中文字幕久久久久| 国产亚洲精品久久久com| 可以在线观看毛片的网站| 亚洲av成人不卡在线观看播放网| 国产黄色小视频在线观看| 免费观看精品视频网站| 91av网一区二区| 久久精品91无色码中文字幕| a级一级毛片免费在线观看| 村上凉子中文字幕在线| 亚洲国产精品成人综合色| 99精品在免费线老司机午夜| 此物有八面人人有两片| 精品免费久久久久久久清纯| 九色国产91popny在线| 亚洲美女黄片视频| 精品久久久久久久末码| 日韩大尺度精品在线看网址| 好男人电影高清在线观看| 18禁裸乳无遮挡免费网站照片| 国产探花在线观看一区二区| 精品久久久久久,| 日韩欧美精品免费久久 | 久久人妻av系列| 女同久久另类99精品国产91| 久久久久久久精品吃奶| 国产又黄又爽又无遮挡在线| 一区二区三区高清视频在线| 深爱激情五月婷婷| 狂野欧美白嫩少妇大欣赏| 欧美日韩瑟瑟在线播放| 亚洲国产欧洲综合997久久,| 国内精品久久久久久久电影| 日韩免费av在线播放| 国产精品女同一区二区软件 | 国产私拍福利视频在线观看| 少妇人妻精品综合一区二区 | 级片在线观看| 午夜福利高清视频| 成人午夜高清在线视频| av专区在线播放| 日本 av在线| 国产欧美日韩一区二区精品| 国产伦在线观看视频一区| 欧美bdsm另类| 国产成人aa在线观看| 757午夜福利合集在线观看| 欧美bdsm另类| 国产亚洲精品久久久com| 亚洲国产欧美网| 亚洲专区中文字幕在线| 69人妻影院| 精品久久久久久久人妻蜜臀av| 精华霜和精华液先用哪个| 18美女黄网站色大片免费观看| 人妻丰满熟妇av一区二区三区| 国产精品影院久久| 日韩高清综合在线| 午夜视频国产福利| www日本黄色视频网| 久久九九热精品免费| 亚洲一区二区三区不卡视频| 精品国产美女av久久久久小说| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式 | 久久久色成人| 国产伦在线观看视频一区| 91麻豆精品激情在线观看国产| 女警被强在线播放| 国产精品一区二区三区四区免费观看 | 97碰自拍视频| 嫩草影院入口| 欧美高清成人免费视频www| 国产伦精品一区二区三区视频9 | 人妻丰满熟妇av一区二区三区| 国产精品久久久人人做人人爽| 亚洲人与动物交配视频| 在线观看66精品国产| 丰满乱子伦码专区| 精品人妻一区二区三区麻豆 | 99热精品在线国产| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 99国产精品一区二区蜜桃av| 国产精品野战在线观看| 午夜a级毛片| 国产欧美日韩一区二区精品| 国产免费一级a男人的天堂| 午夜福利视频1000在线观看| 国产熟女xx| 国产精品一及| 首页视频小说图片口味搜索| 全区人妻精品视频| 日本免费a在线| 尤物成人国产欧美一区二区三区| 久久亚洲精品不卡| 亚洲人与动物交配视频| av天堂中文字幕网| 我要搜黄色片| 人妻丰满熟妇av一区二区三区| 色吧在线观看| 国产精品久久久久久久久免 | 国产久久久一区二区三区| 亚洲国产精品sss在线观看| 99久久精品国产亚洲精品| 国产亚洲精品av在线| a级毛片a级免费在线| 校园春色视频在线观看| 91久久精品电影网| 午夜福利在线观看吧| 国产亚洲精品综合一区在线观看| 国产高清视频在线观看网站| 老司机福利观看| 欧美三级亚洲精品| 搡老妇女老女人老熟妇| 免费人成在线观看视频色| 日本 av在线| 国产精品久久电影中文字幕| 久久婷婷人人爽人人干人人爱| 欧美极品一区二区三区四区| 亚洲片人在线观看| 国内久久婷婷六月综合欲色啪| 国产极品精品免费视频能看的| 青草久久国产| 中文亚洲av片在线观看爽| 亚洲精品久久国产高清桃花| 一区二区三区激情视频| 精品一区二区三区视频在线 | 国产中年淑女户外野战色| 国产精品综合久久久久久久免费| 中文字幕精品亚洲无线码一区| 午夜福利高清视频| 老司机午夜十八禁免费视频| 人人妻人人澡欧美一区二区| 99久久综合精品五月天人人| 搡女人真爽免费视频火全软件 | 亚洲成人精品中文字幕电影| 色视频www国产| 成人永久免费在线观看视频| 亚洲av一区综合| 国产精品野战在线观看| 深夜精品福利| svipshipincom国产片| 97碰自拍视频| 国产精品嫩草影院av在线观看 | 中文字幕人妻熟人妻熟丝袜美 | 久久香蕉国产精品| 国产日本99.免费观看| 亚洲成人免费电影在线观看| 丝袜美腿在线中文| 在线观看舔阴道视频| 亚洲18禁久久av| 亚洲精华国产精华精| 国产成人a区在线观看| 在线视频色国产色| 美女cb高潮喷水在线观看| 欧美不卡视频在线免费观看| 99国产综合亚洲精品| 啦啦啦韩国在线观看视频| 他把我摸到了高潮在线观看| 97超级碰碰碰精品色视频在线观看| 午夜福利在线在线| 熟女人妻精品中文字幕| 亚洲中文字幕日韩| 日韩av在线大香蕉| 欧美最新免费一区二区三区 | 一进一出好大好爽视频| 欧美成人免费av一区二区三区| 国产爱豆传媒在线观看| 国产麻豆成人av免费视频| 亚洲 欧美 日韩 在线 免费| 男女午夜视频在线观看| 国产精品久久久久久人妻精品电影| 18+在线观看网站| 他把我摸到了高潮在线观看| 亚洲av成人精品一区久久| 男女做爰动态图高潮gif福利片| 国产精品亚洲av一区麻豆| 在线观看av片永久免费下载| 久久午夜亚洲精品久久| 久久久久国产精品人妻aⅴ院| 国产精品久久久久久人妻精品电影| 国产亚洲精品一区二区www| 九色成人免费人妻av| 日本一本二区三区精品| 男女做爰动态图高潮gif福利片| 97超级碰碰碰精品色视频在线观看| 午夜影院日韩av| 日韩有码中文字幕| 亚洲中文字幕一区二区三区有码在线看| 69av精品久久久久久| 18+在线观看网站| 欧美av亚洲av综合av国产av| 成人性生交大片免费视频hd| 身体一侧抽搐| 久久久久九九精品影院| 小说图片视频综合网站| 久久天躁狠狠躁夜夜2o2o| 欧美日韩瑟瑟在线播放| 一个人免费在线观看电影| 亚洲国产日韩欧美精品在线观看 | 校园春色视频在线观看| 国产精品免费一区二区三区在线| av专区在线播放| 深爱激情五月婷婷| 波野结衣二区三区在线 | 三级毛片av免费| 99视频精品全部免费 在线| 国产高清有码在线观看视频| 久久精品亚洲精品国产色婷小说| 国产精品爽爽va在线观看网站| 国产激情欧美一区二区| 国产欧美日韩精品亚洲av| 欧美国产日韩亚洲一区| 久久精品夜夜夜夜夜久久蜜豆| 国产精品亚洲av一区麻豆| av天堂在线播放| 99在线视频只有这里精品首页| 亚洲电影在线观看av| 久久6这里有精品| 亚洲中文日韩欧美视频| 琪琪午夜伦伦电影理论片6080| 人妻丰满熟妇av一区二区三区| 免费av毛片视频| 欧美另类亚洲清纯唯美| 国产不卡一卡二| 国产视频内射| 亚洲av成人不卡在线观看播放网| 热99re8久久精品国产| 欧美精品啪啪一区二区三区| 日韩欧美免费精品| 久久精品国产亚洲av涩爱 | 亚洲国产日韩欧美精品在线观看 | 亚洲在线观看片| 九九久久精品国产亚洲av麻豆| av片东京热男人的天堂| 欧美又色又爽又黄视频| 51国产日韩欧美| 中文亚洲av片在线观看爽| 最近最新中文字幕大全免费视频| 99久久九九国产精品国产免费| 免费观看的影片在线观看| 在线观看免费视频日本深夜| 国模一区二区三区四区视频| 国产亚洲精品一区二区www| 免费av毛片视频| 夜夜夜夜夜久久久久| 国产亚洲av嫩草精品影院| 嫁个100分男人电影在线观看| 国产激情欧美一区二区| 男女下面进入的视频免费午夜| 好男人在线观看高清免费视频| 国语自产精品视频在线第100页| 亚洲av一区综合| 少妇的丰满在线观看| 欧美性感艳星| 男女之事视频高清在线观看| 最近最新免费中文字幕在线| 久久国产精品人妻蜜桃| 久久久久国产精品人妻aⅴ院| 国产欧美日韩一区二区三| 99国产精品一区二区蜜桃av| 久久草成人影院| 丁香六月欧美| 男女午夜视频在线观看| 啪啪无遮挡十八禁网站| 一卡2卡三卡四卡精品乱码亚洲| 白带黄色成豆腐渣| 国产一区二区在线观看日韩 | 最后的刺客免费高清国语| 亚洲中文字幕日韩| а√天堂www在线а√下载| 亚洲在线自拍视频| 美女高潮喷水抽搐中文字幕| 成人18禁在线播放| 中文字幕人成人乱码亚洲影| 美女高潮的动态| 日韩欧美 国产精品| 午夜福利在线观看吧| 日本熟妇午夜| 黄色视频,在线免费观看| 日日夜夜操网爽| 国产成年人精品一区二区| 国产高清videossex| 美女免费视频网站| 美女高潮的动态| 亚洲精品在线观看二区| 精品熟女少妇八av免费久了| 成年女人看的毛片在线观看| 一区二区三区免费毛片| 亚洲色图av天堂| 免费人成视频x8x8入口观看| 精品久久久久久久人妻蜜臀av| 久久精品国产亚洲av香蕉五月| 又爽又黄无遮挡网站| 成人av一区二区三区在线看| 18禁裸乳无遮挡免费网站照片| 久久久久国产精品人妻aⅴ院| 日本免费a在线| 午夜免费男女啪啪视频观看 | 琪琪午夜伦伦电影理论片6080| 精品国产三级普通话版| 在线观看av片永久免费下载| 两人在一起打扑克的视频| 亚洲自拍偷在线| 国产亚洲欧美98| 国产精品香港三级国产av潘金莲| 亚洲第一电影网av| 色综合亚洲欧美另类图片| a级一级毛片免费在线观看| 国产一区二区三区在线臀色熟女| 国产亚洲精品一区二区www| 国产成人影院久久av| 国产野战对白在线观看| 久久人妻av系列| 国产一区二区三区在线臀色熟女| 全区人妻精品视频| 亚洲欧美激情综合另类| 高清日韩中文字幕在线| 亚洲中文日韩欧美视频| 欧美成人a在线观看| 中文字幕av成人在线电影| 婷婷丁香在线五月| 哪里可以看免费的av片| 婷婷亚洲欧美| 美女cb高潮喷水在线观看| 久久人人精品亚洲av| 国产高清videossex| 乱人视频在线观看| 国产亚洲精品久久久久久毛片| 97超级碰碰碰精品色视频在线观看| 精品国产美女av久久久久小说| 免费av观看视频| 狠狠狠狠99中文字幕| 麻豆久久精品国产亚洲av| 少妇的丰满在线观看| 天堂动漫精品| 又黄又爽又免费观看的视频| 精品熟女少妇八av免费久了| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕久久专区| 久久久久久久精品吃奶| 中文字幕av成人在线电影| 中文字幕人妻熟人妻熟丝袜美 | 成人高潮视频无遮挡免费网站| 国产亚洲av嫩草精品影院| 日本五十路高清| 久久人妻av系列| 两个人视频免费观看高清| 亚洲电影在线观看av| 久久草成人影院| 九色成人免费人妻av| 国产精品国产高清国产av| 女警被强在线播放| 精品国产美女av久久久久小说| 内射极品少妇av片p| 欧美色视频一区免费| 最新中文字幕久久久久| 人妻丰满熟妇av一区二区三区| 久久6这里有精品| 大型黄色视频在线免费观看| 黄色丝袜av网址大全| 欧美一区二区亚洲| 12—13女人毛片做爰片一| 国产精品久久电影中文字幕| 亚洲精品粉嫩美女一区| 免费看光身美女| 操出白浆在线播放| 丰满的人妻完整版| 国产成+人综合+亚洲专区| 久久久国产成人免费| 久久久久久久久大av|