——以鹽城越冬丹頂鶴及其生境的變化為例"/>
  • <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    耦合種群動態(tài)的生境格局變化分析粒度與景觀因子選擇
    ——以鹽城越冬丹頂鶴及其生境的變化為例

    2016-07-21 01:11:02歐維新甘玉婷婷
    生態(tài)學報 2016年10期

    歐維新,甘玉婷婷

    1 南京農(nóng)業(yè)大學 土地管理學院,南京 210095 2 農(nóng)村土地資源利用與整治國家地方聯(lián)合工程研究中心,南京 210095

    ?

    耦合種群動態(tài)的生境格局變化分析粒度與景觀因子選擇
    ——以鹽城越冬丹頂鶴及其生境的變化為例

    歐維新1,2,*,甘玉婷婷1

    1 南京農(nóng)業(yè)大學 土地管理學院,南京210095 2 農(nóng)村土地資源利用與整治國家地方聯(lián)合工程研究中心,南京210095

    摘要:選擇鹽城珍禽國家級自然保護區(qū),以丹頂鶴越冬生境景觀與種群動態(tài)這一格局-過程關系為切入點,嘗試從丹頂鶴最小存活面積特征與種群動態(tài)視角探討一下最佳粒度、景觀格局變化分析的指數(shù)遴選方法。結(jié)果發(fā)現(xiàn):根據(jù)丹頂鶴最小生存面積確定200m為最大轉(zhuǎn)換粒度,綜合景觀指數(shù)隨空間粒度變化的規(guī)律和粒度轉(zhuǎn)換精度損失評價的結(jié)果,確定最佳分析粒度為70m。在眾多景觀因子中,運用Spearman秩相關分析,再利用逐步回歸分析建立起景觀與丹頂鶴數(shù)量指標間的聯(lián)系,最終篩選出反映景觀面積(CA)和空間格局(IJI,ENN_MN)的3種影響顯著景觀因子,其解釋貢獻率(累計R2)達到70.5%,且其所代表的生境景觀的組成和結(jié)構(gòu)信息較為客觀地反映出研究區(qū)丹頂鶴種群動態(tài)分布的顯示狀態(tài)。與純粹的景觀指數(shù)分析相比,這種方法更能反應景觀格局演變特征的特定生態(tài)學意義。

    關鍵詞:粒度效應;景觀指數(shù);生態(tài)過程;景觀格局變化;越冬丹頂鶴生境

    基于景觀生態(tài)學理論的景觀格局動態(tài)變化分析是當前研究珍稀物種生境保護的主要方法之一[1- 4]。景觀指數(shù)以其具有反映景觀各要素的結(jié)構(gòu)和空間配置信息而得以廣泛應用[5]。在用于真實刻畫景觀格局變化及其生態(tài)過程,景觀指數(shù)的通常需考慮3個方面,即選擇合適的研究分析尺度,如一定空間幅度下如何選擇與特定生態(tài)過程關聯(lián)的景觀粒度水平;指數(shù)間需滿足統(tǒng)計上的相互獨立性;所選指數(shù)對所研究的生態(tài)過程具有很好的敏感性或指示性。

    粒度是指最小可辨識單元所代表的特征長度、面積或體積(如樣方、像元)[6]。進行景觀格局與生態(tài)過程分析時,如果粒度過小,會造成區(qū)域空間信息數(shù)據(jù)量過大而掩蓋一些重要信息;如果粒度過大,則又造成關鍵細節(jié)信息的缺失[7],已有的研究表明,隨著粒度的粗化,不同的景觀類型會有不同的響應結(jié)果,景觀格局的最佳分析粒度指的是針對某一景觀類型分析的最合適尺度[8]。故最佳粒度因能真實地反映生境類型景觀的生態(tài)特征而顯得尤為重要。由于景觀指數(shù)具有明顯的尺度依賴性[9],粒度效應研究具有確定景觀指數(shù)隨空間粒度變化臨界閾值、辨明景觀指數(shù)的反應敏感性等作用[10- 11],因而粒度效應分析可以作為遴選景觀最佳分析粒度與敏感性指數(shù)的可靠方法。由于景觀指數(shù)眾多且多數(shù)極具高度相關性,隨意地指數(shù)選擇易造成指數(shù)間的冗余。因此,景觀格局指數(shù)選擇時的相關性研究已成為必須的環(huán)節(jié)。目前對該方面的研究已有不少,但更多地僅限于研究景觀指數(shù)本身的相關性考察[12- 13],結(jié)合特定生態(tài)過程進行指數(shù)的遴選研究相對較少。

    本文選擇鹽城珍禽國家級自然保護區(qū),以丹頂鶴越冬生境及其種群動態(tài)這一格局-過程關系為切入點,嘗試探討一下最佳粒度、景觀格局變化分析的指數(shù)遴選方法。

    1研究區(qū)域

    1.1研究區(qū)域概況

    本文的研究區(qū)為鹽城珍禽國家級自然保護區(qū)的核心區(qū)與緩沖區(qū),由3個越冬丹頂鶴熱點小區(qū)組成(圖1),該區(qū)地處江蘇中部沿海,太平洋西海岸,北至黃沙港,南至四卯酉河,西接海堤,為淤漲型潮灘濕地,總面積有87×104hm2。其中核心區(qū)為保存最為完好的原生灘涂濕地,其南北兩側(cè)的緩沖區(qū)有較多的人類開發(fā)活動,保護與開發(fā)的矛盾最為突出。研究區(qū)擁有蘆葦濕地、堿蓬濕地、米草濕地、淤泥光灘、養(yǎng)殖水域等濕地類型,每年有數(shù)百萬計水禽在此棲息或越冬,來鹽城越冬的丹頂鶴數(shù)量占世界種群的一半以上[14],其中,80%左右的丹頂鶴集中在研究區(qū)越冬棲息。研究區(qū)因保留了大片原生生境使其成為丹頂鶴越冬棲息的關鍵地段,而人類圍墾開發(fā)活動使緩沖區(qū)生境景觀變化劇烈進而對丹頂鶴棲息條件造成很大影響。因此,該區(qū)景觀尺度變化、景觀格局變化等特征在響應丹頂鶴的種群動態(tài)具有典型性。

    圖1 研究區(qū)位置示意圖Fig.1 Location and scope of the study area

    1.2數(shù)據(jù)來源與處理

    利用1991—1992年、1995年、1997年、1999—2006年、2010年軌道號為119/037、120/036共計13期TM影像作為基礎數(shù)據(jù),運用ENVI4.8軟件進行幾何校正、影像拼接、影像裁剪等預處理后,得到研究區(qū)遙感圖。采用分區(qū)分類的方法依據(jù)海岸帶將研究區(qū)劃分為陸地和海灘部分,對于陸地部分,本研究采用最大似然法,對影像進行監(jiān)督分類,提取研究區(qū)土地利用類型特征,在ARCGIS中將研究區(qū)劃分為9類,分別是耕地、建筑用地、米草濕地、河流水域、淤泥光灘、蘆葦白茅濕地、堿蓬濕地、水產(chǎn)養(yǎng)殖塘、鹽田。根據(jù)丹頂鶴在鹽城越冬的習性特征(即它們通常喜歡在蘆葦和堿蓬地棲息,偶爾會到鹽田和水塘去覓食,但會回避淤泥光灘、互花米草、道路、河流和居民點等一些人類活動頻繁的地方),結(jié)合丹頂鶴越冬生境適宜性特征及已有成果[15- 16],本文將研究區(qū)地類進一步歸并為越冬丹頂鶴的適宜生境(包括蘆葦白茅濕地、堿蓬濕地)、補充生境(包括水產(chǎn)養(yǎng)殖塘、鹽田)、非生境(其他利用與覆被類型)三大景觀類型[17]。

    1991—2006年間研究區(qū)所含的3個丹頂鶴越冬棲息熱點地區(qū)(核心區(qū)及射陽蘆葦基地、射陽鹽場及灘涂、四卯酉王港灘涂)中的丹頂鶴數(shù)據(jù),來源于江蘇鹽城國家級珍禽自然保護區(qū)管理局、中國林科院鳥類環(huán)志中心共同開展的丹頂鶴越冬數(shù)量監(jiān)測區(qū)調(diào)查結(jié)果[18]。本文只選擇了與可獲的遙感影像數(shù)據(jù)相應年份的12期丹頂鶴數(shù)據(jù)。

    2研究方法

    本研究旨在結(jié)合丹頂鶴最小存活面積及種群動態(tài)等生態(tài)特征來考察丹頂鶴越冬生境分析的最佳粒度及適宜景觀指數(shù)。為此,研究擬采用粒度效應分析與粒度轉(zhuǎn)化的精度損失分析來確定最佳粒度;適宜景觀指數(shù)則結(jié)合景觀因子的獨立性分析方法與基于種群動態(tài)與格局因子的逐步回歸分析方法來遴選。

    2.1最佳粒度的確定方法

    根據(jù)文獻可知,7.8 km2的生境區(qū)域是丹頂鶴生存的最小需求面積[19- 20],即丹頂鶴單個生存的生境斑塊的直徑必須達到200m。為獲得最佳分析粒度,利用重采樣的方法,將2006年景觀類型圖進行柵格像元的聚合,柵格單元的賦值采用最大面積值法(RMA)。以30—200m為粒度轉(zhuǎn)換的考察區(qū)間,每10m為粒度增幅,劃分并生成18幅不同粒度水平的柵格圖。

    由于area/density/edge等方面的指數(shù)隨粒度的變化比較敏感,且有明顯拐點,其粒度效應方程通常呈現(xiàn)簡單線性或?qū)?shù)函數(shù)[21],因而選取平均斑塊面積(AREA_MN)、斑塊數(shù)量(NP)、邊界總長度(TE)、邊界密度(ED)、周長面積比分維數(shù)(PAFPAC)、面積加權平均分維數(shù)6個代表指數(shù),計算并比較各指數(shù)的粒度拐點。同時,考慮到矢量數(shù)據(jù)的柵格化過程中存在著面積等方面的精度損失[22],本文以景觀斑塊面積與斑塊數(shù)量為對象,分析比較不同粒度水平下的整體精度損失。

    不同粒度水平下的精度損失計算公式如下[23]:

    (1)

    (2)

    式中,E為某一評價指標損失量,L表示損失精度,Agi為i類景觀類型柵格數(shù)量或面積,Abi為i類景觀類型矢量斑塊數(shù)量和面積,n為景觀類型個數(shù)。

    2.2景觀因子篩選方法

    從景觀面積-周長-密度、形狀、隔離與破碎化程度方面選擇19種常用景觀指數(shù),為丹頂鶴生境景觀格局演變分析提供關鍵性的景觀指標,運用相關分析進行因子篩選:

    (1)用spearman等級相關分析對所有指數(shù)進行兩兩之間的相關分析,檢測指數(shù)間的獨立性,選擇獨立性較強的指數(shù);

    (2)為遴選對丹頂鶴分布有重要影響的關鍵生境景觀因子,采用逐步回歸方法建立12年3個熱點地區(qū)的丹頂鶴分布數(shù)據(jù)[15]與景觀因子的聯(lián)系,剔除獨立性差且和數(shù)量變化相關性小的指數(shù),獲得反映丹頂鶴生境格局變化的主要指數(shù)。需要說明一下的是,為了消除世界種群數(shù)量動態(tài)對鹽城丹頂鶴種群數(shù)量的影響,本文將所涉3個熱點小區(qū)每年來越冬的丹頂鶴絕對數(shù)量轉(zhuǎn)化為相對數(shù)量,即用當年每個小區(qū)丹頂鶴數(shù)量與3個小區(qū)丹頂鶴總量的比值來表示。以上分析均在SPSS軟件中進行。

    3結(jié)果與分析

    3.1丹頂鶴越冬生境景觀格局最佳粒度的確定

    3.1.1景觀指數(shù)粒度效應的分析

    根據(jù)粒度效應圖(圖2)可知,隨著粒度變大,6個景觀指數(shù)的響應曲線大致分為兩種情況。

    第一種情況景觀指數(shù)隨粒度變化時拐點明顯,在粒度從40m到100m的變化過程中,平均斑塊面積出現(xiàn)了3次明顯的波動變化,分別是40—60m時先增大后減?。?0—80m時先急劇上升再緩慢下降;80—100m時先上升再急劇下降,后兩次的變化幅度較大。而超過100m后,平均斑塊面積逐漸平緩下降,變化幅度越來越小,最終趨于穩(wěn)定,不再出現(xiàn)突變現(xiàn)象。這是由于粒度達到100m后,景觀融合程度加強,各斑塊面積也越來越均勻。在這過程中,變化明顯的拐點為40、60、70、90m,對應的尺度域為(40,60)、(60,70)、(70,90);與平均斑塊面積的變化相比,斑塊數(shù)量呈現(xiàn)相反的變化,而其波動幅度更小一些。在整個過程中,出現(xiàn)了兩次較明顯的上下波動,分別是在粒度為50—70m和80—100m時,其轉(zhuǎn)折點為60、70、90m,之后出現(xiàn)波動幅度較小,其轉(zhuǎn)折點為140、180m。則對應的尺度域為(60,70)、(70,90)、(90,140)、(140,180);周長面積比分維數(shù)和面積加權分維數(shù)隨粒度的變化整體都是呈現(xiàn)速率逐漸增大的上升而后速率較小的下降趨勢,前者幾個明顯的轉(zhuǎn)折點出現(xiàn)在粒度為90、100、130、140、180m,對應的尺度域為(90,100)、(100,130)、(130,140)、(140,180)。面積加權分維數(shù)的轉(zhuǎn)折點出現(xiàn)在60、70、120、130、170、180m,對應的尺度域為(60,70)、(70,120)、(120,130)、(130,170)、(170,180)。

    表1 用于相關分析的初選景觀指數(shù)表

    圖2 景觀指數(shù)粒度效應圖Fig.2 Scale effect curves of landscape metrics

    第二種情況隨粒度變化有規(guī)律,但拐點并不明顯,主要是隨著粒度的增大,總體都呈現(xiàn)遞減的趨勢的邊界密度和邊界總長度。遞減的速率小于斑塊數(shù)和平均斑塊面積遞減的速率。這說明粒度增大后,柵格單元的擴大使得一些優(yōu)勢地類受到道路、河流等的切割作用減小,形狀更為規(guī)則,同時斑塊邊界的復雜程度的改變比面積數(shù)量的改變要慢。在遞減過程中,邊界密度出現(xiàn)了波動幅度較小的情況,是在粒度為70—100m時,其拐點為70、90m,其尺度域為(70,90)。邊界總長度無數(shù)值上下波動的現(xiàn)象,則沒有粒度轉(zhuǎn)折點。

    根據(jù)這些指數(shù)的粒度分異規(guī)律,綜合分析可知變化明顯的轉(zhuǎn)折點多集中在60、70、90、140、180m,所以對于整個研究區(qū)景觀,尺度域是(60,70)、(70,90)、(90,140)、(140,180)。相關研究表明,景觀指數(shù)隨粒度變化的第一尺度域(前兩個拐點)是景觀格局中最具有比例尺的信息特征,也是粒度選擇的適宜取值范圍[23],因此本研究區(qū)內(nèi)適宜取值范圍為(60,70)。

    3.1.2粒度轉(zhuǎn)化精度損失分析

    利用粒度大小作為橫坐標,分別以損失面積、數(shù)量及損失精度為縱坐標得到精度損失曲線圖。景觀整體的面積損失在30—50m,70—110m變化平緩,而在50—70、110—200m起伏較大;斑塊 數(shù)量損失在30—100m和150 —200m變化平緩,在110 —150m起伏較大,確定較小精度損失的粒度范圍為(30,70)。

    圖3 不同粒度下各指標信息損失分析結(jié)果Fig.3 The result of accuracy loss evaluating under different grain sizes

    綜上兩種方法得到的粒度域,取交集得到(60,70)為最佳分析粒度范圍,考慮到選擇較大的粒度可以適當避免冗余的工作量,因此本文最終確定70m為鹽城丹頂鶴自然保護區(qū)的最佳分析粒度。

    3.2景觀指數(shù)選擇

    3.2.1基于面積相關的指數(shù)剔除

    在進行景觀格局分析中,斑塊面積已被公認為是最為主要的測度因子,因此首先進行景觀指數(shù)與適宜生境面積的Spearman等級相關分析可以盡可能多地剔除冗余的指數(shù)。本文以適宜生境為對象,該類生境的相關景觀指數(shù)相關分析結(jié)果表明(表2),LPI、ED、LSI、AREA_MN、SHAPE_MN、PLADJ、DIVISION、COHESION、SPLIT、AI等10個與面積變化相關性明顯顯著,故首先剔除這些指數(shù)。

    3.2.2基于獨立性檢驗的指數(shù)剔除

    對第一步篩選出來的8個指標,計算兩種生境類型8對景觀指數(shù)兩兩間的Spearman相關系數(shù),進一步從獨立性檢驗上來遴選,并得到9個相對獨立的指標(表3),即PARA_MN_M、ENN_MN_M、IJI_M、CONNECT_M、PD_S、GYRATE_MN_S、ENN_MN_S、IJI_S、CONNECT_S。

    表2 景觀指數(shù)與面積相關性情況表

    *P<0.05;**P<0.01

    表3 9個景觀指數(shù)的Spearman秩相關系數(shù)矩陣

    U1、U2、U3、U4、U5、U6、U7、U8、U9分別表示PARA_MN_M、ENN_MN_M、IJI_M、CONNECT_M、PD_S、GYRATE_MN_S、ENN_MN_S、IJI_S、CONNECT_S; 其中“_M”表示適宜生境;“_S”表示補充生境

    3.2.3基于丹頂鶴種群動態(tài)與格局變化關系的指數(shù)遴選

    以12期3個熱點小區(qū)越冬丹頂鶴的相對數(shù)量為因變量,以上述分析獲得的9個指數(shù)和生境面積CA(包括適宜生境面積CA_M與補充生境面積CA_S)為自變量,進行逐步回歸分析,得到影響丹頂鶴數(shù)量動態(tài)變化的4個顯著因子(表4)。

    從結(jié)果上看,兩類生境景觀中反映景觀面積(CA)和空間格局(IJI,ENN_MN)的3種景觀指數(shù)在解釋丹頂鶴種群變化上具有顯著影響效果,解釋貢獻率(累計R2)達到70.5%。即這4個景觀指數(shù)所代表的景觀格局變化信息與越冬丹頂鶴種群動態(tài)分布特征顯著相關,且越冬丹頂鶴種群動態(tài)及分布結(jié)果有70.5%的信息可由這4個景觀指數(shù)來解釋。因此在分析丹頂鶴種群動態(tài)分析方面,這幾個指數(shù)將是很好的預測變量。

    表4 逐步回歸分析結(jié)果

    生境面積是影響種群變化的最根本的原因,在丹頂鶴越冬適宜生境大面積減少的同時,被圍墾開發(fā)利用的大量補充生境(養(yǎng)殖魚塘、鹽田)成為影響丹頂鶴數(shù)量種群數(shù)量變化的關鍵因素[1- 2]。表4中引入的4個因子中有3個都是與補充生境有關的。其中,CA_S與種群數(shù)量呈顯著負相關,說明大規(guī)模的補充生境面積的增加盡管能為丹頂鶴補充一定的食物和水源,但因其提供的隱蔽條件不夠好,再加上養(yǎng)殖、曬鹽過程中所帶來的人為干擾會對丹頂鶴的集聚造成不利影響。IJI_S與種群數(shù)量顯著正相關、ENN_MN_S與種群數(shù)量顯著正相關,前者說明被其他多種生境類型斑塊(尤其是適宜生境)所毗連、包圍的補充生境斑塊能為丹頂鶴提供安全棲息場所,后者則指出離散分布的補充生境格局有利于丹頂鶴種群的集聚,這兩個指數(shù)進一步從生境景觀空間格局上反映了集中連片的補充生境不利于丹頂鶴棲息和集聚。IJI_M與種群數(shù)量顯著正相關,說明現(xiàn)有適宜生境因為面積損失與破碎化的不斷加劇過程中,其鄰近補充生境能為丹頂鶴提供食物、水源等資源條件而有利于物種的集聚。因此,從引入的景觀因子的生態(tài)學意義來看,其所代表的生境景觀的組成和結(jié)構(gòu)特征也很好地解釋了丹頂鶴種群動態(tài)這一生態(tài)過程。

    4結(jié)論與討論

    (1)本文根據(jù)丹頂鶴最小存活面積確定200m為最大轉(zhuǎn)換粒度,綜合景觀指數(shù)隨空間粒度變化的規(guī)律和粒度轉(zhuǎn)換精度損失評價的結(jié)果,得到研究區(qū)域適宜粒度域為60—70m,并確定最佳分析粒度為70m。尺度選擇是景觀格局研究分析的基礎[22- 23],在適宜粒度下進行景觀格局分析能夠比較準確地表達研究區(qū)的景觀特征概況。目前,通過景觀指數(shù)粒度效應曲線分析確定最佳景觀格局分析粒度的研究較為常見[24- 27],但考慮結(jié)合粒度轉(zhuǎn)換中發(fā)生的某些信息損失評價方法確定最佳粒度并不多見。結(jié)合已有關于粒度效應的一些研究,本研究在粒度效應分析上,主要選用了隨粒度變化具有明顯規(guī)律、易解釋、算式簡便的面積-周長-密度等景觀指數(shù),可以較快地識別出粒度域,簡化了分析粒度效應的程序。另外,粒度轉(zhuǎn)換的精度損失研究中,但大多都只考慮轉(zhuǎn)換過程中的景觀面積損失,本文還增加了斑塊數(shù)量評價。同時,本研究以丹頂鶴最小存活面積這一生態(tài)特征來確定最大轉(zhuǎn)換粒度,這比當前大部分粒度效應研究中單純考慮粒度大小的物理變化更具有生態(tài)學意義,也減少了粒度轉(zhuǎn)換尺度選擇的隨意性。

    (2)丹頂鶴種群動態(tài)分布是生境景觀格局的變化的直接響應,因此遴選并明確影響丹頂鶴數(shù)量分布的關鍵景觀因子,有助于更好地理解景觀格局與生態(tài)過程的關系。本文在眾多景觀因子中,運用Spearman秩相關分析,從整體上揭示景觀指數(shù)的信息重疊;再利用逐步回歸分析建立起景觀與丹頂鶴數(shù)量指標間的聯(lián)系,最終篩選出關鍵的景觀因子。與純粹的景觀指數(shù)分析相比,這種方法更能反應景觀格局的演變特征的特定生態(tài)學意義。

    (3)本文最終篩選出關鍵的4個影響顯著景觀因子,分別是兩類生境景觀中反應景觀面積(CA)和空間格局(IJI,ENN_MN)的3種景觀指數(shù),解釋貢獻率(累計R2)達到70.5%,且其所代表的生境景觀的組成和結(jié)構(gòu)信息較為客觀地反映了研究區(qū)丹頂鶴種群動態(tài)分布的顯示狀態(tài)。同時研究也表明:在丹頂鶴越冬適宜生境大面積減少和破碎化背景下,大量被圍墾開發(fā)利用的補充生境(養(yǎng)殖魚塘、鹽田)已經(jīng)成為影響丹頂鶴種群數(shù)量變化的關鍵因素;盡管補充生境能為丹頂鶴越冬棲息提供一定的食物和水源條件,但集中連片的補充生境不利于丹頂鶴棲息和集聚。

    (4)從本研究的實證來看,因為數(shù)據(jù)可獲得性的限制,還有一些有待進一步深化完善的地方,比如,研究中采用的最大轉(zhuǎn)換粒度是根據(jù)丹頂鶴繁殖地的最小存活面積,但研究區(qū)是丹頂鶴用于覓食和夜棲的越冬地,它們在鹽城的繁殖領域性不明顯,但由于當前仍缺乏這方面的研究結(jié)果,因而采用這個數(shù)據(jù)來替代。但應該指出的是,本研究嘗試如何把生態(tài)過程或生態(tài)特征與粒度、格局等更好地結(jié)合起來,更多地是從方法上來探討景觀格局分析中粒度和景觀指數(shù)選擇問題,這對今后的景觀格局分析研究還是具有一定的理論參考價值。

    參考文獻(References):

    [1]張曼胤. 江蘇鹽城濱海濕地景觀變化及其對丹頂鶴生境的影響. 長春:東北師范大學,2008.

    [2]孫賢斌, 劉紅玉. 江蘇鹽城海濱區(qū)域丹頂鶴適宜越冬生境變化. 生態(tài)學雜志, 2011, 30(4): 694- 699.

    [3]Cao M C, Liu G H. Habitat suitability change of red-crowned crane in Yellow River Delta Nature Reserve. Journal of Forestry Research, 2008, 19(2): 141- 147.

    [4]Horne B V. Density as a misleading indicator of habitat quality. The Journal of Wildlife Management, 1983, 47(4): 893- 901.

    [5]張金屯, 邱揚, 鄭鳳英. 景觀格局的數(shù)量研究方法. 山地學報, 2000, 18(4): 346- 352.

    [6]Jenerette G D, Wu J G. On the definitions of scale. Bulletin of the Ecological Society of America, 2000, 81(4): 104- 105.

    [7]Qi Y, Wu J G. Effects of changing spatial resolution on the results of landscape pattern analysis using spatial autocorrelation indices. Landscape Ecology, 1996, 11(1): 39- 49.

    [8]Wu J G. Effects of changing scale on landscape pattern analysis: scaling relations. Landscape Ecology, 2004, 19(2): 125- 138.

    [9]Wu J G, Jelinski D E, Luck M, Tueller P T. Multiscale analysis of landscape heterogeneity: scale variance and pattern metrics. Geographic Information Science, 2000, 6(1): 6- 19.

    [10]張東水, 蘭樟仁, 丘榮祖. 閩江口濕地遙感影像最佳景觀觀察尺度的選擇. 遙感信息, 2006, (4): 30- 32.

    [11]張慶印, 樊軍. 高精度遙感影像下農(nóng)牧交錯帶小流域景觀特征的粒度效應. 生態(tài)學報, 2013, 33(24): 7739- 7747.

    [12]甄 霖, 謝高地, 楊 麗, 成升魁. 涇河流域分縣景觀格局特征及相關性. 生態(tài)學報, 2005, 25(12): 3343- 3353.

    [13]布仁倉, 胡遠滿, 常禹, 李秀珍, 賀紅士. 景觀指數(shù)之間的相關分析. 生態(tài)學報, 2005, 25(10): 2764- 2775.

    [14]馬志軍, 李文軍, 王子健, 湯鴻霄. 中國鹽城生物圈保護區(qū)丹頂鶴(Grusjaponensis)生境的變化及其適應性. AMBIO-人類環(huán)境雜志, 1998, (6): 461- 464.

    [15]胡遠滿, 舒瑩, 李秀珍, 王凌, 李玉祥, 楊玉成. 遼寧雙臺河口自然保護區(qū)丹頂鶴繁殖生境變化及其繁殖容量分析. 生態(tài)學雜志, 2004, 23(5): 7- 12.

    [16]呂士成. 鹽城沿海灘涂丹頂鶴的分布現(xiàn)狀及其趨勢分析. 生態(tài)科學, 2008, 27(3): 154- 158.

    [17]歐維新, 逄謙, 甘玉婷婷. 鹽城海濱濕地資源利用變化及其對丹頂鶴越冬生境的影響. 中國人口·資源與環(huán)境, 2014, 24(7): 30- 36.

    [18]呂士成. 丹頂鶴越冬期在人工濕地的夜棲數(shù)量分布. 野生動物雜志, 2007, 28(2): 11- 13.

    [19]Su L Y. Comparative Feeding Ecology of the Red-Crowned and White-Naped Cranes. Columbia, Missouri: University of Missouri, 1993.

    [20]Csuti B, Kiester A R. Hierchical Gap Analysis for Identifying Priority Areas for Biodiversity // Gap Analysis—A Landscape Approach to Biodiversity Planning. Bethesda, Maryland: American of Society for Photogrammetry and Sensing, 1996: 23- 58.

    [22]趙文武, 傅伯杰, 陳利項. 尺度推繹研究中的幾點基本問題. 地球科學進展, 2002, 17(6): 905- 911.

    [23]楊存建, 劉紀遠, 張增祥, 王長有. 土地利用數(shù)據(jù)尺度轉(zhuǎn)換的精度損失分析. 山地學報, 2001, 19(3): 258- 264.

    [24]趙文武, 傅伯杰, 陳利頂. 景觀指數(shù)的粒度變化效應. 第四紀研究, 2003, 23(3): 326- 333.

    [25]Li H, Wu J G. Use and misuse of landscape indices. Landscape Ecology, 2004, 19(4): 389- 399.

    [26]Griffith J A, Martinko E A, Price K P. Landscape structure analysis of Kansas at three scales. Landscape and Urban Planning, 2000, 52(1): 45- 61.

    [27]于磊, 趙彥偉, 張遠, 李敏. 基于最佳分析粒度的大遼河流域濕地景觀格局分析. 環(huán)境科學學報, 2011, 31(4): 873- 879.

    Grain size and landscape indices selection by coupling population dynamics and habitat pattern analysis: a case study of wintering red-crowned crane and its habitat in Yancheng

    OU Weixin1,2,*, GAN Yutingting1

    1CollegeofLandManagement,NanjingAgriculturalUniversity,Nanjing210095,China2National&LocalJointEngineering,ResearchCenterforRuralLandResourcesUseandConsolidation,Nanjing210095,China

    Abstract:Landscape pattern is scale-dependent. Thus, understanding landscape structure and functioning requires multiscale information. Scaling functions are the most precise and concise methods for explicitly quantifying multiscale characteristics. If landscape indices are ecologically relevant and reflect important attributes of spatial pattern, they can functionally link the dynamics of ecological processes to landscape structure. Therefore, the selection of appropriate scale (e.g., grain size) and landscape indices are critical to landscape pattern analysis. The major objective of this study was to explore how optimal grain size and landscape indices can be selected for landscape pattern analysis, to improve our understanding and prediction of ecological processes. For this, we conducted a case study on wintering Red-crowned crane and its habitat in Yancheng. To obtain the optimal grain size for landscape pattern analysis of the wintering habitat of the Red-crowned crane, we followed two steps. The first step was scale (grain) effect analysis, to find a grain size that the test metrics could detect in case of any variation (sensitivity). In this step, the minimum survival area for the species was considered for identifying the maximum grain size for scaling. The second step was accuracy test by evaluating the loss of landscape area and patch numbers in each grain size level. The optimal grain size was obtained by integrated analysis of the results from the above two steps. In landscape indices selection, 19 landscape metrics (computed at the optimal grain size) were subjected to Spearman rank correlation analysis to assess the independence. Then, step-wise regressions were performed to evaluate the effects of the spatial attributes (landscape indices with higher independence) at three hotspots on the abundance of cranes (based on population dynamics of the past 12 years[1]). Variables with significance level above 90% were selected as the optimal landscape indices for pattern analysis. The results showed the following: (1) on the basis of the minimum survival area of the red-crowned crane, 200 m was the highest obtained grain size, and 70 m was the optimal grain size identified by integrating the results of scale (grain) effect and accuracy after assessing the grain conversion area loss. (2) Effects of landscape patterns at the hotspots were analyzed using 8 pairs of landscape indices (results of metrics for the most suitable and supplement habitat types for cranes were computed) as the independent variables. Three landscape indices (CA, IJI, ENN_MN) of two habitat types were found to be significant (R2 = 70.5%). (3) The selected optimal landscape indices of the supplement habitat showed positive effect for isolation and negative effect for area size on population abundance. Our results imply that supplement habitat may provide complementary resource sites for cranes, but continued species concentration at these habitats may negatively affect crane abundance and distribution due to induced human disturbance. Our combined findings on optimal grain size and landscape indices proved to be satisfactory for landscape pattern analysis, suggesting that our approaches are reasonable. Further, compared to landscape pattern analysis by using random grain size and simple landscape metrics, the results obtained using our approaches, as shown in this case study, are ecological relevant.

    Key Words:grain effect; metrics; ecological process; landscape pattern change; wintering red-crowned crane habitat

    基金項目:江蘇省青藍工程項目

    收稿日期:2014- 12- 10; 網(wǎng)絡出版日期:2015- 10- 10

    *通訊作者

    Corresponding author.E-mail: wxounj@163.com

    DOI:10.5846/stxb201412102450

    歐維新,甘玉婷婷.耦合種群動態(tài)的生境格局變化分析粒度與景觀因子選擇——以鹽城越冬丹頂鶴及其生境的變化為例.生態(tài)學報,2016,36(10):2996- 3004.

    Ou W X, Gan Y T T.Grain size and landscape indices selection by coupling population dynamics and habitat pattern analysis: a case study of wintering red-crowned crane and its habitat in Yancheng.Acta Ecologica Sinica,2016,36(10):2996- 3004.

    亚洲国产最新在线播放| 啦啦啦韩国在线观看视频| 国产精品av视频在线免费观看| 91精品一卡2卡3卡4卡| 蜜桃久久精品国产亚洲av| 最近中文字幕2019免费版| 麻豆av噜噜一区二区三区| 男人爽女人下面视频在线观看| 在线 av 中文字幕| 好男人视频免费观看在线| 99久国产av精品国产电影| 国产人妻一区二区三区在| 日韩一区二区三区影片| 99热6这里只有精品| 国产精品一二三区在线看| 中文字幕亚洲精品专区| 亚洲图色成人| 黄色欧美视频在线观看| 亚洲人成网站在线观看播放| 亚洲精品一二三| 久久鲁丝午夜福利片| 在线a可以看的网站| 日韩欧美 国产精品| 国产精品福利在线免费观看| 免费看a级黄色片| 2021天堂中文幕一二区在线观| 中文天堂在线官网| 午夜日本视频在线| 国产精品精品国产色婷婷| 久久精品夜色国产| 国产伦在线观看视频一区| 亚洲精品日韩在线中文字幕| 超碰av人人做人人爽久久| 久久久久性生活片| 99re6热这里在线精品视频| 免费高清在线观看视频在线观看| 观看美女的网站| 国产亚洲一区二区精品| 亚洲在久久综合| av在线观看视频网站免费| 成人特级av手机在线观看| 18禁裸乳无遮挡免费网站照片| 99热这里只有精品一区| 国产成人精品福利久久| 亚洲国产精品专区欧美| 日本一二三区视频观看| 欧美潮喷喷水| 亚洲精品,欧美精品| 舔av片在线| 午夜福利视频1000在线观看| 亚洲欧美一区二区三区国产| 最近手机中文字幕大全| 成人性生交大片免费视频hd| 成人二区视频| 成人鲁丝片一二三区免费| 亚洲图色成人| 国产日韩欧美在线精品| 国产伦精品一区二区三区四那| 日韩在线高清观看一区二区三区| 97超碰精品成人国产| 久久久久久久久久人人人人人人| 国产永久视频网站| 免费播放大片免费观看视频在线观看| 有码 亚洲区| 一夜夜www| 欧美xxxx黑人xx丫x性爽| 秋霞在线观看毛片| 丝瓜视频免费看黄片| 伊人久久国产一区二区| 亚洲精品一区蜜桃| 秋霞伦理黄片| 亚洲欧美日韩东京热| 夜夜看夜夜爽夜夜摸| 国内少妇人妻偷人精品xxx网站| 国产亚洲5aaaaa淫片| 亚洲国产精品专区欧美| 我要看日韩黄色一级片| 欧美日韩国产mv在线观看视频 | 国产精品麻豆人妻色哟哟久久 | 噜噜噜噜噜久久久久久91| 美女内射精品一级片tv| 亚洲va在线va天堂va国产| 狂野欧美白嫩少妇大欣赏| 老师上课跳d突然被开到最大视频| 一个人看视频在线观看www免费| 欧美精品一区二区大全| 美女主播在线视频| 日韩大片免费观看网站| 午夜福利视频精品| 一级黄片播放器| 午夜精品一区二区三区免费看| 一个人看视频在线观看www免费| 久久精品人妻少妇| 热99在线观看视频| 国产精品久久久久久精品电影| 国产成人福利小说| 校园人妻丝袜中文字幕| 午夜老司机福利剧场| 色视频www国产| 国产在线一区二区三区精| 少妇的逼水好多| 亚洲性久久影院| 99re6热这里在线精品视频| h日本视频在线播放| 三级男女做爰猛烈吃奶摸视频| 欧美xxⅹ黑人| 欧美日韩在线观看h| 久久久久久伊人网av| 丰满乱子伦码专区| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 欧美丝袜亚洲另类| 日本一本二区三区精品| 联通29元200g的流量卡| 久久精品国产亚洲av涩爱| 亚洲怡红院男人天堂| 亚洲av电影在线观看一区二区三区 | 又爽又黄a免费视频| 久久久久久久久久人人人人人人| 久久久精品94久久精品| 久久久久久久国产电影| 永久免费av网站大全| 搡老乐熟女国产| 亚洲精品中文字幕在线视频 | 成人漫画全彩无遮挡| 男人狂女人下面高潮的视频| 九草在线视频观看| av免费在线看不卡| 午夜激情福利司机影院| 久久久久国产网址| av又黄又爽大尺度在线免费看| 国产成人aa在线观看| 国产伦在线观看视频一区| 久久久久精品久久久久真实原创| 在线 av 中文字幕| 日韩欧美一区视频在线观看 | 特级一级黄色大片| 国产成人精品福利久久| 蜜臀久久99精品久久宅男| 精品一区二区免费观看| 成人漫画全彩无遮挡| 国产精品一区www在线观看| 高清毛片免费看| 国产成人一区二区在线| 亚洲精品aⅴ在线观看| 亚洲综合色惰| 搡老妇女老女人老熟妇| 一级黄片播放器| 精品国内亚洲2022精品成人| 欧美日韩综合久久久久久| 日韩中字成人| 99九九线精品视频在线观看视频| 亚洲av电影不卡..在线观看| 最近中文字幕2019免费版| 日韩av免费高清视频| 久久人人爽人人片av| 亚洲欧美日韩东京热| 国产精品一及| 欧美3d第一页| 国产精品一区二区性色av| 久久久精品欧美日韩精品| 晚上一个人看的免费电影| 全区人妻精品视频| 91久久精品国产一区二区三区| 久久久久久久大尺度免费视频| 大陆偷拍与自拍| 久久人人爽人人爽人人片va| 99热这里只有是精品在线观看| 国产精品一区二区三区四区免费观看| 国国产精品蜜臀av免费| 日本三级黄在线观看| 永久免费av网站大全| 视频中文字幕在线观看| 男女那种视频在线观看| 亚洲自拍偷在线| videossex国产| 青春草视频在线免费观看| 久久99蜜桃精品久久| 97热精品久久久久久| 久久精品熟女亚洲av麻豆精品 | 国产淫片久久久久久久久| kizo精华| 亚洲综合色惰| 亚洲精品456在线播放app| h日本视频在线播放| 有码 亚洲区| 嘟嘟电影网在线观看| 成人亚洲欧美一区二区av| 777米奇影视久久| 成人鲁丝片一二三区免费| 午夜免费观看性视频| 久久热精品热| 国产av国产精品国产| 美女高潮的动态| 少妇高潮的动态图| 淫秽高清视频在线观看| 成人美女网站在线观看视频| 男人舔女人下体高潮全视频| 亚洲精品成人av观看孕妇| 国产黄片美女视频| 国产精品精品国产色婷婷| 国产在视频线在精品| 亚洲精品视频女| 亚洲精品影视一区二区三区av| 99热这里只有是精品50| 久久精品国产自在天天线| 在线 av 中文字幕| 在线a可以看的网站| 亚洲国产欧美人成| 久久久久久久亚洲中文字幕| av国产免费在线观看| 久久久久网色| 成人亚洲欧美一区二区av| 亚洲最大成人av| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 大香蕉97超碰在线| 免费无遮挡裸体视频| 精品久久国产蜜桃| 久久久a久久爽久久v久久| 久久99蜜桃精品久久| 美女cb高潮喷水在线观看| 久久久久免费精品人妻一区二区| 亚洲av成人av| 丝瓜视频免费看黄片| 日韩av在线免费看完整版不卡| 六月丁香七月| 熟妇人妻久久中文字幕3abv| 久久久久久久午夜电影| 亚洲av成人精品一区久久| 日本午夜av视频| 国产淫片久久久久久久久| 国产精品.久久久| 久久精品综合一区二区三区| 亚洲在线观看片| 美女大奶头视频| 国产精品嫩草影院av在线观看| 精品少妇黑人巨大在线播放| 久久久久久久午夜电影| 成人综合一区亚洲| 丰满人妻一区二区三区视频av| 国产69精品久久久久777片| 国产精品久久久久久久电影| av在线天堂中文字幕| 男女那种视频在线观看| 亚洲怡红院男人天堂| 麻豆精品久久久久久蜜桃| 天堂中文最新版在线下载 | 毛片女人毛片| 国产白丝娇喘喷水9色精品| 噜噜噜噜噜久久久久久91| 十八禁国产超污无遮挡网站| 精品国产三级普通话版| 国产精品蜜桃在线观看| 久久久午夜欧美精品| 插阴视频在线观看视频| 少妇猛男粗大的猛烈进出视频 | 男的添女的下面高潮视频| 亚洲高清免费不卡视频| 日韩国内少妇激情av| 午夜激情久久久久久久| 老女人水多毛片| 国产 一区精品| 春色校园在线视频观看| 亚洲经典国产精华液单| 国产高清有码在线观看视频| av在线亚洲专区| 国产精品人妻久久久影院| 晚上一个人看的免费电影| 男人舔奶头视频| 一本久久精品| 黄色一级大片看看| 日本av手机在线免费观看| 秋霞伦理黄片| 国产国拍精品亚洲av在线观看| 亚洲图色成人| 精品国产三级普通话版| 成人亚洲精品av一区二区| 婷婷色麻豆天堂久久| 高清欧美精品videossex| 少妇熟女aⅴ在线视频| 国产精品日韩av在线免费观看| 国产亚洲av片在线观看秒播厂 | 99热6这里只有精品| 永久网站在线| 舔av片在线| 女的被弄到高潮叫床怎么办| 色5月婷婷丁香| 少妇人妻一区二区三区视频| 视频中文字幕在线观看| 中文字幕亚洲精品专区| 国产高清不卡午夜福利| 国产精品人妻久久久影院| 国产美女午夜福利| 成人午夜精彩视频在线观看| 国产女主播在线喷水免费视频网站 | 亚洲精品影视一区二区三区av| 国产亚洲av片在线观看秒播厂 | 亚洲欧洲日产国产| 五月天丁香电影| 亚洲人成网站在线观看播放| 看非洲黑人一级黄片| 国产免费视频播放在线视频 | 亚洲国产欧美在线一区| 中文字幕亚洲精品专区| 婷婷色综合大香蕉| 亚洲成人一二三区av| 22中文网久久字幕| av女优亚洲男人天堂| 国产男女超爽视频在线观看| 亚洲欧美日韩卡通动漫| 人妻制服诱惑在线中文字幕| 国产黄色小视频在线观看| 免费av观看视频| 日韩av免费高清视频| 国产精品麻豆人妻色哟哟久久 | 男女那种视频在线观看| 97超视频在线观看视频| 一级爰片在线观看| 美女黄网站色视频| 十八禁国产超污无遮挡网站| 日韩三级伦理在线观看| 成人美女网站在线观看视频| 中文精品一卡2卡3卡4更新| 久久精品国产鲁丝片午夜精品| 美女大奶头视频| 国产精品一二三区在线看| 日韩av在线免费看完整版不卡| 在线 av 中文字幕| 欧美精品一区二区大全| 国产成人aa在线观看| 三级国产精品欧美在线观看| 男人爽女人下面视频在线观看| 夜夜看夜夜爽夜夜摸| 欧美日韩国产mv在线观看视频 | 国产精品99久久久久久久久| 国产综合精华液| 亚洲欧美成人精品一区二区| 一级二级三级毛片免费看| 免费观看精品视频网站| 色哟哟·www| 中文字幕久久专区| 亚洲内射少妇av| 国产成人a区在线观看| 亚洲伊人久久精品综合| 国产大屁股一区二区在线视频| 国产精品久久久久久精品电影| 亚洲不卡免费看| 婷婷色综合大香蕉| 亚洲美女搞黄在线观看| 国产成人91sexporn| 特大巨黑吊av在线直播| eeuss影院久久| 欧美 日韩 精品 国产| 91av网一区二区| 亚洲一区高清亚洲精品| 三级国产精品欧美在线观看| 噜噜噜噜噜久久久久久91| 久久久久久久国产电影| 日韩人妻高清精品专区| 视频中文字幕在线观看| 久久鲁丝午夜福利片| 久久99热这里只频精品6学生| 特级一级黄色大片| 日日干狠狠操夜夜爽| 亚洲最大成人中文| 天天躁日日操中文字幕| 国产男女超爽视频在线观看| 男插女下体视频免费在线播放| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 亚洲第一区二区三区不卡| 国产亚洲精品久久久com| 亚洲欧美精品自产自拍| 精品人妻熟女av久视频| 亚洲精品久久久久久婷婷小说| 最近2019中文字幕mv第一页| 久久99精品国语久久久| 久久草成人影院| 哪个播放器可以免费观看大片| 婷婷色av中文字幕| 日韩精品有码人妻一区| av黄色大香蕉| 亚洲国产欧美在线一区| 午夜视频国产福利| 免费大片18禁| 国内精品一区二区在线观看| 久久韩国三级中文字幕| 熟女电影av网| 在线免费十八禁| 纵有疾风起免费观看全集完整版 | 欧美区成人在线视频| 亚洲精品国产av成人精品| av.在线天堂| 午夜老司机福利剧场| 国产伦一二天堂av在线观看| 七月丁香在线播放| 人妻少妇偷人精品九色| 国产午夜精品一二区理论片| 26uuu在线亚洲综合色| 亚洲成人久久爱视频| 网址你懂的国产日韩在线| 亚洲精品色激情综合| 高清毛片免费看| 人人妻人人澡欧美一区二区| 99久国产av精品国产电影| 国产精品1区2区在线观看.| 舔av片在线| 精品久久久久久久久久久久久| 亚洲第一区二区三区不卡| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 免费不卡的大黄色大毛片视频在线观看 | 老女人水多毛片| 亚洲欧美精品专区久久| 亚洲欧美中文字幕日韩二区| 免费观看av网站的网址| 搡老妇女老女人老熟妇| 极品少妇高潮喷水抽搐| 麻豆久久精品国产亚洲av| 美女cb高潮喷水在线观看| 伦理电影大哥的女人| 激情五月婷婷亚洲| 亚洲精品,欧美精品| 午夜福利在线观看吧| 久久亚洲国产成人精品v| 亚洲av成人av| 精品人妻一区二区三区麻豆| 黄片无遮挡物在线观看| 男插女下体视频免费在线播放| 亚洲,欧美,日韩| 五月玫瑰六月丁香| 成人二区视频| 国产一级毛片七仙女欲春2| 亚洲av电影在线观看一区二区三区 | 亚洲内射少妇av| 久久精品国产亚洲网站| av.在线天堂| 久久6这里有精品| 99久久精品一区二区三区| 日韩制服骚丝袜av| 免费看日本二区| 1000部很黄的大片| 免费观看的影片在线观看| 人妻制服诱惑在线中文字幕| or卡值多少钱| 一级av片app| 精品一区二区三区人妻视频| 超碰av人人做人人爽久久| 真实男女啪啪啪动态图| 人人妻人人澡欧美一区二区| 午夜福利视频1000在线观看| 欧美人与善性xxx| 天堂√8在线中文| 亚洲人与动物交配视频| 色网站视频免费| 亚洲国产欧美人成| 男女那种视频在线观看| 日本wwww免费看| 久久精品国产亚洲网站| 欧美日韩一区二区视频在线观看视频在线 | 午夜福利在线在线| 久久久久久久亚洲中文字幕| 深夜a级毛片| 国产美女午夜福利| 日本黄大片高清| 99热6这里只有精品| 婷婷色综合大香蕉| 一级爰片在线观看| 国产高清国产精品国产三级 | 青春草视频在线免费观看| 国产高清国产精品国产三级 | 久久亚洲国产成人精品v| 亚洲精品久久午夜乱码| 日本免费在线观看一区| 久久精品国产亚洲网站| 亚洲怡红院男人天堂| 日韩电影二区| 69人妻影院| 欧美最新免费一区二区三区| 好男人在线观看高清免费视频| 免费看a级黄色片| 亚洲精品乱码久久久v下载方式| 亚洲av福利一区| 91狼人影院| 亚洲高清免费不卡视频| 天天躁夜夜躁狠狠久久av| 超碰97精品在线观看| 久久久久久久久大av| 免费人成在线观看视频色| 一个人免费在线观看电影| 色综合色国产| 日韩亚洲欧美综合| 国产一区二区三区av在线| 中文字幕人妻熟人妻熟丝袜美| 男女边吃奶边做爰视频| 亚洲最大成人中文| 国产精品嫩草影院av在线观看| 自拍偷自拍亚洲精品老妇| 亚洲图色成人| a级毛片免费高清观看在线播放| 国产在线男女| 黑人高潮一二区| 久久久久性生活片| 看十八女毛片水多多多| 国产成人精品婷婷| 成人亚洲精品av一区二区| 黄片无遮挡物在线观看| 建设人人有责人人尽责人人享有的 | 亚洲精品亚洲一区二区| 99热6这里只有精品| 婷婷色麻豆天堂久久| 国产成人aa在线观看| 日韩人妻高清精品专区| 毛片女人毛片| 日本-黄色视频高清免费观看| 亚洲国产精品成人久久小说| 国内精品美女久久久久久| 精品人妻熟女av久视频| 日韩av在线免费看完整版不卡| 久久久久精品性色| 亚洲精品亚洲一区二区| 高清欧美精品videossex| 看免费成人av毛片| 国产亚洲最大av| 免费观看精品视频网站| 国产精品爽爽va在线观看网站| 麻豆成人av视频| 久久6这里有精品| 黄色日韩在线| 2021少妇久久久久久久久久久| 亚洲成人av在线免费| 亚洲精品影视一区二区三区av| 女人久久www免费人成看片| 日日摸夜夜添夜夜爱| 中国美白少妇内射xxxbb| 男女边吃奶边做爰视频| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 欧美激情国产日韩精品一区| 国精品久久久久久国模美| 伦理电影大哥的女人| 最后的刺客免费高清国语| 亚洲av日韩在线播放| 非洲黑人性xxxx精品又粗又长| 蜜臀久久99精品久久宅男| 晚上一个人看的免费电影| 亚洲美女视频黄频| 国产高清三级在线| 亚洲av免费在线观看| 免费观看精品视频网站| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 久久久久精品久久久久真实原创| 晚上一个人看的免费电影| 免费高清在线观看视频在线观看| 亚洲精品中文字幕在线视频 | 中国美白少妇内射xxxbb| 久久久久久久久久黄片| 丝袜喷水一区| 国产 亚洲一区二区三区 | 亚洲精品456在线播放app| 日本色播在线视频| 久久草成人影院| 日韩伦理黄色片| 日韩欧美国产在线观看| 国产精品人妻久久久影院| 国产伦精品一区二区三区四那| 亚洲精品国产成人久久av| 人人妻人人看人人澡| 国产精品女同一区二区软件| 精品一区二区三区视频在线| 国产白丝娇喘喷水9色精品| 能在线免费看毛片的网站| 一级毛片黄色毛片免费观看视频| 亚洲va在线va天堂va国产| 成人毛片60女人毛片免费| 秋霞在线观看毛片| 亚洲欧美日韩无卡精品| 99久久人妻综合| 亚洲熟妇中文字幕五十中出| 久久久久久久久中文| 青春草亚洲视频在线观看| 天天一区二区日本电影三级| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av涩爱| 亚洲aⅴ乱码一区二区在线播放| .国产精品久久| 久久午夜福利片| 日本黄大片高清| 亚洲欧美中文字幕日韩二区| 日韩av在线免费看完整版不卡| 麻豆精品久久久久久蜜桃| 麻豆久久精品国产亚洲av| 亚洲欧美中文字幕日韩二区| 人妻夜夜爽99麻豆av| 五月天丁香电影| 亚洲av国产av综合av卡| 欧美丝袜亚洲另类| 美女大奶头视频| 国产在线一区二区三区精| xxx大片免费视频| 免费看不卡的av| 搡老妇女老女人老熟妇| 国产高清不卡午夜福利| 久久久久精品性色| 汤姆久久久久久久影院中文字幕 | 国产探花极品一区二区| 大片免费播放器 马上看| 国产午夜福利久久久久久| 在线免费观看不下载黄p国产| 欧美日韩精品成人综合77777| 免费黄频网站在线观看国产| 精品久久久久久久久av| 国产成人一区二区在线|