• <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
    女人被躁到高潮嗷嗷叫费观| 99精国产麻豆久久婷婷| 国产精品三级大全| 国产不卡av网站在线观看| 青春草国产在线视频| 亚洲国产av影院在线观看| 飞空精品影院首页| 国产福利在线免费观看视频| 麻豆乱淫一区二区| av在线播放精品| freevideosex欧美| 亚洲三级黄色毛片| 热99国产精品久久久久久7| 色婷婷av一区二区三区视频| 亚洲国产精品一区二区三区在线| 波多野结衣av一区二区av| 综合色丁香网| 亚洲伊人久久精品综合| 亚洲激情五月婷婷啪啪| 九草在线视频观看| a 毛片基地| 成人18禁高潮啪啪吃奶动态图| 一区二区三区精品91| 看十八女毛片水多多多| 18禁动态无遮挡网站| 亚洲成人手机| 看非洲黑人一级黄片| 久久久久国产网址| 少妇熟女欧美另类| xxx大片免费视频| 五月开心婷婷网| 亚洲第一av免费看| 亚洲av综合色区一区| 国产日韩一区二区三区精品不卡| 国产精品成人在线| 香蕉丝袜av| 老熟女久久久| 国产免费一区二区三区四区乱码| 少妇的丰满在线观看| 国产毛片在线视频| 亚洲精品国产av蜜桃| 搡老乐熟女国产| 免费观看a级毛片全部| 亚洲,欧美,日韩| 韩国av在线不卡| 久久97久久精品| 日韩不卡一区二区三区视频在线| 午夜福利在线观看免费完整高清在| 亚洲第一av免费看| 美女国产高潮福利片在线看| 中文字幕最新亚洲高清| 亚洲精品视频女| 国产淫语在线视频| 麻豆乱淫一区二区| 国产精品久久久久成人av| 久久久国产一区二区| 日韩av在线免费看完整版不卡| 日韩欧美一区视频在线观看| 最近的中文字幕免费完整| 日韩成人av中文字幕在线观看| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 午夜福利乱码中文字幕| 欧美精品亚洲一区二区| 国产精品女同一区二区软件| freevideosex欧美| 欧美日韩综合久久久久久| 亚洲精品国产av蜜桃| 97精品久久久久久久久久精品| 一区二区三区四区激情视频| 亚洲,一卡二卡三卡| 超碰成人久久| 久久久精品区二区三区| 97在线人人人人妻| 国产淫语在线视频| 欧美人与善性xxx| 成人毛片60女人毛片免费| 日韩不卡一区二区三区视频在线| av片东京热男人的天堂| 伦理电影免费视频| 成人国产麻豆网| 一本色道久久久久久精品综合| 久久久久人妻精品一区果冻| 男女边吃奶边做爰视频| 亚洲精品久久午夜乱码| 欧美精品av麻豆av| 人人妻人人添人人爽欧美一区卜| 97在线视频观看| 国产欧美日韩综合在线一区二区| 日本-黄色视频高清免费观看| 高清视频免费观看一区二区| 日本91视频免费播放| 亚洲精品第二区| 日韩av免费高清视频| 人人妻人人澡人人看| 日韩人妻精品一区2区三区| 丝袜美腿诱惑在线| 亚洲,欧美精品.| 欧美成人午夜免费资源| 在线免费观看不下载黄p国产| 免费观看性生交大片5| 性少妇av在线| 在线观看免费日韩欧美大片| 国产高清不卡午夜福利| 亚洲精华国产精华液的使用体验| 超碰97精品在线观看| 757午夜福利合集在线观看| 在线观看免费高清a一片| 国产成人啪精品午夜网站| 99精品久久久久人妻精品| 精品福利永久在线观看| www.999成人在线观看| 黄色女人牲交| 国产成人欧美| 黄色 视频免费看| 亚洲自偷自拍图片 自拍| 窝窝影院91人妻| 天堂影院成人在线观看| 大码成人一级视频| 国产欧美日韩综合在线一区二区| 久久久国产成人免费| 亚洲欧美一区二区三区黑人| 黄色视频,在线免费观看| 岛国视频午夜一区免费看| 精品人妻1区二区| 无限看片的www在线观看| 精品一区二区三区av网在线观看| 成年女人毛片免费观看观看9| 久久人妻福利社区极品人妻图片| 亚洲av成人av| 男女做爰动态图高潮gif福利片 | 亚洲全国av大片| 国产成人免费无遮挡视频| 午夜久久久在线观看| 欧美日韩福利视频一区二区| 高清毛片免费观看视频网站 | 亚洲国产欧美一区二区综合| 女性生殖器流出的白浆| 欧美日韩亚洲综合一区二区三区_| 深夜精品福利| 国产av一区二区精品久久| 国产精品免费视频内射| 91在线观看av| 国产99久久九九免费精品| 真人做人爱边吃奶动态| 夜夜躁狠狠躁天天躁| 国产高清videossex| 久久热在线av| 久久午夜综合久久蜜桃| 亚洲av成人一区二区三| 免费在线观看视频国产中文字幕亚洲| 欧美乱码精品一区二区三区| 水蜜桃什么品种好| 在线观看免费视频日本深夜| 亚洲三区欧美一区| 亚洲午夜精品一区,二区,三区| 欧美精品啪啪一区二区三区| 激情在线观看视频在线高清| 欧美精品亚洲一区二区| 一级a爱视频在线免费观看| 一级毛片女人18水好多| 午夜成年电影在线免费观看| 久久久久国产一级毛片高清牌| 国产激情久久老熟女| 国产区一区二久久| 大香蕉久久成人网| 免费看十八禁软件| 国产精华一区二区三区| 久99久视频精品免费| 国产黄色免费在线视频| 国产精品一区二区在线不卡| 高清在线国产一区| 夜夜躁狠狠躁天天躁| 中文字幕色久视频| 成人18禁在线播放| 国产麻豆69| 三级毛片av免费| 久久精品成人免费网站| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美一区二区三区黑人| 国产午夜精品久久久久久| 岛国视频午夜一区免费看| 99精品欧美一区二区三区四区| 午夜福利免费观看在线| 国内毛片毛片毛片毛片毛片| 免费观看人在逋| 国产亚洲精品久久久久5区| 欧美日本中文国产一区发布| 90打野战视频偷拍视频| 天天添夜夜摸| 在线天堂中文资源库| 久久久久亚洲av毛片大全| 老司机靠b影院| 丝袜在线中文字幕| 亚洲欧美精品综合一区二区三区| 成年版毛片免费区| 热99国产精品久久久久久7| 不卡av一区二区三区| 久久香蕉精品热| 国产又爽黄色视频| 国产一区二区在线av高清观看| 美女高潮到喷水免费观看| 午夜福利在线观看吧| 亚洲一区高清亚洲精品| 久久影院123| 女性被躁到高潮视频| 99久久久亚洲精品蜜臀av| 一a级毛片在线观看| 亚洲av美国av| 精品国产乱码久久久久久男人| 国产精品成人在线| 亚洲一区高清亚洲精品| 淫妇啪啪啪对白视频| 一边摸一边做爽爽视频免费| 亚洲精品一区av在线观看| 久久精品国产清高在天天线| 亚洲性夜色夜夜综合| 国产精品偷伦视频观看了| 欧美成人免费av一区二区三区| 一区二区三区国产精品乱码| 午夜视频精品福利| 午夜免费观看网址| 成在线人永久免费视频| 午夜免费激情av| 夫妻午夜视频| 亚洲在线自拍视频| 在线十欧美十亚洲十日本专区| 国产99白浆流出| 国产一区在线观看成人免费| 亚洲熟妇中文字幕五十中出 | 亚洲成a人片在线一区二区| 91老司机精品| 国产精品免费视频内射| 国产又爽黄色视频| 亚洲狠狠婷婷综合久久图片| 国产精华一区二区三区| 国产成人av教育| 亚洲男人的天堂狠狠| 亚洲欧美日韩高清在线视频| 99久久99久久久精品蜜桃| 久久精品国产综合久久久| 亚洲午夜理论影院| 又黄又爽又免费观看的视频| 亚洲国产精品合色在线| 村上凉子中文字幕在线| 国产高清激情床上av| 久久久久久久午夜电影 | 国产精品久久电影中文字幕| 亚洲欧美日韩高清在线视频| 亚洲一区二区三区不卡视频| 黑人巨大精品欧美一区二区mp4| 欧美一区二区精品小视频在线| 天堂俺去俺来也www色官网| 天天躁夜夜躁狠狠躁躁| 老熟妇仑乱视频hdxx| 精品一区二区三卡| 19禁男女啪啪无遮挡网站| 国产精品爽爽va在线观看网站 | 亚洲国产欧美一区二区综合| 色精品久久人妻99蜜桃| 在线观看免费午夜福利视频| 色婷婷av一区二区三区视频| av片东京热男人的天堂| 99久久综合精品五月天人人| 淫妇啪啪啪对白视频| 国产激情欧美一区二区| 亚洲精品一卡2卡三卡4卡5卡| 一边摸一边做爽爽视频免费| 搡老岳熟女国产| 女生性感内裤真人,穿戴方法视频| 99riav亚洲国产免费| 精品国产一区二区久久| 高清在线国产一区| 看免费av毛片| 黑丝袜美女国产一区| 久久精品成人免费网站| 精品久久久久久电影网| 淫妇啪啪啪对白视频| 午夜激情av网站| 欧美日韩av久久| 757午夜福利合集在线观看| 亚洲色图av天堂| 老熟妇仑乱视频hdxx| a级片在线免费高清观看视频| 51午夜福利影视在线观看| 国产成人av教育| 黄色丝袜av网址大全| 看片在线看免费视频| 999久久久精品免费观看国产| 国产亚洲精品第一综合不卡| 成人亚洲精品一区在线观看| 十八禁人妻一区二区| 99精品久久久久人妻精品| 亚洲第一欧美日韩一区二区三区| 国产aⅴ精品一区二区三区波| 国产免费男女视频| 最好的美女福利视频网| 亚洲中文av在线| 中亚洲国语对白在线视频| 欧美成人免费av一区二区三区| 丝袜美足系列| 久久欧美精品欧美久久欧美| 亚洲片人在线观看| 日韩视频一区二区在线观看| 99riav亚洲国产免费| 91国产中文字幕| av网站在线播放免费| 久久热在线av| 在线观看66精品国产| 天天影视国产精品| 成人亚洲精品一区在线观看| 亚洲五月天丁香| 午夜精品在线福利| 男女高潮啪啪啪动态图| 欧美精品啪啪一区二区三区| 999精品在线视频| 日本欧美视频一区| 国产免费现黄频在线看| 国产精品久久视频播放| x7x7x7水蜜桃| xxx96com| 午夜两性在线视频| 国产av一区二区精品久久| 亚洲国产欧美一区二区综合| 免费不卡黄色视频| 免费在线观看日本一区| 夜夜看夜夜爽夜夜摸 | 99久久人妻综合| 一进一出好大好爽视频| 亚洲精品美女久久久久99蜜臀| 亚洲七黄色美女视频| 国产熟女午夜一区二区三区| 国产单亲对白刺激| 日韩免费高清中文字幕av| 欧美日韩亚洲高清精品| 久久精品aⅴ一区二区三区四区| 天天添夜夜摸| 国产熟女午夜一区二区三区| 亚洲精品久久午夜乱码| 亚洲va日本ⅴa欧美va伊人久久| 国产深夜福利视频在线观看| 亚洲成av片中文字幕在线观看| 精品免费久久久久久久清纯| 久久久久国产一级毛片高清牌| 窝窝影院91人妻| 老汉色∧v一级毛片| 亚洲久久久国产精品| 欧美乱妇无乱码| 50天的宝宝边吃奶边哭怎么回事| 美国免费a级毛片| 欧美一区二区精品小视频在线| 深夜精品福利| 在线观看午夜福利视频| 伦理电影免费视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美国产一区二区入口| 亚洲熟妇中文字幕五十中出 | 黑人猛操日本美女一级片| 一级毛片高清免费大全| 日韩免费av在线播放| 午夜福利在线免费观看网站| 精品无人区乱码1区二区| 丁香六月欧美| 免费在线观看影片大全网站| 岛国视频午夜一区免费看| 国产伦人伦偷精品视频| 日韩一卡2卡3卡4卡2021年| 91大片在线观看| 一本大道久久a久久精品| 欧美日韩瑟瑟在线播放| 动漫黄色视频在线观看| 午夜免费观看网址| 欧美日韩视频精品一区| 亚洲中文日韩欧美视频| 99久久久亚洲精品蜜臀av| 亚洲男人天堂网一区| 真人一进一出gif抽搐免费| av国产精品久久久久影院| 精品少妇一区二区三区视频日本电影| 在线观看免费视频网站a站| 国内毛片毛片毛片毛片毛片| 亚洲人成伊人成综合网2020| 人人澡人人妻人| 久热爱精品视频在线9| 天天躁夜夜躁狠狠躁躁| 久久久久久久午夜电影 | 一级黄色大片毛片| 美女高潮到喷水免费观看| 99久久综合精品五月天人人| 成年版毛片免费区| 精品国产美女av久久久久小说| 欧美人与性动交α欧美精品济南到| 在线国产一区二区在线| 国产不卡一卡二| 女性被躁到高潮视频| 欧美日韩中文字幕国产精品一区二区三区 | 99热国产这里只有精品6| 色播在线永久视频| 男人舔女人下体高潮全视频| 亚洲五月婷婷丁香| 一区二区三区国产精品乱码| 韩国精品一区二区三区| 日日干狠狠操夜夜爽| e午夜精品久久久久久久| 午夜激情av网站| 亚洲中文av在线| 久久草成人影院| 午夜激情av网站| 亚洲,欧美精品.| 夜夜爽天天搞| 一级片'在线观看视频| 国产精品国产高清国产av| 国产av一区二区精品久久| 后天国语完整版免费观看| 99久久99久久久精品蜜桃| 男人舔女人下体高潮全视频| 国产黄色免费在线视频| 人妻丰满熟妇av一区二区三区| 国产亚洲精品第一综合不卡| 一级片免费观看大全| 国产aⅴ精品一区二区三区波| 国产视频一区二区在线看| 一二三四社区在线视频社区8| 久久久精品国产亚洲av高清涩受| 亚洲av五月六月丁香网| 国产一卡二卡三卡精品| 久久香蕉激情| 岛国在线观看网站| 精品一区二区三卡| 男女下面插进去视频免费观看| 国产单亲对白刺激| 老司机在亚洲福利影院| 精品一区二区三卡| 两性夫妻黄色片| netflix在线观看网站| 天天影视国产精品| 嫁个100分男人电影在线观看| 99国产精品99久久久久| 国产av在哪里看| 免费女性裸体啪啪无遮挡网站| 国产精品 欧美亚洲| 国产一区二区三区视频了| 中文字幕另类日韩欧美亚洲嫩草| 国产av一区在线观看免费| 久久婷婷成人综合色麻豆| 在线看a的网站| 亚洲欧洲精品一区二区精品久久久| 久久久久久久久久久久大奶| 久久精品91无色码中文字幕| 日日摸夜夜添夜夜添小说| 精品国产一区二区三区四区第35| 夜夜爽天天搞| 免费在线观看日本一区| 精品人妻1区二区| 丝袜人妻中文字幕| 免费在线观看日本一区| 国产三级在线视频| 男人的好看免费观看在线视频 | aaaaa片日本免费| 超碰成人久久| 91精品三级在线观看| 成人国语在线视频| 久久香蕉精品热| 日韩视频一区二区在线观看| 多毛熟女@视频| 久久国产精品影院| 天堂√8在线中文| 女人精品久久久久毛片| 欧美日韩av久久| 19禁男女啪啪无遮挡网站| 成人亚洲精品av一区二区 | 国产有黄有色有爽视频| 成人18禁在线播放| 国产成人欧美| 99国产精品一区二区蜜桃av| 国产在线精品亚洲第一网站| 午夜免费成人在线视频| 性色av乱码一区二区三区2| 丰满迷人的少妇在线观看| 涩涩av久久男人的天堂| 国产成人影院久久av| 国产1区2区3区精品| 国产真人三级小视频在线观看| 日韩高清综合在线| 9191精品国产免费久久| 婷婷六月久久综合丁香| 国产亚洲欧美98| 18禁观看日本| 黄片大片在线免费观看| 不卡av一区二区三区| 99在线人妻在线中文字幕| 欧美中文综合在线视频| 一进一出抽搐gif免费好疼 | 欧美大码av| 久久人人精品亚洲av| 久9热在线精品视频| 一二三四在线观看免费中文在| 90打野战视频偷拍视频| 老司机福利观看| 99久久久亚洲精品蜜臀av| 母亲3免费完整高清在线观看| 一级毛片高清免费大全| av视频免费观看在线观看| 精品久久久久久电影网| 纯流量卡能插随身wifi吗| 亚洲精华国产精华精| 久久精品亚洲精品国产色婷小说| 国产欧美日韩精品亚洲av| 在线播放国产精品三级| 久久国产乱子伦精品免费另类| 久9热在线精品视频| av免费在线观看网站| 琪琪午夜伦伦电影理论片6080| 国产aⅴ精品一区二区三区波| 无遮挡黄片免费观看| 99在线视频只有这里精品首页| 亚洲av五月六月丁香网| 99精国产麻豆久久婷婷| 我的亚洲天堂| 免费少妇av软件| 美女 人体艺术 gogo| 亚洲欧美精品综合一区二区三区| 俄罗斯特黄特色一大片| 久久国产亚洲av麻豆专区| 亚洲男人天堂网一区| www.熟女人妻精品国产| 男女下面进入的视频免费午夜 | 久久久久久大精品| 久久影院123| 午夜免费成人在线视频| 69精品国产乱码久久久| 一区二区三区激情视频| 免费人成视频x8x8入口观看| 在线国产一区二区在线| 精品一品国产午夜福利视频| 热99re8久久精品国产| 亚洲三区欧美一区| 久久久久久久久久久久大奶| av电影中文网址| 久久久国产精品麻豆| 亚洲 欧美一区二区三区| 亚洲第一青青草原| 免费av中文字幕在线| 午夜福利欧美成人| 97超级碰碰碰精品色视频在线观看| 999久久久精品免费观看国产| 18美女黄网站色大片免费观看| 欧美乱色亚洲激情| 嫩草影院精品99| 黄片播放在线免费| 在线观看免费视频网站a站| 欧美日韩亚洲国产一区二区在线观看| 精品人妻1区二区| 操出白浆在线播放| 亚洲欧美精品综合一区二区三区| 免费女性裸体啪啪无遮挡网站| 免费一级毛片在线播放高清视频 | 身体一侧抽搐| 国产精品偷伦视频观看了| 亚洲精品成人av观看孕妇| 这个男人来自地球电影免费观看| 久久国产精品影院| 午夜福利在线观看吧| 曰老女人黄片| 最近最新中文字幕大全免费视频| 成人永久免费在线观看视频| 国产一区二区三区在线臀色熟女 | 国产在线观看jvid| 夜夜夜夜夜久久久久| 成人三级黄色视频| 国产亚洲欧美精品永久| 一a级毛片在线观看| 免费av中文字幕在线| 男女下面进入的视频免费午夜 | 天堂影院成人在线观看| 日韩精品青青久久久久久| bbb黄色大片| 丝袜在线中文字幕| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼 | 精品福利观看| 老鸭窝网址在线观看| 曰老女人黄片| 国产成人精品久久二区二区91| 欧美精品啪啪一区二区三区| 18禁美女被吸乳视频| 国产亚洲精品一区二区www| 天天添夜夜摸| 91九色精品人成在线观看| 90打野战视频偷拍视频| 丝袜美足系列| 久久中文字幕人妻熟女| 成人永久免费在线观看视频| 亚洲九九香蕉| 一本综合久久免费| 午夜福利在线免费观看网站| 亚洲中文日韩欧美视频| 91成年电影在线观看| 美国免费a级毛片| 国内毛片毛片毛片毛片毛片| 色婷婷久久久亚洲欧美| 久久久久久久久久久久大奶| 一二三四在线观看免费中文在| 亚洲人成电影观看| 欧美日韩亚洲高清精品| 热99re8久久精品国产| av天堂在线播放| 久久久国产成人精品二区 | 中文欧美无线码| 免费在线观看影片大全网站| 国产成人一区二区三区免费视频网站| 婷婷精品国产亚洲av在线| 在线av久久热| 在线观看免费视频网站a站| 欧美大码av|