999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

北洛河流域植物多樣性地理格局與環(huán)境關(guān)系

2021-09-23 12:40:12田起隆許小明王浩嘉雷斯越易海杰鄒亞東王妙倩張曉萍1
中國環(huán)境科學(xué) 2021年9期
關(guān)鍵詞:物種植物

田起隆,許小明,呂 渡,王浩嘉,雷斯越,易海杰,賀 潔,何 亮,薛 帆,鄒亞東,王妙倩,張曉萍1,*

北洛河流域植物多樣性地理格局與環(huán)境關(guān)系

田起隆1,2,許小明3,呂 渡1,2,王浩嘉3,雷斯越3,易海杰1,2,賀 潔3,何 亮3,薛 帆3,鄒亞東3,王妙倩3,張曉萍1,3*

(1.中國科學(xué)院水利部水土保持研究所,黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點(diǎn)實(shí)驗(yàn)室,陜西 楊凌 712100;2.中國科學(xué)院大學(xué),北京 100049;3.西北農(nóng)林科技大學(xué)水土保持研究所,陜西 楊凌 712100)

通過野外調(diào)查和歷史文獻(xiàn)及標(biāo)本統(tǒng)計(jì)分析,確定并分析北洛河流域種子植物信息數(shù)據(jù).結(jié)果表明:該區(qū)共有種子植物123科581屬1671種,占黃土高原種子植物的83.67%、67.24%、51.83%,集中分布于大型科(>50種)、較大科(21~50種)和單種屬(1種)、小型屬(2~5種)內(nèi).11個(gè)干流流經(jīng)的地理單元植物相似性聚類分為5組,較早分離的是黃龍縣,最大類因子午嶺山脈連接而包含4個(gè)地理單元(富縣、合水縣、宜君縣、黃陵縣).富縣、合水縣和黃龍縣植物豐富度SD值最高,洛川縣、華池縣和甘泉縣SD值最低.該區(qū)地理成分在屬級水平上有15個(gè)分布區(qū)類型16個(gè)變型,以溫帶分布及其變型為主(287屬,占49.40%),并與熱帶分布及其變型聯(lián)系緊密(115屬,占19.79%),植物豐富度SD值與溫帶成分顯著正相關(guān);黃陵縣、洛川縣、黃龍縣/值最高,熱帶性質(zhì)最強(qiáng),華池縣、定邊縣/值最低,熱帶性質(zhì)最弱./值受海拔影響,隨溫度和降水從東南向西北遞減.該區(qū)植物起源古老,分化程度高,物種多樣性豐富,地理成分復(fù)雜,是重要的植物種質(zhì)資源庫.植物物種多樣性地帶性格局過渡性明顯,體現(xiàn)出這些類群的生態(tài)位需求.需加強(qiáng)保護(hù)該區(qū)原生環(huán)境和植物,合理開發(fā)利用,可選擇6大優(yōu)勢科(菊科、禾本科、薔薇科、豆科、唇形科和毛茛科)中適地適樹適草的物種(溫帶性質(zhì)),進(jìn)行植被恢復(fù)和水土流失治理.

物種多樣性;科屬組成;地理成分;豐富度;溫帶性質(zhì)

植物物種多樣性是生物多樣性的重要組成部分,是植物在一定自然地理環(huán)境,特別是自然歷史等綜合條件作用下長期適應(yīng)進(jìn)化的結(jié)果,蘊(yùn)涵大量生物與非生物環(huán)境相互作用的生態(tài)信息,其地理格局更是自然地理環(huán)境的反應(yīng)和環(huán)境變遷的鑒證或依據(jù)[1-3].北洛河流域地處黃河中游,地形地貌復(fù)雜,溝壑縱橫,具有黃土高原典型的丘陵溝壑、高塬溝壑和土石山等地貌特征;跨越溫帶半干旱和大陸性季風(fēng)氣候帶;生境破碎化嚴(yán)重、植被類型典型多樣、物種豐富[4-5].但該區(qū)年降水時(shí)空異質(zhì)性強(qiáng)且多暴雨,破壞性大,加之生態(tài)環(huán)境脆弱、水土流失嚴(yán)重,導(dǎo)致該地區(qū)災(zāi)害頻繁、生產(chǎn)落后,人民貧困[6].

近年來學(xué)者們對水土流失現(xiàn)狀及防治做了大量的研究,在此過程中深刻地認(rèn)識到植被對水土流失和氣候的調(diào)控作用.它能削減降雨能量、增加入滲,減緩徑流,減輕和控制侵蝕敏感因子,提高土壤抗沖性[7-13];也能增加降水,提高濕度,減小溫差,影響水文循環(huán),進(jìn)而改善區(qū)域氣候和生態(tài)環(huán)境,因此對植被的研究很有必要[14-19].植被治理措施中適地適樹適草的概念成為熱點(diǎn),植物適合種在什么區(qū)域抑或是什么區(qū)域適合種植什么植物來防治水土流失、與環(huán)境協(xié)調(diào)發(fā)展、發(fā)揮功能特性這都體現(xiàn)著該區(qū)植物多樣性研究的重要性.目前,大多研究主要集中于植被演替下群落多樣性和子午嶺、黃龍山植物區(qū)系等方面[20-22],關(guān)于環(huán)境因子驅(qū)動(dòng)下區(qū)域尺度植物多樣性之間的關(guān)系,需進(jìn)一步深入研究.這對理解植物的起源、遷移和分布等方面具有重要意義.

為此,對流域內(nèi)植物多樣性進(jìn)行分析,揭示不同地區(qū)植物的科屬組成和地理成分,對比分析各地理單元植物多樣性分布格局與環(huán)境因子關(guān)系,為植物資源的合理利用、生態(tài)環(huán)境建設(shè)、生物多樣性保護(hù)及黃土高原水土流失治理提供科學(xué)依據(jù).

1 材料與方法

1.1 研究區(qū)概況

