趙靜，楊政，張欣，范靜兒，郭曉成

摘要：【目的】由于石榴凍害的綜合氣候指標(biāo)研究相對較少，這就造成石榴氣象服務(wù)存在一定局限性。且石榴凍害指標(biāo)分級制訂存在差異，缺乏統(tǒng)一標(biāo)準(zhǔn)，不利于石榴凍害監(jiān)測預(yù)報(bào)預(yù)警和災(zāi)情評估。因此，對臨潼石榴凍害綜合氣候指標(biāo)進(jìn)行分析，為石榴凍害的監(jiān)測預(yù)報(bào)預(yù)警服務(wù)和災(zāi)情評估提供技術(shù)支撐?！痉椒ā客ㄟ^實(shí)地走訪、資料查閱等方式調(diào)查石榴凍害災(zāi)情，建立臨潼石榴凍害綜合氣候指標(biāo)分級標(biāo)準(zhǔn)?；?991—2020年37個(gè)凍害實(shí)例，運(yùn)用主成分分析、K-means聚類分析，對臨潼石榴凍害綜合氣候指標(biāo)進(jìn)行研究分析。【結(jié)果】極端最低氣溫、低溫持續(xù)日數(shù)、過程有害積寒與石榴凍害密切相關(guān)，通過對石榴凍害數(shù)據(jù)標(biāo)準(zhǔn)化處理，構(gòu)建石榴凍害綜合氣候等級指標(biāo)（Z），包括：輕度，-7.56≤Z<-6.77;中度，-8.36≤Z<-7.56;重度，Z<-8.36。通過與典型年對比驗(yàn)證，并進(jìn)行智能網(wǎng)格預(yù)報(bào)應(yīng)用，發(fā)現(xiàn)實(shí)際發(fā)生的凍害情況與本文建立的臨潼石榴凍害綜合氣候等級指標(biāo)吻合較好，為氣象服務(wù)石榴高質(zhì)量發(fā)展提供模式支撐。【結(jié)論】通過典型年對比驗(yàn)證及智能網(wǎng)格預(yù)報(bào)應(yīng)用，石榴凍害綜合氣候指標(biāo)能較好的預(yù)報(bào)預(yù)測石榴凍害風(fēng)險(xiǎn)，為氣象服務(wù)石榴高質(zhì)量發(fā)展提供模式支撐。但暫未對臨潼石榴遭受凍害后的植物生理指標(biāo)進(jìn)行測定，其他氣象因子對石榴凍害的影響程度及其指標(biāo)建立和等級劃分等工作，還有待今后進(jìn)一步研究。

關(guān)鍵詞：石榴;凍害;氣候指標(biāo);等級劃分;智能網(wǎng)格預(yù)報(bào)應(yīng)用

文章編號：2096-8108（2022）03-0029-05中圖分類號：S665.4文獻(xiàn)標(biāo)識碼：A

Classification and Application of Comprehensive Climatic Index of Pomegranate Freezing Injury

ZHAO Jing1，YANG Zheng1， ZHANG Xin1，F(xiàn)AN Jinger1，GUO Xiaocheng2

（1.Meteorological Bureau of Lintong County， Xi an 710600， China;

2. Fruit Industry Technology Promotion Center of Xian City， Xian 710600，China）

Abstract：【Objective】The research on comprehensive climatic index of pomegranate freezing injury is relatively few， which results in some limitations of pomegranate meteorological service. In addition， there are differences in classification of pomegranate freezing injury index and lack of unified standard， which is not conducive to monitor， forecast the assessment of pomegranate freezing injury. Therefore， the comprehensive climate index of pomegranate freezing injury in Lintong was analyzed to provide technical support for monitoring， forecasting and the service of pomegranate freezing injury assessment.【Methods】The comprehensive climate index classification standard of pomegranate freezing injury in Lintong was established by field visit and data reference. Based on 37 freezing injury cases from 1991 to 2020， principal component analysis and K-means clustering analysis were used to study and analyze comprehensive climatic indexes of pomegranate freezing damage in Lintong. 【Results】 The extreme minimum temperature， the number of days of low temperature， and the process harmful cold accumulation were closely related to pomegranate freezing injury. Through the standardized processing of pomegranate freezing injury data， the comprehensive climate grade index （Z） of pomegranate freezing injury was constructed， including mild， -7.56≤Z<-6.77; medium， -8.36≤Z<-7.56; sever Z<-8.36. 【Conclusion】Through the verification of many years and the application of intelligent grid prediction， the comprehensive climate index of pomegranate freezing injury could better predict the risk of pomegranate freezing injury， and provided model support for the high-quality development of pomegranate meteorological services. However， the physiological indexes of pomegranate after freezing injury in Lintong had not been determined. The influence of other meteorological factors on the freezing injury of pomegranate and the establishment and classification of indexes needed to be further studied in the future.18C28418-DC83-4079-8820-5389436BE7FA

