白昕欣, 黃蕭霖, 秦孟晟, 張岳芳, 陳留根, 郭 智, 郝 璐**

(1.南京信息工程大學(xué)江蘇省農(nóng)業(yè)氣象重點實驗室 南京 210044; 2.江蘇省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與環(huán)境研究所 南京 210014)

亞熱帶地區(qū)水稻田地表反照率變化特征*

白昕欣1, 黃蕭霖1, 秦孟晟1, 張岳芳2, 陳留根2, 郭 智2, 郝 璐1**

(1.南京信息工程大學(xué)江蘇省農(nóng)業(yè)氣象重點實驗室 南京 210044; 2.江蘇省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與環(huán)境研究所 南京 210014)

利用地面實測資料研究稻田地表反照率, 一方面可以更好地刻畫以稻田為主要土地利用方式的流域地氣之間的能量分配過程; 另一方面, 可以為陸面模式提供更為準(zhǔn)確的參數(shù)值, 以及為遙感反演的地表反照率提供驗證, 從而為更好地解釋土地利用/覆被變化對全球氣候變化的影響機制提供參考。本文利用江蘇省農(nóng)業(yè)科學(xué)院溧水試驗基地四分量儀測得的2016年稻田地表反照率數(shù)據(jù), 分析了稻田地表反照率特征, 并結(jié)合同期觀測的太陽短波輻射、溫度、濕度、風(fēng)速、風(fēng)向等氣象數(shù)據(jù), 進(jìn)行相關(guān)性分析, 識別影響稻田地表反照率的主要氣象因子, 為進(jìn)一步量化地表反照率與溫度及濕度等的參數(shù)化關(guān)系提供參考。結(jié)果表明: 晴天稻田地表反照率整體上呈“U”型分布, 中午較低, 下午和上午較高。晴天稻田地表反照率在一天內(nèi)的變化呈不對稱特性,其不對稱性主要是由露水和風(fēng)速、風(fēng)向引起。太陽高度角較小時, 露水的散射作用使得上午時分的地表反照率值較下午高; 而太陽高度角較大時, 西南風(fēng)促使作物葉面傾斜, 從而使得下午的地表反照率值較上午高。稻田晴天地表反照率值較陰雨天高。地表反照率在晴天與出射短波輻射相關(guān)系數(shù)最高(0.670,P＜0.01), 在陰天與相對濕度之間的相關(guān)程度最高(-0.480,P＜0.05)。在整個觀測期間, 稻田生長季內(nèi)地表反照率呈現(xiàn)先升高后降低的趨勢, 地表反照率最高值出現(xiàn)在灌漿期到成熟期之間, 插秧到分蘗期之間最低, 其中灌漿期地表反照率與太陽短波輻射及濕度間的相關(guān)程度較高, 并且均通過了P＜0.01顯著性檢驗。分蘗期和拔節(jié)期是水稻生長季內(nèi)地表反照率變化較快的兩個生育期, 并受氣象因素的顯著影響。

水稻田; 地表反照率; 生長季; 不對稱; 氣象因子

地表反照率控制著到達(dá)地球表面的凈輻射能量[1],反映地氣之間的能量分配狀況, 是地球氣候系統(tǒng)的關(guān)鍵變量, 對地表和大氣之間的水熱交換具有重要的意義[2-4]。當(dāng)入射到地表的太陽短波輻射一定時, 地表反照率越大, 說明地表反射的太陽輻射越多, 地表吸收的能量較少; 而地表反照率越小, 說明地表反射的太陽輻射越少, 地表吸收的能量較多, 則會引起一定的增溫效應(yīng)[5]。地表反照率影響地表溫度、蒸散、光合作用和呼吸作用等生物物理及化學(xué)過程,因此地表反照率直接或間接地影響著全球及區(qū)域氣候[5]。

高分辨率的地表反照率參數(shù)在監(jiān)測降水以及農(nóng)作物災(zāi)害等突發(fā)自然過程, 作物種植、作物收獲等地區(qū)尺度的人類活動過程[6-8], 自然生態(tài)系統(tǒng)的輻射能量平衡過程[9]和城市熱島[10]等中具有重要作用。作為許多陸面氣候模式和地表能量平衡方程中的一個重要輸入?yún)?shù), 地表反照率是影響地球氣候系統(tǒng)的關(guān)鍵變量[11]。

地表反照率的大小受多種因素的共同影響。地表反照率隨太陽高度角的增大而減小, 隨土壤濕度的增加而減小[12]。陳云浩等[13]利用NOAA衛(wèi)星資料對我國西北地區(qū)的地表反照率進(jìn)行反演發(fā)現(xiàn), 在植被覆蓋度和土壤濕度狀況的共同影響下, 我國西北地區(qū)地表反照率呈現(xiàn)明顯的地帶性分布。楊娟等[10]利用MODIS地表雙向反照率產(chǎn)品MOD43B3研究了北京地區(qū)地表反照率時空分布及變化特征, 發(fā)現(xiàn)山區(qū)的地表反照率最低, 其次為北京城區(qū), 而平原區(qū)和永定河流域的地表反照率較高。蔡福等[14]用NCEP再分析資料, 發(fā)現(xiàn)一年中夏季7月的地表反照率最低, 春秋季的地表反照率較高, 而冬季 1月的地表反照率最高。遙感方法是獲取大區(qū)域地表反照率的有效方法。但是, 遙感數(shù)據(jù)反映的是地表綜合信息, 與其分辨率有關(guān), 所以遙感反演的地表反照率需要實測資料的驗證。另一方面, 當(dāng)前陸面過程模式中, 地表反照率作為診斷變量由其他參數(shù)導(dǎo)出,在時間和空間上存在較大的不確定性[2]。因此, 考慮下墊面的空間異質(zhì)性, 盡可能地開展各種不同下墊面地表反照率觀測試驗具有重要意義。