北洛河流域地處黃河中游,是黃河的二級支流、渭河的一級支流,干流發(fā)源于陜西省定邊縣白于山郝莊梁,流經(jīng)吳起、志丹、甘泉、富縣、洛川等縣,在大荔縣韋林鎮(zhèn)的倉西村匯入渭河.地理坐標(biāo)為東經(jīng)107°34′45′′~110°01′37′′,北緯34°57′13′′~ 37°17′03′′,是陜西省縱跨緯度最大的河流,河道全長680km,流域面積為2.69萬 km2.屬溫帶半干旱和大陸性季風(fēng)氣候,四季分明,春季干燥、夏季旱澇相間、秋季濕潤、冬季寒冷干燥.多年平均降水量514mm,由南向北、自東向西遞減,降雨年內(nèi)分布不均,集中在5~9月份,且多暴雨,破壞性大,水土流失極為嚴(yán)重.該區(qū)地形地貌復(fù)雜多變,地形破碎、溝壑縱橫、山高坡陡,具有典型的丘陵溝壑地貌,可分為黃土丘陵溝壑區(qū)、黃土高塬溝壑區(qū)、土石山區(qū)、渭河、洛河沖積平原區(qū)及風(fēng)沙區(qū).主要土壤類型為黃綿土、黑壚土和灰褐土.

圖1 北洛河流域主要地貌特征及涉及縣域

1.風(fēng)沙丘陵區(qū);2.丘陵溝壑區(qū);3.土石山區(qū);4.高塬溝壑區(qū);5.階地平原區(qū)

1.2 研究方法

