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

    未來RCPs情景下湖北省中稻高溫熱害風險研究

    2022-04-02 02:28:22劉可群劉凱文秦鵬程趙小芳楊青青
    湖北農(nóng)業(yè)科學 2022年5期
    關鍵詞:中稻增幅湖北省

    葉 佩 ,劉可群 ,劉凱文 ,3,秦鵬程 ,趙小芳 ,楊青青

    (1.荊州農(nóng)業(yè)氣象試驗站,湖北 荊州 434000;2.武漢區(qū)域氣候中心,武漢 430074;3.長江大學農(nóng)學院,湖北 荊州 434025)

    近年來,全球氣候變暖對農(nóng)業(yè)生產(chǎn)和發(fā)展的不利影響日趨顯著,糧食安全生產(chǎn)問題面臨更加嚴峻的挑戰(zhàn)[1]。水稻(Oryza sativaL.)是中國重要的三大主要糧食作物之一,其種植面積占世界種植總面積的26.9%,產(chǎn)量占世界稻谷總產(chǎn)量的33.9%[2]。在中國有65%以上的人口以水稻為主食,水稻的安全生產(chǎn)對保障中國糧食安全、促進農(nóng)戶增收以及推動農(nóng)村經(jīng)濟發(fā)展有著舉足輕重的作用[3]。湖北省作為華中地區(qū)主要水稻種植區(qū)之一,水稻種植面積超過該省糧食作物種植總面積的50%以上,產(chǎn)量占全省糧食總產(chǎn)的70%左右,為該省的第一大主糧作物,其中,水稻種植類型以中稻為主[4]。夏季為水稻生長的關鍵時期,此時受副熱帶高壓影響,水稻生育期內(nèi)出現(xiàn)高溫熱害的可能性大[5]。隨著全球氣候變暖進程的加劇,湖北省中稻高溫熱害發(fā)生的強度和頻率呈增強、增多的趨勢[6],在抽穗揚花期水稻遭遇極端高溫容易造成花粉發(fā)育不良,使水稻空秕率增加,結實率下降,最終導致減產(chǎn),影響水稻的產(chǎn)量和效益[7,8]。劉文茹等[9]運用 BCC-CSM1-1 氣候系統(tǒng)模式輸出的典型濃度情景RCPs 數(shù)據(jù),發(fā)現(xiàn)RCP4.5 情景下長江中下游地區(qū)夏季氣溫變化速率較大,增幅每10 年達1.2 ℃。因此,研究未來氣候變化情景下高溫熱害對湖北省中稻的影響及其風險變化十分有必要。

    高溫熱害作為影響水稻安全生產(chǎn)的重要環(huán)境生態(tài)因子,未來氣候變化背景下該因子對水稻的影響及風險研究日益受到學者的關注。李琪等[10]評估了未來氣候變化情景下孕穗-抽穗期不同強度和持續(xù)時間高溫熱害對水稻產(chǎn)量的影響,王連喜等[11]分析了未來氣候變化情景下孕穗-抽穗期水稻對高溫的適應性。宋瑞明等[12]基于CanESM2 模式SDSM統(tǒng)計降尺度后的RCPs 情景數(shù)據(jù)預估未來2021—2070 年氣候變化條件下江蘇水稻高溫熱害變化規(guī)律。萬素琴等[13]將區(qū)域氣候模式PRECIS 構建的氣候變化情景與水稻生產(chǎn)模型ORYZA2000 相結合,研究了未來氣候變化對華中地區(qū)中稻產(chǎn)量的影響。熊偉等[14]預估 2021—2050 年 RCP2.6 和 RCP8.5 氣候情景下,全國水稻關鍵生育期高溫日數(shù)、高溫積溫均呈增加趨勢,水稻高溫熱害風險也將進一步加大。目前大部分研究的氣候情景數(shù)據(jù)鮮有訂正,加上未來氣候情景下的湖北省中稻高溫熱害風險也未見報道,因此,本研究采用方差訂正后的RCPs 數(shù)據(jù)來探究未來氣候情景下的湖北省中稻高溫熱害風險時空分布特征,以期為該區(qū)域高溫熱害風險預警、水稻生產(chǎn)布局、防災減災等管理決策提供理論依據(jù)。

    1 資料來源與研究方法

    1.1 資料來源

    農(nóng)業(yè)氣象資料選用湖北省18 個農(nóng)業(yè)氣象觀測站(圖1)中稻抽穗開花期觀測數(shù)據(jù),時間為1981—2016 年;氣象資料為湖北省 1960—2019 年 76 個氣象站點逐日氣象數(shù)據(jù),包括最高氣溫、平均氣溫,數(shù)據(jù)來源于湖北省氣象局氣候檔案室。氣候模式數(shù)據(jù)來自 ISI-MIP(Inter-Sectoral Impact Model Intercomparison Project)推薦使用的5 個全球海氣耦合模式中對最高氣溫模擬性能最好的HadGem2-ES(Had-Gem)模式下的 1961—2080 年數(shù)據(jù),其中 1961—1990 年為基準時段,1991—2005 年為驗證時段,2021—2080 年為未來時段,該數(shù)據(jù)基于WATCH 再分析資料和概率分布統(tǒng)計偏差訂正方法進行了降尺度處理[15],空間分辨率為 0.5°×0.5°,研究路徑為穩(wěn)定路徑 RCP4.5 和高端路徑 RCP8.5[9]。

    圖1 湖北省18 個農(nóng)業(yè)氣象觀測站點分布

    1.2 研究方法

    1.2.1 研究時段 由表1 可知,7 月下旬至8 月下旬為湖北省中稻抽穗開花期,是中稻最易受高溫危害的時段,因此,將7 月21 至8 月31 日作為中稻高溫熱害風險研究時段。

    表1 湖北省18 個農(nóng)業(yè)氣象站點中稻抽穗期

    1.2.2 高溫熱害指標 高溫危害積溫同時包含高溫強度和持續(xù)時間的特征,且能定量描述高溫熱害影響程度,應用廣泛[16],因此,本研究以水稻抽穗開花期超過35 ℃以上的高溫危害積溫評價高溫發(fā)生強度,根據(jù)田間試驗[16]與實際調(diào)查[7],高溫危害積溫達到 16、32、48 ℃·d 分別為輕度、中度、重度水稻熱害指標。高溫熱害發(fā)生頻率為歷年出現(xiàn)高溫的總年數(shù)與參與統(tǒng)計的總年數(shù)之比。

    1)高溫危害積溫。

    2)高溫發(fā)生頻率。

    式中,HDD 為高溫危害積溫,DTi是水稻生育期內(nèi)逐日高溫危害積溫,Tmax為逐日最高氣溫,T0為35,n為研究時段水稻生育期總天數(shù);i為水稻生育期內(nèi)第i天,i=1,2,…,n。P為高溫熱害發(fā)生頻率;M為該站發(fā)生某等級高溫熱害年數(shù);N為全部總年數(shù)。

    1.2.3 RCPs 數(shù)據(jù)訂正 RCPs 情景下的全球氣候模式(GCM)數(shù)據(jù)的模擬效果較好[17],但與實測資料還存在一定的誤差,降尺度和方差訂正可在一定程度上彌補氣候模式數(shù)據(jù)分辨率精度不夠的缺陷,提高氣候模式數(shù)據(jù)的準確度[18]。

    首先運用雙線性內(nèi)插法將格點數(shù)據(jù)插值到觀測站點F(X,Y)上,假設模擬格點數(shù)據(jù)的經(jīng)緯度坐標為F11(X1,Y1)、F12(X1,Y2)、F21(X2,Y1)、F22(X2,Y2),4個格點模擬的最高氣溫分別為T11、T12、T21、T22,首先在X方向進行線性插值,然后在Y方向進行線性插值,得到雙線性插值結果。

    其次利用1961—1990 年校準時段模式數(shù)據(jù)與實測氣象站點數(shù)據(jù)通過傳遞函數(shù)系數(shù)和方差信息建立方差訂正模型,具體如下。

    第一步,將站點校準年1961—1990 年觀測值、歷史模擬數(shù)據(jù) 1961—1990 年、未來 2021—2050 和2051—2080 年需訂正模擬數(shù)據(jù),轉化為同期等長的序列其中i=1,…,n,根據(jù)公式(5)至公式(8)計算旬均值,并以基準時期模擬數(shù)據(jù)旬均值序列(Xs)作為自變量,同一旬的觀測數(shù)據(jù)旬均值作為因變量(Xo)建立傳遞函數(shù)。

    式中,k為某年某旬天數(shù),取實際值8 或9 或10或11,用最小二乘法計算系數(shù)a、b、c的值。

    1.2.4 高溫熱害風險評估模型 采用各站點不同等級高溫熱害發(fā)生強度和頻率建立高溫熱害風險評估模型。

    式中,G 表示高溫熱害風險;n表示不同等級高溫熱害,n=1,2,3 分別代表輕度、中度及重度高溫熱害;Wi表示不同等級高溫熱害權重;輕度的權重為0.1,中度為0.2,重度為0.3;Pi表示研究時段不同等級高溫熱害發(fā)生頻率。

    1.2.5 模型訂正 以武漢站1 月上旬逐日最高氣溫訂正為例,步驟如下。①利用基準時段1961—1990年模擬數(shù)據(jù)和觀測資料,通過公式(8)、公式(9)分別計算逐年1 月上旬旬均值,得到容量為30 的模擬和觀測序列;②以觀測序列為因變量,模擬序列為自變量代入式(10)進行回歸分析,得到1 月上旬最高氣溫回歸系數(shù)a、b、c(表2);③將回歸系數(shù)帶入式(11),建立傳遞函數(shù),對該旬逐日模擬結果進行均值訂正;④計算旬方差訂正因子λ并代入式(12)進行方差訂正。該旬逐日最高氣溫訂正模型如下。

    表2 武漢站逐旬最高氣溫傳遞函數(shù)系數(shù)

    1.2.6 模型驗證 利用1991—2005 年氣候模擬數(shù)據(jù)插值到氣象站點形成待訂正數(shù)據(jù),根據(jù)各站點最高氣溫的回歸系數(shù)和方差訂正因子,建立訂正模型,得到1991—2005 年訂正后逐日模擬數(shù)據(jù)。通過觀測數(shù)據(jù)和訂正前后模擬數(shù)據(jù)之間的逐月相對誤差(RE)和均方根誤差(RMSE)分析最高氣溫的訂正效果。其中,RE用于表示模式值偏離觀測值的實際大小,數(shù)值越小,偏離觀測值的程度越?。籖MSE反映模擬值相對于觀測值的精確程度,數(shù)值越小,精確度越高。

    式中,simi為訂正后逐月模擬值,obsi為逐月觀測值,n為樣本容量。

    分別計算訂正后各站逐月最高氣溫與觀測逐月最高氣溫的RE(圖2a)和RMSE(圖2b)。方差訂正后數(shù)據(jù)與觀測值的RE和RMSE較訂正前明顯減小,且 訂 正 后 7—8 月 最 高 氣 溫RE為 0.04、RMSE在1.5 ℃范圍內(nèi),使用方差訂正對未來情景下的模式數(shù)據(jù)進行訂正可有效降低模擬數(shù)據(jù)與觀測值的偏差。

    圖2 1991—2005 年訂正前后最高氣溫與觀測值的相對誤差(a)和均方根誤差(b)

    2 結果與分析

    2.1 RCPs 情景高溫熱害強度變化特征

    2.1.1 2021—2050 年高溫熱害強度變化特征 由圖 3 可知,2021—2050 年不同 RCPs 情景下中稻高溫熱害強度增減不一,尤以8.5 情景變幅更明顯,其中鄂西南和鄂東南增幅最大,鄂東北和江漢平原次之,在鄂西北出現(xiàn)遞減趨勢。RCR4.5 情景下增幅高值區(qū)集中在巴東、興山、秭歸等地,增幅分別為19.07、23.22、18.75 ℃·d;其次為鄂東南的赤壁、咸寧、通城、陽新等站點,其增幅分別為12.39、11.65、10.07、10.25 ℃·d;鄂西北和鄂東北的大部地區(qū)增幅最小,其增幅維持在10 ℃·d 以下(圖3a)。RCP8.5 情景下高溫危害積溫強度增幅較大的區(qū)域主要集中在鄂西南的巴東-興山-秭歸及鄂東南的咸寧、赤壁一帶,其中巴東、秭歸、興山、恩施等站點增幅超過20 ℃·d;增幅較小的區(qū)域主要位于江漢平原的鐘祥、天門以及鄂東北部的隨州、麻城一帶,增幅均不超過5 ℃·d;增幅最小的區(qū)域主要位于鄂西北部的鄖西、鄖縣、竹山、老河口、谷城一帶,高溫熱害強度有下降趨勢(圖3b)。

    圖3 RCPs情景下2021—2050 年高溫危害積溫較基準期的空間變化特征

    2.1.2 2051—2080 年高溫積溫強度變化特征 由圖 4 可知,2051—2080 年不同 RCPs 情景下中稻高溫熱害強度增量呈現(xiàn)南高北低的空間分布,RCP4.5 情景下,增幅高值區(qū)為巴東-恩施-興山等站點,增幅均超過了20 ℃·d,其次為赤壁-崇陽-咸寧-通山一帶,增幅均超過了15 ℃·d,鄂東北和江漢平原南部高溫危害積溫增幅不超過10 ℃·d,鄂西北和鄂東北北部高溫危害積溫增幅不超過5 ℃·d,在鄖西、鄖縣、竹山等站點高溫危害積溫有降低趨勢。RCP8.5情景下高溫危害積溫變幅范圍更大,空間分布變幅顯著,其中在鄂西北部高溫危害積溫有顯著降低趨勢,其余各地區(qū)為顯著增加趨勢。整體上湖北省各區(qū)域高溫熱害強度呈增加的趨勢,RCP8.5 情景高溫熱害強度較RCP4.5 增減幅度大,RCP8.5 情景較RCP4.5 波動劇烈。與 2021—2050 年相比,2051—2080 年高溫危害積溫增幅更大。

    圖4 RCPs情景下2051—2080 年高溫危害積溫較基準期的空間變化特征

    2.2 RCPs 情景下高溫熱害頻率變化特征

    2.2.1 2021—2050年高溫熱害頻率變化特征 RCP4.5、RCP8.5 情景下2021—2050 年湖北省各區(qū)域不同等級高溫熱害發(fā)生頻率增減不一,RCP8.5 情景增減變化更為顯著(圖5)。鄂西南為重度高溫熱害頻率增幅最大的地區(qū),RCP4.5、RCP8.5 情景下重度高溫熱害頻率增幅分別達7.84%、12.35%。鄂東南主要表現(xiàn)為中度等級的高溫熱害頻率增加,RCP4.5、RCP8.5 情景下中度高溫熱害頻率增幅分別達到12.98%、14.79%。鄂東北和江漢平原主要表現(xiàn)為輕度等級的高溫熱害頻率大幅增加,RCP4.5、RCP8.5情景下輕度增幅分別為14.29%、23.81% 以及16.90%、22.62%。鄂西北高溫熱害發(fā)生頻率有降低趨勢,RCP4.5 情景下輕度、中度、重度、全部高溫熱害頻率較基準期的變化分別為-1.11%、0.67%、0、-0.44%,RCP8.5 情 景 下 分 別 為 -0.44%、0、0、-0.44%,屬全省高溫熱害頻率低變化區(qū)。2 種RCP情景下湖北省全區(qū)2021—2050 年各等級高溫熱害發(fā)生頻率相對基準年呈現(xiàn)增加趨勢,主要以輕度增幅為主,RCP4.5 和RCP8.5 情景下全省各等級高溫熱害發(fā)生頻率變幅分別為7.26%、5.28%、3.27%、15.81%以及10.94%、5.04%、5.20%、21.18%。

    圖5 RCPs情景下2021—2050 年湖北省各區(qū)域不同等級高溫熱害頻率較基準期變化特征

    2.2.2 2051—2080 年高溫熱害頻率變化特征 由圖 6 可見,與 2021—2050 年相比,2051—2080 年高溫熱害發(fā)生頻率和程度變幅更大,尤以RCP8.5 顯著。其中鄂西南在RCP4.5 情景下以輕度高溫熱害頻率增幅為主,為10.78%,在RCP8.5 情景下以重度高溫熱害頻率增幅為主,達到27.06%。鄂東南在RCP4.5 情景下以中度高溫熱害頻率增加最顯著,為11.51%,在RCP8.5 情景下以重度高溫熱害頻率增幅為主,為19.47%。鄂東北和江漢平原為以輕度高溫熱害頻率增幅最為顯著,RCP4.5 和RCP8.5 情景下增幅為16.43%、29.52%以及16.67%、21.67%。鄂西北地區(qū)高溫熱害發(fā)生頻率有降低的趨勢,在RCP8.5情景下輕度高溫熱害頻率減幅為5.33%。對湖北省全區(qū)而言,2051—2080 年中稻高溫熱害發(fā)生頻率有顯著增加的趨勢,RCP4.5 和RCP8.5 情景下不同程度的高溫熱害增幅分別為9.10%、5.84%、4.96%、19.91% 以及 9.58%、6.78%、11.64%、28.01%,其中RCP4.5 情景下主要以輕度為主,RCP8.5 情景下以重度為主,表明RCP8.5 情景下極端高溫熱害發(fā)生更頻繁。

    圖6 RCPs情景下2051—2080 年湖北省各區(qū)域不同等級高溫熱害頻率較基準期變化特征

    2.3 RCPs 情景下高溫熱害風險指數(shù)

    2.3.1 RCPs 情景下2021—2050 年高溫熱害風險指數(shù) 由圖7 可見,2021—2050 年高溫熱害風險呈現(xiàn)由北向南遞增,尤其以RCP8.5 顯著,鄂西南低山、鄂東南等地高溫熱害風險增幅較高;其次為鄂東北和江漢平原大部地區(qū),在鄂西北地區(qū)高溫熱害風險為降低趨勢。這種高溫熱害風險變化趨勢與近年來湖北省各地區(qū)出現(xiàn)極端高溫的強度、頻率變化及對區(qū)域性增溫的貢獻有關[19]。RCP4.5 情景下高溫熱害風險增幅較高的地區(qū)主要為鄂東南的咸寧、赤壁以及鄂西南的巴東等地,其風險超過0.12;風險變幅較小的地區(qū)主要位于鄂西北和江漢平原大部,其中鄖西、鄖縣、竹山等站點風險為0。RCP8.5 情景下高溫熱害風險增幅較高的地區(qū)分布在鄂西南的巴東、長陽以及鄂東南的咸寧、赤壁、通城、崇陽等地;風險增幅較低的地區(qū)則主要分布在鄂西北一帶,其中在鄖西、鄖縣、老河口、竹山等站風險變化為負值,該地區(qū)風險有降低趨勢。未來30 年湖北省中稻高溫熱害風險呈增加的趨勢,RCP4.5、RCP8.5 情景下的風險平均增幅分別為0.05 和0.06,RCP8.5 情景下的風險變化相對于RCP4.5 情景更劇烈。

    圖7 RCPs情景下2021—2050 年湖北省中稻高溫熱害風險指數(shù)

    2.3.2 RCPs 情景下2051—2080 年高溫熱害風險指數(shù) 由圖8 可知,2051—2080 年的高溫熱害風險變化較2021—2050 年更加顯著,在RCP8.5 情景下高溫熱害風險增幅呈明顯帶狀的南高北低趨勢。其中鄂西南和鄂東南風險增幅更大,鄂西北減幅也更明顯。鄂東南和鄂西南高溫熱害強度及發(fā)生頻率增加顯著,導致高溫熱害危險性大幅增加,鄂西北地區(qū)高溫熱害強度及頻率均有降低趨勢,該地區(qū)高溫熱害危險性有所緩解。

    圖8 RCPs情景下2051—2080 年湖北省中稻高溫熱害風險指數(shù)

    3 小結與討論

    3.1 小結

    1)RCPs 情景數(shù)據(jù)經(jīng)過方差訂正,能有效降低模擬數(shù)據(jù)與觀測值的相對誤差和均方根誤差,訂正后的RCPs 數(shù)據(jù)對未來氣候變化趨勢模擬更具有代表性。

    2)未來 RCPs 情景下,2051—2080 年湖北省中稻高溫熱害強度、頻率及風險較基準時段變幅大于2021—2050 年,尤以 RCP8.5 情景顯著。

    3)2 種RCPs 情景下,中稻高溫熱害頻率、強度增加,導致湖北省中稻高溫熱害風險整體呈增加的趨勢。未來湖北省中稻高溫熱害風險增幅最大的地區(qū)位于鄂東南大部和鄂西南低山地區(qū),該地區(qū)高溫熱害主要以中度和重度為主,其次為江漢平原和鄂東北地區(qū),該地區(qū)高溫熱害主要以輕度為主,鄂西北地區(qū)有降低的趨勢,其高溫熱害頻率、強度有減小的趨勢。RCP8.5 情景下的風險變幅高于RCP4.5 情景,RCP8.5 情景下的不穩(wěn)定性較高。

    3.2 討論

    本研究通過訂正后的RCP 數(shù)據(jù)來研究未來氣候變化情景下中稻高溫熱害風險變化,可提高預測風險的精確度,為農(nóng)業(yè)生產(chǎn)防災減災提供理論依據(jù)。但未考慮未來水稻品種更新、田間管理制度以及氣候變化對生育期的改變,使得未來情景下高溫熱害的統(tǒng)計結果存在一定不確定性[20];每種氣候模式都有其擅長和不擅長模擬的區(qū)域及氣候要素,因此,本研究選擇最高氣溫模擬性能最好的(HadGem)模式下模擬的HadGem2-ES 數(shù)據(jù),但研究認為多模式集合的概率分布是降低氣候變化影響評估的不確定性的有效方法,接下來則更需要研究多模式結合作物生長模型模擬中稻生產(chǎn)面臨的高溫熱害風險趨勢[21]。未來 RCPs 情景下,2051—2080 年湖北省中稻高溫熱害強度、頻率和風險均明顯高于2021—2050 年,說明隨著未來氣候變化加劇,將進一步危害該地區(qū)中稻安全生產(chǎn)[22]。在空間上,鄂西南低山和鄂東南高溫熱害強度、頻率和風險較基準時段的增幅最大,該地區(qū)不利于中稻的安全種植,可通過調(diào)整其生產(chǎn)布局[23]和播期[24]以及篩選培育抗逆品種來緩解高溫熱害;鄂東北和江漢平原次之,在鄂西北出現(xiàn)遞減趨勢,未來該地區(qū)在水資源充足的情況下可擴大中稻的種植面積。

    猜你喜歡
    中稻增幅湖北省
    2023年全省一般公共預算收支情況一覽表
    山西財稅(2024年1期)2024-03-15 06:29:00
    上半年全省一般公共預算收支情況一覽表
    山西財稅(2022年7期)2022-08-18 01:27:48
    2016年湘審中稻品種(上)
    2016年汽車產(chǎn)銷雙雙超過2800萬輛增幅比上年較快提升
    汽車與安全(2017年2期)2017-03-14 02:56:27
    內(nèi)蒙古自治區(qū)2015年審定通過的水稻新品種
    中國稻米(2016年1期)2016-06-29 07:59:53
    江西省2014年審定通過的水稻新品種(下)
    中國稻米(2016年1期)2016-06-29 07:59:38
    2015年湘審中稻品種(下
    春風化雨 深入淺出——記湖北省優(yōu)秀理論宣講員龍會忠
    學習月刊(2015年19期)2015-07-09 03:38:34
    湖北省2015年3月生產(chǎn)資料價格
    儀器儀表業(yè)今年產(chǎn)銷增幅或達15%
    自動化博覽(2014年4期)2014-02-28 22:31:11
    亚洲成国产人片在线观看| 91午夜精品亚洲一区二区三区| 在线观看www视频免费| 啦啦啦视频在线资源免费观看| 性少妇av在线| 电影成人av| 伦理电影大哥的女人| 亚洲人成网站在线观看播放| 亚洲美女搞黄在线观看| 曰老女人黄片| 亚洲成人手机| 亚洲欧洲日产国产| 26uuu在线亚洲综合色| 男女下面插进去视频免费观看| 国产毛片在线视频| 十八禁高潮呻吟视频| 午夜福利影视在线免费观看| 一级片'在线观看视频| 汤姆久久久久久久影院中文字幕| 久久久国产精品麻豆| 热99久久久久精品小说推荐| 天天躁夜夜躁狠狠躁躁| 亚洲欧美日韩另类电影网站| 王馨瑶露胸无遮挡在线观看| 一区二区三区四区激情视频| 欧美另类一区| 日韩中字成人| 成人毛片60女人毛片免费| 欧美97在线视频| 国产午夜精品一二区理论片| 如何舔出高潮| 一级黄片播放器| 观看美女的网站| 国产1区2区3区精品| av免费在线看不卡| 青草久久国产| 一级片免费观看大全| 亚洲成色77777| 国产一区二区激情短视频 | 中文欧美无线码| 久久久久久久久久人人人人人人| 亚洲人成网站在线观看播放| 男女午夜视频在线观看| 久久久久久免费高清国产稀缺| av电影中文网址| 王馨瑶露胸无遮挡在线观看| 国产精品欧美亚洲77777| 亚洲欧美成人综合另类久久久| 伊人亚洲综合成人网| 看十八女毛片水多多多| 亚洲第一av免费看| 亚洲精品国产av蜜桃| 最近的中文字幕免费完整| 性色avwww在线观看| 国产精品女同一区二区软件| 亚洲成国产人片在线观看| 亚洲欧洲国产日韩| 91国产中文字幕| 久久99一区二区三区| 亚洲第一青青草原| 国产不卡av网站在线观看| 波多野结衣av一区二区av| 巨乳人妻的诱惑在线观看| 国产在线一区二区三区精| 亚洲三区欧美一区| 97精品久久久久久久久久精品| 中文字幕精品免费在线观看视频| 激情五月婷婷亚洲| 777米奇影视久久| 桃花免费在线播放| 日韩一卡2卡3卡4卡2021年| 制服诱惑二区| 2022亚洲国产成人精品| 日本-黄色视频高清免费观看| 在线天堂最新版资源| 色网站视频免费| 亚洲中文av在线| 日韩,欧美,国产一区二区三区| 亚洲第一青青草原| 99国产精品免费福利视频| 1024香蕉在线观看| 欧美最新免费一区二区三区| 国产精品av久久久久免费| 亚洲国产日韩一区二区| 精品国产一区二区三区久久久樱花| 亚洲第一区二区三区不卡| 成人免费观看视频高清| 国产毛片在线视频| 久久久久久久亚洲中文字幕| 免费黄网站久久成人精品| 国产日韩欧美视频二区| 久久国内精品自在自线图片| 高清不卡的av网站| 丰满迷人的少妇在线观看| 熟女少妇亚洲综合色aaa.| 一区二区av电影网| 亚洲av中文av极速乱| 欧美激情高清一区二区三区 | 十分钟在线观看高清视频www| 午夜福利,免费看| 亚洲一区中文字幕在线| 日本vs欧美在线观看视频| 亚洲欧美成人精品一区二区| 热99国产精品久久久久久7| 日韩不卡一区二区三区视频在线| a级片在线免费高清观看视频| 9191精品国产免费久久| 2022亚洲国产成人精品| 国产亚洲av片在线观看秒播厂| 精品一区在线观看国产| 亚洲经典国产精华液单| 国产乱来视频区| 国产成人免费观看mmmm| 亚洲男人天堂网一区| 国产黄色视频一区二区在线观看| 大片免费播放器 马上看| 成人影院久久| 啦啦啦在线观看免费高清www| 国产极品粉嫩免费观看在线| 成人手机av| 一级片免费观看大全| 久久精品人人爽人人爽视色| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲三级黄色毛片| 九色亚洲精品在线播放| 咕卡用的链子| 啦啦啦视频在线资源免费观看| 久久午夜综合久久蜜桃| 亚洲精品国产一区二区精华液| 美女国产视频在线观看| 国产熟女午夜一区二区三区| 亚洲人成电影观看| 亚洲精品aⅴ在线观看| 大陆偷拍与自拍| 国产在线免费精品| 深夜精品福利| 午夜免费男女啪啪视频观看| 欧美97在线视频| 亚洲欧美色中文字幕在线| 亚洲精品一区蜜桃| 亚洲精品国产av蜜桃| 亚洲av综合色区一区| 九色亚洲精品在线播放| www日本在线高清视频| 亚洲国产最新在线播放| 卡戴珊不雅视频在线播放| 9色porny在线观看| 国产人伦9x9x在线观看 | 成年av动漫网址| 女人被躁到高潮嗷嗷叫费观| 午夜久久久在线观看| 欧美精品一区二区大全| 欧美97在线视频| 热99国产精品久久久久久7| 国产成人精品婷婷| 色网站视频免费| 新久久久久国产一级毛片| 美女福利国产在线| 免费在线观看黄色视频的| 国产在线视频一区二区| 纯流量卡能插随身wifi吗| 亚洲,一卡二卡三卡| 黄片无遮挡物在线观看| 婷婷色综合大香蕉| 日韩伦理黄色片| 国产国语露脸激情在线看| 少妇被粗大猛烈的视频| 日韩伦理黄色片| 亚洲av.av天堂| 18禁裸乳无遮挡动漫免费视频| 极品人妻少妇av视频| 国产精品亚洲av一区麻豆 | 国产极品天堂在线| 超色免费av| 国产精品成人在线| av免费观看日本| 成人国产麻豆网| 精品国产一区二区久久| 色婷婷av一区二区三区视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av.av天堂| 亚洲国产欧美在线一区| 色哟哟·www| 制服人妻中文乱码| 黄色怎么调成土黄色| 免费观看性生交大片5| av卡一久久| 亚洲精品国产av成人精品| 五月天丁香电影| 国产精品久久久久久精品电影小说| 人妻少妇偷人精品九色| 777久久人妻少妇嫩草av网站| 一边摸一边做爽爽视频免费| 天堂俺去俺来也www色官网| 婷婷成人精品国产| 亚洲国产成人一精品久久久| 熟妇人妻不卡中文字幕| 久久久国产欧美日韩av| 欧美精品高潮呻吟av久久| 男人添女人高潮全过程视频| 久久婷婷青草| 日韩在线高清观看一区二区三区| 在线精品无人区一区二区三| 黄色一级大片看看| 香蕉丝袜av| 在线看a的网站| 汤姆久久久久久久影院中文字幕| 欧美精品一区二区免费开放| 久久久久久伊人网av| 一级片'在线观看视频| 成人影院久久| 久久97久久精品| 午夜91福利影院| 99久国产av精品国产电影| 国产黄色视频一区二区在线观看| 久久久久久免费高清国产稀缺| 国产欧美日韩一区二区三区在线| 久久婷婷青草| 91精品三级在线观看| 亚洲,欧美,日韩| 国产白丝娇喘喷水9色精品| 永久网站在线| 亚洲精品自拍成人| 欧美日韩一级在线毛片| 国产成人91sexporn| xxx大片免费视频| 激情五月婷婷亚洲| 人人妻人人添人人爽欧美一区卜| 大码成人一级视频| 久久99精品国语久久久| 美女高潮到喷水免费观看| 国产色婷婷99| 女的被弄到高潮叫床怎么办| 成年人午夜在线观看视频| 不卡av一区二区三区| 国产精品一区二区在线不卡| 中文欧美无线码| 国产有黄有色有爽视频| 日韩精品免费视频一区二区三区| 麻豆乱淫一区二区| 性色avwww在线观看| 欧美精品av麻豆av| 久久久精品国产亚洲av高清涩受| 成年女人毛片免费观看观看9 | 五月伊人婷婷丁香| 日韩中字成人| 成年女人毛片免费观看观看9 | 欧美日韩亚洲国产一区二区在线观看 | 久久久久国产网址| 五月开心婷婷网| av.在线天堂| 欧美亚洲 丝袜 人妻 在线| 精品国产乱码久久久久久小说| 久久久精品国产亚洲av高清涩受| 精品久久蜜臀av无| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利视频精品| 少妇熟女欧美另类| 极品少妇高潮喷水抽搐| 宅男免费午夜| 国产精品一区二区在线观看99| 亚洲国产精品一区二区三区在线| 亚洲精品,欧美精品| 日韩中字成人| 午夜激情久久久久久久| 国精品久久久久久国模美| 亚洲四区av| 成人二区视频| 国产成人精品一,二区| 欧美激情 高清一区二区三区| 久久久久久久久久人人人人人人| 亚洲精品日韩在线中文字幕| 精品国产乱码久久久久久小说| 久久久久久伊人网av| 青春草亚洲视频在线观看| 国产免费视频播放在线视频| 国产精品99久久99久久久不卡 | 欧美日韩精品网址| 亚洲国产成人一精品久久久| 日韩一区二区视频免费看| 黄色一级大片看看| 男女免费视频国产| 街头女战士在线观看网站| av片东京热男人的天堂| 搡老乐熟女国产| 在线 av 中文字幕| 国产精品蜜桃在线观看| 青春草视频在线免费观看| 久久久久国产网址| 母亲3免费完整高清在线观看 | 99久久综合免费| 女性被躁到高潮视频| 一区福利在线观看| 中文天堂在线官网| 欧美日韩亚洲国产一区二区在线观看 | 免费黄频网站在线观看国产| 午夜免费男女啪啪视频观看| 色网站视频免费| 只有这里有精品99| 亚洲精品乱久久久久久| 观看美女的网站| 亚洲精品在线美女| 高清视频免费观看一区二区| 亚洲精品,欧美精品| 欧美日韩一级在线毛片| av网站在线播放免费| 久久久久久久精品精品| 下体分泌物呈黄色| 亚洲美女搞黄在线观看| 国产精品国产av在线观看| 叶爱在线成人免费视频播放| 如日韩欧美国产精品一区二区三区| 中文精品一卡2卡3卡4更新| 久久国产亚洲av麻豆专区| 最近2019中文字幕mv第一页| 亚洲国产毛片av蜜桃av| 日本黄色日本黄色录像| 久久99一区二区三区| 亚洲欧美成人综合另类久久久| 国产高清国产精品国产三级| 高清欧美精品videossex| 久久精品久久久久久久性| videos熟女内射| 99国产精品免费福利视频| 成人手机av| 日韩av免费高清视频| av电影中文网址| 国产精品免费视频内射| 一级a爱视频在线免费观看| 精品亚洲乱码少妇综合久久| 国产伦理片在线播放av一区| 成年动漫av网址| 亚洲,欧美,日韩| 少妇被粗大的猛进出69影院| 亚洲欧美清纯卡通| 夜夜骑夜夜射夜夜干| 精品亚洲乱码少妇综合久久| 亚洲,欧美,日韩| 欧美精品一区二区大全| 亚洲av男天堂| 王馨瑶露胸无遮挡在线观看| 蜜桃国产av成人99| 国产精品不卡视频一区二区| 大话2 男鬼变身卡| 国产一区二区在线观看av| 国产精品一区二区在线不卡| 女性被躁到高潮视频| 自线自在国产av| 看免费成人av毛片| a级毛片在线看网站| 各种免费的搞黄视频| 人人妻人人添人人爽欧美一区卜| 曰老女人黄片| 国产野战对白在线观看| 久久精品aⅴ一区二区三区四区 | 一区二区三区精品91| 黑人猛操日本美女一级片| 欧美精品亚洲一区二区| 91成人精品电影| 国产av码专区亚洲av| 少妇的丰满在线观看| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 亚洲激情五月婷婷啪啪| 精品一区在线观看国产| 久久久久国产精品人妻一区二区| 日韩电影二区| 天天操日日干夜夜撸| 亚洲天堂av无毛| 卡戴珊不雅视频在线播放| 天天躁夜夜躁狠狠躁躁| 国产探花极品一区二区| 色哟哟·www| 日本爱情动作片www.在线观看| 国产精品久久久久久精品古装| 人妻一区二区av| 18在线观看网站| 高清av免费在线| 午夜91福利影院| 999精品在线视频| 国产精品久久久av美女十八| 如日韩欧美国产精品一区二区三区| 久久久久网色| 五月开心婷婷网| 丰满乱子伦码专区| 国产乱来视频区| 欧美日韩综合久久久久久| 在线观看国产h片| 国产成人欧美| 天天躁夜夜躁狠狠久久av| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 久久这里只有精品19| 久久精品久久久久久久性| 肉色欧美久久久久久久蜜桃| 国产精品女同一区二区软件| 亚洲国产成人一精品久久久| 国产黄色免费在线视频| 性色avwww在线观看| 亚洲在久久综合| 欧美国产精品一级二级三级| 纯流量卡能插随身wifi吗| 最近最新中文字幕大全免费视频 | 亚洲成av片中文字幕在线观看 | av国产久精品久网站免费入址| 欧美在线黄色| a级毛片在线看网站| 高清不卡的av网站| av在线老鸭窝| 香蕉国产在线看| 丝袜美足系列| 国产亚洲精品第一综合不卡| 国产精品国产三级专区第一集| 卡戴珊不雅视频在线播放| 午夜福利乱码中文字幕| 久久精品久久久久久噜噜老黄| 曰老女人黄片| 久久鲁丝午夜福利片| 亚洲av福利一区| 99九九在线精品视频| 久久国产精品男人的天堂亚洲| 天天躁夜夜躁狠狠久久av| 久久青草综合色| 国产亚洲精品第一综合不卡| 高清黄色对白视频在线免费看| 久久97久久精品| 高清在线视频一区二区三区| 久久久亚洲精品成人影院| 国产成人免费无遮挡视频| 女人高潮潮喷娇喘18禁视频| 欧美精品人与动牲交sv欧美| 免费观看无遮挡的男女| 亚洲少妇的诱惑av| av网站免费在线观看视频| 观看美女的网站| 久久婷婷青草| 国产极品粉嫩免费观看在线| 中文字幕av电影在线播放| 香蕉精品网在线| 最近中文字幕2019免费版| 多毛熟女@视频| 999久久久国产精品视频| 国产精品 国内视频| av国产精品久久久久影院| 国产xxxxx性猛交| 国产97色在线日韩免费| 国产精品蜜桃在线观看| 91国产中文字幕| 最近2019中文字幕mv第一页| 久久久久国产网址| av一本久久久久| 亚洲国产欧美日韩在线播放| www.精华液| 夫妻性生交免费视频一级片| 色哟哟·www| 免费黄网站久久成人精品| 精品福利永久在线观看| 9色porny在线观看| 十八禁高潮呻吟视频| 久久久久久久久久久久大奶| 欧美日韩视频高清一区二区三区二| 国产精品蜜桃在线观看| 久久久久久久久久久久大奶| 久久午夜福利片| 99re6热这里在线精品视频| 色哟哟·www| 欧美成人午夜免费资源| 婷婷色麻豆天堂久久| 一区二区三区四区激情视频| 老女人水多毛片| 国产一区二区三区av在线| 欧美日韩国产mv在线观看视频| 久久久精品免费免费高清| 最黄视频免费看| 五月天丁香电影| 国产探花极品一区二区| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 日本免费在线观看一区| 一二三四中文在线观看免费高清| 国产淫语在线视频| 久久久久精品久久久久真实原创| 亚洲国产色片| 精品酒店卫生间| 亚洲欧洲精品一区二区精品久久久 | 黄色怎么调成土黄色| 精品99又大又爽又粗少妇毛片| 亚洲欧美精品自产自拍| 久久久久人妻精品一区果冻| 少妇的逼水好多| 日韩三级伦理在线观看| 亚洲国产看品久久| 777米奇影视久久| 国产乱人偷精品视频| 国产极品粉嫩免费观看在线| 三上悠亚av全集在线观看| 国产熟女午夜一区二区三区| 久久久久精品性色| 午夜福利在线免费观看网站| 观看美女的网站| 少妇的逼水好多| 自拍欧美九色日韩亚洲蝌蚪91| 国产野战对白在线观看| 少妇人妻久久综合中文| 久久精品国产亚洲av天美| 一区在线观看完整版| 国产日韩欧美视频二区| 香蕉丝袜av| 亚洲国产欧美网| 亚洲第一青青草原| 国产精品人妻久久久影院| 国产精品久久久久久精品古装| 伊人亚洲综合成人网| 亚洲四区av| 成年人免费黄色播放视频| 国产男女内射视频| 美女视频免费永久观看网站| 国产成人精品久久二区二区91 | 国产成人精品久久久久久| 黄色毛片三级朝国网站| 国产精品二区激情视频| 国产一区亚洲一区在线观看| 高清在线视频一区二区三区| 亚洲av成人精品一二三区| 精品99又大又爽又粗少妇毛片| 日本午夜av视频| 三级国产精品片| 亚洲国产精品一区二区三区在线| 一区二区三区乱码不卡18| 日韩熟女老妇一区二区性免费视频| 9色porny在线观看| 亚洲久久久国产精品| 成人漫画全彩无遮挡| 欧美变态另类bdsm刘玥| 免费av中文字幕在线| 高清在线视频一区二区三区| 少妇人妻 视频| 中文字幕制服av| 国产精品成人在线| 日韩一区二区三区影片| 亚洲精品aⅴ在线观看| videossex国产| 国产精品二区激情视频| www.熟女人妻精品国产| 久热这里只有精品99| 桃花免费在线播放| 中文字幕制服av| 国产精品成人在线| 久久久久国产精品人妻一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产成人免费无遮挡视频| 午夜免费男女啪啪视频观看| 久久久国产一区二区| 91成人精品电影| 亚洲综合色惰| 国产爽快片一区二区三区| 亚洲熟女精品中文字幕| 韩国av在线不卡| 人成视频在线观看免费观看| 国语对白做爰xxxⅹ性视频网站| 大香蕉久久网| 最近中文字幕高清免费大全6| 免费看不卡的av| 欧美精品av麻豆av| 1024香蕉在线观看| 国产高清国产精品国产三级| 久久av网站| 免费久久久久久久精品成人欧美视频| 一级黄片播放器| 午夜日本视频在线| 建设人人有责人人尽责人人享有的| 超色免费av| 一级a爱视频在线免费观看| 国产成人精品无人区| 国产精品久久久久成人av| 中文字幕另类日韩欧美亚洲嫩草| 国语对白做爰xxxⅹ性视频网站| 国产免费现黄频在线看| 人人澡人人妻人| 国产伦理片在线播放av一区| 国产日韩一区二区三区精品不卡| 亚洲人成网站在线观看播放| 国产精品人妻久久久影院| 99久久人妻综合| 中文欧美无线码| 啦啦啦啦在线视频资源| 成人手机av| 亚洲国产欧美日韩在线播放| 一区在线观看完整版| 老熟女久久久| 亚洲av男天堂| 精品一区二区三卡| 亚洲国产精品一区三区| 欧美日本中文国产一区发布| 校园人妻丝袜中文字幕| 女性生殖器流出的白浆| 精品国产国语对白av| 久久女婷五月综合色啪小说| 毛片一级片免费看久久久久| 成人免费观看视频高清| 国产av国产精品国产| www.av在线官网国产| 午夜福利视频精品| 国产男女超爽视频在线观看| 免费在线观看视频国产中文字幕亚洲 | 又大又黄又爽视频免费| 亚洲精品中文字幕在线视频| 亚洲天堂av无毛| 欧美亚洲日本最大视频资源| 亚洲美女黄色视频免费看| 亚洲成人一二三区av| 欧美精品一区二区免费开放| 丝袜在线中文字幕| av免费观看日本| 久热久热在线精品观看|