華北平原冬小麥產(chǎn)量變異的氣象影響因子分析

吳冰潔，王 靖，唐建昭，王 娜, 3，徐 琳，白慧卿，鄭雋卿，王 娜，李 揚(yáng)

?

華北平原冬小麥產(chǎn)量變異的氣象影響因子分析

吳冰潔1, 2，王 靖1**，唐建昭1，王 娜1, 3，徐 琳1，白慧卿1，鄭雋卿1，王 娜1，李 揚(yáng)1

（1. 中國(guó)農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，北京 100193；2. 成都錦天聯(lián)華科技有限責(zé)任公司，成都 610072；3. 山東省農(nóng)業(yè)科學(xué)院作物研究所，濟(jì)南 250100）

采用1988?2015年華北平原冬小麥種植區(qū)46個(gè)市的統(tǒng)計(jì)產(chǎn)量和相應(yīng)46個(gè)氣象站點(diǎn)的逐日氣象資料，通過(guò)Logistic曲線(xiàn)和雙曲線(xiàn)方法分離出氣象產(chǎn)量，并構(gòu)建氣象產(chǎn)量與生長(zhǎng)季主要?dú)庀笠蜃拥亩嘣y(tǒng)計(jì)關(guān)系，以明確華北平原冬小麥產(chǎn)量變異的氣象影響因子。結(jié)果表明：（1）1988?2015年華北平原冬小麥產(chǎn)量在3200～6800kg·hm?2，中部地區(qū)產(chǎn)量最高，南部地區(qū)產(chǎn)量的變異高于中部和北部地區(qū)。（2）生長(zhǎng)季日照時(shí)數(shù)、溫度和降水平均值的年際變化影響了17%～78%的氣象產(chǎn)量的變異，其中54%的地區(qū)達(dá)到顯著水平（P＜0.05）。影響程度較高的地區(qū)主要分布在河北南部、山東西部和河南的東北部地區(qū)。（3）播種?返青階段的降水顯著影響產(chǎn)量變異，降水量每增加1%，天津、駐馬店及山東西北部等地產(chǎn)量將上升13～74kg·hm?2，而河北北部、河南南部、山東南部等地產(chǎn)量將下降16～80kg·hm?2。返青?成熟階段對(duì)產(chǎn)量變異影響較大的因子為最低氣溫，平均最低氣溫每上升1?C，天津和石家莊、山東東部和西部及河南東部等地產(chǎn)量將增加50～295kg·hm?2，而北京、唐山和棗莊等地將減少76～124kg·hm?2。總體來(lái)看，溫度對(duì)華北平原冬小麥產(chǎn)量變異影響范圍更廣且更加顯著，但氣象因子對(duì)產(chǎn)量變異的影響受局地品種和管理措施等影響呈現(xiàn)較大的空間差異。

產(chǎn)量變率；氣象產(chǎn)量；日照時(shí)數(shù)；溫度；降水

華北平原是中國(guó)重要的小麥產(chǎn)區(qū)，小麥總產(chǎn)和種植面積占全國(guó)小麥總產(chǎn)和種植面積的50%左右[1]，提高該區(qū)域小麥總產(chǎn)對(duì)保障中國(guó)糧食安全具有重要意義。增加糧食總產(chǎn)的主要途徑，一是擴(kuò)大糧食播種面積，二是提高糧食單產(chǎn)水平[2]。然而中國(guó)人均耕地面積少，后備耕地資源緊缺，種植面積進(jìn)一步擴(kuò)大的可能性低[3]，因此，提升糧食單產(chǎn)是增加區(qū)域糧食總產(chǎn)的有效途徑。小麥單產(chǎn)的提高主要集中在品種和栽培措施的改進(jìn)[4?5]以及對(duì)氣候變化的趨利避害方面[6?7]。分析氣象因子對(duì)小麥產(chǎn)量形成的影響對(duì)適應(yīng)氣候變化和防災(zāi)減災(zāi)具有重要意義。Fang等[8]通過(guò)人工增溫試驗(yàn)，表明小麥產(chǎn)量隨著溫度升高而增加。鄭娜等[9]發(fā)現(xiàn)生長(zhǎng)季較高的光照和溫度日較差對(duì)冬小麥生長(zhǎng)發(fā)育和產(chǎn)量形成有利，特別是越冬前較高的日照時(shí)數(shù)和溫度日較差對(duì)冬小麥產(chǎn)量有顯著的正影響。閆素輝等[10?11]認(rèn)為生育中后期弱光脅迫對(duì)小麥產(chǎn)量形成有不利影響，弱光環(huán)境影響植株蒸騰速率和氣孔導(dǎo)度，導(dǎo)致光合產(chǎn)物合成受抑制、小麥灌漿速率下降，從而顯著降低小麥產(chǎn)量。譚凱炎等[12]指出拔節(jié)?成熟階段溫度過(guò)高不利于穗花發(fā)育、籽粒形成和灌漿，尤其是花后高溫顯著降低冬小麥產(chǎn)量，因?yàn)闇囟壬咄鶗?huì)伴隨著干熱風(fēng)和高溫逼熟等災(zāi)害天氣的影響而導(dǎo)致小麥減產(chǎn)[13]。生育前期降水增加一般有利于小麥產(chǎn)量提高，如返青?拔節(jié)期降水有利于小麥的生長(zhǎng)發(fā)育，拔節(jié)?抽穗是冬小麥生長(zhǎng)需水盛期，降水量過(guò)少會(huì)對(duì)產(chǎn)量產(chǎn)生不利影響[14]，而后期尤其是開(kāi)花期和成熟期降水過(guò)多會(huì)導(dǎo)致小麥減產(chǎn)[15]。前人研究表明冬小麥干旱高風(fēng)險(xiǎn)區(qū)主要分布在河北省環(huán)渤海地區(qū)及山東省的長(zhǎng)島[16?17]。王春乙等[17]發(fā)現(xiàn)河南大部分地區(qū)冬小麥生育后期會(huì)受到干旱的影響。以上研究表明，氣象因子對(duì)作物產(chǎn)量的影響隨地區(qū)和冬小麥生育階段不同而存在顯著差異，就華北平原冬小麥而言，應(yīng)針對(duì)其不同的生育階段展開(kāi)研究。