張亞峰等[15]研究了干旱區(qū)春小麥(Triticum aestivum)發(fā)育期反照率變化特征, 發(fā)現(xiàn)晴天地表反照率在每日內(nèi)的變化呈“U”型分布, 而其日變化的不對稱性主要是受上午附著在春小麥葉表面的露珠影響。余予等[16]研究了安徽省壽縣地區(qū)小麥田和水稻(Oryza sativa)田的地表反照率變化, 發(fā)現(xiàn)麥田和稻田整個生育期內(nèi)的地表反照率具有大致相同的變化趨勢。盡管前人對農(nóng)田地表反照率有了一定的研究,但是對不同氣候年型下的農(nóng)田地表反照率, 尤其是亞熱帶地區(qū)水稻田地表反照率的相關(guān)研究較少。本研究觀測了 2016年偏暖濕氣候年型下的稻田地表反照率, 分析了稻田地表反照率變化特征。

本文以亞熱帶地區(qū)秦淮河流域水稻田為研究對象, 利用江蘇省農(nóng)業(yè)科學(xué)院溧水試驗基地的四分量儀, 分析了不同天氣狀況下以及生長季內(nèi)稻田地表反照率變化特征, 并結(jié)合溫度、濕度及風(fēng)速、風(fēng)向等氣象數(shù)據(jù), 分析了其可能影響因素, 以期為陸面模式提供相關(guān)參數(shù)以及為遙感反演的地表反照率提供地面驗證。

1 研究方法

1.1 研究區(qū)概況

研究站點位于江蘇省農(nóng)業(yè)科學(xué)院溧水試驗基地,該基地位于溧水縣白馬鎮(zhèn)(119.2°E, 31.6°N, 海拔高度38 m), 總面積為 0.8 km2, 站點四周均為水稻田, 水源來自于白馬湖水庫?；镜孛差愋蜑榍鹆? 土質(zhì)類型為白漿土?；貙儆谇鼗春恿饔? 為亞熱帶季風(fēng)氣候, 雨量充沛, 年降水量為 1 048 mm, 四季分明, 年平均溫度15.4 ℃。根據(jù)溧水站多年氣象觀測數(shù)據(jù)(1985—2016年, 包括氣溫、降水量以及相對濕度)分析, 溧水2016年5月和6月的平均氣溫較歷年偏低, 而其他月份的平均氣溫較歷年偏高。2016年水稻生長季平均氣溫略高于歷年同期值。除8月份外,2016年水稻生長季內(nèi)各月的相對濕度和降水量均大于歷年同期值?？傮w上, 2016年水稻生長季屬于偏暖濕年份(表1)。觀測站點種植的水稻品種為‘南粳9108’, 2016年的插秧時間為6月17日, 6月27日進(jìn)入分蘗期, 7月15日進(jìn)入拔節(jié)期, 8月12日進(jìn)入抽穗期, 8月15日進(jìn)入灌漿期, 10月10日進(jìn)入成熟期, 10月25日收獲。

表1 研究站點(溧水)2016年氣象條件與歷年平均值(1985—2015年)比較Table 1 Monthly means of meteorological factors in 2016 and averages from 1985-2015 in the study area of Lishui

1.2 觀測儀器與數(shù)據(jù)預(yù)處理

本文的地表反照率和氣象數(shù)據(jù)來自基地安裝的四分量凈輻射傳感器(CNR4)、雨量筒(TE525MM)、紅外溫度傳感器(109-L50)以及風(fēng)速風(fēng)向傳感器等常規(guī)氣象觀測儀器。由CNR4可分別測得向上和向下的短波輻射及向上和向下的長波輻射量, 短波輻射光譜范圍為 0.31～2.80 μm, 長波輻射光譜范圍為4.20～42.00 μm。四分量凈輻射傳感器、紅外溫度傳感器安裝在稻田中央距地面2 m高度處, 風(fēng)速風(fēng)向傳感器安裝在稻田中央距地面 3 m高度處, 雨量筒水平安裝在稻田中央距地面 70 cm高度處, 數(shù)據(jù)無重復(fù)。在數(shù)據(jù)分析過程中, 剔除缺測值, 選擇白天8: 00—17: 00的半小時反照率數(shù)據(jù), 求得白天的反照率日均值, 分析稻田生長季和生長季內(nèi)水稻不同生育階段以及非生長季的反照率變化特征。選擇 4—11月之間幾個典型晴天的白天地表反照率數(shù)據(jù),分析稻田反照率的日變化特征, 再結(jié)合溫度、相對濕度、風(fēng)速風(fēng)向等氣象數(shù)據(jù)分析晴天白天反照率的不對稱現(xiàn)象。根據(jù)晴天、陰天和雨天3種不同天氣狀況下的白天地表反照率數(shù)據(jù), 分析稻田反照率在不同天氣狀況下的日變化特征。

采用作圖法分析典型晴天和不同天氣條件下以及整個生長季內(nèi)稻田地表反照率的變化特征, 用SPSS軟件的典型相關(guān)分析方法分析稻田地表反照率同短波輻射、溫度及相對濕度之間的相關(guān)關(guān)系。用MATLAB軟件處理研究站點的歷年氣象條件數(shù)據(jù)。

2 結(jié)果與分析

2.1 晴天稻田地表反照率的日變化特征

2.1.1 水稻不同發(fā)育期典型晴天的稻田地表反照率

