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

    基于正交試驗的沉積物-水界面營養(yǎng)鹽交換通量研究——以海州灣海洋牧場為例

    2017-11-23 01:41:32張俊波
    中國環(huán)境科學 2017年11期
    關(guān)鍵詞:粉砂營養(yǎng)鹽沉積物

    張 碩,方 鑫,黃 宏,張 虎,張俊波,5,6

    ?

    基于正交試驗的沉積物-水界面營養(yǎng)鹽交換通量研究——以海州灣海洋牧場為例

    張 碩1,2,方 鑫1,黃 宏3,張 虎4,張俊波1,5,6*

    (1.上海海洋大學海洋科學學院,上海201306;2.上海海洋大學大洋漁業(yè)資源可持續(xù)開發(fā)省部共建教育部重點實驗室,上海201306;3.上海海洋大學海洋生態(tài)與環(huán)境學院,上海201306;4.江蘇省海洋水產(chǎn)研究所,江蘇南通226007;5.上海海洋大學國家遠洋漁業(yè)工程技術(shù)研究中心,上海201306;6.上海海洋大學水產(chǎn)科學國家級實驗教學示范中心,上海201306)

    為了研究沉積物-水界面營養(yǎng)鹽交換特性,采用正交試驗設(shè)計,對海州灣海洋牧場2016年5月3個站位沉積物-水界面營養(yǎng)鹽交換通量進行了研究,分析了沉積物類型、溫度、DO和pH值對沉積物-水界面營養(yǎng)鹽交換通量的影響.結(jié)果表明:影響營養(yǎng)鹽交換通量的因素間存在主次關(guān)系,影響NH4+-N交換的因素依次:DO>溫度>沉積物類型;影響NO3-+NO2--N交換的因素依次為:沉積物類型>DO>溫度;影響PO43--P交換的因素依次為:DO>沉積物類型>溫度;影響SiO32--Si交換的因素依次為:溫度>pH值.因素之間的交互作用對各項營養(yǎng)鹽交換通量有著重要的影響,在建立沉積物-水界面營養(yǎng)鹽交換模型時,更應著重考慮交互作用的影響效果.本研究實驗得出的沉積物類型、溫度、溶解氧和pH值這4種因素對海州灣營養(yǎng)鹽交換影響的結(jié)果與海州灣往年調(diào)查研究結(jié)果大致相符.

    沉積物-水界面;營養(yǎng)鹽;交換通量;正交試驗;交互作用

    沉積物和水體之間的相互作用(營養(yǎng)物質(zhì)的循環(huán)、遷移和存儲)對近岸淺海生態(tài)環(huán)境有重要的影響[1].沉積物-水界面是近岸淺海區(qū)域中非常重要的界面之一,沉積物-水界面營養(yǎng)鹽的交換對維持海洋初級生產(chǎn)力和水體營養(yǎng)鹽含量平衡有著顯著的影響[2].控制沉積物-水界面間營養(yǎng)鹽遷移的因素主要包括擴散、吸附-解吸、沉淀(礦化)、溶解和有機質(zhì)分解等,沉積物中有機物通過微生物的作用轉(zhuǎn)化成間隙水中的營養(yǎng)鹽,再通過吸附-解吸、自生礦物沉淀、微生物的吸收利用及擴散等作用從沉積物中發(fā)生轉(zhuǎn)移[3].在這些過程中,溫度、溶解氧、pH值、鹽度、生物擾動和沉積物自身理化性質(zhì)等因素也會促進或抑制沉積物-水界面間營養(yǎng)的交換[4].沉積物的粒度是影響TP(總磷)、TN(總氮)含量的重要因素,沉積物顆粒越細,TP、TN含量越高,沉積物粒度還制約著沉積物中TP、TN的垂直分布,其含量隨深度增加而逐漸降低[5].有研究表明,無機磷占TP的50%~90%,無機磷通過厭氧釋放轉(zhuǎn)化為PO43-[6-7],有機氮占TN的70%~90%,大部分有機氮會在微生物的作用下發(fā)生礦化作用,其主要產(chǎn)物是NH4+-N[8-9],后經(jīng)硝化作用氧化為NO3--N或NO2--N,或通過擴散作用滲入間隙水中,從而改變N、P營養(yǎng)鹽在沉積物-水界面的濃度梯度,發(fā)生交換作用[10-13].溫度對硝化和反硝化作用有顯著的影響,并認為硝化反應的微生物適宜溫度為20~30℃.溫度在5~30℃之間時,溫度每上升10℃,細菌的增長速率增大1倍[14],在溫度適宜的情況下,硝化作用發(fā)生后改變其原有的濃度梯度,加快了其交換速率,同時pH值升高也會加快NO3-的交換速率.SiO32-的交換速率雖然與溫度呈正相關(guān)的關(guān)系,但同時還會受到沉積物粒度、鹽度等因素的影響[15].溫度的升高和較低的含氧量會加快PO43-從有機質(zhì)中的釋放,沉積物的類型也會影響PO43-的釋放過程[16],在厭氧條件下,PO43-的釋放速率則是好氧條件下的2~8倍[17].此外,溫度和DO之間呈顯著負相關(guān)[18],pH值與海水溫度存在顯著的相關(guān)性[19],從另一個側(cè)面影響著營養(yǎng)鹽的交換.

    海州灣位于蘇魯交界處,該海域自2002年開始實施人工魚礁項目,經(jīng)過10多年的建設(shè),人工魚礁對海州灣漁場漁業(yè)資源的恢復起到了重要作用.從近些年的水質(zhì)生態(tài)調(diào)查報告結(jié)果看[20-21],海州灣的生態(tài)環(huán)境得到顯著的改善,海洋牧場的功能逐漸得到體現(xiàn).而海洋牧場對環(huán)境的調(diào)控能力取決于該海域營養(yǎng)物質(zhì)的交換與補充.因此,開展對海洋沉積物和海水中營養(yǎng)鹽交換的特性研究顯得尤為重要.海州灣是一個半開敞型海灣[22],灣內(nèi)有多條入海河流,會導致海水的pH值產(chǎn)生波動,且海州灣氣候四季分明,有研究表明海州灣C、N元素的賦存受到沉積物粒度的制約[23].目前,多數(shù)學者僅就單個因素對營養(yǎng)鹽交換的影響規(guī)律進行了研究,董慧[4]和Zhen[24]通過單因素的控制實驗發(fā)現(xiàn),NH4+-N、NO3--N和PO43--P在貧氧環(huán)境下的交換通量要高于富氧條件下,隨著溫度的升高,NH4+-N、NO3--N和PO43--P的交換通量增大.但影響沉積物-水界面營養(yǎng)鹽交換的因素之間并不是相互獨立的,它們會共同作用于交換過程.故本研究采用統(tǒng)計學中的正交試驗法設(shè)計實驗室培養(yǎng)實驗,從多因素(沉積物粒徑、溫度、溶解氧、pH值)角度,開展沉積物-水界面營養(yǎng)鹽交換特性研究,從而為該海域建立營養(yǎng)鹽交換模型提供基本參數(shù),也為海域水環(huán)境修復提供理論依據(jù).

    1 材料與方法

    1.1 研究站點與樣品采集

    為了解海州灣海洋牧場沉積物物理性質(zhì),本研究于2015年5月,由近及遠在海州灣魚礁區(qū)選取6個站位(RA1、RA2、RA3、RA4、RA5、RA6),對照區(qū)選擇2個站位(CA1、CA2)進行現(xiàn)場采樣調(diào)查,具體站位如圖1所示,在每個站位使用抓斗式采泥器采集4組表層沉積物.根據(jù)2015年采樣得到的沉積物類型結(jié)果,又于2016年5月,從8個站位中選取3個站位(CA1、RA1、RA4)進行現(xiàn)場調(diào)查采樣,在每個站位使用柱狀采泥器(內(nèi)置可替換的PVC管)采集站位處的27根無擾動的表層沉積物(共計81根),每根PVC管直徑為5cm,長度為35cm,且在采樣時保證每根PVC管內(nèi)的沉積物長度大于20cm,如若未達到20cm則重新采樣,把PVC管的兩頭用橡膠塞子密封,放入加冰的便攜式冷藏柜中冷藏并帶回實驗室.同時采集3桶(75L)沉積物上方的底層海水,放入冷藏柜中帶回實驗室.

    圖1 海州灣海洋牧場采樣站位

    1.2 正交試驗設(shè)計

    根據(jù)2015年5月調(diào)查采樣得到的結(jié)果, CA1、RA1、RA2站位的沉積物類型分別為粉砂質(zhì)黏土(TY)、黏土質(zhì)粉砂(YT)、粉砂質(zhì)砂(TS),平均粒徑依次增大,水深分別為12.6m、14.1m、14.5m. 908近海專項調(diào)查數(shù)據(jù)顯示:海州灣海域全年表層水溫變化在31.00~5.70℃之間,平均溫度為17.58℃,底層水溫在30.00~5.62℃之間,平均溫度為16.78℃;表層溶解氧含量介于3.59~ 11.03mg/L,平均含量為7.33mg/L,底層溶解氧含量介于3.29~8.67mg/L,平均含量為7.03mg/L;表層pH值介于7.45~8.96,平均值為8.02,底層pH值介于7.81~8.34,平均值為8.09(908近海專項調(diào)查數(shù)據(jù)).從近些年海州灣的相關(guān)研究來看,海州灣海域全年平均海表溫16.76℃,平均溫度變化范圍為5.51~28.52℃[23].海州灣平均溶解氧濃度為6.67mg/L,其具體變化范圍為3.77~11.8mg/ L[25-26],且受多條入海河流影響導致海水的pH值會上下波動,其pH平均值為8.11.除了徑流、降雨、大氣交換等因素的影響,還會受到生物生長的影響,6月和10月是海州灣赤潮的多發(fā)期,部分海域的pH值會達到8.87[27].結(jié)合908近海專項調(diào)查數(shù)據(jù)與相關(guān)研究數(shù)據(jù),采用正交設(shè)計研究沉積物類型、溫度、DO和pH值對沉積物-水界面營養(yǎng)鹽交換通量的影響,試驗共設(shè)定上述4個因素(A、B、C、D),每個因素設(shè)定3個水平,在實驗中不僅考慮單個因素的影響,而且考慮因素之間的相互作用,選取L27(313)正交設(shè)計表,具體試驗設(shè)計見表1,每組實驗設(shè)置3個平行.

    表1 正交試驗設(shè)計方案

    注:“貧氧”條件是向海水中充氮氣,DO值為(3±0.5)mg/L;“自然狀態(tài)”表示自然狀態(tài)下的DO濃度(6.5±0.5)mg/L;“富氧”條件是向海水中充空氣,DO值為(10±0.5)mg/L.

    1.3 實驗方法

    將采集的PVC管中的沉積物樣品解凍,并將其中的沉積物樣品小心推入直徑為5cm,高度為50cm的有機玻璃培養(yǎng)管中,培養(yǎng)管中沉積物的高度約為18cm,將采集的上覆水倒入2L的燒杯中(上覆水的體積約為600mL),用50%的HCL和20%的NaOH調(diào)整至所需的pH值,并充分攪勻,在沉積物上方加入約25cm的調(diào)整過pH值的上覆水,加水過程中注意避免攪動表層的沉積物.按表1調(diào)整每根培養(yǎng)管中的pH值,在需要充氮氣和充氧氣的培養(yǎng)管中放入連接著氣泵的起泡石,調(diào)節(jié)起泡石的高度,避免攪動表層沉積物,并把這些培養(yǎng)管分別放入10,20,30℃的避光培養(yǎng)箱.分別在培養(yǎng)0,2,4,8,12,24,36,48h后采集上覆水50mL,并加入等體積對應站位的海水.采集的上覆水采用0.45μm醋酸纖維膜過濾,并加入三氯甲烷冷藏保存[28],以備分析.

    1.4 測定及計算方法

    沉積物粒徑采用馬爾文Mastersizer 2000激光粒度儀測定,參照謝帕德三角圖進行分類[29].培養(yǎng)后采集的水樣采用Clever Chem380全自動間斷分析儀(Dechem-Tech)測定,各項營養(yǎng)鹽的具體方法為:NH4+的測定采用苯酚-次氯酸鹽比色法;NO3-的測定采用鎘柱還原法;NO2-的測定采用重氮-偶氮法;PO43-的測定為鉬藍分光光度法;SiO32-的測定為硅鉬藍法;沉積物-水界面營養(yǎng)鹽交換通量通過公式(1)計算計算得出[30].

    式中:為沉積物-水界面營養(yǎng)鹽的交換通量, mmol/(m2·d);()為上覆水營養(yǎng)鹽質(zhì)量的變化值, mmol;為培養(yǎng)柱的截面積,m2; Δ為培養(yǎng)時間的變化量,d.

    1.5 數(shù)據(jù)統(tǒng)計與分析

    由RGui 3.3.1軟件對營養(yǎng)鹽的結(jié)果進行多因素方差分析,主效應圖由Origin 9.1繪制,交互作用圖由Minitab 16.1.0繪制.

    2 結(jié)果與討論

    2.1 不同因素對沉積物-水界面營養(yǎng)鹽交換通量的影響

    營養(yǎng)鹽的交換通量方差分析結(jié)果如表2所示,在影響營養(yǎng)鹽交換通量的各因素之間存在著交互作用的影響,且這些因素對各營養(yǎng)鹽交換通量的影響程度有所不同,影響程度的顯著性差異較大.其中DO對NH4+-N交換通量有極顯著的影響(<0.01),沉積物類型、溫度及沉積物類型和DO的交互作用對NH4+-N交換通量有顯著的影響(<0.05).通過比較表2中顯著性值的大小可以得出各個因素或交互作用對NH4+-N交換通量影響的主次關(guān)系:DO(C)>(A′C)>溫度(B)>沉積物類型(A)>pH值(D)>(C′D)>(A′D)>(B′D)> (B′C)>(A′B).通過比較因素獨立作用的顯著性及其主效應圖(圖2),得出NH4+-N交換通量隨著沉積物類型的平均粒徑增大而增大,隨著DO的升高而減小,當溫度為20℃時,其交換通量最大.DO是影響NH4+-N交換通量最顯著的因素,在有氧條件下NH4+可通過硝化作用氧化為NO2-,并能被繼續(xù)氧化為NO3-,在缺氧條件下, NO3-通過反硝化作用,NO2-和NH4+通過厭氧氨氧化逐步還原生成 N2,或NO3-通過異化硝酸鹽還原作用還原為NH4+[31].當水體中的DO濃度較大時,沉積物表層(氧滲透層)處于氧化環(huán)境,硝化作用更容易發(fā)生,沉積物間隙水中的NH4+-N經(jīng)過硝化細菌作用轉(zhuǎn)化為NO3-,導致表層沉積物間隙水中的NH4+-N濃度降低.此外,沉積物中NO3-在缺氧環(huán)境下會發(fā)生異化硝酸鹽還原作用,提高沉積物中NH4+-N的濃度[31-33].實驗中,當水體中的DO濃度降低時,沉積物表層的氧滲透層發(fā)生異化硝酸鹽還原作用,改變了原有的濃度梯度,濃度梯度的改變直接影響物質(zhì)的交換強度[34],從而增大了NH4+-N的交換通量;溫度對NH4+-N交換通量存在顯著的影響,溫度在5~30℃之間時,溫度每上升10℃,細菌的增長速率增大1倍[14].溫度對參與硝化和反硝化作用的微生物活性有顯著的影響,并認為硝化反應的微生物適宜溫度為20~30℃,反硝化的適宜溫度為25℃左右[35-37].實驗發(fā)現(xiàn)在20℃時,NH4+-N平均交換通量達到最大值,這可能是由于20℃時微生物的活性較強,沉積物間隙水中大量NH4+-N發(fā)生硝化作用,導致水體中的NH4+-N濃度大于沉積物間隙水中的濃度,從而發(fā)生交換作用.但是隨著溫度的升高,硝化作用減弱,從而減低了NH4+-N交換通量;沉積物組成粒徑較大的黏土質(zhì)粉砂(YT)和粉砂質(zhì)砂(TS)NH4+-N平均交換通量比粒徑較小的粉砂質(zhì)黏土(TY)的交換通量大,沉積物的粒度大小不僅可以改變沉積物中氧氣的分布,而且沉積物的粒徑越大,微生物能更容易接觸到顆粒的表面,因而能培養(yǎng)更多的微生物,使得更多的微生物參與到硝化反應中,從而加快了NH4+-N的硝化過程[38-40].

    如表2所示,沉積物類型及沉積物類型和DO的交互作用對NO3-+NO2--N交換通量有極顯著的影響(<0.01),溫度、DO、沉積物類型和溫度的交互作用及沉積物類型和pH的交互作用對NO3-+NO2--N交換通量有顯著的影響(< 0.05).通過比較表2中顯著性值的大小可以得出各個因素或交互作用對NO3-+NO2--N交換通量影響的主次關(guān)系:沉積物類型(A)>(A′C)>DO (C)>溫度(B)>(A′B)>(A′D)>pH值(D)>(B′D)> (B′C)>(C′D).通過比較因素獨立作用的顯著性及其主效應圖(圖2),得出當沉積物類型為黏土質(zhì)粉砂(YT)時,NO3-+NO2--N交換通量最大,NO3-+ NO2--N的交換通量隨著DO的升高而減小,當溫度為10℃,其交換通量最大.沉積物類型是影響NO3-+NO2--N交換通量最顯著的因素,沉積物組成粒徑較大的黏土質(zhì)粉砂(YT)和粉砂質(zhì)砂(TS) 中NO3-+NO2—N的平均交換通量高于粒徑較小的粉砂質(zhì)黏土(TY)的交換通量大,與NH4+-N不同的是黏土質(zhì)粉砂(YT)的交換通量比粉砂質(zhì)砂(TS)和粉砂質(zhì)黏土(TY)高,由于沉積物類型為黏土質(zhì)粉砂(YT)和粉砂質(zhì)砂(TS)時,NH4+-N表現(xiàn)出較高的交換通量,大量NH4+-N在細菌的作用下發(fā)生硝化作用轉(zhuǎn)化成為NO3--N和NO2--N[41].當沉積物表層(氧滲透層)的NO3--N濃度大于沉積物間隙水中NO3--N的濃度時,NO3--N會向沉積物間隙水擴散,在擴散過程中有一部分NO3-- N會發(fā)生反硝化作用,導致交換通量的減小.DO對NO3-+NO2--N的交換通量也有顯著的影響, DO的濃度直接關(guān)系到水體的氧化還原環(huán)境.當水體為氧化環(huán)境時,氧滲透層的NH4+-N更容易發(fā)生硝化作用,NO3-+NO2--N的交換通量應該比還原條件下更大.但在本實驗中,NO3-+NO2--N的交換通量隨溶解氧的升高而降低.原因可能是在還原條件下,沉積物對NO3--N有較強的吸附性,從而抑制NO3-+NO2--N的交換,導致其交換通量減小[42-44].但當溫度超過25℃時,硝化反應速率降低[45],致使NO3-+NO2--N的交換通量增大,從而出現(xiàn)轉(zhuǎn)折.

    如表2所示,DO對PO43--P交換通量有極顯著的影響(<0.01),沉積物類型、溫度、溫度和DO的交互作用對PO43--P交換通量有顯著的影響(<0.05).通過比較表2中顯著性值的大小可以得出各個因素或交互作用對PO43--P交換通量影響的主次關(guān)系:DO(C)>沉積物類型(A)>(B′C)>溫度(B)>(A′B)>(B′D)>(A′D)>(C′D)>(A′C)>pH值(D).通過比較因素獨立作用的顯著性及其主效應圖(圖2),得出PO43--P交換通量隨著沉積物類型的平均粒徑增大而增大,隨著DO的升高而減小,隨著溫度的升高而增大.DO是影響PO43--P交換通量最顯著的因素, P不會直接參與氧化還原過程,Fe(Ⅲ)與PO43--P雖然形成難溶的鐵結(jié)合態(tài)P,當沉積物表層(氧滲透層)含氧量較高時,Fe(Ⅲ)會限制PO43--P向水體的遷移[46],而當沉積物表層(氧滲透層)含氧量較低時,Fe(Ⅲ)被還原成Fe(Ⅱ),PO43--P從結(jié)合態(tài)P中釋放出來[32,47-48],從而導致了PO43--P的交換通量隨DO的降低而增大;沉積物的組成對PO43--P的交換通量也產(chǎn)生了顯著的影響,水體中的水生動物的排泄、浮游動植物的尸體以及外源磷沉降都會導致沉積物中有機磷含量的增高,微生物能分解有機質(zhì),促進沉積物中營養(yǎng)鹽向水層的釋放[31],較大的粒徑可以培養(yǎng)更多的微生物對有機質(zhì)進行分解;溫度對PO43--P的影響表現(xiàn)在兩個方面,沉積物-水界面上的吸附反應通常是一個放熱過程,當溫度升高時,磷酸鹽會發(fā)生解吸作用[49],溫度的升高還會提高微生物的活性,從而使PO43--P的交換通量增大.

    沉積物類型和溫度的交互作用對SiO32--Si交換通量有極顯著的影響(<0.01),溫度、pH值及沉積物類型和DO的交互作用對SiO32--Si交換通量有顯著的影響(<0.05).通過比較表2中顯著性值的大小可以得出各個因素或交互作用對SiO32--Si交換通量影響的主次關(guān)系:(A′B)> (A′C)>溫度(B)>pH值(D)>沉積物類型(A)>(A′D)>(B′C)>DO(C)>(C′D)>(B′D).通過比較因素獨立作用的顯著性及其主效應圖(圖2),得出SiO32--Si交換通量隨著溫度的升高而增大,隨著pH的升高而增大.溫度對SiO32--Si的交換通量有顯著的影響,溫度會顯著影響生物硅的溶解度和溶解動力學過程[50-52].有研究表明,生物硅的溶解度和溫度呈線性相關(guān),其關(guān)系為eq=23.8+ 936,式中:eq是生物硅的溶解度,是反應溫度,溫度變化范圍為4.5~28℃[53].生物硅的大量溶解,導致溶解態(tài)硅濃度上升,從而增大了SiO32--Si的交換通量;雖然pH值在海洋環(huán)境中相對穩(wěn)定,但在河口地區(qū),pH值的變化相對明顯[54].生物硅的溶解動力在自然水體中都隨pH值的增大而升高,而較高的pH值可以使硅醇鍵更容易斷裂[55].研究發(fā)現(xiàn)當pH值從6.1上升至8.1時,二氧化硅的溶解速度上升了2倍[56],溶解速度的增大是導致硅酸鹽交換通量變大的主要因素.

    表2 正交設(shè)計方差分析

    注:*表示有顯著性差異(<0.05);**表示有極顯著性差異(<0.01).′表示2個因素的交互作用.

    圖2 營養(yǎng)鹽交換通量主效應Fig.2 Main effect plot of nutrients exchange flux

    2.2 各因素交互作用對沉積物-水界面營養(yǎng)鹽交換通量的影響

    由表2和圖3可知,只有(A′C)的交互作用對NH4+-N的交換通量產(chǎn)生了顯著的影響,沉積物類型和DO在各自不同水平上都會有不同的效果.在單因素條件下,NH4+-N的交換通量隨沉積物組成平均粒徑的增大而減小,隨DO的增大而減小,但 (A′C)存在對NH4+-N交換通量的交互作用,在平均粒徑較小的沉積物類型和低DO水平共同作用下,NH4+-N的交換通量較大,在沉積物為黏土質(zhì)粉砂(YT)和DO為自然水平的共同作用下,NH4+-N的交換通量較小.因此,在研究NH4+-N的交換通量時,需要對沉積物類型和DO的交互作用進行考慮.

    圖3 NH4+-N的交互作用

    圖4 NO3-+NO2--N的交互作用

    由表2和圖4可知,(A′B)、(A′C)、(A′D)的交互作用都對NO3-+NO2--N的交換通量產(chǎn)生了顯著的影響,在單因素條件下,沉積物平均粒徑較大時,NO3-+NO2--N的交換通量較大,20°C時NO3-+NO2--N的交換通量最大,隨DO的升高而降低,且pH值對NO3-+NO2--N的交換通量并無顯著的影響,但(A′B)、(A′C)、(A′D)存在對NO3-+NO2--N交換通量的交互作用,在平均粒徑較大的沉積物類型和低溫度水平共同作用下,NO3-+NO2--N的交換通量較大.在平均粒徑較小的沉積物類型和低DO水平共同作用下, NO3-+NO2--N的交換通量較大.在沉積物為粉砂質(zhì)砂(TS)和pH值為9的共同作用下,NO3-+ NO2--N的交換通量最大.因此,在研究NO3-+ NO2--N的交換通量時,需要對沉積物類型、DO、溫度和pH值的交互作用進行考慮.

    圖5 PO43--P的交互作用

    圖6 SiO32--Si的交互作用

    由表2和圖5可知,只有(B′C)的交互作用對PO43--P的交換通量產(chǎn)生了顯著的影響,在單因素條件下,PO43--P的交換通量隨溫度的上升而增大,隨DO的增大而減小,但(B′C)存在對PO43--P交換通量的交互作用,在低溫和低DO水平的共同作用下,PO43--P的交換通量較小.因此,在研究PO43--P的交換通量時,需要對溫度和DO的交互作用進行考慮.

    由表2和圖6可知,(A′B)和(A′C)的交互作用對SiO32--Si的交換通量產(chǎn)生了顯著的影響,在單因素條件下,沉積物類型和DO濃度并未對SiO32--Si的交換通量產(chǎn)生顯著的影響,而SiO32-- Si的交換通量隨溫度的升高而增大,但 (A′B)和(A′C)存在對SiO32--Si交換通量的交互作用,在黏土質(zhì)粉砂(YT)的沉積物類型和高溫的共同作用下,SiO32--Si的交換通量最大,在黏土質(zhì)粉砂(YT)的沉積物類型和低溶解氧水平的共同作用下,SiO32--Si的交換通量最小.因此,在研究SiO32--Si的交換通量時,需要對沉積物類型和溫度的交互作用及沉積物類型和DO的交互作用進行考慮.

    2.3 正交試驗結(jié)果與往年海州灣實測數(shù)據(jù)比較

    根據(jù)往年海州灣實地采樣,通過實驗室培養(yǎng)得到其交換通量的數(shù)據(jù)(其中2014年部分交換通量數(shù)據(jù)引自高春梅等[57]的研究),并結(jié)合海州灣理化參數(shù)進行分析,如表3所示.各項營養(yǎng)鹽沉積物-水界面的交換通量隨沉積物類型、溫度、DO等的變化趨勢與本文的研究結(jié)果大致相同, DO含量較小的夏季,NH4+-N的交換通量比春季和秋季大.秋季的水溫條件應該更適合NH4+-N的交換,但還是比夏季的交換通量小,其原因可能是夏季雨量充沛,處于泄洪時期,河流攜帶大量營養(yǎng)鹽進入水體.粒徑較小的黏土質(zhì)粉砂NH4+-N平均交換通量也比粒徑較大的粉砂質(zhì)砂大;沉積物為黏土質(zhì)粉砂時,NO3-+NO2--N的交換通量最大. DO含量較小的夏季, NO3-+NO2--N的交換通量比春季和秋季大; DO含量較小的夏季, PO43--P的交換通量比春季和秋季大,且隨著沉積物粒徑的增大而增大.秋季的溫度應更適合PO43--P的交換,但還是比春季的交換通量小,其原因可能是秋季浮游植物豐度較大[58-60],并認為硅藻沉積與沉積物中P 的埋藏有密切聯(lián)系[61],導致其通量減小;SiO32--Si的交換通量隨著夏季溫度升高而增大.

    表3 海州灣往年營養(yǎng)鹽交換通量及理化參數(shù)

    注:表中數(shù)值為正值,表示營養(yǎng)鹽由沉積物向上覆水遷移;負值表示營養(yǎng)鹽由上覆水向沉積物遷移.

    3 結(jié)論

    3.1 影響NH4+-N交換的因素依次為:DO>溫度>沉積物類型;影響NO3-+NO2--N交換的因素依次為:沉積物類型>DO>溫度;影響PO43--P交換的因素依次為: DO>沉積物類型>溫度;影響SiO32--Si交換的因素依次為:溫度>pH值.

    3.2 在分析各因素對營養(yǎng)鹽交換通量的影響和建立沉積物-水界面模型時,應著重考慮其交互作用的影響效果.

    3.3 基于正交實驗所得的沉積物類型、溫度、DO和pH值這4種因素對沉積物-水界面營養(yǎng)鹽交換通量的影響結(jié)果與往年海州灣現(xiàn)場調(diào)查結(jié)果相比較,結(jié)果顯示本研究結(jié)果與往年海州灣現(xiàn)場調(diào)查的營養(yǎng)鹽交換通量變化趨勢具有一致性.

    [1] Reay W G. Sediment-water column oxygen and nutrient fluxes in nearshore environments of the lower Delmarva Peninsula, USA [J]. Marine Ecology Progress, 1995,118(1-3):215-227.

    [2] Callender E, Hammond D E. Nutrient exchange across the sediment-water interface in the Potomac River estuary [J]. Estuarine Coastal & Shelf Science, 1982,15(4):395-413.

    [3] 宋金明.中國近海沉積物-海水界面化學過程與生源物質(zhì)循環(huán)研究[J]. 海洋科學, 2000,24(2):56.

    [4] 董 慧.河口區(qū)沉積物--水界面營養(yǎng)鹽交換通量研究——以李村河為例[D]. 青島:中國海洋大學, 2012.

    [5] 李玲玲.黃河口濕地沉積物中營養(yǎng)鹽分布及交換通量的研究[D]. 青島:中國海洋大學, 2010.

    [6] 宋國棟,劉素美,張國玲.黃東海表層沉積物中磷的分布特征[J]. 環(huán)境科學, 2014,35(1):157-162.

    [7] 李玲玲.黃河口濕地沉積物中營養(yǎng)鹽分布及交換通量的研究[D]. 青島:中國海洋大學, 2010.

    [8] Hopkinson C S. Nutrient regeneration in shallow-water sediments of the estuarine plume region of the nearshore Georgia Bight, USA [J]. Deep Sea Research Part B Oceanographic Literature Review, 1987,94(1):127-142.

    [9] Mortimer R J G, Krom M D, Watson P G, et al. Sediment-Water Exchange of Nutrients in the Intertidal Zone of the Humber Estuary, UK [J]. Marine Pollution Bulletin, 1999,37(3-7):261- 279.

    [10] Nedwell D B, Trimmer M. Nitrogen fluxes through the upper estuary of the Great Ouse, England: the role of the bottom sediments [J]. Marine Ecology Progress, 1998,163(1):109-124.

    [11] Rysgaard S, Risgaard-Petersen N, Peter S N, et al. Oxygen regulation of nitrification and denitrification in sediments [J]. Limnology and Oceanography, 1994,39(7):1643-1652.

    [12] Conley D J, Stockenberg A, Carman R, et al. Sediment-water Nutrient Fluxes in the Gulf of Finland, Baltic Sea [J]. Estuarine Coastal & Shelf Science, 1997,45(5):591-598.

    [13] Rysgaard S. Seasonal variation in nitrification and denitrification in estuarine sediment colonized by benthic microalgae and bioturbating infauna [J]. Marine Ecology Progress Series, 1995, 126(4):154-5.

    [14] Antoniou P, Hamilton J, Koopman B, et al. Effect of temperature and pH on the effective maximum specific growth rate of nitrifying bacteria [J]. Water Research, 1990,24(1):97-101.

    [15] 石 峰.營養(yǎng)鹽在東海沉積物-海水界面交換速率和交換通量的研究[D]. 青島:中國海洋大學, 2003.

    [16] Holdren G C, Armstrong D E. Factors affecting phosphorus release from intact lake sediment cores [J]. Environmental Science & Technology, 1979,14(1):79-87.

    [17] 范成新,相崎守弘.好氧和厭氧條件對霞浦湖沉積物-水界面氮磷交換的影響[J]. 湖泊科學, 1997,9(4):337-342.

    [18] 黃 爽.東海赤潮高發(fā)區(qū)顆粒有機物的來源、分布、分解及其環(huán)境效應[D]. 青島:中國海洋大學, 2012.

    [19] Gieskes J M. Effect of Temperature on the pH of Seawater [J]. Limnology and Oceanography, 1970,14(2):679-685.

    [20] 章守宇,張煥君,焦俊鵬,等.海州灣人工魚礁海域生態(tài)環(huán)境的變化[J]. 水產(chǎn)學報, 2006,30(4):475-480.

    [21] 朱 珠,張 碩.2012年海州灣海洋牧場示范區(qū)水質(zhì)狀況分析與評價[C]. 杭州:中國水產(chǎn)學會學術(shù)年會, 2015.

    [22] 張存勇.連云港近岸海域海洋工程對生態(tài)環(huán)境的影響及其研究[D]. 青島:中國海洋大學, 2006.

    [23] 謝 冕.海州灣南部近岸海域氮、磷營養(yǎng)鹽變化規(guī)律及營養(yǎng)鹽限制狀況[D]. 青島:國家海洋局第一海洋研究所, 2013.

    [24] Zhen S. Phosphorus speciation and effects of environmental factors on release of phosphorus from sediments obtained from Taihu Lake, Tien Lake, and East Lake [J]. Toxicological & Environmental Chemistry, 2015,97(3/4):335-348.

    [25] 張 旭.連云港海州灣近岸海域溶解氧含量及飽和度特征分析[J]. 中國水運:理論版, 2007,5(10):69-70.

    [26] 方南娟,梅肖樂.海州灣近岸海域水質(zhì)現(xiàn)狀與趨勢分析[J]. 水產(chǎn)養(yǎng)殖, 2013,34(1):28-31.

    [27] 周德山.海州灣海域赤潮形成的環(huán)境因子研究[D]. 蘇州:蘇州大學, 2008.

    [28] 徐善良,沈 勤,嚴小軍,等.水樣中氮磷營養(yǎng)鹽的短期保存技術(shù)研究[C]. 福州:中國科學技術(shù)協(xié)會年會, 2010:52-59.

    [29] 趙東波.常用沉積物粒度分類命名方法探討[J]. 海洋地質(zhì)前沿, 2009,25(8):41-44.

    [30] Michaud E, Sundby B, Desrosiers G, et al. The functional group approach to bioturbation: II. The effects of the Macomabalthica, community on fluxes of nutrients and dissolved organic carbon across the sediment-water interface [J]. Journal of Experimental Marine Biology & Ecology, 2006,337(2):178-189.

    [31] Francis C A, Beman J M, Kuypers M M. New processes and players in the nitrogen cycle: the microbial ecology of anaerobic and archaeal ammonia oxidation [J]. Isme Journal, 2007,1(1): 19-27.

    [32] Strous M, Kuenen J G, Jetten M S M. Key Physiology of Anaerobic Ammonium Oxidation [J]. Applied & Environmental Microbiology, 1999,65(7):3248-3250.

    [33] Thamdrup B, Dalsgaard T. Thamdrup B, et al. Production of N2through anaerobic ammonium oxidation coupled to nitrate reduction in marine sediments [J]. Applied & Environmental Microbiology, 2002,68(3):1312-1318.

    [34] Syers J K, Harris R F, Armstrong D E. Phosphate chemistry in lake sediments [J]. Journal of Environmental Quality, 1973,2(1): 1-14.

    [35] Banerjee A, Elefsiniotis P, Tuhtar D. The effect of addition of potato-processing wastewater on the acidogenesis of primary sludge under varied hydraulic retention time and temperature [J]. Journal of Biotechnology, 1999,72(3):203-212.

    [36] 馬 娟,彭永臻,王 麗,等.溫度對反硝化過程的影響以及pH值變化規(guī)律[J]. 中國環(huán)境科學, 2008,28(11):1004-1008.

    [37] 尚會來,彭永臻,張靜蓉,等.溫度對短程硝化反硝化的影響[J]. 環(huán)境科學學報, 2009,29(3):516-520.

    [38] Vlaeminck S E, Terada A, Smets B F, et al. Aggregate Size and Architecture Determine Microbial Activity Balance for One-Stage Partial Nitritation and Anammox [J]. Applied & Environmental Microbiology, 2009,76(3):900-909.

    [39] Volcke E I, Picioreanu C, De B B, et al. Effect of granule size on autotrophic nitrogen removal in a granular sludge reactor [J]. Environmental Technology, 2010,31(11):1271-1280.

    [40] 張 艷,白相東,袁四化,等.顆粒粒徑對三氮硝化-反硝化影響特征研究[J]. 科學技術(shù)與工程, 2015,15(18):243-246.

    [41] Herbert R A. Nitrogen cycling in coastal marine ecosystems [J]. FEMS Microbiology Reviews, 1999,23(5):563-590.

    [42] 戚曉紅,劉素美,張 經(jīng).東、黃海沉積物-水界面營養(yǎng)鹽交換速率的研究[J]. 海洋科學, 2006,30(3):9-15.

    [43] Cerco C F. Measured and modelled effects of temperature, dissolved oxygen and nutrient concentration on sediment-water nutrient exchange [J]. Hydrobiologia, 1989,174(3):185-194.

    [44] 張 輝.黃東海沉積物中營養(yǎng)鹽分布及交換通量研究[D]. 青島:中國海洋大學, 2009.

    [45] Fux C, Boehler M, Huber P, et al. Biological treatment of ammonium-rich wastewater by partial nitritation and subsequent anaerobic ammonium oxidation (anammox) in a pilot plant. [J]. Journal of Biotechnology, 2002,99(3):295-306.

    [46] Jensen H S, Mortensen P B, Rasmussen E, et al. Phosphorus cycling in a coastal marine sediment, Aarhus Bay, Denmark [J]. Limnology and Oceanography, 1995,40(5):908-917.

    [47] Andersen F ?, Ring P. Comparison of phosphorus release from littoral and profundal sediments in a shallow, eutrophic lake [J]. Hydrobiologia, 1999,408-409:175-183.

    [48] Anschutz P, Chaillou G, Lecroart P. Phosphorus diagenesis in sediment of the ThauLagoon [J]. Estuarine Coastal & Shelf Science, 2007,72(3):447-456.

    [49] 趙志梅.渤海灣沉積物磷形態(tài)及營養(yǎng)鹽在沉積物-水界面交換的研究[D]. 咸陽:西北農(nóng)林科技大學, 2005.

    [50] Lawson D S, Hurd D C, Pankratz H S. Silica dissolution rates of decomposing phytoplankton assemblages at various temperatures [J]. American Journal of Science, 1978,278(10):1373-1393.

    [51] Kamatani A. Dissolution rates of silica from diatoms decomposing at various temperatures [J]. Marine Biology, 1982, 68(1):91-96.

    [52] Cappellen P V, Qiu L. Biogenic silica dissolution in sediments of the Southern Ocean. I. Solubility [J]. Deep Sea Research Part II Topical Studies in Oceanography, 1997,44(5):1109-1128.

    [53] Loucaides, Socratis, Cappellen P V, et al. Dissolution of biogenic silica from land to ocean: Role of salinity and pH [J]. Limnology & Oceanography, 2008,53(4):1614-1621.

    [54] Jourabchi P, Meile C, Pasion L R, et al. Quantitative interpretation of pore water O2and pH distributions in deep-sea sediments [J]. Geochimica et Cosmochimica Acta, 2008,72(72):1350-1364.

    [55] Dove P M, Elston S F. Dissolution kinetics of quartz in sodium chloride solutions: Analysis of existing data and a rate model for 25℃[J]. Geochimica et Cosmochimica Acta, 1992,56(12):4147- 4156.

    [56] Loucaides S. Dissolution of biogenic silica: Roles of pH, salinity, pressure, electrical charging and reverse weathering [J]. Journal of Histochemistry & Cytochemistry Official Journal of the Histochemistry Society, 2009,22(12):1092-1104.

    [57] 高春梅,鄭伊汝,張 碩.海州灣海洋牧場沉積物-水界面營養(yǎng)鹽交換通量的研究[J]. 大連海洋大學學報, 2016,31(1):95-102.

    [58] 楊曉改.海州灣及其鄰近海域浮游生物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系[D]. 青島:中國海洋大學, 2015.

    [59] 孔 明,張 路,尹洪斌,等.藍藻暴發(fā)對巢湖表層沉積物氮磷及形態(tài)分布的影響[J]. 中國環(huán)境科學, 2014,34(5):1285-1292.

    [60] 張玉鳳,田 金,楊 爽,等.大連灣海域營養(yǎng)鹽時空分布、結(jié)構(gòu)特征及其生態(tài)響應[J]. 中國環(huán)境科學, 2015,35(1):236-243.

    [61] Diaz J, Ingall E, Benitez-Nelson C, et al. Marine polyphosphate: a key player in geologic phosphorus sequestration [J]. Science, 2008,320(5876):652-655.

    Exchange fluxes of nutrients at the sediment-water interface based on orthogonal experimental design--Taking the marine ranching area of Haizhou Bay as an example.

    ZHANG Shuo1,2, FANG Xin1, HUANG Hong3, ZHANG Hu4, ZHANG Jun-bo1,5,6*

    (1.College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China;2.Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai Ocean University, Shanghai 201306, China;3.College of Marine Ecology and Environment, Shanghai Ocean University, Shanghai 201306, China;4.Maine Fisheries Research Institution of Jiangsu, Nantong 226007, China;5.National Engineering Research Center for Oceanic Fisheries, Shanghai Ocean University, Shanghai 201306, China;6.National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai 201306, China)., 2017,37(11):4266~4276

    The exchange fluxes of nutrients between sediment and water interface were studied in Haizhou Bay, where 3 stations were sampled inMay 2016 to further investigate the characteristics of nutrients fluxes exchange based on orthogonal experiment method. The effects of sediment type, temperature, DO and pH on the exchange fluxes of nutrients between sediment and water interface were analyzed. The primary and secondary relations of various factors above affecting nutrient fluxes could be observed. With respect to NH4+-N, its influencing factors were shown as follows, DO>temperature>sediment type; for NO3-+NO2--N, sediment type>DO>temperature; for PO43--P, DO>sediment type>temperature; for SiO32--Si, temperature>pH. The interaction among these factors was found as an important effect on the fluxes of nutrients, which should be taken into account in the establishment of the sediment-water interface nutrients exchange model. Results in this study on the effect of sediment types, temperatures, dissolved oxygen, and pH on the nutrient exchange fluxes in Haizhou Bay were basically consistent with the investigation data of Haizhou Bay in former years.

    sediment-water interface;nutrinets;exchange flux;orthogonal experiment;interaction

    X131.2

    A

    1000-6923(2017)11-4266-11

    張 碩(1976-),男,天津人,副教授,博士,主要從事海洋生態(tài)修復方面研究.發(fā)表論文40余篇.

    2017-04-07

    海州灣海洋牧場示范項目(D-8006-12-0018,D8006-15-8014);公益性行業(yè)(農(nóng)業(yè))科研專項(201003068,201303047);國家高技術(shù)研究發(fā)展計劃項目(2006A A100303)

    * 責任作者, 副教授, jb_zhang@shou.edu.cn

    猜你喜歡
    粉砂營養(yǎng)鹽沉積物
    晚更新世以來南黃海陸架沉積物源分析
    海洋通報(2022年2期)2022-06-30 06:07:04
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對磷的吸收及釋放研究進展
    典型粉砂地層盾構(gòu)選型及施工參數(shù)研究
    涼水河子河營養(yǎng)鹽滯留能力評估
    躍進總干渠粉砂地基上節(jié)制閘抗滑穩(wěn)定性分析
    瓊東海域冬季、夏季營養(yǎng)鹽結(jié)構(gòu)特征及其對浮游植物生長的影響
    2012年冬季南海西北部營養(yǎng)鹽分布及結(jié)構(gòu)特征
    討論用ICP-AES測定土壤和沉積物時鈦對鈷的干擾
    太湖草源性“湖泛”水域沉積物營養(yǎng)鹽釋放估算
    湖泊科學(2014年2期)2014-02-27 10:29:43
    狂野欧美白嫩少妇大欣赏| 精品一区二区三区视频在线| 超碰97精品在线观看| 亚洲第一av免费看| 最近中文字幕2019免费版| 欧美亚洲 丝袜 人妻 在线| 黑人高潮一二区| 国产成人a∨麻豆精品| 超碰97精品在线观看| 蜜桃在线观看..| 天天操日日干夜夜撸| 色婷婷av一区二区三区视频| 在线观看免费高清a一片| 免费看av在线观看网站| 午夜免费男女啪啪视频观看| 免费高清在线观看视频在线观看| 在线观看国产h片| 国产美女午夜福利| 大又大粗又爽又黄少妇毛片口| 久久久久久久久久人人人人人人| 十分钟在线观看高清视频www | 精品人妻熟女av久视频| 国产成人精品一,二区| 久久国产精品大桥未久av | 高清在线视频一区二区三区| 女性生殖器流出的白浆| 亚洲美女黄色视频免费看| 亚洲精品国产色婷婷电影| 日韩亚洲欧美综合| av在线app专区| 伦理电影大哥的女人| 日本av免费视频播放| 噜噜噜噜噜久久久久久91| 精品人妻熟女av久视频| 国产伦精品一区二区三区四那| 大话2 男鬼变身卡| 各种免费的搞黄视频| 午夜91福利影院| 九九爱精品视频在线观看| 日韩,欧美,国产一区二区三区| 久久99蜜桃精品久久| 久久ye,这里只有精品| 国产精品99久久99久久久不卡 | 欧美xxⅹ黑人| 久久人人爽人人爽人人片va| 精品人妻熟女av久视频| 国产在线男女| 蜜臀久久99精品久久宅男| a级毛片免费高清观看在线播放| av女优亚洲男人天堂| 久久久久久久亚洲中文字幕| 午夜福利网站1000一区二区三区| 欧美少妇被猛烈插入视频| 欧美日韩av久久| 男人添女人高潮全过程视频| a级毛片在线看网站| 一边亲一边摸免费视频| 黑人高潮一二区| 日韩伦理黄色片| 精品一区二区免费观看| 97在线视频观看| 久久久久久久精品精品| 王馨瑶露胸无遮挡在线观看| 亚洲国产最新在线播放| 亚洲av电影在线观看一区二区三区| 国产精品女同一区二区软件| av国产久精品久网站免费入址| 中文乱码字字幕精品一区二区三区| 久久人人爽人人片av| 一级av片app| 97在线人人人人妻| 久久婷婷青草| 国产一区二区三区综合在线观看 | 中文乱码字字幕精品一区二区三区| 99热这里只有是精品50| 我的老师免费观看完整版| 日韩大片免费观看网站| 99久久精品国产国产毛片| 一级爰片在线观看| 人妻夜夜爽99麻豆av| 九色成人免费人妻av| 丝袜脚勾引网站| 99久久精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 一级a做视频免费观看| 久久久久精品性色| 97超碰精品成人国产| 国产极品天堂在线| 免费大片18禁| 18禁动态无遮挡网站| 国产精品.久久久| 亚洲伊人久久精品综合| 我要看黄色一级片免费的| 国产高清三级在线| 亚洲精品aⅴ在线观看| 日本猛色少妇xxxxx猛交久久| 国产淫片久久久久久久久| 久久婷婷青草| 69精品国产乱码久久久| 91精品伊人久久大香线蕉| 99久久人妻综合| 99精国产麻豆久久婷婷| 黑人高潮一二区| 简卡轻食公司| 国产欧美亚洲国产| 精品国产一区二区久久| 国产精品久久久久成人av| 中文字幕久久专区| 男女边摸边吃奶| 久久人人爽av亚洲精品天堂| 久久精品国产a三级三级三级| 最近中文字幕2019免费版| xxx大片免费视频| 亚洲欧美中文字幕日韩二区| 精品少妇黑人巨大在线播放| 丰满人妻一区二区三区视频av| 国产片特级美女逼逼视频| 日韩一区二区视频免费看| 男人舔奶头视频| 美女xxoo啪啪120秒动态图| 18禁动态无遮挡网站| 久久久久久久久久人人人人人人| 久久久精品94久久精品| 热re99久久精品国产66热6| 亚洲欧美日韩东京热| 亚洲av成人精品一二三区| 国语对白做爰xxxⅹ性视频网站| 久久国产精品男人的天堂亚洲 | 午夜福利在线观看免费完整高清在| 国产亚洲欧美精品永久| 麻豆精品久久久久久蜜桃| 精品少妇黑人巨大在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 人人妻人人澡人人看| 一本久久精品| 午夜免费男女啪啪视频观看| 精品久久久久久电影网| 老司机影院成人| 国产精品国产av在线观看| 如日韩欧美国产精品一区二区三区 | 国产欧美亚洲国产| 一级毛片 在线播放| 成年人午夜在线观看视频| 人人妻人人看人人澡| 国产熟女午夜一区二区三区 | av黄色大香蕉| 91午夜精品亚洲一区二区三区| 免费大片18禁| 乱人伦中国视频| 天堂中文最新版在线下载| 国产亚洲最大av| 嘟嘟电影网在线观看| 久久久久久久久久人人人人人人| 国产日韩一区二区三区精品不卡 | 精品久久久噜噜| 久久久午夜欧美精品| 最近中文字幕高清免费大全6| 国产深夜福利视频在线观看| 我的老师免费观看完整版| 18禁在线播放成人免费| 新久久久久国产一级毛片| 久久免费观看电影| 国产高清有码在线观看视频| 高清在线视频一区二区三区| 久久久久久久久久久久大奶| 日本爱情动作片www.在线观看| 老熟女久久久| av在线观看视频网站免费| 久久国产亚洲av麻豆专区| 久久久久久久亚洲中文字幕| 这个男人来自地球电影免费观看 | 国产欧美日韩一区二区三区在线 | 久久久久视频综合| 欧美3d第一页| 交换朋友夫妻互换小说| 日日爽夜夜爽网站| 中文字幕久久专区| 国产熟女午夜一区二区三区 | 国产日韩欧美在线精品| 国产精品成人在线| 国语对白做爰xxxⅹ性视频网站| 人妻 亚洲 视频| 日本午夜av视频| 极品教师在线视频| 亚洲国产精品999| 亚洲欧美精品自产自拍| 免费高清在线观看视频在线观看| 少妇被粗大猛烈的视频| 久久影院123| 精品人妻熟女av久视频| 大又大粗又爽又黄少妇毛片口| 欧美区成人在线视频| 欧美bdsm另类| 视频中文字幕在线观看| 久久精品国产亚洲网站| 国产av一区二区精品久久| 欧美xxⅹ黑人| 中国三级夫妇交换| 久久97久久精品| 久久99热6这里只有精品| 在线观看免费视频网站a站| 日本-黄色视频高清免费观看| 亚洲欧洲国产日韩| 一级毛片我不卡| 蜜桃久久精品国产亚洲av| 亚洲av日韩在线播放| 黄色怎么调成土黄色| 新久久久久国产一级毛片| 这个男人来自地球电影免费观看 | 久久精品国产亚洲av天美| 我的老师免费观看完整版| 在线观看美女被高潮喷水网站| 精品亚洲成a人片在线观看| 亚洲精品视频女| 久久国产亚洲av麻豆专区| 久久狼人影院| 日韩制服骚丝袜av| 99热这里只有是精品在线观看| 看十八女毛片水多多多| xxx大片免费视频| 亚洲成人一二三区av| 99视频精品全部免费 在线| 久久ye,这里只有精品| 亚洲国产精品国产精品| 日韩欧美 国产精品| 亚洲色图综合在线观看| 久久这里有精品视频免费| 3wmmmm亚洲av在线观看| 国产片特级美女逼逼视频| 人妻制服诱惑在线中文字幕| 国产亚洲91精品色在线| 成年av动漫网址| 日韩欧美 国产精品| 中文字幕av电影在线播放| 青春草视频在线免费观看| 欧美日韩精品成人综合77777| 大香蕉97超碰在线| 九草在线视频观看| 国产精品一区二区在线不卡| 久久婷婷青草| 黑人高潮一二区| 中文字幕亚洲精品专区| 日韩欧美 国产精品| 91精品伊人久久大香线蕉| 免费看不卡的av| 日本vs欧美在线观看视频 | 国产成人91sexporn| 免费在线观看成人毛片| 夜夜看夜夜爽夜夜摸| 黄色怎么调成土黄色| 一级毛片aaaaaa免费看小| 少妇人妻 视频| 秋霞在线观看毛片| 国产免费一级a男人的天堂| 黄色毛片三级朝国网站 | 久久韩国三级中文字幕| 97超碰精品成人国产| 亚洲精品色激情综合| 国产欧美日韩一区二区三区在线 | 中国国产av一级| 又黄又爽又刺激的免费视频.| 国产熟女午夜一区二区三区 | 国产成人午夜福利电影在线观看| 在线观看三级黄色| 乱人伦中国视频| 中文欧美无线码| 简卡轻食公司| 三级国产精品欧美在线观看| 精品少妇内射三级| 一个人看视频在线观看www免费| 一级a做视频免费观看| 又爽又黄a免费视频| 中文欧美无线码| 亚洲欧美中文字幕日韩二区| 国产免费福利视频在线观看| 最近手机中文字幕大全| 成人免费观看视频高清| 亚洲国产精品一区三区| a级一级毛片免费在线观看| 99热网站在线观看| 欧美日韩亚洲高清精品| freevideosex欧美| 人妻少妇偷人精品九色| 夜夜看夜夜爽夜夜摸| 中文字幕av电影在线播放| 亚洲精品一二三| 国产国拍精品亚洲av在线观看| 亚洲精品一区蜜桃| 哪个播放器可以免费观看大片| 亚洲人与动物交配视频| 亚洲国产精品一区二区三区在线| 性色avwww在线观看| 九草在线视频观看| 人人妻人人添人人爽欧美一区卜| 亚洲精华国产精华液的使用体验| 欧美另类一区| 久久这里有精品视频免费| 亚洲一级一片aⅴ在线观看| 三上悠亚av全集在线观看 | 精品亚洲乱码少妇综合久久| 亚洲av成人精品一二三区| 成人二区视频| 国产精品成人在线| 国产日韩欧美视频二区| 精品人妻偷拍中文字幕| 五月天丁香电影| 精品国产露脸久久av麻豆| 91在线精品国自产拍蜜月| 国产精品国产av在线观看| av又黄又爽大尺度在线免费看| 欧美成人午夜免费资源| 国产熟女午夜一区二区三区 | 亚洲,欧美,日韩| 欧美少妇被猛烈插入视频| 久久久久久久亚洲中文字幕| 各种免费的搞黄视频| 亚洲国产精品国产精品| 日本wwww免费看| 久热久热在线精品观看| 国产 一区精品| 国产精品国产三级国产专区5o| 交换朋友夫妻互换小说| 成年人午夜在线观看视频| 你懂的网址亚洲精品在线观看| 一二三四中文在线观看免费高清| 九草在线视频观看| 国精品久久久久久国模美| 欧美人与善性xxx| 中国国产av一级| 午夜视频国产福利| 有码 亚洲区| 少妇精品久久久久久久| 99久久中文字幕三级久久日本| 国产爽快片一区二区三区| 亚洲天堂av无毛| 狂野欧美白嫩少妇大欣赏| 我的老师免费观看完整版| 一本大道久久a久久精品| 欧美最新免费一区二区三区| 好男人视频免费观看在线| 五月开心婷婷网| 亚洲精品国产成人久久av| 久久狼人影院| 免费观看a级毛片全部| 亚洲av免费高清在线观看| 免费少妇av软件| 国产成人精品无人区| 卡戴珊不雅视频在线播放| 久久鲁丝午夜福利片| 欧美日韩视频高清一区二区三区二| 蜜桃在线观看..| 丰满乱子伦码专区| 国产av一区二区精品久久| 有码 亚洲区| 免费人成在线观看视频色| 777米奇影视久久| 少妇精品久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 久久影院123| 免费看不卡的av| 亚洲精品国产成人久久av| 成人午夜精彩视频在线观看| 久久ye,这里只有精品| 久久毛片免费看一区二区三区| 97在线视频观看| 人妻 亚洲 视频| 国产精品欧美亚洲77777| 国产伦精品一区二区三区四那| 免费久久久久久久精品成人欧美视频 | 精品午夜福利在线看| 久久久久久久亚洲中文字幕| 久久久国产一区二区| 日韩电影二区| 多毛熟女@视频| av专区在线播放| 国产亚洲午夜精品一区二区久久| 午夜免费观看性视频| 大片电影免费在线观看免费| 久久婷婷青草| 多毛熟女@视频| 另类亚洲欧美激情| 一级二级三级毛片免费看| 国产精品人妻久久久久久| a级毛片免费高清观看在线播放| 少妇的逼水好多| 狂野欧美激情性xxxx在线观看| kizo精华| 亚洲精品久久午夜乱码| 国产精品一区二区三区四区免费观看| 王馨瑶露胸无遮挡在线观看| 免费观看av网站的网址| 下体分泌物呈黄色| 亚洲真实伦在线观看| 在线观看www视频免费| 在线观看免费日韩欧美大片 | 99国产精品免费福利视频| 天美传媒精品一区二区| 大香蕉久久网| 午夜福利影视在线免费观看| 高清黄色对白视频在线免费看 | 哪个播放器可以免费观看大片| 欧美成人精品欧美一级黄| 最黄视频免费看| 国产91av在线免费观看| 99热网站在线观看| 亚洲三级黄色毛片| 噜噜噜噜噜久久久久久91| 午夜福利,免费看| 嫩草影院新地址| 精品酒店卫生间| 国产精品久久久久久av不卡| 在线天堂最新版资源| 日韩亚洲欧美综合| 最新的欧美精品一区二区| 麻豆精品久久久久久蜜桃| 国产 一区精品| 久久久久久人妻| 国产伦精品一区二区三区视频9| 十八禁高潮呻吟视频 | 欧美老熟妇乱子伦牲交| 久久97久久精品| .国产精品久久| 春色校园在线视频观看| 男的添女的下面高潮视频| 久久人人爽人人片av| 深夜a级毛片| 熟妇人妻不卡中文字幕| 我的女老师完整版在线观看| 黄色欧美视频在线观看| 青春草亚洲视频在线观看| 亚洲久久久国产精品| 91久久精品国产一区二区成人| av天堂中文字幕网| 亚洲国产色片| 黄片无遮挡物在线观看| 自线自在国产av| 一本—道久久a久久精品蜜桃钙片| 在线观看一区二区三区激情| 亚洲欧洲日产国产| 亚洲在久久综合| 免费观看性生交大片5| 成人免费观看视频高清| 春色校园在线视频观看| 我的老师免费观看完整版| 国产片特级美女逼逼视频| 纵有疾风起免费观看全集完整版| 国产精品99久久99久久久不卡 | 免费看光身美女| av国产精品久久久久影院| 日本av免费视频播放| 天堂俺去俺来也www色官网| 欧美精品一区二区免费开放| 女的被弄到高潮叫床怎么办| 99热国产这里只有精品6| 欧美成人精品欧美一级黄| 精品国产乱码久久久久久小说| 国产精品久久久久久精品古装| 成人黄色视频免费在线看| 国产爽快片一区二区三区| 亚洲av男天堂| 欧美一级a爱片免费观看看| 国产国拍精品亚洲av在线观看| 在线观看免费视频网站a站| 丰满乱子伦码专区| 不卡视频在线观看欧美| 美女福利国产在线| 肉色欧美久久久久久久蜜桃| 久久久久精品性色| 日产精品乱码卡一卡2卡三| 麻豆精品久久久久久蜜桃| 欧美精品一区二区免费开放| 久久久久久久国产电影| 日本vs欧美在线观看视频 | av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 激情五月婷婷亚洲| 99re6热这里在线精品视频| 啦啦啦视频在线资源免费观看| 国产白丝娇喘喷水9色精品| 欧美日韩国产mv在线观看视频| 欧美三级亚洲精品| 欧美精品亚洲一区二区| 欧美精品一区二区免费开放| 久久久久久久久久久丰满| 亚洲内射少妇av| 亚洲在久久综合| 插阴视频在线观看视频| 91午夜精品亚洲一区二区三区| 久热久热在线精品观看| 亚洲自偷自拍三级| 麻豆成人av视频| 午夜91福利影院| 日韩人妻高清精品专区| 午夜精品国产一区二区电影| 啦啦啦在线观看免费高清www| 亚洲国产精品999| 好男人视频免费观看在线| 欧美精品人与动牲交sv欧美| 亚洲精品456在线播放app| 中国三级夫妇交换| 中文资源天堂在线| 最后的刺客免费高清国语| 男女啪啪激烈高潮av片| 亚洲美女搞黄在线观看| 一级片'在线观看视频| 亚洲久久久国产精品| 一本久久精品| 亚洲精品,欧美精品| 国产乱人偷精品视频| 精品久久久久久久久av| 国产成人精品福利久久| 亚洲国产精品成人久久小说| 国产亚洲5aaaaa淫片| 久久精品国产鲁丝片午夜精品| 欧美国产精品一级二级三级 | 成年人午夜在线观看视频| 国产一区二区三区av在线| 亚洲精品国产av蜜桃| 我的女老师完整版在线观看| 91精品一卡2卡3卡4卡| 国产精品久久久久久av不卡| a级毛片在线看网站| 成年av动漫网址| 色哟哟·www| 五月开心婷婷网| 简卡轻食公司| 亚洲欧洲国产日韩| 午夜视频国产福利| 国产免费一区二区三区四区乱码| 这个男人来自地球电影免费观看 | av线在线观看网站| 99热国产这里只有精品6| 久热久热在线精品观看| 亚洲熟女精品中文字幕| 国产精品99久久99久久久不卡 | 26uuu在线亚洲综合色| 成人美女网站在线观看视频| 大片电影免费在线观看免费| 天堂俺去俺来也www色官网| 久久久a久久爽久久v久久| 少妇人妻久久综合中文| 久久久久久久久久久丰满| 亚洲国产av新网站| 亚洲国产成人一精品久久久| 观看av在线不卡| 久久精品国产亚洲av天美| 久久久久久伊人网av| av在线老鸭窝| 国产欧美另类精品又又久久亚洲欧美| 日韩视频在线欧美| 欧美少妇被猛烈插入视频| 日日摸夜夜添夜夜添av毛片| 日本av手机在线免费观看| 黄色一级大片看看| 成人无遮挡网站| 亚洲电影在线观看av| 下体分泌物呈黄色| 日产精品乱码卡一卡2卡三| 男人狂女人下面高潮的视频| 又大又黄又爽视频免费| 国产在线一区二区三区精| 黑人高潮一二区| 久久久精品免费免费高清| 国产日韩一区二区三区精品不卡 | 国产91av在线免费观看| 丰满乱子伦码专区| 国产av精品麻豆| 国产视频内射| 97超视频在线观看视频| 国产在线视频一区二区| 精品人妻熟女av久视频| 精品一区二区免费观看| 国产毛片在线视频| 99久国产av精品国产电影| 最黄视频免费看| 午夜视频国产福利| 国精品久久久久久国模美| 人体艺术视频欧美日本| 五月天丁香电影| 亚洲,欧美,日韩| 亚洲国产av新网站| 91成人精品电影| 亚洲精品色激情综合| 欧美97在线视频| 国产精品一区二区在线不卡| 成年人免费黄色播放视频 | 人妻人人澡人人爽人人| 午夜福利在线观看免费完整高清在| 最后的刺客免费高清国语| 亚洲精品乱码久久久v下载方式| 最近最新中文字幕免费大全7| 秋霞伦理黄片| 黄色怎么调成土黄色| 成人亚洲欧美一区二区av| 国产精品熟女久久久久浪| 汤姆久久久久久久影院中文字幕| 国产精品蜜桃在线观看| av免费在线看不卡| 在线观看免费视频网站a站| 一级二级三级毛片免费看| av有码第一页| 男人和女人高潮做爰伦理| 91午夜精品亚洲一区二区三区| 男人狂女人下面高潮的视频| 久久久久久伊人网av| 久久久精品免费免费高清| 日本爱情动作片www.在线观看| 熟女电影av网| 亚洲精品第二区| 人人妻人人澡人人看| 日韩电影二区| 最近中文字幕高清免费大全6| 在线天堂最新版资源| 免费黄色在线免费观看| 美女主播在线视频|