Keywords：pomegranate; freezing injury; climate indicator; grade division; application of intelligent grid forecast

石榴作為陜西臨潼特色優(yōu)勢產(chǎn)業(yè)，一直以來受到政府有關(guān)部門及社會各界的高度重視和大力支持，臨潼現(xiàn)有石榴種植面積5 333余hm2，在驪山北麓形成了一條東西長約30 km，南北寬約5～10 km的石榴林帶，全區(qū)10個(gè)街辦的2萬余農(nóng)戶從中受益，農(nóng)戶年均石榴產(chǎn)業(yè)收入達(dá)萬元左右，占當(dāng)?shù)毓r(nóng)年收入的約80%。在當(dāng)前全球變暖的背景下，北半球區(qū)域性極端低溫事件頻繁發(fā)生[1-2]，石榴遭受凍害的幾率明顯增加，對石榴生產(chǎn)造成不同程度的損失。曾有研究引起石榴凍害發(fā)生的氣象條件，并提出相應(yīng)的預(yù)防措施[3-10]，但針對石榴凍害的綜合氣候指標(biāo)研究相對較少，這就造成石榴氣象服務(wù)存在一定局限性。且石榴凍害指標(biāo)分級制訂存在差異，缺乏統(tǒng)一標(biāo)準(zhǔn)，不利于石榴凍害監(jiān)測預(yù)報(bào)預(yù)警和災(zāi)情評估[11]。因此，本文在前人研究成果的基礎(chǔ)上，以臨潼區(qū)自然條件下發(fā)生的凍害為基礎(chǔ)，結(jié)合實(shí)地調(diào)查、翻閱資料所得到的低溫過程為研究對象，確定臨潼石榴凍害臨界氣溫，并對凍害實(shí)例進(jìn)行致災(zāi)因子分析，運(yùn)用聚類分析等方法，構(gòu)建臨潼石榴精細(xì)化凍害綜合氣候指標(biāo)，為石榴凍害的監(jiān)測預(yù)報(bào)預(yù)警服務(wù)和災(zāi)情評估提供技術(shù)支撐。

1數(shù)據(jù)來源

1996-2020年臨潼區(qū)石榴產(chǎn)量數(shù)據(jù)，1991-2020年冬季（當(dāng)年12月至次年2月）臨潼國家氣象觀測站的逐日平均氣溫、逐日最低氣溫和過程降溫幅度等數(shù)據(jù)。

2研究方法

2.1計(jì)算過程有害積寒

采用文獻(xiàn)[12-13]提出的過程有害積寒概念，作為臨潼石榴凍害過程中的氣象因子進(jìn)行深入分析。過程有害積寒是指，果樹在低溫過程中逐時(shí)低于臨界受害氣溫的累積值，其計(jì)算公式可表示為：X=14∑Nn=1（Tc-Tmin）2/（Tm-Tmin）（1）式（1）中，X代表過程有害積寒（單位：℃·d），N代表低溫過程持續(xù)日數(shù)（單位：d），Tc代表凍害臨界氣溫（單位：℃），Tmin代表日最低氣溫（單位：℃），Tmin≤Tc，Tm代表日平均氣溫（單位：℃）。

2.2致災(zāi)因子標(biāo)準(zhǔn)化

導(dǎo)致臨潼石榴凍害的致災(zāi)因子較多，為了增強(qiáng)數(shù)據(jù)可比性，消除量綱影響，需要對各致災(zāi)因子進(jìn)行標(biāo)準(zhǔn)化處理[14-17]，公式為：Xi=X′i-∑ni=1（X′i-）2n（2）處理后的標(biāo)準(zhǔn)化致災(zāi)因子最小值為0，最大值為1，其與原數(shù)值相比，只是數(shù)據(jù)形式的變換，將因子處理成量綱為1，保留了因子的統(tǒng)計(jì)學(xué)和生物學(xué)意義[18]，且不影響數(shù)據(jù)研究及分析工作。

