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

    興凱湖春季水體懸浮顆粒物和CDOM吸收特*

    2015-05-10 01:27:02李思佳宋開山陳智文邵田田穆光熠
    湖泊科學(xué) 2015年5期
    關(guān)鍵詞:興凱湖吸收光譜點(diǎn)位

    李思佳,宋開山,陳智文,趙 瑩,邵田田,穆光熠,3,關(guān) 瑩

    (1:吉林師范大學(xué)生態(tài)環(huán)境研究所,四平 136000)

    (2:中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所,長(zhǎng)春 130102)

    (3:東北師范大學(xué)生命科學(xué)學(xué)院,長(zhǎng)春 130000)

    興凱湖春季水體懸浮顆粒物和CDOM吸收特*

    李思佳1,2,宋開山2,陳智文1,趙 瑩2,邵田田2,穆光熠2,3,關(guān) 瑩2**

    (1:吉林師范大學(xué)生態(tài)環(huán)境研究所,四平 136000)

    (2:中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所,長(zhǎng)春 130102)

    (3:東北師范大學(xué)生命科學(xué)學(xué)院,長(zhǎng)春 130000)

    為了分析興凱湖水體光學(xué)活性物質(zhì)的吸收特性、來源和空間分布以及對(duì)400~700nm范圍內(nèi)總吸收的貢獻(xiàn),于2013年5月對(duì)該水體進(jìn)行野外實(shí)驗(yàn),對(duì)水體中浮游藻類、非藻類顆粒物和有色可溶性有機(jī)物的吸收特性和水質(zhì)參數(shù)進(jìn)行測(cè)定. 結(jié)果表明:總懸浮顆粒物的吸收光譜與非藻類顆粒物相似,色素顆粒物含量較少且單一,非藻類顆粒物在總懸浮顆粒物吸收中占主導(dǎo)地位,其貢獻(xiàn)率始終在50%以上.CDOM吸收曲線的擬合函數(shù)斜率值Sg均高于其它水體.440nm處總懸浮顆粒物和非藻類顆粒物的吸收系數(shù)ap(440)、ad(440)與總懸浮顆粒物、無機(jī)懸浮顆粒物和有機(jī)懸浮顆粒物濃度相關(guān)性均較好,與葉綠素a(Chl.a)濃度的相關(guān)性較差.興凱湖與其它Ⅱ類水體的差異性表現(xiàn)在440nm處CDOM吸收系數(shù)ag(440)與Chl.a濃度、溶解性有機(jī)碳(DOC)濃度均無顯著相關(guān)性,說明DOC以無色部分為主.總體上,大興凱湖各吸收系數(shù)和水質(zhì)參數(shù)均值均低于小興凱湖,后者水質(zhì)受農(nóng)耕區(qū)退水及周圍漁業(yè)、旅游業(yè)的影響較大.

    興凱湖;吸收光譜;有色可溶性有機(jī)物;非藻類顆粒物;浮游植物

    利用水質(zhì)參數(shù)和固有光學(xué)特性與表光光學(xué)特性之間的關(guān)系構(gòu)建的水質(zhì)遙感估算分析或半分析模型,除對(duì)傳感器的穩(wěn)定性和輻射定標(biāo)有較高的要求外,還需要高精度的固有光學(xué)、表觀光學(xué)參數(shù)[1-8]作為支撐.相對(duì)于Ⅰ類水體,Ⅱ類水體生物光學(xué)模型受陸源影響較大,浮游藻類、非藻類顆粒物和有色可溶性有機(jī)物(CDOM)來源相對(duì)獨(dú)立,非藻類顆粒物和CDOM的光譜吸收均在藍(lán)光波段重疊,甚至成為水體光譜的主導(dǎo)組分[3],從而對(duì)目標(biāo)物質(zhì)光譜特征的提取造成干擾,增加了反演難度.由于Ⅱ類水體所處的地理環(huán)境特點(diǎn)及區(qū)域發(fā)展經(jīng)濟(jì)狀況不同,光活性物質(zhì)組成具有多樣性和時(shí)空差異性,模型依賴的固有光學(xué)量呈現(xiàn)區(qū)域變化特點(diǎn).對(duì)水體類型、季相和區(qū)域差異的模型系數(shù)不可統(tǒng)一參數(shù)化,需要大量不同水體研究的實(shí)驗(yàn)數(shù)據(jù)加以支撐、驗(yàn)證和調(diào)整[9].固有光學(xué)量作為水體生物光學(xué)算法構(gòu)建的重要參數(shù),其特性的研究對(duì)于提高水質(zhì)參數(shù)反演精度,實(shí)現(xiàn)區(qū)域水色遙感算法和生物光學(xué)模型參數(shù)化具有重要理論意義[10-13].因此,研究引起水體光學(xué)活性物質(zhì)和固有光學(xué)量變化機(jī)理,確定不同水體光學(xué)活性物質(zhì)的貢獻(xiàn)率尤為重要.

    吸收系數(shù)是固有光學(xué)量之一,國(guó)外已有較多針對(duì)Ⅱ類水體組分吸收特性的研究[14-15],國(guó)內(nèi)相關(guān)研究?jī)H局限于太湖[16]、巢湖[17]等富營(yíng)養(yǎng)化水體.作為中俄界湖的興凱湖,是世界上少有的面積廣闊、生態(tài)完好、物種多樣的濕地之一[18],在涵養(yǎng)水源、補(bǔ)給地下水、降解污染和維持水平衡方面具有重要作用.近幾年來農(nóng)業(yè)面源污染以及旅游業(yè)的發(fā)展導(dǎo)致我國(guó)部分湖水污染加重,已有的吸收光譜模型是否能應(yīng)用于興凱湖水體亟待研究.興凱湖獨(dú)特的地理位置及社會(huì)經(jīng)濟(jì)狀況使了解整個(gè)興凱湖湖泊生態(tài)系統(tǒng)變化受到限制,通過分析興凱湖我國(guó)部分的固有光學(xué)特性,有助于提升對(duì)整個(gè)興凱湖光學(xué)特性的認(rèn)識(shí),并為進(jìn)一步對(duì)興凱湖實(shí)施全面遙感機(jī)理監(jiān)測(cè)提供參考依據(jù).

    1 研究區(qū)概況

    興凱湖(45°01′00″~45°34′30″N,131°58′30″~133°07′30″E)地處我國(guó)黑龍江省東南部和俄羅斯遠(yuǎn)東濱海地區(qū),為造山運(yùn)動(dòng)地殼陷落形成的構(gòu)造湖,由大、小興凱湖組成.小興凱湖是興凱湖湖進(jìn)再湖退殘留而成,兩湖中間隔著一條10m左右高的沙崗,并以湖崗上的第一、二泄洪閘和新開流水道相通,湖底高程為58~69m.小興凱湖總面積為170km2,平均水深1.8m,最深為3.5m;大興凱湖是中俄兩國(guó)的界湖,略呈橢圓形,總面積為4380km2,其中北部1240km2屬我國(guó)境內(nèi),平均水深3.5m,最深處達(dá)10m.該區(qū)域?qū)儆跍貛Т箨懶约撅L(fēng)氣候,春、秋季風(fēng)大降水少,夏季溫?zé)岫嘤?,冬季漫長(zhǎng)干燥.興凱湖主要依賴地表徑流、湖面降水補(bǔ)給以及湖東北部大片沼澤濕地地下徑流補(bǔ)充,匯入河流在我國(guó)境內(nèi)有8條,流入大、小興凱湖的各有4條,東南側(cè)以松阿察河湖口與俄羅斯相連,松阿察河是興凱湖的唯一泄流通道,匯入烏蘇里江.小興凱湖屬草原湖泊且湖水較淺,湖床具有平緩斜坡,沿岸附近地勢(shì)低平多沼澤、蘆葦,湖底屬淤泥底質(zhì),為國(guó)有漁場(chǎng);大興凱湖湖水開闊,以砂粒底質(zhì)為主.

    2 材料與方法

    2.1 野外采樣

    于2012年5月18日和20日分別在小興凱湖、大興凱湖各采集10個(gè)樣點(diǎn)(1#~10#樣點(diǎn)位于小興凱湖,11#~20#樣點(diǎn)位于大興凱湖,見圖1),并用手持式GPS進(jìn)行定位導(dǎo)向.水樣采集使用有機(jī)玻璃采水器,采樣前先將采樣器和聚四氟乙烯水樣瓶潤(rùn)洗,采樣深度設(shè)置為表層0~50cm,每個(gè)點(diǎn)位采集水樣約1L,采樣后立即用黑色塑料袋罩住,避免陽光直射,放入車載保溫箱.

    圖1 興凱湖采樣點(diǎn)分布Fig.1 Distribution of sampling sites in Lake Xingkai

    2.2 參數(shù)測(cè)定

    2.2.1 葉綠素a(Chl.a)和懸浮物濃度的測(cè)定 為防止酸化引起色素溶解,采樣后立即完成水質(zhì)參數(shù)的測(cè)定.采用0.45μm Whatman GF/F膜過濾150~200ml水樣(根據(jù)水中Chl.a濃度決定過濾水樣體積),將濾膜剪碎放入15ml離心管中,加入10ml 90%的丙酮溶液,充分混合均勻后避光冷藏24~48h,對(duì)其離心后取上清液倒入1cm石英比色皿中,以90%的丙酮溶液為空白參考,置于UV-2600島津-紫外可見分光光度計(jì)測(cè)定波長(zhǎng)為630、647、664和750nm的吸光度,計(jì)算Chl.a濃度[19].總懸浮顆粒物(TSM)質(zhì)量、有機(jī)顆粒物(OSM)質(zhì)量和無機(jī)顆粒物(ISM)質(zhì)量采用烘燒、稱重的方法測(cè)定.溶解性有機(jī)碳(DOC)濃度采用TOC儀測(cè)定.不同水質(zhì)參數(shù)濃度見圖2.

    2.2.2 吸光系數(shù)的測(cè)定 采用定量濾膜技術(shù)(QFT)[20-21]測(cè)定總懸浮顆粒物的吸收系數(shù)ap(λ)、非藻類顆粒物的吸收系數(shù)ad(λ)、浮游植物色素顆粒物的吸收系數(shù)aph(λ)和CDOM的吸收系數(shù)ag(λ).用0.47μm Whatman GF/F濾膜在真空下過濾50~400ml水樣,利用島津-紫外可見分光光度計(jì)UV-2006在380~800nm范圍測(cè)定濾膜上顆粒物的吸光度,每隔1nm測(cè)定1個(gè)吸光度.為消除儀器間的差異,用各波段的吸光度減去750nm波長(zhǎng)處的吸光度以校正,并與經(jīng)蒸餾水潤(rùn)濕的同樣濕潤(rùn)程度的空白濾膜作對(duì)比.ap(λ)的計(jì)算公式為:

    (1)

    式中,V為過濾水樣的體積,S為過濾留在濾膜上的有效面積,ODs(λ)為濾膜上懸浮顆粒物的吸光度.

    采用500μl 0.1%的NaClO溶液漂白總懸浮顆粒物濾膜15min,利用分光光度計(jì)測(cè)定ad(λ),方法同上.根據(jù)線性疊加原理,aph(λ)計(jì)算公式為:

    aph(λ)=ap(λ)-ad(λ)

    (2)

    利用0.22μm微孔聚碳酸酯膜過濾的水樣測(cè)定CDOM的吸光度,方法同上.按照式(3)、(4)計(jì)算ag(λ):

    (3)

    (4)

    圖2 興凱湖所有點(diǎn)位各組分濃度Fig.2 Component concentrations of all the sampling sites in Lake Xingkai

    3 結(jié)果與分析

    3.1 TSM的吸收光譜特征

    興凱湖TSM吸收光譜特征與非藻類顆粒物吸收光譜特征相似(圖3a),結(jié)合已有的研究結(jié)果[22-23]表明,在葉綠素a濃度較低的情況下,TSM吸收表現(xiàn)出非藻類顆粒物吸收光譜特征,且不同點(diǎn)位ap(λ)變化較大,具有空間異質(zhì)性.由于水體中個(gè)別點(diǎn)位Chl.a濃度較低,光譜曲線在675nm左右出現(xiàn)不明顯的吸收峰;浮游植物色素顆粒物另一個(gè)特征吸收波段440nm未出現(xiàn)明顯的吸收峰,與同季節(jié)其他研究結(jié)果相同[22],非藻類顆粒物在藍(lán)光波段的強(qiáng)吸收使其在光譜曲線中占優(yōu)勢(shì).興凱湖1#~20#點(diǎn)位的TSM以ISM為主,OSM所占比例較小(圖2).分別對(duì)ap(440)、ap(675)與TSM、ISM、OSM、Chl.a濃度進(jìn)行相關(guān)分析(表1),ap(440)與TSM、ISM和OSM濃度具有極顯著的相關(guān)性(P<0.01),與Chl.a濃度的相關(guān)性較差,ap(675)與Chl.a、OSM濃度有極顯著的相關(guān)性(P<0.01),這與其他研究結(jié)果一致[24],進(jìn)一步說明春季興凱湖為非藻類顆粒物占主導(dǎo)的典型Ⅱ類水體,主要由于春季東北氣溫相對(duì)較低,藻類正處于生長(zhǎng)期,同時(shí)受流域內(nèi)人類活動(dòng)和流域內(nèi)農(nóng)業(yè)污染負(fù)荷增加的影響,非藻類顆粒物濃度較高[25].1#~10#、11#~20#點(diǎn)位ap(440)均值分別為21.29±10.27m-1、12.55±3.56m-1,總體上小興凱湖TSM濃度高于大興凱湖,主要受人類活動(dòng)的影響[18,26].

    表1 440、675nm組分吸收系數(shù)、各水色因子濃度之間的線性關(guān)系

    **表示在0.01水平上顯著相關(guān);*表示在0.05水平上顯著相關(guān).

    3.2 非藻類顆粒物的吸收光譜特征

    非藻類顆粒物包括水體中礦物沉積、非生命有機(jī)碎屑,如排泄物與藻類細(xì)胞降解物,以及活的非色素有機(jī)體[27].興凱湖非藻類顆粒物的光譜吸收曲線呈現(xiàn)指數(shù)衰減規(guī)律和空間異質(zhì)性(圖3b),ad(440)均值為14.67±7.19m-1.由于Ⅰ類水體中非藻類顆粒物多來自浮游植物的降解,ad(440)與Chl.a濃度具有較高的相關(guān)性,而Ⅱ類水體ad(440)通常表現(xiàn)出與懸浮物濃度較高的相關(guān)性[28].而興凱湖1#~20#點(diǎn)位的ad(440)、ad(675)與Chl.a濃度無顯著相關(guān)性,與TSM、ISM、OSM濃度相關(guān)性均較好(表1),ad(440)與顆粒物濃度呈非線性相關(guān),與其他研究結(jié)果一致[29-30],說明非藻類顆粒物主要為陸源性.進(jìn)一步進(jìn)行曲線擬合發(fā)現(xiàn),ad(440)與TSM濃度的擬合效果相對(duì)較好,其次是ISM和OSM濃度,擬合關(guān)系均普遍較差(圖3d、e、f),這與太湖2007年6月的研究結(jié)果一致[23],這可能是由于樣本數(shù)相對(duì)較少,需在后續(xù)研究中增加樣本數(shù).春季隨著水溫升高、光照時(shí)數(shù)增長(zhǎng),入湖河流徑流量增加且攜帶流域內(nèi)大量泥沙、農(nóng)耕區(qū)農(nóng)田退水和生活污水等無機(jī)顆粒物.在該區(qū)自然降水少、春季多大風(fēng)的條件下,湖心區(qū)水面開闊,底泥容易發(fā)生再懸浮,非藻類顆粒物在TSM中占優(yōu)勢(shì).同時(shí),東北地區(qū)5月氣溫相對(duì)南方仍較低,為春、夏交替季節(jié),藻類正處于生長(zhǎng)期且降解產(chǎn)物少,一方面在湖水冬季長(zhǎng)期低溫和冰封的條件下,本身換水周期長(zhǎng),少量藻類降解產(chǎn)物等有機(jī)顆粒物積累且含量穩(wěn)定;另一方面氣溫的變化導(dǎo)致湖水水流垂直對(duì)流擾動(dòng)加強(qiáng),也是非藻類顆粒物為主導(dǎo)組分的原因之一.

    總體而言,1#~10#點(diǎn)位TSM、OSM、ISM、DOC和Chl.a濃度均值均高于11#~20#點(diǎn)位,ad(440)最大值出現(xiàn)在6#點(diǎn)位,最小值出現(xiàn)在19#點(diǎn)位(圖2),大、小興凱湖ad(440)無顯著差異(圖3b).小興凱湖西北部為“857”、“856”等國(guó)營(yíng)農(nóng)場(chǎng),生活污水、農(nóng)田退水通過河道和排干流入小興凱湖,再加上小興凱湖魚類人工養(yǎng)殖及周圍旅游業(yè)的發(fā)展,都對(duì)水體產(chǎn)生了嚴(yán)重污染.另外,穆棱河流域已形成沖溝較為嚴(yán)重的水土流失區(qū)[22],也給小興凱湖帶來大量無機(jī)物.1#~10#點(diǎn)位的ad(440)最小值出現(xiàn)在1#點(diǎn)位,湖崗在大興凱湖一側(cè)形成沙灘景觀,在小興凱湖一側(cè)多沼澤,1#點(diǎn)位距各入湖河流及排干口較遠(yuǎn);11#~20#點(diǎn)位的ad(440)最大值出現(xiàn)在17#點(diǎn)位,靠近新開流,其附近有興凱湖水產(chǎn)養(yǎng)殖場(chǎng)、長(zhǎng)林子捕魚組,最小值出現(xiàn)在19#點(diǎn)位,靠近一閘口.1#~20#點(diǎn)位非藻類顆粒物在400~700nm按指數(shù)衰減規(guī)律擬合斜率Sd值的變化范圍為0.0043~0.0053nm-1,平均值為0.0048±0.0003nm-1,決定系數(shù)R2均在0.96以上.

    圖3 興凱湖水體ap(λ)、ad(λ)和aph(λ)的吸收光譜(灰色代表1#~10#樣點(diǎn),黑色代表11#~20#樣點(diǎn),下同);ad(440)與TSM、ISM、OSM濃度的相關(guān)性Fig.3 Absorption spectra of ap(λ), ad(λ),aph(λ);Correlation coefficients between ad(440) and TSM, ISM, OSM concentrations of Lake Xingkai

    3.3 浮游植物色素顆粒物的吸收光譜特征

    小興凱湖西北部為農(nóng)業(yè)耕作區(qū),尤其在近幾年的旅游業(yè)及漁業(yè)發(fā)展的帶動(dòng)下,沿岸及入湖河流攜帶大量營(yíng)養(yǎng)物質(zhì),促進(jìn)浮游植物的生長(zhǎng),水生植物系統(tǒng)較大興凱湖發(fā)達(dá).根據(jù)圖2,5#點(diǎn)位未表現(xiàn)出藻類吸收特征,營(yíng)養(yǎng)鹽濃度低.興凱湖Chl.a濃度最大值出現(xiàn)在泄洪一閘附近的1#點(diǎn)位,營(yíng)養(yǎng)鹽的聚積使浮游植物增加;最小值出現(xiàn)在11#點(diǎn)位,位于新開流與一閘之間,與小湖沒有閘口,營(yíng)養(yǎng)鹽流動(dòng)較弱.1#~10#點(diǎn)位Chl.a濃度最小值出現(xiàn)在8#點(diǎn)位,距離農(nóng)業(yè)排干口較遠(yuǎn).11#~20#點(diǎn)位Chl.a濃度最大值出現(xiàn)在20#點(diǎn)位,位于白泡子排干口處,農(nóng)田退水和生活污水的排入帶來較多營(yíng)養(yǎng)鹽,浮游植物大量生長(zhǎng).

    圖4 aph(440)、aph(675)與Chl.a濃度的相關(guān)性Fig.4 Correlation coefficients between aph(440), aph(675) and Chl.a concentration

    3.4 CDOM的吸收光譜特征

    3.4.1 CDOM的吸收 CDOM包括DOC中有色部分的60%以上,占到可溶性有機(jī)物質(zhì)的10%~90%,CDOM和DOC無色部分均來源于地表徑流的陸源輸入和浮游植物降解[42-43].由于CDOM成分復(fù)雜,難以確認(rèn)其濃度,常用355、375、400、440nm等波長(zhǎng)處的吸收系數(shù)來表征CDOM的濃度[44].由于CDOM吸收集中在紫外短波段,本文對(duì)興凱湖水體CDOM 250~700nm波段的吸收特性進(jìn)行分析,發(fā)現(xiàn)其呈指數(shù)衰減(圖5a).ag(440)介于0.37~2.17m-1之間,均值為1.19±0.54m-1,與4、6月份太湖、巢湖[17]樣點(diǎn)間差異較大.ag(440)吸收系數(shù)與Chl.a、DOC濃度均無顯著相關(guān)性(表1),與前人研究結(jié)果有所差異[16],表明興凱湖水體DOC中無色部分所占比例高.由于興凱湖春季藻類濃度低,無色DOC的主要來源不是浮游植物降解產(chǎn)物,而是河流的陸源性輸入,隨著東北地區(qū)5月氣溫、光照時(shí)數(shù)增加,紫外輻射的增強(qiáng)也促進(jìn)了CDOM的光化學(xué)降解作用,提高無色DOC的比例.興凱湖1#~10#、11#~20#點(diǎn)位的ag(440)分別為1.34±0.44、1.04±0.6m-1,最大值出現(xiàn)在13#點(diǎn)位,最小值出現(xiàn)在11#和12#點(diǎn)位.13#點(diǎn)位位于龍王廟和松阿察河出口處,松阿察河為興凱湖唯一出口,大量陸源輸入的CDOM積聚在此;而11#點(diǎn)位Chl.a濃度最小,其與大小湖的閘口距離較遠(yuǎn);12#點(diǎn)位在當(dāng)壁鎮(zhèn)附近,無農(nóng)田污水排干口,ag(440)較低.1#~10#點(diǎn)位東側(cè)點(diǎn)位ag(440)高于西側(cè),最大值出現(xiàn)在3#點(diǎn)位,最小值出現(xiàn)在9#和10#點(diǎn)位(圖6),主要是因?yàn)樾∨d凱湖東側(cè)不僅有農(nóng)田排干入湖口,而且流域內(nèi)流經(jīng)虎林市的東地河攜帶大量城市污水注入此湖.

    CDOM的分子量M(ag(250)/ag(365))(簡(jiǎn)寫為M250/365)可反映腐殖酸和富里酸在CDOM中的比例,M值越大,腐殖酸所占比例越低[45].興凱湖M250/365變化范圍為5.55~8.99(均值為7.38±0.99),均值低于太湖梅梁灣冬季水體(8.09±0.56)[16]和巢湖(8.66±0.088)[17].由于CDOM的陸源為河流攜帶有機(jī)成分的輸入,使腐殖酸的比例偏大[46].興凱湖穆棱河流域等集水區(qū)多為農(nóng)業(yè)耕作區(qū),大量輸入有機(jī)成分且春季浮游植物處于生長(zhǎng)期,色素濃度較太湖和巢湖低,表現(xiàn)為興凱湖水體腐殖酸所占比例相對(duì)較高.

    圖5 CDOM的光譜吸收曲線(a);各波段CDOM吸收光譜擬合波段的斜率Sg(b)Fig.5 Absorption spectra of CDOM(a); Shape factor Sg of CDOM absorption spectra from bonds(b)

    3.4.2 CDOM吸收曲線的擬合函數(shù)斜率值Sg與特征波長(zhǎng)吸收系數(shù)的關(guān)系Sg值反映水體中CDOM含量的高低,其值的大小與波段選取的范圍和CDOM組成有關(guān),光譜微分能夠很好地描述光譜斜率隨波段的變化情況[17].利用最小二乘法,計(jì)算CDOM吸收曲線在250~500nm波段的擬合函數(shù)斜率Sg,其均值為0.3161±0.1278nm-1.總體上高于4月太湖Sg值(0.0171~0.0032nm-1)[47],這一方面可能是由于參考波長(zhǎng)和擬合波段選擇的不同導(dǎo)致計(jì)算出的Sg值在不同水體或同一水體間存在差異[46],另一方面可能是因?yàn)椴煌wCDOM組成及來源不同.本研究選取250~290、290~320、320~400、400~500nm波段計(jì)算指數(shù)函數(shù)斜率值Sg(圖5b),不同波段內(nèi)Sg值差異較大,且Sg值隨波長(zhǎng)的增加而減小,與前人的研究結(jié)論一致[48-50].將250~290、290~320、320~400和400~500nmSg值與ag(440)進(jìn)行回歸分析,其相關(guān)系數(shù)分別為0.54、0.57、0.79和0.91,選擇相關(guān)性較好的320~400、400~500nm處進(jìn)行函數(shù)擬合(圖7),發(fā)現(xiàn)均存在較好的二次函數(shù)關(guān)系(R2=0.62,R2=0.83).針對(duì)CDOM的遙感反演,如QAA算法中光譜斜率S值應(yīng)相應(yīng)增大.

    圖6 興凱湖所有點(diǎn)位ag(440)值的空間分布Fig.6 Distribution of ag(440) of all the sampling sites in Lake Xingkai

    圖7 兩種波段下ag(440)與Sg的關(guān)系(a:320~400nm; b:400~500nm)Fig.7 Relationship between ag(440) and Sg under two bands(a: 320~400nm; b: 400~500nm)

    3.5 各組分對(duì)水體總吸收系數(shù)的貢獻(xiàn)

    水體各組分的吸收貢獻(xiàn)率可衡量太陽輻射的有效利用率和水體組分對(duì)光的衰減程度,根據(jù)線性疊加原理,水體總吸收系數(shù)是各個(gè)水色因子吸收系數(shù)的和:

    a(λ)=ad(λ)+aph(λ)+ag(λ)+aw(λ)

    (5)

    式中,aw(λ)為純水的吸收系數(shù).由于純水的吸收系數(shù)為常數(shù)[51],本文以非藻類顆粒物、色素顆粒物和CDOM的吸收系數(shù)為總吸收系數(shù),研究在PAR波段內(nèi)分析不同組分吸收系數(shù)的貢獻(xiàn)率(圖8a).總體上,興凱湖1#~20#點(diǎn)位水體各組分吸收系數(shù)貢獻(xiàn)率為非藻類顆粒物>色素顆粒物>CDOM,CDOM的貢獻(xiàn)主要集中于紫外短波波段,非藻類顆粒物的貢獻(xiàn)率始終大于50%,表明興凱湖水體非藻類顆粒物的吸收和散射是PAR波段光衰減的主要因子.約550nm處,由于葉綠素a的吸收,非藻類顆粒物的貢獻(xiàn)曲線逐漸下降,直至670nm附近到達(dá)最低谷后逐漸上升,色素顆粒物的貢獻(xiàn)曲線與非藻類顆粒物的貢獻(xiàn)曲線走勢(shì)相反,表明春季興凱湖5月水體色素含量單一并以葉綠素a為主.非藻類顆粒物和CDOM的吸收集中于短波段并呈指數(shù)遞減,其貢獻(xiàn)率曲線在約675nm后出現(xiàn)不穩(wěn)定的增幅,可能是儀器的噪聲所致.

    選取各點(diǎn)位PAR藍(lán)光(400~500nm)、綠光(500~600nm)、紅光(600~700nm)的特征波長(zhǎng)440、550、675nm(圖8b、c、d),進(jìn)行各組分貢獻(xiàn)率的研究.藍(lán)光波段小興凱湖色素貢獻(xiàn)大于CDOM,主要受農(nóng)業(yè)退水和生活污水的影響;大興凱湖靠近松阿察河的13#點(diǎn)位CDOM的貢獻(xiàn)大于色素.綠光波段小興凱湖點(diǎn)位色素貢獻(xiàn)大于CDOM,大興凱湖點(diǎn)位CDOM貢獻(xiàn)大于小湖;紅光波段處色素貢獻(xiàn)率比例大,主要是因?yàn)榉窃孱愵w粒和CDOM貢獻(xiàn)集中于短波段,點(diǎn)位1#、18#、19#、20#的色素貢獻(xiàn)率相對(duì)最大,其中18#、19#、20#點(diǎn)位分別位于一閘口、二閘口和白泡子排干附近,陸源輸入較多營(yíng)養(yǎng)鹽.

    圖8 興凱湖所有點(diǎn)位各組分平均吸收系數(shù)貢獻(xiàn)率(a);400~500nm(b)、500~600nm(c)和600~700nm(d)波段各組分平均吸收系數(shù)貢獻(xiàn)率Fig.8 The average absorption contribution(a); The average absorption contribution over the range of 400-500nm(b), 500-600nm(c) and 600-700nm(d) in all the sampling sites of Lake Xingkai

    4 結(jié)論

    [1] Zibordi G, Berthon JF, Mélin Fetal. Validation of satellite ocean color primary products at optically complex coastal sites: Northern Adriatic Sea, Northern Baltic Proper and Gulf of Finland.RemoteSensingofEnvironment, 2009, 113: 2574-2591.

    [2] Zhang Y, Liu M, Qin Betal. Modeling remote-sensing reflectance and retrieving chlorophyll-a concentration in extremely turbid case-2 waters(Lake Taihu China).IEEETransactionsonGeoscienceandRemoteSensing, 2009, 47(7): 1937-1948.

    [3] Zimba PV, Gitelson AA. Remote estimation of chlorophyll concentration in hyper-eutrophic aquatic systems: Model tuning and accuracy optimization.ContinentalShelfResearch,Aquaculture, 2006, 256: 272-286.

    [4] Gitelson AA, Schalles JF, Hladik CM. Remote chlorophyll-a retrieval in turbid, productive estuaries:Chesapeake Bay case study.RemoteSensingofEnvironment, 2007, 109(4): 464-472.

    [5] Gitelson AA, Gurlin D, Moses WJetal. A bio-optical algorithm for the remote estimation of chlorophyll-a concentration in case 2 waters.EnvironmentalResearchLetters, 2009, 4(4): 045003.

    [6] Gons HJ, Rijkeboer M, Ruddick KG. A chlorophyll-retrieval algorithm for satellite imagery(medium resolution imaging spectrometer) of inland and coastal waters.JournalofPlanktonResearch, 2002, 24(9): 947-951.

    [7] Brando VE, Dekker AG. Satellite hyperspectral remote sensing for estimating estuarine and coastal water quality.IEEETransactionsonGeoscienceandRemoteSensing, 2003, 41(6): 1378-1387.

    [8] Song K, Li L, Tedesco LPetal. Remote estimation of chlorophyll-a in turbid inland waters: Three-band model versus GA-PLS model.RemoteSensingofEnvironment, 2013, 136: 342-357.

    [9] 馬榮華,唐軍武.湖泊水色遙感參數(shù)獲取與算法分析.水科學(xué)進(jìn)展,2006, 17(5): 720-726.

    [10] Siegel DA, Michaels AF. Quantification of non-algal light attenuation in the Sargasso Sea: Implications for biogeochemistry and remote sensing.DeepSeaResearchPartII: Topical Studies in Oceanography, 1996, 43(2): 321-345.

    [11] Zhang Y, Zhang B, Wang Xetal. A study of absorption characteristics of chromophoric dissolved organic matter and particles in Lake Taihu, China.Hydrobiologia, 2007, 592(1): 105-120.

    [12] D’Sa EJ, Miller RL. Bio-optical properties in waters influenced by the Mississippi River during low flow conditions.RemoteSensingofEnvironment, 2003, 84(4): 538-549.

    [13] Lee ZP, Carder KL. Absorption spectrum of phytoplankton pigments derived from hyperspectral remote-sensing reflectance.RemoteSensingofEnvironment, 2004, 89(3): 361-368.

    [14] Lakhtaki A, Varadan VK, Varadan VV. Scattering and absorption characteristics of lossy dielectric,chiral,nonspherical objects.AppliedOptics, 1985, 24(23): 4146-4154.

    [15] Darecki M, Weeks A, Sagan Setal. Optical characteristics of two contrasting case 2 waters and their influence on remote sensing algorithms.ContinentalShelfResearch, 2003, 23: 237-250.

    [16] 張運(yùn)林,秦伯強(qiáng),楊龍?jiān)?太湖梅梁灣水體懸浮顆粒物和CDOM的吸收特性.生態(tài)學(xué)報(bào),2006, 26(12): 3970-3979.

    [17] 施 坤,李云梅,王 橋等.太湖、巢湖水體CDOM吸收特性和組成的異同.環(huán)境科學(xué),2010,31(5): 1183-1191.

    [18] 劉東興.淺析興凱湖濕地保護(hù).東北水利水電,2009,6: 50-52.

    [19] Arar EJ, Collins GB. Method 445.0: In vitro determination of chlorophyll a and pheophytin a in marine and freshwater algae by fluorescence. United States Environmental Protection Agency, Office of Research and Development, National Exposure Research Laboratory, 1997.

    [20] Bricaud A, Stramski D. Spectral absorption coefficients of living phytoplankton and non-algal biogenous matter: A comparison between the Peru upwelling area and the Sargasso Sea.LimnologyandOceanography, 1990, 35: 562-582.

    [21] Mitchell BG, Kahru M, Wieland Jetal. Determination of spectral absorption coefficients of particles, dissolved material and phytoplankton for discrete water samples.OceanOpticsProtocolsforSatelliteOceanColorSensorValidation.Revision, 2002, 3: 231-257.

    [22] 張 紅,黃家柱,李云梅等.滇池水體光學(xué)物質(zhì)的吸收特性研究.環(huán)境科學(xué),2011,32(2): 452-463.

    [23] 馬榮華,段洪濤,唐軍武等.湖泊水環(huán)境遙感.北京:科學(xué)出版社,2010.

    [24] 王遠(yuǎn)東,劉殿偉,宋開山等.吉林省查干湖水體懸浮顆粒物光譜吸收特性研究.光譜學(xué)與光譜分析,2011,31(1): 162-167.

    [25] 李 方,徐京萍,何艷芬等.長(zhǎng)春市石頭口門水庫顆粒物光譜吸收特性.湖泊科學(xué),2009,21(2): 280-287. DOI 10.18307/2009.0218.

    [26] 樸德雄,王鳳昆.興凱湖水環(huán)境狀況及其保護(hù)對(duì)策.湖泊科學(xué),2011,23(2): 196-202. DOI 10.18307/2011.0206.

    [27] 柳彩霞,郭子祺,張寶鋼等.太湖流域昆承湖春季顆粒物和有色可溶性有機(jī)物吸收特性.湖泊科學(xué),2011,23(5): 773-782. DOI 10.18307/2011.0516.

    [28] Babin M, Stramski D, Ferrari GMetal. Variations in the light absorption coefficients of phytoplankton,nonalgal particles, and dissolved organic matter in coastal waters around Europe.JournalofGeophysicalResearch, 2003, 108(C7):1-4.

    [29] Roesler CS, Perry MJ, Carder KL. Modeling in situ phytoplankton absorption from total absorption spectra in productive inland marine waters.LimnologyandOceanography, 1989, 34(8): 1510-1523.

    [30] Loisel H, Morel A. Light scattering and chlorophyll concentration in case 1 waters: Areexamination.LimnologyandOceanography, 1998, 43(5): 847-858.

    [31] 樂成峰,李云梅,查 勇等.太湖梅梁灣水體組分吸收特性季節(jié)差異分析.環(huán)境科學(xué),2008,29(9): 2448-2455.

    [32] 張亭祿,徐青娜,陳樹果等.海水吸收光譜測(cè)量方法比較——定量濾膜技術(shù)和現(xiàn)場(chǎng)測(cè)量方法.中國(guó)海洋大學(xué)學(xué)報(bào),2011,41(7/8): 164-169.

    [33] 王桂芬,曹文熙,許大志等.南海北部水體浮游植物比吸收系數(shù)的變化.熱帶海洋學(xué)報(bào),2005,24(5): 1-10.

    [34] 徐京萍.石頭口門水庫水體光學(xué)特性及水色物質(zhì)遙感反演研究[學(xué)位論文].長(zhǎng)春:中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所,2009.

    [35] 袁宇翔,于洪賢,姜 明.小興凱湖浮游植物群落結(jié)構(gòu)及多樣性.濕地科學(xué),2013,11(2): 151-157.

    [36] Bricaud A, Babin M, Morel Aetal. Variability in the chlorophyll specific absorption coefficients of natural phytoplankton: Analysis and parameterization.JournalofGeophysicalResearch, 1995, 100(C7): 13321-13332.

    [37] Kirk JTO. Light and photosynthesis in aquatic ecosystems. Cambridge: Cambridge University Press, 1994: 1-431.

    [39] Binding CE, Jerome JH, Bukata RPetal. Spectral absorption properties of dissolved and particulate matter in Lake Erie.RemoteSensingofEnvironment, 2008, 112(4): 1702-1711.

    [40] 孫德勇,李云梅,王 橋等.太湖水體中懸浮顆粒物的比吸收光譜特性及其參數(shù)化分析.環(huán)境科學(xué)學(xué)報(bào),2009,29(9): 1820-1829.

    [41] 邢小罡,趙冬至,劉玉光等.渤海非色素顆粒物和黃色物質(zhì)的吸收特性研究.海洋環(huán)境科學(xué),2008,27(6): 595-598.

    [42] Song K, Liu D, Li Letal. Spectral absorption properties of colored dissolved organic matter(CDOM) and total suspended matter(TSM) of inland waters.ProceedingsoftheInternationalSocietyforOpticalEngineering, 2010, 7811: 78110B.

    [43] Zhang YL, Yin Y, Feng LQetal. Characterizing chromophoric dissolved organic matter in Lake Tianmuhu and its catchment basin using excitation emission matrix fluorescence and parallel factor analysis.WaterResearch, 2011, 45(16): 5110-5122.

    [44] Del Castillo CE, Coble PG. Seasonal variability of the colored dissolved organic matter during the 1994-95 NE and SW monsoons in the Arabian Sea.DeepSeaResearchPartII: Topical Studies in Oceanography, 2000, 47(7): 1563-1579.

    [45] De Haan H. Solar UV-light penetration and photodegradation of humic substances in peaty lake water.LimnologyandOceanography, 1993, 38(5): 1072-1076.

    [46] Ma R. Tang J, Dai Jetal. Absorption and scattering properties of water body in Taihu Lake, China: absorption.InternationalJournalofRemoteSensing, 2006, 27(19): 4277-4304.

    [47] Kowalczuk PA. Stedmon C, Markager S. Modeling absorption by CDOM in the Baltic Sea from salinity and chlorophyll.MarineChemistry, 2006, 101(1): 1-11.

    [48] Spencer RGM, Butler KD, Aiken GR. Dissolved organic carbon and chromophoric dissolved organic matter properties of rivers in the USA.JournalofGeophysicalResearch, 2012, 117: G3.

    [49] Zhang Y, Qin B, Zhu Getal. Chromophoric dissolved organic matter(CDOM) absorption characteristics in relation to fluorescence in Lake Taihu, China, a large shallow subtropical lake.Hydrobiologia, 2007, 581(1): 43-52.

    [50] Zhang Y, Qin B, Zhang Letal. Spectral absorption and fluorescence of chromophoric dissolved organic matter in shallow lakes in the middle and lower reaches of the Yangtze River.JournalofFreshwaterEcology, 2005, 20(3): 451-459.

    [51] 劉忠華,李云梅,呂 恒等.太湖春季水體固有光學(xué)特性及其對(duì)遙感反射率變化的影響.生態(tài)學(xué)報(bào),2012,32(2): 438-447.

    Absorption characteristics of particulates and CDOM in spring in the Lake Xingkai

    LI Sijia1,2, SONG Kaishan2, CHEN Zhiwen1, ZHAO Ying2, SHAO Tiantian2, MU Guangyi2,3& GUAN Ying2

    (1:InstituteofEcologicalEnvironment,JilinNormalUniversity,Siping136000,P.R.China)

    (2:NortheastInstituteofGeographyandAgriculturalEcology,ChineseAcademyofSciences,Changchun130102,P.R.China)

    (3:LifeScienceCollege,NortheastNormalUniversity,Changchun130000,P.R.China)

    Field surveys were carried out in the Lake Xingkai(including both the Great Lake Xingkai and the Small Lake Xingkai) in May 2013, and water samples were collected for examination of the absorption features. In the laboratory, tripton and chromophoric dissolved organic matter(CDOM) absorption of 20 water samples were measured through quantitative filter technique. Absorption coefficients were used to analyse the absorption characteristics, source, spatial distribution and relative contribution of optical active constituents over the range of PAR(400-700nm). The results showed that the absorption spectra of tripton are similar to nonalgal particles, indicating less phytoplankton presence in the water column. The total particulate absorption is dominated by the nonalgal particles, of which the contribution rate is up to 50% of the total absorption, andSgis higher than that in water bodies being investigated elsewhere. Positive correlation of the particles was observed withinap(440),ad(440) and total suspended matter, inorganic suspended matter and organic suspended matter, respectively, while weak correlation was recorded in terms of chlorophyll-a(Chl.a) compared with results reported in literatures. In the current study, low correlations of the particles were observed withinag(440), Chl.a and dissolve dissolved organic carbon(DOC). This is probably due to the fact that more colorless DOCs are present in the Lake Xingkai. Absorption coefficient and water quality parameters showed insignificant difference between the Great Lake Xingkai and the Small Lake Xingkai, though the latter is greatly impacted by farmland irrigation discharge, fisheries and tourism sewage.

    Lake Xingkai; absorption spectra; CDOM; nonalgal particles; phytoplankton

    *性國(guó)家自然科學(xué)基金項(xiàng)目(44171293)和興凱湖生態(tài)環(huán)境保護(hù)基線調(diào)查項(xiàng)目(JXZFCG2013-XXK21-03)聯(lián)合資助.2014-05-26收稿;2014-12-18收修改稿.李思佳(1988~),女,碩士研究生;E-mail: jlsdlsj@126.com.

    J.LakeSci.(湖泊科學(xué)), 2015, 27(5): 941-952

    DOI 10.18307/2015.0522

    ?2015 byJournalofLakeSciences

    **通信作者;E-mail: songks@iga.ac.cn.

    猜你喜歡
    興凱湖吸收光譜點(diǎn)位
    1年期國(guó)債收益率跌至1%
    銀行家(2025年1期)2025-02-08 00:00:00
    興凱湖霧凇景觀
    機(jī)器人快速示教方法及示教點(diǎn)位姿變換的研究
    穆棱河-興凱湖平原土地質(zhì)量地球化學(xué)評(píng)價(jià)
    原子吸收光譜分析的干擾與消除應(yīng)用研究
    淺析原子吸收光譜法在土壤環(huán)境監(jiān)測(cè)中的應(yīng)用
    垂直面內(nèi)建立基線的特殊點(diǎn)位高程測(cè)量法
    美哉,興凱湖霧凇
    茶油氧化過程中紫外吸收光譜特性
    那棵樹
    雪花(2015年6期)2015-12-10 23:48:35
    国产探花极品一区二区| 有码 亚洲区| 亚洲欧美日韩卡通动漫| 欧美日韩精品成人综合77777| 性色av一级| 中文字幕亚洲精品专区| 欧美日韩精品成人综合77777| 男女无遮挡免费网站观看| 亚洲欧美中文字幕日韩二区| 女人精品久久久久毛片| 国产男人的电影天堂91| 最新中文字幕久久久久| 国产色婷婷99| 亚洲精品久久午夜乱码| 黄色配什么色好看| 免费在线观看成人毛片| 啦啦啦中文免费视频观看日本| 国产黄色视频一区二区在线观看| 欧美日韩视频精品一区| 99热这里只有是精品50| 熟女人妻精品中文字幕| 精品久久国产蜜桃| 成人黄色视频免费在线看| 免费黄色在线免费观看| 午夜视频国产福利| 蜜桃久久精品国产亚洲av| 免费看日本二区| 嫩草影院入口| 一级二级三级毛片免费看| 久久久久久久精品精品| 久久青草综合色| 国产精品国产三级专区第一集| 国产淫语在线视频| 久久久久久久精品精品| 欧美+日韩+精品| 丝袜在线中文字幕| www.色视频.com| 国产片特级美女逼逼视频| 欧美xxⅹ黑人| 午夜免费观看性视频| 少妇人妻 视频| 色94色欧美一区二区| 另类精品久久| 九九在线视频观看精品| 成人黄色视频免费在线看| www.av在线官网国产| 日日撸夜夜添| 国产永久视频网站| 欧美三级亚洲精品| 大又大粗又爽又黄少妇毛片口| 晚上一个人看的免费电影| 亚洲av成人精品一区久久| 国产一区有黄有色的免费视频| 久久久久久久久久成人| 久久久久久久久久成人| 天天操日日干夜夜撸| 久久久久国产网址| 精品人妻一区二区三区麻豆| 黑人猛操日本美女一级片| 欧美人与善性xxx| 妹子高潮喷水视频| 日韩 亚洲 欧美在线| 欧美日韩av久久| 少妇猛男粗大的猛烈进出视频| 亚洲欧美成人综合另类久久久| 国产91av在线免费观看| 国产精品国产三级国产专区5o| 国产精品久久久久久久久免| 伦理电影免费视频| 精品少妇黑人巨大在线播放| 国产成人精品一,二区| 精品一品国产午夜福利视频| 国产伦精品一区二区三区四那| 美女视频免费永久观看网站| 人妻少妇偷人精品九色| 2021少妇久久久久久久久久久| 男人和女人高潮做爰伦理| 黄色毛片三级朝国网站 | 两个人免费观看高清视频 | 精品久久久久久久久av| 亚洲av国产av综合av卡| av福利片在线观看| 午夜福利网站1000一区二区三区| 午夜91福利影院| 亚洲国产精品专区欧美| 亚洲欧美一区二区三区黑人 | 久久午夜福利片| 日韩在线高清观看一区二区三区| 热re99久久精品国产66热6| 亚洲图色成人| 精华霜和精华液先用哪个| 嘟嘟电影网在线观看| 一级a做视频免费观看| 久久亚洲国产成人精品v| 亚洲人成网站在线观看播放| 韩国高清视频一区二区三区| 麻豆精品久久久久久蜜桃| 久久久久人妻精品一区果冻| 99久久精品热视频| 亚洲第一区二区三区不卡| 青春草国产在线视频| 亚洲精品国产色婷婷电影| 在线观看三级黄色| 中国三级夫妇交换| 久久热精品热| 中文在线观看免费www的网站| 91精品国产国语对白视频| 一级片'在线观看视频| 热99国产精品久久久久久7| 成人漫画全彩无遮挡| 久久久久视频综合| 亚洲精品aⅴ在线观看| 国内精品宾馆在线| 国产熟女午夜一区二区三区 | 插逼视频在线观看| av又黄又爽大尺度在线免费看| 我要看日韩黄色一级片| 国产中年淑女户外野战色| 少妇人妻精品综合一区二区| av国产精品久久久久影院| 色视频在线一区二区三区| 9色porny在线观看| a级毛片在线看网站| 日韩av免费高清视频| 午夜激情福利司机影院| 日本黄色日本黄色录像| 自拍偷自拍亚洲精品老妇| 亚洲av国产av综合av卡| 街头女战士在线观看网站| 大码成人一级视频| 丝袜在线中文字幕| 精品国产一区二区三区久久久樱花| 国产精品一二三区在线看| 国产精品久久久久久精品电影小说| 亚洲成人一二三区av| 国产91av在线免费观看| 成人18禁高潮啪啪吃奶动态图 | 久久狼人影院| 国产无遮挡羞羞视频在线观看| 91成人精品电影| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久精品古装| 99热这里只有精品一区| 综合色丁香网| 国产亚洲最大av| 青春草视频在线免费观看| 欧美日韩亚洲高清精品| 久久久久久久亚洲中文字幕| 国产日韩欧美亚洲二区| 国产伦在线观看视频一区| 国产黄片美女视频| a级毛片在线看网站| 伦理电影免费视频| 涩涩av久久男人的天堂| 亚洲激情五月婷婷啪啪| 久久青草综合色| 精品亚洲成a人片在线观看| av福利片在线观看| 我要看日韩黄色一级片| 自线自在国产av| 在线看a的网站| 偷拍熟女少妇极品色| 久久精品国产亚洲av天美| 午夜日本视频在线| 男女免费视频国产| 国产精品一区二区在线不卡| 新久久久久国产一级毛片| 永久免费av网站大全| 国产高清国产精品国产三级| av黄色大香蕉| 男人舔奶头视频| 精品少妇久久久久久888优播| 建设人人有责人人尽责人人享有的| 十八禁网站网址无遮挡 | 丝袜喷水一区| 国产成人精品福利久久| 亚洲av欧美aⅴ国产| 亚洲精品乱码久久久久久按摩| 一级,二级,三级黄色视频| 性高湖久久久久久久久免费观看| 亚洲av福利一区| 日本色播在线视频| 亚洲精品乱码久久久久久按摩| 我的女老师完整版在线观看| 青春草亚洲视频在线观看| 嫩草影院新地址| 一级av片app| av有码第一页| 成人黄色视频免费在线看| 国产欧美另类精品又又久久亚洲欧美| 日本猛色少妇xxxxx猛交久久| 日产精品乱码卡一卡2卡三| 成人国产麻豆网| 18禁在线无遮挡免费观看视频| 亚洲人与动物交配视频| 亚洲成人一二三区av| 中文字幕久久专区| 国产白丝娇喘喷水9色精品| 亚洲欧洲国产日韩| 好男人视频免费观看在线| 狠狠精品人妻久久久久久综合| 人妻 亚洲 视频| 女性生殖器流出的白浆| 国产av精品麻豆| 欧美精品一区二区大全| 国产精品国产三级专区第一集| 国产精品免费大片| 国产黄片视频在线免费观看| 亚洲精品成人av观看孕妇| 亚洲av在线观看美女高潮| 成年美女黄网站色视频大全免费 | 嫩草影院入口| 久久亚洲国产成人精品v| 69精品国产乱码久久久| 国产精品久久久久成人av| 另类亚洲欧美激情| 在线观看三级黄色| 亚洲精品一区蜜桃| 中文欧美无线码| 春色校园在线视频观看| 国产精品一区二区三区四区免费观看| 亚洲精品国产色婷婷电影| 久久久久久久久久久丰满| 交换朋友夫妻互换小说| 人妻人人澡人人爽人人| 久久精品国产亚洲av涩爱| 亚洲国产欧美在线一区| 免费看av在线观看网站| 交换朋友夫妻互换小说| 王馨瑶露胸无遮挡在线观看| 午夜福利,免费看| av免费在线看不卡| 妹子高潮喷水视频| 久久免费观看电影| 婷婷色综合www| 亚洲久久久国产精品| 自拍欧美九色日韩亚洲蝌蚪91 | av免费在线看不卡| 国产欧美日韩一区二区三区在线 | 成人国产av品久久久| 久久韩国三级中文字幕| 少妇的逼水好多| 91aial.com中文字幕在线观看| 午夜福利影视在线免费观看| 亚洲中文av在线| 日本猛色少妇xxxxx猛交久久| 熟妇人妻不卡中文字幕| 久久99蜜桃精品久久| 日韩免费高清中文字幕av| 午夜免费观看性视频| 欧美日韩精品成人综合77777| 国产 精品1| 成人免费观看视频高清| 成人18禁高潮啪啪吃奶动态图 | 少妇被粗大的猛进出69影院 | 色94色欧美一区二区| 搡老乐熟女国产| 在线观看人妻少妇| 99久久精品热视频| 下体分泌物呈黄色| 国产在线男女| 香蕉精品网在线| 亚洲精品国产av成人精品| 人妻人人澡人人爽人人| 精品人妻一区二区三区麻豆| 成人美女网站在线观看视频| 午夜影院在线不卡| 日本91视频免费播放| 大又大粗又爽又黄少妇毛片口| 欧美+日韩+精品| av又黄又爽大尺度在线免费看| 在线亚洲精品国产二区图片欧美 | 国产免费福利视频在线观看| 国精品久久久久久国模美| 赤兔流量卡办理| 亚洲美女视频黄频| 免费看日本二区| 精品人妻偷拍中文字幕| 91精品伊人久久大香线蕉| 久久毛片免费看一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品国产精品| 99热国产这里只有精品6| 麻豆成人午夜福利视频| 中文字幕人妻丝袜制服| 性色avwww在线观看| 欧美激情国产日韩精品一区| 夫妻午夜视频| 99热这里只有是精品50| 亚洲国产av新网站| 精品一品国产午夜福利视频| 日本午夜av视频| 亚洲人与动物交配视频| 天堂中文最新版在线下载| 国产黄片视频在线免费观看| 99re6热这里在线精品视频| 一级a做视频免费观看| 在线观看人妻少妇| 桃花免费在线播放| 美女视频免费永久观看网站| 日本av免费视频播放| 欧美日韩av久久| 亚洲电影在线观看av| 免费观看无遮挡的男女| 黄色配什么色好看| 天美传媒精品一区二区| 亚洲国产精品999| 免费av中文字幕在线| 精品亚洲成国产av| 我的女老师完整版在线观看| 涩涩av久久男人的天堂| 七月丁香在线播放| 久热久热在线精品观看| 亚洲欧美精品专区久久| 欧美3d第一页| 国产一区二区三区av在线| 国产在视频线精品| 国产男女超爽视频在线观看| 91精品伊人久久大香线蕉| 亚洲在久久综合| 青春草亚洲视频在线观看| 久久这里有精品视频免费| 日韩中文字幕视频在线看片| 新久久久久国产一级毛片| 视频区图区小说| 热re99久久精品国产66热6| 男人舔奶头视频| av国产久精品久网站免费入址| 内射极品少妇av片p| 亚洲国产精品一区三区| 亚洲国产精品999| 精品一区在线观看国产| 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| 久久狼人影院| 亚洲欧洲日产国产| 精品一品国产午夜福利视频| 国国产精品蜜臀av免费| 精品视频人人做人人爽| 国产精品欧美亚洲77777| 美女中出高潮动态图| 一级毛片 在线播放| 成年av动漫网址| 少妇人妻一区二区三区视频| 亚洲av综合色区一区| 中文字幕人妻丝袜制服| 观看美女的网站| 国产无遮挡羞羞视频在线观看| 午夜免费男女啪啪视频观看| 九九在线视频观看精品| 国产高清有码在线观看视频| 欧美日本中文国产一区发布| 校园人妻丝袜中文字幕| www.av在线官网国产| 久久99热这里只频精品6学生| a级片在线免费高清观看视频| 亚洲国产欧美在线一区| 国产一区有黄有色的免费视频| 极品人妻少妇av视频| 久久精品国产亚洲av天美| 国内揄拍国产精品人妻在线| 久久精品国产鲁丝片午夜精品| 国内揄拍国产精品人妻在线| 国产在线一区二区三区精| 精品国产一区二区久久| 精品亚洲成a人片在线观看| 美女国产视频在线观看| 国产日韩欧美亚洲二区| 亚洲精品色激情综合| 91成人精品电影| 99久久中文字幕三级久久日本| 国产在线一区二区三区精| 日本wwww免费看| 国产91av在线免费观看| 亚洲国产精品专区欧美| 又黄又爽又刺激的免费视频.| 亚洲av.av天堂| 丁香六月天网| 亚州av有码| 在线观看人妻少妇| 午夜激情久久久久久久| 亚洲精华国产精华液的使用体验| 99热国产这里只有精品6| 久久人人爽人人片av| av网站免费在线观看视频| 精品人妻熟女av久视频| 蜜臀久久99精品久久宅男| 久久精品夜色国产| av播播在线观看一区| 国产精品久久久久久精品古装| 亚洲熟女精品中文字幕| 国产无遮挡羞羞视频在线观看| 国产 一区精品| 男人狂女人下面高潮的视频| 人人妻人人澡人人看| 国产精品蜜桃在线观看| 精品一区二区免费观看| 欧美精品高潮呻吟av久久| 爱豆传媒免费全集在线观看| 精品久久久噜噜| 日韩大片免费观看网站| 精品视频人人做人人爽| 我要看日韩黄色一级片| 亚洲欧美日韩东京热| 啦啦啦中文免费视频观看日本| 岛国毛片在线播放| 精品久久久精品久久久| 纵有疾风起免费观看全集完整版| 麻豆成人午夜福利视频| 国产黄色免费在线视频| 热re99久久精品国产66热6| 亚洲精品日韩av片在线观看| 欧美成人午夜免费资源| 一本色道久久久久久精品综合| 久久国产精品男人的天堂亚洲 | 亚洲精品,欧美精品| 美女视频免费永久观看网站| 亚洲国产av新网站| 久久国产精品男人的天堂亚洲 | 日本vs欧美在线观看视频 | 97在线视频观看| 国产av一区二区精品久久| 亚洲av电影在线观看一区二区三区| 熟女av电影| 人人澡人人妻人| 十八禁网站网址无遮挡 | 人妻夜夜爽99麻豆av| 男人爽女人下面视频在线观看| 男女边摸边吃奶| 久久精品久久久久久久性| 精品久久久久久久久亚洲| 国产深夜福利视频在线观看| 五月玫瑰六月丁香| 国产成人午夜福利电影在线观看| 在线观看美女被高潮喷水网站| 黄色一级大片看看| 国语对白做爰xxxⅹ性视频网站| 你懂的网址亚洲精品在线观看| 午夜精品国产一区二区电影| 亚洲精品视频女| 麻豆乱淫一区二区| 深夜a级毛片| 99热这里只有是精品在线观看| 男的添女的下面高潮视频| 久久精品国产亚洲网站| 日韩伦理黄色片| 国产爽快片一区二区三区| 亚洲av.av天堂| 大又大粗又爽又黄少妇毛片口| 激情五月婷婷亚洲| 亚洲国产日韩一区二区| 好男人视频免费观看在线| 成年美女黄网站色视频大全免费 | 青春草视频在线免费观看| 大片免费播放器 马上看| 免费观看在线日韩| 亚洲人与动物交配视频| 69精品国产乱码久久久| 国产亚洲av片在线观看秒播厂| 少妇的逼好多水| 熟女av电影| 一二三四中文在线观看免费高清| 熟女电影av网| 欧美日韩av久久| 国产精品福利在线免费观看| 成年av动漫网址| 亚洲精品色激情综合| 久久狼人影院| 一级二级三级毛片免费看| 伦理电影大哥的女人| 六月丁香七月| 大话2 男鬼变身卡| 丰满饥渴人妻一区二区三| 26uuu在线亚洲综合色| 免费黄频网站在线观看国产| 黑人高潮一二区| 日韩精品免费视频一区二区三区 | 美女大奶头黄色视频| 另类精品久久| 亚洲伊人久久精品综合| 国产伦精品一区二区三区四那| 观看美女的网站| 亚洲成人av在线免费| 国产精品一区www在线观看| 欧美日韩在线观看h| 久久久久久人妻| 亚洲欧美成人综合另类久久久| 丝瓜视频免费看黄片| 国产精品人妻久久久影院| 国产在线男女| 久久国产亚洲av麻豆专区| 寂寞人妻少妇视频99o| 午夜激情福利司机影院| 视频区图区小说| 成人综合一区亚洲| 人人妻人人添人人爽欧美一区卜| 老司机影院毛片| 日韩亚洲欧美综合| 久久免费观看电影| 久久久久久人妻| 久久人人爽人人爽人人片va| 内射极品少妇av片p| 人妻人人澡人人爽人人| 国产精品国产三级专区第一集| 蜜桃久久精品国产亚洲av| 97超视频在线观看视频| 日日摸夜夜添夜夜添av毛片| 国产成人a∨麻豆精品| 两个人的视频大全免费| 妹子高潮喷水视频| 女人精品久久久久毛片| 国产精品不卡视频一区二区| 亚洲成人手机| av网站免费在线观看视频| 香蕉精品网在线| 日韩中文字幕视频在线看片| 高清视频免费观看一区二区| 亚洲精品国产av成人精品| 日产精品乱码卡一卡2卡三| 亚洲av在线观看美女高潮| 97在线视频观看| 大香蕉久久网| 五月天丁香电影| 成人黄色视频免费在线看| 亚洲国产毛片av蜜桃av| 综合色丁香网| 国产免费福利视频在线观看| 男人狂女人下面高潮的视频| 精品午夜福利在线看| 亚洲久久久国产精品| 免费观看的影片在线观看| 男女边摸边吃奶| 欧美成人午夜免费资源| 欧美另类一区| 国产精品.久久久| 精品人妻熟女av久视频| 成人国产麻豆网| 午夜免费观看性视频| 黄色毛片三级朝国网站 | 91精品国产九色| 丰满人妻一区二区三区视频av| 尾随美女入室| 一区二区三区乱码不卡18| 免费高清在线观看视频在线观看| 精品99又大又爽又粗少妇毛片| 亚洲av福利一区| 欧美日本中文国产一区发布| 国产淫片久久久久久久久| 亚洲熟女精品中文字幕| 精品一区二区三卡| 免费观看无遮挡的男女| 天天躁夜夜躁狠狠久久av| 亚洲精品乱码久久久v下载方式| 欧美国产精品一级二级三级 | 亚洲国产色片| 3wmmmm亚洲av在线观看| 日韩欧美一区视频在线观看 | 国产成人精品久久久久久| 色视频www国产| 日韩三级伦理在线观看| 欧美日本中文国产一区发布| 欧美日韩av久久| 搡女人真爽免费视频火全软件| 秋霞伦理黄片| 美女福利国产在线| 精品人妻偷拍中文字幕| 亚洲精品一区蜜桃| 亚洲国产精品一区三区| 天堂中文最新版在线下载| 国产精品久久久久久久久免| 精品久久久久久久久av| 人人妻人人澡人人爽人人夜夜| 免费黄色在线免费观看| 七月丁香在线播放| 在线精品无人区一区二区三| 亚洲欧美清纯卡通| 黄片无遮挡物在线观看| 中文字幕久久专区| 午夜福利视频精品| 亚洲精品久久久久久婷婷小说| 一个人免费看片子| 在线观看美女被高潮喷水网站| 亚洲av日韩在线播放| 亚洲人成网站在线播| 国产精品欧美亚洲77777| 色网站视频免费| 国产无遮挡羞羞视频在线观看| 777米奇影视久久| 欧美日本中文国产一区发布| 国产美女午夜福利| 一级毛片 在线播放| 精品人妻熟女毛片av久久网站| 国产精品嫩草影院av在线观看| 国产熟女午夜一区二区三区 | 边亲边吃奶的免费视频| av线在线观看网站| 久久午夜综合久久蜜桃| 亚洲欧美精品自产自拍| 成年av动漫网址| 在线观看国产h片| 国产成人精品久久久久久| 中文字幕免费在线视频6| 狂野欧美激情性xxxx在线观看| 人体艺术视频欧美日本| 丝袜在线中文字幕| 免费av中文字幕在线| 一级毛片黄色毛片免费观看视频| 超碰97精品在线观看| 色哟哟·www| 老司机影院毛片| 国产精品伦人一区二区| 97在线视频观看| av.在线天堂| 久久久久久人妻| 9色porny在线观看|