當(dāng)前對(duì)作物產(chǎn)量和氣象因子的關(guān)系研究主要基于站點(diǎn)或者省級(jí)和國(guó)家尺度，如Lobell等[18?19]在全球尺度上的研究指出，1980年以來(lái)的溫度升高導(dǎo)致全球小麥產(chǎn)量降低。You等[20]評(píng)估了氣象因子對(duì)中國(guó)小麥產(chǎn)量的影響，發(fā)現(xiàn)小麥生長(zhǎng)季溫度每升高1℃，產(chǎn)量下降3%～10%。Tao等[21]基于省級(jí)尺度研究表明，在中國(guó)小麥的主產(chǎn)省份中，黑龍江、遼寧和新疆生長(zhǎng)季溫度日較差與小麥產(chǎn)量呈顯著負(fù)相關(guān)，但與廣東的小麥產(chǎn)量呈正相關(guān)。但在市級(jí)尺度上對(duì)產(chǎn)量和氣象因子關(guān)系鮮有研究報(bào)道，而中國(guó)詳細(xì)完整的小麥產(chǎn)量統(tǒng)計(jì)年鑒資料的基本單元為市級(jí)尺度，因此，研究市級(jí)尺度產(chǎn)量和氣候的關(guān)系對(duì)采取措施應(yīng)對(duì)氣候變化和抵御不利氣象災(zāi)害有重要的意義。本研究以華北平原冬小麥生產(chǎn)區(qū)市級(jí)小麥產(chǎn)量為研究對(duì)象，分析不同區(qū)域影響小麥產(chǎn)量變異的氣象因子，以期為提高華北平原冬小麥的氣候適應(yīng)性和抵御農(nóng)業(yè)氣象災(zāi)害提供科學(xué)依據(jù)。

1 資料與方法

1.1 研究區(qū)域和數(shù)據(jù)來(lái)源

華北平原主要包括河北省、河南省、山東省、北京市和天津市，該區(qū)屬半濕潤(rùn)暖溫帶，冬小麥生長(zhǎng)季總輻射為1920～2760MJ·m?2，總?cè)照諘r(shí)數(shù)為1400～2000h，總降水量為100～300mm。土壤類(lèi)型主要為潮土和褐土，土層深厚，適宜耕作，種植制度以冬小麥?夏玉米輪作為主。

華北平原46個(gè)市的1988?2015年逐年冬小麥單產(chǎn)數(shù)據(jù)來(lái)源于中國(guó)種植業(yè)信息網(wǎng)和《河北農(nóng)村統(tǒng)計(jì)年鑒》、《河北經(jīng)濟(jì)統(tǒng)計(jì)年鑒》、《河南農(nóng)村統(tǒng)計(jì)年鑒》、《河南統(tǒng)計(jì)年鑒》、《山東統(tǒng)計(jì)年鑒》、《北京統(tǒng)計(jì)年鑒》和《天津統(tǒng)計(jì)年鑒》。其中河北省各市統(tǒng)計(jì)產(chǎn)量時(shí)間序列為1993?2015年，山東省日照市和萊蕪市、河南省濟(jì)源市的產(chǎn)量統(tǒng)計(jì)年份分別是1990?2015年、1993?2015年和1996?2015年，其余市的產(chǎn)量數(shù)據(jù)年份均為1988?2015年。

華北平原各市對(duì)應(yīng)的氣象站點(diǎn)數(shù)據(jù)來(lái)源于中國(guó)氣象局國(guó)家氣象信息中心地面氣象觀(guān)測(cè)資料（圖1），包括逐日日照時(shí)數(shù)（h）、日平均氣溫（?C）、日最高氣溫（?C）、日最低氣溫（?C）和日降水量（mm）。個(gè)別站點(diǎn)缺少日照時(shí)數(shù)觀(guān)測(cè)數(shù)據(jù)，則使用鄰近站點(diǎn)日照時(shí)數(shù)代替。將華北平原冬小麥生長(zhǎng)季劃分為播種、越冬、返青、拔節(jié)、開(kāi)花和成熟6個(gè)發(fā)育階段，不同種植區(qū)域冬小麥各生育階段長(zhǎng)度根據(jù)1992? 2013年華北平原各種植區(qū)域內(nèi)43個(gè)農(nóng)業(yè)氣象觀(guān)測(cè)站實(shí)際觀(guān)測(cè)值的平均值確定。

圖1 華北平原冬小麥種植區(qū)和46個(gè)氣象站點(diǎn)分布

1.2 研究方法

1.2.1 氣象產(chǎn)量分離

糧食生產(chǎn)受自然環(huán)境因子和社會(huì)生產(chǎn)力水平的綜合影響，在長(zhǎng)時(shí)間序列的作物產(chǎn)量與氣象因子關(guān)系的統(tǒng)計(jì)研究中，通常將實(shí)際糧食單產(chǎn)（Y）分解為趨勢(shì)產(chǎn)量（Yt）、氣象產(chǎn)量（Yw）和隨機(jī)誤差產(chǎn)量（ε），即

受冬小麥品種、田間管理措施、地理環(huán)境和地方經(jīng)濟(jì)變革等因素影響，各地冬小麥實(shí)際產(chǎn)量隨時(shí)間變化的趨勢(shì)會(huì)有所不同，因此，分別采用Logistic曲線(xiàn)和雙曲線(xiàn)法擬合華北平原各市冬小麥趨勢(shì)產(chǎn)量，進(jìn)而從實(shí)際產(chǎn)量中分離出氣象產(chǎn)量并與氣象因子進(jìn)行回歸，以回歸方程解釋率（決定系數(shù)R2）高為標(biāo)準(zhǔn)選擇相應(yīng)的趨勢(shì)產(chǎn)量擬合方法，以提高氣象產(chǎn)量分離的準(zhǔn)確性和合理性。

1.2.2 Logistic曲線(xiàn)

Logistic曲線(xiàn)方程[22]的形式為

式中，a、b為待定參數(shù)，k為實(shí)際產(chǎn)量的最大值。對(duì)式（3）兩邊取自然對(duì)數(shù)可得

1.2.3 雙曲線(xiàn)

雙曲線(xiàn)（Hyperbola）形式為

式中，a、b為模型擬合參數(shù)。

表1列出了最終選擇的各市氣象產(chǎn)量分離方法。

1.2.4 產(chǎn)量?氣象因子關(guān)系模型

采用多元線(xiàn)性回歸分析確定冬小麥氣象產(chǎn)量與生育期主要?dú)庀笠蜃樱ㄗ罡邷囟?、最低溫度、平均溫度、溫度日較差、日照時(shí)數(shù)和降水）平均值以及各階段（播種?越冬、越冬?返青、返青?拔節(jié)、拔節(jié)?開(kāi)花、開(kāi)花?成熟）氣象因子平均值的關(guān)系，即

式中，xj為第j個(gè)氣象因子的值，y為冬小麥氣象產(chǎn)量，a0是回歸常數(shù)，aj是偏回歸系數(shù)，反映在其它所有自變量不變的情況下，某一自變量單位變化對(duì)產(chǎn)量的影響程度。

表2列出了研究區(qū)各省冬小麥生育階段起始時(shí)間。