2.3確定凍害臨界氣溫

趙靜等[19-21]針對臨潼石榴受凍特點(diǎn)做出研究，并提出-10 ℃是臨潼石榴受凍的臨界溫度，因此，本文確定-10 ℃為臨潼石榴主栽產(chǎn)品凍害的臨界氣溫，即當(dāng)最低氣溫≤-10 ℃時(shí)，石榴凍害過程開始。

3結(jié)果分析

3.1確定凍害致災(zāi)因子

本文主要分析了凍害過程極端最低氣溫、持續(xù)日數(shù)、有害積寒、日最低氣溫降溫幅度等[22-24]以氣溫作為主要分析對象的氣象因子。通過各因子與凍害實(shí)例進(jìn)行相關(guān)分析，結(jié)果（見表1）表明：過程極端最低氣溫（X1）、≤-10 ℃低溫持續(xù)日數(shù)（X2）、日最低氣溫≤-10 ℃的過程有害積寒（X3）、≤-10 ℃過程降溫幅度（X4）等與石榴凍害均為顯著相關(guān)關(guān)系[25]。由于研究選取的4個(gè)致災(zāi)因子是基于臨潼石榴災(zāi)情調(diào)查分析的結(jié)果，因此，各個(gè)致災(zāi)因子具有一定的代表性。

3.2計(jì)算臨潼石榴凍害綜合氣候指標(biāo)

通過計(jì)算各個(gè)致災(zāi)因子間的相關(guān)性發(fā)現(xiàn)，4個(gè)致災(zāi)因子中除了X4以外的其余致災(zāi)因子，兩兩間均呈顯著相關(guān)，且達(dá)到0.05顯著性水平（見表1），這表明X1、X2、X3之間互有影響，并不獨(dú)立，適合運(yùn)用主成分分析法進(jìn)行數(shù)據(jù)分析。

因此，采用主成分分析方法時(shí)，對X1、X2、X3這3個(gè)致災(zāi)因子進(jìn)行綜合簡化并構(gòu)建公式。將各致災(zāi)因子進(jìn)行標(biāo)準(zhǔn)化處理，運(yùn)用R統(tǒng)計(jì)軟件中的主成分分析函數(shù)，對標(biāo)準(zhǔn)化后的石榴凍害3個(gè)致災(zāi)因子X1、X2、X3序列進(jìn)行分析（見表2），提取第1個(gè)主成分，其累計(jì)方差貢獻(xiàn)率已經(jīng)達(dá)到84.6%，說明第1個(gè)主成分能夠充分代表數(shù)據(jù)間的波動原因。

通過求解正交矩陣Q（見表3），可構(gòu)建臨潼石榴凍害綜合氣候指標(biāo)Z，表示為：

Z=0.585X1-0.560X2-0.586X3（3）

因?yàn)橹鞒煞直磉_(dá)式中系數(shù)的正負(fù)本身并沒有實(shí)際意義，有意義的是不同主成分之間統(tǒng)一變量的系數(shù)大小：由式（3）可知，在第1個(gè)主成分中，3個(gè)變量的系數(shù)比較接近，這意味著第1個(gè)主成分是對3個(gè)變量的綜合測度，即體現(xiàn)3個(gè)致災(zāi)因子的信息。

凍害綜合氣候指標(biāo)Z與冬季極端最低氣溫X1呈正相關(guān)關(guān)系，而與X2、X3呈負(fù)相關(guān)關(guān)系，即表明冬季極端最低氣溫越低，≤-10 ℃低溫持續(xù)日數(shù)越多、≤-10 ℃有害積寒越大，則凍害綜合氣候指標(biāo)Z越小，這與臨潼石榴實(shí)際發(fā)生的凍害過程一致。

3.3凍害綜合氣候指標(biāo)分級

