• <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富集與降解中的研究
    在线精品无人区一区二区三| 久久久国产欧美日韩av| 黄色毛片三级朝国网站| 亚洲内射少妇av| 99热网站在线观看| 男的添女的下面高潮视频| 美女中出高潮动态图| 欧美日韩av久久| 国产又色又爽无遮挡免| www.精华液| 日韩一本色道免费dvd| 欧美黄色片欧美黄色片| 黑人巨大精品欧美一区二区蜜桃| 亚洲av免费高清在线观看| 欧美日韩视频高清一区二区三区二| 亚洲欧美一区二区三区黑人 | 亚洲精华国产精华液的使用体验| 国产成人aa在线观看| 免费av中文字幕在线| 日产精品乱码卡一卡2卡三| 一级毛片黄色毛片免费观看视频| 国产伦理片在线播放av一区| 免费黄频网站在线观看国产| 女的被弄到高潮叫床怎么办| 麻豆精品久久久久久蜜桃| 啦啦啦在线免费观看视频4| 日韩中文字幕视频在线看片| 亚洲av男天堂| 国产极品天堂在线| 欧美黄色片欧美黄色片| 97在线人人人人妻| 日韩 亚洲 欧美在线| 一区二区三区乱码不卡18| 久久韩国三级中文字幕| 欧美国产精品一级二级三级| 日韩中字成人| 超色免费av| 亚洲精品自拍成人| 精品国产露脸久久av麻豆| 欧美av亚洲av综合av国产av | 中文精品一卡2卡3卡4更新| 亚洲美女黄色视频免费看| 久久精品国产亚洲av高清一级| 亚洲第一av免费看| 久久久久国产精品人妻一区二区| 高清不卡的av网站| 欧美日韩视频高清一区二区三区二| 亚洲精品自拍成人| 免费观看在线日韩| 免费大片黄手机在线观看| 久久久久久人妻| 在现免费观看毛片| av在线播放精品| 色吧在线观看| 最黄视频免费看| 成人黄色视频免费在线看| 在线天堂最新版资源| 九九爱精品视频在线观看| 水蜜桃什么品种好| 亚洲四区av| 性色avwww在线观看| 精品一区在线观看国产| av免费观看日本| 美女大奶头黄色视频| 老司机影院毛片| 成年人免费黄色播放视频| 亚洲视频免费观看视频| 国产免费又黄又爽又色| 人妻人人澡人人爽人人| 亚洲成人手机| 欧美精品国产亚洲| 美女主播在线视频| 久久免费观看电影| 精品国产一区二区三区久久久樱花| 免费看不卡的av| 一区二区av电影网| 国产精品蜜桃在线观看| 国产日韩欧美视频二区| 国产精品免费视频内射| 成年人午夜在线观看视频| 最近2019中文字幕mv第一页| a级毛片在线看网站| 国产高清不卡午夜福利| 大陆偷拍与自拍| 满18在线观看网站| 纵有疾风起免费观看全集完整版| 国产精品.久久久| 久久久久久久久久久免费av| 99国产综合亚洲精品| 大香蕉久久网| 国产伦理片在线播放av一区| 久久久欧美国产精品| 亚洲男人天堂网一区| 成年动漫av网址| 久久亚洲国产成人精品v| 欧美人与善性xxx| 国产在线免费精品| 考比视频在线观看| 精品人妻在线不人妻| 国产午夜精品一二区理论片| 美女午夜性视频免费| 国产一区有黄有色的免费视频| 免费黄频网站在线观看国产| 午夜免费观看性视频| 少妇的逼水好多| 亚洲精品一二三| 亚洲第一青青草原| 久久久精品免费免费高清| 亚洲国产精品国产精品| 精品一区二区三卡| 亚洲精品一二三| 成人午夜精彩视频在线观看| 日韩在线高清观看一区二区三区| 日韩,欧美,国产一区二区三区| 国产又色又爽无遮挡免| 高清视频免费观看一区二区| 亚洲久久久国产精品| h视频一区二区三区| 午夜影院在线不卡| 蜜桃国产av成人99| 妹子高潮喷水视频| 国产精品蜜桃在线观看| 久久久久久久亚洲中文字幕| 久久久国产欧美日韩av| 日韩视频在线欧美| 国产伦理片在线播放av一区| 免费高清在线观看日韩| 女性生殖器流出的白浆| 亚洲成av片中文字幕在线观看 | 久久久久久人妻| 国产精品久久久久久久久免| 国产又爽黄色视频| 久久人人97超碰香蕉20202| 十八禁高潮呻吟视频| 涩涩av久久男人的天堂| www.自偷自拍.com| av线在线观看网站| 黑人巨大精品欧美一区二区蜜桃| 伊人久久大香线蕉亚洲五| 久久久久久人妻| 国产欧美日韩一区二区三区在线| www.自偷自拍.com| 日韩一区二区视频免费看| 久久精品国产a三级三级三级| 久热久热在线精品观看| 性色av一级| 丰满饥渴人妻一区二区三| 韩国高清视频一区二区三区| 激情视频va一区二区三区| 国产精品国产三级国产专区5o| 日韩一区二区三区影片| 曰老女人黄片| 一区二区三区激情视频| 中文乱码字字幕精品一区二区三区| 精品少妇内射三级| 91精品国产国语对白视频| 卡戴珊不雅视频在线播放| 91在线精品国自产拍蜜月| 日韩欧美精品免费久久| 午夜91福利影院| 亚洲精品日本国产第一区| 欧美 亚洲 国产 日韩一| 精品少妇一区二区三区视频日本电影 | 我的亚洲天堂| 精品第一国产精品| 午夜老司机福利剧场| 午夜激情av网站| 九色亚洲精品在线播放| 水蜜桃什么品种好| 纯流量卡能插随身wifi吗| 国产乱人偷精品视频| 天天影视国产精品| 亚洲三级黄色毛片| 国产精品久久久av美女十八| 亚洲国产毛片av蜜桃av| 午夜免费观看性视频| 久久精品国产亚洲av涩爱| 亚洲成色77777| 十八禁网站网址无遮挡| 久久午夜福利片| av卡一久久| 国产片特级美女逼逼视频| 久久久久久免费高清国产稀缺| 超碰97精品在线观看| 午夜福利视频在线观看免费| 岛国毛片在线播放| 国产精品国产三级专区第一集| 十分钟在线观看高清视频www| 中文字幕最新亚洲高清| 美国免费a级毛片| 亚洲精品成人av观看孕妇| 日日啪夜夜爽| 18+在线观看网站| 欧美人与善性xxx| 丝袜美腿诱惑在线| 日本av免费视频播放| 国产成人免费观看mmmm| 久久久久国产一级毛片高清牌| 国产一区二区 视频在线| 性色avwww在线观看| 久久久精品区二区三区| 久久精品国产鲁丝片午夜精品| av在线app专区| 国产精品久久久av美女十八| 亚洲国产精品999| 日本-黄色视频高清免费观看| 美女主播在线视频| 午夜久久久在线观看| 在线观看美女被高潮喷水网站| 国产免费视频播放在线视频| 日本欧美国产在线视频| 成人二区视频| 成人午夜精彩视频在线观看| 国产精品.久久久| 日本午夜av视频| 人人澡人人妻人| 国产免费一区二区三区四区乱码| 五月天丁香电影| 三级国产精品片| 欧美 日韩 精品 国产| 国产亚洲av片在线观看秒播厂| 日韩欧美一区视频在线观看| 亚洲人成77777在线视频| 亚洲欧美一区二区三区黑人 | 亚洲精品国产一区二区精华液| 国产人伦9x9x在线观看 | 成人午夜精彩视频在线观看| 久久久欧美国产精品| 中文字幕另类日韩欧美亚洲嫩草| 91国产中文字幕| 亚洲精品乱久久久久久| 少妇猛男粗大的猛烈进出视频| 一区二区三区激情视频| 黑人巨大精品欧美一区二区蜜桃| 午夜影院在线不卡| 国产伦理片在线播放av一区| videosex国产| 天天操日日干夜夜撸| 制服诱惑二区| 国产成人免费观看mmmm| 一级毛片电影观看| 免费观看性生交大片5| 人妻少妇偷人精品九色| 中文字幕制服av| 欧美+日韩+精品| 久久免费观看电影| 亚洲精品乱久久久久久| av卡一久久| 99久久精品国产国产毛片| 97精品久久久久久久久久精品| 亚洲伊人色综图| 久久这里有精品视频免费| 久久人妻熟女aⅴ| 亚洲综合色惰| 丝袜美足系列| 精品少妇一区二区三区视频日本电影 | 99精国产麻豆久久婷婷| 26uuu在线亚洲综合色| 国产人伦9x9x在线观看 | 精品国产一区二区三区四区第35| 另类亚洲欧美激情| 国产欧美亚洲国产| 亚洲精品,欧美精品| 夫妻午夜视频| 97人妻天天添夜夜摸| 国产日韩欧美亚洲二区| 国产av精品麻豆| 亚洲人成电影观看| 久久久久精品久久久久真实原创| 中文字幕色久视频| 亚洲精品av麻豆狂野| 侵犯人妻中文字幕一二三四区| 亚洲欧美色中文字幕在线| 亚洲国产精品999| 纵有疾风起免费观看全集完整版| 国产一区亚洲一区在线观看| 9191精品国产免费久久| 亚洲五月色婷婷综合| 久久久久久免费高清国产稀缺| 一区二区日韩欧美中文字幕| 国产成人精品福利久久| 国产成人精品久久久久久| 日韩精品有码人妻一区| 国产一区二区 视频在线| 精品少妇一区二区三区视频日本电影 | 一区福利在线观看| 亚洲经典国产精华液单| 一区二区三区乱码不卡18| 日日爽夜夜爽网站| 麻豆精品久久久久久蜜桃| 欧美在线黄色| 大香蕉久久网| 成人手机av| 精品一区在线观看国产| 我要看黄色一级片免费的| 国产在线一区二区三区精| 国产精品av久久久久免费| 国产精品一国产av| 国产精品女同一区二区软件| 热99国产精品久久久久久7| 在线亚洲精品国产二区图片欧美| 久久狼人影院| 国产av精品麻豆| 免费在线观看黄色视频的| 美女xxoo啪啪120秒动态图| av女优亚洲男人天堂| 国产精品久久久av美女十八| 亚洲在久久综合| 男女免费视频国产| 国产成人精品婷婷| 啦啦啦中文免费视频观看日本| 99久久精品国产国产毛片| 狠狠精品人妻久久久久久综合| 男女无遮挡免费网站观看| 国产男女超爽视频在线观看| 如日韩欧美国产精品一区二区三区| 一边亲一边摸免费视频| 亚洲婷婷狠狠爱综合网| 久久久久国产精品人妻一区二区| 男女午夜视频在线观看| 不卡视频在线观看欧美| 日韩av在线免费看完整版不卡| 久久久久视频综合| 国产1区2区3区精品| 建设人人有责人人尽责人人享有的| 久久狼人影院| 中文字幕人妻丝袜制服| 日韩伦理黄色片| 人人妻人人添人人爽欧美一区卜| 水蜜桃什么品种好| 免费少妇av软件| 老司机影院毛片| 又大又黄又爽视频免费| 国产综合精华液| 一级黄片播放器| 久久这里有精品视频免费| 丝袜美腿诱惑在线| 国产成人91sexporn| 午夜免费鲁丝| 久久精品夜色国产| 精品一品国产午夜福利视频| 亚洲在久久综合| 只有这里有精品99| 中文字幕色久视频| 国产精品熟女久久久久浪| 久久久国产欧美日韩av| 亚洲第一区二区三区不卡| 国产一区二区在线观看av| 免费av中文字幕在线| 国产片特级美女逼逼视频| 午夜影院在线不卡| 久久免费观看电影| 两个人看的免费小视频| 黄网站色视频无遮挡免费观看| 日本91视频免费播放| 性高湖久久久久久久久免费观看| 婷婷色麻豆天堂久久| 色婷婷av一区二区三区视频| 97人妻天天添夜夜摸| 亚洲男人天堂网一区| 如日韩欧美国产精品一区二区三区| 色视频在线一区二区三区| 少妇人妻精品综合一区二区| 精品一区二区三卡| 国产福利在线免费观看视频| 成年女人毛片免费观看观看9 | av免费观看日本| av国产精品久久久久影院| 丰满饥渴人妻一区二区三| 亚洲欧美成人精品一区二区| 麻豆av在线久日| 欧美另类一区| 成年人午夜在线观看视频| 美女中出高潮动态图| 黄频高清免费视频| 免费av中文字幕在线| 建设人人有责人人尽责人人享有的| 久久久久国产网址| 日韩中文字幕视频在线看片| 水蜜桃什么品种好| 国产成人午夜福利电影在线观看| 亚洲美女黄色视频免费看| 精品少妇内射三级| 欧美成人精品欧美一级黄| 国产男女超爽视频在线观看| av片东京热男人的天堂| 99久国产av精品国产电影| av视频免费观看在线观看| av在线观看视频网站免费| 久久综合国产亚洲精品| 亚洲国产成人一精品久久久| 最近中文字幕2019免费版| 777米奇影视久久| 欧美成人午夜精品| 狠狠婷婷综合久久久久久88av| 日本vs欧美在线观看视频| 高清欧美精品videossex| 免费少妇av软件| 1024视频免费在线观看| 午夜免费鲁丝| 亚洲av在线观看美女高潮| av.在线天堂| 国产精品 欧美亚洲| 麻豆乱淫一区二区| 欧美人与善性xxx| 日本vs欧美在线观看视频| 91精品国产国语对白视频| 激情视频va一区二区三区| 午夜福利在线免费观看网站| 少妇人妻久久综合中文| 人妻少妇偷人精品九色| 亚洲国产精品999| 国产欧美日韩综合在线一区二区| www.精华液| 女人被躁到高潮嗷嗷叫费观| 国产探花极品一区二区| 婷婷色综合大香蕉| 少妇猛男粗大的猛烈进出视频| 伊人久久国产一区二区| 亚洲欧美精品自产自拍| 天美传媒精品一区二区| 涩涩av久久男人的天堂| 午夜日本视频在线| 丝袜美足系列| 国产成人精品一,二区| 日本av免费视频播放| 午夜91福利影院| 国产不卡av网站在线观看| 久久精品久久久久久噜噜老黄| 热re99久久国产66热| 国产一区有黄有色的免费视频| 满18在线观看网站| 国产亚洲最大av| 欧美成人午夜免费资源| 精品福利永久在线观看| 国产免费视频播放在线视频| 午夜免费观看性视频| 如日韩欧美国产精品一区二区三区| 91在线精品国自产拍蜜月| 女人高潮潮喷娇喘18禁视频| 视频在线观看一区二区三区| 国产免费一区二区三区四区乱码| 热99久久久久精品小说推荐| 搡女人真爽免费视频火全软件| 色播在线永久视频| 国产一级毛片在线| 岛国毛片在线播放| 七月丁香在线播放| 春色校园在线视频观看| 伦精品一区二区三区| 中国三级夫妇交换| 少妇人妻精品综合一区二区| 深夜精品福利| 精品酒店卫生间| 亚洲国产av影院在线观看| 九色亚洲精品在线播放| 成人亚洲精品一区在线观看| 激情五月婷婷亚洲| 午夜福利视频在线观看免费| 久久精品亚洲av国产电影网| 在线观看免费高清a一片| 久久精品aⅴ一区二区三区四区 | 69精品国产乱码久久久| 老司机影院成人| 国产免费现黄频在线看| av一本久久久久| 巨乳人妻的诱惑在线观看| 最黄视频免费看| 少妇人妻精品综合一区二区| 最近最新中文字幕大全免费视频 | 黄色一级大片看看| av免费在线看不卡| 韩国高清视频一区二区三区| 丝袜美腿诱惑在线| 中国国产av一级| 国产成人精品一,二区| 欧美bdsm另类| 日韩一区二区三区影片| 免费女性裸体啪啪无遮挡网站| 日韩成人av中文字幕在线观看| 一个人免费看片子| 青草久久国产| 免费日韩欧美在线观看| 国产精品二区激情视频| 日日撸夜夜添| 精品国产超薄肉色丝袜足j| 在线天堂最新版资源| 美女xxoo啪啪120秒动态图| 老女人水多毛片| 欧美日韩成人在线一区二区| 国产精品久久久久成人av| 亚洲精品美女久久av网站| 国产精品欧美亚洲77777| 亚洲精品一区蜜桃| 亚洲精品久久久久久婷婷小说| 在线观看免费高清a一片| 亚洲,欧美,日韩| 熟女少妇亚洲综合色aaa.| 伊人久久国产一区二区| 99香蕉大伊视频| 一本久久精品| 久久婷婷青草| 少妇被粗大猛烈的视频| 亚洲精品成人av观看孕妇| 男的添女的下面高潮视频| av国产久精品久网站免费入址| 蜜桃在线观看..| 国产精品 欧美亚洲| 成年女人在线观看亚洲视频| 欧美精品国产亚洲| 国产日韩欧美在线精品| 国产精品秋霞免费鲁丝片| 纯流量卡能插随身wifi吗| 久久久久精品久久久久真实原创| 啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| 日日爽夜夜爽网站| 国产精品二区激情视频| 日韩大片免费观看网站| 国产乱人偷精品视频| 日韩中字成人| 人妻少妇偷人精品九色| 国产成人精品久久久久久| 9191精品国产免费久久| 精品酒店卫生间| 我的亚洲天堂| 尾随美女入室| 精品午夜福利在线看| 男人舔女人的私密视频| av在线观看视频网站免费| 女性被躁到高潮视频| 香蕉丝袜av| 日韩欧美精品免费久久| 亚洲中文av在线| 91在线精品国自产拍蜜月| 十八禁高潮呻吟视频| 久久99精品国语久久久| 国产片特级美女逼逼视频| 国产熟女欧美一区二区| 中文字幕最新亚洲高清| 日韩视频在线欧美| 国产精品秋霞免费鲁丝片| 久久久久久久大尺度免费视频| 91aial.com中文字幕在线观看| 亚洲精品成人av观看孕妇| 久久久亚洲精品成人影院| 精品国产一区二区三区久久久樱花| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品视频女| 天堂俺去俺来也www色官网| 这个男人来自地球电影免费观看 | 亚洲欧美精品综合一区二区三区 | 一级毛片 在线播放| 女的被弄到高潮叫床怎么办| 久久99蜜桃精品久久| 亚洲五月色婷婷综合| 午夜福利在线观看免费完整高清在| 亚洲四区av| 最近最新中文字幕大全免费视频 | 肉色欧美久久久久久久蜜桃| 国语对白做爰xxxⅹ性视频网站| 叶爱在线成人免费视频播放| 亚洲 欧美一区二区三区| 亚洲第一青青草原| 成人18禁高潮啪啪吃奶动态图| 伦精品一区二区三区| 中文字幕最新亚洲高清| 蜜桃在线观看..| 国产日韩一区二区三区精品不卡| 亚洲精华国产精华液的使用体验| 精品一区在线观看国产| 免费看不卡的av| 丰满迷人的少妇在线观看| 久久精品aⅴ一区二区三区四区 | 精品少妇内射三级| 国产在线免费精品| 国产人伦9x9x在线观看 | 亚洲av男天堂| 91成人精品电影| 桃花免费在线播放| 老熟女久久久| 国产精品秋霞免费鲁丝片| 日韩制服丝袜自拍偷拍| 精品久久久久久电影网| 国产深夜福利视频在线观看| 亚洲欧美成人综合另类久久久| 在线亚洲精品国产二区图片欧美| 免费播放大片免费观看视频在线观看| 丝瓜视频免费看黄片| 精品一区二区三区四区五区乱码 | videosex国产| 午夜91福利影院| 肉色欧美久久久久久久蜜桃| 精品少妇久久久久久888优播| 久久久久精品人妻al黑| av国产精品久久久久影院| 九草在线视频观看| 国产成人a∨麻豆精品| 精品国产一区二区三区久久久樱花| 99久久精品国产国产毛片| 午夜福利网站1000一区二区三区| 亚洲国产色片| 日韩av在线免费看完整版不卡| 少妇 在线观看| 久久精品亚洲av国产电影网| 欧美亚洲日本最大视频资源| 少妇被粗大猛烈的视频| 乱人伦中国视频| 日韩三级伦理在线观看| 精品视频人人做人人爽| 91久久精品国产一区二区三区| kizo精华| 啦啦啦中文免费视频观看日本| 日韩熟女老妇一区二区性免费视频| av国产久精品久网站免费入址|