從水稻播種之前到收獲的不同階段, 各選取一個典型晴天分析白天的地表反照率變化(圖1)。所選擇的晴天日期分別是4月29日(插秧前期)、5月16日(插秧前期)、6月17日(插秧期)、7月22日(拔節(jié)期)、8月 15日(灌漿期)、9月2日(灌漿期)、10月9日(成熟期)和11月2日(收獲后期)。

從圖 1可以看出, 稻田日地表反照率整體呈現(xiàn)中午小、上午和下午較大的“U”型分布特征, 與張亞峰等[15]的研究結(jié)果相同。在所選取的幾個典型晴天中, 地表反照率最低值均出現(xiàn)在中午左右, 而最高值在上午和下午均有分布。日反照率的最高值為0.23, 出現(xiàn)在8月15日和9月2日的上午; 最低值為0.07, 出現(xiàn)在6月17日接近中午的時候。水稻的晴天日地表反照率呈向上的拋物線主要受太陽高度角影響。地表反照率隨太陽高度角的增加而減小,中午太陽高度角較大, 上午和下午太陽高度角較小,因此一天中中午的地表反照率較低, 上午和下午的地表反照率較高。

圖1 水稻不同生育期典型晴天日(月-日)的稻田地表反照率Fig.1 Surface albedo of paddy fields in typical sunny days(month-day) in different growth periods of rice

隨著水稻生長, 地表反照率的日變化特征變得更加明顯。水稻插秧前期, 日地表反照率變化幅度較小, 最大值和最小值之間差異較小; 水稻插秧之后, 日地表反照率波動較大, 最大值與最小值之間差異明顯增大。典型晴天稻田地表反照率在中午較小, 上午和下午較大, 同太陽高度角呈相反的變化趨勢, 不同天氣條件下, 地表反照率也不同, 因此水稻生長季內(nèi)的地表反照率變化受太陽高度角和氣象條件的共同影響。

2.1.2 晴天地表反照率不對稱特征及其成因

晴天水稻田的地表反照率日變化呈現(xiàn)不對稱現(xiàn)象[15], 不對稱性由上午與下午地表反照率的差值占二者均值的百分?jǐn)?shù)表示。選取4月29日(插秧前期)、5月 16日(插秧前期)、6月 17日(插秧期)、7月 22日(拔節(jié)期)、8 月 15 日(灌漿期)、9 月 2 日(灌漿期)、10月9日(成熟期)和11月2日(收獲后期)的典型晴天分析晴天稻田地表反照率的不對稱特征(圖2), 選取 4月 29日(插秧前期)、5月 16日(插秧前期)、6月 17日(插秧期)、7月 22日(拔節(jié)期)、8月 15日(灌漿期)和9月2日(灌漿期)的典型晴天分析晴天稻田地表反照率的不對稱成因(圖3、4、5)。

地表反照率半小時均值日變化呈不對稱性的可能原因: 一是植被葉面受到風(fēng)速和風(fēng)向的影響[17-19],二是附著在植被上面的露水的散射作用使得反照率增加[20-22]。由圖2可以看出, 當(dāng)太陽高度角小于60°時, 一般上午的地表反照率大于下午; 當(dāng)太陽高度角大于60°時, 下午的地表反照率值大于上午。利用典型晴天的氣象數(shù)據(jù), 計算了各時間點的平均露點溫度。通過露點溫度和空氣實際溫度對比(圖3a)可以看出, 夜晚實際氣溫和露點溫度較為接近, 有利于露水的產(chǎn)生; 白天實際氣度和露點溫度差異較大, 不利于露水的產(chǎn)生。由于夜晚產(chǎn)生的露水的散射作用使上午的地表反照率增加, 從而導(dǎo)致了上午的地表反照率大于下午的不對稱現(xiàn)象。

由典型晴天的平均入射短波輻射和出射短波輻射圖(圖 3b)可以看出, 晴天上午的入射短波輻射明顯高于下午, 而上午的出射短波輻射和下午相當(dāng),說明這一時期的地表反照率主要受到下墊面狀況的影響。由典型晴天的風(fēng)速圖(圖4)可以看出, 上午和下午的風(fēng)速較小, 中午前后的風(fēng)速較大, 說明太陽高度角較大(中午前后)時地表反照率不對稱性可能是受風(fēng)的影響。典型晴天的風(fēng)玫瑰圖(圖5)表明, 在觀測區(qū)域, 約 85%的風(fēng)速大于 2 m·s-1, 其主導(dǎo)風(fēng)向為西南風(fēng)。植被葉面受風(fēng)向的影響, 葉傾角發(fā)生改變, 由于地表反照率隨太陽高度角的增加而減小,上午葉面相對太陽的高度角增大, 反照率減小; 下午葉面相對太陽的高度角減小, 地表反照率增大,因此產(chǎn)生了下午地表反照率大于上午的反照率不對稱現(xiàn)象。

圖2 稻田不同時期典型晴天地表反照率的不對稱性Fig.2 Surface albedo asymmetry of paddy fields in typical sunny days in different growth periods of rice

2.2 不同天氣狀況下的地表反照率

選取 10月 9日(晴天)、10月 11日(陰天)和 10月7日(雨天)代表3種不同天氣狀況, 地表反照率調(diào)查結(jié)果(圖 6)表明, 稻田地表反照率整體上晴天＞陰天＞雨天。晴天和陰天的地表反照率不對稱性較為明顯, 均為上午較高, 下午較低; 而雨天的地表反照率較為對稱, 這可能與晴天和陰天上午地表附著較多的露水有關(guān)。晴天地表接收的太陽短波輻射較大,土壤濕度較小, 因此地表反照率較高。受太陽高度角的影響, 地表反照率的日變化較為明顯。陰天到達(dá)地表的太陽短波輻射較少土壤濕度較大, 地表反照率較小, 反照率的日變化幅度最小。雨天條件下, 到達(dá)地表的太陽輻短波射最小, 土壤濕度較大, 反照率也相應(yīng)減小, 地表反照率的日變化較幅度較小。

