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

    秦淮河流域五種參考作物蒸散量估算方法的比較及改進*

    2016-11-16 09:02:29秦孟晟施婷婷
    中國農業(yè)氣象 2016年4期
    關鍵詞:秦淮河適用性修正

    秦孟晟,郝 璐,施婷婷,孫 磊,孫 閣

    ?

    秦淮河流域五種參考作物蒸散量估算方法的比較及改進*

    秦孟晟,郝 璐**,施婷婷,孫 磊,孫 閣

    (南京信息工程大學生態(tài)氣象環(huán)境研究中心,南京 210044)

    基于秦淮河流域內部及周邊共7個氣象站2000-2013年的逐日氣象資料,使用FAO-56 Penman- Monteith、Irmak-Allen、Makkink、Turc、Jensen-Haise和Hargreaves共6種方法估算各站點逐日參考作物蒸散量(ET0)。以FAO-56 Penman-Monteith結果為標準,修正其余5種方法估算公式的原始經驗系數,并通過平均絕對誤差、平均相對誤差、相關系數等精度評價指標和Wilcoxon非參數檢驗法,分別從年、月尺度對比分析5種方法修正前后的估算結果,旨在獲得一種適于秦淮河流域的數據要求低,估算過程簡單,精度較高的ET0估算方法。分別以5種方法的ET0日值為自變量,P-M法ET0日值為因變量,建立逐月線性回歸方程,尋找經驗系數的修正倍數,對5種方法經驗系數進行逐月修正。結果表明,使用原始經驗系數時,年尺度上,Irmak-Allen、Makkink、Turc法存在較大誤差,Hargreaves法相關性較差,均不適于秦淮河流域;月尺度上,Irmak-Allen法在5-8月,Turc在9-11月,Hargreaves法在4月及9-11月適用性較好,其余月份誤差較大,Makkink和J-H法分別在1-12月和3-11月存在顯著差異,故5種方法均不能代替P-M法在年內12個月使用。使用修正后經驗系數,年尺度上Makkink法適用性最好,平均絕對誤差和平均相對誤差分別為14.9mm·a-1和1.4%,相關系數為0.89,無顯著差異,其次為Turc法,I-A法估算結果仍存在顯著差異,Hargreaves法相關性仍較差;月尺度上,從估算精度考慮,Turc和Makkink法搭配使用,4-10月推薦使用Turc法,其平均絕對誤差為2.1~6.1mm·mon-1,平均相對誤差為2.9%~4.3%,無顯著差異,月平均相對誤差波動較小,穩(wěn)定性好,1-3月和11-12月推薦使用Makkink法,其平均絕對誤差為1.2~2.4mm·mon-1,平均相對誤差為3.2%~5.7%,無顯著差異,月平均相對誤差波動較小,穩(wěn)定性好,從時間連續(xù)性考慮,推薦使用Hargreaves法,其平均絕對誤差為1.9~10.4mm·mon-1,平均相對誤差為3.9%~9.2%,無顯著差異,月平均相對誤差波動較小,穩(wěn)定性好。

    秦淮河流域;Penman-Monteith;輻射法;溫度法;蒸散

    蒸散是地表水循環(huán)的關鍵環(huán)節(jié)[1],也是確定作物需水量的重要基礎,一般通過參考作物蒸散量(Reference Evapotranspiration,ET0)乘以作物系數得到,因此,ET0不僅是準確估算作物需水量的關鍵因素,也是指導合理灌溉的重要參數之一[2]。準確估算參考作物蒸散量對農業(yè)水分利用以及全面分析水循環(huán)等具有重要意義[3]。

    依據估算原理及數據需要,ET0估算方法可劃分為綜合法、輻射法、溫度法和蒸發(fā)皿法[4]。一般認為綜合法精度最高,聯(lián)合國糧農組織(FAO)將FAO-56 Penman-Monteith(P-M)法作為標準方法向全球推廣[5-7]。但P-M法需要數據要求高,估算過程復雜,在資料缺失地區(qū)難以使用?;趩我粴庀笠氐妮椛浞?、溫度法需要數據少,估算過程簡單,但其原始經驗系數k均有適用區(qū)域,在適用區(qū)域以外使用需對k進行修正[8]。一些研究以P-M法估算結果為標準,針對輻射法和溫度法,結合研究區(qū)的實際情況,對不同方法的k值進行修正,使其在研究區(qū)ET0估算中獲得更好的適用效果。如張曉琳等[9]研究表明,經修正后的Hargreaves法在漢江流域適用性最好,左德鵬等[10]研究表明,經修正后的Romanenko法在石羊河流域使用性最好,趙璐等[11]研究表明,經修正后的Irmak-Alleen法在川中丘陵有良好的適用效果。

    目前國內有關ET0估算方法對比及修正的研究大多數集中在北方干旱半干旱地區(qū),其中多數研究僅限于溫度法,對輻射法和綜合法研究較少,且基本是在年尺度上進行修正,未考慮研究區(qū)不同月份的氣候差異的影響,因此,逐月修正k值對提高估算方法的適用性具有重要意義。秦淮河流域位于亞熱帶季風氣候區(qū),具有氣候濕潤,四季分明,雨熱同季,光照充足的特點。多年平均降水量為1048mm,降水充沛,但年內分配不均[12]。流域內頻發(fā)的水旱災害對生態(tài)環(huán)境與水資源的合理配置有很大影響,尋找一個數據要求低,過程簡單且精度較高的ET0估算方法對實現流域內水資源的科學管理和農業(yè)合理灌溉具有重要意義。本文以秦淮河流域歷史氣象資料為依據,估算2000-2013年秦淮河流域逐日ET0,以P-M法估算結果為標準,針對基于輻射和氣溫的5種方法ET0估算結果進行對比分析,探討各方法在秦淮河流域的適用性。結合流域不同月份的實際情況對5種方法原始經驗系數k進行修正,使其在秦淮河流域ET0估算過程中有更好的適用效果。

    1 資料與方法

    1.1 數據

    采用秦淮河流域及其周邊共7個氣象站2000-2013年逐日地面實測氣象資料,包括日平均氣溫(℃)、日最低氣溫 (℃)、日最高氣溫(℃)、10m處平均風速(m·s-1)、日照時數(h)、日平均相對濕度(%)等,數據由中國氣象科學數據共享服務網以及江蘇省氣象局提供。研究區(qū)范圍及站點分布見圖1。

    1.2 參考作物蒸散量估算方法

    選取較常用的4種綜合法(FAO-56 Penman- Monteith[13]、Turc、Jensen-Haise、Irmak-Allen),以及基于氣溫的Hargreaves法和基于輻射的Makkink法。各公式中的經驗系數k初始值選擇公式提出時的原始值。

    (1)Irmak-Allen法(I-A)是根據美國東部濕潤地區(qū)數據由FAO推導得出,所需數據種類較少,公式為

    式中,Rn為地表凈輻射(MJ·m-2·d-1);T為日平均氣溫(℃);原始經驗系數k=0.489[14]。

    (2)Makkink法(Ma)由Penman法簡化而來,公式為

    式中,Δ為飽和水汽壓-氣溫關系斜率(kPa·℃-1);γ為干濕計常數(kPa·℃-1);Rs為太陽輻射(MJ·m-2·d-1);為蒸發(fā)潛熱,2.45MJ·kg-1;原始經驗系數k=0.61[15]。

    (3)Turc法(Tu)基于歐洲西部區(qū)域提出,被推薦作為濕潤地區(qū)的估算方法,公式為

    式中,n為日照時數(h);N為可照時數(h);Ra為地外輻射(MJ·m-2·d-1);T為日平均氣溫;原始經驗系數k=0.013[16]。

    (4)Jensen-Haise法(J-H)是通過觀測取樣土壤的蒸散而得到的ET0估算方法,公式為

    式中變量與Turc法的變量一致,原始經驗系數k=0.87[17]。

    (5)Hargreaves法(Har)是在美國西部較干旱的氣候條件下建立的,在基于溫度的方法中有一定優(yōu)勢,公式為

    式中,Ra為地外輻射(MJ·m-2·d-1);Tmax和Tmin分別為日最高氣溫和日最低氣溫;原始經驗系數k=0.0023[18]。

    1.3 分析方法

    以P-M法估算結果為標準,對其余5種方法的估算結果進行比較分析。采用平均絕對誤差MAE(mm)、平均相對誤差MRE(%)和相關系數r作為精度評價的指標,各指標的計算方法參閱文獻[8]。通過非參數檢驗法Wilcoxon秩和檢驗分析不同方法估算結果與P-M法是否有顯著性差異[19]。J-H法在冬季(12、1、2月)估算結果出現負值,故年尺度上不分析J-H法的適用性,月尺度上僅分析J-H法在3-11月的適用性。MAE在年尺度上以50mm·a-1為參照標準:大于50mm·a-1為適用性差,需修正k值,小于50mm·a-1,表示適用性較好,月尺度參照標準為10mm·mon-1。MRE在年尺度上以10%為參照標準:大于10%,為適用性差,需修正k值,小于10%,表示適用性較好,月尺度參照標準亦為10%。為區(qū)分各估算方法月MRE的離散程度和波動性,引入變異系數CV(%),CV小于10%、10%~40%、40%~100%和大于100%分別表示弱變異、低等變異、中等變異和強變異。CV值越大,說明月MRE波動越大,估算方法在月尺度上穩(wěn)定性越差。r在0~0.3、0.3~0.5、0.5~0.8和0.8~1相關性分別為弱、低、中和高。Wilcoxon秩和檢驗置信區(qū)間選擇99%,即α為0.01,當Wilcoxon檢驗P值小于0.01,說明兩者存在顯著差異,需修正k值,P值大于0.01,說明兩者無顯著差異。

    考慮到秦淮河流域不同月份氣候差異的影響,以需要修正的估算方法ET0日值為自變量,P-M法ET0日值為因變量,建立逐月線性回歸方程,尋找逐月k值的修正倍數,對相應公式的k值進行逐月修正。

    2 結果與分析

    2.1 基于初始經驗系數的參考作物蒸散量估算結果分析

    2.1.1 ET0年值

    由圖2可以看出,利用各方法中原始經驗系數k估算的2000-2013年歷年ET0值與P-M法估算結果均有明顯差異,其中I-A法持續(xù)明顯偏大;Hargreaves法的估算結果差異相對較小,但具有階段性特點,2000-2011年略高,而2012-2013年略低;其余3種方法(Ture,Makkink,J-H)估算結果明顯偏小。從變化趨勢看,Makkink法估算結果的年際變化趨勢與P-M法最接近,除第一個谷值外,其余峰值和谷值均同時達到;I-A法估算結果與P-M法相比波動性較小,但也能同時達到峰值;Turc和Hargreaves法估算結果的年際變化與P-M法差異較大,峰谷值年份不同。

    從表1中的誤差指標看,Makkink法與P-M法估算結果的MAE和MRE均最大,為222.9mm·a-1和21.1%,但r為0.89。說明兩者雖誤差大,但相關性高。Hargreaves法與P-M法估算結果的MAE和MRE均最小,為39.8mm·a-1和3.7%,但r僅0.42。說明兩者雖誤差最小,但相關性低。Turc法和I-A法Wilcoxon P值小于0.01,與P-M法估算結果存在顯著差異。綜上所述,利用原始K值估算歷年ET0年值時,I-A、Makkink和Turc法誤差較大,且估算結果與P-M法存在顯著差異;Hargreaves法與P-M法r小于0.5,相關性較差,均不能代替P-M法在秦淮河流域使用。

    表1 基于原始經驗系數的4種方法與P-M法估算的ET0年值比較的統(tǒng)計特征

    注:MAE為平均絕對誤差;MRE為平均相對誤差;r為相關系數。**表示P<0.01;Wilcoxon P值<0.05,為有顯著差異,Wilcoxon P值>0.05,為無顯著差異。下同。

    Note:MAE is mean absolute error; MRE is mean relative error; r is correlation coefficient.**indicates P<0.01;Wilcoxon P<0.05,significant difference,Wilcoxon P>0.05,no significant difference. The same as below.

    2.1.2 ET0月值

    由圖3可以看出,利用各方法中原始經驗系數k估算的2000-2013年多年平均逐月ET0值與P-M法估算結果均有明顯差異。與年尺度結果一致,與P-M法估算結果相比較,I-A法估算的ET0值持續(xù)明顯偏大;Hargreaves法的估算結果差異相對較小,3-6月及10-11月略高,7-8月及12月略低,其余月份估算結果接近;Turc法在1-8月及12月偏小,9-11月接近;J-H法7-8月略小,其余月份明顯偏?。籑akkink法全年持續(xù)明顯偏小。從變化趨勢看,所有方法年內變化趨勢基本一致,均表現為1-5月逐漸增大,除J-H法外的5種方法均在6月出現小幅降低,7月所有方法均達到最大,8-12月均逐漸降低。

    從誤差指標看(表2),與P-M法相比,I-A法估算的1-4月及9-12月ET0值其MAE均大于10mm·mon-1,MRE均大于10%,且存在顯著差異,適用性較差,5-8月適用性較好;Makkink法估算的1-12月ET0值其MRE均大于10%且存在顯著差異,適用性較差;Turc法1-8月和12月存在顯著差異,適用性較差,9-11月適用性較好;J-H法在3-11月均存在顯著差異,適用性較差;Hargreaves法估算的1-12月ET0值其MRE雖均小于10%,但一些月份(1-3月、5-8月、12月)估算結果MAE較大,超過10mm·mon-1,故Hargreaves法在4月和9-11月適用性較好,其余月份適用性較差。

    秦淮河流域不同月份氣候差異較大,ET0主導因子不同,因此有必要討論各方法估算ET0月值的穩(wěn)定性。從表2中MRE指標可以看出,Makkink法和Hargreaves法年內MRE變化幅度較小,其余方法年初和年末月份MRE較大,年中月份MRE較小。對比各方法逐月MRE可見,Makkink法在4-10月MRE較大,其余方法較小,而I-A法和Turc法在1-3月和11-12月MRE遠大于Makkink法和Hargreaves法。因此,分別計算1-12月、1-3月和11-12月、4-10月3個時段CV值分析月MRE波動狀況,進而分析不同方法在不同時段估算ET0月值的穩(wěn)定性。J-H法僅分析4-10月這一時段。由圖4可以看出,Makkink法月MRE的CV值在3個時段內均最小,月MRE誤差波動較小,受不同月份氣候變化影響?。黄浯螢镠argreaves法,CV值在3個時段內均處于10%~40%,為低等變異,月MRE誤差波動較??;I-A法在1-12月、4-10月的CV值超過40%,為中等變異,在1-3月和11-12月這一時段為低等變異;Turc法在1-12月CV值超過100%,為強變異,誤差波動大,估算結果受不同月份氣候變化影響較大,4-10月時段為低等變異,但CV值接近40%,1-3月和11-12月時段為中等變異;J-H法在4-10月CV值大于40%,為中等變異。

    表2 基于原始經驗系數的5種方法與P-M法估算ET0月值比較的統(tǒng)計特征

    綜上所述,利用原始經驗系數估算ET0月值時,I-A法在5-8月適用性較好,其余月份MAE大于10mm·mon-1,MRE大于10%且估算結果與P-M法存在顯著差異,月MRE年內波動較大,穩(wěn)定性較差;Makkink法月MRE年內波動小但均大于10%,估算結果與P-M法存在顯著差異;Turc法在9-11月適用性較好,其余月份存在顯著差異,年內月MRE的CV值超過100%,穩(wěn)定性差;J-H法在3-11月估算結果與P-M法存在顯著差異;Hargreaves法4月及9-11月適用效果較好,其余月份MAE大于10mm·mon-1。月尺度上雖部分方法在一些月份適用效果較好,但5種方法均不能代替P-M法估算全年ET0月值。

    2.2 參考作物蒸散量估算公式中原始經驗系數的修正

    根據上文所述,利用原始k值估算歷年ET0年值時,I-A、Makkink和Turc法誤差較大,且估算結果與P-M法存在顯著差異;Hargreaves法相關性較差,均不能代替P-M法在秦淮河流域使用。利用原始k值估算ET0月值時,I-A法在5-8月,Turc法在9-11月,Hargreaves法在4月及9-11月適用性較好,其余月份誤差較大且存在顯著差異,Makkink和J-H法分別在1-12月和3-11月存在顯著差異,5種方法均不能代替P-M法估算全年ET0月值。

    造成上述現象的主要原因是各方法原始經驗系數k是基于各自特定的氣候條件獲取的,并不適用于秦淮河流域地區(qū)。因此,有必要對各公式中原始經驗系數k進行修正,以增加各方法的適用性。秦淮河流域不同月份氣候差異較大[12],ET0主導因子不同,以各方法逐月ET0日值為自變量,P-M法相應月份ET0日值為因變量,建立各方法逐月線性回歸方程,尋找各方法k值的逐月修正倍數,對相應公式的k值進行修正,使MAE、MRE趨向0,r趨向1,Wilcoxon P值大于0.01,月MRE的CV值趨向0。I-A、Makkink、Turc、J-H和Hargreaves法的原始經驗系數分別為0.489、0.61、0.013、0.87和0.0023。修正后各方法逐月經驗系數見表3。

    表3 修正后5種方法逐月經驗系數

    2.3 基于修正后經驗系數的參考作物蒸散量估算結果分析

    2.3.1 ET0年值

    由圖5可以看出,利用各方法修正后經驗系數k估算的2000-2013年歷年ET0值與P-M法估算結果差異變小,其中I-A法2000-2011年持續(xù)偏大,2013年略低;Makkink法2000-2009年差異較小,在2000-2001及2004年略高,其余年份略低,而2010-2013年明顯偏低;Turc法2000-2004及2006-2009年偏高,其余年份偏低,Hargreaves法2000-2011年偏高,2012-2013年偏低。從變化趨勢看,修正后各方法年際變化趨勢和峰谷值出現年份無明顯變化,主要原因是逐月修正k值對ET0年際變化趨勢無較大影響。

    由表4可見,從誤差指標看,與P-M法相比,修正后 Hargreaves法MAE、MRE變?yōu)樽畲?,但分別小于50mm·a-1和10%,r為0.45,誤差較小,相關性仍最差,估算結果與P-M法變?yōu)榇嬖陲@著差異;修正后I-A法MAE、MRE略低于Hargreaves法,r增大為0.71,誤差較小,相關性較好,估算結果與P-M法仍有顯著差異。修正后Makkink法的MAE、MRE由最大變?yōu)樽钚?,r值仍為0.89,估算結果與P-M法變?yōu)闊o顯著差異;修正后Turc法MAE、MRE明顯減小,略高于Makkink法,r增大為0.74,估算結果與P-M法變?yōu)闊o顯著差異。綜上所述,利用修正后k值估算歷年ET0年值時,Makkink法誤差最小,相關性最好,無顯著差異,適用效果最好,其次為Turc法,I-A法估算結果與P-M法有顯著差異,Hargreaves 相關性較差。

    表4 基于修正后經驗系數的4種方法與P-M法估算的ET0年值比較的統(tǒng)計特征

    2.3.2 ET0月值

    由圖6可以看出,利用各方法修正后經驗系數k估算的2000-2013年多年平均逐月ET0值與P-M法估算結果差異明顯減小,其中I-A和Hargreaves法估算的1-12月ET0值略高,Makkink和Turc法與P-M法1-12月ET0值基本相同,J-H法3-5月略低,6月略高,7-11月與P-M法基本相同。從變化趨勢看,各方法與修正前(圖3)基本一致,1-5月逐漸增大,除J-H法外的5種方法在6月均出現小幅降低,7月所有方法達到最大,8-12月逐漸降低。

    從誤差指標看,修正后5種方法誤差均明顯減小,普遍無顯著性差異(表5)。與P-M法相比,修正后I-A法估算的1-12月ET0值其MAE小于10mm·mon-1,最高為8月的8.2mm·mon-1,MRE除2月和12月外均小于10%,其中最大為1月的7.2%,1-12月估算結果與P-M法無顯著差異;修正后Makkink法估算的1-12月ET0值其MAE和MRE均小于10mm·mon-1和10%,最高為8月的6.7mm·mon-1和12月的5.7%,4月、6-7月及9-10月值與P-M法有顯著差異;修正后Turc法估算的1-12月ET0值其MAE小于10mm·mon-1,最高為8月的6.1mm·mon-1,MRE除1-2、12月外均小于10%,最高為3月的4.5%,1-12月估算結果與P-M法無顯著差異;修正后J-H法估算的3-11月ET0值其MAE小于10mm·mon-1,最高為3月的9.1mm·mon-1,MRE除3月外均低于10%,最高為11月的8.7%,1-12月估算結果與P-M法無顯著差異;修正后Hargreaves法估算的8月ET0值其MAE為10.4mm·mon-1,其余月份低于10mm·mon-1,1-12月MRE低于10%,最高為12月的9.2%;1-12月估算結果與P-M無顯著差異。

    表5 基于修正后經驗系數的5種方法與P-M法估算ET0月值比較的統(tǒng)計特征

    對比圖7和圖4可以看出,修正后Makkink法估算月ET0值時MRE的CV值在1-12月、1-3月和11-12月、4-10月3個時段均小幅上升,但仍最小,為低等變異,穩(wěn)定性較好;修正后Turc法在3個時段內的CV值均明顯下降,其中4-10月的CV值略高于Makkink法和J-H法,為20.35%,另兩個時段內CV值高于40%,為中等變異;修正后I-A法在1-12月和4-10月兩個時段CV值均明顯下降,另一個時段小幅上升,4-10月CV值為20.63%,為低等變異,其余兩個時段為中等變異;J-H法在4-10月由中等變異減弱為低等變異;修正后Hargreaves法在3個時段內CV值均無明顯變化,仍為低等變異,穩(wěn)定較好。

    綜上所述,利用修正后經驗系數估算ET0月值時,Turc法在4-10月適用效果最好,誤差最小,穩(wěn)定性僅次于Makkink法和J-H法,估算結果與P-M法無顯著差異;其次為I-A法,在4-10月僅誤差略高于Turc法,穩(wěn)定性較好,無顯著差異;Makkink法在1-3、11-12月時段適用效果最好,誤差最小,穩(wěn)定性最好,估算結果與P-M法無顯著差異;J-H法10-11月適用效果較好,3月MRE大于10%,4-9月存在顯著差異;基于溫度的Hargreave法在1-12月有較好的適用性,僅8月MAE略高,為10.4mm·mon-1,其余誤差指標均較低,估算結果與P-M法無顯著差異。月尺度上,從估算精度考慮,4-10月推薦使用Turc法,1-3月及11-12月推薦使用Makkink法;從時間連續(xù)性考慮,推薦使用Hargreaves法。

    3 結論與討論

    (1)修正前,與P-M法估算結果相比,年尺度上Irmak-Allen、Makkink、Turc法存在較大誤差和顯著性差異,Hargreaves法相關性較差,均不適用于秦淮河流域;月尺度上I-A法在5-8月,Turc在9-11月,Hargreaves法在4月及9-11月適用性較好,Makkink和J-H法分別在1-12月和3-11月存在顯著差異,5種方法均不能代替P-M法在年內12個月使用。其原因主要是各方法原始經驗系數k是在特定氣候條件下獲取的,不適于秦淮河流域。

    (2)利用回歸分析方法對各種方法k值進行逐月修正,以提高5種方法在秦淮河流域的適用性。修正后,年尺度上Makkink法適用性最好,MAE和MRE為14.9mm·a-1和1.4%,r為0.89,無顯著差異,其次為Turc法,I-A法估算結果與P-M法存在顯著差異,Hargreaves法相關性較差;月尺度上從估算精度考慮,4-10月使用Turc法,MAE為2.1~6.1mm·mon-1,MRE為2.9%~4.3%,與P-M法無顯著差異,月MRE波動較小,穩(wěn)定性好,1-3月及11-12月使用Makkink法,MAE為1.2~2.4mm·mon-1,MRE為3.2%~5.7%,與P-M法無顯著差異,月MRE波動較小,穩(wěn)定性好,從時間連續(xù)性考慮,使用Hargreaves法,僅8月MAE略高,為10.4mm·mon-1,其余月份為1.9~9.2mm·mon-1,MRE為3.9%~9.2%,與P-M法無顯著差異,月MRE波動較小,穩(wěn)定性好。

    國內外學者對不同ET0估算方法的適用性比較研究較多,但專門針對經驗系數修正的研究報道卻很少,且以往研究未考慮不同月份氣候差異帶來的影響。劉倪等[20]提出按不同月份對k進行修正,但未付諸實踐。趙捷等[8]指出,修正后Makkink法在黑河流域年尺度和月尺度的適用性最好,年尺度與本文結果一致,月尺度有所不同,本文認為Makkink法在4-10月適用效果較差,主要原因是趙捷等僅對多年平均ET0月值進行對比分析,未逐月討論修正后Makkink法適用性。張曉琳等[9]指出,修正后Hargreaves法在漢江流域適用性最好,左德鵬等[10]指出修正后Hargreaves法在西北干旱地區(qū)適用性較好,估算精度僅次于修正后Romanenko法,與本文研究結論不一致,其原因可能有以下兩點:(1)Hargreaves法是建立在美國西部較干旱氣候條件下,左德鵬等研究區(qū)域氣候條件與之相似,而秦淮河流域屬亞熱帶季風氣候,溫暖濕潤。(2)Hargreaves法在溫度法中有一定優(yōu)勢,上述兩文只限于溫度法的討論。Xystrakis等[21]指出輻射法和綜合法總體而言適用性好于溫度法,本文中Makkink為輻射法,Turc為綜合法。趙璐等[11]在川中丘陵地區(qū)對I-A法進行逐月修正,取得了良好的適用效果。但其只使用相對誤差和回歸方程決定系數兩項指標檢驗適用性,考慮不夠全面,本文綜合考慮多項指標,認為修正后I-A法在冬季適用性較差??梢姡壳癊T0估算方法在秦淮河流域的修正和使用仍有很大不確定性,需進一步完善。

    本文分別給出多種方法逐月修正后的公式,相比P-M法,數據要求低,估算過程簡單,能夠在年尺度和月尺度上較為準確地估算ET0值,為秦淮河流域內農業(yè)灌溉的合理規(guī)劃和水資源的科學管理提供科學依據。但修正方法存在局限性,文章通過建立回歸方程尋找修正倍數的方法進行經驗系數修正,未考慮深層次的物理機制,在進一步的研究中,需結合ET0對Tmin、Tmax、RH等氣候要素的敏感性分析修正經驗系數[22]。文中修正后公式還存在地域局限性,若要在其它濕潤地區(qū)使用,需綜合考慮不同地區(qū)的氣候差異,結合當地站點資料進行分析后再使用。

    參考文獻References

    [1]Trambauer P,Dutra E,Maskey S,et al.Comparison of different evaporation estimates over the African continent[J]. Hydrology and Earth System Science,2014,18(1):193-212.

    [2]陳超,龐艷梅,潘學標,等.四川地區(qū)參考作物蒸散量的變化特征及氣候影響因素分析[J].中國農業(yè)氣象,2011,32(1):35-40.

    Chen C,Pang Y M,Pan X B,et al.Analysison change of reference crop evapotranspiration and climatic influence factors in Sichuan[J].Chinese Journal of Agrometeorology, 2011,32 (1):35-40.(in Chinese)

    [3]李春強,洪克勤,李保國,等.河北省近35年(1965-2009)參考作物蒸散量的時空變化[J].中國農業(yè)氣象,2008,29(4): 414-419.

    Li C Q,Hong K Q,Li B G,et al.Analysis of temporal-spatial variations of reference evapotranspiration during 1965 and 1999 in Hebei Province[J].Chinese Journal of Agrometeorology, 2008,29 (4):414-419.(in Chinese)

    [4]曹金峰,李玉中,劉曉英,等.四種參考作物蒸散量綜合法的比較[J].中國農業(yè)氣象,2015,36(4):428-436.

    Cao J F,Li Y Z,Liu X Y,et al.Comparison of four combinationmethods for reference crop evapotranspiration[J]. Chinese Journal of Agrometeorology,2015,36(4):428-436.(in Chinese)

    [5]范伶俐,張光亞.廣東實際蒸散發(fā)和潛在蒸散發(fā)的關系研究[J].廣東海洋大學學報,2013,33 (3):71-77.

    Fan L L,Zhang G Y.Relationship between actual evapotranspiration and potential evapotranspiration in Guangdong[J].Journal of Guangdong Ocean University,2013,33(3):71-77.(in Chinese)

    [6]謝賢群,王菱.中國北方近50年潛在蒸發(fā)的變化[J].自然資源學報,2007,22(5):683-691.

    Xie X Q,Wang L.Changes of potential evaporation in Northern China over the past 50 years[J].Journal of Natural Resources,2007,22(5):683-691.(in Chinese)

    [7]Thompson J R,Green A J,Kingston D G.Potential evapotranspiration-related uncertainty in climate change impacts on river flow:an assessment for the Mekong River basin[J].Journal of Hydrology,2014,510:259-279.

    [8]趙捷,徐宗學,左德鵬,等.基于輻射的潛在蒸散發(fā)量估算方法在黑河流域的適用性分析[J].干旱區(qū)資源與環(huán)境,2013, 27(10):107 -114.

    Zhao J,Xu Z X,Zuo D P,et al.Radiation:based methods to estimate potential evapotranspiration in the Heihe River Basin[J].Journal of Arid Land Resources and Environment, 2013,27(10):107-144.(in Chinese)

    [9]張曉琳,熊立華,林琳,等.五種參考作物蒸散量發(fā)公式在漢江流域的應用[J].干旱區(qū)地理,2012,35(2):229-237.

    Zhang X L,Xiong L H,Lin L,et al.Application of five potential evapotranspiration equations in Hanjiang Basin[J]. Arid Land Geography,2012,35(2):229-237.(in Chinese)

    [10]左德鵬,徐宗學,劉兆飛.基于氣溫的潛在蒸散發(fā)量估算方法在我國西北干旱地區(qū)的應用比較[J].干旱區(qū)資源與環(huán)境,2009,23(10):123-131.

    Zuo D P,Xu Z X,Liu Z F.Comparison of temperature- based methods for estimating potential evapotranspiration in the arid region of Northwest China[J].Journal of Arid Land Resources and Environment,2009,23(10):123-131.(in Chinese)

    [11]趙璐,梁川,崔寧博,等.不同ET0計算方法在川中丘陵地區(qū)的比較及改進[J].農業(yè)工程學報,2012,28(24): 92-98.

    Zhao L,Liang C,Cui N B,et al.Comparison and improvement of different calculation methods for ET0in hilly area of central Sichuan Basin[J].Transactions of the CSAE,2012, 28(24):92-98.(in Chinese)

    [12]芮菡藝.秦淮河流域土地利用變化及其水文效應[D].南京:南京大學,2012:1-65.

    Rui H Y.Land use changes and its hydrological responses in the Qinhuai River basin[D].Nanjing:Nanjing University, 2012:1-65.(in Chinese)

    [13]Allen R G,Pereira L S,Raes D,et al.Crop evapotranspiration guidelines for computing crop water requirements[M].FAO Irrigation and Drainage Paper,1998:1-300.

    [14]Irmak S,Irmak A,Allen R G,et al.Solar and net radiation- ased equations to estimate reference evapotranspiration in humid climates[J].Journal of Irrigation and Drainage Engineering-ASCE,129(5):336 -347.

    [15]Makkink G F.Testing the Penman formula by means of lysimeters[J].Journal of the Institution of Water Engineers, 957,(11):277-288.

    [16]Turc L.Estimation of irrigation water requirements,potential evapotranspiration:a simple climatic formula evolved up to date[J].Advances in Agronomy,1961,12:13-49.

    [17]Jensen M E,Haise H R.Estimating evapotranspiration from solar radiation[J].Journal of Irrigation and Drainage Engineering-ASCE,1963,93(IR3):15-41.

    [18]Hargreaves L G,Hargreaves G H,Riley J P.Irrigation water requirements for Senegal river basin[J].Journal of Irrigation and Drainage Engineering-ASCE,1985,111 (3):265-275.

    [19]王聲鋒,段愛旺,張展羽,等.半干旱地區(qū)不同水文年Hargreaves和P-M公式的對比分析[J].農業(yè)工程學報, 2008,24(7):29-33.

    Wang S F,Duan A W,Zhang Z Y,et al.Comparison and analysis of Hargreaves equation and Penman-Monteith equation during the different hydrological years in the semi-arid region[J].Transactions of the CSAE,2008,24(7): 9-33.(in Chinese)

    [20]劉倪,夏偉,吳曉蔚,等.幾種參考作物蒸散量計算方法的比較[J].河北科技大學學報,2009,30(1):18-24.

    Liu N,Xia W,Wu X W,et al.Comparative study on calculation methods of reference crop evapotranspiration[J]. Journal of Hebei University of Science and Technology,2009,30 (1):18-24.(in Chinese)

    [21]Xystrakis F,Matzarakis A.Evaluation of 13 empirical reference potential evapotranspiration equations on the island of crete in southern Greece[J].Journal of Irrigation and Drainage Engineering-ASCE,2011,137(4): 211-222.

    [22]王曉東,馬曉群,許瑩,等.淮河流域參考作物蒸散量變化特征及主要氣象因子的貢獻分析[J].中國農業(yè)氣象,2013, 34(6):661-667.

    Wang X D,Ma X Q,Xu Y,et al.Variation of reference crop evapotranspiration and contribution of main factors in the Huaihe basin[J].Chinese Journal of Agrometeorology,2013, 34(6):661-667.(in Chinese)

    Comparison and Modification of Five Crop Reference Evapotranspiration Models for Qinhuai River Basin

    QIN Meng-sheng, HAO Lu, SHI Ting-ting, SUN Lei, SUN Ge

    (International Center For Ecology,Meteorology and Environment,Nanjing University of Information Science and Technology,Nanjing 210044,China)

    The daily reference crop evapotranspiration (ET0) was estimated using the FAO-56 Penman-Monteith and five other methods (Irmak-Allen, Makkink, Turc, Jensen-Haise and Hargreaves) and meteorological data from seven sites located in inside and surrounding areas of the Qinhuai River Basin for the period of 2000-2013.Taken the FAO-56 Penman-Montieth method as a reference, the original empirical coefficients of five other methods were calibrated. The results were analyzed with mean absolute error (MAE), mean relative error (MRE), correlation coefficient (r) and non-parametric Wilcoxon test at an annual and a monthly scale respectively. We aimed to obtain one method which requires less data with high accuracy for the Qinhuai River Basin. With daily results of five methods as independent variable and daily results of FAO-56 Penman-Montieth method as dependent variable, monthly linear regression equations were established. Monthly correction parameters could be found based on these equations. This research indicated that, Irmak-Allen, Makkink, Turc and Hargreaves methods were not applicable at an annual scale when the original empirical coefficients were used. At a monthly scale, when original empirical coefficients were used, large biases and significant differences were found for most of the methods in most months except May to August for Irmak-Allen method, September to November for Turc method, April and September to November for Hargreaves method. After model calibration, the Makkink method performed best followed by Turc method at annual scale. The MAE,MRE and r was 14.9mm·y-1,1.42% and 0.89 respectively. No significant difference existed between the results of Makkink and the FAO-56 Penman-Montieth method. In contrast, there were significant differences between the results of Irmak-Allen and FAO-56 Penman-Montieth method. The Hargreaves method was still not applicable due to a poor correlation. At a monthly scale, considering the accuracy of estimation, Makkink and Turc methods were used by collocation. Turc method was recommended to use from April to October, MAE and MRE was 2.1-6.1mm·mon-1and 2.9%-4.3%, Makkink method was recommended to use in the period from January to May and November to December, MAE and MRE was 1.2-2.4mm·mon-1and 3.2%-5.7%. No significant difference existed and coefficient of variation of monthly MRE was small in the two periods for each method. Considering continuity of time, Hargreaves method was recommended to use from January to December, MAE and MRE was 1.9-10.4mm·mon-1and 4.0%-9.2%.No significant difference existed and coefficient of variation of monthly MRE was small from January to December for Hargreaves method.

    Qinhuai River Basin; Penman-Monteith; Radiation-based method; Temperature-based method

    10.3969/j.issn.1000-6362.2016.04.002

    2015-10-29

    通訊作者。E-mail:hl_haolu@163.com

    國家自然科學基金(41571026;71373130)

    秦孟晟(1991-),博士生,主要研究方向為應用氣象。E-mial:qinmengsheng1@163.com

    猜你喜歡
    秦淮河適用性修正
    Some new thoughts of definitions of terms of sedimentary facies: Based on Miall's paper(1985)
    修正這一天
    快樂語文(2021年35期)2022-01-18 06:05:30
    強調簡潔和適用性 MICHI by Rotel X5/X3合并功放
    合同解釋、合同補充與合同修正
    法律方法(2019年4期)2019-11-16 01:07:28
    環(huán)保技術在土木工程領域的適用性探討
    小議階級分析理論的適用性
    夜航
    軟件修正
    秦淮河水冷,戰(zhàn)事幾回傷
    中國三峽(2016年11期)2017-01-15 14:00:08
    南京城與秦淮河
    中國三峽(2016年11期)2017-01-15 14:00:04
    女人高潮潮喷娇喘18禁视频| 免费在线观看黄色视频的| 午夜久久久在线观看| netflix在线观看网站| 热re99久久精品国产66热6| 丝瓜视频免费看黄片| 一进一出好大好爽视频| 最近最新中文字幕大全免费视频| 别揉我奶头~嗯~啊~动态视频| 欧美中文综合在线视频| 国产精品98久久久久久宅男小说| 日韩 欧美 亚洲 中文字幕| 捣出白浆h1v1| 一级,二级,三级黄色视频| 欧美一级毛片孕妇| 亚洲自偷自拍图片 自拍| 亚洲全国av大片| 人人妻人人爽人人添夜夜欢视频| 母亲3免费完整高清在线观看| 久久久精品国产亚洲av高清涩受| 国产高清国产精品国产三级| 9热在线视频观看99| 国产精品麻豆人妻色哟哟久久| 19禁男女啪啪无遮挡网站| a在线观看视频网站| 国产亚洲午夜精品一区二区久久| 超碰成人久久| 久久精品熟女亚洲av麻豆精品| 欧美激情极品国产一区二区三区| 香蕉国产在线看| 黑人猛操日本美女一级片| 一区二区三区国产精品乱码| 99久久人妻综合| 久久久国产精品麻豆| 后天国语完整版免费观看| av天堂久久9| 亚洲色图av天堂| 日韩视频在线欧美| 电影成人av| 国产亚洲欧美精品永久| 国产精品久久久久成人av| 12—13女人毛片做爰片一| 欧美精品一区二区免费开放| 欧美中文综合在线视频| 亚洲专区国产一区二区| 日韩欧美三级三区| 美女主播在线视频| 亚洲精品av麻豆狂野| 国产亚洲一区二区精品| 国产精品二区激情视频| 免费在线观看黄色视频的| 啦啦啦在线免费观看视频4| 热re99久久精品国产66热6| 久久精品aⅴ一区二区三区四区| 色94色欧美一区二区| 最黄视频免费看| 国产精品一区二区在线不卡| 变态另类成人亚洲欧美熟女 | av片东京热男人的天堂| 精品国产国语对白av| 丝瓜视频免费看黄片| 国产欧美日韩一区二区三| 18禁美女被吸乳视频| 国产区一区二久久| 激情在线观看视频在线高清 | 欧美日本中文国产一区发布| 黄色片一级片一级黄色片| 亚洲精品av麻豆狂野| 亚洲专区中文字幕在线| 国内毛片毛片毛片毛片毛片| 久久久国产一区二区| av电影中文网址| 一个人免费看片子| 亚洲欧美一区二区三区黑人| 久久久久视频综合| 五月天丁香电影| 国产亚洲精品第一综合不卡| cao死你这个sao货| 国产日韩一区二区三区精品不卡| 国产无遮挡羞羞视频在线观看| kizo精华| 成人黄色视频免费在线看| 久久热在线av| 纯流量卡能插随身wifi吗| 亚洲久久久国产精品| 丝袜人妻中文字幕| 亚洲 欧美一区二区三区| 纯流量卡能插随身wifi吗| 日本黄色日本黄色录像| 18禁黄网站禁片午夜丰满| 丁香六月欧美| 波多野结衣av一区二区av| 久久99一区二区三区| 免费久久久久久久精品成人欧美视频| 亚洲午夜理论影院| 免费日韩欧美在线观看| 精品久久久久久电影网| 18在线观看网站| 欧美中文综合在线视频| 不卡av一区二区三区| 日韩欧美国产一区二区入口| 亚洲成人国产一区在线观看| 精品国产乱码久久久久久男人| 亚洲一区二区三区欧美精品| 久久久国产成人免费| 黑人猛操日本美女一级片| 大片免费播放器 马上看| 国产亚洲av高清不卡| 天天影视国产精品| 国产精品.久久久| 国产成人影院久久av| 国产单亲对白刺激| 久久 成人 亚洲| 一边摸一边抽搐一进一出视频| 亚洲国产欧美日韩在线播放| 久久亚洲精品不卡| 极品教师在线免费播放| 欧美日韩视频精品一区| 国产黄频视频在线观看| 国产一区二区在线观看av| 久久久久精品国产欧美久久久| 黄色毛片三级朝国网站| 不卡av一区二区三区| 在线观看www视频免费| 男女午夜视频在线观看| videosex国产| 欧美乱妇无乱码| www.熟女人妻精品国产| 国产激情久久老熟女| 久久久国产成人免费| 日本av手机在线免费观看| 亚洲中文字幕日韩| 久久亚洲精品不卡| 亚洲av国产av综合av卡| 婷婷丁香在线五月| 中文字幕精品免费在线观看视频| 涩涩av久久男人的天堂| 久久久国产精品麻豆| 亚洲熟女精品中文字幕| 亚洲成国产人片在线观看| 菩萨蛮人人尽说江南好唐韦庄| 夜夜夜夜夜久久久久| 欧美大码av| 美国免费a级毛片| 国产一区二区三区视频了| 亚洲av成人不卡在线观看播放网| 韩国精品一区二区三区| 久热这里只有精品99| 亚洲av日韩精品久久久久久密| 视频区欧美日本亚洲| 69av精品久久久久久 | 黄色毛片三级朝国网站| 欧美成人午夜精品| 亚洲九九香蕉| 黄片播放在线免费| 波多野结衣av一区二区av| 1024视频免费在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产又爽黄色视频| av福利片在线| 国产精品98久久久久久宅男小说| 高清视频免费观看一区二区| 国产精品免费视频内射| 欧美一级毛片孕妇| 99久久99久久久精品蜜桃| 91成人精品电影| 亚洲国产精品一区二区三区在线| 久久久精品国产亚洲av高清涩受| 久久人妻av系列| 国产精品欧美亚洲77777| 精品第一国产精品| 欧美精品亚洲一区二区| 日本五十路高清| 大香蕉久久网| a级毛片在线看网站| 久久久久国产一级毛片高清牌| 9色porny在线观看| 一级黄色大片毛片| 精品一区二区三卡| 久久久水蜜桃国产精品网| 国产成人精品在线电影| 美女主播在线视频| 亚洲情色 制服丝袜| 97人妻天天添夜夜摸| 日韩成人在线观看一区二区三区| 黑丝袜美女国产一区| 亚洲一区二区三区欧美精品| 国产97色在线日韩免费| 亚洲精品久久午夜乱码| 久久 成人 亚洲| 天堂俺去俺来也www色官网| h视频一区二区三区| 99国产综合亚洲精品| 视频区图区小说| 日本av手机在线免费观看| 国产亚洲av高清不卡| tube8黄色片| 考比视频在线观看| 国产精品自产拍在线观看55亚洲 | 一个人免费在线观看的高清视频| 男女下面插进去视频免费观看| 国产精品成人在线| 女人精品久久久久毛片| 人人妻人人添人人爽欧美一区卜| 中文字幕精品免费在线观看视频| 在线亚洲精品国产二区图片欧美| 一个人免费看片子| 一本一本久久a久久精品综合妖精| 国产免费福利视频在线观看| 黄片大片在线免费观看| 久久久精品国产亚洲av高清涩受| 国产精品香港三级国产av潘金莲| a级片在线免费高清观看视频| 色在线成人网| 国产精品国产av在线观看| 国产精品九九99| 嫁个100分男人电影在线观看| 男人操女人黄网站| 正在播放国产对白刺激| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文字幕日韩| 久久毛片免费看一区二区三区| 国产成人精品久久二区二区91| 丰满迷人的少妇在线观看| 国产欧美日韩一区二区三区在线| 男女边摸边吃奶| 免费一级毛片在线播放高清视频 | 夜夜夜夜夜久久久久| 中文亚洲av片在线观看爽 | 精品熟女少妇八av免费久了| 精品少妇一区二区三区视频日本电影| 欧美日韩亚洲综合一区二区三区_| 久久久久精品国产欧美久久久| 自线自在国产av| 国产黄色免费在线视频| 中文字幕精品免费在线观看视频| 亚洲全国av大片| 亚洲九九香蕉| 99国产精品99久久久久| 精品久久久久久久毛片微露脸| 日韩大码丰满熟妇| 欧美性长视频在线观看| 国产一区二区三区综合在线观看| 欧美一级毛片孕妇| 免费看a级黄色片| 中文欧美无线码| 五月天丁香电影| 国产91精品成人一区二区三区 | 久久狼人影院| 亚洲成a人片在线一区二区| 757午夜福利合集在线观看| 国产一区二区激情短视频| svipshipincom国产片| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩综合在线一区二区| 激情在线观看视频在线高清 | 亚洲色图 男人天堂 中文字幕| 亚洲欧美一区二区三区久久| 美女高潮到喷水免费观看| 午夜成年电影在线免费观看| tube8黄色片| 这个男人来自地球电影免费观看| 亚洲国产欧美网| 黄片小视频在线播放| 日韩中文字幕视频在线看片| 久久精品国产a三级三级三级| 最近最新免费中文字幕在线| 一本大道久久a久久精品| 天堂动漫精品| 国产一区二区 视频在线| 欧美在线黄色| 中文字幕色久视频| 国产在线一区二区三区精| 人成视频在线观看免费观看| av不卡在线播放| 亚洲第一青青草原| 国产精品 欧美亚洲| av天堂在线播放| 亚洲天堂av无毛| av不卡在线播放| 亚洲精品美女久久av网站| 夫妻午夜视频| 在线观看一区二区三区激情| 99热国产这里只有精品6| 曰老女人黄片| 国产成人精品久久二区二区91| 中文字幕av电影在线播放| 亚洲精品一二三| 一本—道久久a久久精品蜜桃钙片| 中文亚洲av片在线观看爽 | 成人国产av品久久久| 国产黄色免费在线视频| 岛国在线观看网站| 亚洲精品中文字幕一二三四区 | 免费高清在线观看日韩| 久久青草综合色| 搡老岳熟女国产| 国产成人一区二区三区免费视频网站| 亚洲成人免费电影在线观看| 电影成人av| 不卡av一区二区三区| 亚洲人成77777在线视频| 免费在线观看日本一区| 国产欧美日韩精品亚洲av| 久久午夜综合久久蜜桃| 在线观看免费视频网站a站| 国产一区二区在线观看av| 亚洲精品国产区一区二| 久久午夜亚洲精品久久| 欧美日韩黄片免| 国产黄色免费在线视频| 757午夜福利合集在线观看| 成在线人永久免费视频| 9191精品国产免费久久| 精品久久久久久久毛片微露脸| 黑人欧美特级aaaaaa片| 波多野结衣av一区二区av| 在线观看免费日韩欧美大片| 一个人免费在线观看的高清视频| 人人妻人人澡人人爽人人夜夜| 日韩人妻精品一区2区三区| 国产精品 欧美亚洲| 热99国产精品久久久久久7| 亚洲美女黄片视频| 亚洲国产精品一区二区三区在线| 免费不卡黄色视频| 精品久久久精品久久久| 国产日韩欧美在线精品| 夜夜爽天天搞| 国产日韩一区二区三区精品不卡| 999久久久精品免费观看国产| 精品人妻在线不人妻| 久9热在线精品视频| 午夜日韩欧美国产| 午夜福利乱码中文字幕| 色精品久久人妻99蜜桃| 免费在线观看日本一区| 在线观看人妻少妇| 亚洲七黄色美女视频| 久久影院123| 黑人巨大精品欧美一区二区蜜桃| 日本a在线网址| 美女福利国产在线| 亚洲人成77777在线视频| 久久 成人 亚洲| 深夜精品福利| 精品少妇久久久久久888优播| 热99久久久久精品小说推荐| 91大片在线观看| 两个人免费观看高清视频| 欧美成人免费av一区二区三区 | 国产日韩欧美视频二区| 欧美亚洲 丝袜 人妻 在线| 色老头精品视频在线观看| 51午夜福利影视在线观看| 美国免费a级毛片| 日韩视频一区二区在线观看| 黑人操中国人逼视频| 性高湖久久久久久久久免费观看| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 国产精品久久电影中文字幕 | 天天躁狠狠躁夜夜躁狠狠躁| 日韩中文字幕视频在线看片| 久久狼人影院| av一本久久久久| 精品视频人人做人人爽| 老司机福利观看| 欧美精品高潮呻吟av久久| 国产精品亚洲av一区麻豆| 精品亚洲乱码少妇综合久久| 老司机午夜十八禁免费视频| videos熟女内射| 国产麻豆69| 人人妻,人人澡人人爽秒播| 女警被强在线播放| 午夜免费成人在线视频| 日韩中文字幕视频在线看片| 日韩精品免费视频一区二区三区| 久久久国产一区二区| 少妇精品久久久久久久| 首页视频小说图片口味搜索| 午夜福利欧美成人| 国产亚洲午夜精品一区二区久久| 精品人妻在线不人妻| svipshipincom国产片| 少妇被粗大的猛进出69影院| 久久人妻av系列| 男女床上黄色一级片免费看| 黄色毛片三级朝国网站| 天天添夜夜摸| 热99国产精品久久久久久7| 精品国产一区二区三区久久久樱花| 热re99久久国产66热| 午夜福利免费观看在线| 一本久久精品| 国产成人精品无人区| 亚洲va日本ⅴa欧美va伊人久久| 免费观看人在逋| 热99国产精品久久久久久7| 免费观看人在逋| 欧美国产精品va在线观看不卡| 亚洲九九香蕉| 中文字幕av电影在线播放| 国产精品二区激情视频| 美女福利国产在线| 狠狠婷婷综合久久久久久88av| 91老司机精品| 少妇被粗大的猛进出69影院| 亚洲精品乱久久久久久| av不卡在线播放| 法律面前人人平等表现在哪些方面| 欧美性长视频在线观看| 亚洲,欧美精品.| 国产单亲对白刺激| 国产欧美日韩一区二区三区在线| 久久久久久久精品吃奶| 一区二区三区国产精品乱码| 久久天堂一区二区三区四区| 少妇的丰满在线观看| 好男人电影高清在线观看| 久久亚洲真实| 中文字幕精品免费在线观看视频| 国产精品免费一区二区三区在线 | 丁香欧美五月| 大香蕉久久成人网| 亚洲自偷自拍图片 自拍| 亚洲久久久国产精品| 夜夜夜夜夜久久久久| 变态另类成人亚洲欧美熟女 | 国产成人精品久久二区二区91| 久久青草综合色| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 久热爱精品视频在线9| 久久人妻熟女aⅴ| 国产精品电影一区二区三区 | 国产日韩欧美在线精品| 精品一区二区三区四区五区乱码| 在线观看免费视频日本深夜| 大型av网站在线播放| 国产成人精品久久二区二区免费| 少妇被粗大的猛进出69影院| 老司机午夜福利在线观看视频 | 大陆偷拍与自拍| 欧美乱码精品一区二区三区| 日本五十路高清| 俄罗斯特黄特色一大片| 亚洲 欧美一区二区三区| 久久精品国产亚洲av香蕉五月 | 国产区一区二久久| 亚洲专区中文字幕在线| 欧美日韩成人在线一区二区| 国产精品久久久久成人av| 高清视频免费观看一区二区| 两人在一起打扑克的视频| 美女视频免费永久观看网站| 少妇裸体淫交视频免费看高清 | 人人妻人人添人人爽欧美一区卜| 99re在线观看精品视频| 国产成人免费观看mmmm| 999久久久精品免费观看国产| 大型av网站在线播放| 咕卡用的链子| 美女高潮到喷水免费观看| 色视频在线一区二区三区| 天天影视国产精品| 国产97色在线日韩免费| 99久久人妻综合| 18禁黄网站禁片午夜丰满| 美女视频免费永久观看网站| 亚洲欧美一区二区三区黑人| 大片免费播放器 马上看| 日韩欧美一区二区三区在线观看 | 黑丝袜美女国产一区| 成人18禁高潮啪啪吃奶动态图| 侵犯人妻中文字幕一二三四区| av有码第一页| 国产在线精品亚洲第一网站| 亚洲av成人不卡在线观看播放网| 天天操日日干夜夜撸| 岛国在线观看网站| 国产福利在线免费观看视频| 超碰成人久久| 亚洲精品中文字幕一二三四区 | 久久人妻福利社区极品人妻图片| 久久影院123| 欧美在线黄色| 自拍欧美九色日韩亚洲蝌蚪91| 日韩成人在线观看一区二区三区| 国产精品成人在线| 亚洲自偷自拍图片 自拍| 国产精品九九99| 国产亚洲午夜精品一区二区久久| 超碰97精品在线观看| 色94色欧美一区二区| 黄色丝袜av网址大全| av天堂在线播放| 日本一区二区免费在线视频| 热re99久久精品国产66热6| 精品国产乱码久久久久久男人| 久久久久久久久免费视频了| 欧美午夜高清在线| 欧美另类亚洲清纯唯美| 亚洲国产毛片av蜜桃av| 日本wwww免费看| 亚洲精品中文字幕在线视频| 久久性视频一级片| 超碰97精品在线观看| 亚洲av国产av综合av卡| 久久久久网色| 在线看a的网站| 搡老熟女国产l中国老女人| 国产男女超爽视频在线观看| 人人妻人人添人人爽欧美一区卜| 在线永久观看黄色视频| 国产精品 欧美亚洲| 日韩欧美一区二区三区在线观看 | 久久性视频一级片| 国产无遮挡羞羞视频在线观看| 亚洲人成电影观看| 一区二区三区激情视频| 成年人黄色毛片网站| 不卡一级毛片| 69精品国产乱码久久久| 欧美精品高潮呻吟av久久| 国产在线精品亚洲第一网站| 欧美激情高清一区二区三区| 女同久久另类99精品国产91| 99热网站在线观看| 1024香蕉在线观看| 搡老熟女国产l中国老女人| a在线观看视频网站| 十八禁高潮呻吟视频| 亚洲精品国产色婷婷电影| 日韩三级视频一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 又大又爽又粗| 精品人妻1区二区| 久久ye,这里只有精品| 成年人免费黄色播放视频| 国产精品熟女久久久久浪| 免费人妻精品一区二区三区视频| 国产深夜福利视频在线观看| 亚洲综合色网址| 捣出白浆h1v1| 五月天丁香电影| a在线观看视频网站| videos熟女内射| 亚洲精品国产色婷婷电影| 性高湖久久久久久久久免费观看| 国产成人一区二区三区免费视频网站| 日韩视频在线欧美| 国精品久久久久久国模美| 老司机影院毛片| 一区二区三区乱码不卡18| 亚洲九九香蕉| 精品国产乱码久久久久久小说| 精品少妇久久久久久888优播| 99精品欧美一区二区三区四区| 高清在线国产一区| kizo精华| 精品国产乱子伦一区二区三区| 99re6热这里在线精品视频| 最近最新中文字幕大全电影3 | 久久久久久人人人人人| 精品国产乱码久久久久久小说| 久久久久久人人人人人| 亚洲人成电影免费在线| 午夜免费成人在线视频| 中文字幕色久视频| 人妻 亚洲 视频| 成人影院久久| 午夜福利,免费看| 男男h啪啪无遮挡| 亚洲性夜色夜夜综合| 老司机亚洲免费影院| 男女边摸边吃奶| 亚洲熟女毛片儿| 亚洲人成电影免费在线| 国产男靠女视频免费网站| 免费久久久久久久精品成人欧美视频| 日韩精品免费视频一区二区三区| 多毛熟女@视频| 久久中文看片网| 久久 成人 亚洲| 亚洲九九香蕉| 黄色 视频免费看| 欧美另类亚洲清纯唯美| 亚洲人成77777在线视频| 动漫黄色视频在线观看| 亚洲国产欧美在线一区| 十八禁网站网址无遮挡| 国产成人精品久久二区二区免费| 国产精品.久久久| 高清欧美精品videossex| 在线看a的网站| 手机成人av网站| 久久九九热精品免费| 老熟妇仑乱视频hdxx| 亚洲人成伊人成综合网2020| 国产成人系列免费观看| 欧美国产精品一级二级三级| 高清视频免费观看一区二区| 青青草视频在线视频观看| 日本黄色日本黄色录像| 久久 成人 亚洲| 大陆偷拍与自拍| 又黄又粗又硬又大视频| av超薄肉色丝袜交足视频| 国产日韩欧美视频二区| 国产精品美女特级片免费视频播放器 | 日韩欧美国产一区二区入口| 国产亚洲欧美精品永久| av有码第一页|