K-means聚類分析是經(jīng)典的聚類算法，常應(yīng)用在氣象數(shù)據(jù)分析和農(nóng)業(yè)氣象災(zāi)害等級劃分等方面[26-29]。本文采用R統(tǒng)計(jì)軟件中的K-means聚類分析函數(shù)，通過分析37個(gè)凍害等級對應(yīng)的凍害綜合氣候指標(biāo)，并獲取各凍害等級對應(yīng)的Z值，研究中選擇凍害綜合氣候指標(biāo)作為聚類分析的變量，1996—2020年臨潼區(qū)石榴產(chǎn)量數(shù)據(jù)作為標(biāo)識變量，聚類數(shù)分為3類，初始類中心按照系統(tǒng)指定方式確定。結(jié)果表明：3個(gè)聚類中心點(diǎn)經(jīng)過2次迭代后，得到最終的各聚類中心分別為-6.77，-7.56，-8.36（見圖1）。由此可知，臨潼石榴的輕度、中度、重度凍害綜合氣候指標(biāo)分級分別為-7.56≤Z<-6.77，-8.36≤Z<-7.56，Z<-8.36，在業(yè)務(wù)應(yīng)用中設(shè)置相應(yīng)風(fēng)險(xiǎn)等級顏色，輕度用藍(lán)色表示、中度用黃色表示、重度用紅色表示（見表4）。18C28418-DC83-4079-8820-5389436BE7FA

3.4典型年對比驗(yàn)證及網(wǎng)格預(yù)報(bào)應(yīng)用

3.4.1典型年份對比驗(yàn)證

通過結(jié)合臨潼石榴畝（667 m2）產(chǎn)數(shù)據(jù)與凍害綜合氣候指標(biāo)對應(yīng)關(guān)系發(fā)現(xiàn)，凍害綜合氣候指標(biāo)及其分級與實(shí)際發(fā)生的凍害情況較吻合，即可用凍害綜合氣候指標(biāo)分析臨潼石榴遭受凍害的程度，并對畝（667 m2）產(chǎn)記錄未知年份和歷史災(zāi)情進(jìn)行驗(yàn)證。

1）計(jì)算1991-1993年凍害綜合氣候指標(biāo)Z發(fā)現(xiàn)，1991年Z值為-14.71，達(dá)重度凍害;1992年Z值為-9.50，達(dá)重度凍害;1993年Z值為-8.78，達(dá)重度凍害。根據(jù)查閱資料和走訪調(diào)查，1991年臨潼石榴遭受了歷史上毀滅性災(zāi)害，驗(yàn)證了臨潼石榴在相應(yīng)年份遭受的凍害等級。

2）2020年凍害綜合氣候指標(biāo)Z值為-8.44，達(dá)重度凍害。根據(jù)園區(qū)調(diào)查，由于2020年冬季石榴樹體受凍，造成石榴果量減少，減產(chǎn)明顯，通過z值驗(yàn)證。

3.4.2精細(xì)化網(wǎng)格預(yù)報(bào)應(yīng)用

為更好將指標(biāo)及分級應(yīng)用于氣象服務(wù)，選取2020年12月27-31日降溫過程為對象模擬進(jìn)行石榴精細(xì)化網(wǎng)格預(yù)報(bào)服務(wù)。

應(yīng)用個(gè)例：臨潼區(qū)氣象臺發(fā)布重要天氣報(bào)告，受高空偏南暖濕氣流和北方冷空氣南下影響，預(yù)計(jì)12月27-31日將有一次降雪、吹風(fēng)、強(qiáng)降溫天氣過程。27日晚-28日城區(qū)有小雪，南部山區(qū)小到中雪，28-30日有5級左右偏北風(fēng)，陣風(fēng)可達(dá)6級，30-31日最低氣溫北部平原-10～-8 ℃，南部山區(qū)-12～-10 ℃。

精細(xì)化網(wǎng)格預(yù)報(bào)30日、31日2：00后均出現(xiàn)可能影響石榴生長的凍害溫度，根據(jù)公式（3）計(jì)算Z值為：城區(qū)-6.97，為輕度風(fēng)險(xiǎn);山區(qū)-8.73，為重度風(fēng)險(xiǎn)（見表5）。

根據(jù)表5，基于臨潼GIS地圖制作臨潼區(qū)石榴產(chǎn)區(qū)受凍風(fēng)險(xiǎn)圖用于氣象服務(wù)，為石榴生產(chǎn)管理提供保障服務(wù)（見圖2）。

4討論與結(jié)論

