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

    遼河水體光學(xué)吸收特性的季節(jié)變化

    2016-04-11 03:31:01邵田田宋開山
    生態(tài)學(xué)報(bào) 2016年7期
    關(guān)鍵詞:季節(jié)變化遼河顆粒物

    邵田田,宋開山,丁 智,趙 瑩,王 銘,杜 嘉,張 柏

    1 中國科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所, 長春 130102 2 中國科學(xué)院大學(xué), 北京 100049 3 中國科學(xué)院地理科學(xué)與資源研究所, 北京 100101

    ?

    遼河水體光學(xué)吸收特性的季節(jié)變化

    邵田田1,2,宋開山1,*,丁智3,趙瑩1,王銘1,杜嘉1,張柏1

    1 中國科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所, 長春130102 2 中國科學(xué)院大學(xué), 北京100049 3 中國科學(xué)院地理科學(xué)與資源研究所, 北京100101

    摘要:吸收特性是水體重要的光學(xué)特性,也是建立生物光學(xué)模型的基本參數(shù)。利用2013年5月和9月遼河流域水體實(shí)測數(shù)據(jù)對總懸浮顆粒物、浮游植物、非藻類顆粒物以及有色溶解有機(jī)物(CDOM)的吸收特性、季節(jié)變化進(jìn)行了研究。研究表明,總顆粒物的吸收光譜曲線與非藻類顆粒物的吸收光譜基本保持一致并且與非藻類顆粒物的吸收[ad(λ)]以及浮游植物的吸收[aph(λ)]都呈現(xiàn)顯著相關(guān)。兩個(gè)季節(jié),ad(440)均是總吸收的主要組成部分,CDOM的吸收[aCDOM(440)]則大于aph(440)的貢獻(xiàn)率,但9月ad(440)的吸收所占比例較5月高。遼河水體非藻類顆粒物中礦物顆粒物的含量較其他一般水體要高。通過分析5月及9月的aph(440)/aph(675)發(fā)現(xiàn)浮游植物中輔助色素與葉綠素a濃度(Chla)的組成季節(jié)性變化不大,但其空間特征表現(xiàn)較為明顯。另外,CDOM的吸收特性與非藻類顆粒物的吸收相似,呈現(xiàn)出較大的空間差異性;遼河SCDOM基本小于其它河流、湖泊、河口等水體,其CDOM主要由大分子物質(zhì)組成。

    關(guān)鍵詞:遼河;吸收特性;顆粒物;CDOM;季節(jié)變化

    水體光學(xué)活性物質(zhì)(OpticallyActiveConstituents,OACs)的固有光學(xué)特性(InherentOpticalProperties,IOPs)具有不受光照條件影響的優(yōu)點(diǎn),是基于水體輻射傳輸理論的生物光學(xué)模型反演各種水質(zhì)參數(shù)的基礎(chǔ),具有嚴(yán)格的物理意義,在水色遙感中有著重要作用[1-2]。

    內(nèi)陸水體的吸收系數(shù)和后向散射系數(shù)是水體的兩個(gè)主要固有光學(xué)參量[3],其中固有光學(xué)吸收特性受4種物質(zhì)影響:純水、浮游植物、非藻類顆粒物和有色溶解有機(jī)物(ChromophoricDissolvedOrganicMatter,CDOM)。在指定波段,水體的總吸收系數(shù)等于各組分在該波段的吸收系數(shù)之和[4],即a=aw+aph+ad+aCDOM,其中aw、aph、ad和aCDOM分別代表4種相應(yīng)物質(zhì)的吸收系數(shù)。純水對光的衰減是所有水體都共有的[5],因此水體固有光學(xué)吸收特性主要受制于非藻類顆粒物、浮游植物和CDOM的光譜吸收,而作為水體的固有光學(xué)特性,三者是水色遙感的生物-光學(xué)模式、光輻射傳輸、初級生產(chǎn)力估算和生態(tài)過程模擬的基本參數(shù)[5]。不同于Ⅰ類大洋水體,在近岸的Ⅱ類水體和內(nèi)陸的湖泊水體中,由于陸源的非藻類顆粒物和CDOM的大量輸入,導(dǎo)致水體的固有光學(xué)特性不再僅由浮游植物所決定[6-7]。大量研究都比較集中于海洋Ⅰ類水體、國外內(nèi)陸水體以及國內(nèi)富營養(yǎng)化較為嚴(yán)重的長江中下游的內(nèi)陸湖泊[5, 8-9],少量研究集中于石頭口門水庫[10]、官廳水庫[7]、千島湖[11]等污染較輕或者很純凈的水體,而對于河流水體固有光學(xué)吸收特性的研究相對較少[12]。河流作為連接大洋水體和內(nèi)陸水體的紐帶,其固有光學(xué)特性不同于內(nèi)陸湖泊,也不同于近岸II類水體。

    遼河水系流經(jīng)沈陽、撫順、鞍山、本溪等大中型工業(yè)城市,隨著工農(nóng)業(yè)和大城市的迅速發(fā)展,人為因素對遼河流域水環(huán)境的影響越來越大,所產(chǎn)生的危害也逐年增加,使遼河處于較為嚴(yán)重的污染狀態(tài)[13- 14],從而影響水體的固有光學(xué)吸收特性。因此作為東北地區(qū)南部向海輸入的主要途徑,河流中的總懸浮顆粒物吸收(非藻類顆粒物和浮游植物)、CDOM的吸收將會(huì)影響遼河口以及渤海灣水體的固有特性,為更好的全面理解內(nèi)陸水體提供支持?jǐn)?shù)據(jù),也為生物光學(xué)模型的發(fā)展和完善提供重要借鑒。

    1材料與方法

    1.1研究區(qū)概況

    遼河流域位于中國東北地區(qū)南部(116° 30′—125° 47′E, 38° 43′—45°N),流域面積為22.9×104km2。流域由2個(gè)獨(dú)立水系所組成:一為東、西遼河,于福德店匯流為遼河干流,經(jīng)雙臺子河由盤山入海,全長1.39×103km;另一水系為渾河、太子河于三岔河匯合后經(jīng)大遼河由營口入海,全長415.4km。該區(qū)氣候?qū)儆诎霛駶?半干旱的溫帶氣候區(qū),溫度由南向北降低(4—7 °C),而降水量則自東向西遞減(350—1000mm)。遼河流域包括東、西遼河流域、遼河干流流域以及渾太河流域。而本文的研究區(qū)同時(shí)還包括了大凌河流域,主要是因?yàn)榇罅韬釉谶|東灣注入渤海,與遼河流域相接,聯(lián)系密切。因此,在本文中所涉及到的遼河流域均是指包括大凌河在內(nèi)的區(qū)域。大凌河位于遼寧省西部,是遼寧省西部最大的河流,全長397km,流域面積2.35×104km2。

    1.2野外樣品的采集

    于2013年5月(春季)、2013年9月(秋季)對遼河流域主要河流相同位置點(diǎn)進(jìn)行重復(fù)采樣,每次采集40組,共計(jì)80組數(shù)據(jù)進(jìn)行遼河水體固有光學(xué)吸收特性及季節(jié)性變化的分析,具體采樣點(diǎn)位置如表1所示。

    在采樣過程中,每個(gè)樣點(diǎn)采集2.5L水樣并立刻放在車載冰箱里冷藏,送至實(shí)驗(yàn)室測量其光學(xué)參數(shù)等,整個(gè)測試過程需要2—3d。

    1.3固有光學(xué)特性物質(zhì)吸收的測定

    總懸浮物的吸收系數(shù)采用定量濾膜技術(shù)(QFT:quantitativefilttertechnique)進(jìn)行測定[15],用直徑25mm的GF/F(Whatman)濾膜過濾一定體積的水樣,在UV-2600紫外分光光度計(jì)下測定其380—800nm吸光度,并將各波段的吸光度減去750nm波長處的吸光度來消除儀器間的差異,然后根據(jù)式(1)計(jì)算總懸浮物的吸收系數(shù):

    (1)

    式中,ap(λ)為λnm處總懸浮物的吸收系數(shù),V為過濾的水樣的體積(m3),S為過濾留在濾膜上的有效面積(m2)。

    非藻類顆粒物的吸收是利用UV-2600紫外分光光度計(jì)測得非藻類顆粒物的吸光度,然后根據(jù)式(1)計(jì)算其吸收系數(shù)。在非藻類顆粒物的吸收測定之前,需要用0.1%的次氯酸鈉浸泡濾膜,將濾膜上的色素溶解,并使用空白濾膜作為參比,測得其吸光度。浮游植物的吸收系數(shù)aph(λ)則是通過總懸浮物吸收系數(shù)減去非藻類顆粒物吸收系數(shù)所得到:

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

    (2)

    CDOM吸收系數(shù)的測定是通過0.22μm的微孔聚碳酸酯膜過濾水樣后,利用UV-2600紫外分光光度計(jì)測得200—800nm的吸光度,然后根據(jù)式(3)計(jì)算各波長的吸收系數(shù)。

    aCDOM(λ′)=2.303OD(λ)/l

    (3)

    式中,a(λ′)為未校正的吸收系數(shù)(m-1);D(λ)為吸光度;l為光程路徑(m)。為消除過濾液中殘留細(xì)小顆粒物的散射,本文用750nm處吸收系數(shù)進(jìn)行校正[16],得到校正后的吸收系數(shù)(m-1):

    aCDOM(λ)=aCDOM(λ′)-aCDOM(750)(λ/750)

    (4)

    1.4其他參數(shù)的獲取

    葉綠素a(Chlorophylla,Chla)濃度通過90%的丙酮溶液萃取,并用UV-2600紫外分光光度計(jì)分別測得630、647、664和750nm處的吸光度,計(jì)算出葉綠素a的濃度[17]??倯腋☆w粒物濃度(TSM)、無機(jī)顆粒物濃度(ISM)、有機(jī)顆粒物濃度(OSM)的測定采用稱重法獲得。溶解有機(jī)碳(DOC)濃度的測試基于總有機(jī)碳分析儀(日本島津,TOC-VCPN)獲得。

    1.5S值的獲取

    SCDOM是CDOM吸收曲線的斜率,可以反映CDOM組成分子的大小,表征CDOM的來源[18-19]。對于SCDOM的獲得主要基于吸收波譜波段以及擬合方法的選擇[20-21]。本研究選擇常用波段范圍280—400nm、280—320nm以及350—400nm,根據(jù)式(5)采用非線性擬合方法進(jìn)行S值的求算:

    aCDOM(λ)=aCDOM(λ0)exp[SCDOM(λ0-λ)]

    (5)

    式中,aCDOM(λ)是CDOM在特定波長的吸收系數(shù)(m-1);aCDOM(λ0)是參考波段(通常選擇440nm)的吸收系數(shù)(m-1)。

    另外,對于II類水體來說,可以用式(6)表示非藻類顆粒物在波長λ處的吸收系數(shù)[7]。

    ad(λ)=ad(λ0)exp[Sd(λ0-λ)]

    (6)

    其中,ad(λ)表示非藻類顆粒物在波長λ的吸收系數(shù)(m-1);λ0是參照波長,一般取440nm;Sd是指數(shù)函數(shù)曲線斜率參數(shù)(μm-1)。本文采用最小二乘法擬合,建立非藻類顆粒物在400—700nm吸收模型。

    2結(jié)果與分析

    2.1總懸浮顆粒物吸收特征

    春、秋季遼河采樣點(diǎn)的總懸浮顆粒物吸收系數(shù)(ap(λ))、非藻類顆粒物吸收系數(shù)(ad(λ))、浮游植物吸收系數(shù)(aph(λ)的吸收光譜曲線如圖1所示。

    同時(shí)對比圖1和表2可以發(fā)現(xiàn),同一季節(jié)不同位置點(diǎn)ap(λ)的差異很大;不同季節(jié)相同采樣點(diǎn)ap(λ)差異性也比較大,9月的吸收明顯高于5月。但是從曲線的形狀來看,不同季節(jié)所有采樣點(diǎn)的吸收系數(shù)大致隨波長的增加呈現(xiàn)逐漸降低的趨勢,但在675nm附近存在不太明顯的葉綠素a的吸收峰,這與浮游植物在該波段的特征吸收有關(guān);而在另一個(gè)浮游植物特征吸收波段440nm則觀察不到吸收峰,這主要是因?yàn)榉窃孱愵w粒物的吸收在總懸浮顆粒物的吸收中占有絕對優(yōu)勢(圖1),除675nm附近的浮游植物強(qiáng)吸收波段外,非藻類顆粒物吸收在其他波段的吸收要遠(yuǎn)遠(yuǎn)大于浮游植物,致使浮游植物的吸收特征在總懸浮物的吸收中表現(xiàn)不明顯。河流的流動(dòng)性影響沉積物發(fā)生再懸浮,使得本文樣品的總懸浮顆粒物的吸收曲線更接近于非藻類顆粒物的吸收,而這種顆粒物吸收類型的曲線經(jīng)常出現(xiàn)在無機(jī)顆粒物濃度高的河流、河口、沿岸帶等Ⅱ類水體[22-23]。

    本文將ap(440)、ap(675)、ad(440)、ad(675)、aph(400)、aph(675)以及aCDOM(440)與水質(zhì)參數(shù)TSM、ISM、OSM、Chla和DOC濃度進(jìn)行相關(guān)分析(表3)。可以看出,5月和9月,總懸浮顆粒物的吸收與非藻類顆粒物的吸收都具有很強(qiáng)的相關(guān)性,特別是在440nm,相關(guān)系數(shù)均達(dá)0.99;ap(440)與aph(440)雖呈顯著性相關(guān),但相關(guān)性較弱;而ap(675)與aph(675)相關(guān)系數(shù)分別為0.81(5月)和0.55(9月),呈顯著相關(guān);另外,ap(440)和ap(675)與TSM呈現(xiàn)顯著相關(guān),ap(675)與Chla也表現(xiàn)出顯著正相關(guān)??倯腋☆w粒物的吸收與ISM、OSM在9月份表現(xiàn)出顯著正相關(guān),而5月份不存在相關(guān)性。另外,兩個(gè)季節(jié),總顆粒物的吸收與DOC未表現(xiàn)出相關(guān)性。

    2.2非藻類顆粒物的吸收特征

    水體中非藻類懸浮顆粒物主要包括礦物沉積物、非活性有機(jī)體(如排泄物、退化的藻類細(xì)胞等)以及非藻類活性有機(jī)體[20],因此其吸收特性主要由礦物顆粒物與有機(jī)物所決定。遼河水體ad(λ)隨著波長的增加呈指數(shù)遞減(圖1),ad(440)和ad(675)的變化范圍呈現(xiàn)出較大的空間差異性(表2)。

    表2不同季節(jié)總懸浮物[totalsuspendedparticulate, ap(λ)]、非藻類顆粒物[non-algalparticles, ad(λ)]以及浮游植物[phytoplankton, aph(λ)]的統(tǒng)計(jì)分析

    Table2Statisticalanalysisofabsorptionoftotalsuspendedparticulate[ap(λ)],non-algalparticles[ad(λ)]andphytoplankton[aph(λ)]indifferentmonths

    ap(λ): 總懸浮物Totalsuspendedparticulate;ad(λ): 非藻類顆粒物non-algalparticlesabsorption;aph(λ): 浮游植物phytoplankton

    許多實(shí)測數(shù)據(jù)模擬表明,不同波段的非藻類吸收系數(shù)之間存在一定的關(guān)系。我們選擇在水色遙感中較為重要的440nm為參考波段,采用最小二乘法擬合,建立非藻類顆粒物的吸收模型(圖2)。可以看出,兩期數(shù)據(jù)典型樣點(diǎn)的模擬值與實(shí)測值之間的吻合度非常好,非線性擬合的決定系數(shù)均大于0.99,所有樣點(diǎn)決定系數(shù)均大于0.95。遼河水體5月份ad(λ)在400—700nm按指數(shù)衰減規(guī)律擬合得到斜率Sd值的變化范圍為10.40—13.24μm-1,均值為11.14μm-1;而9月份ad(λ)的Sd變化范圍8.51—13.06μm-1,均值為11.26μm-1。

    與國內(nèi)其它研究相比,遼河Sd高于官廳水庫(均值6.7μm-1)[7],與巢湖(均值11.9μm-1)[24]、太湖(均值11.7μm-1)[5]、三峽水庫(均值11.5μm-1)[25]、石頭口門水庫(均值11.78μm-1,11.38μm-1)[10]的Sd值相近,而小于滇池(均值12.5μm-1)[26]和珠江口(均值12.0μm-1)[23]。Babin等[27]認(rèn)為,光譜斜率Sd的大小可能與非藻類顆粒物中所含的有機(jī)物與礦物顆粒的比例有關(guān);而Binding等[20]則進(jìn)一步確定Sd的大小與礦物顆粒物濃度有關(guān),斜率越小,其濃度越高。由此可見,遼河水體非藻類顆粒物中礦物顆粒物的含量較其他一般水體要高。

    已有研究表明,對于大洋I類水體而言,由于其非藻類顆粒物主要來源于浮游植物的降解[23],ad(λ)與Chla濃度存在較好的相關(guān)關(guān)系。而對于遼河流動(dòng)水體來說,ad(λ)與Chla濃度幾乎不存在相關(guān)性,與ISM、OSM以及DOC的相關(guān)性較差,但與TSM存在顯著相關(guān)性(表3),結(jié)果與昆承湖的研究類似[28]。主要是因?yàn)閼腋∧嗌车汝懺葱詿o機(jī)質(zhì)輸入是遼河非藻類顆粒物的主要來源,而浮游植物的降解作用所占比例較小。另外本文對特定波段(400nm, 440nm, 510nm, 550nm, 675nm)的非藻類顆粒物的比吸收系數(shù)[ad(λ):TSM]進(jìn)行計(jì)算[23],結(jié)果如表4??梢钥闯?,ad(λ):TSM隨波長的增加而減小。遼河水體ad(λ):TSM的比值表現(xiàn)出較大的季節(jié)波動(dòng)性,9月份明顯高于5月份的均值,主要原因在于9月份ad(λ)較5月份大很多(表2,圖1),而TSM較為相近,導(dǎo)致兩者的比值差距較大。另外對比前人研究發(fā)現(xiàn),5月份遼河水體ad(λ):TSM近似于美國Erie湖的結(jié)果[23],小于太湖;而9月份較太湖高出很多。

    2.3浮游植物吸收特征

    遼河水體浮游植物的吸收曲線變化較大(圖1),反映了浮游植物吸收呈現(xiàn)較大差別的空間特征。另外,aph(λ)主要受水體中浮游藻類的濃度及種群組成的影響,隨葉綠素濃度的變化而變化。浮游植物在藍(lán)光波段的440nm附近和紅光波段的675nm附近存在兩個(gè)吸收峰(圖1),是浮游植物的兩個(gè)特征吸收波段,675nm的吸收是葉綠素a作用的結(jié)果,而440nm的吸收則是葉綠素a和其他輔助色素的共同作用結(jié)果。因此藍(lán)光波段與紅光波段的比值[aph(440)/aph(675)]能在一定程度上反映輔助色素與葉綠素a濃度的組成比例[29]。本文5月aph(440)/aph(675)的變化范圍為0.14—5.35m-1(均值為2.08m-1),9月aph(440)/aph(675)的變化范圍0.10—9.46m-1(均值為2.56m-1),說明浮游植物中輔助色素與葉綠素a濃度的組成季節(jié)性變化不大。

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

    ap(440), ap(675): 總懸浮物在440和675nm的吸收absorptionoftotalsuspendedparticulateat440and675nm; ad(440), ad(675): 非藻類顆粒物在440和675nm的吸收absorptionofnon-algalparticlesat440and675nm; aph(440), aph(675): 浮游植物在440和675nm的吸收absorptionofphytoplanktonat440and675nm; aCDOM(440): 有色溶解有機(jī)物在440nm的吸收absorptionofchromophoricdissolvedorganicmatter(CDOM);TSM: 總懸浮顆粒物濃度TotalSuspendedMatter;ISM: 無機(jī)顆粒物濃度InorganicSuspendedMatter;OSM: 有機(jī)顆粒物濃度OrganicSuspendedMatter;Chla: 葉綠素a濃度Chlorophylla;DOC: 溶解有機(jī)碳濃度DissolvedOrganicCarbon

    5月及9月aph(440)、aph(675)的變化范圍及均值如表2所示,5月Chla的變化范圍為0.001—21.42μg/L,9月Chla的變化范圍為2.94—138.84μg/L。而據(jù)表3可知,遼河葉綠素a濃度與浮游植物吸收系數(shù)具有顯著正相關(guān)關(guān)系。前人研究結(jié)果表明浮游植物吸收系數(shù)與葉綠素a濃度存在線性或非線性關(guān)系,本文研究也證實(shí)了這種關(guān)系的存在。

    2.4有色溶解有機(jī)物的吸收特征

    CDOM的光譜吸收特征與非藻類顆粒物的吸收相似,在700nm之后基本為0,而在短波波段(280—500nm)呈指數(shù)增長(圖3),并呈現(xiàn)出較大差異性。5月aCDOM(440)吸收范圍為0.37—2.21m-1,均值為0.94m-1,而9月份的變化范圍為0.33—3.09m-1(均值為1.03m-1),兩期數(shù)據(jù)aCDOM(440)的最大值均出現(xiàn)在s6(新開河,隸屬西遼河),而且濃度最高的幾個(gè)采樣點(diǎn)均位于西遼河流域,表現(xiàn)出較大的空間差異性。西遼河樣點(diǎn)濃度較高的主要原因是西遼河流域較為干旱,徑流量和降雨量較少,從而使河流中CDOM的濃度較高[30]。

    CDOM是溶解有機(jī)物中有顏色的部分,因而在很多研究中都可以發(fā)現(xiàn)CDOM在特定波段的吸收系數(shù)與DOC濃度存在非常顯著的正相關(guān)[17, 19]。本文將aCDOM(440)與DOC濃度進(jìn)行相關(guān)分析,發(fā)現(xiàn)兩者在5月和9月都存在顯著正相關(guān),分別為0.55和0.90(表3)。DOC組成可以分為有色部分和無色部分,而當(dāng)DOC中無色的部分為常數(shù)或變化不大時(shí),CDOM吸收系數(shù)則會(huì)與DOC濃度存在非常好的相關(guān)性[31]。遼河作為流動(dòng)水體,攜帶更多陸源的DOC,因而含有較多的有色DOC,無色DOC的含量相對較低,從而使遼河水體兩個(gè)季節(jié)CDOM吸收系數(shù)會(huì)與DOC濃度存在較強(qiáng)的相關(guān)性。而9月兩者的相關(guān)性較5月更高,原因在于2013年7、8、9月,遼河流域的強(qiáng)降雨導(dǎo)致河流徑流量增加,陸源DOC含量更高,因而9月相關(guān)性更強(qiáng)。

    另外,本文將aCDOM(440)與其他組分及濃度進(jìn)行相關(guān)分析發(fā)現(xiàn),除與DOC具有顯著正相關(guān)外,與其他所有參量幾乎都不存在相關(guān)性(表3)。aCDOM(440)與Chla之間較差的相關(guān)性說明遼河水體中自生CDOM的影響較小,而與總懸浮物的相關(guān)性說明外源物質(zhì)對CDOM的影響也是有限,從而說明河流CDOM的來源及作用較為復(fù)雜。

    SCDOM能夠反映水體中CDOM分子量的大小與來源,因此本文利用Matlab,基于最小二乘法對280—400nm、280—320nm以及350—400nm的吸收曲線進(jìn)行擬合求取斜率值,如圖4所示。可以看出,SCDOM隨采樣位置點(diǎn)的空間變化而發(fā)生變化,但是5月與9月的S280—400、S280—320以及S350—400值隨采樣點(diǎn)空間位置的變化較為一致,這說明遼河水體SCDOM受季節(jié)變化影響不大,跟采樣點(diǎn)所處的地理位置與環(huán)境條件有關(guān)。

    本文對三個(gè)波段指數(shù)函數(shù)斜率值進(jìn)行統(tǒng)計(jì)(表5)并與前人的研究進(jìn)行比較,發(fā)現(xiàn)結(jié)果與河流較為相近,而低于湖泊、河口的SCDOM。如Spencer等[32]對美國河流的CDOM進(jìn)行擬合發(fā)現(xiàn),大多數(shù)的河流S350—400值在0.0165—0.0183nm-1內(nèi)變化,而所有河流的最大值為0.020nm-1,最小值為0.0160nm-1,稍大于遼河水體;Chen等[12]發(fā)現(xiàn)長江S350-400值在0.0121—0.0236nm-1之間變化,均值為0.0148nm-1,低于遼河水體;Zhang等[33]發(fā)現(xiàn)長江中下游淺型湖泊S280—400的均值為0.02nm-1,大于本文5月份和9月份的值;劉忠華等[9]通過對太湖春季水體的研究發(fā)現(xiàn)其S280—320在0.0170—0.0268nm-1之間變化,均值為0.0203nm-1,明顯大于遼河。因此,相對于其他河流與湖泊,遼河CDOM組成中大分子較多,這與遼河流域水土流失嚴(yán)重,水體攜帶較多的陸源DOM關(guān)系密切。

    2.5光學(xué)活性物質(zhì)對總吸收的貢獻(xiàn)

    決定內(nèi)陸水體光吸收特性的物質(zhì)主要有純水、浮游植物、非藻類顆粒物、CDOM四種。在給定波段,內(nèi)陸水體總的吸收系數(shù)等于各組分在該波段各組分的吸收系數(shù)之和,其中,純水的吸收是固定的[34]。本文對440nm處各種組分對總吸收的貢獻(xiàn)進(jìn)行分析,由于純水在440nm處的吸收非常小,僅有0.00635m-1,因此文中所指的總吸收將忽略純水的吸收(圖5)。

    可以看出,兩個(gè)月份非藻類顆粒物的吸收在總吸收中占據(jù)主導(dǎo)地位。5月有38個(gè)采樣點(diǎn)的非藻類顆粒物吸收的貢獻(xiàn)率在40%以上,32個(gè)樣點(diǎn)的貢獻(xiàn)率在50%以上,平均貢獻(xiàn)率為67.46%;浮游植物的吸收對總吸收的貢獻(xiàn)率較非藻類顆粒物的吸收低很多,平均貢獻(xiàn)率僅為14.45%;CDOM的吸收對總吸收的貢獻(xiàn)遠(yuǎn)小于非藻類顆粒物,但大于浮游植物的吸收,平均貢獻(xiàn)率為18.09%。相對于5月,9月的非藻類顆粒物的吸收的貢獻(xiàn)率更高,除有一個(gè)采樣點(diǎn)的貢獻(xiàn)率為32.01%外,其余所有樣點(diǎn)的貢獻(xiàn)率都在50%以上,最高達(dá)到96.09%;浮游植物的吸收貢獻(xiàn)率基本上低于20%,平均貢獻(xiàn)率為9.54%;CDOM的平均貢獻(xiàn)率10.27%,略大于浮游植物的吸收。由此可以看出,影響遼河春、秋季水體吸收特性的主要因素是非藻類顆粒物,而浮游植物和CDOM的影響相對較小,主要原因是遼河春季水體中葉綠素濃度相對較低,非藻類顆粒物是總懸浮物的主要組成部分;而由于7—9月降雨量的急劇增加,使地表徑流增加,從而使得非藻類顆粒物的含量較高,其吸收貢獻(xiàn)率較5月有所增加。

    3結(jié)論

    本文通過遼河水體不同月份固有光學(xué)特性的分析,得到如下結(jié)論:

    (1)遼河總顆粒物的吸收光譜曲線基本與浮游植物吸收光譜一致,并且具有很強(qiáng)的空間差異性;總顆粒物的吸收與非藻類顆粒物的吸收以及浮游植物的吸收都呈現(xiàn)顯著相關(guān)。兩個(gè)季節(jié),非藻類顆粒物在440nm的吸收是總吸收的主要組成部分,而aCDOM(440)則大于aph(440)的吸收。9月非藻類顆粒物的吸收所占比例較5月高。

    (2)通過非藻類顆粒物吸收斜率的對比分析發(fā)現(xiàn)遼河水體非藻類顆粒物中礦物顆粒物的含量較其他一般水體要高。

    (3)兩季節(jié)的aph(440)/aph(675)的變化范圍比較相近(5月均值為2.08m-1,9月為2.56m-1),說明浮游植物輔助色素與葉綠素a濃度的組成季節(jié)性變化不大;而同一季節(jié)隨采樣點(diǎn)變化較大,說明其空間特征表現(xiàn)較為明顯。

    (4)遼河CDOM吸收斜率(SCDOM)基本小于其他水體(河流、湖泊、河口),說明其CDOM的組成主要是大分子物質(zhì),與河流攜帶較多的陸源DOM關(guān)系密切。

    致謝:感謝中國科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所水環(huán)境遙感學(xué)科組碩士研究生關(guān)瑩、劉蕾、李思佳在數(shù)據(jù)采集與室內(nèi)分析中給予的幫助;感謝馬建行師弟在數(shù)據(jù)分析中給予的幫助和支持。

    參考文獻(xiàn)(References):

    [1]MaRH,TangJW,DaiJF.Bio-opticalmodelwithoptimalparametersuitableforTaihuLakeinwatercolourremotesensing.InternationalJournalofRemoteSensing, 2006, 27(19): 4305- 4328.

    [2]Dall′OlmoG,GitelsonA.Effectofbio-opticalparametervariabilityontheremoteestimationofchlorophyll-aconcentrationinturbidproductivewaters:experimentalresults.AppliedOptics, 2005, 44(3): 412- 422.

    [3]孫德勇, 李云梅, 王橋, 樂成峰, 黃昌春, 王利珍. 內(nèi)陸湖泊水體固有光學(xué)特性的典型季節(jié)差異. 應(yīng)用生態(tài)學(xué)報(bào), 2008, 19(5): 1117- 1124.

    [4]KirkJTO.LightandPhotosynthesisinAquaticEcosystems.Britain:CambridgeUniversityPress, 1994.

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

    [6]IOCCG.IOCCGReportNumber5:RemoteSensingofInherentOpticalProperties:Fundamentals,TestsofAlgorithms,andApplications.Dartmouth,NS,Canada:IOCCG, 2006.

    [7]雷霞, 郭子祺, 田野, 謝飛, 秦靜欣. 官廳水庫秋季懸浮顆粒物和CDOM吸收特性. 湖泊科學(xué), 2013, 25(6): 883- 891.

    [8]金鑫, 李云梅, 王橋, 劉忠華, 王彥飛, 張紅, 尹斌, 徐祎凡, 徐昕. 巢湖水體漫衰減系數(shù)特性及其影響因子分析. 環(huán)境科學(xué), 2011, 32(6): 1581- 1590.

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

    [10]李方, 徐京萍, 何艷芬, 張柏, 劉殿偉, 宋開山, 王宗明, 章光新. 長春市石頭口門水庫顆粒物光譜吸收特性. 湖泊科學(xué), 2009, 21(2): 280- 287.

    [11]王明珠, 張運(yùn)林, 施坤, 高贇, 劉剛, 蔣浩. 千島湖水體懸浮顆粒物吸收特性及其典型季節(jié)差異. 環(huán)境科學(xué), 2014, 35(7): 2528- 2538.

    [12]ChenH,MengW,ZhengBH,WangCY,AnLH.Opticalsignaturesofdissolvedorganicmatterinthewatershedofagloballylargeriver(YangtzeRiver,China).Limnologica-EcologyandManagementofInlandWaters, 2013, 43(6): 482- 491.

    [13]惠秀娟, 楊濤, 李法云, 胡成, 王彤, 徐成斌, 侯偉, 馬放. 遼寧省遼河水生態(tài)系統(tǒng)健康評價(jià). 應(yīng)用生態(tài)學(xué)報(bào), 2011, 22(1): 181- 188.

    [14]王雪蕾, 蔡明勇, 鐘部卿, 姚延娟, 殷守敬, 吳迪. 遼河流域非點(diǎn)源污染空間特征遙感解析. 環(huán)境科學(xué), 2013, 34(10): 3788- 3796.

    [15]MitchellBG.Algorithmsfordeterminingtheabsorptioncoefficientforaquaticparticulatesusingthequantitativefiltertechnique//Proc.SPIE1302,OceanOpticsX, 137.Orlando,FL,UnitedStates, 1990: 137- 148.

    [16]BricaudA,MorelA,PrieurL.Absorptionbydissolvedorganicmatterofthesea(yellowsubstance)intheUVandvisibledomains.LimnologyandOceanography, 1981, 26(1): 43- 53.

    [17]SongKS,WangZM,BlackwellJ,ZhangB,LiF,ZhangYZ,JiangGJ.WaterqualitymonitoringusingLandsatThemateMapperdatawithempiricalalgorithmsinChaganLake,China.JournalofAppliedRemoteSensing, 2011, 5(1): 053506.

    [18]HelmsJR,StubbinsA,RitchieJD,MinorEC,KieberDJ,MopperK.Absorptionspectralslopesandsloperatiosasindicatorsofmolecularweight,source,andphotobleachingofchromophoricdissolvedorganicmatter.LimnologyandOceanography, 2008, 53(3): 955- 969.

    [19]BloughNV,GreenSA.Spectroscopiccharacterizationandremotesensingofnon-livingorganicmatter//ZeppRG,SonntagC,eds.TheRoleofNonlivingOrganicMatterintheEarth′sCarbonCycle.Chichester:JohnWiley&Sons, 1995: 23- 45.

    [20]BindingCE,JohnHJ,BukataRP,BootyWG.SpectralabsorptionpropertiesofdissolvedandparticulatematterinLakeErie.RemoteSensingofEnvironment, 2008, 112(4): 1702- 1711.

    [21]AstorecaR,RousseauV,LancelotC.Coloureddissolvedorganicmatter(CDOM)inSouthernNorthSeawaters:Opticalcharacterizationandpossibleorigin.Estuarine,CoastalandShelfScience, 2009, 85(4): 633- 640.

    [22]VaehaetaloAV,WetzelRG,PaerlHW.LightabsorptionbyphytoplanktonandchromophoricdissolvedorganicmatterinthedrainagebasinandestuaryoftheNeuseRiver,NorthCarolina(U.S.A).FreshwaterBiology, 2005, 50(3): 477- 493.

    [23]曹文熙, 楊躍忠, 許曉強(qiáng), 黃良民, 張建林. 珠江口懸浮顆粒物的吸收光譜及其區(qū)域模式. 科學(xué)通報(bào), 2003, 48(17): 1876- 1882.

    [24]戴永寧, 李素菊, 王學(xué)軍. 巢湖水體固有光學(xué)特性研究. 環(huán)境科學(xué)研究, 2008, 21(5): 173- 177.

    [25]蓋利亞, 劉正軍, 張繼賢. 三峽壩區(qū)水體吸收系數(shù)的特征研究. 遙感學(xué)報(bào), 2010, 14(2): 313- 332.

    [26]張紅, 黃家柱, 李云梅, 呂恒, 吳傳慶, 金鑫, 王彥飛, 尹斌. 滇池水體光學(xué)物質(zhì)的吸收特性研究. 環(huán)境科學(xué), 2011, 32(2): 452- 463.

    [27]BabinM,StramskiD,FerrariGM,ClaustreH,BricaudA,ObolenskyG,HoepffnerN.Variationsinthelightabsorptioncoefficientsofphytoplankton,nonalgalparticles,anddissolvedorganicmatterincoastalwatersaroundEurope.JournalofGeophysicalResearch:Oceans(1978—2012), 2003, 108(C7),doi: 10.1029/2001JC000882.

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

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

    [30]SobekS,TranvikLJ,PrairieYT,KortelainenP,ColeJJ.Patternsandregulationofdissolvedorganiccarbon:Ananalysisof7,500widelydistributedlakes.LimnologyandOceanography, 2007, 52(3): 1208- 1219.

    [31]ChenZQ,LiY,PanJM.DistributionsofcoloreddissolvedorganicmatteranddissolvedorganiccarboninthePearlRiverEstuary,China.ContinentalShelfResearch, 2004, 24(16): 1845- 1856.

    [32]SpencerRGM,ButlerKD,AikenGR.DissolvedorganiccarbonandchromophoricdissolvedorganicmatterpropertiesofriversintheUSA.JournalofGeophysicalResearch, 2012, 117(G3):G03001,doi: 10.1029/2011JG001928.

    [33]ZhangYL,QinBQ,ZhangL,ZhuGW,ChenWM.SpectralabsorptionandfluorescenceofchromophoricdissolvedorganicmatterinshallowlakesinthemiddleandlowerreachesoftheYangtzeRiver.JournalofFreshwaterEcology, 2005, 20(3): 451- 459.

    [34]PopeRM,FryES.Absorptionspectrum(380—700nm)ofpurewater. (Ⅱ).Integratingcavitymeasurements.AppliedOptics, 1997, 36(33): 8710- 8723.

    AbsorptioncharacteristicsandseasonalvariationsofopticallyactivewaterconstituentsfromLiaoheRiver

    SHAOTiantian1, 2,SONGKaishan1,*,DINGZhi3,ZHAOYing1,WANGMing1,DUJia1,ZHANGBai1

    1 Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun 130102, China 2 University of Chinese Academy of Sciences, Beijing 100049, China 3 Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China

    KeyWords:LiaoheRiver;absorptioncharacteristics;particulates;CDOM;seasonalvariations

    Abstract:Theabsorptioncharacteristicsofopticallyactiveconstituents(OACs)inwatercolumnareimportantopticalpropertiesandkeyparametersforbio-opticalmodeling.Inherentopticalproperties(IOPs)ofOACsaremainlydeterminedbythreecomponents,i.e.,thespectralabsorptionofnon-algalparticles,phytoplanktonandchromophoricdissolvedorganicmatter(CDOM).Spectralabsorptionproperties,seasonalandspatialvariabilityofparticulates(totalsuspendedmatters(TSM),includingphytoplanktonandnon-algalparticles)andCDOMabsorptioncharacteristicwereanalyzedfromwatersamplescollectedinMayandSeptember2013,respectively.CompositionsofCDOMandnon-algalparticleswerealsoexaminedatthesametime.Theresultsdemonstratedthatmostabsorptioncurvesofthetotalsuspendedparticulate[ap(λ)]weresimilartonon-algalparticlesabsorption[ad(λ)],showinggreatvariationsalongtheLiaoheRiverobservedatthesametime.Significantcorrelationswereobtainedbetweenap(λ)andad(λ),especiallythecorrelationsbetweenap(440)andad(440) (R=0.99bothinMayandSeptember).Meanwhile,strongrelationshipsbetweenap(λ)andtheabsorptionofphytoplankton[aph(λ)]wereexhibitedat675nm,withcorrelationcoefficient(R)=0.81inMayandR=0.55inSeptember,respectively.Moreover,closerelationshipsbetweenapandTSMatboth440and675nmwererevealedwhilesignificantcorrelationswereobtainedbetweenap(675)andchlorophylla(Chla)forsamplescollectedinbothseasons.Accordingtoourinvestigation,nosignificantrelationshipswereexhibitedbetweenap(λ)anddissolvedorganiccarbon(DOC).Acloserelationshipsbetweenap(λ)andinorganicsuspendedparticulate(ISM)andorganicsuspendedparticulate(OSM)wererecordedinSeptember,whileitwasnotthecaseforsamplesinMay.SignificantcorrelationswerealsoyieldedbetweenCDOMandDOCconcentrationforsamplescollectedinSeptember(R=0.90)at440nm,whichwashigherthanthatinMay(R=0.55).Moreover,contentofmineralparticulatesinLiaoheRiverwashigherthanmostotherwaterbodiesbasedontheanalysisofslopesofnon-algalparticlesabsorptions.ThecompositionofauxiliarypigmentsandChlaconcentrationinphytoplanktonshowedsmallseasonalvariationsbutgreatspatialvariabilityasindicatedbytheratioofaph(440)/aph(675).Furthermore,absorptioncharacteristicsofCDOMexhibitedsimilarcurveswithnon-algalparticles,howevershowedgreatvariabilityamongdifferentsamplingsites.Inaddition,CDOMabsorptionslopes(SCDOM)oftheLiaoheRiver,thoughexhibitinggreatlyspatialvariations,weresmallerthanotherwaterbodies,whichindicatedthattheCDOMcompositionofLiaoheRiverwasincliningtohavehighmolecularweight.Ourinvestigationindicatedthatthetotalparticulateabsorptionwasdominatedbythenon-algalparticlesat440nmbothinMayandSeptember.Theaveragecontributionrateofnon-algalparticlesis67.46%inMay,whilefor97.5%ofwatersamples,thecontributionrateofnon-algalparticlesreachesupto50%inSeptember.Consequently,contributionrateofnon-algalparticlesat440nminSeptemberwashigherthanthatinMayforthehighprecipitationwhilecontributionrateofaCDOM(440)waslitterhigherthanthatofaph(440)inbothmonths.

    基金項(xiàng)目:國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目(2013CB430401); 國家自然科學(xué)基金項(xiàng)目(41471290)

    收稿日期:2014- 09- 10; 網(wǎng)絡(luò)出版日期:2015- 04- 16

    *通訊作者

    Correspondingauthor.E-mail:songks@iga.ac.cn

    DOI:10.5846/stxb201409101787

    邵田田,宋開山,丁智,趙瑩,王銘,杜嘉,張柏.遼河水體光學(xué)吸收特性的季節(jié)變化.生態(tài)學(xué)報(bào),2016,36(7):1861- 1871.

    ShaoTT,SongKS,DingZ,ZhaoY,WangM,DuJ,ZhangB.AbsorptioncharacteristicsandseasonalvariationsofopticallyactivewaterconstituentsfromLiaoheRiver.ActaEcologicaSinica,2016,36(7):1861- 1871.

    猜你喜歡
    季節(jié)變化遼河顆粒物
    遼河口
    傾聽
    遼河(2020年8期)2020-09-02 14:25:25
    遼河文訊
    遼河(2018年5期)2018-11-15 03:39:26
    南平市細(xì)顆粒物潛在來源分析
    基于ELM算法的遼河沈陽段水污染評價(jià)
    金華市小型封閉水體氮污染季節(jié)變化規(guī)律及其原因分析
    綠色科技(2016年24期)2017-03-30 17:24:11
    基于GIS的東川地區(qū)植被指數(shù)季節(jié)變化研究
    商(2016年22期)2016-07-08 14:05:14
    季節(jié)變化對蛋雞養(yǎng)殖的危害及其避抗措施
    錯(cuò)流旋轉(zhuǎn)填料床脫除細(xì)顆粒物研究
    不同類型人工濕地凈化水質(zhì)季節(jié)變化分析
    日韩伦理黄色片| 免费av中文字幕在线| 男女之事视频高清在线观看 | 国产一区亚洲一区在线观看| 黄频高清免费视频| 国产精品香港三级国产av潘金莲 | 女人爽到高潮嗷嗷叫在线视频| 啦啦啦在线观看免费高清www| 街头女战士在线观看网站| 免费女性裸体啪啪无遮挡网站| 亚洲人成网站在线观看播放| 成人三级做爰电影| 日韩伦理黄色片| 在线观看免费午夜福利视频| 国产精品二区激情视频| 一级黄片播放器| 精品一区在线观看国产| 乱人伦中国视频| 男女高潮啪啪啪动态图| 天堂俺去俺来也www色官网| 一级毛片电影观看| 麻豆av在线久日| 亚洲,欧美,日韩| 欧美日韩亚洲高清精品| 色播在线永久视频| 国产精品人妻久久久影院| 中文欧美无线码| 天天添夜夜摸| 精品国产一区二区三区四区第35| 国产亚洲av高清不卡| 国语对白做爰xxxⅹ性视频网站| 精品久久蜜臀av无| 亚洲七黄色美女视频| 宅男免费午夜| 国产在线一区二区三区精| 国产精品一区二区精品视频观看| 亚洲色图综合在线观看| 亚洲国产欧美在线一区| 视频在线观看一区二区三区| 99精国产麻豆久久婷婷| 亚洲成人av在线免费| 日韩电影二区| 亚洲,欧美,日韩| 免费高清在线观看日韩| 在现免费观看毛片| 51午夜福利影视在线观看| 狠狠婷婷综合久久久久久88av| 国产乱人偷精品视频| 搡老岳熟女国产| 一区在线观看完整版| 国产伦理片在线播放av一区| 青青草视频在线视频观看| 伦理电影大哥的女人| 高清视频免费观看一区二区| 人人妻人人澡人人爽人人夜夜| 免费观看人在逋| 日韩一本色道免费dvd| 十分钟在线观看高清视频www| 欧美在线一区亚洲| 黄片无遮挡物在线观看| 国产激情久久老熟女| 亚洲精品av麻豆狂野| 高清欧美精品videossex| 久热爱精品视频在线9| 国产xxxxx性猛交| 精品酒店卫生间| 夫妻性生交免费视频一级片| 欧美日本中文国产一区发布| 精品久久久精品久久久| 欧美久久黑人一区二区| 亚洲一级一片aⅴ在线观看| 69精品国产乱码久久久| 欧美日韩视频精品一区| 欧美xxⅹ黑人| av有码第一页| 亚洲视频免费观看视频| 丝袜喷水一区| 飞空精品影院首页| 少妇人妻精品综合一区二区| 一级毛片我不卡| 日韩成人av中文字幕在线观看| 亚洲av福利一区| 99热国产这里只有精品6| 人成视频在线观看免费观看| 国产精品二区激情视频| 秋霞伦理黄片| 女人高潮潮喷娇喘18禁视频| 亚洲图色成人| 国产日韩欧美亚洲二区| 久久久国产欧美日韩av| 国产精品久久久av美女十八| 一区福利在线观看| 亚洲第一青青草原| 久久国产精品大桥未久av| 美女大奶头黄色视频| 电影成人av| 香蕉国产在线看| 国产成人91sexporn| 婷婷色麻豆天堂久久| 国产成人精品久久二区二区91 | 国产福利在线免费观看视频| 欧美精品一区二区免费开放| 一本色道久久久久久精品综合| 国产精品人妻久久久影院| 亚洲色图 男人天堂 中文字幕| 亚洲一区二区三区欧美精品| 中文字幕另类日韩欧美亚洲嫩草| 久久久久精品人妻al黑| 久久久精品免费免费高清| 国产精品麻豆人妻色哟哟久久| 三上悠亚av全集在线观看| 亚洲成人一二三区av| 亚洲欧美一区二区三区国产| 超碰成人久久| 王馨瑶露胸无遮挡在线观看| 青青草视频在线视频观看| 激情五月婷婷亚洲| 久久久久精品国产欧美久久久 | 日日爽夜夜爽网站| 精品第一国产精品| 国产成人啪精品午夜网站| 18禁裸乳无遮挡动漫免费视频| 色婷婷av一区二区三区视频| 69精品国产乱码久久久| 国产 一区精品| 老鸭窝网址在线观看| 男女免费视频国产| 亚洲av电影在线进入| 看免费成人av毛片| 国产极品天堂在线| 一级毛片黄色毛片免费观看视频| 精品人妻熟女毛片av久久网站| 宅男免费午夜| 国产av一区二区精品久久| avwww免费| 国产又爽黄色视频| 亚洲成人免费av在线播放| 国产亚洲精品第一综合不卡| 国产一级毛片在线| 久久人人爽人人片av| 狂野欧美激情性bbbbbb| 黄色视频在线播放观看不卡| 老司机影院成人| 成年人午夜在线观看视频| 欧美 亚洲 国产 日韩一| 9热在线视频观看99| 色婷婷久久久亚洲欧美| 老司机影院成人| 最近2019中文字幕mv第一页| 国产精品久久久久久精品电影小说| 亚洲欧美中文字幕日韩二区| 水蜜桃什么品种好| a 毛片基地| 久热爱精品视频在线9| 国产精品久久久久久精品电影小说| 日本欧美国产在线视频| 欧美国产精品va在线观看不卡| www日本在线高清视频| 色综合欧美亚洲国产小说| 最新在线观看一区二区三区 | 最近手机中文字幕大全| 久久久精品国产亚洲av高清涩受| 美女主播在线视频| 亚洲天堂av无毛| 国产国语露脸激情在线看| 亚洲美女视频黄频| 极品人妻少妇av视频| 七月丁香在线播放| 一级毛片 在线播放| 日韩大片免费观看网站| 爱豆传媒免费全集在线观看| 大码成人一级视频| av网站在线播放免费| 熟女少妇亚洲综合色aaa.| 激情五月婷婷亚洲| www日本在线高清视频| 男女午夜视频在线观看| 国产老妇伦熟女老妇高清| 激情五月婷婷亚洲| 欧美精品av麻豆av| 日韩免费高清中文字幕av| 久久精品久久久久久久性| 99久久人妻综合| 免费高清在线观看日韩| 少妇人妻久久综合中文| 亚洲婷婷狠狠爱综合网| 高清在线视频一区二区三区| 欧美日韩综合久久久久久| 久久久久视频综合| 国产精品久久久久久精品电影小说| 免费看不卡的av| 欧美精品一区二区免费开放| 啦啦啦中文免费视频观看日本| 99久久人妻综合| 欧美激情高清一区二区三区 | 久久人人爽人人片av| 精品少妇黑人巨大在线播放| 国产极品粉嫩免费观看在线| 欧美日韩国产mv在线观看视频| 91成人精品电影| 99国产综合亚洲精品| 观看av在线不卡| 一级毛片电影观看| 日韩 亚洲 欧美在线| 少妇精品久久久久久久| 99精国产麻豆久久婷婷| 可以免费在线观看a视频的电影网站 | 制服人妻中文乱码| 亚洲精品视频女| 国产精品 国内视频| 午夜福利网站1000一区二区三区| 亚洲精品自拍成人| 亚洲一卡2卡3卡4卡5卡精品中文| 日本vs欧美在线观看视频| xxx大片免费视频| 久久久久精品国产欧美久久久 | 久久久久人妻精品一区果冻| 日本av手机在线免费观看| 午夜91福利影院| 免费黄色在线免费观看| 久久久久久久国产电影| 久久人人97超碰香蕉20202| tube8黄色片| 午夜免费男女啪啪视频观看| 久久久久精品人妻al黑| 天天添夜夜摸| 亚洲av欧美aⅴ国产| 亚洲熟女精品中文字幕| 中国三级夫妇交换| www.熟女人妻精品国产| 亚洲色图综合在线观看| 欧美人与性动交α欧美软件| 老司机影院毛片| 叶爱在线成人免费视频播放| 美女大奶头黄色视频| 69精品国产乱码久久久| 亚洲免费av在线视频| 丰满乱子伦码专区| 亚洲情色 制服丝袜| 高清欧美精品videossex| 日韩中文字幕欧美一区二区 | 精品卡一卡二卡四卡免费| 丝袜脚勾引网站| 女人被躁到高潮嗷嗷叫费观| 99久久人妻综合| av有码第一页| 日韩一卡2卡3卡4卡2021年| 国语对白做爰xxxⅹ性视频网站| 久久久久久人妻| 巨乳人妻的诱惑在线观看| 国产一卡二卡三卡精品 | 一级毛片 在线播放| 国产精品久久久久成人av| 美女中出高潮动态图| 久久天躁狠狠躁夜夜2o2o | 少妇猛男粗大的猛烈进出视频| 国产乱人偷精品视频| 老司机影院毛片| 99久久精品国产亚洲精品| 午夜久久久在线观看| 天美传媒精品一区二区| 国语对白做爰xxxⅹ性视频网站| 9热在线视频观看99| 黄色 视频免费看| 国产精品久久久av美女十八| 黄片小视频在线播放| 免费高清在线观看日韩| 日韩人妻精品一区2区三区| 婷婷色综合大香蕉| 国产男人的电影天堂91| 我的亚洲天堂| 欧美 亚洲 国产 日韩一| 一区二区三区激情视频| 伊人久久国产一区二区| 国产午夜精品一二区理论片| 成人亚洲精品一区在线观看| av网站免费在线观看视频| 99香蕉大伊视频| 熟女少妇亚洲综合色aaa.| 久久97久久精品| 男人爽女人下面视频在线观看| 国产一级毛片在线| 国产一区有黄有色的免费视频| 深夜精品福利| 国产免费福利视频在线观看| 如日韩欧美国产精品一区二区三区| 一本色道久久久久久精品综合| 久久热在线av| 又黄又粗又硬又大视频| 亚洲精品国产区一区二| 国产xxxxx性猛交| 在线天堂中文资源库| 久久99一区二区三区| 啦啦啦在线观看免费高清www| 中文乱码字字幕精品一区二区三区| 哪个播放器可以免费观看大片| 狂野欧美激情性xxxx| 91精品国产国语对白视频| 中文字幕亚洲精品专区| 在线天堂最新版资源| 免费不卡黄色视频| 高清欧美精品videossex| 成人手机av| 人人妻人人澡人人看| www.精华液| netflix在线观看网站| 麻豆精品久久久久久蜜桃| 日本91视频免费播放| 十八禁网站网址无遮挡| 国产精品久久久久久人妻精品电影 | 激情视频va一区二区三区| 久久国产精品男人的天堂亚洲| 黄色一级大片看看| 午夜日本视频在线| 桃花免费在线播放| 久久精品久久精品一区二区三区| 久久精品亚洲熟妇少妇任你| 欧美日韩国产mv在线观看视频| 欧美成人午夜精品| 国产又爽黄色视频| 国产成人精品在线电影| 国产一区二区三区av在线| 国产成人啪精品午夜网站| 国语对白做爰xxxⅹ性视频网站| 黄色 视频免费看| 国产成人午夜福利电影在线观看| 国产伦人伦偷精品视频| 又大又黄又爽视频免费| 9色porny在线观看| 午夜激情av网站| 免费黄频网站在线观看国产| 国产一区二区激情短视频 | 最近2019中文字幕mv第一页| 午夜福利视频在线观看免费| 国产成人精品福利久久| 日韩av在线免费看完整版不卡| 成人国语在线视频| 99久久综合免费| 欧美日韩亚洲高清精品| 宅男免费午夜| 国产人伦9x9x在线观看| 亚洲欧美一区二区三区久久| 国产精品免费视频内射| 亚洲国产看品久久| 色网站视频免费| 久久99精品国语久久久| 欧美日韩福利视频一区二区| 国产精品三级大全| 色婷婷久久久亚洲欧美| 丁香六月天网| 亚洲人成77777在线视频| 波多野结衣一区麻豆| 成人18禁高潮啪啪吃奶动态图| 水蜜桃什么品种好| 国产亚洲av片在线观看秒播厂| www.av在线官网国产| 中国三级夫妇交换| 在线 av 中文字幕| 欧美人与善性xxx| bbb黄色大片| 亚洲国产精品一区三区| 久久久国产欧美日韩av| 亚洲成人av在线免费| 看免费av毛片| 1024香蕉在线观看| 色综合欧美亚洲国产小说| 男人舔女人的私密视频| 亚洲第一青青草原| 99热国产这里只有精品6| 久久久久网色| 两个人免费观看高清视频| 街头女战士在线观看网站| 亚洲欧美精品综合一区二区三区| av.在线天堂| 久久久久国产一级毛片高清牌| 热99国产精品久久久久久7| 美女中出高潮动态图| 亚洲国产欧美一区二区综合| 久久久久精品国产欧美久久久 | 亚洲,欧美精品.| 操美女的视频在线观看| 亚洲三区欧美一区| 久久精品国产亚洲av高清一级| 亚洲国产欧美在线一区| 久久久久人妻精品一区果冻| 午夜老司机福利片| 久久国产亚洲av麻豆专区| 久久精品aⅴ一区二区三区四区| av福利片在线| 不卡av一区二区三区| 成年人免费黄色播放视频| www.精华液| 日本91视频免费播放| 欧美日韩成人在线一区二区| 亚洲欧美清纯卡通| 成年人免费黄色播放视频| 男女高潮啪啪啪动态图| 在线观看人妻少妇| 美女主播在线视频| 在线 av 中文字幕| 啦啦啦 在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 老汉色∧v一级毛片| 在线天堂最新版资源| 国产日韩一区二区三区精品不卡| 又粗又硬又长又爽又黄的视频| 黄色毛片三级朝国网站| 新久久久久国产一级毛片| 亚洲一码二码三码区别大吗| 精品一区在线观看国产| 国产精品久久久久久精品古装| 亚洲av日韩精品久久久久久密 | a级片在线免费高清观看视频| 国产麻豆69| 秋霞在线观看毛片| 人人妻人人澡人人爽人人夜夜| 激情五月婷婷亚洲| videos熟女内射| 夫妻性生交免费视频一级片| 免费观看av网站的网址| 国产99久久九九免费精品| 成人午夜精彩视频在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲一区二区三区欧美精品| 看免费av毛片| 午夜福利视频在线观看免费| 亚洲精品视频女| 日日爽夜夜爽网站| 欧美精品高潮呻吟av久久| 麻豆av在线久日| 一二三四中文在线观看免费高清| 精品亚洲成国产av| 午夜日本视频在线| 少妇猛男粗大的猛烈进出视频| 日韩视频在线欧美| 91国产中文字幕| 午夜91福利影院| 久久久精品免费免费高清| 别揉我奶头~嗯~啊~动态视频 | 午夜福利,免费看| 久久这里只有精品19| 男人爽女人下面视频在线观看| 欧美黄色片欧美黄色片| 亚洲av成人精品一二三区| 亚洲精品日本国产第一区| 亚洲成av片中文字幕在线观看| 国产一区二区在线观看av| 久久久久精品性色| 国产欧美亚洲国产| 在线观看三级黄色| 午夜免费观看性视频| 国产xxxxx性猛交| 国产黄频视频在线观看| 侵犯人妻中文字幕一二三四区| av有码第一页| 久久久久国产一级毛片高清牌| 国产片特级美女逼逼视频| 国产伦理片在线播放av一区| 精品国产一区二区久久| 最近手机中文字幕大全| 日韩欧美一区视频在线观看| 香蕉丝袜av| 国产高清国产精品国产三级| 欧美成人精品欧美一级黄| 久久人妻熟女aⅴ| 人成视频在线观看免费观看| 国产视频首页在线观看| 亚洲欧美色中文字幕在线| 亚洲免费av在线视频| 亚洲精品一区蜜桃| bbb黄色大片| 亚洲熟女毛片儿| 这个男人来自地球电影免费观看 | 1024香蕉在线观看| 最新在线观看一区二区三区 | 高清欧美精品videossex| 国产成人精品福利久久| 两个人看的免费小视频| 亚洲国产精品999| 2018国产大陆天天弄谢| 免费少妇av软件| 久久性视频一级片| 中国三级夫妇交换| 亚洲国产欧美日韩在线播放| 满18在线观看网站| av卡一久久| 久久久久久久久久久免费av| 夫妻午夜视频| 久久国产精品男人的天堂亚洲| 2018国产大陆天天弄谢| 久热这里只有精品99| 大片电影免费在线观看免费| 日韩欧美一区视频在线观看| 一级毛片 在线播放| 十八禁高潮呻吟视频| 日韩精品有码人妻一区| 最近2019中文字幕mv第一页| 老鸭窝网址在线观看| 久久久久精品性色| 最近最新中文字幕大全免费视频 | 丝袜美足系列| 一级毛片黄色毛片免费观看视频| 亚洲精品国产av成人精品| 国产精品.久久久| 亚洲色图综合在线观看| 美女高潮到喷水免费观看| 国产成人a∨麻豆精品| 国产一区二区三区综合在线观看| 最近中文字幕高清免费大全6| 熟女av电影| 亚洲情色 制服丝袜| 18禁国产床啪视频网站| 久久久久精品国产欧美久久久 | 最近最新中文字幕大全免费视频 | 亚洲 欧美一区二区三区| 成年动漫av网址| 久久久久久人妻| 纵有疾风起免费观看全集完整版| 91精品伊人久久大香线蕉| 亚洲久久久国产精品| 精品一区在线观看国产| 国产精品久久久人人做人人爽| 亚洲国产中文字幕在线视频| 欧美av亚洲av综合av国产av | 久久99精品国语久久久| 精品一区二区三区av网在线观看 | 丝袜在线中文字幕| 亚洲色图 男人天堂 中文字幕| 亚洲一区中文字幕在线| 免费观看av网站的网址| 久久韩国三级中文字幕| avwww免费| 男女边吃奶边做爰视频| 欧美黑人精品巨大| 多毛熟女@视频| 秋霞伦理黄片| 久久久久视频综合| 男女无遮挡免费网站观看| 日本wwww免费看| 欧美精品高潮呻吟av久久| 菩萨蛮人人尽说江南好唐韦庄| 老熟女久久久| 午夜av观看不卡| 国产 精品1| 美女脱内裤让男人舔精品视频| 亚洲精品久久午夜乱码| 精品人妻在线不人妻| 可以免费在线观看a视频的电影网站 | 999久久久国产精品视频| 亚洲精品久久久久久婷婷小说| 一本大道久久a久久精品| 看十八女毛片水多多多| 夫妻午夜视频| 亚洲第一av免费看| 日韩人妻精品一区2区三区| 亚洲国产精品一区三区| 精品一品国产午夜福利视频| 国产成人91sexporn| 18禁观看日本| 亚洲人成电影观看| 九色亚洲精品在线播放| 最近最新中文字幕大全免费视频 | 免费观看性生交大片5| 18禁国产床啪视频网站| 自拍欧美九色日韩亚洲蝌蚪91| 国产乱人偷精品视频| 亚洲天堂av无毛| 久久精品久久久久久久性| 美女大奶头黄色视频| www日本在线高清视频| 天天添夜夜摸| 日日爽夜夜爽网站| 人妻一区二区av| 亚洲av电影在线观看一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲一级一片aⅴ在线观看| 女人被躁到高潮嗷嗷叫费观| 香蕉国产在线看| 狂野欧美激情性bbbbbb| av网站免费在线观看视频| 天天躁日日躁夜夜躁夜夜| 国产男女超爽视频在线观看| 久久热在线av| 亚洲国产精品999| 日韩av免费高清视频| 精品国产一区二区三区四区第35| 国产97色在线日韩免费| 十八禁高潮呻吟视频| 亚洲国产av影院在线观看| 欧美日韩亚洲高清精品| 大片免费播放器 马上看| 国产在线免费精品| 中文字幕av电影在线播放| 国产欧美亚洲国产| 青青草视频在线视频观看| 亚洲,欧美,日韩| 久久影院123| 免费在线观看完整版高清| 亚洲国产精品一区三区| 国产精品国产av在线观看| √禁漫天堂资源中文www| 男女无遮挡免费网站观看| 国产精品免费大片| 日本欧美视频一区| 啦啦啦在线观看免费高清www| 少妇精品久久久久久久| av又黄又爽大尺度在线免费看| 国产野战对白在线观看| 久久久久视频综合| 欧美精品高潮呻吟av久久| 国产男人的电影天堂91| 国产又爽黄色视频| 午夜久久久在线观看| h视频一区二区三区| 伦理电影免费视频|