• <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)主要是指浮游動物通過攝食作用影響浮游植物生物量的累積。上下行效應(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)主要受浮游動物的下行控制作用;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ā)生的重要因素。而浮游動物對浮游植物的下行控制研究目前主要集中在研究區(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]、浮游動物[20]、漁業(yè)資源[21-22]等與環(huán)境因子間關(guān)系的研究中引入了GAM模型分析方法。

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

    2 材料與方法

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

    長江口及鄰近海域現(xiàn)場調(diào)查站位分布如圖1所示。共執(zhí)行3個航次,調(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自動分析儀測定(其中DIN為NO3-N、NO2-N和NH4-N總和),懸浮物采用重量法測定,浮游動物生物量指淺水Ⅰ型浮游生物網(wǎng)采樣品的濕重生物量(含膠質(zhì)浮游動物)。

    圖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和浮游動物生物量(ZB)。由于葉綠素值屬于高度右偏態(tài)分布,分析時(shí)取log10(Chl)轉(zhuǎn)換以達(dá)到近似正態(tài)分布[27],同理對于浮游動物生物量和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è)形成兩個高值區(qū)(>5.0μg/L,圖2b);秋季,葉綠素濃度普遍較低,僅在3個站位稍高于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)顯示:春季,與整個調(diào)查區(qū)相比,葉綠素高值區(qū)的溫度平均值相同,鹽度稍高,TSS顯著低值,而營養(yǎng)鹽(DIN、PO4-P和Si O3-Si)也相對較低,且PO4-P和Si O3-Si濃度在3個季節(jié)中最低;夏季,葉綠素高值區(qū)的溫度相對整個調(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

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

    3.2 GAM模型結(jié)果

    GAM模型結(jié)果(見表2)顯示:溫度(SST)、鹽度(SSS)、懸浮物(TSS)、營養(yǎng)鹽(DIN、PO4和Si O3)和浮游動物生物量(ZB)變化對葉綠素時(shí)空變化的貢獻(xiàn)合計(jì)可達(dá)74.6%(模型1),依次去除統(tǒng)計(jì)學(xué)意義上不顯著(p>0.05)的TSS(模型2)、浮游動物生物量(模型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、浮游動物生物量和溫度三因子的直接相關(guān)性不顯著。

    圖3 2010—2011年度長江口及鄰近海域浮游動物生物量平面分布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á)到一個峰值,然后降低且在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ò)散與混合、水平移動、沉降損失與被攝食消耗等因素[29]。浮游植物的生長率主要受營養(yǎng)鹽和光照的上行效應(yīng)影響,而被攝食消耗則以浮游動物的下行效應(yīng)為主。本文觀測結(jié)果發(fā)現(xiàn),長江口及鄰近海域3個季節(jié)葉綠素高值區(qū)營養(yǎng)鹽相對全區(qū)平均值較低,夏春兩季葉綠素高值區(qū)(也是全年葉綠素高值)的水體TSS相對較低,夏季葉綠素高值區(qū)鹽度較低,而溫度在葉綠素高值區(qū)差別不大;浮游動物生物量高值區(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)鹽和光照條件等因子的共同作用),而浮游動物生物量和溫度等因子直接影響不顯著。

    圖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)儲存能力,屬于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),說明較好的光照條件是浮游植物形成高生物量的一個重要條件。李雁賓[37]通過環(huán)境因子變異系數(shù)研究指出光照對硅藻赤潮的季節(jié)性生消影響最大,光照條件的增大有助于硅藻在與甲藻的競爭中占優(yōu);孫百曄[12]利用生長光照效應(yīng)分析及模擬證實(shí)了海水中的光照條件是限制赤潮進(jìn)一步向近岸方向擴(kuò)展的關(guān)鍵因素。但GAM模型分析結(jié)果顯示,TSS與葉綠素變化之間的直接相關(guān)性不顯著,說明光照條件對長江口及鄰近海域浮游植物現(xiàn)存量的控制作用不起決定性作用。

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

    因此,在長江口及鄰近海域,營養(yǎng)鹽的上行效應(yīng)是浮游植物現(xiàn)存量的主要控制作用,而光照條件和浮游動物生物量與浮游植物現(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]楊東方,王凡,高振會,等.長江口理化因子影響初級生產(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ā)期東海微型浮游動物攝食研究[J].應(yīng)用生態(tài)學(xué)報(bào),2003,14(7):1073-1080.

    [14]孫軍,宋書群.東海春季水華期浮游植物生長與微型浮游動物攝食[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)儲存能力對比[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| 人成视频在线观看免费观看| 69精品国产乱码久久久| 日本一区二区免费在线视频| 久久影院123| 高清不卡的av网站| 精品福利永久在线观看| 亚洲精品第二区| 国产黄频视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 校园人妻丝袜中文字幕| 97人妻天天添夜夜摸| 亚洲三区欧美一区| 亚洲五月色婷婷综合| 女的被弄到高潮叫床怎么办| 97在线人人人人妻| 精品亚洲成国产av| 天堂中文最新版在线下载| 午夜福利网站1000一区二区三区| 黄片无遮挡物在线观看| 欧美精品av麻豆av| 色综合欧美亚洲国产小说| 又大又爽又粗| 亚洲国产毛片av蜜桃av| 曰老女人黄片| 妹子高潮喷水视频| 老汉色av国产亚洲站长工具| 欧美激情高清一区二区三区 | 亚洲人成77777在线视频| 国产成人精品无人区| 亚洲av男天堂| 丁香六月欧美| 超碰97精品在线观看| 国产精品国产av在线观看| 美女扒开内裤让男人捅视频| 大香蕉久久成人网| 国产男人的电影天堂91| 女性生殖器流出的白浆| 日韩精品有码人妻一区| 色视频在线一区二区三区| 久久99热这里只频精品6学生| 啦啦啦中文免费视频观看日本| 色播在线永久视频| 精品一区在线观看国产| 这个男人来自地球电影免费观看 | 韩国高清视频一区二区三区| 久久人人97超碰香蕉20202| 国产在线视频一区二区| 这个男人来自地球电影免费观看 | 在线亚洲精品国产二区图片欧美| 精品一品国产午夜福利视频| 波多野结衣av一区二区av| 午夜福利在线免费观看网站| a级毛片在线看网站| 日本一区二区免费在线视频| 久久影院123| 女人爽到高潮嗷嗷叫在线视频| 国产免费视频播放在线视频| 另类精品久久| 少妇被粗大的猛进出69影院| 亚洲一区二区三区欧美精品| 欧美国产精品一级二级三级| 欧美精品亚洲一区二区| 精品一区二区三卡| 日韩av不卡免费在线播放| 日本一区二区免费在线视频| 1024视频免费在线观看| 一级毛片 在线播放| 久久 成人 亚洲| 97人妻天天添夜夜摸| 人人妻,人人澡人人爽秒播 | 中文字幕亚洲精品专区| 国产精品久久久av美女十八| 国产福利在线免费观看视频| bbb黄色大片| 亚洲一级一片aⅴ在线观看| 精品第一国产精品| 又粗又硬又长又爽又黄的视频| 69精品国产乱码久久久| 国产片内射在线| 亚洲国产精品一区三区| 欧美国产精品va在线观看不卡| 国产日韩欧美视频二区| 日韩,欧美,国产一区二区三区| 最近最新中文字幕大全免费视频 | 男人爽女人下面视频在线观看| 99久久综合免费| 美女福利国产在线| 老鸭窝网址在线观看| 亚洲国产日韩一区二区| 国产毛片在线视频| 国产淫语在线视频| 国产av一区二区精品久久| 欧美日韩亚洲综合一区二区三区_| 观看av在线不卡| 多毛熟女@视频| 久久精品国产综合久久久| xxxhd国产人妻xxx| 亚洲国产欧美日韩在线播放| 1024香蕉在线观看| 男女免费视频国产| 观看av在线不卡| 久久久久人妻精品一区果冻| 欧美日韩精品网址| 考比视频在线观看| 日韩人妻精品一区2区三区| 午夜老司机福利片| 美国免费a级毛片| av女优亚洲男人天堂| 少妇人妻精品综合一区二区| 国产伦理片在线播放av一区| 咕卡用的链子| 色播在线永久视频| 久久99热这里只频精品6学生| 免费在线观看视频国产中文字幕亚洲 | 久久久精品区二区三区| 高清视频免费观看一区二区| 亚洲一区二区三区欧美精品| 在线观看国产h片| 最近手机中文字幕大全| 中文欧美无线码| 美女主播在线视频| 九色亚洲精品在线播放| 人妻人人澡人人爽人人| 国产精品免费视频内射| 丁香六月天网| 久久久欧美国产精品| 久久婷婷青草| 成年人免费黄色播放视频| 国产伦人伦偷精品视频| 日韩av在线免费看完整版不卡| 亚洲人成电影观看| 街头女战士在线观看网站| 最新的欧美精品一区二区| 日韩大片免费观看网站| 日韩av免费高清视频| 男女国产视频网站| 日日撸夜夜添| 少妇人妻 视频| 日韩,欧美,国产一区二区三区| 国产欧美日韩综合在线一区二区| 少妇猛男粗大的猛烈进出视频| 亚洲五月色婷婷综合| 亚洲在久久综合| 免费高清在线观看视频在线观看| 丰满迷人的少妇在线观看| 国产精品久久久久久精品古装| a级片在线免费高清观看视频| 国产精品一区二区在线观看99| 99久久精品国产亚洲精品| 最近中文字幕高清免费大全6| 国产老妇伦熟女老妇高清| av视频免费观看在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久久久精品国产欧美久久久 | 亚洲在久久综合| 久久久久久久久久久免费av| 看免费成人av毛片| 日韩中文字幕欧美一区二区 | 国产精品久久久久久精品电影小说| 制服人妻中文乱码| 最近中文字幕高清免费大全6| 午夜日韩欧美国产| 欧美国产精品va在线观看不卡| 国产日韩欧美视频二区| 国产一卡二卡三卡精品 | 另类亚洲欧美激情| 久久久久久久大尺度免费视频| 日本91视频免费播放| avwww免费| 亚洲成人国产一区在线观看 | 中文字幕av电影在线播放| 欧美人与性动交α欧美软件| 9热在线视频观看99| 女人爽到高潮嗷嗷叫在线视频| 9热在线视频观看99| 久久狼人影院| 免费人妻精品一区二区三区视频| 亚洲精品av麻豆狂野| 啦啦啦在线观看免费高清www| 最近的中文字幕免费完整| 欧美97在线视频| 免费人妻精品一区二区三区视频| av在线播放精品| 久久久久精品性色| 亚洲国产精品999| 亚洲精品自拍成人| 亚洲成人手机| 黄色 视频免费看| 欧美日韩亚洲综合一区二区三区_| 日韩制服骚丝袜av| 亚洲精品av麻豆狂野| 午夜福利视频在线观看免费| 亚洲av电影在线进入| 婷婷色综合大香蕉| 一边摸一边抽搐一进一出视频| 亚洲精品成人av观看孕妇| avwww免费| 国产高清不卡午夜福利| 国产高清国产精品国产三级| 久久人人爽人人片av| 国产成人精品在线电影| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美精品综合一区二区三区| 亚洲国产精品999| 69精品国产乱码久久久| 少妇人妻精品综合一区二区| 亚洲精品一二三| 日韩一卡2卡3卡4卡2021年| 女人爽到高潮嗷嗷叫在线视频| 中文字幕亚洲精品专区| 成人亚洲欧美一区二区av| 国产麻豆69| 老司机在亚洲福利影院| 国产毛片在线视频| 国产精品无大码| 免费av中文字幕在线| 黄色视频在线播放观看不卡| 18禁国产床啪视频网站| 国产成人欧美在线观看 | 精品亚洲乱码少妇综合久久| 亚洲成人av在线免费| 欧美精品av麻豆av| 乱人伦中国视频| 亚洲国产欧美一区二区综合| 午夜久久久在线观看| 亚洲欧美中文字幕日韩二区| 看十八女毛片水多多多| 国产精品亚洲av一区麻豆 | 久久鲁丝午夜福利片| 91aial.com中文字幕在线观看| 一区二区三区乱码不卡18| 日本91视频免费播放| 国产精品av久久久久免费| 国产一区二区 视频在线| 啦啦啦啦在线视频资源| a级片在线免费高清观看视频| 汤姆久久久久久久影院中文字幕| 黄色怎么调成土黄色| 无遮挡黄片免费观看| 欧美xxⅹ黑人| 少妇人妻 视频| 亚洲精品成人av观看孕妇| 18禁国产床啪视频网站| 日韩中文字幕视频在线看片| 一级a爱视频在线免费观看| 成人亚洲精品一区在线观看| 免费久久久久久久精品成人欧美视频| 水蜜桃什么品种好| 国产av码专区亚洲av| 18禁裸乳无遮挡动漫免费视频| 国产1区2区3区精品| 亚洲精品国产av蜜桃| 在线天堂最新版资源| svipshipincom国产片| 精品少妇一区二区三区视频日本电影 | 欧美xxⅹ黑人| 久久久久久久久免费视频了| 男人添女人高潮全过程视频| 亚洲天堂av无毛| 亚洲四区av| 成人午夜精彩视频在线观看| 亚洲熟女精品中文字幕| 久久久久精品性色| 国产欧美亚洲国产| 超色免费av| 亚洲综合精品二区| 亚洲欧美成人精品一区二区| 国产亚洲欧美精品永久| 国产精品一二三区在线看| 欧美黑人欧美精品刺激| √禁漫天堂资源中文www| 另类精品久久| 青春草亚洲视频在线观看| 在线免费观看不下载黄p国产| 女性生殖器流出的白浆| 国产成人免费观看mmmm| 中文字幕精品免费在线观看视频| 国产精品久久久人人做人人爽| 只有这里有精品99| 五月天丁香电影| 亚洲国产精品999| 麻豆av在线久日| 午夜福利影视在线免费观看| kizo精华| 国产精品亚洲av一区麻豆 | 满18在线观看网站| 精品第一国产精品| 久久狼人影院| 午夜老司机福利片| 国产片特级美女逼逼视频| 国产一区二区激情短视频 | 90打野战视频偷拍视频| 国产精品人妻久久久影院| 国产欧美日韩一区二区三区在线| 精品免费久久久久久久清纯 | 一区二区av电影网| 国产老妇伦熟女老妇高清| 欧美老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 侵犯人妻中文字幕一二三四区| 精品国产一区二区三区久久久樱花| 亚洲综合色网址| 国产日韩欧美视频二区| 精品国产一区二区三区久久久樱花| 国产男女内射视频| 亚洲精品在线美女| av网站免费在线观看视频| 欧美精品av麻豆av| 国产一区二区 视频在线| av在线老鸭窝| 制服诱惑二区| 青青草视频在线视频观看| 乱人伦中国视频| 男人操女人黄网站| 99久国产av精品国产电影| 美女午夜性视频免费| 人人妻人人澡人人看| 亚洲欧美清纯卡通| 亚洲国产精品999| 国产精品欧美亚洲77777| 在线观看www视频免费| 亚洲精品aⅴ在线观看| 我的亚洲天堂| 午夜免费男女啪啪视频观看| 国产精品秋霞免费鲁丝片| 国产成人欧美| 男的添女的下面高潮视频| 久久久精品94久久精品| 国产免费又黄又爽又色| 国产极品粉嫩免费观看在线| 综合色丁香网| 亚洲少妇的诱惑av| 在线看a的网站| 两个人看的免费小视频| 亚洲少妇的诱惑av| 老鸭窝网址在线观看| 久久99热这里只频精品6学生| 亚洲视频免费观看视频| 中国三级夫妇交换| 黄色 视频免费看| 国产精品av久久久久免费| 国产97色在线日韩免费| 老熟女久久久| 男女床上黄色一级片免费看| 精品第一国产精品| 日韩视频在线欧美| 99精品久久久久人妻精品| 久久久久久免费高清国产稀缺| 在线天堂最新版资源| 国产免费一区二区三区四区乱码| 在线观看免费视频网站a站| 久久这里只有精品19| 国产精品国产三级国产专区5o| 欧美精品一区二区大全| 亚洲成人手机| 尾随美女入室| 成人18禁高潮啪啪吃奶动态图| 极品少妇高潮喷水抽搐| 十分钟在线观看高清视频www| 天堂俺去俺来也www色官网| 免费观看人在逋| 免费人妻精品一区二区三区视频| 精品国产国语对白av| 亚洲成人国产一区在线观看 | 多毛熟女@视频| 曰老女人黄片| 国产av码专区亚洲av| 日韩一卡2卡3卡4卡2021年| 久久久久久人妻| 亚洲成色77777| 看免费av毛片| 狠狠精品人妻久久久久久综合| 成人手机av| 女人爽到高潮嗷嗷叫在线视频| 久久女婷五月综合色啪小说| 制服人妻中文乱码| 在线免费观看不下载黄p国产| 亚洲精品aⅴ在线观看| 国产99久久九九免费精品| 精品一品国产午夜福利视频| 亚洲伊人久久精品综合| 少妇人妻精品综合一区二区| 婷婷色麻豆天堂久久| 国产精品一区二区在线不卡| 国产又色又爽无遮挡免| 又大又爽又粗| 欧美日韩国产mv在线观看视频| 999久久久国产精品视频| 18在线观看网站| 国产免费一区二区三区四区乱码| 女性生殖器流出的白浆| 成人漫画全彩无遮挡| 一区二区日韩欧美中文字幕| 色婷婷av一区二区三区视频| 在线亚洲精品国产二区图片欧美| 日韩av在线免费看完整版不卡| 免费在线观看视频国产中文字幕亚洲 | 韩国av在线不卡| 成年人午夜在线观看视频| 欧美av亚洲av综合av国产av | 午夜福利一区二区在线看| 赤兔流量卡办理| 亚洲av中文av极速乱| 操美女的视频在线观看| 欧美精品av麻豆av| 看免费av毛片| 日本欧美视频一区| videosex国产| 一区在线观看完整版| svipshipincom国产片| 国产1区2区3区精品| 国产片特级美女逼逼视频| 国产乱人偷精品视频| 亚洲七黄色美女视频| 久久热在线av| 亚洲欧美成人精品一区二区| 亚洲精品国产一区二区精华液| 国产深夜福利视频在线观看| 亚洲熟女精品中文字幕| 如何舔出高潮| 久久免费观看电影| 免费观看性生交大片5| 久久久国产精品麻豆| 久久ye,这里只有精品| 精品一品国产午夜福利视频| 免费不卡黄色视频| 青春草国产在线视频| 亚洲美女黄色视频免费看| 9热在线视频观看99| 亚洲av国产av综合av卡| 欧美日韩av久久| a级毛片在线看网站| 热99国产精品久久久久久7| 考比视频在线观看| 亚洲成国产人片在线观看| 国产乱来视频区| 国产精品国产三级专区第一集| av女优亚洲男人天堂| 久久精品国产亚洲av涩爱| 熟女少妇亚洲综合色aaa.| 丝袜喷水一区| 国产精品秋霞免费鲁丝片| 一边摸一边做爽爽视频免费| 男女免费视频国产| 久久久久国产精品人妻一区二区| 国产免费福利视频在线观看| 欧美精品亚洲一区二区| 交换朋友夫妻互换小说| 午夜福利影视在线免费观看| 制服丝袜香蕉在线| 无限看片的www在线观看| 亚洲成色77777| 黑人欧美特级aaaaaa片| 亚洲中文av在线| 国产一区二区 视频在线| 一本一本久久a久久精品综合妖精| 极品人妻少妇av视频| 精品国产一区二区三区四区第35| 日日撸夜夜添| 最近的中文字幕免费完整| 中文天堂在线官网| 国产精品久久久人人做人人爽| 精品一区二区免费观看| 在线亚洲精品国产二区图片欧美| 少妇的丰满在线观看| 熟妇人妻不卡中文字幕| 精品久久蜜臀av无| 男女午夜视频在线观看| 可以免费在线观看a视频的电影网站 | 飞空精品影院首页| 男女边摸边吃奶| 亚洲成人一二三区av| 色精品久久人妻99蜜桃| 美女视频免费永久观看网站| 亚洲人成电影观看| 免费黄网站久久成人精品| 婷婷色综合大香蕉| 日本91视频免费播放| 视频区图区小说| 黑人巨大精品欧美一区二区蜜桃| 国产精品一国产av| 99精品久久久久人妻精品| 亚洲五月色婷婷综合| 久久久久久久久久久久大奶| 国产精品蜜桃在线观看| 丝袜人妻中文字幕| 天堂俺去俺来也www色官网| 免费在线观看黄色视频的| 2021少妇久久久久久久久久久| av网站免费在线观看视频| 国产精品久久久久久精品电影小说| av在线观看视频网站免费| 日韩一区二区三区影片| 亚洲第一区二区三区不卡| 精品国产国语对白av| 纵有疾风起免费观看全集完整版| 狂野欧美激情性bbbbbb| 少妇被粗大猛烈的视频| 国产爽快片一区二区三区| 成年女人毛片免费观看观看9 | 青春草国产在线视频| 美女脱内裤让男人舔精品视频| 99久久99久久久精品蜜桃| 女性生殖器流出的白浆| 久久99热这里只频精品6学生| 看免费成人av毛片| 亚洲精品久久午夜乱码| 五月开心婷婷网| av卡一久久| 美女脱内裤让男人舔精品视频| 欧美日韩视频高清一区二区三区二| 尾随美女入室| xxxhd国产人妻xxx| 午夜福利在线免费观看网站| 国产熟女欧美一区二区| 王馨瑶露胸无遮挡在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产精品偷伦视频观看了| 日韩中文字幕欧美一区二区 | 亚洲成人国产一区在线观看 | 街头女战士在线观看网站| 亚洲av成人精品一二三区| 纯流量卡能插随身wifi吗| 视频在线观看一区二区三区| 免费黄色在线免费观看| 91国产中文字幕| 久久精品久久久久久久性| 成人亚洲欧美一区二区av| 久久韩国三级中文字幕| 亚洲av欧美aⅴ国产| 秋霞在线观看毛片| 亚洲国产av新网站| 午夜精品国产一区二区电影| 国产亚洲欧美精品永久| 亚洲欧美清纯卡通| 久久99精品国语久久久| 成人影院久久| 99精国产麻豆久久婷婷| 亚洲成人手机| 亚洲欧美一区二区三区国产| 考比视频在线观看| 中文字幕最新亚洲高清| 午夜福利,免费看| 日日撸夜夜添| 成年av动漫网址| 一区在线观看完整版| 黄色视频在线播放观看不卡| 国产精品 国内视频| 国产深夜福利视频在线观看| 久久精品熟女亚洲av麻豆精品| 国产熟女欧美一区二区| 日韩伦理黄色片| 亚洲欧美日韩另类电影网站| 久久av网站| 久久免费观看电影| 亚洲av电影在线观看一区二区三区| 中文字幕亚洲精品专区| 人成视频在线观看免费观看| 男男h啪啪无遮挡| 亚洲精品美女久久av网站| 欧美日韩成人在线一区二区| 久久青草综合色| 日韩 欧美 亚洲 中文字幕| 99久久精品国产亚洲精品| 如何舔出高潮| 午夜日本视频在线| av女优亚洲男人天堂| videosex国产| 99久久精品国产亚洲精品| 欧美 日韩 精品 国产| 国产精品久久久久久人妻精品电影 | 久久久国产精品麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 最近的中文字幕免费完整| 亚洲在久久综合| 啦啦啦在线观看免费高清www| xxxhd国产人妻xxx| 伦理电影大哥的女人| 国产1区2区3区精品| 黄频高清免费视频| 精品国产国语对白av| 国产毛片在线视频| 国产精品 欧美亚洲| 免费黄色在线免费观看| 纵有疾风起免费观看全集完整版| 一本久久精品| 看十八女毛片水多多多| 精品卡一卡二卡四卡免费| 99国产精品免费福利视频| 日本vs欧美在线观看视频| 美女国产高潮福利片在线看| 在线精品无人区一区二区三| 久久久久人妻精品一区果冻| 青春草视频在线免费观看| 黄片无遮挡物在线观看| 毛片一级片免费看久久久久| 熟女av电影| 大片电影免费在线观看免费| avwww免费| 国产av精品麻豆| av在线播放精品| avwww免费| 自拍欧美九色日韩亚洲蝌蚪91| 自线自在国产av| 午夜av观看不卡| 热re99久久国产66热| 99re6热这里在线精品视频| 操美女的视频在线观看| 一区二区三区激情视频| 欧美精品一区二区大全| 欧美另类一区| 久久国产精品男人的天堂亚洲|