孫恩虹 張建平 張凱 童建川 蘇家奎 葉坤

摘要：【目的】研究關(guān)鍵生育期氣象條件對青花椒品質(zhì)的影響，并構(gòu)建青花椒氣候品質(zhì)評估模型，為青花椒氣候品質(zhì)評估及區(qū)域優(yōu)勢氣候資源開發(fā)利用等提供技術(shù)支撐?！痉椒ā恳郧嗷ń菲贩N九葉青為試驗材料，于2015—2020年分別在重慶地區(qū)11個花椒主產(chǎn)區(qū)選擇代表性花椒園3～5個取樣，測定青花椒品質(zhì)指標;采用相關(guān)分析法篩選出影響青花椒揮發(fā)油、麻味素和硒含量等品質(zhì)指標的主要氣象因子，建立氣象因子與青花椒品質(zhì)指標的回歸模型;利用模型推算出不同花椒品質(zhì)等級相對應(yīng)的氣象因子閾值，利用層次分析法構(gòu)建青花椒氣候品質(zhì)評估模型;以未參與建模的品質(zhì)檢測數(shù)據(jù)驗證模型?！窘Y(jié)果】影響青花椒氣候品質(zhì)優(yōu)劣的關(guān)鍵氣象因子為果實膨大期的日照時數(shù)（X1）、累計降水量（X2）、氣溫日較差（X3）、≥85%高濕日數(shù)（X4），采收期的日照時數(shù)（X5）、累計降水量（X6）、氣溫日較差（X7）和雨日數(shù)（X8），采用層次分析法確定評估模型中各氣象因子的權(quán)重系數(shù)分別為0.19、0.06、0.17、0.13、0.07、0.19、0.06和0.13，青花椒氣候品質(zhì)評估模型為：ICQ=0.19X1+0.06X2+0.17X3+0.13X4+0.07X5+0.19X6+0.06X7+0.13X8。青花椒氣候品質(zhì)劃分為特優(yōu)、優(yōu)、良和一般4個等級，評估指數(shù)（ICQ）相應(yīng)閾值為ICQ≥0.85、0.75≤ICQ<0.85、0.65≤ICQ<0.75、ICQ<0.65。基于評估模型，對2020年重慶江津幾江地區(qū)青花椒氣候品質(zhì)進行評估分析，評估指數(shù)為0.91，判定為特優(yōu)等級，與花椒園現(xiàn)場取樣檢測結(jié)果相符。【結(jié)論】研究建立的評估方法和指標能客觀反映青花椒的品質(zhì)等級，可應(yīng)用于青花椒生產(chǎn)的氣象條件適宜度服務(wù)和花椒采收期預(yù)報服務(wù)。

關(guān)鍵詞： 青花椒;氣候品質(zhì);氣象因子;評估模型

中圖分類號： S573;S162.5? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2021）07-1942-10

Establishment and validation of elevation model of climate quality for Zanthoxylum schinifolium Sieb.et Zucc

SUN En-hong1，2， ZHANG Jian-ping3*， ZHANG Kai1，TONG Jian-chuan4，

SU Jia-kui4， YE Kun1

（1Jiangjin Modern Agrometeorology Experimental Station of Chongqing， Chongqing? 402260， China;

2Jiangjin Meteorological Bureau of Chongqing，? Chongqing? 402260， China; 3Chongqing Institute

of Meteorological Sciences， Chongqing? 401147，? China; 4Jiangjin Agriculture and Rural

Affairs Committee of Chongqing， Chongqing? 402260， China）

Abstract：【Objective】The effects of meteorological conditions during the key growth period on quality of Zanthoxylum schinifolium Sieb.et Zucc were researched and the evaluation model of climate quality for Z. schinifolium was established， in order to provide technical support for climate quality evaluation for Z. schinifolium and exploitation and utilization of regional climate resources. 【Method】Took Jiuyeqing Z. schinifolium as material，samples were collected at major producing districts in Chongqing during 2015-2020. Three to five orchards were selected in each district. Quality indexes were tested for Z. schinifolium. Correlation analysis were used to screen climate factors that influenced volatile oil， amide compounds and selenium of Z. schinifolium. Regression model was established by climate factor and quality index. The threshold of climate factors corresponding to the Z. schinifolium quality and classification criteria for climatic conditions of Z. schinifolium were also determined by regression model. Analytical hierarchy process（AHP） were used to establish the evaluation model. To predict the effects of the evaluation model by using the other samples that were not included in the model development. 【Result】There were eight key climate factors that reflected the climate quality of Z. schinifolium， which were sunshine duration（X1）， cumulative precipitation（X2）， daily temperature range（X3）， ≥85% humidity days du-ring the fruit formation period（X4） and sunshine duration（X5）， cumulative precipitation（X6）， daily temperature range（X7）， rainy days（X8） during picking period. The weight coefficients of climate factors were 0.19， 0.06， 0.17， 0.13， 0.07， 0.19， 0.06，0.13，respectively. The climate quality evaluation model of Z. schinifolium was ICQ=0.19X1+0.06X2+0.17X3+0.13X4+0.07X5+0.19X6+0.06X7+0.13X8. The climate quality of Z. schinifolium were divided into 4 grades： extra excellent， excellent，good and normal. The thresholds of evaluation indexes（ICQ） wereICQ≥0.85，0.75≤ICQ<0.85，0.65≤ICQ<0.75 and ICQ<0.65. Based on the model，the comprehensive index of climate quality of Z. schinifolium picked in Jijiang，Jiangjin District，Chongqing in 2020 was 0.91，which was judged as extra excellent and consistent with actual conditions in sampling sites. 【Conclusion】The established evaluation methods and indicators can objectively reflect the quality grade of Z. schinifolium，they can be applied to serve in climate conditions suitability of Z. schinifolium production and harvest forecast of Z. schinifolium.

Key words： Zanthoxylum schinifolium Sieb.et Zucc; climate quality; climate factor; evaluation model

Foundation item： Innovation and Development Project of China Meterological Administration（CXFZ202 1J68）

0 引言

