• <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
    女性生殖器流出的白浆| 观看美女的网站| 国内少妇人妻偷人精品xxx网站| 亚洲av二区三区四区| 国产伦理片在线播放av一区| 国产精品麻豆人妻色哟哟久久| 久久久a久久爽久久v久久| 国产视频内射| 啦啦啦啦在线视频资源| 日本黄色片子视频| 国产精品久久久久久精品古装| 色视频www国产| 热re99久久精品国产66热6| 嘟嘟电影网在线观看| 秋霞在线观看毛片| 建设人人有责人人尽责人人享有的 | 在线观看av片永久免费下载| 国产白丝娇喘喷水9色精品| 在线观看一区二区三区激情| 日韩一区二区视频免费看| 免费看日本二区| 国产91av在线免费观看| 在线观看人妻少妇| 晚上一个人看的免费电影| 少妇精品久久久久久久| 成人美女网站在线观看视频| 精品一区二区三卡| videos熟女内射| 日韩视频在线欧美| 午夜视频国产福利| 亚洲欧美成人精品一区二区| 国产毛片在线视频| 国产伦在线观看视频一区| 精品少妇黑人巨大在线播放| 97超碰精品成人国产| 色5月婷婷丁香| 国产在线视频一区二区| 国产白丝娇喘喷水9色精品| 黑丝袜美女国产一区| 中文天堂在线官网| 欧美三级亚洲精品| 国产精品一区二区性色av| 国产伦在线观看视频一区| 亚洲中文av在线| 亚洲第一区二区三区不卡| 亚洲一区二区三区欧美精品| av在线观看视频网站免费| 日本午夜av视频| 亚洲国产av新网站| 亚洲一级一片aⅴ在线观看| 国产成人免费观看mmmm| 777米奇影视久久| 啦啦啦在线观看免费高清www| 欧美精品国产亚洲| 免费久久久久久久精品成人欧美视频 | 80岁老熟妇乱子伦牲交| 久久精品国产亚洲av涩爱| 久久久久久久久久久免费av| 欧美日韩综合久久久久久| 边亲边吃奶的免费视频| 只有这里有精品99| 亚洲成人中文字幕在线播放| 美女国产视频在线观看| 亚洲欧洲日产国产| 一级二级三级毛片免费看| 亚洲,一卡二卡三卡| av国产久精品久网站免费入址| 国产在视频线精品| 少妇 在线观看| 亚洲精品成人av观看孕妇| 亚洲av在线观看美女高潮| 联通29元200g的流量卡| 赤兔流量卡办理| 看免费成人av毛片| 黄色日韩在线| 在线看a的网站| av在线观看视频网站免费| 美女国产视频在线观看| 97在线视频观看| 精品一区在线观看国产| 免费大片18禁| 国产精品嫩草影院av在线观看| 国产精品一二三区在线看| 欧美成人午夜免费资源| 男人添女人高潮全过程视频| 国产成人一区二区在线| 高清视频免费观看一区二区| 亚洲成色77777| 亚洲av国产av综合av卡| 九九久久精品国产亚洲av麻豆| 777米奇影视久久| 黄色怎么调成土黄色| 欧美成人一区二区免费高清观看| 少妇 在线观看| 亚洲图色成人| a级一级毛片免费在线观看| 中文资源天堂在线| 久久国产精品大桥未久av | 永久网站在线| 亚洲精品视频女| 亚洲av成人精品一二三区| 久久久国产一区二区| 天堂8中文在线网| 91精品伊人久久大香线蕉| av女优亚洲男人天堂| 欧美日韩亚洲高清精品| 精品亚洲乱码少妇综合久久| h视频一区二区三区| 伦理电影大哥的女人| 国产女主播在线喷水免费视频网站| av播播在线观看一区| 中国三级夫妇交换| 美女xxoo啪啪120秒动态图| 最后的刺客免费高清国语| 少妇熟女欧美另类| 亚洲图色成人| 中国美白少妇内射xxxbb| 身体一侧抽搐| 日韩伦理黄色片| 少妇丰满av| 亚洲精品日韩av片在线观看| 日韩精品有码人妻一区| 国产成人免费观看mmmm| 国产综合精华液| 国内精品宾馆在线| 狠狠精品人妻久久久久久综合| 啦啦啦啦在线视频资源| 精品久久久噜噜| 成人免费观看视频高清| av在线观看视频网站免费| 亚洲av福利一区| 国产成人免费无遮挡视频| 国产美女午夜福利| 亚洲图色成人| 国产伦理片在线播放av一区| 国产成人精品久久久久久| 国产精品无大码| 高清av免费在线| 亚洲aⅴ乱码一区二区在线播放| 男人爽女人下面视频在线观看| 女的被弄到高潮叫床怎么办| 欧美97在线视频| h日本视频在线播放| 中文天堂在线官网| 久久久久久伊人网av| 色视频在线一区二区三区| 午夜激情福利司机影院| 亚洲图色成人| 亚洲欧美日韩东京热| freevideosex欧美| 精品久久久久久电影网| a级毛片免费高清观看在线播放| 97热精品久久久久久| 啦啦啦啦在线视频资源| 国产欧美日韩精品一区二区| 久久女婷五月综合色啪小说| 只有这里有精品99| 高清av免费在线| 18禁裸乳无遮挡免费网站照片| 欧美 日韩 精品 国产| 国产午夜精品一二区理论片| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 国产精品欧美亚洲77777| 一区二区三区乱码不卡18| 欧美丝袜亚洲另类| 国产淫语在线视频| 亚洲电影在线观看av| 亚洲av成人精品一二三区| 毛片一级片免费看久久久久| 亚洲国产成人一精品久久久| 亚洲国产精品999| 男女啪啪激烈高潮av片| 久久精品久久久久久噜噜老黄| 黄色怎么调成土黄色| 亚洲婷婷狠狠爱综合网| 简卡轻食公司| 国产一级毛片在线| 精品亚洲成a人片在线观看 | 丰满乱子伦码专区| 久久久久精品性色| 久久久a久久爽久久v久久| 久久国产精品大桥未久av | 少妇人妻久久综合中文| 精品少妇久久久久久888优播| 亚洲欧美一区二区三区国产| 在线观看人妻少妇| 大香蕉久久网| 日本与韩国留学比较| 高清av免费在线| 一区二区三区乱码不卡18| 亚洲av不卡在线观看| 精品一品国产午夜福利视频| 高清午夜精品一区二区三区| 国产成人91sexporn| 国产精品女同一区二区软件| 久久精品夜色国产| 亚洲国产毛片av蜜桃av| 亚洲成色77777| 欧美日韩精品成人综合77777| 久久人人爽av亚洲精品天堂 | 国产黄频视频在线观看| 国产视频首页在线观看| av不卡在线播放| 老女人水多毛片| 国产精品一区二区在线观看99| 精品久久久久久电影网| 免费观看无遮挡的男女| 午夜精品国产一区二区电影| 人妻制服诱惑在线中文字幕| 精品人妻视频免费看| 青春草亚洲视频在线观看| 狂野欧美白嫩少妇大欣赏| 男人添女人高潮全过程视频| 欧美日韩综合久久久久久| 十分钟在线观看高清视频www | av又黄又爽大尺度在线免费看| 高清欧美精品videossex| 国产一区有黄有色的免费视频| 五月开心婷婷网| 在线观看一区二区三区| 精品少妇久久久久久888优播| 三级国产精品片| 国产精品国产av在线观看| 久久热精品热| 国产精品免费大片| 男人和女人高潮做爰伦理| 亚洲第一av免费看| 高清黄色对白视频在线免费看 | 一级毛片电影观看| 少妇人妻精品综合一区二区| 又黄又爽又刺激的免费视频.| 亚洲成人av在线免费| 最黄视频免费看| 国产黄片视频在线免费观看| 成年免费大片在线观看| 午夜免费男女啪啪视频观看| 日韩免费高清中文字幕av| 一级黄片播放器| 夜夜爽夜夜爽视频| 国产成人精品久久久久久| 亚洲精品乱久久久久久| 亚洲欧美清纯卡通| 精品亚洲乱码少妇综合久久| 国产精品.久久久| 又大又黄又爽视频免费| 91在线精品国自产拍蜜月| 啦啦啦啦在线视频资源| 黄色配什么色好看| 看非洲黑人一级黄片| 欧美日韩国产mv在线观看视频 | 欧美日韩精品成人综合77777| 中文天堂在线官网| 久久影院123| a级毛色黄片| 18禁动态无遮挡网站| 91狼人影院| 3wmmmm亚洲av在线观看| 色哟哟·www| 亚洲人成网站在线播| 人妻夜夜爽99麻豆av| 人人妻人人添人人爽欧美一区卜 | 午夜激情福利司机影院| 国产男女超爽视频在线观看| 一级a做视频免费观看| 如何舔出高潮| 亚洲成人av在线免费| 26uuu在线亚洲综合色| 婷婷色综合www| 久久亚洲国产成人精品v| 美女主播在线视频| 亚洲va在线va天堂va国产| 狠狠精品人妻久久久久久综合| 国产熟女欧美一区二区| 成人国产av品久久久| 一二三四中文在线观看免费高清| 在线观看一区二区三区激情| 国产 一区精品| 国语对白做爰xxxⅹ性视频网站| 国产淫语在线视频| 黄片wwwwww| 一本一本综合久久| 亚洲久久久国产精品| 亚洲自偷自拍三级| 婷婷色av中文字幕| 九九在线视频观看精品| 丰满乱子伦码专区| 免费黄网站久久成人精品| 22中文网久久字幕| 国产精品精品国产色婷婷| 精品亚洲成国产av| 国产黄色免费在线视频| 亚洲,一卡二卡三卡| 国产免费又黄又爽又色| 这个男人来自地球电影免费观看 | 青春草视频在线免费观看| 在线观看av片永久免费下载| 丰满迷人的少妇在线观看| 国产乱人偷精品视频| 欧美精品一区二区免费开放| 丰满少妇做爰视频| 久久毛片免费看一区二区三区| 日本色播在线视频| 亚洲av免费高清在线观看| 十八禁网站网址无遮挡 | 黄色怎么调成土黄色| 国产日韩欧美在线精品| 亚洲美女黄色视频免费看| 国模一区二区三区四区视频| 久久久精品免费免费高清| 国产精品秋霞免费鲁丝片| 国产精品久久久久久精品电影小说 | 国产精品久久久久久久电影| 免费大片18禁| 一二三四中文在线观看免费高清| 中文字幕免费在线视频6| 久久久久久久久久久免费av| 日韩三级伦理在线观看| 97精品久久久久久久久久精品| 人人妻人人澡人人爽人人夜夜| 亚州av有码| 99久国产av精品国产电影| 国产精品不卡视频一区二区| 久久久久精品久久久久真实原创| 狂野欧美白嫩少妇大欣赏| 成人亚洲精品一区在线观看 | 亚洲av电影在线观看一区二区三区| 亚洲欧美一区二区三区国产| 大香蕉久久网| 多毛熟女@视频| av卡一久久| 亚洲丝袜综合中文字幕| 国产午夜精品久久久久久一区二区三区| 欧美成人a在线观看| 亚洲人与动物交配视频| 国产男女超爽视频在线观看| 伦理电影免费视频| 午夜精品国产一区二区电影| 国产乱人视频| 嫩草影院入口| 亚洲精华国产精华液的使用体验| 大码成人一级视频| 青春草国产在线视频| 永久网站在线| 一级二级三级毛片免费看| 最近最新中文字幕大全电影3| 久久久久精品久久久久真实原创| 五月天丁香电影| 亚洲久久久国产精品| 国产亚洲欧美精品永久| 婷婷色av中文字幕| 国产精品一区二区性色av| 欧美日韩视频高清一区二区三区二| 亚洲国产毛片av蜜桃av| 五月伊人婷婷丁香| 在线 av 中文字幕| 国产欧美亚洲国产| 国产精品免费大片| videos熟女内射| 五月伊人婷婷丁香| 简卡轻食公司| 伦理电影大哥的女人| 日韩中文字幕视频在线看片 | 日本欧美视频一区| 亚洲精品久久久久久婷婷小说| 美女内射精品一级片tv| 在线观看一区二区三区激情| 久久精品久久久久久噜噜老黄| 久久人人爽人人片av| 寂寞人妻少妇视频99o| 亚洲色图综合在线观看| 80岁老熟妇乱子伦牲交| 亚洲欧洲国产日韩| 能在线免费看毛片的网站| 秋霞伦理黄片| 国内少妇人妻偷人精品xxx网站| 午夜免费鲁丝| 国产v大片淫在线免费观看| 国产av码专区亚洲av| 亚洲av福利一区| 女的被弄到高潮叫床怎么办| 日韩电影二区| 久久99精品国语久久久| 亚洲国产最新在线播放| 18+在线观看网站| 午夜福利在线在线| av一本久久久久| 午夜福利视频精品| 久久精品国产亚洲av天美| 男的添女的下面高潮视频| 亚洲国产精品一区三区| 国产成人aa在线观看| 亚洲色图av天堂| 精品一区二区免费观看| 99久国产av精品国产电影| 蜜桃久久精品国产亚洲av| 女性生殖器流出的白浆| 精品久久国产蜜桃| 亚洲欧美中文字幕日韩二区| 在线观看免费视频网站a站| 天天躁夜夜躁狠狠久久av| 天美传媒精品一区二区| 国产日韩欧美亚洲二区| 夜夜看夜夜爽夜夜摸| 国产精品欧美亚洲77777| 日韩一区二区三区影片| 久久av网站| 少妇被粗大猛烈的视频| 在线免费十八禁| av天堂中文字幕网| 欧美性感艳星| 久久精品熟女亚洲av麻豆精品| 五月伊人婷婷丁香| 亚洲欧美日韩另类电影网站 | 国产精品偷伦视频观看了| 亚洲欧美日韩另类电影网站 | 男女边吃奶边做爰视频| av在线app专区| 欧美日韩视频精品一区| 欧美成人a在线观看| 赤兔流量卡办理| 久久久午夜欧美精品| 五月伊人婷婷丁香| 青青草视频在线视频观看| 国产精品伦人一区二区| 岛国毛片在线播放| a级毛片免费高清观看在线播放| 国产爽快片一区二区三区| 国产精品偷伦视频观看了| 国产淫片久久久久久久久| 国产精品一区二区三区四区免费观看| 狂野欧美激情性xxxx在线观看| 成人综合一区亚洲| 亚洲美女视频黄频| 超碰97精品在线观看| 亚洲第一区二区三区不卡| 亚州av有码| 性色av一级| 亚洲欧美精品自产自拍| 欧美日韩视频高清一区二区三区二| 午夜免费观看性视频| 日日撸夜夜添| 建设人人有责人人尽责人人享有的 | 国产综合精华液| 九九爱精品视频在线观看| 免费观看在线日韩| 日本av手机在线免费观看| 日韩一区二区视频免费看| 777米奇影视久久| 丝瓜视频免费看黄片| 欧美xxⅹ黑人| 亚洲,一卡二卡三卡| 国产伦理片在线播放av一区| 亚洲内射少妇av| 午夜免费男女啪啪视频观看| 久久精品久久精品一区二区三区| 国产熟女欧美一区二区| 少妇人妻精品综合一区二区| 国产一区二区在线观看日韩| 涩涩av久久男人的天堂| 国国产精品蜜臀av免费| 久久鲁丝午夜福利片| 18禁裸乳无遮挡免费网站照片| 亚洲av日韩在线播放| 亚洲精品第二区| 黄片wwwwww| 观看免费一级毛片| 国产视频首页在线观看| 三级经典国产精品| 在现免费观看毛片| 欧美精品一区二区免费开放| 97在线人人人人妻| 一级a做视频免费观看| 亚洲,欧美,日韩| 一区二区三区免费毛片| 22中文网久久字幕| 中文字幕亚洲精品专区| 爱豆传媒免费全集在线观看| 热99国产精品久久久久久7| 成年美女黄网站色视频大全免费 | 国产免费视频播放在线视频| 亚洲av日韩在线播放| 亚洲欧美日韩另类电影网站 | 欧美xxxx性猛交bbbb| 国产 一区精品| 精品亚洲成a人片在线观看 | 女人久久www免费人成看片| 国产亚洲精品久久久com| 日韩av不卡免费在线播放| 美女视频免费永久观看网站| 各种免费的搞黄视频| 人妻 亚洲 视频| 久久人人爽av亚洲精品天堂 | 最近2019中文字幕mv第一页| 精品久久久久久久久亚洲| 蜜桃亚洲精品一区二区三区| 哪个播放器可以免费观看大片| 成年免费大片在线观看| 国产黄频视频在线观看| 亚洲精品456在线播放app| 欧美日韩视频精品一区| 特大巨黑吊av在线直播| 深夜a级毛片| 一级毛片 在线播放| 国产片特级美女逼逼视频| 一级毛片黄色毛片免费观看视频| 午夜福利视频精品| 看免费成人av毛片| 观看美女的网站| 欧美激情国产日韩精品一区| 偷拍熟女少妇极品色| 日本黄色日本黄色录像| 久久久色成人| 日本av免费视频播放| 精品亚洲成国产av| 美女xxoo啪啪120秒动态图| 国产av精品麻豆| 久久精品夜色国产| 菩萨蛮人人尽说江南好唐韦庄| 80岁老熟妇乱子伦牲交| 嫩草影院入口| 亚洲精品日韩在线中文字幕| 丝瓜视频免费看黄片| 97在线视频观看| 国产成人精品福利久久| 精品人妻熟女av久视频| 26uuu在线亚洲综合色| 免费人成在线观看视频色| 观看美女的网站| 成年美女黄网站色视频大全免费 | av一本久久久久| 国产毛片在线视频| 成人影院久久| 亚洲无线观看免费| 超碰av人人做人人爽久久| 免费大片黄手机在线观看| 大片免费播放器 马上看| 国产熟女欧美一区二区| 亚洲av.av天堂| 国产精品久久久久久av不卡| 成人高潮视频无遮挡免费网站| 成人漫画全彩无遮挡| 在线播放无遮挡| 在线观看美女被高潮喷水网站| 欧美+日韩+精品| 国精品久久久久久国模美| 三级经典国产精品| 51国产日韩欧美| www.色视频.com| 久久亚洲国产成人精品v| 亚洲综合精品二区| 国产真实伦视频高清在线观看| 成人美女网站在线观看视频| 亚洲,一卡二卡三卡| 国内少妇人妻偷人精品xxx网站| 亚洲美女搞黄在线观看| 97在线视频观看| 在线精品无人区一区二区三 | 国产色爽女视频免费观看| 一级a做视频免费观看| 另类亚洲欧美激情| 日本爱情动作片www.在线观看| 国产白丝娇喘喷水9色精品| 亚洲精品国产av成人精品| 国内少妇人妻偷人精品xxx网站| a级毛色黄片| 久久久久精品久久久久真实原创| 尤物成人国产欧美一区二区三区| 我的女老师完整版在线观看| 成人美女网站在线观看视频| 丰满乱子伦码专区| 人体艺术视频欧美日本| 国产精品久久久久久精品古装| 啦啦啦在线观看免费高清www| 色综合色国产| 亚洲人与动物交配视频| av免费在线看不卡| 高清视频免费观看一区二区| 亚洲av免费高清在线观看| 精华霜和精华液先用哪个| av黄色大香蕉| 九色成人免费人妻av| 久久亚洲国产成人精品v| 美女cb高潮喷水在线观看| 边亲边吃奶的免费视频| 高清欧美精品videossex| 伦理电影大哥的女人| 日本欧美国产在线视频| freevideosex欧美| 极品教师在线视频| 日韩不卡一区二区三区视频在线| 精品亚洲成国产av| 日本猛色少妇xxxxx猛交久久| 国产伦理片在线播放av一区| 91精品国产国语对白视频| 精品国产一区二区三区久久久樱花 | 国产精品嫩草影院av在线观看| 免费高清在线观看视频在线观看| 亚洲精品国产av蜜桃| 国产欧美另类精品又又久久亚洲欧美| 久久婷婷青草| 校园人妻丝袜中文字幕| 最近中文字幕2019免费版| 麻豆成人午夜福利视频| 亚洲熟女精品中文字幕| 国产大屁股一区二区在线视频| 国产男女内射视频| 国产一区二区三区综合在线观看 | 涩涩av久久男人的天堂| 国产亚洲欧美精品永久| 亚洲丝袜综合中文字幕| 久久99热这里只频精品6学生| 波野结衣二区三区在线| 在线 av 中文字幕| 日韩不卡一区二区三区视频在线|