李 楠，薛曉萍

(山東省氣候中心，山東 濟南 250031)

引 言

設(shè)施農(nóng)業(yè)小氣候條件既與外界大氣候條件不同，又與其密切相關(guān)[1-3]。為清楚地了解設(shè)施內(nèi)小氣候變化規(guī)律，近幾年山東、天津、河北等氣象部門先后在具有當(dāng)?shù)卮硇缘娜展鉁厥覂?nèi)安裝了自動觀測設(shè)備，其中，山東組建了標(biāo)準(zhǔn)化的設(shè)施農(nóng)業(yè)小氣候觀測站網(wǎng)，為開展日光溫室氣象服務(wù)提供了觀測資料基礎(chǔ)保障。山東在設(shè)施農(nóng)業(yè)小氣候預(yù)報以及基于小氣候預(yù)報結(jié)果的低溫災(zāi)害影響預(yù)報方面已形成了較為成熟的技術(shù)標(biāo)準(zhǔn)和業(yè)務(wù)體系，但在小氣候監(jiān)測的定量評價方面還有所欠缺。如何定量評價日光溫室內(nèi)光、溫、濕等小氣候要素對作物生長的適宜程度，提出更有針對性的管理措施，改善溫室小氣候環(huán)境，促進蔬菜生長發(fā)育，是設(shè)施農(nóng)業(yè)氣象服務(wù)關(guān)注的重要內(nèi)容之一[4-6]。

本文擬在三要素小氣候適宜度模型基礎(chǔ)上，增加作物根系活動層且變化相對較明顯的10 cm地溫要素，構(gòu)建基于日光溫室內(nèi)氣溫、空氣相對濕度、太陽輻射和10 cm地溫的四要素小氣候適宜度定量評價模型，豐富模型的物理意義，以期進一步提高設(shè)施農(nóng)業(yè)氣象服務(wù)的技術(shù)水平以及為決策服務(wù)提供依據(jù)。

1 資料和方法

1.1 資 料

所用資料為山東壽光、寒亭設(shè)施農(nóng)業(yè)小氣候自動觀測站日光溫室內(nèi)小氣候要素的逐小時、逐日資料，觀測站安裝在當(dāng)?shù)卦O(shè)施農(nóng)業(yè)實際生產(chǎn)中具有代表性的日光溫室內(nèi)，坐北朝南，東西長80 m，南北寬12 m，脊高4.5 m，墻體為土墻，墻體厚1.5 m。種植作物為番茄，10月開始移栽，12月開始進入采收期。

日光溫室內(nèi)小氣候自動觀測站按照國家標(biāo)準(zhǔn)《設(shè)施農(nóng)業(yè)小氣候觀測規(guī)范 日光溫室和塑料大棚》(GB/T 38757—2020)[21]的要求安裝，山東地區(qū)日光溫室主要生產(chǎn)服務(wù)季是當(dāng)年11月至次年4月，因此日光溫室小氣候數(shù)據(jù)以生產(chǎn)服務(wù)季為周期。

在適宜度定量評價模型的檢驗過程中所用的番茄生長發(fā)育農(nóng)業(yè)氣象指標(biāo)由山東省農(nóng)業(yè)科學(xué)院蔬菜花卉研究所提供；番茄生長狀況和產(chǎn)量資料來自小氣候觀測站所在日光溫室內(nèi)的番茄實際觀測值。

1.2 方 法

1.2.1 溫室內(nèi)氣溫適宜度評價模型

日光溫室內(nèi)逐小時氣溫適宜度定量評價模型公式[2]如下：

(1)

式中：T(ti)為日光溫室內(nèi)第i時刻的氣溫適宜度；ti為設(shè)施內(nèi)第i時刻的氣溫；tl、th和to(℃)分別為作物生長發(fā)育所需要的最低氣溫、最高氣溫和最適宜氣溫；B=(th-to)/(to-tl)。

日光溫室逐日氣溫適宜度定量評價公式[2]如下：

(2)

式中：Td為日光溫室內(nèi)日氣溫的適宜度；T(ti)為日光溫室內(nèi)第i時刻的氣溫適宜度指數(shù)；n為一日內(nèi)日光溫室內(nèi)逐時氣溫有觀測記錄的次數(shù)。

1.2.2 溫室內(nèi)空氣相對濕度評價模型

日光溫室內(nèi)逐小時空氣相對濕度的適宜度定量評價模型公式[2]如下：

(3)

式中：U(ui)為日光溫室內(nèi)第i時刻的空氣相對濕度適宜度；ui為日光溫室內(nèi)第i時刻的空氣相對濕度；ul、uh和uo(%)分別為作物生長發(fā)育所需要的最小、最大和最適宜空氣相對濕度；D=(uh-uo)/(uo-ul)。

