楊素苗 董宇航 馬筱建 石海強(qiáng) 杜紀(jì)壯 姚樹然 段鵬偉

摘? ? 要：為探明氣象因子對果實品質(zhì)的影響，確定適宜‘富士蘋果栽培的氣候條件范圍，應(yīng)用主成分分析、多元線性回歸及線性規(guī)劃對河北省44個蘋果園及當(dāng)?shù)貧庀笠蜃舆M(jìn)行分析。結(jié)果表明，氣候因子對果實的著色指數(shù)、裂果指數(shù)、質(zhì)構(gòu)特性影響顯著，對單果質(zhì)量、可溶性固形物含量、可滴定酸含量沒有顯著影響;適宜‘富士蘋果栽培的氣候條件為生長期日均氣溫17.05～18.42 ℃，平均日較差9.94～14.61 ℃，日均降水量1.29～3.09 mm，日均濕度55.01%～71.55%，日均日照時數(shù)5.72～5.92 h，及年日均氣溫8.38～15.60 ℃，平均日較差9.59～14.26 ℃，日均降水量0.79～1.93 mm，日均濕度50.15%～65.67%，日均日照時數(shù)5.19～7.48 h。綜上，河北省生長期內(nèi)的日較差、降水量和濕度，以及周年的溫度、日較差、降雨量、濕度和光照時數(shù)均適宜‘富士蘋果栽培，部分地區(qū)生長期內(nèi)過高的日照時數(shù)和溫度是影響果實品質(zhì)的關(guān)鍵氣象因子。

關(guān)鍵詞：蘋果;果實品質(zhì);主成分回歸;線性規(guī)劃;氣象因子

中圖分類號：S162.5+5? ? ? ? ?文獻(xiàn)標(biāo)識碼：A? ? ? ? ? DOI 編碼：10.3969/j.issn.1006-6500.2021.03.014

The Effects of Climate Factors on Fruit Quality of 'Fuji' Apple in Hebei Province

YANG Sumiao1， DONG Yuhang2， MA Xiaojian1， SHI Haiqiang1， DU Jizhuang1， YAO Shuran2， DUAN Pengwei1

（1.Shijiazhuang Institute of Fruit Trees， Hebei Academy of Agriculture and Forestry Sciences， Shijiazhuang， Hebei 050000， China; 2.Meteorological Institute of Hebei Province， Shijiazhuang，? Hebei? 050000， China）

Abstract： In order to ascertain the impact of meteorological factors on fruit quality and determine the suitable range of climatic conditions for 'Fuji' apple cultivation， 44 apple orchards and local meteorological factors in Hebei Province were analyzed through the principal component analysis， multiple linear regression and linear programming. The results showed that climate factors had significant effects on fruit color index， cracking index， and texture characteristics， but had no significant effects on fruit weight， soluble solid content， and titratable acid content. The suitable climatic conditions for 'Fuji' apple cultivation were as follows： the average daily temperature was 17.05-18.42℃， the average diurnal temperature range was 9.94-14.61 ℃， the average daily precipitation was 1.29-3.09 mm， the average daily humidity was 55.01%-71.55%， the average daily sunshine hours was 5.72-5.92 h during the growth period;the average daily temperature was 8.38 -15.60 ℃， the average daily temperature difference was 9.59-14.26 ℃， the average daily precipitation was 0.79-1.93 mm， the average daily humidity was 50.15%-65.67%， and the average daily sunshine hours was 5.19-7.48 h in the whole year. In summary， the diurnal temperature range， precipitation and humidity during the growth period， as well as the annual temperature， diurnal temperature range， rainfall， humidity and light hours are suitable for the cultivation of 'Fuji' apples in Hebei Province， and the key meteorological factors affecting fruit quality are the too long sunshine hours and the too high temperature during the growth period in some areas.

Key words： apple; fruit quality; principal component regression; linear programming; climate factor