1.2.1 數(shù)據(jù)獲取 本研究選取流域中上游典型地貌植被區(qū)干流所流經(jīng)的定邊縣、吳起縣、合水縣、富縣、黃陵縣等中國11個(gè)縣級地理樣本單元(圖1),于2016~2020年,根據(jù)不同地貌植被類型劃分,選擇代表性林分類型,設(shè)置面積為20m×20m(喬木)、10m×10m(灌木)和1m×1m(草叢)樣方的標(biāo)準(zhǔn)樣地,進(jìn)行野外植物群落調(diào)查,記錄其中出現(xiàn)的物種名稱[23].并根據(jù)國家標(biāo)本信息基礎(chǔ)設(shè)施(http://www. nsii.org.cn/)標(biāo)本館標(biāo)本和大量的分布數(shù)據(jù),結(jié)合西北農(nóng)林科技大學(xué)生命科學(xué)學(xué)院植物標(biāo)本館(WUK)、中國科學(xué)院植物研究所標(biāo)本館(PE)等19個(gè)標(biāo)本館記錄,參照《中國植物志》、《Flora of China》、《黃土高原植物志》、《甘肅植物志》、《陜西維管植物名錄》、《子午嶺木本植物志》、《陜西子午嶺植物圖鑒》、《西北農(nóng)牧交錯(cuò)帶常見植物圖譜》和各地《縣志》等書籍及相關(guān)文獻(xiàn),運(yùn)用數(shù)理統(tǒng)計(jì)方法、構(gòu)建研究區(qū)種子植物信息數(shù)據(jù)庫[24].經(jīng)過中國數(shù)字植物標(biāo)本館(http://www.cvhac.cn/)和The Plant List (http: //www.theplantlist.org/)對物種拉丁名反復(fù)核實(shí)(避免因科屬調(diào)整和一種多名等因素造成的誤差),去除栽培引種,確定該區(qū)植物名錄,編制每個(gè)縣的種子植物名錄,確定各縣種子植物組成及地理成分.

年均溫、年均降水量和各地理單元的海拔平均值數(shù)據(jù)來源于WorldClim數(shù)據(jù)庫(2.1版, http://www. worldclim.org)

1.2.2 統(tǒng)計(jì)分析

(1)采用S?renson系數(shù),進(jìn)行等級聚類分析該區(qū)各縣植物組成的相似性,計(jì)算公式為:

=[2/(+)]×100% (1)

式中:為各縣物種相似性系數(shù);分別為兩地全部物種數(shù);為兩地共有物種數(shù)[25].

(2)植物種系分化度分析,用來評判該區(qū)域植物豐富度的數(shù)值[26].某一區(qū)域植物中,每個(gè)科所含屬越多,每個(gè)屬所含種越多,被視為植物種系分化度越大,物種組成越豐富.用以下公式計(jì)算得出[27].

SD=2/1+3/1 (2)

式中:123,分別表示植物物種的科、屬、種的數(shù)量;SD值表示該植物的種系分化度.

(3)/值為植物區(qū)系中熱帶成分()屬與溫帶成分()屬的比值,是衡量植物區(qū)系性質(zhì)的一個(gè)指標(biāo),/值越高植物區(qū)系的熱帶性質(zhì)越強(qiáng).是各種熱帶分布類型屬(2~7)之和,是各種溫帶分布類型屬(8~14)之和[28].

數(shù)據(jù)處理、分析采用Excel 2010、SPSS 26.0,地圖繪制由ArcGIS 10.5完成.

2 結(jié)果與分析

2.1 植物種類組成及相似性

據(jù)統(tǒng)計(jì),北洛河流域種子植物共有123科581屬1671種,其中裸子植物4科8屬13種,分別占該區(qū)種子植物總科、屬、種數(shù)的3.25%、1.38%、0.78%;被子植物119科573屬1658種,分別占該區(qū)種子植物總科、屬、種數(shù)的96.75%、98.62%、99.22%,在數(shù)量上占絕對優(yōu)勢.被子植物中雙子葉植物有97科460屬1368種,占被子植物總數(shù)的82.51%;單子葉植物有22科113屬290種,占17.49%,雙子葉植物優(yōu)勢明顯.

圖2 北洛河流域種子植物相似性聚類樹狀圖

G表示分組類型;Ⅰ.定邊縣;Ⅱ.吳起縣;Ⅲ.華池縣;Ⅳ.志丹縣;Ⅴ.合水縣;Ⅵ.富縣;Ⅶ.甘泉縣;Ⅷ.黃陵縣;Ⅸ.宜君縣;Ⅹ.洛川縣;Ⅺ.黃龍縣

基于不同縣級地理樣本單元種子植物相似性S?renson系數(shù)構(gòu)建物種相似性歐式距離矩陣,采用Ward法進(jìn)行聚類分析(圖2),可以反應(yīng)出各地多樣性特征間的聯(lián)系,理解其發(fā)展演化過程,為判斷區(qū)域植物過渡屬性提供科學(xué)依據(jù)[29].11個(gè)干流流經(jīng)的地理單元植物聚類分為5組,即G1組:Ⅺ黃龍縣;G2組:Ⅲ華池縣、Ⅱ吳起縣;G3組:Ⅶ 甘泉縣、Ⅳ志丹縣;G4組:Ⅹ洛川縣、Ⅰ定邊縣;G5組:Ⅵ富縣、Ⅴ合水縣、Ⅸ宜君縣、Ⅷ黃陵縣.聚類分析樹狀圖整體反映該區(qū)植物物種相似性變化的地理分布格局.

2.2 植物科屬構(gòu)成多樣性

科作為植物分類學(xué)中較自然的演化單位,是揭示該區(qū)植物多樣性特征的重要途徑[30].研究區(qū)內(nèi)野生種子植物共123科,對所含的屬、種數(shù)進(jìn)行統(tǒng)計(jì)(表1),結(jié)果顯示,其中>50種的大型科有6科,為菊科Asteraceae(59屬204種)、禾本科Poaceae(52屬127種)、薔薇科Rosaceae(24屬120種)、豆科Fabaceae(34屬115種)、唇形科Lamiaceae(26屬65種)和毛茛科Ranunculaceae(12屬61種),占該區(qū)種子植物總科數(shù)的4.88%(比例最低),但占總屬數(shù)、總種數(shù)的35.63%、41.41%,是該區(qū)的基礎(chǔ)和主體,表明植物種類集中在>50種的科中,優(yōu)勢現(xiàn)象十分顯著,無論是在屬級還是種級水平上均是該區(qū)的優(yōu)勢科.含21~50種的較大科有14科130屬403種,分別占該區(qū)總科、屬、種數(shù)的11.38%、22.37%、24.12%;含11~20種的中型科有19科,含83屬270種,分別占該區(qū)總科、屬、種數(shù)的15.45%、14.29%、16.16%;含2~10種的寡種科有57科134屬279種,分別占該區(qū)總科、屬、種數(shù)的46.34%、23.06%、16.70%;單種科有27科,分別占該區(qū)總科、屬、種數(shù)的21.95%、4.65%、1.61%.

表1 北洛河流域種子植物科的組成

屬是反映植物系統(tǒng)發(fā)育過程和多樣性特征的重要研究對象[31].對本區(qū)野生種子植物581屬所含的種數(shù)進(jìn)行統(tǒng)計(jì)(表2),結(jié)果顯示,主要有單種屬(含1種)、小型屬(2~5種)、中型屬(6~9種)、較大屬(10~14種)和大型屬(315種)5個(gè)類型.該區(qū)含315種的大型屬有10屬229種,為蒿屬(49種)、薹草屬(24種)、委陵菜屬(23種)、柳屬(22種)、堇菜屬(22種)、鐵線蓮屬(21種)、錦雞兒屬(18種)、黃芪屬(18種)、萹蓄屬(17種)和披堿草屬(15種),是該區(qū)的優(yōu)勢屬,分別占該區(qū)總屬數(shù)、總種數(shù)的1.72%、13.71%;含10~14種的較大屬有20屬234種,分別占總屬數(shù)總種數(shù)的3.44%、14.00%;含6~9種的中型屬有37屬270種,分別占總屬數(shù)總種數(shù)的6.37%、16.16%;含2~5種的小型屬有225屬649種,分別占總屬數(shù)總種數(shù)的38.73%、38.84%;單種屬有289屬,占總屬數(shù)的49.74%,在各類型屬中所占比例最高,其所含種數(shù)占總種數(shù)的17.29%,僅次于小型屬類型所含的種數(shù).因此,屬型的統(tǒng)計(jì)分析表明,該區(qū)88.47%的植物種類集中于小型屬和單種屬,它們在該區(qū)植物組成中占絕對優(yōu)勢.也可以說單種屬構(gòu)成了該區(qū)植物屬的主體,而小型屬所包含的植物種類構(gòu)成了該區(qū)植物種的主體.

表2 北洛河流域種子植物屬的組成

2.3 植物屬的分布區(qū)類型

根據(jù)吳征鎰[32]對中國種子植物屬的分布區(qū)類型劃分,本研究區(qū)種子植物共581屬,可以將其劃分為15個(gè)分布區(qū)類型和16個(gè)變型.其中世界分布型有77屬,占區(qū)域內(nèi)分布屬的13.25%;北溫帶分布有140屬,占總屬數(shù)的24.1%,其變型北溫帶和南溫帶(全溫帶)間斷分布有34屬、變型歐亞和南美溫帶間斷分布有4屬,分別占總屬數(shù)的5.85%、0.69%;泛熱帶分布有58屬,占總屬數(shù)的9.98%,其變型熱帶亞洲、非洲和中、南美洲間斷分布僅1種,占總屬數(shù)的0.17%;舊世界溫帶分布有48屬,占總屬數(shù)的8.26%.其有2個(gè)變型,地中海、西亞(或中亞)和東亞間斷分布有8屬、歐亞和南部非洲(有時(shí)也在大洋洲)間斷分布有5屬,分別占總屬數(shù)的1.38%、0.86%;東亞和北美洲間斷分布有30屬,占總屬數(shù)的5.17%;東亞分布有27屬,占總屬數(shù)的4.65%.其有2個(gè)變型,中國-喜馬拉雅有7屬、中國-日本分布有15屬,分別占總屬數(shù)的1.21%、2.58%;地中海、西亞至中亞分布有20屬,占總屬數(shù)的3.44%.其有3個(gè)變型,地中海區(qū)至西亞或中亞和墨西哥或古巴間斷分布、地中海至溫帶-熱帶亞洲均為2屬,占總屬數(shù)的0.34%,地中海至中亞和南非洲、大洋洲間斷分布僅有1屬;溫帶亞洲分布有18屬,占總屬數(shù)的3.1%;熱帶亞洲至熱帶大洋洲分布、中國特有分布有13屬,占總屬數(shù)的2.24%;舊世界熱帶分布、熱帶亞洲至熱帶非洲分布有12屬,占總屬數(shù)的2.07%;熱帶亞洲(印度一馬來西亞)分布有11屬,占總屬數(shù)的1.9%.有2個(gè)變型,越南(或中南半島)至華南(或西南)分布、爪哇(或蘇門達(dá)臘)、喜馬拉雅間斷或星散分布到華南、西南均只有1屬;中亞分布有6屬,占總屬數(shù)的1.03%.有3變型,中亞東部(亞洲中部)分布有6屬,西亞至西喜馬拉雅和西藏分布有2屬,中亞至喜馬拉雅和我國西南分布有1屬;熱帶亞洲和熱帶美洲間斷分布有5屬,占總屬數(shù)的0.86%;大洋洲和南美洲間斷分布有2屬,占總屬數(shù)的0.34%;熱帶亞洲、非洲(或東非、馬達(dá)加斯加)和大洋洲間斷分布只有1屬,僅占屬數(shù)的0.17%(表3).

表3 北洛河流域種子植物屬的分布區(qū)類型統(tǒng)計(jì)

除世界分布屬外,研究區(qū)內(nèi)植物屬的分布型包括熱帶成分、溫帶成分、古地中海成分、東亞成分及中國特有成分.該區(qū)溫帶成分4分布型4變型,共計(jì)287屬,占總屬的49.40%,溫帶成分優(yōu)勢較明顯;熱帶成分6分布型4變型,有115屬,占總屬的19.79%,在區(qū)域內(nèi)占有一定比例;東亞成分有1個(gè)分布型2個(gè)分布變型,49屬,占總屬的8.43%;古地中海成分2分布型6變型,共40屬,占總屬的6.88%;中國特有成分只有13屬,僅占2.24%.

2.4 植物多樣性與地理成分相關(guān)性

該區(qū)涉及11個(gè)縣,各縣科屬組成為(圖3):洛川縣67科177屬258種,植物種系分化度SD值為6.49,/值為0.29;華池縣82科237屬363種,SD值為7.32,/值為0.18;甘泉縣86科256屬406種,SD值為7.70,/值為0.25;宜君縣97科304屬559種,SD值為8.90,/值為0.27;黃陵縣97科323屬568種,SD值為9.19,/值為0.31;志丹縣81科288屬512種,SD值為9.88,/值為0.24;吳起縣79科270屬523種,SD值為10.04,/值為0.22;定邊縣68科229屬459種,SD值為10.12,/值為0.19;黃龍縣94科341屬693種,SD值為11.00,/值為0.28;合水縣93科314屬715種,SD值為11.06,/值為0.21;富縣108科398屬801種,SD值為11.10,/值為0.26.

該區(qū)富縣、合水縣和黃龍縣植物豐富度最高, SD值均在11以上;洛川縣、華池縣和甘泉縣的豐富度最低,SD值均在8以下;宜君縣、黃陵縣、志丹縣、吳起縣和定邊縣植物豐富度居中,但定邊縣僅459種,SD卻達(dá)到10.12,這說明定邊縣與含有500種以上的宜君縣、黃陵縣、志丹縣、吳起縣相比物種豐富度高.黃陵縣、洛川縣、黃龍縣熱帶性質(zhì)最強(qiáng),華池縣、定邊縣熱帶性質(zhì)最低,這種分布模式是氣候與環(huán)境共同作用的結(jié)果.

圖3 北洛河流域種子植物多樣性特征

SD.植物區(qū)系種系分化度;R/T.熱帶成分/溫帶成分;Ⅰ.定邊縣;Ⅱ.吳起縣;Ⅲ.華池縣;Ⅳ.志丹縣;Ⅴ.合水縣;Ⅵ.富縣;Ⅶ.甘泉縣;Ⅷ.黃陵縣;Ⅸ.宜君縣;Ⅹ.洛川縣;Ⅺ.黃龍縣

表4 北洛河流域種子植物SD值與分布區(qū)類型關(guān)系

注:*表示在0.05水平上顯著相關(guān).

將11個(gè)縣級地理樣本單元SD值與/值,值(熱帶成分)與值(溫帶成分)做相關(guān)分析發(fā)現(xiàn):本研究區(qū)植物豐富度與溫帶成分顯著正相關(guān),熱帶成分與/值顯著正相關(guān)(表4).這說明溫帶成分的多樣性分布格局是導(dǎo)致流域內(nèi)不同區(qū)域種系分化的基礎(chǔ).熱帶成分的滲入是導(dǎo)致各區(qū)域植物多樣性差異的原因.溫帶成分更適合在此區(qū)域內(nèi)生長發(fā)育,保持物種延續(xù).所以,研究/值與環(huán)境因子關(guān)系有利于揭示該區(qū)植物多樣性分布格局,該區(qū)域內(nèi)植被恢復(fù)應(yīng)優(yōu)先選用溫帶性質(zhì)植物.

2.5 R/T值與環(huán)境因子關(guān)系

通過將各地理單元/值與環(huán)境因子回歸分析(圖4)發(fā)現(xiàn),/值隨海拔升高呈現(xiàn)降低趨勢圖4(a),這表明海拔是限制熱帶成分?jǐn)U散分布的重要因素./值隨年均溫升高呈增大趨勢,這與/和年均降水量變化趨勢相同.這說明溫度和降水是影響熱帶成分植物遷移的直接生態(tài)因子[圖4(b)、圖4(c)].這些因子的共同作用解釋了,本流域植物區(qū)系呈現(xiàn)過渡性特征明顯的原因,對于區(qū)域植被恢復(fù)和生態(tài)建設(shè)具有重要意義.

圖4 北洛河流域R/T值與環(huán)境因子關(guān)系

3 討論

本研究發(fā)現(xiàn),北洛河流域種子植物共有123科581屬1671種,占黃土高原種子植物的83.67%、67.24%、51.83%[33],占中國種子植物的36.50%、18.16%和6.13%[34].該區(qū)具有很高的物種多樣性保護(hù)價(jià)值,是黃土高原代表性的植物種質(zhì)資源庫,對地區(qū)生態(tài)建設(shè)和經(jīng)濟(jì)發(fā)展具有重要作用.其中單種科和寡種科共有84科,占總科數(shù)的68.29%(比例最高),是科組成的主體,但種數(shù)僅占到總種數(shù)的18.31%;大型科和較大科共20科337屬1095種,只占總科數(shù)的16.26%,屬、種卻占到58%、65.53%,這說明屬和種趨向于集中在有限的少數(shù)科內(nèi);其中優(yōu)勢科為菊科、禾本科、薔薇科、豆科、唇形科和毛茛科.該區(qū)含315種的大型屬有10屬229種,優(yōu)勢屬,分別為:蒿屬、薹草屬、委陵菜屬、柳屬、堇菜屬、鐵線蓮屬、錦雞兒屬、黃芪屬、萹蓄屬和披堿草屬.單種屬與小型屬合計(jì)514屬(含938種),分別占總屬數(shù)、總種數(shù)的88.47%、56.13%,該區(qū)主要的優(yōu)勢種基本上均來自該類型,是在該區(qū)植物中占絕對優(yōu)勢,表現(xiàn)出該區(qū)保存有較多的殘遺屬和新建屬,有利于古老殘遺植物類群的保存與新建類群的分化遷移定居[35].這說明本區(qū)植物起源古老,分化程度高,組成復(fù)雜[36-37].綜上,可從6大優(yōu)勢科(菊科、禾本科、薔薇科、豆科、唇形科和毛茛科)中選擇適地適樹適草的物種,進(jìn)行當(dāng)?shù)刂脖换謴?fù)及建設(shè).

