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

空間信息技術在河口潮間帶棲息地分類研究中的應用

2016-03-04 08:31:56張婷婷侯俊利
海洋漁業 2016年4期
關鍵詞:物種分類特征

張婷婷,趙 峰,侯俊利,高 宇,莊 平

(中國水產科學研究院東海水產研究所,農業部東海與遠洋漁業資源開發利用重點實驗室,農業部東海與長江口漁業資源環境科學觀測實驗站,上海 200090)

·綜述·

空間信息技術在河口潮間帶棲息地分類研究中的應用

張婷婷,趙 峰,侯俊利,高 宇,莊 平

(中國水產科學研究院東海水產研究所,農業部東海與遠洋漁業資源開發利用重點實驗室,農業部東海與長江口漁業資源環境科學觀測實驗站,上海 200090)

研究河口潮間帶底棲生境(intertidal benthic habitat)變化,能為河口濕地生物的保護、底棲環境質量的評價、受損生境的生態恢復提供的科學依據和管理建議。本文從關注空間和時間尺度的研究出發,對河口潮間帶棲息地監測、分類與評價方面的研究進展進行了歸納總結,分為以下六個方面進行論述:河口潮間帶底棲生境的時空新特征、河口潮間帶底棲生境特征參數、河口潮間帶底棲生境特征獲取方式、數據處理、河口潮間帶底棲生境分類模型選擇、多時空尺度分析。河口潮間帶底棲生境特征的分類方式由于考慮生態過程和格局中所涉及的空間生態學問題(邊界分割、尺度效應、分類精度、動態變化等),模擬底棲物種及群落分布更為便捷,最終使物種/群落棲息地分類結果更加精確、真實、全面。最后,本文提出利用SIT對河口潮間帶底棲生境進行分類將對底棲生態系統評價及研究生態功能機制帶來新的視角,同時也需注意一些問題。例如,潮間帶生境特征參數需要選擇合適的統計方法在不同的尺度上進行甄別和篩選,以期有效提高指示參數在研究尺度上的解釋度。

空間信息技術;河口潮間帶;底棲生境;分類研究

潮間帶棲息地(intertidal habitat)是河口生態系統的重要組成部分,其與陸地、海洋生態系統的營養物質交換頻繁、環境復雜多變,具有支持植被、提供動物棲息所,維持生物多樣性的生態功能[1]。在全球氣候變化、上游水沙變化和外來物種入侵等自然因素以及促淤圈圍、深水航道工程、水庫建設等人為因素的共同作用下,高度敏感、脆弱的河口潮間帶棲息地是受污染物和人類活動影響最早和最直接的生態區,近幾十年來發生了巨大變化,其生態系統的功能和服務遭到嚴重破壞[2]。作為與河口潮間帶理化環境密切相關的底棲動物群落[3],其分布能有效反映河口生態系統的健康狀況[4],是指示河口潮間帶生境變化最直接的指示因子。研究和監測河口潮間帶底棲生境(intertidal benthic habitat)變化,能為河口濕地生物的保護、底棲環境質量的評價、受損生境的生態恢復提供科學依據和管理建議。

隨著空間信息技術(spatial information technologies,SIT)的發展,SIT技術逐漸囊括地理信息系統技術(geographic information system,GIS)、遙感(remote sensing,RS)、空間模型(spatial modeling)、空間統計學(spatial statistics)[5],利用SIT技術能獲取棲息地的空間信息及動態的時間信息,具有多樣性與多維性,使科學家們能得到更為精確、全面、系統的河口潮間帶底棲生境監測與評價結果[6]。

本文從關注空間和時間尺度的研究出發,以滿足宏觀河口管理的需求為目標,在現行SIT技術快速發展的背景下,對河口潮間帶棲息地監測、分類與評價方面的研究進展進行了歸納總結,分為以下六個方面進行論述:河口潮間帶底棲生境的時空新特征、河口潮間帶底棲生境特征參數、河口潮間帶底棲生境特征獲取方式、數據處理、河口潮間帶底棲生境分類模型選擇、多時空尺度分析。以期構建大尺度上河口底棲物理生境分類指標體系和分類方法,為河口底棲生境分類和基于生境的生態保護及水生態修復措施等管理策略的制訂提供依據。

1 河口潮間帶底棲生境的時空新特征

近20年來,廣大科研工作者已經在河口潮間帶生境變化與底棲動物群落的結構、數量動態、時空分布格局、次級生產力、食性和營養級、重金屬富集毒理等方面的關系,進行了廣泛而深入的研究。在這些研究中,研究者面臨的首要問題是對潮間帶底棲生境進行劃分和描述。

通常,“生境”或“棲息地”(habitat)的定義為維持特定物種或群落生存所需的空間和其中全部物理、化學、其它生物因子的總和[7-8]。對棲息地定義的理解,導致劃分棲息地類型的研究大多直接依據環境中的生物或非生物因子。研究者對河口潮間帶底棲生物棲息地的劃分,往往也是根據不同的生物及非生物因子分類規則,進行比較研究,分類規則大體分為三類[9]:一、根據某一主要生態或環境因子分類,如鹽沼植物生境(海草、紅樹林)[10-11]、潮溝生境[12-13]、鹽度梯度生境[14]、潮上中下帶生境[15];二、綜合一組河口潮灘濕地中差異性最大的環境因子[16],或者一組對生物群落影響最大的環境因子進行分類,如構成生境適宜性指數(habitat suitability index,HSI)[17]、環境指數(index of environmental condition,IEC)[18]、生態整合度指數(index of ecological integrity,IEI)[19]進行生境分類;三、劃分為人工干擾生境與自然生境,如大型航道工程生境[20]、土方工程生境[21]、圍墾生境[22]。

然而,這三種分類方式都有一定的局限性。在小區域范圍內進行短期靜態研究時,或影響目標物種生存的棲息地特征因子較為簡單的情況下,以上三種方法是合適的[23]。但是河口是河流與海洋的結合地段,受河流終段和近海水域相互影響,因此,河口潮間帶是一類空間高度異質化的復雜體,包括光灘、潮汐、鹽沼植被等,其生物和非生物因子隨著不同的生境有著很大的變化,如植物、水文、地形、土壤、鹽度、浮游生物等。另外,受到潮汐、風暴及人類活動的影響,這些物理生化環境因子會在很小的空間范圍內(有時小于100 m)和時間尺度上發生變化[24]。因此,對如此諸多、動態變化的環境因子進行大規模采集存在人力、物力的限制,傳統的棲息地定義及隨之產生的分類規則亦不能全面滿足監測和評價的需求,亟需從時空角度上考慮棲息地新特征,確定新的分類規則和采樣手段,對復雜多變的河口潮間帶底棲生境進行分類,用以滿足長期、便捷、更加精確監測其動態變化的需求[8]。