【研究意義】青花椒（Zanthoxylum schinifolium Sieb. et Zucc）屬蕓香科（Rutaceae）花椒屬（Zanthoxylum L.）落葉植物，因果實成熟后為青色而得名，是著名的香料和油料樹種，廣泛栽培于我國西南地區(qū)（姚佳和蒲彪，2010;程敏等，2016）。目前，西南地區(qū)主栽的九葉青花椒被全國多個地區(qū)引種，已成為鄉(xiāng)村振興的重要支柱產(chǎn)業(yè)之一。但花椒品質(zhì)優(yōu)劣與氣象條件密不可分，即便是同種花椒品種，由于其生長時的光照時間和降水情況不同，其揮發(fā)油和麻味素等品質(zhì)成分也不相同（韓若梅等，2017）。因此，開展青花椒氣候品質(zhì)研究對提高區(qū)域氣候資源利用，提升農(nóng)產(chǎn)品經(jīng)濟效益具有重要作用。【前人研究進展】目前，國內(nèi)已有大量學(xué)者采用多元回歸和層次分析等方法針對葡萄（張曉煜等，2007）、茶葉（金志鳳等，2015）、沙田柚（梁敏妍等，2019）、柑橘（謝遠玉等，2019）和油茶（江滌非等，2021）等特色經(jīng)濟作物進行了農(nóng)產(chǎn)品氣候品質(zhì)評估研究。在花椒品質(zhì)與氣候關(guān)系方面也有學(xué)者開展了基礎(chǔ)研究。李智杰等（2006）應(yīng)用因子分析法對全國10個花椒主產(chǎn)地的環(huán)境因子進行分析，發(fā)現(xiàn)花椒產(chǎn)地的溫度狀況和降水量等氣象條件對花椒品質(zhì)的影響較大。孫玉蓮等（2014）采用數(shù)理統(tǒng)計方法，利用歷史花椒觀測資料，研究甘肅臨夏地區(qū)影響花椒生長的關(guān)鍵生育期及主要氣象因子，建立臨夏地區(qū)花椒各生育期時空動態(tài)氣候預(yù)測模型。杜晉城等（2015）調(diào)查不同海拔靈山正路椒的物候期、品質(zhì)及產(chǎn)量變化，發(fā)現(xiàn)花椒物候期隨海拔升高而延遲，品質(zhì)和產(chǎn)量在1800～2600 m內(nèi)與海拔呈正相關(guān)。劉麗等（2017）基于花椒農(nóng)業(yè)氣象平行觀測試驗數(shù)據(jù)，確定了武都花椒生長關(guān)鍵時期和影響其生長的關(guān)鍵氣象要素，在確定適宜性區(qū)劃指標的基礎(chǔ)上，結(jié)合GIS技術(shù)建立武都山區(qū)花椒氣候適宜性區(qū)劃。廖要明和楊文平（2019）對鳳縣花椒生產(chǎn)的氣候條件和主要氣象災(zāi)害進行系統(tǒng)分析，構(gòu)建花椒生產(chǎn)的氣候影響評價指標體系，實現(xiàn)對區(qū)域花椒生產(chǎn)的氣候生態(tài)適應(yīng)性定量評估。張群等（2020）采用化學(xué)分析方法測定獅子頭、美鳳椒和秦安1號花椒葉片的黃酮含量，發(fā)現(xiàn)花椒葉總黃酮存在顯著的季節(jié)變化和年際變化，同時花椒葉總黃酮含量與日照百分率、降水量和空氣濕度顯著相關(guān)，與當(dāng)月的日照百分率呈負相關(guān)，與上月的降水量和當(dāng)月的空氣濕度呈正相關(guān)。可見，建立氣象條件與果實品質(zhì)指標的量化關(guān)系是開展氣候品質(zhì)評估工作的基礎(chǔ)?！颈狙芯壳腥朦c】目前對影響花椒產(chǎn)量和品質(zhì)的主要氣象條件的響應(yīng)關(guān)系研究雖有涉及，但研究甚少且深度不夠，特別是西南地區(qū)針對因年際間氣候條件差異引起的花椒品質(zhì)優(yōu)劣的定量評估鮮有報道?！緮M解決的關(guān)鍵問題】基于重慶青花椒主產(chǎn)區(qū)（縣）不同年份的青花椒品質(zhì)檢測資料和同期氣象資料，利用相關(guān)分析法、層次分析法和加權(quán)求和法研究青花椒主要品質(zhì)指標與氣象要素的關(guān)系，構(gòu)建青花椒氣候品質(zhì)評估模型，以期更好地為青花椒氣候品質(zhì)評估及區(qū)域優(yōu)勢氣候資源開發(fā)利用等提供技術(shù)支撐。

1 材料與方法

1. 1 研究資料

1. 1. 1 青花椒品質(zhì)資料 供試青花椒品種為九葉青，該品種最早發(fā)現(xiàn)于重慶市江津區(qū)先鋒鎮(zhèn)，是目前我國栽培面積最大的青花椒品種和最重要的花椒資源之一（徐潔，2007）。青花椒品質(zhì)生化指標主要包括揮發(fā)油和麻味素（余曉琴等，2009b;薛小輝和蒲彪，2013）等。其中，揮發(fā)油指青花椒果皮中具有芳香氣味、可隨水蒸氣蒸餾出來且與水不相混溶的揮發(fā)性油狀物質(zhì)的總和（唐課文等，2003），是青花椒香氣的主要成分;麻味素由一類不飽和脂肪族酰胺組成，是青花椒麻味的主要來源（劉雄等，2004）;此外，由于江津土壤富硒，所產(chǎn)青花椒果實中還富含對人體健康有益的硒元素，因此將青花椒硒含量作為一項重要的花椒品質(zhì)指標。

本研究所用青花椒品質(zhì)數(shù)據(jù)來自2015—2020年重慶地區(qū)11個花椒主產(chǎn)區(qū)（縣）（江津、酉陽、豐都、奉節(jié)、南川、銅梁、永川、璧山、合川、潼南和涪陵）的檢測資料。每個區(qū)（縣）選擇代表性花椒園3～5個，為保證樣本的均一性，同一椒園按東、南、西、北4個方位各選定花椒樹2株，共8株，每株花椒樹上選取不同方位的果穗，每個樣本3個重復(fù)。供試九葉青花椒樹樹齡均為5～6年，管理水平中上，平均產(chǎn)量11250～15000 kg/ha，屬當(dāng)?shù)卣．a(chǎn)量水平。青花椒揮發(fā)油含量參考《中華人民共和國藥典法》（一部）附錄中揮發(fā)油測定法（甲法）測定;麻味素含量參考DB 50/T 321—2009《花椒麻味物質(zhì)的檢測方法 高效液相色譜法》測定;硒含量參考GB 5009.93—2017《食品安全國家標準 食品中硒的測定》中的原子熒光光譜法測定。此外，青花椒氣候品質(zhì)等級標準還參考GB/T? 30391—2013《花椒》、GB/T 30385—2013《香辛料和調(diào)味品》、LY/T 1652—2005《花椒質(zhì)量等級》、LY/T? 2042—2012《九葉青花椒豐產(chǎn)栽培技術(shù)規(guī)程》、QX/T 486—2019《農(nóng)產(chǎn)品氣候品質(zhì)認證技術(shù)規(guī)范》和DBS 50/003—2014《食品安全地方標準 保鮮花椒》等6項標準。

