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

    渾太河春季不同水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系

    2018-03-07 08:49:56孟云飛李晨張吉吳丹楊璐殷旭旺徐宗學(xué)張遠(yuǎn)
    關(guān)鍵詞:功能區(qū)流域動(dòng)物

    孟云飛,李晨,張吉,吳丹,楊璐,殷旭旺,徐宗學(xué),張遠(yuǎn)

    (1.大連海洋大學(xué)水產(chǎn)與生命學(xué)院,遼寧省水生生物學(xué)重點(diǎn)實(shí)驗(yàn)室,遼寧大連116023;2.北京師范大學(xué)水科學(xué)研究院,北京100875;3.中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評(píng)估國家重點(diǎn)實(shí)驗(yàn)室,流域水生態(tài)保護(hù)技術(shù)研究室,北京100012)

    水生態(tài)功能區(qū)是通過對(duì)環(huán)境要素 (植被、土壤、地理、氣候),水生生境類型,水生生物分布格局等因素的分析[1],根據(jù)生態(tài)系統(tǒng)不同尺度上的空間特征與維持生態(tài)系統(tǒng)完整性的特點(diǎn),將陸地和水體相似的區(qū)域劃分而成的地理單元[1-3]。水生態(tài)功能區(qū)是基于流域水生態(tài)系統(tǒng),為了統(tǒng)一管理水生態(tài)系統(tǒng)而提出的全新概念[4-8],目的是保護(hù)水生態(tài)環(huán)境的完整性、社會(huì)經(jīng)濟(jì)發(fā)展的可持續(xù)性,并根據(jù)水生態(tài)功能區(qū)的概念和特點(diǎn),以水生態(tài)功能區(qū)為基本單元,對(duì)生態(tài)系統(tǒng)進(jìn)行健康評(píng)價(jià)和不同水生態(tài)環(huán)境的保護(hù)和修復(fù)。

    大型底棲動(dòng)物 (macroinvertebrate)如環(huán)節(jié)動(dòng)物、線形動(dòng)物、軟體動(dòng)物和甲殼動(dòng)物[9],是河流生態(tài)系統(tǒng)的重要組成部分,可推動(dòng)生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量流動(dòng)[10-11],對(duì)生態(tài)系統(tǒng)的構(gòu)成具有重要的生物學(xué)意義。大型底棲動(dòng)物不僅具有較強(qiáng)的區(qū)域性,而且對(duì)生態(tài)系統(tǒng)的變化比較敏感,不同種類的大型底棲動(dòng)物對(duì)環(huán)境的適應(yīng)以及對(duì)環(huán)境變化的耐受程度和敏感程度不同,所以大型底棲動(dòng)物對(duì)環(huán)境狀況具有指示作用[12-13]。研究不同水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)與環(huán)境因子的關(guān)系,有利于深入了解河流生態(tài)系統(tǒng)中不同水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)特點(diǎn)。

    水生態(tài)功能區(qū)可以較好地反映流域水生態(tài)系統(tǒng)特征。研究表明,太子河不同水生態(tài)功能區(qū)內(nèi)的土地利用類型、人為干擾程度、氣候特征等均存在差異,從而導(dǎo)致水環(huán)境因子,以及水生生物與環(huán)境因子的響應(yīng)關(guān)系均存在明顯差異[2]。本研究中,以渾太河流域?yàn)檠芯繀^(qū)域,基于渾太河流域不同水生態(tài)功能分區(qū)的不同水生態(tài)特征,采用生態(tài)學(xué)方法,對(duì)不同水生態(tài)功能區(qū)大型底棲動(dòng)物的群落結(jié)構(gòu)、環(huán)境因子的比較和對(duì)應(yīng)區(qū)域大型底棲動(dòng)物群落結(jié)構(gòu)與環(huán)境因子的相關(guān)性進(jìn)行了分析,旨在為渾太河水生生態(tài)功能區(qū)健康評(píng)估提供參考。

    1 材料與方法

    1.1 材料

    1.1.1 研究區(qū)域概況 渾太河流域 (40.45°~42.30°N, 122.00°~125.30°E) 位于中國東北部,屬于濕潤半濕潤暖溫帶季風(fēng)氣候,流域面積為2.73×104km2, 集水面積為 1.15×104km2, 年平均溫度為9℃,年平均降水量為686.4 mm。渾太河起源于長白山,最終流入渤海,由渾河和太子河兩條河流組成,是遼河最大的支流。土地覆蓋主要由森林組成,由于經(jīng)濟(jì)的發(fā)展,森林面積目前正越來越多地轉(zhuǎn)化為工農(nóng)業(yè)用地,淡水生態(tài)系統(tǒng)的惡化越來越嚴(yán)重,因此,渾太河流域生態(tài)狀況受到越來越多的關(guān)注。

    張楠等[1]的研究中確定渾太河分為3個(gè)水生態(tài)功能分區(qū),分別為平原農(nóng)業(yè)水生態(tài)功能Ⅰ區(qū)、丘陵森林水生態(tài)功能Ⅱ區(qū)和山地森林水生態(tài)功能Ⅲ區(qū)。水生態(tài)功能Ⅰ區(qū)平均海拔為29 m,地貌類型主要是平原,年降水量較低,由于農(nóng)田、工業(yè)和城市分布,該區(qū)植被覆蓋率低且受人為干擾比較大。水生態(tài)功能Ⅱ區(qū)是平原向山地過度區(qū),平均海拔為245 m,地貌類型主要是丘陵,年降水量高于水生態(tài)功能Ⅰ區(qū),人類活動(dòng)相對(duì)較少。水生態(tài)功能Ⅲ區(qū)平均海拔為349 m,地貌類型主要是山地,年降水量較高,森林覆蓋率高且受人為干擾少,生境多樣化較高。

    1.1.2 采樣時(shí)間與站點(diǎn)設(shè)置 2012年春季在渾太河干流和支流共設(shè)置65個(gè)采樣點(diǎn) (圖1),其中在渾河設(shè)置22個(gè)采樣點(diǎn),太子河設(shè)置43個(gè)采樣點(diǎn),即水生態(tài)功能Ⅰ區(qū)設(shè)置30個(gè)采樣點(diǎn) (H11~H16,H18~H22, H25, T24~T41), 水生態(tài)功能Ⅱ區(qū)設(shè)置16個(gè)采樣點(diǎn) (H1~H10, T1~T4, T42~T43),水生態(tài)功能Ⅲ區(qū)設(shè)置19個(gè)采樣點(diǎn) (T5~T23)。

    1.2 方法

    1.2.1 底棲樣品采集、分類與鑒定 在所選定采樣點(diǎn)的100 m范圍之內(nèi),使用索伯網(wǎng) (surber,網(wǎng)口為30 cm×30 cm,網(wǎng)孔徑為500 μm)隨機(jī)采集2個(gè)平行樣本。水面較大的樣點(diǎn)于近岸處取樣。將網(wǎng)內(nèi)所有采集物經(jīng)過60目網(wǎng)篩轉(zhuǎn)入到30 mL塑料瓶中,加入體積分?jǐn)?shù)為90%的酒精溶液進(jìn)行保存。在實(shí)驗(yàn)室采用人工挑揀的方法將底棲動(dòng)物樣品轉(zhuǎn)入到100 mL的塑料瓶中,同時(shí)加入體積分?jǐn)?shù)為95%的酒精溶液保存待檢。在顯微鏡或解剖鏡下對(duì)底棲動(dòng)物進(jìn)行分類、計(jì)數(shù),利用解剖鏡和顯微鏡進(jìn)行分類鑒定,樣品盡量鑒定到屬或者種。

    圖1 渾太河流域采樣點(diǎn)位圖Fig.1 Location sketch of sample sites in the Huntai River

    1.2.2 環(huán)境因子的測定 采用MAGELLAN全球定位系統(tǒng) (eXplorist-200)現(xiàn)場測定海拔和經(jīng)緯度;采用YSI-85型便攜式水質(zhì)分析儀現(xiàn)場測定水體的pH、水溫 (WT)、電導(dǎo)率 (COND)、溶解氧(DO)、總?cè)芙夤腆w (TDS);使用流速儀(FP111)現(xiàn)場測定水深 (WD)、流速 (CV);采用皮尺測量河寬 (CW);根據(jù)水深、流速和河寬計(jì)算流量 (FLO);參照 《國家地面水環(huán)境質(zhì)量標(biāo)準(zhǔn) (GB 3838-2002)》在實(shí)驗(yàn)室測定總氮 (TN)、總磷 (TP)、懸浮物 (SES)、高錳酸鉀鹽指數(shù)(CODMn)、鹽度 (SAL)、硬度 (HAR)等水質(zhì)指標(biāo)。

    1.3 數(shù)據(jù)處理

    采用Pcord 5.0軟件對(duì)3個(gè)水生態(tài)功能區(qū)的指示種進(jìn)行分析;使用 Biodiversity pro軟件計(jì)算Shannon-Wiener多樣性指數(shù) (H′)、Pielou均勻度指數(shù) (J)和物種豐富度指數(shù) (D);采用SPSS 20.0軟件進(jìn)行單因素方差分析 (One-way ANOVA),分析比較3個(gè)水生態(tài)功能區(qū)各指數(shù)間的差異性和環(huán)境因子的差異性,并進(jìn)行Tukey-Kramer檢驗(yàn)。采用Origin 8.5軟件對(duì)存在顯著差異的環(huán)境因子進(jìn)行箱體圖比較分析,以便能更直觀地看到不同水生態(tài)功能區(qū)環(huán)境因子的差異。使用CANOCO 4.5軟件對(duì)各樣點(diǎn)大型底棲動(dòng)物密度數(shù)據(jù)進(jìn)行趨勢對(duì)應(yīng)分析 (Detrended correspondence analysis, DCA),當(dāng)SD(梯度)值>2時(shí),選擇典范對(duì)應(yīng)分析 (Canonical correspondence analysis, CCA), 分析大型底棲動(dòng)物群落結(jié)構(gòu)的主要驅(qū)動(dòng)因子,否則采用冗余度分析 (Redundancy analysis,RDA)。利用 SPSS 20.0軟件對(duì)3個(gè)水生態(tài)功能區(qū)大型底棲動(dòng)物Shannon-Wiener指數(shù)、均勻度指數(shù)和豐富度指數(shù)與環(huán)境因子進(jìn)行Pearson相關(guān)系數(shù)分析。

    在進(jìn)行CCA分析時(shí),對(duì)大型底棲動(dòng)物數(shù)據(jù)和環(huán)境因子數(shù)據(jù) (除pH外)進(jìn)行 [lg(x+1)]轉(zhuǎn)換;并對(duì)數(shù)據(jù)進(jìn)行999次蒙特卡洛循環(huán)檢驗(yàn)(Monte Carlo test,MCT),以判定顯著影響大型底棲動(dòng)物群落空間分布特征的環(huán)境因子。

    2 結(jié)果與分析

    2.1 不同水生態(tài)功能區(qū)大型底棲動(dòng)物的群落特征

    2012年春季在渾太河流域65個(gè)采樣點(diǎn)共采集到大型底棲動(dòng)物71屬 (種),隸屬7綱13目35科。其中,優(yōu)勢種有三帶環(huán)足搖蚊Cricotopus trifasciatus、雙線環(huán)足搖蚊 Cricotopus bicinctus、三輪環(huán)足搖蚊Cricotopus triannulatus、霍甫水絲蚓Limnodrilus hoffmeisteri、克拉伯水絲蚓Limnodrilus claparedianus;水生昆蟲51種 (占總種數(shù)的71.83%),軟體動(dòng)物10種 (占總種數(shù)的14.08%),環(huán)節(jié)動(dòng)物7種 (占總種數(shù)的9.86%),甲殼動(dòng)物3種 (占總種數(shù)的4.23%)(表1)。3個(gè)水生態(tài)功能區(qū)大型底棲動(dòng)物的分布情況顯示:水生態(tài)功能Ⅰ區(qū)采集到的大型底棲動(dòng)物42屬 (種),隸屬7綱13目21科,其中水生昆蟲25種 (占59.52%),軟體動(dòng)物8種(占19.05%),環(huán)節(jié)動(dòng)物6種 (占14.29%),甲殼動(dòng)物3種 (占7.14%);水生態(tài)功能Ⅱ區(qū)采集到的大型底棲動(dòng)物35屬 (種),隸屬5綱9目18科,其中水生昆蟲27種 (占77.14%),軟體動(dòng)物4種(占11.43%),環(huán)節(jié)動(dòng)物4種 (占11.43%),未發(fā)現(xiàn)甲殼動(dòng)物;水生態(tài)功能Ⅲ區(qū)采集到的大型底棲動(dòng)物47屬 (種),隸屬5綱9目26科,其中水生昆蟲37種 (占 78.72%),軟體動(dòng)物 4種 (占8.51%),環(huán)節(jié)動(dòng)物5種 (占13.51%),甲殼動(dòng)物1種 (占2.70%)(表1)。比較分析可知,水生態(tài)功能Ⅰ區(qū)大型底棲動(dòng)物主要由水生昆蟲、軟體動(dòng)物、環(huán)節(jié)動(dòng)物組成,甲殼動(dòng)物較少;水生態(tài)功能Ⅱ區(qū)大型底棲動(dòng)物主要由水生昆蟲組成,軟體動(dòng)物和環(huán)節(jié)動(dòng)物相對(duì)較少;水生態(tài)功能Ⅲ區(qū)大型底棲動(dòng)物主要由水生昆蟲組成,軟體動(dòng)物和環(huán)節(jié)動(dòng)物相對(duì)較少,甲殼動(dòng)物更少。

    表1 渾太河流域3個(gè)水生態(tài)功能區(qū)大型底棲動(dòng)物的群落組成Tab.1 Community composition of macrobenthic fauna in three freshwater ecosystem function regions in Huntai River Basin

    對(duì)渾太河流域3個(gè)水生態(tài)功能區(qū)指示種的分析結(jié)果如表2所示。結(jié)果表明:水生態(tài)功能Ⅰ區(qū)有2個(gè)指示種 (屬),分別為梯形多足搖蚊Polypedilum scalaenum和中華新米蝦Neocaridina denticulata sinensis;水生態(tài)功能Ⅱ區(qū)有2個(gè)指示種 (屬),分別為韋氏四節(jié)蜉Acentrella sp.和三帶環(huán)足搖蚊Cricotopus trifasciatus;水生態(tài)功能Ⅲ區(qū)有4個(gè)指示種(屬),分別為韋氏四節(jié)蜉Heterocleon amplum、動(dòng)蜉屬 Cinygma sp.、黑翅蜉 Ephemera nigropter、稻田蛆Dicranota sp.。

    從表3可見:Shannon-Wiener多樣性指數(shù)和物種豐富度在3個(gè)不同的水生態(tài)功能區(qū)間有極顯著性差異 (P<0.01),水生態(tài)功能Ⅲ區(qū)Shannon-Wiener多樣性指數(shù) (1.45)和物種豐富度 (5.58)最高,水生態(tài)功能Ⅱ區(qū)的Shannon-Wiener多樣性指數(shù)最低,為0.74,水生態(tài)功能Ⅰ區(qū)的物種豐富度最低,為2.60;3個(gè)水生態(tài)功能區(qū)的Pielou均勻度指數(shù)有顯著性差異 (P<0.05),Pielou均勻度指數(shù)在水生態(tài)功能Ⅰ區(qū)最高,在水生態(tài)功能Ⅱ區(qū)最低,分別為0.64和0.38。這表明,不同的水生態(tài)功能區(qū)大型底棲動(dòng)物指數(shù)有顯著性差異 (P<0.05)。

    2.2 不同水生態(tài)功能區(qū)環(huán)境因子的比較

    不同水生態(tài)功能區(qū)環(huán)境因子的比較見表4。環(huán)境因子單因素方差分析表明:pH、懸浮物、草地面積在3個(gè)水生態(tài)功能區(qū)間無顯著性差異 (P>0.05);水溫、河寬、流速、流量、總磷、森林面積、水田面積在3個(gè)水生態(tài)功能區(qū)間有顯著性差異(P<0.05);海拔、電導(dǎo)率、鹽度、溶解氧、總?cè)芙夤腆w、水深、總氮、高錳酸鉀鹽指數(shù)、硬度、緩沖區(qū)流域總面積、建設(shè)用地面積、旱田面積在3個(gè)水生態(tài)功能區(qū)間有極顯著性差異 (P<0.01);其中水生態(tài)功能Ⅰ區(qū)電導(dǎo)率、鹽度、總?cè)芙夤腆w、水深、總氮、硬度、緩沖區(qū)流域總面積、建設(shè)用地面積、旱田面積的平均值高于水生態(tài)功能Ⅱ、Ⅲ區(qū),水生態(tài)功能Ⅱ區(qū)溶解氧的平均值最高,水生態(tài)功能Ⅲ區(qū)海拔的平均值最高。

    表2 渾太河流域3個(gè)水生態(tài)功能區(qū)的指示種Tab.2 Indicator species in three freshwater ecosystem function regions in Huntai River Basin

    表3 渾太河流域3個(gè)水生態(tài)功能區(qū)大型底棲動(dòng)物指數(shù)比較Tab.3 Comparison of macrobenthic indices in three freshwater ecosystem function regions in Huntai River Basin

    利用Origin 85箱體圖對(duì)單因素方差分析中3個(gè)水生態(tài)功能區(qū)差異極顯著的環(huán)境因子做進(jìn)一步比較分析。從圖2可見:3個(gè)水生態(tài)功能區(qū)的海拔、總?cè)芙夤腆w、高錳酸鉀鹽指數(shù)的平均值差異性比較明顯;而水生態(tài)功能Ⅱ區(qū)和水生態(tài)功能Ⅲ區(qū)的電導(dǎo)率、鹽度、溶解氧、水深、總氮、硬度緩沖區(qū)流域總面積、建設(shè)用地面積、旱田面積的平均值差異性不明顯,但水生態(tài)功能Ⅰ區(qū)與水生態(tài)功能Ⅱ區(qū)、Ⅲ區(qū)平均值差異性相比較明顯。

    2.3 不同水生態(tài)功能區(qū)大型底棲動(dòng)物與環(huán)境因子的關(guān)系

    選擇典范對(duì)應(yīng)分析 (CCA)結(jié)果顯示:渾太河流域3個(gè)水生態(tài)功能區(qū)內(nèi)對(duì)大型底棲動(dòng)物群落分布影響的環(huán)境因子具有差異性。水生態(tài)功能Ⅰ區(qū)主要的驅(qū)動(dòng)因子是海拔、水深、總氮;水生態(tài)功能Ⅱ區(qū)主要的驅(qū)動(dòng)因子是高錳酸鉀鹽指數(shù)、鹽度;水生態(tài)功能Ⅲ區(qū)主要的驅(qū)動(dòng)因子是溶解氧、電導(dǎo)率、鹽度 (圖3)。

    皮爾森相關(guān)系數(shù)分析結(jié)果 (表5)顯示:水生態(tài)功能Ⅰ區(qū)海拔與Shannon-Wiener多樣性指數(shù)呈顯著正相關(guān) (P<0.05),總氮與Shannon-Wiener多樣性指數(shù)呈顯著負(fù)相關(guān) (P<0.05),高錳酸鉀鹽與Pielou均勻度指數(shù)呈顯著負(fù)相關(guān) (P<0.05),海拔與物種豐富度呈極顯著正相關(guān) (P<0.01),總氮與物種豐富度呈極顯著負(fù)相關(guān) (P<0.01),溶解氧與物種豐富度呈顯著正相關(guān) (P<0.05);水生態(tài)功能Ⅱ區(qū)環(huán)境因子與Shannon-Wiener指數(shù)不相關(guān)(P>0.05),海拔與Pielou均勻度指數(shù)呈顯著負(fù)相關(guān) (P<0.05),電導(dǎo)率、鹽度、總?cè)芙夤腆w與Pielou均勻度指數(shù)呈顯著正相關(guān) (P<0.05),海拔與物種豐富度呈顯著正相關(guān) (P<0.05);水生態(tài)功能Ⅲ區(qū)環(huán)境因子與這些指數(shù)相關(guān)性不明顯 (P>0.05)。據(jù)此表明,渾太河流域3個(gè)水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)與環(huán)境因子的相關(guān)性具有差異性。

    表4 不同水生態(tài)功能區(qū)的環(huán)境因子比較Tab.4 Comparison of environmental factors in different freshwater ecosystem function regions

    圖2 渾太河流域3個(gè)水生態(tài)功能區(qū)環(huán)境特征的比較Fig.2 Environmental characteristics in three freshwater ecosystem function regions in Huntai River Basin

    圖3 不同水生態(tài)功能區(qū)環(huán)境因子與大型底棲動(dòng)物群落的典范對(duì)應(yīng)分析結(jié)果Fig.3 Canonical correspondence analysis(CCA)of species-environmental relatioship in different freshwater ecosystem function regions

    表5 皮爾森相關(guān)系數(shù)分析3個(gè)水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)與環(huán)境因子相關(guān)性的結(jié)果Tab.5 Pearson correlation coefficient analysis of macrobenthic community structure and environmental factors related to the results in three aquatic ecological functional regions

    3 討論

    3.1 不同水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)的差異

    在渾太河流域共采集鑒定大型底棲動(dòng)物71屬(種),隸屬7綱13目35科,其中,水生昆蟲51種,占71.83%。這與渠曉東等[15]對(duì)太子河流域大型底棲動(dòng)物調(diào)查結(jié)果 (大型底棲動(dòng)物共鑒定131種,主要為水生昆蟲)相一致。渾太河流域不同水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)表現(xiàn)出明顯的空間異質(zhì)性,3個(gè)水生態(tài)功能區(qū)的大型底棲動(dòng)物物種組成比例存在明顯差異,指示種有顯著性差異,Shannon-Wiener多樣性指數(shù)、Pielou均勻度和物種豐富度在各個(gè)水生態(tài)功能區(qū)間有顯著性差異,其中物種豐富度呈現(xiàn)遞增的趨勢。這可能是因?yàn)樗飳?duì)不同棲息地的適應(yīng)機(jī)制不同,造成3個(gè)水生態(tài)功能區(qū)大型底棲動(dòng)物群落結(jié)構(gòu)的差異[16]。3個(gè)水生態(tài)功能區(qū)由平原區(qū)逐步過渡到山地森林區(qū),生境由簡單到復(fù)雜,人為干擾程度逐漸減小等一系列因素造成物種豐富度呈現(xiàn)出遞增的趨勢;水生態(tài)功能Ⅰ區(qū)和Ⅱ區(qū)為平原和丘陵區(qū),建設(shè)用地和耕地所占面積相對(duì)水生態(tài)功能Ⅲ區(qū)比例較大,而Shannon-Wiener多樣性指數(shù)低于水生態(tài)功能Ⅲ區(qū),流域內(nèi)建設(shè)用地和耕地比例的增加會(huì)加劇水體惡化,從而導(dǎo)致大型底棲動(dòng)物多樣性下降[17-18];不同的水生態(tài)功能區(qū)河流底質(zhì)類型、理化因子也不盡相同,不同的河流流域生境也造成大型底棲動(dòng)物群落結(jié)構(gòu)的不同,因?yàn)榱饔虻踪|(zhì)不僅是大型底棲動(dòng)物應(yīng)對(duì)生物和非生物擾動(dòng)的緩沖區(qū)[19],而且是大型底棲動(dòng)物的棲息和產(chǎn)卵場所[20],所以底質(zhì)對(duì)大型底棲動(dòng)物的分布具有重要作用,不同的水生態(tài)功能區(qū)底質(zhì)類型不同,使得大型底棲動(dòng)物的群落結(jié)構(gòu)具有顯著差異。由于3個(gè)水生態(tài)功能區(qū)人為干擾程度和土地利用類型等均有差異,造成3個(gè)功能區(qū)的指示種呈現(xiàn)顯著差異。3個(gè)水生態(tài)功能區(qū)的Pielou均勻度指數(shù)、Shannon-Wiener多樣性指數(shù)和Species richness指數(shù)有顯著或極顯著性差異,這與不同水生態(tài)功能區(qū)的城鎮(zhèn)化程度、人為干擾活動(dòng)、污染狀況、土地利用類型等均有關(guān)系,這些都直接或間接地影響水生生物的生長,對(duì)太子河、渭河的有關(guān)研究中也得到了證實(shí)[21-27]。

    3.2 不同水生態(tài)功能區(qū)環(huán)境因子的差異

    渾太河流域不同水生態(tài)功能區(qū)間環(huán)境因子存在較大的差異,這與王云濤等[28]得出的太子河流域不同水生態(tài)區(qū)間環(huán)境特征存在較大差異相一致。本研究中發(fā)現(xiàn):3個(gè)水生態(tài)功能區(qū)間的水溫、河寬、流速、流量、總磷、森林面積、水田面積等環(huán)境因子存在顯著性差異,水生態(tài)功能Ⅲ區(qū)到Ⅰ區(qū)地貌類型從山地到平原,河道由窄變寬,水溫也越來越高;海拔、電導(dǎo)率、鹽度、溶解氧、總?cè)芙夤腆w、水深、總氮、高錳酸鉀鹽指數(shù)、硬度、緩沖區(qū)流域總面積、建設(shè)用地面積、旱田面積等環(huán)境因子在3個(gè)水生態(tài)功能區(qū)間存在極顯著性差異。在水生態(tài)功能Ⅰ區(qū)、Ⅱ區(qū)為低海拔區(qū),適宜居住和開展工農(nóng)業(yè)生產(chǎn),然而建筑用地面積、水田和旱田面積較高,人為干擾較大,水體污染比較嚴(yán)重;水生態(tài)功能Ⅲ區(qū)主要是山地森林,不利于人類活動(dòng),土地利用類型就與其他兩個(gè)區(qū)域有較大差異,導(dǎo)致水體理化因子出現(xiàn)梯度變化,水生態(tài)功能Ⅰ區(qū)、Ⅱ區(qū)、Ⅲ區(qū)的海拔依次增加,而水溫、電導(dǎo)率、鹽度、總?cè)芙夤腆w、河寬、流量、總氮、高錳酸鉀鹽指數(shù)則呈現(xiàn)依次降低的趨勢。有研究表明,電導(dǎo)率與耕地面積和建設(shè)面積呈線性相關(guān)[29],水生態(tài)功能Ⅰ區(qū)、Ⅱ區(qū)、Ⅲ區(qū)的自然地貌分別為平原區(qū) (平均海拔為29.23 m)、丘陵區(qū) (平均海拔為245 m)、山地區(qū) (平均海拔為349.6 m),由于自然地貌的不同,不同水生態(tài)功能區(qū)的城鎮(zhèn)化程度、土地利用類型也不同,城鎮(zhèn)化對(duì)河流生態(tài)系統(tǒng)的影響通過多種因子共同作用,導(dǎo)致環(huán)境因素存在顯著性差異[30-31]。

    3.3 不同水生態(tài)功能區(qū)大型底棲動(dòng)物群落及其驅(qū)動(dòng)因子的比較

    渾太河流域大型底棲動(dòng)物群落結(jié)構(gòu)受到環(huán)境因子的影響,且不同水生態(tài)功能區(qū)影響大型底棲動(dòng)物群落結(jié)構(gòu)的環(huán)境因子差異顯著:海拔、水深、總氮對(duì)水生態(tài)功能Ⅰ區(qū)內(nèi)大型底棲動(dòng)物群落結(jié)構(gòu)影響較大;高錳酸鉀鹽指數(shù)、鹽度對(duì)水生態(tài)功能Ⅱ區(qū)內(nèi)大型底棲動(dòng)物群落結(jié)構(gòu)影響較大;溶解氧、電導(dǎo)率、鹽度對(duì)水生態(tài)功能Ⅲ區(qū)內(nèi)大型底棲動(dòng)物群落結(jié)構(gòu)影響較大。Marshall等[32]認(rèn)為,環(huán)境因子、生物之間的相互作用以及生物因子是影響生物群落結(jié)構(gòu)的3個(gè)要素。有相關(guān)研究顯示,緯度、海拔對(duì)大型底棲動(dòng)物群落結(jié)構(gòu)有顯著影響[33],物種豐富度隨海拔的升高而降低,導(dǎo)致這種機(jī)制形成的原因是隨著海拔的變化,環(huán)境因子、土地利用類型均有差異[34-36],從而影響了大底棲動(dòng)物群落結(jié)構(gòu)的空間異質(zhì)性。眾多研究發(fā)現(xiàn),海拔差異會(huì)導(dǎo)致土地利用、人為活動(dòng)等因素的差異,一系列因素對(duì)底棲動(dòng)物群落結(jié)構(gòu)具有顯著的影響[37]。研究表明,總氮是河流生態(tài)系統(tǒng)重要的水化指標(biāo),能顯著影響大型底棲動(dòng)物的群落結(jié)構(gòu)[38-39]。同時(shí),不同的水生態(tài)功能區(qū)大型底棲動(dòng)物指數(shù)與環(huán)境因子的相關(guān)性差異也十分明顯,有研究顯示,由于不同水生態(tài)功能區(qū)水質(zhì)、水文、棲息地質(zhì)量等的差異,造成對(duì)大型底棲動(dòng)物群落結(jié)構(gòu)有顯著相關(guān)性的環(huán)境因子也不盡相同[40-42],以致于不同水生態(tài)功能區(qū)對(duì)大型底棲動(dòng)物的群落結(jié)構(gòu)影響相關(guān)的環(huán)境因子存在顯著差異。不同水生態(tài)功能區(qū)具有不同的河流底質(zhì)類型,不同的底質(zhì)類型對(duì)水生生物群落結(jié)構(gòu)的影響也不盡相同。在對(duì)渭河和渾河的相關(guān)研究中發(fā)現(xiàn),不同的底質(zhì)類型,著生藻的群落結(jié)構(gòu)具有明顯差異[43-44]?;诓煌鷳B(tài)功能分區(qū)進(jìn)行水生生物群落結(jié)構(gòu)及其環(huán)境因子的相關(guān)性差異的比較,有利于在流域尺度上對(duì)不同水生態(tài)功能區(qū)水生生物群落結(jié)構(gòu)與環(huán)境因子,以及水生生物與環(huán)境因子相關(guān)性的進(jìn)一步認(rèn)識(shí),從水生生物層面上更深入地理解水生態(tài)功能區(qū),以便對(duì)水生態(tài)健康的調(diào)查評(píng)估提供指導(dǎo)。

    [1] 張楠,張遠(yuǎn),孔維靜,等.太子河流域水生態(tài)功能Ⅱ級(jí)區(qū)的劃分[J].環(huán)境科學(xué)研究,2013,26(5):472-479.

    [2] 孟偉,張遠(yuǎn),張楠,等.流域水生態(tài)功能區(qū)概念、特點(diǎn)與實(shí)施策略[J].環(huán)境科學(xué)研究,2013,26(5):465-471.

    [3] 孟偉,張遠(yuǎn),張楠,等.流域水生態(tài)功能分區(qū)與質(zhì)量目標(biāo)管理技術(shù)研究的若干問題[J].環(huán)境科學(xué)學(xué)報(bào),2011,31(7):1345-1351.

    [4] 高永年,高俊峰.太湖流域水生態(tài)功能分區(qū)[J].地理研究,2010,29(1):111-117.

    [5] Omernik J M,Bailey G B.Distinguishing between watersheds and ecoregions[J].Journal of the American Water Resources Association,1997,33(5):935-949.

    [6] 孟偉,張遠(yuǎn),鄭丙輝.水生態(tài)區(qū)劃方法及其在中國的應(yīng)用前景[J].水科學(xué)進(jìn)展,2007,18(2):293-300.

    [7] 馬溪平,周世嘉,張遠(yuǎn),等.流域水生態(tài)功能分區(qū)方法與指標(biāo)體系探討[J].環(huán)境科學(xué)與管理,2010,35(12):59-64,70.

    [8] 孫小銀,周啟星.中國水生態(tài)分區(qū)初探[J].環(huán)境科學(xué)學(xué)報(bào),2010,30(2):415-423.

    [9] Barbour M T,Gerritsen J,Snyder B D,et al.Rapid Bioassessment Protocols for Use in Streams and Wadeable Rivers: Periphyton,Benthic Macroinvertebrates and Fish[M].2nd ed.Washington DC:U.S.Environmental Protecting Agency,Office of Water,1999.

    [10] Covich A P,Palmer M A,Crowl T A.The role of benthic invertebrate species in freshwater ecosystems:zoobenthic species influence energy flows and nutrient cycling[J].Bio Science,1999,49(2):119-127.

    [11] Devine J A,Vanni M J.Spatial and seasonal variation in nutrient excretion by benthic invertebrates in a eutrophic reservoir[J].Freshwater Biology,2002,47(6):1107-1121.

    [12] Pedersen M L,F(xiàn)riberg N,Skriver J,et al.Restoration of Skjern River and its valley-short-term effects on river habitats,macrophytes and macroinvertebrates[J].Ecological Engineering,2007,30(2):145-156.

    [13] Brown S C,Smith K,Batzer D.Macroinvertebrate responses to wetland restoration in northern New York[J].Environmental Entomology,1997,26(5):1016-1024.

    [14] 李艷利,李艷粉,徐宗學(xué).影響渾太河流域大型底棲動(dòng)物群落結(jié)構(gòu)的環(huán)境因子分析[J].環(huán)境科學(xué),2015,36(1):94-106.

    [15] 渠曉東,張遠(yuǎn),馬淑芹,等.太子河流域大型底棲動(dòng)物群落結(jié)構(gòu)空間分布特征[J].環(huán)境科學(xué)研究,2013,26(5):509-515.

    [16] 閔文武,王培培,李麗娟,等.渭河流域浮游植物功能群與環(huán)境因子的關(guān)系[J].環(huán)境科學(xué)研究,2015,28(9):1397-1406.

    [17] Beauger A,Lair N,Reyes-Marchant P,et al.The distribution of macroinvertebrate assemblages in a reach of the River Allier(France),in relation to riverbed characteristics[J].Hydrobiologia,2006,571(1):63-76.

    [18] Jun Y C,Kim N Y,Kwon S J,et al.Effects of land use on benthic macroinvertebratecommunities:comparison oftwomountain streams in Korea[J].Annales de Limnologie-International Journal of Limnology,2011,47(S1):S35-S49.

    [19] Beisel J N,Usseglio-Polatera P,Thomas S,et al.Stream community structure in relation to spatial variation:the influence of mesohabitat characteristics[J].Hydrobiologia,1998,389(1-3):73-88.

    [20] Jiang X M,Xiong J,Qiu J W,et al.Structure of macroinvertebrate communities in relation to environmental variables in a subtropical Asian river system[J].International Review of Hydrobiology,2010,95(1):42-57.

    [21] 殷旭旺,張遠(yuǎn),渠曉東,等.太子河著生藻類群落結(jié)構(gòu)空間分布特征[J].環(huán)境科學(xué)研究,2013,26(5):502-508.

    [22] 丁森,張遠(yuǎn),渠曉東,等.影響太子河流域魚類空間分布的不同尺度環(huán)境因子分析[J].環(huán)境科學(xué),2012,33(7):2272-2280.

    [23] 李麗娟,金文,王博涵,等.太子河河岸帶土地利用類型與硅藻群落結(jié)構(gòu)的關(guān)系[J].環(huán)境科學(xué)研究,2015,28(11):1662-1669.

    [24] 殷旭旺,渠曉東,李慶南,等.基于著生藻類的太子河流域水生態(tài)系統(tǒng)健康評(píng)價(jià)[J].生態(tài)學(xué)報(bào),2012,32(6):1677-1691.

    [25] 劉麟菲,譚冰冰,殷旭旺,等.在種與屬兩個(gè)級(jí)別評(píng)價(jià)太子河硅藻群落與環(huán)境因子的關(guān)系[J].生態(tài)學(xué)報(bào),2014,34(22):6 613-6621.

    [26] 劉麟菲,宋佳,王博涵,等.渭河流域硅藻群落特征及水生態(tài)健康評(píng)價(jià)[J].環(huán)境科學(xué)研究,2015,28(10):1560-1569.

    [27] 李麗娟,張吉,吳丹,等.太子河流域魚類功能群結(jié)構(gòu)與多樣性對(duì)土地利用類型的響應(yīng)[J].生態(tài)學(xué)報(bào),2017,37(20):6863-6874.

    [28] 王云濤,張遠(yuǎn),高欣,等.太子河流域不同水生態(tài)區(qū)魚類群落分布與環(huán)境因子的關(guān)聯(lián)性[J].環(huán)境科學(xué)研究,2016,29(2):192-201.

    [29] Bu H M,Meng W,Zhang Y,et al.Relationships between land use patterns and water quality in the Taizi River Basin,China[J].Ecological Indicators,2014(41):187-197.

    [30] Sprague L A,Zuellig R E,Dupree J A.Effects of urban development on stream ecosystems along the front range of the Rocky Mountains,Colorado and Wyoming[R].USGS,2006:134-145.

    [31] Helms B S,F(xiàn)eminella J W,Pan S F.Detection of biotic responses to urbanization using fish assemblages from small streams of western Georgia,USA[J].Urban Ecosystems,2005,8(1):39-57.

    [32] Marshall S,Elliott M.Environmental influences on the fish assemblage of the Humber estuary,U.K[J].Estuarine,Coastal and Shelf Science,1998,46(2):175-184.

    [33] Jacobsen D,Schultz R,Encalada A.Structure and diversity of stream invertebrate assemblages:the influence of temperature with altitude and latitude[J].Freshwater Biology,1997,38(2):247-261.

    [34] Rundle S D,Jenkins A,Ormerod S J.Macroinvertebrate communities in streams in the Himalaya,Nepal[J].Freshwater Biology,1993,30(1):169-180.

    [35] Ao M,Alfred J R B,Gupta A.Studies on some lotic systems in the north-eastern hill region of India[J].Limnologica,1984,15:135-141.

    [36] Suren A M.Macroinvertebrate communities of streams in western Nepal:effects of altitude and land use[J].Freshwater Biology,1994,32(2):323-336.

    [37] 張勇,劉朔孺,于海燕,等.錢塘江中游流域不同空間尺度環(huán)境因子對(duì)底棲動(dòng)物群落的影響[J].生態(tài)學(xué)報(bào),2012,32(14):4309-4317.

    [38] 吳東浩,于海燕,吳海燕,等.基于大型底棲無脊椎動(dòng)物確定河流營養(yǎng)鹽濃度閾值——以西苕溪上游流域?yàn)槔齕J].應(yīng)用生態(tài)學(xué)報(bào),2010,21(2):483-488.

    [39] 謝旭,俞存根,鄭基,等.魚山漁場披山海域蝦類群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系[J].大連海洋大學(xué)學(xué)報(bào),2017,32(6):713-723.

    [40] Gao X,Niu C J,Chen Y S,et al.Spatial heterogeneity of stream environmental conditions and macroinvertebrates community in an agriculture dominated watershed and management implications for a large river(the Liao River,China)basin[J].Environmental Monitoring and Assessment,2014,186(4):2375-2391.

    [41] Qu X D,Tang T,Xie Z C,et al.Distribution of the macroinvertebrate communities in the Xiangxi River system and relationships with environmental factors[J].Journal of Freshwater Ecology,2005,20(2):233-238.

    [42] Jowett I G.Hydraulic constraints on habitat suitability for benthic invertebrates in gravel-bed rivers[J].River Research and Applications,2003,19(5-6):495-507.

    [43] 殷旭旺,徐宗學(xué),鄢娜,等.渭河流域河流著生藻類的群落結(jié)構(gòu)與生物完整性研究[J].環(huán)境科學(xué)學(xué)報(bào),2013,33(2):518-527.

    [44] 項(xiàng)珍龍,楊璐,李晨,等.渾河流域底質(zhì)類型對(duì)硅藻群落分布特征的影響[J].環(huán)境科學(xué)學(xué)報(bào),doi:10.13671/j.hjkxxb.2017.0320.

    猜你喜歡
    功能區(qū)流域動(dòng)物
    壓油溝小流域
    堡子溝流域綜合治理
    羅堰小流域
    燒結(jié)混合機(jī)各功能區(qū)的研究與優(yōu)化
    山東冶金(2019年3期)2019-07-10 00:53:54
    水利規(guī)劃與設(shè)計(jì)(2018年1期)2018-01-31 01:53:37
    衡水市不同功能區(qū)土壤營養(yǎng)元素的比較
    最萌動(dòng)物榜
    動(dòng)物可笑堂
    水功能區(qū)納污能力計(jì)算及污染物總量控制
    動(dòng)物捉迷藏
    校园春色视频在线观看| 大陆偷拍与自拍| 国产91精品成人一区二区三区| 色在线成人网| 波多野结衣高清无吗| 欧美日韩黄片免| 欧美成人午夜精品| 国产精品乱码一区二三区的特点 | 一区在线观看完整版| a级毛片在线看网站| 最近最新中文字幕大全免费视频| 日本a在线网址| 亚洲五月天丁香| 91国产中文字幕| 夜夜爽天天搞| 桃红色精品国产亚洲av| 免费看十八禁软件| 国产精华一区二区三区| 一进一出抽搐动态| 99精品久久久久人妻精品| 一区在线观看完整版| 久久这里只有精品19| 狂野欧美激情性xxxx| 中文字幕最新亚洲高清| 国产精品一区二区免费欧美| 黄色怎么调成土黄色| 香蕉丝袜av| 久久狼人影院| 国产高清激情床上av| 国产免费av片在线观看野外av| 久久天堂一区二区三区四区| 亚洲激情在线av| 日韩精品青青久久久久久| 国产野战对白在线观看| 国产在线观看jvid| 精品久久久久久,| 在线观看免费视频网站a站| 久久天堂一区二区三区四区| 亚洲精品粉嫩美女一区| 99国产精品一区二区蜜桃av| 人人妻人人爽人人添夜夜欢视频| 亚洲精品美女久久av网站| 国内毛片毛片毛片毛片毛片| 久久久久久久久久久久大奶| 国产精华一区二区三区| 人人妻人人澡人人看| 国产成人欧美在线观看| 久久国产精品人妻蜜桃| 麻豆一二三区av精品| 精品国产乱子伦一区二区三区| 日本三级黄在线观看| 国产熟女午夜一区二区三区| 久久午夜亚洲精品久久| 国产精品亚洲av一区麻豆| 久久人妻av系列| 好看av亚洲va欧美ⅴa在| 久热爱精品视频在线9| 人人妻,人人澡人人爽秒播| 丰满的人妻完整版| 亚洲人成77777在线视频| 精品福利永久在线观看| 亚洲色图 男人天堂 中文字幕| 人人妻人人澡人人看| 欧美另类亚洲清纯唯美| 黄色a级毛片大全视频| 国产成人免费无遮挡视频| 女性被躁到高潮视频| 黄色成人免费大全| 一进一出抽搐gif免费好疼 | 久久久久久久久久久久大奶| 欧美日韩福利视频一区二区| 午夜老司机福利片| 亚洲avbb在线观看| 看片在线看免费视频| 成年版毛片免费区| 人人妻,人人澡人人爽秒播| 91麻豆精品激情在线观看国产 | 91九色精品人成在线观看| 热99国产精品久久久久久7| 激情视频va一区二区三区| 亚洲 国产 在线| 国产高清激情床上av| videosex国产| 一个人观看的视频www高清免费观看 | 69精品国产乱码久久久| 精品欧美一区二区三区在线| 中文字幕av电影在线播放| 热re99久久国产66热| 日本黄色视频三级网站网址| 纯流量卡能插随身wifi吗| 色在线成人网| 亚洲视频免费观看视频| 90打野战视频偷拍视频| 日韩欧美免费精品| 波多野结衣一区麻豆| 桃色一区二区三区在线观看| 中文字幕av电影在线播放| 一边摸一边抽搐一进一出视频| 精品日产1卡2卡| 一级片'在线观看视频| 中文字幕av电影在线播放| 国产成+人综合+亚洲专区| 999精品在线视频| 午夜福利一区二区在线看| 黑丝袜美女国产一区| 欧美乱妇无乱码| 黄色丝袜av网址大全| 老司机深夜福利视频在线观看| 亚洲成人久久性| 他把我摸到了高潮在线观看| 高清黄色对白视频在线免费看| 两人在一起打扑克的视频| 国产1区2区3区精品| 日韩中文字幕欧美一区二区| 亚洲少妇的诱惑av| 久久中文看片网| 中文字幕另类日韩欧美亚洲嫩草| 亚洲一区二区三区欧美精品| 国产精品亚洲av一区麻豆| 看免费av毛片| 手机成人av网站| 12—13女人毛片做爰片一| 男女做爰动态图高潮gif福利片 | 老司机靠b影院| 久久精品aⅴ一区二区三区四区| 黑丝袜美女国产一区| 国产极品粉嫩免费观看在线| 成人三级做爰电影| 久久久久国产精品人妻aⅴ院| 日韩精品青青久久久久久| 久久久水蜜桃国产精品网| 国产成人精品久久二区二区91| 丝袜人妻中文字幕| aaaaa片日本免费| 淫妇啪啪啪对白视频| 亚洲专区国产一区二区| 日韩精品免费视频一区二区三区| 久久中文字幕人妻熟女| 中国美女看黄片| 国产亚洲欧美98| 多毛熟女@视频| 精品国产国语对白av| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美最黄视频在线播放免费 | 他把我摸到了高潮在线观看| 成人亚洲精品一区在线观看| 嫁个100分男人电影在线观看| 大码成人一级视频| 色综合婷婷激情| 久久久久久免费高清国产稀缺| 一区二区日韩欧美中文字幕| 精品免费久久久久久久清纯| 久久精品国产99精品国产亚洲性色 | 天堂中文最新版在线下载| 免费在线观看亚洲国产| 国产激情欧美一区二区| 777久久人妻少妇嫩草av网站| 亚洲国产欧美日韩在线播放| 老汉色∧v一级毛片| 天天躁夜夜躁狠狠躁躁| 国产高清激情床上av| 中文字幕人妻丝袜制服| 欧美日韩一级在线毛片| 最新美女视频免费是黄的| www日本在线高清视频| 乱人伦中国视频| 一二三四社区在线视频社区8| 久久久久国内视频| 成熟少妇高潮喷水视频| 好看av亚洲va欧美ⅴa在| 日韩免费av在线播放| 亚洲,欧美精品.| 日韩人妻精品一区2区三区| 大码成人一级视频| 欧美精品啪啪一区二区三区| 日本欧美视频一区| 久久亚洲真实| 咕卡用的链子| 日韩欧美一区二区三区在线观看| 在线观看免费午夜福利视频| 少妇 在线观看| 亚洲中文av在线| 一级片'在线观看视频| 亚洲一区二区三区色噜噜 | 亚洲精品久久午夜乱码| 两个人免费观看高清视频| 在线观看免费午夜福利视频| 又紧又爽又黄一区二区| 色综合站精品国产| 久久人妻av系列| 欧美日本亚洲视频在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜激情av网站| 国产午夜精品久久久久久| 大码成人一级视频| 国产一区二区三区视频了| 国产有黄有色有爽视频| 我的亚洲天堂| 国产高清激情床上av| 国产精品美女特级片免费视频播放器 | 欧美成人午夜精品| 亚洲av第一区精品v没综合| 我的亚洲天堂| 午夜福利一区二区在线看| 国产国语露脸激情在线看| 欧美激情久久久久久爽电影 | 欧美日韩福利视频一区二区| 亚洲专区国产一区二区| 一夜夜www| 91老司机精品| av视频免费观看在线观看| bbb黄色大片| 成年版毛片免费区| 亚洲性夜色夜夜综合| 91大片在线观看| 色在线成人网| 真人一进一出gif抽搐免费| 又黄又粗又硬又大视频| 制服人妻中文乱码| 欧美在线黄色| 桃色一区二区三区在线观看| 少妇粗大呻吟视频| 午夜免费激情av| 又黄又粗又硬又大视频| 女性生殖器流出的白浆| 淫妇啪啪啪对白视频| 亚洲第一av免费看| 97超级碰碰碰精品色视频在线观看| 亚洲专区国产一区二区| 欧美成人性av电影在线观看| 精品国产国语对白av| 黄片小视频在线播放| 久久久国产成人精品二区 | 欧美丝袜亚洲另类 | 亚洲精品一二三| av片东京热男人的天堂| 纯流量卡能插随身wifi吗| 亚洲国产欧美一区二区综合| 亚洲专区字幕在线| 国产免费现黄频在线看| 高潮久久久久久久久久久不卡| 精品国产一区二区三区四区第35| 亚洲免费av在线视频| 亚洲精品国产区一区二| www.www免费av| 日韩免费高清中文字幕av| 国产精品二区激情视频| 国内久久婷婷六月综合欲色啪| 日本五十路高清| 国产精品偷伦视频观看了| 不卡av一区二区三区| 久久久水蜜桃国产精品网| 日本三级黄在线观看| 91成人精品电影| 欧美人与性动交α欧美精品济南到| 亚洲午夜精品一区,二区,三区| 90打野战视频偷拍视频| 99精品在免费线老司机午夜| 国产成年人精品一区二区 | 亚洲av成人av| 久久热在线av| 男女做爰动态图高潮gif福利片 | 人妻丰满熟妇av一区二区三区| 成人三级黄色视频| 老司机午夜十八禁免费视频| 大型黄色视频在线免费观看| 久久性视频一级片| 国产亚洲av高清不卡| 18禁美女被吸乳视频| 成人精品一区二区免费| 亚洲精品国产精品久久久不卡| 日韩欧美在线二视频| 免费看a级黄色片| 真人一进一出gif抽搐免费| 国产精品99久久99久久久不卡| 国产精品av久久久久免费| 看免费av毛片| 99国产精品免费福利视频| √禁漫天堂资源中文www| 国产亚洲精品第一综合不卡| 老司机午夜十八禁免费视频| 国产深夜福利视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲中文字幕日韩| 国产成人影院久久av| 欧美精品啪啪一区二区三区| 久久国产精品男人的天堂亚洲| 亚洲狠狠婷婷综合久久图片| 两性夫妻黄色片| 黄片播放在线免费| 精品国产美女av久久久久小说| 十分钟在线观看高清视频www| 欧美日韩亚洲高清精品| 香蕉国产在线看| 老熟妇乱子伦视频在线观看| 满18在线观看网站| 欧美精品啪啪一区二区三区| 51午夜福利影视在线观看| 亚洲中文日韩欧美视频| 欧美激情极品国产一区二区三区| 亚洲色图综合在线观看| 在线免费观看的www视频| 亚洲国产中文字幕在线视频| 午夜a级毛片| 欧美中文日本在线观看视频| 午夜影院日韩av| 18禁黄网站禁片午夜丰满| 精品一区二区三区视频在线观看免费 | 国产亚洲欧美精品永久| 在线观看一区二区三区激情| 成人国语在线视频| 波多野结衣av一区二区av| 精品无人区乱码1区二区| 老司机靠b影院| 大香蕉久久成人网| 夜夜看夜夜爽夜夜摸 | 一本大道久久a久久精品| 久久久国产精品麻豆| 脱女人内裤的视频| 国产午夜精品久久久久久| 国产成人av激情在线播放| 他把我摸到了高潮在线观看| 正在播放国产对白刺激| 一进一出好大好爽视频| www.熟女人妻精品国产| av视频免费观看在线观看| 亚洲avbb在线观看| 亚洲精品久久成人aⅴ小说| 久久久久国产一级毛片高清牌| 女人高潮潮喷娇喘18禁视频| 最新在线观看一区二区三区| 在线视频色国产色| 久久精品国产99精品国产亚洲性色 | 精品国产亚洲在线| 成人手机av| 久久国产精品影院| 18禁裸乳无遮挡免费网站照片 | 高清在线国产一区| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 免费在线观看完整版高清| 天堂影院成人在线观看| 两个人免费观看高清视频| 精品久久久久久成人av| 久久精品国产亚洲av香蕉五月| 校园春色视频在线观看| 日韩欧美一区视频在线观看| 国产精品 欧美亚洲| avwww免费| 精品日产1卡2卡| 巨乳人妻的诱惑在线观看| 在线av久久热| 啦啦啦免费观看视频1| 欧美中文日本在线观看视频| 国产一区在线观看成人免费| 午夜成年电影在线免费观看| 午夜免费观看网址| 久久人人爽av亚洲精品天堂| 国产欧美日韩精品亚洲av| 人人澡人人妻人| 丰满的人妻完整版| 精品久久久久久久久久免费视频 | www.精华液| 一a级毛片在线观看| 人妻久久中文字幕网| 久久久国产欧美日韩av| 成人三级做爰电影| 少妇被粗大的猛进出69影院| 两个人看的免费小视频| 老汉色av国产亚洲站长工具| 国产成人免费无遮挡视频| 女性生殖器流出的白浆| 咕卡用的链子| а√天堂www在线а√下载| 在线十欧美十亚洲十日本专区| 欧美激情高清一区二区三区| 色老头精品视频在线观看| 69av精品久久久久久| 欧美激情久久久久久爽电影 | 欧美日韩精品网址| 久久久久国内视频| 欧美亚洲日本最大视频资源| 免费在线观看黄色视频的| 好男人电影高清在线观看| 午夜老司机福利片| 久久国产乱子伦精品免费另类| 国产精品亚洲av一区麻豆| 波多野结衣高清无吗| 日本wwww免费看| 日本五十路高清| 最新在线观看一区二区三区| 日韩 欧美 亚洲 中文字幕| 成人亚洲精品av一区二区 | 国产99久久九九免费精品| 欧美激情极品国产一区二区三区| 亚洲欧美一区二区三区久久| 51午夜福利影视在线观看| 岛国在线观看网站| 色精品久久人妻99蜜桃| 999久久久精品免费观看国产| 国产欧美日韩一区二区精品| 99久久人妻综合| 十分钟在线观看高清视频www| 亚洲国产精品sss在线观看 | www.自偷自拍.com| 亚洲一区高清亚洲精品| 男女下面进入的视频免费午夜 | 亚洲aⅴ乱码一区二区在线播放 | 18禁黄网站禁片午夜丰满| 国产精品影院久久| 91麻豆av在线| 一进一出好大好爽视频| 免费在线观看日本一区| 国产一区二区三区在线臀色熟女 | 日韩高清综合在线| 免费在线观看影片大全网站| 免费不卡黄色视频| 一级a爱片免费观看的视频| 欧美另类亚洲清纯唯美| 欧美黄色片欧美黄色片| 美女福利国产在线| 一夜夜www| 水蜜桃什么品种好| 两性夫妻黄色片| 亚洲 欧美 日韩 在线 免费| 亚洲人成电影免费在线| 男女午夜视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲男人天堂网一区| 亚洲 欧美一区二区三区| 国产免费现黄频在线看| 欧美精品啪啪一区二区三区| 国产成人影院久久av| 免费在线观看影片大全网站| av片东京热男人的天堂| 国产高清激情床上av| 国产日韩一区二区三区精品不卡| 国产欧美日韩精品亚洲av| 色综合婷婷激情| 天堂中文最新版在线下载| 欧美精品亚洲一区二区| 亚洲色图 男人天堂 中文字幕| 好男人电影高清在线观看| 国产精品亚洲av一区麻豆| 精品一品国产午夜福利视频| 淫妇啪啪啪对白视频| 午夜福利免费观看在线| a级毛片在线看网站| 一二三四在线观看免费中文在| 在线观看免费日韩欧美大片| 韩国精品一区二区三区| 成人三级做爰电影| 在线视频色国产色| 国产麻豆69| 国产精品综合久久久久久久免费 | 亚洲av美国av| 美女扒开内裤让男人捅视频| 日本五十路高清| 99久久国产精品久久久| 999久久久精品免费观看国产| 久久精品成人免费网站| 美女国产高潮福利片在线看| 91av网站免费观看| 韩国av一区二区三区四区| 女性生殖器流出的白浆| 亚洲 欧美一区二区三区| 深夜精品福利| 久久午夜综合久久蜜桃| 在线永久观看黄色视频| 国产亚洲欧美精品永久| 亚洲第一av免费看| 在线av久久热| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 纯流量卡能插随身wifi吗| 欧美人与性动交α欧美精品济南到| 不卡一级毛片| 桃色一区二区三区在线观看| 在线永久观看黄色视频| 一级毛片高清免费大全| 男女下面插进去视频免费观看| 日本黄色视频三级网站网址| 亚洲成a人片在线一区二区| 久久精品aⅴ一区二区三区四区| 人人妻人人添人人爽欧美一区卜| 叶爱在线成人免费视频播放| 黑人巨大精品欧美一区二区mp4| 精品熟女少妇八av免费久了| 日本欧美视频一区| 黄色a级毛片大全视频| 亚洲性夜色夜夜综合| 国产精品一区二区精品视频观看| 99香蕉大伊视频| 国产在线观看jvid| 人人澡人人妻人| av在线天堂中文字幕 | 黄频高清免费视频| 激情视频va一区二区三区| 极品人妻少妇av视频| 日韩欧美一区视频在线观看| 看片在线看免费视频| 大型av网站在线播放| 妹子高潮喷水视频| 久久精品国产99精品国产亚洲性色 | 两性午夜刺激爽爽歪歪视频在线观看 | 国产91精品成人一区二区三区| 法律面前人人平等表现在哪些方面| 另类亚洲欧美激情| 亚洲国产欧美网| a级毛片在线看网站| 18禁观看日本| 午夜久久久在线观看| 国产一区二区三区综合在线观看| 热99国产精品久久久久久7| 国产精品成人在线| 欧美av亚洲av综合av国产av| 久久人妻av系列| 大陆偷拍与自拍| 女同久久另类99精品国产91| 精品久久久久久电影网| 欧美午夜高清在线| 韩国av一区二区三区四区| 亚洲一区二区三区欧美精品| 婷婷精品国产亚洲av在线| 国产片内射在线| 久久久久久大精品| 侵犯人妻中文字幕一二三四区| 视频区图区小说| 国产高清激情床上av| 久久 成人 亚洲| 啪啪无遮挡十八禁网站| 999久久久国产精品视频| www.精华液| 欧美日韩中文字幕国产精品一区二区三区 | 久久国产亚洲av麻豆专区| 欧美成人性av电影在线观看| 少妇粗大呻吟视频| 精品久久久久久成人av| 女性生殖器流出的白浆| 在线观看66精品国产| 国产麻豆69| 久久天堂一区二区三区四区| 国产熟女xx| 黄片大片在线免费观看| 亚洲美女黄片视频| 午夜免费鲁丝| 99国产精品一区二区蜜桃av| 天天躁夜夜躁狠狠躁躁| 丰满人妻熟妇乱又伦精品不卡| 90打野战视频偷拍视频| 欧美成人午夜精品| 亚洲成国产人片在线观看| 亚洲久久久国产精品| 99国产精品一区二区三区| 久久精品成人免费网站| 国产无遮挡羞羞视频在线观看| 夜夜爽天天搞| 国产麻豆69| 免费看a级黄色片| 亚洲国产中文字幕在线视频| 免费在线观看完整版高清| 91精品国产国语对白视频| 精品国产乱码久久久久久男人| 免费在线观看视频国产中文字幕亚洲| 午夜亚洲福利在线播放| 美女国产高潮福利片在线看| 18美女黄网站色大片免费观看| 国产亚洲精品久久久久久毛片| 国产av又大| 黑人巨大精品欧美一区二区蜜桃| 欧美中文日本在线观看视频| av欧美777| 免费在线观看日本一区| 嫩草影院精品99| 国产精品国产高清国产av| 久久久久久久久久久久大奶| 一区二区三区激情视频| 99精品久久久久人妻精品| 国产亚洲精品第一综合不卡| 亚洲精品国产一区二区精华液| av视频免费观看在线观看| 亚洲欧美日韩高清在线视频| 视频区欧美日本亚洲| 国产免费男女视频| 亚洲一区二区三区欧美精品| 黄色 视频免费看| 免费不卡黄色视频| 一区福利在线观看| 国产精品一区二区免费欧美| 国产精品日韩av在线免费观看 | 黄色怎么调成土黄色| 色综合站精品国产| 不卡av一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| av视频免费观看在线观看| 久久精品亚洲精品国产色婷小说| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久久人人做人人爽| 在线观看舔阴道视频| 亚洲色图 男人天堂 中文字幕| 国产主播在线观看一区二区| 欧美日韩国产mv在线观看视频| 久久午夜亚洲精品久久| 一区二区日韩欧美中文字幕| 日韩免费高清中文字幕av| 欧美日韩一级在线毛片| 免费av毛片视频| 91九色精品人成在线观看| 国产熟女午夜一区二区三区| 他把我摸到了高潮在线观看| 99国产精品一区二区三区| 国产成人av教育| 日本精品一区二区三区蜜桃| 亚洲色图av天堂| 久久久国产成人精品二区 |