近年來,隨著SIT技術的發展以及對空間生態學概念理解的加深,可以利用GIS對棲息地類型特征或空間單元(斑塊)賦予空間變量特征,記錄其在范圍、面積、形狀、斑塊位置、與鄰近斑塊空間關系的變化等空間結構和紋理信息[8]。因此,基于棲息地斑塊的內在屬性和空間性質[25],可以將棲息地特征重新劃分為三類:(1)屬性特征(typology):指生物和非生物因子特征,如地形、地貌、植物群落組成、植物優勢度、植物演替階段、其它動物群落組成等。通過棲息地斑塊屬性特征的相似性,或基于此運算得到的棲息地指數確定不同棲息地類型(以往的三類分類研究屬于基于該類特征)。(2)紋理特征(texture):關注每一種棲息地類型的斑塊數目、大小及其頻度分布。當每種棲息地類型中斑塊數量和面積大小計算出來時,即可對紋理特征進行量化。(3)結構特征(structure):由棲息地斑塊的空間結構(范圍、面積、形狀)、布局特征以及斑塊之間的空間地理關系(斑塊位置、與鄰近斑塊空間關系的變化)確定。空間結構可以由多種空間統計方法或空間指數描述(如破碎化指數、景觀指數等)。棲息地斑塊空間性質、紋理和結構特征能反映多種環境因子下不同時空尺度上生境的綜合狀況。通過描述棲息地的紋理和結構特征這類新屬性,研究其與生物種群、群落之間的關系,能從與棲息地水平相關的多尺度條件下,描述和理解生物多樣性狀態和變化機制。

2 河口潮間帶底棲生境特征參數

一般河口潮間帶底棲生境的狀態是通過對上述的屬性特征(物理-化學性質、物種組成和/或種群相對豐度)以及紋理、結構空間特征(斑塊面積、分布、數目、破碎化程度等)進行采集而評估的。除了生境斑塊的內在特征外,特征信息采集通常還涉及一些額外信息,如水文、生物群落生理生化指標、以及潛在的驅動因子和人類脅迫因子(土地利用方式等)。表1根據以往文獻總結了這些生境特征數據組及額外信息的類型以及獲取方式[1,26-28]。

3 河口潮間帶底棲生境特征獲取方式

根據主要的數據來源,河口潮間帶底棲生境特征獲取方式主要分為野外采樣和輔以地面驗證數據的遙感監測方式。

3.1 野外采集生物及底泥信息

野外采樣是指基于野外調查和測量,通常采用固定的樣點或樣條,在小尺度上探查棲息地中物種組成或相對豐度。在未使用SIT技術之前,一般通過對某些區域實施抽樣采樣的調查方法,通過先驗的隨機法或后驗的空間模型法(外推法),推測其它區域物種組成或相對豐度。野外采樣能有效記錄局部和區域尺度上底棲物種空間分布結構以及其環境決定因子空間變化。然而,野外采樣的方法十分耗時耗力,同時還需專業的分類學家參與。因此,野外采樣制圖,很少能覆蓋完整的研究區域,并且小尺度上的異質性格局通常被忽略。單單通過野外采樣,不可能設計出棲息地高分辨率的采樣方案,更不能實現對棲息地進行長期監測的目標[1]。

隨著遙感技術的發展,潮間帶棲息地的大多環境屬性特征能得到有效地大范圍、連續精確監測[8](見3.2),但其底棲環境的大部分生物特征和底泥特征不能直接被遙感技術測量。因此,為了闡明底棲生物群落及其生境之間的關系,利用原位采樣的方法,采集底棲環境的生物和底泥特征必不可少(表1)。

3.2 遙感監測生境環境屬性特征

利用遙感技術能對河口潮間帶的多種環境因子進行識別,從而實現對棲息地多種環境因子大規模連續監測的目標。近幾十年來,已有研究者利用多種遙感源對河口棲息地環境特征進行監測[1],例如,航拍圖像、衛星影像(潮間帶最常用的衛星是Quickbird、IKONOS、SPOT、Landsat TM/ETM),高光譜圖像、聲學遙感。河口潮間帶的多個主要環境特征已得到了成功地提取和反演,如高程[29]、植被[30]、土地利用[31]、潮溝[32]、土壤[33](表1)。

表1 河口潮間帶底棲生境特征參數及其主要獲取方式Tab.1 Characteristics of intertidal benthic habitat in estuary and their acquisition methods

與野外采樣相比,遙感有明顯的優勢。例如避免對河口潮間帶造成物理侵入或人為的干擾、不需要太多的生態專業知識[8]、可以得到連續覆蓋的大尺度環境數據層,尤其是高分辨率數據,能在更小的尺度上辨識和監測高度破碎化的潮間帶底棲生境系統[34]。

在大部分的研究中,河口潮間帶底棲生境圖是通過遙感環境數據和原位采樣數據之間的關系推譯出來的。在這種方式下,河口潮間帶底棲生境分布估測圖最終是由完整的環境數據作為指示因子所繪制成圖。

3.3 Fragstats軟件提取生境空間特征

棲息地斑塊的空間特征與種群和群落的多個生態過程和分布格局密切相關,如物種關系、繁殖、分散、干擾等。ZIV[35]認為棲息地類型越多越復雜,即多樣性(diversity)程度高,將為更多的生物種群創造更多的生存機會及更大的棲息地選擇權;棲息地斑塊大小(patch size)是決定種群數量及滅絕速率的重要因子,斑塊面積越大,可支持的種群數量越多、滅絕速率越低;棲息地斑塊性(patchiness)越大,即同一棲息地類型中斑塊數量越多,一個群落中個體在整個區域生存機會越大,滅絕速率越低,種群的穩定性越高。

