• <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)于無窮多個素因子的問題
    靈感的土壤
    影響因子
    影響因子
    eeuss影院久久| 麻豆av噜噜一区二区三区| 久久99热这里只有精品18| 天天躁夜夜躁狠狠久久av| www.av在线官网国产| 国产午夜精品论理片| 亚洲美女搞黄在线观看| 久久久久性生活片| 91久久精品国产一区二区成人| 免费不卡的大黄色大毛片视频在线观看 | 亚洲色图av天堂| 91精品一卡2卡3卡4卡| 国产色爽女视频免费观看| 欧美成人精品欧美一级黄| 午夜精品在线福利| 免费搜索国产男女视频| 日韩强制内射视频| 日日撸夜夜添| 26uuu在线亚洲综合色| 日日啪夜夜撸| 国产精品久久久久久久电影| 少妇熟女aⅴ在线视频| 亚洲人成网站在线观看播放| 国产综合懂色| 99热这里只有精品一区| 国产视频内射| 国产老妇伦熟女老妇高清| 一本久久中文字幕| 综合色丁香网| 亚洲欧洲日产国产| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品一区www在线观看| 免费电影在线观看免费观看| 亚洲成人中文字幕在线播放| 我的老师免费观看完整版| 熟妇人妻久久中文字幕3abv| 国产综合懂色| 麻豆av噜噜一区二区三区| 成人二区视频| 热99re8久久精品国产| 岛国在线免费视频观看| 天堂√8在线中文| 亚洲欧美成人综合另类久久久 | 波野结衣二区三区在线| 国产精品人妻久久久久久| 色尼玛亚洲综合影院| 国产精品嫩草影院av在线观看| 成人漫画全彩无遮挡| 久久精品国产清高在天天线| 国产爱豆传媒在线观看| 亚洲色图av天堂| 日韩成人伦理影院| 欧美潮喷喷水| 高清在线视频一区二区三区 | 亚洲无线在线观看| 91麻豆精品激情在线观看国产| 此物有八面人人有两片| 国产熟女欧美一区二区| 国产探花极品一区二区| 久久久久久伊人网av| 国产精品综合久久久久久久免费| 免费av不卡在线播放| av在线亚洲专区| 精品国内亚洲2022精品成人| 一本精品99久久精品77| 国产v大片淫在线免费观看| 51国产日韩欧美| 在线观看午夜福利视频| 超碰av人人做人人爽久久| 国产亚洲91精品色在线| 日日撸夜夜添| 亚洲欧洲国产日韩| 日韩欧美精品v在线| 日本五十路高清| 国产高潮美女av| 日本黄色视频三级网站网址| 一个人看视频在线观看www免费| 少妇的逼水好多| 国产色爽女视频免费观看| 直男gayav资源| 女的被弄到高潮叫床怎么办| 乱系列少妇在线播放| 插阴视频在线观看视频| 毛片一级片免费看久久久久| 日本爱情动作片www.在线观看| 久久久久久大精品| av又黄又爽大尺度在线免费看 | 日韩国内少妇激情av| a级一级毛片免费在线观看| 亚洲欧美日韩高清在线视频| 老女人水多毛片| 在线播放国产精品三级| 人妻久久中文字幕网| 中出人妻视频一区二区| 啦啦啦啦在线视频资源| 麻豆乱淫一区二区| 乱系列少妇在线播放| 在线免费十八禁| 日本欧美国产在线视频| 春色校园在线视频观看| 在线观看午夜福利视频| 美女内射精品一级片tv| 久久韩国三级中文字幕| 嫩草影院新地址| a级毛片免费高清观看在线播放| 国产精品三级大全| 亚洲精品色激情综合| 在线免费观看的www视频| 99久久精品热视频| 日韩精品有码人妻一区| 亚洲精品日韩av片在线观看| 国产亚洲5aaaaa淫片| av在线观看视频网站免费| 亚洲精品乱码久久久久久按摩| 国产高潮美女av| av在线蜜桃| 两个人的视频大全免费| 身体一侧抽搐| 国产精品久久久久久久电影| 网址你懂的国产日韩在线| 少妇熟女aⅴ在线视频| 不卡视频在线观看欧美| 国产精品99久久久久久久久| 成人永久免费在线观看视频| 国产精品不卡视频一区二区| 亚洲性久久影院| 成人美女网站在线观看视频| 色视频www国产| 乱系列少妇在线播放| 好男人视频免费观看在线| 最近手机中文字幕大全| 麻豆成人午夜福利视频| 国产黄色小视频在线观看| 亚洲自偷自拍三级| 日韩中字成人| 波多野结衣高清无吗| 国产女主播在线喷水免费视频网站 | 又黄又爽又刺激的免费视频.| 人妻久久中文字幕网| 国产麻豆成人av免费视频| 91久久精品国产一区二区成人| 国产亚洲91精品色在线| 91av网一区二区| ponron亚洲| 成年女人看的毛片在线观看| 看十八女毛片水多多多| 搞女人的毛片| 久久久久久久久中文| 黄色欧美视频在线观看| 特大巨黑吊av在线直播| 非洲黑人性xxxx精品又粗又长| 日韩亚洲欧美综合| 一级毛片久久久久久久久女| 精华霜和精华液先用哪个| АⅤ资源中文在线天堂| 观看美女的网站| 亚洲最大成人手机在线| a级毛色黄片| 亚洲av不卡在线观看| 国产 一区 欧美 日韩| 99热精品在线国产| 尾随美女入室| 国产毛片a区久久久久| 国产高清不卡午夜福利| av在线蜜桃| 亚洲乱码一区二区免费版| 在线观看一区二区三区| 日韩欧美三级三区| 久久久久久久午夜电影| а√天堂www在线а√下载| 18禁裸乳无遮挡免费网站照片| 精品人妻熟女av久视频| 久久精品久久久久久久性| 蜜桃久久精品国产亚洲av| 午夜精品一区二区三区免费看| 在线国产一区二区在线| 熟妇人妻久久中文字幕3abv| 亚洲成人av在线免费| 欧美精品一区二区大全| 哪个播放器可以免费观看大片| 久久99热这里只有精品18| 天堂√8在线中文| 91在线精品国自产拍蜜月| 国产亚洲精品久久久com| 干丝袜人妻中文字幕| 亚洲图色成人| 国产激情偷乱视频一区二区| 熟妇人妻久久中文字幕3abv| 三级经典国产精品| 波多野结衣高清无吗| 国产精品一及| 国产精品精品国产色婷婷| 色播亚洲综合网| 久久久国产成人免费| 99热这里只有精品一区| 亚洲av中文字字幕乱码综合| 三级经典国产精品| 中文欧美无线码| 亚洲国产欧美在线一区| 色播亚洲综合网| 久久久色成人| 色视频www国产| 国产在视频线在精品| 免费人成在线观看视频色| 国产午夜精品论理片| .国产精品久久| 久久精品人妻少妇| 精华霜和精华液先用哪个| 99久久中文字幕三级久久日本| 欧美性猛交黑人性爽| 亚洲aⅴ乱码一区二区在线播放| 国产精华一区二区三区| 麻豆精品久久久久久蜜桃| 内地一区二区视频在线| 美女被艹到高潮喷水动态| 亚洲最大成人中文| 在线天堂最新版资源| 国产乱人视频| 午夜a级毛片| 亚洲内射少妇av| а√天堂www在线а√下载| 最近手机中文字幕大全| 午夜激情福利司机影院| 欧洲精品卡2卡3卡4卡5卡区| 免费人成在线观看视频色| 久久久久国产网址| 日本免费一区二区三区高清不卡| 黄片无遮挡物在线观看| 国产精品,欧美在线| 最近视频中文字幕2019在线8| 久久午夜福利片| 热99re8久久精品国产| 亚洲国产精品国产精品| kizo精华| 中文亚洲av片在线观看爽| 久久久久久久久大av| 国产精品日韩av在线免费观看| 亚洲精品粉嫩美女一区| 夜夜爽天天搞| 午夜激情福利司机影院| 99久久人妻综合| 久久这里有精品视频免费| 一个人看的www免费观看视频| 日本免费一区二区三区高清不卡| 看非洲黑人一级黄片| 九色成人免费人妻av| 岛国在线免费视频观看| 精品人妻熟女av久视频| 亚洲不卡免费看| 亚洲丝袜综合中文字幕| 国产精品国产高清国产av| 中文字幕av成人在线电影| 免费av毛片视频| 日韩成人av中文字幕在线观看| 久久午夜亚洲精品久久| 91精品国产九色| 看片在线看免费视频| 欧美日韩综合久久久久久| 欧美日本亚洲视频在线播放| 亚洲自偷自拍三级| 69av精品久久久久久| 少妇人妻精品综合一区二区 | 最近最新中文字幕大全电影3| 欧美人与善性xxx| 亚洲av不卡在线观看| 成人高潮视频无遮挡免费网站| 哪个播放器可以免费观看大片| 亚洲性久久影院| 黄色日韩在线| 久久99蜜桃精品久久| 国产高清视频在线观看网站| 亚洲av一区综合| 麻豆成人午夜福利视频| 天堂影院成人在线观看| 国产伦一二天堂av在线观看| 美女脱内裤让男人舔精品视频 | 有码 亚洲区| 一级毛片我不卡| 精品日产1卡2卡| 日韩在线高清观看一区二区三区| 全区人妻精品视频| 草草在线视频免费看| 超碰av人人做人人爽久久| 熟女人妻精品中文字幕| 最近最新中文字幕大全电影3| 岛国毛片在线播放| 精品人妻偷拍中文字幕| 黑人高潮一二区| 午夜福利在线观看吧| 特大巨黑吊av在线直播| 欧美高清成人免费视频www| 国产精品一区二区性色av| 精品国产三级普通话版| 一级毛片久久久久久久久女| 麻豆成人午夜福利视频| av在线老鸭窝| 人人妻人人澡人人爽人人夜夜 | 天美传媒精品一区二区| 成年免费大片在线观看| av专区在线播放| 寂寞人妻少妇视频99o| 久久这里只有精品中国| 亚洲欧美日韩卡通动漫| 久久久国产成人精品二区| 国产在线男女| 在线观看av片永久免费下载| 精品午夜福利在线看| 好男人视频免费观看在线| 国产精品一二三区在线看| 国产探花在线观看一区二区| 精品一区二区三区人妻视频| av免费在线看不卡| 毛片一级片免费看久久久久| 欧美3d第一页| av在线亚洲专区| 22中文网久久字幕| 亚洲国产欧洲综合997久久,| 插阴视频在线观看视频| 国产av一区在线观看免费| 亚洲av二区三区四区| 国产一区二区在线观看日韩| 亚洲精品影视一区二区三区av| 中文字幕久久专区| 国产黄片美女视频| 我的老师免费观看完整版| 欧美bdsm另类| 精品久久国产蜜桃| 色播亚洲综合网| 麻豆精品久久久久久蜜桃| 精品久久久久久久久久免费视频| 少妇熟女aⅴ在线视频| 大型黄色视频在线免费观看| 搡女人真爽免费视频火全软件| 99在线人妻在线中文字幕| 久久6这里有精品| 久久人人爽人人片av| 国产综合懂色| 两性午夜刺激爽爽歪歪视频在线观看| 成人毛片60女人毛片免费| 18禁黄网站禁片免费观看直播| 国产精品三级大全| 禁无遮挡网站| 中文字幕免费在线视频6| kizo精华| 午夜福利高清视频| 色综合亚洲欧美另类图片| 一级二级三级毛片免费看| 国产精品久久久久久精品电影小说 | 九九爱精品视频在线观看| 久久久久久久久久成人| 成人毛片60女人毛片免费| 久久久成人免费电影| 欧美性感艳星| 最近2019中文字幕mv第一页| 国产在视频线在精品| 亚洲国产精品合色在线| 国产精品福利在线免费观看| 亚洲内射少妇av| 给我免费播放毛片高清在线观看| 亚洲图色成人| 久久鲁丝午夜福利片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲七黄色美女视频| 午夜精品国产一区二区电影 | 夜夜夜夜夜久久久久| 内射极品少妇av片p| av免费观看日本| 少妇猛男粗大的猛烈进出视频 | 晚上一个人看的免费电影| 日本熟妇午夜| 亚洲欧洲国产日韩| 亚洲va在线va天堂va国产| 天堂网av新在线| 一级毛片aaaaaa免费看小| 国产女主播在线喷水免费视频网站 | 69人妻影院| 精品99又大又爽又粗少妇毛片| 久久久久国产网址| 国产伦精品一区二区三区视频9| 久久午夜福利片| 国产在线精品亚洲第一网站| 悠悠久久av| 久久久久久久午夜电影| 啦啦啦啦在线视频资源| 亚洲色图av天堂| 亚洲一级一片aⅴ在线观看| 人妻少妇偷人精品九色| 91精品一卡2卡3卡4卡| av免费观看日本| 天天躁日日操中文字幕| 99久久精品一区二区三区| 日产精品乱码卡一卡2卡三| 日韩av在线大香蕉| 精品久久久久久成人av| 日本免费一区二区三区高清不卡| 男女视频在线观看网站免费| 国产av麻豆久久久久久久| 少妇被粗大猛烈的视频| 久99久视频精品免费| 成人亚洲欧美一区二区av| 亚洲内射少妇av| 亚洲国产精品久久男人天堂| 99视频精品全部免费 在线| 国产精品一区二区三区四区久久| 日韩 亚洲 欧美在线| 精品少妇黑人巨大在线播放 | 欧美日韩在线观看h| 天天躁日日操中文字幕| 久久久精品大字幕| 神马国产精品三级电影在线观看| 69人妻影院| 日韩欧美国产在线观看| 一进一出抽搐gif免费好疼| 日韩强制内射视频| 午夜福利成人在线免费观看| 精品不卡国产一区二区三区| 亚洲18禁久久av| 日本撒尿小便嘘嘘汇集6| 色综合色国产| av又黄又爽大尺度在线免费看 | 色哟哟·www| 男人狂女人下面高潮的视频| 熟女电影av网| 99在线人妻在线中文字幕| 国产v大片淫在线免费观看| 精品不卡国产一区二区三区| 精品久久国产蜜桃| 一级av片app| 国产精品一区二区在线观看99 | 国内精品宾馆在线| 18+在线观看网站| 亚洲欧美成人综合另类久久久 | 高清毛片免费观看视频网站| 免费人成在线观看视频色| 日韩欧美精品v在线| 一级毛片我不卡| 免费一级毛片在线播放高清视频| 亚洲乱码一区二区免费版| 18禁在线无遮挡免费观看视频| 亚洲丝袜综合中文字幕| 18禁黄网站禁片免费观看直播| 午夜免费激情av| 久久亚洲国产成人精品v| 色综合亚洲欧美另类图片| 国产黄色视频一区二区在线观看 | 欧美日韩国产亚洲二区| 国产成人91sexporn| 欧美xxxx性猛交bbbb| 亚洲成av人片在线播放无| av免费观看日本| 99riav亚洲国产免费| 啦啦啦韩国在线观看视频| 亚洲真实伦在线观看| 欧美zozozo另类| 国产高清三级在线| 成人特级黄色片久久久久久久| 日日撸夜夜添| 熟女人妻精品中文字幕| 国产伦理片在线播放av一区 | 亚洲内射少妇av| 干丝袜人妻中文字幕| 在线a可以看的网站| 麻豆国产97在线/欧美| 九九爱精品视频在线观看| 九九久久精品国产亚洲av麻豆| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 亚洲成人av在线免费| 青春草亚洲视频在线观看| 精品人妻一区二区三区麻豆| 亚洲乱码一区二区免费版| 亚洲四区av| 亚洲国产精品合色在线| 亚洲内射少妇av| 欧美高清性xxxxhd video| 欧美+亚洲+日韩+国产| 在线观看午夜福利视频| 中文资源天堂在线| 丰满人妻一区二区三区视频av| 婷婷色av中文字幕| 91久久精品电影网| 一个人观看的视频www高清免费观看| 少妇的逼水好多| 美女黄网站色视频| 国产一级毛片在线| 天天躁夜夜躁狠狠久久av| 亚洲婷婷狠狠爱综合网| 亚洲无线观看免费| 中文精品一卡2卡3卡4更新| 欧美另类亚洲清纯唯美| 日本免费一区二区三区高清不卡| 嫩草影院入口| 国产精品不卡视频一区二区| 亚洲18禁久久av| 精品久久久久久成人av| 波多野结衣高清无吗| 国内精品宾馆在线| 国内精品一区二区在线观看| 亚洲欧洲日产国产| 别揉我奶头 嗯啊视频| 三级国产精品欧美在线观看| 久久午夜福利片| 伊人久久精品亚洲午夜| 深夜精品福利| 精品久久久久久久久久久久久| 亚洲乱码一区二区免费版| 亚洲成a人片在线一区二区| 欧美一区二区精品小视频在线| 国产亚洲精品久久久com| 午夜免费男女啪啪视频观看| 国产国拍精品亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 26uuu在线亚洲综合色| 久久精品国产99精品国产亚洲性色| 1000部很黄的大片| 99久国产av精品国产电影| 波野结衣二区三区在线| 国产成人a∨麻豆精品| 亚洲欧美日韩东京热| 欧美性猛交╳xxx乱大交人| 少妇的逼水好多| 亚洲精品粉嫩美女一区| 亚洲在线观看片| 亚洲美女视频黄频| 一级av片app| 日韩精品青青久久久久久| 国内精品宾馆在线| av福利片在线观看| 女人十人毛片免费观看3o分钟| 久久亚洲国产成人精品v| 精品久久久久久久久av| 国产成人a区在线观看| 丰满的人妻完整版| 晚上一个人看的免费电影| 国产av在哪里看| 日韩 亚洲 欧美在线| 亚洲一区二区三区色噜噜| 一级毛片我不卡| 久久久久久久久久久免费av| 一个人免费在线观看电影| 欧美高清成人免费视频www| 国产国拍精品亚洲av在线观看| 国产 一区 欧美 日韩| 菩萨蛮人人尽说江南好唐韦庄 | 日韩强制内射视频| 国产激情偷乱视频一区二区| 国产高清视频在线观看网站| 国内久久婷婷六月综合欲色啪| 国产不卡一卡二| 蜜桃亚洲精品一区二区三区| 亚洲欧美精品专区久久| 人妻系列 视频| 波多野结衣高清无吗| 亚洲欧美中文字幕日韩二区| 你懂的网址亚洲精品在线观看 | 在线观看美女被高潮喷水网站| 欧洲精品卡2卡3卡4卡5卡区| 久久热精品热| 国产精品一区www在线观看| 人妻制服诱惑在线中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | a级毛片免费高清观看在线播放| 免费人成在线观看视频色| 波多野结衣高清作品| 久久久久久久亚洲中文字幕| av免费在线看不卡| 亚洲天堂国产精品一区在线| 国产精品av视频在线免费观看| 精品99又大又爽又粗少妇毛片| 哪个播放器可以免费观看大片| 嘟嘟电影网在线观看| 18禁黄网站禁片免费观看直播| 又爽又黄无遮挡网站| 亚洲欧美日韩无卡精品| 女人被狂操c到高潮| 色哟哟哟哟哟哟| 国产片特级美女逼逼视频| av在线蜜桃| 99在线视频只有这里精品首页| 悠悠久久av| 国产精品一区www在线观看| 午夜福利高清视频| 国产真实乱freesex| 一进一出抽搐gif免费好疼| 午夜爱爱视频在线播放| 国产乱人视频| 久久人人爽人人片av| 蜜臀久久99精品久久宅男| 国产精品一二三区在线看| 中文欧美无线码| 中出人妻视频一区二区| eeuss影院久久| 一卡2卡三卡四卡精品乱码亚洲| 五月伊人婷婷丁香| 99热这里只有精品一区| 九九久久精品国产亚洲av麻豆| av黄色大香蕉| 日韩三级伦理在线观看| 联通29元200g的流量卡| 国产精品日韩av在线免费观看| 人人妻人人看人人澡| 欧美性猛交╳xxx乱大交人| 成人永久免费在线观看视频| 嫩草影院新地址| 欧美不卡视频在线免费观看| 好男人视频免费观看在线| 亚洲人成网站在线播放欧美日韩| 国产精华一区二区三区| 1000部很黄的大片| 美女 人体艺术 gogo| 亚洲欧美精品专区久久| 如何舔出高潮| 国内揄拍国产精品人妻在线| 精品人妻一区二区三区麻豆| 国产精品人妻久久久影院|