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

    海南東寨港紅樹林急速退化的空間分布特征及影響因素分析

    2016-01-25 09:17:01孫艷偉廖寶文王榮麗丁冬靜鐘才榮中國林業(yè)科學(xué)研究院熱帶林業(yè)研究所廣東廣州5050海南東寨港國家級(jí)自然保護(hù)區(qū)管理局海南???79
    關(guān)鍵詞:空間分布紅樹林影響因素

    孫艷偉,廖寶文,管 偉,王榮麗,丁冬靜,鐘才榮(中國林業(yè)科學(xué)研究院熱帶林業(yè)研究所,廣東廣州5050; 海南東寨港國家級(jí)自然保護(hù)區(qū)管理局,海南???79)

    海南東寨港紅樹林急速退化的空間分布特征及影響因素分析

    孫艷偉1,廖寶文1,管偉1,王榮麗1,丁冬靜1,鐘才榮2
    (1中國林業(yè)科學(xué)研究院熱帶林業(yè)研究所,廣東廣州510520; 2海南東寨港國家級(jí)自然保護(hù)區(qū)管理局,海南海口571129)

    摘要:【目的】對(duì)我國一些紅樹林濕地生態(tài)系統(tǒng)近期出現(xiàn)的急速退化現(xiàn)象作相關(guān)分析.【方法】以海南東寨港整個(gè)港灣紅樹林區(qū)為研究典型區(qū)域,應(yīng)用3S技術(shù),分析近30年間東寨港紅樹林濕地景觀格局的演變過程;同時(shí)采用高分辨率影像數(shù)據(jù),聚焦分析了演豐東河沿岸、塔市及三江3個(gè)片區(qū)紅樹林群落退化特征及空間分布規(guī)律;結(jié)合大量野外樣方調(diào)查數(shù)據(jù),系統(tǒng)分析了近海岸紅樹林生態(tài)系統(tǒng)退化的主要影響因素.【結(jié)果和結(jié)論】從1999—2013年,東寨港紅樹林面積由1 709. 4 hm2減少至1 679. 5 hm2,年均減少2. 1 hm2.而沿海岸線2 km緩沖區(qū)范圍內(nèi),周邊養(yǎng)殖塘面積則由1987年的59. 1 hm2增至2013年的1 986. 9 hm2.截止2013年,東寨港紅樹林枯死面積約為4 hm2.從空間分布上來看,紅樹林枯死群落主要分布于地形低洼積水處或污染物傳輸通道兩側(cè).東寨港河流型紅樹林群落退化最為嚴(yán)重,集中分布于演豐東河中下游河段;同時(shí),近海灣的前沿地段也有退化群落零星分布.處于地帶性演替后期的紅樹林群落,退化程度較高.陸源污染物輸入強(qiáng)度是影響紅樹林退化群落分布的重要因素.

    關(guān)鍵詞:東寨港;紅樹林;退化特征;空間分布;影響因素

    孫艷偉,廖寶文,管偉,等.海南東寨港紅樹林急速退化的空間分布特征及影響因素分析[J].華南農(nóng)業(yè)大學(xué)學(xué)報(bào),2015,36(6):111-118.

    優(yōu)先出版時(shí)間:2015-10-16

    優(yōu)先出版網(wǎng)址: http: / /www.cnki.net/kcms/detail/44.1110.s.20151016.1631.034.html

    長期以來,地處熱帶和亞熱帶海岸潮間帶的紅樹林群落正經(jīng)歷著嚴(yán)重的退化過程[1].在發(fā)展中國家,紅樹林的退化速度甚至超過很多其他熱帶樹種[2].人類活動(dòng)對(duì)紅樹林生態(tài)系統(tǒng)服務(wù)功能的影響可分為積極影響和消極影響.積極影響主要有生態(tài)系統(tǒng)管理、生態(tài)工程、生態(tài)恢復(fù)與重建、生態(tài)評(píng)價(jià)與規(guī)劃等;消極影響包括圍墾、砍伐、采礦、挖捕動(dòng)物、海水養(yǎng)殖、污水排放、旅游等[3].我國沿海區(qū)域海岸帶占國土面積13%,承載著約40%的全國人口,創(chuàng)造了約60%的GDP[4].沿海區(qū)域經(jīng)濟(jì)的高速發(fā)展和快速的城市化建設(shè)直接吞噬和影響著濱海紅樹林濕地生態(tài)系統(tǒng),土地利用的變化和陸源污染物的大量輸入使紅樹林面積不斷縮減、生態(tài)功能急劇衰退[5].因此,研究濱海紅樹林生態(tài)系統(tǒng)的退化過程及其機(jī)制對(duì)紅樹林濕地資源的保護(hù)和恢復(fù)具有重要的科學(xué)意義和現(xiàn)實(shí)意義.

    海南東寨港是我國最典型、最原始的天然紅樹林分布區(qū)[6].近年來,東寨港紅樹林區(qū)出現(xiàn)一定范圍的海洋污損動(dòng)物——團(tuán)水虱(Sphaeromadae)危害,已造成了紅樹林死亡現(xiàn)象,引起了政府和社會(huì)的高度關(guān)注[7].影響團(tuán)水虱分布的因素復(fù)雜多樣,既包括紅樹林的群落特征、立地條件(地形因子、水質(zhì)因子),也包括陸源污染物的輸入等人類活動(dòng)的擾動(dòng).而目前學(xué)者對(duì)影響團(tuán)水虱地理分布的相關(guān)環(huán)境因子尚未形成確定的結(jié)論.Brook[8]認(rèn)為佛羅里達(dá)灣北部紅樹林中團(tuán)水虱的分布與溫度、鹽度、pH等水質(zhì)因素沒有顯著性關(guān)系,并認(rèn)為其他鉆孔動(dòng)物對(duì)團(tuán)水虱的分布起一定限制作用.邱勇等[9]認(rèn)為東寨港紅樹林區(qū)水體中總氮、總磷含量和浮游生物量是影響光背團(tuán)水虱分布的主要水質(zhì)因子,并且網(wǎng)紋藤壺是光背團(tuán)水虱在垂直分布上的限制因素.本文借助5期Landsat衛(wèi)星影像及近期高分辨率航拍數(shù)據(jù),結(jié)合實(shí)地踏勘調(diào)查,對(duì)整個(gè)東寨港灣紅樹林群落的退化特征及其空間分布規(guī)律進(jìn)行了系統(tǒng)的調(diào)查研究,同時(shí)還進(jìn)一步深入分析了地形以及陸源污染物的輸入強(qiáng)度等環(huán)境因子對(duì)紅樹林群落退化的影響.

    1 材料與方法

    1.1研究區(qū)概況

    海南東寨港紅樹林區(qū)(指東寨港整個(gè)港灣區(qū)域)位于海口市美蘭區(qū)東北部(110°30'~110°37' E,19°51'~20°01' N),目前擁有紅樹林總面積約為1 679. 5 hm2.該區(qū)域?qū)俚湫偷臒釒Ъ撅L(fēng)海洋性氣候,年平均氣溫23. 3~23. 8℃,極端最高氣溫38. 9℃,極端最低氣溫2. 6℃;年均日照2 000 h以上;年均相對(duì)濕度85%;年均降水量1 676. 4 mm;潮型為不規(guī)則的半日潮[6].東寨港紅樹林樹種主要有:白骨壤Avicennia marina、桐花樹Aegiceras cornicu latum、木欖Bruguiera gymnorrhiza、秋茄Kandelia obovata、角果木Ceriops tagal、海蓮B.sexangula、海漆Excoecaria agallocha、欖李Lumnitzera racemosa、紅海欖Rhizophorasty losa、老鼠簕Acanthus ilicifolius等.

    近半個(gè)世紀(jì)以來,在經(jīng)濟(jì)利益的驅(qū)使下,整個(gè)東寨港圍墾、毀林建塘、過度旅游開發(fā)等現(xiàn)象不斷涌現(xiàn),致使紅樹林群落呈現(xiàn)不斷退化的趨勢.據(jù)不完全統(tǒng)計(jì),20世紀(jì)60年代至今,由于圍墾,海南東寨港紅樹林面積減少近50%[10].近幾年,團(tuán)水虱爆發(fā)是導(dǎo)致紅樹林死亡退化的直接原因.盡管保護(hù)區(qū)管理局采取了人工種植、關(guān)閉養(yǎng)鴨塘、禁止圍墾和新建蝦塘等多項(xiàng)措施,并取得了一定成效,但紅樹林退化現(xiàn)象并沒有得到有效控制.

    1.2研究方法

    1.2.1紅樹林植被提取本研究所選用的數(shù)據(jù)主要來源于5期Landsat衛(wèi)星影像,主要包括LT51230461987253BJC00、elp123r046 _ 7t19991224、LT51230462004044BJC00、LT51230462009281BJC00、LC81230462013180LGN00;結(jié)合研究區(qū)1∶10 000地形圖,在ENVI5. 0平臺(tái)下,采用二次多項(xiàng)式糾正法對(duì)覆蓋研究區(qū)的5期遙感影像進(jìn)行幾何校正(誤差不超過0. 5個(gè)像元)、圖像增強(qiáng)和區(qū)域剪切等處理.根據(jù)野外現(xiàn)場調(diào)查及專家經(jīng)驗(yàn),采用目視解譯法,提取研究區(qū)海陸邊界2 km緩沖區(qū)范圍以內(nèi)的水產(chǎn)養(yǎng)殖塘區(qū).實(shí)地踏勘調(diào)查發(fā)現(xiàn),東寨港紅樹林主要分布于距海岸線1 km緩沖區(qū)范圍內(nèi).基于此,本文參考Liu等[11]提出的基于輔助數(shù)據(jù)的專家決策樹紅樹林分類思路,采用分類回歸樹分析(CART),整合研究區(qū)DEM、海陸分界線以及NDVI等輔助數(shù)據(jù),分別構(gòu)建了5期東寨港紅樹林分類的專家決策樹,結(jié)果顯示紅樹林分類精度達(dá)到90%.

    1.2.2退化紅樹林群落提取為準(zhǔn)確提取枯死紅樹林群落的分布及面積,采用2012年2月24日SPOT5遙感影像及近期高分航拍影像為數(shù)據(jù)源,經(jīng)幾何校正、區(qū)域裁切后,結(jié)合專家目視解譯,獲取東寨港枯死紅樹林群落的空間分布及面積.

    1.2.3紅樹林退化樣方調(diào)查考慮到紅樹林團(tuán)水虱危害發(fā)生的立地條件及群落結(jié)構(gòu)的差異性,分別在東寨港紅樹林區(qū)的塔市、演豐東河及三江等3個(gè)片區(qū)典型區(qū)域共設(shè)置28個(gè)固定樣方,樣方尺寸為10 m×15 m.于2013年7月和10月對(duì)整個(gè)東寨港紅樹林區(qū)的退化情況進(jìn)行了現(xiàn)場踏勘調(diào)查,調(diào)查記錄每個(gè)樣方的直角坐標(biāo)系、不同樹種個(gè)體數(shù)、郁閉度及樹均高,同時(shí)記錄樣方對(duì)角線上每株樹木的種類、胸徑、冠幅、蛀蝕高度以及蛀蝕等級(jí)等參數(shù),并計(jì)算每個(gè)樣方的死亡率和退化率.調(diào)查樣方的布設(shè)情況詳見圖1.依據(jù)樹木主干、氣生根的蛀蝕程度以及枝葉凋落程度將樹種退化程度劃分為6個(gè)等級(jí),退化等級(jí)由低到高依次為健康、1級(jí)、2級(jí)、3級(jí)、4級(jí)和5級(jí),退化等級(jí)越高表明樹木的蛀蝕程度越嚴(yán)重.

    圖1 東寨港紅樹林群落調(diào)查樣方位置示意圖Fig.1 The location map of sampling plots of mangrove forest in Dongzhaigang Harbor

    1.2.4相關(guān)環(huán)境因子的獲取選取高程、坡度、坡向等地形因子以及陸源污染物輸入強(qiáng)度作為紅樹林退化的影響因子.其中,地形因子數(shù)值應(yīng)用ArcGIS軟件的空間分析功能提取每個(gè)樣地所處空間網(wǎng)格的柵格值.而陸源污染物輸入強(qiáng)度主要考慮臨近海岸線2 km緩沖區(qū)范圍內(nèi)的陸源總N(TN)、總P(TP)等營養(yǎng)物質(zhì)的年輸出負(fù)荷量.首先,應(yīng)用ArcGIS軟件將研究區(qū)劃分為2 km×2 km網(wǎng)格,然后將2013年土地利用圖層與網(wǎng)格進(jìn)行空間關(guān)聯(lián),統(tǒng)計(jì)每個(gè)網(wǎng)格內(nèi)各土地利用類型的面積,最后采用經(jīng)典的Johns輸出系數(shù)法估算每個(gè)網(wǎng)格的陸源面源污染負(fù)荷量,輸出系數(shù)參見表1.

    表1 研究區(qū)N、P輸出系數(shù)1)Tab.1  Nitrogen and phosphorus export coefficients in the study area t·km-2

    演豐東河為內(nèi)陸淡水河流,自西向東匯入東寨港.為分析陸源污染負(fù)荷輸入對(duì)河流型紅樹林退化的影響,本文采用SWAT模型將演豐東河流域劃分為15個(gè)子流域單元,同時(shí)利用Johns輸出系數(shù)法估算每個(gè)子流域陸源面源污染負(fù)荷輸出量,主要考慮TN、TP等營養(yǎng)物質(zhì)負(fù)荷量.

    基于樣方調(diào)查及相關(guān)環(huán)境因子提取,應(yīng)用Spearman相關(guān)性分析方法,定量分析紅樹林群落退化程度(樣方的死亡率、退化率)與立地條件(包括Shannon-Wiener指數(shù)[14]、地表高程、坡度及坡向等)、陸源污染物輸入強(qiáng)度(TN、TP等陸源營養(yǎng)物質(zhì)輸入負(fù)荷量等)等環(huán)境因子的關(guān)聯(lián)性.

    2 結(jié)果與分析

    2.1近30年紅樹林面積動(dòng)態(tài)變化

    利用1987、1999、2004、2009和2013年等5個(gè)時(shí)相的Landsat衛(wèi)星影像,提取得到東寨港近30年來紅樹林空間分布特征,面積動(dòng)態(tài)變化情況見表2和圖2.從統(tǒng)計(jì)結(jié)果來看,近30年來東寨港紅樹林面積呈現(xiàn)先增加后縮減的趨勢,同時(shí)斑塊破碎化不斷加劇.自1986年東寨港升級(jí)為國家級(jí)紅樹林自然保護(hù)區(qū)以后,隨著管理力度的加強(qiáng),紅樹林得到保護(hù)和恢復(fù).至1999年為止,整個(gè)港灣紅樹林面積達(dá)到1 709. 4 hm2,與1987年相比面積增加了11. 18%.但自此以后,在經(jīng)濟(jì)利益的驅(qū)使下,港灣內(nèi)圍墾、毀林建塘、過度旅游開發(fā)等現(xiàn)象不斷涌現(xiàn),紅樹林灘涂區(qū)域不斷被擠占和破壞,致使紅樹林面積不斷縮減.自2009年以來,局部區(qū)域團(tuán)水虱危害嚴(yán)重,進(jìn)一步加劇了紅樹林群落的退化.截止2013年,東寨港紅樹林面積約為1679. 5hm2,與1999年相比減少了29. 9 hm2,年均減少2. 1 hm2.另外,遙感影像判讀結(jié)果顯示,沿海岸線2 km緩沖區(qū)范圍內(nèi),周邊養(yǎng)殖塘面積由1987年的59. 1 hm2增至2013年的1 986. 9 hm2,增加了33倍.

    由圖2可見,紅樹林斑塊面積的縮減主要表現(xiàn)為斑塊邊緣萎縮和潮溝切割所致的斑塊內(nèi)部退化.在海岸帶開發(fā)過程中,對(duì)河流-潮汐通道截彎取直以及旅游船舶引起的頻繁巨浪都造成紅樹植物的死亡和退化[15].岸邊畜禽養(yǎng)殖、陸源污染物的排放以及東寨港的旅游開發(fā)等人類開發(fā)活動(dòng)改變了紅樹林生態(tài)系統(tǒng)的生境、結(jié)構(gòu)和生物地化循環(huán),從而引起了生態(tài)系統(tǒng)服務(wù)功能的降低.潮溝是海陸物質(zhì)交換的主要通道[16].自1987年以來,紅樹林斑塊內(nèi)部的潮溝呈增多的趨勢,直接導(dǎo)致同一斑塊割裂為幾個(gè),破碎化加劇.

    表2 1987—2013年東寨港紅樹林面積變化情況Tab.2 Changes of mangrove area in Dongzhaigang Harbor from 1987 to 2013

    圖2 1987—2013年東寨港紅樹林分布動(dòng)態(tài)變化Fig.2 Dynamic changes of mangrove distribution in Dongzhaigang Harbor from 1987 to 2013

    2.2紅樹林群落退化特征分析

    2.2.1紅樹林退化群落的空間分布特征自2009年以來,團(tuán)水虱大量爆發(fā),導(dǎo)致部分區(qū)域紅樹林大面積枯死,紅樹林濕地生態(tài)系統(tǒng)退化嚴(yán)重.東寨港紅樹林可分為河流型紅樹林群落(演豐東河片區(qū))和前沿型紅樹林群落(塔市、三江片區(qū)) 2個(gè)類型.經(jīng)過全面實(shí)地踏勘調(diào)查及近期高分遙感影像分析顯示,東寨港紅樹林區(qū)范圍內(nèi)團(tuán)水虱危害地理空間分布廣泛.從東寨港河流型紅樹林分布來看,河流中下游紅樹林群落退化嚴(yán)重,同時(shí),近海灣的前沿地段退化群落有零星分布.由團(tuán)水虱蛀蝕形成大小不等的林窗斑塊共27塊,總枯死面積約為4 hm2,年均枯死面積0. 8 hm2(圖3).從空間分布來看,塔市紅樹林片區(qū)出現(xiàn)4個(gè)斑塊(0. 738 hm2),分布于靠近陸地一側(cè),主要危害樹種為白骨壤;演豐東河紅樹林片區(qū)出現(xiàn)21個(gè)斑塊(3. 17 hm2),分布于河流的中下游區(qū)域,主要危害樹種為木欖、海蓮、尖瓣海蓮等處于演替后期的樹種;三江紅樹林片區(qū)出現(xiàn)2個(gè)斑塊(0. 087 hm2),分布于潮灘前沿,主要危害樹種為秋茄.以上結(jié)果說明,對(duì)于河流型紅樹林群落,處于地帶性演替后期的樹種是團(tuán)水虱危害的主要對(duì)象.

    圖3 2013年東寨港紅樹林枯死斑塊分布Fig.3 The distribution map of death patches of mangroves in Dongzhaigang Harbor in 2013

    2.2.2不同樹種的退化特征典型樣方調(diào)查顯示,所有調(diào)查樣方的退化率約為16%,但河流型紅樹林群落與前沿型紅樹林群落不同樹種的退化程度呈現(xiàn)出差異性(表3).對(duì)于河流型紅樹林群落,演豐東河片區(qū)不同樹種的退化程度由重到輕依次為:木欖(28. 8%)>尖瓣海蓮(10. 7%)>海蓮(10. 4%)>秋茄>紅海欖;對(duì)于前沿型紅樹林群落,塔市片區(qū)白骨壤群落和三江片區(qū)秋茄的退化率都為15. 0%.分析表明,團(tuán)水虱對(duì)紅樹林不同樹種的蛀蝕具有選擇性.

    表3 各樣地的群落退化特征1)Tab.3 Community degradation characteristics of sample plots

    2.2.3群落特征與退化特征的相關(guān)性紅樹林群落特征在一定程度上反映了紅樹林的生長狀況,是紅樹林群落退化的直觀反映.紅樹林群落的退化特征可采用死亡率和退化率2個(gè)指標(biāo)來表達(dá).典型樣方的Spearman相關(guān)性分析顯示,樣方的死亡率與種群密度和群落多樣性指數(shù)(Shannon-Wiener指數(shù))呈顯著正相關(guān),而與郁閉度、平均胸徑及平均冠幅相關(guān)性不顯著;樣方的退化率則與郁閉度呈顯著負(fù)相關(guān),而與其他指標(biāo)相關(guān)性不顯著(表4).以上結(jié)果表明,處于地帶性演替后期的紅樹林群落更易遭受團(tuán)水虱的蛀蝕而死亡,而郁閉度是群落退化程度的直觀反映.

    表4 各樣地的群落特征與退化特征相關(guān)性分析1)Tab.4 The correlation analyses of community characteristics and degradation characteristics of sample plots

    2.3紅樹林退化程度與環(huán)境因子的關(guān)系

    海灣是人類活動(dòng)最為密集的區(qū)域之一.隨著人類對(duì)海灣開發(fā)進(jìn)程的加快,人類活動(dòng)對(duì)紅樹林濕地生態(tài)系統(tǒng)的擾動(dòng)也日益加劇.紅樹植物的生長不僅受到灘涂高程、地形以及淹水時(shí)間等生境立地條件的影響,同時(shí)也受到人類活動(dòng)產(chǎn)生的陸源污染物輸入的影響.

    2.3.1立地條件與退化程度的關(guān)系應(yīng)用GIS軟件的空間分析功能,采用DEM數(shù)據(jù)分別提取了典型樣地所處地理網(wǎng)格的高程、坡度以及坡向等地形因子的數(shù)值.相關(guān)性分析結(jié)果表明,樣地的死亡率與地表高程呈顯著負(fù)相關(guān)關(guān)系,而與坡度及坡向相關(guān)關(guān)系不顯著(表5).在現(xiàn)場踏勘調(diào)查過程也發(fā)現(xiàn),地形低洼處、小排污口以及潮溝兩側(cè)紅樹被團(tuán)水虱蛀蝕的程度更為嚴(yán)重.由此可見,污染物傳輸通道的邊緣以及地形低洼的沉積區(qū)更易受到團(tuán)水虱的危害,這與范航清等[16]的研究結(jié)論一致.

    灘涂表面的相對(duì)高程是影響紅樹植物分布的關(guān)鍵生境因子,直接導(dǎo)致紅樹植物淹沒深度和水淹時(shí)間的差異.演豐東河片區(qū)紅樹林群落以木欖、海蓮等演替后期樹種為優(yōu)勢種,而這些樹種更適宜在較高潮灘的生境生長[17].放養(yǎng)鴨群對(duì)表層土壤的松動(dòng)、高位養(yǎng)殖塘排污口排放污水的沖刷等都會(huì)導(dǎo)致局部地形的降低,退潮時(shí)凹陷的地表形成積水,致使紅樹林水淹時(shí)間變長,長時(shí)間的淹水不僅超出了紅樹植物的忍耐程度,而且也給團(tuán)水虱提供更多的濾食時(shí)間,導(dǎo)致團(tuán)水虱數(shù)量的增長,從而更快更多地破壞紅樹植物的呼吸根及根莖[7].隨著根系不斷被蛀蝕和腐爛,地表又會(huì)進(jìn)一步發(fā)生沉降,最終導(dǎo)致紅樹林退化程度不斷加重.此外,典型樣地調(diào)查結(jié)果還表明,紅樹植物的枯亡與所處地理位置的坡度和坡向的關(guān)聯(lián)不大.

    表5 各樣地的地形因子與退化特征相關(guān)性分析1)Tab.5  The correlation analyses of topographical factors and community degradation characteristics of sample plots

    2.3.2陸源污染物輸入強(qiáng)度與退化程度的關(guān)系陸源污染物的輸入已經(jīng)成為近岸海域污染的主要來源.其中,總N、總P等陸源營養(yǎng)物質(zhì)的大量輸入是導(dǎo)致海水富營養(yǎng)化的主要原因,水體中總氮、總磷含量和浮游生物量也是影響光背團(tuán)水虱Sphraeroma retrolaevis分布的主要水質(zhì)因子[9].對(duì)于東寨港地區(qū),陸源營養(yǎng)物質(zhì)一部分來源于海灣周邊大面積的耕地、林地等產(chǎn)生的面源污染,而另一個(gè)主要來源是養(yǎng)殖塘換水和滲漏輸出的污染物,并占到營養(yǎng)物質(zhì)輸出總量的四成以上.近年來東寨港周邊養(yǎng)殖塘呈現(xiàn)快速增加的趨勢,大量營養(yǎng)物質(zhì)的輸出正不斷污染臨近海域.水體的富營養(yǎng)化也進(jìn)一步加劇了沿岸紅樹林的退化.

    為分析陸源污染物的輸入對(duì)紅樹林退化的影響,本文將陸源TN、TP輸入強(qiáng)度圖層與紅樹林枯死斑塊圖層進(jìn)行疊加,可以直觀表達(dá)陸源污染物的輸入強(qiáng)度與紅樹林退化斑塊之間的空間匹配關(guān)系.如圖4所示,從空間位置關(guān)系上來看,紅樹林枯死斑塊與陸源TN、TP輸入強(qiáng)度高值區(qū)具有很好的對(duì)應(yīng)關(guān)系.整體上來看,紅樹林枯死斑塊所對(duì)應(yīng)的陸源TN輸入強(qiáng)度介于1. 1~3. 6 t·km-2之間,高于此區(qū)域TN平均輸入強(qiáng)度0. 9 t·km-2;而陸源TP輸入強(qiáng)度介于0. 24~0. 35 t·km-2之間,高于此區(qū)域TN平均輸入強(qiáng)度0. 2 t·km-2.

    圖5為演豐東河流域陸源污染物輸入強(qiáng)度與紅樹林枯死斑塊的空間關(guān)系.由圖5可知,紅樹林枯死斑塊主要分布于演豐東河中下游區(qū)域,與陸源TN、 TP輸入強(qiáng)度高值區(qū)具有較好的一致性.上述結(jié)果表明,周邊陸源營養(yǎng)物質(zhì)的輸入強(qiáng)度是影響紅樹林退化群落分布的重要因素.

    圖4 東寨港陸源污染物輸入強(qiáng)度與紅樹林枯死斑塊的空間關(guān)系Fig.4 The spatial corresponding relationship between terrestrial pollutant loading and death patch in Dongzhaigang Harbor

    圖5 演豐東河流域陸源污染物輸入強(qiáng)度與紅樹林枯死斑塊的空間關(guān)系Fig.5  The spatial corresponding relationship between terrestrial pollutant loading and death patch in the Yanfengdong River watershed

    3 討論

    近海岸帶的紅樹林生態(tài)系統(tǒng)具有極其重要的生態(tài)和環(huán)境價(jià)值.但大量研究表明氣候變化、病蟲害等自然因素以及無序的海岸帶開發(fā)已經(jīng)造成紅樹林資源的大面積縮減.從宏觀大尺度的視角,定量分析典型紅樹林群落的急速退化特征及其影響因素,可為我國紅樹林資源的保護(hù)及生態(tài)恢復(fù)提供科學(xué)依據(jù).海南東寨港是我國最典型、最原始的天然紅樹林分布區(qū),成為紅樹林生態(tài)系統(tǒng)保護(hù)和恢復(fù)研究的理想場所.通過多期遙感影像的判讀分析顯示,近30年來,東寨港紅樹林面積經(jīng)歷了先增長后減小的變化趨勢,紅樹林斑塊亦呈現(xiàn)破碎化趨勢.人類活動(dòng)對(duì)海岸帶擾動(dòng)強(qiáng)度的增加是導(dǎo)致紅樹林斑塊面積縮減及景觀破碎化最為直接的原因.自1980年東寨港紅樹林自然保護(hù)區(qū)成立以來,保護(hù)區(qū)管理局采取了多項(xiàng)切實(shí)可行的紅樹林保護(hù)和恢復(fù)措施,對(duì)保護(hù)當(dāng)?shù)丶t樹林生態(tài)系統(tǒng)的多樣性、完整性起到了至關(guān)重要的作用.

    近年來,東寨港紅樹林區(qū)出現(xiàn)一定范圍的海洋污損動(dòng)物——團(tuán)水虱危害,已造成了近4 hm2紅樹林枯死現(xiàn)象,引起了政府和社會(huì)的高度關(guān)注[7].從空間分布上來看,紅樹林枯死群落主要分布于地形低洼積水處或污染物傳輸通道兩側(cè);河流型紅樹林生態(tài)系統(tǒng)受到團(tuán)水虱的危害最為嚴(yán)重,同時(shí)近海灣的前沿地段也出現(xiàn)一定程度的危害現(xiàn)象.實(shí)地樣方調(diào)查數(shù)據(jù)顯示,團(tuán)水虱危害呈現(xiàn)明顯的“種群選擇性”,即木欖、尖瓣海蓮等處于地帶性演替后期的紅樹林群落更易遭受團(tuán)水虱的蛀蝕而死亡.結(jié)合GIS空間對(duì)應(yīng)關(guān)系分析發(fā)現(xiàn),沿海陸源氮磷營養(yǎng)物質(zhì)的輸入強(qiáng)度是影響紅樹林退化群落分布的重要因素.東寨港紅樹林區(qū)周邊近1 986. 9 hm2的養(yǎng)殖塘是陸源污染物的主要來源,每年向近海排放大量氮磷營養(yǎng)物質(zhì)以及消毒劑等污染物,對(duì)紅樹林生態(tài)系統(tǒng)的健康構(gòu)成嚴(yán)重威脅.因此,位于地形低洼的積水區(qū)、潮溝兩側(cè)以及污水排放口等特殊地理位置的紅樹林群落應(yīng)該是今后紅樹林退化的重點(diǎn)防控區(qū)域.

    參考文獻(xiàn):

    [1]ALONGI D.Present state and future of the World’s mangrove forests[J].Envir Conserv,2002,29 (3):331-349.

    [2]RAKOTOMAVO A,F(xiàn)ROMARD F.Dynamics of mangrove forests in the Mangoky River delta,Madagascar,under the influence of natural and human factors[J].Forest Ecol Management,2010,259 (6):1161-1169.

    [3]伍淑婕,梁士楚.人類活動(dòng)對(duì)紅樹林生態(tài)系統(tǒng)服務(wù)功能的影響[J].海洋環(huán)境科學(xué),2008,27 (5):537-542.

    [4]趙曉濤,楊威,周丹,等.影響我國河口地區(qū)可持續(xù)發(fā)展的五大問題[J].海洋開發(fā)與管理,2008,25 (3):91-93.

    [5]REDDI E U B,RAMAN A V,SATYANARAYANA B,et al.Degradation of mangrove ecosystem due to hinterland farm practices: A case for Coringa,East Coast of India [J].南京林業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版,2003,27 (2):1-6.

    [6]廖寶文,李玫,陳玉軍,等.海南東寨港紅樹林生態(tài)系統(tǒng)研究[M].青島:中國海洋大學(xué)出版社,2007: 2-31.

    [7]徐蒂,廖寶文,朱寧華,等.海南東寨港紅樹林退化原因初探[J].生態(tài)科學(xué),2014,33 (2):294-300.

    [8]BROOK R A.Discovery of Sphaeroma terebrans,a woodboring isopod,in the red mangrove,Rhizophora mangle,habitat of Northern Florida Bay[J].AMBIO: J Human Envir,2004,33 (3):171-173.

    [9]邱勇,李俊,黃勃,等.影響東寨港紅樹林中光背團(tuán)水虱分布的生態(tài)因子研究[J].海洋科學(xué),2013,37 (4):21-25.

    [10]王胤,左平,黃仲琪,等.海南東寨港紅樹林濕地面積變化及其驅(qū)動(dòng)力分析[J].四川環(huán)境,2006,25 (3):44-49.

    [11]LIU K,LI X,SHI X,et al.Monitoring mangrove forest changes using remote sensing and GIS data with decisiontree learning[J].Wetlands,2008,28 (2):336-346.

    [12]吳哲,陳歆,劉貝貝,等.基于InVEST模型的海南島氮磷營養(yǎng)物質(zhì)負(fù)荷的風(fēng)險(xiǎn)評(píng)估[J].熱帶作物學(xué)報(bào),2013,34 (9):1791-1797.

    [13]李廷友,林振山.海水圍塘混合養(yǎng)殖生態(tài)系統(tǒng)氮磷平衡的研究[J].井岡山大學(xué)學(xué)報(bào):自然科學(xué)版,2010,31 (2):32-35.

    [14]王兵,鄭秋紅,郭浩.基于Shannon-Wiener指數(shù)的中國森林物種多樣性保育價(jià)值評(píng)估方法[J].林業(yè)科學(xué)研究,2008,21(2):268-274.

    [15]薛春汀.人類活動(dòng)對(duì)密克羅尼西亞聯(lián)邦庫賽埃島紅樹林海岸的影響[J].海洋湖沼通報(bào),2002,(2):17-23.

    [16]范航清,劉文愛,鐘才榮,等.中國紅樹林蛀木團(tuán)水虱危害分析研究[J].廣西科學(xué),2014,21 (2):140-152.

    [17]廖寶文,李玫,陳玉軍,等.中國紅樹林恢復(fù)與重建技術(shù)[M].北京:科學(xué)出版社,2010: 120-122.

    【責(zé)任編輯李曉卉】

    Spatial patterns and environmental factors of rapidly degraded mangroves at Dongzhaigang Harbor in Hainan

    SUN Yanwei1,LIAO Baowen1,GUAN Wei1,WANG Rongli1,DING Dongjing1,ZHONG Cairong2
    (1 Research Institute of Tropical Forestry,Chinese Academy of Forestry,Guangzhou 510520,China; 2 Hainan Dongzhaigang National Natural Reserve Authority,Haikou 571129,China)

    Abstract:【Objective】The purpose of this study was to analyze mangrove degradation in some coastlands of southern China.【Method】Mangrove of Dongzhaigang Harbor in Hainan Province was investigated as a case study.Landscape pattern dynamics of mangrove wetland from 1987 to 2013 were analyzed at Dongzhaigang Harbor by interpreting five phase remote sensing images.The degradation characteristics of mangrove forest and its spatial patterns were also further studied using SPOT5 and aerial imagery.Furthermore,the influencing factors of mangrove degradation such as topographical factors and terrestrial contaminant loading were analyzed by surveying numerous sampling plots.【Result and conclusion】It was found that mangrove area decreased from 1 709. 4 hm2in 1999 to 1 679. 5 hm2in 2013,while the area of ponds in the 2 km coastal buffers increased from 59. 1 hm2in 1987 to 1 986. 9 hm2in 2013.The dead area of mangrove patches had reached 4 hm2by the end of 2013.Spatially,mangrove degradation mainly located in low-lying land and both sides of contaminant transport channels.The degradation degree of River-type mangrove communities was the most serious,whose major distribution was in the middle and lower rea-book=112,ebook=384ches of Yanfengdong River.Meanwhile,the degraded mangrove communities also distributed sporadically in frontier areas of Dongzhaigang Harbor.The degradation degree was higher for mangrove communities in the climax stage of succession.The distribution of degraded mangrove communities was mainly determined by terrestrial pollutant loading.

    Key words:Dongzhaigang Harbor; mangrove forest; degradation characteristic; spatial pattern; influencing factor

    基金項(xiàng)目:科技部林業(yè)公益性行業(yè)科研專項(xiàng)(201404305) ;國家自然科學(xué)基金(41176084) ;國家林業(yè)局海南東寨港紅樹林生態(tài)站項(xiàng)目

    作者簡介:孫艷偉(1983—),男,助理研究員,博士,E-mail: sunyanwei2002@ 126.com;通信作者:廖寶文(1963—),男,研究員,博士,E-mail: baowenliao@ ritf.ac.cn

    收稿日期:2014-11-09

    文章編號(hào):1001-411X(2015) 06-0111-08

    文獻(xiàn)標(biāo)志碼:A

    中圖分類號(hào):Q149

    猜你喜歡
    空間分布紅樹林影響因素
    藏著寶藏的紅樹林
    神奇的紅樹林
    走過紅樹林
    歌海(2018年4期)2018-05-14 12:46:15
    白龜山濕地重金屬元素分布特征及其來源分析
    綠色科技(2016年20期)2016-12-27 17:34:13
    基于GIS技術(shù)的福建省柳葉白前資源適宜性空間分布研究
    江蘇省臭氧污染變化特征
    科技視界(2016年18期)2016-11-03 23:51:58
    鐵路客流時(shí)空分布研究綜述
    環(huán)衛(wèi)工人生存狀況的調(diào)查分析
    中國市場(2016年35期)2016-10-19 02:30:10
    農(nóng)業(yè)生產(chǎn)性服務(wù)業(yè)需求影響因素分析
    商(2016年27期)2016-10-17 07:09:07
    村級(jí)發(fā)展互助資金組織的運(yùn)行效率研究
    商(2016年27期)2016-10-17 04:40:12
    啪啪无遮挡十八禁网站| 宅男免费午夜| 亚洲久久久国产精品| 女生性感内裤真人,穿戴方法视频| 热re99久久国产66热| 久久久久久亚洲精品国产蜜桃av| www.999成人在线观看| 国产蜜桃级精品一区二区三区| 韩国精品一区二区三区| 亚洲伊人色综图| www.自偷自拍.com| 99久久人妻综合| 涩涩av久久男人的天堂| 亚洲第一av免费看| 国产免费av片在线观看野外av| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久久av美女十八| 免费少妇av软件| 黄色视频不卡| 久久精品亚洲av国产电影网| 日本三级黄在线观看| 超碰成人久久| 乱人伦中国视频| 亚洲成av片中文字幕在线观看| 18禁美女被吸乳视频| 免费不卡黄色视频| 两性午夜刺激爽爽歪歪视频在线观看 | xxxhd国产人妻xxx| 麻豆成人av在线观看| 国产又爽黄色视频| 丝袜美腿诱惑在线| 免费在线观看亚洲国产| 日韩国内少妇激情av| 99国产极品粉嫩在线观看| 免费高清视频大片| 午夜精品久久久久久毛片777| 欧美中文综合在线视频| 国产免费男女视频| 99re在线观看精品视频| 亚洲av熟女| 免费在线观看亚洲国产| 久久精品国产综合久久久| 国产精品综合久久久久久久免费 | 在线观看免费视频网站a站| 91av网站免费观看| 国产精品久久电影中文字幕| 亚洲熟女毛片儿| 1024香蕉在线观看| 国产三级黄色录像| 国产精华一区二区三区| 久久热在线av| 黄网站色视频无遮挡免费观看| 欧美在线黄色| 高清黄色对白视频在线免费看| 日本vs欧美在线观看视频| 在线观看免费高清a一片| 国产一区二区三区视频了| 宅男免费午夜| 99在线视频只有这里精品首页| 级片在线观看| 亚洲精品中文字幕在线视频| 欧美老熟妇乱子伦牲交| 一区二区三区精品91| 级片在线观看| 午夜视频精品福利| 久久天堂一区二区三区四区| 欧美日韩乱码在线| 精品午夜福利视频在线观看一区| 老司机亚洲免费影院| 在线观看一区二区三区激情| 亚洲一码二码三码区别大吗| 久久热在线av| 日韩精品青青久久久久久| 亚洲人成77777在线视频| aaaaa片日本免费| 好看av亚洲va欧美ⅴa在| 男女午夜视频在线观看| 51午夜福利影视在线观看| 久久草成人影院| 免费日韩欧美在线观看| 亚洲精品一二三| 男人的好看免费观看在线视频 | 欧美黑人欧美精品刺激| 国产欧美日韩综合在线一区二区| 国产1区2区3区精品| 亚洲全国av大片| 国产黄a三级三级三级人| 99热国产这里只有精品6| 可以在线观看毛片的网站| 午夜精品国产一区二区电影| 亚洲性夜色夜夜综合| e午夜精品久久久久久久| 一本大道久久a久久精品| 搡老岳熟女国产| 99国产精品一区二区蜜桃av| 国产精品一区二区免费欧美| 精品人妻1区二区| 久久中文字幕人妻熟女| 欧美日韩亚洲国产一区二区在线观看| 亚洲午夜精品一区,二区,三区| 久久国产精品男人的天堂亚洲| 亚洲国产精品一区二区三区在线| 校园春色视频在线观看| 黄色 视频免费看| 日本黄色日本黄色录像| 桃色一区二区三区在线观看| 9热在线视频观看99| 午夜福利一区二区在线看| 丝袜人妻中文字幕| 99精品在免费线老司机午夜| 亚洲国产精品sss在线观看 | 色婷婷久久久亚洲欧美| 日韩精品免费视频一区二区三区| 夜夜爽天天搞| av国产精品久久久久影院| 亚洲专区中文字幕在线| 69av精品久久久久久| 亚洲精品国产色婷婷电影| 久久国产精品影院| 国产免费现黄频在线看| 一级片'在线观看视频| 一级a爱视频在线免费观看| 免费在线观看完整版高清| 91国产中文字幕| 精品国产乱子伦一区二区三区| 久久久国产成人精品二区 | 丰满饥渴人妻一区二区三| 又黄又爽又免费观看的视频| 不卡一级毛片| 一级a爱片免费观看的视频| 免费高清视频大片| a级毛片在线看网站| av片东京热男人的天堂| 国产aⅴ精品一区二区三区波| 久久国产亚洲av麻豆专区| 一级a爱视频在线免费观看| 女人高潮潮喷娇喘18禁视频| 国产精品国产av在线观看| 国产有黄有色有爽视频| 欧美日韩一级在线毛片| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品综合一区在线观看 | 日韩大码丰满熟妇| 亚洲欧美一区二区三区久久| 美女午夜性视频免费| 免费看十八禁软件| 色播在线永久视频| 1024香蕉在线观看| 欧美中文日本在线观看视频| 一二三四在线观看免费中文在| 欧美午夜高清在线| 十分钟在线观看高清视频www| 国产伦一二天堂av在线观看| 国产蜜桃级精品一区二区三区| 久久久久国内视频| 欧美乱色亚洲激情| 99国产精品免费福利视频| 国产麻豆69| 日韩国内少妇激情av| 亚洲激情在线av| 99香蕉大伊视频| 国产伦一二天堂av在线观看| 国产蜜桃级精品一区二区三区| 中文欧美无线码| 午夜91福利影院| 在线观看免费视频网站a站| 如日韩欧美国产精品一区二区三区| 黄片播放在线免费| 亚洲成人精品中文字幕电影 | 日韩免费高清中文字幕av| 啦啦啦免费观看视频1| 久久性视频一级片| 欧美日韩瑟瑟在线播放| 久久久国产精品麻豆| 可以免费在线观看a视频的电影网站| 婷婷丁香在线五月| 欧美最黄视频在线播放免费 | 国产精品久久久久久人妻精品电影| 欧美成狂野欧美在线观看| 午夜免费成人在线视频| 国产精品日韩av在线免费观看 | 精品第一国产精品| 国产精品秋霞免费鲁丝片| 黄频高清免费视频| 视频区欧美日本亚洲| 中文字幕人妻丝袜一区二区| 91成人精品电影| 久久久久久人人人人人| 亚洲三区欧美一区| 欧美成人午夜精品| 欧美 亚洲 国产 日韩一| 精品国产乱码久久久久久男人| 伦理电影免费视频| 黑人操中国人逼视频| 一进一出好大好爽视频| 一级,二级,三级黄色视频| 日韩有码中文字幕| 久久国产精品人妻蜜桃| 91精品国产国语对白视频| 久久久精品欧美日韩精品| 午夜老司机福利片| 免费在线观看日本一区| 韩国精品一区二区三区| 久久精品国产亚洲av高清一级| av网站在线播放免费| 91老司机精品| 欧美国产精品va在线观看不卡| 国产一区二区三区在线臀色熟女 | 法律面前人人平等表现在哪些方面| 亚洲欧美一区二区三区久久| 亚洲精品一卡2卡三卡4卡5卡| 国产不卡一卡二| 久久精品aⅴ一区二区三区四区| 久热这里只有精品99| 久久久水蜜桃国产精品网| 亚洲欧美一区二区三区黑人| 如日韩欧美国产精品一区二区三区| 精品一区二区三区四区五区乱码| 亚洲欧洲精品一区二区精品久久久| 国产av又大| 欧美一级毛片孕妇| 丰满人妻熟妇乱又伦精品不卡| 国产精品秋霞免费鲁丝片| 极品人妻少妇av视频| 男人操女人黄网站| 青草久久国产| 久久精品成人免费网站| 在线观看午夜福利视频| 亚洲aⅴ乱码一区二区在线播放 | 天天躁夜夜躁狠狠躁躁| 涩涩av久久男人的天堂| 国产精品98久久久久久宅男小说| 久久精品国产99精品国产亚洲性色 | 亚洲一区二区三区不卡视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产在线观看jvid| 怎么达到女性高潮| 亚洲欧美激情综合另类| 精品久久蜜臀av无| 精品久久久久久电影网| 嫩草影院精品99| 亚洲国产精品999在线| 国产视频一区二区在线看| 女人高潮潮喷娇喘18禁视频| 在线观看免费视频日本深夜| 在线十欧美十亚洲十日本专区| 高潮久久久久久久久久久不卡| 亚洲成a人片在线一区二区| 男人操女人黄网站| 女人精品久久久久毛片| 亚洲免费av在线视频| 天堂动漫精品| 在线永久观看黄色视频| 一进一出好大好爽视频| 啦啦啦 在线观看视频| 欧美日本中文国产一区发布| 亚洲精品美女久久av网站| 亚洲av美国av| 超色免费av| 亚洲欧美一区二区三区黑人| 亚洲激情在线av| 韩国精品一区二区三区| 国产亚洲精品一区二区www| 别揉我奶头~嗯~啊~动态视频| 国产又色又爽无遮挡免费看| 色播在线永久视频| 国产精品偷伦视频观看了| 国产不卡一卡二| 久久国产精品影院| 黄色怎么调成土黄色| 极品人妻少妇av视频| 日日干狠狠操夜夜爽| 啦啦啦免费观看视频1| 一二三四在线观看免费中文在| 99久久综合精品五月天人人| 中文字幕另类日韩欧美亚洲嫩草| 精品国产一区二区三区四区第35| 精品一区二区三区视频在线观看免费 | 女警被强在线播放| 日日干狠狠操夜夜爽| 夜夜爽天天搞| 国产成人影院久久av| videosex国产| 免费女性裸体啪啪无遮挡网站| 国产精品香港三级国产av潘金莲| 久久久久亚洲av毛片大全| 亚洲人成网站在线播放欧美日韩| 亚洲va日本ⅴa欧美va伊人久久| 亚洲自偷自拍图片 自拍| 在线国产一区二区在线| 国产1区2区3区精品| 日韩大尺度精品在线看网址 | 51午夜福利影视在线观看| 国产亚洲欧美在线一区二区| 国产主播在线观看一区二区| 亚洲成人免费av在线播放| 国产一区二区激情短视频| 久久性视频一级片| 熟女少妇亚洲综合色aaa.| 日韩一卡2卡3卡4卡2021年| 最近最新免费中文字幕在线| 国产精品 国内视频| 一夜夜www| 一个人观看的视频www高清免费观看 | 国产又色又爽无遮挡免费看| 亚洲少妇的诱惑av| 精品电影一区二区在线| 涩涩av久久男人的天堂| 在线永久观看黄色视频| 日韩欧美免费精品| 日韩人妻精品一区2区三区| 精品国产一区二区久久| www.熟女人妻精品国产| 国产精品电影一区二区三区| 久久香蕉精品热| 另类亚洲欧美激情| 午夜a级毛片| 国产蜜桃级精品一区二区三区| www.999成人在线观看| 亚洲黑人精品在线| 中文欧美无线码| 色婷婷av一区二区三区视频| 亚洲第一欧美日韩一区二区三区| 男人操女人黄网站| 久久久精品国产亚洲av高清涩受| 1024视频免费在线观看| av福利片在线| 久久香蕉激情| 一个人观看的视频www高清免费观看 | www.精华液| 中文字幕精品免费在线观看视频| 又紧又爽又黄一区二区| 一a级毛片在线观看| 黄色视频,在线免费观看| 视频区图区小说| 国产在线精品亚洲第一网站| 亚洲av五月六月丁香网| 日本黄色视频三级网站网址| 免费少妇av软件| 国产精品av久久久久免费| 老司机午夜十八禁免费视频| 一级毛片女人18水好多| 激情在线观看视频在线高清| 日韩欧美国产一区二区入口| 国产一区二区激情短视频| 一边摸一边抽搐一进一小说| 亚洲avbb在线观看| 人妻丰满熟妇av一区二区三区| 日日爽夜夜爽网站| 又黄又爽又免费观看的视频| 国产精品久久久av美女十八| 国产精品日韩av在线免费观看 | 女性生殖器流出的白浆| 18美女黄网站色大片免费观看| 99久久精品国产亚洲精品| 两个人看的免费小视频| 亚洲自偷自拍图片 自拍| 国产精品1区2区在线观看.| 国产精品影院久久| 天堂动漫精品| 一边摸一边抽搐一进一小说| 国产人伦9x9x在线观看| 一级毛片女人18水好多| 老司机深夜福利视频在线观看| 亚洲一区中文字幕在线| 国产成人系列免费观看| 欧美日韩乱码在线| √禁漫天堂资源中文www| 国产伦一二天堂av在线观看| 亚洲自拍偷在线| 成在线人永久免费视频| 美女扒开内裤让男人捅视频| 国产精品免费视频内射| 19禁男女啪啪无遮挡网站| 桃色一区二区三区在线观看| 国产熟女xx| 涩涩av久久男人的天堂| 在线看a的网站| 在线av久久热| 亚洲欧美日韩另类电影网站| 在线永久观看黄色视频| 免费搜索国产男女视频| 视频区图区小说| 99国产精品99久久久久| 亚洲伊人色综图| 亚洲人成77777在线视频| 日本wwww免费看| 中文字幕精品免费在线观看视频| 悠悠久久av| 久久香蕉国产精品| 操美女的视频在线观看| 在线观看舔阴道视频| 啦啦啦免费观看视频1| 桃红色精品国产亚洲av| 精品国产乱子伦一区二区三区| 曰老女人黄片| av中文乱码字幕在线| 国产精品爽爽va在线观看网站 | 免费搜索国产男女视频| 成年女人毛片免费观看观看9| 黄色毛片三级朝国网站| 悠悠久久av| 免费高清视频大片| 久久这里只有精品19| 老司机深夜福利视频在线观看| 亚洲精品一区av在线观看| 欧美成人性av电影在线观看| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 午夜福利一区二区在线看| 一区在线观看完整版| 国产高清国产精品国产三级| 变态另类成人亚洲欧美熟女 | 美女国产高潮福利片在线看| 国产高清激情床上av| 国产精品偷伦视频观看了| 看片在线看免费视频| 新久久久久国产一级毛片| 国产成人精品无人区| 黄片播放在线免费| 精品电影一区二区在线| 在线看a的网站| 国内久久婷婷六月综合欲色啪| 国产精品久久久久成人av| 啦啦啦免费观看视频1| 视频在线观看一区二区三区| 91精品国产国语对白视频| 亚洲欧美日韩无卡精品| 一级片免费观看大全| 最近最新中文字幕大全电影3 | 亚洲美女黄片视频| 精品福利永久在线观看| 天堂影院成人在线观看| 日韩免费av在线播放| 叶爱在线成人免费视频播放| 亚洲国产中文字幕在线视频| 999精品在线视频| 亚洲av电影在线进入| 女人精品久久久久毛片| 国产av精品麻豆| 精品久久久久久电影网| 女性被躁到高潮视频| 国产91精品成人一区二区三区| 日韩人妻精品一区2区三区| 亚洲专区中文字幕在线| 国产蜜桃级精品一区二区三区| 久久精品91无色码中文字幕| 日日爽夜夜爽网站| 亚洲国产精品sss在线观看 | 69精品国产乱码久久久| 又黄又粗又硬又大视频| 一二三四社区在线视频社区8| 视频区图区小说| 9热在线视频观看99| 亚洲熟妇中文字幕五十中出 | 久热爱精品视频在线9| 国产精品九九99| 久久影院123| 国产野战对白在线观看| 日韩高清综合在线| 欧美人与性动交α欧美软件| 国产主播在线观看一区二区| 一区在线观看完整版| 欧美日韩av久久| 国产黄a三级三级三级人| 久久久国产成人精品二区 | 黑丝袜美女国产一区| 精品久久久久久成人av| 欧美亚洲日本最大视频资源| 成在线人永久免费视频| 国产aⅴ精品一区二区三区波| 纯流量卡能插随身wifi吗| 亚洲av第一区精品v没综合| 免费观看精品视频网站| 男女午夜视频在线观看| 中文欧美无线码| 激情在线观看视频在线高清| 美女高潮到喷水免费观看| 日本三级黄在线观看| 好男人电影高清在线观看| a级毛片在线看网站| 色播在线永久视频| 亚洲成人国产一区在线观看| 首页视频小说图片口味搜索| 国产高清视频在线播放一区| 男女下面插进去视频免费观看| 久久午夜综合久久蜜桃| 伊人久久大香线蕉亚洲五| 在线观看午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看 | 精品一区二区三区av网在线观看| 久久婷婷成人综合色麻豆| 久久精品成人免费网站| 免费在线观看日本一区| 99re在线观看精品视频| tocl精华| 欧美日韩瑟瑟在线播放| 美女高潮喷水抽搐中文字幕| 成人三级做爰电影| 亚洲国产欧美日韩在线播放| 51午夜福利影视在线观看| 亚洲色图综合在线观看| 在线观看免费视频网站a站| 国产欧美日韩综合在线一区二区| 美女福利国产在线| 久久天堂一区二区三区四区| 免费看十八禁软件| 50天的宝宝边吃奶边哭怎么回事| 亚洲全国av大片| 中文字幕另类日韩欧美亚洲嫩草| 亚洲第一欧美日韩一区二区三区| 国产视频一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片女人18水好多| 一区福利在线观看| 国产成人av教育| 女人被躁到高潮嗷嗷叫费观| 这个男人来自地球电影免费观看| 久久青草综合色| 免费高清视频大片| 丰满迷人的少妇在线观看| 日韩精品免费视频一区二区三区| 一进一出抽搐动态| 成年版毛片免费区| 国产精品爽爽va在线观看网站 | 国产高清videossex| 神马国产精品三级电影在线观看 | 久久99一区二区三区| 精品一区二区三区av网在线观看| av中文乱码字幕在线| 欧美人与性动交α欧美精品济南到| 精品卡一卡二卡四卡免费| 亚洲第一青青草原| 高清av免费在线| 日本五十路高清| 国产欧美日韩一区二区精品| 亚洲色图av天堂| 成人亚洲精品一区在线观看| 精品一区二区三区av网在线观看| 最好的美女福利视频网| 国产熟女午夜一区二区三区| 欧美日韩精品网址| 国产av精品麻豆| 黄色怎么调成土黄色| 国产91精品成人一区二区三区| 久久香蕉激情| 亚洲人成网站在线播放欧美日韩| 亚洲人成77777在线视频| 欧美久久黑人一区二区| 少妇被粗大的猛进出69影院| 欧美中文日本在线观看视频| 久久人人97超碰香蕉20202| 真人一进一出gif抽搐免费| 亚洲国产精品合色在线| 国产有黄有色有爽视频| 脱女人内裤的视频| 岛国视频午夜一区免费看| 国产精品 欧美亚洲| 亚洲成人免费av在线播放| 一进一出抽搐动态| 每晚都被弄得嗷嗷叫到高潮| 桃红色精品国产亚洲av| 水蜜桃什么品种好| 夜夜爽天天搞| 亚洲国产精品合色在线| 高清黄色对白视频在线免费看| 日日爽夜夜爽网站| 深夜精品福利| 天堂影院成人在线观看| 日本 av在线| 国产亚洲欧美精品永久| 搡老岳熟女国产| 成人国语在线视频| 人人妻,人人澡人人爽秒播| 无人区码免费观看不卡| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费日韩欧美大片| 夜夜看夜夜爽夜夜摸 | 香蕉久久夜色| 91成年电影在线观看| 久久国产亚洲av麻豆专区| 免费在线观看完整版高清| 久久国产精品人妻蜜桃| 男女做爰动态图高潮gif福利片 | 久久久水蜜桃国产精品网| 欧美大码av| 在线永久观看黄色视频| 国产av一区在线观看免费| 午夜福利欧美成人| www国产在线视频色| 欧美在线黄色| 大码成人一级视频| 在线视频色国产色| 亚洲欧美精品综合久久99| 亚洲国产毛片av蜜桃av| 久久久久国产一级毛片高清牌| 欧美亚洲日本最大视频资源| 在线观看免费日韩欧美大片| 久久精品91蜜桃| 十八禁人妻一区二区| 国产精品99久久99久久久不卡| 国产精品av久久久久免费| 国产熟女xx| 天天添夜夜摸| 美国免费a级毛片| 亚洲狠狠婷婷综合久久图片| 在线观看免费高清a一片| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲专区中文字幕在线| 在线观看日韩欧美| 神马国产精品三级电影在线观看 | 午夜精品在线福利| 免费在线观看影片大全网站| 免费在线观看黄色视频的| 91字幕亚洲|