遙感易于辨識各種棲息地斑塊的空間特征,如紋理和結構特征,然后通過簡單的運算得到景觀指數(landscape metrics)進行量化。利用Fragstats軟件工具可以便捷地通過遙感專題圖計算河口潮間帶生境景觀指數,以代表生境在不同尺度上的景觀結構、形狀、屏障、連接度、邊緣、多樣性等生境特征(表1)。通過量化的景觀指數,對棲息地進行分類,例如分為破碎性生境、連續性生境、邊緣銳化生境、邊緣平滑生境、多屏障生境、多廊道生境。目前,研究者已經就應用景觀特征描述河口潮間帶生境格局和動態變化方面進行了相關研究[36-38],并通過景觀指數定量分析景觀特征對大型底棲動物群落結構和分布的影響,證實了河口潮間帶生境斑塊大小、斑塊數目、斑塊組成和斑塊連接度對大型底棲動物物種多樣性、物種優勢度、物種豐度有顯著影響[39-41]。例如,DE JUAN等[40]通過提取新西蘭八大河口區六大生境類型的生境破碎化指數、斑塊大小和數目等景觀指數,將其與物種多樣性指數進行相關分析。結果表明,斑塊大小對河口大型底棲動物物種多樣性有顯著影響。BEATRIZ等[41]深入研究了巴西圣保羅北部海岸生境破碎化后,邊緣景觀、斑塊大小、斑塊組成對潮下帶底棲動物群落的影響。研究表明,生境破碎化后盡管小斑塊的邊緣效應增加了物種的豐度,但物種優勢度和多樣性大大降低。THRUSH等[42]通過研究生境斑塊數目、斑塊大小、斑塊連接度等指數與大型底棲動物群落豐度、總密度和種群組成的關系。結果表明,物理環境上的景觀連接度對大型底棲動物群落結構的恢復起了重要作用。

4 數據處理

隨著SIT技術的發展,一方面關于棲息地狀態的數據來源更加豐富,另一方面,處理這些數據的手段也更為多樣。因此,研究者可以更加全面考慮生態過程和格局中所涉及的空間生態學問題,最終使物種/群落棲息地分類結果更加精確、真實、全面。

一個全面的數據處理程序應該包括以下步驟[43]:一、整合相關數據,使之可用于同一分析過程中,包括多源遙感數據融合、遙感數據與野外采樣數據整合;二、評價數據的充分性(數據精確性、全面性、完整性),包括生物數據的零值處理、遙感專題圖分類精度優化、生態特征變量代表層次選擇;三、處理高度相關的環境變量——空間自相關分析;四、棲息地連續特征變量的邊界分割。

4.1 數據整合

由于棲息地特征獲取方式多樣,數據來源多元化,在棲息地分類制圖前,需要對元數據進行數據整合。通常有二類數據整合方式:一、不同來源的遙感數據融合。不同分辨率的原始遙感數據代表的尺度不同,例如小尺度數據無法簡單代表大尺度數據[8],因此,由空間、光譜分辨率差異導致的尺度問題對棲息地分類有較大影響。已有大量文獻對不同遙感來源的數據融合方法進行了綜述[1,44],主要有象元(pixel)、特征(feature)和解譯(interpretation)水平的融合;二、遙感數據與野外采樣數據整合[44]。這類整合的優勢為能將棲息地特征進行互補,利用高精度的野外采樣數據用以支持遙感數據或對遙感解譯圖像進行驗證,或提供遙感所不能獲取的底棲環境額外信息,如土壤信息等。另外,遙感能用作野外采樣的空間信息的補充數據或校正依據,提供在野外采樣中難以探查的棲息地空間格局信息,如破碎化、連接度等。但野外采樣數據大多為位點上的離散數據,需要在GIS平臺上對其進行柵格化或矢量化后,與連續的遙感數據進行整合。

4.2 零值處理

在對底棲種群或群落的野外采樣中常常面臨大量的零值數據,影響隨后的科學分析。零值數據的來源主要有三類[45]:一、真實零值,表示研究區內真的不存在某一物種;二、假零值,也稱為偽零值,這類零值是對物種較低的采樣概率所導致;三、采樣干擾零值,已知某一物種常出現在某地,但由于侵入式的采樣手段導致其驚嚇躲避從而無法記錄采樣。目前已有很多方法用于處理這三類采樣時產生的零值[46]。最簡單的方法是在數據處理前剔除零值。但簡單地去除零值常常會對結果產生影響,使分析結果出現偏差[46-47]。另一種方法是把所有的觀測值加一個正值的常量,接著進行對數函數(Log)轉化。這一方法要求人為選擇一個常量,這個常量的選擇不當將造成估計值偏差[47]。有研究對此方法進行改進,利用觀測值的標準分布概率函數(Delta對數正態或Gamma函數)得到一個二次型分布正值代替該常數。此法稱為Delta法[48],由于這一方法中常量是非零值觀測值的產物,大大減少了估值偏差。

4.3 遙感專題圖分類精度優化

遙感解譯得到土地利用/土地覆蓋專題圖是運算景觀指數的基礎,遙感影像專題圖分辨率(由土地覆蓋分類策略決定)代表地理空間信息的詳細程度,對景觀元素(斑塊等)分類及景觀指數分類有重要影響[49]。改變分類圖像的專題分辨率同樣將導致景觀類別數目和空間格局的變化,從而導致景觀指數的變化,許多研究結果證明了專題圖分辨率對大多數景觀指數有顯著的影響[49-51]。解決這一問題的方法是在允許的條件下,盡可能使用高空間分辨率的遙感影像進行解譯得到土地利用/土地覆蓋專題圖。高空間分辨率影像具有豐富的空間信息。因此在高分辨率影像上,同一地物表現出更大的差異,得到的景觀特征更加豐富。但傳統的基于像元光譜信息的分類方法會由于同一地物的光譜異質性造成嚴重的“椒鹽效應”(在一個小區域內,同種地物被分類成不同的類別,原本均一的地塊被“打碎”的效應),目前取而代之的是面向對象分類方法對高空間分辨率遙感影像進行分類。該分類法能夠結合形狀、紋理、拓撲、層次結構等空間信息,較其它空間分辨率遙感影像,可將河口潮間帶植被分類精度提高10%~25%[30]。

4.4 生境特征變量層級選擇

隨著SIT技術和空間生態學的發展,生態學家已經認識到物種和群落的分布依賴于尺度變化和在特定尺度上的棲息地異質性。在不同的尺度上,影響生物種群或群落結構和分布格局的關鍵棲息地特征不盡相同,往往從大尺度到小尺度上,這些關鍵性的棲息地特征之間具有層級性(hierarchical)關系[52]。MACEDO等[52]列舉了在河口流域尺度和局域尺度影響底棲動物群落結構與分布的棲息地特征因子,如局域尺度上重要棲息地特征為潮溝、底泥、水流、水化學和物理指標、水生植被、人為捕撈因子;而在流域尺度上的重要棲息地特征為地形、地貌、氣候、自然植被、草地、農用地、城市、人口密度、道路。因此,在選擇和處理影響生物種群或群落結構和分布的棲息地所有潛在特征因子時,需要根據不同的研究目的,選擇代表不同尺度下合適等級的生境特征類別。或者預先通過尺度效應分析,明確不同尺度下影響生物種群或群落結構和分布的特定棲息地特征,再行選擇分析。

4.5 空間自相關分析

