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

    基于云模型的冬小麥氣候適宜度評(píng)價(jià)方法*
    ——以安徽省宿州市為例

    2016-09-10 06:15:19德周文鱗楊霏云
    關(guān)鍵詞:宿州市日照時(shí)數(shù)冬小麥

    李 德周文鱗楊霏云

    (1. 安徽省宿州市氣象局 宿州 234000; 2. 安徽省合肥市氣象局 合肥 236001;3. 中國(guó)氣象局氣象干部培訓(xùn)中心 北京 100081)

    基于云模型的冬小麥氣候適宜度評(píng)價(jià)方法*
    ——以安徽省宿州市為例

    李德1周文鱗2楊霏云3

    (1. 安徽省宿州市氣象局 宿州 234000; 2. 安徽省合肥市氣象局 合肥 236001;3. 中國(guó)氣象局氣象干部培訓(xùn)中心 北京 100081)

    為建立冬小麥氣候適宜度量化評(píng)價(jià)方法,本文基于云模型理論,依據(jù)光、溫和水界限指標(biāo),采用“3En”法則確定云參數(shù),建立日照、氣溫和降水對(duì)小麥生長(zhǎng)影響的云模型。運(yùn)用積分回歸法,確定權(quán)重系數(shù),采用加權(quán)綜合法和幾何平均法,確定不同生育期和全生育期氣候適宜度,利用1954—2013年安徽省宿州市各縣(區(qū))冬小麥單產(chǎn)和 1995—2013年觀測(cè)地段產(chǎn)量因素等資料進(jìn)行檢驗(yàn)。結(jié)果表明,日照適宜度可用左半云,氣溫和降水適宜度可用梯形云來(lái)表達(dá)。計(jì)算的冬小麥全生育期氣候適宜度,與宿州市各縣(區(qū))冬小麥氣候產(chǎn)量呈顯著或極顯著正相關(guān)關(guān)系; 與觀測(cè)地段的冬小麥氣候產(chǎn)量、千粒重、每穗籽粒數(shù)和乳熟期株高呈顯著正相關(guān),相關(guān)系數(shù)分別為0.588 0(P<0.01)、0.756 1(P<0.01)、0.670 7(P<0.01)和0.464 3(P<0.05)。返青—拔節(jié)期、抽穗—乳熟期2個(gè)時(shí)期的氣候適宜度與單位面積穗數(shù)、每穗籽粒數(shù)的相關(guān)系數(shù)分別為0.558 9(P<0.05)、0.710 7(P<0.01)和0.736 1(P<0.01)、0.744 2(P<0.01),拔節(jié)—抽穗期氣候適宜度與單位面積穗數(shù)的相關(guān)系數(shù)為0.649 8(P<0.01)。1954—2013年宿州市日照與降水適宜度以每10年0.005和0.008的速度降低,氣溫適宜度以每10年0.028的速率升高。研究結(jié)果可作為評(píng)價(jià)宿州冬小麥對(duì)氣候條件的適應(yīng)性及制定相應(yīng)策略的參考依據(jù)。

    冬小麥 云模型 “3En”法則 氣候適宜度 氣候產(chǎn)量 觀測(cè)產(chǎn)量 氣候因子 生育期 宿州

    農(nóng)業(yè)氣候適宜性評(píng)價(jià)問(wèn)題,是農(nóng)業(yè)氣象學(xué)的主要研究?jī)?nèi)容之一[1-2],特別是在全球氣候變暖背景下,客觀定量評(píng)價(jià)氣候條件對(duì)農(nóng)作物生長(zhǎng)發(fā)育的影響,為農(nóng)業(yè)生產(chǎn)采取相應(yīng)措施應(yīng)對(duì)或適應(yīng)氣候變化具有重要意義[3-5]。如何客觀定量評(píng)價(jià)氣候條件對(duì)農(nóng)作物生長(zhǎng)發(fā)育影響的優(yōu)劣程度,國(guó)內(nèi)外學(xué)者進(jìn)行了大量研究[1]。在評(píng)價(jià)方法上,一是采用某時(shí)段若干氣候要素特征值的比較法來(lái)研判氣候條件對(duì)農(nóng)作物生長(zhǎng)發(fā)育適宜程度的傳統(tǒng)方法[3,6-7],二是針對(duì)傳統(tǒng)方法不能客觀反映農(nóng)作物與氣候條件之間相互作用與影響過(guò)程具有不確定性的特點(diǎn),基于模糊數(shù)學(xué)理論,提出用氣候適宜度的概念[8-9]來(lái)評(píng)價(jià)光溫水各要素及其綜合影響,并在作物氣候適宜性評(píng)價(jià)[10-12]、作物產(chǎn)量評(píng)估[13-14]與發(fā)育期預(yù)測(cè)[15]、農(nóng)業(yè)生態(tài)氣候區(qū)劃、種植制度與種植區(qū)域變化[16-18]和不同時(shí)間尺度氣候適宜度模型及其相關(guān)參數(shù)客觀化等方面進(jìn)行了研究[2,19]。然而,自氣候適宜度概念提出以來(lái),一直存在著生物學(xué)意義不清的缺陷。同時(shí),一旦用一個(gè)隸屬度函數(shù)來(lái)描述氣候要素對(duì)農(nóng)作物影響的模糊性,其模糊概念就被強(qiáng)行納入到精確計(jì)算之中而不再有模糊性。隸屬函數(shù)即轉(zhuǎn)化為僅是物理意義上的確定關(guān)系而非統(tǒng)計(jì)學(xué)上的隨機(jī)性。其實(shí),就連模糊集理論創(chuàng)始人Zadeh本人也只是用定性推理方法近似指定隸屬度函數(shù)[20]。因此,建立能反映農(nóng)作物-氣候相互作用過(guò)程中的隨機(jī)性和模糊性的氣候適宜度評(píng)價(jià)模型[1,8],是農(nóng)業(yè)氣象學(xué)需要解決的重大問(wèn)題之一。另一方面,由于農(nóng)作物生長(zhǎng)環(huán)境的復(fù)雜性和多變性,至今人們所能掌握的僅是某些關(guān)鍵期或時(shí)段的一些定性或離散化定量指標(biāo),很難掌握精細(xì)到每一時(shí)刻(即時(shí)間間隔足夠小)氣候條件對(duì)農(nóng)作物影響的精確指標(biāo),因而也就不能確定氣候條件對(duì)其影響的精確物理關(guān)系。因此,氣候適宜度評(píng)價(jià)模型構(gòu)建需要另辟蹊徑。

    云模型是一種定性知識(shí)描述與定性概念及其定量數(shù)值之間的轉(zhuǎn)換模型[20-21],目前已在圖像分類[22]、土地生態(tài)評(píng)價(jià)[23]、防洪減災(zāi)能力評(píng)估[24]、水庫(kù)誘發(fā)地震風(fēng)險(xiǎn)[25]和麥田噴藥氣象適宜等級(jí)評(píng)定[26]等方面有了應(yīng)用,但在作物氣候適宜度評(píng)價(jià)方面應(yīng)用較少。本文以位于黃淮冬麥區(qū)淮北平原中部的安徽省宿州市冬小麥生長(zhǎng)期間光、溫、水逐日氣候要素測(cè)量值,結(jié)合不同生育時(shí)段的光、溫、水氣候指標(biāo),基于云模型理論的特點(diǎn)與優(yōu)勢(shì)[20-21],確定不同生育期日照、溫度和降水云模型參數(shù),結(jié)合計(jì)算機(jī)仿真技術(shù),得到不同時(shí)期氣候要素影響的云隸屬度,用于評(píng)價(jià)氣候條件對(duì)冬小麥生長(zhǎng)發(fā)育影響的優(yōu)劣程度,從而為科學(xué)鑒定冬小麥對(duì)氣候條件的適應(yīng)性和制定相應(yīng)種植策略提供依據(jù)。

    1 研究區(qū)概況與研究方法

    1.1 研究地概況與資料來(lái)源

    安徽省宿州市位于淮北平原中部,介于 33.3°~34.7°N、116.2°~118.2°E之間,屬亞熱帶與暖溫帶氣候過(guò)渡區(qū),其下轄泗縣、靈璧縣、蕭縣、碭山縣和埇橋區(qū),常年冬小麥種植面積在40萬(wàn)hm2左右,冬小麥一般在秋季10月中、下旬播種,翌年6月初成熟。冬小麥生育期間多年平均降水量為300 mm左右、≥0 ℃積溫為2 300 ℃·d、日照時(shí)數(shù)約1 300 h,基本能滿足冬小麥生長(zhǎng)需求。

    1995—2013年冬小麥生育狀況等資料,取自宿州市農(nóng)業(yè)氣象試驗(yàn)站觀測(cè)地段觀測(cè)資料。觀測(cè)地段位于宿州市中部埇橋區(qū)內(nèi),在當(dāng)?shù)貙儆谥械确柿λ?,種植的品種: 1995—2009年為‘皖麥 19’、2010—2013年為‘皖麥52’。觀測(cè)地段冬小麥生育狀況等觀測(cè)方法與標(biāo)準(zhǔn),按照《農(nóng)業(yè)氣象觀測(cè)規(guī)范》執(zhí)行[27]。以1995—2013年度冬小麥各發(fā)育期的算術(shù)平均日期作為宿州市各縣(區(qū))冬小麥生育期常年平均日期(表1)。

    表1 1995—2013年宿州市冬小麥主要生育期平均日期Table 1 Average date of primary growth stages of winter wheat from 1955 to 2013 in Suzhou City

    冬小麥生育期氣候資料為1954—2013年宿州市各縣(區(qū))冬小麥生長(zhǎng)期間逐日日照、氣溫和降水量測(cè)量值,取自各縣(區(qū))國(guó)家氣象觀測(cè)站。統(tǒng)計(jì) 1954—2013年冬小麥各生育期逐日氣溫的平均值、日照時(shí)數(shù)和降水量累計(jì)值,作為各年度冬小麥各生育期的平均溫度、日照和降水量。按照世界氣象組織(WMO)規(guī)定,取某氣象要素1981—2010年的平均值作為冬小麥各生育期某要素的氣候平均值。

    1954—2013年度宿州市及泗縣、靈璧、蕭縣、碭山縣和埇橋區(qū)冬小麥逐年單產(chǎn)資料,取自宿州市統(tǒng)計(jì)局整編統(tǒng)計(jì)年鑒。

    1.2 研究方法

    1.2.1 云模型理論原理

    云模型由概率論和模糊數(shù)學(xué)演化發(fā)展而來(lái),能較好刻畫事件發(fā)生的模糊性和隨機(jī)性[20-21]。其定義為: 設(shè)Ω是一個(gè)精確數(shù)值表示的定量論域,C是Ω上的定性概念,即一個(gè)描述性的語(yǔ)言值或指標(biāo),對(duì)于任意一個(gè)論域中的元素 x,都存在一個(gè)有穩(wěn)定傾向的隨機(jī)數(shù)μ(x)∈[0,1],稱之為x對(duì)C的隸屬度,則x在論域Ω上的分布稱為云模型(cloud model),簡(jiǎn)稱為云(cloud),每個(gè)[x,μ(x)]稱為一個(gè)云滴[20]。

    云模型由期望(Ex)、熵(En)和超熵(He)3個(gè)數(shù)字特征或參數(shù)來(lái)表征[20-21]。其中: Ex標(biāo)定了云對(duì)象在論域中的位置,即云的重心位置,它完全屬于該集合。En是概念模糊度的度量,其大小直接決定了在論域中可被某一概念所接受的元素?cái)?shù),即亦此亦彼性的裕度。He也稱為熵的熵,是En的不確定性度量。

    正態(tài)云模型是在正態(tài)分布和鐘形隸屬函數(shù)的基礎(chǔ)上發(fā)展起來(lái)。李德毅等[20-21]已證明其具有一定普適性。以正態(tài)云模型為基礎(chǔ)可以進(jìn)行各種云運(yùn)算,獲得到左、右半云、組合云等等。當(dāng)期望值(Ex)是一個(gè)區(qū)間,即 a<xc時(shí),即論域中的元素值x不大于c或不小于c時(shí),μ(x)≡1,正態(tài)云即轉(zhuǎn)化為右半云或左半云[20]。

    同時(shí),李德毅等[20-21]還給出了正態(tài)云的生成算法,即:

    ③令[xi,μi(x)]為云滴,其隸屬度:

    ④重復(fù)步驟①~③,直到產(chǎn)生N個(gè)云滴為止,即生成云圖。

    1.2.2 基于云模型的冬小麥氣候適宜度評(píng)價(jià)方法

    光、溫、水對(duì)冬小麥生長(zhǎng)發(fā)育的影響,具有一定的隨機(jī)性和模糊性[1-2]。隨機(jī)性表現(xiàn)同一氣候要素的相同量值對(duì)冬小麥不同生育時(shí)期,或不同年份某氣候要素同一量值對(duì)冬小麥同一生育期生長(zhǎng)發(fā)育的影響是不確定的。模糊性表現(xiàn)為各氣候要素對(duì)冬小麥生長(zhǎng)發(fā)育影響的優(yōu)劣程度或適宜程度是難以確定的量化指標(biāo),但是這種影響又是客觀存在的。因此,依據(jù)云模型理論將氣候條件對(duì)冬小麥生長(zhǎng)發(fā)育影響的云模型進(jìn)行定義:

    假設(shè)Ω是一地某氣候要素可能出現(xiàn)情形的論域,C是論域Ω內(nèi)某氣候要素影響冬小麥生長(zhǎng)發(fā)育的描述值,對(duì)于論域Ω中任一元素(即量值x)其對(duì)冬小麥生長(zhǎng)發(fā)育的影響程度,都存在一個(gè)有穩(wěn)定度傾向的隨機(jī)數(shù)μ(x)∈[0,1],稱之為某氣候要素在某一數(shù)值x時(shí)對(duì)冬小麥生長(zhǎng)發(fā)育影響程度 C的隸屬度,那么,該氣候要素 x在其所有可能取值的論域Ω上的分布稱之為云模型,簡(jiǎn)稱云,每個(gè)[x,μ(x)]稱之為一個(gè)云滴,云滴的分布具有隨機(jī)性和模糊性。

    因此,基于云模型理論的冬小麥氣候適宜度評(píng)價(jià)方法,可歸結(jié)為光、溫、水各要素影響冬小麥生長(zhǎng)發(fā)育的云模型構(gòu)建。而氣候要素云模型構(gòu)建,主要是依據(jù)氣候要素影響冬小麥生長(zhǎng)發(fā)育的一些離散化的關(guān)鍵指標(biāo),如最適宜溫度、最高溫度、最低溫度、最適宜需水量等等,運(yùn)用云模型參數(shù)構(gòu)建方法,確定云模型的參數(shù)得到云模型。然后,再由氣候要素的實(shí)測(cè)值或預(yù)報(bào)值驅(qū)動(dòng)云模型,得到氣候要素對(duì)冬小麥生長(zhǎng)發(fā)育影響的隸屬度。依據(jù)隸屬度的大小,實(shí)現(xiàn)氣候條件對(duì)冬小麥生長(zhǎng)發(fā)育影響優(yōu)劣程度的量化評(píng)價(jià)。

    對(duì)于具有無(wú)確定度信息(即沒(méi)有條件約束)的數(shù)據(jù),如某一時(shí)段年降水量、年蒸散量等數(shù)據(jù),可以采用逆向正態(tài)云算法[20-21]求取云模型的 3個(gè)數(shù)字特征。對(duì)于有確定度(即有約束性)的問(wèn)題,如當(dāng)日照時(shí)數(shù)高于某一界限指標(biāo)時(shí),或氣溫、降水量在一定范圍內(nèi)時(shí),對(duì)冬小麥生長(zhǎng)發(fā)育的適宜程度最大等,其云模型的3個(gè)參數(shù),一般地,Ex 即取其約束性的描述值或界限值[21,26]。En可按照“3En”法則[20-21],取Ex距離上限(a)或下限(b)值的1/3來(lái)表示[20,28],即當(dāng)影響關(guān)系為左半云時(shí),En=(Ex-a)/3; 為右半云時(shí),En=(b-Ex)/3。而He,一般根據(jù)量值的波動(dòng)情況,并結(jié)合經(jīng)驗(yàn)確定[21,26]。

    1.2.3 冬小麥不同生育時(shí)期氣候適宜度

    冬小麥生長(zhǎng)發(fā)育過(guò)程中,某一生育時(shí)期的光照適宜度[μ(S)、溫度適宜度μ(T)、降水適宜度μ(R)相互作用決定一個(gè)生育時(shí)期的氣候適宜度[2,29]。即:

    式中: μj(S,T,R)為冬小麥第j個(gè)生育時(shí)期的氣候適宜度,μj(S)、μj(T)、μj(R)分別為第j個(gè)生育時(shí)期的日照適宜度、溫度適宜度、降水適宜度。

    1.2.4 冬小麥全生育期氣候適宜度

    冬小麥各個(gè)生育時(shí)期的生理生態(tài)特征不同,對(duì)氣候條件的需求也不同。同時(shí),各生育時(shí)期氣候因子對(duì)冬小麥生長(zhǎng)發(fā)育及產(chǎn)量形成的影響程度亦有差異。因此,為客觀反映不同時(shí)期氣候因子對(duì)冬小麥生長(zhǎng)發(fā)育的影響強(qiáng)度,計(jì)算全生育期氣候適宜度,必須為各生育期氣候適宜度設(shè)定權(quán)重。本文采用加權(quán)法[19]確定冬小麥全生育的氣候適宜度。權(quán)重系數(shù)采用一元積分回歸法確定,即分別計(jì)算逐生育期日照時(shí)數(shù)、氣溫和降水量對(duì)冬小麥產(chǎn)量的積分效應(yīng)系數(shù)(asj、atj、arj),每個(gè)生育期的積分效應(yīng)系數(shù)的絕對(duì)值除以所有生育期的積分效應(yīng)系數(shù)的絕對(duì)值的總和,作為每個(gè)生育期日照時(shí)數(shù)、氣溫和降水量適宜度的權(quán)重系數(shù)。以各要素逐生育期的隸屬度乘以對(duì)應(yīng)的權(quán)重系數(shù),然后相加,得到冬小麥全生育期日照時(shí)數(shù)、氣溫和降水量的適宜度。即:

    式中: bsj、btj、brj分別為第j個(gè)生育期日照時(shí)數(shù)、氣溫和降水量隸屬度的權(quán)重系數(shù),μ(S)、μ(T)、μ(R)分別為冬小麥全生育期日照時(shí)數(shù)、氣溫和降水量適宜度。

    仿照式(2),采用幾何平均法,得到冬小麥全生育氣候適宜度。即

    1.3 模型檢驗(yàn)與應(yīng)用

    為檢驗(yàn)本文建立的冬小麥氣候適宜度云模型能否客觀反映氣候條件對(duì)作物生長(zhǎng)發(fā)育的影響,將1954—2013年安徽省宿州市埇橋區(qū)及泗縣、靈璧、蕭縣和碭山縣冬小麥各年度的單產(chǎn)資料,先分解為趨勢(shì)產(chǎn)量和氣候產(chǎn)量,前者的變化主要是由社會(huì)經(jīng)濟(jì)因子造成,后者則主要取決于氣候因子。本文采用正交多項(xiàng)式,結(jié)合3年滑動(dòng)平均法,擬合得到趨勢(shì)產(chǎn)量,與實(shí)際產(chǎn)量對(duì)比,得到1954—2013年的氣候產(chǎn)量。經(jīng)檢驗(yàn),各地氣候產(chǎn)量符合正態(tài)分布。再利用1954—2013年宿州市埇橋區(qū)和泗縣、靈璧、蕭縣和碭山縣國(guó)家氣象觀測(cè)站氣象要素資料,計(jì)算各縣(區(qū))逐年冬小麥氣候適宜度,與對(duì)應(yīng)的冬小麥氣候產(chǎn)量作相關(guān)分析。

    同時(shí),利用1995—2013年宿州市農(nóng)業(yè)氣象試驗(yàn)站冬小麥觀測(cè)地段的逐年單產(chǎn)(kg·hm-2)、千粒重(g)、每穗籽粒重(g)、單位面積穗數(shù)(×104個(gè)·hm-2)和乳熟期植株高度(cm)等資料,分析與其氣候適宜度的關(guān)系,檢驗(yàn)氣候適宜度對(duì)冬小麥關(guān)鍵生長(zhǎng)要素的解釋能力。單產(chǎn)和部分產(chǎn)量因素和植株高度測(cè)定方法,均按照文獻(xiàn)[27]執(zhí)行。

    2 結(jié)果與分析

    2.1 日照適宜度云模型

    日照條件由少到多的變化,對(duì)冬小麥生長(zhǎng)發(fā)育的影響程度表現(xiàn)為由“不適宜”到“適宜”,在“適宜”程度時(shí)的日照時(shí)數(shù)稱為臨界值。由于在一地的長(zhǎng)期栽培、選育,促使農(nóng)作物會(huì)對(duì)日照條件具有很強(qiáng)的生態(tài)適應(yīng)性[30-31],而一地日照時(shí)數(shù)的氣候平均值,反映一定時(shí)期日照條件的平均狀態(tài)。因此,本文取冬小麥各生育期間日照時(shí)數(shù)的氣候平均值,作為日照時(shí)數(shù)適宜度最大時(shí)的臨界值。表2為宿州市埇橋區(qū)冬小麥不同生育期間日照時(shí)數(shù)的氣候平均值。

    一般地,當(dāng)日照時(shí)數(shù)超過(guò)臨界值時(shí),日照適宜度[μ(x)]最大,即μ(x)≡1。低于臨界值時(shí),光照條件不足,會(huì)對(duì)冬小麥的正常生長(zhǎng)發(fā)育產(chǎn)生影響,日照適宜度隨之降低,且偏少氣候平均值越多,日照適宜度越小,直至衰減為μ(x)=0。同時(shí),其在衰減過(guò)程中,其變化具有明顯的隨機(jī)性和模糊性。因此日照適宜度選用左半云模型來(lái)表達(dá)。

    表2 1981—2010年宿州市埇橋區(qū)冬小麥不同生育期日照時(shí)數(shù)氣候平均值Table 2 Climatic average of sunshine hours of winter wheat at different growth stages from 1981 to 2010 in Yongqiao District of Suzhou City

    根據(jù)1.2節(jié)方法,得到埇橋區(qū)冬小麥不同生育期日照時(shí)數(shù)適宜度的云模型參數(shù)(表3)。

    表3 宿州市埇橋區(qū)冬小麥不同生育期日照適宜度云模型參數(shù)Table 3 Cloud model parameters for the sunshine hours suitability of winter wheat at different growth stages in Yongqiao District of Suzhou City

    利用MATLAB語(yǔ)言,根據(jù)表3中各云模型參數(shù),采用逆向云生成法[20],取 500個(gè)云滴,通過(guò)計(jì)算機(jī)仿真,生成埇橋區(qū)冬小麥各生育期日照時(shí)數(shù)適宜度[μj(S)]云模型(圖1),j=1,2,…,7,即表1中的7個(gè)生育時(shí)期。

    由圖1和表3中Ex可以看出,由于不同生育期日照時(shí)數(shù)臨界值不同,其云模型的Ex存在一定差異,但在日照時(shí)數(shù)少于Ex時(shí),其對(duì)冬小麥生育影響的隸屬度(即云滴)呈霧化狀態(tài),較好地刻畫了日照條件低于臨界值時(shí),對(duì)冬小麥生育影響的不確定性,而非傳統(tǒng)上建立的日照隸屬度函數(shù),只是一條平滑的曲線[10-13],僅僅反映的是確定關(guān)系。

    圖1 宿州市埇橋區(qū)冬小麥不同生育期日照適宜度云模型Fig. 1 Cloud models for the sunshine hours suitability of winter wheat at different growth stages in Yongqiao District of Suzhou Citya: 播種—三葉期; b: 三葉期—越冬開(kāi)始; c: 越冬期; d: 返青—拔節(jié)期; e: 拔節(jié)—抽穗期; f: 抽穗—乳熟期; g: 乳熟—成熟期。a:seeding to trileaf stage; b: trileaf stage to wintering start; c: wintering stage; d: reviving to jointing stage; e: jointing to heading stage; f:heading to milky ripe stage; g: milky ripe to mature stage.

    利用式(1),當(dāng)已知某生育時(shí)期的日照時(shí)數(shù)(x),可以仿真得到日照時(shí)數(shù)為某一量值時(shí)的一組云滴[x,μ(x)],即同一日照時(shí)數(shù)在云模型中生成的云滴是隨機(jī)的,或者說(shuō)同一變量x可以隨機(jī)產(chǎn)生一組不同隸屬度μ(x)。

    本文在計(jì)算1954—2013年各縣(區(qū))冬小麥不同生育期日照時(shí)數(shù)所產(chǎn)生的云滴時(shí),均取 5個(gè)隸屬度的平均值,作為某一年不同生育期日照時(shí)數(shù)的適宜度,即:

    式中: i=1,2,…,5,即云滴[x,μ(x)]的個(gè)數(shù)。

    2.2 氣溫適宜度云模型

    冬小麥不同生育時(shí)期對(duì)溫度條件的要求不同,且存在最低、最高和最適宜溫度,其中最適宜溫度常常具有一定的幅度[10,31]。當(dāng)溫度處于最適宜范圍內(nèi)時(shí),冬小麥生長(zhǎng)良好,其溫度適宜度[μ(x)]最大,即μ(x)≡1。當(dāng)出現(xiàn)低于或高于冬小麥生長(zhǎng)發(fā)育的最適宜下限和上限溫度時(shí),即溫度處于最適宜溫度下限值到當(dāng)?shù)囟嗄昶骄畹蜌鉁刂?,或處于最適宜溫度上限到當(dāng)?shù)囟嗄昶骄罡邷囟戎g時(shí),溫度適宜度會(huì)減小,減小規(guī)律具有隨機(jī)性和模糊性。因此,溫度對(duì)冬小麥生長(zhǎng)發(fā)育的影響,具有梯形云特征[20]。本文采用梯形云來(lái)刻畫。綜合利用前人的研究成果[10-11,13],并結(jié)合宿州市多年農(nóng)業(yè)氣象服務(wù)經(jīng)驗(yàn),得到冬小麥各生育時(shí)期最低溫度、最高溫度和最適宜溫度指標(biāo)(表4)。

    表4 宿州市冬小麥不同生育期溫度指標(biāo)Table 4 Temperature indexes of winter wheat at different growth stages in Suzhou ℃

    根據(jù)1.2節(jié)方法,得到埇橋區(qū)冬小麥不同生育時(shí)期溫度適宜度的云模型參數(shù)(表5)。

    表5 宿州市埇橋區(qū)冬小麥不同生育時(shí)期溫度適宜度云模型參數(shù)Table 5 Cloud model parameters for the temperature suitability of winter wheat at different growth stages in Yongqiao District ofSuzhou City

    利用MATLAB語(yǔ)言,根據(jù)表5中各云模型參數(shù),采用逆向云生成法,取500個(gè)云滴,通過(guò)計(jì)算機(jī)仿真,生成埇橋區(qū)冬小麥各生育期氣溫適宜度[μj(T)]云模型(圖2)。

    由圖 2可以看出,不同生育期溫度處于適宜范圍內(nèi)時(shí),其隸屬度恒等于1。當(dāng)出現(xiàn)低于或高于適宜的界限溫度時(shí),溫度隸屬度呈霧化狀衰減,即溫度的適宜度越來(lái)越小,直至降低為0。這與溫度對(duì)冬小麥生育影響的生物學(xué)規(guī)律相一致。

    同理,利用式(1),取 1954—2013年各縣(區(qū))某一年某生育期氣溫?cái)?shù)據(jù),可以得到 5個(gè)云滴時(shí)的氣溫隸屬度,并由其均值作為不同年份不同生育期內(nèi)的氣溫適宜度[μj(T)]。

    2.3 降水適宜度云模型

    不同作物或者同一作物不同品種的不同生育期,對(duì)水分的需求量存在差異。本文按照侯英雨等[2]方法,即冬小麥某一生育時(shí)期(或階段)的需水量(ETp)是指在最適宜的土壤水分和肥力條件下,且田間正常發(fā)育、無(wú)病害并達(dá)到高產(chǎn)水平的特定作物的農(nóng)田蒸散量。其計(jì)算公式如下:

    ET0為計(jì)算時(shí)段內(nèi)參考作物蒸散量,采用國(guó)際糧農(nóng)組織(FAO)推薦的 Penman-Monteith模型(1998版)計(jì)算[32-33]。其中,Kc為相應(yīng)時(shí)段的作物系數(shù),本文綜合文獻(xiàn)[34-37]獲得(表6)。

    圖2 宿州市埇橋區(qū)冬小麥不同生育期溫度適宜度云模型Fig. 2 Cloud models for temperature suitability of winter wheat at different growth stages in Yongqiao District of Suzhou Citya: 播種—三葉期; b: 三葉期—越冬開(kāi)始; c: 越冬期; d: 返青—拔節(jié)期; e: 拔節(jié)—抽穗期; f: 抽穗—乳熟期; g: 乳熟—成熟期。a:seeding to trileaf stage; b: trileaf to wintering start; c: wintering stage; d: reviving to jointing stage; e: jointing to heading stage; f: heading to milky ripe stage; g: milky ripe to mature stage.

    表6 宿州市冬小麥各月作物系數(shù)(Kc)Table 6 Crop coefficients of winter wheat in each month in Suzhou City(Kc)

    淮北平原屬雨養(yǎng)冬麥區(qū),大氣降水是冬小麥生長(zhǎng)發(fā)育所需水分的主要來(lái)源[6,35,37]。一般地,降水量過(guò)多或過(guò)少,對(duì)冬小麥的生長(zhǎng)發(fā)育均不利,即冬小麥生長(zhǎng)發(fā)育對(duì)降水量的需求有一個(gè)適宜范圍。在適宜降水量范圍之內(nèi),冬小麥生長(zhǎng)發(fā)育良好,氣候適宜度可達(dá)到最大值。當(dāng)降水量低于或超過(guò)適宜降水量下限或上限時(shí),降水適宜度會(huì)隨之減小。因此,降水量對(duì)冬小麥生長(zhǎng)發(fā)育的影響與溫度條件的影響類似,可以利用梯形云來(lái)刻畫。

    利用文獻(xiàn)[19]的成果,結(jié)合淮北平原氣候生態(tài)特點(diǎn),在拔節(jié)開(kāi)始以前和進(jìn)入拔節(jié)期以后,確定當(dāng)降水量(R)滿足 0.65ETp≤R≤1.20ETp時(shí),氣候適宜度[(x)]最大,即μ(x)≡1。返青—拔節(jié)期間,當(dāng)降水量滿足0.75ETp≤R≤1.50 ETp時(shí),氣候適宜度[μ(x)]最大,即μ(x)≡1。當(dāng)R<0.65ETp或<0.75ETp或R>1.20ETp或>1.50ETp時(shí),降水量對(duì)冬小麥生長(zhǎng)發(fā)育均不利,氣候適宜度會(huì)隨之減小,直至降至0。而R=0.65ETp或0.75ETp和R=1.20ETp或1.50ETp,分別為相應(yīng)生育時(shí)期適宜降水量的下限和上限指標(biāo)。降水量為0和1960—2013年間降水量的最大值,分別作為降水量發(fā)生可能量值的限定指標(biāo)(表7)。

    表7 宿州市冬小麥不同生育期水分指標(biāo)Table 7 Water indexes of winter wheat at different growth stages in Suzhou City mm

    依據(jù)1.2節(jié)方法,得到埇橋區(qū)冬小麥不同生育時(shí)期降水適宜度的云模型參數(shù)(表8)。

    利用MATLAB語(yǔ)言,根據(jù)表8中各云模型參數(shù),采用逆向云生成法,取 500個(gè)云滴,通過(guò)計(jì)算機(jī)仿真,生成埇橋區(qū)各生育時(shí)期降水量適宜度[μj(R)]云模型(圖3)。

    表8 宿州市埇橋區(qū)冬小麥不同生育期降水適宜度云模型參數(shù)Table 8 Cloud model parameters for the precipitation suitability of winter wheat at different growth stages in Yongqiao District of Suzhou City

    圖3 宿州埇橋區(qū)冬小麥不同生育期降水適宜度云模型Fig. 3 Cloud models for precipitation suitability of winter wheat at different growth stages in Yongqiao District of Suzhou Citya: 播種—三葉期; b: 三葉期—越冬開(kāi)始; c: 越冬期; d: 返青—拔節(jié)期; e: 拔節(jié)—抽穗期; f: 抽穗—乳熟期; g: 乳熟—成熟期。a:seeding to trileaf stage; b: trileaf to wintering start; c: wintering stage; d: reviving to jointing stage; e: jointing to heading stage; f: heading to milky ripe stage; g: milky ripe to mature stage.

    由圖 3可以看出,不同生育期降水量在適宜范圍內(nèi)時(shí),其隸屬度恒等于1。當(dāng)出現(xiàn)低于或高于適宜的界限降水量時(shí),降水隸屬度呈霧化狀衰減,即降水的適宜度越來(lái)越小,直至降低為0。這與降水對(duì)冬小麥生育影響的生物學(xué)規(guī)律一致。

    同時(shí),由圖 3可以發(fā)現(xiàn),播種—三葉期和乳熟—成熟期2個(gè)云圖,左半云較右半云衰減趨勢(shì)明顯,即當(dāng)降水量低于適宜值下限指標(biāo)時(shí),這 2個(gè)時(shí)期降水適宜度的降低速度,明顯快于降水量多于適宜值上限時(shí)的降低速度,表明降水偏少較偏多時(shí)的不利影響為大。這應(yīng)是宿州市小麥生產(chǎn)需要注意的。

    同理,利用式(1),取 1954—2013年間各縣(區(qū))某一年某生育期內(nèi)的降水量數(shù)據(jù),可以得到 5個(gè)云滴的降水隸屬度,并由其均值作為某一年不同生育期的降水適宜度μj(R)。

    2.4 冬小麥氣候適宜度變化特征

    2.4.1 冬小麥全生育期內(nèi)光、溫、水適宜度變化

    將1954—2013年宿州市埇橋區(qū)、泗縣、靈璧、蕭縣和碭山縣各年冬小麥氣候適宜度進(jìn)行平均,得到宿州市各年氣候適宜度(圖4)。

    由圖4可以看出,1954—2013年間,宿州市冬小麥日照適宜度和溫度適宜度最高,平均值在 0.9以上; 降水適宜度最小,平均值在0.6以下,且降水適宜度的年際間波動(dòng)幅度明顯高于日照和溫度適宜度,說(shuō)明宿州市光、熱資源一般能夠滿足冬小麥生長(zhǎng)發(fā)育需求,但降水條件是冬小麥產(chǎn)量形成的限制因素。經(jīng)相關(guān)統(tǒng)計(jì)分析,1954—2013年間,日照適宜度和降水適宜度呈下降趨勢(shì),10年下降速率分別為0.005和0.008。氣溫適宜度呈上升趨勢(shì),10年上升速率為0.028,且通過(guò)P<0.01的F檢驗(yàn),表明氣溫對(duì)宿州冬小麥生長(zhǎng)發(fā)育影響的正效應(yīng)在增加,降水與日照時(shí)數(shù)的負(fù)效應(yīng)在增大。

    2.4.2 冬小麥全生育期氣候適宜度變化特征

    圖5是宿州市冬小麥全生育期氣候適宜度的年際變化及其變化趨勢(shì)。

    由圖5可以看出,1954—2013年間,宿州市冬小麥氣候適宜度在 0.7左右,但總體變化呈增加趨勢(shì),10年上升速率為0.002。從年際間變化看,1954—1960年為氣候適宜度波動(dòng)較大階段。自 1960年開(kāi)始到1974年,氣候適宜度變化平穩(wěn),且多為0.7。其后,自1975年以來(lái),氣候適宜度又進(jìn)入波動(dòng)較大階段,除1990—1994年為一個(gè)連續(xù)在0.7以上的時(shí)段外,其他年際間的氣候適宜度基本上是在0.7的均值線附近作上下波動(dòng)。但進(jìn)入 2003年后,氣候適宜度總體上呈降低趨勢(shì)。由是可見(jiàn),在全球氣候變化影響下,宿州氣候?qū)Χ←溕L(zhǎng)發(fā)育產(chǎn)生的負(fù)效應(yīng)正在增大,這在一定程度上增加了本區(qū)域冬小麥生產(chǎn)的氣候風(fēng)險(xiǎn)。

    圖4 宿州市冬小麥全生育期光、溫、水適宜度年際變化(1954—2013年)Fig. 4 Inter-annual variations of climatic suiabilities of annual sunshine hours,temperature and precipitation during winter wheat growth period in Suzhou City from 1954 to 2013

    圖5 宿州市冬小麥全生育期氣候適宜度年際變化(1954—2013年)Fig. 5 Inter-annual variation of climatic suitability of winter wheat during whole growth period in Suzhou City from 1954 to 2013

    2.5 氣候適宜度與氣候產(chǎn)量的關(guān)系

    根據(jù)本文確定的方法和云模型參數(shù),計(jì)算得到宿州市埇橋區(qū)及泗縣、靈璧、蕭縣和碭山縣冬小麥生育期間的氣候適宜度,并計(jì)算與各自氣候產(chǎn)量的相關(guān)系數(shù)(樣本數(shù)量均為60個(gè))。

    結(jié)果表明,埇橋區(qū)、泗縣、靈璧、蕭縣和碭山縣氣候適宜度與氣候產(chǎn)量之間呈顯著正相關(guān)關(guān)系,相關(guān)系數(shù)分別為0.314 1(P<0.05)、0.332 7(P<0.01)、0.321 1(P<0.05)、0.273 2(P<0.05)和0.332 7(P<0.01),其中灌溉條件相對(duì)較弱的泗縣和碭山縣相關(guān)系數(shù)最大,通過(guò)P<0.01的顯著性檢驗(yàn)。說(shuō)明本文設(shè)計(jì)的冬小麥氣候適宜度模型能客觀反映宿州市冬小麥的氣候適宜性水平及其動(dòng)態(tài)變化。

    2.6 氣候適宜度與產(chǎn)量因素的關(guān)系

    利用 1995—2013年宿州市農(nóng)業(yè)氣象試驗(yàn)站冬小麥觀測(cè)地段的監(jiān)測(cè)數(shù)據(jù),選取地段單產(chǎn)、單位面積穗數(shù)、每穗籽粒數(shù)、每穗籽粒重、千粒重和乳熟期植株高度等,檢驗(yàn)其與氣候適宜度之間的相關(guān)關(guān)系。采取正交多項(xiàng)式方法,對(duì)地段1995—2013年的單產(chǎn)資料進(jìn)行分離,得到地段相應(yīng)年度的氣候產(chǎn)量。分別計(jì)算1995—2013年冬小麥全生育期氣候適宜度與地段氣候產(chǎn)量、千粒重、每穗籽粒數(shù)、乳熟期植株高度之間的相關(guān)關(guān)系(樣本數(shù)量均為19個(gè))。

    結(jié)果表明,冬小麥全生育期氣候適宜度與氣候產(chǎn)量、千粒重、每穗籽粒數(shù)、乳熟期植株高度之間呈顯著或極顯著正相關(guān)關(guān)系(圖6),相關(guān)系數(shù)分別為0.588 0(P<0.01)、0.756 1(P<0.01)、0.670 7(P<0.01)、0.464 3(P<0.05)??梢?jiàn),冬小麥全生育期氣候適宜度對(duì)冬小麥生長(zhǎng)發(fā)育狀況和產(chǎn)量高低具有很強(qiáng)的解釋能力。即氣候適宜度越大,產(chǎn)量水平越高、千粒重越大、穗粒數(shù)越多。

    同時(shí),全生育期氣候適宜度與各參量的回歸方程(圖6),經(jīng)F檢驗(yàn),各方程均通過(guò)P<0.01及以上水平顯著性檢驗(yàn),即可以由不同生育期的氣候適宜度,反演得到相應(yīng)的產(chǎn)量因素值,用于生長(zhǎng)狀況評(píng)估。

    圖6 宿州市冬小麥全生育期氣候適宜度與氣候產(chǎn)量、千粒重、穗粒數(shù)和株高的相關(guān)關(guān)系(1995—2013年)Fig. 6 Correlation between the climatic suitability of winter wheat during whole growth period and climatic yield,1000-grain weight,grain number per spike and plant height in Suzhou City from 1995 to 2013

    另外,考慮到不同產(chǎn)量構(gòu)成因素形成的時(shí)期不同,對(duì)氣候條件的需求也不同。返青—拔節(jié)期、拔節(jié)—抽穗期、抽穗—乳熟期分別是冬小麥穗數(shù)、單穗籽粒重形成的關(guān)鍵期[38],因此,分別檢驗(yàn) 3個(gè)時(shí)期的氣候適宜度與單位面積穗數(shù)、單株成穗數(shù)和每穗籽粒重的相關(guān)關(guān)系。結(jié)果表明,各相關(guān)關(guān)系分別達(dá)到顯著或極顯著水平,其中: 返青—拔節(jié)期氣候適宜度與單位面積穗數(shù)、每穗籽粒數(shù)的相關(guān)系數(shù)分別為0.558 9(P<0.05)、0.710 7(P<0.01),拔節(jié)—抽穗期氣候適宜度與單位面積穗數(shù)的相關(guān)系數(shù)為0.649 8(P<0.01),抽穗—乳熟期氣候適宜度與單株成穗數(shù)、每穗籽粒重的相關(guān)系數(shù)分別為0.736 1(P<0.01)、0.744 2(P<0.01)??梢?jiàn),隨著生育時(shí)期的細(xì)分,氣候適宜度與不同生育期所需達(dá)到的產(chǎn)量因素目標(biāo)之間的正相關(guān)性在增大,進(jìn)一步說(shuō)明,本文建立的氣候適宜度對(duì)不同氣候條件影響冬小麥生長(zhǎng)發(fā)育的程度具有很強(qiáng)的解釋能力。

    3 結(jié)論

    1)本文基于云模型理論,建立了日照時(shí)數(shù)適宜度評(píng)價(jià)的左半云模型、溫度與降水適宜度評(píng)價(jià)的梯形云模型。模型生物學(xué)意義明確,且兼顧了氣候條件對(duì)農(nóng)作物生育影響優(yōu)劣程度的隨機(jī)性和模糊性。經(jīng)統(tǒng)計(jì)檢驗(yàn),計(jì)算得到的氣候適宜度,與冬小麥單產(chǎn)及其部分產(chǎn)量構(gòu)成因素之間具有顯著的正相關(guān)關(guān)系,分別通過(guò)顯著(P<0.05)或極顯著(P<0.01)水平檢驗(yàn); 不同關(guān)鍵生育期的氣候適宜度與其培育目標(biāo)之間呈顯著正相關(guān)關(guān)系,尤以抽穗—乳熟期氣候適宜度與單株成穗數(shù)、每穗籽粒重的相關(guān)系數(shù)最高,分別為0.736 1(P<0.01)和0.744 2(P<0.01)。因此,所建模型可以用于解釋氣候條件對(duì)冬小麥生育影響的優(yōu)劣程度。

    2)1954—2013年宿州市冬小麥生育期間各要素氣候適宜度存在明顯年際變化,尤以降水適宜度年際間波動(dòng)幅度明顯,且多年平均值在0.6以下,而日照與溫度適宜度在0.9以上。60年間日照與降水適宜度呈下降趨勢(shì),每 10年的下降速率分別為 0.005 和 0.008。氣溫適宜度呈增加趨勢(shì),每 10年的上升速率為 0.028。但進(jìn)入 2003年以來(lái),綜合氣候適宜度呈下降趨勢(shì),氣候條件對(duì)宿州冬小麥生育產(chǎn)生的負(fù)效應(yīng)在增大。這與同為黃海冬麥區(qū)的河南省、河北省冬小麥氣候適宜度變化特征基本一致[11,19,39]。

    3)本文構(gòu)建氣候適宜度評(píng)價(jià)云模型時(shí),涉及的日照時(shí)數(shù)、溫度和降水量等指標(biāo),多是從已有文獻(xiàn)并結(jié)合宿州市多年氣象服務(wù)經(jīng)驗(yàn)獲取的。其中,日照時(shí)數(shù)臨界指標(biāo)是依據(jù)作物在一個(gè)地區(qū)長(zhǎng)期栽培、選育后對(duì)其環(huán)境具有較強(qiáng)適應(yīng)性[30-31]的生態(tài)學(xué)原理確定的,而不同于一些文獻(xiàn)選取可照時(shí)數(shù)的 70%作為臨界值[10-11],能更好地反映作物對(duì)一地氣候條件的適應(yīng)性,這為本文創(chuàng)新之處。同時(shí),本文在確定降水量指標(biāo)時(shí),雖然利用權(quán)重系數(shù)間接考慮了不同生育期需水的動(dòng)態(tài)特性,但未能考慮降水的后效性和有效性[2]問(wèn)題,是所建模型的不足之處,是今后需要加以改進(jìn)之處。

    4)對(duì)于定性概念值轉(zhuǎn)換成云參數(shù),目前尚無(wú)完整的理論來(lái)證明何種方法確定En較為科學(xué)[20]。根據(jù)研究問(wèn)題實(shí)際,本文結(jié)合經(jīng)驗(yàn)和云模型的“3En”法則[20-21]確定云參數(shù)。如何確定更客觀的 En、He將是今后工作的重點(diǎn)。另外,本文構(gòu)建的氣候適宜度評(píng)價(jià)云模型,是依據(jù)淮北平原宿州氣候特點(diǎn)和相關(guān)指標(biāo)來(lái)確定的,在參考應(yīng)用時(shí),應(yīng)結(jié)合各地特點(diǎn)分別建立適用于區(qū)域特點(diǎn)的云模型。

    References

    [1] 羅懷良,陳國(guó)階,朱波. 農(nóng)業(yè)生態(tài)氣候適宜度研究進(jìn)展[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃,2004,25(1): 28-32 Luo H L,Chen G J,Zhu B. Review on suitability of agro-climate[J]. Journal of China Agricultural Resources and Regional Planning,2004,25(1): 28-32

    [2] 侯英雨,張艷紅,王良宇,等. 東北地區(qū)春玉米氣候適宜度模型[J]. 應(yīng)用生態(tài)學(xué)報(bào),2013,24(11): 3207-3212 Hou Y Y,Zhang Y H,Wang L Y,et al. Climatic suitability model for spring maize in Northeast China[J]. Chinese Journal of Applied Ecology,2013,24(11): 3207-3212

    [3] Tao F L,Yokozawa M,Xu Y L,et al. Climate changes and trends in phenology and yields of field crops in China,1981-2000[J]. Agricultural and Forest Meteorology,2006,138(1/4): 82-92

    [4] 劉志娟,楊曉光,王文峰,等. 氣候變化背景下我國(guó)東北三省農(nóng)業(yè)氣候資源變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào),2009,20(9):2199-2206 Liu Z J,Yang X G,Wang W F,et al. Characteristics of agricultural climate resources in three provinces of Northeast China under global climate change[J]. Chinese Journal of Applied Ecology,2009,20(9): 2199-2206

    [5] 謝立勇,李悅,徐玉秀,等. 氣候變化對(duì)農(nóng)業(yè)生產(chǎn)與糧食安全影響的新認(rèn)知[J]. 氣候變化研究進(jìn)展,2014,10(4):235-239 Xie L Y,Li Y,Xu Y X,et al. Updated understanding on the impacts of climate change on food production and food security[J]. Progressus Inquisitiones de Mutatione Climatis,2014,10(4): 235-239

    [6] 馬曉群,張宏群,吳文玉,等. 安徽省冬小麥品種生態(tài)氣候適宜性分析和精細(xì)化區(qū)劃[J]. 中國(guó)農(nóng)業(yè)氣象,2012,33(1):86-92 Ma X Q,Zhang H Q,Wu W Y,et al. Analyzing and zoning of the eco-climate suitability on winter wheat varieties in Anhui Province[J]. Chinese Journal of Agrometeorology,2012,33(1): 86-92

    [7] 李德,楊太明,劉瑞娜,等. 安徽省設(shè)施農(nóng)業(yè)冬季低溫風(fēng)險(xiǎn)分析和區(qū)劃[J]. 中國(guó)農(nóng)業(yè)氣象,2013,34(6): 703-709 Li D,Yang T M,Liu R N,et al. Low temperature risk division in winter for protected agriculture in Anhui Province[J]. Chinese Journal of Agrometeorology,2013,34(6): 703-709

    [8] 顧恒岳,艾南山. 農(nóng)業(yè)氣候系統(tǒng)及其動(dòng)態(tài)模型[J]. 大自然探索,1984(1): 43-56 Gu H Y,Ai N S. The dynamic models of agro-climate system[J]. Discovery of Nature,1984(1): 43-56

    [9] 徐紅,李寶鴻. 山東省農(nóng)業(yè)生態(tài)氣候資源適宜度分析[J].山東師大學(xué)報(bào): 自然科學(xué)版,1993,8(4): 41-46 Xu H,Li B H. Analysis of agro-ecolmate resources in Shandong Province with suitability degree[J]. Journal of Shandong Normal University: Natural Sciences Edition,1993,8(4):41-46

    [10] 黃璜. 中國(guó)紅黃壤地區(qū)作物生產(chǎn)的氣候生態(tài)適應(yīng)性研究[J].自然資源學(xué)報(bào),1996,11(4): 340-346 Huang H. A study on the climatic ecology adaptability of the crop production in the red and yellow soils region of China[J]. Journal of Natural Resources,1996,11(4): 340-346

    [11] 趙峰,千懷遂,焦士興. 農(nóng)作物氣候適宜度模型研究——以河南省冬小麥為例[J]. 資源科學(xué),2003,25(6): 77-82 Zhao F,Qian H S,Jiao S X. Climatic suitability model of crop:A case study of winter wheat in Henan Province[J]. Resources Science,2003,25(6): 77-82

    [12] 任玉玉,千懷遂. 河南省棉花氣候適宜度變化趨勢(shì)分析[J].應(yīng)用氣象學(xué)報(bào),2006,17(1): 87-93 Ren Y Y,Qian H S. Climatic suitability of cotton and its changes in Henan Province[J]. Journal of Applied Meteorological Science,2006,17(1): 87-93

    [13] 魏瑞江,宋迎波,王鑫. 基于氣候適宜度的玉米產(chǎn)量動(dòng)態(tài)預(yù)報(bào)方法[J]. 應(yīng)用氣象學(xué)報(bào),2009,20(5): 622-627 Wei R J,Song Y B,Wang X. Method for dynamic forecast of corn yield based on climatic suitability[J]. Journal of Applied Meteorological Science,2009,20(5): 622-627

    [14] 易雪,王建林,宋迎波. 氣候適宜指數(shù)在早稻產(chǎn)量動(dòng)態(tài)預(yù)報(bào)上的應(yīng)用[J]. 氣象,2010,36(6): 85-89 Yi X,Wang J L,Song Y B. Application of climatic suitability index to dynamical prediction of early rice yield[J]. Meteorological Monthly,2010,36(6): 85-89

    [15] 李昊宇,王建林,鄭昌玲,等. 氣候適宜度在華北冬小麥發(fā)育期預(yù)報(bào)中的應(yīng)用[J]. 氣象,2012,38(12): 1554-1559 Li H Y,Wang J L,Zheng C L,et al. The development period prediction of winter wheat based on climatic suitability inNorth China[J]. Meteorological Monthly,2012,38(12):1554-1559

    [16] 馬興祥,鄧振鏞,李棟梁,等. 甘肅省春小麥生態(tài)氣候適宜度在適生種植區(qū)劃中的應(yīng)用[J]. 應(yīng)用氣象學(xué)報(bào),2005,16(6):820-827 Ma X X,Deng Z Y,Li D L,et al. Study on eco-climate applicability of spring wheat for condign planting division in Gansu Province[J]. Journal of Applied Meteorological Science,2005,16(6): 820-827

    [17] 姚小英,蒲金涌,姚茹莘,等. 氣候暖干化背景下甘肅旱作區(qū)玉米氣候適宜性變化[J]. 地理學(xué)報(bào),2011,66(1): 59-67 Yao X Y,Pu J Y,Yao R X,et al. Variation of climate suitability of maize in arid area in Gansu under the condition of climate dry-warming[J]. Acta Geographica Sinica,2011,66(1):59-67

    [18] 段海來(lái),千懷遂,李明霞,等. 中國(guó)亞熱帶地區(qū)柑桔的氣候適宜性[J]. 應(yīng)用生態(tài)學(xué)報(bào),2010,21(8): 1915-1925 Duan H L,Qian H S,Li M X,et al. Climatic suitability of citrus in subtropical China[J]. Chinese Journal of Applied Ecology,2010,21(8): 1915-1925

    [19] 魏瑞江,張文宗,康西言,等. 河北省冬小麥氣候適宜度動(dòng)態(tài)模型的建立及應(yīng)用[J]. 干旱地區(qū)農(nóng)業(yè)研究,2007,25(6):5-9 Wei R J,Zhang W Z,Kang X Y,et al. Application and establishment of climatic suitability dynamic model of winter wheat in Hebei Province[J]. Agricultural Research in the Arid Areas,2007,25(6): 5-9

    [20] 李德毅,杜鹢. 不確定性人工智能[M]. 北京: 國(guó)防工業(yè)出版社,2005: 131-165 Li D Y,Du Y. Artificial Intelligent with Uncertainty[M]. Beijing: National Defence Industry Press,2005: 131-165

    [21] 劉常昱,李德毅,潘莉莉. 基于云模型的不確定性知識(shí)表示[J]. 計(jì)算機(jī)工程與應(yīng)用,2004,40(2): 32-35 Liu C Y,Li D Y,Pan L L. Uncertain knowledge representation based on cloud model[J]. Computer Engineering and Applications,2004,40(2): 32-35

    [22] 李剛,萬(wàn)幼川. 基于高維云模型和 RBF神經(jīng)網(wǎng)絡(luò)的遙感影像不確定性分類方法[J]. 測(cè)繪科學(xué),2012,37(1): 115-118 Li G,Wan Y C. Uncertainty classification method of remote sensing image based on high-dimensional cloud model and RBF neural network[J]. Science of Surveying Mapping,2012,37(1): 115-118

    [23] 張楊,嚴(yán)金明,江平,等. 基于正態(tài)云模型的湖北省土地資源生態(tài)安全評(píng)價(jià)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(22): 252-258 Zhang Y,Yan J M,Jiang P,et al. Normal cloud model based evaluation of land resources ecological security in Hubei Province[J]. Transactions of the Chinese Society of Agricultural Engineering,2013,29(22): 252-258

    [24] 孫鴻鵠,程先富,倪玲,等. 基于云模型和熵權(quán)法的巢湖流域防洪減災(zāi)能力評(píng)估[J]. 災(zāi)害學(xué),2015,30(1): 222-227 Sun H H,Cheng X F,Ni L,et al. Capacity evaluation of flood disaster prevention and reduction in Chaohu Basin based on cloud model and entropy weight method[J]. Journal of Catastrophology,2015,30(1): 222-227

    [25] 張秋文,章永志,鐘鳴. 基于云模型的水庫(kù)誘發(fā)地震風(fēng)險(xiǎn)多級(jí)模糊綜合評(píng)價(jià)[J]. 水利學(xué)報(bào),2014,45(1): 87-95 Zhang Q W,Zhang Y Z,Zhong M. A cloud model based approach for multi-hierarchy fuzzy comprehensive evaluation of reservoir-induced seismic risk[J]. Journal of Hydraulic Engineering,2014,45(1): 87-95

    [26] 李德,楊霏云,孫義,等. 基于云模型的淮北平原麥田噴藥(肥)氣象適宜度評(píng)價(jià)[J]. 生態(tài)學(xué)雜志,2016,35(1): 259-268 Li D,Yang F Y,Sun Y,et al. Meteorological suitability evaluation based on cloud model for spraying pesticide/ fertilizer in wheat fields of Huaibei Plain[J]. Chinese Journal of Ecology,2016,35(1): 259-268

    [27] 國(guó)家氣象局. 農(nóng)業(yè)氣象觀測(cè)規(guī)范(上卷)[M]. 北京: 氣象出版社,1993: 10-31 China Meteorological Administration. Observation Specification of Agro-meteorology[M]. Beijing: China Meteorological Press,1993: 10-31

    [28] 張國(guó)英,沙云,劉旭紅,等. 高維云模型及其在多屬性評(píng)價(jià)中的應(yīng)用[J]. 北京理工大學(xué)學(xué)報(bào),2004,24(12): 1065-1069 Zhang G Y,Sha Y,Liu X H,et al. High dimensional cloud model and its application in multiple attribute evaluation[J]. Transactions of Beijing Institute of Technology,2004,24(12):1065-1069

    [29] 馬樹(shù)慶. 吉林省農(nóng)業(yè)氣候研究[M]. 北京: 氣象出版社,1994: 33 Ma S Q. Agro-climate Research in Jilin Province[M]. Beijing:China Meteorological Press,1994: 33

    [30] 韓湘玲. 作物生態(tài)學(xué)[M]. 北京: 氣象出版社,1991: 15-18 Han X L. Crop Ecology[M]. Beijing: China Meteorological Press,1991: 15-18

    [31] 龔紹先. 糧食作物與氣象[M]. 北京: 北京農(nóng)業(yè)大學(xué)出版社,1987: 35-45 Gong S X. Crop and Meteorology[M]. Beijing: Agriculture University Press,1987: 35-45

    [32] Allen R G,Smith M,Pereira L S,et al. An update for the calculation of reference evapotranspiration[J]. ICID Bulletin,1994,43(2): 35-92

    [33] 劉鈺. Pereira L S,Teixeira J L,等. 參照騰發(fā)量的新定義及計(jì)算方法對(duì)比[J]. 水利學(xué)報(bào),1997(6): 27-33 Liu Y,Pereira L S,Teixeira J L,et al. Update definition and computation of reference evapotranspiration comparison with former method[J]. Journal of Hydraulic Engineering,1997(6):27-33

    [34] 陳曉藝,馬曉群,孫秀邦. 安徽省冬小麥發(fā)育期農(nóng)業(yè)干旱發(fā)生風(fēng)險(xiǎn)分析[J]. 中國(guó)農(nóng)業(yè)氣象,2008,29(4): 472-476 Chen X Y,Ma X Q,Sun X B. Risk analysis of agricultural drought for winter wheat during growing period in Anhui Province[J]. Chinese Journal of Agrometeorology,2008,29(4): 472-476

    [35] 馬曉群,姚筠,許瑩. 安徽省農(nóng)作物干旱損失動(dòng)態(tài)評(píng)估模型及其試用[J]. 災(zāi)害學(xué),2010,25(1): 13-17 Ma X Q,Yao Y,Xu Y. A model for dynamic assessment ofcrop yield losses from drought and its tryout in Anhui Province[J]. Journal of Catastrophology,2010,25(1): 13-17

    [36] 丁德峻,張旭暉. 江蘇淮北地區(qū)主要作物需水量的初步研究[J]. 自然資源,1994(3): 40-46 Ding D J,Zhang X H. A preliminary study on water requirement of major crops in Huaibei area of Jiangsu Province[J]. Resources Science,1994(3): 40-46

    [37] 王曉東,馬曉群,許瑩,等. 淮河流域主要農(nóng)作物全生育期水分盈虧時(shí)空變化分析[J]. 資源科學(xué),2013,35(3): 665-672 Wang X D,Ma X Q,Xu Y,et al. Temporal analysis of the crop water surplus deficit index for the whole growth period in the Huaihe Basin[J]. Resources Science,2013,35(3):665-672

    [38] 胡承霖. 安徽麥作學(xué)[M]. 合肥: 安徽科學(xué)技術(shù)出版社,2009: 1-14 Hu C L. Anhui Wheat Crop Cultivation[M]. Hefei: Anhui Science and Technology Press,2009: 1-14

    [39] 千懷遂,焦士興,趙峰. 河南省冬小麥氣候適宜性變化研究[J]. 生態(tài)學(xué)雜志,2005,24(5): 503-507 Qian H S,Jiao S X,Zhao F. Climate suitability change of winter wheat in Henan Province[J]. Chinese Journal of Ecology,2005,24(5): 503-507

    Evaluation of climate suitability of winter wheat based on cloud model analysis — A case study of Suzhou, Anhui Province*

    LI De1,ZHOU Wenlin2,YANG Feiyun3
    (1. Suzhou Meteorological Bureau in Anhui Province,Suzhou 234000,China; 2. Hefei Meteorological Bureau in Anhui Province,Hefei 236001,China; 3. China Meteorological Administration Training Center,Beijing 100081,China)

    Based on three limiting indexes(light,temperature and water),an evaluation system was set up along with its parameters for normal cloud models and the “3En” rule to evaluate the suitability of winter wheat to climatic factors. The weighted comprehensive assessment and geometric mean method were used to determine the climate suitability of winter wheat for the whole growth period,while the integral regression method was used to determine the weight coefficients for each growth stage. The calculated values were analyzed in relation to winter wheat yield in every district of Suzhou City of Anhui Province from 1954 to 2013 and the actual observed yield and its’ components from 1995 to 2013. The results showed that the trapezium cloud model was suitable for air temperature and precipitation suitability evaluation,while the left-half cloud model was suitable for sunshine suitability analysis. There was significantly positive correlation between the climate suitability of winter wheat for the whole growth period and climate-driven yield,1000-grain weight,kernel number per ear,and plant height at milk stage in the study area. The corresponding cor-relation coefficients were 0.588 0(P < 0.01),0.756 1(P < 0.01),0.670 7(P < 0.01) and 0.464 3(P < 0.05),respectively. The correlation coefficients between the climate suitability and panicle per unit area of winter wheat were 0.558 9(P < 0.05),0.649 8(P < 0.01)and 0.736 1(P < 0.01) at returning green-jointing stage,jointing-heading stage and heading-milk ripe stage,respectively. There was a significantly positive correlation between climate suitability and kernel number per ear — 0.710 7(P < 0.01) and 0.744 2(P < 0.01) at returning green-jointing stage and heading-milk ripe stage,respectively. The suitability of sunshine and precipitation decreased at the rates of 0.005 and 0.008 per decade. However,the suitability of temperature in the study area increased at the rate of 0.028 per decade in 1954-2013. The study laid the scientific reference for evaluating adaptability and developing response strategies of winter wheat to climatic conditions.

    Winter wheat; Cloud model; ‘3En’ rule; Climate suitability; Climatic yield; Observed yield; Climate factor;Growth period; Suzhou city

    Dec. 14,2015; accepted Feb. 19,2016

    S11+9; S162.5+3

    A

    1671-3990(2016)07-0944-13

    10.13930/j.cnki.cjea.151317

    * 安徽省氣象新技術(shù)集成項(xiàng)目(AHXJ201207)和國(guó)家公益性(氣象)行業(yè)專項(xiàng)項(xiàng)目(GYHY201006027)資助

    李德,主要從事應(yīng)用氣象業(yè)務(wù)與研究工作。E-mail: szlide@sohu.com

    2015-12-14 接受日期: 2016-02-19

    * This study was supported by the Project of New Meteorological Technology Integration in Anhui Province(AHXJ201207),and the Special Fund for Meteorological Scientific Research in the Public Interest(GYHY201006027).

    Corresponding author,LI De,E-mail: szlide@sohu.com

    猜你喜歡
    宿州市日照時(shí)數(shù)冬小麥
    福州市近70年日照變化趨勢(shì)分析
    宿州市昆侖種業(yè)有限公司
    西昌近60年日照時(shí)數(shù)的變化特征分析
    1961~2020年曲麻萊縣日照時(shí)數(shù)變化特征
    拼 搏
    寶藏(2018年10期)2018-10-18 02:29:20
    攬 月
    寶藏(2018年10期)2018-10-18 02:29:20
    1980年~2017年大冶市日照時(shí)數(shù)變化特征分析
    甘肅冬小麥田
    宿州市2000國(guó)家坐標(biāo)系轉(zhuǎn)換方法的研究
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    国产在线免费精品| 免费高清在线观看视频在线观看| 天天操日日干夜夜撸| 国产在线视频一区二区| 久久女婷五月综合色啪小说| av网站免费在线观看视频| 亚洲国产av影院在线观看| 黑人欧美特级aaaaaa片| 午夜日韩欧美国产| 国产精品久久久久久精品电影小说| a 毛片基地| 欧美人与性动交α欧美软件| 欧美成人午夜精品| 在线永久观看黄色视频| 亚洲av日韩在线播放| 国产麻豆69| 国产精品自产拍在线观看55亚洲 | 国产日韩欧美亚洲二区| 19禁男女啪啪无遮挡网站| 国产精品 国内视频| 99久久国产精品久久久| 看免费av毛片| 女人爽到高潮嗷嗷叫在线视频| 不卡一级毛片| 一二三四社区在线视频社区8| 母亲3免费完整高清在线观看| 免费不卡黄色视频| 咕卡用的链子| 十八禁高潮呻吟视频| 啦啦啦 在线观看视频| 超色免费av| 午夜久久久在线观看| 亚洲精品国产av成人精品| 久久久久久人人人人人| 亚洲精品国产精品久久久不卡| 十八禁人妻一区二区| 熟女少妇亚洲综合色aaa.| 天天躁日日躁夜夜躁夜夜| 国产欧美日韩一区二区三区在线| 在线天堂中文资源库| 超色免费av| 欧美少妇被猛烈插入视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美老熟妇乱子伦牲交| 国产成人精品久久二区二区91| 夜夜夜夜夜久久久久| 啦啦啦中文免费视频观看日本| 97人妻天天添夜夜摸| 极品少妇高潮喷水抽搐| 国产精品.久久久| 嫁个100分男人电影在线观看| 老司机午夜福利在线观看视频 | 国产高清videossex| 亚洲精品国产av成人精品| 男人舔女人的私密视频| 午夜福利,免费看| 国产亚洲精品久久久久5区| 欧美日韩一级在线毛片| 99久久国产精品久久久| 欧美精品人与动牲交sv欧美| 免费在线观看视频国产中文字幕亚洲 | 人人妻人人澡人人看| 久久久久视频综合| 一本—道久久a久久精品蜜桃钙片| 欧美人与性动交α欧美精品济南到| 中亚洲国语对白在线视频| 欧美精品一区二区大全| 亚洲国产欧美一区二区综合| 美女大奶头黄色视频| 免费看十八禁软件| 国产精品偷伦视频观看了| 电影成人av| 美女大奶头黄色视频| 亚洲精品自拍成人| 最新的欧美精品一区二区| 中亚洲国语对白在线视频| 中文字幕色久视频| 又黄又粗又硬又大视频| 极品少妇高潮喷水抽搐| 国产一区二区三区在线臀色熟女 | 在线观看免费高清a一片| 婷婷色av中文字幕| 女人久久www免费人成看片| 国产熟女午夜一区二区三区| av国产精品久久久久影院| 婷婷色av中文字幕| 90打野战视频偷拍视频| 99热网站在线观看| 国产在视频线精品| 国产免费av片在线观看野外av| 母亲3免费完整高清在线观看| 午夜久久久在线观看| 国产精品久久久久久精品电影小说| 国产xxxxx性猛交| 少妇裸体淫交视频免费看高清 | 老熟妇仑乱视频hdxx| 日韩,欧美,国产一区二区三区| 成人手机av| 色视频在线一区二区三区| 日韩一区二区三区影片| 欧美激情极品国产一区二区三区| 日韩制服骚丝袜av| av不卡在线播放| 丝袜喷水一区| 亚洲一区二区三区欧美精品| 黑人欧美特级aaaaaa片| 日韩制服丝袜自拍偷拍| 免费观看人在逋| 亚洲成国产人片在线观看| 欧美人与性动交α欧美软件| videos熟女内射| 精品乱码久久久久久99久播| 久久久国产成人免费| 亚洲色图综合在线观看| 老熟妇仑乱视频hdxx| 性高湖久久久久久久久免费观看| 亚洲免费av在线视频| netflix在线观看网站| 午夜免费成人在线视频| 永久免费av网站大全| av网站在线播放免费| 97人妻天天添夜夜摸| 另类亚洲欧美激情| 如日韩欧美国产精品一区二区三区| 国产精品久久久久久精品电影小说| 精品亚洲成国产av| cao死你这个sao货| 欧美久久黑人一区二区| 12—13女人毛片做爰片一| 十八禁网站免费在线| 欧美日韩一级在线毛片| 欧美日韩国产mv在线观看视频| 亚洲av成人不卡在线观看播放网 | 久久精品国产综合久久久| av线在线观看网站| 精品国产乱码久久久久久男人| 女人精品久久久久毛片| 日韩 亚洲 欧美在线| 午夜福利一区二区在线看| 欧美精品一区二区大全| 亚洲欧美日韩高清在线视频 | 久久久精品免费免费高清| av超薄肉色丝袜交足视频| 高清在线国产一区| 另类精品久久| 亚洲成av片中文字幕在线观看| 成人黄色视频免费在线看| 窝窝影院91人妻| 黑人巨大精品欧美一区二区蜜桃| 国产福利在线免费观看视频| 欧美另类亚洲清纯唯美| 国产精品久久久久久精品古装| 99国产精品免费福利视频| 别揉我奶头~嗯~啊~动态视频 | 母亲3免费完整高清在线观看| 黑人欧美特级aaaaaa片| 国产精品麻豆人妻色哟哟久久| 天天操日日干夜夜撸| 欧美久久黑人一区二区| 国产精品久久久av美女十八| 国产一区二区三区在线臀色熟女 | 久久久久精品国产欧美久久久 | 久久天躁狠狠躁夜夜2o2o| 我的亚洲天堂| 青春草视频在线免费观看| 十八禁人妻一区二区| 一级毛片精品| 美女午夜性视频免费| 久久精品成人免费网站| 日日爽夜夜爽网站| 久久国产精品男人的天堂亚洲| 在线看a的网站| 少妇人妻久久综合中文| 日日摸夜夜添夜夜添小说| 久久久久久免费高清国产稀缺| 国产三级黄色录像| 色婷婷av一区二区三区视频| 啦啦啦视频在线资源免费观看| 丝瓜视频免费看黄片| 亚洲欧美一区二区三区黑人| av国产精品久久久久影院| 中文字幕制服av| 岛国毛片在线播放| 又大又爽又粗| 桃花免费在线播放| 国产精品亚洲av一区麻豆| 一级毛片电影观看| 一级黄色大片毛片| 乱人伦中国视频| 国产精品自产拍在线观看55亚洲 | 男人添女人高潮全过程视频| 一个人免费在线观看的高清视频 | 一级毛片精品| 国产免费av片在线观看野外av| 黄片大片在线免费观看| 日本av免费视频播放| 97精品久久久久久久久久精品| 亚洲中文字幕日韩| 国产亚洲av高清不卡| 纯流量卡能插随身wifi吗| 老熟妇仑乱视频hdxx| 免费黄频网站在线观看国产| 久久精品aⅴ一区二区三区四区| 国产在线视频一区二区| 国内毛片毛片毛片毛片毛片| 免费观看a级毛片全部| 亚洲国产欧美在线一区| 91老司机精品| 又黄又粗又硬又大视频| av不卡在线播放| 波多野结衣av一区二区av| 日韩,欧美,国产一区二区三区| 久久精品成人免费网站| 色综合欧美亚洲国产小说| 亚洲va日本ⅴa欧美va伊人久久 | 日韩欧美免费精品| 成人三级做爰电影| 亚洲精品中文字幕在线视频| 亚洲精品粉嫩美女一区| 日韩中文字幕视频在线看片| 久久久久国产精品人妻一区二区| 国产一区二区三区在线臀色熟女 | 男人舔女人的私密视频| 国产一区二区三区在线臀色熟女 | 免费日韩欧美在线观看| 久久毛片免费看一区二区三区| 欧美人与性动交α欧美精品济南到| 9191精品国产免费久久| 亚洲成人免费av在线播放| 免费高清在线观看视频在线观看| 亚洲国产欧美网| 另类亚洲欧美激情| 国产成人一区二区三区免费视频网站| 国产精品国产三级国产专区5o| 国产欧美日韩精品亚洲av| 色综合欧美亚洲国产小说| 狂野欧美激情性bbbbbb| 高清视频免费观看一区二区| 亚洲精品国产一区二区精华液| 一级a爱视频在线免费观看| 成人免费观看视频高清| 成人国产av品久久久| 青草久久国产| 国产av又大| 免费少妇av软件| av不卡在线播放| 99热国产这里只有精品6| 国产一区二区激情短视频 | 亚洲成人手机| 亚洲欧洲日产国产| 99热网站在线观看| 90打野战视频偷拍视频| 免费黄频网站在线观看国产| 黄色视频在线播放观看不卡| 精品国产乱子伦一区二区三区 | 亚洲精品久久久久久婷婷小说| 亚洲九九香蕉| 宅男免费午夜| 国产日韩欧美在线精品| 在线亚洲精品国产二区图片欧美| 亚洲精品国产av成人精品| 欧美大码av| 欧美一级毛片孕妇| 男人操女人黄网站| 无遮挡黄片免费观看| 99久久国产精品久久久| 免费在线观看日本一区| 日韩欧美一区视频在线观看| 亚洲免费av在线视频| 少妇粗大呻吟视频| 中文字幕高清在线视频| 久久影院123| av片东京热男人的天堂| 亚洲三区欧美一区| 91精品三级在线观看| 国产男女内射视频| 久久人妻福利社区极品人妻图片| 欧美av亚洲av综合av国产av| 丁香六月欧美| 久久久精品国产亚洲av高清涩受| 天堂中文最新版在线下载| av欧美777| 欧美精品亚洲一区二区| 91精品国产国语对白视频| e午夜精品久久久久久久| 久久免费观看电影| 一二三四社区在线视频社区8| 欧美精品高潮呻吟av久久| 色播在线永久视频| 日韩欧美一区视频在线观看| 久久人妻福利社区极品人妻图片| 69av精品久久久久久 | 两个人免费观看高清视频| 我要看黄色一级片免费的| 欧美中文综合在线视频| 一区二区av电影网| 国产黄色免费在线视频| 成年美女黄网站色视频大全免费| 在线天堂中文资源库| 高清视频免费观看一区二区| 国产又爽黄色视频| 国产高清视频在线播放一区 | kizo精华| 日韩视频在线欧美| 极品少妇高潮喷水抽搐| 成人av一区二区三区在线看 | 日韩视频在线欧美| 国产极品粉嫩免费观看在线| 美女主播在线视频| 成人黄色视频免费在线看| 色精品久久人妻99蜜桃| 精品国产乱码久久久久久男人| 岛国毛片在线播放| 国产欧美日韩精品亚洲av| 亚洲国产毛片av蜜桃av| 欧美+亚洲+日韩+国产| 波多野结衣av一区二区av| 亚洲精品中文字幕在线视频| 美女国产高潮福利片在线看| 国产在视频线精品| 精品国产乱子伦一区二区三区 | 国产成人av教育| 欧美激情高清一区二区三区| 高清视频免费观看一区二区| 国产有黄有色有爽视频| 美女午夜性视频免费| 狠狠精品人妻久久久久久综合| 91精品伊人久久大香线蕉| 国产成人av激情在线播放| 国产亚洲午夜精品一区二区久久| 人人妻人人澡人人看| 亚洲精品日韩在线中文字幕| 亚洲性夜色夜夜综合| 亚洲男人天堂网一区| 精品视频人人做人人爽| 丝袜在线中文字幕| 99热全是精品| 欧美人与性动交α欧美软件| 丁香六月欧美| 亚洲精品久久午夜乱码| 伦理电影免费视频| 国产精品亚洲av一区麻豆| 中文字幕人妻熟女乱码| 99热网站在线观看| 精品久久蜜臀av无| 亚洲第一欧美日韩一区二区三区 | 9色porny在线观看| 亚洲情色 制服丝袜| 我的亚洲天堂| 亚洲av片天天在线观看| 女人被躁到高潮嗷嗷叫费观| 国产真人三级小视频在线观看| 亚洲av日韩精品久久久久久密| 久久天躁狠狠躁夜夜2o2o| 一本—道久久a久久精品蜜桃钙片| 久久久久久久大尺度免费视频| 国产深夜福利视频在线观看| 欧美变态另类bdsm刘玥| 久久国产精品大桥未久av| 91国产中文字幕| 人妻久久中文字幕网| 丁香六月欧美| 午夜91福利影院| 日韩欧美一区二区三区在线观看 | 国产成人av教育| 免费av中文字幕在线| 美女国产高潮福利片在线看| 欧美日韩亚洲高清精品| 精品一品国产午夜福利视频| 日韩视频一区二区在线观看| 最近最新免费中文字幕在线| 大片电影免费在线观看免费| 亚洲自偷自拍图片 自拍| 男人爽女人下面视频在线观看| 99久久人妻综合| 午夜福利免费观看在线| 国产色视频综合| 国产人伦9x9x在线观看| 咕卡用的链子| 婷婷成人精品国产| 国产欧美日韩精品亚洲av| 最新的欧美精品一区二区| 成人国语在线视频| 国产1区2区3区精品| 黑人巨大精品欧美一区二区mp4| 少妇被粗大的猛进出69影院| 日韩欧美免费精品| 精品一区二区三区四区五区乱码| 一级毛片女人18水好多| 精品少妇一区二区三区视频日本电影| 久久精品aⅴ一区二区三区四区| 国产有黄有色有爽视频| 午夜老司机福利片| 亚洲av电影在线观看一区二区三区| 久久久久网色| 国产一区有黄有色的免费视频| 一二三四社区在线视频社区8| 精品免费久久久久久久清纯 | 大香蕉久久成人网| 精品第一国产精品| 亚洲欧美激情在线| 巨乳人妻的诱惑在线观看| 国产欧美亚洲国产| videosex国产| 亚洲国产看品久久| 免费少妇av软件| 99国产综合亚洲精品| 日韩有码中文字幕| 可以免费在线观看a视频的电影网站| 亚洲色图综合在线观看| 亚洲欧美日韩另类电影网站| 他把我摸到了高潮在线观看 | 国产区一区二久久| 高潮久久久久久久久久久不卡| a 毛片基地| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av日韩精品久久久久久密| 午夜福利在线观看吧| 一级片'在线观看视频| 精品卡一卡二卡四卡免费| 欧美日韩黄片免| 亚洲第一av免费看| 亚洲一码二码三码区别大吗| 国产成人一区二区三区免费视频网站| 中文字幕高清在线视频| 黄色a级毛片大全视频| 99热全是精品| 热99久久久久精品小说推荐| 欧美精品一区二区免费开放| 久久影院123| 黑人巨大精品欧美一区二区mp4| 日韩视频在线欧美| 人成视频在线观看免费观看| 国产片内射在线| 少妇 在线观看| 日韩欧美一区二区三区在线观看 | 欧美日韩亚洲综合一区二区三区_| 免费高清在线观看视频在线观看| 国产97色在线日韩免费| 少妇猛男粗大的猛烈进出视频| av在线app专区| 美女脱内裤让男人舔精品视频| bbb黄色大片| 啦啦啦视频在线资源免费观看| 欧美亚洲日本最大视频资源| 亚洲男人天堂网一区| 一区二区av电影网| 99九九在线精品视频| 激情视频va一区二区三区| 久久狼人影院| 欧美日韩亚洲综合一区二区三区_| 亚洲国产av影院在线观看| 99国产精品99久久久久| 纯流量卡能插随身wifi吗| 亚洲欧美一区二区三区黑人| 各种免费的搞黄视频| 午夜91福利影院| 日韩精品免费视频一区二区三区| tocl精华| xxxhd国产人妻xxx| 少妇粗大呻吟视频| 亚洲性夜色夜夜综合| 亚洲色图综合在线观看| 在线观看人妻少妇| 如日韩欧美国产精品一区二区三区| 97人妻天天添夜夜摸| 亚洲成人免费av在线播放| 国产免费福利视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 建设人人有责人人尽责人人享有的| 久久毛片免费看一区二区三区| 最新的欧美精品一区二区| 高潮久久久久久久久久久不卡| 丁香六月天网| 永久免费av网站大全| 一级毛片电影观看| 啦啦啦啦在线视频资源| 男女之事视频高清在线观看| 欧美一级毛片孕妇| 水蜜桃什么品种好| 80岁老熟妇乱子伦牲交| 热99re8久久精品国产| 一本大道久久a久久精品| 亚洲国产看品久久| 91av网站免费观看| 中文精品一卡2卡3卡4更新| 黄色视频在线播放观看不卡| 国产淫语在线视频| 亚洲国产精品一区二区三区在线| 9色porny在线观看| 成年av动漫网址| 香蕉丝袜av| 成年女人毛片免费观看观看9 | 一级毛片精品| 嫩草影视91久久| 欧美激情高清一区二区三区| √禁漫天堂资源中文www| 老司机在亚洲福利影院| 色94色欧美一区二区| 国产成人免费观看mmmm| 97精品久久久久久久久久精品| 大陆偷拍与自拍| 成人三级做爰电影| 美女视频免费永久观看网站| 丝袜美腿诱惑在线| 午夜免费鲁丝| 国产日韩欧美亚洲二区| 欧美激情高清一区二区三区| 色视频在线一区二区三区| 国产在线一区二区三区精| 51午夜福利影视在线观看| 在线av久久热| 天堂俺去俺来也www色官网| 亚洲欧洲精品一区二区精品久久久| 黑人巨大精品欧美一区二区蜜桃| 一级,二级,三级黄色视频| 欧美日韩国产mv在线观看视频| 欧美97在线视频| 欧美日韩成人在线一区二区| 搡老熟女国产l中国老女人| 巨乳人妻的诱惑在线观看| 美女大奶头黄色视频| 人人澡人人妻人| 美女中出高潮动态图| 久久性视频一级片| 我要看黄色一级片免费的| 亚洲av日韩精品久久久久久密| 一本—道久久a久久精品蜜桃钙片| 久久久久网色| 男女床上黄色一级片免费看| svipshipincom国产片| 妹子高潮喷水视频| 性高湖久久久久久久久免费观看| 国产欧美日韩一区二区精品| 男人爽女人下面视频在线观看| 99国产精品免费福利视频| 色婷婷av一区二区三区视频| 亚洲精品中文字幕在线视频| 免费高清在线观看日韩| 日本猛色少妇xxxxx猛交久久| 色婷婷久久久亚洲欧美| 亚洲国产日韩一区二区| 亚洲综合色网址| 久久久久久免费高清国产稀缺| 国产免费福利视频在线观看| 中国国产av一级| 亚洲全国av大片| 免费观看av网站的网址| 亚洲熟女精品中文字幕| 国产片内射在线| 国产老妇伦熟女老妇高清| 国产亚洲一区二区精品| av在线老鸭窝| 人妻一区二区av| 中国美女看黄片| 日日夜夜操网爽| 精品人妻1区二区| 69精品国产乱码久久久| a级毛片黄视频| 成年人免费黄色播放视频| 2018国产大陆天天弄谢| 999精品在线视频| 亚洲欧美一区二区三区黑人| 国产一区二区 视频在线| 巨乳人妻的诱惑在线观看| 中文字幕人妻熟女乱码| a级毛片黄视频| 国产欧美日韩精品亚洲av| 黑人欧美特级aaaaaa片| 女人高潮潮喷娇喘18禁视频| 一区二区三区精品91| 我要看黄色一级片免费的| 男女床上黄色一级片免费看| av又黄又爽大尺度在线免费看| 精品亚洲成国产av| 成年人黄色毛片网站| 在线 av 中文字幕| 黑丝袜美女国产一区| 一本—道久久a久久精品蜜桃钙片| 老熟妇乱子伦视频在线观看 | 亚洲精品国产区一区二| 精品免费久久久久久久清纯 | 50天的宝宝边吃奶边哭怎么回事| 1024香蕉在线观看| 欧美一级毛片孕妇| 悠悠久久av| 日本av免费视频播放| 国产精品久久久久久精品电影小说| 免费高清在线观看日韩| 精品国产乱子伦一区二区三区 | 91老司机精品| 日韩电影二区| av网站在线播放免费| 免费少妇av软件| av网站在线播放免费| 精品人妻在线不人妻| 丰满人妻熟妇乱又伦精品不卡| 国产野战对白在线观看| 欧美在线一区亚洲| 国产在视频线精品| 午夜福利在线免费观看网站| 91老司机精品| 欧美老熟妇乱子伦牲交| 欧美久久黑人一区二区| 亚洲精品美女久久av网站| 亚洲色图综合在线观看| 国产在线免费精品| 国产一区二区 视频在线| 人妻一区二区av| svipshipincom国产片| 久久久欧美国产精品| 欧美中文综合在线视频| 巨乳人妻的诱惑在线观看| 夜夜夜夜夜久久久久| 欧美 日韩 精品 国产|