根據(jù)不同縣級地理樣本單元種子植物相似性聚類分析可知,11個(gè)干流流經(jīng)的地理單元種子植物多樣性可分為5大類:第一大類黃龍縣,位于流域的東南端,在相似性聚類分析中與其他縣較早分離,因山脈、河流或峽谷阻隔而產(chǎn)生地理隔離,使個(gè)體之間遺傳信息不能交換,經(jīng)長期演化和分化易形成新物種[38];第二大類將華池縣及吳起縣聚類,表明這兩個(gè)區(qū)域大部分都處于溫帶森林草原亞地帶,地形地貌及氣候環(huán)境特征相似,因此兩個(gè)區(qū)域的植物組成相似性較高;第三大類將甘泉縣及志丹縣聚類,均處于溫帶森林草原向暖溫帶北部落葉櫟林亞地帶的過渡地帶[39];第四大類將洛川縣及定邊縣聚類,但相似性系數(shù)僅32.9%(<50%),表明二者間的聯(lián)系較弱[25].這由于兩地農(nóng)耕歷史悠久,受人類活動(dòng)影響,使該地區(qū)物種數(shù)較低易于聚類.第五大類將富縣、合水縣、宜君縣和黃陵縣聚類,都屬于土石山區(qū),有子午嶺連接且地理位置更為接近,共有種的比例很高,因此植物類型相似度很高[37].聚類結(jié)果符合張振萬等[40]關(guān)于陜北黃土高原植物分區(qū)界線的初步劃分.