1. 1. 2 氣象資料 氣象資料使用花椒樣本對應(yīng)的區(qū)域氣象站數(shù)據(jù)，包括2015—2020年花椒關(guān)鍵生育期內(nèi)的氣溫、降水量、空氣濕度、雨日數(shù)和日照時數(shù)。評估所用氣象資料為江津區(qū)2020年區(qū)域站氣象觀測資料。由于影響青花椒品質(zhì)的溫度、水分、光照和濕度等參與建模的氣象要素的量級不同，本研究參照常規(guī)農(nóng)業(yè)氣象條件定量化等級評價標準（孫秀邦等，2017），建模前對氣象數(shù)據(jù)進行無量綱化處理，將各氣象因子統(tǒng)一劃分為4個等級，分別賦予0.4～1.0的數(shù)值。4個等級賦值標準如下：

Mi=[1.0? ?Ti01≤Xi≤Ti020.8? ?Ti11≤XiTi22] （1）

式中，Mi為影響青花椒品質(zhì)的第i個氣候品質(zhì)指標;Xi為氣象因子的實測值;Ti01和Ti02分別為青花椒氣候品質(zhì)特優(yōu)時氣象因子的下限值和上限值;Ti11和Ti12分別為青花椒氣候品質(zhì)優(yōu)時氣象因子的下限值和上限值;Ti21和Ti22分別為青花椒氣候品質(zhì)良時氣象因子的下限值和上限值;如果氣象因子低于Ti21或高于Ti22，則青花椒氣候品質(zhì)一般。

1. 2 研究方法

1. 2. 1 青花椒氣候品質(zhì)等級劃分 結(jié)合區(qū)域氣候特征及青花椒的果實特性（趙志峰，2018;孫恩虹等，2019），參考LY/T 1652—2005《花椒質(zhì)量等級》中花椒品質(zhì)分級標準、DBS 50/003—2014《食品安全地方標準 保鮮花椒》、DB 50/T 705—2016《富硒農(nóng)產(chǎn)品》、相關(guān)文獻（余曉琴等，2009b）以及花椒生產(chǎn)相關(guān)的氣象和農(nóng)業(yè)專家依據(jù)實際生產(chǎn)中氣象因子對花椒生長的影響提出的青花椒分級意見（馬力文等，2018），將青花椒揮發(fā)油、麻味素和硒含量分為4個等級（表1）。

1. 2. 2 青花椒氣候品質(zhì)評估模型的構(gòu)建方法 每個區(qū)（縣）選擇7～8個樣本共78份樣本用于建立模型。首先，運用Pearson簡單相關(guān)系數(shù)（黃娟等，2018）將青花椒揮發(fā)油、麻味素和硒含量等檢測數(shù)據(jù)分別與青花椒品質(zhì)形成關(guān)鍵生育期即果實膨大期（4月上旬—5月下旬）和采收期（6月上旬—7月中旬）內(nèi)的日照時數(shù)、降水量、氣溫日較差、≥85%高濕日數(shù)和降水日數(shù)等氣象因子分別進行相關(guān)分析，篩選出通過顯著性檢驗的氣象因子。然后，分別研究揮發(fā)油、麻味素、硒含量與相關(guān)顯著的氣象因子的線性、非線性關(guān)系（馬力文等，2018），選擇通過R檢驗且對青花椒品質(zhì)形成具有生物學(xué)意義的氣象因子，建立不同品質(zhì)指標與氣象因子的關(guān)系模型，并利用模型推算出影響青花椒品質(zhì)各氣象因子的閾值，進而劃分青花椒氣候品質(zhì)等級標準?；谇嗷ń饭麑嵣L發(fā)育過程中各氣象因子對品質(zhì)的影響程度，利用層次分析法定量計算，確定各相關(guān)性顯著氣象因子的權(quán)重系數(shù)、果實膨大期氣候品質(zhì)指數(shù)（m1）和采收期氣候品質(zhì)指數(shù)（m2），分別由相應(yīng)關(guān)鍵品質(zhì)因素的加權(quán)和構(gòu)成，即式（2）：

mi=[m1=j=14ajyjm2=j=58bjyj] （2）

式中，aj、bj分別為果實膨大期氣候品質(zhì)指數(shù)和采收期氣候品質(zhì)指數(shù)各關(guān)鍵氣象因子的權(quán)重系數(shù)，由層次分析法（AHP）確定（劉遵春等，2006）。yj為各關(guān)鍵氣象因子。

最后，按照QX/T 486—2019《農(nóng)產(chǎn)品氣候品質(zhì)認證技術(shù)規(guī)范》，采用加權(quán)指數(shù)求和方法，建立青花椒氣候品質(zhì)評估模型為：

ICQ=[i=1naimi] （3）

式中，ICQ為氣候品質(zhì)評估指數(shù);n為氣候品質(zhì)評估指標的數(shù)量;αi為氣候品質(zhì)評估指標的權(quán)重。

根據(jù)氣候品質(zhì)指數(shù)計算結(jié)果，結(jié)合青花椒生產(chǎn)實際，將青花椒氣候品質(zhì)評估標準統(tǒng)一劃分為4個等級，分別為特優(yōu)（1級）、優(yōu)（2級）、良（3級）和一般（4級）。

1. 2. 3 青花椒氣候品質(zhì)評估模型檢驗 為檢驗該評估方法的適應(yīng)性，在青花椒主產(chǎn)區(qū)（縣）各選擇2～3個未參與建模的樣本，共29個樣本用于模型檢驗。對花椒檢測樣本及對應(yīng)氣象站點的氣象數(shù)據(jù)進行評估檢驗，將計算的氣候品質(zhì)評估指數(shù)按照分類閾值進行品質(zhì)等級歸類，并與各年份對應(yīng)取樣地點的實測值對應(yīng)的實際等級對比驗證其合理性。

1. 3 統(tǒng)計分析

利用SAS v8對數(shù)據(jù)進行相關(guān)分析和回歸分析，采用yaahp進行層次分析，用SigmaPlot 14.0作圖。

2 結(jié)果與分析

2. 1 青花椒品質(zhì)指標

由圖1可知，青花椒揮發(fā)油含量在5.6～7.1 mL/100 g，可見，青花椒揮發(fā)油指標高于優(yōu)等標準;麻味素含量為14.8～18.8 mg/g，高于四川和云南等地的青花椒（余曉琴等，2009b）;依據(jù)DB 50/T 705—2016《富硒農(nóng)產(chǎn)品》硒含量在0.020～1.000 mg/kg有益于人體健康，2015—2020年青花椒果實中硒含量為0.050～0.120 mg/kg，硒含量符合富硒農(nóng)產(chǎn)品相關(guān)標準。青花椒品質(zhì)數(shù)據(jù)分析結(jié)果表明，青花椒風(fēng)味和營養(yǎng)俱佳，與實際情況相吻合。但由于年際間氣象條件差別較大，導(dǎo)致青花椒品質(zhì)年際間存在明顯差異，2020年揮發(fā)油含量最高，2018年麻味素含量最高，2017年硒含量最高。

