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

    金沙江中下游浮游植物功能群與環(huán)境因子的關(guān)系

    2023-04-29 00:50:35龔進(jìn)玲朱挺兵楊德國(guó)胡飛飛杜紅春李學(xué)梅孟子豪吳興兵
    水生態(tài)學(xué)雜志 2023年3期
    關(guān)鍵詞:浮游植物環(huán)境因子金沙江

    龔進(jìn)玲 朱挺兵 楊德國(guó) 胡飛飛 杜紅春 李學(xué)梅 孟子豪 吳興兵

    摘要:分析不同季節(jié)浮游植物群落結(jié)構(gòu)組成,明確主要環(huán)境影響因子,為金沙江干流水生生物保護(hù)和生態(tài)修復(fù)提供理論依據(jù)。于2019年春夏秋三季對(duì)金沙江中下游(河流水體和庫(kù)區(qū))浮游植物群落組成、功能群以及多樣性指數(shù)進(jìn)行調(diào)查,并分析了浮游植物功能群與環(huán)境因子的關(guān)系。結(jié)果表明:金沙江中下游共鑒定出浮游植物8門(mén)32科54屬96種,以硅藻門(mén)(39.58%)為主,其次是綠藻門(mén)(31.25%),河流水體與庫(kù)區(qū)在浮游植物群落組成上差異不顯著;Shannon-Wiener指數(shù)、Margalef指數(shù)和Pielou指數(shù)平均值分別為2.42、1.81和0.92,單因素方差分析表明,3項(xiàng)多樣性指數(shù)在各采樣點(diǎn)之間、各季節(jié)之間無(wú)顯著差異;浮游植物共劃分出24種功能群,其中8種(D、P、X1、J、C、LM、F、G)為該區(qū)域的優(yōu)勢(shì)功能群。河流水體出現(xiàn)頻率大于30%的浮游植物功能群有12種,主要優(yōu)勢(shì)功能群以D和P為主;庫(kù)區(qū)有15種,主要優(yōu)勢(shì)功能群以D、J、P和X1為主。Pearson相關(guān)性分析和冗余分析(RDA)結(jié)果顯示,溫度、濁度、海拔、氨氮和硝態(tài)氮是影響金沙江中下游流域浮游植物功能群分布格局的主要因子。

    關(guān)鍵詞:浮游植物;功能群;環(huán)境因子;金沙江

    中圖分類(lèi)號(hào):S91? ? ? ? 文獻(xiàn)標(biāo)志碼:A? ? ? ? 文章編號(hào):1674-3075(2023)03-0086-08

    浮游植物作為水體中最重要的初級(jí)生產(chǎn)者,其生產(chǎn)力幾乎占全球總初級(jí)生產(chǎn)力的50%(Field et al,1998),是水生態(tài)系統(tǒng)中物質(zhì)和能量流動(dòng)的驅(qū)動(dòng)因子。浮游植物群落結(jié)構(gòu)與環(huán)境因子關(guān)系密切,能夠敏感地反映水生態(tài)系統(tǒng)對(duì)水環(huán)境污染和生態(tài)破壞的響應(yīng),是水質(zhì)評(píng)價(jià)的重要指示生物。由于傳統(tǒng)的藻類(lèi)分類(lèi)法未考慮藻類(lèi)對(duì)環(huán)境響應(yīng)的物種特異性和環(huán)境特異性,導(dǎo)致水質(zhì)評(píng)價(jià)結(jié)果存在差異(Reynolds et al,2002)。為了進(jìn)一步探究藻類(lèi)群落結(jié)構(gòu)變化與環(huán)境因子之間的關(guān)系,Reynolds等(2002)提出了對(duì)淡水浮游植物進(jìn)行功能群劃分,將相同生態(tài)環(huán)境下能共同生長(zhǎng)的浮游植物歸為同一功能群,以便最大限度揭示浮游植物與環(huán)境因子的動(dòng)態(tài)關(guān)系。Reynolds等(2002)和Padisák等(2009)確定了40組不同的功能群(Functional Group,F(xiàn)G),其中有39組功能群應(yīng)用廣泛。近年來(lái),浮游植物功能群分類(lèi)方法逐漸成為國(guó)內(nèi)外學(xué)者研究的熱點(diǎn),被廣泛應(yīng)用于群落生態(tài)學(xué)研究領(lǐng)域(王汨等,2020;Ma et al,2019;趙耿楠等,2021)。

    金沙江屬長(zhǎng)江上游江段,起于青海省玉樹(shù)州直門(mén)達(dá),止于四川省宜賓市長(zhǎng)江岷江交匯處,是長(zhǎng)江的重要組成部分。受大規(guī)模水資源開(kāi)發(fā)利用影響,金沙江流域生態(tài)系統(tǒng)面臨著諸多問(wèn)題,如梯級(jí)開(kāi)發(fā)對(duì)流水生境的破壞、沿江排污、過(guò)度捕撈等。提高對(duì)金沙江水環(huán)境保護(hù)的關(guān)注度,對(duì)維持整個(gè)長(zhǎng)江水生態(tài)系統(tǒng)健康有著重要的意義(和曉榮等,2015)。目前關(guān)于金沙江浮游植物的研究主要集中于長(zhǎng)江源區(qū)、金沙江干支流及其附屬水庫(kù)(崔玉香等,2020;魏志兵等,2020;馬寶珊等,2020;陳儉勇等,2021)。金沙江中下游是水電開(kāi)發(fā)的重要區(qū)域,有關(guān)中下游浮游植物群落結(jié)構(gòu)及其與水環(huán)境關(guān)系的研究較少,本研究以金沙江中下游(河流水體和庫(kù)區(qū))為研究對(duì)象,分析不同季節(jié)浮游植物群落結(jié)構(gòu)組成,劃分功能群通過(guò)Pearson和RDA分析其與水環(huán)境因子之間的相關(guān)性,明確影響該區(qū)域浮游植物功能群的主要環(huán)境因子,為金沙江干流水生生物保護(hù)和生態(tài)修復(fù)提供理論依據(jù)。

    1? ?材料與方法

    1.1? ?采樣點(diǎn)布設(shè)

    參考《水域生態(tài)系統(tǒng)觀測(cè)規(guī)范》(蔡慶華等,2007)并根據(jù)實(shí)際地形地貌結(jié)合采樣的可行性,在金沙江中下游設(shè)置了11個(gè)采樣點(diǎn),并通過(guò)水流速度將采樣點(diǎn)劃分為河流水體(>0.3 m/s)和庫(kù)區(qū)(0~0.03 m/s)。其中河流水體7個(gè),分別為金安橋(JA)、龍開(kāi)口(LK)、雅江匯口上游(YS)、雅江匯口下游(YX)、皎平渡(JP)、巧家(QJ)和水富(SF),庫(kù)區(qū)4個(gè),分別為魯?shù)乩↙D)、觀音巖(GY)、溪洛渡(XL)和新市(XS)。具體采樣斷面見(jiàn)圖1。

    1.2? ?樣品采集與處理

    于2019年春季(4月)、夏季(8月)和秋季(11月)對(duì)金沙江中下游浮游植物和水環(huán)境參數(shù)進(jìn)行調(diào)查。2019年冬季因新冠疫情未能開(kāi)展調(diào)查。

    浮游植物的定性、定量樣品采集及保存方法參考文獻(xiàn)(章宗涉和黃祥飛,1991)。用1 L采水器在相同水層取水樣10 L,混合均勻后取1 L水樣作為定量樣本,加入魯哥試劑固定,靜置48 h后濃縮至30 mL備檢。采集的樣品通過(guò)奧林巴斯顯微鏡(CX41)進(jìn)行觀察,浮游植物種類(lèi)鑒定和計(jì)數(shù)參考相關(guān)資料(胡鴻鈞和魏印心,2006;章宗涉和黃祥飛,1991)。

    采集樣本的同時(shí)現(xiàn)場(chǎng)對(duì)各采樣斷面的環(huán)境參數(shù)進(jìn)行檢測(cè),包括海拔(EL)、透明度(SD)、濁度(TUR)、溶氧(DO)、pH、和水溫(WT)等參數(shù);實(shí)驗(yàn)室內(nèi)分別對(duì)總磷(TP)、總氮(TN)、氨氮(NH4+-N)、亞硝態(tài)氮(NO2--N)、硝態(tài)氮(NO3--N)和磷酸鹽(PO43--P)等指標(biāo)進(jìn)行檢測(cè),所有指標(biāo)的測(cè)定參照《水與廢水監(jiān)測(cè)分析方法》(國(guó)家環(huán)保局,1989)中的方法進(jìn)行。

    1.3? ?數(shù)據(jù)處理與分析

    浮游植物群落優(yōu)勢(shì)種由Mcnaughton優(yōu)勢(shì)種指數(shù)(Y)確定優(yōu)勢(shì)種(Y[≥]0.02);浮游植物群落多樣性由Shannon-Wiener多樣性指數(shù)(H')、Margalef種類(lèi)豐富度指數(shù)(D)和Pielou均勻度指數(shù)(J)進(jìn)行分析(Belaoussoff et al,2003;Lampitt et al,1993),公式如下:

    式中:fi是第i物種的出現(xiàn)頻率,Pi=ni/N,是第i物種個(gè)數(shù)在總個(gè)數(shù)的比例,S為種類(lèi)數(shù),N為總密度。

    參照Reynolds等(2002)和Padisák等(2009)對(duì)浮游植物功能群劃分的標(biāo)準(zhǔn),對(duì)具有相同適應(yīng)性特征、易于在相同生境條件下共存的浮游植物種類(lèi)進(jìn)行功能分組。

    利用SPSS 20.0對(duì)數(shù)據(jù)進(jìn)行單因素方差分析(One-way ANOVA),采用Duncans多重比較進(jìn)行差異顯著性檢驗(yàn);采用Pearson相關(guān)性分析對(duì)浮游生物群落結(jié)構(gòu)與環(huán)境因子進(jìn)行分析;采用 CANOCO 5. 0 軟件對(duì)浮游生物和水環(huán)境因子進(jìn)行冗余分析(Redundancy Analysis,RDA)。

    2? ?結(jié)果與分析

    2.1? ?種類(lèi)組成

    共檢測(cè)到浮游植物8門(mén)32科54屬96種。種類(lèi)組成以硅藻門(mén)(38種)、綠藻門(mén)(30種)和藍(lán)藻門(mén)(13種)為主,共占總種類(lèi)數(shù)的84.37%,裸藻門(mén)(7種)、甲藻門(mén)(3種)、黃藻門(mén)(3種)、隱藻門(mén)(1種)和金藻門(mén)(1種)物種數(shù)較少,占總種類(lèi)數(shù)的15.63%(圖2)。

    從空間分布上看,各采樣點(diǎn)浮游植物種類(lèi)較為豐富,巧家和新市最高,為44種,河流水體的龍開(kāi)口最少,為36種。河流水體的7個(gè)采樣點(diǎn)物種平均數(shù)為40種,略低于4個(gè)庫(kù)區(qū)平均數(shù)42種。從季節(jié)上看,春秋兩季種類(lèi)數(shù)最多,各70種;夏季較少,為56種(圖3)。浮游植物種類(lèi)組成均以硅藻門(mén)為主,占比多超過(guò)50%。

    2.2? ?浮游植物功能群

    2019年金沙江中下游浮游植物可劃分為24個(gè)功能群(表1):D、C、P、E、F、G、H1、J、K、Lo、LM、MP、N、S2、TC、W1、W2、X1、X2、X3、Y。河流水體共鑒定出24個(gè)功能群,其中功能群D和P的出現(xiàn)頻率最高,為100%,在所有采樣點(diǎn)全年均有出現(xiàn),出現(xiàn)頻率大于30%的功能群有12種;而庫(kù)區(qū)共鑒定出23個(gè)功能群,其中功能群D的出現(xiàn)率為100%,出現(xiàn)頻率大于30%的功能群有15種。

    2019年金沙江中下游流域浮游植物功能群相對(duì)生物量情況見(jiàn)圖4。按采樣點(diǎn)相對(duì)生物量>20%的功能群定義為優(yōu)勢(shì)功能群,共有8個(gè),分別為D、P、X1、J、C、LM、F和G。從空間分布上看,河流水體主要以D和P為主,所有采樣點(diǎn)相對(duì)生物量平均后,在優(yōu)勢(shì)功能群中分別占比22.83%和31.26%,相對(duì)比例最高的點(diǎn)位分別在巧家(34.05%)和金安橋(44.46%)。庫(kù)區(qū)D、J、P和X1優(yōu)勢(shì)度相對(duì)較大,在優(yōu)勢(shì)功能群中分別占比21.23%、16.75%、15.38%和12.89%,各點(diǎn)位中最高相對(duì)比例分別在魯?shù)乩?4.99%)、新市(33.88%)、觀音巖(18.39%)和魯?shù)乩?0.07%)。從季節(jié)上看,春季以D+P+C功能群類(lèi)為主導(dǎo),相對(duì)比例之和為51.79%;夏季C功能群消退,J功能群大量生長(zhǎng)繁殖,D+P+J功能類(lèi)群逐漸發(fā)展并占優(yōu)勢(shì);秋季X1功能群取代J功能群,以D+P+X1為主導(dǎo),相對(duì)比例之和為64.32%。綜上所述,金沙江不同采樣點(diǎn)的浮游植物功能群存在一定差異,D和P功能群是金沙江中下游各樣點(diǎn)常年共有的優(yōu)勢(shì)種功能群,這些功能群主要分布在渾濁和中富營(yíng)養(yǎng)的水體中。

    2.3? ?浮游植物多樣性變化

    金沙江中下游浮游植物Shannon-Wiener多樣性指數(shù)(H')為1.49~2.84,均值2.42。Pielou均勻度指數(shù)(J)的變化范圍為0.83~0.98,均值為0.92。Marglef豐富度指數(shù)(D)為0.83~2.62,均值1.81。金沙江中下游河流水體的3項(xiàng)指數(shù)平均值較庫(kù)區(qū)低(圖5)。通過(guò)單因素方差分析比較金沙江中下游時(shí)空差異,結(jié)果表明,3類(lèi)多樣性指數(shù)在各采樣點(diǎn)之間、各季節(jié)之間無(wú)顯著差異。

    2.4? ?浮游植物功能群與環(huán)境因子關(guān)系

    2.4.1? ?Pearson相關(guān)性分析? ?金沙江中下游8個(gè)優(yōu)勢(shì)功能群生物量與水環(huán)境因子的相關(guān)性分析結(jié)果表明(表2),D功能群與硝態(tài)氮呈極顯著負(fù)相關(guān)(R2=-0.476,P<0.01),與海拔呈極顯著正相關(guān)(R2=0.448,P<0.01);X1功能群與氨氮呈顯著負(fù)相關(guān)(R2=-0.379,P<0.05);F功能群與濁度呈極顯著負(fù)相關(guān)(R2=-0.480,P<0.01);G功能群與總氮呈顯著正相關(guān)(R2=0.426,P<0.05);LM功能群與氨氮呈極顯著正相關(guān)(R2=0.536,P<0.01)。總體來(lái)說(shuō),功能群浮游植物優(yōu)勢(shì)功能群生物量與濁度、海拔、氨氮、硝態(tài)氮和總氮等環(huán)境因子的關(guān)系密切。

    2.4.2? ?RDA分析? ?對(duì)12項(xiàng)水環(huán)境因子進(jìn)行主成分分析,提取特征值大于1的主成分,共有5個(gè),累積貢獻(xiàn)率76.74%。選取載荷量超過(guò)0.4的指標(biāo)作為主要影響因子,篩選出5項(xiàng)環(huán)境因子指標(biāo),分別是NH4+-N(0.633)、TN(0.626)、WT(0.614)、TUR(0.611)和PO43--P(0.415)。對(duì)浮游植物優(yōu)勢(shì)功能群與5項(xiàng)水環(huán)境因子進(jìn)行DCA分析,結(jié)果顯示最大的梯度長(zhǎng)度為1.8,小于3,因此選擇對(duì)浮游植物優(yōu)勢(shì)功能群與水環(huán)境因子做RDA分析。分析結(jié)果表明,軸一與軸二共解釋了優(yōu)勢(shì)功能群數(shù)據(jù)累積變化率19.6%(表3),優(yōu)勢(shì)功能群與環(huán)境因子的相關(guān)性分別為0.662和0.550,前兩軸累積方差為79.0%,表明排序能夠較好的反映浮游植物優(yōu)勢(shì)功能群與環(huán)境因子的關(guān)系。

    RDA分析排序(圖6)顯示,對(duì)優(yōu)勢(shì)功能群影響較大的環(huán)境因子是水溫和氨氮。其中排序軸1與氨氮(0.4863)的正相關(guān)性較大,其次是水溫(0.4419)。排序軸2與水溫正相關(guān)性最大(0.3477),與氨氮呈負(fù)相關(guān)較大(-0.2720)。優(yōu)勢(shì)功能群在各象限均存在,其中LM、P與氨氮、濁度及總氮呈正相關(guān);C與濁度呈正相關(guān);G、J、D與水溫呈較強(qiáng)正相關(guān)性;F、X1與所有水環(huán)境因子呈負(fù)相關(guān)。

    3? ?討論

    3.1? ?浮游植物功能群及其群落結(jié)構(gòu)特征

    2019年金沙江中下游共發(fā)現(xiàn)浮游植物96種,與魏志兵等(2020)的調(diào)查結(jié)果(100種)接近,但低于龔?fù)⒌堑龋?014)在金沙江水富段的調(diào)查結(jié)果(124種),造成這種差異可能與采樣點(diǎn)不完全一致、調(diào)查時(shí)間相隔久遠(yuǎn)等因素有關(guān)。本研究中,浮游植物種類(lèi)組成主要集中在硅藻、綠藻和藍(lán)藻,這與往年關(guān)于長(zhǎng)江干流的研究結(jié)果相似(魏志兵等,2020;王岳等,2015)。金沙江中下游河流水體與庫(kù)區(qū)的浮游植物種類(lèi)數(shù)差異不大,4個(gè)庫(kù)區(qū)平均數(shù)略高于河流水體,這表明金沙江中下游梯級(jí)水電建設(shè)對(duì)浮游植物的影響是一個(gè)相對(duì)緩慢長(zhǎng)期的過(guò)程。同時(shí)季節(jié)上差異明顯,春秋兩季種類(lèi)數(shù)較夏季豐富,可能是夏季處于豐水期,河流內(nèi)水體流速較快,含沙量高且透明度低,使水體中藻類(lèi)減少。

    通過(guò)Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和Marglef豐富度指數(shù)對(duì)水質(zhì)級(jí)別進(jìn)行劃分,結(jié)果顯示金沙江中下游的水體基本處于中污染,這與魏志兵等(2020)2017-2018年對(duì)金沙江干流流域研究結(jié)果一致,說(shuō)明水質(zhì)并未進(jìn)一步惡化。

    金沙江中下游共劃分24種浮游植物功能群,其中D、P、X1、J、C、LM、F和G為優(yōu)勢(shì)功能群。河流水體出現(xiàn)的24種功能群中只有12種的出現(xiàn)頻率超過(guò)30%,其中以適應(yīng)透明度低、濁度高的功能群D和P為主要優(yōu)勢(shì)功能群。金沙江干流含沙量高,光照強(qiáng)度弱的特點(diǎn)剛好適合這些類(lèi)群的生長(zhǎng)繁殖(王岳等,2015)。庫(kù)區(qū)共鑒定23種功能群,其中出現(xiàn)頻率大于30%的有15種。其種數(shù)稍低于河流水體,但從出現(xiàn)頻率上看,明顯高于河流水體,這可能是大壩攔截增加了水庫(kù)水體滯留時(shí)間和營(yíng)養(yǎng)鹽累計(jì)效率,使水庫(kù)產(chǎn)生熱力學(xué)分層,從而促進(jìn)浮游植物快速繁殖(申恒倫等,2012)。庫(kù)區(qū)主要演替功能群為D、J、P和X1,這些功能群均適應(yīng)中-富營(yíng)養(yǎng)水體,本研究中水體偏堿性,營(yíng)養(yǎng)鹽含量高,有利于這些類(lèi)群的生長(zhǎng)。浮游植物功能群在3個(gè)季節(jié)上呈P+D+C→D+P+J→P+D+X1的變化趨勢(shì),D和P功能群占絕對(duì)優(yōu)勢(shì),所指示的生境特征表明金沙江中下游為渾濁的富營(yíng)養(yǎng)型水體。春季水溫低,光照強(qiáng)度弱,適宜C功能群生長(zhǎng)繁殖。夏季C功能群逐漸消退,J功能群成為優(yōu)勢(shì)功能群。金沙江干流7-8月降雨豐沛,水體混合均勻,水體中TN和NH4+-N的含量較高,豐富的營(yíng)養(yǎng)元素有利于J功能群的生長(zhǎng)繁殖(張俊芳等,2021)。X1功能群對(duì)水體混合敏感,適應(yīng)生境為營(yíng)養(yǎng)型-超富營(yíng)養(yǎng)型水體。秋季降雨量小,水體穩(wěn)定性較好,適合其生長(zhǎng)繁殖。

    3.2? ?浮游植物群落結(jié)構(gòu)與水環(huán)境因子的關(guān)系

    浮游植物功能群與水體生境特征相對(duì)應(yīng),溫度、光照、水體動(dòng)力學(xué)特征、營(yíng)養(yǎng)鹽狀況、浮游動(dòng)物覓食壓力等因素能在特定生境中影響浮游植物功能群的組成(黃國(guó)佳等,2015)。研究表明,不同種類(lèi)的浮游植物適宜生長(zhǎng)的溫度不同,而硅藻一般適宜較低的溫度(張歡等,2017)。這與金沙江中下游浮游植物優(yōu)勢(shì)功能群多為硅藻的現(xiàn)象一致。RDA排序顯示,水溫是影響浮游浮游植物功能群較顯著的環(huán)境變量,這與漢江(杜紅春,2020)、湯浦水庫(kù)(張俊芳等,2021)等研究結(jié)果相似。與水溫直接關(guān)聯(lián)的海拔對(duì)浮游植物生長(zhǎng)影響較大。Pearson相關(guān)性分析中優(yōu)勢(shì)功能群D與海拔呈極顯著正相關(guān),表明海拔是影響浮游植物功能群的關(guān)鍵因子之一。研究表明,水體中顆粒物質(zhì)對(duì)光的吸收和反射會(huì)導(dǎo)致進(jìn)入水體的光照發(fā)生衰減,浮游植物等的生產(chǎn)力將會(huì)將低(Obrador & Pretus,2008)。金沙江中下游被水庫(kù)分割成幾段不連續(xù)的河流水體,自然河道中含沙量較大,加上該區(qū)域內(nèi)采砂活動(dòng)多,造就了高濁度的生態(tài)環(huán)境。相對(duì)靜止的庫(kù)區(qū)在流水的作用下發(fā)生上下層交換,導(dǎo)致顆粒物增加,使浮游植物的生長(zhǎng)受到一定程度的限制。Pearson相關(guān)性分析表明,優(yōu)勢(shì)功能群F與濁度呈極顯著負(fù)相關(guān),與同為長(zhǎng)江水系的西江(陶峰等,2019)浮游植物的研究相一致,說(shuō)明濁度在一定程度上影響浮游植物功能群。氮、磷是浮游植物生長(zhǎng)必須的營(yíng)養(yǎng)鹽,對(duì)浮游植物的密度、生物量和種類(lèi)組成有顯著的影響(劉林峰等,2018)。RDA排序圖中,第一軸、第二軸與氨氮相關(guān)性較大,表明氮對(duì)金沙江中下游浮游植物功能群的影響較大。相關(guān)分析中,D、LM、P優(yōu)勢(shì)功能群與氨氮呈正相關(guān),而部分功能群與硝態(tài)氮呈負(fù)相關(guān)。金沙江中下游隨著海拔降低,城市人口活動(dòng)較密集,農(nóng)業(yè)污染和生活廢水逐漸增多,導(dǎo)致水體營(yíng)養(yǎng)鹽濃度升高。有研究表明,大氣沉降也是水體氮的最主要來(lái)源(Pitk?nen et al,2001)。降雨的沖刷和攪拌提高了水體中氮的濃度,因此氮的影響比較復(fù)雜。

    綜上所述,浮游植物功能群演替受不同環(huán)境因子影響。總體而言,水溫、濁度、海拔、氨氮、硝態(tài)氮是金沙江中下游浮游植物功能群演替的關(guān)鍵驅(qū)動(dòng)因子。

    參考文獻(xiàn)

    蔡慶華,曹明,黃祥飛,2007. 水域生態(tài)系統(tǒng)觀測(cè)規(guī)范[M].北京:中國(guó)環(huán)境科學(xué)出版社.

    陳儉勇,孫婷婷,李紅濤,等,2021. 金沙江上游蘇洼龍水電站建設(shè)期影響區(qū)域浮游植物群落結(jié)構(gòu)與水質(zhì)現(xiàn)狀分析[J].淡水漁業(yè),51(2):87-97.

    崔玉香,張靜,柴元冰,等,2020. 長(zhǎng)江源區(qū)浮游植物群落特征研究[J].環(huán)境科學(xué)與技術(shù),43(1):20-25.

    杜紅春,2020. 漢江干流浮游生物群落結(jié)構(gòu)和功能群特征及水質(zhì)評(píng)價(jià)[D]. 武漢:華中農(nóng)業(yè)大學(xué).

    國(guó)家環(huán)保局,1989. 水和廢水監(jiān)測(cè)分析方法[M].北京:中國(guó)環(huán)境科學(xué)出版社.

    龔?fù)⒌?,楊偉階,何滔,等,2014. 金沙江水富段水生生物的群落結(jié)構(gòu)及水質(zhì)監(jiān)測(cè)[J].淡水漁業(yè), (4):25-34,22.

    和曉榮,王秀英,和武,2015. 金沙江中游庫(kù)區(qū)水污染綜合防治措施研究[J].水利水電技術(shù), 506(12):76-81.

    胡鴻鈞,魏印心,2006. 中國(guó)淡水藻類(lèi)——系統(tǒng)、分類(lèi)及生態(tài)[M].北京:科學(xué)出版社.

    黃國(guó)佳,李秋華,陳椽,等,2015. 貴州高原三板溪水庫(kù)浮游植物功能群時(shí)空分布特征[J].環(huán)境科學(xué)學(xué)報(bào), 35(2):418-428.

    劉林峰,周先華,高健,等,2018. 神農(nóng)架大九湖濕地浮游植物群落結(jié)構(gòu)特征及營(yíng)養(yǎng)狀態(tài)評(píng)價(jià)[J]. 湖泊科學(xué),30(2):417-430.

    馬寶珊,徐濱,魏開(kāi)金,等,2020. 安寧河中游浮游植物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系[J].生態(tài)學(xué)雜志,39(10):3332-3341.

    申恒倫,蔡慶華,邵美玲,等,2012. 三峽水庫(kù)香溪河流域梯級(jí)水庫(kù)浮游植物群落結(jié)構(gòu)特征[J].湖泊科學(xué),24(2):197-205.

    陶峰,劉凱,藺丹清,等,2019. 安慶市西江浮游植物群落多樣性時(shí)空格局[J].大連海洋大學(xué)學(xué)報(bào), 34(2):90-97.

    王汨,楊柏賀,馬思琦,等,2020. 北運(yùn)河水系浮游植物功能群與水體環(huán)境因子的關(guān)系[J].水產(chǎn)學(xué)雜志,33(2):70-76.

    魏志兵,何勇鳳,龔進(jìn)玲,等,2020. 金沙江干流浮游植物群落結(jié)構(gòu)特征及其時(shí)空變化[J].長(zhǎng)江流域資源與環(huán)境,29(6):1356-1365.

    王岳,夏爽,裴國(guó)鳳,等,2015. 長(zhǎng)江干流湖北段浮游藻類(lèi)群落結(jié)構(gòu)特征[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 31(6):916-922.

    趙耿楠,潘保柱,丁一桐,等,2021. 渭河干流和秦嶺北麓典型支流浮游植物功能群特征及水質(zhì)評(píng)價(jià)[J].生態(tài)學(xué)報(bào),41(8):3226-3237.

    章宗涉,黃祥飛,1991. 淡水浮游生物研究方法[M].北京:科學(xué)出版社.

    張俊芳,胡曉紅,馬沛明,等,2021. 湯浦水庫(kù)浮游植物功能群季節(jié)演替及關(guān)鍵驅(qū)動(dòng)因子[J].水生態(tài)學(xué)雜志,42(3):55-62.

    張歡,張佳磊,劉德富,等,2017. 三峽水庫(kù)水溫對(duì)浮游植物群落演替和多樣性的影響[J]. 環(huán)境工程技術(shù)學(xué)報(bào),7(2): 134-139.

    Belaoussoff S,Kevan P G,Murphy S,et al,2003. Assessing tillage disturbance on assemblages of ground beetles(Coleoptera:carabidae)by using a range of ecological indices[J].Biodiversity & Conservation,12(5):851-882.

    Field C B,Behrenfeld M J,Randerson J T,et al,1998. Primary production of the biosphere: integrating terrestrial and oceanic components[J]. Science,281(5374):237-240.

    Pitk?nen H, Lehtoranta J, R?ike A,2001. Internal Nutrient Fluxes Counteract Decreases in External Load: The Case of the Estuarial Eastern Gulf of Finland,Baltic Sea[J]. Ambio A Journal of the Human Environment,30(4):195-201.

    Lampitt R S,Wishner K F,Turley C M,et al,1993. Marine snow studies in the Northeast Atlantic Ocean: distribution,composition and role as a food source for migrating plankton[J]. Marine Biology,116(4):689-702.

    Ma C,Chula M P,Yu H,et al,2019. Spatial and temporal variation of phytoplankton functional groups in extremely alkaline Dali Nur Lake,North China[J]. Journal of Freshwater Ecology,34(1):91-105.

    Obrador B,Pretus J L,2008. Light regime and components of turbidity in a Mediterranean coastal lagoon[J]. Estuarine Coastal & Shelf Science,77(1):123-133.

    Padisák J,Crossetti L O,Naselli-Flores L,2009. Use and misuse in the application of the phytoplankton functional classification: a critical review with updates[J]. Hydrobiologia, 621(1):1-19.

    Reynolds C S,Vera H,Carla K,et al,2002. Towards a functional classification of the freshwater phytoplankton[J]. Journal of Plankton Research,(5):417-428.

    (責(zé)任編輯? ?鄭金秀)

    Payment Preference and Latent Class Identification of Residents for Ecological Compensation in the Xijiang River Basin

    CHENG Peng1, TANG Hou‐tian2, DONG Yue1, JIANG Ping1, LIU Yao‐lin1

    (1. School of Resource and Environmental Sciences, Wuhan University, Wuhan? ?430079, P.R. China;

    2. School of Public Administration, Central South University, Changsha? ?410083, P.R. China)

    Abstract: In this study, the Xijiang River basin was selected for a case study to explore differences in the payment preference of residents and factors influencing the improvement of watershed ecology. Our aim was to provide information and theoretical basis for the government to formulate and implement differentiated protection policies for the ecological environment of the watershed. Guangxi Province and Guangdong Province are the major areas of the Xijiang River basin, and Laibin and Wuzhou (Guangxi Province) and Zhaoqing and Foshan (Guangdong Province) were selected for study. The assessment framework for ecological compensation standards was based on data on the natural landscape proportion, tourism index, forest coverage, water quality, biodiversity and water quantity in the Xijiang River basin. The multinomial Logit model and latent class model were used to analyze the payment preferences of residents in the four cities and identified the latent preference classes of residents. Results show that: (1) The respondents had significantly different payment preferences for ecological attributes of the river basin, and they attached importance to water quality, forest coverage, and tourism index, but attached little importance to the proportion of natural landscape and biodiversity. (2) The overall maximum willingness of respondents to pay for improving ecosystem services in the Xijiang River basin followed the order: Wuzhou [549.35 yuan / (household · year)] > Foshan [507.60 yuan / (household · year)] > Zhaoqing [452.89 yuan / (household · year)] > Laibin [330.27 yuan / (household · year)]. The compensation surplus of Wuzhou and Foshan in the developed areas of Guangxi and Guangdong was higher than in the less developed areas of Laibin and Zhaoqing, indicating a strong correlation between the level of regional economic development and resident willingness to pay for ecological attributes. (3) The estimation results of the latent class model confirm that the residents of the four cities had significant preference heterogeneity. Laibin, Zhaoqing, and Foshan were divided into three latent classes, respectively, and Wuzhou was divided into two latent classes. In the future, the government should develop a range of ecological compensation based on the actual situation in different areas, implement a differentiated watershed ecological compensation policy, and use a variety of measures to motivate residents to participate in ecological and environmental protection.

    Key words:ecosystem service; watershed ecological compensation; preference heterogeneity; choice experiment method; latent class model; Xijiang River basin

    收稿日期:2021-08-20? ? ? 修回日期: 2023-05-23

    基金項(xiàng)目:農(nóng)業(yè)財(cái)政專(zhuān)項(xiàng)“長(zhǎng)江漁業(yè)資源與環(huán)境調(diào)查”;財(cái)政部和農(nóng)業(yè)農(nóng)村部:現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系專(zhuān)項(xiàng)資金(CARS-46);中國(guó)水產(chǎn)科學(xué)研究院院級(jí)創(chuàng)新團(tuán)隊(duì)項(xiàng)目(2020TD57)。

    作者簡(jiǎn)介:龔進(jìn)玲,1990年生,女,碩士,主要從事水生生物學(xué)研究。E-mail:jlgong123@163.com

    通信簡(jiǎn)介:楊德國(guó),1964年生,男,研究員。E-mail:yangdg@yfi.ac.cn

    猜你喜歡
    浮游植物環(huán)境因子金沙江
    金沙江(女聲獨(dú)唱)
    開(kāi)掛的“金沙江”
    行走金沙江
    塔里木河阿拉爾段及多浪水庫(kù)浮游生物的調(diào)查研究
    阿克蘇地區(qū)杏李與歐洲李氣體交換特征對(duì)比研究
    桉樹(shù)人工林胸徑變化特征與環(huán)境因子的關(guān)系研究
    渭河陜西段浮游植物群落結(jié)構(gòu)及水質(zhì)評(píng)價(jià)
    濟(jì)南流域春季浮游植物多樣性的研究
    秦皇島海域夜光藻種群密度與環(huán)境因子的關(guān)系
    峽江樞紐庫(kù)區(qū)水生生物調(diào)查及水質(zhì)評(píng)價(jià)
    亚洲中文av在线| 在线观看人妻少妇| 免费高清在线观看视频在线观看| 新久久久久国产一级毛片| 欧美人与性动交α欧美精品济南到| 人妻 亚洲 视频| 又紧又爽又黄一区二区| 91老司机精品| 国产精品亚洲av一区麻豆| 国产精品成人在线| 国产欧美日韩一区二区三 | 人人妻人人添人人爽欧美一区卜| 国产免费福利视频在线观看| 国产精品九九99| 亚洲人成77777在线视频| 无遮挡黄片免费观看| av在线app专区| 国产在线免费精品| 高清黄色对白视频在线免费看| 永久免费av网站大全| 视频在线观看一区二区三区| 国产熟女午夜一区二区三区| 午夜久久久在线观看| 中文字幕人妻熟女乱码| 亚洲精品日韩在线中文字幕| 亚洲国产精品一区二区三区在线| 七月丁香在线播放| 黄片播放在线免费| 手机成人av网站| 建设人人有责人人尽责人人享有的| 国产真人三级小视频在线观看| 国产av一区二区精品久久| www.av在线官网国产| 国产福利在线免费观看视频| 丝袜在线中文字幕| 成人午夜精彩视频在线观看| 丝袜美腿诱惑在线| 亚洲五月婷婷丁香| 中文字幕色久视频| 校园人妻丝袜中文字幕| 一二三四在线观看免费中文在| 国产亚洲av高清不卡| 日韩 欧美 亚洲 中文字幕| 成年人免费黄色播放视频| 欧美激情高清一区二区三区| 69精品国产乱码久久久| 青春草亚洲视频在线观看| 脱女人内裤的视频| 悠悠久久av| 天天影视国产精品| 成人亚洲精品一区在线观看| 久久免费观看电影| 日韩制服骚丝袜av| 操出白浆在线播放| 可以免费在线观看a视频的电影网站| 伦理电影免费视频| 亚洲,一卡二卡三卡| 久久久久网色| 国产无遮挡羞羞视频在线观看| 赤兔流量卡办理| 看十八女毛片水多多多| 日韩大码丰满熟妇| 国产亚洲午夜精品一区二区久久| 中文字幕高清在线视频| 在线观看人妻少妇| 日本91视频免费播放| 亚洲国产精品国产精品| 亚洲av美国av| 一区二区三区精品91| 国产97色在线日韩免费| 夜夜骑夜夜射夜夜干| 美女国产高潮福利片在线看| 亚洲av美国av| 国产一区二区在线观看av| 久久影院123| 一边摸一边做爽爽视频免费| 国产亚洲av高清不卡| 丝袜脚勾引网站| 男人舔女人的私密视频| 日本黄色日本黄色录像| 亚洲伊人久久精品综合| 青春草亚洲视频在线观看| 国产精品久久久久成人av| 亚洲三区欧美一区| 另类亚洲欧美激情| 欧美xxⅹ黑人| 久久久久精品国产欧美久久久 | 美女主播在线视频| 黄色a级毛片大全视频| 在线 av 中文字幕| 精品亚洲成a人片在线观看| 日韩 欧美 亚洲 中文字幕| av有码第一页| 男人操女人黄网站| 在线观看国产h片| 夜夜骑夜夜射夜夜干| 91九色精品人成在线观看| 亚洲少妇的诱惑av| 亚洲人成77777在线视频| 久久免费观看电影| 亚洲国产中文字幕在线视频| 亚洲av电影在线观看一区二区三区| 欧美精品啪啪一区二区三区 | 在线观看免费日韩欧美大片| 亚洲精品久久午夜乱码| 久久久久国产精品人妻一区二区| 在线观看一区二区三区激情| 国产精品.久久久| 电影成人av| 99香蕉大伊视频| 欧美精品av麻豆av| 2021少妇久久久久久久久久久| 亚洲国产欧美日韩在线播放| 亚洲人成网站在线观看播放| 性色av乱码一区二区三区2| 亚洲av男天堂| 一级毛片女人18水好多 | 交换朋友夫妻互换小说| 国产精品亚洲av一区麻豆| 搡老乐熟女国产| 丰满人妻熟妇乱又伦精品不卡| 精品国产乱码久久久久久小说| 久久性视频一级片| 91国产中文字幕| 免费人妻精品一区二区三区视频| 视频在线观看一区二区三区| 国产精品 国内视频| 大陆偷拍与自拍| 桃花免费在线播放| 欧美激情极品国产一区二区三区| 国产三级黄色录像| 中文字幕人妻丝袜一区二区| 热99久久久久精品小说推荐| 欧美日韩亚洲综合一区二区三区_| 欧美乱码精品一区二区三区| 免费在线观看黄色视频的| 一级片免费观看大全| 成在线人永久免费视频| 黄色一级大片看看| 久久毛片免费看一区二区三区| 精品人妻1区二区| 国产一级毛片在线| 日韩制服丝袜自拍偷拍| www.自偷自拍.com| 尾随美女入室| 一区二区三区精品91| 国产免费视频播放在线视频| 欧美日韩黄片免| 精品一区二区三卡| 成人免费观看视频高清| 如日韩欧美国产精品一区二区三区| 国产成人av教育| 一级片'在线观看视频| 黄色怎么调成土黄色| 欧美日韩av久久| a级片在线免费高清观看视频| 日本a在线网址| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久精品精品| 一边摸一边做爽爽视频免费| 成人午夜精彩视频在线观看| 亚洲熟女精品中文字幕| 亚洲欧美日韩高清在线视频 | 精品人妻熟女毛片av久久网站| 你懂的网址亚洲精品在线观看| √禁漫天堂资源中文www| 超碰成人久久| 在线观看人妻少妇| 天天添夜夜摸| 99热网站在线观看| 丰满饥渴人妻一区二区三| 一区福利在线观看| 国产精品 国内视频| 免费在线观看黄色视频的| 一本久久精品| 亚洲伊人色综图| 九草在线视频观看| 一区二区av电影网| 美国免费a级毛片| 午夜福利视频在线观看免费| 99精国产麻豆久久婷婷| 国产成人91sexporn| 国产1区2区3区精品| 蜜桃国产av成人99| 国产免费现黄频在线看| 日韩伦理黄色片| 一区二区三区乱码不卡18| 亚洲 国产 在线| 捣出白浆h1v1| 国产精品秋霞免费鲁丝片| 日日爽夜夜爽网站| 在线观看一区二区三区激情| 欧美在线一区亚洲| 黑人猛操日本美女一级片| 丝袜美足系列| 久久精品久久精品一区二区三区| 午夜免费成人在线视频| 在线精品无人区一区二区三| 欧美日韩福利视频一区二区| 亚洲av片天天在线观看| 国产亚洲精品第一综合不卡| 汤姆久久久久久久影院中文字幕| 久久影院123| 真人做人爱边吃奶动态| 亚洲欧洲日产国产| 19禁男女啪啪无遮挡网站| 乱人伦中国视频| 男人添女人高潮全过程视频| 91麻豆精品激情在线观看国产 | 亚洲欧洲国产日韩| 高清av免费在线| 久久ye,这里只有精品| 丰满人妻熟妇乱又伦精品不卡| 久久人人爽av亚洲精品天堂| 亚洲欧美一区二区三区久久| 国产精品一区二区精品视频观看| 一级a爱视频在线免费观看| 女人高潮潮喷娇喘18禁视频| 亚洲中文av在线| 精品国产一区二区久久| 人妻人人澡人人爽人人| 亚洲成av片中文字幕在线观看| 国产一区二区三区av在线| 欧美日韩视频高清一区二区三区二| 久久久精品区二区三区| 国产在线免费精品| 欧美人与性动交α欧美精品济南到| 国产1区2区3区精品| 日韩制服丝袜自拍偷拍| 国产成人精品无人区| 国产人伦9x9x在线观看| 亚洲精品国产av蜜桃| 国产一卡二卡三卡精品| 国产在线免费精品| 在线av久久热| 三上悠亚av全集在线观看| 欧美激情极品国产一区二区三区| av片东京热男人的天堂| 欧美人与性动交α欧美软件| videosex国产| 成人三级做爰电影| 久久精品久久久久久久性| 80岁老熟妇乱子伦牲交| av网站免费在线观看视频| 精品一区在线观看国产| 老司机靠b影院| 性高湖久久久久久久久免费观看| 视频区图区小说| 亚洲五月色婷婷综合| 男女下面插进去视频免费观看| 亚洲综合色网址| 少妇粗大呻吟视频| 男人添女人高潮全过程视频| 国产成人av教育| 一本综合久久免费| 欧美精品高潮呻吟av久久| 后天国语完整版免费观看| 热re99久久国产66热| 91成人精品电影| 成在线人永久免费视频| 七月丁香在线播放| 免费在线观看影片大全网站 | 亚洲美女黄色视频免费看| 免费av中文字幕在线| 激情视频va一区二区三区| 久久精品国产亚洲av涩爱| 日韩伦理黄色片| 婷婷色综合www| 日本欧美视频一区| 欧美在线一区亚洲| 咕卡用的链子| 国产av一区二区精品久久| 欧美精品一区二区免费开放| 九色亚洲精品在线播放| 另类亚洲欧美激情| 在线观看免费日韩欧美大片| 99久久99久久久精品蜜桃| 久久精品国产a三级三级三级| 久久精品人人爽人人爽视色| 两人在一起打扑克的视频| 2021少妇久久久久久久久久久| 久久久精品免费免费高清| 免费看十八禁软件| 婷婷色综合www| 人妻 亚洲 视频| 后天国语完整版免费观看| 成人国产av品久久久| 欧美老熟妇乱子伦牲交| 日韩伦理黄色片| 又黄又粗又硬又大视频| 日本欧美视频一区| 亚洲av男天堂| 黄网站色视频无遮挡免费观看| 国产一区二区激情短视频 | 丝袜美腿诱惑在线| 91精品伊人久久大香线蕉| 在线观看一区二区三区激情| 亚洲av国产av综合av卡| 精品高清国产在线一区| 啦啦啦 在线观看视频| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 一本一本久久a久久精品综合妖精| 麻豆国产av国片精品| 亚洲av电影在线观看一区二区三区| 在线观看www视频免费| 国产淫语在线视频| 亚洲一区二区三区欧美精品| 国产免费一区二区三区四区乱码| 日本欧美视频一区| 精品欧美一区二区三区在线| 亚洲国产欧美日韩在线播放| 永久免费av网站大全| 免费av中文字幕在线| 国产爽快片一区二区三区| 黄频高清免费视频| 校园人妻丝袜中文字幕| 久久久久精品人妻al黑| 在线av久久热| 久久久久久免费高清国产稀缺| 国产日韩欧美在线精品| 满18在线观看网站| av国产久精品久网站免费入址| 精品一区二区三卡| 欧美日韩视频精品一区| 人妻人人澡人人爽人人| 成人亚洲欧美一区二区av| 啦啦啦在线观看免费高清www| 999精品在线视频| 精品国产一区二区三区久久久樱花| 午夜激情久久久久久久| 国产亚洲一区二区精品| 大香蕉久久网| 久久ye,这里只有精品| 丝瓜视频免费看黄片| 美女大奶头黄色视频| 国产一区二区 视频在线| 亚洲国产毛片av蜜桃av| 九草在线视频观看| 一级毛片我不卡| 啦啦啦视频在线资源免费观看| 又紧又爽又黄一区二区| 国产一区二区激情短视频 | 国产精品秋霞免费鲁丝片| 曰老女人黄片| 夫妻性生交免费视频一级片| 亚洲av电影在线观看一区二区三区| 在线观看www视频免费| 精品少妇黑人巨大在线播放| 亚洲欧洲精品一区二区精品久久久| 国产高清视频在线播放一区 | 十八禁网站网址无遮挡| 中文精品一卡2卡3卡4更新| 日韩一区二区三区影片| 视频区欧美日本亚洲| 精品久久蜜臀av无| 飞空精品影院首页| 国产精品 国内视频| 在线观看免费高清a一片| 美女中出高潮动态图| 欧美激情极品国产一区二区三区| 又粗又硬又长又爽又黄的视频| 久久人人爽人人片av| 国产欧美亚洲国产| 各种免费的搞黄视频| 啦啦啦啦在线视频资源| 中文字幕人妻丝袜制服| 黄片播放在线免费| 男人舔女人的私密视频| 最近中文字幕2019免费版| 欧美日韩福利视频一区二区| 91精品伊人久久大香线蕉| 久久久久久久大尺度免费视频| 黄片小视频在线播放| videos熟女内射| 999久久久国产精品视频| 亚洲欧美精品自产自拍| 亚洲国产精品一区二区三区在线| 一级毛片电影观看| 三上悠亚av全集在线观看| 最近最新中文字幕大全免费视频 | 丰满饥渴人妻一区二区三| 叶爱在线成人免费视频播放| 爱豆传媒免费全集在线观看| 欧美成人午夜精品| 色网站视频免费| av天堂久久9| 国产一区有黄有色的免费视频| 两性夫妻黄色片| 色婷婷久久久亚洲欧美| 午夜影院在线不卡| 性色av乱码一区二区三区2| 国产一区亚洲一区在线观看| 国产精品久久久久成人av| 国产日韩欧美在线精品| 人成视频在线观看免费观看| 日本猛色少妇xxxxx猛交久久| a 毛片基地| 国产成人系列免费观看| 午夜福利免费观看在线| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色a级毛片大全视频| 9热在线视频观看99| 久久久久久久国产电影| 日韩大码丰满熟妇| 又黄又粗又硬又大视频| 最近最新中文字幕大全免费视频 | 黑人巨大精品欧美一区二区蜜桃| 最新的欧美精品一区二区| 亚洲精品av麻豆狂野| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 一个人免费看片子| 99香蕉大伊视频| 99国产精品一区二区三区| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 亚洲,欧美精品.| 操美女的视频在线观看| 亚洲精品第二区| 国产色视频综合| 亚洲午夜精品一区,二区,三区| videosex国产| 国产片内射在线| 一本一本久久a久久精品综合妖精| 男女边吃奶边做爰视频| 搡老岳熟女国产| 亚洲美女黄色视频免费看| 亚洲国产欧美网| 日韩熟女老妇一区二区性免费视频| 两个人看的免费小视频| 一区福利在线观看| 爱豆传媒免费全集在线观看| 悠悠久久av| 国产麻豆69| 精品国产超薄肉色丝袜足j| 色精品久久人妻99蜜桃| 人人妻人人添人人爽欧美一区卜| 丝袜脚勾引网站| 中文字幕另类日韩欧美亚洲嫩草| 色视频在线一区二区三区| 精品人妻在线不人妻| 国产欧美日韩综合在线一区二区| 亚洲自偷自拍图片 自拍| 欧美黄色片欧美黄色片| 国产精品秋霞免费鲁丝片| 欧美性长视频在线观看| 国产成人一区二区在线| 日韩av免费高清视频| 99久久人妻综合| 欧美日韩视频精品一区| 日韩 亚洲 欧美在线| 日本91视频免费播放| 国产成人影院久久av| 2021少妇久久久久久久久久久| 日韩熟女老妇一区二区性免费视频| 国产精品欧美亚洲77777| tube8黄色片| 亚洲人成电影免费在线| 久久精品久久久久久久性| 亚洲国产欧美一区二区综合| 老司机深夜福利视频在线观看 | 日韩精品免费视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色av中文字幕| 电影成人av| 欧美国产精品一级二级三级| 晚上一个人看的免费电影| 只有这里有精品99| 亚洲精品一卡2卡三卡4卡5卡 | 久久久久久免费高清国产稀缺| 午夜av观看不卡| 国产日韩欧美视频二区| 视频区欧美日本亚洲| 精品久久蜜臀av无| 国产亚洲午夜精品一区二区久久| 女人久久www免费人成看片| e午夜精品久久久久久久| 丰满迷人的少妇在线观看| 精品国产超薄肉色丝袜足j| 99精国产麻豆久久婷婷| 国产成人欧美在线观看 | 天堂中文最新版在线下载| 亚洲欧洲日产国产| 免费在线观看视频国产中文字幕亚洲 | 制服人妻中文乱码| 欧美日韩亚洲高清精品| 久久狼人影院| 久久性视频一级片| 亚洲精品国产一区二区精华液| 五月天丁香电影| 国产欧美日韩精品亚洲av| 久久99精品国语久久久| 午夜免费鲁丝| 黑人巨大精品欧美一区二区蜜桃| 香蕉丝袜av| 热99久久久久精品小说推荐| 久久久久久免费高清国产稀缺| 黄片小视频在线播放| 国产深夜福利视频在线观看| 国产精品av久久久久免费| 两性夫妻黄色片| av电影中文网址| 欧美97在线视频| 成年人黄色毛片网站| 欧美人与性动交α欧美软件| 亚洲精品一二三| 欧美人与性动交α欧美精品济南到| 女人爽到高潮嗷嗷叫在线视频| 午夜福利,免费看| 国产精品一区二区在线不卡| 狠狠精品人妻久久久久久综合| 男女下面插进去视频免费观看| 久久精品久久久久久久性| 日韩大码丰满熟妇| 99热全是精品| 亚洲av欧美aⅴ国产| 欧美日韩黄片免| 亚洲人成电影观看| 久久久久国产精品人妻一区二区| 男女下面插进去视频免费观看| 精品久久久精品久久久| 久久av网站| 色播在线永久视频| 波野结衣二区三区在线| 9191精品国产免费久久| 国产欧美日韩精品亚洲av| 在线观看国产h片| 看免费成人av毛片| 久久久久久久精品精品| 精品少妇黑人巨大在线播放| 色视频在线一区二区三区| 熟女av电影| 亚洲精品av麻豆狂野| 国产免费现黄频在线看| 久久精品国产综合久久久| videosex国产| 少妇猛男粗大的猛烈进出视频| 久久久精品免费免费高清| 中文精品一卡2卡3卡4更新| 十八禁高潮呻吟视频| 国产色视频综合| 久久久久国产精品人妻一区二区| 波野结衣二区三区在线| 在线观看国产h片| av在线老鸭窝| 少妇人妻久久综合中文| av视频免费观看在线观看| 亚洲综合色网址| 国产真人三级小视频在线观看| 国产成人av激情在线播放| 99国产综合亚洲精品| 久久人人爽人人片av| e午夜精品久久久久久久| 欧美老熟妇乱子伦牲交| 2021少妇久久久久久久久久久| 午夜精品国产一区二区电影| 一级片'在线观看视频| 十八禁网站网址无遮挡| 中文乱码字字幕精品一区二区三区| 久9热在线精品视频| 少妇精品久久久久久久| 久久久久久久久免费视频了| 51午夜福利影视在线观看| 亚洲欧美一区二区三区国产| 制服诱惑二区| 国产在线观看jvid| 色视频在线一区二区三区| 亚洲av日韩精品久久久久久密 | 美国免费a级毛片| 久久99精品国语久久久| 天天影视国产精品| 啦啦啦视频在线资源免费观看| 一级片免费观看大全| 亚洲精品一区蜜桃| 久久久久久久久免费视频了| 亚洲欧美一区二区三区黑人| 老汉色∧v一级毛片| xxx大片免费视频| 亚洲,欧美,日韩| 男人爽女人下面视频在线观看| 女警被强在线播放| kizo精华| 亚洲精品国产av成人精品| 青草久久国产| 在线观看人妻少妇| 啦啦啦中文免费视频观看日本| 日本a在线网址| 黄色 视频免费看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美一区二区三区黑人| 中文字幕亚洲精品专区| 欧美中文综合在线视频| 国产成人精品在线电影| 啦啦啦啦在线视频资源| 高清av免费在线| av有码第一页| 国产日韩欧美视频二区| xxx大片免费视频| 久久久国产精品麻豆| 99热全是精品| 自拍欧美九色日韩亚洲蝌蚪91| videosex国产| 在线观看www视频免费| 亚洲精品国产一区二区精华液| 免费高清在线观看日韩| 亚洲精品国产区一区二| 夫妻午夜视频| 国产有黄有色有爽视频| 国产一级毛片在线| 美女扒开内裤让男人捅视频| 黑丝袜美女国产一区| 黄色 视频免费看| 精品一区在线观看国产| 精品免费久久久久久久清纯 |