該區(qū)種子植物地理成分分析表明,從屬的分布區(qū)類型來看,該區(qū)含有15個(gè)分布區(qū)類型和16個(gè)變型,溫帶性質(zhì)明顯(特別是“北溫帶分布”成分),熱帶成分占有一定比例,有部分東亞成分和古地中海成分,其他成分很少.植物豐富度與溫帶成分有顯著正相關(guān)性,這說明該區(qū)種子植物地理成分類型符合該區(qū)所處地理位置為溫帶的特征.黃陵縣、洛川縣、黃龍縣、宜君縣、富縣和甘泉縣/值均在0.25以上,熱帶性質(zhì)高,這與其均處于暖溫帶落葉櫟林亞地帶相一致;志丹縣、吳起縣、定邊縣和華池縣/值均在0.25以下,熱帶性質(zhì)偏低,這與該區(qū)都處于溫帶相一致;而合水縣屬于暖溫帶北部落葉櫟林亞地帶,其/值卻低于志丹縣、吳起縣,這很有可能受到土壤、氣候或降水等因素的影響,這值得更深一步的研究.同時(shí),/值這種受海拔影響,隨溫度和降水大體上從東南向西北遞減的特征,能夠反映植物區(qū)系的地帶性格局,這種分布模式反應(yīng)出這些類群的生態(tài)位需求,可間接支持多樣性格局假說中的“寒冷忍耐假說”相關(guān)論斷,即冬季低溫是熱帶起源種北遷(或垂直向上遷移)的主要限制[41-43].這與其他學(xué)者的研究結(jié)果相一致[44-46],反映出黃土高原暖溫帶氣候—溫帶季風(fēng)性半濕潤氣候—半干旱氣候特點(diǎn).這表明本流域多種地理成分交錯(cuò),新老混雜,過渡特征明顯,廣布的溫帶性質(zhì)物種具有更好的適應(yīng)性,這種分布模式是植物長期適應(yīng)進(jìn)化的結(jié)果,體現(xiàn)出這些類群的生態(tài)位需求.

4 結(jié)論

4.1 黃土高原典型地貌植被區(qū)的北洛河流域有種子植物123科581屬1671種,裸子植物4科8屬13種,被子植物119科573屬1658種.大型科和較大科共20科337屬1095種是該區(qū)植物的主體,單種屬與小型屬共514屬938種,是該區(qū)植物的主體成分,植物起源古老,分化程度高,物種多樣性豐富,組成復(fù)雜,是重要的植物種質(zhì)資源庫.

4.2 該區(qū)植物相似性聚類可分為5組不同地理單元的植物組成,屬的分布區(qū)類型多樣,有15個(gè)分布區(qū)類型和16個(gè)變型,以溫帶分布為主,地理成分交錯(cuò),新老混雜.植物豐富度SD值與溫帶成分有顯著正相關(guān)性,/值受海拔影響,隨溫度和降水從東南向西北呈遞減的過渡性地帶特征.

4.3 該區(qū)進(jìn)行植被恢復(fù)和水土流失治理可從這6大優(yōu)勢科(菊科、禾本科、薔薇科、豆科、唇形科和毛茛科)中選擇適地適樹適草的溫帶性質(zhì)物種.

[1] 王荷生.植物區(qū)系地理[M]. 北京:科學(xué)出版社, 1992.

Wang H S. Floristic Geography [M]. Beijing:Science Press, 1992.

[2] 馬克平.試論生物多樣性的概念[J]. 生物多樣性, 1993,(1):20-22.

Ma K P. On the concept of biodiversity [J]. Biodiversity Science, 1993,(1):20-22.

[3] 吳征鎰,周浙昆,孫 航,等.種子植物分布區(qū)類型及其起源和分化[M]. 昆明:云南科技出版社, 2006.

Wu Z Y, Zhou Z K, Li D Z, et al. The areal-types of seed plant and their origin and differerntiation [M]. Kunming: Yunnan Science and Technology Press, 2006.

[4] 冉大川,劉 斌,王 宏,等.黃河中游典型支流水土保持措施減洪減沙作用研究[M]. 鄭州:黃河水利出版社, 2006.

Ran D C, Liu B, Wang H, et al. Study on flood and sediment reduction by soil and water conservation measures in typical tributaries of the Middle Yellow River [M]. Zhenzhou: Yellow River Water Conservancy Press, 2006.

[5] 陳 妮,李譚寶,張曉萍,等.北洛河流域植被覆蓋度時(shí)空變化的遙感動(dòng)態(tài)分析[J]. 水土保持通報(bào), 2013,33(3):206-210,327.

Chen N, Li T B, Zhang X P, et al. Spatiotemporal variations of vegetation coverage in Beiluo River watershed based on remote sensing data analysis [J]. Bulletin of Soil and Water Conservation, 2013,33(3):206-210,327.

[6] 唐克麗.黃河泥沙與黃土高原水土流失綜合治理問題[J]. 中國水土保持, 1985,(12):12-14.