生物數據(如豐富度)以及生境特征數據在地理空間上具有相關性,即空間格局本身不獨立,存在由相鄰個體變量的相互作用而產生的空間依賴性,或稱為空間自相關性[53-54]。若不能處理好數據的這一特征,將對模型參數估值及其不確定性造成影響,如棲息地模型的結果將產生空間相關的殘差、有偏差的誤差變量、錯誤的顯著水平和過高的回歸系數[55],最終導致不正確的統計推斷。因此,在選擇棲息地模型進行分類研究時,需要選擇能處理空間自相關性的模型。目前較為常用的、整合考慮了空間自相關性的種群棲息地模型包括[43]:自回歸方法、基于克里金插值(Kriging)的地統計法、廣義線性混合模型,廣義評估函數、地理權重回歸法等。

4.6 棲息地連續特征變量的邊界分割

隨著SIT技術的發展,高度異質性的潮間帶底棲生境的大部分環境特征能被遙感所探查,所解譯的連續環境數據層能在GIS系統中進行空間分析和處理。為了能清晰地表示生物及其棲息地的空間分布格局,在對棲息地不同類型進行劃分時需要確定相鄰類別之間的邊界,這意味著需要對連續的環境數據層進行分割(segmentation)。以往,一般相鄰的棲息地類型之間都是由人為劃分的清晰邊界來分割。然而,在現實情況下,河口棲息地特征之間的邊界并非總是清晰明確的(例如巖礁和周圍軟泥的相交地帶),尤其是底棲生境的環境特征經常呈梯度式變化。棲息地特征之間的分割邊界可以用更為真實的梯度方式來表達。新的空間統計方法,可以實現利用空間連續的環境數據組以連續或離散的方式描述和預測河口潮間帶底棲生物格局。這些方法更容易構建潮間帶底棲群落形成的真實情況,并避免了以往僵化固定的離散邊界。常用的邊界分割方法[7]有專家解譯、主成分分析、神經網絡法、貝葉斯決策法則、監督分類等。

5 河口潮間帶底棲生境分類模型選擇

針對不同研究對象,需要選擇合適的模型,對棲息地生境特征與生物分布數據進行擬合并模擬分類。已有的河口潮間帶底棲生境分類研究可以分為兩類[7]:單一底棲物種棲息地分類(single benthic species habitatmapping),底棲群落棲息地分類(benthic communitymapping)。

5.1 單一底棲物種棲息地分類

單一底棲物種棲息地分類所得到的棲息地類型特征必須本質上與目標物種支持環境狀況相匹配。因此理論上,生成底棲生境地圖最簡單的方法是將環境數據組作為指示因子與物種分布數據進行模擬,確定目標物種的棲息地。這是一種先整合、再預測的非監督分類策略,采取的是一種從上至下的方法,在現有的研究中使用最為廣泛。

隨著GIS和RS的快速發展,研究人員更易獲取更高質量和詳細的環境數據層,這些環境數據層能夠在質量上、空間范圍或分辨率上滿足目標物種的估測。在應用非監督分類策略時,這些連續的環境數據層預先通過專家解譯或目標分類的方法分割(見4.6)為棲息地“空間單元”,然后再利用采樣點上獲取的物種的生物數據(正常值/缺省值)進行建模分析。建模過程通常利用的是最簡單的相關統計方法。根據生物與環境數據組之間的相關關系,對生物物種的出現值(分布值)進行外推(圖1-A)。

圖1 底棲棲息地分類策略Fig.1 Strategies of benthic habitatm apping

另一種分類策略是先預測,再整合的監督分類策略,這是一種從下至上的方法,較前一種方法更加復雜和有針對性。在此,地面采樣所獲得的目標物種生物數據(出現值)將作為環境數據預測值的函數,通過物種分布模型(species distribution modelling,SDM)[43],并生成一張物種分布預測圖。在物種分布預測圖的基礎上,原位采樣所獲得的生物/地理地面驗證數據被用作“空間單元”(此處的空間單元為棲息地類型)分割的標準,將物種分布預測圖劃分物種棲息地分類圖(圖1-B)。SDM是一類能聯系物種分布數據(在已知地點的出現值/豐富度)和這些地點環境或/及空間特征信息的數學模型。相比傳統的多元線性回歸模型等[43],SDM能在時間和空間上進行外推,利用不充足的數據,實質上能等同表達物種所處環境的物種擬合模型,SDM模型能提供有用的生態學解釋和顯示強有力的預測能力。

一般情況下SDM方法都是類似回歸的,假設一個物種的出現值或豐富度可以通過環境因子的加性結合模擬,或通過表示環境因子之間關系的參數所模擬。SDM模型最常用的模型方法為廣義加性模型(generalized addictive model,GAM)[43]。由于其能描述非線性的生態響應(二次方、冪方或其它參數轉換),在擬合生態過程中提供了除線性關系以外更為真實的靈活度。另外,利用機器學習和數據挖掘的模擬方法[43],也是近來 SDM發展的熱點。如神經網絡法(artificial neural networks,ANNs)、多變量適應回歸曲線法(multivariate adaptive regression splines)、分類和回歸樹(classification and regression trees)、最大熵值模型(maximum entropy models)等。這些模型中有一些能提供很好控制的變量選擇和系數估計,另一些能自動探測和擬合預測因子之間的關系。為評價模型的模擬結果及避免模型偏差,一般使用CCI值、kappa系數、AIC(akaike information criterion,AIC)值用于表示模擬結果(生物群落分布、棲息地類型)和實測數據之間的擬合程度。已有研究表明,SDM可以成功地用于單一底棲物種棲息地分類研究[57-58],它們預測行為的最終結果能優于大部分傳統的方法[56]。

5.2 底棲群落棲息地分類

群落棲息地分類步驟與單一物種棲息地分類類似,同樣有兩種策略(圖1),如非監督分類策略中,群落生物數據通常以組別的形式存在(由多元統計、聚類分析等方法得到),進而與預先分割的環境空間單位進行比較分析,并對群落類別進行外推。

監督分類策略中,進行群落棲息地分類時,理論上可以通過每次生成一種物種棲息地分布圖,經過多次后最終將這些單一物種棲息地分布圖整合疊加為群落棲息地分布圖。

另一種監督分類常用的方法是,首先根據多元統計分析、專家判斷、或者已有的生物區系分類系統,對原位采樣生物數據進行群落分類。然后,利用這些群落類別以監督分類的形式對連續覆蓋的環境數據層進行分割,構成不同的群落棲息地類型。這類群落棲息地分類還可以引入一個群落生境(biotope)的概念,即“一個非生物環境組成的棲息地及其相關的生物群落的結合”[59]。這一術語能在空間上有效表達群落信息。從近年來大量棲息地分類研究來看,亟需發展基于群落生境概念的棲息地分類策略[43],這樣的棲息地分類規則具有層級結構,如高層級由大范圍的非生物因子(地理、地質因子)構成,低層級由在生物區系水平的物種信息構成,這樣所形成的棲息地分類系統既不會太粗略,又不會太過于細致(過多種群的細節信息),不用太過考慮集群結構的季節性變化,更有效地達到和實施保護目的[59]。

