• <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)度算法
    国产熟女欧美一区二区| 亚洲五月色婷婷综合| 久久人人97超碰香蕉20202| 成人午夜精彩视频在线观看| 国产精品 国内视频| 日韩大片免费观看网站| 中文字幕高清在线视频| 韩国精品一区二区三区| 国产成人精品福利久久| 免费在线观看黄色视频的| 男女高潮啪啪啪动态图| 少妇人妻精品综合一区二区| 久久国产精品男人的天堂亚洲| 国产成人精品在线电影| 欧美亚洲 丝袜 人妻 在线| 午夜久久久在线观看| 中文字幕精品免费在线观看视频| 国产一区有黄有色的免费视频| 午夜福利免费观看在线| 日本vs欧美在线观看视频| 91精品伊人久久大香线蕉| 建设人人有责人人尽责人人享有的| 午夜影院在线不卡| 18禁国产床啪视频网站| 精品一品国产午夜福利视频| 久久久久国产一级毛片高清牌| 一区在线观看完整版| 日本午夜av视频| 性少妇av在线| 一级毛片电影观看| 午夜日韩欧美国产| 在线观看免费高清a一片| 最近中文字幕高清免费大全6| 免费观看性生交大片5| 美女福利国产在线| 女性被躁到高潮视频| 亚洲av福利一区| 国产精品三级大全| 王馨瑶露胸无遮挡在线观看| 狂野欧美激情性bbbbbb| 国产一区有黄有色的免费视频| 成年动漫av网址| 深夜精品福利| 国产伦理片在线播放av一区| 国产成人精品在线电影| 观看av在线不卡| 亚洲欧洲国产日韩| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看国产h片| 国产在视频线精品| 国产老妇伦熟女老妇高清| 一二三四在线观看免费中文在| 男女午夜视频在线观看| 亚洲成人国产一区在线观看 | 曰老女人黄片| 在线观看免费视频网站a站| 黄频高清免费视频| 欧美国产精品一级二级三级| a 毛片基地| 国产淫语在线视频| 汤姆久久久久久久影院中文字幕| 日韩成人av中文字幕在线观看| 久久久久久免费高清国产稀缺| 亚洲精品一二三| 亚洲一区二区三区欧美精品| 欧美日韩视频高清一区二区三区二| 一边摸一边做爽爽视频免费| 亚洲精品美女久久av网站| 亚洲色图 男人天堂 中文字幕| 精品国产一区二区久久| 热re99久久国产66热| 久久久久久人妻| 午夜福利免费观看在线| 大码成人一级视频| 2021少妇久久久久久久久久久| 国产精品久久久av美女十八| 又大又爽又粗| 黑人猛操日本美女一级片| 午夜激情av网站| 日韩欧美一区视频在线观看| 两性夫妻黄色片| 久久人人爽av亚洲精品天堂| 蜜桃国产av成人99| 美女脱内裤让男人舔精品视频| 色吧在线观看| 下体分泌物呈黄色| 人人妻人人爽人人添夜夜欢视频| 亚洲天堂av无毛| 国产99久久九九免费精品| 久久久精品区二区三区| 亚洲av电影在线进入| 国产麻豆69| 99九九在线精品视频| 亚洲伊人色综图| 大陆偷拍与自拍| 色吧在线观看| 久久人人爽人人片av| 午夜福利,免费看| 成人黄色视频免费在线看| av在线播放精品| 亚洲专区中文字幕在线 | 一本久久精品| 午夜福利乱码中文字幕| 免费人妻精品一区二区三区视频| 日日爽夜夜爽网站| 亚洲精品一区蜜桃| 久久天躁狠狠躁夜夜2o2o | 国产精品国产三级国产专区5o| 人人妻人人爽人人添夜夜欢视频| 国产成人免费观看mmmm| 欧美精品亚洲一区二区| 亚洲精品一二三| 一区在线观看完整版| 国产成人a∨麻豆精品| 日日啪夜夜爽| 高清黄色对白视频在线免费看| 婷婷色综合大香蕉| 亚洲欧洲精品一区二区精品久久久 | 水蜜桃什么品种好| 久久ye,这里只有精品| 亚洲美女视频黄频| 看免费av毛片| 国产精品免费视频内射| 波多野结衣av一区二区av| 另类亚洲欧美激情| 19禁男女啪啪无遮挡网站| 99国产精品免费福利视频| 亚洲欧美精品自产自拍| 亚洲天堂av无毛| 中文精品一卡2卡3卡4更新| 在线观看国产h片| 不卡av一区二区三区| 国产精品久久久久久久久免| 秋霞伦理黄片| 成人亚洲精品一区在线观看| 两性夫妻黄色片| 国产亚洲午夜精品一区二区久久| 一二三四中文在线观看免费高清| 国产黄色视频一区二区在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲成人国产一区在线观看 | 中文天堂在线官网| 最近中文字幕高清免费大全6| 男人舔女人的私密视频| 精品福利永久在线观看| 国产日韩一区二区三区精品不卡| 新久久久久国产一级毛片| 欧美在线一区亚洲| 亚洲精品av麻豆狂野| 91成人精品电影| 日韩一本色道免费dvd| 校园人妻丝袜中文字幕| 涩涩av久久男人的天堂| 丝袜人妻中文字幕| 精品久久久精品久久久| 人人妻人人爽人人添夜夜欢视频| 少妇被粗大猛烈的视频| 国产片特级美女逼逼视频| 国产不卡av网站在线观看| 亚洲天堂av无毛| 午夜福利影视在线免费观看| 欧美人与善性xxx| 日韩不卡一区二区三区视频在线| 97人妻天天添夜夜摸| 久久av网站| 精品国产一区二区久久| 建设人人有责人人尽责人人享有的| 看十八女毛片水多多多| 制服人妻中文乱码| 国产精品麻豆人妻色哟哟久久| 久久久亚洲精品成人影院| 一区二区三区四区激情视频| 巨乳人妻的诱惑在线观看| 大香蕉久久成人网| 免费日韩欧美在线观看| www.av在线官网国产| 中国三级夫妇交换| 国产精品久久久久久精品古装| 天天添夜夜摸| 欧美精品av麻豆av| 国产日韩欧美亚洲二区| 街头女战士在线观看网站| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说| av免费观看日本| 我要看黄色一级片免费的| 国产午夜精品一二区理论片| 亚洲av男天堂| 欧美 日韩 精品 国产| 成年美女黄网站色视频大全免费| 啦啦啦中文免费视频观看日本| 欧美精品高潮呻吟av久久| 一区二区三区精品91| 午夜91福利影院| 国产av一区二区精品久久| 欧美日韩一级在线毛片| 久久久久国产精品人妻一区二区| 欧美在线黄色| 美女脱内裤让男人舔精品视频| 男人操女人黄网站| 在线观看三级黄色| 咕卡用的链子| 国产亚洲av片在线观看秒播厂| 男人爽女人下面视频在线观看| 男女无遮挡免费网站观看| 这个男人来自地球电影免费观看 | 久久久精品区二区三区| 日韩视频在线欧美| 人体艺术视频欧美日本| 久久这里只有精品19| 毛片一级片免费看久久久久| 欧美 日韩 精品 国产| 欧美激情 高清一区二区三区| 亚洲av国产av综合av卡| 国产亚洲欧美精品永久| 最近手机中文字幕大全| av线在线观看网站| 大香蕉久久网| 国产黄频视频在线观看| 精品一区二区三区av网在线观看 | 天天操日日干夜夜撸| 国产99久久九九免费精品| 精品一区二区三卡| 国产一区二区在线观看av| a级毛片黄视频| 19禁男女啪啪无遮挡网站| 国产精品免费大片| tube8黄色片| 黄片无遮挡物在线观看| 97精品久久久久久久久久精品| 久久久国产精品麻豆| av一本久久久久| 国产乱人偷精品视频| 汤姆久久久久久久影院中文字幕| 在现免费观看毛片| 亚洲成色77777| 在线观看国产h片| 女人高潮潮喷娇喘18禁视频| 一区在线观看完整版| av女优亚洲男人天堂| 女性生殖器流出的白浆| 久久久久久久精品精品| 午夜免费男女啪啪视频观看| 又粗又硬又长又爽又黄的视频| 久久久久久人人人人人| 日韩大码丰满熟妇| 777久久人妻少妇嫩草av网站| 久久久久久久久久久免费av| 人人澡人人妻人| 久久国产精品男人的天堂亚洲| 少妇被粗大猛烈的视频| 国产无遮挡羞羞视频在线观看| 国产福利在线免费观看视频| 午夜91福利影院| 一本一本久久a久久精品综合妖精| 国产成人欧美在线观看 | 成人国产av品久久久| 亚洲美女搞黄在线观看| 天天操日日干夜夜撸| svipshipincom国产片| 91精品伊人久久大香线蕉| www.自偷自拍.com| 人人妻人人澡人人爽人人夜夜| 亚洲av成人精品一二三区| 美国免费a级毛片| 搡老岳熟女国产| av网站免费在线观看视频| 男女高潮啪啪啪动态图| 亚洲欧美中文字幕日韩二区| 1024香蕉在线观看| 久久性视频一级片| av卡一久久| tube8黄色片| 男的添女的下面高潮视频| 激情五月婷婷亚洲| 久久精品久久久久久久性| 在线天堂中文资源库| 亚洲欧洲国产日韩| 欧美成人精品欧美一级黄| 久久久久久久久久久久大奶| 日本欧美国产在线视频| 精品国产露脸久久av麻豆| 黄色怎么调成土黄色| 人妻 亚洲 视频| 最新在线观看一区二区三区 | 少妇精品久久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 妹子高潮喷水视频| 国产精品国产三级专区第一集| 国产老妇伦熟女老妇高清| 国产高清不卡午夜福利| 精品卡一卡二卡四卡免费| 亚洲精品美女久久久久99蜜臀 | 亚洲欧美一区二区三区国产| 中文天堂在线官网| 欧美 亚洲 国产 日韩一| 人人妻人人澡人人爽人人夜夜| 久久久久久久久免费视频了| 亚洲国产精品999| 亚洲免费av在线视频| 狂野欧美激情性bbbbbb| 欧美变态另类bdsm刘玥| 久久天躁狠狠躁夜夜2o2o | a级毛片黄视频| 如何舔出高潮| 久久精品国产亚洲av涩爱| 少妇被粗大猛烈的视频| 2021少妇久久久久久久久久久| 日韩欧美一区视频在线观看| 国产福利在线免费观看视频| 国产精品 欧美亚洲| 狂野欧美激情性bbbbbb| 亚洲,欧美精品.| 少妇人妻 视频| 免费观看性生交大片5| 韩国av在线不卡| 女性生殖器流出的白浆| 成人国产av品久久久| 成年女人毛片免费观看观看9 | 亚洲国产欧美一区二区综合| 亚洲精品国产av成人精品| 国产一区二区激情短视频 | 欧美少妇被猛烈插入视频| 波多野结衣av一区二区av| 亚洲av成人精品一二三区| 国产精品久久久久久久久免| 观看av在线不卡| 性色av一级| 亚洲精品视频女| 丝袜脚勾引网站| 欧美日韩一区二区视频在线观看视频在线| 卡戴珊不雅视频在线播放| 十八禁网站网址无遮挡| 精品少妇久久久久久888优播| 亚洲人成网站在线观看播放| 国产视频首页在线观看| 国产国语露脸激情在线看| 亚洲三区欧美一区| 久久性视频一级片| www日本在线高清视频| 久久久久精品人妻al黑| 久久精品久久久久久噜噜老黄| 久久精品亚洲熟妇少妇任你| 日本欧美视频一区| 久久精品亚洲av国产电影网| 国产免费一区二区三区四区乱码| 18禁国产床啪视频网站| 国产97色在线日韩免费| 亚洲精品在线美女| 纵有疾风起免费观看全集完整版| 丁香六月欧美| 久久精品aⅴ一区二区三区四区| 日本色播在线视频| 亚洲精品av麻豆狂野| 啦啦啦啦在线视频资源| 色网站视频免费| av一本久久久久| 久久精品aⅴ一区二区三区四区| av网站在线播放免费| 精品亚洲成a人片在线观看| 久久这里只有精品19| 自线自在国产av| 欧美日韩一级在线毛片| 亚洲美女搞黄在线观看| 国产成人91sexporn| av线在线观看网站| 极品人妻少妇av视频| 日韩制服丝袜自拍偷拍| 操美女的视频在线观看| 国产一区有黄有色的免费视频| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 国产成人免费观看mmmm| 搡老岳熟女国产| 大码成人一级视频| 精品国产一区二区三区四区第35| 欧美中文综合在线视频| 1024香蕉在线观看| 日韩成人av中文字幕在线观看| 在线精品无人区一区二区三| 中文字幕av电影在线播放| 岛国毛片在线播放| 国产成人a∨麻豆精品| 亚洲av日韩精品久久久久久密 | 热99国产精品久久久久久7| 黄频高清免费视频| 亚洲欧美中文字幕日韩二区| 久久免费观看电影| 亚洲熟女毛片儿| 街头女战士在线观看网站| 另类精品久久| 观看av在线不卡| 三上悠亚av全集在线观看| 又大又黄又爽视频免费| 97在线人人人人妻| 自拍欧美九色日韩亚洲蝌蚪91| 日韩 欧美 亚洲 中文字幕| 2018国产大陆天天弄谢| 青春草视频在线免费观看| 日本猛色少妇xxxxx猛交久久| 午夜福利乱码中文字幕| 亚洲国产精品一区三区| 日本黄色日本黄色录像| 精品少妇内射三级| 男女之事视频高清在线观看 | 建设人人有责人人尽责人人享有的| 亚洲精品,欧美精品| 亚洲精品国产区一区二| 亚洲欧洲国产日韩| 亚洲国产毛片av蜜桃av| 国产精品 欧美亚洲| 1024香蕉在线观看| 啦啦啦中文免费视频观看日本| 亚洲中文av在线| 国产日韩一区二区三区精品不卡| 精品国产超薄肉色丝袜足j| 欧美国产精品va在线观看不卡| 亚洲成人国产一区在线观看 | 日本vs欧美在线观看视频| 激情视频va一区二区三区| 日韩人妻精品一区2区三区| 嫩草影院入口| 亚洲欧美成人精品一区二区| 熟女少妇亚洲综合色aaa.| 99久久人妻综合| 一区二区三区精品91| 午夜激情av网站| 精品福利永久在线观看| 女性被躁到高潮视频| 男女之事视频高清在线观看 | 亚洲精品国产一区二区精华液| 丝袜喷水一区| 国产一区二区在线观看av| 亚洲精品国产av蜜桃| 欧美亚洲 丝袜 人妻 在线| 97人妻天天添夜夜摸| av.在线天堂| 啦啦啦视频在线资源免费观看| 日韩欧美一区视频在线观看| 日韩制服骚丝袜av| 秋霞在线观看毛片| 亚洲欧美一区二区三区黑人| xxxhd国产人妻xxx| 成人免费观看视频高清| 欧美 亚洲 国产 日韩一| 老鸭窝网址在线观看| 黄网站色视频无遮挡免费观看| 免费观看a级毛片全部| 亚洲av成人不卡在线观看播放网 | 国产一区二区三区av在线| 午夜精品国产一区二区电影| 一本色道久久久久久精品综合| 亚洲国产看品久久| 亚洲情色 制服丝袜| 亚洲中文av在线| 亚洲欧美清纯卡通| 老鸭窝网址在线观看| 啦啦啦啦在线视频资源| 啦啦啦视频在线资源免费观看| 欧美最新免费一区二区三区| 9热在线视频观看99| 亚洲精品,欧美精品| 男人舔女人的私密视频| 中文精品一卡2卡3卡4更新| kizo精华| 午夜日韩欧美国产| 男女下面插进去视频免费观看| 亚洲熟女毛片儿| 日本一区二区免费在线视频| 国产成人系列免费观看| 成人18禁高潮啪啪吃奶动态图| av有码第一页| 久久久久网色| 亚洲熟女精品中文字幕| 日本91视频免费播放| 男女高潮啪啪啪动态图| 大香蕉久久网| 在线 av 中文字幕| 如日韩欧美国产精品一区二区三区| 99久久人妻综合| 高清av免费在线| 人妻一区二区av| 热re99久久国产66热| 国产福利在线免费观看视频| 超碰97精品在线观看| 色婷婷av一区二区三区视频| 欧美变态另类bdsm刘玥| 卡戴珊不雅视频在线播放| 男女免费视频国产| 丰满饥渴人妻一区二区三| 久久精品久久久久久噜噜老黄| 一级爰片在线观看| 久久久久久久久免费视频了| 国产精品av久久久久免费| 一级a爱视频在线免费观看| 欧美日韩视频精品一区| av在线app专区| 男女国产视频网站| 国产一区二区激情短视频 | 亚洲欧美一区二区三区久久| 亚洲国产欧美一区二区综合| 亚洲精品国产av蜜桃| 亚洲av福利一区| 男女床上黄色一级片免费看| 亚洲精品国产av成人精品| 亚洲av国产av综合av卡| 亚洲av在线观看美女高潮| 亚洲欧洲日产国产| 成年人免费黄色播放视频| 一本大道久久a久久精品| 国产精品国产三级专区第一集| www.自偷自拍.com| 国产精品 欧美亚洲| 国产成人精品福利久久| 亚洲欧美中文字幕日韩二区| 看非洲黑人一级黄片| 日韩视频在线欧美| 亚洲精品久久午夜乱码| 亚洲av日韩在线播放| 免费看av在线观看网站| 少妇的丰满在线观看| 亚洲精品久久成人aⅴ小说| 我的亚洲天堂| 制服丝袜香蕉在线| 蜜桃国产av成人99| av又黄又爽大尺度在线免费看| 国产成人一区二区在线| 最近手机中文字幕大全| 中文欧美无线码| 亚洲美女搞黄在线观看| 久久久精品94久久精品| 免费高清在线观看日韩| 久久人人爽人人片av| 性少妇av在线| a级片在线免费高清观看视频| 精品卡一卡二卡四卡免费| 视频在线观看一区二区三区| 精品亚洲乱码少妇综合久久| 在线观看国产h片| 美女脱内裤让男人舔精品视频| 麻豆乱淫一区二区| 最黄视频免费看| 日韩电影二区| 午夜福利视频精品| 亚洲精品自拍成人| 十八禁高潮呻吟视频| 国产激情久久老熟女| 亚洲视频免费观看视频| 色综合欧美亚洲国产小说| 在线观看国产h片| 国产又色又爽无遮挡免| 欧美精品av麻豆av| 国产午夜精品一二区理论片| 免费久久久久久久精品成人欧美视频| 赤兔流量卡办理| 校园人妻丝袜中文字幕| 黑丝袜美女国产一区| av电影中文网址| 制服丝袜香蕉在线| 韩国av在线不卡| 看非洲黑人一级黄片| 亚洲av综合色区一区| 大香蕉久久成人网| 18禁动态无遮挡网站| bbb黄色大片| 国产色婷婷99| 在线观看免费高清a一片| 亚洲精品自拍成人| 日本91视频免费播放| 日本av手机在线免费观看| 午夜福利网站1000一区二区三区| 欧美日韩亚洲高清精品| 免费av中文字幕在线| 精品国产一区二区三区久久久樱花| 在线观看www视频免费| 777米奇影视久久| 五月开心婷婷网| 97精品久久久久久久久久精品| 国产成人免费观看mmmm| 99国产精品免费福利视频| 久久久久精品性色| 亚洲国产av新网站| 亚洲第一青青草原| 欧美 亚洲 国产 日韩一| 欧美日韩一级在线毛片| 欧美变态另类bdsm刘玥| 一本久久精品| 99香蕉大伊视频| 国产精品国产av在线观看| 一本久久精品| 观看av在线不卡| 国产午夜精品一二区理论片| 中文字幕精品免费在线观看视频| 别揉我奶头~嗯~啊~动态视频 | 亚洲成人一二三区av| 国产男女内射视频| 欧美激情极品国产一区二区三区| 少妇人妻精品综合一区二区| 亚洲精品乱久久久久久| 岛国毛片在线播放| 国产亚洲av高清不卡| av又黄又爽大尺度在线免费看| 五月开心婷婷网| 亚洲三区欧美一区| 国产亚洲欧美精品永久| 午夜福利在线免费观看网站| 妹子高潮喷水视频| www.精华液| 在线天堂最新版资源| 在线亚洲精品国产二区图片欧美| 国产人伦9x9x在线观看| 久久精品国产亚洲av高清一级| 自线自在国产av| 黄片播放在线免费| 欧美在线一区亚洲|