• <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
    從江加榜梯田
    女生性感内裤真人,穿戴方法视频| 波多野结衣av一区二区av| 免费一级毛片在线播放高清视频 | 亚洲精品一二三| 窝窝影院91人妻| 亚洲精华国产精华精| 国产野战对白在线观看| 男男h啪啪无遮挡| 十八禁人妻一区二区| 成人三级黄色视频| 男人舔女人下体高潮全视频| 亚洲av成人一区二区三| 69精品国产乱码久久久| av视频免费观看在线观看| 黄色视频不卡| 妹子高潮喷水视频| 国产aⅴ精品一区二区三区波| 黄色视频,在线免费观看| 自线自在国产av| 人妻丰满熟妇av一区二区三区| 久久亚洲真实| 免费少妇av软件| 亚洲午夜理论影院| 两性夫妻黄色片| 80岁老熟妇乱子伦牲交| 欧美 亚洲 国产 日韩一| 国产精品美女特级片免费视频播放器 | 国产精品乱码一区二三区的特点 | 国产又爽黄色视频| 免费日韩欧美在线观看| 热99re8久久精品国产| 日韩三级视频一区二区三区| 亚洲精品在线美女| 国产精品一区二区免费欧美| 中亚洲国语对白在线视频| 91成人精品电影| 亚洲av成人一区二区三| 国产精品秋霞免费鲁丝片| 欧美日本中文国产一区发布| 国产精品免费视频内射| 大型av网站在线播放| 女人被躁到高潮嗷嗷叫费观| 亚洲欧洲精品一区二区精品久久久| 亚洲黑人精品在线| 午夜免费成人在线视频| 久久久久久久久中文| 一个人观看的视频www高清免费观看 | 国产成人av激情在线播放| 又黄又爽又免费观看的视频| aaaaa片日本免费| 五月开心婷婷网| 少妇粗大呻吟视频| 精品一区二区三区av网在线观看| 高清毛片免费观看视频网站 | 在线国产一区二区在线| 欧美日韩福利视频一区二区| 少妇被粗大的猛进出69影院| 99国产精品一区二区三区| 久久久久久亚洲精品国产蜜桃av| 久久中文字幕人妻熟女| 97碰自拍视频| 色综合婷婷激情| 亚洲成人免费av在线播放| 国产精品一区二区在线不卡| 中国美女看黄片| 久久中文看片网| 国产无遮挡羞羞视频在线观看| 久久午夜综合久久蜜桃| 国产成人av教育| 久久中文字幕人妻熟女| 亚洲成人精品中文字幕电影 | 成人黄色视频免费在线看| 91精品三级在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产激情欧美一区二区| 一进一出抽搐动态| 日本wwww免费看| 成熟少妇高潮喷水视频| 多毛熟女@视频| 搡老乐熟女国产| 91在线观看av| 女人爽到高潮嗷嗷叫在线视频| 国产高清国产精品国产三级| av国产精品久久久久影院| 757午夜福利合集在线观看| 午夜免费成人在线视频| 一个人免费在线观看的高清视频| 欧美乱妇无乱码| 视频区图区小说| 亚洲男人的天堂狠狠| xxxhd国产人妻xxx| 亚洲国产精品合色在线| 好男人电影高清在线观看| 亚洲,欧美精品.| 黑丝袜美女国产一区| 在线观看免费日韩欧美大片| 中出人妻视频一区二区| 无限看片的www在线观看| 久久久久亚洲av毛片大全| 深夜精品福利| 国产成人精品无人区| 首页视频小说图片口味搜索| 91在线观看av| 亚洲精品一区av在线观看| 嫩草影院精品99| 视频区图区小说| 两个人看的免费小视频| 一二三四在线观看免费中文在| 国产精品香港三级国产av潘金莲| 亚洲熟妇中文字幕五十中出 | 国产成人av激情在线播放| 人人澡人人妻人| 一级毛片精品| 亚洲精品成人av观看孕妇| 波多野结衣av一区二区av| 女性生殖器流出的白浆| 老汉色av国产亚洲站长工具| 午夜福利免费观看在线| 51午夜福利影视在线观看| 女警被强在线播放| 97超级碰碰碰精品色视频在线观看| 99re在线观看精品视频| 国产又色又爽无遮挡免费看| 国产成年人精品一区二区 | 久久久久精品国产欧美久久久| 欧美乱色亚洲激情| 久久久久久免费高清国产稀缺| 99re在线观看精品视频| 美女高潮到喷水免费观看| 麻豆av在线久日| www国产在线视频色| 成人精品一区二区免费| 久久婷婷成人综合色麻豆| 他把我摸到了高潮在线观看| 国产成人精品无人区| 在线观看免费午夜福利视频| 91九色精品人成在线观看| 黄片播放在线免费| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久人妻精品电影| 丰满人妻熟妇乱又伦精品不卡| 热re99久久精品国产66热6| 国产精品日韩av在线免费观看 | 香蕉久久夜色| 亚洲精品久久成人aⅴ小说| 嫁个100分男人电影在线观看| 亚洲 欧美一区二区三区| 12—13女人毛片做爰片一| 久久久精品欧美日韩精品| 国产亚洲欧美精品永久| 亚洲人成网站在线播放欧美日韩| 国产精品免费一区二区三区在线| 性少妇av在线| 欧美成人午夜精品| 中文字幕人妻丝袜制服| 亚洲在线自拍视频| 日韩欧美国产一区二区入口| 最近最新免费中文字幕在线| 91精品国产国语对白视频| 无遮挡黄片免费观看| 麻豆av在线久日| 亚洲欧美日韩高清在线视频| 啦啦啦免费观看视频1| 一区二区三区精品91| 99久久99久久久精品蜜桃| 99re在线观看精品视频| bbb黄色大片| 欧美乱码精品一区二区三区| 国产精品免费一区二区三区在线| 国产麻豆69| 欧美精品啪啪一区二区三区| 人人妻人人爽人人添夜夜欢视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品在线美女| 免费高清在线观看日韩| 亚洲男人的天堂狠狠| 老司机靠b影院| 国产精品av久久久久免费| 欧美另类亚洲清纯唯美| 好男人电影高清在线观看| 免费人成视频x8x8入口观看| 亚洲国产欧美网| 一a级毛片在线观看| 欧美乱妇无乱码| av电影中文网址| 欧美午夜高清在线| 两人在一起打扑克的视频| 国产精品一区二区精品视频观看| 亚洲狠狠婷婷综合久久图片| 久久久久精品国产欧美久久久| 伊人久久大香线蕉亚洲五| 大型av网站在线播放| 看免费av毛片| 午夜福利一区二区在线看| 黄色片一级片一级黄色片| 9热在线视频观看99| 啦啦啦 在线观看视频| 精品日产1卡2卡| 国产99白浆流出| 午夜亚洲福利在线播放| 91老司机精品| 在线观看一区二区三区激情| 国产精品美女特级片免费视频播放器 | 欧美日韩瑟瑟在线播放| 90打野战视频偷拍视频| 国产av一区在线观看免费| 自线自在国产av| 久久天躁狠狠躁夜夜2o2o| 丝袜美足系列| 两个人免费观看高清视频| 视频区图区小说| 一级a爱视频在线免费观看| 国产亚洲欧美在线一区二区| 成人18禁高潮啪啪吃奶动态图| 91成人精品电影| 国产色视频综合| 亚洲自偷自拍图片 自拍| 亚洲熟妇中文字幕五十中出 | 亚洲成av片中文字幕在线观看| 欧美丝袜亚洲另类 | 美女高潮喷水抽搐中文字幕| 免费久久久久久久精品成人欧美视频| 久久久久亚洲av毛片大全| www.www免费av| 国产伦一二天堂av在线观看| 女生性感内裤真人,穿戴方法视频| 法律面前人人平等表现在哪些方面| 18禁观看日本| 久久久久九九精品影院| 国产精品国产高清国产av| 国产午夜精品久久久久久| 国产精品久久电影中文字幕| 亚洲国产毛片av蜜桃av| 久久久国产成人免费| 美女国产高潮福利片在线看| 免费观看人在逋| 黄色怎么调成土黄色| 在线观看一区二区三区| 在线av久久热| 亚洲成人久久性| 美女高潮喷水抽搐中文字幕| 久久久久国内视频| 日本三级黄在线观看| 手机成人av网站| 丰满饥渴人妻一区二区三| 国产精品国产av在线观看| 国产av精品麻豆| 国产成人av教育| 99国产精品99久久久久| 乱人伦中国视频| 嫩草影院精品99| 女人被狂操c到高潮| 亚洲专区字幕在线| 美女大奶头视频| 亚洲片人在线观看| 精品一品国产午夜福利视频| 久久精品人人爽人人爽视色| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区在线不卡| 不卡一级毛片| 19禁男女啪啪无遮挡网站| 视频区图区小说| 国产熟女xx| 男女下面进入的视频免费午夜 | 天堂影院成人在线观看| 动漫黄色视频在线观看| 国产黄a三级三级三级人| 日本wwww免费看| 欧美 亚洲 国产 日韩一| 两个人免费观看高清视频| 男人的好看免费观看在线视频 | 欧美日韩视频精品一区| 国产成人av激情在线播放| 久久久久久人人人人人| 身体一侧抽搐| 久久精品亚洲熟妇少妇任你| 国产无遮挡羞羞视频在线观看| 十八禁网站免费在线| 久久香蕉激情| 国产亚洲精品综合一区在线观看 | 婷婷六月久久综合丁香| 如日韩欧美国产精品一区二区三区| 免费观看精品视频网站| av天堂在线播放| 久久99一区二区三区| 国产高清国产精品国产三级| 国产精品综合久久久久久久免费 | 1024香蕉在线观看| 亚洲人成电影免费在线| 日韩中文字幕欧美一区二区| 日韩大码丰满熟妇| 亚洲欧美日韩另类电影网站| 欧美成人性av电影在线观看| 91九色精品人成在线观看| 啦啦啦 在线观看视频| 最新美女视频免费是黄的| 男女下面进入的视频免费午夜 | 日日干狠狠操夜夜爽| 精品卡一卡二卡四卡免费| 国产精品一区二区三区四区久久 | 久久性视频一级片| 日韩中文字幕欧美一区二区| 十八禁人妻一区二区| 97人妻天天添夜夜摸| 中文字幕人妻丝袜制服| 亚洲av片天天在线观看| 免费高清在线观看日韩| 黑人操中国人逼视频| 亚洲熟妇熟女久久| 国产精品一区二区在线不卡| 国产91精品成人一区二区三区| 高清av免费在线| 国产亚洲精品第一综合不卡| 国产97色在线日韩免费| 波多野结衣高清无吗| 国产av又大| 在线观看66精品国产| 欧美在线一区亚洲| 国产乱人伦免费视频| 日日干狠狠操夜夜爽| 午夜福利影视在线免费观看| 叶爱在线成人免费视频播放| 日韩免费高清中文字幕av| 91成年电影在线观看| 免费在线观看日本一区| 国产欧美日韩精品亚洲av| xxx96com| 99国产精品免费福利视频| 中文字幕最新亚洲高清| 国产精品久久久人人做人人爽| 校园春色视频在线观看| 国产不卡一卡二| 国产单亲对白刺激| 在线观看一区二区三区激情| 老司机福利观看| 国产成人系列免费观看| 人妻久久中文字幕网| 丁香欧美五月| 精品久久蜜臀av无| 一a级毛片在线观看| 婷婷丁香在线五月| 在线观看免费视频网站a站| 欧美成人免费av一区二区三区| 国产男靠女视频免费网站| 一进一出抽搐动态| 欧美不卡视频在线免费观看 | 免费在线观看日本一区| 亚洲欧美日韩另类电影网站| 亚洲av美国av| 性少妇av在线| 老司机在亚洲福利影院| 亚洲中文av在线| 婷婷丁香在线五月| 久久婷婷成人综合色麻豆| 最近最新中文字幕大全免费视频| 脱女人内裤的视频| 在线观看一区二区三区激情| 精品免费久久久久久久清纯| 精品国产乱码久久久久久男人| 精品日产1卡2卡| 美女国产高潮福利片在线看| 99精品欧美一区二区三区四区| 国产成人av教育| 精品电影一区二区在线| 黄色片一级片一级黄色片| 大型黄色视频在线免费观看| 免费不卡黄色视频| 久久久国产欧美日韩av| 人妻久久中文字幕网| 国产真人三级小视频在线观看| 久久久国产成人精品二区 | 日韩有码中文字幕| 纯流量卡能插随身wifi吗| 国产欧美日韩一区二区精品| 欧美日韩瑟瑟在线播放| 身体一侧抽搐| 亚洲成人久久性| 中文字幕另类日韩欧美亚洲嫩草| 久9热在线精品视频| 一本综合久久免费| 日韩精品免费视频一区二区三区| 日韩人妻精品一区2区三区| 国产熟女xx| 长腿黑丝高跟| 少妇的丰满在线观看| av免费在线观看网站| 久久精品人人爽人人爽视色| 国产单亲对白刺激| 欧美日韩av久久| 国产免费现黄频在线看| 免费av中文字幕在线| 久久人人精品亚洲av| 成人黄色视频免费在线看| 欧美不卡视频在线免费观看 | 日韩av在线大香蕉| 国产精品野战在线观看 | 国产高清国产精品国产三级| 村上凉子中文字幕在线| 99在线视频只有这里精品首页| 国产精品国产高清国产av| 国产成人精品在线电影| 久久人人爽av亚洲精品天堂| 婷婷六月久久综合丁香| 久久草成人影院| 男人舔女人的私密视频| 又大又爽又粗| 国产精品1区2区在线观看.| 成在线人永久免费视频| 一级毛片女人18水好多| 久久精品国产亚洲av香蕉五月| 亚洲免费av在线视频| 成年版毛片免费区| 国产精品综合久久久久久久免费 | 国产视频一区二区在线看| 亚洲久久久国产精品| 国产精品美女特级片免费视频播放器 | 久久天躁狠狠躁夜夜2o2o| 老司机亚洲免费影院| 91av网站免费观看| 精品国产乱子伦一区二区三区| 两个人免费观看高清视频| 中文字幕色久视频| 欧美日韩亚洲高清精品| av网站免费在线观看视频| 久久国产精品人妻蜜桃| 在线观看免费视频网站a站| av福利片在线| 国产视频一区二区在线看| 欧美成人性av电影在线观看| 亚洲av成人一区二区三| 一级a爱片免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 日韩精品中文字幕看吧| 国产一区二区在线av高清观看| 国产av一区在线观看免费| 久久久久久人人人人人| 日日干狠狠操夜夜爽| 国产亚洲精品综合一区在线观看 | 久久香蕉激情| 又黄又粗又硬又大视频| 国产蜜桃级精品一区二区三区| 极品人妻少妇av视频| 欧美精品亚洲一区二区| 久久久水蜜桃国产精品网| 成人三级做爰电影| 18禁观看日本| 国产成人影院久久av| 国产精品1区2区在线观看.| 国产精品一区二区在线不卡| 黑人巨大精品欧美一区二区蜜桃| 夜夜爽天天搞| 欧美国产精品va在线观看不卡| av在线天堂中文字幕 | 无限看片的www在线观看| 久热这里只有精品99| 精品熟女少妇八av免费久了| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品一二三| 亚洲精品国产一区二区精华液| 好男人电影高清在线观看| 51午夜福利影视在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产又色又爽无遮挡免费看| 欧美精品一区二区免费开放| 久热爱精品视频在线9| 亚洲aⅴ乱码一区二区在线播放 | 亚洲片人在线观看| 午夜精品在线福利| 两性午夜刺激爽爽歪歪视频在线观看 | 99精国产麻豆久久婷婷| 国产成人啪精品午夜网站| 亚洲精品中文字幕在线视频| www.www免费av| 美女国产高潮福利片在线看| 久久久久久大精品| 中文字幕av电影在线播放| 亚洲久久久国产精品| 又紧又爽又黄一区二区| 91精品三级在线观看| 黑人欧美特级aaaaaa片| 国产精品免费一区二区三区在线| 国产97色在线日韩免费| 97碰自拍视频| 精品一区二区三区四区五区乱码| 又黄又爽又免费观看的视频| 亚洲av第一区精品v没综合| 免费在线观看视频国产中文字幕亚洲| 99国产综合亚洲精品| 亚洲av成人一区二区三| 黑丝袜美女国产一区| 国产熟女xx| 国内毛片毛片毛片毛片毛片| 咕卡用的链子| 91av网站免费观看| 不卡av一区二区三区| 久久99一区二区三区| 老司机亚洲免费影院| 51午夜福利影视在线观看| 色婷婷久久久亚洲欧美| 日韩av在线大香蕉| 国产xxxxx性猛交| 很黄的视频免费| 亚洲成人免费电影在线观看| 一个人观看的视频www高清免费观看 | 黄片播放在线免费| 99久久人妻综合| 日本一区二区免费在线视频| 亚洲精品在线观看二区| 在线观看午夜福利视频| 国产亚洲精品一区二区www| 欧美午夜高清在线| 高清在线国产一区| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av在线| 日本撒尿小便嘘嘘汇集6| 国产男靠女视频免费网站| 女人精品久久久久毛片| 好男人电影高清在线观看| 国产成人av教育| 波多野结衣高清无吗| 18禁黄网站禁片午夜丰满| 男人舔女人的私密视频| 欧美日韩福利视频一区二区| 天天躁夜夜躁狠狠躁躁| 极品人妻少妇av视频| av中文乱码字幕在线| 桃色一区二区三区在线观看| 国产野战对白在线观看| 国产97色在线日韩免费| 亚洲自拍偷在线| 久久草成人影院| 一区福利在线观看| 成人18禁在线播放| 精品乱码久久久久久99久播| 精品一区二区三区av网在线观看| 天堂√8在线中文| 亚洲中文日韩欧美视频| 露出奶头的视频| 搡老熟女国产l中国老女人| 欧美日韩亚洲综合一区二区三区_| 中文字幕精品免费在线观看视频| 亚洲av成人av| av在线天堂中文字幕 | 在线观看66精品国产| 日韩高清综合在线| 亚洲欧美精品综合一区二区三区| 国产av精品麻豆| 亚洲九九香蕉| 久久 成人 亚洲| 免费高清在线观看日韩| 免费在线观看视频国产中文字幕亚洲| 日本a在线网址| 一本综合久久免费| 91av网站免费观看| 热99国产精品久久久久久7| 大型av网站在线播放| 一二三四在线观看免费中文在| 久久久久九九精品影院| 欧美亚洲日本最大视频资源| 国产国语露脸激情在线看| 自线自在国产av| 欧美中文综合在线视频| 露出奶头的视频| 亚洲三区欧美一区| 大码成人一级视频| 美女高潮到喷水免费观看| 变态另类成人亚洲欧美熟女 | 日日摸夜夜添夜夜添小说| 亚洲黑人精品在线| 亚洲精品在线观看二区| 欧美久久黑人一区二区| а√天堂www在线а√下载| 久久久久久人人人人人| www.熟女人妻精品国产| 久久中文字幕人妻熟女| 麻豆久久精品国产亚洲av | 国产一区二区三区综合在线观看| 老熟妇乱子伦视频在线观看| 国产欧美日韩综合在线一区二区| 精品国产美女av久久久久小说| 一本综合久久免费| 欧美人与性动交α欧美精品济南到| 国产又色又爽无遮挡免费看| 在线国产一区二区在线| 中文欧美无线码| 在线观看一区二区三区| 9191精品国产免费久久| 久久午夜亚洲精品久久| 精品国产一区二区三区四区第35| 一级毛片高清免费大全| 日韩大码丰满熟妇| 巨乳人妻的诱惑在线观看| 真人做人爱边吃奶动态| 桃红色精品国产亚洲av| 精品一区二区三卡| 麻豆一二三区av精品| 久久久久九九精品影院| 99re在线观看精品视频| 亚洲av成人一区二区三| 制服诱惑二区| 两个人免费观看高清视频| 日韩免费高清中文字幕av| av福利片在线| 18禁国产床啪视频网站| 不卡av一区二区三区| 日韩高清综合在线| 国产精品综合久久久久久久免费 | 成年人黄色毛片网站| 男人舔女人下体高潮全视频| 午夜91福利影院| 真人一进一出gif抽搐免费| 国产av一区二区精品久久| 欧美人与性动交α欧美软件| 国产国语露脸激情在线看|