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

    基于SPEI的新疆干旱變化特征分析

    2022-06-16 03:27:55阿迪拉?阿布都熱合曼昝梅
    安徽農(nóng)業(yè)科學 2022年11期
    關(guān)鍵詞:干旱

    阿迪拉?阿布都熱合曼 昝梅

    摘要 以新疆55個國家氣象站1961—2015年逐日氣溫和降水數(shù)據(jù),采用季尺度下的標準化降水蒸散指數(shù)(SPEI),以國家劃分的氣象干旱等級為判別標準,分析新疆干旱過程。結(jié)果表明,近55年新疆輕旱次數(shù)呈西北多東南少的趨勢分布;中旱次數(shù)呈西北和東北多、南部偏少的趨勢;重旱次數(shù)呈準噶爾盆地、吐魯番盆地和塔里木盆地東部干旱次數(shù)偏高,其余地區(qū)較低的分布趨勢;特旱次數(shù)呈南多北少的分布趨勢。近55年新疆干旱影響范圍及干旱強度極差較大,且兩者呈正相關(guān)。農(nóng)作物平均受災率與影響范圍年際變化具有較高的一致性。通過REOF時空分布將新疆分成3個模態(tài):第一模態(tài)呈北疆與南疆相反格局的特征;第二模態(tài)高值區(qū)位于北疆北部,負值區(qū)位于新疆西南部;第三模態(tài)表現(xiàn)為新疆東西部呈相反格局的狀態(tài)。各模態(tài)干旱強度隨時間變化趨勢也不同。

    關(guān)鍵詞 標準化降水蒸散指數(shù)(SPEI);干旱;REOF;時空分布特征

    中圖分類號 P429? 文獻標識碼 A? 文章編號 0517-6611(2022)11-0178-06

    doi:10.3969/j.issn.0517-6611.2022.11.046

    開放科學(資源服務)標識碼(OSID):

    Analysis of Drought Change Characteristics in Xinjiang Based on SPEI

    Adila·Abudureheman1, ZAN Mei1,2

    (1.School of Geography and Tourism, Xinjiang Normal University, Urumqi,Xinjiang 830054;2.Xinjiang Key Laboratory of Lake Environment and Resources in Arid Areas, Urumqi, Xinjiang 830054)

    Abstract Based on the daily temperature and precipitation data from 55 national meteorological stations in Xinjiang from 1961 to 2015,the standardized precipitation evapotranspiration index (SPEI) at the seasonal scale was used to analyze the drought process in Xinjiang with the meteorological drought grades classified by the country as the criterion.The results showed that the number of light droughts in Xinjiang in the past 55 years showed a trend of more in the northwest and less in the southeast;the number of moderate droughts was more in the northwest and northeast and less in the south.The number of severe drought was higher in the Junggar Basin, Turpan Basin and eastern Tarim Basin, and lower in other areas. The number of special drought was more in the south and less in the north.The range and intensity of drought in Xinjiang in the last 55 years were significantly different, and the two were positively correlated. The average disaster rate of crops had a high consistency with the interannual variation of the affected area.Xinjiang was divided into three modes based on the spatial and temporal distribution of REOF. The first mode was the opposite pattern of northern Xinjiang and southern Xinjiang. The high value region of the second mode was located in the north of northern Xinjiang, and the negative value region was located in the southwest of Xinjiang. The third mode showed that the pattern of east and west of Xinjiang was opposite. The variation trend of drought intensity of each mode was also different with time.

    Key words Standardized precipitation evapotranspiration index (SPEI);Drought;REOF;Temporal and spatial distribution characteristics

    干旱是自然界中頻發(fā)的一種災害現(xiàn)象,其影響范圍廣,持續(xù)時間長,造成的災難性后果嚴重影響到我國農(nóng)業(yè)生產(chǎn)、糧食安全與社會經(jīng)濟發(fā)展[1-2]。IPCC報告中指出,自20世紀以來全球普遍升溫,蒸發(fā)量大幅度增加,干旱趨勢增強[3]。我國干旱、半干旱區(qū)開展不同尺度的區(qū)域干旱檢測與風險評估有利于改善地區(qū)干旱現(xiàn)象、提高地區(qū)生態(tài)系統(tǒng)的自我恢復能力。干旱指數(shù)是開展區(qū)域干旱檢測與風險評估的有效手段之一[4]。

    鑒于造成干旱現(xiàn)象的因素錯綜復雜,為了更好地研究干旱,也為了能更準確地表示干旱的嚴重程度,各國學者都提出了專門用于干旱監(jiān)測的指標。當前關(guān)于干旱指標的研究較多,其中應用較廣泛的氣象干旱指標有標準化降水指數(shù)(SPI)、標準化降水蒸散指數(shù)(SPEI)、氣象干旱綜合指數(shù)(MCI)、帕默爾干旱指數(shù)(PDSI)等,各指標都有各自不同的優(yōu)劣勢[5-9]。Vicente-Serrano等[10]提出的SPEI指數(shù),是在SPI指數(shù)的基礎上,結(jié)合SPI、PDSI指數(shù)的優(yōu)點,綜合考慮了降水的蒸散發(fā)影響,且從多時間尺度角度比較分析,適用性更佳,逐步成為氣象干旱監(jiān)測中較理想的工具之一。目前利用不同干旱指標對新疆進行干旱變化特征的研究較多,但利用季尺度數(shù)據(jù)對研究區(qū)進行干旱時空變化特征及干旱災情的研究較少,3個月時間尺度可反映季節(jié)干旱狀況,12和24個月可分別反映年際、年代際干旱狀況。新疆作為我國國土空間面積最大的區(qū)域,地處西北邊陲,氣候干燥、降水稀少且水資源十分匱乏,干旱現(xiàn)象頻發(fā),嚴重制約了當?shù)厣鐣?jīng)濟發(fā)展,加大了與我國內(nèi)陸地區(qū)的發(fā)展差距,基于此,筆者利用3個月尺度的SPEI指數(shù),定量分析新疆干旱現(xiàn)象發(fā)生的時空演化規(guī)律,厘清新疆發(fā)生干旱的基本特征,為新疆的干旱監(jiān)控、預警、防范與減災工作提供理論借鑒和科學依據(jù)。

    1 資料與方法

    1.1 研究區(qū)概況

    新疆維吾爾自治區(qū)地處73°29′54″~96°23′03″E、34°20′11″~ 49°10′55″N,位于我國西北。全疆總面積為166×104 km2,占全國總面積的1/6[11]。新疆屬于典型的溫帶大陸性氣候,年均降水量遠低于沿海地區(qū)。山脈與盆地相間排列構(gòu)成了新疆“三山夾兩盆”的地貌特征,由北至南分別是阿爾泰山、準噶爾盆地、天山山脈、塔里木盆地和昆侖山。以天山山脈為界形成了南北疆差異顯著的氣候特征。在全球氣候變暖的背景下,加劇了新疆的干旱特征,進而造成干旱災害頻繁發(fā)生,嚴重威脅著新疆社會經(jīng)濟發(fā)展。

    1.2 數(shù)據(jù)來源

    該研究的氣象站觀測數(shù)據(jù)來源于中國氣象科學數(shù)據(jù)中心(http://data.cma.cn/site/index.html),選取了新疆55個國家氣象站點的平均氣溫和降水等數(shù)據(jù)。數(shù)據(jù)來源準確、可靠性高、研究期數(shù)據(jù)完整。其中SPEI數(shù)據(jù)集及各站點經(jīng)緯度、海拔等數(shù)據(jù)來源于相關(guān)文獻[12]。該研究選用的新疆區(qū)域氣象站點分布如圖1所示。

    1.3 研究方法

    1.3.1

    標準化降水蒸散指數(shù)(SPEI)的計算方法[13]。

    (1)利用Thornthwaite方法,綜合考慮氣溫因素的影響來計算潛在蒸散量(PET),公式如下:

    PET=16.0× 10TiHa (1)

    式中, Ti為月平均氣溫,H為年熱量總和,a為常數(shù)。

    (2)計算降水與潛在蒸散量的差值。計算公式如下:

    Di = Pi -PET i (2)

    式中, Pi 為月降水量,PET i 為月潛在蒸散量。

    (3)對數(shù)據(jù)序列進行正態(tài)化處理,由于 Di 序列可能存有負值,采用log-logistic概率分布對 Di 進行擬合,并求累計函數(shù)。計算公式如下:

    F(x)=1+αx-γβ-1(3)

    式中,α、β、γ 分別為尺度參數(shù)、形狀參數(shù)、Origin參數(shù),均用線性矩方法擬合得到。

    (4)利用正態(tài)逆變換求出SPEI。計算公式如下:

    W =-2ln( P )(4)

    SPEI= W-C0-C1W+C2W21+d1W+d2W2+d3W3(5)

    式中,C0=2.515 517,C1=0.802 853,C2=0.010 328,d1=1.432 788,d2=0.189 269,d3=0.001 308。當P>0.5時,P值取1-P [14]。

    1.3.2 干旱強度的運算。

    據(jù)國家氣候中心最新發(fā)布的氣象干旱等級(GB/T 20481—2006)[15]中對干旱過程及強度的定義:當SPEI指數(shù)持續(xù)10 d為輕旱以上等級時,可被判定為一次干旱過程。干旱強度定義為一次干旱過程累計的SPEI指數(shù)之和,其值越小,干旱過程越強。SPEI指數(shù)具有多時間尺度的特征,SPEI-3即季尺度的數(shù)值,能夠較明確地表示季節(jié)性土壤濕度的演變特征。統(tǒng)計新疆地區(qū)1961—2015年SPEI-3的數(shù)值,結(jié)合新疆特有地形地貌條件,以及前人的研究也表明SPEI 指數(shù)具有較好的區(qū)域適用性,因此SPEI 指數(shù)作為干旱評估研究標準是合理的,結(jié)果也是可參考的。SPEI指數(shù)的干旱等級劃分見表1。

    1.3.3 干旱影響范圍的確定。

    基于站點的SPEI對干旱影響范圍進行計算。SPEI值達到輕旱及以上等級的臺站數(shù)( Ni)與總站數(shù)(N )之比,即SPEI值≤-0.5時就表示有干旱的發(fā)生,作為不同級別干旱影響范圍的指標,其反映干旱面積大小。將干旱影響范圍超過90%的定義為大范圍干旱。計算公式為:

    Ri=NiN×100%(6)

    1.3.4 正交函數(shù)方法。經(jīng)驗正交函數(shù)(EOF)和旋轉(zhuǎn)經(jīng)驗正交函數(shù)(REOF)是分析氣象要素時空分布特征時廣泛使用的統(tǒng)計方法[16-19]。正交函數(shù)方法能將同時隨時間及空間變化的某要素場進行分解,將其分解為僅隨空間變化的空間模態(tài)和僅隨時間變化的時間系數(shù)的線性組合。采用EOF方法,根據(jù)EOF原理進行分解,對任意標準化氣象矩陣 X 可展開下列矩陣運算[20]:

    mXn=(mVk×λk)×1λk×kTn=mLk×kFn(7)

    式中,mLk是空間特征向量矩陣,kTn是時間系數(shù),mVk是空間特征向量矩陣,kFn是主成分陣,λk為mXn的相關(guān)系數(shù)陣mRn的特征值。

    REOF基于EOF對載荷特征向量場進行方差極大旋轉(zhuǎn)變換,使旋轉(zhuǎn)空間向量場上各空間點旋轉(zhuǎn)一個或幾個,與EOF方法相比,它解決了不能明確表示不同地理區(qū)域特征的局限。旋轉(zhuǎn)后所得的典型空間分布更好地反映了不同地理區(qū)域的變化情況,也反映了不同區(qū)域的相關(guān)分布狀況。該研究對干旱過程強度進行REOF分析,所采用的空間插值方法為克里金插值。

    2 結(jié)果與分析

    2.1 新疆不同等級干旱次數(shù)空間分布

    從新疆不同等級年干旱次數(shù)空間分布(圖2)可以看出,1961—2015年各氣象站點輕旱等級年平均干旱次數(shù)大體呈西北多東南少的趨勢分布,發(fā)生輕旱次數(shù)最少的站點為十三間房,為274次。哈密地區(qū)西部及吐魯番市輕旱次數(shù)偏少,小于299次;研究區(qū)南部焉耆、拜城,北部的布爾津、鞏留等16個站點輕旱次數(shù)在298~311次;研究區(qū)西北、西南及東部的吉木乃、塔城、烏魯木齊、莎車等23個站點發(fā)生輕旱次數(shù)在311~323次;輕旱次數(shù)在323~335次的12個站點多數(shù)分布在天山山脈及阿爾泰山脈周邊,其余少數(shù)站點全疆零星散布。

    研究區(qū)各站點中旱等級年干旱次數(shù)的分布趨勢并不明顯,大體呈西北和東北多、南部偏少的趨勢。發(fā)生中旱次數(shù)最少的站點為巴音郭楞蒙古自治州西部的輪臺站(48次);發(fā)生中旱次數(shù)最多的站點為伊犁州的新源站(89次);發(fā)生中旱次數(shù)較少(48~64次)的站點大多分布在南疆地區(qū),站點數(shù)為13個;中旱次數(shù)在64~80次的站點均勻分布在全疆各地,發(fā)生中旱次數(shù)較多的站點在北疆的布爾津、新源以及南疆的拜城等7個站點。

    研究區(qū)各站點重旱等級年干旱次數(shù)空間分布明顯,準噶爾盆地、吐魯番盆地周邊及塔里木盆地東部地區(qū)干旱次數(shù)偏高,其余地區(qū)較低。發(fā)生重旱次數(shù)最多的站點為哈密地區(qū)的十三間房站點(62次),最少的站點為和靜站(25次)。發(fā)生重旱次數(shù)在25~39次的站點數(shù)為38個,發(fā)生重旱次數(shù)在39~61次的站點數(shù)為15個。

    研究區(qū)各站點特旱等級年干旱次數(shù)呈南多北少的趨勢分布,發(fā)生干旱次數(shù)最多的站點為莎車站(29次),最少的站點為十三間房(5次)。吉木乃、博樂、巴音布魯克、塔什庫爾干等21 個站點發(fā)生干旱次數(shù)為5~14次;發(fā)生干旱次數(shù)在14~19次的站點數(shù)為17個;發(fā)生重旱次數(shù)較多的站點集中在新疆南部地區(qū),有15個站點發(fā)生干旱次數(shù)在19~28次。

    2.2 新疆干旱影響范圍年際變化特征及干旱災情分析

    根據(jù)氣象干旱過程及干旱強度的定義,經(jīng)站點比例計算得到近55年新疆地區(qū)干旱影響范圍。由圖3a可知,研究區(qū)干旱影響范圍( Ri )極差較大,超過90%;影響較小的干旱范圍在1987、1992和1993年出現(xiàn),分別為0、3.6%、0;影響較大的干旱范圍出現(xiàn)在1997和2008年,分別為98.1%、92.7%;近55年平均干旱影響范圍為46.3%。從平均累計強度(SPEI)分析干旱過程的年際特征,由圖3b可知,干旱強度的極差也較大,干旱強度最弱的年份出現(xiàn)在1987年,累計強度為529.5;最強的干旱強度年出現(xiàn)在1997年,為-625.3。干旱影響范圍與區(qū)域平均累計強度年際變化特征彼此一致,說明大多數(shù)情況下干旱強度及其影響范圍呈現(xiàn)正相關(guān),干旱強度越強,影響范圍越廣,干旱強度越弱其影響范圍也越小。

    旱災對新疆農(nóng)業(yè)發(fā)展影響較大。氣象、農(nóng)業(yè)及水文干旱為干旱的3種類型,其中氣象干旱與農(nóng)業(yè)干旱的關(guān)聯(lián)性較高,從某種程度上來說,農(nóng)業(yè)干旱是氣象干旱的進一步發(fā)展。受旱面積和受災率作為農(nóng)業(yè)干旱指標,能夠較全面地反映農(nóng)業(yè)受災情況[21]。因資料年代受限,在此僅討論1990—2008年全疆受旱面積情況[22]。從全疆1990—2008年災情平均情況(圖4)來看,2008年農(nóng)作物年均受旱面積最大,超過80×104 hm2;2002年全疆農(nóng)作物受旱面積最小,在2×104 hm2以下,其余年份受旱面積均在2×104~50×104 hm2。由于新疆面積廣,氣候條件特殊,因而導致年份間受災面積相差較大,該研究以農(nóng)作物受旱面積占受災面積的百分率來代表當年農(nóng)業(yè)干旱受災率,分析干旱對農(nóng)業(yè)生產(chǎn)產(chǎn)生的影響。從1990—2008年農(nóng)作物平均受災率(圖4)來看,2004、2006、2008年均超過了50%,其對應的干旱影響范圍分別為50.9%、76.3%、92.7%,這與上述的影響范圍年際變化特征具有較高的一致性;平均受災率較低的年份出現(xiàn)在1992、1993、1998、1999、2002和2003年,這些年份干旱影響范圍分別為3.6%、0、30.9%、67.2%、45.4%、14.5%,其中1992、1993和2002年的平均受災率與干旱影響范圍年際變化特征一致性較高,部分年份一致性較低。SPEI指數(shù)能夠較好地反映不同時間尺度干旱的評估,但由于SPEI主要考慮的是氣象因子,沒有考慮到其他人為因素的影響,因而可能導致了部分年份一致性較低的情況。

    2.3 新疆干旱強度時空特征

    2.3.1 干旱強度空間模態(tài)。

    根據(jù)該研究定義得到近55年來新疆出現(xiàn)的所有干旱過程,計算新疆55個氣象站點逐年所有干旱過程的累計強度,將MCI指數(shù)矩陣先進行EOF分解,然后再對前5個分量進行REOF處理并分析,得到REOF結(jié)果(表2)。從表2可看出,REOF前5個模態(tài)累計方差貢獻率達到了66.37%。

    圖5為對新疆SPEI指數(shù)進行REOF分析的前3個空間模態(tài)的分布情況,從圖中可看出,第一模態(tài)表現(xiàn)出以北天山為分界線,新疆北部與南部呈現(xiàn)相反格局的空間分布特征;

    REOF特征向量值最小值為-0.015,方差貢獻率為24.76%,高值區(qū)位于東疆地區(qū)和南疆地區(qū),說明基于季尺度的SPEI來說,南疆為新疆干旱指數(shù)變化敏感的中心地區(qū)。第二模態(tài)表現(xiàn)為正中心區(qū)在北疆阿勒泰、塔城地區(qū)以及準噶爾盆地北部,位于研究區(qū)東南方向的庫木庫勒盆地為過渡區(qū),負值中心在新疆西南以及東部部分地區(qū),方差貢獻率為14.24%。第三模態(tài)表現(xiàn)為新疆西部與東部呈現(xiàn)相反格局的模態(tài),高值區(qū)在新疆西部,其位于天山西部高值區(qū)的伊犁河谷地區(qū)的水汽主要來自大西洋,而南疆地區(qū)近些年來的濕潤化也較為明顯,據(jù)前人研究,由于新疆獨特的地形作用,北冰洋和西伯利亞沼澤水汽通過阿爾泰山和天山的缺口,導致水汽從中灌入,進而使南疆西部出現(xiàn)較濕潤的氣候。

    2.3.2 各模態(tài)時間系數(shù)變化。

    據(jù)REOF分解出來的3個空間模態(tài)對應的時間序列,結(jié)合“2.2”干旱影響范圍年際變化特征的內(nèi)容,分析近55年來新疆各模態(tài)干旱強度隨時間的變化特征。

    REOF1突出表現(xiàn)了南疆的正值區(qū),對應的時間序列(圖6a)在1961—1965年處于正負值間變化的趨勢;1966—1996年基本上處于在正值區(qū)域飄浮變化的趨勢,說明在這段時間新疆進入了相對干旱緩解的階段。其中達到正值區(qū)的典型年份是1972、1976、1987、1992和1996年;1997—2015年處于在負值區(qū)域變化的趨勢,這說明這段時間新疆變干趨勢明顯,典型干旱年份為2001、2007—2009、2013年,從圖5a、圖6a可看出,南北疆處于變化相反的格局,自20世紀90年代末期,北疆處于較干旱的變化趨勢。

    REOF2空間分布特征(圖5b)與REOF1相差較大,第二模態(tài)高值區(qū)位于北疆北部地區(qū),對應的時間序列(圖6b)在1961—1967年處于正負值間變化不定的趨勢;1968—1973年處于在正值間變化的趨勢,其典型年份為1969年,說明在這段時間干旱有變緩趨勢;1974—2001年處于長期在正負值間變換的階段,北疆西北部、南疆除東南部外的大部分地區(qū)都處于干旱狀態(tài),典型干旱年份為1982、1997年。

    REOF3空間模態(tài)(圖5c)表現(xiàn)為新疆東西部呈相反的格局,高值區(qū)主要位于北疆西部、南疆西南部。由時間序列變化趨勢(圖6c)可知,1961—1980年處于正負值間振蕩變化的趨勢;1981—2015年,除少數(shù)年份外,大部分年份處于正值區(qū)域,典型年份為1987年。這說明20世紀80年代后,干旱趨勢得到了明顯緩解,且造成了東干西濕的局面,近些年來西北地區(qū)的濕潤化也確實得到了前人的驗證[23-24]。

    3 討論與結(jié)論

    標準化降水蒸散指數(shù)(SPEI)為主要考慮蒸散發(fā)影響的氣象干旱指數(shù),新疆地形條件復雜、氣候多變,該研究僅從氣象角度來分析干旱特征,并未從灌溉或社會經(jīng)濟因素考慮。因此該研究的干旱影響范圍與農(nóng)作物實際受旱情況一致性有待探討??傮w上,通過

    SPEI分析得出,1961—2015年新疆不同等級年干旱次數(shù)空間分布明顯;通過計算SPEI干旱強度,較好地反映了近55年來新疆干旱影響范圍年際變化特征與干旱災情之間的關(guān)系;基于SPEI指數(shù),用REOF方法將新疆分成了3個模態(tài),并與其對應的時間系數(shù)進行了對比討論。具體得出以下結(jié)論:

    (1)近55年新疆輕旱次數(shù)呈西北多東南少的趨勢分布;中旱次數(shù)分布趨勢不明顯,大體呈西北和東北多、南部偏少的趨勢;重旱次數(shù)空間分布明顯,準噶爾盆地、吐魯番盆地及塔里木盆地東部干旱次數(shù)偏高,其余地區(qū)較低;特旱次數(shù)呈南多北少的分布趨勢。

    (2)近55年新疆干旱影響范圍及干旱強度極差較大,且兩者呈正相關(guān)。從1990—2008年災情情況來看,2008年農(nóng)作物受旱面積最大,2002年受旱面積最小。農(nóng)作物平均受災率與影響范圍年際變化具有較高的一致性。

    (3)通過REOF將新疆分成3個模態(tài),第一模態(tài)呈北疆與南疆相反格局的特征,其時間系數(shù)表明,自20世紀90年代末期以來,北疆處于較干旱的變化趨勢;第二模態(tài)高值區(qū)位于北疆北部,負值區(qū)位于新疆西南部,時間系數(shù)表明其經(jīng)歷了干旱緩解期到長期干旱的變化;第三模態(tài)表現(xiàn)為新疆東西部呈相反格局的狀態(tài),其時間系數(shù)表明20世紀80年代后期干旱趨勢得到了明顯的緩解。

    參考文獻

    [1] 沈國強,鄭海峰,雷振鋒.基于SPEI指數(shù)的1961—2014年東北地區(qū)氣象干旱時空特征研究[J].生態(tài)學報,2017,37(17):5882-5893.

    [2] 張棋,許德合,丁嚴.基于SPEI和時空立方體的中國近40年干旱時空模式挖掘[J].干旱地區(qū)農(nóng)業(yè)研究,2021,39(3):194-201.

    [3] 滕懷頤,馮克鵬.基于SPEI干旱指數(shù)的西北地區(qū)干旱時空分布[J].農(nóng)業(yè)與技術(shù),2021,41(8):887-893.

    [4] 楊金虎,張強,王勁松,等.近60a來中國西南春季持續(xù)性干旱異常特征分析[J].干旱區(qū)地理,2015,38(2):215-222.

    [5] 李治國,朱玲玲,張延偉,等.基于SPI指數(shù)的近55年河南省干旱時空變化特征[J].江蘇農(nóng)業(yè)科學,2018,46(10):237-242.

    [6] 謝坤,丁明軍,張華敏,等.基于SPEI指數(shù)的鄱陽湖流域近60 a干旱時空變化特征[J].長江流域資源與環(huán)境,2021,30(8):2015-2027.

    [7] 成青燕,高曉清,林紓,等.基于MCI指標的甘肅省近50年干旱特征分析[J].干旱地區(qū)農(nóng)業(yè)研究,2017,35(1):211-218.

    [8] 宋琳琳,張強,任余龍,等.PDSI及sc_PDSI干旱指數(shù)在中國西南地區(qū)適用性分析[J].中國沙漠,2021,41(2):242-251.

    [9] 沈彥軍,李紅軍,雷玉平.干旱指數(shù)應用研究綜述[J].南水北調(diào)與水利科技,2013,11(4):128-133,186.

    [10]? VICENTE-SERRANO S M,BEGUERA S,LPEZ-MORENO J I,et al.A multiscalar drought index sensitive to global warming:The standardized precipitation evapotranspiration index[J].Journal of climate,2010,23(7):1696-1718.

    [11] 潘素敏,張明軍,王圣杰.新疆夏季云下二次蒸發(fā)對雨滴穩(wěn)定同位素影響的定量研究[J].干旱區(qū)地理,2018,41(3):488-498.

    [12] 姚俊強,毛煒嶧,胡文峰,等.1961—2015年新疆區(qū)域SPEI干旱指數(shù)數(shù)據(jù)集[J].中國國家冰川凍土沙漠科學數(shù)據(jù),2019,4(3):112-121.

    [13] 張菡,張喜亮,李金建,等.基于SPEI的四川省盆地區(qū)季節(jié)性干旱時空變化特征分析[J].干旱地區(qū)農(nóng)業(yè)研究,2018,36(5):242-250,256.

    [14] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局,中國國家標準化管理委員會.氣象干旱等級:GB/T 20481—2017[S].北京:中國標準出版社,2017.

    [15] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局,中國國家標準化管理委員會.氣象干旱等級:GB/T 20481—2006[S].北京:中國標準出版社,2006.

    [16] NORTH G R,BELL T L,CAHALAN R F,et al.Sampling errors in the estimation of empirical orthogonal functions[J].Monthly weather review,1982,110(7):699-706.

    [17] 施能.氣象科研與預報中的多元分析方法[M].北京:氣象出版社,1995.

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

    [19] 程麟鈞,王帥,宮正宇,等.中國臭氧濃度的時空變化特征及分區(qū)[J].中國環(huán)境科學,2017,37(11):4003-4012.

    [20] 李家啟,申雙和,劉俊,等.重慶地區(qū)地閃活動空間分布及其趨勢分析[J].西南大學學報(自然科學版),2012,34(5):77-84.

    [21] 趙海燕,張強,高歌,等.中國1951—2007年農(nóng)業(yè)干旱的特征分析[J].自然災害學報,2010,19(4):201-206.

    [22] 新疆維吾爾自治區(qū)統(tǒng)計局,國家統(tǒng)計局新疆調(diào)查總隊.新疆統(tǒng)計年鑒[M].北京:中國統(tǒng)計出版社,2018:33-155.

    [23] 姚旭陽,張明軍,張宇,等.中國西北地區(qū)氣候轉(zhuǎn)型的新認識[J/OL].干旱區(qū)地理,2021-10-18[2021-11-22].http://kns.cnki.net/kcms/detail/65.1103.X.20211018.1255.004.html.

    [24] 李明,孫洪泉,蘇志誠.中國西北氣候干濕變化研究進展[J].地理研究,2021,40(4):1180-1194.

    猜你喜歡
    干旱
    探討夏季極端天氣對園林植物生長的影響及應對措施
    區(qū)域動態(tài)氣象干旱強度指數(shù)與應用探討
    科學與財富(2017年3期)2017-03-15 17:06:45
    基于距平的白城地區(qū)干旱時間分布特征分析
    臨夏地區(qū)干旱特征及干濕氣候區(qū)劃
    干旱鹽堿共脅迫下玉米miR398的表達
    夏季高溫干旱時節(jié)高山蔬菜種植管理策略
    小說里的世界 虛擬里的真實
    基于多源衛(wèi)星遙感的長江流域旱情監(jiān)測研究
    模擬干旱對草地早熟禾種子萌發(fā)的影響
    基于SPI指數(shù)的農(nóng)作物生長期干旱時間變化研究
    国产亚洲欧美在线一区二区| 亚洲成av片中文字幕在线观看| 国产精品国产高清国产av | 亚洲第一青青草原| 久久久久网色| 最新美女视频免费是黄的| 免费人妻精品一区二区三区视频| 精品一区二区三卡| 老熟妇乱子伦视频在线观看| av天堂久久9| 午夜精品久久久久久毛片777| 久久久水蜜桃国产精品网| 欧美 日韩 精品 国产| 99九九在线精品视频| 免费不卡黄色视频| 大陆偷拍与自拍| 久热这里只有精品99| 国产成人欧美| 热re99久久国产66热| 免费久久久久久久精品成人欧美视频| 日本av免费视频播放| 女性生殖器流出的白浆| 亚洲国产毛片av蜜桃av| 久热爱精品视频在线9| 国产精品一区二区在线不卡| 如日韩欧美国产精品一区二区三区| 亚洲熟妇熟女久久| 国产男靠女视频免费网站| 欧美 亚洲 国产 日韩一| √禁漫天堂资源中文www| 视频区图区小说| 亚洲国产精品一区二区三区在线| 免费黄频网站在线观看国产| 老司机深夜福利视频在线观看| 久久ye,这里只有精品| 最黄视频免费看| 天天添夜夜摸| 精品一品国产午夜福利视频| 成人特级黄色片久久久久久久 | 欧美 亚洲 国产 日韩一| 最新美女视频免费是黄的| 99久久精品国产亚洲精品| 51午夜福利影视在线观看| 人人妻人人爽人人添夜夜欢视频| 成人国产一区最新在线观看| 国产精品久久久久成人av| 成人三级做爰电影| a在线观看视频网站| 亚洲中文av在线| 欧美变态另类bdsm刘玥| 不卡一级毛片| 日日摸夜夜添夜夜添小说| 日本一区二区免费在线视频| 亚洲一码二码三码区别大吗| 国产三级黄色录像| 一区福利在线观看| 国产野战对白在线观看| 国产一区二区三区视频了| 热99国产精品久久久久久7| 久久久久久久久免费视频了| 久久久精品国产亚洲av高清涩受| 啪啪无遮挡十八禁网站| 19禁男女啪啪无遮挡网站| 日韩 欧美 亚洲 中文字幕| 在线亚洲精品国产二区图片欧美| 夜夜骑夜夜射夜夜干| 看免费av毛片| bbb黄色大片| 一区二区三区精品91| 肉色欧美久久久久久久蜜桃| 久久久欧美国产精品| 亚洲欧美色中文字幕在线| 国产av一区二区精品久久| 悠悠久久av| 久久精品亚洲精品国产色婷小说| 夫妻午夜视频| 国产成人av教育| 一级片免费观看大全| 90打野战视频偷拍视频| 一本一本久久a久久精品综合妖精| 国产男靠女视频免费网站| 999精品在线视频| 一级黄色大片毛片| 久久国产亚洲av麻豆专区| 多毛熟女@视频| 国产精品秋霞免费鲁丝片| 超色免费av| 男女免费视频国产| 午夜福利欧美成人| 999精品在线视频| 777久久人妻少妇嫩草av网站| 9色porny在线观看| 狂野欧美激情性xxxx| 黄色片一级片一级黄色片| 久久免费观看电影| 老司机影院毛片| 午夜福利免费观看在线| 久久久欧美国产精品| 欧美av亚洲av综合av国产av| 999久久久国产精品视频| 成人特级黄色片久久久久久久 | 另类精品久久| 国产视频一区二区在线看| 在线观看舔阴道视频| 中文字幕av电影在线播放| 大片电影免费在线观看免费| 午夜老司机福利片| 桃红色精品国产亚洲av| 18禁美女被吸乳视频| 人妻一区二区av| 少妇精品久久久久久久| 日韩欧美国产一区二区入口| 麻豆乱淫一区二区| 在线 av 中文字幕| 一区二区三区国产精品乱码| 黄片播放在线免费| 久久久精品94久久精品| 国产高清videossex| 色播在线永久视频| 一本大道久久a久久精品| 精品福利观看| 欧美激情久久久久久爽电影 | 久久久久网色| 在线观看免费视频日本深夜| 女性被躁到高潮视频| 欧美精品亚洲一区二区| 亚洲精品在线美女| 欧美日韩黄片免| www.自偷自拍.com| 亚洲 国产 在线| 日韩欧美国产一区二区入口| 亚洲欧美色中文字幕在线| av欧美777| 国产免费av片在线观看野外av| 99久久人妻综合| 999久久久国产精品视频| 另类精品久久| 一边摸一边抽搐一进一出视频| 黑人操中国人逼视频| 少妇猛男粗大的猛烈进出视频| 欧美精品av麻豆av| 精品久久蜜臀av无| 亚洲成国产人片在线观看| 熟女少妇亚洲综合色aaa.| av又黄又爽大尺度在线免费看| 欧美日韩视频精品一区| 视频在线观看一区二区三区| 亚洲七黄色美女视频| 精品免费久久久久久久清纯 | av不卡在线播放| 无限看片的www在线观看| 一级毛片精品| 亚洲欧洲精品一区二区精品久久久| 桃花免费在线播放| 69av精品久久久久久 | 日韩成人在线观看一区二区三区| 黄色成人免费大全| 九色亚洲精品在线播放| 午夜激情av网站| 五月天丁香电影| 免费日韩欧美在线观看| 一边摸一边抽搐一进一出视频| 九色亚洲精品在线播放| 少妇 在线观看| 亚洲午夜精品一区,二区,三区| 一本综合久久免费| 精品国产国语对白av| 一区二区三区国产精品乱码| 在线观看舔阴道视频| 久久精品亚洲熟妇少妇任你| 亚洲av片天天在线观看| 脱女人内裤的视频| 国产人伦9x9x在线观看| 久久ye,这里只有精品| 久久ye,这里只有精品| 高清在线国产一区| 搡老乐熟女国产| 高潮久久久久久久久久久不卡| 91老司机精品| 成人影院久久| 正在播放国产对白刺激| 成人18禁高潮啪啪吃奶动态图| 我的亚洲天堂| 怎么达到女性高潮| 美女福利国产在线| 香蕉久久夜色| av免费在线观看网站| 亚洲中文字幕日韩| 国产精品98久久久久久宅男小说| av不卡在线播放| 久久精品亚洲精品国产色婷小说| 一本综合久久免费| 欧美性长视频在线观看| 老司机深夜福利视频在线观看| 中文字幕高清在线视频| 不卡av一区二区三区| 水蜜桃什么品种好| 国产视频一区二区在线看| 每晚都被弄得嗷嗷叫到高潮| 精品亚洲成a人片在线观看| 国产又爽黄色视频| 欧美黄色片欧美黄色片| 国产三级黄色录像| 欧美日本中文国产一区发布| 捣出白浆h1v1| 久久中文字幕一级| 男女免费视频国产| 超色免费av| 国产av又大| 国产区一区二久久| 国产aⅴ精品一区二区三区波| 啦啦啦在线免费观看视频4| 亚洲av成人一区二区三| 99国产精品一区二区三区| 久久国产精品影院| 日韩中文字幕欧美一区二区| 韩国精品一区二区三区| 在线观看一区二区三区激情| 一进一出好大好爽视频| 亚洲精品国产精品久久久不卡| 香蕉丝袜av| 中亚洲国语对白在线视频| 久久香蕉激情| 亚洲成a人片在线一区二区| 18在线观看网站| 精品国产一区二区久久| 久久久久久亚洲精品国产蜜桃av| 国产精品亚洲av一区麻豆| 欧美成人免费av一区二区三区 | 日韩 欧美 亚洲 中文字幕| 亚洲av美国av| 免费观看av网站的网址| 丝袜人妻中文字幕| 超色免费av| 色婷婷av一区二区三区视频| 18禁裸乳无遮挡动漫免费视频| www.熟女人妻精品国产| 日本av手机在线免费观看| 欧美乱码精品一区二区三区| 热99国产精品久久久久久7| 中文字幕精品免费在线观看视频| 在线观看66精品国产| 国产视频一区二区在线看| 黄网站色视频无遮挡免费观看| 女人高潮潮喷娇喘18禁视频| av视频免费观看在线观看| 精品亚洲成a人片在线观看| 久久精品熟女亚洲av麻豆精品| 91字幕亚洲| 十八禁网站网址无遮挡| 国产精品 国内视频| 国产精品久久久久久人妻精品电影 | 国产亚洲精品第一综合不卡| 日本欧美视频一区| 欧美乱码精品一区二区三区| 美女国产高潮福利片在线看| 中文亚洲av片在线观看爽 | 九色亚洲精品在线播放| 老司机靠b影院| 欧美在线黄色| 99精品欧美一区二区三区四区| 久久久久久久精品吃奶| 无人区码免费观看不卡 | 国产日韩欧美亚洲二区| 久久精品aⅴ一区二区三区四区| 岛国毛片在线播放| 亚洲精品中文字幕一二三四区 | 精品久久久精品久久久| 欧美国产精品va在线观看不卡| 国产精品一区二区精品视频观看| 极品少妇高潮喷水抽搐| 蜜桃在线观看..| 18禁美女被吸乳视频| 欧美日韩一级在线毛片| 久久午夜综合久久蜜桃| 老司机深夜福利视频在线观看| 国产精品久久久久久精品古装| 久久久久久久国产电影| 成年动漫av网址| 成人免费观看视频高清| 男女下面插进去视频免费观看| 别揉我奶头~嗯~啊~动态视频| 精品欧美一区二区三区在线| 成人精品一区二区免费| 777久久人妻少妇嫩草av网站| 天天躁夜夜躁狠狠躁躁| av视频免费观看在线观看| 国产午夜精品久久久久久| 久久中文字幕人妻熟女| 亚洲av国产av综合av卡| 成人国产av品久久久| 亚洲色图综合在线观看| 色婷婷久久久亚洲欧美| 久久久久久久大尺度免费视频| 丁香欧美五月| av有码第一页| 日本av免费视频播放| 婷婷丁香在线五月| 99精品在免费线老司机午夜| 老司机在亚洲福利影院| 亚洲精品国产区一区二| 搡老岳熟女国产| 日本黄色日本黄色录像| 国产精品成人在线| 久久久久国产一级毛片高清牌| 国精品久久久久久国模美| 国产精品一区二区精品视频观看| 国产野战对白在线观看| 国产一区二区三区综合在线观看| 日本a在线网址| 黄色片一级片一级黄色片| 成人18禁在线播放| 一本综合久久免费| 国产欧美日韩一区二区三区在线| 1024香蕉在线观看| 欧美精品一区二区大全| 精品欧美一区二区三区在线| netflix在线观看网站| 婷婷成人精品国产| 亚洲,欧美精品.| 电影成人av| 大码成人一级视频| 日日夜夜操网爽| 国产精品久久久久久精品电影小说| 热re99久久精品国产66热6| 午夜福利一区二区在线看| 精品国产一区二区三区久久久樱花| 美女主播在线视频| 99精品在免费线老司机午夜| 国产日韩欧美视频二区| av在线播放免费不卡| 午夜福利在线免费观看网站| 久9热在线精品视频| 制服人妻中文乱码| 老司机亚洲免费影院| 精品亚洲乱码少妇综合久久| 精品少妇黑人巨大在线播放| 亚洲精华国产精华精| 成人国语在线视频| 天堂8中文在线网| 欧美性长视频在线观看| 日韩欧美国产一区二区入口| 成人免费观看视频高清| 69精品国产乱码久久久| 亚洲情色 制服丝袜| 亚洲五月色婷婷综合| 久久毛片免费看一区二区三区| 天天操日日干夜夜撸| 熟女少妇亚洲综合色aaa.| 欧美 亚洲 国产 日韩一| 午夜精品国产一区二区电影| 一个人免费在线观看的高清视频| 性少妇av在线| 欧美日韩国产mv在线观看视频| 大型黄色视频在线免费观看| 亚洲精品国产精品久久久不卡| 中文字幕人妻熟女乱码| 免费在线观看日本一区| 中文亚洲av片在线观看爽 | 亚洲中文字幕日韩| 亚洲精品久久成人aⅴ小说| 91成年电影在线观看| 亚洲av日韩精品久久久久久密| 又大又爽又粗| 免费在线观看黄色视频的| 国产99久久九九免费精品| 18禁观看日本| 精品国内亚洲2022精品成人 | 免费人妻精品一区二区三区视频| 国产精品一区二区免费欧美| 亚洲专区国产一区二区| 一进一出好大好爽视频| 老鸭窝网址在线观看| 国产单亲对白刺激| 一区二区三区乱码不卡18| 新久久久久国产一级毛片| 亚洲九九香蕉| 中文字幕色久视频| 在线 av 中文字幕| 久热这里只有精品99| 久久精品成人免费网站| www.精华液| 亚洲三区欧美一区| 少妇被粗大的猛进出69影院| 嫁个100分男人电影在线观看| 两个人看的免费小视频| 丰满饥渴人妻一区二区三| 欧美黑人欧美精品刺激| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦视频在线资源免费观看| 亚洲 国产 在线| 国产男女超爽视频在线观看| 成年女人毛片免费观看观看9 | 国产黄频视频在线观看| 久久青草综合色| 亚洲av成人不卡在线观看播放网| 极品教师在线免费播放| 亚洲第一欧美日韩一区二区三区 | videos熟女内射| 满18在线观看网站| 亚洲欧洲日产国产| 黑人巨大精品欧美一区二区mp4| 另类亚洲欧美激情| 黄片大片在线免费观看| 怎么达到女性高潮| 日韩欧美国产一区二区入口| 天天添夜夜摸| 777米奇影视久久| 麻豆成人av在线观看| 天堂动漫精品| 黑人猛操日本美女一级片| 国产精品久久久av美女十八| 菩萨蛮人人尽说江南好唐韦庄| 丝袜在线中文字幕| 大香蕉久久成人网| 手机成人av网站| 国产男女超爽视频在线观看| 久久人人爽av亚洲精品天堂| 日韩欧美一区视频在线观看| 亚洲成国产人片在线观看| 亚洲全国av大片| 日韩大码丰满熟妇| 美女国产高潮福利片在线看| 亚洲熟妇熟女久久| 亚洲精品国产一区二区精华液| 99国产极品粉嫩在线观看| 亚洲av国产av综合av卡| 天天躁日日躁夜夜躁夜夜| 91麻豆精品激情在线观看国产 | 老汉色∧v一级毛片| 黄频高清免费视频| 久久 成人 亚洲| 久久久久精品人妻al黑| 成人影院久久| 久热这里只有精品99| 国产97色在线日韩免费| 国产亚洲av高清不卡| 91精品三级在线观看| 精品午夜福利视频在线观看一区 | 自线自在国产av| 久久精品国产亚洲av高清一级| 久久久久视频综合| 黄色视频在线播放观看不卡| 欧美午夜高清在线| 99精品久久久久人妻精品| 亚洲av片天天在线观看| 午夜福利欧美成人| 男人操女人黄网站| 91老司机精品| 新久久久久国产一级毛片| 人成视频在线观看免费观看| 国产日韩一区二区三区精品不卡| 亚洲伊人久久精品综合| www.熟女人妻精品国产| 国产男靠女视频免费网站| 国内毛片毛片毛片毛片毛片| 熟女少妇亚洲综合色aaa.| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 国产精品秋霞免费鲁丝片| 亚洲自偷自拍图片 自拍| 久久精品国产亚洲av香蕉五月 | 岛国毛片在线播放| 亚洲精品乱久久久久久| 国产精品亚洲一级av第二区| 香蕉丝袜av| 亚洲精品av麻豆狂野| 欧美午夜高清在线| 免费av中文字幕在线| 国产成人系列免费观看| 女警被强在线播放| 国产高清视频在线播放一区| 成年人黄色毛片网站| 日韩一区二区三区影片| 大型av网站在线播放| 国产亚洲精品一区二区www | 亚洲中文av在线| 青青草视频在线视频观看| 99国产极品粉嫩在线观看| 国产av又大| 一区二区三区乱码不卡18| 亚洲精品美女久久av网站| 欧美一级毛片孕妇| 一个人免费看片子| 三级毛片av免费| 国产aⅴ精品一区二区三区波| 999精品在线视频| 久久久久国内视频| 久久久久久久大尺度免费视频| 亚洲免费av在线视频| 欧美黑人精品巨大| 一级毛片女人18水好多| 丝袜美腿诱惑在线| 国产精品免费视频内射| 国产极品粉嫩免费观看在线| 免费少妇av软件| 欧美精品av麻豆av| 成人免费观看视频高清| 久久久久网色| 成人手机av| 乱人伦中国视频| 亚洲国产中文字幕在线视频| 亚洲男人天堂网一区| 精品卡一卡二卡四卡免费| 天堂中文最新版在线下载| 一边摸一边做爽爽视频免费| 国产极品粉嫩免费观看在线| 啦啦啦免费观看视频1| 最近最新中文字幕大全电影3 | 亚洲性夜色夜夜综合| 757午夜福利合集在线观看| 高清毛片免费观看视频网站 | 亚洲情色 制服丝袜| 国产成人精品在线电影| 亚洲色图av天堂| 一区在线观看完整版| 天天躁狠狠躁夜夜躁狠狠躁| 国产aⅴ精品一区二区三区波| 90打野战视频偷拍视频| 9热在线视频观看99| 69精品国产乱码久久久| 亚洲成人免费av在线播放| 欧美精品av麻豆av| av有码第一页| 97人妻天天添夜夜摸| 免费人妻精品一区二区三区视频| 精品人妻1区二区| 欧美中文综合在线视频| 久久久久精品人妻al黑| 欧美精品啪啪一区二区三区| 制服诱惑二区| 丁香欧美五月| 欧美大码av| 成人亚洲精品一区在线观看| 亚洲一区中文字幕在线| 别揉我奶头~嗯~啊~动态视频| 丁香六月天网| 国产在线精品亚洲第一网站| 俄罗斯特黄特色一大片| 欧美成人午夜精品| 一区二区av电影网| 国产av一区二区精品久久| 美女福利国产在线| 国产精品 国内视频| 午夜精品国产一区二区电影| 午夜久久久在线观看| 19禁男女啪啪无遮挡网站| 在线十欧美十亚洲十日本专区| 亚洲五月婷婷丁香| 男女午夜视频在线观看| 国产人伦9x9x在线观看| 黑丝袜美女国产一区| 国产视频一区二区在线看| 无人区码免费观看不卡 | 看免费av毛片| 免费观看av网站的网址| 久久婷婷成人综合色麻豆| 日韩中文字幕视频在线看片| 深夜精品福利| 欧美国产精品va在线观看不卡| 成年人黄色毛片网站| 亚洲第一av免费看| 久久99热这里只频精品6学生| 高清黄色对白视频在线免费看| 成人av一区二区三区在线看| 久久久久国产一级毛片高清牌| 亚洲精品美女久久久久99蜜臀| 亚洲va日本ⅴa欧美va伊人久久| 一区在线观看完整版| 男女无遮挡免费网站观看| 国精品久久久久久国模美| 亚洲黑人精品在线| 日韩制服丝袜自拍偷拍| 久久午夜综合久久蜜桃| 中文字幕人妻丝袜制服| 亚洲精品中文字幕在线视频| 高清av免费在线| 五月开心婷婷网| 国产精品99久久99久久久不卡| 精品一区二区三区四区五区乱码| 蜜桃在线观看..| 香蕉久久夜色| 久久午夜亚洲精品久久| 在线观看人妻少妇| 我的亚洲天堂| 最近最新中文字幕大全免费视频| 久久久久久久大尺度免费视频| 欧美国产精品一级二级三级| 国产精品.久久久| 国产亚洲av高清不卡| 国产97色在线日韩免费| 国产欧美日韩精品亚洲av| 曰老女人黄片| 极品人妻少妇av视频| 美女高潮喷水抽搐中文字幕| 欧美国产精品va在线观看不卡| 免费在线观看黄色视频的| 久久精品国产亚洲av香蕉五月 | 日本五十路高清| 亚洲情色 制服丝袜| 大型av网站在线播放| 日韩制服丝袜自拍偷拍| 丰满迷人的少妇在线观看| 精品免费久久久久久久清纯 | 日韩熟女老妇一区二区性免费视频| 久热爱精品视频在线9| 成人国语在线视频| 老汉色av国产亚洲站长工具| 欧美日韩av久久| 性少妇av在线| 久久久久精品人妻al黑| 欧美激情久久久久久爽电影 | tube8黄色片| 精品少妇久久久久久888优播| 啪啪无遮挡十八禁网站|