中國是蘋果生產(chǎn)大國，種植面積和產(chǎn)量均居世界第一?！皇渴俏覈闹髟云贩N，約占蘋果種植品種的70%[1]，但通過對主產(chǎn)區(qū)‘富士蘋果的抽樣檢測發(fā)現(xiàn)，我國蘋果的理化品質(zhì)呈下降趨勢[2-5]。果實品質(zhì)受氣候條件影響較大，Sugiura等[6]發(fā)現(xiàn)氣候變暖導(dǎo)致‘富士和‘津輕蘋果花期提前、果實成熟度變高，果實硬度和酸含量降低，可溶性固形物含量提高。申順吏等[7]認(rèn)為溫度類因子是蘋果產(chǎn)量的主要影響因素;Syukur等[8]研究發(fā)現(xiàn)，隨著溫度的升高，結(jié)果數(shù)量、果實產(chǎn)量、果品級別、果實硬度等均降低。童盼盼等[9]研究認(rèn)為，果實品質(zhì)和糖心形成受光照強(qiáng)度和晝夜溫差的影響，光照強(qiáng)度與晝夜溫差大有利于果實品質(zhì)提高和糖心形成。張磊等[10]研究結(jié)果表明，氣溫、降水、日照時數(shù)和相對濕度是影響高酸蘋果品質(zhì)最主要的氣象因子。魏欽平等[11-12]研究認(rèn)為單果質(zhì)量和果形指數(shù)的相對光照強(qiáng)度值要求較低，果實著色、可溶性固形物和可溶性糖對相對光照強(qiáng)度的需求均在80%以上。李興軍等[13]研究發(fā)現(xiàn)蘋果果實生育期的相對濕度不宜超過70%;段曉鳳等[14]認(rèn)為“富士”蘋果的相對濕度要求為65.1%～73.3%。綜上，同一氣候因子可能影響果實的多個品質(zhì)，同一品質(zhì)指標(biāo)也可能同時受多個氣候因子影響。為了探明氣象因子對果實品質(zhì)的影響，明確適宜‘富士蘋果栽培的氣象因子范圍，以我國環(huán)渤海蘋果優(yōu)勢產(chǎn)區(qū)的河北省44個果園的‘富士蘋果為試材，結(jié)合所屬縣的氣象數(shù)據(jù)，分析了氣象因子對果實品質(zhì)的影響。

1 材料和方法

1.1 試驗材料

‘富士蘋果果實選自河北省33個蘋果種植縣共44個‘富士蘋果園，南起邯鄲市磁縣，北至承德市平泉市，北緯36°41′至41°20′，東經(jīng)113°78′至119°27′，果園海拔高度范圍為-5～700 m。2019年河北省年平均氣溫12.9 ℃，年平均降水量437.6 mm，年平均日照時數(shù)2 373.3 h。

自2019年10—11月，分別于當(dāng)?shù)毓麑嵶钸m采收期，在每個果園隨機(jī)選擇長勢中庸的成齡長枝富士蘋果樹3株，每株樹從4個方向隨機(jī)摘取10個果實，即每個果園共取30個果實，取樣后4 h內(nèi)保存于果品保鮮庫中。

1.2 試驗方法

1.2.1 果實品質(zhì)指標(biāo) 對果實單果質(zhì)量、著色指數(shù)、裂果指數(shù)、可溶性固形物含量、可滴定酸含量、脆裂性、膠粘性、咀嚼性共8項指標(biāo)進(jìn)行測定。

單果質(zhì)量，根據(jù)NY/T 2316-2013蘋果品質(zhì)指標(biāo)評價規(guī)范稱量。著色指數(shù)，目測著色部分在整個果實的面積占比。裂果指數(shù)，目測裂果部分在整個果實的面積占比?？扇苄怨绦挝?，采用PAL-1型數(shù)字糖度計測定?？傻味ㄋ岷?，采用GMK-835F型水果酸度計測定。脆裂性是指樣品折斷時的力，膠粘性是指半固體食品吞咽前破碎它需要的能量，咀嚼性是指咀嚼固體樣品時需要的能量，感官指果實在牙齒的咀嚼過程中對外力的持續(xù)抵抗作用;脆裂性、膠粘性和咀嚼性用美國FTC質(zhì)構(gòu)儀TPA（texture profile analysis）法測定，方法：延果梗將蘋果切為兩瓣，使用14 mm的打孔器在蘋果的4個角取樣，用切分寬度4.5 mm的雙刀切取居中部位作為試樣。質(zhì)構(gòu)儀感應(yīng)元量程250 N，圓柱形探頭，回升高度30 mm，形變百分量60，檢測速度60 mm·min-1，最小力0.1 N。