表1 各市氣象產(chǎn)量的分離方法

表2 各省冬小麥的生育階段日期

Note: E- is the first ten-day of a month; M- is the middle ten-day of a month; L- is the last ten-day of a month.

1.2.5 氣象因子對(duì)氣象產(chǎn)量的決定系數(shù)（R2）

式中，yi為第i年的氣象產(chǎn)量，yavg為氣象產(chǎn)量的n年平均值，xi為第i年的氣象因子，xavg為氣象因子的n年平均值，n=1，…，28。決定系數(shù)R2取值越接近1，表明氣象因子對(duì)氣象產(chǎn)量變異的解釋率越高，對(duì)氣象產(chǎn)量變異的影響越大。

1.3 數(shù)據(jù)處理

用各站實(shí)際產(chǎn)量和各氣象因子時(shí)間序列資料進(jìn)行線(xiàn)性回歸計(jì)算，得到相應(yīng)的變化傾向率和變異系數(shù)。

產(chǎn)量?jī)A向率y的計(jì)算式為

式中，t為年份，a×10 為傾向率。

變異系數(shù)（CV）的計(jì)算式為

利用ArcGIS10.2軟件表達(dá)研究因子的時(shí)空分布特征，SPSS19.0進(jìn)行統(tǒng)計(jì)分析。

2 結(jié)果與分析

2.1 華北平原冬小麥產(chǎn)量的時(shí)空變異分析

由圖2a可見(jiàn)，華北平原冬小麥實(shí)際產(chǎn)量在3200～6800kg·hm?2，平均為5200kg·hm?2。高產(chǎn)區(qū)主要分布在河北、山東以及河南的東北部地區(qū)，低產(chǎn)區(qū)分布在河南西部山區(qū)和東南部。由圖2b可見(jiàn)，各市冬小麥產(chǎn)量線(xiàn)性變化的傾向率在?601～1957kg·hm?2·10a?1，96%的市為顯著增加（P＜0.05），僅北京和山東萊蕪兩個(gè)地區(qū)為顯著下降（P＜0.05）。產(chǎn)量增加較多的區(qū)域主要分布在河北南部和河南東部地區(qū)，其中河南商丘、周口和駐馬店三市產(chǎn)量增加最顯著，每10a增加1500kg·hm?2以上。圖2c顯示，1988?2015年華北平原冬小麥實(shí)際產(chǎn)量年際變異系數(shù)在7%～32%，其中南部小麥產(chǎn)量變異最大。綜合來(lái)看，華北北部地區(qū)平均產(chǎn)量較高，增速較慢，年際變異較低，而華北南部平均產(chǎn)量較低，增速較快，年際變異較高。

圖2 1988?2015年華北平原冬小麥實(shí)際產(chǎn)量平均值（a）、傾向率（b）和變異系數(shù)（c）的空間分布

2.2 生長(zhǎng)季氣象因子對(duì)冬小麥氣象產(chǎn)量變異的貢獻(xiàn)分析

圖3表明，華北平原冬小麥生長(zhǎng)季主要?dú)庀笠蜃悠骄档淖兓軌蚪忉寶庀螽a(chǎn)量變異的17%～78%，其中54%的市達(dá)到顯著水平（P＜0.05），生長(zhǎng)季主要?dú)庀笠蜃拥钠骄底儺悓?duì)氣象產(chǎn)量變異解釋程度較高的區(qū)域主要集中在河北的邢臺(tái)和邯鄲，山東西南部的聊城、濟(jì)寧和泰安以及河南東北部的開(kāi)封、商丘等地區(qū)。從單個(gè)氣象因子來(lái)看，圖4a顯示生長(zhǎng)季平均氣溫對(duì)氣象產(chǎn)量變異的解釋率在0～56%，26%的市達(dá)到顯著水平（P＜0.05），主要集中在河北北部和河南東部。圖4b表明生長(zhǎng)季平均最高氣溫對(duì)氣象產(chǎn)量變異解釋率在0～23%，15%的地區(qū)達(dá)到顯著水平（P＜0.05），主要集中在河南北部和河北北部。圖4c顯示生長(zhǎng)季平均最低氣溫對(duì)氣象產(chǎn)量變異的解釋率在0～67%，47%的市達(dá)到顯著水平（P＜0.05），集中在河南北部和河北北部。圖4d顯示生長(zhǎng)季降水對(duì)氣象產(chǎn)量變異的解釋率為0～24%，僅15%的地區(qū)達(dá)到顯著水平（P＜0.05），解釋率較高的區(qū)域主要分布在天津和河南北部。圖4e顯示生長(zhǎng)季日照時(shí)數(shù)對(duì)氣象產(chǎn)量變異的解釋率在0～51%，22%的地區(qū)達(dá)到顯著水平（P＜0.05），解釋率較高的區(qū)域分布在山東南部、河南東北部和河北西南部。圖4f表明溫度日較差對(duì)氣象產(chǎn)量變異的解釋率在0～52%，37%的市達(dá)到顯著水平（P＜0.05），主要分布于河北保定、山東東南部以及河南東部。

圖3 冬小麥氣象產(chǎn)量與生長(zhǎng)季氣象因子平均值間多元線(xiàn)性回歸方程的決定系數(shù)（1988?2015年）

2.3 不同生育階段氣象因子對(duì)冬小麥氣象產(chǎn)量變異的貢獻(xiàn)分析