在潮間帶區域,棲息地以梯度式變化環境特征為主,并且在其非生物環境因子劇烈變化(潮水波動)、小尺度環境異質性極高的區域,群落格局容易被掩蓋,這一分類策略由于可以采用群落生境進行群落棲息地分類,其層級結構可以實現在不同空間尺度上識別底棲群落格局,但目前在潮間帶的研究還鮮有報道。

6 多時空尺度分析

隨著對物種和群落的棲息分布存在尺度效應這一認識的發展,對于復雜多變的河口潮間帶生態系統,在單一空間尺度分析棲息地特征與底棲生物群落分布之間的關系,已經不能滿足生態管理和保護的要求,空間層次法(spatial hierarchical approach)中多空間尺度分析(multiple scales analysis)可以用于解決這一問題[60]。

多空間尺度分析是根據系統的層級組織特征,將大尺度由上自下分為多個層級(level)小尺度的方式,在各個層級的空間尺度下,測試變量敏感性的方法[60]。采用多空間尺度的分析方法避免了人為隨意確定研究物種或群落區域尺度的問題。可量化、多層次的多個空間尺度范圍使分析群落與環境之間的關系具有可比性,可以明確群落在不同空間尺度上格局與過程中的變化和不確定性,得到更為全面和精確的結果和評價。

尤其是利用景觀指數作為棲息地空間特征的指示參數時,必須進行尺度效應分析。已有大量研究圍繞景觀指數在不同尺度上與大型底棲動物群落的關系展開,證明了在不同尺度上,潮間帶底棲動物群落對生境景觀特征的響應程度不同。THRUSH等[61]分別在樣點和景觀尺度上分析了新西蘭北島卡沃灣的景觀特征與底棲生物群落多樣性的相關關系。結果表明,在樣點尺度上,局部的景觀連接度對物種多樣性影響顯著,但影響程度較弱;在更廣的景觀尺度上,生境的連接度對物種多樣性增加效果更為顯著且有力。POTTER等[62]分別在流域尺度和河岸區域尺度上對北卡羅來納州73個站點提取景觀特征,并研究其與底棲動物群落的關系。結果表明,流域尺度上提取的景觀特征(解釋度為69.5%~75.4%)比河岸區域尺度上提取的景觀特征(解釋度為57.4%~65.2%)能更好地解釋大型底棲動物群落結構的變化。綜上所述,量化景觀指數,分析棲息地空間特征與底棲動物群落的關系必須在多時空尺度下進行。

另一個方面,大多重要的大型底棲物種為運動型(如蟹類、蝦類、魚類),分布范圍廣,且在不同的生活史階段(如覓食、避難、繁殖和育幼等)以不同的方式利用潮間帶棲息地資源。為了生成合理、準確的棲息地地圖,亟需深刻理解這些重要的底棲物種生活史階段棲息地。這些地圖需要進行多時間尺度上(月變化、季節變化、年變化)的動態分析,用于反映由季節或生物個體自身遺傳引起的棲息地偏好性在生活史各個階段的時空動態變化。然而,已有一些研究率先開展了對潮間帶底棲生物(蟹[63]、蝦[64]等)的長期動態監測,探索其棲息地時空動態選擇及變化規律。

7 小結與展望

SIT技術和空間生態學的發展,使獲取棲息地的空間信息及動態的時間信息具有多樣性與多維性。例如,棲息地類型特征或空間單元(斑塊)被賦予空間變量特征;通過遙感技術,大多潮間帶環境屬性特征能得到有效地大范圍、連續精確監測,河口潮間帶底棲生境特征能得到更為全面、系統、精確地描述。另一方面,數據來源豐富多樣,數據處理手段更為多樣,分類方式可以更加全面考慮生態過程和格局中所涉及的空間生態學問題(邊界分割、尺度效應、分類精度、動態變化等問題)。最后,空間模型和群落生境概念的發展使得底棲物種及群落分布的模擬更為便捷,最終使物種/群落棲息地分類結果更加精確、真實、全面。

潮間帶底棲棲息地分類對潮間帶底棲生態系統進行生態評價十分重要。例如,在對潮間帶底棲棲息地進行精確分類的基礎上,可以選擇與大型底棲生物群落棲息地質量密切相關的主要環境因子作為評價指標。應用棲息地質量量化模型對影響因子進行比較并賦予權重和分值以建立潮間帶底棲棲息地質量評估體系。

另外,在潮間帶底棲棲息地分類的基礎上,還可以結合潮間帶棲息地多樣性、底棲生物多樣性與食物網結構等關系,研究河口潮間帶底棲棲息地隨著環境梯度、地理空間上的變化,在局部和區域尺度上其生態功能的變化機制。

顯而易見,利用SIT對河口潮間帶底棲生境進行分類將對底棲生態系統評價及研究生態功能機制帶來新的視角。然而,利用SIT技術目前也存在一些問題。例如,利用景觀指數作為棲息地空間特征的指示參數時,不是所有的參數都具備顯著性的指示作用,需要選擇合適的統計方法在不同的尺度上進行甄別和篩選,才能有效提高指示參數在研究尺度上的解釋度。另外,由于新發展的物種分布模型等空間統計方法較傳統的數學方法更為復雜,現階段在數據處理和模型模擬過程中,這些優勢明顯的方法在河口潮間帶底棲生態系統監測和評價研究中還未能大范圍進行應用,亟待發展和推廣。

參考文獻:

[1] GODET L,FOURNIER J,TOUPOINT N,et al.Mapping and monitoring intertidal benthic habitats:a review of techniques and a proposal for a new visual methodology for the European coasts[J].Progress in Physical Geography,2009,33(3):378-402.

[2] ZHARIKOV Y,SKILLETER G A,LONERAGAN N R,et al.Mapping and characterising subtropical estuarine landscapes using aerial photography and GIS for potential application in wildlife conservation and management[J].Biological Conservation 2005,125(1):87-100.

[3] MCLUSKY D S.The estuarine ecosystem[M].New York:Chapman&Hall,1989.

[4] 沈新強,陳亞瞿,羅民波,等.長江口底棲生物修復的初步研究[J].農業環境科學學報,2006,25(2):373-376.

SHEN X Q,CHEN Y Q,LUO M B,et al.Preliminary study on the restoration of benthos in the Changjiang Estruary[J].Journal of Agro-Environment Science,2006,25(2):373-376.