1.2.2 氣象因子主成分分析 主成分回歸是一種處理多重共線性的新參數(shù)方法，其基本思想是首先將原始變量轉(zhuǎn)化為若干主成分，再建立主成分的回歸方程，然后再基于主成分與原始變量之間的對應(yīng)關(guān)系，得到原自變量回歸模型的估計方程。選擇生長期內(nèi)（2019年4—10月）日均氣溫、平均日較差、日均降水、日均濕度、日均日照時數(shù)和一年內(nèi)（2019年1—12月）日均氣溫、平均日較差、日均降水、日均日照時數(shù)，共10項指標(biāo)（依次記做X1，X2……X10）進(jìn)行主成分分析。

氣象數(shù)據(jù)來自果園自建氣象站或距其最近的鄉(xiāng)鎮(zhèn)（縣城）氣象站。

1.2.3 多元線性回歸分析 對測定的8項果實品質(zhì)指標(biāo)（依次記做Y1，Y2……Y8）和氣象因子主成分進(jìn)行回歸分析，并將該回歸方程中的主成分值還原為X表達(dá)。

1.2.4 線性規(guī)劃 根據(jù)Y與X之間的關(guān)系，求得氣象因子的最優(yōu)解，得到有利于‘富士蘋果品質(zhì)形成的適宜氣候條件。

1.3 數(shù)據(jù)分析

采用Microsoft Excel 2019進(jìn)行數(shù)據(jù)統(tǒng)計與作圖，采用SPSS 22.0進(jìn)行主成分回歸分析，采用Lingo 12.0進(jìn)行線性規(guī)劃。

2 結(jié)果與分析

2.1 果實品質(zhì)指標(biāo)分析

從表1可知，蘋果單果質(zhì)量介于321.93 g和143.60 g之間，均值為247.53 g，極差為178.33 g;蘋果著色指數(shù)均值為82.08%，超過均值的果園共有26個，占59.09%，達(dá)到100%的果園有4個;裂果方面，平均裂果指數(shù)為12.47%，無裂果的果園有13個，占29.55%，裂果指數(shù)不超過10%的果園有24個，占54.55%;果實可溶性固形物含量的均值為13.59%，最高值為17.44%，最低值為10.51%，極差為6.93%;果實可滴定酸含量的均值為0.21%;脆裂性平均值為61.47 N，最大值為81.64 N，最小值為47.20 N，極差為34.44 N;膠粘性平均值為8.44 N，最大值和最小值分別為13.46 N和6.00 N，極差為7.46 N;咀嚼性介于22.48 J和8.94 J之間，平均值為14.81 J，極差為13.54 J。7個指標(biāo)的變異系數(shù)表現(xiàn)為裂果指數(shù)>咀嚼性>膠粘性>著色指數(shù)>可滴定酸含量>單果質(zhì)量>脆裂性>可溶性固形物含量。

2.2 氣象因子數(shù)據(jù)及相關(guān)性分析

由表2可知，生長季的平均氣溫最大值和最小值分別為17.05 ℃和23.28 ℃，日較差均值為11.61 ℃，年平均氣溫的最大值和最小值分別為8.38 ℃和15.6 ℃，日較差均值為13.55 ℃。生長季總降水量均值為402.32 mm，年總降水量均值為423.4 mm，生長季降雨量占全年降雨量的95.02%。年降雨量最小為288.35 mm，年降雨量最大為704.45 mm，前者與后者的比值為41%，可見年降雨量的區(qū)域差異較大?？諝鉂穸鹊木?0.15%～71.55%，年日照時數(shù)均值1 894.35～2 730.20 h。不同氣象因子的變異系數(shù)表現(xiàn)為全年日均降水量>生長期日均降水>全年日均氣溫>全年平均日較差>生長期平均日較差>全年日均日照時數(shù)>全年日均濕度>生長期日均氣溫>生長期日均濕度>生長期日均日照時數(shù)。