1）影響石榴凍害的主要?dú)庀笠蜃訛闃O端最低氣溫、低溫持續(xù)日數(shù)和低溫過程有害積寒。將3個(gè)致災(zāi)因子進(jìn)行標(biāo)準(zhǔn)化處理后，得到的石榴凍害綜合氣候指標(biāo)，包括：輕度，-7.56≤Z<-6.77;中度，-8.36≤Z<-7.56;重度，Z<-8.36。

2）通過典型年對比驗(yàn)證及智能網(wǎng)格預(yù)報(bào)應(yīng)用，石榴凍害綜合氣候指標(biāo)能較好的預(yù)報(bào)預(yù)測石榴凍害風(fēng)險(xiǎn)，為氣象服務(wù)石榴高質(zhì)量發(fā)展提供模式支撐。

3）由于引起果樹凍害的因子復(fù)雜，且凍害程度還和植物學(xué)因子的有關(guān)[30-32]。本文目前只針對影響石榴凍害重要的氣象因子—?dú)鉁剡M(jìn)行了分析，暫未對臨潼石榴遭受凍害后的植物生理指標(biāo)進(jìn)行測定。其他氣象因子對石榴凍害的影響程度及其指標(biāo)建立和等級劃分等工作，還有待今后進(jìn)一步研究。

4）隨著國內(nèi)外石榴品種軟籽化發(fā)展趨勢，此模式將根據(jù)臨潼主栽品種的變化和服務(wù)需要及時(shí)修訂和完善。

參考文獻(xiàn)

[1]龔志強(qiáng)，王曉娟，崔冬林，等.區(qū)域性極端低溫事件的識別及其變化特征[J].應(yīng)用氣象學(xué)報(bào)，2012，23（2）：195-204.

[2]李清泉，孫丞虎，袁媛，等.近20年我國氣候監(jiān)測診斷業(yè)務(wù)技術(shù)的主要進(jìn)展[J].應(yīng)用氣象學(xué)報(bào)，2013，24（6）：666-676.

[3]馮玉增，李戰(zhàn)鴻，趙艷麗，等.石榴凍害與低溫程度的關(guān)系[J].河北果樹，2003（1）：14-17.

[4]謝建華.石榴樹凍害類型及影響因素[J].陜西農(nóng)業(yè)科學(xué)，2012（6）：267.

[5]李劍雄.北京房山地區(qū)石榴凍害調(diào)查分析[J].現(xiàn)代園藝，2018（8）：7-8.

[6]劉昭武，張其忠，蔡永偉.濱州市石榴生長的氣象條件分析[J].時(shí)代農(nóng)機(jī)，2015（1）：124-125.

[7]李曦.豫西地區(qū)軟籽石榴凍害成因及防凍措施[J].現(xiàn)代農(nóng)業(yè)科技，2018（14）：94-97.

[8]劉文獻(xiàn)，趙丹，王坤宇等.石榴凍害的發(fā)生與預(yù)防[J].現(xiàn)代農(nóng)村科技，2010（3）：21.

[9]畢潤霞，馬敏，羅華，等.棗莊市石榴凍害預(yù)防及凍后補(bǔ)救措施[J].現(xiàn)代農(nóng)業(yè)科技，2018（17）：90-92.

[10]劉英勝.石榴凍害的原因及預(yù)防措施[J].河北果樹，2017（2）：31-32.

[11]楊凱，陳彬彬，陳惠，等.福建省臺灣青棗寒害綜合氣候指標(biāo)與等級劃分[J].應(yīng)用氣象學(xué)報(bào)，2020，31（4）：427-434.

[12]杜堯東，李春梅，毛慧琴.廣東省香蕉與荔枝寒害致災(zāi)因子和綜合氣候指標(biāo)研究[J].生態(tài)學(xué)雜志，2006，25（2）：225-230.

[13]柏秦鳳.華南寒害致災(zāi)氣候因子及綜合指數(shù)研究[M].北京：中國氣象科學(xué)研究院，2008：4-5.

[14]閔晶晶，孫景榮，劉還珠，等.一種改進(jìn)的BP算法及在降水預(yù)報(bào)中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào)，2010，21（1）：55-62.

[15]陳家金，李麗純，林晶，等.福建省枇杷氣象災(zāi)害綜合風(fēng)險(xiǎn)評估[J].應(yīng)用氣象學(xué)報(bào)，2014，25（2）：232-241.

