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

    黃土高原不同植被覆蓋土地利用等對(duì)蒸散發(fā)量的影響

    2021-12-17 00:46李燕張凱瑞
    人民黃河 2021年12期
    關(guān)鍵詞:黃土高原土地利用

    李燕 張凱瑞

    摘 要:為實(shí)現(xiàn)對(duì)黃土高原蒸散發(fā)量的動(dòng)態(tài)連續(xù)監(jiān)測(cè),基于蒸散發(fā)量(ET)、歸一化植被指數(shù)(NDVI)、土地利用、氣溫和降水量數(shù)據(jù),采用SEN趨勢(shì)分析、Hurst指數(shù)及偏相關(guān)分析法,探究了黃土高原2001—2020年ET時(shí)空分布特征、變化趨勢(shì)、未來持續(xù)狀態(tài)及降水、氣溫、不同植被覆蓋度、不同坡度、不同土地利用下的ET變化規(guī)律。結(jié)果表明:①黃土高原ET空間分布呈從北向南增大的趨勢(shì),高蒸發(fā)區(qū)集中在渭河平原、六盤山和祁連山區(qū)域,ET隨坡度、植被覆蓋度的增大而增大;②占黃土高原總面積92.02%的區(qū)域ET表現(xiàn)為顯著增大趨勢(shì)(顯著性水平α<0.01),增速為2~5 mm/a,草地和農(nóng)用地的ET增速較大,水體、濕地、建筑用地的ET增速較小;③ET與降水和氣溫均以正相關(guān)關(guān)系為主。

    關(guān)鍵詞:蒸散發(fā);植被覆蓋;土地利用;黃土高原

    中圖分類號(hào):S161.4;TP79

    文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2021.12.014

    引用格式:李燕,張凱瑞.黃土高原不同植被覆蓋/土地利用等對(duì)蒸散發(fā)量的影響[J].人民黃河,2021,43(12):68-73.

    Abstract: In order to achieve the dynamic and continuous monitoring of evapotranspiration in the Loess Plateau, based on evapotranspiration (ET), normalized vegetation index (NDVI), land use data, temperature and precipitation data, this paper adopted SEN trend analysis, R/S analysis and partial correlation analysis and studied the temporal and spatial distribution characteristics, change trend, future continuous state, precipitation and air quality of ET under different vegetation coverage, slope and land use in the Loess Plateau from 2001 to 2020. The results show that: a) the spatial distribution of ET in the Loess Plateau shows an increasing trend from north to south in the past 20 years. b) ET shows a significant increasing trend in a large area during 20 years, and the increasing rate is concentrated in 2-5 mm/a, accounting for 92.02% of the total area of the Loess Plateau. The ETs of grassland and agricultural land increase more, while the ETs of water, wetland and construction land increase less. c) In the future, ET will decrease in most areas of the Loess Plateau, and the persistence in most pixels of Inner Mongolia, Gansu and Qinghai is stronger than that in other areas. d) The relationship between ET and precipitation/air temperature is mainly positive correlation.

    Key words: evapotranspiration; vegetation cover; land use; Loess Plateau

    蒸散發(fā)主要由土壤蒸發(fā)和植物蒸騰組成,將部分地表水分輸送到大氣中[1-3],是干旱與半干旱區(qū)水量消耗的主要途徑[4-7]。黃土高原地處半濕潤(rùn)、半干旱、干旱區(qū)過渡地帶,受西風(fēng)環(huán)流、高原季風(fēng)和東亞季風(fēng)環(huán)流共同影響,干旱事件頻繁發(fā)生[8-12]。因此,充分認(rèn)識(shí)蒸散發(fā)量(ET)和植被覆蓋/土地利用的時(shí)空變化規(guī)律及其響應(yīng)關(guān)系,不僅有利于理解ET對(duì)水循環(huán)的影響機(jī)制,也可為黃土高原水源涵養(yǎng)和水資源合理開發(fā)利用提供理論參考。

    近年來,專家們從不同角度出發(fā)不斷對(duì)蒸散發(fā)開展研究,提出很多種蒸散發(fā)量的計(jì)算方法[13-16]。如Bowen[13]提出了波文比-能量平衡法、Wilm等[14]提出了空氣動(dòng)力學(xué)方法、Penman[15]提出了Penman算法,這些算法雖然能夠獲取較高精度的數(shù)據(jù),但區(qū)域代表性差[16-17]。隨著遙感技術(shù)發(fā)展,美國(guó)國(guó)家航空航天局(NASA)于2011年基于Penman-Monteith遙感模型和MODIS數(shù)據(jù)研發(fā)了全球陸地蒸散發(fā)數(shù)據(jù)產(chǎn)品(MOD16),該產(chǎn)品獲取的蒸散發(fā)信息空間覆蓋范圍廣、數(shù)據(jù)分辨率高,精度能夠滿足大尺度的研究需求[18-20]。目前,基于MOD16數(shù)據(jù)研究ET時(shí)空變化及不同植被覆蓋/土地利用下ET的差異性成為熱點(diǎn)之一[21-25]。如董晴晴等[24]基于MOD16研究了2000—2013年渭河流域ET的空間格局,發(fā)現(xiàn)ET呈緩慢增大趨勢(shì),但年際波動(dòng)性較大,且空間異質(zhì)性較強(qiáng);張翔等[25]基于MOD16研究了2008—2017年呼倫貝爾地區(qū)植被的ET,發(fā)現(xiàn)其整體呈增大趨勢(shì),但不同區(qū)域的ET差異明顯,年際波動(dòng)性較大;谷佳賀等[26]基于MOD16、NDVI和不同土地利用數(shù)據(jù)對(duì)黃河流域ET時(shí)空分布及不同植被覆蓋下ET變化規(guī)律研究發(fā)現(xiàn),年均ET呈東南大、西北小的空間分布格局,與不同植被覆蓋和土地利用區(qū)域有很大關(guān)系;Wang等[27]研究發(fā)現(xiàn)低植被覆蓋時(shí),植被變化對(duì)ET產(chǎn)生負(fù)面影響,隨著植被覆蓋度增大,這種負(fù)效應(yīng)逐漸變?yōu)檎?yīng)。目前對(duì)蒸散發(fā)的研究雖較多,但多采用單一時(shí)刻遙感數(shù)據(jù)對(duì)某一地區(qū)ET進(jìn)行研究[16-19],很少探討長(zhǎng)時(shí)間序列、不同植被覆蓋下的ET變化特征和未來可能持續(xù)狀態(tài)。

    筆者基于2001—2020年黃土高原的MODIS ET、NDVI數(shù)據(jù),采用Sen、Hurst指數(shù)和偏相關(guān)方法對(duì)植被覆蓋度/不同土地利用類型的ET時(shí)空分布特征、變化趨勢(shì)及氣候變化對(duì)蒸散發(fā)的影響進(jìn)行分析,以期掌握退耕還林(還草)等生態(tài)工程措施下黃土高原蒸散發(fā)、植被覆蓋的時(shí)空變化特征,為該地區(qū)生態(tài)環(huán)境保護(hù)提供科學(xué)依據(jù)。

    1 研究區(qū)概況與研究方法

    1.1 研究區(qū)概況

    黃土高原地處中國(guó)西北地區(qū),總面積64.9萬km2,包括山西、內(nèi)蒙古、陜西、河南、寧夏、甘肅、青海7個(gè)?。▍^(qū))。降水量少而不均,從西北向東南逐漸增大,年降水量300~800 mm,降水主要集中在7—9月;年均氣溫為3.6~14.3 ℃,冬季嚴(yán)寒、夏季暖熱。黃土質(zhì)地疏松,具有濕陷性,經(jīng)人為擾動(dòng)后,極易發(fā)生風(fēng)力和水力侵蝕[1]。區(qū)域內(nèi)地貌類型多樣,植被覆蓋度呈東南高、西北低的空間分布格局,主要土地類型為草地、林地和農(nóng)用地[9-10](見圖1)。

    1.2 數(shù)據(jù)來源與處理

    (1)NDVI數(shù)據(jù):來源于Google Earth Engine(GEE,https://code.earthengine.google.com/)平臺(tái)提供的2001—2020年MOD13Q1產(chǎn)品的NDVI數(shù)據(jù),空間分辨率為250 m,時(shí)間分辨率為16 d,經(jīng)過幾何校正和大氣校正,是標(biāo)準(zhǔn)3級(jí)產(chǎn)品數(shù)據(jù)。使用Hants方法對(duì)NDVI時(shí)序數(shù)據(jù)進(jìn)行去噪,并將年均NDVI小于0.05的像元排除在外[28]。

    (2)ET數(shù)據(jù):來源于GEE平臺(tái)提供的2001—2020年的MOD16A2產(chǎn)品,空間分辨率為500 m,時(shí)間分辨率為8 d。該數(shù)據(jù)是利用物理意義明確且不受地表溫度影響的Penman-Monteith模型計(jì)算得到的,很多學(xué)者研究發(fā)現(xiàn)MOD16雖然在某些地區(qū)存在高估或低估現(xiàn)象,但其精度能夠滿足研究大尺度的應(yīng)用需求[29-31]。

    (3)土地利用數(shù)據(jù):來源于GEE平臺(tái)提供的MCD12Q1產(chǎn)品,是根據(jù)生物圈計(jì)劃(IGBP)進(jìn)行分類的全球產(chǎn)品,包括17種土地利用類型。高艷妮等[32]曾以野外調(diào)查獲取的草地樣點(diǎn)數(shù)據(jù)對(duì)不同來源的土地利用類型數(shù)據(jù)進(jìn)行精度評(píng)價(jià),發(fā)現(xiàn)MCD12Q1的采樣點(diǎn)驗(yàn)證精度為83.74%,故該分類精度能夠滿足本研究的要求。

    (4)氣溫、降水量數(shù)據(jù):采用中國(guó)氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/data/cdcindex)提供的月數(shù)據(jù),時(shí)間跨度為2001—2020年,得到研究區(qū)及周圍的130個(gè)氣象站(見圖2)的氣溫、降水量數(shù)據(jù),定義12月至次年2月為冬季、3—5月為春季、6—8月為夏季、9—11月為秋季,利用氣象插值軟件ANUSPLINE將季節(jié)性氣溫、降水量數(shù)據(jù)插值為500 m的柵格數(shù)據(jù)集。

    2 結(jié)果與分析

    2.1 ET與NDVI/土地利用空間分布特征分析

    利用均值法計(jì)算得到2001—2020年黃土高原ET與NDVI的多年平均值(見圖3),可以看出:ET空間分布上呈從北向南增大的趨勢(shì),其中高蒸發(fā)區(qū)集中在渭河流域、六盤山和祁連山山脈,而鄂爾多斯高原、寧夏高原和后套平原地帶的蒸散發(fā)較低;NDVI的空間分布與ET空間分布具有相似性,在ET較小的區(qū)域NDVI值較小(植被覆蓋度低),ET較大的區(qū)域NDVI大。結(jié)合圖1、圖2可以發(fā)現(xiàn),ET較大的區(qū)域主要土地利用類型為農(nóng)用地和建筑用地,且坡度主要為0°~2°;ET較小的內(nèi)蒙古、寧夏和甘肅的土地利用類型主要為草地和林地,坡度主要為0°~2°;黃土高原南部地區(qū)ET較大,坡度一般大于2°。因此,黃土高原的ET與NDVI/土地利用類型和坡度具有一定相關(guān)性,隨坡度增大而增大。

    2.2 NDVI/ET空間變化趨勢(shì)分析

    采用Sen和M-K檢驗(yàn)方法得出的2001—2020年黃土高原ET變化速率空間分布見圖4(a)??梢钥闯觯篍T在占黃土高原總面積92.02%的區(qū)域上表現(xiàn)為顯著增大趨勢(shì)(顯著性水平α<0.01),增速為2~5 mm/a,ET增速大于5 mm/a的區(qū)域集中在陜西省洛河、甘肅隴東和祁連山地區(qū),而在鄂爾多斯和太行山南坡的ET增速為2~3 mm/a;ET呈減小趨勢(shì)的區(qū)域零星分布在渭河平原地區(qū),減小速率為0~1 mm/a。從不同土地利用類型、不同坡度的ET變化統(tǒng)計(jì)數(shù)據(jù)發(fā)現(xiàn),草地和農(nóng)用地的ET增速較大,90%以上的面積ET增速在2 mm/a以上,草地ET增速較大的為2°~15°坡度的區(qū)域,增速基本在3 mm/a以上,而農(nóng)用地ET增速較大的為6°~15°坡度的區(qū)域,見圖4(b)。

    基于Sen和M-K檢驗(yàn)方法得出的2001—2020年黃土高原NDVI變化趨勢(shì)空間分布見圖5(a)。NDVI整體呈顯著增大趨勢(shì)(顯著性水平α<0.01),尤其在呂梁山脈、山西高原及甘肅隴東地區(qū)的增長(zhǎng)趨勢(shì)較大,NDVI呈減小趨勢(shì)的區(qū)域零星分布在陜西省的渭河平原地區(qū)。農(nóng)用地的NDVI增長(zhǎng)趨勢(shì)最明顯,尤其在25°以上坡度的區(qū)域。草地和林地及建筑用地的NDVI增長(zhǎng)趨勢(shì)較為接近,2°~6°坡度的區(qū)域NDVI增長(zhǎng)趨勢(shì)最大,見圖5(b)。

    2.3 ET與NDVI的空間響應(yīng)關(guān)系

    利用2001—2020年黃土高原ET與NDVI的長(zhǎng)時(shí)間序列數(shù)據(jù),計(jì)算得到二者的相關(guān)系數(shù)(見圖6)。占總面積93.89%的區(qū)域ET與NDVI成正相關(guān)關(guān)系,相關(guān)系數(shù)為0.2~0.8;相關(guān)性較好的區(qū)域集中在呂梁山脈和六盤山地帶,相關(guān)系數(shù)在0.6以上。ET與NDVI成負(fù)相關(guān)關(guān)系的區(qū)域集中在祁連山、渭河平原、后套平原和太行山脈南坡,相關(guān)系數(shù)為-0.4~0。草地的ET與NDVI的相關(guān)系數(shù)較其他幾種土地利用類型的相關(guān)系數(shù)大,坡度為2°~6°的草地相關(guān)系數(shù)在0.53以上,坡度為0°~2°的草地相關(guān)系數(shù)為0.47(坡度在25°以上的草地ET與NDVI的相關(guān)系數(shù)小,為0.18);其次為建筑用地,ET與NDVI的相關(guān)系數(shù)隨著坡度的減小而減小,坡度為15°~25°的建筑用地ET與NDVI相關(guān)系數(shù)為0.48;農(nóng)用地的ET與NDVI相關(guān)系數(shù)在不同坡度上的差異較小,相關(guān)系數(shù)均在0.2以上。

    2.4 ET與氣溫、降水的空間響應(yīng)關(guān)系

    在控制氣溫的條件下ET與降水的相關(guān)關(guān)系為正、負(fù)相關(guān)并存,空間分布上由北向南相關(guān)系數(shù)逐漸減?。ㄒ妶D7(a))。占黃土高原總面積90.41%的區(qū)域ET與降水以正相關(guān)為主,相關(guān)系數(shù)為0~0.8,相關(guān)系數(shù)在0.8以上的區(qū)域集中在毛烏素沙地;ET與降水成負(fù)相關(guān)關(guān)系的區(qū)域集中在海拔較高的祁連山、太行山脈、后套平原地區(qū)。在控制降水的條件下,ET與氣溫為正、負(fù)相關(guān)關(guān)系并存,空間分布上表現(xiàn)為中部相關(guān)系數(shù)大、南北兩邊相關(guān)系數(shù)?。ㄒ妶D7(b))。

    毛烏素沙地的ET與氣溫以正相關(guān)為主,相關(guān)系數(shù)在0.6以上;黃土高原南部地區(qū)氣溫與ET相關(guān)系數(shù)為0~0.4;ET與氣溫成負(fù)相關(guān)關(guān)系的地區(qū)集中在鄂爾多斯高原東北部、山西高原,相關(guān)系數(shù)在-0.2以下。

    3 討 論

    本文通過研究2001—2020年黃土高原不同土地利用類型和不同植被覆蓋下的蒸散發(fā)變化趨勢(shì),發(fā)現(xiàn)其蒸散發(fā)量表現(xiàn)為顯著增大趨勢(shì),其中2003年和2016年氣溫異常大于其他年份,導(dǎo)致ET也異常大于其他年份。ET空間上表現(xiàn)為北小南大、西小東大的態(tài)勢(shì),主要原因是西部和北部地勢(shì)高、氣溫低、降水量小、地表干燥及高寒地區(qū)植被生長(zhǎng)季較短,蒸散發(fā)主要靠土壤淺層水分的微弱蒸發(fā),因此蒸散發(fā)量小;黃土高原南部和東部地區(qū)夏季降水量較大,太陽(yáng)輻射強(qiáng),且該地區(qū)的林地、草地、農(nóng)作物等覆蓋度高,植被蒸騰和土壤水分蒸發(fā)作用較強(qiáng),因此蒸散發(fā)量較大。

    2001—2020年黃土高原的洛河流域、甘肅隴東和祁連山周圍蒸散發(fā)量呈顯著增大趨勢(shì)(顯著性水平α<0.01),原因一方面可能與近20 a黃土高原土地利用類型發(fā)生顯著變化及全球氣候變化有關(guān)。有研究[6,36-38]表明,不同土地利用類型對(duì)ET的貢獻(xiàn)率不同,沼澤地主要覆蓋淺水及水生植物,因此沼澤地水分的直接蒸發(fā)比草地、林地等的蒸騰作用強(qiáng)烈;而林地由于林木根系可以吸收土壤深層水分用于蒸騰,因此其ET比其他植被類型大。邵薇薇等[37]研究表明,陸面植被覆蓋率增大導(dǎo)致蒸散發(fā)量增大的同時(shí),會(huì)使陸面水熱通量發(fā)生變化;Piao等[38]研究表明,溫暖的春季促使植被物候提前及植被生長(zhǎng)加速,影響地表蒸散發(fā)量的變化,特別是林地面積的增大加大了蒸散發(fā)量。ET增大的另一方面原因可能是黃土高原東南部降水量嚴(yán)重受到東亞夏季風(fēng)的制約,自20世紀(jì)70年代以來,東亞夏季風(fēng)強(qiáng)度總體表現(xiàn)為顯著減弱趨勢(shì),氣候變暖導(dǎo)致植被蒸騰作用加快,因此蒸散發(fā)增大。而本研究也發(fā)現(xiàn)黃土高原ET與氣溫和降水呈顯著正相關(guān)關(guān)系,其中降水對(duì)ET的影響范圍較廣、影響程度較大,空間分布上存在明顯的由北向南ET與降水相關(guān)系數(shù)逐漸減小的趨勢(shì),尤其在毛烏素沙地。

    4 結(jié) 論

    本文基于NDVI、ET及氣溫、降水量數(shù)據(jù),分析了2001—2020年黃土高原ET與NDVI、不同土地利用類型、氣溫、降水的關(guān)系、時(shí)空格局和變化趨勢(shì)及氣候因子對(duì)ET的影響,得出結(jié)論如下。

    (1)黃土高原ET空間分布呈從北向南增大的趨勢(shì),且隨坡度的增大而增大。

    (2)占黃土高原總面積92.02%的區(qū)域ET表現(xiàn)為顯著增大趨勢(shì)(顯著性水平α<0.01),增速為2~5 mm/a;草地和農(nóng)用地的ET增速較大,林地ET增速隨坡度的增大而減小,建筑用地的ET隨坡度增大而增大;ET隨植被覆蓋度的增大而增大。

    (3)占黃土高原總面積90.41%的區(qū)域ET與降水成正相關(guān)關(guān)系,相關(guān)系數(shù)為0~0.8。ET與氣溫以正相關(guān)關(guān)系為主,空間分布表現(xiàn)為中部相關(guān)系數(shù)大、南北兩邊相關(guān)系數(shù)小。

    參考文獻(xiàn):

    [1] MONTEITH J L. Evaporation and Environment[J]. Symposia of the Society for Experimental Biology,1965,19:205-234.

    [2] XIONG Y J, ZHAO S H, TIAN F, et al. An Evapotranspiration Product for Arid Regions Based on the Three-Temperature Model and Thermal Remote Sensing[J]. Journal of Hydrology,2015,530:392-404.

    [3] 周志鵬,孫文義,穆興民,等.2001—2017年黃土高原實(shí)際蒸散發(fā)的時(shí)空格局[J].人民黃河,2019,41(6):76-80,84.

    [4] NISHIDA K, NEMANI R R, GLASSY J M, et al. Development of an Evapotranspiration Index from Aqua/MODIS for Monitoring Surface Moisture Status[J]. IEEE Transactions on Geoscience & Remote Sensing,2003,41(2):493-501.

    [5] DONG Q, ZHAN C, WANG H, et al. A Review on Evapotranspiration Data Assimilation Based on Hydrological Models[J]. Journal of Geographical Sciences,2016,26(2):230-242.

    [6] 張戴煒,王春穎,韓宇平,等.基于溫度的ET_0模型改進(jìn)及在海河流域的適用性[J].人民黃河,2021,43(3):155-160.

    [7] 鄧興耀,劉洋,劉志輝,等.中國(guó)西北干旱區(qū)蒸散發(fā)時(shí)空動(dòng)態(tài)特征[J].生態(tài)學(xué)報(bào),2017,37(9):2994-3008.

    [8] 周忠學(xué).陜北黃土高原土地利用變化與社會(huì)經(jīng)濟(jì)發(fā)展關(guān)系及效應(yīng)評(píng)價(jià)[D].西安:陜西師范大學(xué),2007:21-37.

    [9] 孫藝杰,劉憲鋒,任志遠(yuǎn),等.1960—2016年黃土高原多尺度干旱特征及影響因素[J].地理研究,2020,39(4):591-601.

    [10] 楊濤,陸桂華,李會(huì)會(huì),等.氣候變化下水文極端事件變化預(yù)測(cè)研究進(jìn)展[J].水科學(xué)進(jìn)展,2011,22(2):279-286.

    [11] 秦大河,丁一匯,王紹武,等.中國(guó)西部生態(tài)環(huán)境變化與對(duì)策建議[J].地球科學(xué)進(jìn)展,2002,17(3):314-319.

    [12] 姚玉璧,張強(qiáng),李耀輝,等.干旱災(zāi)害風(fēng)險(xiǎn)評(píng)估技術(shù)及其科學(xué)問題與展望[J].資源科學(xué),2013,35(9):1884-1897.

    [13] BOWEN I S. The Ration of Heat Losses by Conduction and Evaporation from Any Water Surface[J]. Physical Review,1926,27(6):779-798.

    [14] WILM H G,THORNTHWAITE C W. Report of the Committee on Transpiration and Evaporation[J]. Transactions of the American Sections of the American Geophysical Union,1944,25(5):683-693.

    [15] PENMAN H L. Natural Evaporation from Open Water, Bare Soil and Grass[J]. Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences,1948,193(1032):120-145.

    [16] 韓松俊,張寶忠.基于Penman方法和互補(bǔ)原理的蒸散發(fā)研究歷程與展望[J].水利學(xué)報(bào),2018,49(9):1158-1168.

    [17] KONG J, HU Y, YANG L, et al. Estimation of Evapotranspiration for the Blown-Sand Region in the Ordos Basin Based on the SEBAL Model[J]. International Journal of Remote Sensing,2019,40(5-6):1945-1965.

    [18] MU Q, ZHAO M, STEVEN W. Improvements to a MODIS Global Terrestrial Evapotranspiration Algorithm[J]. Remote Sensing of Environment,2011,115(8):1781-1800.

    [19] MU Q, HEINSCH, F A, ZHAO M, et al. Development of a Global Evapotranspiration Algorithm Based on MODIS and Global Meteorology Data[J]. Remote Sensing of Environment,2007,111(4):519-536.

    [20] 徐麗娟,張文鴿,李皓冰,等.鄂爾多斯市潛在蒸散發(fā)時(shí)空變化分析[J].人民黃河,2018,40(9):53-57.

    [21] 喻元,白建軍,王建博,等.基于MOD16的關(guān)中地區(qū)實(shí)際蒸散發(fā)時(shí)空特征分析[J].干旱地區(qū)農(nóng)業(yè)研究,2015,33(3):245-253.

    [22] 何慧娟,卓靜,董金芳,等.基于MOD16監(jiān)測(cè)陜西省地表蒸散變化[J].干旱區(qū)地理,2015,38(5):960-967.

    [23] 溫媛媛,趙軍,王炎強(qiáng),等.基于MOD16的山西省地表蒸散發(fā)時(shí)空變化特征分析[J].地理科學(xué)進(jìn)展,2020,39(2):255-264.

    [24] 董晴晴,占車生,王會(huì)肖,等.2000年以來的渭河流域?qū)嶋H蒸散發(fā)時(shí)空格局分析[J].干旱區(qū)地理,2016,39(2):327-335.

    [25] 張翔,朱曉昱,沈貝貝,等.呼倫貝爾植被蒸散發(fā)時(shí)空變化研究[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃,2020,41(4):308-316.

    [26] 谷佳賀,薛華柱,董國(guó)濤,等.黃河流域NDVI/土地利用對(duì)蒸散發(fā)時(shí)空變化的影響[J].干旱區(qū)地理,2021,44(1):158-167.

    [27] WANG Y, LIU Y B, JIN J X. Contrast Effects of Vegetation Cover Change on Evapotranspiration During a Revegetation Period in the Poyang Lake Basin, China[J].Forests,2018,9(4):217-230.

    [28] 秦格霞,吳靜,李純斌,等.中國(guó)北方草地植被物候變化及其對(duì)氣候變化的響應(yīng)[J].應(yīng)用生態(tài)學(xué)報(bào),2019,30(12):4099-4107.

    [29] 姜艷陽(yáng),王文,周正昊.MODIS MOD16蒸散發(fā)產(chǎn)品在中國(guó)流域的質(zhì)量評(píng)估[J].自然資源學(xué)報(bào),2017,32(3):517-528.

    [30] 李佳,辛?xí)灾?,彭志晴,?地表蒸散發(fā)遙感產(chǎn)品比較與分析[J].遙感技術(shù)與應(yīng)用,2021,36(1):103-120.

    [31] 晉銳,李新,馬明國(guó),等.陸地定量遙感產(chǎn)品的真實(shí)性檢驗(yàn)關(guān)鍵技術(shù)與試驗(yàn)驗(yàn)證[J].地球科學(xué)進(jìn)展,2017,32(6):630-642.

    [32] 高艷妮,郭艷芳,王維,等.不同土地利用/覆蓋數(shù)據(jù)在中國(guó)北方草地的精度評(píng)價(jià)[J].生態(tài)學(xué)雜志,2019,38(1):283-293.

    [33] 汪攀,劉毅敏.Sen’s 斜率估計(jì)與Mann-Kendall法在設(shè)備運(yùn)行趨勢(shì)分析中的應(yīng)用[J].武漢科技大學(xué)學(xué)報(bào),2014,37(6):454-457,472.

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

    [35] 徐建華.計(jì)量地理學(xué)[M].北京:高等教育出版社,2014:109-200.

    [36] 胡光印,金會(huì)軍,董治寶,等.黃河源區(qū)土地利用/覆蓋變化(LUCC)研究[J].冰川凍土,2014,36(3):573-581.

    [37] 邵薇薇,楊大文,孫福寶,等.黃土高原地區(qū)植被與水循環(huán)的關(guān)系[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,49(12):1958-1962.

    [38] PIAO S, FANG J, ZHOU L, et al. Variations in Satellite-Derived Phenology in Chinas Temperate Vegetation[J].Global Change Biology, 2006,12(4):672-685.

    【責(zé)任編輯 張華興】

    猜你喜歡
    黃土高原土地利用
    共享社會(huì)經(jīng)濟(jì)路徑在土地利用、能源與碳排放研究的應(yīng)用
    我國(guó)黃土高原區(qū)域地質(zhì)地貌特征及風(fēng)土人情
    黃土高坡上的江南小景
    黃河哺育了黃土高原
    西咸新區(qū)國(guó)土資源“一張圖”平臺(tái)建設(shè)初探
    北京市土地利用現(xiàn)狀分析
    景頗族土地利用方式及其特征探析
    宜君縣黃土高原流域治理實(shí)施的必要性與建設(shè)條件
    小石獅——黃土高原農(nóng)家的“守護(hù)神”
    亚洲成人av在线免费| 久久久久视频综合| 久久韩国三级中文字幕| a级毛色黄片| 涩涩av久久男人的天堂| 国产男女超爽视频在线观看| 国产伦精品一区二区三区视频9| 色视频在线一区二区三区| 一本一本综合久久| 国产爽快片一区二区三区| 秋霞伦理黄片| 男的添女的下面高潮视频| 成人国语在线视频| 丰满少妇做爰视频| 极品少妇高潮喷水抽搐| 久久久久国产网址| 日韩中字成人| 少妇猛男粗大的猛烈进出视频| 天堂俺去俺来也www色官网| 夜夜爽夜夜爽视频| 午夜福利网站1000一区二区三区| 桃花免费在线播放| 免费观看的影片在线观看| 九九爱精品视频在线观看| 亚洲在久久综合| 亚洲精品久久久久久婷婷小说| 国产欧美另类精品又又久久亚洲欧美| 王馨瑶露胸无遮挡在线观看| 亚洲av在线观看美女高潮| 亚洲精品美女久久av网站| 免费人妻精品一区二区三区视频| 成人漫画全彩无遮挡| 欧美精品亚洲一区二区| 亚洲人成网站在线观看播放| 最新的欧美精品一区二区| 2018国产大陆天天弄谢| 九色亚洲精品在线播放| 国产又色又爽无遮挡免| 大香蕉久久成人网| 伦理电影免费视频| 日韩 亚洲 欧美在线| 国产日韩欧美在线精品| 欧美精品一区二区大全| 91久久精品国产一区二区三区| 国产片特级美女逼逼视频| tube8黄色片| 人体艺术视频欧美日本| 99热网站在线观看| av不卡在线播放| 中文字幕人妻熟人妻熟丝袜美| 99re6热这里在线精品视频| 亚洲精品日韩在线中文字幕| 青青草视频在线视频观看| 日韩一区二区三区影片| 亚洲性久久影院| 一区二区三区精品91| 激情五月婷婷亚洲| 天堂俺去俺来也www色官网| 97超碰精品成人国产| 亚洲精品亚洲一区二区| 99九九线精品视频在线观看视频| 国产亚洲一区二区精品| 黄色配什么色好看| 九九久久精品国产亚洲av麻豆| 少妇猛男粗大的猛烈进出视频| 18禁动态无遮挡网站| 日本黄色日本黄色录像| 精品国产一区二区久久| 久久国内精品自在自线图片| 久久久a久久爽久久v久久| 夜夜看夜夜爽夜夜摸| 99久久中文字幕三级久久日本| 99久久综合免费| 内地一区二区视频在线| 中文字幕人妻丝袜制服| √禁漫天堂资源中文www| 中文乱码字字幕精品一区二区三区| 天天影视国产精品| 这个男人来自地球电影免费观看 | 伦精品一区二区三区| 最新的欧美精品一区二区| 亚洲欧洲精品一区二区精品久久久 | 最近2019中文字幕mv第一页| 搡女人真爽免费视频火全软件| 国产高清有码在线观看视频| 国产乱人偷精品视频| 国产av精品麻豆| 一级毛片 在线播放| 人妻人人澡人人爽人人| 免费人成在线观看视频色| 国产女主播在线喷水免费视频网站| 777米奇影视久久| 日韩亚洲欧美综合| 亚洲精品久久成人aⅴ小说 | 亚洲四区av| av在线老鸭窝| 亚洲激情五月婷婷啪啪| 天美传媒精品一区二区| 乱人伦中国视频| 最近中文字幕2019免费版| 久久久久久久久大av| 欧美日韩国产mv在线观看视频| 精品亚洲成国产av| 伊人亚洲综合成人网| 精品久久蜜臀av无| a级毛片黄视频| 一级毛片aaaaaa免费看小| 交换朋友夫妻互换小说| 18禁在线播放成人免费| 亚洲少妇的诱惑av| 亚洲欧美精品自产自拍| av天堂久久9| 韩国高清视频一区二区三区| 亚洲av二区三区四区| 女人精品久久久久毛片| 51国产日韩欧美| 18+在线观看网站| 一区二区三区精品91| 日韩成人av中文字幕在线观看| 亚洲国产欧美日韩在线播放| 成人免费观看视频高清| 国产国拍精品亚洲av在线观看| 日韩 亚洲 欧美在线| 尾随美女入室| 日本色播在线视频| 午夜老司机福利剧场| 在线观看一区二区三区激情| 我的老师免费观看完整版| 精品久久国产蜜桃| 国产一级毛片在线| 欧美日韩一区二区视频在线观看视频在线| 又粗又硬又长又爽又黄的视频| 国产 精品1| 成人午夜精彩视频在线观看| 欧美bdsm另类| 精品国产国语对白av| 免费黄频网站在线观看国产| 久久久亚洲精品成人影院| 久久久久人妻精品一区果冻| 欧美亚洲日本最大视频资源| 女人久久www免费人成看片| 精品久久久久久电影网| 99热6这里只有精品| 精品少妇黑人巨大在线播放| 一本—道久久a久久精品蜜桃钙片| 在现免费观看毛片| 在线播放无遮挡| 一区在线观看完整版| videossex国产| 亚洲国产毛片av蜜桃av| 国产午夜精品久久久久久一区二区三区| 久久97久久精品| 涩涩av久久男人的天堂| 亚洲精品日韩在线中文字幕| 欧美日韩国产mv在线观看视频| 亚洲三级黄色毛片| 国产一区二区三区av在线| 黄色毛片三级朝国网站| 极品少妇高潮喷水抽搐| 最新的欧美精品一区二区| av在线播放精品| 高清在线视频一区二区三区| 永久免费av网站大全| 日韩强制内射视频| 精品亚洲乱码少妇综合久久| 视频在线观看一区二区三区| 亚洲高清免费不卡视频| 日韩在线高清观看一区二区三区| 成人国产av品久久久| 国产熟女欧美一区二区| 国产亚洲午夜精品一区二区久久| 欧美亚洲日本最大视频资源| 边亲边吃奶的免费视频| 日日摸夜夜添夜夜爱| 国产亚洲av片在线观看秒播厂| 美女大奶头黄色视频| 亚洲av电影在线观看一区二区三区| 在线亚洲精品国产二区图片欧美 | 亚洲人与动物交配视频| 久久精品人人爽人人爽视色| 久久午夜综合久久蜜桃| 伊人久久国产一区二区| 国产精品偷伦视频观看了| 久久韩国三级中文字幕| 久久婷婷青草| 欧美bdsm另类| 亚洲精品第二区| 伦理电影大哥的女人| 亚洲精品色激情综合| 少妇熟女欧美另类| 99视频精品全部免费 在线| 欧美国产精品一级二级三级| 99热这里只有精品一区| 91午夜精品亚洲一区二区三区| 久久97久久精品| 国产成人精品一,二区| 欧美亚洲 丝袜 人妻 在线| 国产成人免费无遮挡视频| 精品一区二区三卡| 精品国产一区二区久久| 国产69精品久久久久777片| 亚洲av二区三区四区| 精品久久久久久久久av| 女的被弄到高潮叫床怎么办| 国产成人精品无人区| 天美传媒精品一区二区| 国产精品麻豆人妻色哟哟久久| 汤姆久久久久久久影院中文字幕| 日本欧美国产在线视频| 伊人久久精品亚洲午夜| 国产精品三级大全| 中文字幕人妻熟人妻熟丝袜美| 这个男人来自地球电影免费观看 | 永久免费av网站大全| 精品少妇久久久久久888优播| 精品亚洲成国产av| 欧美bdsm另类| a级毛色黄片| 国产av码专区亚洲av| 国产精品一区二区在线不卡| 又黄又爽又刺激的免费视频.| 日产精品乱码卡一卡2卡三| 在线观看免费高清a一片| 欧美97在线视频| 欧美日韩精品成人综合77777| 日本欧美国产在线视频| 久久精品人人爽人人爽视色| av国产精品久久久久影院| 亚洲欧美色中文字幕在线| 日本与韩国留学比较| 亚洲精品国产av成人精品| 日韩在线高清观看一区二区三区| 在线观看三级黄色| 成年av动漫网址| 母亲3免费完整高清在线观看 | 一区二区av电影网| 丰满饥渴人妻一区二区三| 99国产综合亚洲精品| 在线观看国产h片| 午夜福利视频在线观看免费| 3wmmmm亚洲av在线观看| 少妇熟女欧美另类| 边亲边吃奶的免费视频| 一区二区三区免费毛片| 欧美激情 高清一区二区三区| 一级毛片 在线播放| 另类精品久久| 满18在线观看网站| 最后的刺客免费高清国语| 国产又色又爽无遮挡免| 国产精品国产三级国产av玫瑰| 久久人人爽人人片av| 亚洲人成77777在线视频| 国产视频首页在线观看| 亚洲美女搞黄在线观看| 欧美日韩在线观看h| 一本大道久久a久久精品| 国产探花极品一区二区| 国产片内射在线| 精品亚洲乱码少妇综合久久| 久久精品国产自在天天线| 国产精品一区www在线观看| 波野结衣二区三区在线| 欧美性感艳星| 久热这里只有精品99| 亚洲伊人久久精品综合| 日本黄色片子视频| 22中文网久久字幕| 亚洲av免费高清在线观看| 岛国毛片在线播放| 青春草亚洲视频在线观看| 久久久久久久久久久丰满| 汤姆久久久久久久影院中文字幕| 亚洲国产精品成人久久小说| 插逼视频在线观看| 国产免费一区二区三区四区乱码| 女人久久www免费人成看片| 成人18禁高潮啪啪吃奶动态图 | 草草在线视频免费看| 久久99热6这里只有精品| 亚洲精品aⅴ在线观看| 99久久综合免费| 99国产精品免费福利视频| 一级毛片电影观看| 99久久精品一区二区三区| 日韩大片免费观看网站| av在线app专区| 男人操女人黄网站| 婷婷色av中文字幕| 日产精品乱码卡一卡2卡三| 国产老妇伦熟女老妇高清| 欧美日韩综合久久久久久| 久久97久久精品| 熟女av电影| 国产精品三级大全| av专区在线播放| 最近最新中文字幕免费大全7| 国产欧美日韩一区二区三区在线 | av黄色大香蕉| 精品一区二区免费观看| 高清毛片免费看| 免费观看a级毛片全部| 亚洲欧美日韩另类电影网站| 日韩av在线免费看完整版不卡| 国产亚洲精品久久久com| 久久久国产精品麻豆| 亚洲情色 制服丝袜| 午夜视频国产福利| 欧美少妇被猛烈插入视频| 黄色欧美视频在线观看| 精品午夜福利在线看| 乱人伦中国视频| 欧美bdsm另类| 熟妇人妻不卡中文字幕| 视频区图区小说| 亚洲成人一二三区av| 男女边摸边吃奶| a级毛片黄视频| 肉色欧美久久久久久久蜜桃| 欧美日韩国产mv在线观看视频| 亚洲中文av在线| 久久综合国产亚洲精品| 欧美+日韩+精品| 男女高潮啪啪啪动态图| 久热这里只有精品99| 自拍欧美九色日韩亚洲蝌蚪91| av黄色大香蕉| 成人国产av品久久久| 欧美bdsm另类| 99热这里只有精品一区| 精品国产乱码久久久久久小说| 建设人人有责人人尽责人人享有的| 中文字幕人妻丝袜制服| 少妇猛男粗大的猛烈进出视频| 亚洲欧美色中文字幕在线| 精品久久久久久久久亚洲| xxxhd国产人妻xxx| 熟女av电影| 亚洲欧洲精品一区二区精品久久久 | 韩国高清视频一区二区三区| 国产免费福利视频在线观看| 一级毛片我不卡| 中国美白少妇内射xxxbb| 国产视频首页在线观看| 欧美97在线视频| 久久久久视频综合| 国产深夜福利视频在线观看| 69精品国产乱码久久久| 亚洲五月色婷婷综合| 亚洲精品乱码久久久久久按摩| 日韩欧美精品免费久久| 你懂的网址亚洲精品在线观看| 久久97久久精品| 久久精品国产亚洲av涩爱| √禁漫天堂资源中文www| 美女国产高潮福利片在线看| av.在线天堂| 多毛熟女@视频| 涩涩av久久男人的天堂| 精品国产乱码久久久久久小说| 51国产日韩欧美| 大陆偷拍与自拍| 久久99热这里只频精品6学生| 中国美白少妇内射xxxbb| 亚洲av福利一区| 日日啪夜夜爽| 婷婷色综合大香蕉| 亚洲图色成人| 国产一区二区在线观看日韩| 蜜桃久久精品国产亚洲av| freevideosex欧美| 国产精品一二三区在线看| 亚洲精品视频女| 亚洲国产最新在线播放| 中文字幕人妻丝袜制服| 女性被躁到高潮视频| 欧美日韩综合久久久久久| 性色avwww在线观看| 欧美激情国产日韩精品一区| 免费观看无遮挡的男女| 国产伦理片在线播放av一区| 综合色丁香网| 欧美日本中文国产一区发布| 久久综合国产亚洲精品| 免费大片黄手机在线观看| 国产精品欧美亚洲77777| 亚洲人成网站在线观看播放| 狂野欧美激情性xxxx在线观看| av在线老鸭窝| 九色成人免费人妻av| 精品人妻一区二区三区麻豆| 99久久中文字幕三级久久日本| 中文字幕制服av| 看非洲黑人一级黄片| av国产久精品久网站免费入址| .国产精品久久| 亚洲国产最新在线播放| 男女免费视频国产| 嫩草影院入口| 熟女电影av网| 欧美 亚洲 国产 日韩一| 国产白丝娇喘喷水9色精品| 国产成人a∨麻豆精品| 免费看光身美女| 啦啦啦中文免费视频观看日本| 另类亚洲欧美激情| 日韩成人av中文字幕在线观看| 亚洲人成网站在线观看播放| 免费大片18禁| 国产色婷婷99| 亚洲综合色网址| 日韩欧美一区视频在线观看| 国产成人91sexporn| 亚洲少妇的诱惑av| 国产一区亚洲一区在线观看| 丁香六月天网| 亚洲欧洲日产国产| 国产淫语在线视频| 日韩强制内射视频| 日本黄色日本黄色录像| 最近中文字幕高清免费大全6| 国产在线免费精品| 久久精品久久精品一区二区三区| 国产成人精品福利久久| 日韩成人av中文字幕在线观看| 精品酒店卫生间| 亚洲成人手机| 在线看a的网站| 最近2019中文字幕mv第一页| 制服诱惑二区| 精品亚洲乱码少妇综合久久| 国产国拍精品亚洲av在线观看| 亚洲人成网站在线观看播放| 91在线精品国自产拍蜜月| 精品酒店卫生间| 午夜老司机福利剧场| 久久久久久久久久人人人人人人| 男人爽女人下面视频在线观看| a级毛片在线看网站| 天堂俺去俺来也www色官网| av免费观看日本| 国产深夜福利视频在线观看| 欧美精品国产亚洲| 尾随美女入室| 美女内射精品一级片tv| 亚洲欧美精品自产自拍| 久久精品国产亚洲网站| 国产极品粉嫩免费观看在线 | 日韩av不卡免费在线播放| 国产精品一区二区三区四区免费观看| 亚洲,一卡二卡三卡| 日本黄色日本黄色录像| 精品少妇久久久久久888优播| 波野结衣二区三区在线| 99九九在线精品视频| 狠狠婷婷综合久久久久久88av| 丝袜脚勾引网站| 亚洲av日韩在线播放| 久久精品国产亚洲网站| 少妇丰满av| 特大巨黑吊av在线直播| 亚州av有码| 新久久久久国产一级毛片| 九色成人免费人妻av| 亚洲性久久影院| 国产日韩一区二区三区精品不卡 | 免费播放大片免费观看视频在线观看| 又黄又爽又刺激的免费视频.| 亚洲国产精品专区欧美| 免费人成在线观看视频色| 啦啦啦在线观看免费高清www| 热99久久久久精品小说推荐| 精品久久国产蜜桃| 人妻制服诱惑在线中文字幕| 国产在视频线精品| 永久网站在线| 啦啦啦啦在线视频资源| 国产极品粉嫩免费观看在线 | 99久久中文字幕三级久久日本| 亚洲精品一区蜜桃| 日本猛色少妇xxxxx猛交久久| 男人添女人高潮全过程视频| 97在线视频观看| 欧美日本中文国产一区发布| 国产亚洲精品久久久com| 这个男人来自地球电影免费观看 | 99视频精品全部免费 在线| 国产精品久久久久久精品电影小说| 欧美xxⅹ黑人| 在线观看www视频免费| 成人国产麻豆网| 中文字幕免费在线视频6| 亚洲精品乱久久久久久| 亚洲欧美色中文字幕在线| 制服人妻中文乱码| 成人午夜精彩视频在线观看| 国产免费视频播放在线视频| 国产有黄有色有爽视频| 国产高清有码在线观看视频| 久久毛片免费看一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲精品一二三| 青春草亚洲视频在线观看| 午夜免费观看性视频| 天堂8中文在线网| 欧美+日韩+精品| 亚洲综合精品二区| 久久久久久久久久人人人人人人| 少妇 在线观看| 国产极品粉嫩免费观看在线 | 国产白丝娇喘喷水9色精品| 免费日韩欧美在线观看| 成人黄色视频免费在线看| 在线观看人妻少妇| 国产男人的电影天堂91| 国产一区二区三区综合在线观看 | 性色av一级| 极品人妻少妇av视频| 日日啪夜夜爽| 搡老乐熟女国产| videosex国产| 最近手机中文字幕大全| 伊人亚洲综合成人网| 最近2019中文字幕mv第一页| 中文欧美无线码| 妹子高潮喷水视频| 少妇人妻 视频| 日韩欧美一区视频在线观看| 午夜福利视频精品| 日韩一本色道免费dvd| a级毛色黄片| 国产成人免费观看mmmm| 色5月婷婷丁香| 啦啦啦中文免费视频观看日本| 男女免费视频国产| 国产精品久久久久久久久免| 久久久久精品久久久久真实原创| 日韩一区二区视频免费看| 精品久久久久久久久av| 看十八女毛片水多多多| 亚洲人与动物交配视频| 如日韩欧美国产精品一区二区三区 | 黄色欧美视频在线观看| 国产精品人妻久久久影院| 午夜福利,免费看| 男男h啪啪无遮挡| 国产av码专区亚洲av| 久久久精品区二区三区| 少妇被粗大猛烈的视频| 男人添女人高潮全过程视频| 一二三四中文在线观看免费高清| 伊人亚洲综合成人网| 国产69精品久久久久777片| 亚洲精品乱久久久久久| 精品一品国产午夜福利视频| 日日摸夜夜添夜夜爱| 久久热精品热| 啦啦啦啦在线视频资源| 亚洲av福利一区| 免费少妇av软件| 亚洲精品国产色婷婷电影| 人人妻人人澡人人爽人人夜夜| 亚洲精品成人av观看孕妇| 国产一区二区三区av在线| 欧美一级a爱片免费观看看| 青春草亚洲视频在线观看| 少妇高潮的动态图| 久久久久国产精品人妻一区二区| 成人手机av| 一个人免费看片子| 国产精品三级大全| av视频免费观看在线观看| av国产久精品久网站免费入址| 亚洲第一区二区三区不卡| 女的被弄到高潮叫床怎么办| 亚洲伊人久久精品综合| 亚洲性久久影院| 在线 av 中文字幕| 18在线观看网站| 国产高清三级在线| 国产黄片视频在线免费观看| 久久国产精品大桥未久av| 国产69精品久久久久777片| 在线观看三级黄色| 久久亚洲国产成人精品v| 国内精品宾馆在线| 美女主播在线视频| 亚洲少妇的诱惑av| 天美传媒精品一区二区| 国产永久视频网站| 欧美丝袜亚洲另类| 九色成人免费人妻av| 国产男女超爽视频在线观看| 韩国高清视频一区二区三区| 精品一区二区免费观看| 伦精品一区二区三区| 香蕉精品网在线| 日韩大片免费观看网站| 看免费成人av毛片| 好男人视频免费观看在线| 国产永久视频网站| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧洲国产日韩| 国产毛片在线视频| 老司机亚洲免费影院| 最近中文字幕2019免费版| 香蕉精品网在线| 国产亚洲最大av| 亚洲欧美一区二区三区黑人 | 九色亚洲精品在线播放| 一级毛片 在线播放| 黄色配什么色好看| 人妻系列 视频| 国产国拍精品亚洲av在线观看| 成人亚洲欧美一区二区av|