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

    基于生態(tài)位和模糊數(shù)學的冬小麥適宜性評價

    2016-09-21 02:34:15李陽煦王培法王曉英羅陽歡
    生態(tài)學報 2016年14期
    關鍵詞:冬小麥生育期氣候

    王 麗,李陽煦,王培法,王曉英,羅陽歡,吳 浩

    1 南京信息工程大學地理與遙感學院,南京 210044  2 福州市勘測院數(shù)字城市研究中心,福州 350003

    ?

    基于生態(tài)位和模糊數(shù)學的冬小麥適宜性評價

    王麗1,*,李陽煦2,王培法1,王曉英1,羅陽歡1,吳浩1

    1 南京信息工程大學地理與遙感學院,南京2100442 福州市勘測院數(shù)字城市研究中心,福州350003

    為定量評價氣象、土壤等要素對作物生長的影響,利用江蘇省及周邊40個氣象站1980—2010年日氣象資料,選取影響冬小麥品種生長發(fā)育的主要氣象(如氣溫、降水、日照時數(shù)等)和土壤(如土壤厚度、有機質(zhì)和pH值等)生態(tài)因子,基于生態(tài)位理論和模糊數(shù)學的方法,依據(jù)冬小麥生長對生態(tài)各因子的響應關系建立各生態(tài)因子適宜度模型,借助GIS空間插值和空間分析,計算江蘇省冬小麥品種的種植適宜度并進行等級劃分。結果表明:從單要素看,研究區(qū)氣溫適宜度和日照時數(shù)適宜度都由北往南逐漸降低,降水適宜度由南往北逐漸降低;大部分地區(qū)土壤厚度和土壤有機質(zhì)適宜,但大部分地區(qū)的pH值適宜度較低;從氣候適宜度看,江蘇省氣候適宜度的范圍為0.68—0.81,中北部大部分區(qū)域適宜度大于0.73,全省的氣候適宜度均適合種植冬小麥,大體上呈現(xiàn)由北往南逐漸減小的變化;從土壤適宜度看,全省中部、南部與西北部大部分區(qū)域適宜度大于0.70,適合種植冬小麥。從綜合適宜度看,全省大部分地區(qū)適宜度為0.50—0.86,適宜于冬小麥生長。綜合考慮農(nóng)業(yè)氣候資源和土壤資源對江蘇省冬小麥品種種植的綜合適宜性進行評價,評價結果為充分利用江蘇省農(nóng)業(yè)生態(tài)資源、指導及科學制定冬小麥品種區(qū)域種植規(guī)劃提供科學依據(jù);按作物品種分生育期多角度的精細化研究方法和建立的各因子的適宜度模型可為今后作物區(qū)域適宜性評價提供一種新的思路,對同類研究具有一定的借鑒作用。

    土壤;地理信息系統(tǒng);作物;模糊數(shù)學;氣候;生態(tài)適宜度;冬小麥品種

    綜合考察光、溫、水等農(nóng)業(yè)氣候資源和土壤資源同作物生長發(fā)育和產(chǎn)量形成所必需條件的匹配性,是作物適宜性評價及作物適宜種植區(qū)劃的主要手段,也是農(nóng)業(yè)氣象學的重要研究內(nèi)容。小麥作為中國三大主要糧食之一,國內(nèi)外學者對冬小麥和春小麥適宜性評價做了大量研究。如Tuan 等[1],Mendas 等[2],吳洪顏等[3],千懷遂等[4],蒲金涌等[5],Wang Dacheng等[6],王琦等[7]對冬小麥做了適宜性的研究,李自珍等[8]針對春小麥生態(tài)位適宜度模型進行了研究。近年來,冬小麥適宜性研究由定性向客觀定量化方向發(fā)展,特別在指標的選擇上,由傳統(tǒng)的定性方法轉向集成多種地理因子的綜合環(huán)境指標定量數(shù)理統(tǒng)計方法。如錢存鳴等[9]從江蘇不同生態(tài)區(qū)選擇7個地點,根據(jù)它們的氣候條件、土壤類型和品質(zhì)特點,將全省分成3個品質(zhì)區(qū)域。王龍俊等[10]根據(jù)生態(tài)環(huán)境因子、土壤、小麥消費習慣等將江蘇省分為4個一級種植區(qū),12個二級種植區(qū),并提出框架性區(qū)劃方案。以上研究應用少數(shù)的氣候資料觀測站點,并不能精確反映整個江蘇的小麥氣候生態(tài)狀況。一些學者將GIS技術應用到農(nóng)作物精細化評價中,如代立芹等[11]利用模糊數(shù)據(jù)方法,建立河北省冬小麥各生育期氣候適宜度評價模型,并分析了30 a中各生育期氣溫、降水、日照和氣候適宜度的時空變化特征。馬曉群等[12]對安徽省冬小麥品種做了生態(tài)氣候適宜性分析和精細化區(qū)劃。這些研究定量到了空間面上,但主要考慮了氣候?qū)Χ←溕L的影響,以氣候適宜性評價為主。也有少數(shù)學者綜合考慮了氣候、土壤等多種資源參與評價,如李奇峰等[13]利用GIS技術構建評價單元,選取氣候、土壤、地形等關鍵的作物生態(tài)環(huán)境因子,采用空間聚類分析法構建作物生態(tài)適宜性評價體系。陸洲等[14]采用生態(tài)距離和二次閾值分級方法解決作物生態(tài)適宜性評價中分級缺乏物理意義的問題,并以北京的冬小麥為例做了驗證。上述研究為深入開展冬小麥適宜性評價提供了較好的研究思路。

    江蘇省是中國重要的冬小麥產(chǎn)區(qū)之一,而影響冬小麥生長的自然資源要素不僅包括氣候因素、還有土壤以及地形等因素,因此正確認識江蘇省冬小麥種植的自然資源狀況就需要盡可能全面地評價江蘇省自然資源各因子與冬小麥生長需求的匹配程度。目前從模糊數(shù)學角度全面考慮氣候、土壤中多個生態(tài)因子進行江蘇省冬小麥種植適宜性評價的研究較少,從冬小麥具體品種層次上的研究報道也較少。本文根據(jù)模糊數(shù)學的方法和生態(tài)適宜度理論,從冬小麥品種揚麥16號的各生育期的氣候指標出發(fā),綜合考慮農(nóng)業(yè)氣候資源和土壤資源對江蘇省冬小麥品種的綜合適宜性量化評價,實現(xiàn)區(qū)域空間上冬小麥品種的適宜性區(qū)劃,以期為充分利用江蘇省生態(tài)資源、指導及科學制定冬小麥品種區(qū)域種植規(guī)劃提供依據(jù)。

    1 研究數(shù)據(jù)及其處理

    1.1研究數(shù)據(jù)

    由國家氣候中心提供的江蘇省及周邊40個氣象站1980—2010年日觀測資料,包括平均氣溫、降水量、日照時數(shù)數(shù)據(jù);江蘇省土壤資料如土層厚度、土壤pH值和土壤有機質(zhì)含量,由南京土壤所提供;從STRM網(wǎng)站下載90 m空間分辨率的數(shù)字高程模型DEM(digital elevation model);江蘇省冬小麥揚麥16號品種的物候數(shù)據(jù)及產(chǎn)量數(shù)據(jù),由江蘇省農(nóng)科院提供。

    1.2數(shù)據(jù)處理

    表1給出了冬小麥品種揚麥16的7個生育期的劃分標準。首先,將江蘇省及周邊40個氣象站30 a的氣溫、降水、日照時數(shù)等氣象數(shù)據(jù)按照表1的生育期整理,并計算各生育期30 a氣候平均值。其次,為了保證站點氣象資料空間化的精度,對7個生育期分別采用反距離加權平均法(invers distance weight,IDW)、普通克里金(Ordinary Kriging,OK)、協(xié)同克里金(Co-Kriging,CK)以及徑向基函數(shù)(radial basis function,RBF)中的規(guī)則樣條、張力樣條和薄板樣條6種空間插值方法實現(xiàn)空間化,并對結果進行交叉檢驗,使用平均絕對誤差(mean absolute error)和均方根誤差(root mean square error)作為評價插值精度的依據(jù)。通過7個生育期6種插值方法的結果對比可得,氣溫最佳的插值方法為CK,降水最佳的插值方法為OK,日照時數(shù)插值最佳的方法為RBF的規(guī)則樣條插值。

    表1 揚麥16號小麥品種生育期

    2 研究方法

    2.1適宜性評價模型

    采用多準則綜合評價模型計算每個評價單元的綜合適宜度,作物適宜性評價模型[15]為:

    (1)

    式中,IFI為綜合適宜度;Fi為第i個評價指標的隸屬度;Wi為第i個評價指標的權重。

    [16],將江蘇省冬小麥的種植適宜性劃分為最適宜(>0.75)、適宜(≥0.5—0.75)、次適宜(≥0.25—0.5)和不適宜(<0.25)4個等級。

    2.2權重Wi設置

    作物的生態(tài)適應性指數(shù)是由光照、氣溫、降水、土壤類型、有機質(zhì)、pH值、氮磷鉀含量等因素共同決定的。對于某種生物,各因子的重要性不同。根據(jù)最小因素定律,生物的生長發(fā)育是由數(shù)量最缺的因子決定的[17]??紤]到不同因子在生態(tài)適應評價中決定強度的差異,本文采用考慮限制因子的相對權重法[18]確定各評價因子的權重。具體權重公式(2)如下:

    (2)

    式中,Wi,含義同上;Ui為評價單元一個因子的隸屬度;Um為評價單元中隸屬度最小的因子的隸屬度。該公式表明在評價單元格中尋找隸屬度最小的因子,若是小于0.2明顯的限制因子,則賦予較大的相對權重,否則由其他因子的限制性共同決定的。

    2.3模糊數(shù)學隸屬函數(shù)

    作物對某環(huán)境因子的適應有很適宜、適宜或不適宜之分,這些都為定性描述,而利用模糊數(shù)學隸屬函數(shù)的方法可以將各因子對作物的適宜程度定量化,此隸屬度取值為0—1之間,隸屬度值越大表明作物在該區(qū)域的生態(tài)適應性越高。依據(jù)生態(tài)位適宜度理論模型,作物對生態(tài)資源環(huán)境的要求主要有3類型:適宜區(qū)間型(如氣溫,pH值)、越多越好型(如土層厚度)、越低越好型(如土壤污染程度)。

    2.3.1氣溫隸屬函數(shù)

    氣溫存在適宜區(qū)間,作物光合作用有最低、最適和最高氣溫。對同一種作物來說,輻射強度的變化,會導致其同化率的變化,從而引起光合速率的氣溫曲線的變化。同樣,不同環(huán)境氣溫下光合速率反應曲線也是不同的。在正常的環(huán)境氣溫下,凈光合速率的最適氣溫不是一點,而是一個較寬的氣溫區(qū)域。作物通過光合作用,完成作物生長發(fā)育和產(chǎn)量形成的全過程。氣溫對光合作用強度的影響有兩種效應:一是氣溫增高時光化學過程加快而使總光合作用強度加強,二是氣溫增高時呼吸消耗增加。因此凈光合產(chǎn)物在初期隨氣溫增加而加強,而在超過最適氣溫以后,凈光合產(chǎn)物則隨氣溫增加而降低。以往學者建立的線形等模型并不能充分表達氣溫對作物的影響變化特點,根據(jù)Larcher和Mirza等氣溫與作物光合作用的關系曲線,氣溫與光合速率反應曲線近似呈S型[19-20],結合研究地區(qū)實際情況,據(jù)此建立如式(3)的氣溫隸屬函數(shù)ST:

    (3)

    式中,T是某一生育期內(nèi)的實際日平均氣溫(℃);TL、TH、T2、T1分別為作物在該時段內(nèi)生長發(fā)育的下限氣溫、上限氣溫和最適氣溫的上下限(℃)。此處參數(shù)依據(jù)參考資料[21-22]與江蘇省農(nóng)業(yè)科學院提供品種小麥栽培種植資料獲取。

    2.3.2日照隸屬函數(shù)

    作物利用太陽輻射及其自身具有的光合作用功能,將太陽能轉化為有機干物質(zhì)。本文根據(jù)Camilo L.Medina[23]、黃璜[24]的研究,把日照時數(shù)達到可照時數(shù)的70%,即日照百分率70%以上的光照條件作為冬小麥的適宜狀態(tài)。日照時數(shù)的隸屬函數(shù)SSH為[25]:

    (4)

    式中,SH為實際日照時數(shù)(h);SH0為日照百分率為 70%的日照時數(shù)(h);b為系數(shù)。

    2.3.3降水隸屬函數(shù)

    評價降水對作物生長的影響,按照作物產(chǎn)量與降水量關系曲線圖[26],建立作物的降水隸屬函數(shù)SR,計算公式為:

    (5)

    式中,R為某生育期內(nèi)的平均降水量(mm);RL、RH、R2、R1分別為作物在該時段內(nèi)生長發(fā)育的下限、上限和最適需水量的上下限(mm)。k=1,n=2。此處參數(shù)依據(jù)參考資料[21-22]與江蘇省農(nóng)業(yè)科學院提供品種小麥栽培種植資料獲取。

    2.3.4土壤pH值隸屬函數(shù)

    建立土壤pH值的隸屬度函數(shù)SpH,如式(6):

    (6)

    式中,pH為實際土壤pH值;pHSL和 pHSH分別為適宜于作物生長的下限值和上限值;pHOL和 pHOH分別為作物最佳生長的下限值和上限值,此處參數(shù)依據(jù)參考資料[21-22],分別取6.8和7.0。

    2.3.5土壤有機質(zhì)隸屬函數(shù)

    土壤有機質(zhì)隸屬度函數(shù)SOM[27]:

    (7)

    式中,OM為實際土壤有機質(zhì)含量(%);OMSL和OMOL分別為作物生長的土壤有機質(zhì)含量適宜下限值及最佳下限值(%),最佳下限值取1.3%。由于土壤有機質(zhì)因子在最不利情況下,作物仍有可能生存,因此這里將生態(tài)位適宜度模型中的最低適應值由 0 換為a,a∈ [ 0,1)。

    2.3.6土層厚度隸屬函數(shù)

    土層厚度隸屬度函數(shù)SSD[28]:

    (8)

    式中,SD為實際土層厚度(cm);SDOL作物生長的土層厚度下限,取20cm。

    3 結果與分析

    3.1各要素生態(tài)適宜度評價

    根據(jù)以上氣溫、降水量和日照時數(shù)等各隸屬函數(shù)及生態(tài)指標,可計算得到不同生育期各因子的適宜度。由氣候因子各生育期的適宜度與相應生育期的權重乘積累和計算得到氣候因子整個生育期的適宜度。冬小麥各生育期內(nèi)與生育期間氣象要素的權重結合文獻[11]和江蘇省農(nóng)科院品種小麥專家意見,得出整個生育期上各生育期的相對權重(表2),通過空間分析的柵格計算得到整個生育期的各氣象因子適宜度(圖1)。因為土壤各因子隨生育期的變化很小,可忽略,即只算整個生育期的土壤各因子的適宜度。

    表2 冬小麥各因子各生育期權重

    氣溫適宜度的范圍為0.54—0.81(圖1),由北往南逐漸降低,但在最南端如昆山、蘇州等地適宜度較高;同一緯度上,東部地區(qū)的適宜度略高于西部地區(qū);最適宜區(qū)域在江蘇省北部,適宜度較低的區(qū)域在江蘇省南部。由于播種期昆山、吳江地區(qū)適宜度為1,與周邊適宜度相差較大;拔節(jié)期南京、溧水、溧陽、宜興和吳江一帶氣溫適宜度為1,與周邊相差較大,使得將各生育期的氣溫取平均值

    后,南京與吳江地區(qū)顯得異常。

    圖1 江蘇省各因子的適宜度分布圖Fig.1 Distribution of suitability of ecological factors

    降水適宜度的范圍為0.25—0.63(圖1),基本上呈現(xiàn)由南部向北部逐漸降低;高郵、東臺區(qū)域略低于其南北兩側;最適宜區(qū)域在江蘇省南部,如無錫、蘇州等地;適宜度較低的區(qū)域在江蘇省北部,如徐州、連云港。

    日照時數(shù)適宜度的范圍為0.53—0.91(圖1),江蘇省大部分地區(qū)的日照時數(shù)都適宜冬小麥的生長。由北往南逐漸降低;同一緯度上,大部分西部地區(qū)的適宜度略低于東部地區(qū);最適宜區(qū)域在江蘇省北部,如連云港等;適宜度較低的區(qū)域在江蘇省南部,如蘇州等地。

    經(jīng)資料查閱,冬小麥的耕作層深度一般為20cm以上,土壤厚度的最低下限為40cm,計算的土壤厚度適宜度范圍較大,見圖1。除去河流、湖泊,全省土壤厚度的適宜度大部分為1,即厚度大于40cm;江蘇省的土壤厚度均適宜冬小麥的生長。

    江蘇省土壤pH值的適宜度的范圍見圖1,適宜度大于0.6的區(qū)域占全省的61%。北部、蘇中的西部和東部、蘇南部分地區(qū)pH值適宜度較高,如徐州、宿遷、連云港、蘇州等地;中部與蘇南東部的部分地區(qū)適宜度為0,土壤酸性超出冬小麥適宜范圍。全省有大約50%的面積pH值適宜度大于0.8。

    土壤有機質(zhì)含量的值達到1.3%為適宜[29],根據(jù)有機質(zhì)隸屬函數(shù),計算的江蘇省有機質(zhì)適宜度的范圍見圖1,適宜度大于0.8的區(qū)域占全省64%,大于0.6的區(qū)域占全省75%。大部分地區(qū)有機質(zhì)含量較高,中、東南部如蘇州、南通、泰州等地有機質(zhì)適宜度最高;北部地區(qū)適宜度低于南部,但基本大于0.2,屬于適宜范圍。所以,除水域外,江蘇省大部分地區(qū)的土壤有機質(zhì)適宜冬小麥生長。

    3.2氣候生態(tài)適宜度評價和土壤生態(tài)適宜度評價

    3.2.1冬小麥氣候生態(tài)適宜度

    單從氣候的角度考查江蘇省冬小麥品種的適宜性,可用氣候適宜度來定量評價。利用式(1),首先由各生育期的氣溫、降水、日照時數(shù)適宜度和冬小麥生育期內(nèi)權重計算得到各生育期氣候適宜度,再將不同生育期氣候適宜度和相應的生育期間權重乘積求和,即得到冬小麥整個生育期的氣候適宜度范圍為0.68—0.81(圖2),江蘇省中北部大部分區(qū)域適宜度大于0.73,全省的氣候適宜度均適合種植冬小麥,大體上呈現(xiàn)由北往南逐漸減小的變化。鎮(zhèn)江、南通一帶適宜度略低于其南部與北部,但都大于0.60。

    3.2.2冬小麥土壤生態(tài)適宜度

    單從土壤的角度考查江蘇省冬小麥品種的適宜性,可用土壤適宜度來定量評價。利用式(1),綜合土層厚度、土壤pH值、土壤有機質(zhì)單因子的適宜度,其中3因子權重賦值依據(jù)文獻[30],分別為0.29、0.19和0.52計算得到江蘇省冬小麥土壤適宜度分布圖3。由圖3看出江蘇省土壤適宜度為0—1,中部、南部與西北部等大部分區(qū)域土壤適宜度大于0.70,適合種植冬小麥,其中0.7以上區(qū)域占全省的82%,適宜度大于0.89的區(qū)域占全省的35.7%,主要分布在長江、淮河沿岸、太倉等區(qū)域??臻g分布圖上呈現(xiàn)南部優(yōu)于北部的趨勢。

    圖2  江蘇省氣候適宜度Fig.2 Climatic suitability of Jiangsu Province

    圖3 江蘇省土壤適宜度Fig.3 Soil suitability of Jiangsu Province

    3.3冬小麥綜合適宜度評價

    圖4 江蘇省冬小麥綜合適宜度分布圖 Fig.4 Comprehensive suitability for winter wheat in Jiangsu Province

    冬小麥的生長是氣候、土壤等多個環(huán)境因素的綜合作用,因此全面考慮氣候、土壤等因素的影響更為合理。根據(jù)以上氣溫、降水、光照、土壤類型、土壤有機質(zhì)、土壤pH值各因子適宜度,通過權重公式(2)可得各因子的權重圖層,利用式(1)計算得到冬小麥綜合適宜度空間分布圖4。

    從圖4中顏色深淺變化表示種植冬小麥適宜程度的高低??梢钥闯鲈诟鳉夂蛞蜃泳C合作用下江蘇省冬小麥種植適宜度呈現(xiàn)北高南低。依據(jù)文獻[16]將適宜性分為最適宜、適宜、次適宜和不適宜4個等級,同時參考研究區(qū)土地利用圖和冬小麥的產(chǎn)量等調(diào)查資料,做出分級調(diào)整,由冬小麥產(chǎn)量數(shù)據(jù)將最適宜區(qū)為適宜度0.7以上,適宜度0.65為適宜區(qū)和次適宜區(qū)的分界點,江河湖泊等水域都低于0.45,將適宜度低于0.45的劃為不適宜區(qū)。由圖4可看出江蘇省大部分區(qū)域適宜度為0.50—0.86,適宜種植冬小麥。最適宜的區(qū)域主要在江蘇北部和中部區(qū)域,南部區(qū)域主要為適宜區(qū),不適宜的區(qū)域主要分布在蘇南的太湖、滆湖,蘇中的高郵湖以及蘇北的洪澤湖、長江等水域。該分級結果和吳洪顏等[3]、金善寶[17]學者研究結論相一致。

    4 結論與討論

    本文利用多種插值方法對比選取最佳方法實現(xiàn)氣象要素的空間化,利用模糊數(shù)學建立適宜度模型,分析江蘇省冬小麥的氣候、土壤適宜度情況,并將各要素的適宜度賦予權重,進行綜合分析,得到江蘇省冬小麥適宜種植區(qū)劃。結果表明:

    1)江蘇省氣溫適宜度為0.54—0.81,日照適宜度為0.53—0.91,降水適宜度為0.25—0.63,總體上氣候條件較適宜冬小麥生長。氣溫適宜度總體上呈從北往南逐漸降低的趨勢;降水適宜度總體上呈現(xiàn)從南往北逐漸降低的趨勢;日照時數(shù)適宜度從北往南逐漸降低。

    2)江蘇省大部分區(qū)域的土壤厚度能滿足冬小麥的生長,適宜度為1;土壤有機質(zhì)基本滿足冬小麥生長,適宜度大于0.8的區(qū)域占全省64%,大于0.6的區(qū)域占全省75%;土壤pH值適宜度較適宜冬小麥生長,適宜度大于0.6的區(qū)域占全省的61%。

    3)江蘇省冬小麥整個生育期的氣候適宜度范圍為0.68—0.81,中北部大部分區(qū)域適宜度大于0.73,全省的氣候適宜度均適合種植冬小麥,大體上呈現(xiàn)由北往南逐漸減小的變化。江蘇省土壤適宜度為0—1,中部、南部與西北部等大部分區(qū)域土壤適宜度大于0.7,適合種植冬小麥,其中0.7以上區(qū)域占全省的82%,適宜度大于0.89的區(qū)域占全省的35.7%,主要分布在長江、淮河沿岸、太倉等區(qū)域??臻g分布圖上呈現(xiàn)南部優(yōu)于北部的趨勢。

    4)江蘇省冬小麥綜合適宜度看,整個江蘇省基本都適宜種植冬小麥,空間上看北部與中部最適宜冬小麥生長,南部較適宜冬小麥生長,不適宜區(qū)域主要分布在江河湖泊等水體區(qū)域。

    本文研究中以冬小麥優(yōu)質(zhì)品種揚麥16號為例,該品種屬于春性中熟小麥品種,在各生態(tài)要素的隸屬度計算時,由于缺少該品種準確的三基點資料,根據(jù)專家意見,各氣候要素的三基點參考中熟小麥的三基點資料,土壤pH值、土層厚度以及土壤有機質(zhì)含量等三基點的取值參考冬小麥大類的三基點資料。這可能與該品種的實際三基點取值有偏差。

    致謝:江蘇省農(nóng)業(yè)科學院提供冬小麥資料,范曉梅、孫德勇老師幫助寫作,特此致謝。

    參考文獻(References):

    [1]Nguyen Thanh Tuan, Qiu J J, Ann Verdoodt, Li H, Eric Van Ranst. Temperature and precipitation suitability evaluation for the winter wheat and summer maize cropping system in the Huang-Huai-Hai plain of China. Agricultural Sciences in China, 2011, 10(2): 275-288.

    [2]Abdelkader Mendasa, Amina Delali. Integration of multicriteria decision analysis in GIS to develop land suitability for agriculture: Application to durum wheat cultivation in the region of Mleta in Algeria. Computers and Electronics in Agriculture, 2012, 83: 117-126.

    [3]吳洪顏,商兆堂,程婷,黃文杰,陸志剛.基于聚類分析的江蘇冬小麥農(nóng)業(yè)氣候區(qū)劃研究.中國農(nóng)學通報,2012,28(26):119-124.

    [4]千懷遂,焦士興,趙峰.河南省冬小麥氣候適宜性變化研究.生態(tài)學雜志,2005,24(5):503-507.

    [5]蒲金涌,張存杰.甘肅省冬小麥水分適應性動態(tài)變化研究.資源科學,2008,30(9):1297-1402.

    [6]Wang D C, Li C J, Song X Y, Wang J H, Yang X D, Huang W J, Wang J Y, Zhou J H. Assessment of land suitability potentials for selecting winter wheat cultivation areas in Beijing, China, using RS and GIS. Agricultural Sciences in China, 2011, 10(9): 1419-1430.

    [7]王琦,郭建茂,鄭騰飛,施俊怡.基于作物模型的冬小麥氣候適宜度算法研究.中國農(nóng)學通報,2013,29(32):374-378.

    [8]李自珍,李文龍,馬智慧,楊雨華.作物生態(tài)位適宜度模型及其在春小麥生長系統(tǒng)中的應用.蘭州大學學報(自然科學版),2010,46(2):45-50.

    [9]錢存鳴,馬兆私,周朝飛,葛永福,盛培英,陳志德,姚國才,柏貴華.江蘇省小麥品質(zhì)區(qū)劃研究.江蘇農(nóng)業(yè)科學,1990(6):8-11.

    [10]王龍俊,陳榮振,朱新開,楊力,姜東,陳維新,李旭.江蘇省小麥區(qū)劃研究初報.江蘇農(nóng)業(yè)科學,2002(2):15-18.

    [11]代立芹,李春強,魏瑞江,于長文.河北省冬小麥氣候適宜度及其時空變化特征分析.中國農(nóng)業(yè)氣象,2011,32(3):399-406.

    [12]馬曉群,張宏群,吳文玉,戚尚恩,陳曉藝.安徽省冬小麥品種生態(tài)氣候適宜性分析和精細化區(qū)劃.中國農(nóng)業(yè)氣象,2012,33(1):86-92.

    [13]李奇峰,劉曦,孔箐鋅,陸洲,潘瑜春,秦向陽.基于空間聚類的作物生態(tài)適宜性評價方法研究——以北京市為例.中國農(nóng)學通報,2011,27(11):262-265.

    [14]陸洲,秦向陽,李奇峰,于瑩,臧辰龍,淮賀舉.作物生態(tài)適宜性定量化評價方法及通用工具.農(nóng)業(yè)工程學報,2012, 28(20):195-202.

    [15]Sudabe Jafari, Narges Zaredar. Land suitability analysis using multi attribute decision making approach. International Journal of Environmental Science and Development, 2010, 1(5): 441-445.

    [16]Madrau S, Zucca C, Urgeghe A M, Julitta F, Previtali F. Land suitability for crop options evaluation in areas affected by desertification: the case study of Feriana in Tunisia. Land Degradation and Desertification: Assessment, Mitigation and Remediation, 2010: 179-193.

    [17]金善寶.中國小麥學.北京:農(nóng)業(yè)出版社,1996.

    [18]曹衛(wèi)星,朱艷.作物管理知識模型.北京:中國農(nóng)業(yè)出版社,2004.

    [19]張靜.作物——地域多種組合中作物生態(tài)適宜性評價與權重配置方法的研究[D].南京:南京農(nóng)業(yè)大學,2005.

    [20]Mirza Hasanuzzaman, Kamrun Nahar, Masayuki Fujita. Extreme temperature responses, oxidative stress and antioxidant defense in plants. 2013 [2013-3-13]. http://dx.doi.org/10.5772/54833.

    [21]胡立勇,丁艷鋒.作物栽培學.北京:高等教育出版社,2008.

    [22]楊文鈺,屠乃美.作物栽培學各論.南方本(第二版).北京:中國農(nóng)業(yè)出版社,2011.

    [23]Camilo L Medina, Rogéria P Souza, Eduardo C Machado, Rafael V Ribeiro, José A B Silva. Photosynthetic response of citrus grown under reflective aluminized polypropylene shading nets. Scientia Horticulturae, 2002, 96(1/4): 115-125.

    [24]黃璜. 中國紅黃壤地區(qū)作物生產(chǎn)的氣候生態(tài)適應性研究. 自然資源學報, 1996, 11(4): 340-346.

    [25]賴純佳, 千懷遂, 段海來, 宋秋洪, 俞芬, 張靜芬, 張軼秀. 基于數(shù)據(jù)處理及圖件的小麥-水稻種植制度的氣候風險評估. 農(nóng)業(yè)工程學報, 2011, 27(2): 230-236.

    [26]徐學選, 高鵬, 蔣定生. 延安降水對農(nóng)作物生長適宜性的模糊分析. 水土保持研究, 2000, 7(2): 73-76, 118-118.

    [27]歐陽志云, 王如松, 符貴南. 生態(tài)位適宜度模型及其在土地利用適宜性評價中的應用. 生態(tài)學報, 1996, 16(2): 113-120.

    [28]單英杰. 基于GIS和模型的種植設計系統(tǒng)研究[D]. 南京: 南京農(nóng)業(yè)大學, 2007.

    [29]山東省小麥高產(chǎn)栽培研究協(xié)作組. 山東省冬小麥穩(wěn)產(chǎn)千斤的生育特點及栽培規(guī)范. 山東農(nóng)業(yè)科學, 1981, (2): 1-7.

    [30]于婧, 聶艷, 周勇, 何佑勇. 生態(tài)位適宜度方法在基于GIS的耕地多宜性評價中的應用. 土壤學報, 2006, 43(2): 190-196.

    Assessment of ecological suitability of winter wheat in Jiangsu Province based on the niche-fitness theory and fuzzy mathematics

    WANG Li1,*, LI Yangxu2, WANG Peifa1, WANG Xiaoying1, LUO Yanghuan1, WU Hao1

    1SchoolofRemoteSensing,NanjingUniversityofInformationScienceandTechnology,Nanjing210044,China2TheResearchCenterofDigitalCity,FuzhouInvestigationandSurveyingInstitute,Fuzhou350003,China

    In order to quantitatively evaluate the influence of weather and soil on the growth of wheat crops, we use the day-to-day meteorological data from Jiangsu Province during 19802010, and develop an ecological suitability model of winter wheat on the basis of the niche suitability theory and fuzzy mathematics. The growth suitability of winter wheat in Jiangsu Province is evaluated and the regions in the province are delineated for different degrees of suitability with the aid of GIS spatial interpolation and analysis technologies; the different degrees of suitability assigned included the most suitable, more suitable, suitable, and not suitable. The major ecological factors affecting winter wheat growth such as temperature, precipitation, sunshine, soil depth, organic matter, and pH value were considered during the analysis. The analysis process includes four parts. First, the analysis was performed taking into account every growth period of winter wheat, to gain spatial distribution of meteorological stations during 19802010, using the most suitable method from the six spatial interpolation methods, including IDW, Kriging, and spline. Second, the analysis was conducted using the fuzzy mathematical models and considering the obtained ecoclimatic indices for the crop′s growth period, and the suitability of temperature, precipitation, and sunlight in different growth periods. Third, the weight of each meteorological factor at the different stages, the suitability of temperature, precipitation, and sunlight for the entire life period were calculated for analysis. Lastly, the comprehensive adaptability model for winter wheat was constructed using the relative weight method which determine the weight of all factors, each high suitability degree factor multiplied to reduce its weight while each weak factor multiplied to increase its weight. The integrated suitability values of winter wheat were calculated using this comprehensive adaptability model taking into consideration the previous results that took into account the temperature, precipitation, sunshine, soil thickness, soil pH, and soil organic matter concentration. The conclusions of our research are: 1) The suitability of temperature and sunshine hours gradually decreases from north to south within the study area, while the suitability of rainfall gradually decreases from south to north. 2) The suitability of soil thickness and soil organic matter are higher for crop growth in most of the study area, while the suitability of the soil pH value is low. For climatic suitability, the suitability value of Jiangsu Province varied from 0.68 to 0.81, where the suitability in middle and north area exceeds 0.73. 3) The climate of Jiangsu Province is suitable for the growth of winter wheat, and the suitability value decreases from north to south. Considering soil suitability, the value for middle, south and northwest parts of study area are more than 0.7, which indicates that they are suitable for the growth of winter wheat. 4) The comprehensive suitability degree for most part of research area varies from 0.5 to 0.86, which is considered suitable for the growth of winter wheat. Our research can provide reasonable guidance to the scientific plans for the cultivation of winter wheat in Jiangsu Province. Furthermore, the developed suitability model that considers crop variety and its growth stages provides a new method to the crop ecological suitability evaluation, and can be used in similar researches.

    soil; geographic information systems; crops; fuzzy mathematics; climate; ecological suitability; winter wheat varieties

    國家自然科學基金項目(41175077);黔科合重大專項([2011]6003);2014年度大學生實踐創(chuàng)新訓練計劃江蘇省級重點項目(201410300017Z)

    2014-12-02; 網(wǎng)絡出版日期:2015-10-30

    Corresponding author.E-mail: xuer217@163.com

    10.5846/stxb201412022394

    王麗,李陽煦,王培法,王曉英,羅陽歡,吳浩.基于生態(tài)位和模糊數(shù)學的冬小麥適宜性評價.生態(tài)學報,2016,36(14):4465-4474.

    Wang L, Li Y X, Wang P F, Wang X Y, Luo Y H, Wu H.Assessment of ecological suitability of winter wheat in Jiangsu Province based on the niche-fitness theory and fuzzy mathematics.Acta Ecologica Sinica,2016,36(14):4465-4474.

    猜你喜歡
    冬小麥生育期氣候
    瞧,氣候大不同
    氣候變暖會怎樣?
    甘肅冬小麥田
    植物保護(2017年1期)2017-02-13 06:44:34
    基于作物生育期的潛在蒸散的時空演變特征及R/S 分析
    西藏科技(2015年5期)2015-09-26 11:55:25
    冬小麥和春小麥
    中學生(2015年4期)2015-08-31 02:53:50
    立冬
    2013-2014年度二二二團冬小麥各生育期氣象條件分析
    阿拉爾市2012年棉花生育期氣象條件分析
    冬小麥——新冬18號
    冬小麥—新冬41號
    国产精品一区二区在线观看99| 亚洲av日韩在线播放| 亚洲综合精品二区| 在线免费十八禁| 日韩免费高清中文字幕av| 免费观看无遮挡的男女| 观看美女的网站| 日本av手机在线免费观看| 日日摸夜夜添夜夜爱| 91久久精品国产一区二区成人| 国产av一区二区精品久久 | 一区二区三区四区激情视频| 1000部很黄的大片| 黄色配什么色好看| 日日摸夜夜添夜夜爱| 又大又黄又爽视频免费| 欧美日韩一区二区视频在线观看视频在线| 免费看不卡的av| 少妇熟女欧美另类| 精品亚洲成国产av| 高清视频免费观看一区二区| 国产成人精品久久久久久| av网站免费在线观看视频| 国产黄频视频在线观看| 欧美日韩精品成人综合77777| 亚洲怡红院男人天堂| 夫妻午夜视频| av国产久精品久网站免费入址| 99久久精品国产国产毛片| 一区二区三区四区激情视频| 18禁裸乳无遮挡免费网站照片| 国产精品一区www在线观看| 色视频www国产| 亚洲精品成人av观看孕妇| 亚洲av不卡在线观看| 久久久久久久亚洲中文字幕| 成人高潮视频无遮挡免费网站| h日本视频在线播放| 纵有疾风起免费观看全集完整版| 精品久久久久久久久亚洲| 精品一区二区免费观看| 亚洲综合色惰| 最近中文字幕高清免费大全6| 欧美日韩国产mv在线观看视频 | 干丝袜人妻中文字幕| 妹子高潮喷水视频| 国产精品熟女久久久久浪| 一区二区av电影网| 午夜福利高清视频| 网址你懂的国产日韩在线| 欧美高清成人免费视频www| 亚洲三级黄色毛片| 国产精品麻豆人妻色哟哟久久| 亚洲人成网站高清观看| 久久精品国产a三级三级三级| 国产日韩欧美在线精品| 国产欧美亚洲国产| 亚洲av成人精品一二三区| 一级a做视频免费观看| 天美传媒精品一区二区| 亚洲精品日本国产第一区| 如何舔出高潮| 人妻 亚洲 视频| 亚洲中文av在线| 日本黄大片高清| 亚洲精品,欧美精品| 国产伦精品一区二区三区视频9| 啦啦啦在线观看免费高清www| 欧美少妇被猛烈插入视频| 国产亚洲午夜精品一区二区久久| 久久久精品94久久精品| 成人无遮挡网站| 色视频www国产| 亚洲国产毛片av蜜桃av| 日韩,欧美,国产一区二区三区| 国产成人精品久久久久久| 日本午夜av视频| 久久久久视频综合| 在线观看国产h片| 高清欧美精品videossex| 免费看日本二区| 日韩亚洲欧美综合| 国产乱来视频区| 嫩草影院入口| 中国国产av一级| 精品一区在线观看国产| 亚洲熟女精品中文字幕| 日本av免费视频播放| 成年美女黄网站色视频大全免费 | 天美传媒精品一区二区| 免费在线观看成人毛片| 日日啪夜夜爽| 久久久久国产网址| 精品人妻视频免费看| 国产精品三级大全| 亚洲精品久久午夜乱码| 少妇人妻精品综合一区二区| 国产日韩欧美在线精品| 久久国产精品大桥未久av | 天堂8中文在线网| 久久久午夜欧美精品| 国产成人一区二区在线| 最近最新中文字幕大全电影3| av在线观看视频网站免费| 久久女婷五月综合色啪小说| 不卡视频在线观看欧美| 亚洲丝袜综合中文字幕| 国产淫片久久久久久久久| 亚洲精品国产av蜜桃| 久热这里只有精品99| 久久女婷五月综合色啪小说| 精品人妻视频免费看| 少妇的逼好多水| 男人舔奶头视频| 亚洲自偷自拍三级| 天美传媒精品一区二区| 久久久久网色| 亚洲精品一区蜜桃| 一个人看的www免费观看视频| 日产精品乱码卡一卡2卡三| xxx大片免费视频| 亚洲aⅴ乱码一区二区在线播放| 欧美xxxx性猛交bbbb| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美成人精品一区二区| 亚洲精品aⅴ在线观看| 22中文网久久字幕| 日本wwww免费看| 久久ye,这里只有精品| 中文字幕精品免费在线观看视频 | 亚洲美女视频黄频| 国产男女内射视频| 亚洲国产最新在线播放| 久久久久久人妻| 97超碰精品成人国产| 国产在线视频一区二区| 街头女战士在线观看网站| 亚洲精品色激情综合| 啦啦啦视频在线资源免费观看| 人妻少妇偷人精品九色| 亚洲丝袜综合中文字幕| 少妇被粗大猛烈的视频| 男女无遮挡免费网站观看| 水蜜桃什么品种好| 黄色怎么调成土黄色| 欧美bdsm另类| 精品熟女少妇av免费看| 日本免费在线观看一区| 亚洲国产日韩一区二区| 少妇的逼水好多| 国产高清三级在线| tube8黄色片| 久久亚洲国产成人精品v| 天天躁夜夜躁狠狠久久av| 亚洲国产欧美人成| 欧美成人精品欧美一级黄| 亚洲国产日韩一区二区| 国产国拍精品亚洲av在线观看| 简卡轻食公司| 久久久久久久久久成人| 街头女战士在线观看网站| 欧美一区二区亚洲| 最近最新中文字幕大全电影3| 身体一侧抽搐| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品电影小说 | 免费看日本二区| 午夜精品国产一区二区电影| 建设人人有责人人尽责人人享有的 | 色5月婷婷丁香| av国产免费在线观看| 亚洲成人手机| 赤兔流量卡办理| 五月开心婷婷网| av福利片在线观看| 日韩一区二区视频免费看| 亚洲精品一区蜜桃| 蜜桃久久精品国产亚洲av| 极品教师在线视频| 中国美白少妇内射xxxbb| 美女主播在线视频| 亚洲欧美中文字幕日韩二区| 夜夜爽夜夜爽视频| 国产永久视频网站| 欧美成人a在线观看| 国产日韩欧美在线精品| 赤兔流量卡办理| 小蜜桃在线观看免费完整版高清| 亚洲图色成人| av黄色大香蕉| 在线观看一区二区三区激情| 国产精品不卡视频一区二区| 国产伦理片在线播放av一区| 国产欧美日韩精品一区二区| 老熟女久久久| 两个人的视频大全免费| 亚洲色图综合在线观看| 成人一区二区视频在线观看| 免费播放大片免费观看视频在线观看| 亚洲欧美一区二区三区国产| 久久久久久久久久久免费av| 成人特级av手机在线观看| 久久久久久久大尺度免费视频| 在线天堂最新版资源| 一级毛片电影观看| 精品亚洲成a人片在线观看 | 免费黄网站久久成人精品| 欧美精品国产亚洲| 免费看日本二区| 亚洲欧美一区二区三区黑人 | 青春草国产在线视频| 免费黄频网站在线观看国产| 国产毛片在线视频| 欧美精品人与动牲交sv欧美| 综合色丁香网| 大又大粗又爽又黄少妇毛片口| 国产午夜精品久久久久久一区二区三区| 天天躁日日操中文字幕| 97超碰精品成人国产| 在线观看一区二区三区| 久久 成人 亚洲| 一级黄片播放器| 国产色爽女视频免费观看| 国产黄片美女视频| 亚洲av电影在线观看一区二区三区| 99久久人妻综合| 日韩人妻高清精品专区| 一区二区三区乱码不卡18| 男人添女人高潮全过程视频| 欧美bdsm另类| 人妻系列 视频| 欧美日韩综合久久久久久| 久久影院123| 国产日韩欧美亚洲二区| av播播在线观看一区| 国产免费一区二区三区四区乱码| 尾随美女入室| 少妇的逼好多水| 一级av片app| 欧美日韩一区二区视频在线观看视频在线| av国产精品久久久久影院| 亚洲丝袜综合中文字幕| 久久久久视频综合| 高清毛片免费看| 日韩一本色道免费dvd| 国产成人a∨麻豆精品| 国产成人午夜福利电影在线观看| 国产亚洲91精品色在线| 直男gayav资源| 男人舔奶头视频| 97热精品久久久久久| 欧美老熟妇乱子伦牲交| 夫妻午夜视频| 欧美97在线视频| 少妇被粗大猛烈的视频| 超碰97精品在线观看| 亚洲精品aⅴ在线观看| 免费观看无遮挡的男女| 精品少妇黑人巨大在线播放| 成人特级av手机在线观看| 精品亚洲成国产av| 久久久久精品久久久久真实原创| 国产亚洲一区二区精品| 国内少妇人妻偷人精品xxx网站| 好男人视频免费观看在线| 人人妻人人看人人澡| 久久久久久伊人网av| 日本黄色片子视频| 汤姆久久久久久久影院中文字幕| 亚洲国产色片| 国产白丝娇喘喷水9色精品| 97在线视频观看| 97在线人人人人妻| 26uuu在线亚洲综合色| 亚洲国产精品成人久久小说| 夜夜看夜夜爽夜夜摸| 视频中文字幕在线观看| 国产伦精品一区二区三区四那| 日日撸夜夜添| 欧美高清性xxxxhd video| 22中文网久久字幕| 亚洲va在线va天堂va国产| 国产毛片在线视频| 国产亚洲最大av| 九九爱精品视频在线观看| 午夜福利在线在线| 欧美xxxx性猛交bbbb| .国产精品久久| 成人漫画全彩无遮挡| 又黄又爽又刺激的免费视频.| 精华霜和精华液先用哪个| 在线 av 中文字幕| 肉色欧美久久久久久久蜜桃| 午夜福利视频精品| 人妻制服诱惑在线中文字幕| 丰满迷人的少妇在线观看| 久久精品国产鲁丝片午夜精品| 免费人妻精品一区二区三区视频| 亚洲人与动物交配视频| 永久免费av网站大全| 美女福利国产在线 | 3wmmmm亚洲av在线观看| 久久精品夜色国产| 插逼视频在线观看| av卡一久久| 人人妻人人爽人人添夜夜欢视频 | 观看美女的网站| 久热这里只有精品99| 亚洲国产精品专区欧美| 欧美成人午夜免费资源| 一区在线观看完整版| 久久久久性生活片| 国产精品一区二区性色av| videossex国产| 精品酒店卫生间| 国产亚洲午夜精品一区二区久久| 伊人久久国产一区二区| av国产免费在线观看| 建设人人有责人人尽责人人享有的 | 少妇人妻久久综合中文| 久热久热在线精品观看| 这个男人来自地球电影免费观看 | 欧美xxxx黑人xx丫x性爽| 在线亚洲精品国产二区图片欧美 | 晚上一个人看的免费电影| 日韩欧美 国产精品| 最后的刺客免费高清国语| 国产探花极品一区二区| 91久久精品国产一区二区三区| 极品少妇高潮喷水抽搐| 日韩精品有码人妻一区| 国产深夜福利视频在线观看| 一区在线观看完整版| 日韩一区二区三区影片| 这个男人来自地球电影免费观看 | 91精品国产国语对白视频| 小蜜桃在线观看免费完整版高清| av又黄又爽大尺度在线免费看| 直男gayav资源| 黄片wwwwww| 少妇 在线观看| 国产女主播在线喷水免费视频网站| 国产免费视频播放在线视频| 中国美白少妇内射xxxbb| 少妇的逼好多水| 成人综合一区亚洲| 国产精品成人在线| 成人影院久久| 亚洲国产精品成人久久小说| 国产成人精品福利久久| 男女国产视频网站| 久热这里只有精品99| 日本猛色少妇xxxxx猛交久久| 精品久久久久久久久亚洲| 水蜜桃什么品种好| 日韩av在线免费看完整版不卡| 毛片女人毛片| av.在线天堂| 能在线免费看毛片的网站| 免费看光身美女| a级毛片免费高清观看在线播放| 青春草视频在线免费观看| 欧美精品国产亚洲| 久久婷婷青草| 欧美精品国产亚洲| 久久婷婷青草| 成人特级av手机在线观看| 国产av一区二区精品久久 | 中文字幕免费在线视频6| 亚洲av中文av极速乱| 国产精品久久久久久精品电影小说 | 在线观看av片永久免费下载| 日韩大片免费观看网站| 久久久久精品性色| 全区人妻精品视频| 国产在视频线精品| 色网站视频免费| 国产成人a∨麻豆精品| 女人久久www免费人成看片| 天美传媒精品一区二区| 内射极品少妇av片p| 久久综合国产亚洲精品| 人妻一区二区av| 国产精品女同一区二区软件| kizo精华| 色视频在线一区二区三区| 国产成人精品一,二区| 日韩人妻高清精品专区| 日韩中字成人| 久久99精品国语久久久| a级毛片免费高清观看在线播放| 晚上一个人看的免费电影| 国产女主播在线喷水免费视频网站| 亚洲国产最新在线播放| 国产伦精品一区二区三区四那| 少妇丰满av| 熟女人妻精品中文字幕| 伊人久久精品亚洲午夜| 亚洲精品乱久久久久久| 韩国高清视频一区二区三区| 亚洲av二区三区四区| 精品国产三级普通话版| 2021少妇久久久久久久久久久| 97超碰精品成人国产| 国产精品av视频在线免费观看| 国产亚洲精品久久久com| 免费少妇av软件| 18禁裸乳无遮挡动漫免费视频| 国产深夜福利视频在线观看| 亚州av有码| 亚洲国产毛片av蜜桃av| xxx大片免费视频| h视频一区二区三区| 啦啦啦在线观看免费高清www| 一级爰片在线观看| 国产一区亚洲一区在线观看| 秋霞伦理黄片| 免费av中文字幕在线| 国产av一区二区精品久久 | 久久人人爽人人片av| 在线观看av片永久免费下载| 国产成人aa在线观看| av在线蜜桃| 乱系列少妇在线播放| 91久久精品电影网| 免费黄色在线免费观看| 久久久午夜欧美精品| 99热6这里只有精品| 午夜福利在线在线| 国产无遮挡羞羞视频在线观看| 一区二区三区精品91| 美女高潮的动态| av国产免费在线观看| 久久久久久人妻| 久久久久久伊人网av| 精品亚洲成国产av| 日本猛色少妇xxxxx猛交久久| 97超碰精品成人国产| 18禁在线播放成人免费| 日韩国内少妇激情av| 91精品国产国语对白视频| 国产黄片美女视频| 亚洲综合精品二区| 两个人的视频大全免费| 国产精品国产三级国产专区5o| 全区人妻精品视频| 美女xxoo啪啪120秒动态图| 午夜精品国产一区二区电影| 国产淫片久久久久久久久| 久久久欧美国产精品| 高清日韩中文字幕在线| 韩国av在线不卡| 国产高清不卡午夜福利| av又黄又爽大尺度在线免费看| 国产精品一区www在线观看| 一个人看视频在线观看www免费| 日韩三级伦理在线观看| 国产精品久久久久久久久免| 亚洲av欧美aⅴ国产| 午夜福利影视在线免费观看| 熟女av电影| 一本—道久久a久久精品蜜桃钙片| 久久久久久九九精品二区国产| 亚洲国产最新在线播放| 欧美xxⅹ黑人| 国产淫片久久久久久久久| 一区二区三区精品91| 中文字幕制服av| 精品国产乱码久久久久久小说| 99久久人妻综合| 欧美精品国产亚洲| kizo精华| 午夜福利视频精品| 国产精品一区二区性色av| 日韩 亚洲 欧美在线| 亚洲美女黄色视频免费看| 亚洲美女搞黄在线观看| 日韩三级伦理在线观看| 亚洲国产精品999| 777米奇影视久久| 性色av一级| 国产大屁股一区二区在线视频| 国产一区二区三区av在线| 五月天丁香电影| 精品熟女少妇av免费看| 丝袜脚勾引网站| 男人舔奶头视频| 精品99又大又爽又粗少妇毛片| 国产伦精品一区二区三区视频9| av黄色大香蕉| 黑人猛操日本美女一级片| 久久精品国产亚洲av涩爱| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久电影| 午夜激情福利司机影院| 亚洲欧美日韩无卡精品| 亚洲精品一区蜜桃| 夫妻午夜视频| 一个人看视频在线观看www免费| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美精品自产自拍| 午夜激情福利司机影院| 搡女人真爽免费视频火全软件| 精品视频人人做人人爽| 久久99热这里只频精品6学生| 肉色欧美久久久久久久蜜桃| 久久精品夜色国产| 又大又黄又爽视频免费| 日韩欧美一区视频在线观看 | 视频区图区小说| 免费黄色在线免费观看| 高清午夜精品一区二区三区| 亚洲精品国产av成人精品| 1000部很黄的大片| 啦啦啦啦在线视频资源| av免费在线看不卡| av视频免费观看在线观看| 校园人妻丝袜中文字幕| av.在线天堂| 人妻系列 视频| 日本wwww免费看| h视频一区二区三区| 久久99热这里只频精品6学生| 全区人妻精品视频| 欧美3d第一页| 国产伦精品一区二区三区视频9| 亚洲国产av新网站| 国产乱人偷精品视频| 亚洲国产欧美人成| 一个人看视频在线观看www免费| 亚洲综合精品二区| 91aial.com中文字幕在线观看| 91精品国产九色| 街头女战士在线观看网站| 日本黄色片子视频| 国产精品一及| 国产av一区二区精品久久 | 亚洲av综合色区一区| 精品久久久久久电影网| 久久久欧美国产精品| 日日摸夜夜添夜夜爱| 国产成人精品久久久久久| 啦啦啦啦在线视频资源| 亚洲av男天堂| 人人妻人人看人人澡| 欧美日韩国产mv在线观看视频 | 中文字幕亚洲精品专区| 国产精品人妻久久久影院| 亚洲,欧美,日韩| 美女cb高潮喷水在线观看| 久久精品人妻少妇| 国产精品久久久久久久久免| xxx大片免费视频| 人人妻人人爽人人添夜夜欢视频 | 日韩一区二区三区影片| 久久久久久久久大av| 日韩欧美一区视频在线观看 | 亚洲精品国产av蜜桃| 最近最新中文字幕免费大全7| 欧美日本视频| 国产黄片视频在线免费观看| 在线天堂最新版资源| 色视频www国产| 菩萨蛮人人尽说江南好唐韦庄| 国产午夜精品一二区理论片| 王馨瑶露胸无遮挡在线观看| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲av天美| 亚洲精品国产成人久久av| 性高湖久久久久久久久免费观看| 99热6这里只有精品| 亚洲国产av新网站| 日日摸夜夜添夜夜添av毛片| 亚洲婷婷狠狠爱综合网| 国产极品天堂在线| 久久久久久久久久久丰满| 欧美少妇被猛烈插入视频| 久久女婷五月综合色啪小说| 韩国av在线不卡| 日韩大片免费观看网站| 在线精品无人区一区二区三 | 午夜福利在线在线| 97在线人人人人妻| 亚洲国产精品一区三区| 成年美女黄网站色视频大全免费 | 亚洲国产最新在线播放| 自拍偷自拍亚洲精品老妇| 亚洲不卡免费看| 美女xxoo啪啪120秒动态图| 熟女电影av网| 中文字幕亚洲精品专区| 一级毛片 在线播放| 国产淫语在线视频| 欧美另类一区| 卡戴珊不雅视频在线播放| 国产又色又爽无遮挡免| 国产亚洲最大av| 日韩在线高清观看一区二区三区| 精品久久国产蜜桃| 亚洲国产日韩一区二区| 日韩成人av中文字幕在线观看| 国产一区二区三区综合在线观看 | 校园人妻丝袜中文字幕| 三级经典国产精品| 大香蕉97超碰在线| av在线蜜桃| 最近的中文字幕免费完整| freevideosex欧美| 免费黄网站久久成人精品| 久久久久久久久大av| 国产成人精品一,二区| 国产精品女同一区二区软件| 免费观看av网站的网址| 久久久久久久大尺度免费视频| 日本猛色少妇xxxxx猛交久久|