2. 2 影響青花椒品質(zhì)關(guān)鍵氣象因子的選取

利用青花椒揮發(fā)油、麻味素、硒含量品質(zhì)檢測資料和各年份對應(yīng)站點氣象資料進行相關(guān)分析，篩選出相關(guān)性通過0.05顯著性檢驗的8個氣象因子，其中青花椒果實膨大期累計降水量對揮發(fā)油和硒含量均有顯著（P<0.05，下同）或極顯著（P<0.01，下同）影響。如表2表示，青花椒揮發(fā)油含量與果實膨大期的日照時數(shù)、累計降水量和氣溫日較差呈極顯著或顯著正相關(guān)，即果實膨大期充足的日照、充沛的降水與較大的氣溫日較差均有利于揮發(fā)油的積累;揮發(fā)油含量與采收期的累計降水量呈極顯著負相關(guān)，表明過多的雨水會導(dǎo)致青花椒含水量過高，且降水天氣影響青花椒采收，導(dǎo)致油腺破裂，色澤變暗，嚴重影響花椒品質(zhì)。青花椒麻味素含量與果實膨大期≥85%高濕日數(shù)呈極顯著負相關(guān)，高濕日數(shù)多會導(dǎo)致青花椒銹病和斑點病等病害發(fā)生發(fā)展，影響青花椒外觀及麻味素等物質(zhì)含量，與生產(chǎn)實際相符;此外，麻味素含量與采收期的雨日數(shù)呈極顯著負相關(guān)，與采收期的日照時數(shù)和氣溫日較差呈極顯著正相關(guān)，采收前期青花椒未完全成熟，充足的光熱資源可促進麻味素等物質(zhì)的進一步積累。青花椒硒含量與果實膨大期的累計降水量呈極顯著負相關(guān)，即青花椒硒含量積累主要在果實膨大期，降水過多可能影響青花椒對硒的吸收，與潘文杰等（2006）的研究結(jié)果一致。植物對礦質(zhì)元素的需求量與其生物量密切相關(guān)，且植物對硒的吸收有一定飽和現(xiàn)象（管文文等，2018），花椒果實膨大期果實快速膨大，營養(yǎng)物質(zhì)大量積累，所以采收期植物硒的吸收變化不大。因此，關(guān)鍵生育期內(nèi)充足的光照、較大的氣溫日較差、適宜的降水，以及較少的≥85%高濕日數(shù)和雨日數(shù)對青花椒最終優(yōu)良品質(zhì)的形成較有利。

2. 3 青花椒各品質(zhì)與氣象條件的關(guān)系

2. 3. 1 揮發(fā)油與氣象條件的關(guān)系 將各取樣點揮發(fā)油檢測數(shù)據(jù)與其果實膨大期的日照時數(shù)、累計降水量、氣溫日較差以及采收期的累計降水量等通過顯著性檢驗的氣象因子進行回歸分析，結(jié)果見表3。根據(jù)擬合方程，揮發(fā)油含量與果實膨大期的日照時數(shù)和氣溫日較差均呈線性關(guān)系，日照時數(shù)增加50 h，揮發(fā)油含量約增加1.20 mL/100 g，日較差提高1 ℃，揮發(fā)油含量增加0.97 mL/100 g。揮發(fā)油含量與果實膨大期的累計降水量呈對數(shù)關(guān)系，揮發(fā)油含量隨累計降水量的增加而升高，當(dāng)降水量由200 mm增加到300 mm時，揮發(fā)油含量增加1.16 mL/100 g。青花椒揮發(fā)油含量與采收期的累計降水量呈負線性關(guān)系，降水量增加100 mL，揮發(fā)油含量降低1.80 mL/100 g。因此，揮發(fā)油的積累需要果實膨大期有充足的溫光水資源，增強青花椒光合作用，提高光合產(chǎn)物凈積累，采收期的降水量不宜過多。

2. 3. 2 麻味素與氣象條件的關(guān)系 將各取樣點麻味素檢測數(shù)據(jù)與其果實膨大期的≥85%高濕日數(shù)，采收期日照時數(shù)、氣溫日較差和雨日數(shù)等通過顯著性檢驗的氣象因子進行回歸分析，結(jié)果見表4。由表4可看出，麻味素含量與果實膨大期的≥85%高濕日數(shù)呈負線性關(guān)系，≥85%高濕日數(shù)每增加10 d，麻味素含量約降低3.16 mg/g。青花椒麻味素含量與采收期的日照時數(shù)和氣溫日較差均呈線性關(guān)系，日照時數(shù)增加20 h，麻味素含量約增加1 mg/g;氣溫日較差提高1 ℃，麻味素含量增加2.491 mg/g。青花椒麻味素含量與采收期的雨日數(shù)呈負線性關(guān)系，雨日數(shù)增加2 d，麻味素含量降低1 mg/g。因此，青花椒采收期更需要晴好少雨天氣。

2. 3. 3 硒與氣象條件的關(guān)系 選擇與硒含量極顯著相關(guān)的果實膨大期累計降水量與硒含量進行回歸分析，結(jié)果見表5。由表5可看出，青花椒果實中的硒含量與果實膨大期的累計降水量呈對數(shù)關(guān)系，當(dāng)降水量由359.30 mm增加到391.69 mm時，硒含量減少0.015 mg/kg，即果實膨大期雨水過多不利于青花椒硒含量的積累和轉(zhuǎn)化。

2. 4 青花椒氣候品質(zhì)等級劃分

2. 4. 1 青花椒氣候品質(zhì)評估模型中氣象指標的分級賦值 根據(jù)青花椒品質(zhì)等級標準（表1），利用青花椒揮發(fā)油、麻味素和硒含量等品質(zhì)指標與氣象因素的回歸方程（表3～表5），反推出與青花椒品質(zhì)相對應(yīng)的氣象因子閾值（表6）。依據(jù)表6中氣象指標的等級，分別賦予0.4～1.0的數(shù)值，將氣象數(shù)據(jù)進行歸一化處理。由表6可見，優(yōu)等級的青花椒一般要求果實膨大期氣溫日較差在7.7 ℃以上為佳，日照時數(shù)在147.0 h以上，利于揮發(fā)油和麻味素積累，品質(zhì)更好。此時，葉片光合作用旺盛，葉面水分蒸發(fā)量大，需要充足的水分;相反，雨水過多，日照過少，濕度過大，則會導(dǎo)致硒積累量減少，麻味降低，因此該階段高濕日數(shù)少于23 d，累計降水量126.90～391.69 mm，更有利于青花椒優(yōu)良品質(zhì)形成。青花椒一般在6月初陸續(xù)進入采收期，采收前期麻味素含量得到進一步積累，保障青花椒優(yōu)良品質(zhì)的采收期氣候條件要求是多晴朗天氣，日照時數(shù)大于等于142.0 h，雨日少于30 d，晝夜溫差大于6.7 ℃。

