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

    基于MaxEnt模型的當(dāng)歸全球生態(tài)適宜區(qū)和生態(tài)特征研究

    2017-09-12 13:00:38張東方孫成忠謝彩香
    生態(tài)學(xué)報 2017年15期
    關(guān)鍵詞:概率因子土壤

    張東方, 張 琴, 郭 杰,孫成忠,吳 杰,聶 祥,謝彩香,*

    1 中國醫(yī)學(xué)科學(xué)院, 北京協(xié)和醫(yī)學(xué)院 藥用植物研究所,北京 100193 2 江蘇師范大學(xué),徐州 221116 3 中國測繪科學(xué)研究院,北京 100039 4 中國中醫(yī)科學(xué)院中藥研究所,北京 100029 5 貴州省畢節(jié)市氣象局,畢節(jié) 551714

    基于MaxEnt模型的當(dāng)歸全球生態(tài)適宜區(qū)和生態(tài)特征研究

    張東方1,2, 張 琴1, 郭 杰1,孫成忠3,吳 杰4,聶 祥5,謝彩香1,*

    1 中國醫(yī)學(xué)科學(xué)院, 北京協(xié)和醫(yī)學(xué)院 藥用植物研究所,北京 100193 2 江蘇師范大學(xué),徐州 221116 3 中國測繪科學(xué)研究院,北京 100039 4 中國中醫(yī)科學(xué)院中藥研究所,北京 100029 5 貴州省畢節(jié)市氣象局,畢節(jié) 551714

    研究藥用植物當(dāng)歸在全球范圍內(nèi)的生態(tài)適宜區(qū)和生態(tài)特征,根據(jù)109個當(dāng)歸樣本點和37個生態(tài)因子,利用最大信息熵模型(MaxEnt)和地理信息系統(tǒng)(Geographic information System, GIS)對其進行全球適宜性區(qū)劃。結(jié)果表明,當(dāng)歸生態(tài)適宜區(qū)主要集中在北緯20—50°范圍內(nèi)的北美洲、歐洲、亞洲和南緯15—35°范圍內(nèi)的南美洲和非洲,共約593.07萬km2。其中生態(tài)相似度最高的區(qū)域主要分布在中國甘肅省南部、四川省、西藏自治區(qū)東部、云南省北部、貴州省、陜西省西南部等地區(qū),可見中國十分適合當(dāng)歸的種植,是當(dāng)歸的主產(chǎn)區(qū)。影響當(dāng)歸地理分布的主要氣候因子適宜值范圍分別為:最暖季降水量為300—700 mm,最冷季平均溫為-3—7℃,降水的季節(jié)性數(shù)值為70—98;影響當(dāng)歸地理分布的主要土壤因子適宜值范圍分別為:碳酸鈣含量為50%—60%,粘土比重為17%—24%。研究結(jié)果為當(dāng)歸人工引種栽培及精細(xì)化種植管理提供參考。

    當(dāng)歸;最大信息熵模型;地理信息系統(tǒng);生態(tài)適宜性;生態(tài)特征

    當(dāng)歸為傘形科(Umbelliferae)當(dāng)歸屬(Angelica)植物當(dāng)歸(Angelicasinensis(Oliv.)Diels.)的干燥根,是臨床常用中藥,素有“十方九歸”之稱[1],現(xiàn)代藥理研究表明:當(dāng)歸在心臟保護、胃腸保護、抗輻射和抗腫瘤以及鎮(zhèn)痛等方面均有較好療效[2- 5]。當(dāng)歸作為傳統(tǒng)出口大宗藥材,不僅用于配方,而且在傳統(tǒng)膳食及保健品中均有較多應(yīng)用,全年銷量達到13000—15000 t[6]。

    當(dāng)歸屬植物主要分布于北溫帶地區(qū),北美和東亞為其世界分布中心,東亞則以中國種類最豐,在日本、朝鮮和越南也有當(dāng)歸栽培品[7]。我國當(dāng)歸產(chǎn)地廣泛分布于甘肅、云南、貴州、四川、陜西、青海、湖南、湖北等地[8- 9]。目前當(dāng)歸主要來源于栽培,雖然有著悠久的栽培歷史及廣泛的應(yīng)用前景,但是當(dāng)歸存在種質(zhì)資源退化、栽培技術(shù)落后、質(zhì)量下降等問題,而且栽培過程中還存在著嚴(yán)重的連作障礙,導(dǎo)致當(dāng)歸可栽培土地逐漸減少,其產(chǎn)量遠(yuǎn)遠(yuǎn)跟不上市場需求,嚴(yán)重影響著當(dāng)歸的現(xiàn)代化發(fā)展。

    氣候和土壤是決定物種分布的最主要因素,近年來已經(jīng)有很多國內(nèi)外學(xué)者研究了環(huán)境因子對植物潛在分布的影響[10- 11],但對當(dāng)歸的研究目前主要集中在植物學(xué)、栽培學(xué)、藥理藥效、臨床應(yīng)用和化學(xué)成分提取等方面,當(dāng)歸在全球的資源分布及生態(tài)適宜性研究鮮有報道。為了避免盲目引種擴種和連作障礙等問題帶來的經(jīng)濟損失,研究當(dāng)歸全球適宜性區(qū)劃及其生態(tài)特征,對當(dāng)歸科學(xué)引種栽培及精細(xì)化種植管理具有重大意義。本文利用GIS技術(shù)和最大信息熵模型,通過數(shù)據(jù)庫和文獻檢索,收集當(dāng)歸分布信息,結(jié)合氣候、地形和土壤等相關(guān)生態(tài)因子,對當(dāng)歸地理分布進行區(qū)劃,找出影響當(dāng)歸分布的主要生態(tài)因子以及最適合生長的區(qū)域,為當(dāng)歸人工引種栽培及選址提供參考,以提高當(dāng)歸質(zhì)量和產(chǎn)量,實現(xiàn)中藥產(chǎn)業(yè)可持續(xù)發(fā)展。

    1 材料

    1.1 最大信息熵模型

    最大信息熵模型MaxEnt軟件(版本為3.3.3 k)來源于官網(wǎng)(http://www.cs.princeton.edu/—schapire/maxtent),MaxEnt模型是一種基于最大熵理論而提出的生態(tài)位模型[12-13],是把研究區(qū)所有柵格單元作為最大熵的可能分布空間,將已知物種分布點的柵格單元作為樣點,根據(jù)樣點單元的環(huán)境變量得出約束條件,尋找此約束條件下的最大熵的可能分布[14-15],具有操作簡單、運算速度快、預(yù)測結(jié)果好等優(yōu)點[16-19]。

    1.2 當(dāng)歸分布數(shù)據(jù)

    根據(jù)全球生物多樣性信息數(shù)據(jù)庫網(wǎng)絡(luò)(GBIF,http://www.gbif.org/)和中國數(shù)字植物標(biāo)本館(CVH,http://www.cvh.org.cn/)以及文獻檢索得到全球109個當(dāng)歸樣本點,主要分布在亞洲、美洲等(如圖1)。

    圖1 當(dāng)歸樣本分布圖Fig.1 Sample distribution of Angelica sinensis

    1.2 環(huán)境數(shù)據(jù)

    氣候數(shù)據(jù)和高程數(shù)據(jù)來源于世界氣候網(wǎng)(http://www.worldclim.org/),坐標(biāo)系為WGS84,柵格大小約為25km2。氣候數(shù)據(jù)為1950年至2000年監(jiān)測數(shù)據(jù)的平均值,共19個氣候因子,圖層中的溫度數(shù)值(℃)為實際數(shù)值×10。地形數(shù)據(jù)包括高程、坡度和坡向數(shù)據(jù),其中坡度和坡向數(shù)據(jù)利用ArcGIS軟件對高程數(shù)據(jù)進行表面分析得到。土壤數(shù)據(jù)來自世界土壤數(shù)據(jù)庫(http://webarchive.iiasa.ac.at/Research / LUC/External-World-soil-database/ HTML/),坐標(biāo)系為WGS84,柵格大小約為1km2,共15個土壤因子,則全部生態(tài)因子共37個。地圖資料來源于國家基礎(chǔ)地理信息系統(tǒng)(http://nfgis.nsdi.gov. cn/)。

    2 方法

    2.1 數(shù)據(jù)預(yù)處理

    采用ENVI4.8和ArcGIS 10.0對土壤數(shù)據(jù)進行重采樣,得到與氣候數(shù)據(jù)相同的柵格大小,同時對土壤數(shù)據(jù)進行格式轉(zhuǎn)換,得到與氣候數(shù)據(jù)和海拔數(shù)據(jù)相同的數(shù)據(jù)格式(.asc),可直接加載于MaxEnt軟件。當(dāng)歸樣本點數(shù)據(jù)文件保存為csv格式。

    2.2 MaxEnt模型參數(shù)設(shè)置

    將分布點數(shù)據(jù)(csv格式)和生態(tài)因子數(shù)據(jù)集導(dǎo)入MaxEnt軟件中,設(shè)置參數(shù)建模運行,具體參數(shù)設(shè)置為:設(shè)置刀切法檢驗權(quán)重,MaxEnt自定義設(shè)置ROC評價曲線,其他參數(shù)設(shè)置為軟件默認(rèn)值。

    2.3 區(qū)劃結(jié)果可信度和準(zhǔn)確度

    本實驗應(yīng)用ROC曲線和AUC值對預(yù)測結(jié)果進行精度評價。ROC曲線分析法在物種潛在分布預(yù)測模型評價中得到了廣泛的應(yīng)用,是目前認(rèn)可度較高的診斷試驗評價指標(biāo)[20-22]。通過計算ROC接受曲線下方的面積得到AUC值,AUC的數(shù)值范圍為0.5—1,值越大表示預(yù)測越精確。

    2.4 當(dāng)歸生態(tài)適宜區(qū)劃

    利用ArcGIS 10.0軟件對MaxEnt模型計算的結(jié)果進行疊加分析和地圖制作,繪制出基于主要生態(tài)因子的當(dāng)歸生態(tài)適宜區(qū)劃全球圖,采用人工分級方法基于生態(tài)相似度劃分出不同等級分布區(qū)域,其生態(tài)相似度共劃分為<10%,10%—30%,30%—50%,50%—70%和>70% 5個等級,并確定當(dāng)歸生態(tài)適宜分布區(qū)與行政歸屬范圍。通過MaxEnt模型運行結(jié)果中的響應(yīng)曲線得出生態(tài)因子適宜值范圍。

    圖2 基于MaxEnt模型預(yù)測當(dāng)歸全球分布的ROC曲線 Fig.2 The ROC curve predicting the distribution of Angelica sinensis based on MaxEnt

    3 結(jié)果與分析

    3.1 區(qū)劃結(jié)果可信度和準(zhǔn)確度

    AUC指的是ROC曲線與橫坐標(biāo)圍成的區(qū)域面積,其大小反映了模型模擬效果的準(zhǔn)確性,數(shù)值越大,模型模擬的準(zhǔn)確性越高。如圖2所示,訓(xùn)練集AUC值達到0.989,測試集AUC值達到了0.967,表明模型模擬效果達到極高的水平,由模型運算得出的當(dāng)歸生態(tài)適宜區(qū)劃具有較高的可信度和準(zhǔn)確度。

    3.2 生態(tài)因子貢獻率

    MaxEnt模型自動輸出各個生態(tài)因子貢獻率,其中貢獻率較大的幾個生態(tài)因子從大到小依次是最暖季降水量、最冷季平均溫、碳酸鈣含量、降水的季節(jié)性、坡度、氣溫的季節(jié)性、粘土比重、海拔、土壤中粘粒組的陽離子交換能力、氣溫年范圍等,其累計貢獻率為92%,見表1。

    3.3 當(dāng)歸分布區(qū)氣候特征

    基于各氣候因子響應(yīng)曲線(圖3)可知,當(dāng)存在概率大于0.5時,其對應(yīng)的生態(tài)因子值比較適合當(dāng)歸的生長[23]。

    (1)最暖季降水量(貢獻率26.5%)在約250mm以下時,當(dāng)歸的存在概率極低,在0.5以下。隨著最暖季降水量變大,當(dāng)歸存在概率迅速提高,約為650mm時達到最大。隨后,隨著最暖季降水量變大,當(dāng)歸存在概率降低,當(dāng)最暖季降水量達到700mm時,當(dāng)歸存在概率降到0.4以下,隨著降水量的增大,存在概率略有增加,但小于0.5。以存在概率0.5以上為適宜范圍,因此,當(dāng)歸最暖季降水量的適宜范圍約為300—700 mm。

    表1 生態(tài)因子貢獻率

    (2)最冷季平均溫(貢獻率15.4%)在約-11℃以下時,當(dāng)歸存在概率很低,在0.1以下。隨著溫度升高,存在概率變大,約7℃時達到最大,隨后隨著溫度升高,當(dāng)歸存在概率降低,在11℃時降到0.1以下。根據(jù)存在概率0.5以上為適宜范圍,當(dāng)歸最冷季平均溫的適宜范圍約為-3—7 ℃。

    (3)降水的季節(jié)性(貢獻率8.6%)在約55以下時,當(dāng)歸存在概率低于0.3,隨著其上升,存在概率不斷變大,當(dāng)存在概率為0.5時,降水的季節(jié)性為70,降水的季節(jié)性為80時存在概率最大,然后不斷下降。因此存在概率為0.5以上時,當(dāng)歸分布區(qū)降水的季節(jié)性范圍約為70—98。

    圖3 主要氣候因子響應(yīng)曲線:最暖季降水量、最冷季平均溫和降水的季節(jié)性Fig.3 Response curves of three main climate factors: Precipitation of warmest quarter,Mean temperature of coldest quarter and Coefficient of variation of precipitation seasonality

    3.4 當(dāng)歸分布區(qū)土壤特征

    基于各土壤因子響應(yīng)曲線(圖4)可知適合當(dāng)歸種植的土壤因子范圍。

    圖4 主要土壤因子響應(yīng)曲線:碳酸鈣含量和粘土比重Fig.4 Response curves of two main soil factors:the calcium carbonate (lime) content in top- and subsoil and Percentage clay respectively in the top- and subsoil

    (1)碳酸鈣含量(9.5%),當(dāng)碳酸鈣含量約為55%時,當(dāng)歸的存在概率為0.5,隨著碳酸鈣含量增加當(dāng)歸的存在概率也不斷增加,當(dāng)碳酸鈣含量約為60%時,存在概率達到最大,當(dāng)碳酸鈣含量繼續(xù)增加時,存在概率反而先下降后上升然后下降,但是根據(jù)實際情況可知,碳酸鈣的含量范圍一般為0—1,因此碳酸鈣含量的適宜性范圍值約為55%—60%。

    (2)粘土比重(6.3%),當(dāng)粘土比重約為14%時,存在概率為0.3,粘土比重約為17%時,存在概率為0.5,隨著粘土比重加大,存在概率達到最大后變小,當(dāng)粘土比重變?yōu)?4%時,存在概率再次降至0.5,因此粘土比重的適宜范圍值約為17%—24%。

    3.5 當(dāng)歸生態(tài)區(qū)劃

    基于MaxEnt模型獲得的當(dāng)歸生態(tài)區(qū)劃結(jié)果,其存在概率取值在0到0.929之間,對其進行歸一化處理,將存在概率值歸一化到0到1之間,根據(jù)專家經(jīng)驗法確定70%以上為相似度最高等級,以20%間隔劃分出小于10%,10%至30%,30%至50%,50%至70%和大于70%,共5個等級。

    3.5.1 當(dāng)歸全球區(qū)劃

    圖5 當(dāng)歸全球生態(tài)區(qū)劃圖Fig.5 Global ecological division of Angelica sinensis

    當(dāng)歸潛在生態(tài)適宜區(qū)域主要分布在北半球北緯20—50°范圍內(nèi)的北美洲,歐洲,亞洲以及南半球南緯15—35°范圍內(nèi)的南美洲和非洲(圖5),共約593.07萬km2。北美洲的生態(tài)適宜區(qū)域主要位于美國中部的南達科他州、內(nèi)布拉斯加州、奧克拉荷馬州、堪薩斯州、德克薩斯州、新墨西哥洲等地以及美國西海岸的加利福尼亞州、華盛頓州、俄勒岡州;另外,墨西哥的奇瓦瓦州等地區(qū)也有少量分布。北美洲區(qū)域中當(dāng)歸的生態(tài)相似度主要在10%至30%之間,共約47.96萬km2。歐洲的生態(tài)適宜區(qū)域主要包括羅馬尼亞、烏克蘭西部、波蘭南部、俄羅斯西南部土耳其等地區(qū),生態(tài)相似度主要在10%至30%之間,共約27.22萬km2。南美洲的生態(tài)適宜區(qū)域主要分布在秘魯、智利、玻利維亞和阿根廷,生態(tài)相似度主要有10%至30%和30%至50%兩個等級,面積分別為9.01萬km2和2.05萬km2。非洲的生態(tài)適宜區(qū)域主要分布在贊比亞、津巴布韋、馬拉維科摩羅、南非和萊托索,生態(tài)相似度主要在10%到30%之間,共約28.23萬km2。亞洲的生態(tài)適宜區(qū)域主要分布在印度南部、阿富汗東部、巴基斯坦東南部、尼泊爾、不丹、中國、韓國、朝鮮、日本,其生態(tài)相似度在10%至30%,30%至50%,50%至70%,大于70% 4個等級均有分布,面積分別為207.58萬km2、132.73萬km2、82.45萬km2、44.31萬km2,共467.07萬km2,見表2。

    表2 當(dāng)歸不同生態(tài)相似度區(qū)域分布

    3.5.2 當(dāng)歸中國區(qū)劃

    利用ArcGIS 10.0對MaxEnt模型計算的結(jié)果與中國行政分區(qū)矢量圖進行疊加分析,得到當(dāng)歸在中國生態(tài)適宜區(qū)劃圖(圖6)。當(dāng)歸生態(tài)適宜區(qū)域在中國分布很廣泛,其生態(tài)相似度在10%—30%,30%—50%,50%—70%,>70% 4個等級均有分布,面積分別為138.2萬km2、114.03萬km2、77.56萬km2、41.08萬km2,共370.87萬km2(圖7)。其中生態(tài)相似度70%以上的區(qū)域主要集中在甘肅省南部、四川省、西藏自治區(qū)東部、云南省北部、貴州省、陜西省西南部??梢?中國十分適合當(dāng)歸的種植,是當(dāng)歸的主產(chǎn)區(qū)。

    圖6 當(dāng)歸中國生態(tài)區(qū)劃圖Fig.6 Ecological division of Angelica sinensis in China

    圖7 當(dāng)歸不同生態(tài)相似度區(qū)域中國面積分布圖Fig.7 Ecological similarity area of Angelica sinensis in China

    4 討論

    中藥材品質(zhì)和藥效與其生長的地域環(huán)境密切相關(guān)[24- 25],當(dāng)歸作為重要的中藥材之一,其產(chǎn)業(yè)發(fā)展一直受到密切關(guān)注。由于連作障礙等問題存在,當(dāng)歸的生長面積不斷縮減,迫切需要在世界范圍內(nèi)研究當(dāng)歸生態(tài)特征及其生態(tài)適宜區(qū),實現(xiàn)當(dāng)歸精細(xì)化種植管理及科學(xué)引種栽培,本文首次基于最大信息熵模型研究當(dāng)歸在全球范圍內(nèi)的適宜區(qū)及其生態(tài)特征。

    4.1 當(dāng)歸生態(tài)特征

    影響當(dāng)歸地理分布的生態(tài)因子主要有最暖季降水量、最冷季平均溫、碳酸鈣含量、降水的季節(jié)性、坡度、氣溫的季節(jié)性、粘土比重、海拔、土壤中粘粒組的陽離子交換能力和氣溫年范圍。最暖季降水量和最冷季平均溫對當(dāng)歸的地理分布貢獻率較大,最暖季是當(dāng)歸生長的旺盛期,表明當(dāng)歸生長期的降水量對當(dāng)歸的生長至關(guān)重要,最暖季降水量范圍的適宜值范圍為300—700 mm,這個降水量條件保證了水分含量但卻不會過分潮濕,這與目前對當(dāng)歸與道地產(chǎn)區(qū)環(huán)境適宜性的研究結(jié)果相一致,當(dāng)歸喜濕潤地帶,不耐旱不耐澇,水分含量過多易導(dǎo)致當(dāng)歸發(fā)生根腐病[26]。最冷季平均溫則是對當(dāng)歸的分布有著限制性影響,最冷季的極限溫度決定了當(dāng)歸是否能夠存活,本實驗結(jié)果表明當(dāng)歸分布區(qū)最冷季平均溫的適宜范圍為-3—7℃,而據(jù)相關(guān)文獻報道[27],適宜當(dāng)歸生長地區(qū)的年均溫約為2—6℃,研究結(jié)果和之前報道基本一致。

    影響當(dāng)歸地理分布的土壤因子主要有碳酸鈣含量、粘土比重。土壤碳酸鈣含量對當(dāng)歸生長有重要影響,表明當(dāng)歸對土壤的無機元素含量要求較高,相關(guān)文獻表明,當(dāng)歸原產(chǎn)區(qū)土壤中含有豐富的氮、磷、鉀、鎂等無機元素[28],因此在栽培過程中要注意土壤中元素含量和組成。粘土比重反映了土壤的物理性質(zhì),粘土比重小,土壤疏松透氣,排水性良好。由此可見,當(dāng)歸需要生活在無機元素含量豐富,疏松透氣排水性好的土壤中,這與文獻報道中對當(dāng)歸栽培區(qū)土壤研究結(jié)果相符合[29],因此在生態(tài)相似度相對低的地區(qū)引種栽培過程時,可以通過改善土壤質(zhì)地,合理施肥,以提高當(dāng)歸生長質(zhì)量。

    4.2 當(dāng)歸潛在分布區(qū)

    當(dāng)歸潛在分布區(qū)主要集中在北緯20—50°之間,緯度可能是限制其分布的一個因素,但在相同緯度的美洲地區(qū),當(dāng)歸分布極少,這可能是與不同地區(qū)海拔有關(guān),因為亞洲地區(qū)的當(dāng)歸主要是分布在海拔相對較高的地區(qū),這與當(dāng)歸喜低溫、高寒氣候的生物學(xué)特征相一致[30]。當(dāng)歸在美洲和歐洲有部分相似度為10%—30%的潛在分布區(qū)域,但目前關(guān)于這些地區(qū)的當(dāng)歸分布狀況很少有報導(dǎo),因此,在其他條件允許的情況下,這些地區(qū)可以成為當(dāng)歸合適的引種栽培區(qū)。分析結(jié)果顯示中國地區(qū)的四川省和甘肅省是當(dāng)歸的主產(chǎn)區(qū),在甘肅省南部和四川省的大片地區(qū)其生態(tài)相似度都大于70%,屬于生態(tài)相似度最高等級;同時,西藏自治區(qū)東部、云南省北部、貴州省、陜西省西南部等地區(qū)也有部分相似度大于70%的區(qū)域,另外,湖南、重慶、湖北等地區(qū)有生態(tài)相似度50%—70%的區(qū)域,目前這些地區(qū)雖有當(dāng)歸分布,但是規(guī)模較小,可以考慮擴種。此外,青海、河南、安徽、江蘇、浙江、江西、福建、湖南等地存在著大面積相似度為10%—30%的分布區(qū),因此,將這些地區(qū)的局部條件加以控制,可以會成為當(dāng)歸的生產(chǎn)基地,這與相關(guān)文獻中對當(dāng)歸適宜性分布研究結(jié)果[11]較一致。本文當(dāng)歸適宜區(qū)分析結(jié)果與第三次中藥資源普查結(jié)果相一致,當(dāng)歸資源生產(chǎn)區(qū)域均包含在本文生態(tài)相似度較高的潛在分布區(qū)域中,且潛在分布區(qū)域面積比第三次資源普查的當(dāng)歸真實分布區(qū)更大,因為模型預(yù)測物種的潛在分布,理論上大于物種的實際分布。同時,本文中不同地區(qū)的生態(tài)相似度與當(dāng)歸品質(zhì)評價結(jié)果相一致,即生態(tài)相似度較高的地區(qū)當(dāng)歸品質(zhì)較高[31- 32],因此,本文研究結(jié)果具有一定的科學(xué)指導(dǎo)意義。

    本文基于MaxEnt軟件和一定生態(tài)數(shù)據(jù)對全球當(dāng)歸地理分布進行預(yù)測,旨在為當(dāng)歸引種栽培提供科學(xué)參考依據(jù)。但由于影響藥用植物分布的因素有很多,且每種模型有其優(yōu)缺點,因此在擴種引種的具體操作中,為兼顧經(jīng)濟,種植地等因素,應(yīng)先進行試驗種植,方可大規(guī)模引種。

    [1] 嚴(yán)輝, 段金廒, 宋秉生, 王年鶴, 何子清, 孫成忠. 我國當(dāng)歸藥材生產(chǎn)現(xiàn)狀與分析. 中國現(xiàn)代中藥, 2009, 11(4): 12- 17.

    [2] 宋秋月, 付迎波, 劉江, 鄭丹, 韓力, 黃學(xué)石. 當(dāng)歸的化學(xué)成分研究. 中草藥, 2011, 42(10): 1900- 1904.

    [3] 李曦, 張麗宏, 王曉曉, 楊雯, 金玉青, 呂光華. 當(dāng)歸化學(xué)成分及藥理作用研究進展. 中藥材, 2013, 36(6): 1023- 1028.

    [4] Huang C Y, Kuo W W, Kuo C H, Tsai F J, Liu P Y, Hsieh D J Y. Protective effect ofDanggui(Radix Angelicae Sinensis) on angiotensin II-induced apoptosis in H9c2 cardiomyoblast cells. BMC Complementary and Alternative Medicine, 2014, 14: 358- 358.

    [5] Liu S P, Dong W G, Wu D F, Luo H S, Yu J P. Protective effect ofangelicasinensispolysaccharide on experimental immunological colon injury in rats. World Journal of Gastroenterology, 2003, 9(12): 2786- 2790.

    [6] 孫紅梅, 張本剛, 齊耀東, 張昭, 郭增祥. 當(dāng)歸藥材資源調(diào)查與分析. 中國農(nóng)學(xué)通報, 2009, 25(23): 437- 441.

    [7] 晉玲, 吳迪, 崔治家, 朱田田, 張弦飛, 金堯, 張裴斯. 當(dāng)歸藥材資源種類及分布研究. 中獸醫(yī)醫(yī)藥雜志, 2013, 32(1): 74- 77.

    [8] 張尚智, 賀莉萍, 韓黎明, 楊文璽, 劉玲玲. 野生當(dāng)歸及其當(dāng)歸植物資源的研究. 中國現(xiàn)代中藥, 2012, 14(4): 33- 36.

    [9] 吳寶華, 隋立軍, 鐘心堯, 酒連娣, 張友民. 當(dāng)歸屬植物資源及開發(fā)利用. 安徽農(nóng)業(yè)科學(xué), 2014, 42(3): 702- 703.

    [10] Engler R, Randin C F, Thuiller W, Dullinger S, Zimmermann N E, Araújo M B, Pearman P B, Le Lay G, Piedallu C, Albert C H, Choler P, Coldea G, De Lamo X, Dirnb?ck T, Gégout J C, Gómez-García D, Grytnes J A, Heegaard E, H?istad F, Nogués-Bravo D, Normand S, PucaM, Sebastià M T, Stanisci A, Theurillat J P, Trivedi M R, Vittoz P, Guisan A. 21st century climate change threatens mountain flora unequally across Europe. Global Change Biology, 2011, 17(7): 2330- 2341.

    [11] 邵璞, 曾曉東. 氣候年際變率對全球植被平均分布的影響. 生態(tài)學(xué)報, 2011, 31(6): 1494- 1505.

    [12] 馬松梅, 張明理, 張宏祥, 孟宏虎, 陳曦. 利用最大熵模型和規(guī)則集遺傳算法模型預(yù)測孑遺植物裸果木的潛在地理分布及格局. 植物生態(tài)學(xué)報, 2010, 34(11): 1327- 1335.

    [13] 劉恩斌, 周國模, 葛宏立. 基于最大熵原理的浙江毛竹胸徑分布及測量不確定度評定. 生態(tài)學(xué)報, 2009, 29(1): 86- 91.

    [14] Phillips S J, Anderson R P, Schapire R E. Maximum entropy modeling of species geographic distributions. Ecological Modelling, 2006, 190(3/4): 231- 259.

    [15] Phillips S J, Dudík M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography, 2008, 31(2): 161- 175.

    [16] Phillips S J, Dudík M, Schapire R E. A maximum entropy approach to species distribution modeling//Proceedings of the Twenty-First International Conference on Machine Learning. New York, NY, USA: ACM, 2004: 655- 662.

    [17] 胡忠俊, 張鐿鋰, 于海彬. 基于MaxEnt模型和GIS的青藏高原紫花針茅分布格局模擬. 應(yīng)用生態(tài)學(xué)報, 2015, 26(2): 505- 511.

    [18] Merow C, Silander J A Jr. A comparison of Maxlike and Maxent for modelling species distributions. Methods in Ecology and Evolution, 2014, 5(3): 215- 225.

    [19] 劉振生, 高惠, 滕麗微, 蘇云, 王曉勤, 孔芳毅. 基于MAXENT模型的賀蘭山巖羊生境適宜性評價. 生態(tài)學(xué)報, 2013, 33(22): 7243- 7249.

    [20] 王運生, 謝丙炎, 萬方浩, 肖啟明, 戴良英. ROC曲線分析在評價入侵物種分布模型中的應(yīng)用. 生物多樣性, 2007, 15(4): 365- 372.

    [21] 李振華, 王艷芳, 伊勒泰, 張艾華, 任凱, 張春紅, 畢雅瓊, 李旻輝. 阿拉善地區(qū)荒漠肉蓯蓉生態(tài)適宜性區(qū)劃研究. 中國中藥雜志, 2015, 40(5): 785- 792.

    [22] Wang Y S, Xie B Y, Wan F H, Xiao Q M, Dai L Y. Application of ROC curve analysis in evaluating the performance of alien species′ potential distribution models. Biodiversity Science, 2007, 15(4): 365- 372. (本條文獻為第20條文獻英文信息,請核對)

    [23] 麻亞鴻. 基于最大熵模型(MaxEnt)和地理信息系統(tǒng)(ArcGIS)預(yù)測蘚類植物的地理分布范圍——以廣西花坪自然保護區(qū)為例[D]. 上海: 上海師范大學(xué), 2013: 112- 113.

    [24] 王歡, 李慧, 曾凡琳, 謝彩香. 黃花蒿空間分布及全球潛在氣候適宜區(qū). 中藥材, 2015, 38(3): 460- 466.

    [25] 邵璽文, 韓梅, 韓忠明, 孔偉偉, 楊利民. 不同生境條件下黃芩光合日變化與環(huán)境因子的關(guān)系. 生態(tài)學(xué)報, 2009, 29(3): 1470- 1477.

    [26] 趙銳明, 陳垣, 郭鳳霞, 汪洋, 張尚智, 王瑞琪. 甘肅岷縣野生當(dāng)歸資源分布特點及其與栽培當(dāng)歸生長特性的比較研究. 草業(yè)學(xué)報, 2014, 23(2): 29- 37.

    [27] 趙慶芳, 曾小英, 馬瑞君. 當(dāng)歸栽培研究進展. 甘肅科技縱橫, 2002, 31(6): 37- 37.

    [28] 趙楊景, 陳四保, 高光耀, 馮毓秀, 楊世林, 徐麗珍, 杜力軍, 胡世林, 馮雪峰. 道地與非道地當(dāng)歸栽培土壤的理化性質(zhì). 中國中藥雜志, 2002, 27(1): 19- 22.

    [29] 姚玉璧, 馬鵬里, 張秀云. 道地與近道地當(dāng)歸栽培氣候生態(tài)與土壤環(huán)境區(qū)劃. 中國農(nóng)學(xué)通報, 2011,27(27): 156- 160.

    [30] 許世峰. 地膜當(dāng)歸栽培技術(shù)要點. 甘肅科技, 2013, 29(7): 139- 140.

    [31] 張宏意, 湯瑞平, 廖文波, 畢文鋼. 不同產(chǎn)地當(dāng)歸多糖含量的比較//廣東省植物學(xué)會學(xué)術(shù)研討會(廣東省植物學(xué)會2009年度青年學(xué)術(shù)研討會暨第四屆廣東植物學(xué)研究生學(xué)術(shù)交流會)論文集. 廣州: 中國植物學(xué)會,廣東省植物學(xué)會, 2009.

    [32] 楊樹聲. 中國當(dāng)歸藥材的質(zhì)量與安全性評價. 甘肅科技, 2015, 31(18): 74- 76.

    ResearchontheglobalecologicalsuitabilityandcharacteristicsofregionswithAngelicasinensisbasedontheMaxEntmodel

    ZHANG Dongfang1,2, ZHANG Qin1, GUO Jie1, SUN Chengzhong3,WU Jie4,NIE Xiang5,XIE Caixiang1,*

    1InstituteofMedicinalPlantDevelopment,PekingUnionMedicalCollege,ChineseAcademyofMedicalSciences,Beijing100193,China2JiangsuNormalUniversity,Xuzhou221116,China3ChineseAcademyofSurveying&Mapping,Beijing100039,China4InstituteofChineseMateriaChinaAcademyofChineseMedicalSciences,Beijing100029,China5MeteorologicalBureauofBiJieofGuiZhouProvince,Bijie551714,China

    In order to investigate the global ecological suitability and characteristics of regions for the medicinal plantAngelicasinensis, a maximum entropy model (MaxEnt) and geographic information system (GIS) technology were used to investigate global ecology suitability regions based on 109 sample points ofA.sinensisand 37 ecological factors. The results showed that global ecological regions that are suitable forA.sinensisare located mainly in the range of 20—50°N and 15—35°S, which includes sections of North America, Europe, Asia, South America and Africa. The total area of suitable regions found was about 5930700 km2. The regions with the most similarity were found to be mainly in the south of Gansu province, Sichuan, the east of the Tibet Autonomous Region, the north of Yunnan province, Guizhou, the southwest of Shanxi province, and other nearby regions. These results signify that China is very suitable for the cultivation and production ofA.sinensis. The climatic factors that were found to affect the distribution ofA.sinensiswere mainly precipitation during warmest quarter, mean temperature of the coldest quarter, and the coefficient of variation of precipitation seasonality, with the most suitable value ranges being 300—700 mm, -3—7℃ and 70—98 respectively. The soil factors affecting the distribution ofA.sinensiswere found to be mainly the calcium carbonate (lime) content and percentage clay, the value ranges of which were found to be 60%—0% and 17%—24%, respectively, in the ecologically suitable regions. These results provide a scientific reference for the refinement of the artificial cultivation ofA.sinensis.

    Angelicasinensis; maximum entropy model; geographic information system; ecological suitability; ecological characteristics

    國家自然科學(xué)基金項目(81473304,81130069);國家科技支撐計劃課題(2015BAI05B01);國家科技支撐計劃(2011BAI07B08);重大新藥創(chuàng)制國家科技重大專項(2014ZX09304307001)

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

    日期:2017- 03- 22

    10.5846/stxb201605030837

    *通訊作者Corresponding author.E-mail: caixiangxie@163.com

    張東方, 張琴, 郭杰,孫成忠,吳杰,聶祥,謝彩香.基于MaxEnt模型的當(dāng)歸全球生態(tài)適宜區(qū)和生態(tài)特征研究.生態(tài)學(xué)報,2017,37(15):5111- 5120.

    Zhang D F, Zhang Q, Guo J, Sun C Z,Wu J,Nie X,Xie C X.Research on the global ecological suitability and characteristics of regions withAngelicasinensisbased on the MaxEnt model.Acta Ecologica Sinica,2017,37(15):5111- 5120.

    猜你喜歡
    概率因子土壤
    第6講 “統(tǒng)計與概率”復(fù)習(xí)精講
    第6講 “統(tǒng)計與概率”復(fù)習(xí)精講
    因子von Neumann代數(shù)上的非線性ξ-Jordan*-三重可導(dǎo)映射
    土壤
    概率與統(tǒng)計(一)
    概率與統(tǒng)計(二)
    一些關(guān)于無窮多個素因子的問題
    靈感的土壤
    影響因子
    影響因子
    国产免费现黄频在线看| 亚洲情色 制服丝袜| 精品亚洲成国产av| 日韩大片免费观看网站| 高清视频免费观看一区二区| 婷婷色av中文字幕| 日本vs欧美在线观看视频| 成人综合一区亚洲| 国产伦精品一区二区三区视频9| 亚洲丝袜综合中文字幕| 午夜福利影视在线免费观看| 亚洲国产精品一区二区三区在线| 91成人精品电影| 综合色丁香网| 欧美日韩综合久久久久久| 婷婷成人精品国产| 久久ye,这里只有精品| 男女边摸边吃奶| 秋霞在线观看毛片| 五月玫瑰六月丁香| 天堂中文最新版在线下载| 精品亚洲成a人片在线观看| 免费av中文字幕在线| 18在线观看网站| 日韩精品免费视频一区二区三区 | 成人国产麻豆网| 色5月婷婷丁香| 亚洲精品第二区| 欧美少妇被猛烈插入视频| 精品人妻偷拍中文字幕| 在线观看免费视频网站a站| 国产无遮挡羞羞视频在线观看| 激情五月婷婷亚洲| 久久久久久久久久成人| 春色校园在线视频观看| 亚洲综合色惰| 成年人免费黄色播放视频| 99热6这里只有精品| 91久久精品国产一区二区三区| 国产永久视频网站| 狠狠精品人妻久久久久久综合| 美女福利国产在线| 国产极品天堂在线| 国产精品国产av在线观看| 97精品久久久久久久久久精品| 亚洲av.av天堂| 18+在线观看网站| 一本一本综合久久| 日韩精品有码人妻一区| 国产成人一区二区在线| 午夜激情福利司机影院| av天堂久久9| 亚洲精品久久久久久婷婷小说| 色5月婷婷丁香| 午夜精品国产一区二区电影| xxxhd国产人妻xxx| 久久久久久伊人网av| 黑丝袜美女国产一区| 另类亚洲欧美激情| 国产精品99久久99久久久不卡 | 日日摸夜夜添夜夜添av毛片| 国产亚洲精品久久久com| 久久久久人妻精品一区果冻| 精品久久久久久久久亚洲| 日韩欧美一区视频在线观看| av电影中文网址| 亚洲精品色激情综合| 三级国产精品欧美在线观看| 久久狼人影院| 久久99热6这里只有精品| 青春草亚洲视频在线观看| 国产极品粉嫩免费观看在线 | 丰满乱子伦码专区| 全区人妻精品视频| 国产高清三级在线| 欧美性感艳星| 只有这里有精品99| 乱码一卡2卡4卡精品| 最近的中文字幕免费完整| 亚洲内射少妇av| 热re99久久精品国产66热6| 一区二区三区精品91| 成人午夜精彩视频在线观看| 日韩 亚洲 欧美在线| 青春草亚洲视频在线观看| 国产探花极品一区二区| 丝袜脚勾引网站| 久久综合国产亚洲精品| 91aial.com中文字幕在线观看| 精品人妻熟女毛片av久久网站| av线在线观看网站| 丝袜美足系列| 日韩av免费高清视频| 69精品国产乱码久久久| 欧美+日韩+精品| 国产色婷婷99| 久久久久国产精品人妻一区二区| 日韩av免费高清视频| 91国产中文字幕| 中国美白少妇内射xxxbb| 狠狠精品人妻久久久久久综合| 日日爽夜夜爽网站| 日本-黄色视频高清免费观看| 免费黄色在线免费观看| 一区二区三区四区激情视频| 搡女人真爽免费视频火全软件| 美女脱内裤让男人舔精品视频| 久久久久久久亚洲中文字幕| 久久久久精品性色| 黄片无遮挡物在线观看| 在线精品无人区一区二区三| 麻豆精品久久久久久蜜桃| 久久99热这里只频精品6学生| 一本—道久久a久久精品蜜桃钙片| 69精品国产乱码久久久| 国产精品成人在线| 插逼视频在线观看| 久久精品国产亚洲av涩爱| 欧美97在线视频| 亚洲精品中文字幕在线视频| 国产欧美日韩综合在线一区二区| 午夜福利在线观看免费完整高清在| 另类亚洲欧美激情| 中文字幕人妻熟人妻熟丝袜美| a 毛片基地| 亚洲精品国产av成人精品| 免费观看性生交大片5| 少妇人妻 视频| 日韩一本色道免费dvd| 日韩在线高清观看一区二区三区| 国产亚洲一区二区精品| 男男h啪啪无遮挡| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品偷伦视频观看了| 欧美最新免费一区二区三区| 高清视频免费观看一区二区| 久久精品熟女亚洲av麻豆精品| 久久久久久久国产电影| 夜夜骑夜夜射夜夜干| a级毛片免费高清观看在线播放| 精品一区二区免费观看| 国产精品偷伦视频观看了| 中文字幕人妻丝袜制服| 久久精品国产亚洲av涩爱| 精品一区二区三卡| 18禁观看日本| 在线观看免费高清a一片| 简卡轻食公司| 黄色怎么调成土黄色| 高清午夜精品一区二区三区| 精品久久久久久久久av| 亚洲精品美女久久av网站| 99热这里只有是精品在线观看| 在线观看美女被高潮喷水网站| 色94色欧美一区二区| 国产色婷婷99| 99热这里只有精品一区| 日韩亚洲欧美综合| 好男人视频免费观看在线| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 国产精品.久久久| 精品99又大又爽又粗少妇毛片| 女性被躁到高潮视频| 一边亲一边摸免费视频| 精品国产乱码久久久久久小说| 欧美人与性动交α欧美精品济南到 | 一区二区日韩欧美中文字幕 | 视频中文字幕在线观看| 免费久久久久久久精品成人欧美视频 | 亚洲婷婷狠狠爱综合网| 国产精品女同一区二区软件| 欧美日韩视频高清一区二区三区二| av天堂久久9| 黄片播放在线免费| 国产av国产精品国产| 国产成人a∨麻豆精品| 日韩欧美精品免费久久| 女性被躁到高潮视频| 久久久久国产精品人妻一区二区| 亚洲欧美日韩另类电影网站| 一边摸一边做爽爽视频免费| 妹子高潮喷水视频| 久久久久久久大尺度免费视频| 久久久久久久久久久丰满| 亚洲少妇的诱惑av| 黄色欧美视频在线观看| 99久久精品一区二区三区| 精品人妻在线不人妻| 国产高清不卡午夜福利| a级毛片在线看网站| 久久久久久久久久久久大奶| 免费高清在线观看视频在线观看| 少妇被粗大猛烈的视频| 桃花免费在线播放| 国产男人的电影天堂91| 少妇的逼好多水| 在线观看国产h片| 久久午夜综合久久蜜桃| av国产精品久久久久影院| 18禁在线无遮挡免费观看视频| 美女福利国产在线| 亚洲欧美清纯卡通| 亚洲av免费高清在线观看| 曰老女人黄片| 性色av一级| 在线亚洲精品国产二区图片欧美 | 母亲3免费完整高清在线观看 | 又大又黄又爽视频免费| 日韩中字成人| 日本vs欧美在线观看视频| 高清av免费在线| 欧美另类一区| 国产成人a∨麻豆精品| 大香蕉97超碰在线| 免费日韩欧美在线观看| 亚洲精品第二区| 亚洲精品中文字幕在线视频| 国产一级毛片在线| 在线观看人妻少妇| 久久毛片免费看一区二区三区| 亚洲无线观看免费| 免费av不卡在线播放| 夜夜看夜夜爽夜夜摸| 91久久精品电影网| 亚洲精品自拍成人| 国产精品偷伦视频观看了| 在线观看免费日韩欧美大片 | 午夜激情av网站| 国产视频内射| 伊人久久国产一区二区| videos熟女内射| 久久毛片免费看一区二区三区| 男女啪啪激烈高潮av片| 91久久精品国产一区二区三区| 九色亚洲精品在线播放| 人妻一区二区av| av福利片在线| 99热网站在线观看| 一边亲一边摸免费视频| 老司机亚洲免费影院| 国产成人aa在线观看| 在线观看www视频免费| 欧美日韩精品成人综合77777| 久久99精品国语久久久| 国产极品粉嫩免费观看在线 | 色视频在线一区二区三区| 人人妻人人添人人爽欧美一区卜| 成人影院久久| 美女内射精品一级片tv| 国产亚洲精品久久久com| 97精品久久久久久久久久精品| 欧美 日韩 精品 国产| tube8黄色片| 免费人成在线观看视频色| 尾随美女入室| 狠狠婷婷综合久久久久久88av| 日韩一本色道免费dvd| 人人妻人人澡人人看| av一本久久久久| 国模一区二区三区四区视频| 大香蕉97超碰在线| a级毛色黄片| tube8黄色片| 在线精品无人区一区二区三| 免费高清在线观看日韩| 天天操日日干夜夜撸| 国产精品偷伦视频观看了| 99热网站在线观看| 精品久久国产蜜桃| 久久综合国产亚洲精品| 18禁在线无遮挡免费观看视频| 国产精品99久久99久久久不卡 | 国国产精品蜜臀av免费| 精品午夜福利在线看| 久久97久久精品| 在现免费观看毛片| 韩国高清视频一区二区三区| 18禁观看日本| av专区在线播放| 女性生殖器流出的白浆| 国产精品一国产av| xxx大片免费视频| 热re99久久国产66热| 性色av一级| 国产淫语在线视频| 久久精品人人爽人人爽视色| 22中文网久久字幕| 一区二区av电影网| 欧美性感艳星| 五月天丁香电影| 免费黄网站久久成人精品| 亚洲无线观看免费| 国产成人精品无人区| 日韩伦理黄色片| 久久ye,这里只有精品| 国内精品宾馆在线| 久久久久久久国产电影| 久热这里只有精品99| 大陆偷拍与自拍| 欧美人与性动交α欧美精品济南到 | 婷婷色麻豆天堂久久| 中国三级夫妇交换| 王馨瑶露胸无遮挡在线观看| 老女人水多毛片| 久久久亚洲精品成人影院| 国产男人的电影天堂91| 91精品伊人久久大香线蕉| 黄色视频在线播放观看不卡| 国产在线一区二区三区精| 十八禁网站网址无遮挡| 高清黄色对白视频在线免费看| 国产男女内射视频| 美女视频免费永久观看网站| 18+在线观看网站| 国产成人精品无人区| 各种免费的搞黄视频| 午夜福利视频精品| 国产精品国产三级国产专区5o| 国产免费又黄又爽又色| 大陆偷拍与自拍| 国产伦精品一区二区三区视频9| 日韩成人av中文字幕在线观看| 久久久久国产精品人妻一区二区| av在线老鸭窝| 国产高清不卡午夜福利| 草草在线视频免费看| 国产免费视频播放在线视频| 国产69精品久久久久777片| 午夜免费男女啪啪视频观看| 国产av精品麻豆| 欧美bdsm另类| 色视频在线一区二区三区| 国产精品国产三级国产专区5o| 蜜桃在线观看..| 日产精品乱码卡一卡2卡三| a级毛片免费高清观看在线播放| 国产老妇伦熟女老妇高清| 国产日韩欧美亚洲二区| 不卡视频在线观看欧美| 国产精品成人在线| 99九九线精品视频在线观看视频| 精品一区二区三区视频在线| 乱码一卡2卡4卡精品| 亚洲成色77777| 成人国产麻豆网| 国产精品女同一区二区软件| 亚洲精品久久午夜乱码| av卡一久久| 午夜91福利影院| 成年人午夜在线观看视频| 九色亚洲精品在线播放| 麻豆乱淫一区二区| 亚洲欧美日韩另类电影网站| 亚洲国产毛片av蜜桃av| 丝袜喷水一区| 秋霞伦理黄片| 日日摸夜夜添夜夜添av毛片| 欧美+日韩+精品| 老司机影院毛片| 高清黄色对白视频在线免费看| a级毛片在线看网站| 国产女主播在线喷水免费视频网站| 久久热精品热| 黄片无遮挡物在线观看| 少妇人妻 视频| 极品人妻少妇av视频| 久久人人爽人人爽人人片va| 日韩av免费高清视频| 观看美女的网站| 一级毛片我不卡| 亚洲精品国产av蜜桃| 精品酒店卫生间| 最黄视频免费看| 久久 成人 亚洲| 亚洲精品一区蜜桃| 国产白丝娇喘喷水9色精品| 美女主播在线视频| 国产免费视频播放在线视频| 亚洲av不卡在线观看| 岛国毛片在线播放| 亚洲欧洲国产日韩| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av天美| 久久国产精品大桥未久av| 中国三级夫妇交换| 亚洲欧美清纯卡通| 精品视频人人做人人爽| 男女啪啪激烈高潮av片| 亚洲国产精品一区三区| 亚洲国产欧美日韩在线播放| 久久久国产精品麻豆| 99久久中文字幕三级久久日本| 婷婷色综合www| 日韩精品免费视频一区二区三区 | 国产男女超爽视频在线观看| 99久久中文字幕三级久久日本| 人妻少妇偷人精品九色| 久久精品国产亚洲av涩爱| 久久影院123| 青春草国产在线视频| 国产不卡av网站在线观看| 国产精品偷伦视频观看了| 国产av国产精品国产| 国产亚洲精品第一综合不卡 | 麻豆成人av视频| 欧美日韩视频精品一区| 日日摸夜夜添夜夜爱| 中文字幕人妻熟人妻熟丝袜美| 全区人妻精品视频| 久久久久视频综合| 下体分泌物呈黄色| 久久亚洲国产成人精品v| 欧美日韩精品成人综合77777| 男男h啪啪无遮挡| 亚洲国产毛片av蜜桃av| 伊人亚洲综合成人网| 青春草亚洲视频在线观看| 亚洲国产欧美日韩在线播放| 国产精品.久久久| 免费不卡的大黄色大毛片视频在线观看| 女性被躁到高潮视频| h视频一区二区三区| 免费人妻精品一区二区三区视频| 美女xxoo啪啪120秒动态图| 女的被弄到高潮叫床怎么办| 99精国产麻豆久久婷婷| 中文精品一卡2卡3卡4更新| 下体分泌物呈黄色| 国产又色又爽无遮挡免| 国产熟女午夜一区二区三区 | 久久久久久久久久成人| 观看av在线不卡| 精品国产国语对白av| 国产日韩欧美在线精品| 91久久精品电影网| 成人漫画全彩无遮挡| 中文字幕免费在线视频6| 日本-黄色视频高清免费观看| 亚洲欧洲精品一区二区精品久久久 | 日韩成人伦理影院| 看非洲黑人一级黄片| 国产亚洲精品久久久com| 国产精品不卡视频一区二区| 一区二区av电影网| 国产精品嫩草影院av在线观看| 18+在线观看网站| 视频区图区小说| 人妻少妇偷人精品九色| 国产一区二区三区综合在线观看 | 免费观看a级毛片全部| 午夜91福利影院| 国产午夜精品一二区理论片| 国产欧美日韩综合在线一区二区| 日韩亚洲欧美综合| 亚洲精品美女久久av网站| 欧美精品国产亚洲| 如何舔出高潮| 满18在线观看网站| 熟妇人妻不卡中文字幕| 在现免费观看毛片| 免费观看av网站的网址| 丰满少妇做爰视频| 九草在线视频观看| 91在线精品国自产拍蜜月| 亚洲性久久影院| 日韩中文字幕视频在线看片| 少妇熟女欧美另类| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久| 香蕉精品网在线| 两个人的视频大全免费| 一区二区三区四区激情视频| 欧美一级a爱片免费观看看| 只有这里有精品99| 亚洲综合精品二区| 国产不卡av网站在线观看| 久久久久精品久久久久真实原创| av不卡在线播放| 国产精品99久久久久久久久| 中国三级夫妇交换| 亚洲精品中文字幕在线视频| 高清黄色对白视频在线免费看| 日本爱情动作片www.在线观看| 亚洲国产毛片av蜜桃av| 一区在线观看完整版| 亚洲精品国产av成人精品| a级片在线免费高清观看视频| 一区二区av电影网| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 青春草亚洲视频在线观看| 丁香六月天网| 日本av免费视频播放| 99久久精品国产国产毛片| 亚洲国产日韩一区二区| 国产精品三级大全| 97在线人人人人妻| 日产精品乱码卡一卡2卡三| 18禁在线无遮挡免费观看视频| 少妇被粗大的猛进出69影院 | 国产一区二区三区综合在线观看 | 一本色道久久久久久精品综合| 韩国高清视频一区二区三区| 在线看a的网站| 欧美日韩国产mv在线观看视频| 日韩欧美一区视频在线观看| 我的女老师完整版在线观看| 国产乱来视频区| 在线看a的网站| 日韩欧美精品免费久久| 亚洲在久久综合| 大片电影免费在线观看免费| 熟妇人妻不卡中文字幕| 国产精品偷伦视频观看了| 色婷婷久久久亚洲欧美| 国产成人精品在线电影| 人人妻人人添人人爽欧美一区卜| 国产av一区二区精品久久| 菩萨蛮人人尽说江南好唐韦庄| 国产精品一二三区在线看| 中文字幕人妻丝袜制服| 91久久精品国产一区二区成人| 国产不卡av网站在线观看| 国产亚洲欧美精品永久| 新久久久久国产一级毛片| 69精品国产乱码久久久| 久久99热6这里只有精品| 欧美精品一区二区大全| 午夜福利影视在线免费观看| 女人久久www免费人成看片| 青春草亚洲视频在线观看| 日韩精品免费视频一区二区三区 | 五月开心婷婷网| 黑人欧美特级aaaaaa片| 人妻 亚洲 视频| 国产一区亚洲一区在线观看| 免费黄网站久久成人精品| 日本av免费视频播放| 一级黄片播放器| 久久免费观看电影| 国产精品熟女久久久久浪| 九九在线视频观看精品| 久久久久久久精品精品| 色5月婷婷丁香| 狂野欧美激情性xxxx在线观看| 久久国产精品男人的天堂亚洲 | 久久久久精品久久久久真实原创| 亚洲精品中文字幕在线视频| 日韩av不卡免费在线播放| 亚洲精品日本国产第一区| 夫妻性生交免费视频一级片| 最近中文字幕高清免费大全6| 免费看不卡的av| 免费av中文字幕在线| 日本欧美视频一区| 亚洲av男天堂| 纯流量卡能插随身wifi吗| 哪个播放器可以免费观看大片| 99热这里只有是精品在线观看| 22中文网久久字幕| 久久国内精品自在自线图片| 久久国产精品大桥未久av| 卡戴珊不雅视频在线播放| 亚洲在久久综合| 成人漫画全彩无遮挡| 在线播放无遮挡| 中文字幕最新亚洲高清| h视频一区二区三区| 啦啦啦在线观看免费高清www| 国产乱人偷精品视频| 又大又黄又爽视频免费| 香蕉精品网在线| 男人操女人黄网站| 一边摸一边做爽爽视频免费| 大话2 男鬼变身卡| 性色avwww在线观看| 99精国产麻豆久久婷婷| 嫩草影院入口| 国产成人一区二区在线| 久久久久人妻精品一区果冻| 蜜桃在线观看..| 国产高清国产精品国产三级| 精品人妻在线不人妻| 国产亚洲午夜精品一区二区久久| 亚洲三级黄色毛片| 精品久久蜜臀av无| 一级爰片在线观看| 十分钟在线观看高清视频www| 中文字幕人妻熟人妻熟丝袜美| 中国三级夫妇交换| 午夜影院在线不卡| 中文精品一卡2卡3卡4更新| 久久久久久久亚洲中文字幕| xxx大片免费视频| 国产高清三级在线| 亚洲精品色激情综合| 91成人精品电影| xxx大片免费视频| 精品一区二区三卡| 男女边吃奶边做爰视频| 蜜桃在线观看..| 国产日韩一区二区三区精品不卡 | 精品国产露脸久久av麻豆| 中文字幕精品免费在线观看视频 | 你懂的网址亚洲精品在线观看| 高清欧美精品videossex| 天堂中文最新版在线下载| 日韩亚洲欧美综合| av网站免费在线观看视频| 精品一区二区三区视频在线| 国产永久视频网站| 少妇熟女欧美另类| 欧美日韩在线观看h| 18禁动态无遮挡网站|