由表3和表4可知，華北平原冬小麥在播種－越冬階段對(duì)氣象產(chǎn)量變異影響最大的氣象因子主要為降水和氣溫日較差。該生育階段內(nèi)降水每增加1%，天津、河北秦皇島、河南駐馬店及山東日照、煙臺(tái)和濟(jì)南等地小麥氣象產(chǎn)量增加13～72kg·hm?2，而北京、河北廊坊，山東濱州，河南南部的濟(jì)源、漯河、周口和信陽(yáng)小麥氣象產(chǎn)量降低23～80kg·hm?2。該階段氣溫日較差每升高1?C，天津、河北保定、山東菏澤和濮陽(yáng)、河南漯河和三門(mén)峽氣象產(chǎn)量增加94～279kg·hm?2，而河北滄州和邢臺(tái)、山東濟(jì)南和濟(jì)寧及河南南陽(yáng)和北部地區(qū)氣象產(chǎn)量減少54～173kg·hm?2。表4顯示，在冬小麥越冬?返青階段對(duì)氣象產(chǎn)量變異影響最主要的因子是降水量，降水量每增加1%，北京、河北邯鄲、河南鶴壁、濟(jì)源和信陽(yáng)及山東中部和南部等地氣象產(chǎn)量減少16～62kg·hm?2；由表3、表4可得，在返青?拔節(jié)階段影響氣象產(chǎn)量變異的主要因子為最低氣溫和日照時(shí)數(shù)，該階段平均最低氣溫每上升1?C，天津、河北石家莊、河南濟(jì)源、商丘和駐馬店及山東濮陽(yáng)、日照和青島氣象產(chǎn)量增加66～254kg·hm?2；該階段平均日照時(shí)數(shù)每增加1h，河北邢臺(tái)及山東淄博冬小麥氣象產(chǎn)量增加108～158kg·hm?2，而河北邯鄲、山東濱州和煙臺(tái)及河南西南部地區(qū)等地氣象產(chǎn)量減少108～215kg·hm?2。表3表明在拔節(jié)?開(kāi)花階段對(duì)氣象產(chǎn)量變異影響最大的因子主要為最低氣溫和氣溫日較差，該階段平均最低氣溫每升高1?C，河北石家莊及河南焦作、南陽(yáng)和周口氣象產(chǎn)量增加50～142kg·hm?2，而河北唐山和滄州及山東棗莊氣象產(chǎn)量減少76～124kg·hm?2；該階段氣溫日較差每上升1?C，河北衡水和廊坊冬小麥氣象產(chǎn)量增加68～98kg·hm?2，而河北秦皇島和滄州、河南鶴壁及山東德州和東部地區(qū)冬小麥氣象產(chǎn)量減少120～335kg·hm?2。表3顯示開(kāi)花?成熟階段對(duì)冬小麥產(chǎn)量變異影響最大的主要因子是最低氣溫，該階段平均最低氣溫每升高1?C，河南東部和山東西部等地氣象產(chǎn)量增加102～295kg·hm?2，僅北京小麥氣象產(chǎn)量減少104kg·hm?2。

表3 各市不同生長(zhǎng)階段溫度的單位變化（+1℃）對(duì)氣象產(chǎn)量的影響程度（kg·hm?2）

注：每個(gè)市名后括號(hào)中的數(shù)值代表該階段溫度的單位變化引起的氣象產(chǎn)量的變化量即方程6中的偏回歸系數(shù)。+表示增加，?表示減少。下同。

Note: The number in the parenthesis after each city name represents the change in meteorological yield induced by the unit change in temperature, i.e. the partial regression coefficient of Eq.6. “+” means increase, and “–”means decrease. The same as below.

表4 各市不同生長(zhǎng)階段降水的單位變化（+1%）和平均日照時(shí)數(shù)的單位變化（+1 h）對(duì)氣象產(chǎn)量的影響程度（kg·hm?2）

注：每個(gè)市名后括號(hào)中的數(shù)值代表該階段降水和日照的單位變化引起的氣象產(chǎn)量的變化量即方程6中的偏回歸系數(shù)。

Note: The number in the parenthesis after each city name represents the change in meteorological yield induced by the unit change in precipitation and sunshine hours, i.e. the partial regression coefficient of Eq.6.

3 結(jié)論與討論

華北平原1988?2015年冬小麥實(shí)際產(chǎn)量變化在3200～6800kg·hm?2，整體呈現(xiàn)增加趨勢(shì)，主要原因包括化肥施用量的增加和優(yōu)良品種的培育和推廣以及灌溉設(shè)施的建設(shè)等[23]。華北平原冬小麥實(shí)際產(chǎn)量中部地區(qū)高于北部和南部，可能原因是良好的耕作管理措施和土壤性質(zhì)，且中部地區(qū)有完善的灌溉系統(tǒng)，因此是傳統(tǒng)小麥高產(chǎn)區(qū)[24?25]。主要的氣象因素是華北南部地區(qū)生長(zhǎng)季總輻射顯著低于北部和中部地區(qū)[26]，并且南部地區(qū)氣溫高加之氣溫日較差小，小麥呼吸作用強(qiáng)度較大，光合產(chǎn)物更多用于呼吸消耗，導(dǎo)致產(chǎn)量低于中部和北部。

已有較多研究探討了從實(shí)際產(chǎn)量中分離氣象產(chǎn)量的方法[27?28]，但各方法均存在一定的適用范圍和局限性，本研究采用傳統(tǒng)農(nóng)業(yè)氣象學(xué)氣象產(chǎn)量分離方法，該方法已經(jīng)廣泛用于氣象因子變化對(duì)作物產(chǎn)量變化的影響分析中[13,27,29]。分離氣象產(chǎn)量時(shí)，對(duì)趨勢(shì)產(chǎn)量的模擬可以采用多種方法[29?30]，但各方法均存在一定局限性。為便于分析氣象因子變化對(duì)產(chǎn)量變化的影響，本研究主要選擇了常用的6個(gè)氣象因子，結(jié)果表明這些氣象因子對(duì)氣象產(chǎn)量變異影響程度表現(xiàn)為最低溫度＞溫度日較差＞平均氣溫＞日照時(shí)數(shù)＞降水＞最高溫度，反映了華北平原冬小麥產(chǎn)量一定程度上仍受到了低溫的限制。趙凱娜等[32]在河南地區(qū)的研究結(jié)果也表明生育期氣候要素對(duì)冬小麥產(chǎn)量變異的影響程度為最低溫度＞平均氣溫＞最高溫度。肖登攀等[14]的研究表明華北冬小麥整個(gè)生育期除駐馬店站點(diǎn)外，產(chǎn)量與溫度均呈正相關(guān)關(guān)系，而本研究表明不同的生育階段溫度變化對(duì)冬小麥的影響不同，反映了站點(diǎn)尺度與市級(jí)尺度的研究結(jié)果不同，Tao等[21]研究也表明，國(guó)家或者省級(jí)尺度的氣候?產(chǎn)量關(guān)系不能應(yīng)用到市級(jí)尺度上，這反映了氣候?qū)Ξa(chǎn)量影響的復(fù)雜性。本研究與趙凱娜等[31]的研究都表明，相比于其它氣象因子，華北平原冬小麥生長(zhǎng)季內(nèi)降水對(duì)小麥產(chǎn)量影響較小，主要是因?yàn)樵摰貐^(qū)冬小麥普遍以灌溉為主[32]，降水的變化對(duì)產(chǎn)量變異影響較小。同時(shí)，通過(guò)將研究結(jié)果與《中國(guó)氣象災(zāi)害大典》紀(jì)錄進(jìn)行對(duì)比，發(fā)現(xiàn)氣象產(chǎn)量的高低與實(shí)際觀(guān)測(cè)一致，如北京 1987年11月寒潮大風(fēng)降溫，大部分麥田半數(shù)以上葉片被凍傷。河北省1993年滄州、邢臺(tái)連續(xù)22個(gè)月受旱，小麥顯著減產(chǎn)；1994年1月降溫顯著，全省小麥?zhǔn)軆觥Ｉ綎|省1987年12月?1988年4月，小麥?zhǔn)芎递^重的有濰坊、煙臺(tái)、青島、威海，1993年4月青島、淄博部分縣、德州武城出現(xiàn)低溫凍害，導(dǎo)致冬小麥產(chǎn)量降低。河南省1990年淮北地區(qū)小麥凍害嚴(yán)重，1995年1?2月全省大部基本無(wú)降水，嚴(yán)重干旱造成小麥減產(chǎn)。分離的氣象產(chǎn)量在這些年份較低，表明氣象產(chǎn)量的分離結(jié)果具有合理性。