2.3 水稻不同發(fā)育期內(nèi)的稻田地表反照率

由圖 7可以看出, 從 6月17日水稻插秧到 10月25日水稻收獲, 稻田的地表反照率整體上呈現(xiàn)先升高后降低的趨勢, 與Raja等[23]和Zhang等[24]的研究結(jié)果相同。整個生長季內(nèi), 地表反照率最高值出現(xiàn)在灌漿期, 插秧及成熟期的地表反照率最低。從6月27日開始水稻進(jìn)入分蘗期, 地表反照率呈現(xiàn)緩慢波動上升的趨勢。7月15日水稻進(jìn)入拔節(jié)期, 水稻生長旺盛, 地表反照率迅速上升。8月12日水稻進(jìn)入抽穗期, 地表反照率較大。8月 15日開始, 水稻進(jìn)入灌漿期, 地表反照率呈現(xiàn)波動下降的趨勢。10月 10日開始, 水稻進(jìn)入成熟期, 由于水稻停止生長作物開始枯萎, 地表反照率迅速降低。由圖 8可以看出, 6月底到7月上旬(分蘗期)之間以及9月初到10月上旬(灌漿期)之間降水量較多, 并且這兩個時期的地表反照率波動較大; 而水稻生長季內(nèi)的其他時期, 降水量較小且地表反照率的波動性不明顯,因此, 降水的增多可能是導(dǎo)致分蘗期和灌漿期地表反照率波動較大的一個共同原因。

2.4 稻田反照率與溫度、濕度及短波輻射之間的關(guān)系

2.4.1 不同天氣狀況下稻田反照率與溫度、濕度及短波輻射之間的關(guān)系

圖3 稻田不同時期典型晴天的平均露點溫度與實際溫度(a)和平均短波輻射量(b)的日變化特征Fig.3 Daily dynamics of averages of dew point temperature and actual temperature (a) and average shortwave radiation (RS) (b)of paddy fields in typical sunny days in different growth periods of rice

圖4 稻田不同時期典型晴天平均風(fēng)速日變化Fig.4 Daily dynamics of average wind speed of paddy fields in typical sunny days in different growth periods of rice

圖5 稻田不同時期典型晴天的風(fēng)玫瑰圖Fig.5 Wind rose map of paddy fields in typical sunny days in different growth periods of rice

圖6 不同天氣狀況下的稻田地表反照率Fig.6 Surface albedo of paddy fields under different weather conditions

相關(guān)性分析表明(表2), 晴天地表反照率與地表2 m處的溫度、濕度及短波輻射之間的相關(guān)程度較高, 其中與出射短波輻射顯著相關(guān)(0.670,P＜0.01)。陰天地表反照率與相對濕度之間的相關(guān)程度最高。雨天地表反照率與三者之間的相關(guān)程度均比較低。

2.4.2 不同發(fā)育期稻田反照率與溫度、濕度及短波輻射之間的關(guān)系

1)稻田反照率與溫度及短波輻射變化

在水稻生長初期, 稻田反照率與溫度存在明顯的相反趨勢(圖9)。4月20日到7月上旬地表反照率減小, 地表吸收的太陽輻射增加, 地表溫度上升。9月底, 地表反照率開始升高, 地表吸收的太陽輻射減少, 地表溫度下降。其他階段, 植被生長較為旺盛,地表反照率與地表溫度的關(guān)系較不明顯。

圖7 水稻生長季內(nèi)的稻田地表反照率(Ⅰ為水稻插秧到分蘗期, Ⅱ為水稻分蘗期, Ⅲ為水稻拔節(jié)期, Ⅳ為水稻抽穗期,Ⅴ為水稻灌漿期, Ⅵ為水稻成熟期)Fig.7 Surface albedo of paddy fields during rice growing season (Ⅰ is from transplant stage to tillering stage, Ⅱ is tillering stage,Ⅲ is jointing stage, Ⅳ is heading stage, Ⅴ is filling stage, and Ⅵ is maturity stage)

圖8 水稻生長季內(nèi)的降水量Fig.8 Precipitation during rice growing season

表2 不同天氣狀況下稻田反照率與溫度、濕度及短波輻射之間的相關(guān)關(guān)系Table 2 Correlativity between surface albedo and temperature, humidity and shortwave radiation of paddy field in different weather conditions

圖9 觀測期內(nèi)的稻田地表反照率與溫度的變化Fig.9 Changes of surface albedo and temperature in paddy field during the observation period

地表的入射短波輻射與出射短波輻射在時間上的變化趨勢基本保持一致, 而地表反照率的變化趨勢略滯后于入射和出射短波輻射的變化趨勢(圖 10)。在稻田未整地之前, 地表反照率同入射短波輻射的變化趨勢最為一致, 這一時期下墊面狀況基本保持不變,而入射到地表的短波輻射主要受太陽高度角和天氣狀況影響; 開始整地插秧后, 下墊面改變對出射短波輻射的影響較大, 地表反照率相對入射短波輻射滯后現(xiàn)象較為明顯, 其中水稻生長較為旺盛時期這種滯后現(xiàn)象最為明顯, 說明這一階段, 除了太陽高度角和天氣狀況的影響, 作物的生長狀況對地表反照率的影響也較大。

2)不同發(fā)育期稻田反照率與溫度、濕度及短波輻射之間的相關(guān)關(guān)系