Tang K L. Yellow River sediment and comprehensive control of soil and water loss in Loess Plateau [J]. Soil and Water Conservation in China, 1985,(12):12-14.

[7] 朱顯謨.黃土地區(qū)植被因素對于水土流失的影響[J]. 土壤學(xué)報(bào), 1960,(2):110-121.

Zhu X M. Influence of vegetation factors on soil erosion in loess area [J]. Acta pedologica Sinica, 1960,(2):110-121.

[8] 余新曉.森林植被減弱降雨侵蝕能量的數(shù)理分析[J]. 水土保持學(xué)報(bào), 1988,(2):24-30.

Yu X X. Mathematical analysis of forest vegetation reducing rainfall erosion energy [J]. Journal of Soil and Water Conservation, 1988,(2): 24-30.

[9] 焦菊英,王萬忠,李 靖,等.黃土丘陵溝壑區(qū)水土保持人工林減蝕效應(yīng)研究[J]. 林業(yè)科學(xué), 2002,38(5):87-94.

Jiao J Y, Wang W Z, Li J, et al. The soil erosion reducing benefit of soil and water conservation of plantation on the loess hilly and gully region [J]. Scientia Silvae Sinicae, 2002,38(5):87-94.

[10] 李占斌,朱冰冰,李 鵬.土壤侵蝕與水土保持研究進(jìn)展[J]. 土壤學(xué)報(bào), 2008,(5):802-809.

Li Z B, Zhu B B, Li P. Advancement in study on soil erosion and soil and water conservation [J]. Acta pedologica Sinica, 2008,(5):802-809.

[11] Ran L, Lu X, Xu J. Effects of vegetation restoration on Soil conservation and sediment loads in China: A critical review [J]. Critical Reviews in Environmental Science & Technology, 2013, 43(13):1384-1415.

[12] Chen H, Zhang X P, Muratjan A, et al. Effects of vegetation and rainfall types on surface runoff and soil erosion on steep slopes on the Loess Plateau, China [J]. CATENA, 2018,170:141-149.

[13] Wang Y, Feng X, Zhao J, et al. Quantifying the effects of vegetation restorations on the soil erosion export and nutrient loss on the Loess Plateau [J]. Frontiers in Plant Science, 2020,11.

[14] Zhang L, Dawes W R, Walker G R. Response of mean annual evapotranspiration to vegetation changes at catchment scale [J]. Water Resources Research, 2001,37(3):701-708.

[15] 徐憲立,馬克明,傅伯杰,等.植被與水土流失關(guān)系研究進(jìn)展[J]. 生態(tài)學(xué)報(bào), 2006,(9):3137-3143.

Xu X L, Ma K M, Fu B J, et al. Research review of the relationship between vegetation and soil loss [J]. Acta Ecologica Sinica, 2006,(9): 3137-3143.

[16] Makarieva A M , Gorshkov V G , Li B L. Precipitation on land versus distance from the ocean: Evidence for a forest pump of atmospheric moisture [J]. Ecological Complexity, 2009,6(3):302-307.

[17] 劉國彬,上官周平,姚文藝,等.黃土高原生態(tài)工程的生態(tài)成效[J]. 中國科學(xué)院院刊, 2017,32(1):11-19.

Liu G B, Shang G Z P, Yao W Y, et al. Ecological effects of soil conservation in Loss Plateau [J]. Bulletin of Chinese Academy of Sciences, 2017,32(1):11-19.

[18] Wu G L, Liu Y F, Zeng C, et al. Trade-off between vegetation type, soil erosion control and surface water in global semi-arid regions: A meta-analysis [J]. Journal of Applied Ecology, 2020,57(5).

[19] 張 琨,呂一河,傅伯杰,等.黃土高原植被覆蓋變化對生態(tài)系統(tǒng)服務(wù)影響及其閾值[J]. 地理學(xué)報(bào), 2020,75(5):949-960.

Zhang K, Lu Y H, Fu B J, et al. The effects of vegetation coverage changes on ecosystem service and their threshold in the Loess Plateau [J]. Acta Geographica Sinica, 2020,75(5):949-960.

[20] 李 軻,趙勇鋼,劉小芳,等.黃土丘陵區(qū)典型草地演替中植物群落特征與土壤儲水量關(guān)系[J]. 水土保持通報(bào), 2019,39(1):55-60.

Li K, Zhao Y G, Liu X F, et al.Relationship between plant community characteristics and soil water storage along typical grassland succession in loss hilly region [J].Bulletin of Soil and Water Conservation, 2019,39(1):55-60.

[21] 李登武.陜北黃土高原維管植物區(qū)系地理研究[D]. 楊凌:西北農(nóng)林科技大學(xué), 2008.

Li D W. The floristic geography of vascular plant in the Loess Plateau of Northern Shanxi Provinece [D]. Yangling: Northwest Agriculture and Forestry University, 2008.

[22] 王凱博.子午嶺植被演替過程中物種多樣性研究[D]. 楊凌:西北農(nóng)林科技大學(xué), 2008.

Wang K B, Studies on the species diversity in vegetation succession in ziwuling area [D]. Yangling: Northwest Agriculture and Forestry University, 2008.

[23] 方精云,朱江玲,郭兆迪,等.植物群落清查的主要內(nèi)容、方法和技術(shù)規(guī)范[J]. 生物多樣性, 2009,17(6):533-548.

Fang J Y, Zhu J L, Guo Z D, et al. Methods andprotocolsforplant community inventory [J]. Biodiversity Science, 2009,17(6):533-548.

[24] Qian H, Deng T, Jin Y, et al. Phylogenetic dispersion and diversity in regional assemblages of seed plants in China [J]. Proceedings of the National Academy of Sciences, 2019,116(46):23192-23201.

[25] 張鐿鋰,張雪梅.植物區(qū)系地理研究中的重要參數(shù)——相似性系數(shù)[J]. 地理研究, 1998,17(1):59-63.

Zhang Y L, Zhang X M. An important parameter in the study of floristic geography - Similarity coefficient [J]. Geographical Research, 1998,17(1):59-63.

[26] 白 龍,段博文,陳 曦,等.遼寧省西部低山丘陵區(qū)草地類型分布及植物區(qū)系特征[J]. 草地學(xué)報(bào), 2020,28(6):1726-1735.

Bai L, Duan B W, Chen X,etc. .Distribution of Grassland types and Floristic characteristics in the low mountains and hilly of Western Liaoning Province [J]. Acta Agrestia Sinica, 2020,28(6):1726-1735.

