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

    福建敖江口表層沉積硅藻空間分布特征

    2021-12-14 03:29:10張紀(jì)暉周成旭李冬玲林忠洲邱悅沙龍濱
    關(guān)鍵詞:河口區(qū)環(huán)境變量江口

    張紀(jì)暉, 周成旭, 李冬玲, 林忠洲, 邱悅, 沙龍濱*

    福建敖江口表層沉積硅藻空間分布特征

    張紀(jì)暉1a, 周成旭1b, 李冬玲1a, 林忠洲2, 邱悅1a, 沙龍濱1a*

    (1. 寧波大學(xué), a. 地理與空間信息技術(shù)系; b. 食品與藥學(xué)學(xué)院, 浙江 寧波 315832; 2. 寧波海洋研究院, 浙江 寧波 315000)

    為探討福建敖江口表層沉積硅藻空間分布特征,2019年7月對(duì)敖江口進(jìn)行表層沉積硅藻的采樣調(diào)查,并研究了硅藻與環(huán)境因子的關(guān)系。結(jié)果表明,從13個(gè)站點(diǎn)共檢出硅藻114種,隸屬于39屬。硅藻豐度具有從河口向外海先減少后增加的變化特征。聚類分析表明,主要硅藻屬種可劃分為淡水種組合帶和半咸水-咸水種組合帶。典型對(duì)應(yīng)分析(CCA)表明,表層海水溫度和鹽度是影響硅藻分布的主要環(huán)境變量。河口區(qū)鹽度相對(duì)較低,主要分布淡水種雙面曲殼藻、優(yōu)美曲殼藻和顆粒溝鏈藻;外海區(qū)北部鹽度相對(duì)較高,主要分布咸水種流水雙菱藻;外海區(qū)南部鹽度相對(duì)較低,主要分布淡水種優(yōu)美曲殼藻和顆粒溝鏈藻;灘涂的沉積物較粗,主要分布優(yōu)美曲殼藻。此外,河口區(qū)和外海區(qū)南部可能存在一定的水體污染。

    福建;敖江;硅藻;環(huán)境;典型對(duì)應(yīng)分析;聚類分析

    近岸海域是航道交通、生產(chǎn)生活的重點(diǎn)區(qū)域。隨著社會(huì)經(jīng)濟(jì)的發(fā)展和海岸海洋資源開發(fā)力度的不斷加大,人類活動(dòng)對(duì)近岸海域的生態(tài)環(huán)境產(chǎn)生著越來越重要的影響:不僅包括海水養(yǎng)殖,海洋捕撈等社會(huì)經(jīng)濟(jì)活動(dòng),還包括伴隨這些活動(dòng)而產(chǎn)生的各類環(huán)境污染。為研究日益復(fù)雜化的近岸海域水體環(huán)境,在研究中往往需要尋找一種對(duì)環(huán)境變量敏感的物種,通過研究該物種的時(shí)空分布變化特征,從而指示海域的環(huán)境變化。

    硅藻是一種個(gè)體微小的單細(xì)胞生物,分布范圍極其廣泛,在陸地湖泊、近岸海域以及海洋中均有發(fā)現(xiàn)[1]。硅藻對(duì)環(huán)境因子(如溫度、鹽度、酸堿度、營(yíng)養(yǎng)鹽含量等)的變化靈敏[2],能對(duì)水體環(huán)境變化做出迅速的反應(yīng),因而其組合結(jié)構(gòu)就是對(duì)當(dāng)前水體環(huán)境特點(diǎn)最直接的體現(xiàn)[3],硅藻已成為廣泛用于水體環(huán)境監(jiān)測(cè)的指示類物種[4–5]。

    敖江口是福建敖江干流的入???,位于福建連江縣,坐落于我國(guó)浙閩東南沿海地區(qū),是我國(guó)海洋經(jīng)濟(jì)最為發(fā)達(dá)的區(qū)域之一。連江縣生產(chǎn)生活主要以第一產(chǎn)業(yè)為主,以海洋經(jīng)濟(jì)為傳統(tǒng)支柱產(chǎn)業(yè),海水養(yǎng)殖業(yè)尤為發(fā)達(dá),以養(yǎng)殖鮑魚、縊蟶和海帶等水產(chǎn)為名。發(fā)達(dá)的海水養(yǎng)殖業(yè)勢(shì)必會(huì)對(duì)敖江口水體環(huán)境產(chǎn)生影響,然而目前使用硅藻作為指標(biāo)對(duì)該區(qū)域水體環(huán)境進(jìn)行研究較少,目前已有采用其他指標(biāo)對(duì)敖江口周邊和敖江流域水體環(huán)境狀況進(jìn)行研究的報(bào)道,阮金山等[6]對(duì)連江縣東南部養(yǎng)殖水體中的重金屬和有機(jī)污染物含量進(jìn)行了研究;張玉珍等[7]對(duì)敖江全流域底棲動(dòng)物和魚類群落結(jié)構(gòu)進(jìn)行了研究;張丹丹等[8]對(duì)敖江下游河段的抗生素抗性基因分布特征進(jìn)行了研究。

    本文通過在福建敖江口采集表層沉積物,調(diào)查河口、灘涂和外海區(qū)表層沉積硅藻分布情況,探討了敖江口表層沉積硅藻的空間分布特征,并對(duì)主要硅藻屬種與表層海水溫度、鹽度、pH等環(huán)境因子間的關(guān)系進(jìn)行分析,了解該海域河口、灘涂和外海區(qū)影響硅藻分布的主要環(huán)境因素。研究結(jié)果將為了解敖江口硅藻對(duì)環(huán)境的響應(yīng)機(jī)制、該區(qū)域生態(tài)系統(tǒng)的初級(jí)生產(chǎn)力變化,以及人類活動(dòng)對(duì)該區(qū)域生態(tài)環(huán)境的綜合影響研究提供依據(jù)。

    1 材料和方法

    1.1 樣品和環(huán)境變量的采集

    表層沉積物采樣工作于2019年7月完成。使用抓斗在13個(gè)采樣點(diǎn)抓取底部表層沉積物,樣品密封保存送至實(shí)驗(yàn)室。使用HORIBA設(shè)備收集采樣點(diǎn)表層海水溫度、鹽度、pH、濁度等環(huán)境信息。采樣站點(diǎn)位置見圖1,其中外海站點(diǎn)4個(gè)(WH1、WH2、WH3、WH4);灘涂站點(diǎn)6個(gè)(TT1、TT2、TT3、TT4、TT5、TT6);河口站點(diǎn)3個(gè)(HD1、HD2、HD4)。

    1.2 方法

    硅藻鑒定 稱取約1 g烘干樣品放入燒杯, 加入10%的稀鹽酸攪拌均勻,待充分反應(yīng)后加入30%的雙氧水去除有機(jī)質(zhì),待完全反應(yīng)并出現(xiàn)明顯分層后,加入去離子水靜置約10 h,抽去上層清液, 重復(fù)洗酸3次;將樣品倒入比色管中,加入明膠溶液和去離子水定容搖勻,轉(zhuǎn)移至固定好蓋玻片的培養(yǎng)皿中靜置24 h后引流;將處理好的蓋玻片用Naphrax膠(dn=1.73)固定于載玻片上。在徠卡光學(xué)顯微鏡下進(jìn)行鑒定,放大倍率為1 000×油鏡,隨機(jī)逐行統(tǒng)計(jì)硅藻屬種與數(shù)量,對(duì)于不完整的硅藻,若是超過一半則視為整體處理。參考Krammer等[9–12]和金德祥等[13]的方法進(jìn)行硅藻鑒定,大部分硅藻鑒定到種,少數(shù)不確定的鑒定到屬。

    圖1 敖江口采樣站點(diǎn)分布

    粒度分析 稱取烘干樣品0.2 g置于干凈燒杯中,先加入10 mL 10%的雙氧水溶液去除有機(jī)質(zhì),然后加入10 mL 10%的HCl溶液去除碳酸鹽;水浴加熱2 h后,燒杯中加滿去離子水靜置12 h。待樣品顆粒盡量分散后抽取上層清液,加入2滴1 mol/L的六偏磷酸鈉[(NaPO3)6]溶液,放入超聲波儀中震蕩分散15 min后在Beckman Coulter LS13320型激光粒度儀上進(jìn)行粒度測(cè)試,使用自帶軟件LS13320對(duì)數(shù)據(jù)進(jìn)行處理分析。依照???洛克三角分級(jí)和命名法分為3個(gè)粒度,即粘土(<4m)、粉砂(4~63m)和砂(>63m)。

    1.3 數(shù)據(jù)處理

    硅藻豐度計(jì)算 絕對(duì)豐度=/,式中,為絕對(duì)豐度(cell/g),為硅藻總個(gè)數(shù)(cell),為樣品干質(zhì)量(g)。在實(shí)際操作中,參照梁鈺瑩等[14]的方法,上式可轉(zhuǎn)化為=(×)/(×),式中,為顯微鏡視野內(nèi)硅藻個(gè)數(shù)(cell),為顯微鏡視野面積(cm2),為培養(yǎng)皿面積(cm2)。硅藻豐度的空間分布特征采用Surfer軟件分析。

    聚類分析 聚類分析是一種生態(tài)學(xué)研究的常用方法。以樣點(diǎn)間的相似性為依據(jù)劃分為不同的組,從而找出樣點(diǎn)間的共性與規(guī)律。根據(jù)樣點(diǎn)特征,當(dāng)樣點(diǎn)本身具有空間或時(shí)間的連續(xù)屬性時(shí),應(yīng)當(dāng)采用順序聚類的方式;反之當(dāng)樣點(diǎn)本身不具有空間或時(shí)間連續(xù)性時(shí),應(yīng)當(dāng)采用無序聚類的方式。本文采用Tilia軟件的無序聚類方式。

    相關(guān)性分析 硅藻分布不僅受生境的影響,還受到外在因素的制約,如局地海流、水體營(yíng)養(yǎng)化等。而通過相關(guān)性分析,如典型對(duì)應(yīng)分析(canonical correspondence analysis, CCA),可以將硅藻屬種與環(huán)境變量的相互關(guān)系通過二維圖的形式來表現(xiàn)出來,從而得出硅藻屬種與各環(huán)境變量間的關(guān)系[15]。CCA分析是一種有效的直接梯度分析法,將環(huán)境變量、物種數(shù)據(jù)和站點(diǎn)同時(shí)顯示在一個(gè)低維空間中,可以非常直觀地展現(xiàn)出環(huán)境變量與物種之間的關(guān)系特征。CCA分析的假設(shè)前提是物種與環(huán)境的關(guān)系呈現(xiàn)非線性的單峰響應(yīng)特點(diǎn)。因此,在進(jìn)行CCA分析之前,應(yīng)當(dāng)首先進(jìn)行去趨勢(shì)分析(detrended corre- spondence analysis, DCA),得到物種數(shù)據(jù)的梯度長(zhǎng)度, 當(dāng)梯度長(zhǎng)度>4時(shí),表明物種數(shù)據(jù)具有單峰分布的特點(diǎn),適用于CCA分析;若梯度長(zhǎng)度<3時(shí)物種數(shù)據(jù)更適用于冗余分析(redundancy analysis, RDA)等線性模型。本文采用CANOCO軟件進(jìn)行相關(guān)性分析。

    2 結(jié)果和分析

    2.1 環(huán)境變量特征

    從表1可知,研究區(qū)域內(nèi)表層海水溫度存在一定的波動(dòng),為26.00℃~28.07℃,最高值出現(xiàn)在TT1站點(diǎn),最低值出現(xiàn)在WH1站點(diǎn)。鹽度的變化很大(4.50‰~26.90‰),主要是由于敖江徑流輸入的影響,從而導(dǎo)致河口區(qū)的鹽度明顯低于外海和灘涂區(qū)。沉積物粒度存在顯著差異,河口區(qū)與外海區(qū)的沉積物多為粘土質(zhì)粉砂,而灘涂區(qū)的沉積物多為砂。研究區(qū)pH變化不大(7.86~8.26),溶解氧含量為5.86~7.47 mg/L,濁度為30.0~237.0 NTU,濁度的高值分別出現(xiàn)在HD1、HD2、HD4和TT6站點(diǎn)。

    2.2 主要硅藻屬種

    從沉積物中共鑒定出硅藻114種,分屬39屬,其中河口區(qū)站點(diǎn)的硅藻屬種最多,平均54種;灘涂區(qū)站點(diǎn)最少,平均22種。硅藻屬種最多的是HD4站點(diǎn),位于河口區(qū),共鑒定出59種;最少的是TT1站點(diǎn),位于灘涂區(qū),共鑒定出13種。

    敖江口表層沉積硅藻常見屬為小環(huán)藻屬()、曲殼藻屬()、輻環(huán)藻屬()、圓篩藻屬()和溝鏈藻屬();主要硅藻有雙面曲殼藻()、優(yōu)美曲殼藻()、諾氏輻環(huán)藻()、簡(jiǎn)單雙眉藻()、顆粒溝鏈藻()、條紋小環(huán)藻()、具槽帕拉藻()、流水雙菱藻()、菱形海線藻()和卵形摺盤藻()。

    表1 研究區(qū)站點(diǎn)的環(huán)境因子

    2.3 硅藻豐度特征

    將敖江口各站點(diǎn)的硅藻豐度導(dǎo)入Surfer軟件中,使用克里金插值法進(jìn)行插值計(jì)算,得到敖江口表層硅藻的空間分布情況(圖2)。硅藻在河口區(qū)和外海區(qū)分布相對(duì)集中,硅藻豐度由河口向外海呈現(xiàn)先減少后增加的特征。最大硅藻豐度位于HD2站點(diǎn),最小硅藻豐度位于TT2站點(diǎn)。硅藻在灘涂區(qū)豐度最小,平均豐度為5 136.67 cell/g;在河口區(qū)豐度最大,平均豐度為67 579.35 cell/g。

    2.4 硅藻分布特征和聚類分析

    采用Tilia軟件進(jìn)行聚類分析結(jié)果表明(圖3), 條紋小環(huán)藻、具槽直鏈藻、菱形海線藻和卵形摺盤藻在除TT1以外的所有站點(diǎn)均有分布,是研究區(qū)內(nèi)的優(yōu)勢(shì)種。根據(jù)聚類分析結(jié)果,可以將敖江口表層沉積硅藻劃分為2個(gè)組合帶:I帶有HD1、HD2和HD4共3個(gè)站點(diǎn),顆粒溝鏈藻和雙面曲殼藻在這3個(gè)站點(diǎn)含量較高,在其他站點(diǎn)幾乎沒有發(fā)現(xiàn);雖然優(yōu)美曲殼藻在這3個(gè)站點(diǎn)也有發(fā)現(xiàn),但占比相對(duì)較少,其主要在TT1站點(diǎn)有大量分布,因而I帶的主要硅藻組合為顆粒溝鏈藻-雙面曲殼藻。II帶有WH2、WH3、WH4、TT5和TT6共5個(gè)站點(diǎn),流水雙菱藻在這些站點(diǎn)的含量較高,在其他站點(diǎn)較少發(fā)現(xiàn);條紋小環(huán)藻在這些站點(diǎn)中的含量明顯高于I帶站點(diǎn),具槽直鏈藻、菱形海線藻和卵形摺盤藻的含量則變化不大,因而II帶主要硅藻組合為流水雙菱藻-條紋小環(huán)藻。

    圖2 敖江口硅藻豐度(cell/g)分布插值圖

    圖3 敖江口硅藻的聚類分析。A. bia: 雙面曲殼藻; A. del: 優(yōu)美曲殼藻; A. nor: 諾氏輻環(huán)藻; A. exi: 簡(jiǎn)單雙眉藻; A. gra: 顆粒溝鏈藻; C. sch: 叔曼小環(huán)藻; C. str: 條紋小環(huán)藻; E. gib: 駝峰窗紋藻; N. cus: 急尖舟形藻; P. sul: 具槽帕拉藻; S. flu: 流水雙菱藻; T. nit: 菱形海線藻; T. coc:卵形摺盤藻。圖5同。

    從圖3還可看出,WH1和TT1是特殊站點(diǎn), WH1站點(diǎn)中諾氏輻環(huán)藻的含量很高,優(yōu)美曲殼藻、顆粒溝鏈藻和叔曼小環(huán)藻()僅有少量,未檢出流水雙菱藻。相比于其他外海區(qū)站點(diǎn), WH1站點(diǎn)的硅藻分布特征更接近河口區(qū)站點(diǎn)。TT1站點(diǎn)的硅藻分布特征最為特殊,未見優(yōu)勢(shì)種(如條紋小環(huán)藻)以及I帶和II帶的主要硅藻(如顆粒溝鏈藻和流水雙菱藻),其含量最高的是優(yōu)美曲殼藻和簡(jiǎn)單雙眉藻,還有少量的駝峰窗紋藻()和急尖舟形藻()。

    2.5 相關(guān)性分析

    首先對(duì)敖江口表層沉積硅藻進(jìn)行DCA分析, 得到物種的梯度長(zhǎng)度為6.006,表明硅藻屬種在坐標(biāo)軸上具有非線性單峰分布的特點(diǎn),說明CCA分析更適用于敖江口表層沉積硅藻-環(huán)境變量的相關(guān)性分析。硅藻與環(huán)境變量的CCA分析結(jié)果表明,前兩軸特征值分別為0.837和0.567,可以解釋53.09%的硅藻-環(huán)境變量關(guān)系的變化特征,研究區(qū)硅藻屬種和環(huán)境變量之間具有很高的相關(guān)性。通過膨脹系數(shù)(VIF≥20)以及蒙特卡羅篩選檢驗(yàn)(999次循環(huán))后,最終篩選出5個(gè)環(huán)境變量:表層海水溫度、鹽度、濁度、pH和溶解氧含量。將這5個(gè)環(huán)境變量與硅藻屬種再次進(jìn)行CCA分析,結(jié)果表明前兩軸特征值分別為0.634和0.481,共解釋了70.16%的硅藻-環(huán)境變量關(guān)系變化特征(表2),比未篩選前有很大提高。此外,表層海水溫度與鹽度的特征值分別占5個(gè)環(huán)境變量的39.88%和30.28%,進(jìn)一步說明表層海水溫度與鹽度是影響敖江口表層沉積硅藻分布的關(guān)鍵環(huán)境變量。

    對(duì)站點(diǎn)-環(huán)境變量間的關(guān)系進(jìn)行CCA分析,環(huán)境變量用矢量表示,箭頭指向表示其極值,具有相似的硅藻屬種和環(huán)境變量的站點(diǎn)在圖中位置相近[15], 站點(diǎn)的投影在環(huán)境變量正延長(zhǎng)線上表示其與環(huán)境變量呈正相關(guān)性,反之呈負(fù)相關(guān)性[16]。從圖4可見, HD1、HD2、HD4和WH1站點(diǎn)與鹽度呈負(fù)相關(guān), WH3、WH4、TT5、TT6、WH2和TT1站點(diǎn)與鹽度呈正相關(guān);HD1、HD2、TT1和WH2站點(diǎn)與表層海水溫度呈正相關(guān),WH1、WH3、WH4、TT5、TT6和HD4站點(diǎn)與表層海水溫度呈負(fù)相關(guān)。從分析結(jié)果可知, 河口區(qū)水體鹽度相對(duì)較低,海水溫度相對(duì)較高;外海區(qū)海水溫度總體上相對(duì)較低,然而南部(WH1站點(diǎn))和北部(WH2、WH3和WH4站點(diǎn))的鹽度存在差異,表現(xiàn)為南部比北部低。灘涂區(qū)TT1站點(diǎn)的海水溫度和鹽度較高,而TT5和TT6站點(diǎn)的水體環(huán)境則更接近于外海區(qū)北部的WH3和WH4站點(diǎn)。

    表2 硅藻屬種-環(huán)境變量的CCA分析

    圖4 敖江口站點(diǎn)與環(huán)境變量的CCA關(guān)系圖。Turb: 濁度; Temp: 表層海水溫度; DO: 溶解氧; Sali: 鹽度。下圖同。

    對(duì)主要硅藻屬種-環(huán)境變量的關(guān)系進(jìn)行CCA分析,環(huán)境變量使用矢量表示,針對(duì)某一具體的環(huán)境變量,可以將該變量沿矢量反方向延長(zhǎng),從而將其近似地看作這一環(huán)境變量的“數(shù)軸”,硅藻屬種與環(huán)境變量的關(guān)系體現(xiàn)在硅藻投影到環(huán)境變量“數(shù)軸”上的得分情況。從圖5可見,大部分硅藻與鹽度呈正相關(guān),只有叔曼小環(huán)藻、諾氏輻環(huán)藻、顆粒溝鏈藻和雙面曲殼藻與鹽度呈負(fù)相關(guān)。在與表層海水溫度的相關(guān)性上,硅藻屬種的分布比較均衡。

    3 結(jié)論和討論

    河口與近岸海域是一個(gè)復(fù)雜的水體環(huán)境,作為河水、海水和人類活動(dòng)的交集地,往往會(huì)受到三者的共同影響, 因而,河口與近岸海域中的硅藻分布特征往往是溫、鹽、沉積環(huán)境和人類活動(dòng)等多重因素共同作用的結(jié)果[17]。福建敖江口是典型的潮間帶海域,通過對(duì)其表層沉積物中的硅藻進(jìn)行鑒定,發(fā)現(xiàn)條紋小環(huán)藻、具槽直鏈藻、菱形海線藻和卵形摺盤藻是該海域的優(yōu)勢(shì)種(圖3),其中條紋小環(huán)藻、具槽直鏈藻和菱形海線藻都是我國(guó)近岸海域最為常見的廣布種[18–21],卵形摺盤藻則是福建沿海的常見種[13]。結(jié)合CCA分析表明,由于敖江徑流和外海的共同影響,表層海水溫度和鹽度是影響敖江口表層沉積硅藻自河口向外海分布呈現(xiàn)空間差異的最重要環(huán)境變量。此外,從沉積物粒度、硅藻豐度及硅藻組成來看,河口、灘涂和外海沉積環(huán)境的差異和人類活動(dòng)因素,也對(duì)敖江口的表層沉積硅藻分布存在一定影響。

    圖5 敖江口硅藻與環(huán)境變量關(guān)系圖

    敖江徑流的輸入使鹽度成為硅藻分布最關(guān)鍵的影響因子。作為受敖江徑流影響最大的區(qū)域,河口區(qū)具有低鹽度、低表層海水溫度的特征,硅藻組合屬于I帶,主要分布有淡水種雙面曲殼藻、顆粒溝鏈藻和優(yōu)美曲殼藻[18,22–25]。此外,河口區(qū)也是敖江口硅藻豐度最高的區(qū)域(圖2)。河口區(qū)的濁度較高,推測(cè)可能是人類活動(dòng)引起的水體富營(yíng)養(yǎng)化。顆粒溝鏈藻對(duì)水體富營(yíng)養(yǎng)化程度具有指示作用,常在中度營(yíng)養(yǎng)鹽的水體中占據(jù)優(yōu)勢(shì)地位[23],該種在河口區(qū)較為豐富,指示了河口區(qū)的水體富營(yíng)養(yǎng)化較為嚴(yán)重。

    外海區(qū)遠(yuǎn)離敖江徑流,受人類活動(dòng)影響小,總體上海水溫度較低,然而在鹽度上卻存在空間差異,表現(xiàn)為外海區(qū)北部水體鹽度高于南部水體。聚類分析結(jié)果表明,外海區(qū)北部WH2、WH3和WH4站點(diǎn)硅藻組合屬于II帶,優(yōu)勢(shì)種為流水雙菱藻,是典型的咸水種[26–27],表明水體鹽度較高;但外海區(qū)南部水體中流水雙菱藻的含量很低,淡水種優(yōu)美曲殼藻和顆粒溝鏈藻的含量卻相對(duì)較高(圖3),推測(cè)這片區(qū)域可能受到了來自敖江南部的閩江淡水影響。此外,諾氏輻環(huán)藻在這片區(qū)域的含量較高(圖3),其被認(rèn)為是一種水體富營(yíng)養(yǎng)化指示種[23],因而外海區(qū)南部的水體中可能也存在一定的水體污染。

    河口區(qū)與外海區(qū)的主要沉積物都是粘土質(zhì)粉砂,而灘涂區(qū)的主要沉積物則是砂(表1)。這種沉積物粒度的空間差異性可能導(dǎo)致灘涂區(qū)硅藻分布特征有別于河口和外海區(qū)。聚類結(jié)果表明TT1站點(diǎn)的優(yōu)美曲殼藻和簡(jiǎn)單雙眉藻的含量非常高(圖3)。有研究表明, 優(yōu)美曲殼藻對(duì)沉積物粒度具有指示作用[19]。因此,在灘涂區(qū),沉積物粒度對(duì)硅藻分布的影響較為顯著。

    綜上所述,通過對(duì)福建敖江口表層沉積硅藻的采樣調(diào)查與結(jié)果分析,發(fā)現(xiàn)敖江口表層沉積硅藻在空間分布上具有顯著的差異性,表現(xiàn)為河口區(qū)、外海區(qū)北部、外海區(qū)南部和灘涂區(qū)4個(gè)區(qū)域主要硅藻屬種的差異性。相關(guān)性分析表明,表層沉積硅藻的空間分布受到鹽度和海水溫度的顯著影響。此外,研究區(qū)域中沉積物粒度和人為活動(dòng)因素對(duì)表層沉積硅藻的空間分布也具有一定的影響。

    [1] ZHANG J P. Study on sedimentary diatoms of late Holocene from middle offshore sea area of Fujian Province [D]. Xiamen: Xiamen University, 2009: 1–2.

    張金鵬. 福建中部近岸海域晚全新世沉積硅藻研究 [D]. 廈門: 廈門大學(xué), 2009: 1–2.

    [2] ANDERSON N J, VOS P. Learning from the past: Diatoms as palaeoecological indicators of changes in marine environments [J]. Neth J Aquat Ecol, 1992, 26(1): 19–30. doi: 10.1007/BF02298025.

    [3] QIU L C, WEI G F, LI X J, et al. Species diversity and temporal-spatial distribution of benthic diatoms in Jianjiang River, Guangdong Province [J]. J Trop Subtrop Bot, 2016, 24(2): 197–207. doi: 10.11926/j.issn. 1005-3395.2016.02.011.

    邱禮城, 韋桂峰, 李湘姣, 等. 廣東鑒江水系底棲硅藻多樣性與時(shí)空分布特征 [J]. 熱帶亞熱帶植物學(xué)報(bào), 2016, 24(2): 197–207. doi: 10.11926/j.issn.1005-3395.2016.02.011.

    [4] HE Q, WEI G F, LIU J, et al. Species diversity and spatial distribution of benthic diatom in Zeng River, Guangdong Province, southern China [J]. J Trop Subtrop Bot, 2011, 19(3): 245–253. doi: 10.3969/j.issn. 1005-3395.2011.03.009.

    何琦, 韋桂峰, 劉靜, 等. 增江底棲硅藻種類多樣性及分布特征 [J].熱帶亞熱帶植物學(xué)報(bào), 2011, 19(3): 245–253. doi: 10.3969/j.issn. 1005-3395.2011.03.009.

    [5] XU C C, YANG X Y, LU L, et al. Comparison of benthic diatom assemblages from mangrove and reed wetland in Qi’ao Island, Zhuhai, Guangdong Province, China [J]. J Trop Subtrop Bot, 2010, 18(3): 217– 223. doi: 10.3969/j.issn.1005-3395.2010.03.001.

    徐承超, 楊曉瑩, 盧亮, 等. 珠海淇澳島紅樹林和蘆葦濕地底棲硅藻群落比較 [J]. 熱帶亞熱帶植物學(xué)報(bào), 2010, 18(3): 217–223. doi: 10.3969/j.issn.1005-3395.2010.03.001.

    [6] RUAN J S, ZHONG S L, LIN H X, et al. Quality security assessment of shellfish-culture Lianjiang southeast sea area, Fujian [J]. J Fujian Fish, 2011, 33(2): 10–17.

    阮金山, 鐘碩良, 林后祥, 等. 福建連江縣東南部海域養(yǎng)殖貝類質(zhì)量安全評(píng)估 [J]. 福建水產(chǎn), 2011, 33(2): 10–17.

    [7] ZHANG Y Z, HUANG W D, WANG Z Y, et al. Evaluation of aquatic ecosystem health in Aojiang Basin, Fujian province [J]. J Lake Sci, 2015, 27(6): 1079–1086. doi: 10.18307/2015.0612.

    張玉珍, 黃文丹, 王智苑, 等. 福建敖江流域水域生態(tài)系統(tǒng)健康評(píng)估 [J]. 湖泊科學(xué), 2015, 27(6): 1079–1086. doi: 10.18307/2015.0612.

    [8] ZHANG D D, GUO Y P, REN H Y, et al. Characteristics of antibiotic resistance genes in downstream areas of the Aojiang River, Fujian Province [J]. Environ Sci, 2018, 39(6): 2600–2606. doi: 10.13227∕j. hjkx.201710121.

    張丹丹, 郭亞平, 任紅云, 等. 福建省敖江下游抗生素抗性基因分布特征 [J]. 環(huán)境科學(xué), 2018, 39(6): 2600–2606. doi: 10.13227∕j.hjkx. 201710121.

    [9] KRAMMER K, LANGE-BERTALOT H. Bacillariophyceae. teil 1: Naviculaceae [M]// ETTL H, G?RTNER G, GERLOFF J, et al. Süss- wasserflora von Mitteleuropa, 2(1). New York: Gustav Fischer Verlag, Stuttgart, 1986: 480–750.

    [10] KRAMMER K, LANGE-BERTALOT H. Bacillariophyceae. teil 2: Bacillariaceae, epithemiaceae, surirellaceae [M]// ETTL H, G?RTNER G, GERLOFF J, et al. Süsswasserflora von Mitteleuropa, 2(2). New York: Gustav Fischer Verlag, Stuttgart, 1988: 300–421.

    [11] KRAMMER K, LANGE-BERTALOT H. Bacillariophyceae. teil 3: Centrales, fragilariaceae, eunotiaceae [M]// ETTL H, G?RTNER G, GERLOFF J, et al. Süsswasserflora von Mitteleuropa, 2(3). Jena: Gustav Fischer Verlag, Stuttgart, 1991: 320–325.

    [12] KRAMMER K, LANGE-BERTALOT H. Bacillariophyceae. teil 4: Achnanthaceae. kritische erg?nzungen zu(Lineolatae) und[M]// ETTL H, G?RTNER G, GERLOFF J, et al. Süss- wasserflora von Mitteleuropa, 2(4). Jena: Gustav Fischer Verlag, Stutt- gart, 1991: 276–422.

    [13] JIN D X, CHEN Z D, LIN J M, et al. Benthic Diatoms in the Oceans of China, Volume 1 [M]. Beijing: China Ocean Press, 1982: 102–322.

    金德祥, 程兆第, 林均民, 等. 中國(guó)海洋底棲硅藻類(上卷) [M]. 北京: 海洋出版社, 1982: 102–322.

    [14] LIANG Y Y, LI D L, SHA L B, et al. Early to mid-Holocene sedimentary environment evolution in Xianghu Lake Area, Zhejiang Province, based on a diatom record [J]. Quat Sci, 2018, 38(4): 842–853. doi: 10.11928/j.issn.1001-7410.2018.04.04.

    梁鈺瑩, 李冬玲, 沙龍濱, 等. 浙江湘湖早-中全新世的硅藻記錄及沉積環(huán)境演變 [J]. 第四紀(jì)研究, 2018, 38(4): 842–853. doi: 10. 11928/j.issn.1001-7410.2018.04.04.

    [15] SHA L B, LIU Y G, LI D L, et al. Diatoms from the surface sediments off west Greenland and their relationship to sea-ice cover [J]. Acta Micropalaeontol Sin, 2012, 29(4): 321–332.

    沙龍濱, 劉焱光, 李冬玲, 等. 格陵蘭西部海域表層沉積硅藻分布與海冰覆蓋率關(guān)系探討 [J]. 微體古生物學(xué)報(bào), 2012, 29(4): 321–332.

    [16] SHEN L N, CHEN M, LAN B B, et al. Diatom distribution as an environmental indicator in surface sediments of the West Philippine Basin [J]. Chin J Ocean Limnol, 2017, 35(2): 431–443. doi: 10.1007/ s00343-016-5306-8.

    [17] WANG S S, LIU D Y, WANG Y Y, et al. Temporal and spatial variations of benthic diatom communities at three estuaries in the Bohai Sea [J]. Haiyang Xuebao, 2020, 42(8): 101–114. doi: 10.3969/j. issn.0253-4193.2020.08.009.

    王珊珊, 劉東艷, 王玉鈺, 等. 渤海3個(gè)河口區(qū)底棲硅藻群落的時(shí)空變化特征 [J]. 海洋學(xué)報(bào), 2020, 42(8): 101–114. doi: 10.3969/j.issn. 0253-4193.2020.08.009.

    [18] HUANG Y. Distribution of the surface sediment diatoms in the outer bay of Qinzhou bay of Guangxi [J]. Mar Sci, 2017, 41(1): 96–103. doi: doi: 10.11759/hykx20150927001.黃玥. 廣西欽州灣外灣表層沉積硅藻分布特征 [J]. 海洋科學(xué), 2017, 41(1): 96–103. doi: 10.11759/hykx20150927001.

    [19] SHANG Z W, TIAN L Z, WANG H, et al. Diatom assemblages from surficial sediments in north-central Bohai Bay and their implications for environments [J]. Geol China, 2012, 39(4): 1099–1107. doi: 10. 3969/j.issn.1000-3657.2012.04.027.

    商志文, 田立柱, 王宏, 等. 渤海灣中北部表層沉積硅藻分布及環(huán)境指示意義 [J]. 中國(guó)地質(zhì), 2012, 39(4): 1099–1107. doi: 10.3969/j. issn.1000-3657.2012.04.027.

    [20] HUANG Y, HUANG Y H. Characterastics of surface sediments diatom distribution in Zhenzhu Bay of Guangxi [J]. Adv Marine Sci, 2016, 34 (3): 411–4 20. doi: 10.3969/j.issn.1671-6647.0000.00.011.

    黃玥, 黃元輝. 廣西珍珠灣表層沉積硅藻分布特征 [J]. 海洋科學(xué)進(jìn)展, 2016, 34(3): 411–420. doi: 10.3969/j.issn.1671-6647.0000.00.011.

    [21] HUANG Y H, HUANG Y, JIANG H. Surface water temperature changes in the northern South China Sea over last ca. 15 000 years: Evidence from marine diatoms [J]. Mar Geol Quat Geol, 2007, 27(5): 65–74. doi: 10. 16562/j.cnki.0256-1492.2007.05.002.

    黃元輝, 黃玥, 蔣輝. 南海北部15kaBP以來表層海水溫度變化: 來自海洋硅藻的記錄 [J]. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2007, 27(5): 65–74. doi: 10.16562/j.cnki.0256-1492.2007.05.002.

    [22] GAO Y H, CHEN C P, LI Y, et al. A preliminary study on nanodiatoms from Holocene sediments in the Jiulong River and Dongyu Island of Fujian Province, China [J]. Acta Oceanol Sin, 2003, 25(6): 89–98. doi: 10.3321/j.issn:0253-4193.2003.06.011.

    高亞輝, 陳長(zhǎng)平, 李揚(yáng), 等. 福建九龍江和東嶼地區(qū)全新世微型硅藻的初步研究 [J]. 海洋學(xué)報(bào), 2003, 25(6): 89–98. doi: 10.3321/j.issn: 0253-4193.2003.06.011.

    [23] ZHAO L. Taxonomic and ecological studies on common epipsammic diatoms in the coasts of China [D]. Xiamen: Xiamen University, 2016: 85.

    趙龍. 中國(guó)沿海常見沙生境硅藻的分類學(xué)與生態(tài)學(xué)研究 [D]. 廈門: 廈門大學(xué), 2016: 85.

    [24] WU Y N, LI D L, SHA L B, et al. A diatom record of rapid sea-level rise during early to mid-Holocene in Baima lake area, Zhejiang Province [J]. Acta Palaeontol Sin, 2016, 55(4): 495–507. doi: 10.19800/j.cnki. aps.2016.04.010.

    吳燕妮, 李冬玲, 沙龍濱, 等. 浙江白馬湖地區(qū)硅藻記錄的早-中全新世海平面上升事件 [J]. 古生物學(xué)報(bào), 2016, 55(4): 495–507. doi: 10.19800/j.cnki.aps.2016.04.010.

    [25] CAO Q. Seasonal change of diatoms and climate response in recent 800 years of Douhu, Fujian [D]. Fuzhou: Fujian Normal University, 2018: 21.

    曹奇. 福建斗湖硅藻季節(jié)變化及近800年來的氣候響應(yīng) [D]. 福州: 福建師范大學(xué), 2018: 21.

    [26] CHEN C. Distribution of diatoms in surface sediment in Chinese inshore waters and their environmental significance [D]. Qingdao: Ocean University of China, 2013: 90.

    陳淳. 中國(guó)近岸海域表層沉積硅藻的地理分布特征及其環(huán)境意義 [D]. 青島: 中國(guó)海洋大學(xué), 2013: 90.

    [27] LAN B B. Diatoms and their palaeoenvironments from the cores of Xijiang Delta in China [D]. Xiamen: Third Institute of Oceanography, SOA, 2007: 47.

    蘭彬斌. 西江三角洲巖芯中的硅藻及其記錄的古環(huán)境 [D]. 廈門: 國(guó)家海洋局第三海洋研究所, 2007: 47.

    Diatom Distribution in Surface Sediments of Aojiang River Estuary in Fujian, China

    ZHANG Jihui1a, ZHOU Chengxu1b, LI Dongling1a, LIN Zhongzhou2, QIU Yue1a, SHA Longbin1a*

    (1a. Department of Geography and Spatial Information Techniques; 1b. College of Food and Pharmaceutical Sciences, Ningbo University,Ningbo 315211, Zhejiang, China; 2. Ningbo Institute of Oceanography,Ningbo 315000, Zhejiang, China)

    To understand the spatial distribution characteristic of diatoms in Aojiang River estuary area, Fujian, surface sediment samples from 13 sites were collected in July 2019, and the relationship between distribution of diatoms and environmental factors was studied. The results showed that there were 114 species of diatoms belonging to 39 genera identified from sediment samples. The abundance of diatoms from estuarine to offshore was firstly decreased and then increased. Cluster analysis showed that the main diatom species could be divided into fresh water assemblage zone and brackish-brackish water assemblage zone. Surface water temperature and salinity were the main environmental variables that affect the distribution of diatoms by canonical correspondence analysis (CCA). The salinity in the estuary area was relatively low, freshwater diatoms, such as,andwere the dominant species. The salinity of northern offshore area of Aojiang River was high, marine diatomwas dominant. However, the salinity of southern offshore area was low, the abundance of freshwater diatomwas higher than that in northern area. In addition, it is speculated that there may be some water pollution in estuary and the southern offshore area.

    Aojiang River; Diatom; Environment; Canonical correspondence analysis; Cluster analysis

    10.11926/jtsb.4375

    2021-01-08

    2021-04-30

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2018YFD0900702); 國(guó)家自然科學(xué)基金項(xiàng)目(41776193, 41876215); 浙江省自然科學(xué)基金項(xiàng)目(Y19D010011)資助

    This work was supported by the Project for Key Research & Development of China (Grant No. 2018YFD0900702), the National Natural Science Foundation of China (Grant No. 41776193, 41876215), and the Project for Natural Science in Zhejiang (Grant No. Y19D010011).

    張紀(jì)暉(1998~ ), 男, 主要從事海洋環(huán)境研究。E-mail: zhangjihui0@163.com

    通信作者 Corresponding author. E-mail: shalongbin@nbu.edu.cn

    猜你喜歡
    河口區(qū)環(huán)境變量江口
    基于最大熵模型的云南思茅松潛在分布區(qū)
    鹿寨江口寶塔與柳宗元詩(shī)
    文史春秋(2022年4期)2022-06-16 07:13:26
    漢江河口區(qū)四明灘段四明上護(hù)岸坍岸修復(fù)技術(shù)
    從桌面右鍵菜單調(diào)用環(huán)境變量選項(xiàng)
    “江口沉銀”的時(shí)代價(jià)值探析
    徹底弄懂Windows 10環(huán)境變量
    “傳承—信仰—未來”教育主題的設(shè)計(jì)者
    ——記山東省東營(yíng)市河口區(qū)河安小學(xué)校長(zhǎng)劉波
    教書育人(2017年23期)2017-09-26 06:47:50
    基于三階段DEA—Malmquist模型的中國(guó)省域城鎮(zhèn)化效率測(cè)度及其收斂分析
    閩江下游及河口區(qū)表層沉積物多種微量元素特征
    泉州灣洛陽江口2種紅樹林生境大型底棲動(dòng)物群落多樣性比較
    美女主播在线视频| 国模一区二区三区四区视频| 精华霜和精华液先用哪个| 精品人妻偷拍中文字幕| 成人亚洲精品一区在线观看 | 久久久成人免费电影| 日本午夜av视频| 欧美日韩视频精品一区| 亚洲av中文av极速乱| 男女啪啪激烈高潮av片| 麻豆久久精品国产亚洲av| 色综合色国产| 国产成人一区二区在线| 另类亚洲欧美激情| 男女那种视频在线观看| 最近最新中文字幕大全电影3| 亚洲成人一二三区av| 建设人人有责人人尽责人人享有的 | 秋霞在线观看毛片| 在线免费十八禁| 精品人妻熟女av久视频| 中文乱码字字幕精品一区二区三区| 亚洲欧美中文字幕日韩二区| 美女cb高潮喷水在线观看| 女人十人毛片免费观看3o分钟| 性插视频无遮挡在线免费观看| 欧美性猛交╳xxx乱大交人| 看非洲黑人一级黄片| 哪个播放器可以免费观看大片| 久久国产乱子免费精品| 国产色爽女视频免费观看| 伊人久久国产一区二区| 久久久a久久爽久久v久久| 小蜜桃在线观看免费完整版高清| 久久久久国产网址| 人体艺术视频欧美日本| 成人毛片a级毛片在线播放| 91午夜精品亚洲一区二区三区| 成人亚洲欧美一区二区av| 你懂的网址亚洲精品在线观看| 白带黄色成豆腐渣| 内射极品少妇av片p| 色综合色国产| 全区人妻精品视频| av网站免费在线观看视频| 午夜福利在线在线| 免费高清在线观看视频在线观看| 欧美日韩在线观看h| 日韩精品有码人妻一区| 最新中文字幕久久久久| 国产 一区精品| 美女国产视频在线观看| 亚洲国产欧美在线一区| 夫妻午夜视频| 久久久久久国产a免费观看| 九九久久精品国产亚洲av麻豆| 国产精品秋霞免费鲁丝片| 精品久久久久久电影网| 高清日韩中文字幕在线| 亚洲av免费在线观看| 七月丁香在线播放| 亚洲怡红院男人天堂| 久久精品久久久久久噜噜老黄| 婷婷色麻豆天堂久久| 国产亚洲5aaaaa淫片| 日本色播在线视频| 久久久久精品久久久久真实原创| 国产精品av视频在线免费观看| 国产精品国产三级国产av玫瑰| 国产男人的电影天堂91| 国产成年人精品一区二区| 2018国产大陆天天弄谢| 激情 狠狠 欧美| 国产av国产精品国产| 亚洲国产精品成人综合色| 毛片一级片免费看久久久久| 校园人妻丝袜中文字幕| 午夜福利视频精品| 亚洲欧美中文字幕日韩二区| 美女高潮的动态| 亚洲不卡免费看| 日日摸夜夜添夜夜爱| 美女cb高潮喷水在线观看| 亚洲精品第二区| 久久久久久久国产电影| 国产探花在线观看一区二区| 国产乱人偷精品视频| 下体分泌物呈黄色| 国产美女午夜福利| 日韩强制内射视频| 亚洲精品成人久久久久久| 1000部很黄的大片| 国产精品99久久久久久久久| 国产 一区 欧美 日韩| 午夜激情福利司机影院| 嫩草影院新地址| 免费观看a级毛片全部| 麻豆成人av视频| 极品少妇高潮喷水抽搐| 午夜爱爱视频在线播放| 国产亚洲最大av| 国产精品秋霞免费鲁丝片| 制服丝袜香蕉在线| 色5月婷婷丁香| 国产在线一区二区三区精| 人妻少妇偷人精品九色| 免费大片18禁| 国产精品一及| 男女边摸边吃奶| 伦理电影大哥的女人| 大又大粗又爽又黄少妇毛片口| 久久精品熟女亚洲av麻豆精品| 在线精品无人区一区二区三 | 美女被艹到高潮喷水动态| 欧美丝袜亚洲另类| 亚洲国产成人一精品久久久| 精品久久久噜噜| 成人免费观看视频高清| 麻豆久久精品国产亚洲av| 国产精品成人在线| 亚洲av免费在线观看| 天美传媒精品一区二区| 99久久精品国产国产毛片| 色视频www国产| 久久久久久久久久久丰满| 80岁老熟妇乱子伦牲交| 女人久久www免费人成看片| 少妇人妻一区二区三区视频| 久久久国产一区二区| 欧美zozozo另类| 一级黄片播放器| 在线观看人妻少妇| 国产毛片在线视频| 永久网站在线| 欧美bdsm另类| 久久久国产一区二区| 校园人妻丝袜中文字幕| 好男人在线观看高清免费视频| 香蕉精品网在线| 欧美高清成人免费视频www| 1000部很黄的大片| www.色视频.com| 亚洲精品久久午夜乱码| 一级黄片播放器| 精品视频人人做人人爽| 永久网站在线| 色吧在线观看| 男男h啪啪无遮挡| 国产亚洲一区二区精品| 特大巨黑吊av在线直播| 成人综合一区亚洲| 国产黄片美女视频| 国产精品偷伦视频观看了| 精品午夜福利在线看| 亚洲欧美日韩卡通动漫| 最近的中文字幕免费完整| 亚洲图色成人| 欧美人与善性xxx| 成人毛片60女人毛片免费| 只有这里有精品99| 黑人高潮一二区| 日本爱情动作片www.在线观看| 久久久色成人| 亚洲成人一二三区av| 极品少妇高潮喷水抽搐| 在线精品无人区一区二区三 | 夜夜看夜夜爽夜夜摸| 日本午夜av视频| 搡老乐熟女国产| 久久国内精品自在自线图片| 欧美成人精品欧美一级黄| 日韩一区二区视频免费看| 久久久久久久久久久免费av| 五月天丁香电影| 欧美xxxx黑人xx丫x性爽| 简卡轻食公司| 欧美成人精品欧美一级黄| 中国美白少妇内射xxxbb| 亚洲第一区二区三区不卡| 3wmmmm亚洲av在线观看| 激情 狠狠 欧美| 欧美 日韩 精品 国产| 啦啦啦啦在线视频资源| 成人综合一区亚洲| 99久久九九国产精品国产免费| 亚洲欧美精品自产自拍| 婷婷色综合大香蕉| 欧美性感艳星| 亚洲精品日本国产第一区| 成人鲁丝片一二三区免费| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久久电影| 国产黄频视频在线观看| 一边亲一边摸免费视频| 亚洲aⅴ乱码一区二区在线播放| 久久久久国产精品人妻一区二区| 国产成人a区在线观看| 看十八女毛片水多多多| 神马国产精品三级电影在线观看| 高清av免费在线| 一区二区三区乱码不卡18| 欧美成人精品欧美一级黄| 亚洲精华国产精华液的使用体验| 插逼视频在线观看| 自拍偷自拍亚洲精品老妇| 搡女人真爽免费视频火全软件| 久久综合国产亚洲精品| 少妇熟女欧美另类| 国产高清三级在线| 国产精品一区二区三区四区免费观看| 亚洲精品视频女| 爱豆传媒免费全集在线观看| 一本久久精品| 乱码一卡2卡4卡精品| 少妇高潮的动态图| 亚洲熟女精品中文字幕| 18禁在线播放成人免费| 亚洲国产色片| 国产精品麻豆人妻色哟哟久久| 白带黄色成豆腐渣| av播播在线观看一区| av一本久久久久| 麻豆成人午夜福利视频| 亚洲不卡免费看| 亚洲最大成人中文| 久久精品人妻少妇| 婷婷色麻豆天堂久久| 亚州av有码| 亚洲高清免费不卡视频| 亚洲人成网站在线播| 十八禁网站网址无遮挡 | 日韩强制内射视频| 免费在线观看成人毛片| 午夜免费观看性视频| 亚洲国产高清在线一区二区三| 五月伊人婷婷丁香| 国产黄片视频在线免费观看| 色视频www国产| 国产伦在线观看视频一区| 久久久亚洲精品成人影院| 国产国拍精品亚洲av在线观看| 日韩亚洲欧美综合| 国产成人a区在线观看| 成人亚洲精品一区在线观看 | 亚洲性久久影院| 2022亚洲国产成人精品| 欧美少妇被猛烈插入视频| 好男人在线观看高清免费视频| 午夜精品一区二区三区免费看| av天堂中文字幕网| 五月伊人婷婷丁香| 亚洲欧美日韩卡通动漫| 高清毛片免费看| 国产一区二区三区综合在线观看 | 一级二级三级毛片免费看| 水蜜桃什么品种好| 国产欧美日韩精品一区二区| 人妻夜夜爽99麻豆av| 久久人人爽人人爽人人片va| 亚洲综合精品二区| 直男gayav资源| 伦理电影大哥的女人| 欧美bdsm另类| 一二三四中文在线观看免费高清| 女人被狂操c到高潮| 搞女人的毛片| 五月伊人婷婷丁香| 超碰av人人做人人爽久久| 美女高潮的动态| 麻豆成人av视频| 国产精品福利在线免费观看| 欧美性猛交╳xxx乱大交人| 国产精品一二三区在线看| 免费看不卡的av| 久久久午夜欧美精品| 丰满乱子伦码专区| 国产精品一区www在线观看| 联通29元200g的流量卡| 国内精品美女久久久久久| 只有这里有精品99| 国产精品国产三级国产av玫瑰| 天天躁日日操中文字幕| 秋霞在线观看毛片| 晚上一个人看的免费电影| 一级爰片在线观看| av天堂中文字幕网| 久久久久久九九精品二区国产| 2018国产大陆天天弄谢| 男人和女人高潮做爰伦理| 在线观看一区二区三区激情| 免费观看在线日韩| 日韩大片免费观看网站| 禁无遮挡网站| 99热这里只有是精品50| 日本免费在线观看一区| 国产精品人妻久久久影院| 午夜老司机福利剧场| av在线app专区| 免费观看av网站的网址| 97超视频在线观看视频| 日韩亚洲欧美综合| 国产又色又爽无遮挡免| av在线观看视频网站免费| 国产精品av视频在线免费观看| 在线观看美女被高潮喷水网站| 又爽又黄无遮挡网站| 日韩,欧美,国产一区二区三区| 黄片无遮挡物在线观看| 欧美一区二区亚洲| 亚洲成人一二三区av| 国产精品一区二区三区四区免费观看| 日韩精品有码人妻一区| 国产成人a∨麻豆精品| 国产在线一区二区三区精| 99精国产麻豆久久婷婷| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美清纯卡通| av在线老鸭窝| 欧美精品人与动牲交sv欧美| 大片电影免费在线观看免费| 黄色视频在线播放观看不卡| 秋霞伦理黄片| 国产精品无大码| 国模一区二区三区四区视频| 国产一区有黄有色的免费视频| 人妻夜夜爽99麻豆av| 狂野欧美激情性xxxx在线观看| 免费看av在线观看网站| 成人免费观看视频高清| 午夜激情久久久久久久| 精品久久国产蜜桃| 久久精品国产自在天天线| 自拍偷自拍亚洲精品老妇| 黄片无遮挡物在线观看| 免费大片黄手机在线观看| 久久精品国产鲁丝片午夜精品| 久久6这里有精品| 日韩制服骚丝袜av| 少妇高潮的动态图| 日韩国内少妇激情av| 神马国产精品三级电影在线观看| 国产成人a∨麻豆精品| 乱码一卡2卡4卡精品| 麻豆久久精品国产亚洲av| av在线天堂中文字幕| 伊人久久国产一区二区| 国产成人精品久久久久久| 成人漫画全彩无遮挡| 极品教师在线视频| 人人妻人人爽人人添夜夜欢视频 | 欧美xxxx性猛交bbbb| 国产成人freesex在线| 高清日韩中文字幕在线| 51国产日韩欧美| 久久精品国产亚洲av涩爱| 我的女老师完整版在线观看| 在线免费十八禁| 欧美老熟妇乱子伦牲交| 少妇人妻 视频| av女优亚洲男人天堂| 熟女av电影| 亚洲怡红院男人天堂| 亚洲av免费高清在线观看| 亚洲,欧美,日韩| 亚州av有码| 啦啦啦啦在线视频资源| 欧美三级亚洲精品| 我要看日韩黄色一级片| 在线亚洲精品国产二区图片欧美 | 日本与韩国留学比较| 一区二区三区四区激情视频| 最近最新中文字幕大全电影3| 国产亚洲5aaaaa淫片| 91精品伊人久久大香线蕉| 久久精品熟女亚洲av麻豆精品| 国产黄频视频在线观看| 成人无遮挡网站| 国产成人免费观看mmmm| 亚洲va在线va天堂va国产| 少妇的逼水好多| 99热网站在线观看| 亚洲精品成人av观看孕妇| 亚洲经典国产精华液单| 亚洲精品国产色婷婷电影| 亚洲不卡免费看| 国产精品成人在线| 午夜福利视频1000在线观看| 欧美日本视频| 最近最新中文字幕免费大全7| 成年免费大片在线观看| 你懂的网址亚洲精品在线观看| 国产视频首页在线观看| 2021少妇久久久久久久久久久| 人妻 亚洲 视频| 国产黄频视频在线观看| 毛片女人毛片| 特级一级黄色大片| 五月玫瑰六月丁香| 亚洲精品456在线播放app| 在线天堂最新版资源| 超碰av人人做人人爽久久| 国产精品麻豆人妻色哟哟久久| 国产一区二区三区av在线| 麻豆成人午夜福利视频| 高清午夜精品一区二区三区| 少妇被粗大猛烈的视频| 国产精品国产三级国产av玫瑰| 狂野欧美激情性bbbbbb| 午夜福利高清视频| 日韩成人av中文字幕在线观看| 18+在线观看网站| 日本wwww免费看| 香蕉精品网在线| 2021天堂中文幕一二区在线观| 国产精品99久久99久久久不卡 | 可以在线观看毛片的网站| 99热网站在线观看| 国产久久久一区二区三区| 亚洲av欧美aⅴ国产| 在现免费观看毛片| 插逼视频在线观看| 国产成人a∨麻豆精品| 久久鲁丝午夜福利片| 青春草亚洲视频在线观看| 精品久久久噜噜| 亚洲av.av天堂| 国产v大片淫在线免费观看| 亚洲第一区二区三区不卡| 久久久精品欧美日韩精品| 嘟嘟电影网在线观看| 国产91av在线免费观看| 亚洲欧美精品自产自拍| 久久99热这里只频精品6学生| 国产淫语在线视频| 亚洲国产av新网站| 韩国高清视频一区二区三区| 欧美人与善性xxx| 国产黄a三级三级三级人| 狠狠精品人妻久久久久久综合| 日韩欧美 国产精品| 欧美成人一区二区免费高清观看| 尾随美女入室| 国产精品人妻久久久影院| 国产爽快片一区二区三区| 免费高清在线观看视频在线观看| 欧美成人一区二区免费高清观看| 国产成人一区二区在线| 久久人人爽人人爽人人片va| 亚洲不卡免费看| 亚洲精品乱码久久久久久按摩| 亚洲四区av| 亚洲av二区三区四区| 免费大片18禁| 黄色欧美视频在线观看| 免费av毛片视频| 丝瓜视频免费看黄片| 亚洲成人中文字幕在线播放| 一级黄片播放器| 色网站视频免费| 亚洲欧美日韩卡通动漫| 国产成人一区二区在线| 精品久久久精品久久久| 男插女下体视频免费在线播放| a级一级毛片免费在线观看| 精品一区二区免费观看| 成人二区视频| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 最近中文字幕2019免费版| 一区二区三区四区激情视频| 中国美白少妇内射xxxbb| 一级a做视频免费观看| 观看免费一级毛片| 精品国产乱码久久久久久小说| a级毛色黄片| 欧美xxⅹ黑人| 王馨瑶露胸无遮挡在线观看| 看免费成人av毛片| 久久综合国产亚洲精品| 久久精品熟女亚洲av麻豆精品| 中国三级夫妇交换| 亚洲精品国产av蜜桃| 午夜福利在线在线| 久久久久精品性色| 欧美三级亚洲精品| av国产久精品久网站免费入址| 看免费成人av毛片| 免费大片18禁| kizo精华| 日韩强制内射视频| 久久久精品免费免费高清| xxx大片免费视频| 国产综合懂色| 免费av毛片视频| 久久久久久久国产电影| 深夜a级毛片| 国产亚洲一区二区精品| freevideosex欧美| 亚洲精品第二区| 精品人妻偷拍中文字幕| 久久久色成人| 全区人妻精品视频| 亚洲高清免费不卡视频| 九草在线视频观看| 最近中文字幕高清免费大全6| 91狼人影院| 亚洲av一区综合| 久久亚洲国产成人精品v| 噜噜噜噜噜久久久久久91| 看黄色毛片网站| 男女边吃奶边做爰视频| 国产伦精品一区二区三区四那| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 亚洲精品乱久久久久久| av黄色大香蕉| 国产成人免费观看mmmm| 久久精品国产鲁丝片午夜精品| 男男h啪啪无遮挡| 免费看光身美女| 亚洲av一区综合| 搞女人的毛片| 在线观看国产h片| 国产精品不卡视频一区二区| 搡老乐熟女国产| 国产精品人妻久久久久久| 久久久精品94久久精品| 国产欧美日韩精品一区二区| 精品久久国产蜜桃| 人妻 亚洲 视频| 久久精品综合一区二区三区| 国产一区有黄有色的免费视频| 最近最新中文字幕免费大全7| 国产大屁股一区二区在线视频| 亚洲国产精品专区欧美| 97精品久久久久久久久久精品| 色网站视频免费| 最近手机中文字幕大全| 精品国产一区二区三区久久久樱花 | av播播在线观看一区| 18禁动态无遮挡网站| 少妇人妻 视频| 尤物成人国产欧美一区二区三区| 日韩制服骚丝袜av| av在线app专区| 五月玫瑰六月丁香| 亚洲三级黄色毛片| 男女边吃奶边做爰视频| 80岁老熟妇乱子伦牲交| 嫩草影院精品99| 久久久久久久久久成人| 最近中文字幕高清免费大全6| 伦精品一区二区三区| av国产精品久久久久影院| 黄色一级大片看看| 色播亚洲综合网| 色网站视频免费| 婷婷色综合www| 永久免费av网站大全| 午夜爱爱视频在线播放| 三级国产精品欧美在线观看| 在现免费观看毛片| 久久热精品热| 建设人人有责人人尽责人人享有的 | 亚洲精品国产av蜜桃| 日韩不卡一区二区三区视频在线| 久久精品国产鲁丝片午夜精品| 各种免费的搞黄视频| 国产色婷婷99| 日本猛色少妇xxxxx猛交久久| 亚洲精华国产精华液的使用体验| av在线天堂中文字幕| 欧美+日韩+精品| 午夜福利视频精品| 一级毛片黄色毛片免费观看视频| 春色校园在线视频观看| 我的老师免费观看完整版| 水蜜桃什么品种好| 亚洲av中文av极速乱| 建设人人有责人人尽责人人享有的 | 亚洲av在线观看美女高潮| 久久久久久伊人网av| 成人欧美大片| 一级毛片久久久久久久久女| 色综合色国产| 婷婷色综合www| 永久免费av网站大全| 狂野欧美激情性bbbbbb| 国产精品国产三级专区第一集| 一级av片app| 少妇丰满av| 亚洲精品一区蜜桃| 免费电影在线观看免费观看| 精品久久久噜噜| 国产高清国产精品国产三级 | 人人妻人人澡人人爽人人夜夜| 亚洲高清免费不卡视频| 亚洲欧美精品自产自拍| 久久精品夜色国产| 国产午夜精品久久久久久一区二区三区| 国产一级毛片在线| 亚洲伊人久久精品综合| 联通29元200g的流量卡| 成人漫画全彩无遮挡| 亚洲色图av天堂| 男人和女人高潮做爰伦理| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美中文字幕日韩二区| 只有这里有精品99| 丝瓜视频免费看黄片| 久久久国产一区二区| 97超碰精品成人国产| 亚洲色图综合在线观看| 国产精品久久久久久久电影| 国产精品99久久99久久久不卡 |