單成鋼,王憲昌,張教洪,馬炳鵬,李桂祥,倪大鵬,朱彥威,王志芬
(1.山東省農(nóng)業(yè)科學(xué)院農(nóng)產(chǎn)品研究所/藥用植物研究中心/農(nóng)業(yè)農(nóng)村部新食品資源加工重點(diǎn)實(shí)驗(yàn)室,山東 濟(jì)南 250100;2.沂南駿達(dá)農(nóng)業(yè)發(fā)展有限公司,山東 沂南 276314;3.山東省果樹研究所,山東 泰安 271000)
光照是植物進(jìn)行光合作用的基礎(chǔ),影響著植物光合作用過程中的同化力形成、酶活化、氣孔開放等。遮蔭是使植物接收的光量子密度受到限制的環(huán)境或人為因子[1],一般會(huì)導(dǎo)致植物的日照時(shí)間減少或者照射強(qiáng)度降低。遮蔭一方面使植物生長環(huán)境的小氣候產(chǎn)生變化,另一方面對植物的生長、發(fā)育、生理生化特性以及礦質(zhì)營養(yǎng)代謝產(chǎn)生影響,這種影響既有有利的一面也有不利的一面,而且影響程度隨植物自身的差異性而有所不同[2]。
近年來,隨著對清潔能源需求增加和光伏扶貧工程的推進(jìn),我國的光伏產(chǎn)業(yè)得到迅猛發(fā)展。截至2016年底,我國光伏發(fā)電新增裝機(jī)容量3 454萬千瓦,累計(jì)裝機(jī)容量7 742萬千瓦,新增和累計(jì)裝機(jī)容量均為全球第一[3]。為了合理利用土地和集約利用空間,催生出了光伏與農(nóng)業(yè)相結(jié)合的產(chǎn)業(yè)形態(tài)——光伏農(nóng)業(yè),包括光伏大棚、露地光伏農(nóng)業(yè)等形式,無論哪種形式均是為了能夠在滿足光伏發(fā)電的同時(shí)開展有效的農(nóng)業(yè)生產(chǎn)。但在這種復(fù)合模式下植物生長環(huán)境如光照、溫度、濕度等將產(chǎn)生什么變化,這種變化對生長期間的植物有多大影響,目前鮮見研究報(bào)道[4,5]。
本文著重研究了露地光伏-丹參生長模式下的光照變化,以及這種變化對丹參生長的影響,以期為光伏條件下開展合理的農(nóng)業(yè)生產(chǎn)提供科學(xué)依據(jù)。
試驗(yàn)在沂南縣岸堤鎮(zhèn)興旺莊村瑞諾多晶硅光伏電站進(jìn)行,多晶硅陣列寬3 m,骨架間距6.6 m,丹參品種為魯?shù)?號,種植時(shí)間是2017年3月20日。采用大壟雙行種植,每一陣列行間種植4壟丹參,每壟種植2行,壟寬0.9 m,壟高0.3 m,統(tǒng)一常規(guī)水肥管理。光伏板與種植模式布局見圖1。
注:γ=25°;h1=1.4 m; H=2.9 m;h2=2.2 m; d=6.6 m;e=1.5 m; f=0.9 m。
試驗(yàn)地點(diǎn)的太陽高度角和方位角可以通過計(jì)算得出[6]。
1.2.1 太陽赤緯的計(jì)算 太陽赤緯隨季節(jié)而變化,春分和秋分時(shí)太陽赤緯為0°,夏至和冬至分別為23°23′。其計(jì)算式[7]如式(1)。
δ=0.3730+23.2567sinθ+0.1149sin2θ-0.1717sin3θ-0.7580cosθ+0.3650cos2θ+0.0201cos3θ
(1)
式中:δ為太陽赤緯;θ為日角,即θ=2πt/365.2422,其中t=N-N0,N為積日,即日期在年內(nèi)的順序號,N0=79.6764+0.2422×(年份-1985)-INT[(年份-1985)/4]。
1.2.2 太陽高度角的計(jì)算 太陽高度角是指某地太陽光線與通過該地與地心相連的地表切面的夾角,簡單來說即太陽光的入射方向與地平面之間的夾角,用α表示,計(jì)算式[8]如式(2)。
sinα=sinφsinδ+cosφcosδcosω
(2)
式中:δ為太陽赤緯;φ為觀測地地理緯度;ω為地方時(shí)角,ω=15°(t-12),午后為正,午前為負(fù)。
1.2.3 太陽方位角的計(jì)算 太陽方位角是以目標(biāo)物的北方向?yàn)槠鹗挤较颍蕴柟獾娜肷浞较驗(yàn)榻K止方向,按順時(shí)針方向所測量的角度,用β表示。中國區(qū)域太陽方位角范圍在90°~270°,正午為180°左右。其計(jì)算公式[8]如式(3)。
cosβ=(sinαsinφ-sinδ)cosαcosφ
(3)
式中:δ為太陽赤緯;φ為觀測地地理緯度;α為太陽高度角。
1.3.1 光伏板在地面陰影長度 即光伏板對陽光遮擋所形成的陰影長度。其計(jì)算公式如式(4)。
(4)
其中:L為光伏板地面陰影長度;H為光伏板高度;α為太陽高度角。
1.3.2 有效陰影長度 即光伏板在南北(與行向垂直)方向上的陰影長度。此陰影長度對位于行間的作物形成有效的遮蔭。其計(jì)算公式如式(5)。
(5)
其中:L′為有效陰影長度;L為光伏板地面陰影長度;β為太陽方位角。
1.4.1 株高測定 開花期每小區(qū)選取代表性植株10株,測量株高。
1.4.2 凈光合速率測定 利用 LI-6400XT 便攜式光合作用測定系統(tǒng),在丹參開花期(6月15日)8—16時(shí)每間隔2小時(shí)測定一次。每處理選擇長勢具有代表性植株10株,在每株相同部位選取3個(gè)葉片,測定復(fù)葉的中央小葉,取3個(gè)讀數(shù)的平均值作為該株的測定結(jié)果。
1.4.3 SPAD值測定 使用SPAD-502型葉綠素儀進(jìn)行測定。 選擇健康成熟葉片的中央小葉, 分別在小葉葉基、葉中、葉尖處測得SPAD值, 求平均作為每片葉的值。每株測定3個(gè)小葉,求平均作為每個(gè)植株的值,每處理測定10株。
試驗(yàn)數(shù)據(jù)采用Microsoft Excel 軟件進(jìn)行處理,采用SPSS統(tǒng)計(jì)軟件進(jìn)行差異顯著性分析。
由表1可見,在6月中旬,太陽高度角從6時(shí)的11.93°逐步升高至12時(shí)的77.53°,而后逐步降低,至18時(shí)降至14.93°。6月中旬日出時(shí)太陽位于東北方向,角度為69.7°,至9時(shí)左右太陽方位角達(dá)到90°,即位于正東方,至12時(shí)方位角為171.65°,接近正南方,而后逐步向西移動(dòng),16時(shí)左右到達(dá)正西方,而后繼續(xù)向西北方向移動(dòng),18時(shí)與正北方夾角為288.28°。
由于受到太陽高度角和方位角的共同作用,有效陰影長度也呈現(xiàn)有規(guī)律的變化,日出時(shí)因?yàn)樘栁挥跂|北方向,陰影來自于北側(cè)光伏陣列的遮擋,太陽高度角小,陰影比較長,約4.76 m;隨著太陽的南移和升高,陰影長度迅速縮短,至9時(shí),太陽方位到達(dá)正東方附近,有效陰影長度僅為0.17 m;至正午時(shí)分,太陽到達(dá)正南方,太陽高度角77.53°,有效陰影長度為0.63 m;正午之后的有效陰影變化趨勢與上午正好對稱,先由大變小,再由小變大,見表1。
表1 2017年6月15日太陽高度角、方位角和有效陰影長度的變化
以表1的時(shí)刻為橫坐標(biāo),有效陰影長度為縱坐標(biāo)繪制圖2,可見圖中存在4個(gè)單調(diào)區(qū)間,分別對應(yīng)方位角0°~90°,90°~180°,180°~270°和270°~360°四個(gè)區(qū)間。
圖2 有效陰影長度隨時(shí)間變化
在每一個(gè)單調(diào)區(qū)間內(nèi),以有效陰影長度為自變量,以對應(yīng)時(shí)間點(diǎn)為因變量,進(jìn)行曲線擬合,可得每個(gè)區(qū)間的擬合曲線,其中,對應(yīng)時(shí)間點(diǎn)換算為數(shù)字形式,例如6時(shí)30分,換算為6.50,6時(shí)15分換算為6.25。所得擬合曲線見表2。
表2 每個(gè)單調(diào)區(qū)間的擬合曲線及決定系數(shù)
利用表2中曲線1可以計(jì)算得到在太陽方位角0°~90°區(qū)間,即日出后不久,每個(gè)處理丹參壟的日照時(shí)間點(diǎn),并由此計(jì)算出該時(shí)段的日照時(shí)數(shù);同理,利用曲線4計(jì)算可得在方位角270°~360°區(qū)間,即日落前不久,每一壟丹參的日照時(shí)數(shù),這兩個(gè)時(shí)段的陰影來自丹參種植區(qū)北側(cè)光伏板的遮擋。而在太陽方位角90°~180°和180°~270°區(qū)間,陰影來自丹參種植區(qū)南側(cè)光伏板的遮擋,從表1的數(shù)據(jù)分析出這兩個(gè)時(shí)段遮蔭范圍較少,主要是對R1壟的遮擋,分別分析R1壟內(nèi)兩行,遠(yuǎn)離光伏板的一行始終不被遮蔭,而靠近光伏板的一行經(jīng)計(jì)算上午和下午分別被遮陰2小時(shí)和1.5小時(shí),因此,R1壟的光照時(shí)間采用兩行光照時(shí)間的平均值。由此可得6月15日各處理全天日照時(shí)數(shù),見表3。
考慮到丹參的整個(gè)生育期,分別計(jì)算2017年3月20日(春分日)、4月20日、5月20日、6月21日(夏至日)、7月20日、8月20日、9月23日(秋分日)、10月20日、11月20日每個(gè)遮蔭處理的日照時(shí)數(shù)。
表3 各處理丹參在一天中日照時(shí)數(shù) (h)
注:表中距離為每一丹參壟南側(cè)邊界與相鄰北側(cè)光伏板高點(diǎn)的地面投影間的水平距離,表示該壟接受來自北側(cè)光照的臨界點(diǎn)。
由圖3可見,在無遮蔭條件下(CK),日照時(shí)數(shù)隨季節(jié)呈現(xiàn)有規(guī)律的變化,從春分日的12.00 h逐漸增至夏至日的17.15 h,特別是5月下旬至6月下旬,增加幅度最大,日照時(shí)數(shù)增加了25%;夏至日后,日照時(shí)數(shù)迅速下降,與之前的上升趨勢幾乎呈對稱狀態(tài),至秋分日又回落到12 h;之后繼續(xù)下降,11月20日降至9.75 h。與對照相比,各遮蔭處理的日照時(shí)數(shù)都受到不同程度的影響。R1是受遮蔭最為嚴(yán)重的處理,春分日的光照時(shí)間不足15 min,之后逐漸增大,至夏至日達(dá)到最大值10.43 h,之后又逐漸減少,秋分日之后基本接收不到太陽直射,全天處于遮蔭狀態(tài)。R2在3月20日日照時(shí)數(shù)僅為30 min,在隨后的1個(gè)月內(nèi)迅速達(dá)到12 h,一直持續(xù)到8月20日,日照時(shí)數(shù)均保持在12 h左右;但到了9月23日(秋分日),日照時(shí)數(shù)又迅速回落到30 min,之后的生育期內(nèi)都處于遮蔭狀態(tài)。R3自春分日開始一直到秋分日,日照情況相對比較平穩(wěn),保持在12 h左右,6月20日到8月20日略低于12 h,主要因?yàn)楸眰?cè)光伏板的遮擋,秋分日之后到10月20日,陽光逐步被完全遮擋。R4與R3情況相似,但8月20日前遮蔭較R3更為嚴(yán)重,而陽光直射可以一直持續(xù)到10月中下旬,10月20日尚有將近7 h的光照,這是因?yàn)镽4雖然更容易被北側(cè)光伏板遮擋,但距離南側(cè)光伏板最遠(yuǎn),遮光最為輕微。
為了更加清晰地比較分析5個(gè)處理生育期內(nèi)的總?cè)照諘r(shí)數(shù),可以通過積分分別計(jì)算5條曲線與橫軸之間構(gòu)成圖形的面積來求取,結(jié)果見圖4。
圖3 丹參生育期內(nèi)的日照時(shí)數(shù)變化
圖4 各處理生育期內(nèi)的總?cè)照諘r(shí)數(shù)
四個(gè)處理與對照的總?cè)照諘r(shí)數(shù)存在顯著差異,隨著與南側(cè)光伏板的距離增大總?cè)照諘r(shí)數(shù)逐漸增大,呈線性變化趨勢,其決定系數(shù)R2=0.972,表明與光伏板的距離與日照時(shí)數(shù)間存在顯著相關(guān)性。對照的總?cè)照諘r(shí)數(shù)最大,為1 045 h,R1的日照時(shí)數(shù)最小,為400 h,僅占對照的38%;R2總?cè)照諘r(shí)數(shù)明顯高于R1,約為610 h,高出R1約52%,占對照的58%;R3和R4雖然從4月至8月每日的日照時(shí)數(shù)略低于R2,但由于他們在生育期前期和后期的光照時(shí)間明顯高于R2,光照天數(shù)明顯多于R2,因此,總體來看R3和R4的總?cè)照諘r(shí)數(shù)分別高出R2約19%和32%,占對照的69%和77%。
由圖5可以看出,遮蔭對丹參株高有著顯著的影響,R1遮光最嚴(yán)重,株高平均65 cm,高出對照約18%,而R2平均株高60 cm,高出對照9%。在0.05水平下R1、R2差異顯著,且均顯著高于其他處理,R3、R4與對照沒有顯著差異。
由圖6可知,8—12時(shí),各處理的凈光合速率均呈現(xiàn)上升趨勢,12時(shí)達(dá)最大值;12—16時(shí),4個(gè)遮蔭處理均呈下降趨勢,而對照的凈光合速率14—16時(shí)又升高。遮蔭降低了丹參的光合速率。8時(shí)對照的凈光合速率顯著高于其他處理;10時(shí),R3、R4的光合速率迅速升高,但仍略低于對照,差異不顯著,而R1和R2顯著低于對照,分別比對照低78%和71%;12時(shí),R1和R2的凈光合速率也表現(xiàn)出顯著差異,但仍明顯低于對照和另外兩個(gè)處理,分別比對照低69%和37%;14時(shí)各處理的表現(xiàn)與10時(shí)相似,16時(shí)各遮蔭處理的光合速率均降至極低狀態(tài),顯著低于對照。
圖5 遮蔭對株高的影響
圖6 遮蔭對光合速率的影響
由圖7可知,遮蔭顯著降低丹參葉片的SPAD值,遮蔭越重,SPAD值越低。R1和R2的SPAD值顯著低于對照,分別比對照低7.6%和4.9%;R3和R4略低于對照,分別比對照低1.9%和2.1%,與對照間差異不顯著。
如表4所示,CK根鮮重、根干重、根條數(shù)、根總長、產(chǎn)量值最大,R3的根鮮重與CK差異顯著,其他指標(biāo)R3、R4均與CK差異不顯著,其產(chǎn)量分別比對照低11.97%和3.61%。除根總長、根直徑外,R2的其他指標(biāo)均與對照差異顯著,其中產(chǎn)量約低于對照30%,根干重、根條數(shù)、折干率分別比對照低31%、30%、14%。R1的所有根部性狀指標(biāo)與對照均有顯著差異,其中產(chǎn)量約降低50%,根鮮重、干重、根條數(shù)、根總長、根直徑、折干率分別比對照低35%、49%、42%、20%、10%、25%。
圖7 遮蔭對SPAD值的影響
表4 各處理丹參根部性狀和產(chǎn)量
注:同列不同小寫字母表示處理間差異達(dá)0.05顯著性水平。
研究表明,利用天文學(xué)公式計(jì)算太陽方位角和高度角,對樹木樹冠進(jìn)行模型化分析,計(jì)算得出不同經(jīng)緯度某一時(shí)刻樹木的投影面積,可以確定果樹合理的種植間距[11,12]。建筑物的遮蔭效應(yīng),夏季以建筑物北側(cè)日照條件最好,冬季以南側(cè)日照條件最好,全年綜合情況則以南側(cè)最優(yōu),北側(cè)最差[13]。本研究在計(jì)算得到太陽方位角和高度角的基礎(chǔ)上,推算出光伏板在一天之中以及丹參整個(gè)生長季的遮蔭狀況,進(jìn)一步得出距離南側(cè)光伏板由近及遠(yuǎn)的丹參壟,在生育期內(nèi)接收的日照總時(shí)數(shù)分別為對照的38%、58%、69%和77%,距離和遮蔭程度顯著相關(guān),這一結(jié)論與建筑物的遮蔭情況相一致[13]。
除了光環(huán)境的改變,光伏板遮蔭也能夠改變植被的其他生長條件。能夠有效降低土壤pH值和電導(dǎo)率,導(dǎo)致土壤水分的入滲量和蒸發(fā)量都比較小,水分上下流通不暢,不利于水分通道的形成,可能導(dǎo)致土壤容重增加。光伏板遮蔭作用為更多耐蔭性植物生長提供了可能,從而增加了物種豐富度、均勻度和優(yōu)勢度[14]。光伏電廠干擾下高寒荒漠草原區(qū)植被群落的植物種數(shù)、蓋度、地上生物量、豐富度指數(shù)、土壤含水量、有機(jī)質(zhì)和全氮含量升高[15]。
這些條件的改變使正處于生長期間的植物受到顯著影響。例如:遮蔭顯著降低甘草株高,影響光合作用,造成顯著減產(chǎn)[5]。不同強(qiáng)度的遮光可顯著減弱作物的光合作用,使單葉凈光合速率下降[16],氣孔關(guān)閉[17],從而導(dǎo)致干物質(zhì)積累量和產(chǎn)量下降;不同品種對遮蔭的適應(yīng)能力存在差異[4,18-21。遮蔭會(huì)使喜光植物葉片變薄,顯著降低其葉片SPAD值、葉綠素a/b值和最大光合速率,導(dǎo)致CO2同化速率和氮的吸收率降低,從而延緩其生長[23-25]。而耐陰或喜陰植物在遮蔭情況下,葉綠素質(zhì)量分?jǐn)?shù)均較對照極顯著升高,可明顯增強(qiáng)光系統(tǒng)對低密度光量子的吸收、傳遞和轉(zhuǎn)化[26-28]。不同遮蔭條件對不同地被植物根鮮重、分株數(shù)有不同程度的影響:金葉過路黃植株根鮮重、分株數(shù)均在27%的遮光度下最高;射干植株鮮重及高度隨遮光度增大呈下降趨勢,根鮮重及分枝數(shù)變化緩慢[29]。
丹參為喜光植物,具有較高的光飽和點(diǎn)和比較低的光補(bǔ)償點(diǎn),對遮蔭環(huán)境具有一定的適應(yīng)能力[30,31]。本研究結(jié)果也顯示,生育期內(nèi)總遮光度62%,丹參仍能夠完成正常的生命周期,并且能夠達(dá)到對照產(chǎn)量的50%,凈光合速率、株高、SPAD值以及根部性狀與對照均有顯著差異;遮光42%,約能達(dá)到對照產(chǎn)量的70%;遮光度<31%,產(chǎn)量基本不受影響。
本研究為光伏農(nóng)業(yè)模式提供了一種遮蔭程度的計(jì)算方法,可以根據(jù)光伏板的高度、布局等準(zhǔn)確推算出每日和整個(gè)生育期的日照時(shí)數(shù),這為光伏農(nóng)業(yè)電站的設(shè)計(jì)、種植品種的選擇和搭配提供了參考依據(jù)。在本研究的光伏-丹參生產(chǎn)模式下,只要保證全生育期累計(jì)日照時(shí)數(shù)720 h以上,就可以基本保證其產(chǎn)量;即便全生育期日照時(shí)數(shù)只有400 h,也可以帶來正常產(chǎn)量50%的收益。這對保障光伏農(nóng)業(yè)電站的整體效益具有指導(dǎo)價(jià)值。