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

    象山港浮游動(dòng)物β多樣性及其成分變化的環(huán)境因子解釋

    2017-10-13 12:29:13戴美霞朱藝峰毛碩乾
    生態(tài)學(xué)報(bào) 2017年17期
    關(guān)鍵詞:象山港嵌套周轉(zhuǎn)

    戴美霞, 朱藝峰, 2, 3, *,林 霞, 3,毛碩乾,2

    1 寧波大學(xué)海洋學(xué)院, 寧波 315211 2 寧波海洋研究院, 寧波 315832 3 寧波大學(xué)應(yīng)用海洋生物技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室, 寧波 315211

    象山港浮游動(dòng)物β多樣性及其成分變化的環(huán)境因子解釋

    戴美霞1, 朱藝峰1, 2, 3, *,林 霞1, 3,毛碩乾1,2

    1 寧波大學(xué)海洋學(xué)院, 寧波 315211 2 寧波海洋研究院, 寧波 315832 3 寧波大學(xué)應(yīng)用海洋生物技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室, 寧波 315211

    根據(jù)象山港24個(gè)站位浮游動(dòng)物樣品和配套環(huán)境數(shù)據(jù),采用R語(yǔ)言的gdm工具包對(duì)浮游動(dòng)物β多樣性進(jìn)行了廣義非相似性模型(GDM)分析,并利用betapart工具包對(duì)β多樣性進(jìn)行了成分(周轉(zhuǎn)和嵌套性)分解,探討了環(huán)境因子與浮游動(dòng)物β多樣性及其成分間的關(guān)系。GDM模型分析結(jié)果表明,有10個(gè)環(huán)境變量(表層水溫、溶解氧、水深、透明度、pH、葉綠素a、地理距離、電導(dǎo)率、鹽度、懸浮顆粒物)對(duì)浮游動(dòng)物β多樣性有影響,解釋了GDM模型偏差比例的75.2%。在這10個(gè)變量中,水溫、溶解氧、水深是驅(qū)動(dòng)β多樣性變化重要因子,3個(gè)變量累計(jì)相對(duì)貢獻(xiàn)比例占63.9%(以GDM模型偏差的可解釋比例作為100%),其他7個(gè)變量占36.1%;3個(gè)變量中,又以水溫最重要,相對(duì)貢獻(xiàn)比例占到38.4%。從各預(yù)測(cè)變量的影響梯度看,地理距離、pH和鹽度分別在約高于25 km、7.8和25,水溫、葉綠素a含量分別在約低于22 ℃和0.5 μg/L時(shí),隨著變量梯度的增加,β多樣性增大;而溶解氧、電導(dǎo)率、透明度和水深則隨著梯度增加β多樣性一直增加,懸浮顆粒物含量對(duì)β多樣性幾乎無(wú)影響。據(jù)β多樣性成分分解結(jié)果,象山港浮游動(dòng)物在時(shí)空變化上總體以周轉(zhuǎn)為主,嵌套性很低。在象山港,浮游動(dòng)物嵌套性主要發(fā)生于大型浮游動(dòng)物和幼體類群,尤其是幼體類群,嵌套性在時(shí)空上幾乎都高于周轉(zhuǎn)。進(jìn)一步的Pearson相關(guān)性分析表明,與大型浮游動(dòng)物嵌套性顯著相關(guān)的環(huán)境因子是水溫和溶解氧,而與幼體類嵌套性顯著相關(guān)的環(huán)境因子除水溫和溶解氧外,還有電導(dǎo)率和流速。

    象山港; 浮游動(dòng)物;β多樣性成分; 環(huán)境因子; 廣義非相似性模型(GDM)

    Abstract: To explore the relationships between environmental factors and zooplanktonic beta diversity—and its component changes—zooplankton samples and corresponding environmental data were collected from 24 stations in Xiangshan Bay. Zooplanktonic beta diversity was analyzed by a generalized dissimilarity model (GDM) with the gdm package in the R language, and was further separated into two components—turnover and nestedness—using the betapart package. Analytical results with GDM showed that 10 environmental variables (surface water temperature, dissolved oxygen, water depth, transparency, pH, chlorophyll a, geographical distance, conductivity, salinity, and suspended particulate matter) affected zooplanktonic beta diversity and could explain 75.2% of the GDM deviation. Among the 10 variables, water temperature, dissolved oxygen, and water depth were the important driving factors for beta diversity changes, which explained 63.9% of the cumulative relative contribution (taking the interpretable proportion of the GDM deviation as 100%), while the other 7 factors accounted for 36.1%. Furthermore, water temperature was the most important driving factor among these three factors, which accounted for 38.4% of the cumulative relative contribution. According to the gradient effects of predictors, beta diversity increased with the increasing gradients of factors when the geographical distance, pH, salinity, water temperature, and chlorophyll a were above ca. 25 km, 7.8, 25, and under 22°C and 0.5 μg/L respectively. Moreover, zooplanktonic beta diversity always increased with increasing dissolved oxygen, conductivity, transparency, and water depth. However, no obvious effect was found with the content change of suspended particulate matter. According to the results of partitioning beta diversity into turnover and nestedness components, turnover played a dominant role in the spatiotemporal changes of zooplankton, while the nestedness had small effects in Xiangshan Bay. Furthermore, the nestedness mainly occurred in the groups of macro-zooplankton and lava, especially in the latter group, and the spatiotemporal nestedness was almost always higher than turnover in the two groups. Further analyses of Pearson′s correlation demonstrated that the water temperature and dissolved oxygen were significantly correlated with the nestedness of macro-zooplankton, and that the water temperature, dissolved oxygen, and two other factors—conductivity and velocity—were significantly associated with the nestedness in the lava group.

    KeyWords: Xiangshan Bay; zooplankton; component of beta diversity; environmental factor; Generalised Dissimilarity Modelling

    生物多樣性度量是生態(tài)學(xué)重要問題之一,β多樣性是除α多樣性外的一個(gè)關(guān)鍵指標(biāo),它可定義為不同時(shí)空群落間的物種成分變化。自從Whittaker提出β多樣性概念后,β多樣性已長(zhǎng)期應(yīng)用于陸地、水生動(dòng)植物生態(tài)學(xué)研究和生物多樣性保護(hù)規(guī)劃中[1- 2],因此,有關(guān)β多樣性度量方法已不斷被關(guān)注和拓展[3- 4]。

    在眾多的β多樣性測(cè)量方法中,以平均非相似性的單指數(shù)測(cè)量法[5]、多元排序的空間梯度長(zhǎng)度法和方差分配法最常見[6]。單指數(shù)測(cè)量法很難定量分辨出不同環(huán)境梯度和地理距離分別對(duì)β多樣性的影響;梯度長(zhǎng)度法則主要基于線性模型,也難以反映生態(tài)學(xué)中最常見的2種非線性關(guān)系,即物種成分非相似性隨環(huán)境因子和空間距離增加呈曲線關(guān)系[7- 8],以及物種周轉(zhuǎn)速率與環(huán)境因子梯度間呈曲線關(guān)系[9];方差分配法盡管可以定量探索不同預(yù)測(cè)變量對(duì)β多樣性的影響,但只能對(duì)少于4個(gè)變量(組)進(jìn)行分析,因此,許多影響因子需要合并為變量組[10]。為解決這些問題,廣義非相似性模型(GDM,Generalised Dissimilarity Modelling)已得到開發(fā)并在不同領(lǐng)域得到大量應(yīng)用[11],其實(shí)質(zhì)是矩陣回歸,把各個(gè)預(yù)測(cè)變量通過(guò)樣條函數(shù)轉(zhuǎn)換為統(tǒng)一單位的生物成分非相似性,于是,可用不同預(yù)測(cè)變量、包括經(jīng)緯度變量來(lái)預(yù)測(cè)區(qū)域β多樣性變化,因此,GDM結(jié)果可全面代表研究區(qū)域的β多樣性格局[12]。

    β多樣性格局生成的實(shí)質(zhì)來(lái)自2個(gè)不同生態(tài)現(xiàn)象—地點(diǎn)間的物種替換和豐富度差異,即β多樣性可進(jìn)一步分為周轉(zhuǎn)和嵌套性2個(gè)成分[13]。所謂周轉(zhuǎn),指物種被替換(站位間有不同物種),而嵌套性則必須具備物種豐富度差異的子集特性[14],即群落間不僅存在豐富度差異,而且1個(gè)群落必須是另1個(gè)群落的子集(站位間無(wú)新物種),或者說(shuō),1個(gè)群落是另1個(gè)高豐富度群落非隨機(jī)物種丟失的結(jié)果,這種物種丟失的非隨機(jī)過(guò)程可由任何因素產(chǎn)生[15]。然而,廣泛應(yīng)用的Whittakerβ多樣性無(wú)法分辨這2個(gè)成分[16],而區(qū)分這2個(gè)成分可為分析群落形成機(jī)制提供一個(gè)更強(qiáng)大的工具。如2個(gè)研究區(qū)域有相等的β多樣性值,似乎β多樣性的格局是相似的,但β多樣性的形成機(jī)制可完全不同,因1個(gè)區(qū)域可來(lái)自是高周轉(zhuǎn)、低嵌套性,而另1個(gè)區(qū)域則可以是低周轉(zhuǎn)、高嵌套性的結(jié)果[17]。此外,區(qū)分這2個(gè)成分對(duì)了解生物地理學(xué)、生態(tài)學(xué)和規(guī)劃保護(hù)問題至關(guān)重要,因嵌套性地點(diǎn)間存在豐富度子集特性,因此,只保護(hù)豐富度最高地點(diǎn)即可,而周轉(zhuǎn)為主的區(qū)域則需要保護(hù)更多的地點(diǎn)[15]。

    在這2個(gè)β多樣性成分中,對(duì)嵌套性的認(rèn)識(shí)實(shí)際上比β多樣性概念還要早[15]。嵌套性是二進(jìn)制生態(tài)數(shù)據(jù)矩陣的屬性,可用于探索物種對(duì)生境片段化、對(duì)特定環(huán)境棲息地的適應(yīng)性等[18]。其定量方法研究較多,早期主要采用矩陣溫度及其改進(jìn)方法,隨后又提出并建議使用NODF(基于重疊和下降填充的嵌套性)方法[19]。然而,Baselga[15]近年用只含嵌套性、無(wú)物種周轉(zhuǎn)的β多樣性案例,指出了NODF方法不能準(zhǔn)確反映Whittakerβ多樣性變化,并提出了新的嵌套性計(jì)算方法。由于β多樣性成分分解新方法剛出現(xiàn)不久,目前知道這2者成分的相對(duì)重要性甚少[2]。

    訖今為止,國(guó)內(nèi)β多樣性研究大多針對(duì)陸地植物,對(duì)海洋生物定量研究很少。已有綜述指出[20],國(guó)內(nèi)針對(duì)浮游動(dòng)物幾乎只有α多樣性研究,探討浮游動(dòng)物β多樣性與環(huán)境因子的定量關(guān)系目前也缺乏分析方法,因而是一個(gè)薄弱環(huán)節(jié)。為此,本文利用先前已報(bào)道的不同網(wǎng)目采集的群落數(shù)據(jù)和同步環(huán)境監(jiān)測(cè)數(shù)據(jù)[21],進(jìn)一步挖掘象山港浮游動(dòng)物β多樣性及其成分變化,并探索環(huán)境因子的驅(qū)動(dòng)作用。目標(biāo)是探明象山港浮游動(dòng)物:(1)β多樣性的驅(qū)動(dòng)因子及其梯度影響;(2)β多樣性周轉(zhuǎn)與嵌套性成分的時(shí)空格局特征;(3)β多樣性成分與環(huán)境因子的相關(guān)性,以期為象山港海洋牧場(chǎng)建設(shè)和漁業(yè)管理提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 調(diào)查站位與采樣

    有關(guān)調(diào)查站位設(shè)置細(xì)節(jié)、室內(nèi)樣品處理、種類鑒定方法等見先前報(bào)道的群落比較研究[21]。采樣時(shí)間為2010年3月、7月、9月和12月。按象山港空間水文特征,把狹長(zhǎng)海灣劃分為6個(gè)區(qū)(I—VI),其中,I區(qū)為杭州灣沖淡水影響區(qū),并受區(qū)域圍墾影響;II區(qū)是與東海水交換的主要區(qū)域;III區(qū)為狹灣入口,受港口水交換影響,還受碼頭運(yùn)輸和大橋建設(shè)影響;IV區(qū)為狹灣區(qū),其潮汐混合強(qiáng)度較III區(qū)更弱,此外還受烏沙山電廠影響;V區(qū)主要受養(yǎng)殖、國(guó)華電廠、顏公河影響;VI區(qū)則受養(yǎng)殖、國(guó)華電廠、圍墾、鳧溪河影響。

    在各區(qū)劃分經(jīng)緯度網(wǎng)格(每格4.7 km × 5.6 km),各月分別在各區(qū)隨機(jī)取得一個(gè)經(jīng)緯度網(wǎng)格(去除圍墾區(qū)、養(yǎng)殖區(qū)等受限區(qū)域),并在網(wǎng)格內(nèi)任意取1個(gè)站位,以獲得更多含不同經(jīng)緯度的空間環(huán)境信息(圖1)。采樣時(shí),同時(shí)用淺水I、II型浮游生物網(wǎng)(網(wǎng)目分別為505 μm和160 μm)由底到表垂直拖網(wǎng)2次。全年共得24個(gè)站位樣品。

    圖1 采樣站位Fig.1 Sampling stations

    1.2 環(huán)境因子監(jiān)測(cè)

    現(xiàn)場(chǎng)用水質(zhì)儀(Horiba U-5000,日本)同步測(cè)定表層水0.5 m處的溫度(SWT)、鹽度(SLT)、溶解氧(DO)和pH;透明度(Trs)用圓盤法、水深(Dep)用繩長(zhǎng)法、流速(CV)用智能流速儀(LGY-II,南京)測(cè)定。水樣在室內(nèi)用電導(dǎo)率儀測(cè)定電導(dǎo)率(CDT),懸浮顆粒物(SPM)按重量法測(cè)定[22],葉綠素a (Chla,Chlorophyll a)熒光法測(cè)定[23]。

    1.3 β多樣性分析

    1.3.1 GDM分析

    GDM建模時(shí),以站位間所有物種Bray-Curtis非相似性作為響應(yīng)變量,以多項(xiàng)式樣條函數(shù)(polynomial spline function)變換的環(huán)境變量差值作為預(yù)測(cè)變量,再用指數(shù)函數(shù)連接響應(yīng)變量和預(yù)測(cè)變量,變換后得[11,24]:

    式中,i,j代表不同站位,dij為站位間群落成分非相似性;b為回歸常數(shù);k為預(yù)測(cè)變量(k=1, 2, …,p),xki、xkj分別為i,j站位的k變量監(jiān)測(cè)值。fk(xk)代表經(jīng)樣條擬合的轉(zhuǎn)換函數(shù),fk(xk)=a×I-spline,其中,a為常數(shù),I-spline是樣條函數(shù),該樣條函數(shù)形式為平滑單調(diào)增加、以符合隨站位間生態(tài)距離越遠(yuǎn)、群落成分差異性越大的生態(tài)特點(diǎn)[7]。為防止模型過(guò)度擬合,采用各個(gè)預(yù)測(cè)變量分布的0、50%和100%分位數(shù)作為樣條函數(shù)分段擬合節(jié)點(diǎn)(knot),并使用gdm工具包缺省劃分的3個(gè)樣條函數(shù)[11]。有關(guān)變量的樣條函數(shù)轉(zhuǎn)換細(xì)節(jié)參照文獻(xiàn)[25]。

    GDM擬合后,由于預(yù)測(cè)變量已通過(guò)樣條函數(shù)統(tǒng)一轉(zhuǎn)換為無(wú)量綱的生物成分非相似性單位,故不同預(yù)測(cè)變量對(duì)β多樣性的相對(duì)貢獻(xiàn)大小,可通過(guò)軟件輸出的偏響應(yīng)圖(Partial response graph)的曲線高低幅度直接進(jìn)行比較;而各預(yù)測(cè)變量在圖中的形狀變化可指示預(yù)測(cè)變量不同梯度對(duì)β多樣性影響[8,11]。GDM運(yùn)行后,同樣也會(huì)自動(dòng)去除對(duì)GDM無(wú)貢獻(xiàn)的預(yù)測(cè)變量,此時(shí),該變量的偏響應(yīng)圖不會(huì)輸出[26]。

    1.3.2 β多樣性的成分分解

    不難計(jì)算,當(dāng)2個(gè)站位有相同物種數(shù)時(shí),因豐富度相同,嵌套性不發(fā)生,即b=c時(shí),βsor=βsim,說(shuō)明β多樣性差異完全由空間周轉(zhuǎn)產(chǎn)生;當(dāng)2個(gè)站位有不相同物種數(shù)時(shí),因b≠c,則βsor≠βsim。很明顯,由于βsor和βsim在無(wú)嵌套性時(shí)是相等的,因此,它們的差值則代表了嵌套性(βnes)成分的測(cè)量,即:βnes=βsor-βsim。當(dāng)計(jì)算多站位情況時(shí),同樣可按照該公式推導(dǎo)計(jì)算,但需要考慮多站位間物種的共有信息,各成分詳細(xì)推導(dǎo)過(guò)程見文獻(xiàn)[15]。

    1.4 數(shù)據(jù)分析

    1.4.1 物種豐度數(shù)據(jù)處理

    I型、II型拖網(wǎng)樣品的2重復(fù)浮游動(dòng)物豐度分別取平均。對(duì)I型和II型網(wǎng)中的相同物種,取高豐度數(shù)據(jù)作為該物種有效采樣數(shù)據(jù)。如象山港主要優(yōu)勢(shì)種針刺擬哲水蚤(Paracalanusaculeatus)在I型網(wǎng)中的豐度僅為0.34個(gè)/m3,而在II型網(wǎng)高達(dá)1169個(gè)/m3[21],顯然,使用I型網(wǎng)或兩網(wǎng)平均值數(shù)據(jù),都不能代表該種在象山港中的豐度,而用II型網(wǎng)的高豐度數(shù)據(jù)更具代表性。

    1.4.2 統(tǒng)計(jì)分析

    GDM建模時(shí),使用浮游動(dòng)物的豐度數(shù)據(jù)。經(jīng)緯度采用UTM投影(WGS 1984 UTM Zone 51N),用以探討地理距離上的β多樣性差異,在GDM擬合時(shí),配套的經(jīng)緯度作為1個(gè)預(yù)測(cè)變量[11]。

    為探討β多樣性的周轉(zhuǎn)與嵌套性成分變化,物種豐度數(shù)據(jù)需經(jīng)P/A(Presence/Absence)轉(zhuǎn)換。由于象山港以橈足類和中小型浮游動(dòng)物為主[21],為探索不同類群可能存在的β多樣性成分上差異,還進(jìn)一步把浮游動(dòng)物分為橈足類、幼體類,以及大型浮游動(dòng)物類(毛顎類、水母類、端足類、漣蟲類、櫻蝦類、磷蝦類、糠蝦類、多毛類、仔魚類)3個(gè)類群分別進(jìn)行探討。β多樣性成分與環(huán)境因子間的Pearson相關(guān)系數(shù)顯著性經(jīng)Mantel permutation排列法檢驗(yàn),排列999次。以上所有統(tǒng)計(jì)分析經(jīng)R語(yǔ)言(版本2.15.3)編程,主要加載包為vegan、gdm[11]和betapart[28]。

    2 結(jié)果

    2.1 環(huán)境因子監(jiān)測(cè)結(jié)果

    不同環(huán)境因子的時(shí)空監(jiān)測(cè)結(jié)果列于表1。由表1可知,不同月份鹽度(SLT)監(jiān)測(cè)值變化不大。流速(CV)差異與采樣時(shí)的潮汐有關(guān),監(jiān)測(cè)點(diǎn)的水深(Dep)除與潮汐變化有關(guān)外,還與采樣點(diǎn)地形有關(guān)。懸浮顆粒物(SPM)9月明顯低于其他各月,故透明度(Trs)也以9月最大,且明顯高于其他月份。葉綠素a(Chla)以9月最高,12月最低;電導(dǎo)率(CDT)以12月最低,而pH低值主要發(fā)生在7月和9月。象山港不同季節(jié)水溫(SWT)變化明顯,以9月最高,3月最低;溶解氧(DO)變化趨勢(shì)與水溫相反。從空間上分析,港口到港底,鹽度、流速、pH基本呈下降態(tài)勢(shì),而葉綠素a、透明度、水溫增加。水深在港口區(qū)較淺,懸浮顆粒物主要分布在港口區(qū)與狹灣區(qū),電導(dǎo)率與溶解氧不同區(qū)域有所波動(dòng)。

    表1 環(huán)境因子的時(shí)空特征

    2.2 象山港不同浮游動(dòng)物類群時(shí)空差異

    由表2可知,橈足類豐富度(種類數(shù))4個(gè)月份總體相似,大型浮游動(dòng)物豐富度以7月和9月較高,幼體類也以7月和9月較多??臻g上,橈足類豐富度以港口區(qū)最高,港底區(qū)最低;大型浮游動(dòng)物在港口區(qū)和狹灣區(qū)相似,都高于港底區(qū);而幼體類空間分布較相似。

    2.3 浮游動(dòng)物β多樣性的GDM擬合及預(yù)測(cè)變量貢獻(xiàn)

    據(jù)GDM輸出結(jié)果,Null模型的總偏差量為29.40,而GDM預(yù)測(cè)變量模型產(chǎn)生的偏差量為7.29,故預(yù)測(cè)變量可解釋浮游動(dòng)物β多樣性比例為75.2% (即100%-7.29/29.4×100%),未能解釋比例占24.8%。

    在構(gòu)建GDM模型的11個(gè)變量中,未檢出流速對(duì)β多樣性的影響(表3)。從各變量相對(duì)貢獻(xiàn)百分比上看,列于前3位變量分別是表層水溫、溶解氧和水深,累計(jì)貢獻(xiàn)63.9%。其他變量累計(jì)占36.1%,其中,經(jīng)緯度產(chǎn)生的地理距離對(duì)β多樣性有貢獻(xiàn),比例占6.3%,幾乎與葉綠素a和pH的相對(duì)貢獻(xiàn)比例相當(dāng),也高于鹽度的相對(duì)貢獻(xiàn),指示在象山港β多樣性研究時(shí),需要考慮采樣點(diǎn)間空間距離產(chǎn)生的影響。

    表2 浮游動(dòng)物類群豐富度時(shí)空差異

    表3 預(yù)測(cè)變量對(duì)浮游動(dòng)物群落β多樣性的相對(duì)貢獻(xiàn)/%

    2.4 浮游動(dòng)物β多樣性的預(yù)測(cè)變量梯度影響

    預(yù)測(cè)變量對(duì)β多樣性的梯度影響見圖2。各變量的梯度影響大體可分為4種類型,第1種是變量?jī)H在高梯度區(qū)有影響,如來(lái)自經(jīng)緯度的地理距離(Geo)、pH和鹽度,分別在約高于25 km、7.8和25時(shí)會(huì)促進(jìn)β多樣性增加;第2種是僅在低梯度區(qū)有影響,如水溫和葉綠素a分別在約低于22 ℃和0.5 μg/L時(shí),β多樣性才隨梯度值的增加而增大;第3種是幾乎隨著變量值的增加一直增加,如溶解氧、電導(dǎo)率、透明度和水深;第4種是影響極低,如懸浮顆粒物含量幾乎不會(huì)對(duì)β多樣性產(chǎn)生影響。

    圖2 預(yù)測(cè)變量梯度對(duì)浮游動(dòng)物β多樣性影響的GDM偏響應(yīng)圖Fig.2 Partial response graph of GDM for the gradient effects of predictors on zooplanktonic β diversityGeo:地理距離Geographical distance/m; SLT:鹽度Salinity; SWT:表層水溫Surface water temperature/℃; Chla:葉綠素a Chlorophyll a/(μg/L); DO:溶解氧Dissolved oxygen/(mg/L); CDT:電導(dǎo)率Conductivity/(ms/cm); Trs:透明度Transparency/m; Dep:水深Depth/m; SPM:懸浮顆粒物Suspended particulate matter/(mg/L); x軸的變量單位為原始單位; y軸指示I-spline轉(zhuǎn)換函數(shù), 單位為群落非相似性連接單位-ln(1-dij)

    2.5 象山港浮游動(dòng)物β多樣性成分的時(shí)空特征

    象山港浮游動(dòng)物β多樣性及其周轉(zhuǎn)、嵌套性成分時(shí)空差異列于表4。對(duì)總浮游動(dòng)物,總β多樣性在7月份最低,主要原因是其周轉(zhuǎn)最低,指示該月份物種替換較少,即7月份群落成分比其他月份更穩(wěn)定;對(duì)比各月份的周轉(zhuǎn)與嵌套性結(jié)果,周轉(zhuǎn)均比嵌套性值高出2倍以上,顯示各月份不同地點(diǎn)間浮游動(dòng)物總體上都以物種替換為主。空間上,從港口區(qū)到港底區(qū)總β多樣性呈下降趨勢(shì),但周轉(zhuǎn)卻增加,顯然,嵌套性值從港口到港底區(qū)下降是導(dǎo)致空間上總β多樣性下降的主因。

    表4 浮游動(dòng)物β多樣性及其成分時(shí)空差異

    據(jù)表4還可知,3個(gè)類群浮游動(dòng)物的β多樣性及成分時(shí)空變化也明顯不同。首先,對(duì)橈足類,不同時(shí)空上的總β多樣性主要由周轉(zhuǎn)引起,嵌套性較小。在7月份,β多樣性較低,而在其他月份更相似,主要原因是橈足類周轉(zhuǎn)較低;從港口區(qū)到港底區(qū)β多樣性只略有下降,說(shuō)明橈足類在象山港周轉(zhuǎn)和嵌套性都較穩(wěn)定。其次,大型浮游動(dòng)物在7月份總β多樣性也最低,其他月份間相近,但7月份β多樣性中的周轉(zhuǎn)與嵌套性值比較接近,反映站位間的物種嵌套性主要發(fā)生在7月;空間上,港口區(qū)到港底區(qū)大型浮游動(dòng)物總β多樣性增加,這主要由物種周轉(zhuǎn)增加、以及嵌套性值僅在港口區(qū)和狹灣區(qū)較高引起。最后,對(duì)幼體類,除7月份總β多樣性最低、周轉(zhuǎn)與嵌套性相似外,其他月份的嵌套性都高于周轉(zhuǎn),指示其他月份幼體類β多樣性中的嵌套性比周轉(zhuǎn)起更重要作用;空間上,港口區(qū)幼體類的β多樣性比港底區(qū)高,主要原因是周轉(zhuǎn)和嵌套性在港口區(qū)都高于港底區(qū),而且嵌套性都高于周轉(zhuǎn)。

    2.6 不同類群浮游動(dòng)物β多樣性成分與環(huán)境因子相關(guān)性

    不同類群浮游動(dòng)物的β多樣性成分與環(huán)境因子間的相關(guān)性、及其Mantel檢驗(yàn)見表5。對(duì)橈足類,經(jīng)緯度、水深、水溫、溶解氧與周轉(zhuǎn)存在顯著相關(guān)性,但未檢出各因子對(duì)嵌套性有顯著相關(guān)。對(duì)大型浮游動(dòng)物,經(jīng)緯度、鹽度與物種周轉(zhuǎn)有顯著相關(guān)性,水溫和溶解氧則與嵌套性顯著相關(guān)。在幼體類,僅pH值與周轉(zhuǎn)顯著相關(guān),而流速、水溫、電導(dǎo)率、溶解氧與嵌套性顯著相關(guān)。對(duì)總浮游動(dòng)物而言,影響浮游動(dòng)物周轉(zhuǎn)的顯著相關(guān)因子分別是經(jīng)緯度、水深、pH、水溫、電導(dǎo)率和溶解氧,而影響嵌套性的顯著相關(guān)因子是流速、水溫和溶解氧。

    3 討論

    3.1 驅(qū)動(dòng)象山港浮游動(dòng)物β多樣性變化的主要環(huán)境因子

    據(jù)本研究結(jié)果,GDM根據(jù)樣條函數(shù)的分段曲線擬合功能,不僅可分辨出不同環(huán)境因子、地理距離分別產(chǎn)生的β多樣性影響,而且還能識(shí)別各變量梯度的影響速率。經(jīng)GDM擬合,10個(gè)預(yù)測(cè)變量解釋了浮游動(dòng)物β多樣性的75.2%,該解釋比例高于大多數(shù)水生生物、陸地生物β多樣性研究結(jié)果[24,29]。這說(shuō)明監(jiān)測(cè)的預(yù)測(cè)變量可較好地反映象山港浮游動(dòng)物β多樣性變化,也進(jìn)一步證實(shí):由于大多數(shù)海洋生物調(diào)查的物種數(shù)據(jù)偏少,往往很難建立代表性的物種模型,但GDM使用生物成分非相似性,故對(duì)數(shù)據(jù)點(diǎn)很少時(shí)特別有用,并可廣泛用于物種分布分析[30]。

    表5 浮游動(dòng)物類群β多樣性成分Pearson相關(guān)性與Mantel檢驗(yàn)

    Lon: 經(jīng)度; Lat: 緯度; 其他簡(jiǎn)寫變量見表1; *P<0.05;**P<0.01

    在監(jiān)測(cè)的變量中,水溫、溶解氧、水深是驅(qū)動(dòng)β多樣性變化3個(gè)最主要因素。除水深與浮游動(dòng)物空間容納量有關(guān)外,由于溶解氧含量往往伴隨季節(jié)水溫變化呈極強(qiáng)的負(fù)相關(guān)[31],說(shuō)明溶解氧對(duì)β多樣性的驅(qū)動(dòng)強(qiáng)烈受季節(jié)水溫影響,且當(dāng)水溫約在22℃以下時(shí),β多樣性隨水溫上升而快速增加,該水溫上升條件明顯發(fā)生于3月到7月份。據(jù)已有研究,象山港水溫從2月到8月為上升周期,隨水溫增加浮游動(dòng)物種類會(huì)快速增加[32]。從整個(gè)象山港浮游動(dòng)物優(yōu)勢(shì)種更替看,3月到7月種類更替率高達(dá)77.8%,而高溫期的7月到9月更替率極低(16.7%)[21]。這些結(jié)果指示,隨著水溫上升,β多樣性隨著浮游動(dòng)物更替率增加而增加,到高溫期,伴隨溶解氧的降低,β多樣性幾乎不再增加。

    除上述3個(gè)主要因素外,鹽度往往是驅(qū)動(dòng)浮游動(dòng)物空間分布的另一個(gè)重要因素[20],然而,本研究發(fā)現(xiàn)鹽度對(duì)β多樣性影響較低,主要原因是象山港鹽度梯度周年變化不大;再者,由于象山港水交換率低,外海高鹽水所攜帶的浮游動(dòng)物影響區(qū)域有限,因高鹽水充其量也只能到達(dá)象山港約1/3處[33]。本研究還發(fā)現(xiàn),往往作為站位描述的經(jīng)緯度(地理距離)卻對(duì)浮游動(dòng)物β多樣性有影響,其影響比例高于鹽度。當(dāng)站位間距離大約超過(guò)25 km時(shí),β多樣性影響開始增大,事實(shí)上,這反映了生物成分非相似性隨地理距離增加而增大特點(diǎn)[1, 7],這種非相似性差異顯然不是人類活動(dòng)產(chǎn)生,因此,在探討人類活動(dòng)對(duì)β多樣性影響的大尺度研究時(shí),地理距離尺度本身存在的β多樣性差異必須關(guān)注,否則所識(shí)別的驅(qū)動(dòng)因子也會(huì)不同。如對(duì)小型底棲動(dòng)物線蟲的4個(gè)尺度大小研究時(shí)發(fā)現(xiàn),在宏觀尺度上(>1500 km),物種周轉(zhuǎn)由食物資源的數(shù)量驅(qū)動(dòng),而在小尺度上(m),則由食物資源的質(zhì)量控制[4]。此外,由地理距離尺度產(chǎn)生的β多樣性差異甚至還會(huì)放大人類活動(dòng)的影響因子解釋,如在澳大利亞沿岸大型海藻β多樣性研究中,當(dāng)添加經(jīng)緯度變量時(shí),對(duì)GDM模型解釋量增加了不到1%[8];而在墨西哥,對(duì)節(jié)肢動(dòng)物分布的貢獻(xiàn)則以地理距離最重要[9],顯然,此時(shí)若不考慮地理距離影響,則該影響會(huì)計(jì)算到人類活動(dòng)的影響中,從而放大人類活動(dòng)的影響。

    3.2 不同類群浮游動(dòng)物周轉(zhuǎn)、嵌套性與環(huán)境因子關(guān)系

    本研究結(jié)果顯示,象山港浮游動(dòng)物β多樣性以物種周轉(zhuǎn)為主。對(duì)大型浮游動(dòng)物,經(jīng)緯度差異和鹽度是驅(qū)動(dòng)其周轉(zhuǎn)變化的顯著影響因子。在象山港,港口與港底區(qū)水環(huán)境很不同,早被認(rèn)為是2個(gè)不同的水體,各自以不同的浮游動(dòng)物優(yōu)勢(shì)種在港中部交匯混合[32],因此,經(jīng)緯度引起的地理距離相距越遠(yuǎn)物種差異越大、周轉(zhuǎn)越高。而鹽度越高周轉(zhuǎn)加快,這明顯與外海高鹽水?dāng)y帶大型浮游動(dòng)物、從而增加物種豐富度有關(guān)。如象山港肥胖箭蟲(Sagittaenflata)的出現(xiàn)可指示其來(lái)源于暖水性外海種[33],主要大型浮游動(dòng)物-百陶箭蟲(Sagittabedoti)也主要分布于鹽度梯度較大、外海高鹽水與沿岸水的交匯區(qū),雙生水母(Diphyeschamissonis)也同樣與表層鹽度有較強(qiáng)相關(guān)性[34]。對(duì)幼體類,物種周轉(zhuǎn)則顯著相關(guān)于pH值。從監(jiān)測(cè)結(jié)果看,象山港pH值以狹灣區(qū)和港底區(qū)最低,且發(fā)生在7月和9月份,主要原因是6—9月份為當(dāng)?shù)刎S水期,由于上游化工廠多,徑流會(huì)帶入酸化水,而酸化水會(huì)影響浮游動(dòng)物幼體發(fā)育,甚至死亡[35],從而降低浮游動(dòng)物的豐富度。而對(duì)橈足類,經(jīng)緯度、水深、水溫和溶解氧均與周轉(zhuǎn)有顯著相關(guān)性,由于象山港主要類群是橈足類,豐度幾乎占79%以上(不包括橈足類幼體)[21],因此,這些因素的解釋類似于總β多樣性,這已在3.1節(jié)中論述。

    從象山港浮游動(dòng)物的嵌套性結(jié)果看,除橈足類嵌套性極低外,大型浮游動(dòng)物、幼體類都存在不同程度的嵌套性。從概念上看,嵌套性指不同站位間的浮游動(dòng)物種類發(fā)生丟失,而不是指種類發(fā)生更替,因此,這明顯與種類的擴(kuò)散分布狀態(tài)有關(guān)。對(duì)大型浮游動(dòng)物而言,嵌套性主要發(fā)生在港口區(qū)和狹灣區(qū),驅(qū)動(dòng)因子主要是水溫和溶解氧。原因是港口區(qū)易受外海影響,水交換率高溶解氧較大,通過(guò)水交換也會(huì)帶來(lái)東海區(qū)大型浮游動(dòng)物,但水交換影響通常只侵入到象山港約1/3處[33],因此,從港口到港中部,隨著水動(dòng)力的減弱,主要大型浮游動(dòng)物的毛顎類、水母類、糠蝦類種類數(shù)分布幾乎呈遞減趨勢(shì)[21]。而這些優(yōu)勢(shì)類群在象山港海域主要發(fā)生在7—9月的高溫期,此時(shí),隨著臺(tái)灣暖流靠近長(zhǎng)江口,帶來(lái)了毛顎類[32]、水母類[34]、糠蝦類、磷蝦類、十足類,以及仔魚等象山港優(yōu)勢(shì)類群[21]。在幼體類中,嵌套性在時(shí)空上大多高于周轉(zhuǎn),空間上也以港口區(qū)嵌套性更高,驅(qū)動(dòng)幼體類嵌套性的主要因子除與大型浮游動(dòng)物相似的水溫、溶解氧影響外,還有流速和電導(dǎo)率。眾所周知,幼體類的運(yùn)動(dòng)能力比大型浮游動(dòng)物弱的多,流速更易改變幼體類的擴(kuò)散分布,這是驅(qū)動(dòng)幼體類嵌套性顯著變化的一個(gè)影響因素,其中,港口區(qū)流速受潮汐影響更大,而港底區(qū)水動(dòng)力則很弱,使得港口區(qū)幼體類分布受影響更大。而電導(dǎo)率驅(qū)動(dòng)浮游動(dòng)物嵌套性變化目前不清楚,據(jù)報(bào)道,電導(dǎo)率與水體中可溶性營(yíng)養(yǎng)物有關(guān)[36],因此,有關(guān)電導(dǎo)率與幼體嵌套性關(guān)系是否與水體中營(yíng)養(yǎng)物分布有關(guān)有待深入研究。

    [1] 陳圣賓, 歐陽(yáng)志云, 徐衛(wèi)華, 肖燚. Beta多樣性研究進(jìn)展. 生物多樣性, 2010, 18(4): 323- 335.

    [3] Legendre P, De Cáceres M. Beta diversity as the variance of community data: dissimilarity coefficients and partitioning. Ecology Letters, 2013, 16(8): 951- 963.

    [4] Gambi C, Pusceddu A, Benedetti-Cecchi L, Danovaro R. Species richness, species turnover and functional diversity in nematodes of the deep Mediterranean Sea: searching for drivers at different spatial scales. Global Ecology and Biogeography, 2014, 23(1): 24- 39.

    [5] Koleff P, Gaston K J, Lennon J J. Measuring beta diversity for presence-absence data. Journal of Animal Ecology, 2003, 72(3): 367- 382.

    [6] Borcard D, Legendre P, Drapeau P. Partialling out the spatial component of ecological variation. Ecology, 1992, 73(3): 1045- 1055.

    [7] Nekola J C, White P S. The distance decay of similarity in biogeography and ecology. Journal of Biogeography, 1999, 26(4): 867- 878.

    [8] Leaper R, Hill N A, Edgar G J, Ellis N, Lawrence E, Pitcher C R, Barrett N S, Thomson R. Predictions of beta diversity for reef macroalgae across southeastern Australia. Ecosphere, 2011, 2(7): 1- 18.

    [9] Bell K L, Heard T A, Manion G, Ferrier S, van Klinken R D. The role of geography and environment in species turnover: phytophagous arthropods on a Neotropical legume. Journal of Biogeography, 2013, 40(9): 1755- 1766.

    [10] Borcard D, Gillet F, Legendre P.數(shù)量生態(tài)學(xué)——R語(yǔ)言的應(yīng)用. 賴江山, 譯. 北京: 高等教育出版社, 2014: 137- 201.

    [11] Ferrier S, Manion G, Elith J, Richardson K J. Using generalized dissimilarity modelling to analyse and predict patterns of beta diversity in regional biodiversity assessment. Diversity and Distributions, 2007, 13(3): 252- 264.

    [12] Rosauer D F, Ferrier S, Manion G, Laffan S, Williams K. Nice weather for frogs-using environmental data to model phylogenetic turnover // Lees B G, Laffan S W, eds. 10th International Conference on GeoComputation. Australia: University of New South Wales, 2009.

    [13] Baselga A. The relationship between species replacement, dissimilarity derived from nestedness, and nestedness. Global Ecology and Biogeography, 2012, 21(12): 1223- 1232.

    [14] Ulrich W, Gotelli N J. Null model analysis of species nestedness patterns. Ecology, 2007, 88(7): 1824- 1831.

    [15] Baselga A. Partitioning the turnover and nestedness components of beta diversity. Global Ecology and Biogeography, 2010, 19(1): 134- 143.

    [16] Harrison S, Ross S J, Lawton J H. Beta diversity on geographic gradients in Britain. Journal of Animal Ecology, 1992, 61(1): 151- 158.

    [17] Baselga A, Gómez-Rodríguez C, Lobo J M. Historical legacies in world amphibian diversity revealed by the turnover and nestedness components of beta diversity. PLoS One, 2012, 7(2): e32341.

    [18] Gururaja K V S, Ramachandra T V. Nestedness pattern in freshwater fishes of the Western Ghats: an indication of stream islands along riverscapes. Current Science, 2008, 95(12): 1707- 1714.

    [19] Almeida-Neto M, Guimar?es P, Guimar?es P R, Loyola R D, Ulrich W. A consistent metric for nestedness analysis in ecological systems: reconciling concept and measurement. Oikos, 2008, 117(8): 1227- 1239.

    [20] 徐兆禮. 中國(guó)近海浮游動(dòng)物多樣性研究的過(guò)去和未來(lái). 生物多樣性, 2011, 19(6): 635- 645.

    [21] 朱藝峰, 王銀, 林霞, 吳燕萍, 趙金龍, 嚴(yán)小軍. 象山港兩種網(wǎng)目網(wǎng)采浮游動(dòng)物群落比較. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(8): 2277- 2286.

    [22] 中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局, 中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì). GB/T 17378.4—2007 海洋監(jiān)測(cè)規(guī)范 第4部分: 海水分析. 北京: 中國(guó)標(biāo)準(zhǔn)出版社, 2008.

    [23] 中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局, 中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì). GB/T 12763.6—2007 海洋調(diào)查規(guī)范 第6部分: 海洋生物調(diào)查. 北京: 中國(guó)標(biāo)準(zhǔn)出版社, 2008.

    [24] Overton J M, Barker G M, Price R. Estimating and conserving patterns of invertebrate diversity: a test case of New Zealand land snails. Diversity and Distributions, 2009, 15(5): 731- 741.

    [25] Manion G. A technique for constructing monotonic regression splines to enable non-linear transformation of GIS rasters // Anderssen R S, Braddock R D, Newham L T H, eds. 18th World IMACS congress and MODSIM09 international congress on modelling and simulation. Australia: The Modelling and Simulation Society of Australia and New Zealand Inc. and the International Association for Mathematics and Computers in Simulation, 2009: 13- 17.

    [26] Valdujo P H, Carnaval A C O Q, Graham C H. Environmental correlates of anuran beta diversity in the Brazilian Cerrado. Ecography, 2013, 36(6): 708- 717.

    [27] Lennon J J, Koleff P, Greenwood J J D, Gaston K J. The geographical structure of British bird distributions: diversity, spatial turnover and scale. Journal of Animal Ecology, 2001, 70(6): 966- 979.

    [28] Baselga A, Orme C D L. Betapart: an R package for the study of beta diversity. Methods in Ecology and Evolution, 2012, 3(5): 808- 812.

    [29] Koubbi P, Moteki M, Duhamel G, Goarant A, Hulley P A, O′Driscoll R, Ishimaru T, Pruvost P, Tavernier E, Hosie G. Ecoregionalization of myctophid fish in the Indian sector of the Southern Ocean: Results from generalized dissimilarity models. Deep Sea Research II: Topical Studies in Oceanography, 2011, 58(1/2): 170- 180.

    [30] Burgman M A, Lindenmayer D B, Elith J. Managing landscapes for conservation under uncertainty. Ecology, 2005, 86(8): 2007- 2017.

    [31] Wu M L, Wang Y S, Sun C C, Wang H L, Dong J D, Han S H. Identification of anthropogenic effects and seasonality on water quality in Daya Bay, South China Sea. Journal of Environmental Management, 2009, 90(10): 3082- 3090.

    [32] 柏懷萍, 陳亞瞿, 宋海棠, 邊元慈. 象山港浮游動(dòng)物調(diào)查報(bào)告. 海洋漁業(yè), 1984, (6): 249- 253.

    [33] 王曉波, 邱武生, 秦銘俐, 魏永杰. 象山港浮游動(dòng)物生態(tài)群落分布的研究. 海洋環(huán)境科學(xué), 2009, 28(Z1): 62- 64.

    [34] 陳小慶, 陳斌, 黃備, 王婕妤, 鄭基, 寧平, 俞存根. 夏季舟山漁場(chǎng)及鄰近海域浮游動(dòng)物群落結(jié)構(gòu)特征分析. 動(dòng)物學(xué)研究, 2010, 31(1): 99- 107.

    [35] 趙信國(guó), 劉廣緒. 海洋酸化對(duì)海洋無(wú)脊椎動(dòng)物的影響研究進(jìn)展. 生態(tài)學(xué)報(bào), 2015, 35(7): 2388- 2398.

    [36] Walker C E, Pan Y D. Using diatom assemblages to assess urban stream conditions. Hydrobiologia, 2006, 516(1): 179- 189.

    InterpretationofenvironmentalfactorsaffectingzooplanktonicbetadiversityanditscomponentsinXiangshanBay

    DAI Meixia1, ZHU Yifeng1, 2, 3,*, LIN Xia1,3, MAO Shuoqian1,2

    1SchoolofMarineSciences,NingboUniversity,Ningbo315211,China2NingboInstituteofOceanography,Ningbo315832,China3KeyLaboratoryofAppliedMarineBiotechnology,theMinistryofEducation,NingboUniversity,Ningbo315211,China

    海洋公益性行業(yè)科研專項(xiàng)(201105009- 3);國(guó)家科技支撐項(xiàng)目(2011BAD13B08);浙江省自然科學(xué)基金(Y5080274, Y5100369);寧波大學(xué)水產(chǎn)養(yǎng)殖浙江省重中之重一級(jí)學(xué)科開放基金資助(421500052)

    2016- 06- 06; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 04- 24

    10.5846/stxb201606061083

    *通訊作者Corresponding author.E-mail: zhuyifeng@nbu.edu.cn

    戴美霞, 朱藝峰,林霞,毛碩乾.象山港浮游動(dòng)物β多樣性及其成分變化的環(huán)境因子解釋.生態(tài)學(xué)報(bào),2017,37(17):5780- 5789.

    Dai M X, Zhu Y F, Lin X, Mao S Q.Interpretation of environmental factors affecting zooplanktonic beta diversity and its components in Xiangshan Bay.Acta Ecologica Sinica,2017,37(17):5780- 5789.

    猜你喜歡
    象山港嵌套周轉(zhuǎn)
    例析“立幾”與“解幾”的嵌套問題
    象山港周邊主要入海污染物特征研究
    基于嵌套Logit模型的競(jìng)爭(zhēng)性選址問題研究
    關(guān)于壓縮貨車周轉(zhuǎn)時(shí)間的探討
    基于SolidWorks周轉(zhuǎn)輪系裝配與運(yùn)動(dòng)仿真
    象山港海域生態(tài)分區(qū)研究
    象山港海灣生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估
    周轉(zhuǎn)性材料租賃參考價(jià)格
    C35復(fù)合雙摻高性能混凝土在象山港跨海大橋中的應(yīng)用
    一種基于區(qū)分服務(wù)的嵌套隊(duì)列調(diào)度算法
    avwww免费| 日韩视频一区二区在线观看| 色综合欧美亚洲国产小说| 久久久久久久国产电影| 国产精品久久久久久人妻精品电影 | 国产欧美日韩精品亚洲av| 乱人伦中国视频| 亚洲精品中文字幕一二三四区 | 亚洲全国av大片| 久久免费观看电影| 一区在线观看完整版| 窝窝影院91人妻| 99riav亚洲国产免费| 日韩欧美一区视频在线观看| 日韩一区二区三区影片| 欧美国产精品va在线观看不卡| 国产激情久久老熟女| 亚洲精华国产精华精| 午夜福利,免费看| 国产无遮挡羞羞视频在线观看| 日日夜夜操网爽| 18禁美女被吸乳视频| 国产欧美日韩精品亚洲av| 色在线成人网| 国产精品亚洲一级av第二区| 国产成人免费无遮挡视频| 久久久国产成人免费| 老司机影院毛片| 国产精品久久久人人做人人爽| 国产伦理片在线播放av一区| 色婷婷久久久亚洲欧美| www.精华液| 69精品国产乱码久久久| 欧美精品高潮呻吟av久久| 曰老女人黄片| 黄色视频在线播放观看不卡| 国产成人精品久久二区二区免费| 免费不卡黄色视频| 麻豆av在线久日| 欧美日本中文国产一区发布| 久久热在线av| 久久婷婷成人综合色麻豆| 在线看a的网站| 一区二区三区乱码不卡18| 19禁男女啪啪无遮挡网站| 99九九在线精品视频| 香蕉国产在线看| 国产三级黄色录像| 女人爽到高潮嗷嗷叫在线视频| 国产av精品麻豆| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产精品一区二区三区| 国产一区二区在线观看av| 亚洲精品美女久久av网站| 女人高潮潮喷娇喘18禁视频| 99久久99久久久精品蜜桃| 欧美成人免费av一区二区三区 | 久久中文字幕一级| 蜜桃在线观看..| 欧美日韩亚洲高清精品| 亚洲中文日韩欧美视频| 国产片内射在线| 亚洲国产av影院在线观看| 久久香蕉激情| 国产欧美日韩一区二区精品| 色94色欧美一区二区| 国产亚洲精品第一综合不卡| 国产精品98久久久久久宅男小说| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久精品电影小说| 精品国产乱子伦一区二区三区| 成人国语在线视频| 电影成人av| 十八禁网站网址无遮挡| 久久亚洲精品不卡| av天堂在线播放| 精品一区二区三卡| 国产精品影院久久| 亚洲成人免费电影在线观看| 不卡av一区二区三区| svipshipincom国产片| 午夜精品久久久久久毛片777| 日韩大码丰满熟妇| 亚洲色图 男人天堂 中文字幕| 日韩欧美免费精品| 十八禁网站免费在线| 亚洲欧美色中文字幕在线| 欧美日韩黄片免| 老司机亚洲免费影院| 国产成+人综合+亚洲专区| 热99国产精品久久久久久7| 国产av国产精品国产| 午夜成年电影在线免费观看| 日本av免费视频播放| 男女边摸边吃奶| 亚洲,欧美精品.| 国产精品麻豆人妻色哟哟久久| 久久精品熟女亚洲av麻豆精品| 女人被躁到高潮嗷嗷叫费观| 欧美精品啪啪一区二区三区| 午夜福利视频精品| 中文字幕最新亚洲高清| 久久久精品国产亚洲av高清涩受| 日韩中文字幕视频在线看片| 国产麻豆69| 最新的欧美精品一区二区| 精品人妻熟女毛片av久久网站| 在线观看免费视频网站a站| 一边摸一边做爽爽视频免费| 精品高清国产在线一区| 久久精品亚洲精品国产色婷小说| 在线 av 中文字幕| 超碰成人久久| 天堂8中文在线网| 婷婷丁香在线五月| 久久人妻福利社区极品人妻图片| 午夜福利影视在线免费观看| 手机成人av网站| 蜜桃国产av成人99| 国产精品香港三级国产av潘金莲| 香蕉丝袜av| 成人国产av品久久久| 久久中文字幕一级| 男女边摸边吃奶| 精品熟女少妇八av免费久了| 亚洲精品国产一区二区精华液| 黄色怎么调成土黄色| 国产成人免费无遮挡视频| 免费一级毛片在线播放高清视频 | 自线自在国产av| 日韩成人在线观看一区二区三区| a级片在线免费高清观看视频| 91麻豆精品激情在线观看国产 | 国产亚洲精品久久久久5区| 黄频高清免费视频| 国产成人精品无人区| 最新的欧美精品一区二区| 亚洲精品美女久久av网站| 国产精品自产拍在线观看55亚洲 | 一本一本久久a久久精品综合妖精| 蜜桃在线观看..| 99在线人妻在线中文字幕 | 下体分泌物呈黄色| 亚洲精品中文字幕在线视频| 99re6热这里在线精品视频| 别揉我奶头~嗯~啊~动态视频| 久久人人爽av亚洲精品天堂| 精品人妻1区二区| 大香蕉久久网| 国产欧美日韩综合在线一区二区| 久久精品国产综合久久久| 伊人久久大香线蕉亚洲五| 国产免费现黄频在线看| 老司机午夜十八禁免费视频| 搡老熟女国产l中国老女人| tube8黄色片| 久久亚洲真实| 日日夜夜操网爽| 精品人妻1区二区| 在线看a的网站| 免费av中文字幕在线| 一级毛片精品| 一边摸一边抽搐一进一出视频| 成年人午夜在线观看视频| 国产精品国产高清国产av | 高清av免费在线| 日本av手机在线免费观看| 亚洲av片天天在线观看| av不卡在线播放| 国产精品一区二区精品视频观看| 日本五十路高清| 一区二区日韩欧美中文字幕| 亚洲久久久国产精品| 国产精品自产拍在线观看55亚洲 | 丝袜人妻中文字幕| 大片电影免费在线观看免费| 男女高潮啪啪啪动态图| 大片免费播放器 马上看| 一边摸一边抽搐一进一出视频| 老司机深夜福利视频在线观看| 女同久久另类99精品国产91| 欧美 日韩 精品 国产| a级片在线免费高清观看视频| 国产三级黄色录像| 国产亚洲精品第一综合不卡| 亚洲午夜理论影院| 亚洲午夜精品一区,二区,三区| 老汉色∧v一级毛片| 久久精品aⅴ一区二区三区四区| 狠狠婷婷综合久久久久久88av| 亚洲五月婷婷丁香| 免费不卡黄色视频| 欧美黄色淫秽网站| 操美女的视频在线观看| 丁香六月欧美| 亚洲自偷自拍图片 自拍| 超碰97精品在线观看| 日本欧美视频一区| 大陆偷拍与自拍| 欧美日韩视频精品一区| 老司机午夜福利在线观看视频 | 国产精品久久久久成人av| 一本综合久久免费| 色婷婷久久久亚洲欧美| 黄色 视频免费看| 日韩欧美一区二区三区在线观看 | 丝袜在线中文字幕| 久久精品国产a三级三级三级| 日韩制服丝袜自拍偷拍| 午夜成年电影在线免费观看| 丁香六月欧美| 婷婷成人精品国产| 午夜福利,免费看| 亚洲国产精品一区二区三区在线| 俄罗斯特黄特色一大片| 亚洲av国产av综合av卡| 免费观看a级毛片全部| 成人三级做爰电影| 国产成+人综合+亚洲专区| 日韩视频在线欧美| 免费女性裸体啪啪无遮挡网站| 91av网站免费观看| 欧美在线一区亚洲| 国产aⅴ精品一区二区三区波| 免费在线观看视频国产中文字幕亚洲| 久久久久网色| 在线观看www视频免费| 黄色视频不卡| 嫁个100分男人电影在线观看| 欧美一级毛片孕妇| 成年版毛片免费区| 涩涩av久久男人的天堂| 国产真人三级小视频在线观看| 亚洲久久久国产精品| 成年人免费黄色播放视频| bbb黄色大片| 亚洲伊人久久精品综合| 搡老乐熟女国产| 国产精品1区2区在线观看. | 一区二区av电影网| 精品第一国产精品| 老熟妇乱子伦视频在线观看| 好男人电影高清在线观看| kizo精华| 欧美日韩成人在线一区二区| 午夜久久久在线观看| 丰满迷人的少妇在线观看| 久久精品亚洲熟妇少妇任你| 精品亚洲成a人片在线观看| 久久久欧美国产精品| 精品国产乱码久久久久久男人| 国内毛片毛片毛片毛片毛片| 中亚洲国语对白在线视频| 亚洲情色 制服丝袜| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩一区二区三| 少妇裸体淫交视频免费看高清 | a级片在线免费高清观看视频| 精品少妇黑人巨大在线播放| tocl精华| 国产精品一区二区免费欧美| 男男h啪啪无遮挡| 免费日韩欧美在线观看| 他把我摸到了高潮在线观看 | 亚洲综合色网址| 欧美激情久久久久久爽电影 | 香蕉久久夜色| 黑人巨大精品欧美一区二区mp4| 丝瓜视频免费看黄片| 国产又色又爽无遮挡免费看| 欧美日本中文国产一区发布| 一夜夜www| 国产极品粉嫩免费观看在线| √禁漫天堂资源中文www| 精品卡一卡二卡四卡免费| netflix在线观看网站| 日韩欧美国产一区二区入口| 深夜精品福利| 汤姆久久久久久久影院中文字幕| 91大片在线观看| 欧美性长视频在线观看| 亚洲伊人色综图| av在线播放免费不卡| 亚洲精品一卡2卡三卡4卡5卡| 黄色视频在线播放观看不卡| 欧美亚洲日本最大视频资源| 亚洲男人天堂网一区| 日本精品一区二区三区蜜桃| 在线观看一区二区三区激情| 精品人妻1区二区| 国产精品麻豆人妻色哟哟久久| 欧美精品高潮呻吟av久久| 国产欧美日韩一区二区三区在线| 在线观看www视频免费| 一边摸一边抽搐一进一出视频| 亚洲午夜精品一区,二区,三区| 少妇猛男粗大的猛烈进出视频| 亚洲va日本ⅴa欧美va伊人久久| 女同久久另类99精品国产91| 国产精品久久久人人做人人爽| 国产av精品麻豆| 午夜精品久久久久久毛片777| 最新美女视频免费是黄的| 黑人操中国人逼视频| 午夜福利影视在线免费观看| 久久人人爽av亚洲精品天堂| 亚洲午夜理论影院| 丁香六月欧美| 午夜日韩欧美国产| 久久国产精品人妻蜜桃| av不卡在线播放| 成人国产一区最新在线观看| 欧美老熟妇乱子伦牲交| 亚洲av日韩精品久久久久久密| 久久久久久久久久久久大奶| 老司机午夜十八禁免费视频| 国产精品熟女久久久久浪| 免费看a级黄色片| 天天添夜夜摸| 国产成人免费观看mmmm| 亚洲午夜理论影院| 久久人人97超碰香蕉20202| 午夜福利免费观看在线| 国产男女内射视频| 欧美日韩av久久| videos熟女内射| 在线观看免费日韩欧美大片| 欧美 亚洲 国产 日韩一| 1024香蕉在线观看| 免费久久久久久久精品成人欧美视频| 午夜精品国产一区二区电影| 久久久国产成人免费| 免费日韩欧美在线观看| 天堂动漫精品| 欧美一级毛片孕妇| 无遮挡黄片免费观看| 国产黄频视频在线观看| 九色亚洲精品在线播放| 日韩视频一区二区在线观看| 最近最新中文字幕大全电影3 | 国产成人啪精品午夜网站| 久久久久久久大尺度免费视频| 日韩大码丰满熟妇| 中文字幕精品免费在线观看视频| av有码第一页| 亚洲精品国产色婷婷电影| 妹子高潮喷水视频| 一个人免费在线观看的高清视频| 国内毛片毛片毛片毛片毛片| 亚洲成a人片在线一区二区| 亚洲精品国产色婷婷电影| 丝袜美足系列| 日本五十路高清| 精品久久久久久久毛片微露脸| 香蕉国产在线看| 日本av手机在线免费观看| 久久久久久久大尺度免费视频| 啦啦啦 在线观看视频| 亚洲天堂av无毛| 黄色a级毛片大全视频| 乱人伦中国视频| 亚洲av电影在线进入| 99精品久久久久人妻精品| 老司机福利观看| 亚洲熟女毛片儿| 久久久欧美国产精品| 汤姆久久久久久久影院中文字幕| 91精品国产国语对白视频| 亚洲精品国产区一区二| 精品乱码久久久久久99久播| 一级片免费观看大全| √禁漫天堂资源中文www| 女性被躁到高潮视频| 国产不卡一卡二| 亚洲av欧美aⅴ国产| 99香蕉大伊视频| 久久久久网色| 亚洲熟妇熟女久久| 国精品久久久久久国模美| 成年动漫av网址| 69精品国产乱码久久久| 一二三四在线观看免费中文在| a级片在线免费高清观看视频| www.熟女人妻精品国产| 巨乳人妻的诱惑在线观看| 狂野欧美激情性xxxx| 国产成人欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品成人av观看孕妇| 母亲3免费完整高清在线观看| 精品少妇内射三级| 90打野战视频偷拍视频| 脱女人内裤的视频| 国产欧美日韩精品亚洲av| 高清av免费在线| 欧美+亚洲+日韩+国产| bbb黄色大片| 亚洲三区欧美一区| 国产av又大| 欧美精品啪啪一区二区三区| 国产福利在线免费观看视频| netflix在线观看网站| 久久国产精品大桥未久av| 正在播放国产对白刺激| 久久中文看片网| 中文字幕色久视频| 少妇粗大呻吟视频| 人人妻人人澡人人爽人人夜夜| 在线观看免费午夜福利视频| av不卡在线播放| 男女无遮挡免费网站观看| 男女边摸边吃奶| 亚洲色图综合在线观看| 久久久久久免费高清国产稀缺| 啦啦啦免费观看视频1| 亚洲精品自拍成人| 国产人伦9x9x在线观看| 美女国产高潮福利片在线看| 欧美精品高潮呻吟av久久| 国产免费现黄频在线看| 美女主播在线视频| 窝窝影院91人妻| av有码第一页| 欧美一级毛片孕妇| 黄网站色视频无遮挡免费观看| 国产三级黄色录像| 日韩大片免费观看网站| 亚洲欧美一区二区三区久久| av线在线观看网站| 叶爱在线成人免费视频播放| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 少妇被粗大的猛进出69影院| 汤姆久久久久久久影院中文字幕| 成年人黄色毛片网站| 一进一出好大好爽视频| 欧美日韩成人在线一区二区| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产成人一精品久久久| 男女免费视频国产| 我的亚洲天堂| 丁香欧美五月| 国产熟女午夜一区二区三区| 国产av国产精品国产| 精品视频人人做人人爽| 亚洲人成电影免费在线| 丝袜美足系列| 热99久久久久精品小说推荐| 日韩一卡2卡3卡4卡2021年| 美国免费a级毛片| 久久av网站| 精品久久久久久电影网| 成人特级黄色片久久久久久久 | 午夜免费成人在线视频| 51午夜福利影视在线观看| 丰满人妻熟妇乱又伦精品不卡| tube8黄色片| 亚洲人成电影观看| 啦啦啦视频在线资源免费观看| 久久久欧美国产精品| 午夜老司机福利片| 精品少妇久久久久久888优播| 一边摸一边抽搐一进一出视频| 少妇被粗大的猛进出69影院| 中亚洲国语对白在线视频| 精品久久久精品久久久| 在线观看一区二区三区激情| 国产成人影院久久av| 99九九在线精品视频| 少妇精品久久久久久久| 久久亚洲真实| 色尼玛亚洲综合影院| 日本五十路高清| 亚洲精华国产精华精| 久久久久精品国产欧美久久久| 亚洲美女黄片视频| 欧美精品一区二区大全| 亚洲欧洲精品一区二区精品久久久| 18禁国产床啪视频网站| 黑人欧美特级aaaaaa片| 精品久久久久久久毛片微露脸| 久久青草综合色| 欧美黄色片欧美黄色片| 久久久国产一区二区| 精品亚洲成国产av| 免费av中文字幕在线| 久久狼人影院| 国产一区二区 视频在线| 蜜桃在线观看..| 青青草视频在线视频观看| 性高湖久久久久久久久免费观看| 久久人妻av系列| 欧美黑人精品巨大| 成年女人毛片免费观看观看9 | 亚洲精品av麻豆狂野| 叶爱在线成人免费视频播放| 国产精品免费视频内射| 又大又爽又粗| 丝袜在线中文字幕| 午夜精品久久久久久毛片777| videos熟女内射| 国产成人免费观看mmmm| 欧美变态另类bdsm刘玥| 搡老乐熟女国产| 亚洲中文av在线| 久久国产亚洲av麻豆专区| 亚洲,欧美精品.| 91大片在线观看| 人人妻人人爽人人添夜夜欢视频| 国产在线精品亚洲第一网站| 久久av网站| 亚洲午夜理论影院| 99精国产麻豆久久婷婷| 激情在线观看视频在线高清 | 制服诱惑二区| 一本大道久久a久久精品| 最近最新免费中文字幕在线| 精品少妇一区二区三区视频日本电影| 人人妻人人澡人人看| 十八禁高潮呻吟视频| 亚洲精品国产一区二区精华液| 精品一区二区三卡| 我的亚洲天堂| 免费观看人在逋| 亚洲精品av麻豆狂野| 国产亚洲欧美在线一区二区| 啦啦啦在线免费观看视频4| 日本黄色日本黄色录像| 热99国产精品久久久久久7| 在线观看免费日韩欧美大片| 老熟妇仑乱视频hdxx| 久久精品国产99精品国产亚洲性色 | www.熟女人妻精品国产| 亚洲精品国产精品久久久不卡| 中文字幕av电影在线播放| 91av网站免费观看| 精品福利观看| 久久精品成人免费网站| 老司机靠b影院| 成人国产av品久久久| 如日韩欧美国产精品一区二区三区| 老汉色∧v一级毛片| 成人国产av品久久久| 99热网站在线观看| 国产一区有黄有色的免费视频| 久久国产精品人妻蜜桃| 人人妻人人爽人人添夜夜欢视频| 成人永久免费在线观看视频 | 一二三四在线观看免费中文在| 中国美女看黄片| www.自偷自拍.com| 亚洲人成电影免费在线| av在线播放免费不卡| 亚洲欧洲日产国产| 日本精品一区二区三区蜜桃| 久久久精品94久久精品| 国产亚洲精品第一综合不卡| e午夜精品久久久久久久| 精品国产超薄肉色丝袜足j| 桃红色精品国产亚洲av| 精品久久久精品久久久| 日本精品一区二区三区蜜桃| 99热国产这里只有精品6| 黄色 视频免费看| 纵有疾风起免费观看全集完整版| 亚洲午夜精品一区,二区,三区| av超薄肉色丝袜交足视频| 国产成+人综合+亚洲专区| 日韩三级视频一区二区三区| 国产精品电影一区二区三区 | 狠狠精品人妻久久久久久综合| 亚洲av国产av综合av卡| 动漫黄色视频在线观看| 亚洲av第一区精品v没综合| 精品第一国产精品| 国产精品久久久av美女十八| 日本av免费视频播放| 欧美日韩亚洲高清精品| 精品福利永久在线观看| 一级黄色大片毛片| 国产精品av久久久久免费| 久久久国产一区二区| 国产成人系列免费观看| 99精品久久久久人妻精品| 久久久久久久久免费视频了| 国产av一区二区精品久久| 国产精品亚洲av一区麻豆| tube8黄色片| aaaaa片日本免费| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩福利视频一区二区| 国产精品98久久久久久宅男小说| 色精品久久人妻99蜜桃| 欧美人与性动交α欧美软件| 欧美激情 高清一区二区三区| 黄色怎么调成土黄色| 国产精品国产av在线观看| 国产高清videossex| 国产不卡一卡二| 桃花免费在线播放| 涩涩av久久男人的天堂| 女警被强在线播放| 99九九在线精品视频| 91成人精品电影| a在线观看视频网站| 麻豆乱淫一区二区| 人人妻人人澡人人看| 成年版毛片免费区| 一边摸一边抽搐一进一小说 | 1024视频免费在线观看| 女人爽到高潮嗷嗷叫在线视频| av有码第一页| 99精品在免费线老司机午夜| 亚洲久久久国产精品| av网站在线播放免费|