圖10 觀測期內(nèi)的短波輻射量與地表反照率Fig.10 Changes of shortwave radiation (Rs) and surface albedo of paddy fields during the observation period

相關(guān)性分析表明, 整個觀測期內(nèi)地表反照率與溫度及短波輻射之間呈正相關(guān)關(guān)系, 與地表2 m處相對濕度呈負(fù)相關(guān)關(guān)系。水稻生長季內(nèi)地表反照率與濕度及短波輻射之間的相關(guān)程度較高。水稻不同生育期內(nèi)地表反照率與溫度、濕度及短波輻射之間的相關(guān)程度差異較大, 其中灌漿期地表反照率與太陽短波輻射及相對濕度間的相關(guān)程度較高, 并且均通過了P＜0.01顯著性檢驗。相對于水稻生育期, 稻田非生育期內(nèi)(水稻插秧之前和收獲之后)的地表反照率與地表2 m處的溫度相關(guān)程度較高, 與濕度及短波輻射之間的相關(guān)程度較低(表2)。

3 討論

3.1 稻田與其他農(nóng)田的地表反照率對比

已有的相關(guān)研究結(jié)果表明[25-27], 不同下墊面地表反照率差異較大, 其中雪地的地表反照率最大,可以接近0.9左右, 其次為植被覆蓋較少的荒漠、戈壁、退化的草地及裸地等, 在植被覆蓋較多的地區(qū)如農(nóng)田和森林, 地表反照率較小, 森林的地表反照率最小, 農(nóng)田下墊面中因種植作物類型的差異[28],地表反照率也有所不同。干旱區(qū)的麥田反照率大于半濕潤區(qū), 這可能是由于半濕潤區(qū)土壤濕度較大。對比小麥生育期內(nèi)反照率[15-16]和本文稻田生育期內(nèi)反照率的觀測結(jié)果, 發(fā)現(xiàn)麥田和稻田整個生育期內(nèi)的反照率變化規(guī)律大致相同, 均為先升高再減小。余予等[16]對比了小麥田和水稻田的反照率變化, 發(fā)現(xiàn)生長初期小麥田的反照率比水稻田高 0.02, 生長中期小麥田和水稻田的反照率相當(dāng), 孕穗期后水稻田的反照率高于小麥田的反照率, 完全成熟時水稻田的反照率比小麥田反照率高約 0.04。本文的水稻生育期內(nèi)地表反照率同余予等[16]的變化趨勢基本相似, 都是在抽穗期和灌漿期較大, 從插秧到分蘗期和成熟期較小。但是, 在同一生育期內(nèi), 本文的地表反照率較小, 這可能是由于溧水實驗基地位于秦淮河流域內(nèi), 周圍水庫較多, 水分充裕, 土壤濕度較大的緣故。

3.2 影響地表反照率的因子

Liu等[25]發(fā)現(xiàn)玉米(Zea mays)農(nóng)田地表反照率與太陽高度角和土壤濕度呈負(fù)相關(guān)關(guān)系, 與植被指數(shù)呈正相關(guān)關(guān)系。鄭志遠(yuǎn)等[29]發(fā)現(xiàn)裸土地表反照率隨太陽高度角及土壤濕度的增加而減小。張樂樂等[30]研究發(fā)現(xiàn), 土壤水分增加, 地表反照率明顯減小,而且這種關(guān)系還受到植被生長周期的影響。Usowicz等[31]發(fā)現(xiàn)隨著生物炭用量的增加, 地表反照率明顯減小, 并且降雨過后, 休耕地的地表反照率降低,而草地則相反。李慧芳等[32]研究了內(nèi)蒙古錫林浩特市地表反照率, 均發(fā)現(xiàn)地表反照率與歸一化植被指數(shù)(normalized difference vegetation index, NDVI)呈負(fù)相關(guān)關(guān)系。Bsaibes等[33]研究發(fā)現(xiàn), 葉面積指數(shù)(leaf area index, LAI)同地表反照率之間的關(guān)系受土地覆被類型的影響, 且 LAI較小時, 地表反照率的差異較大, 而隨著 LAI的增加, 地表反照率趨于穩(wěn)定。由此可見, 地表反照率的大小受到多種因素的共同影響, 如下墊面的顏色、植被覆蓋率、土壤溫度、土壤濕度、太陽高度角和大氣透過率等。

我們的結(jié)果表明, 典型晴天狀況下, 稻田地表反照率與太陽高度角有關(guān), 白天中午地表反照率最低, 上午和下午地表反照率較大, 因此地表反照率隨太陽高度角的增加而減小。同時, 晴天地表反照率每日內(nèi)的變化特征還受到當(dāng)天晝夜溫差及風(fēng)的影響, 太陽高度角小于60°時, 晝夜溫差導(dǎo)致上午植被表面附著較多露水, 同一太陽高度角下上午地表反照率較下午大; 而太陽高度角小于60°時, 風(fēng)使得植被葉面傾斜, 同一太陽高度角下下午地表反照率較上午大。地表反照率也與天氣狀況有關(guān), 晴天大氣透過率較大, 地表反照率也較大, 而陰雨天大氣透過率較小, 地表反照率也較小。地表反照率還與水稻的生長狀況有關(guān)。水稻生長初期, 地表反照率隨水稻的生長而增加, 在水稻接近成熟時, 稻田地表反照率變化較為平緩。當(dāng)水稻成熟時, 地表反照率開始下降, 水稻收獲之后, 地表反照率進(jìn)一步降低。對于土壤濕度及植被覆蓋度等因子對地表反照率的影響研究, 將是未來需要開展的重點。