[5] HESS G R,CHESHIRE H M.Spatial information technologies-integrating the tools with the curricula[J].Journal of Forestry,2002,100(1):29-34.

[6] 張婷婷,張 濤,侯俊利,等.空間信息技術在漁業資源及生態環境監測與評價中的應用[J].海洋漁業,2014,36(3):272-281.

ZHANG T T,ZHANG T,HOU J L,et al.Application of spatial information technologies in monitoring and assessment of fishery resources and ecological environment[J].Marine Fisheries,2014,36(3):272-281.

[7] BROWN C J,SMITH S J,LAWTON P,et al.Benthic habitat mapping:A review of progress towards improved understanding of the spatial ecology of the seafloor using acoustic techniques[J].Estuarine Coastal and Shelf Science,2011,92(3):502-520.

[8] LENGYEL S,KOBLER A,KUTNAR L,et al.A review and a framework for the integration of biodiversity monitoring at the habitat level[J].Biodiversity and Conservation,2008,17(14):3341-3356.

[9] SHOKRI M R,GLADSTONE W.Limitations of habitats as biodiversity surrogates for conservation planning in estuaries[J].EnvironmentalMonitoring and Assessment,2013,185(4):3477-3492.

[10] 童春富,章飛軍,陸健健.長江口海三棱藨草帶生長季大型底棲動物群落變化特征[J].動物學研究,2007,6(6):640-646.

TONG C F,ZHANG F J,LU J J.Variation characteristics of the macrobenthic fauna community in the Scirpusmariqueter zone of the Yangtze Estuary during the growing seasons[J].Zoological Research,2007,6(6):640-646.

[11] 謝志發,何文珊,劉文亮,等.不同發育時間的互花米草鹽沼對大型底棲動物群落的影響[J].生態學雜志,2008,27(1):63-67.

XIE Z F,HE W S,LIU W L,et al.Influence of Spartina alterniflora salt marsh at its different development stages on macrobenthos[J].Chinese Journal of Ecology,2008,27(1):63-67.

[12] RODRIGUESA M,QUINTINO V,SAMPAIO L,et al.Benthic biodiversity patterns in Ria de Aveiro,Western Portugal: Environmental-biological relationships[J].Estuarine Coastal and Shelf Science,2011,95(2-3):338-348.

[13] 宋慈玉,儲忝江,盛 強,等.長江口鹽沼分級潮溝系統中大型底棲動物群落結構特征[J].復旦學報(自然科學版),2011,50(3):253-259.

SONG C Y,CHU T J,SHENG Q,et al.Macrobenthos communities in salt marshes of the Yangtze River Esturay along a creek-order gradient[J].Journal of Fudan University(Natural Science),2011,50(3):253-259.

[14] BARNESR S K.Distribution patterns of macrobenthic biodiversity in the intertidal seagrass beds of an estuarine system,and their conservation significance[J].Biodiversity and Conservation,2013,22(2):357-372.

[15] 全為民,趙云龍,朱江興,等.上海市潮灘濕地大型底棲動物的空間分布格局 [J].生態學報,2008,28(10):5179-5187.

QUANW M,ZHAO Y L,ZHU J X,et al.The spatial pattern ofmacrozoobenthic communities in the tidal wetlands of Shanghai City[J].Acta Ecologica Sinica,2008,28(10):5179-5187.

[16] GAEVER S V,OLU K,DERYCKE S,et al.Metazoan meiofaunal communities at cold seeps along the Norwegian margin:Influence of habitat heterogeneity and evidence for connection with shallow-water habitats[J].Deep-Sea Research Part I-Oceanographic Research Papers,2009,56(5):772-785.

[17] DEGRAER S,VERFAILLIE E,WILLEMSW,et al.Habitat suitability modelling as a mapping tool for macrobenthic communities:An example from the Belgian part of the North Sea[J].Continental Shelf Research,2008,28(3):369-379.

[18] ENGLE V D,SUMMERS J K,GASTON G R.A benthic index of environmental-condition of Gulf of Mexico estuaries[J].Estuaries,1994,17(2):372-384.

[19] KASHIAN D R,BURTON T M.A comparison of macroinvertebrates of two Great Lakes coastal wetlands:Testing potential metrics for an index of ecological integrity[J].Journal of Great Lakes Research,2000,26(4):460-481.

[20] LERCARID,DEFEO O.Variation of a sandy beach macrobenthic community along a human-induced environmental gradient[J].Estuarine Coastal and Shelf Science,2003,58(3):17-24.

[21] 章飛軍,童春富,謝志發,等.長江口潮間帶大型底棲動物群落演替[J].生態學報,2007,27(12):4944-4952.

ZHANG F J,TONG C F,XIE Z F,et al.The recolonisation progress of intertidal benthic fauna community in the Changjiang Estuary[J].Acta Ecologica Sinica,2007,27(12):4944-4952.

[22] NASER H A.Effects of reclamation on macrobenthic assemblages in the coastline of the Arabian Gulf:A microcosm experimental approach[J].Marine Pollution Bulletin,2011,62(3):520-524.

[23] DIAZR J,SOLANM,VALENTE RM.A review of approaches for classifying benthic habitats and evaluating habitat quality[J].Journal of Environmental Management,2004,73(3):165-181.

[24] BENDELL L I,WAN P C Y.Application of aerial photography in combination with GIS for coastal management at small spatial scales:a case study of shellfish aquaculture[J].Journal of Coastal Conservation,2011,15(4):417-431.

[25] FIRBANK L G,BARR C J,BUNCE R G H,et al.Assessing stock and change in land cover and biodiversity in GB:an introduction to Countryside Survey 2000[J].Journal of Environmental Management,2003,67(3):207-218.

[26] BRUNS D A.Macroinvertebrate response to land cover,habitat,and water chemistry in a miningimpacted river ecosystem:A GISwatershed analysis[J].Aquatic Sciences,2005,67(4):403-423.

[27] JENSEN O P,CHRISTMAN M C,MILLER T J. Landscape-based geostatistics:a case study of the distribution of blue crab in Chesapeake Bay[J].Environmetrics,2006,17(6):605-621.

[28] RAKOCINSKIC F.Evaluatingmacrobenthic process indicators in relation to organic enrichment and hypoxia[J].Ecological Indicators,2012,13(1):1-12.

[29] 韓 震,惲才興,戴志軍,等.淤泥質潮灘高程及沖淤變化遙感定量反演方法研究——以長江口崇明東灘為例[J].海洋湖沼通報,2009,11(1):12-18.

