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

    全天候區(qū)域地表蒸散發(fā)反演
    ——以黑河流域?yàn)槔?/h1>
    2019-05-10 02:53:30廖前瑜段四波韓曉靜
    中國農(nóng)業(yè)信息 2019年1期
    關(guān)鍵詞:全天候黑河植被指數(shù)

    廖前瑜 ,任 超 ,冷 佩 ,段四波,韓曉靜

    (1. 桂林理工大學(xué)測繪地理信息學(xué)院,廣西桂林541004;2. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京100081)

    0 引言

    蒸散發(fā)是地表水分循環(huán)與能量平衡的重要組成分量,反映了大氣與地球表面的水分和能量交換[1]。由于全球氣候變化和人類對水資源的過度利用,導(dǎo)致了可利用水資源短缺和相關(guān)生態(tài)問題。因此,精確估算蒸散發(fā)在時(shí)空尺度上的變化,對農(nóng)田灌溉管理、作物需水量估算、水資源優(yōu)化配置等領(lǐng)域的研究具有重要意義。傳統(tǒng)的蒸散發(fā)估算大多基于氣象觀測的單點(diǎn)計(jì)算,盡管能得到相對準(zhǔn)確的蒸散發(fā)量,但無法實(shí)現(xiàn)非均勻下墊面的蒸散發(fā)研究,而遙感對于區(qū)域蒸散發(fā)估算具有無可替代的優(yōu)勢[2-3]。自20 世紀(jì)80 年代以來,眾多學(xué)者圍繞地表蒸散發(fā)遙感反演開展了大量工作,形成了許多切實(shí)可行的算法。這些算法大致可以分為基于能量平衡的方法[4-7],基于水量平衡的方法[8]以及基于遙感地表溫度—植被指數(shù)特征空間的方法[9-13]。其中,基于遙感地表溫度—植被指數(shù)特征空間的方法近年來得到了最為廣泛的應(yīng)用,主要得益于其簡單的原理和輸入以及它能夠避免用于計(jì)算水和熱傳遞的空氣動(dòng)力學(xué)阻抗和表面阻抗的復(fù)雜參數(shù)化過程。應(yīng)用地表溫度—植被指數(shù)特征空間方法反演地表蒸散發(fā)的一個(gè)重要前提是要確定特征空間的干濕邊。目前,確定地表溫度—植被指數(shù)特征空間干濕邊主要有基于傳統(tǒng)的空間統(tǒng)計(jì)的方法[14-15]和基于能量平衡理論干濕邊計(jì)算方法[13,16]。其中,基于傳統(tǒng)的空間統(tǒng)計(jì)的方法通常具有一定主觀性,且要求研究區(qū)氣象條件一致,導(dǎo)致其在區(qū)域尺度上難以應(yīng)用?;谀芰科胶獾姆椒▌t需要輸入大量的參數(shù),包括一些遙感難以獲取的參數(shù)如風(fēng)速、空氣動(dòng)力學(xué)阻抗及冠層阻抗等,中間參量過多會(huì)導(dǎo)致誤差的累積和傳遞。盡管基于能量平衡估算干濕邊的方法需要大量的輸入?yún)?shù),但因其具有良好的普適性和物理基礎(chǔ)受到了研究學(xué)者們廣泛的關(guān)注。

    最近Leng 等[16]提出了逐像元確定地表溫度—植被指數(shù)特征空間理論干濕邊的思想,避免了傳統(tǒng)特征空間方法在區(qū)域尺度應(yīng)用的局限。由于在有云情況下無法獲取有效的植被指數(shù)和地表溫度參數(shù),Leng 等[16]結(jié)合經(jīng)典的Penman-Monteith 模型,發(fā)展了全天候地表蒸散發(fā)反演方法。然而,該研究僅在站點(diǎn)尺度進(jìn)行了分析。在區(qū)域尺度上,個(gè)別站點(diǎn)觀測的氣象數(shù)據(jù)顯然無法很好地描述氣象因子的空間異質(zhì)性。因此,文章在Leng等[16]研究的基礎(chǔ)上,基于目前常用的衛(wèi)星陸表遙感產(chǎn)品和較高空間分辨率的格網(wǎng)氣象數(shù)據(jù),深入探索區(qū)域尺度全天候蒸散發(fā)的反演。

    1 數(shù)據(jù)

    1.1 研究區(qū)與地面觀測數(shù)據(jù)

    黑河流域(97°E~102°E,38°N~43°N)面積約為 14.29 萬 km2,位于我國西部甘肅省境內(nèi),流域氣候主要受到中高緯度西風(fēng)帶環(huán)流和極地冷氣團(tuán)的影響。該流域從上游到中、下游,以水為紐帶,形成了“冰雪/凍土—森林—草甸—人工/天然綠洲—荒漠—湖泊”的多元自然景觀。流域內(nèi)寒區(qū)和干旱區(qū)并存,地表異質(zhì)性顯著,尤其是山區(qū)冰凍圈和極端干旱的河流尾閭地區(qū)形成了鮮明對比。因此,準(zhǔn)確獲取流域內(nèi)復(fù)雜環(huán)境狀況下的地表蒸散發(fā)數(shù)據(jù),對于深刻理解區(qū)域水循環(huán)和地氣之間能量交換,具有重要意義。此外,該文選擇黑河流域作為研究區(qū)的另一個(gè)重要原因是由中國科學(xué)院西北生態(tài)資源環(huán)境研究院主導(dǎo)開展一系列長期的生態(tài)水文實(shí)驗(yàn),積累了豐富的水文氣象觀測數(shù)據(jù),為研究工作的開展奠定了良好的基礎(chǔ)。

    該研究反演區(qū)域尺度全天候地表蒸散發(fā),涉及的地面觀測數(shù)據(jù)包括自動(dòng)氣象站氣象觀測數(shù)據(jù)和基于渦動(dòng)的蒸散發(fā)觀測數(shù)據(jù)。由中國科學(xué)院西北生態(tài)資源環(huán)境研究院主導(dǎo)的“黑河流域生態(tài)—水文過程綜合遙感觀測聯(lián)合試驗(yàn)”于2012 年5 月底至9 月中旬在黑河流域中游地區(qū)布置了多套渦動(dòng)相關(guān)儀和自動(dòng)氣象站,同步獲取蒸散發(fā)與氣象數(shù)據(jù),為該文全天候蒸散的反演和驗(yàn)證提供了優(yōu)質(zhì)數(shù)據(jù)源。這些觀測數(shù)據(jù)可從寒區(qū)旱區(qū)科學(xué)數(shù)據(jù)中心(http://westdc.westgis.ac.cn/)免費(fèi)下載??紤]到不同站點(diǎn)觀測時(shí)期的差異,該文選取研究區(qū)內(nèi)2012 年6 月1 日至9 月15 日的數(shù)據(jù)開展研究。根據(jù)獲取的驗(yàn)證數(shù)據(jù)選擇了4個(gè)具有代表性站點(diǎn)來探索區(qū)域地表全天候蒸散發(fā)反演,表1 給出了這4 個(gè)站點(diǎn)的編碼名稱、地表類型、高程、測量高度及經(jīng)緯度。

    圖1 研究區(qū)域的地表高程Fig.1 DEM in the study area

    表1 黑河生態(tài)水文遙感試驗(yàn)通量觀測站點(diǎn)信息Table 1 Information of Heihe ecological hydrological remote sensing test flux observation sites

    1.2 遙感數(shù)據(jù)

    與地面觀測數(shù)據(jù)對應(yīng),該文使用2012 年6 月1 日至9 月15 日的MODIS(Moderateresolution Imaging Spectroradiometer)陸表參數(shù)產(chǎn)品作為研究的遙感數(shù)據(jù)源。MODIS 是搭載在Terra 和Aqua 衛(wèi)星上的一個(gè)重要傳感器,具有36 個(gè)波段的觀測通道,覆蓋中國全境,具有較高的空間分辨率和時(shí)間分辨率。目前,基于MODIS 觀測數(shù)據(jù),已經(jīng)發(fā)布了多種陸表和大氣產(chǎn)品數(shù)據(jù),在全球范圍內(nèi)得到了廣泛應(yīng)用。該文使用的MODIS 產(chǎn)品及其參數(shù)如表2 所示。

    表2 MODIS 數(shù)據(jù)詳情表Table 2 MODIS data detail

    MODIS 產(chǎn)品數(shù)據(jù)來源于LAADS DAAC(Level-1 and Atmosphere Archive and Distribution System Distributed Active Archive Center)(https://ladsweb.modaps.eosdis.nasa.gov/search/)網(wǎng)站。該文使用MRT 和MRTSwath 工具對不同級(jí)別的MODIS 數(shù)據(jù)進(jìn)行重投影等預(yù)處理。研究中MOD021KM、MOD03 和MOD05_L2 產(chǎn)品被用來反演地表短波凈輻射,其中MOD021KM 產(chǎn)品空間分辨率為250 m 和500 m 的大氣頂部輻射亮度和反射率數(shù)據(jù)被整合成空間分辨率為1 km 的數(shù)據(jù);MOD35_L2 云掩膜產(chǎn)品被用來判斷MODIS 圖像的每個(gè)像元是否受到云的干擾;MOD09A1、MOD11A1、MOD13A2 和MOD15A2 作為全天候地表蒸散發(fā)反演的參數(shù)輸入,其中MOD09A1 和MOD15A2 在8 d 時(shí)間內(nèi)被認(rèn)為是不變的,MOD13A2 歸一化植被指數(shù)通過hants 濾波后插值得到日尺度的植被指數(shù)[17]。

    1.3 格網(wǎng)氣象數(shù)據(jù)

    研究所用的格網(wǎng)氣象數(shù)據(jù)為中國氣象局陸面數(shù)據(jù)同化系統(tǒng)CLDAS(CMA Land Data Assimilation System)的Verson 1.0 產(chǎn)品(簡稱CLDAS-v1.0),數(shù)據(jù)下載于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)。CLDAS-v1.0 產(chǎn)品覆蓋范圍為東亞區(qū)域(0°~60°N,70°E~150°E),空間分辨率為0.0625°,時(shí)間分辨率為1 h。研究中用到的CLDAS-v1.0 氣象產(chǎn)品要素包括短波輻射、氣溫、風(fēng)速和比濕。這些數(shù)據(jù)被重采樣到MODIS 數(shù)據(jù)1 km空間分辨率。以2012 年6 月1 日為例,圖2 顯示了北京時(shí)間上午11:00 黑河流域的短波輻射、氣溫、風(fēng)速和比濕。

    2 方法

    該文在Leng 等[16]建立的全天候地表蒸散發(fā)反演方法基礎(chǔ)上,深入探索區(qū)域全天候地表蒸散發(fā)的反演。在該反演方法中,所有的像元根據(jù)是否存在有效的地表溫度數(shù)據(jù)被劃分為晴空和有云像元。其中,晴空像元的地表蒸散發(fā)利用逐像元地表溫度—植被指數(shù)特征空間方法進(jìn)行反演,有云像元的地表蒸散發(fā)則基于Penman-Monteith 模型計(jì)算。基于MODIS 陸表產(chǎn)品和CLDAS 格網(wǎng)氣象數(shù)據(jù),探索該全天候蒸散發(fā)反演方法在區(qū)域尺度上的可行性。圖3 為區(qū)域尺度全天候蒸散發(fā)反演的技術(shù)流程圖。

    圖2 2012 年6 月1 日CLDAS 產(chǎn)品:(a)短波輻射;(b)氣溫;(c)風(fēng)速;(d)比濕Fig.2 CLDAS products on June 1,2012

    圖3 全天候蒸散發(fā)反演方法流程圖Fig.3 Flow chart of the all-sky evapotranspiration retrieval method

    2.1 LST-NDVI 特征空間估算晴天像元蒸散發(fā)

    區(qū)別于傳統(tǒng)基于空間信息的“地表溫度—植被指數(shù)”特征空間,Leng 等[16]發(fā)展了獨(dú)立于氣象和下墊面條件的逐像元兩階段“地表溫度—植被指數(shù)”特征空間新模式。圖4 為示意圖。該模式假設(shè)對于每一個(gè)晴空像元(如圖1 中的像元C),都對應(yīng)存在一個(gè)虛擬的特征空間,該虛擬特征空間由 4 個(gè)極端條件(Ts,max,Ts,min,Tv,max,Tv,min)對應(yīng)的地表溫度構(gòu)成。其中,Ts,max和Tv,min兩點(diǎn)構(gòu)成的斜線把梯形特征空間劃分為兩個(gè)三角形,分別為由 Ts,max,Ts,min,Tv,min3 點(diǎn)構(gòu)成的下三角形,由 Ts,max,Tv,max,Tv,min3 點(diǎn)構(gòu)成的上三角形。在實(shí)際應(yīng)用中,只要已知這4 個(gè)極端地表情形(即干燥裸土、受水分脅迫的全植被覆蓋、供水良好的全植被覆蓋和濕潤裸土)對應(yīng)的地表溫度,便可確定像元C 的虛擬特征空間,從而能夠根據(jù)像元C 在該虛擬特征空間中所處位置,估算蒸散發(fā)。該特征空間新模式的一個(gè)顯著優(yōu)點(diǎn)是在實(shí)際應(yīng)用中能夠避免傳統(tǒng)基于空間信息確定干邊和濕邊時(shí)對氣象和下墊面條件的苛刻要求,從而能夠直接反演大范圍蒸散發(fā)。

    圖4 LST-NDVI 特征空間示意圖Fig.4 Schematic diagram of LST-NDVI feature space

    LST-NDVI 特征空間假設(shè)NDVI 為某一值時(shí),蒸散發(fā)隨LST 增大而線性遞增。當(dāng)某像元的LST 和NDVI 落在下三角形中,土壤根層水分充足,植被蒸騰接近于潛在最大值,而表層土壤水分蒸散由像元處于三角形中的位置決定。在該情況下,土壤蒸發(fā)(ETs)和植被蒸騰(ETv)可以表示為:

    式(1)~(2)中,Es,p和Ev,p分別為最大土壤蒸發(fā)量和最大植被蒸騰量。AB和AC為長度(圖4)。

    類似地,當(dāng)某像元的LST 和NDVI 落在上三角形中,則表層土壤已干,土壤蒸發(fā)為0,植被開始受到脅迫,植被蒸騰量取決于該像元處于三角形中的位置。在該情況下,土壤蒸發(fā)(ETs)和植被蒸騰(ETv)可以表示為:

    其中 Ts,max,Ts,min,Tv,max和 Tv,min分別是干燥裸土,飽和裸土,水分脅迫全植被和充分灌溉的全植被的理論地表溫度(K),如圖4 所示。Moran 等[18]基于能量平衡原理,假設(shè)4 種極端條件給出了4 個(gè)理論地表溫度的計(jì)算公式:

    式(7)中,Rn是地表凈輻射(W/m2),rs是表面阻抗。ra是空氣動(dòng)力學(xué)阻抗,γ為干濕表常數(shù)(kPa/℃),Δ 為飽和水汽壓隨空氣溫度變化的斜率(kPa/℃),VPD為參考高度的空氣飽和水汽壓差(kPa),G是土壤熱通量(W/m2),Cp為空氣定壓比熱(J/(kg℃))。其中,rs和Rn參數(shù)的計(jì)算可參考 Leng 等[16]和 Todorovic 等[20]的文獻(xiàn)。

    式(5)中,Rn,s和Rn,v分別為土壤表面凈輻射和植被表面凈輻射(W/m2),G是土壤熱通量(W/m2),Ta是空氣溫度(K)。Cv為空氣定容比熱(J/(℃m3)),γ為干濕表常數(shù)(kPa/℃),Δ 為飽和水汽壓隨空氣溫度變化的斜率(kPa/℃),VPD為參考高度的空氣飽和水汽壓差(kPa)。rcx是葉面氣孔幾乎完全關(guān)閉情況的冠層阻抗,設(shè)定為2 000 s/m;rcp是充足水分條件下的植被冠層阻抗,為100 與葉面積指數(shù)的比值;ra是空氣動(dòng)力學(xué)阻抗。

    地表蒸散發(fā)可以表示為植被蒸騰和土壤蒸發(fā)的和:

    2.2 Penman-Monteith 公式估算有云條件像元蒸散發(fā)

    在有云條件下,無法有效獲取LST 和NDVI 參數(shù),因此使用經(jīng)典的Penman-Monteith模型[19]估算地表蒸散發(fā):

    3 結(jié)果與分析

    3.1 CLDAS 氣象數(shù)據(jù)分析

    該文使用所選4 個(gè)站點(diǎn)自動(dòng)氣象站實(shí)測的短波輻射、風(fēng)速、氣溫、比濕對CLDAS 格網(wǎng)數(shù)據(jù)產(chǎn)品重采樣為1 km 空間分辨率后的格網(wǎng)數(shù)據(jù)進(jìn)行分析,結(jié)果如圖5 所示。

    圖5 CLDAS 數(shù)據(jù)產(chǎn)品與實(shí)測數(shù)據(jù)散點(diǎn)圖:(a)短波輻射;(b)風(fēng)速;(c)氣溫;(d)比濕Fig.5 Comparison of CLDAS products and measured data

    從圖5 可以看出CLDAS 短波輻射、風(fēng)速、氣溫、比濕數(shù)據(jù)誤差較大,尤其是短波輻射和風(fēng)速產(chǎn)品。這可能是由于CLDAS 氣象數(shù)據(jù)在空間上受到氣象站點(diǎn)疏密的影響,導(dǎo)致了在缺少氣象站點(diǎn)的西北地區(qū)數(shù)據(jù)精度下降。Vinukollu 等[21]使用相同的一套驅(qū)動(dòng)數(shù)據(jù)對比驗(yàn)證了P-M 模型[22]、P-T 模型[23]、表面能量平衡模型[5],發(fā)現(xiàn)不同模型之間蒸散發(fā)反演的差異主要取決于凈輻射數(shù)據(jù)的選擇??紤]到短波輻射是地表蒸散發(fā)的能量來源,其精度高低不僅限制了地表蒸散發(fā)的反演結(jié)果,還是地表蒸散發(fā)的主要影響因子?;诖?,該文采用Tang 等[24]的方法利用MODIS 數(shù)據(jù)估算短波凈輻射轉(zhuǎn)換為短波輻射來代替CLDAS 中的短波輻射,作為全天候蒸散發(fā)反演的輸入?yún)?shù)。

    為了驗(yàn)證MODIS 估算的短波凈輻射在黑河地區(qū)的可靠性,該文選取了黑河流域4個(gè)站點(diǎn)的2012 年6 月1 日至9 月15 日的地面有效測量數(shù)據(jù)進(jìn)行對比。其中,地面觀測的短波凈輻射采樣頻率為10 min,MODIS 數(shù)據(jù)反演的短波凈輻射是分辨率為1 km 的瞬時(shí)值。另外,根據(jù)MOD35_L2 云產(chǎn)品判斷每個(gè)像元是否會(huì)受到云的影響。MOD35_L2 中認(rèn)為是clear 和probably clear 的像元,判斷為晴空的像元;MOD35_L2 中認(rèn)為是uncertain和cloudy 的像元,判斷為有云的像元。圖6 是該4 個(gè)站點(diǎn)的地表短波凈輻射觀測值和MODIS 估算的短波凈輻射散點(diǎn)圖,從圖6 可以看出,晴空條件下,基于MODIS 反演的短波凈輻射的均方根誤差小于30 W/m2;有云時(shí),均方根誤差在70 W/m2左右。相比于利用CLDAS 數(shù)據(jù)獲取的短波輻射精度,基于MODIS 反演的短波凈輻射顯然具有更高的精度。

    圖6 估算的地表短波凈輻射與實(shí)測的地表短波凈輻射散點(diǎn)圖(a)1 號(hào)站;(b)7 號(hào)站;(c)10 號(hào)站;(d)17 號(hào)站Fig.6 Comparison of estimated net surface short-wave radiation and measured data

    3.2 全天候區(qū)域地表蒸散發(fā)反演與分析

    為了驗(yàn)證Leng 等[16]方法的可靠性,該文首先使用單點(diǎn)的實(shí)測氣象數(shù)據(jù)和MODIS數(shù)據(jù)反演地表蒸散發(fā),并利用黑河流域4 個(gè)站點(diǎn)的2012 年6 月1 日至9 月15 日的地面有效測量蒸散發(fā)數(shù)據(jù)進(jìn)行對比。驗(yàn)證結(jié)果如圖7 所示。使用單點(diǎn)的氣象數(shù)據(jù)和MODIS數(shù)據(jù)反演的地表蒸散發(fā)RMSE 值在70 W/m2以內(nèi),具有較高的精度。這表明了利用Leng 等[16]的方法反演全天候地表蒸散發(fā)是可行的。在此基礎(chǔ)之上,利用CLDAS 短波輻射、風(fēng)速、比濕、氣溫的產(chǎn)品與MODIS 估算了黑河流域的全天候地表蒸散發(fā),并將反演結(jié)果與地面實(shí)測數(shù)據(jù)進(jìn)行了比較,4 個(gè)站點(diǎn)的平均偏差值為30.8 W/m2,平均均方根誤差值為111.9 W/m2,誤差相對較大。采用Tang 等[24]的方法轉(zhuǎn)換得到短波輻射作為全天候蒸散發(fā)反演的輸入?yún)?shù)后,精度有了明顯的提高,其RMSE 在80 W/m2以內(nèi)。

    圖7 不同驅(qū)動(dòng)數(shù)據(jù)估算的蒸散發(fā)和實(shí)測的蒸散發(fā)散點(diǎn)圖(a)1 號(hào)站、(b)7 號(hào)站、(c)10 號(hào)站、(d)17 號(hào)站Fig.7 Comparison of estimated evapotranspiration and measured data

    基于能量平衡逐像元計(jì)算干濕邊的方法具有良好的普適性和物理基礎(chǔ)。然而大量的輸入?yún)?shù)造成了誤差的積累和傳遞。通過提高關(guān)鍵輸入?yún)?shù)的精度,能夠顯著地提高蒸散發(fā)的反演精度。

    4 結(jié)論

    該文以黑河流域?yàn)檠芯繀^(qū),基于遙感數(shù)據(jù)和較高空間分辨率的格網(wǎng)氣象數(shù)據(jù),實(shí)現(xiàn)了區(qū)域尺度地表蒸散發(fā)反演。在驗(yàn)證分析了CLDAS 數(shù)據(jù)精度的基礎(chǔ)之上,利用MODIS數(shù)據(jù)估算的短波輻射代替CLDAS 短波輻射數(shù)據(jù)作為蒸散發(fā)反演模型的輸入?yún)?shù),最后利用2012 年6 月1 日至9 月15 日的CLDAS 數(shù)據(jù)和MODIS 數(shù)據(jù)反演了黑河流域的全天候地表蒸散發(fā),并對結(jié)果進(jìn)行精度檢驗(yàn)。研究表明,基于MODIS 遙感數(shù)據(jù)估算得到精度較高的短波輻射替代CLDAS 中的短波輻射產(chǎn)品,能顯著提高黑河地區(qū)蒸散發(fā)的反演精度。

    雖然通過該文方法得到了精度較高的區(qū)域尺度全天候地表蒸散發(fā),但研究仍然存在以下不足。首先是驗(yàn)證站點(diǎn)較少,且相對空間距離較近,不利于對地表蒸散發(fā)的全面評(píng)價(jià);其次是由于蒸散發(fā)反演和驗(yàn)證過程中涉及不同尺度的數(shù)據(jù)源,該研究沒有充分考慮尺度效應(yīng)對反演結(jié)果可能造成的影響。

    猜你喜歡
    全天候黑河植被指數(shù)
    復(fù)合式雷達(dá)全天候感知數(shù)字孿生系統(tǒng)
    構(gòu)筑全方位全天候全覆蓋預(yù)警信息網(wǎng)
    黑河的“護(hù)衛(wèi)隊(duì)”
    到張掖看黑河
    文學(xué)港(2019年5期)2019-05-24 14:19:42
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    淺述全天候高速公路便攜收費(fèi)終端
    張掖黑河濕地國家級(jí)自然保護(hù)區(qū)
    24時(shí)間萬全體制全天候安全走行
    主要植被指數(shù)在生態(tài)環(huán)評(píng)中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29

    午夜激情福利司机影院| 国产综合精华液| 久久午夜福利片| 青春草国产在线视频| 联通29元200g的流量卡| 看十八女毛片水多多多| 日韩制服骚丝袜av| 超碰av人人做人人爽久久| 久99久视频精品免费| 色尼玛亚洲综合影院| 亚洲,欧美,日韩| 欧美xxⅹ黑人| 又粗又硬又长又爽又黄的视频| 性色avwww在线观看| 亚洲美女搞黄在线观看| 成年版毛片免费区| 国产伦精品一区二区三区视频9| 天美传媒精品一区二区| 最近的中文字幕免费完整| 内地一区二区视频在线| 欧美最新免费一区二区三区| 一级毛片 在线播放| 天天躁日日操中文字幕| 极品教师在线视频| 九草在线视频观看| 神马国产精品三级电影在线观看| 日本午夜av视频| 亚洲最大成人av| 在线观看美女被高潮喷水网站| 午夜爱爱视频在线播放| 亚洲人成网站高清观看| 波多野结衣巨乳人妻| 午夜免费男女啪啪视频观看| xxx大片免费视频| 亚洲国产精品专区欧美| 日韩三级伦理在线观看| 免费电影在线观看免费观看| 少妇的逼好多水| 麻豆国产97在线/欧美| 国产综合精华液| 午夜老司机福利剧场| 在线免费十八禁| 超碰97精品在线观看| 色尼玛亚洲综合影院| 久久精品熟女亚洲av麻豆精品 | 丝袜喷水一区| 国产午夜精品一二区理论片| 91午夜精品亚洲一区二区三区| 日韩av免费高清视频| 一本一本综合久久| 久久久精品免费免费高清| 纵有疾风起免费观看全集完整版 | 国产精品伦人一区二区| 色综合色国产| 最近最新中文字幕免费大全7| 一区二区三区免费毛片| 精品久久久噜噜| 欧美xxⅹ黑人| 亚洲最大成人av| 麻豆乱淫一区二区| 黄色配什么色好看| 日韩中字成人| 亚洲国产最新在线播放| 又黄又爽又刺激的免费视频.| 高清视频免费观看一区二区 | 美女被艹到高潮喷水动态| 麻豆精品久久久久久蜜桃| 99热这里只有是精品在线观看| 水蜜桃什么品种好| 欧美 日韩 精品 国产| 精品国内亚洲2022精品成人| 国产毛片a区久久久久| 1000部很黄的大片| 亚洲熟女精品中文字幕| 丰满人妻一区二区三区视频av| or卡值多少钱| 97热精品久久久久久| 干丝袜人妻中文字幕| 亚洲三级黄色毛片| 看十八女毛片水多多多| 亚洲精品视频女| 亚洲av免费在线观看| 嫩草影院新地址| 欧美另类一区| 国产又色又爽无遮挡免| 国产伦一二天堂av在线观看| 色5月婷婷丁香| 欧美丝袜亚洲另类| av黄色大香蕉| 国产激情偷乱视频一区二区| 免费少妇av软件| 久久这里有精品视频免费| 国产午夜精品论理片| 麻豆久久精品国产亚洲av| 大话2 男鬼变身卡| 亚洲av免费高清在线观看| 亚洲av电影在线观看一区二区三区 | 欧美一级a爱片免费观看看| a级一级毛片免费在线观看| 天堂av国产一区二区熟女人妻| 日日摸夜夜添夜夜爱| 插阴视频在线观看视频| 久久久久精品性色| 成人无遮挡网站| 麻豆av噜噜一区二区三区| 免费看日本二区| 亚洲天堂国产精品一区在线| 啦啦啦啦在线视频资源| 久久久久久久久久久丰满| 人人妻人人澡人人爽人人夜夜 | 国产极品天堂在线| 久久午夜福利片| 日本一二三区视频观看| 亚洲av不卡在线观看| 国产亚洲精品av在线| 美女脱内裤让男人舔精品视频| 男女边摸边吃奶| 舔av片在线| 久久国内精品自在自线图片| 久久久a久久爽久久v久久| 99视频精品全部免费 在线| 亚洲av成人精品一二三区| 欧美另类一区| 高清欧美精品videossex| 天堂影院成人在线观看| 人体艺术视频欧美日本| 男女下面进入的视频免费午夜| 综合色av麻豆| 精品人妻熟女av久视频| 少妇的逼水好多| 国产一级毛片七仙女欲春2| 日韩欧美精品v在线| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久久久久久av| 日本与韩国留学比较| 亚洲欧美成人精品一区二区| 成人高潮视频无遮挡免费网站| 欧美最新免费一区二区三区| 边亲边吃奶的免费视频| 久久久久久久久大av| av在线播放精品| 久久人人爽人人片av| 国产欧美日韩精品一区二区| .国产精品久久| 国产男人的电影天堂91| 久久久久性生活片| 欧美不卡视频在线免费观看| 97在线视频观看| 日本黄色片子视频| 日本免费在线观看一区| 国产精品精品国产色婷婷| 狠狠精品人妻久久久久久综合| 精品不卡国产一区二区三区| 中国国产av一级| 久久久久久久国产电影| ponron亚洲| 18禁在线播放成人免费| 一级爰片在线观看| 精品人妻视频免费看| av卡一久久| 精品久久久噜噜| 久久99精品国语久久久| 2021少妇久久久久久久久久久| 日韩欧美一区视频在线观看 | 51国产日韩欧美| 精品少妇黑人巨大在线播放| 国产亚洲91精品色在线| 成人亚洲精品av一区二区| 亚洲第一区二区三区不卡| 亚洲最大成人av| 内地一区二区视频在线| 边亲边吃奶的免费视频| 六月丁香七月| 看黄色毛片网站| 五月伊人婷婷丁香| 国产精品美女特级片免费视频播放器| av天堂中文字幕网| 色哟哟·www| 亚洲四区av| 亚洲精品乱久久久久久| 免费播放大片免费观看视频在线观看| 丰满乱子伦码专区| 男人和女人高潮做爰伦理| 午夜精品一区二区三区免费看| 波野结衣二区三区在线| 亚洲欧美日韩无卡精品| or卡值多少钱| 欧美xxxx性猛交bbbb| 亚洲人成网站在线观看播放| 纵有疾风起免费观看全集完整版 | 亚洲美女搞黄在线观看| 日韩一本色道免费dvd| 69人妻影院| 国产一区有黄有色的免费视频 | 日本欧美国产在线视频| 黄色配什么色好看| 亚洲欧洲日产国产| 亚洲av不卡在线观看| 国产成人精品久久久久久| 18禁在线播放成人免费| 老师上课跳d突然被开到最大视频| 午夜日本视频在线| 亚洲av在线观看美女高潮| 大香蕉久久网| 国产精品国产三级专区第一集| 亚洲综合色惰| 亚洲最大成人av| 午夜免费激情av| 免费黄频网站在线观看国产| 全区人妻精品视频| 亚洲av成人精品一二三区| 日日摸夜夜添夜夜爱| 国产91av在线免费观看| 身体一侧抽搐| 天美传媒精品一区二区| 少妇的逼水好多| 日韩制服骚丝袜av| 国产单亲对白刺激| 日韩精品青青久久久久久| 天天躁日日操中文字幕| 麻豆成人av视频| 亚州av有码| 超碰97精品在线观看| 永久网站在线| av国产久精品久网站免费入址| 免费av观看视频| 搡女人真爽免费视频火全软件| 久久久久精品性色| 三级经典国产精品| 免费av毛片视频| 国产一级毛片七仙女欲春2| 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 亚洲欧美成人综合另类久久久| 乱人视频在线观看| 久久精品夜色国产| 免费av观看视频| 蜜臀久久99精品久久宅男| 在线观看美女被高潮喷水网站| 婷婷色av中文字幕| 久久久亚洲精品成人影院| 国产免费视频播放在线视频 | 国产精品av视频在线免费观看| 最近的中文字幕免费完整| 最近中文字幕2019免费版| 一级二级三级毛片免费看| 99热这里只有是精品50| 97人妻精品一区二区三区麻豆| 亚洲av免费高清在线观看| 国产成人精品福利久久| 在线观看美女被高潮喷水网站| 国产色婷婷99| 97超视频在线观看视频| 青春草国产在线视频| 国产成人精品福利久久| 成人一区二区视频在线观看| a级毛色黄片| 青春草视频在线免费观看| 国内精品一区二区在线观看| 男女啪啪激烈高潮av片| 亚洲欧洲国产日韩| 成人午夜高清在线视频| 国产精品伦人一区二区| 成人高潮视频无遮挡免费网站| 男人舔女人下体高潮全视频| 能在线免费看毛片的网站| 中文字幕av在线有码专区| 亚洲精品乱码久久久久久按摩| 超碰av人人做人人爽久久| 成人美女网站在线观看视频| 久久精品久久久久久噜噜老黄| 国产麻豆成人av免费视频| av福利片在线观看| 少妇猛男粗大的猛烈进出视频 | 全区人妻精品视频| 国模一区二区三区四区视频| 少妇高潮的动态图| 两个人视频免费观看高清| 亚洲av中文字字幕乱码综合| 欧美精品国产亚洲| 免费黄频网站在线观看国产| 亚洲伊人久久精品综合| 欧美高清性xxxxhd video| 亚洲真实伦在线观看| 久久久色成人| 在线天堂最新版资源| 69人妻影院| 欧美变态另类bdsm刘玥| 一级a做视频免费观看| 在线观看人妻少妇| 国产成人免费观看mmmm| 伊人久久国产一区二区| 亚洲国产av新网站| 嘟嘟电影网在线观看| 久久精品人妻少妇| 简卡轻食公司| 大话2 男鬼变身卡| 亚洲av免费高清在线观看| 国产男女超爽视频在线观看| 丝袜美腿在线中文| 国产毛片a区久久久久| 狂野欧美激情性xxxx在线观看| 国产成人精品一,二区| 最新中文字幕久久久久| 大香蕉97超碰在线| 全区人妻精品视频| 午夜免费激情av| 成人亚洲精品一区在线观看 | 国产精品国产三级专区第一集| 亚洲精品国产av成人精品| 国产精品不卡视频一区二区| 欧美最新免费一区二区三区| 少妇的逼水好多| 国产免费福利视频在线观看| 久久久成人免费电影| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 亚洲欧美精品专区久久| 嫩草影院新地址| 亚洲av电影不卡..在线观看| 国产有黄有色有爽视频| 国产高清国产精品国产三级 | 最新中文字幕久久久久| 成人av在线播放网站| 国产欧美日韩精品一区二区| 五月伊人婷婷丁香| 99久久精品国产国产毛片| 天天一区二区日本电影三级| 久久久久久久久久黄片| 欧美3d第一页| 国产成人精品福利久久| 国产亚洲午夜精品一区二区久久 | 精品人妻一区二区三区麻豆| 国产爱豆传媒在线观看| 色播亚洲综合网| 国产一区有黄有色的免费视频 | 我要看日韩黄色一级片| 狠狠精品人妻久久久久久综合| 永久免费av网站大全| 欧美日韩视频高清一区二区三区二| 中文字幕人妻熟人妻熟丝袜美| 2022亚洲国产成人精品| 亚洲国产最新在线播放| 免费av观看视频| 夫妻午夜视频| 亚洲无线观看免费| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 少妇猛男粗大的猛烈进出视频 | 日韩精品有码人妻一区| 亚洲美女搞黄在线观看| 日韩不卡一区二区三区视频在线| 亚洲国产精品国产精品| 国产精品一及| 亚洲精品一区蜜桃| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最近最新中文字幕大全电影3| www.色视频.com| 久久久久久久久久久丰满| 国产精品av视频在线免费观看| 内地一区二区视频在线| 久久久久久伊人网av| 国产久久久一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品456在线播放app| 国产精品福利在线免费观看| 一级av片app| 亚洲熟女精品中文字幕| 亚洲精品国产av蜜桃| av免费观看日本| 欧美成人精品欧美一级黄| 日本色播在线视频| 18禁动态无遮挡网站| 3wmmmm亚洲av在线观看| 国产在线一区二区三区精| 亚洲av成人精品一二三区| 亚洲不卡免费看| 久久这里只有精品中国| 国产亚洲av片在线观看秒播厂 | 成年av动漫网址| 日韩精品有码人妻一区| 天堂影院成人在线观看| 七月丁香在线播放| 欧美日韩亚洲高清精品| 国产午夜福利久久久久久| av在线天堂中文字幕| 蜜桃久久精品国产亚洲av| av免费观看日本| 久久久色成人| 波多野结衣巨乳人妻| 国产精品久久久久久精品电影| 性色avwww在线观看| 亚洲精品一区蜜桃| 国产 一区精品| 97精品久久久久久久久久精品| 亚洲av日韩在线播放| 国模一区二区三区四区视频| 亚洲,欧美,日韩| 日日撸夜夜添| 能在线免费观看的黄片| 午夜福利视频精品| 日日干狠狠操夜夜爽| 18+在线观看网站| 国产精品三级大全| 91狼人影院| 国产黄片视频在线免费观看| 伦精品一区二区三区| 免费观看性生交大片5| 精品酒店卫生间| 亚洲精品影视一区二区三区av| 黄色日韩在线| 亚洲最大成人手机在线| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 久久99蜜桃精品久久| 国产v大片淫在线免费观看| eeuss影院久久| 中文字幕免费在线视频6| 热99在线观看视频| 十八禁国产超污无遮挡网站| 97在线视频观看| 尾随美女入室| 一个人观看的视频www高清免费观看| 日韩国内少妇激情av| 激情五月婷婷亚洲| 免费观看性生交大片5| 久热久热在线精品观看| 亚洲精品久久久久久婷婷小说| 毛片女人毛片| 午夜精品一区二区三区免费看| 精品一区二区免费观看| 亚洲美女搞黄在线观看| 亚洲欧洲国产日韩| 欧美高清性xxxxhd video| 亚洲精品日韩av片在线观看| 精品一区二区三区视频在线| 亚洲国产高清在线一区二区三| 天天躁日日操中文字幕| 国产午夜精品论理片| 国产av不卡久久| 极品少妇高潮喷水抽搐| 亚洲自拍偷在线| 久久国内精品自在自线图片| 久久久久久久久久成人| 久久人人爽人人爽人人片va| 精品久久久久久成人av| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| av.在线天堂| 欧美成人a在线观看| 免费看光身美女| 美女xxoo啪啪120秒动态图| 日韩成人伦理影院| 一级a做视频免费观看| 国产淫语在线视频| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久久成人| 亚洲综合色惰| 日韩欧美一区视频在线观看 | 日韩国内少妇激情av| 大片免费播放器 马上看| 国产有黄有色有爽视频| 三级毛片av免费| 少妇人妻一区二区三区视频| 久久久久久久久久久丰满| 777米奇影视久久| 午夜免费男女啪啪视频观看| 亚洲aⅴ乱码一区二区在线播放| 精品国产露脸久久av麻豆 | 精品国产一区二区三区久久久樱花 | 亚洲av国产av综合av卡| 成人亚洲精品av一区二区| 色吧在线观看| 欧美高清成人免费视频www| 久久97久久精品| 最近中文字幕2019免费版| 国产精品久久久久久精品电影小说 | 欧美高清成人免费视频www| 国产精品久久久久久久电影| 一本一本综合久久| 汤姆久久久久久久影院中文字幕 | 亚洲真实伦在线观看| 国产91av在线免费观看| 黑人高潮一二区| kizo精华| 亚洲精品影视一区二区三区av| 人妻夜夜爽99麻豆av| 久久精品久久久久久久性| 日韩精品有码人妻一区| 日韩电影二区| 日本爱情动作片www.在线观看| 插逼视频在线观看| 狂野欧美白嫩少妇大欣赏| 黄片wwwwww| 国产色爽女视频免费观看| 日韩强制内射视频| 国产精品1区2区在线观看.| 在线免费十八禁| 亚洲国产精品sss在线观看| 亚洲精品一区蜜桃| 丰满乱子伦码专区| 网址你懂的国产日韩在线| 女人十人毛片免费观看3o分钟| 日本-黄色视频高清免费观看| 亚洲人与动物交配视频| 免费播放大片免费观看视频在线观看| 日韩电影二区| 免费观看无遮挡的男女| 欧美精品一区二区大全| 极品教师在线视频| 亚洲av成人av| 亚洲,欧美,日韩| 亚洲综合精品二区| 婷婷色麻豆天堂久久| 亚洲人成网站在线播| 亚洲18禁久久av| 亚洲性久久影院| 亚洲综合精品二区| 亚洲乱码一区二区免费版| 亚洲精品一区蜜桃| 国产精品爽爽va在线观看网站| 一个人看的www免费观看视频| 免费观看的影片在线观看| 久久99精品国语久久久| 国产亚洲av片在线观看秒播厂 | 午夜精品一区二区三区免费看| 深爱激情五月婷婷| 美女脱内裤让男人舔精品视频| 97精品久久久久久久久久精品| 久久精品夜夜夜夜夜久久蜜豆| 亚洲丝袜综合中文字幕| 色播亚洲综合网| 国产永久视频网站| 欧美成人午夜免费资源| 久久久精品94久久精品| 看免费成人av毛片| 好男人在线观看高清免费视频| 日本与韩国留学比较| 国产精品久久久久久久久免| 搡女人真爽免费视频火全软件| 国产精品麻豆人妻色哟哟久久 | 最近中文字幕2019免费版| 久久精品国产亚洲网站| 国产成人免费观看mmmm| 高清视频免费观看一区二区 | 亚洲乱码一区二区免费版| 免费黄频网站在线观看国产| 日韩成人av中文字幕在线观看| 国产午夜福利久久久久久| 国产一级毛片在线| 国产亚洲91精品色在线| 少妇高潮的动态图| 三级毛片av免费| 亚洲av男天堂| 日韩在线高清观看一区二区三区| 亚洲精品国产av成人精品| 亚洲成人精品中文字幕电影| 亚洲在线观看片| 欧美zozozo另类| 天天躁夜夜躁狠狠久久av| 亚洲av在线观看美女高潮| 嫩草影院精品99| 欧美日韩一区二区视频在线观看视频在线 | 最后的刺客免费高清国语| 免费看a级黄色片| 国产高清三级在线| 国产高清有码在线观看视频| 22中文网久久字幕| 伊人久久国产一区二区| av在线天堂中文字幕| 亚洲欧美精品自产自拍| 狠狠精品人妻久久久久久综合| 久久久精品94久久精品| 在线天堂最新版资源| 91精品国产九色| 一级a做视频免费观看| 波野结衣二区三区在线| 国产亚洲一区二区精品| 中文精品一卡2卡3卡4更新| 国产精品.久久久| 免费在线观看成人毛片| 亚洲精品色激情综合| 黄色一级大片看看| 人妻少妇偷人精品九色| 久久久久久伊人网av| 亚洲在久久综合| 久99久视频精品免费| 国产精品久久久久久精品电影| 欧美+日韩+精品| 国产人妻一区二区三区在| 精品国产三级普通话版| 久久久久免费精品人妻一区二区| 99久国产av精品国产电影| 免费播放大片免费观看视频在线观看| 97超碰精品成人国产| 日本一二三区视频观看| 婷婷色综合大香蕉| 在线观看人妻少妇| 91久久精品国产一区二区三区| 国产又色又爽无遮挡免| 少妇高潮的动态图| 日韩欧美精品免费久久| 亚洲欧美精品自产自拍| 亚洲欧美一区二区三区国产| 国产av在哪里看| 国产单亲对白刺激| 国产 亚洲一区二区三区 | 狠狠精品人妻久久久久久综合| 国产av不卡久久| 免费大片黄手机在线观看| 一本一本综合久久| 五月伊人婷婷丁香| 亚洲精品乱码久久久v下载方式| 日日撸夜夜添| 中国美白少妇内射xxxbb|