2. 4. 2 氣候品質(zhì)構(gòu)成因素權(quán)重系數(shù)的確定

2. 4. 2. 1 氣候品質(zhì)綜合指數(shù)的權(quán)重系數(shù) 由相關(guān)分析篩選出影響青花椒品質(zhì)的8個氣象因子，結(jié)合生產(chǎn)實踐構(gòu)造層次結(jié)構(gòu)關(guān)系。根據(jù)青花椒生長發(fā)育特性，即品質(zhì)指標主要在果實膨大期和采收期形成，再結(jié)合專家經(jīng)驗（張玉榮等，2021），利用yaahp層次分析法，確定花椒果實膨大期和采收期氣象指數(shù)的權(quán)重分別為0.55和0.45，如式（4）所示。

ci=[0.55、0.45]? ?（4）

2. 4. 2. 2 氣候品質(zhì)指數(shù)構(gòu)成因素的權(quán)重系數(shù) 根據(jù)關(guān)鍵生育期各氣象要素對青花椒揮發(fā)油、麻味素和硒含量的影響程度，結(jié)合專家知識，對果實膨大期的4個氣象因素構(gòu)造判斷矩陣。模型判斷矩陣的一致性檢驗結(jié)果<0.1，通過一致性檢驗，說明本研究所構(gòu)造的判斷矩陣中各影響因子的相互關(guān)系較一致。經(jīng)計算，果實膨大期日照時數(shù)、累計降水量、氣溫日較差和≥85%高濕日數(shù)的權(quán)重向量如式（5）。

aj=[0.35、0.11、0.31、0.24]? （5）

同理，構(gòu)造采收期氣象因素判斷矩陣，模型判斷矩陣的一致性檢驗結(jié)果<0.1，通過一致性檢驗。經(jīng)計算，果實采收期日照時數(shù)、累計降水量、氣溫日較差和雨日數(shù)的權(quán)重向量如式（6）。

bj=[0.16、0.42、0.13、0.29] （6）

綜合式（4）、（5）和（6），得到青花椒氣候品質(zhì)評估模型為：ICQ=0.19X1+0.06X2+0.17X3+0.13X4+0.07X5+0.19X6+0.06X7+0.13X8? ? ? ? ? ? ? ? ? ? ? ? ? ? （7）

式中，X1～X4分別為花椒果實膨大期的日照時數(shù)、累計降水量、氣溫日較差和≥85%高濕日數(shù)，X5～X8分別為花椒采收期的日照時數(shù)、累計降水量、氣溫日較差和雨日數(shù)。

2. 5 青花椒氣候品質(zhì)評估模型的驗證

為檢驗?zāi)Ｐ偷脑u估效果，選用2015—2020年未參與建模的29份青花椒檢測樣本，按照青花椒品質(zhì)分類閾值，分別將實測的每項品質(zhì)要素按照4個等級歸類（表1）。同時，對青花椒樣本對應(yīng)站點的氣象因子進行分級賦值（表6），將氣象數(shù)據(jù)進行歸一化處理;再利用模型估算青花椒氣候品質(zhì)評估指數(shù)，并進行氣候品質(zhì)評估指數(shù)等級分類（表7）。結(jié)果表明，揮發(fā)油、麻味素和硒含量實測結(jié)果的等級劃分與氣候品質(zhì)評估等級變化趨勢基本一致，有2個樣本將2級估算為1級，有1個樣本將1級估算為2級（圖2），評估模型預(yù)測的準確率為89.7%?？傮w來看，該評估方法和指標能客觀反映青花椒的品質(zhì)等級，對品質(zhì)等級的評定效果較好。

2. 6 青花椒氣候品質(zhì)模型的實際應(yīng)用效果分析

應(yīng)用構(gòu)建的青花椒氣候品質(zhì)評估模型，結(jié)合區(qū)域站2020年4月6日—7月15日的逐日氣溫（最高、最低氣溫）、日照時數(shù)、空氣相對濕度、降水量和雨日數(shù)等觀測資料，分析評估2020年青花椒氣候品質(zhì)評估指數(shù)。以花椒種植區(qū)位上有代表性的江津幾江地區(qū)為例，用江津幾江站2020年氣象觀測資料計算，得到江津幾江地區(qū)青花椒氣象條件適宜度評估指標，結(jié)果見表8。

將表8各氣候品質(zhì)指標代入氣候品質(zhì)評估模型，得到江津幾江地區(qū)2020年青花椒氣候品質(zhì)評估指數(shù)（ICQ）為：ICQ=1.0×0.19+1.0×0.06+1.0×0.17+1.0×0.13+0.8×0.07+0.8×0.19+0.8×0.06+0.8×0.13=0.91。? 同理，分別計算江津其他種植區(qū)域的青花椒氣候品質(zhì)評估指數(shù)，經(jīng)計算得到氣候品質(zhì)評估指數(shù)的平均值為0.91，根據(jù)青花椒氣候品質(zhì)評估指數(shù)等級劃分標準，得到2020年青花椒氣候品質(zhì)等級為特優(yōu)。2020年4月6日—5月31日為青花椒果實膨大期，氣溫日較差為8.7～11.0 ℃，較大的晝夜溫差有利于花椒果實的快速膨大及揮發(fā)油、麻味素和硒的積累與轉(zhuǎn)化;日照時數(shù)274.1 h，較充足，對果實的膨大和品質(zhì)形成較有利;累計降水量171.5～333.1 mm，≥85%高濕日數(shù)8 d，雨量充沛但并不過多，雨日數(shù)適中，總體有利于花椒果實品質(zhì)的形成，但不會導(dǎo)致花椒銹病和斑點病等病害發(fā)生，從而影響花椒果實外觀及揮發(fā)油和麻味素的含量。2020年6月1日—7月15日為采收期，日照時數(shù)137.2 h，氣溫日較差6.4～8.8 ℃，對花椒光合產(chǎn)物后期積累無不利影響;累計降水量352.6～635.1 mm，雨日數(shù)24～31 d，降水量偏多，對花椒的采收略有影響，但對青花椒果實品質(zhì)的形成無明顯影響。綜合上述分析結(jié)果，結(jié)合氣象資料、現(xiàn)場實地調(diào)查和企業(yè)品質(zhì)檢查結(jié)果，整體認為2020年青花椒生長期間的氣候條件有利于其品質(zhì)的形成，計算結(jié)果與實際情況較一致。