本研究表明在冬小麥播種?越冬階段影響氣象產(chǎn)量變異的主要?dú)庀笠蜃訛榻邓蜏囟热蛰^差，該階段降水在華北大部分地區(qū)對(duì)產(chǎn)量表現(xiàn)為正影響，主要原因是降水增加能夠促進(jìn)小麥生長(zhǎng)和產(chǎn)生更多分蘗[33]。該階段氣溫日較差大對(duì)河南小麥生長(zhǎng)產(chǎn)生了負(fù)效應(yīng)，主要是因?yàn)槿照諘r(shí)數(shù)高，導(dǎo)致作物蒸散增加，土壤墑情變差，影響了作物出苗和苗期生長(zhǎng)。何守法等[34?35]研究也表明，豫中地區(qū)從播種?越冬前期，日照時(shí)數(shù)對(duì)小麥生長(zhǎng)以負(fù)影響為主。越冬?返青階段影響小麥氣象產(chǎn)量變異的主要因子是降水，華北平原大部分區(qū)域該階段降水的增加對(duì)小麥生長(zhǎng)有不利的影響，這個(gè)結(jié)果與張怡等[36]的研究一致，主要因?yàn)樵撾A段小麥停止生長(zhǎng)，蒸散量低，而降水過(guò)程常伴隨著大幅降溫，影響冬小麥的安全越冬，從而影響產(chǎn)量[37]。返青?拔節(jié)階段最低氣溫對(duì)氣象產(chǎn)量的影響均為正效應(yīng)，羅蔣梅等[38?39]研究也表明該階段較高溫度對(duì)小麥產(chǎn)量有益，最低溫度較低會(huì)導(dǎo)致冬小麥遭受晚霜凍害。日照時(shí)數(shù)對(duì)產(chǎn)量影響也存在不同的效應(yīng)，該期小麥需水量增加，日照時(shí)數(shù)低會(huì)減少土壤蒸發(fā)，降低水分脅迫，然而在降水豐沛的地區(qū)，日照時(shí)數(shù)高會(huì)使小麥積累更多干物質(zhì)，促進(jìn)小麥增產(chǎn)。拔節(jié)?開(kāi)花階段對(duì)產(chǎn)量影響最大的是最低氣溫和氣溫日較差，兩因子在不同區(qū)域影響不同。前人研究[34,37]表明，日照時(shí)數(shù)在拔節(jié)?開(kāi)花階段主要是負(fù)影響，因?yàn)槿照諘r(shí)數(shù)高導(dǎo)致作物蒸散高，加劇了春旱的影響。開(kāi)花?成熟階段影響氣象產(chǎn)量變異的關(guān)鍵因子為最低氣溫，最低氣溫除北京外其余均對(duì)小麥產(chǎn)量有正影響，灌漿期較高的最低溫度有利于小麥充分灌漿[40]。

華北北部冬小麥產(chǎn)量的主要?dú)夂蛳拗埔蜃訛榘喂?jié)?開(kāi)花階段的低溫。該區(qū)域要做好防凍保溫措施，可通過(guò)調(diào)整播種期[41]和采用覆膜[42]等方法避免冬小麥遭受凍害。華北中部地區(qū)冬小麥產(chǎn)量的主要?dú)夂蛳拗埔蛩貫椴シN?越冬階段的干旱，越冬?返青階段較多的降水，返青?開(kāi)花期的光照不足以及灌漿期的低溫。華北東部冬小麥產(chǎn)量的主要?dú)夂蛳拗埔蛩貫椴シN?越冬階段的干旱，越冬?返青階段較多的降水，返青?拔節(jié)階段的低溫。中部和東部地區(qū)應(yīng)在冬小麥播種?越冬期增加灌溉量，為防止作物蒸散過(guò)高，可采取一些保水措施，從而增加土壤墑情，保證小麥出苗率。華北南部冬小麥產(chǎn)量的主要?dú)夂蛳拗埔蛩貫椴シN?返青階段較多的降水，返青?拔節(jié)階段的低溫以及灌漿期的低溫。南部地區(qū)可以通過(guò)采用抗低溫的品種[43]，同時(shí)充分利用小麥生長(zhǎng)階段光熱資源充足的優(yōu)勢(shì)進(jìn)行播量和群體的調(diào)控管理來(lái)實(shí)現(xiàn)高產(chǎn)。

大量研究表明，氣候變化對(duì)作物產(chǎn)量有重要影響，但本研究的主要目標(biāo)是揭示影響產(chǎn)量變異的主要?dú)庀笠蜃?，而產(chǎn)量變異的影響因子主要是氣候變率而不是氣候變化。因此，在分離氣象產(chǎn)量時(shí)未考慮氣候變化對(duì)趨勢(shì)產(chǎn)量的可能影響[44?45]。1988?2015年，39%的市氣溫呈顯著上升趨勢(shì)，而降水量?jī)H在開(kāi)封市顯著減少，17%的市日照時(shí)數(shù)顯著減少，其變化特點(diǎn)與本研究氣象產(chǎn)量變化影響因子分析的結(jié)果保持一致。分析氣象因子與氣象產(chǎn)量的關(guān)系時(shí)采用多元線(xiàn)性逐步回歸方法，該方法能很好地排除自相關(guān)性較高的變量，且對(duì)最后影響氣象產(chǎn)量變異結(jié)果中的氣象因子進(jìn)行了相關(guān)性分析，將其中相關(guān)的因子剔除后再與氣象產(chǎn)量作回歸分析，結(jié)果與之前的結(jié)論保持一致。本研究在各個(gè)市直接選用單個(gè)氣象站點(diǎn)來(lái)代表整個(gè)市級(jí)氣象數(shù)據(jù)，而沒(méi)有使用空間插值再分區(qū)統(tǒng)計(jì)等方法，主要原因是華北平原地形平坦，氣象要素分布均勻，這46個(gè)氣象站點(diǎn)一定程度上代表了該市冬小麥種植區(qū)的主要?dú)夂驐l件，因此，利用本市的氣象數(shù)據(jù)和產(chǎn)量作回歸分析具有科學(xué)性，前人研究也有采用此類(lèi)方法[34,37,46]。其次，當(dāng)空間插值的數(shù)據(jù)點(diǎn)不均勻分布時(shí)，插值后進(jìn)行平均難以很好地反映研究要素的空間分布特征[47]，且采用插值后進(jìn)行平均得到的結(jié)果要低于當(dāng)前研究結(jié)果的精度。進(jìn)一步研究應(yīng)分析影響華北平原實(shí)際產(chǎn)量變異的自然和社會(huì)經(jīng)濟(jì)等因子，為提高華北平原冬小麥產(chǎn)量和資源利用效率提供科學(xué)參考。

