• <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)管理指標體系構建
    遼河流域水文特性分析
    亚洲欧美一区二区三区黑人| 嫁个100分男人电影在线观看| 男人舔女人的私密视频| 最近最新中文字幕大全电影3| 日本a在线网址| 欧美黄色片欧美黄色片| av在线蜜桃| 欧美日韩瑟瑟在线播放| 国产精品 欧美亚洲| 69av精品久久久久久| 久久久久久久久久黄片| 午夜视频精品福利| 久久国产精品影院| 一本综合久久免费| 亚洲avbb在线观看| 香蕉国产在线看| 19禁男女啪啪无遮挡网站| 亚洲美女黄片视频| 特级一级黄色大片| 久久久久久久久免费视频了| 黄色片一级片一级黄色片| 一级a爱片免费观看的视频| 搞女人的毛片| 亚洲在线观看片| 国内精品久久久久久久电影| 禁无遮挡网站| 亚洲欧美一区二区三区黑人| 九九热线精品视视频播放| 又粗又爽又猛毛片免费看| 黄色成人免费大全| 99久久99久久久精品蜜桃| 国产精品永久免费网站| 日韩欧美免费精品| 久久天堂一区二区三区四区| 12—13女人毛片做爰片一| 欧美日韩国产亚洲二区| xxx96com| 欧美中文综合在线视频| 午夜激情欧美在线| 精品久久久久久久末码| 国内毛片毛片毛片毛片毛片| 波多野结衣高清作品| 色噜噜av男人的天堂激情| 女人高潮潮喷娇喘18禁视频| 老熟妇仑乱视频hdxx| 国产乱人伦免费视频| 久久久久免费精品人妻一区二区| 国产av一区在线观看免费| 性欧美人与动物交配| 中文字幕精品亚洲无线码一区| 国产亚洲精品久久久com| 色播亚洲综合网| 欧美成人性av电影在线观看| 琪琪午夜伦伦电影理论片6080| 精品国产三级普通话版| 久久久久久人人人人人| 亚洲熟女毛片儿| 在线观看66精品国产| 成人鲁丝片一二三区免费| 欧美色视频一区免费| 国产精品精品国产色婷婷| 亚洲欧美日韩无卡精品| 伊人久久大香线蕉亚洲五| 免费在线观看影片大全网站| 午夜日韩欧美国产| 国产aⅴ精品一区二区三区波| 亚洲国产中文字幕在线视频| 国产麻豆成人av免费视频| 久久午夜亚洲精品久久| 亚洲精华国产精华精| 99国产精品一区二区蜜桃av| 国产精品,欧美在线| 亚洲电影在线观看av| 一本综合久久免费| 9191精品国产免费久久| 看免费av毛片| 国产精品一区二区免费欧美| 老司机午夜十八禁免费视频| 成熟少妇高潮喷水视频| 我要搜黄色片| 国产伦精品一区二区三区视频9 | 他把我摸到了高潮在线观看| av国产免费在线观看| 亚洲真实伦在线观看| 亚洲 欧美 日韩 在线 免费| 岛国在线观看网站| 欧美不卡视频在线免费观看| 又黄又爽又免费观看的视频| 手机成人av网站| h日本视频在线播放| 欧美日韩国产亚洲二区| 精品国产乱子伦一区二区三区| 中文在线观看免费www的网站| 亚洲人成伊人成综合网2020| 国产精品美女特级片免费视频播放器 | 每晚都被弄得嗷嗷叫到高潮| 淫妇啪啪啪对白视频| 久久久久国产一级毛片高清牌| 性色avwww在线观看| 成人特级av手机在线观看| 久久婷婷人人爽人人干人人爱| 欧美+亚洲+日韩+国产| 亚洲人成网站在线播放欧美日韩| 99国产精品一区二区三区| 波多野结衣巨乳人妻| 精品久久久久久久久久免费视频| xxx96com| 在线播放国产精品三级| 最新在线观看一区二区三区| 国产黄色小视频在线观看| 九色成人免费人妻av| 99国产极品粉嫩在线观看| 免费av不卡在线播放| 欧美色欧美亚洲另类二区| 午夜精品在线福利| 夜夜爽天天搞| 日韩欧美免费精品| 亚洲熟妇熟女久久| 他把我摸到了高潮在线观看| 在线观看66精品国产| 国产乱人伦免费视频| 国内精品一区二区在线观看| 免费在线观看日本一区| 91在线观看av| 黑人巨大精品欧美一区二区mp4| 亚洲国产高清在线一区二区三| 精品人妻1区二区| 国产精品香港三级国产av潘金莲| 久久久久久久久久黄片| 国产精品,欧美在线| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久av美女十八| 亚洲精品在线观看二区| 色老头精品视频在线观看| 日日摸夜夜添夜夜添小说| 高潮久久久久久久久久久不卡| 精品久久久久久久久久久久久| 老司机福利观看| 色视频www国产| 国产美女午夜福利| 国产精华一区二区三区| 久久久久久久精品吃奶| 搡老妇女老女人老熟妇| АⅤ资源中文在线天堂| 黄片小视频在线播放| 丰满人妻一区二区三区视频av | 精品无人区乱码1区二区| 亚洲成人久久爱视频| 可以在线观看的亚洲视频| 亚洲成a人片在线一区二区| 色尼玛亚洲综合影院| 国内精品久久久久精免费| 两性午夜刺激爽爽歪歪视频在线观看| 激情在线观看视频在线高清| 亚洲色图av天堂| 久久久久国产精品人妻aⅴ院| 亚洲一区高清亚洲精品| 亚洲精华国产精华精| 成在线人永久免费视频| 19禁男女啪啪无遮挡网站| www日本在线高清视频| 国产伦精品一区二区三区四那| av黄色大香蕉| 夜夜夜夜夜久久久久| 成年人黄色毛片网站| www.www免费av| 中文字幕熟女人妻在线| 国产午夜精品久久久久久| 国产午夜福利久久久久久| 高潮久久久久久久久久久不卡| 精品久久久久久久久久久久久| 精品午夜福利视频在线观看一区| 99久久综合精品五月天人人| 美女高潮的动态| 亚洲成av人片免费观看| a在线观看视频网站| 色精品久久人妻99蜜桃| 日韩欧美免费精品| 欧美黑人欧美精品刺激| 日本五十路高清| 又粗又爽又猛毛片免费看| www.自偷自拍.com| 国产爱豆传媒在线观看| 一进一出抽搐gif免费好疼| 精品国产超薄肉色丝袜足j| 国产精品国产高清国产av| 亚洲第一电影网av| 桃色一区二区三区在线观看| 免费在线观看成人毛片| 国产真实乱freesex| 在线观看日韩欧美| 欧美在线一区亚洲| 久久久成人免费电影| 男女下面进入的视频免费午夜| 欧美日本亚洲视频在线播放| 91av网一区二区| 色在线成人网| 女人高潮潮喷娇喘18禁视频| 天天躁狠狠躁夜夜躁狠狠躁| 美女大奶头视频| 丁香欧美五月| 最好的美女福利视频网| 精品久久久久久,| 欧美日本亚洲视频在线播放| 国内揄拍国产精品人妻在线| 亚洲一区二区三区色噜噜| 国产精品久久久久久久电影 | 免费看a级黄色片| 我的老师免费观看完整版| 亚洲欧美精品综合一区二区三区| 亚洲一区二区三区不卡视频| 免费一级毛片在线播放高清视频| 亚洲欧美日韩高清在线视频| 日日夜夜操网爽| 久久久久免费精品人妻一区二区| 欧美不卡视频在线免费观看| 国产精品一区二区三区四区免费观看 | 人人妻,人人澡人人爽秒播| 男人舔女人的私密视频| 身体一侧抽搐| 每晚都被弄得嗷嗷叫到高潮| 国产极品精品免费视频能看的| 亚洲精品在线美女| 国语自产精品视频在线第100页| 国产亚洲精品一区二区www| 欧美精品啪啪一区二区三区| 1000部很黄的大片| 日本黄大片高清| 精品人妻1区二区| 熟妇人妻久久中文字幕3abv| 757午夜福利合集在线观看| 国产视频一区二区在线看| 亚洲av成人一区二区三| 亚洲国产欧美网| 国产日本99.免费观看| 两人在一起打扑克的视频| 欧美3d第一页| 成人亚洲精品av一区二区| 日本一本二区三区精品| 99精品在免费线老司机午夜| 欧美一级a爱片免费观看看| 国产亚洲精品综合一区在线观看| 亚洲欧美日韩东京热| 国产午夜精品论理片| 精品不卡国产一区二区三区| 少妇人妻一区二区三区视频| 天天躁狠狠躁夜夜躁狠狠躁| 18禁国产床啪视频网站| 国产精品98久久久久久宅男小说| 国产毛片a区久久久久| 麻豆成人av在线观看| 窝窝影院91人妻| 亚洲无线观看免费| 99久久综合精品五月天人人| svipshipincom国产片| 精品国内亚洲2022精品成人| 国产精品久久电影中文字幕| 最新美女视频免费是黄的| 国产高清三级在线| 人人妻人人看人人澡| 黄频高清免费视频| 嫁个100分男人电影在线观看| 亚洲人成网站高清观看| 久久性视频一级片| 老司机午夜十八禁免费视频| 国产人伦9x9x在线观看| 亚洲国产欧美网| 精品不卡国产一区二区三区| 最新中文字幕久久久久 | 亚洲国产日韩欧美精品在线观看 | 日韩av在线大香蕉| 99在线视频只有这里精品首页| 人妻久久中文字幕网| 一级作爱视频免费观看| 成年人黄色毛片网站| 97超级碰碰碰精品色视频在线观看| 国产激情欧美一区二区| 国产精品野战在线观看| 国产成+人综合+亚洲专区| 日本与韩国留学比较| 99精品欧美一区二区三区四区| 亚洲av熟女| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品综合一区在线观看| 午夜福利在线观看免费完整高清在 | 成人av一区二区三区在线看| 久久婷婷人人爽人人干人人爱| 免费看美女性在线毛片视频| 国产精品久久久久久人妻精品电影| 亚洲人成电影免费在线| 精品午夜福利视频在线观看一区| 日本黄色片子视频| 美女扒开内裤让男人捅视频| 一级毛片高清免费大全| 搡老熟女国产l中国老女人| 亚洲中文字幕一区二区三区有码在线看 | 人人妻人人看人人澡| 不卡一级毛片| 搡老岳熟女国产| 淫秽高清视频在线观看| 草草在线视频免费看| 1024香蕉在线观看| 亚洲av熟女| 一本一本综合久久| 国产人伦9x9x在线观看| 日韩欧美在线乱码| 久久香蕉精品热| 老熟妇乱子伦视频在线观看| 国产99白浆流出| 久久人人精品亚洲av| 国产真人三级小视频在线观看| 精品乱码久久久久久99久播| 欧美性猛交╳xxx乱大交人| 国产成人精品无人区| 国产精华一区二区三区| 一进一出抽搐动态| 成年免费大片在线观看| 欧美成人免费av一区二区三区| 美女被艹到高潮喷水动态| 色老头精品视频在线观看| 国产精品亚洲一级av第二区| 一进一出抽搐gif免费好疼| 999久久久国产精品视频| av欧美777| 男女床上黄色一级片免费看| www日本黄色视频网| 欧美又色又爽又黄视频| 国产成人系列免费观看| 国产乱人视频| 在线视频色国产色| 国产一区二区在线av高清观看| 欧美日韩中文字幕国产精品一区二区三区| 听说在线观看完整版免费高清| 极品教师在线免费播放| 嫁个100分男人电影在线观看| 免费看日本二区| 一级毛片女人18水好多| 天堂网av新在线| 国产午夜精品论理片| 中文资源天堂在线| 悠悠久久av| 少妇人妻一区二区三区视频| 国产一区在线观看成人免费| 国产三级黄色录像| 久久久久久久精品吃奶| av片东京热男人的天堂| 国产综合懂色| 又黄又爽又免费观看的视频| 男人舔女人的私密视频| 免费av毛片视频| 18禁美女被吸乳视频| 欧美日韩一级在线毛片| 亚洲人成网站高清观看| 成人国产综合亚洲| 日韩欧美精品v在线| 中文字幕人妻丝袜一区二区| 91av网站免费观看| 欧美日韩福利视频一区二区| 一本精品99久久精品77| 欧美不卡视频在线免费观看| 一本精品99久久精品77| 中文字幕最新亚洲高清| 国产精品电影一区二区三区| 变态另类成人亚洲欧美熟女| 精华霜和精华液先用哪个| 免费大片18禁| 久久精品91无色码中文字幕| 精品久久久久久成人av| 亚洲无线观看免费| 国产三级中文精品| 国产成年人精品一区二区| 丝袜人妻中文字幕| 中出人妻视频一区二区| 精品久久久久久久人妻蜜臀av| 欧美激情久久久久久爽电影| 美女大奶头视频| 黄片小视频在线播放| 男女视频在线观看网站免费| 欧美色视频一区免费| 亚洲天堂国产精品一区在线| 黑人巨大精品欧美一区二区mp4| 一a级毛片在线观看| 免费观看人在逋| 97超级碰碰碰精品色视频在线观看| 可以在线观看毛片的网站| 国产精品一区二区精品视频观看| 在线视频色国产色| av在线天堂中文字幕| 一级作爱视频免费观看| 日本一本二区三区精品| 亚洲精品一卡2卡三卡4卡5卡| 啦啦啦免费观看视频1| 黄色视频,在线免费观看| 精品熟女少妇八av免费久了| 欧美一级a爱片免费观看看| 一区二区三区激情视频| 日本 av在线| 日韩av在线大香蕉| 午夜久久久久精精品| 国产毛片a区久久久久| 亚洲精品一区av在线观看| 色综合站精品国产| 亚洲欧美一区二区三区黑人| 中亚洲国语对白在线视频| 精品久久久久久成人av| 日日摸夜夜添夜夜添小说| 免费在线观看日本一区| 搡老熟女国产l中国老女人| 久久国产精品影院| 亚洲 国产 在线| a级毛片a级免费在线| 淫秽高清视频在线观看| 日本一二三区视频观看| 校园春色视频在线观看| 亚洲人成网站高清观看| 色综合婷婷激情| 久久久久久久久久黄片| 亚洲精品久久国产高清桃花| 日本一本二区三区精品| 91在线观看av| 999久久久精品免费观看国产| 国产精品一区二区三区四区久久| 国产精品日韩av在线免费观看| 久久九九热精品免费| 精品久久久久久久人妻蜜臀av| 国产精品香港三级国产av潘金莲| 少妇裸体淫交视频免费看高清| 亚洲人成电影免费在线| 99热精品在线国产| 日韩欧美国产一区二区入口| 极品教师在线免费播放| 夜夜夜夜夜久久久久| 9191精品国产免费久久| 老司机在亚洲福利影院| 九九热线精品视视频播放| 午夜免费观看网址| 午夜福利在线在线| 亚洲狠狠婷婷综合久久图片| 又大又爽又粗| 法律面前人人平等表现在哪些方面| 亚洲专区国产一区二区| 国产高清激情床上av| 亚洲国产欧美一区二区综合| 欧美高清成人免费视频www| 村上凉子中文字幕在线| 99国产综合亚洲精品| 亚洲人成网站高清观看| 国产高清视频在线播放一区| 噜噜噜噜噜久久久久久91| 精品日产1卡2卡| 老汉色∧v一级毛片| 天堂√8在线中文| 国产欧美日韩一区二区精品| 午夜久久久久精精品| 国产极品精品免费视频能看的| 黄片小视频在线播放| 成人亚洲精品av一区二区| 两个人看的免费小视频| 成人永久免费在线观看视频| 国产欧美日韩一区二区三| 亚洲av熟女| 日本在线视频免费播放| 特级一级黄色大片| 久久久国产欧美日韩av| 亚洲九九香蕉| 淫妇啪啪啪对白视频| 999精品在线视频| 国产成人福利小说| 亚洲精品色激情综合| 亚洲狠狠婷婷综合久久图片| 日本撒尿小便嘘嘘汇集6| 99国产精品99久久久久| 亚洲国产精品999在线| 男女视频在线观看网站免费| 99久久精品一区二区三区| 12—13女人毛片做爰片一| 男插女下体视频免费在线播放| 国产精品电影一区二区三区| 成人av一区二区三区在线看| 叶爱在线成人免费视频播放| 国产精品永久免费网站| 99久久久亚洲精品蜜臀av| 亚洲九九香蕉| 成人18禁在线播放| 欧美三级亚洲精品| 88av欧美| 亚洲人成电影免费在线| 亚洲av中文字字幕乱码综合| 久久久久亚洲av毛片大全| 男女视频在线观看网站免费| 99精品在免费线老司机午夜| 精品人妻1区二区| 色播亚洲综合网| 婷婷丁香在线五月| 国产精品99久久久久久久久| 国产一区二区在线观看日韩 | 老司机午夜福利在线观看视频| 亚洲色图av天堂| 亚洲精品色激情综合| aaaaa片日本免费| 免费看光身美女| 国产精品一及| 在线免费观看不下载黄p国产 | 一二三四社区在线视频社区8| 亚洲一区高清亚洲精品| 99视频精品全部免费 在线 | 亚洲色图 男人天堂 中文字幕| 成年女人毛片免费观看观看9| 小蜜桃在线观看免费完整版高清| 男女午夜视频在线观看| a级毛片a级免费在线| 人妻久久中文字幕网| 欧美日韩瑟瑟在线播放| 免费电影在线观看免费观看| 欧美成狂野欧美在线观看| 九色国产91popny在线| 一个人看视频在线观看www免费 | 国产高清三级在线| 国产精品影院久久| 日本免费a在线| 精品一区二区三区四区五区乱码| 熟妇人妻久久中文字幕3abv| 午夜福利免费观看在线| 成人欧美大片| 久久久水蜜桃国产精品网| 国产伦一二天堂av在线观看| 精品免费久久久久久久清纯| 国产精品99久久久久久久久| 欧美成人一区二区免费高清观看 | 久久中文字幕人妻熟女| 波多野结衣巨乳人妻| 老鸭窝网址在线观看| 久久久久国产精品人妻aⅴ院| 最近视频中文字幕2019在线8| 国产精品99久久99久久久不卡| 国产精品久久视频播放| 亚洲熟女毛片儿| 听说在线观看完整版免费高清| 小说图片视频综合网站| 日韩高清综合在线| 国产午夜福利久久久久久| 一夜夜www| 老司机在亚洲福利影院| 两个人视频免费观看高清| 99久久99久久久精品蜜桃| 亚洲精品久久国产高清桃花| 国产视频一区二区在线看| 国模一区二区三区四区视频 | 免费观看的影片在线观看| 成人特级黄色片久久久久久久| 俺也久久电影网| 色在线成人网| 午夜福利成人在线免费观看| 99久久精品热视频| 亚洲av成人一区二区三| 久久久精品欧美日韩精品| 亚洲午夜理论影院| 高清毛片免费观看视频网站| 精品国产亚洲在线| 午夜福利视频1000在线观看| 亚洲av成人不卡在线观看播放网| 午夜激情欧美在线| www.自偷自拍.com| 国内久久婷婷六月综合欲色啪| 丰满的人妻完整版| 国产综合懂色| 日本a在线网址| 一级黄色大片毛片| 色播亚洲综合网| 亚洲av成人不卡在线观看播放网| 日韩有码中文字幕| 丁香欧美五月| 淫妇啪啪啪对白视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看成人毛片| 一区福利在线观看| 美女扒开内裤让男人捅视频| 精品电影一区二区在线| xxxwww97欧美| 黄色女人牲交| 中文字幕久久专区| 91久久精品国产一区二区成人 | 岛国在线观看网站| 最近在线观看免费完整版| 亚洲一区二区三区色噜噜| 校园春色视频在线观看| 18禁黄网站禁片免费观看直播| 国产精品野战在线观看| 最近最新免费中文字幕在线| av天堂中文字幕网| e午夜精品久久久久久久| 午夜福利欧美成人| 亚洲午夜理论影院| 亚洲九九香蕉| av福利片在线观看| 嫁个100分男人电影在线观看| 精华霜和精华液先用哪个| 99热这里只有精品一区 | 成年免费大片在线观看| 成年女人看的毛片在线观看| 天天躁日日操中文字幕| 三级男女做爰猛烈吃奶摸视频| 亚洲成a人片在线一区二区| 2021天堂中文幕一二区在线观| 亚洲av成人精品一区久久| 亚洲熟妇中文字幕五十中出| 国产精品一区二区三区四区免费观看 | 亚洲国产精品合色在线| 亚洲人与动物交配视频| 欧美日本亚洲视频在线播放| 脱女人内裤的视频| 久久久国产精品麻豆| 亚洲熟女毛片儿| bbb黄色大片| 深夜精品福利|