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

    1960-2020年黃河流域氣候干濕狀況時(shí)空分異及變化趨勢(shì)*

    2022-03-22 10:35:36高秉麗靳甜甜張?jiān)葡?/span>朱月華
    中國(guó)農(nóng)業(yè)氣象 2022年3期
    關(guān)鍵詞:黃河流域狀況時(shí)空

    鞏 杰,高秉麗,李 焱,靳甜甜,張?jiān)葡?,朱月華

    1960-2020年黃河流域氣候干濕狀況時(shí)空分異及變化趨勢(shì)*

    鞏杰,高秉麗,李 焱,靳甜甜,張?jiān)葡?,朱月華

    (蘭州大學(xué)資源環(huán)境學(xué)院/西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,蘭州 730000)

    基于濕潤(rùn)指數(shù),利用黃河流域及周邊1960-2020年113個(gè)氣象站點(diǎn)日值氣候數(shù)據(jù),通過(guò)氣候傾向率、Morlet小波、簡(jiǎn)單相關(guān)分析和Hurst指數(shù)等方法,分析近61a黃河流域干濕狀況時(shí)空演變特征,預(yù)估黃河流域未來(lái)干濕變化趨勢(shì),以期為流域防旱減災(zāi)、合理利用氣候資源提供科學(xué)依據(jù)。結(jié)果表明:(1)空間上,黃河流域年和季節(jié)濕潤(rùn)指數(shù)呈現(xiàn)為東南部高、西北部低。流域年、春秋兩季以變干為主,濕潤(rùn)指數(shù)呈顯著下降的站點(diǎn)占總站點(diǎn)的比例分別為5.31%、7.96%、6.19%;流域夏、冬季以變濕為主,濕潤(rùn)指數(shù)呈顯著上升的站點(diǎn)占比分別為10.62%、13.27%。(2)年際變化上,1960-2020年黃河流域年和季節(jié)濕潤(rùn)指數(shù)均無(wú)顯著變化趨勢(shì)。流域全年、春、夏、秋、冬季分別處于半濕潤(rùn)、半干旱、半濕潤(rùn)、半濕潤(rùn)和半干旱狀態(tài)。黃河流域年和季節(jié)濕潤(rùn)指數(shù)存在多時(shí)間尺度的周期振蕩特征,第一主周期分別為3a、26a、14a、3a和15a。(3)降水量、相對(duì)濕度與濕潤(rùn)指數(shù)呈極顯著正相關(guān),日照時(shí)數(shù)與濕潤(rùn)指數(shù)呈極顯著負(fù)相關(guān),均為黃河流域濕潤(rùn)指數(shù)的主要驅(qū)動(dòng)因素。(4)黃河流域年和季節(jié)濕潤(rùn)指數(shù)的Hurst指數(shù)均大于0.50,表明流域年、春季和秋季氣候?qū)⒊掷m(xù)變干,夏季和冬季氣候?qū)⒊掷m(xù)變濕。

    濕潤(rùn)指數(shù);未來(lái)趨勢(shì);時(shí)空分異;Morlet小波分析;黃河流域

    《氣候變化與土地特別報(bào)告》指出,近百年來(lái)全球氣候日益變暖,2006?2015年全球陸地溫度在1850?1900年基礎(chǔ)上增加了1.53℃[1]。全球氣候變暖及其變化的時(shí)空差異影響著全球和區(qū)域降水、蒸散發(fā)等,進(jìn)而導(dǎo)致全球和區(qū)域氣候的干濕狀況及其空間格局發(fā)生變化,影響農(nóng)牧業(yè)生產(chǎn)和生態(tài)環(huán)境[2?3]。全球和區(qū)域尺度干濕變化特征認(rèn)知是應(yīng)對(duì)和適應(yīng)氣候變化的關(guān)鍵參考因素,備受關(guān)注[4]。干濕狀況可反映地氣系統(tǒng)中水分收支平衡變化[5]。濕潤(rùn)指數(shù)同時(shí)考慮了大氣降水和潛在蒸散發(fā),計(jì)算方便簡(jiǎn)單,且所需資料較少,是目前常用的氣候干濕變化表征指標(biāo),能客觀反映區(qū)域干濕狀況[6],可用于流域水資源綜合管理應(yīng)用[7]。近年來(lái),一些學(xué)者利用濕潤(rùn)指數(shù)開(kāi)展了中國(guó)氣候干濕狀況時(shí)空演變特征研究。如Li等[8]基于濕潤(rùn)指數(shù)發(fā)現(xiàn)中國(guó)淮河流域在年際尺度上呈不顯著變濕趨勢(shì)。Liu等[9]采用濕潤(rùn)指數(shù)發(fā)現(xiàn)1961?2016年中國(guó)西北地區(qū)、青藏高原及華南東部地區(qū)呈濕潤(rùn)化趨勢(shì),而華南西部、華北和東北地區(qū)呈干旱化趨勢(shì)。賈艷青等[10]發(fā)現(xiàn)1960?2016年中國(guó)西北地區(qū)、青藏高原及內(nèi)蒙古地區(qū)呈現(xiàn)濕潤(rùn)化,而華北地區(qū)在變干??梢?jiàn),不同區(qū)域濕潤(rùn)指數(shù)對(duì)氣候變化的響應(yīng)存在一定差異。

    推斷和預(yù)測(cè)未來(lái)氣候變化是氣候研究中的前沿課題。目前,未來(lái)氣候變化預(yù)測(cè)方法主要有兩類:一是從氣候動(dòng)力學(xué)模型出發(fā),二是基于氣候資料和氣候代用資料,利用不同的統(tǒng)計(jì)分析方法推斷未來(lái)氣候變化趨勢(shì)[11]。其中,統(tǒng)計(jì)分析方法中的Hurst指數(shù)可以較好地揭示氣候要素時(shí)序變化趨勢(shì),被廣泛用于不同區(qū)域氣候要素未來(lái)趨勢(shì)預(yù)測(cè)[12]。史雯雨等[13]基于Hurst指數(shù)開(kāi)展了金沙江流域未來(lái)氣溫變化趨勢(shì)預(yù)估;侯青青等[14]利用Hurst指數(shù)預(yù)測(cè)了黃土高原標(biāo)準(zhǔn)化降水蒸散指數(shù)(Standardized Precipitation Evapotranspiration Index,SPEI)的未來(lái)變化趨勢(shì),發(fā)現(xiàn)黃土高原未來(lái)干旱化趨勢(shì)加重。

    黃河流域作為氣候變化敏感區(qū)之一,氣候時(shí)空分異大,加之受人類活動(dòng)等影響,流域降水和蒸散發(fā)的變化趨勢(shì)明顯,流域氣候干濕狀況隨之而變。近年來(lái),國(guó)內(nèi)許多學(xué)者針對(duì)黃河流域開(kāi)展了干旱時(shí)空變化特征方面研究。Wang等[15]采用SPEI發(fā)現(xiàn)黃河流域春、夏、秋季干旱程度均有所加重;Huang等[16]采用綜合干旱指數(shù)指出黃河流域中游旱情最為嚴(yán)重;徐冬平等[17]發(fā)現(xiàn)黃河流域自北向南、自西向東水分虧缺量逐漸遞減;周帥等[18]采用標(biāo)準(zhǔn)化降水指數(shù)(Standardized Precipitation Index,SPI)揭示了黃河流域四季干旱化程度加劇。但已有研究對(duì)黃河流域氣候干濕狀況時(shí)空變化及格局認(rèn)識(shí)還不夠系統(tǒng)深入,對(duì)流域未來(lái)氣候變化趨勢(shì)預(yù)測(cè)報(bào)道相對(duì)較少[19]。因此,本研究基于1960?2020年黃河流域及周邊113個(gè)氣象站點(diǎn)資料,開(kāi)展黃河流域氣候干濕時(shí)空變化特征分析,預(yù)測(cè)流域氣候干濕狀況未來(lái)趨勢(shì),以期為流域水資源利用、農(nóng)牧業(yè)生產(chǎn)及氣候變化應(yīng)對(duì)等提供科學(xué)依據(jù)。

    1 資料與方法

    1.1 研究區(qū)概況

    黃河流域(32°10'?41°50'N、95°53'?119°05'E)總面積約79.5×104km2,從西向東橫跨青藏高原、內(nèi)蒙古高原、黃土高原和黃淮海平原,包括青海、四川、甘肅、寧夏、內(nèi)蒙古、陜西、山西、河南及山東9省區(qū)(圖1),選取各省區(qū)在黃河流域的區(qū)域作為研究區(qū)[20]。黃河流域地形地貌差異顯著,地勢(shì)西高東低,海拔介于?9~6250m。全流域多年平均降水量440mm,多年平均氣溫7℃,降水和氣溫均呈由南向北、由東向西遞減的特征[21]。

    圖1 黃河流域及相關(guān)氣象站點(diǎn)分布

    1.2 數(shù)據(jù)來(lái)源

    選取黃河流域及周邊氣象站點(diǎn)1960?2020年的逐日氣象數(shù)據(jù),數(shù)據(jù)來(lái)源于中國(guó)氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/),包括最高氣溫、平均氣溫、最低氣溫、平均風(fēng)速、日照時(shí)數(shù)、相對(duì)濕度和降水量等。受建站時(shí)間、儀器損壞等影響,部分站點(diǎn)存在數(shù)據(jù)缺測(cè)(失)。研究中設(shè)定:1月內(nèi)日值缺測(cè)多于7d則為缺測(cè)月,1年內(nèi)月值缺測(cè)多于1月則為缺測(cè)年[22]。剔除缺測(cè)月、年值及觀測(cè)年限不足61a(1960?2020年)的站點(diǎn)后,以數(shù)據(jù)較齊全的113個(gè)站點(diǎn)為基礎(chǔ)數(shù)據(jù),利用相鄰日該氣象要素的平均值對(duì)零星缺測(cè)值進(jìn)行插補(bǔ)[23]。以3?5月為春季,6?8月為夏季,9?11月為秋季,12月?翌年2月為冬季進(jìn)行分析。

    1.3 研究方法

    1.3.1濕潤(rùn)指數(shù)

    濕潤(rùn)指數(shù)是區(qū)域降水量和潛在蒸散發(fā)的比值,可客觀反映區(qū)域水熱平衡狀況,是一個(gè)理想的能夠表征地表干濕狀況的指標(biāo)[24]。計(jì)算式為

    式中,HI為濕潤(rùn)指數(shù);R為降水量(mm);PET為潛在蒸散發(fā)(mm)。各站點(diǎn)的潛在蒸散發(fā)(PET)采用Penman-Monteith模型計(jì)算,即

    1.3.2 Hurst指數(shù)

    累積離差為

    極差為

    標(biāo)準(zhǔn)差為

    1.3.3 其它方法

    2 結(jié)果與分析

    2.1 黃河流域氣候干濕狀況的空間分布格局

    2.1.1 年/季濕潤(rùn)指數(shù)空間分布

    1960?2020年黃河流域多年平均濕潤(rùn)指數(shù)介于0.15~1.31,總體呈現(xiàn)自東南向西北逐漸遞減的態(tài)勢(shì)(圖2a)。處于干旱狀態(tài)的區(qū)域占流域總面積的3.39%,主要為內(nèi)蒙古西部和寧夏東北部。處于半干旱狀態(tài)的區(qū)域占比為39.54%,主要為青海北部、甘肅北部、寧夏中北部、陜西西北部及內(nèi)蒙古大部。處于半濕潤(rùn)狀態(tài)的區(qū)域占比為56.55%,主要位于流域西南部及東南部;處于濕潤(rùn)狀態(tài)的區(qū)域占比僅為0.52%,主要分布于黃河下游的山東片區(qū)(圖2a)。由圖2b可知,黃河流域年降水量自東南向西北遞減,與流域年濕潤(rùn)指數(shù)空間分布格局相似,表明流域年濕潤(rùn)指數(shù)空間變化對(duì)流域年降水量變化極為敏感[6],且流域年濕潤(rùn)指數(shù)0.2、0.5等值線分別與200mm、400mm等降水量線基本吻合(圖2b)。黃河流域中東部為潛在蒸散發(fā)高值區(qū),潛在蒸散發(fā)低值區(qū)主要位于流域西南部,潛在蒸散發(fā)與濕潤(rùn)指數(shù)的空間分布相似性較弱,可見(jiàn),相較于潛在蒸散發(fā),黃河流域濕潤(rùn)指數(shù)對(duì)降水量的敏感性更高(圖2c)。

    圖2 黃河流域年濕潤(rùn)指數(shù)、降水量和潛在蒸散發(fā)的空間分布(1960-2020年)

    黃河流域季節(jié)濕潤(rùn)指數(shù)與年濕潤(rùn)指數(shù)的空間分布格局相似,流域各季濕潤(rùn)指數(shù)均為自東南向西北遞減(圖3)。黃河流域春、夏、秋、冬季濕潤(rùn)指數(shù)范圍分別介于0.07~0.92、0.26~2.98、0.19~1.86及0.027~0.48,平均濕潤(rùn)指數(shù)分別為0.34、0.86、0.77及0.19。春季流域處于干旱狀態(tài)的區(qū)域占比為20.82%,主要位于內(nèi)蒙古西部、寧夏北部、甘肅東北角等;處于半干旱狀態(tài)的區(qū)域占比為63.08%,主要分布于東西走向的流域中部地帶及流域東北部等;處于半濕潤(rùn)狀態(tài)的區(qū)域占比為16.10%,主要為四川大部、甘肅西南部、陜西西南部和河南西南部等(圖3a)。夏季流域處于半干旱、半濕潤(rùn)和濕潤(rùn)狀態(tài)的區(qū)域占比分別為14.47%、58.29%和27.24%(圖 3b)。夏季處于半干旱狀態(tài)的區(qū)域主要位于內(nèi)蒙古西部、寧夏北部、甘肅東北部等;處于半濕潤(rùn)狀態(tài)的片區(qū)主要有青海北部、甘肅局部、寧夏南部、陜西大部、山西中西部、內(nèi)蒙古中東部及河南西北部等;處于濕潤(rùn)狀態(tài)的區(qū)域主要為流域西南部及流域東部等。秋季流域處于半濕潤(rùn)狀態(tài)的區(qū)域占比超過(guò)流域面積的一半,主要位于青海大部、甘肅局部、寧夏南部、陜西北部、內(nèi)蒙古東部和山西北部等(圖 3c)。冬季流域處于干旱狀態(tài)的區(qū)域占比為59.84%,主要位于青海和甘肅大部、寧夏、內(nèi)蒙古中西部及陜西西北部等(圖3d)。綜上,黃河流域年、夏季和秋季處于半濕潤(rùn)狀態(tài)的區(qū)域面積占比最大,春季時(shí)流域多為半干旱狀態(tài),而冬季多為干旱狀態(tài)(圖3)。

    2.1.2 年/季濕潤(rùn)指數(shù)變化趨勢(shì)

    圖3 黃河流域季節(jié)濕潤(rùn)指數(shù)空間分布(1960?2020年)

    Fig. 3 Spatial distribution of seasonal humidity index in the Yellow River Basin(1960?2020)

    圖4 1960-2020年黃河流域年/季濕潤(rùn)指數(shù)變化傾向率[(10a)?1]的空間分布

    Fig. 4 Spatial distribution of annual and seasonal climatic variables [(10y)?1] of humidity index in the Yellow River Basin in 1960-2020

    2.2 黃河流域氣候干濕狀況時(shí)間變化特征

    2.2.1 流域平均濕潤(rùn)指數(shù)變化趨勢(shì)

    由圖5可見(jiàn),1960?2020年黃河流域夏季和秋季濕潤(rùn)指數(shù)偏高,而春季和冬季的濕潤(rùn)指數(shù)偏低,年濕潤(rùn)指數(shù)則介于其間,年和四季濕潤(rùn)指數(shù)均無(wú)顯著變化趨勢(shì)。研究期內(nèi)黃河流域平均濕潤(rùn)指數(shù)為0.57(半濕潤(rùn)狀態(tài)),年濕潤(rùn)指數(shù)最大值出現(xiàn)在1964年(0.92);年濕潤(rùn)指數(shù)最小值出現(xiàn)在1986年(0.41)。春、夏、秋和冬季平均濕潤(rùn)指數(shù)分別為0.36、0.92、0.80和0.21,表明四季分別處于半干旱、半濕潤(rùn)、半濕潤(rùn)和半干旱狀態(tài)(圖5)。

    圖5 1960?2020年黃河流域平均年/季濕潤(rùn)指數(shù)變化過(guò)程

    2.2.2 流域平均濕潤(rùn)指數(shù)變化周期

    由圖6可見(jiàn),1960?2020年黃河流域年際和季節(jié)濕潤(rùn)指數(shù)存在多重時(shí)間尺度的振蕩和周期變化。黃河流域年濕潤(rùn)指數(shù)在2~5a、6~11a和20~32a變化較為密集(圖6a),流域存在一個(gè)3a的主周期和7a、9a、14a、25a的次周期,其中3a經(jīng)歷了密集的“干?濕”變化循環(huán)(圖6f)。流域春季濕潤(rùn)指數(shù)在2~14a、18~32a的時(shí)段較為密集(圖6b),存在26a、4a、6a和12a四個(gè)周期,其中26a經(jīng)歷了“濕?干?濕?干?濕”的交替循環(huán)(圖6g)。流域夏季濕潤(rùn)指數(shù)在2~5a、6~9a和11~18a的時(shí)間尺度上較為密集(圖6c);夏季濕潤(rùn)指數(shù)在14a時(shí)間周期的小波方差值最大(即主周期),并在該周期經(jīng)歷了4次“濕?干”交替循環(huán),此外在3a、8a時(shí)間周期也存在極值(即次周期)(圖6h)。流域秋季濕潤(rùn)指數(shù)在2~5a、6~12a和17~26a的時(shí)間尺度上較為密集(圖6d),小波方差圖顯示3a、9a、18a、25a處存在極值,可推斷流域秋季濕潤(rùn)指數(shù)存在一個(gè)3a的主周期和9a、18a、25a的次周期(圖6i)。黃河流域冬季濕潤(rùn)指數(shù)在2~7a、10~22a的時(shí)間尺度上較為密集(圖6e);由圖6j可見(jiàn),冬季存在一個(gè)15a的主周期和3a、4a、6a的次周期,其中15a經(jīng)歷了4次“干?濕”變化循環(huán)。綜上,黃河流域年濕潤(rùn)指數(shù)存在2~5a、6~11a、20~32a的周期變化,流域年和四季濕潤(rùn)指數(shù)分別存在3a、26a、14a、3a和15a的第一主周期(圖6)。

    圖6 黃河流域年/季濕潤(rùn)指數(shù)小波變換系數(shù)實(shí)部等值線圖和小波方差圖

    2.3 黃河流域氣候干濕狀況的影響因素分析

    2.4 黃河流域未來(lái)干濕狀況趨勢(shì)預(yù)測(cè)

    表1 1960?2020年黃河流域濕潤(rùn)指數(shù)與各氣象因子的相關(guān)系數(shù)

    注:*和**分別表示相關(guān)系數(shù)通過(guò)0.05、0.01水平的顯著性檢驗(yàn)。

    圖7 黃河流域年/季濕潤(rùn)指數(shù)序列的重標(biāo)極差(R/S)分析結(jié)果

    3 結(jié)論與討論

    3.1 結(jié)論

    (1)空間上,1960?2020年黃河流域年和四季濕潤(rùn)指數(shù)呈現(xiàn)東南高西北低的分布特征。流域全年、春、秋季以變干為主,濕潤(rùn)指數(shù)呈顯著下降的站點(diǎn)占比分別為5.31%、7.96%、6.19%;流域夏、冬季以變濕為主,濕潤(rùn)指數(shù)呈顯著上升的站點(diǎn)占比分別為10.62%和13.27%。

    (2)年際變化上,1960?2020年黃河流域年和季節(jié)濕潤(rùn)指數(shù)均無(wú)顯著變化趨勢(shì)。流域全年、春、夏、秋和冬季分別處于半濕潤(rùn)、半干旱、半濕潤(rùn)、半濕潤(rùn)和半干旱狀態(tài)。黃河流域濕潤(rùn)指數(shù)在年際和季節(jié)變化上具有多個(gè)時(shí)間尺度的周期振蕩特征,第一主周期分別為3a、26a、14a、3a和15a。

    (3)黃河流域年和各季濕潤(rùn)指數(shù)均與降水量、相對(duì)濕度呈極顯著正相關(guān),與日照時(shí)數(shù)呈極顯著負(fù)相關(guān),降水量、相對(duì)濕度和日照時(shí)數(shù)是影響黃河流域濕潤(rùn)指數(shù)變化的主導(dǎo)因素,最高氣溫、平均氣溫、最低氣溫和平均風(fēng)速對(duì)濕潤(rùn)指數(shù)影響較小。

    (4)黃河流域年和季節(jié)干濕狀況的未來(lái)變化趨勢(shì)與過(guò)去61a干濕狀況變化趨勢(shì)有著良好的自相似性,即黃河流域年、春季和秋季在未來(lái)一段時(shí)間內(nèi)氣候?qū)⒊掷m(xù)變干,其中秋季呈較強(qiáng)持續(xù)性變干;夏季和冬季流域氣候?qū)⒊掷m(xù)變濕。

    3.2 討論

    氣候干濕狀況是表征區(qū)域氣候變化的重要指標(biāo)[31]。黃河流域年和四季濕潤(rùn)指數(shù)空間分布不均,受地形、海陸距離、不同氣象因子等的影響,流域年和四季濕潤(rùn)指數(shù)自東南向西北內(nèi)陸逐漸遞減,這與劉勤等[32]的研究結(jié)果一致。黃河流域春季以變干為主,這可能是由于春季降水較少而春播和植被返青需水耗水增加,且隨著氣溫回升等導(dǎo)致潛在蒸散發(fā)增大[33]。流域秋季以變干為主,這可能是由于該季節(jié)降水減少趨勢(shì)大于潛在蒸散發(fā)減少趨勢(shì),且秋季干旱頻率較高,干旱歷時(shí)長(zhǎng)及干旱烈度大[18]。夏季流域以變濕為主即旱情減緩,可能是由于該季節(jié)氣溫升高且降水豐沛,降水量增加在很大程度上消減了潛在蒸散發(fā)增大等造成的影響。流域冬季受干冷氣團(tuán)控制,溫度極低,潛在蒸散發(fā)小,且伴隨降水的增加將導(dǎo)致流域以變濕為主。Hurst指數(shù)表明,黃河流域年、春季和秋季在未來(lái)一段時(shí)間內(nèi)氣候?qū)⒊掷m(xù)變干,而夏季和冬季氣候?qū)⒊掷m(xù)變濕。Huang等[16]研究也表明黃河流域未來(lái)氣溫將明顯升高,可能導(dǎo)致流域水資源減少和旱情加重。綜合考慮黃河流域濕潤(rùn)指數(shù)時(shí)空變化特征、未來(lái)趨勢(shì)及農(nóng)牧業(yè)與社會(huì)經(jīng)濟(jì)發(fā)展態(tài)勢(shì),黃河流域氣候變化應(yīng)對(duì)和水資源利用的對(duì)策建議主要有,黃河流域絕大部分地區(qū)年濕潤(rùn)指數(shù)小于1,季節(jié)性干旱特征明顯,生態(tài)環(huán)境脆弱,需遵循開(kāi)發(fā)利用和建設(shè)保護(hù)并舉原則;春季和秋季是黃河流域抗旱減災(zāi)的重點(diǎn)時(shí)段,采取節(jié)水灌溉、種植耐寒(旱)作物、加強(qiáng)農(nóng)田水利基本建設(shè)、水資源的科學(xué)規(guī)劃與合理利用等措施減緩黃河流域未來(lái)干旱趨勢(shì)。同時(shí),黃河流域不同區(qū)域干濕狀況變化趨勢(shì)不一致,應(yīng)因地制宜,進(jìn)行水資源合理利用和旱澇防災(zāi)減災(zāi)等。此外,本研究側(cè)重于黃河流域干濕狀況的趨勢(shì)預(yù)測(cè),對(duì)流域未來(lái)干濕變化趨勢(shì)的空間差異分析不夠。今后應(yīng)進(jìn)一步分析黃河流域未來(lái)干濕變化趨勢(shì)的空間差異及其成因。

    本研究利用濕潤(rùn)指數(shù)分析了近61a黃河流域干濕狀況時(shí)空變化,可為黃河流域農(nóng)牧業(yè)生產(chǎn)、水資源合理利用及防災(zāi)減災(zāi)等提供一定科學(xué)支撐。但濕潤(rùn)指數(shù)僅是一個(gè)氣象指標(biāo),一個(gè)區(qū)域干濕狀況研究不僅要考慮氣象要素的變化,還應(yīng)考慮地形地勢(shì)、人類活動(dòng)、作物種類、土壤性質(zhì)等非氣象要素的變化。因此,氣候變化背景下的濕潤(rùn)指數(shù)、地形地貌、人類活動(dòng)等對(duì)黃河流域干濕狀況的共同影響是未來(lái)亟待解決的議題之一。

    [1] 賈根鎖.IPCC《氣候變化與土地特別報(bào)告》對(duì)陸氣相互作用的新認(rèn)知[J].氣候變化研究進(jìn)展,2020,16(1):9-16.

    Jia G S.New understanding of land-climate interactions from IPCC special report on climate change and land[J]. Climate Change Research,2020,16(1):9-16.(in Chinese)

    [2] 陳潔,劉玉潔,潘韜,等.1961-2010年中國(guó)降水時(shí)空變化特征及對(duì)地表干濕狀況影響[J].自然資源學(xué)報(bào),2019,34(11): 2440-2453.

    Chen J,Liu Y J,Pan T,et al.Spatiotemporal variation of precipitation in China and its impact on surface dry-wet conditions during 1961-2010[J].Journal of Natural Resources, 2019,34(11):2440-2453.(in Chinese)

    [3] 王慶莉,韓玉江,郭斌,等.近57年甘孜州氣候干濕狀況的時(shí)空演變[J].中國(guó)農(nóng)業(yè)氣象,2019,40(7):435-443.

    Wang Q L,Han Y J,Guo B,et al.Temporal and spatial evolution of climate dry and wet conditions in Ganzi in the past 57 years[J].Chinese Journal of Agrometeorology, 2019, 40(7):435-443.(in Chinese)

    [4] 劉珂,姜大膀.RCP4.5情景下中國(guó)未來(lái)干濕變化預(yù)估[J].大氣科學(xué),2015,39(3):489-502.

    Liu K,Jiang D B.Projected changes in the dry/wet climate of China under the RCP4.5 scenario[J].Chinese Journal of Atmospheric Sciences,2015,39(3):489-502.(in Chinese)

    [5] 張波,曹華,吳戰(zhàn)平,等.貴州氣候干濕狀況及影響因素分析[J].氣象與環(huán)境科學(xué),2021,44(3):39-46.

    Zhang B,Cao H,Wu Z P,et al.Analysis on the status of dry-wet climate and the influence factors in Guizhou[J]. Meteorological and Environmental Sciences,2021,44(3): 39-46.(in Chinese)

    [6] 王允,劉普幸,曹立國(guó),等.基于濕潤(rùn)指數(shù)的1960-2011年中國(guó)西南地區(qū)地表干濕變化特征[J].自然資源學(xué)報(bào),2014, 29(5):830-838.

    Wang Y,Liu P X,Cao L G,et al.Characteristics of southwestern China dry-wet condition based on wetness index in 1960-2011[J].Journal of Natural Resources,2014, 29(5):830-838.(in Chinese)

    [7] Zhang K X,Pan S M,Zhang W,et al.Influence of climate change on reference evapotranspiration and aridity index and their temporal-spatial variations in the Yellow River Basin,China,from 1961 to 2012[J].Quaternary International, 2015,380-381:75-82.

    [8] Li M,Chu R H,Islam ARMT,et al.Attribution analysis of long-term trends of aridity index in the Huai River Basin,Eastern China[J].Sustainability,2020,12(5):1743.

    [9] Liu L B,Wang Y,You N S,et al.Changes in aridity and its driving factors in China during 1961-2016[J].International Journal of Climatology,2019,39(1):50-60.

    [10] 賈艷青,張勃.近57年中國(guó)北方氣候干濕變化及與太平洋年代際振蕩的關(guān)系[J].土壤學(xué)報(bào),2019,56(5):1085-1097.

    Jia Y Q,Zhang B.Relationship of dry-wet climate changes in Northern China in the past 57 years with Pacific Decadal Oscillation(PDO)[J].Acta Pedologica Sinica,2019,56(5): 1085-1097.(in Chinese)

    [11] 馮新靈,羅隆誠(chéng),馮自立.中國(guó)近50年降水變化趨勢(shì)及突變的Hurst指數(shù)試驗(yàn)[J].干旱區(qū)地理,2009, 32(6):859-866.

    Feng X L,Luo L C,Feng Z L.Hurst index experiment on precipitation change trend and mutation of China in the near 50 years[J].Arid Land Geography,2009,32(6):859- 866.(in Chinese)

    [12] 姜麗霞,王晾晾,呂佳佳,等.基于Hurst指數(shù)的黑龍江省作物生長(zhǎng)季降水趨勢(shì)研究[J].氣象與環(huán)境學(xué)報(bào), 2020,36(2): 70-77.

    Jiang L X,Wang L L,Lv J J,et al.A study on rainfall trend during the crop growing season in Heilongjiang province based on Hurst exponent method[J].Journal of Meteorology and Environment,2020,36(2):70-77.(in Chinese)

    [13] 史雯雨,楊勝勇,李增永,等.近57年金沙江流域氣溫變化特征及未來(lái)趨勢(shì)預(yù)估[J].水土保持研究,2021, 28(1):211- 217.

    Shi W Y, Yang S Y, Li Z Y, et al. Variation characteristics and the future trend estimation of temperature in Chinsha River Basin over the past 57 years[J]. Research of Soil and Water Conservation,2021,28(1): 211-217. (in Chinese)

    [14] 侯青青,裴婷婷,陳英,等.1986-2019年黃土高原干旱變化特征及趨勢(shì)[J].應(yīng)用生態(tài)學(xué)報(bào),2021,32(2):649- 660.

    Hou Q Q,Pei T T,Chen Y,et al.Variations of drought and its trend in the Loess Plateau from 1986 to 2019[J].Chinese Journal of Applied Ecology,2021,32(2):649-660.(in Chinese)

    [15] Wang F,Wang Z M,Yang H B,et al.Copula-based droughtanalysis using Standardized Precipitation Evapotranspiration Index:a case study in the Yellow River Basin,China[J]. Water,2019,11(6):1298.

    [16] Huang S Z,Chang J X,Leng G Y,et al.Integrated index for drought assessment based on variable fuzzy set theory:a case study in the Yellow River Basin,China[J].Journal of Hydrology,2015,527:608-618.

    [17] 徐冬平,卜金洋,楊軍輝.黃河流域水分盈虧量時(shí)空變化及其影響因素[J].干旱區(qū)研究,2016,33(4):739-746.

    Xu D P,Bu J Y,Yang J H.Spatiotemporal change of water budget and its affecting factors in the Yellow River Basin[J].Arid Zone Research,2016,33(4):739-746.(in Chinese)

    [18] 周帥,王義民,暢建霞,等.黃河流域干旱時(shí)空演變的空間格局研究[J].水利學(xué)報(bào),2019,50(10):1231-1241.

    Zhou S,Wang Y M,Chang J X,et al.Research on spatio-temporal evolution of drought patterns in the Yellow River Basin[J].Journal of Hydraulic Engineering,2019, 50(10):1231-1241.(in Chinese)

    [19] 黎云云.氣候和土地利用變化下流域干旱評(píng)估-傳播-驅(qū)動(dòng)-預(yù)測(cè)研究[D].西安:西安理工大學(xué),2018.

    Li Y Y.Research on the drought assessment-propagation- driving-prediction under the climate and land use land cover change scenarios[D].Xi’an:Xi’an University of Technology,2018.(in Chinese)

    [20] 張志強(qiáng),劉歡,左其亭,等.2000-2019年黃河流域植被覆蓋度時(shí)空變化[J].資源科學(xué),2021,43(4):849-858.

    Zhang Z Q,Liu H,Zuo Q T,et al.Spatiotemporal change of fractional vegetation cover in the Yellow River Basin during 2000-2019[J].Resources Science,2021,43(4):849- 858.(in Chinese)

    [21] 楊澤康,田佳,李萬(wàn)源,等.黃河流域生態(tài)環(huán)境質(zhì)量時(shí)空格局與演變趨勢(shì)[J].生態(tài)學(xué)報(bào),2021,41(19): 7627-7636.

    Yang Z K,Tian J,Li W Y,et al.Spatio-temporal pattern and evolution trend of ecological environment quality in the Yellow River Basin[J].Acta Ecologica Sinica,2021,41(19): 7627-7636.(in Chinese)

    [22] 范帥邦,肖春柳,曹永強(qiáng),等.1964-2019年遼寧省平均風(fēng)速時(shí)空演變特征及其影響因素[J].地理科學(xué), 2021,41(4): 717-727.

    Fan S B,Xiao C L,Cao Y Q,et al.Multiscale analysis of wind speed and its influencing factors in Liaoning province from 1964 to 2019[J].Scientia Geographica Sinica,2021, 41(4):717-727.(in Chinese)

    [23] 吳霞,王培娟,霍治國(guó),等.1961-2015年中國(guó)潛在蒸散時(shí)空變化特征與成因[J].資源科學(xué),2017,39(5): 964-977.

    Wu X,Wang P J,Huo Z G,et al.Spatio-temporal distribution characteristics of potential evapotranspiration and impact factors in China from 1961 to 2015[J].Resources Science, 2017,39(5):964-977.(in Chinese)

    [24] 胡琦,董蓓,潘學(xué)標(biāo),等.不同時(shí)間尺度下華北平原干濕氣候時(shí)空變化及成因分析[J].中國(guó)農(nóng)業(yè)氣象, 2017,38(5): 267-277.

    Hu Q,Dong B,Pan X B,et al.Spatiotemporal variation and causes analysis of dry-wet climate at different time scales in North China Plain[J].Chinese Journal of Agrometeorology, 2017,38(5):267-277.(in Chinese)

    [25] Allen R G,Pereira L S,Raes D,et al.Crop evapotranspiration: guidelines for computing crop water requirements[M] Rome:United Nations Food and Agriculture Organization, 1998.

    [26] Garcia M,Raes D,Allen R,et al.Dynamics of reference evapotranspiration in the Bolivian highlands (Altiplano)[J]. Agricultural and Forest Meteorology,2004,125(1-2):1-82.

    [27] 劉鈺,蔡林根.參照騰發(fā)量的新定義及計(jì)算方法對(duì)比[J].水利學(xué)報(bào),1997(6):28-34.

    Liu Y,Cai L G.Update definition and computation of reference evapotranspiration comparison with former method[J].Journal of Hydraulic Engineering,1997(6):28- 34.(in Chinese)

    [28] 杜軍,周刊社,高佳佳,等.未來(lái)氣候變暖對(duì)羌塘自然保護(hù)區(qū)高寒草地牧草青草期的影響[J].中國(guó)農(nóng)業(yè)氣象, 2021, 42(10):845-858.

    Du J,Zhou K S,Gao J J,et al.Influence of climate warming in future on the green grass season for alpine grassland in Qiangtang Nature Reserve of Tibet[J].Chinese Journal of Agrometeorology,2021,42(10): 845-858.(in Chinese)

    [29] 王炎強(qiáng),趙軍,李忠勤,等.1977-2017年薩吾爾山冰川變化及其對(duì)氣候變化的響應(yīng)[J].自然資源學(xué)報(bào),2019,34(4): 802-814.

    Wang Y Q,Zhao J,Li Z Q,et al.Glacier changes in the Sawuer Mountain during 1977-2017 and their response to climate change[J].Journal of Natural Resources,2019,34(4): 802-814.(in Chinese)

    [30] 張景華,封志明,姜魯光,等.瀾滄江流域植被NDVI與氣候因子的相關(guān)性分析[J].自然資源學(xué)報(bào),2015, 30(9):1425- 1435.

    Zhang J H,Feng Z M,Jiang L G,et al.Analysis of the correlation between NDVI and climate factors in the Lancang River Basin[J].Journal of Natural Resources,2015, 30(9):1425-1435.(in Chinese)

    [31] 李秀芬,郭昭濱,趙慧穎,等.大興安嶺氣候干濕變化及對(duì)森林火災(zāi)的影響[J].應(yīng)用氣象學(xué)報(bào),2018, 29(5):619-629.

    Li X F,Guo Z B,Zhao H Y,et al.Change of dry and wet climate and its influence on forest fires in the Great Xing’an Mountains[J].Journal of Applied Meteorological Science,2018,29(5):619-629.(in Chinese)

    [32] 劉勤,嚴(yán)昌榮,何文清.黃河流域干旱時(shí)空變化特征及其氣候要素敏感性分析[J].中國(guó)農(nóng)業(yè)氣象,2016,37(6): 623- 632.

    Liu Q,Yan C R,He W Q.Drought variation and its sensitivity coefficients to climatic factors in the Yellow River Basin[J].Chinese Journal of Agrometeorology, 2016,37(6):623-632.(in Chinese)

    [33] 史建國(guó),嚴(yán)昌榮,何文清,等.黃河流域水分虧缺時(shí)空格局變化研究[J].自然資源學(xué)報(bào),2008(1):113-119.

    Shi J G, Yan C R, He W Q, et al. Study on spatial and temporal variation of water surplus and deficiency in Yellow River Basin[J]. Journal of Natural Resources, 2008(1):113-119.(in Chinese)

    Spatiotemporal Variation of Climate Dry-Wet Condition and Its Potential Trend in the Yellow River Basin from 1960 to 2020

    GONG Jie, GAO Bing-li, LI Yan, JIN Tian-tian, ZHANG Yun-xia, ZHU Yue-hua

    (College of Earth and Environmental Sciences/Key Laboratory of Western China’s Environmental Systems, Ministry of Education, Lanzhou University, Lanzhou 730000, China)

    Based on the humidity index, using the daily climate data of 113 meteorological stations in and adjacent the Yellow River Basin(YRB) from 1960 to 2020, climate trend rate, Morlet wavelet, simple correlation analysis, and Hurst index were used to analysis the spatiotemporal characteristics of dry-wet condition in the YRB in the latest 61 years and to predict the future trend of dry-wet changes in the YRB, to provide scientific basis for drought prevention and mitigation and rational utilization of climate resources in the YRB. The results showed that: (1) among all the geographic zones, the annual and seasonal humidity index in the YRB was high in the southeastern YRB and low in the northwestern YRB. Spatially, the climate of the annual, spring and autumn in the YRB was mainly dry, the proportions of stations with significant decrease of humidity index of the total stations were 5.31%, 7.96% and 6.19%, respectively.The climate of summer and winter was mainly wet, the proportions of stations with significant increase of humidity index were 10.62% and 13.27%, respectively. (2) As for the interannual change, the annual and seasonal humidity index of the YRB from 1960 to 2020 showed no significant change. Semi-humid, semi-drought, semi-humid, semi-humid, semi-drought were found in entire year, spring, summer, autumn, and winter, respectively. The annual and seasonal humidity index of the YRB had periodic oscillation characteristics at multiple time scales, and the first main period was 3 years, 26 years, 14 years, 3 years and 15 years, respectively. (3) The change of humidity index in the YRB was extremely significantly positively correlated with precipitation and relative humidity. There was a extremely significant negative correlation between sunshine duration and humidity index. Precipitation, relative humidity, and sunshine duration were the main factors affecting the change of humidity index in the YRB. (4) The Hurst index of annual and seasonal humidity index in the YRB were larger than 0.50, indicating that the climate of annual, spring, and autumn of the YRB would continue to be drying, the summer and winter climate would continue to be wetting.

    Humidity index; Future trend; Spatiotemporal variation; Morlet wavelet analysis; Yellow River Basin

    鞏杰,高秉麗,李焱,等.1960-2020年黃河流域氣候干濕狀況時(shí)空分異及變化趨勢(shì)[J].中國(guó)農(nóng)業(yè)氣象,2022,43(3):165-176

    10.3969/j.issn.1000-6362.2022.03.001

    2021?09?02

    國(guó)家自然科學(xué)基金重大項(xiàng)目(41991231)

    鞏杰,博士,教授,主要從事景觀生態(tài)學(xué)、生態(tài)系統(tǒng)服務(wù)、環(huán)境遙感與生態(tài)評(píng)價(jià)等方面的科研和教學(xué)工作,E-mail: jgong@lzu.edu.cn

    猜你喜歡
    黃河流域狀況時(shí)空
    生態(tài)環(huán)境部啟動(dòng)新一年度黃河流域“清廢行動(dòng)”
    跨越時(shí)空的相遇
    聲敏感患者的焦慮抑郁狀況調(diào)查
    鏡中的時(shí)空穿梭
    黃河流域燦爛的齊家文化——陶器
    收藏界(2019年2期)2019-10-12 08:26:10
    2019年中國(guó)國(guó)際收支狀況依然會(huì)保持穩(wěn)健
    增強(qiáng)大局意識(shí) 提升黃河流域生態(tài)保護(hù)發(fā)展水平
    玩一次時(shí)空大“穿越”
    第五節(jié) 2015年法學(xué)專業(yè)就業(yè)狀況
    三十六計(jì)之順手牽羊
    午夜免费鲁丝| 亚洲中文字幕日韩| 一区二区三区乱码不卡18| 亚洲国产中文字幕在线视频| 别揉我奶头~嗯~啊~动态视频 | 午夜免费男女啪啪视频观看| 嫩草影视91久久| 久久久久国产精品人妻一区二区| 国产成人精品久久二区二区免费| 国产日韩欧美亚洲二区| 国产av国产精品国产| www.精华液| 少妇人妻 视频| 一区二区三区精品91| 久久中文字幕一级| 国产一卡二卡三卡精品| 精品人妻一区二区三区麻豆| 啦啦啦啦在线视频资源| 99国产精品一区二区三区| 精品第一国产精品| 一本—道久久a久久精品蜜桃钙片| 日本av免费视频播放| 亚洲国产精品一区二区三区在线| 十八禁网站网址无遮挡| 欧美97在线视频| 91成人精品电影| 午夜福利视频在线观看免费| 老司机深夜福利视频在线观看 | 久久久国产一区二区| av福利片在线| 又大又爽又粗| 性色av乱码一区二区三区2| 99re6热这里在线精品视频| 久热这里只有精品99| 日本五十路高清| 中文字幕亚洲精品专区| 久久精品亚洲熟妇少妇任你| 精品国产乱码久久久久久小说| 一级,二级,三级黄色视频| 亚洲精品美女久久久久99蜜臀 | 99国产综合亚洲精品| 男人舔女人的私密视频| 在线观看www视频免费| 777久久人妻少妇嫩草av网站| 色播在线永久视频| 国产精品一区二区在线不卡| 欧美日韩亚洲国产一区二区在线观看 | www日本在线高清视频| 啦啦啦在线免费观看视频4| 国产片内射在线| 中文字幕色久视频| 欧美日韩福利视频一区二区| 1024视频免费在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产欧美日韩精品亚洲av| 一区二区日韩欧美中文字幕| 少妇精品久久久久久久| 啦啦啦啦在线视频资源| 久热爱精品视频在线9| 亚洲欧美清纯卡通| 王馨瑶露胸无遮挡在线观看| 黄色a级毛片大全视频| 欧美国产精品一级二级三级| 久久久久久亚洲精品国产蜜桃av| 美女主播在线视频| 国产精品偷伦视频观看了| 男女之事视频高清在线观看 | 亚洲激情五月婷婷啪啪| 欧美激情高清一区二区三区| 大香蕉久久网| 宅男免费午夜| 成人午夜精彩视频在线观看| 搡老岳熟女国产| 嫩草影视91久久| 在线亚洲精品国产二区图片欧美| 秋霞在线观看毛片| 手机成人av网站| 夜夜骑夜夜射夜夜干| 男男h啪啪无遮挡| 亚洲色图综合在线观看| 久久综合国产亚洲精品| av天堂久久9| 久久ye,这里只有精品| av网站在线播放免费| 国产97色在线日韩免费| 激情五月婷婷亚洲| 操美女的视频在线观看| 91国产中文字幕| 亚洲人成电影观看| 国产黄色免费在线视频| cao死你这个sao货| cao死你这个sao货| 久久国产精品大桥未久av| 成人亚洲精品一区在线观看| 亚洲国产日韩一区二区| 亚洲国产精品999| 亚洲人成网站在线观看播放| 久热这里只有精品99| 亚洲国产中文字幕在线视频| 女警被强在线播放| 99热国产这里只有精品6| 观看av在线不卡| 国产精品成人在线| 国产精品偷伦视频观看了| 欧美亚洲 丝袜 人妻 在线| 国产精品 国内视频| 久久久亚洲精品成人影院| 成年美女黄网站色视频大全免费| 欧美激情高清一区二区三区| 成人亚洲欧美一区二区av| 岛国毛片在线播放| 国产亚洲欧美在线一区二区| 国产一区二区激情短视频 | 考比视频在线观看| 蜜桃国产av成人99| 久久久久久久国产电影| 国产精品一二三区在线看| 一级a爱视频在线免费观看| 又黄又粗又硬又大视频| 超碰成人久久| 一本大道久久a久久精品| 9191精品国产免费久久| 真人做人爱边吃奶动态| 亚洲av片天天在线观看| 欧美久久黑人一区二区| 两个人看的免费小视频| 91老司机精品| 日日夜夜操网爽| 脱女人内裤的视频| 亚洲欧美日韩高清在线视频 | 欧美变态另类bdsm刘玥| 宅男免费午夜| 精品少妇黑人巨大在线播放| 啦啦啦 在线观看视频| 欧美激情极品国产一区二区三区| 日韩,欧美,国产一区二区三区| 一二三四社区在线视频社区8| 尾随美女入室| 看免费成人av毛片| 大香蕉久久网| 精品久久久久久久毛片微露脸 | 精品亚洲成a人片在线观看| 欧美精品高潮呻吟av久久| 天天影视国产精品| 日韩精品免费视频一区二区三区| 午夜影院在线不卡| 亚洲人成网站在线观看播放| 午夜福利,免费看| av有码第一页| 色网站视频免费| 大型av网站在线播放| 婷婷色麻豆天堂久久| 久久九九热精品免费| 久久精品久久久久久噜噜老黄| 两人在一起打扑克的视频| 91麻豆av在线| 久久午夜综合久久蜜桃| 国产av国产精品国产| 亚洲欧美一区二区三区久久| 又黄又粗又硬又大视频| 午夜激情av网站| 免费看av在线观看网站| 男人舔女人的私密视频| 波多野结衣av一区二区av| 一区二区日韩欧美中文字幕| 国产成人一区二区在线| 色94色欧美一区二区| 免费观看av网站的网址| 亚洲黑人精品在线| 中文欧美无线码| 极品人妻少妇av视频| 美女扒开内裤让男人捅视频| 国产人伦9x9x在线观看| 啦啦啦视频在线资源免费观看| 美国免费a级毛片| 下体分泌物呈黄色| 国产无遮挡羞羞视频在线观看| 丝瓜视频免费看黄片| 欧美精品高潮呻吟av久久| 国产欧美日韩一区二区三 | 欧美成人精品欧美一级黄| 男女午夜视频在线观看| 欧美97在线视频| 国产成人一区二区在线| 亚洲人成77777在线视频| 纵有疾风起免费观看全集完整版| 91国产中文字幕| 色婷婷av一区二区三区视频| 婷婷丁香在线五月| 精品第一国产精品| 久久精品国产综合久久久| 黄色视频在线播放观看不卡| 免费在线观看影片大全网站 | 久久女婷五月综合色啪小说| 国产成人精品无人区| 亚洲精品自拍成人| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品免费视频内射| 久久久久久免费高清国产稀缺| 老熟女久久久| 午夜福利一区二区在线看| www.999成人在线观看| 欧美日韩亚洲综合一区二区三区_| 精品福利永久在线观看| 这个男人来自地球电影免费观看| av福利片在线| 一边亲一边摸免费视频| 脱女人内裤的视频| 久久精品亚洲熟妇少妇任你| 性色av一级| 丝袜喷水一区| 深夜精品福利| 人体艺术视频欧美日本| 这个男人来自地球电影免费观看| bbb黄色大片| 91麻豆av在线| 国产成人免费无遮挡视频| 少妇粗大呻吟视频| 少妇被粗大的猛进出69影院| 亚洲视频免费观看视频| 在线精品无人区一区二区三| 天堂中文最新版在线下载| 日本av免费视频播放| 可以免费在线观看a视频的电影网站| 欧美精品av麻豆av| 人妻人人澡人人爽人人| 一区二区三区激情视频| 午夜免费成人在线视频| 尾随美女入室| 久久人妻熟女aⅴ| 亚洲天堂av无毛| 亚洲精品国产一区二区精华液| 亚洲,一卡二卡三卡| 色婷婷久久久亚洲欧美| 美女视频免费永久观看网站| 美女中出高潮动态图| 女性生殖器流出的白浆| 麻豆av在线久日| 色婷婷av一区二区三区视频| 亚洲国产精品一区三区| 2021少妇久久久久久久久久久| 黄色 视频免费看| avwww免费| 精品久久久久久久毛片微露脸 | 精品视频人人做人人爽| 亚洲中文av在线| 又大又爽又粗| 狂野欧美激情性xxxx| 亚洲,一卡二卡三卡| 一二三四在线观看免费中文在| 伊人久久大香线蕉亚洲五| 考比视频在线观看| 久久中文字幕一级| 免费观看a级毛片全部| 午夜久久久在线观看| 亚洲国产最新在线播放| 亚洲av综合色区一区| 日韩电影二区| 亚洲人成77777在线视频| 91国产中文字幕| 亚洲精品国产一区二区精华液| 大片免费播放器 马上看| 久久久精品区二区三区| 日韩 欧美 亚洲 中文字幕| 欧美激情极品国产一区二区三区| 日韩人妻精品一区2区三区| 桃花免费在线播放| 日日爽夜夜爽网站| 国产精品一二三区在线看| 日韩 欧美 亚洲 中文字幕| 波野结衣二区三区在线| 久久九九热精品免费| 亚洲国产中文字幕在线视频| 天天躁日日躁夜夜躁夜夜| 国产成人欧美| 91麻豆精品激情在线观看国产 | 曰老女人黄片| 亚洲成人免费av在线播放| www.av在线官网国产| 十八禁人妻一区二区| 一级片'在线观看视频| 色播在线永久视频| 每晚都被弄得嗷嗷叫到高潮| 18禁国产床啪视频网站| 欧美精品一区二区大全| 在线av久久热| 精品少妇内射三级| 飞空精品影院首页| 尾随美女入室| 1024视频免费在线观看| 又粗又硬又长又爽又黄的视频| 在线观看免费高清a一片| 亚洲欧美激情在线| 久久久精品免费免费高清| 国产免费福利视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产精品一区二区免费欧美 | 亚洲午夜精品一区,二区,三区| 91成人精品电影| 亚洲人成77777在线视频| 亚洲国产av影院在线观看| av不卡在线播放| 少妇被粗大的猛进出69影院| 操出白浆在线播放| 亚洲精品日本国产第一区| 欧美黑人欧美精品刺激| 成人手机av| 成人影院久久| 亚洲专区国产一区二区| av欧美777| 亚洲av日韩在线播放| 少妇人妻久久综合中文| 亚洲九九香蕉| 别揉我奶头~嗯~啊~动态视频 | 伦理电影免费视频| 老司机在亚洲福利影院| 美女中出高潮动态图| 51午夜福利影视在线观看| 亚洲伊人久久精品综合| 国产成人a∨麻豆精品| 999精品在线视频| 高清视频免费观看一区二区| 欧美亚洲 丝袜 人妻 在线| 国产高清videossex| 久久久久国产精品人妻一区二区| 一区福利在线观看| 永久免费av网站大全| 777米奇影视久久| 国产精品熟女久久久久浪| 在线观看人妻少妇| 人体艺术视频欧美日本| 操美女的视频在线观看| 99九九在线精品视频| 极品少妇高潮喷水抽搐| 久久久久久久久久久久大奶| 亚洲精品av麻豆狂野| 亚洲av在线观看美女高潮| 亚洲男人天堂网一区| 在线观看免费日韩欧美大片| 青春草亚洲视频在线观看| 99热全是精品| 91麻豆精品激情在线观看国产 | 超碰97精品在线观看| 亚洲精品一二三| 亚洲人成电影观看| 精品第一国产精品| 精品国产超薄肉色丝袜足j| 亚洲av美国av| 波多野结衣一区麻豆| 亚洲色图综合在线观看| 亚洲精品国产av蜜桃| 日韩熟女老妇一区二区性免费视频| 日本午夜av视频| 国产成人免费无遮挡视频| 自线自在国产av| 黑丝袜美女国产一区| 久久久久久亚洲精品国产蜜桃av| 婷婷色综合www| 国产日韩欧美在线精品| 纵有疾风起免费观看全集完整版| 国产男女内射视频| 久热爱精品视频在线9| 亚洲成人国产一区在线观看 | 一二三四在线观看免费中文在| 欧美黑人精品巨大| 亚洲精品日韩在线中文字幕| 两个人看的免费小视频| 丝袜在线中文字幕| 亚洲av日韩精品久久久久久密 | av在线播放精品| 成人国产一区最新在线观看 | 啦啦啦在线免费观看视频4| 久久人人97超碰香蕉20202| 成人亚洲欧美一区二区av| 一本色道久久久久久精品综合| 在现免费观看毛片| 两个人免费观看高清视频| 亚洲人成电影观看| 精品少妇久久久久久888优播| 国产爽快片一区二区三区| 五月开心婷婷网| 两人在一起打扑克的视频| 免费日韩欧美在线观看| 午夜免费观看性视频| 丝袜人妻中文字幕| 看免费av毛片| 伊人久久大香线蕉亚洲五| 精品国产一区二区三区久久久樱花| 久久久久视频综合| 国产精品99久久99久久久不卡| 亚洲成人国产一区在线观看 | 麻豆乱淫一区二区| 日本黄色日本黄色录像| 日韩制服骚丝袜av| 国产精品熟女久久久久浪| 国产精品久久久人人做人人爽| 青草久久国产| 午夜福利视频精品| 亚洲欧美色中文字幕在线| 久久久久久久久久久久大奶| 欧美久久黑人一区二区| 亚洲专区中文字幕在线| 少妇裸体淫交视频免费看高清 | 视频区图区小说| 亚洲欧美清纯卡通| 国产av一区二区精品久久| 欧美日韩视频精品一区| 国产免费一区二区三区四区乱码| 国产精品久久久av美女十八| 久久天躁狠狠躁夜夜2o2o | 80岁老熟妇乱子伦牲交| 超色免费av| 亚洲,一卡二卡三卡| 一边摸一边做爽爽视频免费| 国产伦理片在线播放av一区| 别揉我奶头~嗯~啊~动态视频 | 欧美乱码精品一区二区三区| 青春草视频在线免费观看| 青草久久国产| 丰满少妇做爰视频| 欧美变态另类bdsm刘玥| av福利片在线| www日本在线高清视频| 人妻一区二区av| 色婷婷av一区二区三区视频| 九色亚洲精品在线播放| 成年人免费黄色播放视频| 女人精品久久久久毛片| 97精品久久久久久久久久精品| 成年人免费黄色播放视频| 久久久久国产一级毛片高清牌| 99热全是精品| 老司机影院毛片| 国产片内射在线| 丝瓜视频免费看黄片| 老鸭窝网址在线观看| 啦啦啦在线免费观看视频4| 国产又爽黄色视频| 久久中文字幕一级| 免费久久久久久久精品成人欧美视频| 在线观看www视频免费| 美女脱内裤让男人舔精品视频| 国产成人啪精品午夜网站| 欧美成人精品欧美一级黄| 日本a在线网址| 我的亚洲天堂| 欧美精品高潮呻吟av久久| 青草久久国产| 欧美 亚洲 国产 日韩一| 国产成人一区二区在线| 久久久久久亚洲精品国产蜜桃av| 国产亚洲av高清不卡| 黑丝袜美女国产一区| 久久ye,这里只有精品| 老熟女久久久| 亚洲精品日韩在线中文字幕| 亚洲av片天天在线观看| 国产精品国产av在线观看| 精品久久久久久电影网| 侵犯人妻中文字幕一二三四区| 激情视频va一区二区三区| 香蕉丝袜av| 亚洲情色 制服丝袜| videos熟女内射| 久久99精品国语久久久| 婷婷色综合www| 丁香六月天网| 成在线人永久免费视频| 国产一区二区在线观看av| 91精品国产国语对白视频| 欧美 日韩 精品 国产| 在线精品无人区一区二区三| 久久性视频一级片| 久久影院123| 性色av一级| 中文字幕人妻丝袜一区二区| 大香蕉久久成人网| 久久热在线av| 久久精品久久久久久久性| 欧美激情 高清一区二区三区| a级毛片在线看网站| 黑人猛操日本美女一级片| 欧美精品亚洲一区二区| 一区二区日韩欧美中文字幕| 国产精品偷伦视频观看了| 美女福利国产在线| 国产精品香港三级国产av潘金莲 | av天堂久久9| 免费在线观看影片大全网站 | 欧美在线一区亚洲| 免费久久久久久久精品成人欧美视频| av不卡在线播放| 一本一本久久a久久精品综合妖精| 国产精品成人在线| 精品一区二区三区av网在线观看 | 亚洲国产精品成人久久小说| 午夜福利一区二区在线看| 国产成人91sexporn| 亚洲av综合色区一区| 亚洲美女黄色视频免费看| 欧美性长视频在线观看| www.av在线官网国产| 久久久精品94久久精品| 亚洲伊人色综图| 如日韩欧美国产精品一区二区三区| 免费日韩欧美在线观看| videosex国产| 人人澡人人妻人| 久久亚洲精品不卡| 99久久精品国产亚洲精品| 久久久久久免费高清国产稀缺| 纵有疾风起免费观看全集完整版| 午夜日韩欧美国产| 中文字幕人妻丝袜一区二区| 国产男女内射视频| 国产亚洲精品第一综合不卡| 日本午夜av视频| 久久人人爽人人片av| 亚洲欧美一区二区三区黑人| 精品高清国产在线一区| 深夜精品福利| 欧美少妇被猛烈插入视频| av在线app专区| 99香蕉大伊视频| 后天国语完整版免费观看| 国产精品免费视频内射| 国产亚洲精品久久久久5区| 久久国产亚洲av麻豆专区| 久久av网站| 免费av中文字幕在线| 欧美精品高潮呻吟av久久| 国产精品秋霞免费鲁丝片| 超碰97精品在线观看| 久久精品成人免费网站| 国产老妇伦熟女老妇高清| 色网站视频免费| 中文字幕av电影在线播放| 亚洲男人天堂网一区| 免费av中文字幕在线| 日本欧美视频一区| av片东京热男人的天堂| 三上悠亚av全集在线观看| 成人手机av| 麻豆乱淫一区二区| 欧美精品人与动牲交sv欧美| 国产精品久久久久成人av| 免费高清在线观看日韩| 高清不卡的av网站| 99re6热这里在线精品视频| 伦理电影免费视频| 午夜福利乱码中文字幕| 大码成人一级视频| 久久久国产欧美日韩av| 精品久久久久久电影网| 天天躁夜夜躁狠狠躁躁| 亚洲熟女精品中文字幕| 建设人人有责人人尽责人人享有的| 欧美精品人与动牲交sv欧美| 日本欧美国产在线视频| 人体艺术视频欧美日本| 大片免费播放器 马上看| 国产精品九九99| av国产精品久久久久影院| 一本一本久久a久久精品综合妖精| 久久久精品94久久精品| av欧美777| 国产男女内射视频| 青春草视频在线免费观看| av又黄又爽大尺度在线免费看| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 成人国产一区最新在线观看 | 麻豆乱淫一区二区| 亚洲av国产av综合av卡| 国产精品一区二区精品视频观看| 永久免费av网站大全| www日本在线高清视频| 日本a在线网址| 国产精品熟女久久久久浪| 成年美女黄网站色视频大全免费| av电影中文网址| 国产有黄有色有爽视频| 天天躁日日躁夜夜躁夜夜| 老司机影院毛片| 亚洲国产最新在线播放| 成在线人永久免费视频| 欧美乱码精品一区二区三区| 最近中文字幕2019免费版| 欧美日韩国产mv在线观看视频| 国产男女超爽视频在线观看| 亚洲av片天天在线观看| 老司机靠b影院| 国产一区二区三区av在线| 亚洲av综合色区一区| 成人黄色视频免费在线看| 免费在线观看影片大全网站 | 制服人妻中文乱码| 黄片小视频在线播放| 久久综合国产亚洲精品| 高清视频免费观看一区二区| 黄片播放在线免费| 国产在线视频一区二区| 一区二区三区四区激情视频| 国产精品九九99| 91精品伊人久久大香线蕉| 成人亚洲精品一区在线观看| 精品卡一卡二卡四卡免费| 国产av精品麻豆| 超碰成人久久| 交换朋友夫妻互换小说| 狂野欧美激情性bbbbbb| 最新的欧美精品一区二区| 性色av一级| 免费观看a级毛片全部| 国产爽快片一区二区三区|