表3 水稻不同發(fā)育期地表反照率與溫度、濕度及短波輻射之間的相關(guān)關(guān)系Table 3 Correlativity between surface albedo and temperature, humidity and shortwave radiation of paddy fields in different growth periods of rice

4 結(jié)論

1)晴天稻田地表反照率整體上呈“U”型分布, 中午較低, 下午和上午較高。隨著水稻的生長, 地表反照率的日變化特征變得更為明顯。晴天稻田地表反照率的不對稱性主要是由露水和風(fēng)速風(fēng)向引起, 太陽高度角較小時, 露水的散射作用使得上午的地表反照率值較下午高, 而太陽高度角較大時, 西南風(fēng)促使作物葉面傾斜, 從而使得下午的地表反照率值較上午高。

2)不同天氣狀況下的地表反照率不同, 稻田晴天地表反照率值較陰雨天高。晴天地表反照率變化趨勢較大, 陰天和雨天變化趨勢較小。晴天和陰天時的地表反照率的不對稱性較為明顯。地表反照率在晴天與出射短波輻射相關(guān)系數(shù)最高(0.67,P＜0.01),在陰天與相對濕度之間的相關(guān)程度最高(-0.48,P＜0.05)。

3)生長季內(nèi)地表反照率呈現(xiàn)先升高后降低的趨勢, 并會受到降水的影響。整個生長季內(nèi), 地表反照率最高值出現(xiàn)在灌漿期到成熟期之間, 插秧到分蘗期之間最低。水稻生長季內(nèi)地表反照率與濕度及短波輻射之間的相關(guān)程度較高, 但是不同生育期內(nèi)地表反照率與溫度、濕度及短波輻射之間的相關(guān)程度差異較大, 其中灌漿期地表反照率與太陽短波輻射及濕度間的相關(guān)程度較高, 并且均通過了P＜0.01顯著性檢驗。

References

[1]Russell M J, Nunez M, Chladil M A, et al.Conversion of nadir, narrowband reflectance in red and near-infrared channels to hemispherical surface albedo[J].Remote Sensing of Environment, 1997, 61(1): 16–23

[2]肖登攀, 陶福祿, Moiwo Juana P.全球變化下地表反照率研究進(jìn)展[J].地球科學(xué)進(jìn)展, 2011, 26(11): 1217–1224

Xiao D P, Tao F L, Moiwo Juana P.Research progress on surface albedo under global change[J].Advances in Earth Science, 2011, 26(11): 1217–1224

[3]Betts R A.Offset of the potential carbon sink from boreal forestation by decreases in surface albedo[J].Nature, 2000,408(6809): 187–190

[4]蔣熹.冰雪反照率研究進(jìn)展[J].冰川凍土, 2006, 28(5):728–738

Jiang X.Progress in the research of snow and ice albedo[J].Journal of Glaciology and Geocryology, 2006, 28(5):728–738

[5]Schwaiger H P, Bird D N.Integration of albedo effects caused by land use change into the climate balance: Should we still account in greenhouse gas units?[J].Forest Ecology and Management, 2010, 260(3): 278–286

[6]趙東妮, 王艷華, 任傳友, 等.3種水稻趨勢產(chǎn)量擬合方法的比較分析[J].中國生態(tài)農(nóng)業(yè)學(xué)報, 2017, 25(3): 345–355

Zhao D N, Wang Y H, Ren C Y, et al.Comparative analysis of three fitting methods of rice trend yield[J].Chinese Journal of Eco-Agriculture, 2017, 25(3): 345–355

[7]戴彤, 王靖, 赫迪, 等.1961—2010年氣候變化對西南冬小麥潛在和雨養(yǎng)產(chǎn)量影響的模擬分析[J].中國生態(tài)農(nóng)業(yè)學(xué)報,2016, 24(3): 293–305

Dai T, Wang J, He D, et al.Impact simulation of climate change on potential and rainfed field yields of winter wheat in Southwest China from 1961 to 2010[J].Chinese Journal of Eco-Agriculture, 2016, 24(3): 293–305

[8]Gao B, Jia L, Wang T X.Derivation of land surface albedo at high resolution by combining HJ-1A/B reflectance observations with MODIS BRDF products[J].Remote Sensing, 2014,6(9): 8966–8985

[9]Wang R, Chen J M, Pavlic G, et al.Improving winter leaf area index estimation in coniferous forests and its significance in estimating the land surface albedo[J].ISPRS Journal of Photogrammetry and Remote Sensing, 2016, 119: 32–48

[10]楊娟, 陳洪濱, 王開存, 等.利用 MODIS衛(wèi)星資料分析北京地區(qū)地表反照率時空分布及變化特征[J].遙感技術(shù)與應(yīng)用, 2006, 21(5): 403–406

Yang J, Chen H B, Wang K C, et al.Analysis of the surface albedo distribution and variation in Beijing region by using the MODIS data[J].Remote Sensing Technology and Application, 2006, 21(5): 403–406

[11]張超, 宋海清, 吳國周, 等.NCAR/CLM系列陸面模式對內(nèi)蒙古地表溫度的模擬評估[J].中國生態(tài)農(nóng)業(yè)學(xué)報, 2017,25(2): 299–308

Zhang C, Song H Q, Wu G Z, et al.Simulation and comparative analysis of surface temperature over Inner Mongolia using four NCAR community land models[J].Chinese Journal of Eco-Agriculture, 2017, 25(2): 299–308

[12]Dexter R H, Bailey W G, McArthur L J B.Albedo of wheat during a growing season: Diurnal symmetry and asymmetry[C]// 26thConference on Agricultural and Forest Meteorology.British Columbia, Canada: American Meteorological Society, 2004

