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

    基于MERIS數(shù)據(jù)的滇池葉綠素濃度時空變化(2003-2009年)及趨勢*

    2011-12-11 08:01:18朱晶晶王勝強(qiáng)松下文經(jīng)
    湖泊科學(xué) 2011年4期
    關(guān)鍵詞:滇池氣溶膠反射率

    朱晶晶,陳 晉**,王勝強(qiáng),楊 偉,,松下文經(jīng)

    (1:地表過程與資源生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,北京師范大學(xué),北京100875)

    (2:日本筑波大學(xué)地球科學(xué)系,筑波305-8572,日本)

    基于MERIS數(shù)據(jù)的滇池葉綠素濃度時空變化(2003-2009年)及趨勢*

    朱晶晶1,陳 晉1**,王勝強(qiáng)1,楊 偉1,2,松下文經(jīng)2

    (1:地表過程與資源生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,北京師范大學(xué),北京100875)

    (2:日本筑波大學(xué)地球科學(xué)系,筑波305-8572,日本)

    基于多時相MERIS數(shù)據(jù),本文對滇池葉綠素a濃度的時空變化趨勢進(jìn)行了研究.以野外實(shí)測數(shù)據(jù)為基礎(chǔ),對應(yīng)用較好的三種葉綠素a濃度反演模型進(jìn)行了驗(yàn)證比較,通過精度評價和誤差分析選擇最優(yōu)的三波段模型;將其應(yīng)用到經(jīng)過幾何糾正和大氣糾正等預(yù)處理后的MERIS數(shù)據(jù)系列,得到2003-2009年時間序列下的57幅滇池葉綠素a濃度分布圖.分析結(jié)果表明,在2003-2009年間,滇池葉綠素a濃度總體呈周期性波動緩慢上升趨勢,每年2-3月濃度最低,4-11月較快上升達(dá)到峰值,12-1月開始回落;空間分布上,葉綠素a濃度高值區(qū)多出現(xiàn)在滇池北部、岸邊帶等靠近城鎮(zhèn)的人口聚集區(qū)域,且濃度越高的季節(jié),空間分布差異越顯著.

    葉綠素a;時空變化;滇池;MERIS數(shù)據(jù);三波段模型

    近年來隨著各種污染物質(zhì)的排放,內(nèi)陸湖泊水體富營養(yǎng)化現(xiàn)象日益嚴(yán)重,直接影響到當(dāng)?shù)厝嗣竦纳詈蜕a(chǎn)[1].富營養(yǎng)化水體的典型表現(xiàn)是藻類物質(zhì)的大量繁殖,而葉綠素a作為藻類的重要組成部分,其濃度常被作為表征水體富營養(yǎng)化程度的重要指標(biāo).傳統(tǒng)水質(zhì)監(jiān)測通過采集水樣、過濾、萃取以及分光光度計(jì)分析來測定葉綠素a濃度,雖然技術(shù)成熟、相對精度較高,但人力、物力消耗大,難以實(shí)現(xiàn)大區(qū)域的實(shí)時監(jiān)測[2].遙感技術(shù)作為一種區(qū)域性的調(diào)查手段,具有監(jiān)測范圍廣、頻率高、成本低的優(yōu)勢,為實(shí)時的大范圍水質(zhì)動態(tài)監(jiān)測提供了有力工具[3-4].

    遙感手段監(jiān)測水質(zhì)一般通過建立水體反射光譜與葉綠素a濃度之間的關(guān)系模型,從而反演葉綠素a濃度.就Ⅱ類水體而言,水體反射光譜除與葉綠素a濃度相關(guān)外,還受到水中其他組分如黃質(zhì)、非浮游藻類懸浮物等的影響,使得葉綠素a濃度的反演精度相比Ⅰ類水體(海洋等)具有更大的不確定性[5].許多學(xué)者對此作了大量研究,建立了各種模型來提高反演精度,如Neville和Gower提出的葉綠素?zé)晒飧叨确P停?]、形式簡單的兩波段模型[7]、Gitelson等提出的三波段模型[8]等,但這些算法往往具有一定的區(qū)域依存性,其參數(shù)需要根據(jù)研究水域的區(qū)域特征進(jìn)行修正.遙感手段監(jiān)測水質(zhì)的另一個問題是現(xiàn)有傳感器的時空分辨率難以滿足實(shí)際工作的要求,如:LansatTM/ETM等數(shù)據(jù)空間分辨率較高(30m),但時間分辨率偏低(16d);而MODIS雖然可以每日獲得數(shù)據(jù),但其空間分辨率(250-1000m)又不能精細(xì)地刻畫水質(zhì)的空間差異.為此,歐共體專門設(shè)計(jì)了面向水色遙感的傳感器MERIS(Medium Resolution Imaging Spectrometer,搭載于環(huán)境遙感衛(wèi)星ENVISAT上),使其同時具有較高的重訪周期(3d)和空間分辨率(300m),極大地提高了Ⅱ類水體水質(zhì)時空動態(tài)監(jiān)測的能力[9].但遺憾的是,MERIS數(shù)據(jù)在我國內(nèi)陸湖泊的應(yīng)用仍十分缺乏.

    滇池作為云南省昆明市的重要水資源之一,近年來富營養(yǎng)化程度日趨嚴(yán)重.受數(shù)據(jù)和經(jīng)費(fèi)條件的限制,其葉綠素a濃度的時空動態(tài)規(guī)律和發(fā)展趨勢的掌握仍存在較大不確定性,直接影響到各種環(huán)境治理措施和政策的實(shí)施.因此,本文試圖通過野外實(shí)測數(shù)據(jù)構(gòu)建一種適合滇池區(qū)域精度較高的葉綠素a反演模型,在此基礎(chǔ)上,利用多時相MERIS數(shù)據(jù)實(shí)現(xiàn)滇池多年多季相葉綠素a濃度分布制圖,并分析葉綠素a濃度的時空變化規(guī)律和發(fā)展趨勢.

    1 研究區(qū)域與數(shù)據(jù)

    滇池(24°50'N,102°41'E)位于云南省昆明市南端,面積約 300km2,平均水深為 4.3m,最深處為 11.3m.近年來,該湖4-11月均出現(xiàn)全湖性藍(lán)藻水華暴發(fā),富營養(yǎng)化非常嚴(yán)重[10].

    為獲取實(shí)地水體反射率光譜和葉綠素a濃度數(shù)據(jù),我們在2007年10月23日、2009年3月12日和2009年8月1-7日對滇池進(jìn)行了三次野外觀測,采樣點(diǎn)數(shù)目分別是3、3、47個,共計(jì)53個.觀測項(xiàng)目包括水體光譜(光譜儀采集)和葉綠素a濃度(水樣實(shí)驗(yàn)室內(nèi)測定).同時收集了2003-2009年間所有覆蓋滇池且基本無云的MERIS1B影像(共57幅).其中對應(yīng)于實(shí)地采樣日期有兩幅影像,分別是2007年10月24日(軌道29531)和2009年3月13日(軌道36774),但影像獲取時間均晚一天.整體而言,受夏季云覆蓋影響,5-8月影像數(shù)據(jù)較少,而秋冬季數(shù)據(jù)(11-4月)較完備,其中2007年數(shù)據(jù)最全.

    2 數(shù)據(jù)預(yù)處理

    2.1 實(shí)測數(shù)據(jù)預(yù)處理

    水體光譜采集時間為當(dāng)?shù)貢r間10:00-14:00點(diǎn),使用美國ASD FieldSpec手持式光譜儀,采樣間隔1nm,波譜范圍325-1075nm.測量參數(shù)包括波長為λ時的離水幅亮度Lu(λ)、總的下行輻照度Ed(λ)、天空光漫射幅亮度Lsky(λ).恰好處于水面以上的離水遙感反射率Rrs(λ)(sr-1),按照下式計(jì)算[11-12]:

    式中,Cal(λ)是波長為λ時的標(biāo)準(zhǔn)白板的反射比;F(θ)是太陽天頂角為θ時的天空漫射光反射比,可通過氣-水界面的菲涅爾反射系數(shù)公式計(jì)算得到[13].計(jì)算結(jié)果如圖1所示.為方便與MERIS數(shù)據(jù)對比,需將離水遙感反射率Rrs(λ)轉(zhuǎn)換為水體反射率R(λ),可根據(jù)水面BRDF特征對Rrs(λ)進(jìn)行半球空間積分得到[13].

    葉綠素a濃度在實(shí)驗(yàn)室內(nèi)用分光光度計(jì)方法測量[14].首先用100%甲醇萃取葉綠素a(置于4℃黑暗環(huán)境下24h);將萃取后的葉綠素a置于UV-1700分光光度計(jì)上測量750、663、645、630nm的吸光度,則葉綠素a濃度[Chl.a](μg/L)可通過下式計(jì)算得到:

    式中,V甲醇是萃取葉綠素a的甲醇體積(ml);ODλ是UV測得的波長為λ(nm)時的吸光度;V水樣是水樣體積(ml).滇池53個樣點(diǎn)的葉綠素a濃度統(tǒng)計(jì)結(jié)果是:平均 值 80.72μg/L,最 小 值 22.98μg/L,最 大 值318.60μg/L,標(biāo)準(zhǔn)偏差 46.34μg/L,變異系數(shù)約為57.41%.

    2.2 MERIS數(shù)據(jù)預(yù)處理

    MERIS傳感器1B產(chǎn)品包含幾何糾正后的大氣層頂幅亮度和輔助數(shù)據(jù)(經(jīng)緯度,太陽天頂角、方位角,觀測天頂角、方位角,DEM高程,風(fēng)速,臭氧濃度,相對濕度等)[9].首先,利用 MERIS數(shù)據(jù)處理軟件BEAM將大氣層頂幅亮度轉(zhuǎn)換為大氣頂層反射率,并進(jìn)行投影轉(zhuǎn)換;然后,選擇一幅影像作為基準(zhǔn),對其他影像進(jìn)行Image to Image的幾何配準(zhǔn).

    圖1 實(shí)測水體遙感反射率光譜Fig.1 Measurements of remote-sensing reflectance just above the water surface

    一般而言,水色遙感器接收到的總輻射中反映水體生物光學(xué)特征的信號僅占5%-15%,必須通過大氣糾正得到較準(zhǔn)確的水體反射率影像[15].我們采用暗水體法對滇池進(jìn)行大氣糾正,該方法由Wang和Gordon于1994年提出,主要針對Ⅰ類水體[16].它假設(shè)近紅外波段(>700nm)水體輻亮度近似為0,衛(wèi)星所接收到該波段范圍的輻射全部來自大氣效應(yīng);由此可推算出當(dāng)?shù)氐拇髿鈪?shù)和氣溶膠糾正因子,進(jìn)而對可見光波段的輻射值進(jìn)行校正.滇池作為Ⅱ類水體,該假設(shè)無法成立;但滇池附近的撫仙湖為清潔的Ⅰ類水體,其大氣參數(shù)可通過上述方法求得.因而在假設(shè)大氣狀況在一定區(qū)域內(nèi)穩(wěn)定的前提下,先按照暗水體法求得撫仙湖的氣溶膠糾正因子,并外推至滇池區(qū)域,從而實(shí)現(xiàn)對滇池的大氣糾正.暗水體法基本原理如下:

    到達(dá)衛(wèi)星傳感器的總輻射能量,即波長為λ的大氣層頂輻射值Lt(λ)可描述為:

    其中,Lr(λ)為大氣瑞利散射,La(λ)為氣溶膠散射,Lra(λ)為氣溶膠與大氣分子之間的多次散射,Lg(λ)為直射太陽光的鏡面反射,Lb(λ)為來自水體底部的反射,Lw(λ)為來自水體的輻射,t為水面到衛(wèi)星傳感器之間的大氣衰減系數(shù).由于傳感器在設(shè)計(jì)時一般都會避開鏡面反射角度,因此Lg(λ)項(xiàng)可以忽略;滇池的水深狀況可保證可見光不會到達(dá)水體底部,則Lb(λ)項(xiàng)可以忽略;若只考慮氣溶膠的單次散射效應(yīng),則可把La(λ)和 Lra(λ)合寫成 Las(λ).因而(3)式可簡化為:

    將幅亮度L轉(zhuǎn)換成反射率ρ后,上式轉(zhuǎn)變?yōu)?

    式中,瑞利散射反射率ρa(bǔ)s(λ)和大氣衰減系數(shù)t可參照文獻(xiàn)的方法計(jì)算得到[17-18];氣溶膠反射率的計(jì)算分步驟進(jìn)行,首先引進(jìn)波長 λi、λj的氣溶膠糾正因子 ε(λi,λj):

    式中,ωa是氣溶膠單次散射反照度,τa是氣溶膠光學(xué)厚度,pa是氣溶膠散射相函數(shù),i、j表示波長λ的波段序號,θ0、θ分別是太陽、衛(wèi)星天頂角.一般情況下,氣溶膠散射相位函數(shù)與波長無關(guān),氣溶膠單次散射反照度是波長的弱函數(shù),而氣溶膠光學(xué)厚度隨波長而變[19]:

    式中,α為混濁度因子;β是埃斯屈朗指數(shù).由此(6)式可簡化為:

    式中,n為待定常數(shù),表示氣溶膠散射對波長的依賴性,可通過實(shí)際數(shù)據(jù)計(jì)算得到.

    實(shí)際操作中,針對撫仙湖這一清潔水體,在假設(shè)近紅外波段水體輻亮度近似為0的前提下,選擇MERIS 12、13 波段(778.75、865nm),利用式(5)、(8)計(jì)算這兩個波段下的氣溶膠糾正因子 ε(778.75,865),進(jìn)而求得待定常數(shù)n.然后以較長的13波段作為參考,計(jì)算出波長為λ時的氣溶膠糾正因子ε(λ,865),并將ε(λ,865)外推應(yīng)用到滇池水域,依照式(6)計(jì)算出滇池的氣溶膠反射率.然后利用式(5)從衛(wèi)星信號中分離出滇池的水體反射率.

    2007年10月24日和2009年3月13日滇池實(shí)測光譜數(shù)據(jù)中有6個樣點(diǎn)的采樣時間與MERIS衛(wèi)星過境時刻基本一致(日期相隔一天).對所對應(yīng)的MERIS數(shù)據(jù)進(jìn)行大氣糾正得到水體反射率,并與實(shí)測光譜相比較(圖2).大氣糾正后的水體反射率光譜在形狀和數(shù)值上均與實(shí)測值很接近,各樣點(diǎn)的平均相對誤差(MRE)均小于30%,6個樣點(diǎn)的總體MRE為19.67%,在可以接受的精度范圍內(nèi).

    3 葉綠素a濃度反演模型建立與驗(yàn)證

    3.1 模型構(gòu)建

    葉綠素a濃度的反演模型通過實(shí)測光譜數(shù)據(jù)和葉綠素a濃度建立回歸關(guān)系實(shí)現(xiàn).首先將實(shí)測光譜采樣到MERIS波段;然后選擇三種常用且精度較高的葉綠素a濃度反演模型(葉綠素?zé)晒飧叨确P?、兩波段模型和三波段模?,分別建立其模型指數(shù)與葉綠素a濃度的回歸關(guān)系,選擇精度較高的模型作為本研究區(qū)葉綠素a濃度的反演模型.現(xiàn)地采樣點(diǎn)共計(jì)53個,隨機(jī)選取27個用于模型構(gòu)建,剩余26個進(jìn)行模型驗(yàn)證.本文所用3種模型的一般表達(dá)式如下:

    (1)熒光高度法模型:Index1=FLH(λ2)λ1-λ3,其中 λ2是熒光峰位置,λ1、λ3是其兩邊肩部位置,三者分別對應(yīng) MERIS 的 9、8、10 波段(708.75、681.25、753.75nm).

    (2)兩波段模型:Index2=R(λ1)/R(λ2),其中 λ1、λ2對應(yīng) MERIS 的9、7 波段(708.75、665nm).

    (3)三波段模型:Index3=[R-1(λ1)-R-1(λ2)]×R(λ3),其中 λ1、λ2、λ3對應(yīng) MERIS 的 7、9、10 波段(665、708.75、753.75nm).

    3種模型的回歸關(guān)系式、決定系數(shù)R2和均方根誤差(RMSE)見表1.三波段模型精度相對于其他兩個模型更高,因此將其選作滇池葉綠素a濃度的反演模型,以[Chl.a]表示葉綠素a濃度(μg/L),R(λ)表示波長為λ時的水體反射率,則該模型可表示為:

    該模型的 R2為0.86,RMSE 為10.73μg/L,n 為53.

    表1 熒光高度法、兩波段和三波段模型在滇池區(qū)域的比較*Tab.1 Comparison of FLH,Two-band and Three-band model in Lake Dianchi

    3.2 模型驗(yàn)證

    利用滇池26個驗(yàn)證樣點(diǎn)的實(shí)測數(shù)據(jù)對滇池三波段模型(式9)進(jìn)行檢驗(yàn),結(jié)果表明總體RMSE為13.86μg/L(圖3a),但在葉綠素a濃度高值區(qū),模型有低估的趨勢.此外,我們也利用MERIS遙感數(shù)據(jù)對模型進(jìn)行了驗(yàn)證,以保證模型在遙感數(shù)據(jù)上應(yīng)用的有效性.圖4對比了2007年10月24日和2009年3月13日6個樣點(diǎn)的葉綠素a濃度實(shí)測數(shù)據(jù)和MERIS數(shù)據(jù)反演結(jié)果,反演結(jié)果與實(shí)測濃度基本趨勢一致,總體RMSE為14.01μg/L,略大于直接使用實(shí)測光譜數(shù)據(jù)的誤差(圖3b),但與滇池葉綠素a濃度范圍20-320μg/L相比,模型精度在可以接受的范圍內(nèi).

    圖2 實(shí)測光譜與大氣糾正前后的光譜比較Fig.2 Comparison between measured and atmospheric-corrected water reflectance

    4 結(jié)果分析

    將所建立的葉綠素a反演三波段模型(公式9)應(yīng)用到滇池MERIS系列數(shù)據(jù)上,計(jì)算得到2003-2009年57幅滇池葉綠素a濃度分布圖.具體流程為:1)對MERIS圖像進(jìn)行幾何糾正、大氣糾正等預(yù)處理;2)根據(jù)每一個像元的水體反射率計(jì)算三波段指數(shù),代入公式(9)計(jì)算得到葉綠素a濃度.

    圖3 實(shí)測葉綠素a濃度與三波段模型預(yù)測值(a)和基于MERIS數(shù)據(jù)的模型反演值(b)比較Fig.3 Relationship between measured and estimated Chlorophyll-a concentration by three-band model with in-situ spectra(a)and by three-band model with MERIS data(b)

    圖4 2003-2009年間滇池葉綠素a濃度變化Fig.4 Variation of chlorophyll-a concentration from 2003 to 2009 in Lake Dianchi

    4.1 滇池葉綠素a濃度的時間變化趨勢

    2003-2009年滇池全湖葉綠素a濃度平均值變化范圍約為 10-130μg/L,其中,2003、2007、2009年為相對高值年,2004、2005年為相對低值年(圖4).不同年份同一月份葉綠素a濃度的變化情況見圖5,受數(shù)據(jù)時間覆蓋的限制,這里只分析11-4月.總體而言,濃度越高,年際差異越顯著,較高濃度月份(11、12、1月)呈現(xiàn)周期性波動中上升趨勢.已有研究表明,藻類生長狀況是湖泊營養(yǎng)條件、溫度、光照、風(fēng)及風(fēng)浪等因素綜合作用的結(jié)果[20-24].下面我們主要從營養(yǎng)條件和氣象因素角度分析滇池葉綠素a年際變化的原因.

    昆明市環(huán)境監(jiān)測中心的監(jiān)測數(shù)據(jù)表明[20],近二十年來,滇池的總氮、銨氮和總磷年均濃度呈波動式上升趨勢,而氮、磷是藻類生長的重要營養(yǎng)物質(zhì),這是導(dǎo)致葉綠素a濃度呈現(xiàn)上升趨勢的根本原因.而年際波動則主要受水溫及湖水量變化的影響,一般而言,水溫較高有利于藻類生長繁殖,而降水導(dǎo)致的湖水量增加則可以稀釋營養(yǎng)物質(zhì)及藻類濃度.中國氣象局提供的數(shù)據(jù)表明:滇池流域年平均氣溫分別約16.50℃(2003年)、15.67℃(2004年)、16.68℃(2005 年)、16.38℃(2006 年)、15.64℃(2007 年)、15.42℃(2008 年)和 16.57℃(2009 年),年降水量分別為833mm(2003 年)、1094mm(2004 年)、976mm(2005 年)、993mm(2006 年)、933mm(2007 年)、978mm(2008)和566mm(2009年).可以看出,2003-2009年間氣溫與降水量的匹配模式存在一定的差異,高溫少雨年(如2003、2009年)一般對應(yīng)于葉綠素a濃度高值年,而低溫多雨年(如2004年)一般對應(yīng)于葉綠素a濃度低值年.但這里2007年異常高值和2005年相對低值卻不能很好地從氣象要素的變化角度解釋.2005年的低值估計(jì)與該年度昆明市實(shí)施了入湖河流截污整治、湖濱帶生態(tài)恢復(fù)等多項(xiàng)工程有關(guān),而2007年異常高值可能與該年度污染事故等原因?qū)е氯牒⒘琢棵驮鲇嘘P(guān).由于缺乏具體監(jiān)測數(shù)據(jù),上述推論仍需進(jìn)一步考證.

    2007年作為高值年且數(shù)據(jù)較完備,故進(jìn)一步分析了其年內(nèi)變化規(guī)律.滇池葉綠素a濃度在一年內(nèi)呈現(xiàn)“降-升-降”的季節(jié)變化模式(圖6).其中,2、3月葉綠素a濃度最低,全湖平均值約為15μg/L;4-6月開始緩慢上升,達(dá)到40μg/L左右;9-11月達(dá)到較高值,約為103μg/L;12月開始下降,至2月降至最低濃度.這與觀測到的4-11月份間滇池處于藻類暴發(fā)期的事實(shí)基本一致[10].

    圖5 2003-2009年間滇池葉綠素a濃度不同月份變化Fig.5 Variation of chlorophyll-a concentration by the month from 2003 to 2009 in Lake Dianchi

    圖6 2007年滇池葉綠素a濃度變化Fig.6 Variation of chlorophyll-a concentration in 2007,Lake Dianchi

    上述季節(jié)變化模式與滇池氣候特點(diǎn)密切相關(guān).滇池地處亞熱帶高原,屬冬暖夏涼性季風(fēng)氣候,氣溫年較差小.低溫期出現(xiàn)在12-2月,月平均氣溫小于10℃;高溫期出現(xiàn)在5-9月,月平均氣溫約為20℃.但降水量年內(nèi)分布不均,6-8月雨季集中了全年60%以上的降雨,而冬、春季降水較少.這種溫度與降水量的匹配模式?jīng)Q定了葉綠素a濃度的年內(nèi)變化規(guī)律.高溫少雨月份(9-11月)不僅溫度適宜藻類生長繁殖,而且湖水量較少也有助于提高葉綠素a濃度;而對于低溫少雨月份(12-3月),低溫則作為主要限制因子抑制了藻類生長繁殖.夏季月份(6-8月)雖然溫度適宜藻類生長繁殖,但降水導(dǎo)致的湖水量增加卻稀釋了葉綠素a濃度.總體而言,滇池葉綠素a濃度季節(jié)變化受溫度的影響更為顯著.此外,各月葉綠素a濃度也與前期月份藻類生長累積有一定關(guān)系,導(dǎo)致葉綠素a濃度變化具有漸變和時間后延的特點(diǎn).

    4.2 滇池葉綠素a濃度的空間變化趨勢

    滇池2007年除5、7、8月外所有月份葉綠素a濃度的空間分布情況說明其具有明顯的空間分異和季節(jié)變化特點(diǎn)(圖7).2-3月,滇池處于全年水溫最低的時期,藻類生長不活躍,全湖除最北端的草海外大部分區(qū)域的葉綠素a濃度值均低于50μg/L(深藍(lán)色、藍(lán)色表示),空間差異不顯著.4月開始,伴隨春季氣溫回升,藻類的生長開始活躍,全湖葉綠素a濃度逐漸升高,高值區(qū)由草海沿北部近岸水域逐漸向中部擴(kuò)大.6月全湖葉綠素a濃度有所降低,這可能與降水帶來的湖水量增加稀釋了葉綠素a濃度有關(guān).9-12月全湖整體濃度可超過100μg/L,且呈現(xiàn)由北向南、由西向東遞減的趨勢,中部、南部的個別區(qū)域也出現(xiàn)點(diǎn)狀高濃度區(qū).我們進(jìn)一步分析其他年份葉綠素a濃度的空間分布特點(diǎn),所得到的規(guī)律與2007年基本相似,不同年份僅存在濃度高低差異.

    圖7 2007年滇池葉綠素a濃度空間分布Fig.7 Spatial distribution of chlorophyll-a concentration in 2007,Lake Dianchi

    上述空間分布特點(diǎn)與滇池周邊污染源分布密切相關(guān),北部的草海由于靠近昆明市,接受了大量的生活和工農(nóng)業(yè)污水,氮、磷等營養(yǎng)物質(zhì)含量和水溫均偏高,且草海與南部的滇池外海區(qū)域連接的河口較小,水量交流緩慢,污染的凈化過程遲緩[10,24],導(dǎo)致其一直是全湖氮、磷濃度最高的區(qū)域.而北部近岸帶(如海埂等)由于靠近昆明市,附近聚集著大量的旅游度假中心及村鎮(zhèn),是污染物排放的重要地點(diǎn),氮、磷等營養(yǎng)物質(zhì)含量也較高,一旦溫度回升就會導(dǎo)致藻類大量繁殖[10,24].西岸??诤湍隙说睦リ査蛴捎诟浇S的排放濃度較高,也出現(xiàn)點(diǎn)狀分布的高值區(qū).其他濃度較高的地點(diǎn),則是由于季風(fēng)作用導(dǎo)致藻類聚集形成的.理論上來說,一定強(qiáng)度的風(fēng)浪可抑制藻類的堆積,但微風(fēng)則相反.滇池流域雖處于多風(fēng)地帶,但大風(fēng)期多集中于冬春季,主導(dǎo)風(fēng)向?yàn)槲髂巷L(fēng),9-12月風(fēng)速較小,有利于藻類聚集于水體表層,并向西南部漂移,因而在湖濱地帶以及風(fēng)力較小的灣內(nèi)易成為葉綠素a濃度的較高值點(diǎn)[24].

    5 結(jié)論

    本文基于野外實(shí)測數(shù)據(jù)構(gòu)建了滇池葉綠素a濃度反演的三波段模型,并用實(shí)測數(shù)據(jù)和MERIS衛(wèi)星數(shù)據(jù)進(jìn)行了模型驗(yàn)證;利用該模型和預(yù)處理后的MERIS水體反射率數(shù)據(jù),對2003-2009年滇池葉綠素a濃度分布進(jìn)行了制圖分析,得到如下結(jié)論:

    (1)2003-2009年滇池全湖葉綠素a濃度呈周期性波動上升趨勢,具有明顯的年際變化和季節(jié)變化特點(diǎn).每年2、3月濃度最低,全湖平均值約為15μg/L,且不同年份變化不大;4-6月濃度開始上升,達(dá)到40μg/L左右;而9-11月濃度則達(dá)到最高值(超過60μg/L),最高值月份年際差異明顯;12月濃度開始下降,但依然處于較高水平,直至2月才降至較低濃度.

    (2)滇池葉綠素a濃度空間分布亦具有明顯季節(jié)變化特征.2-3月,全湖大部分區(qū)域的葉綠素a濃度值較低,只有最北端的草海水域濃度較高.4-6月開始升高,且高值區(qū)由草海沿北部近岸水域逐漸向中部擴(kuò)大.9-12月,全湖整體濃度較高,呈現(xiàn)由北向南、由東向西遞減的趨勢,中部、南部的個別區(qū)域也出現(xiàn)點(diǎn)狀高濃度區(qū).

    (3)滇池葉綠素a濃度的時空變化是湖泊營養(yǎng)條件、溫度、降水、風(fēng)及藻類種類等因素共同作用的結(jié)果.年際和年內(nèi)波動主要受溫度、降水等氣象因素控制;空間分布差異則主要受氮、磷等營養(yǎng)物質(zhì)排放源的分布以及風(fēng)速、風(fēng)向的影響.

    本文尚有以下不足:首先,MERIS數(shù)據(jù)時間覆蓋范圍不夠,多集中于秋冬春季節(jié),導(dǎo)致夏季規(guī)律分析受到限制;其次,三波段模型估測葉綠素a濃度時存在低估情形.因此,提高三波段模型的估算精度以及考慮用其他衛(wèi)星數(shù)據(jù)來增強(qiáng)數(shù)據(jù)的時間覆蓋率,將有助于全面、準(zhǔn)確地掌握滇池葉綠素a濃度的時空動態(tài)規(guī)律和發(fā)展趨勢.

    [1]Ayres W,Busia A,Dinar A et al.Integrated lake and reservoir management:World band approach and experience.Washington DC:World Bank,1996.

    [2]Lorenzen CJ.Determination of chlorophyll and Pheo-Pigments:Spectrophotometric Equations.Limnology and Oceanography,1967,12(2):343-346.

    [3]Dall'Olmo G,Gitelson AA,Rundquist DC.Towards a unified approach for remote estimation of chlrophyll-a in both terrstrial vegetatioin and turbid productive waters. GeophysicalReseach Letters, 2003, 30(doi:10.1029/2003GL018065).

    [4]Kutser T,Metsamaa L,Strombeck N et al.Monitoring cyanobacterial blooms by satellite remote sensing.Estuarine Coastal and Shelf Science,2003,67:303-312.

    [5]Gitelson AA,Schalles JF,Rundquist DC et al.Comparative reflectance properties of algal cultures with manipulated densities.Journal of Applied Phycology,1999,11(4):345-354.

    [6]Neville RA,Gower JFR.Passive remote sensing of phytoplankton via chlorophyll fluorescence.Journal of Geophysical Research,1977,82(24):3487-3493.

    [7]Gitelson AA,Keydan GP.Remote sensing of inland waters quality Measurements in the visible spectrum.Acta Hydrophysi-ca,1990,34:5-27.

    [8]Gitelson AA,Dall'Olmo G,Moses W et al.A simple semi-analytical model for remote estimation of chlorophyll-a in turbid waters:Validation.Remote Sensing of Environment,2008,112(9):3582-3593.

    [9]European Space Agency.MERIS Product Handbook,2006:10-28.

    [10]Jin X,Wang L,He L.Experience and lessons learned brief for Lake Dianchi.2003:159-178.

    [11]Lee Z,Carder KL,Mobley CD et al.Hyperspectral remote sensing for shallow waters.2.Deriving bottom depths and water properties by optimization.Applied Optics,1999,38(18):3831-3843.

    [12]Mobley CD.Estimation of the remote-sensing reflectance from above-surface measurements.Applied Optics,1999,38(36):7442-7455.

    [13]Dekker AG.Analytical algorithms for lake water TSM estimation for retrospective analyses of TM and SPOT sensor data.International Journal of Remote Sensing,2002,23(1):15-35.

    [14]SCOR-UNESCO.Determination of photosynthetic pigment in seawater.Monographs on oceanographic methodology.France:In Paris,1966.

    [15]Huot JP,Tait H,Rast M et al.The optical imaging instruments and their applications:AATSR and MERIS.ESA Bulletin,2001,106:56-66.

    [16]Gordon HR,Wang M.Retrieval of water-leaving radiance and aerosol optical thickness over the oceanswith Sea WiFS:a preliminary algorithm.Applied Optics,1994,33(3):443-452.

    [17]Selby JEA,Kneizys FX,Chetwynd JHJ et al.Atmospheric transmittance/radiance:computer code LOWTRAN 4.AFGLTR-78-0053,1978.

    [18]Gordon HR,Brown JW,Evans RH.Exact Rayleigh scattering calculations for use with the Nimbus-7 Coastal Zone Color Scanner.Applied Optics,1988,27(5):862-871.

    [19]?ngstr?m AK.On the atmospheric transmission of Sun radiation and on the dust on the air.Geografiska Annaler,1929,11:156-166.

    [20]王紅梅,陳 燕.滇池近20a富營養(yǎng)化變化趨勢及原因分析.環(huán)境科學(xué)導(dǎo)刊,2009,28(3):57-60.

    [21]Chen Y,Qin B,Teubner K et al.Long-term dynamics of phytoplankton assemblages:Microcystis-domination in Lake Taihu,a large shallow lake in China.Journal of Plankton Research,2003,25(1):445-453.

    [22]李 原,張 梅,王若南.滇池的水華藍(lán)藻的時空變化.云南大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,27(3):272-276.

    [23]萬 能,宋立榮,王若南等.滇池藻類生物量時空分布及其影響因子.水生生物學(xué)報(bào),2008,32(2):184-188.

    [24]Xing K,Guo H,Sun Y et al.Assessment of the spatial-temporal eutrophic chatacter in the Lake Dianchi.Journal of Geographical Sciences,2005,15(1):37-43.

    Spatial-temporal variation of chlorophyll-a concentration in Lake Dianchi from 2003 to 2009 and trend analysis based on MERIS data

    ZHU Jingjing1,CHEN Jin1,WANG Shengqiang1,YANG Wei1,2& MATSUSHITA Bunkei2
    (1:State Key Laboratory of Earth Surface Processes and Resource Ecology,Beijing Normal University,Beijing 100875,P.R.China)
    (2:Department of Geography,Tsukuba University,Tsukuba 305-8572,Japan)

    Spatial-temporal variation trends of Chlorophyll-a concentration in Lake Dianchi were analyzed based on MERIS data.First,three popular retrieval algorithms of Chlorophyll-a concentration were calibrated by field measurements.It was found that three-band model was the best one with higher estimation accuracy.Then,the optimized three-band model was applied to 57 MERIS images from 2003 to 2009 which were pre-processed through geometric correction,atmospheric correction,etc.At last,a series of maps about Chlorophyll-a concentration distribution of Lake Dianchi were produced.The result showed that Chlorophyll-a concentration was rising slowly with a periodic fluctuation in Lake Dianchi.During a year,Chlorophyll-a concentration showed a decreasing-increasing-decreasing pattern,where the minimum appeared in February and March,and the maximum occurred within September and November.Regarding to spatial distribution,Chlorophyll-a concentration was higher at the edge and the north of Lake Dianchi,which are close to the Kunming city and small towns.In addition,the higher Chlorophyll-a concentration was,the more obvious the spatial variation was.

    Chlorophyll-a;spatial-temporal variation;Lake Dianchi;MERIS data;three-band model

    * 國家自然科學(xué)基金項(xiàng)目(40871162)資助.2010-09-17收稿;2010-12-06收修改稿.朱晶晶,女,1986年生,碩士研究生;E-mail:jingjingzhu.3@gmail.com.

    ** 通訊作者;E-mail:chenjin@ires.cn.

    猜你喜歡
    滇池氣溶膠反射率
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    具有顏色恒常性的光譜反射率重建
    在滇池東岸看西山
    文苑(2020年8期)2020-09-09 09:30:56
    氣溶膠傳播之謎
    氣溶膠中210Po測定的不確定度評定
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    四川盆地秋季氣溶膠與云的相關(guān)分析
    滇池濕地,美回來了
    大氣氣溶膠成核監(jiān)測
    成年人午夜在线观看视频| 少妇人妻久久综合中文| 亚洲成人精品中文字幕电影| 欧美精品一区二区大全| 欧美区成人在线视频| 国产精品一区二区在线观看99| a级一级毛片免费在线观看| 成人综合一区亚洲| 97在线人人人人妻| 香蕉精品网在线| 三级国产精品片| 丰满少妇做爰视频| 精品午夜福利在线看| 国产一区二区在线观看日韩| 自拍偷自拍亚洲精品老妇| 深夜a级毛片| 日本黄大片高清| 亚洲国产精品999| 欧美97在线视频| 国产色爽女视频免费观看| 久久精品国产a三级三级三级| 性色av一级| 国产黄色视频一区二区在线观看| 国产精品嫩草影院av在线观看| 美女cb高潮喷水在线观看| 亚洲最大成人手机在线| 性色avwww在线观看| 欧美日韩视频高清一区二区三区二| 欧美另类一区| 午夜激情久久久久久久| 久久久国产一区二区| 各种免费的搞黄视频| 亚洲自偷自拍三级| 久久久久久久久大av| 婷婷色综合大香蕉| 国产一区二区三区av在线| av播播在线观看一区| 国产精品无大码| 亚洲熟女精品中文字幕| 欧美日韩视频精品一区| 日本与韩国留学比较| av播播在线观看一区| 日日啪夜夜撸| 国产淫片久久久久久久久| 亚洲婷婷狠狠爱综合网| 美女内射精品一级片tv| videossex国产| 一区二区三区免费毛片| 国模一区二区三区四区视频| 18禁裸乳无遮挡免费网站照片| 午夜福利视频精品| 久久韩国三级中文字幕| 在线观看免费高清a一片| 91精品伊人久久大香线蕉| 在线观看国产h片| 少妇裸体淫交视频免费看高清| 超碰av人人做人人爽久久| 日本爱情动作片www.在线观看| 亚洲,一卡二卡三卡| 欧美成人a在线观看| 大片免费播放器 马上看| 能在线免费看毛片的网站| 国模一区二区三区四区视频| 国内精品宾馆在线| 亚洲成人久久爱视频| 亚洲欧美日韩另类电影网站 | 黄色日韩在线| 亚洲经典国产精华液单| 亚洲成人av在线免费| 夫妻性生交免费视频一级片| 久久久久久久亚洲中文字幕| 狠狠精品人妻久久久久久综合| 丝袜脚勾引网站| 一级黄片播放器| 久久久久久久大尺度免费视频| 亚洲久久久久久中文字幕| 亚洲精品一区蜜桃| 观看免费一级毛片| 黄色欧美视频在线观看| 亚洲精品一二三| 国产精品.久久久| av.在线天堂| 男的添女的下面高潮视频| 麻豆久久精品国产亚洲av| 最近2019中文字幕mv第一页| 亚洲精品国产av蜜桃| 亚洲图色成人| 成年女人在线观看亚洲视频 | 日韩视频在线欧美| 尤物成人国产欧美一区二区三区| 涩涩av久久男人的天堂| 日韩一区二区视频免费看| 亚洲性久久影院| 国产白丝娇喘喷水9色精品| 成年女人看的毛片在线观看| 久久久久久九九精品二区国产| 蜜臀久久99精品久久宅男| 日本色播在线视频| 精品国产三级普通话版| 成人国产av品久久久| 婷婷色麻豆天堂久久| 日韩欧美精品v在线| 九九久久精品国产亚洲av麻豆| 国产精品爽爽va在线观看网站| 中国三级夫妇交换| 日韩av免费高清视频| 国产高清有码在线观看视频| 最后的刺客免费高清国语| 色婷婷久久久亚洲欧美| 久久精品久久久久久久性| 国产久久久一区二区三区| 午夜免费观看性视频| 免费看光身美女| 插逼视频在线观看| 欧美亚洲 丝袜 人妻 在线| 看非洲黑人一级黄片| 一级毛片 在线播放| 男女国产视频网站| av在线播放精品| 成人无遮挡网站| 亚洲精品成人av观看孕妇| 晚上一个人看的免费电影| videos熟女内射| 国产精品人妻久久久久久| 女人久久www免费人成看片| 日本-黄色视频高清免费观看| 自拍偷自拍亚洲精品老妇| 我的女老师完整版在线观看| 中国国产av一级| 一区二区三区乱码不卡18| 国产精品国产三级国产专区5o| 亚洲久久久久久中文字幕| 成人午夜精彩视频在线观看| 九草在线视频观看| 九九在线视频观看精品| 亚洲精品自拍成人| 国产在线一区二区三区精| 亚洲美女搞黄在线观看| tube8黄色片| 国产精品99久久久久久久久| 夜夜看夜夜爽夜夜摸| 国产黄色免费在线视频| 少妇丰满av| 精品午夜福利在线看| 69人妻影院| 亚洲国产色片| 国产精品久久久久久精品古装| 中文字幕人妻熟人妻熟丝袜美| 18+在线观看网站| 国产有黄有色有爽视频| 婷婷色麻豆天堂久久| 国产成人91sexporn| 激情 狠狠 欧美| 亚洲va在线va天堂va国产| 亚洲精品第二区| 亚洲人成网站高清观看| 国产毛片a区久久久久| 午夜福利网站1000一区二区三区| av国产免费在线观看| 精品久久国产蜜桃| 久久韩国三级中文字幕| 久久久久精品性色| 99热这里只有精品一区| 国产精品国产av在线观看| 亚洲成人av在线免费| 只有这里有精品99| 黑人高潮一二区| 91精品伊人久久大香线蕉| 最近2019中文字幕mv第一页| 国产精品麻豆人妻色哟哟久久| 成人高潮视频无遮挡免费网站| 男人舔奶头视频| 99久久中文字幕三级久久日本| 69人妻影院| 三级国产精品片| 精品一区二区三卡| 国产黄片美女视频| av播播在线观看一区| 国产男女内射视频| 91精品伊人久久大香线蕉| 久久久久性生活片| 黄色欧美视频在线观看| 欧美最新免费一区二区三区| 亚洲精品久久久久久婷婷小说| 黑人高潮一二区| 久久精品久久精品一区二区三区| 日韩 亚洲 欧美在线| 国产亚洲av片在线观看秒播厂| 国产白丝娇喘喷水9色精品| 午夜福利在线在线| 亚洲自偷自拍三级| 丰满少妇做爰视频| 久热久热在线精品观看| 国产精品不卡视频一区二区| 人体艺术视频欧美日本| 少妇人妻 视频| 亚洲精品国产色婷婷电影| 一本一本综合久久| 午夜精品国产一区二区电影 | 亚洲av二区三区四区| 国产真实伦视频高清在线观看| 丰满乱子伦码专区| 国产亚洲av嫩草精品影院| 嘟嘟电影网在线观看| 男女无遮挡免费网站观看| 精品国产露脸久久av麻豆| 日本午夜av视频| 亚洲国产欧美在线一区| 女人十人毛片免费观看3o分钟| 亚洲成人久久爱视频| 国产精品麻豆人妻色哟哟久久| 精品酒店卫生间| 国产午夜精品久久久久久一区二区三区| 国产av不卡久久| av国产久精品久网站免费入址| 国产v大片淫在线免费观看| 日韩成人伦理影院| 久久热精品热| 日韩三级伦理在线观看| 日本爱情动作片www.在线观看| 99久久精品国产国产毛片| 日本黄色片子视频| 亚洲av中文字字幕乱码综合| 日本免费在线观看一区| 久久精品国产自在天天线| 狂野欧美激情性xxxx在线观看| 高清欧美精品videossex| 国产91av在线免费观看| 久久久久国产网址| 亚洲欧美精品专区久久| 国产视频首页在线观看| 亚洲欧美中文字幕日韩二区| 久久久久久久久久久丰满| 伦理电影大哥的女人| 深夜a级毛片| 久久影院123| 亚洲成人久久爱视频| 久久久久久久亚洲中文字幕| 网址你懂的国产日韩在线| 亚洲精品日韩av片在线观看| 成人特级av手机在线观看| 亚洲欧美中文字幕日韩二区| 97在线人人人人妻| 亚洲成色77777| 又大又黄又爽视频免费| 在线观看美女被高潮喷水网站| 国产乱来视频区| 成人午夜精彩视频在线观看| 免费观看在线日韩| 国产精品一及| 亚洲精品影视一区二区三区av| 免费看av在线观看网站| 欧美老熟妇乱子伦牲交| 亚洲精品aⅴ在线观看| 美女国产视频在线观看| 免费观看性生交大片5| 干丝袜人妻中文字幕| 18禁在线无遮挡免费观看视频| 在线免费十八禁| 老女人水多毛片| 美女脱内裤让男人舔精品视频| 菩萨蛮人人尽说江南好唐韦庄| 啦啦啦中文免费视频观看日本| 简卡轻食公司| 国产免费又黄又爽又色| 国产男女超爽视频在线观看| 五月伊人婷婷丁香| 美女xxoo啪啪120秒动态图| 亚洲人与动物交配视频| 高清在线视频一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲精品,欧美精品| 久久久久久国产a免费观看| 毛片女人毛片| 国产亚洲av片在线观看秒播厂| 欧美日韩综合久久久久久| 草草在线视频免费看| 色播亚洲综合网| 婷婷色av中文字幕| 欧美潮喷喷水| 亚洲电影在线观看av| 最近最新中文字幕免费大全7| 久久6这里有精品| 欧美丝袜亚洲另类| 18+在线观看网站| 日韩一本色道免费dvd| 国产一区二区亚洲精品在线观看| 一区二区三区乱码不卡18| 国产日韩欧美亚洲二区| 国产精品一区二区在线观看99| 男女下面进入的视频免费午夜| 一二三四中文在线观看免费高清| 18禁动态无遮挡网站| 十八禁网站网址无遮挡 | 麻豆国产97在线/欧美| 国产综合精华液| 熟女人妻精品中文字幕| 国产伦精品一区二区三区四那| 男的添女的下面高潮视频| 久久久久久国产a免费观看| 国产淫片久久久久久久久| 男女无遮挡免费网站观看| 亚洲av日韩在线播放| 日韩av免费高清视频| 色婷婷久久久亚洲欧美| 国产欧美日韩精品一区二区| 免费看a级黄色片| 免费观看无遮挡的男女| 内地一区二区视频在线| 伦理电影大哥的女人| 菩萨蛮人人尽说江南好唐韦庄| 一本久久精品| 国产探花在线观看一区二区| 日韩av不卡免费在线播放| 久久影院123| 欧美成人精品欧美一级黄| 网址你懂的国产日韩在线| 99热网站在线观看| 成年人午夜在线观看视频| 97超碰精品成人国产| 亚洲天堂av无毛| 天天躁日日操中文字幕| 日本一本二区三区精品| 国产爱豆传媒在线观看| 嫩草影院精品99| 干丝袜人妻中文字幕| 看十八女毛片水多多多| 日日撸夜夜添| 午夜爱爱视频在线播放| 看黄色毛片网站| 亚洲欧美日韩无卡精品| 69av精品久久久久久| 精品人妻一区二区三区麻豆| av福利片在线观看| 亚洲av免费在线观看| 久久久久国产精品人妻一区二区| 国产伦理片在线播放av一区| av国产免费在线观看| 在线观看免费高清a一片| 国产大屁股一区二区在线视频| 成人毛片a级毛片在线播放| 中国国产av一级| 亚洲精品国产av蜜桃| 国产高清不卡午夜福利| 我的女老师完整版在线观看| 久久精品夜色国产| 国产伦理片在线播放av一区| 日韩伦理黄色片| 久久国内精品自在自线图片| 人人妻人人爽人人添夜夜欢视频 | 亚洲国产欧美人成| 麻豆精品久久久久久蜜桃| 国产探花极品一区二区| www.色视频.com| 黄色视频在线播放观看不卡| 特级一级黄色大片| 少妇人妻 视频| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 亚洲精品,欧美精品| 美女内射精品一级片tv| 成人欧美大片| 国产精品福利在线免费观看| 久久韩国三级中文字幕| 联通29元200g的流量卡| 亚洲熟女精品中文字幕| 王馨瑶露胸无遮挡在线观看| 99热这里只有是精品50| 欧美亚洲 丝袜 人妻 在线| 欧美zozozo另类| 午夜精品一区二区三区免费看| 亚洲av中文av极速乱| 午夜免费鲁丝| 国产色爽女视频免费观看| 日产精品乱码卡一卡2卡三| 亚洲国产精品成人综合色| 夜夜爽夜夜爽视频| 午夜日本视频在线| 女人十人毛片免费观看3o分钟| 亚洲人成网站高清观看| 亚洲av中文av极速乱| 久久午夜福利片| 波多野结衣巨乳人妻| 九九久久精品国产亚洲av麻豆| 亚洲欧美一区二区三区国产| 伊人久久精品亚洲午夜| 久久久久久久久大av| 日韩 亚洲 欧美在线| 成人高潮视频无遮挡免费网站| 亚洲国产av新网站| 欧美xxxx性猛交bbbb| 免费看光身美女| av女优亚洲男人天堂| 尾随美女入室| 亚洲美女视频黄频| 禁无遮挡网站| 尾随美女入室| 国产成年人精品一区二区| 国产精品国产三级国产专区5o| 九草在线视频观看| 亚洲国产欧美在线一区| 国产69精品久久久久777片| 国产老妇女一区| 激情 狠狠 欧美| 99久久人妻综合| 日韩av在线免费看完整版不卡| av卡一久久| 国产又色又爽无遮挡免| 亚洲av男天堂| 久久久久性生活片| 网址你懂的国产日韩在线| 一级毛片电影观看| 免费黄网站久久成人精品| 亚洲av日韩在线播放| 在线观看美女被高潮喷水网站| 一级爰片在线观看| 蜜臀久久99精品久久宅男| 97热精品久久久久久| 99久久精品一区二区三区| 久久精品国产自在天天线| 黄色视频在线播放观看不卡| 女人久久www免费人成看片| 九色成人免费人妻av| 中文字幕免费在线视频6| 纵有疾风起免费观看全集完整版| 亚洲国产欧美人成| 熟妇人妻不卡中文字幕| 在线看a的网站| 天美传媒精品一区二区| 亚洲天堂国产精品一区在线| 久久精品人妻少妇| 日日摸夜夜添夜夜爱| 亚洲av中文字字幕乱码综合| 亚洲av成人精品一二三区| 99热这里只有是精品在线观看| 99精国产麻豆久久婷婷| 欧美国产精品一级二级三级 | 久久99热这里只有精品18| 亚洲精品第二区| 肉色欧美久久久久久久蜜桃 | 熟女av电影| 亚洲国产欧美在线一区| 国产精品蜜桃在线观看| 久久久亚洲精品成人影院| 丰满少妇做爰视频| 尾随美女入室| 中文精品一卡2卡3卡4更新| 极品教师在线视频| 91午夜精品亚洲一区二区三区| 亚洲高清免费不卡视频| 人人妻人人看人人澡| 97超视频在线观看视频| 如何舔出高潮| 成人毛片a级毛片在线播放| 午夜老司机福利剧场| 尾随美女入室| 五月玫瑰六月丁香| 香蕉精品网在线| 日韩一本色道免费dvd| 成人无遮挡网站| 人人妻人人爽人人添夜夜欢视频 | 性色av一级| 国产免费一级a男人的天堂| 久久6这里有精品| 中文资源天堂在线| 久久久国产一区二区| 亚洲不卡免费看| 午夜福利在线在线| 色婷婷久久久亚洲欧美| 99久久九九国产精品国产免费| 丰满少妇做爰视频| 免费黄频网站在线观看国产| 欧美日韩精品成人综合77777| 国产美女午夜福利| 肉色欧美久久久久久久蜜桃 | 特级一级黄色大片| 一级a做视频免费观看| 国产在线一区二区三区精| 国产精品无大码| 99热6这里只有精品| 国产免费一级a男人的天堂| 日韩欧美精品免费久久| 久久久久久久久久久免费av| 精品一区在线观看国产| 观看免费一级毛片| 日本午夜av视频| 久久久久性生活片| 永久免费av网站大全| 色播亚洲综合网| 在线精品无人区一区二区三 | a级毛片免费高清观看在线播放| 97在线人人人人妻| 26uuu在线亚洲综合色| 久久精品国产亚洲网站| 两个人的视频大全免费| 禁无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 一区二区三区四区激情视频| 国产欧美日韩精品一区二区| 一本色道久久久久久精品综合| 蜜桃久久精品国产亚洲av| 国产午夜精品久久久久久一区二区三区| 黄色视频在线播放观看不卡| 天堂中文最新版在线下载 | 国产精品蜜桃在线观看| 一区二区三区精品91| 亚洲精品国产av成人精品| 欧美最新免费一区二区三区| 91在线精品国自产拍蜜月| 国产伦理片在线播放av一区| 日本猛色少妇xxxxx猛交久久| 国产日韩欧美在线精品| 亚洲av.av天堂| 直男gayav资源| 欧美日韩精品成人综合77777| 色视频www国产| 亚洲一区二区三区欧美精品 | 在线看a的网站| 一级毛片我不卡| 亚洲成人av在线免费| 亚洲高清免费不卡视频| 国产成人精品福利久久| 午夜亚洲福利在线播放| 亚洲性久久影院| 亚洲最大成人av| 99久久精品一区二区三区| 日韩av免费高清视频| 日韩视频在线欧美| 99久久中文字幕三级久久日本| 中文字幕人妻熟人妻熟丝袜美| 日韩不卡一区二区三区视频在线| 如何舔出高潮| 五月伊人婷婷丁香| 一区二区三区四区激情视频| 80岁老熟妇乱子伦牲交| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 久久久久久久精品精品| av.在线天堂| 免费看av在线观看网站| 亚州av有码| 尾随美女入室| 国产亚洲精品久久久com| 日韩人妻高清精品专区| 国产久久久一区二区三区| 成人特级av手机在线观看| .国产精品久久| 99热全是精品| 国产精品国产三级国产av玫瑰| 麻豆成人av视频| 最近手机中文字幕大全| 亚洲欧美精品专区久久| 免费少妇av软件| 久久97久久精品| av网站免费在线观看视频| 蜜桃久久精品国产亚洲av| 亚洲av福利一区| 十八禁网站网址无遮挡 | 在线观看免费高清a一片| 22中文网久久字幕| 日本三级黄在线观看| 国产高清国产精品国产三级 | 在线观看三级黄色| 免费大片黄手机在线观看| 丝袜脚勾引网站| 草草在线视频免费看| 亚洲精品一二三| 美女cb高潮喷水在线观看| 亚洲精品日韩av片在线观看| 成年女人在线观看亚洲视频 | 交换朋友夫妻互换小说| 欧美成人午夜免费资源| 在线观看国产h片| 天堂中文最新版在线下载 | 欧美亚洲 丝袜 人妻 在线| 丝袜喷水一区| 午夜免费鲁丝| 久久久久久九九精品二区国产| 99久久人妻综合| 18禁在线播放成人免费| 国产精品一区二区性色av| 蜜桃久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 肉色欧美久久久久久久蜜桃 | 国产一区亚洲一区在线观看| 国产欧美亚洲国产| 国产精品99久久99久久久不卡 | 国产在线男女| 好男人视频免费观看在线| 国产精品成人在线| 青春草视频在线免费观看| 嫩草影院入口| 国产精品福利在线免费观看| 中文字幕久久专区| 亚洲精品影视一区二区三区av| 2021天堂中文幕一二区在线观| 夜夜爽夜夜爽视频| 性插视频无遮挡在线免费观看| 自拍偷自拍亚洲精品老妇| 日韩视频在线欧美| 亚洲精品影视一区二区三区av| 国产精品成人在线| 久久精品熟女亚洲av麻豆精品| 王馨瑶露胸无遮挡在线观看| 欧美日韩在线观看h| 男女无遮挡免费网站观看| 国产免费一区二区三区四区乱码| 亚洲成人精品中文字幕电影| 日韩一区二区视频免费看| 亚洲综合精品二区| 看非洲黑人一级黄片| 亚洲精品一二三| 国产老妇伦熟女老妇高清| 精品一区二区免费观看| 十八禁网站网址无遮挡 |