3 討論

花椒果實生長發(fā)育過程中，氣溫、日照時數(shù)和空氣濕度等氣象要素對其生長發(fā)育具有顯著影響（樊惠芳等，2004;張奇等，2018）。分析重慶青花椒主產(chǎn)區(qū)連續(xù)6年的青花椒品質(zhì)數(shù)據(jù)與對應(yīng)氣象要素的相關(guān)性，發(fā)現(xiàn)青花椒不同品質(zhì)指標受氣象要素影響的關(guān)鍵時期不同，其中，揮發(fā)油在果實膨大期達到顯著影響的氣象因子多于采收期，表明前期氣象條件對揮發(fā)油影響更大，溫光資源不足和水分過少會影響花椒揮發(fā)油的積累和轉(zhuǎn)化;麻味素在采收期達到顯著影響的氣象因子多于果實膨大期，特別是日照時數(shù)和氣溫日較差對花椒麻味素形成仍有一定影響，與余曉琴等（2009a）研究發(fā)現(xiàn)花椒麻味素含量的峰值出現(xiàn)時間晚于揮發(fā)油的結(jié)論相符，即花椒麻味素在采收期還要完成進一步的積累，后期陰雨寡照等氣象條件會影響麻味素含量;硒含量的關(guān)鍵期是果實膨大期，與潘文杰等（2016）研究發(fā)現(xiàn)作物硒含量積累主要集中在生長發(fā)育中期的結(jié)論一致。

不同氣象要素對花椒果實品質(zhì)影響程度不同，其中青花椒果實膨大期的日照時數(shù)和氣溫日較差，采收期累計降水量對青花椒品質(zhì)影響最大，其權(quán)重分別為0.19、0.17和0.19，其次是果實膨大期的高濕日數(shù)、采收期的雨日數(shù)，其權(quán)重系數(shù)均為0.13;果實膨大期的累計降水量，采收期的日照時數(shù)和氣溫日較差對青花椒品質(zhì)的影響相對較小，其權(quán)重分別為0.06、0.07和0.06。在果實膨大期，不僅光照和氣溫日較差影響花椒品質(zhì)形成，而且降水的多寡和≥85%高濕日數(shù)也會影響花椒營養(yǎng)物質(zhì)的積累，降水過少不利于花椒光合作用形成揮發(fā)油;相反降水過多，影響花椒根系活力且土壤有效硒有一定損失，不利于花椒對硒的吸收和積累（潘文杰等，2016）。

本研究基于花椒氣候品質(zhì)指標，采取易于獲取資料的氣象要素，利用層次分析法和加權(quán)求和法，同時結(jié)合專家經(jīng)驗與花椒生產(chǎn)實際，從而構(gòu)建青花椒氣候品質(zhì)評估模型，可反映年際間氣候變化對花椒品質(zhì)產(chǎn)生的影響，突破了傳統(tǒng)定性描述花椒品質(zhì)的方法，實現(xiàn)了從氣象學(xué)角度對花椒品質(zhì)等級的量化評估，可為年份花椒評級提供參考依據(jù)，有助于提升花椒氣候品質(zhì)和品牌內(nèi)涵。青花椒氣候品質(zhì)評估模型是結(jié)合青花椒生長的地理氣候特點和實際農(nóng)業(yè)生產(chǎn)需求而設(shè)計，引入了花椒果實硒含量作為花椒品質(zhì)指標之一，在作物氣候品質(zhì)評估方面具有一定的創(chuàng)新性，應(yīng)在實踐中進一步檢驗、修正、補充和完善。

4 結(jié)論

氣候品質(zhì)評估結(jié)果將青花椒劃分為特優(yōu)、優(yōu)、良好、一般4個等級，符合當(dāng)?shù)厣a(chǎn)實際，對椒農(nóng)提高青花椒的商品價值具有積極意義。同時，明確了適宜花椒品質(zhì)形成各氣象指標的適宜度值，研究結(jié)果不僅可應(yīng)用于青花椒氣候品質(zhì)的評估，還可進一步開展研究訂正指標，運用于青花椒生產(chǎn)的氣象條件適宜度服務(wù)和花椒采收期預(yù)報服務(wù)。

參考文獻：

程敏，高陽華，楊世琦，葉勤玉，齊會娟，趙彤. 2016. 重慶市花椒氣候生態(tài)區(qū)劃研究[J]. 西南師范大學(xué)學(xué)報（自然科學(xué)版），41（6）：53-59. doi：10.13718/j.cnki.xsxb.2016.06.010. [Cheng M，Gao Y H，Yang S Q，Ye Q Y，Qi H J，Zhao T. 2016. On climatic ecological regionalization of Zanthoxylum bungeanum in Chonging[J]. Journal of Southwest China Normal University（Natural Science Edition），41（6）： 53-59.]

杜晉城，代仕高，王澤亮，李德榮，馬文寶，陳炙，羅成榮. 2015. 不同海拔高度對靈山正路椒產(chǎn)量和品質(zhì)的影響[J]. 經(jīng)濟林研究，33（2）： 97-100. doi： l0.14067/j.cnki. 1003-8981.2015.02.017. [Du J C，Dai S G，Wang Z L，Li D R，Ma W B，Chen Z，Luo C Y. 2015. Effects of diffe-rent altitudes on yield and quality of Zanthoxylum bungeanum ‘Lingshan[J]. Nonwood Forest Research，33（2）： 97-100.]

樊惠芳，郭旭新，程斌，王穩(wěn)江，徐鳳蘭. 2004. 渭北旱原花椒產(chǎn)量與水分關(guān)系的研究[J]. 西北林學(xué)院學(xué)報，19（4）：32-34. [Fan H F，Guo X X，Cheng B，Wang W J，Xu F L. 2004. Studies on the relationship between the yield and moisture of Chinese Prickly Ash in Weibei arid area[J]. Journal of Northwest Forestry University，19（4）：32-34.]

管文文，戴其根，張洪程，尹雪斌. 2018. 水稻不同生育期對硒的吸收、轉(zhuǎn)運及累積規(guī)律[J]. 土壤，50（6）： 1150-1154. doi：10.13758/j.cnki.tr.2018.06.014. [Guan W W，Dai Q G，Zhang H C，Yin X B. 2018. Accumulation and distribution of selenium in rice （Oryza sativa） at different growth stages[J]. Soil，50（6）：1150-1154.]