HAN Z,YUN C X,DAIZ J,et al.Remote sensing quantitative inversion research of the silt tidal flat elevation and deposition and erosion—Chongming Dongtan as the example[J].Transactions of Oceanology and Limnology,2009,11(1):12-18.

[30] OUYANG Z T,GAO Y,XIE X,et al.Spectral discrimination of the invasive plant Spartina alterniflora at multiple phenological stages in a saltmarsh wetland[J].Plos One,2013,8(6):e67315.

[31] ZHANG M Q,GUO H Q,XIE X,et al.Identification of land-cover characteristics using MODIS time series data:An application in the Yangtze River Estuary[J].Plos One,2013,8(7):e70079.

[32] 郭永飛.長江口水邊線和潮溝遙感信息提取研究[D].上海:上海海洋大學,2011.

GUO Y F.The study on waterside line and tidal channel remote sensing information extraction in the Yangzte River Estuary[D].Shanghai:Shanghai Ocean University,2011.

[33] 彭 飛.長江口南匯潮灘土壤濕度遙感研究[D].上海:上海海洋大學,2012.

PENG F.Soil moiture quantitative study of the Nanhui tidal flat in the Yangzte River Esturay[D].Shanghai:ShanghaiOcean University,2012.

[34] FREITASR,RODRIGUESA M and QUINTINO V.Benthic biotopes remote sensing using acoustics[J].Journal of Experimental Marine Biology and Ecology,2003,285:339-353.

[35] ZIV Y.The effect of habitatheterogeneity on species diversity patterns:a community-level approach using an object-oriented landscape simulation model(SHALOM)[J].Ecological Modelling,1998,111(2-3):135-170.

[36] 丁 駿.長江口中華鱘棲息地景觀格局變化及其適宜性評價[D].上海:上海海洋大學.2012.

DING J.The Chinese sturgeon habitat landscapepattern change in the Yangtze Estuary and its habitat suitability assessment[D].Shanghai:Shanghai Ocean University,2012.

[37] 平 原.南匯邊灘濕地景觀格局分析和生態系統健康評價[D].上海:華東師范大學,2010.

PING Y.The landscape pattern analysis and assessment of wetland ecosystem health for Nanhui Biantan,Shanghai[D].Shanghai:East China Normal University,2010.

[38] 孫永光.長江口不同年限圍墾區景觀結構與功能分異[D].上海:華東師范大學,2011.

SUN Y G.Differentiation of landscape structure and function during different reclamation time at the Changjiang River Estuary[D].Shanghai:East China Normal University,2011.

[39] PIERRI-DAUNT A B,TANAKA M O.Assessing habitat fragmentation on marine epifaunal macroinvertebrate communities:an experimental approach[J].Landscape Ecology,2014,29(1):17-28.

[40] JUAN SD,HEWITT J.Relative importance of local biotic and environmental factors versus regional factors in driving macrobenthic species richness in intertidal areas[J].Marine Ecology Progress Series,2011,423(1):117-129.

[41] BEATRIZ P,TANAKA M O.Assessing habitat fragmentation on marine epifaunal macroinvertebrate communities:an experimental approach[J].Landscape Ecology,2014,29(1):17-28.

[42] THRUSH SF,HALLIDAY J,HEWITT JE,et al.The effects of habitat loss,fragmentation,and community homogenization on resilience in estuaries[J].Ecological Applications,2008,18(1):12-21.

[43] ELITH J,LEATHWICK J R.Species Distribution Models:Ecological Explanation and Prediction Across Space and Time[J].Annual Review of Ecology Evolution and Systematics,2009,40(40):677-697.

[44] WANG K,FRANKLIN S E,GUO X L,et al.Problems in remote sensing of landscapes and habitats[J].Progress in Physical Geography,2009,33(6):747-768.

[45] LECOMTE JB,BENOIT H P,ETIENNE M P,etal.Modeling the habitat associations and spatial distribution of benthic macroinvertebrates:A hierarchical Bayesianmodel for zero-inflated biomass data[J].Ecological Modelling,2013,265(1767):74-84.

[46] MARTIN TG,WINTLE B A,RHODES JR,et al.Zero tolerance ecology:improving ecological inference by modelling the source of zero observations[J].Ecology Letters,2005,8(11):1235-1246.

[47] MAUNDER M N,PUNT A E.Standardizing catch and effort data:a review of recent approaches[J].Fisheries Research,2004,70(2-3):141-159.

[48] STEFANSSON G.Analysis of groundfish survey abundance data:Combining the GLM and delta approaches[J].Ices Journal of Marine Science,1996,53(3):577-588.

[49] BUYANTUYEV A,WU J G.Effects of thematic resolution on landscape pattern analysis[J].Landscape Ecology,2007,22(1):7-13.

[50] HUANG C,GEIGER E L,KUPFER JA.Sensitivity of landscape metrics to classification scheme[J].International Journal of Remote Sensing,2006,27(14):2927-2948.

[51] CASTILLA G,LARKIN K,LINKE J,et al.The impact of thematic resolution on the patch-mosaic model of natural landscapes[J].Landscape Ecology,2009,24(1):15-23.

[52] MACEDO D R,HUGHESR M,LIGEIRO R,et al.The relative influence of catchment and site variables on fish and macroinvertebrate richness in cerrado biome streams[J].Landscape Ecology,2014,29(6):1001-1016.

[53] KOENIG W D.Spatial autocorrelation of ecological phenomena[J].Trends in Ecology&Evolution,1999,14(1):22-26.

[54] LICHSTEIN JW,SIMONSTR,SHRINER SA,et al.Spatial autocorrelation and autoregressivemodels in ecology[J].Ecological Monographs,2002,72(3):445-463.

[55] KISSLING W D,CARL G.Spatial autocorrelation and the selection of simultaneous autoregressive models[J].Global Ecology and Biogeography,2008,17(1):59-71.

[56] ELITH J,GRAHAM C H,ANDERSON R P,et al.Novel methods improve prediction of species'distributions from occurrence data[J].Ecography,2006,29(2):129-151.

[57] GALPARSORO I,BORJA A,BALD J,et al.Predicting suitable habitat for the European lobster(Homarus gammarus),on the Basque continental shelf(Bay of Biscay),using Ecological-NicheFactor Analysis[J].Ecological Modelling,2009,220(4):556-567.

[58] IAMPIETRO P J,YOUNG M A,KVITEK R G.Multivariate Prediction of Rockfish Habitat Suitability in Cordell Bank National Marine Sanctuary and Del Monte Shalebeds,California,USA[J].Marine Geodesy,2008,31(4):359-371.

[59] CONNOR DW,ALLEN JH,GOLDING N,et al.The Marine Habitat Classification for Britain and Ireland Version 04.05.[Z].2004.www.jncc.gov.uk/MarineHabitatClassification(internet version).

