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

    基于國家標(biāo)本資源共享平臺數(shù)據(jù)的中國杜鵑花科植物標(biāo)本采集地理偏差分析

    2018-05-20 06:56:48張玉雪任敏肖翠
    關(guān)鍵詞:氣候因子杜鵑花差值

    張玉雪,任敏,肖翠

    1.中國科學(xué)院遺傳與發(fā)育所農(nóng)業(yè)資源研究中心,河北 石家莊 050021

    2.安徽農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,安徽 合肥 230036

    3.中國科學(xué)院植物研究所,北京 100093

    1 背景

    1.1 研究背景

    杜鵑花科 (Ericaceae) 有 967 種,多為落葉灌木、半常綠灌木、常綠灌木、常綠小喬木,分布遍及歐洲、亞洲、北美洲和大洋洲,其中亞洲是杜鵑花科植物的集中分布地區(qū)[1]。據(jù)記載,杜鵑花科植物起源時間約為晚白堊紀(jì)至早第三紀(jì)的過渡期,中國西南至中國中部最有可能是杜鵑花屬植物的起源地[2]。我國有杜鵑花科 15 屬,約 757 種,分布于全國各地,主產(chǎn)地在西南部山區(qū),尤以西藏、四川、云南三省區(qū)相鄰地區(qū)為盛[3]。目前對于杜鵑花的資源保護(hù)利用并不成熟,由于人為的過度采伐促使野生杜鵑花資源受到嚴(yán)重破壞[4]。

    杜鵑花科植物在我國西藏、四川、云南等地區(qū)分布廣泛,在區(qū)域或全國尺度已進(jìn)行了對杜鵑花屬區(qū)系組成、地理分布格局的研究。青藏高原的川西地區(qū)與藏東南地區(qū)杜鵑花屬植物類群及其分布格局的異同,揭示了我國杜鵑花屬植物地理區(qū)系及其垂直分布特征由東至西的趨異化過程與性質(zhì)[5]。莊平分析探討了我國杜鵑花屬植物地理分布格局,喜瑪拉雅山、狹義橫斷山、川西山地既是杜鵑花屬的集中分布區(qū),也是某些類群擴(kuò)散、遷移的地理屏障,對我國現(xiàn)代杜鵑花分布區(qū)及其分布格局的形成影響重大[6]。

    物種豐富度是生物多樣性最基本的衡量指標(biāo)之一,目前關(guān)于物種豐富度格局的研究先后有學(xué)者提出了氣候、地質(zhì)歷史過程、環(huán)境異質(zhì)性以及隨機(jī)過程等眾多假說機(jī)制[7-13]。盡管物種豐富度格局受到多種環(huán)境因子的共同影響,但氣候?qū)ξ锓N豐富度格局的影響被認(rèn)為是最主要的決定因子[14-15]。其中水熱動態(tài)假說最受關(guān)注,其認(rèn)為水分和能量共同決定植物物種豐富度格局[16-17],該觀點(diǎn)在中國黑戈壁[18]、內(nèi)蒙古草原[19]等地區(qū)研究中已得到驗(yàn)證。另外,氣候穩(wěn)定性假說也是討論較多的,該假說認(rèn)為氣候穩(wěn)定性越高的地區(qū),物種分布區(qū)受氣候變化影響越低,動植物的種類越豐富[20]。

    1.2 研究意義

    植物標(biāo)本是植物學(xué)家長期從事科研活動的積累和人類自然遺產(chǎn)的永久記錄之一,是研究物種的分布及其歷史、現(xiàn)狀、系統(tǒng)演化的重要資料。針對標(biāo)本數(shù)據(jù)的地理采集偏差檢驗(yàn)和評估的研究有部分學(xué)者已經(jīng)著手,如張大才等收集了有關(guān)數(shù)據(jù)庫和標(biāo)本館的標(biāo)本資料分析了橫斷山區(qū)樹線以上區(qū)域種子植物標(biāo)本的采集現(xiàn)狀和物種豐富度[21],前人對杜鵑花科植物多是在杜鵑花屬水平以下進(jìn)行其多樣性、區(qū)系與地理分布的研究[2,5,6],但是基于國家標(biāo)本資源共享平臺 (NSII) 的數(shù)據(jù)對中國杜鵑花科植物標(biāo)本采集偏差的分析很少且在杜鵑花科水平上很少見到關(guān)于其多樣性空間分布格局與氣候因子關(guān)系的研究。本研究擬以國家標(biāo)本資源共享平臺收錄的中國杜鵑花科植物的植物標(biāo)本采集記錄為研究對象,選取縣級行政單位為地理單元,進(jìn)行我國杜鵑花科植物采集現(xiàn)狀以及偏差分析;通過對該科植物地理分布與氣候因子關(guān)系的探討,結(jié)合 ArcGIS10.2 軟件和 SPSS 軟件開展杜鵑花科植物在我國的適生區(qū)分布預(yù)測,擬探討以下問題:(1) 各氣候因子變量對杜鵑花科物種豐富度的影響;(2) 水分、能量對杜鵑花科物種豐富度的影響。本研究結(jié)果為今后開展杜鵑花科植物區(qū)系的調(diào)查及生物多樣性的保護(hù)和資源開發(fā)利用提供理論參考。

    2 材料與方法

    2.1 數(shù)據(jù)的篩選與校對

    本研究中選取的杜鵑花科植物數(shù)字化標(biāo)本信息全部來源于國家標(biāo)本資源共享平臺 (NSII) (National Specimen Information Infrastructure,” (http://www.nsii.org.cn/)。經(jīng)初步篩選后得到 149,559 份杜鵑花科的標(biāo)本信息,為進(jìn)一步提高標(biāo)本數(shù)據(jù)質(zhì)量,現(xiàn)對標(biāo)本數(shù)據(jù)進(jìn)行以下的校對工作。

    步驟如下:(1) 植物拉丁名稱校對。利用 (http://tnrs.iplantcollaborative.org/TN RSapp.html 拉丁文名稱核對網(wǎng)站對標(biāo)本拉丁名稱的標(biāo)準(zhǔn)化與錯誤數(shù)據(jù)進(jìn)行剔除,對標(biāo)本拉丁名稱進(jìn)行查詢,將準(zhǔn)確率得分在 0.9 以下的標(biāo)本數(shù)據(jù)剔除,準(zhǔn)確率得分在 0.9 以上的標(biāo)本數(shù)據(jù)接收為正確的拉丁名稱。(2) 剔除數(shù)據(jù)缺失的采集地點(diǎn) (香港、澳門特行政區(qū)、臺灣省以及南海諸島數(shù)據(jù)不全) 以及在中國行政區(qū)以外的采集記錄。(3) 對于采集地點(diǎn)在自然保護(hù)區(qū)內(nèi)的標(biāo)本記錄,根據(jù)電子地圖查詢自然保護(hù)區(qū)所在地,若自然保護(hù)區(qū)完整地在某個縣級地區(qū)內(nèi),將保護(hù)區(qū)名更改為該縣名。(4) 縣級行政單位名稱校對。根據(jù)中國縣級行政區(qū)劃地圖確定縣名。部分舊縣名經(jīng)過行政區(qū)劃的合并和重新劃分,在查找其歷史沿革后,確定最新的縣名。最終經(jīng)過整理后共獲取了約 124, 940 條地標(biāo)化至縣級的杜鵑花科植物標(biāo)本數(shù)據(jù),共包含有杜鵑花科植物 1611 種。

    2.2 采集完整程度的評估方法

    標(biāo)本采集偏差指在植物標(biāo)本采集過程中受地理?xiàng)l件,天氣,交通,采集人的主觀意識等條件的影響導(dǎo)致采集時采集地區(qū)分布不均勻,部分地區(qū)采集遺漏,或采集數(shù)量過少導(dǎo)致數(shù)據(jù)不完整。本研究共使用兩種方法來評估各個縣杜鵑花科植物在該地的標(biāo)本采集完整程度。第一種為差值法,第二種為斜率法。

    2.2.1 差值評估法

    該方法是由 Chao 提出了一種基于評估稀有種數(shù)量的群落豐富度估算方法,常稱為 Chao1 估計(jì)量。Chao1 的估計(jì)以已有取樣中物種個體數(shù)量的分布頻率為基礎(chǔ),當(dāng)取樣中每個物種都是兩個或兩個以上的個體時,Chao1 認(rèn)為該地區(qū)的取樣比較完全;而當(dāng)取樣中每個物種只有一個個體時,Chao1 認(rèn)為該地區(qū)仍有新物種有待發(fā)現(xiàn)。

    其計(jì)算方程為:

    其中,Sobs 是采集到的物種數(shù),a 是取樣中只有一個個體的物種數(shù),b 是取樣中只有兩個個體的物種數(shù),S*是 Chao1 估計(jì)的物種數(shù)[11]。根據(jù)該公式可得知每個縣中稀有種數(shù)量越多,該地區(qū)新物種未被發(fā)現(xiàn)的概率越大,估計(jì)豐富度越高;而每個縣中稀有種數(shù)量越少,該地區(qū)新物種未被發(fā)現(xiàn)的概率越小,估計(jì)豐富度越小。由于稀有種的存在和野外調(diào)查的限制,植物群落中記錄在案的物種數(shù)量一般少于真實(shí)的物種數(shù)量。并經(jīng) Shen 等和 Chao 等證實(shí),對于大部分?jǐn)?shù)據(jù)集來說,Chao1 估計(jì)量是一種計(jì)算物種豐富度的良好指標(biāo),陽文靜使用 Jackknife 和 Bootstrap 等模型計(jì)算出的估計(jì)結(jié)果和 Chao1 估計(jì)量的結(jié)果基本相似,說明 Chao1 估計(jì)量具科學(xué)性[22]。

    本研究選取物種豐富度 Chao1 估計(jì)量來計(jì)算中國各個縣的杜鵑花科植物物種豐富度,并用 Chao1 估計(jì)量與記錄在案的物種豐富度之間的差值來評估一個縣的采集完整程度。在評估時,本研究中使用 Spade 軟件進(jìn)行計(jì)算,并在計(jì)算過程中本研究將稀有種個數(shù)定為 10。最終計(jì)算結(jié)果差值越大表明稀有種出現(xiàn)的概率越大,該縣中稀有種數(shù)量越多,新物種未被發(fā)現(xiàn)的概率越大;反之越小。

    2.2.2 斜率評估法

    該方法是由李巧等提出的利用物種累積曲線 (species accumulate on curve, SAC) 末端百分之十的斜率的大小作為標(biāo)本采集完整程度的標(biāo)準(zhǔn)[23]。物種累積曲線用于描述隨機(jī)抽取采集記錄情況下,隨著抽樣強(qiáng)度增加,物種豐富度隨之積累的情況,是理解調(diào)查樣地物種組成和預(yù)測物種豐富度的有效工具,在生物多樣性和群落調(diào)查中,被廣泛用于抽樣量充分性的判斷。

    本研究假設(shè),對于每一個縣級行政單位,其取樣完整度越低對應(yīng)的物種累積曲線越趨近于直線,而取樣完整度越高對應(yīng)的物種累積曲線的彎曲程度越高。其中物種累積曲線的彎曲程度體現(xiàn)了物種的不完整程度,物種積累曲線的彎曲程度可以用物種積累曲線尾端 10% 的斜率表示[24]。

    本研究中使用“R”軟件計(jì)算物種累積曲線尾端 10% 部分的斜率,表示物種各積累曲線的彎曲程度,代表某個縣采集數(shù)據(jù)的不完整程度。當(dāng)斜率值趨近于 0 時,說明取樣較為完整;當(dāng)斜率近于 1 時,說明取樣較為不完整。該斜率值還表示如果采集強(qiáng)度增加,將會有新物種被發(fā)現(xiàn)。

    2.3 氣候因子變量數(shù)據(jù)篩選

    19 個氣候因子從世界氣候數(shù)據(jù)庫下載 (http://www.world-climates.com),并將其進(jìn)一步歸類為水分因子和能量因子。應(yīng)用 ArcGIS10.2 獲得與中國分布區(qū)相關(guān)的 19 個氣候因子數(shù)據(jù)變量,GIS 分析底層是杜鵑花科標(biāo)本匹配的每個縣經(jīng)緯度數(shù)據(jù) (表 1)。

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

    得到氣候因子數(shù)據(jù)及杜鵑花科植物的 3 個物種豐富度數(shù)據(jù),在 SPSS 22.0 進(jìn)行下一步分析,包括:

    (1) 氣候因子的 Spearman 秩相關(guān)檢驗(yàn)

    影響杜鵑花科分布的氣候因子有 19 個,但由于變量的多重共線性會導(dǎo)致物種分布模型過度擬合[23-24],所以需要將高度相關(guān)的氣候因子剔除掉。利用 Spearman 秩相關(guān)檢驗(yàn)在 spss 22.0 中篩選出相關(guān)性最小的變量 (Spearman's Rho b 0.75),相關(guān)變量的 Spearman Rho 大于 0.75 時只保留一個生態(tài)邏輯上最有意義的氣候因子。

    表1 影響杜鵑花科植物空間分布的 19 個氣候因子Table 1 19 climatic factors affecting the spatial distribution of rhododendron plants

    (2) 逐步回歸分析

    將氣候因子分為能量和水分兩組,通過逐步回歸分析篩選出對 3 個物種豐富度最重要的氣候因子變量,同時記錄每組的相關(guān)性指數(shù) R square。逐步回歸的思路是將物種豐富度作為因變量,將氣候因子作為自變量逐個添加到回歸模型中,每添加一個解釋變量后都要進(jìn)行 F 檢驗(yàn),并對已經(jīng)選入的解釋變量逐個進(jìn)行 t 檢驗(yàn),如果原來引入的解釋變量由于后面解釋變量的加入變得不再顯著時,則將其刪除。以確保每次加入新的變量之前回歸方程中只包含顯著性變量。

    (3) 一元線性回歸分析

    探究不同氣候因子變量與杜鵑花科物種豐富度之間關(guān)系。

    (4) 分析氣候因子、物種豐富度與緯度之間相關(guān)性

    通過顯著性水平檢驗(yàn)發(fā)現(xiàn)氣候因子、物種豐富度與緯度之間用一元二次回歸解釋時相關(guān)性較高,所以利用一元二次回歸分析來探討氣候因子、物種豐富度與緯度之間相關(guān)性。

    3 結(jié)果與分析

    3.1 我國杜鵑花科植物標(biāo)本采集現(xiàn)狀

    從采集數(shù)目上來看,杜鵑花科植物標(biāo)本采集主要集中于長江以南地區(qū),尤其西南地區(qū)分布最多,西北與東北地區(qū)分布較少。四川最多,云南次之,其標(biāo)本數(shù)據(jù)均過萬,第三為西藏,占我國杜鵑花科植物標(biāo)本數(shù)目一半以上,見表 1。

    由圖一 A 可以看出長江以北地區(qū),如黑龍江省、河北省、山西省、山東省、江蘇省等省份的標(biāo)本采集地區(qū)零星分布,其中標(biāo)本采集數(shù)目最少的是山東省,僅有四份。黑龍江共 273 分標(biāo)本,主要集中在大興安嶺地區(qū)的三個縣:呼瑪縣、漠河縣、塔河縣;河北省標(biāo)本共 361 份,主要集中于張家口市內(nèi)恒山、太行山、燕山三山交匯之處。江蘇省僅 239 份標(biāo)本,主要集中于無錫市與宜興市多丘陵地區(qū)。綜上,杜鵑花科植物標(biāo)本集中于我國西南地區(qū):青藏高原東南部、橫斷山區(qū)、四川平原、云貴高原、大巴山區(qū)等,該地區(qū)有著豐富多樣的小生境類型,是我國的植物特有種中心,稀有種種類豐富,四川盆地的標(biāo)本分布較少,可能是由于其熱量比周邊山區(qū)高,而降雨則低于周邊山區(qū),不適宜杜鵑花科植物生長,且耕地連片,植被破壞較嚴(yán)重,故分布較少。其中分子生物學(xué)證據(jù)也表示,這些地區(qū)地形的復(fù)雜性造成了植物種間隔離,加速了物種的分化。

    3.2 物種豐富度估計(jì)值與觀察值的差值分析

    由圖 1B、C 中可以看出物種豐富度記錄值和 Chao1 估計(jì)值最大的均為云南省的貢山獨(dú)龍族怒族自治縣、德欽縣和西藏自治區(qū)的墨脫縣。物種豐富度記錄值較高的地點(diǎn)多集中于橫斷山區(qū)、云貴高原、青藏高原西部以及四川盆地。物種豐富度記錄值低的地區(qū)主要是新疆維吾爾族自治區(qū)、藏北高原、內(nèi)蒙古中部、大興安嶺山區(qū)、巴顏喀拉山脈、云貴高原東北部。同一個省內(nèi),各縣采集強(qiáng)度參差不齊。

    圖1 A.標(biāo)本采集數(shù)量;B.原始豐富度 (物種數(shù));C.豐富度估計(jì)值;D.殘差 (豐富度估計(jì)值-實(shí)際值);E.物種積累曲線尾端 10% 斜率Fig.1 A: Number of collected specimens; B: Original richness; C: Richness estimates; D: Richness residue (estimated value minus actual value); E: Gradient of the last 10% section of the accumulated species curve

    圖 1D 差值圖中可直觀看出差值的范圍是 0 至 84 種。差值的大小代表采集的完整程度,現(xiàn)將如圖所示的差值分為五個區(qū)間:0-3.7,3.7-9.4,9.4-18.9,18.9-33.0,33.0-84.2,并定義五個區(qū)間分別代表:極完整、較完整、中等完整、較不完整、極不完整。結(jié)果表明標(biāo)本數(shù)目較多的橫斷山區(qū)有存在采集不完整的縣,差值法得出的極不完整的縣如表 4。

    根據(jù)差值法,可由圖二看出,采集完整程度與采集數(shù)目并沒有較大關(guān)聯(lián)性,如采集數(shù)目最多的峨眉山市 Chao1 估計(jì)值與物種豐富度記錄值的差值為 34.6,屬于采集極不完整的縣;仁化縣標(biāo)本采集數(shù)僅 72 份,物種豐富度觀察值為 15,差值僅 0.9,屬于采集極完整的縣。說明差值法在一定程度上反映了采集完整程度,但其可靠性受標(biāo)本基數(shù)、稀有種數(shù)目的設(shè)定影響。

    3.3 物種積累曲線尾端 10% 斜率分析

    以 0.05 為閾值,在有標(biāo)本采集記錄的縣級行政單位中,標(biāo)本采集不完整程度指數(shù) (SAC slope) 在 0.012500-0.818182 之間變化,有 73 個縣 (全部縣的 0.12%) 的斜率值 (SAC slope) 小于或等于 0.05,可以認(rèn)為是名錄較為完整的縣。根據(jù)斜率法計(jì)算得出完整性較高的縣主要集中于橫斷山區(qū)、云貴高原、藏南谷地、大興安嶺山區(qū)以及大巴山區(qū),與標(biāo)本采集數(shù)目多的地區(qū)相符。對于采集完整的縣,Chao1 估計(jì)值與采集豐富度的差值取值范圍是 0-40.8,而采集不完整的縣,Chao1 估計(jì)值與采集豐富度的差值取值范圍是 0-84.2,差值的取值范圍縮??;斜率值最高的二十個縣標(biāo)本采集數(shù)量均低于 50 份 (表 5),以云縣為例,云縣位于云南省西南地區(qū),地處橫斷山系縱谷區(qū)南部,地貌類型豐富,為亞熱帶季風(fēng)氣候,自然資源豐富,適宜杜鵑花科植物生長,但標(biāo)本數(shù)量只有 15 份,這說明該縣確實(shí)采集不完整,需加強(qiáng)標(biāo)本采集工作。故經(jīng)過斜率法對標(biāo)本進(jìn)行質(zhì)量評估后,使用采集完整性高的數(shù)據(jù)在一定程度上可以提高數(shù)據(jù)分析的科學(xué)性。

    表2 Spear man 秩相關(guān)檢驗(yàn)出的氣候因子變量Table 2 19 climatic factors affecting the spatial distribution of rhododendron plants

    表3 標(biāo)本分布總數(shù)及比例Table 3 Distribution of specimen and the proportion

    3.4 氣候因子對杜鵑花科物種豐富度的影響

    由 Spearman 秩相關(guān)檢驗(yàn) (表 2) 以及一元線性回歸結(jié)果顯示,影響杜鵑花科物種豐富度的氣候因子主要包括晝夜溫差與年溫差比值 (bio3)、年溫度變化范圍 (bio7)、年平均雨量 (bio12) 和最暖季度平均雨量 (bio18) 等 4 個氣候因子 (圖 3)。研究分析發(fā)現(xiàn),觀察物種豐富度、估計(jì)物種豐富度、估計(jì)與觀察物種豐富度的差值與年溫度變化范圍 (bio7) 均呈顯著的負(fù)相關(guān),即隨著年溫度變化范圍的升高,3 個物種豐富度呈遞減趨勢 (圖 3B,F(xiàn),J) ;而 3 個物種豐富度與晝夜溫差與年溫差比值 (bio3)、年平均雨量 (bio12)、最暖季度平均雨量 (bio18) 等氣候因子均呈顯著的正相關(guān) (圖 3A,C,D,E,G,H,I,K,L)。其中年溫度變化范圍 (bio7) 對 3 個物種豐富度提供的解釋率最高 (圖 3B,F(xiàn),J),這說明反映溫度季節(jié)性分異 (氣候穩(wěn)定性) 的年溫度變化范圍 (bio7) 在很大程度上決定了杜鵑花科物種豐富度的分布格局。通過估計(jì)與觀察物種豐富度的差值與 4 個氣候因子的線性關(guān)系可以看出,在雨水充足、氣溫穩(wěn)定且較高的生境中杜鵑花科新物種未被發(fā)現(xiàn)的概率較大 (圖 3I,J,K,L),綜合可以看出杜鵑花科喜愛雨水充足、氣溫穩(wěn)定且較高的氣候環(huán)境。

    表4 極不完整縣 (差值法)Table 4 Counties of incomplete collection (difference method)

    圖2 差值與標(biāo)本采集數(shù)目的相關(guān)性Fig.2 Relevance of differentials and specimen quantities

    圖3 杜鵑花科物種豐富度與氣候因子的關(guān)系Fig.3 Relevance of rhododendron species richness and climatic factors

    圖3 杜鵑花科物種豐富度與氣候因子的關(guān)系Fig.3 Relevance of rhododendron species richness and climatic factors

    表5 斜率最高的 20 個縣及其標(biāo)本數(shù)目Table 5 20 counties with the highest gradients and respective specimen quantities

    表6 能量與水分對杜鵑花科物種豐富度的解釋 (R square)Table 6 Implication of energy and moisture on rhododendron species richness (R square)

    3.5 能量、水分對杜鵑花科物種豐富度的影響

    多元逐步回歸分析表明,進(jìn)入回歸方程的能量氣候因子為晝夜溫差與年溫差比值 (bio3)、年溫度變化范圍 (bio7) (表 6)。由分析可知,能量因子對 3 個物種豐富度的解釋率均高于水分因子,并且能量、水分對 3 個物種豐富度的單獨(dú)解釋率均小于二者共同作用對 3 個物種豐富度的解釋率,表明相對于水分因子對杜鵑花科物種豐富度的影響,能量因子對杜鵑花科物種豐富度的影響更為顯著,并且能量水分共同作用的影響最為顯著。

    4 結(jié)論與討論

    經(jīng)以上兩種評估方法,發(fā)現(xiàn)我國杜鵑花科地理采集偏差十分嚴(yán)重,標(biāo)本采集工作仍需更進(jìn)。

    植物標(biāo)本的采集工作是一個長期的歷史過程。我國國土面積遼闊,植被類型多樣,生境類型復(fù)雜,地形地貌復(fù)雜多變,系統(tǒng)化的進(jìn)行標(biāo)本采集幾乎難以實(shí)現(xiàn)。杜鵑花科植物物種豐富度與標(biāo)本采集數(shù)目最高的區(qū)域位于橫斷山區(qū)。橫斷山脈地跨川滇藏三省,獨(dú)特的地貌和高海拔低緯度的特點(diǎn),使得該地區(qū)自然條件獨(dú)具一格,生物區(qū)系絢麗多彩,如轎子山自然保護(hù)區(qū)、高黎貢山自然保護(hù)區(qū)、點(diǎn)蒼山自然保護(hù)區(qū)等[25]。杜鵑花科植物分為落葉與常綠兩大類,常綠類居多且喜高山深山。故高山及亞高山是杜鵑花科植物的分布中心。

    標(biāo)本采集的完整性使用兩種方法進(jìn)行評估,本研究中兩種評估方法得出的結(jié)果相似性一致 (圖 1D;圖 1E)。僅有不到1%的縣級行政單位被評估為采樣完全,表明未來野外調(diào)查的任務(wù)仍然十分艱巨。我國西南三省是世界杜鵑花科植物的發(fā)源地與現(xiàn)代分布中心,特別是常綠種類,由于未合理開發(fā)利用以及管理不當(dāng),我國對于杜鵑花科植物引種馴化程度不高,甚至有部分高山資源遭到破壞[4]。因此,針對專科專屬的普查十分重要,并對于建立杜鵑花科自然保護(hù)區(qū)有重要意義。

    根據(jù)滇西北縣域生物多樣性本底調(diào)查與評估,盡管中國已在生物多樣性調(diào)查方面做了大量的工作,也初步在大尺度上掌握了其基本組成。但是中國生物多樣性本底調(diào)查工作仍十分薄弱,物種分布數(shù)據(jù)掌握不夠詳盡,如部分縣級區(qū)的物種分布數(shù)據(jù)僅與物種數(shù)相當(dāng),館藏標(biāo)本的代表性較差,各地區(qū)、各類群的采集程度十分不均衡,而取樣的偏差直接影響到生物多樣性豐富度統(tǒng)計(jì)的準(zhǔn)確性[26-27]。如今越來越多的學(xué)者在使用標(biāo)本數(shù)據(jù)前進(jìn)行質(zhì)量評估,如覃海寧,楊永在中國高等植物受威脅物種名錄的確立過程中就進(jìn)行了標(biāo)本質(zhì)量評估[28]。

    研究結(jié)果表明,在影響杜鵑花科 3 個物種豐富度的 4 個氣候因子中,年溫度變化范圍 (bio7) 顯著影響物種豐富度 (圖 3)。如上所述,年溫度變化范圍 (bio7) 表示的是研究縣內(nèi)溫度的季節(jié)性分異,即該縣氣候的穩(wěn)定性。根據(jù)氣候穩(wěn)定性假說,在氣候穩(wěn)定性越高的地區(qū),物種分布區(qū)受氣候變化影響越低,動植物的種類越豐富[20],因此年溫度變化范圍小的地區(qū)杜鵑花科物種豐富度可能高于年溫度變化范圍大的地區(qū)。根據(jù) Rapoport 法則,從南到北隨著緯度的升高,年溫度變化范圍 (bio7) 呈上升趨勢,所以物種豐富度呈下降趨勢[29],這在一定程度上解釋了本研究中杜鵑花科物種豐富度從南至北的遞減趨勢。在神農(nóng)架、云南獨(dú)龍江等地區(qū)植物物種垂直分布格局的研究中,Rapoport 法則均得到驗(yàn)證[30-31]。在研究中國種子植物物種多樣性大尺度分布格局中也發(fā)現(xiàn),年溫度變化范圍在一定程度上決定了物種豐富度的空間分布格局[31]。

    水熱動態(tài)假說在研究影響區(qū)域植物物種豐富度格局中最受關(guān)注,該假說認(rèn)為水分和能量共同決定植物物種豐富度分布格局[16-17]。該觀點(diǎn)在中國黑戈壁[18]、內(nèi)蒙古草原[19]等地區(qū)研究中已得到驗(yàn)證。水分不僅是植物生物化學(xué)反應(yīng)過程中不可缺少的溶劑,更是光合作用中光反應(yīng)階段的關(guān)鍵反應(yīng)物,水分和能量共同控制著植物的生長,進(jìn)而影響植物多樣性格局[19]。本文對能量、水分與物種豐富度的回歸結(jié)果顯示,能量和水分共同作用對杜鵑花科 3 個物種豐富度的解釋率均高于能量、水分單獨(dú)對杜鵑花科 3 個物種豐富度的解釋率,這表明了能量、水分對 3 個物種豐富度的單獨(dú)解釋力較弱,唯有二者共同作用才能對杜鵑花科物種豐富度有著強(qiáng)有力的解釋力,該結(jié)果與上述研究結(jié)果一致。

    猜你喜歡
    氣候因子杜鵑花差值
    青海省果洛州高原鼠兔宜生區(qū)分布及其對氣候變化的響應(yīng)
    杜鵑花開
    差值法巧求剛體轉(zhuǎn)動慣量
    待到杜鵑花開時
    心聲歌刊(2020年4期)2020-09-07 06:37:14
    杜鵑花
    哦,杜鵑花!
    枳殼及其炮制品色差值與化學(xué)成分的相關(guān)性
    中成藥(2017年6期)2017-06-13 07:30:35
    四川水稻產(chǎn)量及其構(gòu)成要素對不同生育期氣候因子的響應(yīng)分析
    天山北坡NDVI對氣候因子響應(yīng)的敏感性分析
    基于區(qū)域最大值與平均值差值的動態(tài)背光調(diào)整
    中亚洲国语对白在线视频| 亚洲一区二区三区色噜噜| 欧美乱色亚洲激情| 色综合欧美亚洲国产小说| 国产真实乱freesex| 可以在线观看的亚洲视频| 熟女人妻精品中文字幕| 欧美zozozo另类| 久久九九热精品免费| 亚洲av免费高清在线观看| bbb黄色大片| 欧美3d第一页| 亚洲av电影不卡..在线观看| 亚洲精品456在线播放app | 国产精品一区二区三区四区久久| 国产一区二区在线av高清观看| 亚洲激情在线av| 琪琪午夜伦伦电影理论片6080| 听说在线观看完整版免费高清| 中亚洲国语对白在线视频| 日本一二三区视频观看| 99久久精品国产亚洲精品| www.熟女人妻精品国产| 内地一区二区视频在线| 亚洲三级黄色毛片| 日韩精品中文字幕看吧| 久久久久免费精品人妻一区二区| 欧美色欧美亚洲另类二区| 国产精品久久久久久人妻精品电影| 亚洲成av人片免费观看| 亚洲五月婷婷丁香| 99热这里只有是精品在线观看 | 人妻夜夜爽99麻豆av| 人人妻人人澡欧美一区二区| 男人舔女人下体高潮全视频| 欧美一区二区国产精品久久精品| 亚洲七黄色美女视频| 色精品久久人妻99蜜桃| 1024手机看黄色片| 日韩 亚洲 欧美在线| 精品午夜福利视频在线观看一区| 日日摸夜夜添夜夜添小说| 在线国产一区二区在线| 国产视频一区二区在线看| 日本黄大片高清| 1000部很黄的大片| www日本黄色视频网| 国产精品av视频在线免费观看| 十八禁网站免费在线| 久久6这里有精品| 51国产日韩欧美| 精品午夜福利视频在线观看一区| 国产成人a区在线观看| 99久久99久久久精品蜜桃| 亚洲色图av天堂| 此物有八面人人有两片| 99久久精品热视频| 舔av片在线| 男插女下体视频免费在线播放| 精品一区二区三区视频在线| 欧美丝袜亚洲另类 | 欧美潮喷喷水| 亚洲经典国产精华液单 | 午夜福利18| 婷婷色综合大香蕉| 亚洲18禁久久av| 两个人视频免费观看高清| av女优亚洲男人天堂| 国产综合懂色| 色5月婷婷丁香| 黄色丝袜av网址大全| 18+在线观看网站| 熟女电影av网| 亚洲国产高清在线一区二区三| 亚洲精品亚洲一区二区| 一级a爱片免费观看的视频| 在线观看一区二区三区| 男人舔女人下体高潮全视频| 观看美女的网站| 亚洲欧美清纯卡通| 女人被狂操c到高潮| 欧美最黄视频在线播放免费| 欧美高清成人免费视频www| 激情在线观看视频在线高清| 人人妻,人人澡人人爽秒播| 草草在线视频免费看| a级毛片a级免费在线| a级毛片免费高清观看在线播放| 亚洲无线在线观看| 天天躁日日操中文字幕| 中国美女看黄片| 中出人妻视频一区二区| 国产亚洲精品av在线| 人妻夜夜爽99麻豆av| 日本 av在线| 国产一区二区三区在线臀色熟女| 国产精品久久视频播放| 桃色一区二区三区在线观看| 免费av观看视频| 人妻夜夜爽99麻豆av| 亚洲av美国av| 国产伦人伦偷精品视频| 男插女下体视频免费在线播放| 婷婷丁香在线五月| 亚洲av中文字字幕乱码综合| 午夜福利在线观看免费完整高清在 | 美女cb高潮喷水在线观看| 婷婷六月久久综合丁香| 国产成年人精品一区二区| 日韩免费av在线播放| 国产高清激情床上av| 国产黄色小视频在线观看| 午夜激情福利司机影院| 在线国产一区二区在线| 欧美一区二区精品小视频在线| 国产高清激情床上av| 黄片小视频在线播放| 欧美黑人巨大hd| 老司机福利观看| 99热这里只有是精品50| 欧美+日韩+精品| 一区二区三区激情视频| 丰满的人妻完整版| 久久久精品欧美日韩精品| 我要搜黄色片| 1024手机看黄色片| 欧美成人免费av一区二区三区| 伊人久久精品亚洲午夜| 国产高清视频在线播放一区| 欧美精品国产亚洲| 色播亚洲综合网| 欧美乱妇无乱码| 狂野欧美白嫩少妇大欣赏| 天堂影院成人在线观看| 亚洲国产欧洲综合997久久,| 搡女人真爽免费视频火全软件 | 舔av片在线| 国产真实乱freesex| 高清在线国产一区| 国产精品亚洲美女久久久| 国内少妇人妻偷人精品xxx网站| 色在线成人网| 婷婷色综合大香蕉| 欧美性猛交黑人性爽| 国产真实乱freesex| 亚洲在线观看片| 婷婷精品国产亚洲av| 久久精品91蜜桃| 成人三级黄色视频| 国产免费av片在线观看野外av| 51午夜福利影视在线观看| 男人的好看免费观看在线视频| 久久久久精品国产欧美久久久| 免费av观看视频| 在线观看av片永久免费下载| www.色视频.com| av视频在线观看入口| 午夜福利成人在线免费观看| 搞女人的毛片| 人人妻,人人澡人人爽秒播| 亚洲精品日韩av片在线观看| 国产精品一区二区三区四区久久| 最近中文字幕高清免费大全6 | 18美女黄网站色大片免费观看| 中出人妻视频一区二区| 91九色精品人成在线观看| 神马国产精品三级电影在线观看| 国产三级黄色录像| 欧美色视频一区免费| 老女人水多毛片| 国产极品精品免费视频能看的| 此物有八面人人有两片| 亚洲国产欧洲综合997久久,| 18禁裸乳无遮挡免费网站照片| 国产综合懂色| 老熟妇乱子伦视频在线观看| 在线观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 亚洲一区高清亚洲精品| 蜜桃亚洲精品一区二区三区| 久久国产精品人妻蜜桃| 中文资源天堂在线| .国产精品久久| 国产主播在线观看一区二区| 欧美性猛交╳xxx乱大交人| 中国美女看黄片| 色视频www国产| a级一级毛片免费在线观看| 在线a可以看的网站| 国产三级在线视频| 日本免费a在线| 很黄的视频免费| 老司机深夜福利视频在线观看| 美女大奶头视频| 夜夜爽天天搞| 亚洲国产精品久久男人天堂| 一级黄色大片毛片| 五月玫瑰六月丁香| 狂野欧美白嫩少妇大欣赏| 赤兔流量卡办理| 老熟妇仑乱视频hdxx| 亚洲国产精品sss在线观看| 亚洲美女视频黄频| 亚洲精华国产精华精| 国产精品永久免费网站| 变态另类成人亚洲欧美熟女| 国产不卡一卡二| 69人妻影院| 免费黄网站久久成人精品 | 人人妻人人澡欧美一区二区| 亚洲熟妇中文字幕五十中出| 午夜福利在线观看吧| 乱码一卡2卡4卡精品| 日本黄色视频三级网站网址| 男女下面进入的视频免费午夜| 免费大片18禁| 很黄的视频免费| 国产精品亚洲美女久久久| 少妇丰满av| 日韩高清综合在线| 日韩欧美国产在线观看| 免费在线观看影片大全网站| 久久精品综合一区二区三区| 麻豆av噜噜一区二区三区| 看片在线看免费视频| 全区人妻精品视频| 长腿黑丝高跟| 禁无遮挡网站| 国产激情偷乱视频一区二区| 中出人妻视频一区二区| 亚洲av中文字字幕乱码综合| 嫩草影视91久久| 在线观看美女被高潮喷水网站 | 天堂影院成人在线观看| 免费在线观看影片大全网站| 亚洲最大成人av| 中出人妻视频一区二区| 久久久久免费精品人妻一区二区| 男人舔奶头视频| 日本精品一区二区三区蜜桃| 一级黄片播放器| 国产精品亚洲美女久久久| 亚洲成人中文字幕在线播放| 久久6这里有精品| 高清在线国产一区| 亚洲成人免费电影在线观看| 我的女老师完整版在线观看| 欧美在线一区亚洲| 无人区码免费观看不卡| 国产成+人综合+亚洲专区| www日本黄色视频网| 国产精华一区二区三区| 亚洲av中文字字幕乱码综合| 人人妻人人澡欧美一区二区| 免费无遮挡裸体视频| 他把我摸到了高潮在线观看| 亚洲成av人片免费观看| 欧美高清成人免费视频www| eeuss影院久久| 亚洲欧美清纯卡通| 天堂av国产一区二区熟女人妻| 国产精品影院久久| 国产av不卡久久| 97热精品久久久久久| 首页视频小说图片口味搜索| 老熟妇乱子伦视频在线观看| 黄片小视频在线播放| 国产中年淑女户外野战色| 国内精品一区二区在线观看| 欧美黑人欧美精品刺激| 可以在线观看的亚洲视频| 国内精品一区二区在线观看| 国产私拍福利视频在线观看| 一进一出抽搐动态| 欧美绝顶高潮抽搐喷水| 高清日韩中文字幕在线| 成人av在线播放网站| 一个人看视频在线观看www免费| 男人舔女人下体高潮全视频| 亚洲精品粉嫩美女一区| 亚洲成人久久爱视频| 成人国产综合亚洲| 高清日韩中文字幕在线| 国产淫片久久久久久久久 | 人妻制服诱惑在线中文字幕| 3wmmmm亚洲av在线观看| 午夜福利在线观看吧| 99精品在免费线老司机午夜| 午夜福利18| 午夜福利在线观看吧| 精品久久久久久久久久免费视频| 嫩草影视91久久| 久久久久久久午夜电影| 在线看三级毛片| 99久久九九国产精品国产免费| 男女床上黄色一级片免费看| 很黄的视频免费| 国产精品1区2区在线观看.| 国产精品,欧美在线| 欧美色欧美亚洲另类二区| 亚洲av二区三区四区| 国产欧美日韩精品亚洲av| 国产成人欧美在线观看| 九九在线视频观看精品| www.999成人在线观看| 日韩有码中文字幕| 搡老熟女国产l中国老女人| 午夜免费成人在线视频| 嫁个100分男人电影在线观看| 国产精品不卡视频一区二区 | 99久久成人亚洲精品观看| 一区二区三区高清视频在线| 免费大片18禁| 国产色爽女视频免费观看| 一区二区三区高清视频在线| 久久精品综合一区二区三区| 欧美黄色淫秽网站| 在线观看免费视频日本深夜| 九九热线精品视视频播放| 亚洲第一欧美日韩一区二区三区| 三级国产精品欧美在线观看| 在线观看美女被高潮喷水网站 | 欧美一区二区精品小视频在线| 91久久精品国产一区二区成人| 久久久久久久久中文| 久久精品夜夜夜夜夜久久蜜豆| 最后的刺客免费高清国语| 最近最新中文字幕大全电影3| 两个人的视频大全免费| 超碰av人人做人人爽久久| www.熟女人妻精品国产| av在线蜜桃| 一夜夜www| 精品一区二区三区av网在线观看| 日韩人妻高清精品专区| 欧美zozozo另类| 中文字幕久久专区| 2021天堂中文幕一二区在线观| 97热精品久久久久久| 成人精品一区二区免费| 欧美bdsm另类| 午夜精品在线福利| netflix在线观看网站| 国产免费av片在线观看野外av| 熟女人妻精品中文字幕| 婷婷亚洲欧美| 九九热线精品视视频播放| 内地一区二区视频在线| 精品一区二区三区人妻视频| 亚洲精品粉嫩美女一区| 在线看三级毛片| 国产精品1区2区在线观看.| 亚洲乱码一区二区免费版| 黄片小视频在线播放| 简卡轻食公司| 日日摸夜夜添夜夜添av毛片 | 亚洲成人中文字幕在线播放| 免费观看人在逋| 国产视频内射| 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 亚洲无线在线观看| 日韩欧美精品免费久久 | 国内少妇人妻偷人精品xxx网站| 色噜噜av男人的天堂激情| 亚洲熟妇熟女久久| 精品免费久久久久久久清纯| 女人被狂操c到高潮| 可以在线观看毛片的网站| 婷婷亚洲欧美| 一级a爱片免费观看的视频| 精品乱码久久久久久99久播| 狂野欧美白嫩少妇大欣赏| 国产综合懂色| 人人妻,人人澡人人爽秒播| 亚洲国产精品成人综合色| 亚洲三级黄色毛片| 精品久久久久久,| 免费观看的影片在线观看| 久久久久久久久久成人| 男女床上黄色一级片免费看| 国产精品久久久久久久久免 | 国产精品日韩av在线免费观看| 欧美高清性xxxxhd video| 国产高潮美女av| 毛片一级片免费看久久久久 | 欧美国产日韩亚洲一区| 村上凉子中文字幕在线| 国产成人欧美在线观看| 亚洲片人在线观看| 久久人妻av系列| 美女黄网站色视频| 亚洲第一电影网av| 午夜福利在线观看免费完整高清在 | 97超视频在线观看视频| 一本一本综合久久| 18+在线观看网站| 国内精品美女久久久久久| 亚洲熟妇中文字幕五十中出| 国产精华一区二区三区| 欧美最黄视频在线播放免费| 成人特级黄色片久久久久久久| 我的女老师完整版在线观看| 狂野欧美白嫩少妇大欣赏| 成人美女网站在线观看视频| 午夜两性在线视频| 在线观看一区二区三区| 欧美成人一区二区免费高清观看| 中文字幕免费在线视频6| 两个人的视频大全免费| 精品午夜福利在线看| 精品日产1卡2卡| 欧美一区二区精品小视频在线| 亚洲精品影视一区二区三区av| 亚洲国产日韩欧美精品在线观看| 亚洲精品色激情综合| 蜜桃亚洲精品一区二区三区| 神马国产精品三级电影在线观看| 可以在线观看的亚洲视频| 亚洲国产高清在线一区二区三| 一夜夜www| 亚洲中文日韩欧美视频| 日本精品一区二区三区蜜桃| 欧美高清性xxxxhd video| 色在线成人网| bbb黄色大片| 欧美三级亚洲精品| 色播亚洲综合网| 老司机深夜福利视频在线观看| 综合色av麻豆| 熟妇人妻久久中文字幕3abv| 老鸭窝网址在线观看| 午夜两性在线视频| 18禁黄网站禁片午夜丰满| 日韩欧美国产一区二区入口| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| 久久精品国产99精品国产亚洲性色| 国产乱人视频| av女优亚洲男人天堂| 九九在线视频观看精品| 精品久久久久久久末码| 国产乱人伦免费视频| 午夜亚洲福利在线播放| 亚洲无线观看免费| 我要看日韩黄色一级片| 午夜免费男女啪啪视频观看 | 99久久九九国产精品国产免费| 偷拍熟女少妇极品色| 丁香六月欧美| 床上黄色一级片| 久久国产精品人妻蜜桃| 久久中文看片网| 色尼玛亚洲综合影院| eeuss影院久久| 亚洲一区高清亚洲精品| 国产白丝娇喘喷水9色精品| 免费在线观看影片大全网站| 淫秽高清视频在线观看| 亚洲美女搞黄在线观看 | 一级av片app| 好男人电影高清在线观看| 男人和女人高潮做爰伦理| 成年版毛片免费区| 嫩草影院新地址| 91九色精品人成在线观看| 琪琪午夜伦伦电影理论片6080| 淫妇啪啪啪对白视频| 国产精品,欧美在线| 国产久久久一区二区三区| 丁香欧美五月| 夜夜夜夜夜久久久久| 男人舔奶头视频| bbb黄色大片| 国内少妇人妻偷人精品xxx网站| 亚洲最大成人手机在线| av福利片在线观看| 毛片一级片免费看久久久久 | 欧美一区二区精品小视频在线| 夜夜夜夜夜久久久久| 久久这里只有精品中国| 欧美国产日韩亚洲一区| 国产伦精品一区二区三区视频9| 国产精品久久电影中文字幕| 男女那种视频在线观看| 国产成+人综合+亚洲专区| 18禁在线播放成人免费| 免费av观看视频| a级毛片免费高清观看在线播放| 免费一级毛片在线播放高清视频| 亚洲无线在线观看| netflix在线观看网站| 超碰av人人做人人爽久久| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩卡通动漫| 国产精品1区2区在线观看.| 久久久久九九精品影院| 1000部很黄的大片| 精品一区二区三区视频在线| 日韩精品青青久久久久久| 日本黄色片子视频| 欧美在线黄色| 老司机福利观看| 亚洲自拍偷在线| 夜夜夜夜夜久久久久| 国产午夜精品久久久久久一区二区三区 | 成熟少妇高潮喷水视频| 婷婷六月久久综合丁香| 真实男女啪啪啪动态图| 国产精品一及| 国产又黄又爽又无遮挡在线| 久久中文看片网| 亚洲三级黄色毛片| 色综合站精品国产| 91在线观看av| 欧美性感艳星| 色综合婷婷激情| 亚洲av电影在线进入| 看片在线看免费视频| 欧美精品国产亚洲| 又爽又黄a免费视频| 国产欧美日韩一区二区三| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产亚洲av香蕉五月| 久久久久精品国产欧美久久久| www.色视频.com| 国产 一区 欧美 日韩| 久久中文看片网| 亚洲成人免费电影在线观看| 国产成人影院久久av| 日本黄色视频三级网站网址| 给我免费播放毛片高清在线观看| 国产一级毛片七仙女欲春2| 久久午夜亚洲精品久久| 亚洲中文字幕日韩| 一本久久中文字幕| 亚洲中文字幕日韩| 一进一出抽搐动态| 久久99热这里只有精品18| 中文字幕av在线有码专区| 久久久久久久精品吃奶| 日本免费一区二区三区高清不卡| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久大精品| 亚洲经典国产精华液单 | 村上凉子中文字幕在线| 男女那种视频在线观看| 国产伦人伦偷精品视频| 免费在线观看影片大全网站| 亚洲avbb在线观看| a级毛片免费高清观看在线播放| 国产精品1区2区在线观看.| 美女 人体艺术 gogo| 精品一区二区免费观看| 中文资源天堂在线| 成年人黄色毛片网站| 丁香欧美五月| 可以在线观看毛片的网站| 啪啪无遮挡十八禁网站| 国产乱人视频| 国产在视频线在精品| 美女xxoo啪啪120秒动态图 | 美女高潮的动态| 亚洲18禁久久av| 亚洲国产精品合色在线| 18美女黄网站色大片免费观看| 亚洲第一区二区三区不卡| www.999成人在线观看| 一区二区三区免费毛片| 国产毛片a区久久久久| 亚洲精华国产精华精| 国产精品1区2区在线观看.| 免费观看的影片在线观看| 人妻制服诱惑在线中文字幕| 嫩草影院新地址| 成人高潮视频无遮挡免费网站| 一个人免费在线观看电影| 国产一区二区三区视频了| 亚洲av不卡在线观看| 国产精品日韩av在线免费观看| 麻豆成人av在线观看| 成人鲁丝片一二三区免费| 天堂√8在线中文| 国产毛片a区久久久久| 狂野欧美白嫩少妇大欣赏| 久久精品国产清高在天天线| 亚洲无线观看免费| 1024手机看黄色片| 草草在线视频免费看| 好看av亚洲va欧美ⅴa在| 日本一二三区视频观看| 少妇被粗大猛烈的视频| 午夜日韩欧美国产| 午夜福利视频1000在线观看| 美女免费视频网站| 小蜜桃在线观看免费完整版高清| 日韩欧美精品免费久久 | 欧美色欧美亚洲另类二区| 波多野结衣高清作品| 伦理电影大哥的女人| 免费在线观看成人毛片| 五月玫瑰六月丁香| 内射极品少妇av片p| 久久久久九九精品影院| 97碰自拍视频| 99热这里只有精品一区| 长腿黑丝高跟| 国产免费男女视频| 十八禁人妻一区二区| 国产乱人视频| 成人亚洲精品av一区二区| 亚洲自偷自拍三级| 色5月婷婷丁香| www.www免费av| 男人舔奶头视频|