• <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电影在线播放| 亚洲欧美日韩另类电影网站| 久久久久久久亚洲中文字幕| 亚洲av二区三区四区| 日本与韩国留学比较| 国产男人的电影天堂91| 能在线免费看毛片的网站| 久久99热6这里只有精品| 国产成人精品一,二区| av福利片在线观看| 免费观看a级毛片全部| 少妇人妻久久综合中文| 只有这里有精品99| 国产乱人偷精品视频| 亚洲av综合色区一区| 内地一区二区视频在线| 国产深夜福利视频在线观看| 亚洲人成网站在线播| 国产欧美日韩一区二区三区在线 | 99久久中文字幕三级久久日本| 菩萨蛮人人尽说江南好唐韦庄| 少妇人妻 视频| 少妇丰满av| 久久久久精品久久久久真实原创| 日本午夜av视频| 26uuu在线亚洲综合色| 少妇被粗大猛烈的视频| 日本色播在线视频| 国产高清不卡午夜福利| 黄色日韩在线| 日本午夜av视频| 日本av手机在线免费观看| 韩国av在线不卡| 秋霞伦理黄片| 午夜福利在线观看免费完整高清在| 99久久精品热视频| 亚洲综合精品二区| 美女中出高潮动态图| 欧美xxⅹ黑人| 国产精品国产三级国产av玫瑰| 成人无遮挡网站| 欧美成人午夜免费资源| 亚洲成色77777| 久久久久精品久久久久真实原创| 91精品国产国语对白视频| 国产一级毛片在线| h日本视频在线播放| 久久久久精品久久久久真实原创| 午夜福利视频精品| 精品亚洲成a人片在线观看| 国产高清不卡午夜福利| 大片电影免费在线观看免费| 五月天丁香电影| 国产一区有黄有色的免费视频| 九九在线视频观看精品| 中文精品一卡2卡3卡4更新| 国语对白做爰xxxⅹ性视频网站| 亚洲精品一区蜜桃| 美女脱内裤让男人舔精品视频| 国产成人免费观看mmmm| 国产在视频线精品| 午夜免费观看性视频| 边亲边吃奶的免费视频| 一边亲一边摸免费视频| 99热全是精品| 国产精品一区二区三区四区免费观看| 99久久中文字幕三级久久日本| 亚洲综合精品二区| 80岁老熟妇乱子伦牲交| 久久99蜜桃精品久久| 综合色丁香网| 久久精品国产a三级三级三级| 亚洲中文av在线| 亚洲欧美成人综合另类久久久| 黄色毛片三级朝国网站 | 2022亚洲国产成人精品| 亚洲精品aⅴ在线观看| 少妇精品久久久久久久| 精品亚洲成国产av| 久久这里有精品视频免费| 精品酒店卫生间| 夜夜看夜夜爽夜夜摸| 肉色欧美久久久久久久蜜桃| 日韩av在线免费看完整版不卡| 如日韩欧美国产精品一区二区三区 | 亚州av有码| 国产视频首页在线观看| 日本av手机在线免费观看| av国产久精品久网站免费入址| 免费观看的影片在线观看| 国产亚洲精品久久久com| 成年人午夜在线观看视频| 色视频www国产| 亚洲精品乱码久久久v下载方式| 99久久精品一区二区三区| 国产免费福利视频在线观看| av一本久久久久| 国产日韩欧美视频二区| 中文资源天堂在线| 久久精品久久久久久噜噜老黄| 热re99久久国产66热| 亚洲美女搞黄在线观看| 精品国产一区二区三区久久久樱花| 亚洲国产精品专区欧美| 一本久久精品| 人人妻人人看人人澡| 久久99蜜桃精品久久| 亚洲色图综合在线观看| 欧美精品高潮呻吟av久久| 亚洲第一区二区三区不卡| 高清av免费在线| 精品人妻熟女av久视频| 人妻制服诱惑在线中文字幕| 日韩伦理黄色片| 99九九在线精品视频 | 午夜久久久在线观看| 亚洲人与动物交配视频| av专区在线播放| 亚洲av国产av综合av卡| 久久 成人 亚洲| 亚洲精品一区蜜桃| 男女边摸边吃奶| 在线观看国产h片| 女的被弄到高潮叫床怎么办| av天堂中文字幕网| av一本久久久久| 国产精品久久久久久精品古装| 久久久久国产网址| 亚洲自偷自拍三级| 国模一区二区三区四区视频| 亚洲人成网站在线播| 尾随美女入室| 欧美亚洲 丝袜 人妻 在线| 日韩中文字幕视频在线看片| 精品卡一卡二卡四卡免费| 久久精品久久久久久久性| 人人妻人人爽人人添夜夜欢视频 | 国产综合精华液| 啦啦啦啦在线视频资源| 大又大粗又爽又黄少妇毛片口| 欧美xxⅹ黑人| 一级片'在线观看视频| 国产免费一区二区三区四区乱码| 九九爱精品视频在线观看| 国产熟女欧美一区二区| 在线免费观看不下载黄p国产| 黄色欧美视频在线观看| 成人亚洲欧美一区二区av| 久久热精品热| 亚洲怡红院男人天堂| 国产日韩一区二区三区精品不卡 | 久久久久久伊人网av| 国产成人a∨麻豆精品| 丝袜在线中文字幕| 国产精品不卡视频一区二区| 成年美女黄网站色视频大全免费 | 国内少妇人妻偷人精品xxx网站| 国产黄色免费在线视频| 七月丁香在线播放| 日韩欧美精品免费久久| 狠狠精品人妻久久久久久综合| 成人午夜精彩视频在线观看| 伊人久久精品亚洲午夜| 国产av码专区亚洲av| 亚洲欧美精品专区久久| 久久久久视频综合| 99九九线精品视频在线观看视频| 插阴视频在线观看视频| 成人二区视频| 两个人免费观看高清视频 | 久久99热这里只频精品6学生| 男人爽女人下面视频在线观看| 特大巨黑吊av在线直播| 欧美高清成人免费视频www| 99视频精品全部免费 在线| 日韩不卡一区二区三区视频在线| 国产免费一区二区三区四区乱码| 少妇的逼水好多| 亚洲精品日韩av片在线观看| 最后的刺客免费高清国语| 赤兔流量卡办理| 国产极品粉嫩免费观看在线 | 99热6这里只有精品| 久久久久久久久久久免费av| 国产男女内射视频| 色视频www国产| 国产成人午夜福利电影在线观看| 国产女主播在线喷水免费视频网站| 国产日韩欧美亚洲二区| 国产精品不卡视频一区二区| 欧美日韩视频精品一区| 日本爱情动作片www.在线观看| 国产精品伦人一区二区| 亚洲在久久综合| av播播在线观看一区| 国产精品一区二区在线不卡| 老司机影院成人| 免费少妇av软件| 天堂中文最新版在线下载| 九色成人免费人妻av| 少妇人妻精品综合一区二区| 欧美激情国产日韩精品一区| 久久久久精品久久久久真实原创| 18禁在线无遮挡免费观看视频| 久久久久久久久久久丰满| 九九在线视频观看精品| 亚洲国产欧美在线一区| 国产精品99久久99久久久不卡 | 国产成人免费观看mmmm| 日韩免费高清中文字幕av| 亚洲图色成人| 国产探花极品一区二区| 日本黄大片高清| 欧美bdsm另类| 精品少妇内射三级| 国产无遮挡羞羞视频在线观看| 成人免费观看视频高清| 街头女战士在线观看网站| 久久久久久久精品精品| 久久精品国产a三级三级三级| av.在线天堂| 岛国毛片在线播放| 国产精品国产三级国产av玫瑰| 老司机亚洲免费影院| 人人妻人人看人人澡| av免费观看日本| 国精品久久久久久国模美| 乱码一卡2卡4卡精品| 久久久久久久久久人人人人人人| 女人精品久久久久毛片| 能在线免费看毛片的网站| 精品亚洲成a人片在线观看| 寂寞人妻少妇视频99o| 午夜久久久在线观看| 韩国av在线不卡| 日韩欧美一区视频在线观看 | 日日撸夜夜添| 中文乱码字字幕精品一区二区三区| 80岁老熟妇乱子伦牲交| 最新的欧美精品一区二区| 狂野欧美激情性xxxx在线观看| 国产中年淑女户外野战色| 熟女av电影| 老女人水多毛片| 国产在线免费精品| 日韩制服骚丝袜av| 国产精品久久久久成人av| 成人美女网站在线观看视频| 国产一区二区三区综合在线观看 | 春色校园在线视频观看| av黄色大香蕉| 黄色一级大片看看| 自拍偷自拍亚洲精品老妇| 久久ye,这里只有精品| 国产老妇伦熟女老妇高清| 日韩电影二区| 深夜a级毛片| 日韩成人av中文字幕在线观看| 成人二区视频| 色94色欧美一区二区| 国产一区二区在线观看日韩| 91久久精品国产一区二区成人| 91成人精品电影| 成年美女黄网站色视频大全免费 | 亚洲综合色惰| 久久97久久精品| 91精品国产九色| 一级毛片aaaaaa免费看小| 一级毛片黄色毛片免费观看视频| 国产精品国产三级国产av玫瑰| 简卡轻食公司| 国产在线视频一区二区| 多毛熟女@视频| 国产亚洲av片在线观看秒播厂| 精华霜和精华液先用哪个| 国产淫片久久久久久久久| 日韩熟女老妇一区二区性免费视频| 国国产精品蜜臀av免费| 男人添女人高潮全过程视频| 十分钟在线观看高清视频www | 日韩av不卡免费在线播放| 国模一区二区三区四区视频| 精品酒店卫生间| 成人漫画全彩无遮挡| av在线老鸭窝| 69精品国产乱码久久久| 亚洲成人一二三区av| 交换朋友夫妻互换小说| 我要看黄色一级片免费的| 成年人免费黄色播放视频 | 乱码一卡2卡4卡精品| 国产熟女午夜一区二区三区 | 精品一品国产午夜福利视频| 日韩三级伦理在线观看| 国产色爽女视频免费观看| 国产欧美日韩综合在线一区二区 | 国产精品久久久久久久电影| 热re99久久国产66热| 久久精品国产亚洲av涩爱| 91久久精品国产一区二区三区| 亚洲国产毛片av蜜桃av| 菩萨蛮人人尽说江南好唐韦庄| 一级av片app| 亚洲av国产av综合av卡| 少妇 在线观看| 你懂的网址亚洲精品在线观看| 日日摸夜夜添夜夜爱| 午夜福利,免费看| 免费av中文字幕在线| 欧美一级a爱片免费观看看| 美女国产视频在线观看| 国产深夜福利视频在线观看| 最新中文字幕久久久久| 91精品国产九色| 自拍偷自拍亚洲精品老妇| 欧美精品亚洲一区二区| 在线观看av片永久免费下载| 99久国产av精品国产电影| 亚洲av二区三区四区| 免费观看性生交大片5| 97超碰精品成人国产| 又粗又硬又长又爽又黄的视频| 久久精品国产鲁丝片午夜精品| 国产一区二区在线观看日韩| 欧美日韩亚洲高清精品| 国产一区二区在线观看日韩| 欧美日韩亚洲高清精品| 久久精品国产自在天天线| 在线观看www视频免费| 亚洲精品乱码久久久久久按摩| 大香蕉97超碰在线| 大话2 男鬼变身卡| 国产淫片久久久久久久久| 免费观看在线日韩| 在现免费观看毛片| 国产女主播在线喷水免费视频网站| 免费观看的影片在线观看| 国产高清不卡午夜福利| 一区二区三区乱码不卡18| .国产精品久久| 久久 成人 亚洲| 成年人免费黄色播放视频 | 日韩欧美一区视频在线观看 | 久久久久久久久久成人| 我要看日韩黄色一级片| 伦精品一区二区三区| 亚洲精品国产成人久久av| h视频一区二区三区| 久久久久视频综合| 欧美97在线视频| 最近最新中文字幕免费大全7| 97在线人人人人妻| 一本—道久久a久久精品蜜桃钙片| 97在线人人人人妻| 国产精品国产三级国产av玫瑰| 国产日韩欧美亚洲二区| 日韩三级伦理在线观看| 男人添女人高潮全过程视频| 日本猛色少妇xxxxx猛交久久| 国产精品麻豆人妻色哟哟久久| 国产伦精品一区二区三区视频9| 99热这里只有精品一区| 26uuu在线亚洲综合色| 特大巨黑吊av在线直播| 亚洲av成人精品一二三区| 99热这里只有精品一区| 亚洲综合色惰| 一级爰片在线观看| 又粗又硬又长又爽又黄的视频| 日本欧美视频一区| 男女边吃奶边做爰视频| 97超碰精品成人国产| 噜噜噜噜噜久久久久久91| 久久久欧美国产精品| av一本久久久久| 国产精品人妻久久久久久| 噜噜噜噜噜久久久久久91| 蜜桃在线观看..| 青春草亚洲视频在线观看| 熟女av电影| 人体艺术视频欧美日本| 国产 一区精品| 三级国产精品欧美在线观看| 亚洲精品国产成人久久av| 伊人久久国产一区二区| 欧美精品一区二区大全| 成人漫画全彩无遮挡| 午夜精品国产一区二区电影| 亚洲经典国产精华液单| 男女国产视频网站| 香蕉精品网在线| 亚洲经典国产精华液单| 五月玫瑰六月丁香| 国产爽快片一区二区三区| 午夜日本视频在线| 最后的刺客免费高清国语| 精品国产一区二区三区久久久樱花| 一级毛片我不卡| 在线天堂最新版资源| 2021少妇久久久久久久久久久| 精品卡一卡二卡四卡免费| 国产精品.久久久| 嫩草影院入口| 插逼视频在线观看| 下体分泌物呈黄色| 只有这里有精品99| 美女视频免费永久观看网站| 久久久久精品久久久久真实原创| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品乱久久久久久| 午夜影院在线不卡| 一级毛片电影观看| 久久久久久久久久成人| av卡一久久| 国产精品一区二区在线观看99| 大片免费播放器 马上看| 精品国产一区二区三区久久久樱花| 卡戴珊不雅视频在线播放| 国产精品一区二区在线观看99| 国模一区二区三区四区视频| 寂寞人妻少妇视频99o| 熟女人妻精品中文字幕| 国产精品偷伦视频观看了| 国产精品久久久久久久久免| 亚洲美女视频黄频| 一级爰片在线观看| 精品一区在线观看国产| 日本黄色日本黄色录像| 午夜91福利影院| 寂寞人妻少妇视频99o| 国产美女午夜福利| 如日韩欧美国产精品一区二区三区 | 亚洲精品乱久久久久久| 国产伦精品一区二区三区四那| 我要看日韩黄色一级片| 黄色一级大片看看| 春色校园在线视频观看| 亚洲激情五月婷婷啪啪| 欧美成人精品欧美一级黄| 国产极品粉嫩免费观看在线 | 国产欧美日韩精品一区二区| 精品少妇内射三级| 插逼视频在线观看| 亚洲人与动物交配视频| 免费观看a级毛片全部| 99热这里只有精品一区| 日韩三级伦理在线观看| 国产成人午夜福利电影在线观看| 一区二区三区四区激情视频| 一级毛片黄色毛片免费观看视频| 日韩 亚洲 欧美在线| 国产真实伦视频高清在线观看| 亚洲国产成人一精品久久久| 免费不卡的大黄色大毛片视频在线观看| 在线观看免费视频网站a站| 一级二级三级毛片免费看| 国产成人精品福利久久| 国内精品宾馆在线| 久久99精品国语久久久| 三上悠亚av全集在线观看 | a级一级毛片免费在线观看| 亚洲国产精品国产精品| 少妇人妻 视频| 如何舔出高潮| 看非洲黑人一级黄片| 亚洲精品国产av蜜桃| 久久久久久久久久久久大奶| 亚洲国产精品成人久久小说| 中文天堂在线官网| 男人爽女人下面视频在线观看| 我的老师免费观看完整版| 女性被躁到高潮视频| 色网站视频免费| 国语对白做爰xxxⅹ性视频网站| 全区人妻精品视频| 男男h啪啪无遮挡| 最近中文字幕2019免费版| 一区二区三区乱码不卡18| 九草在线视频观看| 男男h啪啪无遮挡| 精品酒店卫生间| 精品少妇内射三级| .国产精品久久| 午夜91福利影院| 欧美日韩一区二区视频在线观看视频在线| 一级毛片我不卡| 亚洲经典国产精华液单| av卡一久久| .国产精品久久| 丰满少妇做爰视频| .国产精品久久| 六月丁香七月| 老司机影院毛片| 热re99久久国产66热| 国产男女超爽视频在线观看| 久久 成人 亚洲| 女性生殖器流出的白浆| 国产av国产精品国产| 大香蕉久久网| 国产精品熟女久久久久浪| 久久99精品国语久久久| 精品人妻一区二区三区麻豆| 日韩强制内射视频| 欧美最新免费一区二区三区| a 毛片基地| 久久国产乱子免费精品| 在线观看一区二区三区激情| 激情五月婷婷亚洲| 亚洲经典国产精华液单| 狂野欧美白嫩少妇大欣赏| 久久久久精品久久久久真实原创| 免费人成在线观看视频色| 久久青草综合色| 中文字幕人妻丝袜制服| 国产乱来视频区| 亚洲怡红院男人天堂| 人妻 亚洲 视频| av在线播放精品| 久久99一区二区三区| 亚洲精品456在线播放app| 爱豆传媒免费全集在线观看| 国产毛片在线视频| 中文字幕精品免费在线观看视频 | 国产精品一区二区三区四区免费观看| 亚洲精品日本国产第一区| 日本av免费视频播放| 少妇精品久久久久久久| 国产日韩一区二区三区精品不卡 | 高清黄色对白视频在线免费看 | 亚洲精品乱久久久久久| 人人妻人人爽人人添夜夜欢视频 | 九草在线视频观看| 下体分泌物呈黄色| 免费播放大片免费观看视频在线观看| 99久久精品一区二区三区| 91精品国产国语对白视频| 91久久精品电影网| 蜜桃久久精品国产亚洲av| √禁漫天堂资源中文www| 性高湖久久久久久久久免费观看| av女优亚洲男人天堂| 国产一区二区在线观看日韩| 成年女人在线观看亚洲视频| 卡戴珊不雅视频在线播放| 久久热精品热| 久久精品国产a三级三级三级| 丰满乱子伦码专区| 成人亚洲欧美一区二区av| 国产日韩欧美视频二区| 国产精品国产三级国产av玫瑰| 国产伦在线观看视频一区| 春色校园在线视频观看| 如日韩欧美国产精品一区二区三区 | 一级黄片播放器| 亚洲国产av新网站| 国产午夜精品一二区理论片| av国产久精品久网站免费入址| 2022亚洲国产成人精品| 一个人免费看片子| 人人澡人人妻人| 亚洲国产精品一区二区三区在线| 极品教师在线视频| 啦啦啦在线观看免费高清www| 黄色配什么色好看| 国产精品99久久99久久久不卡 | 国产免费视频播放在线视频| 男女啪啪激烈高潮av片| 99久久人妻综合| av一本久久久久| 99视频精品全部免费 在线| 在线观看www视频免费| 777米奇影视久久| 精品酒店卫生间| 少妇 在线观看| 我的女老师完整版在线观看| 久久久久久久久久成人| 亚洲国产色片| 国产一区二区在线观看av| 18+在线观看网站| 99久久精品国产国产毛片| 深夜a级毛片| 成年女人在线观看亚洲视频| 精品卡一卡二卡四卡免费| 国产黄频视频在线观看| 国产色爽女视频免费观看| 亚洲色图综合在线观看| 男女边摸边吃奶| 在线观看人妻少妇| 日本爱情动作片www.在线观看| 伊人亚洲综合成人网| 丝瓜视频免费看黄片| 欧美性感艳星| 边亲边吃奶的免费视频| 九九爱精品视频在线观看| 青春草亚洲视频在线观看| 国产成人精品久久久久久| 国产爽快片一区二区三区| 久久久久人妻精品一区果冻| 亚洲欧美精品自产自拍| 男女国产视频网站| 免费观看av网站的网址| 在线亚洲精品国产二区图片欧美 | 久久久国产一区二区| 99re6热这里在线精品视频| 一区二区三区乱码不卡18| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久久久免| 一本色道久久久久久精品综合| 午夜免费男女啪啪视频观看| 国产免费视频播放在线视频| 我的老师免费观看完整版| 国产精品伦人一区二区| 女人精品久久久久毛片|