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

    氣候變化對(duì)闊葉紅松林潛在地理分布區(qū)的影響

    2017-02-08 06:13:58馬芳芳周旺明于大炮代力民
    生態(tài)學(xué)報(bào) 2017年2期
    關(guān)鍵詞:紅松林分布區(qū)闊葉

    賈 翔,馬芳芳,周旺明,周 莉,于大炮, 秦 靜, 代力民,*

    1 中國(guó)科學(xué)院沈陽(yáng)應(yīng)用生態(tài)研究所,森林生態(tài)與管理重點(diǎn)實(shí)驗(yàn)室,沈陽(yáng) 110016 2 中國(guó)科學(xué)院大學(xué),北京 100049 3 遼寧省林業(yè)種苗管理總站,沈陽(yáng) 110036

    氣候變化對(duì)闊葉紅松林潛在地理分布區(qū)的影響

    賈 翔1,2,馬芳芳1,2,周旺明1,周 莉1,于大炮1, 秦 靜3, 代力民1,*

    1 中國(guó)科學(xué)院沈陽(yáng)應(yīng)用生態(tài)研究所,森林生態(tài)與管理重點(diǎn)實(shí)驗(yàn)室,沈陽(yáng) 110016 2 中國(guó)科學(xué)院大學(xué),北京 100049 3 遼寧省林業(yè)種苗管理總站,沈陽(yáng) 110036

    物種地理分布主要取決于它對(duì)氣候、地形等環(huán)境因子的適應(yīng)性?;?2個(gè)環(huán)境因子和闊葉紅松林的4類(lèi)主要建群樹(shù)種——紅松、紫椴、水曲柳和蒙古櫟的地理分布數(shù)據(jù),采用最大熵模型模擬了闊葉紅松林的潛在分布區(qū)域,并分析決定闊葉紅松林地理分布的主要?dú)夂蚝偷匦我蜃?最后利用政府間氣候變化專(zhuān)門(mén)委員會(huì)(IPCC)發(fā)布的3種排放場(chǎng)景(SRES-A2、SRES-A1B、SRES-B1)下2020、2050、2080年的氣候數(shù)據(jù)預(yù)測(cè)闊葉紅松林的未來(lái)潛在分布區(qū)。結(jié)果表明:各樹(shù)種的受試者工作特征曲線下面積(AUC值)都大于0.8,說(shuō)明模型有很好的預(yù)測(cè)能力;影響闊葉紅松林分布的主導(dǎo)環(huán)境因子是年降雨量、季節(jié)性降雨量、海拔、年平均溫度、最濕季度的平均溫度。在基準(zhǔn)氣候條件下,闊葉紅松林的高度適宜分布區(qū)主要分布在長(zhǎng)白山和小興安嶺地區(qū),占研究區(qū)總面積的11.69%,低度適宜區(qū)面積、不適宜區(qū)面積分別占研究區(qū)總面積的23%和65.31%。模型預(yù)測(cè)結(jié)果顯示,未來(lái)在A2、A1B和B1氣候情景下,闊葉紅松林高度適宜區(qū)的南界與北界都向北移動(dòng),其面積有縮減的趨勢(shì),而低度適宜區(qū)的面積有增加的趨勢(shì)。

    闊葉紅松林;氣候變化;最大熵模型;潛在地理分布;氣候變化情景

    分析物種-環(huán)境之間的關(guān)系已經(jīng)成為生態(tài)學(xué)、生物地理學(xué)中的一個(gè)焦點(diǎn)問(wèn)題[1],物種在各種各樣的生物因子和非生物因子的共同作用下,都有其暫時(shí)的地理分布范圍。物種地理分布主要取決于它對(duì)氣候、地形等環(huán)境因子的適應(yīng)性,特別是森林群落優(yōu)勢(shì)物種的地理分布與氣候有著密切的關(guān)系,其地理分布受到氣候變化的深刻影響[2]。而很多大氣環(huán)流模型預(yù)測(cè)研究表明,到21世紀(jì)末,地表平均溫度將升高1.1—6.4℃[3],因此,全球氣候變暖將會(huì)對(duì)物種地理分布產(chǎn)生極大的影響。

    關(guān)于物種地理分布的研究方法已有很多,其中比較熱門(mén)的研究方法是物種分布模型(Species Distribution Models, SDMs)。物種分布模型主要是利用物種的分布數(shù)據(jù)與環(huán)境數(shù)據(jù),依據(jù)特定的算法估計(jì)物種的生態(tài)位,并投影到景觀中,以概率的形式反映物種對(duì)生境的偏好程度[4]。在所有物種分布模型中,最大熵模型(Maxent)相比而言預(yù)測(cè)效果最好[5-6]。

    闊葉紅松林樹(shù)種資源豐富,紅松(Pinuskoraiensis)為主要成林樹(shù)種,并與其他針闊葉樹(shù)種如魚(yú)鱗云衫(Piceajezoensis)、臭冷杉(Abiesnephrolepis)、紫椴(Tiliaamurensis)、水曲柳(Fraxinusmandshurica)和蒙古櫟(Quercusmongolica)等形成針闊混交林[7],是東北東部山區(qū)典型的地帶性植被,也是第四紀(jì)冰川期后保存完好的原始裸子植物群落[8],其分布中心在東北長(zhǎng)白山、張廣才嶺、小興安嶺和完達(dá)山山地。由于歷史上的過(guò)度采伐,闊葉紅松林遭到嚴(yán)重破壞,其分布面積急劇萎縮。為了實(shí)現(xiàn)闊葉紅松林的生態(tài)恢復(fù)及可持續(xù)經(jīng)營(yíng),研究闊葉紅松林在未來(lái)氣候變化情況下的潛在分布區(qū)顯得至關(guān)重要。目前,關(guān)于闊葉紅松林地理分布的研究還局限于單一的樹(shù)種[9-10],而單一樹(shù)種的分布區(qū)域不能直接代替闊葉紅松林群落的分布區(qū)域,并且這些研究?jī)H僅探討了為數(shù)不多的環(huán)境因子,如溫度累積指標(biāo)生長(zhǎng)度·日(GDD)和水分指標(biāo)可能蒸散率(PER)對(duì)其產(chǎn)生的影響[11]。本文以闊葉紅松林的4類(lèi)主要建群樹(shù)種——紅松、紫椴、水曲柳和蒙古櫟為研究對(duì)象,通過(guò)最大熵模型分析決定闊葉紅松林地理分布的主要?dú)夂蚝偷匦我蜃?并分別模擬4類(lèi)樹(shù)種的潛在分布區(qū),根據(jù)它們的共同分布區(qū)域來(lái)決定闊葉紅松林的分布區(qū),最后利用氣候變化下的未來(lái)環(huán)境數(shù)據(jù)預(yù)測(cè)闊葉紅松林的未來(lái)潛在分布區(qū),為闊葉紅松林的生態(tài)恢復(fù)和經(jīng)營(yíng)管理提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域

    闊葉紅松林主要分布于我國(guó)東北地區(qū),其地處亞歐大陸東緣,地理坐標(biāo)38°43′—53°23′N(xiāo),118°50′—135°05′E,包括遼寧、吉林、黑龍江3省以及內(nèi)蒙古自治區(qū)東部的呼倫貝爾市、興安盟、通遼市、赤峰市,南北長(zhǎng)1600多km,東西寬1400多km[12],南臨黃海和渤海,東北與俄羅斯接壤,東南以鴨綠江為界與朝鮮民主主義人民共和國(guó)隔江相望,西部與內(nèi)蒙古自治區(qū)接壤[13]。

    1.2 數(shù)據(jù)采集

    1.2.1 樹(shù)種地理分布數(shù)據(jù)

    通過(guò)中國(guó)數(shù)字植物標(biāo)本館(http://www.cvh.org.cn/)、教學(xué)標(biāo)本資源共享平臺(tái)(http://mnh.scu.edu.cn/)、中國(guó)自然保護(hù)區(qū)標(biāo)本資源共享平臺(tái)(http://www.papc.cn/)等數(shù)據(jù)庫(kù)獲得標(biāo)本采集地信息,同時(shí)查閱《中國(guó)植物志》和東三省植物志以及相關(guān)研究文獻(xiàn)確定其現(xiàn)在分布點(diǎn)。去除模糊記錄的分布點(diǎn)信息,對(duì)具有詳細(xì)信息的分布點(diǎn),利用百度拾取坐標(biāo)系統(tǒng)(http://api.map.baidu.com/lbsapi/getpoint/index.html)確定其分布點(diǎn)的經(jīng)緯度。另外,為提高預(yù)測(cè)的可靠性,剔除引種栽培的樹(shù)種分布點(diǎn)。

    1.2.2 氣候因子

    氣候基準(zhǔn)數(shù)據(jù)來(lái)源于WorldClim數(shù)據(jù)庫(kù)(http://www.worldclim.org),該數(shù)據(jù)庫(kù)收集了1950—2000年全球各地氣象站的每月氣象數(shù)據(jù),采用插值法生成全球氣候柵格數(shù)據(jù),將2000年作為基準(zhǔn)年。本文基于ArcGIS 10.0平臺(tái),提取了東北地區(qū)的19個(gè)氣候因子(Bio1—Bio19)(表1),空間分辨率為30″(約1km2)。徐影等(2002)對(duì)國(guó)際上較流行的5種大氣環(huán)流模式(HadCM、GFDL、ECHAM、CSIRO以及CGCM)進(jìn)行比較研究,研究結(jié)果表明ECHAM4和HadCM2兩個(gè)模式對(duì)東亞和中國(guó)地區(qū)的氣候模擬效果最好[14],本文采用的未來(lái)氣候數(shù)據(jù)由大氣環(huán)流模型HadCM模擬,版本為HadCM3。研究涉及的3種排放場(chǎng)景來(lái)自《IPCC排放情景特別報(bào)告》,分別代表了碳排放高速增長(zhǎng)(SRES-A2)、中速增長(zhǎng)(SRES-A1B)和低速增長(zhǎng)(SRES-B1)。其中A2情景是描述了一個(gè)非均衡的世界:各地域間生產(chǎn)力方式的趨同異常緩慢,由此導(dǎo)致人口持續(xù)增長(zhǎng),經(jīng)濟(jì)發(fā)展主要面向區(qū)域,人均經(jīng)濟(jì)增長(zhǎng)和技術(shù)變化是不連續(xù)的,并低于其他情景的發(fā)展速度;A1情景描述的世界是:經(jīng)濟(jì)增長(zhǎng)迅速、全球人口峰值將出現(xiàn)在本世紀(jì)中葉,新的更高效的技術(shù)被迅速引進(jìn),A1情景進(jìn)一步劃分為3組情景,而A1B是能源使用較為平衡的一種類(lèi)型;B1情景是描述了一個(gè)趨同的世界:全球人口數(shù)量與A1相同,但經(jīng)濟(jì)結(jié)構(gòu)向服務(wù)和信息經(jīng)濟(jì)方向迅速調(diào)整,伴之以材料密集程度的下降,以及清潔和資源高效技術(shù)的引進(jìn)[15]。以上3種排放場(chǎng)景、3個(gè)未來(lái)時(shí)段(2020、2050、2080年)的未來(lái)氣候數(shù)據(jù)由國(guó)際熱帶農(nóng)業(yè)中心(CIAT, http://www.ccafs-climate.org)提供,空間分辨率為30″(1km2)。

    表1 研究采用的環(huán)境因子

    * 經(jīng)篩選后輸入模型的因子

    1.2.3 地形因子和底圖

    地形數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http://www.gscloud.cn)。下載分辨率為90m的數(shù)字高程地圖(DEM),從中提取海拔、坡度、坡向(表1),而后進(jìn)行重采樣生成分辨率為30″的柵格數(shù)據(jù)圖層。中國(guó)地圖和中國(guó)省級(jí)行政區(qū)劃圖來(lái)自國(guó)家基礎(chǔ)地理信息系統(tǒng)網(wǎng)站(http://nfgis.nsdi.gov.cn/)。

    1.3 數(shù)據(jù)分析

    1.3.1 樹(shù)種地理分布區(qū)

    將樹(shù)種分布點(diǎn)整理到EXCEL中,共獲取紅松分布點(diǎn)167個(gè),蒙古櫟分布點(diǎn)210個(gè),紫椴分布點(diǎn)177個(gè),水曲柳分布點(diǎn)106個(gè)(圖1)。陳新美等對(duì)Maxent模型預(yù)測(cè)物種分布所需樣本量的研究表明,樣本量大小對(duì)Maxent模型預(yù)測(cè)物種空間分布的精度影響不大[16],所以本文數(shù)據(jù)個(gè)數(shù)可以進(jìn)行模擬。

    圖1 紅松、蒙古櫟、紫椴和水曲柳分布點(diǎn)的地理位置Fig.1 Geographic locations of Pinus koraiensis, Quercus mongolica, Tilia amurensis and Fraxinus mandshurica

    1.3.2 環(huán)境因子的篩選

    由于各環(huán)境因子之間具有一定的相關(guān)性,如果直接應(yīng)用于模型,可能會(huì)產(chǎn)生過(guò)度擬合現(xiàn)象,所以對(duì)環(huán)境因子進(jìn)行相關(guān)性分析后才可以應(yīng)用于物種分布模型[17]。參照Yang等篩選環(huán)境因子的方法[1],對(duì)22個(gè)環(huán)境因子圖層進(jìn)行多重共線性分析(SPSS 19.0)來(lái)檢驗(yàn)圖層之間的相關(guān)性,若兩個(gè)環(huán)境因子之間的Pearson指數(shù)|r|≥0.8,那么只能有一個(gè)因子選入模型,最后篩選得到10個(gè)環(huán)境因子,包括7個(gè)氣候因子(Bio1、Bio 2、Bio 3、Bio 4、Bio 8、Bio 12和Bio 15)和3個(gè)地形因子(海拔、坡度、坡向)(表1)。

    1.3.3 物種分布模型

    Maxent模型是物種分布模型中表現(xiàn)較好的一種模型,它基于可以免費(fèi)使用的Maxent軟件V3.3.3.k(http://www.cs.princeton.edu/—schapire/maxent/)。Maxent軟件運(yùn)行輸入數(shù)據(jù)包括研究區(qū)域的一組環(huán)境圖層和一個(gè)物種在該區(qū)域的分布數(shù)據(jù)。參照Moreno等研究方法[18],隨機(jī)選取25%的分布點(diǎn)作為測(cè)試集(testing data),重復(fù)運(yùn)算10次進(jìn)行建模,其他參數(shù)為默認(rèn)值。采用受試者工作特征曲線(receiver operating characteristic curve, ROC)下面積(Area under curve, AUC)來(lái)評(píng)價(jià)模型的擬合程度。選擇AUC值最大的一次運(yùn)算結(jié)果進(jìn)行物種地理分布的預(yù)測(cè),AUC值越大表明擬合越好,目前以AUC值0.5—0.6為較差,0.6—0.7為一般,0.7—0.8為較準(zhǔn)確,0.8—0.9為很準(zhǔn)確,0.9—1為極準(zhǔn)確。在環(huán)境參數(shù)設(shè)置中開(kāi)啟刀切法(Jackknife)來(lái)評(píng)價(jià)各個(gè)環(huán)境變量的權(quán)重,結(jié)合每個(gè)環(huán)境因子的貢獻(xiàn)百分比來(lái)判定主導(dǎo)因子。模型的輸出格式為ASCLL格式文件,利用ASCII to Raster轉(zhuǎn)換為柵格數(shù)據(jù),按柵格數(shù)值的大小將研究區(qū)域分為3個(gè)適宜等級(jí),具體參考政府間氣候變化專(zhuān)門(mén)委員會(huì)(IPCC)關(guān)于評(píng)估可能性的劃分標(biāo)準(zhǔn)[19]:P<0.05為不適宜區(qū);0.05≤P<0.33為低度適宜區(qū);P≥0.33為高度適宜區(qū)。本文將4類(lèi)樹(shù)種共同具有的高度適宜區(qū)作為闊葉紅松林的高度適宜區(qū),用4類(lèi)樹(shù)種的不適宜區(qū)求并集作為闊葉紅松林的不適宜區(qū),其余作為低度適宜區(qū)。

    2 結(jié)果與分析

    2.1 模型模擬的準(zhǔn)確性評(píng)價(jià)

    目前,ROC曲線分析法已經(jīng)廣泛應(yīng)用于準(zhǔn)確評(píng)價(jià)物種分布模型的預(yù)測(cè)能力。本文對(duì)4類(lèi)樹(shù)種的分布進(jìn)行模擬,各樹(shù)種的AUC值如下(表2)。其中紅松表現(xiàn)最好,其次是水曲柳、紫椴、蒙古櫟。當(dāng)Test AUC>0.8時(shí),Maxent模型具有很好的預(yù)測(cè)能力,模擬結(jié)果保持穩(wěn)定,可以預(yù)測(cè)物種分布區(qū)域[20]。

    2.2 影響闊葉紅松林分布的主導(dǎo)環(huán)境因子

    通過(guò)對(duì)4類(lèi)樹(shù)種的環(huán)境因子貢獻(xiàn)百分比疊加分析(圖2),可知影響闊葉紅松林分布的主導(dǎo)環(huán)境因子是年降雨量(Bio12)、季節(jié)性降雨量(Bio15)、海拔(ELE)、年平均溫度(Bio1)和最濕季度的平均溫度(Bio8),對(duì)于紅松、蒙古櫟、水曲柳,它們的累計(jì)貢獻(xiàn)率都大于85%,對(duì)于紫椴,它們的累計(jì)貢獻(xiàn)率大于80%,其中年降雨量為最主要的因子?;贘ackknife模塊的輸出結(jié)果也能很好地反映這些因子對(duì)各樹(shù)種分布的重要性。

    表2 各樹(shù)種模型模擬的AUC值

    2.3 闊葉紅松林適宜分布區(qū)的空間變化趨勢(shì)分析

    2.3.1 闊葉紅松林潛在分布區(qū)的模擬

    運(yùn)行Maxent模型后輸出4類(lèi)樹(shù)種的潛在分布區(qū),利用Reclassy工具劃分等級(jí),然后疊加得到闊葉紅松林分布圖(圖3)。闊葉紅松林的高度適宜分布區(qū)主要分布在長(zhǎng)白山和小興安嶺地區(qū),東起吉林琿春市、黑龍江東寧縣,北界位于黑龍江伊春市、鶴崗市、蘿北縣,向南經(jīng)小興安嶺、張廣才嶺、大青山、吉林哈達(dá)嶺、遼寧本溪,到遼寧岫巖滿族自治縣、鳳城市北部、桓仁滿族自治縣。闊葉紅松林的低度適宜區(qū)主要分布在東北地區(qū)的東南面,北界位于黑龍江黑河市,向南延伸至遼寧普蘭店市、瓦房店市。闊葉紅松林的高度適宜區(qū)面積占研究區(qū)總面積的11.69%,低度適宜區(qū)面積占研究區(qū)總面積的23%,不適宜區(qū)面積占研究區(qū)總面積的65.31%(表3)。

    圖2 各樹(shù)種的環(huán)境因子貢獻(xiàn)百分比 Fig.2 Percentage contribution of environmental variables to the species distribution models

    圖3 闊葉紅松林的潛在分布圖Fig.3 Potential distribution of broadleaved Pinus koraiensis forest

    2.3.2 未來(lái)闊葉紅松林潛在分布區(qū)的預(yù)測(cè)

    由基準(zhǔn)環(huán)境數(shù)據(jù)與物種分布數(shù)據(jù)建立的模型,預(yù)測(cè)效果良好。在此基礎(chǔ)上,將物種和環(huán)境因子之間的聯(lián)系投影到未來(lái)氣候情景中用于預(yù)測(cè)4類(lèi)樹(shù)種的未來(lái)分布,得到闊葉紅松林在3種氣候情景下,2020、2050、2080年的預(yù)測(cè)分布圖(圖4)。

    在A2、A1B、B1氣候情景下,闊葉紅松林高度適宜區(qū)將主要集中分布于長(zhǎng)白山自然保護(hù)區(qū)和小興安嶺部分地區(qū)。其分布南界與北界都向北移動(dòng),在A2、A1B情景下的2080年尤為明顯。

    在A2氣候情景下,闊葉紅松林的高度適宜區(qū)面積波動(dòng)不定,由原來(lái)的11.69%驟減到2020年的4.77%,到2050年增加到5.66%,而后2080年又減少為0.60%;闊葉紅松林的低度適宜區(qū)面積從2000年的23%呈現(xiàn)增加的趨勢(shì),到2080年為37.18%(圖4,表3)。

    圖4 不同氣候情景下,闊葉紅松林的潛在分布預(yù)測(cè)圖Fig.4 Predicted potential distribution of broadleaved Pinus koraiensis forest under different climate change scenariosⅠ: A2氣候情景;Ⅱ: A1B氣候情景;Ⅲ: B1氣候情景;1: 21世紀(jì)20年代;2: 21世紀(jì)50年代;3: 21世紀(jì)80年代。如:Ⅱ—2,指A1B氣候情景下21世紀(jì)50年代

    在A1B氣候情景下,闊葉紅松林的高度適宜區(qū)面積表現(xiàn)為減少趨勢(shì),由基準(zhǔn)年的11.69%驟減到2020年的4.71%,到2050年減小到2.17%,而2080年縮少至0.46%;闊葉紅松林的低度適宜區(qū)面積從2000年的23%增加到2020年的28.01%,至此接近穩(wěn)定。

    表3 闊葉紅松林的不同適宜等級(jí)面積百分比

    而在B1氣候情景下,闊葉紅松林的高度適宜區(qū)面積開(kāi)始稍有增加而后減少,由11.69%增長(zhǎng)至2020年的13.14%,但在2050和2080年其面積分別降為6.81%和4.36%;闊葉紅松林的低度適宜區(qū)面積從2000年(23%)到2020年(23.54%)幾乎沒(méi)有變化,但是到2050年增加為30.10%,而后保持穩(wěn)定。

    3 討論

    已有研究表明,限制植物地理分布的主要因子是植物的耐寒性、完成生活史所需的熱量供應(yīng)和可利用水[21]。殷曉潔等研究表明影響蒙古櫟的最主要因子是年均降水量[10],本文得到與其相同的結(jié)論。孫曉紅研究表明影響紅松地理分布的最主要?dú)夂蛞蜃邮亲罾湓伦顭嵩缕骄鶞囟炔?其次是年降水量[9],與本研究結(jié)果有一定的差異。其主要原因是不同的研究所采用的環(huán)境因子不同,本文增加了地形因子作為環(huán)境因子,此外,由于環(huán)境因子的篩選過(guò)程不同,也會(huì)將不同的環(huán)境因子選入模型。

    對(duì)東北地區(qū)氣溫和降水變化的研究表明,在A2、A1B、B1三種排放情景下未來(lái)100年的氣溫和降水總體趨勢(shì)均呈逐漸增加的趨勢(shì)。此外,在A2情景下氣溫和降水的增幅最高,其次是A1B情景,變化幅度最小的是B1排放情景[22]。有研究表明紅松地理分布區(qū)面積與降水量呈正相關(guān)關(guān)系,與年平均溫度呈負(fù)相關(guān)關(guān)系[9],而全球氣候變化將導(dǎo)致東北地區(qū)溫度明顯升高,溫度帶北移,降水量有所增加,但抵消不了由溫度升高造成的蒸散量的增加,因此全球氣候變化將使東北東部山區(qū)面臨暖干的氣候變化趨勢(shì),所以實(shí)際上氣候的暖干化促使闊葉紅松林分布區(qū)面積有減少的趨勢(shì),生態(tài)適宜性顯著下降[11]。程肖俠等研究表明,在氣候變暖降水增加的背景下,紅松有成為大興安嶺地區(qū)森林主要樹(shù)種的趨勢(shì)[23],而本文中在3種氣候情景下,低度適宜區(qū)都有顯著的北移,到2020年,其分布區(qū)域擴(kuò)展到大興安嶺地區(qū),在A2、A1B氣候情景下到2080年,大興安嶺的部分地區(qū)也將會(huì)成為闊葉紅松林的高度適宜區(qū)。但對(duì)于小興安嶺地區(qū)的研究卻有所不同,當(dāng)未來(lái)年降水增加10%左右、年氣溫增幅大于5℃時(shí),林窗模型的模擬結(jié)果表明闊葉紅松林將被蒙古櫟、紫椴和裂葉榆為主組成的闊葉林所取代[24],在CGCM2 情景下,也能夠得到相同的結(jié)論[25],而本文選取SRES-A2、SRES-A1B、SRES-B1三種氣候情景,結(jié)果有很大的差別。原因是物種分布模型目前有很多不確定性,首先,在模型選取上,現(xiàn)在有很多種類(lèi)的物種分布模型,不同的物種分布模型預(yù)測(cè)結(jié)果有所差異,其次,環(huán)境因子的選取也有很大變化,選取不同的氣候情景,預(yù)測(cè)結(jié)果會(huì)有所不同[26]。此外,在A2氣候情景下,從2020年到2050年闊葉紅松林的高度適宜區(qū)面積增加,其原因是闊葉紅松林的適宜性等級(jí)分布是采取4類(lèi)樹(shù)種的疊加方式得到的,因此可能會(huì)造成累計(jì)誤差。

    物種分布區(qū)是物種生態(tài)與進(jìn)化歷史的復(fù)雜表達(dá)形式,在不同的時(shí)空尺度上受到多種因素、不同強(qiáng)度的控制[27-28]。通常認(rèn)為決定物種分布區(qū)共有四種因素:環(huán)境因素、生物因素、物種的擴(kuò)散能力、物種適應(yīng)新環(huán)境的進(jìn)化能力[29-31]。而本文僅僅考慮環(huán)境因素中的部分環(huán)境因子作為預(yù)測(cè)指標(biāo)進(jìn)行預(yù)測(cè),與實(shí)際還有一定的差距,還需要做進(jìn)一步的研究。

    4 結(jié)論

    目前,闊葉紅松林主要分布于長(zhǎng)白山和小興安嶺地區(qū)。但是在A2、A1B、B1氣候情景下預(yù)測(cè)闊葉紅松林的分布,結(jié)果有很大的差異,其分布界線北移、高度適宜分布區(qū)域面積減小,其中長(zhǎng)白山地區(qū)的闊葉紅松林高度適宜分布區(qū)更集中于長(zhǎng)白山自然保護(hù)區(qū)。此外,氣候變化程度越嚴(yán)重,對(duì)闊葉紅松林的分布影響也越大,尤其是在A2、A1B氣候情景下,到2080年闊葉紅松林的高度適宜區(qū)面積所占研究區(qū)域面積不足1%,說(shuō)明如果不加以保護(hù)環(huán)境而使氣候得以改善,那么闊葉紅松林在未來(lái)有可能逐漸消失于東北地區(qū)。而影響闊葉紅松林分布的主導(dǎo)環(huán)境因子主要是年降雨量(Bio12)、季節(jié)性降雨量(Bio15)、海拔(ELE)、年平均溫度(Bio1)和最濕季度的平均溫度(Bio8)。

    [1] Yang X Q, Kushwaha S P S, Saran S, Xu J C, Roy P S. Maxent modeling for predicting the potential distribution of medicinal plant,JusticiaadhatodaL. in Lesser Himalayan foothills. Ecological Engineering, 2013, 51: 83-87.

    [2] 張新時(shí), 周廣勝, 高瓊, 倪健, 唐海萍. 中國(guó)全球變化與陸地生態(tài)系統(tǒng)關(guān)系研究. 地學(xué)前緣, 1997, 4(1/2): 137-144.

    [3] IPCC(Intergovernmental Panel on Climate Change). Climate change 2007: Impacts, adaptation and vulnerability//Parry M L, Canziani O F, Palutikof J P, van der Linden P J, Hanson C E, eds. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press, 2007.

    [4] 李國(guó)慶, 劉長(zhǎng)成, 劉玉國(guó), 楊軍, 張新時(shí), 郭柯. 物種分布模型理論研究進(jìn)展. 生態(tài)學(xué)報(bào), 2013, 33(16): 4827-4835.

    [5] Elith J, Graham C H, Anderson R P, Dudík M, Ferrier S, Guisan A, Hijmans R J, Huettmann F, Leathwick J R, Lehmann A, Li J, Lohmann L G, Loiselle B A, Manion G, Moritz C, Nakamura M, Nakazawa Y, Overton J M M, Peterson A T, Phillips S J, Richardson K, Scachetti-Pereira R, Schapire R E, Soberón J, Williams S, Wisz M S, Zimmermann N E. Novel methods improve prediction of species′ distributions from occurrence data. Ecography, 2006, 29(2): 129-151.

    [6] Elith J, Graham C H. Do they? How do they? Why do they differ? On finding reasons for differing performances of species distribution models. Ecography, 2009, 32(1): 66-77.

    [7] 代力民, 谷會(huì)巖, 邵國(guó)凡, 王慶禮. 中國(guó)長(zhǎng)白山闊葉紅松林. 沈陽(yáng): 遼寧科學(xué)技術(shù)出版社, 2004.

    [8] 代力民, 邵國(guó)凡. 森林經(jīng)營(yíng)決策——理論與實(shí)踐. 沈陽(yáng): 遼寧科學(xué)技術(shù)出版社, 2005.

    [9] 孫曉紅. 中國(guó)紅松地理分布與氣候的關(guān)系研究[D]. 北京: 中國(guó)科學(xué)院研究生院, 2012.

    [10] 殷曉潔, 周廣勝, 隋興華, 何奇瑾, 李榮平. 蒙古櫟地理分布的主導(dǎo)氣候因子及其閾值. 生態(tài)學(xué)報(bào), 2013, 33(1): 103-109.

    [11] 吳正方. 東北闊葉紅松林分布區(qū)生態(tài)氣候適宜性及全球氣候變化影響評(píng)價(jià). 應(yīng)用生態(tài)學(xué)報(bào), 2003, 14(5): 771-775.

    [12] 王紹強(qiáng), 周成虎, 劉紀(jì)遠(yuǎn), 李克讓, 楊曉梅. 東北地區(qū)陸地碳循環(huán)平衡模擬分析. 地理學(xué)報(bào), 2001, 56(4): 390-400.

    [13] 王新闖. 東北地區(qū)森林生態(tài)系統(tǒng)碳儲(chǔ)量空間分布格局[D]. 沈陽(yáng): 中國(guó)科學(xué)院研究生院, 2011.

    [14] 徐 影, 丁一匯, 趙宗慈. 近30年人類(lèi)活動(dòng)對(duì)東亞地區(qū)氣候變化影響的檢測(cè)與評(píng)估. 應(yīng)用氣象學(xué)報(bào), 2002, 13(5): 513-525.

    [15] 秦大河, 陳振林, 羅勇, 丁一匯, 戴曉蘇, 任賈文, 翟盤(pán)茂, 張小曳, 趙宗慈, 張德二, 高學(xué)杰, 沈永平. 氣候變化科學(xué)的最新認(rèn)知. 氣候變化研究進(jìn)展, 2007, 3(2): 63-73.

    [16] 陳新美, 雷淵才, 張雄清, 賈宏炎. 樣本量對(duì)Maxent模型預(yù)測(cè)物種分布精度和穩(wěn)定性的影響. 林業(yè)科學(xué), 2012, 48(1): 53-59.

    [17] Sillero N. What does ecological modelling model? A proposed classification of ecological niche models based on their underlying methods. Ecological Modelling, 2011, 222(8): 1343-1346.

    [18] Moreno R, Zamora R, Molina J R, Vasquez A, Herrera M A. Predictive modeling of microhabitats for endemic birds in South Chilean temperate forests using Maximum entropy (Maxent). Ecological Informatics, 2011, 6(6): 364-370.

    [19] IPCC(Intergovernmental Panel on Climate Change). Climate change 2007: The physical science basis//Solomon S Qin D, Manning M, Chen Z, Marquis M, Averyt K B, Tignor M, Miller H L, eds. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press, 2007.

    [20] Fitzpatrick M C, Gove A D, Sanders N J, Dunn R R. Climate change, plant migration, and range collapse in a global biodiversity hotspot: theBanksia(Proteaceae) of Western Australia. Global Change Biology, 2008, 14(6): 1337-1352.

    [21] Woodward F I. Climate and Plant Distribution. Cambridge: Cambridge University Press, 1987.

    [22] 王翼, 婁德君, 曲金華, 張凱. IPCC-AR4模式資料對(duì)東北地區(qū)氣候及可利用水資源的預(yù)估研究. 自然資源學(xué)報(bào), 2009, 24(9): 1647-1656.

    [23] 程肖俠, 延曉冬. 氣候變化對(duì)中國(guó)大興安嶺森林演替動(dòng)態(tài)的影響. 生態(tài)學(xué)雜志, 2007, 26(8): 1277-1284.

    [24] 鄧慧平, 吳正方, 周道瑋. 全球氣候變化對(duì)小興安嶺闊葉紅松林影響的動(dòng)態(tài)模擬研究. 應(yīng)用生態(tài)學(xué)報(bào), 2000, 11(1): 43-46.

    [25] 周丹卉, 賀紅士, 李秀珍, 周春華, 王緒高, 陳宏偉. 小興安嶺不同年齡林分對(duì)氣候變化的潛在響應(yīng). 北京林業(yè)大學(xué)學(xué)報(bào), 2007, 29(4): 110-117.

    [26] Xu D Y, Yan H. A study of the impacts of climate change on the geographic distribution ofPinuskoraiensisin China. Environment International, 2001, 27(2/3): 201-205.

    [27] Brown J H. Macroecology. Chicago: University of Chicago Press, 1995.

    [28] Franklin J. Mapping Species Distributions: Spatial Inference and Prediction. Cambridge: Cambridge University Press, 2009.

    [29] Soberon J, Peterson A T. Interpretation of models of fundamental ecological niches and species′ distributional areas. Biodiversity Informatics, 2005, 2: 1-10.

    [30] Soberón J, Nakamura M. Niches and distributional areas: concepts, methods, and assumptions. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(S2): 19644-19650.

    [31] Soberón J M. Niche and area of distribution modeling: a population ecology perspective. Ecography, 2010, 33(1): 159-167.

    Impacts of climate change on the potential geographical distribution of broadleaved Korean pine (Pinuskoraiensis) forests

    JIA Xiang1,2, MA Fangfang1,2, ZHOU Wangming1, ZHOU Li1, YU Dapao1, QIN Jing3, DAI Limin1,*

    1KeyLaboratoryofForestEcologyandManagement,InstituteofAppliedEcology,ChineseAcademyofSciences,Shenyang110016,China2UniversityofChineseAcademyofSciences,Beijing100049,China3SeedsandSeedlingsofWoodsManagementStation,Shenyang110036,China

    Species-environment relationships have always been a central issue in ecology and biogeography. The adaptation of species to a changing natural environment, including topographical and climatic factors, determines its geographic distribution. Accumulating evidence has revealed that Northeast China has experienced the greatest temperature increase since the 1950s, which in turn affected species distribution in this region. Broadleaved Korea pine (Pinuskoraiensis) mixed forest (BKF) is the native forest type in the southern part of Northeast China. The distribution area of this species has shrunk substantially due to historical overexploitation. It is important for forest managers to be able to predict the potential geographic distribution of BKF based on species-environment relationships. Although the distributions of individual species in BKF such asP.koraiensisandQuercusmongolicahave been reported previously, few studies have focused on the potential geographic distribution of BKF. In this study, four dominant tree species -P.koraiensis,Q.mongolica,TiliaamurensisandFraxinusmandshurica, which together account for more than 80% of the growing stock in primary BKF, were selected to represent this forest type. Nineteen climatic and three topographic variables in Northeast China that are considered to be most likely influences on the geographic distribution of tree species were selected as environmental factors. To identify the major climatic and topographic factors controlling BKF distribution and simulate the potential geographic distribution of BKF under current climatic condition, the geographic distribution records of the dominant tree species, together with the environmental factors, were used in the Maxent model. The future geographic distributions of BKF were consequently predicted for the 2020s, 2050s, and 2080s, based on three kinds of climate change scenarios (SRES-A2, SRES-A1B, SRES-B1) published by the Intergovernmental Panel on Climate Change (IPCC). For the four dominant tree species, the simulations showed that the area under the curve indexes (AUC) were 0.925, 0.890, 0.859, and 0.847, respectively. All these values exceeded 0.8, which indicates that the models had a good predictive performance. The major environmental factors affecting the distribution of BKF included annual precipitation, precipitation seasonality, elevation, annual mean temperature and mean temperature of wettest quarter. For the entire region, 11.69% of the total area was identified to be of high suitability for BKF distribution, 23% was of low suitability, and 65.31% of the area was unsuitable. Under the A2, A1B, and B1 scenarios, the model predicted that both the southern and northern boundary of the high suitability area for BKF will shift northward. Overall, the high suitability area in this region was predicted to decrease, with the extent of the decrease depending on the severity of climate change. For example, under the A2 and A1B scenarios, the high suitability area in Northeast China will be less than 1% of the region by the 2080s. Overall, the results indicate that if no effective measures are taken to mitigate climate change, there is a great possibility that BKF will disappear from Northeast China.

    broadleaved Korea pine mixed forest; climate change; Maxent model; potential geographic distribution; climate change scenarios

    科技基礎(chǔ)性工作專(zhuān)項(xiàng)資助項(xiàng)目(2015FY210200-9);中國(guó)科學(xué)院特色研究所項(xiàng)目資助項(xiàng)目(Y5YZX151YD)

    2015-08-10;

    日期:2016-06-13

    10.5846/stxb201508101680

    * 通訊作者Corresponding author.E-mail: lmdai@iae.ac.cn

    賈翔,馬芳芳,周旺明,周莉,于大炮, 秦靜, 代力民.氣候變化對(duì)闊葉紅松林潛在地理分布區(qū)的影響.生態(tài)學(xué)報(bào),2017,37(2):464-473.

    Jia X, Ma F F, Zhou W M, Zhou L, Yu D P, Qin J, Dai L M.Impacts of climate change on the potential geographical distribution of broadleaved Korean pine (Pinuskoraiensis) forests.Acta Ecologica Sinica,2017,37(2):464-473.

    猜你喜歡
    紅松林分布區(qū)闊葉
    ◆ 園林苗木
    沒(méi)有紅松的紅松林
    恩平市植物組成及區(qū)系特征研究
    廣東地區(qū)夾竹桃科植物區(qū)系分析
    紅松林采伐更新的確定方式
    物種分布區(qū)特征及其對(duì)生物多樣性保育的意義
    炎黃地理(2019年5期)2019-10-20 01:59:19
    東折棱河楓樺紅松林與柞樹(shù)紅松林結(jié)構(gòu)特征對(duì)比分析
    斜紋夜蛾(Spodopteralitura)在我國(guó)的時(shí)空分布概述
    園林苗木
    高見(jiàn)的詩(shī)
    亚洲av福利一区| 在线精品无人区一区二区三 | 国模一区二区三区四区视频| 我要看日韩黄色一级片| 综合色丁香网| 午夜免费鲁丝| 哪个播放器可以免费观看大片| 精品人妻偷拍中文字幕| 国产乱人视频| 久久99蜜桃精品久久| 日韩精品有码人妻一区| 我的女老师完整版在线观看| 日韩av免费高清视频| 在线观看美女被高潮喷水网站| 色播亚洲综合网| 成人无遮挡网站| 王馨瑶露胸无遮挡在线观看| 午夜老司机福利剧场| 亚洲精品第二区| 成年免费大片在线观看| 国产精品蜜桃在线观看| 久久99热这里只频精品6学生| 亚洲精品日韩在线中文字幕| 国产av码专区亚洲av| 日韩成人伦理影院| 另类亚洲欧美激情| 内射极品少妇av片p| 白带黄色成豆腐渣| 亚洲色图av天堂| 两个人的视频大全免费| 成人黄色视频免费在线看| 亚洲欧洲日产国产| 王馨瑶露胸无遮挡在线观看| 亚洲欧美一区二区三区国产| 男人狂女人下面高潮的视频| 新久久久久国产一级毛片| 精品久久久久久久人妻蜜臀av| 亚洲国产精品999| 你懂的网址亚洲精品在线观看| 婷婷色麻豆天堂久久| 99久久精品热视频| 国产亚洲午夜精品一区二区久久 | 亚洲精品乱码久久久v下载方式| 男女无遮挡免费网站观看| 特级一级黄色大片| 全区人妻精品视频| 久久99蜜桃精品久久| 啦啦啦在线观看免费高清www| 久久人人爽人人爽人人片va| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 边亲边吃奶的免费视频| 大码成人一级视频| 久久精品夜色国产| 日韩国内少妇激情av| 国产精品福利在线免费观看| 视频中文字幕在线观看| 熟妇人妻不卡中文字幕| 国产av码专区亚洲av| 一级毛片 在线播放| 日韩欧美精品免费久久| 丝袜喷水一区| 午夜福利网站1000一区二区三区| 在线观看国产h片| 午夜激情福利司机影院| 精品少妇黑人巨大在线播放| 日韩不卡一区二区三区视频在线| 午夜福利视频1000在线观看| 精品一区二区三卡| 男女啪啪激烈高潮av片| 国产亚洲最大av| 秋霞在线观看毛片| 久久鲁丝午夜福利片| av免费观看日本| 久久99热这里只频精品6学生| 18禁动态无遮挡网站| 人妻制服诱惑在线中文字幕| 亚洲av男天堂| 99久久精品国产国产毛片| 成人美女网站在线观看视频| 日韩人妻高清精品专区| 免费大片18禁| 精品99又大又爽又粗少妇毛片| 看免费成人av毛片| 高清日韩中文字幕在线| 亚洲国产精品成人久久小说| 国产毛片在线视频| 久久国产乱子免费精品| 欧美成人一区二区免费高清观看| 免费少妇av软件| 成人毛片a级毛片在线播放| 国产成人精品婷婷| 高清欧美精品videossex| 麻豆精品久久久久久蜜桃| 涩涩av久久男人的天堂| 午夜精品一区二区三区免费看| 久久6这里有精品| av女优亚洲男人天堂| 欧美xxxx性猛交bbbb| 青春草视频在线免费观看| 久久精品国产鲁丝片午夜精品| 秋霞伦理黄片| 成人亚洲欧美一区二区av| 精品人妻视频免费看| 在线观看国产h片| 亚洲国产色片| 成人美女网站在线观看视频| 小蜜桃在线观看免费完整版高清| 国产探花极品一区二区| 国产熟女欧美一区二区| av播播在线观看一区| 七月丁香在线播放| 深爱激情五月婷婷| 好男人在线观看高清免费视频| 如何舔出高潮| 精品久久久久久久人妻蜜臀av| 三级男女做爰猛烈吃奶摸视频| 免费看a级黄色片| 精品久久久久久电影网| 青春草亚洲视频在线观看| 午夜福利视频1000在线观看| 久久精品熟女亚洲av麻豆精品| 大片免费播放器 马上看| 国产精品蜜桃在线观看| 精品国产一区二区三区久久久樱花 | 中文字幕免费在线视频6| 最近中文字幕2019免费版| 国产成年人精品一区二区| 人妻少妇偷人精品九色| 国产欧美亚洲国产| 久久人人爽人人爽人人片va| 如何舔出高潮| 边亲边吃奶的免费视频| 我要看日韩黄色一级片| av在线app专区| 亚洲欧美精品自产自拍| 色视频在线一区二区三区| 蜜桃久久精品国产亚洲av| 九九在线视频观看精品| 国产视频首页在线观看| 亚洲精品一区蜜桃| 日韩av在线免费看完整版不卡| 99九九线精品视频在线观看视频| 欧美高清性xxxxhd video| 国产av码专区亚洲av| 九色成人免费人妻av| 99热全是精品| 熟女电影av网| 18禁在线播放成人免费| 国产大屁股一区二区在线视频| 中文天堂在线官网| 亚洲av一区综合| 内地一区二区视频在线| 国内少妇人妻偷人精品xxx网站| 国产精品精品国产色婷婷| 亚洲av不卡在线观看| 高清午夜精品一区二区三区| 在线观看美女被高潮喷水网站| 99久久人妻综合| 成人毛片a级毛片在线播放| 18禁在线播放成人免费| 伊人久久国产一区二区| 亚洲真实伦在线观看| 国产免费视频播放在线视频| 久久久a久久爽久久v久久| 天天躁夜夜躁狠狠久久av| 嫩草影院新地址| 新久久久久国产一级毛片| 午夜福利高清视频| 午夜福利视频精品| 国产在线男女| 男女下面进入的视频免费午夜| 美女cb高潮喷水在线观看| 80岁老熟妇乱子伦牲交| 精品一区二区免费观看| 三级国产精品片| 国产白丝娇喘喷水9色精品| 日韩三级伦理在线观看| 欧美日韩亚洲高清精品| 精品熟女少妇av免费看| 亚洲色图综合在线观看| 久热这里只有精品99| 美女主播在线视频| 一本一本综合久久| 免费大片黄手机在线观看| 亚洲一级一片aⅴ在线观看| 狂野欧美白嫩少妇大欣赏| 欧美激情久久久久久爽电影| 久久精品综合一区二区三区| 欧美 日韩 精品 国产| 免费观看av网站的网址| 欧美zozozo另类| 男女无遮挡免费网站观看| 97人妻精品一区二区三区麻豆| 亚洲va在线va天堂va国产| 小蜜桃在线观看免费完整版高清| 免费电影在线观看免费观看| 国产精品蜜桃在线观看| 亚洲av日韩在线播放| 国产综合精华液| 亚洲天堂国产精品一区在线| 自拍偷自拍亚洲精品老妇| 97热精品久久久久久| 久久久欧美国产精品| 日韩制服骚丝袜av| 超碰97精品在线观看| 人体艺术视频欧美日本| 日韩成人av中文字幕在线观看| 国产精品麻豆人妻色哟哟久久| 免费高清在线观看视频在线观看| 久久久色成人| 波野结衣二区三区在线| 综合色丁香网| 人人妻人人看人人澡| 国产精品av视频在线免费观看| 亚洲成人久久爱视频| 国产白丝娇喘喷水9色精品| 97在线视频观看| 内地一区二区视频在线| 精品人妻视频免费看| 国产精品嫩草影院av在线观看| 国产成人一区二区在线| 夫妻午夜视频| 精品久久久噜噜| 黄片wwwwww| 亚洲国产精品成人久久小说| 国产av码专区亚洲av| 一级二级三级毛片免费看| 亚洲天堂国产精品一区在线| 最新中文字幕久久久久| 美女高潮的动态| 欧美xxⅹ黑人| 性色avwww在线观看| 日韩在线高清观看一区二区三区| 丰满人妻一区二区三区视频av| 国产69精品久久久久777片| 在线a可以看的网站| 18禁在线播放成人免费| 日韩欧美一区视频在线观看 | 日日啪夜夜撸| 色婷婷久久久亚洲欧美| 国产精品.久久久| 丰满乱子伦码专区| 免费少妇av软件| 久久久久精品久久久久真实原创| 国产69精品久久久久777片| 国内精品美女久久久久久| 一区二区三区四区激情视频| 91狼人影院| av又黄又爽大尺度在线免费看| 欧美成人精品欧美一级黄| 欧美精品一区二区大全| 亚洲精品国产av成人精品| 特级一级黄色大片| 听说在线观看完整版免费高清| 一级爰片在线观看| 精品国产三级普通话版| 香蕉精品网在线| 男女边摸边吃奶| 美女cb高潮喷水在线观看| 看免费成人av毛片| 久久久久久久精品精品| 国产伦理片在线播放av一区| 97在线视频观看| 国产综合懂色| 国产精品秋霞免费鲁丝片| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 最后的刺客免费高清国语| 国产伦在线观看视频一区| 国产国拍精品亚洲av在线观看| 亚洲av不卡在线观看| freevideosex欧美| 少妇人妻精品综合一区二区| 亚洲精品aⅴ在线观看| 18禁在线无遮挡免费观看视频| 麻豆成人av视频| 一级a做视频免费观看| 国产乱人偷精品视频| 少妇人妻一区二区三区视频| 三级国产精品片| 午夜福利网站1000一区二区三区| 久久精品久久精品一区二区三区| 高清视频免费观看一区二区| 听说在线观看完整版免费高清| 国产探花在线观看一区二区| 一本色道久久久久久精品综合| 国产v大片淫在线免费观看| 亚洲精品视频女| 毛片女人毛片| 熟女电影av网| 久久久久九九精品影院| 色网站视频免费| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 九色成人免费人妻av| 国产 一区 欧美 日韩| 亚洲欧美日韩无卡精品| 人妻系列 视频| 国产精品一区二区性色av| 最后的刺客免费高清国语| 大香蕉久久网| 18禁在线播放成人免费| 国产 一区 欧美 日韩| 国产黄频视频在线观看| 少妇熟女欧美另类| 久久久久久九九精品二区国产| 免费观看的影片在线观看| 网址你懂的国产日韩在线| 深爱激情五月婷婷| 亚洲无线观看免费| 国产高清国产精品国产三级 | 日韩av在线免费看完整版不卡| 亚洲最大成人手机在线| 欧美三级亚洲精品| 欧美日韩国产mv在线观看视频 | 国产精品av视频在线免费观看| 男人狂女人下面高潮的视频| 亚洲精品视频女| 在线免费十八禁| 夜夜看夜夜爽夜夜摸| 六月丁香七月| .国产精品久久| 91在线精品国自产拍蜜月| 免费av观看视频| 欧美精品国产亚洲| 天堂网av新在线| 亚洲av.av天堂| 中文字幕制服av| 亚洲熟女精品中文字幕| 久久人人爽人人片av| 22中文网久久字幕| 简卡轻食公司| 亚洲精品自拍成人| 在线免费观看不下载黄p国产| 我的老师免费观看完整版| 18禁裸乳无遮挡免费网站照片| 看免费成人av毛片| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 69人妻影院| 中国美白少妇内射xxxbb| 精品一区在线观看国产| 91精品一卡2卡3卡4卡| 久久久久久国产a免费观看| 亚洲三级黄色毛片| 日韩av免费高清视频| 亚洲美女搞黄在线观看| 男的添女的下面高潮视频| 能在线免费看毛片的网站| 久久精品国产亚洲av涩爱| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产免费福利视频在线观看| 美女cb高潮喷水在线观看| 最近中文字幕2019免费版| 欧美xxxx黑人xx丫x性爽| 亚洲成色77777| 欧美激情在线99| 中国美白少妇内射xxxbb| 免费电影在线观看免费观看| 国产亚洲一区二区精品| 永久免费av网站大全| 高清欧美精品videossex| 国产av码专区亚洲av| 国语对白做爰xxxⅹ性视频网站| 免费av不卡在线播放| 欧美成人a在线观看| 成人亚洲精品一区在线观看 | 亚洲精品日本国产第一区| 国产av码专区亚洲av| 亚洲综合精品二区| 国内精品美女久久久久久| 色婷婷久久久亚洲欧美| 久久久久久九九精品二区国产| 国产 精品1| 午夜福利高清视频| 久久人人爽人人片av| 蜜臀久久99精品久久宅男| 99热6这里只有精品| 久久久午夜欧美精品| 欧美xxxx性猛交bbbb| 久久久亚洲精品成人影院| 国产精品99久久99久久久不卡 | 熟妇人妻不卡中文字幕| 亚洲综合色惰| 成人美女网站在线观看视频| 欧美国产精品一级二级三级 | 亚洲人成网站高清观看| 2018国产大陆天天弄谢| 日韩一区二区视频免费看| 久久久久久久午夜电影| 国产成人午夜福利电影在线观看| 我要看日韩黄色一级片| 久久精品国产亚洲av涩爱| 精品少妇久久久久久888优播| 内射极品少妇av片p| 久久精品久久精品一区二区三区| .国产精品久久| 欧美日韩国产mv在线观看视频 | 精品午夜福利在线看| 男人添女人高潮全过程视频| 日韩强制内射视频| 特级一级黄色大片| 少妇猛男粗大的猛烈进出视频 | 日韩 亚洲 欧美在线| 久久精品人妻少妇| 免费看av在线观看网站| 又大又黄又爽视频免费| av在线老鸭窝| 国产av不卡久久| 国产欧美日韩精品一区二区| 我要看日韩黄色一级片| 欧美激情久久久久久爽电影| 18禁裸乳无遮挡免费网站照片| 欧美三级亚洲精品| 久久久久国产精品人妻一区二区| 国产精品偷伦视频观看了| 一级av片app| 99热国产这里只有精品6| 九九爱精品视频在线观看| 亚洲美女视频黄频| 国产色爽女视频免费观看| 啦啦啦在线观看免费高清www| 美女被艹到高潮喷水动态| 国产视频首页在线观看| 不卡视频在线观看欧美| 性插视频无遮挡在线免费观看| 在线a可以看的网站| 亚洲四区av| 国产女主播在线喷水免费视频网站| 欧美xxxx性猛交bbbb| 国产精品偷伦视频观看了| 如何舔出高潮| 久久久久久伊人网av| 真实男女啪啪啪动态图| 国产亚洲午夜精品一区二区久久 | 男女国产视频网站| 2021天堂中文幕一二区在线观| 日韩 亚洲 欧美在线| 日产精品乱码卡一卡2卡三| 男人狂女人下面高潮的视频| 精品一区二区免费观看| 青春草视频在线免费观看| 激情五月婷婷亚洲| 欧美激情国产日韩精品一区| 大片免费播放器 马上看| 肉色欧美久久久久久久蜜桃 | 精品久久久久久久久亚洲| 国产爱豆传媒在线观看| 秋霞在线观看毛片| 免费播放大片免费观看视频在线观看| tube8黄色片| 久久久久久国产a免费观看| 日日摸夜夜添夜夜添av毛片| 国产av国产精品国产| 夫妻性生交免费视频一级片| 国产片特级美女逼逼视频| 亚洲精品亚洲一区二区| 特大巨黑吊av在线直播| 免费av毛片视频| 欧美97在线视频| 久久久精品欧美日韩精品| 亚洲最大成人中文| 日韩欧美精品免费久久| 男女啪啪激烈高潮av片| 欧美极品一区二区三区四区| 大码成人一级视频| 99久久精品一区二区三区| 成人漫画全彩无遮挡| 国产精品一区二区三区四区免费观看| 亚洲最大成人手机在线| 免费高清在线观看视频在线观看| 亚洲不卡免费看| 亚洲国产欧美在线一区| 国产毛片在线视频| 丰满少妇做爰视频| 天天躁夜夜躁狠狠久久av| 天堂俺去俺来也www色官网| 精品久久久久久久久亚洲| 涩涩av久久男人的天堂| a级毛片免费高清观看在线播放| 亚洲欧美日韩东京热| 欧美日韩国产mv在线观看视频 | 欧美另类一区| 人妻制服诱惑在线中文字幕| 国产亚洲精品久久久com| 国产精品国产三级专区第一集| 亚洲人成网站在线播| 国产男人的电影天堂91| 99re6热这里在线精品视频| 美女内射精品一级片tv| 91久久精品电影网| 香蕉精品网在线| 各种免费的搞黄视频| 国产在线男女| 亚洲国产精品成人综合色| 亚洲电影在线观看av| 亚洲精品成人av观看孕妇| 婷婷色综合大香蕉| 国产色婷婷99| 综合色丁香网| 亚洲第一区二区三区不卡| 青春草国产在线视频| 国内揄拍国产精品人妻在线| 亚洲三级黄色毛片| 啦啦啦中文免费视频观看日本| 国产av码专区亚洲av| 少妇人妻精品综合一区二区| 亚洲欧美成人综合另类久久久| 青春草视频在线免费观看| 在线观看美女被高潮喷水网站| 国产高清不卡午夜福利| 免费黄色在线免费观看| 99久久精品热视频| 赤兔流量卡办理| av在线播放精品| 三级国产精品欧美在线观看| 男人狂女人下面高潮的视频| 亚洲精品乱久久久久久| 日韩在线高清观看一区二区三区| 好男人视频免费观看在线| 精品人妻偷拍中文字幕| 欧美成人午夜免费资源| 久久久久久国产a免费观看| 黄色一级大片看看| 亚洲国产精品成人久久小说| av国产免费在线观看| 一级爰片在线观看| 午夜免费鲁丝| 我的老师免费观看完整版| 如何舔出高潮| 久久久国产一区二区| 日韩国内少妇激情av| 日韩欧美精品免费久久| 蜜臀久久99精品久久宅男| 亚洲精品一区蜜桃| 久久久精品欧美日韩精品| 久久久久久久久久人人人人人人| 一级黄片播放器| 欧美日韩精品成人综合77777| 丝瓜视频免费看黄片| 亚洲一级一片aⅴ在线观看| 又爽又黄无遮挡网站| 日本一本二区三区精品| 久久99热这里只有精品18| 99久久人妻综合| 七月丁香在线播放| 新久久久久国产一级毛片| 国产精品久久久久久精品电影| 久久久久久久久久久丰满| 国产伦精品一区二区三区视频9| 国产综合懂色| 亚洲自偷自拍三级| 免费人成在线观看视频色| 久久久久久久午夜电影| 色5月婷婷丁香| 偷拍熟女少妇极品色| 伦精品一区二区三区| 欧美zozozo另类| 美女脱内裤让男人舔精品视频| 国产亚洲av嫩草精品影院| 三级国产精品欧美在线观看| 一级黄片播放器| 亚洲欧美成人综合另类久久久| 国产精品麻豆人妻色哟哟久久| 日韩av免费高清视频| 欧美日本视频| 啦啦啦中文免费视频观看日本| 国产成人a区在线观看| 成人亚洲欧美一区二区av| 国产在视频线精品| 午夜福利在线观看免费完整高清在| 免费观看的影片在线观看| 亚洲精品亚洲一区二区| 啦啦啦啦在线视频资源| av免费观看日本| 内地一区二区视频在线| 99热网站在线观看| 欧美一级a爱片免费观看看| 日韩av免费高清视频| 18+在线观看网站| 午夜亚洲福利在线播放| 国产高清不卡午夜福利| 国产片特级美女逼逼视频| av在线观看视频网站免费| 久久99热6这里只有精品| 亚洲欧洲国产日韩| 真实男女啪啪啪动态图| 激情五月婷婷亚洲| 日韩国内少妇激情av| 精品久久久精品久久久| 成人毛片a级毛片在线播放| 国产精品熟女久久久久浪| 男女边摸边吃奶| videossex国产| 免费看不卡的av| 免费观看无遮挡的男女| 99九九线精品视频在线观看视频| h日本视频在线播放| 日韩亚洲欧美综合| freevideosex欧美| 国产精品蜜桃在线观看| 成人国产av品久久久| 亚洲欧美精品自产自拍| 51国产日韩欧美| 国产精品国产三级专区第一集| 丝瓜视频免费看黄片| 一级片'在线观看视频| 欧美+日韩+精品| 2021少妇久久久久久久久久久| 国产午夜精品久久久久久一区二区三区| 亚洲美女搞黄在线观看| 一级片'在线观看视频| 在线精品无人区一区二区三 | 亚洲综合精品二区|