[27] 張光富.安徽板橋自然保護(hù)區(qū)植物多樣性[M]. 南京:南京師范大學(xué)出版社, 2001.

Zhang G F. Plant diversity in Banqiao Nature Reserve [M]. Nanjing: Nanjing Normal University Press, 2001.

[28] 董厚德.遼寧東部的主要植被類型及其分布[J]. 植物生態(tài)學(xué)與地植物學(xué)叢刊, 1981,(4):241-257.

Dong H D. Main vegetation types and their distribution in eastern Liaoning [J]. Series of Plant Ecology and Geobotany, 1981,(4):241- 257.

[29] 武吉華,張 坤,江 源,等.植物地理學(xué)[M]. 北京:高等教育出版社, 2004.

Wu J H, Zhang K, Jiang Y, et al. Phytogeography [M]. Beijing: Higher Education Press, 2004.

[30] 廖文波,張宏達(dá),仲銘錦.廣東亞熱帶植物區(qū)系表征科的區(qū)系地理學(xué)分析[J]. 生態(tài)科學(xué),1994,(1):47-55.

Liao W B, Zhang H D, Zhong M J.Floristic geographical analysis of the subtropical flora of Guangdong Province [J]. Ecological Science, 1994,(1):47-55.

[31] Cox C B, Moore P D.生物地理學(xué)[M]. 趙鐵橋譯.北京:高等教育出版社, 2007.

Cox C B, Moore P D. Biogeography [M]. Zhao Tieqiao trans. Beijing: Higher Education Press, 2007.

[32] 吳征鎰.中國種子植物屬的分布區(qū)類型[J]. 云南植物研究, 1991, (S4):1-139.

Wu Z Y. Distribution type of seed plant genera in China [J]. Yunnan Plant Research, 1991,(S4):1-139.

[33] 張文輝,李登武,劉國彬,等.黃土高原地區(qū)種子植物區(qū)系特征[J]. 植物研究, 2002,(3):373-379.

Zhang W H, Li D W, Liu G B,etc. Floristic characteristics of seed plants in the Loess Plateau [J]. Plant research, 2002,(3):373-379.

[34] 李錫文.中國種子植物區(qū)系統(tǒng)計(jì)分析[J]. 云南植物研究, 1996,18(4): 363-384.

Li X W. Floristic statistics and analyses of seed plants from China [J]. Yunnan Plant Research, 1996,18(4):363-384.

[35] 張文輝,康永祥,李 紅.陜西木本植物區(qū)系研究[J]. 植物研究, 1999,19(4):374-384.

Zhang W H, Kang Y X, Li H. Study on the woody plant flora of Shaanxi [J]. Plant Research, 1999,19(4):374-384.

[36] 李登武,張文輝,茍小東,等.陜北黃土高原植物區(qū)系的演變及新生代植物區(qū)系分析(英文) [J]. 草地學(xué)報(bào), 2011,19(2):181-190.

Li D W, Zhang W H, Gou X D. The evolution of the flora of the Loess Plateau in northern Shaanxi and the analysis of the flora of the Cenozoic [J]. Acta Agrestia Sinica, 2011,19(2):181-190.

[37] 左家哺.中國種子植物區(qū)系定量化研究V.區(qū)系相似性[J]. 熱帶亞熱帶植物學(xué)報(bào), 1996,(3):18-25.

Zuo J F. Quantitative study on seed plant flora of china v. Floristic similarity [J]. Journal of Tropical and Subtropical Botany, 1996, (3):18-25.

[38] 王義鳳,肖向明.黃土高原地區(qū)主要植被類型的氣候梯度分布[J]. 植物學(xué)報(bào), 1993,35(4):291-299.

Wang Y F, Xiao X M. Climatic gradiant of main vegetation types in the loess plateau region [J].Acta Botanica Sinica, 1993,35(4):291- 299.

[39] 張希彪,上官周平.黃土高原子午嶺種子植物區(qū)系特征研究[J]. 生態(tài)學(xué)雜志, 2005,24(8):872-877.

Zhang X B, Shang G Z P. Floristic characteristics of seed plants in Ziwuling of Loess Plateau [J]. Chinese Journal of Ecology, 2005,24(8): 872-877.

[40] 張振萬,張春林,陳彥生,等.陜北黃土高原植物區(qū)系分區(qū)界線初探[J]. 西北植物學(xué)報(bào), 1992,(4):53-58.

Zhang Z W, Zhang C L, Chen Y S. A Preliminary study on the border line of the Floristic regionalization in the loess plateau of northern shaanxi [J]. Acta BotanyNorthwest, 1992,(4):53-58.

[41] Currie D J, Francis A P. Regional versus climatic effect on taxon richness in angiosperms: Reply to Qian and Ricklefs [J]. The American Naturalist, 2004,163(5):780-785.

[42] 田義超,黃遠(yuǎn)林,張 強(qiáng),等.北部灣南流江流域土地覆蓋及生物多樣性模擬[J]. 中國環(huán)境科學(xué), 2020,v.40(3):394-408.

Tian Y C, Huang Y L, Zhang Q, et al. Simulation of land cover and biodiversity in Nanliu River Basin in Beibu Gulf [J]. China Environmental Science, 2020,v.40(3):394-408.

[43] 陳模舜,金則新,柯世省,等.極瀕危物種天臺鵝耳櫪群落特征及其與環(huán)境因子的關(guān)系[J]. 林業(yè)科學(xué), 2020,56(9):1-11.

Chen M S, Jin Z X, Ke S X, et al. Community characteristics and their relations with environmental variables of critically endangered species[J]. Scientia Silvae Sinicae, 2002,38(5):87-94.

[44] 唐志堯,方精云.植物物種多樣性的垂直分布格局[J]. 生物多樣性, 2004,(1):20-28.

Tang Z Y, Fang J Y. A review on the elevational patterns of plant species diversity [J]. Biodiversity Science, 2004,(1):20-28.

[45] 趙鳴飛,王宇航,邢開雄,等.黃土高原山地森林群落植物區(qū)系特征與地理格局[J]. 地理學(xué)報(bào), 2014,69(7):916-925.

Zhao M F, Wang Y H, Xing K X, et al. Floristic traits and geographic patterns of natural montane forests at community level on the Loess Plateau, North China [J]. Acta Geographica Sinica, 2014,69(7):916- 925.

[46] 董雪蕊,張 紅,張明罡.基于系統(tǒng)發(fā)育的黃土高原地區(qū)木本植物多樣性及特有性格局[J]. 生物多樣性, 2019,27(12):1269–1278.

