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

    綜合干旱指數(shù)構(gòu)建及其在不同草原類型中的應(yīng)用

    2021-11-26 06:04:32趙水霞王文君吳英杰
    關(guān)鍵詞:旱情葉綠素降水

    趙水霞,王文君,吳英杰,李 瑋,全 強(qiáng)

    綜合干旱指數(shù)構(gòu)建及其在不同草原類型中的應(yīng)用

    趙水霞,王文君※,吳英杰,李 瑋,全 強(qiáng)

    (1. 中國(guó)水利水電科學(xué)研究院內(nèi)蒙古陰山北麓荒漠草原生態(tài)水文野外科學(xué)觀測(cè)研究站,北京 100038;2.水利部牧區(qū)水利科學(xué)研究所,呼和浩特 010020)

    干旱指數(shù)是評(píng)估旱情等級(jí)、制定防旱減災(zāi)對(duì)策的重要指標(biāo)。該研究從干旱發(fā)生的物理機(jī)制出發(fā),基于CRITIC(Criteria Importance Through Inter-criteria Correlation)客觀賦權(quán)法,綜合考慮氣象標(biāo)準(zhǔn)化降水蒸散指數(shù)(Standard Precipitation Evapotranspiration Index,SPEI)、植被供水指數(shù)(Vegetation Supply Water Index,VSWI)和日光誘導(dǎo)葉綠素?zé)晒庵笖?shù)(Solar-induced Chlorophyll Fluorescence,SIF)在干旱發(fā)展過(guò)程中的作用機(jī)制,構(gòu)建融合多源信息的綜合干旱指數(shù),探究其在不同草原類型干旱中的適用性和敏感性,進(jìn)而分析干旱的時(shí)空分布及演化特征。結(jié)果表明:綜合干旱指數(shù)既能敏銳捕捉氣象干旱的早期開(kāi)始,亦能從水文和牧業(yè)干旱角度綜合反映旱情的持續(xù)過(guò)程和時(shí)間,相比單一干旱指數(shù),具有綜合表征多尺度干旱演化特征的優(yōu)勢(shì);2007—2018年,內(nèi)蒙古錫林郭勒草原干旱頻率呈波動(dòng)中下降趨勢(shì),干旱強(qiáng)度主要以輕中旱為主,重旱及以上干旱集中在春秋兩季;空間分布上從西北向東南逐漸降低,呈現(xiàn)荒漠草原>典型草原>草甸草原,典型草原中北部旱情嚴(yán)重的特征。該研究對(duì)蒙古高原不同草原類型干旱事件精準(zhǔn)監(jiān)測(cè)及旱災(zāi)防御具有重要意義。

    干旱;降水;綜合干旱指數(shù);葉綠素?zé)晒庵笖?shù);適用性;錫林郭勒草原

    0 引 言

    干旱作為全球最為突出的自然災(zāi)害之一,具有持續(xù)時(shí)間長(zhǎng)、影響范圍廣及旱災(zāi)損失大等特點(diǎn),全球氣候變暖趨勢(shì)下,干旱災(zāi)害風(fēng)險(xiǎn)日益嚴(yán)峻,嚴(yán)重影響和制約著生態(tài)環(huán)境安全、農(nóng)牧業(yè)高質(zhì)量發(fā)展及人類生活水平提高,使得實(shí)現(xiàn)干旱的精準(zhǔn)化監(jiān)測(cè),變被動(dòng)減災(zāi)為主動(dòng)防災(zāi)是國(guó)內(nèi)外亟待解決的關(guān)鍵科學(xué)問(wèn)題,具有重要的現(xiàn)實(shí)意義[1-3]。錫林郭勒草原位于蒙古高原東南部,是內(nèi)蒙古天然草原的主體,在植被類型上具有一定的完整性和典型性,而頻發(fā)的干旱事件嚴(yán)重制約了當(dāng)?shù)厣鐣?huì)經(jīng)濟(jì)的可持續(xù)發(fā)展,使得實(shí)現(xiàn)旱情分布及演化趨勢(shì)的精準(zhǔn)化監(jiān)測(cè)成為“十四五”期間亟待解決的關(guān)鍵科學(xué)問(wèn)題之一[4-5]。

    適用的干旱指數(shù)是確保旱情精準(zhǔn)監(jiān)測(cè)的基礎(chǔ),受干旱復(fù)雜性及自然水循環(huán)過(guò)程相互關(guān)聯(lián)的影響,單一干旱指數(shù)往往難以準(zhǔn)確反映干旱的多類型及多尺度特征,使得構(gòu)建基于指標(biāo)-影響關(guān)聯(lián)分析的多源信息綜合干旱指數(shù)成為目前研究的難點(diǎn)問(wèn)題及熱點(diǎn)趨勢(shì)[4, 6-7]。綜合干旱指數(shù)的構(gòu)建分為考慮氣象、水文和農(nóng)業(yè)致旱因子的綜合干旱指數(shù)和多種單一干旱指數(shù)組合的綜合干旱指數(shù),構(gòu)建方法多為權(quán)重組合、多變量聯(lián)合分布和機(jī)器學(xué)習(xí),單一干旱指數(shù)較少時(shí)常選用結(jié)構(gòu)簡(jiǎn)單且易于理解的權(quán)重組合方法,如CRITIC(Criteria Importance Through Inter-criteria Correlation客觀賦權(quán)法[8]?;谥潞狄蜃拥木C合干旱指數(shù)多采用線性回歸及數(shù)據(jù)挖掘等手段,綜合考慮氣象降水盈虧、土壤水分供給及植被生長(zhǎng)條件的相互作用,進(jìn)而構(gòu)建適用于不同區(qū)域下墊面的綜合干旱監(jiān)測(cè)模型[9-10];以草地生態(tài)系統(tǒng)為例,降水與蒸散失衡造成土壤水分虧缺是干旱發(fā)生的主要誘因,植被受干旱水分脅迫影響,光合作用和物質(zhì)代謝受阻,進(jìn)而降低葉片生長(zhǎng)速率和植被蒸騰效應(yīng),導(dǎo)致植被生物量銳減、生產(chǎn)力下降及草畜不平衡。傳統(tǒng)的地面數(shù)據(jù)點(diǎn)狀監(jiān)測(cè),難以準(zhǔn)確刻畫(huà)大尺度旱情空間分布,隨著遙感監(jiān)測(cè)技術(shù)的不斷發(fā)展,越來(lái)越多的遙感干旱指數(shù)被用于綜合干旱指數(shù)的構(gòu)建及旱情的精準(zhǔn)監(jiān)測(cè),如歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)、溫度植被干旱指數(shù)(Temperature Vegetation Drought Index,TVDI)及土壤濕度監(jiān)測(cè)指數(shù)(Soil Moisture Monitoring Index,SMMI)等[4,11-13]。近年來(lái),新型日光誘導(dǎo)葉綠素?zé)晒猓⊿olar-induced Chlorophyll Fluorescence,SIF)遙感指數(shù)被發(fā)現(xiàn)對(duì)水分脅迫非常敏感,在碳氮循環(huán)、初級(jí)生產(chǎn)力估算及植被脅迫等領(lǐng)域被廣泛應(yīng)用,其原理是植被在吸收太陽(yáng)能進(jìn)行光合作用時(shí)會(huì)釋放出一種長(zhǎng)波信號(hào),該信號(hào)為自然狀態(tài)下的葉綠素?zé)晒猓渥鳛楣夂献饔玫母碑a(chǎn)物對(duì)水分脅迫的敏感性顯著高于主要通過(guò)反射率估算的植被指數(shù)(歸一化植被指數(shù)、增強(qiáng)型植被指數(shù))[14-18],在有效監(jiān)測(cè)植被瞬時(shí)光合作用、受脅迫狀態(tài)及物候監(jiān)測(cè)等方面具有明顯優(yōu)勢(shì),一定程度上緩減了作物對(duì)干旱脅迫的滯后效應(yīng),提高了干旱監(jiān)測(cè)的準(zhǔn)確度。

    融合多源地面及遙感信息的綜合干旱指數(shù)能較準(zhǔn)確及時(shí)地監(jiān)測(cè)干旱過(guò)程、客觀全面地評(píng)估干旱程度,但從綜合干旱指數(shù)的應(yīng)用與推廣技術(shù)上來(lái)看,綜合干旱指數(shù)的區(qū)域適用性和差異性問(wèn)題仍然存在[8,19-21],有必要從干旱機(jī)制及致旱因子相互耦合作用出發(fā),進(jìn)一步構(gòu)建適用于不同下墊面及不同氣候條件的有效綜合干旱指數(shù)及旱情等級(jí)評(píng)估標(biāo)準(zhǔn)。本研究采用氣象學(xué)、生態(tài)學(xué)及土壤學(xué)等多學(xué)科交叉理論,在考慮表征氣象干旱的標(biāo)準(zhǔn)化降水蒸散指數(shù)(Standard Precipitation Evapotranspiration Index,SPEI)和土壤含水率的植被供水指數(shù)(Vegetation Supply Water Index,VSWI)基礎(chǔ)上,嘗試引入植被葉綠素?zé)晒庵笖?shù)(Solar-induced Chlorophyll Fluorescence,SIF),利用CRITIC 客觀賦權(quán)及回歸分析等方法,構(gòu)建適用于錫林郭勒不同類型草原旱情監(jiān)測(cè)的綜合干旱指數(shù),以期提高綜合干旱指數(shù)的針對(duì)性,實(shí)現(xiàn)旱情的精準(zhǔn)化監(jiān)測(cè)與評(píng)估,為氣候變化下牧區(qū)防旱減災(zāi)及可持續(xù)發(fā)展提供科學(xué)依據(jù)。

    1 研究區(qū)概況與數(shù)據(jù)來(lái)源

    1.1 研究區(qū)概況

    錫林郭勒草原作為內(nèi)蒙古北方地區(qū)重要的生態(tài)屏障及畜牧業(yè)生產(chǎn)基地,受氣候變化和持續(xù)性干旱的影響,區(qū)域防旱減災(zāi)任務(wù)面臨嚴(yán)峻挑戰(zhàn)。研究區(qū)位于內(nèi)蒙古中東部,地勢(shì)呈現(xiàn)東南低山丘陵、西北地形平坦的特點(diǎn),如圖1所示。區(qū)域氣候?qū)儆跍貛Т箨懶约撅L(fēng)氣候,具有明顯的季節(jié)性差異和地域差異,降水量由東南向西北逐漸減少,多年平均降水量為136~376 mm,年均氣溫為0.88~5.21 ℃;受降水分布格局及水分差異的影響,錫林郭勒草原類型從東北向西南涵蓋了東南部草甸草原(多倫、西烏珠穆沁旗)、中部典型草原(東烏珠穆沁旗、那仁寶力格、錫林浩特和阿巴嘎旗)和西部荒漠草原(二連浩特、蘇尼特左旗)[22-23],為綜合干旱指數(shù)的研究提供了較全面的區(qū)域條件。

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

    研究區(qū)氣象數(shù)據(jù)和土壤相對(duì)濕度數(shù)據(jù)資料來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://data.cma.cn)。選取錫林郭勒草原時(shí)間序列較長(zhǎng)且數(shù)據(jù)較為齊全的8個(gè)氣象站點(diǎn)資料,通過(guò)數(shù)據(jù)整編獲取2007—2018年各站點(diǎn)日均降水、溫度、日照時(shí)數(shù)和相對(duì)濕度等數(shù)據(jù),用于標(biāo)準(zhǔn)化降水蒸散指數(shù)的計(jì)算和綜合干旱指數(shù)的構(gòu)建。

    2 研究方法

    2.1 單一干旱指數(shù)計(jì)算

    1)標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)

    SPEI是一種綜合考慮降水和蒸散的標(biāo)準(zhǔn)化干旱指數(shù),其絕對(duì)值的大小反映了干旱或濕潤(rùn)的嚴(yán)重程度,假設(shè)降水與蒸散的差值服從Log-logistic概率分布,對(duì)其進(jìn)行標(biāo)準(zhǔn)化處理,得到如下方程[24]:

    式中常數(shù)項(xiàng)0、1、2、1、2、3分別為2.515 517、0.802 853、0.010 328、1.432 788、0.189 269、0.001 308;為累積概率加權(quán)矩;累積概率為對(duì)降水和潛在蒸散量差值數(shù)據(jù)系列的正態(tài)化處理結(jié)果,其中,潛在蒸散量的計(jì)算采用彭曼公式[22]。

    通過(guò)計(jì)算各站點(diǎn)日值SPEI指數(shù),進(jìn)一步均值化得到研究區(qū)時(shí)間尺度統(tǒng)一的每16 d一次SPEI數(shù)值。

    2)植被供水指數(shù)(VSWI)

    植被供水指數(shù)(VSWI)以地面植被指數(shù)和地表溫度為指標(biāo),用于間接反映土壤含水率及作物受干旱影響狀況[25]。當(dāng)植被受到干旱脅迫時(shí),葉面氣孔為減少水分損失產(chǎn)生一定閉合,進(jìn)而導(dǎo)致冠層溫度升高,植被指數(shù)降低,如下所示:

    式中T為地表溫度,℃;NDVI為歸一化植被指數(shù);分別采用中分辨率成像光譜儀MODIS傳感器通道31(T)、通道1和通道2(NDVI)數(shù)據(jù)進(jìn)行反演。

    3)日光誘導(dǎo)葉綠素?zé)晒庵笖?shù)(SIF)

    植被在吸收太陽(yáng)輻射通量后,在波長(zhǎng)600~800 nm范圍內(nèi)發(fā)射出長(zhǎng)波信號(hào),稱為日光誘導(dǎo)葉綠素?zé)晒?,其?huì)受水分虧缺影響呈現(xiàn)不同程度下降,可用于監(jiān)測(cè)植被的生理狀態(tài)及水分脅迫狀況[14-15]。目前可用于探測(cè)熒光的衛(wèi)星主要有GOSAT(日本)、MetOp(GOME-2)、OCO-2(美國(guó))、TANSAT(中國(guó))和Sentinel-5P(TROPOMI)等[26-27],Li等[28]采用Cubist模型樹(shù)算法,基于OCO-2 SIF、MODIS EVI及氣象數(shù)據(jù),構(gòu)建了時(shí)間分辨率為8 d,空間分辨率為0.05°的葉綠素?zé)晒猱a(chǎn)品(Global‘OCO-2’ SIF,GOSIF),該產(chǎn)品相比OCO-2 SIF和MODIS GPP而言,可更加準(zhǔn)確地估算植被初級(jí)生產(chǎn)力,反映區(qū)域氣候變化特征[29]。

    本研究采用2007—2018年錫林郭勒草原的GOSIF數(shù)據(jù)(https://globalecology.unh.edu/),在ARCGIS數(shù)據(jù)空間分析系統(tǒng)中對(duì)不同類型草原的葉綠素?zé)晒鈹?shù)據(jù)進(jìn)行提取和插值分析。

    2.2 綜合干旱指數(shù)(SSV)構(gòu)建方法

    采用CRITIC 賦權(quán)法進(jìn)行綜合干旱指數(shù)(SSV)的構(gòu)建。CRITIC賦權(quán)法是一種客觀賦權(quán)方法,將SPEI、VSWI和SIF作為綜合干旱指數(shù)的輸入指標(biāo),利用指標(biāo)間的沖突性和信息量大小確定其權(quán)重值。兩個(gè)指標(biāo)間的相關(guān)系數(shù)越高,代表沖突性越??;單個(gè)指標(biāo)的標(biāo)準(zhǔn)離差越大,其包含的信息量越大。該方法結(jié)合了相關(guān)性權(quán)重和信息量權(quán)重,在評(píng)價(jià)因子權(quán)重中具有顯著的優(yōu)越性[12]。

    采用離差標(biāo)準(zhǔn)化法對(duì)數(shù)據(jù)進(jìn)行歸一化處理,保證基本度量單位統(tǒng)一且服從0~1之間的統(tǒng)計(jì)概率分布。

    式中XD分別為數(shù)據(jù)歸一化處理前后的值;min和max分別為樣本的最小值和最大值。

    指標(biāo)間的沖突性T和單個(gè)指標(biāo)的信息量C表示為

    式中r為指標(biāo)和指標(biāo)之間的相關(guān)系數(shù);δ為第個(gè)指標(biāo)的標(biāo)準(zhǔn)差,為一個(gè)評(píng)價(jià)指標(biāo)的評(píng)價(jià)數(shù)量。C越大,則該指標(biāo)的重要性越大。設(shè)W為權(quán)重系數(shù),則

    式中為綜合干旱指標(biāo)中單一指標(biāo)數(shù)量。

    綜上所述,綜合干旱指數(shù)(SSV)可表示為

    式中D分別為SPEI、SIF和VSWI的歸一化數(shù)值。

    3 結(jié)果與分析

    3.1 基于綜合干旱指數(shù)干旱等級(jí)劃分

    基于數(shù)據(jù)正態(tài)標(biāo)準(zhǔn)化處理得到的綜合干旱指數(shù)能綜合表征降水與蒸散的水量平衡過(guò)程、土壤水分脅迫強(qiáng)度及植被敏感度狀態(tài)。針對(duì)不同草原下墊面植被類型、植被蓋度差異及其對(duì)降水和溫度的敏感性差異,建立適用于不同季節(jié)和不同草原類型的綜合干旱指數(shù)旱情等級(jí)評(píng)價(jià)標(biāo)準(zhǔn),對(duì)旱情發(fā)生發(fā)展過(guò)程及演化趨勢(shì)的精準(zhǔn)定量化監(jiān)測(cè)具有重要意義。根據(jù)CRITIC客觀賦權(quán)法各致旱因子指標(biāo)間的沖突性和對(duì)比強(qiáng)度,不同草原類型各指標(biāo)權(quán)重系數(shù)如表1所示。

    錫林郭勒不同草原類型降水和土壤相對(duì)濕度是表征氣象干旱和農(nóng)業(yè)干旱的重要指標(biāo),且土壤相對(duì)濕度信息是采取應(yīng)對(duì)干旱措施的重要前提條件。本研究參照土壤相對(duì)濕度干旱等級(jí)劃分標(biāo)準(zhǔn)方法[30],通過(guò)建立與綜合干旱指數(shù)的相關(guān)性模型,獲得基于綜合干旱指數(shù)的植被生長(zhǎng)期干旱等級(jí)劃分標(biāo)準(zhǔn),如表2所示。

    表1 CRITIC客觀賦權(quán)法各指標(biāo)權(quán)重系數(shù)

    注:SIF為日光誘導(dǎo)葉綠素?zé)晒庵笖?shù);SPEI為標(biāo)準(zhǔn)化降水蒸散指數(shù);VSWI為植被供水指數(shù)。

    Note: SIF is solar-induced chlorophyll fluorescence; SPEI is standard precipitation evapotranspiration index; VSWI is vegetation supply water index.

    表2 綜合干旱指數(shù)SSV等級(jí)劃分

    3.2 基于綜合干旱指數(shù)的干旱時(shí)空分布特征分析

    3.2.1 錫林郭勒草原干旱的年際變化特征

    圖2為2007—2018年基于綜合干旱指數(shù)的錫林郭勒草原不同程度干旱發(fā)生頻率,可以發(fā)現(xiàn):1)2007年以來(lái),研究區(qū)不同程度的干旱每年均有發(fā)生,主要以輕旱和中旱為主;2)2009—2012年,干旱發(fā)生頻率呈一定下降趨勢(shì),旱情有所緩減,但極端干旱事件增多,重旱和特旱發(fā)生頻次增加;3)2012年以后,干旱頻率顯著增加且旱情程度有所加重,特旱頻率最高年份出現(xiàn)在2016年,當(dāng)年干旱發(fā)生頻率達(dá)到82.3%。整體而言,21世紀(jì),錫林郭勒盟草原旱情呈波動(dòng)中略有下降特征。

    經(jīng)查閱《中國(guó)水旱災(zāi)害公報(bào)》《中國(guó)氣象災(zāi)害年鑒》《內(nèi)蒙古水旱災(zāi)害》及錫林郭勒盟防汛辦旱情統(tǒng)計(jì)數(shù)據(jù),錫林郭勒草原2007—2018年旱情主要呈現(xiàn)以下特征:2007年春夏連旱,旱情等級(jí)為重旱;2012年和2013年多地為無(wú)旱和輕旱;2014年8月上旬,錫林郭勒草原大部分地區(qū)出現(xiàn)輕度干旱,且西北部干旱嚴(yán)重;2014—2016年全區(qū)發(fā)生連續(xù)3a嚴(yán)重干旱,錫林郭勒草原大部分地區(qū)承受嚴(yán)重災(zāi)損,受2014和2015年持續(xù)性干旱影響,2016年錫林郭勒草原旱情嚴(yán)重,受旱面積達(dá)到1 064.24萬(wàn)hm2,直接經(jīng)濟(jì)損失8.27億元;2017年4—5 月阿巴嘎旗、蘇尼特左旗及東烏珠穆沁旗等地發(fā)生干旱。上述旱情結(jié)果與基于綜合干旱指數(shù)評(píng)判的干旱頻率和干旱程度相契合,一定程度上揭示了綜合干旱指數(shù)的可靠性。

    錫林郭勒草原是內(nèi)蒙古自治區(qū)主要的畜牧業(yè)和飼草料生產(chǎn)基地,將內(nèi)蒙古自治區(qū)2007—2018年農(nóng)作物因旱受災(zāi)面積、受災(zāi)人口、因旱飲水困難牲畜及旱災(zāi)直接經(jīng)濟(jì)損失情況分別與錫林郭勒草原干旱頻率做相關(guān)性分析,并將其與不同程度的干旱頻率做皮爾遜顯著性檢驗(yàn),如圖3和表3所示??梢园l(fā)現(xiàn),年際干旱發(fā)生頻率與農(nóng)作物受旱面積、因旱飲水困難人畜及直接經(jīng)濟(jì)損失等均呈明顯正相關(guān);因研究時(shí)段錫林郭勒草原主要以輕中旱為主,不同干旱等級(jí)中,中旱頻率與農(nóng)作物受旱面積和直接經(jīng)濟(jì)損失呈顯著性相關(guān),顯著性水平分別達(dá)到0.05和0.01,輕旱頻率與受災(zāi)人口呈顯著性相關(guān),相關(guān)性達(dá)0.05顯著性水平。

    圖3 錫林郭勒草原干旱發(fā)生頻率與內(nèi)蒙古地區(qū)旱災(zāi)損失相關(guān)性分析

    表3 內(nèi)蒙古自治區(qū)旱災(zāi)損失與錫林郭勒草原旱情等級(jí)的顯著性檢驗(yàn)

    注:* 在 0.05 級(jí)別(雙尾)相關(guān)性顯著;** 在 0.01 級(jí)別(雙尾)相關(guān)性顯著。下同。

    Note: * significant correlation at 0.05 level (two tailed); ** significant correlation at 0.01 level (two tailed). The same below.

    3.2.2 錫林郭勒草原干旱的季節(jié)性時(shí)空分布特征

    對(duì)錫林郭勒不同草原類型2007—2018年春、夏和秋3個(gè)季節(jié)干旱特征及季節(jié)性干旱頻率時(shí)空分布特征進(jìn)行分析,如圖4和圖5所示??梢钥闯?,錫林郭勒草原輕中度干旱發(fā)生頻率從強(qiáng)到弱表現(xiàn)為秋季、春季、夏季,空間分布由西北向東南逐漸下降,呈明顯的條帶狀,全區(qū)干旱嚴(yán)重區(qū)域主要集中在西北部;重旱及以上干旱主要發(fā)生在春秋兩季,夏季受降水量集中且雨強(qiáng)較大等影響,重度以上干旱基本沒(méi)有,重旱嚴(yán)重區(qū)主要集中在錫林郭勒典型草原中北部阿巴嘎旗和那仁寶力格。按草原類型分析,荒漠草原的季節(jié)性干旱頻率顯著高于典型草原和草甸草原,但典型草原輕度干旱發(fā)生頻繁,荒漠草原中度及以上干旱相對(duì)頻發(fā),典型草原秋季輕中旱和春季重特旱頻率與荒漠草原相當(dāng)。

    3.3 綜合干旱指數(shù)與單一干旱指數(shù)的敏感性分析

    為直觀反映綜合干旱指數(shù)對(duì)干旱的敏感性,選取能充分體現(xiàn)降水盈虧程度的降水距平百分率(Precipitation Anomaly Percentage)氣象干旱指數(shù),將其與葉綠素?zé)晒庵笖?shù)(SIF)、標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)、植被供水指數(shù)(VSWI)和綜合干旱指數(shù)(SSV)繪制在隨時(shí)間變化的同一坐標(biāo)系統(tǒng)中,如圖6所示。

    從圖6可知,單一干旱指數(shù)隨降水距平百分率指數(shù)的變化趨勢(shì)均比較顯著,雨量越充沛,干旱指數(shù)越大,干旱程度越小。其中,標(biāo)準(zhǔn)化降水蒸散指數(shù)屬于氣象干旱指數(shù),一定程度上對(duì)干旱的早期開(kāi)始過(guò)程比較敏感,但其波動(dòng)范圍較大,適用于短歷時(shí)旱情監(jiān)測(cè),對(duì)長(zhǎng)序列旱情的精準(zhǔn)監(jiān)測(cè)有失偏頗;以最大和最小降水距平百分率指數(shù)為參照,發(fā)現(xiàn)植被供水指數(shù)和葉綠素?zé)晒庵笖?shù)因反映的是土壤水分供給情況及植被光合作用效率,相比降水距平百分率指數(shù)而言,最大存在1~2個(gè)月的滯后期,而綜合干旱指數(shù)SSV在耦合土壤水分和植被有效光合效率的基礎(chǔ)上,對(duì)降水虧缺的短歷時(shí)變化和長(zhǎng)序列演變響應(yīng)均比較及時(shí),相比單一干旱指數(shù)而言,不僅能捕捉氣象干旱的早期開(kāi)始,亦能從多要素反映農(nóng)牧業(yè)干旱的持續(xù)過(guò)程,一定程度上彌補(bǔ)了單一干旱指數(shù)對(duì)降水響應(yīng)不及時(shí)的缺陷,進(jìn)一步驗(yàn)證了綜合干旱指數(shù)的有效性和可靠性。

    選取2016年和2017年為典型年,對(duì)比分析綜合干旱指數(shù)相比降水距平百分率指數(shù)和土壤相對(duì)濕潤(rùn)度指數(shù)(0~20 cm)在監(jiān)測(cè)精度上的改進(jìn)。2016年,錫林郭勒草原大部分地區(qū)均有不同程度的干旱,主要集中在西烏珠穆沁旗、東烏珠穆沁旗和蘇尼特右旗,受災(zāi)人口分別為4.107 7萬(wàn)、1.679 4萬(wàn)和1.351 7萬(wàn)人,農(nóng)作物受旱面積為220.68萬(wàn)、325.19萬(wàn)、184.38 萬(wàn)hm2,直接經(jīng)濟(jì)損失為148 959萬(wàn)、9 750 萬(wàn)、260萬(wàn)元,相比較而言,西烏珠穆沁旗干旱較東烏珠穆沁旗嚴(yán)重,且蘇尼特右旗干旱特征突出。采用2016年錫林郭勒草原綜合干旱指數(shù)和降水距平百分率指數(shù)分析干旱時(shí)空分布(圖7),可以發(fā)現(xiàn)綜合干旱指數(shù)與實(shí)際旱情更加吻合,且基于降水距平百分率指數(shù)的草甸草原旱情較綜合干旱指數(shù)評(píng)估的旱情偏重,荒漠草原則偏輕。

    受降雨持續(xù)偏少、氣溫異常偏高及大風(fēng)天氣較多的影響,錫林郭勒草原阿巴嘎旗于2017年4月1日—9月7日發(fā)生持續(xù)性干旱,直接經(jīng)濟(jì)損失達(dá)17 000萬(wàn)元;5月1日—8月29日和9月7日,東烏珠穆沁旗和蘇尼特左旗發(fā)生旱災(zāi),直接經(jīng)濟(jì)損失分別為59 000萬(wàn)、11 270萬(wàn)元;4月10和5月26日,錫林浩特市和西烏珠穆沁旗分別開(kāi)始干旱,且錫林浩特市為春夏連旱,直接經(jīng)濟(jì)損失分別為6 213萬(wàn)、28 298萬(wàn)元。以2017年旱災(zāi)直接經(jīng)濟(jì)損失為參照,干旱程度表現(xiàn)為東烏珠穆沁旗>西烏珠穆沁旗>阿巴嘎旗>蘇尼特左旗>錫林浩特。表4為2017年4 —9 月,錫林郭勒草原干旱地區(qū)基于綜合干旱指數(shù)和土壤相對(duì)濕度指數(shù)每16 d的旱情等級(jí)結(jié)果,分析發(fā)現(xiàn):綜合干旱指數(shù)監(jiān)測(cè)旱情等級(jí)相比土壤相對(duì)濕度指數(shù)偏重,尤其對(duì)荒漠草原;綜合干旱指數(shù)與土壤相對(duì)濕度指數(shù)相比實(shí)際旱情均存在一定偏差,但干旱程度還與農(nóng)作物受災(zāi)面積、受災(zāi)人口等有關(guān),無(wú)法單純依靠直接經(jīng)濟(jì)損失判別,且綜合干旱指數(shù)監(jiān)測(cè)干旱的開(kāi)始與持續(xù)時(shí)間,相比土壤相對(duì)濕度指數(shù)而言,與實(shí)際旱情更吻合。

    表4 基于綜合干旱指數(shù)與土壤相對(duì)濕度指數(shù)的旱情對(duì)比分析

    3.4 綜合干旱指數(shù)在不同類型草原中的適用強(qiáng)度分析

    不同類型草原受植被蓋度、植被光合效率及土壤下滲率等不一的影響,對(duì)降水和土壤水分的時(shí)滯效應(yīng)存在差異。本研究分析不同草原類型綜合干旱指數(shù)對(duì)降水和土壤相對(duì)濕度的響應(yīng)規(guī)律有利于精準(zhǔn)描述旱情演化特征。錫林郭勒草甸草原、典型草原和荒漠草原的綜合干旱指數(shù)與降水距平百分率指數(shù)相關(guān)性如圖8所示;不同草原類型綜合干旱指數(shù)與不同深度土壤相對(duì)濕度指數(shù)的相關(guān)顯著性檢驗(yàn)如表5所示。

    表5 不同深度土壤相對(duì)濕度指數(shù)與綜合干旱指數(shù)的相關(guān)性顯著檢驗(yàn)

    從綜合干旱指數(shù)與降水距平百分率指數(shù)的相關(guān)性可以發(fā)現(xiàn),草甸草原和典型草原長(zhǎng)序列綜合干旱指數(shù)與降水距平百分率指數(shù)的變化趨勢(shì)和相關(guān)性較荒漠草原顯著;荒漠草原因生物多樣性較少、生產(chǎn)力較低、土壤養(yǎng)分貧瘠且水分虧缺嚴(yán)重等特性,對(duì)一次強(qiáng)降水變化的響應(yīng)更為敏感,對(duì)基于葉綠素?zé)晒獾木C合干旱指數(shù)響應(yīng)不顯著,但不同草原類型綜合干旱指數(shù)與土壤相對(duì)濕度指數(shù)顯著相關(guān)(<0.01),且 0~20 cm深度的土壤相對(duì)濕度與綜合干旱指數(shù)相關(guān)性最顯著。因此,綜合干旱指數(shù)適用于錫林郭勒不同草原類型的旱情監(jiān)測(cè),但荒漠草原相對(duì)草甸草原和典型草原的適用強(qiáng)度偏弱。

    4 討 論

    干旱的精準(zhǔn)化監(jiān)測(cè)有利于準(zhǔn)確把握旱情動(dòng)態(tài)、評(píng)估區(qū)域干旱狀況,提升農(nóng)牧業(yè)生產(chǎn)力水平。傳統(tǒng)的干旱監(jiān)測(cè)包括地面降水、氣溫和蒸散等氣象要素監(jiān)測(cè)和衛(wèi)星遙感大尺度干旱指數(shù)反演,將地面點(diǎn)狀與遙感大尺度數(shù)據(jù)相結(jié)合,能有效反應(yīng)干旱時(shí)空尺度的異質(zhì)性。本研究采用權(quán)重組合方法,從致旱因子相互響應(yīng)機(jī)制出發(fā),融合地面降水和蒸散發(fā)、遙感植被光合效率及土壤水分等水循環(huán)過(guò)程因子,構(gòu)建適用于不同下墊面條件的綜合干旱指數(shù),從氣象、水文及農(nóng)業(yè)角度綜合刻畫(huà)干旱,進(jìn)而實(shí)現(xiàn)旱情的精準(zhǔn)捕捉。

    降水是干旱的主要驅(qū)動(dòng)因子,一般而言,氣象干旱指數(shù)能及時(shí)反映降水的虧缺狀況,進(jìn)而對(duì)干旱進(jìn)行早期預(yù)警,而農(nóng)業(yè)干旱和水文干旱則存在一定的滯后性。本研究融合新型葉綠素?zé)晒庵笖?shù)一定程度上彌補(bǔ)了農(nóng)業(yè)干旱滯后性影響,與單一干旱指數(shù)和氣象干旱指數(shù)Pa對(duì)比,融合多源信息的綜合干旱指數(shù)不僅能有效捕捉干旱的早期開(kāi)始,亦能多方面反映干旱的持續(xù)過(guò)程及旱情等級(jí)。以錫林郭勒草原為例,2007—2018年,基于綜合干旱指數(shù)的錫林郭勒草原旱情主要以輕中旱為主,秋季和春季的干旱強(qiáng)度較強(qiáng),對(duì)人畜飲水困難及農(nóng)作物受災(zāi)等產(chǎn)生直接影響,這與杜波波等[31-32]結(jié)論相契合,進(jìn)一步揭示了綜合干旱指數(shù)的可靠性。西部二連浩特和蘇尼特左旗荒漠草原土壤貧瘠,保水性差,植被類型多為戈壁針茅和石生針茅,與中東部典型草原和草甸草原對(duì)比,其對(duì)降水過(guò)程的響應(yīng)更加敏感;荒漠草原牧草對(duì)一次有效降水過(guò)程(≥10 mm)的生理響應(yīng)短于草甸草原,一般為7~12 d[33],草原類型越干旱,其植被蓋度對(duì)降水的滯后效應(yīng)和累積效應(yīng)越明顯[34],除此之外,荒漠草原受長(zhǎng)序列蒸散發(fā)強(qiáng)耗損及植被蓋度對(duì)降水滯后效應(yīng)等影響,綜合干旱指數(shù)與降水距平百分率指數(shù)和土壤相對(duì)濕度指數(shù)的相關(guān)系數(shù)相對(duì)較低,但正相關(guān)性顯著,一定程度上揭示了不同草原類型對(duì)氣候因子響應(yīng)程度不一,只有綜合氣象、土壤及植被等多要素因子,才能精準(zhǔn)刻畫(huà)旱情變化,進(jìn)而評(píng)估旱情等級(jí)。

    5 結(jié) 論

    本研究基于錫林郭勒草原2007—2018年每16 d一次的標(biāo)準(zhǔn)化降水蒸散指數(shù)、植被供水指數(shù)及葉綠素?zé)晒庵笖?shù),構(gòu)建了適用于不用類型草原和植被不同生長(zhǎng)季的綜合干旱指數(shù),從干旱指數(shù)敏感性及可靠性等方面分析了綜合干旱指數(shù)的適用性,在此基礎(chǔ)上,探究了錫林郭勒草原干旱的年際及季節(jié)性變化特征,得出以下主要結(jié)論:

    1)綜合干旱指數(shù)從干旱機(jī)理出發(fā),對(duì)干旱信息進(jìn)行了最優(yōu)組合,輸入信息易于獲取和計(jì)算,相比單一干旱指數(shù)在旱情的定性和定量分析中具有較大優(yōu)越性,融合葉綠素?zé)晒庵笖?shù)的綜合干旱指數(shù)在錫林郭勒不同草原類型干旱分析中具有適用性和可靠性。

    2)2007年以來(lái),錫林郭勒草原主要以輕旱和中旱為主,呈現(xiàn)由西向東逐漸下降的條帶狀,重旱及以上干旱集中在春季和秋季,空間上以中部典型草原和荒漠草原過(guò)渡帶頻發(fā)為主。

    3)荒漠草原對(duì)降水過(guò)程的響應(yīng)比較敏感,植被蓋度對(duì)降水的滯后效應(yīng)和累積效應(yīng)顯著,使用單一干旱指數(shù)監(jiān)測(cè)其旱情,一定程度上會(huì)造成實(shí)際情況的夸大和不實(shí),而綜合水熱、土壤及植被立地條件等探究干旱過(guò)程,有利于實(shí)現(xiàn)精準(zhǔn)監(jiān)測(cè)。

    [1] 張強(qiáng),韓蘭英,張立陽(yáng),等. 論氣候變暖背景下干旱和干旱災(zāi)害風(fēng)險(xiǎn)特征與管理策略[J]. 地球科學(xué)進(jìn)展,2014,29(1):80-91.

    Zhang Qiang, Han Lanying, Zhang Liyang, et al. Analysis on the character and management strategy of drought disaster and risk under the climatic warming[J]. Advances in Earth Science, 2014, 29(1): 80-91. (in Chinese with English abstract)

    [2] King A D, Pitman A J, Henley B J, et al. The role of climate variability in Australian drought[J]. Nature Climate Change, 2020, 10(3): 177-179.

    [3] 吳循,周青. 氣候變暖對(duì)陸地生態(tài)系統(tǒng)的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2008,16(1):223-228.

    Wu Xun, Zhou Qing. Effect of climate change and global warming on terrestrial ecosystems[J]. Chinese Journal of Eco-Agriculture, 2008, 16(1): 223-228. (in Chinese with English abstract)

    [4] 屈艷萍,呂娟,蘇志誠(chéng),等. 抗旱減災(zāi)研究綜述及展望[J]. 水利學(xué)報(bào),2018,49(1):115-125.

    Qu Yanping, Lv Juan, Su Zhicheng, et al. Research review and perspective of drought mitigation[J]. Journal of Hydraulic Engineering, 2018, 49(1): 115-125. (in Chinese with English abstract)

    [5] 許世衛(wèi). 中國(guó)農(nóng)業(yè)監(jiān)測(cè)預(yù)警的研究進(jìn)展與展望[J]. 農(nóng)學(xué)學(xué)報(bào),2018,8(1):197-202.

    Xu Shiwei. Review on research in agricultural monitoring and early warning in China[J]. Journal of Agriculture, 2018, 8(1): 197-202. (in Chinese with English abstract)

    [6] Hao Z, Singh V P. Drought characterization from a multivariate perspective: A review[J]. Journal of Hydrology, 2015, 527: 668-678.

    [7] 趙福年,王潤(rùn)元,王鶯,等. 干旱過(guò)程、時(shí)空尺度及干旱指數(shù)構(gòu)建機(jī)制的探討[J]. 災(zāi)害學(xué),2018,33(4):32-39.

    Zhao Funian, Wang Runyuan, Wang Ying, et al. Discussion of process, temporal and spatial scale for drought and establishment of drought indices[J]. Journal of Catastrophology, 2018, 33(4): 32-39. (in Chinese with English abstract)

    [8] 吳志勇,程丹丹,何海,等. 綜合干旱指數(shù)研究進(jìn)展[J]. 水資源保護(hù),2021,37(1):36-45.

    Wu Zhiyong, Cheng Dandan, He Hai, et al. Research progress of composite drought index[J]. Water Resources Protection, 2021, 37(1): 36-45. (in Chinese with English abstract)

    [9] 溫慶志,孫鵬,張強(qiáng),等. 基于多源遙感數(shù)據(jù)的農(nóng)業(yè)干旱監(jiān)測(cè)模型構(gòu)建及應(yīng)用[J]. 生態(tài)學(xué)報(bào),2019,39(20):7757-7770.

    Wen Qingzhi, Sun Peng, Zhang Qiang, et al. An integrated agricultural drought monitoring model based on multi-source remote sensing data: model development and application[J]. Acta Ecologica Sinica, 2019, 39(20): 7757-7770. (in Chinese with English abstract)

    [10] 杜靈通,田慶久,王磊,等. 基于多源遙感數(shù)據(jù)的綜合干旱監(jiān)測(cè)模型構(gòu)建[J]. 農(nóng)業(yè)工程學(xué)報(bào),2014,30(9):126-132.

    Du Lingtong, Tian Qingjiu, Wang Lei, et al. A synthesized drought monitoring model based on multi-source remote sensing data[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(9): 126-132. (in Chinese with English abstract)

    [11] 李軍,吳旭樹(shù),王兆禮,等. 基于新型綜合干旱指數(shù)的珠江流域未來(lái)干旱變化特征研究[J]. 水利學(xué)報(bào),2021,52(4):486-497.

    Li Jun, Wu Xushu, Wang Zhaoli, et al. Changes of drought characteristics in future in Pearl River Basin describing by a new comprehensive standardized drought index[J]. Journal of Hydraulic Engineering, 2021, 52(4): 486-497. (in Chinese with English abstract)

    [12] 王璽圳,粟曉玲,張更喜. 綜合干旱指數(shù)的構(gòu)建及其在涇惠渠灌區(qū)的應(yīng)用[J]. 干旱地區(qū)農(nóng)業(yè)研究,2019,37(4):223-230.

    Wang Xizhen, Su Xiaoling, Zhang Gengxi. Construction of comprehensive drought index and its application in Jinghui irrigation area[J]. Agricultural Research in the Arid Areas, 2019, 37(4): 223-230. (in Chinese with English abstract)

    [13] 王思楠,李瑞平,李夏子. 基于綜合干旱指數(shù)的毛烏素沙地腹部土壤水分反演及分布[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(13):113-121.

    Wang Sinan, Li Ruiping, Li Xiazi. Inversion and distribution of soil moisture in belly of Maowusu sandy land based on comprehensive drought index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(13): 113-121. (in Chinese with English abstract)

    [14] 章釗穎,王松寒,邱博,等. 日光誘導(dǎo)葉綠素?zé)晒膺b感反演及碳循環(huán)應(yīng)用進(jìn)展[J]. 遙感學(xué)報(bào),2019,23(1):37-52.

    Zhang Zhaoying, Wang Songhan, Qiu Bo, et al. Retrieval of sun-induced chlorophyll fluorescence and advancements in carbon cycle application[J]. Journal of Remote Sensing, 2019, 23(1): 37-52. (in Chinese with English abstract)

    [15] 王雅楠,韋瑾,湯旭光,等. 應(yīng)用葉綠素?zé)晒夤浪阒脖豢偝跫?jí)生產(chǎn)力研究進(jìn)展[J]. 遙感技術(shù)與應(yīng)用,2020,35(5):975-989.

    Wang Yanan, Wei Jin, Tang Xuguang, et al. Progress of using the chlorophyll fluorescence to estimate terrestrial gross primary production[J]. Remote Sensing Technology and Application, 2020, 35(5): 975-989. (in Chinese with English abstract)

    [16] 劉雷震,武建軍,周洪奎,等. 葉綠素?zé)晒饧捌湓谒置{迫監(jiān)測(cè)中的研究進(jìn)展[J]. 光譜學(xué)與光譜分析,2017,37(9):2780-2787.

    Liu Leizhen, Wu Jianjun, Zhou Hongkui, et al. Chlorophyll fluorescence and its progress in detecting water stress[J]. Spectroscopy and Spectral Analysis, 2017, 37(9): 2780-2787. (in Chinese with English abstract)

    [17] Wang H, Xiao J. Improving the capability of the SCOPE model for simulating solar-induced fluorescence and gross primary production using data from OCO-2 and flux towers[J]. Remote Sensing, 2021, 13(4): 794.

    [18] 錢(qián)新. 基于葉綠素?zé)晒膺b感的植被干旱脅迫指數(shù)構(gòu)建[D]. 南京:南京大學(xué),2019.

    Qian Xin. A Vegetation Drought Stress Index Based on Solar-Induced Chlorophyll Fluorescence[D]. Najing: Najing University, 2019. (in Chinese with English abstract)

    [19] 陳家寧,孫懷衛(wèi),王建鵬,等. 綜合氣象干旱指數(shù)改進(jìn)及其適用性分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(16):71-77.

    Chen Jianing, Sun Huaiwei, Wang Jianpeng, et al. Improvement of comprehensive meteorological drought index and its applicability analysis[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(16): 71-77. (in Chinese with English abstract)

    [20] 江笑薇,白建軍,劉憲峰. 基于多源信息的綜合干旱監(jiān)測(cè)研究進(jìn)展與展望[J]. 地球科學(xué)進(jìn)展,2019,34(3):275-287.

    Jiang Xiaowei, Bai Jianjun, Liu Xianfeng. Research progress and prospect of integrated drought monitoring based on multisource information[J]. Advances in Earth Science, 2019, 34(3): 275-287. (in Chinese with English abstract)

    [21] 張建平,劉宗元,王靖,等. 西南地區(qū)綜合干旱監(jiān)測(cè)模型構(gòu)建與驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(5):102-107.

    Zhang Jianping, Liu Zongyuan, Wang Jing, et al. Construction and validation of comprehensive drought monitoring model in Southwest China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(5): 102-107. (in Chinese with English abstract)

    [22] 賈若楠,杜鑫,李強(qiáng)子,等. 近15年錫林郭勒盟植被變化時(shí)空特征及其對(duì)氣候的響應(yīng)[J]. 中國(guó)水土保持科學(xué),2016,14(5):47-56.

    Jia Ruonan, Du Xin, Li Qiangzi, et al. Spatiotemporal characteristics of vegetation variation in Xilin Gol League in recent 15 years and its response to climate[J]. Science of Soil and Water Conservation, 2016, 14(5): 47-56. (in Chinese with English abstract)

    [23] 趙水霞,王文君,吳英杰,等. 近59a錫林郭勒草原旱災(zāi)驅(qū)動(dòng)氣候因子分析[J]. 干旱區(qū)研究,2021,38(3):785-793.

    Zhao Shuixia, Wang Wenjun, Wu Yingjie, et al. Analysis of drought-driving climatic factors of Xilin Gol grassland in the past 59 years[J]. Arid Zone Research, 2021, 38(3): 785-793. (in Chinese with English abstract)

    [24] Vicente-Serrano S M, Begueria S, Lopez-Moreno J I. A multiscale drought index sensitive to global warming: The standardized precipitation evapotranspiration index[J]. Journal of Climate, 2010, 23(7): 1696-1718.

    [25] Javed T, Zhang J, Bhattarai N, et al. Drought characterization across agricultural regions of China using standardized precipitation and vegetation water supply indices[J]. Journal of Cleaner Production, 2021, 313:127866. doi:10.1016/j. jclepro.2021.127866

    [26] 張永江,劉良云,侯名語(yǔ),等. 植物葉綠素?zé)晒膺b感研究進(jìn)展[J]. 遙感學(xué)報(bào),2009,13(5):963-978.

    Zhang Yongjiang, Liu Liangyun, Hou Mingyu, et al. Progress in remote sensing of vegetation chlorophyll fluorescence[J]. National Remote Sensing Bulletin, 2009, 13(5): 963-978. (in Chinese with English abstract)

    [27] 張立福,王思恒,黃長(zhǎng)平. 太陽(yáng)誘導(dǎo)葉綠素?zé)晒獾男l(wèi)星遙感反演方法[J]. 遙感學(xué)報(bào),2018,22(1):1-12.

    Zhang Lifu, Wang Siheng, Huang Changping. Top-of-atmosphere hyperspectral remote sensing of solar-induced chlorophyll fluorescence: A review of methods[J]. National Remote Sensing Bulletin, 2018, 22(1): 1-12. (in Chinese with English abstract)

    [28] Li X, Xiao J F. A global, 0.05-degree product of solar-induced chlorophyll fluorescence derived from OCO-2, MODIS, and reanalysis data[J]. Remote Sensing, 2019, 11(5), 517. doi: 10. 3390/rs11050517.

    [29] Qiu R, Han G, X Ma, et al. A comparison of OCO-2 SIF, MODIS GPP, and GOSIF data from gross primary production (GPP) estimation and seasonal cycles in North America[J]. Remote Sensing, 2020, 12(2): 258.

    [30] 崔園園,張強(qiáng),李威,等. CLDAS融合土壤相對(duì)濕度產(chǎn)品適用性評(píng)估及在氣象干旱監(jiān)測(cè)中的應(yīng)用[J]. 海洋氣象學(xué)報(bào),2020,40(4):105-113.

    Cui Yuanyuan, Zhang Qiang, Li Wei, et al. Applicability evaluation of CLDAS merged soil relative moisture product and its application in meteorological drought monitoring[J]. Journal of Marine Meteorology, 2020, 40(4): 105-113. (in Chinese with English abstract)

    [31] 杜波波,阿拉騰圖婭,包剛. 2002-2016年錫林郭勒草原干旱時(shí)空特征[J].水土保持研究,2019,26(4):190-194,202.

    Du Bobo, Alatengtuya, Bao Gang. Temporal and spatial characteristics of drought in the Xilingol grassland from 2002 to 2016[J]. Research of Soil and Water Conservation, 2019, 26(4): 190-194, 202. (in Chinese with English abstract)

    [32] 張巧鳳,劉桂香,于紅博,等. 錫林郭勒草原土壤含水量遙感反演模型及干旱監(jiān)測(cè)[J]. 草業(yè)學(xué)報(bào),2017,26(11):1-11.

    Zhang Qiaofeng, Liu Guixiang, Yu Hongbo, et al. Soil moisture modelling and drought monitoring using remote sensing in Xilingol grassland[J]. Acta Prataculturae Sinica, 2017, 26(11): 1-11. (in Chinese with English abstract)

    [33] 劉志剛,劉麗萍,游曉勇,等. 錫林郭勒草原氣候變化與干旱特征[J]. 內(nèi)蒙古氣象,2008(1):17-18.

    Liu Zhigang, Liu Liping, You Xiaoyong, et al. Characteristics of climate change and drought in Xilingol Grassland[J]. Meteorology Journal of Inner Mongolia, 2008(1): 17-18. (in Chinese with English abstract)

    [34] 王舉鳳,何亮,陸紹娟,等. 內(nèi)蒙古不同類型草原光合植被覆蓋度對(duì)降水變化的響應(yīng)[J]. 生態(tài)學(xué)報(bào),2020,40(16):5620-5629.

    Wang Jufeng, He Liang, Lu Shaojuan, et al. Photosynthetic vegetation cover response to precipitation on the Inner Mongolian Steppe[J]. Acta Ecologica Sinica, 2020, 40(16): 5620-5629. (in Chinese with English abstract)

    Construction and application of comprehensive drought index in different steppe types

    Zhao Shuixia, Wang Wenjun※, Wu Yingjie, Li Wei, Quan Qiang

    (1.,,100038,;2.,,010020,)

    The drought index is important to the drought grade evaluation, drought prevention, and mitigation measures formulating.Based on the physical mechanism of drought, the objective weighting method of Criteria Importance Through Inter-criteria Correlation (CRITIC) was used to construct a comprehensive drought index with multi-source information. The mechanism of Standardized Precipitation Evapotranspiration Index (SPEI), Vegetation Supply Water Index (VSWI), and Solar-induced Chlorophyll Fluorescence index (SIF) were comprehensively considered in the drought development process. The vegetation types, coverage, and their sensitivity to precipitation and temperature of different steppes were also been considered in the comprehensive drought index construction. Taking Xinlingol steppe in Inner Mongolia as an example, the applicability and sensitivity of comprehensive drought index in different steppe types were explored compared with the precipitation anomaly percentage index of meteorological drought and the Relative Soil Moisture Index (RSM) of agricultural drought. The temporal and spatial distribution and evolution characteristics of drought were further analyzed in desert, typical and meadow steppes during 2007—2018. The results showed that: 1) The comprehensive drought index can not only sensitively capture the early stage of drought but also comprehensively reflect the duration and time of drought from the perspective of hydrology and pastoral area drought. The comprehensive drought index makes up for the limitation of drought information revealed by a single drought index to a certain extent, which can effectively identify the temporal and spatial distribution characteristics of drought, and improve the pertinence and accuracy of drought monitoring. 2) The temporal and spatial distribution characteristics of drought in Xinlingol steppe evaluated by the comprehensive drought index were more consistent with the actual drought condition compared with precipitation anomaly percentage index and RSM, and the correlation between 0-20 cm RSM and comprehensive drought index was the most significant. The drought frequency of the Xilingol steppe during 2007—2018 showed a fluctuating downward trend, which had a significant positive correlation with the drought crop area, affected population, and direct economic loss. 3) The drought intensity of the Xinlingol steppe was mainly light and moderate drought, in which the severe and extreme drought mainly happened in spring and autumn. For different steppe types of Xilingol, the drought frequency of desert steppe was significantly higher than that of typical steppe and meadow steppe, and the spatial distribution of drought shows a gradually decreasing trend from northwest to southeast and serious drought in the middle and north part of the typical steppe. This study carried out an optimal combination of drought information, in which the input information is easy to obtain and calculate, and the application of solar-induced chlorophyll fluorescence index effectively revealed the response of vegetation to drought stress.The application intensity of the comprehensive drought index in the desert steppe is relatively weak compared with meadow and typical steppe because the soil and biodiversity situation of the desert steppe is poor, the vegetation is sensitive to a precipitation process, and the cumulative effect is significant. This study is of great significance to the drought events monitoring and prevention in different steppe types on the Mongolian Plateau.

    drought; precipitation; comprehensive drought index; chlorophyll fluorescence index; applicability; Xilingol steppe

    趙水霞,王文君,吳英杰,等. 綜合干旱指數(shù)構(gòu)建及其在不同草原類型中的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(16):99-107.doi:10.11975/j.issn.1002-6819.2021.16.013 http://www.tcsae.org

    Zhao Shuixia, Wang Wenjun, Wu Yingjie, et al. Construction and application of comprehensive drought index in different steppe types[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(16): 99-107. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2021.16.013 http://www.tcsae.org

    2021-06-22

    2021-08-10

    內(nèi)蒙古自治區(qū)科技計(jì)劃重大專項(xiàng)及科技計(jì)劃項(xiàng)目(2020ZD0020,2021GG0050);中國(guó)水科院“五大人才”創(chuàng)新團(tuán)隊(duì)項(xiàng)目和基本科研專項(xiàng)(MK0145B022021;MK2020J12)

    趙水霞,博士,工程師,研究方向?yàn)楦珊蛋敫珊祬^(qū)生態(tài)水文。Email:zhaosx@iwhr.com

    王文君,高級(jí)工程師,研究方向?yàn)楦珊蛋敫珊祬^(qū)生態(tài)水文。Email:wangwj@iwhr.com

    10.11975/j.issn.1002-6819.2021.16.013

    P426.616

    A

    1002-6819(2021)-16-0099-09

    猜你喜歡
    旱情葉綠素降水
    人工增雨解旱情
    基于不同旱情指數(shù)的石羊河流域春旱監(jiān)測(cè)研究
    黑龍江省玉米生長(zhǎng)季自然降水與有效降水對(duì)比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    提取葉綠素
    桃樹(shù)葉綠素含量與SPAD值呈極顯著正相關(guān)
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    葉綠素家族概述
    降水現(xiàn)象儀模擬軟件設(shè)計(jì)與實(shí)現(xiàn)
    云南省旱情信息系統(tǒng)升級(jí)及完善
    水利信息化(2015年5期)2015-12-21 12:54:40
    ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
    亚洲成av人片在线播放无| 亚洲国产精品成人综合色| 日韩欧美免费精品| 99热全是精品| 真人做人爱边吃奶动态| 中文字幕av成人在线电影| 亚州av有码| 舔av片在线| 亚洲欧美日韩高清专用| 精品久久久久久久久久免费视频| 成人亚洲欧美一区二区av| 精品久久久久久久久亚洲| 黄色视频,在线免费观看| 亚洲国产精品成人久久小说 | 久久久久久久亚洲中文字幕| 在线天堂最新版资源| 日本 av在线| 国产精品精品国产色婷婷| 午夜爱爱视频在线播放| 欧美极品一区二区三区四区| 精品少妇黑人巨大在线播放 | 亚洲第一区二区三区不卡| 男女做爰动态图高潮gif福利片| 综合色丁香网| 成人二区视频| 日韩国内少妇激情av| 婷婷精品国产亚洲av在线| 男女边吃奶边做爰视频| 中文字幕精品亚洲无线码一区| 欧美一级a爱片免费观看看| 欧美日韩综合久久久久久| avwww免费| 国产高清有码在线观看视频| 国产久久久一区二区三区| 久久久久久久久中文| 亚洲人成网站高清观看| 天堂√8在线中文| 国产乱人视频| 九九在线视频观看精品| 欧美一区二区精品小视频在线| 色综合亚洲欧美另类图片| 听说在线观看完整版免费高清| 午夜老司机福利剧场| 午夜福利视频1000在线观看| 久久99热6这里只有精品| 国内久久婷婷六月综合欲色啪| 日韩av在线大香蕉| 欧美一区二区国产精品久久精品| 久久久久久伊人网av| 色综合站精品国产| 看非洲黑人一级黄片| 黑人高潮一二区| 天堂√8在线中文| 在线观看一区二区三区| 日韩制服骚丝袜av| 亚洲熟妇中文字幕五十中出| 97在线视频观看| 一级毛片电影观看 | 日本黄色片子视频| 欧美成人免费av一区二区三区| 亚洲欧美日韩高清在线视频| 久久鲁丝午夜福利片| 国内精品美女久久久久久| 欧美激情国产日韩精品一区| 国产 一区 欧美 日韩| 校园人妻丝袜中文字幕| 亚洲精品国产成人久久av| 亚洲精品影视一区二区三区av| 国产欧美日韩一区二区精品| 午夜精品在线福利| 欧美成人免费av一区二区三区| 一级毛片久久久久久久久女| 日本一二三区视频观看| 国产精品久久久久久久久免| 99热6这里只有精品| 国产色爽女视频免费观看| 色综合色国产| 免费看a级黄色片| 亚洲自偷自拍三级| 日日干狠狠操夜夜爽| 亚洲成人久久爱视频| 中文字幕久久专区| 一级黄色大片毛片| 菩萨蛮人人尽说江南好唐韦庄 | 成年av动漫网址| 午夜精品一区二区三区免费看| 夜夜看夜夜爽夜夜摸| 国内揄拍国产精品人妻在线| 黑人高潮一二区| 91久久精品国产一区二区成人| 99在线视频只有这里精品首页| 亚洲美女搞黄在线观看 | 亚洲内射少妇av| 日本欧美国产在线视频| 尤物成人国产欧美一区二区三区| 一个人看视频在线观看www免费| 在线a可以看的网站| 九九久久精品国产亚洲av麻豆| 精品午夜福利在线看| 日韩欧美精品免费久久| 精品99又大又爽又粗少妇毛片| 欧美日本视频| 日韩欧美 国产精品| 国产高清视频在线观看网站| 亚洲精品影视一区二区三区av| 精品日产1卡2卡| 久久久欧美国产精品| 一区二区三区高清视频在线| 免费高清视频大片| 免费人成视频x8x8入口观看| 久久精品国产99精品国产亚洲性色| 99热精品在线国产| 全区人妻精品视频| 午夜a级毛片| 校园人妻丝袜中文字幕| 国产av不卡久久| 国产亚洲av嫩草精品影院| 99热这里只有是精品在线观看| 大又大粗又爽又黄少妇毛片口| 网址你懂的国产日韩在线| 日韩欧美精品v在线| 久久久精品欧美日韩精品| 午夜免费激情av| 午夜久久久久精精品| .国产精品久久| 婷婷六月久久综合丁香| 夜夜爽天天搞| 国产淫片久久久久久久久| 色视频www国产| 中文字幕免费在线视频6| 国产伦在线观看视频一区| 亚洲精品亚洲一区二区| 亚洲内射少妇av| 欧美色视频一区免费| 老师上课跳d突然被开到最大视频| 97超碰精品成人国产| 国产av麻豆久久久久久久| 免费观看在线日韩| 午夜影院日韩av| 美女大奶头视频| 日韩三级伦理在线观看| 国产av不卡久久| 久久韩国三级中文字幕| 97超级碰碰碰精品色视频在线观看| 天堂网av新在线| 97超碰精品成人国产| 99久久久亚洲精品蜜臀av| 免费av毛片视频| 精品免费久久久久久久清纯| 亚洲人与动物交配视频| 99精品在免费线老司机午夜| 成人特级黄色片久久久久久久| 3wmmmm亚洲av在线观看| 成人漫画全彩无遮挡| 别揉我奶头~嗯~啊~动态视频| 国产成人a区在线观看| 欧美不卡视频在线免费观看| 午夜影院日韩av| 日韩三级伦理在线观看| 免费搜索国产男女视频| 内地一区二区视频在线| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久久久亚洲av鲁大| 免费不卡的大黄色大毛片视频在线观看 | aaaaa片日本免费| av专区在线播放| 久久精品人妻少妇| 亚洲综合色惰| 亚洲七黄色美女视频| 国产69精品久久久久777片| 亚洲图色成人| 国产伦精品一区二区三区四那| 精品久久久噜噜| 亚洲av中文字字幕乱码综合| 亚洲色图av天堂| av黄色大香蕉| 欧美丝袜亚洲另类| 老司机午夜福利在线观看视频| 亚洲性久久影院| 亚洲真实伦在线观看| 深爱激情五月婷婷| 在线免费十八禁| 中文字幕人妻熟人妻熟丝袜美| 亚洲最大成人中文| 国产爱豆传媒在线观看| 日韩精品有码人妻一区| 欧美激情在线99| 成人亚洲精品av一区二区| 亚洲性久久影院| 国产av在哪里看| 深爱激情五月婷婷| 97超级碰碰碰精品色视频在线观看| 又爽又黄a免费视频| 欧美不卡视频在线免费观看| 欧美日韩在线观看h| 国产一区二区亚洲精品在线观看| 日本一二三区视频观看| 国产探花在线观看一区二区| 成人毛片a级毛片在线播放| 精品午夜福利在线看| 麻豆国产av国片精品| 国产成人a区在线观看| .国产精品久久| 人人妻,人人澡人人爽秒播| 久久久久久大精品| 三级毛片av免费| 欧美色视频一区免费| 天堂动漫精品| 欧美在线一区亚洲| 熟女人妻精品中文字幕| 99热这里只有是精品50| 成人av一区二区三区在线看| 亚洲av中文字字幕乱码综合| 成年免费大片在线观看| 黄色日韩在线| 在线观看一区二区三区| 校园人妻丝袜中文字幕| 欧美中文日本在线观看视频| 欧美激情国产日韩精品一区| 成人三级黄色视频| 精品久久久久久久末码| 精品久久久久久久久久免费视频| 三级男女做爰猛烈吃奶摸视频| 午夜日韩欧美国产| 欧美+亚洲+日韩+国产| 亚洲无线在线观看| 乱系列少妇在线播放| 亚洲国产日韩欧美精品在线观看| 婷婷亚洲欧美| 日日啪夜夜撸| 变态另类丝袜制服| 免费av毛片视频| 简卡轻食公司| 久久久久免费精品人妻一区二区| 97碰自拍视频| 色综合站精品国产| 久久久国产成人精品二区| 欧美潮喷喷水| 99国产极品粉嫩在线观看| 国国产精品蜜臀av免费| 成人午夜高清在线视频| 国产91av在线免费观看| 伊人久久精品亚洲午夜| av在线蜜桃| 欧美+亚洲+日韩+国产| 中国国产av一级| 18禁在线无遮挡免费观看视频 | 赤兔流量卡办理| 成人高潮视频无遮挡免费网站| 国产女主播在线喷水免费视频网站 | 九九在线视频观看精品| 日韩一本色道免费dvd| 国产精品,欧美在线| 成人午夜高清在线视频| 久99久视频精品免费| 国产亚洲精品av在线| 欧美又色又爽又黄视频| 日韩制服骚丝袜av| 亚洲av免费高清在线观看| 青春草视频在线免费观看| 欧美又色又爽又黄视频| 国产人妻一区二区三区在| 成年版毛片免费区| 激情 狠狠 欧美| 亚洲在线自拍视频| eeuss影院久久| 可以在线观看的亚洲视频| 变态另类成人亚洲欧美熟女| 亚洲aⅴ乱码一区二区在线播放| 中文字幕av在线有码专区| 国产精品,欧美在线| 99久国产av精品| 天天躁夜夜躁狠狠久久av| 成人无遮挡网站| 国产精品亚洲美女久久久| 国产高清视频在线观看网站| 国产综合懂色| www.色视频.com| 欧美色欧美亚洲另类二区| 三级国产精品欧美在线观看| 最近手机中文字幕大全| 成人亚洲欧美一区二区av| 美女xxoo啪啪120秒动态图| 国产 一区 欧美 日韩| 看非洲黑人一级黄片| 欧美色欧美亚洲另类二区| 两个人的视频大全免费| 午夜精品在线福利| 在线天堂最新版资源| 69人妻影院| 小蜜桃在线观看免费完整版高清| 国产亚洲91精品色在线| 亚洲av成人精品一区久久| 色噜噜av男人的天堂激情| 国产日本99.免费观看| 久久精品91蜜桃| 亚洲精品影视一区二区三区av| 国产免费男女视频| 嫩草影院精品99| 青春草视频在线免费观看| 久久精品国产清高在天天线| 亚洲美女搞黄在线观看 | 狂野欧美激情性xxxx在线观看| 久久久久久久久久久丰满| 国产一区二区三区av在线 | 欧美一级a爱片免费观看看| 亚洲精品成人久久久久久| 三级毛片av免费| 永久网站在线| 亚洲精品国产av成人精品 | 狠狠狠狠99中文字幕| 国产伦精品一区二区三区四那| 日韩欧美三级三区| 老女人水多毛片| 久久精品国产自在天天线| 午夜激情欧美在线| 久久久久久九九精品二区国产| 国产精品精品国产色婷婷| 国产伦精品一区二区三区四那| 亚洲av美国av| 国产色爽女视频免费观看| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 国产不卡一卡二| 亚洲一区高清亚洲精品| 夜夜看夜夜爽夜夜摸| 岛国在线免费视频观看| 免费看光身美女| 舔av片在线| 精品一区二区三区av网在线观看| 久久久成人免费电影| 国产精品久久久久久久久免| 久久国产乱子免费精品| 99久国产av精品| 无遮挡黄片免费观看| 午夜激情欧美在线| 99国产极品粉嫩在线观看| 国产成人aa在线观看| 欧美日韩综合久久久久久| 色播亚洲综合网| 国产精品乱码一区二三区的特点| 亚洲av第一区精品v没综合| 尾随美女入室| 国产精品女同一区二区软件| av在线亚洲专区| 露出奶头的视频| 国产欧美日韩精品一区二区| 精品久久久久久久久av| 亚洲av中文字字幕乱码综合| 九九热线精品视视频播放| 青春草视频在线免费观看| 3wmmmm亚洲av在线观看| 免费观看人在逋| 一级黄片播放器| 久久精品国产亚洲网站| 久久亚洲国产成人精品v| 国产一区二区三区av在线 | 中文资源天堂在线| 两个人的视频大全免费| 狂野欧美激情性xxxx在线观看| 又爽又黄a免费视频| 日本成人三级电影网站| or卡值多少钱| 日韩亚洲欧美综合| 大香蕉久久网| 久久久久久久午夜电影| 久久人人爽人人片av| 亚洲美女搞黄在线观看 | 91麻豆精品激情在线观看国产| 久久久精品94久久精品| 国产久久久一区二区三区| 在线观看av片永久免费下载| 18禁在线播放成人免费| 欧美成人精品欧美一级黄| 欧美+亚洲+日韩+国产| 六月丁香七月| 18+在线观看网站| www日本黄色视频网| 18禁在线无遮挡免费观看视频 | 日日撸夜夜添| 午夜激情福利司机影院| 免费观看人在逋| 日日摸夜夜添夜夜添小说| 国内揄拍国产精品人妻在线| 在线看三级毛片| 日本黄色视频三级网站网址| 69人妻影院| 欧美潮喷喷水| 大型黄色视频在线免费观看| 成人三级黄色视频| 欧美在线一区亚洲| 日韩欧美国产在线观看| 日韩一本色道免费dvd| 国产欧美日韩精品亚洲av| 中出人妻视频一区二区| 久久精品国产亚洲av涩爱 | 一个人看视频在线观看www免费| 国产午夜福利久久久久久| 毛片一级片免费看久久久久| 91在线观看av| 精品久久久久久久久亚洲| 日韩欧美精品免费久久| 国产成人freesex在线 | 国内少妇人妻偷人精品xxx网站| 嫩草影院新地址| 99精品在免费线老司机午夜| 又粗又爽又猛毛片免费看| 色吧在线观看| 禁无遮挡网站| 不卡视频在线观看欧美| 国产一区二区在线av高清观看| 国产麻豆成人av免费视频| 精品一区二区三区av网在线观看| 看免费成人av毛片| 亚洲精品粉嫩美女一区| 干丝袜人妻中文字幕| 中文在线观看免费www的网站| 国产精品综合久久久久久久免费| 欧美最黄视频在线播放免费| 成年版毛片免费区| 国产精品久久久久久精品电影| 精品无人区乱码1区二区| 中文亚洲av片在线观看爽| 99热只有精品国产| 啦啦啦啦在线视频资源| 在线天堂最新版资源| 国产成人福利小说| 九九爱精品视频在线观看| 黄色日韩在线| 国产精品三级大全| a级毛色黄片| 此物有八面人人有两片| 毛片女人毛片| 精品不卡国产一区二区三区| 一级毛片电影观看 | 国产午夜精品论理片| 香蕉av资源在线| 久久精品国产鲁丝片午夜精品| 男女下面进入的视频免费午夜| 免费看日本二区| 午夜福利在线在线| 午夜精品在线福利| 在线观看66精品国产| 午夜福利在线观看吧| 最近在线观看免费完整版| 久久久久久久久久成人| 国产精品爽爽va在线观看网站| 啦啦啦观看免费观看视频高清| 亚州av有码| 91精品国产九色| 国产精品综合久久久久久久免费| 久久九九热精品免费| 在线观看66精品国产| 成人午夜高清在线视频| 日日摸夜夜添夜夜添av毛片| 长腿黑丝高跟| 久久国内精品自在自线图片| 亚洲三级黄色毛片| 丰满的人妻完整版| 精品午夜福利视频在线观看一区| 男人舔奶头视频| 国产人妻一区二区三区在| 你懂的网址亚洲精品在线观看 | av黄色大香蕉| 国产高清激情床上av| 91精品国产九色| 国产 一区精品| 男插女下体视频免费在线播放| av天堂中文字幕网| 国模一区二区三区四区视频| 欧美日韩综合久久久久久| 插阴视频在线观看视频| 一级毛片久久久久久久久女| 欧美+日韩+精品| 亚洲婷婷狠狠爱综合网| 国产亚洲精品久久久久久毛片| 熟女人妻精品中文字幕| 日韩中字成人| 国产乱人偷精品视频| 亚洲18禁久久av| 美女被艹到高潮喷水动态| 成人漫画全彩无遮挡| 亚洲欧美日韩东京热| 国产成人a∨麻豆精品| 三级国产精品欧美在线观看| 99热这里只有是精品在线观看| 热99在线观看视频| 国产欧美日韩精品一区二区| 嫩草影院新地址| 国产熟女欧美一区二区| 亚洲人与动物交配视频| 日韩精品有码人妻一区| 欧美潮喷喷水| 亚洲欧美日韩卡通动漫| 人妻夜夜爽99麻豆av| 干丝袜人妻中文字幕| 香蕉av资源在线| 赤兔流量卡办理| 狠狠狠狠99中文字幕| 国产精品亚洲一级av第二区| 成人性生交大片免费视频hd| 美女大奶头视频| 男女那种视频在线观看| 乱人视频在线观看| 国产精品亚洲一级av第二区| 免费在线观看成人毛片| 搡女人真爽免费视频火全软件 | 一a级毛片在线观看| 亚洲人成网站高清观看| 超碰av人人做人人爽久久| 亚洲精品一卡2卡三卡4卡5卡| 简卡轻食公司| 男人舔女人下体高潮全视频| 亚洲人成网站高清观看| 日韩成人av中文字幕在线观看 | 国产人妻一区二区三区在| 国产 一区精品| 国产精品一区二区免费欧美| 黄色视频,在线免费观看| 欧美日韩综合久久久久久| 最近2019中文字幕mv第一页| 久久欧美精品欧美久久欧美| 欧洲精品卡2卡3卡4卡5卡区| 香蕉av资源在线| av在线老鸭窝| 在线观看免费视频日本深夜| 亚洲精品亚洲一区二区| av女优亚洲男人天堂| 日本免费一区二区三区高清不卡| 在线免费十八禁| 日韩欧美 国产精品| 亚洲av免费高清在线观看| 亚洲av五月六月丁香网| 国产精品女同一区二区软件| 日韩在线高清观看一区二区三区| 色综合站精品国产| 午夜日韩欧美国产| 嫩草影院精品99| 在线播放无遮挡| 国产黄色小视频在线观看| 内射极品少妇av片p| 色播亚洲综合网| 小说图片视频综合网站| 高清毛片免费观看视频网站| 亚洲人成网站在线播放欧美日韩| 插阴视频在线观看视频| 日本成人三级电影网站| 国产亚洲精品久久久com| 人人妻人人澡人人爽人人夜夜 | 久久人人爽人人爽人人片va| 国产av麻豆久久久久久久| 国国产精品蜜臀av免费| 中文字幕av在线有码专区| 欧美性猛交╳xxx乱大交人| 尾随美女入室| 天美传媒精品一区二区| 欧美绝顶高潮抽搐喷水| 嫩草影院入口| 你懂的网址亚洲精品在线观看 | 日韩高清综合在线| 亚洲成人精品中文字幕电影| 精品一区二区三区av网在线观看| 两个人视频免费观看高清| 国产精品综合久久久久久久免费| 高清毛片免费看| 国产一区二区亚洲精品在线观看| 变态另类成人亚洲欧美熟女| 桃色一区二区三区在线观看| 亚洲欧美成人综合另类久久久 | 最近最新中文字幕大全电影3| 国产欧美日韩精品亚洲av| 久久鲁丝午夜福利片| 成人av在线播放网站| 久久婷婷人人爽人人干人人爱| 午夜免费男女啪啪视频观看 | 成年女人看的毛片在线观看| 日本黄色片子视频| 美女免费视频网站| 成人综合一区亚洲| 最近中文字幕高清免费大全6| 老熟妇仑乱视频hdxx| 亚洲第一电影网av| 欧美性感艳星| 寂寞人妻少妇视频99o| 中文字幕免费在线视频6| 露出奶头的视频| 亚洲电影在线观看av| 人妻久久中文字幕网| 成年女人永久免费观看视频| 波野结衣二区三区在线| 久久精品国产亚洲网站| 久久精品国产清高在天天线| 国产伦精品一区二区三区四那| 亚洲国产精品成人综合色| 99国产极品粉嫩在线观看| 天堂网av新在线| 亚洲av一区综合| 国产精品一区二区性色av| 国产久久久一区二区三区| 深爱激情五月婷婷| 人妻夜夜爽99麻豆av| 国产精品永久免费网站| 舔av片在线| 午夜福利在线观看免费完整高清在 | 女人被狂操c到高潮| 99在线人妻在线中文字幕| 一级黄色大片毛片| 蜜桃亚洲精品一区二区三区| 久久久久久久午夜电影| 日本黄色视频三级网站网址| 中文字幕免费在线视频6| 成人漫画全彩无遮挡| 日韩欧美 国产精品| 日韩欧美国产在线观看| 岛国在线免费视频观看|