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

    遼河流域大型底棲動物耐污值

    2015-01-19 05:48:36渠曉東劉思思
    生態(tài)學報 2015年14期
    關鍵詞:遼河流域樣點級別

    趙 瑞, 高 欣, 丁 森, 張 遠,* , 渠曉東, 劉思思

    1 中國環(huán)境科學研究院, 環(huán)境基準與風險評估國家重點實驗室, 北京 100012 2 中國環(huán)境科學研究院, 流域水生態(tài)研究室, 北京 100012 3 中國水利水電科學研究院, 水環(huán)境研究所, 北京 100038

    遼河流域大型底棲動物耐污值

    趙 瑞1,2, 高 欣1,2, 丁 森1,2, 張 遠1,2,*, 渠曉東3, 劉思思1,2

    1 中國環(huán)境科學研究院, 環(huán)境基準與風險評估國家重點實驗室, 北京 100012 2 中國環(huán)境科學研究院, 流域水生態(tài)研究室, 北京 100012 3 中國水利水電科學研究院, 水環(huán)境研究所, 北京 100038

    大型底棲動物耐污值是構建河流生物評價指標的基礎。受自然地理區(qū)域影響,同一分類單元的耐污值會出現地域上的差異,因此計算和修訂適用于某一研究區(qū)域的耐污值是準確評價該區(qū)域河流健康狀況的前提。目前尚未有針對我國北方河流大型底棲動物耐污值計算和修訂的報道,基于2009年春季、秋季和2010年夏季在遼河流域308個樣點的大型底棲動物數據,采用Simpson多樣性指數的水質分級標準對所有樣點的水質等級進行了劃分,并利用累積分位數法計算獲得了遼河流域大型底棲動物195個分類單元的耐污值;同時參考專家經驗和相關文獻核定了其余71個分類單元的耐污值。依據耐污值敏感特征劃分標準,遼河流域大型底棲動物優(yōu)勢類群以敏感種和一般耐污種為主,兩者幾乎占到總分類單元數的84%,耐污類群種類較少,僅占16%。其中敏感類群有齒蛉科Corydalidae,襀科Perlidae,原石蛾科Rhyacophilidae,細裳蜉科Leptophlebiidae等;一般耐污類群有紋石蛾科Hydropsychidae,四節(jié)蜉科Baetidae,扁蜷螺科Planorbidae等;耐污類群有顫蚓科Tubificidae、石蛭科Erpobdellidae等。

    大型底棲動物; 耐污值; 遼河流域

    耐污值(TV)是指生物在時間和空間上對環(huán)境壓力的耐受性,反映了生物在環(huán)境壓力干擾下其生存和繁殖相對能力的高低[1]。在河流生態(tài)學領域,耐污值最早是作為河流生物監(jiān)測指標的重要組分而被提出[2-3],主要應用于大型底棲動物監(jiān)測評價體系,所涉及指標如Hilsenhoff生物指數(HBI)[4]、北卡羅來納州生物指數(NCBI)[5]、敏感性/耐污性的豐富度綜合指數[6]等。隨著耐污值計算方法的不斷改進[5-7],研究者發(fā)現不同地理區(qū)域大型底棲動物耐污值存在差異。對于全球分布的大型底棲動物來說,受氣候、地理特征等環(huán)境因子影響,同一科所包含的屬或同一屬所包含的種在大尺度地理區(qū)域間會有所不同[5],進而會影響到同一科/屬的耐污值發(fā)生變化。美國國家環(huán)境保護局(USEPA,United States Environmental Protection Agency)在綜合考慮地理區(qū)域影響的基礎上,分別建立了美國西北部、中西部、東南部、中西部偏北和沿大西洋中部海岸區(qū)5個區(qū)域的大型底棲動物耐污值體系[6]。

    目前,我國大型底棲動物耐污值的研究報道多集中于中南部地區(qū),如:江西廬山、安徽南部山區(qū)、秦淮河和長江三角洲等[8-9],而我國北方地區(qū)的相關報道較少[10]。此外,以往的耐污值研究過于依賴專家經驗,或參與構建耐污值的樣本數量過少。針對上述大型底棲動物耐污值研究中的不足,本研究基于2009—2010年遼河全流域大型底棲動物基礎數據,初步推導了遼河流域大型底棲動物主要分類單元耐污值,為整個遼河流域水質生物評價提供理論依據和數據支持,也為我國東北地區(qū)大型底棲動物耐污值數據庫的建立提供數據支撐。

    1 材料和方法

    1.1 研究區(qū)域概況

    遼河流域位于中國東北地區(qū)南部(116°30′E—125°47′E,38°43′N—45°00′N),河流全長1345 km,流域面積約22.0×104km2。遼河流域年降水量 350—1000 mm,降水分布極不均勻,降水量自東南向西北逐漸減少,且多集中于 6—9 月,約占全年降水量的 70%,屬溫帶、暖溫帶半濕潤大陸性季風氣候。遼河流域上游為山丘區(qū),植被差,覆蓋率在30%以下,水土流失嚴重,是風沙干旱嚴重區(qū),東部為丘陵區(qū),森林茂密,植被覆蓋率高,中下游為平原區(qū)。遼河流域土地利用類型以農田為主(38.29%),草地(23.65%)和林地(24.3%)次之,裸地(6.45%),城市用地(4.43%)和水田(4.4%)所占比例最少,其中草地和林地主要分布在遼河流域的東南、西南和西北的一級河流區(qū),這些區(qū)域的棲境狀況也最佳。遼河流域各支流源頭區(qū)域的水質和棲境狀況良好,中下游河段COD和BOD5含量較高,受工業(yè)污染影響較大[11]。

    1.2 大型底棲動物樣品采集

    2009年春季、秋季和2010年夏季,在遼河流域共設置308個樣點(圖1),開展大型底棲動物野外監(jiān)測。利用索伯網(面積0.3 m×0.3 m,孔徑0.5 mm),選擇不同生境類型的采樣斷面收集3個平行樣品,將樣品沖洗并用60目篩網過濾,放入500 mL廣口瓶,加75%酒精保存。在實驗室條件下,利用體式顯微鏡(Olympus SZ61)將絕大多數物種鑒定到種或屬,搖蚊鑒定到亞科,寡毛類和線蟲鑒定到綱。

    圖1 遼河流域樣點圖

    1.3 水質級別劃分

    我國已有的多樣性指數水質劃分標準是黃玉瑤等[12]提出的,王備新等[9]在計算我國東部地區(qū)大型無脊椎動物耐污值時,根據Lenat[5]提出的五級水質劃分標準,對我國的Shannon-Wiener多樣性指數(H′)水質劃分標準進行了調整,即:H′≥4,極清潔;H′≥3,清潔;H′≥2,輕污;H′≥1,中污;H′<1,重污。本研究利用H′進行水質劃分發(fā)現樣點在各水質級別分布不均勻,而用修訂過的Simpson多樣性指數(D′)劃分的結果較好,故本研究采用D′劃分水質級別并賦分[13](表1)。

    表1 水質級別劃分標準及賦分

    1.4 耐污值計算方法

    1.4.1 分類單元和豐富度值的轉化

    依據大型底棲動物分類單元在不同水質級別中出現的個體數將其分為4個類型,并賦予4個不同的豐富度值(表2)。

    表2 個體數和豐富度值間的轉化方法

    1.4.2 初始耐污值(PTV)確定方法

    PTV的確定包含以下幾個關鍵步驟:

    (1)為水質級別賦分,最清潔賦1分;清潔賦2分;輕污賦3分;中污賦4分;重污賦5分(表1);

    (2)以Hydropsycheorientalis為例,按照表2的方法將其在各樣點的個體數轉化為豐富度值;

    (3)計算H.orientalis在各水質級別中的總豐富度,同時計算相應的平均豐富度值如其在最清潔水體中總豐富度值為144,樣本數22,平均豐富度值為144/22=6.55(表3);

    (4)計算H.orientalis在5個水質級別的平均豐富度的累積百分位數(表3)

    1)累加H.orientalis的5個水質級別的平均豐富度值得總平均豐富度值

    6.55+5.05+2.61+1.40+0.73=16.34

    2)計算H.orientalis在每個水質級別中的累積百分位數

    最清潔

    6.55/16.34=40%

    清潔

    (6.55+5.05)/16.34=45%

    輕污

    (6.55+5.05+2.61)/16.34=87%

    中污

    (6.55+5.05+2.61+1.40)/16.34=96%

    重污

    (6.55+5.05+2.61+1.40+0.73)/16.34=100%

    統(tǒng)計H.orientalis在所有樣點出現的次數為98。

    表3 H. orientalis在不同水質級別下的平均豐富度值和累積百分位數

    參考Lenat[5]的方法,確定計算H.orientalis初始耐污值的累積百分位數具體原則如下:分類單元被采集到的次數大于10,以75%分位數值作為初始耐污值; 5—10次,以50%分位數值為初始耐污值;少于5次,不參加統(tǒng)計,結合文獻資料,采用專家觀點的方法確定耐污值。其計算公式如下:

    PTV=L1+(A-N1)/(N2-N1)×(L2-L1)

    式中,L1為上個水質級別的分值;L2為下個水質級別的分值;N1為L1所對應的百分位數;N2為L2所對應的百分位數。A為75%或50%。

    依據以上原則,H.orientalis的出現次數為98,應采用75%分位數對應值作為初始耐污值,從表3可知,75%分位數位于清潔(45%)和輕污(87%)之間,清潔的分值為2,輕污的分值為3,將L1=2,L2=3,N1=45%,N2=87%,代入上式計算其初始耐污值PTV為2.3。

    1.4.3 耐污值(TV)的確定

    采用上述方法確定的初始耐污值的分布范圍一般位于1—4.5。根據公式:TV=2×(1.43×PTV-1.43)將初始耐污值轉換為目前通用的耐污值(分布范圍為0—10)。通過上述計算獲得H.orientalis的TV為3.6。

    1.5 耐污特征的劃分

    根據大型底棲動物耐污值的高低,可將其分為3類[14]:TV≤3,敏感類群(Intolerant group);3—7,一般耐污類群(Intermediate group);TV≥7,耐污類群(Tolerant group)。

    2 結果

    2.1 大型底棲動物群落結構

    在遼河流域共調查采集大型底棲動物236屬(種),共計403039頭,隸屬于4門12綱31目102科。其中節(jié)肢動物門(Arthropoda)8目178種,占底棲動物總數的60%,軟體動物門(Mollusca)5目32種,占11%, 環(huán)節(jié)動物門(Annelida)4目25種,占8%,扁形動物門(Platyhelminthes)1目1種,占2%。

    2.2 水質級別劃分

    利用修訂過的Simpson指數(D′)對遼河流域308個樣點進行水質級別劃分,其中水質級別屬于污染較輕的樣點(最清潔、清潔、輕污)占54.5%,污染較重點位(中污、重污)樣點占到45.5%,樣點比例表明適合耐污值的計算(表3)。

    2.3 大型底棲動物主要分類單元耐污值

    本研究獲得了遼河流域195個大型底棲動物分類單元的耐污值(附錄I),包括昆蟲綱廣翅目1科1屬,蜻蜓目4科5屬,襀翅目3科5屬,毛翅目10科13屬,蜉蝣目8科17屬,鞘翅目6科5屬,半翅目2科1屬,雙翅目12科3亞科18屬,彈尾目1科,端足目1科1屬,十足目1科1屬,基眼目1科3屬,中腹足目4科5屬,蚌目1科1屬,顫蚓目3科4屬,蛭綱3科5屬,蛛形綱1科1屬,渦蟲綱1科1屬。采用專家觀點的方法給出了71科、屬級的耐污值。

    遼河流域耐污值較高的類群有顫蚓科Tubificidae(5—9.9)和石蛭科Erpobdellidae(7.8—8.6)。敏感類群有齒蛉科Corydalidae(2.1),襀科Perlidae(1.3—2.1),原石蛾科Rhyacophilidae(0.6—3),細裳蜉科Leptophlebiidae(0.4—1.7)等,一般耐污類群有紋石蛾科Hydropsychidae(3.7—5.3),四節(jié)蜉科Baetidae(3.4—6.2),扁蜷螺科Planorbidae(3.7—5.5)等。大型底棲動物以敏感種和中度耐污種為主,兩者占到總分類單元個數的84%,而耐污類群種類較少,僅占16%。

    3 討論

    選取適合的方法和標準對樣點水質級別進行劃分是準確計算大型底棲動物耐污值的前提。Lenat[5]采用EPT分類單元數作為水質級別的劃分標準,其將EPT分類單元數大于41的點位水質級別劃分為清潔。本研究發(fā)現遼河流域所有樣點EPT值相對較低,即便水質條件和棲境狀況最好的點位,其EPT分類單元數也不超過27。王備新等[9]報道了我國東部地區(qū)祁門歷溪雙河口、七陽坑和唐云里的EPT值分別為18、14和17;本研究組在廣西紅水河上游長洲河的EPT值為16(未發(fā)表數據)。因此,不難推測美國與中國本底EPT物種豐富度的差異是導致Lenat的水質級別劃分標準不適用于本研究區(qū)域的原因之一。此外,美國水質級別劃分標準是以不同水生態(tài)區(qū)域為基礎,這需要大量樣本數支持,但目前較難實現大樣本數據的獲取,利用EPT進行水質級別劃分對于小樣本數據并不適用。因此,在中國開展長期的水生生物監(jiān)測可更準確地開展水質生物評價,也是與國際已有研究成果進行比較的基礎。

    在對我國東部河流底棲動物耐污值的研究中,王備新等[9]改用多樣性指數進行水質級別劃分,并將以往的4級標準改為5級,以便于與美國分類單元耐污值進行對比。同樣,本研究將Simpson多樣性指數4級劃分標準改為5級,對樣點進行水質級別劃分,所計算遼河流域大型底棲動物耐污值比較準確可信。但本研究所采用的方法仍存在不足,例如只采用Simpson多樣性指數進行樣點的水質級別劃分,而沒有考慮到河流生態(tài)特征的縱向差異性,用一套標準來劃分源頭區(qū)與下游區(qū)的水質等級并不合理。美國對此是按照水生態(tài)區(qū)為基礎,制定了不同生態(tài)區(qū)的河流水質級別劃分標準[5]。因此,中國亟需建立基于水生態(tài)區(qū)的河流管理模式,進而建立合理的基于生物學的河流水質級別劃分標準。

    耐污值作為反映生物耐受外界環(huán)境干擾的一種生物屬性,有一定的穩(wěn)定性,它是與外界環(huán)境特征長期相互作用的結果。王備新等[9]研究發(fā)現生活在不同地理區(qū)的同一生物類群,假如其環(huán)境特征差異較大,會導致其耐污值發(fā)生變化。Lenat在對美國北卡州大型底棲動物耐污值的研究中,考慮到了不同類型河流中大型底棲動物在分類單元組成及耐污性方面的差異,分別計算了不同類型河流大型底棲動物的耐污值[5]。在此基礎上,美國EPA也建立了5個不同地區(qū)的大型底棲動物分類單元耐污值體系[6]。Bressler等[1]在對密西西比河大型底棲動物耐污值的研究中也提到了地區(qū)差異性的問題,底質組成、水文情勢等都可能是引起大型底棲動物耐污值變化的因素。例如,某物種分別出現在沙質底的河流以及受人為干擾后變?yōu)樯迟|底的河流,后者由于存在人類干擾,河流的健康狀況有所下降,物種耐受能力相對升高,因此該物種在前者的耐污值就會低于后者。本研究所計算的耐污值與我國東部地區(qū)間存在一定差異,這與南北方地貌、氣候、水文條件密切相關。比如某些冷水性物種,在南方僅分布在高山溪流區(qū),而在北方地區(qū)分布較廣,從而表現出耐受性相對較強。例如直突搖蚊亞科大多數物種屬冷水性[15],在遼河流域TV是6.9,而在我國東部地區(qū)TV是4.0。底棲動物的分布有一定的地域性,同一科屬所包含的種類不同,每個種類對污染的耐受能力不同,導致不同地區(qū)同一科屬級分類單元耐污值存在差異;即便是不同地理區(qū)域的同一物種,其耐污能力可能隨環(huán)境不同而發(fā)生變化,因此耐污值的計算和應用應具有區(qū)域性的特點[1,5,9]。本研究中大型底棲動物的耐污值均是根據2009年和2010年兩年在遼河流域所獲數據而計算的,故推薦于我國東北地區(qū)使用。

    目前計算耐污值的方法日趨完善,除Lenat[5]的計算方法外,Bressler[1]還采用了PCA對密西西比河大型底棲動物耐污值進行了研究,并與之前北卡州Lenat的計算結果進行了比較,最終修訂了部分物種的耐污值。同時,美國EPA探討了目前國際上普遍使用的PCA[1]、EPT[5]、GAM[16]以及預測模型[17]等方法的適用性,提出使用這些方法的前提要有大量的環(huán)境和生物數據作為基礎,且能夠表征該地區(qū)河流的干擾梯度。此外,PCA、GAM和預測模型法在計算大型底棲動物耐污值時準確性更強,因此這些方法在我國未來大型底棲動物耐污值研究中應予以重視。

    [1] Bressler D W, Stribling J B, Paul M J, Hicks M B. Stressor tolerance values for benthic macroinvertebrates in Mississippi. Hydrobiologia, 2006, 573(1): 155-172.

    [2] Chutter F M. An empirical biotic index of the quality of water in South African streams and rivers. Water Research, 1972, 6(1): 19-30.

    [3] Winget R N, Mangum F A. Biotic condition index-integrated biological, physical, and chemical stream parameters for management: Ogden Utah. U. S. Forest Service, Intermountain Region, 1979: 50-50.

    [4] Hilsenhoff W L. An improved biotic index of organic stream pollution. Great Lakes Entomologist, 1987, 20(1): 31-39.

    [5] Lenat D R. A biotic index for the southeastern United States: derivation and list of tolerance values, with criteria for assigning water-quality rating. Journal of the North American Benthological Society, 1993, 12(3): 279-290.

    [6] Barbour M T, Gerritsen J, Snyder B D, Stribling J B. Rapid Bioassessment Protocols for Use in Streams and Wadeable Rivers: Periphyton, Benthic Macroinvertebrates and Fish. Second Edition. Washingtion, D.C.: U. S. Environmental Protection Agency, Office of Water, EPA841-B-9-9002, 1999: B-1-B-50.

    [7] Hilsenhoff W L. Rapid field assessment of organic pollution with a family-level biotic index. Journal of the North American Benthological Society, 1988, 7(1): 65-68.

    [8] 王建國, 黃恢柏, 楊明旭, 唐振華, 趙鳳霞. 廬山地區(qū)底棲大型無脊椎動物耐污值與水質生物學評價. 應用與環(huán)境生物學報, 2003, 9(3): 279-284.

    [9] 王備新, 楊蓮芳. 我國東部底棲無脊椎動物主要分類單元耐污值. 生態(tài)學報, 2004, 24(12): 2768-2775.

    [10] 邢樹威, 王俊才, 丁振軍, 姜永偉. 遼寧省大型底棲無脊椎動物耐污值及水質評價. 環(huán)境保護科學, 2013, 39(3): 29-33.

    [11] Gao X, Niu C J, Chen Y S, Yin X W. 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. Environmental Monitoring and Assessment, 2013, 186(4): 2375-2391.

    [12] 黃玉瑤, 滕德興,趙忠憲. 應用大型無脊椎動物群落結構特征及其多樣性指數監(jiān)測薊運河污染. 動物學集刊, 1982, 2: 133-144.

    [13] Washington H G. Diversity, biotic, similarity indices: A review with special relevance to aquatic ecosystems. Water Research, 1984, 18(6): 653-694.

    [14] Maxted J R, Barbour M T, Gerritsen J, Poretti V, Primrose N, Silvia A, Penrose D, Renfrow, R. Assessment framework for mid-Atlantic coastal plain streams using benthic macroinvertebrates. Journal of the North American Benthological Society, 2000, 19(1): 128-144.

    [15] Armitage P D, Cranston P S, Pinde L C V. The Chironomidae: Biology and Ecology of Non-Biting Midges. London: Chapman & Hall, 1995: 62-68.

    [16] U. S. Environmental Protection Agency. The evaluation of methods for creating defensible, repeatable, objective and accurate tolerance values for aquatic Taxa. U S Environmental Protection Agency Office of Research and Accurate Tolerance Values for Aquatic Taxa. Washington, DC: EPA, 2006:1-50.

    [17] Hawkins C P, Norris R H, Hogue J N, Feminella J W. Development and evaluation of predictive models for measuring the biological integrity of streams. Ecological Applications, 2000, 10(5): 1456-1477.

    附錄 I 遼河流域大型底棲無脊椎動物分類單元耐污值與同一分類單元在中國東部地區(qū)、美國 5 個地區(qū)耐污值的比較( LR,遼河;EC,中國東部地區(qū);NC,北卡州;WI,威斯康星州;OH,俄亥俄州;DI,愛達荷州;MACS,大西洋海岸中部溪流工作組)

    Appendix I Tolerance values of macroinvertebrate taxa in Liao River, and the value of the same taxon in the Eastern area of China and NC(North Carolina, Southeast, WI (Wisconsin,Upper Midwest), OH (Ohio, Midwest), ID(Idaho DEP, Northwest), and MACS (Mid-Atlantic Coastal Streams

    Tolerance values of macroinvertebrate taxa in Liao River basin

    ZHAO Rui1,2, GAO Xin1,2, DING Sen1,2, ZHANG Yuan1,2,*, QU Xiaodong3, LIU Sisi1,2

    1StateKeyLaboratoryofEnvironmentalCriteriaandRiskAssessment,ChineseResearchAcademyofEnvironmentalSciences,Beijing100012,China2LaboratoryofRiverineEcologicalConservationandTechnology,ChineseResearchAcademyofEnvironmentalSciences,Beijing100012,China3DepartmentofWaterEnvironment,ChinaInstituteofWaterResourcesandHydropowerResearch,Beijing100038,China

    Tolerance represents the relative capacity of aquatic organisms to survive and reproduce under different levels of environmental stress. Tolerance values (TVs) of macroinvertebrates are broadly used for calculating metrics applied in bioassessment. Variability among TVs may be related to regional differences in macroinvertebrate sensitivity to pollution; whether the calculations and amendments of TVs are appropriate for one region is the basis for accurate evaluation of river health condition. Two methods, namely, expert opinion and statistical analysis, have been used to derive the TVs of macroinvertebrate taxa. The TVs currently used in bioassessment of northern rivers in China were acquired from relevant literature and expert opinion. The bioassessment results, however, will be more accurate if the TVs are derived from statistical analysis. The objective of the current study was to derive the TVs appropriate for macroinvertebrate taxa in northeastern rivers in China. Macroinvertebrate data from 308 samples were collected in spring and autumn in 2009 and summer in 2010 in the Liao River basin. Using the revised Simpson index, we assigned the water quality into five categories as: ≥ 10, excellent; ≥ 6, good; ≥ 3, good-fair; ≥ 2, fair; and < 2, poor. Different scores were marked to different water quality: excellent, 1; good, 2; good-fair, 3; fair, 4 and poor, 5. Macroinvertebrates were grouped and coded based on abundance at each site: < 1 ind./sample = 0; 1—2 ind./sample = 1; 3—9 ind./sample = 3; and ≥ 10 ind./sample = 10. The average abundance values for each taxon were calculated for each of the five water quality classes. The average abundance values were converted into a cumulative percentile: the 75thpercentile value was used for a taxon that had > 10 records and the 50thpercentile value was used for a taxon that had 5—10 records. This yielded preliminary tolerance values of 0.0—4.5, which were transformed to final tolerance values on a 0—10 scale. A total of 195 macroinvertebrate taxa TVs were calculated, including one genus in one family for each Megaloptera, Amphipoda, Decapoda, Basommatophora, Unionoida, Tubificida, Hirudinea, Turbellaria, and Arachnida; one family of Collembola; two families (including one taxon at genus-level) of Hemiptera; five genera in three families of Plecoptera; five genera in four families of Odonata; five genera in four families of Mesogastropoda; six families (including five taxa at genus-level) of Coleoptera; 13 genera in 10 families of Trichoptera; 17 genera in eight families of Ephemeroptera; and 18 genera and three sub-families in 12 families of Diptera. The other 71 TVs were derived from expert opinion and relevant literature. TVs were classified according to sensitivity as: TV ≤ 3, intolerant; 3 < TV < 7, intermediate; and TV ≥ 7, tolerant. The dominant macroinvertebrate families in the Liao River basin were divided into the three groups as follows: intolerant group, including Corydalidae (2.1), Perlidae (1.3—2.1), Rhyacophilidae (0.6—3.0), and Leptophlebiidae (0.4—1.7); intermediate group, including Hydropsychidae (3.7—5.3), Baetidae (3.4—6.2), and Planorbidae (3.7—5.5); and tolerant group, including Tubificidae (5.0—9.9) and Erpobdellidae (7.8—8.6). The intolerant and intermediate groups comprised nearly 84% of the total taxa, while the tolerant group comprised only 16% of the total taxa. We suggest that TVs can be used to assess the water quality of streams and rivers in northern China.

    macroinvertebrates; tolerance value; Liao River basin

    分類單元TaxaLRECNCWIIDOHMACS節(jié)肢動物門Arthropoda昆蟲綱Insecta廣翅目Megaloptera齒蛉科Corydalidae2.13.80黃石蛉Protohermesgrandis2.1星齒嶺屬Protohermesspp.2.13.8泥蛉科Sialidae4.5*4泥蛉屬Sialisspp.4.5*蜻蜓目Odonata蟌科Coenagrionidae7.4996.19蔚藍細蟌Cercionsexlineatum7.5*Enllagmacivile5*春蜓科Gomphidae1.52.71施春蜓屬Sieboldiusspp.3.5*小葉春蜓屬Gomphidiaspp.5.2新月戴春蜓Davidiuslunatus2.2擴腹春蜓屬Stylurusspp.1長唇春蜓Anisogomphusmaacki0.7葉春蜓屬Ictinogomphusspp.3.5*蜻科Libellulidae9.48.5色蟌科Claopterygidae0.72.4日本色蟌Calopteryxjaponica0.7扇蟌科Platicnymidae4.5*9扇蟌屬platycnemisspp.4.5*襀翅目Plecoptera襀科Perlidae1.81.211石蠅屬Neoperlaspp.2.11.613.1大山石蠅屬Oyamiaspp.1.3網襀科Perlodidae1.42石蠅屬Isoperlaspp.0Megrcysochracea1.9Eccopturaspp.3.5*綠襀科Chloroperlidae211蘇瓦襀屬Suwalliaspp.201毛翅目Trichoptera紋石蛾科Hydropsychidae53.744短脈紋石蛾屬Cheumatopsychespp.4.53.86.5552.95東方側枝紋石蛾Hydropsycheorientalis3.3缺尾高原紋石蛾Hydropsychekozha-ntschikovi4.6紋石蛾Hydropsycheyaeyamensis3.1紋石蛾屬Hydropsychespp.4.344旋刺枝紋石蛾Hydropsychenevae3.7Arctopsychespp.41角石蛾科Stenopsychidae1.74.7條紋角石蛾Stenopsychemarmorata1.7角石蛾屬Stenopsychespp.1.74.7原石蛾科Rhyacophilidae1.810原石蛾屬Rhyacophilaspp.1.80

    國家水體污染控制與治理科技重大專項(2012ZX07501-001-04); 中歐環(huán)境可持續(xù)發(fā)展計劃(DCI-ASIE/2013/323-261); 國家自然科學基金(2010KYYW06); 國家自然科學基金(41401066)

    2013-11-08;

    2014-09-09

    10.5846/stxb201311082699

    *通訊作者Corresponding author.E-mail: zhangyuan@craes.org.cn

    趙瑞, 高欣, 丁森, 張遠, 渠曉東, 劉思思.遼河流域大型底棲動物耐污值.生態(tài)學報,2015,35(14):4797-4809.

    Zhao R, Gao X, Ding S, Zhang Y, Qu X D, Liu S S. Tolerance values of macroinvertebrate taxa in Liao River basin.Acta Ecologica Sinica,2015,35(14):4797-4809.

    猜你喜歡
    遼河流域樣點級別
    小麥條銹病田間為害損失的初步分析
    湖北植保(2022年4期)2022-08-23 10:51:52
    痘痘分級別,輕重不一樣
    基于空間模擬退火算法的最優(yōu)土壤采樣尺度選擇研究①
    土壤(2021年1期)2021-03-23 07:29:06
    邁向UHD HDR的“水晶” 十萬元級別的SIM2 CRYSTAL4 UHD
    新年導購手冊之兩萬元以下級別好物推薦
    基于分融策略的土壤采樣設計方法*
    土壤學報(2017年5期)2017-11-01 09:21:27
    你是什么級別的
    中學生(2017年13期)2017-06-15 12:57:48
    遼寧省遼河流域石油煉制排放限值的制定
    遼河流域水生態(tài)管理指標體系構建
    遼河流域水文特性分析
    久久中文字幕一级| 另类亚洲欧美激情| 免费久久久久久久精品成人欧美视频| 久久av网站| 中文字幕亚洲精品专区| 一级毛片我不卡| 亚洲国产精品一区三区| 99热国产这里只有精品6| 人体艺术视频欧美日本| 国产又色又爽无遮挡免| 欧美+亚洲+日韩+国产| 亚洲精品日本国产第一区| 丝袜脚勾引网站| 欧美久久黑人一区二区| 啦啦啦在线观看免费高清www| 狂野欧美激情性xxxx| 男女国产视频网站| 天天躁日日躁夜夜躁夜夜| 日本欧美国产在线视频| 午夜两性在线视频| 十八禁高潮呻吟视频| 少妇精品久久久久久久| 高清视频免费观看一区二区| 夫妻性生交免费视频一级片| 一本一本久久a久久精品综合妖精| 性高湖久久久久久久久免费观看| 大话2 男鬼变身卡| 免费在线观看黄色视频的| 日韩av不卡免费在线播放| 日韩av不卡免费在线播放| 中文字幕人妻熟女乱码| 久久久久国产一级毛片高清牌| 国产亚洲午夜精品一区二区久久| 亚洲色图综合在线观看| 五月天丁香电影| 国产精品熟女久久久久浪| 免费日韩欧美在线观看| 另类精品久久| 1024香蕉在线观看| 中文字幕人妻丝袜制服| 午夜av观看不卡| 黄色片一级片一级黄色片| 国产亚洲一区二区精品| 国产熟女欧美一区二区| 午夜免费男女啪啪视频观看| 无遮挡黄片免费观看| 久久99热这里只频精品6学生| 国产国语露脸激情在线看| 十八禁网站网址无遮挡| 国产真人三级小视频在线观看| 国产精品免费视频内射| 精品国产国语对白av| 男女边吃奶边做爰视频| 首页视频小说图片口味搜索 | 免费女性裸体啪啪无遮挡网站| 亚洲久久久国产精品| 亚洲视频免费观看视频| 成在线人永久免费视频| 18禁裸乳无遮挡动漫免费视频| 好男人视频免费观看在线| 中文字幕制服av| 国产欧美日韩精品亚洲av| 国产黄频视频在线观看| 亚洲色图综合在线观看| 亚洲专区中文字幕在线| 精品国产一区二区三区久久久樱花| 免费在线观看黄色视频的| 校园人妻丝袜中文字幕| 欧美日韩亚洲高清精品| 亚洲成人免费av在线播放| 黄色一级大片看看| 19禁男女啪啪无遮挡网站| 亚洲国产精品一区三区| 久久精品久久久久久噜噜老黄| 午夜福利一区二区在线看| 蜜桃国产av成人99| 9191精品国产免费久久| 欧美大码av| 啦啦啦在线观看免费高清www| 亚洲熟女精品中文字幕| 青草久久国产| 男女免费视频国产| 免费看不卡的av| 亚洲人成网站在线观看播放| 欧美亚洲日本最大视频资源| 丰满饥渴人妻一区二区三| 欧美日韩成人在线一区二区| 国产男女超爽视频在线观看| 国产在线免费精品| 巨乳人妻的诱惑在线观看| 男人舔女人的私密视频| 欧美日韩av久久| 一区二区三区乱码不卡18| 大片电影免费在线观看免费| 七月丁香在线播放| 日本av免费视频播放| 国产免费福利视频在线观看| 永久免费av网站大全| 国产亚洲欧美在线一区二区| 欧美激情极品国产一区二区三区| 91精品国产国语对白视频| 两个人看的免费小视频| 高清av免费在线| 国产精品免费大片| 国产精品国产三级专区第一集| 亚洲欧美日韩高清在线视频 | 久久亚洲精品不卡| 日本欧美国产在线视频| 国产一区二区 视频在线| 国产精品三级大全| 国产国语露脸激情在线看| 亚洲熟女精品中文字幕| 只有这里有精品99| 国产高清不卡午夜福利| 成人免费观看视频高清| 另类亚洲欧美激情| 在线观看免费视频网站a站| 亚洲国产精品国产精品| 纯流量卡能插随身wifi吗| 晚上一个人看的免费电影| 亚洲中文日韩欧美视频| 又大又爽又粗| 乱人伦中国视频| 美女扒开内裤让男人捅视频| 欧美日韩亚洲国产一区二区在线观看 | 日韩制服骚丝袜av| 日韩人妻精品一区2区三区| 麻豆av在线久日| 波多野结衣一区麻豆| 最黄视频免费看| 国产福利在线免费观看视频| 国产精品久久久久久人妻精品电影 | 免费女性裸体啪啪无遮挡网站| 国产欧美亚洲国产| 亚洲欧美精品自产自拍| 亚洲av男天堂| av电影中文网址| 天天躁夜夜躁狠狠躁躁| 国产一区二区激情短视频 | 这个男人来自地球电影免费观看| 精品一区二区三区四区五区乱码 | 国产高清视频在线播放一区 | 国产伦理片在线播放av一区| 久久免费观看电影| 国产人伦9x9x在线观看| 精品福利观看| 国产精品久久久人人做人人爽| 交换朋友夫妻互换小说| 一边亲一边摸免费视频| 一本大道久久a久久精品| 久久人人97超碰香蕉20202| 日韩,欧美,国产一区二区三区| 大香蕉久久成人网| 成人影院久久| 亚洲五月色婷婷综合| 日本91视频免费播放| 国产熟女欧美一区二区| 国产成人91sexporn| 亚洲精品久久午夜乱码| 一本—道久久a久久精品蜜桃钙片| 女人精品久久久久毛片| 国产99久久九九免费精品| 一边摸一边抽搐一进一出视频| 亚洲av成人不卡在线观看播放网 | 亚洲精品乱久久久久久| 脱女人内裤的视频| 国产亚洲一区二区精品| 国产精品久久久人人做人人爽| 国产福利在线免费观看视频| 最黄视频免费看| 精品一区在线观看国产| 亚洲欧美中文字幕日韩二区| videos熟女内射| 大香蕉久久网| 中文欧美无线码| 制服人妻中文乱码| 国产精品99久久99久久久不卡| 国产野战对白在线观看| 美女高潮到喷水免费观看| 国产色视频综合| 国产日韩一区二区三区精品不卡| 欧美日韩av久久| 国产xxxxx性猛交| 免费观看av网站的网址| 女人精品久久久久毛片| 亚洲国产最新在线播放| 成人国产一区最新在线观看 | 狠狠精品人妻久久久久久综合| 国产又色又爽无遮挡免| 性高湖久久久久久久久免费观看| 欧美精品一区二区免费开放| 搡老乐熟女国产| 97人妻天天添夜夜摸| 在线观看免费视频网站a站| 亚洲av日韩在线播放| 日本午夜av视频| 19禁男女啪啪无遮挡网站| 视频在线观看一区二区三区| 欧美黄色淫秽网站| 50天的宝宝边吃奶边哭怎么回事| 午夜日韩欧美国产| 亚洲欧美日韩另类电影网站| 制服人妻中文乱码| 精品福利观看| 国产免费又黄又爽又色| 亚洲人成77777在线视频| 精品熟女少妇八av免费久了| av在线播放精品| 丝袜在线中文字幕| 久久毛片免费看一区二区三区| 99国产精品一区二区蜜桃av | 久久精品人人爽人人爽视色| 麻豆av在线久日| 亚洲情色 制服丝袜| 国产在线一区二区三区精| 亚洲国产精品一区二区三区在线| 亚洲精品久久午夜乱码| 国产精品国产av在线观看| 国产国语露脸激情在线看| 国产精品一区二区精品视频观看| 精品国产一区二区三区久久久樱花| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一区蜜桃| 国产一区二区三区av在线| 2021少妇久久久久久久久久久| 午夜福利,免费看| 亚洲人成电影免费在线| 精品少妇一区二区三区视频日本电影| 欧美 亚洲 国产 日韩一| 国产在线一区二区三区精| 汤姆久久久久久久影院中文字幕| 中文乱码字字幕精品一区二区三区| 看免费成人av毛片| 亚洲欧洲精品一区二区精品久久久| 日本五十路高清| 丝袜人妻中文字幕| 欧美国产精品va在线观看不卡| 久久鲁丝午夜福利片| 中文精品一卡2卡3卡4更新| 黄色 视频免费看| 日本色播在线视频| 亚洲精品一卡2卡三卡4卡5卡 | cao死你这个sao货| av国产久精品久网站免费入址| 久久久欧美国产精品| 十分钟在线观看高清视频www| 黄网站色视频无遮挡免费观看| 丝袜美腿诱惑在线| 丝袜美足系列| 99re6热这里在线精品视频| 成人三级做爰电影| 亚洲精品一区蜜桃| 国产一区二区三区av在线| www日本在线高清视频| 久久久久精品人妻al黑| 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费av在线播放| 波多野结衣一区麻豆| 老司机在亚洲福利影院| 老司机影院毛片| 精品一区二区三区av网在线观看 | 亚洲国产最新在线播放| 女性生殖器流出的白浆| 欧美精品av麻豆av| 国产又爽黄色视频| 老司机午夜十八禁免费视频| 999精品在线视频| 大型av网站在线播放| 一级毛片女人18水好多 | 日韩制服骚丝袜av| 91精品伊人久久大香线蕉| 一级黄色大片毛片| 国产免费又黄又爽又色| 国产av精品麻豆| 国产亚洲欧美在线一区二区| 操出白浆在线播放| 亚洲精品av麻豆狂野| 久久精品国产a三级三级三级| 熟女av电影| 纵有疾风起免费观看全集完整版| 中文精品一卡2卡3卡4更新| 午夜免费鲁丝| 黄网站色视频无遮挡免费观看| 人人妻人人爽人人添夜夜欢视频| 亚洲色图 男人天堂 中文字幕| 啦啦啦视频在线资源免费观看| 老熟女久久久| 大话2 男鬼变身卡| 成人影院久久| 99久久人妻综合| 亚洲国产精品一区二区三区在线| 少妇猛男粗大的猛烈进出视频| 国产一卡二卡三卡精品| 久久久久精品国产欧美久久久 | 久久久亚洲精品成人影院| 国产免费一区二区三区四区乱码| 搡老乐熟女国产| 啦啦啦中文免费视频观看日本| 晚上一个人看的免费电影| 岛国毛片在线播放| 久久这里只有精品19| 最近手机中文字幕大全| 亚洲精品久久成人aⅴ小说| 又大又爽又粗| 日本av手机在线免费观看| 乱人伦中国视频| 99国产精品一区二区蜜桃av | 91精品伊人久久大香线蕉| 搡老岳熟女国产| 日韩大码丰满熟妇| 久久国产亚洲av麻豆专区| tube8黄色片| 十八禁人妻一区二区| 亚洲成人手机| 国产主播在线观看一区二区 | 精品亚洲乱码少妇综合久久| 99香蕉大伊视频| 国产亚洲一区二区精品| 在线观看免费日韩欧美大片| 亚洲精品成人av观看孕妇| 交换朋友夫妻互换小说| 欧美日韩av久久| 午夜福利在线免费观看网站| 美女扒开内裤让男人捅视频| svipshipincom国产片| 国产高清videossex| 日韩 亚洲 欧美在线| 美女中出高潮动态图| 男女床上黄色一级片免费看| 久久精品亚洲av国产电影网| 午夜福利视频精品| av在线老鸭窝| 免费高清在线观看视频在线观看| 亚洲成人手机| 成年动漫av网址| 一区二区三区精品91| 美女福利国产在线| 国产又爽黄色视频| 成人手机av| 好男人视频免费观看在线| 日韩 欧美 亚洲 中文字幕| 亚洲色图 男人天堂 中文字幕| av电影中文网址| 亚洲,欧美精品.| 十分钟在线观看高清视频www| 一级黄色大片毛片| 国产黄色免费在线视频| 久久久久网色| 国产成人精品无人区| av电影中文网址| 欧美日韩综合久久久久久| 久久 成人 亚洲| www.999成人在线观看| 欧美精品av麻豆av| 国产午夜精品一二区理论片| 精品一区在线观看国产| 亚洲第一青青草原| 亚洲国产精品一区三区| 中文字幕人妻丝袜一区二区| av国产精品久久久久影院| 在线精品无人区一区二区三| 黄色怎么调成土黄色| 亚洲视频免费观看视频| 久久综合国产亚洲精品| 欧美人与善性xxx| av天堂久久9| 国产欧美日韩一区二区三 | 两人在一起打扑克的视频| av天堂在线播放| 国产成人91sexporn| 国产精品久久久久久精品古装| 亚洲精品久久久久久婷婷小说| 欧美性长视频在线观看| 国产成人精品久久二区二区91| www日本在线高清视频| 亚洲国产精品一区二区三区在线| 午夜激情av网站| 母亲3免费完整高清在线观看| 人妻人人澡人人爽人人| 男女之事视频高清在线观看 | 国产成人精品在线电影| 美国免费a级毛片| 国产av精品麻豆| 日韩免费高清中文字幕av| 色94色欧美一区二区| 啦啦啦中文免费视频观看日本| 久久鲁丝午夜福利片| 日日爽夜夜爽网站| 国产高清不卡午夜福利| 国产在视频线精品| 99国产精品一区二区蜜桃av | 观看av在线不卡| 中文字幕制服av| 女警被强在线播放| 国产野战对白在线观看| 一级片'在线观看视频| 欧美97在线视频| 婷婷丁香在线五月| 亚洲人成网站在线观看播放| 久久久久精品国产欧美久久久 | 少妇精品久久久久久久| 欧美人与性动交α欧美精品济南到| 国产黄色视频一区二区在线观看| 视频区欧美日本亚洲| 午夜激情久久久久久久| 国产老妇伦熟女老妇高清| 亚洲精品一二三| 十分钟在线观看高清视频www| 热re99久久精品国产66热6| 制服诱惑二区| 不卡av一区二区三区| 欧美日本中文国产一区发布| www.熟女人妻精品国产| 女人高潮潮喷娇喘18禁视频| 在线观看免费日韩欧美大片| 美女国产高潮福利片在线看| videosex国产| 两个人看的免费小视频| 日韩精品免费视频一区二区三区| 99热国产这里只有精品6| 亚洲欧美一区二区三区国产| 男女之事视频高清在线观看 | 丝瓜视频免费看黄片| 国产精品偷伦视频观看了| 老司机亚洲免费影院| 精品国产一区二区三区久久久樱花| 50天的宝宝边吃奶边哭怎么回事| 一区福利在线观看| 桃花免费在线播放| 纯流量卡能插随身wifi吗| 视频在线观看一区二区三区| a级毛片黄视频| 亚洲国产成人一精品久久久| 一级毛片我不卡| 91麻豆精品激情在线观看国产 | 赤兔流量卡办理| av又黄又爽大尺度在线免费看| www.自偷自拍.com| 国产一区二区三区av在线| 久久久久国产精品人妻一区二区| 黄频高清免费视频| 大片免费播放器 马上看| 国产精品二区激情视频| 亚洲国产中文字幕在线视频| 国产亚洲av片在线观看秒播厂| 免费不卡黄色视频| 男女无遮挡免费网站观看| 中文字幕最新亚洲高清| 国产精品.久久久| 国产在视频线精品| 麻豆av在线久日| 老汉色∧v一级毛片| 亚洲国产av新网站| 青春草视频在线免费观看| www.熟女人妻精品国产| 欧美激情极品国产一区二区三区| 我要看黄色一级片免费的| 久久免费观看电影| 午夜视频精品福利| 狂野欧美激情性xxxx| 亚洲精品久久成人aⅴ小说| 精品第一国产精品| 女性生殖器流出的白浆| 女性被躁到高潮视频| 十八禁人妻一区二区| 国产欧美日韩一区二区三区在线| 成年女人毛片免费观看观看9 | 777米奇影视久久| 国产av一区二区精品久久| 高潮久久久久久久久久久不卡| 精品久久久久久久毛片微露脸 | 亚洲精品中文字幕在线视频| 桃花免费在线播放| 国产精品久久久人人做人人爽| xxxhd国产人妻xxx| 人体艺术视频欧美日本| 国产精品一区二区精品视频观看| 亚洲天堂av无毛| 久久精品国产a三级三级三级| 中文字幕色久视频| kizo精华| 丝袜美足系列| 人成视频在线观看免费观看| 亚洲精品在线美女| 最新在线观看一区二区三区 | xxx大片免费视频| 欧美老熟妇乱子伦牲交| 国产免费一区二区三区四区乱码| 一级片免费观看大全| 老司机靠b影院| 欧美黄色淫秽网站| 宅男免费午夜| 咕卡用的链子| 一级,二级,三级黄色视频| 少妇粗大呻吟视频| 亚洲精品乱久久久久久| 尾随美女入室| 一级毛片女人18水好多 | 成人三级做爰电影| 黄色毛片三级朝国网站| 亚洲成人国产一区在线观看 | 日韩av免费高清视频| 国产野战对白在线观看| 视频区图区小说| 久久久久久久国产电影| 国产一区二区三区av在线| 亚洲精品在线美女| 久久精品亚洲av国产电影网| 啦啦啦啦在线视频资源| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美精品综合一区二区三区| 国产有黄有色有爽视频| 国产一区二区三区av在线| 美国免费a级毛片| 国产精品久久久久久精品古装| 国产精品成人在线| 久久久久久免费高清国产稀缺| 久久久亚洲精品成人影院| 日韩 亚洲 欧美在线| 国产成人精品在线电影| 男女之事视频高清在线观看 | 欧美乱码精品一区二区三区| 亚洲人成77777在线视频| 丝袜美足系列| 国产男女超爽视频在线观看| 精品国产国语对白av| 大陆偷拍与自拍| 亚洲av片天天在线观看| 国产一区二区激情短视频 | 精品少妇久久久久久888优播| 精品久久久精品久久久| 黄色视频不卡| 美女福利国产在线| 国产在线免费精品| 中文精品一卡2卡3卡4更新| 欧美黄色淫秽网站| 好男人电影高清在线观看| 婷婷色综合大香蕉| 男人舔女人的私密视频| 美女主播在线视频| www.自偷自拍.com| 七月丁香在线播放| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 亚洲国产看品久久| 大片电影免费在线观看免费| 午夜免费观看性视频| 中文字幕最新亚洲高清| 人妻一区二区av| 91老司机精品| 国产在线视频一区二区| 国产成人精品久久二区二区免费| 精品国产一区二区三区久久久樱花| 亚洲精品国产av成人精品| 肉色欧美久久久久久久蜜桃| 天天躁日日躁夜夜躁夜夜| 日韩熟女老妇一区二区性免费视频| 精品国产国语对白av| 久久人妻熟女aⅴ| 下体分泌物呈黄色| 久久这里只有精品19| 18禁国产床啪视频网站| 国产精品久久久久成人av| 国产精品成人在线| av在线老鸭窝| 一区二区三区四区激情视频| 老司机影院毛片| 在现免费观看毛片| 最新的欧美精品一区二区| 中文字幕亚洲精品专区| 精品国产乱码久久久久久小说| 精品亚洲成国产av| 日本一区二区免费在线视频| 丁香六月天网| 国产激情久久老熟女| 男女无遮挡免费网站观看| 亚洲国产精品一区三区| 女性生殖器流出的白浆| 99国产精品一区二区三区| 只有这里有精品99| 天天躁夜夜躁狠狠躁躁| 免费不卡黄色视频| 免费少妇av软件| 性少妇av在线| 日本91视频免费播放| 欧美中文综合在线视频| 蜜桃在线观看..| 99国产精品99久久久久| 中文字幕人妻丝袜制服| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品久久久av美女十八| 国产精品 欧美亚洲| 精品少妇黑人巨大在线播放| 精品少妇久久久久久888优播| 在线精品无人区一区二区三| 91精品伊人久久大香线蕉| 51午夜福利影视在线观看| 久久人人爽人人片av| 中文字幕人妻熟女乱码| 国产野战对白在线观看| 国产精品欧美亚洲77777| 尾随美女入室| 日日摸夜夜添夜夜爱| 美女高潮到喷水免费观看| 一本久久精品| 最近最新中文字幕大全免费视频 | 欧美日韩成人在线一区二区| 亚洲 国产 在线| 国产成人欧美在线观看 | 亚洲精品久久久久久婷婷小说| 中文字幕亚洲精品专区| 可以免费在线观看a视频的电影网站| 久久天躁狠狠躁夜夜2o2o | 50天的宝宝边吃奶边哭怎么回事| 国产一区二区激情短视频 |