[1] 中華人民共和國(guó)國(guó)家統(tǒng)計(jì)局.中國(guó)統(tǒng)計(jì)年鑒2013[M]. 北京:中國(guó)統(tǒng)計(jì)出版社,2013. National Bureau of Statistics of the People’s Republic of China. China Statistical Yearbook 2013[M].Beijing:China Statistics Press,2013.(in Chinese)

[2] 楊曉光,劉志娟.作物產(chǎn)量差研究進(jìn)展[J].中國(guó)農(nóng)業(yè)科學(xué),2014,47(14):2731-2741.Yang X G,Liu Z J.Advances in research on crop yield gaps[J].Scientia Agricultura Sinica,2014,47(14):2731-2741.(in Chinese)

[3] 劉巽浩,胡躍高,陳阜.糧食:中國(guó)農(nóng)業(yè)的永恒難題與主題[J].農(nóng)業(yè)現(xiàn)代化研究,2010,31(4):385-391.Liu X H,Hu Y G,Chen F.The problem of cereals: eternal theme and difficult problem of Chinese agriculture[J].Research of Agricultural Modernization,2010,31(4):385-391.(in Chinese)

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

[5] 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.

[6] Schapendonk A,Xu H Y,Putten P,et al.Heat-shock effects on photosynthesis and sink-source dynamics in wheat(Triticum aestivum L.)[J].NJAS-Wageningen Journal of Life Sciences,2008,55(1):37-54.

[7] 田展,梁卓然,史軍,等.近50年氣候變化對(duì)中國(guó)小麥生產(chǎn)潛力的影響分析[J].中國(guó)農(nóng)學(xué)通報(bào),2013,29(9):61-69.Tian Z,Liang Z R,Shi J,et al.Analysis of impact on China wheat potential productivity of climate change during 1961-2010[J].Chinese Agricultural Science Bulleti,2013,29(9):61-69.(in Chinese)

[8] Fang S B,Cammarano D,Zhou G S,et al. Effects of increased day and night temperature with supplemental infrared heating on winter wheat growth in North China[J].European Journal of Agronomy,2015,64:67-77.

[9] 鄭娜,劉秀位,王錫平.利用產(chǎn)量差距方法進(jìn)行華北冬小麥產(chǎn)量受氣候影響規(guī)律的分析:以河北省欒城縣為例[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2014,22(2):234-240. Zheng N,Liu X W,Wang X P.Determining the effects of climate on winter wheat yield in Northern China via yield gap analysis:a case study of Luancheng county,Hebei province[J].Chinese Journal of Eco-Agriculture,2014,22(2):234-240.(in Chinese)

[10]閆素輝,李文陽(yáng),楊安中,等.弱光對(duì)小麥花后旗葉光合及籽粒灌漿的影響[J].麥類(lèi)作物學(xué)報(bào),2011,31(1):77-81.Yan S H,Li W Y,Yang A Z,et al.Effects of weak light at grain filling stage on photosynthetic characteristics and grain filling of winter wheat[J].Journal of Triticeae Crops,2011,31(1):77-81.(in Chinese)

[11]張黎萍,荊奇,戴廷波,等.溫度和光照強(qiáng)度對(duì)不同品質(zhì)類(lèi)型小麥旗葉光合特性和衰老的影響[J].應(yīng)用生態(tài)學(xué)報(bào),2008,19(2):311-316.Zhang L P,Jing Q,Dai T B,et al.Effects of temperature and illumination on flag leaf photosynthetic characteristics and senescence of wheat cultivars with different grain quality[J].The Journal of Applied Ecology,2008,19(2):311-316.(in Chinese)

[12]譚凱炎,楊曉光,任三學(xué),等.高溫脅迫對(duì)華北地區(qū)冬小麥灌漿及產(chǎn)量的影響[J].生態(tài)學(xué)報(bào),2015,35(19):6355-6361. Tan K Y,Yang X G,Ren S X,et al.Impact of high temperature stress at the grain-filling stage on winter wheat yield[J].Acta Ecologica Sinica,2015,35(19):6355-6361.(in Chinese)

[13]許世衛(wèi),喻聞,劉淑云,等.冬小麥單產(chǎn)與氣象因素關(guān)系研究:以山東省微山縣為例[J].農(nóng)學(xué)學(xué)報(bào),2013,3(10):45-51. Xu S W,Yu W,Liu S Y,et al.The impact of meteorological factors on the winter wheat yield:the case study from Weishan county, Shandong province[J].Journal of Agricul- ture,2013,3(10):45-51.(in Chinese)

[14]肖登攀,陶福祿,沈彥俊,等.華北平原冬小麥對(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)

[14]Asseng S,Jamieson P D,Kimball B,et al.Simulated wheat growth affected by rising temperature, increased water deficit and elevated atmospheric CO2[J].Field Crops Resear- ch,2004,85(2-3):85-102.

[15]吳東麗,王春乙,薛紅喜,等.華北地區(qū)冬小麥干旱風(fēng)險(xiǎn)區(qū)劃[J].生態(tài)學(xué)報(bào),2011,31(3):760-769. Wu D L,Wang C Y,Xue H X,et al.The drought risk zoning of winter wheat in North China[J].Acta EcologicaSinica,2011,31(3):760-769.(in Chinese)

[16]王素艷,霍治國(guó),李世奎,等.北方冬小麥干旱災(zāi)損風(fēng)險(xiǎn)區(qū)劃[J].作物學(xué)報(bào),2005,31(3):267-274. Wang S Y,Huo Z G,Li S K,et al.Risk regionalization of winter wheat loss caused by drought in North of China[J].Acta Agronomica Sinica,2005,31(3):267-274.(in Chinese)

[17]王春乙,張玉靜,張繼權(quán).華北地區(qū)冬小麥主要?dú)庀鬄?zāi)害風(fēng)險(xiǎn)評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(s1):203-213. Wang C Y,Zhang Y J,Zhang J Q.Risk assessment of main meteorological disasters of winter wheat in North China[J].Transactions of the CSAE,2016,32(s1):203-213.(in Chinese)