[13]陳云浩, 李曉兵, 謝鋒.我國西北地區(qū)地表反照率的遙感研究[J].地理科學(xué), 2001, 21(4): 327–333

Chen Y H, Li X B, Xie F.Study on surface albedo distribution over Northwest China using remote sensing technique[J].Scientia Geographica Sinica, 2001, 21(4): 327–333

[14]蔡福, 祝青林, 何洪林, 等.中國月平均地表反照率的估算及其時空分布[J].資源科學(xué), 2005, 27(1): 114–120

Cai F, Zhu Q L, He H L, et al.Estimation and spatio-temporal distribution of monthly mean surface albedo in China[J].Resources Science, 2005, 27(1): 114–120

[15]張亞峰, 王新平, 潘顏霞, 等.干旱區(qū)春小麥發(fā)育期反照率變化及其不對稱性[J].農(nóng)業(yè)工程學(xué)報, 2011, 27(12):120–125

Zhang Y F, Wang X P, Pan Y X, et al.Variations of surface albedo and its diurnal asymmetry at different growing stages of spring wheat in arid areas[J].Transactions of the CSAE,2011, 27(12): 120–125

[16]余予, 李揚云, 童應(yīng)祥, 等.壽縣地區(qū)小麥和水稻田地表反照率觀測分析[J].氣候與環(huán)境研究, 2009, 14(6): 639–645

Yu Y, Li Y Y, Tong Y X, et al.Observation and analysis of surface albedo of wheat and rice fields in Shouxian region[J].Climatic and Environment Research, 2009, 14(6):639–645

[17]Song J.Diurnal asymmetry in surface albedo[J].Agricultural and Forest Meteorology, 1998, 92(3): 181–189

[18]Ripley E A, Redmann R E.Grassland[M]//Monteith J L.Vegetation and the Atmosphere.London: Academic Press,1976

[19]Grant I F, Prata A J, Cechet R P.The impact of the diurnal variation of albedo on the remote sensing of the daily mean albedo of grassland[J].Journal of Applied Meteorology, 2000,39(2): 231–244

[20]Minnis P, Mayor S, Smith W L, et al.Asymmetry in the diurnal variation of surface albedo[J].IEEE Transactions on Geoscience and Remote Sensing, 1997, 35(4): 879–890

[21]王勝, 張強, 張鴻.稀疏植被地表反照率及土壤熱傳導(dǎo)率特征研究[J].中國沙漠, 2008, 28(1): 119–124

Wang S, Zhang Q, Zhang H.Characteristics of surface albedo and soil heat conductivity in sparse vegetation site[J].Journal of Desert Research, 2008, 28(1): 119–124

[22]文麗娟, 呂世華, 陳世強, 等.干旱區(qū)綠洲地表反照率不對稱觀測研究[J].太陽能學(xué)報, 2009, 30(7): 953–956

Wen L J, Lü S H, Chen S Q, et al.Observation study on diurnal asymmetry in surface albedo of oasis in arid region[J].Acta Energiae Solaris Sinica, 2009, 30(7): 953–956

[23]Raja R, Nayak A K, Rao K S, et al.Effect of fly ash deposition on photosynthesis, growth and yield of rice[J].Bulletin of Environmental Contamination and Toxicology, 2014, 93(1):106–112

[24]Zhang Y F, Wang X P, Pan Y X, et al.Diurnal and seasonal variations of surface albedo in a spring wheat field of arid lands of Northwestern China[J].International Journal of Biometeorology, 2013, 57(1): 67–73

[25]Liu H Z, Tu G, Dong W J.Three-year changes of surface albedo of degraded grassland and cropland surfaces in a semiarid area[J].Chinese Science Bulletin, 2008, 53(8):1246–1254

[26]孫俊, 胡澤勇, 荀學(xué)義, 等.黑河中上游不同下墊面反照率特征及其影響因子分析[J].高原氣象, 2011, 30(3): 607–613

Sun J, Hu Z Y, Xun X Y, et al.Albedo characteristics in different underlying surfaces in the mid and upper-reaches of HEIFE and its impact factor analysis[J].Plateau Meteorology,2011, 30(3): 607–613

[27]陳琛.淮河流域農(nóng)田生態(tài)系統(tǒng)能量平衡與閉合研究[D].合肥: 安徽農(nóng)業(yè)大學(xué), 2012: 1–71

Chen C.Observational study of surface energy balance and energy closure in the farmland ecosystems in Huaihe River Basin[D].Hefei: Anhui Agricultural University, 2012: 1–71

[28]張雯, 謝輝, 張平, 等.扁桃||冬小麥間作模式下樹冠結(jié)構(gòu)對間作區(qū)域光環(huán)境的影響[J].中國生態(tài)農(nóng)業(yè)學(xué)報, 2016,24(6): 753–761

Zhang W, Xie H, Zhang P, et al.Effect of tree canopy structure on light condition in almond-winter wheat intercropping systems[J].Chinese Journal of Eco-Agriculture, 2016, 24(6):753–761

[29]鄭志遠(yuǎn), 韋志剛, 李振朝, 等.敦煌荒漠戈壁地區(qū)裸土地表反照率參數(shù)化研究[J].大氣科學(xué), 2014, 38(2): 297–308

Zheng Z Y, Wei Z G, Li Z C, et al.Study of parameterization of surface albedo of bare soil over the ggobi desert in the Dunhuang Region[J].Chinese Journal of Atmospheric Sciences, 2014, 38(2): 297–308

[30]張樂樂, 趙林, 李韌, 等.青藏高原唐古拉地區(qū)暖季土壤水分對地表反照率及其土壤熱參數(shù)的影響[J].冰川凍土,2016, 38(2): 351–358

