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

    水足跡分析中國耕地水資源短缺時空格局及驅動機制

    2019-11-08 01:24:28操信春吳夢洋郭相平王衛(wèi)光
    農業(yè)工程學報 2019年18期
    關鍵詞:省區(qū)足跡耕地

    操信春,劉 喆,吳夢洋,郭相平,王衛(wèi)光

    水足跡分析中國耕地水資源短缺時空格局及驅動機制

    操信春1,2,劉 喆2,吳夢洋2,郭相平2,王衛(wèi)光1※

    (1. 河海大學水文水資源與水利工程科學國家重點實驗室,南京 210098;2. 河海大學農業(yè)工程學院,南京 210098)

    為全面評估區(qū)域農業(yè)水資源供需關系,基于水足跡理論構建了耕地水資源短缺指數(arable land water scarcity index,AWSI)。在分析1999-2014年中國AWSI時空分布格局的基礎上,借助偏最小二乘法揭示了AWSI的主控因子。結果顯示:中國AWSI的年均值約為0.413,總體上處于高度水資源壓力狀態(tài),且有隨時間加劇的趨勢;各年份AWSI以華北平原為中心向外遞減式擴散;面臨極高水資源壓力(AWSI>0.800)的省區(qū)均分布在北方地區(qū),長江以南省區(qū)均面臨中度水資源壓力(0.100

    水資源;偏最小二乘回歸;農業(yè);高效用水;藍綠水;灰水足跡;評價指標

    0 引 言

    日益嚴重的水短缺問題已經被視為全球的系統(tǒng)性風險[1]。作為水資源消耗最大的部門,農業(yè)高效用水不僅直接決定食物生產也與生態(tài)安全密切相關,緩解農業(yè)水短缺是促進區(qū)域水資源持續(xù)利用的重要措施[2]。農業(yè)生產過程涉及復雜的水文循環(huán)和面源污染源的遷移過程,消耗水資源同時對環(huán)境造成不利影響。因此,水資源利用的效率和可持續(xù)性是農業(yè)生產與水資源之間關系評價的主要內容。前者關注所占用水資源在農業(yè)生產系統(tǒng)的表現,而后者還涉及區(qū)域水資源承載能力[3-6]?;谒Y源收支對比的水短缺指標構建是進行區(qū)域用水可持續(xù)性評價最常用的手段[7-8]。隨著研究視角的演變,水資源收支的核算手段已由傳統(tǒng)水資源取用發(fā)展到同時衡量水量和水質的水足跡方法[9]。具體到農作物水足跡為作物生長過程中所直接和間接需求的水資源量,包括藍水、綠水及灰水足跡[10-11]。其中藍水和綠水足跡分別為灌溉水和土壤存儲有效降水的消耗數量,灰水足跡則為農業(yè)生產過程所排放污染物的稀釋水量[12]。水足跡理論的核心內涵為廣義水資源與真實水耗用,不僅核算作物生長對藍綠水資源的消耗,還量化了其對水環(huán)境的影響。在不同時空尺度作物水足跡評估的基礎上[13-14],學者開始建立相應區(qū)域水資源短缺評價指標,如水短缺足跡[15-16]、用水效果[17],水質短缺指數[18]以及水足跡稀缺性[19]等。然而這些方法仍大多以藍水為評估對象,未將占主導地位的綠水考慮為水資源短缺評價的參與因素,也沒有很好地將廣義水資源和真實水消耗相統(tǒng)一。本文擬在水足跡理論視角下,建立耕地水資源短缺指數(arable land water scarcity index,AWSI);在核算AWSI時空分布格局的基礎上識別其主要影響因子;并以實現區(qū)域農業(yè)水資源高效可持續(xù)利用目標,展開農業(yè)水足跡調控策略和水資源短缺緩解措施方面的探討,以期在發(fā)展水資源管理相關科學理論的同時,為區(qū)域農業(yè)用水相關宏觀策略的制定提供參考。

    1 研究方法與數據來源

    1.1 耕地水資源短缺指數

    本文基于水足跡理論視角下的廣義水資源與真實水耗用內涵,構建如下耕地水資源短缺評價指標:

    AWSI=CWF/AWR (1)

    式中AWSI為耕地水資源短缺指數,無量綱;CWF為區(qū)域農作物水足跡總量,m3;AWR為區(qū)域農業(yè)廣義水資源量,即可以提供作為農作物生產的藍水和綠水水資源量,m3。

    CWF=CWFb+CWFg+CWFgrey(2)

    AWR=AWRb+AWRg(3)

    式中CWFb、CWFg及CWFgrey分別為作物的藍水、綠水和灰水足跡,m3;AWRg為耕地土壤可儲存的有效降水總量,m3。參照《水足跡評價手冊》,CWFb與CWFg基于作物需水量估算[17],CWFgrey的計算方法為

    CWFgrey=/(max?nat) (4)

    式中為化肥施用量,kg;為淋溶率;max為水體最大容許濃度,kg/m3;nat為自然本底濃度,kg/m3。AWRb與AWRg分別為農業(yè)可利用藍水(常規(guī)水資源)和綠水資源量。由于數據缺乏,按區(qū)域藍水資源供給于不同行業(yè)可能性均等之假設,計算AWRb為

    AWRb=TWR·AWU/TWU (5)

    式中TWR為區(qū)域常規(guī)水資源總量,m3;TWU與AWU分別為總用水量與農業(yè)用水量,m3。

    AWRg=0.1·P(6)

    式中為耕地面積,ha;P為有效降水,mm;P在收集降水量的基礎上以旬為步長,采用下式對地區(qū)和土壤類型均有較廣泛適用性的方法計算[20]:

    基于廣義農業(yè)水資源與水足跡理論構建AWSI,可全面揭示區(qū)域在滿足農業(yè)生產過程對水資源及水環(huán)境需求上的能力,受區(qū)域降雨、灌溉設施、管理水平、耕作措施及耕地率等影響。高AWSI值意味著區(qū)域的農業(yè)生產用水可持續(xù)不足,區(qū)域同時面臨較嚴峻的水資源及水環(huán)境壓力。將AWSI指標進行等級劃分可直觀判斷特定區(qū)域的水資源短缺程度并為進行區(qū)域間對比提供便捷。當前缺乏基于水足跡視角,同時考慮藍、綠和灰水足跡的區(qū)域水資源短缺指標構建及等級劃分。Raskin等[21]提出的以傳統(tǒng)水資源(藍水)為對象區(qū)域水短缺(water shortage index,WSI)的等級劃分方法已被學術界認可和廣泛采用[22-23]。由于AWSI與WSI構建的原理相同,即同為區(qū)域水資源耗用量與水資源可利用量之間比值的核算。因此,參照WSI,AWSI等級閾值及其所對應的水短缺程度如表1所示。

    表1 基于耕地水資源短缺指數的耕地水資源短缺等級劃分

    1.2 空間自相關和偏最小二乘回歸

    全局空間自相關可以描述某種地理屬性在整個區(qū)域(全國)的空間特征,并判定該屬性值在空間上是否表現聚集特征。具體通過計算全局Moran’s的來揭示AWSI的總體空間關聯(lián)或差異特征。在給定顯著水平(0.01)時,若Moran’s值顯著(即-score大于2.58)為正,表示AWSI相當的省區(qū)在空間聚集,否則各省區(qū)與其周邊地區(qū)AWSI總體上差異較大[24]。

    在量化空間格局的基礎上,借助偏最小二乘回歸(partial least-squares regression,PLSR)來分析AWSI的驅動機制。PLSR在進行因變量對自變量的回歸建模中,能夠有效地解決多元回歸分析中的變量多重共線性問題[25]。采用成分提取的方法,用少數幾個綜合指標來代表原有的眾多因素,綜合后的新指標為能綜合反映原有眾多因素的大部分信息的原有因素主成分。PLSR權重是原始變量所定義分數的線性組合,它們可用來描述預測變量和PLSR主成分之間的定量關系。參數變量投影重要性指標(variable importance for the projection,VIP)可以用來表示自變量對因變量預測的重要程度,該指標值大的因素在解釋AWSI時相關性較高。同時,利用回歸系數來量化PLSR模型中每個因素對AWSI影響的方向與強度。選擇11個因素作為PLSR空間分析的響應變量,在進行建模之前將驅動因素進行初步相關性分析,結果見表2。

    表2 AWSI驅動因子的相關矩陣

    注:**表示在0.01水平上顯著相關;*表示在0.05水平上顯著;相對濕度(RH)、平均氣溫(AT)、降水量()、日照時數(SH)、地均化肥(CF)、地均農藥(PP)、地均機械總動力(MP)、灌溉率(IR)、灌溉效率(IE)、糧食面積比例(GA)以及人均GDP(GP),下同。

    Note: **Correlation is significant at 1%, *Correlation is significant at 5%. Abbreviations: relative humidity (RH), average temperature (AT), precipitation (), sunshine hours (SH), chemical fertilizer per area (CF), pesticides per area (PP), power of machinery per area (MP), irrigation rate (IR), irrigation efficiency (IE), proportion of grain area (GA) and per capita GDP (GP), the same as below.

    1.3 數據來源

    本文的研究區(qū)域為中國大陸31個省級行政單元,研究時段為1999-2014年。各省區(qū)歷年水資源總量、用水總量及農業(yè)用水量來自1999-2014年的《中國水資源公報》;各省區(qū)歷年耕地面積、農作物播種面積和產量、化肥施用量、人口、GDP、有效灌溉面積、農業(yè)機械總動力等數據來自《中國統(tǒng)計年鑒》。各省區(qū)主要氣象站點的氣象數據來自中國氣象數據平臺(http://data.cma.cn);作物系數參考文獻[26-27]和水利部灌溉試驗總站收集資料;灌溉效率由水資源公報及各省大型灌區(qū)測算值對有效灌溉面積的加權得到。

    2 結果與分析

    2.1 中國耕地水資源短缺指數

    1999-2014年中國年均AWR與CWF分別約為2 539.8與1 041.7 Gm3,因而形成的AWSI多年平均值為0.413,說明中國耕地上總體處于高度水資源壓力狀態(tài),面臨較嚴重的水短缺。圖1給出了歷年AWR、CWF以及AWSI的全國值,以判別中國農業(yè)廣義水資源收支及水短缺狀況隨時間變化規(guī)律。

    圖1 1999-2014中國農業(yè)廣義水資源收支及耕地水資源短缺指數

    圖1顯示,全國可用于農業(yè)生產的廣義水資源量AWR在年際間呈小幅波動,無明顯的變化趨勢,各值介于2 200~2 900 Gm3之間。這主要與歷年氣候條件、耕地面積總量、取水總量及用水組成等因素總體上的相對穩(wěn)定有關。研究時段內,全國耕地面積始終維持在1.27×108hm2左右;水資源總量和農業(yè)用水量也分別在2.7×1012m3和3 700×108m3左右小幅變化。雖然初始出現回落現象,CWF總體上呈隨時間穩(wěn)步增長的趨勢。作物耗水和水分生產力變化不大的前提下,農業(yè)生產規(guī)模的擴張和種植結構的改變是CWF穩(wěn)步增長的主要驅動力。15 a間,農作物總播種面積由1.56×108hm2擴大到1.65×108hm2,與此同時高耗水的糧食作物產量由約5.08×108t增長到6.07×108t左右。AWR與CWF的數值及變化趨勢共同決定AWSI的表現:前4 a穩(wěn)定在0.320左右并面臨中度水資源壓力,2003年之后持續(xù)上升,2006年超過0.400落入高度水資源壓力區(qū)間并于2011年達到研究時段的峰值0.534;受全國水資源總量波動導致AWRblue年際變化加大的影響,近5 a AWSI以0.476為均值出現較大幅度的波動,不僅說明中國耕地水短缺加劇,也揭示其不確定性在增強。

    2.2 AWSI的時空分布

    圖2為AWSI歷年全局Moran’s及檢驗值Z-Score。由圖可以看出,Moran’s在所有的年份均大于0,同時相應的Z-Score值大于顯著水平為0.01時的臨界值(2.58)。說明中國耕地水資源短缺指數表現為正自相關關系,相似值之間在空間上呈相互聚集的總體格局。AWSI的決定因素涉及自然氣候、水資源供需以及農業(yè)生產特征等與作物-水資源關系相聯(lián)系的諸多方面。而以上各因素在時間上表現無突變現象的同時在地理上基本上均表現出臨近區(qū)域的類似性[28-29],這應該是歷年耕地水資源短缺指數在空間上呈聚集現象的主要原因。

    圖2 1999-2014中國AWSI全局Moran’s I指數及其檢驗

    根據AWSI的空間分布圖可以分別判別出較大值和較低值所具體聚集的區(qū)域。觀察各年份空間分布狀況發(fā)現,雖然數值存在差異,但是各省AWSI不同年份在全國相對大小基本穩(wěn)定,空間上的聚集特征在時間上未發(fā)現較大變化。因此給出1999、2004、2009及2014年4個年份的AWSI的空間分布圖,以觀察其在空間上的分布和聚集特征,如圖3所示。同時,計算出各省區(qū)AWSI的年際平均值、極差、標準差、年均變化率及變異系數等主要統(tǒng)計量,為明晰省區(qū)耕地水資源短缺及其隨時間變化情況,結果列于表3。

    圖3顯示,不同年份AWSI的總體空間分布格局與全局空間自相關的計算結果相符,表現為明顯而穩(wěn)定的空間聚集現象。AWSI較大省區(qū)聚集于黃淮海平原,并隨著時間推移有向西北和東北方向延伸的趨勢;青藏地區(qū)及其周邊的西南省區(qū)是AWSI較低的區(qū)域;長江右岸各省的耕地水短缺程度位于全國中間位置,也呈現隨時間惡化的趨勢。AWSI在局部空間的分布狀況直接由水資源稟賦與農業(yè)生產規(guī)模決定,這些條件在地理上聚集是形成各年份相似且局部聚集空間格局的原因。結合圖3與表2來看,山東省1999-2014年間的AWSI為2.056,居全國之首;河北、寧夏、北京、河南、天津及上海的AWSI值也明顯高于極高水資源壓力閾值的0.800;以上省區(qū)為中國耕地農業(yè)生產水資源利用最不可持續(xù)的區(qū)域。江蘇與遼寧也屬于極高水資源壓力的范疇,農業(yè)生產用水過程面臨較大挑戰(zhàn)。10個省區(qū)的AWSI介于0.400~0.800之間,面臨高度水資源壓力,這些省區(qū)分布在極高水資源壓力區(qū)域的邊緣;表現為低度水資源壓力和中度水資源壓力(0.100

    圖3 典型年份耕地水資源短缺指數AWSI的空間分布

    表3 各省區(qū)AWSI在1999-2014年間AWSI的主要統(tǒng)計量

    不同年份間的變化方面,極差與標準差的空間分布與耕地水資源短缺指數相似,基本上表現為AWSI值越大的省區(qū)年際間絕對變化越大,山東、河北、寧夏、北京與河南的極差超過1.000,而青海與西藏低于0.100;西藏標準差幾乎為0,與年均值、極差共同揭示該區(qū)不僅AWSI低,且在年際間維持穩(wěn)定;其他省區(qū)間標準差介于0.008(青海)~0.635(寧夏)之間,差異較大;省區(qū)間的AWSI值及其年際間絕對變化較大,然而除西藏與寧夏外,變異系數基本上介于0.100~0.300之間,年際間相對變化較為接近。西藏與寧夏的AWSI表現不同于其他大部分省區(qū)主要是由其水資源和農業(yè)生產特征決定的。雖然耕地上的有效降水量不足400 mm,但由于河川徑流量大,西藏的年均AWR達到了400 Gm3,為所有省區(qū)之最;然而由于農業(yè)生產規(guī)模小,西藏年均CWF僅為1.55 Gm3,全國最??;歷年CWF占AWR的比例均極小,這也是AWSI小且波動微弱的原因。寧夏的降水稀少與農業(yè)水資源缺口大而造成嚴重的耕地水資源短缺,同時,AWR的年際變化和CWF的大幅增長共同決定了AWSI的年際間的變異系數最大。表3還顯示,所有省區(qū)AWSI整體上均呈隨時間增大的趨勢,說明近15 a中國耕地水短缺全面加劇,農業(yè)用水效率提升是當前面臨重要需求。由農業(yè)生產規(guī)模擴張導致的農作物水足跡增大(圖1)是AWSI隨時間增長的主要原動力。大部分省區(qū)AWSI的年均增長率不超過10.0%,新疆的增速達到了26.9%,明顯高于其他省區(qū)。新疆AWR保持在100 Gm3,無明顯變化趨勢,而CWF由21.51 Gm3一路上升到2014年的75.75 Gm3。這也使得新疆由低度~中度水短資源壓力地區(qū)轉換為高度~極高水資源壓力地區(qū)。

    2.3 AWSI驅動要素解析

    利用PLSR對AWSI的主控因子判別的分析結果見表4、圖4及圖5。表4所示,PLSR模型提取出的2個偏最小二乘回歸主成分分別解釋ASWI空間變異的58.9%和13.8%,而增加更多成分到ASWI偏最小二乘回歸模型中,未能提高對ASWI空間變異的貢獻率。圖4顯示,第1成分中農業(yè)機械總動力(MP)對ASWI大小有正作用,降水量()對ASWI大小呈負作用。第2成分中糧食面積比例(GA)對ASWI大小有正向影響,日照時數(SH)對ASWI有負向影響,表現為不顯著。

    表4 ASWI的PLSR模型

    一般情況下,VIP值在1之上的變量為因變量的重要影響因素。AWSI的PLSR模型中(圖5)3個變量的VIP值大于1,分別為農業(yè)機械總動力(MP,VIP=1.761,回歸系數=0.457),降水量(,VIP=1.192,回歸系數=?0.312),人均GDP(GP,VIP=1.170,回歸系數=0.030)。其中農業(yè)機械總動力和人均GDP的增大會導致ASWI增大,降水量則和ASWI呈顯著負相關。與此相對,VIP較小的因素包括相對濕度(RH,VIP=0.729)、平均氣溫(AT,VIP=0.321)、地均化肥(CF,VIP=0.822)及地均農藥(PP,VIP=0.458)。

    圖4 AWSI偏最小二乘模型第一、二成分權重

    圖5 ASWI偏最小二乘模型中各因子變量投影重要性指標VIP及回歸系數

    降水量為與農作物類型和水資源稟賦的空間分布有直接聯(lián)系的氣象因子。中國的降水量大體呈現東南向西北遞減的空間格局,而日照時數正好基本相反。降水量與區(qū)域農業(yè)廣義水資源密不可分,其直接決定耕地可利用綠水資源的同時很大程度上影響藍水資源的豐欠。農業(yè)機械總動力、糧食面積比重以及人均GDP均為能綜合反映區(qū)域農業(yè)生產條件與用水水平的社會經濟因子。農業(yè)機械總動力最主要部分體現于農田耕作機械和灌排機械,二者分別對區(qū)域農業(yè)用水狀況緊密相連的同時影響區(qū)域的耕地水短缺狀況。耕作機械動力投入量較大的區(qū)域大多土地肥沃平整、農業(yè)基礎設施完善,如位于水資源相對匱缺地區(qū)的華北平原;由于主要作物生育期降水有限,這些地區(qū)農業(yè)生產對灌溉水及排灌機械的需求較大,如華北各省市農業(yè)機械總動力高達10 kW/hm2,而長江流域糧食主產區(qū)僅為3 kW/hm2左右;農業(yè)機械投入和農業(yè)生產效率的提高會促使灌溉面積和農業(yè)生產規(guī)模的擴大,從而形成更大的農作物水足跡和更嚴峻的耕地水短缺。由于糧食作物耗水量和單位面積水足跡值均高于其他作物,糧食面積比重較大省區(qū)農作物水足跡總量也趨于較大,進而導致耕地水資源短缺指數的增大。人均GDP衡量區(qū)域經濟發(fā)展程度也側面反映行業(yè)用水結構并影響耕地的水短缺狀況:由于工業(yè)和服務業(yè)發(fā)展水平和水資源需求均較高,人均GDP越大的省區(qū)農業(yè)用水比例相對低,這會造成當地分配于農業(yè)領域的可利用水資源量較小,進而增大耕地水資源短缺指數。氣象要素相對濕度(RH)與平均氣溫(AT)雖影響作物耗水過程和藍綠水足跡,但是與作物灰水足跡以及AWR無直接關聯(lián),這可能是二者VIP值較低的原因。農業(yè)生產投入地均化肥(CH)及地均農藥(PP)在區(qū)域間的差別不如其他社會經濟因子明顯,且當前灰水足跡估算過程中難以考慮區(qū)域間化肥有效利用程度的差別。這是他們在空間上對AWSI的影響相對較小主要原因。此外,土壤特性對于作物生產水足跡相關,從而可能成為AWSI的驅動要素。因數據來源缺陷,該要素在PLSR分析未予考慮,這是未來區(qū)域農業(yè)水資源短缺驅動機制分析中需要注意的問題。

    3 結 論

    1)水足跡視角下,耕地水資源短缺指數(arable land water scarcity index,AWSI)能夠全面量化區(qū)域廣義水資源可持續(xù)利用程度。15 a間全國AWSI為0.413,農業(yè)生產領域一直處于隨時間加劇趨勢的高度水資源壓力狀態(tài),將持續(xù)面臨嚴峻的資源和環(huán)境兩方面導致的水短缺。

    2)耕地水資源短缺指數同時存在時間和空間上的不平衡性,華北平原是耕地水短缺問題最為嚴重的區(qū)域,長江以南和東北的糧食主產區(qū)也面臨中度以上水資源壓力,西北旱區(qū)的耕地水短缺加劇趨勢最為明顯。

    3)降水是對AWSI影響最為明顯的氣象因子,而農業(yè)機械總動力、糧食面積比例以及人均GDP對AWSI表現出了較明顯的影響。在經濟發(fā)展和農業(yè)生產水平較高的糧食主產區(qū)應基于全區(qū)域的農產品產出保障和水足跡總量控制來進行農業(yè)水資源管理策略制定,以實現水資源的可持續(xù)利用。

    [1] 姜秋香,周智美,王子龍,等. 基于水土資源耦合的水資源短缺風險評價及優(yōu)化[J]. 農業(yè)工程學報,2017,33(12):136-143.

    Jiang Qiuxiang, Zhou Zhimei, Wang Zilong, et al. Risk assessment and optimization of water resources shortage based on water and land resources coupling[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 136-143. (in Chinese with English abstract)

    [2] Cao X, Wu M, Zheng Y, et al. Can China achieve food security through the development of irrigation?[J]. Regional Environmental Change, 2018, 18(2): 465-475.

    [3] Velasco M, Aznar S, Belmonte U, et al. Advances in water use efficiency in agriculture: A bibliometric analysis[J]. Water, 2018, 10(4): 377.

    [4] Wang Y, Zhang Y, Zhang R, et al. Reduced irrigation increases the water use efficiency and productivity of winter wheat-summer maize rotation on the North China Plain[J]. Science of the Total Environment, 2018, 618: 112-120.

    [5] 李慧,周維博,莊妍,等. 延安市農業(yè)水土資源匹配及承載力[J]. 農業(yè)工程學報,2016,32(5):156-162.

    Li Hui, Zhou Weibo, Zhuang Yan, et al. Agricultural water and soil resources matching patterns and carrying capacity in Yan’an City[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(5): 156-162. (in Chinese with English abstract)

    [6] Pedro M, Solera A, Ferrer J, et al. A review of water scarcity and drought indexes in water resources planning and management[J]. Journal of Hydrology, 2015, 527: 482-493.

    [7] Li J, Liu Z, He C, et al. Water shortages raised a legitimate concern over the sustainable development of the drylands of northern China: Evidence from the water stress index[J]. Science of The Total Environment, 2017, 590: 739-750.

    [8] Sun S, Wang Y, Engel B, et al. Effects of virtual water flow on regional water resource stress: A case study of grain in China[J]. Science of the Total Environment, 2016, 550: 871-879.

    [9] Hoekstra A. Water footprint assessment: Evolvement of a new research field[J]. Water Resources Management, 2017, 31(10): 3061-3081.

    [10] Xinchun C, Mengyang W, Rui S, et al. Water footprint assessment for crop production based on field measurements: A case study of irrigated paddy rice in East China[J]. Science of the Total Environment, 2018, 610: 84-93.

    [11] Chu Y, Shen Y, Yuan Z. Water footprint of crop production for different crop structures in the Hebei southern plain, North China[J]. Hydrology and Earth System Sciences, 2017, 21(6): 3061-3069.

    [12] Zhang Y, Huang K, Yu Y, et al. Mapping of water footprint research: A bibliometric analysis during 2006-2015[J]. Journal of Cleaner Production, 2017, 149: 70-79.

    [13] Chukalla A, Krol M, Hoekstra A. Grey water footprint reduction in irrigated crop production: Effect of nitrogen application rate, nitrogen form, tillage practice and irrigation strategy[J]. Hydrology and Earth System Sciences, 2018, 22(6): 3245-3259.

    [14] Dambrosio E, Degirolamo A M, Rulli M C. Assessing sustainability of agriculture through water footprint analysis and in-stream monitoring activities[J]. Journal of Cleaner Production, 2018, 200: 454-470.

    [15] Zhuo L, Mekonnen M, Hoekstra A, et al. Inter-and intra-annual variation of water footprint of crops and blue water scarcity in the Yellow River basin (1961-2009)[J]. Advances in Water Resources, 2016, 87: 29-41.

    [16] Veettil A, Mishra A. Water security assessment using blue and green water footprint concepts[J]. Journal of Hydrology, 2016, 542: 589-602.

    [17] 操信春,任杰,吳夢洋,等. 基于水足跡的中國農業(yè)用水效果評價[J]. 農業(yè)工程學報,2018,34(5):1-8.

    Cao Xinchun, Ren Jie, Wu Mengyang, et al. Assessing agricultural water use effect of China based on water footprint framework[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(5): 1-8. (in Chinese with English abstract)

    [18] Liu J, Liu Q, Yang H. Assessing water scarcity by simultaneously considering environmental flow requirements, water quantity, and water quality[J]. Ecological Indicators, 2016, 60: 434-441.

    [19] Cao X, Huang X, Huang H, et al. Changes and driving mechanism of water footprint scarcity in crop production: A study of Jiangsu Province, China[J]. Ecological Indicators, 2018, 95: 444-454.

    [20] 操信春,邵光成,王小軍,等. 中國農業(yè)廣義水資源利用系數及時空格局分析[J]. 水科學進展,2017,28(1):14-21.

    Cao Xinchun, Shao Guangcheng, Wang Xiaojun, et al. Generalized water efficiency and strategic implications for food security and water management: A case study of grain production in China[J]. Advances in Water Science, 2017, 28(1): 14-21. (in Chinese with English abstract)

    [21] Raskin P, Gleick P, Kirshen P, et al. Water Futures: Assessment of long-range patterns and prospects[R]. Stockholm, Sweden: Stockholm Environment Institute, 1997.

    [22] Munia H, Guillaume J, Mirumachi N, et al. Water stress in global transboundary river basins: Significance of upstream water use on downstream stress[J]. Environmental Research Letters, 2016, 1(1): 014002.

    [23] Xinchun C, Mengyang W, Xiangping G, et al. Assessing water scarcity in agricultural production system based on the generalized water resources and water footprint framework[J]. Science of the Total Environment, 2017, 609: 587-597.

    [24] Xia J, Wang L, Yu J, et al. Impact of environmental factors on water quality at multiple spatial scales and its spatial variation in Huai River Basin, China[J]. Science China Earth Sciences, 2018, 61(1): 82-92.

    [25] Ye X, Wu Z, Wang Z, et al. Seasonal prediction of the Yangtze River runoff using a partial least squares regression Model[J]. Atmosphere-Ocean, 2018, 56(2): 117-128.

    [26] 陳玉民,郭國雙,王廣興,等. 中國主要作物需水量與灌溉[M]. 北京:中國水利電力出版社,1995.

    [27] 段愛旺,孫景生,劉鈺,等. 北方地區(qū)主要農作物灌溉用水定額[M]. 北京:中國農業(yè)科學技術出版社,2004.

    [28] 孫才志,陳麗新,劉玉玉. 中國農作物綠水占用指數估算及時空差異分析[J]. 水科學進展,2010,21(5):637-643.

    Sun Caizhi, Chen Lixin, Liu Yuyu. Spatial and temporal variation of crops green water occupancy index in China[J]. Advances in Water Science, 2010, 21(5): 637-643. (in Chinese with English abstract)

    [29] 操信春,吳普特,王玉寶,等. 中國灌區(qū)水分生產率及其時空差異分析[J]. 農業(yè)工程學報,2012,28(13):1-7.

    Cao Xinchun, Wu Pute, Wang Yubao, et al. Analysis on temporal and spatial differences of water productivity in irrigation districts in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(13): 1-7. (in Chinese with English abstract)

    Temporal-spatial distribution and driving mechanism of arable land water scarcity index in China from water footprint perspective

    Cao Xinchun1,2, Liu Zhe2, Wu Mengyang2, Guo Xiangping2, Wang Weiguang1※

    (1.,210098,; 2.,,210098,)

    Efficient water use in agriculture production system is widely accepted as an important foundation of regional water resources management, water shortage alleviation and environmental sustainability. The arable land water scarcity index (AWSI) to describe relationship between crop production and potentially water resources was established based on water footprint framework in current study. AWSI was defined as the ratio of total water footprint in regional crops cultivation to available agricultural water resources, including blue and green water. AWSI in 31 provinces, municipalities and autonomous regions of China from 1999 to 2014 was calculated. Then, the spatial-temporal pattern and driving mechanism in the observed period were explored with the help of the methods of spatial autocorrelation analysis. A total of 10 potential factors such as relative humidity (RH), average temperature (AT), precipitation (), sunshine hours (SH), chemical fertilizer per area (CF), pesticides per area (PP), power of machinery per area (MP), irrigation rate (IR), irrigation efficiency (IE), proportion of grain area (GA) and per capita GDP (GP) were selected in driving mechanism assessment. Given that the high co-dependence of these potential factors, partial least squares regression (PLSR) was used to elucidate the linkages between the ASWI and the selected factors. The results showed that, annual value of AWSI in China was estimated about 0.413, and the country faced high water stress during the studied period; AWSI in almost all of the provinces, municipalities and autonomous regions increased over time, indicating that water scarcity in agricultural production system of China was intensifying. Spatial autocorrelation analysis showed that the global Moran’s I was higher than 0 in all the calculated years, implying provinces, municipalities and autonomous regions with similar AWSI presented an obvious aggregation characteristic in agriculture production of China. Provinces with high AWSI was in the North China Plain and all the regions facing extremely high water stress (AWSI>0.800) were distributed in the north of China; most of the provinces located in south of the Yangtze River were classified as moderate water stress (0.100

    water resources; partial least-squares regression (PLSR); agriculture; efficient water use; blue and green water; grey water footprint; evaluating indicator

    10.11975/j.issn.1002-6819.2019.18.012

    TV213.9; S279

    A

    1002-6819(2019)-18-0094-07

    操信春,劉 喆,吳夢洋,郭相平,王衛(wèi)光. 水足跡分析中國耕地水資源短缺時空格局及驅動機制[J]. 農業(yè)工程學報,2019,35(18):94-100.doi:10.11975/j.issn.1002-6819.2019.18.012 http://www.tcsae.org

    Cao Xinchun, Liu Zhe, Wu Mengyang, Guo Xiangping, Wang Weiguang. Temporal-spatial distribution and driving mechanism of arable land water scarcity index in China from water footprint perspective[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(18): 94-100. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.18.012 http://www.tcsae.org

    2019-03-12

    2019-08-10

    國家自然科學基金(51979074、51609065);國家重點研發(fā)計劃(2018YFF0215702);江蘇省社會科學基金(17GLC013);中國博士后科學基金資助項目(2018T110436,2017M611681)

    操信春,博士,教授,主要從事農業(yè)水資源高效利用研究。Email:caoxinchun@hhu.edu.cn。

    王衛(wèi)光,博士,教授,主要從事水資源演變與高效利用研究。Email:wangweiguang@hhu.edu.cn。

    中國農業(yè)工程學會會員:操信春(E040000462A)

    猜你喜歡
    省區(qū)足跡耕地
    自然資源部:加強黑土耕地保護
    我國將加快制定耕地保護法
    我校成功舉辦第四屆沿黃九省區(qū)黃河論壇
    保護耕地
    北京測繪(2021年12期)2022-01-22 03:33:36
    新增200億元列入耕地地力保護補貼支出
    成長足跡
    足跡
    心聲歌刊(2018年6期)2018-01-24 00:56:12
    真抓實干,為這26省區(qū)點贊!
    足跡
    社會與公益(2016年2期)2016-04-13 02:49:06
    我國沿海省區(qū)海洋產業(yè)集聚水平比較研究
    欧美一级a爱片免费观看看 | 久久久久久久久免费视频了| 日本撒尿小便嘘嘘汇集6| 50天的宝宝边吃奶边哭怎么回事| 不卡一级毛片| АⅤ资源中文在线天堂| 99香蕉大伊视频| 亚洲天堂国产精品一区在线| 久久国产精品影院| 日韩精品免费视频一区二区三区| 国产成人精品久久二区二区91| 精品熟女少妇八av免费久了| 黄色a级毛片大全视频| 精品不卡国产一区二区三区| 黄片小视频在线播放| 亚洲五月天丁香| 中文字幕色久视频| www日本在线高清视频| 神马国产精品三级电影在线观看 | 国产成人系列免费观看| 欧美成人性av电影在线观看| 色av中文字幕| 不卡av一区二区三区| 男人舔女人的私密视频| 亚洲av电影在线进入| 麻豆国产av国片精品| 国产精品一区二区免费欧美| 国产高清视频在线播放一区| 亚洲无线在线观看| 国产成人一区二区三区免费视频网站| 国内久久婷婷六月综合欲色啪| 久久久久国产一级毛片高清牌| 亚洲成人久久性| 久久精品国产综合久久久| 啦啦啦观看免费观看视频高清 | 亚洲国产欧美日韩在线播放| 免费无遮挡裸体视频| 日韩精品中文字幕看吧| 99国产精品99久久久久| 亚洲色图综合在线观看| 日日夜夜操网爽| 色综合婷婷激情| 好男人电影高清在线观看| 久久香蕉国产精品| 亚洲免费av在线视频| 波多野结衣巨乳人妻| 天天躁夜夜躁狠狠躁躁| 亚洲狠狠婷婷综合久久图片| 中出人妻视频一区二区| 天天一区二区日本电影三级 | 国产亚洲av嫩草精品影院| 级片在线观看| 亚洲中文字幕一区二区三区有码在线看 | 男人舔女人下体高潮全视频| 在线观看午夜福利视频| 午夜久久久在线观看| 亚洲午夜精品一区,二区,三区| 国产一区在线观看成人免费| 久久中文字幕人妻熟女| 精品午夜福利视频在线观看一区| 麻豆成人av在线观看| 亚洲视频免费观看视频| 香蕉国产在线看| av天堂久久9| 侵犯人妻中文字幕一二三四区| 国产主播在线观看一区二区| 成人手机av| 精品一区二区三区av网在线观看| 天天躁夜夜躁狠狠躁躁| 亚洲九九香蕉| 女人精品久久久久毛片| 精品国产一区二区久久| 久久精品国产亚洲av香蕉五月| av天堂在线播放| av欧美777| 日本欧美视频一区| 欧美在线黄色| 成人免费观看视频高清| 久久人妻福利社区极品人妻图片| 免费高清在线观看日韩| 午夜亚洲福利在线播放| 青草久久国产| 99久久综合精品五月天人人| 亚洲av片天天在线观看| 国语自产精品视频在线第100页| 亚洲一区二区三区色噜噜| 精品电影一区二区在线| 国产精品久久久久久亚洲av鲁大| 国产真人三级小视频在线观看| 午夜亚洲福利在线播放| 成人亚洲精品av一区二区| 亚洲av五月六月丁香网| 他把我摸到了高潮在线观看| 99香蕉大伊视频| 一级毛片高清免费大全| 欧美av亚洲av综合av国产av| 最近最新免费中文字幕在线| 国产激情欧美一区二区| 搡老岳熟女国产| av网站免费在线观看视频| 午夜福利成人在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩精品亚洲av| 淫秽高清视频在线观看| 欧美另类亚洲清纯唯美| 亚洲成人久久性| 欧美日韩中文字幕国产精品一区二区三区 | 青草久久国产| 热re99久久国产66热| 国产三级在线视频| 亚洲精品久久国产高清桃花| 99国产精品99久久久久| 国产三级在线视频| 日韩高清综合在线| 精品卡一卡二卡四卡免费| 女人高潮潮喷娇喘18禁视频| 国产伦人伦偷精品视频| 黑人欧美特级aaaaaa片| 九色亚洲精品在线播放| 又紧又爽又黄一区二区| 午夜福利视频1000在线观看 | 国产人伦9x9x在线观看| 国产精品亚洲av一区麻豆| 在线观看日韩欧美| 老司机午夜十八禁免费视频| 成人亚洲精品av一区二区| 熟女少妇亚洲综合色aaa.| 成人18禁高潮啪啪吃奶动态图| 久久人妻av系列| 美国免费a级毛片| 久久香蕉精品热| 12—13女人毛片做爰片一| av在线播放免费不卡| 久久精品aⅴ一区二区三区四区| 在线天堂中文资源库| 大陆偷拍与自拍| www.www免费av| 国产激情欧美一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费激情av| 日日摸夜夜添夜夜添小说| 午夜福利高清视频| 国产不卡一卡二| 国产91精品成人一区二区三区| 韩国av一区二区三区四区| 成人国产综合亚洲| 亚洲中文字幕一区二区三区有码在线看 | 日本撒尿小便嘘嘘汇集6| 欧美一级毛片孕妇| 在线播放国产精品三级| 国产亚洲精品av在线| 国产精品乱码一区二三区的特点 | 中亚洲国语对白在线视频| 免费在线观看黄色视频的| 可以免费在线观看a视频的电影网站| av网站免费在线观看视频| 日本撒尿小便嘘嘘汇集6| 一级毛片女人18水好多| 亚洲午夜精品一区,二区,三区| 最好的美女福利视频网| 亚洲第一av免费看| av在线天堂中文字幕| 69精品国产乱码久久久| 久久热在线av| 一级a爱视频在线免费观看| 搡老岳熟女国产| 亚洲狠狠婷婷综合久久图片| 日本一区二区免费在线视频| 欧美一级a爱片免费观看看 | 国产91精品成人一区二区三区| АⅤ资源中文在线天堂| 男女下面插进去视频免费观看| 如日韩欧美国产精品一区二区三区| 午夜成年电影在线免费观看| 欧美最黄视频在线播放免费| 亚洲一区中文字幕在线| 男人舔女人的私密视频| 国产精品,欧美在线| 亚洲九九香蕉| 欧美人与性动交α欧美精品济南到| 久久久久久亚洲精品国产蜜桃av| 99国产精品99久久久久| 丰满的人妻完整版| 日本vs欧美在线观看视频| 国产精品一区二区精品视频观看| 婷婷六月久久综合丁香| 国产精品美女特级片免费视频播放器 | 88av欧美| 国产成人欧美在线观看| 日本免费一区二区三区高清不卡 | 国产一区二区三区视频了| 国产精品av久久久久免费| 国产高清videossex| 黑人巨大精品欧美一区二区蜜桃| 国产精华一区二区三区| 国产伦一二天堂av在线观看| АⅤ资源中文在线天堂| 欧美另类亚洲清纯唯美| ponron亚洲| 国产成人精品无人区| 脱女人内裤的视频| 国产精品久久久人人做人人爽| 波多野结衣一区麻豆| 一级毛片女人18水好多| 99在线视频只有这里精品首页| 欧美成狂野欧美在线观看| 高潮久久久久久久久久久不卡| 午夜福利18| 一本久久中文字幕| 欧美激情久久久久久爽电影 | 亚洲天堂国产精品一区在线| 国产真人三级小视频在线观看| 波多野结衣一区麻豆| 国产精品乱码一区二三区的特点 | 国产免费男女视频| av有码第一页| 丰满人妻熟妇乱又伦精品不卡| 乱人伦中国视频| 亚洲精品久久国产高清桃花| 国产伦一二天堂av在线观看| 色综合站精品国产| 欧美一级a爱片免费观看看 | 国产xxxxx性猛交| 国产熟女午夜一区二区三区| 国产亚洲av高清不卡| 日韩欧美免费精品| 色综合婷婷激情| 国产91精品成人一区二区三区| 最新在线观看一区二区三区| 色哟哟哟哟哟哟| 日韩av在线大香蕉| 一本综合久久免费| 亚洲av日韩精品久久久久久密| 午夜福利成人在线免费观看| 亚洲美女黄片视频| 国产成人av教育| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久人妻精品电影| 真人做人爱边吃奶动态| 国产精品一区二区三区四区久久 | 正在播放国产对白刺激| 午夜久久久久精精品| 亚洲一区中文字幕在线| 亚洲欧美日韩无卡精品| 91成年电影在线观看| aaaaa片日本免费| 成人特级黄色片久久久久久久| 91av网站免费观看| 亚洲av成人一区二区三| 国产精品九九99| 天堂影院成人在线观看| 久久亚洲精品不卡| 国产97色在线日韩免费| 精品国产超薄肉色丝袜足j| 黄色视频不卡| e午夜精品久久久久久久| 欧美乱码精品一区二区三区| 97碰自拍视频| 亚洲成av片中文字幕在线观看| 淫秽高清视频在线观看| 免费在线观看黄色视频的| 多毛熟女@视频| 亚洲va日本ⅴa欧美va伊人久久| 老熟妇仑乱视频hdxx| 他把我摸到了高潮在线观看| 丝袜美腿诱惑在线| 国产欧美日韩一区二区精品| 成人18禁在线播放| 在线十欧美十亚洲十日本专区| 老熟妇仑乱视频hdxx| 极品教师在线免费播放| 波多野结衣巨乳人妻| 免费在线观看完整版高清| 亚洲av电影不卡..在线观看| 搡老岳熟女国产| 三级毛片av免费| 激情视频va一区二区三区| 久久久久久久精品吃奶| 欧美国产精品va在线观看不卡| 午夜福利免费观看在线| 啦啦啦韩国在线观看视频| 精品人妻在线不人妻| 国产精品久久电影中文字幕| 美女免费视频网站| 久9热在线精品视频| 侵犯人妻中文字幕一二三四区| 搡老妇女老女人老熟妇| 咕卡用的链子| 大陆偷拍与自拍| 国产激情欧美一区二区| 国产成年人精品一区二区| √禁漫天堂资源中文www| 波多野结衣av一区二区av| 国产麻豆成人av免费视频| 久久久久久免费高清国产稀缺| 欧美在线一区亚洲| 色老头精品视频在线观看| 亚洲九九香蕉| 国产伦一二天堂av在线观看| 亚洲av日韩精品久久久久久密| 我的亚洲天堂| 99国产精品免费福利视频| 亚洲avbb在线观看| 美女免费视频网站| 久久久久国产一级毛片高清牌| 乱人伦中国视频| 在线观看日韩欧美| 国产精品综合久久久久久久免费 | 国产伦一二天堂av在线观看| 黄色视频,在线免费观看| 亚洲精品av麻豆狂野| 午夜福利,免费看| 亚洲第一电影网av| 色综合亚洲欧美另类图片| 欧美日韩亚洲国产一区二区在线观看| 精品国产一区二区三区四区第35| 成人三级做爰电影| 国产不卡一卡二| 久久久国产成人精品二区| 亚洲伊人色综图| 国产精品美女特级片免费视频播放器 | 香蕉丝袜av| 热re99久久国产66热| 欧美日本中文国产一区发布| 亚洲色图综合在线观看| 国产一区二区激情短视频| 午夜免费鲁丝| 日韩三级视频一区二区三区| 亚洲 欧美一区二区三区| 国产91精品成人一区二区三区| www.精华液| 亚洲国产看品久久| 亚洲男人的天堂狠狠| 免费少妇av软件| 欧美中文日本在线观看视频| 亚洲第一av免费看| 成人亚洲精品av一区二区| 丁香欧美五月| 欧美另类亚洲清纯唯美| 国产高清激情床上av| 欧美日韩亚洲国产一区二区在线观看| 国产不卡一卡二| 大码成人一级视频| 午夜a级毛片| 欧美另类亚洲清纯唯美| 日韩欧美一区视频在线观看| 在线播放国产精品三级| 亚洲成人国产一区在线观看| 国产成人av激情在线播放| 午夜福利,免费看| 国产日韩一区二区三区精品不卡| 天堂动漫精品| 成人欧美大片| 国产精品爽爽va在线观看网站 | 免费在线观看视频国产中文字幕亚洲| 51午夜福利影视在线观看| 老司机深夜福利视频在线观看| 久久中文字幕一级| 亚洲最大成人中文| 久久久久久久久久久久大奶| 女同久久另类99精品国产91| 18禁裸乳无遮挡免费网站照片 | 免费看a级黄色片| 国产精品自产拍在线观看55亚洲| 中文字幕高清在线视频| 叶爱在线成人免费视频播放| 欧美精品啪啪一区二区三区| 日本 欧美在线| 免费无遮挡裸体视频| 免费搜索国产男女视频| 亚洲五月色婷婷综合| 久久精品成人免费网站| 国产乱人伦免费视频| 人成视频在线观看免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久国产高清桃花| 激情视频va一区二区三区| 给我免费播放毛片高清在线观看| 欧美一区二区精品小视频在线| 中文字幕人妻熟女乱码| 99国产极品粉嫩在线观看| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品第一综合不卡| 色av中文字幕| 一区福利在线观看| 制服人妻中文乱码| 亚洲久久久国产精品| 看免费av毛片| 嫩草影院精品99| 国产精品亚洲美女久久久| 一夜夜www| 激情在线观看视频在线高清| 日韩高清综合在线| 老司机靠b影院| 国产欧美日韩综合在线一区二区| 男女午夜视频在线观看| 成人亚洲精品av一区二区| 国产高清激情床上av| 国产片内射在线| 国产精品,欧美在线| 我的亚洲天堂| 国产成人精品在线电影| 久久精品亚洲熟妇少妇任你| 中国美女看黄片| 亚洲电影在线观看av| 久9热在线精品视频| 大型黄色视频在线免费观看| 成人国语在线视频| 一a级毛片在线观看| 国产亚洲欧美精品永久| 18禁国产床啪视频网站| 国内精品久久久久精免费| 午夜a级毛片| 亚洲成a人片在线一区二区| 亚洲国产毛片av蜜桃av| 亚洲第一欧美日韩一区二区三区| 午夜免费成人在线视频| 精品一区二区三区四区五区乱码| 啦啦啦免费观看视频1| 日本三级黄在线观看| 国产精品免费一区二区三区在线| 日日摸夜夜添夜夜添小说| 亚洲狠狠婷婷综合久久图片| 岛国在线观看网站| 久久久久九九精品影院| 成人av一区二区三区在线看| 夜夜夜夜夜久久久久| 美女午夜性视频免费| 国产av一区在线观看免费| 国产成年人精品一区二区| 黄色毛片三级朝国网站| 免费女性裸体啪啪无遮挡网站| 一边摸一边抽搐一进一出视频| 亚洲精品久久国产高清桃花| 亚洲最大成人中文| 亚洲av成人一区二区三| 无遮挡黄片免费观看| 欧美成狂野欧美在线观看| 日韩国内少妇激情av| 国产精品九九99| 欧美成人免费av一区二区三区| 久久精品91蜜桃| av网站免费在线观看视频| 国产亚洲精品av在线| 欧美黄色片欧美黄色片| 精品欧美一区二区三区在线| 午夜免费成人在线视频| 欧美成狂野欧美在线观看| 亚洲成人免费电影在线观看| 免费人成视频x8x8入口观看| 欧美绝顶高潮抽搐喷水| 欧美日韩中文字幕国产精品一区二区三区 | 91精品国产国语对白视频| av在线天堂中文字幕| 成人国产综合亚洲| 成年人黄色毛片网站| 97碰自拍视频| 美女扒开内裤让男人捅视频| 亚洲久久久国产精品| 亚洲五月色婷婷综合| 成年人黄色毛片网站| 久久欧美精品欧美久久欧美| 免费久久久久久久精品成人欧美视频| 亚洲成av片中文字幕在线观看| 丝袜美足系列| 一个人观看的视频www高清免费观看 | 亚洲国产精品999在线| 成人18禁在线播放| 女生性感内裤真人,穿戴方法视频| 国产精品亚洲美女久久久| 日本 欧美在线| 亚洲一区二区三区色噜噜| 日本vs欧美在线观看视频| 男女下面进入的视频免费午夜 | 亚洲国产日韩欧美精品在线观看 | 欧美中文综合在线视频| 一二三四社区在线视频社区8| netflix在线观看网站| 亚洲黑人精品在线| 少妇粗大呻吟视频| 可以免费在线观看a视频的电影网站| 少妇被粗大的猛进出69影院| 91精品国产国语对白视频| 久久精品国产亚洲av高清一级| 亚洲三区欧美一区| 日本撒尿小便嘘嘘汇集6| 亚洲精品久久成人aⅴ小说| 国内精品久久久久精免费| 国产精品av久久久久免费| 国产精品免费视频内射| 国产麻豆成人av免费视频| 精品国产一区二区三区四区第35| 夜夜躁狠狠躁天天躁| 丝袜美腿诱惑在线| 一边摸一边做爽爽视频免费| 精品无人区乱码1区二区| 午夜a级毛片| 看免费av毛片| 两个人视频免费观看高清| 在线观看www视频免费| 99久久精品国产亚洲精品| 久久婷婷人人爽人人干人人爱 | 最近最新中文字幕大全免费视频| 国产99白浆流出| 在线av久久热| 亚洲男人天堂网一区| 在线国产一区二区在线| 一级,二级,三级黄色视频| 18禁观看日本| 国产精品久久电影中文字幕| 国产三级黄色录像| 精品福利观看| 99久久国产精品久久久| 久久婷婷人人爽人人干人人爱 | 欧美在线一区亚洲| 每晚都被弄得嗷嗷叫到高潮| 精品日产1卡2卡| 亚洲av成人av| 亚洲国产精品合色在线| 别揉我奶头~嗯~啊~动态视频| 免费少妇av软件| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久5区| 激情视频va一区二区三区| 国产97色在线日韩免费| 精品乱码久久久久久99久播| 久久精品国产综合久久久| 又大又爽又粗| 国产亚洲精品第一综合不卡| 9色porny在线观看| 少妇裸体淫交视频免费看高清 | 精品国内亚洲2022精品成人| 亚洲av熟女| 亚洲精品中文字幕在线视频| 亚洲欧美精品综合久久99| 一二三四社区在线视频社区8| 欧美性长视频在线观看| 久久香蕉激情| 一区福利在线观看| 男男h啪啪无遮挡| 最近最新免费中文字幕在线| 亚洲精品在线观看二区| 国产aⅴ精品一区二区三区波| 叶爱在线成人免费视频播放| 亚洲精品在线美女| 国产精品秋霞免费鲁丝片| 国产精品久久久人人做人人爽| 精品卡一卡二卡四卡免费| 这个男人来自地球电影免费观看| 男女午夜视频在线观看| 久久久久久久久免费视频了| 亚洲欧美一区二区三区黑人| 老熟妇仑乱视频hdxx| 欧美中文综合在线视频| 精品久久久久久久毛片微露脸| 女同久久另类99精品国产91| 国产97色在线日韩免费| 中文字幕人妻丝袜一区二区| netflix在线观看网站| 免费观看人在逋| 看黄色毛片网站| 亚洲色图av天堂| 狠狠狠狠99中文字幕| 曰老女人黄片| 老鸭窝网址在线观看| 黄色成人免费大全| 青草久久国产| 亚洲全国av大片| 成人国语在线视频| avwww免费| 久久午夜综合久久蜜桃| АⅤ资源中文在线天堂| 99久久99久久久精品蜜桃| 亚洲欧美一区二区三区黑人| 亚洲专区国产一区二区| 精品欧美国产一区二区三| 嫩草影视91久久| 狂野欧美激情性xxxx| 在线永久观看黄色视频| 免费观看精品视频网站| 亚洲性夜色夜夜综合| 亚洲一区中文字幕在线| 一本久久中文字幕| 天天躁夜夜躁狠狠躁躁| 国产单亲对白刺激| 他把我摸到了高潮在线观看| 看免费av毛片| 叶爱在线成人免费视频播放| 老司机在亚洲福利影院| 如日韩欧美国产精品一区二区三区| 久久久久亚洲av毛片大全| 777久久人妻少妇嫩草av网站| 视频在线观看一区二区三区| 国产午夜福利久久久久久| 深夜精品福利| 少妇的丰满在线观看| 欧美不卡视频在线免费观看 | 久久精品国产综合久久久| 夜夜爽天天搞| 每晚都被弄得嗷嗷叫到高潮| x7x7x7水蜜桃| 精品一品国产午夜福利视频| 午夜久久久久精精品| 在线天堂中文资源库| 欧美日韩瑟瑟在线播放| 久久精品成人免费网站| 久久人人97超碰香蕉20202| 亚洲av电影不卡..在线观看| 午夜福利,免费看| 黑人操中国人逼视频| 亚洲国产精品成人综合色| 欧美亚洲日本最大视频资源| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 男男h啪啪无遮挡|