邱璟怡,郭建茂*,周文艷,朱敬宜
(1.南京信息工程大學(xué) 無(wú)錫研究院,江蘇 南京 210044;2.國(guó)家氣候中心,北京 100081;3.南京信息工程大學(xué)江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室 應(yīng)用氣象學(xué)院,江蘇 南京 210044)
陸面過(guò)程也稱(chēng)為陸-氣相互作用,是模擬地表?xiàng)l件變化和研究陸-氣相互作用的重要工具。陸氣相互作用特征對(duì)不同的下墊面特性十分敏感[1],不同地表狀況條件下,陸氣之間物質(zhì)和能量交換情況存在差異,這在很大程度上決定了局地氣候差異,對(duì)區(qū)域和全球氣候有重要影響[2]。針對(duì)地表覆蓋數(shù)據(jù)對(duì)陸面溫度模擬影響的一些研究表明,植被PFT在陸地表面的覆蓋范圍在所有覆蓋類(lèi)型中最大,因而其影響最為顯著[3]。大部分學(xué)者認(rèn)為,植物群落的降溫效果與葉面積指數(shù)呈正相關(guān)關(guān)系[4-5],不同的植被葉面積指數(shù)規(guī)律不同,也造成了對(duì)陸面溫度不同程度的影響,如闊葉林的降溫范圍和降溫幅度比草地要大[6];稻田冠層高度50 cm處溫度(T50)與氣溫(T)的關(guān)系近似有T50=0.767T1.076[7];棉花、玉米、大豆旺盛生長(zhǎng)期覆膜可以使葉面積指數(shù)增大,降溫幅度也分別增大1.1 ℃、0.9 ℃和0.7 ℃[8]。因此,對(duì)植被生長(zhǎng)進(jìn)行更好的參數(shù)化描述、在陸面過(guò)程模式中包含人類(lèi)活動(dòng)和生態(tài)系統(tǒng)過(guò)程并實(shí)現(xiàn)高分辨率模擬是新一代陸面過(guò)程模式發(fā)展的戰(zhàn)略方向[2]。
本研究在春玉米一年一熟主要種植制度東北平原、以冬小麥-夏玉米一年兩熟輪作為主要種植制度的華北平原,在陸面模式中利用具體的作物代替“抽象作物”,引入發(fā)育期的概念,并在不同發(fā)育階段配置更符合實(shí)際的分配系數(shù),使其在理論上更符合實(shí)際,從而改進(jìn)陸面模式的模擬效果。
陸面模式發(fā)展至今經(jīng)歷了4個(gè)階段:第一代“水桶模式”不考慮土壤和植被種類(lèi)[9];第二代“土壤大氣植被傳輸模型”[9]發(fā)現(xiàn),植物蒸散的重要作用,開(kāi)始引入植被方案;第三代“大氣植被相互作用模型”,考慮了生化過(guò)程,使得植物能響應(yīng)氣候變化而動(dòng)態(tài)生長(zhǎng)[10];第四代陸面過(guò)程模式正在發(fā)展中,引入大量的子模式,旨在對(duì)過(guò)程加入更精細(xì)化的描述,使模擬過(guò)程更加完備。
國(guó)家氣候中心大氣-植被相互作用模式(BCC_AVIM),是我國(guó)研發(fā)的國(guó)家氣候中心氣候系統(tǒng)模式(BCC_CSM)的陸面過(guò)程分量模式。BCC_ AVIM1.0在美國(guó)NCAR-CLM3的基礎(chǔ)上,融合了動(dòng)態(tài)植被模式(AVIM2)的生物化學(xué)模塊[11-12],它能對(duì)植被光合作用固定CO2、植被生長(zhǎng)及凋落、土壤呼吸釋放CO2返回大氣的陸地碳循環(huán)等過(guò)程進(jìn)行合理刻畫(huà)[13],已具備對(duì)動(dòng)態(tài)植被和陸面碳循環(huán)過(guò)程的模擬能力。
BCC_AVIM2.0在原有基礎(chǔ)上進(jìn)一步作出了重要改進(jìn)。除了在湖泊、凍土、積雪等方面的改進(jìn)外,在作物生長(zhǎng)方面,實(shí)現(xiàn)了利用動(dòng)態(tài)的預(yù)測(cè)植被生長(zhǎng)的物候方案描述植被動(dòng)態(tài)生長(zhǎng)過(guò)程,更新描述通過(guò)植被冠層的太陽(yáng)輻射傳遞方案,并在植被模塊中加入了稻田參數(shù)化方案[14]。
BCC_AVIM2.0融合動(dòng)態(tài)植被模式、更新動(dòng)態(tài)預(yù)測(cè)模型,對(duì)生長(zhǎng)環(huán)境相對(duì)特殊、水表環(huán)境顯著的稻田采用了新參數(shù)化方案,其結(jié)果都表明了改進(jìn)植被的重要性[13]。在重要的農(nóng)業(yè)發(fā)達(dá)地區(qū),擬進(jìn)一步地將下墊面植被覆蓋類(lèi)型進(jìn)行細(xì)分,并且依據(jù)不同作物的特點(diǎn)進(jìn)行更符合實(shí)際的動(dòng)態(tài)模擬,以提升陸面模式的模擬效果。
本文利用BCC_AVIM2.0進(jìn)行陸面過(guò)程模擬改進(jìn)。在BCC_AVIM2.0的支持文件中,除水稻外,其他作物都采用“crop1”,默認(rèn)為同一種作物作為“第一種作物”,是一種不加以區(qū)分的抽象的“混合作物”,此外其不足之處還有3點(diǎn):(1)用來(lái)控制生長(zhǎng)的模式分4種(MAX.GROWTH、NORMAL GROWTH、LEAF FALL、NO LEAF);這是更接近于物候的表述,而非發(fā)育期的表述,過(guò)于粗糙,不能較好地體現(xiàn)某一種作物生長(zhǎng)發(fā)育的特性,也難以刻畫(huà)作物隨著發(fā)育期階段的不同生長(zhǎng)情況,模式模擬的作物生長(zhǎng)與實(shí)際情況有較大的出入。例如,農(nóng)業(yè)上抽穗開(kāi)花期是一個(gè)重要的發(fā)育期,抽穗開(kāi)花前后作物有明顯的區(qū)別,BCC_AVIM中未考慮開(kāi)花結(jié)果等過(guò)程;事實(shí)上,還有很多生長(zhǎng)發(fā)育參數(shù)也受到發(fā)育期的控制。所以,在陸面過(guò)程改進(jìn)中,首先需引入具體作物,再引入該作物的發(fā)育期算法,使模擬的作物發(fā)育過(guò)程更接近實(shí)際。(2)原模式對(duì)凈光合產(chǎn)物分配系數(shù)設(shè)定為四階段式固定參數(shù),與發(fā)育期不對(duì)應(yīng);改進(jìn)方案引入發(fā)育期的概念并由播種期和氣象要素來(lái)動(dòng)態(tài)表征發(fā)育期,然后依據(jù)發(fā)育期動(dòng)態(tài)采用來(lái)源于大量實(shí)測(cè)結(jié)果分配系數(shù),給每種具體作物匹配相應(yīng)的分配系數(shù),分配方式更符合作物的實(shí)際。(3)原模式中分配系數(shù)在不同階段點(diǎn)間呈明顯的階梯狀突變,改進(jìn)模式中分配系數(shù)采用了與生育期階段有關(guān)的插值函數(shù),使之日變化漸進(jìn)平緩。
東北是中國(guó)重要的商品糧生產(chǎn)基地[15],也是我國(guó)最大的玉米產(chǎn)區(qū),其產(chǎn)量占全國(guó)總產(chǎn)量1/3左右[16]。華北平原是我國(guó)重要的糧食主產(chǎn)區(qū),現(xiàn)有耕地面積32.6×106hm2,約占全國(guó)耕地面積的40%[17]。由于東北和華北在農(nóng)業(yè)生產(chǎn)中的重要地位,以及其農(nóng)業(yè)特色下墊面對(duì)陸面模式的重要影響,考慮首先在東北和華北地區(qū)選取研究區(qū)并進(jìn)行改進(jìn)。
在東北地區(qū)玉米主產(chǎn)區(qū)的模擬中,將“Crop”改為東北春玉米,并采用東北春玉米的作物生長(zhǎng)發(fā)育和分配參數(shù);在華北地區(qū)模擬中,將“Crop”改為華北冬小麥和華北夏玉米,以期達(dá)到“冬小麥-夏玉米”輪作代替“Crop”單季生長(zhǎng)的目的。
引入具體作物的發(fā)育期和分配系數(shù)等參數(shù),在計(jì)算上只增加了積溫計(jì)算和插值計(jì)算,因而,理論上并不會(huì)造成陸面模式BCC_AVIM的運(yùn)行速度減慢,實(shí)際運(yùn)行測(cè)試表明也確實(shí)如此。
采用具體的作物代替原模式中的抽象“crop”,并引入發(fā)育期控制和按照發(fā)育期采用光合產(chǎn)物分配系數(shù)后,在華北研究區(qū)(110°~120°E,31°~40°N)和東北研究區(qū)(120°~125°E,40°~50°N)分別進(jìn)行控制試驗(yàn),2000年的葉面積指數(shù)模擬結(jié)果如圖1所示。基于GLASS(global land surface satellite)葉面積指數(shù)產(chǎn)品[18]數(shù)據(jù),得到了2000年?yáng)|北、華北研究區(qū)內(nèi)葉面積指數(shù)全年變化情況(圖2),下面將以此作為實(shí)際葉面積指數(shù)的參考數(shù)據(jù)。
圖1 2000年?yáng)|北、華北研究區(qū)內(nèi)葉面積指數(shù)模擬結(jié)果
圖2 2000年?yáng)|北、華北研究區(qū)內(nèi)葉面積指數(shù)全年變化情況
在東北研究區(qū)內(nèi)引入發(fā)育期期前后,BCC_ AVIM模擬的葉面積指數(shù)均呈現(xiàn)單峰結(jié)果(圖1),均在8月達(dá)到峰值,符合玉米生長(zhǎng)的特點(diǎn);原模式模擬出的峰值葉面積指數(shù)達(dá)到2.94,而引入發(fā)育期后模式得到的結(jié)果為2.73,GLASS[18]給出的實(shí)際葉面積指數(shù)峰值為2.82。實(shí)際葉面積指數(shù)數(shù)值低于原模式的峰值結(jié)果,且更接近于利用實(shí)際作物引入發(fā)育期后的模擬結(jié)果,因此,從逐月的時(shí)間變化來(lái)看,引入發(fā)育期后,BCC_AVIM模式對(duì)于東北地區(qū)葉面積指數(shù)的年變化特征,有更好的模擬效果,該改進(jìn)方式在東北地區(qū)較為成功,值得進(jìn)一步地進(jìn)行深入研究,擴(kuò)大使用區(qū)域。
2.3.1 葉面積指數(shù)的時(shí)間變化特點(diǎn)
在華北研究區(qū),葉面積指數(shù)隨時(shí)間變化曲線具有明顯的雙峰(圖2),與該地區(qū)“冬小麥-夏玉米”的一年兩熟的種植制度相符,而引入發(fā)育期后模擬結(jié)果(圖1)仍未明顯地體現(xiàn)兩季發(fā)育的特點(diǎn)。引入發(fā)育期后的模擬結(jié)果在1—2月與原模式模擬結(jié)果一致,3—5月和12月略高于原模式(圖3),而6—11月則明顯偏低。
在3—5月引入發(fā)育期后的模擬結(jié)果稍有偏高,恰好反映了引入發(fā)育期控制的模式在春季冬小麥越冬之后進(jìn)入快速生長(zhǎng)時(shí)期,對(duì)模式葉面積指數(shù)部分存在改進(jìn)作用,但該改進(jìn)效果并未達(dá)到預(yù)期,4月中旬達(dá)到的第一輪峰值沒(méi)有得到體現(xiàn)。另外,華北研究區(qū)7—8月的實(shí)際葉面積指數(shù)峰值可以達(dá)到2.74,而原模式模擬結(jié)果為2.69,引入發(fā)育期后的模式模擬結(jié)果僅為2.46,存在一定偏差。
圖3 2000年3月、4月、5月和12月華北研究區(qū)內(nèi)葉面積指數(shù)模擬結(jié)果對(duì)比
經(jīng)過(guò)研究分析,造成引入發(fā)育期后輸出結(jié)果仍不甚理想的原因可能在于BCC_AVIM2.0模式在計(jì)算光合作用時(shí),采用物候控制,認(rèn)為冬季溫度低于閾值時(shí)不生長(zhǎng),且春季沒(méi)有在上年作物越冬的基礎(chǔ)上生長(zhǎng),而是作為另一種作物才開(kāi)始發(fā)芽;之后5—6月沒(méi)有反映冬小麥成熟收割農(nóng)事活動(dòng)以及夏玉米播種初期的低葉面積指數(shù)的生長(zhǎng)狀況,所以沒(méi)有反映出“冬小麥-夏玉米”一年兩熟的特點(diǎn)。
華北實(shí)際種植作物復(fù)雜,除了參與模擬的冬小麥、夏玉米外,還有棉花、花生、高粱、谷子、大豆、煙草、甘薯以及多種蔬菜等其他作物;華北地區(qū)實(shí)際種植方式有間作、套作等,合理的播種方式可以增大葉面積指數(shù)、提高產(chǎn)量。作物品種和種植方式共同導(dǎo)致了引入發(fā)育期后模式葉面積指數(shù)的峰值偏低。
事實(shí)上,3—5月引入發(fā)育期后的模擬結(jié)果略高于原模式模擬結(jié)果,已反映了模式改進(jìn)在春季冬小麥快速生長(zhǎng)時(shí)期,對(duì)模式葉面積指數(shù)存在值得肯定的改進(jìn)作用,但對(duì)冬小麥光合作用低溫控制物候沒(méi)有改進(jìn),使得改進(jìn)效果大打折扣。
2.3.2 葉面積指數(shù)的空間分布
將引入發(fā)育期后模擬的葉面積指數(shù)結(jié)果減去原模擬的葉面積指數(shù)結(jié)果,得到的2種方案在華北研究區(qū)內(nèi)的模擬結(jié)果差值(圖4)。根據(jù)其空間分布特點(diǎn),重點(diǎn)研究華北地區(qū)造成引入發(fā)育期模式后模擬結(jié)果在4—5月偏高和在7—9月偏低的原因。
圖4 2000年華北研究區(qū)內(nèi)引入發(fā)育期后葉面積指數(shù)模擬結(jié)果與原模擬結(jié)果差值
4—5月引入發(fā)育期后模擬葉面積指數(shù)結(jié)果偏高于原模式的位置主要分布在34°N以北的區(qū)域,且其分布方式有較明顯的緯向分布特點(diǎn)。由于在實(shí)際生產(chǎn)過(guò)程中,冬小麥順利越冬后,生長(zhǎng)速度會(huì)比較快,因而此時(shí)的作物出現(xiàn)生長(zhǎng)速度較快的特點(diǎn),引入發(fā)育期后的模式葉面積指數(shù)在此時(shí)略高于原模式,更符合實(shí)際情況。
6月起,冬小麥?zhǔn)崭詈?,在華北研究區(qū)引入發(fā)育期后模式模擬出的葉面積指數(shù)低于原模式,且偏低區(qū)域位置自南而北移動(dòng)。另外,由于原模式設(shè)定同一溫度閾值決定物候,造成的結(jié)果是江蘇、河南緯度偏南的作物先成熟,往北逐漸推遲;而事實(shí)上,由于品種本地適應(yīng)的原因,以及人為控制播種、收割時(shí)間,實(shí)際作物生長(zhǎng)對(duì)南北溫差的差異性表象并沒(méi)有那么明顯。因此,引入發(fā)育期前后模式之差具有緯向分布特點(diǎn),說(shuō)明了引入發(fā)育期控制更好地反映了生長(zhǎng)實(shí)際,訂正了原模式南部過(guò)早的成熟特征,更符合實(shí)際中作物生長(zhǎng)的特點(diǎn)[19]。
有實(shí)驗(yàn)結(jié)果表明,BCC_AVIM2.0模式性能穩(wěn)定,對(duì)全球植被的季節(jié)變化、植被凈初級(jí)生產(chǎn)力等模擬較好[20],可合理再現(xiàn)全球和局地的地面溫度變化特征[11]。因此,在實(shí)現(xiàn)用具體的作物改進(jìn)華北研究區(qū)內(nèi)作物生長(zhǎng)情況后,對(duì)引入發(fā)育期后模式的整體凈初級(jí)生產(chǎn)力(NPP)變化特征、不同器官的生物量(葉、莖、根的NPP)和2 m溫度(距離地面2 m高度處的溫度)的改進(jìn)情況進(jìn)行分析,討論引入發(fā)育期對(duì)以上重要輸出結(jié)果的影響。
2.4.1 引入發(fā)育期前后NPP變化特點(diǎn)
由于受到BCC_AVIM2.0模式自身物候限制,研究NPP時(shí)發(fā)現(xiàn),春季,NPP有分散的偏大偏小區(qū)域;夏季,引入發(fā)育期后的模式在作物占比在70%以上的區(qū)域表現(xiàn)出偏大的特點(diǎn),而在作物較少的區(qū)域偏小;冬季,華北研究區(qū)內(nèi)模擬結(jié)果無(wú)明顯變化。
有研究表明,小麥-玉米生育期的總NPP遠(yuǎn)高于其他農(nóng)業(yè)生態(tài)系統(tǒng),因而在夏季,引入發(fā)育期后NPP的模擬結(jié)果在作物占比較大的區(qū)域表現(xiàn)出偏大的特點(diǎn),展現(xiàn)了模式較好的改進(jìn)效果;而10月的異常偏小,考慮是受到了模式本身物候的限制,夏玉米播種后生長(zhǎng)緩慢。綜合來(lái)說(shuō),華北地區(qū)應(yīng)當(dāng)有較高的NPP,而引入發(fā)育期后在夏季的作物區(qū)體現(xiàn)了這一特點(diǎn),引入發(fā)育期對(duì)反映華北研究區(qū)內(nèi)實(shí)際的NPP有一定的作用。
2.4.2 引入發(fā)育期前后葉、莖、根的NPP
利用BCC_AVIM輸出的葉、莖、根的NPP,分別在華北研究區(qū)內(nèi)進(jìn)行了各月的平均(圖5),發(fā)現(xiàn)葉、莖、根的NPP在1—5月并無(wú)明顯差異;由于在調(diào)整模式的過(guò)程中,增加了“開(kāi)花期”這一特殊發(fā)育階段,由于開(kāi)花期前為營(yíng)養(yǎng)生長(zhǎng)階段,開(kāi)花后進(jìn)入生殖生長(zhǎng)階段,在本文中基于原模式的框架并沒(méi)有專(zhuān)門(mén)的種子部分的專(zhuān)門(mén)儲(chǔ)藏器官來(lái)接受并形成產(chǎn)量,所以把它加入到莖部,因而葉、根的NPP在6—10月表現(xiàn)出引入發(fā)育期后的模式低于原模式的特點(diǎn),更符合實(shí)際;莖在6—11月均出現(xiàn)引入發(fā)育期后偏大的現(xiàn)象,這個(gè)特點(diǎn)在東北研究區(qū)中也有較好的體現(xiàn),這種偏大更符合實(shí)際夏玉米干物質(zhì)積累的特點(diǎn);仔細(xì)觀察根的NPP可以發(fā)現(xiàn),原模式為單峰生長(zhǎng),而引入發(fā)育期后的模式出現(xiàn)了弱雙峰的特點(diǎn),這更符合玉米生長(zhǎng)時(shí)夏季根系分配系數(shù)減小的實(shí)際生長(zhǎng)情況,也可能是華北研究區(qū)內(nèi)可以實(shí)現(xiàn)兩季輪作的反映,根的弱雙峰特點(diǎn)更確認(rèn)了引入發(fā)育期后模式可以更好模擬華北研究區(qū)內(nèi)兩季輪作的假設(shè)。
圖5 2000年華北研究區(qū)內(nèi)各月根、莖、葉的凈初級(jí)生產(chǎn)力
2.4.3 引入發(fā)育期前后2 m溫度的變化特點(diǎn)
圖6給出了華北研究區(qū)內(nèi)引入發(fā)育期后2 m溫度模擬結(jié)果減去原模擬結(jié)果的差值空間分布,2 m溫度與葉面積指數(shù)存在一定的負(fù)相關(guān)關(guān)系。植被對(duì)地表溫度具有調(diào)控作用,增加植被覆蓋可使地表溫度變幅減小。因此,冬季葉面積指數(shù)偏大的位置2 m溫度偏高,夏季葉面積指數(shù)偏小的區(qū)域2 m溫度偏高,這符合植被對(duì)溫度的調(diào)控作用。并且,這種差異在9月份達(dá)到最大,這與引入發(fā)育期后葉面積指數(shù)模擬結(jié)果在9月份偏小最明顯的情況是一致的。
從整個(gè)華北研究區(qū)來(lái)看,在模式的溫度等原始輸入場(chǎng)無(wú)變化的情況下,由于對(duì)模式中作物生長(zhǎng)過(guò)程和分配系數(shù)的調(diào)整,引入發(fā)育期后的模式運(yùn)行出不同于原模式的葉面積指數(shù)結(jié)果,并且該結(jié)果加入了對(duì)2 m溫度的計(jì)算中,產(chǎn)生了符合實(shí)際的反饋機(jī)制。
在BCC_AVIM2.0中引入具體作物發(fā)育期和隨發(fā)育期動(dòng)態(tài)變化的分配系數(shù)后,東北研究區(qū)內(nèi)有很好的模擬效果;而華北研究區(qū)內(nèi),模式實(shí)現(xiàn)了作物跟隨發(fā)育期生長(zhǎng),并反饋于其他輸出要素,但受其他因素限制,葉面積指數(shù)未表現(xiàn)出明顯的雙峰特點(diǎn)。引入發(fā)育期后葉面積指數(shù)模擬結(jié)果在3—5月略高于原模式,反映了冬小麥越冬后快速生長(zhǎng)的特點(diǎn),其效果值得肯定;引入發(fā)育期后作物生長(zhǎng)的南北差異變小,是對(duì)原模式受溫度控制為主的一種訂正,更加符合實(shí)際華北地區(qū)作物生長(zhǎng)的特點(diǎn)。在夏季,引入發(fā)育期后的NPP模式結(jié)果在作物占比較大的區(qū)域表現(xiàn)出偏大的特點(diǎn),更符合實(shí)際,展現(xiàn)了模式較好的改進(jìn)效果。引入發(fā)育期后模式輸出的葉、根的NPP在6—10月明顯低于原模式,而莖的NPP在6—11月偏大,表示“開(kāi)花期”這一特殊發(fā)育階段成功加入了作物生長(zhǎng)的過(guò)程中;且根的NPP出現(xiàn)了弱雙峰的特點(diǎn),表示引入發(fā)育期后模式可以更好模擬華北研究區(qū)內(nèi)兩季輪作。在原始輸入場(chǎng)無(wú)變化的情況下,引入發(fā)育期后模式模擬的作物生長(zhǎng)過(guò)程發(fā)生了變化,并且該變化對(duì)2 m溫度的計(jì)算產(chǎn)生了符合實(shí)際的反饋,成功作用于整個(gè)模式模擬的過(guò)程中。
圖6 2000年華北研究區(qū)內(nèi)引入發(fā)育期后2 m溫度模擬結(jié)果與原模擬結(jié)果差值
當(dāng)前的BCC_AVIM2.0已經(jīng)可以實(shí)現(xiàn)對(duì)陸面過(guò)程很好的模擬效果,并且廣泛應(yīng)用于相關(guān)領(lǐng)域。由于各領(lǐng)域日益增長(zhǎng)的對(duì)模式的需求,需進(jìn)一步提高模式的精細(xì)程度。本研究已實(shí)現(xiàn)在華北研究區(qū)內(nèi)加入冬小麥、夏玉米2種作物,并引入發(fā)育期,使作物跟隨各發(fā)育期不同的分配系數(shù)動(dòng)態(tài)生長(zhǎng)。由研究結(jié)果可見(jiàn),改進(jìn)已初見(jiàn)成效,輸出結(jié)果可以作用于模式其他變量。下一步仍需繼續(xù)對(duì)BCC_AVIM2.0模式進(jìn)行改進(jìn)研究,實(shí)現(xiàn)作物具體化,增加作物種類(lèi)和熟制,調(diào)整模式的物候設(shè)定,通過(guò)空間位置設(shè)置具體閾值等方法,改進(jìn)越冬作物的生長(zhǎng);繼續(xù)調(diào)整和改進(jìn)光合產(chǎn)物分配系數(shù),讓分配系數(shù)由作物決定并具有地方特色,進(jìn)而使作物生長(zhǎng)更加符合實(shí)際,得到更好的模擬結(jié)果。