[60] NODA T.Spatialhierarchicalapproach in community ecology:a way beyond high context-dependency and low predictability in local phenomena[J].Population Ecology,2004,46(2):105-117.

[61] THRUSH S F,HEWITT JE,LOHRER A M,et al.When small changes matter:the role of cross-scale interactions between habitat and ecological connectivity in recovery[J].Ecological Applications,2013,23(1):226-238.

[62] POTTER K M,CUBBAGE FW,SCHABERG R H.Multiple-scale landscape predictors of benthic macroinvertebrate community structure in North Carolina[J].Landscape and Urban Planning,2005,71(2-4):77-90.

[63] YEUNG C,MCCONNAUGHEY R A.Using acoustic backscatter from a sidescan sonar to explain fish and invertebrate distributions:a case study in Bristol Bay,Alaska[J].Ices Journal of Marine Science,2008,65(2):242-254.

[64] CALLAWAY A,SMYTH J,BROWN C J,et al.The impact of scour processes on a smothered reef system in the Irish Sea[J].Estuarine Coastal and Shelf Science,2009,84(3):409-418.

A review of application of spatial technology in mapping intertidal habitat in estuary

ZHANG Ting-ting,ZHAO Feng,HOU Jun-li,GAO Yu,ZHUANG Ping
(Key Laboratory of East China Sea&Oceanic Fishery Resources Exploitation and Utilization,East China Sea fisheries Research Institute,Chinese Academy of Fishery Sciences,Shanghai200090,China)

Mapping the changes of intertidal benthic habitat(IBH)in estuary provides important scientific proofs and management suggestions formany areas,e.g.,the biological conservation of wetland in estuary,the assessment of benthic habitat quality and ecological recovery of damaged habitat.With the development of spatial information technologies(SIT),scientists have the ability to acquire more explicit,holistic and systematic monitoring and assessment for intertidal benthic habitat in estuary.From the spatial context,this article summarized six aspects of progress in mapping intertidal habitat in estuary recently:(1)new spatialtemporal characteristics of IBH definition in esturary;(2)description of characteristics of IBH in estuary;(3)acquisition methods of IBH characteristics data in estuary;(4)data processing;(5)mapping IBH with differentmodels based on species/community;(6)multiple temporal-spatial scales analysis.Based on these progresses,the IBH in estuary can be fully,systematically and explicitly characterized.By consideration of issues in spatial ecology(boundary segmentation,scales effects,classification accuracy,dynamics change,etc.),the distributionmodelling of benthic species/community and habitatmapping becomemore convenient with more spatial explicitly,realistic and holistic results.Finally,this paper proposed that mapping with spatial information technologieswould bring new prosepects for studies of ecosystem evaluation and ecological function in intertial habitats.Nevertheless,there are some issues that should be carefully deltwith inmapping with spatial information technologies.For example,characteristics of intertidal benthic habitat in estuary should be seriously chosen and identified as the most effective indicators with proper scientific statisitics at different scales.Presently,the spatial information technologies have not been widely applied in mapping intertidal habitat in estuary,which needs urgent development and promotion.

spatial information technologies;fisheries resources;ecological monitoring;ecological assessment;application prospect

S 931

A

1004-2490(2016)04-0424-13

2015-05-26

國家自然科學基金(31400410)

張婷婷(1985-),女,助理研究員,博士,主要從事河口生態監測和評價科研工作。E-mail:zhangtt@ecsf.ac.cn

猜你喜歡
物種分類特征
吃光入侵物種真的是解決之道嗎?
英語世界(2023年10期)2023-11-17 09:18:18
分類算一算
如何表達“特征”
不忠誠的四個特征
當代陜西(2019年10期)2019-06-03 10:12:04
分類討論求坐標
回首2018,這些新物種值得關注
電咖再造新物種
汽車觀察(2018年10期)2018-11-06 07:05:26
數據分析中的分類討論
教你一招:數的分類
抓住特征巧觀察
主站蜘蛛池模板: 精品成人一区二区| 精品人妻一区无码视频| 国产91特黄特色A级毛片| 超碰aⅴ人人做人人爽欧美| 国产亚洲欧美在线人成aaaa| 麻豆国产精品视频| 成人另类稀缺在线观看| 伊人久久综在合线亚洲2019| jizz亚洲高清在线观看| 亚洲网综合| 亚洲欧洲一区二区三区| 97精品久久久大香线焦| 男女男精品视频| 欧美日韩国产综合视频在线观看| 在线综合亚洲欧美网站| 国产亚洲男人的天堂在线观看| 国产美女一级毛片| 成人亚洲天堂| 伊人无码视屏| 免费观看国产小粉嫩喷水| 日韩国产高清无码| 亚洲最大在线观看| 欧美一级色视频| 四虎成人免费毛片| 国产成人夜色91| 国产精品丝袜在线| 国内精品久久九九国产精品| 色综合天天娱乐综合网| 久久精品国产91久久综合麻豆自制| 国产成年女人特黄特色毛片免 | 欧美特黄一级大黄录像| 五月婷婷激情四射| 免费一级毛片| 高清码无在线看| 四虎影视8848永久精品| 天天色天天操综合网| 亚洲人人视频| 国产青青操| 91 九色视频丝袜| 宅男噜噜噜66国产在线观看| 一级高清毛片免费a级高清毛片| 蜜芽国产尤物av尤物在线看| 福利一区在线| 白浆视频在线观看| 麻豆AV网站免费进入| 国产大片黄在线观看| 欧美成一级| 狠狠ⅴ日韩v欧美v天堂| 国产91麻豆免费观看| 欧美成人精品一区二区| 欧美三级自拍| 午夜精品一区二区蜜桃| 91成人精品视频| 欧美19综合中文字幕| 欧美日本一区二区三区免费| 激情在线网| 国产精品久久久久久影院| 综合网久久| 亚洲成人手机在线| 黄色网在线| 91精品人妻一区二区| 三上悠亚精品二区在线观看| 日韩高清成人| 18禁色诱爆乳网站| 久视频免费精品6| 美女内射视频WWW网站午夜 | 亚洲精品卡2卡3卡4卡5卡区| 欧美午夜精品| 亚洲 欧美 中文 AⅤ在线视频| 国产婬乱a一级毛片多女| 国产手机在线ΑⅤ片无码观看| 2020国产精品视频| 视频一区亚洲| 国产在线精品99一区不卡| 日本在线免费网站| 国产午夜无码片在线观看网站| 亚洲高清中文字幕在线看不卡| 国产99在线观看| 亚洲男人天堂2020| 青青青视频91在线 | 色亚洲成人| 91精品最新国内在线播放|