Dong X Y, Zhang H, Zhang M G, et al. Explaining the diversity and endemic patterns based on phylogenetic approach for woody plants of the Loess Plateau [J]. Biodiversity Science, 2019,27(12):1269–1278.

Relationship between geographical pattern of plant diversity and environmental factors in Beiluo River Basin.

TIAN Qi-long1,2, XU Xiao-ming3, LYU Du1,2, WANG Hao-jia3, LEI Si-yue3, YI Hai-jie1,2, HE Jie3, HE Liang3, XUE Fan3, ZHOU Ya-dong3, WANG Miao-qian3, ZHANG Xiao-ping1,3*

(1.State Key Laboratory of Soil Erosion and Dryland Farming on Loess Plateau, Institute of Soil and Water Conservationg, Chinese Academy of Science and Ministry of Water Resources, Yangling 712100, China;2.University of Chinese Academy of Sciences, Beijing 10049, China;3.Institute of Soil and Water Conservation, Northwest A&F University, Yangling 712100, China)., 2021,41(9):4378~4387

The diversity of seed plants in Beiluo River Basin was identified through field investigation and statistical analysis of historical documents. Result: There were 1671 species of seed plants belonging to 581genera and 123 families in this area, accounting for 83.67%, 67.24% and 51.83% of the total seed plants in the Loess Plateau. of all the detected 1671 species, the majority (>65%) belonged to large (>50 species) and larger family (21~50 species), and more than half (~56%) belonged to single species genus and small genus (2~5 species). According to cluster analysis of the plant similarity for the 11main streams, 5groups were obtained, with Huanglong County as the firstly separated category, and the largest category including four geographical units, i.e. Fu County, Heshui County, Yijun County and Huangling County, which were all affected by the connection of Ziwuling mountains. The SD values of floristic richness of Fu County, Heshui County and Huanglong County were highest, while those of Luochuan County, Huachi County and Ganquan County were the lowest. There were 15 areal types and 16forms of geographical elements at the genus level in this area, which were dominated by the temperate distribution and its forms (287genera, accounting for 49.40%) and closely related to the tropical distribution and its forms (115 genera, accounting for 19.79%). There was a significant positive correlation between SD value of plant richness and temperate components. Huangling County, Luochuan County and Huanglong County had the highest/value and the strongest tropical property, while Huachi County and Dingbian County had the lowest/value and the weakest tropical property. The/value was affected by altitude and had a similar distribution pattern with temperature and precipitation decreasing from the southeast to the northwest. It is concluded that the plants in this area are important plant germplasm resources, because they have an ancient origin, a high degree of differentiation, rich species diversity and complex geographical elements. The transition pattern of plant species diversity zone is significant, which reflects the niche demands of these groups. Species of the six dominant families (Asteraceae, Poaceae, Rosaceae, Fabaceae, Lamiaceae and Ranunculaceae) that are suitable for the land and trees and grasses (temperate in nature) can be selected for vegetation restoration and soil erosion control.

species diversity;family and genus composition;geographical component;richness;temperate properties

X173

A

1000-6923(2021)09-4378-10

田起隆(1992-),男,陜西西安人,博士研究生,主要從事水土保持植被研究.發(fā)表論文2篇.

2021-01-26

國家自然科學(xué)基金資助項(xiàng)目(41877083,41440012)

* 責(zé)任作者, 研究員, zhangxp@ms.iswc.ac.cn

猜你喜歡
物種植物
物種大偵探
物種大偵探
吃光入侵物種真的是解決之道嗎?
英語世界(2023年10期)2023-11-17 09:18:18
回首2018,這些新物種值得關(guān)注
電咖再造新物種
汽車觀察(2018年10期)2018-11-06 07:05:26
植物的防身術(shù)
把植物做成藥
哦,不怕,不怕
將植物穿身上
植物罷工啦?
主站蜘蛛池模板: 久久性妇女精品免费| 免费在线成人网| 一本大道香蕉中文日本不卡高清二区| 亚洲乱伦视频| 国产乱子伦精品视频| 熟妇人妻无乱码中文字幕真矢织江 | 国产精品乱偷免费视频| 国产福利免费视频| 波多野结衣视频网站| 成人在线观看不卡| 久久久精品无码一区二区三区| 国产jizzjizz视频| 国产午夜无码专区喷水| 欧美成人国产| 四虎成人精品在永久免费| 极品国产在线| 黄色成年视频| 国产精品亚洲一区二区三区z | 国产情精品嫩草影院88av| 欧美黄色a| 日韩精品一区二区三区免费在线观看| 青青操视频在线| 日韩免费成人| 日本人妻丰满熟妇区| 欧美精品成人一区二区在线观看| 无码区日韩专区免费系列| 在线一级毛片| 91在线无码精品秘九色APP| 色国产视频| 亚洲色图综合在线| 国产情侣一区| 99久久精品视香蕉蕉| 无码AV日韩一二三区| 色噜噜狠狠狠综合曰曰曰| 国产成+人+综合+亚洲欧美| 手机精品福利在线观看| 久久国产精品影院| 亚洲高清在线播放| 四虎影视库国产精品一区| 精品国产亚洲人成在线| 国产一区二区丝袜高跟鞋| 夜夜拍夜夜爽| 亚洲欧美天堂网| 精品少妇三级亚洲| a毛片在线播放| 欧美日韩一区二区三区四区在线观看| 丁香五月婷婷激情基地| 亚洲欧洲日本在线| 色偷偷男人的天堂亚洲av| 久久美女精品| 99热这里只有精品2| 免费中文字幕一级毛片| 狠狠色香婷婷久久亚洲精品| 正在播放久久| av免费在线观看美女叉开腿| 久久这里只精品热免费99| 狠狠操夜夜爽| 玖玖精品视频在线观看| 老司国产精品视频91| 国产黑人在线| 亚洲精品视频网| 日本福利视频网站| 欧美三级不卡在线观看视频| 狠狠色成人综合首页| 88av在线看| 久久综合伊人 六十路| 99久久人妻精品免费二区| 日韩欧美中文亚洲高清在线| 青青草原国产免费av观看| 成年人视频一区二区| 十八禁美女裸体网站| 国内黄色精品| 国产美女主播一级成人毛片| 美女视频黄频a免费高清不卡| 天天综合天天综合| 国产精品亚洲αv天堂无码| 国产男人的天堂| 最新午夜男女福利片视频| 波多野结衣视频一区二区| 成人在线观看一区| 欧美精品另类| 日韩精品一区二区三区swag|