• <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號
    免费在线观看视频国产中文字幕亚洲 | 欧美日韩精品网址| 亚洲第一av免费看| 少妇裸体淫交视频免费看高清 | 可以免费在线观看a视频的电影网站| 精品少妇久久久久久888优播| 曰老女人黄片| 欧美在线黄色| 大话2 男鬼变身卡| 色综合欧美亚洲国产小说| 国产免费一区二区三区四区乱码| 国产无遮挡羞羞视频在线观看| 欧美精品高潮呻吟av久久| www.999成人在线观看| 亚洲精品自拍成人| 丰满迷人的少妇在线观看| 免费日韩欧美在线观看| 波多野结衣av一区二区av| 人人妻人人澡人人爽人人夜夜| 亚洲av日韩在线播放| 1024视频免费在线观看| 又大又黄又爽视频免费| 久久久久久久久免费视频了| 中文字幕精品免费在线观看视频| 国产精品久久久久成人av| 中文乱码字字幕精品一区二区三区| 欧美成人午夜精品| 亚洲精品第二区| 久久久久国产精品人妻一区二区| 妹子高潮喷水视频| 人人妻,人人澡人人爽秒播 | 一级片免费观看大全| 十八禁人妻一区二区| 亚洲五月婷婷丁香| 亚洲美女黄色视频免费看| 人成视频在线观看免费观看| 纵有疾风起免费观看全集完整版| 狠狠婷婷综合久久久久久88av| 欧美激情高清一区二区三区| 国产成人精品无人区| 国产色视频综合| 91国产中文字幕| 日韩制服骚丝袜av| 久久久久国产一级毛片高清牌| 国产97色在线日韩免费| 日韩中文字幕视频在线看片| 欧美国产精品一级二级三级| av国产精品久久久久影院| 91精品国产国语对白视频| 美女视频免费永久观看网站| 久久久久精品人妻al黑| 国产欧美日韩综合在线一区二区| 日日摸夜夜添夜夜爱| 免费看av在线观看网站| 视频区欧美日本亚洲| 又大又黄又爽视频免费| 国产1区2区3区精品| 老司机影院成人| 国产日韩欧美视频二区| 捣出白浆h1v1| 性少妇av在线| 日本av手机在线免费观看| 色播在线永久视频| 又黄又粗又硬又大视频| 男人爽女人下面视频在线观看| 婷婷色av中文字幕| 亚洲第一青青草原| 18禁黄网站禁片午夜丰满| 亚洲精品av麻豆狂野| 老汉色∧v一级毛片| 亚洲男人天堂网一区| 欧美精品一区二区免费开放| 真人做人爱边吃奶动态| 精品亚洲成国产av| 国产成人精品久久久久久| 国产精品人妻久久久影院| av不卡在线播放| 国产在线观看jvid| 亚洲精品日韩在线中文字幕| 黄色怎么调成土黄色| 欧美日韩亚洲高清精品| 精品久久久久久久毛片微露脸 | 国产成人欧美| 国产精品.久久久| 国产成人av教育| 午夜福利视频精品| 亚洲av成人不卡在线观看播放网 | 亚洲熟女精品中文字幕| 国产男人的电影天堂91| 黄色怎么调成土黄色| 国产免费一区二区三区四区乱码| 国产精品国产三级国产专区5o| 国产在线一区二区三区精| 亚洲av美国av| 成人黄色视频免费在线看| 午夜免费鲁丝| 国产野战对白在线观看| 男女边摸边吃奶| 夫妻午夜视频| 婷婷成人精品国产| 看十八女毛片水多多多| 精品久久久久久电影网| a 毛片基地| 国产不卡av网站在线观看| 男的添女的下面高潮视频| 日本色播在线视频| 欧美日韩一级在线毛片| 一级毛片电影观看| 美女福利国产在线| 又紧又爽又黄一区二区| 精品亚洲成国产av| 啦啦啦 在线观看视频| 免费看十八禁软件| 黄色毛片三级朝国网站| 日本欧美视频一区| 青草久久国产| 午夜日韩欧美国产| 亚洲第一av免费看| 国产精品一国产av| 日韩免费高清中文字幕av| 老熟女久久久| www.av在线官网国产| 一区在线观看完整版| 18禁裸乳无遮挡动漫免费视频| 激情五月婷婷亚洲| 精品少妇久久久久久888优播| 亚洲欧美中文字幕日韩二区| 自线自在国产av| 久久亚洲国产成人精品v| 欧美人与性动交α欧美精品济南到| 丝袜美足系列| av欧美777| 国产精品偷伦视频观看了| 丁香六月欧美| 国产一区有黄有色的免费视频| 国产成人精品久久二区二区91| 久久久精品94久久精品| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 国产人伦9x9x在线观看| 天天躁夜夜躁狠狠躁躁| 男女免费视频国产| 黑人欧美特级aaaaaa片| 国产深夜福利视频在线观看| 美女大奶头黄色视频| 欧美国产精品一级二级三级| 伊人久久大香线蕉亚洲五| 中文字幕另类日韩欧美亚洲嫩草| 激情五月婷婷亚洲| 乱人伦中国视频| www.av在线官网国产| 久久久久久久久久久久大奶| 国产极品粉嫩免费观看在线| 多毛熟女@视频| 夜夜骑夜夜射夜夜干| 悠悠久久av| 18禁国产床啪视频网站| 国产福利在线免费观看视频| 国产高清不卡午夜福利| 亚洲激情五月婷婷啪啪| 两人在一起打扑克的视频| 一级毛片女人18水好多 | 免费观看av网站的网址| 亚洲欧美成人综合另类久久久| 亚洲av片天天在线观看| 国产真人三级小视频在线观看| 汤姆久久久久久久影院中文字幕| 热re99久久国产66热| 亚洲欧美色中文字幕在线| 一级,二级,三级黄色视频| 黄片播放在线免费| 伊人久久大香线蕉亚洲五| 欧美激情高清一区二区三区| 伦理电影免费视频| 国产爽快片一区二区三区| 叶爱在线成人免费视频播放| 国产精品一区二区免费欧美 | 日韩熟女老妇一区二区性免费视频| 日韩大片免费观看网站| 亚洲精品国产区一区二| 久久人妻熟女aⅴ| 丰满人妻熟妇乱又伦精品不卡| 在线观看www视频免费| 女人久久www免费人成看片| 亚洲av男天堂| 秋霞在线观看毛片| 亚洲人成77777在线视频| 在线av久久热| 热re99久久精品国产66热6| 免费高清在线观看日韩| 亚洲精品国产色婷婷电影| 超碰97精品在线观看| 午夜福利免费观看在线| 国产高清videossex| 亚洲av综合色区一区| 90打野战视频偷拍视频| 多毛熟女@视频| 美女视频免费永久观看网站| 国产成人一区二区在线| 国产熟女午夜一区二区三区| 免费日韩欧美在线观看| 国产伦理片在线播放av一区| 老司机午夜十八禁免费视频| 一区二区三区激情视频| 精品福利永久在线观看| 免费一级毛片在线播放高清视频 | 国产精品久久久久久精品古装| 亚洲av电影在线进入| 91麻豆av在线| 999久久久国产精品视频| 久久精品国产亚洲av高清一级| 亚洲av电影在线观看一区二区三区| 亚洲专区中文字幕在线| 天天操日日干夜夜撸| 午夜福利影视在线免费观看| 一本一本久久a久久精品综合妖精| 日本vs欧美在线观看视频| 久久久久国产精品人妻一区二区| 国产有黄有色有爽视频| av国产精品久久久久影院| 搡老岳熟女国产| 大香蕉久久成人网| 我要看黄色一级片免费的| 蜜桃在线观看..| 亚洲av成人精品一二三区| 亚洲精品日本国产第一区| 精品人妻熟女毛片av久久网站| 爱豆传媒免费全集在线观看| 亚洲 国产 在线| 51午夜福利影视在线观看| 欧美日韩视频高清一区二区三区二| 老熟女久久久| 两性夫妻黄色片| 我要看黄色一级片免费的| 久久亚洲国产成人精品v| 两人在一起打扑克的视频| 亚洲成人免费电影在线观看 | 国产精品99久久99久久久不卡| www.av在线官网国产| 又粗又硬又长又爽又黄的视频| 免费在线观看影片大全网站 | 人妻人人澡人人爽人人| 国产成人一区二区在线| √禁漫天堂资源中文www| 国产爽快片一区二区三区| 亚洲一区二区三区欧美精品| 9191精品国产免费久久| 日本黄色日本黄色录像| 欧美精品av麻豆av| 免费少妇av软件| 99精品久久久久人妻精品| 国产精品免费视频内射| 人体艺术视频欧美日本| 欧美日韩亚洲国产一区二区在线观看 | 一个人免费看片子| 99香蕉大伊视频| 欧美精品亚洲一区二区| 丝袜美腿诱惑在线| 好男人视频免费观看在线| 亚洲成人国产一区在线观看 | 亚洲成人国产一区在线观看 | 亚洲,欧美精品.| 国产成人欧美在线观看 | 成年女人毛片免费观看观看9 | 亚洲精品久久久久久婷婷小说| av电影中文网址| 一二三四社区在线视频社区8| 久久久久国产一级毛片高清牌| 欧美日韩亚洲高清精品| 成人手机av| 欧美人与性动交α欧美软件| 大香蕉久久网| 亚洲欧美日韩高清在线视频 | 亚洲精品一区蜜桃| 亚洲色图 男人天堂 中文字幕| 一本久久精品| 看十八女毛片水多多多| 国产高清国产精品国产三级| 日本欧美视频一区| 精品久久久久久电影网| 亚洲精品国产av蜜桃| 亚洲国产精品一区二区三区在线| 国精品久久久久久国模美| 国产欧美日韩一区二区三 | 免费av中文字幕在线| av线在线观看网站| 亚洲五月婷婷丁香| a级毛片黄视频| 久久精品人人爽人人爽视色| 黄片播放在线免费| 国产精品久久久人人做人人爽| 日本午夜av视频| 国产1区2区3区精品| 韩国精品一区二区三区| 久久亚洲国产成人精品v| 久久精品亚洲av国产电影网| 青青草视频在线视频观看| 婷婷色av中文字幕| 亚洲av日韩精品久久久久久密 | 王馨瑶露胸无遮挡在线观看| 欧美亚洲日本最大视频资源| 久久av网站| 欧美日韩成人在线一区二区| xxxhd国产人妻xxx| 好男人电影高清在线观看| 久久久国产精品麻豆| 操美女的视频在线观看| 天天躁日日躁夜夜躁夜夜| a级片在线免费高清观看视频| 亚洲精品乱久久久久久| 精品亚洲乱码少妇综合久久| 一区在线观看完整版| 亚洲精品第二区| 久久99一区二区三区| 亚洲成人手机| 19禁男女啪啪无遮挡网站| 一级毛片电影观看| 久久精品aⅴ一区二区三区四区| 老司机午夜十八禁免费视频| 亚洲av成人不卡在线观看播放网 | 黄色毛片三级朝国网站| 尾随美女入室| 91老司机精品| 亚洲精品久久成人aⅴ小说| 丝袜人妻中文字幕| 丁香六月欧美| 日本欧美视频一区| 18禁黄网站禁片午夜丰满| 另类精品久久| 国产亚洲欧美在线一区二区| 中文精品一卡2卡3卡4更新| 国语对白做爰xxxⅹ性视频网站| 亚洲少妇的诱惑av| 日韩av不卡免费在线播放| 午夜久久久在线观看| 久久久久网色| 久久天堂一区二区三区四区| 国产精品成人在线| 美女福利国产在线| 青春草视频在线免费观看| 一级片免费观看大全| 欧美日韩国产mv在线观看视频| 欧美人与善性xxx| 一级毛片电影观看| 午夜免费观看性视频| 亚洲成国产人片在线观看| 国产爽快片一区二区三区| 超色免费av| 午夜激情av网站| 亚洲九九香蕉| 久久久精品94久久精品| 少妇粗大呻吟视频| 亚洲国产精品一区二区三区在线| 国产免费视频播放在线视频| 成年人免费黄色播放视频| 亚洲情色 制服丝袜| 80岁老熟妇乱子伦牲交| 免费久久久久久久精品成人欧美视频| 国产成人系列免费观看| 狂野欧美激情性xxxx| 女性生殖器流出的白浆| 中文字幕最新亚洲高清| a级片在线免费高清观看视频| 久久ye,这里只有精品| 两个人免费观看高清视频| 精品一区二区三区四区五区乱码 | 亚洲黑人精品在线| 丝袜脚勾引网站| 免费女性裸体啪啪无遮挡网站| 午夜日韩欧美国产| 波野结衣二区三区在线| 欧美乱码精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 美女脱内裤让男人舔精品视频| 精品第一国产精品| 亚洲av美国av| 91老司机精品| 国产欧美亚洲国产| 亚洲 欧美一区二区三区| 精品人妻熟女毛片av久久网站| 国产日韩一区二区三区精品不卡| 免费看不卡的av| 多毛熟女@视频| 精品视频人人做人人爽| 美女扒开内裤让男人捅视频| 午夜福利乱码中文字幕| 亚洲欧洲精品一区二区精品久久久| 日韩电影二区| 国产精品欧美亚洲77777| 国产xxxxx性猛交| 少妇的丰满在线观看| 日本91视频免费播放| 操美女的视频在线观看| 国产免费又黄又爽又色| 熟女少妇亚洲综合色aaa.| 国产又爽黄色视频| 精品欧美一区二区三区在线| 亚洲国产看品久久| 丰满人妻熟妇乱又伦精品不卡| 91麻豆精品激情在线观看国产 | 国产精品成人在线| 黄色 视频免费看| 日本一区二区免费在线视频| av天堂在线播放| 18禁裸乳无遮挡动漫免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产男人的电影天堂91| 蜜桃在线观看..| 亚洲国产精品一区三区| 亚洲成av片中文字幕在线观看| 人人妻人人添人人爽欧美一区卜| 欧美人与善性xxx| 久久99热这里只频精品6学生| 九草在线视频观看| 亚洲欧美一区二区三区久久| 久久久国产精品麻豆| 久久人人爽av亚洲精品天堂| 国产精品香港三级国产av潘金莲 | 99精品久久久久人妻精品| 午夜视频精品福利| 精品一区二区三区四区五区乱码 | 99精品久久久久人妻精品| 中文字幕制服av| 国产欧美日韩综合在线一区二区| 国产在线一区二区三区精| 亚洲伊人久久精品综合| 精品人妻一区二区三区麻豆| 美女扒开内裤让男人捅视频| 另类亚洲欧美激情| 亚洲av电影在线观看一区二区三区| 熟女av电影| 久久久久久久久免费视频了| 性高湖久久久久久久久免费观看| 免费看不卡的av| 多毛熟女@视频| 另类亚洲欧美激情| 亚洲精品久久午夜乱码| 人人妻人人澡人人爽人人夜夜| 美国免费a级毛片| 肉色欧美久久久久久久蜜桃| 免费日韩欧美在线观看| 国产福利在线免费观看视频| 久久人人爽人人片av| 一区在线观看完整版| 国产欧美日韩一区二区三区在线| 日韩一卡2卡3卡4卡2021年| 国产野战对白在线观看| 啦啦啦在线观看免费高清www| 久久精品熟女亚洲av麻豆精品| 国产精品一区二区在线不卡| 十八禁网站网址无遮挡| 丝袜美足系列| 黄频高清免费视频| 下体分泌物呈黄色| 欧美精品一区二区免费开放| 亚洲国产精品国产精品| 91麻豆av在线| 国产精品久久久人人做人人爽| 久久久久久久久免费视频了| 色94色欧美一区二区| 一本大道久久a久久精品| 在线看a的网站| 亚洲五月色婷婷综合| 汤姆久久久久久久影院中文字幕| 中文字幕制服av| 首页视频小说图片口味搜索 | 久久人人爽av亚洲精品天堂| 黄色一级大片看看| 午夜福利影视在线免费观看| 久久精品成人免费网站| 亚洲国产成人一精品久久久| 女人精品久久久久毛片| 日韩精品免费视频一区二区三区| 青青草视频在线视频观看| 视频区欧美日本亚洲| av电影中文网址| 国产成人精品久久二区二区91| 九色亚洲精品在线播放| 久9热在线精品视频| 亚洲少妇的诱惑av| 2018国产大陆天天弄谢| av视频免费观看在线观看| 亚洲精品国产色婷婷电影| 久久鲁丝午夜福利片| 一级毛片我不卡| 男女下面插进去视频免费观看| 777久久人妻少妇嫩草av网站| 18在线观看网站| 飞空精品影院首页| 一级毛片 在线播放| 80岁老熟妇乱子伦牲交| 一级毛片女人18水好多 | 亚洲精品成人av观看孕妇| 欧美人与性动交α欧美精品济南到| 十八禁高潮呻吟视频| 亚洲精品国产一区二区精华液| 91精品三级在线观看| 大片免费播放器 马上看| 国产福利在线免费观看视频| 久久久久久免费高清国产稀缺| 免费不卡黄色视频| 久久亚洲精品不卡| 色视频在线一区二区三区| 国产有黄有色有爽视频| 一级黄片播放器| 久久精品久久久久久噜噜老黄| 午夜影院在线不卡| 国产精品一区二区免费欧美 | 日韩大片免费观看网站| 国产精品国产av在线观看| 在线观看免费高清a一片| 免费女性裸体啪啪无遮挡网站| 亚洲精品成人av观看孕妇| bbb黄色大片| 国产xxxxx性猛交| 久久综合国产亚洲精品| 免费高清在线观看日韩| 一二三四在线观看免费中文在| 美女主播在线视频| 久久狼人影院| 看免费av毛片| 精品免费久久久久久久清纯 | 一级毛片 在线播放| avwww免费| 久久国产精品人妻蜜桃| 五月开心婷婷网| 在现免费观看毛片| 中文字幕人妻熟女乱码| 在线观看一区二区三区激情| 欧美大码av| 午夜福利免费观看在线| 侵犯人妻中文字幕一二三四区| 国产国语露脸激情在线看| 三上悠亚av全集在线观看| 后天国语完整版免费观看| 啦啦啦啦在线视频资源| 1024视频免费在线观看| 欧美精品亚洲一区二区| 男女下面插进去视频免费观看| 老司机深夜福利视频在线观看 | 一级毛片电影观看| 大香蕉久久网| 日本欧美国产在线视频| 亚洲熟女精品中文字幕| 日本欧美视频一区| 脱女人内裤的视频| 波野结衣二区三区在线| 国产成人精品久久二区二区免费| 国产精品一区二区在线观看99| 国产成人免费无遮挡视频| 精品人妻熟女毛片av久久网站| 香蕉国产在线看| 亚洲熟女毛片儿| a级毛片在线看网站| 一级a爱视频在线免费观看| 亚洲av国产av综合av卡| 亚洲欧美激情在线| 国产黄色免费在线视频| 亚洲欧美激情在线| 久久99精品国语久久久| 久久人人爽人人片av| 婷婷色麻豆天堂久久| 亚洲专区国产一区二区| 国产精品国产三级专区第一集| 日韩中文字幕欧美一区二区 | 午夜激情av网站| 色精品久久人妻99蜜桃| 一二三四在线观看免费中文在| 亚洲精品日韩在线中文字幕| 狂野欧美激情性bbbbbb| 丝袜人妻中文字幕| 久久久久久免费高清国产稀缺| 男人操女人黄网站| 欧美日韩亚洲高清精品| av福利片在线| 人人妻,人人澡人人爽秒播 | 久久精品国产综合久久久| 国产伦人伦偷精品视频| 亚洲自偷自拍图片 自拍| 欧美激情高清一区二区三区| 欧美国产精品va在线观看不卡| 精品免费久久久久久久清纯 | 欧美精品高潮呻吟av久久| 精品欧美一区二区三区在线| 精品视频人人做人人爽| 午夜久久久在线观看| 国产精品久久久人人做人人爽| 国产精品欧美亚洲77777| 少妇被粗大的猛进出69影院| 青春草亚洲视频在线观看| 赤兔流量卡办理| 搡老乐熟女国产| 欧美黑人欧美精品刺激| 久久国产亚洲av麻豆专区| 性高湖久久久久久久久免费观看| 国产1区2区3区精品| 亚洲成色77777| 咕卡用的链子| 日韩人妻精品一区2区三区| 国产精品免费视频内射| 亚洲精品久久久久久婷婷小说| 久久精品久久久久久久性| 久久久精品国产亚洲av高清涩受| 亚洲视频免费观看视频| 中文字幕精品免费在线观看视频| 69精品国产乱码久久久| 国产精品国产av在线观看| 最近手机中文字幕大全| 日韩人妻精品一区2区三区| 天堂俺去俺来也www色官网| 别揉我奶头~嗯~啊~动态视频 | 青春草视频在线免费观看| 精品一区二区三区av网在线观看 | 亚洲精品国产区一区二|