由表3可知，全年和生長季的氣溫、日較差、日照時數(shù)3者之間相關(guān)性均為極顯著。年降水量與年日均濕度、年日均光照時數(shù)相關(guān)性極顯著，生長季日均氣溫與日照時數(shù)同樣呈極顯著相關(guān)。通常認(rèn)為，二因子相關(guān)系數(shù)大于0.7時，可能存在共線性，因此建立回歸方程需要對氣象因子數(shù)據(jù)進(jìn)行偏回歸處理。

2.3 氣象因子的主成分提取

通過對生長期內(nèi)日均氣溫、平均日較差、日均降水、日均濕度、日均日照時數(shù)和一年內(nèi)日均氣溫、平均日較差、日均降水、日均日照時數(shù)等氣象因子進(jìn)行主成分轉(zhuǎn)化后，得到10個主成分因子及其對應(yīng)的氣象因子特征值（表4）?？芍?，前5個主成分的累計貢獻(xiàn)率達(dá)到99.09%，滿足統(tǒng)計學(xué)分析要求。

2.4 果實品質(zhì)與氣象因子主成分的回歸分析

建立果實各品質(zhì)指標(biāo)與5個主成分的回歸方程，再根據(jù)表4中主成分與氣象因子的關(guān)系式，將主成分轉(zhuǎn)換為氣象因子，得到果實單果質(zhì)量、著色指數(shù)、裂果指數(shù)、可溶性固形物含量、可滴定酸含量、脆裂性、膠粘性和咀嚼性與氣象因子的回歸關(guān)系（表5）。根據(jù)顯著相關(guān)性判定，果實的著色指數(shù)、裂果指數(shù)、脆裂性、膠粘性、咀嚼性回歸方程均具有統(tǒng)計學(xué)意義，即上述5項指標(biāo)與氣象因子相關(guān)性顯著。通過對擬合值和實際值進(jìn)行對比，見圖1—圖5，可以看出大部分果園的實際值與擬合值基本一致。

2.5 確定氣象因子適宜范圍

根據(jù)氣象因子與品質(zhì)的線性規(guī)劃方程組，求解氣象因子的范圍，既要保證氣象因子在測定范圍之內(nèi)，又要確保果實品質(zhì)指標(biāo)達(dá)到平均值以上水平，即果實品質(zhì)因素的選擇為果實著色指數(shù)≥82.08%、裂果指數(shù)≤12.47%、脆裂性≥61.47 N、膠粘性≥8.44 N、咀嚼性≥14.81 J。通過Lingo軟件依次進(jìn)行計算，結(jié)果見表6，可知氣象因子的適宜范圍：生長期的日均氣溫17.05～18.42 ℃，平均日較差9.94～14.61 ℃，日均降水量1.29～3.09 mm，日均濕度55.01%～71.55%，日均日照時數(shù)5.72～5.92 h;年日均氣溫8.38～15.60 ℃，平均日較差9.59～14.26 ℃，日均降水量0.79～1.93 mm，日均濕度50.15%～65.67%，日均日照時數(shù)5.19～7.48 h。

3 結(jié)論與討論

主成分回歸法降維優(yōu)勢明顯，具有多重共線性的自變量經(jīng)過主成分轉(zhuǎn)換后，消除了共線性的危害[15-16]。本文中氣候因子選擇了全年和生長期兩類時期的各5項指標(biāo)，部分氣象因子之間相關(guān)性顯著，通過主成分提取及回歸分析，降低了模型擬合的誤差，在一定程度上會使得回歸方程更加可靠。使用Lingo軟件進(jìn)行線性規(guī)劃，在求解氣象因子范圍時，將所有相關(guān)品質(zhì)的條件及氣象因子范圍都列為限制條件，然后分別求得所有自變量之和的最大值和最小值，其結(jié)果值即氣象因子的優(yōu)化范圍。張強(qiáng)[1]在求解氣象因子范圍時，依次對每項品質(zhì)進(jìn)行單獨求解，然后對所有求得的氣象因子值進(jìn)行匯總來確定范圍，這種方法可能會導(dǎo)致求得的氣象因子只有一個最大值或者最小值，無法確定區(qū)間。

