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

    基于NEX-GDDP數(shù)據(jù)集的青藏高原牧區(qū)雪災(zāi)風(fēng)險預(yù)估*

    2022-06-27 08:46:06陳虹舉楊建平丁永建賀青山王彥霞葛秋伶
    災(zāi)害學(xué) 2022年2期
    關(guān)鍵詞:雪深雪災(zāi)牧區(qū)

    陳虹舉,楊建平,丁永建,2,4,賀青山,3,冀 欽,王彥霞,3,唐 凡,3,葛秋伶,3

    (1.中國科學(xué)院西北生態(tài)環(huán)境資源研究院 冰凍圈科學(xué)國家重點實驗室,甘肅 蘭州 730000;2.中國科學(xué)院西北生態(tài)環(huán)境資源研究院 內(nèi)陸河流域生態(tài)水文重點實驗室,甘肅 蘭州 730000;3.中國科學(xué)院大學(xué),北京100049;4.中國科學(xué)院大學(xué) 資源與環(huán)境學(xué)院,北京 100049)

    青藏高原是我國高寒畜牧區(qū)。高亢的地勢、嚴寒的氣候、冰川、凍土、積雪等共存的生態(tài)環(huán)境形成了較為獨特的牧業(yè)經(jīng)濟的本底資源,作為區(qū)域經(jīng)濟體的重要組成部分,也使其直接暴露于冰凍圈環(huán)境之中。青藏高原高寒草地生長季節(jié)短,草地生產(chǎn)力與承載力較低[1],加之牧民世代放牧為生,受教育水平低[2],致使畜牧業(yè)經(jīng)濟這一承險體脆弱且適應(yīng)能力較低,極易受到積雪災(zāi)害影響。該青藏高原高寒牧區(qū)雪災(zāi)過程表現(xiàn)為:積雪過厚,牧草被大雪掩埋,靠牧草為生的家畜因吃不到草,凍餓而死[3]。因而雪災(zāi)風(fēng)險評估成為青藏高原災(zāi)害防治中必不可少的一項基礎(chǔ)性工作,特別是準確的預(yù)估未來青藏高原雪災(zāi)風(fēng)險對當?shù)卣蜕鐣茖W(xué)準確地指定防災(zāi)措施、及時開展雪災(zāi)緊急援救以及災(zāi)害恢復(fù)有著重要指導(dǎo)作用。

    國內(nèi)外學(xué)者針對牧區(qū)雪災(zāi)研究可以分為以下幾個方面:①雪災(zāi)危險性研究,KlUVER[4]利用歷史氣象資料分析美國本土降雪頻率與趨勢;李海紅等[5]2006在中國牧區(qū)雪災(zāi)等級指標研究中,依據(jù)積雪掩埋牧草程度、積雪持續(xù)日數(shù)和積雪面積三個方面,將雪災(zāi)災(zāi)情分為輕災(zāi)、中災(zāi)、重災(zāi)和特大災(zāi)害四級;巴桑等[6]基于多種源數(shù)據(jù)對西藏地區(qū)積雪變化趨勢進行分析;②承災(zāi)體對雪災(zāi)脆弱性研究,解偉等[7]通過積雪深度與牲畜死亡率之間關(guān)系,厘定了雪災(zāi)保險率;③雪災(zāi)的適應(yīng)性與災(zāi)后救援研究,趙富昌等[8]通過對青藏高原歷史雪災(zāi)分析,提出如何合理載青藏高原雪災(zāi)高風(fēng)險區(qū)建立草料儲備庫;趙霞[9]則利用GIS方法對錫林郭勒盟雪災(zāi)救助區(qū)進行區(qū)劃;④雪災(zāi)預(yù)警模型研究,王瑋[10]結(jié)合RS和GIS對青藏高原積雪動態(tài)變化進行研究,并構(gòu)建了雪災(zāi)危害等級預(yù)警模型;⑤建立雪災(zāi)風(fēng)險指標模型并對雪災(zāi)區(qū)劃分析,何永清等[11]通過選取牲畜數(shù)量、雪災(zāi)次數(shù)、牧草面積、人均GDP以及牧草產(chǎn)量基礎(chǔ)上對青海省雪災(zāi)風(fēng)險評價并對風(fēng)險區(qū)進行區(qū)劃。TACHIIRI等[12]在考慮前1年NDVI和社會普查牲畜死亡量基礎(chǔ)上,通過樹形結(jié)構(gòu)模型對蒙古雪災(zāi)風(fēng)險進行研究??傮w來說,牧區(qū)雪災(zāi)作為多因子氣象災(zāi)害,危險性程受多種因素影響,如積雪深度、積雪持續(xù)時間、積雪覆蓋范圍、積雪時候的氣溫以及風(fēng)速等,因此其識別難度大。雪災(zāi)脆弱性研究中只考慮雪災(zāi)危險性的一種因素與牲畜損失建立關(guān)系,其結(jié)果并不能夠反映牲畜脆弱性。且目前對雪災(zāi)的研究是集中于歷史時段以及未來短期預(yù)警研究,對未來牧區(qū)長期的雪災(zāi)風(fēng)險定量預(yù)估涉及較少。

    因此,此次研究針對自然狀況以及傳統(tǒng)雪災(zāi)防御措施下的牲畜受災(zāi)狀況,定量預(yù)估未來RCPs情景下牧區(qū)雪災(zāi)風(fēng)險。首先選取雪災(zāi)期間積雪平均深度、雪災(zāi)持續(xù)時間、雪災(zāi)期間的日最高氣溫均值以及雪災(zāi)期間日最低氣溫均值四個因子對雪災(zāi)危險度進行識別。其次建立牲畜死亡率和雪災(zāi)危險指數(shù)的脆弱性曲線關(guān)系。以NEX-GDDP數(shù)據(jù)、中國逐日雪深模擬預(yù)估數(shù)據(jù)集、牲畜統(tǒng)計數(shù)據(jù)以及1980年以來青藏高原草地生產(chǎn)力數(shù)據(jù)為數(shù)據(jù)基礎(chǔ),運用“風(fēng)險=危險性×暴露度×脆弱性”方法[13],對青藏高原牧區(qū)未來雪災(zāi)風(fēng)險進行預(yù)估。其中NEX-GDDP數(shù)據(jù)是由美國航空航天局2015年發(fā)布的首套基于CMIP5的多模式高分辨率數(shù)據(jù)集NASA Earth Exchange/Global Daily Downscaled Projections(NEX-GDDP),它利用統(tǒng)計降尺度方法,將21個CMIP5模式模擬得到的歷史時期(1950—2005年)和兩種情景下(RCP4.5和RCP8.5)預(yù)估時段(2006—2100年)的逐日降水、最高氣溫和最低氣溫數(shù)據(jù)轉(zhuǎn)換成0.25°×0.25°空間分辨率[14]。

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)及其來源

    中國逐日雪深模擬預(yù)估數(shù)據(jù)集:該數(shù)據(jù)是以NEX-GDDP數(shù)據(jù)集為依托,利用人工神經(jīng)網(wǎng)絡(luò)模型模擬獲取,包括基準時段(1986—2005年)、未來預(yù)估數(shù)據(jù)(2016—2065年)兩個時段,RCP4.5和RCP8.5兩種氣候情景的逐日雪深數(shù)據(jù),空間分辨率為0.25°×0.25°。數(shù)據(jù)來源于國家青藏高原科學(xué)數(shù)據(jù)中心(DOI:10.11888/Snow.tpdc.271636.CSTR:18406.11.Snow.tpdc.271636)。

    NEX-GDDP數(shù)據(jù)集:該數(shù)據(jù)包括基準時段(1986—2005年)和未來預(yù)估數(shù)據(jù)(2006—2100年), RCP4.5和RCP8.5兩種情景的逐日降水、最高氣溫和最低氣溫預(yù)估數(shù)據(jù),空間分辨率為0.25°×0.25°,數(shù)據(jù)地址為https://nex.nasa.gov/nex/projects/1356/。

    1980年以來青藏高原草地生產(chǎn)力數(shù)據(jù)(莫興國):該數(shù)據(jù)是將NPP換算為干物質(zhì),再根據(jù)根冠比估算了干草產(chǎn)量,空間分辨率1 km×1 km,數(shù)據(jù)來源于國家青藏高原科學(xué)數(shù)據(jù)中心DOI:10.11888/Ecolo.tpdc.270430。

    統(tǒng)計年鑒數(shù)據(jù):甘肅、青海、云南、新疆、四川以及西藏地區(qū)市級統(tǒng)計年鑒,時序1986—2005年,此次研究摘錄了縣級行政單元內(nèi)牧區(qū)各種牲畜數(shù)量數(shù)據(jù)。

    中國雪深長時間序列數(shù)據(jù)集:該數(shù)據(jù)集是以遙感數(shù)據(jù)為基礎(chǔ)反演獲取,其包含1978年10月24日到2012年12月31日逐日的中國范圍的積雪厚度分布數(shù)據(jù),其空間分辨率為0.25°×0.25°。數(shù)據(jù)來源于國家青藏高原科學(xué)數(shù)據(jù)中心[15](2019)(DOI: 10.3972/westdc.001.2014.db)。

    氣象災(zāi)害統(tǒng)計數(shù)據(jù):氣象災(zāi)害統(tǒng)計數(shù)據(jù)來自于各省的氣象災(zāi)害大典[16-19],主要記錄了2000年以前發(fā)生的氣象災(zāi)害資料,文中統(tǒng)計了雪災(zāi)導(dǎo)致青藏高原各地區(qū)牲畜死亡數(shù)量。

    站點氣象數(shù)據(jù):文中使用了青藏高原范圍內(nèi)的地面氣象臺站數(shù)據(jù),其包括日最高氣溫、日最低氣溫和積雪深度數(shù)據(jù),來源于國家氣象科學(xué)數(shù)據(jù)中心(http://www.cma.gov.cn/),時間序列為1986—2005年。

    1.2 研究方法

    1.2.1 中國逐日雪深模擬預(yù)估數(shù)據(jù)集對青藏高原降雪模擬能力評估方法

    中國逐日雪深模擬預(yù)估數(shù)據(jù)中不同模式雪深的模擬結(jié)果在同一區(qū)域精度差異較大,為了明晰此次模擬數(shù)據(jù)和遙感模擬逐日雪深數(shù)據(jù)差異,文中在青藏高原區(qū)域?qū)δM的不同模式下雪深數(shù)據(jù)與站點觀測雪深數(shù)據(jù)進行比較分析,同時為了比較用模式模擬的雪深數(shù)據(jù)和應(yīng)用廣泛的遙感模擬雪深數(shù)據(jù)差異,因此文中加入了中國雪深長時間序列數(shù)據(jù)集進行對比。評估數(shù)據(jù)所選擇時間段為1986—2005年。具體過程如下:

    以青藏高原內(nèi)所有氣象站點所在位置為基準,分別統(tǒng)計三種積雪深度數(shù)據(jù)的站點年累積積雪深度。然后用氣象站數(shù)據(jù)作為觀測值,分析其余兩種雪深數(shù)據(jù)的模擬能力,評估中采用觀測值和兩種雪深數(shù)據(jù)的中心化均方根誤差和相關(guān)系數(shù)作為統(tǒng)計指標,依據(jù)這兩個指標的統(tǒng)計結(jié)果,評估兩種雪深數(shù)據(jù)對青藏高原積雪的刻畫能力,并評選出中國逐日雪深模擬預(yù)估數(shù)據(jù)集中最優(yōu)模式。為了確定最優(yōu)模式可靠性,進一步統(tǒng)計最優(yōu)模式、站點雪深數(shù)據(jù)以及遙感反演雪深數(shù)據(jù)的所有站點逐日雪深概率密度函數(shù),并進行比較分析。

    (1)

    (2)

    (3)

    (4)

    1.2.2 雪災(zāi)危險性指數(shù)的計算方法

    雪災(zāi)危險性是指不同等級降雪事件發(fā)生可能性。牧區(qū)雪災(zāi)作為多因子氣象自然災(zāi)害,其危險性受多種因素影響,如積雪深度、積雪持續(xù)時間以及積雪時段內(nèi)氣溫。對雪災(zāi)危險性識別,不僅要了解不同氣象因子強弱等級,而且要明晰不同因子權(quán)重。因此,首先通過文獻查閱、牧區(qū)雪災(zāi)劃分等級標準[3]以及專家咨詢,分析不同因子對雪災(zāi)的重要性,并進行比較,在此基礎(chǔ)上,結(jié)合層次分析法的判斷矩陣構(gòu)造,設(shè)定了不同氣象因子的權(quán)重(表1)。得到雪災(zāi)危險性模型:

    H=0.449 2×X1+0.301 1×X2-0.100 2×X3-0.149 5×X4。

    (5)

    式中:H為雪災(zāi)危險性指數(shù),X1為積雪期間積雪平均深度,X2積雪的持續(xù)時間,X3為積雪期間的日最高氣溫均值,X4為積雪期間日最低氣溫均值。其中H值越高,雪災(zāi)危險性越大。在確定雪災(zāi)危險性指數(shù)的計算方法的基礎(chǔ)上,再利用中國逐日雪深模擬預(yù)估數(shù)據(jù)集和NEX-GDDP數(shù)據(jù)集對青藏高原區(qū)域內(nèi)的雪災(zāi)危險性指數(shù)進行識別。

    1.2.3 牲畜空間暴露量網(wǎng)格化方法

    牲畜暴露量是指在孕災(zāi)環(huán)境中可能受雪災(zāi)影響并產(chǎn)生損失的牲畜數(shù)量及其空間分布。在以往自然災(zāi)害損失評估中,承載體的空間展布大多是以行政單元進行,認為承載體在行政區(qū)內(nèi)是均勻分布,而忽略了一些因素下的承載體分布的不勻承性,如本次研究牲畜數(shù)量單位網(wǎng)格數(shù)量多少與各個縣的草地面積以及各個網(wǎng)格產(chǎn)草量有著密切的關(guān)系,因此本次研究利用GIS技術(shù),以柵格數(shù)據(jù)形式將牲畜的空間分布進行展示。為了與危險度網(wǎng)格單元一致,同樣選擇0.25°×0.25°網(wǎng)格單元。具體過程如下。

    實際的牲畜空間分布:

    (6)

    式中:ADs為縣域內(nèi)所包含的某一柵格單元內(nèi)的牲畜數(shù)量;GW柵格單元內(nèi)產(chǎn)草量;SGW為縣域內(nèi)所有柵格單元總產(chǎn)草量;SAD為此縣牲畜實際擁有總數(shù)量。

    “以草定畜”牲畜空間分布:

    (7)

    TLADs=TSADs×scale。

    (8)

    式中:TSADs為某一柵格單元內(nèi)羊的承載數(shù)量;GW為柵格單元內(nèi)產(chǎn)草量;scale為此縣牛和羊數(shù)量的比值;TLADs大型牲畜數(shù)量。其中羊的食草量為1.38 kg/d,大型牲畜食草量為4.5 kg/d,1年總天數(shù)為365 d,草地利用率為0.46[20]。

    1.2.4 草原雪災(zāi)牲畜損失率曲線建立

    脆弱性是指一定致災(zāi)因子強度下,承載體可能遭受損失的程度。而在青藏高原區(qū)域內(nèi)不同雪災(zāi)發(fā)生地牲畜總數(shù)量、人口、社會經(jīng)濟發(fā)展水平不一,如若用絕對牲畜受雪災(zāi)損失數(shù)量來評估雪災(zāi),其結(jié)果必定不能夠真實的情況,因此此次研究選擇相對損失指標“損失率”來構(gòu)建雪災(zāi)的損失曲線。其具體過程如下:

    (9)

    式中:DR為牲畜死亡率;DsAn為草原雪災(zāi)牲畜損失數(shù)量;Aan為當年牲畜總數(shù)。

    在統(tǒng)計歷史時期不同雪災(zāi)發(fā)生情況下牲畜死亡率同時,在根據(jù)危險性計算公式(1),計算出發(fā)生雪災(zāi)時期的雪災(zāi)危險性指數(shù),然后建立雪災(zāi)危險性程度和牲畜的死亡率的相關(guān)系曲線。

    DR=f(H)。

    (10)

    式中:DR為雪災(zāi)期間牲畜損失率;H為發(fā)生雪災(zāi)期間雪災(zāi)危險性指數(shù);f為相關(guān)函數(shù)。

    1.2.5 青藏高原牧區(qū)雪災(zāi)定量預(yù)估方法

    研究方法采用了2004年聯(lián)合國開發(fā)計劃署(UNDP)提出的風(fēng)險研究方法,對未來兩個時期(2016—2035年和2046—2065年),兩種氣候情景(RCP4.5和RCP8.5)下牲畜受雪災(zāi)風(fēng)險進行預(yù)估。具體表達如下:

    風(fēng)險=危險性×暴露度×脆弱性。

    (11)

    1.3 技術(shù)路線

    牧區(qū)雪災(zāi)從形成機理角度來看,主要分為雪災(zāi)危險性(即積雪深度、積雪持續(xù)時間以及降雪時段內(nèi)溫度等多因子綜合形成的致災(zāi)因子危險性)、牧區(qū)牲畜暴露度(即承災(zāi)體在時間和空間可能與致災(zāi)因子的接觸、重合量)和脆弱性(即牧區(qū)牲畜雪災(zāi)的脆弱性曲線,不同致災(zāi)因子強度對承災(zāi)體作用造成損失的相關(guān)曲線)三個部分。其中牧區(qū)雪災(zāi)危險性指危險性事件發(fā)生規(guī)模、程度及其可能性,然而雪災(zāi)危險性是多因子綜合而成,相比于單因子災(zāi)害識別難度較大,例如:地震是單因子災(zāi)害,可直接利用地震釋放能量大小量化地震危險性大小(地震等級),雪災(zāi)作為多因子氣象災(zāi)害,牧區(qū)雪災(zāi)危險不能僅根據(jù)積雪深度或者積雪時長來定義雪災(zāi)危險性強度,危險性識別應(yīng)該考慮多種氣象因子影響;牧區(qū)牲畜雪災(zāi)暴露度則是指在承載體在可能發(fā)生雪災(zāi)危險性區(qū)域內(nèi)的暴露量以及空間分布(青藏高原范圍內(nèi)都有發(fā)生雪災(zāi)危險可能,因而本文以整個青藏高原內(nèi)牲畜暴露量作為研究對象);牧區(qū)牲畜雪災(zāi)脆弱性是指一定雪災(zāi)危險強度下,青藏高原內(nèi)某一區(qū)域牲畜可能遭受損失的多少。在總結(jié)前人研究基礎(chǔ)上,結(jié)合青藏高原牧區(qū)雪災(zāi)實際情況,確定了本文研究的研究方案(圖1)。

    表1 不同氣象因子對牧區(qū)雪災(zāi)危險性的重要程度對比矩陣

    2 結(jié)果分析

    2.1 中國逐日雪深模擬預(yù)估數(shù)據(jù)對青藏高原降雪模擬能力評估結(jié)果

    在青藏高原范圍內(nèi),遙感模擬雪深與站點實測的均方根誤差遠大于模式模擬的雪深與站點實測的均方根誤差,就相關(guān)系數(shù)而言,CESM1-BGC模式模擬的積雪深度數(shù)據(jù)與站點實測值的相關(guān)性系數(shù)大小接近于遙感模擬雪深的數(shù)據(jù)與站點實測值的相關(guān)系數(shù),整體來說CESM1-BGC模式模擬的積雪深度數(shù)據(jù)最接近于站點積雪深度觀測值(表2和圖2)。

    為了驗證用ESM1-BGC模式模擬的雪深數(shù)據(jù)與站點實測數(shù)據(jù)、微波模擬的雪深數(shù)據(jù)空間差異,為此基于站點點位對三種雪深數(shù)據(jù)的年累積積雪時間和年累計積雪深度進行統(tǒng)計并進行空間比較,結(jié)果顯示(圖3):無論是累積積雪天數(shù),還是累積積雪深度,在青藏高原中部、南部、北部以及東部邊緣區(qū)域此次模擬結(jié)果比遙感模擬的結(jié)果更加接近于站點實測值,而在青藏高原西部以及西南區(qū)域,此次模擬的逐日雪深數(shù)據(jù)與站點觀測結(jié)果差異比較大,整體看ESM1-BGC模式模擬的雪深數(shù)據(jù)更接近于站點實測數(shù)據(jù)。為了進一步明晰ESM1-BGC模式模擬的日雪深數(shù)據(jù)分布合理性,統(tǒng)計了概率密度函數(shù),統(tǒng)計結(jié)果顯示:CESM1-BGC模式模擬的逐日雪深概率密度函和遙感模擬的逐日雪深概率密度函相比較,CESM1-BGC模式模擬的逐日雪深概率密度函分布更接近于站點實測值(圖4)。

    圖1 研究思路與技術(shù)路線圖

    圖2 幾種雪深數(shù)據(jù)年累積積雪深度比較

    表2 青藏高原地區(qū)21個模式下模擬積雪深度數(shù)據(jù)和站點積雪深度的比較

    圖3 青藏高原區(qū)域基于站點點位三種逐日雪深數(shù)據(jù)空間差異比較(審圖號:GS(2016)1552號,底圖無修改,下同)

    圖4 青藏高原區(qū)域內(nèi)全部氣象站點點位逐日雪深數(shù)據(jù)概率密度函數(shù)

    圖5 青藏高原區(qū)域內(nèi)雪災(zāi)危險性指數(shù)時空分布

    表3 青藏高原區(qū)雪災(zāi)危險性指數(shù)統(tǒng)計結(jié)果

    圖6 牲畜空間暴露量分布圖

    2.2 雪災(zāi)危險性分析

    利用公式(1),運用CESM1-BGC模式下逐日雪深數(shù)據(jù)對雪災(zāi)危險性進行識別,結(jié)果顯示從空間分布來看,雪災(zāi)危險性指數(shù)較高的區(qū)域主要分布在藏北高原、岡底斯山脈沿線、昆侖山脈西段沿線、祁連山脈沿線、三江源區(qū)域和橫斷山脈區(qū)域(圖5)。危險指數(shù)強度均值來看,相比較于歷史時期,RCP4.5情景下,未來遠期雪災(zāi)危險性呈現(xiàn)出增加趨勢,其余情況下則減弱;發(fā)生雪災(zāi)危險性區(qū)域面積占比變化來看,相比于歷史時期,RCP4.5情景下未來近期、未來遠期和RCP8.5情景下未來近期、未來遠期發(fā)生雪災(zāi)危險性的范圍減少6%、11%、6%和14%,縮減區(qū)域處于青藏高原中部和祁連山地區(qū)。在未來RCP8.5情景下,發(fā)生雪災(zāi)危險性區(qū)域面積占比減少都比較明顯,但是發(fā)生區(qū)的雪災(zāi)危險度強度減弱并不明顯,有的時段甚至可能增強(圖5和表3)。

    2.3 牲畜空間暴露量分布

    根據(jù)青藏高原草地生產(chǎn)力空間分布不均勻,對牲畜數(shù)量進行空間網(wǎng)格化展布,結(jié)果(圖6和表4)顯示: 1986—2005年實際牲畜量遠高于2000年草原可載畜量。說明2000年前后,青藏高原牧區(qū)處于嚴重的超載狀態(tài),從空間分布來看,超載區(qū)域主要分布于以拉薩市為中心的周邊區(qū)域、以西寧為中心周邊區(qū)域以及橫斷山脈區(qū)域。2017年草原可載畜量相比于2000年增加11%,其原因是青藏高原2017年牧草產(chǎn)量增加的原因。

    表4 三種情況青藏高原下牲畜空間年總暴露量(萬只)

    2.4 牲畜損失率和雪災(zāi)危險性指數(shù)之間關(guān)系

    青藏高原牧區(qū)牲畜雪災(zāi)損失定量研究關(guān)鍵是明晰雪災(zāi)危險性強弱和牧區(qū)牲畜死亡量多少之間的定量關(guān)系,通過災(zāi)害大典記錄的雪災(zāi)造成牲畜死亡量、雪災(zāi)發(fā)生時間、雪災(zāi)發(fā)生地牲畜暴露量、雪災(zāi)發(fā)生時氣溫、雪災(zāi)發(fā)生時積雪持續(xù)時長和最大積雪深度數(shù)據(jù),首先利用公式(1)計算出雪災(zāi)發(fā)生時的危險性指數(shù),然后計算出此雪災(zāi)造成牲畜死亡率,最后以雪災(zāi)危險性指數(shù)作為自變量,雪災(zāi)造成牲畜死亡率作為因變量,構(gòu)建青藏高原牧區(qū)牲畜在不同雪災(zāi)危險性強度下的損失率曲線。擬合結(jié)果如(圖7)所示。

    圖7 青藏高原牲畜損失率曲線

    曲線方程如公式(12)所示,該方程相關(guān)性為0.568 9,表明可以用該損失率曲線對不同雪災(zāi)危險指數(shù)下的損失率估算。

    y=0.000 04x1.891 9。

    (12)

    2.5 青藏高原牲畜受雪災(zāi)影響年均死亡量預(yù)估結(jié)果

    雪災(zāi)風(fēng)險預(yù)估所用數(shù)據(jù)來自于此次研究所得的雪災(zāi)危險時空分布數(shù)據(jù)、牲畜空間分布數(shù)據(jù)和牲畜草原雪災(zāi)損失率三類。其中牲畜空間分布數(shù)據(jù)則選用了“以草定畜”牲畜空間分布數(shù)據(jù),這是由于本次研究更多側(cè)重于自然狀況下牲畜受雪災(zāi)損失,但實際牲畜空間分布數(shù)量遠遠高于“以草定畜”牲畜數(shù)量,而牲畜超載量較大區(qū)域主要分布于農(nóng)牧交錯區(qū)域,經(jīng)過文獻查閱[16]以及現(xiàn)場調(diào)研,發(fā)現(xiàn)在農(nóng)牧交錯區(qū)域超載牲畜主要由于圈養(yǎng)牲畜量較大,往往雪災(zāi)對圈養(yǎng)牲畜影響并不大。因此2000年草原載畜量計算牲畜總量用來表示歷史時期青藏高原牲畜暴露量,2017年草原載畜量計算牲畜總量用作為未來青藏高原牲畜暴露量。

    預(yù)估結(jié)果顯示:空間分布來看,牲畜死亡量的空間分布與雪災(zāi)危險性指數(shù)的空間分布有著相似性,死亡量較大區(qū)域也分布在藏北高原、岡底斯山脈沿線、昆侖山脈西段沿線、祁連山脈沿線、三江源區(qū)域和橫斷山脈區(qū)域(圖8)。就牲畜死亡量來看,在未來RCPs情景下,牲畜受雪災(zāi)影響年死亡量呈現(xiàn)出減少趨勢,減少量較多地區(qū)是青藏高原中部區(qū)域以及三江源區(qū)域。相比較與歷史時期,RCP4.5情景下未來近期、未來遠期和RCP8.5情景下未來近期、未來遠期牲畜受雪災(zāi)影響,牲畜年均損失量分別減少了10.2%、10.3%、28.3%和45.5%。綜合來看,未來RCP4.5和RCP8.5情景下,發(fā)生雪災(zāi)危險性區(qū)域面積占比變化與青藏高原牲畜受雪災(zāi)影響死亡量變化呈現(xiàn)出較強的一致性,因此可以得出未來青藏高原牲畜受雪災(zāi)死亡減少主要受到了發(fā)生雪災(zāi)危險性區(qū)域面積占比減少影響(表3、圖8和表5)。

    2.6風(fēng)險預(yù)估結(jié)果不確定性分析

    (1)雪災(zāi)危險性不確定性分析。在此次研究中,運用模式數(shù)據(jù)識別了青藏高原區(qū)域內(nèi)雪災(zāi)危險性,而模式數(shù)據(jù)本身就存在有一定的不確定性[21],因此,對用于識別雪災(zāi)危險度的逐日雪深數(shù)據(jù)進行驗證,既ESM1-BGC模式數(shù)據(jù)是否符合實際情況?

    此次研究對雪災(zāi)危險性識別是基于NEX-GDDP數(shù)據(jù)集和中國逐日雪深模擬預(yù)估數(shù)據(jù)集,其中中國逐日雪深模擬預(yù)估數(shù)據(jù)集是以NEX-GDDP數(shù)據(jù)集為依托,利用人工神經(jīng)網(wǎng)絡(luò)模型模擬獲取,模型模擬過程中輸入變量為逐日降水數(shù)據(jù)、逐日最高氣溫和逐日最低氣溫,因此模擬雪深數(shù)據(jù)準確性很大程度上依賴于NEX-GDDP數(shù)據(jù)的準確性,為了研究青藏高原降雪量模擬的不確定性,統(tǒng)計了青藏高原區(qū)域內(nèi)CESM1-BGC模式模擬的降雪時期內(nèi)降水量年均值,然后與實際降水量進行比較(青藏高原降雪時段主要為冬春兩季)(圖9和圖10),結(jié)果顯示:青藏高原西部和西南部模式模擬降水量值要遠遠高于站點實測降水值。因此在此區(qū)域模擬雪深誤差較大原因是由于NEX-GDDP數(shù)據(jù)中CESM1-BGC模式下對降雪時段內(nèi)降水誤差較大。此次NEX-GDDP數(shù)據(jù)對該區(qū)域降水的模擬準確度不夠,導(dǎo)致該區(qū)域內(nèi)逐日雪深模擬精度較低。因此這種模式數(shù)據(jù)本身存在不確定性很大程度上影響本次風(fēng)險預(yù)估結(jié)果的不確定性。隨著CMIP6數(shù)據(jù)的發(fā)布,LUN等[22]通過對比CMIP6模式和CMIP5模式對青藏高原模擬能力發(fā)現(xiàn)CMIP6模式具有較好的降水模擬能力。其次CMIP6包含最新的共享社會經(jīng)濟路徑與典型濃度路徑組合情景(Shared Socioeconomic Pathways and the Representative Concentration Pathways, SSP)融入了社會經(jīng)濟發(fā)展的影響,將提供更可靠的氣候變化可能結(jié)果。然而,目前CMIP6數(shù)據(jù)的模式分辨率較小且不同模式分辨率差異較大,增加數(shù)據(jù)利用難度。因此,隨著CMIP6數(shù)據(jù)分辨率提高、統(tǒng)一,以后青藏高原牧區(qū)雪災(zāi)風(fēng)險研究中利用CMIP6數(shù)據(jù)將提高預(yù)估精度。

    (2)氣候變化可能對青藏高原牧區(qū)雪災(zāi)影響。青藏牧區(qū)雪災(zāi)形成過程是積雪過厚,牧草被大雪掩埋,靠牧草為生的家畜因吃不到草,凍餓而死。積雪的累積持續(xù)時間、最大降雪厚度以及積雪時期的氣溫將直接影響雪災(zāi)的嚴重程度,LEATHERS對積雪持續(xù)的時間進行了特征分析結(jié)果表明積雪持續(xù)的時間與降雪和溫度有著高度的相關(guān)性[24]。降雪量變化不大情況下溫度越高積雪累積持續(xù)時間越少。在全球變暖大背景下,IPCC第四次評估報告顯示1956—2005全球地表氣溫上升趨勢為0.13℃/10 a,而中國發(fā)布《第二次國家氣候被評估報告》顯示1951—2009年中國氣溫上升趨勢為0.23℃/10 a。相比于全球尺度和整個中國區(qū)域尺度,青藏高原作為氣候變化敏感區(qū)域,1960—2010年氣溫上升趨勢為0.37℃/10 a[25],明顯高于全球尺度和中國區(qū)域尺度氣溫上升,同時也有研究表明北半球積雪區(qū)內(nèi),青藏高原積雪范圍減少最為顯著[26],綜合以上幾點,青藏高原積雪持續(xù)時間和積雪范圍減少明顯,積雪期氣溫升高,因此區(qū)域內(nèi)雪災(zāi)發(fā)生的危險性降低。在未來,特別是RCP8.5情景下,發(fā)生雪災(zāi)危險性程度降低可能更加明顯。

    圖8 青藏高原牲畜受雪災(zāi)影響死亡量空間分布

    表5青藏高原牲畜受雪災(zāi)影響總的年均死亡量

    圖9 青藏高原區(qū)域內(nèi)利用站點實測的春冬兩季多年年均降水量[23]

    圖10 青藏高原區(qū)域內(nèi)CESM1-BGC模式模擬的降雪時期內(nèi)降水量多年年均值

    3 結(jié)論

    此次研究選取了歷史基準時段(1986—2005年)、未來近期(2016—2035年)和未來遠期(2046—2065年)三個時段,RCP4.5和RCP8.5兩種氣候情景,分別對青藏高原牧區(qū)雪災(zāi)危險、牲畜空間暴露量、牲畜的脆弱性以及青藏高原牲畜損失量進行研究,研究結(jié)果表明:

    (1)青藏高原區(qū)域內(nèi),所有模擬雪深數(shù)據(jù)模式中,CESM1-BGC模式下模擬的青藏高原雪深精度最高,文章中選擇CESM1-BGC模式模擬的逐日雪深數(shù)據(jù)用來識別雪深危險性??傮w來說,相比較于歷史時期,在未來RCP4.5和RCP8.5情景下,青藏高原發(fā)生雪災(zāi)危險性的區(qū)域縮小,RCP4.5情景下未來近期、未來遠期和RCP8.5情景下未來近期、未來遠期發(fā)生雪災(zāi)危險性的范圍較少6%、11%、6%和14%;但是發(fā)生區(qū)的雪災(zāi)危險強度有的時段可能增強;空間分布來看,危險性指數(shù)較高的區(qū)域主要分布在藏北高原、岡底斯山脈沿線、昆侖山脈西段沿線、祁連山脈沿線、三江源區(qū)域和橫斷山脈山脈區(qū)域。

    (2)利用“以草定畜”,計算結(jié)果顯示2017年青藏牧區(qū)草地載畜量相比較2000年載畜量增加11%,在未來有可能青藏高原載畜量將近一步增加。實際的牲畜與用“以草定畜”的牲畜比較來看,青藏高原區(qū)域內(nèi)實際牲畜保有量遠遠超出了青藏高原載畜量,從空間分布來看,超載區(qū)域主要分布于以拉薩和西寧為中心周邊區(qū)域以及橫斷山脈區(qū)域。

    (3)從空間分布來看牲畜死亡量較大區(qū)域與雪災(zāi)危險性指數(shù)較高的區(qū)域有著相似性,主要也分布在藏北高原、岡底斯山脈沿線、昆侖山脈西段沿線、祁連山脈沿線、三江源區(qū)域和橫斷山脈區(qū)域。相比較歷史時期,RCP4.5情景下未來近期、未來遠期和RCP8.5情景下未來近期、未來遠期牲畜受雪災(zāi)損失量分別減少了10.2%、10.3%、28.3%和45.5%,未來發(fā)生雪災(zāi)危險性面積減少是造成牲畜損失量減小主要原因,其中在RCP8.5情景下,損失最低。

    此次研究針對自然狀況以及傳統(tǒng)雪災(zāi)防御措施下的牲畜受災(zāi)狀況,并未考慮雪災(zāi)防御和災(zāi)后救援能力提升,經(jīng)過多年發(fā)展,牧區(qū)雪災(zāi)風(fēng)險預(yù)估已經(jīng)取得了豐碩成果[3-12,27,28],如:雪災(zāi)預(yù)警系統(tǒng)根據(jù)氣象數(shù)據(jù)對即將發(fā)生雪災(zāi)進行預(yù)報,雪災(zāi)救助區(qū)劃分更加合理,便于開展牧區(qū)雪災(zāi)救助。各種措施的實施,整體提升了對牧區(qū)雪災(zāi)的適應(yīng)能力。其次,青藏高原作為我國重要畜牧業(yè)基地,隨著牧草制度的發(fā)展,出現(xiàn)了牧草租借和流轉(zhuǎn),這將進一步實現(xiàn)畜牧業(yè)規(guī)?;a(chǎn),同時會有效推動牧區(qū)雪災(zāi)保險的普及。另外,隨著教育水平提升,牧業(yè)人口綜合素質(zhì)將得到全面提升,對雪災(zāi)的認識程度更加深刻??傮w來說,未來牧區(qū)規(guī)劃更加合理、雪災(zāi)預(yù)告更加精確、雪災(zāi)救助手段更加科學(xué)、救助設(shè)備更加先進以及牧業(yè)從業(yè)人員更加專業(yè),青藏高原牧區(qū)雪災(zāi)造成損失將進一步減少。

    猜你喜歡
    雪深雪災(zāi)牧區(qū)
    青海海西地區(qū)近41年雪災(zāi)時空分布特征
    試析藏北牧區(qū)的“鹽語”民俗①
    藏北牧區(qū)圈窩飼草種植及應(yīng)用前景
    一種基于頻率與相位算法的雷達雪深探測方法
    高原冬季雪深與重慶夏季降水的年際關(guān)系研究
    鐵路防災(zāi)雪深圖像采集的設(shè)計和實現(xiàn)
    歷史時期以來山西省雪災(zāi)特征與氣候變化
    高寒陰濕牧區(qū)羊痘病防治
    近51 a青藏高原雪災(zāi)時空分布特征
    2000-2010年新疆雪災(zāi)時空自相關(guān)分析
    免费观看精品视频网站| 中文字幕制服av| 不卡一级毛片| 欧美黄色片欧美黄色片| 欧美精品人与动牲交sv欧美| 欧美另类亚洲清纯唯美| 久久精品亚洲熟妇少妇任你| 一级a爱片免费观看的视频| 国产精品二区激情视频| 亚洲av欧美aⅴ国产| 很黄的视频免费| 亚洲欧美日韩高清在线视频| 精品国产超薄肉色丝袜足j| 国产高清激情床上av| 久久久国产精品麻豆| 国产av又大| 欧美大码av| 亚洲午夜精品一区,二区,三区| videos熟女内射| 五月开心婷婷网| 如日韩欧美国产精品一区二区三区| 在线播放国产精品三级| 搡老熟女国产l中国老女人| 曰老女人黄片| 国产乱人伦免费视频| 在线十欧美十亚洲十日本专区| 欧美最黄视频在线播放免费 | 天天操日日干夜夜撸| 中文字幕另类日韩欧美亚洲嫩草| 欧美人与性动交α欧美精品济南到| 正在播放国产对白刺激| 99国产精品一区二区三区| 中文字幕制服av| 亚洲精品美女久久av网站| 日韩大码丰满熟妇| 久久久久久久久免费视频了| 中文亚洲av片在线观看爽 | 亚洲在线自拍视频| 久久久久久人人人人人| 中文字幕人妻丝袜一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 王馨瑶露胸无遮挡在线观看| 乱人伦中国视频| 成年人免费黄色播放视频| 高清毛片免费观看视频网站 | 在线播放国产精品三级| 国产成人精品久久二区二区免费| 999久久久国产精品视频| 亚洲av成人av| 如日韩欧美国产精品一区二区三区| 丰满迷人的少妇在线观看| 一二三四在线观看免费中文在| 精品一区二区三卡| 亚洲七黄色美女视频| 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美网| 18禁观看日本| 香蕉丝袜av| 人人妻人人澡人人看| 精品国内亚洲2022精品成人 | 亚洲五月天丁香| 欧美中文综合在线视频| 亚洲三区欧美一区| 亚洲精品国产精品久久久不卡| 欧美精品一区二区免费开放| 啦啦啦在线免费观看视频4| 在线观看www视频免费| a级片在线免费高清观看视频| 麻豆av在线久日| 午夜成年电影在线免费观看| 欧美日韩精品网址| 99热只有精品国产| 99热国产这里只有精品6| 精品一区二区三区av网在线观看| 国产无遮挡羞羞视频在线观看| 99热只有精品国产| 亚洲色图 男人天堂 中文字幕| 我的亚洲天堂| 久久精品国产a三级三级三级| videosex国产| 狂野欧美激情性xxxx| 中文字幕精品免费在线观看视频| 亚洲美女黄片视频| 欧美日韩乱码在线| 搡老乐熟女国产| 变态另类成人亚洲欧美熟女 | 亚洲国产毛片av蜜桃av| 中亚洲国语对白在线视频| 久久久久精品国产欧美久久久| 一二三四在线观看免费中文在| 在线永久观看黄色视频| 久久香蕉激情| 在线观看舔阴道视频| 新久久久久国产一级毛片| 国产精品.久久久| 黄色片一级片一级黄色片| 日日爽夜夜爽网站| 丝瓜视频免费看黄片| 免费高清在线观看日韩| 久久久久久亚洲精品国产蜜桃av| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成a人片在线一区二区| netflix在线观看网站| 十八禁高潮呻吟视频| 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色 | 欧美激情高清一区二区三区| 国产精品欧美亚洲77777| 又紧又爽又黄一区二区| 亚洲,欧美精品.| 亚洲色图 男人天堂 中文字幕| 啦啦啦视频在线资源免费观看| 国产在线一区二区三区精| 国产精品av久久久久免费| 欧美一级毛片孕妇| 国产亚洲av高清不卡| 十八禁网站免费在线| 免费观看a级毛片全部| 精品高清国产在线一区| tube8黄色片| 亚洲黑人精品在线| 午夜免费鲁丝| 亚洲伊人色综图| 欧美日韩乱码在线| 久久国产精品人妻蜜桃| 国产深夜福利视频在线观看| 国产精品 国内视频| 9色porny在线观看| 国产深夜福利视频在线观看| 亚洲av熟女| 免费av中文字幕在线| 国产亚洲精品久久久久久毛片 | 亚洲专区字幕在线| 黄片小视频在线播放| 久9热在线精品视频| 在线观看免费高清a一片| 午夜免费成人在线视频| 电影成人av| 两个人看的免费小视频| 女人高潮潮喷娇喘18禁视频| 中亚洲国语对白在线视频| 人人妻人人添人人爽欧美一区卜| 色在线成人网| 欧美亚洲 丝袜 人妻 在线| 高潮久久久久久久久久久不卡| 精品电影一区二区在线| 国产精品久久久久久人妻精品电影| 亚洲一区高清亚洲精品| 91精品国产国语对白视频| 又黄又爽又免费观看的视频| 悠悠久久av| 韩国精品一区二区三区| 婷婷精品国产亚洲av在线 | 久久热在线av| 看黄色毛片网站| 国产视频一区二区在线看| 成人国产一区最新在线观看| 一区福利在线观看| 日韩三级视频一区二区三区| 亚洲精华国产精华精| 操美女的视频在线观看| 亚洲欧洲精品一区二区精品久久久| 色播在线永久视频| 黄色 视频免费看| 日韩欧美三级三区| av超薄肉色丝袜交足视频| 成熟少妇高潮喷水视频| 精品福利观看| 90打野战视频偷拍视频| 久久影院123| 欧美日韩亚洲高清精品| 69av精品久久久久久| 国产精品一区二区在线不卡| 国产日韩欧美亚洲二区| 黄色毛片三级朝国网站| 日韩成人在线观看一区二区三区| 中亚洲国语对白在线视频| 免费在线观看日本一区| 精品一品国产午夜福利视频| 欧美日韩福利视频一区二区| 国产精品98久久久久久宅男小说| 日本精品一区二区三区蜜桃| 亚洲精品乱久久久久久| 老熟妇仑乱视频hdxx| 免费观看人在逋| 国产精品欧美亚洲77777| 黄色 视频免费看| 免费女性裸体啪啪无遮挡网站| 久久午夜综合久久蜜桃| 色在线成人网| 欧美日韩一级在线毛片| 在线十欧美十亚洲十日本专区| 人妻一区二区av| 精品久久久久久电影网| 中文亚洲av片在线观看爽 | 九色亚洲精品在线播放| 国产精品一区二区在线观看99| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲五月色婷婷综合| 亚洲精品国产精品久久久不卡| 婷婷丁香在线五月| 精品电影一区二区在线| 桃红色精品国产亚洲av| 又黄又爽又免费观看的视频| 国产一区二区三区综合在线观看| 美女扒开内裤让男人捅视频| 午夜免费成人在线视频| 久9热在线精品视频| 免费日韩欧美在线观看| 久久精品人人爽人人爽视色| 成人18禁高潮啪啪吃奶动态图| 欧美+亚洲+日韩+国产| 欧美在线一区亚洲| 色播在线永久视频| 亚洲精品美女久久av网站| 亚洲综合色网址| 不卡一级毛片| 久久精品成人免费网站| 在线十欧美十亚洲十日本专区| 久久久久精品人妻al黑| 午夜久久久在线观看| 欧美 日韩 精品 国产| 黄频高清免费视频| 亚洲国产欧美网| 亚洲精品一二三| 日本vs欧美在线观看视频| 亚洲欧美激情在线| 69av精品久久久久久| 国产1区2区3区精品| av在线播放免费不卡| 一区在线观看完整版| 国产精品国产高清国产av | 高清欧美精品videossex| 亚洲精品一二三| 女人被狂操c到高潮| 亚洲熟女毛片儿| 欧美日韩亚洲综合一区二区三区_| 国产aⅴ精品一区二区三区波| 日本五十路高清| 大香蕉久久成人网| 色老头精品视频在线观看| 国产亚洲欧美在线一区二区| 高清在线国产一区| 久久狼人影院| 夜夜夜夜夜久久久久| 亚洲熟女毛片儿| 日本a在线网址| 看免费av毛片| 人妻一区二区av| 成熟少妇高潮喷水视频| 免费黄频网站在线观看国产| 久9热在线精品视频| 午夜激情av网站| 国产精品自产拍在线观看55亚洲 | 大陆偷拍与自拍| 黄片播放在线免费| 国产精品久久久人人做人人爽| 性少妇av在线| 好男人电影高清在线观看| av超薄肉色丝袜交足视频| 韩国av一区二区三区四区| 亚洲,欧美精品.| 精品午夜福利视频在线观看一区| 精品亚洲成a人片在线观看| 一区二区三区精品91| 亚洲国产精品合色在线| 国产有黄有色有爽视频| 亚洲午夜理论影院| 日韩三级视频一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 少妇被粗大的猛进出69影院| 亚洲avbb在线观看| 日本欧美视频一区| 免费不卡黄色视频| www.自偷自拍.com| 18禁黄网站禁片午夜丰满| 国产一区在线观看成人免费| 久久中文看片网| 国产精品香港三级国产av潘金莲| 视频在线观看一区二区三区| 免费一级毛片在线播放高清视频 | e午夜精品久久久久久久| 老司机亚洲免费影院| 欧美精品高潮呻吟av久久| 交换朋友夫妻互换小说| 国产精品一区二区在线不卡| 高清av免费在线| 亚洲专区国产一区二区| 岛国毛片在线播放| 久久久久国产一级毛片高清牌| 亚洲五月色婷婷综合| 久久国产乱子伦精品免费另类| 少妇裸体淫交视频免费看高清 | 亚洲中文字幕日韩| 亚洲av美国av| 亚洲人成电影免费在线| 精品久久久精品久久久| 欧美黑人精品巨大| 在线观看免费日韩欧美大片| 丝袜美足系列| 精品国产一区二区三区四区第35| 视频在线观看一区二区三区| 免费观看精品视频网站| 91成人精品电影| 中国美女看黄片| 99精品在免费线老司机午夜| 亚洲伊人色综图| 老鸭窝网址在线观看| 淫妇啪啪啪对白视频| 午夜影院日韩av| 在线观看午夜福利视频| 色播在线永久视频| 午夜福利一区二区在线看| 十八禁人妻一区二区| 色尼玛亚洲综合影院| 很黄的视频免费| 亚洲中文日韩欧美视频| 国产日韩欧美亚洲二区| 中文字幕人妻丝袜制服| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久久久精品古装| 91成年电影在线观看| 美女扒开内裤让男人捅视频| 一级毛片女人18水好多| 国产亚洲精品第一综合不卡| 欧美精品人与动牲交sv欧美| 窝窝影院91人妻| 天天躁狠狠躁夜夜躁狠狠躁| 视频在线观看一区二区三区| 国产1区2区3区精品| 国产真人三级小视频在线观看| 十八禁人妻一区二区| 国产成人免费无遮挡视频| 亚洲第一青青草原| 搡老乐熟女国产| av视频免费观看在线观看| 欧美另类亚洲清纯唯美| av在线播放免费不卡| 中出人妻视频一区二区| 欧美av亚洲av综合av国产av| 欧美一级毛片孕妇| 成年动漫av网址| 久久久久视频综合| 日韩成人在线观看一区二区三区| 中出人妻视频一区二区| 久久狼人影院| 精品无人区乱码1区二区| 女性被躁到高潮视频| 午夜免费鲁丝| 91字幕亚洲| 亚洲av欧美aⅴ国产| av国产精品久久久久影院| 极品少妇高潮喷水抽搐| 久久人妻熟女aⅴ| 高清在线国产一区| 韩国精品一区二区三区| 老鸭窝网址在线观看| 窝窝影院91人妻| 久久精品熟女亚洲av麻豆精品| 久久性视频一级片| 咕卡用的链子| 亚洲成人国产一区在线观看| 操美女的视频在线观看| 麻豆乱淫一区二区| 久久人妻福利社区极品人妻图片| 制服人妻中文乱码| 亚洲成人免费av在线播放| 在线观看免费日韩欧美大片| 亚洲专区字幕在线| www.999成人在线观看| 亚洲精品国产一区二区精华液| 中文字幕制服av| 精品福利永久在线观看| 老司机亚洲免费影院| 亚洲七黄色美女视频| 色婷婷久久久亚洲欧美| 久久亚洲真实| 侵犯人妻中文字幕一二三四区| 老鸭窝网址在线观看| 亚洲欧美精品综合一区二区三区| 一边摸一边抽搐一进一小说 | 老司机影院毛片| 91在线观看av| 欧美日韩乱码在线| 女性被躁到高潮视频| 亚洲aⅴ乱码一区二区在线播放 | 国产又色又爽无遮挡免费看| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美一区二区三区在线观看 | 精品一区二区三卡| 男女高潮啪啪啪动态图| 日本精品一区二区三区蜜桃| 国产成人免费无遮挡视频| 色94色欧美一区二区| 91在线观看av| 国产成人精品久久二区二区91| 欧美日韩av久久| 中亚洲国语对白在线视频| 欧美成人午夜精品| 免费观看a级毛片全部| 亚洲专区中文字幕在线| 欧美黄色淫秽网站| 亚洲国产看品久久| 在线播放国产精品三级| 精品人妻熟女毛片av久久网站| 怎么达到女性高潮| 老司机影院毛片| 欧美乱色亚洲激情| 亚洲五月天丁香| 一边摸一边做爽爽视频免费| 国产欧美日韩一区二区精品| 丝瓜视频免费看黄片| 国产区一区二久久| 日韩 欧美 亚洲 中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 一本一本久久a久久精品综合妖精| 久久人人爽av亚洲精品天堂| videosex国产| 两性午夜刺激爽爽歪歪视频在线观看 | 天堂中文最新版在线下载| 咕卡用的链子| 久久久久精品人妻al黑| 高清在线国产一区| 多毛熟女@视频| 亚洲专区国产一区二区| 黄频高清免费视频| 亚洲精品av麻豆狂野| 一级a爱片免费观看的视频| 成年人黄色毛片网站| 一边摸一边抽搐一进一出视频| 久久人妻熟女aⅴ| 王馨瑶露胸无遮挡在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 王馨瑶露胸无遮挡在线观看| 亚洲中文av在线| 亚洲 欧美一区二区三区| 午夜免费观看网址| 熟女少妇亚洲综合色aaa.| 巨乳人妻的诱惑在线观看| 亚洲第一欧美日韩一区二区三区| 国产日韩欧美亚洲二区| 一级作爱视频免费观看| 老司机靠b影院| 香蕉国产在线看| 国产精品98久久久久久宅男小说| 亚洲精品中文字幕一二三四区| 久久热在线av| 欧美老熟妇乱子伦牲交| 男女高潮啪啪啪动态图| 欧美在线黄色| 在线av久久热| 国产精品99久久99久久久不卡| 国产精品一区二区在线观看99| а√天堂www在线а√下载 | 少妇粗大呻吟视频| 视频区欧美日本亚洲| 欧美国产精品一级二级三级| 自线自在国产av| 制服诱惑二区| 亚洲五月色婷婷综合| 最近最新中文字幕大全免费视频| 人妻丰满熟妇av一区二区三区 | 嫁个100分男人电影在线观看| 又黄又粗又硬又大视频| 俄罗斯特黄特色一大片| 韩国av一区二区三区四区| а√天堂www在线а√下载 | aaaaa片日本免费| 国产一区在线观看成人免费| 国产在线精品亚洲第一网站| 亚洲人成77777在线视频| 欧美人与性动交α欧美软件| 亚洲精品自拍成人| 好男人电影高清在线观看| 大码成人一级视频| 1024香蕉在线观看| 精品人妻熟女毛片av久久网站| 99久久99久久久精品蜜桃| 99国产综合亚洲精品| 免费日韩欧美在线观看| 欧美乱妇无乱码| 欧美+亚洲+日韩+国产| 老司机午夜福利在线观看视频| 日韩 欧美 亚洲 中文字幕| 精品第一国产精品| 国产av精品麻豆| 国产精品国产高清国产av | 精品亚洲成a人片在线观看| 日韩欧美在线二视频 | 老司机深夜福利视频在线观看| 亚洲一区二区三区欧美精品| 好男人电影高清在线观看| 色综合欧美亚洲国产小说| 久久天堂一区二区三区四区| 男男h啪啪无遮挡| 亚洲av日韩在线播放| 麻豆乱淫一区二区| 日韩一卡2卡3卡4卡2021年| 精品一品国产午夜福利视频| 国产激情欧美一区二区| 满18在线观看网站| 精品国产超薄肉色丝袜足j| 日本撒尿小便嘘嘘汇集6| 久久精品国产亚洲av香蕉五月 | 午夜亚洲福利在线播放| 久久久精品区二区三区| 十分钟在线观看高清视频www| 少妇被粗大的猛进出69影院| 亚洲中文日韩欧美视频| 亚洲五月天丁香| 999久久久国产精品视频| 人妻久久中文字幕网| 在线免费观看的www视频| 久久精品国产综合久久久| 久久天躁狠狠躁夜夜2o2o| 身体一侧抽搐| 18禁裸乳无遮挡动漫免费视频| 91麻豆精品激情在线观看国产 | 极品少妇高潮喷水抽搐| 国产高清国产精品国产三级| 国产精品自产拍在线观看55亚洲 | 亚洲五月婷婷丁香| 国产亚洲精品第一综合不卡| 精品第一国产精品| 91av网站免费观看| 欧美日韩黄片免| 十分钟在线观看高清视频www| avwww免费| 怎么达到女性高潮| 最近最新中文字幕大全免费视频| 悠悠久久av| 老司机福利观看| 99在线人妻在线中文字幕 | 如日韩欧美国产精品一区二区三区| 村上凉子中文字幕在线| av天堂久久9| 黑丝袜美女国产一区| 亚洲精品一卡2卡三卡4卡5卡| 99久久国产精品久久久| 亚洲精品久久成人aⅴ小说| 看片在线看免费视频| 亚洲av日韩在线播放| 成人av一区二区三区在线看| 亚洲精品国产一区二区精华液| 亚洲精品国产精品久久久不卡| 曰老女人黄片| 午夜老司机福利片| 一进一出好大好爽视频| 天天躁日日躁夜夜躁夜夜| 老鸭窝网址在线观看| 丝瓜视频免费看黄片| 正在播放国产对白刺激| 亚洲九九香蕉| 妹子高潮喷水视频| 精品久久蜜臀av无| 91精品三级在线观看| 757午夜福利合集在线观看| 午夜精品久久久久久毛片777| 日韩视频一区二区在线观看| 在线观看舔阴道视频| 国产成人精品久久二区二区91| 成人18禁高潮啪啪吃奶动态图| 一级a爱视频在线免费观看| 精品免费久久久久久久清纯 | www.999成人在线观看| 一级毛片高清免费大全| 午夜亚洲福利在线播放| 一级黄色大片毛片| 久久九九热精品免费| 国产精品一区二区精品视频观看| 精品免费久久久久久久清纯 | 男女下面插进去视频免费观看| 看片在线看免费视频| 国产精品九九99| 好看av亚洲va欧美ⅴa在| 亚洲七黄色美女视频| 国产视频一区二区在线看| 日韩有码中文字幕| 精品久久久久久,| 多毛熟女@视频| 天堂√8在线中文| 91九色精品人成在线观看| 视频在线观看一区二区三区| 国产欧美日韩一区二区三区在线| 久久久国产成人精品二区 | 亚洲精品粉嫩美女一区| 午夜日韩欧美国产| 99国产精品一区二区蜜桃av | 在线观看日韩欧美| 亚洲成人免费av在线播放| 黄色丝袜av网址大全| 亚洲精品美女久久久久99蜜臀| 黄色毛片三级朝国网站| 久久久久久久精品吃奶| 建设人人有责人人尽责人人享有的| 一区在线观看完整版| 丝瓜视频免费看黄片| 高清av免费在线| 日韩三级视频一区二区三区| 亚洲自偷自拍图片 自拍| a在线观看视频网站| 午夜精品国产一区二区电影| 少妇的丰满在线观看| 国产91精品成人一区二区三区| 一区二区三区国产精品乱码| 久久人妻熟女aⅴ| 黑人巨大精品欧美一区二区mp4| 91国产中文字幕| 色在线成人网| 黑人欧美特级aaaaaa片| 99国产综合亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 国产极品粉嫩免费观看在线| av在线播放免费不卡| 两个人看的免费小视频|