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

    基于層次分析法和GIS的秦嶺地區(qū)魔芋潛在分布研究

    2015-02-05 01:23:30衛(wèi)海燕郭彥龍
    生態(tài)學報 2015年21期
    關鍵詞:種植區(qū)魔芋秦嶺

    高 蓓, 衛(wèi)海燕,郭彥龍, 顧 蔚

    1 陜西師范大學旅游與環(huán)境學院,西安 710062 2 陜西師范大學生命科學學院,西安 710062 3 陜西省農(nóng)業(yè)遙感信息中心, 西安 710015

    基于層次分析法和GIS的秦嶺地區(qū)魔芋潛在分布研究

    高 蓓1,3, 衛(wèi)海燕1,*,郭彥龍1, 顧 蔚2

    1 陜西師范大學旅游與環(huán)境學院,西安 710062 2 陜西師范大學生命科學學院,西安 710062 3 陜西省農(nóng)業(yè)遙感信息中心, 西安 710015

    魔芋(Amorphophallusrivieri)為我國傳統(tǒng)的食用與藥用植物,主產(chǎn)于秦嶺以南地區(qū)。依據(jù)陜西秦嶺地區(qū)28個魔芋分布點信息,利用秦嶺地區(qū)45 個氣象臺站1961—2010 年氣候資料、土壤數(shù)據(jù)(分辨率1 km)和DEM高程數(shù)據(jù)(分辨率30 m),結(jié)合前人研究,通過魔芋產(chǎn)量與環(huán)境指標的相關性分析,獲取相關性顯著的20 個評價指標,包括氣候指標13 個、土壤指標4 個和地形指標3 個,運用GIS技術(shù)和多元回歸模型對氣候指標進行柵格化處理,基于層次分析法和加權(quán)平均法獲得評價指標權(quán)重,建立陜西秦嶺地區(qū)魔芋潛在種植分布模型,確定魔芋潛在種植的空間分布。結(jié)果顯示:陜西秦嶺地區(qū)魔芋最適宜種植區(qū)面積1 214.42 km2,占可種植區(qū)面積的10.18%;適宜種植區(qū)面積2 015.60 km2,占可種植區(qū)面積的16.90%;次適宜種植區(qū)面積3 115.03 km2,占可種植區(qū)面積的26.12%;不適宜種植區(qū)面積5 580.02 km2占可種植區(qū)面積的46.80%。適宜魔芋潛在種植區(qū)域主要分布在陜西漢中中南部、安康中南部以及商洛東南部。

    層次分析法;GIS;潛在種植分布;魔芋;陜西秦嶺地區(qū)

    魔芋(Amorphophallusrivieri)又名蒟蒻,主產(chǎn)于秦嶺以南的山區(qū)或高原地區(qū),是我國傳統(tǒng)食用與藥用植物[1-2]。陜西秦嶺地區(qū)是我國魔芋種植的傳統(tǒng)區(qū)域,是中國南北氣候的分界線,也是生態(tài)自然環(huán)境保存較好的區(qū)域[3- 5]??茖W的劃分農(nóng)作物種植區(qū)對該區(qū)域生態(tài)平衡的維持、生態(tài)系統(tǒng)穩(wěn)定性的維護、可持續(xù)發(fā)展中生態(tài)與經(jīng)濟關系的協(xié)調(diào)都有重要影響。目前國內(nèi)外對魔芋研究主要集中于種質(zhì)資源[6-7]、遺傳多樣性和基因挖掘[8]、育種[9]、栽培[10]、病害[11-12]、食品[13-14]和藥用價值等方面[15- 18],隨著人們對魔芋功效的逐漸了解,市場對魔芋的廣泛需求,推動了魔芋種植面積的不斷擴大,經(jīng)濟效益愈加明顯,魔芋種植成為耕地緊張、耕作條件差的地區(qū)農(nóng)民脫貧致富的一條路徑[19]。由于缺少魔芋生存環(huán)境適宜性空間分布信息,導致魔芋產(chǎn)業(yè)盲目發(fā)展,不利于產(chǎn)業(yè)可持續(xù)發(fā)展。

    物種生存環(huán)境適宜性研究是物種保護的重要環(huán)節(jié),對物種資源利用有著重要意義[20]。層次分析法(AHP)是一種定性和定量相結(jié)合、系統(tǒng)化、層次化的分析方法,為多目標、多準則以及內(nèi)容復雜問題提供簡便實用的決策方法[21]。由于具有處理復雜決策問題的實用性和有效性特點,該方法已在不同領域得到應用[22-23],為物種適宜性研究提供了借鑒。

    對魔芋的空間分布研究主要集中在氣候指標對魔芋生長的影響、氣候種植區(qū)劃以及魔芋種質(zhì)資源分布等方面[24- 26],其對魔芋的適宜性研究多以數(shù)理統(tǒng)計為主,或以行政區(qū)劃為基本評價單元,精細化程度不夠,而根據(jù)魔芋產(chǎn)量和生長期數(shù)據(jù),綜合考慮氣候、土壤和地形的影響,來探討魔芋種植適宜性評價和合理種植區(qū)化研究未見報道。本文利用陜西秦嶺地區(qū)DEM數(shù)據(jù)、土壤數(shù)據(jù)和氣象臺站資料,以魔芋生長密切相關的生態(tài)要素作為評價指標,采用層次分析法和加權(quán)平均法[27],應用ArcGIS軟件,對陜西秦嶺地區(qū)魔芋潛在種植分布進行區(qū)劃,為客觀調(diào)控和指導魔芋適宜性規(guī)模種植提供基礎資料,為當?shù)胤N植規(guī)劃和經(jīng)濟發(fā)展提供參考依據(jù)。

    1 數(shù)據(jù)來源與研究方法

    1.1 研究對象

    魔芋性喜溫涼,生長于海拔250—2500 m山間,一般于3月中旬—4月上旬發(fā)芽,10月中旬—11月收獲。生長條件為年平均氣溫11 ℃以上,≥10 ℃積溫2900—8000 ℃·d,7—8月月平均氣溫12—26 ℃、月平均最高氣溫15—35 ℃,塊莖膨大期氣象條件要求為7—9月日溫差≥10 ℃,年降雨量800 mm以上,6—9月月降雨量100—200 mm,7—8月月平均相對濕度76%—90%,無霜期220 d以上[28]。影響魔芋生長的主要氣象災害是春季倒春寒、夏季干旱與延遲性冷害和秋季連陰雨。

    適宜魔芋生長的主要土壤條件為土層深度達30 cm以上,有機質(zhì)含量豐富、土層深厚、土壤疏松肥沃的壤土,種植適宜的pH值6.5—7.0[24]。適宜魔芋生長的地形特征為濕度較高的稍傾斜的背風地帶,具有灌溉條件且夏季暴雨不致土壤嚴重沖刷的地段[28]。

    1.2 數(shù)據(jù)來源

    本文的魔芋分布點位置與環(huán)境數(shù)據(jù)來源于西南大學中國魔芋中心“魔芋種質(zhì)資源的調(diào)查研究和保存”項目、《中國魔芋種植區(qū)劃》[26]和中國數(shù)字植物標本館(http://www.cvh.org.cn/)以及陜西省安康市魔芋產(chǎn)業(yè)技術(shù)研究中心的陜西省魔芋調(diào)查報告及實驗數(shù)據(jù),其中從中國數(shù)字植物標本館獲得研究區(qū)20個分布點信息,由相關文獻獲得研究區(qū)12個分布點信息,由“魔芋種質(zhì)資源的調(diào)查研究和保存”項目獲得研究區(qū)8個分布點信息,由陜西省魔芋調(diào)查報告及實驗數(shù)據(jù)獲得研究區(qū)22個分布點61個樣方數(shù)據(jù),剔除重復點數(shù)據(jù),共獲得陜西秦嶺地區(qū)28個分布點67個樣方數(shù)據(jù)。陜西省秦嶺地區(qū)魔芋產(chǎn)量面積數(shù)據(jù)來源于陜西農(nóng)業(yè)廳。

    生態(tài)環(huán)境適宜性評價因子的選擇是評價的關鍵[29],依據(jù)魔芋產(chǎn)量與環(huán)境條件指標的相關分析以及魔芋生長特點,選擇了13個氣候、4個土壤和3個地形因子構(gòu)建陜西秦嶺地區(qū)魔芋潛在種植分布的環(huán)境評價指標(表1),其中氣象數(shù)據(jù)來自陜西省46個氣象站1961—2010年逐日觀測數(shù)據(jù),土壤數(shù)據(jù)和地形數(shù)據(jù)來自中科院國際科學數(shù)據(jù)服務平臺(ISDSP,http://datamirror.csdb,cn),土壤柵格數(shù)據(jù)的分辨率為1 km,高程數(shù)字模型(DEM)的分辨率為30 m。陜西省1∶10萬土地利用數(shù)據(jù)來自于中國西部環(huán)境與生態(tài)科學數(shù)據(jù)中心(http://westdc.westgis.ac.cn)。

    表1 陜西秦嶺地區(qū)魔芋潛在種植分布模型的評價指標Table 1 Evaluation factors of potential distribution of A.rivieri in Qinling mountains of Shaanxi province

    1.3 研究方法

    1.3.1 候數(shù)據(jù)小網(wǎng)格化

    參考山地無氣象站地區(qū)氣候要素的推算模擬方法,選用回歸統(tǒng)計方法和對比分析方法[30]。利用研究區(qū)46個氣象站的經(jīng)度、緯度和海拔作為自變量,SPSS進行多元回歸分析,建立氣象要素空間分布模型(表2)。運用ArcGIS 9.3的空間插值模塊,通過對比Simple kriging、Original kriging、IDW和RBF 4種插值方法,選擇能夠較好反映局地氣候垂直變化的IDW法進行空間插值。由于受地形起伏變化的影響,各氣候指標的總體擬合精度不夠高,為提高擬合精度,對氣候指標的殘差部分進行空間內(nèi)插,以訂正氣候指標。將氣候指標多元回歸分析柵格圖和殘差的柵格分布圖進行圖像疊加,最終生成氣侯指標空間分布柵格圖像。

    1.3.2 評價指標的確定

    影響魔芋生長和分布的因素很多,依據(jù)陜西秦嶺地區(qū)各魔芋種植縣的魔芋產(chǎn)量、魔芋生長期觀測數(shù)據(jù)以及當?shù)厣鷳B(tài)環(huán)境資料,采用定性與定量結(jié)合的方法,遵循客觀性、區(qū)域差異性、主導性和可操作性等原則進行指標的篩選[31]。

    表2 陜西秦嶺地區(qū)氣候指標的空間推算模型Table 2 The small grid reckon model of climate factors in Qinling mountains of Shaanxi provice

    降雨量(mm)、氣溫(℃)、積溫(℃·d)、無霜期(d)和相對濕度(%);**表示通過了0.01的顯著性水平檢驗

    圖1 陜西秦嶺地區(qū)魔芋潛在種植分布研究結(jié)構(gòu)Fig.1 Structure of A.rivieri Potential distribution evaluation in Qinling mountains of Shaanxi province

    1.3.3 綜合評價體系構(gòu)建

    構(gòu)建各評價指標權(quán)重[27],以確定魔芋潛在種植分布的權(quán)重評價體系,即影響魔芋生長的3個準則層評價因素(氣候、土壤和地形),20個評價指標(表1)對魔芋潛在種植分布的影響作用是不同的,這就需要根據(jù)每個評價指標影響作用的重要程度分別賦予不同的權(quán)重(表3)。邀請12位在魔芋栽培、病害、種質(zhì)資源、魔芋產(chǎn)品開發(fā)和土壤領域有豐富知識經(jīng)驗的專家和技術(shù)人員參加問卷填寫,將課題研究的目的告訴專家,請專家根據(jù)多年的實踐經(jīng)驗填寫調(diào)查問卷。首先分別在各層中評價因子之間作兩兩重要性定量判斷比較;其次根據(jù)判斷矩陣計算上一層次指標與本層次相關因素之間的重要性次序的權(quán)重值(權(quán)向量);然后進行層次單排序和層次總排序,并對AHP構(gòu)造的判斷矩陣進行一致性檢驗,即CI(Consistencyindex)與判斷矩陣的平均隨機一致性指標RI(Randomindex)的比值CR(Consistencyratio)進行一致性檢驗[27]。

    經(jīng)計算得出CI=0.06,RI=1.19,CR=CI/RI=0.049<0.10,層次排序的結(jié)果具有滿意的一致性。

    表3 魔芋潛在種植分布研究評價指標的AHP權(quán)重Table 3 Index weights of analytic hierarchy process of the function of potential distribution suitability of A.rivieri

    1.3.4 基于GIS綜合評價模型的建立

    采用加權(quán)平均法,通過模糊隸屬函數(shù)標準化的指標值與層次分析法求出的每個評價指標的權(quán)重值[27,32],利用公式(1)計算每個評價柵格的綜合分。

    (1)

    式中,S為魔芋潛在種植分布適宜度;χi為二級指標的標準化值;ωi為二級指標權(quán)重;ωj為級指標權(quán)重;m為一級指標的個數(shù);n為二級指標的個數(shù)。

    1.3.5 模型檢驗

    為了驗證模型的適用性,選擇陜西秦嶺地區(qū)采樣地的6 個樣本,利用均方根誤差[32](RMSE),檢驗魔芋種植區(qū)劃評價模型的適用性。

    (2)

    式中,N為樣本數(shù),Pi是第i個預測值,Qi代表第i個觀察值。均方根誤差表示的預測值與實測值相應的接近程度,RMSE越小,模型預測值與實際值越接近,模型效果越好。

    本模型的RMSE檢驗結(jié)果為7.8%,RMSE值小于10%,表明模擬值與實際觀測值一致性很高,模擬效果可以接受[29],證明該模型適用于陜西秦嶺地區(qū)魔芋潛在空間分布研究。同時,對秦嶺外的6個采樣地數(shù)據(jù)進行檢驗,RMSE值為6.1%,表明模型也適用于研究陜西秦嶺地區(qū)以外的魔芋潛在空間分布。

    2 結(jié)果

    圖2 陜西秦嶺地區(qū)魔芋潛在種植分布適宜性區(qū)劃 Fig.2 Potential distribution of evaluation on A.ecosuitability in Qinling mountains of Shaanxi provice

    采用分布指標疊加法來進行區(qū)域劃分,運用ArcGIS空間分析模塊,將20 個評價指標按照魔芋潛在種植適宜性評價指標的評分標準(表4)賦予相應的值。利用GIS綜合評價模型(公式1),計算出每個柵格的潛在種植適宜性綜合評價值,并將柵格分值分為4個等級:潛在種植適宜性評價指數(shù)(S)6.00—7.31為最適宜區(qū),5.00—6.00為適宜區(qū),4.00—5.00為次適宜區(qū),2.49—4.00為不適宜區(qū)??紤]到魔芋種植有一定的局限性,參考陜西省2004年土地利用現(xiàn)狀圖,在魔芋適宜性區(qū)劃時將居民區(qū)、水體、林地、灌叢和沙地等不適合種植區(qū)域剔除掉,得到陜西秦嶺地區(qū)魔芋潛在種植適宜性區(qū)劃圖和空間分布圖(圖2和表5)。該圖符合實際,可以為政府部門決策和農(nóng)民種植生產(chǎn)提供參考。

    表4 魔芋潛在種植分布評價指標的評分標準Table 4 Appraisal criterion of all division factors of potential suitability of A.rivieri

    表5 陜西秦嶺地區(qū)魔芋潛在種植適宜性分布等級和面積Table 5 Area and potential distribution of evaluation on A.ecosuitability in Qinling mountains of Shaanxi provice

    陜西省秦嶺地區(qū)魔芋最適宜種植面積為1 214.42 km2,占陜西省秦嶺地區(qū)可種植區(qū)面積的10.18%,主要分布在漢中的勉縣東南部、漢臺區(qū)中部、城固中部、洋縣中西部、寧強中部、西鄉(xiāng)中部局地、鎮(zhèn)巴,安康的紫陽、嵐皋、石泉南部,商洛的鎮(zhèn)安東部、山陽、柞水東部、商州中部。這些地區(qū)海拔700—1 100 m,熱量充足,年平均氣溫均在14.0 ℃以上,年降水量為700—1 100 mm,7—8月平均相對濕度為80%左右,7—8月最高平均氣溫為29.0—30.5 ℃,7—9月氣溫日溫差≥10 ℃日數(shù)在38 d以上。該區(qū)域生態(tài)環(huán)境條件較優(yōu),較適于魔芋生長,也是目前陜西魔芋的主要種植區(qū),最適宜建立規(guī)模經(jīng)營的魔芋基地,應在目前種植基礎上逐漸擴大規(guī)模。

    陜西秦嶺地區(qū)魔芋適宜種植面積有2 015.60 km2,占可種植區(qū)面積的16.90%,主要分布在漢中的西鄉(xiāng)中部、漢臺區(qū)大部、勉縣南部、寧強大部,安康的漢濱區(qū)、漢陰、旬陽、平利,商洛的鎮(zhèn)安大部、商州區(qū)、洛南,西安的藍田、周至,寶雞的眉縣,渭南的華縣、華陰。這些地區(qū)海拔400—700 m中低山區(qū),區(qū)域土壤主要是沙壤土和砂質(zhì)土,光熱資源豐富,雨水較充沛,生態(tài)適宜性條件可以滿足魔芋種植的要求,主要表現(xiàn)為絕大多數(shù)評價指標符合魔芋種植的生長條件,但有個別區(qū)劃指標不符合要求,主要表現(xiàn)夏季高溫干旱時有發(fā)生,這些地區(qū)可作為陜西秦嶺地區(qū)魔芋種植的推廣發(fā)展區(qū),建議在這些區(qū)域,一方面培育適合該區(qū)域種植的品種或改變種植條件,另一方面在專家和技術(shù)人員的指導下,有選擇的種植魔芋,以增加當?shù)剞r(nóng)業(yè)產(chǎn)業(yè)發(fā)展,增加農(nóng)民收入。

    陜西秦嶺地區(qū)魔芋次適宜種植面積為3 115.03 km2,占可種植區(qū)面積的26.12%,主要分布在寶雞北部,渭南的渭城區(qū),安康和商洛也有一部分。這些區(qū)域氣候溫和,降雨較多,但是雨熱分布不均,不能滿足魔芋整個生長期的要求,即使種植也僅是能夠成活,不能促進地下塊莖的生長,不推薦這些地區(qū)種植魔芋。而魔芋不適宜種植面積為5 580.02 km2,占可種植區(qū)面積的46.80%,主要分布在西安、渭南、咸陽和寶雞東部等地區(qū)。這些地區(qū)為海拔500 m以下的平原地區(qū),年平均氣溫均在11—13 ℃,年降水量500—650 mm,7—8月平均相對濕度70—75%,7—9月氣溫日溫差≥10 ℃日數(shù)在22—30 d,冬季和春季易發(fā)生干旱,夏季高溫頻繁發(fā)生,秋季連陰雨也時有發(fā)生,這些地區(qū)不宜發(fā)展魔芋產(chǎn)業(yè)。

    3 討論

    3.1 模型建立適用性

    本文針對魔芋生長條件的要求,用AHP對魔芋潛在種植適宜性分布評價指標進行選擇,綜合考慮氣候、土壤和地形要素,選擇20個評價指標,運用ArcGIS軟件,結(jié)合基礎地理數(shù)據(jù)和土壤數(shù)據(jù),同時考慮到無氣象站觀測點對潛在分布準確性的影響,對基礎氣象觀測數(shù)據(jù)進行數(shù)字化處理,建立氣象數(shù)據(jù)模型,最終生成分辨率為100 m的柵格數(shù)據(jù),建立了陜西秦嶺地區(qū)魔芋潛在種植分布適宜性區(qū)劃模型,模擬準確達到好的程度(RMSE值小于10%)。

    AHP構(gòu)造的判斷矩陣一致性檢驗結(jié)果表明,一級評價指標層和二級評價指標層中各具體參評指標對于目標層——陜西秦嶺地區(qū)魔芋潛在分布區(qū)研究的相對重要性的CR檢驗顯示CR=0.049<0.1,表明層次結(jié)構(gòu)的判斷矩陣的一致性是可以接受的。層次分析模型的構(gòu)造是合理的,可以運用于魔芋潛在種植區(qū)分布。

    3.2 魔芋的潛在分布

    植物分布的適宜生態(tài)環(huán)境首先要為植物生長提供適當?shù)乃疅釛l件,即適當?shù)臏囟群退謼l件,Trisuratet等[33]研究表明植物生境中春季和秋季的水熱條件是影響其生長的主導氣候因子,本文中魔芋潛在種植分布適宜性評價體系權(quán)重(表3)表明,氣候要素權(quán)重為0.57、土壤要素權(quán)重為0.26,地形要素權(quán)重為0.17,說明在氣候要素對魔芋種植和產(chǎn)量影響作用大,土壤要素次之,地形指標對于魔芋生長的作用在大尺度上有限。對于二級評價指標層,無霜期(0.0798)、年平均氣溫(0.0741)、土壤pH值(0.0728)和耕層厚度(0.0702)所占權(quán)重較大,其中無霜期決定魔芋的生長季長短,作為喜溫作物,當?shù)孛鏈囟冉抵? ℃或0 ℃以下時,魔芋就會遭受霜凍的危害;年平均氣溫表征種植區(qū)的基本熱量情況,在秦嶺地區(qū)氣溫與降水密切相關,魔芋是喜溫涼陰濕的物種,氣溫和降水是控制著其生長的主要指標。土壤酸堿度是土壤的重要的基本性質(zhì)之一,在實際種植過程中魔芋必需的營養(yǎng)元素的有效性與土壤的pH值有關;耕層厚度在一定程度上決定了土壤可供養(yǎng)分的情況,是決定魔芋產(chǎn)量的重要因素。

    魔芋適宜生長在溫涼陰濕的環(huán)境中,對氣溫降水有一定要求(表4)。春季,只要氣溫能夠達到閾值(10 ℃左右),魔芋能萌芽[2]。魔芋種植最適宜的年無霜期為大于260 d,年平均氣溫為14—20 ℃,≥10 ℃積溫為≥4000 ℃,7—8月平均最高氣溫為20—30 ℃。水分條件是影響魔芋潛在種植分布的另一項要素,最適宜魔芋種植的年降雨量是1 200—1 800 mm,7—8月平均相對濕度是80%—90%,所占權(quán)重分別為0.0570和0.0399,表明水分在魔芋生長過程中是制約其分布的指標,分布在較濕潤的環(huán)境。魔芋最適宜的耕層厚度是40—60 cm,土壤類型為沙壤土,土壤pH值6.5—7.2,土壤有機質(zhì)含量≥20.0 kg/m3。在相同緯度上,海拔為700—1 000 m,坡度為3—10°半陰坡是魔芋生長的最佳地形。

    表6 魔芋不同等級潛在種植區(qū)面積Table 6 The areas of different levels of potential distribution of evaluation on A.ecosuitability in different cities

    本研究構(gòu)建的模型模擬出陜西秦嶺地區(qū)魔芋潛在適宜區(qū),與調(diào)查資料結(jié)合魔芋生長的生態(tài)條件繪制的中國魔芋分布[23]基本一致(圖2和表5)。最適宜區(qū)和適宜區(qū)面積為3 230.02 km2,占陜西秦嶺地可種植區(qū)面積的27.09%,主要分布在漢中、安康、商洛等地,分別為687.36、612.70 km2和463.12 km2(表6),這些地區(qū)是陜西魔芋種植的傳統(tǒng)產(chǎn)區(qū),因此對魔芋種植和野生資源的保護應有目的選擇最適宜區(qū)或適宜區(qū)加以管理,同時注意病害對魔芋生產(chǎn)的影響,特別在魔芋育種選址和建設時,應充分考慮其適宜區(qū),優(yōu)選最適宜區(qū)。

    4 結(jié)論

    本文依據(jù)生態(tài)環(huán)境數(shù)據(jù),運用AHP和ArcGIS空間分析技術(shù),分析了魔芋在陜西秦嶺地區(qū)的潛在種植分布和主要影響指標,模擬出魔芋的潛在種植地理分布。研究結(jié)果顯示魔芋在研究區(qū)適宜生長(最適宜和適宜區(qū))的面積僅為3 230.02 km2,占研究區(qū)總面積的3.99%,其生長區(qū)域狹窄,對生態(tài)環(huán)境要求較嚴格。同時由于全球氣候變化、受人類活動影響以及魔芋自身病害特點(自然條件下種植1a需停種2—3a的種植方式),導致魔芋的分布區(qū)面積減小,確定適宜種植區(qū)對魔芋原產(chǎn)地保護具有重要意義。本文綜合區(qū)劃的陜西秦嶺地區(qū)魔芋潛在種植空間分布,為魔芋生產(chǎn)基地的選址和建設提供數(shù)據(jù)支撐,為宏觀調(diào)控和指導魔芋的規(guī)模化種植提供參考依據(jù)。

    [1] 龐杰, 張盛林, 劉佩瑛, 張興國. 中國魔芋資源的研究. 資源科學, 2001, 23(5): 87- 89.

    [2] 劉佩英. 魔芋學. 北京: 中國農(nóng)業(yè)出版社, 2004: 15- 22.

    [3] 白紅英, 馬新萍, 高翔, 侯欽磊. 基于DEM的秦嶺山地1月氣溫及0℃等溫線變化. 地理學報, 2012, 67(11): 1143- 1450.

    [4] 康慕誼, 朱源. 秦嶺山地生態(tài)分界線的論證. 生態(tài)學報, 2007, 27(7): 2274- 2284.

    [5] 王昌海, 溫亞利, 李強, 胡崇德, 司開創(chuàng). 秦嶺自然保護區(qū)群的社會效益計量研究. 中國人口·資源與環(huán)境, 2011, 21(7): 113- 121.

    [6] 牛義, 張盛林, 王志敏, 李川, 孫遠航. 中國魔芋資源的研究與利用. 西南農(nóng)業(yè)大學學報: 自然科學版, 2005, 27(5): 634- 638.

    [7] 宣慢. 魔芋種質(zhì)資源形態(tài)多樣性與ISSR分析 [D]. 重慶: 西南農(nóng)業(yè)大學, 2011: 1- 17.

    [8] 任盤宇, 潘明清. 云南南部5種魔芋屬植物居群遺傳結(jié)構(gòu)的ISSR分析. 武漢大學學報: 理學版, 2013, 59(1): 99- 104.

    [9] 牛義, 張大學, 劉海利, 王啟軍, 劉紅艷, 張盛林. 魔芋新品種渝魔1號的選育. 中國蔬菜, 2010, (2): 88- 90.

    [10] 張紅驥, 邵梅, 杜鵬, 盧俊, 何霞紅, 于德才, 朱有勇. 云南省魔芋與玉米多樣性栽培控制魔芋軟腐病. 生態(tài)學雜志, 2012, 31(2): 332- 336.

    [11] 雷珍珍, 葉晶龍, 程海麗, 陳云, 望彗星, 許克靜, 樂超銀. 花魔芋抗軟腐病植株的鑒定及其抗性機理的初步研究. 植物學報, 2013, 48(3): 295- 302.

    [12] 何斐, 段佳麗, 羅寶芳, 張忠良, 劉列平, 張小明, 薛泉宏. 嵐皋魔芋軟腐病病原細菌生物多樣性研究. 西北農(nóng)林科技大學學報: 自然科學版, 2013, 41(7): 91- 98.

    [13] 劉艷, 郭華春, 張雅瓊, 李婉琳. 魔芋與玉米間作群體中魔芋植株生長及葡苷聚糖含量變化的研究. 西南農(nóng)業(yè)學報, 2013, 26(3): 1120- 1125.

    [14] 彭述輝, 龐杰. 魔芋葡甘聚糖在環(huán)保、醫(yī)藥及食品中的應用. 長江蔬菜, 2006, (12): 30- 33.

    [15] Kraemer W J, Vingren J L, Silvestre R, Spiering B A, Hatfield D L, Ho J Y, Fragala M S, Maresh C M, Volek J S. Effect of adding exercise to a diet containing glucomannan. Metabolism Clinical and Experimental, 2007, 56(8): 1149- 1158.

    [16] Keithley J, Swanson B. Glucomannan and obesity: a critical review. Alternatives Therapies in Health and Medicine, 2005, 11(6): 30- 34.

    [17] Alonso S M, Teijeiro O D, Remunán L C, Alonso M J. Glucomannan, a promising polysaccharide for biopharmaceutical purposes. European Journal of Pharmaceutics and Biopharmaceutics, 2009, 72(2): 453- 462.

    [18] Onishi N, Kawamoto S, Ueda K, Yamanaka Y, Katayama A, Suzuki H, Aki T, Hashimoto K, Hide M, Ono K. Dietary pulverized konjac glucomannan prevents the development of allergic rhinitis- like symptoms and IgE response in mice, Bioscience, Biotechnology and Biochemistry, 2007, 71(10): 2551- 2556.

    [19] 謝松峰, 崔鳴. 陜南貧困山區(qū)魔芋產(chǎn)業(yè)發(fā)展現(xiàn)狀和問題及對策建議. 農(nóng)業(yè)現(xiàn)代化研究, 2012, 33(6): 270- 278.

    [20] 朱耿平, 劉國卿, 卜文俊, 高玉葆. 生態(tài)位模型的基本原理及其在生物多樣性保護中的應用. 生物多樣性, 2013, 21(1): 90- 98.

    [21] 陳海生, 劉國順, 劉大雙, 陳偉強. GIS支持下的河南省煙草生態(tài)適宜性綜合評價. 中國農(nóng)業(yè)科學, 2009, 42(7): 2425- 2433.

    [22] 王青, 戴思蘭, 何晶, 季玉山, 王朔. 灰色關聯(lián)法和層次分析法在盆栽多頭小菊株系選擇中的應用. 中國農(nóng)業(yè)科學, 2012, 45(17): 3653- 3660.

    [23] 翟治芬, 王蘭英, 孫敏章, 嚴昌榮, 梁瑜珊. 基于AHP與Rough Set的農(nóng)業(yè)節(jié)水技術(shù)綜合評價. 生態(tài)學報, 2012, 32(3): 931- 941.

    [24] 童碧慶, 吳俊銘. 貴州山地氣候條件對魔芋生長的影響及種植區(qū)劃. 中國農(nóng)業(yè)氣象, 2004, 25(4): 35- 38.

    [25] 朱伯華, 楊玲, 謝國生, 李晏斌, 周洪彪, 費甫華. 湖北省魔芋種植區(qū)劃的系統(tǒng)分析. 湖北農(nóng)業(yè)科學, 2005, (3): 11- 14.

    [26] 楊代明, 劉佩瑛. 中國魔芋種植區(qū)劃. 西南農(nóng)業(yè)大學學報: 自然科學版, 1990, 12(1): 1- 7.

    [27] Lu C Y, Gu W, Dai A H, Wei H Y. Assessing habitat suitability based on geographic information system (GIS) and fuzzy: a case study ofSchisandrasphenantheraRehd. et Wils. In Qinling Mountains, China. Ecological Modelling, 2012, 242: 105- 115.

    [28] 朱義, 張盛林. 魔芋防病與高效栽培技術(shù). 武漢: 中國三峽出版社, 2008: 13- 30.

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

    [30] 傅抱璞, 虞靜明, 盧其堯. 山地氣候資源與開發(fā)利用. 南京: 南京大學出版社, 1996: 87- 103.

    [31] 沈漢, 鄒國元. 菜地土壤評價中參評因素的選定與分級指標的劃分. 土壤通報, 2004, 35(5): 553- 557.

    [32] 郭彥龍, 顧蔚, 路春燕, 衛(wèi)海燕. 基于五味子甲素與乙素的秦嶺地區(qū)野生華中五味子的潛在空間分布. 植物學報, 2013, 48(4): 411- 422.

    [33] Trisurat Y, Shrestha R P, Kjelgren R. Plant species vulnerability to climate change in Peninsular Thailand. Applied Geography, 2011, 31(3): 1106- 1114.

    Potential distribution ofAmorphophallusrivieriin the Qinling Mountains based on analytic hierarchy process and geographic information system

    GAO Bei1,3, WEI Haiyan1,*, GUO Yanlong1, GU Wei2

    1CollegeofTourismandEnvironment,ShaanxiNormalUniversity,Xi′an710062,China

    2CollegeofLifeSciences,ShaanxiNormalUniversity,Xi′an710062,China

    3RemoteSensingInformationCenterforAgricultureofShaanxiProvince,Xi′an710015,China

    Amorphophallusrivieri(the corpse flower) is a traditional edible and medicinal plant in China. This species is distributed in the south part of the Qinling Mountains, China. We assimilated data aboutA.riviericultivation, environmental information from 28 sampling sites in the Qinling Mountains, climate data from 45 weather stations in the Qinling Mountains from 1961 to 2010, soil data with 1 km × 1 km spatial resolution and DEM data with 30 m × 30 m spatial resolution in the Qinling Mountains,A.rivieridata collected throughout China, and a specific report onA.rivieriin Shaanxi Province. We obtained 20 assessment factors that were significantly correlated when evaluatingA.rivieriyield against environmental factors. The key environmental factors affecting the distribution ofA.riviericultivation included 13 dominant climate factors, 4 dominant soil factors, and 3 dominant topographical factors. These dominant factors were 1) Frost-free duration (D), 2) Annual average temperature (Tn), 3) Annual total active temperature(≥10°C)(T≥10djw), 4) Monthly mean maximum temperature from July to August (T78zg), 5) Annual precipitation (Pn), 6) Monthly mean daily temperature range from July to September (T79gc), 7) Monthly mean temperature from May to October (T510p), 8) Monthly mean temperature from July to August (T78p), 9) Monthly mean relative air humidity from July to August (Q78), 10) June precipitation (P6), 11) July precipitation (P7), 12) August precipitation (P8),13) September precipitation (P9), 14) Topsoil depth (H), 15) Topsoil pH(H2O)(pH), 16) topsoil texture classification (C), 17) Topsoil organic matter (O), 18) Aspect (A), 19) Slope (S), 20) Altitude elevation (h). Using Geographic Information System (GIS) and a multivariate regression model, the climate factors were rasterized. Then, we used fuzzy mathematics analysis, analytic hierarchy process (AHP), and the weighted means method to set up the subjection function and determine the weight of each factor. We set up a model of ecological suitability forA.rivieriin the Qinling Mountains of Shaanxi Province, and determined the spatial distribution of suitable planting areas for this species. The root-mean-square error (RMSE) was used to evaluate the accuracy of the model predictions. The RMSE value reached 7.8%, which indicated that the predictive accuracy of the model was “Excellent.” The results showed that the ecological planting suitability model identified a relationship between potentialA.riviericultivation distribution and the environmental factors. Highly suitable, moderately suitable, marginally suitable, and unsuitable planting areas of 1 214.42 km2, 2 015.60 km2, and 3 115.03 km2, and 5 580.02km2were identified forA.rivieri, respectively. The potential suitable planting areas forA.rivieriwere mainly distributed in the south central part of Hanzhong district, the south central part of Ankang district, and the southeast part of Shanluo district. This information on the potential suitable planting area ofA.rivieriis valuable for providing baseline data, scientific information, and a research platform for understanding the ecology, geography, and environmental science of this important medicinal species.

    AHP; GIS; potential planting distribution;Amorphophallusrivieri; Qinling mountains of Shaanxi Province

    國家自然科學基金資助項目(31070293);國家“十一五”科技支撐計劃項目(2006BAI06A13-06)

    2013- 09- 12;

    日期:2015- 04- 14

    10.5846/stxb201309122261

    *通訊作者Corresponding author.E-mail: weihy@snnu.edu.cn

    高蓓,衛(wèi)海燕,郭彥龍,顧蔚.基于層次分析法和GIS的秦嶺地區(qū)魔芋潛在分布研究 .生態(tài)學報,2015,35(21):7108- 7116.

    Gao B, Wei H Y, Guo Y L, Gu W.Potential distribution ofAmorphophallusrivieriin the Qinling Mountains based on analytic hierarchy process and geographic information system.Acta Ecologica Sinica,2015,35(21):7108- 7116.

    猜你喜歡
    種植區(qū)魔芋秦嶺
    暑期秦嶺游
    天天愛科學(2022年4期)2022-05-23 12:41:46
    洞穿秦嶺
    當代陜西(2022年5期)2022-04-19 12:10:28
    C市主要草莓種植區(qū)土壤重金屬鎘、鉛現(xiàn)狀調(diào)查
    不同種植區(qū)隴東苜蓿營養(yǎng)價值的比較研究
    中國飼料(2021年17期)2021-11-02 08:15:24
    魔芋產(chǎn)業(yè)化,打開致富新大門
    草莓種植區(qū)土壤中典型有機氮化合物的分布及來源
    好忙好忙的秦嶺
    巨魔芋
    魔芋腐乳發(fā)酵工藝優(yōu)化
    中國釀造(2016年12期)2016-03-01 03:08:27
    河北昌黎縣葡萄種植區(qū)農(nóng)業(yè)地球化學特征
    亚洲aⅴ乱码一区二区在线播放| 天天躁日日操中文字幕| 天天一区二区日本电影三级| 老熟妇乱子伦视频在线观看| 国产亚洲精品av在线| 久久国产精品人妻蜜桃| 国产色婷婷99| 欧美高清成人免费视频www| 国产av麻豆久久久久久久| 十八禁国产超污无遮挡网站| 亚洲在线自拍视频| www.熟女人妻精品国产| 亚洲av电影不卡..在线观看| 99久久成人亚洲精品观看| 性欧美人与动物交配| 国产亚洲精品久久久com| 亚洲18禁久久av| 精华霜和精华液先用哪个| 亚洲专区中文字幕在线| 亚洲美女黄片视频| 搡老熟女国产l中国老女人| 精品国内亚洲2022精品成人| 亚洲国产精品合色在线| 午夜福利高清视频| 3wmmmm亚洲av在线观看| 美女免费视频网站| 日本精品一区二区三区蜜桃| 床上黄色一级片| 美女高潮喷水抽搐中文字幕| 99国产综合亚洲精品| 99久久无色码亚洲精品果冻| 国产精华一区二区三区| 久99久视频精品免费| 国产精品美女特级片免费视频播放器| 一区二区三区免费毛片| 国产成人aa在线观看| 一卡2卡三卡四卡精品乱码亚洲| 能在线免费观看的黄片| 热99re8久久精品国产| 一级黄片播放器| 精品无人区乱码1区二区| 性欧美人与动物交配| 老司机午夜福利在线观看视频| 久久九九热精品免费| 51国产日韩欧美| 午夜精品在线福利| 最近最新免费中文字幕在线| 变态另类丝袜制服| 午夜福利在线观看免费完整高清在 | 波多野结衣高清作品| 日本熟妇午夜| 久久精品国产亚洲av天美| 两人在一起打扑克的视频| 精品人妻一区二区三区麻豆 | 大型黄色视频在线免费观看| 国产免费av片在线观看野外av| 宅男免费午夜| 欧美一区二区国产精品久久精品| 色哟哟哟哟哟哟| 91字幕亚洲| 日韩精品青青久久久久久| 男女下面进入的视频免费午夜| 亚洲熟妇中文字幕五十中出| 男女下面进入的视频免费午夜| 中文字幕精品亚洲无线码一区| 18+在线观看网站| av女优亚洲男人天堂| 久久午夜亚洲精品久久| 久久久色成人| 搡老妇女老女人老熟妇| 久久性视频一级片| 国产一区二区在线av高清观看| 免费在线观看日本一区| 中文字幕av在线有码专区| 免费看光身美女| 麻豆成人av在线观看| 又粗又爽又猛毛片免费看| 三级男女做爰猛烈吃奶摸视频| 成年女人看的毛片在线观看| ponron亚洲| 此物有八面人人有两片| 色综合亚洲欧美另类图片| 我要看日韩黄色一级片| 狠狠狠狠99中文字幕| 18禁黄网站禁片午夜丰满| 99久久成人亚洲精品观看| 99精品在免费线老司机午夜| 看十八女毛片水多多多| 国产精品人妻久久久久久| 黄色配什么色好看| 少妇熟女aⅴ在线视频| 国产精品99久久久久久久久| 网址你懂的国产日韩在线| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av天美| 精品日产1卡2卡| 亚洲精华国产精华精| 成人av在线播放网站| 国产v大片淫在线免费观看| 精品久久久久久久人妻蜜臀av| 一个人看视频在线观看www免费| 伦理电影大哥的女人| 国产一级毛片七仙女欲春2| 亚洲av成人不卡在线观看播放网| 精华霜和精华液先用哪个| 他把我摸到了高潮在线观看| 五月伊人婷婷丁香| 色吧在线观看| 88av欧美| 99国产精品一区二区蜜桃av| 在线观看美女被高潮喷水网站 | 一二三四社区在线视频社区8| 欧美性感艳星| 少妇被粗大猛烈的视频| 午夜福利18| 麻豆成人av在线观看| 日韩欧美国产一区二区入口| 91av网一区二区| 午夜影院日韩av| 3wmmmm亚洲av在线观看| 国产精品,欧美在线| 亚洲av美国av| 亚洲性夜色夜夜综合| 淫秽高清视频在线观看| 日本a在线网址| 国产精华一区二区三区| x7x7x7水蜜桃| 狂野欧美白嫩少妇大欣赏| 欧美xxxx性猛交bbbb| 久久久国产成人精品二区| 亚洲av不卡在线观看| 精品日产1卡2卡| 久久伊人香网站| 看免费av毛片| 成人国产综合亚洲| 最近最新中文字幕大全电影3| 99热这里只有是精品50| 成人欧美大片| 亚洲中文字幕日韩| 午夜免费成人在线视频| 国产野战对白在线观看| 国产免费男女视频| 成人av一区二区三区在线看| 婷婷色综合大香蕉| 好男人电影高清在线观看| 中文字幕av在线有码专区| 亚洲色图av天堂| 久久国产精品影院| 伦理电影大哥的女人| 亚洲专区国产一区二区| 日韩欧美精品v在线| 我的女老师完整版在线观看| 淫秽高清视频在线观看| 国产精华一区二区三区| 内地一区二区视频在线| a级一级毛片免费在线观看| 日韩精品中文字幕看吧| 最好的美女福利视频网| 亚洲18禁久久av| 男女视频在线观看网站免费| 一进一出抽搐动态| 乱码一卡2卡4卡精品| 欧美绝顶高潮抽搐喷水| 免费在线观看影片大全网站| 麻豆久久精品国产亚洲av| 色综合亚洲欧美另类图片| 男女之事视频高清在线观看| 国产亚洲精品久久久com| 蜜桃久久精品国产亚洲av| 欧美高清成人免费视频www| 一个人看的www免费观看视频| 无遮挡黄片免费观看| 色综合站精品国产| 欧美成狂野欧美在线观看| 亚洲国产精品久久男人天堂| 在线十欧美十亚洲十日本专区| 99久久久亚洲精品蜜臀av| 宅男免费午夜| 久久亚洲真实| 国产单亲对白刺激| 在线观看舔阴道视频| 日韩欧美国产一区二区入口| 国产精品日韩av在线免费观看| 久9热在线精品视频| 亚洲不卡免费看| 日韩欧美在线乱码| 91久久精品电影网| 男女下面进入的视频免费午夜| 深爱激情五月婷婷| 熟女人妻精品中文字幕| 亚洲欧美激情综合另类| 日韩大尺度精品在线看网址| 在线观看66精品国产| 久久婷婷人人爽人人干人人爱| 一夜夜www| 在线免费观看不下载黄p国产 | 尤物成人国产欧美一区二区三区| 伊人久久精品亚洲午夜| 亚洲专区中文字幕在线| www.999成人在线观看| 欧美在线一区亚洲| 桃色一区二区三区在线观看| 久久久久久久久大av| 亚洲熟妇熟女久久| 黄片小视频在线播放| 国产精品精品国产色婷婷| 青草久久国产| 免费看日本二区| 国产精品不卡视频一区二区 | 天天躁日日操中文字幕| 国产v大片淫在线免费观看| 国产伦在线观看视频一区| 婷婷精品国产亚洲av在线| 久久精品综合一区二区三区| 国产亚洲av嫩草精品影院| 欧美另类亚洲清纯唯美| 伦理电影大哥的女人| 国产午夜精品久久久久久一区二区三区 | 日本精品一区二区三区蜜桃| 青草久久国产| 天堂网av新在线| 九九在线视频观看精品| 久久国产乱子伦精品免费另类| 亚洲性夜色夜夜综合| 亚洲国产色片| 乱码一卡2卡4卡精品| 亚洲色图av天堂| 人人妻人人澡欧美一区二区| 国产91精品成人一区二区三区| 丝袜美腿在线中文| 美女高潮的动态| 成人欧美大片| 成人特级黄色片久久久久久久| 久久这里只有精品中国| 亚洲最大成人av| 麻豆成人午夜福利视频| 国产又黄又爽又无遮挡在线| 久久久久久久久久成人| 好男人电影高清在线观看| 丁香欧美五月| 成年女人永久免费观看视频| 欧美成人一区二区免费高清观看| 夜夜夜夜夜久久久久| 成人精品一区二区免费| 日本黄大片高清| 日韩欧美一区二区三区在线观看| 欧美在线一区亚洲| 色尼玛亚洲综合影院| 少妇高潮的动态图| 村上凉子中文字幕在线| 亚洲无线在线观看| 99久久精品国产亚洲精品| 欧美+日韩+精品| 国产野战对白在线观看| 少妇裸体淫交视频免费看高清| 欧美黄色淫秽网站| 国产午夜精品论理片| 亚洲精品一区av在线观看| 午夜亚洲福利在线播放| 91麻豆av在线| 天天一区二区日本电影三级| 亚洲片人在线观看| 俺也久久电影网| 国产av麻豆久久久久久久| 三级国产精品欧美在线观看| АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 亚洲欧美日韩东京热| 亚洲国产精品999在线| 天美传媒精品一区二区| 亚洲一区高清亚洲精品| 韩国av一区二区三区四区| 久久国产乱子伦精品免费另类| 欧美激情在线99| 亚洲男人的天堂狠狠| 午夜两性在线视频| 男女床上黄色一级片免费看| 成人性生交大片免费视频hd| 欧美一区二区亚洲| 亚洲av熟女| 国内精品久久久久久久电影| 国产激情偷乱视频一区二区| h日本视频在线播放| 国产精品av视频在线免费观看| 亚洲精品影视一区二区三区av| 非洲黑人性xxxx精品又粗又长| 日韩欧美在线二视频| 少妇高潮的动态图| 老熟妇仑乱视频hdxx| 欧美性猛交╳xxx乱大交人| 午夜精品久久久久久毛片777| 蜜桃亚洲精品一区二区三区| 欧美日韩黄片免| www.999成人在线观看| 麻豆国产av国片精品| av专区在线播放| ponron亚洲| 制服丝袜大香蕉在线| 蜜桃久久精品国产亚洲av| 亚洲最大成人中文| 少妇丰满av| 欧美性感艳星| 日本撒尿小便嘘嘘汇集6| 久久久久国产精品人妻aⅴ院| 久久精品夜夜夜夜夜久久蜜豆| 欧美zozozo另类| 校园春色视频在线观看| 亚洲中文字幕日韩| 九九热线精品视视频播放| 一级黄片播放器| 国产熟女xx| 久久99热这里只有精品18| 一本久久中文字幕| 波多野结衣巨乳人妻| eeuss影院久久| 亚洲精品影视一区二区三区av| 欧美国产日韩亚洲一区| 午夜福利在线观看免费完整高清在 | 国产高清视频在线播放一区| 日韩欧美三级三区| 精品一区二区三区视频在线| 国产精品电影一区二区三区| 国产三级在线视频| 免费搜索国产男女视频| 精品久久国产蜜桃| 18禁裸乳无遮挡免费网站照片| 精品一区二区三区av网在线观看| av在线观看视频网站免费| 国产一区二区在线av高清观看| 女人被狂操c到高潮| 又紧又爽又黄一区二区| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 欧美激情久久久久久爽电影| 国产蜜桃级精品一区二区三区| 18禁黄网站禁片午夜丰满| 欧美绝顶高潮抽搐喷水| 国产91精品成人一区二区三区| 精品熟女少妇八av免费久了| 中文资源天堂在线| 亚洲片人在线观看| 国产亚洲av嫩草精品影院| 亚洲国产日韩欧美精品在线观看| 国产一区二区在线观看日韩| 毛片一级片免费看久久久久 | 久久久久性生活片| 欧美日韩国产亚洲二区| av天堂在线播放| 久久天躁狠狠躁夜夜2o2o| 国产在线男女| av中文乱码字幕在线| 国模一区二区三区四区视频| 观看免费一级毛片| 国产精品综合久久久久久久免费| 国产高清有码在线观看视频| 久久伊人香网站| 欧美一区二区精品小视频在线| 午夜亚洲福利在线播放| 真实男女啪啪啪动态图| 麻豆成人av在线观看| 噜噜噜噜噜久久久久久91| 欧美黄色淫秽网站| 在线观看免费视频日本深夜| 一进一出好大好爽视频| av女优亚洲男人天堂| 色av中文字幕| 在线十欧美十亚洲十日本专区| 日韩欧美一区二区三区在线观看| 国产精品自产拍在线观看55亚洲| 国产精品伦人一区二区| 99久久九九国产精品国产免费| 91麻豆精品激情在线观看国产| 在线观看午夜福利视频| 99热这里只有是精品50| 国产午夜福利久久久久久| 18禁在线播放成人免费| 午夜福利成人在线免费观看| 欧美激情久久久久久爽电影| 在线免费观看的www视频| 简卡轻食公司| 欧美日本视频| 欧美日本亚洲视频在线播放| 亚洲精品久久国产高清桃花| 久久久久久久久中文| 婷婷六月久久综合丁香| 熟女人妻精品中文字幕| 欧美高清成人免费视频www| 欧美在线黄色| 国产精品美女特级片免费视频播放器| 国产精品人妻久久久久久| 亚洲国产精品sss在线观看| 久久久久国内视频| 一本综合久久免费| 亚洲国产欧洲综合997久久,| 午夜激情欧美在线| 欧美最新免费一区二区三区 | 欧美黄色淫秽网站| 永久网站在线| a级毛片a级免费在线| 亚洲天堂国产精品一区在线| 亚洲成人精品中文字幕电影| 18禁裸乳无遮挡免费网站照片| 亚洲国产日韩欧美精品在线观看| 草草在线视频免费看| 狂野欧美白嫩少妇大欣赏| 757午夜福利合集在线观看| 亚州av有码| 中文字幕精品亚洲无线码一区| 亚洲av二区三区四区| 久久人人爽人人爽人人片va | 免费在线观看影片大全网站| 99久久无色码亚洲精品果冻| 成人精品一区二区免费| 日韩欧美一区二区三区在线观看| 97超级碰碰碰精品色视频在线观看| 国产一级毛片七仙女欲春2| 精品免费久久久久久久清纯| 中文字幕高清在线视频| 久久中文看片网| 色av中文字幕| 97碰自拍视频| 亚洲最大成人手机在线| 赤兔流量卡办理| 男人舔女人下体高潮全视频| 亚洲成人久久爱视频| 欧美区成人在线视频| 99国产综合亚洲精品| 啦啦啦观看免费观看视频高清| 一a级毛片在线观看| 国产中年淑女户外野战色| 日本熟妇午夜| 国产一区二区激情短视频| 亚洲欧美激情综合另类| 亚洲精品在线美女| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清在线视频| 精华霜和精华液先用哪个| 一进一出抽搐动态| 亚洲成人免费电影在线观看| 亚洲av第一区精品v没综合| 欧美潮喷喷水| 亚洲真实伦在线观看| 男人的好看免费观看在线视频| 窝窝影院91人妻| 国产免费av片在线观看野外av| 两人在一起打扑克的视频| 国产精品1区2区在线观看.| 欧美一区二区国产精品久久精品| 欧美日韩亚洲国产一区二区在线观看| 日本免费一区二区三区高清不卡| 日韩高清综合在线| 亚洲一区二区三区不卡视频| 国产午夜精品久久久久久一区二区三区 | 黄色配什么色好看| 国产午夜福利久久久久久| 一个人观看的视频www高清免费观看| 校园春色视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 免费人成视频x8x8入口观看| a级一级毛片免费在线观看| 亚洲美女视频黄频| 久久久国产成人免费| 91狼人影院| 黄色女人牲交| 麻豆av噜噜一区二区三区| 91在线观看av| 国产高清视频在线观看网站| 成年女人毛片免费观看观看9| 12—13女人毛片做爰片一| 亚洲精品久久国产高清桃花| 人妻制服诱惑在线中文字幕| 有码 亚洲区| 别揉我奶头~嗯~啊~动态视频| 看黄色毛片网站| 哪里可以看免费的av片| 精品人妻一区二区三区麻豆 | 俺也久久电影网| 免费看光身美女| 亚洲av免费在线观看| 亚洲一区高清亚洲精品| 三级毛片av免费| 看黄色毛片网站| 有码 亚洲区| 免费在线观看亚洲国产| 亚洲国产精品合色在线| 亚洲av一区综合| 99热精品在线国产| 中文字幕人成人乱码亚洲影| 亚洲国产色片| 国产黄片美女视频| 搡老妇女老女人老熟妇| 成人欧美大片| 十八禁网站免费在线| 老司机福利观看| 日本免费一区二区三区高清不卡| 一边摸一边抽搐一进一小说| 亚洲最大成人中文| 色哟哟·www| 乱码一卡2卡4卡精品| 99国产极品粉嫩在线观看| 美女cb高潮喷水在线观看| netflix在线观看网站| 深夜精品福利| 级片在线观看| 亚洲av成人不卡在线观看播放网| 亚洲黑人精品在线| 国产精品国产高清国产av| 亚洲内射少妇av| 又爽又黄a免费视频| 搡老岳熟女国产| 麻豆av噜噜一区二区三区| 中文字幕av成人在线电影| 欧美性猛交黑人性爽| 久久人妻av系列| 99久久精品热视频| 国产精品久久久久久久久免 | 男人舔女人下体高潮全视频| xxxwww97欧美| 在线天堂最新版资源| 黄色女人牲交| 18禁黄网站禁片午夜丰满| 国产日本99.免费观看| 成人av一区二区三区在线看| 嫩草影院入口| 欧美极品一区二区三区四区| 精品不卡国产一区二区三区| 18美女黄网站色大片免费观看| 黄色配什么色好看| 国产av不卡久久| 级片在线观看| 久久精品国产亚洲av香蕉五月| 国产精品伦人一区二区| 午夜福利高清视频| 欧美成狂野欧美在线观看| 激情在线观看视频在线高清| 久久久久久国产a免费观看| 高潮久久久久久久久久久不卡| 色综合亚洲欧美另类图片| 亚洲男人的天堂狠狠| 五月伊人婷婷丁香| 别揉我奶头~嗯~啊~动态视频| .国产精品久久| 美女黄网站色视频| 色综合欧美亚洲国产小说| 亚洲精品一区av在线观看| 亚洲欧美激情综合另类| 亚洲 国产 在线| 久久亚洲真实| 亚洲最大成人中文| 最近中文字幕高清免费大全6 | 国内久久婷婷六月综合欲色啪| 最好的美女福利视频网| 中文字幕精品亚洲无线码一区| 国产精品久久久久久久电影| 嫩草影院新地址| 欧美日本视频| 琪琪午夜伦伦电影理论片6080| 中文字幕av在线有码专区| 在线观看一区二区三区| 午夜福利成人在线免费观看| 最新中文字幕久久久久| 久久亚洲真实| 国产精品1区2区在线观看.| 91在线精品国自产拍蜜月| 精品久久久久久成人av| 精品久久久久久,| 亚洲国产色片| 午夜两性在线视频| 高清毛片免费观看视频网站| 精品福利观看| 精品免费久久久久久久清纯| 日韩国内少妇激情av| 久久婷婷人人爽人人干人人爱| 精品日产1卡2卡| eeuss影院久久| 51国产日韩欧美| 中文字幕av成人在线电影| 色综合站精品国产| 午夜视频国产福利| 国产成人福利小说| 高清在线国产一区| 深夜精品福利| 午夜精品一区二区三区免费看| 日韩欧美三级三区| 亚洲av成人精品一区久久| 国产精品av视频在线免费观看| 麻豆av噜噜一区二区三区| 久久欧美精品欧美久久欧美| 99在线视频只有这里精品首页| 亚洲中文字幕一区二区三区有码在线看| 欧美xxxx黑人xx丫x性爽| 久久国产精品影院| 国产精品嫩草影院av在线观看 | 国产亚洲精品av在线| 亚洲狠狠婷婷综合久久图片| 日韩国内少妇激情av| 国产伦精品一区二区三区四那| 成人国产综合亚洲| 亚洲七黄色美女视频| 国产精品人妻久久久久久| 88av欧美| 国产在线男女| 哪里可以看免费的av片| 九色国产91popny在线| 高清日韩中文字幕在线| 国产综合懂色| 高清日韩中文字幕在线| 国产成人啪精品午夜网站| 国产蜜桃级精品一区二区三区| 最近最新中文字幕大全电影3| 婷婷色综合大香蕉| 成人亚洲精品av一区二区| 美女免费视频网站| 国产精品亚洲av一区麻豆| 久久欧美精品欧美久久欧美|