[18]Lobell D B,Schlenker W,Costa R.Climate trends and global crop production since 1980[J].Science,2011,333(6042):616-620.

[19]Lobell D B,Field C B.Global scale climate-crop yield relationships and the impacts of recent warming[J].Environ- mental Research Letters,2007,2(1):1-7.

[20]You L Z,Rosegrant M W,Wood S,et al.Impact of growing season temperature on wheat productivity in China[J].Agricultural and Forest Meteorology,2009,149(6):1009-1014.

[21]Tao F L,Yokozawa M,Liu J Y,et al.Climate-crop yield relationships at provincial scales in China and the impacts of recent climate trends[J].Climate Research,2008,38(1):83-94.

[22]楊繼武.農(nóng)業(yè)氣象預(yù)報(bào)和情報(bào)[M].北京:氣象出版社,1994. Yang J W.Agricultural meteorological forecast and informa- tion[M].Beijing:China Meteorology Press,1994.(in Chinese)

[23]傅友蘭,屈寶香,李剛,等.中國(guó)糧食產(chǎn)量變化與年際波動(dòng)特點(diǎn)分析[J].中國(guó)食物與營(yíng)養(yǎng),2012,18(1):16-19. Fu Y L,Qu B X,Li G,et al.Changes of food production and the characteristics of inter-annual fluctuations in China[J].Food and Nutrition in China,2012,18(1):16-19.(in Chinese)

[24]Lu C,Fan L.Winter wheat yield potentials and yield gaps in the North China Plain[J].Field Crops Research,2013,143(1):98-105.

[25]莫興國(guó),林忠輝,劉蘇峽.黃淮海地區(qū)冬小麥生產(chǎn)力時(shí)空變化及其驅(qū)動(dòng)機(jī)制分析[J].自然資源學(xué)報(bào),2006,21(3):449-457.Mo X G,Lin Z H,Liu S X.Spatial-temporal evolution and driving forces of winter wheat productivity in the Huang-Huai-Hai Region[J].Journal of Natural Resources,2006,21(3):449-457.(in Chinese)

[26]李克南.華北地區(qū)冬小麥-夏玉米作物生產(chǎn)體系產(chǎn)量差特征解析[D].北京:中國(guó)農(nóng)業(yè)大學(xué),2014. Li K N.Yield gap analysis focused on winter wheat and summer maize rotation in the North China Plain[D].Beijing:China Agricultural University,2014.(in Chinese)

[27]房世波.分離趨勢(shì)產(chǎn)量和氣候產(chǎn)量的方法探討[J].自然災(zāi)害學(xué)報(bào),2011,20(6):13-18. Fang S B.Exploration of method for discrimination between trend crop yield and climatic fluctuant yield[J].Journal of Natural Disasters,2011,20(6):13-18.(in Chinese)

[28]王媛,方修琦,徐錟.氣候變化背景下“氣候產(chǎn)量”計(jì)算方法的探討[J].自然資源學(xué)報(bào),2004,19(4):531-536. Wang Y,Fang X Q,Xu T.A method for calculating the climatic yield of grain under climate change[J].Journal of Natural Resources,2004,19(4):531-536.(in Chinese)

[29]王桂芝,陸金帥,陳克垚,等.基于HP濾波的氣候產(chǎn)量分離方法探討[J].中國(guó)農(nóng)業(yè)氣象,2014,35(2):195-199. Wang G Z,Lu J S,Chen K Y,et al.Exploration of method in separating climatic output based on HP filter[J].Chinese Journal of Agrometeorology,2014,35(2):195-199.(in Chinese)

[30]趙東妮,王艷華,任傳友,等.3種水稻趨勢(shì)產(chǎn)量擬合方法的比較分析[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2017,25(3):345-355.Zhao D N,Wang Y H,Ren C Y,et al.Comparative analysis of three fitting methods of rice trend yield[J].Chinese Journal of Eco-Agriculture,2017,25(3):345-355.(in Chinese)

[31]趙凱娜,寧曉菊,秦耀辰,等.縣域冬小麥生育期氣候要素與產(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)

[32]劉濤,周廣勝,譚凱炎,等.華北地區(qū)冬小麥灌溉制度及其環(huán)境效應(yīng)研究進(jìn)展[J].生態(tài)學(xué)報(bào),2016,36(19):5979-5986.Liu T,Zhou G S,Tan K Y,et al.Review on research of irrigation regime and its environmental effect in winter wheat field of North China Plain[J].Acta Ecologica Sinica,2016,36(19):5979-5986.(in Chinese)

[33]柳芳,黎貞發(fā).降水量和積溫變化對(duì)天津冬小麥產(chǎn)量的影響[J].中國(guó)農(nóng)業(yè)氣象,2010,31(3):431-435.Liu F,Li Z F.Impacts of precipitation and accumulated temperature on winter wheat yield in tianjin[J].Chinese Jour- nal of Agrometeorology,2010,31(3):431-435.(in Chinese)

[34]何守法,董中東,詹克慧,等.豫中補(bǔ)灌區(qū)主要?dú)夂蛞蜃訉?duì)小麥產(chǎn)量的影響[J].河南農(nóng)業(yè)科學(xué),2009,19(7):47-49.He S F,Dong Z D,Zhan K H,et al.Influence of main climatic factors on yield of wheat in irrigation area of Henan province[J].Henan Agricultural Sciences,2009,19(7):47-49.(in Chinese)

[35]趙孟偉.氣象因素對(duì)開(kāi)封地區(qū)冬小麥產(chǎn)量影響的通徑分析[J].華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,36(4):26-31.Zhao M W.Path analysis on the influence of meteorological factors on winter wheat yield in the Kaifeng area[J].Journal of North China University of Water Resources & Electric Power,2015,36(4):26-31.(in Chinese)

[36]張怡,史本林.氣候變化對(duì)豫東冬小麥產(chǎn)區(qū)的綜合影響[J].江蘇農(nóng)業(yè)科學(xué),2015,43(2):336-339.Zhang Y,Shi B L.Integrated effects of climate change on winter wheat production in east Henan[J].Jiangsu Agricul- tural Science,2015,43(2):336-339.(in Chinese)

[37]劉淑云,谷衛(wèi)剛,尚明華,等.魯西北生育期降水與小麥產(chǎn)量相關(guān)關(guān)系回歸模型研究[J].作物雜志,2013,24 (6):90-94. Liu S Y,Gu W G,Shang M H,et al.Optimum regression model between rainfall in wheat growth period and wheat yield[J].Crops,2013,24(6):90-94.(in Chinese)

