• <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ī)碳、氮變化的影響
    神奇的海島
    免费黄网站久久成人精品 | 免费看a级黄色片| 搡女人真爽免费视频火全软件 | 国产精品av视频在线免费观看| 波多野结衣高清作品| 亚洲aⅴ乱码一区二区在线播放| 免费在线观看成人毛片| 女人被狂操c到高潮| 久久久国产成人精品二区| 如何舔出高潮| 51午夜福利影视在线观看| 日韩欧美一区二区三区在线观看| 内射极品少妇av片p| 国产午夜精品久久久久久一区二区三区 | 国产精华一区二区三区| 亚洲熟妇熟女久久| 国产精品免费一区二区三区在线| 在线观看一区二区三区| 一级黄色大片毛片| 国产成人a区在线观看| 两人在一起打扑克的视频| 18禁裸乳无遮挡免费网站照片| 在线观看美女被高潮喷水网站 | 欧美bdsm另类| 日韩中字成人| 国产精品一及| a级毛片免费高清观看在线播放| 99国产综合亚洲精品| 观看美女的网站| 亚洲精品影视一区二区三区av| 国产淫片久久久久久久久 | 很黄的视频免费| 精品午夜福利视频在线观看一区| 国产在线精品亚洲第一网站| 搞女人的毛片| 国产高潮美女av| 日韩精品青青久久久久久| 久久久久久九九精品二区国产| 久久人人精品亚洲av| 久久久久亚洲av毛片大全| 日韩中字成人| 色哟哟·www| av女优亚洲男人天堂| 免费av不卡在线播放| 老熟妇仑乱视频hdxx| 久久久成人免费电影| 91久久精品国产一区二区成人| 亚洲精品456在线播放app | 12—13女人毛片做爰片一| 亚洲熟妇中文字幕五十中出| 日韩 亚洲 欧美在线| 99热精品在线国产| 日韩精品青青久久久久久| 90打野战视频偷拍视频| 亚洲一区高清亚洲精品| 激情在线观看视频在线高清| 男女做爰动态图高潮gif福利片| av视频在线观看入口| 成年女人毛片免费观看观看9| 好看av亚洲va欧美ⅴa在| av在线观看视频网站免费| 淫秽高清视频在线观看| 老司机午夜十八禁免费视频| 伊人久久精品亚洲午夜| 在线观看一区二区三区| 久久九九热精品免费| 久久国产精品人妻蜜桃| 动漫黄色视频在线观看| 99热只有精品国产| 一进一出抽搐动态| 国产成+人综合+亚洲专区| 国产男靠女视频免费网站| 国产精品女同一区二区软件 | 国产精品99久久久久久久久| 一区福利在线观看| 亚洲黑人精品在线| 国产精品久久电影中文字幕| 国产亚洲精品久久久久久毛片| 香蕉av资源在线| www.www免费av| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 丰满人妻一区二区三区视频av| 十八禁人妻一区二区| 亚洲天堂国产精品一区在线| 老熟妇仑乱视频hdxx| 色综合站精品国产| 免费在线观看亚洲国产| 欧美激情久久久久久爽电影| 美女黄网站色视频| 美女 人体艺术 gogo| 欧美高清成人免费视频www| 午夜免费成人在线视频| 中亚洲国语对白在线视频| 尤物成人国产欧美一区二区三区| 亚洲国产高清在线一区二区三| 亚洲av美国av| 欧美色视频一区免费| 99精品在免费线老司机午夜| 成人国产一区最新在线观看| 一夜夜www| 午夜福利成人在线免费观看| 一区二区三区激情视频| 国产高清视频在线播放一区| 黄色女人牲交| 亚洲国产色片| 久久久久久久久久成人| 久久这里只有精品中国| 久久久久九九精品影院| 波多野结衣巨乳人妻| 亚洲美女黄片视频| 岛国在线免费视频观看| av在线老鸭窝| ponron亚洲| 国产精品精品国产色婷婷| 国产精品亚洲av一区麻豆| 一个人看视频在线观看www免费| av国产免费在线观看| 毛片一级片免费看久久久久 | 国产在线男女| 欧美日韩黄片免| АⅤ资源中文在线天堂| 精品无人区乱码1区二区| 国产精品人妻久久久久久| 人人妻人人澡欧美一区二区| 毛片一级片免费看久久久久 | 国产精品不卡视频一区二区 | 国产 一区 欧美 日韩| 日本黄大片高清| 色哟哟·www| 免费人成视频x8x8入口观看| 欧美中文日本在线观看视频| 国内精品一区二区在线观看| 欧美又色又爽又黄视频| 欧美激情在线99| 高清在线国产一区| 国产精品久久久久久亚洲av鲁大| 亚洲黑人精品在线| 看片在线看免费视频| 精品久久久久久,| 亚洲av熟女| 亚洲人成网站在线播| 韩国av一区二区三区四区| 搡老岳熟女国产| 一二三四社区在线视频社区8| 日韩 亚洲 欧美在线| 搞女人的毛片| 十八禁网站免费在线| 一进一出抽搐动态| .国产精品久久| 欧美黑人欧美精品刺激| 少妇人妻精品综合一区二区 | 亚洲av电影在线进入| 精品午夜福利视频在线观看一区| 免费在线观看影片大全网站| 97超视频在线观看视频| 亚洲av电影在线进入| .国产精品久久| 日韩人妻高清精品专区| 日韩欧美一区二区三区在线观看| 国产精品免费一区二区三区在线| 色哟哟哟哟哟哟| 亚洲18禁久久av| 久久天躁狠狠躁夜夜2o2o| 欧美另类亚洲清纯唯美| 国产探花极品一区二区| 最好的美女福利视频网| 舔av片在线| 精品熟女少妇八av免费久了| 久久国产精品人妻蜜桃| 麻豆成人av在线观看| 日本 av在线| 天美传媒精品一区二区| 桃色一区二区三区在线观看| 麻豆一二三区av精品| 久久精品91蜜桃| 免费av毛片视频| 最近中文字幕高清免费大全6 | 日本免费一区二区三区高清不卡| 国产主播在线观看一区二区| 日韩欧美免费精品| 别揉我奶头~嗯~啊~动态视频| 亚洲av不卡在线观看| 午夜日韩欧美国产| 12—13女人毛片做爰片一| 国产91精品成人一区二区三区| 亚洲七黄色美女视频| 亚洲一区高清亚洲精品| 在线国产一区二区在线| av在线蜜桃| 国产成年人精品一区二区| 免费观看的影片在线观看| 国产久久久一区二区三区| 永久网站在线| 免费搜索国产男女视频| 中文字幕av在线有码专区| 久久精品国产自在天天线| 毛片女人毛片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av.av天堂| 欧美xxxx性猛交bbbb| 日韩精品中文字幕看吧| 九九热线精品视视频播放| 欧美极品一区二区三区四区| 国产精品一区二区性色av| 日本免费a在线| 老女人水多毛片| 人人妻,人人澡人人爽秒播| 久久午夜福利片| 真人做人爱边吃奶动态| 我的女老师完整版在线观看| 久久草成人影院| 国产免费一级a男人的天堂| 亚洲av一区综合| 老女人水多毛片| 村上凉子中文字幕在线| 日韩 亚洲 欧美在线| 国产黄片美女视频| 黄色视频,在线免费观看| 国产av一区在线观看免费| av在线观看视频网站免费| 久久久久九九精品影院| 白带黄色成豆腐渣| 精品久久国产蜜桃| 欧美高清性xxxxhd video| 午夜福利高清视频| 高清毛片免费观看视频网站| 色5月婷婷丁香| 国产精品免费一区二区三区在线| 亚洲avbb在线观看| 99国产极品粉嫩在线观看| 男人的好看免费观看在线视频| 中文资源天堂在线| 国产私拍福利视频在线观看| 色视频www国产| 亚洲欧美清纯卡通| 国产精品一区二区三区四区免费观看 | av福利片在线观看| 床上黄色一级片| 国产v大片淫在线免费观看| 日韩中文字幕欧美一区二区| 国产精品爽爽va在线观看网站| 蜜桃久久精品国产亚洲av| 日韩人妻高清精品专区| 1024手机看黄色片| 国产毛片a区久久久久| 亚洲av成人av| 亚洲天堂国产精品一区在线| 久久亚洲精品不卡| 首页视频小说图片口味搜索| 欧美日韩黄片免| 一夜夜www| 日日摸夜夜添夜夜添小说| 国产高清激情床上av| 国产成人aa在线观看| 国语自产精品视频在线第100页| 欧美一区二区国产精品久久精品| 亚洲不卡免费看| 自拍偷自拍亚洲精品老妇| 国模一区二区三区四区视频| 亚洲,欧美,日韩| 免费高清视频大片| 国产 一区 欧美 日韩| 亚洲精品成人久久久久久| 亚洲人成电影免费在线| 亚洲中文字幕日韩| 精品久久久久久久末码| 久久久久国产精品人妻aⅴ院| 色哟哟·www| 国产综合懂色| 一区福利在线观看| 桃色一区二区三区在线观看| 12—13女人毛片做爰片一| 一区二区三区免费毛片| 成年版毛片免费区| 国产伦人伦偷精品视频| 国产av麻豆久久久久久久| 国产私拍福利视频在线观看| 色尼玛亚洲综合影院| 免费看a级黄色片| 欧美成人一区二区免费高清观看| 一本精品99久久精品77| 天天一区二区日本电影三级| 国产三级中文精品| АⅤ资源中文在线天堂| 欧美丝袜亚洲另类 | h日本视频在线播放| 免费高清视频大片| 老司机福利观看| 国产视频一区二区在线看| 精品乱码久久久久久99久播| 亚洲精品久久国产高清桃花| 久久久精品大字幕| 国产精品综合久久久久久久免费| 噜噜噜噜噜久久久久久91| 色视频www国产| 白带黄色成豆腐渣| 亚洲av免费在线观看| 亚洲人成电影免费在线| 91在线观看av| 亚洲国产精品sss在线观看| 99久久无色码亚洲精品果冻| av黄色大香蕉| 老熟妇仑乱视频hdxx| 久久久精品欧美日韩精品| 老司机福利观看| 亚洲 欧美 日韩 在线 免费| 久久久久国内视频| 在线观看免费视频日本深夜| 亚洲精品影视一区二区三区av| 日本熟妇午夜| 性色av乱码一区二区三区2| 国产伦精品一区二区三区四那| 国产精品一及| 久久国产精品人妻蜜桃| 亚洲av美国av| 真实男女啪啪啪动态图| 亚洲国产精品999在线| 亚洲国产精品成人综合色| 乱码一卡2卡4卡精品| 欧美+亚洲+日韩+国产| 禁无遮挡网站| 欧美日韩黄片免| 午夜免费激情av| 国产黄片美女视频| 久久午夜福利片| 亚洲色图av天堂| 两个人视频免费观看高清| 国产精品不卡视频一区二区 | 欧美日韩瑟瑟在线播放| 亚洲狠狠婷婷综合久久图片| 看黄色毛片网站| 免费高清视频大片| 可以在线观看毛片的网站| 久久久久久久久久黄片| 十八禁国产超污无遮挡网站| 精品乱码久久久久久99久播| 乱码一卡2卡4卡精品| 久久久成人免费电影| 国产成人啪精品午夜网站| 少妇被粗大猛烈的视频| 一卡2卡三卡四卡精品乱码亚洲| 美女大奶头视频| 久久久久性生活片| 日韩欧美国产一区二区入口| 免费在线观看亚洲国产| 亚洲美女黄片视频| 亚洲欧美精品综合久久99| 亚洲精品一区av在线观看| 亚洲成人久久爱视频| 99久久久亚洲精品蜜臀av| 日韩亚洲欧美综合| 精品人妻偷拍中文字幕| 可以在线观看的亚洲视频| 色哟哟哟哟哟哟| 国产精品嫩草影院av在线观看 | 丁香六月欧美| 中文字幕人成人乱码亚洲影| 97热精品久久久久久| 亚洲午夜理论影院| 97超视频在线观看视频| 一级av片app| 老熟妇仑乱视频hdxx| 午夜福利在线观看吧| 国产伦精品一区二区三区四那| 香蕉av资源在线| 欧美中文日本在线观看视频| 日本黄大片高清| 看十八女毛片水多多多| 国产一区二区激情短视频| 别揉我奶头 嗯啊视频| 99久久久亚洲精品蜜臀av| 91在线精品国自产拍蜜月| 99热精品在线国产| 热99在线观看视频| 国产精品美女特级片免费视频播放器| 黄色一级大片看看| a级毛片免费高清观看在线播放| 内地一区二区视频在线| 成人毛片a级毛片在线播放| 日韩av在线大香蕉| 中文字幕久久专区| 久久国产乱子免费精品| 好男人在线观看高清免费视频| 亚洲欧美日韩东京热| 久久精品影院6| 别揉我奶头 嗯啊视频| 亚洲真实伦在线观看| 男人的好看免费观看在线视频| 人妻制服诱惑在线中文字幕| 国产av不卡久久| 欧美色欧美亚洲另类二区| 国产乱人伦免费视频| 免费观看精品视频网站| 波野结衣二区三区在线| 日韩中字成人| 国产成人aa在线观看| 99视频精品全部免费 在线| 欧美一区二区亚洲| 亚洲人成网站高清观看| 国产成人啪精品午夜网站| 黄色日韩在线| 欧美一区二区国产精品久久精品| 亚洲精品乱码久久久v下载方式| 国产精品嫩草影院av在线观看 | 麻豆av噜噜一区二区三区| 日韩中文字幕欧美一区二区| 听说在线观看完整版免费高清| 亚洲国产精品久久男人天堂| 久久热精品热| 亚洲av免费高清在线观看| 狠狠狠狠99中文字幕| 亚洲男人的天堂狠狠| 亚洲中文字幕一区二区三区有码在线看| 嫩草影院精品99| 级片在线观看| 午夜a级毛片| 亚洲成a人片在线一区二区| 国产极品精品免费视频能看的| 亚洲欧美日韩东京热| 精品久久久久久久久av| 亚洲国产精品久久男人天堂| 又粗又爽又猛毛片免费看| 男插女下体视频免费在线播放| 国产伦一二天堂av在线观看| 草草在线视频免费看| 亚洲人成伊人成综合网2020| 国产极品精品免费视频能看的| 日韩av在线大香蕉| 好看av亚洲va欧美ⅴa在| 亚洲内射少妇av| 真实男女啪啪啪动态图| 在线观看一区二区三区| 麻豆av噜噜一区二区三区| 久久久久久九九精品二区国产| 18禁裸乳无遮挡免费网站照片| 激情在线观看视频在线高清| 欧美性猛交╳xxx乱大交人| 国产在线精品亚洲第一网站| 在线观看一区二区三区| 黄色丝袜av网址大全| 亚洲,欧美精品.| 国产欧美日韩精品一区二区| 极品教师在线免费播放| 婷婷精品国产亚洲av| 国产成人影院久久av| 国产一区二区三区视频了| 午夜激情欧美在线| 五月伊人婷婷丁香| 国产色爽女视频免费观看| 精品久久久久久久久av| 窝窝影院91人妻| 欧美性猛交╳xxx乱大交人| 中文字幕av在线有码专区| 床上黄色一级片| aaaaa片日本免费| 女同久久另类99精品国产91| 男女之事视频高清在线观看| 午夜视频国产福利| 精品午夜福利在线看| 日韩欧美在线二视频| 美女免费视频网站| 99热这里只有是精品在线观看 | 最后的刺客免费高清国语| 亚洲成av人片免费观看| 国产精品99久久久久久久久| 国产av不卡久久| 老女人水多毛片| 宅男免费午夜| 国产高清三级在线| 久久精品国产清高在天天线| 精品久久久久久成人av| 国产高清激情床上av| 亚洲美女搞黄在线观看 | 午夜久久久久精精品| 3wmmmm亚洲av在线观看| 欧美黄色淫秽网站| 亚洲自拍偷在线| 少妇人妻精品综合一区二区 | 日韩欧美在线乱码| 国产成年人精品一区二区| 亚洲美女搞黄在线观看 | 国产欧美日韩精品一区二区| 久久中文看片网| 国产成人啪精品午夜网站| 18禁黄网站禁片免费观看直播| 97超级碰碰碰精品色视频在线观看| 成人毛片a级毛片在线播放| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 九色成人免费人妻av| 午夜精品一区二区三区免费看| 天堂动漫精品| 国产熟女xx| 久久亚洲真实| 亚洲成人久久性| 国产成人欧美在线观看| 亚洲综合色惰| 欧美高清成人免费视频www| 免费搜索国产男女视频| 国产爱豆传媒在线观看| 校园春色视频在线观看| www.999成人在线观看| 欧美日本视频| 久久天躁狠狠躁夜夜2o2o| 久久国产精品影院| 亚洲人与动物交配视频| 亚洲片人在线观看| 欧美最黄视频在线播放免费| 成人午夜高清在线视频| 日韩欧美 国产精品| 99久国产av精品| 在线免费观看的www视频| 中文字幕av在线有码专区| 欧美黑人巨大hd| 亚洲专区中文字幕在线| 日韩欧美在线二视频| 成年女人永久免费观看视频| 国产白丝娇喘喷水9色精品| 久久性视频一级片| ponron亚洲| 一卡2卡三卡四卡精品乱码亚洲| 久久精品影院6| 欧美精品啪啪一区二区三区| 精品一区二区三区视频在线| 深夜a级毛片| 九色国产91popny在线| 噜噜噜噜噜久久久久久91| 变态另类成人亚洲欧美熟女| x7x7x7水蜜桃| 免费电影在线观看免费观看| 91在线观看av| 亚洲五月婷婷丁香| 亚洲国产精品sss在线观看| 亚洲在线观看片| 免费观看人在逋| 亚洲中文字幕日韩| 又爽又黄无遮挡网站| 最新中文字幕久久久久| 精品久久久久久久久久免费视频| 搡老妇女老女人老熟妇| 99久久久亚洲精品蜜臀av| 欧美日韩亚洲国产一区二区在线观看| 欧美精品国产亚洲| 99久久九九国产精品国产免费| 可以在线观看的亚洲视频| 国产精品av视频在线免费观看| 波多野结衣巨乳人妻| 一区二区三区免费毛片| 自拍偷自拍亚洲精品老妇| 男女那种视频在线观看| 在线观看美女被高潮喷水网站 | 国产精品久久视频播放| 人人妻,人人澡人人爽秒播| 中文在线观看免费www的网站| 嫩草影院精品99| 此物有八面人人有两片| 亚洲男人的天堂狠狠| a级毛片免费高清观看在线播放| 亚洲国产精品合色在线| 欧美日韩福利视频一区二区| 最新中文字幕久久久久| 日韩欧美在线乱码| 午夜福利视频1000在线观看| av国产免费在线观看| av专区在线播放| 久久午夜福利片| 亚洲美女黄片视频| 脱女人内裤的视频| 亚洲精品影视一区二区三区av| 午夜福利在线在线| 午夜福利成人在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 精品福利观看| 久久99热6这里只有精品| 精品午夜福利在线看| 激情在线观看视频在线高清| 国产精品美女特级片免费视频播放器| 丰满的人妻完整版| 好看av亚洲va欧美ⅴa在| 真人一进一出gif抽搐免费| av欧美777| 欧美日本视频| 国产欧美日韩一区二区精品| 草草在线视频免费看| 91狼人影院| 一级作爱视频免费观看| 日本 av在线| 国产熟女xx| 1024手机看黄色片| 黄色女人牲交| 内射极品少妇av片p| 日韩有码中文字幕| 国产91精品成人一区二区三区| 日本熟妇午夜| 少妇人妻精品综合一区二区 | 成熟少妇高潮喷水视频| a级毛片a级免费在线| 三级男女做爰猛烈吃奶摸视频| 午夜精品久久久久久毛片777| 亚洲成人精品中文字幕电影| 亚洲在线自拍视频| 男女视频在线观看网站免费| 日韩精品青青久久久久久| 午夜福利成人在线免费观看| 欧美黑人欧美精品刺激| 国产探花极品一区二区| 热99re8久久精品国产| 欧美黑人欧美精品刺激| 网址你懂的国产日韩在线| 欧美3d第一页| 色综合站精品国产| 永久网站在线|