氣象因子在蘋果果實品質(zhì)形成中發(fā)揮著重要作用[17-20]。河北省地區(qū)周年的氣象因子均適宜‘富士果實品質(zhì)形成。生長期的日較差、降水量和濕度都在適宜范圍之內(nèi)。生長期的溫度和日照時數(shù)上限高于優(yōu)化后的氣象因子范圍，這主要是由于日照時數(shù)增加和溫度升高會導(dǎo)致裂果指數(shù)增加，從而降低果實的外觀品質(zhì)。蘋果種植的氣候條件沒有標(biāo)準(zhǔn)的限定范圍，本文中氣象和果實品質(zhì)數(shù)據(jù)結(jié)果可為‘富士蘋果的栽培提供一定參考。

本研究僅對‘富士蘋果品質(zhì)而言，其他栽培品種因限于觀測資料的局限性，本文未曾討論。本文雖只運用了一年的數(shù)據(jù)，但果園取樣遍布了河北省絕大部分的蘋果種植區(qū)，基本能說明該地區(qū)的氣象因子對果品影響。在氣象因素中選擇了全年和生長期共10項指標(biāo)，回歸方程的擬合值與實際值基本一致，而對果實某單一品質(zhì)影響的氣候因子指標(biāo)可能更少，回歸方程可能會更加簡化，如果實著色的影響光照條件關(guān)鍵在于轉(zhuǎn)色期的兩周之內(nèi)[21-22]，果實裂果可能在于某一日的極端高溫或日較差[23]，更加精準(zhǔn)影響因子的確定需要下一步試驗研究與驗證。

綜上所述，氣象因子對‘富士蘋果的著色指數(shù)、裂果指數(shù)和質(zhì)構(gòu)特性影響顯著，對果實的單果質(zhì)量、可溶性固形物含量和可滴定酸含量影響不顯著。河北省地區(qū)適宜‘富士蘋果生長的氣象條件范圍為生長期的日均氣溫17.05～18.42 ℃，平均日較差9.94～14.61 ℃，日均降水量1.29～3.09 mm，日均濕度55.01%～71.55%，日均日照時數(shù)5.72～5.92 h;年日均氣溫8.38～15.60 ℃，平均日較差9.59～14.26 ℃，日均降水量0.79～1.93 mm，日均濕度50.15%～65.67%，日均日照時數(shù)5.19～7.48 h。

參考文獻(xiàn)：

[1] 張強(qiáng). 蘋果果實品質(zhì)與土壤養(yǎng)分和氣象因子關(guān)系的研究[D]. 北京： 中國農(nóng)業(yè)大學(xué)， 2018.

[2] 聶繼云， 李靜， 楊振鋒， 等. 蘋果品質(zhì)和質(zhì)量安全問題與對策[J]. 中國果樹， 2007（3）： 60-62.

[3]生吉萍， 王健健， 冷傳祝. 蘋果及其加工制品的質(zhì)量安全及應(yīng)對策略[J]. 食品科學(xué)技術(shù)學(xué)報， 2015， 33（4）： 11-15.

[4] 楊振鋒， 叢佩華， 李靜， 等. 蘋果內(nèi)在品質(zhì)分析[J]. 中國果樹， 2002（2）： 16-17.

[5] 楊振鋒， 聶繼云， 李靜， 等. 富士蘋果理化品質(zhì)分析[J]. 北方果樹， 2007（1）： 11-12.

[6] SUGIURA T， OGAWA H， FUKUDA N， et al. Changes in the taste and textural attributes of apples in response to climate change[J]. Scientific Reports， 2013， 3（2）： 2418.

