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

    長江口及鄰近海域浮游植物現(xiàn)存量的上下行控制分析

    2014-10-27 05:37:34宋洪軍張學(xué)雷王保棟孫霞王曉辛明
    海洋學(xué)報(bào) 2014年8期
    關(guān)鍵詞:模型

    宋洪軍,張學(xué)雷*,王保棟,孫霞,王曉,辛明

    (1.國家海洋局 第一海洋研究所,山東 青島266061;2.海洋生態(tài)環(huán)境科學(xué)與工程國家海洋局重點(diǎn)實(shí)驗(yàn)室,山東 青島266061)

    1 引言

    根據(jù)營養(yǎng)級聯(lián)關(guān)系(Cascading Trophic Interactions)[1],浮游植物現(xiàn)存量主要受營養(yǎng)鹽水平與光照條件的上行控制(Bottom-up Control)和藻類攝食者的下行控制(Top-down Control)兩效應(yīng)影響[2-3]。上行效應(yīng)一般是指浮游植物生長受環(huán)境條件影響,下行效應(yīng)主要是指浮游動(dòng)物通過攝食作用影響浮游植物生物量的累積。上下行效應(yīng)對浮游植物現(xiàn)存量控制的相對作用在不同水域和季節(jié)是不同的,例如,Sommer[4]現(xiàn)場培養(yǎng)實(shí)驗(yàn)研究指出在寡營養(yǎng)的紅海北部浮游植物生物量主要受上行效應(yīng)控制,下行控制作用相對較弱;Goericke[5]研究指出阿拉伯海雨季的高營養(yǎng)鹽低葉綠素現(xiàn)象(“high-nutrient,low-chlorophyll”,HNLC)主要受浮游動(dòng)物的下行控制作用;Thompson等[6]在澳大利亞Huon河口浮游植物水華過程研究發(fā)現(xiàn),冬季該區(qū)域硅藻的優(yōu)勢主要受光照條件的上行效應(yīng)影響,夏季甲藻水華則與潛在的氮限制有關(guān)。

    近年來,國內(nèi)學(xué)者關(guān)于長江口水域浮游植物的營養(yǎng)鹽限制研究已開展了較多的工作,并形成了一定的科學(xué)認(rèn)識。研究結(jié)果顯示,長江口及鄰近水域溶解無機(jī)氮(DIN)和硅酸鹽(Si O3-Si)的濃度均較高(年平均值分別約為20μmol/L和15μmol/L),而磷酸鹽(PO4-P)濃度較低(年平均值約0.5μmol/L)[7];硅氮比值(Si/DIN)大于1或小于1的情況均有,但四季的Si O3-Si濃度遠(yuǎn)高于2μmol/L的閾值,因此長江口及鄰近海域不構(gòu)成硅限制[8];氮磷比值(DIN/P)變化范圍較大,一般高于Redfield比值(16∶1),但PO4-P濃度基本也都高于0.2μmol/L,在浮游植物生長的限制性閾值(0.1μmol/L)之上[9];但長期變化來看,DIN濃度的持續(xù)升高以及PO4-P濃度的持續(xù)低值致使DIN/P的升高趨勢顯著[10],導(dǎo)致長江口及鄰近海域浮游植物生長存在潛在的磷限制。光照條件對長江口浮游植物生長的影響也已引起大家重視,例如,方濤[11]通過現(xiàn)場培養(yǎng)實(shí)驗(yàn)研究了不同光照條件下長江口水域浮游植物生長和營養(yǎng)鹽吸收作用機(jī)制以及不同浮游生物類群對碳貢獻(xiàn)率的變化;孫百曄[12]通過對長江口及鄰近海域浮游植物生物量平面分布、垂直分布和年代變化光照效應(yīng)的綜合分析,證實(shí)光照是影響該海域赤潮發(fā)生的重要因素。而浮游動(dòng)物對浮游植物的下行控制研究目前主要集中在研究區(qū)域水華過程中的攝食實(shí)驗(yàn)方面[13-14],對浮游植物現(xiàn)存量方面的影響分析尚未見報(bào)道。因此,針對控制長江口及鄰近海域浮游植物現(xiàn)存量的上下行效應(yīng)仍需開展綜合研究。

    廣義相加模型(Generalized Additive Model,GAM)是一種基于數(shù)據(jù)分析的數(shù)學(xué)模型,數(shù)據(jù)決定著響應(yīng)變量與預(yù)測因子間的關(guān)系,而不是根據(jù)影響機(jī)制所假設(shè)的某種參數(shù)關(guān)系[15]。該模型能有效揭示響應(yīng)變量與預(yù)測因子間非線性和非單調(diào)潛在聯(lián)系,在生態(tài)學(xué)和海洋學(xué)研究中已經(jīng)廣泛應(yīng)用。國際上已有利用GAM 模型在湖泊[16]、海灣[17]與海峽[18]等區(qū)域研究葉綠素與不同環(huán)境因子間關(guān)系的報(bào)道,國內(nèi)研究則在異養(yǎng)細(xì)菌[19]、浮游動(dòng)物[20]、漁業(yè)資源[21-22]等與環(huán)境因子間關(guān)系的研究中引入了GAM模型分析方法。

    本文利用2010-2011年度的調(diào)查資料研究了長江口及鄰近海域營養(yǎng)鹽、光照(懸浮物)和浮游動(dòng)物對浮游植物現(xiàn)存量(以葉綠素a濃度表征)空間分布和季節(jié)變化的影響,以GAM模型為工具量化分析了影響葉綠素時(shí)空變化的主要環(huán)境因子,探討了上下行效應(yīng)對該區(qū)域浮游植物現(xiàn)存量的控制,以期豐富對長江口及鄰近海域浮游植物生態(tài)學(xué)的認(rèn)識,為生態(tài)系統(tǒng)動(dòng)力學(xué)研究提供數(shù)據(jù)支撐。

    2 材料與方法

    2.1 數(shù)據(jù)來源

    長江口及鄰近海域現(xiàn)場調(diào)查站位分布如圖1所示。共執(zhí)行3個(gè)航次,調(diào)查時(shí)間分別為2010年8月(夏季)、2010年11月(秋季)和2011年5月(春季)。樣品采集、處理和測定參照《海洋調(diào)查規(guī)范》[23-24]和《海洋監(jiān)測規(guī)范》[25]執(zhí)行,其中葉綠素a(表征浮游植物現(xiàn)存量)根據(jù)萃取熒光法以Turner-Designs 700型熒光計(jì)測定,溫鹽由RBR(XR-620)儀器測定,營養(yǎng)鹽使用SKALAR自動(dòng)分析儀測定(其中DIN為NO3-N、NO2-N和NH4-N總和),懸浮物采用重量法測定,浮游動(dòng)物生物量指淺水Ⅰ型浮游生物網(wǎng)采樣品的濕重生物量(含膠質(zhì)浮游動(dòng)物)。

    圖1 2010—2011年長江口及鄰近海域調(diào)查采樣站位分布(n=24)Fig.1 Sampling stations off the Changjiang Estuary in 2010—2011(n=24)

    2.2 GAM 模型

    利用R程序包mgcv中的gam函數(shù)[26]模擬研究海域葉綠素對各環(huán)境因子的響應(yīng)。環(huán)境因子包括溫度(表層,SST)、鹽度(表層,SSS)、懸浮物(TSS)、DIN、PO4-P、Si O3-Si和浮游動(dòng)物生物量(ZB)。由于葉綠素值屬于高度右偏態(tài)分布,分析時(shí)取log10(Chl)轉(zhuǎn)換以達(dá)到近似正態(tài)分布[27],同理對于浮游動(dòng)物生物量和TSS也采取對數(shù)轉(zhuǎn)換。分析時(shí)采取步退式模型選擇方法(根據(jù)統(tǒng)計(jì)學(xué)p值和R2值以及模型廣義交叉驗(yàn)證分值即GCV值,逐步減少變量),得到最適模型及其響應(yīng)因子(即影響葉綠素變化的主要環(huán)境因子)。

    3 結(jié)果

    3.1 現(xiàn)場觀測結(jié)果

    2010—2011年度三航次葉綠素觀測值在0.28~32.01μg/L之間,平均值2.75μg/L。季節(jié)比較以夏季葉綠素值最高,春季次之,秋季最低,平均值依次為:5.54μg/L、1.92μg/L、0.78μg/L。春季,葉綠素高值區(qū)主要分布在122.5°E以東海域(>3.0μg/L,圖2a);夏季,在長江口門外南北兩側(cè)形成兩個(gè)高值區(qū)(>5.0μg/L,圖2b);秋季,葉綠素濃度普遍較低,僅在3個(gè)站位稍高于1.0μg/L(圖2c)。葉綠素高值區(qū)的季節(jié)分布與孫軍等[28]報(bào)道相似。

    圖2 2010—2011年度長江口及鄰近海域葉綠素平面分布Fig.2 Spatial distribution of chlorophyll concentration off the Changjiang Estuary in 2010—2011

    同航次環(huán)境因子的觀測值(表1)顯示:春季,與整個(gè)調(diào)查區(qū)相比,葉綠素高值區(qū)的溫度平均值相同,鹽度稍高,TSS顯著低值,而營養(yǎng)鹽(DIN、PO4-P和Si O3-Si)也相對較低,且PO4-P和Si O3-Si濃度在3個(gè)季節(jié)中最低;夏季,葉綠素高值區(qū)的溫度相對整個(gè)調(diào)查區(qū)幾乎無變化,鹽度稍低(不同于春季),TSS顯著低值,而3項(xiàng)營養(yǎng)鹽也低于全調(diào)查區(qū);秋季葉綠素高值區(qū)溫度、鹽度和TSS的平均值較之全調(diào)查區(qū)稍高,而3項(xiàng)營養(yǎng)鹽在葉綠素高值區(qū)的平均值相對較低(這與春夏兩季類似)。DIN/P在全調(diào)查區(qū)春夏秋三季節(jié)依次為37.4、34.7、19.7,而在葉綠素高值區(qū)依次為43.3、32.9、19.2;Si/DIN在全調(diào)查區(qū)春夏秋三季節(jié)依次為0.44、0.55、0.76,而在葉綠素高值區(qū)依次為0.43、0.61、0.76。這說明葉綠素高值區(qū)的DIN/P在春季相對全調(diào)查區(qū)較高,而在夏季和秋季稍低,Si/DIN僅在夏季相對較高,春季和秋季則與全調(diào)查區(qū)基本相同。

    表1 2010—2011年長江口及鄰近海域環(huán)境因子觀測值Tab.1 Variations of environmental factors off the Changjiang Estuary in 2010—2011

    浮游動(dòng)物濕重生物量三航次觀測值在1.25~2 563.79 mg/m3之間,平均值為362.62 mg/m3。季節(jié)比較來看,春季浮游動(dòng)物生物量最高,夏季次之,而秋季最低,平均值依次為:571.54 mg/m3、371.31 mg/m3、145.02 mg/m3。春季(如圖3a所示),浮游動(dòng)物高值區(qū)(>1 000 mg/m3)分布在調(diào)查區(qū)中央的B4和C4站位及西北部的A2站位,其中B4站位亦為葉綠素高值區(qū),C4站位靠近葉綠素高值區(qū);夏季(如圖3b所示),浮游動(dòng)物高值區(qū)(>1 000 mg/m3)分布在調(diào)查區(qū)中央的B3和C3站位,其中C3站位處于葉綠素高值區(qū)內(nèi),而B3站位處于兩個(gè)葉綠素高值區(qū)之間;秋季(圖3c),浮游動(dòng)物相對高值區(qū)(>300 mg/m3)處于近岸的B1和C1站位以及中央的C4站位,亦與葉綠素高值區(qū)的C2和C5站位相鄰,且浮游動(dòng)物生物量在3個(gè)季節(jié)中的相對低值與葉綠素趨勢一致。

    3.2 GAM模型結(jié)果

    GAM模型結(jié)果(見表2)顯示:溫度(SST)、鹽度(SSS)、懸浮物(TSS)、營養(yǎng)鹽(DIN、PO4和Si O3)和浮游動(dòng)物生物量(ZB)變化對葉綠素時(shí)空變化的貢獻(xiàn)合計(jì)可達(dá)74.6%(模型1),依次去除統(tǒng)計(jì)學(xué)意義上不顯著(p>0.05)的TSS(模型2)、浮游動(dòng)物生物量(模型3)和溫度(模型4)后擬合結(jié)果僅有微小變化,其中GCV值僅下降0.001,而R2最大降低不足0.03;去除統(tǒng)計(jì)學(xué)意義上顯著(p<0.05)的鹽度后R2顯著降低為0.535,GCV值也升至0.099(模型5),繼續(xù)去除DIN后,R2顯著降低為0.428,GCV值也升至0.118(模型6)。這表明上述因素中影響長江口及鄰近海域葉綠素變化的主要因子為鹽度和營養(yǎng)鹽(DIN、PO4-P和Si O3-Si三者共同作用),而葉綠素變化與TSS、浮游動(dòng)物生物量和溫度三因子的直接相關(guān)性不顯著。

    圖3 2010—2011年度長江口及鄰近海域浮游動(dòng)物生物量平面分布Fig.3 Spatial distribution of zooplankton biomass off the Changjiang Estuary in 2010—2011

    表2 GAM模型統(tǒng)計(jì)結(jié)果(樣本數(shù)n=72)Tab.2 Statistical results of generalized additive models(samples number:n=72)

    圖4示出了主要影響因子(鹽度、DIN、PO4-P和Si O3-Si)與葉綠素變化間的關(guān)系:葉綠素(對數(shù)值)變化隨著鹽度和PO4-P升高呈逐漸下降的趨勢;DIN濃度較低時(shí)對應(yīng)的葉綠素變化無明顯規(guī)律,在DIN濃度約32μmol/L時(shí)葉綠素達(dá)到一個(gè)峰值,然后降低且在DIN約42μmol/L時(shí)至葉綠素最低,隨著DIN的繼續(xù)升高,葉綠素亦有升高趨勢,但由于數(shù)據(jù)點(diǎn)較少以及置信區(qū)間較為離散導(dǎo)致結(jié)果可信度較低;Si O3-Si與葉綠素變化間的明顯規(guī)律為在Si O3-Si濃度約22μmol/L時(shí)葉綠素表現(xiàn)為一高值,Si O3-Si較低時(shí)葉綠素值較低,而Si O3-Si較高時(shí)葉綠素值變化較大。此外,受長江沖淡水的影響,研究區(qū)域的鹽度與DIN、PO4-P、Si O3-Si等營養(yǎng)鹽及TSS之間存在顯著的相關(guān)關(guān)系(p<0.001,見圖5)。這屬于GAM模型中的變量間共相關(guān)問題,從直接聯(lián)系上分析應(yīng)該去除鹽度因子,但去除鹽度后模型結(jié)果變化較大(如表2模型5所示)。這說明鹽度對葉綠素變化的影響可能體現(xiàn)了營養(yǎng)鹽和光照條件等因子的作用。

    圖4 長江口及鄰近海域葉綠素變化主要影響因子分析Fig.4 Partial effects of the main environmental factors on the chlorophyll variability off the Changjiang Estuary

    4 討論

    水體葉綠素濃度主要取決于浮游植物生長率、垂向擴(kuò)散與混合、水平移動(dòng)、沉降損失與被攝食消耗等因素[29]。浮游植物的生長率主要受營養(yǎng)鹽和光照的上行效應(yīng)影響,而被攝食消耗則以浮游動(dòng)物的下行效應(yīng)為主。本文觀測結(jié)果發(fā)現(xiàn),長江口及鄰近海域3個(gè)季節(jié)葉綠素高值區(qū)營養(yǎng)鹽相對全區(qū)平均值較低,夏春兩季葉綠素高值區(qū)(也是全年葉綠素高值)的水體TSS相對較低,夏季葉綠素高值區(qū)鹽度較低,而溫度在葉綠素高值區(qū)差別不大;浮游動(dòng)物生物量高值區(qū)與葉綠素高值區(qū)未體現(xiàn)出較好的對應(yīng)關(guān)系,但兩者基本都有交叉或位置鄰近。GAM模型結(jié)果進(jìn)一步顯示,營養(yǎng)鹽是影響長江口及鄰近海域葉綠素變化的主要環(huán)境因子,鹽度在解釋葉綠素變化趨勢上亦起了較大作用(由于鹽度與營養(yǎng)鹽和TSS顯著相關(guān),導(dǎo)致鹽度體現(xiàn)了營養(yǎng)鹽和光照條件等因子的共同作用),而浮游動(dòng)物生物量和溫度等因子直接影響不顯著。

    圖5 長江口及鄰近海域鹽度與營養(yǎng)鹽和懸浮物之間的相關(guān)性(n=72)Fig.5 The correlations bet ween surface salinity and nutrients &total suspended solids off the Changjiang Estuary(n=72)

    長江沖淡水低鹽、高渾濁度和高營養(yǎng)鹽的特性,導(dǎo)致了長江口及鄰近海域形成了鹽度、TSS、營養(yǎng)鹽等一系列梯度變化,而光照條件和營養(yǎng)鹽含量達(dá)到最佳平衡時(shí)往往形成浮游植物高生物量區(qū)[30-31]。韓秀榮[7]利用2002—2007年的調(diào)查資料統(tǒng)計(jì)得知葉綠素高值發(fā)生時(shí)(最大值可達(dá)70μg/L)的DIN、PO4-P以及濁度均維持在較低的水平,而Si O3-Si濃度約在20μmol/L、鹽度約30 PSU以及溫度約為20℃時(shí)葉綠素濃度最高。本文GAM模型中的PO4-P和Si O3-Si結(jié)果(PO4-P濃度較低,Si O3-Si濃度約22μmol/L時(shí)葉綠素表現(xiàn)為一高值)與其相一致,但與DIN和鹽度模擬結(jié)果存在差異(圖4),這可能與兩者研究區(qū)域和時(shí)間上的不同(韓秀榮[7]研究區(qū)域往南延伸,研究時(shí)間以春季為主)有關(guān)。本文引入的GAM模型分析方法在量化上具有較大的優(yōu)勢,在甄別影響浮游植物現(xiàn)存量的多環(huán)境因子不同作用方面具有很好的效果。

    長江口及鄰近海域近年來已成為赤潮頻發(fā)區(qū),春季赤潮藻類群以甲藻為主(米氏凱倫藻、鏈狀亞歷山大藻、具齒原甲藻等)[32],這與水體 PO4-P含量較低以及較高的DIN/P有關(guān)。李瑞香等[33]的船基圍隔實(shí)驗(yàn)結(jié)果表明低磷條件下具齒原甲藻生長比中肋骨條藻生長迅速且持續(xù)時(shí)間長;孫軍等[34]研究表明米氏凱倫藻在氮磷比為80∶1的條件下比生長率最快;呂頌輝和李英[35]研究結(jié)果顯示具齒原甲藻和鏈狀亞歷山大藻較之中肋骨條藻和尖刺擬菱形藻具有較高的細(xì)胞氮磷營養(yǎng)儲(chǔ)存能力,屬于k-策略者,有助于甲藻赤潮的長時(shí)間持續(xù)。本文研究結(jié)果亦表明春季葉綠素高值區(qū)在空間和季節(jié)比較上均具有較低的PO4-P含量和較高的DIN/P比值,為該海域春季甲藻赤潮頻發(fā)提供了有利條件。另一方面,夏季葉綠素高值區(qū)的Si/DIN在空間比較上比值較高(見3.1部分),這與夏季藻華種類主要是硅藻(優(yōu)勢種為中肋骨條藻和細(xì)長翼鼻狀藻[36])相一致。長期變化來看(近50年),DIN的升高和Si O3-Si的降低導(dǎo)致DIN/P和DIN/Si比值均呈升高的趨勢,也是近年來赤潮優(yōu)勢種由硅藻演變?yōu)榧自宓闹饕颍?0]。

    長江口及鄰近海域水體的高渾濁度使光照條件成為影響浮游植物生長的重要環(huán)境因子。本文觀測結(jié)果顯示,春、夏兩季葉綠素高值區(qū)的TSS濃度在時(shí)空比較上均顯著低值(見表1),說明較好的光照條件是浮游植物形成高生物量的一個(gè)重要條件。李雁賓[37]通過環(huán)境因子變異系數(shù)研究指出光照對硅藻赤潮的季節(jié)性生消影響最大,光照條件的增大有助于硅藻在與甲藻的競爭中占優(yōu);孫百曄[12]利用生長光照效應(yīng)分析及模擬證實(shí)了海水中的光照條件是限制赤潮進(jìn)一步向近岸方向擴(kuò)展的關(guān)鍵因素。但GAM模型分析結(jié)果顯示,TSS與葉綠素變化之間的直接相關(guān)性不顯著,說明光照條件對長江口及鄰近海域浮游植物現(xiàn)存量的控制作用不起決定性作用。

    本文觀測結(jié)果和模型分析顯示,長江口及鄰近海域浮游動(dòng)物生物量和葉綠素濃度之間相關(guān)性不顯著,這表明浮游動(dòng)物生物量在影響浮游植物現(xiàn)存量時(shí)空變化方面不是決定性因子?,F(xiàn)場實(shí)驗(yàn)結(jié)果顯示,微型浮游動(dòng)物在靠岸站位對浮游植物現(xiàn)存量的攝食壓力相對較高[13],而在水華爆發(fā)前后以及水華區(qū)和非水華區(qū)微型浮游動(dòng)物比攝食率各有不同[14]。另一方面,營養(yǎng)傳遞導(dǎo)致的浮游動(dòng)物生物量升高或下降相對浮游植物現(xiàn)存量變化存在一定的滯后性,往往導(dǎo)致兩者之間直接相關(guān)性不顯著。高的浮游植物現(xiàn)存量可為浮游動(dòng)物提供足夠食物,而浮游動(dòng)物生物量升高會(huì)反過來加劇對浮游植物的攝食,使浮游植物現(xiàn)存量下降;另一方面,浮游植物低生物量時(shí),可能是環(huán)境條件不適宜浮游生物的生長(此時(shí)浮游動(dòng)物生物量也低),也可能是較高的浮游動(dòng)物攝食作用所導(dǎo)致。此外,浮游動(dòng)物相對浮游植物一般具有較強(qiáng)的水平或垂直移動(dòng)能力,也可能影響兩者之間的直接聯(lián)系。因此,就單一或某幾個(gè)航次結(jié)果來看,很難找出浮游動(dòng)物和浮游植物現(xiàn)存量之間較好的相關(guān)關(guān)系。

    因此,在長江口及鄰近海域,營養(yǎng)鹽的上行效應(yīng)是浮游植物現(xiàn)存量的主要控制作用,而光照條件和浮游動(dòng)物生物量與浮游植物現(xiàn)存量雖然在時(shí)空分布上存在一定的聯(lián)系,但非決定性控制因素。

    [1]Carpenter S R,Kitchell J F,Hodgson J R.Cascading trophic interactions and lake productivity[J].BioScience,1985,35(10):634-639.

    [2]Mc Queen D J,Post J R,Mills E L.Trophic relationships in freshwater pelagic ecosystems[J].Canadian Journal of Fisheries and Aquatic Sciences,1986,43(8):1571-1581.

    [3]Leh man J T.Interacting growth and loss rates:the balance of top-down and bottom-up controls in plankton communities[J].Limnology and Oceanography,1991,36(8):1546-1554.

    [4]Sommer U.Scarcity of medium-sized phytoplankton in the northern Red Sea explained by strong bottom-up and weak top-down control[J].Marine Ecology Progress Series,2000,197:19-25.

    [5]Goericke R.Top-down control of phytoplankton biomass and community structure in the monsoonal Arabian Sea[J].Li mnology and Oceanography,2002,47(5):1307-1323.

    [6]Thompson P A,Bonham P I,Swadling km.Phytoplankton blooms in the Huon Estuary,Tasmania:top-down or bottom-up control?[J].Journal of Plankton Research,2008,30(7):735-753.

    [7]韓秀榮.長江口及鄰近海域浮游植物生長的多環(huán)境效應(yīng)因子影響解析研究[D].青島:中國海洋大學(xué),2009.

    [8]王保棟.長江口及鄰近海域富營養(yǎng)化狀況及其生態(tài)效應(yīng)[D].青島:中國海洋大學(xué),2006.

    [9]楊東方,王凡,高振會(huì),等.長江口理化因子影響初級生產(chǎn)力的探索:Ⅱ.磷不是長江口浮游植物生長的限制因子[J].海洋科學(xué)進(jìn)展,2006,24(1):97-107.

    [10]王江濤,曹婧.長江口海域近50a來營養(yǎng)鹽的變化及其對浮游植物群落演替的影響[J].海洋環(huán)境科學(xué),2012,31(3):310-315.

    [11]方濤.光照和N、P營養(yǎng)鹽的共同作用對長江口浮游植物生長的影響[D].上海:華東師范大學(xué),2008.

    [12]孫百曄.長江口及鄰近海域浮游植物生長的光照效應(yīng)研究[D].青島:中國海洋大學(xué),2008.

    [13]孫軍,劉東艷,王宗靈,等.春季赤潮頻發(fā)期東海微型浮游動(dòng)物攝食研究[J].應(yīng)用生態(tài)學(xué)報(bào),2003,14(7):1073-1080.

    [14]孫軍,宋書群.東海春季水華期浮游植物生長與微型浮游動(dòng)物攝食[J].生態(tài)學(xué)報(bào),2009,29(12):6429-6438.

    [15]Hastie T,Tibshirani R.Generalized additive models[J].Statistical Science,1986,1(3):297-310.

    [16]Lamon III E C,Reckhow K H,Havens K E.Using generalized additive models for prediction of chlorophyllain Lake Okeechobee,F(xiàn)lorida[J].Lakes &Reservoirs:Research and Management,1996,2(1/2):37-46.

    [17]Raitsos D E,Korres G,Triantafyllou G,et al.Assessing chlorophyll variability in relation to the environmental regime in Pagasitikos Gulf,Greece[J].Journal of Marine Systems,2012,94:S16-S22.

    [18]Tew-Kai E,Marsac F.Patterns of variability of sea surface chlorophyll in the Mozambique Channel:A quantitative approach[J].Journal of Marine Systems,2009,77(1/2):77-88.

    [19]Chen B Z,Liu H B,Huang B Q.Environmental controlling mechanisms on bacterial abundance in the South China Sea inferred from generalized additive models(GAMs)[J].Journal of Sea Research,2012,72:69-76.

    [20]Chen B Z,Landry M R,Huang B Q,et al.Does warming enhance the effect of microzooplankton grazing on marine phytoplankton in the ocean?[J].Li mnology and Oceanography,2012,57(2):519-526.

    [21]鄭波,陳新軍,李綱.GL M和GA M模型研究東黃海鮐資源漁場與環(huán)境因子的關(guān)系[J].水產(chǎn)學(xué)報(bào),2008,32(3):379-386.

    [22]朱國平.基于廣義可加模型研究時(shí)間和環(huán)境因子對南極半島北部南極磷蝦漁場的影響[J].水產(chǎn)學(xué)報(bào),2012,36(12):1863-1871.

    [23]GB/T 12763.6-2007,海洋調(diào)查規(guī)范 第6部分:海洋生物調(diào)查[S].北京:中國標(biāo)準(zhǔn)出版社,2008.

    [24]GB/T 12763.4-2007,海洋調(diào)查規(guī)范 第4部分:海水化學(xué)要素調(diào)查[S].北京:中國標(biāo)準(zhǔn)出版社,2008.

    [25]GB 17378.4-2007,海洋監(jiān)測規(guī)范 第4部分:海水分析[S].北京:中國標(biāo)準(zhǔn)出版社,2008.

    [26]Wood S N.Stable and efficient multiple smoothing parameter estimation for generalized additive models[J].Journal of the American Statistical Association,2004,99(467):673-686.

    [27]Campbell J W.The lognor mal distribution as a model for bio-optical variability in the sea[J].Journal of Geophysical Research:Oceans,1995,100(C7):13237-13254.

    [28]孫軍,宋書群,欒青杉.長江口水域浮游植物生物量及結(jié)構(gòu)特點(diǎn)[M]//俞志明,沈志良,陳亞瞿,等.長江口水域富營養(yǎng)化.北京:科學(xué)出版社,2011:184-203.

    [29]Banse K.Grazing and zooplankton production as key controls of phytoplankton production in the open ocean[J].Oceanography,1994,7(1):13-20.

    [30]寧修仁,史君賢,蔡昱明,等.長江口和杭州灣海域生物生產(chǎn)力鋒面及其生態(tài)學(xué)效應(yīng)[J].海洋學(xué)報(bào),2004,26(6):96-106.

    [31]孫軍,宋書群,欒青杉.長江口水域浮游植物關(guān)鍵過程及其對富營養(yǎng)化的影響[M]//俞志明,沈志良,陳亞瞿,等.長江口水域富營養(yǎng)化.北京:科學(xué)出版社,2011:336-349.

    [32]周名江,朱明遠(yuǎn).“我國近海有害赤潮發(fā)生的生態(tài)學(xué)、海洋學(xué)機(jī)制及預(yù)測防治”研究進(jìn)展[J].地球科學(xué)進(jìn)展,2006,21(7):673-679.

    [33]李瑞香,朱明遠(yuǎn),王宗靈,等.東海兩種赤潮生物種間競爭的圍隔實(shí)驗(yàn)[J].應(yīng)用生態(tài)學(xué)報(bào),2003,14(7):1049-1054.

    [34]孫軍,劉東艷,陳宗濤,等.不同氮磷比率對青島大扁藻、新月柱鞘藻和米氏凱倫藻生長影響及其生存策略研究[J].應(yīng)用生態(tài)學(xué)報(bào),2004,15(11):2122-2126.

    [35]呂頌輝,李英.我國東海4種赤潮藻的細(xì)胞氮磷營養(yǎng)儲(chǔ)存能力對比[J].過程工程學(xué)報(bào),2006,6(3):439-444.

    [36]孫軍,宋書群,欒青杉.長江口水域浮游植物群集組成結(jié)構(gòu)與特點(diǎn)[M]//俞志明,沈志良,陳亞瞿,等.長江口水域富營養(yǎng)化.北京:科學(xué)出版社,2011:144-184.

    [37]李雁賓.長江口及鄰近海域季節(jié)性赤潮生消過程控制機(jī)理研究[D].青島:中國海洋大學(xué),2008.

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    电影成人av| 波多野结衣av一区二区av| 国产精品久久久久久久久免| 亚洲欧美中文字幕日韩二区| 丰满少妇做爰视频| 好男人视频免费观看在线| 日日摸夜夜添夜夜爱| 大香蕉久久成人网| 天堂中文最新版在线下载| 一级毛片黄色毛片免费观看视频| 亚洲成人av在线免费| 久久久久精品人妻al黑| 亚洲精品美女久久久久99蜜臀 | 最新的欧美精品一区二区| 交换朋友夫妻互换小说| 麻豆av在线久日| 欧美人与性动交α欧美精品济南到 | 欧美黄色片欧美黄色片| 国产白丝娇喘喷水9色精品| 久久人妻熟女aⅴ| 边亲边吃奶的免费视频| 国产精品成人在线| av在线观看视频网站免费| 青草久久国产| 男人操女人黄网站| 久久精品国产自在天天线| 精品国产一区二区久久| 亚洲国产av影院在线观看| 菩萨蛮人人尽说江南好唐韦庄| 性少妇av在线| 各种免费的搞黄视频| 99热网站在线观看| 成人国语在线视频| 视频区图区小说| 日产精品乱码卡一卡2卡三| 亚洲人成电影观看| 成人毛片a级毛片在线播放| 欧美精品人与动牲交sv欧美| 亚洲av国产av综合av卡| 天堂中文最新版在线下载| 性高湖久久久久久久久免费观看| 熟女av电影| 国产精品香港三级国产av潘金莲 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美色中文字幕在线| 日韩 亚洲 欧美在线| 人妻人人澡人人爽人人| 亚洲欧美清纯卡通| 在线 av 中文字幕| 午夜激情av网站| 韩国av在线不卡| 国产精品一二三区在线看| 成人二区视频| 啦啦啦在线免费观看视频4| 国产亚洲最大av| 2018国产大陆天天弄谢| 日韩av免费高清视频| 一级毛片电影观看| 亚洲情色 制服丝袜| 黄色怎么调成土黄色| 久久国产亚洲av麻豆专区| 卡戴珊不雅视频在线播放| 精品卡一卡二卡四卡免费| 亚洲,一卡二卡三卡| 国产免费视频播放在线视频| 99九九在线精品视频| 欧美人与性动交α欧美软件| 欧美成人午夜精品| 性色av一级| 女人被躁到高潮嗷嗷叫费观| 丝瓜视频免费看黄片| av在线老鸭窝| 亚洲精品久久成人aⅴ小说| a级片在线免费高清观看视频| 十八禁网站网址无遮挡| 国产免费现黄频在线看| 丝袜美腿诱惑在线| 成人毛片a级毛片在线播放| 青草久久国产| 亚洲伊人久久精品综合| xxx大片免费视频| 欧美精品av麻豆av| 亚洲色图综合在线观看| 国语对白做爰xxxⅹ性视频网站| 国产老妇伦熟女老妇高清| 少妇熟女欧美另类| 欧美 日韩 精品 国产| 免费黄频网站在线观看国产| 高清在线视频一区二区三区| 精品视频人人做人人爽| 亚洲第一青青草原| 久久精品国产亚洲av涩爱| 人妻一区二区av| 久久久久久久久久人人人人人人| 亚洲欧美一区二区三区黑人 | 国产1区2区3区精品| 精品少妇久久久久久888优播| 啦啦啦视频在线资源免费观看| 久久精品国产综合久久久| 老司机影院成人| av女优亚洲男人天堂| 另类精品久久| 中文字幕亚洲精品专区| 国产日韩欧美亚洲二区| 一个人免费看片子| 少妇被粗大的猛进出69影院| 中文字幕人妻丝袜一区二区 | av在线老鸭窝| 99re6热这里在线精品视频| 午夜福利在线免费观看网站| 男女下面插进去视频免费观看| 亚洲av中文av极速乱| 26uuu在线亚洲综合色| 国产成人免费观看mmmm| av网站在线播放免费| 黄片播放在线免费| 精品卡一卡二卡四卡免费| 人妻少妇偷人精品九色| 免费看av在线观看网站| 新久久久久国产一级毛片| 蜜桃国产av成人99| 香蕉精品网在线| 国产乱人偷精品视频| 久久精品aⅴ一区二区三区四区 | 永久网站在线| 母亲3免费完整高清在线观看 | 成人国语在线视频| 丰满乱子伦码专区| 在线观看美女被高潮喷水网站| 午夜日韩欧美国产| 777久久人妻少妇嫩草av网站| 日韩一区二区三区影片| 亚洲伊人色综图| 亚洲久久久国产精品| 边亲边吃奶的免费视频| 国产成人91sexporn| 黑人欧美特级aaaaaa片| 大片免费播放器 马上看| 久久午夜福利片| 中文乱码字字幕精品一区二区三区| 一区二区三区激情视频| 亚洲成人一二三区av| 日本91视频免费播放| 18+在线观看网站| 老司机影院毛片| 永久免费av网站大全| 亚洲熟女精品中文字幕| 一级片'在线观看视频| 国产精品免费大片| 美女xxoo啪啪120秒动态图| 嫩草影院入口| 精品少妇一区二区三区视频日本电影 | 免费高清在线观看日韩| 久久久亚洲精品成人影院| 国产精品秋霞免费鲁丝片| 午夜91福利影院| a级毛片黄视频| 三上悠亚av全集在线观看| 亚洲第一av免费看| 久久精品国产亚洲av涩爱| 国产福利在线免费观看视频| 一区在线观看完整版| 亚洲一级一片aⅴ在线观看| 国产综合精华液| www日本在线高清视频| 国产精品久久久久久久久免| 国产免费又黄又爽又色| 亚洲成av片中文字幕在线观看 | 免费播放大片免费观看视频在线观看| 一级片'在线观看视频| 少妇猛男粗大的猛烈进出视频| 国产黄色视频一区二区在线观看| av线在线观看网站| 宅男免费午夜| 免费观看av网站的网址| 亚洲天堂av无毛| 国产激情久久老熟女| 大片免费播放器 马上看| 麻豆乱淫一区二区| 国产亚洲av片在线观看秒播厂| 亚洲精品一二三| 捣出白浆h1v1| 久久久亚洲精品成人影院| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产色婷婷电影| 一本久久精品| 成人二区视频| av在线app专区| 不卡视频在线观看欧美| 亚洲激情五月婷婷啪啪| 亚洲精品乱久久久久久| 免费观看性生交大片5| 国产成人aa在线观看| 国产在线一区二区三区精| 99国产精品免费福利视频| 久久99一区二区三区| 亚洲美女视频黄频| 国产亚洲最大av| 超碰成人久久| 999久久久国产精品视频| 免费黄频网站在线观看国产| 亚洲男人天堂网一区| 亚洲av电影在线进入| 成人毛片a级毛片在线播放| 如日韩欧美国产精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽av亚洲精品天堂| 中文乱码字字幕精品一区二区三区| 亚洲国产成人一精品久久久| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜制服| 一区二区三区精品91| 久久久精品国产亚洲av高清涩受| 欧美精品国产亚洲| 欧美日韩国产mv在线观看视频| 97精品久久久久久久久久精品| 久久久久久人妻| 免费少妇av软件| 久久精品国产亚洲av涩爱| 午夜日本视频在线| 欧美国产精品一级二级三级| 中文字幕精品免费在线观看视频| 天天躁夜夜躁狠狠躁躁| 日韩视频在线欧美| 久久久久久久久久久免费av| 国产精品无大码| 成人国产麻豆网| 久久热在线av| 国产 精品1| 国产精品麻豆人妻色哟哟久久| av有码第一页| 波多野结衣一区麻豆| 男女边吃奶边做爰视频| 国语对白做爰xxxⅹ性视频网站| 熟妇人妻不卡中文字幕| 日韩在线高清观看一区二区三区| av.在线天堂| 日本vs欧美在线观看视频| 王馨瑶露胸无遮挡在线观看| av线在线观看网站| 欧美最新免费一区二区三区| av片东京热男人的天堂| 王馨瑶露胸无遮挡在线观看| 新久久久久国产一级毛片| 亚洲欧美精品自产自拍| 在线观看三级黄色| 9热在线视频观看99| 久久精品aⅴ一区二区三区四区 | 男女免费视频国产| 精品久久久久久电影网| 黄色毛片三级朝国网站| 国产精品久久久久久精品古装| 两个人免费观看高清视频| 99精国产麻豆久久婷婷| 亚洲美女视频黄频| 亚洲色图综合在线观看| 热re99久久精品国产66热6| 久久久精品免费免费高清| 人成视频在线观看免费观看| 久久久久视频综合| 久久精品国产综合久久久| 国产欧美日韩一区二区三区在线| 深夜精品福利| 91午夜精品亚洲一区二区三区| 久久久久久伊人网av| 午夜91福利影院| 国产成人精品一,二区| 成人亚洲欧美一区二区av| 久久免费观看电影| 日本猛色少妇xxxxx猛交久久| a级毛片黄视频| 男人添女人高潮全过程视频| 免费观看av网站的网址| 在线观看免费高清a一片| 在线 av 中文字幕| 成年av动漫网址| 国产精品.久久久| av在线观看视频网站免费| 国产日韩一区二区三区精品不卡| 曰老女人黄片| 18禁观看日本| 国产精品.久久久| 26uuu在线亚洲综合色| 亚洲国产精品一区三区| 免费看av在线观看网站| 九九爱精品视频在线观看| 1024香蕉在线观看| 99国产综合亚洲精品| 日韩精品有码人妻一区| 免费观看性生交大片5| 久热久热在线精品观看| 十八禁网站网址无遮挡| 欧美人与性动交α欧美精品济南到 | 伊人亚洲综合成人网| 亚洲综合色网址| 一级毛片我不卡| 丁香六月天网| 欧美激情高清一区二区三区 | 久久免费观看电影| 亚洲一级一片aⅴ在线观看| 2021少妇久久久久久久久久久| 亚洲欧美色中文字幕在线| 亚洲国产精品一区二区三区在线| 大码成人一级视频| 亚洲欧美清纯卡通| 少妇精品久久久久久久| 最近中文字幕2019免费版| 成人国产av品久久久| 色婷婷久久久亚洲欧美| 波多野结衣av一区二区av| 最近最新中文字幕大全免费视频 | 免费黄网站久久成人精品| 国产探花极品一区二区| 哪个播放器可以免费观看大片| 久久青草综合色| 女人被躁到高潮嗷嗷叫费观| 精品人妻熟女毛片av久久网站| 国产成人精品无人区| 国产成人91sexporn| 欧美日韩视频高清一区二区三区二| 国产片特级美女逼逼视频| 午夜福利乱码中文字幕| 亚洲av欧美aⅴ国产| 91久久精品国产一区二区三区| 中文乱码字字幕精品一区二区三区| 亚洲三级黄色毛片| 亚洲国产精品一区二区三区在线| 国产精品三级大全| 水蜜桃什么品种好| 欧美激情 高清一区二区三区| 美女大奶头黄色视频| 国产精品欧美亚洲77777| 亚洲一区中文字幕在线| 午夜免费观看性视频| 日韩一本色道免费dvd| 中文字幕最新亚洲高清| 免费看av在线观看网站| 性高湖久久久久久久久免费观看| 国产精品成人在线| av一本久久久久| 免费看av在线观看网站| 黄色配什么色好看| 亚洲国产av新网站| 日韩欧美精品免费久久| 精品亚洲成a人片在线观看| 最近最新中文字幕大全免费视频 | 又大又黄又爽视频免费| 国产免费视频播放在线视频| 咕卡用的链子| 欧美日韩视频精品一区| 又粗又硬又长又爽又黄的视频| 亚洲国产欧美在线一区| 久久久久人妻精品一区果冻| 91精品伊人久久大香线蕉| 婷婷色av中文字幕| 伊人久久国产一区二区| 国产在线视频一区二区| 午夜日本视频在线| 亚洲国产av新网站| 综合色丁香网| 亚洲国产色片| 1024香蕉在线观看| 99国产综合亚洲精品| 26uuu在线亚洲综合色| 亚洲国产精品一区三区| 黄色视频在线播放观看不卡| 久久久a久久爽久久v久久| 国产精品.久久久| 青草久久国产| 男女高潮啪啪啪动态图| 超碰成人久久| 十八禁高潮呻吟视频| 久久女婷五月综合色啪小说| 久久99精品国语久久久| 一级毛片电影观看| 亚洲国产精品成人久久小说| 波多野结衣av一区二区av| 亚洲经典国产精华液单| 国产成人免费观看mmmm| 性色avwww在线观看| 久久人人97超碰香蕉20202| 精品少妇黑人巨大在线播放| 亚洲精品久久久久久婷婷小说| 男女高潮啪啪啪动态图| a 毛片基地| 国产色婷婷99| 一级毛片黄色毛片免费观看视频| 成人黄色视频免费在线看| 高清不卡的av网站| 亚洲精品国产av蜜桃| 精品人妻熟女毛片av久久网站| 免费观看在线日韩| av免费在线看不卡| 丰满少妇做爰视频| www.自偷自拍.com| 国产精品二区激情视频| 精品福利永久在线观看| 国产精品欧美亚洲77777| 亚洲精品第二区| 亚洲精华国产精华液的使用体验| 一级,二级,三级黄色视频| 国产人伦9x9x在线观看 | 人体艺术视频欧美日本| 99久久中文字幕三级久久日本| 99久久综合免费| 精品99又大又爽又粗少妇毛片| 久久久久国产网址| 少妇人妻 视频| 午夜日韩欧美国产| 伊人久久国产一区二区| 99精国产麻豆久久婷婷| 色吧在线观看| 国产精品熟女久久久久浪| 一个人免费看片子| 国产精品.久久久| 亚洲国产av影院在线观看| 国语对白做爰xxxⅹ性视频网站| 精品人妻一区二区三区麻豆| 丝袜美腿诱惑在线| 一本—道久久a久久精品蜜桃钙片| 99久久精品国产国产毛片| 成人午夜精彩视频在线观看| 国产精品偷伦视频观看了| 久久99热这里只频精品6学生| 纯流量卡能插随身wifi吗| 欧美最新免费一区二区三区| 各种免费的搞黄视频| 一级a爱视频在线免费观看| 免费在线观看完整版高清| 26uuu在线亚洲综合色| 精品国产国语对白av| 精品视频人人做人人爽| 最新的欧美精品一区二区| 亚洲精品久久午夜乱码| 精品亚洲乱码少妇综合久久| 日韩av不卡免费在线播放| 人人妻人人添人人爽欧美一区卜| 国产在线一区二区三区精| 久久精品国产亚洲av高清一级| 精品一区在线观看国产| 国产精品蜜桃在线观看| 亚洲经典国产精华液单| 啦啦啦在线观看免费高清www| 亚洲精品成人av观看孕妇| 欧美激情极品国产一区二区三区| 桃花免费在线播放| 一级a爱视频在线免费观看| 欧美最新免费一区二区三区| 亚洲内射少妇av| 久久久久久免费高清国产稀缺| 9热在线视频观看99| tube8黄色片| 视频在线观看一区二区三区| 黄色毛片三级朝国网站| 久久精品国产亚洲av天美| 国产熟女欧美一区二区| 熟女少妇亚洲综合色aaa.| 超碰成人久久| 亚洲视频免费观看视频| 亚洲欧美成人综合另类久久久| 一本大道久久a久久精品| 亚洲欧美成人精品一区二区| 青草久久国产| 丝袜美足系列| 亚洲美女搞黄在线观看| 亚洲欧美精品自产自拍| 国产日韩一区二区三区精品不卡| 久久精品国产鲁丝片午夜精品| 午夜福利网站1000一区二区三区| 久久久久久久精品精品| 精品国产一区二区三区四区第35| 超色免费av| 欧美 亚洲 国产 日韩一| 在线精品无人区一区二区三| 亚洲国产欧美在线一区| 亚洲国产欧美日韩在线播放| av网站免费在线观看视频| 精品国产国语对白av| 国产又色又爽无遮挡免| 免费高清在线观看视频在线观看| 久久精品久久精品一区二区三区| 捣出白浆h1v1| 春色校园在线视频观看| 亚洲av国产av综合av卡| 蜜桃国产av成人99| 成人黄色视频免费在线看| 亚洲伊人色综图| 精品久久蜜臀av无| 免费看不卡的av| 国产av国产精品国产| 日本免费在线观看一区| 亚洲一区中文字幕在线| 超色免费av| av线在线观看网站| 男女边吃奶边做爰视频| 九草在线视频观看| 天天躁日日躁夜夜躁夜夜| 美女脱内裤让男人舔精品视频| 不卡av一区二区三区| 18在线观看网站| 九九爱精品视频在线观看| 久久久久视频综合| 亚洲欧美成人精品一区二区| 亚洲国产欧美网| 久久99一区二区三区| 亚洲五月色婷婷综合| 午夜av观看不卡| 亚洲国产精品国产精品| 多毛熟女@视频| 日本黄色日本黄色录像| 99久久精品国产国产毛片| 国产黄色免费在线视频| av线在线观看网站| 又粗又硬又长又爽又黄的视频| 国产乱人偷精品视频| 丝袜人妻中文字幕| 99久久综合免费| 日韩在线高清观看一区二区三区| 亚洲av福利一区| 亚洲成人一二三区av| 久久精品国产亚洲av高清一级| 大片免费播放器 马上看| 国产不卡av网站在线观看| 女人精品久久久久毛片| 18禁动态无遮挡网站| 寂寞人妻少妇视频99o| 国产免费视频播放在线视频| 亚洲av在线观看美女高潮| 超色免费av| 在线 av 中文字幕| 国产成人精品在线电影| 亚洲欧洲精品一区二区精品久久久 | 亚洲伊人色综图| 日本午夜av视频| 国产精品国产三级国产专区5o| 人人妻人人澡人人爽人人夜夜| 亚洲av电影在线观看一区二区三区| 老女人水多毛片| 99热国产这里只有精品6| 成年人午夜在线观看视频| videosex国产| 亚洲精品久久午夜乱码| 午夜免费男女啪啪视频观看| av又黄又爽大尺度在线免费看| 飞空精品影院首页| 精品福利永久在线观看| 天天躁日日躁夜夜躁夜夜| 日韩,欧美,国产一区二区三区| √禁漫天堂资源中文www| 老女人水多毛片| 久久精品国产自在天天线| a级毛片黄视频| 午夜福利在线观看免费完整高清在| 美女福利国产在线| 久久久久精品人妻al黑| 国产视频首页在线观看| 亚洲国产看品久久| 精品一区二区免费观看| 亚洲精品国产av成人精品| 久久这里有精品视频免费| 18禁裸乳无遮挡动漫免费视频| av女优亚洲男人天堂| 99国产精品免费福利视频| 黑人欧美特级aaaaaa片| 亚洲精品一二三| 国产视频首页在线观看| 亚洲精品久久成人aⅴ小说| 成人国语在线视频| 欧美精品高潮呻吟av久久| 日韩av免费高清视频| 亚洲第一青青草原| 午夜免费鲁丝| 性高湖久久久久久久久免费观看| 亚洲婷婷狠狠爱综合网| 中文字幕另类日韩欧美亚洲嫩草| 老司机影院成人| 亚洲精品自拍成人| 丝袜在线中文字幕| 日韩伦理黄色片| 亚洲伊人色综图| 欧美精品国产亚洲| 免费看不卡的av| 久久97久久精品| 亚洲国产日韩一区二区| 啦啦啦视频在线资源免费观看| 性高湖久久久久久久久免费观看| av国产久精品久网站免费入址| 熟女电影av网| 丝袜脚勾引网站| 成人毛片a级毛片在线播放| a 毛片基地| 丝袜人妻中文字幕| 欧美亚洲日本最大视频资源| 精品国产乱码久久久久久小说| 丝瓜视频免费看黄片| 七月丁香在线播放| 午夜免费观看性视频| 美女xxoo啪啪120秒动态图| 精品人妻熟女毛片av久久网站| 精品少妇内射三级| 人体艺术视频欧美日本| 99久国产av精品国产电影| 母亲3免费完整高清在线观看 | 精品亚洲乱码少妇综合久久| 精品久久蜜臀av无| 成年动漫av网址| 国产精品一区二区在线不卡| 国产探花极品一区二区| 亚洲精品av麻豆狂野| 久久女婷五月综合色啪小说| 日本欧美国产在线视频| 男女边摸边吃奶| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕免费大全7| 亚洲av在线观看美女高潮|