Zhang L L, Zhao L, Li R, et al.Investigating the influence of soil moisture on albedo and soil thermodynamic parameters during the warm season in Tanggula Range, Tibetan Plateau[J].Journal of Glaciology and Geocryology, 2016, 38(2):351–358

[31]Usowicz B, Lipiec J, Mateusz ?ukowski M, et al.The effect of biochar application on thermal properties and albedo of loess soil under grassland and fallow[J].Soil and Tillage Research, 2016, 164: 45–51

[32]李慧芳, 袁占良, 余濤, 等.HJ-1/CCD 地表反照率估算及其與NDVI關(guān)系分析[J].遙感信息, 2012, 27(4): 16–21

Li H F, Yuan Z L, Yu T, et al.Surface albedo estimating based on HJ-1/CCD and relationship analysis between albedo and NDVI[J].Remote Sensing Information, 2012, 27(4): 16–21

[33]Bsaibes A, Courault D, Baret F, et al.Albedo and LAI estimates from FORMOSAT-2 data for crop monitoring[J].Remote Sensing of Environment, 2009, 113(4): 716–729

Characteristics of surface albedo in subtropical paddy rice fields*

BAI Xinxin1, HUANG Xiaolin1, QIN Mengsheng1, ZHANG Yuefang2,CHEN Liugen2, GUO Zhi2, HAO Lu1**
(1.Jiangsu Key Laboratory of Agricultural Meteorology, Nanjing University of Information Science & Technology, Nanjing 210044,China; 2.Institute of Agricultural Resources and Environment, Jiangsu Academy of Agricultural Sciences, Nanjing 210014, China)

As a key input parameter of numerical climate models and surface energy balance equations, surface albedo affects climatic systems on the earth.Data observed on rice paddy surface albedo not only better depict energy distribution between ground and atmosphere in paddy fields, but also provide more accurate parameters values for land surface models.It also provides authentication for remote sensing inversions of surface albedo and therefore better explains the effects of land use/cover changemechanism on global climate change.Using CNR4 measured surface albedo in paddy fields at Lishui Experiment Station in 2016,we analyzed the characteristics of surface albedo of paddy field.In the study, we combined the analysis of surface albedo with temperature, humidity, wind speed, wind direction and other meteorological data.Some other possible factors affecting surface albedo in paddy fields were also analyzed.The correlation coefficient between surface albedo and temperature reflected the effect of land surface albedo on temperature.The correlation between surface albedo and meteorological factors (e.g., humidity and solar radiation) was used to study the relationship among surface albedo, temperature and humidity.The results showed that surface albedo had a U-type of distribution for sunny days, which was lowest at noon, highest in the afternoon and morning, but also affected by precipitation.Surface albedo on sunny days was asymmetrical, which was mainly caused by dew, wind speed and wind direction.As actual air temperature was close to dew point temperature at night, more dew settled on leaves in the morning.When solar altitude angle was small, the scattering effect of dew resulted in higher surface albedo in the morning than in the afternoon.When solar altitude angle was high, dew on leaf surfaces gradually evaporated.Southwest winds influenced the tilting of crop leaves, resulting in higher surface albedo in the afternoon than in the morning.Surface albedo was higher on sunny days than on cloudy and rainy days.The correlation coefficient between surface albedo and outgoing shortwave radiation was highest (0.670,P＜ 0.01) on sunny days, while the correlation coefficient between surface albedo and relative humidity was highest (-0.480,P＜0.05) on cloudy days.The trend in surface albedo during the growing season initially increased and then decreased.Surface albedo was highest during grain-filling to mature stage of rice and was lowest during transplanting to tillering stage.The correlation coefficient among surface albedo, humidity and shortwave radiation for rice growing season was high.Also the correlation coefficients between surface albedo and temperature, and then between humidity and shortwave radiation were not the same for different growth stages.During grain-filling period, the correlation among surface albedo, solar shortwave radiation and humidity was more significant, all atP＜ 0.01 significance level.Surface albedo changed greatly at tillering and jointing stages, and were significantly affected by climatic factors.

Mar.31, 2017; accepted Jun.30, 2017

Rice paddy field; Surface albedo; Growing season; Asymmetrical; Meteorological factors

P422

A

1671-3990(2017)11-1707-11

10.13930/j.cnki.cjea.170276

白昕欣, 黃蕭霖, 秦孟晟, 張岳芳, 陳留根, 郭智, 郝璐.亞熱帶地區(qū)水稻田地表反照率變化特征[J].中國生態(tài)農(nóng)業(yè)學(xué)報, 2017, 25(11): 1707-1717

Bai X X, Huang X L, Qin M S, Zhang Y F, Chen L G, Guo Z, Hao L.Characteristics of surface albedo in subtropical paddy rice fields[J].Chinese Journal of Eco-Agriculture, 2017, 25(11): 1707-1717

* 國家自然科學(xué)基金項目(41571026, 71373130)和南京信息工程大學(xué)本科生優(yōu)秀畢業(yè)論文(設(shè)計)支持計劃項目資助

** 通訊作者: 郝璐, 主要研究方向為應(yīng)用氣象。E-mail: hl_haolu@163.com

白昕欣, 主要研究方向為應(yīng)用氣象。E-mail: baiinin@163.com

2017-03-31 接受日期: 2017-06-30

* This study was supported by the National Natural Science Foundation of China (41571026, 71373130) and the Outstanding Undergraduate Thesis (Design) Support Program of Nanjing University of Information Science & Technology.

** Corresponding author, E-mail: hl_haolu@163.com