日光溫室逐日空氣相對濕度適宜度定量評價公式[2]如下：

(4)

式中：Ud為逐日日光溫室內(nèi)空氣相對濕度適宜度指數(shù)；U(ui)為日光溫室內(nèi)第i時刻的空氣相對濕度適宜度指數(shù)；m為日光溫室內(nèi)08:00—20:00逐時空氣相對濕度有觀測記錄的次數(shù)。

1.2.3 溫室內(nèi)太陽輻射適宜度評價模型

因為日光溫室內(nèi)的最大太陽輻射輻照度直接影響溫室內(nèi)的最低氣溫和最高氣溫，進而影響日光溫室內(nèi)作物生長，因此，在日光溫室黃瓜適宜度評價模型研究成果[2]基礎(chǔ)上，將原有太陽輻射適宜度評價模型的最適宜及最不適宜兩個等級細(xì)化，改寫模型公式如下：

(5)

式中：Q(qi)為日光溫室內(nèi)第i時刻的太陽輻射適宜度；qi為日光溫室內(nèi)第i時刻的太陽總輻射輻照度；ql、qh和qo(W·m-2)分別為可達到溫室作物所處發(fā)育階段生長所需最低、最高和最適宜氣溫的日最大總輻射輻照度。

日光溫室逐日太陽輻射適宜度指數(shù)評價公式如下：

(6)

式中：Qd為逐日的日光溫室內(nèi)太陽輻射適宜度指數(shù)；Q(qj)為日光溫室內(nèi)揭簾至蓋簾時間內(nèi)第j時刻的太陽輻射適宜度；p為日光溫室內(nèi)揭簾到蓋簾時間內(nèi)逐時太陽輻射有觀測記錄次數(shù)。

1.2.4 溫室內(nèi)10 cm地溫適宜度評價模型

通過對2017年11月至2018年4月寒亭站和壽光站日光溫室內(nèi)逐日10 cm地溫與氣溫進行相關(guān)性分析，兩者顯著相關(guān)，擬合方程決定系數(shù)R2為0.82(圖1)。因此，與溫室內(nèi)氣溫適宜度評價模型相似，構(gòu)建日光溫室內(nèi)10 cm地溫適宜度評價模型如下：

(7)

式中：T(dti)為日光溫室內(nèi)第i時刻10 cm地溫適宜度；dti、dtl、dth和dto(℃)分別為日光溫室內(nèi)第i時刻10 cm地溫、作物生長發(fā)育所需要的最低10 cm地溫、最高10 cm地溫和最適宜10 cm地溫；P=(dth-dto)/(dto-dtl)。T(dti)是在0～1變化的拋物線函數(shù)(圖2)，反映10 cm地溫條件從不適宜到適宜再到不適宜的連續(xù)變化過程。逐日10 cm地溫適宜度定量評價公式如下：

(8)

式中：Dd為日光溫室內(nèi)10 cm地溫逐日適宜度指數(shù)；T(dti)為第i時刻的10 cm地溫適宜度；n為一日內(nèi)日光溫室內(nèi)逐時10 cm地溫有觀測記錄次數(shù)。

圖1 2017年11月至2018年4月寒亭和壽光站日光溫室內(nèi)10 cm地溫與氣溫散點圖Fig.1 The scatter plot between 10 cm ground temperature and air temperature in heliogreenhouses at Hanting and Shouguang station from November 2017 to April 2018

圖2 日光溫室番茄小氣候適宜度指數(shù)對10 cm地溫的反應(yīng)過程Fig.2 Reaction process of microclimate suitability of tomato in heliogreenhouse to 10 cm ground temperature

1.2.5 溫室內(nèi)綜合小氣候要素適宜度評價模型

將魏瑞江等[2]基于溫度、空氣相對濕度和輻射輻照度3個因子的等權(quán)重小氣候要素適宜度模型改進為包含10 cm地溫在內(nèi)的4個因子的等權(quán)重小氣候要素適宜度模型，公式如下：

(9)

式中：Sd為日光溫室內(nèi)小氣候要素綜合日適宜度指數(shù)；Td、Ud、Dd和Qd分別為溫室內(nèi)氣溫、空氣相對濕度、10 cm地溫及太陽輻射逐日適宜度指數(shù)。

1.2.6 最優(yōu)分割法

采用最優(yōu)分割法[22]劃分小氣候要素適宜度評價等級，最優(yōu)分割法是在對n個樣品分段并使組內(nèi)離差平方和最小的分割方法。