[16]何斌.陜西省農(nóng)業(yè)干旱風(fēng)險(xiǎn)評估方法研究[D].西安：西安理工大學(xué)，2017：20-21.

[17]王永利，侯瓊，苗百嶺，等.內(nèi)蒙古馬鈴薯干旱風(fēng)險(xiǎn)區(qū)劃[J].應(yīng)用氣象學(xué)報(bào)，2017，28（4）：504-512.

[18]李登科，權(quán)文婷，謝飛舟.陜西省小麥條銹病氣象預(yù)報(bào)模型[J].干旱氣象，2017，35（1）：128-133.

[19]趙靜，王昱，何倩茹.石榴生長主要物候期氣象服務(wù)指標(biāo)分析[J].資源與環(huán)境科學(xué)，2015（19）：263-264.

[20]李艷麗，趙靜，劉映寧，等.低溫雨雪冰凍災(zāi)害對我省石榴影響分析[J].陜西農(nóng)業(yè)科學(xué)，2009（1）：116-118.18C28418-DC83-4079-8820-5389436BE7FA

[21]趙靜，賈毅萍，柏永耀.石榴受凍年份的氣候特點(diǎn)[J].山西果樹，2007（6）：24-25.

[22]楊凱，陳彬彬，陳惠，等.基于寒害過程的福建芒果種植氣候風(fēng)險(xiǎn)區(qū)劃[J].中國農(nóng)業(yè)氣象，2019，40（11）：723-732.

[23]楊凱，林晶，陳惠，等.福建枇杷低溫害臨界溫度和綜合氣候指標(biāo)[J].中國農(nóng)業(yè)氣象，2013，34（3）：468-473.

[24]陳家金，王加義，黃川容，等.福建省引種臺灣青棗的寒凍害風(fēng)險(xiǎn)分析與區(qū)劃[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2013，21（12）：1537-1544.

[25]劉靜，馬力文，張曉煜，等.春小麥干熱風(fēng)災(zāi)害監(jiān)測指標(biāo)與損失評估模型方法探討———以寧夏引黃灌區(qū)為例[J].應(yīng)用氣象學(xué)報(bào)，2004，15（2）：217-225.

[26]張洪兵，賈來喜，李璐.R寶典[M].北京：電子工業(yè)出版社，2007：353-364.

[27]任義方，趙艷霞，王春乙.河南省冬小麥干旱保險(xiǎn)風(fēng)險(xiǎn)評估與區(qū)劃[J].應(yīng)用氣象學(xué)報(bào)，2011，22（5）：537-548.

[28]田俊，崔海建.江西省雙季早稻灌漿乳熟期高溫?zé)岷τ绊懺u估[J].中國農(nóng)業(yè)氣象，2015，36（1）：67-73.

[29]劉瑞娜，楊太明，陳鵬，等.安徽省油菜花期連陰雨災(zāi)害損失評估指標(biāo)[J].中國農(nóng)業(yè)氣象，2016，37（4）：471-478.

[30]陳由強(qiáng)，葉冰瑩，高一平，等.低溫脅迫下枇杷幼葉細(xì)胞內(nèi)Ca2+水平及細(xì)胞超微結(jié)構(gòu)變化的研究[J].武漢植物學(xué)研究，2000，18（2）：138-142.

[31]馮美利.杧果葉片組織的細(xì)胞結(jié)構(gòu)與性的初步研究[J].熱帶農(nóng)業(yè)科技，2003，26（1）：5-7.

[32]胡春霞，王麗，湯杰.低溫對南果梨的生理生化指標(biāo)的影響[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào)，2009，40（3）：349-352.

收稿日期：2022-03-10

基金項(xiàng)目：西安市科技局“科技+”行動計(jì)劃農(nóng)業(yè)技術(shù)研發(fā)項(xiàng)目（項(xiàng)目編號：XA2020-NYJSYF-0001）

第一作者簡介：趙靜（1974-），女，陜西白水人，大學(xué)本科學(xué)歷，副研級高工，主要從事農(nóng)業(yè)氣象服務(wù)研究。電話：13609180183;E-mail：359289309@qq.com18C28418-DC83-4079-8820-5389436BE7FA