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

    冬小麥觀測(cè)產(chǎn)量與統(tǒng)計(jì)產(chǎn)量的差異性分析*

    2021-03-10 03:03:40孟翠麗宋迎波
    中國(guó)農(nóng)業(yè)氣象 2021年2期
    關(guān)鍵詞:冬小麥平均值站點(diǎn)

    劉 維,孟翠麗,宋迎波**

    冬小麥觀測(cè)產(chǎn)量與統(tǒng)計(jì)產(chǎn)量的差異性分析*

    劉 維1,孟翠麗2,宋迎波1**

    (1.國(guó)家氣象中心,北京 100081;2.武漢農(nóng)業(yè)氣象試驗(yàn)站,武漢 430040)

    分析對(duì)比全國(guó)123個(gè)冬小麥農(nóng)業(yè)氣象觀測(cè)站1991?2017年觀測(cè)產(chǎn)量與所在縣統(tǒng)計(jì)產(chǎn)量的年代際、變異系數(shù)和傾向率的差異;利用2006?2010年各縣冬小麥種植面積平均值省內(nèi)占比作為權(quán)重因子,將縣觀測(cè)產(chǎn)量與縣統(tǒng)計(jì)產(chǎn)量加權(quán)集成為省級(jí)尺度觀測(cè)產(chǎn)量和省級(jí)尺度統(tǒng)計(jì)產(chǎn)量,并與統(tǒng)計(jì)局發(fā)布的省公布產(chǎn)量進(jìn)行對(duì)比,分析省級(jí)尺度三種不同產(chǎn)量的年代際變化和傾向率差異。結(jié)果表明:(1)縣觀測(cè)產(chǎn)量和縣統(tǒng)計(jì)產(chǎn)量均表現(xiàn)為高產(chǎn)縣數(shù)量增幅明顯,低產(chǎn)縣數(shù)量減幅明顯;21世紀(jì)10年代兩者均為高產(chǎn)年代,21世紀(jì)00年代兩者差值達(dá)到峰值。(2)縣尺度變異系數(shù)觀測(cè)產(chǎn)量離散程度高于縣統(tǒng)計(jì)產(chǎn)量,49個(gè)縣觀測(cè)產(chǎn)量的變異系數(shù)小于0.20,僅8個(gè)站大于0.40;72個(gè)縣統(tǒng)計(jì)產(chǎn)量變異系數(shù)小于0.20,僅9個(gè)站點(diǎn)變異系數(shù)大于0.30。(3)各縣觀測(cè)產(chǎn)量中有73個(gè)站點(diǎn)傾向率呈現(xiàn)顯著增加趨勢(shì),大多集中在河北、河南、山東、江蘇、安徽等冬小麥主產(chǎn)?。?00個(gè)縣統(tǒng)計(jì)產(chǎn)量呈顯著增加趨勢(shì)。(4)21世紀(jì)00年代為各省觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量的高產(chǎn)年代,20世紀(jì)90年代為低產(chǎn)年代;山東、安徽、河北、江蘇、陜西和山西等省10a觀測(cè)產(chǎn)量的平均值在各個(gè)年代均高于省統(tǒng)計(jì)產(chǎn)量平均值。(5)除新疆和山西外其余省份省級(jí)尺度觀測(cè)產(chǎn)量?jī)A向率均通過(guò)顯著性檢驗(yàn);所有省份省統(tǒng)計(jì)產(chǎn)量和公布產(chǎn)量?jī)A向率均通過(guò)顯著性檢驗(yàn),且產(chǎn)量增幅均為正值??傮w來(lái)說(shuō),基于觀測(cè)產(chǎn)量的冬小麥產(chǎn)量序列可以為產(chǎn)量預(yù)報(bào)提供新的數(shù)據(jù)來(lái)源。

    冬小麥;農(nóng)試站;觀測(cè)產(chǎn)量;統(tǒng)計(jì)產(chǎn)量;面積權(quán)重集成

    作物產(chǎn)量預(yù)報(bào)是現(xiàn)代農(nóng)業(yè)氣象業(yè)務(wù)服務(wù)與科研開(kāi)發(fā)的重要工作之一,國(guó)家級(jí)和省級(jí)產(chǎn)量預(yù)報(bào)產(chǎn)品更是各級(jí)政府部門制定農(nóng)業(yè)生產(chǎn)政策重要的參考依據(jù),是氣象為農(nóng)服務(wù)的重要組成部分[1?2]。隨著氣象部門精密監(jiān)測(cè)、精準(zhǔn)預(yù)報(bào)、精細(xì)服務(wù)要求的提出,精細(xì)化服務(wù)成為農(nóng)業(yè)氣象業(yè)務(wù)服務(wù)的必然要求,隨著農(nóng)作物長(zhǎng)勢(shì)監(jiān)測(cè)評(píng)估與預(yù)報(bào)、農(nóng)業(yè)氣象災(zāi)害監(jiān)測(cè)和評(píng)估、農(nóng)用天氣預(yù)報(bào)等農(nóng)業(yè)氣象服務(wù)進(jìn)入公里級(jí)別[2?3],開(kāi)展地市級(jí)和縣級(jí)作物產(chǎn)量預(yù)報(bào),拓展產(chǎn)量預(yù)報(bào)服務(wù)范圍,提升作物產(chǎn)量預(yù)報(bào)的精細(xì)化和精準(zhǔn)度是未來(lái)農(nóng)業(yè)氣象業(yè)務(wù)服務(wù)的發(fā)展趨勢(shì)。目前產(chǎn)量預(yù)報(bào)業(yè)務(wù)應(yīng)用中的產(chǎn)量資料唯一來(lái)源為統(tǒng)計(jì)部門的公布產(chǎn)量,但公布產(chǎn)量獲取存在較長(zhǎng)的滯后性,特別是縣級(jí)尺度的產(chǎn)量序列,獲取難度更大。而全國(guó)氣象部門的農(nóng)業(yè)氣象觀測(cè)體系中有653個(gè)農(nóng)業(yè)氣象觀測(cè)站、70個(gè)農(nóng)業(yè)氣象試驗(yàn)站,可實(shí)時(shí)監(jiān)測(cè)作物不同生育期生長(zhǎng)狀況,且各觀測(cè)站逐年上報(bào)觀測(cè)地段觀測(cè)的單位面積作物產(chǎn)量數(shù)據(jù)(觀測(cè)產(chǎn)量)與所在縣統(tǒng)計(jì)局公布的單位面積的作物產(chǎn)量資料(統(tǒng)計(jì)產(chǎn)量),這也為氣象部門開(kāi)展基于站點(diǎn)觀測(cè)產(chǎn)量的產(chǎn)量預(yù)報(bào)提供了數(shù)據(jù)支撐。當(dāng)前多數(shù)研究集中在利用農(nóng)試站生育期觀測(cè)資料進(jìn)行作物模型模擬[4?7]、災(zāi)害評(píng)估[8?9]、氣候變化對(duì)作物生長(zhǎng)及種植制度的影響[10?13]等方面;吳冰潔等[14]研究了氣象因子對(duì)華北平原冬小麥?zhǔn)屑?jí)產(chǎn)量的影響;趙凱娜等[15]利用灰色關(guān)聯(lián)度分析了河南省縣域冬小麥產(chǎn)量與生育期氣象要素的關(guān)系。當(dāng)前開(kāi)展縣級(jí)尺度冬小麥產(chǎn)量的研究就少,對(duì)農(nóng)試站觀測(cè)產(chǎn)量本身的可用性及與統(tǒng)計(jì)產(chǎn)量的差異研究也未有相關(guān)報(bào)道。

    冬小麥作為最主要的夏收糧食作物,在國(guó)家糧食安全中占有重要地位,加之氣候變化對(duì)冬小麥種植制度及產(chǎn)量的影響在加大[16?17],冬小麥產(chǎn)量預(yù)報(bào)工作顯得尤其重要。當(dāng)前全國(guó)氣象部門共有221個(gè)氣象站能提供冬小麥區(qū)域觀測(cè)產(chǎn)量和當(dāng)?shù)亟y(tǒng)計(jì)產(chǎn)量,因此本研究從歷年冬小麥觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量入手,分析兩者變化規(guī)律,并利用種植面積權(quán)重占比將縣尺度產(chǎn)量集成為省級(jí)尺度產(chǎn)量,并對(duì)集成的省級(jí)產(chǎn)量進(jìn)行對(duì)比,明確觀測(cè)產(chǎn)量的可用性與可預(yù)報(bào)性,為開(kāi)展不同空間尺度的冬小麥預(yù)報(bào)提供產(chǎn)量資料基礎(chǔ),以提高冬小麥產(chǎn)量預(yù)報(bào)的精細(xì)化和準(zhǔn)確率,滿足現(xiàn)代農(nóng)業(yè)氣象服務(wù)新需求。

    1 材料與方法

    1.1 研究區(qū)域及產(chǎn)量數(shù)據(jù)

    研究區(qū)域包括河南、山東、安徽、新疆等北方10個(gè)省、自治區(qū),根據(jù)2013?2017年統(tǒng)計(jì)產(chǎn)量數(shù)據(jù),該區(qū)域冬小麥5a平均總產(chǎn)量占全國(guó)總產(chǎn)量的90%以上,其中河南、山東、安徽、河北、江蘇5省約占82%,是最重要的小麥生產(chǎn)地。冬小麥單位面積實(shí)測(cè)產(chǎn)量數(shù)據(jù)為全國(guó)冬小麥農(nóng)業(yè)氣象觀測(cè)站1991?2017年的產(chǎn)量觀測(cè)數(shù)據(jù)(新疆和江蘇為1995?2017年,省級(jí)產(chǎn)量同),剔除觀測(cè)年份少于10a的站點(diǎn)后剩余123個(gè)站點(diǎn);各站點(diǎn)冬小麥實(shí)測(cè)產(chǎn)量數(shù)據(jù)來(lái)源于國(guó)家氣象信息中心。

    縣尺度單位面積產(chǎn)量資料包括兩類,一是直接使用123個(gè)站點(diǎn)的實(shí)測(cè)產(chǎn)量觀測(cè)數(shù)據(jù)作為站點(diǎn)所在縣的產(chǎn)量,稱為縣觀測(cè)產(chǎn)量;二是直接使用123個(gè)站點(diǎn)所在縣對(duì)應(yīng)的統(tǒng)計(jì)產(chǎn)量資料,稱為縣統(tǒng)計(jì)產(chǎn)量,該資料來(lái)源于各農(nóng)試站提供的當(dāng)年縣級(jí)統(tǒng)計(jì)產(chǎn)量。

    省級(jí)產(chǎn)量有三類,一是直接使用國(guó)家統(tǒng)計(jì)局公布的各省平均單產(chǎn),稱為省公布產(chǎn)量,來(lái)源于各省統(tǒng)計(jì)年鑒;二是由各縣產(chǎn)量集成計(jì)算,即以省內(nèi)涉及的各觀測(cè)站所在縣5a(2006?2010年)冬小麥種植面積平均值占比(占省內(nèi)所有觀測(cè)站所在縣種植面積之和)作為權(quán)重因子,對(duì)各站單產(chǎn)進(jìn)行加權(quán)集成得到所在省的冬小麥平均單產(chǎn),由各縣統(tǒng)計(jì)產(chǎn)量集成的省級(jí)產(chǎn)量稱為省統(tǒng)計(jì)產(chǎn)量,由各縣所轄農(nóng)試站觀測(cè)產(chǎn)量集成的稱為省觀測(cè)產(chǎn)量;各省觀測(cè)站站點(diǎn)信息及面積權(quán)重因子等資料見(jiàn)表1。

    表1 各省觀測(cè)站和縣數(shù)量及站點(diǎn)面積權(quán)重因子

    1.2 研究方法

    1.2.1 產(chǎn)量?jī)A向率

    對(duì)產(chǎn)量序列與年份進(jìn)行線性回歸分析[18?19],回歸系數(shù)b為產(chǎn)量線性傾向率,b>0表明產(chǎn)量隨時(shí)間呈增加趨勢(shì),b<0表明產(chǎn)量隨時(shí)間呈減少趨勢(shì),并對(duì)相關(guān)系數(shù)進(jìn)行顯著性檢驗(yàn)。

    1.2.2 產(chǎn)量系列變異系數(shù)

    用變異系數(shù)來(lái)反映縣觀測(cè)產(chǎn)量與縣統(tǒng)計(jì)產(chǎn)量各自間的離散程度[20],即

    式中,CVi為變異系數(shù),CVi越大表示產(chǎn)量越離散,SDi為縣觀測(cè)產(chǎn)量或縣統(tǒng)計(jì)產(chǎn)量的標(biāo)準(zhǔn)差,Meani為產(chǎn)量平均值,i為觀測(cè)產(chǎn)量或統(tǒng)計(jì)產(chǎn)量。

    2 結(jié)果與分析

    2.1 縣級(jí)尺度觀測(cè)產(chǎn)量與統(tǒng)計(jì)產(chǎn)量的差異

    2.1.1 平均值差異

    由圖1a可見(jiàn),1991?2000年10a間,河北、山東、江蘇、新疆、河南東部等地縣觀測(cè)產(chǎn)量平均值大部在4000kg·hm?2以上;山西、陜西、湖北、河南西部多在4000kg·hm?2以下;共有18個(gè)縣10a平均產(chǎn)量超過(guò)6000kg·hm?2,38個(gè)縣平均產(chǎn)量不足4000kg·hm?2,其中僅12個(gè)縣不足3000kg·hm?2。而從縣統(tǒng)計(jì)產(chǎn)量10a均值來(lái)看,山東、江蘇、新疆、河南東部多在4000kg·hm?2以上,僅7個(gè)縣10a平均產(chǎn)量超過(guò)6000kg·hm?2,55個(gè)縣統(tǒng)計(jì)產(chǎn)量平均值不足4000kg·hm?2,其中有26個(gè)不足3000kg·hm?2。從兩者差值來(lái)看,共有84個(gè)縣觀測(cè)產(chǎn)量高于縣統(tǒng)計(jì)產(chǎn)量,其中有36個(gè)縣超過(guò)1000kg·hm?2,尤以甘肅、陜西兩省突出;僅有39個(gè)縣觀測(cè)產(chǎn)量低于縣統(tǒng)計(jì)產(chǎn)量。

    由圖1b可見(jiàn),2001?2010年10a間,大部省份10a縣觀測(cè)產(chǎn)量平均值都在4000kg·hm?2以上,僅17個(gè)縣平均產(chǎn)量不足4000kg·hm?2,集中在山西、陜西和湖北;而有38個(gè)縣超過(guò)6000kg·hm?2。而從縣統(tǒng)計(jì)產(chǎn)量10a均值來(lái)看,共有42個(gè)縣統(tǒng)計(jì)產(chǎn)量不足4000kg·hm?2,且大部分集中在甘肅、陜西、山西、湖北等省,共有19個(gè)縣統(tǒng)計(jì)產(chǎn)量高于6000kg·hm?2。從兩者差值來(lái)看,共有91個(gè)縣觀測(cè)產(chǎn)量高于縣統(tǒng)計(jì)產(chǎn)量,其中有52個(gè)超過(guò)1000kg·hm?2;僅32個(gè)縣觀測(cè)產(chǎn)量低于縣統(tǒng)計(jì)產(chǎn)量。

    圖1 1991?2000年(a)、2001?2010年(b)和2011?2017(c)年縣尺度觀測(cè)產(chǎn)量(1)、統(tǒng)計(jì)產(chǎn)量(2)和兩者之差(3)

    由圖1c可見(jiàn),2011?2017年7a間,大部省份7a縣觀測(cè)產(chǎn)量平均值都在4000kg·hm?2以上,僅10個(gè)縣平均產(chǎn)量不足4000kg·hm?2,集中在山西、湖北;而有65個(gè)縣超過(guò)6000kg·hm?2,且大部分集中在河南、山東、河北、安徽及新疆。而從縣統(tǒng)計(jì)產(chǎn)量7a均值來(lái)看,共有22個(gè)縣統(tǒng)計(jì)產(chǎn)量不足4000kg·hm?2,45個(gè)縣統(tǒng)計(jì)產(chǎn)量高于6000kg·hm?2。從兩者差值來(lái)看,共有87個(gè)縣觀測(cè)產(chǎn)量高于統(tǒng)計(jì)產(chǎn)量,其中46個(gè)超過(guò)1000kg·hm?2;僅有36個(gè)縣觀測(cè)產(chǎn)量低于縣統(tǒng)計(jì)產(chǎn)量。

    可見(jiàn),對(duì)于縣觀測(cè)產(chǎn)量來(lái)說(shuō),超過(guò)6000kg·hm?2的縣個(gè)數(shù)從20世紀(jì)90年代的18個(gè)增至21世紀(jì)10年代的65個(gè),而不足4000kg·hm?2的個(gè)數(shù)從38個(gè)減至10個(gè)。對(duì)于縣統(tǒng)計(jì)產(chǎn)量來(lái)說(shuō),不足4000kg·hm?2的縣個(gè)數(shù)從20世紀(jì)90年代的55個(gè)減至21世紀(jì)00年代的42個(gè)再到10年代的22個(gè),而超過(guò)6000kg·hm?2的個(gè)數(shù)從7個(gè)增至45個(gè)。整體上看,縣觀測(cè)產(chǎn)量和縣統(tǒng)計(jì)產(chǎn)量均增幅明顯,21世紀(jì)10年代均為高產(chǎn)年代,但兩者差值在00年代達(dá)到峰值,2010年以后兩者差值減小,但仍高于20世紀(jì)90年代。

    2.1.2 變異系數(shù)差異

    對(duì)1991?2017年縣觀測(cè)產(chǎn)量變異系數(shù)的分析表明(圖2a),49個(gè)縣觀測(cè)產(chǎn)量變異系數(shù)小于0.20,46個(gè)變異系數(shù)在0.20~0.30,陜西、山西、甘肅等地共有8個(gè)站大于0.40。對(duì)縣統(tǒng)計(jì)產(chǎn)量變異系數(shù)的分析表明(圖2b),14個(gè)縣統(tǒng)計(jì)產(chǎn)量變異系數(shù)小于0.10,58個(gè)縣變異系數(shù)在0.10~0.20,僅9個(gè)站點(diǎn)變異系數(shù)大于0.30,表明統(tǒng)計(jì)產(chǎn)量年際間差異較小。從兩者對(duì)比來(lái)看,新疆和山東兩省區(qū)內(nèi)所有縣變異系數(shù)都小于0.30,觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量穩(wěn)定性最高;僅甘肅通渭、陜西永壽和河北黃驊三地觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量變異系數(shù)均大于0.3。

    整體上看,縣尺度變異系數(shù)觀測(cè)產(chǎn)量離散程度要高于縣統(tǒng)計(jì)產(chǎn)量,因?yàn)榭h觀測(cè)產(chǎn)量是農(nóng)試站地塊實(shí)際測(cè)得的產(chǎn)量,代表的范圍有限,且易受到管理措施和氣象災(zāi)害的影響;而縣統(tǒng)計(jì)產(chǎn)量則是全縣的平均單產(chǎn),代表全縣的產(chǎn)量水平,且當(dāng)某地出現(xiàn)災(zāi)害導(dǎo)致減產(chǎn)時(shí),全縣其余地塊可能出現(xiàn)不受災(zāi)或者增產(chǎn)的可能性,縣內(nèi)不同地塊存在補(bǔ)償性,這就導(dǎo)致縣統(tǒng)計(jì)產(chǎn)量相對(duì)穩(wěn)定,所以觀測(cè)產(chǎn)量波動(dòng)大,統(tǒng)計(jì)產(chǎn)量波動(dòng)較小。一般情況在作物品種和種植方式較為穩(wěn)定以及農(nóng)業(yè)措施較為完善的情形下,縣單產(chǎn)的年際間變化可以認(rèn)為是由于氣象條件波動(dòng)造成的,一定程度上也為利用觀測(cè)產(chǎn)量進(jìn)行產(chǎn)量預(yù)報(bào)提供了客觀條件。

    2.1.3 線性傾向率差異

    由圖3可見(jiàn),1991?2017年(圖3a),各縣觀測(cè)產(chǎn)量中有73個(gè)站點(diǎn)線性變化呈現(xiàn)顯著增加趨勢(shì),大多集中在河北、河南、山東、江蘇、安徽等主產(chǎn)??;僅湖北麻城、新疆莎車、甘肅文縣和山東福山4個(gè)縣呈現(xiàn)顯著減少趨勢(shì),其余46個(gè)站點(diǎn)線性變化不顯著;傾向率高于200kg·hm?2·a?1的5個(gè)站全部通過(guò)顯著性檢驗(yàn),表明產(chǎn)量增幅明顯。而各站所在縣的統(tǒng)計(jì)產(chǎn)量(圖3b)有100個(gè)縣線性變化呈現(xiàn)顯著增加趨勢(shì);僅湖北隨州和谷城、山東淄博、新疆若羌和山西臨猗5個(gè)站點(diǎn)呈現(xiàn)顯著減少趨勢(shì),僅18個(gè)站點(diǎn)線性變化趨勢(shì)不顯著。總共72個(gè)縣觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量的線性變化趨勢(shì)均通過(guò)顯著性檢驗(yàn),其中63個(gè)縣都呈現(xiàn)顯著增加趨勢(shì),僅6個(gè)縣線性變化趨勢(shì)相反。縣觀測(cè)產(chǎn)量?jī)A向率在?164~266kg·hm?2·a?1,波動(dòng)幅度較大;而縣統(tǒng)計(jì)產(chǎn)量?jī)A向率波動(dòng)較小,產(chǎn)量增幅小于縣觀測(cè)產(chǎn)量,且符合產(chǎn)量逐年增加的趨勢(shì)。

    圖2 1991?2017年縣觀測(cè)產(chǎn)量(a)和統(tǒng)計(jì)產(chǎn)量(b)的變異系數(shù)

    注:畫黑色圈的站點(diǎn)通過(guò)0.05水平顯著性檢驗(yàn)。

    Note: The stations with black cycle mean P<0.05.

    2.2 省級(jí)尺度產(chǎn)量的差異

    2.2.1 各年代平均值的差異

    結(jié)合站點(diǎn)面積權(quán)重計(jì)算各省級(jí)年尺度觀測(cè)產(chǎn)量與統(tǒng)計(jì)產(chǎn)量,在此基礎(chǔ)上計(jì)算各省每10a平均值,由表2可以看出,對(duì)省觀測(cè)產(chǎn)量來(lái)說(shuō),除新疆和山西外,其余8省10a均值呈現(xiàn)增加的趨勢(shì),2011?2017年7a間平均值最高,山東和安徽增幅超過(guò)2000kg·hm?2;新疆為2001?2010年10a平均值最高,山西為最低。對(duì)省統(tǒng)計(jì)產(chǎn)量來(lái)說(shuō),各省10a平均值均為2011?2017年7a最高,河南和安徽增幅超過(guò)2000kg·hm?2,湖北和山西在2001?2010年10a平均值最低。對(duì)于省公布產(chǎn)量,除了湖北省2001?2010年10a均值最低外,其余各省均為2011?2017年7a間平均值最高,1991?2000年10a平均值最小,呈現(xiàn)增產(chǎn)的趨勢(shì)。

    從省觀測(cè)產(chǎn)量與省統(tǒng)計(jì)產(chǎn)量?jī)烧卟钪祦?lái)看,1991?2000年10a間,僅湖北省觀測(cè)產(chǎn)量低于省統(tǒng)計(jì)產(chǎn)量,山東省兩者最接近,為203kg·hm?2,甘肅差值最高,達(dá)2000kg·hm?2;2001?2010年10a間,河南、新疆、湖北三省兩者最為接近,其差值不足25kg·hm?2,但河南為負(fù)值,甘肅差值最高,達(dá)2367kg·hm?2;2011?2017年7a間,安徽兩者差值最低,為99kg·hm?2,河南和新疆觀測(cè)產(chǎn)量小于統(tǒng)計(jì)產(chǎn)量。從省觀測(cè)產(chǎn)量與省公布產(chǎn)量?jī)烧卟钪祦?lái)看,除新疆外,其余各省7a觀測(cè)產(chǎn)量均高于省公布產(chǎn)量,尤其是安徽、江蘇、陜西三省各個(gè)年代差值均高于1000kg·hm?2,甘肅高于2000kg·hm?2;河南在2000年以后省觀測(cè)產(chǎn)量與公布產(chǎn)量較為接近。

    總體來(lái)看,與站點(diǎn)尺度每10a均值一致,2001?2010年為各省觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量的高產(chǎn)年代,1991?2000年為低產(chǎn)年代;山東、安徽、河北、江蘇、陜西和山西等省10a觀測(cè)產(chǎn)量均值在各個(gè)年代均高于省統(tǒng)計(jì)產(chǎn)量均值;除山西省觀測(cè)產(chǎn)量與統(tǒng)計(jì)產(chǎn)量差值各年代維持在1000kg·hm?2以外,其余省份各年代間差值變化各不相同。省公布產(chǎn)量則呈現(xiàn)年代際增產(chǎn)的趨勢(shì),這也與農(nóng)業(yè)技術(shù)管理水平的提升保持一致;但總體上低于省觀測(cè)產(chǎn)量和省統(tǒng)計(jì)產(chǎn)量。

    表2 省級(jí)尺度觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量以及公布產(chǎn)量每10a均值(kg·hm?2)

    注:Aveobs、Avestaand Aveann分別為省觀測(cè)產(chǎn)量、省統(tǒng)計(jì)產(chǎn)量和省公布產(chǎn)量每10a平均值。

    Notes: Aveobs, Avestaand Aveannmean average of observed yield, statistical yield and announced yield.

    2.2.2 省級(jí)變異系數(shù)的差異

    從各省所有縣觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量的標(biāo)準(zhǔn)差和變異系數(shù)的平均值來(lái)看(表3),各省內(nèi)所有縣觀測(cè)產(chǎn)量標(biāo)準(zhǔn)差的平均值都在1000kg·hm?2以上,安徽省標(biāo)準(zhǔn)差的平均值最大;各省縣統(tǒng)計(jì)產(chǎn)量標(biāo)準(zhǔn)差的平均值在571~1285kg·hm?2,新疆最小,安徽最大;且觀測(cè)產(chǎn)量標(biāo)準(zhǔn)差平均值都高于統(tǒng)計(jì)產(chǎn)量標(biāo)準(zhǔn)差平均值。從兩者變異系數(shù)平均值來(lái)看,山東、江蘇、新疆三省區(qū)觀測(cè)產(chǎn)量變異系數(shù)平均值均小于0.20;而山東、江蘇、新疆、湖北四省區(qū)統(tǒng)計(jì)產(chǎn)量變異系數(shù)平均值小于0.20。除河北和甘肅兩省變異系數(shù)平均值相同外,其余各省觀測(cè)產(chǎn)量變異系數(shù)均高于統(tǒng)計(jì)產(chǎn)量變異系數(shù)平均值??傮w來(lái)看,省內(nèi)平均變異系數(shù)觀測(cè)產(chǎn)量離散程度高于統(tǒng)計(jì)產(chǎn)量,這也與縣尺度變異系數(shù)變化結(jié)果相吻合。

    表3 各省內(nèi)縣觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量標(biāo)準(zhǔn)差以及變異系數(shù)平均值對(duì)比

    注:SDobs和SDsta分別為省觀測(cè)產(chǎn)量和省統(tǒng)計(jì)產(chǎn)量平均值(kg·hm?2);CVobs和CVsta分別為省觀測(cè)產(chǎn)量和省統(tǒng)計(jì)產(chǎn)量變異系數(shù)。

    Notes:SDobsand SDstamean average of standard deviation of observed yield and statistical yield(kg·ha?1); CVobsand CVstamean average of coefficient of variation of observed yield and statistical yield.

    2.2.3 省級(jí)尺度產(chǎn)量線性傾向率的差異

    由省級(jí)尺度不同產(chǎn)量?jī)A向率可知(表4),除新疆和山西省觀測(cè)產(chǎn)量?jī)A向率未通過(guò)顯著性檢驗(yàn)外,其余各省傾向率都通過(guò)顯著性檢驗(yàn),且新疆傾向率為負(fù)值;從傾向率數(shù)值來(lái)看,安徽、山東、河北、河南各省觀測(cè)產(chǎn)量增幅明顯。省統(tǒng)計(jì)產(chǎn)量和公布產(chǎn)量?jī)A向率均通過(guò)顯著性檢驗(yàn),且產(chǎn)量增幅均為正值;山東、河北、湖北、陜西各省兩者傾向率都低于觀測(cè)產(chǎn)量的傾向率;僅安徽省觀測(cè)、統(tǒng)計(jì)和公布產(chǎn)量?jī)A向率均高于100kg·hm?2·a?1。新疆因地域遼闊,各地小麥生產(chǎn)氣候條件差異很大,導(dǎo)致觀測(cè)產(chǎn)量?jī)A向率不顯著;山西由于冬小麥主產(chǎn)區(qū)位于汾河中下游且一般都有灌溉條件,導(dǎo)致觀測(cè)產(chǎn)量相對(duì)穩(wěn)定變化不明顯,因此觀測(cè)產(chǎn)量?jī)A向率不顯著。省公布產(chǎn)量由于考慮全省的冬小麥生產(chǎn)情況,受氣象條件影響相對(duì)較小,往往隨著品種改良及農(nóng)業(yè)技術(shù)水平提升而呈現(xiàn)增產(chǎn)的趨勢(shì)。

    表4 各省觀測(cè)產(chǎn)量、統(tǒng)計(jì)產(chǎn)量和公布產(chǎn)量?jī)A向率

    注:Tenobs、Tensta和Tenann分別為省觀測(cè)產(chǎn)量、統(tǒng)計(jì)產(chǎn)量和公布產(chǎn)量?jī)A向率,*表示通過(guò)0.05水平的顯著性檢驗(yàn)。

    Notes: Tenobs, Tenstaand Tenannmean the tendency ratio of observed yield, statistical yield and announced yield.*is P<0.05.

    3 結(jié)論與討論

    3.1 結(jié)論

    (1)縣觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量均表現(xiàn)出高產(chǎn)縣(單產(chǎn)高于6000kg·hm?2)數(shù)量增多、低產(chǎn)縣(單產(chǎn)低于4000kg·hm?2)數(shù)量減少的趨勢(shì);21世紀(jì)10年代兩者產(chǎn)量均為高產(chǎn)年代。

    (2)縣觀測(cè)產(chǎn)量變異系數(shù)離散程度高于縣統(tǒng)計(jì)產(chǎn)量,縣觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量變異系數(shù)小于0.20的個(gè)數(shù)分別為49和72;新疆和山東所有縣兩種產(chǎn)量變異系數(shù)都小于0.30,產(chǎn)量穩(wěn)定性最高??h觀測(cè)產(chǎn)量和統(tǒng)計(jì)產(chǎn)量分別有73和100個(gè)站點(diǎn)傾向率呈現(xiàn)顯著增加趨勢(shì),主要集中在河北、河南、山東、江蘇、安徽等冬小麥主產(chǎn)省。

    (3)21世紀(jì)10年代為省觀測(cè)產(chǎn)量和省統(tǒng)計(jì)產(chǎn)量的高產(chǎn)年代,山東、安徽、河北、江蘇、陜西和山西六省10a省觀測(cè)產(chǎn)量的平均值在所有年代均高于省統(tǒng)計(jì)產(chǎn)量平均值。

    (4)除河北和甘肅外其余各省觀測(cè)產(chǎn)量變異系數(shù)平均值均高于統(tǒng)計(jì)產(chǎn)量的平均值,除新疆和山西外其余各省觀測(cè)產(chǎn)量?jī)A向率均通過(guò)顯著性檢驗(yàn),所有省統(tǒng)計(jì)和公布產(chǎn)量?jī)A向率均通過(guò)顯著性檢驗(yàn),且產(chǎn)量增幅均為正值。

    3.2 討論

    縣尺度作物產(chǎn)量受到農(nóng)業(yè)技術(shù)措施[11,21?22]、生育期氣象條件變化[23?24]、統(tǒng)計(jì)因素等影響導(dǎo)致波動(dòng),而統(tǒng)計(jì)產(chǎn)量變異性小于觀測(cè)產(chǎn)量的主要原因在于:(1)縣觀測(cè)點(diǎn)以外地區(qū)的災(zāi)害有無(wú)及發(fā)生類型與特點(diǎn)不同,全縣大災(zāi)年中也有不受災(zāi)或?yàn)?zāi)輕的增產(chǎn)地塊,豐收年也有少數(shù)受災(zāi)減產(chǎn)地塊,存在相互補(bǔ)償;而觀測(cè)地段的災(zāi)害發(fā)生較為單一。(2)不同農(nóng)田對(duì)單一災(zāi)害大面積發(fā)生時(shí)的響應(yīng)不同,例如干旱年低濕地增產(chǎn),多雨年高崗地增產(chǎn);秋涼年早播有利,秋暖年晚播有利;不同地塊之間存在相互補(bǔ)償導(dǎo)致縣統(tǒng)計(jì)產(chǎn)量較為穩(wěn)定,而觀測(cè)地段不存在相互補(bǔ)償現(xiàn)象。(3)部分地區(qū)出于政策原因往往在減產(chǎn)年少報(bào)減產(chǎn),而增產(chǎn)年少報(bào)增產(chǎn)。

    利用站點(diǎn)觀測(cè)資料進(jìn)行省級(jí)集成時(shí),僅選擇了觀測(cè)站點(diǎn)所在縣的近5a冬小麥種植面積平均值作為面積權(quán)重,由于各地冬小麥種植面積的不斷變化,加之冬小麥地段觀測(cè)站點(diǎn)個(gè)數(shù)遠(yuǎn)少于主產(chǎn)縣個(gè)數(shù),用觀測(cè)站點(diǎn)面積進(jìn)行加權(quán)集成也有一定的局限性;如何在沒(méi)有冬小麥全部縣級(jí)作物種植面積的情形下,進(jìn)行產(chǎn)量的升尺度集成也是值得探討的話題。當(dāng)前利用衛(wèi)星遙感技術(shù)提取作物種植面積相關(guān)技術(shù)也較為成熟[25?28],國(guó)家級(jí)農(nóng)業(yè)氣象業(yè)務(wù)中冬小麥遙感提取最小尺度僅做到地級(jí)市一級(jí),且針對(duì)的是北方冬小麥主產(chǎn)區(qū),要做到縣尺度提取面積尚有困難,在現(xiàn)有條件下,利用觀測(cè)站所在縣5a種植面積加權(quán)集成到省級(jí)尺度相對(duì)容易實(shí)現(xiàn)。

    由于部分農(nóng)業(yè)氣象觀測(cè)站存在遷站、觀測(cè)作物種類改變等客觀原因,以及農(nóng)業(yè)氣象觀測(cè)站的種植管理水平、災(zāi)害發(fā)生情況不同于大田,加之最終測(cè)產(chǎn)中存在的人工統(tǒng)計(jì)誤差等因素,會(huì)導(dǎo)致觀測(cè)產(chǎn)量與統(tǒng)計(jì)產(chǎn)量之間存在一定的差異,因此,基于觀測(cè)產(chǎn)量資料的產(chǎn)量預(yù)報(bào)應(yīng)該主要針對(duì)氣象產(chǎn)量的增減趨勢(shì)開(kāi)展預(yù)報(bào)服務(wù),判別當(dāng)年氣象條件對(duì)作物產(chǎn)量的影響是增產(chǎn)還是減產(chǎn),利用氣象產(chǎn)量的增減幅度,結(jié)合統(tǒng)計(jì)模型預(yù)報(bào)最終產(chǎn)量。

    站點(diǎn)觀測(cè)產(chǎn)量最直接反應(yīng)了當(dāng)?shù)貧庀髼l件的利弊影響,因此數(shù)據(jù)更具代表性和客觀性;同時(shí),目前產(chǎn)量預(yù)報(bào)應(yīng)用的產(chǎn)量資料唯一來(lái)源為統(tǒng)計(jì)局公布產(chǎn)量,但公布產(chǎn)量獲取存在較長(zhǎng)的滯后性,特別是縣級(jí)尺度的產(chǎn)量序列,獲取難度更大。因此,利用各站點(diǎn)觀測(cè)產(chǎn)量建立氣象部門的產(chǎn)量資料序列,并在此基礎(chǔ)上建立縣級(jí)和省級(jí)作物產(chǎn)量預(yù)報(bào)模型,可以補(bǔ)充和驗(yàn)證利用公布數(shù)據(jù)建立的產(chǎn)量預(yù)報(bào)模型結(jié)論,同時(shí)氣象部門應(yīng)用新數(shù)據(jù)源開(kāi)展的作物產(chǎn)量預(yù)報(bào)結(jié)論也可為統(tǒng)計(jì)和農(nóng)業(yè)部門提供參考。

    [1] 徐敏,徐經(jīng)緯,高萍,等.不同統(tǒng)計(jì)模型在冬小麥產(chǎn)量預(yù)報(bào)中的預(yù)報(bào)能力評(píng)估:以江蘇麥區(qū)為例[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2020,28(3): 438-447.

    Xu M,Xu J W,Gao P,et al.Evaluation of winter wheat yield prediction ability of different statistical models:a case study of Jiangsu wheat-growing region[J].Chinese Journal of Eco-Agriculture,2020,28(3): 438-447.(in Chinese)

    [2] 侯英雨,張蕾,吳門新,等.國(guó)家級(jí)現(xiàn)代農(nóng)業(yè)氣象業(yè)務(wù)技術(shù)進(jìn)展[J].應(yīng)用氣象學(xué)報(bào),2018,29(6):641-656.

    Hou Y Y,Zhang L,Wu M X,et al.Advances of modern agrometeorological service and technology in China[J]. Journal of Applied Meteorological Science,2018,29(6): 641-656.(in Chinese)

    [3] 吳門新,莊立偉,侯英雨,等.中國(guó)農(nóng)業(yè)氣象業(yè)務(wù)系統(tǒng)(CAgMSS)設(shè)計(jì)與實(shí)現(xiàn)[J].應(yīng)用氣象學(xué)報(bào),2019,30(5): 513-527.

    Wu M X,Zhuang L W,Hou Y Y,et al.The design and implementation of China agricultural meteorological service system(CAgMSS)[J].Journal of Applied Meteorological Science,2019,30(5):513-527.(in Chinese)

    [4] 李秀芬,馬樹(shù)慶,宮麗娟,等.基于WOFOST的東北地區(qū)玉米生育期氣象條件適宜度評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)氣象,2013, 34(1):43-49.

    Li X F,Ma S Q,Gong L J,et al.Evaluation of meteorological suitability degree during maize growth period based on WOFOST in northeast China[J].Chinese Journal of Agrometeorology,2013,34(1):43-49.(in Chinese)

    [5] 李秀芬,馬樹(shù)慶,趙慧穎,等.基于WOFOST模型的內(nèi)蒙古河套灌區(qū)玉米低溫冷害評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)氣象,2016, 37(3):352-360.

    Li X F,Ma S Q,Zhao H Y,et al.Evaluation on maize chilling injury based on WOFOST model in Hetao irrigation region in Inner Mongolia[J].Chinese Journal of Agrometeorology, 2016,37(3):352-360.(in Chinese)

    [6] 張陽(yáng),王連喜,李琪,等.基于WOFOST 模型的吉林省中西部春玉米灌溉模擬[J].中國(guó)農(nóng)業(yè)氣象,2018,39(6):411-420.

    Zhang Y,Wang L X,Li Q,et al.Irrigation simulation of spring maize in central and western of Jilin province based on WOFOST model[J].Chinese Journal of Agrometeorology, 2018,39(6):411-420.(in Chinese)

    [7] 欒青,郭建平,馬雅麗,等.基于線性生長(zhǎng)假設(shè)的作物積溫模型穩(wěn)定性比較[J].中國(guó)農(nóng)業(yè)氣象,2020,41(11):695-706.

    Luan Q,Guo J P,Ma Y L, et al.Comparison of model’s stability about integrated temperature based on linear hypotheses[J].Chinese Journal of Agrometeorology,2020, 41(11):695-706.(in Chinese)

    [8] 張軼,劉布春,楊曉娟,等.基于農(nóng)作物災(zāi)情的長(zhǎng)江中下游地區(qū)糧食產(chǎn)量損失評(píng)估[J].中國(guó)農(nóng)業(yè)氣象,2018,39(4): 280-291.

    Zhang Y,Liu B C,Yang X J,et al.Grain yield loss evaluation based on agro-meteorological disaster exposure in the Middle-Lower Yangtze Plain[J].Chinese Journal of Agrometeorology,2018,39(4):280-291.(in Chinese)

    [9] 萬(wàn)能涵,楊曉光,劉志娟,等.氣候變化背景下中國(guó)主要作物農(nóng)業(yè)氣象災(zāi)害時(shí)空分布特征Ⅲ:華北地區(qū)夏玉米干旱[J].中國(guó)農(nóng)業(yè)氣象,2018,39(4):209-219.

    Wan N H,Yang X G,Liu Z J,et al.Temporal and spatial variations of agro-meteorological disasters of main crops in China in a changing climate(Ⅲ):drought of summer maize in North China Plain[J].Chinese Journal of Agrometeorology, 2018,39(4):209-219.(in Chinese)

    [10] 譚凱炎,鄔定榮,趙花榮.氣候變暖背景下華北平原冬小麥生育期溫度條件變化趨勢(shì)分析[J].中國(guó)農(nóng)業(yè)氣象,2017, 38(6):333-341.

    Tan K Y,Wu D R,Zhao H R.Trend analysis of temperature conditions over different growth periods of winter wheat under climate warming in North China Plain[J].Chinese Journal of Agrometeorology,2017,38(6):333-341.(in Chinese)

    [11] 李勤英,姚鳳梅,張佳華,等.不同農(nóng)藝措施對(duì)縮小冬小麥產(chǎn)量差和提高氮肥利用率的評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)氣象, 2018,39(6):370-379.

    Li Q Y,Yao F M,Zhang J H,et al.Evaluation of different agronomic measures on narrowing the yield gap and improving nitrogen use efficiency of winter wheat[J]. Chinese Journal of Agrometeorology,2018,39(6):370-379. (in Chinese)

    [12] 馬雪晴,胡琦,潘學(xué)標(biāo),等.1961-2015年華北平原夏玉米生長(zhǎng)季氣候年型及其影響分析[J].中國(guó)農(nóng)業(yè)氣象,2019, 40(2):65-75.

    Ma X Q,Hu Q,Pan X B,et al.Analysis of annual climate types and its impact on summer maize in the North China Plain over the period 1961-2015[J].Chinese Journal of Agrometeorology,2019,40(2):65-75.(in Chinese)

    [13] 劉思華,李晶,黃晚華,等.湖南省春玉米光溫生產(chǎn)潛力和產(chǎn)量差特征及限制因子分析[J].中國(guó)農(nóng)業(yè)氣象,2020, 41(2):94-101.

    Liu S H,Li J,Huang W H,et al.Characteristics and limiting factors of light-temperature potential productivity and yield gap of spring maize in Hunan province[J].Chinese Journal of Agrometeorology,2020,41(2):94-101.(in Chinese)

    [14] 吳冰潔,王靖,唐建昭,等.華北平原冬小麥產(chǎn)量變異的氣象影響因子分析[J].中國(guó)農(nóng)業(yè)氣象,2018,39(10):623-635.

    Wu B J,Wang J,Tang J Z,et al.Meteorological influencing factors on variation in winter wheat yield in the North China Plain[J].Chinese Journal of Agrometeorology,2018, 39(10):623-635.(in Chinese)

    [15] 趙凱娜,寧曉菊,秦耀辰,等.縣域冬小麥生育期氣候要素與產(chǎn)量的灰色關(guān)聯(lián)度分析:以河南省為例[J].中國(guó)農(nóng)業(yè)氣象,2017,38(11):729-737.

    Zhao K N,Ning X J,Qin Y C,et al.Grey correlation analysis between climatic elements during winter wheat growing period and yield in counties:a case study of Henan Province[J].Chinese Journal of Agrometeorology,2017, 38(11):729-737.(in Chinese)

    [16] Asseng S,Ewert F,Rosenzweig C,et al.Uncertainty in simulating wheat yields under climate change[J].Nature Climate Change,20133(9):827-832.

    [17] Tao F,Zhang Z,Xiao D,et al.Responses of wheat growth and yield to climate change in different climate zones of China,1981-2009[J].Agricultural and Forest Meteorology, 2014,189:91-104.

    [18] 崔日鮮.山東省太陽(yáng)總輻射的時(shí)空變化特征分析[J].自然資源學(xué)報(bào),2014,29(10):1780-1791.

    Cui R X.The analysis of spatiotemporal variation characteristics of global solar radiation in Shandong province[J].Journal of Natural Resources,2014,29(10): 1780-1791.(in Chinese)

    [19] 王秀芬,楊艷昭,尤飛.近30年來(lái)黑龍江省氣候變化特征分析[J].中國(guó)農(nóng)業(yè)氣象, 2011,32(S1):28-32.

    Wang X F,Yang Y Z,You F.Analysis on characteristics of climate change in recent 30 years in Heilongjiang province[J].Chinese Journal of Agrometeorology,2011, 32(S1):28-32.(in Chinese)

    [20] 錢永蘭,毛留喜,周廣勝.全球主要糧食作物產(chǎn)量變化及其氣象災(zāi)害風(fēng)險(xiǎn)評(píng)估[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(1):226-235.

    Qian Y L,Mao L X,Zhou G S.Changes in global main crop yields and its meteorological risk assessment[J]. Transactions of the CSAE,2016,32(1):226-235.(in Chinese)

    [21] Wang F,He Z,Sayre K,et al.Wheat cropping systems and technologies in China[J].Field Crops Research,2009,111(3): 181-188.

    [22] Zhang P,Wei T,Wang H X,et al.Effects of straw mulch on soil water and winter wheat production in dryland farming[J].Scientific Reports,2015(5):10725.

    [23] 肖登攀,陶福祿,沈彥俊,等.華北平原冬小麥對(duì)過(guò)去30年氣候變化響應(yīng)的敏感性研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2014,22(4):430-438.

    Xiao D P,Tao F L,Shen Y J,et al.Sensitivity of response of winter wheat to climate change in the North China Plain in the last three decades[J].Chinese Journal of Eco-Agriculture, 2014,22(4):430-438.(in Chinese)

    [24] 郝立生,吳雁,王榮英.海河低平原春季氣候變化對(duì)冬小麥產(chǎn)量的影響[J].氣象與減災(zāi)研究,2007,30(4):20-24.

    Hao L S,Wu Y,Wang R Y.Impact of spring climate change on wheat yield in lower Haihe Plain[J].Meteorological and Disaster Reduction Research,2007,30(4):20-24.(in Chinese)

    [25] 錢永蘭,楊邦杰,焦險(xiǎn)峰.基于遙感抽樣的國(guó)家尺度農(nóng)作物面積統(tǒng)計(jì)方法評(píng)估[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(11):180- 187.

    Qian Y L,Yang B J,Jiao X F.Accuracy assessment on the crop area estimating method based on RS sampling at national scale[J].Transact ions of the CSAE,2007,23(11): 180-187.(in Chinese)

    [26] 陳鵬飛,楊飛,杜佳.基于環(huán)境減災(zāi)衛(wèi)星時(shí)序歸一化植被指數(shù)的冬小麥產(chǎn)量估測(cè)[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(11): 124-131.

    Chen P F,Yang F,Du J.Yield forecasting for winter wheat using time series NDVI from HJ satellite[J].Transactions of the CSAE,2013,29(11):124-131.(in Chinese)

    [27] 黃健熙,侯矞焯,蘇偉,等.基于GF-1 WFV 數(shù)據(jù)的玉米與大豆種植面積提取方法[J].農(nóng)業(yè)工程學(xué)報(bào),2017,33(7): 164-170.

    Huang J X,Hou Y Z,Su W,et al.Mapping corn and soybean cropped area with GF-1 WFV data[J]. Transactions of the CSAE,2017,33(7):164-170.(in Chinese)

    [28] 張超,童亮,劉哲.基于多時(shí)相GF-1 WFV 和高分紋理的制種玉米田識(shí)別[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,50(2):163-168.

    Zhang C,Tong L,Liu Z,et al.Identification method of seed maize plot based on multi-temporal GF-1 WFV and Kompsat-3 texture[J].Transactions of the Chinese Society for Agricultural Machinery,2019,50(2):163-168.(in Chinese)

    Studies on the Difference of Observed Yield and Statistical Yield of Winter Wheat

    LIU Wei1, MENG Cui-li2, SONG Ying-bo1

    (1. National Meteorological Center, Beijing 100081, China; 2. Argo-meteorological Station of Wuhan, Wuhan 430040)

    The difference of interdecadal variations, coefficient of variation and tendency ratio between the observed yield of winter wheat from 123 agrometeorological observation stations and the statistical yield of winter wheat at county level where the observation station was located from 1991 to 2017. The proportion of average winter wheat planting area in each county in five years(2006?2010) was used as the weight factor to integrate the observed yield and statistical yield at province level, at the same time using the announced yield at province level from National Bureau of Statistics. The interdecadal variations and tendency ration of three different yields at provincial level were compared and analyzed. The results showed that:(1) the number of high yield counties increased significantly, and low yield counties decreased significantly in both observed yield and statistical yield counties. The two yield were both high yield years in the 2010s, and the difference between the two reached peak value in the 2000s. (2) The coefficient of variation of observed yield at the county scale was higher than the statistical yield. The coefficient of variation of statistical yield in 49 counties were less than 0.20 and only 8 were greater than 0.40, while 72 statistical yield counties were less than 0.20 and only 9 were greater than 0.30. The coefficient of variation of statistical yield in all counties in Xinjiang and Shandong provinces were less than 0.30. (3) The tendency ratio of 73 observed yield counties showed a significant increase mostly concentrated in the major producing provinces such as Hebei, Henan, Shandong, Jiangsu, and Anhui; and 100 statistical yield counties showed the same significant increase. The tendency ratio of observed and statistical yield in 72 counties passed the significance test at the same time. (4) The 2000s were the high yield years for both observed and statistical yield at provinces level and 1990s were the low yield years. The average of the observed yield in every 10 years was higher than the average of the statistical yield in Shandong, Anhui, Hebei, Jiangsu, Shaanxi and Shanxi province. (5) Eight provinces had passed the significant test on tendency ration of the observed yield at the provincial level except for Xinjiang and Shanxi province. While tendency ration of the statistical yield and the announced yield in all provinces had passed the significant test and the yield growth was positive. In general, the winter wheat yield series based on the observed yield could provide a new data source for yield forecast.

    Winter wheat; Agrometeorological observation stations; Observed yield; Statistical yield ; Area weight factor integration

    10.3969/j.issn.1000-6362.2021.02.004

    劉維,孟翠麗,宋迎波.冬小麥觀測(cè)產(chǎn)量與統(tǒng)計(jì)產(chǎn)量的差異性分析[J].中國(guó)農(nóng)業(yè)氣象,2021,42(2):123-133

    2020?07?07

    國(guó)家氣象中心預(yù)報(bào)員專項(xiàng)(Y201912);2019年國(guó)內(nèi)外作物產(chǎn)量預(yù)報(bào)專項(xiàng);2020年國(guó)內(nèi)外作物產(chǎn)量預(yù)報(bào)專項(xiàng)

    宋迎波,研究員,研究方向?yàn)樽魑锂a(chǎn)量預(yù)報(bào),E-mail:songyb@cma.gov.cn

    劉維,E-mail:rainvswindvs@163.com

    猜你喜歡
    冬小麥平均值站點(diǎn)
    “平均值代換”法在數(shù)學(xué)解題中的應(yīng)用
    基于Web站點(diǎn)的SQL注入分析與防范
    電子制作(2019年14期)2019-08-20 05:43:42
    2017~2018年冬季西北地區(qū)某站點(diǎn)流感流行特征分析
    首屆歐洲自行車共享站點(diǎn)協(xié)商會(huì)召開(kāi)
    甘肅冬小麥田
    怕被人認(rèn)出
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    平面圖形中構(gòu)造調(diào)和平均值幾例
    基于電流平均值的改進(jìn)無(wú)功檢測(cè)法
    冬小麥——新冬18號(hào)
    av免费在线看不卡| 如何舔出高潮| 免费少妇av软件| 日日摸夜夜添夜夜添av毛片| 18+在线观看网站| 欧美成人精品欧美一级黄| 最近2019中文字幕mv第一页| 国产无遮挡羞羞视频在线观看| 国产亚洲av片在线观看秒播厂| 亚洲人成网站在线观看播放| 五月开心婷婷网| 你懂的网址亚洲精品在线观看| 91午夜精品亚洲一区二区三区| 一区二区三区精品91| 夫妻午夜视频| 国产精品久久久久成人av| 伊人久久国产一区二区| 欧美+日韩+精品| 亚洲av中文av极速乱| 日韩av免费高清视频| 精品卡一卡二卡四卡免费| 亚洲国产最新在线播放| 国产欧美亚洲国产| 久久鲁丝午夜福利片| 免费大片黄手机在线观看| 老司机影院毛片| videosex国产| 人体艺术视频欧美日本| 九九在线视频观看精品| 国产精品国产三级专区第一集| 欧美日韩视频精品一区| 国产欧美日韩综合在线一区二区| 免费观看a级毛片全部| 一区在线观看完整版| 国产一区二区在线观看日韩| 亚洲av.av天堂| 亚洲av福利一区| 美女xxoo啪啪120秒动态图| 春色校园在线视频观看| 99热全是精品| 自线自在国产av| 精品亚洲乱码少妇综合久久| 在线亚洲精品国产二区图片欧美 | a 毛片基地| 精品国产国语对白av| 一级黄片播放器| 久久这里有精品视频免费| 我的老师免费观看完整版| 爱豆传媒免费全集在线观看| 99久久精品国产国产毛片| 欧美成人午夜免费资源| 国产午夜精品一二区理论片| 人妻制服诱惑在线中文字幕| 亚洲国产日韩一区二区| 男女免费视频国产| 视频在线观看一区二区三区| 亚洲欧美成人综合另类久久久| 免费播放大片免费观看视频在线观看| 色婷婷av一区二区三区视频| 国产一区亚洲一区在线观看| 9色porny在线观看| 嫩草影院入口| 亚洲av免费高清在线观看| 少妇人妻精品综合一区二区| 亚洲美女搞黄在线观看| 九草在线视频观看| 精品国产乱码久久久久久小说| 久久精品久久精品一区二区三区| 亚洲国产欧美日韩在线播放| 国产欧美亚洲国产| 又大又黄又爽视频免费| 黄色怎么调成土黄色| 老司机亚洲免费影院| 另类精品久久| 亚洲av综合色区一区| 丝瓜视频免费看黄片| 久久久国产精品麻豆| 国产伦精品一区二区三区视频9| 亚洲成人手机| 亚洲国产精品国产精品| 亚洲成人av在线免费| 女人精品久久久久毛片| 午夜老司机福利剧场| 国产高清三级在线| 大话2 男鬼变身卡| 97超视频在线观看视频| 91久久精品国产一区二区成人| 色婷婷久久久亚洲欧美| 欧美精品亚洲一区二区| a级毛色黄片| 免费少妇av软件| 天天操日日干夜夜撸| 夫妻午夜视频| 三级国产精品片| 欧美人与善性xxx| av播播在线观看一区| 精品久久国产蜜桃| 天美传媒精品一区二区| av卡一久久| 国产黄色视频一区二区在线观看| 日韩一区二区视频免费看| 制服丝袜香蕉在线| 亚洲欧洲国产日韩| 一边亲一边摸免费视频| 国产综合精华液| 美女内射精品一级片tv| av在线app专区| 午夜福利在线观看免费完整高清在| 国产欧美另类精品又又久久亚洲欧美| 国产精品久久久久久精品电影小说| 日韩精品免费视频一区二区三区 | 国产黄色视频一区二区在线观看| 亚洲av不卡在线观看| 亚洲欧美色中文字幕在线| 亚洲经典国产精华液单| 日韩av不卡免费在线播放| 亚洲av成人精品一二三区| 日韩三级伦理在线观看| 国产日韩一区二区三区精品不卡 | 老司机亚洲免费影院| 亚洲精品成人av观看孕妇| 中文字幕免费在线视频6| 久久精品熟女亚洲av麻豆精品| 22中文网久久字幕| 99九九线精品视频在线观看视频| 精品久久蜜臀av无| 精品人妻熟女av久视频| 国产精品不卡视频一区二区| 日韩伦理黄色片| 国产在视频线精品| 精品视频人人做人人爽| 欧美亚洲 丝袜 人妻 在线| 天美传媒精品一区二区| 久久精品国产鲁丝片午夜精品| 精品熟女少妇av免费看| 多毛熟女@视频| a级毛片黄视频| 特大巨黑吊av在线直播| 久久国产精品男人的天堂亚洲 | 免费少妇av软件| 男的添女的下面高潮视频| 搡老乐熟女国产| 欧美成人精品欧美一级黄| 久久国产精品大桥未久av| av又黄又爽大尺度在线免费看| 狂野欧美白嫩少妇大欣赏| 伦精品一区二区三区| 草草在线视频免费看| 热99国产精品久久久久久7| 成年女人在线观看亚洲视频| 亚洲精品,欧美精品| 91精品伊人久久大香线蕉| 天天影视国产精品| 免费人妻精品一区二区三区视频| 91精品国产国语对白视频| 一个人免费看片子| 日韩av在线免费看完整版不卡| 黑人欧美特级aaaaaa片| 丝瓜视频免费看黄片| 国产成人精品无人区| 制服人妻中文乱码| 亚洲不卡免费看| 成人毛片60女人毛片免费| 99九九线精品视频在线观看视频| 丰满饥渴人妻一区二区三| 日韩欧美一区视频在线观看| 成人毛片a级毛片在线播放| 亚洲国产精品专区欧美| 18禁在线播放成人免费| 成年人午夜在线观看视频| 高清毛片免费看| av.在线天堂| 男女高潮啪啪啪动态图| 午夜精品国产一区二区电影| 免费高清在线观看日韩| 久久久久久久久大av| 午夜福利视频精品| 国产 精品1| 久久久a久久爽久久v久久| 精品少妇黑人巨大在线播放| 七月丁香在线播放| 亚洲欧美一区二区三区国产| 国产成人精品福利久久| 国产欧美另类精品又又久久亚洲欧美| 日韩 亚洲 欧美在线| 80岁老熟妇乱子伦牲交| 亚洲国产色片| 简卡轻食公司| 成人二区视频| 色5月婷婷丁香| 亚洲精品日韩av片在线观看| 亚洲欧美日韩卡通动漫| 久久久久久人妻| 大香蕉久久网| 国产一区二区三区av在线| 少妇的逼好多水| a 毛片基地| 美女中出高潮动态图| 纵有疾风起免费观看全集完整版| 亚洲综合色网址| 简卡轻食公司| 久久这里有精品视频免费| 22中文网久久字幕| 伦理电影免费视频| 少妇人妻精品综合一区二区| 18+在线观看网站| 91国产中文字幕| 午夜免费鲁丝| 欧美日韩亚洲高清精品| 国产成人精品福利久久| 一区二区av电影网| 亚洲婷婷狠狠爱综合网| 日日啪夜夜爽| 成人亚洲欧美一区二区av| 五月玫瑰六月丁香| a级毛片黄视频| av在线app专区| 五月开心婷婷网| 少妇丰满av| 亚洲婷婷狠狠爱综合网| 丝瓜视频免费看黄片| 成人免费观看视频高清| freevideosex欧美| 一级a做视频免费观看| 99久国产av精品国产电影| videos熟女内射| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区| 国产日韩欧美亚洲二区| 欧美日韩在线观看h| 五月天丁香电影| 久久久久国产网址| 多毛熟女@视频| 亚洲精品一区蜜桃| 一区二区三区精品91| 久久国内精品自在自线图片| 久久久久国产网址| 亚洲精品国产色婷婷电影| 成年av动漫网址| 亚洲精品亚洲一区二区| 夜夜骑夜夜射夜夜干| 精品少妇黑人巨大在线播放| 国产精品免费大片| 亚洲少妇的诱惑av| 亚洲美女黄色视频免费看| 2021少妇久久久久久久久久久| 国产成人freesex在线| 一级片'在线观看视频| 国产欧美日韩一区二区三区在线 | 亚洲色图综合在线观看| 赤兔流量卡办理| 美女国产视频在线观看| 另类精品久久| 母亲3免费完整高清在线观看 | 亚洲一级一片aⅴ在线观看| 国产一级毛片在线| 国产精品 国内视频| 一级毛片我不卡| 在线精品无人区一区二区三| 高清av免费在线| 伊人久久国产一区二区| 哪个播放器可以免费观看大片| a 毛片基地| 九九在线视频观看精品| 亚洲精品视频女| 成人免费观看视频高清| 日韩av不卡免费在线播放| 国产精品国产三级国产av玫瑰| 久久午夜福利片| 狂野欧美白嫩少妇大欣赏| 日韩电影二区| 全区人妻精品视频| 亚洲综合色惰| 青春草亚洲视频在线观看| 久久精品夜色国产| 熟女电影av网| 啦啦啦在线观看免费高清www| 国产视频内射| 母亲3免费完整高清在线观看 | 美女脱内裤让男人舔精品视频| 五月伊人婷婷丁香| 狠狠婷婷综合久久久久久88av| 亚洲av成人精品一二三区| 国产黄片视频在线免费观看| 麻豆精品久久久久久蜜桃| 亚洲精品一二三| 人人妻人人澡人人看| 男男h啪啪无遮挡| 精品亚洲成a人片在线观看| 我的老师免费观看完整版| 国产日韩一区二区三区精品不卡 | 天天躁夜夜躁狠狠久久av| 性色avwww在线观看| 日本黄色日本黄色录像| 亚洲国产精品一区三区| 在线观看一区二区三区激情| 日本与韩国留学比较| 国产精品欧美亚洲77777| 两个人的视频大全免费| 精品久久蜜臀av无| 精品一区二区三区视频在线| 午夜免费观看性视频| 黄片无遮挡物在线观看| 国产精品久久久久久av不卡| 中文字幕久久专区| 黄片无遮挡物在线观看| 美女内射精品一级片tv| 青春草视频在线免费观看| 2022亚洲国产成人精品| 久久综合国产亚洲精品| 丝袜在线中文字幕| 国产在线免费精品| 亚洲国产精品专区欧美| 欧美成人精品欧美一级黄| 久久热精品热| 精品少妇内射三级| 99久久精品一区二区三区| 国产黄频视频在线观看| 中文字幕av电影在线播放| av线在线观看网站| 夜夜看夜夜爽夜夜摸| 日本av手机在线免费观看| 久久毛片免费看一区二区三区| 久久人人爽人人爽人人片va| 99久久人妻综合| 亚洲精品久久久久久婷婷小说| 18禁在线播放成人免费| 亚洲美女搞黄在线观看| 永久网站在线| 国产综合精华液| av在线app专区| 一级黄片播放器| 国产乱人偷精品视频| 色5月婷婷丁香| 街头女战士在线观看网站| 91精品伊人久久大香线蕉| 999精品在线视频| 美女中出高潮动态图| 99热这里只有精品一区| 久久久久久久久久人人人人人人| 97超视频在线观看视频| 欧美 亚洲 国产 日韩一| 成人综合一区亚洲| 母亲3免费完整高清在线观看 | 边亲边吃奶的免费视频| 只有这里有精品99| 免费观看无遮挡的男女| 国产高清有码在线观看视频| 国产成人av激情在线播放 | 久久精品国产鲁丝片午夜精品| 性高湖久久久久久久久免费观看| 晚上一个人看的免费电影| 午夜福利,免费看| 天天躁夜夜躁狠狠久久av| 亚洲内射少妇av| 一级黄片播放器| 人人妻人人澡人人爽人人夜夜| 天美传媒精品一区二区| 伊人久久精品亚洲午夜| 男女边吃奶边做爰视频| 大香蕉97超碰在线| 国产成人av激情在线播放 | av在线app专区| 国产黄色免费在线视频| 黄色毛片三级朝国网站| 亚洲性久久影院| 久久影院123| a级片在线免费高清观看视频| √禁漫天堂资源中文www| 日韩一区二区视频免费看| 色哟哟·www| 欧美一级a爱片免费观看看| 日韩中字成人| 国产精品久久久久久精品电影小说| 人妻一区二区av| 日韩制服骚丝袜av| 亚洲国产精品成人久久小说| 国产精品.久久久| a级毛片免费高清观看在线播放| 考比视频在线观看| av天堂久久9| 日日摸夜夜添夜夜添av毛片| 欧美日韩视频高清一区二区三区二| 伦精品一区二区三区| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 美女国产视频在线观看| 麻豆乱淫一区二区| 成人手机av| videossex国产| 夫妻性生交免费视频一级片| 最后的刺客免费高清国语| 性高湖久久久久久久久免费观看| 丰满迷人的少妇在线观看| videos熟女内射| 国产精品成人在线| 在现免费观看毛片| 九草在线视频观看| 色哟哟·www| 午夜精品国产一区二区电影| freevideosex欧美| 美女cb高潮喷水在线观看| 午夜视频国产福利| 看非洲黑人一级黄片| 秋霞在线观看毛片| 全区人妻精品视频| 久久99热这里只频精品6学生| 91成人精品电影| 亚洲熟女精品中文字幕| 最近2019中文字幕mv第一页| 久久久午夜欧美精品| 在线观看www视频免费| 国产 一区精品| 纯流量卡能插随身wifi吗| 亚洲美女视频黄频| 又黄又爽又刺激的免费视频.| 久久人人爽人人爽人人片va| 久久毛片免费看一区二区三区| 亚洲国产精品一区三区| 国产欧美亚洲国产| 啦啦啦视频在线资源免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费大片18禁| 美女中出高潮动态图| 午夜激情久久久久久久| 国产又色又爽无遮挡免| 国产成人av激情在线播放 | 99久久中文字幕三级久久日本| 精品久久国产蜜桃| 国产精品免费大片| 国产黄色视频一区二区在线观看| 人妻一区二区av| 婷婷色综合大香蕉| av国产久精品久网站免费入址| a级片在线免费高清观看视频| 久久精品国产自在天天线| 97精品久久久久久久久久精品| 亚洲国产精品一区二区三区在线| 9色porny在线观看| 中文字幕久久专区| 女人精品久久久久毛片| 一区二区三区精品91| 人妻人人澡人人爽人人| 亚洲久久久国产精品| 欧美精品国产亚洲| 熟妇人妻不卡中文字幕| 日产精品乱码卡一卡2卡三| 国产亚洲精品久久久com| 伦理电影大哥的女人| 丰满迷人的少妇在线观看| 国产欧美亚洲国产| 高清欧美精品videossex| 久久热精品热| 久久精品久久久久久久性| 视频在线观看一区二区三区| 人妻少妇偷人精品九色| 精品熟女少妇av免费看| 久久99热这里只频精品6学生| h视频一区二区三区| 国产一区二区在线观看av| 99热这里只有精品一区| 黑人高潮一二区| 久久毛片免费看一区二区三区| 视频中文字幕在线观看| 精品国产一区二区三区久久久樱花| 亚洲国产精品成人久久小说| 91aial.com中文字幕在线观看| 国产亚洲精品久久久com| 国产免费一级a男人的天堂| 人妻少妇偷人精品九色| 国产免费一级a男人的天堂| 欧美 亚洲 国产 日韩一| 久久久国产精品麻豆| 青春草亚洲视频在线观看| 男女无遮挡免费网站观看| 在线观看三级黄色| 国产伦理片在线播放av一区| 丝袜美足系列| 伊人亚洲综合成人网| 免费黄色在线免费观看| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区国产| 制服丝袜香蕉在线| 精品少妇内射三级| 国产一区二区三区综合在线观看 | 在线观看三级黄色| 日本欧美视频一区| 精品久久蜜臀av无| 日韩视频在线欧美| 大码成人一级视频| 国产日韩欧美亚洲二区| 最黄视频免费看| 亚洲av综合色区一区| 国产成人精品久久久久久| 亚洲激情五月婷婷啪啪| videossex国产| 精品酒店卫生间| 久久久久久久精品精品| 97超碰精品成人国产| 亚洲精品国产色婷婷电影| 国产精品人妻久久久久久| 黑人欧美特级aaaaaa片| 国产成人freesex在线| 国产黄片视频在线免费观看| 岛国毛片在线播放| 精品卡一卡二卡四卡免费| 高清午夜精品一区二区三区| 国产精品熟女久久久久浪| 国产精品 国内视频| 人妻系列 视频| 日本av免费视频播放| 女的被弄到高潮叫床怎么办| 又大又黄又爽视频免费| av在线app专区| 一本—道久久a久久精品蜜桃钙片| 国产日韩欧美在线精品| 色哟哟·www| 中国美白少妇内射xxxbb| 青春草亚洲视频在线观看| 国产在线视频一区二区| 午夜福利视频精品| 国产男女超爽视频在线观看| 精品国产国语对白av| 女性生殖器流出的白浆| 91精品伊人久久大香线蕉| 插逼视频在线观看| av电影中文网址| 伦理电影大哥的女人| 亚洲国产精品一区三区| 亚洲精华国产精华液的使用体验| 亚洲欧洲日产国产| 另类亚洲欧美激情| 久久久久人妻精品一区果冻| 国产淫语在线视频| 日日撸夜夜添| 国产在线一区二区三区精| 国产熟女午夜一区二区三区 | 伦理电影大哥的女人| 国产免费视频播放在线视频| 肉色欧美久久久久久久蜜桃| 国产无遮挡羞羞视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 色5月婷婷丁香| 欧美 亚洲 国产 日韩一| 午夜av观看不卡| 黄色欧美视频在线观看| 国精品久久久久久国模美| 国产成人精品无人区| 三上悠亚av全集在线观看| 国产精品三级大全| 欧美激情国产日韩精品一区| 国产免费视频播放在线视频| 高清午夜精品一区二区三区| 人妻人人澡人人爽人人| 人人妻人人澡人人爽人人夜夜| 18禁在线无遮挡免费观看视频| 国产免费现黄频在线看| 国产免费视频播放在线视频| 人人妻人人澡人人爽人人夜夜| av电影中文网址| 亚洲欧美精品自产自拍| 曰老女人黄片| 少妇人妻久久综合中文| 欧美97在线视频| 亚洲精品色激情综合| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产av新网站| 久久久久久久久久久丰满| 26uuu在线亚洲综合色| 丝袜在线中文字幕| 欧美日韩在线观看h| 91精品国产九色| 国产 一区精品| 国产欧美日韩一区二区三区在线 | 狂野欧美白嫩少妇大欣赏| 这个男人来自地球电影免费观看 | 国产白丝娇喘喷水9色精品| 99热6这里只有精品| 日韩大片免费观看网站| av专区在线播放| 天堂8中文在线网| 亚洲经典国产精华液单| 人妻人人澡人人爽人人| 母亲3免费完整高清在线观看 | videosex国产| 国产av国产精品国产| 国产精品一区二区在线观看99| 老司机影院成人| 黄片无遮挡物在线观看| 中文字幕人妻熟人妻熟丝袜美| 欧美国产精品一级二级三级| 人妻一区二区av| 美女内射精品一级片tv| 两个人的视频大全免费| 最后的刺客免费高清国语| 亚洲精品成人av观看孕妇| 22中文网久久字幕| 狂野欧美白嫩少妇大欣赏| 制服诱惑二区| 欧美精品高潮呻吟av久久| 99九九线精品视频在线观看视频| 黑丝袜美女国产一区| 日韩 亚洲 欧美在线| 国产午夜精品久久久久久一区二区三区| 99热这里只有精品一区| 97超视频在线观看视频| 男人添女人高潮全过程视频| 18禁在线播放成人免费| 亚洲经典国产精华液单| 久久久久网色| 欧美精品一区二区免费开放| 精品亚洲乱码少妇综合久久| 国产男女内射视频| 亚洲久久久国产精品| av电影中文网址|