2 結(jié)果分析

利用2017年和2018年日光溫室生產(chǎn)服務(wù)季內(nèi)小氣候要素逐時觀測資料及番茄產(chǎn)量資料進行分析檢驗，兩年產(chǎn)量分別為135 686.6和149 253.7 kg·hm-2，在品種特性、肥水管理、土壤質(zhì)地和栽培管理水平一致條件下，產(chǎn)量取決于設(shè)施內(nèi)的小氣候條件。為了對比小氣候適宜度定量評價模型改進前與改進后的效果，應(yīng)用魏瑞江等[2]的三要素適宜度模型與本文提出的四要素適宜度模型分別計算適宜度指數(shù)。產(chǎn)量采用10 d平均值，綜合小氣候適宜度指數(shù)亦用與之對應(yīng)的10 d平均值。

圖3為2017、2018年日光溫室生產(chǎn)服務(wù)季每10 d平均番茄產(chǎn)量分別與三要素、四要素綜合小氣候適宜度指數(shù)的散點圖?？梢钥闯觯?017年和2018年每10 d番茄產(chǎn)量與三要素綜合小氣候適宜度指數(shù)的決定系數(shù)分別為0.59和0.69，均通過0.01的顯著性檢驗；每10 d番茄產(chǎn)量隨著四要素綜合小氣候適宜度指數(shù)的增加而增加，2017年和2018年兩者的決定系數(shù)R2分別為0.93和0.90，均通過0.01的顯著性檢驗，且決定系數(shù)明顯優(yōu)于三要素綜合模型。

日光溫室內(nèi)綜合小氣候適宜度指數(shù)越高，溫室內(nèi)小氣候條件越適宜于番茄生長發(fā)育，其產(chǎn)量越高。綜上所述，改進后的四要素適宜度模型評價效果優(yōu)于三要素適宜度模型。

圖3 2017年(a、c)、2018年(b、d)日光溫室生產(chǎn)服務(wù)季每10 d平均番茄產(chǎn)量與三要素(a、b)及四要素(c、d)綜合小氣候適宜度指數(shù)散點圖Fig.3 The scatter plots between average tomato yield and integrated microclimate suitability index of three (a, b) and four (c, d) elements for every 10 days in heliogreenhouses during production seasons in 2017 (a, c) and 2018 (b, d)

3 模型的應(yīng)用

3.1 適宜度評價等級指標(biāo)

通過計算2012—2017年生產(chǎn)服務(wù)季寒亭日光溫室內(nèi)番茄的氣溫、空氣相對濕度、太陽輻射、10 cm地溫及小氣候綜合適宜度指數(shù)，發(fā)現(xiàn)日光溫室內(nèi)逐日氣溫適宜度指數(shù)最小值為0.0，最大值為0.85；日空氣相對濕度適宜度指數(shù)最小值為0.0，最大值為0.80；日太陽輻射適宜度指數(shù)最小值為0，最大值為1.0；日10 cm地溫適宜度指數(shù)最小值為0，最大值為1.0；日綜合適宜度指數(shù)最小值為0.0，最大值為0.78。利用最優(yōu)分割法，結(jié)合作物實際生長狀況，將各要素適宜度劃分為3個等級，如表1所示。

表1 適宜度評價等級指標(biāo)Tab.1 Grade indexes of suitability evaluation

3.2 小氣候適宜度評價

3.2.1 各要素適宜度

表2列出2018年11月至2019年4月寒亭日光溫室番茄小氣候要素不同等級適宜度各月出現(xiàn)日數(shù)占比?？梢钥闯觯瑲鉁剡m宜度為“適宜”等級的日數(shù)3月占比最高，為71%，11月、12月和1月占比最低；“不適宜”等級的日數(shù)1月占比最高，為55%，其他月份均在10%以下；生產(chǎn)服務(wù)季內(nèi)氣溫要素較適宜于番茄生長。空氣相對濕度適宜度未出現(xiàn)“適宜”等級；“不適宜”等級的日數(shù)11月至次年1月占比最高，為80%以上，2月開始適宜度明顯提高，大部分日數(shù)處于“較適宜”等級。日光溫室太陽輻射適宜度為“適宜”的等級4月和3月占比最高，分別為50%和35%；“不適宜”等級的日數(shù)1月占比最高為52%，其次11月為37%，可見各地日光溫室內(nèi)接受到的太陽輻射對番茄生長發(fā)育最適宜的月份是4月，其次為3月，整個冬季光照條件均不適宜番茄生長。10 cm地溫適宜度為“適宜”等級的日數(shù)3、4月占比最高，分別為97%和100%，表明進入春季后地溫回升明顯，完全可以提供番茄生長所需熱量，冬季的12月和1月10 cm地溫適宜度為“不適宜”等級的日數(shù)占比達100%，其次為11月；2月適宜度各等級日數(shù)占比差別不大，其中“適宜”等級占比為57%，“較適宜”占比為14%；整體看2月地溫條件適宜。

