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

    廟島群島南五島生態(tài)系統(tǒng)凈初級生產(chǎn)力空間分布及其影響因子

    2015-02-13 01:09:21石洪華王曉麗豐愛平
    生態(tài)學(xué)報 2015年24期
    關(guān)鍵詞:坡向海島生產(chǎn)力

    池 源,石洪華,*,王曉麗,李 捷,豐愛平

    1國家海洋局第一海洋研究所, 青島 266061 2天津理工大學(xué)環(huán)境科學(xué)與安全工程學(xué)院,天津 300384 3青島理工大學(xué)環(huán)境與市政工程學(xué)院, 青島 266033

    廟島群島南五島生態(tài)系統(tǒng)凈初級生產(chǎn)力空間分布及其影響因子

    池 源1,石洪華1,*,王曉麗2,李 捷3,豐愛平1

    1國家海洋局第一海洋研究所, 青島 266061 2天津理工大學(xué)環(huán)境科學(xué)與安全工程學(xué)院,天津 300384 3青島理工大學(xué)環(huán)境與市政工程學(xué)院, 青島 266033

    凈初級生產(chǎn)力(NPP)估算對于海島碳源/匯研究具有重要意義。以廟島群島南五島為例,結(jié)合CASA模型和區(qū)域特征構(gòu)建NPP估算模型,借助RS和GIS技術(shù)進(jìn)行NPP估算,進(jìn)而分析南五島NPP空間分布特征及其影響因子。結(jié)果表明:南五島NPP總量為11043.52t C/a,平均密度為340.19g C m-2a-1,處于全國平均水平,高于同緯度的西部地區(qū),低于東部沿海大陸地區(qū);夏季NPP總量占全年的80%左右,春季和秋季分別占11%和7%,冬季僅占1.3%;不同海島的NPP平均密度由大到小依次為大黑山島、北長山島、廟島、南長山島和小黑山島,各島NPP平均密度與建設(shè)用地比例呈明顯負(fù)相關(guān);不同地表覆蓋類型的NPP平均密度由大到小依次為闊葉林、針葉林、農(nóng)田、草地、建設(shè)用地和裸地,林地具有較高的NPP值,說明南五島的人工林建設(shè)具有重要生態(tài)作用;NDVI和地表覆蓋類型是NPP最主要的影響參數(shù),地形參數(shù)通過影響NDVI和地表覆蓋類型間接作用于NPP結(jié)果;NPP與土壤pH、有效磷、全磷、全鉀呈顯著負(fù)相關(guān),與全氮、總碳、總有機(jī)碳呈顯著正相關(guān),與含水量、速效鉀和含鹽量之間相關(guān)關(guān)系不明顯。

    海島生態(tài)系統(tǒng);凈初級生產(chǎn)力;廟島群島南五島;空間分布;影響因子

    凈初級生產(chǎn)力(NPP)是植物光合作用產(chǎn)生有機(jī)質(zhì)總量扣除自養(yǎng)呼吸后的剩余部分,是地表碳循環(huán)的重要組成[1]。NPP不僅能夠直接反映植被群落在自然環(huán)境條件下的生產(chǎn)能力,表征生態(tài)系統(tǒng)的質(zhì)量狀況,而且是判定生態(tài)系統(tǒng)碳源/匯和調(diào)節(jié)生態(tài)過程的主要因子[2]。同時,3S技術(shù)的發(fā)展為NPP研究提供了新的思路和方法,從遙感影像中獲取地表覆蓋類型并提取歸一化植被指數(shù)(NDVI),結(jié)合氣象、地形等數(shù)據(jù)對區(qū)域NPP進(jìn)行估算,具有較高的準(zhǔn)確度和較強(qiáng)的可操作性,在草地、森林、濕地、農(nóng)田等不同生態(tài)系統(tǒng)以及不同時空尺度的NPP研究中得到了廣泛的應(yīng)用[3-7],然而鮮有關(guān)于海島生態(tài)系統(tǒng)NPP的研究。

    海島兼具陸、海雙重特征,是一種特殊的生態(tài)系統(tǒng)[8]。由于海島位置的特殊性、面積的有限性以及空間的隔離性,海島生態(tài)系統(tǒng)更容易受到風(fēng)暴潮、干旱等自然災(zāi)害的干擾,并且具有明顯的獨(dú)立性和脆弱性[9]。同時,海島森林、灌草、土壤等是生物多樣性維持、水源涵養(yǎng)、水土保持等生態(tài)功能發(fā)揮的基本載體[10]。全球海島面積約占全球陸地面積的1/15,在碳循環(huán)中扮演著非常重要的角色[11],海島NPP研究具有重要意義。

    本文以廟島群島南五島為例,構(gòu)建NPP估算模型,在現(xiàn)場調(diào)查的基礎(chǔ)上,采用不同季節(jié)遙感影像數(shù)據(jù),結(jié)合氣象資料和地形數(shù)據(jù),對南五島NPP的空間分布及其影響因子進(jìn)行研究,以期為南五島生態(tài)系統(tǒng)碳源/匯研究提供依據(jù)。

    1 研究區(qū)概況

    廟島群島位于山東半島北側(cè),黃、渤海交會處,是我國北方海島分布的典型區(qū)域,也是長島國家級自然保護(hù)區(qū)的重要依托,按照地理位置可分為南、中、北3個島群。其中南部島群的海島分布最為集中,距離大陸最近,是山東省長島縣的人口、經(jīng)濟(jì)和文化中心,其主體為南長山島、北長山島、廟島、大黑山島和小黑山島5個有居民海島,合稱“南五島”。該區(qū)域?qū)儆跂|亞季風(fēng)氣候區(qū),年均氣溫12.0℃,1月平均氣溫-1.6℃,7月平均氣溫24.5℃;年均降水量約537mm,降水多集中在6—9月;日照較為充足,年均日照時數(shù)2612h。南五島地勢起伏明顯,山勢大致呈南北走向,最高點海拔約為189m;土壤主要有棕壤、褐土、潮土三大類,以棕壤土分布面積最大,土層厚度約為30cm,多砂礫,土質(zhì)較差[12]。南五島現(xiàn)有林木主要為人工林,以黑松和刺槐為主要優(yōu)勢樹種,林下分布著紫穗槐、檉柳等灌木和各類草本植物。

    2 材料與方法

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

    2.1.1 遙感影像

    采用LANDSAT8衛(wèi)星2013年4月21日、8月11日、11月15日和2014年1月2日(代表不同季節(jié))4個時相南五島所在區(qū)域30m分辨率的無云影像。

    利用ENVI4.7軟件對影像進(jìn)行裁切、輻射定標(biāo)、波段運(yùn)算得到NDVI值(圖1)?;?013年8月遙感影像,通過ArcGIS10.0軟件進(jìn)行人機(jī)交互解譯,將南五島地表覆蓋分為針葉林、闊葉林、草地、農(nóng)田、建設(shè)用地和裸地6類(圖2),結(jié)合現(xiàn)場實地調(diào)研、Google Earth和相關(guān)的圖集資料進(jìn)行解譯精度驗證,解譯精度為92.8%,能夠滿足本次研究的需要。

    圖1 不同季節(jié)南五島NDVI分布Fig.1 NDVI Distribution of five southern islands in different seasons

    圖2 南五島地表覆蓋類型Fig.2 Land cover types of five southern islands

    2.1.2 地形

    采用2011年公布的Aster GDEM第二版DEM數(shù)據(jù),垂直分辨率20m,水平分辨率30m;通過ArcGIS10.0由DEM數(shù)據(jù)提取出高程、坡度和坡向(圖3)。

    2.1.3 氣象和土壤

    降雨量、氣溫、日照時數(shù)、相對濕度來自長島縣氣象站監(jiān)測數(shù)據(jù);太陽總輻射來自煙臺福山氣象站監(jiān)測的多年平均數(shù)據(jù)。

    2012年夏季進(jìn)行土壤現(xiàn)場采樣并檢測,獲得南五島36個樣方的土壤理化性質(zhì)數(shù)據(jù)(圖4)。采用多點混合取土樣方法,在每個樣方內(nèi)分別選取3個取樣點,采用標(biāo)準(zhǔn)環(huán)刀(100cm3)進(jìn)行取樣,將3個取樣點的土壤均勻混合,作為該樣方的土壤樣品。樣品采集后,帶回實驗室,除去其中草根、大石礫等雜質(zhì),過2mm鋼篩并磨細(xì),在105℃高溫下烘至恒重,稱量并計算土壤含水率;pH值采用電位法進(jìn)行測定,有效磷采用鹽酸-氟化銨法進(jìn)行測定,速效鉀采用乙酸銨提取法進(jìn)行測定,全磷采用碳酸鈉熔融法進(jìn)行測定,全鉀采用氫氧化鈉熔融法進(jìn)行測定,含鹽量采用電導(dǎo)法進(jìn)行測定,全氮采用開氏法進(jìn)行測定,總碳采用高溫灼燒法進(jìn)行測定,總有機(jī)碳采用重鉻酸鉀容量法進(jìn)行測定。

    圖3 南五島高程、坡度和坡向Fig.3 Altitude, slope and aspect of five southern islands

    圖4 土壤調(diào)查樣方位置Fig.4 Locations of soil sample sites

    2.2 NPP模型構(gòu)建與計算

    根據(jù)CASA模型[13],構(gòu)建本次NPP遙感估算總模型如下:

    NPP(x, t)=APAR(x, t)×ξ(x, t)

    (1)

    APAR(x, t)=PAR(x, t)×FPAR(x, t)

    (2)

    ξ(x, t)=ft(t)×fw(t)×ξmax

    (3)

    式中,NPP(x,t)為x點t月凈初級生產(chǎn)力;APAR(x,t)為x點t月吸收的光合有效輻射(MJ m-2月-1);ξ(x,t) 為x點t月的實際光能利用率(g C/MJ);PAR(x,t)為x點t月的光合有效輻射(MJ m-2月-1);FPAR(x,t)為x點t月光合有效輻射吸收比例(%);ft(t)和fw(t)分別為研究區(qū)t月的氣溫脅迫因子和水分脅迫因子(%);ξmax為植被最大光能利用率(g C/MJ)。

    2.2.1 PAR計算

    考慮到南五島高程起伏明顯,地形遮蔽對于不同位置能夠接受到的太陽輻射量有著直接影響[14],將地形作為PAR計算的重要因子,構(gòu)建公式如下:

    PAR(x, t)=SOL(t) ×50%×[0.4+0.6×d(x, t)]

    (4)

    式中,SOL(t)為t月太陽總輻射量(MJ m-2月-1);50%表示植被能利用的太陽有效輻射占太陽總輻射的比例,由文獻(xiàn)[5]得出;0.4和0.6分別為區(qū)域太陽散射輻射和直接輻射占太陽總輻射的多年平均比例,由文獻(xiàn)[15-16]得出;d(x,t)為x點t月的太陽輻射地形影響因子,由下式計算得出:

    (5)

    式中,θ(t)為t月遙感影像獲得當(dāng)天南五島所在緯度的正午太陽高度角,Hillshade(x,t)為x點t月的遮蔽度,無量綱,Hillshade(max,t)和Hillshade(min,t)分別為t月遮蔽度的最大值和最小值。Hillshade(x,t)由下式求得[17]:

    (6)

    其中,Slope(x)為點x的坡度,Azimuth為太陽方位角,這里均取180°,Aspect(x)為x點的坡向。

    2.2.2 FPAR計算

    FPAR與NDVI存在明顯的線性關(guān)系[18],可由下式得出:

    FPAR(x,t)1=(NDVI(x,t)-NDVImin)/(NDVImax-NDVImin) ×(FPARmax-FRARmin)+ FPARmin

    (7)

    式中:NDVI(x,t)為x點t月的NDVI值;為了剔除異常值,削弱極值的影響,NDVImax和NDVImin分別取全部月份NDVI值的95%和5%百分位值,F(xiàn)PARmax和FPARmin分別取0.95和0.001[13]。

    同時,研究發(fā)現(xiàn)FPAR與比植被系數(shù)(SR)也具有明顯的線性相關(guān)[19],可由下式得出:

    FPAR(x,t)2=(SR(x,t)-SRmin)/(SRmax-SRmin)×(FPARmax-FPARmin)+FPARmin

    (8)

    式中,SR(x,t)為x點t月的比植被系數(shù),由式9得出。SRmax和SRmin分別取SR 值的95%和5%百分位值。

    (9)

    式7和式8是基于遙感影像像元NDVI值的線性公式,具有跨尺度的特點,能夠運(yùn)用到本次的研究中[13]。同時,由于式7的計算值往往比FPAR實測值高,而式8計算值比實測值低[5],本文同時結(jié)合兩種方法進(jìn)行計算:

    FPAR(x,t)=[FPAR(x,t)1+FPAR(x,t)2]/2

    (10)

    2.2.3 ft和fw計算

    ft由以下方法得出:

    ft(t)=ft(t)1×ft(t)2

    (11)

    ft(t)1反映在不同的最適氣溫情況下植物內(nèi)在的生化作用對光合的限制從而帶來的對光能利用率的影響,由下式求得[19]:

    ft(t)1=0.8+0.02×Topt-0.0005×Topt2

    (12)

    式中,Topt為最適氣溫,取NDVI平均值最高月份的月平均氣溫。當(dāng)月平均氣溫小于或等于-10℃時,ft(t)1取0。

    ft(t)2表示氣溫與最適氣溫偏離時光能利用率減小的趨勢,由下式求得[13,19]:

    (13)

    式中,T為當(dāng)月平均氣溫。當(dāng)某月平均氣溫T比最適氣溫Topt高10℃或低13℃時,該月的ft(t)2值等于月平均氣溫T為最適氣溫Topt時ft(t)2值的一半。

    fw反映了植物所能利用的有效水分條件對光能利用率的影響,由以下公式計算[20]:

    fw(t)=0.5+0.5×E/Ep

    (14)

    式中,E為區(qū)域?qū)嶋H蒸散量,Ep為區(qū)域潛在蒸散量。

    E根據(jù)周廣勝和張新時[21]建立的區(qū)域?qū)嶋H蒸散模型求?。?/p>

    (15)

    式中,r為降水量,Rn為凈輻射量。

    Rn參考《噴灌工程設(shè)計手冊》[22],由下式求得:

    Rn=Rn1-Rn2

    (16)

    Rn1=(1-a)(0.25+0.5n/N)Ra

    (17)

    (18)

    式中,Rn為凈輻射(MJ m-2d-1);a為反射率,取23%;n為實際日照實數(shù),N為該緯度最大日照時數(shù);Ra為大氣頂層的太陽輻射,可查表獲得;σ為斯蒂芬-玻爾茲曼常數(shù),取2×10-9mm d-1K-4,Tk為絕對溫度表示的該月平均溫度(K);ed為水汽壓,可由相對濕度求得。

    Ep由周廣勝、張新時[23]提出的Ep-Rn關(guān)系式求得

    (19)

    2.2.4 ξmax的獲取

    最大光能利用率(ξmax)的取值對NPP結(jié)果有著直接的影響,其具體的取值根據(jù)不同植被類型而有所差異。在諸多關(guān)于ξmax的研究中,Running[4]和朱文泉[5]的研究成果在國內(nèi)外NPP模擬中得到了廣泛的應(yīng)用。前者以生態(tài)生理過程模型模擬了全球10種植被類型的ξmax,但其對于中國的植被而言偏高[24];后者基于誤差最小原則,采用NPP實測數(shù)據(jù)對中國各類植被ξmax進(jìn)行模擬,但由于分辨率過低、混合像元等問題導(dǎo)致在較小空間尺度研究中具有一定的局限性,主要表現(xiàn)為模擬值較實際值偏小[6,25]。本次研究的南五島屬于小空間尺度的研究,同時結(jié)合Running和朱文泉的研究,得到南五島各類植被的ξmax取值(表1)。此外,由于本文像元大小為30m,建設(shè)用地和裸地中存在的部分綠色植被無法識別出來,因此賦予最低的ξmax[2,13,19]。

    表1 最大光能利用率取值Table 1 Values of maximum light use efficiency (ξmax)

    2.2.5 NPP計算

    分別以4月、8月、11月、1月代表春、夏、秋、冬四季,計算南五島不同季節(jié)的NPP值,進(jìn)而得到全年NPP估算值。

    2.3 影響因子分析

    分析南五島NPP與影響因子的關(guān)系以探究南五島NPP的空間分布規(guī)律,為了減少不必要的數(shù)據(jù)量,同時一定程度地削弱NPP極值的影響,將計算結(jié)果重采樣為大小為50m×50m的單元進(jìn)行分析。

    2.3.1 對空間分異性參數(shù)的依賴性

    NPP模型中的計算參數(shù)可分為空間同一性參數(shù)和空間分異性參數(shù)兩類,前者包括氣溫、水汽壓、降雨量、太陽總輻射等氣象因子,相同季節(jié)內(nèi)整個研究區(qū)均采用同一數(shù)值;后者為NDVI、高程、坡度、坡向、地表覆蓋類型,隨著點位的不同而具有差異。分析各點位NPP與空間分異性參數(shù)的相關(guān)關(guān)系,評價各參數(shù)對NPP的影響程度大小。

    通過ArcGIS10.0的Fishnet工具建立50×50m網(wǎng)格,將研究區(qū)域劃分為12094個網(wǎng)格,提取網(wǎng)格中心點所在像元的NDVI、高程、坡度、坡向、地表覆蓋類型的參數(shù)數(shù)據(jù)和NPP數(shù)據(jù)。

    (1)數(shù)據(jù)標(biāo)準(zhǔn)化處理

    對NDVI、高程、坡度按照下式進(jìn)行標(biāo)準(zhǔn)化:

    (20)

    式中,CSi,x為參數(shù)i在x點的標(biāo)準(zhǔn)化數(shù)值,Ci,x為參數(shù)i在x點的原始值,Ci,max為參數(shù)i最大原始值,Ci,min為參數(shù)i最小原始值。

    坡向原始值為按0—360順時針增大,0為正北,180為正南,以向陽性為原則,按照下式進(jìn)行標(biāo)準(zhǔn)化:

    (21)

    式中,ASx為x點標(biāo)準(zhǔn)化坡向值,Ax為x點原始坡向值。

    地表覆蓋類型參照2.2.4的方法賦值。

    (2)簡單相關(guān)分析與偏相關(guān)分析

    運(yùn)用SPSS18軟件,對NPP與NDVI、高程、坡度、坡向和地表覆蓋類型之間兩兩的簡單相關(guān)性進(jìn)行分析。由于各參數(shù)之間也可能存在相互影響,在控制其他參數(shù)影響的環(huán)境下,分析NPP與各參數(shù)的偏相關(guān)性。

    2.3.2 與土壤理化性質(zhì)的相關(guān)性

    基于SPSS18軟件,將36個點位的土壤理化性質(zhì)數(shù)據(jù)與該點位所在50m×50m像元的夏季NPP進(jìn)行相關(guān)性分析,土壤理化性質(zhì)包括pH、含水率、有效磷、速效鉀、全磷、全鉀、含鹽量、全氮、總碳、總有機(jī)碳10項因子。

    3 結(jié)果

    經(jīng)計算得出,南五島全年NPP總量為11043.52t C/a,NPP密度介于0—863.20g C m-2a-1之間,平均密度為340.19g C m-2a-1(圖5)。

    圖5 全年南五島NPP分布特征Fig.5 NPP distribution of five southern islands annually

    3.1 不同季節(jié)NPP估算結(jié)果

    夏季NPP總量占全年的80%左右,月均量為2964.91tC,平均密度達(dá)91.33gC m-2月-1;春季、秋季NPP總量分別占比11%和7%,冬季僅占1.3%;不同季節(jié)NPP標(biāo)準(zhǔn)差與平均密度呈正比(表2,圖6)。

    3.2 各島NPP分布特征

    各島NPP總量與其面積大小呈正比,面積最大的南長山島NPP總量達(dá)4127.8t C/a,占37.4%;面積最小的小黑山島為371.07t C/a,僅占3.4%;各島NPP平均密度存在差異,最大值(大黑山島)與最小值(小黑山島)之間差值達(dá)104g C m-2a-1;NPP最高點出現(xiàn)在大黑山島;各島NPP標(biāo)準(zhǔn)差存在差異,其與NPP平均密度、海島面積不存在明顯的關(guān)系(表3)。

    表2 不同季節(jié)南五島NPP估算結(jié)果Table 2 NPP estimation results of five southern islands in different seasons

    圖6 不同季節(jié)南五島NPP分布特征Fig.6 NPP distribution of five southern islands in different seasons

    表3 南五島各島NPP估算結(jié)果Table 3 NPP estimation results of each island

    3.3 不同地表覆蓋類型中NPP分布特征

    不同地表覆蓋中,針葉林NPP總量最大,占70%以上,裸地NPP總量最小,僅占0.15%;除了建設(shè)用地外,其余地表覆蓋類型NPP總量與面積大小呈正比。闊葉林的NPP平均密度最高,其次為針葉林、農(nóng)田、草地、建設(shè)用地和裸地,這與地表覆蓋類型的最大光能利用率有關(guān);NPP標(biāo)準(zhǔn)差也同平均密度大致呈正相關(guān)關(guān)系(表4)。

    表4 不同地表覆蓋類型NPP估算結(jié)果Table 4 Estimation results in different land cover type

    4 討論與結(jié)論

    4.1 討論

    4.1.1 NPP估算方法探討

    在以往基于遙感的NPP估算研究中,多采用1個時相的遙感影像提取NDVI值,并以此代表全年的情況,這在季節(jié)差異性明顯的地區(qū)顯然具有誤差。鑒于此,本文采用4個時相遙感影像以反映不同季節(jié)的實際情況,應(yīng)用于植被狀況隨季節(jié)變化明顯的北方典型海島地區(qū),能夠盡可能地反映地區(qū)實際情況。

    目前關(guān)于NPP的研究多為省域、全國等大尺度的研究[5-6,24,26-29],其光合有效輻射(PAR)基于不同區(qū)域太陽輻射值的差異從而擁有空間分異性,而在現(xiàn)有較小空間尺度的研究中,由于太陽輻射數(shù)據(jù)較為單一且缺乏空間分異性,PAR多作為常量來進(jìn)行計算,這在地形平坦的區(qū)域是可行的[7],但對于南五島這種具有明顯地勢起伏的區(qū)域而言顯然是不科學(xué)的。因此,本文嘗試探討地形對于PAR的影響,在區(qū)分太陽輻射中散射輻射和直接輻射的基礎(chǔ)上,借助GIS工具計算遮蔽度(Hillshade)。該方法可能具有一定的主觀性,但散射輻射和直接輻射的比例以及遮蔽度的計算方法均來源于較為成熟和長期的研究,以此得到地形影響因子作用下具有空間分異性的PAR數(shù)據(jù),是積極的嘗試,符合南五島的實際情況。遮蔽度的計算結(jié)果根據(jù)時間的不同而具有差異,本文依據(jù)格林尼治天文臺觀察的每日正午太陽高度資料,得到4個時相遙感影像獲取當(dāng)日的南五島正午太陽高度角,結(jié)合坡度和坡向數(shù)據(jù)計算得出遮蔽度,該遮蔽度實際上是正午的遮蔽度,以此代表全天的狀況可能具有一定誤差,但鑒于太陽高度數(shù)據(jù)的可獲得性,本文略有遺憾,在接下來的工作中應(yīng)對此開展進(jìn)一步的研究。

    本文NPP模型中的最大光能利用率(ξmax)來源于Running和朱文泉研究,其對南五島而言具有可行性,同時也擁有一定的主觀性,可能會給計算結(jié)果帶來誤差。在下一步的工作中應(yīng)結(jié)合現(xiàn)場監(jiān)測和驗證,對南五島植被最大光能利用率的研究進(jìn)行完善。

    表5 不同地區(qū)NPP估算結(jié)果Table 5 NPP estimation results in different areas

    4.1.2 NPP估算結(jié)果探討

    將本文計算的南五島NPP結(jié)果與采用相同方法的國內(nèi)其他研究結(jié)果進(jìn)行對比(表5)。與全國相比,南五島NPP平均密度和陶波[26]估算的342g C m-2a-1基本一致,略高于朱文泉[5]的324g C m-2a-1,低于顧娟[27]的393.75g C m-2a-1;與各地區(qū)相比,低于廣東省、江蘇省、盤錦市等沿海地區(qū),高于甘肅、青海、西藏。由此可得,南五島NPP平均密度處于全國平均水平,高于同緯度的西部地區(qū),但低于東部大陸沿海地區(qū)。

    本文計算結(jié)果顯示不同季節(jié)南五島NPP具有較大的差異,其中夏季的NPP總量是冬季的60倍以上,因此對于北方地區(qū)的NPP估算,不同季節(jié)的差異性研究尤為關(guān)鍵。南五島各島的NPP平均密度具有一定差異,結(jié)合各島地表覆蓋類型來看,NPP平均密度與各島建設(shè)用地比例(POC)呈明顯的負(fù)相關(guān),線性擬合度達(dá)0.468。城鎮(zhèn)建設(shè)占用土地,使得林地等具有高NPP值的覆蓋類型面積減少,進(jìn)而降低整島的NPP平均密度。這也某種程度上說明城鎮(zhèn)化進(jìn)程不可避免地會帶來海島NPP的減少。

    南五島均為基巖海島,以剝蝕山丘為主要地貌特征,存在水資源匱乏、土層薄等問題,本身植物生長條件較差,20世紀(jì)50年代以來進(jìn)行了廣泛的人工林種植,以黑松和刺槐為主要優(yōu)勢種,目前南五島林地面積占總面積的比例達(dá)55.8%,本文的結(jié)果顯示針葉林和闊葉林的NPP平均密度分別達(dá)446.9g C m-2a-1和559.71g C m-2a-1,略低于全國的平均數(shù)據(jù)469g C m-2a-1和663g C m-2a-1[5]。考慮到海島自身較脆弱的生態(tài)環(huán)境條件,說明黑松和刺槐具有良好的生命力,是南五島生態(tài)建設(shè)的理想物種。

    4.1.3 NPP對空間分異性參數(shù)的依賴性

    NPP與空間分異性參數(shù)相關(guān)性分析結(jié)果見表6。簡單相關(guān)分析顯示,NPP與NDVI、高程、坡度和地表覆蓋類型在0.01水平上顯著正相關(guān),與坡向不具有顯著的相關(guān)性,相關(guān)系數(shù)按由大到小的順序依次為NDVI、地表覆蓋類型、高程、坡度、坡向。偏相關(guān)分析顯示,除高程外各參數(shù)均與NPP在0.01水平上顯著相關(guān),而高程與NPP不具有顯著的相關(guān)性;除坡向外,各參數(shù)的偏相關(guān)系數(shù)同簡單相關(guān)系數(shù)相比均有不同程度的下降,尤其是高程和坡度下降明顯,且坡度由顯著正相關(guān)變?yōu)轱@著負(fù)相關(guān),同時坡向的相關(guān)性呈現(xiàn)明顯的上升。

    表6 NPP與模型參數(shù)相關(guān)分析結(jié)果Table 6 Results of correlation analysis between NPP and model parameters

    **.在0.01水平(雙側(cè))上顯著相關(guān);* 在 0.05水平(雙側(cè))上顯著相關(guān)

    圖7 NPP平均密度與建設(shè)用地比例關(guān)系 Fig.7 Relationship between mean density of NPP and proportion of construction land

    NDVI和地表覆蓋類型是NPP最主要的影響參數(shù),它們各自與NPP顯著正相關(guān),高程和坡度則通過影響NDVI和地表覆蓋類型進(jìn)而作用于NPP。隨著高程和坡度的上升,城鎮(zhèn)建設(shè)和農(nóng)田開墾程度減小,林地面積相應(yīng)增多,從而帶來較高的NPP;但偏相關(guān)分析結(jié)果顯示,在控制其他參數(shù)不變的情況下,高程與NPP關(guān)系不明顯,而坡度與NPP呈現(xiàn)顯著的負(fù)相關(guān),這是因為單從坡度因子來看,坡度較小的地方往往擁有更好的植被生長條件。坡向越接近正南,理論上能夠獲得更多的太陽輻射,從而擁有較高的NPP值,但不同坡向的地表覆蓋類型往往受到人類的影響而不具有規(guī)律性,因此在控制其他參數(shù)的環(huán)境下坡向與NPP具有顯著的正相關(guān)性,但在簡單相關(guān)分析中卻沒有呈現(xiàn)出該特征。

    4.1.4 NPP與土壤理化性質(zhì)的相關(guān)性

    NPP與10項土壤因子之間的簡單相關(guān)性見表7。結(jié)果顯示NPP與土壤pH、有效磷、全磷、全鉀呈顯著負(fù)相關(guān)關(guān)系,其中與前三項在0.01水平上顯著相關(guān),與全鉀在0.05水平上顯著相關(guān);NPP與全氮、總碳、總有機(jī)碳呈顯著正相關(guān)關(guān)系,與全氮在0.05水平上顯著相關(guān),與后二項在0.01水平上顯著相關(guān);NPP與含水量、速效鉀和含鹽量之間不存在顯著相關(guān)性。此外,10項土壤因子之間也存在著不同水平的相關(guān)關(guān)系。

    表7 NPP與土壤理化性質(zhì)的簡單相關(guān)分析結(jié)果Table 7 Results of simple correlation analysis between NPP and soil characters

    Ⅰ: pH;Ⅱ: 含水量Moisture content;Ⅲ: 有效磷Available phosphorus (AP);Ⅳ: 速效鉀Available potassium (AK);Ⅴ: 全磷Total phosphorus (TP);Ⅵ全鉀: Total potassium (TK);Ⅶ: 含鹽量Salinity;Ⅷ: 全氮Total nitrogen (TN);Ⅸ: 總碳Total carbon (TC);Ⅹ: 總有機(jī)碳Total organic carbon (TOC)

    土壤理化性質(zhì)與生態(tài)系統(tǒng)生產(chǎn)力密切相關(guān)[32]。具有較高生產(chǎn)力的群落凋落物往往較多[33],而凋落物是土壤有機(jī)質(zhì)的重要來源,有機(jī)質(zhì)中的主要元素為碳、氧、氫、氮,上述分析顯示了南五島NPP與土壤全氮、總碳、總有機(jī)碳呈明顯的正相關(guān)關(guān)系。南五島的土壤本身較為貧瘠,磷、鉀等營養(yǎng)元素較少,較高的生產(chǎn)力往往表示著群落能夠更加充分地吸收土壤中的營養(yǎng)元素而使其減少[34],這可能是本文NPP與有效磷、全磷、全鉀顯著負(fù)相關(guān)的原因。大多數(shù)土壤養(yǎng)分元素的有效性受土壤pH 值影響[35],南五島的NPP與pH呈現(xiàn)較高的負(fù)相關(guān),說明較低的pH值更適宜南五島植物的生長。土壤水是植被生長發(fā)育所需用水的主要來源,但過高的含水量會導(dǎo)致植被根系的腐爛不利于植被的生長;受到海水的作用,海島土壤的含鹽量往往過高,但南五島并沒有表現(xiàn)出該特征,說明南五島人工林對于過濾海鹽具有重要作用[12]。

    4.2 結(jié)論

    本文依據(jù)CASA模型,結(jié)合區(qū)域?qū)嶋H情況構(gòu)建廟島群島南五島NPP估算模型,以RS和GIS為主要手段,對研究區(qū)不同海島、不同季節(jié)、不同地表覆蓋類的NPP進(jìn)行估算,進(jìn)而分析NPP與模型計算參數(shù)和土壤理化性質(zhì)的關(guān)系。基本結(jié)論如下:

    (1)南五島NPP總量為11043.52t C/a,平均密度為340.19g C m-2a-1,總體上處于全國平均水平,高于同緯度的西部地區(qū),但低于東部沿海大陸地區(qū)。

    (2)南五島夏季NPP總量占全年的80%左右,春季、秋季NPP總量分別占比11%和7%,冬季僅占1.3%,說明研究區(qū)NPP具有明顯的季節(jié)差異。

    (3)南五島各島中,NPP平均密度由大到小依次為大黑山島、北長山島、廟島、南長山島和小黑山島,各島NPP平均密度與建設(shè)用地比例呈明顯的負(fù)相關(guān)。

    (4)不同地表覆蓋類型中,NPP平均密度由大到小依次為闊葉林、針葉林、農(nóng)田、草地、建設(shè)用地和裸地,林地具有較高的NPP,說明南五島的人工林建設(shè)具有重要的生態(tài)作用。

    (5)NDVI和地表覆蓋類型是NPP最主要的影響參數(shù),高程、坡向通過影響NDVI值和地表覆蓋類型間接作用于NPP。

    (6)南五島NPP與土壤pH、有效磷、全磷、全鉀呈顯著負(fù)相關(guān),與全氮、總碳、總有機(jī)碳呈顯著正相關(guān),與含水量、速效鉀和含鹽量之間相關(guān)關(guān)系不明顯。

    在下一步的研究中,應(yīng)重點關(guān)注地形影響因子的修正完善和植被最大光能利用率的現(xiàn)場驗證,以期為南五島碳源/匯研究提供更加充分的技術(shù)支持。

    [1] Lieth H, Whittaker R H.Primary Productivity of the Biosphere.New York: Springer-Verlag, 1975.

    [2] Field C B, Behrenfeld M J, Randerson J T, Falkowski P.Primary production of the biosphere: integrating terrestrial and oceanic components.Science, 1998, 281(5374): 237-240.

    [3] Paruelo J M, Epstei H E, Lauenroth W K, Burke I C.ANPP estimates from NDVI for the central grassland region of the United States.Ecology, 1997, 78(3): 953-958.

    [4] Running S W, Thornton P E, Nemani R, Glassy J M.Global terrestrial gross and net primary productivity from the earth observing system // Sala O E, Jackson R B, Mooney H A, Howarth R W.Methods in Ecosystem Science.New York: Springer-Verlag, 2000: 44-57.

    [5] 朱文泉, 潘耀忠, 張錦水.中國陸地植被凈初級生產(chǎn)力遙感估算.植物生態(tài)學(xué)報, 2007, 31(3): 413-424.

    [6] 龍慧靈, 李曉兵, 王宏, 魏丹丹, 張程.內(nèi)蒙古草原區(qū)植被凈初級生產(chǎn)力及其與氣候的關(guān)系.生態(tài)學(xué)報, 2010, 30(5): 1367-1378.

    [7] 王莉雯, 衛(wèi)亞星.盤錦濕地凈初級生產(chǎn)力時空分布特征.生態(tài)學(xué)報, 2012, 32(19): 6006-6015.

    [8] 石洪華, 鄭偉, 丁德文, 呂吉斌.典型海島生態(tài)系統(tǒng)服務(wù)及價值評估.海洋環(huán)境科學(xué), 2009, 28(6): 743-748.

    [9] 王曉麗, 王嬡, 石洪華, 鄭偉, 周然.海島陸地生態(tài)系統(tǒng)固碳估算方法.生態(tài)學(xué)報, 2014, 34(1): 88-96.

    [10] 王曉麗, 王嬡, 石洪華, 彭士濤, 宮立鵬, 覃雪波.山東省長島縣南長山島黑松和刺槐人工林的碳儲量.應(yīng)用生態(tài)學(xué)報, 2013, 24(5): 1263-1268.

    [11] Paulay G.Biodiversity on oceanic islands: its origin and extinction.American Zoology, 1994, 34(1): 134-144.

    [12] 石洪華, 王曉麗, 王嬡, 劉振英, 麻德明.北長山島森林喬木層碳儲量及其影響因子.生態(tài)學(xué)報, 2013, 33(19): 6363-6372.

    [13] Potter C S, Randerson J T, Field C B, Matson P A, Vitousek P M, Mooney H A, Klooster S A.Terrestrial ecosystem production: A process model based on global satellite and surface data.Global Biogeochemical Cycles, 1993, 7(4): 811-841.

    [14] 袁淑杰, 繆啟龍, 谷曉平, 邱新法.貴州高原起伏地形下太陽直接輻射的精細(xì)分布.自然資源學(xué)報, 2009, 24(8): 1432-1439.

    [15] 左大康, 弓冉.中國太陽直接輻射、散射輻射和太陽總輻射間的關(guān)系.地理學(xué)報, 1962, 28(3): 175-186.

    [16] 馬金玉, 梁宏, 羅勇, 李世奎.中國近50年太陽直接輻射和散射輻射變化趨勢特征.物理學(xué)報, 2011, 60(6): 069601.

    [17] Burrough P A, McDonell R A.Principles of Geographical Information Systems.2nd ed.New York: Oxford University Press, 1998: 190-190.

    [18] Ruimy A, Saugier B, Dedieu G.Methodology for the estimation of terrestrial net primary production from remotely sensed data.Journal of Geophysical Research, 1994, 99(D3): 5263-5283.

    [19] Field C B, Randerson J T, Malmstr?m C M.Global net primary production: combining ecology and remote sensing.Remote Sensing Environment, 1995, 51(1): 74-88.

    [20] 樸世龍, 方精云, 郭慶華.利用CASA模型估算我國植被凈第一性生產(chǎn)力.植物生態(tài)學(xué)報, 2001, 25(5): 603-608.

    [21] 周廣勝, 張新時.自然植被凈第一性生產(chǎn)力模型初探.植物生態(tài)學(xué)報, 1995, 19(3): 193-200.

    [22] 《噴灌工程設(shè)計手冊》編寫組.噴灌工程設(shè)計手冊.北京: 水利電力出版社, 1989: 348-356.

    [23] 周廣勝, 張新時.全球變化的中國氣候——植被分類研究.植物學(xué)報, 1996, 38(1): 8-17.

    [24] 張鐿鋰, 祁威, 周才平, 丁明軍, 劉林山, 高俊剛, 擺萬奇, 王兆鋒, 鄭度.青藏高原高寒草地凈初級生產(chǎn)力(NPP)時空分異.地理學(xué)報, 2013, 68(9): 1197-1211.

    [25] 穆少杰, 李建龍, 周偉, 楊紅飛, 章超斌, 居為民.2001-2010年內(nèi)蒙古植被凈初級生產(chǎn)力的時空格局及其與氣候的關(guān)系.生態(tài)學(xué)報, 2013, 33(12): 3752-3764.

    [26] 陶波, 李克讓, 邵雪梅, 曹明奎.中國陸地凈初級生產(chǎn)力時空特征模擬.地理學(xué)報, 2003, 58(3): 372-380.

    [27] 顧娟, 李新, 黃春林, 張喜風(fēng), 金鑫.2002—2010年中國陸域植被凈初級生產(chǎn)力模擬.蘭州大學(xué)學(xué)報: 自然科學(xué)版, 2013, 49(2): 203-213.

    [28] 羅艷, 王春林.基于MODIS NDVI的廣東省陸地生態(tài)系統(tǒng)凈初級生產(chǎn)力估算.生態(tài)環(huán)境學(xué)報, 2009, 18(4): 1467-1471.

    [29] 王駟鷂, 劉振波.江蘇省植被凈初級生產(chǎn)力時空分布格局研究.南京信息工程大學(xué)學(xué)報: 自然科學(xué)版, 2012, 4(4): 321-325.

    [30] 劉春雨, 董曉峰, 劉英英, 潘竟虎, 車彥軍.甘肅省凈初級生產(chǎn)力時空變化特征.中國人口·資源與環(huán)境, 2014, 24(1): 163-170.

    [31] 衛(wèi)亞星, 王莉雯, 石迎春, 李路子.青海省草地資源凈初級生產(chǎn)力遙感監(jiān)測.地理科學(xué), 2012, 32(5): 621-627.

    [32] 王長庭, 龍瑞軍, 王啟基, 景增春, 尚占環(huán), 丁路明.高寒草甸不同海拔梯度土壤有機(jī)質(zhì)氮磷的分布和生產(chǎn)力變化及其與環(huán)境因子的關(guān)系.草業(yè)學(xué)報, 2005, 14(4): 15-20.

    [33] 馮宗煒, 王效科, 吳剛.中國森林生態(tài)系統(tǒng)的生物量和生產(chǎn)力.北京: 科學(xué)出版社, 1999.

    [34] Elton C S.The Ecology of Invasions by Animals and Plants.Chicago: University of Chicago Press, 2000.

    [35] Khattak R A, Hussain Z.Evaluation of soil fertility status and nutrition of orchards.Soil &Environment, 2007, 26(1): 22-32.

    The spatial distribution and impact factors of net primary productivity in the island ecosystem of five southern islands of Miaodao Archipelago

    CHI Yuan1, SHI Honghua1,*, WANG Xiaoli2, LI Jie3, FENG Aiping1

    1TheFirstInstituteofOceanography,StateOceanicAdministration,Qingdao266061,China2CollegeofEnvironmentalScienceandSafetyEngineering,TianjinUniversityofTechnology,Tianjin300384,China3SchoolofEnvironmentalandMunicipalEngineering,QingdaoTechnologicalUniversity,Qingdao266033,China

    Island ecosystems are characterized by self-containment and vulnerability due to their limited and isolated areas.They are also important for biodiversity maintenance, freshwater conservation, and other basic ecological functions.Moreover, island ecosystems play important roles in global carbon cycle, which makes the estimation of net primary productivity (NPP) more significant for ecological construction and the study of carbon source and sink on islands.The rapid development of RS and GIS technologies provides a practical method for NPP estimation.In this study, five southern islands of Miaodao Archipelago were examined and an NPP estimation model was established based on the Carnegie-Ames-Stanford Approach (CASA) and regional characteristics.To characterize each season, 4temporal remote sensing images were utilized.The normalized difference vegetation index (NDVI) and land cover types were extracted, and these constituted the data used for the NPP estimation model, along with meteorological and topographical data.The NPP values for different seasons, islands, and land cover types were calculated using GIS methods.The linear correlations between NPP values and various parameters including NDVI, altitude, slope, aspect, and land cover types were obtained to determine the importance of each parameter.Finally, the relationships between 10soil environment-chemical properties and NPP at 36sampling sites were analyzed.The total NPP for the five islands was 11043.52t C/a.The mean density of NPP (MDN) was 340.19g C m-2a-1, which was higher than that of western regions at a similar latitude (e.g., Gansu Province and Qinghai Province) and lower than that of coastal zones in the eastern region (e.g., Guangdong Province, Jiangsu Province, and Panjin City in Liaoning Province).NPP in the summer accounted for approximately 80% of the total yearly NPP, followed by spring (11%), autumn (7%), and winter (1.3%), indicating that NPP in the northern area differed greatly among seasons.The total NPP was highest for Nanchangshan Island, followed by Beichangshan Island, Daheishan Island, Miao Island, and Xiaoheishan Island, and the total NPP was positively correlated with island area.MDN was highest for Daheishan Island, followed by Beichangshan Island, Miao Island, Nanchangshan Island, and Xiaoheishan Island, and was negatively related to the proportion of construction land.The land cover types included broad-leaf forest, coniferous forest, farmland, grassland, construction land, and bare land, in descending order based on MDN.The high NPP for forest land indicated that the plantations on the five southern islands were of great ecological significance.NDVI and land cover types were the most important parameters with respect to the spatial distribution of NPP, and topographical parameters indirectly affected NPP via their effects on NDVI and land cover type.NPP had a significant negative relationship with soil pH (P<0.01), available phosphorus (P<0.01), total phosphorus (P<0.01), and total potassium (P<0.05), and a significant positive relationship with total carbon (P<0.01), total organic carbon (P<0.01) and total nitrogen (P<0.05).No significant relationships with water content, available potassium, or salinity were observed.

    island ecosystem;net primary productivity (NPP);five southern islands of Miaodao Archipelago;spatial distribution;impact factors

    科技基礎(chǔ)性工作專項(2012FY112500);海洋公益性行業(yè)科研專項(201305009, 201505012)

    2014-05-04; < class="emphasis_bold">網(wǎng)絡(luò)出版日期:

    日期:2015-05-19

    10.5846/stxb201405040880

    *通訊作者Corresponding author.E-mail: shihonghua@fio.org.cn

    池源,石洪華,王曉麗,李捷,豐愛平.廟島群島南五島生態(tài)系統(tǒng)凈初級生產(chǎn)力空間分布及其影響因子.生態(tài)學(xué)報,2015,35(24):8094-8106.

    Chi Y, Shi H H, Wang X L, Li J, Feng A P.The spatial distribution and impact factors of net primary productivity in the island ecosystem of five southern islands of Miaodao Archipelago.Acta Ecologica Sinica,2015,35(24):8094-8106.

    猜你喜歡
    坡向海島生產(chǎn)力
    冰與火共存的海島
    奧秘(2020年5期)2020-06-30 10:12:10
    DEM地表坡向變率的向量幾何計算法
    在海島度假
    把“風(fēng)景”轉(zhuǎn)化成“生產(chǎn)力
    紅土地(2016年3期)2017-01-15 13:45:17
    青藏高原東緣高寒草甸坡向梯度上植物光合生理特征研究
    土壤與作物(2015年3期)2015-12-08 00:46:55
    人口紅利與提高生產(chǎn)力
    海峽姐妹(2015年9期)2015-02-27 15:13:07
    展望2014:推動科技創(chuàng)新轉(zhuǎn)化為第一生產(chǎn)力
    杭州科技(2014年3期)2014-02-27 15:26:46
    凡事就怕做到極致 一個詞的生產(chǎn)力
    坡向和坡位對小流域梯田土壤有機(jī)碳、氮變化的影響
    神奇的海島
    看十八女毛片水多多多| 色网站视频免费| 97在线人人人人妻| 69人妻影院| 搡女人真爽免费视频火全软件| www.av在线官网国产| 亚洲成人一二三区av| 在线观看免费高清a一片| 99热国产这里只有精品6| 国产探花极品一区二区| 下体分泌物呈黄色| 日韩大片免费观看网站| 国产爱豆传媒在线观看| 成人午夜精彩视频在线观看| 99热全是精品| 69av精品久久久久久| 国产精品秋霞免费鲁丝片| 黄色日韩在线| 日韩电影二区| 久久精品久久久久久久性| 久久99热这里只频精品6学生| 久久久成人免费电影| 欧美三级亚洲精品| 真实男女啪啪啪动态图| 免费观看性生交大片5| 蜜臀久久99精品久久宅男| 国产淫语在线视频| 特大巨黑吊av在线直播| 久久精品久久精品一区二区三区| 在线观看一区二区三区| 成年女人在线观看亚洲视频 | 国产色婷婷99| 精品久久久久久久末码| 草草在线视频免费看| 欧美激情久久久久久爽电影| 亚洲精品一区蜜桃| 亚洲欧美清纯卡通| 99久久精品热视频| 九九在线视频观看精品| 精品亚洲乱码少妇综合久久| 色婷婷久久久亚洲欧美| 一本色道久久久久久精品综合| 国内少妇人妻偷人精品xxx网站| 天堂网av新在线| 日本一二三区视频观看| 国国产精品蜜臀av免费| 亚洲精品乱码久久久久久按摩| 高清视频免费观看一区二区| 日日啪夜夜爽| 亚洲欧美日韩另类电影网站 | 看黄色毛片网站| 一级av片app| 免费观看的影片在线观看| 国产男女内射视频| 久久久亚洲精品成人影院| 久久久久久久久久成人| 亚洲图色成人| 婷婷色综合www| 汤姆久久久久久久影院中文字幕| 成人二区视频| 一区二区av电影网| 久久久久久久久久久丰满| 日韩制服骚丝袜av| 中文在线观看免费www的网站| 成人综合一区亚洲| 高清av免费在线| 看十八女毛片水多多多| 日韩视频在线欧美| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久精品古装| 好男人视频免费观看在线| 国产精品三级大全| 日日撸夜夜添| 国产精品成人在线| 高清午夜精品一区二区三区| 少妇的逼水好多| 国产色爽女视频免费观看| 美女高潮的动态| 国语对白做爰xxxⅹ性视频网站| 日韩av免费高清视频| eeuss影院久久| av在线蜜桃| 女人久久www免费人成看片| 国产伦精品一区二区三区四那| 精品人妻一区二区三区麻豆| 久久久久久久亚洲中文字幕| 日韩国内少妇激情av| 身体一侧抽搐| 高清日韩中文字幕在线| 午夜免费鲁丝| 久久久久久久精品精品| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美 国产精品| 一个人看的www免费观看视频| 一级毛片aaaaaa免费看小| 欧美3d第一页| 国产精品国产三级国产专区5o| 亚洲成人久久爱视频| 国产精品一区www在线观看| 精品久久久精品久久久| 久久韩国三级中文字幕| 久久久久久国产a免费观看| xxx大片免费视频| 国产精品av视频在线免费观看| 精品久久久久久久久亚洲| 亚洲精品乱码久久久久久按摩| 亚洲第一区二区三区不卡| 全区人妻精品视频| 日韩av在线免费看完整版不卡| 人妻一区二区av| 国产精品久久久久久精品古装| 久久久久网色| av国产久精品久网站免费入址| 国内揄拍国产精品人妻在线| 日本黄大片高清| 内射极品少妇av片p| 欧美日本视频| 精品视频人人做人人爽| 国产一区有黄有色的免费视频| 丝瓜视频免费看黄片| 一级av片app| 欧美丝袜亚洲另类| 99久国产av精品国产电影| 精品久久久精品久久久| 女人十人毛片免费观看3o分钟| 激情五月婷婷亚洲| 国产成年人精品一区二区| 最新中文字幕久久久久| 国产极品天堂在线| 身体一侧抽搐| 高清午夜精品一区二区三区| www.av在线官网国产| 人人妻人人爽人人添夜夜欢视频 | 久久精品熟女亚洲av麻豆精品| 亚洲天堂国产精品一区在线| 如何舔出高潮| 亚洲欧美日韩另类电影网站 | 26uuu在线亚洲综合色| 国产永久视频网站| 亚洲欧美精品自产自拍| 在线a可以看的网站| av天堂中文字幕网| 九草在线视频观看| 麻豆成人av视频| 日韩在线高清观看一区二区三区| 亚洲在线观看片| 丝袜喷水一区| 国产精品蜜桃在线观看| 99久久精品国产国产毛片| 女人十人毛片免费观看3o分钟| 久久99蜜桃精品久久| 男女边摸边吃奶| 国内精品美女久久久久久| 免费av不卡在线播放| 国产探花极品一区二区| 国产欧美日韩精品一区二区| 丝袜喷水一区| 亚洲av在线观看美女高潮| 亚洲美女搞黄在线观看| 直男gayav资源| 国产爽快片一区二区三区| 香蕉精品网在线| 蜜桃亚洲精品一区二区三区| 毛片女人毛片| 欧美成人a在线观看| 69人妻影院| 国产黄片视频在线免费观看| 女人十人毛片免费观看3o分钟| 久久6这里有精品| 亚洲av一区综合| 国产黄片视频在线免费观看| 啦啦啦中文免费视频观看日本| 久久99热这里只频精品6学生| 综合色丁香网| 午夜亚洲福利在线播放| 亚洲av一区综合| 中国国产av一级| 国产亚洲最大av| 日韩av免费高清视频| 黑人高潮一二区| 国产精品人妻久久久影院| 久久人人爽人人爽人人片va| 日韩三级伦理在线观看| 爱豆传媒免费全集在线观看| 99久久中文字幕三级久久日本| 草草在线视频免费看| 3wmmmm亚洲av在线观看| 亚洲欧美精品专区久久| 91午夜精品亚洲一区二区三区| 亚洲av欧美aⅴ国产| 国产精品久久久久久久电影| 国产精品一区二区在线观看99| 熟女电影av网| 精品熟女少妇av免费看| 在线a可以看的网站| 九九在线视频观看精品| 成人亚洲欧美一区二区av| 国产男女超爽视频在线观看| 你懂的网址亚洲精品在线观看| 两个人的视频大全免费| 久久精品人妻少妇| 精品久久久久久电影网| 在线观看一区二区三区| 少妇人妻 视频| 国产精品人妻久久久影院| 爱豆传媒免费全集在线观看| 丰满人妻一区二区三区视频av| 日日撸夜夜添| 亚洲精品色激情综合| 亚洲国产欧美在线一区| 国产亚洲午夜精品一区二区久久 | 一级黄片播放器| 中文乱码字字幕精品一区二区三区| 热99国产精品久久久久久7| 国产人妻一区二区三区在| 最后的刺客免费高清国语| 丰满少妇做爰视频| 亚洲高清免费不卡视频| 日日摸夜夜添夜夜爱| 青青草视频在线视频观看| 一级二级三级毛片免费看| 嫩草影院精品99| 欧美一区二区亚洲| av在线天堂中文字幕| 国产精品三级大全| www.色视频.com| av播播在线观看一区| 国产真实伦视频高清在线观看| freevideosex欧美| 久久久久久九九精品二区国产| 一级毛片黄色毛片免费观看视频| 看十八女毛片水多多多| 岛国毛片在线播放| 99re6热这里在线精品视频| 国产毛片在线视频| 黄色配什么色好看| 国产一区二区三区综合在线观看 | 一级毛片黄色毛片免费观看视频| 亚洲色图av天堂| 亚洲在久久综合| 国产精品99久久久久久久久| 欧美日本视频| 在线观看一区二区三区激情| 伊人久久国产一区二区| 亚洲人成网站高清观看| 亚洲va在线va天堂va国产| 男人和女人高潮做爰伦理| 国产在线一区二区三区精| 日韩,欧美,国产一区二区三区| av在线观看视频网站免费| 日韩伦理黄色片| 激情 狠狠 欧美| 国产一级毛片在线| 日韩在线高清观看一区二区三区| av在线观看视频网站免费| 久久人人爽av亚洲精品天堂 | 中文字幕亚洲精品专区| 精品少妇黑人巨大在线播放| 少妇人妻一区二区三区视频| 国产高清国产精品国产三级 | 晚上一个人看的免费电影| 大陆偷拍与自拍| 听说在线观看完整版免费高清| 日本一本二区三区精品| 搡女人真爽免费视频火全软件| 最近的中文字幕免费完整| 久久女婷五月综合色啪小说 | 联通29元200g的流量卡| 七月丁香在线播放| 高清毛片免费看| 亚洲av国产av综合av卡| 好男人在线观看高清免费视频| 亚洲欧美日韩东京热| 欧美老熟妇乱子伦牲交| 精品人妻偷拍中文字幕| 99久久精品国产国产毛片| 欧美人与善性xxx| 又粗又硬又长又爽又黄的视频| 别揉我奶头 嗯啊视频| 亚洲欧美精品专区久久| 三级国产精品片| 国产成人精品福利久久| 色吧在线观看| 国产91av在线免费观看| 91aial.com中文字幕在线观看| 国产精品久久久久久精品电影小说 | 新久久久久国产一级毛片| 午夜福利在线在线| 国产 精品1| 99精国产麻豆久久婷婷| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美成人精品一区二区| 2018国产大陆天天弄谢| 国产高清三级在线| 国产成人91sexporn| 成人鲁丝片一二三区免费| 22中文网久久字幕| 男人爽女人下面视频在线观看| 岛国毛片在线播放| 午夜亚洲福利在线播放| 亚洲精华国产精华液的使用体验| 精品久久久久久久久av| 在线观看一区二区三区激情| 久久久久久伊人网av| 国产精品伦人一区二区| 蜜臀久久99精品久久宅男| 日本欧美国产在线视频| 国产一区二区三区综合在线观看 | 免费黄色在线免费观看| 全区人妻精品视频| 欧美xxxx黑人xx丫x性爽| 国产探花极品一区二区| 久久久久久久精品精品| 亚洲国产精品999| 久久久国产一区二区| 国产乱人偷精品视频| 深夜a级毛片| 国产精品99久久99久久久不卡 | 国产黄a三级三级三级人| 亚洲精品第二区| 亚洲欧美一区二区三区黑人 | 亚洲成人精品中文字幕电影| 最近2019中文字幕mv第一页| 日韩欧美 国产精品| 最后的刺客免费高清国语| 大香蕉久久网| 国产av国产精品国产| 日韩成人伦理影院| 国产高清国产精品国产三级 | 99re6热这里在线精品视频| 亚洲自偷自拍三级| 99久久人妻综合| 国产精品无大码| 精品一区二区三区视频在线| 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人综合色| 国产成人午夜福利电影在线观看| 中文精品一卡2卡3卡4更新| 午夜精品国产一区二区电影 | 亚洲欧美日韩卡通动漫| 亚洲欧美成人精品一区二区| 亚洲国产欧美人成| 成人一区二区视频在线观看| 国产一区亚洲一区在线观看| av在线app专区| 久久99精品国语久久久| 美女内射精品一级片tv| 三级经典国产精品| 久久6这里有精品| 久久精品国产亚洲av涩爱| 亚洲自偷自拍三级| 丝袜喷水一区| 亚洲av男天堂| 99久久人妻综合| 少妇丰满av| 女人被狂操c到高潮| 午夜福利高清视频| 99久国产av精品国产电影| 成人免费观看视频高清| 午夜免费鲁丝| 搞女人的毛片| 中国国产av一级| 精品久久久噜噜| 国产精品av视频在线免费观看| 国产69精品久久久久777片| 中文天堂在线官网| 黄色视频在线播放观看不卡| 人妻夜夜爽99麻豆av| 国语对白做爰xxxⅹ性视频网站| 久久久久国产网址| 少妇裸体淫交视频免费看高清| 久久久精品免费免费高清| 国产av不卡久久| 国产综合懂色| 亚洲欧美成人综合另类久久久| 建设人人有责人人尽责人人享有的 | 新久久久久国产一级毛片| 黄片wwwwww| 久热久热在线精品观看| 大话2 男鬼变身卡| 亚洲欧美清纯卡通| 精品国产露脸久久av麻豆| 欧美zozozo另类| 国产伦理片在线播放av一区| 亚洲精品国产色婷婷电影| 久久精品熟女亚洲av麻豆精品| 精品一区二区免费观看| 直男gayav资源| 日韩中字成人| 晚上一个人看的免费电影| av播播在线观看一区| 水蜜桃什么品种好| 国产爱豆传媒在线观看| 国产爱豆传媒在线观看| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 男人和女人高潮做爰伦理| av天堂中文字幕网| 91精品伊人久久大香线蕉| 不卡视频在线观看欧美| av一本久久久久| 国产精品人妻久久久久久| 激情 狠狠 欧美| av专区在线播放| 一区二区av电影网| 亚洲av在线观看美女高潮| 亚洲av成人精品一二三区| 丝瓜视频免费看黄片| 一区二区av电影网| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 亚洲av成人精品一区久久| 九九久久精品国产亚洲av麻豆| 亚洲av在线观看美女高潮| 国产精品一及| 亚洲在线观看片| 中文精品一卡2卡3卡4更新| 欧美97在线视频| 日韩av免费高清视频| 国产午夜精品一二区理论片| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 国产高清有码在线观看视频| 久久久久久久久久久免费av| 国产在视频线精品| 欧美精品人与动牲交sv欧美| 狂野欧美激情性bbbbbb| 自拍偷自拍亚洲精品老妇| 成人二区视频| 亚洲美女视频黄频| 午夜精品一区二区三区免费看| 亚洲av福利一区| 免费人成在线观看视频色| 色哟哟·www| 久久久欧美国产精品| 亚洲av一区综合| 免费看光身美女| 久久精品国产鲁丝片午夜精品| 亚洲精品日韩av片在线观看| 亚洲人成网站在线播| 久久久久精品久久久久真实原创| 国产大屁股一区二区在线视频| 久久人人爽人人爽人人片va| 欧美变态另类bdsm刘玥| 少妇裸体淫交视频免费看高清| 精品酒店卫生间| 国模一区二区三区四区视频| 寂寞人妻少妇视频99o| 国产精品一区www在线观看| 日韩av在线免费看完整版不卡| 色婷婷久久久亚洲欧美| 午夜福利网站1000一区二区三区| 99久久精品国产国产毛片| av.在线天堂| 熟女人妻精品中文字幕| 自拍偷自拍亚洲精品老妇| 欧美成人精品欧美一级黄| 国产一区二区三区av在线| 精品99又大又爽又粗少妇毛片| 国产成人免费观看mmmm| av线在线观看网站| 黑人高潮一二区| 亚洲在线观看片| 亚洲精品国产色婷婷电影| 亚洲性久久影院| 在线观看美女被高潮喷水网站| 岛国毛片在线播放| 毛片女人毛片| 亚洲精品自拍成人| 久久精品国产亚洲网站| 亚洲真实伦在线观看| 一级黄片播放器| 日本色播在线视频| 婷婷色av中文字幕| 少妇丰满av| 亚洲国产精品国产精品| 一级a做视频免费观看| 青春草视频在线免费观看| 激情五月婷婷亚洲| 白带黄色成豆腐渣| 日本一二三区视频观看| 日韩av不卡免费在线播放| 美女高潮的动态| 看非洲黑人一级黄片| 老女人水多毛片| 亚洲av福利一区| 久久久欧美国产精品| 深爱激情五月婷婷| 国产男女超爽视频在线观看| 国产精品99久久99久久久不卡 | 久久久久久久大尺度免费视频| 国产成人91sexporn| 亚洲在线观看片| av一本久久久久| 久久久久国产网址| 好男人在线观看高清免费视频| 亚洲精品视频女| 香蕉精品网在线| 久久99热这里只有精品18| 天堂网av新在线| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 能在线免费看毛片的网站| 亚洲精品aⅴ在线观看| 国产亚洲一区二区精品| 国产成人91sexporn| 精品视频人人做人人爽| 在线观看美女被高潮喷水网站| 国产黄片美女视频| 亚洲国产色片| 亚洲精品成人久久久久久| av在线app专区| 久久精品综合一区二区三区| 岛国毛片在线播放| 麻豆成人午夜福利视频| 亚洲自拍偷在线| 亚洲国产精品国产精品| 国产亚洲最大av| 伦理电影大哥的女人| 毛片女人毛片| 黄片无遮挡物在线观看| 中文在线观看免费www的网站| 国产视频首页在线观看| 看十八女毛片水多多多| 一本一本综合久久| 亚洲自拍偷在线| 乱码一卡2卡4卡精品| 97在线视频观看| 久久鲁丝午夜福利片| 毛片女人毛片| 亚洲熟女精品中文字幕| 中文字幕免费在线视频6| 日韩制服骚丝袜av| 欧美变态另类bdsm刘玥| 一二三四中文在线观看免费高清| 成人毛片60女人毛片免费| 国产乱来视频区| 一区二区av电影网| 看黄色毛片网站| 少妇猛男粗大的猛烈进出视频 | 久久99热6这里只有精品| 精品久久久久久电影网| 另类亚洲欧美激情| 毛片一级片免费看久久久久| 日韩欧美 国产精品| 晚上一个人看的免费电影| 久久久国产一区二区| 特级一级黄色大片| 亚洲欧美清纯卡通| 日本与韩国留学比较| 亚洲人成网站高清观看| 久久人人爽av亚洲精品天堂 | 亚洲国产欧美在线一区| 欧美+日韩+精品| 亚洲综合精品二区| 日产精品乱码卡一卡2卡三| 新久久久久国产一级毛片| av一本久久久久| 欧美最新免费一区二区三区| 久久久久精品性色| 男女下面进入的视频免费午夜| 久久精品熟女亚洲av麻豆精品| 久久99热这里只频精品6学生| 成年女人看的毛片在线观看| 婷婷色av中文字幕| 2018国产大陆天天弄谢| 国产毛片a区久久久久| 中文乱码字字幕精品一区二区三区| 日本色播在线视频| 可以在线观看毛片的网站| 亚洲精品久久午夜乱码| 寂寞人妻少妇视频99o| 国产一级毛片在线| 精品久久久久久久久av| 国产白丝娇喘喷水9色精品| 春色校园在线视频观看| 国产真实伦视频高清在线观看| 国产高清三级在线| 校园人妻丝袜中文字幕| av.在线天堂| 国产一级毛片在线| 国产成年人精品一区二区| 特级一级黄色大片| 大香蕉97超碰在线| 少妇被粗大猛烈的视频| 可以在线观看毛片的网站| 嫩草影院新地址| 三级国产精品欧美在线观看| 熟妇人妻不卡中文字幕| av一本久久久久| 久久久久久久久久人人人人人人| 亚洲精品成人av观看孕妇| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 少妇丰满av| 国产伦精品一区二区三区四那| av卡一久久| 两个人的视频大全免费| 成人鲁丝片一二三区免费| 777米奇影视久久| 黄色配什么色好看| 97精品久久久久久久久久精品| 免费不卡的大黄色大毛片视频在线观看| 欧美激情国产日韩精品一区| 免费av毛片视频| 亚洲在线观看片| 一二三四中文在线观看免费高清| 中文字幕av成人在线电影| 国产在线一区二区三区精| 久久综合国产亚洲精品| 三级国产精品欧美在线观看| 美女视频免费永久观看网站| 国产男人的电影天堂91| 老司机影院毛片| 黄色视频在线播放观看不卡| 少妇丰满av|