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

    沂河流域1961-2010年降雨侵蝕力時(shí)空分布特征

    2017-08-01 10:00:17李子君林錦闊姜愛(ài)霞
    水土保持通報(bào) 2017年3期
    關(guān)鍵詞:沂河侵蝕性氣象站

    李子君, 林錦闊, 姜愛(ài)霞, 張 晨

    (山東師范大學(xué) 地理與環(huán)境學(xué)院, 山東 濟(jì)南 250014)

    沂河流域1961-2010年降雨侵蝕力時(shí)空分布特征

    李子君, 林錦闊, 姜愛(ài)霞, 張 晨

    (山東師范大學(xué) 地理與環(huán)境學(xué)院, 山東 濟(jì)南 250014)

    [目的] 分析沂河流域近50 a的降雨量和降雨侵蝕力的時(shí)空變化特征,為流域水土流失防治及土地利用合理規(guī)劃等工作提供參考。[方法] 利用沂河流域及周邊12個(gè)氣象站1961—2010年的日降雨數(shù)據(jù),基于日降雨信息的月降雨侵蝕力模型計(jì)算流域多年平均降雨侵蝕力,采用Mann-Kendall非參數(shù)檢驗(yàn)法及析取Kriging內(nèi)插法分析流域降雨量和降雨侵蝕力的時(shí)空變化特征。[結(jié)果] 沂河流域降雨量和降雨侵蝕力空間分布上呈現(xiàn)出由西南向北逐級(jí)遞減的變化趨勢(shì)。多年平均降雨量為789.41 mm,多年平均降雨侵蝕力為2 626.09(MJ·mm)/(hm2·h·a),兩者都在1965年產(chǎn)生突變;降雨量和降雨侵蝕力年內(nèi)分布主要集中在夏季(6—8月),分別占全年比例的63.02%和71.22%,二者最大值都出現(xiàn)在7月,且秋季對(duì)流域多年降雨量的減少趨勢(shì)貢獻(xiàn)最多,夏季的降雨侵蝕力上升幅度最大。[結(jié)論] 沂河流域的降雨量和降雨侵蝕力空間分布趨勢(shì)相似,不同月份的降雨量與降雨侵蝕力差異不同。

    沂河流域; 降雨侵蝕力; 日降雨量; 時(shí)空分布

    文獻(xiàn)參數(shù): 李子君, 林錦闊, 姜愛(ài)霞, 等.沂河流域1961—2010年降雨侵蝕力時(shí)空分布特征[J].水土保持通報(bào),2017,37(3):86-92.DOI:10.13961/j.cnki.stbctb.2017.03.015; Li Zijun, Lin Jinkuo, Jiang Aixia, et al. Spatio-temporal variations characteristics of rainfall erosivity during 1961—2010 in Yihe Basin[J]. Bulletin of Soil and Water Conservation, 2017,37(3):86-92.DOI:10.13961/j.cnki.stbctb.2017.03.015

    降水是引起土壤侵蝕的主要?jiǎng)恿σ蛩豙1]。降雨侵蝕力(R)作為一項(xiàng)重要指標(biāo),不僅被用來(lái)表征降雨對(duì)土壤侵蝕潛在能力的大小,而且是反映水土流失對(duì)全球氣候變化響應(yīng)的最佳參數(shù),成為許多土壤侵蝕模型(如USLE,RUSLE,AGNPS,CREAMS,EPIC等)的主要參數(shù)之一[2]。在相同地表?xiàng)l件下的水蝕區(qū),同一地區(qū)不同年份,或同一年份不同地區(qū)的侵蝕呈現(xiàn)差異,皆源于降雨侵蝕力的不同[3]。因此,研究區(qū)域降雨侵蝕力差異,特別是以流域?yàn)閱卧慕涤昵治g力(R)值的準(zhǔn)確計(jì)算以及時(shí)空分布特征,可以更好的揭示氣候變化與土壤侵蝕的響應(yīng)關(guān)系,為流域水土保持工作提供參考。1958年Wischmeier等[4]首次提出了運(yùn)用次降雨總動(dòng)能E與30 min最大雨強(qiáng)I30的乘積EI30,計(jì)算降雨侵蝕力的經(jīng)典算法。但由于長(zhǎng)序列的降雨資料不易獲取,導(dǎo)致經(jīng)典算法的使用限制較多,因此國(guó)內(nèi)外許多學(xué)者開(kāi)始對(duì)降雨侵蝕力的新算法及基于年、月和日尺度的簡(jiǎn)易方程進(jìn)行了深入分析[1,5-6]。由于日降雨量的資料較易獲取且相比年、月數(shù)據(jù)的降雨侵蝕力估算結(jié)果更為精確,因此多用于分析降雨侵蝕力的空間分布和時(shí)間變化特征[7]。已有諸多研究表明,降雨侵蝕力年際波動(dòng)顯著,年內(nèi)分配不均且主要集中于汛期[8],具有明顯的空間異質(zhì)性[9]且與降雨量空間分布趨勢(shì)基本一致[10]。

    沂河流域?qū)儆诒狈酵潦絽^(qū),是魯中南低山丘陵區(qū)的一部分,是淮河流域的暴雨中心,也是山東省乃至淮河流域水土流失最為嚴(yán)重的區(qū)域[11]。流域土壤表層疏松,土壤涵水能力低,嚴(yán)重的土壤侵蝕狀況影響了區(qū)域的生態(tài)安全和社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展。以往的研究中有學(xué)者指出,魯中南地區(qū)年降雨侵蝕力具有上升趨勢(shì),此區(qū)域不僅包含省級(jí)水土流失重點(diǎn)治理區(qū)為魯中南低山丘陵區(qū),還有國(guó)家級(jí)水土流失重點(diǎn)治理區(qū)為沂蒙山區(qū),因此,加強(qiáng)魯中南地區(qū)的降雨侵蝕力的變化預(yù)測(cè)和趨勢(shì)研究,有助于減少山東省水土流失防治的不確定性和成本[12]。本文擬基于沂河流域內(nèi)部及周邊區(qū)域12個(gè)氣象站點(diǎn)1961—2010年共50 a的日降雨量資料,分析流域降雨量和降雨侵蝕力的空間分布及年際、年內(nèi)變化特征,以期更好地了解流域降雨的潛在侵蝕差異,為流域水土流失防治及土地利用合理規(guī)劃等工作提供參考,同時(shí)為沂蒙山區(qū)生態(tài)安全屏障建設(shè)提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    沂河發(fā)源于山東省沂源縣魯山南麓,流經(jīng)魯中南和蘇北的13個(gè)縣(市),在江蘇新沂駐入駱馬湖,全長(zhǎng)超過(guò)500 km,流域面積17 325 km2,屬淮河流域。本研究中的沂河流域指沂河流域臨沂站以上的部分,流域面積10 034 km2,地理位置位于117°25′—119°49′E,33°30′—36°20′N(xiāo)。沂河流域在地貌上屬構(gòu)造剝蝕堆積平原區(qū)地貌,由于地殼運(yùn)動(dòng)和長(zhǎng)期風(fēng)雨的侵蝕等外力作用,造成流域地表斷裂縱橫、山巒起伏、嶺谷相間的地貌景觀。流域氣候類(lèi)型屬暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候,受大陸氣團(tuán)和海洋性氣團(tuán)交替影響,具有光照充足、四季分明、夏季炎熱、冬季寒冷的特點(diǎn)。多年平均氣溫為13.05~14.3 ℃,南北相差較大。流域降水在年內(nèi)分布極不均勻,多集中于夏秋季節(jié)。整個(gè)流域水系呈極不對(duì)稱(chēng)的扁長(zhǎng)半樹(shù)葉形,較大的一級(jí)支流多從右岸匯入,主要有東汶河、蒙河、祊河。流域土壤類(lèi)型主要有粗骨土、褐土、棕壤、紅黏土、潮土等。植被類(lèi)型主要有兩種:樹(shù)林植被和草甸植被。樹(shù)林植被以油松為代表,一般分布在海拔800 m以上的山腰和其他土壤瘠薄的山脊山坡上,此外還有加拿大楊、二楊、多種歐美楊、旱柳等樹(shù)種。灌草叢以黃荊、酸棗、胡枝子、黃草、白草、結(jié)縷草為主。由于長(zhǎng)期的人為活動(dòng),平地和低矮的丘陵大都開(kāi)為農(nóng)田,農(nóng)作物、多種果樹(shù)等栽培植物。流域以水力侵蝕為主,2010年侵蝕面積(輕度及以上土壤侵蝕面積)為4 252.97 km2,占流域總面積的42.42%。侵蝕嚴(yán)重的區(qū)域大都位于河流中上游地區(qū),主要發(fā)生在海拔600 m以下區(qū)域,其中高程位于200—400 m的區(qū)域侵蝕面積最大。

    2 數(shù)據(jù)來(lái)源與研究方法

    2.1 數(shù)據(jù)來(lái)源及處理

    本文采用1961-2010年50 a間沂河流域內(nèi)及周邊均勻分布的12個(gè)站點(diǎn)的連續(xù)24 h降水量,站點(diǎn)分布均勻,符合區(qū)域降水空間分布特征的分析要求,數(shù)據(jù)來(lái)自于中國(guó)氣象數(shù)據(jù)共享服務(wù)網(wǎng)(http:∥data.cma.cn)和山東省氣象局。臨沂站、莒縣站和沂源站為國(guó)家級(jí)臺(tái)站,其他氣象站為省級(jí)臺(tái)站。但是有些站實(shí)測(cè)降水資料系列相對(duì)較短或缺少記錄(表1)。為了使流域內(nèi)的降水量在不同時(shí)段中具有同一比較基礎(chǔ),以減少或消除因雨量站代表性欠佳而帶來(lái)的偏差,必須利用相關(guān)性較好的雨量站的長(zhǎng)資料系列值進(jìn)行插補(bǔ)、展延。

    利用SPSS軟件對(duì)沂河流域12個(gè)氣象站降水量數(shù)據(jù)進(jìn)行回歸分析后,選取相關(guān)系數(shù)高且具有相對(duì)完整降水資料的站,與需要插補(bǔ)降水資料的站建立線性回歸方程(表2)。經(jīng)檢驗(yàn),所選插補(bǔ)站降水量與需插補(bǔ)站降水量之間的顯著性概率為0.000,小于0.001,說(shuō)明所建線性回歸模型回歸效果極為顯著。

    表1 沂河流域氣象站已有數(shù)據(jù)

    表2 沂河流域氣象站插補(bǔ)、延展情況

    2.2 研究方法

    2.2.1 降雨侵蝕力的計(jì)算 降雨侵蝕力的計(jì)算方法可以分為標(biāo)準(zhǔn)法和模型法。標(biāo)準(zhǔn)法采用次降雨中的最大30 min 降雨強(qiáng)度I30和降雨動(dòng)能E的乘積EI30來(lái)計(jì)算,但需要詳細(xì)的降雨過(guò)程資料,難以實(shí)現(xiàn)。而縱觀各種模型法,雖計(jì)算方法簡(jiǎn)單,但大都未明顯包含雨強(qiáng)的概念,這與降雨侵蝕力提出的初衷相違背[13],影響了計(jì)算的精度。

    本文采用楊軒等人基于日降雨量計(jì)算月降雨侵蝕力的復(fù)合因子模型[14],此模型通過(guò)侵蝕性雨量(Pmer)與侵蝕性降雨日數(shù)(Dmer)之比,很好地將降雨強(qiáng)度反映在模型當(dāng)中,提高了降雨侵蝕力模型的估算精度,其表達(dá)式為:

    (1)

    式中:Rm——月降雨侵蝕力〔(J·mm)/(m2·h)〕;Pm——月降雨量(mm);Pmer——月侵蝕性雨量(mm),即某月內(nèi)日降雨量≥12 mm[2]的降雨量總和;Dmer——某月日降雨量≥12 mm的降雨天數(shù);Dm——某月內(nèi)日降雨量≥0.6 mm的月降雨天數(shù);Zm——某月內(nèi)降雨量最大的1日日降雨量(mm);Z3 m——某月日降雨量≥12 mm時(shí),該月內(nèi)降雨量最大的3日日降雨量之和(mm);Pdmax——某月內(nèi)最大日降雨量(mm)。若要計(jì)算年降雨侵蝕力,則分別計(jì)算一年內(nèi)12個(gè)月份的降雨侵蝕力后累加即可,即:

    (2)

    式中:i——一年內(nèi)的第i月,i=1,2,…,12;Ra——年降雨侵蝕力〔(MJ·mm)/(hm2·h)〕。

    2.2.2 Mann-Kendall秩次相關(guān)檢驗(yàn)法和突變檢驗(yàn)法 用于檢驗(yàn)時(shí)間序列變化的統(tǒng)計(jì)方法很多,其中非參數(shù)檢驗(yàn)法優(yōu)點(diǎn)是不需要樣本遵從一定的分布,也不受少數(shù)異常值的干擾[15],被世界氣象組織推薦并廣泛應(yīng)用于分析降水、徑流、氣溫、水質(zhì)等方面,故文中通過(guò)計(jì)算運(yùn)用Mann-Kendall檢驗(yàn)法對(duì)沂河流域降雨量和降雨侵蝕力進(jìn)行趨勢(shì)檢驗(yàn)和突變分析。

    假定X1,X2,…Xn為時(shí)間序列變量,n為時(shí)間序列的長(zhǎng)度,M-K法定義了統(tǒng)計(jì)量S:

    (3)

    式中:xj,xk——第j年和第k年的相應(yīng)測(cè)量值,且k>j。

    (4)

    式中:Z——正態(tài)分布的統(tǒng)計(jì)量,方差Var(s)=n(n-1)(2n+5)/18。在給定的α置信水平上,如果│Z│≥Z1-a/2,則拒絕原假設(shè),即在α置信水平上,時(shí)間序列數(shù)據(jù)存在明顯的上升或下降趨勢(shì),Z為正值表示增加趨勢(shì),為負(fù)值表示減少趨勢(shì)。取顯著性水平α=0.05水平下,對(duì)應(yīng)Z的絕對(duì)值為1.96。其變化趨勢(shì)的大小用β表示:

    (5)

    若β>0,表示呈上升趨勢(shì);若β<0,表示呈下降趨勢(shì)[16]。

    突變分析主要是檢驗(yàn)時(shí)間序列的不連續(xù)性。設(shè)降雨侵蝕力序列為x1,x2,…,xn,Sk表示第i個(gè)樣本的xi>xj(1≤j≤i)的累計(jì)數(shù),定義統(tǒng)計(jì)量:

    (j=1,2,…,i;k=1,2,…,n)

    (6)

    在時(shí)間序列隨機(jī)獨(dú)立的假定情況下,Sk的均值和方差分別為:

    E[Sk]=k(k-1)/4,

    var[Sk]=k(k-1)(2k+5)/72 (1≤k≤n)

    (7)

    將Sk標(biāo)準(zhǔn)化:

    (8)

    其中UF1=0。給定顯著性水平α,若│UFk│>Ua則表明序列存在明顯的趨勢(shì)變化。所有UFk可組成一條曲線。按時(shí)間序列逆序Xn,Xn-1,…,X1重復(fù)上述公式(7)—(9)過(guò)程,同時(shí)使UBk=-UFk(k=n,n-1,…,1),其中UB1=0。若UBk和UFk的值大于0,則表明序列呈上升趨勢(shì),小于0表明序列呈下降趨勢(shì),當(dāng)他們超過(guò)臨界線時(shí),表明上升或下降趨勢(shì)顯著。2條曲線出現(xiàn)交點(diǎn),且交點(diǎn)在臨界線之間,那么交點(diǎn)對(duì)應(yīng)的時(shí)刻是降雨侵蝕力突變開(kāi)始的年份。根據(jù)上述方法,確定降雨量和降雨侵蝕力突變開(kāi)始的年份[17]。

    2.2.3 滑動(dòng)T檢驗(yàn) 滑動(dòng)T檢驗(yàn)是通過(guò)考察2組樣本平均值的差異是否顯著來(lái)檢驗(yàn)突變。其基本思想是把一氣候序列中的2段子序列均值有無(wú)顯著性差異看做來(lái)自2個(gè)總體均值有無(wú)顯著性差異的問(wèn)題來(lái)檢驗(yàn)。如果兩段子序列的均值差異超過(guò)了一定顯著性水平,可以認(rèn)為均值發(fā)生了突變。但子序列時(shí)段選擇帶有人為性,為避免任意選擇子序列長(zhǎng)度造成突變點(diǎn)的漂移,具體使用時(shí),可反復(fù)變動(dòng)子序列長(zhǎng)度進(jìn)行比較,以提高計(jì)算結(jié)果的可靠性。

    (9)

    方程(9)遵從自由度v=n1+n2-2的t分布[17]。

    3 結(jié)果與分析

    3.1 降雨量和降雨侵蝕力空間分布

    由公式(1)和(2) 計(jì)算出流域各月降雨侵蝕力,經(jīng)累加得出年降雨侵蝕力后,計(jì)算出各氣象站多年平均降雨侵蝕力。由于不同插值方法反映整體趨勢(shì)與局部分布格局有所差異,而基于地統(tǒng)計(jì)學(xué)的析取Kriging內(nèi)插法對(duì)區(qū)域降雨侵蝕力空間插值,可以突出地反映降雨侵蝕力的局部信息,在細(xì)節(jié)上的表現(xiàn)更為靈敏[18]。因此本文利用ArcGIS10.2地統(tǒng)計(jì)學(xué)拓展模塊將各離散測(cè)站的多年平均降雨量和多年平均降雨侵蝕力進(jìn)行析取Kriging插值分析,進(jìn)而得到流域多年平均降雨量等值線和多年平均降雨侵蝕力的空間分布(圖1)。

    圖1 沂河流域降雨量與降雨侵蝕力空間分布

    由圖1可知,沂河流域多年平均降雨量范圍為710.11~986.38mm,平均值為789.41mm,降雨侵蝕力范圍為2 093.70~3 220.27(MJ·mm)/(hm2·h·a),平均值為2 626.09(MJ·mm)/(hm2·h·a),其中費(fèi)縣和平邑縣南部的區(qū)域?yàn)榻涤昵治g力的高值區(qū),達(dá)到2 800(MJ·mm)/(hm2·h·a)以上,而降雨量的高值區(qū)位于平邑縣;流域降雨量和降雨侵蝕力的低值區(qū)位于沂源縣北部區(qū)域,降雨侵蝕力小于2 200(MJ·mm)/(hm2·h·a)??傮w來(lái)看,沂河流域的降雨量和降雨侵蝕力空間分布趨勢(shì)相似,都呈現(xiàn)出由西南向北逐級(jí)遞減的變化規(guī)律。以往并未有學(xué)者對(duì)該流域的降雨侵蝕力進(jìn)行分析,但有學(xué)者對(duì)整個(gè)沂蒙山區(qū)的降雨侵蝕力研究顯示[10],沂蒙山區(qū)的降雨侵蝕力從空間分布上呈現(xiàn)從南向西北、東北2個(gè)方向遞減的趨勢(shì),這與本研究結(jié)論基本一致(圖1)。

    思雨在想怎樣把這兩天發(fā)生的長(zhǎng)發(fā)絲風(fēng)波說(shuō)明白,能讓小姨子理解自己,搬動(dòng)小姨子這個(gè)救兵,幫自己解決這個(gè)難題。

    將12個(gè)氣象站點(diǎn)多年平均降雨量、多年平均侵蝕性降雨量和多年平均降雨侵蝕力繪制成圖2。

    圖2 沂河流域各氣象站點(diǎn)的年均降雨量、 年均侵蝕性降雨量和年均降雨侵蝕力

    由圖2可知,沂河流域12個(gè)點(diǎn)站中年均降雨侵蝕力與年均降雨量、年均侵蝕性降雨量具有一定的協(xié)同性。流域最高降雨侵蝕力出現(xiàn)在棗莊站,最低降雨侵蝕力出現(xiàn)在沂源站。其中平邑站雖然具有最高的年均降雨量,但是由于其年均侵蝕性降雨量較低,因此降雨侵蝕力反而不高。通過(guò)線性回歸法構(gòu)建流域各氣象站多年平均降雨侵蝕力與降雨量的相關(guān)關(guān)系,可以得出年均降雨侵蝕力與年均侵蝕性降雨量的相關(guān)系數(shù)R=0.899(p<0.01)明顯高于年均降雨量的相關(guān)系數(shù)R=0.874(p<0.01)。由此看來(lái),沂河流域的降雨侵蝕力受侵蝕性降雨量的影響更大。

    3.2 降雨量和降雨侵蝕力的年際變化

    從沂河流域降雨量和降雨侵蝕力年際變化趨勢(shì)來(lái)看(圖3),降雨量的變化趨勢(shì)略有下降,而降雨侵蝕力呈上升趨勢(shì)。其中降雨量的最小值出現(xiàn)在2002年,只有471.16mm,最大值出現(xiàn)在1964年,達(dá)到1 183.58mm。降雨侵蝕力的最小值出現(xiàn)在1976年,僅為1 510.85(MJ·mm)/(hm2·h·a),而2000年以后,流域降雨侵蝕力的波動(dòng)的程度較為劇烈,2003年達(dá)到了研究時(shí)段內(nèi)的最大值,且較常年偏高,為4 674.09(MJ·mm)/(hm2·h·a),可以看出,逐年降雨侵蝕力的年際差相差較大,最高值為最低值的3.09倍。流域50a來(lái),降雨侵蝕力與降雨量的變化趨勢(shì)基本一致,但也有部分年份兩者關(guān)系不一致:如2003年雖然沒(méi)有最高的降雨量,但是由于在整個(gè)研究時(shí)段內(nèi)侵蝕性降雨量較高,因此其年降雨侵蝕力的值也是最大,而1964和1974年雖然降雨量和侵蝕性降雨量都很高但是其降雨侵蝕力相對(duì)較低,主要原因是計(jì)算降雨侵蝕力時(shí),采用了極端降雨量(Zm,Z3m)等新指標(biāo),導(dǎo)致在降雨量和侵蝕性降雨量相差不大的情況下,Zm,Z3m較大的月份降雨侵蝕力相對(duì)較小。

    通過(guò)計(jì)算沂河流域及各氣象站降雨量、降雨侵蝕力的傾斜度β值并進(jìn)行Mann-Kendall波動(dòng)趨勢(shì)檢驗(yàn)(表3),可以得出與圖3同樣的規(guī)律。自1961—2010年,流域降雨量略有下降,下降幅度僅為-0.91。各氣象站點(diǎn)的降雨量除新泰站、沂源站、鄒城站、泗水站呈上升趨勢(shì)外,其余各站均呈下降趨勢(shì),其中臨沂站降雨量的下降幅度β值以及Z值最大,但所有氣象站的降雨量變化趨勢(shì)均未超過(guò)0.05的顯著性水平。流域降雨侵蝕力呈現(xiàn)上升趨勢(shì),上升幅度為6.45。除臨沂站、莒縣站和蒼山站的降雨侵蝕力呈下降趨勢(shì)外,其余各站均呈上升趨勢(shì),且沂南站侵蝕力上升幅度最大達(dá)到16.39,但各站的變化趨勢(shì)在0.05顯著性水平上均不顯著。由此說(shuō)明,研究時(shí)段內(nèi),流域降雨侵蝕力存在一定隨機(jī)波動(dòng),單個(gè)氣象站降雨量和降雨侵蝕力的變化對(duì)流域整體的氣候變化趨勢(shì)影響并不大。

    圖3 1961-2010年沂河流域年降雨量、年降雨侵蝕力變化表3 沂河流域及各氣象站降雨量及降雨侵蝕力趨勢(shì)檢驗(yàn)

    氣象站降雨量?jī)A斜度βZ降雨侵蝕力傾斜度βZ 沂河流域-0.91-0.506.450.80新泰站0.610.386.700.70蒙陰站-0.99-0.493.450.43平邑站-1.73-0.847.200.75沂水站-1.66-0.904.390.38沂南站-0.16-0.1316.391.36鄒城站0.030.0312.570.82棗莊站-0.48-0.456.790.65臨沂站-3.18-1.85-8.37-0.90沂源站0.470.227.501.00莒縣站-2.43-1.28-8.45-0.80蒼山站-1.02-0.43-1.15-0.18泗水站0.540.2211.091.19

    注:Z為正態(tài)分布統(tǒng)計(jì)量。下同。

    為分析沂河流域降雨量、降雨侵蝕力的突變狀況,對(duì)流域降雨量、降雨侵蝕力時(shí)間序列進(jìn)行Mann-Kendall突變分析。

    圖4 1961-2010年沂河流域年降雨量突變分析

    圖5 1961-2010年沂河流域年降雨侵蝕力突變分析

    3.3 降雨量和降雨侵蝕力的年內(nèi)變化

    沂河流域季節(jié)降雨量以及降雨侵蝕力時(shí)間變化趨勢(shì)結(jié)果顯示(表4),春、冬季節(jié)降雨量呈現(xiàn)上升趨勢(shì),夏、秋季節(jié)降雨量呈現(xiàn)下降趨勢(shì)。4個(gè)季節(jié)中秋季降雨量的下降趨勢(shì)最大,說(shuō)明秋季對(duì)流域多年降雨量的減少趨勢(shì)貢獻(xiàn)最多。降雨侵蝕力除秋季呈下降趨勢(shì)外,其他季節(jié)均為上升趨勢(shì),且夏季降雨侵蝕力上升幅度最大,Kendall傾斜度為7.21,但在0.05顯著性水平上并不顯著。冬季降雨侵蝕力上升幅度雖然最少僅為0.15,但其Z值卻最大,達(dá)到1.52。有研究認(rèn)為,在全球氣候變暖的大環(huán)境下,部分冬季降雪會(huì)轉(zhuǎn)變?yōu)榻涤?,從而?dǎo)致侵蝕力的升高[19],這與本研究的結(jié)果相吻合。

    表4 沂河流域季節(jié)降雨量、降雨侵蝕 力Kendall傾斜度β及趨勢(shì)檢驗(yàn)

    分別計(jì)算流域50 a來(lái)各月平均降雨量、各月平均降雨侵蝕力占年均降雨量和年均降雨侵蝕力的百分比,并以此繪制降雨量和降雨侵蝕力的年內(nèi)分配圖(圖6)。從圖中可以看出,沂河流域降雨量和降雨侵蝕力年內(nèi)分布規(guī)律基本一致均呈單峰型。降雨量和降雨侵蝕力年內(nèi)分布主要集中在夏季(6—8月),分別占全年比例的63.02%和71.22%,二者最大值都出現(xiàn)在7月,分別為28.92%和31%。流域6—9月的降雨侵蝕力明顯高于降雨量,而其余月份降雨侵蝕力均低于降雨量。雨滴動(dòng)能的大小直接影響擊濺侵蝕的強(qiáng)度,沂河流域夏季暴雨集中,次降雨的雨強(qiáng)大,雨滴動(dòng)能大,而其他月份次降雨的雨強(qiáng)弱,雨滴動(dòng)能小,因此導(dǎo)致不同月份降雨量與降雨侵蝕力的差異不同。

    圖6 1961-2010年沂河流域降雨量 與降雨侵蝕力年內(nèi)分配特征

    4 討論與結(jié)論

    (1) 沂河流域降雨量和降雨侵蝕力總體呈現(xiàn)出由西南向北逐級(jí)遞減的變化趨勢(shì),年均降雨量變化范圍為710.11~986.38 mm,平均值為789.41 mm,降雨侵蝕力范圍為2 093.70~3 220.27 (MJ·mm)/(hm2·h·a),平均值為2 626.09(MJ·mm)/(hm2·h·a),流域年均降雨侵蝕力受侵蝕性降雨的影響更大。其中降雨侵蝕力的高值區(qū)位于費(fèi)縣和平邑縣南部的區(qū)域,達(dá)到2 800(MJ·mm)/(hm2·h·a)以上,降雨量的高值區(qū)位于平邑縣;流域降雨量和降雨侵蝕力的低值區(qū)位于沂源縣北部區(qū)域,降雨侵蝕力小于2 200(MJ·mm)/(hm2·h·a)。

    (2) 50 a來(lái),降雨量變化趨勢(shì)略有下降,下降幅度為-0.91,最小值出現(xiàn)在2002年,只有471.16 mm,最大值出現(xiàn)在1964年,達(dá)到1 183.58 mm。降雨侵蝕力呈上升趨勢(shì),上升幅度為6.45,降雨侵蝕力的最小值出現(xiàn)在1976年,僅為1 510.85(MJ·mm)/(hm2·h·a),2003年達(dá)到了最大值,且較常年偏高,為4 674.09(MJ·mm)/(hm2·h·a),流域降雨量和降雨侵蝕力都在1965年產(chǎn)生突變。研究時(shí)段內(nèi),降雨侵蝕力存在一定隨機(jī)波動(dòng),單個(gè)氣象站降雨量和降雨侵蝕力的顯著變化對(duì)流域整體的氣候變化趨勢(shì)影響并不大。

    (3) 流域降雨量和降雨侵蝕力年內(nèi)分布主要集中在夏季(6—8月),分別占全年比例的63.02%和71.22%,二者最大值都出現(xiàn)在7月,分別為28.92%和31%,且不同月份的降雨量與降雨侵蝕力差異不同。秋季對(duì)流域多年降雨量的減少趨勢(shì)貢獻(xiàn)最多,夏季降雨侵蝕力上升幅度最大,但各季節(jié)的變化趨勢(shì)均不顯著。

    (4) 此外劉正佳等[20]同樣采用基于日降雨量計(jì)算月降雨侵蝕力復(fù)合因子模型,對(duì)沂蒙山區(qū)降雨侵蝕力的時(shí)空分布分析得出,流域6—9月的降雨侵蝕力明顯高于降雨量,而其余月份降雨侵蝕力均低于降雨量。本研究與之相比,得出相同規(guī)律,但是本研究計(jì)算沂河流域降雨侵蝕力多年平均值為2 626.09(MJ·mm)/(hm2·h·a),劉正佳等采用相同的計(jì)算方法得到沂蒙山區(qū)多年平均值為1 777.5(MJ·mm)/(hm2·h·a),而楊紹洋等[10]采用日雨量修正模型得到沂蒙山區(qū)多年平均值為4 133.92(MJ·mm)/(hm2·h·a)。由此可見(jiàn),計(jì)算方法相同,選取不同的區(qū)域尺度和時(shí)間尺度,或選取的研究區(qū)域相同,但采用不同的計(jì)算公式得到的降雨侵蝕力R值均存在較大差異,因此如何采取統(tǒng)一的標(biāo)準(zhǔn)來(lái)更為精確的計(jì)算降雨侵蝕力的大小,在今后的研究中可進(jìn)一步探討。此外,本研究侵蝕性降雨采用的是黃土高原坡面侵蝕的侵蝕性降雨雨量標(biāo)準(zhǔn)(P≥12 mm)[2],此標(biāo)準(zhǔn)在沂河流域的適用性如何,還需進(jìn)一步考證。

    [1] 章文波,付金生.不同類(lèi)型雨量資料估算降雨侵蝕力[J].資源科學(xué),2003,25(1):35-41.

    [2] 謝云,劉寶元,章文波.侵蝕性降雨標(biāo)準(zhǔn)研究[J].水土保持學(xué)報(bào),2000,14(4):6-11.

    [3] 卜兆宏,董勤瑞,周伏建,等.降雨侵蝕力因子新算法的初步研究[J].土壤學(xué)報(bào),1992,29(4):408-418.

    [4] Wischmeier W H, Smith D D. Rainfall energy and its relation-ship to soil loss[J]. Transactions of American Geophysical Union, 1958,39(3):285-291.

    [5] Yu B, Hashim G M, Eusof Z. Estimating the R-factor using limited rainfall data: A case study from Peninsular Malaysia[J].Journal of Soil and Water Conservation, 2001,56(2):101-105.

    [6] 章文波,謝云,劉寶元.利用日雨量計(jì)算降雨侵蝕力的方法研究[J].地理科學(xué),2002,22(6):705-711.

    [7] 賴(lài)成光,陳曉宏,王兆禮,等.珠江流域1960—2012年降雨侵蝕力時(shí)空變化特征[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(8):159-167.

    [8] 劉惠英,任洪玉,張長(zhǎng)偉,等.三峽庫(kù)區(qū)香溪河流域降雨侵蝕力的時(shí)空分布特征[J].中國(guó)水土保持科學(xué),2015,13(3):1-7.

    [9] 范建容,嚴(yán)冬,郭祥.GIS支持下的長(zhǎng)江上游降雨侵蝕力時(shí)空分布特征分析[J].水土保持研究,2010,17(1):92-96.

    [10] 楊韶洋,劉霞,姚孝友,等.沂蒙山區(qū)降雨侵蝕力空間分布推算方法[J].中國(guó)水土保持科學(xué),2015,13(2):1-10.

    [11] 趙善倫,尹民,孫希華.山東省水土流失經(jīng)濟(jì)損失與生態(tài)價(jià)值損失評(píng)估[J].經(jīng)濟(jì)地理,2002,22(5):616-619.

    [12] 馬良,左長(zhǎng)清,尹忠東,等.山東省降雨侵蝕力多年變化特征分析[J].中國(guó)水土保持科學(xué),2010,8(4):79-85.

    [13] 王萬(wàn)忠,焦菊英,郝小品,等.中國(guó)降雨侵蝕力R值的計(jì)算與分布(I)[J].水土保持學(xué)報(bào),1995,9(4):5-18.

    [14] 楊軒,梁音,方繼青,等.基于日降雨信息的月降雨侵蝕力模型[J].土壤學(xué)報(bào),2010,47(2):216-222.

    [15] 任洪玉,劉惠英.三峽庫(kù)區(qū)大寧河流域降雨侵蝕力的時(shí)空變化[J].水土保持通報(bào),2016,36(3):1-7.

    [16] 康淑媛,張勃,柳景峰,等.基于Mann-Kendall法的張掖市降水量時(shí)空分布規(guī)律分析[J].資源科學(xué),2009,31(3):501-508.

    [17] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京:氣象出版社,2007.

    [18] 門(mén)明新,宇振榮,許皞.基于地統(tǒng)計(jì)學(xué)的河北省降雨侵蝕力空間格局研究[J].中國(guó)農(nóng)業(yè)科學(xué),2006,39(11):2270-2277.

    [19] 劉國(guó)彬,梁一民.黃土高原草地植被恢復(fù)與土壤抗沖性形成過(guò)程(I):草地植被恢復(fù)生物量特征[J].水土保持研究,1997,4(S1):102-110,121.

    [20] 劉正佳,劉占仁,王絲絲,等.基于日降雨的沂蒙山區(qū)降雨侵蝕力時(shí)空變化研究[J].水土保持研究,2012,19(2):34-37.

    Spatio-temporal Variations Characteristics of Rainfall Erosivity During 1961-2010 in Yihe Basin

    LI Zijun, LIN Jinkuo, JIANG Aixia, ZHANG Chen

    (SchoolofGeographyandEnvironment,ShandongNormalUniversity,Jinan,Shandong250014,China)

    [Objective] Analyzing the spatio-temporal variations of rainfall and rainfall erosivity of the Yihe Basin, to provide some references for soil and water loss control and land use planning. [Methods] A monthly rainfall erosivity model was adopted to calculate the rainfall erosivity by using daily precipitation data that were collected from 12 weather stations of Yihe Basin and its surrounding regions during 1961—2010. The spatio-temporal variations of rainfall and rainfall erosivity were analyzed by Mann-Kendall nonparametric tests and disjunctive Kriging space interpolation method. [Results] The spatial distributions of rainfall and rainfall erosivity in Yihe Basin showed a decreasing trends from the southwest part to the north part. The yearly averaged rainfall and rainfall erosivity were 789.41 mm and 2 626.09(MJ·mm)/(hm2·h·a). Both of them were observed changed abruptly in 1965. The monthly rainfall and rainfall erosivity in Yihe Basin were mainly distributed in summer(June to August), accounting for 63.02% and 71.22% of the corresponding annual values respectively, and their maximum value were all found in July. The decrease proportion of rainfall in Autumn contributed mostly to the yearly decreasing trend of rainfall. The rainfall erosivity in summer had the greatest contribution. [Conclusion] The spatial distribution of rainfall and rainfall erosivity was similar in Yihe Basin. The monthly rainfall was different from the monthly rainfall erosivity.

    Yihe Basin; rainfall erosivity; daily precipitation; spatio-temporal variation

    2016-09-20

    2016-10-08

    國(guó)家自然科學(xué)基金項(xiàng)目“沂河流域土地覆被變化和蓄水工程水文效應(yīng)的分析與模擬”(41101079); 山東省自然基金(ZR2015DM012); 山東師范大學(xué)大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃國(guó)家級(jí)項(xiàng)目(201610445017)

    李子君(1972—),女(漢族),山東省莒縣人,博士,教授,碩士生導(dǎo)師,主要從事土地利用變化方面的研究。E-mail:lizj@sdjzu.edu.cn。

    姜愛(ài)霞(1961—),女(漢族),山東省威海市人,本科,高級(jí)實(shí)驗(yàn)師,主要從事土壤侵蝕試驗(yàn)方面的研究。E-mail:jaixia79@sina.com。

    A

    1000-288X(2017)03-0086-07

    S157.1, TP79

    猜你喜歡
    沂河侵蝕性氣象站
    青海省1961—2018 年侵蝕性降水時(shí)空變化特征
    珠峰上架起世界最高氣象站
    基于能量的坡面侵蝕性徑流及其水沙傳遞關(guān)系
    沂河濕地生態(tài)系統(tǒng)現(xiàn)狀及生態(tài)恢復(fù)對(duì)策
    園林科技(2021年1期)2022-01-19 03:14:02
    2011—2017年蘇北沿海侵蝕性降雨特征研究
    江蘇水利(2020年5期)2020-06-14 09:08:44
    心靈氣象站
    風(fēng)吹過(guò)沂河淌(組詩(shī))
    米非司酮配伍化療治療侵蝕性葡萄胎的療效觀察
    沂河臨沂站洪水預(yù)報(bào)影響因素分析
    自動(dòng)氣象站應(yīng)該注意的一些防雷問(wèn)題
    青青草视频在线视频观看| 国产欧美日韩一区二区三区在线 | 人妻夜夜爽99麻豆av| 3wmmmm亚洲av在线观看| 视频中文字幕在线观看| a级毛色黄片| 成人一区二区视频在线观看| 欧美日韩综合久久久久久| 中文欧美无线码| 国产男女内射视频| 久久久久精品性色| 国产免费又黄又爽又色| 身体一侧抽搐| 大码成人一级视频| 免费观看的影片在线观看| 亚洲精品456在线播放app| 免费黄频网站在线观看国产| av黄色大香蕉| 啦啦啦中文免费视频观看日本| 精品少妇久久久久久888优播| 一二三四中文在线观看免费高清| 欧美日韩在线观看h| 国产精品av视频在线免费观看| 亚洲av免费在线观看| 少妇的逼好多水| 夫妻午夜视频| 亚洲精品久久午夜乱码| 天天躁日日操中文字幕| 一二三四中文在线观看免费高清| 欧美xxxx黑人xx丫x性爽| 五月天丁香电影| 青春草国产在线视频| 伊人久久国产一区二区| 欧美潮喷喷水| 久久鲁丝午夜福利片| 国产精品一区二区性色av| 神马国产精品三级电影在线观看| 精华霜和精华液先用哪个| 国产午夜精品久久久久久一区二区三区| 婷婷色麻豆天堂久久| 26uuu在线亚洲综合色| 久久99热这里只有精品18| 精品久久久久久久久av| 日韩强制内射视频| 色播亚洲综合网| 日韩亚洲欧美综合| 狂野欧美激情性bbbbbb| 嘟嘟电影网在线观看| 丝袜喷水一区| kizo精华| 久久久久久久久久久免费av| 26uuu在线亚洲综合色| 女人十人毛片免费观看3o分钟| 干丝袜人妻中文字幕| 色5月婷婷丁香| 欧美亚洲 丝袜 人妻 在线| 97在线视频观看| 丝袜脚勾引网站| 国产探花在线观看一区二区| 亚洲精品亚洲一区二区| 久久久久久久精品精品| 男插女下体视频免费在线播放| 午夜激情福利司机影院| 欧美极品一区二区三区四区| 亚洲av在线观看美女高潮| 久久精品熟女亚洲av麻豆精品| 国产又色又爽无遮挡免| 18+在线观看网站| 欧美激情国产日韩精品一区| 黄色日韩在线| 国产高潮美女av| 日韩中字成人| 一级片'在线观看视频| 久久综合国产亚洲精品| 国语对白做爰xxxⅹ性视频网站| 最近2019中文字幕mv第一页| 国产爽快片一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲在久久综合| 肉色欧美久久久久久久蜜桃 | 免费观看无遮挡的男女| 男人爽女人下面视频在线观看| 亚洲不卡免费看| 中国美白少妇内射xxxbb| 成年女人看的毛片在线观看| 午夜福利在线在线| 日韩欧美精品v在线| 日日摸夜夜添夜夜添av毛片| 国内精品美女久久久久久| 最近最新中文字幕免费大全7| 欧美zozozo另类| 99热全是精品| 欧美潮喷喷水| 午夜福利网站1000一区二区三区| 五月玫瑰六月丁香| 男插女下体视频免费在线播放| 99热这里只有精品一区| 国产一区二区三区综合在线观看 | 午夜视频国产福利| 大又大粗又爽又黄少妇毛片口| 伦理电影大哥的女人| 2022亚洲国产成人精品| 亚洲精品乱久久久久久| 狂野欧美激情性xxxx在线观看| 女人十人毛片免费观看3o分钟| 国产人妻一区二区三区在| 别揉我奶头 嗯啊视频| 91精品一卡2卡3卡4卡| 国产男女超爽视频在线观看| 神马国产精品三级电影在线观看| 视频区图区小说| 国产精品爽爽va在线观看网站| 亚洲高清免费不卡视频| 成年免费大片在线观看| 日本一本二区三区精品| 精品酒店卫生间| 国产精品偷伦视频观看了| 欧美性感艳星| 精品人妻一区二区三区麻豆| av网站免费在线观看视频| 91久久精品国产一区二区三区| 建设人人有责人人尽责人人享有的 | 一本久久精品| 亚洲精品一二三| 亚洲成人精品中文字幕电影| 久久6这里有精品| 亚洲图色成人| 大话2 男鬼变身卡| 国产黄片视频在线免费观看| 老师上课跳d突然被开到最大视频| 欧美日韩一区二区视频在线观看视频在线 | 嫩草影院入口| 国产亚洲av片在线观看秒播厂| 亚洲欧美精品自产自拍| 色视频在线一区二区三区| 中文字幕亚洲精品专区| 大话2 男鬼变身卡| 亚洲色图综合在线观看| 精品亚洲乱码少妇综合久久| 亚洲精品乱码久久久久久按摩| 欧美日韩视频高清一区二区三区二| 国产欧美日韩一区二区三区在线 | 少妇的逼水好多| 日日啪夜夜撸| 久久97久久精品| 尤物成人国产欧美一区二区三区| 91午夜精品亚洲一区二区三区| 18禁在线播放成人免费| av黄色大香蕉| 免费高清在线观看视频在线观看| 国模一区二区三区四区视频| 99热这里只有是精品50| 久久99精品国语久久久| 欧美成人一区二区免费高清观看| 国产精品国产av在线观看| 国产 一区精品| 天天躁日日操中文字幕| 欧美高清成人免费视频www| 精华霜和精华液先用哪个| 日韩欧美精品免费久久| 国产精品不卡视频一区二区| 亚洲综合色惰| 精品一区二区三区视频在线| 欧美+日韩+精品| 日韩不卡一区二区三区视频在线| 亚洲最大成人手机在线| 赤兔流量卡办理| 免费观看a级毛片全部| 丰满少妇做爰视频| 午夜免费男女啪啪视频观看| 国产午夜精品久久久久久一区二区三区| av一本久久久久| 青春草国产在线视频| 日日摸夜夜添夜夜添av毛片| 成人毛片60女人毛片免费| 99热国产这里只有精品6| 欧美性猛交╳xxx乱大交人| 麻豆精品久久久久久蜜桃| 国产精品国产av在线观看| 午夜福利在线观看免费完整高清在| 精品一区二区免费观看| 好男人视频免费观看在线| 亚洲欧美中文字幕日韩二区| 狠狠精品人妻久久久久久综合| 国产精品国产av在线观看| 91久久精品电影网| 丰满人妻一区二区三区视频av| 久久97久久精品| 热re99久久精品国产66热6| 日韩一区二区视频免费看| 亚洲国产高清在线一区二区三| 国产伦精品一区二区三区视频9| 在线观看美女被高潮喷水网站| 99热全是精品| 欧美变态另类bdsm刘玥| 国产乱人视频| 男女边吃奶边做爰视频| 亚洲成人久久爱视频| 日韩电影二区| 久久国产乱子免费精品| 国产淫片久久久久久久久| 欧美+日韩+精品| 亚洲精品乱久久久久久| 成人二区视频| 国内少妇人妻偷人精品xxx网站| 在线观看美女被高潮喷水网站| 新久久久久国产一级毛片| 中文字幕av成人在线电影| 少妇猛男粗大的猛烈进出视频 | 亚洲av电影在线观看一区二区三区 | 神马国产精品三级电影在线观看| av卡一久久| 可以在线观看毛片的网站| 亚洲熟女精品中文字幕| 下体分泌物呈黄色| 三级经典国产精品| 少妇裸体淫交视频免费看高清| av在线老鸭窝| 精华霜和精华液先用哪个| 久久久色成人| 99精国产麻豆久久婷婷| 我要看日韩黄色一级片| 天堂网av新在线| 97人妻精品一区二区三区麻豆| 亚洲四区av| 特级一级黄色大片| 亚洲欧洲日产国产| 国产大屁股一区二区在线视频| 亚洲真实伦在线观看| 插阴视频在线观看视频| 亚洲国产日韩一区二区| 男女边吃奶边做爰视频| 日本一本二区三区精品| 搡女人真爽免费视频火全软件| 乱系列少妇在线播放| 亚洲色图av天堂| 丰满乱子伦码专区| 新久久久久国产一级毛片| 2021天堂中文幕一二区在线观| 能在线免费看毛片的网站| 中文资源天堂在线| 男女边吃奶边做爰视频| 伊人久久国产一区二区| 欧美日韩在线观看h| 久久国内精品自在自线图片| 少妇高潮的动态图| 欧美性猛交╳xxx乱大交人| 日本一本二区三区精品| 日韩制服骚丝袜av| 香蕉精品网在线| 午夜福利高清视频| 春色校园在线视频观看| 亚洲国产日韩一区二区| av在线天堂中文字幕| 又黄又爽又刺激的免费视频.| av免费观看日本| 乱系列少妇在线播放| 三级男女做爰猛烈吃奶摸视频| 久久ye,这里只有精品| av在线老鸭窝| 日韩精品有码人妻一区| 亚洲成人一二三区av| 在线看a的网站| 色播亚洲综合网| 亚洲欧美成人精品一区二区| 国产av不卡久久| 亚洲图色成人| 国产一区二区三区av在线| 国产人妻一区二区三区在| 夫妻性生交免费视频一级片| 啦啦啦中文免费视频观看日本| a级一级毛片免费在线观看| 久久99热这里只有精品18| 国产精品一二三区在线看| 午夜福利在线观看免费完整高清在| 免费观看无遮挡的男女| 免费看光身美女| 精品午夜福利在线看| 亚洲久久久久久中文字幕| av在线观看视频网站免费| 日本欧美国产在线视频| 精品国产露脸久久av麻豆| 日韩成人av中文字幕在线观看| 成人欧美大片| 大片电影免费在线观看免费| 又黄又爽又刺激的免费视频.| 日韩伦理黄色片| 国产乱来视频区| 国产免费一级a男人的天堂| 丝袜美腿在线中文| 搞女人的毛片| av在线天堂中文字幕| 午夜免费男女啪啪视频观看| 男女边吃奶边做爰视频| 成人一区二区视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美日韩一区二区视频在线观看视频在线 | 美女高潮的动态| 亚洲av欧美aⅴ国产| 日韩欧美精品免费久久| 欧美最新免费一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 中国美白少妇内射xxxbb| 免费不卡的大黄色大毛片视频在线观看| 在线亚洲精品国产二区图片欧美 | 日本三级黄在线观看| 免费看a级黄色片| 国产精品一区www在线观看| 九九在线视频观看精品| 亚洲婷婷狠狠爱综合网| 伦理电影大哥的女人| 在线观看美女被高潮喷水网站| 日韩欧美一区视频在线观看 | 精品久久久久久久久亚洲| 久久人人爽人人片av| 欧美精品一区二区大全| 日韩亚洲欧美综合| 身体一侧抽搐| 久久久久网色| 国产免费一级a男人的天堂| 丰满少妇做爰视频| 大陆偷拍与自拍| 亚洲,一卡二卡三卡| 伦精品一区二区三区| 国产一区二区亚洲精品在线观看| 亚洲三级黄色毛片| 日韩三级伦理在线观看| 又黄又爽又刺激的免费视频.| 亚洲怡红院男人天堂| 国产欧美日韩一区二区三区在线 | 亚洲欧美清纯卡通| 伊人久久国产一区二区| av在线蜜桃| 少妇裸体淫交视频免费看高清| 国产精品麻豆人妻色哟哟久久| 精品久久久久久久久亚洲| 色视频www国产| 中文字幕av成人在线电影| 在线 av 中文字幕| av在线蜜桃| 99久国产av精品国产电影| 色网站视频免费| 日日摸夜夜添夜夜爱| 高清视频免费观看一区二区| 精品少妇久久久久久888优播| 色5月婷婷丁香| 久久久精品欧美日韩精品| 美女被艹到高潮喷水动态| 国产成人精品婷婷| 中文字幕久久专区| 日本猛色少妇xxxxx猛交久久| 国产色婷婷99| 五月玫瑰六月丁香| 中文字幕久久专区| 亚洲欧洲日产国产| 国产成人午夜福利电影在线观看| 激情五月婷婷亚洲| 亚洲精华国产精华液的使用体验| 舔av片在线| 国产精品人妻久久久久久| 国产一区二区亚洲精品在线观看| 午夜精品一区二区三区免费看| 亚洲av免费在线观看| 最近最新中文字幕大全电影3| 亚洲av男天堂| tube8黄色片| 日韩制服骚丝袜av| 国产 一区 欧美 日韩| 99久久精品热视频| videos熟女内射| 亚洲国产av新网站| 国产大屁股一区二区在线视频| 美女国产视频在线观看| 伊人久久国产一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 深夜a级毛片| 久久久久国产精品人妻一区二区| 婷婷色麻豆天堂久久| 秋霞在线观看毛片| 精品人妻视频免费看| 亚洲精品中文字幕在线视频 | 国产精品99久久久久久久久| 校园人妻丝袜中文字幕| 白带黄色成豆腐渣| 一区二区三区乱码不卡18| 熟女电影av网| 久久久久网色| 亚洲精品日韩在线中文字幕| 干丝袜人妻中文字幕| 精品人妻视频免费看| 嘟嘟电影网在线观看| 欧美日韩精品成人综合77777| 久热久热在线精品观看| 在线观看一区二区三区| 日日啪夜夜爽| 亚洲精品色激情综合| 欧美一级a爱片免费观看看| 国产国拍精品亚洲av在线观看| 成人综合一区亚洲| 成人特级av手机在线观看| 亚洲精华国产精华液的使用体验| 亚洲精品久久午夜乱码| 久久久久久久久久久免费av| 亚洲va在线va天堂va国产| 日本免费在线观看一区| 午夜爱爱视频在线播放| av在线观看视频网站免费| 超碰97精品在线观看| 亚洲成人中文字幕在线播放| 别揉我奶头 嗯啊视频| 亚洲av.av天堂| 国产一区有黄有色的免费视频| 777米奇影视久久| 成人欧美大片| 深爱激情五月婷婷| 亚洲色图av天堂| 精品久久久精品久久久| 两个人的视频大全免费| 久久久久性生活片| 3wmmmm亚洲av在线观看| 中文在线观看免费www的网站| 国产一区亚洲一区在线观看| 精品久久久久久久人妻蜜臀av| 国产视频首页在线观看| 欧美激情国产日韩精品一区| 91精品伊人久久大香线蕉| 欧美日韩综合久久久久久| 99久国产av精品国产电影| 欧美精品一区二区大全| 久久99热6这里只有精品| 亚洲天堂av无毛| 久久久久久久国产电影| 久久久久久久久大av| tube8黄色片| 国产黄a三级三级三级人| 18+在线观看网站| 观看美女的网站| 国产精品.久久久| 18禁裸乳无遮挡免费网站照片| 亚洲欧美精品自产自拍| 久久人人爽人人片av| 少妇人妻一区二区三区视频| 亚洲美女搞黄在线观看| 日韩大片免费观看网站| 国产午夜精品一二区理论片| 国产高清有码在线观看视频| 人体艺术视频欧美日本| 欧美日韩精品成人综合77777| 久久久久久久大尺度免费视频| 久久久成人免费电影| 国产高清三级在线| 久久久久久久午夜电影| 成人高潮视频无遮挡免费网站| 久久午夜福利片| av在线蜜桃| 日韩一区二区三区影片| 3wmmmm亚洲av在线观看| 真实男女啪啪啪动态图| 成人特级av手机在线观看| 99re6热这里在线精品视频| 天天躁日日操中文字幕| 永久免费av网站大全| 麻豆精品久久久久久蜜桃| 日韩伦理黄色片| 在线观看人妻少妇| 国国产精品蜜臀av免费| 看免费成人av毛片| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美一区二区三区黑人 | 两个人的视频大全免费| 欧美丝袜亚洲另类| av在线观看视频网站免费| 一级黄片播放器| 亚洲精品乱码久久久久久按摩| 男人爽女人下面视频在线观看| 欧美xxⅹ黑人| 最近手机中文字幕大全| 最近中文字幕2019免费版| 色播亚洲综合网| 一个人看视频在线观看www免费| 我要看日韩黄色一级片| 美女被艹到高潮喷水动态| 自拍欧美九色日韩亚洲蝌蚪91 | 在线观看三级黄色| 久久精品国产亚洲av天美| 日韩欧美精品v在线| 青春草国产在线视频| 日日啪夜夜撸| 少妇人妻精品综合一区二区| 亚洲av欧美aⅴ国产| 亚洲国产精品成人久久小说| 国产免费福利视频在线观看| 18禁裸乳无遮挡免费网站照片| 99精国产麻豆久久婷婷| av卡一久久| 久久久久久久久久久免费av| 69av精品久久久久久| 日韩国内少妇激情av| 午夜亚洲福利在线播放| 老汉色av国产亚洲站长工具| 青青草视频在线视频观看| 亚洲精品国产区一区二| 日本91视频免费播放| 国产麻豆69| 中文字幕制服av| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品国产三级国产专区5o| 天天躁夜夜躁狠狠久久av| 丝袜脚勾引网站| 精品福利永久在线观看| 97人妻天天添夜夜摸| 一级片免费观看大全| 欧美精品亚洲一区二区| 天天躁夜夜躁狠狠久久av| 国产精品一区二区精品视频观看| 天天操日日干夜夜撸| 曰老女人黄片| av网站免费在线观看视频| 水蜜桃什么品种好| 色精品久久人妻99蜜桃| 少妇人妻精品综合一区二区| 1024视频免费在线观看| 一边摸一边做爽爽视频免费| 十八禁人妻一区二区| 国产亚洲一区二区精品| 制服诱惑二区| 欧美亚洲日本最大视频资源| 男人操女人黄网站| 欧美 日韩 精品 国产| 大码成人一级视频| 捣出白浆h1v1| a级毛片黄视频| 在线观看免费午夜福利视频| h视频一区二区三区| 啦啦啦视频在线资源免费观看| av.在线天堂| 亚洲av综合色区一区| av福利片在线| 高清在线视频一区二区三区| 日韩大片免费观看网站| 日韩精品免费视频一区二区三区| 高清不卡的av网站| 国产女主播在线喷水免费视频网站| 晚上一个人看的免费电影| 午夜激情久久久久久久| 九草在线视频观看| 国产精品99久久99久久久不卡 | 国产亚洲av高清不卡| 日本vs欧美在线观看视频| 1024香蕉在线观看| 咕卡用的链子| 青草久久国产| 久久国产精品男人的天堂亚洲| 国产99久久九九免费精品| 天天操日日干夜夜撸| 成年av动漫网址| 久久久精品94久久精品| 久久国产精品大桥未久av| 亚洲熟女毛片儿| e午夜精品久久久久久久| 国产毛片在线视频| 老司机影院成人| 男女之事视频高清在线观看 | 丁香六月天网| 久久综合国产亚洲精品| 欧美激情高清一区二区三区 | 成年女人毛片免费观看观看9 | 满18在线观看网站| 999久久久国产精品视频| 免费av中文字幕在线| 天天添夜夜摸| 久久久久久免费高清国产稀缺| 国产男人的电影天堂91| 亚洲中文av在线| 少妇 在线观看| 一级毛片 在线播放| 国产精品二区激情视频| 曰老女人黄片| 少妇人妻久久综合中文| kizo精华| 亚洲自偷自拍图片 自拍| 观看av在线不卡| 两个人免费观看高清视频| 亚洲成人手机| 最新的欧美精品一区二区| 国产日韩欧美亚洲二区| 黑人欧美特级aaaaaa片| 男女下面插进去视频免费观看| 中文字幕制服av| √禁漫天堂资源中文www| 国产成人精品在线电影| 久久久精品区二区三区| 一级片'在线观看视频| 国产免费现黄频在线看| 麻豆精品久久久久久蜜桃| 在线观看三级黄色| 午夜福利免费观看在线| 丰满饥渴人妻一区二区三| 午夜激情久久久久久久| 女性生殖器流出的白浆| 青春草国产在线视频| 色综合欧美亚洲国产小说| 国产一区二区激情短视频 | 亚洲综合精品二区| 亚洲精品,欧美精品| 午夜影院在线不卡| 热re99久久国产66热| 久久性视频一级片| 在线观看免费视频网站a站| 日韩欧美精品免费久久| 日本色播在线视频| 另类亚洲欧美激情| 日本欧美国产在线视频| 国产精品秋霞免费鲁丝片| 欧美激情 高清一区二区三区|