表2 2018年11月至2019年4月寒亭日光溫室番茄小氣候要素不同等級適宜度各月出現(xiàn)日數(shù)占比Tab.2 The percentages of appearing days of different elements micro-cliamte suitability of tomato with different grade in each month in Hanting heliogreenhouse from November 2018 to April 2019 單位：%

3.2.2 各要素適宜度變化趨勢

圖4為2018年11月至2019年4月寒亭日光溫室番茄小氣候要素旬平均適宜度指數(shù)變化?？梢钥闯?，日光溫室內(nèi)氣溫、10 cm地溫和太陽輻射適宜度變化均呈不規(guī)則“V”字形，與魏瑞江等[2]分析結(jié)果相似。11月上旬至1月下旬呈下降趨勢，之后開始在波動中上升，各要素適宜度指數(shù)的最低值均出現(xiàn)在1月下旬， 2019年1月下旬外界最低氣溫為-12.3 ℃，且出現(xiàn)3 d以上連續(xù)寡照天氣，日光溫室內(nèi)光照不足，導(dǎo)致溫室氣溫、地溫均降低，其中10 cm地溫降至10.8 ℃，接近番茄生長的地溫下限指標(biāo)；空氣相對濕度前期變化趨于平穩(wěn)，適宜度指數(shù)穩(wěn)定在0.2左右，2月開始由于通風(fēng)時間延長，空氣相對濕度下降，適宜度指數(shù)明顯升高，進入3月大風(fēng)天氣增多，通風(fēng)時間間隔縮減，導(dǎo)致適宜度指數(shù)略有下降。

3.2.3 綜合適宜度

圖5為2018年11月至2019年4月寒亭日光溫室番茄不同等級綜合小氣候適宜度各月出現(xiàn)日數(shù)占比?？梢钥闯?，綜合小氣候適宜度為“適宜”等級的3、4月占比最高，分別為87%和77%，其次為12月，占比為68%，1月最低，占比僅為19%；“不適宜”等級1月占比最高，為61%，4月最低，僅為10%，其次為3月的13%。綜上所述，寒亭日光溫室小氣候條件最適宜番茄生長月份為3月，其次為4月和12月，1月最不適宜。

圖4 2018年11月至2019年4月寒亭日光溫室番茄小氣候要素旬平均適宜度指數(shù)變化Fig.4 The variation of ten-day average suitability index of different elements of tomato in heliogreenhouse in Hanting from November 2018 to April 2019

圖5 2018年11月至2019年4月寒亭日光溫室番茄不同等級綜合小氣候適宜度各月出現(xiàn)日數(shù)占比Fig.5 The percentages of appearing days of comprehensive micro-cliamte suitability of tomato with different grades in Hanting heliogreenhouse in each month from November 2018 to April 2019

4 結(jié)論與討論

(1)引入10 cm地溫，構(gòu)建了日光溫室內(nèi)10 cm地溫適宜度評價模型且改進了太陽輻射適宜度評價模型的適宜度等級指標(biāo)。利用等權(quán)重方法構(gòu)建了四要素綜合小氣候適宜度定量評價模型。寒亭日光溫室內(nèi)番茄實際產(chǎn)量與四要素綜合適宜度指數(shù)顯著相關(guān)，且相關(guān)性優(yōu)于三要素適宜度模型。

(2)在計算2012—2016年寒亭站日光溫室內(nèi)單要素及綜合小氣候適宜度指數(shù)的基礎(chǔ)上，利用最優(yōu)分割法，確定了單要素及綜合適宜度指數(shù)等級評價指標(biāo)。

(3)2018年11月至2019年4月寒亭站日光溫室各要素適宜度指數(shù)的最低值均出現(xiàn)在1月下旬，主要是受到災(zāi)害性天氣的影響；最高值出現(xiàn)在3月。

本研究尚未考慮其他小氣候要素對作物生長發(fā)育的影響，這也是之后需要進一步驗證的內(nèi)容。不同作物、不同品種差異很大，在開展設(shè)施農(nóng)業(yè)小氣候監(jiān)測診斷業(yè)務(wù)工作中，對不同地區(qū)的適宜度評價等級指標(biāo)進行本地化應(yīng)用過程中，可參照本文的指標(biāo)劃分方法。