韓若梅，何盒，申怡飛，張丁玲，劉淑明. 2017. 不同種源大紅袍花椒種實特性及主要化學(xué)成分的地理變化[J]. 西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版），45（6）：69-82. doi： 10.13207/j.cnki.jnwafu.2017.06.010. [Han R M，He H，Shen Y F，Zhang D L，Liu S M. 2017. Geographical varia-tions of seed characteristics and main chemical components of Zanthoxylum bungeanum with different provenances[J]. Journal of Northwest A & F University （Natural Science Edition），45（6）： 69-82.]

黃娟，李新建，吳新國，王雪姣，張仕明，范金霞. 2018. 庫爾勒香梨氣候品質(zhì)評價指標及模型的研究[J]. 沙漠與綠洲氣象，12（3）：87-94. doi：10.12057/j.issn.1002-0799.2018. 03.0013. [Huang J，Li X J，Wu X G，Wang X J，Zhang S M，F(xiàn)an J X. 2018. Research on evaluation index and model of the climatological quality for the korla fragrant pear[J]. Desert and Oasis Meteorology，12（3）：87-94.]

江滌非，郭海峰，卿燃莉，李矛，賈海鷹. 2021. 邵陽縣油茶種植氣候適宜性評價[J]. 江西農(nóng)業(yè)學(xué)報，33（1）：76-81. doi：10.19386 /j.cnki.jxnyxb.2021.01.14. [Jiang D F，Guo H F，Qing R L，Li M，Jia H Y. 2021. Climate suitability evaluation of Camellia in Shaoyang County[J]. Acta Agriculturae Jiangxi，33（1）：76-81.]

金志鳳，王治海，姚益平，李仁忠，王岳飛，盧建偉，葉建剛，徐平. 2015. 浙江省茶葉氣候品質(zhì)等級評價[J]. 生態(tài)學(xué)雜志，34（5）： 1456-1463. [Jin Z F，Wang Z H，Yao Y P，Li R Z，Wang Y F，Lu J W，Ye J G，Xu P. 2015. Evaluation of tea climate quality grade in Zhejiang[J]. Chinese Journal of Ecology，34（5）： 1456-1463.]

李智杰，孔德詳，鄧振義，段宏斌，廉坡. 2006. 不同產(chǎn)地花椒多指標綜合排序及環(huán)境因子影響分析[J]. 西北林學(xué)院學(xué)報，21（5）：74-77. doi：10.3969/j.issn.1001-7461.2006.05. 017. [Li Z J，Sun D X，Deng Z Y，Duan H B，Lian P. 2006. The analysis on the comparison of comprehensive indexes about seeds of bunge prickly ash produced in di-fferent areas and environmental factors[J]. Journal of Northwest Forestry University，21（5）：74-77.]

梁敏妍，黃翠銀，張羽，羅曉玲，王華，史麗，杜堯東. 2019. 廣東仁化縣與廣西容縣沙田柚氣候生態(tài)適應(yīng)性比較[J]. 南方農(nóng)業(yè)學(xué)報，50（11）：2496-2503. doi：10.3969/j.issn. 2095-1191.2019.11.16. [Liang M Y，Huang C Y，Zhang Y，Luo X L，Wang H，Shi L，Du Y D. Climatic and ecological adaptability of Shatian pomelo：A comparative analysis between Renhua of Guangdong and Rongxian of Guangxi[J]. Journal of Southern Agriculture，50（11）：2496-2503.]

廖要明，楊文平. 2019. 鳳縣花椒生產(chǎn)氣候條件和災(zāi)害分析[J]. 中國農(nóng)學(xué)通報，35（23）：84-90. doi：10.11924/j.issn. 1000-6850.casb19030104. [Liao Y M，Yang W P. 2019. Analysis of climate condition and disaster for Chinese prickly ash production in Fengxian County[J]. Chinese Agricultural Science Bulletin ，35（23）：84-90.]

劉麗，張燕，孫暢，蒲靜. 2017. 甘肅武都花椒生長氣象條件分析及氣候適應(yīng)性區(qū)劃研究[J]. 中國農(nóng)學(xué)通報，33（19）： 126-132. doi：10.11924/j.issn.1000-6850.casb17030121. [Liu L，Zhang Y，Sun C，Pu J. 2017. Pepper in Wudu district，Gansu： Meteorological condition analysis and climate adaptability regionalization[J]. Chinese Agricultural Science Bulletin，33（19）： 126-132.]

劉雄，闞建全，付陳梅，陳宗道，聶勛杰. 2004. 花椒麻味成分的提取與分離技術(shù)[J]. 食品與發(fā)酵工業(yè)，30（9）：112-116. doi：10.3321/j.issn：0253-990X.2004.09.028. [Liu X，Kan J Q，F(xiàn)u C M，Chen Z D，Nie X J. 2004. Study on extracting and separating spicy components from Zanthorylum bungeanum[J]. Food and Fermentation Industries，30（9）：112-116.]

劉遵春，包東娥，廖明安. 2006. 層次分析法在金花梨果實品質(zhì)評價上的應(yīng)用[J]. 西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版），34（8）： 125-128. doi：10.3321/j.issn：1671-9387. 2006.08.025. [Liu Z C，Bao D E，Liao M A. 2006. Application of analytic hierarchy process in evaluating Jinhua pear quality[J]. Journal of Northwest Sci-Tech University of Agriculture And Forestry（Natural Science Edition），34（8）： 125-128.]

馬力文，李劍萍，韓穎娟，李萬春. 2018. 賀蘭山東麓‘赤霞珠品質(zhì)形成氣象條件與評級方法研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報，26（3）：453-466. doi： 10.13930/j.cnki.cjea.170579. [Ma L W，Li J P，Han Y J，Li W C. 2018. Meteorological conditions and rating method of quality formation of ‘Cabernet Sauvignon grape in eastern foothills of Helan Mountain[J]. Chinese Journal of Eco-Agriculture，26（3）： 453-466.]

潘文杰，黃建國，姜超. 2006. 烤煙硒積累及其與土壤和氣候關(guān)系[J]. 植物營養(yǎng)與肥料學(xué)報，12（2）：260-265. doi：10. 11674/zwyf.2006.0220. [Pan W J，Huang J G，Jiang C. 2006. Se accumulation in flue-cured tobacco related to soil and climate[J]. Plant Nutrition and Fertilizer Science，12（2）： 260-265.]

孫恩虹，趙潔，韓旭，何永坤. 2019. 不同采收期江津青花椒氣候品質(zhì)變化[J]. 農(nóng)學(xué)學(xué)報，9（4）： 88-91. doi：10.11923/j.issn.2095-4050.cjas18020011. [Sun E H，Zhao J，Han X，He Y K. 2019. Zanthoxylum schinifolium Sieb.et Zucc in Jiangjin：Climate quality change with harvesting time[J]. Journal of Agriculture，9（4）： 88-91.]