[7] 申順吏， 楊俊梅， 鞏在武. 基于灰關(guān)聯(lián)的山西蘋果產(chǎn)量氣候影響因子分析及蘋果產(chǎn)量預(yù)測[J]. 南方農(nóng)業(yè)學(xué)報， 2016， 47（7）：1146-1154.

[8] SITOMPUL S M， SITAMATI， SUGITO Y. Spatial productivity analysis of tropical apple （Malus sylvestris Mill） in relation to temperature with PCRaster[J]. Journal of Agricultural Science and Technology A， 2013， 3： 183-192.

[9] 童盼盼， 張亞若， 湯蕾， 等. 光照和溫度條件對富士蘋果糖心形成的影響[J]. 西北農(nóng)業(yè)學(xué)報， 2020， 29（4）： 579-586.

[10] 張磊， 張曉煜， 段曉鳳， 等. 氣象條件對高酸蘋果若干品質(zhì)因子的影響[J].氣象與環(huán)境科學(xué)， 2018， 41（3）： 51-56.

[11] 魏欽平， 張繼祥， 毛志泉， 等. 蘋果優(yōu)質(zhì)生產(chǎn)的最適氣象因子和氣候區(qū)劃[J]. 應(yīng)用生態(tài)學(xué)報， 2003， 14（5）： 713-716.

[12] 魏欽平， 李嘉瑞， 束懷瑞. 蘋果品質(zhì)與生態(tài)因子關(guān)系的研究[J]. 山東農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)版）， 1998， 29（4）： 532-536.

[13] 李興軍， 張光倫， 劉镕山. 生態(tài)因子對蘋果著紅色的生態(tài)效應(yīng)[J]. 果樹科學(xué)， 2000， 17（2）： 147-150.

[14] 段曉鳳， 張磊， 金飛， 等. 氣象因子對蘋果產(chǎn)量、品質(zhì)的影響研究進(jìn)展[J]. 中國農(nóng)學(xué)通報， 2014， 30（7）： 33-37.

[15] 楊梅， 肖靜， 蔡輝. 多元分析中的多重共線性及其處理方法[J]. 中國衛(wèi)生統(tǒng)計， 2012， 29（4）： 620-624.

[16] SHAH A， CHAUHAN Y， CHAUDHURY B. Principal component analysis based construction and evaluation of cryptocurrency index[J]. Expert Systems With Applications， 2021， 163： 113796.

[17] ORNELAS-PAZ J J， QUINTANA-GALLEGOS B M， ESCALANTE-MINAKATA P， et al. Relationship between the firmness of Golden Delicious apples and the physicochemical characteristics of the fruits and their pectin during development and ripening[J]. Journal of Food Science and Technology， 2018， 55（1）： 33-41.

[18] EL YAACOUBI A， EL JAOUHARI N， BOURIOUG M， et al. Potential vulnerability of Moroccan apple orchard to climate change-induced phenological perturbations： effects on yields and fruit quality[J]. International Journal of Biometeorology， 2020， 64（3）： 377-387.

[19] 孫小諾， 劉悅， 王蓉蓉， 等. 對河北省蘋果品質(zhì)與氣象條件關(guān)系的分析[J]. 現(xiàn)代農(nóng)業(yè)研究， 2019（2）： 33-34.

[20] 羅文文. 不同生境富士蘋果葉片和果實品質(zhì)差異及其與環(huán)境因子關(guān)系研究[D]. 楊凌： 西北農(nóng)林科技大學(xué)， 2014.

[21] 孫文泰， 尹曉寧， 劉興祿， 等. 不同海拔高度對“紅富士”蘋果果實品質(zhì)的影響[J]. 北方園藝， 2013（6）： 12-15.

[22] WERTHEIM S J， WAGENMAKERS P S， BOOTSMA J H， et al. Orchard systems for apple and pear： conditions for success[J]. Acta Horticulturae， 2001， 557： 209-227.

[23] 李肇林. 影響果實外觀的幾個因素[J]. 西北園藝： 綜合， 2019（1）： 34-36.