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

    落葉闊葉林冠層光合有效輻射分量的遙感模擬與分析

    2017-06-22 13:47:48梁守真隋學(xué)艷侯學(xué)會姚慧敏馬萬棟
    生態(tài)學(xué)報 2017年10期
    關(guān)鍵詞:闊葉林植被指數(shù)冠層

    梁守真, 隋學(xué)艷,侯學(xué)會, 王 猛, 姚慧敏,馬萬棟

    1 山東省農(nóng)業(yè)可持續(xù)發(fā)展研究所,濟(jì)南 250100 2 中國科學(xué)院深圳先進(jìn)技術(shù)研究院,深圳 518055 3 環(huán)境保護(hù)部衛(wèi)星環(huán)境應(yīng)用中心,北京 100094

    ?

    落葉闊葉林冠層光合有效輻射分量的遙感模擬與分析

    梁守真1,2,*, 隋學(xué)艷1,侯學(xué)會1, 王 猛1, 姚慧敏1,馬萬棟3

    1 山東省農(nóng)業(yè)可持續(xù)發(fā)展研究所,濟(jì)南 250100 2 中國科學(xué)院深圳先進(jìn)技術(shù)研究院,深圳 518055 3 環(huán)境保護(hù)部衛(wèi)星環(huán)境應(yīng)用中心,北京 100094

    冠層綠色葉片(光合組分)的光合有效輻射分量(綠色FPAR)真實地反映了植被與外界進(jìn)行物質(zhì)和能量交換的能力,獲取冠層光合組分吸收的太陽光合有效輻射,對生態(tài)系統(tǒng)生產(chǎn)力的遙感估算精度的提高具有重要的意義。研究以落葉闊葉林為例,基于SAIL模型模擬森林冠層光合組分和非光合組分吸收的光合有效輻射,研究冠層FPAR變化規(guī)律以及與植被指數(shù)的相關(guān)關(guān)系。結(jié)果表明,冠層結(jié)構(gòu)的改變會影響冠層對PAR的吸收能力,冠層綠色FPAR的大小與植被面積指數(shù)及光合組分面積比相關(guān);在高覆蓋度植被區(qū),冠層綠色FPAR占冠層總FPAR的80%以上,非光合組分的貢獻(xiàn)較小,但在低植被覆蓋區(qū),當(dāng)光合組分和非光合組分面積相同時,綠色FPAR不及冠層總FPAR的50%;相比于NDVI,北方落葉闊葉林冠層EVI與綠色FPAR存在更為顯著的線性相關(guān)關(guān)系(R2>0.99)。

    FPAR; 光合組分; NDVI; EVI; 落葉闊葉林

    光合有效輻射吸收分量FPAR (Fraction of Absorbed Photosynthetically Active Radiation) 是指植被吸收的光合有效輻射PAR (Photosynthetically Active Radiation)在入射PAR中所占的比例。它表征了植被冠層對太陽光能的截獲、吸收能力,是陸地碳循環(huán)研究中的一個重要參數(shù)[1- 2]。植被冠層FPAR主要通過兩種途徑獲?。旱孛娑ㄎ挥^測和基于遙感數(shù)據(jù)的反演[3]。相對于地面定位觀測,遙感可同步獲取不同位置的FPAR,時空性好,是當(dāng)前獲取區(qū)域、全球尺度FPAR 的唯一可行手段。目前FPAR遙感反演研究的重點在于如何提高反演精度,使其能真正反映植被吸收太陽輻射、進(jìn)行光合作用的能力[4]。

    植被冠層由葉片和木本組織組成,按照能否進(jìn)行光合作用,可分為光合組分PAV (Photosynthetically Active Vegetation,綠色葉片)和非光合組分NPV (Non-Photosynthetically Active Vegetation,如樹枝等木本部分)[5- 7]。NPV與PAV無論是組織結(jié)構(gòu)還是輻射特性,都顯著不同[8]。光合組分吸收PAR進(jìn)行光合作用,對生物量的累積起著關(guān)鍵作用;非光合組分也吸收PAR,但它并不進(jìn)行光合作用,或者即使能進(jìn)行光合作用也只是對其呼吸釋放碳的再固定,其光合作用并不增加生態(tài)系統(tǒng)的總初級生產(chǎn)力GPP[9- 12]。森林冠層中正確估算冠層有效的PAR,需要區(qū)分來自不同組分的PAR,提取冠層PAV的FPAR (FPARPAV),即綠色FPAR。在早期的研究中,草地綠色FPAR被描述為綠色生物量占總生物量的比例與冠層總FPAR的乘積,通過計算綠色生物量的比率計算綠色FPAR[13]。但是森林冠層比草地復(fù)雜,冠層NPV和PAV分離困難,因此,其綠色FPAR的研究多采用輻射傳輸模型模擬和實驗室測量的方法[7-8,14]。Xiao等將森林冠層綠色FPAR描述為增強(qiáng)型植被指數(shù)EVI(Enhanced Vegetation Index)的線性函數(shù),但是并未從機(jī)理上給出兩者之間的定量關(guān)系,因此Xiao 等認(rèn)為需要進(jìn)一步研究綠色FPAR 的光譜探測機(jī)理,以更好地定量描述植被指數(shù)與綠色FPAR 之間的關(guān)系[15]。

    輻射傳輸模型可模擬太陽光子在植被內(nèi)部的傳輸過程與作用機(jī)制,它的出現(xiàn)和發(fā)展為冠層的分離提供了可能。本研究以北方落葉闊葉林為研究對象,從輻射傳輸模型SAIL(Scattering by Arbitrary Inclined Leaves)出發(fā),建立冠層光譜和FPAR模擬數(shù)據(jù)集,分析冠層綠色FPAR (FPARPAV)的特征及其與冠層光譜之間的相關(guān)關(guān)系。主要目標(biāo)包含兩個方面,一是確定森林冠層綠色FPAR隨冠層結(jié)構(gòu)的變化特征;二是分析綠色FPAR與冠層光譜以及植被指數(shù)的關(guān)系,探討綠色FPAR (FPARPAV)的機(jī)理和可能的估算方法。

    1 方法

    1.1 模型

    輻射傳輸模型是目前相對成熟的基于物理光學(xué)基礎(chǔ)的模型,它可比較系統(tǒng)、完整地描述植被冠層與入射輻射之間相互作用過程和特征,被廣泛應(yīng)用于遙感FPAR的理論和反演研究[16-17]。SAIL模型是輻射傳輸模型的典型代表,是Verhoef在對Suits 模型進(jìn)行改進(jìn)的基礎(chǔ)上提出的四流輻射傳輸模型[18]。該模型將植被冠層假定為方位隨機(jī)分布、水平均勻無限擴(kuò)展的混濁介質(zhì),只在垂直方向上有限變化。冠層可分為一層或多層,每層中的植被光學(xué)特性參數(shù)和結(jié)構(gòu)特性參數(shù)為常數(shù),并且每層的植被元素視作吸收和散射微粒,考慮了多次散射的作用。該模型可以完美地描述連續(xù)植被冠層的光學(xué)特性,被遙感學(xué)術(shù)界廣泛應(yīng)用[19]。但由于SAIL模型將冠層假設(shè)為混沌介質(zhì),沒有考慮離散冠層的集聚效應(yīng)以及冠層之間的遮蔽效應(yīng)[20],被認(rèn)為在模擬離散的森林冠層時存在一定的局限。離散森林的冠層表現(xiàn)為不同形狀的非均質(zhì)結(jié)構(gòu),具有強(qiáng)烈的三維效應(yīng),需要借助三維模型才能描述復(fù)雜冠層的輻射機(jī)制,但依據(jù)當(dāng)前的條件,建立完善的描述陸地生態(tài)系統(tǒng)結(jié)構(gòu)的三維模型仍不太現(xiàn)實,且三維模型計算更加復(fù)雜,運算耗時長[19]。Gobron在模擬時發(fā)現(xiàn)三維效應(yīng)引起的擾動相當(dāng)穩(wěn)定,并且Widlowski的研究也發(fā)現(xiàn)一維模型幾乎可模擬三維模型生成的反射率[21-22]。因此,復(fù)雜冠層水平一致性的簡化假設(shè)并不會導(dǎo)致一維SAIL模型失效[19]。此外,對于常用的遙感影像,如Landsat/TM、MODIS、MERIS、SPOT/VEGETATION數(shù)據(jù),其空間分辨率通常較低,并無法分辨出單個冠層,在這一分辨率下,將像元看作水平均一的介質(zhì)可能會獲得更好的效果[23- 25]。SAIL模型的分層假設(shè),使其可以計算每個層吸收的PAR[20]。但目前在大多數(shù)研究在應(yīng)用SAIL模型時,整個冠層被視為一個層,組分在層內(nèi)隨機(jī)分布。如果將冠層中的每個組分分離到單獨的層中,分別計算,而不是將其混合為單個層,則可確定每個組分吸收的PAR,從而計算出每個組分的FPAR。

    在本研究中,針葉林冠層被分解為兩層,第一層為PAV(葉片)層,第二層為NPV(樹枝)層。 SAIL模型程序?qū)㈨攲拥南滦腥肷渫?即太陽輻射通量)設(shè)為1.0,冠層的綠色FPAR (FPARPAV)可通過葉片層內(nèi)通量的輸入與輸出平衡來計算:

    FPAR = FLUX3(1)+FLUX2(1)-FLUX1(1)-FLUX3(2)- FLUX2(2)+FLUX1(2)

    (1)

    式中,FLUX1為上行散射通量,FLUX2為下行散射通量,FLUX3為下行直射通量(1.0);“()”中數(shù)字表示層號,1為頂層。

    SAIL模型的輸入?yún)?shù)包括冠層組分光譜(反射率和透射率)、冠層結(jié)構(gòu)(組分傾角分布、組分面積指數(shù))、土壤反射率、太陽天頂角、觀測幾何參數(shù)(觀測天頂角和相對方位角) 以及入射輻射漫輻射比,其輸出包括冠層反射率和每個層的輸入輸出通量。植被面積指數(shù)PAI(Plant Area Index)包括樹枝和葉片面積指數(shù),對于非光合組分,其面積指數(shù)占冠層PAI的比例在0.02—0.5之間變化(冬季除外)[20,26-27],本研究設(shè)置了7種不同的場景,葉片層面積指數(shù)占冠層PAI的比例分別為0.5,0.6,0.7,0.8,0.9,0.95,0.98。各參數(shù)值如表1和圖1所示。

    PAV: Photosynthetically active vegetation; NPV: Non-photosynthetic vegetation

    圖1 SAIL模型采用的冠層組分和土壤光譜 Fig.1 The spectrum of canopy components and soil used in SAIL model

    研究所采用的落葉闊葉林(白橡、白楊和紅榿)冠層組分樹枝和葉片光譜來自NASA ORNL DAAC。其中白楊冠層組分樣本采自于明尼蘇達(dá)超級國家森林公園SNF(Superior National Forest in Minnesota),白橡和紅榿冠層組分來自俄勒岡森林樣帶OTTER(Oregon Transect Ecosystem Research)。研究人員分別在1984年和1990年夏季(生長旺盛期)對森林冠層組分樣本采樣,并在實驗室內(nèi)利用LiCor積分球和Cary- 14光譜儀測量其反射率和透射率[28-29]。本研究對收集的3種闊葉林的高光譜光譜進(jìn)行平均,并以平均的方法計算3個寬波段藍(lán)波段(459—479nm)、紅波段(620—670nm)、近紅外波段(841—876nm)以及PAR波段(400—700nm)的反射率和透射率。土壤光譜數(shù)據(jù)取自約翰·霍普金斯大學(xué) (Johns Hopkins University)采集的淺灰色粉砂壤光譜,其數(shù)據(jù)處理方法與冠層組分光譜相同。

    1.2 植被指數(shù)

    基于植被指數(shù)反演冠層FPAR是遙感領(lǐng)域常用的一種方式,也有研究人員采用實地測量的LAI構(gòu)建反演模型[30- 32]。由于大部分儀器測量并不能區(qū)分葉片和木本組織,測量值多為反映植被冠層總FPAR的量,導(dǎo)致當(dāng)前的大多數(shù)模型是建立在冠層總FPAR和波譜植被指數(shù)之間。在當(dāng)前的植被冠層FPAR遙感研究中,以紅波段、近紅波段等為主構(gòu)建的植被指數(shù)被廣泛應(yīng)用,包括NDVI[33- 35]、EVI[10,36- 38]、DVI[39]等,其中使用最多的是NDVI。在本研究中,NDVI和EVI兩種植被指數(shù)被選用。NDVI和EVI是MODIS植被指數(shù)產(chǎn)品中的兩個指數(shù),NDVI對葉綠素敏感,而EVI對植被冠層結(jié)構(gòu)更加敏感,兩者互為補(bǔ)充[40-41]。其公式分別如下:

    (2)

    (3)

    式中,RNIR,RRed,RBlue分別對應(yīng)MODIS傳感器的近紅外波段(841—876nm)、紅波段(620—670nm)和藍(lán)波段(459—479nm)的反射率。

    2 結(jié)果

    2.1 冠層總FPAR與綠色FPAR

    總FPAR反映了植被冠層吸收的太陽PAR的比例。模擬顯示,冠層總FPAR在分層模擬和與單層模擬結(jié)果一致,FPARPAV和FPARNPV(綠色FPAR)之和與單層的冠層FPAR相等。落葉闊葉林冠層總FPAR隨PAI、PAV的面積比(PAV面積占冠層PAI的比例)的變化而變化(圖2)。在同一PAV面積比條件下,PAI的增加會提高落葉闊葉林冠層的總FPAR。在設(shè)定情景下,冠層總FPAR的最大值為0.958(PAI=7,PAV面積比為0.98)。總FPAR的增加速率隨著PAI的增加而逐漸減小,即在低PAI情況下,冠層總FPAR隨PAI增加迅速,而在高PAI情況下,冠層總FPAR增加緩慢。例如當(dāng)PAI從0.5增加到1時,總FPAR增加值超過0.20;而PAI從6增加到7,冠層總FPAR最大增加量僅為0.012(PAV面積比為0.98時)。當(dāng)PAI<7時,在相同的PAI條件下,PAV面積所占比重越高,冠層總FPAR值越小,而在PAI為7時,冠層總FPAR在PAV為0.98時最大。

    圖2 闊葉林冠層總FPAR和綠色FPAR Fig.2 The change of total FPAR and FPARPAV of canopy

    對于光合組分的FPAR,其變化特征與冠層總FPAR存在一定差異。在同一PAI下,PAV比重越大,FPARPAV越大。 FPARPAV隨PAV的變化幅度與PAI有關(guān)。在PAI為2 時,增幅最為明顯,表現(xiàn)為PAV面積比每增加0.1,FPARPAV增加超過0.03;而PAI為0.1時,PAV面積比增加相同的量,FPARPAV僅增加0.006。在同一PAV面積比下,FPARPAV隨PAI的增加而增加,在PAI為7時,FPARPAV達(dá)到最大值。

    在設(shè)定情景下,FPARPAV占總FPAR的比例在43.07%—99.74%之間波動(圖3)。PAI和PAV面積比的增加都可提高FPARPAV在冠層總FPAR的比重。但PAI的變化會影響FPARPAV比率的變動幅度。在低PAI(0—1)情況下,FPARPAV占總FPAR的比重隨著PAV面積比的增加迅速增加,從45.86%增加到97.56%;而在高PAI(5—7)情況下,當(dāng)PAV面積比從0.5增加到0.98時,FPARPAV的比重僅增加了16.8%,其增幅低于低PAI時增加量。并且在高覆蓋度情況下,FPARPAV占總FPAR的比重超過80%,當(dāng)冠層PAV面積比為0.98時,FPARPAV占總FPAR的99.60%,兩者非常接近。但是在低植被覆蓋度情況下,FPARPAV占總FPAR的比重明顯低于總FPAR,尤其是當(dāng)PAI面積比為0.5時,僅占總FPAR的45.86%。

    圖3 冠層FPARPAV比率與植被面積指數(shù),PAV面積比的關(guān)系Fig.3 The relationship between percent FPARPAV and Plant area index,percent PAV area

    2.2 落葉闊葉林冠層波譜

    根據(jù)圖4,在藍(lán)波段、紅波段和近紅外波段,隨著波長的增加,落葉闊葉林冠層的反射率逐漸升高。冠層在藍(lán)波段有較低的反射率,最高僅為0.074(PAI =0.1,PAV面積比為0.98);而在近紅外波段有高的反射率,超過0.36。在不同的波段冠層反射率對PAI和組分比例的響應(yīng)存在一定的差異。在藍(lán)、紅波段,冠層反射率主要受PAI的影響,組分比例的變動對冠層反射率的影響甚微。隨著PAI的增加,冠層在兩個波段的反射率值逐漸降低,尤其是紅波段反射率下降明顯,從25.04%降低到2.33%。在近紅外波段,冠層反射率對NPV面積比和PAI都非常敏感。在冠層PAI恒定時,隨著PAV比例的升高,冠層反射率升高,這種現(xiàn)象在中等覆蓋度時最為顯著。如PAI=3時,隨這PAV比例從0.5增加到0.98,冠層在近紅外的反射率分別為0.3256、0.3375、0.3495、0.3618、0.3746、0.3814、0.3856。在同一PAV面積比情況下,在PAI為2時最低,PAI為7時最高(0.41)。藍(lán)波段和紅波段與近紅外波段反射率之間的差異隨著PAI的增加逐漸加大,這種變化特征獨立于與冠層組分比例。

    圖4 冠層結(jié)構(gòu)對冠層反射率的影響Fig.4 The influence of canopy structure on canopy reflectance

    2.3 植被指數(shù)與FPAR關(guān)系

    根據(jù)圖5,相同條件下,NDVI總是大于EVI,冠層NDVI和EVI都隨著PAI的增加而升高。當(dāng)PAI較低和較高時,NDVI對PAV的變化不敏感,但在中等覆蓋度條件下(3

    圖5 植被指數(shù)與冠層結(jié)構(gòu)的關(guān)系Fig.5 The relationship between vegetation indices and canopy architecture

    無論NDVI還是EVI,它們與冠層總FPAR都有非常顯著的相關(guān)關(guān)系(R2>0.9)(圖6),但NDVI與總FPAR之間關(guān)系對PAV面積比的變化并不敏感,意味著冠層PAV面積比的波動不會影響NDVI與冠層FPAR的關(guān)系;但在中高覆蓋度條件下,冠層PAV面積比的變動會對 EVI與冠層總FPAR的關(guān)系產(chǎn)生影響。在PAV面積比為0.5時,EVI與冠層總FPAR之間線性方程的R2為0.9978,而在PAV面積比為0.98時,R2為0.9982??傮w上表現(xiàn)為PAV比例越高,R2越大,EVI與冠層總FPAR的關(guān)系越緊密。NDVI、EVI與FPARPAV也表現(xiàn)出了顯著的線性關(guān)系(圖7),但相比于NDVI,EVI與FPARPAV有更高的相關(guān)系數(shù),并且即使在高覆蓋度情況下,線性關(guān)系仍舊顯著,而在此情況下,NDVI與FPARPAV開始表現(xiàn)為非線性關(guān)系。

    圖6 冠層總FPAR與植被指數(shù)的關(guān)系Fig.6 The relationship between total FPAR and vegetation indices

    圖7 冠層FPARPAV與植被指數(shù)的關(guān)系Fig.7 The relationship between FPARPAV and vegetation indices

    3 結(jié)論與討論

    3.1 PAV對冠層FPAR的貢獻(xiàn)

    森林冠層由光合組分PAV和非光合組分NPV組成,冠層總FPAR包含了冠層PAV和NPV的貢獻(xiàn)。PAV吸收的PAV決定了植被冠層可進(jìn)行光合作用的能量, 冠層PAV對冠層總FPAR的貢獻(xiàn)與冠層結(jié)構(gòu)有關(guān)。落葉闊葉林冠層有顯著的季節(jié)變化,冠層PAI和組分比例年內(nèi)波動明顯。在春季,葉片初長,其數(shù)量和葉片面積都較小,此時PAI以及PAV占PAI的比重都較低;夏季是闊葉林的生長旺盛期,冠層PAI和PAV面積高且穩(wěn)定,通常PAI保持在4以上;在秋季,隨著生長季逐漸結(jié)束,落葉闊葉林葉片開始枯落,PAI和PAV逐漸降低[42-43]。冠層結(jié)構(gòu)的這種季節(jié)變動極大地影響了FPARPAV占冠層總FPAR的比例?;谀M結(jié)果發(fā)現(xiàn),在低PAI(0—1)和低PAV面積比(0.5)時, FPARPAV甚至不及冠層總FPAR的一半。因此,在初春和秋末冬初,由于較低的葉面積指數(shù)(落葉林通常在冬季處于休眠狀態(tài),光合作用可忽略)和PAV面積比,冠層吸收的FPAR有很大一部分來自NPV,能有效用于光合作用的能量比重較低。而在闊葉林生長旺盛的夏季,PAI和PAV面積通常很高[19],此時,FPARPAV占冠層總FPAR的比重不低于1%,NPV對冠層FPAR的貢獻(xiàn)低。所以,在春初和秋末冬初,FPARPAV與總FPAR差距較大,以冠層總FPAR估算GPP必將產(chǎn)生較大誤差,導(dǎo)致GPP被高估;但在夏季森林植被茂密的情況下,NPV對冠層吸收的PAR貢獻(xiàn)低,FPARPAV接近總FPAR,以總FPAR代替有效FPAR估算GPP,盡管也會導(dǎo)致結(jié)果偏高但偏差較小。

    3.2 綠色FPAR(FPARPAV)與植被指數(shù)關(guān)系

    冠層結(jié)構(gòu)的變化會對冠層光譜產(chǎn)生影響,進(jìn)而影響到NDVI和EVI。兩種植被指數(shù)都隨著PAI的增加而增加,但NPV的增加會導(dǎo)致植被指數(shù)值降低。根據(jù)光譜模擬結(jié)果可知,隨著PAI的增加,冠層近紅外波段反射率增加而紅波段和藍(lán)波段反射率降低,而兩種植被指數(shù)主要由近紅外波段與兩個可見光波段的差值來決定,因此,它們必然會與PAI有正相關(guān)關(guān)系。NPV和PAV比例的變動對近紅外波段反射率影響明顯,表現(xiàn)為PAV面積比的增加提高近紅外反射率,而紅波段和藍(lán)波段反射率對PAV比例的變化不敏感。植被冠層的這種波譜特征可從冠層組分與光的相互作用來解釋。在可見光的藍(lán)、紅波段,健康植被組分,無論是PAV還是NPV,主要表現(xiàn)為吸收,反射和透射都比較低,導(dǎo)致冠層可見光部分的反射能力受組分比例的影響小。而在近紅外波段,冠層葉片透射性強(qiáng),光子可最大程度上與NPV相互作用,NPV并不產(chǎn)生透射[44-45],這使得近紅外反射率的幅度降低。

    用植被指數(shù)反演冠層FPAR是當(dāng)前FPAR研究的一個重要方法,但是反演值多為冠層總FPAR?;谠O(shè)定情景下的模擬數(shù)據(jù)發(fā)現(xiàn),無論是冠層總FPAR還是FPARPAV,它們與NDVI、EVI之間都存在顯著的線性相關(guān)關(guān)系,似乎基于NDVI和EVI可實現(xiàn)FPARPAV的估算。但是應(yīng)該注意的是,這種關(guān)系是建立在假定NPV面積比例保持恒定不變的基礎(chǔ)上,其可靠性和穩(wěn)定性需要更多的數(shù)據(jù)來證明。現(xiàn)實中,落葉闊葉林PAV占PAI的比重有明顯的季節(jié)性,在春季,它隨著葉片的生長和增加逐漸增加,夏季達(dá)到最低值,秋季由于葉片的逐漸枯落,PAV面積比又開始逐漸降低。對于成熟的落葉闊葉林來說,從春季到夏季,其PAI的增加和降低主要來自葉片的生長,而從夏季到秋季,冠層PAI 的降低同樣主要由葉片的枯落造成[46]。盡管NPV部分也將隨著葉片的生長也會擴(kuò)展,但其生長的速度和增加量難以和葉片相比,因此,將落葉闊葉林的NPV的面積指數(shù)視為定值是合理的。以固定的NPV面積指數(shù)值模擬冠層波譜和FPAR,建立模擬數(shù)據(jù)集。在冬季,落葉闊葉林冠層葉片已經(jīng)枯落,僅有NPV部分,可用該階段的PAI代表冠層NPV面積指數(shù)。根據(jù)黃玫和季勁鈞的研究,落葉闊葉林冬季植被面積指數(shù)大約在0.5[42],基于此我們對數(shù)據(jù)重新模擬,冠層PAI被設(shè)定為1到7變化,其他參數(shù)保持不變。

    圖8 落葉闊葉林冠層植被指數(shù)與綠色FPAR的關(guān)系 Fig.8 The relationship between vegetation indices and canopy FPARPAV

    結(jié)果如圖8所示,無論是NDVI還是EVI,它們與FPARPAV都存在顯著相關(guān)關(guān)系,但EVI與FPARPAV有更高的相關(guān)系數(shù)(R2=0.99),并且在高植被覆蓋情況下,NDVI與FPARPAV的散點圖呈現(xiàn)非線性關(guān)系。因此,相比于NDVI, EVI能更好描述冠層FPARPAV,更適于估算冠層FPARPAV,這種結(jié)果從理論上支持了Xiao的假設(shè)[5]。

    4 結(jié)論

    本研究通過設(shè)置不同情景,應(yīng)用物理模型SAIL模型,采用分層的方法開展冠層FPAR模擬,包括冠層總FPAR、光合組分FPAR以及非光合組分FAPR,并分析它們與NDVI、EVI的相關(guān)關(guān)系。得到以下結(jié)論:

    (1)落葉闊葉林冠層FPAR受冠層結(jié)構(gòu)影響顯著,冠層綠色FPAR隨PAI和PAV面積比的增加而增加。在低植被覆蓋區(qū),冠層綠色FPAR所占比重低,非光合組分的存在對冠層FPAR的影響大于高植被覆蓋區(qū)。

    (2)落葉闊葉林冠層FPAR與植被指數(shù)之間存在顯著的相關(guān)關(guān)系,相比于NDVI, EVI能更好描述冠層綠色FPAR。

    在模擬時,對冠層做了一定的假設(shè),且沒有考慮背景變化的影響。在高覆蓋度地區(qū),光子穿透冠層到達(dá)地面的數(shù)量較少,背景的變化對冠層FPAR和反射光譜影響比較小,但是對于低覆蓋區(qū),太陽光經(jīng)過冠層到達(dá)地面的比例會較高,地面反射率的高低會影響再次進(jìn)入冠層的光子量,因此背景的差異會引起冠層FPAR和光譜的變動。

    致謝:感謝ORNL DAAC提供冠層相關(guān)數(shù)據(jù)。

    [1] Cramer W, Kicklighter D W, Bondeau A, Iii B M, Churkina G, Nemry B, Ruimy A, Schloss A. Comparing global models of terrestrial net primary productivity (NPP): Overview and key results. Global Change Biology, 1999, 5(S1):1- 15.

    [2] Fensholt R, Sandholt I, Rasmussen M S. Evaluation of MODIS LAI, fAPAR and the relation between fAPAR and NDVI in a semi-arid environment using in situ measurements.Remote Sensing of Environment, 2004, 91(3/4): 490- 507.

    [3] 董泰鋒, 蒙繼華, 吳炳方.基于遙感的光合有效輻射吸收比率(FPAR)估算方法綜述.生態(tài)學(xué)報, 2012,32(22): 7190- 7201.

    [4] 侯湖平. 基于遙感的煤礦區(qū)植被凈初級生產(chǎn)力變化的監(jiān)測與評價.徐州:中國礦業(yè)大學(xué),2010:59

    [5] Xiao X M, Hollinger D, Aber J, Goltz M, Davidson E A, Zhang Q Y, Moore III. Satellite-based modeling of gross primary production in an evergreen needleleaf forest. Remote Sensing of Environment, 2004, 89(4): 519- 534.

    [6] Schoettker B, Phinn S, Schmidt M. How does the global Moderate Resolution Imaging Spectroradiometer (MODIS) fraction of photosynthetically active radiation (FPAR) product relate to regionally developed land cover and vegetation products in a semi-arid Australian savanna?. Journal of Applied Remote Sensing, 2010, 4(1):043538.

    [7] Zhang Q Y, Xiao X M, Braswell B, Linder E, Baret F, Moore III B. Estimating light absorption by chlorophyll, leaf and canopy in a deciduous broadleaf forest using MODIS data and a radiative transfer model. Remote Sensing of Environment, 2005, 99(3): 357- 371.

    [8] Serrano L, Gamon J A, Peuelas J. Estimation of canopy photosynthetic and nonphotosynthetic components from spectral transmittance. Ecology, 2000, 81(11): 3149- 3162

    [9] Goward S N, Huemmrich K F. Vegetation canopy PAR absorptance and the normalized difference vegetation index: an assessment using the SAIL model. Remote Sensing of Environment,1992,39(2):119- 140.

    [10] Xiao X M. Light Absorption by leaf chlorophyll and maximum light use efficiency. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(7): 1933- 1935

    [11] 王文杰,祖元剛,王慧梅. 林木非同化器官樹枝(干)光合功能研究進(jìn)展.生態(tài)學(xué)報,2007,27(4):1583- 1595.

    [12] 許慧男.樹枝與葉片綠色組織C4酶相關(guān)光合特性差異.哈爾濱: 東北林業(yè)大學(xué), 2010:1

    [13] Hall F G, Huemmrich K F, Goetz S J,Sellers P J, Nickeson J E.. Satellite remote sensing of surface energy balance: Success, failures, and unresolved issues in FIFE. Journal of Geophysical Research: Atmospheres 1992;97(D17):19061- 19089.

    [14] 高彥華,陳良富,柳欽火, 顧行發(fā), 田國良. 葉綠素吸收的光合有效輻射比率的遙感估算模型研究.遙感學(xué)報, 2006, 10(5): 798- 803

    [15] Xiao X M, Zhang Q Y, Hollinger D, Aber J, Moore B. Modeling gross primary production of an evergreen needleleaf forest using MODIS and climate data. Ecological Applications, 2005, 15(3): 954- 969.

    [16] 李剛.呼倫貝爾溫帶草地FPAR/LAI遙感估算方法研究.北京:中國農(nóng)業(yè)科學(xué)院,2009:7.

    [17] 趙英時. 遙感應(yīng)用分析原理與方法.北京:科學(xué)出版社, 2003:317

    [18] Verhoef W. Light scattering by leaf layers with application to canopy reflectance modeling: the SAIL model. Remote Sensing of Environment, 1984, 16(2): 125- 141.

    [19] Jacquemoud S, Verhoef W, Baret F, Bacour C, Zarco-Tejada P J, Asner G P, Fran?ois C, Ustin S L. PROSPECT + SAIL models: A review of use for vegetation characterization. Remote Sensing of Environment, 2009, 113(S1):S56-S66

    [20] Huemmrich K F, Goward S N. Vegetation canopy PAR absorptance and NDVI: An assessment for ten tree species with the SAIL model. Remote Sensing of Environment, 1997, 61(2): 254- 269.

    [21] Gobron N, Pinty B, Aussedat O, Chen J M, Cohen W B, Fensholt R, Gond V, Huemmrich K F, Lavergne T, Mélin F, Privette J L, Sandholt I, Taberner M, Turner D P, Verstraete M M, Widlowski J L. Evaluation of fraction of absorbed photosynthetically active radiation products for different canopy radiation transfer regimes: Methodology and results using Joint Research Center products derived from SeaWiFS against ground-based estimations. Journal of Geophysical Research: Atmospheres, 2006, 111(D13): D13110

    [22] Widlowski J L, Pinty B, Lavergne T, Verstraete M M, Gobron N. Using 1-D models to interpret the reflectance anisotropy of 3-D canopy targets: Issues and caveats. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(9):2008- 2017.

    [23] Bacour C, Baret F, Béal, D, Weiss M, Pavageau K. Neural network estimation of LAI, fAPAR, fCover and LAI × Cab, from top of canopy MERIS reflectance data: Principles and validation. Remote Sensing of Environment, 2006,105(4):313- 325.

    [24] Weiss M, Baret F, Garrigues S, Lacaze R. LAI and fAPAR CYCLOPES global products derived from VEGETATION. part 2: Validation and comparison with MODIS Collection 4 products. Remote Sensing of Environment, 2007, 110(3): 317- 331.

    [25] 肖艷芳,周德民,趙文吉.輻射傳輸模型多尺度反演植被理化參數(shù)研究進(jìn)展. 生態(tài)學(xué)報, 2013,33(11):3291- 3297.

    [26] Laurent V C E, Verhoef W, Clevers J G P W, Schaepman M E. Inversion of a coupled canopy-atmosphere model using multiangular top-of-atmosphere radiance data: A forest case study.Remote Sensing of Environment, 2011, 115(10): 2603-2612.

    [27] Li W, Fang H L. Estimation of direct, diffuse, and total FPARs from Landsat surface reflectance data and ground-based estimates over six FLUXNET sites. Journal of Geophysical Research: Biogeosciences, 2015,120(1): 2014JG002754.

    [28] Hall F G., Huemmrich K F, Strebel D E, Goetz S J,Nickeson J E, Woods K D. Biophysical, morphological, canopy optical property, and productivity data from the Superior National Forest. NASA, Technical Memorandum TM- 104568. (1992)

    [29] Dungan J. Field Spectroradiometer Calibration Progress Report, Interoffice Memorandum, NASA/Ames Research Center, 5 January 1991

    [30] Ruimy A, Saugier B, Dedieu G. Methodology for the estimation of terrestrial net primary production from remotely sensed data.Journal of Geophysical Research.1994, 99(D3): 5263- 5283

    [31] Wiegand C T, Maas S J, Aase J K, Hatfield J L, Pinter Jr P J, Jackson R D, Kanemasu E T, Lapitan R L. Multisite analyses of spectral-biophysical data for wheat. Remote Sensing of Environment, 1992, 42(1): 1- 21.

    [32] Running S W, Nemani R R, Heinsch F A, Zhao M S, Reeves M, Hashimoto H. A continental zatellite-derived measure of global terrestrial primary production. Bioscicnce,2004,54(6):547- 560.

    [33] Steinberg Daniel C, Goetz Steinberg S J G, Hyer E J. Validation of MODIS F/sub PAR/ products in Boreal forests of Alaska. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(7): 1818- 1828.

    [34] Friedl M A, Davis F W, Michaelsen J, Moritz M A. Scaling and uncertainty in the relationship between the NDVI and land surface biophysical variables: an analysis using a scene simulation model and data from FIFE. Remote Sensing of Environment, 1995, 54(3): 233- 246.

    [35] Myneni R B, Williams D L. On the relationship between FAPAR and NDVI. Remote Sensing of Environment,1994,49(3): 200- 211.

    [36] 陳雪洋,蒙繼華,吳炳方, 朱建軍, 杜鑫, 張飛飛, 紐立明. 基于HJ- 1CCD 的夏玉米FPAR 遙感監(jiān)測模型.農(nóng)業(yè)工程學(xué)報,2010,26 (S1): 241- 245.

    [37] King D A, Turner D P, Ritts W D. Parameterization of a diagnostic carbon cycle model for continental scale application. Remote Sensing of Environment, 2011, 115(7): 1653- 1664

    [38] Nakaji T, Ide R, Oguma H, Saigusa N, Fujinuma Y. Utility of spectral vegetation index for estimation of gross CO2flux under varied sky conditions. Remote Sensing of Environment, 2007, 109(3): 274- 284.

    [39] Roujean J L,Breon F M. Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote Sensing of Environment, 1995, 51(3):375- 384.

    [40] Gitelson A A, Kaufman Y J, Merzlyak M N. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sensing of Environment, 1996, 58(3): 289- 298.

    [41] Huete A R, Liu H Q, Batchily K, van Leeuwen W. A comparison of vegetation indices over a global set of TM images for EOS-MODIS. Remote Sensing of Environment, 1997, 59(3):440- 451.

    [42] 黃玫,季勁鈞.中國區(qū)域植被葉面積指數(shù)時空分布—機(jī)理模型模擬與遙感反演比較.生態(tài)學(xué)報,2010, 30(11): 3057- 3064.

    [43] 柳藝博,居為民,陳鏡明, 朱高龍, 邢白靈, 朱敬芳, 何明珠. 2000—2010年中國森林葉面積指數(shù)時空變化特征.科學(xué)通報,2012,57(16): 1435- 1445.

    [44] van Leeuwen W J D, Huete A R. Effects of standing litter on the biophysical interpretation of plant canopies with spectral indices. Remote Sensing of Environment, 1996, 55(2): 123- 138.

    [45] Asner G P, Wessman C A and Archer S. Scale dependence of absorption of photosynthetically active radiation in terrestrial ecosystems. Ecological Applications,1998, 8(4): 1003- 1021.

    [46] Running S W, Thornton P E, Nemani R, Glassy J M. Global terrestrial gross and net primary productivity from the earth observing system. New York: Springer-Verlag, 2000, 44- 57.

    Simulation and analysis on green fraction of absorbed photosynthetically active radiation of deciduous broadleaved forests canopy through remote sensing model

    LIANG Shouzhen1,2,*, SUI Xueyan1,HOU Xuehui1, WANG Meng1, YAO Huimin1, MA Wandong3

    1ShandongInstituteofAgriculturalSustainableDevelopment,Ji′nan250100,China2ShenzhenInstituteofAdvancedTechnology,ChineseAcademyofSciences,Shenzhen518055,China3SatelliteEnvironmentCenter,MinistryofEnvironmentalProtection,Beijing100094,China

    The fraction of absorbed photosynthetically active radiation (FPAR) is an important biophysical variable, widely used in satellite-based production efficiency models to estimate global primary production(GPP) or net primary production. To the forest, the vegetation canopy is composed primarily of photosynthetically active vegetation (PAV; green leaf) and non-photosynthetic vegetation (NPV; e.g., branches and stems). They absorb photosynthetically active radiation(PAR) from the sun, but only the PAR absorbed by PAV is used for photosynthesis. Remote sensing-driven biogeochemical models that use green FPAR(FPARPAV) in estimating GPP are more likely to be consistent with plant photosynthesis processes. However, there are no field and laboratory experiments to measure green FPAR at the canopy level. In this study, a method based on the radiative transfer model was used, and the scattering by arbitrary inclined leaves (SAIL) model was modified to classify forest canopy, and to model the canopy spectral reflectance and the PAR absorbed by PAV and NPV of three deciduous broadleaved forests in different scenarios (varied plant area index and leaf area index). Green FPAR of the canopy was calculated based on the PAR absorbed by PAV and total PAR, and the characteristics of the green FPAR and the relation between green FPAR and two vegetation indices (normalize difference vegetation index [NDVI] and enhanced vegetation index [EVI]) were analyzed. Our results showed that the variation in canopy structure influenced the canopy absorption of PAR, and that green FPAR was related to the plant area index and the ratio of leaf area index to plant area index. In high-coverage forests, the green FPAR was close to the total FPAR of the vegetation canopy(>80%) and the contribution of NPV to the total FPAR is very low. In contrast, in low-coverage forests or open forests, the difference between green FPAR and the total FPAR was large. Especially when the proportion of the leaf area to plant area low(0.5), the ratio of the green FPAR to the total FPAR was 45.86%. The significant relationships were found between two vegetation index (NDVI and EVI) and green FPAR, but compared with NDVI, EVI was more suitable for describing the variation in the green FPAR.

    FPAR; photosynthetically active vegetation; NDVI; EVI; deciduous broadleaved forests

    國家自然科學(xué)基金(41401407); 山東省農(nóng)業(yè)重大應(yīng)用技術(shù)創(chuàng)新項目(1-05-04)

    2016- 03- 08; 網(wǎng)絡(luò)出版日期:2017- 02- 17

    10.5846/stxb201603080405

    *通訊作者Corresponding author.E-mail: szliang_cas@163.com

    梁守真, 隋學(xué)艷,侯學(xué)會, 王猛, 姚慧敏,馬萬棟.落葉闊葉林冠層光合有效輻射分量的遙感模擬與分析.生態(tài)學(xué)報,2017,37(10):3415- 3424.

    Liang S Z, Sui X Y,Hou X H, Wang M, Yao H M, Ma W D.Simulation and analysis on green fraction of absorbed photosynthetically active radiation of deciduous broadleaved forests canopy through remote sensing model.Acta Ecologica Sinica,2017,37(10):3415- 3424.

    猜你喜歡
    闊葉林植被指數(shù)冠層
    金沙江上游與其他流域典型植被碳含量差異性探究
    基于低空遙感的果樹冠層信息提取方法研究
    基于激光雷達(dá)的樹形靶標(biāo)冠層葉面積探測模型研究
    安徽省淮南森林冠層輻射傳輸過程的特征
    亞熱帶常綠闊葉林的世界之窗 錢江源,探路國家公園
    綠色中國(2019年18期)2020-01-04 01:57:08
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    闊葉林培育及保護(hù)對策
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    主要植被指數(shù)在生態(tài)環(huán)評中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    欧美性长视频在线观看| 国产免费男女视频| 又紧又爽又黄一区二区| 丁香欧美五月| 久久人妻福利社区极品人妻图片| e午夜精品久久久久久久| 亚洲精品久久午夜乱码| 妹子高潮喷水视频| 美女高潮喷水抽搐中文字幕| 欧美亚洲 丝袜 人妻 在线| 亚洲第一欧美日韩一区二区三区| 亚洲国产中文字幕在线视频| 日韩欧美一区视频在线观看| 国产淫语在线视频| 少妇被粗大的猛进出69影院| 91麻豆精品激情在线观看国产 | 纯流量卡能插随身wifi吗| 十八禁高潮呻吟视频| 美女午夜性视频免费| 宅男免费午夜| 成年人免费黄色播放视频| 欧美日韩亚洲综合一区二区三区_| 高清黄色对白视频在线免费看| 这个男人来自地球电影免费观看| 视频区图区小说| 亚洲avbb在线观看| xxx96com| 亚洲国产欧美日韩在线播放| 国产免费现黄频在线看| 日本wwww免费看| 久久国产亚洲av麻豆专区| 黄色 视频免费看| 日韩欧美三级三区| 久久久国产成人免费| 久久久久久久久免费视频了| 999久久久国产精品视频| 久久久久精品人妻al黑| 国产极品粉嫩免费观看在线| 搡老熟女国产l中国老女人| 日韩欧美国产一区二区入口| 精品国产乱码久久久久久男人| 一区在线观看完整版| 亚洲在线自拍视频| 变态另类成人亚洲欧美熟女 | 成人国产一区最新在线观看| 黄色视频,在线免费观看| 亚洲人成伊人成综合网2020| 一区二区三区国产精品乱码| videos熟女内射| 精品国产一区二区三区久久久樱花| 18禁黄网站禁片午夜丰满| 欧美最黄视频在线播放免费 | 波多野结衣av一区二区av| 老司机影院毛片| 99久久精品国产亚洲精品| 一级毛片高清免费大全| 人人妻,人人澡人人爽秒播| 亚洲av日韩精品久久久久久密| 别揉我奶头~嗯~啊~动态视频| 黄色a级毛片大全视频| 丰满饥渴人妻一区二区三| 欧美激情 高清一区二区三区| 18禁裸乳无遮挡免费网站照片 | 啦啦啦免费观看视频1| 精品国产乱码久久久久久男人| 在线av久久热| 成年人午夜在线观看视频| 亚洲国产看品久久| 亚洲国产欧美一区二区综合| 亚洲欧美日韩高清在线视频| 欧美+亚洲+日韩+国产| 人人妻人人澡人人爽人人夜夜| bbb黄色大片| 美女福利国产在线| 国产成人系列免费观看| 亚洲国产精品一区二区三区在线| 正在播放国产对白刺激| 每晚都被弄得嗷嗷叫到高潮| 国产三级黄色录像| 亚洲成国产人片在线观看| 亚洲国产欧美网| 国产一区二区三区视频了| 国产男靠女视频免费网站| 日韩欧美三级三区| 免费在线观看黄色视频的| 免费一级毛片在线播放高清视频 | 在线观看www视频免费| 欧美激情 高清一区二区三区| 欧美成狂野欧美在线观看| 高清av免费在线| 国产精品永久免费网站| 午夜福利在线免费观看网站| 精品免费久久久久久久清纯 | 久久人妻av系列| 亚洲成人手机| 久久国产精品影院| 午夜免费成人在线视频| 在线观看日韩欧美| 丰满人妻熟妇乱又伦精品不卡| 中文字幕高清在线视频| 亚洲精品成人av观看孕妇| 好男人电影高清在线观看| a级片在线免费高清观看视频| 日日爽夜夜爽网站| 精品人妻1区二区| 欧美日韩乱码在线| 制服人妻中文乱码| 乱人伦中国视频| 欧美日韩乱码在线| 久久人妻av系列| 欧美日韩亚洲高清精品| 天堂中文最新版在线下载| 高清av免费在线| 在线观看66精品国产| 亚洲第一青青草原| 两人在一起打扑克的视频| 两性午夜刺激爽爽歪歪视频在线观看 | 一边摸一边做爽爽视频免费| 成人免费观看视频高清| 欧美丝袜亚洲另类 | 欧美日韩国产mv在线观看视频| 久久国产亚洲av麻豆专区| 欧美 日韩 精品 国产| 国产亚洲欧美在线一区二区| 99热网站在线观看| 亚洲自偷自拍图片 自拍| 精品一区二区三卡| 两性夫妻黄色片| 久久国产亚洲av麻豆专区| 午夜影院日韩av| 高清黄色对白视频在线免费看| 高清在线国产一区| 好男人电影高清在线观看| 极品教师在线免费播放| 色综合欧美亚洲国产小说| 变态另类成人亚洲欧美熟女 | 亚洲精品久久成人aⅴ小说| 这个男人来自地球电影免费观看| 午夜日韩欧美国产| 黄色女人牲交| tocl精华| 女性被躁到高潮视频| 日韩欧美一区二区三区在线观看 | 在线观看66精品国产| 亚洲精品美女久久久久99蜜臀| 久久久精品国产亚洲av高清涩受| 久热爱精品视频在线9| 十八禁网站免费在线| aaaaa片日本免费| 老司机亚洲免费影院| 日韩欧美一区二区三区在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 高清视频免费观看一区二区| 亚洲第一欧美日韩一区二区三区| 国产亚洲av高清不卡| 国产一区二区三区综合在线观看| 女人精品久久久久毛片| 999久久久精品免费观看国产| 国产精品九九99| 亚洲成国产人片在线观看| 又紧又爽又黄一区二区| 精品国产乱码久久久久久男人| 国产精品永久免费网站| 新久久久久国产一级毛片| cao死你这个sao货| 精品卡一卡二卡四卡免费| 午夜福利乱码中文字幕| 国产精品秋霞免费鲁丝片| 欧美精品啪啪一区二区三区| 男女床上黄色一级片免费看| 久久久国产一区二区| 一区二区三区激情视频| 亚洲精品中文字幕一二三四区| 欧美精品av麻豆av| 国产精品影院久久| 麻豆av在线久日| 国产激情欧美一区二区| 成年动漫av网址| 成人黄色视频免费在线看| 精品国产国语对白av| 国产淫语在线视频| 日本vs欧美在线观看视频| 国产精品98久久久久久宅男小说| 亚洲精品国产精品久久久不卡| 国产成人欧美在线观看 | 久99久视频精品免费| 国产免费男女视频| 亚洲熟妇熟女久久| 高清在线国产一区| 天天躁夜夜躁狠狠躁躁| а√天堂www在线а√下载 | а√天堂www在线а√下载 | 精品免费久久久久久久清纯 | 久久中文字幕一级| 一级毛片女人18水好多| 国产主播在线观看一区二区| 99久久精品国产亚洲精品| 变态另类成人亚洲欧美熟女 | a级片在线免费高清观看视频| 无人区码免费观看不卡| 激情视频va一区二区三区| 久久久国产精品麻豆| 高清av免费在线| 丝袜在线中文字幕| 精品福利永久在线观看| 日日摸夜夜添夜夜添小说| 精品国产一区二区三区四区第35| 国产精品秋霞免费鲁丝片| 黄色成人免费大全| 国产有黄有色有爽视频| 脱女人内裤的视频| av欧美777| bbb黄色大片| 国产精品免费视频内射| 精品免费久久久久久久清纯 | 国产免费av片在线观看野外av| 搡老岳熟女国产| 下体分泌物呈黄色| 三上悠亚av全集在线观看| 亚洲国产看品久久| 久久人妻熟女aⅴ| aaaaa片日本免费| 亚洲av成人一区二区三| 大型黄色视频在线免费观看| 亚洲色图av天堂| 国产xxxxx性猛交| 国产亚洲精品久久久久5区| 性少妇av在线| 精品国产国语对白av| 91字幕亚洲| 成在线人永久免费视频| 在线av久久热| av不卡在线播放| 久久国产乱子伦精品免费另类| 中文字幕另类日韩欧美亚洲嫩草| 女人精品久久久久毛片| 久久久久久亚洲精品国产蜜桃av| 一边摸一边抽搐一进一小说 | 动漫黄色视频在线观看| 五月开心婷婷网| e午夜精品久久久久久久| 老司机深夜福利视频在线观看| 一a级毛片在线观看| 久久久久国内视频| 国产精品久久电影中文字幕 | 国产精品久久久人人做人人爽| 黄片大片在线免费观看| 一级毛片高清免费大全| 国产真人三级小视频在线观看| 精品久久久久久久毛片微露脸| 久久久久久亚洲精品国产蜜桃av| 久久久久国产精品人妻aⅴ院 | 大片电影免费在线观看免费| 老熟妇乱子伦视频在线观看| 亚洲成人免费av在线播放| 亚洲性夜色夜夜综合| 多毛熟女@视频| av中文乱码字幕在线| 欧美精品一区二区免费开放| 97人妻天天添夜夜摸| 一个人免费在线观看的高清视频| 亚洲色图综合在线观看| 午夜福利,免费看| 99久久综合精品五月天人人| 国产亚洲精品第一综合不卡| 美女 人体艺术 gogo| 亚洲中文字幕日韩| 日韩人妻精品一区2区三区| 精品国产国语对白av| 午夜福利,免费看| 高清av免费在线| 精品人妻1区二区| 麻豆国产av国片精品| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 一级毛片精品| 国产欧美日韩精品亚洲av| 国产精品永久免费网站| 午夜福利乱码中文字幕| 可以免费在线观看a视频的电影网站| 超碰97精品在线观看| 大型av网站在线播放| 午夜91福利影院| 十八禁网站免费在线| www.熟女人妻精品国产| 精品人妻在线不人妻| 日韩一卡2卡3卡4卡2021年| 999久久久精品免费观看国产| 国产精品 国内视频| 一本大道久久a久久精品| 亚洲精品粉嫩美女一区| 国产高清国产精品国产三级| 自线自在国产av| 一级片免费观看大全| 91精品国产国语对白视频| 女人高潮潮喷娇喘18禁视频| tocl精华| 亚洲精品中文字幕在线视频| 中文欧美无线码| 水蜜桃什么品种好| 亚洲黑人精品在线| 久久久久国产精品人妻aⅴ院 | 久久国产精品影院| 大型av网站在线播放| 女人爽到高潮嗷嗷叫在线视频| 大香蕉久久网| 一级毛片精品| 国产一区二区激情短视频| 最近最新中文字幕大全电影3 | 国产97色在线日韩免费| 亚洲伊人色综图| 麻豆国产av国片精品| 黄网站色视频无遮挡免费观看| 母亲3免费完整高清在线观看| 欧美人与性动交α欧美软件| 少妇裸体淫交视频免费看高清 | 国产成人免费观看mmmm| 乱人伦中国视频| 中文字幕人妻熟女乱码| 久久久久久人人人人人| 久久午夜亚洲精品久久| 超碰97精品在线观看| 丝袜人妻中文字幕| 在线观看66精品国产| 国产成人精品无人区| 国产免费男女视频| 性色av乱码一区二区三区2| 热re99久久国产66热| 中文字幕人妻丝袜一区二区| 国产成人av激情在线播放| 久久亚洲精品不卡| 国产精品 欧美亚洲| 国产真人三级小视频在线观看| 免费不卡黄色视频| 久久久久久人人人人人| 国产欧美日韩综合在线一区二区| 精品电影一区二区在线| 麻豆国产av国片精品| 侵犯人妻中文字幕一二三四区| 国产成人精品无人区| 久久精品亚洲熟妇少妇任你| 性色av乱码一区二区三区2| 美国免费a级毛片| 午夜福利免费观看在线| 日日爽夜夜爽网站| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 亚洲精品美女久久久久99蜜臀| 亚洲情色 制服丝袜| 欧美午夜高清在线| 搡老乐熟女国产| 精品人妻在线不人妻| 国产男女内射视频| 一区福利在线观看| 久久久久精品国产欧美久久久| 午夜福利欧美成人| 美女扒开内裤让男人捅视频| 精品一品国产午夜福利视频| 可以免费在线观看a视频的电影网站| 91成人精品电影| 日本撒尿小便嘘嘘汇集6| 伊人久久大香线蕉亚洲五| 在线观看日韩欧美| 国产精品综合久久久久久久免费 | 美女高潮喷水抽搐中文字幕| 国产日韩一区二区三区精品不卡| 欧美黑人欧美精品刺激| 亚洲精品国产色婷婷电影| 悠悠久久av| 久久国产精品男人的天堂亚洲| 中文字幕制服av| 久久国产精品大桥未久av| 国产精品永久免费网站| 国产不卡av网站在线观看| 国产高清videossex| 午夜福利一区二区在线看| 亚洲avbb在线观看| 91麻豆精品激情在线观看国产 | 99久久99久久久精品蜜桃| av片东京热男人的天堂| 极品少妇高潮喷水抽搐| 成年女人毛片免费观看观看9 | 亚洲色图 男人天堂 中文字幕| 国产99白浆流出| 久久午夜亚洲精品久久| 三上悠亚av全集在线观看| 国产欧美日韩综合在线一区二区| 亚洲伊人色综图| 女人精品久久久久毛片| 国产一区在线观看成人免费| 夫妻午夜视频| 美女高潮喷水抽搐中文字幕| 日日爽夜夜爽网站| 99国产综合亚洲精品| 国精品久久久久久国模美| av有码第一页| av电影中文网址| 狠狠婷婷综合久久久久久88av| 亚洲一区二区三区不卡视频| 男女免费视频国产| 亚洲五月色婷婷综合| 男人的好看免费观看在线视频 | 亚洲av美国av| 国产激情欧美一区二区| 国产精品自产拍在线观看55亚洲 | 国产成人精品无人区| 免费日韩欧美在线观看| 久久久久久久午夜电影 | 国产男女内射视频| 涩涩av久久男人的天堂| 欧美人与性动交α欧美精品济南到| 欧美 亚洲 国产 日韩一| 国产主播在线观看一区二区| 99久久国产精品久久久| 国产亚洲av高清不卡| 亚洲av熟女| 中文字幕高清在线视频| 日韩欧美三级三区| 国产aⅴ精品一区二区三区波| 高清视频免费观看一区二区| 真人做人爱边吃奶动态| 久久人妻熟女aⅴ| 欧美日韩精品网址| 在线观看午夜福利视频| 丝袜在线中文字幕| 亚洲欧美激情综合另类| 一a级毛片在线观看| 国产精品免费大片| 99久久99久久久精品蜜桃| 色94色欧美一区二区| 欧美乱色亚洲激情| a级毛片在线看网站| 一边摸一边做爽爽视频免费| 99香蕉大伊视频| 女人高潮潮喷娇喘18禁视频| 熟女少妇亚洲综合色aaa.| 他把我摸到了高潮在线观看| 日日夜夜操网爽| 国产一区二区三区在线臀色熟女 | 久久久国产欧美日韩av| 午夜久久久在线观看| 91麻豆精品激情在线观看国产 | 午夜视频精品福利| 亚洲欧美激情综合另类| 亚洲精品久久成人aⅴ小说| 91大片在线观看| 老汉色av国产亚洲站长工具| 美国免费a级毛片| 男女下面插进去视频免费观看| 国产视频一区二区在线看| 两性夫妻黄色片| 韩国精品一区二区三区| 久久草成人影院| 国产蜜桃级精品一区二区三区 | 俄罗斯特黄特色一大片| 岛国在线观看网站| 大香蕉久久成人网| av线在线观看网站| 亚洲精品国产一区二区精华液| 亚洲第一青青草原| 看黄色毛片网站| 亚洲国产精品一区二区三区在线| 亚洲五月色婷婷综合| 一区在线观看完整版| 亚洲少妇的诱惑av| 高清av免费在线| 人妻一区二区av| 黄色毛片三级朝国网站| 免费观看人在逋| 侵犯人妻中文字幕一二三四区| 欧美日韩视频精品一区| 新久久久久国产一级毛片| 亚洲精品国产一区二区精华液| 人妻一区二区av| 丝袜美足系列| 母亲3免费完整高清在线观看| 国产一区二区三区视频了| 精品无人区乱码1区二区| 欧美精品高潮呻吟av久久| 亚洲欧洲精品一区二区精品久久久| 亚洲av美国av| 国产亚洲精品一区二区www | 丝袜在线中文字幕| 国产成人欧美在线观看 | 亚洲精品av麻豆狂野| 精品久久久久久久毛片微露脸| 欧美日韩瑟瑟在线播放| 亚洲色图 男人天堂 中文字幕| 99国产极品粉嫩在线观看| 国产激情久久老熟女| 精品人妻熟女毛片av久久网站| 在线观看舔阴道视频| 母亲3免费完整高清在线观看| 身体一侧抽搐| 可以免费在线观看a视频的电影网站| 他把我摸到了高潮在线观看| 美女福利国产在线| 老鸭窝网址在线观看| 飞空精品影院首页| 精品欧美一区二区三区在线| 中文字幕人妻丝袜一区二区| 久久精品亚洲精品国产色婷小说| 超色免费av| 99久久人妻综合| 夫妻午夜视频| 女人久久www免费人成看片| 欧美国产精品va在线观看不卡| 久久精品亚洲av国产电影网| 女人高潮潮喷娇喘18禁视频| 亚洲欧洲精品一区二区精品久久久| 女性生殖器流出的白浆| 欧美日韩精品网址| 一进一出抽搐动态| 欧美黑人精品巨大| 91精品国产国语对白视频| 欧美在线一区亚洲| 欧美日韩亚洲高清精品| 欧美另类亚洲清纯唯美| 捣出白浆h1v1| 免费在线观看视频国产中文字幕亚洲| 国产精品成人在线| 91成人精品电影| 最新美女视频免费是黄的| 成在线人永久免费视频| 麻豆成人av在线观看| x7x7x7水蜜桃| 91麻豆精品激情在线观看国产 | 黑人巨大精品欧美一区二区蜜桃| 国产精品欧美亚洲77777| 女性生殖器流出的白浆| 国产精品一区二区精品视频观看| 99久久国产精品久久久| 在线永久观看黄色视频| 成人精品一区二区免费| 免费看十八禁软件| 免费看a级黄色片| 国产高清国产精品国产三级| 精品国产一区二区三区久久久樱花| 一级黄色大片毛片| 免费不卡黄色视频| 免费在线观看完整版高清| 午夜福利影视在线免费观看| 中文欧美无线码| 日韩成人在线观看一区二区三区| 黄色a级毛片大全视频| 精品免费久久久久久久清纯 | 国产成人欧美在线观看 | 美国免费a级毛片| 一区二区三区精品91| 精品亚洲成国产av| 人人妻人人添人人爽欧美一区卜| 在线看a的网站| 12—13女人毛片做爰片一| av网站在线播放免费| 91精品三级在线观看| 色婷婷久久久亚洲欧美| 日本一区二区免费在线视频| 午夜福利一区二区在线看| netflix在线观看网站| 黄片大片在线免费观看| 亚洲国产毛片av蜜桃av| 成熟少妇高潮喷水视频| 天天躁夜夜躁狠狠躁躁| 日本黄色日本黄色录像| 久久青草综合色| 天堂√8在线中文| 一级毛片高清免费大全| 黄色视频,在线免费观看| 窝窝影院91人妻| 亚洲性夜色夜夜综合| 王馨瑶露胸无遮挡在线观看| 久久人妻福利社区极品人妻图片| 欧美日韩国产mv在线观看视频| 少妇裸体淫交视频免费看高清 | 亚洲午夜精品一区,二区,三区| av网站免费在线观看视频| 亚洲成人免费电影在线观看| 久9热在线精品视频| 国产男靠女视频免费网站| 自拍欧美九色日韩亚洲蝌蚪91| 国产主播在线观看一区二区| 久久热在线av| 国产免费现黄频在线看| av不卡在线播放| 在线观看午夜福利视频| 国产成人欧美| 免费不卡黄色视频| 日韩精品免费视频一区二区三区| 亚洲欧洲精品一区二区精品久久久| 亚洲男人天堂网一区| 中文字幕高清在线视频| 日韩欧美免费精品| 精品福利观看| 国产精华一区二区三区| 久久久久久久国产电影| 久久这里只有精品19| 在线观看一区二区三区激情| 亚洲一区高清亚洲精品| 精品一品国产午夜福利视频| av国产精品久久久久影院| 免费女性裸体啪啪无遮挡网站| 日韩中文字幕欧美一区二区| 色婷婷av一区二区三区视频| 国产高清激情床上av| 丁香六月欧美| av网站免费在线观看视频| 久久婷婷成人综合色麻豆| 美女 人体艺术 gogo| tube8黄色片| 人妻久久中文字幕网| 丰满的人妻完整版| 久久国产亚洲av麻豆专区| 欧美亚洲 丝袜 人妻 在线| 欧美在线黄色| 99久久99久久久精品蜜桃|