孫秀邦，陶曙華，田青，曹樹新，劉瑞娜，王周青，陳欣鑫. 2017. 安徽涇縣春茶開采期預(yù)報及氣候品質(zhì)評價[J]. 中國農(nóng)學(xué)通報，33（29）：49-54. [Sun X B，Tao S H，Tian Q，Cao S X，Liu R N，Wang Z Q，Chen X X. 2017. Forecast of first plucking date of spring tea and evaluation of tea climate quality grade in Jingxian，Anhui[J]. Chinese Agri-cultural Science Bulletin，33（29）：49-54.]

孫玉蓮，邊學(xué)軍，韋伯龍，劉洪，張錦泉. 2014. 甘肅臨夏地區(qū)花椒生態(tài)氣候適應(yīng)性分析與產(chǎn)量動態(tài)模型[J]. 西南農(nóng)業(yè)學(xué)報，27（2）： 846-848. [Sun Y L，Bian X J，Wei B L，Liu H，Zhang J Q. 2014. Analysis of eco-climate suitabi-lity of prickly ash and climatic model of yield dynamic in Gansu Linxia[J]. Southwest China Journal of Agricultural Sciences，27（2）： 846-848.]

唐課文，閻建輝，許友，王志剛. 2003. 花椒揮發(fā)油的化學(xué)成分分析[J]. 分析化學(xué)，31（3）：384. doi：10.3321/j.issn：0253-3820.2003.03.037. [Tan K W，Yan J H，Xu Y，Wang Z G. 2003. Analysis of chemical components of volatile oil from Zanthoxylum bungeanum Maxin[J]. Chinese Journal of Analytical Chemis，31（3）：384.]

謝遠玉，王培娟，朱凌金，陳星，黃瑩. 2019. 基于氣象因子的贛南臍橙氣候品質(zhì)指標評價模型[J]. 生態(tài)學(xué)雜志，38（7）：2265-2274. doi：10.13292/j.1000-4890.201907.028. [Xie Y Y，Wang P J，Zhu L J，Chen X，Huang Y. 2019. Climate quality evaluation model for navel orange in Ganzhou[J]. Chinese Journal of Ecology，38（7）：2265-2274.]

徐潔. 2007. 九葉青花椒的起源及遺傳多樣性研究[D]. 重慶： 西南大學(xué). [Xu J. 2007. Studies on the origin and genetic diversity in Zanthoxylum armatum v. novemfolius[D]. Chongqing： Southwest University.]

薛小輝，蒲彪. 2013. 花椒風(fēng)味成分研究與產(chǎn)品開發(fā)現(xiàn)狀[J]. 核農(nóng)學(xué)報，27（11）：1724-1728. doi：10.11869/hnxb.2013. 11.1724. [Xue X H，Pu B. 2013. Research and development situation of Zanthoxylum bungeanum flavor component and products[J]. Journal of Nuclear Agricultural Scien-ces，27（11）：1724-1728.]

姚佳，蒲彪. 2010. 青花椒的研究進展[J]. 中國調(diào)味品，35（6）：35-39. doi：10.3969/j.issn.1000-9973.2010.06.005. [Yao J，Pu B. 2010. Research progress on Zanthoxylum schinifolium Sieb. et Zucc[J]. China Condiment，35（6）：35-39.]

余曉琴，吳素蕊，闞建全，陳宗道. 2009a. 重慶江津花椒不同采收時期的品質(zhì)變化[J]. 食品與發(fā)酵工業(yè)，（35）： 164-167. [Yu X Q，Wu S R，Kan J Q，Chen Z D. 2009a. Study on the quality changes of Jiangjin Zanthoxylum Schinifolium Sieb.et Zucc in different periods[J]. Food and Fermentation Industries，（35）： 164-167.]

余曉琴，鄭顯義，闞建全，郭靜. 2009b. 紅花椒和青花椒主要品質(zhì)特征指標值的評價[J].食品科學(xué)，30（15）： 45-48. doi：10.3321/j.issn：1002-6630.2009.15.009. [Yu X Q，Zheng X Y，Kan J Q，Guo J. 2009b. Evaluation of speci-fic quality of Zanthoxylum bungeanum Maxim and Zanthoxylum schinifolium Sieb. et Zucc[J]. Food Science，30（15）： 45-48.]

張群，陳紅林，鄭濤，蘇柯星，張聰珍，劉淑明. 2020. 花椒葉黃酮含量的變化及其與氣象因子的關(guān)系研究[J]. 西北林學(xué)院學(xué)報，35（3）： 43-47. doi：10.3969j.issn.1001-7461. 2020.03.07. [Zhang Q，Chen H L，Zheng T，Su K X，Zhang C Z，Liu S M. 2020. Seasonal variation of flavonoid content in the leaves of Zanthoxylum bungeanum and its relationship with meteorological factors[J]. Journal of Northwest Forestry University，35（3）： 43-47.]

張奇，任利明，梁宏喆，王海，劉淑明. 2018. 氣象條件對不同品種花椒果實生長發(fā)育的影響[J]. 西北林學(xué)院學(xué)報，33（5）：111-115. doi：10.3969/j.issn.1001-7461.2018.05.17. [Zhang Q，Ren L M，Liang H Z，Wang H，Liu S M. 2018. Effects of meteorological conditions on the growth of different varieties of pepper fruit[J]. Journal of Northwest Forestry University，33（5）：111-115.]

張曉煜，亢艷莉，袁海燕. 2007. 釀酒葡萄品質(zhì)評價及其對氣象條件的響應(yīng)[J]. 生態(tài)學(xué)報，24（2）： 429-436. doi：10. 3321/j.issn：1000-0933.2007.02.039. [Zhang X Y，Kang Y L，Yuan H Y. 2007. The quality evaluation of wine grape and its respond to weather condition[J]. Acta Ecological Sinica，24（2）： 429-436.]

張玉榮，周顯青，彭超. 2021. 不同儲藏年限稻谷的品質(zhì)及蒸谷米加工適應(yīng)性分析[J]. 食品科學(xué)，42（9）：39-45. doi：10.7506/spkx1002-6630-20200422-281. [Zhang Y R，Zhou X Q，Peng C. 2021. Quality and parboiled rice processing suitability of rice with different storage times[J]. Food Science，42（9）：39-45.]

趙志峰. 2018. 花椒中的麻味物質(zhì)[M]. 成都：四川大學(xué)出版社： 1-8. [Zhao Z F. 2018. The numb-taste components of Zanthoxylum[M]. Chengdu： Sichuan University Press：1-8.]

（責(zé)任編輯 麻小燕）