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

    垂向湍流擴(kuò)散和光耦合對(duì)下沉藻增長(zhǎng)的影響——基于內(nèi)陸混濁湖泊(太湖)分析

    2019-02-27 09:12:24趙巧華
    中國(guó)環(huán)境科學(xué) 2019年2期
    關(guān)鍵詞:水柱衰減系數(shù)藻類

    周 妍,趙巧華,劉 鵬

    ?

    垂向湍流擴(kuò)散和光耦合對(duì)下沉藻增長(zhǎng)的影響——基于內(nèi)陸混濁湖泊(太湖)分析

    周 妍,趙巧華*,劉 鵬

    (南京信息工程大學(xué)水文與水資源工程學(xué)院,江蘇 南京 210044)

    在假設(shè)溫度恒定?弱化營(yíng)養(yǎng)鹽限制作用的條件下,利用太湖背景漫射衰減系數(shù)?日變化太陽輻射等數(shù)據(jù),通過下沉藻生長(zhǎng)與水環(huán)境相結(jié)合的耦合模型,模擬下沉藻增長(zhǎng)過程中垂向湍流擴(kuò)散和光之間的耦合.結(jié)果表明:在相對(duì)清潔水體中(背景漫射衰減系數(shù)小于1.1/m),下沉藻類無需垂向湍流擴(kuò)散均可維持增長(zhǎng);混濁水體中(背景漫射衰減系數(shù)介于1.1~3.0/m)下沉藻類增長(zhǎng)需垂向湍流擴(kuò)散維持,且最低垂向湍流擴(kuò)散值隨背景漫射衰減系數(shù)增大而增大,二者間存在指數(shù)函數(shù)關(guān)系;最低垂向湍流擴(kuò)散()?水深()與藻類下沉速度()間的佩克萊數(shù)應(yīng)位于0.38~13.89,否則垂向湍流擴(kuò)散對(duì)比其他因素(藻類沉降和光衰減),對(duì)水柱中下沉藻類的增長(zhǎng)的作用甚小;當(dāng)背景漫射衰減系數(shù)大于3.0/m,水柱平均光能可能難以滿足藻類增長(zhǎng),藻類持續(xù)消亡.該研究有助于厘清氣候變化背景下水生生態(tài)系統(tǒng)中的浮游植物種群演替機(jī)制.

    下沉藻;垂向湍流擴(kuò)散;背景漫射衰減系數(shù);數(shù)值模擬

    近幾十年來,由于內(nèi)陸淡水生態(tài)系統(tǒng)水華事件頻發(fā),使水體環(huán)境惡化極大威脅到了淡水系統(tǒng)的水質(zhì)安全[1].理清藻類水華的發(fā)生機(jī)制可以為水生生態(tài)系統(tǒng)的修復(fù)和治理提供相應(yīng)的理論支持.湖泊中藻類的生長(zhǎng)狀態(tài)是多要素綜合反映,主要包括物理?化學(xué)和生物因素.其中營(yíng)養(yǎng)鹽是藻類生長(zhǎng)的物質(zhì)基礎(chǔ)[2-4];湍流擾動(dòng)是影響藻類垂向分布的關(guān)鍵物理量,而光是架構(gòu)物理過程和生物過程的橋梁.近年來,水體的湍流擾動(dòng)對(duì)浮游植物生長(zhǎng)的影響逐步得到重視:秦伯強(qiáng)等[5]經(jīng)過長(zhǎng)期的觀測(cè)得出適當(dāng)擾動(dòng)有利于太湖藍(lán)藻水華的形成.Zhou等[6]也證實(shí)適當(dāng)湍流擾動(dòng)可促進(jìn)微囊藻生長(zhǎng);Huisman等[7]結(jié)合了光照和水動(dòng)力兩個(gè)條件,提出湍流擴(kuò)散和光照對(duì)藻類的共同影響機(jī)制,從理論上分析藻類生長(zhǎng)中擾動(dòng)的貢獻(xiàn);Cao等[8]以水體流速出發(fā)確定微囊藻所需的極限流速;張卓等[9]提出兩個(gè)主要參數(shù)判斷光限制弱水體中藻類生長(zhǎng)情況.

    自然水體中,不同藻類運(yùn)動(dòng)特征差異明顯:靜水中,藻類的密度大于水體密度,藻類呈現(xiàn)下沉運(yùn)動(dòng),反之則為上升運(yùn)動(dòng).這種在靜水中多表現(xiàn)為沉降運(yùn)動(dòng)的藻類稱為下沉型藻類(以下簡(jiǎn)稱下沉藻).微囊藻可通過偽空胞增大浮力,形成上浮運(yùn)動(dòng)[5],尤其在光限制強(qiáng)水體中多上浮[10].吳生才[11]對(duì)太湖底部柵列藻(綠藻)和微囊藻(藍(lán)藻)的沉降分析中指出太湖柵列藻的沉降普遍強(qiáng)于微囊藻.湍流對(duì)藻類垂向分布的影響是湍流與藻類垂向運(yùn)動(dòng)耦合的結(jié)果[7].光能作為藻類生長(zhǎng)的能量來源之一,在垂直水層上逐漸遞減,因而下沉藻須受到湍流的裹挾作用[12-13],得以在水柱中重新分布,使其獲得更多光能以維持其正向增長(zhǎng)[14].太湖中藻類種群以綠藻居多,硅藻和藍(lán)藻緊隨其次[15];夏季太湖多藍(lán)藻[16],風(fēng)浪擾動(dòng)較強(qiáng)的6月湖心和河道口綠藻占比會(huì)有所增高[17].郭文景等[18]通過自回歸分析了營(yíng)養(yǎng)鹽,溫度和葉綠素等水質(zhì)參數(shù)間的定量關(guān)系;趙巧華等[19]證實(shí)春季太湖水溫較低,藻類生長(zhǎng)對(duì)營(yíng)養(yǎng)鹽的需求較低,藻類生長(zhǎng)主要受溫度和光照影響.太湖富營(yíng)養(yǎng)化的程度比較高,當(dāng)營(yíng)養(yǎng)鹽的物質(zhì)條件比較充足,更加凸顯了光照?擾動(dòng)和溫度對(duì)藻類生長(zhǎng)的影響.可見溫度恒定時(shí),光能作為藻類生長(zhǎng)的重要能量來源,與影響浮游植物垂向分布的湍流間的耦合成為影響下沉藻類生長(zhǎng)的關(guān)鍵物理過程.

    太湖系典型的大型淺水渾濁湖泊,風(fēng)浪作用顯著.以往的研究多集中在恒定光源下通過定性?統(tǒng)計(jì)等方法初步探究湍流?光的耦合對(duì)藻類的生長(zhǎng)影響,在維持下沉藻正向增長(zhǎng)過程中,湍流?光能衰減系數(shù)之間有著明顯的耦合關(guān)系,但其鮮有受到關(guān)注.因此,本文基于實(shí)測(cè)日變化光照強(qiáng)度(非恒量的光源供給),假設(shè)溫度恒定(藻類最適生長(zhǎng)溫度),弱化營(yíng)養(yǎng)鹽對(duì)藻類生長(zhǎng)的限制,擬揭示在大型淺水湖泊中(太湖為例)維持下沉藻臨界(藻類不生長(zhǎng)不消亡)狀況下,垂向湍流擴(kuò)散和背景漫射衰減之間的耦合關(guān)系.

    1 數(shù)據(jù)與方法

    1.1 藻類的生長(zhǎng)擴(kuò)散方程

    本文意理清垂向湍流擴(kuò)散和背景漫射衰減衰減對(duì)下沉藻類生長(zhǎng)影響,故忽略營(yíng)養(yǎng)鹽對(duì)藻類生長(zhǎng)的限制,設(shè)定溫度恒定,僅考慮藻類生長(zhǎng)與光照?藻類自身運(yùn)動(dòng)和水體運(yùn)動(dòng)有關(guān).根據(jù)藻類生長(zhǎng)的因素提出藻類生長(zhǎng)對(duì)流擴(kuò)散方程[7]:

    式中:右邊第一項(xiàng)為藻類的生長(zhǎng).為物種類型;為深度,m;為時(shí)間,h;為某種浮游植物的數(shù)密度, cells/m3;藻類的生產(chǎn)率和消耗率的差值分別為凈生產(chǎn)率(,,)=p(,,)-m;因?qū)雎誀I(yíng)養(yǎng)鹽限制和溫度恒定的假設(shè),生產(chǎn)率為光?水深和時(shí)間決定的函數(shù),由莫諾特模型(2)給定:

    式中:消耗率m是關(guān)于溫度的函數(shù)[20],由于溫度恒定故為常數(shù).式(1)右邊第二項(xiàng)為物種在靜水中的移動(dòng).v為藻類的靜水速度(m/h,正為上浮,負(fù)為下沉).式(1)右第三項(xiàng)為垂向擾動(dòng)引起的擴(kuò)散.為垂向湍流擴(kuò)散系數(shù)(m2/h,以下簡(jiǎn)稱垂向擴(kuò)散系數(shù)).max為藻類種的最大生長(zhǎng)速率;H為藻類的半飽和常數(shù);為光量子通量密度(二者單位均為μmol/(m2·s).光量子隨深度的變化為式(3)所示

    式中:in為進(jìn)入水體的光合有效輻射光量子通量密度;k為浮游植物引起的光衰減,m2/cells;K(/m)為懸浮物等產(chǎn)生的背景漫射衰減系數(shù)(以下簡(jiǎn)稱光衰減系數(shù)).為實(shí)驗(yàn)水深,m;各參數(shù)和出處見表1.

    表1 模型參數(shù)設(shè)置

    注:初始數(shù)值.

    綠藻適宜生長(zhǎng)溫度范圍介于20~25℃[23],而表1中各藻類參數(shù)的水體溫度為18~22℃[21],基本符合該結(jié)論范圍.設(shè)定水深2.8m,分辨率為2.8cm.模擬時(shí)長(zhǎng)為15d,采取迎風(fēng)差分的方式,方法見文獻(xiàn)[24].經(jīng)空間差分處理后得剛性微分方程.用全隱式一階剛性方程的Gear算法[25]對(duì)該微分方程求解.擬定積分步長(zhǎng)3.0×10-6h(約等于0.01s).式(1)的邊界條件中,對(duì)流擴(kuò)散項(xiàng)在=0和=時(shí)為0.

    1.2 水體內(nèi)藻類總生物量和總生物量隨時(shí)間變化

    水柱中藻類總生物量的大小可用(4)描述:

    水柱中藻類總生物量的時(shí)間變化如下:

    本文以上述兩變量來表征水柱中藻類生長(zhǎng)的情況.d/d>0且≠0,水柱中藻類存在正向增長(zhǎng);相反水柱中藻類持續(xù)消耗;=0,水柱中無藻不發(fā)生增長(zhǎng).

    1.3 光量子通量密度和水下光衰減系數(shù)

    入射光上?下行輻照度(u和d)數(shù)據(jù)采集于2017年11月1日7:30~16:30,中國(guó)科學(xué)院太湖湖泊生態(tài)系統(tǒng)觀測(cè)研究站,采集頻率為每30min一次.輻射數(shù)據(jù)采集使用德國(guó)TRIOS公司生產(chǎn)的RAMSES水下光譜儀.波長(zhǎng)為的1mol單色輻射的光量子所具有能量為[26]:

    為阿伏加德羅常數(shù),為光速,為普朗克常數(shù).波長(zhǎng)為的單位能量所具有的光量子數(shù):

    取有效輻射波段(400~700nm)的上?下行輻照度和,可得到光合有效輻射光量子通量密度,

    依時(shí)間變化進(jìn)行插值和擬合,得隨時(shí)間變化的入射光有效輻射光量子通量密度.

    據(jù)水下輻照水深度廓線計(jì)算有效波段內(nèi)0~150cm(N=16)水體的衰減系數(shù)表達(dá)式(9),以確定實(shí)驗(yàn)設(shè)定水下光衰減系數(shù)范圍

    對(duì)實(shí)測(cè)數(shù)據(jù)指數(shù)回歸(Origin86回歸),當(dāng)回歸230.95,深度數(shù)N(=16)33接受值,否則視為無效數(shù)據(jù).

    1.4 光量子通量密度日變化函數(shù)

    根據(jù)計(jì)算得到光量子通量密度隨時(shí)間的變化的插值結(jié)果如圖1.從擬合結(jié)果看,擬合曲線基本符合一天內(nèi)光量子通量密度日變化:日照開始于6~7點(diǎn),符合實(shí)測(cè)期真實(shí)日出時(shí)間.正午入射光量子通量密度達(dá)到最大,約858.2μmol/(m2·s),17點(diǎn)光源為0.一天擬合結(jié)果光暗比近似為10:14,實(shí)驗(yàn)光合反應(yīng)時(shí)間小于消耗時(shí)間.故設(shè)定實(shí)驗(yàn)開始時(shí)間為早8點(diǎn),光量子通量密度約300μmol/(m2·s)開始模擬,光變化周期為24h,整個(gè)水體內(nèi)光源供給條件式:

    表2 光衰減系數(shù)隨時(shí)間的變化

    1.5 水下光衰減系數(shù)設(shè)定

    全天內(nèi)水下光衰減系數(shù)的實(shí)測(cè)擬合結(jié)果如表2.結(jié)合本次實(shí)測(cè)結(jié)果和文獻(xiàn)數(shù)據(jù)[27-28],設(shè)定水下漫射衰減系數(shù)范圍為0.5~7.5/m,間隔0.5/m.

    2 結(jié)果與分析

    2.1 光衰減和垂向擴(kuò)散對(duì)不同深度藻類影響

    2.1.1 藻類生物量變化過程中光衰減和垂向擴(kuò)散的耦合作用 垂向擴(kuò)散改變?cè)孱愒谒械姆植?影響藻類可獲取光能的大小,間接影響藻類的生長(zhǎng).圖2為不同深度(15cm?120cm?270cm分別代表表層,中層和近底層)?不同時(shí)間(1?5?10?15d)藻類生物量數(shù)密度隨光衰減系數(shù)和垂向擴(kuò)散系數(shù)的分布.第1d,擴(kuò)散系數(shù)大于0.01m2/h表層和中層生物量數(shù)密度較高;垂向擴(kuò)散系數(shù)為0.001~0.01m2/h時(shí)近底層生物量數(shù)密度較高.不同水深中生物量數(shù)密度隨光衰減系數(shù)增加而減小,生物量數(shù)密度在中層和近底層差異較顯著.第5d,數(shù)密度高值區(qū)域與第1d類似,近表層數(shù)密度隨光衰減系數(shù)變化差異顯著.第10?15d,垂向擴(kuò)散系數(shù)大于0.01m2/h光衰減系數(shù)小于3.0/m時(shí),表層和中層藻類生物量數(shù)密度較高且隨光衰減系數(shù)變化分布差異明顯;垂向擴(kuò)散系數(shù)介于0.001~ 0.01m2/h、光衰減系數(shù)小于1.5/m時(shí),近底層藻類生物量數(shù)密度較高且分布差異明顯.

    圖2 不同深度藻類生物量數(shù)密度隨垂向擴(kuò)散擾動(dòng)和光衰減系數(shù)的分布(列為15?120?270cm水深;行為1?5?10?15d)

    圖3 藻類凈生長(zhǎng)率隨水深變化

    (a)=0.005m2/h; (b)=0.05m2/h; (c)=5.0m2/h; (d)=100.0m2/h=15d

    垂向擴(kuò)散的存在阻礙了藻類沉降.當(dāng)垂向擴(kuò)散系數(shù)越大,表層和中層存留的藻類生物量數(shù)密度越多;相反垂向擴(kuò)散系數(shù)越小,近底層存留的藻類生物量數(shù)密度越多.然而不論何種深度,光衰減系數(shù)越大,藻類在同一深度可獲得的光能越小.清澈環(huán)境中藻類生數(shù)密度遠(yuǎn)高于渾濁環(huán)境.

    2.1.2 藻類凈生產(chǎn)率變化過程中光衰減和垂向擴(kuò)散的耦合 水柱內(nèi)藻類生長(zhǎng)除了與不同深度藻類數(shù)密度分布有關(guān),還取決于藻類的生長(zhǎng)率大小.圖3給出了不同垂向擴(kuò)散條件下,第15d藻類凈生長(zhǎng)率隨水深的變化.光衰減系數(shù)低值區(qū)域,水柱中任意位置藻類保持正向增長(zhǎng)(圖3a).光衰減系數(shù)越大,光能對(duì)藻類凈生長(zhǎng)率限制越明顯,藻類凈生長(zhǎng)率隨深度遞減率增加.垂向擴(kuò)散系數(shù)越強(qiáng),不同光衰減系數(shù)中藻類垂向生長(zhǎng)率曲線變化幅度較小.光衰減系數(shù)通過削弱水柱中的光能大小限制藻類生長(zhǎng),而垂向擴(kuò)散通過裹挾作用改變?cè)孱愒谒械姆植级趸馑p對(duì)藻類生長(zhǎng)的限制程度.垂向擴(kuò)散和光衰減共同影響了藻類水下凈生長(zhǎng)率的分布情況.

    2.2 水柱中藻類總生物量變化過程中垂向擴(kuò)散和光衰減的耦合分析

    為觀察整個(gè)水柱內(nèi)藻類總生物量的變化,本文根據(jù)水柱內(nèi)藻類總生物量()和其時(shí)間變率()分析藻類生長(zhǎng)過程中垂向擴(kuò)散和光衰減的耦合效應(yīng).圖4表示水柱內(nèi)藻類總生物量時(shí)間變率其臨界值(d/ d=0)隨垂向擴(kuò)散?光衰減系數(shù)的變化.藻類總生物量臨界曲線在第3d后達(dá)到穩(wěn)定,4~15d臨界曲線基本一致.光衰減系數(shù)低于1.1/m的水柱內(nèi),藻類生長(zhǎng)可不依賴垂向擴(kuò)散;當(dāng)光衰減系數(shù)逐漸增大,藻類生長(zhǎng)受光限制的作用逐漸凸顯,需要更大的垂向擴(kuò)散裹夾,維持水柱內(nèi)動(dòng)態(tài)生長(zhǎng)的平衡,而所需垂向擴(kuò)散大小隨光衰減系數(shù)增加而增加;當(dāng)光衰減系數(shù)范圍達(dá)到某一臨界值,不論水柱內(nèi)垂向擴(kuò)散如何增大,水柱內(nèi)藻類總生物量持續(xù)消耗.圖5給出了水柱中藻類總生物量隨光衰減和垂向擴(kuò)散的分布.與圖4總生物量生長(zhǎng)臨界曲線一致,藻類總生物量隨著時(shí)間增加逐漸往低光衰減和高垂向擴(kuò)散區(qū)移動(dòng).當(dāng)垂向擴(kuò)散達(dá)到0.1~ 1.0m2/h時(shí),水柱中藻類總生物量較穩(wěn)定.

    結(jié)合圖4和圖5,水柱內(nèi)藻類正向生長(zhǎng)情況可分兩種:1?光衰減系數(shù)小于1.1/m,無論垂向擴(kuò)散大小,整個(gè)水柱內(nèi)藻類總生物量均可正向增加,小擴(kuò)散下總生物量不如大擴(kuò)散;2?光衰減系數(shù)介于1.1~3.0/m,水柱內(nèi)藻類總生物量正向增長(zhǎng)須依賴一定強(qiáng)弱的垂向擴(kuò)散.藻類總生物量生長(zhǎng)所需垂向擴(kuò)散(臨界值)隨光衰減系數(shù)增大而增大(圖4),變化趨勢(shì)為突增,緩慢增加,突增三階段.對(duì)下沉藻,當(dāng)水體較清澈?垂向擴(kuò)散較弱,藻類多沉底但水柱內(nèi)任意位置凈生長(zhǎng)率恒為正,水柱內(nèi)藻類總生物量維持正向增長(zhǎng);水體較清澈但垂向擴(kuò)散較強(qiáng),雖然藻類在下底層凈生長(zhǎng)率為負(fù),但更多的藻類生物量數(shù)密度維持在中上層使得整個(gè)水柱內(nèi)藻類生物量也可維持正向增長(zhǎng).當(dāng)水柱逐漸渾濁,低擴(kuò)散下整個(gè)水柱無法保持凈生長(zhǎng)率恒為正,積累在下底層藻類逐漸消耗,致使整個(gè)水柱內(nèi)藻類消亡,故此時(shí)需要適當(dāng)?shù)拇瓜驍U(kuò)散使水柱內(nèi)更多的藻類停留在真光層獲取光能從而維系水柱中的藻類生長(zhǎng).

    圖4 水柱中下沉藻類總生物量隨時(shí)間的變化曲線(左邊dW/dt30區(qū)域,右邊相反)

    2.3 垂向擴(kuò)散和光衰減在維持藻類生物量生長(zhǎng)中的耦合

    圖4光衰減系數(shù)介于1.1~3.0/m,整個(gè)水柱內(nèi)藻類總生物量正向增長(zhǎng)所需最低垂向擴(kuò)散隨光衰減系數(shù)增大而增大,而光衰減系數(shù)大于3.0/m,無論垂向擴(kuò)散強(qiáng)度如何,水柱內(nèi)藻類總生物量總在消耗.所以水柱內(nèi)藻類正向增長(zhǎng)垂向擴(kuò)散和光衰減二者共同影響出現(xiàn)在光衰減系數(shù)介于1.1~3.0/m間.藻類總生物量生長(zhǎng)所需垂向擴(kuò)散(臨界值)隨光衰減系數(shù)變化的趨勢(shì)非穩(wěn)定緩慢的增長(zhǎng),而存在突增變化.本文用佩克萊數(shù)Pe[29](Pe?1,沉降作用強(qiáng)):

    藻類自身移動(dòng)時(shí)間尺度(/)和擴(kuò)散影響時(shí)間尺度(2/)之比來判定水柱中沉降和擴(kuò)散對(duì)藻類遷移影響的強(qiáng)弱.圖4中擴(kuò)散擾動(dòng)小于0.01m2/h(Pe= 13.89>>1)時(shí),出現(xiàn)第一次突增.Huisman的模擬[21]Pe313.9(圖4虛線)指出過強(qiáng)的沉降(對(duì)流)對(duì)藻類遷移的影響遠(yuǎn)大于垂向擴(kuò)散,導(dǎo)致后者對(duì)藻類裹挾能力弱難以抵消藻類自身沉降的光能獲取劣勢(shì),故而出現(xiàn)垂向擴(kuò)散急增變化.Pe<1(圖4實(shí)線以上)垂向擴(kuò)散開始占主導(dǎo),但淺水湖泊深度有限,垂向擴(kuò)散增強(qiáng)藻類混合越趨于均勻,當(dāng)垂向擴(kuò)散達(dá)到Pe£0.38[21](圖4點(diǎn)劃線),藻類垂向分布均勻且垂向凈生長(zhǎng)率變化不顯著,藻類生長(zhǎng)情況穩(wěn)定,圖5中垂向擴(kuò)散大于0.4m2/h藻類生物量差異變化不大.伴隨光衰減系數(shù)增大,所需最低垂向擴(kuò)散使得水柱內(nèi)藻類分布越均勻,光的過強(qiáng)衰減對(duì)藻類生長(zhǎng)影響更大,最低垂向擴(kuò)散出現(xiàn)第二次突增變化直至與光衰減系數(shù)垂直.

    Huisman在不考慮光衰減系數(shù)大小下給出最低擾動(dòng)的數(shù)值解[24],該數(shù)值解與藻類沉降速度與入射光強(qiáng)度有關(guān).本文在其基礎(chǔ)上引入光衰減系數(shù),選取最低垂向擴(kuò)散隨光衰減系數(shù)緩慢增加的范圍(0.38£P(guān)e£13.89)探討下沉藻臨界態(tài)下最低垂向擴(kuò)散隨光衰減系數(shù)的變化關(guān)系(式12):

    圖5 水柱內(nèi)藻類總生物量隨光衰減和擴(kuò)散擾動(dòng)分布

    色譜柱最低值未標(biāo)注,因低值并不為零

    圖6 最低垂向擴(kuò)散系數(shù)隨光衰減系數(shù)的擬合(其中 A=v2/4g(Iin))

    式(12)中:為藻類的沉降速率,g(in)為平均入射光源下藻類凈生產(chǎn)率.因生物量時(shí)間變化率一般在3d后差異較小(圖4),取3~15d生長(zhǎng)臨界曲線數(shù)據(jù)做最低垂向擴(kuò)散和光衰減系數(shù)間的擬合(Spss 21.0),結(jié)果如圖6.當(dāng)藻類沉降速度越小,停留在真光層時(shí)間越久,此時(shí)藻類所需的最低垂向擴(kuò)散越小;相反,當(dāng)光衰減系數(shù)不變,入射光源越強(qiáng),水柱內(nèi)同一位置藻類可獲取光能越多,藻類生長(zhǎng)對(duì)擴(kuò)散擾動(dòng)的依賴而降低,所需最低擴(kuò)散擾動(dòng)也有所降低.因光衰減系數(shù)的存在,水下光場(chǎng)呈指數(shù)遞減,當(dāng)水下光衰減劇烈,若僅靠藻類自身的遷移,水柱內(nèi)多數(shù)下沉藻類集中于水柱底部(圖2),在底部藻類因獲取的光能有限,呈現(xiàn)消耗狀態(tài)(圖3).越來越多的藻類沉底消耗,整個(gè)水柱的藻類種群將持續(xù)消耗直至消亡.可見,在光限制水體中,水體的垂向擴(kuò)散是下沉藻維持正向增長(zhǎng)的必要條件, 垂向擴(kuò)散的存在促使多數(shù)沉降藻類維持在真光層區(qū)域,以減少因自身沉降的而帶來的種群消耗.而隨光衰減系數(shù)的增大,藻類種群所需借助的外力作用將越來越強(qiáng).

    3 討論

    光能的獲取是影響藻類生長(zhǎng)的重要因素之一.而其他因素諸如溫度[30-31]?營(yíng)養(yǎng)鹽[32]和水動(dòng)力[33-34]也會(huì)影響藻類的生長(zhǎng).其中水動(dòng)力主要是通過調(diào)整藻類在水體中的位置[7]和營(yíng)養(yǎng)物質(zhì)的釋放[35]影響藻類的生長(zhǎng).

    假設(shè)水體溫度恒定不變?弱化營(yíng)養(yǎng)鹽對(duì)藻類的限制,清澈水柱中光衰減緩慢水柱中平均光能可以滿足具有下沉特性的藻類正向增長(zhǎng);當(dāng)水柱光衰減劇烈,下沉藻類需要一定外力作用維持藻類位于真光層區(qū)域.式(3)中,越接近底層光能越小,根據(jù)式(2),藻類凈生長(zhǎng)率較低,在光衰減過強(qiáng)時(shí)呈消耗狀態(tài)(圖3),這并不利于藻類種群的生長(zhǎng).當(dāng)水柱內(nèi)垂向擴(kuò)散較大,表層和中層的藻類生物量數(shù)密度遠(yuǎn)高于垂向擴(kuò)散較狀態(tài)下生物量數(shù)密度(圖2),下沉藻類依賴垂向擴(kuò)散使更多藻類移動(dòng)至真光層,水柱內(nèi)藻類避免持續(xù)消耗,維持了種群的增長(zhǎng)(圖4和圖5).大多數(shù)值模擬[9,36-37]也指出一些以沉降為主的藻類種需要一定大小的垂向擴(kuò)散才能維持藻類的持續(xù)生長(zhǎng).

    垂向擴(kuò)散是否影響藻類的生長(zhǎng)主要取決于藻類對(duì)水下光能的獲取.本文中背景光衰減系數(shù)較小時(shí)藻類永遠(yuǎn)保持正向增長(zhǎng),與垂向擴(kuò)散強(qiáng)弱無關(guān);背景光衰減系數(shù)增大,藻類依賴垂向擴(kuò)散才能維持正向的增長(zhǎng).但當(dāng)背景光衰減系數(shù)超過一定限度,垂向擴(kuò)散并不能抵消光限制對(duì)藻類增長(zhǎng)影響(圖4),所以水柱內(nèi)藻類總生物量生長(zhǎng)是否依賴于垂向擴(kuò)散與背景光衰減系數(shù)有關(guān).以往的研究多集中一定大小的垂向擴(kuò)散和光衰減是否會(huì)引起藻類生長(zhǎng)的定性研究,缺乏二者之間的定量分析.Huisman等[24]雖給出沉降藻類所需最低垂向擴(kuò)散的數(shù)值表達(dá)式,但該表達(dá)式并未考慮背景光衰減系數(shù)的影響,圖4和圖6體現(xiàn)當(dāng)光衰減系數(shù)增加,藻類所需最低垂向擴(kuò)散會(huì)增大,所以最低垂向擴(kuò)散應(yīng)是藻類速率、入射光強(qiáng)度和背景光衰減系數(shù)共同決定.

    但垂向擴(kuò)散是沒有方向性的, 垂向擴(kuò)散過強(qiáng),延緩下沉藻沉降同時(shí)也會(huì)帶動(dòng)一定數(shù)量的藻類進(jìn)入水柱無光區(qū),這樣反而不利于整個(gè)水柱內(nèi)藻類生長(zhǎng)[7].Huisman等[21]和Portalier等[37]均證實(shí)下沉藻類的正向生長(zhǎng)須存在一個(gè)上限和一個(gè)下限擴(kuò)散值.但本文僅存一個(gè)下限擴(kuò)散值,與二者結(jié)論略不同.在二者研究中藻類移動(dòng)速率較小,水深設(shè)定較大,藻類沉降能力并不凸顯;鑒于本文研究對(duì)象區(qū)域是大型淺水湖泊,水深較淺,從而凸現(xiàn)了藻類沉降效應(yīng).Huisman等[24]也驗(yàn)證當(dāng)藻類移動(dòng)速率相對(duì)于水深較快,上限擴(kuò)散值減小下限擴(kuò)散值增加,導(dǎo)致兩個(gè)擴(kuò)散限度相融合,最后僅存在一條臨界曲線(圖4).所以淺水區(qū)域,沉降能力強(qiáng)的藻類僅存在一個(gè)下限擴(kuò)散值,即可確保整個(gè)水柱內(nèi)藻類種群的正向增長(zhǎng).

    4 結(jié)論

    4.1 在水深2.8m,日變化供給光源光暗比為10:14條件下,水中背景光衰減系數(shù)小于1.1/m,水中光限制作用較小,下沉藻類從表層至底層可保持非消耗狀態(tài),水柱內(nèi)藻類不需要垂向擴(kuò)散可保持種群的正向增長(zhǎng);水中背景光衰減系數(shù)介于1.1~3.0/m時(shí),水柱內(nèi)藻類種群的正向增長(zhǎng)須借助一定大小的垂向擴(kuò)散維持,且垂向擴(kuò)散強(qiáng)弱與光衰減大小有關(guān);水中背景光衰減系數(shù)大于3.0/m,水中光限制過強(qiáng),即便存在高強(qiáng)度的垂向擴(kuò)散,水柱內(nèi)下沉藻類種群依舊無法維持正向增長(zhǎng),生物量逐漸消亡.

    4.2 在不考慮營(yíng)養(yǎng)鹽限制和假設(shè)溫度恒定條件下,下沉藻類生長(zhǎng)主要受生長(zhǎng)、對(duì)流和擴(kuò)散三方面影響.所需最低垂向擴(kuò)散強(qiáng)度應(yīng)與藻類沉降速率,入射光源大小和水下光衰減三方面因素有關(guān).

    即最低垂向擴(kuò)散應(yīng)與藻類沉降速率正向相關(guān),藻類沉降越慢,留在真光層的時(shí)間越久,不需強(qiáng)擴(kuò)散即可正向增長(zhǎng);當(dāng)入射光小,水下藻類獲取光能有限,需要強(qiáng)擴(kuò)散使藻類移動(dòng)至光能較大區(qū)域使藻類得以正向增長(zhǎng);同理水下背景漫射光衰減系數(shù)與擴(kuò)散正向相關(guān).但該擬合結(jié)果須滿足0.38£P(guān)e=/£13.89(代表水深),否則垂向擴(kuò)散對(duì)藻類增長(zhǎng)的影響較小.前者因淺水湖泊中,過高的擴(kuò)散使水柱內(nèi)藻類分布均勻,決定水柱內(nèi)藻類增長(zhǎng)的主要因素為光的衰減;后者則因藻類自身沉降過強(qiáng),過低的擴(kuò)散無法抵消藻類自身的沉降消耗所致.

    [1] 秦伯強(qiáng),王小冬,湯祥明,等.太湖富營(yíng)養(yǎng)化與藍(lán)藻水華引起的飲用水危機(jī)——原因與對(duì)策 [J]. 地球科學(xué)進(jìn)展, 2007,22(9):896-906. Qin B, Wang X, Tang X, et al. Drinking water crisis caused by eutrophication and cyanobacterial bloom in Lake Taihu: cause and measurement [J]. Advances in Earth Science, 2007,22(9):896-906.

    [2] Chen M, Fan M, Liu R, et al. The dynamics of temperature and light on the growth of phytoplankton [J]. Journal of Theoretical Biology, 2015,385:8-19.

    [3] Zhang M, Duan H, Shi X, et al. Contributions of meteorology to the phenology of cyanobacterial blooms: implications for future climate change [J]. Water Research, 2012,46(2):442-452.

    [4] Huisman J, Weissing F J. Light-Limited Growth and Competition for Light in Well-Mixed Aquatic Environments: An Elementary Model [J]. Ecology, 1994,75(2):507-520.

    [5] 秦伯強(qiáng),楊桂軍,馬健榮,等.太湖藍(lán)藻水華“暴發(fā)”的動(dòng)態(tài)特征及其機(jī)制 [J]. 科學(xué)通報(bào), 2016,61(7):759-770. Qin B, Yang G, Ma J, et al. Dynamics of variability and mechanism of harmful cyanobacteria bloom in Lake Taihu, China [J]. Chinese Science Bulletin2016,61(7):759-770.

    [6] Zhou J, Qin B, Han X, et al. Turbulence increases the risk of microcystin exposure in a eutrophic lake (Lake Taihu) during cyanobacterial bloom periods [J]. Harmful Algae, 2016,55:213-220.

    [7] Huisman J, Van Oostveen P, Weissing F J. Critical depth and critical turbulence: Two different mechanisms for the development of phytoplankton blooms [J]. Limnology & Oceanography, 1999,44(7): 1781-1787.

    [8] Cao H, Kong F, Luo L, et al. Effects of Wind and Wind-Induced Waves on Vertical Phytoplankton Distribution and Surface Blooms of Microcystis aeruginosa in Lake Taihu [J]. Journal of Freshwater Ecology, 2006,21(2):231-238.

    [9] 張 卓,宋志堯,黃昌春,等.水藻暴發(fā)的影響因素定量化研究初步 [J]. 環(huán)境科學(xué), 2013,34(7):2603-2610. Zhang Z, Song Z, Huang C, et al. Elementary Quantitative Study on Factors of Phytoplankton Bloom [J]. Environmental Science, 2013, 34(7):2603-2610.

    [10] 俞 茜,陳永燦,劉昭偉.靜止水體中微囊藻屬遷移軌跡的數(shù)值模擬 [J]. 中國(guó)環(huán)境科學(xué), 2017,37(5):1915-1921. Yu Q, Chen Y, Liu Z. Modelling vertical migration trajectory of Microcystis in calm water [J]. China Environmental Science, 2017, 37(5):1915-1921.

    [11] 吳生才.太湖水華藻類越冬生態(tài)機(jī)制和浮游植物的研究 [D]. 南京:中國(guó)科學(xué)院南京地理與湖泊研究所, 2004. Wu S. Study on ecological mechanism of bloom alga overwintering and phytoplankton in Lake Taihu [D]. Nanjing: Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, 2004.

    [12] Macias D, Santana A R, Romero E R, et al. Turbulence as a driver for vertical plankton distribution in the subsurface upper ocean [J]. Scientia Marina, 2013,77(4):541-549.

    [13] 吳挺峰,朱廣偉,秦伯強(qiáng),等.前期風(fēng)場(chǎng)控制的太湖北部湖灣水動(dòng)力及對(duì)藍(lán)藻水華影響 [J]. 湖泊科學(xué), 2012,24(3):409-415. Wu T, Zhu G, Qin B, et al. Prior wind field induced hydrodynamics and its influence on cyanobacterial bloom in northern bays of Lake Taihu, China [J]. Journal of Lake Sciences, 2012,24(3):409-415.

    [14] 俞 茜,陳永燦,劉昭偉,等.藻類光競(jìng)爭(zhēng)模型構(gòu)建及水體紊動(dòng)對(duì)競(jìng)爭(zhēng)的影響 [J]. 中國(guó)環(huán)境科學(xué), 2018,38(7):2665-2670. Yu Q, Chen Y, Liu Z , et al. Development of competition model of algae for light and the impact of turbulence on algal competition [J]. China Environmental Science, 2018,38(7):2665-2670.

    [15] 袁信芳,施華宏,王曉蓉.太湖著生藻類的時(shí)空分布特征 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2006,25(4):1035-1040. Yuan X, Shi H, Wang X. Temporaland spatial distributionsof periphytic algae in Taihu Lake [J]. Journal of Agro-Environment Science, 2006,25(4):1035-1040.

    [16] 虞銳鵬,何恩奇,鈕偉民,等.2010年度太湖水域藻密度和葉綠素時(shí)空分布特征的研究 [C]. 2011中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第一卷): 367-371. Yu R, He E, Niu W, et al. Study on the characteristics of algal density andspatial and temporal distribution in Lake Taihu in 2010 [C]. Proceedings of the 2011 annual conference of the Chinese Society of Environmental Sciences (Volume I):367-371.

    [17] 邢 鵬,孔繁翔,曹煥生,等.太湖浮游細(xì)菌與春末浮游藻類群落結(jié)構(gòu)演替的相關(guān)分析 [J]. 生態(tài)學(xué)報(bào), 2007,27(5):1696-1702. Xing P, Kong F, Cao H, et al. Relationship between bacterioplankton and phytoplankton community dynamics during late spring and early summer in Lake Taihu [J]. Acta Ecologica Sinica, 2007,27(5):1696- 1702.

    [18] 郭文景,符志友,汪 浩,等.水華過程水質(zhì)參數(shù)與浮游植物定量關(guān)系的研究--以太湖梅梁灣為例 [J]. 中國(guó)環(huán)境科學(xué), 2018,38(4): 1517-1525. Guo W, Fu Z, Wang H, et al. The quantitative relation of aquatic parameters and phytoplankton biomass in the process of algal blooms—the case of Meiliang Bay in Taihu Lake [J]. China Environmental Science, 2018,38(4):1517-1525.

    [19] 趙巧華,孫國(guó)棟,王健健,等.水溫?光能對(duì)春季太湖藻類生長(zhǎng)的耦合影響 [J]. 湖泊科學(xué), 2018,30(2):385-393. Zhao Q, Sun G, Wang J, et al. Coupling effect of water temperature and light energy on the algal growth in Lake Taihu [J]. Journal of Lake Sciences, 2018,30(2):385-393.

    [20] 許秋瑾,秦伯強(qiáng),陳偉民,等.太湖藻類生長(zhǎng)模型研究 [J]. 湖泊科學(xué), 2001,13(2):149-157. Xu Q, Qin B, Chen W, et al. Ecological simulation of algae growth in Taihu Lake [J]. Journal of Lake Science, 2001,13(2):149-157.

    [21] Huisman J, Sharples J, Stroom J M, et al. Changes in turbulent mixing shift competition for light between phytoplankton species [J]. Ecology, 2004,85(11):2960-2970.

    [22] 李春華,葉 春,張 詠,等.太湖湖濱帶藻密度與水質(zhì)?風(fēng)作用的分布特征及相關(guān)關(guān)系 [J]. 環(huán)境科學(xué)研究, 2013,26(12):1290-1300. Li C, Ye C, Zhang Y, et al. Temporal and spatial distribution of algal density and its relationship with water quality and wind factor in the littoral zone of Lake Taihu [J]. Research of Environmental Sciences, 2013,26(12):1290-1300.

    [23] Liu X, Lu X, Chen Y. The effects of temperature and nutrient ratios on Microcystis blooms in Lake Taihu, China: An 11-year investigation [J]. Harmful Algae, 2011,10(3):337-343.

    [24] Huisman J, Arrayás M, Ebert U, et al. How do sinking phytoplankton species manage to persist? [J]. American Naturalist, 2002,159(3):245- 254.

    [25] 徐士良. FORTRAN常用算法程序集 [M]. 清華大學(xué)出版社, 1992:321-331. Xu S. FORTRAN common algorithm assembly [M]. Tsinghua University Press, 1992:321-331.

    [26] 周允華,項(xiàng)月琴.光合有效量子通量密度的氣候?qū)W計(jì)算 [J]. 氣象學(xué)報(bào), 1996,4:447-455. Zhou Y, Xiang Y. Climatological estimaton of quantum flux densities [J]. Acta Meteorologica Sinica, 1996,4:447-455.

    [27] 張運(yùn)林,秦伯強(qiáng),陳偉民,等.太湖梅梁灣沿岸帶水體生物學(xué)與光學(xué)特性 [J]. 生態(tài)學(xué)報(bào), 2005,25(3):454-460. Zhang Y, Qin B, Chen W, et al. Biological and optical properties of Meiliang Bay in Lake Taihu [J]. Acta Ecologica Sinica, 2005,25(3): 454-460.

    [28] 張運(yùn)林,秦伯強(qiáng),陳偉民,等.不同風(fēng)浪條件下太湖梅梁灣光合有效輻射的衰減 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2005,16(6):1133-1137. Zhang Y, Qin B, Chen W, et al. Attenuation of photosynthetically available radiation( PAR)in Meiliang Bay under different winds and waves [J]. Chinese Journal of Applied Ecology, 2005,16(6):1133- 1137.

    [29] Visser P M, Ibelings B W, Bormans M, et al. Artificial mixing to control cyanobacterial blooms: a review [J]. Aquatic Ecology, 2016, 50(3):423-441.

    [30] 吳 攀,鄧建明,秦伯強(qiáng),等.水溫和營(yíng)養(yǎng)鹽增加對(duì)太湖冬?春季節(jié)藻類生長(zhǎng)的影響 [J]. 環(huán)境科學(xué)研究, 2013,26(10):1064-1071. Wu P, Deng J, Qin B, et al. Effects of enhanced water temperature and nutrient concentration on algal growth in winter and spring season in Lake Taihu, China [J]. Research of Environmental Sciences, 2013, 26(10):1064-1071.

    [31] 譚 嘯,孔繁翔,于 洋,等.升溫過程對(duì)藻類復(fù)蘇和群落演替的影響 [J]. 中國(guó)環(huán)境科學(xué), 2009,29(6):578-582. Tan X, Kong F, Yu Y, et al. Effects of enhanced temperature on algae recruitment and phytoplankton community succession [J]. China Environmental Science, 2009,29(6):578-582.

    [32] 吳雅麗,許 海,楊桂軍,等.太湖春季藻類生長(zhǎng)的磷營(yíng)養(yǎng)鹽閾值研究 [J]. 中國(guó)環(huán)境科學(xué), 2013,33(9):1622-1629. Wu Y, Xu H, Yang G, et al. Developing the critical phosphorus threshold for spring algal growth in Lake Taihu, China [J]. China Environmental Science, 2013,33(9):1622-1629.

    [33] 陳偉民,陳宇煒,秦伯強(qiáng),等.模擬水動(dòng)力對(duì)湖泊生物群落演替的實(shí)驗(yàn) [J]. 湖泊科學(xué), 2000,12(4):343-352. Chen W, Chen Y, Qin B, et al. Experimental Study on the Biological Community Succession Caused by Water Flow [J]. Journal of Lake Science, 2000,12(4):343-352.

    [34] Huisman J, Codd G A, Paerl H W,et al. Cyanobacterial blooms [J]. Nature Reviews Microbiology, 2018,16(8):471-483.

    [35] 朱廣偉,秦伯強(qiáng),高 光.強(qiáng)弱風(fēng)浪擾動(dòng)下太湖的營(yíng)養(yǎng)鹽垂向分布特征 [J]. 水科學(xué)進(jìn)展, 2004,15(6):775-780. Zhu G, Qin B, Gao G, et al. Vertical distribution of the concentrations of phosphorus and suspended solid in Taihu lake affected by wind-induced wave [J]. Advances in Water Science, 2004,15(6): 775-780.

    [36] Bengfort M, Malchow H. Vertical mixing and hysteresis in the competition of buoyant and non-buoyant plankton prey species in a shallow lake [J]. Ecological Modelling, 2016,323:51-60.

    [37] Portalier S M J, Cherif M, Zhang L,et al. Size-related effects of physical factors on phytoplankton communities [J]. Ecological Modelling, 2016,323:41-50.

    致謝:本文的英文摘要潤(rùn)色工作由唐德才老師協(xié)助完成,在此表示感謝.

    Effect of vertical turbulent diffusion and light coupling on the growth of sinking algae——Analysis based on inland turbid lakes (Lake Taihu).

    ZHOU Yan, ZHAO Qiaohua*, LIU Peng

    (School of Hydrology and Water Resources, Nanjing University of Information Science and Technology, Nanjing 210044, China)., 2019,39(2):792~801

    Based on the assumption of constant temperature and no nutrient salt limitation, a mathematical model coupling sinking algae growth and water environment of Lake Taihu was used to investigate the effect between the vertical turbulent diffusion and background diffused attenuation coefficient(K) on the growth of sinking algae. The simulation results showed that the total biomass growth of the sinking algae could be sustained without vertical turbulent diffusion in the relatively clean water (Kwas less than 1.1/m); in turbid water (Kwas between 1.1~3.0/m), the growth required vertical turbulent diffusion to maintain, and the minimum vertical turbulence diffusion increased with theK. The Peclet number between the minimum vertical turbulence diffusion(), water depth() and the algae sinking speed() should between 0.38~13.89, otherwise the minimum vertical turbulence would have little effect on the growth than other factors (for example, algae sinking speed or light attenuation). When theKwas greater than 3.0/m, the light intensity of the water column was too low to maintain the total biomass growth of algae. This paper helped to clarify the phytoplankton population succession mechanism in aquatic ecosystems under climate change.

    sinking algae;vertical turbulence diffusion;background diffused attenuation coefficient;numerical simulation

    X524

    A

    1000-6923(2019)02-0792-10

    周 妍(1994-),女,安徽安慶人,南京信息工程大學(xué)碩士研究生,主要從事湖泊動(dòng)力與藻類生長(zhǎng)研究.

    2018-07-04

    國(guó)家自然科學(xué)項(xiàng)目(41371222)

    * 責(zé)任作者, 教授, qhzhao@nuist.edu.cn

    猜你喜歡
    水柱衰減系數(shù)藻類
    藻類水華控制技術(shù)及應(yīng)用
    探探鯨的水柱
    細(xì)菌和藻類先移民火星
    軍事文摘(2020年20期)2020-11-16 00:31:40
    Run through the rain
    吃蔬菜有個(gè)“321模式” 三兩葉菜類,二兩其他類,一兩菌藻類
    復(fù)合材料孔隙率的超聲檢測(cè)衰減系數(shù)影響因素
    近岸及內(nèi)陸二類水體漫衰減系數(shù)的遙感反演研究進(jìn)展
    對(duì)《電磁波衰減系數(shù)特性分析》結(jié)果的猜想
    水柱有“魔力”
    浮游藻類在水體PAHs富集與降解中的研究
    天美传媒精品一区二区| 国内精品美女久久久久久| 18禁黄网站禁片免费观看直播| 亚洲国产欧洲综合997久久,| 一个人免费在线观看电影| 麻豆成人午夜福利视频| 成人欧美大片| 一区二区三区高清视频在线| 亚洲av美国av| 久久久国产成人免费| 老女人水多毛片| 国语自产精品视频在线第100页| 国产真实伦视频高清在线观看| 亚洲av美国av| 午夜精品在线福利| 国产成人a区在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲成人中文字幕在线播放| 综合色丁香网| 91精品国产九色| av专区在线播放| 国产精品久久久久久久电影| 美女黄网站色视频| 亚洲欧美清纯卡通| 国产色婷婷99| 久久久久久九九精品二区国产| 亚洲精品色激情综合| 日韩高清综合在线| 天堂影院成人在线观看| 亚洲成人av在线免费| 日韩精品有码人妻一区| 一边摸一边抽搐一进一小说| av黄色大香蕉| 永久网站在线| 有码 亚洲区| 最近在线观看免费完整版| 国产在线精品亚洲第一网站| 亚洲欧美精品自产自拍| 国国产精品蜜臀av免费| 国产视频内射| 99久久九九国产精品国产免费| 少妇人妻一区二区三区视频| 日本a在线网址| 女人被狂操c到高潮| 国产午夜精品论理片| 人人妻人人澡欧美一区二区| 免费大片18禁| 蜜桃久久精品国产亚洲av| 黄色一级大片看看| 菩萨蛮人人尽说江南好唐韦庄 | 午夜视频国产福利| 97碰自拍视频| 色综合站精品国产| 干丝袜人妻中文字幕| 色5月婷婷丁香| 自拍偷自拍亚洲精品老妇| 狂野欧美激情性xxxx在线观看| 熟女人妻精品中文字幕| 美女被艹到高潮喷水动态| 午夜福利高清视频| 亚洲性夜色夜夜综合| 成年免费大片在线观看| 日韩成人伦理影院| 美女大奶头视频| 嫩草影院新地址| 国产精品伦人一区二区| 亚洲国产精品合色在线| 男女做爰动态图高潮gif福利片| 精品久久久久久久末码| 精品人妻一区二区三区麻豆 | 一级毛片久久久久久久久女| 精品日产1卡2卡| 精品少妇黑人巨大在线播放 | 国产精品一区二区性色av| 亚洲自偷自拍三级| 亚洲成人久久爱视频| 麻豆一二三区av精品| 亚洲av中文av极速乱| a级毛片a级免费在线| 中文字幕精品亚洲无线码一区| 国产精品野战在线观看| 国产精品无大码| 亚洲七黄色美女视频| 色av中文字幕| 国产一级毛片七仙女欲春2| 99在线人妻在线中文字幕| 久久99热这里只有精品18| 一级毛片电影观看 | 我要看日韩黄色一级片| 校园春色视频在线观看| 身体一侧抽搐| 免费看美女性在线毛片视频| 成年av动漫网址| 国产高清不卡午夜福利| 有码 亚洲区| 精品无人区乱码1区二区| 中国美白少妇内射xxxbb| 一级a爱片免费观看的视频| 日韩欧美免费精品| 特大巨黑吊av在线直播| 高清午夜精品一区二区三区 | 亚洲国产色片| 中文资源天堂在线| 男女做爰动态图高潮gif福利片| 日韩欧美精品免费久久| 午夜免费男女啪啪视频观看 | 午夜久久久久精精品| 夜夜看夜夜爽夜夜摸| 国产男人的电影天堂91| 国产人妻一区二区三区在| 国产伦精品一区二区三区四那| 国产三级中文精品| 亚洲av熟女| 高清毛片免费看| 国产视频一区二区在线看| 亚州av有码| 成人av在线播放网站| 看黄色毛片网站| 亚洲乱码一区二区免费版| 精品免费久久久久久久清纯| 在线播放国产精品三级| 午夜福利在线观看免费完整高清在 | 国产一区二区在线av高清观看| 少妇人妻精品综合一区二区 | 国产精品亚洲一级av第二区| 九色成人免费人妻av| 久久久久性生活片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av美国av| 女人被狂操c到高潮| 网址你懂的国产日韩在线| 国产精品,欧美在线| 1024手机看黄色片| 免费av观看视频| 亚洲国产精品sss在线观看| 91久久精品国产一区二区成人| 久久久久久久久久久丰满| 亚洲中文字幕一区二区三区有码在线看| 亚洲国产欧美人成| 色av中文字幕| 大型黄色视频在线免费观看| 日日摸夜夜添夜夜爱| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区三区人妻视频| 露出奶头的视频| 女的被弄到高潮叫床怎么办| 成年女人看的毛片在线观看| 舔av片在线| 亚洲最大成人av| 大又大粗又爽又黄少妇毛片口| 色综合色国产| 九色成人免费人妻av| 免费在线观看成人毛片| 亚洲成人av在线免费| 国产一区二区在线av高清观看| 国产亚洲欧美98| 免费无遮挡裸体视频| 久久久久久久久大av| 国产人妻一区二区三区在| 搞女人的毛片| 床上黄色一级片| 性欧美人与动物交配| 在线观看66精品国产| 天堂√8在线中文| 亚洲成a人片在线一区二区| 不卡视频在线观看欧美| 国产成人福利小说| 国产精品久久电影中文字幕| 午夜福利高清视频| 欧美在线一区亚洲| 麻豆一二三区av精品| 日本色播在线视频| 亚洲国产精品成人久久小说 | 国产高清视频在线播放一区| 日日干狠狠操夜夜爽| 亚洲人成网站高清观看| 色噜噜av男人的天堂激情| 久久九九热精品免费| 91在线观看av| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲精品av在线| 亚洲欧美日韩无卡精品| 国产亚洲av嫩草精品影院| 搡老熟女国产l中国老女人| 婷婷六月久久综合丁香| 可以在线观看的亚洲视频| 成人高潮视频无遮挡免费网站| 我的女老师完整版在线观看| a级毛色黄片| av天堂中文字幕网| 深爱激情五月婷婷| 丝袜喷水一区| 亚洲精品一区av在线观看| 亚洲精品影视一区二区三区av| 国产精品电影一区二区三区| 成人国产麻豆网| 成人性生交大片免费视频hd| 三级经典国产精品| 欧美性感艳星| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 成人av在线播放网站| 久久久久国产精品人妻aⅴ院| 国产精品一区二区三区四区久久| 最近手机中文字幕大全| 天天一区二区日本电影三级| 国产一区二区亚洲精品在线观看| 晚上一个人看的免费电影| 国产午夜福利久久久久久| 偷拍熟女少妇极品色| 国产在视频线在精品| 欧美+亚洲+日韩+国产| 性插视频无遮挡在线免费观看| 99精品在免费线老司机午夜| 久久精品国产亚洲av香蕉五月| 人妻少妇偷人精品九色| 九九在线视频观看精品| 舔av片在线| 成熟少妇高潮喷水视频| 国产av在哪里看| 18禁裸乳无遮挡免费网站照片| 欧美最黄视频在线播放免费| 热99在线观看视频| 国产美女午夜福利| 97超碰精品成人国产| 少妇人妻一区二区三区视频| 日韩高清综合在线| 国产色婷婷99| 99国产极品粉嫩在线观看| 六月丁香七月| 欧美激情国产日韩精品一区| 欧美精品国产亚洲| 免费看日本二区| 国产亚洲精品久久久久久毛片| 中国国产av一级| 国产乱人偷精品视频| 97碰自拍视频| 日韩成人av中文字幕在线观看 | 国国产精品蜜臀av免费| 国产综合懂色| 久久精品夜夜夜夜夜久久蜜豆| 别揉我奶头~嗯~啊~动态视频| 伊人久久精品亚洲午夜| a级毛片免费高清观看在线播放| 亚洲美女搞黄在线观看 | 精品日产1卡2卡| 日韩欧美 国产精品| 夜夜爽天天搞| 特级一级黄色大片| 国国产精品蜜臀av免费| 欧美一级a爱片免费观看看| 亚洲国产精品sss在线观看| 1024手机看黄色片| 在线天堂最新版资源| 观看美女的网站| 午夜精品一区二区三区免费看| 国产一区二区三区av在线 | 国产色爽女视频免费观看| 少妇人妻一区二区三区视频| 午夜激情欧美在线| 国产女主播在线喷水免费视频网站 | 真人做人爱边吃奶动态| 国产精品国产三级国产av玫瑰| 成人二区视频| 国产av在哪里看| 男人和女人高潮做爰伦理| 国产精品乱码一区二三区的特点| 亚洲精品国产成人久久av| 老师上课跳d突然被开到最大视频| 免费看a级黄色片| 午夜精品在线福利| 国产视频一区二区在线看| 搡老妇女老女人老熟妇| 午夜日韩欧美国产| 伊人久久精品亚洲午夜| 18禁在线播放成人免费| 亚洲av中文字字幕乱码综合| 亚洲熟妇中文字幕五十中出| 亚洲天堂国产精品一区在线| 看非洲黑人一级黄片| 午夜激情欧美在线| 午夜福利高清视频| 国产片特级美女逼逼视频| 精品福利观看| 亚洲七黄色美女视频| 看黄色毛片网站| 日韩三级伦理在线观看| 级片在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 狠狠狠狠99中文字幕| 真人做人爱边吃奶动态| 床上黄色一级片| 一本精品99久久精品77| 成人国产麻豆网| 欧美丝袜亚洲另类| 中文亚洲av片在线观看爽| 日本爱情动作片www.在线观看 | 一个人观看的视频www高清免费观看| 午夜日韩欧美国产| 欧美一区二区国产精品久久精品| 啦啦啦韩国在线观看视频| 一进一出好大好爽视频| 日本三级黄在线观看| 免费看a级黄色片| 人妻少妇偷人精品九色| 国产女主播在线喷水免费视频网站 | 国产精品99久久久久久久久| 噜噜噜噜噜久久久久久91| 久久久精品大字幕| 亚洲av美国av| 国产高清三级在线| 日韩欧美三级三区| 天堂影院成人在线观看| 99热网站在线观看| 国产av一区在线观看免费| av国产免费在线观看| 日韩 亚洲 欧美在线| 日本-黄色视频高清免费观看| 欧美在线一区亚洲| 少妇高潮的动态图| 搡老妇女老女人老熟妇| 两性午夜刺激爽爽歪歪视频在线观看| 综合色丁香网| 熟妇人妻久久中文字幕3abv| 欧美性猛交╳xxx乱大交人| 成人欧美大片| 国产一区二区激情短视频| 国产综合懂色| 久久久久久久久久黄片| 真实男女啪啪啪动态图| 亚洲一级一片aⅴ在线观看| 别揉我奶头 嗯啊视频| 欧美日本亚洲视频在线播放| 美女高潮的动态| 神马国产精品三级电影在线观看| 简卡轻食公司| 日韩欧美 国产精品| 国产亚洲精品久久久com| 国产精品1区2区在线观看.| 亚洲欧美日韩卡通动漫| 老熟妇乱子伦视频在线观看| 69av精品久久久久久| 国产探花极品一区二区| 日韩成人av中文字幕在线观看 | 日韩欧美精品免费久久| 国产亚洲精品久久久久久毛片| 男女视频在线观看网站免费| 亚洲av五月六月丁香网| 性插视频无遮挡在线免费观看| 黄色配什么色好看| 丝袜美腿在线中文| 黄色日韩在线| 日韩高清综合在线| 欧美三级亚洲精品| 国产欧美日韩一区二区精品| 一区二区三区高清视频在线| 熟女人妻精品中文字幕| 97碰自拍视频| 欧美色欧美亚洲另类二区| 91午夜精品亚洲一区二区三区| 亚洲国产精品成人久久小说 | 麻豆国产97在线/欧美| 精品久久久噜噜| 深夜精品福利| 午夜精品一区二区三区免费看| 免费不卡的大黄色大毛片视频在线观看 | 午夜影院日韩av| 午夜爱爱视频在线播放| 国内少妇人妻偷人精品xxx网站| 久久午夜亚洲精品久久| 国产高清视频在线播放一区| 欧美xxxx性猛交bbbb| 大型黄色视频在线免费观看| 韩国av在线不卡| 亚洲精品在线观看二区| 成年免费大片在线观看| 美女免费视频网站| 中文字幕免费在线视频6| 亚洲精品成人久久久久久| 国产v大片淫在线免费观看| 国产午夜精品久久久久久一区二区三区 | 99热精品在线国产| 国产成人福利小说| 免费av观看视频| 久久久久久国产a免费观看| 国产激情偷乱视频一区二区| 九九热线精品视视频播放| 99热精品在线国产| 日韩制服骚丝袜av| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 美女 人体艺术 gogo| 人人妻人人澡人人爽人人夜夜 | 人妻丰满熟妇av一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 人人妻人人澡欧美一区二区| 极品教师在线视频| 久久久成人免费电影| 国产aⅴ精品一区二区三区波| 黑人高潮一二区| 成人综合一区亚洲| 狂野欧美激情性xxxx在线观看| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 欧美人与善性xxx| 一夜夜www| 成人高潮视频无遮挡免费网站| 色播亚洲综合网| 网址你懂的国产日韩在线| 中国美女看黄片| 日本-黄色视频高清免费观看| 国产国拍精品亚洲av在线观看| 成人一区二区视频在线观看| 亚洲自偷自拍三级| 赤兔流量卡办理| av在线天堂中文字幕| 欧美激情国产日韩精品一区| 亚洲精品影视一区二区三区av| 九九在线视频观看精品| 国产精品一区www在线观看| 日韩av不卡免费在线播放| 桃色一区二区三区在线观看| 男女下面进入的视频免费午夜| 国产乱人偷精品视频| 日韩三级伦理在线观看| 亚洲中文日韩欧美视频| 在线观看美女被高潮喷水网站| 国产成人freesex在线 | 日韩av不卡免费在线播放| 日本撒尿小便嘘嘘汇集6| 一级a爱片免费观看的视频| 亚洲国产精品久久男人天堂| 精品乱码久久久久久99久播| 一级毛片aaaaaa免费看小| 久久久精品欧美日韩精品| 国产成人freesex在线 | 日产精品乱码卡一卡2卡三| 国产精品伦人一区二区| 久久久久久久久大av| 啦啦啦观看免费观看视频高清| 亚洲精品国产成人久久av| 内射极品少妇av片p| 99国产精品一区二区蜜桃av| 国产老妇女一区| 亚洲国产精品成人久久小说 | 免费看美女性在线毛片视频| 亚洲经典国产精华液单| 久久久色成人| 99国产极品粉嫩在线观看| 99久国产av精品| 久久婷婷人人爽人人干人人爱| 国产在线精品亚洲第一网站| or卡值多少钱| 国产白丝娇喘喷水9色精品| 亚洲欧美清纯卡通| 欧美成人一区二区免费高清观看| 国产视频一区二区在线看| 亚洲欧美清纯卡通| 亚洲精品日韩在线中文字幕 | 人妻制服诱惑在线中文字幕| 国产亚洲欧美98| 精品人妻视频免费看| 欧美国产日韩亚洲一区| 国产精品人妻久久久影院| 日日摸夜夜添夜夜添小说| 在线观看美女被高潮喷水网站| 国语自产精品视频在线第100页| 午夜福利视频1000在线观看| 又粗又爽又猛毛片免费看| 欧美中文日本在线观看视频| 免费av毛片视频| 免费人成在线观看视频色| 色综合站精品国产| 国产精品一区二区三区四区免费观看 | 色5月婷婷丁香| 男女视频在线观看网站免费| 日韩欧美 国产精品| 成年免费大片在线观看| 一进一出好大好爽视频| 赤兔流量卡办理| 一区二区三区四区激情视频 | 亚洲国产欧洲综合997久久,| 色5月婷婷丁香| 18禁在线无遮挡免费观看视频 | 欧美丝袜亚洲另类| 久久久成人免费电影| 国产精品电影一区二区三区| 亚洲七黄色美女视频| 国产精品久久久久久久电影| 麻豆乱淫一区二区| 亚洲欧美成人综合另类久久久 | 在线播放国产精品三级| 国产精品日韩av在线免费观看| 国产精品亚洲美女久久久| videossex国产| 在线观看av片永久免费下载| 特大巨黑吊av在线直播| 国产精品一区二区三区四区久久| 亚洲,欧美,日韩| 日韩欧美精品免费久久| 欧美又色又爽又黄视频| 九九爱精品视频在线观看| 免费高清视频大片| 国产精品99久久久久久久久| 美女cb高潮喷水在线观看| 婷婷精品国产亚洲av在线| 日本三级黄在线观看| 国产一区二区激情短视频| 亚洲国产精品合色在线| 精品一区二区三区视频在线| 精品久久久久久久久亚洲| av在线老鸭窝| 综合色丁香网| 国产一级毛片七仙女欲春2| 色视频www国产| 国产精品av视频在线免费观看| 日韩一区二区视频免费看| 一夜夜www| 永久网站在线| 蜜桃亚洲精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品色激情综合| 国产成人一区二区在线| 99久久精品国产国产毛片| 久久人人精品亚洲av| 最近最新中文字幕大全电影3| 听说在线观看完整版免费高清| h日本视频在线播放| 日韩 亚洲 欧美在线| 成人国产麻豆网| 真人做人爱边吃奶动态| 亚洲婷婷狠狠爱综合网| 神马国产精品三级电影在线观看| av黄色大香蕉| 俄罗斯特黄特色一大片| 精品熟女少妇av免费看| 麻豆av噜噜一区二区三区| 精品乱码久久久久久99久播| 欧美+日韩+精品| 美女cb高潮喷水在线观看| 亚洲欧美日韩无卡精品| 人妻丰满熟妇av一区二区三区| 最近中文字幕高清免费大全6| 日本黄大片高清| 如何舔出高潮| 麻豆国产97在线/欧美| 欧美不卡视频在线免费观看| 亚洲人与动物交配视频| 久久这里只有精品中国| 国产亚洲精品综合一区在线观看| 18禁裸乳无遮挡免费网站照片| 人人妻人人澡欧美一区二区| 国产在线精品亚洲第一网站| 69人妻影院| 麻豆久久精品国产亚洲av| 国产v大片淫在线免费观看| 天堂√8在线中文| 99久久九九国产精品国产免费| 日韩 亚洲 欧美在线| www.色视频.com| 亚洲成av人片在线播放无| 97热精品久久久久久| 22中文网久久字幕| 高清毛片免费看| 嫩草影院精品99| 九九在线视频观看精品| 久久精品国产99精品国产亚洲性色| 十八禁网站免费在线| 成人高潮视频无遮挡免费网站| 午夜爱爱视频在线播放| www.色视频.com| 国产成人一区二区在线| 国产欧美日韩一区二区精品| 成人国产麻豆网| 91在线观看av| 国产精品嫩草影院av在线观看| 在线观看一区二区三区| 国产精品人妻久久久久久| 日本一二三区视频观看| www日本黄色视频网| 国产乱人偷精品视频| 欧美性感艳星| 国产免费一级a男人的天堂| 亚洲美女搞黄在线观看 | .国产精品久久| 色视频www国产| 长腿黑丝高跟| 国产精品1区2区在线观看.| 夜夜夜夜夜久久久久| 国产精品一区二区三区四区久久| 欧美不卡视频在线免费观看| 久久人人爽人人片av| 亚洲精品色激情综合| 国产精品国产高清国产av| 国产中年淑女户外野战色| 亚州av有码| 99久久九九国产精品国产免费| 秋霞在线观看毛片| 亚洲av美国av| 白带黄色成豆腐渣| 一级黄片播放器| 久久婷婷人人爽人人干人人爱| 91狼人影院| 久99久视频精品免费| 久久久久久伊人网av| 日本熟妇午夜| 中文字幕精品亚洲无线码一区| 亚洲自偷自拍三级| 一个人看的www免费观看视频| 99热这里只有是精品在线观看| 长腿黑丝高跟| 亚洲成人av在线免费| 如何舔出高潮| 国产一区二区三区在线臀色熟女| 日本黄大片高清|