• <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玫瑰| 久久99热这里只频精品6学生| 国产黄色视频一区二区在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲人成77777在线视频| 久久99精品国语久久久| 欧美性感艳星| 国产黄色视频一区二区在线观看| 日韩一区二区视频免费看| 精品国产露脸久久av麻豆| 日韩中字成人| 各种免费的搞黄视频| 在线观看三级黄色| 精品国产乱码久久久久久小说| 91国产中文字幕| √禁漫天堂资源中文www| 大香蕉97超碰在线| av在线播放免费不卡| 久久人妻福利社区极品人妻图片| 精品熟女少妇八av免费久了| av国产精品久久久久影院| 亚洲成人免费电影在线观看| 欧美老熟妇乱子伦牲交| 国产精品偷伦视频观看了| 国产成人啪精品午夜网站| av不卡在线播放| 久久久久久久国产电影| 成人av一区二区三区在线看| 热re99久久精品国产66热6| 日韩免费av在线播放| 日韩欧美三级三区| 日本av免费视频播放| 国产人伦9x9x在线观看| 欧美成人午夜精品| 高清视频免费观看一区二区| 亚洲五月婷婷丁香| 我的亚洲天堂| 热re99久久精品国产66热6| 国产伦理片在线播放av一区| 国产色视频综合| 成人18禁高潮啪啪吃奶动态图| 色婷婷久久久亚洲欧美| 久久国产精品影院| 窝窝影院91人妻| 精品高清国产在线一区| 成人影院久久| 国产男女内射视频| 国产成人欧美在线观看 | 国产成人欧美在线观看 | 精品一区二区三区视频在线观看免费 | 成人三级做爰电影| 欧美性长视频在线观看| 久久亚洲精品不卡| 汤姆久久久久久久影院中文字幕| 两个人看的免费小视频| 悠悠久久av| 大香蕉久久网| 啪啪无遮挡十八禁网站| 男女免费视频国产| 国产黄色免费在线视频| 久久人人爽av亚洲精品天堂| 少妇粗大呻吟视频| 一个人免费在线观看的高清视频| 18禁裸乳无遮挡动漫免费视频| 国产亚洲精品一区二区www | 狠狠狠狠99中文字幕| 可以免费在线观看a视频的电影网站| av电影中文网址| 国产深夜福利视频在线观看| 亚洲国产看品久久| 日韩欧美一区视频在线观看| 亚洲av第一区精品v没综合| 777米奇影视久久| 国产在视频线精品| 一区福利在线观看| 亚洲五月婷婷丁香| 亚洲色图综合在线观看| 日韩欧美国产一区二区入口| 久久久国产精品麻豆| 欧美日韩一级在线毛片| 真人做人爱边吃奶动态| 亚洲成人免费av在线播放| 99国产综合亚洲精品| 国产区一区二久久| 咕卡用的链子| 欧美亚洲日本最大视频资源| 黄网站色视频无遮挡免费观看| 中文字幕精品免费在线观看视频| 国产精品1区2区在线观看. | 亚洲精品国产精品久久久不卡| 久久久久久免费高清国产稀缺| av视频免费观看在线观看| 51午夜福利影视在线观看| 18禁黄网站禁片午夜丰满| 国产成人啪精品午夜网站| 日本一区二区免费在线视频| 久久香蕉激情| 日本av免费视频播放| 国产97色在线日韩免费| 熟女少妇亚洲综合色aaa.| 国产成人av激情在线播放| 夜夜骑夜夜射夜夜干| 一区二区三区国产精品乱码| 热re99久久国产66热| 一本综合久久免费| 国产精品久久久久成人av| 欧美日韩国产mv在线观看视频| 黄色成人免费大全| 精品高清国产在线一区| 一区二区三区国产精品乱码| 国产又色又爽无遮挡免费看| 桃红色精品国产亚洲av| 热re99久久国产66热| videos熟女内射| 狠狠婷婷综合久久久久久88av| 成年女人毛片免费观看观看9 | 欧美日韩av久久| 亚洲,欧美精品.| 欧美人与性动交α欧美精品济南到| 窝窝影院91人妻| 精品人妻1区二区| 露出奶头的视频| 十分钟在线观看高清视频www| 岛国在线观看网站| 午夜福利,免费看| 女警被强在线播放| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三区视频在线观看免费 | 亚洲av电影在线进入| 老司机影院毛片| 男女免费视频国产| 国产黄频视频在线观看| 久久久国产一区二区| 天天躁日日躁夜夜躁夜夜| 日韩中文字幕欧美一区二区| 久热这里只有精品99| 19禁男女啪啪无遮挡网站| 天天操日日干夜夜撸| 最新在线观看一区二区三区| 麻豆乱淫一区二区| 亚洲性夜色夜夜综合| 人妻久久中文字幕网| 少妇 在线观看| 高清视频免费观看一区二区| 亚洲免费av在线视频| 9热在线视频观看99| 亚洲,欧美精品.| 18在线观看网站| 欧美一级毛片孕妇| 一级毛片精品| 色婷婷久久久亚洲欧美| 日韩制服丝袜自拍偷拍| 久久国产精品男人的天堂亚洲| 久久久久精品人妻al黑| 午夜福利在线免费观看网站| 亚洲精品av麻豆狂野| 黄色怎么调成土黄色| 国产免费视频播放在线视频| 一本综合久久免费| 黄色 视频免费看| 91成人精品电影| 欧美黑人欧美精品刺激| 最近最新中文字幕大全免费视频| 日韩一卡2卡3卡4卡2021年| 欧美国产精品va在线观看不卡| 韩国精品一区二区三区| 18禁国产床啪视频网站| 最近最新中文字幕大全免费视频| 日韩有码中文字幕| 人人妻人人添人人爽欧美一区卜| 天天影视国产精品| 久久午夜亚洲精品久久| 不卡av一区二区三区| 黄片小视频在线播放| 国产国语露脸激情在线看| 国产精品熟女久久久久浪| 宅男免费午夜| 视频在线观看一区二区三区| 在线播放国产精品三级| 亚洲综合色网址| 国产在线观看jvid| 久久久久久久久久久久大奶| 成人影院久久| 999精品在线视频| www.999成人在线观看| 欧美日韩一级在线毛片| 真人做人爱边吃奶动态| 在线观看一区二区三区激情| 亚洲欧美一区二区三区黑人| 国产亚洲欧美精品永久| 欧美精品亚洲一区二区| 他把我摸到了高潮在线观看 | 日本av手机在线免费观看| 丝袜在线中文字幕| 日韩三级视频一区二区三区| 国产av又大| 国产99久久九九免费精品| 中文亚洲av片在线观看爽 | 国产欧美亚洲国产| 一区在线观看完整版| 丁香六月欧美| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费视频网站a站| 在线看a的网站| 亚洲专区国产一区二区| 无限看片的www在线观看| 飞空精品影院首页| 精品人妻熟女毛片av久久网站| 亚洲av国产av综合av卡| 人妻久久中文字幕网| netflix在线观看网站| 国产精品免费视频内射| 巨乳人妻的诱惑在线观看| 午夜福利免费观看在线| 嫁个100分男人电影在线观看| 精品午夜福利视频在线观看一区 | 一本久久精品| av超薄肉色丝袜交足视频| 色婷婷av一区二区三区视频| 最近最新中文字幕大全电影3 | 黑人欧美特级aaaaaa片| tube8黄色片| 中文欧美无线码| 老司机亚洲免费影院| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 超碰97精品在线观看| 亚洲av第一区精品v没综合| 日本av手机在线免费观看| 日韩 欧美 亚洲 中文字幕| 99精品久久久久人妻精品| 下体分泌物呈黄色| 国产有黄有色有爽视频| 99国产精品99久久久久| 欧美黄色片欧美黄色片| 免费av中文字幕在线| 亚洲av成人不卡在线观看播放网| 成年人黄色毛片网站| 欧美激情久久久久久爽电影 | 老司机亚洲免费影院| 久久久久精品人妻al黑| 久久ye,这里只有精品| 亚洲美女黄片视频| 一进一出抽搐动态| 激情视频va一区二区三区| 免费观看av网站的网址| 色播在线永久视频| 亚洲成人国产一区在线观看| 久久久久视频综合| 国产日韩欧美亚洲二区| 波多野结衣一区麻豆| 色老头精品视频在线观看| av有码第一页| 精品少妇久久久久久888优播| 亚洲国产看品久久| 一进一出抽搐动态| 超碰97精品在线观看| cao死你这个sao货| 国产成+人综合+亚洲专区| 捣出白浆h1v1| 成人18禁高潮啪啪吃奶动态图| 亚洲第一av免费看| 一本色道久久久久久精品综合| 成人特级黄色片久久久久久久 | 亚洲av片天天在线观看| 久久久水蜜桃国产精品网| 男女高潮啪啪啪动态图| 99re在线观看精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 免费一级毛片在线播放高清视频 | 日韩免费av在线播放| 午夜福利一区二区在线看| 久久久久久久精品吃奶| 日本撒尿小便嘘嘘汇集6| 国产精品电影一区二区三区 | 天堂俺去俺来也www色官网| av在线播放免费不卡| 亚洲七黄色美女视频| 日韩中文字幕视频在线看片| 女人爽到高潮嗷嗷叫在线视频| 少妇精品久久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 热re99久久精品国产66热6| 精品福利永久在线观看| 欧美精品啪啪一区二区三区| 欧美日韩亚洲综合一区二区三区_| 黄色怎么调成土黄色| 少妇被粗大的猛进出69影院| 午夜久久久在线观看| 女人高潮潮喷娇喘18禁视频| 日本wwww免费看| 亚洲精品国产区一区二| 久久精品亚洲av国产电影网| 一边摸一边做爽爽视频免费| 这个男人来自地球电影免费观看| 久久影院123| 一边摸一边做爽爽视频免费| 国产在线精品亚洲第一网站| 精品一区二区三区视频在线观看免费 | av有码第一页| 日本撒尿小便嘘嘘汇集6| 欧美 亚洲 国产 日韩一| 日韩视频在线欧美| 最黄视频免费看| 成人特级黄色片久久久久久久 | 美国免费a级毛片| 国产欧美日韩精品亚洲av| svipshipincom国产片| 国产在线观看jvid| 中文字幕色久视频| 不卡av一区二区三区| 热re99久久精品国产66热6| 色婷婷久久久亚洲欧美| 亚洲精品一二三| 欧美日韩成人在线一区二区| 啦啦啦在线免费观看视频4| 新久久久久国产一级毛片| 高清黄色对白视频在线免费看| 亚洲欧美一区二区三区久久| av欧美777| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美久久黑人一区二区| www.自偷自拍.com| 精品国产一区二区三区久久久樱花| 亚洲欧洲精品一区二区精品久久久| 国产日韩欧美在线精品| 国产单亲对白刺激| 免费黄频网站在线观看国产| 成人国语在线视频| 亚洲成人免费av在线播放| 日本wwww免费看| 欧美 日韩 精品 国产| 两性夫妻黄色片| 中文亚洲av片在线观看爽 | 亚洲国产av影院在线观看| av不卡在线播放| 亚洲精品一二三| 国产福利在线免费观看视频| 首页视频小说图片口味搜索| 久久人妻福利社区极品人妻图片| 欧美精品一区二区免费开放| 亚洲欧洲日产国产| 下体分泌物呈黄色| 国产伦理片在线播放av一区| 亚洲熟女精品中文字幕| 国精品久久久久久国模美| 中文字幕制服av| 中国美女看黄片| 国产老妇伦熟女老妇高清| 黄频高清免费视频| 欧美性长视频在线观看| 亚洲成人国产一区在线观看| 亚洲伊人久久精品综合| 热99国产精品久久久久久7| h视频一区二区三区| 免费观看a级毛片全部| 午夜福利影视在线免费观看| 欧美精品高潮呻吟av久久| av超薄肉色丝袜交足视频| 真人做人爱边吃奶动态| 日韩大片免费观看网站| 亚洲va日本ⅴa欧美va伊人久久| 久久毛片免费看一区二区三区| 黄色片一级片一级黄色片| 精品国产超薄肉色丝袜足j| 纵有疾风起免费观看全集完整版| 青青草视频在线视频观看| 欧美av亚洲av综合av国产av| av福利片在线| 久久国产精品影院| 男女免费视频国产| 99九九在线精品视频| 日本黄色日本黄色录像| 国产精品免费大片| 丝袜在线中文字幕| a级毛片在线看网站| 天天躁日日躁夜夜躁夜夜| 在线天堂中文资源库| 欧美精品av麻豆av| 欧美精品高潮呻吟av久久| 精品国产乱码久久久久久小说| 黄色视频在线播放观看不卡| 日本a在线网址| 精品少妇久久久久久888优播| 亚洲第一av免费看| bbb黄色大片| av免费在线观看网站| 99国产综合亚洲精品| www.熟女人妻精品国产| 丝袜喷水一区| 国产国语露脸激情在线看| 在线 av 中文字幕| 久久中文字幕一级| 国产片内射在线| 美女午夜性视频免费| 国产成人免费观看mmmm| 建设人人有责人人尽责人人享有的| 国产有黄有色有爽视频| 中文字幕高清在线视频| 人妻一区二区av| www.精华液| 飞空精品影院首页| 久久精品亚洲av国产电影网| 久久久久网色| 亚洲精品粉嫩美女一区| 日本撒尿小便嘘嘘汇集6| 国产亚洲一区二区精品| 美女高潮喷水抽搐中文字幕| 大香蕉久久成人网| 我要看黄色一级片免费的| 欧美黑人欧美精品刺激| 亚洲人成伊人成综合网2020| av免费在线观看网站| 19禁男女啪啪无遮挡网站| 夫妻午夜视频| 每晚都被弄得嗷嗷叫到高潮| 在线看a的网站| 久久影院123| 国产野战对白在线观看| av又黄又爽大尺度在线免费看| 国产又色又爽无遮挡免费看| 久久影院123| 捣出白浆h1v1| 在线观看舔阴道视频| 18禁美女被吸乳视频| 一边摸一边抽搐一进一出视频| 中国美女看黄片| 精品一区二区三区视频在线观看免费 | 欧美乱码精品一区二区三区| 51午夜福利影视在线观看| 曰老女人黄片| 国产黄频视频在线观看| 大型黄色视频在线免费观看| 久久精品国产99精品国产亚洲性色 | 性色av乱码一区二区三区2| 在线十欧美十亚洲十日本专区| 久久亚洲真实| 国产一区二区三区视频了| 午夜老司机福利片| 免费少妇av软件| 成年女人毛片免费观看观看9 | 九色亚洲精品在线播放| 91精品三级在线观看| 午夜激情久久久久久久| 9热在线视频观看99| 亚洲av日韩精品久久久久久密| 久久人人爽av亚洲精品天堂| 欧美老熟妇乱子伦牲交| 丁香六月欧美| 18禁美女被吸乳视频| 国产精品免费视频内射| 中文字幕人妻熟女乱码| 无遮挡黄片免费观看| 美女扒开内裤让男人捅视频| 深夜精品福利| 亚洲欧美一区二区三区久久| 母亲3免费完整高清在线观看| 建设人人有责人人尽责人人享有的| 亚洲一卡2卡3卡4卡5卡精品中文| 精品亚洲成a人片在线观看| 欧美日韩亚洲国产一区二区在线观看 | 日韩欧美国产一区二区入口| 亚洲综合色网址| 欧美老熟妇乱子伦牲交| 国产av国产精品国产| 日韩欧美免费精品| 侵犯人妻中文字幕一二三四区| 欧美日韩视频精品一区| 男女午夜视频在线观看| 正在播放国产对白刺激| 国产三级黄色录像| 欧美精品啪啪一区二区三区| 国产色视频综合| av线在线观看网站| 法律面前人人平等表现在哪些方面| 欧美老熟妇乱子伦牲交| 精品卡一卡二卡四卡免费| 18禁美女被吸乳视频| 亚洲精品成人av观看孕妇| 免费观看a级毛片全部| 涩涩av久久男人的天堂| 国产精品 欧美亚洲| 日本黄色日本黄色录像| 我要看黄色一级片免费的| 久久99一区二区三区| 老司机亚洲免费影院| 性色av乱码一区二区三区2| 黑人操中国人逼视频| 视频区欧美日本亚洲| 国产日韩一区二区三区精品不卡| 搡老岳熟女国产| 国产区一区二久久| 乱人伦中国视频| 国产欧美日韩综合在线一区二区| 麻豆成人av在线观看| 麻豆国产av国片精品| 少妇精品久久久久久久| 国产aⅴ精品一区二区三区波| 国产免费现黄频在线看| 欧美日韩精品网址| 中文字幕人妻熟女乱码| 日韩视频在线欧美| 脱女人内裤的视频| 狠狠狠狠99中文字幕| 日韩一区二区三区影片| 菩萨蛮人人尽说江南好唐韦庄| 一区福利在线观看| 久久国产亚洲av麻豆专区| 国产成人精品久久二区二区91| 妹子高潮喷水视频| 国产成人欧美在线观看 | 久久久久国内视频| 天天躁夜夜躁狠狠躁躁| 国产精品免费大片| 国产精品.久久久| 欧美激情久久久久久爽电影 | 欧美av亚洲av综合av国产av| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 亚洲精品国产色婷婷电影| 欧美黑人精品巨大| 国产精品欧美亚洲77777| 久久狼人影院| av又黄又爽大尺度在线免费看| 亚洲欧美色中文字幕在线| 日韩大片免费观看网站| 国产日韩欧美亚洲二区| 女人精品久久久久毛片| 午夜福利乱码中文字幕| 我的亚洲天堂| 香蕉丝袜av| 亚洲精品在线美女| 少妇的丰满在线观看| 男女床上黄色一级片免费看| 亚洲自偷自拍图片 自拍| 久久九九热精品免费| 免费黄频网站在线观看国产| 国产日韩一区二区三区精品不卡| 国产有黄有色有爽视频| 欧美日韩黄片免| 精品久久蜜臀av无| 成年人免费黄色播放视频| 超碰成人久久| 久久香蕉激情| 精品一品国产午夜福利视频| 亚洲中文日韩欧美视频| 一边摸一边抽搐一进一小说 | 极品教师在线免费播放| 久久国产精品影院| 国产区一区二久久| 99精品欧美一区二区三区四区| 亚洲自偷自拍图片 自拍| 欧美乱妇无乱码| 在线观看一区二区三区激情| 我的亚洲天堂| 国产成人精品久久二区二区91| 中文字幕制服av| 法律面前人人平等表现在哪些方面| 色94色欧美一区二区| 侵犯人妻中文字幕一二三四区| 99热网站在线观看| 国产精品国产av在线观看| 99国产精品99久久久久| 视频在线观看一区二区三区| 99国产极品粉嫩在线观看| 国产高清videossex| 香蕉丝袜av| 蜜桃国产av成人99| 国产男女内射视频| 久久久国产精品麻豆| 一边摸一边做爽爽视频免费| 一本—道久久a久久精品蜜桃钙片| 免费在线观看完整版高清| 99国产精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 男女无遮挡免费网站观看| 日日摸夜夜添夜夜添小说| 久久精品91无色码中文字幕| av天堂在线播放| 国产精品一区二区在线不卡| 国内毛片毛片毛片毛片毛片| 不卡一级毛片| 亚洲avbb在线观看| 一个人免费在线观看的高清视频| 丰满饥渴人妻一区二区三| 婷婷丁香在线五月| 90打野战视频偷拍视频| 亚洲欧美激情在线| 欧美成人免费av一区二区三区 | 亚洲黑人精品在线| 男人舔女人的私密视频| 国产高清国产精品国产三级| 久久中文字幕一级| 久久久久久人人人人人| 久久久精品94久久精品| 亚洲九九香蕉| 一级a爱视频在线免费观看| 香蕉国产在线看| 久久久久久久国产电影| 老司机午夜福利在线观看视频 | 亚洲精品自拍成人| 一边摸一边做爽爽视频免费| 亚洲色图 男人天堂 中文字幕| 亚洲精品在线观看二区| 欧美日韩亚洲国产一区二区在线观看 | 丁香六月欧美| 久久中文字幕一级| 一级a爱视频在线免费观看| 日本av免费视频播放| 国产深夜福利视频在线观看| 两个人看的免费小视频| 夜夜爽天天搞|