[38]羅蔣梅,王建林,申雙和,等.影響冬小麥產(chǎn)量的氣象要素定量評(píng)價(jià)模型[J].南京氣象學(xué)院學(xué)報(bào),2009,32(1):94-99.Luo J M,Wang J L,Shen S H,et al.Quantitative evaluating model for meteorological factor effects on winter wheat yield[J].Journal of Nanjing Institute of Meteorology,2009,32(1):94-99.(in Chinese)

[39]錢(qián)錦霞,郭建平.鄭州地區(qū)冬小麥產(chǎn)量構(gòu)成要素的回歸模型[J].應(yīng)用氣象學(xué)報(bào),2012,23(4):500-504. Qian J X,Guo J P.Regression models of winter wheat yield components in Zhengzhou area[J].Journal of Applied Meteorological Science,2012,23(4):500-504.(in Chinese)

[40]Lobell D B,Sibley A,Ortizmonasterio J I.Extreme heat effects on wheat senescence in India[J].Nature Climate Change,2012,2(3):186-189.

[41]楊桂霞,趙廣才,許軻,等.播期和密度對(duì)冬小麥籽粒產(chǎn)量和營(yíng)養(yǎng)品質(zhì)及生理指標(biāo)的影響[J].麥類(lèi)作物學(xué)報(bào),2010,30(4):687-692. Yang G X,Zhao G C,Xu K,et al.Effect of sowing date and density on grain yield and nutrition quality and physiolo- gical index of winter wheat[J].Journal of Triticeae Crops,2010,30(4):687-692. (in Chinese)

[42]王同花,李援農(nóng).起壟覆膜溝播對(duì)冬小麥土壤水分利用效率及土壤溫度的影響[J].中國(guó)農(nóng)村水利水電,2011,(1):46-49. Wang T H,Li Y N.Effect of furrow planting winter wheat with plastic film mulching on water-use efficiency and soil temperature[J].China Rural Water and Hydropower,2011,(1):46-49.(in Chinese)

[43] 李曉林,白志元,楊子博,等.黃淮麥區(qū)部分主推冬小麥品種越冬及拔節(jié)期的抗寒生理研究[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,41(1):40-48. Li X L,Bai Z Y,Yang Z B,et al.Cold resistant physiology of some main wheat varieties at wintering and jointing stages in Huanghuai area[J].Journal of Northwest A&F University,2013,41(1):40-48.(in Chinese)

[44] 高永剛,顧紅,姬菊枝,等.近43年來(lái)黑龍江氣候變化對(duì)農(nóng)作物產(chǎn)量影響的模擬研究[J].應(yīng)用氣象學(xué)報(bào),2007,18(4):532-538.Gao Y G,Gu H,Ji J Z,et al.Simulation study of climate change impact on crop yield in Heilongjiang Province from 1961to 2003[J].Journal of Natural Resources,2007,18(4):532-538.(in Chinese)

[45] 張明偉,鄧輝,李貴才,等.模型模擬華北地區(qū)氣候變化對(duì)冬小麥產(chǎn)量的影響[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃,2011,32(4):45-49. Zhang M W,Deng H,Li G C,et al.Simulation study on impacts of climate change on the production of winter wheat in North China[J].Chinese Journal of Agricultural Resources and Regional Planning, 2011,32(4):45-49.(in Chinese)

[46] 熊淑萍,盧中民,吳克遠(yuǎn),等.河南省小麥和玉米產(chǎn)量及其有關(guān)氣候因子分析[J].河南農(nóng)業(yè)大學(xué)學(xué)報(bào),2014,48(2):123-129.Xiong S P,Lu Z M,Wu K Y,et al.Analysis on yield of wheat and maize and related climate factorsin Henan Province[J].Journal of Henan Agricultural University,2014,48(2):123-129.(in Chinese)

[47] 朱求安,張萬(wàn)昌,余鈞輝.基于GIS的空間插值方法研究[J].江西師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2004,28(2):183-188. Zhu Q A,Zhang W C,Yu J H.The Spatial Interpolations in GIS[J].Journal of Jiangxi Normal University,2004,28(2):183-188.(in Chinese)

Meteorological Influencing Factors on Variation in Winter Wheat Yield in the North China Plain

WU Bing-jie1, 2, WANG Jing1, TANG Jian-zhao1, WANG Na1, 3, XU Lin1, BAI Hui-qing1, ZHENG Jun-qing1, WANG Na1, LI Yang1

(1. College of Resources and Environmental Science, China Agricultural University, Beijing 100193, China; 2. Chengdu Jintian Lianhua Technology Co., LTD, Chengdu 610072; 3. Crop Research Institute, Shandong Academy of Agricultural Sciences, Jinan 250100)

The meteorological influencing factors on variation in winter wheat yield in the North China Plain (NCP) were investigated by developing yield?meteorological factor relationship based on statistical yield of 46 cities and daily meteorological data from 46 meteorological stations from 1988 to 2015. The results showed that statistical yield of winter wheat ranged from 3200 to 6800kg·ha?1in 1988?2015 with higher yield in the central NCP and higher yield variation in the southern NCP. The variation in growing season main meteorological factors including sunshine hours, temperature and precipitation could account for 17%?78% of variation in the climatic yield of winter wheat with statistical significance at 54% of the cities (P<0.05). The regions with high influencing degree were mainly concentrated in the southern Hebei, southern Shandong and northeastern Henan. Precipitation during the growing period of sowing to turning?green had significant impacts on meteorological yield of winter wheat. Meteorological yield of winter wheat increased by 13?74kg·ha?1in Tianjin, Zhumadian and northwestern Shandong if amount of precipitation during the period increased by 1%, while the meteorological yield in southern Henan and southern Shandong decreased by 16?80kg·ha?1. Minimum temperature during the growing period of turning-green to maturity had significant impacts on variation in meteorological yield of winter wheat. Average minimum temperature during the period increased by 1?C would increase meteorological yield by 50?295kg·ha?1in Tianjin, Shijiazhuang, the eastern and western Shandong, and the western Henan, but would decrease meteorological yield by 76?124kg·ha?1in Beijing, Tangshan and Zaozhuang. In general, the temperature had significant impact on the variati on of winter wheat yield at more sites than sunshine hours and precipitation. However, there was a large spatial difference in meteorological influencing factors on variation in winter wheat yield due to complex interaction of meteorological factors, cultivar and agronomic management practices.

Yield variation; Meteorological yield; Sunshine hours; Temperature; Precipitation

10.3969/j.issn.1000-6362.2018.10.001

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

2018?03?06

。E-mail：wangj@cau.edu.cn

國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2016YFD0300105；2017YFD0300404）

吳冰潔（1994?），女，碩士，主要從事氣候變化對(duì)作物的影響及適應(yīng)研究。E-mail：wubjok@163.com